Primary wave residual static correction method and device, storage medium and electronic equipment

Abstract

The invention provides a residual static correction method and device for first-motion waves, a storage medium and electronic equipment, and the method comprises the steps: constructing a time-distance graph according to the first-motion data corresponding to each first-motion wave, carrying out the meshing of the time-distance graph, and obtaining a grid sequence in the meshed time-distance graph; fitting based on the first arrival data to obtain a first arrival spline interpolation function; traversing the midpoint of each grid in the grid sequence, and constructing an initial node set based on the currently traversed grid and the traversed grid; based on the initial node set and the first arrival spline interpolation function, obtaining a target grid meeting a preset error evaluation strategy from the initial node set; fitting is carried out based on the obtained first arrival data corresponding to the midpoint of the target grid, and a first arrival fitting function is obtained; and based on the first arrival data corresponding to the common shot point, the first arrival data corresponding to the common detection point and the first arrival fitting function, obtaining the residual static correction value of the first arrival wave. The calculation precision can be improved.

Description

Technical Field

[0001] The present invention relates to the field of processing of petroleum geophysical exploration data, and particularly to a residual static correction method, device, storage medium and electronic device for first arrival waves. Background Art

[0002] Obtaining a residual static correction amount based on first arrival data corresponding to each first arrival wave, and performing static correction on seismic exploration data, such as first arrival data, based on the residual static correction amount, plays a crucial role in the process of seismic exploration data processing and directly affects the underground structure morphology and imaging quality. Therefore, improving the accuracy and precision of the residual static correction amount obtained by using the residual static correction method is the key to improving the imaging quality.

[0003] Currently, the residual static correction method is to use a mathematical statistics method to fit the offset and time in the first arrival data in the common midpoint (CMP) domain to obtain a fitting function, and then obtain the residual static correction amount based on the first arrival data corresponding to the common shot point, the first arrival data corresponding to the common receiver point, and the fitting function.

[0004] However, for seismic exploration data in the "double complex" area, due to the rapid change of the underground medium velocity in the "double complex" area, the first arrival data contains more noise and a large time difference. Using all the first arrival data in the "double complex" area for fitting will cause the fitting function not to cover most of the first arrival data, with a large dispersion and low precision. If the first arrival data within a partial offset range is selected for calculation, although the dispersion of the fitting function can be effectively reduced, the reliability of the obtained fitting function is reduced, and local underfitting or overfitting problems are likely to occur. After performing static correction on the first arrival data using the obtained residual static correction amount, serious "nonlinear" and long wavelength problems may occur, and the precision of the obtained residual static correction amount is not high. Summary of the Invention

[0005] In view of this, the present invention provides a residual static correction method, device, storage medium and electronic device for first arrival waves.

[0006] Specifically, the present invention is implemented by the following technical solutions:

[0007] According to a first aspect of the present invention, there is provided a residual static correction method for first arrival waves, the method including:

[0008] Constructing a traveltime graph based on first arrival data corresponding to each first arrival wave, performing grid division on the traveltime graph, and obtaining a grid sequence in the grid-divided traveltime graph;

[0009] Performing fitting based on the first arrival data to obtain a first arrival spline interpolation function;

[0010] Traverse the midpoints of each grid in the grid sequence, and construct an initial node set based on the currently traversed grid and the grids that have been traversed;

[0011] Based on the initial node set and the first-arrival spline interpolation function, obtain the target grid that meets the pre-set error evaluation strategy from the initial node set;

[0012] Perform fitting based on the first-arrival data corresponding to the midpoint of the obtained target grid to obtain the first-arrival fitting function;

[0013] Based on the first-arrival data corresponding to the common shot point, the first-arrival data corresponding to the common receiver point, and the first-arrival fitting function, obtain the residual static correction amount of the first-arrival wave.

[0014] In the method for residual static correction of the first-arrival wave in this embodiment, the target grid is obtained by using the error evaluation strategy, the first-arrival data is optimized, the first-arrival data corresponding to the target grid is used for fitting, and a high-precision first-arrival fitting function is obtained. According to the first-arrival fitting function and the first-arrival data corresponding to the common shot point and the common shoulder wave point, the fitting time difference is calculated, and the residual static correction amount of the first-arrival wave is obtained, which has higher calculation accuracy.

[0015] According to the second aspect of the present invention, there is provided a device for residual static correction of the first-arrival wave, and the device for residual static correction of the first-arrival wave includes:

[0016] A grid sequence acquisition module, configured to construct a travel-time graph based on the first-arrival data corresponding to each first-arrival wave, perform grid division on the travel-time graph, and obtain the grid sequence in the grid-divided travel-time graph;

[0017] A spline interpolation function module, configured to perform fitting based on the first-arrival data to obtain a first-arrival spline interpolation function;

[0018] An initial node construction module, configured to traverse the midpoints of each grid in the grid sequence, and construct an initial node set based on the currently traversed grid and the grids that have been traversed;

[0019] An error evaluation module, configured to obtain a target grid that meets the pre-set error evaluation strategy from the initial node set based on the initial node set and the first-arrival spline interpolation function;

[0020] A fitting function acquisition module, configured to perform fitting based on the first-arrival data corresponding to the midpoint of the obtained target grid to obtain a first-arrival fitting function;

[0021] A residual static correction module, configured to obtain the residual static correction amount of the first-arrival wave based on the first-arrival data corresponding to the common shot point, the first-arrival data corresponding to the common receiver point, and the first-arrival fitting function.

[0022] According to a third aspect of the present invention, there is provided a storage medium having stored thereon a computer program, and when the program is executed by a processor, the steps of the residual static correction method for first arrival waves in any possible implementation manner of the first aspect are implemented.

[0023] According to a fourth aspect of the present invention, there is provided an electronic device, including a memory, a processor, and a computer program stored on the memory and executable on the processor, and when the processor executes the program, the steps of the residual static correction method for first arrival waves in any possible implementation manner of the first aspect are implemented. BRIEF DESCRIPTION OF THE DRAWINGS

[0024] The accompanying drawings herein are incorporated into and constitute a part of this specification, showing embodiments consistent with the present invention, and are used together with the specification to explain the principles of the present invention.

[0025] In order to more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the accompanying drawings required for use in the description of the embodiments or related technologies. Obviously, for those of ordinary skill in the art, other drawings can also be obtained based on these drawings without creative efforts.

[0026] Figure 1 It is a schematic flow chart of a residual static correction method for first arrival waves provided by an embodiment of the present invention;

[0027] Figure 2 It is a schematic diagram of a gridded travel time graph in a residual static correction method for first arrival waves provided by an embodiment of the present invention;

[0028] Figure 3 It is a schematic diagram of a residual static correction device for first arrival waves provided by an embodiment of the present invention;

[0029] Figure 4 It is a schematic structural diagram of an electronic device provided by an embodiment of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0030] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the following will clearly and completely describe the technical solutions in the embodiments of the present invention with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are some, but not all, of the embodiments of the present invention. All other embodiments obtained by those of ordinary skill in the art without creative efforts based on the embodiments of the present invention belong to the scope of protection of the present invention.

[0031] In the related art, due to the rapid change of the underground medium velocity, the first arrival data contains a lot of noise and a large time difference. If only the first arrival data within a partial offset range is selected for calculation, the reliability of the obtained fitting function is reduced, and local underfitting or overfitting problems are likely to occur. After static correction of the first arrival data using the obtained residual static correction amount, serious "nonlinear" and long wavelength problems may be caused, resulting in low calculation accuracy.

[0032] In this embodiment, the first arrival data is fitted by an improved spline curve, the target grid is obtained through an error evaluation strategy, the first arrival data is optimized, and the first arrival data corresponding to the target grid is used for fitting to obtain a high-precision first arrival fitting function. According to the first arrival fitting function and the corresponding first arrival data of the common shot point and the common shoulder wave point, the fitting time difference is calculated, and the residual static correction amount of the first arrival wave is obtained, so as to achieve the purpose of improving the fitting smoothness and accuracy.

[0033] See Figure 1 , an embodiment of the present invention provides a method for residual static correction of first arrival waves, which can be applied to perform residual static correction on the first arrival data corresponding to the first arrival waves. The method may include the following steps:

[0034] S101. Construct a travel time graph based on the first arrival data corresponding to each first arrival wave, grid the travel time graph, and obtain a grid sequence in the gridded travel time graph;

[0035] In this embodiment, as an optional embodiment, a travel time graph is drawn using the first arrival data corresponding to each first arrival wave, and each first arrival data point represents the arrival time of the seismic wave and the corresponding distance.

[0036] In this embodiment, as an optional embodiment, gridding the travel time graph to obtain a grid matrix in the gridded travel time graph includes:

[0037] Dividing the axis where the offset is located in the travel time graph according to a preset offset threshold to obtain a plurality of grids perpendicular to the axis where the time is located in the travel time graph;

[0038] For each grid, obtain the maximum time value and the minimum time value in the corresponding first arrival data within the grid, calculate the average value of the maximum time value and the minimum time value to obtain the time value of the midpoint of the grid, and calculate the average value of the maximum value and the minimum value of the grid on the axis where the offset is located to obtain the offset of the midpoint of the grid.

[0039] In this embodiment, as an optional embodiment, the offset refers to the distance between the shot point and the receiving point. As Figure 2 shown, the vertical columns are grids, and the offset of the midpoint of the grid is the average value of the initial coordinate and the final coordinate of the grid on the abscissa.

[0040] S102. Perform fitting based on the first arrival data to obtain a first arrival spline interpolation function;

[0041] In this embodiment, as an alternative embodiment, perform fitting based on the first arrival data to obtain a first arrival spline interpolation function. Using the first arrival spline interpolation function, the values between known data points can be estimated by interpolation between these points. In this embodiment, as an alternative embodiment, using spline interpolation of the first arrival spline interpolation function, an estimated value of the distance corresponding to this position can be obtained by inputting the time of the first arrival.

[0042] In this embodiment, as an alternative embodiment, the spline interpolation function can be expressed as:

[0043] S(x) = a i + b i (x - x i ) + c i (x - x i ) 2 + d i (x - x i ) 3 , x ∈ [x i , x i+1

[0044] where a i , b i , c i , d i are coefficients, x i is the offset in the first arrival data, and x ∈ [x i , x i+1 indicates interpolation between the first arrival data x i and the first arrival data x i+1 . In this embodiment, each piece of first arrival data corresponds to a node, and S(x) is the time value of the spline interpolation function.

[0045] S103. Traverse the midpoints of each grid in the grid sequence, and construct an initial node set based on the currently traversed grid and the grids that have been traversed;

[0046] In this embodiment, as an alternative embodiment, constructing an initial node set based on the currently traversed grid and the grids that have been traversed includes:

[0047] Extract the midpoint of the first grid in the grid sequence to construct a first node set;

[0048] Extract the midpoint of the first grid and the midpoint of the second grid in the grid sequence to construct a second node set;

[0049] Extract the midpoints of the first grid, the second grid, up to the nth grid in the grid sequence to construct the nth node set, where n is a natural number greater than 2, and the maximum value is the number of grids included in the grid matrix;

[0050] Based on the first node set, the second node set, and the nth node set, obtain the initial node set.

[0051] In this embodiment, as an alternative embodiment, taking the case where there are 4 grids after meshing as an example, the initial node set includes: the first node set composed of the 1st grid, the second node set composed of the 1st grid and the 2nd grid, the third node set composed of the 1st grid to the 3rd grid, and the fourth node set composed of the 1st grid to the 4th grid.

[0052] S104. Based on the initial node set and the first arrival spline interpolation function, obtain a target grid from the initial node set that meets a pre-set error evaluation strategy;

[0053] In this embodiment, as an alternative embodiment, the obtaining a target grid from the initial node set that meets a pre-set error evaluation strategy includes:

[0054] Based on the first arrival spline interpolation function, respectively obtain the error evaluation values of each node set in the initial node set;

[0055] Traverse each error evaluation value, calculate the difference between this error evaluation value and the previous adjacent error evaluation value, and select two node sets with the largest difference that meet the corresponding error evaluation strategy;

[0056] Calculate the difference set of the two node sets to obtain the target grid;

[0057] Delete the target grid from the node set containing the target grid, and delete the node set corresponding to the previous adjacent error evaluation value from the initial node set after deleting the target grid, and execute the step of respectively obtaining the error evaluation values of each node set in the initial node set based on the first arrival spline interpolation function.

[0058] In this embodiment, as an alternative embodiment, the error evaluation strategy includes: the difference is greater than a pre-set difference threshold, and / or the iteration round corresponding to the difference is greater than or equal to a pre-set round threshold, where one iteration round is to execute once the step of respectively obtaining the error evaluation values of each node set in the initial node set based on the first arrival spline interpolation function.

[0059] In this embodiment, as an alternative embodiment, during the process of calculating the error evaluation value, the number of iteration rounds should comprehensively consider the characteristics of the first arrival data, fitting accuracy, and calculation efficiency. As an alternative embodiment, the number of iteration rounds can be preset to 3 times.

[0060] In this embodiment, as an alternative embodiment, taking the initial node set including 5 node sets as an example:

[0061] Using the formula, calculate the error evaluation values of the first node set, the second node set, the third node set, and the fourth node set respectively, and then calculate the first difference between the error evaluation value of the second node set and the error evaluation value of the first node set, the second difference between the error evaluation value of the third node set and the error evaluation value of the second node set, the third difference between the error evaluation value of the fourth node set and the error evaluation value of the third node set, and the fourth difference between the error evaluation value of the fifth node set and the error evaluation value of the fourth node set. Assume that the third difference is the largest and greater than the preset difference threshold, then obtain the difference set between the fourth node set and the third node set, that is, the 4th grid is the target grid;

[0062] In the initial node set of the target grid, the second node set includes the first grid and the second grid, and the third node set also includes the first grid and the second grid. Delete the third node set, and the updated initial node set becomes including 4 node sets. Recalculate the error evaluation value of each node set according to the above method.

[0063] In this embodiment, as an alternative embodiment, the following formula can be used to calculate the error evaluation value:

[0064] E = αE approx + βE smooth

[0065] In the formula, α is the weight parameter of the approximation error, and β is the weight parameter of the smooth error.

[0066] E approx is the approximation error, expressed as:

[0067]

[0068] E smooth is the smooth error, expressed as:

[0069]

[0070] In the formula, y i is the first arrival data, S(x i ) is the time corresponding to the offset of the first arrival data on the interpolation function, and S″(x i ) is the second derivative of the interpolation function at the node x i .

[0071] In this embodiment, as an alternative embodiment, the approximation error weight parameter α and the smooth error weight parameter β in the error evaluation items can be adjusted according to the characteristics of the first arrival data and application requirements. For example, when the velocity changes drastically and there is a velocity inversion, α is taken as 0.6 and β is taken as 0.4; in general cases, α is taken as 0.8 and β is taken as 0.2.

[0072] S105. Fit based on the first arrival data corresponding to the midpoints of the obtained target grids to obtain a first arrival fitting function.

[0073] In this embodiment, as an alternative embodiment, based on the first arrival data corresponding to the midpoints of the obtained target grids, use the steps described in S102 to calculate and obtain a first arrival fitting function. For Figure 2 example, the vertical columns are grids. The offset of the midpoint of the grid is the average of the initial coordinate and the final coordinate of the grid on the abscissa, and the time of the midpoint of the grid is the average of the minimum time and the maximum event corresponding to each first arrival data within the grid.

[0074] S106. Based on the first arrival data corresponding to the common shot points, the first arrival data corresponding to the common geophone points, and the first arrival fitting function, obtain the residual static correction amount of the first arrival wave.

[0075] In this embodiment, the residual static correction amount of the first arrival wave includes: the residual static correction amount of the shot point and the residual static correction amount of the geophone point. For the shot point, use the residual static correction amount of the shot point to correct the first arrival data of the shot point. For the geophone point, use the residual static correction amount of the geophone point to correct the first arrival data of the geophone point. As an alternative embodiment, based on the first arrival data corresponding to the common shot points, the first arrival data corresponding to the common geophone points, and the first arrival fitting function, obtaining the residual static correction amount of the first arrival wave includes:

[0076] Obtain the first time of the common shot point in the first arrival fitting function, and the second time of the first arrival data of the common shot point. Based on the first time and the second time, obtain the residual static correction amount of the shot point in the first arrival data;

[0077] Obtain the third time of the common geophone point in the first arrival fitting function, and the fourth time of the first arrival data of the common geophone point. Based on the third time and the fourth time, obtain the residual static correction amount of the geophone point in the first arrival data.

[0078] In this embodiment, as an alternative embodiment, obtain the offset in the first arrival data of the common shot point, and obtain the time value corresponding to the offset in the first arrival fitting function. This time value is the first time value.

[0079] In this embodiment, as an alternative embodiment, correct the first arrival data based on the residual static correction amount.

[0080] In this embodiment, as an alternative embodiment, the first arrival data corresponding to the common shot point, the first arrival data corresponding to the common receiver point, and the first arrival fitting function are obtained respectively, and the time difference τ between the first arrival data corresponding to the midpoint of each target grid and the fitting curve is calculated respectively.

[0081] In this embodiment, as an alternative embodiment, for each common shot point, the fitting function is used to predict the expected first arrival time when they reach each common receiver point; for each pair of common shot points and common receiver points, the time difference between the actually observed first arrival time and the first arrival time predicted by the fitting function is calculated; for each common shot point, the time differences of all common receiver points are aggregated and the average value is calculated. This average time difference can be regarded as the residual static correction amount of this shot point.

[0082] In this embodiment, as an alternative embodiment, it is assumed that the time difference of the first arrival data of the i-th trace of the S-th shot point is τ i , then there is

[0083]

[0084] If there are m receiving points in the receiver array of the S-th shot point, and the shot point static correction amount τ S of each trace is the same, and the static correction amounts of each receiving point are respectively If the total time differences of all traces of receiving points are added and averaged, that is:

[0085]

[0086] In this embodiment, when analyzing the first arrival data through the fitting model, due to the complexity of the underground medium or other factors, there is a time difference between the actually observed first arrival data and the first arrival time predicted by the fitting function. The residual static correction amount reflects the differences in the first arrival time caused by factors such as changes in the underground structure that the model fails to capture and differences in instrument responses. By calculating the average time difference, an index reflecting the fitting effect on the entire data set can be obtained. Therefore, this average time difference can be regarded as the residual static correction amount.

[0087] In this embodiment, as an alternative embodiment, the time differences of each trace belonging to the same shot point S are added and averaged to obtain the static correction amount of this shot point, and thus the residual static correction amount τ S of this shot point S is obtained. Similarly, at the common receiver point, the residual static correction amount τ R of the receiver point R can be calculated.

[0088] In this embodiment, the residual static correction amount of the shot point is obtained based on common shot points, that is, based on the above formula, the residual static correction amount of the common shot point is obtained as the second time. Then, based on each initial data included in the common shot point, the time difference corresponding to the offset in the first arrival fitting function in each initial data is obtained, and its average value is calculated as the first time. The difference between the first time and the second time is calculated to obtain the residual static correction amount of the shot point.

[0089] In this embodiment, by constructing a travel-time graph according to the first arrival data corresponding to each first arrival wave, meshing the travel-time graph, and obtaining the grid sequence in the meshed travel-time graph; performing fitting based on the first arrival data to obtain a first arrival spline interpolation function; traversing the midpoint of each grid in the grid sequence, and constructing an initial node set based on the currently traversed grid and the traversed grids; based on the initial node set and the first arrival spline interpolation function, obtaining a target grid that meets a pre-set error evaluation strategy from the initial node set; performing fitting based on the first arrival data corresponding to the midpoint of the obtained target grid to obtain a first arrival fitting function; and obtaining the residual static correction amount of the first arrival wave based on the first arrival data corresponding to the common shot point, the first arrival data corresponding to the common receiver point, and the first arrival fitting function. In this way, the target grid is obtained by using the error evaluation strategy, the selected spline nodes are optimized, and the first arrival data corresponding to the target grid is used for fitting, so as to obtain a high-precision first arrival fitting function. According to the first arrival fitting function, the common shot point, and the first arrival data corresponding to the common shoulder wave point, the fitting time difference is calculated, and the residual static correction amount of the first arrival wave is obtained through residual static correction decomposition, which has higher first arrival data fitting accuracy and residual static correction accuracy.

[0090] Based on the same inventive concept, as Figure 3 shown, the embodiment of the present invention further provides a device for residual static correction of the first arrival wave. The device includes:

[0091] A grid sequence acquisition module 301, configured to construct a travel-time graph according to the first arrival data corresponding to each first arrival wave, mesh the travel-time graph, and obtain the grid sequence in the meshed travel-time graph;

[0092] In this embodiment, as an optional embodiment, a travel-time graph is drawn using the first arrival data corresponding to each first arrival wave, where each point represents the time when the seismic wave arrives and the corresponding distance, and the travel-time graph is meshed to obtain the grid matrix in the meshed travel-time graph.

[0093] In this embodiment, as an alternative embodiment, the axis where the offset is located in the traveltime diagram is divided according to a preset offset threshold to obtain a plurality of grids perpendicular to the axis where the time is located in the traveltime diagram; for each grid, the maximum time value and the minimum time value among the corresponding first arrival data within the grid are obtained, the average value of the maximum time value and the minimum time value is calculated to obtain the time value of the midpoint of the grid, and the average value of the maximum value and the minimum value of the grid on the axis where the offset is located is calculated to obtain the offset of the midpoint of the grid.

[0094] The spline interpolation function module 302 is configured to perform fitting based on the first arrival data to obtain a first arrival spline interpolation function;

[0095] In this embodiment, as an alternative embodiment, fitting is performed based on the first arrival data to obtain a first arrival spline interpolation function. Using the first arrival spline interpolation function, the values between known data points can be estimated by interpolation between these points. For example, using the first arrival spline interpolation function, the estimated value of the distance corresponding to this position can be obtained by inputting the arrival time of the first arrival.

[0096] The initial node construction module 303 is configured to traverse the midpoints of each grid in the grid sequence, and construct an initial node set based on the currently traversed grid and the grids that have been traversed;

[0097] In this embodiment, as an alternative embodiment, traverse the midpoints of each grid in the grid sequence, extract the midpoint of the first grid in the grid sequence to construct a first node set; extract the midpoints of the first grid and the second grid in the grid sequence to construct a second node set; extract the midpoints of the first grid, the second grid until the nth grid in the grid sequence to construct an nth node set, where n is a natural number greater than 2, and the maximum value is the number of grids included in the grid matrix; based on the first node set, the second node set, and the nth node set, obtain the initial node set.

[0098] The error evaluation module 304 is configured to obtain a target grid that meets a preset error evaluation strategy from the initial node set based on the initial node set and the first arrival spline interpolation function;

[0099] In this embodiment, as an alternative embodiment, based on the initial node set and the first arrival spline interpolation function, the initial node set is optimized by calculating an error evaluation value to obtain a target grid that meets a preset error evaluation strategy.

[0100] In this embodiment, as an alternative embodiment, the specific steps include: based on the first arrival spline interpolation function, respectively obtaining the error evaluation values of each node set in the initial node set; traversing each error evaluation value, calculating the difference between the error evaluation value and the previous adjacent error evaluation value, and selecting two node sets with the largest difference and meeting the corresponding error evaluation strategy; calculating the difference set of the two node sets to obtain the target grid; deleting the target grid from the node set containing the target grid, and deleting the node set corresponding to the previous adjacent error evaluation value from the initial node set after deleting the target grid, and performing the step of respectively obtaining the error evaluation values of each node set in the initial node set based on the first arrival spline interpolation function.

[0101] In this embodiment, as an alternative embodiment, the error evaluation strategy includes: the difference is greater than a pre-set difference threshold, and / or the iteration round corresponding to the difference is greater than or equal to a pre-set round threshold, where one iteration round is to perform once the step of respectively obtaining the error evaluation values of each node set in the initial node set based on the first arrival spline interpolation function.

[0102] The fitting function acquisition module 305 is configured to perform fitting based on the first arrival data corresponding to the midpoint of the obtained target grid to obtain the first arrival fitting function.

[0103] In this embodiment, as an alternative embodiment, fitting is performed based on the first arrival data corresponding to the midpoint of the target grid to obtain the first arrival fitting function, which is used for the subsequent calculation of the remaining static correction amount.

[0104] The remaining static correction module 306 is configured to obtain the remaining static correction amount of the first arrival wave based on the first arrival data corresponding to the common shot point, the first arrival data corresponding to the common receiver point, and the first arrival fitting function.

[0105] In this embodiment, as an alternative embodiment, the remaining static correction module 306 includes:

[0106] The common shot point remaining static correction calculation unit is configured to obtain the remaining static correction amount of the shot point in the first arrival data.

[0107] The common receiver point remaining static correction calculation unit is configured to obtain the remaining static correction amount of the receiver point in the first arrival data.

[0108] In this embodiment, as an alternative embodiment, obtaining the remaining static correction amount of the shot point in the first arrival data includes: obtaining the first time of the common shot point in the first arrival fitting function and the second time of the first arrival data of the common shot point, and obtaining the remaining static correction amount of the shot point in the first arrival data based on the first time and the second time.

[0109] Obtaining the residual static correction amount of the geophone points in the first arrival data includes: obtaining the third time of the common geophone points in the first arrival fitting function, and the fourth time of the first arrival data of the common geophone points, and obtaining the residual static correction amount of the geophone points in the first arrival data based on the third time and the fourth time.

[0110] In this embodiment, as an optional embodiment, the first arrival data is corrected based on the residual static correction amount.

[0111] Based on the same inventive concept, an embodiment of the present invention further provides a storage medium, on which a computer program is stored, and when the program is executed by a processor, the steps of the residual static correction method of the first arrival wave in any possible implementation manner described above are implemented.

[0112] Optionally, the storage medium may be a non-temporary computer-readable storage medium. For example, the non-temporary computer-readable storage medium may be ROM, random access memory (RAM), CD-ROM, magnetic tape, floppy disk, and optical data storage device, etc.

[0113] Based on the same inventive concept, see Figure 4 , an embodiment of the present invention further provides an electronic device, including a memory 101 (such as a non-volatile memory), a processor 102, and a computer program stored on the memory 101 and executable on the processor 102. When the processor 102 executes the program, the steps of the residual static correction method of the first arrival wave in any possible implementation manner described above are implemented, which is equivalent to the residual static correction device of the first arrival wave as described above. Of course, the processor can also be used to process other data or operations. The electronic device may be a device such as a PC, a server, or a terminal.

[0114] As Figure 3 shown, the electronic device generally may further include: a memory 103, a network interface 104, and an internal bus 105. In addition to these components, other hardware may also be included, which will not be elaborated here.

[0115] It should be noted that the above-mentioned residual static correction device of the first arrival wave can be implemented by software. As a logically meaningful device, it is formed by the processor 102 of the electronic device where it is located reading the computer program instructions stored in the non-volatile memory into the memory 103 and running.

[0116] The embodiments of the subject matter and the functional operations described in this specification can be implemented in digital electronic circuitry, in tangibly embodied computer software or firmware, in computer hardware including the structures disclosed in this specification and their structural equivalents, or in one or more of them in combination. Embodiments of the subject matter described in this specification can be implemented as one or more computer programs, i.e., one or more modules of computer program instructions encoded on a tangible non-transitory program carrier to be executed by, or to control the operation of, data processing apparatus. Alternatively or additionally, the program instructions can be encoded on an artificially generated propagated signal, e.g., a machine-generated electrical, optical, or electromagnetic signal, that is generated to encode and transmit information to the appropriate receiver apparatus for execution by the data processing apparatus. A computer storage medium may be a machine-readable storage device, a machine-readable storage substrate, a random or serial access memory device, or a combination of one or more of them.

[0117] The processes and logical flows described in this specification can be performed by one or more programmable computers executing one or more computer programs to perform the functions by operating on input data and generating output. The processes and logical flows can also be performed by, or the apparatus can be implemented as, special purpose logic circuitry, e.g., an FPGA (field programmable gate array) or an ASIC (application specific integrated circuit).

[0118] Suitable computers for executing computer programs include, by way of example, both general and / or special purpose microprocessors, or any other type of central processing unit. Generally, a central processing unit will receive instructions and data from a read only memory and / or a random access memory. Basic elements of a computer include a central processing unit for performing or executing instructions and one or more memory devices for storing instructions and data. Generally, a computer will also include one or more mass storage devices for storing data, such as magnetic disks, magneto-optical disks, or optical disks, etc., or the computer will be operatively coupled to such mass storage devices to receive data therefrom or to transfer data thereto, or both. However, a computer need not have such devices. In addition, a computer may be embedded in another device, such as a mobile telephone, a personal digital assistant (PDA), a mobile audio or video player, a game console, a global positioning system (GPS) receiver, or a portable storage device such as a universal serial bus (USB) flash drive, to name just a few.

[0119] Computer-readable media suitable for storing computer program instructions and data include all forms of non-volatile memory, media, and memory devices, such as semiconductor memory devices (e.g., EPROM, EEPROM, and flash memory devices), magnetic disks (e.g., internal hard disks or removable disks), magneto-optical disks, and CD-ROM and DVD-ROM disks. The processor and the memory may be supplemented by, or incorporated in, special purpose logic circuitry.

[0120] Although this specification contains many specific implementation details, these should not be construed as limiting the scope of any invention or the scope of what is claimed, but rather as mainly describing the features of specific embodiments of a particular invention. Certain features described in multiple embodiments in this specification may also be implemented in combination in a single embodiment. On the other hand, the various features described in a single embodiment may also be implemented separately in multiple embodiments or in any suitable sub-combination. Additionally, although features may act in certain combinations as described above and even be initially claimed as such, one or more features from a claimed combination may in some cases be removed from that combination, and the claimed combination may be directed to a sub-combination or a variation of a sub-combination.

[0121] Similarly, although operations are depicted in the drawings in a particular order, this should not be understood as requiring that the operations be performed in the particular order shown or sequentially, or that all illustrated operations be performed, to achieve the desired result. In some cases, multitasking and parallel processing may be advantageous.

[0122] Thus, specific embodiments of the subject matter have been described. Other embodiments are within the scope of the appended claims. In some cases, the acts recited in the claims may be performed in a different order and still achieve the desired result. Additionally, the processes depicted in the drawings are not necessarily in the particular order or sequential order shown to achieve the desired result.

[0123] It should be noted that, in this context, relational terms such as "first" and "second" are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any actual relationship or order between these entities or operations. Moreover, the term "comprising", "including" or any other variation thereof is intended to cover non-exclusive inclusion, such that a process, method, article or device comprising a series of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such process, method, article or device. Without further limitation, an element defined by the statement "comprising one..." does not exclude the presence of additional identical elements in the process, method, article or device comprising the element.

[0124] The above are only specific embodiments of the present invention, enabling those skilled in the art to understand or implement the present invention. Various modifications to these embodiments will be obvious to those skilled in the art, and the general principles defined herein can be implemented in other embodiments without departing from the spirit or scope of the present invention. Therefore, the present invention will not be limited to these embodiments shown herein, but rather to the broadest scope consistent with the principles and novel features claimed herein.

Claims

1. A residual static correction method for first arrival waves, characterized in that, Including: Construct a travel-time graph based on the first arrival data corresponding to each first arrival wave, grid the travel-time graph, and obtain a grid sequence in the gridded travel-time graph; Perform fitting based on the first arrival data to obtain a first arrival spline interpolation function; Traverse the midpoints of each grid in the grid sequence, and construct an initial node set based on the currently traversed grid and the grids that have been traversed; Based on the initial node set and the first arrival spline interpolation function, obtain a target grid that meets a pre-set error evaluation strategy from the initial node set; Perform fitting based on the first arrival data corresponding to the midpoint of the obtained target grid to obtain a first arrival fitting function; Based on the first arrival data corresponding to the common shot point, the first arrival data corresponding to the common receiver point, and the first arrival fitting function, obtain the residual static correction amount of the first arrival wave.

2. The method according to claim 1, wherein The constructing an initial node set based on the currently traversed grid and the grids that have been traversed includes: Extract the midpoint of the first grid in the grid sequence to construct a first node set; Extract the midpoints of the first grid and the second grid in the grid sequence to construct a second node set; Extract the midpoints of the first grid, the second grid until the nth grid in the grid sequence to construct an nth node set, where n is a natural number greater than 2, and the maximum value is the number of grids included in the grid matrix; Based on the first node set, the second node set, and the nth node set, obtain the initial node set.

3. The method according to claim 2, wherein The obtaining a target grid that meets a pre-set error evaluation strategy from the initial node set includes: Based on the first arrival spline interpolation function, respectively obtain the error evaluation values of each node set in the initial node set; Traverse each error evaluation value, calculate the difference between this error evaluation value and the previous adjacent error evaluation value, and select two node sets with the largest difference and meeting the corresponding error evaluation strategy; Calculate the difference set of the two node sets to obtain the target grid; Delete the target grid from the node set containing the target grid, delete the node set corresponding to the previous adjacent error evaluation value from the initial node set after deleting the target grid, and execute the step of respectively obtaining the error evaluation values of each node set in the initial node set based on the first arrival spline interpolation function.

4. The method according to claim 3, wherein The error evaluation strategy includes: the difference is greater than a pre-set difference threshold, and / or the iteration round corresponding to the difference is greater than or equal to a pre-set round threshold, where one iteration round is to execute once the step of respectively obtaining the error evaluation values of each node set in the initial node set based on the first arrival spline interpolation function.

5. The method according to any one of claims 1 to 4, characterized in that, The obtaining the residual static correction amount of the first arrival wave based on the first arrival data corresponding to the common shot point, the first arrival data corresponding to the common receiver point, and the first arrival fitting function includes: Obtain the first time of the common shot point in the first arrival fitting function, and the second time of the first arrival data of the common shot point, and based on the first time and the second time, obtain the residual static correction amount of the shot point in the first arrival data; Obtain the third time of the common geophone point in the first arrival fitting function, and the fourth time of the first arrival data of the common geophone point. Based on the third time and the fourth time, obtain the residual static correction amount of the geophone point in the first arrival data.

6. The method according to any one of claims 1 to 4, characterized in that, The method further includes: Correct the first arrival data based on the residual static correction amount.

7. The method according to any one of claims 1 to 4, characterized in that The gridifying the traveltime graph and obtaining the grid matrix in the gridified traveltime graph includes: Dividing the axis where the offset is located in the traveltime graph according to a preset offset threshold to obtain a plurality of grids perpendicular to the axis where the time is located in the traveltime graph; For each grid, obtain the maximum time value and the minimum time value in the corresponding first arrival data within the grid, calculate the average value of the maximum time value and the minimum time value to obtain the time value of the midpoint of the grid, and calculate the average value of the maximum value and the minimum value of the grid on the axis where the offset is located to obtain the offset of the midpoint of the grid.

8. A residual static correction device for first arrival waves, characterized in that, The residual static correction device for the first arrival wave includes: A grid sequence acquisition module, configured to construct a traveltime graph based on the first arrival data corresponding to each first arrival wave, gridify the traveltime graph, and obtain a grid sequence in the gridified traveltime graph; A spline interpolation function module, configured to perform fitting based on the first arrival data to obtain a first arrival spline interpolation function; An initial node construction module, configured to traverse the midpoints of each grid in the grid sequence, and construct an initial node set based on the currently traversed grid and the already traversed grids; An error evaluation module, configured to obtain a target grid that meets a preset error evaluation strategy from the initial node set based on the initial node set and the first arrival spline interpolation function; A fitting function acquisition module, configured to perform fitting based on the first arrival data corresponding to the midpoint of the obtained target grid to obtain a first arrival fitting function; A residual static correction module, configured to obtain the residual static correction amount of the first arrival wave based on the first arrival data corresponding to the common shot point, the first arrival data corresponding to the common geophone point, and the first arrival fitting function.

9. A storage medium, characterized in that, A program or instruction is stored on a storage medium, and when the program or instruction is run by a processor, the steps of the residual static correction method for the first arrival wave according to any one of claims 1 to 7 are implemented.

10. An electronic device, comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, When the processor executes the program, the steps of the residual static correction method for the first arrival wave according to any one of claims 1 to 7 are implemented.

Images