Method, apparatus, and computer device for acquiring seismic exploration information

Abstract

The present application discloses a method, apparatus, and computer device for acquiring seismic exploration information, belonging to the technical field of geophysical exploration. An embodiment of the present application provides a method for acquiring seismic exploration information. Due to the complexity of the underground structure and lithology, the vibration conditions of the mass points at different acquisition times are complex and variable. This method first determines the polarization angle and polarization parameters corresponding to each acquisition time within a preset acquisition time range, and determines the vibration conditions of the mass points at each acquisition time according to the polarization angle and polarization parameters corresponding to each acquisition time. In this way, the true vibration conditions of the mass points underground can be reflected. Therefore, the seismic exploration information can be accurately acquired according to the vibration conditions of the mass points at each acquisition time, thereby improving the accuracy of the acquired seismic exploration information.

Description

Technical Field

[0001] This application relates to the technical field of geophysical exploration. In particular, it relates to a method, device, and computer device for acquiring seismic exploration information. Background Technique

[0002] Currently, during seismic exploration, seismic waves are first generated by a shot point, and then a three-component geophone is used to record the vibration of particles in a three-dimensional space caused by the seismic waves, obtaining three-component seismic data. The three-component seismic data includes seismic data of the X horizontal component, seismic data of the Y horizontal component, and seismic data of the Z vertical component, and the X horizontal component, Y horizontal component, and Z horizontal component are perpendicular to each other in pairs. Then, the three-component seismic data is processed, and seismic exploration information such as underground structural morphology, lithology, and oil and gas distribution is obtained through the processed seismic data.

[0003] In the related art, first, assuming that the formation is a horizontal layered formation and the lithology is an isotropic medium, the radial component R from the shot point to the three-component geophone and the tangential component T perpendicular to the radial component R are determined according to the azimuth of the shot point and the azimuth of the three-component geophone. Then, the seismic data corresponding to the X component and the seismic data corresponding to the Y component are respectively rotationally synthesized to obtain the seismic data corresponding to the radial component R and the seismic data corresponding to the tangential component T. Seismic exploration information is obtained based on the seismic data corresponding to the radial component R, the seismic data corresponding to the tangential component T, and the seismic data corresponding to the Z component.

[0004] However, in the method of the related art, under the assumption that the formation is a horizontal layered formation and the lithology is an isotropic medium, the vibration condition of particles at any moment is the same. Due to the complexity of the underground structural morphology and lithology, the vibration condition of particles caused by seismic waves at the three-component geophone is complex and variable. Therefore, this method cannot reflect the true vibration condition of particles underground, resulting in poor accuracy of the obtained seismic exploration information. Summary of the Invention

[0005] Embodiments of this application provide a method, device, and computer device for acquiring seismic exploration information, which can improve the accuracy of acquiring seismic exploration information. The specific technical solutions are as follows:

[0006] On the one hand, embodiments of this application provide a method for acquiring seismic exploration information, and the method includes:

[0007] Obtain multiple seismic data collected by multiple geophones in a target work area. The seismic data is collected by the geophones within a preset acquisition 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. The first horizontal component, the second horizontal component, and the first vertical component are perpendicular to each other in pairs.

[0008] For each geophone, based on the first seismic data and the second seismic data, determine the complex parameter of particle vibration, the linear parameter of particle vibration, and the cross-correlation parameter at the target acquisition moment. The complex parameter of particle vibration is used to reflect the complexity of particle vibration, the linear parameter of particle vibration is used to reflect the linearity of particle vibration, and the cross-correlation parameter is used to reflect the cross-correlation relationship between the first seismic data and the second seismic data at the target acquisition moment. The target acquisition moment is any acquisition moment within the preset acquisition time range.

[0009] Based on the complex parameter of particle vibration, the linear parameter of particle vibration, and the cross-correlation parameter at the target acquisition moment, determine the polarization parameter at the target acquisition moment and the polarization angle at the target acquisition moment.

[0010] Based on the first seismic data, the second seismic data, the polarization parameter at the target acquisition moment, and the polarization angle, determine the fourth seismic data. The fourth seismic data is the seismic data corresponding to the polarization angle at the target acquisition moment.

[0011] Based on the third seismic data corresponding to each geophone and the fourth seismic data at each acquisition moment, obtain seismic exploration information.

[0012] In a possible implementation manner, the determination process of determining the cross-correlation parameter at the target acquisition moment based on the first seismic data and the second seismic data includes:

[0013] Based on the first seismic data and the second seismic data, determine the frequency distribution range of the seismic wave excited by the shot point.

[0014] Based on the frequency distribution range, determine the time window parameter. The time window parameter is used to reflect the time window width corresponding to determining the cross-correlation parameter at the target acquisition moment.

[0015] Based on the time window parameter and the target acquisition moment, determine the integration time range.

[0016] Based on the integration time range, the first seismic data, and the second seismic data, determine the cross-correlation parameter at the target acquisition moment.

[0017] In another possible implementation, determining the cross-correlation parameter of the target acquisition moment based on the integration time range, the first seismic data, and the second seismic data includes:

[0018] Determine the square root of the energy of the first seismic data within the integration time range to obtain a first amplitude;

[0019] Determine the square root of the energy of the second seismic data within the integration time range to obtain a second amplitude;

[0020] Determine the product of the first amplitude and the second amplitude to obtain a third amplitude;

[0021] Determine the product of the first seismic data and the second seismic data within the integration time range to obtain a first vibration vector, where the first vibration vector includes an amplitude and a vibration direction;

[0022] Determine 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.

[0023] In another possible implementation, determining the polarization parameter of the target acquisition moment and the polarization angle of the target acquisition moment based on the complex particle vibration parameter, the linear particle vibration parameter, and the cross-correlation parameter of the target acquisition moment includes:

[0024] Based on the complex particle vibration parameter, the linear particle vibration parameter, and the cross-correlation parameter of the target acquisition moment, determine the polarization parameter of the target acquisition moment, where the polarization parameter is used to reflect the linear degree of particle vibration;

[0025] When the polarization parameter is not less than a preset polarization parameter, determine the polarization angle of the target acquisition moment based on the first seismic data and the second seismic data;

[0026] When the polarization parameter is less than the preset polarization parameter, determine the polarization angle of the target acquisition moment based on the first seismic data, the second seismic data, and the integration time range, where the integration time range is obtained when determining the cross-correlation parameter.

[0027] In another possible implementation, determining the polarization angle of the target acquisition moment based on the first seismic data and the second seismic data includes:

[0028] Determine the ratio of the second seismic data at the target acquisition moment to the first seismic data at the target acquisition moment to obtain a first ratio;

[0029] Determine the arctangent value of the first ratio to obtain the polarization angle at the target acquisition moment.

[0030] In another possible implementation, determining the polarization angle at the target acquisition moment based on the first seismic data, the second seismic data, and the integration time range includes:

[0031] Determine the square root of the sum of energies of the first seismic data and the second seismic data within the integration time range to obtain the fourth amplitude;

[0032] For each preset polarization angle within the preset polarization angle range, determine the product of the first seismic data at the target integration moment and the cosine value of the preset polarization angle to obtain the fifth seismic data, where the target integration moment is any integration moment within the integration time range;

[0033] Determine the product of the second seismic data at the target integration moment and the sine value of the preset polarization angle to obtain the sixth seismic data;

[0034] Determine the sum value of the fifth seismic data and the sixth seismic data within the integration time range to obtain the fifth amplitude;

[0035] Determine the difference between the fourth amplitude and the fifth amplitude to obtain the sixth amplitude corresponding to the preset polarization angle;

[0036] Take the preset polarization angle corresponding to the minimum sixth amplitude as the polarization angle at the target acquisition moment.

[0037] In another possible implementation, determining the polarization parameter at the target acquisition moment based on the complex particle vibration parameter, the linear particle vibration parameter, and the cross-correlation parameter at the target acquisition moment includes:

[0038] Determine the difference between the cross-correlation parameter at the target acquisition moment and the linear particle vibration parameter to obtain the first difference;

[0039] Determine the product of the complex particle vibration parameter and the first difference to obtain the first product;

[0040] Determine the exponential value with the natural constant as the base and the negative of the first product as the exponent;

[0041] Determine the reciprocal after adding 1 to the exponential value, and take the reciprocal as the polarization parameter.

[0042] In another possible implementation manner, determining, based on the first seismic data, the second seismic data, the polarization parameter and the polarization angle at the target acquisition moment, the third seismic data corresponding to the polarization angle at the target acquisition moment includes:

[0043] Determining the product of the first seismic data at the target acquisition moment and the cosine value of the polarization angle at the target acquisition moment to obtain the seventh seismic data;

[0044] Determining the product of the second seismic data at the target acquisition moment and the sine value of the polarization angle at the target acquisition moment to obtain the eighth seismic data;

[0045] Determining the product of the sum of the seventh seismic data and the eighth seismic data and the polarization parameter at the target acquisition moment to obtain the third seismic data corresponding to the polarization angle at the target acquisition moment.

[0046] On the other hand, an embodiment of the present application provides a device for acquiring seismic exploration information, and the device includes:

[0047] A first acquisition module, configured to acquire a plurality of seismic data collected by a plurality of geophones in a target work area, where the seismic data is collected by the geophones within a preset acquisition time range, and the seismic data includes the first seismic data of the first horizontal component, the second seismic data of the second horizontal component, and the third seismic data of the first vertical component, and the first horizontal component, the second horizontal component, and the first vertical component are perpendicular to each other in pairs;

[0048] A first determination module, configured to, for each geophone, determine a particle vibration complexity parameter, a particle vibration linear parameter, and a cross-correlation parameter at the target acquisition moment based on the first seismic data and the second seismic data, where the particle vibration complexity parameter is used to reflect the complexity degree of particle vibration, the particle vibration linear parameter is used to reflect the linear degree of particle vibration, the cross-correlation parameter is used to reflect the cross-correlation relationship between the first seismic data and the second seismic data at the target acquisition moment, and the target acquisition moment is any acquisition moment within the preset acquisition time range;

[0049] A second determination module, configured to determine the polarization parameter at the target acquisition moment and the polarization angle at the target acquisition moment based on the particle vibration complexity parameter, the particle vibration linear parameter, and the cross-correlation parameter at the target acquisition moment;

[0050] A third determination module, 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 moment, where the fourth seismic data is the seismic data corresponding to the polarization angle at the target acquisition moment;

[0051] A second acquisition module, configured to acquire seismic exploration information based on third seismic data corresponding to each geophone and fourth seismic data at each acquisition moment.

[0052] In a possible implementation manner, the first determination module is configured to: determine a frequency distribution range of seismic waves excited by a shot point 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 determining a cross-correlation parameter at the target acquisition moment; determine an integration time range based on the time window parameter and the target acquisition moment; and determine the cross-correlation parameter at the target acquisition moment based on the integration time range, the first seismic data, and the second seismic data.

[0053] In another possible implementation manner, the first determination module is configured to: determine a square root of energy of the first seismic data within the integration time range to obtain a first amplitude; determine a square root of energy of the second seismic data within 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 within 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 the cross-correlation parameter at the target acquisition moment.

[0054] In another possible implementation manner, the second determination module is configured to: determine a polarization parameter at the target acquisition moment based on the particle vibration complex parameter, the particle vibration linear parameter, and the cross-correlation parameter at the target acquisition moment, where the polarization parameter is used to reflect a linear degree of particle vibration; in a case where the polarization parameter is not less than a preset polarization parameter, determine a polarization angle at the target acquisition moment based on the first seismic data and the second seismic data; and in a case where the polarization parameter is less than the preset polarization parameter, determine the polarization angle at the target acquisition moment based on the first seismic data, the second seismic data, and the integration time range, where the integration time range is obtained when determining the cross-correlation parameter.

[0055] In another possible implementation manner, the second determination module is configured to: determine a ratio of the second seismic data at the target acquisition moment to the first seismic data at the target acquisition moment to obtain a first ratio; and determine an arctangent value of the first ratio to obtain the polarization angle at the target acquisition moment.

[0056] In another possible implementation, the second determination module is configured to determine the square root of the sum of energies of the first seismic data and the second seismic data within the integration time range to obtain a fourth amplitude; for each preset polarization angle within a preset polarization angle range, determine the product of the first seismic data at the target integration moment and the cosine value of the preset polarization angle to obtain a fifth seismic data, where the target integration moment is any integration moment within the integration time range; determine the product of the second seismic data at the target integration moment and the sine value of the preset polarization angle to obtain a sixth seismic data; determine the sum value of the fifth seismic data and the sixth seismic data within the integration time range to obtain a fifth amplitude; determine the difference between the fourth amplitude and the fifth amplitude to obtain a sixth amplitude corresponding to the preset polarization angle; and use the preset polarization angle corresponding to the smallest sixth amplitude as the polarization angle at the target acquisition moment.

[0057] In another possible implementation, the second determination module is configured to determine the difference between the cross-correlation parameter at the target acquisition moment and the linear parameter of particle vibration to obtain a first difference; determine the product of the complex parameter of particle vibration and the first difference to obtain a first product; determine the exponential value with the natural constant as the base and the negative of the first product as the exponent; and determine the reciprocal after adding 1 to the exponential value, and use the reciprocal as the polarization parameter.

[0058] In another possible implementation, the third determination module is configured to determine the product of the first seismic data at the target acquisition moment and the cosine value of the polarization angle at the target acquisition moment to obtain a seventh seismic data; determine the product of the second seismic data at the target acquisition moment and the sine value of the polarization angle at the target acquisition moment to obtain an eighth seismic data; and determine the product of the sum value of the seventh seismic data and the eighth seismic data and the polarization parameter at the target acquisition moment to obtain the third seismic data corresponding to the polarization angle at the target acquisition moment.

[0059] On the other hand, a computer device is provided. The computer device includes a processor and a memory. At least one program code is stored in the memory and is loaded and executed by the processor to implement the operations performed in the method for obtaining seismic exploration information in the embodiments of the present application.

[0060] On the other hand, an embodiment of the present application provides a computer-readable storage medium. At least one program code is stored in the computer-readable storage medium and is loaded and executed by a processor to implement the operations performed in the method for obtaining seismic exploration information in the embodiments of the present application.

[0061] On the other hand, an embodiment of the present application provides a computer program product or a computer program. The computer program product or the computer program includes computer program code, and the computer program code is stored in a computer-readable storage medium. A processor of a computer device reads the computer program code from the computer-readable storage medium, and the processor executes the computer program code to implement the operations performed in the method for obtaining seismic exploration information in the embodiments of the present application.

[0062] The beneficial effects brought by the technical solution provided by the embodiment of the present application are as follows:

[0063] An embodiment of the present application provides a method for obtaining seismic exploration information. Due to the complexity of the underground structure and lithology, the vibration conditions of the particle at different acquisition times are complex and variable. This method first determines the polarization angle and polarization parameters corresponding to each acquisition time within a preset acquisition time range, and determines the vibration condition of the particle at each acquisition time according to the polarization angle and polarization parameters corresponding to each acquisition time. In this way, the true vibration condition of the particle underground can be reflected. Therefore, the seismic exploration information can be accurately obtained according to the vibration condition of the particle at each acquisition time, thereby improving the accuracy of the obtained seismic exploration information. Description of the Drawings

[0064] Figure 1 is a schematic diagram of the vibration trajectories of a particle in the X horizontal component, Y horizontal component, and Z vertical component at different times provided by an embodiment of the present application;

[0065] Figure 2 is a schematic diagram of the vibration trajectories of a particle in the X horizontal component and Y horizontal component caused by random noise provided by an embodiment of the present application;

[0066] Figure 3 is a schematic diagram of the vibration trajectories of a particle in the X horizontal component and Y horizontal component caused by a seismic longitudinal wave provided by an embodiment of the present application;

[0067] Figure 4 is a schematic diagram of the vibration trajectories of a particle in the X horizontal component and Z vertical component caused by a Rayleigh surface wave provided by an embodiment of the present application;

[0068] Figure 5 is a schematic diagram of the projection and angle of the particle vibration on the XOY plane and XOZ plane provided by an embodiment of the present application;

[0069] Figure 6 is a flowchart of a method for obtaining seismic exploration information provided by an embodiment of the present application;

[0070] Figure 7Schematic diagram of the first seismic data of the theoretical synthetic X horizontal component provided by an embodiment of the present application;

[0071] Figure 8 Schematic diagram of the second seismic data of the theoretical synthetic Y horizontal component provided by an embodiment of the present application;

[0072] Figure 9 Is a provided by an embodiment of the present application Figure 7 The first seismic data in Figure 8 Schematic diagram of synthesizing the second seismic data in into the third seismic data;

[0073] Figure 10 Schematic diagram of the synthesized third seismic data provided by an embodiment of the present application;

[0074] Figure 11 Schematic diagram of the third seismic data corresponding to different polarization angles at different acquisition times provided by an embodiment of the present application;

[0075] Figure 12 Schematic diagram of the polarization angle determined from the actually acquired three-component seismic data provided by an embodiment of the present application;

[0076] Figure 13 Schematic diagram of the third seismic data determined from the actually acquired three-component seismic data provided by an embodiment of the present application;

[0077] Figure 14 Schematic diagram of the structure of an acquisition device for seismic exploration information provided by an embodiment of the present application;

[0078] Figure 15 Structure block diagram of a computer device provided by an embodiment of the present application. Specific embodiments

[0079] To make the technical solutions and advantages of the present application clearer, the embodiments of the present application will be further described in detail below.

[0080] For the convenience of elaboration, the particle vibration situation caused by seismic waves will be described here first.

[0081] The seismic wave is generated by the shot point. When the seismic wave propagates to the location of the three-component geophone, it causes the particles to vibrate. When the particles vibrate, they leave the equilibrium position and complete the vibration in three-dimensional space. After the seismic wave passes through the location of the three-component geophone, the particles return to their original equilibrium position. The three-component geophone can record the vibration of the particles caused by the seismic wave in three-dimensional space, comprehensively reflecting the propagation law and characteristics of the seismic wave in the underground medium. The three-component seismic data collected by the three-component geophone can accurately reconstruct the three-dimensional vibration trajectory of the particles. These vibration trajectories are a set of points arranged in time series in three-dimensional space and can be expressed as a function of time. See Figure 1 , Figure 1 which respectively represent the vibration trajectories of the particles on the X horizontal component, Y horizontal component and Z vertical component at different times. Connecting the vibration trajectory points in chronological order shows a curve in three-dimensional space. During different time periods, due to the different types of seismic waves reaching the three-component geophone, the vibration trajectories of the particles are also different. Among them, the main types of seismic waves include P-waves, S-waves, converted waves, Rayleigh surface waves, Love surface waves, random noise, etc. See Figures 2 to 4 , Figure 2 are the vibration trajectories of the particles on the X horizontal component and Y horizontal component caused by random noise, Figure 3 are the vibration trajectories of the particles on the X horizontal component and Y horizontal component caused by P-waves, Figure 4 are the vibration trajectories of the particles on the X horizontal component and Z vertical component caused by Rayleigh surface waves. It can be seen from Figures 2 to 4 that the vibration trajectories of the particles caused by different types of seismic waves are different. The two main characteristics of the vibration trajectory of the particles are the vibration amplitude and the vibration direction. Different types of seismic waves have different polarization characteristics. For example, linear, elliptical and three-dimensional. The polarization characteristics of the particle vibration not only depend on the type of seismic wave, but also are affected by the type of the seismic source and the complexity of the underground medium.

[0082] Although the polarization characteristics of different types of seismic waves are different, the vibration of the particles at the location of the three-component geophone caused by the seismic wave can be determined by a vibration vector whose magnitude and direction both change. The vibration trajectory of the particles can also be represented by different parameters related to the selected coordinate system. For example, in the Cartesian rectangular coordinate system, the projections of the vibration vector A(T) on the X, Y and Z coordinate axes can be respectively denoted as A x (T), A y (T) and A z (T), and the following formula holds:

[0083] A x (T)=|A(T)|sinφ(T)cosω(T)

[0084] Ay A(T) = |A(T)| sin φ(T) sin ω(T)

[0085] A z A(T) = |A(T)| cos φ(T)

[0086] Where, |A(T)| represents the absolute value of the displacement of the particle polarization at the T acquisition moment, ω(T) represents the polarization angle between the particle polarization and the X-axis on the horizontal plane at the T acquisition moment, φ(T) represents the angle between the particle polarization and the vertical plane at the T acquisition moment, A x A(T) represents the displacement of the particle vibration recorded by the sensor of the x horizontal component, A y A(T) represents the displacement of the particle vibration recorded by the sensor of the y horizontal component, A z A(T) represents the displacement of the particle vibration recorded by the sensor of the z vertical component. See Figure 5 , Figure 5 is a schematic diagram of the projection and angle of the particle vibration on the XOY plane and the XOZ plane. Where, A o A(T) is the projection of A(T) on the XOY plane, A1(T) is the projection of A(T) on the XOZ plane, ω(T) is the angle between A o A(T) and the X-axis direction, φ(T) is the angle between A1(T) and the Z-axis direction.

[0087] When the seismic longitudinal wave recorded by the three-component geophone does not interfere with other waves, the particle vibration it causes is linearly polarized. Linear polarization means that the direction of the vibration vector A(T) remains unchanged, while its magnitude can change, that is, the particle vibrates in a straight-line trajectory near the equilibrium position. In a homogeneous isotropic medium, the seismic shear wave is linearly polarized in the plane tangent to the seismic longitudinal wavefront. The Rayleigh surface wave is elliptically polarized, and the polarization plane is vertical. The amplitudes of the horizontal and vertical components of the Rayleigh surface wave vary differently with depth. Different from the Rayleigh surface wave, the Love surface wave is linearly polarized on the horizontal plane perpendicular to the wave propagation direction. Random noise is characterized by the instability of polarization characteristics. Continue to see Figure 2 , which means that the phase difference between its different components in the three-component geophone is random. In a relatively short acquisition time range, different components can be in phase and the recorded shapes are similar, but in a relatively large acquisition time range, the cross-correlation parameter of the seismic data recorded by the sensors of different components of the same random noise in the three-component geophone is approximately zero.

[0088] In order to improve the processing and interpretation efficiency and accuracy of three-component seismic data, the projection of the particle vibration at the position of the three-component geophone on the XOY horizontal plane is defined as the optimal receiving component, that is, the maximum vibration amplitude, and the seismic signal received at this position is the strongest.

[0089] Next, the solution of this application will be introduced in detail:

[0090] An embodiment of this application provides a method for obtaining seismic exploration information, which is executed by a computer device. Refer to Figure 6 , and this method includes:

[0091] Step 601: The computer device acquires a plurality of seismic data collected by a plurality of geophones in the target work area.

[0092] For each geophone, the seismic data is collected by the geophone within a preset acquisition 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. The first horizontal component, the second horizontal component, and the first vertical component are perpendicular to each other in pairs.

[0093] Refer to Figure 7 and Figure 8 , Figure 7 is the first seismic data of the theoretically synthesized X horizontal component, Figure 8 is the second seismic data of the theoretically synthesized Y horizontal component.

[0094] Step 602: For each geophone, the computer device determines a particle vibration complex parameter, a particle vibration linear parameter, and a cross-correlation parameter at the target acquisition moment based on the first seismic data and the second seismic data.

[0095] The particle vibration complex parameter is used to reflect the complexity of particle vibration, the particle vibration linear parameter is used to reflect the linearity of particle vibration, and the cross-correlation parameter is used to reflect the cross-correlation relationship between the first seismic data and the second seismic data at the target acquisition moment. The target acquisition moment is any acquisition moment within the preset acquisition time range.

[0096] Step 602 can be implemented through the following steps (1) to (2), including:

[0097] (1) The computer device determines the cross-correlation parameter at the target acquisition moment based on the first seismic data and the second seismic data.

[0098] Step (1) can be implemented through the following steps (1-1) to (1-4), including:

[0099] (1-1) The computer device determines the frequency distribution range of the seismic waves excited by the shot point based on the first seismic data and the second seismic data.

[0100] Among the seismic waves excited by the shot point, there are P-waves and S-waves. The computer device obtains the frequency distribution ranges of the P-waves and S-waves from the first seismic data and the second seismic data.

[0101] (1-2) The computer device determines the time window parameter based on this frequency distribution range.

[0102] This time window parameter is used to reflect the time window width corresponding to the cross-correlation parameter when determining the target acquisition moment. The computer device can determine the time window parameter based on the frequency distribution range of the seismic longitudinal wave, or can also determine the time window parameter based on the frequency distribution range of the seismic shear wave. Or, this time window parameter can also be half of the time window width, that is, the half time window.

[0103] Here, only the example of the computer device determining the time window parameter based on the frequency distribution range of the seismic shear wave is used for illustration. The computer device determines the dominant frequency of the seismic shear wave based on the frequency distribution range of the seismic shear wave, and takes the reciprocal of this dominant frequency as the time window parameter. Among them, the dominant frequency of the seismic shear wave refers to the frequency with the strongest energy of the seismic shear wave. For example, if the dominant frequency of the seismic shear wave is 10 Hz, then the time window parameter is 0.1. The computer device can also determine the time window parameter by other methods, and no specific limitation is made on this.

[0104] (1-3) The computer device determines the integration time range based on this time window parameter and this target acquisition moment.

[0105] The computer device takes the sum value of this time window parameter and this target acquisition moment as the upper limit value of the integration time range, and takes the difference value between this time window parameter and this target acquisition moment as the lower limit value of the integration time range, and obtains the integration time range.

[0106] (1-4) The computer device determines the cross-correlation parameter of this target acquisition moment based on this integration time range, the first seismic data, and the second seismic data.

[0107] Step (1-4) can be implemented through the following steps (1-4-1) to (1-4-5), including:

[0108] (1-4-1) The computer device determines the square root of the energy of the first seismic data within this integration time range to obtain the first amplitude.

[0109] The first amplitude can be expressed by the following formula: Among them, t0 represents any integration moment within the integration time range from t - n to t + n, t represents this target acquisition moment, n represents the time window parameter, A x (t0) represents the first seismic data at the t0 moment, x represents the first horizontal component, represents the energy of the first seismic data within this integration time range.

[0110] (1-4-2) The computer device determines the square root of the energy of the second seismic data within this integration time range to obtain the second amplitude.

[0111] The second amplitude can be expressed by the following formula: where A y (t0) represents the second seismic data at the integration time of t0, y represents the second horizontal component, represents the energy of the second seismic data within this integration time range.

[0112] (1 - 4 - 3) The computer device determines the product of the first amplitude and the second amplitude to obtain the third amplitude.

[0113] The third amplitude can be expressed by the following formula:

[0114] (1 - 4 - 4) The computer device determines the product of the first seismic data and the second seismic data within this integration time range to obtain the first vibration vector.

[0115] The first vibration vector includes an amplitude and a vibration direction, and the first vibration vector can be expressed by the following formula:

[0116] (1 - 4 - 5) The computer device determines the absolute value of the ratio of the first vibration vector to the third amplitude to obtain the cross - correlation parameter at this target acquisition time.

[0117] The cross - correlation parameter can be expressed by the following formula:

[0118] where c(t) represents the cross - correlation parameter at the acquisition time of t.

[0119] (2) The computer device determines the complex particle vibration parameter and the linear particle vibration parameter based on the first seismic data and the second seismic data.

[0120] Both the complex particle vibration parameter and the linear particle vibration parameter are constants. The computer device can determine the complex particle vibration parameter and the linear particle vibration parameter based on the linear particle vibration situation reflected by the first seismic data and the second seismic data. The more obvious and prominent the linear particle vibration is, the smaller the complex particle vibration parameter is, and the smaller the linear particle vibration parameter is.

[0121] The complex particle vibration parameter is used to control the shape of the polarization parameter. The complex particle vibration parameter can be expressed as d0, d0 ∈ [1, 10000]. In most cases, d0 can take the value of 10. The linear particle vibration parameter serves as the threshold value for controlling the enhancement or suppression of seismic waves with linear polarization characteristics. The linear particle vibration parameter can be expressed as c0, c0 ∈ [0, 1].

[0122] Step 603: The computer device determines the polarization parameter at the target acquisition moment based on the complex parameters of particle vibration, the linear parameters of particle vibration, and the cross-correlation parameter at the target acquisition moment.

[0123] Step 603 can be implemented through the following steps (1) to (4), including:

[0124] (1) The computer device determines the difference between the cross-correlation parameter at the target acquisition moment and the linear parameter of particle vibration to obtain a first difference.

[0125] The first difference can be expressed by the following formula: c(t) - c0, where c(t) represents the cross-correlation parameter at the t acquisition moment, and c0 represents the linear parameter of particle vibration.

[0126] (2) The computer device determines the product of the complex parameter of particle vibration and the first difference to obtain a first product.

[0127] The first product can be expressed by the following formula: d0(c(t) - c0), where d0 represents the complex parameter of particle vibration.

[0128] (3) The computer device determines the exponential value with the natural constant as the base and the negative of the first product as the exponent.

[0129] The exponential value can be expressed by the following formula: where e is the natural constant, and its value is approximately equal to 2.71828.

[0130] (4) The computer device determines the reciprocal after adding 1 to the exponential value, and takes the reciprocal as the polarization parameter at the target acquisition moment.

[0131] The polarization parameter is used to reflect the linear degree of particle vibration and can be expressed by the following formula: where γ(t) represents the polarization parameter at the target acquisition moment.

[0132] After the computer device obtains the polarization parameter at the target acquisition moment, it determines the magnitude relationship between the polarization parameter and the preset polarization parameter. When the polarization parameter is not less than the preset polarization parameter, step 604 is executed; when the polarization parameter is less than the preset polarization parameter, step 605 is executed.

[0133] The preset polarization parameter can be set and changed as needed, and in the embodiments of the present application, no specific limitation is made thereto. For example, the preset polarization parameter is 0.8 or 0.9.

[0134] It should be noted that γ(t) ∈ [0, 1]. When γ(t) = 1, A x (t) and Ay (t). Considering that oil and gas exploration utilizes the information of seismic longitudinal waves and seismic transverse waves to detect underground formations, and both have the characteristic of linear polarization. Therefore, γ(t) is an activation function with the absolute value of the cross-correlation function of A x (t) and A y (t) as the independent variable.

[0135] Step 604: When the polarization parameter at the target acquisition moment is not less than the preset polarization parameter, the computer device determines the polarization angle at the target acquisition moment based on the first seismic data and the second seismic data.

[0136] When the polarization parameter at the target acquisition moment is not less than the preset polarization parameter, it indicates that the particle vibration characteristics are typical linear polarization. The computer device determines the ratio of the second seismic data at the target acquisition moment to the first seismic data at the target acquisition moment, obtains a first ratio, and determines the arctangent value of the first ratio to obtain the polarization angle at the target acquisition moment.

[0137] When the polarization parameter at the target acquisition moment is not less than the preset polarization parameter, the polarization angle can be expressed by the following formula: where ω(t) represents the polarization angle at the acquisition moment t, arctan(.) represents the arctangent function, A y (t) represents the second seismic data at the acquisition moment t, and A x (t) represents the first seismic data at the acquisition moment t.

[0138] Step 605: When the polarization parameter at the target acquisition moment is less than the preset polarization parameter, the computer device determines the polarization angle at the target acquisition moment based on the first seismic data, the second seismic data, and the integration time range.

[0139] The integration time range is obtained when determining the cross-correlation parameter in Step 602. When the polarization parameter is less than the preset polarization parameter, it indicates that the particle vibration may be caused by the superposition of seismic waves of two or more types of linear polarization, or by elliptically polarized Rayleigh surface waves or random noise, etc. The polarization angle needs to be determined by the method of scanning analysis. Correspondingly, Step 605 is implemented through the following steps (1) to (6), including:

[0140] (1) The computer device determines the square root of the sum of the energies of the first seismic data and the second seismic data within the integration time range to obtain a fourth amplitude.

[0141] The fourth amplitude can be expressed by the following formula:

[0142] (2) For any preset polarization angle within the preset polarization angle range, the computer device determines the product of the first seismic data at the target integration moment and the cosine value of the preset polarization angle to obtain the fifth seismic data.

[0143] The target integration moment is any integration moment within the integration time range, and the preset polarization angle range is That is

[0144] The fifth seismic data can be expressed by the following formula: A x (t0)cosω(β), where t0 represents the target integration moment, cosω(β) represents the cosine value of the preset polarization angle, β represents the preset polarization angle, ω represents the polarization angle, and cosω(β) represents the cosine value when ω is β.

[0145] (3) The computer device determines the product of the second seismic data at the target integration moment and the sine value of the preset polarization angle to obtain the sixth seismic data.

[0146] The sixth seismic data can be expressed by the following formula: A y (t0)sinω(β), where sinω(β) represents the sine value when ω is β.

[0147] (4) The computer device determines the sum value of the fifth seismic data and the sixth seismic data within the integration time range to obtain the fifth amplitude.

[0148] The fifth amplitude can be expressed by the following formula:

[0149]

[0150] (5) The computer device determines the difference between the fourth amplitude and the fifth amplitude to obtain the sixth amplitude corresponding to the preset polarization angle.

[0151] The sixth amplitude can be expressed by the following formula:

[0152]

[0153] (6) The computer device takes the preset polarization angle corresponding to the smallest sixth amplitude as the polarization angle at the target acquisition moment.

[0154] The computer device can determine the sixth amplitude corresponding to each preset polarization angle within the preset polarization angle range through the above steps (1) to (5), and the computer device takes the preset polarization angle corresponding to the smallest sixth amplitude as the polarization angle at the target acquisition moment.

[0155] The polarization angle at the target acquisition moment can be expressed as:

[0156]

[0157] See Figure 9 , Figure 9 Based on the polarization angle and polarization parameter at the target acquisition moment, Figure 7 The first seismic data in Figure 8 And the second seismic data in Figure 9 Are combined into the third seismic data, and different colors represent different polarization angles.

[0158] Step 606: The computer device determines the fourth seismic data based on the first seismic data, the second seismic data, the polarization parameter and the polarization angle at the target acquisition moment.

[0159] The fourth seismic data is the seismic data corresponding to the polarization angle at the target acquisition moment. Correspondingly, step 606 can be implemented through the following steps (1) to (3), including:

[0160] (1) The computer device determines the product of the first seismic data at the target acquisition moment and the cosine value of the polarization angle at the target acquisition moment to obtain the seventh seismic data.

[0161] The first seismic data is the data corresponding to the preset acquisition time range, and the computer device obtains the first seismic data at the target acquisition moment.

[0162] The seventh seismic data can be expressed by the following formula: A x (t)cosω(t), where cosω(t) is the cosine value of the polarization angle at the acquisition moment t.

[0163] (2) The computer device determines the product of the second seismic data at the target acquisition moment and the sine value of the polarization angle at the target acquisition moment to obtain the eighth seismic data.

[0164] The second seismic data is the data corresponding to the preset acquisition time range, and the computer device obtains the second seismic data at the target acquisition moment.

[0165] The eighth seismic data can be expressed by the following formula: A y (t)sinω(t), where sinω(t) is the sine value of the polarization angle at the acquisition moment t.

[0166] (3) The computer device determines the product of the sum of the seventh seismic data and the eighth seismic data and the polarization parameter at the target acquisition moment to obtain the fourth seismic data.

[0167] The fourth seismic data can be expressed by the following formula:

[0168] A o (t) = γ(t)(A x(t)cosω(t) + A y (t)sinω(t)), where A o (t) is the fourth seismic data corresponding to the ω polarization angle at the t acquisition moment, and γ(t) is the polarization parameter at the t acquisition moment.

[0169] It should be noted that the method provided in the embodiments of the present application is mainly applied to the processing and interpretation of multi-wave and multi-component seismic data, which can shorten the production cycle of the processing and interpretation work of multi-wave and multi-component seismic data and improve the accuracy of the processing and interpretation results. This method analyzes the three-component seismic data to determine the polarization angle of the particle vibration at the three-component geophone and synthesizes the optimal receiving component to improve the processing and interpretation efficiency and accuracy of the three-component seismic data. Compared with the methods in the related art, this method fully considers the complexity of the underground seismic wave propagation, can improve the detection efficiency and accuracy of complex formations and complex media, and reduce the risks of oil and gas exploration and development.

[0170] See Figure 10 , Figure 10 For synthesizing the third seismic data obtained from the first seismic data in Figure 7 and the second seismic data in Figure 8 , compared with the first seismic data in Figure 7 and the second seismic data in Figure 8 , Figure 10 the third seismic data in

[0171] is more comprehensive and complete, and the third seismic data can more truly and comprehensively reflect the particle vibration situation caused by the seismic wave. Figure 11 , Figure 11 See Figure 11 for the third seismic data corresponding to different polarization angles at different acquisition moments, Figure 12 , Figure 12 which can intuitively express the projection of the particle vibration in the XOY plane, that is, the third seismic data and the polarization angle. See Figure 13 , Figure 13 for the polarization angle determined from the actually acquired three-component seismic data, and different colors represent different polarization angles. See Figure 13 , Figure 13 for the third seismic data determined from the actually acquired three-component seismic data. It can be seen from Figure 12 and Figure 13 that different types of seismic waves have different polarization angles.

[0172] Step 607: The computer device obtains seismic exploration information based on the third seismic data corresponding to each geophone and the fourth seismic data at each acquisition moment.

[0173] The third seismic data is the seismic data of the first vertical component. The computer device obtains seismic exploration information based on the third seismic data of each acquisition moment corresponding to each geophone and the fourth seismic data of each acquisition moment. Among them, for each acquisition moment corresponding to each geophone, the computer device can determine the vibration trajectory of the particle based on the fourth seismic data of this acquisition moment, and determine the vibration trajectory of the particle on the first vertical component based on the third seismic data of this acquisition moment. The seismic exploration information is obtained according to the vibration trajectory at the polarization angle and the vibration trajectory on the first vertical component.

[0174] In the embodiment of the present application, after the computer device obtains the fourth seismic data of each geophone at each acquisition moment, it can correspondingly store the fourth seismic data of each acquisition moment and the polarization angle of each acquisition moment to guide subsequent processing.

[0175] The embodiment of the present application provides a method for obtaining seismic exploration information. Due to the complexity of the underground structure and lithology, the vibration conditions of the particle at different acquisition moments are complex and variable. This method first determines the polarization angle and polarization parameters corresponding to each acquisition moment within the preset acquisition time range, and determines the vibration conditions of the particle at each acquisition moment according to the polarization angle and polarization parameters corresponding to each acquisition moment. This can reflect the real vibration conditions of the particle underground. Therefore, the seismic exploration information can be accurately obtained according to the vibration conditions of the particle at each acquisition moment, thereby improving the accuracy of the obtained seismic exploration information.

[0176] The embodiment of the present application provides a device for obtaining seismic exploration information. Refer to Figure 14 , the device includes:

[0177] The first acquisition module 1401 is used to acquire a plurality of seismic data collected by a plurality of geophones in the target work area. The seismic data is collected by the geophones within the preset acquisition time range, and the seismic data includes the first seismic data of the first horizontal component, the second seismic data of the second horizontal component, and the third seismic data of the first vertical component. The first horizontal component, the second horizontal component, and the first vertical component are perpendicular to each other in pairs;

[0178] The first determination module 1402 is used to, for each geophone, determine the particle vibration complex parameter, the particle vibration linear parameter, and the cross-correlation parameter of the target acquisition moment based on the first seismic data and the second seismic data. The particle vibration complex parameter is used to reflect the complexity of the particle vibration, the particle vibration linear parameter is used to reflect the linearity of the particle vibration, and the cross-correlation parameter is used to reflect the cross-correlation relationship between the first seismic data and the second seismic data at the target acquisition moment. The target acquisition moment is any acquisition moment within the preset acquisition time range;

[0179] The second determination module 1403 is configured to determine the polarization parameter and the polarization angle at the target acquisition moment based on the complex parameters of particle vibration, the linear parameters of particle vibration, and the cross-correlation parameter at the target acquisition moment;

[0180] The third determination module 1404 is configured to determine the fourth seismic data based on the first seismic data, the second seismic data, the polarization parameter and the polarization angle at the target acquisition moment, where the fourth seismic data is the seismic data corresponding to the polarization angle at the target acquisition moment;

[0181] The second acquisition module 1405 is configured to acquire seismic exploration information based on the third seismic data corresponding to each geophone and the fourth seismic data at each acquisition moment.

[0182] In a possible implementation manner, the first determination module 1402 is configured to determine the frequency distribution range of the seismic wave excited by the shot point based on the first seismic data and the second seismic data; determine the time window parameter based on the frequency distribution range, where the time window parameter is used to reflect the time window width corresponding to determining the cross-correlation parameter at the target acquisition moment; determine the integration time range based on the time window parameter and the target acquisition moment; determine the cross-correlation parameter at the target acquisition moment based on the integration time range, the first seismic data, and the second seismic data.

[0183] In another possible implementation manner, the first determination module 1402 is configured to determine the square root of the energy of the first seismic data within the integration time range to obtain the first amplitude; determine the square root of the energy of the second seismic data within the integration time range to obtain the second amplitude; determine the product of the first amplitude and the second amplitude to obtain the third amplitude; determine the product of the first seismic data and the second seismic data within the integration time range to obtain the first vibration vector, where the first vibration vector includes an amplitude and a vibration direction; determine the absolute value of the ratio of the first vibration vector to the third amplitude to obtain the cross-correlation parameter at the target acquisition moment.

[0184] In another possible implementation manner, the second determination module 1403 is configured to determine the polarization parameter at the target acquisition moment based on the complex parameters of particle vibration, the linear parameters of particle vibration, and the cross-correlation parameter at the target acquisition moment, where the polarization parameter is used to reflect the linear degree of particle vibration; in the case where the polarization parameter is not less than the preset polarization parameter, determine the polarization angle at the target acquisition moment based on the first seismic data and the second seismic data; in the case where the polarization parameter is less than the preset polarization parameter, determine the polarization angle at the target acquisition moment based on the first seismic data, the second seismic data, and the integration time range, where the integration time range is obtained when determining the cross-correlation parameter.

[0185] In another possible implementation manner, the second determination module 1403 is configured to determine a ratio of the second seismic data at the target acquisition moment to the first seismic data at the target acquisition moment, to obtain a first ratio; and determine an arctangent value of the first ratio, to obtain a polarization angle at the target acquisition moment.

[0186] In another possible implementation manner, the second determination module 1403 is configured to determine a square root of a sum of energies of the first seismic data and the second seismic data within an integration time range, to obtain a fourth amplitude; for each preset polarization angle within a preset polarization angle range, determine a product of the first seismic data at the target integration moment and a cosine value of the preset polarization angle, to obtain a fifth seismic data, where the target integration moment is any integration moment within the integration time range; determine a product of the second seismic data at the target integration moment 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 within 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 the preset polarization angle corresponding to the minimum sixth amplitude as the polarization angle at the target acquisition moment.

[0187] In another possible implementation manner, the second determination module 1403 is configured to determine a difference between a cross-correlation parameter and a particle vibration linear parameter at the target acquisition moment, to obtain a first difference; determine a product of a particle vibration complex parameter and the first difference, to obtain a first product; determine an exponential value with the natural constant as the base and the negative of the first product as the exponent; and determine a reciprocal after adding 1 to the exponential value, and use the reciprocal as a polarization parameter.

[0188] In another possible implementation manner, the third determination module 1404 is configured to determine a product of the first seismic data at the target acquisition moment and a cosine value of the polarization angle at the target acquisition moment, to obtain a seventh seismic data; determine a product of the second seismic data at the target acquisition moment and a sine value of the polarization angle at the target acquisition moment, to obtain an eighth seismic data; and determine a product of a sum value of the seventh seismic data and the eighth seismic data and the polarization parameter at the target acquisition moment, to obtain a third seismic data corresponding to the polarization angle at the target acquisition moment.

[0189] The embodiment of the present application provides a seismic exploration information acquisition device. Due to the complexity of the underground structure form and lithology, the vibration conditions of particles at different acquisition moments are complex and variable. The device first determines the polarization angle and polarization parameter corresponding to each acquisition moment within a preset acquisition time range, and determines the vibration conditions of particles at each acquisition moment according to the polarization angle and polarization parameter corresponding to each acquisition moment, so as to reflect the real vibration conditions of particles underground. Therefore, the seismic exploration information can be accurately acquired according to the vibration conditions of particles at each acquisition moment, thereby improving the accuracy of the acquired seismic exploration information.

[0190] Figure 15 The structural block diagram of a computer device 1500 provided by an exemplary embodiment of the present application is shown. The computer device 1500 may be a portable mobile computer device, such as: a smart phone, a tablet computer, an MP3 player (Moving Picture Experts Group Audio Layer III), an MP4 (Moving Picture Experts Group Audio Layer IV) player, a notebook computer or a desktop computer. The computer device 1500 may also be referred to by other names such as a user device, a portable computer device, a laptop computer device, a desktop computer device, etc.

[0191] Generally, the computer device 1500 includes: a processor 1501 and a memory 1502.

[0192] The processor 1501 may include one or more processing cores, such as a quad-core processor, an octa-core processor, etc. 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 and a co-processor. The main processor is a processor for processing data in the wake state, also known as the CPU (Central Processing Unit); the co-processor is a low-power processor for processing data in the standby state. In some embodiments, the processor 1501 may be integrated with a GPU (Graphics Processing Unit), and the GPU is responsible for rendering and drawing the content to be displayed on the display screen. In some embodiments, the processor 1501 may further include an AI (Artificial Intelligence) processor, and the AI processor is used to process computing operations related to machine learning.

[0193] The memory 1502 may include one or more computer-readable storage media, which may be non-transitory. The memory 1502 may further include high-speed random access memory, as well as non-volatile memory, such as one or more magnetic disk storage devices and flash memory storage devices. In some embodiments, the non-transitory computer-readable storage media in the memory 1502 is used to store at least one instruction for being executed by the processor 1501 to implement the method for obtaining seismic exploration information provided in the method embodiments of the present application.

[0194] In some embodiments, the computer device 1500 may further optionally include: a peripheral device interface 1503 and at least one peripheral device. The processor 1501, the memory 1502, and the peripheral device interface 1503 may be connected through a bus or signal lines. Each peripheral device may be connected to the peripheral device interface 1503 through a bus, signal lines, or a circuit board. Specifically, the peripheral device includes at least one of a radio frequency circuit 1504, a display screen 1505, a camera assembly 1506, an audio circuit 1507, a positioning assembly 1508, and a power supply 1509.

[0195] The peripheral device interface 1503 may be used to connect at least one peripheral device related to I / O (Input / Output) to the processor 1501 and the memory 1502. In some embodiments, the processor 1501, the memory 1502, and the peripheral device interface 1503 are integrated on the same chip or circuit board; in some other embodiments, any one or two of the processor 1501, the memory 1502, and the peripheral device interface 1503 may be implemented on a separate chip or circuit board, and the present embodiment does not limit this.

[0196] The radio frequency circuit 1504 is used 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 through electromagnetic signals. The radio frequency circuit 1504 converts an electrical signal into an electromagnetic signal for transmission, or converts the received electromagnetic signal into an electrical signal. Optionally, the radio frequency circuit 1504 includes: an antenna system, an RF transceiver, one or more amplifiers, a tuner, an oscillator, a digital signal processor, a codec chipset, a subscriber identity module card, and so on. The radio frequency circuit 1504 can communicate with other computer devices through at least one wireless communication protocol. The wireless communication protocol includes but is not limited to: the World Wide Web, a metropolitan area network, an intranet, generations of mobile communication networks (2G, 3G, 4G, and 5G), a wireless local area network, and / or a WiFi (Wireless Fidelity) network. In some embodiments, the radio frequency circuit 1504 may further include a circuit related to NFC (Near Field Communication), which is not limited in this application.

[0197] The display screen 1505 is used to display a UI (User Interface). The UI may include graphics, text, icons, videos, and any combination thereof. When the display screen 1505 is a touch display screen, the display screen 1505 also has the ability to collect touch signals on or above the surface of the display screen 1505. The touch signal can be input to the processor 1501 as a control signal for processing. At this time, the display screen 1505 can also be used to provide virtual buttons and / or a virtual keyboard, also known as soft buttons and / or a soft keyboard. In some embodiments, there may be one display screen 1505, which is disposed on the front panel of the computer device 1500; in other embodiments, there may be at least two display screens 1505, which are respectively disposed on different surfaces of the computer device 1500 or are in a foldable design; in other embodiments, the display screen 1505 may be a flexible display screen, which is disposed on a curved surface or a foldable surface of the computer device 1500. Even, the display screen 1505 can be set to an irregular non-rectangular shape, that is, a special-shaped screen. The display screen 1505 can be prepared using materials such as LCD (Liquid Crystal Display) and OLED (Organic Light-Emitting Diode).

[0198] The camera component 1506 is used to collect images or videos. Optionally, the camera component 1506 includes a front camera and a rear camera. Generally, the front camera is set on the front panel of the computer device, and the rear camera is set on the back of the computer device. In some embodiments, there are at least two rear cameras, which can be any one of a main camera, a depth camera, a wide-angle camera, and a telephoto camera, so as to realize the function of background blurring by fusing the main camera and the depth camera, panoramic shooting by fusing the main camera and the wide-angle camera, and VR (Virtual Reality) shooting function or other fusion shooting functions. In some embodiments, the camera component 1506 may further include a flash. The flash can be a single-color temperature flash or a dual-color temperature flash. The dual-color temperature flash refers to the combination of a warm light flash and a cold light flash, which can be used for light compensation under different color temperatures.

[0199] The audio circuit 1507 may include a microphone and a speaker. The microphone is used to collect sound waves of the user and the environment, and convert the sound waves into electrical signals and input them to the processor 1501 for processing, or input them to the radio frequency circuit 1504 to achieve voice communication. For the purpose of stereo collection or noise reduction, there may be multiple microphones, which are respectively set at different parts of the computer device 1500. The microphone can also be an array microphone or an omnidirectional collection microphone. The speaker is used to convert the electrical signal from the processor 1501 or the radio frequency circuit 1504 into sound waves. The speaker can be a traditional thin film speaker or a piezoelectric ceramic speaker. When the speaker is a piezoelectric ceramic speaker, it can not only convert the electrical signal into sound waves audible to humans, but also convert the electrical signal into sound waves inaudible to humans for uses such as ranging. In some embodiments, the audio circuit 1507 may further include a headphone jack.

[0200] The positioning component 1508 is used to locate the current geographical location of the computer device 1500 to achieve navigation or LBS (Location Based Service). The positioning component 1508 can be a positioning component based on the US GPS (Global Positioning System), China's Beidou system, or Russia's Galileo system.

[0201] The power supply 1509 is used to supply power to each component in the computer device 1500. The power supply 1509 can be alternating current, direct current, a disposable battery, or a rechargeable battery. When the power supply 1509 includes a rechargeable battery, the rechargeable battery can 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 can also be used to support fast charging technology.

[0202] In some embodiments, the computer device 1500 further 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.

[0203] The acceleration sensor 1511 can detect the magnitude of acceleration on the three coordinate axes of the coordinate system established with the computer device 1500. For example, the acceleration sensor 1511 can be used to detect the components of the gravitational acceleration on the three coordinate axes. The processor 1501 can control the display screen 1505 to display the user interface in a landscape view or a portrait view according to the gravitational acceleration signal collected by the acceleration sensor 1511. The acceleration sensor 1511 can also be used for collecting game or user movement data.

[0204] The gyroscope sensor 1512 can detect the body direction and rotation angle of the computer device 1500. The gyroscope sensor 1512 can cooperate with the acceleration sensor 1511 to collect the 3D actions of the user on the computer device 1500. Based on the data collected by the gyroscope sensor 1512, the processor 1501 can implement the following functions: motion sensing (such as changing the UI according to the user's tilting operation), image stabilization during shooting, game control, and inertial navigation.

[0205] The pressure sensor 1513 can be disposed on the side frame of the computer device 1500 and / or the lower layer of the display screen 1505. When the pressure sensor 1513 is disposed on the side frame of the computer device 1500, it can detect the holding signal of the user on the computer device 1500, and the processor 1501 can perform left - hand / right - hand recognition or quick operation according to the holding signal collected by the pressure sensor 1513. When the pressure sensor 1513 is disposed on the lower layer of the display screen 1505, the processor 1501 can control the operable controls on the UI interface according to the pressure operation of the user on the display screen 1505. The operable controls include at least one of button controls, scroll bar controls, icon controls, and menu controls.

[0206] The fingerprint sensor 1514 is used to collect the user's fingerprint. The processor 1501 identifies the user's identity based on the fingerprint collected by the fingerprint sensor 1514, or the fingerprint sensor 1514 identifies the user's identity based on the collected fingerprint. When the user's identity is identified as a trusted identity, the processor 1501 authorizes the user to perform relevant sensitive operations, which include unlocking the screen, viewing encrypted information, downloading software, making payments, and changing settings, etc. The fingerprint sensor 1514 can be set on the front, back, or side of the computer device 1500. When there are physical buttons or a manufacturer's logo on the computer device 1500, the fingerprint sensor 1514 can be integrated with the physical buttons or the manufacturer's logo.

[0207] The optical sensor 1515 is used to collect the ambient light intensity. In one embodiment, the processor 1501 can control the display brightness of the display screen 1505 according to the ambient light intensity collected by the optical sensor 1515. Specifically, when the ambient light intensity is high, the display brightness of the display screen 1505 is increased; when the ambient light intensity is low, the display brightness of the display screen 1505 is decreased. In another embodiment, the processor 1501 can also dynamically adjust the shooting parameters of the camera module 1506 according to the ambient light intensity collected by the optical sensor 1515.

[0208] The proximity sensor 1516, also known as the distance sensor, is usually set on the front panel of the computer device 1500. The proximity sensor 1516 is used to collect the distance between the user and the front of the computer device 1500. In one embodiment, when the proximity sensor 1516 detects that the distance between the user and the front of the computer device 1500 is gradually decreasing, the processor 1501 controls the display screen 1505 to switch from the lit state to the off state; when the proximity sensor 1516 detects that the distance between the user and the front of the computer device 1500 is gradually increasing, the processor 1501 controls the display screen 1505 to switch from the off state to the lit state.

[0209] Those skilled in the art can understand that Figure 15 the structure shown in does not constitute a limitation on the computer device 1500, and it may include more or fewer components than shown in the figure, or combine some components, or adopt different component arrangements.

[0210] The embodiment of the present application also provides a computer-readable storage medium, in which at least one program code is stored, and the at least one program code is loaded and executed by a processor to implement the operations performed in the method for obtaining seismic exploration information in the embodiment of the present application.

[0211] The embodiments of the present application also provide a computer program product or a computer program. The computer program product or the computer program includes computer program code, and the computer program code is 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, so that the computer device performs the operations executed by the above-mentioned method for obtaining seismic exploration information.

[0212] In some embodiments, the computer program involved in the embodiments of the present application can be deployed to be executed on one computer device, or on multiple computer devices located at one location. Alternatively, it can be executed on multiple computer devices distributed at multiple locations and interconnected through a communication network. The multiple computer devices distributed at multiple locations and interconnected through a communication network can form a blockchain system.

[0213] The above 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 shall be included in the protection scope of the present application.

Claims

1. A method for obtaining seismic exploration information, characterized in that, The method includes: Obtaining multiple seismic data collected by multiple geophones in a target work area, where the seismic data is collected by the geophones within a preset acquisition 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, and the first horizontal component, the second horizontal component, and the first vertical component are perpendicular to each other pairwise; For each geophone, based on the first seismic data and the second seismic data, determining a complex parameter of particle vibration, a linear parameter of particle vibration, and a cross-correlation parameter at a target acquisition moment, where the complex parameter of particle vibration is used to reflect the complexity of particle vibration, the linear parameter of particle vibration is used to reflect the linearity of particle vibration, the cross-correlation parameter is used to reflect the cross-correlation relationship between the first seismic data and the second seismic data at the target acquisition moment, and the target acquisition moment is any acquisition moment within the preset acquisition time range; Based on the complex parameter of particle vibration, the linear parameter of particle vibration, and the cross-correlation parameter at the target acquisition moment, determining a polarization parameter at the target acquisition moment and a polarization angle at the target acquisition moment; Based on the first seismic data, the second seismic data, the polarization parameter at the target acquisition moment, and the polarization angle, determining fourth seismic data, where the fourth seismic data is the seismic data corresponding to the polarization angle at the target acquisition moment; Based on the third seismic data corresponding to each geophone and the fourth seismic data at each acquisition moment, obtaining seismic exploration information.

2. The method according to claim 1, wherein The determination process of determining the cross-correlation parameter at the target acquisition moment based on the first seismic data and the second seismic data includes: Based on the first seismic data and the second seismic data, determining the frequency distribution range of the seismic wave excited by the shot point; Based on the frequency distribution range, determining a time window parameter, where the time window parameter is used to reflect the time window width corresponding to determining the cross-correlation parameter at the target acquisition moment; Based on the time window parameter and the target acquisition moment, determining an integration time range; Based on the integration time range, the first seismic data, and the second seismic data, determining the cross-correlation parameter at the target acquisition moment.

3. The method according to claim 2, characterized in that The determining the cross-correlation parameter at the target acquisition moment 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 within the integration time range to obtain a first vibration vector, where the first vibration vector includes an amplitude and a vibration direction; Determining the absolute value of the ratio of the first vibration vector to the third amplitude to obtain the cross-correlation parameter at the target acquisition moment.

4. The method according to claim 1, wherein Determining the polarization parameter and the polarization angle at the target acquisition time based on the complex parameters of the particle vibration, the linear parameters of the particle vibration, and the cross-correlation parameter at the target acquisition time includes: Based on the complex parameters of the particle vibration, the linear parameters of the particle vibration, and the cross-correlation parameter at the target acquisition time, determining the polarization parameter at the target acquisition time, where the polarization parameter is used to reflect the linear degree of the particle vibration; When the polarization parameter is not less than the preset polarization parameter, determining the polarization angle at the target acquisition time based on the first seismic data and the second seismic data; When the polarization parameter is less than the preset polarization parameter, determining the polarization angle at the target acquisition time based on the first seismic data, the second seismic data, and the integration time range, where the integration time range is obtained when determining the cross-correlation parameter.

5. The method according to claim 4, wherein The determining the polarization angle at 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; Determining the arctangent value of the first ratio to obtain the polarization angle at the target acquisition time.

6. The method according to claim 4, wherein The determining the polarization angle at 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 within the integration time range to obtain a fourth amplitude; For each preset polarization angle within the preset polarization angle range, determining the product of the first seismic data at the target integration time and the cosine value of the preset polarization angle to obtain a fifth seismic data, where the target integration time is any integration time within the integration time range; Determining the product of the second seismic data at the target integration time and the sine value of the preset polarization angle to obtain a sixth seismic data; Determining the sum value of the fifth seismic data and the sixth seismic data within the integration time range to obtain a fifth amplitude; Determining the difference between the fourth amplitude and the fifth amplitude to obtain a sixth amplitude corresponding to the preset polarization angle; Taking the preset polarization angle corresponding to the minimum sixth amplitude as the polarization angle at the target acquisition time.

7. The method according to claim 4, wherein The determining the polarization parameter at the target acquisition time based on the complex parameters of the particle vibration, the linear parameters of the particle vibration, and the cross-correlation parameter at the target acquisition time includes: Determining the difference between the cross-correlation parameter at the target acquisition time and the linear parameters of the particle vibration to obtain a first difference; Determining the product of the complex parameters of the particle vibration and the first difference to obtain a first product; Determining the exponential value with the natural constant as the base and the negative of the first product as the exponent; Determining the reciprocal after adding 1 to the exponential value and taking the reciprocal as the polarization parameter.

8. The method according to claim 1, wherein Determining the third seismic data corresponding to the polarization angle at the target acquisition time based on the first seismic data, the second seismic data, the polarization parameter, and the polarization angle at the target acquisition time includes: Determining the product of the first seismic data at the target acquisition time and the cosine value of the polarization angle at the target acquisition time to obtain the seventh seismic data; Determining the product of the second seismic data at the target acquisition time and the sine value of the polarization angle at the target acquisition time to obtain the eighth seismic data; 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 the third seismic data corresponding to the polarization angle at the target acquisition time.

9. An apparatus for acquiring seismic exploration information, characterized in that The apparatus includes: A first acquisition module, configured to acquire a plurality of seismic data collected by a plurality of geophones in a target work area, where the seismic data is collected by the geophones within a preset acquisition time range, and the seismic data includes the first seismic data of the first horizontal component, the second seismic data of the second horizontal component, and the third seismic data of the first vertical component, and the first horizontal component, the second horizontal component, and the first vertical component are perpendicular to each other pairwise; A first determination module, configured to, for each geophone, determine a particle vibration complexity parameter, a particle vibration linearity parameter, and a cross-correlation parameter at 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 the complexity of particle vibration, the particle vibration linearity parameter is used to reflect the linearity of particle vibration, and the cross-correlation parameter is used to reflect the 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; A second determination module, configured to determine 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; A third determination module, 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 the seismic data corresponding to the polarization angle at the target acquisition time; A second acquisition module, configured to acquire seismic exploration information based on the third seismic data corresponding to each geophone and the fourth seismic data at each acquisition time.

10. A computer device, characterized in that, The computer device includes a processor and a memory, and at least one program code is stored in the memory and is loaded and executed by the processor to implement the method for acquiring seismic exploration information according to any one of claims 1 to 8.

Images