A method, device and storage medium for improving the resolution of seismic data

Abstract

The present invention discloses a method, apparatus and storage medium for improving the resolution of seismic data. The method includes: obtaining the spatio-temporal variable wavelet amplitude spectrum of the current seismic trace; constructing a frequency-variable compression coefficient corresponding to the current seismic trace according to the spatio-temporal variable wavelet amplitude spectrum; calculating the spatio-temporal variable wide-frequency wavelet amplitude spectrum of the current seismic trace according to the frequency-variable compression coefficient; calculating a spatio-temporal variable frequency-up operator corresponding to the current seismic trace based on the spatio-temporal variable wavelet amplitude spectrum and the spatio-temporal variable wide-frequency wavelet amplitude spectrum; and performing frequency-up processing on the current seismic trace by using the spatio-temporal variable frequency-up operator. The frequency-variable wavelet based on time and space extracted by the solution of the present invention is more suitable for the needs of complex structure work areas; the spatio-temporal variable frequency-up operator obtained based on the spatio-temporal variable wavelet amplitude spectrum and the spatio-temporal variable wide-frequency wavelet amplitude spectrum can adaptively determine the corresponding spatio-temporal variable frequency-up operator according to the local spectral characteristics of the seismic data. Even when applied to the seismic data of complex structure work areas, it can effectively improve its resolution.

Description

Technical Field

[0001] The present invention relates to the technical field of oil and gas exploration and data processing, and particularly relates to a method, device, and storage medium for improving the resolution of seismic data. Background Art

[0002] For geophysical exploration, high-resolution seismic data helps to distinguish thinner strata and more precisely describe the characteristics of small fault blocks, small faults, and reservoir edges of oil and gas reservoirs. Moreover, how to improve the resolution of seismic data has always been a difficult point in geophysical exploration.

[0003] In the prior art, a method for improving the resolution of seismic data by designing the compression coefficient of a frequency-variable wavelet is proposed. By performing frequency-up conversion processing on the post-stack seismic data, a good effect of improving the resolution of seismic data can be achieved. However, in the prior art solution, a single wavelet is extracted for the entire work area. According to the spectral characteristics of seismic data, when the subsurface structure of the work area is very complex, the spectral characteristics of seismic data also vary greatly in the time and space domain. That is to say, the method of only extracting a single wavelet in the prior art cannot effectively describe the seismic spectral characteristics of a complex structure work area, resulting in the vertical resolution of the seismic data after frequency-up conversion processing being effectively improved in some areas, while not being significantly improved in other areas. More seriously, the method of only extracting a single wavelet in the prior art also causes a large change in the spatio-temporal relative relationship of the amplitude and energy of seismic data, seriously affecting the amplitude preservation of the seismic data after improving the resolution. Summary of the Invention

[0004] The technical problem to be solved by the present invention is to provide a solution for improving the resolution of seismic data applicable to a complex structure work area, so that not only can the resolution of seismic data be improved, but also the seismic data after improving the resolution has good amplitude preservation.

[0005] To solve the above technical problem, the present invention provides a method for improving the resolution of seismic data, including:

[0006] Obtain the spatio-temporal variable wavelet amplitude spectrum of the current seismic trace;

[0007] Construct the frequency-variable compression coefficient corresponding to the current seismic trace according to the spatio-temporal variable wavelet amplitude spectrum;

[0008] Calculate the spatio-temporal variable wide-frequency wavelet amplitude spectrum of the current seismic trace according to the frequency-variable compression coefficient;

[0009] Calculate the spatio-temporal variable frequency-up conversion operator corresponding to the current seismic trace based on the spatio-temporal variable wavelet amplitude spectrum and the spatio-temporal variable wide-frequency wavelet amplitude spectrum;

[0010] Use the spatio-temporal frequency-up operator to perform frequency-up processing on the current seismic trace.

[0011] Optionally, the obtaining of the spatio-temporal wavelet amplitude spectrum of the current seismic trace includes:

[0012] Determine each surrounding seismic trace corresponding to the current seismic trace from the seismic data of the three-dimensional target work area;

[0013] Obtain the spatio-temporal wavelet amplitude spectrum of the current seismic trace according to the average time-frequency amplitude spectra of each surrounding seismic trace.

[0014] Optionally, the determining of each surrounding seismic trace corresponding to the current seismic trace from the seismic data of the three-dimensional target work area includes:

[0015] Use a spatial sliding window to determine each surrounding seismic trace corresponding to the current seismic trace from the seismic data of the three-dimensional target work area according to the set search radius value.

[0016] Optionally, the obtaining of the spatio-temporal wavelet amplitude spectrum of the current seismic trace according to the average time-frequency amplitude spectra of each surrounding seismic trace includes:

[0017] Use the short-time Fourier transform to calculate the time-frequency spectra of each surrounding seismic trace, and calculate the average time-frequency amplitude spectra of each surrounding seismic trace based on the time-frequency spectra;

[0018] Obtain the spatio-temporal wavelet amplitude spectrum of the current seismic trace according to the average time-frequency amplitude spectra.

[0019] Optionally, the obtaining of the spatio-temporal wavelet amplitude spectrum of the current seismic trace according to the average time-frequency amplitude spectra includes:

[0020] Perform time-frequency domain logarithmic spectrum simulation on the average time-frequency amplitude spectra to obtain the spatio-temporal wavelet amplitude spectrum of the current seismic trace.

[0021] Optionally, the constructing of the frequency-variable compression coefficient corresponding to the current seismic trace according to the spatio-temporal wavelet amplitude spectrum includes:

[0022] Determine the frequency-variable compression coefficient corresponding to the current seismic trace based on the peak frequency in the spatio-temporal wavelet amplitude spectrum.

[0023] Optionally, the determining of the frequency-variable compression coefficient corresponding to the current seismic trace based on the peak frequency in the spatio-temporal wavelet amplitude spectrum includes:

[0024] Determine the frequency-variable compression coefficient a(t,f) corresponding to the current seismic trace according to the following expression,

[0025]

[0026] where f is the frequency of the current seismic trace, f p (t) is the peak frequency of the spatio-temporal variable wavelet amplitude spectrum at the t-th sampling point of the current seismic trace, f Nyq is the Nyquist frequency, f g is a frequency constant with a value range between 1 and 2, a max is the maximum compression coefficient of the frequency-variable wavelet.

[0027] Optionally, the frequency up-conversion processing of the current seismic trace using the spatio-temporal variable frequency up-conversion operator includes:

[0028] Performing time-frequency domain frequency up-conversion processing on the current seismic trace according to the spatio-temporal variable frequency up-conversion operator, and obtaining the time-frequency spectrum after the frequency up-conversion processing;

[0029] Summing the time-frequency spectrum in the time direction, then performing inverse Fourier transform on the summation result, and taking the real part of the inverse Fourier transform as the result of the frequency up-conversion processing.

[0030] To solve the above technical problems, the present invention provides a device for improving the resolution of seismic data, including:

[0031] A spatio-temporal variable wavelet amplitude spectrum acquisition module, configured to acquire the spatio-temporal variable wavelet amplitude spectrum of the current seismic trace;

[0032] A compression coefficient construction module, configured to construct the frequency-variable compression coefficient corresponding to the current seismic trace according to the spatio-temporal variable wavelet amplitude spectrum;

[0033] A wide-band wavelet amplitude spectrum acquisition module, configured to calculate the spatio-temporal variable wide-band wavelet amplitude spectrum of the current seismic trace according to the frequency-variable compression coefficient;

[0034] A frequency up-conversion operator calculation module, configured to calculate the spatio-temporal variable frequency up-conversion operator corresponding to the current seismic trace based on the spatio-temporal variable wavelet amplitude spectrum and the spatio-temporal variable wide-band wavelet amplitude spectrum;

[0035] A frequency up-conversion processing module, configured to perform frequency up-conversion processing on the current seismic trace using the spatio-temporal variable frequency up-conversion operator.

[0036] Optionally, the spatio-temporal variable wavelet amplitude spectrum acquisition module includes a surrounding seismic trace determination unit and a spatio-temporal variable wavelet amplitude spectrum acquisition unit; wherein,

[0037] The surrounding seismic trace determination unit is configured to determine each surrounding seismic trace corresponding to the current seismic trace from the seismic data of the three-dimensional target work area;

[0038] The spatio-temporal variable wavelet amplitude spectrum acquisition unit is configured to acquire the spatio-temporal variable wavelet amplitude spectrum of the current seismic trace according to the average time-frequency amplitude spectrum of each surrounding seismic trace.

[0039] Optionally, the peripheral seismic trace determination unit is specifically configured to determine each peripheral seismic trace corresponding to the current seismic trace from the seismic data of the three-dimensional target work area by using a spatial sliding window according to a set search radius value.

[0040] Optionally, the spatio-temporal variable wavelet amplitude spectrum acquisition unit includes an average time-frequency amplitude spectrum calculation subunit and a spatio-temporal variable wavelet amplitude spectrum acquisition subunit; wherein,

[0041] The average time-frequency amplitude spectrum calculation subunit calculates the time-frequency spectrum of each peripheral seismic trace by using short-time Fourier transform, and calculates the average time-frequency amplitude spectrum of each peripheral seismic trace based on the time-frequency spectrum;

[0042] The spatio-temporal variable wavelet amplitude spectrum acquisition subunit is configured to obtain the spatio-temporal variable wavelet amplitude spectrum of the current seismic trace according to the average time-frequency amplitude spectrum.

[0043] Optionally, the spatio-temporal variable wavelet amplitude spectrum acquisition subunit is specifically configured to perform time-frequency domain logarithmic spectrum simulation on the average time-frequency amplitude spectrum to obtain the spatio-temporal variable wavelet amplitude spectrum of the current seismic trace.

[0044] Optionally, the compression coefficient construction module is configured to determine the frequency-variable compression coefficient corresponding to the current seismic trace based on the peak frequency in the spatio-temporal variable wavelet amplitude spectrum.

[0045] Optionally, the compression coefficient construction module is specifically configured to determine the frequency-variable compression coefficient a(t,f) corresponding to the current seismic trace according to the following expression,

[0046]

[0047] where f is the frequency of the current seismic trace, f p (t) is the peak frequency of the spatio-temporal variable wavelet amplitude spectrum of the t-th sampling point of the current seismic trace, f Nyq is the Nyquist frequency, f g is a frequency constant with a value range of 1 to 2, a max is the maximum compression coefficient of the frequency-variable wavelet.

[0048] Optionally, the frequency up-conversion processing module is specifically configured to perform time-frequency domain frequency up-conversion processing on the current seismic trace according to the spatio-temporal variable frequency up-conversion operator, and obtain the time-frequency spectrum after frequency up-conversion processing; sum the time-frequency spectrum in the time direction, then perform inverse Fourier transform on the summation result, and use the real part of the inverse Fourier transform as the result of the frequency up-conversion processing.

[0049] To solve the above technical problems, the present invention provides a computer device, including a memory, a processor, and a computer program stored on the memory and executable on the processor. When the processor executes the computer program, the above method is implemented.

[0050] In order to solve the above technical problem, the present invention provides a computer-readable storage medium, on which a computer program is stored, and the program implements the above method when executed by a processor.

[0051] Compared with the prior art, one or more embodiments of the above scheme may have the following advantages or beneficial effects:

[0052] When the resolution of seismic data is improved by applying the scheme of the present invention, the space-time variable wavelet amplitude spectrum of the current seismic trace is first obtained; then the frequency-variant compression coefficient corresponding to the current seismic trace is constructed based on the space-time variable wavelet amplitude spectrum, and then the space-time variable wideband frequency wavelet amplitude spectrum of the current seismic trace is calculated based on the frequency-variant compression coefficient; based on the space-time variable wavelet amplitude spectrum and the space-time variable wideband frequency wavelet amplitude spectrum, the space-time variable frequency boosting operator corresponding to the current seismic trace is calculated; and then the space-time variable frequency boosting operator is used to perform frequency boosting processing on the current seismic trace.

[0053] It can be seen that the frequency-variant wavelet based on time and space extracted by the scheme of the present invention is more suitable for the needs of complex structural work areas; the time-space variable frequency-raising operator obtained based on the time-space variable wavelet amplitude spectrum and the time-space variable broadband frequency-raising operator can adaptively determine the corresponding time-space variable frequency-raising operator according to the local frequency spectrum characteristics of the seismic data, so that even the seismic data applied to complex structural work areas can effectively improve its resolution; in addition, according to experiments, it can be seen that this scheme does not destroy the relative relationship of the time-space energy of the collected seismic data, and therefore has a high amplitude preservation property. BRIEF DESCRIPTION OF THE DRAWINGS

[0054] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0055] Figure 1a is the original seismic data;

[0056] Figure 1b Based on a single wavelet pair Figure 1a Seismic data after frequency enhancement processing of the original earthquake;

[0057] Figure 1c To apply the solution of the present invention to Figure 1a Seismic data after frequency enhancement processing of the original earthquake;

[0058] Figure 1d for Figure 1a A partial enlarged view of the original seismic data;

[0059] Figure 1e The local enlarged view of the seismic data after frequency up-conversion processing of the original seismic data in Figure 1a ;

[0060] Figure 1f The amplitude spectrum superposition display chart of the original seismic data and the seismic data after frequency up-conversion processing using the solution of the present invention;

[0061] Figure 2 A flowchart of the method for improving the resolution of seismic data provided by an embodiment of the present invention;

[0062] Figure 3 Another flowchart of the method for improving the resolution of seismic data provided by an embodiment of the present invention;

[0063] Figure 4a The average time-frequency amplitude spectrum corresponding to the xline number in the original seismic data;

[0064] Figure 4b The average time-frequency amplitude spectrum corresponding to the xline number in the seismic data after frequency up-conversion processing of the original seismic data in Figure 4a using the solution of the present invention;

[0065] Figure 5 Another flowchart of the method for improving the resolution of seismic data provided by an embodiment of the present invention;

[0066] Figure 6 Another flowchart of the method for improving the resolution of seismic data provided by an embodiment of the present invention;

[0067] Figure 7 A structural diagram of the device for improving the resolution of seismic data provided by an embodiment of the present invention;

[0068] Figure 8 The internal structural diagram of the frequency up-conversion operator calculation module provided by an embodiment of the present invention;

[0069] Figure 9 A structural diagram of the computer device provided by an embodiment of the present invention. Detailed implementation manners

[0070] Next, the technical solutions in the embodiments of the present invention will be clearly and completely described in conjunction with the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, rather than all of the embodiments. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

[0071] In the above-mentioned prior art, the method of improving the resolution of seismic data by designing the compression coefficient of frequency-variable wavelets extracts a single wavelet for the entire work area. When the subsurface structure of the work area is very complex, the spectral characteristics of seismic data also vary greatly in the time-space domain. That is to say, the solution provided in the prior art cannot effectively describe the seismic spectral characteristics of complex structure work areas, resulting in the effective improvement of the vertical resolution of seismic data in some areas after frequency-up conversion processing, while there is no obvious improvement in some other areas. More seriously, the method of only extracting a single wavelet in the prior art will also cause a large change in the spatio-temporal relative relationship of the amplitude and energy of seismic data (please refer to Figure 1a and Figure 1b ), seriously affecting the amplitude preservation of seismic data after improving the resolution.

[0072] Therefore, in order to improve the resolution of seismic data in complex structure work areas, the present invention provides the following method, device and storage medium for improving the resolution of seismic data.

[0073] First, the method for improving the resolution of seismic data provided by the present invention will be described below.

[0074] Example 1

[0075] As Figure 2 shown, it is a flowchart of a method for improving the resolution of seismic data provided by an embodiment of the present invention, which may include the following steps:

[0076] Step S101: Obtain the spatio-temporal variable wavelet amplitude spectrum of the current seismic trace.

[0077] Step S102: Construct the frequency-variable compression coefficient corresponding to the current seismic trace according to the spatio-temporal variable wavelet amplitude spectrum.

[0078] Step S103: Calculate the spatio-temporal variable wide-frequency wavelet amplitude spectrum of the current seismic trace according to the frequency-variable compression coefficient.

[0079] Step S104: Calculate the spatio-temporal variable frequency-up conversion operator corresponding to the current seismic trace based on the spatio-temporal variable wavelet amplitude spectrum and the spatio-temporal variable wide-frequency wavelet amplitude spectrum.

[0080] Step S105: Perform frequency-up conversion processing on the current seismic trace by using the spatio-temporal variable frequency-up conversion operator.

[0081] It can be seen that the frequency-variable wavelets based on time and space extracted by the solution of the present invention are more suitable for the needs of complex structure work areas; the spatio-temporal variable frequency-up conversion operator obtained based on the spatio-temporal variable wavelet amplitude spectrum and the spatio-temporal variable wide-frequency wavelet amplitude spectrum can adaptively determine the corresponding spatio-temporal variable frequency-up conversion operator according to the local spectral characteristics of seismic data. In this way, even when applied to seismic data in complex structure work areas, the resolution can be effectively improved.

[0082] In addition, please refer to Figure 1c , according to experiments, it can be known that the spatio-temporal energy relative relationship of the collected seismic data is not damaged by this solution, so it has high amplitude preservation. For the method for improving the resolution of seismic data provided by the embodiments of the present invention, due to the adoption of the spatio-temporal variable wavelet processing technology, not only the resolution of the processed data is improved, but also the spatial relative relationship of the amplitude and energy is basically consistent with the characteristics of the original data, and its amplitude preservation is greatly improved. In order to further demonstrate the frequency-upconversion effect of the embodiments of the present invention, Figure 1d and Figure 1e give a partial enlarged comparison diagram of the original seismic data and the seismic data after frequency-upconversion processing. It can be seen that the resolution of the frequency-upconverted data is effectively improved, the number of isophases is significantly increased, and the signal-to-noise ratio and lateral continuity of the profile are both high.

[0083] Figure 1f gives the amplitude spectrum superposition display of the original seismic data and the seismic data after frequency-upconversion processing. The black line is the amplitude spectrum of the original seismic data, and the gray line is the amplitude spectrum of the seismic data after frequency-upconversion processing. It can be seen that the high-frequency region of the amplitude spectrum of the seismic data after frequency-upconversion processing is effectively improved relative to the amplitude spectrum of the original seismic data, while the energy at the low-frequency end is consistent with the original data, verifying that the solution of the present invention can effectively protect the low-frequency information of the original seismic data.

[0084] Example 2

[0085] As Figure 3 shown, it is another flowchart of the method for improving the resolution of seismic data provided by the embodiments of the present invention, which may include the following steps:

[0086] Step S201: Determine each surrounding seismic trace corresponding to the current seismic trace from the seismic data of the three-dimensional target work area.

[0087] In one implementation, each surrounding seismic trace corresponding to the current seismic trace can be determined from the seismic data of the three-dimensional target work area by using a spatial sliding window according to a set search radius value.

[0088] Step S202: Obtain the spatio-temporal variable wavelet amplitude spectrum of the current seismic trace according to the average time-frequency amplitude spectra of each surrounding seismic trace.

[0089] Preferably, to ensure that the phase relationship of the original seismic data is not changed after frequency-upconversion processing, the spatio-temporal variable wavelet extracted by the embodiments of the present invention is a zero-phase wavelet.

[0090] Furthermore, since the spatio-temporal variable wavelet is extracted from the average time-frequency amplitude spectra of each seismic trace within the range covered by a spatial sliding window centered on the current seismic trace (such as a local rectangular window), and this spatial sliding window changes according to the position change of the current seismic trace, the extracted spatio-temporal variable wavelet can have both time-varying and space-varying characteristics. It should also be noted that the spatio-temporal variable wavelet extracted based on the spatial sliding window changes more smoothly, and compared with the conventional single-trace frequency extraction algorithm, it can further improve the lateral continuity of the frequency extraction process.

[0091] In one implementation, the short-time Fourier transform is used to calculate the time-frequency spectra of each surrounding seismic trace, and the average time-frequency amplitude spectrum of each surrounding seismic trace is calculated based on the time-frequency spectra; the spatio-temporal variable wavelet amplitude spectrum of the current seismic trace is obtained according to the average time-frequency amplitude spectrum.

[0092] Specifically, the average time-frequency amplitude spectrum can be subjected to time-frequency domain logarithmic spectrum simulation to obtain the spatio-temporal variable wavelet amplitude spectrum of the current seismic trace.

[0093] The following combines specific examples to illustrate the process of obtaining the spatio-temporal variable wavelet amplitude spectrum.

[0094] (1) Set the search radius α of the spatial sliding window, and extract each seismic trace within the range covered by the spatial sliding window according to the main survey line number (i.e., inline number) and connecting line number (i.e., xline number) of the current seismic trace. For example, if the inline number corresponding to the current seismic trace is i and the xline number is j, then the seismic traces within this spatial sliding window include: each seismic trace with the inline number ranging from (i - α) to (i + α) and the xline number ranging from (j - α) to (j + α).

[0095] (2) Use the short-time Fourier transform, such as the Gabor transform, to calculate the time-frequency spectra of each seismic trace within the spatial sliding window, and calculate the average time-frequency amplitude spectrum of the time-frequency spectra of each seismic trace. For example, calculate the average time-frequency amplitude spectrum corresponding to the current seismic trace with the inline number i and the xline number j according to the following expression

[0096]

[0097] where t is time, f is frequency, and (2α + 1) 2 is the number of seismic traces within the spatial sliding window; sp i+k,j+m (t, f) is the time-frequency spectrum of the seismic trace with the inline number i + k and the xline number j + m, and this time-frequency spectrum can be obtained by the Gabor transform, where the ranges of both k and m are (-α, α).

[0098] (3) Based on the average time-frequency amplitude spectrum, the time-space variable wavelet amplitude spectrum of the current seismic trace is obtained by using the time-frequency domain logarithmic spectrum simulation technique.

[0099] It should be noted that in order to ensure that the time-frequency amplitude spectrum of the extracted time-space variable wavelet is smooth, the logarithmic domain average time-frequency amplitude spectrum can be polynomially fitted. The specific implementation process is as follows: for each time sampling point t, the least squares method is used to optimize the objective function to determine the polynomial coefficients β n , n = 0, 1, 2... N, where N is the pre-given polynomial order. The objective function is as follows,

[0100]

[0101] Correspondingly, the time-space variable wavelet amplitude spectrum WAsp i,j (t, f) corresponding to the t-th sampling point of the current seismic trace with inline number i and xline number j is:

[0102]

[0103] It should be noted that the aforementioned Gabor transform is only a specific algorithm of the short-time Fourier transform. In the embodiments of the present invention, there is no need to limit the specific algorithm of the short-time Fourier transform adopted, and those skilled in the art can make reasonable settings according to the specific circumstances in actual applications.

[0104] Furthermore, according to Figure 4a and Figure 4b it can be known that based on the average time-frequency amplitude spectrum, by using the time-frequency domain logarithmic spectrum simulation technique, the time-frequency amplitude spectrum of the time-space variable wavelet of the current seismic trace can be obtained, as Figure 4b shown. It can be seen that the time-frequency amplitude spectrum characteristics of the time-space variable wavelet of the current seismic trace are similar to the average time-frequency amplitude characteristics, but smoother, so it can meet the needs of time-space variable frequency-up processing.

[0105] Step S203: Construct the frequency-variable compression coefficient corresponding to the current seismic trace according to the time-space variable wavelet amplitude spectrum.

[0106] It should be noted that the role of the constructed frequency-variable compression coefficient factor is to compress the time-space variable wavelet to obtain the time-space variable wide-frequency wavelet.

[0107] Step S204: Calculate the time-space variable wide-frequency wavelet amplitude spectrum of the current seismic trace according to the frequency-variable compression coefficient.

[0108] Step S205: Calculate the time-space variable frequency-up operator corresponding to the current seismic trace based on the time-space variable wavelet amplitude spectrum and the time-space variable wide-frequency wavelet amplitude spectrum.

[0109] Step S206: Use the spatio-temporal frequency-up conversion operator to perform frequency-up conversion processing on the current seismic trace.

[0110] Figure 3 The method embodiments shown have Figure 2 all the beneficial effects of the method embodiments shown. In addition, the spatio-temporal wavelet in the embodiments of the present invention is obtained by extracting based on a spatial sliding window, and its change is more gentle. Compared with the conventional single-trace frequency-up conversion algorithm, it can further improve the lateral continuity of the frequency-up conversion processing. In addition, since the conventional frequency-variable compression coefficient method for improving resolution in the prior art is based on a single wavelet, the corresponding frequency-variable compression coefficient and spatio-temporal variable broadband wavelet are also fixed. However, in the embodiments of the present invention, by combining the time-frequency amplitude spectrum characteristics of the spatio-temporal wavelet, an adaptive low-frequency-preserving spatio-temporal frequency-up conversion operator is proposed on the basis of the conventional frequency-variable compression coefficient, which can effectively improve the resolution of seismic data in complex work areas.

[0111] Example 3

[0112] As Figure 5 shown, it is another flowchart of the method for improving the resolution of seismic data provided by the embodiments of the present invention, which may include the following steps:

[0113] Step S301: Obtain the spatio-temporal wavelet amplitude spectrum of the current seismic trace.

[0114] Step S302: Based on the peak frequency in the spatio-temporal wavelet amplitude spectrum, determine the frequency-variable compression coefficient corresponding to the current seismic trace.

[0115] In one implementation, the frequency-variable compression coefficient a(t, f) corresponding to the current seismic trace can be determined according to the following expression:

[0116]

[0117] where f is the frequency of the current seismic trace, f p (t) is the peak frequency of the spatio-temporal wavelet amplitude spectrum of the t-th sampling point of the current seismic trace, f Nyq is the Nyquist frequency, f g is a frequency constant with a value range between 1 and 2, and a max is the maximum compression coefficient of the frequency-variable wavelet. It should be noted that a max is the maximum compression coefficient of the frequency-variable wavelet, and according to empirical values, it is usually sufficient to take 2.

[0118] It can be seen from the above expression that when the frequency value is between 0 Hz and the peak frequency f pWhen the frequency value is between (t), the determined frequency-variable compression coefficient is 1, indicating that no operation is performed on the frequency axis to ensure the consistency between the low-frequency components of the broadband wavelet and the spatio-temporal variable wavelet. When the frequency value is between the peak frequency f p (t) and the Nyquist frequency f Nyq , the compression factor is a max , indicating that the frequency axis needs to be stretched by a factor of a max to ensure that the broadband wavelet has rich high-frequency information. When the frequency value is between the peak frequency f p (t) and the frequency f p (t)+f g , the frequency compression coefficient gradually changes from 1 to a max to ensure that the broadband wavelet has rich mid-frequency information.

[0119] Step S303: Calculate the spatio-temporal variable broadband wavelet amplitude spectrum of the current seismic trace according to the frequency-variable compression coefficient.

[0120] In one implementation, the time-frequency amplitude spectrum of the spatio-temporal variable broadband wavelet of the current seismic trace can be calculated according to the determined frequency-variable compression coefficient by using the cubic spline interpolation algorithm.

[0121] Step S304: Calculate the spatio-temporal variable frequency-up operator corresponding to the current seismic trace based on the spatio-temporal variable wavelet amplitude spectrum and the spatio-temporal variable broadband wavelet amplitude spectrum.

[0122] Step S305: Perform frequency-up processing on the current seismic trace by using the spatio-temporal variable frequency-up operator.

[0123] It should be noted that Figure 5 Steps S301, S303 to S305 in the method embodiment shown are similar to the steps in the Figure 3 method embodiment shown. For the relevant description, reference can be made to the Figure 3 method embodiment shown, and details are not described here again.

[0124] Figure 5 The method embodiment shown has all the beneficial effects of the Figure 3 method embodiment shown. In addition, since the parameter f p (t) in the embodiment of the present invention changes with time and space, the determined frequency-variable compression coefficient a(t,f) also changes adaptively with time and space. It can be seen that the method provided by the embodiment of the present invention can effectively handle the frequency-up processing of complex structures and protect the low-frequency components of seismic data during the processing.

[0125] Example 4

[0126] Such as Figure 6As shown in the figure, it is another flowchart of the method for improving the resolution of seismic data provided by the embodiments of the present invention, which may include the following steps:

[0127] Step S401: Obtain the spatio-temporal variable wavelet amplitude spectrum of the current seismic trace.

[0128] Step S402: Construct the frequency-variable compression coefficient corresponding to the current seismic trace according to the spatio-temporal variable wavelet amplitude spectrum.

[0129] Step S403: Calculate the spatio-temporal variable wide-frequency wavelet amplitude spectrum of the current seismic trace according to the frequency-variable compression coefficient.

[0130] Step S404: Calculate the spatio-temporal variable frequency-up conversion operator corresponding to the current seismic trace based on the spatio-temporal variable wavelet amplitude spectrum and the spatio-temporal variable wide-frequency wavelet amplitude spectrum.

[0131] Step S405: Perform time-frequency domain frequency-up conversion processing on the current seismic trace according to the spatio-temporal variable frequency-up conversion operator, and obtain the time-frequency spectrum after frequency-up conversion processing;

[0132] Step S406: Sum the time-frequency spectrum in the time direction, then perform inverse Fourier transform on the summation result, and use the real part of the inverse Fourier transform as the result of frequency-up conversion processing.

[0133] The following combines specific examples to illustrate the frequency-up conversion processing process for the current seismic trace.

[0134] (1) For each sampling point t, use the frequency-variable compression coefficient to stretch the frequency coordinate axis to obtain a new frequency coordinate axis, and finally use the cubic spline interpolation algorithm to resample the spatio-temporal variable wavelet from the new frequency coordinate axis to the original frequency coordinate axis:

[0135]

[0136] Among them, is the spatio-temporal variable wide-frequency wavelet amplitude spectrum corresponding to the t-th sampling point of the seismic trace with inline number i and xline number j; F new is the new frequency coordinate axis after stretching; F is the original frequency coordinate axis; F new =F·a(t,f) is the relationship between the new frequency coordinate axis, the original frequency coordinate axis, and the low-frequency-preserving frequency-variable compression coefficient.

[0137] (2) After obtaining the spatio-temporal variable wide-frequency wavelet time-frequency amplitude spectrum, calculate the time-frequency domain frequency-up conversion operator of the local seismic pre-trace according to the following expression:

[0138]

[0139] Among them, H i,j(t,f) is the spatio-temporal variable frequency lifting operator corresponding to the t-th sampling point of the current seismic trace with inline number i and xline number j; μ represents the white noise coefficient, and according to the empirical value, it is usually taken as 0.01; ξ represents the maximum value of the spatio-temporal variable wavelet amplitude spectrum at the t-th sampling point.

[0140] (3) According to the following expression, use the spatio-temporal variable frequency lifting operator to perform frequency lifting processing on the current seismic trace:

[0141]

[0142] where is the time-frequency spectrum after frequency lifting processing of the current seismic trace with inline number i and xline number j.

[0143] (4) Sum the time-frequency spectrum after frequency lifting processing in the time direction, perform inverse Fourier transform to the time domain and take the real part, which is the result of frequency lifting processing of the current trace.

[0144]

[0145] where S i,j is the seismic trace after frequency lifting of the current seismic trace with inline number i and xline number j; real represents taking the real part; ifft represents inverse Fourier transform; T represents the number of sampling points.

[0146] Figure 6 Steps S401 to S404 in the method embodiment shown are similar to Figure 2 Steps S101 to S104 in the method embodiment shown, and reference can be made to Figure 2 The relevant description of the method embodiment shown, and details are not described here again.

[0147] Figure 6 The method embodiment shown has Figure 2 All the beneficial effects of the method embodiment shown. In addition, the embodiment of the present invention calculates the frequency-variable compression coefficient with adaptive low-frequency preservation and the spatio-temporal variable broadband wavelet, so as to obtain the spatio-temporal variable frequency lifting operator, and this spatio-temporal variable frequency lifting operator can perform frequency lifting processing adaptively according to the local spectral characteristics of seismic data. Therefore, when applied to the scenario of seismic data with complex structures, it can not only improve the resolution of seismic data, but also not damage the relative relationship of spatio-temporal energy of the original data, showing high amplitude preservation.

[0148] For the above method embodiment, the device for improving the resolution of seismic data provided by the present invention will be further described.

[0149] Example 5

[0150] As Figure 7As shown in the figure, it is the structural diagram of the device for improving the resolution of seismic data provided by the present invention, including the following modules:

[0151] The spatio-temporal variable wavelet amplitude spectrum acquisition module 510 is used to acquire the spatio-temporal variable wavelet amplitude spectrum of the current seismic trace.

[0152] The compression coefficient construction module 520 is used to construct the frequency-variable compression coefficient corresponding to the current seismic trace according to the spatio-temporal variable wavelet amplitude spectrum.

[0153] The wide-band wavelet amplitude spectrum acquisition module 530 is used to calculate the spatio-temporal variable wide-band wavelet amplitude spectrum of the current seismic trace according to the frequency-variable compression coefficient.

[0154] The frequency up-conversion operator calculation module 540 is used to calculate the spatio-temporal variable frequency up-conversion operator corresponding to the current seismic trace based on the spatio-temporal variable wavelet amplitude spectrum and the spatio-temporal variable wide-band wavelet amplitude spectrum.

[0155] The frequency up-conversion processing module 550 is used to perform frequency up-conversion processing on the current seismic trace by using the spatio-temporal variable frequency up-conversion operator.

[0156] In one implementation, the compression coefficient construction module 520 is used to determine the frequency-variable compression coefficient corresponding to the current seismic trace based on the peak frequency in the spatio-temporal variable wavelet amplitude spectrum.

[0157] Specifically, the compression coefficient construction module 520 is specifically used to determine the frequency-variable compression coefficient a(t, f) corresponding to the current seismic trace according to the following expression:

[0158]

[0159] where f is the frequency of the current seismic trace, f p (t) is the peak frequency of the spatio-temporal variable wavelet amplitude spectrum of the t-th sampling point of the current seismic trace, f Nyq is the Nyquist frequency, f g is a frequency constant with a value range of 1 to 2, a max is the maximum compression coefficient of the frequency-variable wavelet.

[0160] In one implementation, the frequency up-conversion processing module 550 is specifically used to perform time-frequency domain frequency up-conversion processing on the current seismic trace according to the spatio-temporal variable frequency up-conversion operator, and obtain the time-frequency spectrum after frequency up-conversion processing; sum the time-frequency spectrum in the time direction, then perform inverse Fourier transform on the summation result, and use the real part of the inverse Fourier transform as the result of frequency up-conversion processing.

[0161] It can be seen that the frequency-varying wavelet based on time and space extracted by the solution of the present invention is more suitable for the requirements of complex structure work areas; the time-space varying frequency-up operator obtained based on the amplitude spectrum of the time-space varying wavelet and the amplitude spectrum of the time-space varying wide-frequency wavelet can adaptively determine the corresponding time-space varying frequency-up operator according to the local spectral characteristics of the seismic data. In this way, even when applied to the seismic data of complex structure work areas, the resolution can be effectively improved; in addition, according to experiments, the solution of the present invention does not damage the relative relationship of the time-space energy of the collected seismic data, so it has high amplitude preservation performance.

[0162] Example 6

[0163] As Figure 8 shown, the time-space varying wavelet amplitude spectrum acquisition module 510 includes a surrounding seismic trace determination unit 511 and a time-space varying wavelet amplitude spectrum acquisition unit 512.

[0164] Among them, the surrounding seismic trace determination unit 511 is used to determine each surrounding seismic trace corresponding to the current seismic trace from the seismic data of the three-dimensional target work area.

[0165] The time-space varying wavelet amplitude spectrum acquisition unit 512 is used to acquire the time-space varying wavelet amplitude spectrum of the current seismic trace according to the average time-frequency amplitude spectra of the surrounding seismic traces.

[0166] Specifically, the surrounding seismic trace determination unit 511 is specifically used to use a spatial sliding window and according to a set search radius value, determine each surrounding seismic trace corresponding to the current seismic trace from the seismic data of the three-dimensional target work area.

[0167] Specifically, the time-space varying wavelet amplitude spectrum acquisition unit 512 includes an average time-frequency amplitude spectrum calculation subunit and a time-space varying wavelet amplitude spectrum acquisition subunit.

[0168] Among them, the average time-frequency amplitude spectrum calculation subunit is used to calculate the time-frequency spectrum of each surrounding seismic trace by using the short-time Fourier transform, and calculate the average time-frequency amplitude spectrum of each surrounding seismic trace based on the time-frequency spectrum;

[0169] The time-space varying wavelet amplitude spectrum acquisition subunit is used to acquire the time-space varying wavelet amplitude spectrum of the current seismic trace according to the average time-frequency amplitude spectrum.

[0170] In one implementation manner, the time-space varying wavelet amplitude spectrum acquisition subunit is specifically used to perform time-frequency domain logarithmic spectrum simulation on the average time-frequency amplitude spectrum to obtain the time-space varying wavelet amplitude spectrum of the current seismic trace.

[0171] Example 7

[0172] To solve the above technical problems, the present invention provides a computer device, asFigure 9 As shown, it includes a memory 610, a processor 620, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the method described above is implemented.

[0173] The computer device may be a computing device such as a desktop computer, a notebook, a palm computer, and a cloud server. The computer device may include, but is not limited to, a processor 620 and a memory 610. Those skilled in the art can understand that Figure 8 merely examples of computer devices, which do not constitute a limitation on computer devices, may include more or fewer components than shown in the figure, or combine certain components, or different components. For example, the computer device may also include input / output devices, network access devices, buses, etc.

[0174] The so-called processor 620 may be a central processing unit (CPU), or may also be other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The general-purpose processor may be a microprocessor or the processor may also be any conventional processor, etc.

[0175] The memory 610 may be an internal storage unit of the computer device, such as the hard disk or memory of the computer device. The memory 610 may also be an external storage device of the computer device, such as a plug-in hard disk, a smart media card (SMC), a secure digital (SD) card, a flash card, etc. equipped on the computer device. Further, the memory 610 may also include both the internal storage unit and the external storage device of the computer device. The memory 610 is used to store the computer program and other programs and data required by the computer device. The memory 610 may also be used to temporarily store data that has been output or will be output.

[0176] The method implemented when the processor executes the computer program includes the following steps:

[0177] Step S101: Obtain the spatio-temporal variable wavelet amplitude spectrum of the current seismic trace.

[0178] Step S102: Construct the frequency-variable compression coefficient corresponding to the current seismic trace according to the spatio-temporal variable wavelet amplitude spectrum.

[0179] Step S103: Calculate the spatio-temporal variable wide-frequency wavelet amplitude spectrum of the current seismic trace according to the frequency-variable compression coefficient.

[0180] Step S104: Calculate the spatio-temporal variable frequency-up operator corresponding to the current seismic trace based on the spatio-temporal variable wavelet amplitude spectrum and the spatio-temporal variable wide-frequency wavelet amplitude spectrum.

[0181] Step S105: Perform frequency-up processing on the current seismic trace by using the spatio-temporal variable frequency-up operator.

[0182] The method implemented when the processor executes the computer program may include the following steps:

[0183] In one case, the method implemented when the processor executes the computer program may include the following steps:

[0184] Step S201: Determine each surrounding seismic trace corresponding to the current seismic trace from the seismic data of the three-dimensional target work area.

[0185] Step S202: Obtain the spatio-temporal variable wavelet amplitude spectrum of the current seismic trace according to the average time-frequency amplitude spectra of the surrounding seismic traces.

[0186] Step S203: Construct the frequency-variable compression coefficient corresponding to the current seismic trace according to the spatio-temporal variable wavelet amplitude spectrum.

[0187] Step S204: Calculate the spatio-temporal variable wide-frequency wavelet amplitude spectrum of the current seismic trace according to the frequency-variable compression coefficient.

[0188] Step S205: Calculate the spatio-temporal variable frequency-up operator corresponding to the current seismic trace based on the spatio-temporal variable wavelet amplitude spectrum and the spatio-temporal variable wide-frequency wavelet amplitude spectrum.

[0189] Step S206: Perform frequency-up processing on the current seismic trace by using the spatio-temporal variable frequency-up operator.

[0190] In another case, the method implemented when the processor executes the computer program may include the following steps:

[0191] Step S301: Obtain the spatio-temporal variable wavelet amplitude spectrum of the current seismic trace.

[0192] Step S302: Determine the frequency-variable compression coefficient corresponding to the current seismic trace based on the peak frequency in the spatio-temporal variable wavelet amplitude spectrum. Determine the frequency-variable compression coefficient a(t, f) corresponding to the current seismic trace according to the following expression

[0193]

[0194] Among them, f is the frequency of the current seismic trace, f p (t) is the peak frequency of the spatio-temporal variable wavelet amplitude spectrum at the t-th sampling point of the current seismic trace, f Nyq is the Nyquist frequency, f g is a frequency constant with a value range between 1 and 2, a max is the maximum compression coefficient of the frequency-variable wavelet. It should be noted that a max is the maximum compression coefficient of the frequency-variable wavelet, and according to empirical values, it is usually sufficient to take a value of 2.

[0195] Step S303: Calculate the spatio-temporal variable wide-frequency wavelet amplitude spectrum of the current seismic trace according to the frequency-variable compression coefficient.

[0196] Step S304: Calculate the spatio-temporal variable frequency-upconversion operator corresponding to the current seismic trace based on the spatio-temporal variable wavelet amplitude spectrum and the spatio-temporal variable wide-frequency wavelet amplitude spectrum.

[0197] Step S305: Perform frequency-upconversion processing on the current seismic trace using the spatio-temporal variable frequency-upconversion operator.

[0198] In another case, the method implemented by the processor when executing the computer program may include the following steps:

[0199] Step S401: Obtain the spatio-temporal variable wavelet amplitude spectrum of the current seismic trace.

[0200] Step S402: Construct the frequency-variable compression coefficient corresponding to the current seismic trace according to the spatio-temporal variable wavelet amplitude spectrum.

[0201] Step S403: Calculate the spatio-temporal variable wide-frequency wavelet amplitude spectrum of the current seismic trace according to the frequency-variable compression coefficient.

[0202] Step S404: Calculate the spatio-temporal variable frequency-upconversion operator corresponding to the current seismic trace based on the spatio-temporal variable wavelet amplitude spectrum and the spatio-temporal variable wide-frequency wavelet amplitude spectrum.

[0203] Step S405: Perform time-frequency domain frequency-upconversion processing on the current seismic trace according to the spatio-temporal variable frequency-upconversion operator, and obtain the time-frequency spectrum after frequency-upconversion processing;

[0204] Step S406: Sum the time-frequency spectrum in the time direction, then perform inverse Fourier transform on the summation result, and use the real part of the inverse Fourier transform as the result of the frequency-upconversion processing.

[0205] For the specific implementation methods and related descriptions of each step, please refer to the foregoing method embodiments, and they will not be elaborated in this embodiment.

[0206] Example 8

[0207] The embodiments of the present application also provide a computer-readable storage medium. The computer-readable storage medium may be the computer-readable storage medium included in the memory in the above embodiments; or it may be a computer-readable storage medium that exists independently and is not assembled into a computer device. The computer-readable storage medium stores one or more computer programs, and when the programs are executed by a processor, the above-described method is implemented.

[0208] If the integrated module / unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, to implement all or part of the processes in the method embodiments of the present application, it can also be completed by instructing relevant hardware through a computer program. The computer program can be stored in a computer-readable storage medium. When the computer program is executed by a processor, the steps of the above method embodiments can be implemented. Among them, the computer program includes computer program code, and the computer program code can be in the form of source code, object code, executable file or some intermediate form, etc. The computer-readable medium may include: any entity or device capable of carrying the computer program code, recording medium, USB flash drive, mobile hard disk, magnetic disk, optical disc, computer memory 610, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), electrical carrier signal, telecommunication signal, and software distribution medium, etc. It should be noted that the content included in the computer-readable medium can be appropriately increased or decreased according to the requirements of legislation and patent practice in the jurisdiction. For example, in some jurisdictions, according to legislation and patent practice, the computer-readable medium does not include electrical carrier signals and telecommunication signals.

[0209] For the device embodiments, since they are basically similar to the method embodiments, the description is relatively simple. For the relevant parts, please refer to the partial description of the method embodiments.

[0210] Those skilled in the art can clearly understand that, for the convenience and brevity of description, only the division of the above-mentioned functional units and modules is used as an example. In actual applications, the above-mentioned functions can be allocated to different functional units and modules according to needs, that is, the internal structure of the device is divided into different functional units or modules to complete all or part of the functions described above. Each functional unit and module in the embodiment can be integrated into a processing unit, or each unit can exist physically alone, or two or more units can be integrated into one unit. The above-mentioned integrated unit can be implemented in the form of hardware or in the form of a software functional unit. In addition, the specific names of the functional units and modules are only for the convenience of mutual distinction and do not limit the protection scope of this application. The specific working processes of the units and modules in the above system can refer to the corresponding processes in the foregoing method embodiments and will not be elaborated herein.

[0211] It should be noted that in this document, 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 terms "comprise", "include" or any other variant thereof are intended to cover non-exclusive inclusion, so 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 phrase "comprising a..." does not exclude the presence of additional identical elements in the process, method, article or device comprising said element.

[0212] It should be understood that the terms used in the specification of this application are only for the purpose of describing specific embodiments and are not intended to limit this application. As used in the specification of this application and the appended claims, unless the context clearly indicates otherwise, the singular forms "a", "an" and "the" are intended to include the plural forms.

[0213] It should also be understood that the term "and / or" used in the specification of this application and the appended claims refers to any combination and all possible combinations of one or more of the associated listed items, and includes these combinations.

[0214] As used in this specification and the appended claims, the term "if" may be construed, depending on the context, as "when" or "once" or "in response to determining" or "in response to detecting". Similarly, the phrase "if determined" or "if the described condition or event is detected" may be construed, depending on the context, to mean "once determined" or "in response to determining" or "once the described condition or event is detected" or "in response to detecting the described condition or event".

[0215] The above are only the preferred embodiments of the present invention and are not intended to limit the scope of protection of the present invention. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of the present invention are included in the scope of protection of the present invention.

Claims

1. A method for improving the resolution of seismic data, characterized in that Comprising: Obtain the spatio-temporal variable wavelet amplitude spectrum of the current seismic trace; Construct the frequency-variable compression coefficient corresponding to the current seismic trace according to the spatio-temporal variable wavelet amplitude spectrum; Calculate the spatio-temporal variable wide-frequency wavelet amplitude spectrum of the current seismic trace according to the frequency-variable compression coefficient; Based on the spatio-temporal variable wavelet amplitude spectrum and the spatio-temporal variable wide-frequency wavelet amplitude spectrum, calculate the spatio-temporal variable frequency-up operator corresponding to the current seismic trace; Use the spatio-temporal variable frequency-up operator to perform frequency-up processing on the current seismic trace.

2. The method for improving the resolution of seismic data according to claim 1, wherein The obtaining of the spatio-temporal variable wavelet amplitude spectrum of the current seismic trace includes: Determine each surrounding seismic trace corresponding to the current seismic trace from the seismic data of the three-dimensional target work area; Obtain the spatio-temporal variable wavelet amplitude spectrum of the current seismic trace according to the average time-frequency amplitude spectrum of each surrounding seismic trace.

3. The method for improving the resolution of seismic data according to claim 2, wherein, The determining of each surrounding seismic trace corresponding to the current seismic trace from the seismic data of the three-dimensional target work area includes: Use a spatial sliding window according to the set search radius value to determine each surrounding seismic trace corresponding to the current seismic trace from the seismic data of the three-dimensional target work area.

4. The method for improving the resolution of seismic data according to claim 2, wherein The obtaining of the spatio-temporal variable wavelet amplitude spectrum of the current seismic trace according to the average time-frequency amplitude spectrum of each surrounding seismic trace includes: Use the short-time Fourier transform to calculate the time-frequency spectrum of each surrounding seismic trace, and calculate the average time-frequency amplitude spectrum of each surrounding seismic trace based on the time-frequency spectrum; Obtain the spatio-temporal variable wavelet amplitude spectrum of the current seismic trace according to the average time-frequency amplitude spectrum.

5. The method for improving the resolution of seismic data according to claim 4, characterized in that, The obtaining of the spatio-temporal variable wavelet amplitude spectrum of the current seismic trace according to the average time-frequency amplitude spectrum includes: Perform time-frequency domain logarithmic spectrum simulation on the average time-frequency amplitude spectrum to obtain the spatio-temporal variable wavelet amplitude spectrum of the current seismic trace.

6. The method for improving the resolution of seismic data according to claim 1, wherein The constructing of the frequency-variable compression coefficient corresponding to the current seismic trace according to the spatio-temporal variable wavelet amplitude spectrum includes: Based on the peak frequency in the spatio-temporal variable wavelet amplitude spectrum, determine the frequency-variable compression coefficient corresponding to the current seismic trace.

7. The method for improving the resolution of seismic data according to claim 6, characterized in that, The determining of the frequency-variable compression coefficient corresponding to the current seismic trace based on the peak frequency in the spatio-temporal variable wavelet amplitude spectrum includes: Determine the frequency-variable compression coefficient a(t,f) corresponding to the current seismic trace according to the following expression where f is the frequency of the current seismic trace, f p (t) is the peak frequency of the spatio-temporal variable wavelet amplitude spectrum at the t-th sampling point of the current seismic trace, f Nyq is the Nyquist frequency, f g is a frequency constant with a value range between 1 and 2, a max is the maximum compression coefficient of the frequency-variable wavelet.

8. The method for improving the resolution of seismic data according to claim 1, wherein The performing of frequency-up processing on the current seismic trace using the spatio-temporal variable frequency-up operator includes: Perform time-frequency domain frequency-up processing on the current seismic trace according to the spatio-temporal variable frequency-up operator, and obtain the time-frequency spectrum after frequency-up processing; Sum the time-frequency spectrum in the time direction, then perform inverse Fourier transform on the summation result, and use the real part of the inverse Fourier transform as the result of frequency-up processing.

9. A device for improving the resolution of seismic data, characterized in that, Comprising: A spatio-temporal variable wavelet amplitude spectrum obtaining module, configured to obtain the spatio-temporal variable wavelet amplitude spectrum of the current seismic trace; A compression coefficient constructing module, configured to construct the frequency-variable compression coefficient corresponding to the current seismic trace according to the spatio-temporal variable wavelet amplitude spectrum; A wide-frequency wavelet amplitude spectrum obtaining module, configured to calculate the spatio-temporal variable wide-frequency wavelet amplitude spectrum of the current seismic trace according to the frequency-variable compression coefficient; A frequency-up operator calculating module, configured to calculate the spatio-temporal variable frequency-up operator corresponding to the current seismic trace based on the spatio-temporal variable wavelet amplitude spectrum and the spatio-temporal variable wide-frequency wavelet amplitude spectrum; The frequency-up conversion processing module is used to perform frequency-up conversion processing on the current seismic trace by using the spatio-temporal variable frequency-up conversion operator.

10. A computer 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 computer program, the method described in any one of claims 1 to 8 is implemented.

11. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the program is executed by the processor, the method described in any one of claims 1 to 8 is implemented.

Images