CN111323816A - Instantaneous phase gradient attribute extraction method based on ocean broadband seismic data waveform - Google Patents
Instantaneous phase gradient attribute extraction method based on ocean broadband seismic data waveform Download PDFInfo
- Publication number
- CN111323816A CN111323816A CN202010199904.7A CN202010199904A CN111323816A CN 111323816 A CN111323816 A CN 111323816A CN 202010199904 A CN202010199904 A CN 202010199904A CN 111323816 A CN111323816 A CN 111323816A
- Authority
- CN
- China
- Prior art keywords
- instantaneous phase
- seismic data
- phase gradient
- attribute
- broadband
- 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
Links
- 238000000605 extraction Methods 0.000 title claims abstract description 6
- 238000000034 method Methods 0.000 claims abstract description 34
- 238000012545 processing Methods 0.000 claims abstract description 9
- 238000004804 winding Methods 0.000 claims abstract description 9
- 230000009466 transformation Effects 0.000 claims abstract description 5
- 230000001629 suppression Effects 0.000 claims abstract description 4
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 230000000717 retained effect Effects 0.000 claims description 2
- 230000008569 process Effects 0.000 abstract description 3
- 238000004364 calculation method Methods 0.000 description 10
- 238000007796 conventional method Methods 0.000 description 7
- 230000001052 transient effect Effects 0.000 description 7
- 230000008901 benefit Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000011161 development Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000008021 deposition Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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/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
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
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 an instantaneous phase gradient attribute extraction method based on ocean broadband seismic data waveform, which comprises the following steps: acquiring original seismic data through an inclined cable, and performing ghost wave suppression processing to obtain broadband seismic data; carrying out Hilbert transformation on the broadband seismic data, and calculating to obtain an instantaneous phase attribute; the winding phenomenon of the instantaneous phase is kept, and the first derivative related to time is obtained for the instantaneous phase to obtain the instantaneous phase gradient; and identifying the sequence interface by using the instantaneous phase gradient. The invention can be widely applied to marine broadband seismic data, effectively improves the resolution capability of the underground sequence interface and reduces the multi-solution of interpretation. The invention can be widely applied to the marine oil and gas exploration process.
Description
Technical Field
The invention relates to the field of petroleum exploration and development, in particular to an instantaneous phase gradient attribute extraction method based on ocean broadband seismic data waveform characteristics.
Background
In recent years, the technology of collecting and matching processing seismic data of ocean deep water inclined cables is greatly developed, and the frequency band range of the seismic data is greatly widened. The broadband seismic data can improve the imaging precision of an underground complex structure and improve the accuracy of a seismic inversion result, and has remarkable advantages in the aspects of sediment sequence explanation and reservoir parameter prediction. In the practical application process, the potential of the broadband seismic data is far from being thoroughly mined, mainly because the interpretation and attribute calculation method of the broadband seismic data still continues the thought of the conventional narrowband seismic data. In fact, compared with the conventional seismic data, the wavelet morphology of the broadband seismic data has obvious change, but an attribute calculation method specially aiming at the waveform characteristics of the broadband seismic data is lacked.
The instantaneous phase attribute of the conventional seismic signal can effectively reflect internal deposition phenomena such as faults, sand pinch-out, river channels, deposition fans and the like, and is widely applied to exploration and development. There are many methods for calculating the instantaneous phase and instantaneous frequency attributes of seismic signals, such as Hilbert transform, fourier transform, Hilbert-Huang transform, etc. However, the methods all use the arctangent function when calculating the instantaneous phase, and the computer has the phase winding phenomenon when the phase changes from 0-2 pi. Ocean seismic data are generally influenced by ghost waves, and the phase discontinuity characteristics can generate obvious multi-solution due to the side lobe effect of seismic wavelets, so that unwrapping processing is required when the instantaneous frequency is obtained by utilizing the instantaneous phase attribute of conventional data, and the corresponding relation between the phase discontinuity characteristics and the sequence interface is ignored.
Disclosure of Invention
In view of the above problems, an object of the present invention is to provide an attribute calculation method for wavelet phase discontinuity characteristics of broadband seismic data based on marine broadband seismic data waveform characteristics, where the method can be used for accurately identifying a sedimentary sequence interface, and can effectively improve the identification capability of a thin-layer sequence interface and reduce the multi-solution of interpretation results.
In order to achieve the purpose, the invention adopts the following technical scheme:
the invention provides a method for extracting an instantaneous phase gradient attribute based on a marine broadband seismic data waveform, which comprises the following steps of:
1) acquiring broadband seismic data in a work area;
2) carrying out Hilbert transformation on the broadband seismic data in the step 1) to obtain an instantaneous phase attribute theta (t);
3) the winding phenomenon of the instantaneous phase is kept, and the first derivative of the instantaneous phase attribute theta (t) with respect to time is solved to obtain an instantaneous phase gradient omega (t);
4) and (3) utilizing the instantaneous phase gradient omega (t) to identify the sequence interface.
In the step 1), original seismic data of a work area are acquired through the inclined cables, ghost wave suppression processing is carried out, and broadband seismic data are obtained.
Wherein, the method for calculating the instantaneous phase attribute θ (t) in the step 2) is as follows: the signal is complex through Hilbert conversion, and after the signal x (t) is subjected to Hilbert conversion, the signal Z (t) is analyzed to be:
Z(t)=x(t)+iY(t)=A(t)eiθ(t)(1)
the instantaneous phase property θ (t) is found by the imaginary and real parts of the transformed signal, expressed as:
where a represents the amplitude of the signal, i represents the imaginary unit, t represents time, x (t) represents the input signal, y (t) represents the transformed imaginary component, and z (t) represents the analytic signal.
Wherein, the expression of the instantaneous phase gradient ω (t) in the step 3) is:
where ω (t) represents the instantaneous phase gradient, θ (t) represents the instantaneous phase attribute, and t represents time.
When the instantaneous phase gradient ω (t) is obtained in the step 3), the winding information of the phase is retained, and the first derivative of the instantaneous phase is directly obtained.
Wherein, the identification method in the step 4) is as follows: judging the positive polarity and the negative polarity of the seismic data, and if the seismic data is positive polarity, judging that the peak value position of the instantaneous phase gradient attribute corresponds to a negative reflection coefficient, and judging that the peak value interface of the instantaneous phase gradient attribute is an interface with wave impedance from large to small; and if the seismic data is negative, judging that the instantaneous phase gradient attribute peak interface is an interface with wave impedance from small to large if the instantaneous phase gradient attribute peak position corresponds to a positive reflection coefficient.
The identification method in the step 4) judges the positive polarity and the negative polarity of the seismic data on the premise of stable phase of the wideband seismic data of the work area to be identified.
Due to the adoption of the technical scheme, the invention has the following advantages: 1. in the traditional method, although the first derivative of the instantaneous phase is obtained when the instantaneous frequency is calculated, the instantaneous phase is subjected to unwrapping treatment, and the corresponding relation between the discontinuous characteristic of the instantaneous phase and the stratum sequence interface is ignored; the instantaneous phase gradient attribute obtained by the method is changed along with time, the corresponding relation between the instantaneous phase discontinuity characteristic and the sequence interface is reserved, and the sequence interface position can be accurately identified. 2. The phase characteristics of the conventional marine seismic data wavelets are unstable, so that real zero-phase seismic data are often difficult to obtain, and the accurate corresponding relation between the position of phase discontinuity change and the sequence interface cannot be ensured; the invention is based on the broadband seismic data, the wavelet phase characteristics are stable, and the one-to-one correspondence between the instantaneous phase gradient attribute peak position and the sequence interface can be ensured. 3. The conventional marine seismic data are generally influenced by ghost waves, the side lobe effect of the wavelets is serious, and the multi-resolution is strong if the instantaneous phase gradient attribute is calculated based on the conventional seismic data; the broadband seismic data not only broadens low-frequency information, but also broadens high-frequency information, wherein the abundant low-frequency information can effectively improve the signal-to-noise ratio of the seismic data, suppress wavelet side lobes and reduce the multi-solution property of the instantaneous phase gradient property; the high-frequency information in the broadband data can greatly improve the resolution of the instantaneous phase gradient attribute, and greatly improve the description capacity of the thin reservoir and the boundary thereof. 4. Because the method for calculating the instantaneous frequency attribute in the prior art is carried out aiming at the conventional seismic data, the advantages of the broadband seismic data cannot be fully exerted if the method is directly used for the broadband seismic data; the instantaneous phase gradient attribute with the reserved phase winding characteristic is provided for the waveform characteristics of the broadband seismic data, the peak position of the instantaneous phase gradient attribute has a stable corresponding relation with the stratum sequence interface, and the calculation precision is high. The invention can be widely applied to the marine broadband seismic exploration and development process.
Drawings
FIG. 1 is a graph of the calculated instantaneous phase gradient property of a positive polarity Riker wavelet with a dominant frequency of 35 Hz;
FIG. 2 is a graph of calculated instantaneous phase gradient properties of a negative polarity Riker wavelet with a dominant frequency of 35 Hz;
FIG. 3 is a graph of the instantaneous phase gradient property calculated based on positive polarity broadband seismic wavelets in accordance with the present invention;
FIG. 4 is a graph of the instantaneous phase gradient property calculated based on negative polarity broadband seismic wavelets in accordance with the present invention;
FIG. 5 is a schematic diagram illustrating an example of a positive polarity ideal wavelet with a frequency of 0-100 Hz;
FIG. 6 is a schematic diagram illustrating an example of a negative polarity ideal wavelet with a frequency of 0-100 Hz;
FIG. 7 is a diagram illustrating an example of the computation of different wavelet fit seismic signals;
FIG. 8 is a plot of the results of a seismic profile and a transient attribute profile taken somewhere using the conventional method and the method of the present invention, respectively.
Detailed Description
The invention is described in detail below with reference to the figures and examples.
The invention provides an attribute extraction method based on ocean broadband seismic data waveform characteristics, which comprises the following steps:
1) acquiring original seismic data through an inclined cable, and performing ghost wave suppression processing to obtain broadband seismic data;
2) carrying out Hilbert transformation on the broadband seismic data, and calculating to obtain an instantaneous phase attribute theta (t);
3) the winding phenomenon of the instantaneous phase is kept, and the first derivative of the instantaneous phase theta (t) with respect to time is obtained to obtain an instantaneous phase gradient omega (t);
in the above steps 2) and 3), the method for calculating the instantaneous phase attribute θ (t) and the instantaneous phase gradient ω (t) is as follows: the signal is complex through Hilbert transform, and after the signal x (t) is Hilbert transformed, its analytic signal z (t) can be expressed as:
Z(t)=x(t)+iY(t)=A(t)eiθ(t)(1)
the instantaneous phase property θ (t) is found by the imaginary and real parts of the transformed signal, expressed as:
when the instantaneous phase is calculated according to the above formula (2), the arctan function is used,
different from the traditional method for solving the instantaneous frequency attribute, the invention reserves the winding information of the phase position, and directly solves the first derivative of the instantaneous phase position to obtain the instantaneous phase gradient expression:
4) and (3) identifying the sequence interface by using the instantaneous phase gradient attribute profile omega (t), wherein the identification method comprises the following steps: when the phase of the wideband seismic data of the work area to be identified is stable, if the seismic data is positive, the peak position of the instantaneous phase gradient attribute is corresponding to the negative reflection coefficient, and the interface of the peak value of the instantaneous phase gradient attribute is judged to be the interface with the wave impedance from large to small; when the phase of the wideband seismic data of the work area to be identified is stable, if the seismic data is negative, the peak position of the instantaneous phase gradient attribute is corresponding to the positive reflection coefficient, and the interface of the peak value of the instantaneous phase gradient attribute is judged to be the interface from small to large in wave impedance.
Example 1
Fig. 1 and fig. 2 are results of a method without using broadband seismic data, that is, an effect analysis of wavelets with different polarities when the broadband seismic data is not used. FIG. 1 is a graph of calculated instantaneous phase gradient properties of a positive polarity Riker wavelet with a dominant frequency of 35Hz, wherein FIG. 1(a) is a zero phase positive polarity Ricker wavelet with a dominant peak and two significant side lobes; FIG. 1(b) shows the instantaneous phase obtained by Hilbert transform of the wavelet, where a winding phenomenon occurs from pi to-pi in the phase at the position of the corresponding wavelet side lobe, and two discontinuous change points occur; fig. 1(c) shows that the property of the instantaneous phase gradient after the phase wrapping phenomenon is preserved to obtain the first derivative of the instantaneous phase, and two false discontinuous interfaces appear at the positions corresponding to the wavelet sidelobe, which indicates that the wavelet sidelobe effect of the conventional seismic data can generate multi-solution to the calculation result. FIG. 2 is a graph of calculated instantaneous phase gradient property of a negative polarity Riker wavelet with a dominant frequency of 35Hz, wherein FIG. 2(a) is a zero phase negative polarity Ricker wavelet with a dominant peak and two significant side lobes; FIG. 2(b) shows that the phase wrapping phenomenon does not occur at the side lobe position of the corresponding wavelet, and the phase wrapping occurs only at the main peak position of the wavelet, for the instantaneous phase obtained by Hilbert transform of the wavelet; fig. 2(c) shows the transient phase gradient property after the first derivative of the transient phase is obtained by preserving the phase wrap phenomenon, and a discontinuous interface is generated only at the position of the main peak, which illustrates that the transient phase gradient property generates a discontinuous interface only at the position where the phase wrap occurs.
Both fig. 3 and 4 are the results of using the acquisition method of the present invention in conjunction with broadband seismic data. FIG. 3 is an instantaneous phase gradient property calculated based on a positive polarity broadband seismic wavelet of the present invention, wherein FIG. 3(a) is a positive polarity broadband seismic wavelet, the frequency range of the wavelet is 3-100 Hz, and the energy of the wavelet sidelobe has been suppressed; FIG. 3(b) shows that the phase wrapping phenomenon does not occur at the positions of the corresponding wavelet side lobes and the wavelet main peak for the instantaneous phase obtained by Hilbert transform of the wavelet; FIG. 3(c) is the instantaneous phase gradient property after the first derivative is obtained for the instantaneous phase, and since the wavelet phase has no phase wrapping phenomenon, no discontinuous interface is generated in the calculated instantaneous phase gradient property, which illustrates that the instantaneous phase gradient property calculated based on the broadband seismic data does not generate multi-solution compared with FIG. 1 (c);
FIG. 4 is an instantaneous phase gradient attribute calculated based on negative polarity broadband seismic wavelets, wherein FIG. 4(a) is a negative polarity broadband seismic wavelet, the frequency range of the wavelet is 3-100 Hz, and the energy of the wavelet sidelobe has been suppressed; FIG. 4(b) shows the instantaneous phase obtained by Hilbert transform of the wavelet, where a phase wrap-around phenomenon occurs at the position of the main peak of the corresponding wavelet; FIG. 4(c) is the property of the instantaneous phase gradient after the first derivative is taken over the instantaneous phase, resulting in a discontinuous interface only at the location of the corresponding wavelet main peak; this is more accurate than the generation of a continuous interface, which indicates that the acquisition method of the present invention using the combination of broadband seismic data is more accurate.
When the ideal wavelet algorithm with frequency of 0-100 Hz is used, the result is shown in FIG. 5 and FIG. 6, wherein FIG. 5 is the positive polarity ideal wavelet and FIG. 6 is the negative polarity ideal wavelet algorithm. Because the low frequency of the ideal wavelet reaches 0Hz, compared with FIGS. 3-4, the sidelobe of the wavelet is completely suppressed, for the instantaneous phase obtained by Hilbert transformation of the two wavelets, the positive polarity wavelet shown in FIG. 5(b) does not generate the phase wrap-around phenomenon, and no discontinuous interface is generated in the instantaneous phase gradient attribute shown in FIG. 5 (c); phase wrapping occurs at the main peak position of the wavelet of negative polarity shown in FIG. 6(b), and accordingly a discontinuous interface is generated in the transient phase gradient profile shown in FIG. 6 (c).
The method for identifying a gas reservoir using wavelet phase characteristics of the present invention is further described below by way of specific embodiments.
Example 2:
an example of the computation of the fitting of different wavelets to the seismic signal, as shown in FIG. 7, is presented to illustrate the advantage of computing instantaneous phase gradient properties based on broad frequency seismic data. FIG. 7(a) shows the formation reflection coefficients, which correspond to three reflection coefficients at time 250ms, 350ms and 450ms, respectively, and the reflection coefficient values are: -0.5, 0.2, -0.4; the positive polarity wavelets and the reflection coefficients of FIG. 1(a), FIG. 3(a) and FIG. 5(a) are convolved to obtain corresponding seismic trace data and calculate the corresponding instantaneous phase gradient attributes. FIG. 7(b) is a reflected seismic signal synthesized by 35Hz positive polarity rake wavelets, and FIG. 7(c) is an instantaneous frequency attribute calculated by a conventional method, because the conventional method performs unwrapping processing on phase discontinuity characteristics, no discontinuous interface appears in the calculation result, and the resolution is very low; FIG. 7(d) is the result of calculating the instantaneous phase gradient property for the reflection signal of a conventional wavelet fit, and it can be seen that a false discontinuous interface appears at the wavelet side lobe position of the regular reflection coefficient, illustrating that the method produces ambiguity for conventional seismic data; FIG. 7(e) is a 3-100 Hz positive polarity broadband wavelet fitting seismic signal, FIG. 7(f) is a corresponding instantaneous phase gradient attribute result, a discontinuous interface is not generated at a position corresponding to a positive reflection coefficient, an accurate discontinuous interface is only generated at a position corresponding to a negative reflection coefficient, although the discontinuity at the position of the positive reflection coefficient is eliminated, the result is accurate and has no ambiguity problem, and when the thickness of a stratum is thin, the resolution is high; FIGS. 7(g) and 7(h) show the positive polarity ideal wavelet fitting seismic signals at 0-100 Hz and the transient phase gradient properties, which are consistent with the broad frequency wavelets shown in FIGS. 7(e) and 7 (f).
Example 3
The conventional method and the method of the present invention are used to map the seismic profile and the transient attribute profile of a certain location, and the results are shown in fig. 8(a) to 8 (f). Wherein FIG. 8(a) is a seismic section of a conventional horizontal streamer acquisition, after conventional processing; FIG. 8(b) is a cross-sectional view of the instantaneous frequency property calculated by the conventional method on conventional data, which shows that the resolution is very low and the development characteristics of the water channel are difficult to see; FIG. 8(c) is a cross-sectional view of the instantaneous phase gradient property calculated based on conventional data by using the method of the present invention, because ghost wave interference exists in conventional seismic data, the position generated by a discontinuous interface does not have a good correspondence with a stratigraphic sequence interface, and obvious ambiguity exists in the calculation result; FIG. 8(d) is a seismic profile after streamer acquisition and broadband processing; FIG. 8(e) is a cross-sectional view of the instantaneous frequency profile calculated by conventional methods for wideband data, which, although having improved accuracy compared to FIG. 8(b), is still difficult to identify the location of the sequence interface; fig. 8(f) is an instantaneous phase gradient attribute profile calculated for the wideband seismic data according to the present invention, and since ghost waves are removed, the ambiguity problem in fig. 8(c) is solved well, and meanwhile, the energy of the high frequency part of the seismic data is enhanced, the accuracy of the calculation result is greatly improved, and there is a good correspondence relationship with the sequence interface position.
In summary, the present invention is obviously different from the prior art, and mainly includes the following points: 1) the traditional instantaneous frequency attribute calculation method is carried out based on the conventional horizontal streamer acquisition seismic data, the frequency band of the seismic data is narrow, the wavelet sidelobe effect is serious, in order to reduce the multi-solution of the calculation result, the corresponding relation between instantaneous phase discontinuity and a sequence interface is ignored, the phase discontinuity characteristic is used as an abnormal value to be subjected to unwrapping treatment, the resolution ratio of the finally calculated instantaneous frequency attribute is low, the conventional method is directly used for broadband seismic data, and the advantages of the broadband seismic data cannot be fully exerted; 2) the wavelet phase characteristics of the broadband seismic data are stable based on the seismic data acquired by the inclined cable broadband and processed by the broadband, and the instantaneous phase gradient attribute peak positions can be ensured to be in one-to-one correspondence with sequence interfaces; the abundant low-frequency information can effectively improve the signal-to-noise ratio of seismic data, suppress wavelet side lobes and reduce the multi-solution property of the instantaneous phase gradient; the high frequency information in the broadband data can greatly improve the resolution of the instantaneous phase gradient attribute, and greatly improve the description capability of the thin reservoir and the boundary thereof (as shown in fig. 8 (f)).
The above embodiments are only for illustrating the present invention, and the steps may be changed, and on the basis of the technical solution of the present invention, the modification and equivalent changes of the individual steps according to the principle of the present invention should not be excluded from the protection scope of the present invention.
Claims (7)
1. An instantaneous phase gradient attribute extraction method based on ocean broadband seismic data waveforms is characterized by comprising the following steps:
1) acquiring broadband seismic data in a work area;
2) carrying out Hilbert transformation on the broadband seismic data in the step 1) to obtain an instantaneous phase attribute theta (t);
3) the winding phenomenon of the instantaneous phase is kept, and the first derivative of the instantaneous phase theta (t) with respect to time is obtained to obtain an instantaneous phase gradient omega (t);
4) and (3) utilizing the instantaneous phase gradient omega (t) to identify the sequence interface.
2. The method for extracting attributes based on marine broadband seismic data waveform characteristics as claimed in claim 1, wherein the step 1) is to acquire original seismic data in a work area through an inclined cable and perform ghost wave suppression processing to obtain broadband seismic data.
3. The method for extracting attributes based on marine broadband seismic data waveform features according to claim 1, wherein the method for calculating the instantaneous phase attribute θ (t) in the step 2) is as follows: the signal is complex through Hilbert conversion, and after the signal x (t) is subjected to Hilbert conversion, the signal Z (t) is analyzed to be:
Z(t)=x(t)+iY(t)=A(t)eiθ(t)(1)
the instantaneous phase property θ (t) is found by the imaginary and real parts of the transformed signal, expressed as:
where a represents the amplitude of the signal, i represents the imaginary unit, t represents time, x (t) represents the input signal, y (t) represents the transformed imaginary component, and z (t) represents the analytic signal.
4. The method for extracting attributes based on marine broadband seismic data waveform features of claim 3, wherein the expression of the instantaneous phase gradient ω (t) in the step 3) is as follows:
where ω (t) represents the instantaneous phase gradient and θ (t) represents the instantaneous phase property.
5. The method as claimed in claim 4, wherein in the step 3), when the instantaneous phase gradient ω (t) is obtained, the wrapping information of the phase is retained, and the first derivative is directly obtained for the instantaneous phase.
6. The method for extracting attributes based on marine broadband seismic data waveform features according to claim 1, wherein the identification method in the step 4) is as follows: judging the positive polarity and the negative polarity of the seismic data, and if the seismic data is positive polarity, judging that the peak value position of the instantaneous phase gradient attribute corresponds to a negative reflection coefficient, and judging that the peak value interface of the instantaneous phase gradient attribute is an interface with wave impedance from large to small; and if the seismic data is negative, judging that the instantaneous phase gradient attribute peak interface is an interface with wave impedance from small to large if the instantaneous phase gradient attribute peak position corresponds to a positive reflection coefficient.
7. The method as claimed in claim 6, wherein the identification method in step 4) is to determine the positive and negative polarities of the seismic data based on the stable phase of the wideband seismic data of the work area to be identified.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010199904.7A CN111323816B (en) | 2020-03-20 | 2020-03-20 | Instantaneous phase gradient attribute extraction method based on ocean broadband seismic data waveform |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010199904.7A CN111323816B (en) | 2020-03-20 | 2020-03-20 | Instantaneous phase gradient attribute extraction method based on ocean broadband seismic data waveform |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111323816A true CN111323816A (en) | 2020-06-23 |
CN111323816B CN111323816B (en) | 2021-12-14 |
Family
ID=71167625
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010199904.7A Active CN111323816B (en) | 2020-03-20 | 2020-03-20 | Instantaneous phase gradient attribute extraction method based on ocean broadband seismic data waveform |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111323816B (en) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA1221160A (en) * | 1984-09-17 | 1987-04-28 | Clifford L. Shock | Seismic data processing |
US20030018436A1 (en) * | 2001-07-20 | 2003-01-23 | Stark Tracy Joseph | System for multi-dimensional data analysis |
CN102508295A (en) * | 2011-10-28 | 2012-06-20 | 西安交通大学 | Earthquake stratum thickness variation analysis method |
CN105093305A (en) * | 2014-05-15 | 2015-11-25 | 中国石油化工股份有限公司 | Method for determining position of reflection interface of seismic data |
CN105259579A (en) * | 2015-11-10 | 2016-01-20 | 中国地质大学(北京) | A high-amplitude shielding layer rejecting method based on seismic data instantaneous attributes |
US20180088249A1 (en) * | 2015-04-17 | 2018-03-29 | Schlumberger Technology Corporation | Thin Bed Tuning Frequency and Thickness Estimation |
CN107888525A (en) * | 2017-11-10 | 2018-04-06 | 中国电子科技集团公司第四十研究所 | Solution method for winding, the device of a kind of psk signal |
CN108710150A (en) * | 2018-05-22 | 2018-10-26 | 中国海洋石油集团有限公司 | A kind of earthquake random noise minimizing technology based on steady singular spectrum analysis |
CN109765616A (en) * | 2019-01-10 | 2019-05-17 | 中国海洋石油集团有限公司 | A kind of guarantor's width wave field extrapolation bearing calibration and system |
-
2020
- 2020-03-20 CN CN202010199904.7A patent/CN111323816B/en active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA1221160A (en) * | 1984-09-17 | 1987-04-28 | Clifford L. Shock | Seismic data processing |
US20030018436A1 (en) * | 2001-07-20 | 2003-01-23 | Stark Tracy Joseph | System for multi-dimensional data analysis |
CN102508295A (en) * | 2011-10-28 | 2012-06-20 | 西安交通大学 | Earthquake stratum thickness variation analysis method |
CN105093305A (en) * | 2014-05-15 | 2015-11-25 | 中国石油化工股份有限公司 | Method for determining position of reflection interface of seismic data |
US20180088249A1 (en) * | 2015-04-17 | 2018-03-29 | Schlumberger Technology Corporation | Thin Bed Tuning Frequency and Thickness Estimation |
CN105259579A (en) * | 2015-11-10 | 2016-01-20 | 中国地质大学(北京) | A high-amplitude shielding layer rejecting method based on seismic data instantaneous attributes |
CN107888525A (en) * | 2017-11-10 | 2018-04-06 | 中国电子科技集团公司第四十研究所 | Solution method for winding, the device of a kind of psk signal |
CN108710150A (en) * | 2018-05-22 | 2018-10-26 | 中国海洋石油集团有限公司 | A kind of earthquake random noise minimizing technology based on steady singular spectrum analysis |
CN109765616A (en) * | 2019-01-10 | 2019-05-17 | 中国海洋石油集团有限公司 | A kind of guarantor's width wave field extrapolation bearing calibration and system |
Non-Patent Citations (4)
Title |
---|
ALDO VESNAVER: "Instantaneous frequency and phase without unwrapping", 《GEOPHYSICS》 * |
乐友喜等: "用于地震反射界面识别的瞬时相位复合属性", 《物探化探计算技术》 * |
王俊茹: "《物探数字信号分析与处理技术》", 30 September 2009, 北京:地质出版社 * |
石颖等: "地震信号的复地震道分析及应用", 《地球物理学进展》 * |
Also Published As
Publication number | Publication date |
---|---|
CN111323816B (en) | 2021-12-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109633752B (en) | Offshore towing cable data self-adaptive ghost wave compression method based on three-dimensional fast Radon transformation | |
CN109581516B (en) | Denoising method and system for data of curvelet domain statistic adaptive threshold value ground penetrating radar | |
CN107144879A (en) | A kind of seismic wave noise-reduction method combined based on adaptive-filtering with wavelet transformation | |
CN105445801B (en) | A kind of processing method for eliminating 2-d seismic data random noise | |
Zhong‐Lai et al. | A SECOND‐ORDER SYNCHROSQUEEZING S‐TRANSFORM AND ITS APPLICATION IN SEISMIC SPECTRAL DECOMPOSITION | |
CN116520419B (en) | Hot fluid crack channel identification method | |
CN109581481B (en) | Portable high-frequency controllable seismic source seismic signal harmonic interference elimination method | |
CN110749925A (en) | Broadband reverse time migration imaging processing method | |
CN111323816B (en) | Instantaneous phase gradient attribute extraction method based on ocean broadband seismic data waveform | |
CN114089416A (en) | Method for estimating seismic wave attenuation gradient by utilizing Schrodinger equation | |
CN111239814B (en) | Shallow profile data mechanical interference suppression method based on same-phase axis frequency division tracking smoothing | |
CN111880222A (en) | Shearlet transformation-based seismic image minor fault identification enhancement method | |
CN108919345B (en) | Submarine cable land detection noise attenuation method | |
CN116203634A (en) | Ghost wave removing method based on low-rank constraint | |
CN110927781A (en) | Seismic data optimization method and system for high-precision seismic coherence attribute | |
CN114371505A (en) | Multi-wavelet inversion method and system based on seismic frequency division technology | |
CN103901470B (en) | Seabed settled layer reflection strength and time delay estimation method | |
CN107589456B (en) | Method and device for acquiring seismic data and computer readable storage medium | |
CN113777652B (en) | Method and device for extracting seismic source wavelet of marine seismic data and storage medium | |
Wang et al. | Estimate dip by combining gradient structure tensor and multi-window | |
CN114740530B (en) | Medium-high frequency quasi-linear noise suppression method and device based on hyperbolic time window constraint | |
CN113933902B (en) | High-frequency expansion method based on stereoscopic space wavelet | |
CN111880221B (en) | Novel VSP data seismic source wavelet self-adaptive extraction method based on Hilbert transform | |
CN114114422B (en) | Prestack seismic data noise elimination method based on directional multi-scale decomposition | |
CN110542921B (en) | Amplitude transverse change depicting method based on principal component analysis |
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 |