CN103675896B - A kind of diffracted wave and echo method for separate imaging - Google Patents
A kind of diffracted wave and echo method for separate imaging Download PDFInfo
- Publication number
- CN103675896B CN103675896B CN201210316656.5A CN201210316656A CN103675896B CN 103675896 B CN103675896 B CN 103675896B CN 201210316656 A CN201210316656 A CN 201210316656A CN 103675896 B CN103675896 B CN 103675896B
- Authority
- CN
- China
- Prior art keywords
- data
- dimensional
- wave
- road collection
- plane wave
- 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
Links
Landscapes
- Image Processing (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a kind of diffracted wave and echo method for separate imaging, belong to imaging of seismic data process field.The inventive method comprises the following steps: (1) carries out three-dimensional τ p direct transform and obtains plane wave data d (p each the road collection in each liney, px, y, x, τ)；(2) to the plane wave data d (p exported by step (1)y, px, y, x, τ) and use the destructing filtering of plane wave zone to process, obtain the plane wave data body d mainly containing diffracted wave energyr(py, px, y, x, τ)；(3) the main plane wave data body d containing diffracted wave energy obtained from step (2)r(py, px, y, x, τ) in extract three-dimensional τ p numeric field data corresponding to each road collection, Bing Congyuan road concentrates the big gun obtaining each road collection to examine coordinate, then described three-dimensional τ p numeric field data is done the prestack road collection after three-dimensional τ p inverse transformation obtains wave field separation；(4) the prestack road collection after the wave field separation exported by step (3) is carried out preSDM in shot domain and obtains diffracted wave separate imaging result.Present invention achieves diffracted wave signal extraction and little geological anomalous body or the blur-free imaging of tomography.
Description
Technical field
The invention belongs to imaging of seismic data process field, be specifically related to a kind of diffracted wave and be separated into echo
Image space method, for diffracted wave signal extraction and little geological anomalous body or tomography blur-free imaging.
Background technology
More weak owing to earthquake record coming the diffracted wave signal energy of geological anomalous body from childhood, be to its accurate one-tenth
Separate or diffracted wave extraction at imaging advance traveling-wave field as typically requiring.The destructing filter of current existing plane wave zone
Wave diffraction ripple and echo method for separate imaging are a kind of effective methods, and it comprises τ-p conversion, altogether p and puts down
The steps such as face wave datum sorting, plane wave zone destructing filtering, plane wave migration before stack.But plane wave to be made
Migration before stack is suitable with the imaging effect of tradition preSDM in shot domain, needs to decomposite when τ-p converts
Substantial amounts of plane wave data, this brings huge amount of calculation for skew below the most again, even more so that plane
The efficiency of ripple migration before stack is lower than the efficiency of tradition big gun territory skew on the contrary.To this end, can be by after wave field separation
Plane wave data regain prestack road collection by τ-p inverse transformation, then use tradition big gun territory migration before stack real
Existing seismic imaging.Information used by the conventional two dimension positive inverse transformation of τ-p is only the value of offset distance, therefore when twice
Offset distance value is close and time big gun inspection position difference is the biggest, twice geological data in the positive inverse transformation of τ-p by aliasing
Together, the problem that τ-p inverse transformation data signal distortion occurs, affect diffracted wave and echo separate imaging
Effect.
Summary of the invention
The present invention is directed to that diffracted wave signal in prestack migration image is weak, the unsharp problem of imaging, it is provided that a kind of
Diffracted wave and echo method for separate imaging, examine coordinate information according to the big gun in geological data and carry out three-dimensional time
The territory positive inverse transformation of τ-p, deconstructs plane wave filtered plane wave zone diffracted wave data and is intactly converted to fold
Front road collection, and then diffracted wave information is carried out big gun territory prestack migration image, thus avoid the two dimension positive contravariant of τ-p
Change distorted signals problem.
The present invention is achieved by the following technical solutions:
A kind of diffracted wave and echo method for separate imaging, comprise the following steps:
(1) each the road collection in each line is carried out three-dimensional τ-p direct transform and obtain plane wave data
d(py, px, y, x, τ)；
(2) to the plane wave data d (p exported by step (1)y, px, y, x, τ) and use the destructing filter of plane wave zone
Ripple processes, and obtains the plane wave data body d mainly containing diffracted wave energyr(py, px, y, x, τ)；
(3) the main plane wave data body d containing diffracted wave energy obtained from step (2)r(py, px, y, x, τ) in
Extracting the three-dimensional τ-p numeric field data that each road collection is corresponding, Bing Congyuan road concentrates the big gun obtaining each road collection to examine
Coordinate, then does the prestack road after three-dimensional τ-p inverse transformation obtains wave field separation to described three-dimensional τ-p numeric field data
Collection；
(4) the prestack road collection after the wave field separation exported by step (3) is carried out preSDM in shot domain
Obtain diffracted wave separate imaging result (i.e. separated echo and only included the imaging results of diffracted wave).
Described step (1) is specific as follows:
(11) use following three-dimensional τ-p direct transform formula that each the road collection in each line is carried out three
Dimension τ-p direct transform obtains τ-p numeric field data ψ:
Wherein, φiRepresent the i-th track data, hix、hiyIt is the offset distance component of the both direction of the i-th track data,
These two offset distance components are that to be examined coordinate by the big gun of the i-th track data calculated, px、pyRepresent two sides
To ray parameter；
(12) τ-p numeric field data ψ step (11) obtained forms plane wave number in being stored in five dimension data bodies
According to d (py, px, y, x, τ) and (a three-dimensional work area includes that N bar line, every line include M road collection, ψ Shi Yige road
The transformation results of collection, all ψ just constitute five dimension data bodies), wherein, y represents wire size, and x represents Xian Nei road
Collection number.
The plane wave data body d obtained from step (2) described in described step (3)r(py, px, y, x, τ) in take out
Take three-dimensional τ-p numeric field data corresponding to each road collection to be achieved in that by wire size y and road collection x from plane
Wave datum body dr(py, px, y, x, τ) and value obtains the three-dimensional τ-p numeric field data that each road collection is corresponding.
Described three-dimensional τ-p numeric field data is done by described step (3) after three-dimensional τ-p inverse transformation obtains wave field separation
Prestack road collection be achieved in that
After using following three-dimensional τ-p reconstructed formula that described three-dimensional τ-p numeric field data is transformed to wave field separation
Prestack road collection φi(t):
ψ′(py, px, τ) and=ψ (py, px, τ) and * ρ (τ)
Wherein, ρ (τ) is rho filter.
Compared with prior art, the invention has the beneficial effects as follows: the present invention is at the wave field of diffracted wave Yu echo
In separation process, the three-dimensional positive inverse transform method of τ-p is used to achieve prestack road collection without distortions to τ-p territory number
The data reconstruction of prestack road collection is arrived again according to body.By the present invention can effectively realize diffracted wave and echo point
From so that the diffracted wave image formation of little geological anomalous body or tomography becomes apparent from.The present invention may be used for around
Ejected wave signal extraction and little geological anomalous body or the blur-free imaging of tomography.
Accompanying drawing explanation
Fig. 1 is the step block diagram of diffracted wave of the present invention and echo method for separate imaging.
Fig. 2 is the seismic imaging result using the inventive method to obtain in embodiment 1.
Fig. 3 is the seismic imaging result using prior art to obtain in embodiment 1.
Fig. 4-1 is the former prestack road collection in embodiment 2.
Fig. 4-2 is the road collection in embodiment 2 through the two dimension positive inverse transformation of τ-p.
Fig. 4-3 is the road collection in embodiment 2 through the three-dimensional positive inverse transformation of τ-p.
Detailed description of the invention
Below in conjunction with the accompanying drawings the present invention is described in further detail: one diffracted wave of the present invention and reflection wavelength-division
From formation method, as it is shown in figure 1, comprise the following steps:
(1) each the road collection (i.e. prestack road collection in Fig. 1) in each line is carried out three-dimensional time territory
τ-p direct transform:
Owing to seismic data process having some lines, and in every line, include several road collection, therefore adopt
With two stage cycle, outer circulation is line circulation, i.e. y circulates, and interior circulation is the collection circulation of Xian Nei road, i.e. x circulates,
Two stage cycle is nested circulation.
Each road collection (i.e. prestack road collection in Fig. 1) in each line is carried out three-dimensional τ-p direct transform
Specific as follows:
The three-dimensional time territory τ-p direct transform of one road collection data (road collection data will be transformed to three-dimensional τ-p territory
Data) formula is:
Wherein, φiRepresent the i-th track data, hix、hiyIt is the offset distance component of the both direction of the i-th track data,
Namely in both direction, geophone station, to the displacement of shot point, is to be examined coordinate by the big gun of the i-th track data to be calculated
, px、pyRepresent the ray parameter of both direction.By upper formula by the different track data φ in road collectioniCarry out
Superposition obtains referring to three-dimensional τ-p numeric field data ψ；
The process that realizes of the three-dimensional τ-p direct transform of one road collection data includes four circulations, and outermost layer is pyFollow
Ring, the second layer is pxCirculation, third layer is circulation, and innermost layer is τ circulation, and these 4 circulations are belonging to often
Within individual road collection；
The result obtained after three-dimensional τ-p direct transform is stored in five dimension data body d (py, px, y, x, τ) in, wherein, y
Representing wire size, x represents Xian Nei road collection number.Five dimension data body d (py, px, y, x, τ) and it is described plane wave data, its
In contain the wavefield data of multiple plane wave, any one group of (py, px) corresponding 3D data volume represents one
The wavefield data of plane wave data；
(2) to the plane wave data (i.e. d (p exported by step (1)y, px, y, x, τ)) use plane wave zone
Destructing filtering (refer to, M.Turhan Taner, Rock Solid Images, Sergey Fomel .Separation
and imaging of seismic diffractions using plane-wave decomposition[J].Expanded
Abstracts of 76hAnnual International Meeting, SEG.2006) process, it is achieved plane
Wave zone diffracted wave separates with echo, it is thus achieved that the main plane wave data body containing diffracted wave energy
dr(py, px, y, x, τ)；
(3) two stage cycle (circulating for line y outward, Nei Weixiannei road collection x circulates), from five dimensions
According to body dr(py, px, y, x, τ) and middle extraction each road collection (φiRepresent the i-th track data, multiple φiForm a road collection)
Corresponding three-dimensional τ-p numeric field data (i.e. presses wire size y and road collection x from dr(py, px, y, x, τ) and value), and from former
Road concentrates the big gun obtaining each road collection to examine coordinate information, then three-dimensional τ-p numeric field data is three-dimensional τ-p anti-
Conversion, a three-dimensional τ-p numeric field data will be transformed to road collection data, and the formula of employing is:
ψ′(py, px, τ) and=ψ (py, px, τ) and * ρ (τ)
Wherein, the same rho filter used during ρ (τ) is conventional two dimension τ-p inverse transformation.By upper formula by three-dimensional
τ-p numeric field data is overlapped i.e. obtaining the prestack road collection after diffracted wave separates.
Above-mentioned three-dimensional τ-p inverse transformation realizes process and includes four circulations equally, and outermost layer is circulation, and second
Layer is pyCirculation, third layer is pxCirculation, innermost layer is t circulation, and these 4 circulations are belonging in each road collection
Portion；
(4) the prestack road collection after the wave field separation exported by step (3) is carried out preSDM in shot domain
(refer to, RISTOW D, R HL T.Fourier finite-difference migration
[J] .Geophysics, 1994,59 (12): 1881-1893), it is thus achieved that diffracted wave separate imaging result.
Although it is to be appreciated that foregoing invention content is towards three dimensional seismic data, but two dimension data is had equally
There is compatibility.When seismic data is two dimension, pyOnly one of which, its value is equal to 0 μ s/m.
The effect of the inventive method is described by 2 embodiments below:
Embodiment 1: sigsbee model data is carried out diffracted wave separate imaging process.Fig. 2 is to use this
The diffracted wave imaging results that bright method obtains, Fig. 3 is the conventional migration imaging result not doing wave field separation.
Comparison diagram 2 and Fig. 3 is it can be seen that result of implementation of the present invention makes Diffraction Point (as the rectangle frame in figure marks
Shown in) and the imaging of tomography (as shown in the oval mark in figure) become apparent from.
Embodiment 2: with two dimension τ-p just utilizing western certain real data contrast three-dimensional time territory positive inverse transformation of τ-p
The effect of inverse transformation.Fig. 4-1 is former prestack road collection, and Fig. 4-2 is the road collection through the two dimension positive inverse transformation of τ-p,
Fig. 4-3 is to utilize the inventive method through the road collection of the three-dimensional positive inverse transformation of τ-p.Contrast three figures visible, two dimension
There is short flat axle phenomenon in the positive and negative transformation results of τ-p, and Fig. 4-2 and Fig. 4-1 difference are very big, represent this positive contravariant
Data distortion problem is there is in changing；After the three-dimensional positive inverse transformation of τ-p, Fig. 4-3 and Fig. 4-1 is basically identical, table
Show and well recovered by this positive and negative transform data.Therefore, the three-dimensional time territory τ-p in the present invention is just
Inverse transform method is better than in the two dimension positive inverse transform method of τ-p.
Technique scheme is one embodiment of the present invention, for those skilled in the art,
On the basis of the invention discloses application process and principle, it is easy to make various types of improvement or deformation,
Being not limited solely to the method described by the above-mentioned detailed description of the invention of the present invention, the most previously described mode is simply
Preferably, and not restrictive meaning.
Claims (3)
1. a diffracted wave and echo method for separate imaging, it is characterised in that: said method comprising the steps of:
(1) each the road collection in each line is carried out three-dimensional τ-p direct transform and obtain plane wave data d (py,px,y,x,τ)；
(2) to the plane wave data d (p exported by step (1)y,px, y, x, τ) and use the destructing filtering of plane wave zone to process, obtain the plane wave data body d mainly containing diffracted wave energyr(py,px,y,x,τ)；
(3) the main plane wave data body d containing diffracted wave energy obtained from step (2)r(py,px, y, x, τ) in extract three-dimensional τ-p numeric field data corresponding to each road collection, Bing Congyuan road concentrates the big gun obtaining each road collection to examine coordinate, then described three-dimensional τ-p numeric field data is done the prestack road collection after three-dimensional τ-p inverse transformation obtains wave field separation；
(4) the prestack road collection after the wave field separation exported by step (3) is carried out preSDM in shot domain and obtains diffracted wave separate imaging result；
Wherein, described step (1) is specific as follows:
(11) use following three-dimensional τ-p direct transform formula that each the road collection in each line carries out three-dimensional τ-p direct transform and obtain τ-p numeric field data ψ:
Wherein, φiRepresent the i-th track data, hix、hiyBeing the offset distance component of the both direction of the i-th track data, these two offset distance components are that to be examined coordinate by the big gun of the i-th track data calculated, px、pyRepresent the ray parameter of both direction；
(12) τ-p numeric field data ψ step (11) obtained forms plane wave data d (p in being stored in five dimension data bodiesy,px, y, x, τ), wherein, y represents wire size, and x represents Xian Nei road collection number.
Diffracted wave the most according to claim 1 and echo method for separate imaging, it is characterised in that: the plane wave data body d obtained from step (2) described in described step (3)r(py,px, y, x, τ) in extract each road and collect corresponding three-dimensional τ-p numeric field data and be achieved in that by wire size y and the road number of collection x from plane wave data body dr(py,px, y, x, τ) and value obtains the three-dimensional τ-p numeric field data that each road collection is corresponding.
Diffracted wave the most according to claim 2 and echo method for separate imaging, it is characterised in that: described three-dimensional τ-p numeric field data is made the prestack road collection after three-dimensional τ-p inverse transformation obtains wave field separation by described step (3) and is achieved in that
Prestack road collection φ after using following three-dimensional τ-p reconstructed formula that described three-dimensional τ-p numeric field data is transformed to wave field separationi(t):
ψ'(py,px, τ) and=ψ (py,px,τ)*ρ(τ)
Wherein, ρ (τ) is rho filter.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210316656.5A CN103675896B (en) | 2012-08-30 | 2012-08-30 | A kind of diffracted wave and echo method for separate imaging |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210316656.5A CN103675896B (en) | 2012-08-30 | 2012-08-30 | A kind of diffracted wave and echo method for separate imaging |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103675896A CN103675896A (en) | 2014-03-26 |
| CN103675896B true CN103675896B (en) | 2016-08-10 |
Family
ID=50314015
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210316656.5A Active CN103675896B (en) | 2012-08-30 | 2012-08-30 | A kind of diffracted wave and echo method for separate imaging |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103675896B (en) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105093288B (en) * | 2014-05-21 | 2017-09-19 | 中国石油化工股份有限公司 | A kind of diffracted wave separation method based on kinematics wave field attributes |
| 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 |
| CN108693559B (en) * | 2017-04-05 | 2020-04-07 | 中国石油化工股份有限公司 | Seismic wave combined imaging method and system |
| CN109116423B (en) * | 2018-07-11 | 2020-03-17 | 北京奥能恒业能源技术有限公司 | Method and device for suppressing diffracted multiples |
| CN109471172B (en) * | 2018-12-26 | 2020-04-03 | 中国科学院地球化学研究所 | Surface wave purification method and device based on same-phase axis morphological difference |
| CN111025383B (en) * | 2019-11-21 | 2021-09-24 | 徐州工程学院 | Method for qualitatively judging water filling condition of tunnel front karst cave based on diffracted transverse waves |
-
2012
- 2012-08-30 CN CN201210316656.5A patent/CN103675896B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN103675896A (en) | 2014-03-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103675896B (en) | A kind of diffracted wave and echo method for separate imaging | |
| CN102193107B (en) | Method for separating and denoising seismic wave field | |
| CN105510976B (en) | A kind of many subwaves combination adaptive attenuation method | |
| CN102778693B (en) | Diffracted wave separation processing method based on reflection wave layer leveling extraction and elimination | |
| CN106526677B (en) | A kind of wideband reverse-time migration imaging method of marine adaptive compacting ghost reflection | |
| CN104932010B (en) | A kind of diffracted wave separation method based on the sparse Radon transformation of shortcut fringing | |
| CN101382598B (en) | Pressing method for true 3-d seismics data linear noise | |
| CN108549106B (en) | Aliasing noise drawing method and device | |
| CN102053276A (en) | Two-dimensional filtering method for a plurality of gathers of digital seismic signal | |
| CN104024890B (en) | The wave field simultaneously gathered by weighted superposition suppresses 4D noise | |
| CN101893719A (en) | Method for separating and de-noising wave fields of longitudinal and transverse earthquake waves | |
| CN103869364A (en) | Multiple wave suppression method based on dual parabolic Radon transformation | |
| CN107678062A (en) | The integrated forecasting deconvolution of hyperbolic Radon domains and feedback loop methodology multiple suppression model building method | |
| CN107884828A (en) | It is a kind of in spatial frequency domain based on the theoretical terrible ripple drawing method of Green | |
| CN103399348A (en) | Denoising method for seismic signal based on Shearlet transform | |
| CN103675903A (en) | Multi-source seismic wave field separating method by utilizing random denoising technology | |
| CN104914466A (en) | Method for improving seismic data resolution | |
| CN105652322A (en) | T-f-k field polarization filtering method for multi-component seismic data | |
| CN104216010A (en) | Method for increasing quality of seismic data by using harmonic waves of controllable seismic focus | |
| CN103675902A (en) | Optimal direction edge monitoring method | |
| CN107884829A (en) | A kind of method for combining compacting shallow sea OBC Multiple Attenuation in Seismic Data | |
| CN103412325B (en) | The method of the wave field separation of a kind of earthquake P-wave And S and denoising | |
| CN105510975B (en) | Improve the method and device of geological data signal to noise ratio | |
| CN104570114B (en) | A kind of reverse-time migration Noise Elimination method based on wavefield decomposition | |
| CN105301648A (en) | Method of acquiring common reflection surface stacking parameters |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | 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 |