Surface consistency residual static correction method, electronic equipment, storage medium and device

Abstract

The invention discloses a surface consistency residual static correction method, electronic equipment, a storage medium and a device. The method comprises the following steps: acquiring first arrival data of seismic rotation waves; generating an XYO gather based on the first arrival data; establishing a revolution wave earth surface consistency convolution model based on the XYO gather to obtain a seismic trace and a model trace; correlation operation of the seismic trace and the model trace is completed in the XYO trace set, and the time shift amount between the seismic trace and the model trace is obtained; and calculating a shot point residual static correction value and a detection point residual static correction value of the XYO gather based on the time shift amount. According to the method, the characteristics of strong revolution wave energy, high signal-to-noise ratio, no false information interference and sensitivity to near-surface change of the vibroseis are utilized to establish the revolution wave surface consistency convolution model, so that residual static correction calculation of shot points and detection points is completed; according to the method, the problem of systematic static time shifting caused by the near-surface effect in the extremely thick desert area can be effectively solved.

Description

Technical Field

[0001] The present invention belongs to the technical field of ground controllable vibroseis seismic data processing, and more specifically, relates to a surface consistency residual static correction method, electronic equipment, storage medium and device. Background Art

[0002] There are many thick desert exploration areas in western my country (such as the Tarim Basin and the Tengger Desert). Since the surface conditions in these work areas are typical aeolian desert landforms, the thickness of the sand dunes is large. Due to the high efficiency, safety and environmental protection of the vibroseis technology, the vibroseis technology has been applied to such thick desert areas by major domestic oilfield companies. However, the sand dunes are thick and the surface is loose, so the surface absorption attenuation is serious, the energy transmission is limited, and the signal of the exploration target layer cannot be seen in the collected data. Therefore, it is impossible to use the target layer to complete the work required for the remaining static correction (such as the correlation between the model channel and each channel in the CDP channel gather).

[0003] Compared with the extremely weak signal of the target layer, the rotation wave of the controllable source has strong energy, high signal-to-noise ratio, no false information interference (such as no coherent noise and refracted multiple waves), and is sensitive to near-surface changes. Therefore, these characteristics of the rotation wave can be used to complete the residual static correction calculation of the shot point and the detection point.

[0004] The information disclosed in the background technology section of the present invention is only intended to deepen the understanding of the general background technology of the present invention, and should not be regarded as acknowledging or suggesting in any form that the information constitutes the prior art already known to those skilled in the art. Summary of the invention

[0005] The purpose of the present invention is to propose a surface consistency residual static correction method, electronic equipment, storage medium and device, which can utilize the characteristics of high signal-to-noise ratio of rotating waves and sensitivity to near-surface changes to effectively solve the residual static correction of shot points and detection points caused by near-surface effects in thick desert areas.

[0006] To achieve the above-mentioned purpose, the present invention proposes a surface consistency residual static correction method, electronic equipment, storage medium and device

[0007] According to a first aspect of the present invention, a surface consistency residual static correction method is proposed, comprising:

[0008] Obtain the first arrival data of seismic rotation waves;

[0009] Generate an XYO gather based on the first arrival data;

[0010] Among them, X represents the longitudinal direction, Y represents the transverse direction, and O represents the offset;

[0011] Based on the XYO gathers, a surface consistent convolution model of the rotating wave is established to obtain seismic traces and model traces;

[0012] Complete correlation calculations between the seismic trace and the model trace in the XYO trace gather to obtain a time shift between the seismic trace and the model trace;

[0013] The shot point residual static correction amount and the detection point residual static correction amount of the XYO gather are calculated based on the time shift amount.

[0014] Optionally, the obtaining of the first arrival data of the seismic rotation wave specifically includes:

[0015] Applying the field long wavelength static correction and refraction static correction to the shot gather data, flattening the first arrival wave using the existing rotation wave velocity, and intercepting the data of the set time interval up and down along the first arrival wave to obtain the first arrival data.

[0016] Optionally, generating an XYO gather based on the first arrival data comprises:

[0017] The first arrival data are sorted into the XYO domain to generate the XYO gather.

[0018] Optionally, the expression for calculating the residual static correction amount of the shot point and the receiver point of the XYO gather based on the time shift is:

[0019]

[0020] in, is the time shift between the seismic trace and the model trace, is the shot point time shift, i.e. the residual static correction of the shot point, is the detection point time shift, that is, the residual static correction of the detection point, τ q It is the time shift of XYO gather.

[0021] Optionally, the conjugate gradient method or Gauss-Seidel method is used to solve the expression The residual static correction amount of the shot point and the residual static correction amount of the detection point are obtained.

[0022] According to a second aspect of the present invention, a surface consistency residual static correction device is proposed, comprising:

[0023] An acquisition module is used to acquire the first arrival data of the seismic rotation wave after completing the field long wavelength static correction and refraction static correction;

[0024] A generation module, used for generating XYO gathers based on the first arrival data;

[0025] Among them, X represents the longitudinal direction, Y represents the transverse direction, and O represents the offset;

[0026] An establishment module is used to establish a surface consistent convolution model of the rotary wave based on the XYO gather to obtain seismic traces and model traces;

[0027] A correlation operation module, used for completing the correlation operation between the seismic trace and the model trace in the XYO trace gather to obtain the time shift between the seismic trace and the model trace;

[0028] A calculation module is used to calculate the residual static correction amount of the shot point and the residual static correction amount of the detection point of the XYO gather based on the time shift amount.

[0029] Optionally, the expression for calculating the residual static correction amount of the shot point and the receiver point of the XYO gather based on the time shift is:

[0030]

[0031] in, is the time shift between the seismic trace and the model trace, is the time shift of the gun point, is the detection point time shift, τ q It is the time shift of XYO gather.

[0032] Optionally, the conjugate gradient method or Gauss-Seidel method is used to solve the expression The residual static correction amount of the shot point and the residual static correction amount of the detection point are obtained.

[0033] According to a third aspect of the present invention, an electronic device is provided, the electronic device comprising:

[0034] at least one processor; and,

[0035] a memory communicatively connected to the at least one processor; wherein,

[0036] The memory stores instructions that can be executed by the at least one processor, and the instructions are executed by the at least one processor so that the at least one processor can execute the surface consistency residual static correction method described in any one of the first aspects.

[0037] According to a fourth aspect of the present invention, a non-transitory computer-readable storage medium is proposed, which stores computer instructions for causing a computer to execute any surface consistency residual static correction method described in the first aspect.

[0038] The beneficial effects of the present invention are as follows: the present invention utilizes the characteristics of strong energy of the rotation wave of the controllable seismic source, high signal-to-noise ratio, no false information interference, and sensitivity to near-surface changes to establish a rotation wave surface consistency convolution model, and then completes the residual static correction calculation of the shot point and the detection point. The present invention can effectively solve the systematic static time shift problem caused by the near-surface effect in the thick desert area.

[0039] The system of the present invention has other characteristics and advantages, which will be apparent from the drawings incorporated herein and the following detailed description, or will be described in detail in the drawings incorporated herein and the following detailed description, which together serve to explain the specific principles of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0040] The above and other objects, features and advantages of the present invention will become more apparent through a more detailed description of exemplary embodiments of the present invention in conjunction with the accompanying drawings, in which like reference numerals generally represent like components.

[0041] Figure 1 A flow chart showing the steps of a surface consistency residual static correction method according to the present invention.

[0042] Figure 2 A schematic diagram of a three-dimensional common turning point gather CTP according to the present invention is shown.

[0043] Figure 3 A schematic diagram of the common turning point gather CTP ray path according to the present invention is shown.

[0044] Figure 4 A schematic diagram of a surface consistency residual static correction device according to embodiment 2 of the present invention is shown. DETAILED DESCRIPTION

[0045] The present invention will be described in more detail below with reference to the accompanying drawings. Although preferred embodiments of the present invention are shown in the accompanying drawings, it should be understood that the present invention can be implemented in various forms and should not be limited by the embodiments set forth herein. On the contrary, these embodiments are provided to make the present invention more thorough and complete, and to fully convey the scope of the present invention to those skilled in the art.

[0046] like Figure 1 As shown, a surface consistency residual static correction method according to the present invention includes:

[0047] Obtain the first arrival data of seismic rotation waves;

[0048] Generate XYO gathers based on first arrival data;

[0049] Among them, X represents the longitudinal direction, Y represents the transverse direction, and O represents the offset;

[0050] Based on XYO gathers, a surface consistent convolution model of rotating waves is established to obtain seismic traces and model traces.

[0051] Complete the correlation operation between the seismic trace and the model trace in the XYO gather to obtain the time shift between the seismic trace and the model trace;

[0052] The residual statics of the shot points and receiver points of the XYO gather are calculated based on the time shift.

[0053] Specifically, the present invention utilizes the characteristics of the controllable seismic source, such as strong energy of the rotation wave, high signal-to-noise ratio, no false information interference, and sensitivity to near-surface changes, to complete the residual static correction calculation of the shot point and the detection point. First, the first arrival data of the seismic rotation wave is obtained, and then an XYO gather is generated based on the first arrival data. In the western desert area, the velocity and thickness of the low-deceleration zone of the near-surface structure vary greatly, but the underground structure is similar within a certain range. Therefore, we approximately assume that there is a common turning point (common turnpoint), and extract the first arrival data into an XYO gather according to the size of the rotation face element and the size of the offset face element (X represents longitudinal, Y represents lateral, and O represents offset). The schematic diagram of the three-dimensional common turning point gather CTP is shown as follows Figure 2 As shown in the figure, dx is the longitudinal interval, dy is the lateral interval, and df is the offset interval. In an XYO gather, different offsets represent different layers of the underground structure, so an XYO gather can represent the underground structure at its turning point, such as Figure 3 As shown in the figure, based on this characteristic, the propagation effect (Green's function) caused by the background velocity can be obtained in the XYO gather, which can be called the background principal value effect (Green's function). For each trace of the XYO gather, the background principal value effect is removed, and the remaining part is the surface inconsistency factor caused by the near-surface effect. Due to the inhomogeneity of the near-surface, the seismic reflection signal is different at different shot points and receiver points. In order to correct the inhomogeneity of the near-surface, the surface consistency decomposition is usually completed using the target layer reflection signal:

[0054] P ij (t) = S i (t)*R i (t)*Γ ij (t) (1)

[0055] In the above formula, P ij (t) represents the seismic trace received by the jth receiver of the i-th shot point, which is decomposed into three convolution effect components: the convolution effect S of the i-th shot point i (t), the convolution effect R of the jth detection point j(t) and the propagation effect (Green’s function) Γ of the ith shot point exciting the jth receiver point ij (t); i, j represent the unique index numbers of the shot point and the detection point respectively, i=1,…,N s ,j=1,…,N R ; N S ,N R The total number of unique index numbers of the shot points and the detection points respectively; transform formula (1) to the frequency domain:

[0056]

[0057] In the above formula, They are P ij (t),S i (t),R i (t),Γ ij (t), ω is the angular frequency. Similar to the surface consistent convolution model of seismic reflection signals, the surface consistent convolution model can also be established using the backscattered wave, thereby separating the near-surface effects at the shot point and the receiver point. In order to compensate for the propagation effect (Green's function), the backscattered wave surface consistent convolution model needs to be calculated in the XYO gather. In the XYO gather, the backscattered wave surface consistent convolution model can be expressed as:

[0058] Time domain: P b (t)=W(t)*S i (t)*R j (t)*Q(t) (3)

[0059] Frequency domain:

[0060] In the above formula, S i (t),R j (t) Surface consistency effect of shot point and receiver point, YesS i (t),R j (t); W(t) is the average response of the XYO channel gather (also called the model channel), which can be approximated by the superposition channel W′(t) of the XYO channel gather, and its Fourier transform is Q(t) is the residual term related to the underground structure, is the Fourier transform of Q(t), ω is the angular frequency; b is the index number that uniquely corresponds to the index pair of the shot point and the receiver point, (i, j) → b, i = 1, ..., N s ,j=1,…,N R , N S ,N RThe total number of unique index numbers of the shot points and the receiver points respectively; that is, the surface consistent convolution model of the rotating wave is established based on the XYO data set to obtain the seismic trace and the model trace; the correlation operation of the seismic trace and the model trace is completed in the XYO data set to obtain the time shift between the seismic trace and the model trace; the expression (4) is rearranged to obtain the expression on the right side that only contains the near-surface effect:

[0061]

[0062] The expression (5) can be further separated into the expressions of amplitude and phase:

[0063]

[0064]

[0065] Where θ is the phase at each frequency. The residual static correction caused by surface inconsistency is expressed as a phase error in the frequency domain. Therefore, when estimating the surface consistency phase error, it is necessary to assume that the phase error is a constant. Therefore, the phase error in the frequency domain (expression (7)) can be approximated by the time shift constant, that is, the second formula in expression (5) can be written as:

[0066]

[0067] In the formula is the seismic trace P′ b The time shift between (t) and the model track W′(t), and τ q Respectively represent the time shift of the shot point, receiver point and XYO gather, where They are also called the residual corrections of the shot point and the detector point. Through P' b The cross-correlation between W(t) and W′(t) is obtained, and then the residual static correction of the shot point and the residual static correction of the receiver point of the XYO data set are calculated based on the time shift. Expression (8) defines a set of equations for all shot-receiver combinations (i, j). The conjugate gradient method or Gauss-Seidel method can be used to solve expression (8); the calculated residual static correction and They are the subsequent supplements to the field long-wavelength static correction and refraction static correction respectively; the present invention is suitable for the processing of controllable source data in extremely thick desert areas, and makes full use of the characteristics of high signal-to-noise ratio of the return wave and sensitivity to near-surface changes, and effectively solves the systematic static time shift (residual static correction amount) caused by near-surface effects in extremely thick desert areas.

[0068] In one example, obtaining the first arrival data of the seismic rotation wave specifically includes:

[0069] Apply the field long wavelength static correction and refraction static correction to the shot gather data, use the existing rotation wave velocity to flatten the first arrival wave, and intercept the data of the set time interval up and down along the first arrival wave to obtain the first arrival data.

[0070] Specifically, the present invention applies field long-wavelength static correction and refraction static correction to the shot gather data, uses the existing rotation wave velocity to flatten the first arrival wave, and intercepts data of set time intervals up and down along the first arrival wave, for example, intercepts 200-400ms data up and down along the first arrival wave, to obtain the first arrival data.

[0071] In one example, generating an XYO gather based on first arrival data includes:

[0072] Sort the first arrival data into the XYO domain and generate XYO gathers.

[0073] In one example, the expression for calculating the residual static correction of the shot points and receiver points of the XYO gather based on the time shift is:

[0074]

[0075] in, is the time shift between the seismic trace and the model trace, is the shot point time shift, i.e. the residual static correction of the shot point, is the detection point time shift, that is, the residual static correction of the detection point, τ q It is the time shift of XYO gather.

[0076] Specifically, Separate into amplitude and phase expressions:

[0077]

[0078]

[0079] Where θ is the phase at each frequency. The residual static correction caused by surface inconsistency is expressed as a phase error in the frequency domain. Therefore, when estimating the surface consistency phase error, it is necessary to assume that the phase error is a constant. Therefore, the phase error in the frequency domain (expression (7)) can be approximated by the time shift constant, that is, the second formula in expression (5) can be written as:

[0080]

[0081] In the formula is the seismic trace P′ b The time shift between (t) and the model track W′(t), and τ q Respectively represent the time shift of the shot point, receiver point and XYO gather, where They are also called the residual corrections of the shot point and the detector point. Through P' b The cross-correlation between W(t) and W′(t) is obtained.

[0082] In one example, the conjugate gradient method or Gauss-Seidel method is used to solve the expression The residual static correction of the shot point and the residual static correction of the detection point are obtained.

[0083] Specifically, the expression For all the gun inspection combinations (i, j), a system of equations is defined, which can be solved by the conjugate gradient method or the Gauss-Seidel method; the residual static correction calculated is and They are the subsequent supplements to the field long-wavelength static correction and refraction static correction respectively.

[0084] The present invention is further described below in conjunction with the accompanying drawings and specific embodiments, but is not intended to be limiting of the present invention. It should be noted that the embodiments and features in the embodiments of the present invention may be combined with each other without conflict.

[0085] Example 1

[0086] This embodiment provides a surface consistency residual static correction method, including:

[0087] Obtain the first arrival data of the seismic rotation wave, apply the field long wavelength static correction and refraction static correction to the shot gather data, use the existing rotation wave velocity to flatten the first arrival wave, and intercept the data of the set time interval up and down along the first arrival wave to obtain the first arrival data;

[0088] Generate XYO gathers based on first arrival data, sort the first arrival data into XYO domain, and generate XYO gathers; where X represents the longitudinal direction, Y represents the transverse direction, and O represents the offset;

[0089] Based on XYO gathers, a surface consistent convolution model of rotating waves is established to obtain seismic traces and model traces.

[0090] Complete the correlation operation between the seismic trace and the model trace in the XYO gather to obtain the time shift between the seismic trace and the model trace;

[0091] The residual static correction of the shot points and the residual static correction of the receiver points of the XYO gather are calculated based on the time shift. The expression for calculating the residual static correction of the shot points and the receiver points of the XYO gather is:

[0092]

[0093] in, is the time shift between the seismic trace and the model trace, is the shot point time shift, i.e. the residual static correction of the shot point, is the detection point time shift, that is, the residual static correction of the detection point, τ q It is the time shift of XYO gather.

[0094] Example 2

[0095] like Figure 4 As shown, this embodiment provides a surface consistency residual static correction device, including:

[0096] An acquisition module is used to acquire the first arrival data of the seismic rotation wave after completing the field long wavelength static correction and refraction static correction;

[0097] The generation module is used to generate XYO gathers based on the first arrival data;

[0098] Among them, X represents the longitudinal direction, Y represents the transverse direction, and O represents the offset;

[0099] Establish a module for establishing a surface consistent convolution model of the rotary wave based on the XYO gathers to obtain seismic traces and model traces;

[0100] The correlation operation module is used to complete the correlation operation between the seismic trace and the model trace in the XYO trace gather to obtain the time shift between the seismic trace and the model trace;

[0101] The calculation module is used to calculate the residual static correction of the shot point and the residual static correction of the receiver point of the XYO gather based on the time shift.

[0102] Example 3

[0103] This embodiment provides an electronic device, the electronic device comprising:

[0104] at least one processor; and,

[0105] a memory communicatively connected to the at least one processor; wherein,

[0106] The memory stores instructions that can be executed by the at least one processor, and the instructions are executed by the at least one processor so that the at least one processor can perform the surface consistency residual static correction method in Example 1.

[0107] The electronic device according to an embodiment of the present disclosure includes a memory and a processor, and the memory is used to store non-temporary computer-readable instructions. Specifically, the memory may include one or more computer program products, and the computer program product may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may, for example, include a random access memory (RAM) and / or a cache memory (cache), etc. The non-volatile memory may, for example, include a read-only memory (ROM), a hard disk, a flash memory, etc.

[0108] The processor may be a central processing unit (CPU) or other forms of processing units having data processing capabilities and / or instruction execution capabilities, and may control other components in the electronic device to perform desired functions. In one embodiment of the present disclosure, the processor is used to run the computer-readable instructions stored in the memory.

[0109] Those skilled in the art should be able to understand that in order to solve the technical problem of how to obtain a good user experience, the present embodiment may also include well-known structures such as a communication bus and an interface, and these well-known structures should also be included in the protection scope of the present disclosure.

[0110] For detailed description of this embodiment, reference may be made to the corresponding descriptions in the aforementioned embodiments, which will not be repeated here.

[0111] Example 4

[0112] This embodiment provides a non-transitory computer-readable storage medium, which stores computer instructions, and the computer instructions are used to enable a computer to execute the surface consistency residual static correction method in Example 1.

[0113] According to the computer-readable storage medium of the embodiment of the present disclosure, non-transitory computer-readable instructions are stored thereon. When the non-transitory computer-readable instructions are executed by a processor, all or part of the steps of the method of each embodiment of the present disclosure are executed.

[0114] The above-mentioned computer-readable storage media include, but are not limited to: optical storage media (e.g., CD-ROM and DVD), magneto-optical storage media (e.g., MO), magnetic storage media (e.g., magnetic tape or mobile hard disk), media with built-in rewritable non-volatile memory (e.g., memory card) and media with built-in ROM (e.g., ROM box).

[0115] Example 5

[0116] This embodiment provides a surface consistency residual static correction method, including:

[0117] In the western desert area, the velocity and thickness of the low-velocity zone of the near-surface structure vary greatly, but the underground structure is similar within a certain range. Therefore, we approximately assume that there is a common turnpoint, and extract the first arrival data into an XYO gather according to the size of the turnaround bin and the offset bin (X represents the longitudinal direction, Y represents the lateral direction, and O represents the offset). In an XYO gather, different offsets represent different strata of the underground structure, so an XYO gather band can represent the underground structure at its turning point. Based on this characteristic, the propagation effect (Green's function) caused by the background velocity can be obtained in the XYO gather, which can be called the background principal value effect (Green's function). For each track of the XYO gather, the background principal value effect is removed, and the remaining part is the surface inconsistency factor caused by the near-surface effect. Due to the heterogeneity of the near-surface, the seismic reflection signal is different at different shot points and detection points. In order to correct the heterogeneity of the near-surface, the reflection signal of the target layer is usually used to complete the surface consistency decomposition:

[0118] P ij (t) = S i (t)*R i (t)*Γ ij (t) (1)

[0119] In the above formula, P ij (t) represents the seismic trace received by the jth receiver of the i-th shot point, which is decomposed into three convolution effect components: the convolution effect S of the i-th shot point i (t), the convolution effect R of the jth detection point j (t) and the propagation effect (Green’s function) Γ of the ith shot point exciting the jth receiver point ij (t); i, j represent the unique index numbers of the shot point and the detection point respectively, i=1,…,N s ,j=1,…,N R ; N S ,N R The total number of unique index numbers of the shot points and the detection points respectively; transform formula (1) to the frequency domain:

[0120]

[0121] In the above formula, They are P ij (t),S i (t),R i (t),Γ ij(t), ω is the angular frequency. Similar to the surface consistent convolution model of seismic reflection signals, the surface consistent convolution model can also be established using the backscattered wave, thereby separating the near-surface effects at the shot point and the receiver point. In order to compensate for the propagation effect (Green's function), the backscattered wave surface consistent convolution model needs to be calculated in the XYO gather. In the XYO gather, the backscattered wave surface consistent convolution model can be expressed as:

[0122] Time domain: P b (t)=W(t)*S i (t)*R j (t)*Q(t) (3)

[0123] Frequency domain:

[0124] In the above formula, S i (t),R j (t) Surface consistency effect of shot point and receiver point, YesS i (t),R j (t); W(t) is the average response of the XYO channel gather (also called the model channel), which can be approximated by the superposition channel W′(t) of the XYO channel gather, and its Fourier transform is Q(t) is the residual term related to the underground structure, is the Fourier transform of Q(t), ω is the angular frequency; b is the index number that uniquely corresponds to the index pair of the shot point and the receiver point, (i, j) → b, i = 1, ..., N s ,j=1,…,N R , N S ,N R The total number of unique index numbers of the shot points and the receiver points respectively; rearrange expression (4) to obtain the expression on the right side that only contains the near-surface effect:

[0125]

[0126] The expression (5) can be further separated into the expressions of amplitude and phase:

[0127]

[0128]

[0129] Where θ is the phase at each frequency. The residual static correction caused by surface inconsistency is expressed as a phase error in the frequency domain. Therefore, when estimating the surface consistency phase error, it is necessary to assume that the phase error is a constant. Therefore, the phase error in the frequency domain (expression (7)) can be approximated by the time shift constant, that is, the second formula in expression (5) can be written as:

[0130]

[0131] In the formula is the seismic trace P′ b The time shift between (t) and the model track W′(t), and τ q Respectively represent the time shift of the shot point, receiver point and XYO gather, where They are also called the residual corrections of the shot point and the detector point. Through P' b The cross-correlation between W(t) and W′(t) is obtained, and then the residual static correction of the shot point and the residual static correction of the receiver point of the XYO data set are calculated based on the time shift. Expression (8) defines a set of equations for all shot-receiver combinations (i, j). The conjugate gradient method or Gauss-Seidel method can be used to solve expression (8); the calculated residual static correction and They are the subsequent supplements to the field long-wavelength static correction and refraction static correction respectively.

[0132] The specific process of this embodiment is as follows:

[0133] Apply field long wavelength static correction and refraction static correction to the shot gather data;

[0134] The first arrival wave is flattened using the existing rotation wave velocity, and 200-400ms data are intercepted along the first arrival wave;

[0135] Sorting the intercepted data into the XYO domain

[0136] Complete seismic trace P′ in all XYO gathers b (t) is correlated with the superposition channel W′(t) to obtain the correlation channel. The time corresponding to the maximum amplitude on the correlation channel is the time shift to be obtained.

[0137] The conjugate gradient method or Gauss-Seidel method is used to solve equation (8) to obtain the residual static correction of all shot points and receiver points.

[0138] The embodiments of the present invention have been described above, and the above description is exemplary, not exhaustive, and is not limited to the disclosed embodiments. Many modifications and changes will be apparent to those skilled in the art without departing from the scope and spirit of the described embodiments.

Claims

1. A surface consistency residual static correction method, characterized in that: include: Obtain the first arrival data of seismic rotation waves; Generate an XYO gather based on the first arrival data; Among them, X represents the longitudinal direction, Y represents the transverse direction, and O represents the offset; Based on the XYO gathers, a surface consistent convolution model of the rotating wave is established to obtain seismic traces and model traces; Complete correlation calculations between the seismic trace and the model trace in the XYO trace gather to obtain a time shift between the seismic trace and the model trace; The shot point residual static correction amount and the detection point residual static correction amount of the XYO gather are calculated based on the time shift amount.

2. The surface consistency residual static correction method according to claim 1, characterized in that: The obtaining of the first arrival data of the seismic rotation wave specifically includes: Applying the field long wavelength static correction and refraction static correction to the shot gather data, flattening the first arrival wave using the existing rotation wave velocity, and intercepting the data of the set time interval up and down along the first arrival wave to obtain the first arrival data.

3. The surface consistency residual static correction method according to claim 1, characterized in that: Generating an XYO gather based on the first arrival data comprises: The first arrival data are sorted into the XYO domain to generate the XYO gather.

4. The surface consistency residual static correction method according to claim 1, characterized in that: The expression for calculating the residual static correction of the shot point and the receiver point of the XYO gather based on the time shift is: in, is the time shift between the seismic trace and the model trace, is the shot point time shift, i.e. the residual static correction of the shot point, is the detection point time shift, that is, the residual static correction of the detection point, τ q It is the time shift of XYO gather.

5. The surface consistency residual static correction method according to claim 4, characterized in that: Use the conjugate gradient method or Gauss-Seidel method to solve the expression The residual static correction amount of the shot point and the residual static correction amount of the detection point are obtained.

6. An electronic device, characterized in that: The electronic device comprises: 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, and the instructions are executed by the at least one processor so that the at least one processor can execute the surface consistency residual static correction method described in any one of claims 1-5.

7. A non-transitory computer-readable storage medium, characterized in that: The non-transitory computer-readable storage medium stores computer instructions, which are used to enable a computer to execute the surface consistency residual static correction method described in any one of claims 1-5.

8. A surface consistency residual static correction device, characterized in that: include: An acquisition module is used to acquire the first arrival data of the seismic rotation wave after completing the field long wavelength static correction and refraction static correction; A generation module, used for generating XYO gathers based on the first arrival data; Among them, X represents the longitudinal direction, Y represents the transverse direction, and O represents the offset; An establishment module is used to establish a surface consistent convolution model of the rotary wave based on the XYO gather to obtain seismic traces and model traces; A correlation operation module, used for completing the correlation operation between the seismic trace and the model trace in the XYO trace gather to obtain the time shift between the seismic trace and the model trace; A calculation module is used to calculate the residual static correction amount of the shot point and the residual static correction amount of the detection point of the XYO gather based on the time shift amount.

9. The surface consistency residual static correction device according to claim 8, characterized in that: The expression for calculating the residual static correction of the shot point and the receiver point of the XYO gather based on the time shift is: in, is the time shift between the seismic trace and the model trace, is the time shift of the gun point, is the detection point time shift, τ q It is the time shift of XYO gather.

10. The surface consistency residual static correction device according to claim 9, characterized in that: Use the conjugate gradient method or Gauss-Seidel method to solve the expression The residual static correction amount of the shot point and the residual static correction amount of the detection point are obtained.

Images