Static correction method, device and equipment for elastic wave vector seismic exploration data and medium

Abstract

The invention belongs to the technical field of seismic exploration, and particularly discloses an elastic wave vector seismic exploration data static correction method, which comprises the following steps: S1, establishing a longitudinal wave micro-logging surface velocity model; s2, establishing a three-component micro-logging shear wave surface velocity model; s3, performing low-frequency static correction on the longitudinal wave seismic data; s4, performing low-frequency static correction on the transverse wave seismic data; s5, performing low-frequency static correction on the converted wave seismic data; s6, calculating residual static correction values of shot points and detection points of the longitudinal wave field, the converted wave field and the transverse wave field; and S7, calculating total static correction values of shot points and detection points of the longitudinal wave field, the converted wave field and the transverse wave field. According to the method, the static correction problem of various wave field seismic data in elastic wave vector seismic exploration is effectively solved, and the calculation efficiency and the imaging precision are improved. The method is suitable for static correction of seismic exploration data.

Description

Technical Field

[0001] The present invention belongs to the technical field of seismic exploration, and specifically relates to a static correction method, device, equipment and medium for elastic wave vector seismic exploration data. Background Art

[0002] In recent years, as the difficulty of oil and gas exploration has gradually increased, the effective identification and efficient development of oil and gas reservoirs have faced greater difficulties. The traditional single P-wave seismic wavefield exploration method (P-wave source excitation, P-wave component reception) has faced more challenges and limitations. Therefore, the elastic wave vector seismic exploration method (single P-wave source excitation or synchronous excitation of P-wave source and S-wave source, synchronous reception of P-wave and S-wave components) has received increasing attention in the oil and gas exploration industry.

[0003] Elastic wave vector seismic exploration is also known as multi-wave multi-component seismic exploration. The important difference from conventional single P-wave exploration is that elastic wave vector seismic exploration uses multi-wave source excitation and multi-component geophones for reception, and can collect the full-wavefield vector information of underground oil and gas reservoirs, providing the possibility for accurately describing the structural morphology of underground geological bodies, finely predicting fractures, and finely depicting the spatial distribution of oil and gas reservoirs.

[0004] The technical advantages of elastic wave vector seismic exploration mainly include the following aspects: (1) For regions such as gas cloud areas, tight gas reservoirs in clastic rocks, carbonate reservoirs, and shale oil and gas reservoirs, P-wave (P-wave source excitation, P-wave component reception) data has characteristics such as low signal-to-noise ratio and weak energy, while using converted wave (P-wave source excitation, S-wave component reception) and S-wave (S-wave source excitation, S-wave component reception) information can achieve accurate imaging and reservoir characterization; (2) The elastic wave seismic exploration data has rich wavefield information, and can simultaneously obtain P-wave, converted wave, and S-wave wavefield information. In addition to accurately extracting P-wave velocity parameters, it can also more accurately obtain S-wave velocity parameters, and further extract accurate reservoir elastic parameters, such as Poisson's ratio, P-wave to S-wave velocity ratio, Young's modulus, brittleness index, etc., for lithology prediction and reservoir characteristic description; (3) The converted wave data and S-wave data are more sensitive to reservoir fractures, and can better solve the difficult problems such as quantitatively predicting the development characteristics of reservoir fractures and fluid detection using the shear wave splitting characteristics.

[0005] In elastic wave vector seismic exploration, the S-wave static correction problem has always been a bottleneck problem restricting the development of elastic wave vector seismic exploration technology. Affected by the rapid change of the S-wave surface velocity structure and the large thickness of the low-velocity layer near the surface, the difference between the S-wave velocity and the P-wave velocity near the surface is large, the S-wave static correction amount changes complexly, and the S-wave static correction amount at the same geophone position may reach 2-10 times that of the P-wave. The in-phase axis on the converted wave and S-wave imaging profiles is severely distorted, the imaging continuity is poor, and there are serious high- and low-frequency static correction problems.

[0006] Currently, the commonly used methods for solving shear-wave static correction include: P-wave coefficient method for static correction, converted-wave first-arrival picking method for static correction, common-receiver-point stacking method for shear-wave static correction, and surface-wave inversion method for shear-wave static correction. However, the above methods all have defects: the P-wave coefficient method for static correction has relatively cumbersome test steps, low operability, and poor accuracy; for the converted-wave first-arrival picking method for static correction, due to the relatively low signal-to-noise ratio of the actual converted-wave data, the converted-wave first-arrival information is often masked by other wavefield information and noise signals, and the shear-wave first-arrival information is not easy to identify and pick up; the common-receiver-point stacking method for shear-wave static correction is more applicable to multi-wave and multi-component data with simple structure and high signal-to-noise ratio in the work area, while its adaptability is poor for multi-wave and multi-component seismic data with complex structure and low signal-to-noise ratio in the work area; for the surface-wave inversion method for shear-wave static correction, the model inversion accuracy highly depends on the initial layered structure model. To avoid spatial aliasing, the original data acquisition needs to use small trace intervals and high-density reception.

[0007] Based on the above analysis, in elastic-wave vector seismic exploration, there is an urgent need for a method that can solve the shear-wave static correction problem with high efficiency and high accuracy to promote the progress and improvement of elastic-wave vector seismic exploration technology. Summary of the Invention

[0008] The object of the present invention is to provide a method, device, equipment and medium for static correction of elastic-wave vector seismic exploration data to effectively solve the static correction problem of multi-wavefield seismic data in elastic-wave vector seismic exploration and improve the calculation efficiency and imaging accuracy.

[0009] The technical method adopted by the present invention to achieve the above object is as follows:

[0010] A method for static correction of elastic-wave vector seismic exploration data includes the following steps:

[0011] S1. Conduct P-wave micro-log surface surveys in the field seismic work area, interpret the velocity of the P-wave micro-log data, and establish a P-wave micro-log surface velocity model;

[0012] S2. Conduct three-component micro-log surface surveys in the field seismic work area, interpret the velocity of the three-component micro-log data, and establish a three-component micro-log shear-wave surface velocity model;

[0013] S3. Conduct elastic-wave vector seismic exploration in the field seismic work area, pick the first arrivals of the P-wave seismic data obtained from the elastic-wave vector seismic exploration, and use the P-wave micro-log surface velocity model established in step S1 as a constraint. Adopt tomographic inversion technology based on the P-wave first-arrival travel-time information to establish a near-surface P-wave velocity model, calculate the low-frequency static correction amounts of the shot points and receiver points in the P-wave wavefield, and perform low-frequency static correction on the P-wave seismic data through the obtained low-frequency static correction amounts of the shot points and receiver points in the P-wave wavefield;

[0014] S4. Conduct elastic wave vector seismic exploration in the field seismic work area. When the source type is the synchronous excitation of a longitudinal wave source and a transverse wave source, pick the first arrivals of the transverse wave seismic data obtained from the elastic wave vector seismic exploration. Using the three-component micro-logging transverse wave surface velocity model established in step S2 as a constraint, adopt the tomography inversion technology based on the transverse wave first arrival travel time information to establish the first near-surface transverse wave velocity model, and calculate the low-frequency static correction amounts of the shot points and geophone points in the transverse wave field. When the source type is the excitation of a single longitudinal wave source, adopt the acquisition method of small trace interval and high-density reception, optimize the shot-receiver offsets of the longitudinal wave seismic data obtained from the elastic wave vector seismic exploration in step S3, and obtain the Rayleigh surface wave seismic data with obvious linear characteristics. Then, using the three-component micro-logging transverse wave surface velocity model established in step S2 as a constraint, perform inversion according to the Rayleigh surface wave dispersion characteristic curve to establish the second near-surface transverse wave velocity model with the layered structure characteristics of the low-velocity zone, calculate the low-frequency static correction amounts of the shot points and geophone points in the transverse wave field, and perform low-frequency static correction on the transverse wave seismic data through the obtained low-frequency static correction amounts of the shot points and geophone points in the transverse wave field;

[0015] S5. Conduct elastic wave vector seismic exploration in the field seismic work area. Use the acquisition method of longitudinal wave source excitation and transverse wave component reception to obtain converted wave seismic data. Perform low-frequency static correction on the shot point data of the converted wave seismic data through the low-frequency static correction amount of the shot points in the longitudinal wave field obtained in step S3, and perform low-frequency static correction on the geophone point data of the converted wave seismic data through the low-frequency static correction amount of the geophone points in the transverse wave field obtained in step S4;

[0016] S6. Perform normal moveout correction on the low-frequency static corrected longitudinal wave seismic data, converted wave seismic data, and transverse wave seismic data, calculate the surface consistent residual static correction amounts of the longitudinal wave seismic data, converted wave seismic data, and transverse wave seismic data after normal moveout correction, and obtain the residual static correction amounts of the shot points and geophone points in the longitudinal wave field, converted wave field, and transverse wave field;

[0017] S7. According to the low-frequency static correction amounts of the shot points and geophone points in the longitudinal wave field obtained in step S3, the low-frequency static correction amounts of the shot points and geophone points in the transverse wave field obtained in step S4, and the residual static correction amounts of the shot points and geophone points in the longitudinal wave field, converted wave field, and transverse wave field obtained in step S6, calculate the total static correction amounts of the shot points and geophone points in the longitudinal wave field, converted wave field, and transverse wave field.

[0018] It is specified that the specific method for establishing the longitudinal wave micro-logging surface velocity model in step S1 is as follows:

[0019] S11. Pick up the first arrival travel time information of the longitudinal wave seismic wave field at each excitation point in the longitudinal wave micro-logging data;

[0020] S12. For the P-wave micro-logging data of each excitation point, based on the time-depth curve of the P-wave first arrival time varying with depth for different seismic traces, establish a discrete point model of the P-wave layer velocity at the excitation point, which is: ;

[0021] Among them, is the formation depth of the th layer of the P-wave, is the formation depth of the th layer of the P-wave, is the first arrival time of the th layer of the P-wave, is the first arrival time of the +1th layer of the P-wave, is the P-wave velocity of the th layer;

[0022] S13. According to the discrete point models of the P-wave layer velocity at all excitation points in the field seismic work area obtained in steps S11 and S12, then perform three-dimensional spatial uniform sampling interpolation to obtain a P-wave micro-logging surface velocity model with uniform three-dimensional spatial sampling in the entire work area.

[0023] As a limitation: The specific method for establishing the three-component micro-logging shear wave surface velocity model in step S2 is as follows:

[0024] S21. Pick up the first arrival travel time information of the shear wave component seismic wave field at each excitation point in the three-component micro-logging data;

[0025] S22. For the shear wave component micro-logging data of each excitation point, based on the time-depth curve of the shear wave first arrival time varying with depth for different seismic traces, establish a discrete point model of the shear wave layer velocity at the excitation point, which is: ;

[0026] Among them, is the formation depth of the th layer of the shear wave, is the formation depth of the th layer of the shear wave, is the first arrival time of the th layer of the shear wave, is the first arrival time of the +1th layer of the shear wave, is the shear wave velocity of the th layer;

[0027] S23. According to the discrete point models of the shear wave layer velocity at all excitation points in the field seismic work area obtained in steps S21 and S22, then perform three-dimensional spatial uniform sampling interpolation to obtain a three-component micro-logging shear wave surface velocity model with uniform three-dimensional spatial sampling in the entire work area.

[0028] As a limitation: Step S3 is specifically as follows:

[0029] S31. Conduct elastic wave vector seismic exploration in the field seismic work area, and use the acquisition method of P-wave source excitation and P-wave component reception to obtain P-wave seismic data, and pick up the first arrival travel time information of the P-wave seismic wave field at each excitation point;

[0030] S32. Obtain the first arrival travel time information of the P-wave seismic wave field at all excitation points in the field seismic work area according to Step S31;

[0031] S33. With the P-wave micro-logging surface velocity model established in Step S1 as a constraint, use the first arrival travel time information of the P-wave seismic wave field at all excitation points in the field seismic work area, and adopt the tomography inversion technology based on the first arrival travel time information of the P-wave to establish a near-surface P-wave velocity model. The relationship between the first arrival travel time information of the P-wave seismic wave field and the slowness of the P-wave in the formation grid is as follows: ;

[0032] Wherein, is the first arrival travel time of the th P-wave ray from the excitation point S to the receiving point R, is the distance traveled by the P-wave ray in the th layer in the work area, is the slowness of the P-wave in the th layer;

[0033] The relationship between the P-wave slowness and the P-wave velocity is:

[0034] ;

[0035] Wherein, is the P-wave velocity of the th layer;

[0036] The formula for the first arrival travel time information of the P-wave seismic wave field and the slowness of the P-wave in the formation grid can be simplified as:

[0037] ;

[0038] Wherein, is the P-wave first arrival travel time matrix, is the P-wave ray matrix, is the P-wave slowness matrix;

[0039] Adopt the tomography inversion technology based on the first arrival travel time information of the P-wave to establish a near-surface P-wave velocity model, specifically as follows:

[0040] Establish the objective function : ;

[0041] Among them, is the initial P-wave slowness matrix obtained from the P-wave micro-logging surface velocity model established in step S1, is the perturbation amount of the P-wave slowness matrix, is the actually observed P-wave first arrival travel time matrix, is the error between the P-wave ray tracing travel time and the actually observed travel time after the

[0042] The simulated annealing algorithm is used to solve the objective function. The specific calculation method is as follows: For the th iteration calculation, a perturbation amount is randomly generated to form a new state of the objective function. Calculate the change amount caused by the perturbation amount to the objective function ; If , that is, the energy decreases, then this perturbation is accepted; if , that is, the energy increases, then the probability of accepting this perturbation is: ;

[0043] Among them, is the absolute temperature during the th iteration calculation;

[0044] After the th iteration calculation, if the perturbation is accepted, then the parameter value of the objective function is modified; otherwise, the parameter value of the original objective function is still used. Iterate in turn. When iterating to the th time, if the objective function reaches the convergence condition, the loop terminates, the current optimal solution is output, and the calculation ends. At this time, the optimal solution of the P-wave slowness matrix is expressed as: ;

[0045] Among them, is the perturbation amount of the P-wave slowness matrix during the th iteration calculation;

[0046] The finally inverted near-surface P-wave velocity model is: ;

[0047] S34. According to the near-surface P-wave velocity model , calculate the low-frequency static correction amounts of the shot points and geophone points of the P-wave wavefield;

[0048] The calculation formula for the low-frequency static correction amount of the shot points of the P-wave wavefield is as follows: ;

[0049] The calculation formula for the low-frequency static correction amount of the geophone points in the longitudinal wave field is as follows: ;

[0050] Wherein, represents the low-frequency static correction amount of the shot points in the longitudinal wave field, represents the low-frequency static correction amount of the geophone points in the longitudinal wave field, represents the static correction reference surface of the shot points, represents the static correction reference surface of the geophone points, represents the surface elevation of the shot points, represents the surface elevation of the geophone points, represents the velocity of the high-velocity layer in the near-surface longitudinal wave velocity model, represents the velocity of the low-velocity layer in the near-surface longitudinal wave velocity model, represents the thickness of the low-velocity layer at the shot point position, represents the thickness of the low-velocity layer at the geophone point position;

[0051] S35. Perform low-frequency static correction on the longitudinal wave seismic data obtained in step S31 using the low-frequency static correction amounts of the shot points and geophone points in the longitudinal wave field obtained through step S34.

[0052] As a limitation: Step S4 is specifically as follows:

[0053] S41. Conduct elastic wave vector seismic exploration in the field seismic work area;

[0054] When the source type is the synchronous excitation of a longitudinal wave source and a transverse wave source, use the acquisition method of exciting with the transverse wave source and receiving the transverse wave component to obtain transverse wave seismic data, and pick up the first arrival travel time information of the transverse wave seismic wave field at each excitation point;

[0055] According to the obtained first arrival travel time information of the transverse wave seismic wave field at all excitation points in the field seismic work area;

[0056] Constrained by the three-component micro-logging transverse wave surface velocity model established in step S2, use the first arrival travel time information of the transverse wave seismic wave field at all excitation points in the field seismic work area, and adopt the tomographic inversion technology based on the first arrival travel time information of the transverse wave to establish the first near-surface transverse wave velocity model. The first arrival travel time information of the transverse wave seismic wave field and the slowness of the transverse wave in the formation grid satisfy the following relationship: ;

[0057] Wherein, is the first arrival travel time of the th ray of the transverse wave from the excitation point S to the receiving point R, is the distance traveled by the transverse wave ray in the th layer in the work area, is the slowness of the shear wave in the layer;

[0058] The relationship between the shear wave slowness and the shear wave velocity is: ;

[0059] where is the shear wave velocity of the layer;

[0060] The first arrival travel time information of the shear wave seismic wave field and the slowness formula of the shear wave in the formation grid can be simplified to:

[0061] ;

[0062] where is the shear wave first arrival travel time matrix, is the shear wave ray matrix, is the shear wave slowness matrix;

[0063] Adopt the tomographic inversion technology based on the shear wave first arrival travel time information to establish the first near-surface shear wave velocity model, specifically:

[0064] Establish the objective function : ;

[0065] where is the initial shear wave slowness matrix obtained from the three-component micro-log shear wave surface velocity model established in step S2, is the perturbation amount of the shear wave slowness matrix, is the actually observed shear wave first arrival travel time matrix, is the th error between the shear wave ray tracing travel time and the actually observed travel time after iteration;

[0066] Adopt the simulated annealing algorithm to solve the objective function, and the specific calculation method is: for the th iteration calculation, randomly generate a perturbation amount , form a new state of the objective function, and calculate the change amount of the perturbation amount causing the objective function ; if , that is, the energy decreases, then the perturbation is accepted; if , that is, the energy increases, then the probability of accepting the perturbation is: ;

[0067] where is the absolute temperature during the th iteration calculation;

[0068] After the -th iterative calculation, if the perturbation is accepted, the parameter value of the objective function is modified; otherwise, the parameter value of the original objective function is still used. Iterate in turn. When iterating to the -th time, if the objective function reaches the convergence condition, the loop terminates, the current optimal solution is output, and the calculation ends. At this time, the optimal solution of the shear wave slowness matrix The expression is: ;

[0069] where is the perturbation of the shear wave slowness matrix during the -th iterative calculation;

[0070] The first near-surface shear wave velocity model obtained by the final inversion is: ;

[0071] When the source type is a single P-wave source excitation, a acquisition method with a small trace interval and high-density reception is adopted. The P-wave seismic data obtained from elastic wave vector seismic exploration is optimized for the shot-receiver offset, and Rayleigh surface wave seismic data with obvious linear characteristics is obtained. The Rayleigh surface wave dispersion curve is inverted by the damped least squares method to establish the second near-surface shear wave velocity model. The objective function is: ;

[0072] where is the actually observed Rayleigh surface wave dispersion curve, is the three-component micro-log shear wave surface layer velocity model established in step 2, is the perturbation of the second near-surface shear wave velocity model, is the Rayleigh surface wave dispersion curve calculated by forward modeling, is the number of inversion data;

[0073] The damped least squares solution of the objective function is: ;

[0074] where is the Jacobian matrix of, is the transpose matrix of, is the damping factor, is the identity matrix;

[0075] The second near-surface shear wave velocity model obtained by the final inversion is: ;

[0076] S42. Calculate the low-frequency static corrections for the shot points and geophone points of the shear wave field based on the first near-surface shear wave velocity model or the second near-surface shear wave velocity model obtained in step S41. The calculation formula for the low-frequency static correction of the shot points in the shear wave field is as follows: ;

[0077] The calculation formula for the low-frequency static correction of the geophone points in the shear wave field is as follows: ;

[0078] Wherein, represents the low-frequency static correction of the shot points in the shear wave field, represents the low-frequency static correction of the geophone points in the shear wave field, represents the static correction reference surface of the shot points, represents the static correction reference surface of the geophone points, represents the surface elevation of the shot points, represents the surface elevation of the geophone points, represents the velocity of the high-velocity layer of the first near-surface shear wave velocity model or the second near-surface shear wave velocity model, represents the velocity of the low-velocity layer of the first near-surface shear wave velocity model or the second near-surface shear wave velocity model, represents the thickness of the low-velocity layer at the shot point position, represents the thickness of the low-velocity layer at the geophone point position;

[0079] S43. Perform low-frequency static correction on the shear wave seismic data obtained in step S41 using the low-frequency static corrections of the shot points and geophone points of the shear wave field obtained in step S42.

[0080] As a limitation: The specific calculation of the surface-consistent residual static correction of the P-wave seismic data after normal moveout correction in step S6 is as follows: The residual static correction of any seismic trace is expressed as, ;

[0081] Wherein, is the shot point number, is the geophone point number, is the common midpoint number; represents the residual static correction of the shot point number , represents the residual static correction of the geophone point number ; is the structural term, representing the residual moveout generated by the structural undulation at the position of the th common midpoint number;

[0082] The residual static correction amount for any seismic trace Perform residual moveout decomposition, establish a system of linear equations according to the criterion of minimum squared error. Given that the residual static correction amount for any seismic trace is , and the residual static correction amount calculated in each iteration is , establish the objective function , solve the objective function for the solution with the minimum squared error, ;

[0083] Obtain the following partial derivatives,

[0084] , , ; Use an iterative algorithm to solve the system of equations to obtain the residual static correction amount at any shot point position and the residual static correction amount at any geophone point position;

[0085] Use the same method as calculating the surface-consistent residual static correction amount of the P-wave seismic data after normal moveout correction to calculate the surface-consistent residual static correction amount of the converted-wave seismic data after normal moveout correction, and obtain the residual static correction amount at any shot point position and the residual static correction amount at any geophone point position;

[0086] Use the same method as calculating the surface-consistent residual static correction amount of the P-wave seismic data after normal moveout correction to calculate the surface-consistent residual static correction amount of the S-wave seismic data after normal moveout correction, and obtain the residual static correction amount at any shot point position and the residual static correction amount .

[0087] As a limitation: Calculate the total static correction amount of the shot points and geophone points in the P-wave wavefield in step S7: ; ;

[0088] wherein, is the total static correction amount of the shot points in the P-wave wavefield, is the low-frequency static correction amount of the shot points in the P-wave wavefield, is the residual static correction amount of the shot points in the P-wave wavefield;

[0089] is the total static correction amount of the geophone points in the P-wave wavefield, is the low-frequency static correction amount of the geophone points in the P-wave wavefield, is the residual static correction amount of the geophone points in the P-wave wavefield;

[0090] Calculate the total static correction amounts of the shot points and geophone points for the converted wavefield:

[0091] ; ;

[0092] Among them, is the total static correction amount of the shot points for the converted wavefield, is the low-frequency static correction amount of the shot points for the longitudinal wavefield, is the residual static correction amount of the shot points for the converted wavefield;

[0093] is the total static correction amount of the geophone points for the converted wavefield, is the low-frequency static correction amount of the geophone points for the shear wavefield, is the residual static correction amount of the geophone points for the converted wavefield;

[0094] Calculate the total static correction amounts of the shot points and geophone points for the shear wavefield:

[0095] ; ;

[0096] Among them, is the total static correction amount of the shot points for the shear wavefield, is the low-frequency static correction amount of the shot points for the shear wavefield, is the residual static correction amount of the shot points for the shear wavefield;

[0097] is the total static correction amount of the geophone points for the shear wavefield, is the low-frequency static correction amount of the geophone points for the shear wavefield, is the residual static correction amount of the geophone points for the shear wavefield.

[0098] The present invention also discloses a static correction device for elastic wave vector seismic exploration data, including:

[0099] A longitudinal wave micro-log surface velocity model construction module, which is used to conduct a longitudinal wave micro-log surface survey in a field seismic work area, perform velocity interpretation on the longitudinal wave micro-log data, and establish a longitudinal wave micro-log surface velocity model;

[0100] A three-component micro-log shear wave surface velocity model construction module, which is used to conduct a three-component micro-log surface survey in a field seismic work area, perform velocity interpretation on the three-component micro-log data, and establish a three-component micro-log shear wave surface velocity model;

[0101] The longitudinal wave seismic data low-frequency static correction module is used for carrying out elastic wave vector seismic exploration in a field seismic work area, picking the first arrivals of the longitudinal wave seismic data obtained from the elastic wave vector seismic exploration, establishing a near-surface longitudinal wave velocity model by using tomography inversion technology based on the longitudinal wave first arrival travel time information with the established longitudinal wave micro-log surface velocity model as a constraint, calculating the low-frequency static correction amounts of the shot points and geophone points in the longitudinal wave wave field, and performing low-frequency static correction on the longitudinal wave seismic data by using the obtained low-frequency static correction amounts of the shot points and geophone points in the longitudinal wave wave field;

[0102] The shear wave seismic data low-frequency static correction module is used for carrying out elastic wave vector seismic exploration in a field seismic work area. When the source type is longitudinal wave source - shear wave source synchronous excitation, it picks the first arrivals of the shear wave seismic data obtained from the elastic wave vector seismic exploration, establishes a first near-surface shear wave velocity model by using tomography inversion technology based on the shear wave first arrival travel time information with the established three-component micro-log shear wave surface velocity model as a constraint, and calculates the low-frequency static correction amounts of the shot points and geophone points in the shear wave wave field; when the source type is single longitudinal wave source excitation, it adopts a acquisition method with small trace interval and high-density reception, optimizes the shot-receiver offsets of the longitudinal wave seismic data obtained from the elastic wave vector seismic exploration to obtain Rayleigh surface wave seismic data with obvious linear characteristics, then establishes a second near-surface shear wave velocity model with the characteristics of a layered structure in the low-velocity layer by inversion according to the Rayleigh surface wave dispersion characteristic curve with the established three-component micro-log shear wave surface velocity model as a constraint, calculates the low-frequency static correction amounts of the shot points and geophone points in the shear wave wave field, and performs low-frequency static correction on the shear wave seismic data by using the obtained low-frequency static correction amounts of the shot points and geophone points in the shear wave wave field;

[0103] The converted wave seismic data low-frequency static correction module is used for carrying out elastic wave vector seismic exploration in a field seismic work area, obtaining converted wave seismic data by using the acquisition method of longitudinal wave source excitation and shear wave component reception, performing low-frequency static correction on the shot point data of the converted wave seismic data by using the low-frequency static correction amount of the shot point in the longitudinal wave wave field, and performing low-frequency static correction on the geophone point data of the converted wave seismic data by using the low-frequency static correction amount of the geophone point in the shear wave wave field;

[0104] The wave field residual static correction amount calculation module is used for performing normal moveout correction on the longitudinally corrected wave seismic data, converted wave seismic data and shear wave seismic data after low-frequency static correction, calculating the surface consistent residual static correction amounts of the longitudinally corrected wave seismic data, converted wave seismic data and shear wave seismic data after normal moveout correction, and obtaining the residual static correction amounts of the shot points and geophone points in the longitudinal wave wave field, converted wave wave field and shear wave wave field;

[0105] The total static correction amount calculation module for wave fields is used to calculate the total static correction amounts of shot points and geophone points for the longitudinal wave field, converted wave field, and transverse wave field according to the low-frequency static correction amounts of shot points and geophone points in the longitudinal wave field and the transverse wave field, as well as the residual static correction amounts of shot points and geophone points in the longitudinal wave field, converted wave field, and transverse wave field.

[0106] The present invention also discloses an electronic device, including a memory, a processor, and a computer program stored on the memory and capable of running on the processor. When the processor executes the computer program, the above method is implemented.

[0107] The present invention also discloses a computer-readable storage medium. A computer program is stored on the computer-readable storage medium. When the computer program is executed by a processor, the above method is implemented.

[0108] Due to the adoption of the above solution, compared with the prior art, the beneficial effects achieved by the present invention are:

[0109] (1) For a static correction method for elastic wave vector seismic exploration data provided by the present invention, the first arrivals are picked using the shear wave seismic data obtained by elastic wave vector seismic exploration acquisition, and constrained by the three-component micro-log shear wave surface velocity model. Then, a near-surface shear wave inversion model is established using the tomography inversion technology based on the shear wave first arrival travel time information. Compared with traditional methods and technologies, the initial shear wave velocity model established by the present invention makes full use of the three-component micro-log shear wave velocity information, and the initial shear wave velocity model is more accurate; the shear wave velocity model obtained by inversion of the present invention directly comes from real shear wave information, and the inverted shear wave velocity model is more precise, increasing the reliability of the shear wave velocity model inversion result;

[0110] (2) The present invention can simultaneously obtain high-precision near-surface longitudinal wave and shear wave velocity models, and synchronously calculate the static correction amounts of shot points and geophone points for the longitudinal wave, converted wave, and shear wave fields. Compared with traditional methods and technologies, the implementation method of the present invention is more convenient, the calculation efficiency is faster, the method accuracy is higher, and it can effectively solve the static correction problem of seismic data of multiple wave fields in elastic wave vector seismic exploration, improving the popularization and application value of the present invention;

[0111] (3) For a static correction method for elastic wave vector seismic exploration data provided by the present invention, the simulated annealing method based on the non-linear global optimization algorithm is used to solve the tomography inversion problem of the longitudinal wave and shear wave first arrival travel time information. It can not only search in the direction of decreasing the objective function, but also search in the direction of increasing the objective function, will not be trapped in local extrema, can climb out of local extrema, has high search efficiency, can reach the global extremum, and can invert and obtain more accurate near-surface longitudinal wave and shear wave velocity models;

[0112] (4) The present invention also provides corresponding implementation devices, electronic devices, and readable storage media, further making the method more practical, and the devices, electronic devices, and readable storage media have corresponding advantages.

[0113] The present invention is applicable to the static correction of seismic exploration data. BRIEF DESCRIPTION OF THE DRAWINGS

[0114] The present invention will be further described in detail below with reference to the drawings and specific embodiments.

[0115] Figure 1 For the first arrival picking of P-wave seismic shot gather data in Embodiment 1 of the present invention;

[0116] Figure 2 For the near-surface P-wave velocity model established in Embodiment 1 of the present invention;

[0117] Figure 3 For the first arrival picking of S-wave seismic shot gather data in Embodiment 1 of the present invention;

[0118] Figure 4 For the first near-surface S-wave velocity model established in Embodiment 1 of the present invention;

[0119] Figure 5 For the total static correction amount of the geophone points in the P-wave wavefield in Embodiment 1 of the present invention;

[0120] Figure 6 For the total static correction amount of the shot points in the P-wave wavefield in Embodiment 1 of the present invention;

[0121] Figure 7 For the total static correction amount of the geophone points in the S-wave wavefield in Embodiment 1 of the present invention;

[0122] Figure 8 For the total static correction amount of the shot points in the S-wave wavefield in Embodiment 1 of the present invention;

[0123] Figure 9 For the structural block diagram of a static correction device for elastic wave vector seismic exploration data in Embodiment 2 of the present invention;

[0124] Figure 10 For the structural block diagram of an electronic device in Embodiment 2 of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0125] The present invention will be further described below with reference to the embodiments. However, those skilled in the art should understand that the present invention is not limited to the following embodiments, and any improvements and equivalent changes made on the basis of the specific embodiments of the present invention are within the scope of the claims of the present invention.

[0126] Embodiment 1 A Static Correction Method for Elastic Wave Vector Seismic Exploration Data

[0127] An elastic wave vector seismic exploration data static correction method, comprising the following steps:

[0128] S1. Conduct a longitudinal wave micro-logging surface survey in the field seismic work area, perform velocity interpretation on the longitudinal wave micro-logging data, and establish a longitudinal wave micro-logging surface velocity model; specifically including the following steps:

[0129] S11. Pick up the first arrival travel time information of the longitudinal wave seismic wave field at each excitation point in the longitudinal wave micro-logging data;

[0130] S12. For the longitudinal wave micro-logging data at each excitation point, according to the time-depth curve of the first arrival time of the longitudinal wave varying with depth for different seismic channels, establish a discrete point model of the longitudinal wave layer velocity at the excitation point, which is: ;

[0131] Wherein, is the formation depth of the th layer of the longitudinal wave, is the formation depth of the th layer of the longitudinal wave, is the first arrival time of the th layer of the longitudinal wave, is the first arrival time of the +1th layer of the longitudinal wave, is the longitudinal wave velocity of the th layer;

[0132] S13. According to the discrete point models of the longitudinal wave layer velocity at all excitation points in the field seismic work area obtained in steps S11 and S12, then perform three-dimensional spatial uniform sampling interpolation to obtain a longitudinally wave micro-logging surface velocity model with uniform three-dimensional spatial sampling in the entire work area.

[0133] S2. Conduct a three-component micro-logging surface survey in the field seismic work area, perform velocity interpretation on the three-component micro-logging data, and establish a three-component micro-logging shear wave surface velocity model; specifically including the following steps:

[0134] S21. Pick up the first arrival travel time information of the shear wave component seismic wave field at each excitation point in the three-component micro-logging data;

[0135] S22. For the shear wave component micro-logging data at each excitation point, according to the time-depth curve of the first arrival time of the shear wave varying with depth for different seismic channels, establish a discrete point model of the shear wave layer velocity at the excitation point, which is: ;

[0136] Wherein, is the formation depth of the th layer of the shear wave, is the formation depth of the th layer of the shear wave, is the first arrival time of the n-th layer of the shear wave, is the first arrival time of the (n + 1)-th layer of the shear wave, is the shear wave velocity of the n-th layer;

[0137] S23. Obtain the discrete point model of the shear wave layer velocity at all shot points in the field seismic work area according to steps S21 and S22, and then perform three-dimensional spatial uniform sampling interpolation to obtain a three-component micro-log shear wave surface velocity model with uniform three-dimensional spatial sampling in the entire work area.

[0138] S3. Conduct elastic wave vector seismic exploration in the field seismic work area, pick up the first arrivals of the longitudinal wave seismic data obtained from the elastic wave vector seismic exploration, and use the longitudinal wave micro-log surface velocity model established in step S1 as a constraint. Adopt the tomography inversion technology based on the first arrival travel time information of the longitudinal wave to establish a near-surface longitudinal wave velocity model, calculate the low-frequency static correction amounts of the shot points and geophone points in the longitudinal wave field, and perform low-frequency static correction on the longitudinal wave seismic data through the obtained low-frequency static correction amounts of the shot points and geophone points in the longitudinal wave field; specifically, it includes the following steps:

[0139] S31. Conduct elastic wave vector seismic exploration in the field seismic work area, and obtain longitudinal wave seismic data by using the acquisition method of exciting with a longitudinal wave source and receiving longitudinal wave components. Pick up the first arrival travel time information of the longitudinal wave seismic wave field at each shot point, as Figure 1 shown in the first arrival picking of the longitudinal wave seismic shot gather data;

[0140] S32. Obtain the first arrival travel time information of the longitudinal wave seismic wave field at all shot points in the field seismic work area according to step S31;

[0141] S33. Use the longitudinal wave micro-log surface velocity model established in step S1 as a constraint, and utilize the first arrival travel time information of the longitudinal wave seismic wave field at all shot points in the field seismic work area. Adopt the tomography inversion technology based on the first arrival travel time information of the longitudinal wave to establish a near-surface longitudinal wave velocity model. The near-surface longitudinal wave velocity model is as Figure 2 shown. The first arrival travel time information of the longitudinal wave seismic wave field and the slowness of the longitudinal wave in the formation grid satisfy the following relationship: ;

[0142] where, is the first arrival travel time of the i-th ray of the longitudinal wave from the shot point S to the receiving point R, is the distance traveled by the longitudinal wave ray in the n-th layer in the work area, is the slowness of the longitudinal wave in the n-th layer;

[0143] The relationship between the longitudinal wave slowness and the longitudinal wave velocity is as follows: ;

[0144] where, is the longitudinal wave velocity of the th layer;

[0145] The first arrival travel time information of the longitudinal wave seismic wave field and the slowness formula of the longitudinal wave in the formation grid can be simplified as: ;

[0146] where, is the first arrival travel time matrix of the longitudinal wave, is the longitudinal wave ray matrix, is the longitudinal wave slowness matrix;

[0147] Using the tomography inversion technology based on the first arrival travel time information of the longitudinal wave, a near-surface longitudinal wave velocity model is established, specifically:

[0148] Establish an objective function : ;

[0149] where, is the initial longitudinal wave slowness matrix obtained from the longitudinal wave micro-log surface velocity model established in step S1, is the perturbation amount of the longitudinal wave slowness matrix, is the actually observed first arrival travel time matrix of the longitudinal wave, is the th error between the travel time of the longitudinal wave ray tracing and the actually observed travel time after the

[0150] Using the simulated annealing algorithm to solve the objective function, the specific calculation method is: for the th iterative calculation, randomly generate a perturbation amount , form a new state of the objective function, and calculate the change amount of the perturbation amount causing the objective function ; if , that is, the energy decreases, then the perturbation is accepted; if , that is, the energy increases, then the probability of accepting the perturbation is: ;

[0151] where, is the absolute temperature during the th iterative calculation;

[0152] The After the -th iteration calculation, if the perturbation is accepted, the parameter value of the objective function is modified; otherwise, the parameter value of the original objective function is still used, and the iteration loop is carried out in turn. When the -th iteration is reached, if the objective function reaches the convergence condition, the loop terminates, the current optimal solution is output, the calculation ends, and the optimal solution of the P-wave slowness matrix is:

[0153] Among them, is the perturbation of the P-wave slowness matrix during the -th iteration calculation;

[0154] The near-surface P-wave velocity model finally obtained by inversion is: ;

[0155] S34. According to the near-surface P-wave velocity model , calculate the low-frequency static correction amounts of the shot points and geophone points of the P-wave wavefield;

[0156] The calculation formula for the low-frequency static correction amount of the shot points of the P-wave wavefield is as follows: ;

[0157] The calculation formula for the low-frequency static correction amount of the geophone points of the P-wave wavefield is as follows: ;

[0158] Among them, represents the low-frequency static correction amount of the shot points of the P-wave wavefield, represents the low-frequency static correction amount of the geophone points of the P-wave wavefield, represents the static correction datum plane of the shot points, represents the static correction datum plane of the geophone points, represents the surface elevation of the shot points, represents the surface elevation of the geophone points, represents the velocity of the high-velocity layer of the near-surface P-wave velocity model, represents the velocity of the low-velocity layer of the near-surface P-wave velocity model, represents the thickness of the low-velocity layer at the shot point position, represents the thickness of the low-velocity layer at the geophone point position;

[0159] S35. Perform low-frequency static correction on the P-wave seismic data obtained in step S31 using the low-frequency static correction amounts of the shot points and geophone points of the P-wave wavefield obtained in step S34.

[0160] S4. Conduct elastic wave vector seismic exploration in a field seismic work area. When the source type is the synchronous excitation of a longitudinal wave source and a transverse wave source, perform first arrival picking on the transverse wave seismic data obtained from the elastic wave vector seismic exploration. Using the three-component micro-logging transverse wave surface velocity model established in step S2 as a constraint, adopt tomographic inversion technology based on the first arrival travel time information of the transverse wave to establish the first near-surface transverse wave velocity model, and calculate the low-frequency static correction amounts of the shot points and geophone points in the transverse wave field. When the source type is the excitation of a single longitudinal wave source, adopt a acquisition method with a small trace interval and high-density reception, optimize the shot-receiver offsets for the longitudinal wave seismic data obtained from the elastic wave vector seismic exploration in step S31, and the optimized shot-receiver offsets are the longitudinal wave seismic data within the range of 50 m to 2000 m, to obtain Rayleigh surface wave seismic data with obvious linear characteristics. Then, using the three-component micro-logging transverse wave surface velocity model established in step S2 as a constraint, perform inversion according to the Rayleigh surface wave dispersion characteristic curve to establish the second near-surface transverse wave velocity model with the layered structure characteristics of a low-velocity layer, calculate the low-frequency static correction amounts of the shot points and geophone points in the transverse wave field, and perform low-frequency static correction on the transverse wave seismic data through the obtained low-frequency static correction amounts of the shot points and geophone points in the transverse wave field. Specifically, it includes the following steps:

[0161] S41. Conduct elastic wave vector seismic exploration in a field seismic work area;

[0162] When the source type is the synchronous excitation of a longitudinal wave source and a transverse wave source, use the acquisition method of transverse wave source excitation and transverse wave component reception to obtain transverse wave seismic data, and pick up the first arrival travel time information of the transverse wave seismic wave field at each excitation point, as Figure 3 shown in the first arrival picking of the transverse wave seismic shot gather data;

[0163] According to the obtained first arrival travel time information of the transverse wave seismic wave field at all excitation points in the field seismic work area;

[0164] Using the three-component micro-logging transverse wave surface velocity model established in step S2 as a constraint, and using the first arrival travel time information of the transverse wave seismic wave field at all excitation points in the field seismic work area, adopt tomographic inversion technology based on the first arrival travel time information of the transverse wave to establish the first near-surface transverse wave velocity model. The first near-surface transverse wave velocity model is as Figure 4 shown. The first arrival travel time information of the transverse wave seismic wave field and the slowness of the transverse wave in the formation grid satisfy the following relationship: ;

[0165] where, is the first arrival travel time of the th transverse wave ray from the excitation point S to the receiving point R, is the distance traveled by this transverse wave ray in the th layer in the work area, is the slowness of the transverse wave in the th layer;

[0166] The relationship between shear wave slowness and shear wave velocity is as follows: ;

[0167] wherein, is the shear wave velocity of the th layer;

[0168] The first arrival travel time information of the shear wave seismic wave field and the slowness formula of the shear wave in the formation grid can be simplified as: ;

[0169] wherein, is the shear wave first arrival travel time matrix, is the shear wave ray matrix, is the shear wave slowness matrix;

[0170] Adopt the tomography inversion technology based on the shear wave first arrival travel time information to establish the first near-surface shear wave velocity model, specifically:

[0171] Establish the objective function : ;

[0172] wherein, is the initial shear wave slowness matrix obtained from the three-component micro-log shear wave surface velocity model established in step S2, is the perturbation amount of the shear wave slowness matrix, is the actually observed shear wave first arrival travel time matrix, is the th error between the shear wave ray tracing travel time and the actually observed travel time after the

[0173] Adopt the simulated annealing algorithm to solve the objective function. The specific calculation method is: for the th iteration calculation, randomly generate a perturbation amount , form a new state of the objective function, and calculate the change amount of the perturbation amount causing the objective function ; if , that is, the energy decreases, then the perturbation is accepted; if , that is, the energy increases, then the probability of accepting this perturbation is: ;

[0174] wherein, is the absolute temperature during the th iteration calculation;

[0175] The After the -th iteration calculation, if the perturbation is accepted, the parameter value of the objective function is modified; otherwise, the parameter value of the original objective function is still used, and the iteration loop is carried out in turn. When the iteration reaches the -th time, if the objective function reaches the convergence condition, the loop terminates, the current optimal solution is output, the calculation ends, and the optimal solution of the shear wave slowness matrix is:

[0176] where is the perturbation of the shear wave slowness matrix during the -th iteration calculation;

[0177] The finally inverted first near-surface shear wave velocity model is: ;

[0178] When the source type is single P-wave source excitation, the acquisition method of small trace interval and high density reception is adopted. The shot-receiver offset of the P-wave seismic data obtained by elastic wave vector seismic exploration is optimized to obtain Rayleigh surface wave seismic data with obvious linear characteristics. The Rayleigh surface wave dispersion curve is inverted by the damped least squares method to establish the second near-surface shear wave velocity model, and the objective function is: ;

[0179] where is the actually observed Rayleigh surface wave dispersion curve, is the three-component micro-log shear wave surface layer velocity model established in step 2, is the perturbation of the second near-surface shear wave velocity model, is the Rayleigh surface wave dispersion curve calculated by forward modeling, is the number of inversion data;

[0180] The damped least squares solution of the objective function is: ;

[0181] where is the Jacobian matrix of, is the transpose matrix of, is the damping factor, is the identity matrix;

[0182] The finally inverted second near-surface shear wave velocity model is: ;

[0183] S42. Calculate the low-frequency static corrections for the shot points and geophone points of the shear wave field based on the first near-surface shear wave velocity model or the second near-surface shear wave velocity model obtained in step S41. The formula for calculating the low-frequency static correction for the shot points of the shear wave field is as follows: ;

[0184] The formula for calculating the low-frequency static correction for the geophone points of the shear wave field is as follows: ;

[0185] Where, represents the low-frequency static correction for the shot points of the shear wave field, represents the low-frequency static correction for the geophone points of the shear wave field, represents the static correction reference surface for the shot points, represents the static correction reference surface for the geophone points, represents the surface elevation of the shot point, represents the surface elevation of the geophone point, represents the velocity of the high-velocity layer of the first near-surface shear wave velocity model or the second near-surface shear wave velocity model, represents the velocity of the low-velocity layer of the first near-surface shear wave velocity model or the second near-surface shear wave velocity model, represents the thickness of the low-velocity layer at the shot point location, represents the thickness of the low-velocity layer at the geophone point location;

[0186] S43. Perform low-frequency static correction on the shear wave seismic data obtained in step S41 using the low-frequency static corrections for the shot points and geophone points of the shear wave field obtained in step S42.

[0187] S5. Conduct elastic wave vector seismic exploration in the field seismic work area, obtain converted wave seismic data using the acquisition method of exciting with a P-wave source and receiving S-wave components, perform low-frequency static correction on the shot point data of the converted wave seismic data using the low-frequency static correction for the shot points of the P-wave field obtained in step S34, and perform low-frequency static correction on the geophone point data of the converted wave seismic data using the low-frequency static correction for the geophone points of the S-wave field obtained in step S42;

[0188] S6. Perform normal moveout correction on the low-frequency static corrected P-wave seismic data, converted wave seismic data, and S-wave seismic data, calculate the surface consistent residual static corrections for the P-wave seismic data, converted wave seismic data, and S-wave seismic data after normal moveout correction, and obtain the residual static corrections for the shot points and geophone points of the P-wave field, converted wave field, and S-wave field; The specific calculation of the surface consistent residual static correction for the P-wave seismic data after normal moveout correction is: The residual static correction of any seismic trace is expressed as, ;

[0189] Among them, is the shot point number, is the geophone point number, is the common midpoint number; represents the residual static correction of the shot point number , represents the residual static correction of the geophone point number ; is the structural term, representing the residual moveout generated by the structural undulation at the position of the th common midpoint number;

[0190] Perform residual moveout decomposition on the residual static correction of any seismic trace, establish a system of linear equations according to the criterion of the minimum sum of squared errors. Given that the residual static correction of any seismic trace is , the residual static correction calculated in each iteration is , establish the objective function , solve the minimum sum of squared errors solution of the objective function , ; ;

[0191] Obtain the following partial derivatives, , , ;

[0192] Use an iterative algorithm to solve the system of equations to obtain the residual static correction at any shot point position and the residual static correction at any geophone point position;

[0193] Use the same method as calculating the surface consistent residual static correction of P-wave seismic data after normal moveout correction to calculate the surface consistent residual static correction of converted-wave seismic data after normal moveout correction, and obtain the residual static correction at any shot point position and the residual static correction at any geophone point position;

[0194] Use the same method as calculating the surface consistent residual static correction of P-wave seismic data after normal moveout correction to calculate the surface consistent residual static correction of S-wave seismic data after normal moveout correction, and obtain the residual static correction at any shot point position and the residual static correction .

[0195] S7. Calculate the total static corrections for the shot points and geophone points of the P-wave field, S-wave field, and converted wave field based on the low-frequency static corrections of the shot points and geophone points of the P-wave field obtained in step S3, the low-frequency static corrections of the shot points and geophone points of the S-wave field obtained in step S4, and the residual static corrections of the shot points and geophone points of the P-wave field, converted wave field, and S-wave field obtained in step S6; the total static correction of the geophone points of the P-wave field is as shown in Figure 5 , the total static correction of the shot points of the P-wave field is as shown in Figure 6 , the total static correction of the geophone points of the S-wave field is as shown in Figure 7 , the total static correction of the shot points of the S-wave field is as shown in Figure 8 ; the total static correction of the geophone points of the converted wave field is the same as that of the geophone points of the S-wave field, and the total static correction of the shot points of the converted wave field is the same as that of the shot points of the P-wave field;

[0196] Calculate the total static corrections for the shot points and geophone points of the P-wave field: ; ;

[0197] where is the total static correction of the shot points of the P-wave field, is the low-frequency static correction of the shot points of the P-wave field, is the residual static correction of the shot points of the P-wave field;

[0198] is the total static correction of the geophone points of the P-wave field, is the low-frequency static correction of the geophone points of the P-wave field, is the residual static correction of the geophone points of the P-wave field;

[0199] Calculate the total static corrections for the shot points and geophone points of the converted wave field: ; ;

[0200] where is the total static correction of the shot points of the converted wave field, is the low-frequency static correction of the shot points of the P-wave field, is the residual static correction of the shot points of the converted wave field;

[0201] is the total static correction of the geophone points of the converted wave field, is the low-frequency static correction of the geophone points of the S-wave field, is the residual static correction of the geophone points of the converted wave field;

[0202] Calculate the total static corrections for the shot points and geophone points of the S-wave field: ; ;

[0203] Among them, is the total static correction amount of the cross - wave field shot points, is the low - frequency static correction amount of the cross - wave field shot points, is the residual static correction amount of the cross - wave field shot points;

[0204] is the total static correction amount of the cross - wave field receiver points, is the low - frequency static correction amount of the cross - wave field receiver points, is the residual static correction amount of the cross - wave field receiver points.

[0205] Embodiment 2: An elastic - wave vector seismic exploration data static correction device, equipment and medium

[0206] An elastic - wave vector seismic exploration data static correction device, the structure of which is as Figure 9 shown, including:

[0207] A longitudinal - wave micro - logging surface velocity model construction module, which is used to carry out longitudinal - wave micro - logging surface surveys in the field seismic work area, interpret the velocity of the longitudinal - wave micro - logging data, and establish a longitudinal - wave micro - logging surface velocity model;

[0208] A three - component micro - logging shear - wave surface velocity model construction module, which is used to carry out three - component micro - logging surface surveys in the field seismic work area, interpret the velocity of the three - component micro - logging data, and establish a three - component micro - logging shear - wave surface velocity model;

[0209] A longitudinal - wave seismic data low - frequency static correction module, which is used to carry out elastic - wave vector seismic exploration in the field seismic work area, pick the first arrivals of the longitudinal - wave seismic data obtained from the elastic - wave vector seismic exploration, establish a near - surface longitudinal - wave velocity model by using the tomography inversion technology based on the longitudinal - wave first - arrival travel - time information with the established longitudinal - wave micro - logging surface velocity model as a constraint, calculate the low - frequency static correction amounts of the shot points and receiver points of the longitudinal - wave wave field, and perform low - frequency static correction on the longitudinal - wave seismic data through the obtained low - frequency static correction amounts of the shot points and receiver points of the longitudinal - wave wave field;

[0210] The shear wave seismic data low-frequency static correction module is used for carrying out elastic wave vector seismic exploration in a field seismic work area. When the source type is the synchronous excitation of a longitudinal wave source and a shear wave source, it picks the first arrivals of the shear wave seismic data obtained from the elastic wave vector seismic exploration, and uses the three-component micro-log shear wave surface velocity model established as a constraint. By using the tomography inversion technology based on the shear wave first arrival travel time information, it establishes the first near-surface shear wave velocity model and calculates the low-frequency static correction amounts of the shot points and geophone points in the shear wave field. When the source type is the excitation of a single longitudinal wave source, it adopts the acquisition method of small trace spacing and high-density reception, optimizes the shot-receiver offsets of the longitudinal wave seismic data obtained from the elastic wave vector seismic exploration, and obtains the Rayleigh surface wave seismic data with obvious linear characteristics. Then, using the three-component micro-log shear wave surface velocity model established as a constraint, it performs inversion according to the Rayleigh surface wave dispersion characteristic curve, establishes the second near-surface shear wave velocity model with the layered structure characteristics of the low-velocity layer, calculates the low-frequency static correction amounts of the shot points and geophone points in the shear wave field, and performs low-frequency static correction on the shear wave seismic data through the obtained low-frequency static correction amounts of the shot points and geophone points in the shear wave field;

[0211] The converted wave seismic data low-frequency static correction module is used for carrying out elastic wave vector seismic exploration in a field seismic work area. It obtains the converted wave seismic data by using the acquisition method of longitudinal wave source excitation and shear wave component reception, performs low-frequency static correction on the shot point data of the converted wave seismic data through the low-frequency static correction amount of the shot point in the longitudinal wave field, and performs low-frequency static correction on the geophone point data of the converted wave seismic data through the low-frequency static correction amount of the geophone point in the shear wave field;

[0212] The wave field residual static correction amount calculation module is used for performing normal moveout correction on the low-frequency static corrected longitudinal wave seismic data, converted wave seismic data and shear wave seismic data, calculating the surface consistent residual static correction amounts of the longitudinal wave seismic data, converted wave seismic data and shear wave seismic data after normal moveout correction, and obtaining the residual static correction amounts of the shot points and geophone points in the longitudinal wave field, converted wave field and shear wave field;

[0213] The wave field total static correction amount calculation module is used for calculating the total static correction amounts of the shot points and geophone points in the longitudinal wave field, converted wave field and shear wave field according to the low-frequency static correction amounts of the shot points and geophone points in the longitudinal wave field and the low-frequency static correction amounts of the shot points and geophone points in the shear wave field, as well as the residual static correction amounts of the shot points and geophone points in the longitudinal wave field, converted wave field and shear wave field.

[0214] This embodiment also provides an electronic device, the structure of which is as Figure 10 shown, including a memory, a processor, and a computer program stored on the memory and capable of running on the processor. When the processor executes the computer program, it implements the method of Embodiment 1.

[0215] This embodiment also provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the method of Embodiment 1 is implemented.

Claims

1. An elastic wave vector seismic exploration data static correction method, characterized in that Including the following steps: S1. Conduct longitudinal wave micro-logging surface surveys in the field seismic work area, perform velocity interpretation on the longitudinal wave micro-logging data, and establish a longitudinal wave micro-logging surface velocity model; S2. Conduct three-component micro-logging surface surveys in the field seismic work area, perform velocity interpretation on the three-component micro-logging data, and establish a three-component micro-logging shear wave surface velocity model; S3. Conduct elastic wave vector seismic exploration in the field seismic work area, pick the first arrivals of the longitudinal wave seismic data obtained from the elastic wave vector seismic exploration, and using the longitudinal wave micro-logging surface velocity model established in step S1 as a constraint, adopt tomography inversion technology based on longitudinal wave first arrival travel time information to establish a near-surface longitudinal wave velocity model, calculate the low-frequency static correction amounts for the shot points and geophone points of the longitudinal wave field, and perform low-frequency static correction on the longitudinal wave seismic data through the obtained low-frequency static correction amounts for the shot points and geophone points of the longitudinal wave field; S4. Conduct elastic wave vector seismic exploration in the field seismic work area. When the source type is the synchronous excitation of a longitudinal wave source - shear wave source, pick the first arrivals of the shear wave seismic data obtained from the elastic wave vector seismic exploration, and using the three-component micro-logging shear wave surface velocity model established in step S2 as a constraint, adopt tomography inversion technology based on shear wave first arrival travel time information to establish a first near-surface shear wave velocity model, and calculate the low-frequency static correction amounts for the shot points and geophone points of the shear wave field; when the source type is the excitation of a single longitudinal wave source, adopt a acquisition method with a small trace interval and high-density reception, optimize the shot-receiver offsets of the longitudinal wave seismic data obtained from the elastic wave vector seismic exploration in step S3 to obtain Rayleigh surface wave seismic data with obvious linear characteristics, and then using the three-component micro-logging shear wave surface velocity model established in step S2 as a constraint, perform inversion according to the Rayleigh surface wave dispersion characteristic curve to establish a second near-surface shear wave velocity model with a layered structure characteristic of a low-velocity zone, calculate the low-frequency static correction amounts for the shot points and geophone points of the shear wave field, and perform low-frequency static correction on the shear wave seismic data through the obtained low-frequency static correction amounts for the shot points and geophone points of the shear wave field; S5. Conduct elastic wave vector seismic exploration in the field seismic work area, obtain converted wave seismic data using the acquisition method of longitudinal wave source excitation and shear wave component reception, perform low-frequency static correction on the shot point data of the converted wave seismic data through the low-frequency static correction amount of the shot points of the longitudinal wave field obtained in step S3, and perform low-frequency static correction on the geophone point data of the converted wave seismic data through the low-frequency static correction amount of the geophone points of the shear wave field obtained in step S4; S6. Perform normal moveout correction on the longitudinally corrected wave seismic data, converted wave seismic data, and shear wave seismic data after low-frequency static correction, calculate the surface-consistent residual static correction amounts of the longitudinally corrected wave seismic data, converted wave seismic data, and shear wave seismic data after normal moveout correction to obtain the residual static correction amounts for the shot points and geophone points of the longitudinal wave field, converted wave field, and shear wave field; S7. Calculate the total static correction amounts of the shot points and geophone points for the P-wave field, converted-wave field, and S-wave field based on the low-frequency static correction amounts of the shot points and geophone points of the P-wave field obtained in step S3, the low-frequency static correction amounts of the shot points and geophone points of the S-wave field obtained in step S4, and the residual static correction amounts of the shot points and geophone points of the P-wave field, converted-wave field, and S-wave field obtained in step S6.

2. The static correction method for elastic wave vector seismic exploration data according to claim 1, wherein In step S1, the specific method for establishing the P-wave micro-logging surface velocity model is as follows: S11. Pick up the first-arrival travel time information of the P-wave seismic wave field at each excitation point in the P-wave micro-logging data. S12. For the P-wave micro-logging data at each excitation point, establish a discrete point model of the P-wave layer velocity at the excitation point according to the time-depth curve of the first-arrival time of the P-wave varying with depth for different seismic traces, which is: ； Among them, is the formation depth of the th layer of the P-wave, is the formation depth of the th layer of the P-wave, is the first arrival time of the th layer of the P-wave, is the first arrival time of the th + 1 layer of the P-wave, is the P-wave velocity of the th layer; S13. Based on the discrete point models of the P-wave layer velocity at all excitation points in the field seismic work area obtained in steps S11 and S12, perform three-dimensional spatial uniform sampling interpolation to obtain a P-wave micro-logging surface velocity model with uniform three-dimensional spatial sampling in the entire work area.

3. The static correction method for elastic wave vector seismic exploration data according to claim 1, wherein, In step S2, the specific method for establishing the S-wave surface velocity model of the three-component micro-logging is as follows: S21. Pick up the first-arrival travel time information of the S-wave component seismic wave field at each excitation point in the three-component micro-logging data. S22. For the S-wave component micro-logging data at each excitation point, establish a discrete point model of the S-wave layer velocity at the excitation point according to the time-depth curve of the first-arrival time of the S-wave varying with depth for different seismic traces, which is: ； Among them, is the formation depth of the th shear wave layer, is the formation depth of the th shear wave layer, is the first arrival time of the th shear wave layer, is the first arrival time of the th + 1 shear wave layer, is the shear wave velocity of the th layer; S23. Based on the discrete point models of the S-wave layer velocity at all excitation points in the field seismic work area obtained in steps S21 and S22, perform three-dimensional spatial uniform sampling interpolation to obtain a three-component micro-logging S-wave surface velocity model with uniform three-dimensional spatial sampling in the entire work area.

4. The static correction method for elastic wave vector seismic exploration data according to claim 1, characterized in that Step S3 is specifically as follows: S31. Conduct elastic wave vector seismic exploration in the field seismic work area, use the acquisition method of P-wave source excitation and P-wave component reception to obtain P-wave seismic data, and pick up the first-arrival travel time information of the P-wave seismic wave field at each excitation point. S32. Based on the first-arrival travel time information of the P-wave seismic wave field at all excitation points in the field seismic work area obtained in step S31. Using the P-wave micro-logging surface velocity model established in step S1 as a constraint, and using the first-arrival travel time information of the P-wave seismic wave field at all excitation points in the field seismic work area, adopt the tomography inversion technology based on the first-arrival travel time information of the P-wave to establish a near-surface P-wave velocity model. The relationship between the first-arrival travel time information of the P-wave seismic wave field and the slowness of the P-wave in the formation grid is as follows: ； Among them, is the first arrival travel time of the th longitudinal wave ray from the excitation point S to the receiving point R, is the distance traveled by the longitudinal wave ray in the th layer in the work area, is the slowness of the longitudinal wave in the th layer; The relationship formula between the P-wave slowness and the P-wave velocity is: ； Among them, is the longitudinal wave velocity of the The formula of the first-arrival travel time information of the P-wave seismic wave field and the slowness of the P-wave in the formation grid can be simplified as: ； Among them, is the first arrival travel time matrix of the P-wave, is the ray matrix of the P-wave, is the slowness matrix of the P-wave; Adopt the tomography inversion technology based on the first-arrival travel time information of the P-wave to establish a near-surface P-wave velocity model, specifically as follows: Establish the objective function : ； Among them, is the initial P-wave slowness matrix obtained from the P-wave micro-log surface velocity model established in step S1, is the perturbation of the P-wave slowness matrix, is the actually observed P-wave first arrival travel time matrix, is the error between the P-wave ray tracing travel time and the actually observed travel time after the -th iteration; The simulated annealing algorithm is used to solve the objective function. The specific calculation method is as follows: For the th iterative calculation, a perturbation amount is randomly generated to form a new state of the objective function, and the change amount caused by the perturbation amount to the objective function is calculated; if , that is, the energy decreases, then the perturbation is accepted; if , that is, the energy increases, then the probability of accepting the perturbation is: ； Among them, is the absolute temperature during the th iterative calculation; After the th iteration calculation, if the perturbation is accepted, the parameter value of the objective function is modified; otherwise, the parameter value of the original objective function is still used, and the iteration loop is carried out in turn. When it iterates to the th time, if the objective function reaches the convergence condition, the loop terminates, the current optimal solution is output, and the calculation ends. At this time, the optimal solution expression of the P-wave slowness matrix is: ； Among them, is the perturbation of the P-wave slowness matrix during the th iterative calculation; The near-surface P-wave velocity model obtained by the final inversion is as follows: ； S34. Calculate the low-frequency static correction amounts of the shot points and geophone points of the longitudinal wave field according to the near-surface longitudinal wave velocity model , and calculate the low-frequency static correction amounts of the shot points and geophone points of the longitudinal wave field; The calculation formula for the low-frequency static correction amount of the shot point in the P-wave field is as follows: ； The calculation formula for the low-frequency static correction amount of the geophone point in the P-wave field is as follows: ； Among them, represents the low-frequency static correction amount of the shot point for the longitudinal wave field, represents the low-frequency static correction amount of the geophone point for the longitudinal wave field, represents the static correction reference surface of the shot point, represents the static correction reference surface of the geophone point, represents the surface elevation of the shot point, represents the surface elevation of the geophone point, represents the velocity of the high-velocity layer in the near-surface longitudinal wave velocity model, represents the velocity of the low-velocity layer in the near-surface longitudinal wave velocity model, represents the thickness of the low-velocity layer at the shot point location, represents the thickness of the low-velocity layer at the geophone point location; S35. Perform low-frequency static correction on the P-wave seismic data obtained in step S31 using the low-frequency static correction amounts of the shot points and geophone points in the P-wave field obtained through step S34.

5. The static correction method for elastic wave vector seismic exploration data according to claim 1, wherein Step S4 is specifically as follows: S41. Conduct elastic wave vector seismic exploration in the field seismic work area; When the source type is the synchronous excitation of the P-wave source and the S-wave source, S-wave seismic data is obtained by using the acquisition method of S-wave source excitation and S-wave component reception, and the first arrival travel time information of the S-wave seismic wave field at each excitation point is picked up; According to the first arrival travel time information of the S-wave seismic wave field at all excitation points in the field seismic work area obtained; Constrained by the three-component micro-logging S-wave surface velocity model established in step S2, using the first arrival travel time information of the S-wave seismic wave field at all excitation points in the field seismic work area, and adopting the tomography inversion technology based on the first arrival travel time information of the S-wave, a first near-surface S-wave velocity model is established. The relationship between the first arrival travel time information of the S-wave seismic wave field and the slowness of the S-wave in the formation grid is as follows: ； Among them, is the first arrival travel time of the th shear wave ray from the source point S to the receiver point R, is the distance traveled by the shear wave ray in the th layer in the work area, is the slowness of the shear wave in the th layer; The relationship formula between the S-wave slowness and the S-wave velocity is: ； Among them, is the shear wave velocity of the The formula for the first arrival travel time information of the S-wave seismic wave field and the slowness of the S-wave in the formation grid can be simplified as: ； Among them, is the matrix of the first arrival time of shear waves, is the matrix of shear wave rays, is the matrix of shear wave slowness; Adopting the tomography inversion technology based on the first arrival travel time information of the S-wave, a first near-surface S-wave velocity model is established. Specifically: Establish the objective function : ； Among them, is the initial shear-wave slowness matrix obtained from the three-component micro-log shear-wave surface velocity model established in step S2, is the perturbation amount of the shear-wave slowness matrix, is the actually observed shear-wave first arrival travel time matrix, is the error between the shear-wave ray tracing travel time and the actually observed travel time after the The simulated annealing algorithm is used to solve the objective function. The specific calculation method is as follows: For the th iterative calculation, a perturbation is randomly generated to form a new state of the objective function. Calculate the change in the objective function caused by the perturbation ; If , that is, the energy decreases, then the perturbation is accepted; If , that is, the energy increases, then the probability of accepting the perturbation is: ； Among them, is the absolute temperature during the -th iterative calculation; After the th iteration calculation, if the perturbation is accepted, the parameter value of the objective function is modified; otherwise, the parameter value of the original objective function is still used, and the iteration loop is carried out in turn. When iterating to the th time, if the objective function reaches the convergence condition, the loop terminates, the current optimal solution is output, and the calculation ends. At this time, the optimal solution expression of the shear wave slowness matrix is: ； Among them, is the perturbation of the shear wave slowness matrix during the th iterative calculation; The first near-surface shear wave velocity model obtained by the final inversion is as follows: ； When the seismic source type is single P-wave source excitation, a acquisition method with small trace interval and high-density reception is adopted to optimize the offset of the P-wave seismic data obtained from elastic wave vector seismic exploration, and obtain Rayleigh surface wave seismic data with obvious linear characteristics. The damping least squares method is used to invert the Rayleigh surface wave dispersion curve to establish the second near-surface S-wave velocity model, and the objective function is as follows: ； Among them, is the Rayleigh surface wave dispersion curve of actual observation, is the three-component micro-log shear wave surface velocity model established in step 2, is the perturbation of the second near-surface shear wave velocity model, is the Rayleigh surface wave dispersion curve of forward calculation, is the number of inversion data; Objective function The damped least squares solution of is as follows: ； Among them, is the Jacobian matrix of, is the transpose matrix of, is the damping factor, is the identity matrix; The second near-surface shear wave velocity model obtained by the final inversion is as follows: ； S42. According to the first near-surface S-wave velocity model or the second near-surface S-wave velocity model obtained in step S41, calculate the low-frequency static correction amounts of the shot points and receiver points in the S-wave wave field; the calculation formula for the low-frequency static correction amount of the shot points in the S-wave wave field is as follows: ； The calculation formula for the low-frequency static correction amount of the receiver points in the S-wave wave field is as follows: ； Among them, represents the low-frequency static correction amount of the shot point for the shear wave field, represents the low-frequency static correction amount of the geophone point for the shear wave field, represents the static correction reference surface of the shot point, represents the static correction reference surface of the geophone point, represents the surface elevation of the shot point, represents the surface elevation of the geophone point, represents the high-velocity layer velocity of the first near-surface shear wave velocity model or the second near-surface shear wave velocity model, represents the low-velocity layer velocity of the first near-surface shear wave velocity model or the second near-surface shear wave velocity model, represents the low-velocity layer thickness at the shot point location, represents the low-velocity layer thickness at the geophone point location; S43. Perform low-frequency static correction on the S-wave seismic data obtained in step S41 by using the low-frequency static correction amounts of the shot points and receiver points in the S-wave wave field obtained in step S42.

6. The static correction method for elastic wave vector seismic exploration data according to claim 1, characterized in that, In step S6, the surface-consistent residual static correction amount of the P-wave seismic data after normal moveout correction is specifically calculated as follows: the residual static correction amount of any seismic trace is expressed as ； Among them, is the shot point number, is the geophone point number, is the common midpoint number; represents the residual static correction amount of the shot point number , represents the residual static correction amount of the geophone point number ; is the structural term, representing the residual moveout generated by the structural undulation at the position of the th common midpoint number; The residual static correction amount for any seismic trace Perform residual moveout decomposition, establish a system of linear equations according to the criterion of the minimum sum of squared errors. Given that the residual static correction amount for any seismic trace is , and the residual static correction amount calculated in each iteration is , establish the objective function , solve the objective function for the solution with the minimum sum of squared errors ； Obtain the following partial derivatives, ， ， ； An iterative algorithm is used to solve the system of equations to obtain the residual static correction amount at any shot point position and the residual static correction amount at any geophone position ; Calculate the surface-consistent residual static correction amount of the converted-wave seismic data after normal moveout correction by the method of calculating the same surface-consistent residual static correction amount of the P-wave seismic data after normal moveout correction, and obtain the residual static correction amount at any shot point position and the residual static correction amount at any geophone position ; Calculate the surface-consistent residual static correction amount of the shear-wave seismic data after normal moveout correction by using the method of calculating the same surface-consistent residual static correction amount of the longitudinal-wave seismic data after normal moveout correction, and obtain the residual static correction amount at any shot point position and the residual static correction amount at any geophone position .

7. The static correction method for elastic wave vector seismic exploration data according to claim 1, characterized in that Calculate the total static correction amounts of the shot points and receiver points in the P-wave wave field in step S7: ； ； Among them, is the total static correction amount of the vertical wavefield shot points, is the low-frequency static correction amount of the vertical wavefield shot points, is the residual static correction amount of the vertical wavefield shot points; is the total static correction amount of the vertical wavefield geophone points, is the low-frequency static correction amount of the vertical wavefield geophone points, is the residual static correction amount of the vertical wavefield geophone points; Calculate the total static correction amounts of the shot points and receiver points in the converted wave wave field: ； ； Among them, is the total static correction amount of the converted wave field shot point, is the low-frequency static correction amount of the longitudinal wave field shot point, is the residual static correction amount of the converted wave field shot point; For converting the total static correction amount of the converted wave field geophone points, For the low-frequency static correction amount of the geophone points in the shear wave field, For the residual static correction amount of the geophone points in the converted wave field; Calculate the total static correction amounts of the shot points and receiver points in the S-wave wave field: ； ； Among them, is the total static correction amount of the horizontal wavefield shot points, is the low-frequency static correction amount of the horizontal wavefield shot points, is the residual static correction amount of the horizontal wavefield shot points; is the total static correction amount of the geophone points in the shear wave field, is the low-frequency static correction amount of the geophone points in the shear wave field, is the residual static correction amount of the geophone points in the shear wave field.

8. An elastic wave vector seismic exploration data static correction device, characterized in that, Including: A P-wave micro-logging surface velocity model construction module, which is used to carry out P-wave micro-logging surface surveys in the field seismic work area, perform velocity interpretation on the P-wave micro-logging data, and establish a P-wave micro-logging surface velocity model; A three-component micro-logging S-wave surface velocity model construction module, which is used to carry out three-component micro-logging surface surveys in the field seismic work area, perform velocity interpretation on the three-component micro-logging data, and establish a three-component micro-logging S-wave surface velocity model; A P-wave seismic data low-frequency static correction module, which is used to carry out elastic wave vector seismic exploration in the field seismic work area, perform first arrival picking on the P-wave seismic data obtained from the elastic wave vector seismic exploration, establish a near-surface P-wave velocity model by using the established P-wave micro-logging surface velocity model as a constraint, adopt the tomography inversion technology based on the first arrival travel time information of the P-wave, calculate the low-frequency static correction amounts of the shot points and receiver points in the P-wave wave field, and perform low-frequency static correction on the P-wave seismic data by using the obtained low-frequency static correction amounts of the shot points and receiver points in the P-wave wave field; The shear-wave seismic data low-frequency static correction module is used for carrying out elastic wave vector seismic exploration in a field seismic work area. When the source type is the synchronous excitation of a longitudinal wave source and a shear wave source, it picks the first arrivals of the shear-wave seismic data obtained from the elastic wave vector seismic exploration, and, constrained by the three-component micro-log shear-wave surface velocity model established, uses the tomography inversion technology based on the shear-wave first arrival travel time information to establish the first near-surface shear-wave velocity model, and calculates the low-frequency static correction amounts of the shot points and geophone points of the shear-wave wave field; when the source type is the excitation of a single longitudinal wave source, it adopts the acquisition method of a small trace interval and high-density reception, optimizes the shot-receiver offsets of the longitudinal-wave seismic data obtained from the elastic wave vector seismic exploration to obtain Rayleigh surface wave seismic data with obvious linear characteristics, and then, constrained by the three-component micro-log shear-wave surface velocity model established, performs inversion according to the Rayleigh surface wave dispersion characteristic curve to establish the second near-surface shear-wave velocity model with the layered structure characteristics of a low-velocity layer, calculates the low-frequency static correction amounts of the shot points and geophone points of the shear-wave wave field, and performs low-frequency static correction on the shear-wave seismic data through the obtained low-frequency static correction amounts of the shot points and geophone points of the shear-wave wave field; The converted-wave seismic data low-frequency static correction module is used for carrying out elastic wave vector seismic exploration in a field seismic work area, obtains the converted-wave seismic data by using the acquisition method of longitudinal wave source excitation and shear wave component reception, performs low-frequency static correction on the shot point data of the converted-wave seismic data through the low-frequency static correction amount of the shot point of the longitudinal-wave wave field, and performs low-frequency static correction on the geophone point data of the converted-wave seismic data through the low-frequency static correction amount of the geophone point of the shear-wave wave field; The wave field residual static correction amount calculation module is used for performing normal moveout correction on the low-frequency static corrected longitudinal-wave seismic data, converted-wave seismic data and shear-wave seismic data, calculating the surface-consistent residual static correction amounts of the longitudinal-wave seismic data, converted-wave seismic data and shear-wave seismic data after the normal moveout correction, and obtaining the residual static correction amounts of the shot points and geophone points of the longitudinal-wave wave field, converted-wave wave field and shear-wave wave field; The wave field total static correction amount calculation module is used for calculating the total static correction amounts of the shot points and geophone points of the longitudinal-wave wave field, converted-wave wave field and shear-wave wave field according to the low-frequency static correction amounts of the shot points and geophone points of the longitudinal-wave wave field and the shear-wave wave field, and the residual static correction amounts of the shot points and geophone points of the longitudinal-wave wave field, converted-wave wave field and shear-wave wave field; 9. An electronic device, characterized in that, It includes a memory, a processor, and a computer program stored on the memory and capable of running on the processor. When the processor executes the computer program, it implements the method described in any one of claims 1-7.

10. A computer-readable storage medium, characterized in that, A computer program is stored on the computer-readable storage medium. When the computer program is executed by the processor, it implements the method described in any one of claims 1-7.

Images