CN105717542B - 基于三维地震的Wheeler域切片砂体雕刻方法 - Google Patents
基于三维地震的Wheeler域切片砂体雕刻方法 Download PDFInfo
- Publication number
- CN105717542B CN105717542B CN201610220910.XA CN201610220910A CN105717542B CN 105717542 B CN105717542 B CN 105717542B CN 201610220910 A CN201610220910 A CN 201610220910A CN 105717542 B CN105717542 B CN 105717542B
- Authority
- CN
- China
- Prior art keywords
- domain
- wheeler
- seismics
- sand body
- sliced
- 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
- 239000004576 sand Substances 0.000 title claims abstract description 55
- 238000000034 method Methods 0.000 title claims abstract description 26
- 230000008569 process Effects 0.000 title claims abstract description 22
- 238000005553 drilling Methods 0.000 claims abstract description 11
- 230000008859 change Effects 0.000 claims abstract description 10
- 239000010410 layer Substances 0.000 claims description 37
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000008021 deposition Effects 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 3
- 239000002356 single layer Substances 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 description 3
- 230000000704 physical effect Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 2
- 238000012163 sequencing technique Methods 0.000 description 2
- 208000035126 Facies Diseases 0.000 description 1
- 241001074085 Scophthalmus aquosus Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 238000003786 synthesis reaction 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/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- 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/63—Seismic attributes, e.g. amplitude, polarity, instant phase
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
本发明公开了一种基于三维地震的Wheeler域切片砂体雕刻方法，包括以下步骤：井震结合确定层序顶底界面，计算地震同相轴倾角导向体，根据倾角导向是否突变追踪年代地层，年代地层按沉积顺序排列，代表相对地质时间。追踪年代地层的参数为，数据驱动，所有层都连续，最大迭代为2次；突变点打断，迭代10次，建立起年代地层三维层位体。将层位体的所有小层全部拉平，转换到Wheeler域。Wheeler域切片，能够反映砂体的平面分布范围，对砂体进行空间雕刻，Wheeler域地震属性切片，还能反映砂体内部地震属性的平面分布，研究砂体分布特征。该方法可以提高冲积扇砂体、扇三角洲砂体识别的准确性，提高钻井成功率。
Description
技术领域
本发明属于砂体雕刻技术领域，具体地涉及一种基于三维地震的Wheeler域切片砂体雕刻方法。
背景技术
冲积扇砂体、扇三角洲砂体识别的原理是基于砂体的密度和速度与周围的泥岩有较明显的差异，从而每一期的砂体反映在地震资料上有一定的特征，以往通过反演、地震属性分析方法对砂体的识别，高度依赖钻井资料，在离井较远的区域，误差很大，离井越远误差越大，钻井成功率仅为60-70%，因此还有30-40%的钻井成功率有待提高。
对于冲积扇和扇三角洲为储层的油藏，高效井一般只分布在主河道相带内，河道间为致密储层分布区，可作侧向遮挡，为了精细刻画冲积扇和扇三角洲的主河道砂体的分布，迫切需要砂体雕刻技术，以进一步提高钻探成功率。而当前国内外对于岩性横向变化大，储层非均质性强，但主河道内产量高的油藏，还缺乏利用地震资料研究砂体的有效手段。
发明内容
针对上述存在的技术问题，本发明旨在提供一种基于三维地震的Wheeler域切片砂体雕刻方法，可以精细刻画砂体的空间分布范围，以及砂体的储层物性特征，并且提高钻井成功率。
为达到上述目的，本发明的技术方案是：
一种基于三维地震的Wheeler域切片砂体雕刻方法，其特征在于，包括以下步骤：
S01：井震结合确定层序顶底界面；
S02：计算三维地震倾角导向体，通过自动检测三维地震同相轴，计算每个样点的倾角和方位角；
S03：计算三维地震年代地层层位体，所述层位体为一组可连续追踪的层位，按地层沉积顺序排列；
S04：将层位体的所有小层全部拉平，转到Wheeler域，再把地震数据体按小层拉平，生成Wheeler域的地震数据体；
S05：根据井震标定结果显示Wheeler域切片，确定砂体顶底平面分布范围，做出砂体顶底界面构造图；根据Wheeler域地震切片，分析砂体的地震属性特征，研究砂体的储层物性及含油性分布特征。
优选的，所述步骤S01中井震结合确定层序顶底界面，为等时层序界面，对目的层段，通过钻井合成地震记录，以及测井曲线精细层位标定，以单砂体规模的高精度层序格架进行标定。
优选的，所述步骤S02中计算三维地震倾角导向体的算法选择快速傅立叶变换法FFT，计算步长为3线、3道、2个采样点。
优选的，所述步骤S03的年代地层能够精确反映地震同相轴的形态。
优选的，所述步骤S03包括追踪层位体，所述追踪层位体包括：①追踪连续层位体，确保层位闭合；②根据倾角、方位角是否连续变化把连续的层位体打断，层位重合部分，只保留单层；③层位体补洞消除空白洞。
优选的，所述步骤S04中所述Wheeler域数据体在垂向上代表相对地质时间，表示沉积的先后顺序，最下面的地层先沉积，上面的地层后沉积。
与现有技术相比，本发明的有益效果是：
本发明基于三维地震的Wheeler域切片砂体雕刻方法可以实现砂体平面分布的空间形态以及分析砂体的储层物性分布特征的研究，结合冲积扇、扇三角洲所处的大地构造背景、物源方向、沉积中心分布特征，基准面变化规律，砂体的地震属性特点，开展地质综合研究，寻找有利储层发育区，部署勘探开发井位，可以大大提高钻井成功率，成功率由原来的60-70%提高到90%。本发明中的基于三维地震的Wheeler域切片砂体雕刻方法具有广泛的应用前景，且能够大幅度提高钻井成功率。
附图说明
图1为本发明基于三维地震的Wheeler域切片砂体雕刻方法的流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明了，下面结合具体实施方式以及附图，对本发明进一步详细说明。应该理解，这些描述只是示例性的，而并非要限制本发明的范围。此外，在以下说明中，省略了对公知结构和技术的描述，以避免不必要地混淆本发明的概念。
实施例：
图1为本发明的基于三维地震的Wheeler域切片砂体雕刻方法的流程图，如图1所示。
在步骤101，井震结合确定层序顶底界面，为等时层序界面，对目的层段，通过钻井合成地震记录，以及测井曲线精细层位标定，以单砂体规模的高精度层序格架进行标定，为砂体空间雕刻，提供可靠的依据。等时层序界面是计算层位体的终止层，仅仅是限定层位体的计算范围，对层位体的计算结果不产生任何影响。流程进入到步骤102。
在步骤102，所述步骤S02中计算三维地震倾角导向体，算法选择FFT（快速傅立叶变换法），计算地震倾角导向体有多种算法，FFT法更加精准的反映地震同相轴的空间展布形态，计算步长为332（3线、3道、2个采样点），这一组参数能够反映同相轴的变化趋势，又能兼顾砂体横向变化的细节。流程进入到步骤103。
在步骤103，计算三维地震年代地层层位体，所述层位体是一组可连续追踪的层位，按地层沉积顺序排列，层位体能够精确反映地震同相轴的特征。流程进入到步骤104。
在步骤104，追踪层位体的步骤包括：①追踪连续层位体（Continuouse vents）104-1，连续层位体的所有层位都连续分布，层位可重合，目的是确保层位闭合，追踪连续层位体的参数为，起始位置为厚度最大处的第一个采样点，年代地层的间距最大为4ms，检测同相轴的最远距离50道，最大迭代次数2；②把连续层位体打断（Truncate HorizonCube）104-2，根据倾角、方位角是否连续变化把连续的层位体打断，层位重合部分，只保留单层；③层位体补洞（Add iterations to HorizonCube）104-3，连续层位体打断后，层位重合现象被消除，但在剖面上会出现一些没有层位的空白洞，层位体补洞用于消除空白洞，两个年代地层间距最大2ms，超过2ms时增加一层，最大迭代次数10。
在步骤105，将层位体转到Wheeler域，把层位体的小层全部拉平，转到Wheeler域，Wheeler域数据体X、Y方向与常规数据体完全一致，在垂向上不代表深度也不代表时间，而是相对地质时间，表示沉积的先后顺序，最下面的地层先沉积，上面的地层后沉积，再把地震数据体按小层拉平，生成Wheeler域的各种地震属性体，包括常规数据体、反演数据体、相位数据体等。流程进入到步骤106。
在步骤106，根据井震标定结果显示Wheeler域切片，可以确定砂体顶底平面分布范围，做出砂体顶底界面构造图，同时，根据Wheeler域地震属性切片，分析砂体的地震属性特征，研究砂体的储层物性及含油性分布特征。
通过Wheeler域切片对砂体的精细雕刻，研究砂体平面分布形态，分析砂体物性特征，结合冲积扇、扇三角洲所处的大地构造背景，物源方向，沉积中心，基准面变化规律，砂体的地震属性特点，开展地质综合研究，寻找有利储层发育区，部署勘探开发井位，可以大大提高钻井成功率，成功率由原来的60-70%提高到90%。本发明中的基于三维地震的Wheeler域切片砂体雕刻方法具有广泛的应用前景，且能够大幅度提高钻井成功率。
应当理解的是，本发明的上述具体实施方式仅仅用于示例性说明或解释本发明的原理，而不构成对本发明的限制。因此，在不偏离本发明的精神和范围的情况下所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。此外，本发明所附权利要求旨在涵盖落入所附权利要求范围和边界、或者这种范围和边界的等同形式内的全部变化和修改例。
Claims (7)
1.一种基于三维地震的Wheeler域切片砂体雕刻方法，其特征在于，包括以下步骤：
S01：井震结合确定层序顶底界面；
S02：计算三维地震倾角导向体，通过自动检测三维地震同相轴，计算每个样点的倾角和方位角；
S03：计算三维地震年代地层层位体，所述层位体为一组可连续追踪的层位，按地层沉积顺序排列，追踪层位体；
所述追踪层位体包括以下步骤：①追踪连续层位体，确保层位闭合；②根据倾角、方位角是否连续变化把连续的层位体打断，层位重合部分，只保留单层；③层位体补洞消除空白洞；
S04：将层位体的所有小层全部拉平，转到Wheeler域，再把地震数据体按小层拉平，生成Wheeler域的地震数据体；
S05：根据井震标定结果显示Wheeler域切片，确定砂体顶底平面分布范围，做出砂体顶底界面构造图；根据Wheeler域切片，分析砂体的地震属性特征，研究砂体的储层物性及含油性分布特征。
2.根据权利要求1所述的基于三维地震的Wheeler域切片砂体雕刻方法，其特征在于，所述步骤S01中井震结合确定层序顶底界面，为等时层序界面，对目的层段，通过钻井合成地震记录，以及测井曲线精细层位标定，以单砂体规模的高精度层序格架进行标定。
3.根据权利要求1所述的基于三维地震的Wheeler域切片砂体雕刻方法，其特征在于，所述步骤S02中计算三维地震倾角导向体的算法选择快速傅立叶变换法FFT，计算步长为3线、3道、2个采样点。
4.根据权利要求1所述的基于三维地震的Wheeler域切片砂体雕刻方法，其特征在于，所述步骤S03的三维地震年代地层层位体能够精确反映地震同相轴的形态。
5.根据权利要求1所述的基于三维地震的Wheeler域切片砂体雕刻方法，其特征在于，所述追踪连续层位体的参数为，起始位置为厚度最大处的第一个采样点，年代地层的间距最大为4ms，检测同相轴的最远距离50道，最大迭代次数2。
6.根据权利要求1所述的基于三维地震的Wheeler域切片砂体雕刻方法，其特征在于，所述层位体补洞消除空白洞的参数为两个年代地层间距最大2ms，超过2ms时增加一层，最大迭代次数10。
7.根据权利要求1所述的基于三维地震的Wheeler域切片砂体雕刻方法，其特征在于，所述步骤S04中所述Wheeler域的地震数据体在垂向上代表相对地质时间，表示沉积的先后顺序，最下面的地层先沉积，上面的地层后沉积。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610220910.XA CN105717542B (zh) | 2016-04-11 | 2016-04-11 | 基于三维地震的Wheeler域切片砂体雕刻方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610220910.XA CN105717542B (zh) | 2016-04-11 | 2016-04-11 | 基于三维地震的Wheeler域切片砂体雕刻方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105717542A CN105717542A (zh) | 2016-06-29 |
| CN105717542B true CN105717542B (zh) | 2019-02-22 |
Family
ID=56159948
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610220910.XA Active CN105717542B (zh) | 2016-04-11 | 2016-04-11 | 基于三维地震的Wheeler域切片砂体雕刻方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105717542B (zh) |
Families Citing this family (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107976713B (zh) * | 2017-11-08 | 2019-07-09 | 中国石油天然气股份有限公司 | 一种高维地震数据输入下去除沉积背景的方法及装置 |
| CN111273346B (zh) * | 2018-12-04 | 2022-02-01 | 中国石油天然气股份有限公司 | 去除沉积背景的方法、装置、计算机设备及可读存储介质 |
| US11604909B2 (en) | 2019-05-28 | 2023-03-14 | Chevron U.S.A. Inc. | System and method for accelerated computation of subsurface representations |
| CN110244359B (zh) * | 2019-07-03 | 2022-02-25 | 成都理工大学 | 一种基于改进地震切片技术的异常体边缘检测计算方法 |
| US11249220B2 (en) | 2019-08-14 | 2022-02-15 | Chevron U.S.A. Inc. | Correlation matrix for simultaneously correlating multiple wells |
| CN112711069B (zh) * | 2019-10-24 | 2024-03-26 | 中国石油化工股份有限公司 | 基于平原区曲流河河道演化的点沙坝预测方法 |
| CN112904423B (zh) * | 2019-12-03 | 2024-04-30 | 中国石油天然气股份有限公司 | 陆相河道的雕刻方法及装置 |
| US11187826B2 (en) | 2019-12-06 | 2021-11-30 | Chevron U.S.A. Inc. | Characterization of subsurface regions using moving-window based analysis of unsegmented continuous data |
| US11010969B1 (en) | 2019-12-06 | 2021-05-18 | Chevron U.S.A. Inc. | Generation of subsurface representations using layer-space |
| US10984590B1 (en) | 2019-12-06 | 2021-04-20 | Chevron U.S.A. Inc. | Generation of subsurface representations using layer-space |
| US11320566B2 (en) | 2020-01-16 | 2022-05-03 | Chevron U.S.A. Inc. | Multiple well matching within subsurface representation |
| US11263362B2 (en) | 2020-01-16 | 2022-03-01 | Chevron U.S.A. Inc. | Correlation of multiple wells using subsurface representation |
| US11397279B2 (en) | 2020-03-27 | 2022-07-26 | Chevron U.S.A. Inc. | Comparison of wells using a dissimilarity matrix |
| CN113419274B (zh) * | 2021-06-21 | 2022-03-01 | 大庆油田有限责任公司 | 基于高精度层序格架模型的三维地震切片属性体提取方法 |
| CN116413794B (zh) * | 2021-12-31 | 2024-07-23 | 中国石油天然气股份有限公司 | 碳酸盐岩的沉积相类型的预测方法及装置 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| BRPI0822118A2 (pt) * | 2008-01-28 | 2015-06-23 | Landmark Graphics Corp | Método para modelar uma sequência na terra, programa de computação, e, interface gráfica de usuário |
| CN103454678B (zh) * | 2013-08-12 | 2016-02-10 | 中国石油天然气股份有限公司 | 一种地震切片等时性的确定方法及系统 |
| FR3016983A1 (fr) * | 2014-01-28 | 2015-07-31 | IFP Energies Nouvelles | Procede de construction d'un maillage volumique pour la modelisation de structures geologiques |
| CN105372705B (zh) * | 2015-10-27 | 2017-10-27 | 中国石油天然气股份有限公司 | 一种基于多波资料的地层切片方法 |
-
2016
- 2016-04-11 CN CN201610220910.XA patent/CN105717542B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105717542A (zh) | 2016-06-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105717542B (zh) | 基于三维地震的Wheeler域切片砂体雕刻方法 | |
| CN102053270B (zh) | 一种基于沉积地层单元的地震相分析方法 | |
| US7769546B2 (en) | Method for indexing a subsurface volume for the purpose of inferring geologic information | |
| CN103454678B (zh) | 一种地震切片等时性的确定方法及系统 | |
| EP3329307B1 (en) | Assignment of systems tracts | |
| US7769545B2 (en) | Method for determining geological information related to a subsurface volume of interest | |
| CN102200588B (zh) | 地震资料的波形相似性体曲率分析方法 | |
| US20130223187A1 (en) | Geological Structure Contour Modeling and Imaging | |
| NO20160203A1 (en) | Stratigraphic function | |
| GB2506026A (en) | Stratigraphic analysis of seismic data | |
| WO1997033184A1 (en) | Method for geophysical processing and interpretation using seismic trace difference for analysis and display | |
| CN109870719B (zh) | 一种碳酸盐岩致密薄储层的井位布设方法、装置及系统 | |
| EP2525242A2 (en) | Multiscale geologic modeling of a clastic meander belt including asymmetry using multi-point statistics | |
| EP1141750A1 (en) | Hydrocarbon edge detection using seismic amplitude | |
| Burton et al. | Geologically-based permeability anisotropy estimates for tidally-influenced reservoirs using quantitative shale data | |
| CN104636980A (zh) | 针对河道砂油藏类型油气汇集条件的地球物理表征方法 | |
| CN110231652A (zh) | 一种基于密度的含噪声应用空间聚类的地震相提取方法 | |
| CN102096098A (zh) | 地震属性提取方法 | |
| CN105093266A (zh) | 地震地质综合描述地层超覆线方法 | |
| Dando et al. | Complexity in microseismic phase identification: full waveform modelling, traveltime computations and implications for event locations within the Groningen gas field | |
| Anomneze et al. | Application of seismic stratigraphy and structural analysis in the determination of petroleum plays within the Eastern Niger Delta Basin, Nigeria | |
| WO2012060888A1 (en) | System and method for providing a physical property model | |
| Zhao et al. | Application of gradient structure tensor method in CBM fracture identification and sweet spot prediction | |
| CN113640875B (zh) | 一种沉积相识别方法、装置及系统 | |
| Zhai et al. | Quantitative seismic analysis of upper Neogene hierarchical fluvial reservoir architecture, Bohai Bay Basin |
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 |