Disclosure of Invention
The embodiment of the application provides a method, a device and computer equipment for acquiring seismic exploration information, which can improve the accuracy of acquiring the seismic exploration information. The specific technical scheme is as follows:
in one aspect, an embodiment of the present application provides a method for acquiring seismic exploration information, where the method includes:
Acquiring a plurality of seismic data acquired by a plurality of detectors in a target work area, wherein the seismic data are acquired by the detectors in a preset acquisition time range, and the seismic data comprise first seismic data of a first horizontal component, second seismic data of a second horizontal component and third seismic data of a first vertical component, and the first horizontal component, the second horizontal component and the first vertical component are perpendicular to each other;
for each detector, determining a particle vibration complexity parameter, a particle vibration linearity parameter and a cross-correlation parameter of a target acquisition time based on the first seismic data and the second seismic data, wherein the particle vibration complexity parameter is used for reflecting the complexity degree of particle vibration, the particle vibration linearity parameter is used for reflecting the linearity degree of particle vibration, the cross-correlation parameter is used for reflecting the cross-correlation relation of the first seismic data and the second seismic data at the target acquisition time, and the target acquisition time is any acquisition time in the preset acquisition time range;
determining a polarization parameter at the target acquisition time and a polarization angle at the target acquisition time based on the particle vibration complex parameter, the particle vibration linear parameter and the cross-correlation parameter at the target acquisition time;
Determining fourth seismic data based on the first seismic data, the second seismic data, the polarization parameters and the polarization angles at the target acquisition time, wherein the fourth seismic data is the seismic data corresponding to the polarization angles at the target acquisition time;
and acquiring seismic exploration information based on the third seismic data corresponding to each detector and the fourth seismic data at each acquisition time.
In one possible implementation manner, the determining the cross-correlation parameter of the target acquisition time based on the first seismic data and the second seismic data includes:
determining a frequency distribution range of seismic waves excited by shot points based on the first seismic data and the second seismic data;
Based on the frequency distribution range, determining a time window parameter, wherein the time window parameter is used for reflecting the corresponding time window width when the cross-correlation parameter of the target acquisition moment is determined;
determining an integration time range based on the time window parameter and the target acquisition time;
and determining a cross-correlation parameter of the target acquisition moment based on the integration time range, the first seismic data and the second seismic data.
In another possible implementation manner, the determining the cross-correlation parameter of the target acquisition time based on the integration time range, the first seismic data and the second seismic data includes:
Determining the square root of the energy of the first seismic data within the integration time range to obtain a first amplitude;
Determining the square root of the energy of the second seismic data within the integration time range to obtain a second amplitude;
determining the product of the first amplitude and the second amplitude to obtain a third amplitude;
Determining the product of the first seismic data and the second seismic data in the integral time range to obtain a first vibration vector, wherein the first vibration vector comprises an amplitude and a vibration direction;
And determining the absolute value of the ratio of the first vibration vector to the third amplitude to obtain the cross-correlation parameter of the target acquisition moment.
In another possible implementation, the determining the polarization parameter at the target acquisition time and the polarization angle at the target acquisition time based on the particle vibration complexity parameter, the particle vibration linearity parameter, and the cross-correlation parameter at the target acquisition time includes:
determining a polarization parameter of the target acquisition time based on the particle vibration complex parameter, the particle vibration linear parameter and the cross-correlation parameter of the target acquisition time, wherein the polarization parameter is used for reflecting the linear degree of particle vibration;
determining a polarization angle of the target acquisition time based on the first seismic data and the second seismic data under the condition that the polarization parameter is not smaller than a preset polarization parameter;
And under the condition that the polarization parameter is smaller than the preset polarization parameter, determining the polarization angle of the target acquisition moment based on the first seismic data, the second seismic data and an integral time range, wherein the integral time range is obtained when the cross-correlation parameter is determined.
In another possible implementation manner, the determining the polarization angle of the target acquisition time based on the first seismic data and the second seismic data includes:
determining the ratio of the second seismic data at the target acquisition time to the first seismic data at the target acquisition time to obtain a first ratio;
And determining the arctangent value of the first ratio to obtain the polarization angle of the target acquisition time.
In another possible implementation manner, the determining the polarization angle of the target acquisition time based on the first seismic data, the second seismic data and the integration time range includes:
Determining the square root of the sum of the energies of the first seismic data and the second seismic data in the integration time range to obtain a fourth amplitude;
For each preset polarization angle in a preset polarization angle range, determining the product of first seismic data of a target integral moment and a cosine value of the preset polarization angle to obtain fifth seismic data, wherein the target integral moment is any integral moment in the integral time range;
Determining the product of the second seismic data of the target integration moment and the sine value of the preset polarization angle to obtain sixth seismic data;
determining a sum of the fifth seismic data and the sixth seismic data within the integration time range to obtain a fifth amplitude;
Determining the difference value between the fourth amplitude and the fifth amplitude to obtain a sixth amplitude corresponding to the preset polarization angle;
And taking the preset polarization angle corresponding to the minimum sixth amplitude as the polarization angle of the target acquisition moment.
In another possible implementation, the determining the polarization parameter at the target acquisition time based on the particle vibration complexity parameter, the particle vibration linearity parameter, and the cross-correlation parameter at the target acquisition time includes:
determining a difference value between the cross-correlation parameter and the particle vibration linear parameter at the target acquisition time to obtain a first difference value;
Determining the product of the particle vibration complex parameter and the first difference value to obtain a first product;
determining an exponent value based on a natural constant and an exponent based on a negative of the first product;
determining the reciprocal after increasing 1 on the basis of the index value, and taking the reciprocal as the polarization parameter.
In another possible implementation manner, the determining, based on the first seismic data, the second seismic data, the polarization parameter and the polarization angle at the target acquisition time, third seismic data corresponding to the polarization angle at the target acquisition time includes:
Determining the product of the first seismic data at the target acquisition time and the polarization angle cosine value at the target acquisition time to obtain seventh seismic data;
Determining the product of the second seismic data at the target acquisition time and the polarization angle sine value at the target acquisition time to obtain eighth seismic data;
And determining the product of the sum of the seventh seismic data and the eighth seismic data and the polarization parameter at the target acquisition time to obtain third seismic data corresponding to the polarization angle at the target acquisition time.
In another aspect, an embodiment of the present application provides an apparatus for acquiring seismic exploration information, including:
The first acquisition module is used for acquiring a plurality of seismic data acquired by a plurality of detectors in a target work area, wherein the seismic data are acquired by the detectors in a preset acquisition time range, and the seismic data comprise first seismic data of a first horizontal component, second seismic data of a second horizontal component and third seismic data of a first vertical component, and the first horizontal component, the second horizontal component and the first vertical component are perpendicular to each other;
The first determining module is configured to determine, for each detector, a particle vibration complexity parameter, a particle vibration linearity parameter and a cross-correlation parameter of a target acquisition time based on the first seismic data and the second seismic data, where the particle vibration complexity parameter is used to reflect a complexity degree of particle vibration, the particle vibration linearity parameter is used to reflect a linearity degree of particle vibration, the cross-correlation parameter is used to reflect a cross-correlation relationship between the first seismic data and the second seismic data at the target acquisition time, and the target acquisition time is any acquisition time within the preset acquisition time range;
the second determining module is used for determining the polarization parameter at the target acquisition time and the polarization angle at the target acquisition time based on the particle vibration complex parameter, the particle vibration linear parameter and the cross-correlation parameter at the target acquisition time;
The third determining module is used for determining fourth seismic data based on the first seismic data, the second seismic data, the polarization parameters and the polarization angles of the target acquisition time, wherein the fourth seismic data are seismic data corresponding to the polarization angles at the target acquisition time;
the second acquisition module is used for acquiring the seismic exploration information based on the third seismic data corresponding to each detector and the fourth seismic data at each acquisition time.
In one possible implementation manner, the first determining module is configured to determine a frequency distribution range of seismic waves excited by shot points based on the first seismic data and the second seismic data, determine a time window parameter based on the frequency distribution range, where the time window parameter is used to reflect a time window width corresponding to when the cross-correlation parameter of the target acquisition time is determined, determine an integration time range based on the time window parameter and the target acquisition time, and determine the cross-correlation parameter of the target acquisition time based on the integration time range, the first seismic data and the second seismic data.
In another possible implementation manner, the first determining module is configured to determine a square root of energy of the first seismic data in the integration time range to obtain a first amplitude, determine a square root of energy of the second seismic data in the integration time range to obtain a second amplitude, determine a product of the first amplitude and the second amplitude to obtain a third amplitude, determine a product of the first seismic data and the second seismic data in the integration time range to obtain a first vibration vector, where the first vibration vector includes an amplitude and a vibration direction, and determine an absolute value of a ratio of the first vibration vector to the third amplitude to obtain a cross-correlation parameter of the target acquisition time.
In another possible implementation manner, the second determining module is configured to determine a polarization parameter at the target acquisition time based on the particle vibration complex parameter, the particle vibration linearity parameter, and a cross-correlation parameter at the target acquisition time, where the polarization parameter is used to reflect a linearity degree of particle vibration, determine a polarization angle at the target acquisition time based on the first seismic data and the second seismic data if the polarization parameter is not less than a preset polarization parameter, and determine the polarization angle at the target acquisition time based on the first seismic data, the second seismic data, and an integration time range, where the polarization parameter is less than the preset polarization parameter, and the integration time range is obtained when determining the cross-correlation parameter.
In another possible implementation manner, the second determining module is configured to determine a ratio of the second seismic data at the target acquisition time to the first seismic data at the target acquisition time to obtain a first ratio, and determine an arctangent value of the first ratio to obtain a polarization angle at the target acquisition time.
In another possible implementation manner, the second determining module is configured to determine a square root of an energy sum of the first seismic data and the second seismic data in the integration time range to obtain a fourth amplitude, determine, for each preset polarization angle in a preset polarization angle range, a product of a first seismic data at a target integration time and a cosine value of the preset polarization angle to obtain a fifth seismic data, where the target integration time is any integration time in the integration time range, determine a product of a second seismic data at the target integration time and a sine value of the preset polarization angle to obtain a sixth seismic data, determine a sum value of the fifth seismic data and the sixth seismic data in the integration time range to obtain a fifth amplitude, determine a difference value between the fourth amplitude and the fifth amplitude to obtain a sixth amplitude corresponding to the preset polarization angle, and use a preset polarization angle corresponding to the minimum sixth amplitude as the polarization angle at the target acquisition time.
In another possible implementation manner, the second determining module is configured to determine a difference value between the cross-correlation parameter at the target acquisition time and the particle vibration linearity parameter to obtain a first difference value, determine a product of the particle vibration complex parameter and the first difference value to obtain a first product, determine an exponent value with a natural constant as a base and a negative number of the first product as an exponent, determine an inverse number after adding 1 on the basis of the exponent value, and use the inverse number as the polarization parameter.
In another possible implementation manner, the third determining module is configured to determine a product of the first seismic data at the target acquisition time and a cosine value of a polarization angle at the target acquisition time to obtain seventh seismic data, determine a product of the second seismic data at the target acquisition time and a sine value of the polarization angle at the target acquisition time to obtain eighth seismic data, and determine a product of a sum of the seventh seismic data and the eighth seismic data and a polarization parameter at the target acquisition time to obtain third seismic data corresponding to the polarization angle at the target acquisition time.
In another aspect, a computer device is provided, where the computer device includes a processor and a memory, where the memory stores at least one piece of program code that is loaded and executed by the processor to implement operations performed in the seismic survey information acquisition method according to an embodiment of the present application.
In another aspect, embodiments of the present application provide a computer readable storage medium having at least one program code stored therein, the at least one program code being loaded and executed by a processor to implement operations performed in the method for acquiring seismic survey information in embodiments of the present application.
In another aspect, embodiments of the present application provide a computer program product or a computer program comprising computer program code stored in a computer readable storage medium. The processor of the computer device reads the computer program code from the computer readable storage medium and executes the computer program code to perform the operations performed in the seismic survey information acquisition method according to the embodiments of the application.
The technical scheme provided by the embodiment of the application has the beneficial effects that:
The embodiment of the application provides a seismic exploration information acquisition method, which is characterized in that the vibration condition of particles at different acquisition moments is complex and changeable due to the complexity of the underground structural form and lithology, the method firstly determines the polarization angle and the polarization parameter corresponding to each acquisition moment in a preset acquisition time range, and determines the vibration condition of the particles at each acquisition moment according to the polarization angle and the polarization parameter corresponding to each acquisition moment, so that the real vibration condition of the mass points in the underground can be reflected, the seismic exploration information can be accurately acquired according to the vibration condition of the particles at each acquisition moment, and the accuracy of the acquired seismic exploration information is improved.
Detailed Description
In order to make the technical scheme and advantages of the present application more clear, the following further describes the embodiments of the present application in detail.
For ease of explanation, the case of particle vibration caused by seismic waves will be described.
The shot point excitation generates a seismic wave, and when the seismic wave propagates to the position of the three-component detector, the particle is caused to vibrate. When a particle vibrates, the particle leaves the equilibrium position and three-dimensional vibration is completed. After the seismic wave passes through the position of the three-component detector, the particles are restored to the original equilibrium position. The three-component wave detectors can record the vibration condition of the particles in the three-dimensional space, which is caused by the seismic waves, comprehensively reflect the propagation rule and characteristics of the seismic waves in the underground medium, and accurately reconstruct the three-dimensional space vibration tracks of the particles by utilizing the three-component seismic data collected by the three-component wave detectors, wherein the vibration tracks are a set of the points arranged in the three-dimensional space in time sequence and can be expressed as a function of time. Referring to FIG. 1, FIG. 1 shows the trajectories of particles at different times on the X-horizontal, Y-horizontal, and Z-vertical components, respectively. And connecting all vibration track points according to time sequence, and representing the vibration track points as a curve in a three-dimensional space. In different time periods, vibration tracks of particles are different due to different types of seismic waves reaching the three-component detectors, wherein the types of the seismic waves mainly comprise longitudinal waves, transverse waves, converted waves, rayleigh waves, lux waves, random noise and the like. Referring to fig. 2 to 4, fig. 2 is a vibration trace of a particle on an X horizontal component and a Y horizontal component caused by random noise, fig. 3 is a vibration trace of a particle on an X horizontal component and a Y horizontal component caused by a longitudinal wave of an earthquake, and fig. 4 is a vibration trace of a particle on an X horizontal component and a Z vertical component caused by a rayleigh wave, and it can be seen from fig. 2 to 4 that vibration traces of particles caused by different types of seismic waves are different. The two main characteristics of the particle's trajectory are amplitude and direction of vibration, and different types of seismic waves have different polarization characteristics, such as linear, elliptical, and three-dimensional. The polarization characteristics of particle vibrations are not only dependent on the type of seismic wave, but also affected by the type of source and the complexity of the subsurface medium.
Although the polarization characteristics of different types of seismic waves are different, the particle vibrations caused by the seismic waves at the location of the three-component geophone can be determined from the vibration vectors, which vary in size and direction. The trajectory of the particles may also be represented by different parameters related to the selected coordinate system. For example, in a Cartesian coordinate system, the projections of the vibration vector A (T) on the coordinate axes of X, Y and Z can be denoted as Ax(T)、Ay (T) and Az (T), respectively, and the following formulas hold:
Ax(T)=|A(T)|sinφ(T)cosω(T)
Ay(T)=|A(T)|sinφ(T)sinω(T)
Az(T)=|A(T)|cosφ(T)
Where A (T) represents the absolute value of the displacement of the particle polarization at the time of T acquisition, ω (T) represents the polarization angle between the particle polarization at the time of T acquisition and the X-axis, φ (T) represents the angle between the particle polarization at the time of T acquisition and the vertical plane, Ax (T) represents the particle vibration displacement recorded by the sensor for the X-horizontal component, Ay (T) represents the particle vibration displacement recorded by the sensor for the y-horizontal component, and Az (T) represents the particle vibration displacement recorded by the sensor for the z-vertical component. Referring to FIG. 5, FIG. 5 is a schematic diagram of the projection and angle of particle vibrations in the XOY plane and the XOZ plane, where Ao (T) is the projection of A (T) in the XOY plane, A1 (T) is the projection of A (T) in the XOZ plane, ω (T) is the angle of Ao (T) with the X-axis direction, and φ (T) is the angle of A1 (T) with the Z-axis direction.
When the seismic longitudinal wave recorded by the three-component detector does not interfere with other waves, the vibration of the particles caused by the seismic longitudinal wave is linearly polarized, wherein the linear polarization means that the direction of a vibration vector A (T) is unchanged, and the size of the vibration vector A (T) can be changed, namely the particles vibrate in a linear track mode near the balance position. In a homogeneous isotropic medium, the seismic transverse wave is linearly polarized in a plane tangential to the seismic longitudinal wave front. The rayleigh wave is elliptically polarized, the plane of polarization is plumb, and the amplitudes of the horizontal and vertical components of the rayleigh wave vary differently with depth. Unlike rayleigh waves, the lux surface wave is linearly polarized in a horizontal plane perpendicular to the wave propagation direction. Random noise is characterized by instability of the polarization characteristics, which means that its phase difference between the different components of the three-component detector is random, the different components can be in phase and the recorded shape is similar in a shorter acquisition time range, but the cross-correlation parameters of the seismic data recorded by the sensors of the different components of the same random noise in the three-component detector are approximately zero in a larger acquisition time range, as continued reference is made to fig. 2.
In order to improve the processing interpretation efficiency and accuracy of the three-component seismic data, the projection of the particle vibration in the XOY horizontal plane at the position of the three-component detector is defined as the optimal receiving component, namely the maximum vibration amplitude, and the received seismic signal is strongest.
The following will describe the scheme of the present application in detail:
the embodiment of the application provides a method for acquiring seismic exploration information, which is executed by computer equipment, and is shown in fig. 6, and the method comprises the following steps:
step 601, a computer device obtains a plurality of seismic data acquired by a plurality of detectors within a target work area.
For each detector, the seismic data is acquired by the detector within a preset acquisition time range, and the seismic data comprises first seismic data of a first horizontal component, second seismic data of a second horizontal component and third seismic data of a first vertical component, wherein the first horizontal component, the second horizontal component and the first vertical component are perpendicular to each other.
Referring to fig. 7 and 8, fig. 7 is first seismic data of a theoretical synthesized X horizontal component, and fig. 8 is second seismic data of a theoretical synthesized Y horizontal component.
Step 602, for each detector, the computer device determines a particle vibration complex parameter, a particle vibration linearity parameter, and a cross-correlation parameter for a target acquisition time based on the first seismic data and the second seismic data.
The particle vibration complex parameter is used for reflecting the complexity degree of particle vibration, the particle vibration linear parameter is used for reflecting the linearity degree of particle vibration, the cross-correlation parameter is used for reflecting the cross-correlation relation between the first seismic data and the second seismic data at the target acquisition time, and the target acquisition time is any acquisition time within a preset acquisition time range.
Step 602 may be implemented by the following steps (1) to (2), including:
(1) The computer device determines cross-correlation parameters for the target acquisition time based on the first seismic data and the second seismic data.
Step (1) may be achieved by the following steps (1-1) to (1-4), comprising:
(1-1) the computer device determining a frequency distribution range of the shot-excited seismic waves based on the first seismic data and the second seismic data.
The shot point excited seismic waves have seismic longitudinal waves and seismic transverse waves, and the computer equipment acquires the frequency distribution range of the seismic longitudinal waves and the seismic transverse waves from the first seismic data and the second seismic data.
(1-2) The computer device determining the time window parameter based on the frequency distribution range.
The time window parameters are used for reflecting the corresponding time window width when the cross-correlation parameters of the target acquisition time are determined. The computer device may determine the time window parameters based on the frequency distribution range of the seismic longitudinal wave, or may determine the time window parameters based on the frequency distribution range of the seismic transverse wave. Or the time window parameter may also be half the time window width, i.e. half the time window.
The description will be given by taking the example of determining the time window parameter based on the frequency distribution range of the seismic transverse wave by the computer device, and determining the main frequency of the seismic transverse wave based on the frequency distribution range of the seismic transverse wave by the computer device, and taking the inverse of the main frequency as the time window parameter. The main frequency of the transverse wave refers to the frequency with the strongest energy of the transverse wave. For example, if the main frequency of the transverse seismic wave is 10Hz, the time window parameter is 0.1. The computer device may also determine the time window parameters by other methods, which are not particularly limited.
(1-3) The computer device determining an integration time range based on the time window parameter and the target acquisition time instant.
The computer equipment takes the sum value of the time window parameter and the target acquisition time as the upper limit value of the integration time range, and takes the difference value of the time window parameter and the target acquisition time as the lower limit value of the integration time range, so as to obtain the integration time range.
(1-4) The computer device determining a cross-correlation parameter for the target acquisition time based on the integration time range, the first seismic data, and the second seismic data.
Step (1-4) may be achieved by the following steps (1-4-1) to (1-4-5), comprising:
(1-4-1) the computer device determining a square root of the energy of the first seismic data over the integration time range to obtain a first amplitude.
The first amplitude may be expressed by the following formula: Wherein t0 denotes any integration time within the integration time range from t-n to t+n, t denotes the target acquisition time, n denotes the time window parameter, ax(t0) denotes the first seismic data at time t0, x denotes the first horizontal component,Representing the energy of the first seismic data over the integration time range.
(1-4-2) The computer device determining a square root of the energy of the second seismic data over the integration time range to obtain a second amplitude.
The second amplitude may be expressed by the following formula: Where ay(t0) represents the second seismic data at the integration time of t0, y represents the second horizontal component,Representing the energy of the second seismic data over the integration time range.
(1-4-3) The computer device determining a product of the first amplitude and the second amplitude to obtain a third amplitude.
The third amplitude may be expressed by the following formula:
(1-4-4) the computer device determining a product of the first seismic data and the second seismic data over the integration time range to obtain a first vibration vector.
The first vibration vector includes an amplitude and a vibration direction, and the first vibration vector can be expressed by the following formula:
(1-4-5) the computer device determining an absolute value of a ratio of the first vibration vector to the third amplitude to obtain a cross-correlation parameter for the target acquisition time.
The cross-correlation parameter may be expressed by the following formula:
where c (t) represents the cross-correlation parameter at the time of t acquisition.
(2) The computer device determines a particle vibration complex parameter and a particle vibration linearity parameter based on the first seismic data and the second seismic data.
The complex particle vibration parameters and the linear particle vibration parameters are constants, and the computer equipment can determine the complex particle vibration parameters and the linear particle vibration parameters based on the linear particle vibration conditions reflected by the first seismic data and the second seismic data, wherein the more obvious the linear particle vibration is, the more prominent the complex particle vibration parameters are, and the smaller the linear particle vibration parameters are.
The particle vibration complexity parameter is used to control the form of the polarization parameter, and may be represented as d0,d0 e 1,10000, in most cases, d0 may take on a value of 10. The particle vibration linearity parameter, which is used as a threshold value for controlling the enhancement or suppression of seismic waves with linear polarization characteristics, can be expressed as c0,c0 epsilon 0, 1.
Step 603, the computer device determines the polarization parameter at the target acquisition time based on the particle vibration complex parameter, the particle vibration linear parameter, and the cross-correlation parameter at the target acquisition time.
Step 603 may be implemented by the following steps (1) to (4), including:
(1) The computer equipment determines the difference value of the cross-correlation parameter and the particle vibration linear parameter at the target acquisition time to obtain a first difference value.
The first difference may be represented by the following formulas c (t) -c0, where c (t) represents the cross-correlation parameter at the time of t acquisition and c0 represents the particle vibration linearity parameter.
(2) The computer device determines a product of the particle vibration complex parameter and the first difference value to obtain a first product.
The first product may be represented by d0(c(t)-c0, where d0 represents a particle vibration complex parameter.
(3) The computer device determines an exponent value based on the natural constant and an exponent based on a negative number of the first product.
The index value may be expressed by the following formula: where e is a natural constant, which is approximately equal to 2.71828.
(4) The computer device determines the reciprocal after increasing by 1 on the basis of the index value, and takes the reciprocal as the polarization parameter of the target acquisition time.
The polarization parameter is used to reflect the linearity of particle vibration and can be expressed by the following formula:
Wherein γ (t) represents the polarization parameter at the target acquisition time.
After obtaining the polarization parameter at the target acquisition time, the computer device determines the magnitude relation between the polarization parameter and the preset polarization parameter, executes step 604 if the polarization parameter is not less than the preset polarization parameter, and executes step 605 if the polarization parameter is less than the preset polarization parameter.
The preset polarization parameters may be set and changed according to needs, and in the embodiment of the present application, this is not particularly limited. For example, the preset polarization parameter is 0.8 or 0.9.
Note that γ (t) ∈ [0,1], and when γ (t) =1, ax (t) and ay (t) can be reconstructed. Considering that petroleum and natural gas exploration utilizes information of seismic longitudinal waves and seismic transverse waves to detect a subsurface formation, both of which have the characteristic of linear polarization, γ (t) is an activation function having the absolute value of the cross-correlation function of ax (t) and ay (t) as an argument.
Step 604, under the condition that the polarization parameter of the target acquisition time is not smaller than the preset polarization parameter, the computer equipment determines the polarization angle of the target acquisition time based on the first seismic data and the second seismic data.
Under the condition that the polarization parameter at the target acquisition time is not smaller than the preset polarization parameter, the characteristic of particle vibration is typical linear polarization, the computer equipment determines the ratio of the second seismic data at the target acquisition time to the first seismic data at the target acquisition time to obtain a first ratio, and determines the arctangent value of the first ratio to obtain the polarization angle at the target acquisition time.
In the case that the polarization parameter at the target acquisition time is not less than the preset polarization parameter, the polarization angle can be expressed by the following formula: Wherein ω (t) represents the polarization angle at the time of t-acquisition, arctan (°) represents the arctangent function, ay (t) represents the second seismic data at the time of t-acquisition, and ax (t) represents the first seismic data at the time of t-acquisition.
Step 605, under the condition that the polarization parameter of the target acquisition time is smaller than the preset polarization parameter, the computer equipment determines the polarization angle of the target acquisition time based on the first seismic data, the second seismic data and the integration time range.
The integration time range is obtained when the cross-correlation parameter is determined in step 602, and when the polarization parameter is smaller than the preset polarization parameter, it is indicated that the particle vibration may be caused by superposition of two or more types of linearly polarized seismic waves, or caused by an elliptically polarized rayleigh surface wave or random noise, etc., and the polarization angle needs to be determined by a scanning analysis method. Accordingly, step 605 is implemented by the following steps (1) to (6), including:
(1) The computer device determines a square root of an energy sum of the first seismic data and the second seismic data over an integration time range to obtain a fourth amplitude.
The fourth amplitude may be expressed by the following formula:
(2) For any preset polarization angle in the preset polarization angle range, the computer equipment determines the product of the first seismic data of the target integral moment and the cosine value of the preset polarization angle to obtain fifth seismic data.
The target integration time is any integration time within the integration time range, and the preset polarization angle range isI.e. the
The fifth seismic data may be represented by ax(t0) cos ω (β), where t0 represents the target integration time, cos ω (β) represents a preset polarization angle cosine value, β represents a preset polarization angle, ω represents a polarization angle, and cos ω (β) represents a cosine value when ω is β.
(3) The computer equipment determines the product of the second seismic data at the target integration time and the sine value of the preset polarization angle to obtain sixth seismic data.
The sixth seismic data may be represented by ay(t0) sin ω (β), where sin ω (β) represents a sine value when ω is β.
(4) The computer device determines a sum of the fifth seismic data and the sixth seismic data over an integration time range to obtain a fifth amplitude.
The fifth amplitude may be expressed by the following formula:
(5) The computer equipment determines the difference value of the fourth amplitude and the fifth amplitude to obtain a sixth amplitude corresponding to the preset polarization angle.
The sixth amplitude may be expressed by the following formula:
(6) And the computer equipment takes the preset polarization angle corresponding to the minimum sixth amplitude as the polarization angle of the target acquisition moment.
The computer device may determine a sixth amplitude corresponding to each preset polarization angle in the preset polarization angle range through the steps (1) to (5), and the computer device uses the preset polarization angle corresponding to the smallest sixth amplitude as the polarization angle of the target acquisition time.
The polarization angle at the target acquisition time can be expressed as:
Referring to fig. 9, fig. 9 is a schematic diagram of synthesizing the first seismic data in fig. 7 and the second seismic data in fig. 8 into third seismic data based on the polarization angle and the polarization parameter at the target acquisition time, and different colors of fig. 9 represent different polarization angles.
Step 606, the computer device determines fourth seismic data based on the first seismic data, the second seismic data, the polarization parameters and the polarization angle at the target acquisition time.
The fourth seismic data is the seismic data corresponding to the polarization angle at the target acquisition time, and accordingly, step 606 may be implemented by the following steps (1) to (3), including:
(1) The computer equipment determines the product of the first seismic data at the target acquisition time and the polarization angle cosine value at the target acquisition time to obtain seventh seismic data.
The first seismic data are corresponding data in a preset acquisition time range, and the computer equipment acquires the first seismic data of the target acquisition time.
The seventh seismic data may be represented by ax (t) cos ω (t), where cos ω (t) is the polarization angle cosine value at the time of t acquisition.
(2) The computer equipment determines the product of the second seismic data at the target acquisition time and the polarization angle sine value at the target acquisition time to obtain eighth seismic data.
The second seismic data are corresponding data in a preset acquisition time range, and the computer equipment acquires the second seismic data of the target acquisition time.
The eighth seismic data may be represented by ay (t) sin ω (t), where sin ω (t) is the polarization angle sine value at the time of t acquisition.
(3) The computer equipment determines the product of the sum value of the seventh seismic data and the eighth seismic data and the polarization parameter of the target acquisition time to obtain fourth seismic data.
The fourth seismic data may be represented by the following formula:
Ao(t)＝γ(t)(Ax(t)cosω(t)+Ay (t) sin ω (t)), where ao (t) is fourth seismic data corresponding to ω polarization angle at the t acquisition time, and γ (t) is a polarization parameter at the t acquisition time.
It should be noted that, the method provided by the embodiment of the application is mainly applied to processing and interpretation of multi-wave and multi-component seismic data, and can shorten the production period of the processing and interpretation work of the multi-wave and multi-component seismic data and improve the accuracy of the processing and interpretation result. According to the method, the polarization angle of particle vibration at the three-component detector is determined by analyzing the three-component seismic data, and the optimal receiving component is synthesized, so that the processing interpretation efficiency and the processing interpretation precision of the three-component seismic data are improved. Compared with the method in the related art, the method fully considers the complexity of underground seismic wave propagation, can improve the detection efficiency and accuracy of complex stratum complex media, and reduces the risk of petroleum and natural gas exploration and development.
Referring to fig. 10, fig. 10 is a schematic diagram of third seismic data obtained by combining the first seismic data in fig. 7 and the second seismic data in fig. 8, where the third seismic data in fig. 10 is more complete than the first seismic data in fig. 7 and the second seismic data in fig. 8, and the third seismic data can more truly and fully reflect particle vibration caused by seismic waves.
Referring to fig. 11, fig. 11 shows third seismic data corresponding to different polarization angles at different acquisition moments, and fig. 11 can intuitively represent projection of particle vibration on the XOY plane, that is, the third seismic data and the polarization angle. Referring to fig. 12, fig. 12 is a polarization angle determined from three-component seismic data actually acquired, with different colors representing different polarization angles. Referring to fig. 13, fig. 13 is a third seismic data determined from the three-component seismic data actually acquired, and it can be seen from fig. 12 and 13 that the different types of seismic waves have different polarization angles.
In step 607, the computer device obtains seismic exploration information based on the third seismic data corresponding to each detector and the fourth seismic data at each acquisition time.
The third seismic data is seismic data of the first vertical component, and the computer equipment acquires seismic exploration information based on the third seismic data of each acquisition time and the fourth seismic data of each acquisition time corresponding to each detector. For each acquisition time corresponding to each detector, the computer device may determine a vibration trace of the particle based on the fourth seismic data at the acquisition time, and determine a vibration trace of the particle on the first vertical component based on the third seismic data at the acquisition time, and acquire the seismic exploration information according to the vibration trace on the polarization angle and the vibration trace on the first vertical component.
In the embodiment of the application, after the computer equipment acquires the fourth seismic data of each detector at each acquisition time, the fourth seismic data of each acquisition time and the polarization angle of each acquisition time can be correspondingly stored so as to guide subsequent processing.
The embodiment of the application provides a seismic exploration information acquisition method, which is characterized in that the vibration condition of particles at different acquisition moments is complex and changeable due to the complexity of the underground structural form and lithology, the method firstly determines the polarization angle and the polarization parameter corresponding to each acquisition moment in a preset acquisition time range, and determines the vibration condition of the particles at each acquisition moment according to the polarization angle and the polarization parameter corresponding to each acquisition moment, so that the real vibration condition of the mass points in the underground can be reflected, the seismic exploration information can be accurately acquired according to the vibration condition of the particles at each acquisition moment, and the accuracy of the acquired seismic exploration information is improved.
The embodiment of the application provides a device for acquiring seismic exploration information, referring to fig. 14, the device comprises:
A first obtaining module 1401, configured to obtain a plurality of seismic data collected by a plurality of detectors in a target work area, where the seismic data is collected by the detectors in a preset collection time range, and the seismic data includes first seismic data of a first horizontal component, second seismic data of a second horizontal component, and third seismic data of a first vertical component, where the first horizontal component, the second horizontal component, and the first vertical component are perpendicular to each other;
a first determining module 1402, configured to determine, for each detector, a particle vibration complexity parameter, a particle vibration linearity parameter, and a cross-correlation parameter of a target acquisition time based on the first seismic data and the second seismic data, where the particle vibration complexity parameter is used to reflect a complexity degree of particle vibration, the particle vibration linearity parameter is used to reflect a linearity degree of particle vibration, the cross-correlation parameter is used to reflect a cross-correlation relationship between the first seismic data and the second seismic data at the target acquisition time, and the target acquisition time is any acquisition time within a preset acquisition time range;
A second determining module 1403, configured to determine a polarization parameter at the target acquisition time and a polarization angle at the target acquisition time based on the particle vibration complex parameter, the particle vibration linear parameter, and the cross-correlation parameter at the target acquisition time;
a third determining module 1404, configured to determine fourth seismic data based on the first seismic data, the second seismic data, the polarization parameter and the polarization angle at the target acquisition time, where the fourth seismic data is seismic data corresponding to the polarization angle at the target acquisition time;
A second acquiring module 1405, configured to acquire seismic exploration information based on the third seismic data corresponding to each detector and the fourth seismic data at each acquisition time.
In one possible implementation, a first determining module 1402 is configured to determine a frequency distribution range of seismic waves excited by shots based on the first seismic data and the second seismic data, determine a time window parameter based on the frequency distribution range, the time window parameter being configured to reflect a corresponding time window width when the cross-correlation parameter at the target acquisition time is determined, determine an integration time range based on the time window parameter and the target acquisition time, and determine the cross-correlation parameter at the target acquisition time based on the integration time range, the first seismic data, and the second seismic data.
In another possible implementation, the first determining module 1402 is configured to determine a square root of energy of the first seismic data in an integration time range to obtain a first amplitude, determine a square root of energy of the second seismic data in the integration time range to obtain a second amplitude, determine a product of the first amplitude and the second amplitude to obtain a third amplitude, determine a product of the first seismic data and the second seismic data in the integration time range to obtain a first vibration vector, where the first vibration vector includes an amplitude and a vibration direction, and determine an absolute value of a ratio of the first vibration vector to the third amplitude to obtain a cross-correlation parameter of the target acquisition time.
In another possible implementation manner, the second determining module 1403 is configured to determine a polarization parameter at the target acquisition time based on the particle vibration complex parameter, the particle vibration linearity parameter, and the cross-correlation parameter at the target acquisition time, where the polarization parameter is used to reflect the linearity degree of the particle vibration, determine a polarization angle at the target acquisition time based on the first seismic data and the second seismic data if the polarization parameter is not less than a preset polarization parameter, and determine the polarization angle at the target acquisition time based on the first seismic data, the second seismic data, and an integration time range if the polarization parameter is less than the preset polarization parameter, where the integration time range is obtained when determining the cross-correlation parameter.
In another possible implementation manner, the second determining module 1403 is configured to determine a ratio of the second seismic data at the target acquisition time to the first seismic data at the target acquisition time to obtain a first ratio, and determine an arctangent value of the first ratio to obtain the polarization angle at the target acquisition time.
In another possible implementation manner, the second determining module 1403 is configured to determine a square root of an energy sum of the first seismic data and the second seismic data in the integration time range to obtain a fourth amplitude, determine, for each preset polarization angle in the preset polarization angle range, a product of the first seismic data at a target integration time and a cosine value of the preset polarization angle to obtain a fifth seismic data, the target integration time being any integration time in the integration time range, determine a product of the second seismic data at the target integration time and a sine value of the preset polarization angle to obtain a sixth seismic data, determine a sum of the fifth seismic data and the sixth seismic data in the integration time range to obtain a fifth amplitude, determine a difference between the fourth amplitude and the fifth amplitude to obtain a sixth amplitude corresponding to the preset polarization angle, and use a preset polarization angle corresponding to the minimum sixth amplitude as the polarization angle at the target acquisition time.
In another possible implementation, the second determining module 1403 is configured to determine a difference between the cross-correlation parameter and the particle vibration linearity parameter at the target acquisition time to obtain a first difference, determine a product of the particle vibration complex parameter and the first difference to obtain a first product, determine an exponent value based on a natural constant and an exponent value based on a negative of the first product, and determine a reciprocal after adding 1 based on the exponent value, where the reciprocal is used as the polarization parameter.
In another possible implementation manner, the third determining module 1404 is configured to determine a product of the first seismic data at the target acquisition time and a cosine value of a polarization angle at the target acquisition time to obtain seventh seismic data, determine a product of the second seismic data at the target acquisition time and a sine value of the polarization angle at the target acquisition time to obtain eighth seismic data, and determine a product of a sum of the seventh seismic data and the eighth seismic data and a polarization parameter at the target acquisition time to obtain third seismic data corresponding to the polarization angle at the target acquisition time.
The embodiment of the application provides a seismic exploration information acquisition device, which is characterized in that the vibration condition of particles at different acquisition moments is complex and changeable due to the complexity of the underground structural form and lithology, the device firstly determines the polarization angle and the polarization parameter corresponding to each acquisition moment in a preset acquisition time range, and determines the vibration condition of the particles at each acquisition moment according to the polarization angle and the polarization parameter corresponding to each acquisition moment, so that the real vibration condition of the mass points in the underground can be reflected, the seismic exploration information can be accurately acquired according to the vibration condition of the particles at each acquisition moment, and the accuracy of the acquired seismic exploration information is improved.
Fig. 15 shows a block diagram of a computer device 1500 provided in accordance with an exemplary embodiment of the present application. The computer device 1500 may be a portable mobile computer device such as a smart phone, tablet, MP3 player (Moving Picture Experts Group Audio Layer III, MPEG 3), MP4 (Moving Picture Experts Group Audio Layer IV, MPEG 4) player, notebook, or desktop. The computer device 1500 may also be referred to as a user device, portable computer device, laptop computer device, desktop computer device, and the like.
In general, computer device 1500 includes a processor 1501 and memory 1502.
The processor 1501 may include one or more processing cores, such as a 4-core processor, an 8-core processor, or the like. The processor 1501 may be implemented in at least one hardware form of DSP (DIGITAL SIGNAL Processing), FPGA (Field-Programmable gate array), PLA (Programmable Logic Array ). The processor 1501 may also include a main processor, which is a processor for processing data in a wake-up state, also called a CPU (Central Processing Unit ), and a coprocessor, which is a low-power processor for processing data in a standby state. In some embodiments, the processor 1501 may be integrated with a GPU (Graphics Processing Unit, image processor) for rendering and drawing of content to be displayed by the display screen. In some embodiments, the processor 1501 may also include an AI (ARTIFICIAL INTELLIGENCE ) processor for processing computing operations related to machine learning.
Memory 1502 may include one or more computer-readable storage media, which may be non-transitory. Memory 1502 may also include high-speed random access memory, as well as non-volatile memory, such as one or more magnetic disk storage devices, flash memory storage devices. In some embodiments, a non-transitory computer readable storage medium in memory 1502 is used to store at least one instruction for execution by processor 1501 to implement the method of acquiring seismic survey information provided by a method embodiment of the present application.
In some embodiments, computer device 1500 may also optionally include a peripheral device interface 1503 and at least one peripheral device. The processor 1501, memory 1502 and peripheral interface 1503 may be connected by a bus or signal lines. The individual peripheral devices may be connected to the peripheral device interface 1503 via a bus, signal lines, or circuit board. Specifically, the peripheral devices include at least one of radio frequency circuitry 1504, a display screen 1505, a camera assembly 1506, audio circuitry 1507, a positioning assembly 1508, and a power supply 1509.
A peripheral interface 1503 may be used to connect I/O (Input/Output) related at least one peripheral device to the processor 1501 and the memory 1502. In some embodiments, the processor 1501, the memory 1502 and the peripheral interface 1503 are integrated on the same chip or circuit board, and in some other embodiments, either or both of the processor 1501, the memory 1502 and the peripheral interface 1503 may be implemented on separate chips or circuit boards, which is not limited in this embodiment.
The Radio Frequency circuit 1504 is configured to receive and transmit RF (Radio Frequency) signals, also known as electromagnetic signals. The radio frequency circuit 1504 communicates with a communication network and other communication devices via electromagnetic signals. The radio frequency circuit 1504 converts electrical signals to electromagnetic signals for transmission, or converts received electromagnetic signals to electrical signals. Optionally, the radio frequency circuit 1504 includes an antenna system, an RF transceiver, one or more amplifiers, tuners, oscillators, digital signal processors, codec chipsets, subscriber identity module cards, and so forth. The radio frequency circuit 1504 may communicate with other computer devices via at least one wireless communication protocol. The wireless communication protocols include, but are not limited to, the world wide web, metropolitan area networks, intranets, various generations of mobile communication networks (2G, 3G, 4G, and 5G), wireless local area networks, and/or WiFi (WIRELESS FIDELITY ) networks. In some embodiments, the radio frequency circuit 1504 may further include NFC (NEAR FIELD Communication) related circuits, which is not limited by the present application.
Display 1505 is used to display a UI (User Interface). The UI may include graphics, text, icons, video, and any combination thereof. When display screen 1505 is a touch display screen, display screen 1505 also has the ability to collect touch signals at or above the surface of display screen 1505. The touch signal may be input to the processor 1501 as a control signal for processing. At this point, display 1505 may also be used to provide virtual buttons and/or a virtual keyboard, also referred to as soft buttons and/or a soft keyboard. In some embodiments, display screen 1505 may be one, disposed on a front panel of computer device 1500, in other embodiments, display screen 1505 may be at least two, disposed on different surfaces of computer device 1500 or in a folded design, respectively, and in other embodiments, display screen 1505 may be a flexible display screen, disposed on a curved surface or a folded surface of computer device 1500. Even more, the display 1505 may be arranged in a non-rectangular irregular pattern, i.e., a shaped screen. The display screen 1505 may be made of materials such as an LCD (Liquid CRYSTAL DISPLAY) and an OLED (Organic Light-Emitting Diode).
The camera assembly 1506 is used to capture images or video. Optionally, the camera assembly 1506 includes a front camera and a rear camera. Typically, the front camera is disposed on a front panel of the computer device and the rear camera is disposed on a rear surface of the computer device. In some embodiments, the at least two rear cameras are any one of a main camera, a depth camera, a wide-angle camera and a tele camera, so as to realize that the main camera and the depth camera are fused to realize a background blurring function, and the main camera and the wide-angle camera are fused to realize a panoramic shooting and Virtual Reality (VR) shooting function or other fusion shooting functions. In some embodiments, the camera assembly 1506 may also include a flash. The flash lamp can be a single-color temperature flash lamp or a double-color temperature flash lamp. The dual-color temperature flash lamp refers to a combination of a warm light flash lamp and a cold light flash lamp, and can be used for light compensation under different color temperatures.
The audio circuitry 1507 may include a microphone and a speaker. The microphone is used for collecting sound waves of users and the environment, converting the sound waves into electric signals, inputting the electric signals to the processor 1501 for processing, or inputting the electric signals to the radio frequency circuit 1504 for voice communication. The microphone may be provided in a plurality of different locations of the computer device 1500 for stereo acquisition or noise reduction purposes. The microphone may also be an array microphone or an omni-directional pickup microphone. The speaker is used to convert electrical signals from the processor 1501 or the radio frequency circuit 1504 into sound waves. The speaker may be a conventional thin film speaker or a piezoelectric ceramic speaker. When the speaker is a piezoelectric ceramic speaker, not only the electric signal can be converted into a sound wave audible to humans, but also the electric signal can be converted into a sound wave inaudible to humans for ranging and other purposes. In some embodiments, the audio circuit 1507 may also include a headphone jack.
The positioning component 1508 is for positioning a current geographic location of the computer device 1500 to enable navigation or LBS (Location Based Service, location-based services). The positioning component 1508 may be a positioning component based on the United states GPS (Global Positioning System ), the Beidou system of China, or the Galileo system of Russia.
The power supply 1509 is used to power the various components in the computer device 1500. The power supply 1509 may be an alternating current, a direct current, a disposable battery, or a rechargeable battery. When the power supply 1509 includes a rechargeable battery, the rechargeable battery may be a wired rechargeable battery or a wireless rechargeable battery. The wired rechargeable battery is a battery charged through a wired line, and the wireless rechargeable battery is a battery charged through a wireless coil. The rechargeable battery may also be used to support fast charge technology.
In some embodiments, computer device 1500 also includes one or more sensors 1510. The one or more sensors 1510 include, but are not limited to, an acceleration sensor 1511, a gyroscope sensor 1512, a pressure sensor 1513, a fingerprint sensor 1514, an optical sensor 1515, and a proximity sensor 1516.
The acceleration sensor 1511 may detect the magnitudes of accelerations on three coordinate axes of the coordinate system established with the computer device 1500. For example, the acceleration sensor 1511 may be used to detect components of gravitational acceleration in three coordinate axes. The processor 1501 may control the display screen 1505 to display the user interface in a landscape view or a portrait view based on the gravitational acceleration signal acquired by the acceleration sensor 1511. The acceleration sensor 1511 may also be used for the acquisition of motion data of a game or user.
The gyro sensor 1512 may detect a body direction and a rotation angle of the computer apparatus 1500, and the gyro sensor 1512 may collect 3D actions of the user on the computer apparatus 1500 in cooperation with the acceleration sensor 1511. The processor 1501 can implement functions such as motion sensing (e.g., changing a UI according to a tilting operation of a user), image stabilization at photographing, game control, and inertial navigation based on data collected by the gyro sensor 1512.
Pressure sensor 1513 may be disposed on a side frame of computer device 1500 and/or under display screen 1505. When the pressure sensor 1513 is disposed on the side frame of the computer apparatus 1500, a grip signal of the user on the computer apparatus 1500 may be detected, and the processor 1501 performs a left-right hand recognition or a shortcut operation according to the grip signal collected by the pressure sensor 1513. When the pressure sensor 1513 is disposed at the lower layer of the display screen 1505, the processor 1501 realizes control of the operability control on the UI interface according to the pressure operation of the user on the display screen 1505. The operability controls include at least one of a button control, a scroll bar control, an icon control, and a menu control.
The fingerprint sensor 1514 is used to collect a fingerprint of a user, and the processor 1501 identifies the identity of the user based on the fingerprint collected by the fingerprint sensor 1514, or the fingerprint sensor 1514 identifies the identity of the user based on the collected fingerprint. Upon recognizing that the user's identity is a trusted identity, the processor 1501 authorizes the user to perform relevant sensitive operations including unlocking the screen, viewing encrypted information, downloading software, paying for and changing settings, etc. The fingerprint sensor 1514 may be disposed on the front, back, or side of the computer device 1500. When a physical key or vendor Logo is provided on the computer device 1500, the fingerprint sensor 1514 may be integrated with the physical key or vendor Logo.
The optical sensor 1515 is used to collect the ambient light intensity. In one embodiment, processor 1501 may control the display brightness of display screen 1505 based on the intensity of ambient light collected by optical sensor 1515. Specifically, the display luminance of the display screen 1505 is turned up when the ambient light intensity is high, and the display luminance of the display screen 1505 is turned down when the ambient light intensity is low. In another embodiment, the processor 1501 may also dynamically adjust the shooting parameters of the camera assembly 1506 based on the ambient light intensity collected by the optical sensor 1515.
A proximity sensor 1516, also referred to as a distance sensor, is typically provided on the front panel of the computer device 1500. The proximity sensor 1516 is used to capture the distance between the user and the front of the computer device 1500. In one embodiment, the processor 1501 controls the display screen 1505 to switch from the on-screen state to the off-screen state when the proximity sensor 1516 detects a gradual decrease in the distance between the user and the front of the computer device 1500, and the processor 1501 controls the display screen 1505 to switch from the off-screen state to the on-screen state when the proximity sensor 1516 detects a gradual increase in the distance between the user and the front of the computer device 1500.
Those skilled in the art will appreciate that the architecture shown in fig. 15 is not limiting as to the computer device 1500, and may include more or fewer components than shown, or may combine certain components, or employ a different arrangement of components.
The embodiment of the application also provides a computer readable storage medium, and at least one program code is stored in the computer readable storage medium, and the at least one program code is loaded and executed by a processor to realize the operation executed in the seismic exploration information acquisition method in the embodiment of the application.
Embodiments of the present application also provide a computer program product or computer program comprising computer program code stored in a computer readable storage medium. The processor of the computer device reads the computer program code from the computer readable storage medium, and the processor executes the computer program code to cause the computer device to perform the operations performed by the seismic survey information acquisition method described above.
In some embodiments, a computer program according to an embodiment of the present application may be deployed to be executed on one computer device or on multiple computer devices located at one site or on multiple computer devices distributed across multiple sites and interconnected by a communication network, where the multiple computer devices distributed across multiple sites and interconnected by a communication network may constitute a blockchain system.
The foregoing description is only for the convenience of those skilled in the art to understand the technical solution of the present application, and is not intended to limit the present application. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the protection scope of the present application.