CN102455439B - Diffracted wave field separation method based on Kirchhoff integral method - Google Patents

Diffracted wave field separation method based on Kirchhoff integral method Download PDF

Info

Publication number
CN102455439B
CN102455439B CN 201010527864 CN201010527864A CN102455439B CN 102455439 B CN102455439 B CN 102455439B CN 201010527864 CN201010527864 CN 201010527864 CN 201010527864 A CN201010527864 A CN 201010527864A CN 102455439 B CN102455439 B CN 102455439B
Authority
CN
China
Prior art keywords
inclination angle
waves
kirchhoff
diffracted wave
step
Prior art date
Application number
CN 201010527864
Other languages
Chinese (zh)
Other versions
CN102455439A (en
Inventor
白英哲
孙赞东
Original Assignee
中国石油大学(北京)
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 中国石油大学(北京) filed Critical 中国石油大学(北京)
Priority to CN 201010527864 priority Critical patent/CN102455439B/en
Publication of CN102455439A publication Critical patent/CN102455439A/en
Application granted granted Critical
Publication of CN102455439B publication Critical patent/CN102455439B/en

Links

Abstract

The invention provides a diffracted wave field separation method based on a Kirchhoff integral method. In seismic exploration, subsurface information can be inferred through responses of seismic waves to subsurface media; when a subsurface medium catastrophe happens, diffracted waves usually are response to the subsurface medium catastrophe. Main oil gas storage spaces in carbonate rock regions (such as Tarim and Sichuan) are secondary corroded holes and cracks, the sizes of the holes and the cracks are smaller than the wavelength of the seismic waves; the seismic response of the holes and the cracks is mainly diffraction, therefore the research on the diffracted waves has significant theoretical and practical meanings. However, in common deviation, the diffracted waves are converged to be one point on a profile and are submerged in a plurality of events and cannot be recognized. In a common reflection point (CRP) trace gather, since diffraction points and reflection points are the flattening of the events, the diffraction points cannot be recognized in the trace gather. Therefore, the invention provides a method for outputting an inclined angle trace gather through Kirchhoff deviation and hole diffracted waves and interface reflection waves can be rapidly and efficiently separated through obvious difference between the reflection waves and the diffracted waves in the inclined angle trace gather.

Description

Diffracted wave field separation method based on Kirchhoff integral method

Technical field

The present invention relates to the exploration of oil and gas field technical field, belong to the seismic signal analysis category, specifically a kind of seam hole imaging technique based on Kirchhoff integral method.

Background technology

For complex dielectrics, such as the Tarim Basin, oil gas preserves in ancient Cambrian system and ORDOVICIAN CARBONATE stratum, the buried depth on these stratum is below 5000m, primary pore is compacted already, filling totally, reservoir space of today is the secondary pores that will form the day after tomorrow by the effect such as corrosion, fracture.The distribution of these holes is not subjected to the control of sedimentary facies belt, and the scale that distributes and inhomogeneous.Therefore, wanting to seek these hole seam type reservoirs can only be by the means of earthquake.

The size of corrosion hole is often much smaller than the wavelength of seismic event, can think the catastrophe point of lithology.Diffraction will occur in seismic event on the catastrophe point of lithology, as long as detect these diffracted waves, by the seismic migration technology, diffracted wave is restrained, and just can determine the position of corrosion hole.Yet the fact is so not simple, and on the section after the skew, diffracted wave is often converged a point, thereby is submerged between numerous reflection line-ups, is difficult for differentiating.In the CRP gather after common migration before stack (CRP), Diffraction Point and reflection spot show as again evening up of lineups equally, therefore, can not distinguish reflection wave and diffracted wave by the CRP gather (CRP) of migration before stack.

At present also be in the starting stage in the research aspect the diffracted wave field separation both at home and abroad, also do not have the achievement in research of moulding.On-the-spot in exploration, diffracted wave is still the difficult problem of processing.Under many circumstances, owing to there not being effective separation means, diffracted wave can only be suppressed as a kind of interference.But so, underground medium abrupt information has just been lost, and the resolution of seismic prospecting also descends thereupon.Therefore, petroleum prospecting industry is especially faced when having the complicated reservoirs of corrosion hole in the urgent need to a kind of effective diffracted wave field separation method.

Summary of the invention

The invention is intended to provide a kind of by the high s/n ratio seismic data; by kirchhoff skew output inclination angle Dao Ji; obtain the independent imaging results of diffracted wave according to the obvious difference of concentrating diffracted wave and reflection wave in the road, inclination angle; thereby underground geology catastrophe point is carried out Precise imaging, so that better identify Carbonate Reservoir seam reservoir.

The technical scheme that the interval velocity field realization above-mentioned purpose that the high s/n ratio seismic data that this patent use was processed and Accurate velocity analysis obtain is taked is as follows:

Step 1: use the interval velocity model of input to carry out ray tracing, obtain the cautious whilst on tour tabulation of each big gun of ground.The rapid finite difference ray casting that ray tracing can select Vidale (1988) to propose.

Step 2: Kirchhoff method output inclination angle Dao Ji:

At first, read in the input geological data one.And read in respectively the right whilst on tour tabulation of this together corresponding big gun inspection that step 1 obtains;

Secondly, calculate underground each point place's incident ray and reflected ray with the angle of vertical direction by following formula:

tan α s = ∂ τ s / ∂ x ∂ τ s / ∂ z

tan α r = ∂ τ r / ∂ x ∂ τ r / ∂ z

In the formula: τ sBe shot point whilst on tour, τ rBe geophone station whilst on tour, α sFor the angle α r of incident ray and vertical direction is the angle of reflected ray and vertical direction, (x, z) is volume coordinate.

According to geometric relationship, can obtain inclination angle theta, reflection angle With α s, α rBetween relation

θ = α s + α r 2

As can be known by ray tracing, can set up the relation of shot point coordinate and calculation level coordinate and inclination angle theta and geophone station coordinate r by above analysis

s=s(x,z,r,θ)

And then the kirchhoff offset equation can be write as

R ( x , z , θ ) = ∫ Γ u ( s ( x , z , r , θ ) , r , τ s + τ r ) dr

So just realized kirchhoff skew output inclination angle Dao Ji.

Step 3: concentrate in the road, inclination angle that step 2 obtains, reflection wave shows as the curve of a upward opening, and the diffracted wave of observing directly over Diffraction Point shows as a horizontal linear.Therefore, can both be separated by the method for medium filtering.Can carrying out according to following operation of medium filtering:

At first, read in an inclination angle Dao Ji.

Secondly, the amplitude that its adjacent several roads same depth point is got in every each depth point together of concentrating for choosing the road through consists of a sequence, and this sequence is pressed from small to large (from big to small) ordering.

At last, get the intermediate value of above-mentioned sequence as the filtering output of this this depth point, road.

Step 4: road, the inclination angle collection stack of the filtering output that step 3 is obtained just can generate the independent imaging section of diffracted wave.

Above embodiment only is used for explanation the present invention, but not is used for limiting the present invention.

Description of drawings

Horizontal reflection INTERFACE MODEL of Fig. 1 (a);

Fig. 1 (b) is that Fig. 1 (a) is at the inclination angle Dao Ji at the horizontal reflection interface of x=0 place observation;

Inclination angle of Fig. 1 (c) is 25 ° tiltedinterface model;

Fig. 1 (d) is that Fig. 1 (c) is at the inclination angle Dao Ji of the tiltedinterface of x=0 place observation.

The single Diffraction Point model of Fig. 2 (a);

Fig. 2 (b) is the inclination angle Dao Ji of Fig. 2 (a) institute representation model when observing at the x=0 place;

Fig. 2 (c) is the inclination angle Dao Ji of Fig. 2 (a) institute representation model when observing at the x=100m place;

Fig. 2 (d) is the inclination angle Dao Ji of Fig. 2 (a) institute representation model when observing at the x=-100m place;

Each parameter schematic diagram of Fig. 3 ray tracing system.

Fig. 4 multilayered medium diffraction model.5 Diffraction Points are wherein arranged, lay respectively at x=2000m and x=2200m place.

Near in Fig. 5 model shown in Figure 4 the Diffraction Point several inclination angle Dao Ji.The road of black line mark integrates the place as x=2000m.

Fig. 6 (a) is the inclination angle Dao Ji of black line mark among Fig. 5;

Fig. 6 (b) is the diffracted wave after road, inclination angle collection separates among Fig. 6 (a);

Fig. 6 (c) is the reflection wave after road, inclination angle collection separates among Fig. 6 (a).

Skew stacked section before Fig. 7 (a) separates model diffraction shown in Figure 4;

Skew stacked section after Fig. 7 (b) separates model diffraction shown in Figure 4.

Embodiment

By a numerical model working of an invention process is described:

A multilayered medium rate pattern of model (as shown in Figure 4).5 Diffraction Points are wherein arranged, lay respectively at x=2000m and x=2200m place.There are two Diffraction Points at the x=2000m place, and the degree of depth is respectively 3000m and 4400m; There are three Diffraction Points at the x=2200m place, and the degree of depth is respectively 3000m, 4400m and 4800m.

Adopt Kirchhoff just to drill to model, adopt split shooting, 300 roads receive, and totally 300 big guns cover for 150 times, 2 milliseconds of samplings, and record length is 4 seconds.

Then generate several inclination angle Dao Ji (as shown in Figure 5) by the described method of step 2.Concentrating reflecting interface in the road, inclination angle is the curve that is bent upwards, corresponding inclination angle and the degree of depth of its extreme point indication.Diffraction Point is concentrated in the road shown in the black line (x=2000m) and is well evened up, and concentrates on monotonous curve in other roads.

Separate diffracted wave and reflection wave road collection (as shown in Figure 6) afterwards by step 3 again.

Stacked section before and after at last diffracted wave field being separated compares (as shown in Figure 7).Can find out, on the stacked section (Fig. 7 (a)) before diffracted wave separates, the reflection of bed interface is very strong, has covered in the reflection of hole fully; Can't see the reflective information of bed interface on the stacked section of diffracted wave after separating (Fig. 7 (b)), the just reflective information (oval point among the figure) of hole that stays.

List of references

Vidale,J.E.,1988,Finite-difference?traveltime?calculation:Bull.,Seis.Sot.Am.,78,2062-2076.

Vidale,J.E.,1990,3-D?finite-difference?traveltime?calculation:Geophysics,55,523-528.

Claims (1)

1. diffracted wave field separation method based on Kirchhoff integral method, by the high s/n ratio seismic data, by kirchhoff skew output inclination angle Dao Ji, obtain the independent imaging results of diffracted wave according to the obvious difference of concentrating diffracted wave and reflection wave in the road, inclination angle, thereby underground geology catastrophe point is carried out Precise imaging, so that better identify Carbonate Reservoir seam reservoir;
Above-mentioned purpose is realized in the interval velocity field that the high s/n ratio seismic data that described method use was processed and Accurate velocity analysis obtain, and the technical scheme of taking is as follows:
Step 1: use the formation velocity model of input to carry out ray tracing, obtain the cautious whilst on tour tabulation of each big gun of ground; The rapid finite difference ray casting that ray tracing selects Vadile (1988) to propose;
Step 2: Kirchhoff integral method output inclination angle Dao Ji; Described Kirchhoff integral method output road, inclination angle collection is:
At first, read in the input geological data one; And read in respectively the right whilst on tour tabulation of this together corresponding big gun inspection that step 1 obtains;
Secondly, calculate underground each point place's incident ray and reflected ray with the angle of vertical direction by following formula:
In the formula: τ sBe shot point whilst on tour, τ rBe geophone station whilst on tour, α sAngle α for incident ray and vertical direction rBe the angle of reflected ray and vertical direction, (x, z) is volume coordinate;
According to geometric relationship, can obtain inclination angle theta, reflection angle φ and α s, α rBetween relation
As can be known by ray tracing, can set up the relation of shot point coordinate and calculation level coordinate and inclination angle theta and geophone station coordinate r by above analysis
And then the kirchhoff offset equation can be write as
So just realized Kirchhoff skew output inclination angle Dao Ji;
Step 3: concentrate in the road, inclination angle that step 2 obtains, reflection wave shows as the curve of a upward opening, and the diffracted wave of observing directly over Diffraction Point shows as a horizontal linear; Therefore, can both be separated by the method for medium filtering;
Step 4: road, the inclination angle collection stack of the filtering output that step 3 is obtained just can generate the independent imaging section of diffracted wave.
CN 201010527864 2010-11-02 2010-11-02 Diffracted wave field separation method based on Kirchhoff integral method CN102455439B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201010527864 CN102455439B (en) 2010-11-02 2010-11-02 Diffracted wave field separation method based on Kirchhoff integral method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201010527864 CN102455439B (en) 2010-11-02 2010-11-02 Diffracted wave field separation method based on Kirchhoff integral method

Publications (2)

Publication Number Publication Date
CN102455439A CN102455439A (en) 2012-05-16
CN102455439B true CN102455439B (en) 2013-10-23

Family

ID=46038845

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201010527864 CN102455439B (en) 2010-11-02 2010-11-02 Diffracted wave field separation method based on Kirchhoff integral method

Country Status (1)

Country Link
CN (1) CN102455439B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104237940A (en) * 2014-09-29 2014-12-24 中国石油天然气股份有限公司 Diffracted wave imaging method and diffracted wave imaging device based on dynamical features

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102749647B (en) * 2012-07-12 2014-09-17 恒泰艾普石油天然气技术服务股份有限公司 Quantitative method for space of holes in reservoir
CN102854529B (en) * 2012-07-13 2013-10-23 孙赞东 Reflected wave generalized radon spectrum method diffracted wave field separation technology
CN103675897B (en) * 2012-08-30 2017-10-27 中国石油化工股份有限公司 A kind of earthquake diffracted wave method for separate imaging
US9134444B2 (en) * 2012-11-27 2015-09-15 Chevron U.S.A., Inc. System and method for deducing cavern properties
CN103399345B (en) * 2013-08-16 2016-08-24 中国科学院地质与地球物理研究所 The investigation method of a kind of buried hill fissure distribution and device
CN104076395B (en) * 2014-04-16 2016-01-13 孙学凯 Based on minute surface Energy extraction and the formation method of filtering combination
CN103941289B (en) * 2014-05-09 2017-01-25 中国科学院地质与地球物理研究所 Diffraction point imaging method and equipment
CN105093288B (en) * 2014-05-21 2017-09-19 中国石油化工股份有限公司 A kind of diffracted wave separation method based on kinematics wave field attributes
CN104730571A (en) * 2015-03-11 2015-06-24 中国科学院地质与地球物理研究所 Method and device for identifying small-scale geologic body through diffraction refocusing
CN104730572B (en) * 2015-03-11 2016-11-30 中国科学院地质与地球物理研究所 A kind of diffracted wave formation method based on L0 semi-norm and device
CN105785439B (en) * 2016-02-01 2018-05-18 北京中科联衡科技有限公司 The Forecasting Methodology and device of small scale heterogeneous geologic body spatial distribution position
CN106772583B (en) * 2017-01-10 2018-09-04 中国科学院地质与地球物理研究所 A kind of earthquake diffracted wave separation method and device
CN107255832B (en) * 2017-07-10 2019-02-26 中国石油大学(华东) A kind of inversion method of subsurface structure
CN107942374B (en) * 2017-11-15 2019-03-15 中国矿业大学(北京) Diffracted wave field extracting method and device
CN108845354B (en) * 2018-09-26 2019-07-30 西安石油大学 A kind of method of intermediate value resistance filtering separation earthquake diffracted wave
CN109507732A (en) * 2018-11-07 2019-03-22 中石化石油工程技术服务有限公司 A kind of diffracted wave method for separate imaging based on imaging trace gather
CN110118993A (en) * 2019-05-22 2019-08-13 中国矿业大学(北京) Diffraction wave imaging method and device

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101105537A (en) * 2006-07-12 2008-01-16 中国石油集团东方地球物理勘探有限责任公司 High accuracy depth domain prestack earthquake data inversion method

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6904368B2 (en) * 2002-11-12 2005-06-07 Landmark Graphics Corporation Seismic analysis using post-imaging seismic anisotropy corrections

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101105537A (en) * 2006-07-12 2008-01-16 中国石油集团东方地球物理勘探有限责任公司 High accuracy depth domain prestack earthquake data inversion method

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Evgeny Landa, Sergey Fomel, Moshe Reshef.Separation, imaging, and velocity analysis of seismic diffractions using migrated dip-angle gathers.《SEG Las Vegas 2008 Annual Meeting》.2008,2176-2179.
Moshe Reshef,Evgeny Landa.Post-stack velocity analysis in the dip-angle domain using diffractions.《Geophysical Prospecting》.2009,811-821. *
Separation, imaging, and velocity analysis of seismic diffractions using migrated dip-angle gathers;Evgeny Landa, Sergey Fomel, Moshe Reshef;《SEG Las Vegas 2008 Annual Meeting》;20081231;2176-2179 *
钢筋混凝土缺陷的探地雷达检测模拟与成像效果;黄玲等;《物探与化探》;20070430;第31卷(第02期);181-185 *
黄玲等.钢筋混凝土缺陷的探地雷达检测模拟与成像效果.《物探与化探》.2007,第31卷(第02期),181-185.

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104237940A (en) * 2014-09-29 2014-12-24 中国石油天然气股份有限公司 Diffracted wave imaging method and diffracted wave imaging device based on dynamical features
CN104237940B (en) * 2014-09-29 2017-03-08 中国石油天然气股份有限公司 A kind of diffraction wave imaging method based on dynamic characteristic and device

Also Published As

Publication number Publication date
CN102455439A (en) 2012-05-16

Similar Documents

Publication Publication Date Title
Sheriff et al. Exploration seismology
Juhlin et al. 3D baseline seismics at Ketzin, Germany: the CO 2 SINK project
US5671136A (en) Process for seismic imaging measurement and evaluation of three-dimensional subterranean common-impedance objects
Griffiths et al. Applied geophysics for geologists and engineers: the elements of geophysical prospecting
Van der Veen et al. Design and application of a towed land-streamer system for cost-effective 2-D and pseudo-3-D shallow seismic data acquisition
CN100487488C (en) Zero hypocentral distance vertical seismic section compressional-shear wave data depth field corridor stacked section processing method
Kearey et al. An introduction to geophysical exploration
CA2557526C (en) Prediction of shallow drilling hazards using seismic refraction data
US6092025A (en) Hydrocarbon edge detection using seismic amplitude
CN104335072A (en) Integrated passive and active seismic surveying using multiple arrays
Tsuji et al. Methane-hydrate occurrence and distribution in the eastern Nankai Trough, Japan: Findings of the Tokai-oki to Kumano-nada methane-hydrate drilling program
GB2462222A (en) Method for processing acoustic reflections in arrray data to image near borehole geological structure
CN101907725B (en) Crack prediction method and device
Daley et al. Field testing of modular borehole monitoring with simultaneous distributed acoustic sensing and geophone vertical seismic profiles at Citronelle, Alabama
Malehmir et al. 3D seismic reflection imaging of volcanic-hosted massive sulfide deposits: Insights from reprocessing Halfmile Lake data, New Brunswick, Canada
US6055482A (en) Method of seismic signal processing
EP2609450B1 (en) System and method for acquisition and processing of elastic wavefield seismic data
CN102012521B (en) Method for detecting pre-stack cracks in seismic reservoir prediction
CN102466816B (en) Inversion method for stratum elasticity constant parameter of pre-stack seismic data
EP2189818B1 (en) Method of wavefield extrapolation for single-station, dual-sensor towed streamer signals
Pugin et al. Hydrogeological prospecting using P‐and S‐wave landstreamer seismic reflection methods
CN102455439A (en) Diffracted wave field separation method based on Kirchhoff integral method
CA2387760A1 (en) Transfer function method of seismic signal processing and exploration
Urosevic et al. Targeting nickel sulfide deposits from 3D seismicreflection data at Kambalda, Australia
Yordkayhun et al. 3D seismic traveltime tomography imaging of the shallow subsurface at the CO2SINK project site, Ketzin, GermanyCO2SINK 3D seismic traveltime tomography

Legal Events

Date Code Title Description
DD01 Delivery of document by public notice

Addressee: Bai Yingzhe

Document name: Notification to Make Rectification

PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination
GR01 Patent grant
C14 Grant of patent or utility model
CB03 Change of inventor or designer information

Inventor after: Sun Zandong

Inventor after: Liu Zhishui

Inventor after: Bai Yingzhe

Inventor after: Wang Peifu

Inventor after: Liu Lifeng

Inventor after: Tian Jun

Inventor after: Wang Qinghua

Inventor after: Yang Haijun

Inventor after: Liang Xianghao

Inventor after: Wang Zhenbiao

Inventor after: Han Jianfa

Inventor after: Zhang Dong

Inventor before: Bai Yingzhe

Inventor before: Sun Zandong

COR Change of bibliographic data