CN111487675A - Method for generating seismic data high signal-to-noise ratio and high resolution time frequency spectrum - Google Patents

Method for generating seismic data high signal-to-noise ratio and high resolution time frequency spectrum Download PDF

Info

Publication number
CN111487675A
CN111487675A CN202010215366.6A CN202010215366A CN111487675A CN 111487675 A CN111487675 A CN 111487675A CN 202010215366 A CN202010215366 A CN 202010215366A CN 111487675 A CN111487675 A CN 111487675A
Authority
CN
China
Prior art keywords
time
spectrum
signal
seismic
frequency spectrum
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010215366.6A
Other languages
Chinese (zh)
Other versions
CN111487675B (en
Inventor
王仰华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Priority to CN202010215366.6A priority Critical patent/CN111487675B/en
Publication of CN111487675A publication Critical patent/CN111487675A/en
Application granted granted Critical
Publication of CN111487675B publication Critical patent/CN111487675B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/307Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention relates to a method for generating a high-signal-to-noise-ratio high-resolution time frequency spectrum of seismic data, which is realized by multi-channel weighted synthesis along a formation dip angle and a mode of constructing an adaptive mask filter; the method comprises the following steps: estimating dip angles of each time-space position on the seismic section; constructing a plurality of weighted synthesis channels aiming at each time position, and calculating an instantaneous autocorrelation function of an analytic signal of the synthesis channels; calculating standard Weiwei time-frequency distribution and corresponding flat-virtual Weiwei distribution; constructing a time-frequency spectrum self-adaptive mask filter; and realizing mask filtering of standard Weiwei time-frequency distribution to obtain an optimized seismic channel time-frequency spectrum. The method has the characteristics of high signal-to-noise ratio and high resolution, simultaneously has the capability of eliminating cross interference, and has the advantage of clearly separating the spatial distribution of the sandstone reservoir and the spatial distribution of the system containing the coal bed when the actual three-dimensional seismic data is used for reservoir prediction.

Description

Method for generating seismic data high signal-to-noise ratio and high resolution time frequency spectrum
Technical Field
The invention relates to the field of geophysical signal analysis, in particular to a method for generating a seismic data time frequency spectrum with high signal-to-noise ratio and high resolution.
Background
Time-frequency spectrum analysis of seismic data is a key means for reservoir geophysical research. Intuitively, the time-frequency spectrum analysis is to generate corresponding frequency spectrum point by point according to the time sampling sequence of any seismic single channel, and is to convert a one-dimensional seismic channel changing along with time into a two-dimensional spectrum changing along with time and frequency. This two-dimensional spectrum is called the time spectrum.
The Wigner-Ville time frequency distribution method is one of effective methods for calculating the time frequency spectrum of the seismic trace. The implementation process can be divided into three basic steps: calculating analytic signals of the seismic traces, calculating an instantaneous autocorrelation function of the analytic signals, and calculating Fourier transform of the autocorrelation function. Because the analytical signal and the instantaneous autocorrelation function are adopted in the calculation method of the Weiwei time-frequency distribution, the Weiwei time-frequency distribution has the characteristic of high resolution. However, it is because the computation of the instantaneous autocorrelation function leads to severe cross-interference between the components contained in the signal.
In order to suppress the cross-interference between the above-mentioned signal components, researchers have proposed a calculation method called the flat-imaginary wiry distribution. It contains two layers of meaning, one "smooth" refers to the suppression of interference in the time direction, and one "virtual" actually refers to the suppression of interference in the frequency direction. However, the suppression processing in both the time direction and the frequency direction is essentially low-pass filtering, and this has the effect of losing high-frequency components and reducing the resolution of the time spectrum.
Because the time frequency spectrum calculation method is a single-channel analysis method, the generated time frequency spectrum is influenced by the signal-to-noise ratio of seismic data, the signal-to-noise ratio of the time frequency spectrum is low, and the time frequency spectrum between channels lacks space consistency, thereby causing difficulty in the interpretation and application of the seismic time frequency spectrum.
By combining the analysis, the spectrum resolution is high but the cross interference is serious in the standard Weiwei time spectrum, the cross interference is inhibited but the resolution is reduced in the flat-virtual Weiwei time spectrum, and the calculation of the time spectrum is single-channel operation, so that the generated earthquake time spectrum has poor signal-to-noise ratio and is seriously influenced by earthquake data noise, and the application of the time spectrum and the geophysical explanation are difficult.
Disclosure of Invention
The invention aims to overcome the defects of the existing method and provides a method for generating a high-signal-to-noise-ratio high-resolution time spectrum of seismic data. The method for generating the earthquake time frequency spectrum has the advantages of high signal-to-noise ratio and high resolution.
In order to achieve the purpose, the seismic data high signal-to-noise ratio high resolution time frequency spectrum generation method provided by the invention improves the signal-to-noise ratio of the seismic time frequency spectrum by multi-channel weighted synthesis along the formation dip angle and utilizing the spatial continuity of seismic reflection homophase axes; the resolution of the frequency spectrum during earthquake is enhanced by a method of constructing an adaptive mask filter. The method has the advantages of high resolution and high signal-to-noise ratio.
As optimization, the seismic data high signal-to-noise ratio high resolution time frequency spectrum generation method specifically comprises the following steps: (1) estimating dip angles of each time-space position on the seismic section; (2) constructing a plurality of weighted synthesis channels aiming at each time position, and calculating an instantaneous autocorrelation function of an analytic signal of the synthesis channels; (3) calculating standard Weiwei (Wigner-Ville) time-frequency distribution and corresponding flat-virtual Weiwei distribution; (4) constructing a time-frequency spectrum self-adaptive mask filter; (5) and realizing mask filtering of standard Weiwei time-frequency distribution to obtain an optimized seismic channel time-frequency spectrum.
As optimization, the step (1) of estimating the dip angle of each time-space position from the actual seismic section is to construct a nonlinear equation set by adopting a discrete plane wave analytic equation mode, and estimate the time-varying and space-varying stratigraphic dip angles by a calculation method for solving a nonlinear problem.
Estimating the dip angle of each time-space position on the seismic section by adopting a plane wave analytic equation:
Figure 100002_DEST_PATH_IMAGE001
,
in the formulatIs a time sample of the seismic trace,xis the spatial location of the seismic traces,u(t,x) Is the seismic wavefield and σ is the dip of the linear event. Discretizing the plane wave analytic equation and expressing the plane wave analytic equation in a vector matrix form
f(σ) = e T [C(σ)○U]e = 0,
Where U is the data matrix, C (σ) is the matrix composed of the coefficients of the two-dimensional filter in the reverse order, [ C (σ) ○ U]Is a multiplication by elements or a Hadamard product of two matrices, e is a full 1 vector. The matrix [ C ] is realized by introducing all 1 vectors e in the formula(σ)○U]The sum of all elements in. By solving a system of non-linear equationsf(σ) ≈ 0 estimates time-varying and space-varying formation dip angles.
As an optimization, step (2) implements the following three steps at each time position: constructing a composite track by performing multi-channel weighting according to local dip anglesx(t) (ii) a Calculating an analytic signal of the composite tracez(t) (ii) a Calculating the instantaneous autocorrelation function of the analytic signal:k(t,τ) =z(t+τ/2)z*(t-τ/2) in the formulaτIs the lag time of the instantaneous auto-correlation,z*(t) Is analyzing the signalz(t) The complex number of (c) is conjugated.
As optimization, step (3) uses the instantaneous autocorrelation function obtained in step (2) to calculate the standard wiry time-frequency distribution as follows:
Figure 418894DEST_PATH_IMAGE002
,
and the corresponding Ping-Xunwei distribution was calculated as follows:
Figure 100002_DEST_PATH_IMAGE003
in the formulah(τ) Is a symmetric function with time as a variable, which functions as a low-pass filter along the frequency axis,g(t) Is another function of time as a variable, which functions as a low-pass filter along the time axis, andg(0) =h(0)。
as optimization, the mask filter in step (4) is constructed as follows:
Figure 835356DEST_PATH_IMAGE004
,
the middle ratio ofW SP (t,f)|2/ |W(t,f)|2Is the similarity of the two power spectra,is an adjustment parameter that controls the position of the mask filter boundary,ηthen it is the parameter that controls the smoothness of the mask filter boundary, settingηNot less than 2. For the sake of stability of the estimation, int,f) The power spectrum estimate at a point is the sum of the power spectra within a small two-dimensional window centered at that point.
As optimization, the constructed mask filter has the characteristics of the flattest passband and the adjusting parameters thereofEffectively controlling the boundary position of the filter; setting a boundary smoothness parameterη= 4 regulating parameterWhen = 0.05, any imaging spectral point with a power spectrum ratio value less than 20% can be filtered out from the time-frequency spectrum.
As optimization, the implementation process of the step (5) is divided into two steps: distribution of amplitude spectrum of Weiwei time-frequency distributionW(t,f) The l is filtered out and the filter is implemented,
Figure 250550DEST_PATH_IMAGE005
(ii) a The wiry time frequency spectrum is reconstructed,
Figure 932199DEST_PATH_IMAGE006
in the formulaƟ(t,f) The method is an original phase spectrum of Weiwei time-frequency distribution, and finally obtains a seismic channel optimized time-frequency spectrum with high signal-to-noise ratio and high resolution.
After the embodiment is adopted, the invention has the following beneficial effects: (1) the method has the advantages that the method has high resolution, the generated time frequency spectrum has the same high resolution as that of standard Weiwei distribution, and the method has the same capability of eliminating cross interference as that of flat-virtual Weiwei distribution. (2) The method has the advantage of high signal-to-noise ratio, and can clearly separate the spatial distribution of the sandstone reservoir and the spatial distribution of the coal-containing bed when the actual three-dimensional seismic data is used for reservoir prediction.
Drawings
FIG. 1 is a flow chart of a basic implementation of the seismic data high signal-to-noise ratio and high resolution time-frequency spectrum generation method of the present invention; FIG. 2 is a diagram of an actual seismic section (top) and its estimated dip (bottom) according to an embodiment of the present invention.
Fig. 3 is a mask filter proposed by the present invention, having the flattest passband characteristic. The horizontal axis in the diagram is the ratio of two power spectra +W SP (t,f)|2/|W(t,f)|2Calculating the boundary smoothness parameterηThe value is 4, and the adjusting parameters corresponding to the coefficient curves of the three mask filtersThe values are 0.1,0.05,0.01, respectively.
FIG. 4 is a diagram of an embodiment of the present invention applied to a composite trace containing 4 wavelets comparing the standard Weir distribution of the composite trace, the Ping-Vwei distribution of the composite trace, and the time-frequency spectrum generated by the method of the present invention. Fig. 4 shows that the optimized temporal spectrum generated by the method of the present invention has the same high resolution as the standard wiry distribution, while the method has the same cross-interference rejection capability as the flat-imaginary wiry distribution.
FIG. 5 is an embodiment of the present invention applied to actual three-dimensional seismic recordings for reservoir prediction. The top plot is a time slice of a three-dimensional seismic data volume. The spectrum slice of the middle graph reflects the spatial distribution of the coal measure strata of the region, and the primary frequency of the reflection earthquake of the coal measure strata of the region is 15 Hz. The spectral slice of the lower graph is the spatial distribution of the sandstone reservoir in the region, which reflects seismic dominant frequencies of 20 hertz. The figure shows the high signal-to-noise ratio characteristic of the method of the invention, enabling a clear separation of the spatial distribution of the sandstone reservoir and the coal-derived formation.
Detailed Description
The method for generating the seismic data high signal-to-noise ratio and high resolution time frequency spectrum comprises the steps of improving the signal-to-noise ratio of the seismic time frequency spectrum by means of multi-channel weighted synthesis along the dip angle of the stratum and utilizing the spatial continuity of seismic reflection in-phase axes; the resolution of the frequency spectrum during earthquake is enhanced by a method of constructing an adaptive mask filter. The method comprises the following specific steps: (1) estimating dip angles of each time-space position on the seismic section; (2) constructing a plurality of weighted synthesis channels aiming at each time position, and calculating an instantaneous autocorrelation function of an analytic signal of the synthesis channels; (3) calculating standard Weiwei (Wigner-Ville) time-frequency distribution and corresponding flat-virtual Weiwei distribution; (4) constructing a time-frequency spectrum self-adaptive mask filter; (5) and realizing mask filtering of standard Weiwei time-frequency distribution to obtain an optimized seismic channel time-frequency spectrum.
In the step (1), the estimation of the dip angle of each time-space position from the actual seismic section is to construct a nonlinear equation set by discretizing a plane wave analytic equation, and estimate the dip angles of the time-varying and space-varying stratums by a calculation method for solving a nonlinear problem. Estimating the dip angle of each time-space position on the seismic section by adopting a plane wave analytic equation:
Figure 749238DEST_PATH_IMAGE001
in the formulatIs the point in time of the seismic trace,xis the spatial location of the seismic traces,u(t,x) Is the seismic wavefield and σ is the dip of the linear event. Discretizing the plane wave analytic equation and expressing the plane wave analytic equation in a vector matrix formf(σ)=e T [C(σ)○U]e =0, where U is the data matrix and C (σ) is the matrix composed of the coefficients of the two-dimensional filter in the reverse order, with dimensions identical to the data matrix U, [ C (σ) ○ U]Is element-wise multiplication or Hadamard multiplication of two matrices, e being a full "1" vector, introduction of the full "1" vector e in the above formula achieves a pair of matrices [ C (σ) ○ U]The sum of all elements in. By solving a system of non-linear equationsf(σ) ≈ 0 estimates time-varying and space-varying formation dip angles.
The step (2) realizes the following three steps at each time position: constructing a composite track by performing multi-channel weighting according to local dip anglesx(t) (ii) a Calculating an analytic signal of the composite tracez(t) (ii) a Calculating the instantaneous autocorrelation function of the analytic signal:k(t,τ)=z(t+τ/2)z*(t-τ/2) In the formulaτIs the lag time of the autocorrelation function,z*(t) Is analyzing the signalz(t) The complex number of (c) is conjugated.
And (3) calculating the standard Weiwei time-frequency distribution by using the autocorrelation function in the step (2) as follows:
Figure 819146DEST_PATH_IMAGE007
,
and the corresponding Ping-Xunwei distribution was calculated as follows:
Figure 550997DEST_PATH_IMAGE003
,
in the formulah(τ) Is a symmetric function with time as a variable, which functions as a low-pass filter along the frequency axis,g(t) Is another function of time as a variable, which functions as a low-pass filter along the time axis, andg(0) =h(0)。
the mask filter in the step (4) is constructed as follows:
Figure 475484DEST_PATH_IMAGE004
,
the middle ratio ofW SP (t,f)|2/|W(t,f)|2Is the similarity of the two power spectra,is an adjustment parameter that controls the position of the mask filter boundary,ηthen it is the parameter that controls the smoothness of the mask filter boundary, settingηNot less than 2. For the sake of stability, int,f) The amplitude estimate at a point is the sum of the amplitudes within a small two-dimensional window centered at that point. The constructed mask filter has the characteristics of the flattest passband and the adjusting parameterEffectively controlling the boundary position of the filter; setting a boundary smoothness parameterηAdjustment parameter with value 4When the value is 0.05, any imaging spectrum point with the power spectrum ratio of less than 20 percent can be filtered from the time-frequency spectrum.
The implementation process of the step (5) is divided into two steps: distribution of amplitude spectrum of Weiwei time-frequency distributionW(t,f) The l is filtered out and the filter is implemented,
Figure 114407DEST_PATH_IMAGE005
(ii) a The wiry time frequency spectrum is reconstructed,
Figure 590781DEST_PATH_IMAGE006
in the formulaƟ(t,f) The method is an original phase spectrum of Weiwei time-frequency distribution, and finally obtains a seismic channel optimized time-frequency spectrum with high signal-to-noise ratio and high resolution.
In order to make the technical scheme and advantages of the invention more apparent, the following detailed description of the embodiments of the invention is provided in conjunction with the accompanying drawings. The invention provides a method for generating a high-signal-to-noise-ratio high-resolution time spectrum of seismic data. FIG. 1 is a flow chart of a basic implementation of the method of the present invention.
Step S101, estimating the dip angle of each time-space position on the seismic section. Using a plane wave analytic equation:
Figure 918251DEST_PATH_IMAGE008
in the formulatIs the point in time of the seismic trace,xis the spatial location of the seismic traces,u(t,x) Is the seismic wavefield and σ is the dip of the linear event. Discretizing the plane wave analytic equation and expressing the plane wave analytic equation in a vector matrix form
f(σ) = e T [C(σ)○U]e = 0
Where U is the data matrix, C (σ) is the matrix composed of the coefficients of the two-dimensional filter in the reverse order, [ C (σ) ○ U]Is element-wise multiplication or Hadamard product of two matrices, e is a full 1 vector, and introducing the full 1 vector e in the above formula realizes the pair of matrices [ C (sigma) ○ U]The sum of all elements in. By solving a system of non-linear equationsf(σ) ≈ 0 estimates time-varying and space-varying formation dip angles.
FIG. 2 is an actual seismic section (top) and its estimated dip (bottom) of an embodiment of the invention. The arrows in the figure are the estimated local dip directions, coinciding with the true in-phase axes of the actual seismic event of the background.
Step S102, the following three steps are realized at each time position: constructing a composite track by performing multi-channel weighting according to local dip anglesx(t) (ii) a Calculating an analytic signal of the composite tracez(t) (ii) a Calculating the instantaneous autocorrelation function of the analytic signal:k(t,τ)=z(t+τ/2)z*(t-τ/2) In the formulaτIs the lag time of the autocorrelation function,z*(t) Is analyzing the signalz(t) The complex number of (c) is conjugated.
Step S103, calculating the standard Weiwei time-frequency distribution by using the instantaneous autocorrelation function obtained in the step S102 as follows:
Figure 879385DEST_PATH_IMAGE007
,
and the corresponding Ping-Xunwei distribution was calculated as follows:
Figure 44918DEST_PATH_IMAGE003
in the formulah(τ) Is a symmetric function with time as a variable, which functions as a low-pass filter along the frequency axis,g(t) Is another function of time as a variable, which functions as a low-pass filter along the time axis, andg(0) =h(0)。
step S104, the mask filter is constructed as follows:
Figure 692193DEST_PATH_IMAGE004
,
the middle ratio ofW SP (t,f)|2/|W(t,f)|2Is the similarity of the two power spectra,is an adjustment parameter that controls the position of the mask filter boundary,ηthen it is the parameter that controls the smoothness of the mask filter boundary, settingηNot less than 2. For the sake of stability, int,f) The amplitude estimate at a point is the sum of the amplitudes within a small two-dimensional window centered at that point. As shown in FIG. 3, the constructed mask filter has the flattest passband characteristic, which adjusts the parametersEffectively controlling the boundaries of the filter; setting a boundary smoothness parameterηAdjustment parameter with value 4When the value is 0.05, any imaging spectrum point with the power spectrum ratio of less than 20 percent can be filtered from the time-frequency spectrum.
And step S105, realizing mask filtering of the standard Weiwei time-frequency distribution. The implementation process comprises two steps: distribution of amplitude spectrum of Weiwei time-frequency distributionW(t,f) The l is filtered out and the filter is implemented,
Figure 38117DEST_PATH_IMAGE005
(ii) a The wiry time frequency spectrum is reconstructed,
Figure 537363DEST_PATH_IMAGE006
in the formulaƟ(t,f) The method is an original phase spectrum of Weiwei time-frequency distribution, and finally obtains a seismic channel optimized time-frequency spectrum with high signal-to-noise ratio and high resolution.
The invention also provides two embodiments [ fig. 4, fig. 5], which illustrate that the invention has the following beneficial effects:
(1) as shown in FIG. 4, a composite trace contains 4 wavelets, and in its standard Weir distribution, the cross-interference between wavelets is very significant except for the energy blob corresponding to 4 wavelets; on the flat-virtual weiwei distribution of the composite trace, although the cross interference is completely eliminated, 4 effective energy blobs are obviously spread, so that the resolution of the time spectrum is lowered; finally, the time frequency spectrum generated by the method has the high resolution as that of the standard Weiwei distribution, and meanwhile, the method has the capability of eliminating cross interference as that of the flat-virtual Weiwei distribution.
(2) FIG. 5 is an embodiment of the present invention applied to actual three-dimensional seismic recordings for reservoir prediction. The top plot is a time slice of a three-dimensional seismic data volume. The spectral slice of the middle graph reflects the spatial distribution of the coal measure strata of the region, and the primary frequency of the reflection earthquake of the coal measure strata of the region is 15 Hz. Note that the coal measure stratigraphic space variation of the middle graph has similarities with some parts of the upper graph. The spectral slice of the lower graph is the spatial distribution of the sandstone reservoir in the region, and the reflection seismic main frequency of the sandstone reservoir in the region is 20 Hz. As shown in fig. 5, the method has the advantage of high signal-to-noise ratio, and can clearly separate the spatial distribution of the sandstone reservoir and the spatial distribution of the coal-series stratum when the actual three-dimensional seismic data is used for reservoir prediction.

Claims (8)

1. A method for generating a high-signal-to-noise-ratio high-resolution time spectrum of seismic data is characterized in that the signal-to-noise ratio of the seismic time spectrum is improved by multi-channel weighted synthesis along the dip angle of a stratum by utilizing the spatial continuity of seismic reflection event axes; the resolution of the frequency spectrum during earthquake is enhanced by a method of constructing an adaptive mask filter.
2. The method for generating the seismic data time spectrum with high signal-to-noise ratio and high resolution according to claim 1, characterized by the following steps: (1) estimating dip angles of each time-space position on the seismic section; (2) constructing a plurality of weighted synthesis channels aiming at each time position, and calculating an instantaneous autocorrelation function of an analytic signal of the synthesis channels; (3) calculating standard Weiwei (Wigner-Ville) time-frequency distribution and corresponding flat-virtual Weiwei distribution; (4) constructing a time-frequency spectrum self-adaptive mask filter; (5) and realizing mask filtering of standard Weiwei time-frequency distribution to obtain an optimized seismic channel time-frequency spectrum.
3. The method for generating high SNR and high resolution time-frequency spectrum of seismic data as claimed in claim 2, wherein the step (1) of estimating the dip angle of each time-space position from the actual seismic section is to construct a nonlinear equation system by using a discrete plane wave analytic equation, and estimate the time-varying and space-varying stratigraphic dip angles by solving a nonlinear problem.
4. The method for generating a high signal-to-noise ratio and high resolution time spectrum of seismic data as claimed in claim 2, wherein the step (2) implements the following three steps at each time position: building a composite track by multi-track weighting according to local dip anglesx(t) (ii) a Calculating an analytic signal of the composite tracez(t) (ii) a Calculating the instantaneous autocorrelation function of the analytic signal:k(t,τ)=z(t+τ/2)z*(t-τ/2) in the formulaτIs the lag time of the instantaneous auto-correlation,z*(t) Is analyzing the signalz(t) The complex number of (c) is conjugated.
5. The method for generating the seismic data time-frequency spectrum with high signal-to-noise ratio and high resolution according to claim 2, wherein the step (3) uses the instantaneous autocorrelation function obtained in the step (2) to calculate the standard wiry time-frequency distribution as follows:
Figure DEST_PATH_IMAGE001
,
and the corresponding Ping-Xunwei distribution was calculated as follows:
Figure 244160DEST_PATH_IMAGE002
in the formulah(τ) Is a symmetric function with time as a variable, which functions as a low-pass filter along the frequency axis,g(t) Is another function of time as a variable, which functions as a low-pass filter along the time axis, andg(0) =h(0)。
6. the method for generating the seismic data time spectrum with high signal-to-noise ratio and high resolution as claimed in claim 2, wherein the mask filter in the step (4) is constructed as follows:
Figure DEST_PATH_IMAGE003
the middle ratio ofW SP (t,f)|2/ |W(t,f)|2Is the similarity of the two power spectra,is an adjustment parameter that controls the position of the mask filter boundary,ηthen it is the parameter that controls the smoothness of the mask filter boundary, settingη≥ 2。
7. The method of claim 6, wherein the mask filter is constructed to have the flattest passband characteristics and the tuning parameters thereofEffectively controlling the boundary position of the mask filter; setting a boundary smoothness parameterη= 4 regulating parameterWhen = 0.05, any imaging spectral point with a power spectrum ratio value less than 20% can be filtered out from the time-frequency spectrum.
8. The method of claim 2The method for generating the seismic data time spectrum with high signal-to-noise ratio and high resolution is characterized in that the implementation process of the step (5) is divided into two steps: distribution of amplitude spectrum of Weiwei time-frequency distributionW(t,f) The l is filtered out and the filter is implemented,
Figure 30106DEST_PATH_IMAGE004
(ii) a The wiry time frequency spectrum is reconstructed,
Figure DEST_PATH_IMAGE005
in the formulaƟ(t,f) The method is an original phase spectrum of Weiwei time-frequency distribution, and finally obtains a seismic channel optimized time-frequency spectrum with high signal-to-noise ratio and high resolution.
CN202010215366.6A 2020-03-25 2020-03-25 Method for generating seismic data high signal-to-noise ratio and high resolution time frequency spectrum Active CN111487675B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010215366.6A CN111487675B (en) 2020-03-25 2020-03-25 Method for generating seismic data high signal-to-noise ratio and high resolution time frequency spectrum

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010215366.6A CN111487675B (en) 2020-03-25 2020-03-25 Method for generating seismic data high signal-to-noise ratio and high resolution time frequency spectrum

Publications (2)

Publication Number Publication Date
CN111487675A true CN111487675A (en) 2020-08-04
CN111487675B CN111487675B (en) 2021-08-27

Family

ID=71812363

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010215366.6A Active CN111487675B (en) 2020-03-25 2020-03-25 Method for generating seismic data high signal-to-noise ratio and high resolution time frequency spectrum

Country Status (1)

Country Link
CN (1) CN111487675B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112835100A (en) * 2021-01-21 2021-05-25 中国石油天然气股份有限公司 Spectrum decomposition method and device
CN114167499A (en) * 2021-12-02 2022-03-11 大庆油田有限责任公司 Method and device for automatically identifying biological reef, electronic equipment and storage medium

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102305940A (en) * 2011-05-24 2012-01-04 中国石油集团川庆钻探工程有限公司 Method for extracting fluid factor
CN102466819A (en) * 2010-11-03 2012-05-23 中国石油天然气集团公司 Spectrum analysis method of seismic signal and apparatus thereof
US20140297188A1 (en) * 2013-03-29 2014-10-02 Cgg Services Sa Time-frequency representations of seismic traces using wigner-ville distributions
CN105353408A (en) * 2015-11-20 2016-02-24 电子科技大学 Wigner higher-order spectrum seismic signal spectral decomposition method based on matching pursuit
CN105954801A (en) * 2016-06-12 2016-09-21 成都理工大学 Frequency-conversion-attribute-based reservoir permeability evaluation method
CN106483555A (en) * 2016-09-28 2017-03-08 中国地质大学(武汉) The Green's function SPWVD Time-Frequency Analysis Method of ENPEMF data
CN110865410A (en) * 2019-11-29 2020-03-06 山西大学 Method for suppressing random noise of seismic exploration based on NAR-TFPF

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102466819A (en) * 2010-11-03 2012-05-23 中国石油天然气集团公司 Spectrum analysis method of seismic signal and apparatus thereof
CN102305940A (en) * 2011-05-24 2012-01-04 中国石油集团川庆钻探工程有限公司 Method for extracting fluid factor
US20140297188A1 (en) * 2013-03-29 2014-10-02 Cgg Services Sa Time-frequency representations of seismic traces using wigner-ville distributions
CN105353408A (en) * 2015-11-20 2016-02-24 电子科技大学 Wigner higher-order spectrum seismic signal spectral decomposition method based on matching pursuit
CN105954801A (en) * 2016-06-12 2016-09-21 成都理工大学 Frequency-conversion-attribute-based reservoir permeability evaluation method
CN106483555A (en) * 2016-09-28 2017-03-08 中国地质大学(武汉) The Green's function SPWVD Time-Frequency Analysis Method of ENPEMF data
CN110865410A (en) * 2019-11-29 2020-03-06 山西大学 Method for suppressing random noise of seismic exploration based on NAR-TFPF

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
唐俊: "几种时频分析方法的对比分析", 《价值工程》 *
王雨青: "地震信号分数域属性分析及储层流体识别关键技术研究", 《中国博士学位论文全文数据库 基础科学辑》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112835100A (en) * 2021-01-21 2021-05-25 中国石油天然气股份有限公司 Spectrum decomposition method and device
CN114167499A (en) * 2021-12-02 2022-03-11 大庆油田有限责任公司 Method and device for automatically identifying biological reef, electronic equipment and storage medium
CN114167499B (en) * 2021-12-02 2022-08-26 大庆油田有限责任公司 Method and device for automatically identifying biological reef, electronic equipment and storage medium

Also Published As

Publication number Publication date
CN111487675B (en) 2021-08-27

Similar Documents

Publication Publication Date Title
CN111487675B (en) Method for generating seismic data high signal-to-noise ratio and high resolution time frequency spectrum
Cai et al. Seismic data denoising based on mixed time-frequency methods
CN109031422A (en) A kind of seismic signal noise suppressing method based on CEEMDAN and Savitzky-Golay filtering
CN107144879A (en) A kind of seismic wave noise-reduction method combined based on adaptive-filtering with wavelet transformation
CN109946740B (en) Seismic resolution enhancement method based on wide flat spectrum seismic wavelet shaping
CN103645507A (en) A processing method for seismic records
Liu et al. Seismic quality factor estimation using frequency-dependent linear fitting
CN111474582B (en) Precise S transformation method for generating high-precision time frequency spectrum
CN112305612B (en) High-resolution complex spectrum decomposition time-frequency space domain amplitude variation correction method along with offset distance
CN109901224B (en) Method for protecting and suppressing noise of low-frequency signal of seismic data
CN113281809B (en) Spectrum analysis method for seismic signals
CN112255690B (en) Self-adaptive surrounding rock strong reflection separation method based on seismic phase decomposition
CN116859457A (en) Space-time-frequency wave number domain synchronous extrusion superposition tight river reservoir characterization method
CN104932008B (en) The method that the multiple time-frequency spectrum of compensation J conversion improves seismic profile resolution ratio
CN114371505A (en) Multi-wavelet inversion method and system based on seismic frequency division technology
CN112213775B (en) Fidelity frequency-boosting method for high-coverage-frequency pre-stack seismic data
Porsani et al. An adaptive local-slope SVD filtering approach to enhance events on seismic sections
CN114137606A (en) Stable spectrum simulation deconvolution method
CN112925013A (en) Seismic data high-resolution processing method based on full-band continuation fidelity
CN113341462A (en) Zero-phase processing method for marine seismic data
CN110673211A (en) Quality factor modeling method based on logging and seismic data
CN113359187B (en) Wavelet sidelobe elimination method for seismic data
DONG et al. Fast implementation technique of adaptive Kalman filtering deconvolution via dyadic wavelet transform
King et al. The role of interpolation in seismic resolution
Zheng et al. Time-Frequency Domain Deconvolution based on Synchrosqueezing Generalized S Transform

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant