阅读说明：本技术 三维初至波剩余静校正方法及装置 (Three-dimensional first-motion wave residual static correction method and device ) 是由 曹中林 周强 何光明 熊晶璇 李乐 吕文彪 于 2019-10-30 设计创作，主要内容包括：本发明公开了一种三维初至波剩余静校正方法及装置,该方法包括：从炮点道集数据中获取各检波点的与炮点对应的原有初至旅行时间；根据各检波点相对于所述炮点的位置采用曲面拟合方法对各检波点的原有初至旅行时间进行拟合,得到各检波点的拟合初至时间；计算各检波点的拟合初至时间与原有初至旅行时间的差值,并计算所有检波点的该差值的平均值,得到所述炮点的初至剩余静校正量。本发明解决了现有剩余静校正方法结果不够准确的技术问题。(The invention discloses a three-dimensional first-motion wave residual static correction method and a device, wherein the method comprises the following steps: acquiring original first-arrival travel time, corresponding to the shot point, of each wave detection point from shot point gather data; fitting the original first-arrival travel time of each demodulator probe by adopting a surface fitting method according to the position of each demodulator probe relative to the shot point to obtain the fitting first-arrival time of each demodulator probe; and calculating the difference value between the fitting first arrival time of each demodulator probe and the original first arrival travel time, and calculating the average value of the difference values of all the demodulator probes to obtain the first arrival residual static correction value of the shot point. The invention solves the technical problem that the existing residual static correction method is not accurate enough in result.)

1. A three-dimensional first-motion wave residual static correction method is characterized by comprising the following steps:

acquiring original first-arrival travel time, corresponding to the shot point, of each wave detection point from shot point gather data;

fitting the original first-arrival travel time of each demodulator probe by adopting a surface fitting method according to the position of each demodulator probe relative to the shot point to obtain the fitting first-arrival time of each demodulator probe;

and calculating the difference value between the fitting first arrival time of each demodulator probe and the original first arrival travel time, and calculating the average value of the difference values of all the demodulator probes to obtain the first arrival residual static correction value of the shot point.

2. The three-dimensional first-arrival wave residual static correction method according to claim 1, further comprising:

acquiring original first-arrival travel time of the wave detection point corresponding to each shot point from the wave detection point gather data;

fitting original first-arrival travel time of the demodulator probe corresponding to each shot point by adopting a surface fitting method according to the position of each shot point relative to the demodulator probe to obtain the fitting first-arrival time of the demodulator probe corresponding to each shot point;

and calculating the difference value between the fitting first arrival time of the demodulator probe corresponding to each shot point and the original first arrival travel time, and calculating the average value of the difference values of the demodulator probe corresponding to each shot point to obtain the first arrival residual static correction value of the demodulator probe.

3. The three-dimensional first-arrival wave residual static correction method according to claim 1, wherein the fitting of the original first-arrival travel time of each demodulator probe by using a surface fitting method according to the position of each demodulator probe relative to the shot point specifically comprises:

establishing a plane coordinate system by taking the shot as an origin to divide the shot gather data into four quadrants, wherein the X-axis direction of the plane coordinate system is consistent with the line detection direction;

and respectively fitting the original first arrival travel time of each demodulator probe in the quadrant by adopting a surface fitting method according to the coordinates of each demodulator probe in the quadrant aiming at each quadrant to obtain the fitting result of each quadrant.

4. The three-dimensional first-arrival wave residual static correction method according to claim 3, wherein the fitting of the original first-arrival travel time of each demodulator probe in the quadrant by adopting a surface fitting method according to the coordinates of each demodulator probe in the quadrant for each quadrant specifically comprises:

calculating the offset distance between each detection point and the shot point according to the coordinates of each detection point;

establishing a three-dimensional coordinate system, and determining the coordinates of each demodulator probe in the three-dimensional coordinate system according to the offset distance of each demodulator probe, the Y-axis coordinate of each demodulator probe and the original first arrival travel time of each demodulator probe;

and fitting the original first arrival travel time of each wave detection point by adopting a surface fitting method according to the coordinates of each wave detection point in the three-dimensional coordinate system.

5. The three-dimensional first-arrival wave residual static correction method according to claim 2, wherein the fitting of the original first-arrival travel time of each shot point corresponding to each shot point by adopting a surface fitting method according to the position of each shot point relative to the corresponding shot point specifically comprises:

establishing a plane coordinate system by taking the wave detection point as an origin to divide the wave detection point gather data into four quadrants, wherein the X-axis direction of the plane coordinate system is consistent with the direction of the shot line;

and respectively fitting the original first arrival travel time of the demodulator probe corresponding to each shot point by adopting a surface fitting method aiming at each quadrant according to the coordinates of each shot point in the quadrant to obtain the fitting result of each quadrant.

6. The three-dimensional first-arrival wave residual static correction method according to claim 5, wherein the fitting the original first-arrival travel time of the demodulator probe corresponding to each shot point by adopting a surface fitting method according to the coordinates of each shot point in a quadrant for each quadrant specifically comprises:

calculating the offset distance between the detection point and each shot point according to the coordinates of each shot point;

establishing a three-dimensional coordinate system, and determining the coordinates of each shot point in the three-dimensional coordinate system according to the offset distance of each shot point, the Y-axis coordinate of each shot point and the original first arrival travel time of the demodulator probe corresponding to each shot point;

and fitting the original first-arrival travel time of the wave detection point corresponding to each shot point by adopting a surface fitting method according to the coordinates of each shot point in the three-dimensional coordinate system.

7. The three-dimensional first-arrival residual static correction method according to claim 1 or 2, wherein the surface fitting method comprises: a two-dimensional scattered point BSpline curved surface fitting method based on a thin plate.

8. A three-dimensional first-motion wave residual static correction device is characterized by comprising:

the first original first-arrival travel time acquisition unit is used for acquiring original first-arrival travel time, corresponding to the shot point, of each wave detection point from shot point gather data;

the first fitting first-arrival time determining unit is used for fitting the original first-arrival travel time of each demodulator probe by adopting a surface fitting method according to the position of each demodulator probe relative to the shot point to obtain the fitting first-arrival time of each demodulator probe;

and the shot point first arrival residual static correction value calculating unit is used for calculating the difference value between the fitting first arrival time of each demodulator probe and the original first arrival travel time, and calculating the average value of the difference values of all the demodulator probes to obtain the shot point first arrival residual static correction value.

9. The three-dimensional first-arrival wave residual static correction device according to claim 8, further comprising:

a second original first-arrival travel time obtaining unit, configured to obtain original first-arrival travel time of the probe corresponding to each shot point from the probe gather data;

the second fitting first-arrival time determining unit is used for fitting the original first-arrival travel time of the wave detection point corresponding to each shot point by adopting a surface fitting method according to the position of each shot point relative to the wave detection point to obtain the fitting first-arrival time of the wave detection point corresponding to each shot point;

and the demodulator probe first arrival residual static correction value calculating unit is used for calculating the difference value between the fitting first arrival time of the demodulator probe corresponding to each shot point and the original first arrival travel time, and calculating the average value of the difference value of the demodulator probe corresponding to each shot point to obtain the first arrival residual static correction value of the demodulator probe.

10. The three-dimensional first-arrival wave residual static correction device according to claim 8, wherein the first fitted first-arrival time determination unit includes:

the first plane coordinate system establishing module is used for establishing a plane coordinate system by taking the shot as an origin to divide the shot gather data into four quadrants, wherein the X-axis direction of the plane coordinate system is consistent with the line detection direction;

and the first fitting module is used for fitting the original first arrival travel time of each demodulator probe in the quadrant by adopting a surface fitting method according to the coordinates of each demodulator probe in the quadrant aiming at each quadrant respectively to obtain the fitting result of each quadrant.

11. The three-dimensional first-arrival wave residual static correction device according to claim 10, wherein the first fitting module comprises:

the first offset distance calculation submodule is used for calculating the offset distance between each detection point and the shot point according to the coordinates of each detection point;

the first three-dimensional coordinate system establishing submodule is used for establishing a three-dimensional coordinate system and determining the coordinates of each wave detection point in the three-dimensional coordinate system according to the offset distance of each wave detection point, the Y-axis coordinate of each wave detection point and the original first arrival travel time of each wave detection point;

and the first surface fitting submodule is used for fitting the original first arrival travel time of each wave detection point by adopting a surface fitting method according to the coordinate of each wave detection point in the three-dimensional coordinate system.

12. The three-dimensional first-arrival wave residual static correction device according to claim 9, wherein the second fitted first-arrival time determination unit includes:

the second plane coordinate system establishing module is used for establishing a plane coordinate system by taking the demodulator probe as an origin to divide the demodulator probe gather data into four quadrants, wherein the X-axis direction of the plane coordinate system is consistent with the direction of the gun line;

and the second fitting module is used for fitting the original first arrival travel time of the demodulator probe corresponding to each shot point by adopting a surface fitting method according to the coordinates of each shot point in each quadrant respectively to obtain the fitting result of each quadrant.

13. The three-dimensional first-arrival wave residual static correction device according to claim 12, wherein the second fitting module comprises:

the second offset distance calculation submodule is used for calculating the offset distance between the detection point and each shot point according to the coordinates of each shot point;

the second three-dimensional coordinate system establishing submodule is used for establishing a three-dimensional coordinate system and determining the coordinates of each shot point in the three-dimensional coordinate system according to the offset distance of each shot point, the Y-axis coordinate of each shot point and the original first arrival travel time of the wave detection point corresponding to each shot point;

and the second surface fitting submodule is used for fitting the original first arrival travel time of the wave detection point corresponding to each shot point by adopting a surface fitting method according to the coordinate of each shot point in the three-dimensional coordinate system.

14. The three-dimensional first-arrival residual static correction device according to claim 8 or 9, wherein the surface fitting method comprises: a two-dimensional scattered point BSpline curved surface fitting method based on a thin plate.

15. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the steps of the method of any of claims 1 to 7 are implemented when the computer program is executed by the processor.

16. A computer-readable storage medium, in which a computer program is stored which, when being executed in a computer processor, carries out the steps of the method according to any one of claims 1 to 7.

Technical Field

The invention relates to a residual static correction method, in particular to a three-dimensional first-motion wave residual static correction method and device.

Background

Over the years, although seismic data processing workers have greatly developed the automatic residual static correction technology based on reflected waves, when the signal-to-noise ratio of the reflected waves is too low or the residual static correction value is larger than the half period of the reflected waves, the data processing cannot obtain a stable result, the generated error is large, and the requirement of complex mountain three-dimensional low signal-to-noise ratio seismic data processing cannot be met. Therefore, the prior art lacks a more accurate residual static correction method.

Disclosure of Invention

In order to solve at least one of the above technical problems, the present invention provides a three-dimensional first-motion wave residual static correction method and apparatus.

In order to achieve the above object, according to one aspect of the present invention, there is provided a three-dimensional first-arrival residual static correction method, including:

acquiring original first-arrival travel time, corresponding to the shot point, of each wave detection point from shot point gather data;

fitting the original first-arrival travel time of each demodulator probe by adopting a surface fitting method according to the position of each demodulator probe relative to the shot point to obtain the fitting first-arrival time of each demodulator probe;

and calculating the difference value between the fitting first arrival time of each demodulator probe and the original first arrival travel time, and calculating the average value of the difference values of all the demodulator probes to obtain the first arrival residual static correction value of the shot point.

Optionally, the method further includes:

acquiring original first-arrival travel time of the wave detection point corresponding to each shot point from the wave detection point gather data;

fitting original first-arrival travel time of the demodulator probe corresponding to each shot point by adopting a surface fitting method according to the position of each shot point relative to the demodulator probe to obtain the fitting first-arrival time of the demodulator probe corresponding to each shot point;

and calculating the difference value between the fitting first arrival time of the demodulator probe corresponding to each shot point and the original first arrival travel time, and calculating the average value of the difference values of the demodulator probe corresponding to each shot point to obtain the first arrival residual static correction value of the demodulator probe.

Optionally, the fitting the original first arrival travel time of each demodulator probe by using a surface fitting method according to the position of each demodulator probe relative to the shot point specifically includes:

establishing a plane coordinate system by taking the shot as an origin to divide the shot gather data into four quadrants, wherein the X-axis direction of the plane coordinate system is consistent with the line detection direction;

and respectively fitting the original first arrival travel time of each demodulator probe in the quadrant by adopting a surface fitting method according to the coordinates of each demodulator probe in the quadrant aiming at each quadrant to obtain the fitting result of each quadrant.

Optionally, the fitting, by using a surface fitting method, of the original first arrival travel time of each demodulator probe in the quadrant according to the coordinates of each demodulator probe in the quadrant for each quadrant respectively includes:

calculating the offset distance between each detection point and the shot point according to the coordinates of each detection point;

establishing a three-dimensional coordinate system, and determining the coordinates of each demodulator probe in the three-dimensional coordinate system according to the offset distance of each demodulator probe, the Y-axis coordinate of each demodulator probe and the original first arrival travel time of each demodulator probe;

and fitting the original first arrival travel time of each wave detection point by adopting a surface fitting method according to the coordinates of each wave detection point in the three-dimensional coordinate system.

Optionally, the fitting the original first arrival travel time of each shot point corresponding to each shot point by using a surface fitting method according to the position of each shot point relative to the corresponding probe point specifically includes:

establishing a plane coordinate system by taking the wave detection point as an origin to divide the wave detection point gather data into four quadrants, wherein the X-axis direction of the plane coordinate system is consistent with the direction of the shot line;

and respectively fitting the original first arrival travel time of the demodulator probe corresponding to each shot point by adopting a surface fitting method aiming at each quadrant according to the coordinates of each shot point in the quadrant to obtain the fitting result of each quadrant.

Optionally, the fitting, by using a surface fitting method, of the original first arrival travel time of the demodulator probe corresponding to each shot point according to the coordinates of each shot point in the quadrant for each quadrant specifically includes:

calculating the offset distance between the detection point and each shot point according to the coordinates of each shot point;

establishing a three-dimensional coordinate system, and determining the coordinates of each shot point in the three-dimensional coordinate system according to the offset distance of each shot point, the Y-axis coordinate of each shot point and the original first arrival travel time of the demodulator probe corresponding to each shot point;

and fitting the original first-arrival travel time of the wave detection point corresponding to each shot point by adopting a surface fitting method according to the coordinates of each shot point in the three-dimensional coordinate system.

Optionally, the surface fitting method includes: a two-dimensional scattered point BSpline curved surface fitting method based on a thin plate.

In order to achieve the above object, according to another aspect of the present invention, there is provided a three-dimensional first-arrival residual static correction apparatus, comprising:

the first original first-arrival travel time acquisition unit is used for acquiring original first-arrival travel time, corresponding to the shot point, of each wave detection point from shot point gather data;

the first fitting first-arrival time determining unit is used for fitting the original first-arrival travel time of each demodulator probe by adopting a surface fitting method according to the position of each demodulator probe relative to the shot point to obtain the fitting first-arrival time of each demodulator probe;

and the shot point first arrival residual static correction value calculating unit is used for calculating the difference value between the fitting first arrival time of each demodulator probe and the original first arrival travel time, and calculating the average value of the difference values of all the demodulator probes to obtain the shot point first arrival residual static correction value.

Optionally, the apparatus further comprises:

a second original first-arrival travel time obtaining unit, configured to obtain original first-arrival travel time of the probe corresponding to each shot point from the probe gather data;

the second fitting first-arrival time determining unit is used for fitting the original first-arrival travel time of the wave detection point corresponding to each shot point by adopting a surface fitting method according to the position of each shot point relative to the wave detection point to obtain the fitting first-arrival time of the wave detection point corresponding to each shot point;

and the demodulator probe first arrival residual static correction value calculating unit is used for calculating the difference value between the fitting first arrival time of the demodulator probe corresponding to each shot point and the original first arrival travel time, and calculating the average value of the difference value of the demodulator probe corresponding to each shot point to obtain the first arrival residual static correction value of the demodulator probe.

Optionally, the first fitting first arrival time determining unit includes:

the first plane coordinate system establishing module is used for establishing a plane coordinate system by taking the shot as an origin to divide the shot gather data into four quadrants, wherein the X-axis direction of the plane coordinate system is consistent with the line detection direction;

and the first fitting module is used for fitting the original first arrival travel time of each demodulator probe in the quadrant by adopting a surface fitting method according to the coordinates of each demodulator probe in the quadrant aiming at each quadrant respectively to obtain the fitting result of each quadrant.

Optionally, the first fitting module includes:

the first offset distance calculation submodule is used for calculating the offset distance between each detection point and the shot point according to the coordinates of each detection point;

the first three-dimensional coordinate system establishing submodule is used for establishing a three-dimensional coordinate system and determining the coordinates of each wave detection point in the three-dimensional coordinate system according to the offset distance of each wave detection point, the Y-axis coordinate of each wave detection point and the original first arrival travel time of each wave detection point;

and the first surface fitting submodule is used for fitting the original first arrival travel time of each wave detection point by adopting a surface fitting method according to the coordinate of each wave detection point in the three-dimensional coordinate system.

Optionally, the second fitted first arrival time determining unit includes:

the second plane coordinate system establishing module is used for establishing a plane coordinate system by taking the demodulator probe as an origin to divide the demodulator probe gather data into four quadrants, wherein the X-axis direction of the plane coordinate system is consistent with the direction of the gun line;

and the second fitting module is used for fitting the original first arrival travel time of the demodulator probe corresponding to each shot point by adopting a surface fitting method according to the coordinates of each shot point in each quadrant respectively to obtain the fitting result of each quadrant.

Optionally, the second fitting module includes:

the second offset distance calculation submodule is used for calculating the offset distance between the detection point and each shot point according to the coordinates of each shot point;

the second three-dimensional coordinate system establishing submodule is used for establishing a three-dimensional coordinate system and determining the coordinates of each shot point in the three-dimensional coordinate system according to the offset distance of each shot point, the Y-axis coordinate of each shot point and the original first arrival travel time of the wave detection point corresponding to each shot point;

and the second surface fitting submodule is used for fitting the original first arrival travel time of the wave detection point corresponding to each shot point by adopting a surface fitting method according to the coordinate of each shot point in the three-dimensional coordinate system.

Optionally, the surface fitting method includes: a two-dimensional scattered point BSpline curved surface fitting method based on a thin plate.

In order to achieve the above object, according to another aspect of the present invention, there is also provided a computer device, including a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the steps of the three-dimensional first-arrival residual static correction method when executing the computer program.

In order to achieve the above object, according to another aspect of the present invention, there is also provided a computer-readable storage medium storing a computer program which, when executed in a computer processor, implements the steps in the above three-dimensional first arrival residual static correction method.

The invention has the beneficial effects that: the embodiment of the invention carries out residual static correction through the first arrival wave, and has better correction effect and more accurate correction compared with the residual static correction technology based on the reflected wave in the prior art. In addition, the embodiment of the invention performs surface fitting on the original first-arrival travel time data of the three-dimensional wave detection point to obtain smooth first-arrival travel time data, and compared with the conventional curve fitting method, the method has the advantages that the calculated residual static correction value is more accurate, and the integral imaging quality of the three-dimensional low signal-to-noise ratio seismic data is improved.

Drawings

In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to these drawings without creative efforts. In the drawings:

FIG. 1 is a first flowchart of a three-dimensional first-arrival residual static correction method according to an embodiment of the present invention;

FIG. 2 is a second flowchart of the three-dimensional first-arrival residual static correction method according to the embodiment of the present invention;

FIG. 3 is a first flowchart of the present invention for fitting the original first arrival travel time of the demodulator probe;

FIG. 4 is a second flowchart of the present invention for fitting the original first arrival travel time of the demodulator probe;

FIG. 5 is a first flowchart of the present invention for fitting the wave-detection points to the original first-arrival travel time corresponding to each shot point;

FIG. 6 is a second flowchart of the present invention for fitting the wave-detection points to the original first-arrival travel time corresponding to each shot;

FIG. 7 is a first block diagram of a three-dimensional first-arrival residual static correction apparatus according to an embodiment of the present invention;

FIG. 8 is a second block diagram of the three-dimensional first-motion wave residual static correction apparatus according to the embodiment of the present invention;

FIG. 9 is a schematic diagram of establishing a planar coordinate system according to an embodiment of the present invention;

FIG. 10 is a three-dimensional representation of the original first arrival travel time of the probe prior to fitting according to an embodiment of the present invention;

FIG. 11 is a three-dimensional representation of the original first arrival travel time of the fitting demodulator probe according to the embodiment of the invention;

FIG. 12 is a schematic cross-sectional view of a 3D stack of an embodiment of the present invention with only tomographic static correction;

FIG. 13 is a schematic 3D overlay cross-section of an embodiment of the present invention undergoing tomographic and residual static correction processing in accordance with the present invention;

FIG. 14 is a schematic diagram of a computer apparatus according to an embodiment of the present invention.

Detailed Description

In order to make the technical solutions of the present invention better understood, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.

As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.

It should be noted that the terms "comprises" and "comprising," and any variations thereof, in the description and claims of the present invention and the above-described drawings, are intended to cover a non-exclusive inclusion, such that a process, method, system, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those steps or elements expressly listed, but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus.

It should be noted that the embodiments and features of the embodiments may be combined with each other without conflict. The present invention will be described in detail below with reference to the embodiments with reference to the attached drawings.

Fig. 1 is a first flowchart of a three-dimensional first-arrival wave residual static correction method according to an embodiment of the present invention, and as shown in fig. 1, the three-dimensional first-arrival wave residual static correction method according to the present embodiment includes steps S101 to S103.

And step S101, acquiring the original first arrival travel time of each demodulator probe corresponding to the shot point from shot point gather data.

In the embodiment of the invention, the seismic data can be arranged according to the shot points to form shot gather data, and then the original first-arrival travel time of all seismic channels is read on the shot gather data, so that the original first-arrival travel time of each demodulator probe corresponding to the shot points is obtained.

And S102, fitting the original first arrival travel time of each demodulator probe by adopting a surface fitting method according to the position of each demodulator probe relative to the shot point to obtain the fitting first arrival time of each demodulator probe.

In the embodiment of the present invention, the original first arrival travel time data of each probe point is three-dimensional data, and as shown in fig. 10, in the three-dimensional coordinate system, the vertical axis coordinate of each probe point is the original first arrival travel time, and the horizontal axis coordinate of each probe point is determined by the positions of the probe point and the shot point. The original first arrival travel time data of each demodulator probe can be represented by points in a three-dimensional coordinate system.

In the embodiment of the invention, the original first-motion travel time data of each demodulator probe can be fitted by adopting a surface fitting method to obtain the fitting first-motion time of each demodulator probe. In embodiments of the present invention, prior art surface fitting methods, such as cubic splines, may be used for fitting. In the preferred embodiment of the invention, a two-dimensional scattered point BSpline curved surface fitting method based on a thin plate is adopted for fitting to obtain the fitting first-arrival time of each demodulator probe.

And step S103, calculating the difference value between the fitting first arrival time of each demodulator probe and the original first arrival travel time, and calculating the average value of the difference values of all demodulator probes to obtain the first arrival residual static correction value of the shot point.

In the embodiment of the invention, the obtained first arrival residual static correction value of the shot point can be applied to correction of various seismic data and can be used for improving the overall imaging quality of the three-dimensional low signal-to-noise ratio seismic data.

From the above description, it can be seen that the residual static correction is performed by using the first-arrival wave in the embodiments of the present invention, and the correction effect is better and the correction is more accurate compared to the residual static correction technique based on the reflected wave in the prior art. In addition, the embodiment of the invention performs surface fitting on the original first-arrival travel time data of the three-dimensional wave detection point to obtain smooth first-arrival travel time data, and compared with the conventional curve fitting method, the method has the advantages that the calculated residual static correction value is more accurate, and the integral imaging quality of the three-dimensional low signal-to-noise ratio seismic data is improved.

In an embodiment of the present invention, the fitting of the original first arrival travel time of each probe point by using a surface fitting method according to the position of each probe point relative to the shot point in step S102 to obtain the fitting first arrival time of each probe point may specifically include the following steps.

Selecting data of a shot gather, establishing a plane coordinate system by taking the shot as a center, and dividing the data of the shot gather into four quadrants, as shown in fig. 9, wherein a receiving point in fig. 9 is a demodulator probe. As shown in fig. 9, the receiver points in the shot gather data are regularly arranged to form a plurality of detection lines, and in the embodiment of the present invention, the X-axis direction may be the same as the detection line direction when the plane coordinate system is established.

Fitting the original first-arrival travel time of each demodulator probe of each quadrant by adopting a sheet-based two-dimensional scattered point BSpline curved surface fitting method to obtain fitting first-arrival time, wherein the specific implementation formula is as follows:

suppose thatT_ijThe ith shot point corresponds to the jth wave detection point and the first arrival time of the wave detection point respectively,which represents the offset distance, is the distance between the two electrodes,T_ijhas the functional relation:

wherein the content of the first and second substances,

for convenient marking, after determining the number i of the shot points, the functional relation can be simplified into a standard functional equation

Then a rectangular area [ a, b ]]×[c,d]Inner scatter point x_n＝(x_n,1,x_n,2) Has a displacement d_nN1.. N, consider the problem of energy minimization of a two-dimensional sheet:

wherein:

u∈U，u＝u(x₁,x₂)，u_n＝u(x_n)＝u(x_n,1,x_n,2)

u is a function space.

According to the principle of variation, forComprises the following steps:

constructing a discrete problem by dividing a rectangular area [ a, b ] × [ c, d ] into rectangular units,

x_i＝a+(i-1)h_x,i＝1..m₁ (6)

y_j＝c+(j-1)h_y,j＝1..m₂

wherein the content of the first and second substances,representing the step sizes in the x and y directions, respectively.

By b_i(x) Denotes a position in x ═ x_iOne-dimensional BSpline at rectangular grid points (x)_i,y_j) The two-dimensional BSpline of (A) can represent B by a tensor product_k＝b_i(x)b_j(y),k＝m₁×(i-1)+j，i＝1..m₁，j＝1..m₂Make U_m＝span{B₁,B₂,B₃,…,B_mWhere m is (m)₁+3)×(m₂+3)。

Order toThen there is a discrete equation:

namely:

wherein:

B_k,n＝B_k(x_i,y_j),n＝m₁(i-1)+j，in order to ensure that the water-soluble organic acid,

A＝[∫_Ω(B_k,11B_l,11+2B_k,12B_l,12+B_k,22B_l,22)dxdy]_m×m (9)

wherein, B_k,n＝B_k(x_n)，d＝(d₁,…,d_N)^T。

The discrete equation can be expressed as:

the fitting value of the two-dimensional scattered point curved surface can be obtained by solving the discrete equation (11), namely the fitting first arrival time of the wave detection point is obtained.

In an embodiment of the present invention, the method for obtaining the first arrival residual static correction value of the shot point in step S103 may specifically be to use the shot point S_iSubtracting the original first arrival time from the fitting first arrival time corresponding to each detection point, and averaging all the receiving points to obtain a shot point S_iFirst arrival residual static correction amount ST_iThe concrete formula is as follows:

wherein M is the number of detection points, T'_ijAnd (4) corresponding fitting first arrival time of the wave detection points.

Fig. 10 is a three-dimensional schematic diagram of the original first arrival travel time of the demodulator probe before fitting according to the embodiment of the present invention, fig. 11 is a three-dimensional schematic diagram of the fitted first arrival travel time data of the demodulator probe after fitting according to the embodiment of the present invention, and as can be seen from fig. 10 and fig. 11, by the above-mentioned surface fitting method according to the embodiment of the present invention, the original first arrival travel time data of the demodulator probe not only becomes smooth, but also can keep the original form, which is helpful for improving the overall imaging quality of the three-dimensional low signal-to-noise ratio seismic data.

The residual static correction of the embodiment of the invention also comprises the residual static correction of the demodulator probe, and the invention calculates the first arrival residual static correction value of the demodulator probe besides calculating the first arrival residual static correction value of each shot point. As shown in fig. 2, the flow of calculating the first arrival remaining static correction amount of the demodulator probe according to the embodiment of the present invention includes steps S201 to S203.

Step S201, original first arrival travel time of the wave detection point corresponding to each shot point is obtained from the wave detection point gather data.

In the embodiment of the invention, the seismic data can be arranged according to the wave detection points to form wave detection point gather data, and then the original first-arrival travel time of the wave detection points of all seismic channels is read on the wave detection point gather data, so that the original first-arrival travel time of the wave detection points corresponding to each shot point is obtained.

And S202, fitting the original first-arrival travel time of the demodulator probe corresponding to each shot point by adopting a surface fitting method according to the position of each shot point relative to the demodulator probe to obtain the fitting first-arrival time of the demodulator probe corresponding to each shot point.

In the embodiment of the present invention, the method for fitting a curved surface in this step may be the same as that in the embodiment of step S102, and the specific fitting process may refer to the embodiment of step S102, which is not described herein again.

In an embodiment of the present invention, this step may calculate the fitting first arrival time of the demodulator probe corresponding to each shot point by using equations (1) to (11) of the embodiment of step S102 described above.

Step S203, calculating the difference value between the fitting first arrival time of the demodulator probe corresponding to each shot point and the original first arrival travel time, and calculating the average value of the difference value of the demodulator probe corresponding to each shot point to obtain the first arrival residual static correction value of the demodulator probe.

In this embodiment of the present invention, the step may be to apply the detection point R to the data of the detection point gather_jSubtracting the original first arrival time from the fitting first arrival time corresponding to each shot point, and averaging all the shot points to obtain a demodulator probe R_jResidual static correction value RT_jThe specific formula can be:

wherein N is the sum-detection point R_jNumber of corresponding shot points, T'_ijIs a detection point R_jAnd fitting first arrival time corresponding to each shot point.

In the embodiment of the invention, the obtained first arrival residual static correction value of the demodulator probe and the obtained first arrival residual static correction value of the shot point can be applied to the correction of various seismic data and can be used for improving the integral imaging quality of the three-dimensional low signal-to-noise ratio seismic data.

According to the embodiment of the invention, the first arrival residual static correction value of the shot point and the first arrival residual static correction value of the demodulator probe are calculated by a three-dimensional data surface fitting method, and the first arrival residual static correction value can be used for correcting seismic data.

FIG. 12 is a schematic cross-sectional view of a 3D stack of an embodiment of the invention subjected to only tomographic static correction processing, and FIG. 13 is a schematic cross-sectional view of a 3D stack of an embodiment of the invention subjected to tomographic static correction processing and residual static correction processing. As can be seen from fig. 12 and 13, the 3D stacking section after adding the first-arrival residual static correction values of the shot point and the demodulator probe calculated by the present application is improved in terms of signal-to-noise ratio, resolution and continuity of the in-phase axis compared with the 3D stacking section only subjected to tomographic static correction. Especially in the superficial layer (500ms-1500ms), some of the unidentified in-phase axes become quite distinct, and in addition, the in-phase axes in both the middle and deep layers are improved correspondingly.

Fig. 3 is a first flowchart of fitting the original first arrival travel time of the demodulator probe according to the embodiment of the present invention, and as shown in fig. 3, in the embodiment of the present invention, the fitting of the original first arrival travel time of each demodulator probe by using a surface fitting method according to the position of each demodulator probe relative to the shot point in step S102 may specifically include step S301 and step S302.

Step S301, a plane coordinate system is established by taking the shot as an origin to divide the shot gather data into four quadrants, wherein the X-axis direction of the plane coordinate system is consistent with the line detection direction.

In the embodiment of the present invention, a plane coordinate system is established in shot gather data with the shot as a center, and the shot gather data is divided into four quadrants, as shown in fig. 9, where the receiving points in fig. 9 are geophone points. As shown in fig. 9, the receiver points in the shot gather data are regularly arranged to form a plurality of detection lines, and in the embodiment of the present invention, the X-axis direction may be the same as the detection line direction when the plane coordinate system is established.

And step S302, fitting the original first arrival travel time of each demodulator probe in the quadrant by adopting a surface fitting method according to the coordinates of each demodulator probe in the quadrant aiming at each quadrant respectively to obtain the fitting result of each quadrant.

Fig. 4 is a second flowchart of fitting the original first arrival travel time of the demodulator probe according to the embodiment of the present invention, and as shown in fig. 4, the step S302 of fitting the original first arrival travel time of each demodulator probe in a quadrant by using a surface fitting method according to the coordinates of each demodulator probe in the quadrant in each quadrant respectively, may specifically include steps S401 to S403.

Step S401, calculating the offset distance between each detection point and the shot point according to the coordinates of each detection point.

In the embodiment of the present invention, the distance between the pickup point and the shot point, that is, the offset distance, can be calculated from the coordinates of the pickup point and the coordinates (origin) of the shot point.

Step S402, a three-dimensional coordinate system is established, and the coordinates of each wave detection point in the three-dimensional coordinate system are determined according to the offset distance of each wave detection point, the Y-axis coordinates of each wave detection point and the original first arrival travel time of each wave detection point.

In an alternative embodiment of the present invention, in the three-dimensional coordinate system, the offset distance of each demodulator probe and the Y-axis coordinate of each demodulator probe form a horizontal-axis coordinate, and the original first-arrival travel time of each demodulator probe is a vertical-axis coordinate, as shown in fig. 10. The original first arrival travel time data of each demodulator probe can be represented by points in a three-dimensional coordinate system.

And S403, fitting the original first arrival travel time of each demodulator probe by adopting a surface fitting method according to the coordinates of each demodulator probe in the three-dimensional coordinate system.

In an optional embodiment of the invention, a two-dimensional scattered point BSpline curved surface fitting method based on a thin plate can be adopted for fitting, and the fitting first-arrival time of each wave detection point is obtained.

Fig. 5 is a first flowchart of fitting the wave detection points to the original first-arrival travel time corresponding to each shot point according to an embodiment of the present invention, and as shown in fig. 5, in the embodiment of the present invention, the fitting of the original first-arrival travel time corresponding to each shot point of the wave detection points to each shot point by using a surface fitting method according to the position of each shot point relative to the wave detection point in step S202 may specifically include step S501 and step S502.

Step S501, a plane coordinate system is established by taking the demodulator probe as an origin, and the demodulator probe gather data is divided into four quadrants, wherein the X-axis direction of the plane coordinate system is consistent with the direction of the shot line.

In the embodiment of the invention, a plane coordinate system is established in the data of the wave detection point gather by taking the wave detection point as a center, and the data of the wave detection point gather is divided into four quadrants. In the embodiment of the invention, when a plane coordinate system is established, the X-axis direction can be the same as the direction of the shot lines.

And step S502, fitting the original first arrival travel time of the demodulator probe corresponding to each shot point by adopting a surface fitting method according to the coordinates of each shot point in each quadrant respectively to obtain the fitting result of each quadrant.

Fig. 6 is a second flowchart of fitting the wave detection point to the original first-arrival travel time corresponding to each shot point according to the embodiment of the present invention, and as shown in fig. 6, in the embodiment of the present invention, the step S502 of fitting the original first-arrival travel time corresponding to each shot point of the wave detection point to each quadrant by using a surface fitting method according to the coordinates of each shot point in the quadrant may specifically include step S601 and step S602.

Step S601, calculating the offset distance between the detection point and each shot point according to the coordinates of each shot point.

Step S602, a three-dimensional coordinate system is established, and the coordinates of each shot point in the three-dimensional coordinate system are determined according to the offset distance of each shot point, the Y-axis coordinate of each shot point and the original first arrival travel time of the wave detection point corresponding to each shot point.

In an alternative embodiment of the present invention, in the three-dimensional coordinate system, the offset distance of each shot and the Y-axis coordinate of each shot form a horizontal axis coordinate, and the original first arrival travel time of each shot is a vertical axis coordinate, as shown in fig. 10. The original first arrival travel time data for each shot may be represented as a point in a three-dimensional coordinate system.

And step S603, fitting the original first arrival travel time of the demodulator probe corresponding to each shot point by adopting a surface fitting method according to the coordinates of each shot point in the three-dimensional coordinate system.

The method comprises the steps of calculating the first-arrival residual static correction value of the shot point and the first-arrival residual static correction value of the demodulator probe by using a three-dimensional data surface fitting method, wherein the first-arrival residual static correction value is used for correcting seismic data. In addition, by the curved surface fitting method of the embodiment of the invention, the original first arrival travel time data of the detection point is smooth, and the original form can be kept, thereby being beneficial to improving the integral imaging quality of the three-dimensional low signal-to-noise ratio seismic data.

It should be noted that the steps illustrated in the flowcharts of the figures may be performed in a computer system such as a set of computer-executable instructions and that, although a logical order is illustrated in the flowcharts, in some cases, the steps illustrated or described may be performed in an order different than presented herein.

Based on the same inventive concept, the embodiment of the present invention further provides a three-dimensional first-arrival wave residual static correction apparatus, which can be used to implement the three-dimensional first-arrival wave residual static correction method described in the foregoing embodiment, as described in the following embodiment. Because the principle of solving the problem of the three-dimensional first-arrival wave residual static correction device is similar to that of the three-dimensional first-arrival wave residual static correction method, the embodiment of the three-dimensional first-arrival wave residual static correction device can be referred to as the embodiment of the three-dimensional first-arrival wave residual static correction method, and repeated parts are not repeated. As used hereinafter, the term "unit" or "module" may be a combination of software and/or hardware that implements a predetermined function. Although the means described in the embodiments below are preferably implemented in software, an implementation in hardware, or a combination of software and hardware is also possible and contemplated.

Fig. 7 is a first structural block diagram of a three-dimensional first-arrival wave residual static correction device according to an embodiment of the present invention, and as shown in fig. 7, the three-dimensional first-arrival wave residual static correction device according to the embodiment of the present invention includes: the device comprises a first original first arrival travel time acquisition unit 1, a first fitting first arrival time determination unit 2 and a shot point first arrival residual static correction value calculation unit 3.

A first original first-arrival travel time obtaining unit 1, configured to obtain, from the shot gather data, an original first-arrival travel time corresponding to the shot of each probe.

And the first fitting first arrival time determining unit 2 is used for fitting the original first arrival travel time of each demodulator probe by adopting a surface fitting method according to the position of each demodulator probe relative to the shot point to obtain the fitting first arrival time of each demodulator probe.

And the shot point first arrival residual static correction value calculating unit 3 is used for calculating the difference value between the fitting first arrival time of each demodulator probe and the original first arrival travel time, and calculating the average value of the difference values of all the demodulator probes to obtain the shot point first arrival residual static correction value.

Fig. 8 is a second structural block diagram of the three-dimensional first-arrival wave residual static correction device according to the embodiment of the present invention, and as shown in fig. 8, the three-dimensional first-arrival wave residual static correction device according to the embodiment of the present invention further includes: a second original first arrival travel time acquisition unit 4, a second fitting first arrival time determination unit 5 and a demodulator probe first arrival residual static correction value calculation unit 6.

And a second original first-arrival travel time obtaining unit 4, configured to obtain original first-arrival travel time of the probe corresponding to each shot point from the probe gather data.

And the second fitting first arrival time determining unit 5 is used for fitting the original first arrival travel time of each shot point corresponding to the corresponding wave detection point of the wave detection point by adopting a surface fitting method according to the position of each shot point relative to the wave detection point to obtain the fitting first arrival time of each shot point corresponding to the corresponding wave detection point.

And the demodulator probe first arrival residual static correction value calculating unit 6 is used for calculating the difference value between the fitting first arrival time of the demodulator probe corresponding to each shot point and the original first arrival travel time, and calculating the average value of the difference value of the demodulator probe corresponding to each shot point to obtain the first arrival residual static correction value of the demodulator probe.

In an embodiment of the present invention, the first fitted first arrival time determining unit 2 includes:

the first plane coordinate system establishing module is used for establishing a plane coordinate system by taking the shot as an origin to divide the shot gather data into four quadrants, wherein the X-axis direction of the plane coordinate system is consistent with the line detection direction;

and the first fitting module is used for fitting the original first arrival travel time of each demodulator probe in the quadrant by adopting a surface fitting method according to the coordinates of each demodulator probe in the quadrant aiming at each quadrant respectively to obtain the fitting result of each quadrant.

In an embodiment of the invention, the first fitting module comprises:

the first offset distance calculation submodule is used for calculating the offset distance between each detection point and the shot point according to the coordinates of each detection point;

the first three-dimensional coordinate system establishing submodule is used for establishing a three-dimensional coordinate system and determining the coordinates of each wave detection point in the three-dimensional coordinate system according to the offset distance of each wave detection point, the Y-axis coordinate of each wave detection point and the original first arrival travel time of each wave detection point;

and the first surface fitting submodule is used for fitting the original first arrival travel time of each wave detection point by adopting a surface fitting method according to the coordinate of each wave detection point in the three-dimensional coordinate system.

In an embodiment of the present invention, the second fitted first arrival time determining unit 5 includes:

the second plane coordinate system establishing module is used for establishing a plane coordinate system by taking the demodulator probe as an origin to divide the demodulator probe gather data into four quadrants, wherein the X-axis direction of the plane coordinate system is consistent with the direction of the gun line;

and the second fitting module is used for fitting the original first arrival travel time of the demodulator probe corresponding to each shot point by adopting a surface fitting method according to the coordinates of each shot point in each quadrant respectively to obtain the fitting result of each quadrant.

In an embodiment of the invention, the second fitting module comprises:

the second offset distance calculation submodule is used for calculating the offset distance between the detection point and each shot point according to the coordinates of each shot point;

the second three-dimensional coordinate system establishing submodule is used for establishing a three-dimensional coordinate system and determining the coordinates of each shot point in the three-dimensional coordinate system according to the offset distance of each shot point, the Y-axis coordinate of each shot point and the original first arrival travel time of the wave detection point corresponding to each shot point;

and the second surface fitting submodule is used for fitting the original first arrival travel time of the wave detection point corresponding to each shot point by adopting a surface fitting method according to the coordinate of each shot point in the three-dimensional coordinate system.

In an embodiment of the present invention, the surface fitting method includes: a two-dimensional scattered point BSpline curved surface fitting method based on a thin plate.

To achieve the above object, according to another aspect of the present application, there is also provided a computer apparatus. As shown in fig. 14, the computer device comprises a memory, a processor, a communication interface and a communication bus, wherein a computer program that can be run on the processor is stored in the memory, and the steps of the method of the embodiment are realized when the processor executes the computer program.

The processor may be a Central Processing Unit (CPU). The Processor may also be other general purpose processors, Digital Signal Processors (DSPs), Application Specific Integrated Circuits (ASICs), Field Programmable Gate Arrays (FPGAs) or other Programmable logic devices, discrete Gate or transistor logic devices, discrete hardware components, or a combination thereof.

The memory, which is a non-transitory computer readable storage medium, may be used to store non-transitory software programs, non-transitory computer executable programs, and units, such as the corresponding program units in the above-described method embodiments of the present invention. The processor executes various functional applications of the processor and the processing of the work data by executing the non-transitory software programs, instructions and modules stored in the memory, that is, the method in the above method embodiment is realized.

The memory may include a storage program area and a storage data area, wherein the storage program area may store an operating system, an application program required for at least one function; the storage data area may store data created by the processor, and the like. Further, the memory may include high speed random access memory, and may also include non-transitory memory, such as at least one disk storage device, flash memory device, or other non-transitory solid state storage device. In some embodiments, the memory optionally includes memory located remotely from the processor, and such remote memory may be coupled to the processor via a network. Examples of such networks include, but are not limited to, the internet, intranets, local area networks, mobile communication networks, and combinations thereof.

The one or more units are stored in the memory and when executed by the processor perform the method of the above embodiments.

The specific details of the computer device may be understood by referring to the corresponding related descriptions and effects in the above embodiments, and are not described herein again.

In order to achieve the above object, according to another aspect of the present application, there is also provided a computer-readable storage medium storing a computer program which, when executed in a computer processor, implements the steps in the above three-dimensional first arrival residual static correction method. It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by a computer program, which can be stored in a computer-readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. The storage medium may be a magnetic Disk, an optical Disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a Flash Memory (Flash Memory), a Hard Disk (Hard Disk Drive, abbreviated as HDD) or a Solid State Drive (SSD), etc.; the storage medium may also comprise a combination of memories of the kind described above.

It will be apparent to those skilled in the art that the modules or steps of the present invention described above may be implemented by a general purpose computing device, they may be centralized on a single computing device or distributed across a network of multiple computing devices, and they may alternatively be implemented by program code executable by a computing device, such that they may be stored in a storage device and executed by a computing device, or fabricated separately as individual integrated circuit modules, or fabricated as a single integrated circuit module from multiple modules or steps. Thus, the present invention is not limited to any specific combination of hardware and software.

The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.