CN104360394A - Method of pressing ghost waves of marine variable-depth cable earthquake data - Google Patents

Method of pressing ghost waves of marine variable-depth cable earthquake data Download PDF

Info

Publication number
CN104360394A
CN104360394A CN201410717828.9A CN201410717828A CN104360394A CN 104360394 A CN104360394 A CN 104360394A CN 201410717828 A CN201410717828 A CN 201410717828A CN 104360394 A CN104360394 A CN 104360394A
Authority
CN
China
Prior art keywords
data
geological data
frequency spectrum
ripple
ghost
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
CN201410717828.9A
Other languages
Chinese (zh)
Other versions
CN104360394B (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.)
China National Offshore Oil Corp CNOOC
CNOOC Research Institute Co Ltd
Original Assignee
China National Offshore Oil Corp CNOOC
CNOOC Research Institute Co Ltd
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 China National Offshore Oil Corp CNOOC, CNOOC Research Institute Co Ltd filed Critical China National Offshore Oil Corp CNOOC
Priority to CN201410717828.9A priority Critical patent/CN104360394B/en
Publication of CN104360394A publication Critical patent/CN104360394A/en
Application granted granted Critical
Publication of CN104360394B publication Critical patent/CN104360394B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The invention relates to a method of pressing ghost waves of marine variable-depth cable earthquake data. The method comprises the following steps: firstly, transforming variable-depth cable earthquake data to a frequency domain; then, fitting frequency spectrum in the frequency domain, so that the difference between the frequency spectrum of the synthesized data and the frequency spectrum of the actually received data is extremely small, and the ghost wave pressing problem is transformed into a frequency spectrum fitting inversion problem. According to the characteristic of the variable-depth cable data, in an inversion process, a horizontal restrain item is introduced, so that the frequency spectrums on a trapped wave point position supplement with one another horizontally. The method of pressing ghost waves of variable-depth cable data disclosed by the invention is adopted, so that gradation is distinct; the transverse restrain item is adopted and a restrain factor is used to control the transverse restrain item, and therefore, the tapped wave effects of the ghost wave-pressed earthquake data are corrected in the frequency domain very well, the noises can be restrained very well in a time domain, and very good ghost wave pressing effects are achieved. The method disclosed by the invention can be directly used to a treatment process of acquired marine variable-depth cable earthquake data.

Description

A kind of method of marine varying depth cable geological data ghost wave pressure
Technical field
The present invention relates to exploratory engineering of off-shore petroleum/gas reservoir field, particularly about a kind of method of marine varying depth cable geological data ghost wave pressure.
Background technology
Seismic prospecting is the important means of oil-gas exploration, particularly exploratory engineering of off-shore petroleum/gas reservoir, because its drilling cost is high, geologic information and other data few, therefore in seismic prospecting oil-gas exploration at sea, there is prior effect.The geological data gathering wide-band is the pursuit of offshore seismic exploration always, but up to now, due to the restriction of various condition, this remains a unreachable target.The degree of depth at traditional marine acquisition towing cable place is fixing, and namely the geophone station ghost ripple trap frequency obtained is the maximum frequency that can recover.In order to make whole frequency higher, towing cable must be fixed in an as far as possible shallow degree of depth, but the noise of shallow-layer can affect the low-frequency range effective information of data.Towing cable is placed on certain depth to carry out gathering the noise that can reduce record, but geophone station ghost ripple can be introduced, make the spectral amplitude in earthquake frequency range there is trap phenomenon.Therefore, during tradition gathers, need, by close with sea level for trying one's best of towing cable placement, first terrible ripple trap wave point to be positioned at outside the frequency range required for imaging.
In order to overcome the defect of trap, proposing various different acquisition method at present, having comprised upper and lower cable acquisition method, two inspection acquisition method and varying depth cable acquisition method.The collection of varying depth cable is placed on by cable in the oblique line at a non-zero inclination angle, to make the trap wave point of terrible ripple decentralized.The people such as Soubaras re-start research to varying depth cable acquisition technique, and two features are added to it: the section 1) towing cable being optimized for a varying depth, to make the reflection horizon to all degree of depth, trap wave point is all dispersion, for shallow-layer part, the slope of towing cable is comparatively large in front portion, thus trap wave point is enough disperseed.2) a kind of geophone station ghost ripple drawing method is proposed for varying depth towing cable collection.The collection of varying depth cable can provide wider bandwidth for geological data, its object is the as far as possible dark smnr data obtaining maximum possible in low-frequency range by being placed by towing cable, different with the constant depth that tradition gathers, the degree of depth of towing cable changes along with offset distance, and object is to break the restriction of high frequency trap wave point to data bandwidth.The object of varying depth is that trap wave point is disperseed, each wave detector has a different trap frequency, thus make the stack result finally combining different wave detector there is no obvious trap wave point, varying depth cable acquisition technique is as a kind of emerging technology occurred recently, acquire several data volume all over the world after, can find out that varying depth cable data can provide higher-quality data: have higher resolution, imaging is more clear, and low frequency penetrability is better.
At present for the research of varying depth cable image data process, mainly concentrate in the compacting of geophone station ghost ripple, Soubaras (2010,2012) thinks, each geophone station degree of depth is change, and therefore varying depth acquisition technique makes geophone station ghost wavelength-division dispersion at different offset distance place.The characteristic of this dispersion can apply a kind of associating deconvolution technology, suppresses geophone station ghost ripple completely.Dechun Lin, Sablon, (2011) flow processs to varying depth cable data processing such as Yan Gao are optimized analysis, point out the broadband properties of varying depth cable data, and geophone station ghost ripple phenomenon specific to these data, make the process needing to consider constantly frequency band and geophone station ghost ripple in processing procedure, they give one compared with traditional data process in article, the flow process of slightly different varying depth cable process.Rebert, Sablon, Vidal etc. (2012) point out in varying depth cable processing procedure, should be noted that and keep the frequency at data frequency section two ends, and point out to face in varying depth cable processing procedure, as noise compacting, multiple reflection is suppressed, some challenges of velocity modeling and imaging aspect.Although existing terrible attenuation techniques can suppress some impacts of terrible ripple, the impact of terrible ripple can not be eliminated completely.
Summary of the invention
For the problems referred to above, the object of this invention is to provide one and can not only eliminate terrible wave action, and geological data frequency band can be widened, improve the method for the marine varying depth cable geological data ghost wave pressure of seismic data resolution.
For achieving the above object, the present invention takes following technical scheme: a kind of method of marine varying depth cable geological data ghost wave pressure, and it comprises the following steps: 1) take out collection to the geological data that varying depth cable receives, and forms common midpoint gather; 2) common midpoint geological data is carried out normal moveout correction, and the common midpoint geological data after normal moveout correction is carried out Fourier transform, obtain the frequency spectrum of the original earthquake data containing there's something fishy ripple; 3) according to the different depth of geophone station, the frequency spectrum of the original earthquake data containing there's something fishy ripple obtained is simplified, obtains H; 4) least Squares Spectral approximating method is adopted to calculate the frequency spectrum of the geological data after suppressing terrible ripple; 5) frequency spectrum of the geological data after the terrible ripple of compacting is carried out inverse Fourier transform, obtain the geological data after suppressing terrible ripple.
Described step 3) in H (f)=1-e -i2 π f τ, wherein, τ is the whilst on tour of seismic event from wave detector to water meter.
Described step 4) adopt least Squares Spectral approximating method to calculate the terrible ripple of compacting after the frequency spectrum detailed process of geological data be:
Ghost wave pressure system is equivalent to the minimal value solving following formula;
||HS r-S *|| (1)
Wherein, S *the matrix representation of the geological data that known varying depth cable actual acquisition arrives, S rit is the matrix representation of the seismic signal that subsurface reflective is returned;
Formula (1) basis adds the bound term embodying different track data and mutually supplement, then formula (1) becomes:
||HS r-S *||+α||MS r|| (2)
Wherein, || MS r|| be Lateral constrained condition, α is constraint factor;
The basis of formula (2) adds stable functional obtain:
||HS r-S *||+α||MS r||+β||IS r|| (3)
Wherein, β represents stable factor;
The minimal value of solution formula (3), to formula (3) about S rdifferentiate and make derivative equal zero:
H T(HS r-S *)+αM TMS r+βIS r=0 (4)
Arrange:
(H TH+αM TM+βI)S r=H TS *(5)
According to given α and the β parameter of the feature of geological data, α and β is substituted into formula (5) and obtains a system of equations, separate this system of equations and can solve and obtain S r, obtain the frequency spectrum of the geological data after suppressing terrible ripple.
The value of α and β two parameters is most important, α and β determines according to the signal to noise ratio (S/N ratio) of seismic data, if the signal to noise ratio (S/N ratio) of seismic data is b, then α is taken as β is taken as
The present invention is owing to taking above technical scheme, it has the following advantages: 1, the present invention is owing to taking full advantage of the feature of varying depth cable data, prestack CMP road collection (common midpoint gather) is achieved by spectrum matching the compacting of terrible ripple, different according to each seismic trace trap wave point in pressing process, CMP road collection complements each other, therefore the better trap frequency must recovering every track data, has widened seismic signal frequency spectrum, improves the resolution of seismic data.2, the present invention is due to while compensation trap frequency, phase spectrum is corrected, therefore compared with variable step the Method of Deconvolution, not only compensate for notch band preferably in frequency field, also terrible ripple has been suppressed preferably in time domain, and there is no the adjoint noise after variable step deconvolution process, simultaneously abundanter than the result low frequency of variable step the Method of Deconvolution.3, the present invention owing to carrying out terrible wave pressure system on prestack CMP road collection, the still prestack road collection obtained after suppressing terrible ripple, remain AVO (the Amplitude Versus Offset of seismic data, the change of amplitude offset distance) feature, therefore can carry out AVO analysis and prestack inversion on this basis, improve the accuracy of petroleum-gas prediction.4, the present invention is due to by carrying out normal moveout correction to concentrically road collection, then the frequency spectrum of common midpoint gather is calculated, least Squares Spectral approximating method is utilized to obtain the geological data after suppressing terrible ripple, thus reach the object of suppressing terrible ripple, therefore can suppress the terrible ripple in seismologic record preferably, reach with without the on all four effect of terrible Bo CMP road collection.Therefore, the present invention can be widely used in the process of energy exploration exploitation.
Accompanying drawing explanation
Fig. 1 is the schematic flow sheet of the inventive method;
Fig. 2 is the Numerical Experiment result schematic diagram of the embodiment of the present invention, wherein, figure (a) be synthesis without Gui Bo CMP road collection, figure (b) be synthesis containing terrible Bo CMP road collection, figure (c) is that employing the present invention suppresses terrible ripple Hou CMP road collection;
Fig. 3 is the Numerical Experiment spectral contrast schematic diagram of the embodiment of the present invention, and wherein, figure (a) is the spectral amplitude containing Gui Bo CMP road collection of synthesis, and figure (b) is the spectral amplitude adopting the present invention to suppress terrible Bo Hou CMP road collection;
Fig. 4 is the real data process experimental result schematic diagram of the embodiment of the present invention, wherein, figure (a) is raw data CMP road collection, and figure (b) adopts variable step deconvolution process Hou CMP road collection, and figure (c) adopts process Hou CMP road of the present invention collection;
Fig. 5 is the real data process stacked section contrast schematic diagram of the embodiment of the present invention, and wherein, figure (a) is the stacked section after adopting variable step deconvolution process, and figure (b) is the stacked section after adopting process of the present invention.
Embodiment
Below in conjunction with drawings and Examples, the present invention is described in detail.
As shown in Figure 1, the method for marine varying depth cable geological data ghost wave pressure of the present invention, comprises the following steps:
1) collection is taken out to the geological data that varying depth cable receives, form common midpoint gather;
2) common midpoint geological data is carried out normal moveout correction, and the common midpoint geological data after normal moveout correction is carried out Fourier transform, obtain the frequency spectrum of the original earthquake data containing there's something fishy ripple;
3) according to the different depth of geophone station, the frequency spectrum of the original earthquake data containing there's something fishy ripple obtained is simplified, obtains H;
S r(t)=w(t)*r(t) (1)
Wherein, S rt () represents the seismic signal that subsurface reflective is returned, be equivalent to not containing the seismic signal of terrible ripple; R (t) represents the reflection coefficient on stratum, and w (t) represents seismic wavelet.The reflection coefficient of water meter is r s=-1;
Due to the reflection of water meter, seismic signal S (t) of the actual reception of varying depth cable is:
S(t)=S r(t)+r sS r(t-τ)=S r(t)-S r(t-τ) (2)
Wherein, τ is the whilst on tour of seismic event from wave detector to water meter, in frequency field:
S(f)=S r(f)-S r(f)e -i2πfτ(3)
Be reduced to further:
S(f)=S r(f)(1-e -i2πfτ)=H(f)S r(f) (4)
Wherein,
H(f)=1-e -i2πfτ(5)
|H(f)|=2sin(πfτ) (6)
As can be seen from formula (6), the trap effect that the reflection due to water meter causes, the whilst on tour τ from wave detector to water meter is relevant for trap wave point position and seismic event; In varying depth cable data, because the geophone station degree of depth is change, thus the whilst on tour τ of seismic event from wave detector to water meter is also different, and trap wave point position is also different;
Suppose that the degree of depth of geophone station is respectively h j(j=1,2, ‖, n), n represents the number of channels in each common midpoint gather, then:
S j ( f ) = H j ( f ) S j r ( f ) = ( 1 - e - i 2 π fτ J ) S j r ( f ) - - - ( 7 )
Wherein, S jf () represents the frequency spectrum containing there's something fishy Boj road geological data, represent the frequency spectrum of compacting terrible Bo Houj road geological data;
Formula (7) is write as the form of matrix, then:
S=HS r(8)
4) least Squares Spectral approximating method is adopted to calculate the frequency spectrum of the geological data after suppressing terrible ripple;
Suppress terrible ripple and be equivalent to the geological data that known varying depth cable actual acquisition arrives seismic signal Sjr (f) that inverting subsurface reflective is returned, namely terrible wave pressure system is equivalent to the minimal value solving following problem;
||HS r-S *|| (9)(9)
Wherein, || HS r-S *|| be equivalent to carry out matching by the spectrum of geological data of synthesizing and the spectrum of real data;
Because varying depth cable data will be complemented each other by different trap wave points, on formula (9) basis, thus add the bound term embodying different track data and mutually supplement, then formula (9) becomes:
||HS r-S *||+α||MS r|| (10)
Wherein, || MS r|| be Lateral constrained condition, α is constraint factor, and the size of α can affect the power of lateral confinement;
In order to overcome the instability of refutation process when formula (10) is asked minimum, the basis of formula (10) adding stable functional and obtains:
||HS r-S *||+α||MS r||+β||IS r|| (11)
Wherein, β represents stable factor, the minimal value of solution formula (11), and can realize the compacting of varying depth cable data ghost ripple, concrete solution procedure is:
In order to the minimal value of solution formula (11), to formula (11) about S rdifferentiate and make derivative equal zero:
H T(HS r-S *)+αM TMS r+βIS r=0 (12)
Arrange:
(H TH+αM TM+βI)S r=H TS *(13)
According to feature given parameters α and β of geological data, α and β is substituted into formula (13) and obtains a system of equations, separate this system of equations and can solve and obtain S r, obtain the frequency spectrum of the geological data after suppressing terrible ripple;
In practical operation, α and β is two vital parameters, chooses the improper result of calculation that may cause is failed.General α and β can determine according to the signal to noise ratio (S/N ratio) of seismic data, if the signal to noise ratio (S/N ratio) of seismic data is b, then α can be taken as β can be taken as
5) frequency spectrum of the geological data after the terrible ripple of compacting is carried out inverse Fourier transform, obtain the geological data after suppressing terrible ripple.
As shown in Figure 2-5, the varying depth cable data acquisition of field the actual acquisition method of marine varying depth cable geological data ghost wave pressure of the present invention and existing method are carried out Numerical Experiment respectively, and adopt the present invention to process real data, below logarithm value simulated experiment and the result that real data processes is described respectively:
As shown in Figure 2, Numerical Experiment result shows the present invention and has suppressed terrible ripple in seismologic record preferably, and result is with completely the same without terrible Bo CMP road collection.As shown in Figure 3, Numerical Experiment spectral contrast shows the present invention and has recovered the trap effect of terrible ripple preferably.As shown in Figure 4, real data process experimental result display variable step the Method of Deconvolution creates with noise while the terrible ripple of compacting, and the present invention does not produce with noise.As shown in Figure 5, find out that variable step the Method of Deconvolution corrects the trap effect of terrible ripple preferably on spectral amplitude from figure (a), but in time domain, the method can be introduced with noise in the process of the terrible ripple of compacting, reduces the signal to noise ratio (S/N ratio) of geological data; Find out that the present invention not only corrects trap effect preferably on spectral amplitude from figure (b), time domain does not also occur noise, and in processing procedure, have the effect of compacting noise.
In sum, existing mirrored offset associating the Method of Deconvolution is a kind of typical poststack ghost ripple drawing method, this method can suppress terrible ripple preferably for imaging, corrects trap effect, but this method can not be used for the reservoir prediction process of AVO analysis or prestack inversion after treatment.And adopt method of the present invention, prestack road collection carries out terrible wave pressure system, the still prestack road collection obtained after suppressing terrible ripple, the reservoir prediction process such as AVO analysis and prestack inversion can be carried out further.Therefore, can find out that the present invention effectively can suppress terrible ripple on prestack road collection by real data process, compensate terrible ripple and fall into frequently.Relative to conventional variable step deconvolution ghost ripple drawing method, no matter in terrible wave pressure system, or in signal to noise ratio (S/N ratio), be all significantly improved.
The various embodiments described above are only for illustration of the present invention, and wherein each step etc. all can change to some extent, and every equivalents of carrying out on the basis of technical solution of the present invention and improvement, all should not get rid of outside protection scope of the present invention.

Claims (4)

1. a method for marine varying depth cable geological data ghost wave pressure, it comprises the following steps:
1) collection is taken out to the geological data that varying depth cable receives, form common midpoint gather;
2) common midpoint geological data is carried out normal moveout correction, and the common midpoint geological data after normal moveout correction is carried out Fourier transform, obtain the frequency spectrum of the original earthquake data containing there's something fishy ripple;
3) according to the different depth of geophone station, the frequency spectrum of the original earthquake data containing there's something fishy ripple obtained is simplified, obtains H;
4) least Squares Spectral approximating method is adopted to calculate the frequency spectrum of the geological data after suppressing terrible ripple;
5) frequency spectrum of the geological data after the terrible ripple of compacting is carried out inverse Fourier transform, obtain the geological data after suppressing terrible ripple.
2. the method for a kind of marine varying depth cable geological data ghost wave pressure as claimed in claim 1, is characterized in that: described step 3) in H (f)=1-e -i2 π f τ, wherein, τ is the whilst on tour of seismic event from wave detector to water meter.
3. the method for a kind of marine varying depth cable geological data ghost wave pressure as claimed in claim 1 or 2, is characterized in that: described step 4) adopt least Squares Spectral approximating method to calculate the terrible ripple of compacting after the frequency spectrum detailed process of geological data be:
Ghost wave pressure system is equivalent to the minimal value solving following formula;
||HS r-S *|| (1)
Wherein, S *the matrix representation of the geological data that known varying depth cable actual acquisition arrives, S rit is the matrix representation of the seismic signal that subsurface reflective is returned;
Formula (1) basis adds the bound term embodying different track data and mutually supplement, then formula (1) becomes:
||HS r-S *||+α||MS r|| (2)
Wherein, || MS r|| be Lateral constrained condition, α is constraint factor;
The basis of formula (2) adds stable functional obtain:
||HS r-S *||+α||MS r||+β||IS r|| (3)
Wherein, β represents stable factor;
The minimal value of solution formula (3), to formula (3) about S rdifferentiate and make derivative equal zero:
H T(HS r-S *)+αM TMS r+βIS r=0 (4)
Arrange:
(H TH+αM TM+βI)S r=H TS *(5)
According to given α and the β parameter of the feature of geological data, α and β is substituted into formula (5) and obtains a system of equations, separate this system of equations and can solve and obtain S r, obtain the frequency spectrum of the geological data after suppressing terrible ripple.
4. the method for a kind of marine varying depth cable geological data ghost wave pressure as claimed in claim 3, it is characterized in that: the value of α and β two parameters is most important, α and β determines according to the signal to noise ratio (S/N ratio) of seismic data, if the signal to noise ratio (S/N ratio) of seismic data is b, then α is taken as β is taken as
CN201410717828.9A 2014-12-01 2014-12-01 Method of pressing ghost waves of marine variable-depth cable earthquake data Active CN104360394B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410717828.9A CN104360394B (en) 2014-12-01 2014-12-01 Method of pressing ghost waves of marine variable-depth cable earthquake data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410717828.9A CN104360394B (en) 2014-12-01 2014-12-01 Method of pressing ghost waves of marine variable-depth cable earthquake data

Publications (2)

Publication Number Publication Date
CN104360394A true CN104360394A (en) 2015-02-18
CN104360394B CN104360394B (en) 2017-02-22

Family

ID=52527674

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410717828.9A Active CN104360394B (en) 2014-12-01 2014-12-01 Method of pressing ghost waves of marine variable-depth cable earthquake data

Country Status (1)

Country Link
CN (1) CN104360394B (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104977615A (en) * 2015-07-01 2015-10-14 中国海洋石油总公司 Model-statistics-pickup-based multiple suppression method of deep-sea OBC data
CN106199699A (en) * 2016-06-30 2016-12-07 恒泰艾普(北京)能源科技研究院有限公司 Utilize the method that transfer matrix method removes ghost ripple
CN106526677A (en) * 2016-10-26 2017-03-22 中海石油(中国)有限公司 Marine self-adaptive ghost reflection-suppressing broadband reverse time migration imaging method
CN106610507A (en) * 2016-12-01 2017-05-03 中国石油天然气集团公司 Method and device for removing ghost reflection interference
CN106896409A (en) * 2017-03-14 2017-06-27 中国海洋石油总公司 A kind of varying depth cable ghost ripple drawing method based on wave equation boundary values inverting
CN108919357A (en) * 2018-05-16 2018-11-30 中国海洋石油集团有限公司 A kind of ghost reflection drawing method based on frequency spectrum reconfiguration
CN108983284A (en) * 2018-06-22 2018-12-11 中国石油大学(华东) A kind of domain f-p ghost reflection drawing method suitable for marine tiltedly cable data
CN110967735A (en) * 2018-09-28 2020-04-07 中国石油化工股份有限公司 Self-adaptive ghost wave suppression method and system
CN112817047A (en) * 2020-12-31 2021-05-18 北京东方联创地球物理技术有限公司 Ocean earthquake self-adaptive ghost wave removing method and device, electronic equipment and medium
CN113514889A (en) * 2021-07-13 2021-10-19 中山大学 Processing method for improving low-frequency signal energy in ocean deep reflection seismic data

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102169189B (en) * 2011-01-19 2015-01-14 中国海洋石油总公司 Method for eliminating deep water interlamination multiple
CN102636811B (en) * 2012-04-10 2014-01-29 恒泰艾普石油天然气技术服务股份有限公司 Eliminating method of multiple waves in bidimensional seismic data on sea
CN103308943B (en) * 2013-05-10 2015-11-18 中国石油天然气股份有限公司 The method of interlayer multiple attenuation in the process of a kind of oceanic earthquake data and device
CN103308944B (en) * 2013-05-20 2015-10-14 中国石油天然气股份有限公司 The method of terrible wave attenuation in the process of a kind of oceanic earthquake data and device

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
A.KEMAL öZDEMIR等: "Optimized deghosting of over/under towed-streamer data in the presence of noise", 《THE LEADING EDGE》 *
ROBERT SOUBARAS等: "variable-depth streamer acquisition:broadband data for imaging and inversion", 《GEOPHYSICS》 *
张振波 等: "白云凹陷深水复杂构造区斜缆地震资料处理关键技术及应用", 《石油物探》 *
王芳芳 等: "基于逆散射级数法的鬼波压制方法", 《地球物理学报》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104977615A (en) * 2015-07-01 2015-10-14 中国海洋石油总公司 Model-statistics-pickup-based multiple suppression method of deep-sea OBC data
CN106199699A (en) * 2016-06-30 2016-12-07 恒泰艾普(北京)能源科技研究院有限公司 Utilize the method that transfer matrix method removes ghost ripple
CN106526677A (en) * 2016-10-26 2017-03-22 中海石油(中国)有限公司 Marine self-adaptive ghost reflection-suppressing broadband reverse time migration imaging method
CN106610507A (en) * 2016-12-01 2017-05-03 中国石油天然气集团公司 Method and device for removing ghost reflection interference
CN106610507B (en) * 2016-12-01 2018-12-25 中国石油天然气集团公司 A kind of method and apparatus of removal ghosting interference
CN106896409B (en) * 2017-03-14 2018-12-07 中国海洋石油集团有限公司 A kind of varying depth cable ghost reflection drawing method based on wave equation boundary values inverting
CN106896409A (en) * 2017-03-14 2017-06-27 中国海洋石油总公司 A kind of varying depth cable ghost ripple drawing method based on wave equation boundary values inverting
CN108919357A (en) * 2018-05-16 2018-11-30 中国海洋石油集团有限公司 A kind of ghost reflection drawing method based on frequency spectrum reconfiguration
CN108919357B (en) * 2018-05-16 2019-10-11 中国海洋石油集团有限公司 A kind of ghost reflection drawing method based on frequency spectrum reconfiguration
CN108983284A (en) * 2018-06-22 2018-12-11 中国石油大学(华东) A kind of domain f-p ghost reflection drawing method suitable for marine tiltedly cable data
CN110967735A (en) * 2018-09-28 2020-04-07 中国石油化工股份有限公司 Self-adaptive ghost wave suppression method and system
CN112817047A (en) * 2020-12-31 2021-05-18 北京东方联创地球物理技术有限公司 Ocean earthquake self-adaptive ghost wave removing method and device, electronic equipment and medium
CN112817047B (en) * 2020-12-31 2021-10-08 北京东方联创地球物理技术有限公司 Ocean earthquake self-adaptive ghost wave removing method and device, electronic equipment and medium
CN113514889A (en) * 2021-07-13 2021-10-19 中山大学 Processing method for improving low-frequency signal energy in ocean deep reflection seismic data

Also Published As

Publication number Publication date
CN104360394B (en) 2017-02-22

Similar Documents

Publication Publication Date Title
CN104360394A (en) Method of pressing ghost waves of marine variable-depth cable earthquake data
Özdemir et al. Optimized deghosting of over/under towed-streamer data in the presence of noise
CN102305941B (en) Method for determining stratum stack quality factor by direct scanning of prestack time migration
WO2018026875A1 (en) Surface consistent statics solution and amplification correction
Lin et al. Optimizing the processing flow for variable-depth streamer data
CN112305612B (en) High-resolution complex spectrum decomposition time-frequency space domain amplitude variation correction method along with offset distance
CN103913770A (en) Method for processing seismic data based on VSP data
van Gestel et al. Atlantis ocean bottom nodes time-lapse observations
Zhang et al. Interval Q inversion based on zero-offset VSP data and applications
Piazza et al. M-factorial kriging-an efficient aid to noisy seismic data interpretation
Sablon et al. Synchronized Multi-level Source and Variable-depth Streamer-A Combined Ghost-free Solution for Broadband Marine Data
AU2011248910B2 (en) System and method for compensating time and offset varying near-surface effects in seismic data background
Welsh et al. Spectral decomposition of seismic reflection data to detect gas related frequency anomalies
Saunders et al. Improved stratigraphic interpretation using broadband processing–Sergipe Basin, Brazil
Latiff et al. Seismic data analysis to the converted wave acquisition: A case study in offshore Malaysia
Qingfeng et al. Progress and Application of 3DVSP in Shengli Oilfield
van Gestel et al. WATS on NATS time-lapse: Gulf of Mexico case study, King and Horn Mountain fields
Saunders et al. A 2-Stage Approach to Broadband Processing for Improved Stratigraphic Interpretation in the Sergipe Basin, Brazil
Urosevic et al. Acquisition of the Nested 3D seismic survey at Harvey
Sun et al. P-wave azimuthal velocity inversion for fracture detection in an east Texas gas field
Pirogova et al. AVA-study of shallow marine unconsolidated sediments using UHR multichannel seismoacoustic dataset acquired with a deep-towed system
Velasco et al. Near-Surface Characterization Using Point Receiver Data
Lansley Considerations for optimum 3-D survey design, acquisition and processing
Backus et al. OBC sensor response and calibrated reflectivity
Denney Carbonate-Reservoir Characterization: Wide-Azimuth-Seismic Processing

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: 100010 Beijing, Chaoyangmen, North Street, No. 25, No.

Co-patentee after: CNOOC research institute limited liability company

Patentee after: China Offshore Oil Group Co., Ltd.

Address before: 100010 Beijing, Chaoyangmen, North Street, No. 25, No.

Co-patentee before: CNOOC Research Institute

Patentee before: China National Offshore Oil Corporation