Method and device for determining residual static correction value of first-motion wave
阅读说明：本技术 初至波剩余静校正量的确定方法及装置 (Method and device for determining residual static correction value of first-motion wave ) 是由 闫智慧 何永清 王珂 赵薇薇 高子涵 贾恒悦 于 2019-02-27 设计创作，主要内容包括：本申请公开了一种初至波剩余静校正量的确定方法及装置,该方法包括：获取目标区域中炮检对的初至信息；按照激发点和接收点连线的中心点划分炮检对到至少一个集合中；按照炮检距的大小,将每个集合划分为至少两个子集合；叠加基础静校正量和初至时间,得到每个炮检对的叠加值；将每个子集合中炮检对的叠加值进行线性拟合,得到拟合时间；根据拟合时间和炮检对的叠加值,确定每个炮检对的剩余静校正量；根据每个炮检对的剩余静校正量,确定每个炮检对中激发点和接收点的剩余静校正量；根据相同激发点和相同接收点在不同集合中计算得到的剩余静校正量和基础静校正量,确定初至波的剩余静校正量。本申请可以提高短波长初至波静校正量的精度。(The application discloses a method and a device for determining a residual static correction value of a first-motion wave, wherein the method comprises the following steps: acquiring first arrival information of shot-blast inspection pairs in a target area; dividing shot checking pairs into at least one set according to a central point of a connecting line of the excitation point and the receiving point; dividing each set into at least two subsets according to the size of the offset; superposing the basic static correction value and the first arrival time to obtain a superposed value of each shot detection pair; performing linear fitting on the superposition value of the shot-geophone pair in each subset to obtain fitting time; determining the residual static correction value of each shot-checking pair according to the fitting time and the superposition value of the shot-checking pair; determining the residual static correction values of the excitation point and the receiving point in each shot-check pair according to the residual static correction values of each shot-check pair; and determining the residual static correction value of the first arrival wave according to the residual static correction value and the basic static correction value which are obtained by calculating the same excitation point and the same receiving point in different sets. The method and the device can improve the precision of the short-wavelength first-arrival static correction value.)
1. A method for determining a residual static correction amount of a first-arrival wave, the method comprising:
acquiring first-arrival information of shot-geophone pairs in a target area, wherein the first-arrival information comprises a basic static correction value, first-arrival time and shot-geophone offset;
dividing shot-geophone pairs with the central points of connecting lines of the excitation points and the receiving points within a specified range into the same set to obtain at least one set; dividing each set into at least two subsets according to the size of the offset;
superposing the basic static correction value and the first arrival time of each shot-checking pair to obtain a superposed value corresponding to each shot-checking pair;
respectively carrying out linear fitting on the superposed values of the shot-geophone pairs in each subset to obtain fitting time changing along with shot-geophone distances;
respectively determining the residual static correction value of each shot-checking pair according to the fitting time and the superposition value of the shot-checking pair;
determining the residual static correction value of an excitation point and the residual static correction value of a receiving point in each shot-check pair according to the residual static correction value of each shot-check pair;
and determining the residual static correction value of the first arrival wave according to the residual static correction value calculated in different sets by the same excitation point and the same receiving point and the basic static correction value.
2. The method of claim 1, wherein after dividing each set into at least two subsets by offset size, the method further comprises:
judging that the number of shot-check pairs in the subset is greater than a number threshold;
the linear fitting is respectively carried out on the superposed values of the shot-geophone pairs in each subset to obtain the fitting time changing along with the shot-geophone distance, and the method comprises the following steps:
and if the number of the shot-geophone pairs in the subsets is larger than the number threshold, respectively carrying out linear fitting on the superposed values of the shot-geophone pairs in each subset to obtain the fitting time which changes along with the shot-geophone distance.
3. The method of claim 2, wherein determining the remaining static correction for each shot pair separately from the fitting time and the superposition of the shot pairs comprises:
according to the formula dt ═ t2-t1Determining the residual static correction dt of the shot-check pairs in each subset;
wherein, t2As fitting time, t2K and b are fitting parameters obtained by linear fitting, and x is offset; t is t1The superposition value of the shot-check pair.
4. The method of claim 1, wherein determining the remaining static correction for the excitation point and the remaining static correction for the receive point for each shot pair based on the remaining static correction for each shot pair comprises:
and determining the residual static correction value of half of the shot detection pair as the residual static correction value of the excitation point, and determining the residual static correction value of half of the shot detection pair as the residual static correction value of the receiving point.
5. The method according to any one of claims 1 to 4, wherein determining the residual static correction of the first arrival wave based on the residual static corrections calculated in different sets from the same excitation point and the same reception point and the base static correction comprises:
according to the formulaCalculating an average residual static correction value ts of each excitation point;
according to the formula
calculating the sum of the average residual static correction value of the excitation point, the average residual static correction value of the receiving point and the basic static correction value in the same shot-checking pair;
if the sum meets the preset condition, the sum is used as the residual static correction value of the primary wave received by the receiving point in the shot detection pair;
wherein Ns is the total number of residual static correction values calculated by the same excitation point in all sets;ith calculated for the same excitation point1Value of residual static correction amount i11,2, … Ns; nr is the total number of the residual static correction values calculated in all the sets by the same receiving point;i-th calculated for the same receiving point2Value of residual static correction amount i2＝1,2,…Nr。
6. The method according to claim 5, wherein the preset conditions include:
the root mean square of the residual static correction values of all shot-check pairs in all the sets is smaller than an error threshold;
or, after iteration for a specified number of times, calculating to obtain a residual static correction value convergence of the shot-geophone pair;
wherein, iteration is carried out according to the following method:
after the sum of the average residual static correction value of the excitation point, the average residual static correction value of the receiving point and the basic static correction value in the same shot-geophone pair is calculated, the sum is used as the basic static correction value of the shot-geophone pair, and the residual static correction value of the first arrival wave is calculated;
wherein according to the formula
7. An apparatus for determining a residual static correction amount of a first-arrival wave, the apparatus comprising:
the acquisition module is used for acquiring first arrival information of shot-geophone pairs in the target area, wherein the first arrival information comprises a basic static correction value, first arrival time and shot-geophone offset;
the dividing module is used for dividing the shot-geophone pair with the central point of the connecting line of the excitation point and the receiving point acquired by the acquiring module within a specified range into the same set to obtain at least one set; dividing each set into at least two subsets according to the size of the offset;
the superposition module is used for superposing the basic static correction value and the first arrival time of each shot detection pair acquired by the acquisition module to obtain a superposition value corresponding to each shot detection pair;
the fitting module is used for respectively performing linear fitting on the superposed values of the shot-geophone pairs in each subset obtained by superposing the superposition modules to obtain the fitting time changing along with the shot-geophone distance;
the determining module is used for respectively determining the residual static correction value of each shot detection pair according to the fitting time determined by the fitting module and the superposition value of the shot detection pair obtained by superposition of the superposition module;
the determining module is further used for determining the residual static correction value of the excitation point and the residual static correction value of the receiving point in each shot-check pair according to the residual static correction value of each shot-check pair;
and the determining module is further used for determining the residual static correction value of the first-motion wave according to the residual static correction value calculated by the same excitation point and the same receiving point in different sets and the basic static correction value.
8. The apparatus of claim 7, further comprising:
the judging module is used for judging that the number of the shot-geophone pairs in the subset divided by the dividing module is greater than a number threshold value;
a fitting module to:
and when the judging module judges that the number of the shot-geophone pairs in the subsets is larger than the number threshold, respectively carrying out linear fitting on the superposed values of the shot-geophone pairs in each subset to obtain the fitting time which changes along with the shot-geophone distance.
9. The apparatus of claim 8, wherein the means for determining is configured to:
according to the formula dt ═ t2-t1Determining the residual static correction dt of the shot-check pairs in each subset;
wherein, t2As fitting time, t2K and b are fitting parameters obtained by linear fitting, and x is offset; t is t1The superposition value of the shot-check pair.
10. The apparatus of claim 7, wherein the means for determining is configured to:
and determining the residual static correction value of half of the shot detection pair as the residual static correction value of the excitation point, and determining the residual static correction value of half of the shot detection pair as the residual static correction value of the receiving point.
11. The apparatus according to any one of claims 7 to 10, wherein the determining means is configured to:
according to the formulaCalculating an average residual static correction value ts of each excitation point;
according to the formula
calculating the sum of the average residual static correction value of the excitation point, the average residual static correction value of the receiving point and the basic static correction value in the same shot-checking pair;
if the sum meets the preset condition, the sum is used as the residual static correction value of the primary wave received by the receiving point in the shot detection pair;
wherein Ns is the total number of residual static correction values calculated by the same excitation point in all sets;
12. The apparatus of claim 11, wherein the preset condition comprises:
the root mean square of the residual static correction values of all shot-check pairs in all the sets is smaller than an error threshold;
or, after iteration for a specified number of times, calculating to obtain a residual static correction value convergence of the shot-geophone pair;
wherein, iteration is carried out according to the following method:
after the sum of the average residual static correction value of the excitation point, the average residual static correction value of the receiving point and the basic static correction value in the same shot-geophone pair is calculated, the sum is used as the basic static correction value of the shot-geophone pair, and the residual static correction value of the first arrival wave is calculated;
wherein according to the formulaCalculating root mean square rms; n is a radical of0For the total number of all shot pairs in all sets,is the ith3Value of residual static correction value, i, of individual shot-check pairs3＝1,2,…N0。
13. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method of any one of claims 1 to 6 when executing the computer program.
14. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program for executing the method of any one of claims 1 to 6.
Technical Field
The application relates to the technical field of geophysical exploration, in particular to a method and a device for determining a residual static correction value of a first-motion wave.
Background
In seismic exploration, static correction processing is mainly used for solving the problem that reflected wave superposition imaging is influenced due to travel time change when seismic waves travel in a near-surface medium. The first-arrival wave reflects the change of the near-surface stratum, a relatively accurate near-surface model can be inverted by using a first-arrival wave refraction static correction method or a chromatography static correction method, most of the static correction problems of long wavelength and medium and long wavelength can be solved by using the near-surface model, but in the aspect of solving the static correction problem of short wavelength, the error of the residual static correction value of the first-arrival wave of short wavelength determined by using the conventional first-arrival wave static correction method is relatively large, and the precision of the static correction of the first-arrival wave of short wavelength is relatively low. However, the short-wavelength static correction directly affects the imaging quality of the seismic section, and therefore, how to improve the accuracy of the short-wavelength first-arrival static correction and thus the imaging quality of the seismic section becomes a problem to be solved urgently at present.
Disclosure of Invention
The embodiment of the application provides a method for determining a first-arrival wave residual static correction value, which is used for improving the precision of short-wavelength first-arrival wave static correction, and comprises the following steps:
acquiring first-arrival information of shot-geophone pairs in a target area, wherein the first-arrival information comprises a basic static correction value, first-arrival time and shot-geophone offset; dividing shot-geophone pairs with the central points of connecting lines of the excitation points and the receiving points within a specified range into the same set to obtain at least one set; dividing each set into at least two subsets according to the size of the offset; superposing the basic static correction value and the first arrival time of each shot-checking pair to obtain a superposed value corresponding to each shot-checking pair; respectively carrying out linear fitting on the superposed values of the shot-geophone pairs in each subset to obtain fitting time changing along with shot-geophone distances; respectively determining the residual static correction value of each shot-checking pair according to the fitting time and the superposition value of the shot-checking pair; determining the residual static correction value of an excitation point and the residual static correction value of a receiving point in each shot-check pair according to the residual static correction value of each shot-check pair; and determining the residual static correction value of the first arrival wave according to the residual static correction value calculated in different sets by the same excitation point and the same receiving point and the basic static correction value.
The embodiment of the present application further provides a device for determining a first-arrival wave residual static correction value, so as to improve the precision of a short-wavelength first-arrival wave static correction value, where the device includes:
the acquisition module is used for acquiring first arrival information of shot-geophone pairs in the target area, wherein the first arrival information comprises a basic static correction value, first arrival time and shot-geophone offset; the dividing module is used for dividing the shot-geophone pair with the central point of the connecting line of the excitation point and the receiving point acquired by the acquiring module within a specified range into the same set to obtain at least one set; dividing each set into at least two subsets according to the size of the offset; the superposition module is used for superposing the basic static correction value and the first arrival time of each shot detection pair acquired by the acquisition module to obtain a superposition value corresponding to each shot detection pair; the fitting module is used for respectively performing linear fitting on the superposed values of the shot-geophone pairs in each subset obtained by superposing the superposition modules to obtain the fitting time changing along with the shot-geophone distance; the determining module is used for respectively determining the residual static correction value of each shot detection pair according to the fitting time determined by the fitting module and the superposition value of the shot detection pair obtained by superposition of the superposition module; the determining module is further used for determining the residual static correction value of the excitation point and the residual static correction value of the receiving point in each shot-check pair according to the residual static correction value of each shot-check pair; and the determining module is further used for determining the residual static correction value of the first-motion wave according to the residual static correction value calculated by the same excitation point and the same receiving point in different sets and the basic static correction value.
In the embodiment of the application, the shot-geophone pairs are divided into different sets according to the range of the central point of the connecting line of the excitation point and the receiving point to obtain a common central point set, and then the common central point set is further divided according to the shot-geophone distance to obtain a subset of each set. Then, linear fitting is carried out on the superposition values of the basic static correction values of the excitation points and the receiving points and the first arrival time in different subsets to obtain corrected residual static correction values. Therefore, on the basis of the basic static correction value of the original first-arrival wave, the parameter value of the calculated residual static correction value is corrected in a linear fitting mode, so that the accuracy of the calculated residual static correction value is improved, the accuracy of the short-wavelength first-arrival wave static correction is improved, and the imaging quality of the seismic section is improved; and the common central point set is divided into a plurality of subsets for linear fitting, so that the accuracy of the obtained fitting time is improved, and the accuracy of the calculated residual static correction value is improved.
Drawings
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present application, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts. In the drawings:
fig. 1 is a flowchart of a method for determining a residual static correction value of a first-arrival wave according to an embodiment of the present disclosure;
fig. 2 is a schematic diagram of a mesh partitioned according to mesh parameters according to an embodiment of the present application;
fig. 3 is a block diagram of an apparatus for determining a residual static correction amount of a first-arrival wave according to an embodiment of the present application.
Detailed Description
To make the objects, technical solutions and advantages of the embodiments of the present application more apparent, the embodiments of the present application are further described in detail below with reference to the accompanying drawings. The exemplary embodiments and descriptions of the present application are provided herein to explain the present application and not to limit the present application.
The embodiment of the application provides a method for determining a residual static correction value of a first-arrival wave, as shown in fig. 1, the method includes steps 101 to 107:
step 101, first arrival information of shot-geophone pairs in a target area is obtained.
The first arrival information comprises a basic static correction value, first arrival time and offset.
The target area is an area to be subjected to first arrival static correction.
It should be noted that the basic static correction amount includes a basic static correction amount of the excitation point and a basic static correction amount of the receiving point, which may be determined according to the prior art and will not be described herein. The first arrival time and the offset are determined according to actual measurement.
If there is an abnormality in the first arrival time, such as the first arrival time is too long or too short, the first arrival time of the abnormality needs to be eliminated.
102, dividing shot-geophone pairs with the central points of connecting lines of the excitation points and the receiving points within a specified range into the same set to obtain at least one set; and dividing each set into at least two subsets according to the size of the offset.
The specified range is a range divided by a user, the specified range can be one or more, and shot-geophone pairs with the excitation points and the receiving points in the same specified range are divided into the same set.
In the embodiment of the present application, the set obtained by the division is referred to as a first arrival gather G. The first arrival gather G is a common-midpoint domain gather, and the division process can be as follows:
firstly, grid parameters are determined. The grid parameters are generally set at a distance of 0.5 to 2 times the track pitch between two receiving points, for example, if the track pitch between two receiving points is 30 meters, the length of the side of the grid may be 15 meters, 30 meters, or the like.
And secondly, dividing the target area into one or more grids according to the grid parameters.
Specifically, referring to fig. 2, fig. 2 shows a mesh divided by a mesh parameter.
And thirdly, setting a coordinate system, calculating the coordinates of a central point of a connecting line between the excitation point and the receiving point by using the coordinates of the excitation point and the receiving point, and determining the grid where the central point is located.
In fig. 2, point S represents a trigger point, point R represents a receiver point, and after the coordinate system is set, each trigger point and receiver point have corresponding coordinates. On the basis of determining the coordinates of the excitation point and the receiving point, the coordinates of the center point of the excitation point and the center point of the corresponding receiving point can be obtained through calculation according to the coordinates, and the grid where the center point is located is determined according to the coordinates. Referring to fig. 2, the center point of the line connecting the excitation point S and the reception point R is in the grid G.
In the embodiment of the present application, one mesh is taken as one designated range.
And fourthly, respectively calculating the central point of a connecting line of the excitation point and the receiving point in each shot-examination pair, and dividing the shot-examination pairs with the central points in the same grid into the same set, so that a common central point first arrival gather G is formed in the grid.
It should be noted that, since each excitation point is excited, a large number of well-defined reception points are received. Thus, in general, there is a centerpoint in each mesh, that is, each mesh may form a common centerpoint gather.
The offset distances corresponding to the subsets can be manually set, for example, the range of all the offset distances in the set is 500-2000 meters, the step length of the offset distances can be set to 500, the offset distances are divided into 3 sections, namely 500-1000 meters, 1000-1500 meters and 1500-2000 meters, then, the offset pairs with the offset distances of 500-1000 meters are divided into the same subset, the offset pairs with the offset distances of 1000-1500 meters are divided into the same subset, and the offset pairs with the offset distances of 1500-2000 meters are divided into the same subset. The specific size of the offset for dividing the subset is not limited herein.
And 103, superposing the basic static correction value and the first arrival time of each shot detection pair to obtain a superposed value corresponding to each shot detection pair.
Alternatively, it can be represented by the formula t1＝t0+tcCalculating the corresponding superposition value t of each shot detection pair1. Wherein, t0Base static correction for excitation or reception points, tcThe first arrival time corresponding to the excitation point or the receiving point.
And step 104, respectively carrying out linear fitting on the superposed values of the shot-geophone pairs in each subset to obtain the fitting time changing along with the shot-geophone distance.
Optionally, because the amount of data used in the linear fitting is too small, the accuracy of the obtained fitting result is low, and therefore, after each set is divided into at least two subsets according to the offset, it can be further determined that the number of offset pairs in the subsets is greater than a number threshold; and if the number of the shot-geophone pairs in the subsets is larger than the number threshold, respectively carrying out linear fitting on the superposed values of the shot-geophone pairs in each subset to obtain the fitting time which changes along with the shot-geophone distance.
The number threshold may be set by human, and may be, for example, 3 or 5.
Alternatively, it can be represented by the formula t2Kx + b denotes the fitting function obtained after linear fitting, where t2Taking the fitting time as the reference value, k and b as fitting parameters obtained by linear fitting, and x as offset; t is t1Is a superimposed value.
And 105, respectively determining the residual static correction value of each shot-check pair according to the fitting time and the superposition value of the shot-check pairs.
Optionally, according to the formula dt ═ t2-t1The remaining static correction dt for the shot pairs in each subset is determined.
It should be noted that if dt is larger than the error range given by the user, the value is invalid and does not participate in the subsequent calculation.
And step 106, determining the residual static correction value of the excitation point and the residual static correction value of the receiving point in each shot-check pair according to the residual static correction value of each shot-check pair.
Alternatively, for the excitation point and the receiving point included in a shot pair, half of the remaining static corrections for the shot pair may be determined as remaining static corrections for the excitation point and the other half may be determined as remaining static corrections for the receiving point in the shot pair.
Alternatively, it can be expressed by formula
And step 107, determining the residual static correction value of the primary wave according to the residual static correction value calculated in different sets by the same excitation point and the same receiving point and the basic static correction value.
By repeating the steps 103 to 106, the remaining static correction values of the same excitation point or the same receiving point in different sets can be obtained.
Alternatively, it can be according to a formulaCalculating an average residual static correction value ts of each excitation point; according to the formulaCalculating an average residual static correction value tr of each receiving point; calculating the sum of the average residual static correction value of the excitation point, the average residual static correction value of the receiving point and the basic static correction value in the same shot-checking pair; and if the summation meets the preset condition, the summation is used as the residual static correction value of the primary wave received by the receiving point in the shot detection pair.
Wherein Ns is the total number of residual static correction values calculated by the same excitation point in all sets;
Optionally, the preset conditions include the following two conditions:
(1) the root mean square of the remaining static corrections for all shot pairs in all sets is less than the error threshold.
Wherein according to the formula
The error threshold may be set empirically by the user, and the specific value thereof is not limited herein.
(2) And after iteration for a specified number of times, calculating the residual static correction value of the shot-check pair to be converged. Wherein, iteration is carried out according to the following method: after the sum of the average residual static correction value of the excitation point, the average residual static correction value of the receiving point and the basic static correction value in the same shot-checking pair is calculated, the sum is used as the basic static correction value of the shot-checking pair, and the residual static correction value of the primary wave is calculated according to steps 103 to 107.
And if the sum meets the preset condition (2), taking the sum obtained in the last iteration as the residual static correction value of the first-motion wave.
It should be noted that if the calculated sum satisfies any one of the two preset conditions, the sum is used as the remaining static correction amount of the first-motion wave received by the receiving point in the shot detection pair.
In the embodiment of the application, the shot-geophone pairs are divided into different sets according to the range of the central point of the connecting line of the excitation point and the receiving point to obtain a common central point set, and then the common central point set is further divided according to the shot-geophone distance to obtain a subset of each set. Then, linear fitting is carried out on the superposition values of the basic static correction values of the excitation points and the receiving points and the first arrival time in different subsets to obtain corrected residual static correction values. Therefore, on the basis of the basic static correction value of the original first-arrival wave, the parameter value of the calculated residual static correction value is corrected in a linear fitting mode, so that the accuracy of the calculated residual static correction value is improved, the accuracy of the short-wavelength first-arrival wave static correction is improved, and the imaging quality of the seismic section is improved; and the common central point set is divided into a plurality of subsets for linear fitting, so that the accuracy of the obtained fitting time is improved, and the accuracy of the calculated residual static correction value is improved.
The embodiment of the present application further provides a device for determining a residual static correction amount of a first-arrival wave, as shown in fig. 3, the device 300 includes an obtaining module 301, a dividing module 302, a superimposing module 303, a fitting module 304, and a determining module 305.
The obtaining module 301 is configured to obtain first arrival information of a shot-geophone pair in a target region, where the first arrival information includes a basic static correction amount, first arrival time, and a shot-geophone distance.
A dividing module 302, configured to divide the shot-geophone pairs whose central points of the connecting lines of the excitation point and the receiving point obtained by the obtaining module 301 are within a specified range into the same set, so as to obtain at least one set; and dividing each set into at least two subsets according to the size of the offset.
And the superposition module 303 is configured to superpose the basic static correction amount and the first arrival time of each shot detection pair acquired by the acquisition module 301 to obtain a superposition value corresponding to each shot detection pair.
And the fitting module 304 is configured to perform linear fitting on the superposed values of the shot-geophone pairs in each subset obtained by superposing in the superposing module 303, so as to obtain fitting time varying with shot-geophone distances.
And a determining module 305, configured to determine a remaining static correction value of each shot detection pair according to the fitting time determined by the fitting module 304 and a superimposed value of the shot detection pair superimposed by the superimposing module 303.
The determining module 305 is further configured to determine a remaining static correction value of the excitation point and a remaining static correction value of the receiving point in each shot pair according to the remaining static correction value of each shot pair.
The determining module 305 is further configured to determine the remaining static correction value of the first-arrival wave according to the remaining static correction values calculated in different sets by using the same excitation point and the same receiving point, and the base static correction value.
Optionally, the apparatus 300 further comprises:
and the judging module 306 is configured to judge that the number of shot pairs in the subset divided by the dividing module 302 is greater than the number threshold.
A fitting module 304 for: when the judging module 306 judges that the number of the shot-geophone pairs in the subsets is greater than the number threshold, linear fitting is respectively performed on the superposed values of the shot-geophone pairs in each subset, so as to obtain the fitting time varying with the shot-geophone distance.
Optionally, the determining module 305 is configured to:
according to the formula dt ═ t2-t1Determining the residual static correction dt of the shot-check pairs in each subset;
wherein, t2As fitting time, t2K and b are fitting parameters obtained by linear fitting, and x is offset; t is t1The superposition value of the shot-check pair.
Optionally, the determining module 305 is configured to:
and determining the residual static correction value of half of the shot detection pair as the residual static correction value of the excitation point, and determining the residual static correction value of half of the shot detection pair as the residual static correction value of the receiving point.
Optionally, the determining module 305 is configured to:
according to the formulaCalculating an average residual static correction value ts of each excitation point;
according to the formula
calculating the sum of the average residual static correction value of the excitation point, the average residual static correction value of the receiving point and the basic static correction value in the same shot-checking pair;
if the summation meets the preset condition, the summation is used as the residual static correction value of the primary wave received by the receiving point in the shot detection pair;
wherein Ns is the total number of residual static correction values calculated by the same excitation point in all sets;ith calculated for the same excitation point1Value of residual static correction amount i11,2, … Ns; nr is the total number of the residual static correction values calculated in all the sets by the same receiving point;
Optionally, the preset conditions include:
the root mean square of the residual static correction values of all shot-check pairs in all the sets is smaller than an error threshold;
or, after iteration for a specified number of times, calculating to obtain a residual static correction value convergence of the shot-geophone pair;
wherein, iteration is carried out according to the following method:
after the sum of the average residual static correction value of the excitation point, the average residual static correction value of the receiving point and the basic static correction value in the same shot-geophone pair is calculated, the sum is used as the basic static correction value of the shot-geophone pair, and the residual static correction value of the first arrival wave is calculated;
wherein according to the formula
In the embodiment of the application, the shot-geophone pairs are divided into different sets according to the range of the central point of the connecting line of the excitation point and the receiving point to obtain a common central point set, and then the common central point set is further divided according to the shot-geophone distance to obtain a subset of each set. Then, linear fitting is carried out on the superposition values of the basic static correction values of the excitation points and the receiving points and the first arrival time in different subsets to obtain corrected residual static correction values. Therefore, on the basis of the basic static correction value of the original first-arrival wave, the parameter value of the calculated residual static correction value is corrected in a linear fitting mode, so that the accuracy of the calculated residual static correction value is improved, the accuracy of the short-wavelength first-arrival wave static correction is improved, and the imaging quality of the seismic section is improved; and the common central point set is divided into a plurality of subsets for linear fitting, so that the accuracy of the obtained fitting time is improved, and the accuracy of the calculated residual static correction value is improved.
The embodiment of the application also provides computer equipment which comprises a memory, a processor and a computer program which is stored on the memory and can run on the processor, wherein the processor realizes the determination method of the first-arrival-wave residual static correction value when executing the computer program.
An embodiment of the present application also provides a computer-readable storage medium storing a computer program for executing the method for determining the first-arrival-wave residual static correction amount.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application 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.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The above-mentioned embodiments are further described in detail for the purpose of illustrating the invention, and it should be understood that the above-mentioned embodiments are only illustrative of the present invention and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements, etc. made within the spirit and principle of the present invention should be included in the scope of the present invention.