CN106501859B - A kind of moving integration lineups of prestack trace gather even up method - Google Patents
A kind of moving integration lineups of prestack trace gather even up method Download PDFInfo
- Publication number
- CN106501859B CN106501859B CN201611238674.0A CN201611238674A CN106501859B CN 106501859 B CN106501859 B CN 106501859B CN 201611238674 A CN201611238674 A CN 201611238674A CN 106501859 B CN106501859 B CN 106501859B
- Authority
- CN
- China
- Prior art keywords
- trace gather
- lineups
- evening
- amount
- evened
- 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
- 238000000034 method Methods 0.000 title claims abstract description 20
- 230000010354 integration Effects 0.000 title claims abstract description 16
- 230000008569 process Effects 0.000 claims description 3
- 238000012937 correction Methods 0.000 description 5
- 230000006872 improvement Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Traffic Control Systems (AREA)
Abstract
The invention discloses a kind of moving integration lineups of prestack trace gather to even up method, including：Calculate the amount of the evening up DT1 of neighboring track in original trace gather；The amount of evening up DT1 is applied in original trace gather, a part is generated and evens up superposition trace gather；The amount of the evening up DT2 that superposition trace gather calculates neighboring track again is evened up using part；Shortcut " superposition road " is selected, moving integration is implemented to the amount of evening up DT2 and is evened up, entire trace gather is evened up.The present invention can solve the problems, such as the lineups distortion changing of the relative positions caused by conventional Similar operator lineups are evened up and solve two class AVO of lineups reversed polarity.
Description
Technical field
The invention belongs to petroleum exploration field, it is related to a kind of moving integration lineups of prestack trace gather and evens up method.
Background technology
In seismic data process, overlap-add procedure is a kind of effective denoising method.But in actual seismic Data processing,
Since NMO velocity is the NMO algorithms assumed based on flat bed, and it is the interpolation speed of discrete point, it is difficult to accurately calculate；This
Outside, when there are when heterogeneous body or even pre-stack time migration algorithm cannot all obtain the trace gather evened up in seismic raypaths.Cause
This can all have certain residual normal moveout in the data after moving correction, and stacked section can be reduced when doing in-phase stacking
Resolution ratio.
Either the common reflection point superposition in seism processing stage or the AVA and AVAz in seismic data interpretation stage
Analysis, be all by lineups even up after CRP gather based on, lineups injustice can influence the quality of follow-up work,
Even obtain the result of mistake.In processing, lineups are often carried out according to the traveltime-distance equation of CRP gather and even up (dynamic school
Just), although dynamic correction can ensure that CRP gather lineups are evened up substantially, by various factors (such as due to earth's surface elevation
Static correction caused by rising and falling acutely is inaccurate, due to conventional time curve caused by horizontal layer isotropic medium assumed condition
Error in equation) influence, dynamic correction often cannot achieve CRP gather and evens up, and then, reduce common reflection point superposition
Signal-to-noise ratio, and lead to the inaccurate even mistake of elastic parameter analyzed by AVA and AVAz.
It is to obtain the time shift amount of accurate prestack trace gather to carry out the premise that trace gather is evened up, and there has been proposed much about residue
Dynamic bearing calibration, be such as concerned with earlier and high-precision superimposing technique, radial filtering technique, using residual normal moveout as it is local with
The technology etc. that residual quantity is handled when machine static state；Some researchers also proposed solves residual NMO correction using residual velocity analysis
The problem of.Such as by extracting the time-shifting method of sampled point time difference different moments.It passes through to model trace and CMP trace gather subregions
Domain, window carries out cross-correlation calculation when opening, and accurately determines the notable extreme point time shift amount of seismic channel, is stretched to other data
Or it shrinks and moves integrally again.
Mainstream trace gather evens up algorithm there are following some problems at present：
A library track is determined in advance 1. evening up algorithm and generally requiring, then carries out in other roads in trace gather with library track
Comparison, when two seismic channel lineups changing of the relative positions of comparison larger (more than one wavelength), frequently can lead to mistake evens up knot
Fruit, or trace gather can not be evened up.
It evens up algorithm 2. most of and can't solve the problems, such as two class AVO of lineups reversed polarity well.
3. it is poor that general trace gather evens up algorithm lateral stability, or existed and even up phenomenon.
Invention content
Shortcoming present in regarding to the issue above, the present invention provide a kind of moving integration lineups drawing of prestack trace gather
Quadratic method.
To achieve the above object, the present invention provides a kind of moving integration lineups of prestack trace gather and evens up method, including：
Step 1, the amount of the evening up DT1 for calculating neighboring track in original trace gather；
The amount of evening up DT1 is applied in original trace gather by step 2, is generated a part and is evened up superposition trace gather；
Step 3 evens up the amount of the evening up DT2 that superposition trace gather calculates neighboring track again using part；
Step 4, selection shortcut " superposition road ", implement moving integration to the amount of evening up DT2 and even up, even up entire trace gather.
As a further improvement on the present invention, it in the calculating process of step 1, is counted come when controlling using the zero passage of waveform
The length of window.
As a further improvement on the present invention, the step 2 includes：
Step 21 laterally chooses sliding window in original trace gather；
Step 22, using the interim orbit of window as library track, trace gather in window is evened up and is superimposed；
Step 23, stack result are placed on the interim orbit of window, are obtained a part and are evened up superposition trace gather.
Compared with prior art, beneficial effects of the present invention are：
The present invention provides a kind of moving integration lineups of prestack trace gather and evens up method, can solve conventional Similar operator
The lineups distortion changing of the relative positions caused by lineups are evened up and the two class AVO for solving the problems, such as lineups reversed polarity.
Description of the drawings
Fig. 1 is that the moving integration lineups of prestack trace gather disclosed in an embodiment of the present invention even up the flow chart of method；
Fig. 2 is that the present invention is based on the long selection method graphic software platform form figures of the adaptive waveform associated window of zero crossing；
Fig. 3 is that superposition trace gather generating principle graphic software platform form figure is evened up in present invention part；
Fig. 4 is the basic principle graphic software platform form figure of the amount of evening up integral of the invention.
Specific implementation mode
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is
A part of the embodiment of the present invention, instead of all the embodiments.Based on the embodiments of the present invention, ordinary skill people
The every other embodiment that member is obtained without making creative work, shall fall within the protection scope of the present invention.
The present invention is described in further detail below in conjunction with the accompanying drawings：
As shown in Figure 1, the present invention, which provides a kind of moving integration lineups of prestack trace gather, evens up method, including：
Step 1, the amount of the evening up DT1 for calculating neighboring track in original trace gather；
The similitude of neighboring track waveform is compared in step 1.Waveform how is selected also to determine the amount of evening up DT1 result of calculations,
In order to solve the frequency difference due to seismic channel depth layer to trouble caused by window when determining comparison of wave shape, the present invention uses waveform
Zero passage count come the length of window when controlling, can ensure that the waveform of different frequency can obtain proper window in this way
It is long, as shown in Figure 2.
Step 2, in order to improve trace gather signal-to-noise ratio：The amount of evening up DT1 is applied in original trace gather, a part is generated and evens up
It is superimposed trace gather；
In step 2 in order to overcome the amount of evening up caused by trace gather signal-to-noise ratio is too low to calculate stability problem, propose
Superposition algorithm is evened up in part.As shown in figure 3, laterally choosing a small sliding window in original trace gather, the length of window can
It is selected according to actual demand；Using the interim orbit of window as library track, trace gather in window is evened up and is superimposed, superposition knot
Fruit is placed on the interim orbit of window, has thus obtained a part and has evened up superposition trace gather.
Step 3 evens up the amount of the evening up DT2 that superposition trace gather calculates neighboring track again using part.
Step 4, selection shortcut " superposition road ", implement moving integration to the amount of evening up DT2 and even up, even up entire trace gather；Wherein,
The superposition number of neighboring track can be selected according to actual demand.
Step 4 is specially：As shown in figure 4, the amount of evening up of first tr1 and third road tr3 be can by first of tr1 with
What the amount of evening up and second tr2 of second tr2 were summed with the amount of evening up of third road tr3, this recurrence relation makes me
The amount of evening up of neighboring track in trace gather can be utilized to calculate the amounts of evening up of arbitrary twice in trace gather.Two benefits of this way
It is：1. solving the problems, such as lineups Severe distortion；2. solving two class AVO polarity Inversion Problems, (waveform between neighboring track is gradual change
).
The moving integration lineups of prestack trace gather using the present invention even up method, and it is same can to solve conventional Similar operator
The lineups distortion changing of the relative positions caused by phase axis is evened up and the two class AVO for solving the problems, such as lineups reversed polarity.
It these are only the preferred embodiment of the present invention, be not intended to restrict the invention, for those skilled in the art
For member, the invention may be variously modified and varied.Any modification made by all within the spirits and principles of the present invention,
Equivalent replacement, improvement etc., should all be included in the protection scope of the present invention.
Claims (2)
1. a kind of moving integration lineups of prestack trace gather even up method, which is characterized in that including：
Step 1, the amount of the evening up DT1 for calculating neighboring track in original trace gather；
The amount of evening up DT1 is applied in original trace gather by step 2, is generated a part and is evened up superposition trace gather；
The step 2 includes：
Step 21 laterally chooses sliding window in original trace gather；
Step 22, using the interim orbit of window as library track, trace gather in window is evened up and is superimposed；
Step 23, stack result are placed on the interim orbit of window, are obtained a part and are evened up superposition trace gather；
Step 3 evens up the amount of the evening up DT2 that superposition trace gather calculates neighboring track again using part；
Step 4, selection shortcut " superposition road ", implement moving integration to the amount of evening up DT2 and even up, even up entire trace gather.
2. the moving integration lineups of prestack trace gather as described in claim 1 even up method, which is characterized in that in step 1
In calculating process, counted come the length of window when controlling using the zero passage of waveform.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611238674.0A CN106501859B (en) | 2016-12-28 | 2016-12-28 | A kind of moving integration lineups of prestack trace gather even up method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611238674.0A CN106501859B (en) | 2016-12-28 | 2016-12-28 | A kind of moving integration lineups of prestack trace gather even up method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106501859A CN106501859A (en) | 2017-03-15 |
| CN106501859B true CN106501859B (en) | 2018-11-02 |
Family
ID=58334494
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611238674.0A Active CN106501859B (en) | 2016-12-28 | 2016-12-28 | A kind of moving integration lineups of prestack trace gather even up method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106501859B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107144875B (en) * | 2017-04-20 | 2019-03-15 | 中国石油天然气集团公司 | Seismic data corrections method and apparatus |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102565857A (en) * | 2011-12-16 | 2012-07-11 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Automatic remaining dynamic correction method |
| CN103744116A (en) * | 2014-01-23 | 2014-04-23 | 成都晶石石油科技有限公司 | Prestack gather all-time-domain phase-equalization correction method |
| CN105182420A (en) * | 2015-10-13 | 2015-12-23 | 中国石油天然气集团公司 | Dynamic-matching dynamic-correction method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9784866B2 (en) * | 2013-07-28 | 2017-10-10 | Geokinetics Usa, Inc. | Method and apparatus for enhanced monitoring of induced seismicity and vibration using linear low frequency and rotational sensors |
-
2016
- 2016-12-28 CN CN201611238674.0A patent/CN106501859B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102565857A (en) * | 2011-12-16 | 2012-07-11 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Automatic remaining dynamic correction method |
| CN103744116A (en) * | 2014-01-23 | 2014-04-23 | 成都晶石石油科技有限公司 | Prestack gather all-time-domain phase-equalization correction method |
| CN105182420A (en) * | 2015-10-13 | 2015-12-23 | 中国石油天然气集团公司 | Dynamic-matching dynamic-correction method |
Non-Patent Citations (2)
| Title |
|---|
| "一种动态的道集拉平方法研究及应用";周鹏 等;《断块油气田》;20150131;第22卷(第1期);第58-61、66页 * |
| "基于岩性预测的CRP道集优化处理";刘力辉 等;《石油物探》;20130930;第52卷(第5期);第482-488页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106501859A (en) | 2017-03-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| AU2005290207B2 (en) | Method for depth migrating seismic data using pre-stack time migration, demigration, and post-stack depth migration | |
| US9470811B2 (en) | Creating a high resolution velocity model using seismic tomography and impedance inversion | |
| CN105182420B (en) | A kind of Dynamic Matching moves bearing calibration | |
| CN108535775B (en) | Non-stationary seismic data sound impedance inversion method and device | |
| CN102890290B (en) | Prestack depth migration method under a kind of relief surface condition | |
| CN103460074A (en) | Method of wavelet estimation and multiple prediction in full wavefield inversion | |
| CN106226818A (en) | Seismic data processing technique and device | |
| CN110471113A (en) | Bearing calibration, device and storage medium are moved in inverting based on unstable state seismic data | |
| CN103616723A (en) | Amplitude versus offset (AVO) characteristic-based common reflection point (CRP) gather true amplitude recovery method | |
| CN104820242B (en) | A kind of road collection amplitude towards prestack inversion divides compensation method | |
| CN103499836A (en) | High-precision residual static correction method with combination between space variation and a plurality of time windows | |
| CN108957553B (en) | Method and device for dynamic correction without stretching distortion by recursion correction of dynamic correction value | |
| Chen et al. | Ray-parameter based stacking and enhanced pre-conditioning for stable inversion of receiver function data | |
| CN106501859B (en) | A kind of moving integration lineups of prestack trace gather even up method | |
| CN103913767A (en) | Method for solving energy nonuniform problem during multiple three-dimensional seismic data fusion processing before stacking | |
| CN106249292A (en) | A kind of CRP gather optimized treatment method and system | |
| Guo et al. | Becoming effective velocity-model builders and depth imagers, Part 2—The basics of velocity-model building, examples and discussions | |
| CN106324677B (en) | A kind of displacement field remaining dynamic correction method and device | |
| Garabito et al. | Global optimization of the common-offset CRS-attributes: synthetic and field data application | |
| White et al. | Optimised spectral merge of the background model in seismic inversion | |
| CN104635270B (en) | Converted wave staticses method based on composite traces constraint | |
| Filpo et al. | Image-ray concept as the key to 20 years of success of time-to-depth conversion in Petrobras | |
| CN113960660B (en) | Forward simulation based dynamic correction distortion area automatic identification and removal method | |
| CN103852786B (en) | A kind of method and system of the reverse-time migration imaging being applied to land seismic data | |
| Fliedner et al. | Velocity model building by wavefield-continuation imaging in the deepwater Gulf of Mexico |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |