CN113296155A - Seismic reservoir prediction method by synchronous extraction of base stretch frequency modulation - Google Patents

Seismic reservoir prediction method by synchronous extraction of base stretch frequency modulation Download PDF

Info

Publication number
CN113296155A
CN113296155A CN202110558646.1A CN202110558646A CN113296155A CN 113296155 A CN113296155 A CN 113296155A CN 202110558646 A CN202110558646 A CN 202110558646A CN 113296155 A CN113296155 A CN 113296155A
Authority
CN
China
Prior art keywords
frequency
time
seismic
signal
sbct
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.)
Pending
Application number
CN202110558646.1A
Other languages
Chinese (zh)
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.)
Chengdu Univeristy of Technology
Original Assignee
Chengdu Univeristy of Technology
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 Chengdu Univeristy of Technology filed Critical Chengdu Univeristy of Technology
Priority to CN202110558646.1A priority Critical patent/CN113296155A/en
Publication of CN113296155A publication Critical patent/CN113296155A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/301—Analysis for determining seismic cross-sections or geostructures
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G01V2210/624—Reservoir parameters

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 discloses a prediction method for synchronously extracting seismic reservoirs by using fundamental stretch frequency modulation, which comprises the following steps of: s1, inputting an original two-dimensional seismic section signal S (x; t) to be analyzed; s2, decomposing the original two-dimensional seismic section signal S (x; t) by using basis-extension adaptive frequency modulation transform (SBCT) to obtain a time-frequency transform result SBCT (x; f, t); s3, estimating the instantaneous frequency of the x-th seismic signal at each time-frequency position according to the time-frequency transformation result SBCT (x; f, t); s4, according to the synchronous extraction principle, constructing a synchronous extraction operator taking the signal instantaneous frequency as the center of a curve on the time-frequency domain, extracting the original time-frequency spectrum energy to the curve to obtain a new time-frequency coefficient, and removing fuzzy time-frequency energy; s5, taking a module of the time frequency coefficient to obtain a time frequency spectrum after base stretch frequency modulation synchronous extraction transformation. S6, calculating the dominant frequency range of the original seismic signal by adopting Fourier transform, and determining a high frequency value and a low frequency value; then respectively extracting frequency spectrums corresponding to the high-frequency value and the low-frequency value from the time frequency spectrum obtained in the step S5 to obtain two common frequency profiles; and finally, determining whether the reservoirs exist in the seismic section by comparing the attenuation conditions of the seismic signals in the two common frequency sections. The method can obviously improve the time-frequency energy focusing property and improve the identification precision of the seismic reservoir.

Description

Seismic reservoir prediction method by synchronous extraction of base stretch frequency modulation
Technical Field
The invention belongs to the field of signal processing, and particularly relates to a prediction method for synchronously extracting seismic reservoirs through fundamental stretching frequency modulation.
Background
Time-frequency analysis is an important means for earthquake reservoir prediction, can effectively represent the non-stationary characteristics of earthquake signals and reveal the change relation of signal frequency along with time. As an effective means of signal processing, time-frequency analysis clearly describes the close relationship between signal frequency and time, and is therefore widely used in seismic reservoir identification. Many scholars realize oil and gas reservoir identification by performing spectrum analysis on different frequency slices of a geologic body and utilizing the 'absorption and attenuation' law of instantaneous spectral energy of the different frequency slices. In recent years, in order to research complex and variable non-stationary seismic signals and improve the characterization accuracy of the signals, scholars successively put forward a plurality of time-frequency analysis methods with high time-frequency resolution, but most of the time-frequency characterization methods cannot process complex multi-component signals without prior conditions.
The invention provides a novel linear frequency modulation transformation method, namely basis scaling adaptive frequency modulation transformation (SBCT), on the basis of linear frequency modulation transformation (CT), and the method adaptively matches multi-component signals by constructing a basis function. Compared with the existing method, the SBCT method can obtain the time-frequency characterization to realize higher time-frequency energy aggregation. The time-frequency representation results obtained by a Synchronous Extraction Transformation (SET) method after the energy of the time-frequency spectrum is redistributed are more aggregated, and the resolution is greatly improved. From the above analysis, we know that the SET and SBCT have respective advantages in the signal processing field, and thus perform a joint analysis discussion thereof.
The basic idea of the method for predicting the seismic reservoir by the basis-expansion frequency modulation synchronous extraction is to obtain a time-frequency representation result of a signal through SBCT, and then rearrange time-frequency energy to a real instantaneous frequency ridge of the signal so as to realize high aggregation of the time-frequency energy of the signal; the real instantaneous frequency of the signal is estimated on a time-frequency domain of basis-extension adaptive frequency conversion (SBCT), the time-frequency energy which is closely related to the time-frequency characteristic of the signal is extracted from an original time-frequency spectrum by using a frequency fixed point, and a plurality of fuzzy time-frequency energies are removed, so that the time-frequency energy distribution information of the signal is accurately carved, and the positive promotion effect is generated on the improvement of the time-frequency analysis accuracy.
Disclosure of Invention
Aiming at the defects in the prior art, the invention provides a prediction method for synchronously extracting a seismic reservoir by using a fundamental stretch frequency modulation. The method can improve the focusing performance of signal time-frequency representation energy distribution, improve the identification precision of the seismic reservoir and well reconstruct the signal.
In order to achieve the purpose of the invention, the invention adopts the technical scheme that: a method for predicting a base expansion frequency modulation synchronous extraction seismic reservoir stratum comprises the following steps:
s1, inputting an original two-dimensional seismic section signal S (x; t) to be analyzed;
s2, performing basis expansion self-adaptive frequency conversion on the two-dimensional seismic section signals S (x; t) input in the S1 to obtain a time frequency conversion result SBCT (x; f, t);
s3, estimating the instantaneous frequency of the x-th seismic signal at each time-frequency position according to the obtained time-frequency transformation result SBCT (x; f, t);
s4, according to the synchronous extraction principle, a synchronous extraction operator which takes the signal instantaneous frequency as the curve center is constructed on the time-frequency domain to extract the original time-frequency spectrum energy to obtain a new time-frequency coefficient, and fuzzy time-frequency energy is removed;
and S5, performing modulus extraction on the obtained time frequency coefficient to obtain a time frequency spectrum after base stretch frequency modulation synchronous extraction transformation.
S6, calculating the main frequency range of the original seismic signal by adopting Fourier transform, setting 85% of the maximum frequency value in the main frequency range as a high-frequency value, setting the minimum frequency value in the main frequency range as a low-frequency value, and then extracting a frequency spectrum corresponding to the high-frequency value from the time-frequency spectrum obtained in S5 to obtain a common frequency profile of the high-frequency value; and extracting the frequency spectrum corresponding to the low frequency value from the time frequency spectrum obtained in the step S5 to obtain a common frequency section of the low frequency value. The presence of a reservoir in the seismic profile is determined by comparing the attenuation of the seismic signal in two common frequency profiles.
Preferably, the time-frequency transformation result of the fundamental warping adaptive transformation of the original two-dimensional seismic section signal S (x; t) in step S2 is:
Figure BDA0003078080070000031
where t denotes time, f denotes the center of frequency, and h (x; τ) denotes a Gaussian function.
Figure BDA0003078080070000032
Representing a phase function, τ being a time variable, said phase function
Figure BDA0003078080070000033
Comprises the following steps:
Figure BDA0003078080070000034
where k is 1,2, …, n, n is a phase function
Figure BDA0003078080070000035
Order of (a)1,a2,…,an) The adaptive parameter of the phase function can be determined according to the kurtosis theory, namely:
Figure BDA0003078080070000036
where argmax (·) represents the set of parameters when the function takes the maximum value.
Preferably, the instantaneous frequency of the signal at each time-frequency position (t, f) is estimated from the resulting time-frequency transformation results SBCT (x; f, t)
Figure BDA0003078080070000037
Figure BDA0003078080070000038
Wherein,
Figure BDA0003078080070000039
representing the partial derivative of the window function h (x; t) over t,
Figure BDA00030780800700000310
expressed in a window function
Figure BDA00030780800700000311
(ii) the SBCT result of (B);
Figure BDA00030780800700000312
is expressed in function tkSBCT results under h (x; t).
Preferably, the specific method for acquiring the new time-frequency coefficient in step S4 is as follows: according to the synchronous extraction principle, a synchronous extraction operator SEO (x; f, T) taking the signal instantaneous frequency as the center of a curve is constructed on the time-frequency domain and used for extracting the original time-frequency spectrum energy to the curve to obtain a new time-frequency coefficient T (x; f, T), fuzzy time-frequency energy is removed, and the synchronous extraction operator SEO (x; f, T) meets the following requirements:
Figure BDA0003078080070000041
wherein δ (-) is a unit pulse function, and the new time-frequency coefficient T (x; f, T) is:
T(x；f,t)＝SBCT(x；f,t)SEO(x；f,t) (6)
preferably, the step S5 may adopt the following formula to perform inverse transformation on T (x; f, T) in (6) to reconstruct the seismic signal S (x; T);
Figure BDA0003078080070000042
wherein,
Figure BDA0003078080070000043
the method for predicting the seismic reservoir by synchronously extracting the base stretch frequency modulation provided by the invention has the beneficial effects that:
the invention provides a basic telescopic frequency modulation synchronous extraction algorithm by taking an SBCT algorithm as an entry point and considering the problems of unclear identification, poor aggregation and the like in the SET multi-component non-stable non-linear complex signal. Firstly, acquiring a time-frequency representation result of a signal through SBCT, and rearranging time-frequency energy to a real instantaneous frequency ridge line of the signal so as to realize high aggregation of the time-frequency energy of the signal; the real instantaneous frequency of the signal is estimated on the SBCT time-frequency domain, the time-frequency energy which is closely related to the time-frequency characteristics of the signal is extracted from the original time-frequency spectrum by using the frequency fixed point, and a plurality of fuzzy time-frequency energies are removed, so that the time-frequency energy distribution information of the signal is carved at high precision, and the positive promotion effect is generated on improving the time-frequency analysis precision. The method has good effects on reconstruction capability and calculation efficiency, and can remarkably improve time-frequency energy focusing and improve the identification precision of the seismic reservoir.
The idea of the invention is as follows:
firstly, inputting an original two-dimensional seismic section signal s (x; t) to be analyzed;
secondly, performing basis-extension adaptive frequency conversion on the input signal s (x; t), wherein the calculation method comprises the following steps of obtaining a time-frequency conversion result SBCT (x; f, t):
Figure BDA0003078080070000051
wherein t represents time, f represents frequency center, h (x; tau) epsilon L2(R) is expressed as a Gaussian function,
Figure BDA0003078080070000052
representing a phase function, tau being a time variable;
thirdly, the instantaneous frequency of the signal at each time-frequency position (t, f) is estimated according to the obtained time-frequency transformation result SBCT (x; f, t)
Figure BDA0003078080070000053
Fourthly, according to the synchronous extraction principle, a synchronous extraction operator SEO (x; f, T) which takes a signal instantaneous frequency curve as the center is constructed on the time-frequency domain and is used for extracting the original time-frequency spectrum energy to the curve to obtain a new time-frequency coefficient T (x; f, T) and removing the fuzzy time-frequency energy;
fifthly, modulus is taken for the time frequency coefficient obtained in the fourth step, and a time frequency spectrum after base stretch frequency modulation synchronous extraction transformation is obtained.
Sixthly, calculating the main frequency range of the original seismic signal by adopting Fourier transform, setting 85% of the maximum frequency value in the main frequency range as a high-frequency value, setting the minimum frequency value in the main frequency range as a low-frequency value, and then extracting a frequency spectrum corresponding to the high-frequency value from the time-frequency spectrum obtained by S5 to obtain a common frequency profile of the high-frequency value; and extracting the frequency spectrum corresponding to the low frequency value from the time frequency spectrum obtained in the step S5 to obtain a common frequency section of the low frequency value. The presence of a reservoir in the seismic profile is determined by comparing the attenuation of the seismic signal in two common frequency profiles.
The working principle of the invention is as follows: acquiring an original two-dimensional seismic profile signal s (x; t) to be analyzed; decomposing the original two-dimensional seismic section signal s (x; t) by using basis warping adaptive frequency modulation transform (SBCT) to obtain a time frequency spectrum SBCT (x; f, t) of the basis warping adaptive frequency modulation transform of the original two-dimensional seismic section signal s (x; t); estimating the instantaneous frequency of the signal at each time-frequency position according to the time-frequency transformation result SBCT (x; f, t); according to the synchronous extraction principle, a synchronous extraction operator which takes the signal instantaneous frequency as the curve center is constructed on the time-frequency domain and is used for extracting the original time-frequency spectrum energy to obtain a new time-frequency coefficient, and fuzzy time-frequency energy is removed; and taking a module of the time frequency coefficient to obtain a time frequency spectrum after the base stretch frequency modulation synchronous extraction transformation. The method can obviously improve the time-frequency energy focusing property and improve the identification precision of the seismic reservoir.
The invention provides a basic telescopic frequency modulation synchronous extraction algorithm by taking an SBCT algorithm as an entry point and combining the problems of unclear identification, poor aggregation and the like in the SET multi-component non-stable non-linear complex signal. Firstly, acquiring a time-frequency representation result of a signal through SBCT, and rearranging time-frequency energy to a real instantaneous frequency ridge line of the signal so as to realize high aggregation of the time-frequency energy of the signal; the real instantaneous frequency of the signal is estimated on the SBCT time-frequency domain, the time-frequency energy which is closely related to the time-frequency characteristic of the signal is extracted from the original time-frequency spectrum by using the frequency fixed point, a plurality of fuzzy time-frequency energies are removed, the time-frequency energy distribution information of the signal is carved with high precision, and the positive promotion effect is generated on improving the precision of time-frequency analysis. The method has good effects on reconstruction capability and calculation efficiency, and can remarkably improve time-frequency energy focusing and improve the identification precision of the seismic reservoir.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a cross-sectional view of the gas field in the Sichuan basin in the southwest and the middle
FIG. 3 is a low frequency value common frequency section of Zhongjiang gas field processed by the method of the present invention
FIG. 4 is a high-frequency value common-frequency profile of Zhongjiang gas field processed by the method of the present invention
Detailed Description
The invention will be further explained with reference to the drawings.
Example 1: referring to fig. 1, a method for synchronously extracting seismic reservoirs by using basis-stretch frequency modulation comprises the following steps:
s1, inputting an original two-dimensional seismic section signal S (x; t) to be analyzed;
s2, decomposing the original two-dimensional seismic section signal S (x; t) by using the basis stretching frequency modulation adaptive transformation to obtain a basis stretching frequency modulation adaptive transformation result SBCT (x; f, t) of the original two-dimensional seismic section signal S (x; t);
s3, estimating the instantaneous frequency of the signal at each time frequency position according to the time frequency transformation result SBCT (x; f, t);
s4, according to the synchronous extraction principle, a synchronous extraction operator which takes the signal instantaneous frequency as the curve center is constructed on the time-frequency domain for extracting new time-frequency coefficients obtained from the original time-frequency spectrum energy and eliminating fuzzy time-frequency energy;
s5, taking a module of the time frequency coefficient to obtain a time frequency spectrum after base stretch frequency modulation synchronous extraction transformation.
S6, calculating the main frequency range of the original seismic signal by adopting Fourier transform, setting 85% of the maximum frequency value in the main frequency range as a high-frequency value, setting the minimum frequency value in the main frequency range as a low-frequency value, and then extracting a frequency spectrum corresponding to the high-frequency value from the time-frequency spectrum obtained in S5 to obtain a common frequency profile of the high-frequency value; and extracting the frequency spectrum corresponding to the low frequency value from the time frequency spectrum obtained in the step S5 to obtain a common frequency section of the low frequency value. The presence of a reservoir in the seismic profile is determined by comparing the attenuation of the seismic signal in two common frequency profiles.
Preferably, the time-frequency transformation result of the fundamental warping adaptive transformation of the original two-dimensional seismic section signal S (x; t) in step S2 is:
Figure BDA0003078080070000071
wherein t represents time, f represents a frequency center, h (x; τ) represents a Gaussian function,
Figure BDA0003078080070000072
representing a phase function, tau being a time variable;
the phase function
Figure BDA0003078080070000073
Comprises the following steps:
Figure BDA0003078080070000074
where k is 1,2, …, n, n is a phase function
Figure BDA0003078080070000081
Order of (a)1,a2,…,an) The adaptive parameter of the phase function can be determined according to the kurtosis theory, namely:
Figure BDA0003078080070000082
where argmax (·) represents the set of parameters when the function takes the maximum value.
Preferably, in step S3, the instantaneous frequency of the signal at each time-frequency position (t, f) is estimated according to the time-frequency spectrum phase information as:
Figure BDA0003078080070000083
wherein,
Figure BDA0003078080070000084
representing the partial derivative of the window function h (x; t) over t,
Figure BDA0003078080070000085
expressed in a window function
Figure BDA0003078080070000086
(ii) the SBCT result of (B);
Figure BDA0003078080070000087
is expressed in function tkSBCT results under h (x; t).
Preferably, the specific method for acquiring the new time-frequency coefficient in step S4 is as follows: according to the synchronous extraction principle, a synchronous extraction operator SEO (x; f, T) taking the signal instantaneous frequency as the center of a curve is constructed on the time-frequency domain and used for extracting the original time-frequency spectrum energy to the curve to obtain a new time-frequency coefficient T (x; f, T), fuzzy time-frequency energy is removed, and the synchronous extraction operator SEO (x; f, T) meets the following requirements:
Figure BDA0003078080070000088
wherein δ (-) is a unit pulse function, and the new time-frequency coefficient T (x; f, T) is:
T(x；f,t)＝SBCT(x；f,t)SEO(x；f,t) (6)
preferably, the step S5 is to apply the following formula to perform inverse transformation on T (x; f, T) in (6) to reconstruct the seismic signal S (x; T)
Figure BDA0003078080070000089
Wherein,
Figure BDA0003078080070000091
referring to fig. 1 to 4, we take a seismic section as an example, and the two-dimensional seismic section is shown in fig. 2. Fig. 3 to 4 are a low-frequency value common frequency section and a high-frequency value common frequency section extracted from a time spectrum obtained by basis expansion adaptive frequency conversion of a seismic section, respectively. In the figure, the abscissa represents the number of seismic traces, the ordinate represents time, and the right-hand color scale represents energy values. The embodiment proves that the time-frequency resolution of the result graph obtained after the processing of the method is higher, the energy is more concentrated, and the identification precision of the seismic reservoir is improved.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents and improvements made within the spirit and principle of the present invention are intended to be included within the scope of the present invention.

Claims (5)

1. A method for predicting a seismic reservoir by synchronously extracting a base expansion frequency modulation is characterized by comprising the following steps:
s1, inputting an original two-dimensional seismic section signal S (x; t) to be analyzed; wherein x represents a seismic trace and t represents time;
s2, performing basis expansion self-adaptive frequency conversion on the two-dimensional seismic section signals S (x; t) input in the S1 to obtain a time frequency conversion result SBCT (x; f, t), wherein the calculation method comprises the following steps:
Figure FDA0003078080060000011
wherein t represents a time center, f represents a frequency center, h (x; τ) represents a Gaussian function,
Figure FDA0003078080060000012
representing a phase function, tau being a time variable;
s3, estimating the instantaneous frequency of the x-th seismic signal at each time-frequency position (t, f) according to the time-frequency transformation result SBCT (x; f, t) obtained in S2
Figure FDA0003078080060000013
S4, according to the synchronous extraction principle, constructing a synchronous extraction operator SEO (x; f, T) taking a signal instantaneous frequency curve as a center on the time-frequency domain for extracting the original time frequency spectrum energy, obtaining a new time-frequency coefficient T (x; f, T), and removing fuzzy time-frequency energy;
s5, performing modulus extraction on the time-frequency coefficient obtained in the S4 to obtain a time-frequency spectrum after base stretch frequency modulation synchronous extraction transformation;
s6, calculating the main frequency range of the original seismic signal S (x; t) by Fourier transform, setting 85% of the maximum frequency value in the main frequency range as a high-frequency value, setting the minimum frequency value in the main frequency range as a low-frequency value, and then extracting the frequency spectrum corresponding to the high-frequency value from the time-frequency spectrum obtained in S5 to obtain a common frequency profile of the high-frequency value; and extracting the frequency spectrum corresponding to the low frequency value from the time frequency spectrum obtained in the step S5 to obtain a common frequency section of the low frequency value. The presence of a reservoir in the seismic profile is determined by comparing the attenuation of the seismic signal in two common frequency profiles.
2. The method for synchronously extracting the seismic reservoir prediction according to the fundamental tone of claim 1, wherein the time-frequency transformation result of the fundamental tone adaptive transformation of the two-dimensional seismic profile signal S (x; t) in the step S2 is as follows:
Figure FDA0003078080060000021
wherein t represents a time center, f represents a frequency center, h (x; τ) represents a Gaussian function,
Figure FDA0003078080060000022
representing a phase function, τ being a time variable, said phase function
Figure FDA0003078080060000023
Comprises the following steps:
Figure FDA0003078080060000024
where k is 1,2, …, n, n is a phase function
Figure FDA0003078080060000025
Order of (a)1,a2,…,an) The adaptive parameter of the phase function can be determined according to the kurtosis theory, namely:
Figure FDA0003078080060000026
where argmax (·) represents the set of parameters when the function takes the maximum value.
3. The PSFM seismic reservoir prediction method as claimed in claim 1, wherein the calculation of the instantaneous frequency in step S3 comprises estimating the instantaneous frequency of the signal at each time-frequency position (t, f) from the SBCT (x; f, t) result obtained in step S2
Figure FDA0003078080060000027
Figure FDA0003078080060000028
Wherein,
Figure FDA0003078080060000029
representing the partial derivative of the window function h (x; t) over time t,
Figure FDA00030780800600000210
expressed in a window function
Figure FDA00030780800600000211
(ii) the SBCT result of (B);
Figure FDA00030780800600000212
is expressed in function tkSBCT results under h (x; t).
4. The method for synchronously extracting seismic reservoirs according to the fundamental stretch frequency modulation of claim 1, wherein the specific method for obtaining the new time-frequency coefficient T (x; f, T) in the step S4 is as follows: according to the synchronous extraction principle, a synchronous extraction operator SEO (x; f, T) taking the signal instantaneous frequency as the curve center is constructed on the time-frequency domain to extract the original time-frequency spectrum energy to obtain a new time-frequency coefficient T (x; f, T), the fuzzy time-frequency energy is removed, and the synchronous extraction operator SEO (x; f, T) meets the following requirements:
Figure FDA0003078080060000031
wherein δ (-) is a unit pulse function, and the new time-frequency coefficient T (x; f, T) is:
T(x；f,t)＝SBCT(x；f,t)SEO(x；f,t) (7) 。
5. the method for predicting a seismic reservoir according to claim 1, wherein said step S5 is implemented by transforming T (x; f, T) in (7) with the following formula to reconstruct a seismic signal S (x; T);
Figure FDA0003078080060000032
wherein,
Figure FDA0003078080060000033
the invention is suitable for reconstructing seismic signals.
CN202110558646.1A 2021-05-21 2021-05-21 Seismic reservoir prediction method by synchronous extraction of base stretch frequency modulation Pending CN113296155A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110558646.1A CN113296155A (en) 2021-05-21 2021-05-21 Seismic reservoir prediction method by synchronous extraction of base stretch frequency modulation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110558646.1A CN113296155A (en) 2021-05-21 2021-05-21 Seismic reservoir prediction method by synchronous extraction of base stretch frequency modulation

Publications (1)

Publication Number Publication Date
CN113296155A true CN113296155A (en) 2021-08-24

Family

ID=77323730

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110558646.1A Pending CN113296155A (en) 2021-05-21 2021-05-21 Seismic reservoir prediction method by synchronous extraction of base stretch frequency modulation

Country Status (1)

Country Link
CN (1) CN113296155A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114280671A (en) * 2021-09-17 2022-04-05 成都理工大学 Generalized W-transform seismic gas-bearing reservoir characterization method
CN114355442A (en) * 2022-01-12 2022-04-15 成都理工大学 Three-parameter W-transform seismic reservoir identification time-frequency analysis method

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114280671A (en) * 2021-09-17 2022-04-05 成都理工大学 Generalized W-transform seismic gas-bearing reservoir characterization method
CN114355442A (en) * 2022-01-12 2022-04-15 成都理工大学 Three-parameter W-transform seismic reservoir identification time-frequency analysis method

Similar Documents

Publication Publication Date Title
CN106597532B (en) Pre-stack seismic data frequency band expanding method combining well data and horizon data
CN113296154A (en) Method for predicting base expansion frequency modulation synchronous extrusion seismic reservoir
CN106500735B (en) A kind of compressed sensing based FBG signal adaptive restorative procedure
CN109031422A (en) A kind of seismic signal noise suppressing method based on CEEMDAN and Savitzky-Golay filtering
CN113296155A (en) Seismic reservoir prediction method by synchronous extraction of base stretch frequency modulation
CN113378661A (en) Direct current electric energy signal denoising method based on improved wavelet threshold and related detection
CN105353408B (en) A kind of Wigner higher-order spectrum seismic signal spectral factorization methods based on match tracing
CN104182946A (en) Spectral denoising method based on index soft and hard averaging threshold function
CN113851144A (en) Voice signal denoising method based on improved variational modal decomposition and principal component analysis
CN113537102B (en) Feature extraction method of microseismic signals
CN113777650B (en) Sparse time-frequency spectrum decomposition method, device and equipment based on mixed norm and wavelet transform and storage medium
CN103345920B (en) Self-adaptation interpolation weighted spectrum model voice conversion and reconstructing method based on Mel-KSVD sparse representation
Li Speech denoising based on improved discrete wavelet packet decomposition
CN105187341A (en) Stationary wavelet transform denoising method based on cross validation
CN112530449B (en) Speech enhancement method based on bionic wavelet transform
CN111505709B (en) Attenuation qualitative analysis method based on sparse spectral decomposition
CN104778342B (en) A kind of heart sound feature extracting method based on wavelet singular entropy
CN106443771B (en) Improve the method and its velocity inversion method of converted wave seismic data resolution
CN106340304B (en) A kind of online sound enhancement method under the environment suitable for nonstationary noise
CN103915102A (en) Method for noise abatement of LFM underwater sound multi-path signals
CN113095113A (en) Wavelet line spectrum feature extraction method and system for underwater target identification
CN106997766B (en) Homomorphic filtering speech enhancement method based on broadband noise
CN116383605A (en) Vehicle vibration signal denoising method based on wavelet decomposition
CN107025911B (en) Fundamental frequency detection method based on particle group optimizing
CN113066483B (en) Sparse continuous constraint-based method for generating countermeasure network voice enhancement

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