CN103576187A - Method determining multi-wave static correction value employing reflection imaging time difference - Google Patents
Method determining multi-wave static correction value employing reflection imaging time difference Download PDFInfo
- Publication number
- CN103576187A CN103576187A CN201210250364.6A CN201210250364A CN103576187A CN 103576187 A CN103576187 A CN 103576187A CN 201210250364 A CN201210250364 A CN 201210250364A CN 103576187 A CN103576187 A CN 103576187A
- Authority
- CN
- China
- Prior art keywords
- imaging time
- shot
- static correction
- actual reflection
- reflection interval
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention relates to a method determining a multi-wave static correction value employing reflection imaging time difference. Processing on all shot points and detection points on a multi-wave detection line is carried out, common shot point trace superimposition and common detection point trace superimposition are firstly carried out to respectively generate two superimposition trace composite sections, and real time of shallow reflection co-phase shafts of the acquired sections is respectively picked; based on an earth surface consistency static correction principle, static correction mainly influences shallow imaging, so real time can be relatively continuous and smooth after eliminating influence of the static correction, smoothing on the real time is properly carried out to acquire imaging time after smoothing, and the static correction value is acquired employing imaging time difference of the imaging time and the real time. The method avoids effective identification difficulty and picking difficulty during conversion wave first arrival, solves a problem of relatively great residual amount on the conversion wave section caused by severe co-phase shaft jittering.
Description
Technical field
The present invention relates to multi-wave seismic exploration engineering, a kind ofly utilize the method that the catoptric imaging time difference calculates many ripples static correction value.
Background technology
Multi-wave seismic exploration is the effective ways that the much informations such as a comprehensive utilization compressional wave, shear wave, transformed wave carry out meticulous exploration, Direct Prediction of Oil.Although the applicating history of this technology is not long, many challenging achievements in application, have been obtained.Some old oil gas fields enter development late stage mostly, and remaining oil distribution high degree of dispersion is effectively identified thin interbed, minor fault, low amplitude structure, and Accurate Prediction reservoir border and reservoir properties, monitoring fluid all become the key factor of old oil gas field stable yields.Utilize multi-wave seismic technology to overcome the above problems and there is vital role, along with the development of many ripples exploration engineering and raising, the problem that many ripple explorations exist also highlights gradually, particularly shear wave static correction problem does not obtain the breakthrough of essence, and its reason is mainly in many complex areas, and compressional wave static correction problem does not also solve at present completely, shear wave is because corresponding weathering zone is more complicated, data signal to noise ratio (S/N ratio) is lower, adds the impact of the azimuthal anisotropy that may exist, and difficulty is more much bigger than compressional wave static correction.
The horizontal wave propagation velocity of shallow-layer is relatively low, after proofreading and correct, is conducive to the contrast of ripple in length and breadth, is conducive to obtain the dynamic characteristic of more real ripple.But conventional many ripples static correction poor effect based on many ripples near-surface investigation and transformed wave first arrival, on most of many wave profiles seriously, there is larger surplus (being greater than 1/2 cycle of reflection wave) in lineups shake, and many wave tables layer investigation precision is lower; Transformed wave first arrival is difficult to effectively identify and pick up, and the static correction value calculating exists larger error, cannot meet multiwave data processing requirement.
Summary of the invention
The object of the invention is to provide a kind of investigation precision high, and the catoptric imaging time difference of utilizing of sufficient data processing requirements is determined the method for many ripples static correction value.
The present invention realizes by following steps:
1) gather multi-wave seismic data;
2) all big guns on many ripples survey line are carried out to shot record migration trace-stacking, generate shot record migration trace-stacking section;
3) all geophone stations on many ripples survey line are carried out to common receiver trace-stacking, generate common receiver trace-stacking section;
4) on shot record migration trace-stacking section, pick up the actual reflection interval of the reflection horizon lineups that shallow-layer signal to noise ratio (S/N ratio) is high, form shot point actual reflection interval of curve;
Be to have continuity on shot record migration trace-stacking section the actual reflection interval of the reflection horizon lineups that described signal to noise ratio (S/N ratio) is high, the actual reflection interval of traceable reflection horizon lineups.
5) on common receiver trace-stacking section, pick up the reflection horizon lineups actual reflection interval corresponding with shot record migration trace-stacking section, form geophone station actual reflection interval of curve;
6) respectively shot point actual reflection interval of curve and geophone station actual reflection interval of curve are carried out smoothly, make actual reflection interval of continuous imaging, obtain respectively shot point imaging time, geophone station imaging time;
7) utilize shot point imaging time and shot point actual reflection interval of difference, calculate shot statics, computing formula is: shot statics=(actual reflection interval of shot point imaging time-shot point)/2;
8) utilize geophone station imaging time and geophone station actual reflection interval of difference, calculate receiver static correction amount, computing formula is: receiver static correction amount=(actual reflection interval of geophone station imaging time-geophone station)/2.
The present invention utilizes reflected wave information, the continuity of picking up respectively shallow-layer on the folded section of shot record migration trace-stacking and common receiver road is good, the hourage of traceable high s/n ratio lineups, to carrying out hourage smoothly, obtain the correction time of picked up lineups, after correction, these lineups more continuously and convergence.Based on surface-consistent statics principle, now shallow, in, the static correction problem of deep layer all can obtain corresponding solution.This method is avoided transformed wave first arrival and is difficult to the difficult problem of effectively identifying and picking up, and can solve lineups shake on converted wave sectional plane and seriously, have the problem of larger surplus (being greater than 1/2 cycle of reflection wave).
Accompanying drawing explanation
Fig. 1 imaging time difference of the present invention is asked for schematic diagram:
Fig. 2 the present invention utilizes the catoptric imaging time difference to calculate many ripples static correction value process flow diagram.
Embodiment
The present invention is by all shot points and geophone station on many ripples survey line are processed, and first carries out shot record migration trace-stacking and common receiver trace-stacking and generates respectively two stack road composite sections, the real time of picking up respectively gained section superficial reflex lineups; Then based on surface-consistent statics principle, the imaging of static correction major effect shallow-layer, therefore it should be smoother eliminating the rear real time of static correction impact, real time is done and suitably smoothly obtains imaging time (time after level and smooth), then utilize the imaging time difference of imaging time and real time to ask for static correction value, see Fig. 1.
The concrete implementation step of the present invention is as follows:
1) gather multi-wave seismic data;
2) all big guns on many ripples survey line are carried out to shot record migration trace-stacking, generate shot record migration trace-stacking section;
3) all geophone stations on many ripples survey line are carried out to common receiver trace-stacking, generate common receiver trace-stacking section;
4) on shot record migration trace-stacking section, pick up the actual reflection interval of the reflection horizon lineups that shallow-layer signal to noise ratio (S/N ratio) is high, form shot point actual reflection interval of curve;
Be to have continuity on shot record migration trace-stacking section the actual reflection interval of the reflection horizon lineups that described signal to noise ratio (S/N ratio) is high, the actual reflection interval of traceable reflection horizon lineups.
5) on common receiver trace-stacking section, pick up the reflection horizon lineups actual reflection interval corresponding with shot record migration trace-stacking section, form geophone station actual reflection interval of curve;
6) respectively shot point actual reflection interval of curve and geophone station actual reflection interval of curve are carried out smoothly, making originally to have fluctuating in Excel software, actual reflection interval of the curve of the changing of the relative positions becomes and rises and falls more slowly, there is no the changing of the relative positions.Time curve after level and smooth is respectively shot point imaging time, geophone station imaging time;
7) utilize shot point imaging time and shot point actual reflection interval of difference, calculate shot statics, computing formula is: shot statics=(actual reflection interval of shot point imaging time-shot point)/2;
8) utilize geophone station imaging time and geophone station actual reflection interval of difference, calculate receiver static correction amount, computing formula is: receiver static correction amount=(actual reflection interval of geophone station imaging time-geophone station)/2.
Other implementation steps of the present invention are as follows:
Step 1: when carrying out picking up reflection interval, can carry out segmentation to the different reflection horizon of shallow-layer and pick up, not need to follow the trail of same reflection horizon, the lineups that the signal to noise ratio (S/N ratio) of take is relatively high are first-selection.
Step 2: in the time of between reflection time picking, pick up reflection at peak as far as possible, after correction, can guarantee that reflection wave crest is consistent, lineups have a good continuity.
Claims (2)
1. utilize the catoptric imaging time difference to determine a method for many ripples static correction value, feature is to realize by following steps:
1) gather multi-wave seismic data;
2) all big guns on many ripples survey line are carried out to shot record migration trace-stacking, generate shot record migration trace-stacking section;
3) all geophone stations on many ripples survey line are carried out to common receiver trace-stacking, generate common receiver trace-stacking section;
4) on shot record migration trace-stacking section, pick up the actual reflection interval of the reflection horizon lineups that shallow-layer signal to noise ratio (S/N ratio) is high, form shot point actual reflection interval of curve;
5) on common receiver trace-stacking section, pick up the reflection horizon lineups actual reflection interval corresponding with shot record migration trace-stacking section, form geophone station actual reflection interval of curve;
6) respectively shot point actual reflection interval of curve and geophone station actual reflection interval of curve are carried out smoothly, make actual reflection interval of continuous imaging, obtain respectively shot point imaging time, geophone station imaging time;
7) utilize shot point imaging time and shot point actual reflection interval of difference, calculate shot statics, computing formula is: shot statics=(actual reflection interval of shot point imaging time-shot point)/2;
8) utilize geophone station imaging time and geophone station actual reflection interval of difference, calculate receiver static correction amount, computing formula is: receiver static correction amount=(actual reflection interval of geophone station imaging time-geophone station)/2.
2. according to the method for claim 1, feature is to be the actual reflection interval of the reflection horizon lineups that the signal to noise ratio (S/N ratio) described in step 4) is high to have continuity on shot record migration trace-stacking section, the actual reflection interval of traceable reflection horizon lineups.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210250364.6A CN103576187B (en) | 2012-07-19 | 2012-07-19 | A kind of method utilizing the catoptric imaging time difference to determine many ripples static correction value |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210250364.6A CN103576187B (en) | 2012-07-19 | 2012-07-19 | A kind of method utilizing the catoptric imaging time difference to determine many ripples static correction value |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103576187A true CN103576187A (en) | 2014-02-12 |
| CN103576187B CN103576187B (en) | 2016-08-10 |
Family
ID=50048353
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210250364.6A Active CN103576187B (en) | 2012-07-19 | 2012-07-19 | A kind of method utilizing the catoptric imaging time difference to determine many ripples static correction value |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103576187B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111596359A (en) * | 2019-02-21 | 2020-08-28 | 中国石油天然气集团有限公司 | Inclined stratum wide line data correction imaging method and device |
| CN112394412A (en) * | 2020-10-30 | 2021-02-23 | 中国石油天然气集团有限公司 | Shear wave static correction method and device based on shallow layer superposition modeling |
| CN112462426A (en) * | 2020-11-02 | 2021-03-09 | 中国石油天然气集团有限公司 | Shear wave vector static correction method and device |
| CN113189646A (en) * | 2021-04-19 | 2021-07-30 | 滨州学院 | Method for removing dragging type shallow-section stratum abnormal fluctuation |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1365008A (en) * | 2001-01-19 | 2002-08-21 | 中国石油天然气股份有限公司 | Earthquack polyregion interative static correction method |
-
2012
- 2012-07-19 CN CN201210250364.6A patent/CN103576187B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1365008A (en) * | 2001-01-19 | 2002-08-21 | 中国石油天然气股份有限公司 | Earthquack polyregion interative static correction method |
Non-Patent Citations (6)
| Title |
|---|
| FENG FAQUANG ET AL.: "《The multivariate statistic residual static correction method and examples of its application》", 《SEG SAN ANTONIO 2011 ANNUAL MEETING》 * |
| S.CUI ET AL.: "《On-land Converted Wave Static Correction》", 《71ST EAGE CONFERENCE & EXHIBITION》 * |
| 童思友等: "《三维折射波静校正方法及应用》", 《石油地球物理勘探》 * |
| 赵峰等: "《叠加法静校正》", 《石油地球物理勘探》 * |
| 黄明忠等: "《三维炮检域迭代直接静校正技术及应用》", 《勘探地球物理进展》 * |
| 黄明忠等: "《炮检域迭代直接静校正方法及应用效果》", 《勘探地球物理进展》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111596359A (en) * | 2019-02-21 | 2020-08-28 | 中国石油天然气集团有限公司 | Inclined stratum wide line data correction imaging method and device |
| CN112394412A (en) * | 2020-10-30 | 2021-02-23 | 中国石油天然气集团有限公司 | Shear wave static correction method and device based on shallow layer superposition modeling |
| CN112462426A (en) * | 2020-11-02 | 2021-03-09 | 中国石油天然气集团有限公司 | Shear wave vector static correction method and device |
| CN112462426B (en) * | 2020-11-02 | 2024-05-28 | 中国石油天然气集团有限公司 | Transverse wave vector static correction method and device |
| CN113189646A (en) * | 2021-04-19 | 2021-07-30 | 滨州学院 | Method for removing dragging type shallow-section stratum abnormal fluctuation |
| CN113189646B (en) * | 2021-04-19 | 2022-08-30 | 滨州学院 | Method for removing dragging type shallow-section stratum abnormal fluctuation |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103576187B (en) | 2016-08-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102193109B (en) | Direct prestack time migration method for three-dimensional seismic data acquired from irregular surfaces | |
| CN103091710B (en) | A kind of reverse-time migration formation method and device | |
| CN105549087B (en) | A kind of underground coal mine seam seismic exploration when walking and amplitude joint inversion method | |
| CN104570096B (en) | A kind of granite buried hill top surface fault interpretation method | |
| CN104536043A (en) | Depth domain overall velocity model combination method and device | |
| CN102879819B (en) | Seismic data processing method and device for keeping kinematics characteristics of seismic wave field | |
| CN104570079A (en) | Time matching method of longitudinal wave and converted shear wave seismic data | |
| CN104656139A (en) | Static correction method and equipment for converted wave | |
| CN103576187A (en) | Method determining multi-wave static correction value employing reflection imaging time difference | |
| CN102692648A (en) | Refraction wave residual static correction method based on cannon first arrival | |
| CN104090301B (en) | A kind of method for asking for three-D high frequency static correction value | |
| CN102073067A (en) | Method for increasing automatic pickup efficiency of geological data first-motion waves | |
| CN104570073B (en) | A kind of bireflectance seismic imaging method suitable for complicated high-dip structure | |
| CN103630934A (en) | Method for determining large transverse wave static correction value of converted wave demodulator probe | |
| CN107219554A (en) | The automatic obtaining method of the Value of residual static correction of land seismic data | |
| CN104199103A (en) | Static correction method and static correction device for converted wave | |
| CN104977615B (en) | A kind of multiple ripple drawing method of deep water OBC data based on modeling statistics pickup | |
| CN104155694A (en) | Residual static correction method of reflection converted shear wave common geophone stack section | |
| CN105607119A (en) | Near-surface model construction method and static correction value calculation method | |
| CN105116447A (en) | A geological river direction discrimination method based on curvature-abnormal stripes | |
| CN103487834A (en) | Converted wave total inspection wave point overlapping static correction method | |
| CN109100798B (en) | Realize the method, apparatus and processing terminal of refraction multiple wave tomographic inversion | |
| CN104216018A (en) | Surface structure survey method based on small track pitch tomography inversion | |
| CN102798893B (en) | Three-dimensional superimposed tomography inversion static correction method | |
| CN109725354B (en) | Anisotropic speed modeling method and system |
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 | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |