CN109283581A - A kind of reservoir gas-bearing property evaluation method based on the analysis of Depth Domain seismic-wave dispersion - Google Patents
A kind of reservoir gas-bearing property evaluation method based on the analysis of Depth Domain seismic-wave dispersion Download PDFInfo
- Publication number
- CN109283581A CN109283581A CN201811388697.9A CN201811388697A CN109283581A CN 109283581 A CN109283581 A CN 109283581A CN 201811388697 A CN201811388697 A CN 201811388697A CN 109283581 A CN109283581 A CN 109283581A
- Authority
- CN
- China
- Prior art keywords
- depth
- domain
- analysis
- record
- log
- 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
- 239000006185 dispersion Substances 0.000 title claims abstract description 17
- 238000011156 evaluation Methods 0.000 title claims abstract description 11
- 238000005070 sampling Methods 0.000 claims description 11
- 238000012952 Resampling Methods 0.000 claims description 8
- 238000001228 spectrum Methods 0.000 claims description 6
- 238000000034 method Methods 0.000 claims description 5
- 239000002131 composite material Substances 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 3
- 230000007613 environmental effect Effects 0.000 claims description 3
- 239000000284 extract Substances 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 4
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000003786 synthesis reaction Methods 0.000 description 2
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/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
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)
Abstract
The invention discloses a kind of reservoir gas-bearing property evaluation methods based on the analysis of Depth Domain seismic-wave dispersion, belong to oil and gas exploration technical field, comprising the following steps: time domain earthquake record is converted to Depth Domain earthquake record by step 1.；Step 2. carries out deep frequency analysis in Depth Domain, obtains frequency content with the variation of depth；Step 3. calculates the frequency dispersion of each depth, identifies gas-bearing layer.The present invention carries out the dispersion analysis of seismic wave in Depth Domain, then can directly determine the spatial depth position of gas-bearing layer, and the bulk effect of energy quantitative estimation reservoir.
Description
Technical field
The invention belongs to oil and gas exploration technical fields, are related to a kind of storage based on the analysis of Depth Domain seismic-wave dispersion
Layer gassiness evaluation method.
Background technique
The existing reservoir gas-bearing property evaluation method technology based on seismic-wave dispersion analysis carries out (time frequency analysis) in time-domain,
Because of the changeability of geologic medium seismic wave propagation speed, do not known so that the spatial position of the gas-bearing layer of identification determines to exist
Property, while can not also determine bulk effect.
Summary of the invention
The purpose of the present invention is to provide a kind of reservoir gas-bearing property evaluation methods based on the analysis of Depth Domain seismic-wave dispersion.
To realize above-mentioned technical purpose and the technique effect, its technical solution specifically:
Include the following steps；
Time domain earthquake record is converted to Depth Domain earthquake record by step 1.；Layer thickness division is carried out using log and is built
Vertical depth-when relationship seeks the reflection coefficient of frequency domain according to the wave equation derivation formula of multilayer dielectricity reflection coefficient；Clock synchronization
Between domain well bypass road record using homomorphic deconvolution method extract seismic wavelet, obtain the wavelet spectrum of frequency domain；Utilize frequency domain
Convolution model calculates composite traces, obtains time-domain synthetic seismogram as inversefouriertransform, finds out further according to time-depth relationship
The corresponding depth point of each time domain sampling point carries out resampling according to required depth sampling interval, obtains Depth Domain synthesis
Earthquake record；
Step 2. carries out deep frequency analysis in Depth Domain, obtains frequency content with the variation of depth；
Depth Domain earthquake record is pre-processed.The pretreatment of log includes going singular value, environmental correction and adopting again
Sample etc. carries out wavelet transformation and spectrum sigtral response in Depth Domain to pretreated data (signal), obtain frequency content with
The variation of depth；
Step 3. calculates the frequency dispersion of each depth, identifies gas-bearing layer.
Further, when deep and time and depth transfer is during seeking seismic wavelet, and well bypass road is needed to record from Depth Domain
It is transformed into time-domain, the serialization formula of discrete curve can be used, the record value of corresponding log depth point is sought, then weighs
Sampling obtains the record of isochronic interval；Or the layering interval velocity according to each log, point-by-point time value is sought, is then weighed
Sampling just obtains the record of isochronic interval.
Further, in step 1, the log is sound wave, density log curve.
Further, in step 1, the log is equal interval sampling.
The invention has the following advantages: the present invention carries out the dispersion analysis of seismic wave in Depth Domain, then it can directly really
Determine the spatial depth position of gas-bearing layer, and the bulk effect of energy quantitative estimation reservoir.
Specific embodiment
In order to make the objectives, technical solutions, and advantages of the present invention clearer, with reference to embodiments, to the present invention
It is further elaborated.
A kind of reservoir gas-bearing property evaluation method based on the analysis of Depth Domain seismic-wave dispersion, comprising the following steps:
Time domain earthquake record is converted to Depth Domain earthquake record by step 1.；Layer thickness division is carried out using log and is built
Vertical depth-when relationship seeks the reflection coefficient of frequency domain according to the wave equation derivation formula of multilayer dielectricity reflection coefficient；Clock synchronization
Between domain well bypass road record using homomorphic deconvolution method extract seismic wavelet, obtain the wavelet spectrum of frequency domain；Utilize frequency domain
Convolution model calculates composite traces, obtains time-domain synthetic seismogram as inversefouriertransform, finds out further according to time-depth relationship
The corresponding depth point of each time domain sampling point carries out resampling according to required depth sampling interval, obtains Depth Domain synthesis
Earthquake record；
Step 2. carries out deep frequency analysis in Depth Domain, obtains frequency content with the variation of depth；
Depth Domain earthquake record is pre-processed.The pretreatment of log includes going singular value, environmental correction and adopting again
Sample etc. carries out wavelet transformation and spectrum sigtral response in Depth Domain to pretreated data (signal), obtain frequency content with
The variation of depth；
Step 3. calculates the frequency dispersion of each depth, identifies gas-bearing layer.
When deep and time and depth transfer is during seeking seismic wavelet, when well bypass road record being needed to be transformed into from Depth Domain
Between domain, the serialization formula of discrete curve can be used, seek the record value of corresponding log depth point, then resampling obtains
The record of isochronic interval；Or the layering interval velocity according to each log, point-by-point time value is sought, then resampling must
To the record of isochronic interval.
In step 1, the log is sound wave, density log curve.
In step 1, the log is equal interval sampling.
The foregoing is only a preferred embodiment of the present invention, the scope of protection of the present invention is not limited to this, it is any ripe
Know those skilled in the art within the technical scope of the present disclosure, the letter for the technical solution that can be become apparent to
Altered or equivalence replacement are fallen within the protection scope of the present invention.
Claims (4)
1. a kind of reservoir gas-bearing property evaluation method based on the analysis of Depth Domain seismic-wave dispersion, which is characterized in that including following step
Suddenly；
Time domain earthquake record is converted to Depth Domain earthquake record by step 1.；Layer thickness division is carried out using log and is established
Depth-when, relationship according to the wave equation derivation formula of multilayer dielectricity reflection coefficient sought the reflection coefficient of frequency domain；To the time
The well bypass road record in domain extracts seismic wavelet using homomorphic deconvolution method, obtains the wavelet spectrum of frequency domain；Utilize frequency domain pleat
Product module type calculates composite traces, obtains time-domain synthetic seismogram as inversefouriertransform, finds out often further according to time-depth relationship
The corresponding depth point of a time domain sampling point carries out resampling according to required depth sampling interval, obtains Depth Domain synthetically
Shake record；
Step 2. carries out deep frequency analysis in Depth Domain, obtains frequency content with the variation of depth；
Depth Domain earthquake record is pre-processed, and the pretreatment of log includes going singular value, environmental correction and resampling, right
Pretreated data-signal carries out wavelet transformation and spectrum sigtral response in Depth Domain, obtains frequency content with the change of depth
Change；
Step 3. calculates the frequency dispersion of each depth, identifies gas-bearing layer.
2. a kind of reservoir gas-bearing property evaluation method based on the analysis of Depth Domain seismic-wave dispersion as described in claim 1, special
Sign is that when deep and time and depth transfer is during seeking seismic wavelet, when well bypass road record being needed to be transformed into from Depth Domain
Between domain the record value of corresponding log depth point is sought, whens then resampling obtains equal with the serialization formula of discrete curve
The record at interval；Or the layering interval velocity according to each log, point-by-point time value is sought, then resampling just obtains
When interval record.
3. a kind of reservoir gas-bearing property evaluation method based on the analysis of Depth Domain seismic-wave dispersion as described in claim 1, special
Sign is, in step 1, the log is sound wave, density log curve.
4. a kind of reservoir gas-bearing property evaluation method based on the analysis of Depth Domain seismic-wave dispersion as described in claim 1, special
Sign is, in step 1, the log is equal interval sampling.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811388697.9A CN109283581B (en) | 2018-11-20 | 2018-11-20 | Reservoir gas content evaluation method based on depth domain seismic wave frequency dispersion analysis |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811388697.9A CN109283581B (en) | 2018-11-20 | 2018-11-20 | Reservoir gas content evaluation method based on depth domain seismic wave frequency dispersion analysis |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109283581A true CN109283581A (en) | 2019-01-29 |
| CN109283581B CN109283581B (en) | 2020-08-07 |
Family
ID=65176303
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811388697.9A Expired - Fee Related CN109283581B (en) | 2018-11-20 | 2018-11-20 | Reservoir gas content evaluation method based on depth domain seismic wave frequency dispersion analysis |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109283581B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111060970A (en) * | 2019-12-30 | 2020-04-24 | 恒泰艾普（北京）能源科技研究院有限公司 | Reservoir constrained deconvolution method and application |
| CN112099086A (en) * | 2020-09-16 | 2020-12-18 | 中油奥博(成都)科技有限公司 | High-resolution optical fiber borehole seismic data deep frequency analysis method |
| US20220221614A1 (en) * | 2021-01-11 | 2022-07-14 | Shandong University Of Science And Technology | Analysis method, system and storage media of lithological and oil and gas containing properties of reservoirs |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106291689A (en) * | 2016-08-18 | 2017-01-04 | 中国石油天然气集团公司 | A kind of extract the processing method of geological data frequency dispersion attribute, device and prognoses system |
| CN107255831A (en) * | 2017-06-14 | 2017-10-17 | 中国石油化工股份有限公司 | A kind of extracting method of prestack frequency dispersion attribute |
-
2018
- 2018-11-20 CN CN201811388697.9A patent/CN109283581B/en not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106291689A (en) * | 2016-08-18 | 2017-01-04 | 中国石油天然气集团公司 | A kind of extract the processing method of geological data frequency dispersion attribute, device and prognoses system |
| CN107255831A (en) * | 2017-06-14 | 2017-10-17 | 中国石油化工股份有限公司 | A kind of extracting method of prestack frequency dispersion attribute |
Non-Patent Citations (3)
| Title |
|---|
| 徐敬领 等: "利用小波深频分析方法研究沉积储层旋回", 《中国石油大学学报（自然科学版）》 * |
| 杨国权 等: "井间地震资料层位标定方法研究", 《石油物探》 * |
| 胡玮 等: "频散属性在致密薄储层识别中的应用", 《石油物探》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111060970A (en) * | 2019-12-30 | 2020-04-24 | 恒泰艾普（北京）能源科技研究院有限公司 | Reservoir constrained deconvolution method and application |
| CN112099086A (en) * | 2020-09-16 | 2020-12-18 | 中油奥博(成都)科技有限公司 | High-resolution optical fiber borehole seismic data deep frequency analysis method |
| CN112099086B (en) * | 2020-09-16 | 2022-03-29 | 中油奥博(成都)科技有限公司 | High-resolution optical fiber borehole seismic data deep frequency analysis method |
| US20220221614A1 (en) * | 2021-01-11 | 2022-07-14 | Shandong University Of Science And Technology | Analysis method, system and storage media of lithological and oil and gas containing properties of reservoirs |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109283581B (en) | 2020-08-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Lu et al. | Seismic spectral decomposition using deconvolutive short-time Fourier transform spectrogram | |
| CN103842853B (en) | Time shift geological data is corrected for coating and record effect | |
| Yanhu et al. | A method of seismic meme inversion and its application | |
| US7525874B2 (en) | Method for converting seismic data from the time domain to the depth domain | |
| CN109283581A (en) | A kind of reservoir gas-bearing property evaluation method based on the analysis of Depth Domain seismic-wave dispersion | |
| CN104570067B (en) | Phased seismic inversion method in a kind of geophysical exploration | |
| AU2014273165A1 (en) | High resolution estimation of attenuation from vertical seismic profiles | |
| CN105425289B (en) | The method and apparatus for determining low frequency wave impedance | |
| CN102937720A (en) | Method for improving seismic data resolution ratio for well control | |
| CN103852785B (en) | The evaluation method of formation anisotropy | |
| Ringler et al. | A comparison of seismic instrument noise coherence analysis techniques | |
| CN103064112A (en) | Absorption method used for detecting oil-gas possibility of stratum | |
| CN107436451A (en) | A kind of automatic amplitude spectral method for calculating geological data optical cable coupled noise degree of strength | |
| Xue et al. | Q-factor estimation by compensation of amplitude spectra in synchrosqueezed wavelet domain | |
| CN109581500B (en) | Reflection seismic record frequency-variable velocity analysis method | |
| Cheng et al. | A new method for estimating the correlation of seismic waveforms based on the NTFT | |
| Gong et al. | Automatic time picking of microseismic data based on shearlet-AIC algorithm | |
| CN109425892A (en) | The estimation method and system of seismic wavelet | |
| Zhao et al. | A general method for extracting the amplitude spectrum of the seismic wavelet from the seismic traces | |
| Feng et al. | Dispersion calculation method based on S-transform and coordinate rotation for Love channel waves with two components | |
| Qi et al. | Well ties for seismic with severe stratigraphic filtering | |
| CN102353991A (en) | Method for analyzing seismic instantaneous frequency based on physical wavelet matched with seismic wavelet | |
| CN108107482B (en) | Hydrocarbon source rock multiparameter prediction method and system under a kind of condition of less wells | |
| Xu et al. | Time domain elastic-wave full waveform inversion based on first-order approximate instantaneous frequency | |
| Zhang et al. | Fine complex geological structure interpretation based on multiscale seismic dip constraint |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: He Xiaoyan Inventor after: Cao Junxing Inventor after: Wang Shuyi Inventor after: Jiang Xudong Inventor before: Cao Junxing Inventor before: He Xiaoyan Inventor before: Wang Shuyi Inventor before: Jiang Xudong |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200807 |