CN105445802A - High-precision time-depth correction method based on time-frequency spectrum constraints - Google Patents
High-precision time-depth correction method based on time-frequency spectrum constraints Download PDFInfo
- Publication number
- CN105445802A CN105445802A CN201410449078.1A CN201410449078A CN105445802A CN 105445802 A CN105445802 A CN 105445802A CN 201410449078 A CN201410449078 A CN 201410449078A CN 105445802 A CN105445802 A CN 105445802A
- Authority
- CN
- China
- Prior art keywords
- time
- frequency spectrum
- seismic
- depth
- data
- 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.)
- Pending
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a high-precision time-depth correction method based on time-frequency spectrum constraints. The high-precision time-depth correction method comprises the steps of: acquiring seismic data and acoustic logging data used for time-depth correction; acquiring a seismic synthetic recording trace through wavelet extraction and reflection coefficient calculation, and calibrating the seismic synthetic recording trace in a comprehensive calibration system; respectively calculating time-frequency spectra of the seismic and logging data, and displaying synthetic seismic record, cross-well seismic trace data as well as seismic and logging time-frequency spectra in the comprehensive calibration system; comparing a logging time-frequency spectrum interface with a seismic time-frequency spectrum interface and finely adjusting the interfaces under horizon constraint, and meeting the requirement of seismic synthetic record calibration and time-frequency spectrum interface matching; and finally, saving or outputting a time-depth relationship after completing the calibration. The high-precision time-depth correction method utilizes time interface features such as time-frequency spectra of logging and seismic data, organically combines the seismic synthetic record calibration method and the time-frequency spectrum calibration method, and effectively increases logging and seismic time-depth calibration precision.
Description
Technical field
The present invention relates to geophysical survey seismic interpretation and synthetic study field, particularly relate to a kind of high precision time-depth rectify method
.
Background technology
In geophysical survey seismic interpretation and synthetic study, earthquake and well-log information are two kinds of the most frequently used basic datas, and they belong to time domain and Depth Domain usually respectively.Well-log information can obtain comparatively meticulous and reliable geological information, and seismic data has the continuity of horizontal extrapolation, has more spatial feature.Both combination needs time and depth transfer or using in soil deeply realizes.
Time and depth transfer is necessary treatment step geological data being changed into Depth Domain signal from time-domain signal, is to utilize seismic data to carry out one of structure elucidation very the key link.At present conventional time and depth transfer method is dark relation when obtaining by the theogram horizon calibration of routine.Conventional theogram method is the one dimension convolution model under employing HORIZONTAL LAYERED MEDIUM WITH HIGH ACCURACY assumed condition, do not consider the impact of the factors such as transmission loss, multiple reflection, transformed wave, theogram and seismic trace near well is caused to there are differences on across comparison, and due to the impact that composite traces resolution is lower, the time dark relation precision obtained after causing demarcation is not high.Well shake time-frequency spectrum can obtain the isochronic plane of high-resolution, under the constraint of conventional calibration result, and the precision of deeply demarcating when can be improved by the contrast of these isochronic plane therein, thus obtain the higher time dark relation of a precision.
Summary of the invention
The object of the invention is for the deficiencies in the prior art, a kind of high precision time-depth rectify method based on time-frequency spectrum constraint is provided, solves the problem that time dark stated accuracy that is that cause is low because conventional composite traces resolution is lower.
The present invention adopts following technical scheme:
First obtain an initial time dark relation by the horizon calibration of conventional composite traces, axle and geological data can keep a rational relativity in the same way to make the key reflections of synthetic seismic record.Under the constraint of composite traces, adopt the high-precision calibrating method of well shake time-frequency spectrum isochronic plane obtain one more high-precision time dark relation, specifically adopt following steps:
1) obtain geological data for time-depth rectify and sound wave measuring well curve by data acquisition, determine the depth range of well data;
2) obtain synthetic seismic record road by wavelet extraction and reflection coefficient calculating, the composite traces obtained is carried out dark relation when demarcation obtains initial by composite calibration system;
3) time-frequency spectrum of earthquake and log data is calculated respectively, and regulate well logging time-frequency spectrum sampling rate, make earthquake consistent with well logging time-frequency spectrum resolution, synthetic seismic record, excessively well seismic channel data, earthquake and well logging time-frequency spectrum are shown in composite calibration system respectively;
4) under layer restrain, compared in well logging time-frequency spectrum interface and earthquake time-frequency spectrum interface, the corresponding time-frequency spectrum interface of fine setting well logging makes to be consistent with earthquake time-frequency spectrum interface, by stretching or both having been met dark relation when earthquake composite seismogram calibration also meets the high precision of time-frequency spectrum interphase match after compressing adjustment;
5) finally the time dark relation after having demarcated is saved in database or outputs in external file.
The time-frequency spectrum calculating earthquake and log data described in such scheme comprises time-frequency spectrum instrument and generates and time-frequency spectrum generation.Described time-frequency spectrum instrument generation method comprises triangle filtering method, small wave converting method and S-transformation method.Described composite calibration system comprises translation, Tension and Compression subsystem.
The present invention utilizes well to shake the time-frequency spectrum isochronic plane feature of data, synthetic seismic record scaling method and time-frequency spectrum scaling method is organically combined, and effectively improves the precision of well shake time-depth rectify.The time dark relation particularly obtained ensure that the consistance of well shake sedimentary cycle interface feature, contributes to the research work such as well shake isochronal formation cor-relation.
Accompanying drawing explanation
Fig. 1 a kind of implementing procedure figure of the present invention
Dark relation contrast time before and after the demarcation of Fig. 2 time-frequency spectrum
Fig. 3 figure composite traces initial alignment result
Fig. 4 a time-frequency spectrum demarcates the contrast of front well shake time-frequency spectrum
Fig. 4 b time-frequency spectrum demarcates the contrast of rear well shake time-frequency spectrum
Map Table Header information illustrates:
DVST---time dark relation curve road
Sonic---sound wave curve road
Scale---scale road
RC---reflection coefficient road
Wavelet---wavelet road
Synthetic_record---synthetic traces
Seismic---seismic-data traces
Spectrum---time-frequency spectrum road.
Embodiment
Below in conjunction with accompanying drawing, the present invention is described in further detail:
The invention provides a kind of high precision time-depth rectify method and apparatus based on time-frequency spectrum constraint, with reference to Fig. 1, its concrete implementation step is as follows:
(1) obtain geological data for time-depth rectify and sound wave measuring well curve by data acquisition, determine the depth range of well data.
(2) can obtain synthetic seismic record road by wavelet extraction and reflection coefficient calculating, the composite traces obtained is carried out dark relation when demarcation obtains initial by composite calibration system, and Fig. 3 is composite traces initial alignment result.
(3) calculated the time-frequency spectrum calculating earthquake and log data respectively by time-frequency spectrum, and regulate well logging time-frequency spectrum sampling rate, make earthquake consistent with well logging time-frequency spectrum resolution.Synthetic seismic record, excessively well seismic channel data, earthquake and well logging time-frequency spectrum are shown in composite calibration system respectively.Fig. 4 a is well shake time-frequency spectrum display situation before time-frequency is demarcated.
(4) under layer restrain, compared in well logging time-frequency spectrum interface and earthquake time-frequency spectrum interface, the corresponding time-frequency spectrum interface of fine setting well logging makes to be consistent with earthquake time-frequency spectrum interface, both can be met dark relation when earthquake composite seismogram calibration also meets the high precision of time-frequency spectrum interphase match after adjustment.Fig. 4 b is well shake time-frequency spectrum contrast after time-frequency is demarcated, and can find out when composite traces remains unchanged, and well shake time-frequency spectrum isochronic plane relation is more consistent.Fig. 2 show time-frequency spectrum demarcate before and after time dark relation contrast situation, can find out time-frequency spectrum demarcate before and after time dark relation basically identical, dark relation trend when maintaining reason.
(5) export the time dark relation after having demarcated can be saved in database or output in external file finally by data.
Claims (4)
1., based on a high precision time-depth rectify method for time-frequency spectrum constraint, it is characterized in that comprising:
1) obtain geological data for time-depth rectify and sound wave measuring well curve by data acquisition, determine the depth range of well data;
2) obtain synthetic seismic record road by wavelet extraction and reflection coefficient calculating, the composite traces obtained is carried out dark relation when demarcation obtains initial by composite calibration system;
3) time-frequency spectrum of earthquake and log data is calculated respectively, and regulate well logging time-frequency spectrum sampling rate, make earthquake consistent with well logging time-frequency spectrum resolution, synthetic seismic record, excessively well seismic channel data, earthquake and well logging time-frequency spectrum are shown in composite calibration system respectively;
4) under layer restrain, compared in well logging time-frequency spectrum interface and earthquake time-frequency spectrum interface, the corresponding time-frequency spectrum interface of fine setting well logging makes to be consistent with earthquake time-frequency spectrum interface, by stretching or both having been met dark relation when earthquake composite seismogram calibration also meets the high precision of time-frequency spectrum interphase match after compressing adjustment;
5) finally the time dark relation after having demarcated is saved in database or outputs in external file.
2., as claimed in claim 1 based on the high precision time-depth rectify method of time-frequency spectrum constraint, it is characterized in that the time-frequency spectrum of described calculating earthquake and log data comprises time-frequency spectrum instrument and generates and time-frequency spectrum generation.
3., as claimed in claim 2 based on the high precision time-depth rectify method of time-frequency spectrum constraint, it is characterized in that time-frequency spectrum instrument generation method comprises triangle filtering method, small wave converting method and S-transformation method.
4. the high precision time-depth rectify method based on time-frequency spectrum constraint as described in claim 1 or 2,3, is characterized in that described composite calibration system comprises translation, Tension and Compression subsystem.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410449078.1A CN105445802A (en) | 2014-09-04 | 2014-09-04 | High-precision time-depth correction method based on time-frequency spectrum constraints |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410449078.1A CN105445802A (en) | 2014-09-04 | 2014-09-04 | High-precision time-depth correction method based on time-frequency spectrum constraints |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN105445802A true CN105445802A (en) | 2016-03-30 |
Family
ID=55556214
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410449078.1A Pending CN105445802A (en) | 2014-09-04 | 2014-09-04 | High-precision time-depth correction method based on time-frequency spectrum constraints |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105445802A (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106569280A (en) * | 2016-10-28 | 2017-04-19 | 中国石油化工股份有限公司 | Calibration method for improving seismic interpretation precision of fault terrace structures |
| CN106610505A (en) * | 2016-12-29 | 2017-05-03 | 中国石油大学(华东) | Logging and seismic data matching method based on combination of DTW (Dynamic Time Warping) and ABA (Amplitude Balance Adjustment) |
| CN109782343A (en) * | 2018-12-13 | 2019-05-21 | 中国石油天然气集团有限公司 | Stratigraphic cycles analysis method and device |
| WO2020131424A1 (en) * | 2018-12-20 | 2020-06-25 | Chevron U.S.A. Inc. | Methods and systems for calibrating depth in a well to seismic data in a subsurface volume of interest |
| CN112147717A (en) * | 2020-09-27 | 2020-12-29 | 北京中恒利华石油技术研究所 | Pre-stack and post-stack well-seismic combination time depth calibration method |
| CN115096756A (en) * | 2022-07-27 | 2022-09-23 | 浙江双元科技股份有限公司 | Self-calibration surface density detector and calibration method thereof |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101201407A (en) * | 2006-12-12 | 2008-06-18 | 中国石油天然气集团公司 | Relative non-high-frequency leakage equivalent N-drop smooth spectrum analog deconvolution method |
| CN101907725A (en) * | 2010-06-12 | 2010-12-08 | 中国石油天然气股份有限公司 | Crack prediction method and device |
| US20110085415A1 (en) * | 2009-10-09 | 2011-04-14 | Morton Peter J | Method for synchronizing seismic data recorded by two or more separate recording systems |
| CN103558634A (en) * | 2013-10-29 | 2014-02-05 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method and device for obtaining frequency attenuation gradient of seismic data |
| CN103790569A (en) * | 2012-10-26 | 2014-05-14 | 中国石油化工股份有限公司 | Method for correcting acoustic logging information through VSP (vertical seismic profiling) time-depth relationship |
| CN103837893A (en) * | 2014-03-10 | 2014-06-04 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Acoustic logging curve correction method based on VSP data |
-
2014
- 2014-09-04 CN CN201410449078.1A patent/CN105445802A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101201407A (en) * | 2006-12-12 | 2008-06-18 | 中国石油天然气集团公司 | Relative non-high-frequency leakage equivalent N-drop smooth spectrum analog deconvolution method |
| US20110085415A1 (en) * | 2009-10-09 | 2011-04-14 | Morton Peter J | Method for synchronizing seismic data recorded by two or more separate recording systems |
| CN101907725A (en) * | 2010-06-12 | 2010-12-08 | 中国石油天然气股份有限公司 | Crack prediction method and device |
| CN103790569A (en) * | 2012-10-26 | 2014-05-14 | 中国石油化工股份有限公司 | Method for correcting acoustic logging information through VSP (vertical seismic profiling) time-depth relationship |
| CN103558634A (en) * | 2013-10-29 | 2014-02-05 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method and device for obtaining frequency attenuation gradient of seismic data |
| CN103837893A (en) * | 2014-03-10 | 2014-06-04 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Acoustic logging curve correction method based on VSP data |
Non-Patent Citations (2)
| Title |
|---|
| 刘百红等: "时频分析技术及其在储层预测中的应用", 《石油地球物理进展》 * |
| 杨帅等: "精细井震分频标定及其特殊检验方法研究", 《科学技术与工程》 * |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106569280A (en) * | 2016-10-28 | 2017-04-19 | 中国石油化工股份有限公司 | Calibration method for improving seismic interpretation precision of fault terrace structures |
| CN106610505A (en) * | 2016-12-29 | 2017-05-03 | 中国石油大学(华东) | Logging and seismic data matching method based on combination of DTW (Dynamic Time Warping) and ABA (Amplitude Balance Adjustment) |
| CN106610505B (en) * | 2016-12-29 | 2019-03-22 | 中国石油大学(华东) | One kind shaking material matches method based on the united well of DTW and ABA |
| CN109782343A (en) * | 2018-12-13 | 2019-05-21 | 中国石油天然气集团有限公司 | Stratigraphic cycles analysis method and device |
| CN109782343B (en) * | 2018-12-13 | 2020-07-10 | 中国石油天然气集团有限公司 | Stratum convolution analysis method and device |
| WO2020131424A1 (en) * | 2018-12-20 | 2020-06-25 | Chevron U.S.A. Inc. | Methods and systems for calibrating depth in a well to seismic data in a subsurface volume of interest |
| US11275193B2 (en) | 2018-12-20 | 2022-03-15 | Chevron U.S.A. Inc. | Methods and systems for calibrating depth in a well to seismic data in a subsurface volume of interest |
| CN112147717A (en) * | 2020-09-27 | 2020-12-29 | 北京中恒利华石油技术研究所 | Pre-stack and post-stack well-seismic combination time depth calibration method |
| CN112147717B (en) * | 2020-09-27 | 2021-11-16 | 北京中恒利华石油技术研究所 | Pre-stack and post-stack well-seismic combination time depth calibration method |
| CN115096756A (en) * | 2022-07-27 | 2022-09-23 | 浙江双元科技股份有限公司 | Self-calibration surface density detector and calibration method thereof |
| CN115096756B (en) * | 2022-07-27 | 2022-11-22 | 浙江双元科技股份有限公司 | Self-calibration surface density detector and calibration method thereof |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105445802A (en) | High-precision time-depth correction method based on time-frequency spectrum constraints | |
| Munafò et al. | On the relationship between M w and ML for small earthquakes | |
| EP3028071B1 (en) | Method and device for the generation and application of anisotropic elastic parameters in horizontal transverse isotropic (hti) media | |
| CN105588883B (en) | Three-dimensional rock mechanics parameters acquisition methods and system | |
| CN106291677B (en) | A kind of poststack sound impedance inversion method based on match tracing method | |
| Witten et al. | Extended wave-equation imaging conditions for passive seismic data | |
| CN106291682B (en) | A kind of poststack sound impedance inversion method based on base method for tracing | |
| Shao et al. | What the exercise of the SPICE source inversion validation BlindTest 1 did not tell you | |
| CN104122581A (en) | Poststack acoustic wave impedance inversion method | |
| CN102156297A (en) | Fluid substitution method based on sandstone reservoir post-stack seismic data | |
| Stehly et al. | Towards improving ambient noise tomography using simultaneously curvelet denoising filters and SEM simulations of seismic ambient noise | |
| EP3749987A1 (en) | Systems and methods to enhance 3-d prestack seismic data based on non-linear beamforming in the cross-spread domain | |
| CN104007466B (en) | The reservoir that a kind of no restriction from borehole data prestack inversion based on P-wave amplitude realizes and fluid prediction method | |
| CN111208564A (en) | Depth domain horizon calibration method and device | |
| CN105301658A (en) | Method for extracting near-surface deconvolution operator by utilizing record of microlog | |
| CN102914797A (en) | Method and device for acquiring anisotropy coefficient of stratum | |
| CN102998700B (en) | A kind of earth's surface-consistent energy compensation process of vertical seismic profile data | |
| WO2011059977A3 (en) | System and method for analyzing and transforming geophysical and petrophysical data | |
| Westaway | Extrapolation of populations of small earthquakes to predict consequences of low-probability high impact events: The Pohang case study revisited | |
| WO2015124960A1 (en) | Systems and methods for improved inversion analysis of seismic data | |
| CN104345347A (en) | Logging curve recovery method used for compact gas-containing sandstone reservoir prediction | |
| Lin et al. | Semi-automation of borehole seismic travel-time picking by time-frequency analysis | |
| CN104316967B (en) | Large offseting distance NMO velocity Time Displacement Error compensation method | |
| Tan et al. | Combined adaptive multiple subtraction based on event tracing and Wiener filtering | |
| Gei et al. | Effective anisotropy parameters from microseismic reservoir monitoring |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20160330 |
|
| WD01 | Invention patent application deemed withdrawn after publication |