CN101592738B - Method for identifying two-dimensional up-and-down sea-bottom multiple - Google Patents

Method for identifying two-dimensional up-and-down sea-bottom multiple Download PDF

Info

Publication number
CN101592738B
CN101592738B CN2008101141340A CN200810114134A CN101592738B CN 101592738 B CN101592738 B CN 101592738B CN 2008101141340 A CN2008101141340 A CN 2008101141340A CN 200810114134 A CN200810114134 A CN 200810114134A CN 101592738 B CN101592738 B CN 101592738B
Authority
CN
China
Prior art keywords
sea
big gun
wave field
minimum
elevation
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.)
Active
Application number
CN2008101141340A
Other languages
Chinese (zh)
Other versions
CN101592738A (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.)
BGP Inc
Original Assignee
BGP Inc
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 BGP Inc filed Critical BGP Inc
Priority to CN2008101141340A priority Critical patent/CN101592738B/en
Publication of CN101592738A publication Critical patent/CN101592738A/en
Application granted granted Critical
Publication of CN101592738B publication Critical patent/CN101592738B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention relates to a method for identifying two-dimensional up-and-down sea-bottom multiple by wave field continuation, comprising the following steps: adopting seismic data collected by maritime towing line seismic exploration and sea floor altitude figures obtained by field survey and indoor acquisition; starting from the minimum shot size and shot gather, obtaining the maximum sea floor altitude and the minimum sea floor altitude and calculating the difference between the maximum sea floor altitude and the minimum sea floor altitude according to the sea floor altitude to which the array of the shot gather belongs; predicting the sea-bottom multiple according to horizontal or up-and-down wave movement.The advantages of the invention lie in two aspects: first, giving overall consideration to accuracy and efficiency and offering different solutions to different situations at the sea bottom: adopting the phase displacement technology with the highest accuracy and efficiency to the gentle parts and adopting the wave field insertion technology with high accuracy to fluctuant parts; second, introducing the wave field insertion technology to the identification of the up-and-down sea-bottom multiple for the first time.

Description

Method for identifying two-dimensional up-and-down sea-bottom multiple
Technical field
The present invention relates to geophysical exploration technology, belong to the technology category of the regular interference wave of compacting in the seismic data processing procedure, is a kind of method for identifying two-dimensional up-and-down sea-bottom multiple.
Background technology
In the seismic prospecting,, when earthquake data acquisition, can record more intense sea-bottom multiple at sea because the seabed is a strong reflecting interface between fluid and the solid.The existence of sea-bottom multiple can mislead the explanation personnel on the one hand and give the structure elucidation result who makes mistake, and can disturb primary reflection on the other hand, influences the structure imaging of primary reflection and utilizes primary reflection to carry out Reservoir Analysis.Therefore suppress sea-bottom multiple in the method for marine seismic data and be always a gordian technique during method for marine seismic data is handled.
The method of the sea-bottom multiple in the compacting method for marine seismic data generally is divided into two big classes.One class is based on the kinematics difference of sea-bottom multiple and primary reflection, multiple reflection is separated with primary reflection, thereby reach the purpose of multiple suppression; Another kind ofly be based on wave theory simulation or dope the model trace collection of sea-bottom multiple, utilize multiple reflection model trace collection then, adopt adaptive approach that the multiple reflection in the raw data is deducted as multiple reflection.
Method based on the kinematics difference multiple suppression of sea-bottom multiple and primary reflection mainly comprises predictive deconvolution method, F-K conversion filter method, Karhunen-Loeve transformation filter method, Tau-p conversion filter method and Tau-p territory The Method of Deconvolution.Based on wave theory carry out that the multiple reflection forecast method mainly contains that wave equation is just drilled, the wave field extrapolation method, based on the method for WRW model, based on the Born Series Method and the Kirchhoff Series Method of backscattering theory.Multiple reflection self-adaptation subraction comprises that mainly minimum two multipliers are according to matching method, pattern-recongnition method and independent component analysis method.
Wave field extrapolation multiple reflection forecasting techniques can be divided into the sea level by the geological data treated side and to be the wave field extrapolation multiple reflection forecasting techniques of treated side and to be the wave field extrapolation multiple reflection forecasting techniques of treated side with the sea bottom surface.Numerical method by wave field extrapolation can be divided into integral method wave field extrapolation multiple reflection forecasting techniques and phase shift method wave field extrapolation technology.Can be divided into single stage method, two-step approach and three-step approach wave field extrapolation multiple reflection forecasting techniques by the continuation step.
The wave field extrapolation multiple reflection forecasting techniques of traditional with the sea level is treated side is shown in Fig. 1 (a), general adopt integral method earlier geological data from the sea level continuation to the seabed, again continuation to the geological data in seabed from the seabed continuation to the sea level, so just doped sea-bottom multiple.This method has obviously been ignored the reflection that D is ordered, and need use the speed of medium under the seabed during integration.Need two speed when overcoming integration, people generally are approximately a monocline or flat seabed to the fluctuating seabed in an arrangement, shown in Fig. 1 (b).Fig. 1 (a) and (b) just different on the mode of processing seabed, but they need two step continuation to finish the prediction of sea-bottom multiple, and in order to reduce the inconvenience that the two-step approach continuation brings, people had proposed the single stage method continuation afterwards again.This method is that the big gun collection is placed on the mirror image face of sea level with respect to level and smooth sea bottom surface, utilizes upward traveling wave just can obtain the multiple reflection model trace collection of this big gun collection to the big gun collection to the sea level from the continuation of mirror image face like this, sees Fig. 1 (c).
Obviously, sea-bottom multiple Forecasting Methodology in the past based on wave field extrapolation, when acutely rising and falling in the seabed, because can not, therefore there is bigger error between the kinematics character of the sea-bottom multiple of prediction and actual multiple reflection.
Summary of the invention
The present invention seeks under the prerequisite of known submarine elevation and seawater speed, to fluctuating sea-bottom multiple prediction, provide a kind of method for identifying two-dimensional up-and-down sea-bottom multiple of wave field extrapolation of correct realization fluctuating seabed the wave field extrapolation technology introduction of fluctuating interface.
The invention provides following solution:
1) geological data, field survey or the indoor submarine elevation data of asking for of employing marine streamer seismic acquisition;
2) from minimum big gun gunfire collection, read also (can not occur subjective action word) submarine elevation of arranging according to big gun collection place of shot gather data, obtain its maximum submarine elevation and minimum submarine elevation, it is poor to calculate minimum and maximum submarine elevation;
3) when minimum and maximum submarine elevation difference during, by horizontal seabed prediction sea-bottom multiple less than 5 meters:
p n(x, z=0 t) are the n big gun geological data that is positioned at (z=0) on the sea level,
m n(x, z=0 t) are the sea-bottom multiple of determining of the correspondence n big gun that is positioned at (z=0) on the sea level obtained,
Wherein n is a big gun number, and x is the distance of shot point to acceptance point, and z is the elevation coordinate, and t is a time coordinate;
The multiple reflection in the described definite horizontal seabed of step 3) is meant:
(1) to p n(x, z=0 t) carry out two-dimentional fourier transform about x and t, obtain P n(k x, z=0, ω), k wherein xBe respectively horizontal direction circular wavenumber and circular frequency with ω; k z = ω 2 / c 2 - k x 2 Be the vertical direction circular wavenumber, c is the speed of seawater;
(2) by formula M n ( k x , z = 0 , ω ) = P n ( k x , z = 0 , ω ) e - 2 i k z | H ) Obtain the multiple reflection that frequency-wavenumber domain is surveyed, wherein, M n(k x, z=0 is the multiple reflection of the n big gun of frequency-wavenumber domain prediction ω), H is the mean value of minimum and maximum submarine elevation;
(3) to M n(k x, z=0 ω) carries out two-dimentional Fu Shi inverse transformation, obtains the multiple reflection m of the n big gun of time domain prediction n(x, z=0, t).
4) when minimum and maximum submarine elevation difference during, by fluctuating seabed prediction sea-bottom multiple greater than 5 meters:
p n(x, z=0 t) are the n big gun geological data that is positioned on the sea level (z=0),
p n(x, z=Z 0, t) be n big gun seismic wave field of (z=Z0) on bottom surface, shallow sea,
p n(x, z=Z N, t) be the n big gun on bottom surface, deep-sea (seismic wave field of z=ZN,
p n(x, z i, t) be the n big gun at z i=z 0+ (i-1) Δ z (i=1,2 ... the N+1) seismic wave field on the face,
q n(z i, be that the n big gun is at z t) i=z 0+ (i-1) Δ z (i=1,2 ... the N+1) record of bottom, face Shanghai,
m n(x, z=0 t) are the sea-bottom multiple of the prediction of the n big gun on sea level (z=0) obtained,
Wherein n is a big gun number, and x is the distance of shot point to acceptance point, and z is the elevation coordinate, and t is a time coordinate, and Δ z is the step size of depth direction.
The multiple reflection in the described definite fluctuating seabed of step 3) is meant:
(1) to p n(x, z=0 t) carry out fourier transform about x and t, obtain P n(k x, z=0, ω), k wherein xWith ω for being respectively horizontal direction circular wavenumber and circular frequency; k z = ω 2 / c 2 - k x 2 Be the vertical direction circular wavenumber, c is the speed of seawater;
(2) by formula P n ( k x , z = Z 0 , ω ) = P n ( k x , z = 0 , ω ) e - i k z | Z 0 | Obtain the seismic wave field of the frequency-wavenumber domain of (z=Z0) on bottom surface, shallow sea;
(3) from bottom surface, shallow sea (i=1) to bottom surface, deep-sea (i=N+1), by formula P n(k x, z=z I+1, ω)=P n(k x, z=z i, ω) e -2ik Δ zCalculate z i=z 0+ (i-1) the descending wave-wave field on the Δ z face is to P n(k x, z=z I+1, ω) carry out two-dimentional Fu Shi inverse transformation, and the following traveling-wave field q at place, record seabed n(z I+1, t);
(4) from bottom surface, deep-sea (i=N+1) to bottom surface, shallow sea (i=1), by formula p n(x, z=z i, t)=p n(x, z=z i, t)+q (z i, t) carry out wave field and be bumped into, by formula: P n(k x, z=z I-1, ω)=P n(k x, z=z i, ω) e -ik Δ zCarry out the upstream wave field continuation;
(5) by formula M n ( k x , z = 0 , ω ) = P n ( k x , z = Z 0 , ω ) e - i k z | Z 0 | Obtain the sea-bottom multiple of the frequency-wavenumber domain that on the sea level, writes down;
(6) to M n(k x, z=0 ω) carries out two-dimentional Fu Shi inverse transformation, obtains the multiple reflection of the n big gun of time domain prediction, m n(x, z=0, t).
Advantage one of the present invention is to have taken all factors into consideration precision and efficient, and the different situations in seabed are carried out different processing, locates to adopt all the highest phase shift technique of precision and efficient gently in the seabed; At the submarine relief place, adopt the very high wave field of precision to be bumped into technology; The 2nd, first wave field is bumped into the technology sea-bottom multiple that is applied to rise and fall and determines.
Description of drawings
The sea-bottom multiple that Fig. 1 is traditional is determined the method synoptic diagram;
(a) traditional multiple reflection Forecasting Methodology based on wave field extrapolation,
(b) the approximate multiple reflection Forecasting Methodology of inclination or flat sea bottom surface based on wave field extrapolation,
(c) the multiple reflection Forecasting Methodology that begins from the mirror image face based on wave field extrapolation
Fig. 2 the present invention sea-bottom multiple that rises and falls is determined the method synoptic diagram;
(a) according to submarine elevation the seabed is divided into horizontal seabed with the fluctuating seabed and adopt different wave field extrapolation algorithms
(b) corresponding fluctuating benthic division adopts wave field to be bumped into method to carry out wave field extrapolation
The big gun collection record of Fig. 3 input;
The sea-bottom multiple of the prediction of Fig. 4 output.
Specific embodiments
The present invention adopts streamer seismic data, field survey or the indoor submarine elevation data of asking for of marine streamer seismic acquisition, on the basis of taking all factors into consideration computational accuracy and counting yield, utilize phase shift horizontal interface wave field extrapolation method and wave field to be bumped into the sea-bottom multiple that method fluctuating interface wave field continuation method predicts that the seabed produces, obtain the sea-bottom multiple geological data of prediction.
Definition:
Bottom, the shallow sea elevation (Shallow bottom elevation) of SBE or Z0---;
The elevation (Deep bottom elevation) of the bottom, deep-sea of DBE or ZN-;
Poor (the Bottom elevation difference) of the most shallow and deep-sea floor elevation of BED---;
The mean value (Mean of bottom elevation) of the most shallow and deep-sea floor elevation of MBE--;
The most shallow and the deep-sea floor elevation that LBED---user is given poor is used to differentiate (the Limitation of the bottom elevation difference) that flat seabed still is the fluctuating seabed;
The wave field extrapolation step-length that DZ---user is given.
The specific implementation step of fluctuating sea-bottom multiple prediction is as follows:
1, from minimum big gun gunfire collection, reads shot gather data p n(x, z=0, t) (see figure 3) and submarine elevation data.p n(t) for being positioned at the n big gun geological data of (z=0) on the sea level, wherein n is a big gun number for x, z=0, and x is the distance of shot point to acceptance point, and z is the elevation coordinate, and t is a time coordinate;
2, the submarine elevation section of arranging according to big gun collection place obtains the altitude figures SBE and the DBE of the most shallow and bottom, deep-sea.
3, calculate the dark and the most shallow submarine elevation difference BED of place.As shown in Figure 2, when difference of elevation is less than BEI at the bottom of the darkest and shallow sea, the multiple reflection in prediction level seabed according to the following steps.
(1) the mean value MBE of the minimum and maximum submarine elevation of calculating puts H=MBE;
(2) to p n(x, z=0 t) carry out two-dimentional fourier transform about x and t, obtain P n(k x, z=0, ω), k wherein xBe respectively the circular wavenumber and the circular frequency of horizontal direction with ω; k z = ω 2 / c 2 - k x 2 Be the circular wavenumber of vertical direction, c is the speed of seawater;
(3) by formula M n ( k x , z = 0 , ω ) = P n ( k x , z = 0 , ω ) e - 2 i k z | H | Obtain the sea-bottom multiple that frequency-wavenumber domain is surveyed, wherein, M n(k x, z=0 is the multiple reflection of n big gun collection of frequency-wavenumber domain prediction ω), H is the mean value of minimum and maximum submarine elevation;
(4) to M n(k x, z=0 ω) carries out two-dimentional Fu Shi inverse transformation, obtains the multiple reflection m of n big gun collection of time domain prediction n(x, z=0 t), see Fig. 4.
4, when the difference of elevation of the most shallow and bottom, deep-sea during, as shown in Figure 2, adopt wave field to be bumped into the multiple reflection in method prediction fluctuating seabed according to the following steps less than BEI.
(1) to p n(x, z=0 t) carry out fourier transform about x and t, obtain P n(k x, z=0, ω), k wherein xWith ω be circular wavenumber and the circular frequency that is respectively horizontal direction; k z = ω 2 / c 2 - k x 2 Be the circular wavenumber of vertical direction, c is the speed of seawater;
(2) by formula P n ( k x , z = Z 0 , ω ) = P n ( k x , z = 0 , ω ) e - i k z | Z 0 | Obtain the seismic wave field of the frequency-wavenumber domain of (z=Z0) on bottom surface, shallow sea;
(3) from bottom surface, shallow sea (i=1) to bottom surface, deep-sea (i=N+1), by formula P n(k x, z=z I+1, ω)=P n(k x, z=z i, ω) e -2ik Δ zCalculate z i=z 0+ (i-1) the descending wave-wave field on the Δ z face is to P n(k x, z=z I+1, ω) carry out two-dimentional Fu Shi inverse transformation, and the following traveling-wave field q at place, record seabed n(z I+1, t);
(4) from bottom surface, deep-sea (i=N+1) to bottom surface, shallow sea (i=1), by formula p n(x, z=z i, t)=p n(x, z=z i, t)+q (z i, t) carry out wave field and be bumped into, by formula P n(k x, z=z I-1, ω)=P n(k x, z=z i, ω) e -ik Δ zCarry out the upstream wave field continuation;
(5) by formula M n ( k x , z = 0 , ω ) = P n ( k x , z = Z 0 , ω ) e - i k z | Z 0 | Obtain the sea-bottom multiple of the frequency-wavenumber domain that on the sea level, writes down;
(6) to M n(k x, z=0 ω) carries out two-dimentional Fu Shi inverse transformation, obtains the multiple reflection of n big gun collection of time domain prediction, m n(x, z=0 t), see Fig. 4.

Claims (1)

1. a method for identifying two-dimensional up-and-down sea-bottom multiple, realization as follows:
1) geological data, field survey or the indoor submarine elevation data of asking for of employing marine streamer seismic acquisition;
2) from minimum big gun gunfire collection, according to the submarine elevation that arrange at big gun collection place, obtain its maximum submarine elevation and minimum submarine elevation, it is poor to calculate minimum and maximum submarine elevation;
3) when minimum and maximum submarine elevation difference during, by horizontal seabed prediction sea-bottom multiple less than 5 meters:
p n(x, z=0 t) are the n big gun geological data that is positioned at z=0 on the sea level,
m n(x, z=0 t) are the sea-bottom multiple of determining of the correspondence n big gun that is positioned at z=0 on the sea level obtained,
Wherein n is a big gun number, and x is the distance of shot point to acceptance point, and z is the elevation coordinate, and t is a time coordinate;
The multiple reflection in described definite horizontal seabed is meant:
(1) to p n(x, z=0 t) carry out two-dimentional fourier transform about x and t, obtain P n(k x, z=0, ω), k wherein xBe respectively horizontal direction circular wavenumber and circular frequency with ω; Be the vertical direction circular wavenumber, c is the speed of seawater;
(2) by formula
Figure FSB00000376291600012
Obtain the multiple reflection of frequency-wavenumber domain prediction, wherein, M n(k x, z=0 is the multiple reflection of the n big gun of frequency-wavenumber domain prediction ω), H is the mean value of minimum and maximum submarine elevation;
(3) to M n(k x, z=0 ω) carries out two-dimentional Fu Shi inverse transformation, obtains the multiple reflection m of the n big gun of time domain prediction n(x, z=0, t);
4) when minimum and maximum submarine elevation difference during, by fluctuating seabed prediction sea-bottom multiple greater than 5 meters:
p n(x, z=0 t) are the n big gun geological data that is positioned on the z=0 of sea level,
p n(x, z=Z 0, t) be n big gun z=Z on bottom surface, shallow sea 0Seismic wave field,
p n(x, z=Z N, t) be n big gun z=Z on bottom surface, deep-sea NSeismic wave field,
p n(x, z i, t) be the n big gun at z i=z 0+ (i-1) Δ z i=1,2 ... seismic wave field on the N+1 face,
q n(z i, be that the n big gun is at z t) i=z 0+ (i-1) Δ z i=1,2 ... the record of bottom, N+1 face Shanghai,
m n(x, z=0 t) are the sea-bottom multiple of the prediction of the n big gun on the z=0 of sea level obtained,
Wherein n is a big gun number, and x is the distance of shot point to acceptance point, and z is the elevation coordinate, and t is a time coordinate, and Δ z is the step size of depth direction;
The multiple reflection in described definite fluctuating seabed is meant:
(1) to p n(x, z=0 t) carry out fourier transform about x and t, obtain P n(k x, z=0, ω), k wherein xWith ω for being respectively horizontal direction circular wavenumber and circular frequency; Be the vertical direction circular wavenumber, c is the speed of seawater;
(2) by formula
Figure FSB00000376291600022
Obtain z=Z on bottom surface, shallow sea 0The seismic wave field of frequency-wavenumber domain;
(3) from bottom surface, shallow sea i=1 to bottom surface, deep-sea i=N+1, by formula P n(k x, z=z I+1, ω)=P n(k x, z=z i, ω) e -2ik Δ zCalculate z i=z 0+ (i-1) the descending wave-wave field on the Δ z face is to P n(k x, z=z I+1, ω) carry out two-dimentional Fu Shi inverse transformation, and the following traveling-wave field q at place, record seabed n(z I+1, t);
(4) from bottom surface, deep-sea i=N+1 to bottom surface, shallow sea i=1, by formula p n(x, z=z i, t)=p n(x, z=z i, t)+q (z i, t) carry out wave field and be bumped into, by formula P n(k x, z=z I-1, ω)=P n(k x, z=z i, ω) e -ik Δ zCarry out the upstream wave field continuation;
(5) by formula
Figure FSB00000376291600031
Obtain the sea-bottom multiple of the frequency-wavenumber domain that on the sea level, writes down;
(6) to M n(k x, z=0 ω) carries out two-dimentional Fu Shi inverse transformation, obtains the multiple reflection m of the n big gun of time domain prediction n(x, z=0, t).
CN2008101141340A 2008-05-30 2008-05-30 Method for identifying two-dimensional up-and-down sea-bottom multiple Active CN101592738B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2008101141340A CN101592738B (en) 2008-05-30 2008-05-30 Method for identifying two-dimensional up-and-down sea-bottom multiple

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2008101141340A CN101592738B (en) 2008-05-30 2008-05-30 Method for identifying two-dimensional up-and-down sea-bottom multiple

Publications (2)

Publication Number Publication Date
CN101592738A CN101592738A (en) 2009-12-02
CN101592738B true CN101592738B (en) 2011-04-20

Family

ID=41407499

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2008101141340A Active CN101592738B (en) 2008-05-30 2008-05-30 Method for identifying two-dimensional up-and-down sea-bottom multiple

Country Status (1)

Country Link
CN (1) CN101592738B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103576198B (en) * 2012-08-02 2016-04-06 中国石油天然气集团公司 A kind of two-dimentional method for marine seismic data Free Surface multiple reflection Forecasting Methodology
CN106324664B (en) * 2015-06-29 2018-04-13 中国石油化工股份有限公司 A kind of method and apparatus for carrying out surface-related multiple forward modeling
CN105738881A (en) * 2016-03-09 2016-07-06 中国石油大学(华东) Offshore Envisat satellite height finding waveform resetting method based on waveform classification
CN109507722B (en) * 2017-09-15 2020-11-13 中国石油化工股份有限公司 Model and dual-wavefield continuation-based interlayer multiple prediction method and system
CN108957522B (en) * 2018-09-26 2019-03-19 中国海洋大学 A kind of submarine seismograph data multiple wave prediction technique based on Kirchhoff migiation
CN112882094B (en) * 2021-02-25 2022-05-10 中国石油集团东方地球物理勘探有限责任公司 First-arrival wave acquisition method and device, computer equipment and storage medium

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1797037A (en) * 2004-12-29 2006-07-05 中国石油天然气集团公司 Method for carrying out inversion for wave impedance of earthquake wave
CN1797041A (en) * 2004-12-29 2006-07-05 中国石油天然气集团公司 Method for eliminating linear and non-linear interference wave by using filtering operation at deep layer domain
CN1797031A (en) * 2004-12-29 2006-07-05 中国石油天然气集团公司 Method for shifting depth before superposition in combined earthquake data from multiple times of shooting

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1797037A (en) * 2004-12-29 2006-07-05 中国石油天然气集团公司 Method for carrying out inversion for wave impedance of earthquake wave
CN1797041A (en) * 2004-12-29 2006-07-05 中国石油天然气集团公司 Method for eliminating linear and non-linear interference wave by using filtering operation at deep layer domain
CN1797031A (en) * 2004-12-29 2006-07-05 中国石油天然气集团公司 Method for shifting depth before superposition in combined earthquake data from multiple times of shooting

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
耿伟峰等.波场延拓法压制海底多次波研究.《天然气工业》.2007,第27卷期刊第216-218页. *

Also Published As

Publication number Publication date
CN101592738A (en) 2009-12-02

Similar Documents

Publication Publication Date Title
Toomey et al. Tomographic imaging of the shallow crustal structure of the East Pacific Rise at 9° 30′ N
EP3387467B1 (en) Velocity model update with an inversion gradient
US10690792B2 (en) Amplitude-versus-angle analysis for quantitative interpretation
EP2778718B1 (en) Systems and methods for frequency-domain filtering and space-time domain discrimination of seismic data
CN101592738B (en) Method for identifying two-dimensional up-and-down sea-bottom multiple
EP2728385B1 (en) Systems and methods for high-resolution imaging using separated wavefields
CN103576198B (en) A kind of two-dimentional method for marine seismic data Free Surface multiple reflection Forecasting Methodology
CN101620276A (en) Method for attenuation of multiple reflections in seismic data
EP2755057A2 (en) High-fidelity adaptive curvelet domain primary-multiple separation processing of seismic data
US10877172B2 (en) Prediction and subtraction of multiple diffractions
EP2537049A2 (en) Estimating internal multiples in seismic data
EP3353577B1 (en) Determining node depth and water column transit velocity
CN104536045A (en) Ghost wave compression method based on wavelet processing
EP2762926B1 (en) Systems and methods for detecting swell noise in a seismic gather
CN104977615A (en) Model-statistics-pickup-based multiple suppression method of deep-sea OBC data
CN104730572A (en) Diffracted wave imaging method and device based on L0 semi-norm
US20150198729A1 (en) Regularization of spatially aliased seismic data
US9829593B2 (en) Determination of an impulse response at a subsurface image level
CN114994763B (en) Method for separating fluctuation components of continuous submarine seismic wave field in same frequency and same direction
US20160238727A1 (en) Removal of an Estimated Acquisition Effect from a Marine Survey Measurement
Ågesen Processing and interpretation of reflection seismic data from Isfjorden, Svalbard
CN114994763A (en) Same-frequency and same-direction continuous ocean bottom seismic wave field fluctuation component separation method
Wang et al. Link between crustal thickness and Moho transition zone at 9oN East Pacific Rise
Mordret et al. Near Surface Structures and Anisotropy from Cross-correlations of Ambient Seismic Noise at the Valhall Oil Field
Toomey et al. Tomographic imaging of the shallow crustal structure

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
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20091202

Assignee: NATIONAL ENGINEERING RESEARCH CENTER OF OIL AND GAS EXPLORATION SOFTWARE CO., LTD.

Assignor: Dongfang Geophysical Exploration Co., Ltd., China Petrochemical Corp.

Contract record no.: 2013990000780

Denomination of invention: Method for identifying two-dimensional up-and-down sea-bottom multiple

Granted publication date: 20110420

License type: Exclusive License

Record date: 20131121

LICC Enforcement, change and cancellation of record of contracts on the licence for exploitation of a patent or utility model