CN102998700A - Surface consistency compensation method utilizing vertical seismic section data - Google Patents
Surface consistency compensation method utilizing vertical seismic section data Download PDFInfo
- Publication number
- CN102998700A CN102998700A CN2011102651037A CN201110265103A CN102998700A CN 102998700 A CN102998700 A CN 102998700A CN 2011102651037 A CN2011102651037 A CN 2011102651037A CN 201110265103 A CN201110265103 A CN 201110265103A CN 102998700 A CN102998700 A CN 102998700A
- Authority
- CN
- China
- Prior art keywords
- data
- energy
- seismic
- component
- arrival
- 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
Abstract
The invention relates to a surface consistency compensation method utilizing vertical seismic section data in seismic exploration data processing. Noise suppression processing and spherical diffusion compensation are performed to downhole three-component vertical seismic section data, neutron wave starting potion in the seismic data is acquired as first break, a time window is reserved in a vertical component backwards along the first break, quadratic sum of each three-component sample point data in the time window is calculated, the quadratic sum is extracted to serve as first break energy of a downstream wave of a signal, each compensated sample point data is divided by the first break energy of the downstream wave, and the process is repeated until surface consistency compensation of all seismic channels is finished. The first break energy of the downstream wave in an original wave field is mainly used for performing excitation energy and transmission energy consistency compensation, so that high fidelity is achieved. In addition, the surface consistency compensation method is simple to use and high in calculating efficiency.
Description
Technical field
The present invention relates to the Seismic Exploration Data Processing technology, is a kind of method of utilizing the earth's surface-consistent energy compensating of vertical seismic profile data.
Background technology
Wave detector can near the stratum, directly observe descending ripple and upward traveling wave from the stratum during collection of vertical seismic profile (VSP) (VSP) seismic data.Because the impact of other factors such as equipment, finish the VSP data collection of a bite well, need to repeatedly excite, excite the measurement that can only finish a degree of depth section at every turn, when excitation energy not simultaneously, the energy of the wave field integral body between big gun and big gun is not identical yet.For the difference of excitation energy between big gun, mainly utilize the energy of excitation wavelet to carry out excitation energy normalization at present, when not recording excitation wavelet, just there is not the method for good elimination excitation energy.
In addition, the descending ripple that receives for the wave detector of the different degree of depth has different transmission potential losses, this also causes same descending ripple also to exist capacity volume variance between the reflection upward traveling wave at same interface, cause the reflection coefficient that the actual upstream wave field that receives can not the complete reaction stratum, affect the result that seismic prospecting is processed and explained.In processing, seismic data yet there are no at present the method to the transmission potential compensating for loss and damage.
Summary of the invention
The object of the invention provides a kind of first arrival energy that utilizes descending ripple, and is easy to use succinct, the earth's surface-consistent energy compensation process of the vertical seismic profile data that operation efficiency is high.
Concrete steps of the present invention comprise:
1) surface source excites, and adopts the down-hole three-component seismometer to receive and obtains the vertical seismic profiling (VSP) data;
2) to step 1) in the geological data that obtains carry out various noise compression process;
3) to step 2) in the three-component seismic data of each depth track all carry out spherical diffusion compensation;
4) on through the vertical component after the spherical diffusion compensation, pick up the ski-jump of neutron deficiency of geological data as first arrival;
5) in vertical component, open for the moment window along first arrival backward, during to this in window three-component each sampling point data ask quadratic sum, then quadratic sum is extracted square root, as the first arrival energy of this road signal downlink ripple;
Window should just comprise a complete positive secondary lobe and a complete negative secondary lobe in the seismic signal when step 5) described.
6) to step 3) in the data of each sampling point of this road signal after spherical diffusion compensation divided by step 4) in the first arrival energy of descending ripple;
7) repeating step 2)-6), until all seismic traces are all divided by the descending ripple first arrival energy of this trace record, the earth's surface-consistent energy compensating is complete.
The present invention mainly utilizes the first arrival energy of descending ripple in the original wave field, carries out excitation energy and transmission potential uniformity compensation, has very high fidelity, and simultaneously easy to use succinct, operation efficiency is high.
Description of drawings
Fig. 1 actual zero is the VSP record diagram partially;
Fig. 2 excitation wavelet record diagram;
Window is selected synoptic diagram during Fig. 3;
Fig. 4 energy comparison diagram;
After Fig. 5 the present invention compensation zero be the VSP record diagram partially.
Embodiment
Describe the present invention in detail below in conjunction with accompanying drawing.
The present invention is in vertical seismic profile (VSP) (VSP) data, utilizes the first arrival energy of descending ripple, carries out the earth's surface-consistent energy compensation process.
Describe the present invention in detail below in conjunction with accompanying drawing.
Concrete steps of the present invention comprise:
1) adopt surface source to excite, the down-hole three-component seismometer receives and obtains the vertical seismic profiling (VSP) data; And near focus, utilize the three-component seismometer reception to obtain the excitation wavelet record.Fig. 1 is certain actual zero inclined to one side VSP record diagram; Fig. 2 is this zero inclined to one side excitation wavelet record diagram corresponding to each degree of depth of VSP data.
2) to step 1) in the geological data that obtains carry out the noise compression process;
3) to step 2) in the three-component seismic data of each depth track all carry out spherical diffusion compensation;
4) on through the vertical component after the spherical diffusion compensation, pick up the ski-jump of neutron deficiency of geological data as first arrival;
5) in vertical component, open for the moment window along first arrival backward, during to this in window three-component each sampling point data ask quadratic sum, then quadratic sum is extracted square root, as the first arrival energy of this road signal downlink ripple;
Window should just comprise a complete positive secondary lobe and a complete negative secondary lobe in the seismic signal when step 5) described.Window is selected shown in the synoptic diagram during such as Fig. 3, and included sampling point length is window length between two straight lines.
6) to step 3) in the data of each sampling point of this road signal after spherical diffusion compensation divided by step 4) in the first arrival energy of descending ripple; Fig. 4 energy comparison diagram can obviously be seen the excitation wavelet capacity volume variance.
Step 7) repeating step 2)-6), until all seismic traces are all divided by the descending ripple first arrival energy of this trace record, the earth's surface-consistent energy compensating is complete.
Fig. 5 is through zero inclined to one side VSP record diagram behind the earth's surface-consistent energy compensating of the present invention, eliminated the loss of excitation energy difference and transmission potential, the first arrival energy of the descending ripple of each depth track is identical (the descending ripple first arrival energy such as each depth track among Fig. 4 all is 1) all.
Claims (2)
1. the earth's surface-consistent energy compensation process of a vertical seismic profile data, characteristics are that concrete steps are:
1) surface source excites, and adopts the down-hole three-component seismometer to receive and obtains the vertical seismic profiling (VSP) data;
2) to step 1) in the geological data that obtains carry out the noise compression process;
3) three-component seismic data of each depth track is all carried out spherical diffusion compensation;
4) on through the vertical component after the spherical diffusion compensation, pick up the ski-jump of neutron deficiency of geological data as first arrival;
5) in vertical component, open for the moment window along first arrival backward, during to this in window three-component each sampling point data ask quadratic sum, then quadratic sum is extracted square root, as the first arrival energy of this road signal downlink ripple;
6) to step 3) in the data of each sampling point of this road signal after spherical diffusion compensation divided by step 4) in the first arrival energy of descending ripple;
7) repeating step 2)-6), until all seismic traces are all finished the earth's surface-consistent energy compensating divided by the descending ripple first arrival energy of this trace record.
2. method according to claim 1, characteristics are steps 5) window should comprise a complete positive secondary lobe and a complete negative secondary lobe in the seismic signal when described.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110265103.7A CN102998700B (en) | 2011-09-08 | 2011-09-08 | A kind of earth's surface-consistent energy compensation process of vertical seismic profile data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110265103.7A CN102998700B (en) | 2011-09-08 | 2011-09-08 | A kind of earth's surface-consistent energy compensation process of vertical seismic profile data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102998700A true CN102998700A (en) | 2013-03-27 |
| CN102998700B CN102998700B (en) | 2015-09-23 |
Family
ID=47927506
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201110265103.7A Active CN102998700B (en) | 2011-09-08 | 2011-09-08 | A kind of earth's surface-consistent energy compensation process of vertical seismic profile data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102998700B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104375188A (en) * | 2014-11-19 | 2015-02-25 | 中国石油天然气集团公司 | Seismic wave transmission attenuation compensation method and device |
| CN106483563A (en) * | 2015-08-25 | 2017-03-08 | 中国石油天然气股份有限公司 | seismic energy compensation method based on complementary set empirical mode decomposition |
| CN110954940A (en) * | 2018-09-26 | 2020-04-03 | 中国石油化工股份有限公司 | First arrival quality control method based on earth surface consistency model estimation |
| CN111060968A (en) * | 2019-12-24 | 2020-04-24 | 恒泰艾普（北京）能源科技研究院有限公司 | Method for correcting accurate well control speed through spherical diffusion compensation |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4464737A (en) * | 1981-09-28 | 1984-08-07 | Mobil Oil Corporation | Method for migration of seismic reflection waves |
| US20090097356A1 (en) * | 2003-11-14 | 2009-04-16 | Schlumberger Technology Corporation | Processing of combined surface and borehole seismic data |
| CN101630015A (en) * | 2008-07-16 | 2010-01-20 | 中国石油天然气集团公司 | Method for improving precision and efficiency of first-arrival wave pickup |
| CN102176056A (en) * | 2011-02-18 | 2011-09-07 | 中国石油化工股份有限公司 | Seismic effective wave multi-domain energy compensating method based on time-frequency analysis |
-
2011
- 2011-09-08 CN CN201110265103.7A patent/CN102998700B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4464737A (en) * | 1981-09-28 | 1984-08-07 | Mobil Oil Corporation | Method for migration of seismic reflection waves |
| US20090097356A1 (en) * | 2003-11-14 | 2009-04-16 | Schlumberger Technology Corporation | Processing of combined surface and borehole seismic data |
| CN101630015A (en) * | 2008-07-16 | 2010-01-20 | 中国石油天然气集团公司 | Method for improving precision and efficiency of first-arrival wave pickup |
| CN102176056A (en) * | 2011-02-18 | 2011-09-07 | 中国石油化工股份有限公司 | Seismic effective wave multi-domain energy compensating method based on time-frequency analysis |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104375188A (en) * | 2014-11-19 | 2015-02-25 | 中国石油天然气集团公司 | Seismic wave transmission attenuation compensation method and device |
| CN106483563A (en) * | 2015-08-25 | 2017-03-08 | 中国石油天然气股份有限公司 | seismic energy compensation method based on complementary set empirical mode decomposition |
| CN110954940A (en) * | 2018-09-26 | 2020-04-03 | 中国石油化工股份有限公司 | First arrival quality control method based on earth surface consistency model estimation |
| CN111060968A (en) * | 2019-12-24 | 2020-04-24 | 恒泰艾普（北京）能源科技研究院有限公司 | Method for correcting accurate well control speed through spherical diffusion compensation |
| CN111060968B (en) * | 2019-12-24 | 2022-10-21 | 恒泰艾普（北京）能源科技研究院有限公司 | Method for correcting accurate well control speed through spherical diffusion compensation |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102998700B (en) | 2015-09-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102841373B (en) | Microseism positioning method based on azimuth angle constraint | |
| WO2017024523A1 (en) | Inversion method for ray elastic parameter | |
| CN102176054B (en) | Near-surface comprehensive information processing explanation method | |
| CN101419292B (en) | Method for generating transverse wave seismic section by longitudinal wave source multi-component earthquake data | |
| CN103399300B (en) | Wave packet superposition microseism ground location method | |
| CN1873441A (en) | Method for suppressing background noise not related to excitation source in seismic data enhancement | |
| CN103399346B (en) | A kind of well shake associating impedance initial value modeling method | |
| CN103645497B (en) | Emergence angle based multi-component wave field separation method | |
| CN103645503A (en) | Three-dimensional time domain illumination analysis and amplitude compensation method | |
| CN102721977A (en) | High-precision converted wave anisotropy stacking velocity analysis method | |
| CN102073067A (en) | Method for increasing automatic pickup efficiency of geological data first-motion waves | |
| CN102998700B (en) | A kind of earth's surface-consistent energy compensation process of vertical seismic profile data | |
| CN110824563A (en) | Reservoir lithology prediction method based on Xgboost algorithm | |
| CN104090299B (en) | Surface seismic data amplitude compensation method based on VSP primary waves | |
| CN105093319A (en) | Ground micro-seismic static correction method based on three-dimensional seismic data | |
| CN105607119B (en) | Near-surface model construction method and static correction value acquiring method | |
| CN102830170B (en) | Control method and control device for obtaining coal sample transverse wave signal based on ultrasonic test | |
| CN103149588A (en) | Method and system for calculating VTI (Velocity Time Integral) anisotropic parameter by utilizing well seismic calibration | |
| CN105093296A (en) | Method and device for optimizing observation system | |
| CN105301647A (en) | Method for distinguishing calcareous mudstone from sandstone | |
| CN104133244B (en) | A kind of near surface static corrections processing method of focus shear wave exploration | |
| CN105301658A (en) | Method for extracting near-surface deconvolution operator by utilizing record of microlog | |
| CN108594301B (en) | A kind of method and processing terminal of the seismic data fusion with difference characteristic | |
| CN101950032A (en) | Multi-wave exciting method for near surface investigation | |
| CN102914797A (en) | Method and device for acquiring anisotropy coefficient of stratum |
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 |