CN109765621B - 一种沉积相精细刻画的成图方法 - Google Patents
一种沉积相精细刻画的成图方法 Download PDFInfo
- Publication number
- CN109765621B CN109765621B CN201910014721.0A CN201910014721A CN109765621B CN 109765621 B CN109765621 B CN 109765621B CN 201910014721 A CN201910014721 A CN 201910014721A CN 109765621 B CN109765621 B CN 109765621B
- Authority
- CN
- China
- Prior art keywords
- sedimentary
- proportion
- well
- seismic
- thickness
- 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
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种沉积相精细刻画的成图方法。它包括以下步骤：针对目标储层，在地震剖面上根据地震反射轴提取地震属性，形成粗化的沉积相平面图；在井点依据沉积特性和测井形态识别出各个沉积期次，并通过其计算各个沉积期次实际钻井上的比例；对目标储层，利用实际钻井的比例，进行井间插值得到每一个期次的平面厚度比例，按照沉积期次的比例将地震轴进行劈分，提取地震属性，实现粗化的沉积相平面图在垂向储层上的划分；将沉积体各期期次在剖面上和平面上的分布状况进行地震属性提取，并对每种地震属性进行沉积相刻画，雕刻出几个沉积相，即得到沉积相细分的成图。本发明为井网密度稀疏、低于地震分辨率的沉积相精细刻画能应用于油田勘探开发中。
Description
技术领域
本发明涉及一种沉积相精细刻画的成图方法，属于石油天然气开采领域。
背景技术
精细的沉积相平面图是进行油田生产规划、开发过程的动态分析和调整的重要基础。由于沉积体的旋回性，导致沉积相带侧向相变变化大和内部结构复杂，这就必然要求对沉积相的刻画从垂向上和平面上做的精细，从而准确提高井位部署、深挖剩余油和规避含水上升等问题。目前，精细沉积相的描述主要是基于三维地震信息和井信息。然而由于受地震资料分辨能力的制约以及岩性界面的穿时性，传统的地震相分析法不能满足油藏开发地质研究中对储层内部结构精细解释的需要。针对井信息的沉积相刻画，利用多井信息进行沉积相图编制，这种方法无法实现井位之间的砂体预测，只能采用插值，无法准确预测储集体的空间分布形态和规模，难于满足油田开发中后期剩余油挖潜和井位部署的需要。如何精细刻画储层沉积相，以及科学的确定海上油田部井方案，是在油田开发阶段提高油气采收率、井位调整和部署的关键问题。
发明内容
本发明的目的是提供一种沉积相精细刻画的成图方法，本发明为井网密度稀疏、低于地震分辨率的沉积相精细刻画，能广泛应用于油田勘探开发评价领域中。
本发明提供的一种沉积相细分的成图方法，它包括以下步骤：
1)针对目标储层，在地震剖面上根据地震反射轴提取地震属性，形成粗化的沉积相平面图；
2)对步骤1)处理后的所述目标储层，在井点依据沉积特性和测井形态识别出各个沉积期次，并通过对所述沉积期次划分，计算各个所述沉积期次实际钻井上的比例；
3)对所述目标储层，利用2)计算的实际钻井的比例，进行井间插值得到每一个期次的平面厚度比例，按照所述沉积期次的比例将所述地震轴进行劈分，提取地震属性，实现所述粗化的沉积相平面图在垂向储层上的划分；
4)将步骤3)得到的沉积体各期期次在剖面上和平面上的分布状况进行地震属性提取，并对每种地震属性进行沉积相刻画，雕刻出几个沉积相，即得到沉积相细分的成图。
本发明中，步骤1)中的粗化的沉积相进行了精细刻画，从而实现了沉积相的精细刻画，满足了开发中后期井位调整的目的，为剩余油的挖潜和水体的规避起到了至关重要的作用。
上述的方法中，所述目标储层的厚度为30～50m。
上述的方法步骤2)中识别所述沉积期次与划分的步骤如下：井点依据所述沉积特征，整体上期次划分；结合每一口钻井的测井形态，按照测井曲线的回返，针对每口井统一划分相同的数量的期次；
所述地震轴内的沉积期次划分不超过4～5期。
上述的方法步骤2)中对所述沉积期次划分和计算各个所述沉积期次的比例的方法，其步骤如下：整个地震同相轴时间厚度乘以每一个期次的厚度比例，得到每一个期次的时间厚度，生成不同期次的时间顶底层位，从而实现对沉积体进行期次精细划分。
上述的方法步骤2)中对所述沉积期次划分和计算各个所述沉积期次的比例的方法，其步骤如下：a)所述目标储层内有m口井点，划分成n个期次，第n个期次在第m口井的的厚度比例为hnm/Hm，其中，hnm为第n期在m井上的厚度比例，Hm为m井上该目标储层的厚度；
b)在平面上对每一个期次的厚度比例进行平面插值得到每一个期次的厚度比例平面图，每一个期次的厚度比例按照式Ⅰ计算得到；
式Ⅰ中，qn是第n期的厚度比例，wj(j＝1,…,m)是井点的权重系数，hnj/Hj为第n个期次在第j口井的的厚度比例。
本发明具有以下优点：
1、本发明提高了地震的使用精度，利用井资料对其进行合理的解释，从而实现了提高了地震分辨率。沉积体受地震垂向分辨率限制而无法分辨，然而井上的期次可以识别划分，通过对平面上钻井期次进行统计分析，从而实现了对一个地震轴的合理划分。层次清楚，能够对每一目标储层进行精细刻画。
2、充分挖掘地震和钻井信息，使的信息更为可靠。
3、便于操作，可以广泛应用于油田勘探开发评价领域中。
附图说明
图1是本发明某一目标储层的地震属性(a)和沉积相(b)的示意图。
图2是本发明实际钻井储层划分的示意图。
图3是本发明目标储层依据井的比例劈分结果的示意图。
图4是本发明目标储层地震属性结果的示意图，其中图4(a)、(b)分别为前后的地震属性结果。
图5是本发明目标储层沉积相早期和晚期的精细刻画的示意图。
具体实施方式
下述实施例中所使用的实验方法如无特殊说明，均为常规方法。
下面结合附图和实施案例对本发明进行详细的描述。
如图1～图3所示，本发明沉积相精细刻画的成图方法，包括以下步骤：
1)如图1所示，某一目标储层的地震属性平面图。地震剖面上表现为一个地震反射轴，地层厚度较厚，储层厚度约为30-50m，沉积相图较粗不能满足井位调整需要。
2)在井上进行储层期次划分，得到每一个期次的厚度比例平面图，如图2～图3所示。
由于沉积的旋回性，不同井点尽管存在相变，但是每个期次在井点均有所体现，因而可以在井点依据沉积特性和测井形态可以识别出各个期次。通过对实际钻井的期次精细划分，计算各个期次比例，按照这个比例关系将地震轴进行劈分，从而实现垂向储层细分。具体的，包括以下步骤：
①钻井期次识别与划分。
依据沉积特征，整体上期次划分。结合每一口钻井的测井箱形形态变化，按照测井曲线的回返，针对每口井统一划分相同的数量的期次。一般而言针对河流三角洲沉积相，一个地震轴内的沉积期次划分不超过4～5期。
②钻井期次计算方法。
一个地震同向轴的时间厚度大约为60～80ms，以采样率为1ms计算，则样点数一般为60～80个。工区内有m口井，每口井按照步骤①分成n个期次。计算各井不同期次在该井的厚度比例系数，第1期次在第1口井的厚度比例为h11/H1，第n个期次在第1口井的厚度比例为hn1/H1，第1期次在第2口井的厚度比例为h12/H2，第n个期次在第2口井的厚度比例为hn2/H2，第1期次在第m口井的厚度比例为h1m/Hm，第n个期次在第m口井的的厚度比例为hnm/Hm，即第1期次在各井点的厚度比例为
h11/H1，h12/H2，……，h1m/Hm
第2期次在各井点的厚度比例为
h21/H1，h22/H2，……，h2m/Hm
第n期次在各井点的厚度比例为
hn1/H1，hn2/H2，……，hnm/Hm
其中，hnm为第n期在m井上的厚度比例，Hm为m井上该目标储层的厚度。
在平面上对每一个期次的厚度比例进行插值得到每一个期次的厚度比例平面图，即
其中qn是第n期的厚度比例，wj(j＝1，…，m)是井点的权重系数。
3)根据步骤2)②计算的比例劈分地震反射轴，得到对应地震属性所包含地层的各个期次比例关系，如图3所示。
在地震同向反射轴上下地层的约束下，按照钻井的比例和平面插值的比例，对沉积体进行期次划分，从而得到沉积体各期期次在剖面上和平面上的分布状况。
整个地震同相轴时间厚度乘以每一个期次的厚度比例，得到每一个期次的时间厚度，生成不同期次的时间顶底层位，进而可以对细分后的期次进行地震属性提取，能精细刻画小尺度的储层横向变化，满足开发中后期的生产需求。
4)根据步骤3)提取相应地震属性平面图，按照地震属性进行沉积相精细刻画，如图3～图5所示。
对1)中的粗化的沉积相进行了精细刻画，雕刻出几个沉积相，从而实现了沉积相的精细刻画，满足了开发中后期井位调整的目的，为剩余油的挖潜和水体的规避起到了至关重要的作用。
Claims (4)
1.一种沉积相细分的成图方法，它包括以下步骤：
1)针对目标储层，在地震剖面上根据地震反射轴提取地震属性，形成粗化的沉积相平面图；
2)对步骤1)处理后的所述目标储层，在井点依据沉积特性和测井形态识别出各个沉积期次，并通过对所述沉积期次划分，计算各个所述沉积期次实际钻井上的比例；
步骤2)中对所述沉积期次划分和计算各个所述沉积期次的比例的方法，其步骤如下：a)所述目标储层内有m口井点，划分成n个期次，第n个期次在第m口井的的厚度比例为hnm/Hm，其中，hnm为第n期在m井上的厚度比例，Hm为m井上该目标储层的厚度；
b)在平面上对每一个期次的厚度比例进行平面插值得到每一个期次的厚度比例平面图，每一个期次的厚度比例按照式Ⅰ计算得到；
式Ⅰ中，qn是第n期的厚度比例，wj(j＝1,…,m)是井点的权重系数，hnj/Hj为第n个期次在第j口井的的厚度比例；
3)对所述目标储层，利用2)计算的实际钻井的比例，进行井间插值得到每一个期次的平面厚度比例，按照所述沉积期次的比例将所述地震反射轴进行劈分，提取地震属性，实现所述粗化的沉积相平面图在垂向储层上的划分；
4)将步骤3)得到的沉积体各期期次在剖面上和平面上的分布状况进行地震属性提取，并对每种地震属性进行沉积相刻画，雕刻出几个沉积相，即得到沉积相细分的成图。
2.根据权利要求1所述的方法，其特征在于：所述目标储层的厚度为30～50m。
3.根据权利要求1或2所述的方法，其特征在于：步骤2)中识别所述沉积期次与划分的步骤如下：井点依据所述沉积特征，整体上期次划分；结合每一口钻井的测井形态，按照测井曲线的回返，针对每口井统一划分相同的数量的期次；
所述地震反射轴内的沉积期次划分不超过4～5期。
4.根据权利要求1或2所述的方法，其特征在于：步骤2)中对所述沉积期次划分和计算各个所述沉积期次的比例的方法，其步骤如下：整个地震同相轴时间厚度乘以每一个期次的厚度比例，得到每一个期次的时间厚度，生成不同期次的时间顶底层位，从而实现对沉积体进行期次精细划分。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910014721.0A CN109765621B (zh) | 2019-01-08 | 2019-01-08 | 一种沉积相精细刻画的成图方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910014721.0A CN109765621B (zh) | 2019-01-08 | 2019-01-08 | 一种沉积相精细刻画的成图方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109765621A CN109765621A (zh) | 2019-05-17 |
| CN109765621B true CN109765621B (zh) | 2020-06-30 |
Family
ID=66453468
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910014721.0A Active CN109765621B (zh) | 2019-01-08 | 2019-01-08 | 一种沉积相精细刻画的成图方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109765621B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113960694B (zh) * | 2020-07-21 | 2024-02-20 | 中国石油化工股份有限公司 | 一种沉积相识别方法及装置 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104502969A (zh) * | 2014-12-30 | 2015-04-08 | 中国石油化工股份有限公司 | 河道砂岩性油藏的识别方法 |
| CN105629325A (zh) * | 2015-12-23 | 2016-06-01 | 中国石油大学(北京) | 前陆盆地冲积扇精细刻画与预测方法 |
| CN106597543A (zh) * | 2016-11-21 | 2017-04-26 | 中国石油天然气股份有限公司 | 一种地层沉积相划分方法 |
| CN107991716A (zh) * | 2017-10-13 | 2018-05-04 | 中国石油天然气股份有限公司 | 薄层的沉积相图和沉积体厚度的确定方法和装置 |
| CN108008466A (zh) * | 2017-11-07 | 2018-05-08 | 中国石油天然气股份有限公司 | 一种基于砂地比的复合河道沉积期次空间解析方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070276604A1 (en) * | 2006-05-25 | 2007-11-29 | Williams Ralph A | Method of locating oil and gas exploration prospects by data visualization and organization |
| US10120092B2 (en) * | 2016-08-22 | 2018-11-06 | Chevron U.S.A. Inc. | System and method for analysis of depositional settings of subsurface reservoirs |
-
2019
- 2019-01-08 CN CN201910014721.0A patent/CN109765621B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104502969A (zh) * | 2014-12-30 | 2015-04-08 | 中国石油化工股份有限公司 | 河道砂岩性油藏的识别方法 |
| CN105629325A (zh) * | 2015-12-23 | 2016-06-01 | 中国石油大学(北京) | 前陆盆地冲积扇精细刻画与预测方法 |
| CN106597543A (zh) * | 2016-11-21 | 2017-04-26 | 中国石油天然气股份有限公司 | 一种地层沉积相划分方法 |
| CN107991716A (zh) * | 2017-10-13 | 2018-05-04 | 中国石油天然气股份有限公司 | 薄层的沉积相图和沉积体厚度的确定方法和装置 |
| CN108008466A (zh) * | 2017-11-07 | 2018-05-08 | 中国石油天然气股份有限公司 | 一种基于砂地比的复合河道沉积期次空间解析方法 |
Non-Patent Citations (2)
| Title |
|---|
| 川西坳陷陆相致密气藏河道砂岩储层精细刻画技术及其应用;武恒志 等;《石油与天然气地质》;20150430;第36卷(第2期);第230-239页 * |
| 永安镇地区沙四段砂砾岩体沉积特征及分布规律;杜亮慧;《中国优秀硕士论文全文数据库-基础科学辑》;20180615(第06期);摘要及正文第15-16、19-30、60-64 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109765621A (zh) | 2019-05-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101158724B (zh) | 基于偶极小波的储层厚度预测方法 | |
| CN100437147C (zh) | 多参数降维油气水层识别方法 | |
| CN105528522A (zh) | 一种基于空间网格计算连续性油气藏资源量的方法及装置 | |
| CN105629304B (zh) | 一种基于多属性的砂体叠合模式识别方法 | |
| CN110838175A (zh) | 注气开发油藏的地质模型建立方法 | |
| Harris | The role of geology in reservoir simulation studies | |
| CN111239815A (zh) | 基于三维地震属性的砂岩型铀储层成矿沉积要素提取方法 | |
| CN106526697A (zh) | 一种三角洲建设与破坏主要转换面识别方法 | |
| CN109765621B (zh) | 一种沉积相精细刻画的成图方法 | |
| Iltaf et al. | Facies and petrophysical modeling of Triassic Chang 6 tight sandstone reservoir, Heshui oil field, Ordos basin, China | |
| Reda et al. | Hydrocarbon reservoir characterization in the challenging structural setting of southern gulf of Suez: synergistic approach of well log analyses and 2D seismic data interpretation | |
| Sunmonu et al. | Reservoir characterization and by-passed pay analysis of Philus Field in Niger Delta, Nigeria | |
| Chihi et al. | Multiscale modeling of reservoir systems using geostatistical methods | |
| Peijs-van Hilten et al. | Heterogeneity modeling and geopseudo upscaling applied to waterflood performance prediction of an incised valley reservoir: Countess YY Pool, southern Alberta, Canada | |
| Abe et al. | Structure and facies development resulting from neogene gravity tectonics and depositional processes: Application to Afo Field Niger Delta, Nigeria | |
| REN et al. | Geomodeling technology under multifactor control | |
| Zene et al. | Volumetric estimation and OOIP calculation of the Ronier4 block of Ronier oilfield in the Bongor basin, Chad | |
| Zhang et al. | Turbidite facies identification and analysis using geophysical methods in the salt diapir domain of Lower Congo basin | |
| Ivlev | Prediction of geophysical properties of rocks on rare well data and attributes of seismic waves by machine learning methods on the example of the Achimov formation | |
| Furnée | Geo-modeling and fluid flow simulation in paleokarst reservoirs | |
| Oyedele et al. | Reservoir characterization and volumetric analysis of roan field, Niger Delta, Nigeria | |
| Bu et al. | Geological modeling of combined incised channel system in deep-water gravity flow sedimentation: the case of OML130 in Niger delta basin, West Africa | |
| Hamdulla et al. | Building geological model for tertiary reservoir of exploration Ismail oil field, North Iraq | |
| Noad | The use of field analogues in the correlation and static reservoir methodology used in the Tern Field, Northern North Sea, UK | |
| Olowokere | Geostatistical modeling of interval velocity to quantifying hydrocarbon resource in multi-layer reservoir from TMB field, Niger delta |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |