A method for smoothing high-frequency static correction of surface seismic data

Abstract

The application discloses a high-frequency static correction method for smoothing surface seismic data, which comprises the following steps: performing spatial smoothing processing on an original collected surface to obtain a smoothed surface; calculating a static correction amount of the smoothed surface through a datum surface elevation and a replacement velocity; obtaining a total static correction amount, subtracting the static correction amount of the smoothed surface from the total static correction amount to obtain a residual static correction amount corresponding to the smoothed surface; performing large-scale smoothing processing on the residual static correction amount to obtain a static correction low-frequency component of a low-velocity zone; removing the static correction low-frequency component of the low-velocity zone from the residual static correction amount to obtain a remaining high-frequency static correction amount, wherein the high-frequency static correction amount comprises a high-frequency component of the low-velocity zone and a static correction high-frequency component of a real surface; and applying the high-frequency static correction amount to seismic data and first arrival wave travel time data to complete high-frequency static correction. The method can effectively eliminate high-frequency inter-channel time differences caused by a fluctuating surface in original data, thereby improving the stability and precision of subsequent depth domain velocity modeling.

Description

Technical Field

[0001] The embodiments of the present invention relate to the field of oil and gas geophysical exploration engineering technology, and in particular to a high-frequency static correction method for seismic data of smooth ground surfaces. Background Technology

[0002] As seismic data processing gradually moves towards the depth domain, obtaining high-precision depth-domain processed data for interpreters and geologists to determine the geological development of the work area and ultimately deploy well locations becomes crucial. For depth-domain processing, the most important processing steps are depth-domain velocity modeling and depth-domain migration imaging. The accuracy of velocity affects the refinement of migration imaging and the final seismic data processing. For seismic exploration in piedmont areas, velocity modeling is a key processing technique. Generally, depth-domain velocity modeling mainly relies on depth-domain reflection wave grid tomography. However, with the continuous development of oil and gas exploration in my country, the focus of onshore seismic exploration has gradually shifted to complex and undulating surface areas such as mountains, Gobi Desert, Loess Plateau, and deserts in the west. Therefore, near-surface velocity modeling has become a necessary processing step. Near-surface modeling mainly utilizes reflection waves, calculating theoretical models and picking up first-arrival residuals, performing ray tracing at shot-receiver points to establish tomographic matrices, or calculating gradients using backpropagation residuals, iteratively updating the near-surface velocity model. Finally, the inverted near-surface velocity model was used to continue the grid tomographic velocity modeling of the mid-to-deep reflected waves, and the velocity model of the full depth domain was established.

[0003] Since the accuracy of existing depth domain velocity modeling technology is not yet sufficient to fully characterize the high wavenumber information present in actual underground media, it is urgent to correct the high wavenumber information in actual data and then output it to the subsequent depth domain velocity modeling processing stage. Summary of the Invention

[0004] To address the aforementioned technical problems, at least one embodiment of the present invention provides a high-frequency static correction method for smooth surface seismic data, thereby solving the technical difficulties in the depth domain velocity modeling process of conventional seismic exploration.

[0005] In some optional embodiments, the method includes the following steps:

[0006] The original ground surface is spatially smoothed to obtain a smooth ground surface;

[0007] The static correction amount of the smooth ground surface is calculated using the reference surface elevation and the replacement speed;

[0008] Obtain the total static correction amount, and subtract the static correction amount of the smooth ground surface from the total static correction amount to obtain the residual static correction amount corresponding to the smooth ground surface.

[0009] The residual static correction amount is subjected to large-scale smoothing to obtain the static correction low-frequency component of the low-speed band.

[0010] The low-frequency component of the static correction in the low-speed-deceleration zone is removed from the residual static correction amount to obtain the remaining high-frequency static correction amount, which includes the high-frequency component of the low-speed-deceleration zone and the high-frequency component of the static correction of the actual ground surface.

[0011] The high-frequency static correction is applied to the seismic data and the first arrival wave travel time data to obtain the high-frequency statically corrected seismic data.

[0012] In some optional embodiments, the step of spatially smoothing the original sampled surface to obtain a smooth surface includes:

[0013] Based on the grid information defined in the work area and the actual elevation information of the shot detection points in the work area, the original acquisition surface elevation is obtained by interpolation, and the original acquisition surface is smoothed to obtain a smooth ground surface.

[0014] In some optional embodiments, the interpolation yields the original acquisition surface elevation, including:

[0015] The original elevation of the data acquisition surface is obtained by interpolation using the B-spline interpolation algorithm.

[0016] In some optional embodiments, smoothing the original sampling surface to obtain a smooth surface includes:

[0017] The original acquisition surface is smoothed using a Gaussian filter algorithm of a given radius.

[0018] In some optional embodiments, calculating the static correction amount of the smooth ground surface using the reference surface elevation and replacement speed includes:

[0019] The static correction for a smooth surface is calculated using the following formula, based on the reference elevation and replacement rate:

[0020]

[0021] Among them, Elev datum Elev represents the reference surface elevation used in the time domain processing of the work area. surf Vel indicates a smooth surface. replace This indicates the replacement speed information.

[0022] In some alternative embodiments, the total static correction is obtained by processing the reference surface in the time domain.

[0023] In some optional embodiments, the large-scale smoothing of the residual static correction to obtain the static correction low-frequency component in the low-speed band includes:

[0024] Using an IIR-type Gaussian Filter algorithm, the residual static correction is smoothed on a large scale to obtain the static correction low-frequency component of the low-speed band.

[0025] At least one embodiment of the present invention also provides a high-frequency static correction device for seismic data of smooth ground surfaces, characterized in that it comprises:

[0026] The spatial smoothing module is used to perform spatial smoothing on the original collected ground surface to obtain a smooth ground surface.

[0027] The static correction calculation module is used to calculate the static correction of the smooth ground surface using the reference surface elevation and the replacement speed.

[0028] The residual static correction calculation module is used to obtain the total static correction amount, and subtract the static correction amount of the smooth ground surface from the total static correction amount to obtain the residual static correction amount corresponding to the smooth ground surface.

[0029] The low-frequency component calculation module is used to perform large-scale smoothing on the residual static correction amount to obtain the static correction low-frequency component of the low-speed band.

[0030] A high-frequency component calculation module is used to remove the low-frequency static correction component of the low-speed-deceleration band from the residual static correction amount to obtain the remaining high-frequency static correction amount, which includes the high-frequency component of the low-speed-deceleration band and the static correction high-frequency component of the actual ground surface.

[0031] The high-frequency static correction module is used to apply the high-frequency static correction amount to the seismic data and the first arrival wave travel time data to obtain the high-frequency static corrected seismic data.

[0032] At least one embodiment of the present invention also provides an electronic device, characterized in that it comprises:

[0033] At least one processor; and,

[0034] A memory communicatively connected to the at least one processor; wherein,

[0035] The memory stores instructions that can be executed by the at least one processor to enable the at least one processor to perform the high-frequency static correction method for smooth ground surface seismic data as described above.

[0036] At least one embodiment of the present invention also provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the high-frequency static correction method for smooth ground surface seismic data as described above.

[0037] At least one embodiment of the present invention also provides a computer program product, including a computer program that, when executed by a processor, implements the steps of the high-frequency static correction method for smooth ground surface seismic data as described above.

[0038] Compared with the prior art, the high-frequency static correction method for seismic data on smooth ground surfaces provided by the embodiments of the present invention has the following beneficial effects:

[0039] 1) This method improves the accuracy and stability of depth domain processing and avoids the data mismatch problem caused by using the conventional time domain CMP surface scale to the depth domain as the processing reference surface;

[0040] 2) This method directly smooths the ground surface in the depth domain and uses the smoothed surface to obtain high-frequency static correction values. Then, the seismic data is used to directly correct the data, which fundamentally avoids bringing high-frequency information from the data into the modeling process. The resulting modeling results are relatively accurate and stable, and can reflect accurate underground structural information, thus providing accuracy assurance for subsequent steps in seismic exploration depth domain velocity modeling. Attached Figure Description

[0041] One or more embodiments are illustrated by way of example with reference to the accompanying drawings, and these illustrative descriptions do not constitute a limitation on the embodiments.

[0042] Figure 1 This is a flowchart of the steps of the high-frequency static correction method for seismic data on smooth ground surfaces used in Embodiment 1 of the present invention;

[0043] Figure 2 This is the work area grid information map used in Embodiment 2 of the present invention;

[0044] Figure 3 This is the elevation information map of the shot receiver point used in Embodiment 2 of the present invention;

[0045] Figure 4 This is a true surface elevation map obtained by interpolation in Embodiment 2 of the present invention;

[0046] Figure 5 This is a smooth surface curve obtained by smoothing in Embodiment 2 of the present invention;

[0047] Figure 6 This is a low-frequency static correction diagram calculated according to Embodiment 2 of the present invention;

[0048] Figure 7This is a residual static correction diagram calculated according to Embodiment 2 of the present invention;

[0049] Figure 8 This is a low-frequency static correction diagram of the low-speed reduction band calculated according to Embodiment 2 of the present invention;

[0050] Figure 9 This is the high-frequency static correction quantity diagram calculated according to Embodiment 2 of the present invention;

[0051] Figure 10 This is a comparison chart of seismic data before and after the high-frequency static correction output in Embodiment 2 of the present invention. Detailed Implementation

[0052] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the various embodiments of the present invention will be described in detail below with reference to the accompanying drawings. However, those skilled in the art will understand that many technical details are presented in the embodiments of the present invention to facilitate a better understanding of the invention. However, the technical solutions claimed in the present invention can be implemented even without these technical details and various variations and modifications based on the following embodiments. The division of the following embodiments is for ease of description and should not constitute any limitation on the specific implementation of the present invention. The various embodiments can be combined with and referenced by each other without contradiction.

[0053] As mentioned earlier, existing depth-domain velocity modeling techniques cannot fully capture the high-wavenumber information present in actual subsurface media. Therefore, it is necessary to correct the high-wavenumber information in the actual data before outputting it to the subsequent depth-domain velocity modeling processing stage. Theoretically, depth-domain processing should start from the original acquisition surface. The currently recognized depth-domain processing surface is the "simulated surface," which is a slightly smoothed result of the original data acquisition surface. The "simulated surface" is a relatively smooth physical surface that is closer to the real surface after removing high-frequency inter-channel time differences that cannot be distinguished by PSDM and depth-domain velocity analysis, and retaining the low-frequency static correction caused by near-surface velocity anomalies. It is a physically meaningful processing surface that meets the needs of depth-domain velocity analysis and imaging. Through the slightly smoothed surface, the corresponding seismic data and the picked first arrivals should also be corrected to this surface, avoiding the introduction of high-frequency information from the data into the modeling processing stage. The resulting modeling results are relatively accurate and stable, reflecting accurate subsurface structural information and yielding high-precision imaging results.

[0054] The implementation details of the above concept are explained in detail below through examples. The following content is only for the convenience of understanding and is not necessary for implementing this solution.

[0055] Example 1:

[0056] like Figure 1 As shown in the figure, this embodiment provides a high-frequency static correction method for seismic data of smooth ground surfaces. The method mainly includes the following steps:

[0057] S1, Real Surface Spatial Smoothing: Spatial smoothing is performed on the original collected ground surface. An appropriate smoothing radius is selected for smoothing to obtain a "simulated ground surface", that is, a smooth ground surface.

[0058] S2, Calculation of low-frequency components of smooth ground surface: Given a smooth ground surface, the static correction amount of the smooth ground surface can be calculated based on the elevation of the reference plane and the replacement speed.

[0059] S3, Calculation of residual static correction: Obtain the total static correction from the time domain processing reference surface of the work area, and subtract the static correction of the smooth ground surface from the total static correction to obtain the residual static correction corresponding to the smooth ground surface; in essence, the residual static correction includes the low-frequency component of static correction of the low-speed zone, the high-frequency component of static correction of the low-speed zone, and the high-frequency component of static correction of the true ground surface.

[0060] S4, Extraction of static correction low-frequency components in low speed reduction zone: For the residual static correction amount, perform large-scale smoothing on the residual static correction amount (the large scale referred to here is a scale known to those in the art, usually above 2 kilometers, and is not limited here), thus obtaining the static correction low-frequency components in low speed reduction zone.

[0061] S5, High-frequency static correction component extraction: The above residual static correction includes the low-frequency component of the low-speed reduction zone, the high-frequency component of the low-speed reduction zone, and the high-frequency component caused by the real ground surface. Therefore, by removing the low-frequency component of the low-speed reduction zone from the residual static correction, the high-frequency static correction component can be obtained, which is the high-frequency component of the low-speed reduction zone and the high-frequency component caused by the real ground surface.

[0062] S6 applies the high-frequency static correction component to the seismic data and first arrival travel time data to complete the high-frequency static correction and obtain the high-frequency static corrected seismic data for subsequent depth domain velocity modeling.

[0063] In this embodiment, step S1, true surface spatial smoothing, involves statistically analyzing the lead-in information of the original acquired seismic data to obtain the shot point elevation, receiver elevation, and shot point depth, thus acquiring scattered surface elevation information for the entire work area. Based on the processing grid defined for the work area, the true ground surface, i.e., the original acquisition surface, is interpolated, and then this original acquisition surface is smoothed using a Gaussian smoothing method with a given radius. The specific formula is as follows:

[0064] F(x,y)=scatterInterp(g(x,y),x t ,y t ,z t )

[0065] Where F(x,y) represents the interpolated true ground surface, scatterInterp represents the interpolation algorithm (which can be arbitrarily chosen; this embodiment uses B-spline interpolation), and g(x,y) represents the position information of the x and y coordinates in the defined work area grid. t y t , z t This indicates the actual location coordinates and elevation information of the shot receiver.

[0066] F s (x,y)=smooth(F(x,y),r)

[0067] Among them, F s (x, y) represents a small smooth surface obtained by smoothing with radius r. Here, "smooth" represents the smoothing algorithm using the IIR-type Gaussian Filter algorithm. This algorithm is based on a recursive structure, with only 6 multiplication-addition operations per dimension, and is independent of the Gaussian kernel size. This avoids the drawback of traditional algorithms where computational efficiency decreases as the smoothing radius increases. Assuming the input signal is in[n], the intermediate variables, output signals w[n], and out[n] can be calculated using the following formulas:

[0068] w[n] = B·in[n]

[0069] +(b1·w[n-1]+b2·w[n-2]+b3·w[n-3])

[0070] / b0

[0071] out[n] = B·w[n]

[0072] +(b1·out[n+1]+out·w[n+2]+b3

[0073] ·w[n+3]) / b0

[0074] in:

[0075]

[0076] in:

[0077]

[0078] Here, sigma is the smoothing parameter. The larger the sigma, the greater the smoothness.

[0079] In this embodiment, the static correction amount of the smooth surface in step S2 is calculated as follows: Given a smooth surface, the static correction amount can be calculated based on the elevation of the reference surface and the replacement speed. The formula is as follows:

[0080]

[0081] Among them, Elev datum Elev represents the reference surface elevation used in the time domain processing of the work area. surf Vel represents the small, smooth surface obtained in step S1. replace This indicates the replacement speed information.

[0082] In this embodiment, the residual static correction amount in step S3 above is calculated by subtracting the static correction amount for the smooth ground surface from the original static correction amount, using the following formula:

[0083] Static_residual=(Static_datum-Static_low_surf) (6)

[0084] Wherein, Static_residual represents the residual static correction, which essentially includes the low-frequency component of the static correction in the low-velocity zone, the high-frequency component of the static correction in the low-velocity zone, and the high-frequency component of the static correction on the true ground surface. Static_datum represents the total static correction obtained from time-domain processing, and Static_low_surf represents the static correction for the smooth ground surface calculated by S2.

[0085] In this embodiment, the static correction low-frequency component extraction of the low-speed band in step S4 is as follows: the residual static correction amount obtained in S3 is smoothed on a large scale, and the static correction low-frequency component caused by the low-speed band is obtained using the IIR-type Gaussian Filter to which S1 belongs, as shown in the following formula:

[0086] Static_low_weather=Smooth(Static_residual,r)

[0087] Among them, Static_low_weather represents the low-frequency component of the static correction in the low-speed band.

[0088] In this embodiment, the high-frequency static correction quantity extraction in step S5 is as follows: The residual static correction quantity in S3 includes the low-frequency component of the low-speed reduction zone, the high-frequency component of the low-speed reduction zone, and the high-frequency component caused by the real ground surface. Therefore, the high-frequency static correction quantity can be obtained by removing the low-frequency component of the low-speed reduction zone from the residual static correction. The high-frequency static correction quantity includes the high-frequency component of the low-speed reduction zone and the high-frequency component caused by the real ground surface.

[0089] Static_high=Static_residual-Static_low_weather

[0090] Wherein, Static_high represents the final extracted high-frequency static correction amount.

[0091] In this embodiment, step S6 above applies the processed high-frequency static correction amount to the seismic data and first arrival wave travel time data to obtain high-frequency statically corrected seismic data, forming a high-frequency static correction method and technical process for smooth surface seismic data. The final output of this invention is high-frequency static correction information for shot points and receiver points, which can be extended to other work areas with severe surface undulations for application, helping to obtain seismic data more suitable for depth domain velocity modeling.

[0092] Example 2

[0093] The technical solution of the present invention and its beneficial effects will be further illustrated below with a specific example.

[0094] As shown in Figure 1, this embodiment uses a high-frequency static correction method for seismic data on smooth ground surfaces, which includes the following steps:

[0095] S1, true surface space smoothing, based on the grid information defined for the work area and the actual elevation information of the shot receiver points in the work area, such as... Figure 2 and Figure 3 As shown, interpolation yields the true elevation, i.e., the original surface elevation, and smoothing results in a small, smoothed surface, as shown. Figure 4 and Figure 5 As shown.

[0096] S2, Calculation of static correction for smooth ground surface: Based on the smooth ground surface information obtained in S1, as well as the work area's reference surface and replacement speed information, the static correction for the smooth ground surface is calculated as follows: Figure 6 As shown.

[0097] S3, Calculation of residual static correction: Based on the total static correction information of the work area and the static correction of the smooth ground surface obtained in S2, the corresponding residual static correction of the work area is calculated as follows: Figure 7 As shown.

[0098] S4, Extraction of low-frequency components of static correction in the low-speed-deceleration band: Based on the residual static correction obtained in S3, perform large-scale smoothing to obtain the static correction amount in the low-speed-deceleration band, as shown below. Figure 8 As shown.

[0099] S5, High-frequency static correction extraction: Using the residual static correction obtained in S3 and the low-frequency quantity in the low-speed band obtained in S4, the remaining high-frequency static correction is obtained as follows: Figure 9 As shown.

[0100] S6, apply the processed high-frequency static correction to the seismic data and first arrival wave travel time data to obtain the high-frequency statically corrected seismic data as follows: Figure 10 As shown.

[0101] from Figure 8 In terms of application results, the high-frequency static correction method for smooth surface seismic data provided by this invention can effectively eliminate the high inter-channel time difference caused by undulating surface in the original data, thereby improving the stability and accuracy of subsequent depth domain velocity modeling.

[0102] Example 3

[0103] Another embodiment of the present invention relates to a high-frequency static correction device for seismic data of a smooth ground surface, comprising:

[0104] The spatial smoothing module is used to perform spatial smoothing on the original collected ground surface to obtain a smooth ground surface.

[0105] The static correction calculation module is used to calculate the static correction of the smooth ground surface using the reference surface elevation and the replacement speed.

[0106] The residual static correction calculation module is used to obtain the total static correction amount, and subtract the static correction amount of the smooth ground surface from the total static correction amount to obtain the residual static correction amount corresponding to the smooth ground surface.

[0107] The low-frequency component calculation module is used to perform large-scale smoothing on the residual static correction amount to obtain the static correction low-frequency component of the low-speed band.

[0108] A high-frequency component calculation module is used to remove the low-frequency static correction component of the low-speed-deceleration band from the residual static correction amount to obtain the remaining high-frequency static correction amount, which includes the high-frequency component of the low-speed-deceleration band and the static correction high-frequency component of the actual ground surface.

[0109] The high-frequency static correction module is used to apply the high-frequency static correction amount to the seismic data and the first arrival wave travel time data to complete the high-frequency static correction and obtain the high-frequency static corrected seismic data.

[0110] Example 4:

[0111] Another embodiment of the present invention relates to an electronic device, comprising: at least one processor; and a memory communicatively connected to the at least one processor; wherein the memory stores instructions executable by the at least one processor, the instructions being executed by the at least one processor to enable the at least one processor to perform the high-frequency static correction method for smooth ground surface seismic data in the above embodiments.

[0112] The memory and processor are connected via a bus, which can include any number of interconnecting buses and bridges, connecting various circuits of one or more processors and memories. The bus can also connect various other circuits, such as peripheral devices, voltage regulators, and power management circuits, which are well known in the art and will not be described further herein. The bus interface provides an interface between the bus and the transceiver. The transceiver can be a single element or multiple elements, such as multiple receivers and transmitters, providing a unit for communicating with various other devices over a transmission medium. Data processed by the processor is transmitted over the wireless medium via an antenna, which further receives data and transmits it to the processor.

[0113] The processor manages the bus and general processing, and also provides various functions, including timing, peripheral interfaces, voltage regulation, power management, and other control functions. Memory is used to store data used by the processor during operation.

[0114] Example 5:

[0115] Another embodiment of the present invention relates to a computer-readable storage medium storing a computer program. When executed by a processor, the computer program implements the high-frequency static correction method for smooth ground surface seismic data described in the above embodiments.

[0116] That is, those skilled in the art will understand that all or part of the steps in the methods of the above embodiments can be implemented by a program instructing related hardware. This program is stored in a storage medium and includes several instructions to cause a device (which may be a microcontroller, chip, etc.) or processor to execute all or part of the steps of the methods described in the various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0117] Example 6

[0118] Another embodiment of the present invention relates to a computer program product, including a computer program that, when executed by a processor, implements the steps of the high-frequency static correction method for smooth ground surface seismic data described above.

[0119] Those skilled in the art will understand that the above embodiments are specific embodiments for implementing the present invention, and in practical applications, various changes in form and detail may be made without departing from the spirit and scope of the present invention.

Claims

1. A method for high-frequency static correction of seismic data on smooth ground surfaces, characterized in that, include: The original ground surface is spatially smoothed to obtain a smooth ground surface; The static correction amount of the smooth ground surface is calculated using the reference surface elevation and the replacement speed; Obtain the total static correction amount, and subtract the static correction amount of the smooth ground surface from the total static correction amount to obtain the residual static correction amount corresponding to the smooth ground surface. The residual static correction amount is subjected to large-scale smoothing to obtain the static correction low-frequency component of the low-speed band. The low-frequency component of the static correction in the low-speed-deceleration zone is removed from the residual static correction amount to obtain the remaining high-frequency static correction amount, which includes the high-frequency component of the low-speed-deceleration zone and the high-frequency component of the static correction of the actual ground surface. The high-frequency static correction is applied to the seismic data and the first arrival wave travel time data to complete the high-frequency static correction.

2. The high-frequency static correction method for seismic data on smooth ground surfaces according to claim 1, characterized in that, The process of spatially smoothing the original surface of the collected data to obtain a smooth surface includes: Based on the grid information defined in the work area and the actual elevation information of the shot detection points in the work area, the original acquisition surface elevation is obtained by interpolation, and the original acquisition surface is smoothed to obtain a smooth ground surface.

3. The high-frequency static correction method for seismic data on smooth ground surfaces according to claim 2, characterized in that, The interpolation yields the original acquisition surface elevation, including: The original elevation of the data acquisition surface is obtained by interpolation using the B-spline interpolation algorithm.

4. The high-frequency static correction method for seismic data on smooth ground surfaces according to claim 2, characterized in that, The process of smoothing the original sampling surface to obtain a smooth ground surface includes: The original acquisition surface is smoothed using a Gaussian filter algorithm of a given radius.

5. The high-frequency static correction method for seismic data on smooth ground surfaces according to claim 1, characterized in that, The calculation of the static correction amount of the smooth ground surface using the reference surface elevation and replacement speed includes: The static correction for a smooth surface is calculated using the following formula, based on the reference elevation and replacement rate: Among them, Elev datum Elev represents the reference surface elevation used in the time domain processing of the work area. surf Vel represents a smooth surface. replace This indicates the replacement speed information.

6. The high-frequency static correction method for seismic data on smooth ground surfaces according to claim 1, characterized in that, The total static correction is obtained by processing the reference surface in the time domain.

7. The high-frequency static correction method for seismic data on smooth ground surfaces according to claim 1, characterized in that, The large-scale smoothing process performed on the residual static correction amount to obtain the static correction low-frequency component in the low-speed band includes: Using an IIR-type Gaussian Filter algorithm, the residual static correction is smoothed on a large scale to obtain the static correction low-frequency component of the low-speed band.

8. A high-frequency static correction device for seismic data on a smooth ground surface, characterized in that, include: The spatial smoothing module is used to perform spatial smoothing on the original collected ground surface to obtain a smooth ground surface. The static correction calculation module is used to calculate the static correction of the smooth ground surface using the reference surface elevation and the replacement speed. The residual static correction calculation module is used to obtain the total static correction amount, and subtract the static correction amount of the smooth ground surface from the total static correction amount to obtain the residual static correction amount corresponding to the smooth ground surface. The low-frequency component calculation module is used to perform large-scale smoothing on the residual static correction amount to obtain the static correction low-frequency component of the low-speed band. A high-frequency component calculation module is used to remove the low-frequency static correction component of the low-speed-deceleration band from the residual static correction amount to obtain the remaining high-frequency static correction amount, which includes the high-frequency component of the low-speed-deceleration band and the static correction high-frequency component of the actual ground surface. The high-frequency static correction module is used to apply the high-frequency static correction amount to the seismic data and the first arrival wave travel time data to complete the high-frequency static correction.

9. An electronic device, characterized in that, include: At least one processor; as well as, A memory communicatively connected to the at least one processor; wherein, The memory stores instructions that can be executed by the at least one processor to enable the at least one processor to perform the high-frequency static correction method for seismic data of smooth ground surfaces as described in any one of claims 1 to 7.

10. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed by the processor, it implements the high-frequency static correction method for seismic data of smooth ground surfaces as described in any one of claims 1 to 7.

Images