CN105652323B - 一种储层预测方法 - Google Patents
一种储层预测方法 Download PDFInfo
- Publication number
- CN105652323B CN105652323B CN201510490804.9A CN201510490804A CN105652323B CN 105652323 B CN105652323 B CN 105652323B CN 201510490804 A CN201510490804 A CN 201510490804A CN 105652323 B CN105652323 B CN 105652323B
- Authority
- CN
- China
- Prior art keywords
- reservoir
- parameter
- log
- curve
- carried out
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种储层预测方法,通过对研究区内的测井数据进行交汇分析,确定能够区分储层的有效参数;综合利用地震资料和测井数据,对获得的每个有效参数分别进行正演模拟,并对应进行特征响应分析,确定能够区分储层的最优参数;对最优参数进行储层参数反演,根据反演结果确定合理的储层预测结果或将正演模拟结果和反演结果进行对比分析后确定储层预测结果。本发明的方法提高了预测的精度,降低了勘探风险,减少了工作量,是岩性油气藏储层预测的一种有效方案。
Description
技术领域
本发明涉及一种储层预测方法,属于地球物理勘探领域。
背景技术
在地球物理勘探技术中,常规的储层预测流程中需要对测井资料进行分析,对地震资料进行属性提取、地震反演等相关计算,然后对计算得到的属性数据进行分析,识别出储层发育的可能区域和层段,从而完成储层预测工作。
现有的储层预测技术对于特征明显的简单的并且比较厚的储层,效果比较明显,即当交汇图上的储层特征能以独立的区域相应特征反应出来,能够在相对的地震属性数据上圈定出特征区域的分布范围;但对于相对比较薄的储层,在现有储层预测技术中,尽管交汇分析能得到有效参数,但由于分辨率的问题,往往得到的储层预测结果不满意,则需要反复尝试不同的参数反演,操作比较盲目,交汇参数的选择和反演结果没有直观的相对关系,这就大大增加了分析的工作量,而且很难得到满意的结果。因此如何选择最优的反演参数就成为储层预测中比较重要的问题。
发明内容
本发明提供一种储层预测方法,旨在解决当得到的储层预测结果不满意时如何选择最优的反演参数进行储层预测。
为解决上述技术问题,本发明一种储层预测方法的技术方案包含以下步骤:
1)收集研究区内的地震数据以及相应的测井数据;
2)对测井数据进行交汇分析,确定能够区分储层的有效参数;
3)利用测井数据和地震数据,对每一个有效参数分别进行正演模拟,并对应进行特征响应分析,确定能够区分储层的最优参数;
4)利用测井数据和地震数据,对所述最优参数进行储层参数反演,根据反演结果,最终确定储层预测结果。
对所述正演模拟结果和反演结果进行对比分析,最终确定储层预测结果。
步骤1)中所述的测井数据包括测井曲线,测井曲线包括声波曲线、电阻率曲线、自然电位曲线、孔隙度曲线、饱和度曲线、密度曲线、伽马曲线和泥质含量曲线。
步骤2)包含如下步骤:
b1、对所述测井曲线进行横波预测,得到横波曲线;
b2、利用纵波曲线、横波曲线、密度曲线,进行弹性参数计算;
b3、对所述弹性参数进行交汇分析,确定能够区分储层的有效参数。
步骤3)包含如下步骤:
c1、提取地震子波或制作理论地震子波;
c2、利用所述测井曲线制作合成地震记录,并将合成地震记录与过井地震剖面进行对比,确定地震剖面与所述测井曲线之间的时深关系;
c3、导入地震解释层位和断层,依据地震解释层位和断层建立正演模型;
c4、依据建立的正演模型,对每一个有效参数分别进行正演模拟,并对应进行特征响应分析,确定能够区分储层的最优参数。
步骤4)包含如下步骤:
d1、提取地震子波或制作理论地震子波;
d2、利用所述测井曲线制作合成地震记录,并将合成地震记录与过井地震剖面进行对比,确定地震剖面与所述测井曲线之间的时深关系,并精确标定各岩性界面在地震剖面上的反射位置;
d3、根据研究区地震解释层位、断层和主要目的层的沉积相分布特征,建立地质框架模型,根据框架内的地层接触关系是整合、上超、下超或是削截进一步完善地质框架模型,在此地质框架模型基础上,将测井的阻抗信息按照反距离乘方的方式在三维空间内进行内插和外推,建立三维的初始波阻抗约束模型;
d4、对原始地震数据进行波阻抗反演;
d5、在波阻抗反演结果的基础上进行所述最优参数的属性反演;
d6、对属性反演结果进行切片处理及地质解释;
d7、依据经步骤d6处理过的反演结果,最终确定合理的储层预测结果。
本发明的有益效果是:本发明的一种技术方案是利用地震资料和测井数据对交汇分析得到的每一个有效参数分别进行正演模拟,并对应进行特征响应分析,确定能够区分储层的最优参数,然后对最优参数进行参数反演,根据反演结果最终确定储层预测结果,该技术方案解决了当储层预测结果不满意时需要反复尝试不同的参数进行反演的问题,减少了储层预测的工作量,提高了工作效率。
进一步的,对所述正演模拟结果和反演结果进行对比分析后确定储层预测结果,不仅能够解决当储层预测结果不满意时需要反复尝试不同的参数进行反演的问题,还可以得到精度更高、范围更准确地储层预测结果,降低了勘探风险,是岩性油气藏储层预测的一种有效方法。
附图说明
图1为本发明实施例中预测的横波曲线示意图;
图2为本发明实施例中计算的弹性参数示意图;
图3为本发明实施例中测井曲线交汇分析示意图;
图4为本发明实施例中有效参数正演模拟结果示意图;
图5为本发明实施例中最优参数反演结果示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明的技术方案作进一步详细描述。
本发明实现其目的所采取的技术方案是通过对研究区内的测井数据进行交汇分析,确定能够区分储层的有效参数;综合利用地震资料和测井数据,对获得的每个有效参数分别进行正演模拟,并对应进行特征响应分析,确定能够区分储层的最优参数;对最优参数进行储层参数反演,根据反演结果确定合理的储层预测结果或将正演模拟结果和反演结果进行对比分析后确定储层预测结果。
本发明优选的一种实施例的具体步骤如下:
(1)收集研究区内的地震数据以及研究区内相应的测井数据。
收集研究区内的segy地震数据体以及测井数据;测井数据包括测井曲线、井分层及井斜数据,其中测井曲线包括声波曲线、电阻率曲线、自然电位曲线、孔隙度曲线、饱和度曲线、密度曲线、伽马曲线和泥质含量曲线。
(2)对测井数据进行交汇分析,确定能够区分储层的有效参数。
对测井数据进行交汇分析,确定测井数据的交汇图上的储层段是否具有明显特征;将能够表征储层段具有明显特征的参数确定为能够区分储层的有效参数,包含以下步骤:
a.通常情况下,测井曲线没有横波曲线,需要进行横波预测,使用修改的高斯曼方法结合经验公式法预测横波:
其中:ks是饱和岩石的体积模量,kd是干燥岩石格架的体积模量,km是岩石固体基质的体积模量,kf是孔隙流体的体积模量,是孔隙度,μs是饱和岩石的剪切模量,μd是干燥岩石格架的剪切模量,ρs是饱和岩石的密度,预测的横波曲线如图1所示。
b.利用纵波曲线,横波曲线,密度曲线,进行纵波波阻抗,横波波阻抗,拉梅系数,剪切模量,泊松比的计算,如图2所示,计算公式如下:
Ip=ρvp,Is=ρvs,λ=vp 2-2vs 2,μ=vs 2,
其中:Ip是纵波波阻抗,Is是横波波阻抗,λ是拉梅系数,μ是剪切模量,σ是泊松比。
c.对上述计算得到的参数进行交汇分析,确定有效参数,如图3所示。
(3)利用测井数据和地震数据,对每一个有效参数分别进行正演模拟,并对应进行特征响应分析,确定能够区分储层的最优参数,包含以下步骤:
a.提取地震子波或制作理论地震子波Wi(t);
b.根据地震波传播的褶积模型,利用测井曲线制作合成地震记录,并将合成地震记录与过井地震剖面进行对比,确定地震剖面与测井曲线之间的时深关系,其中公式如下:
Si(t)=Wi(t)*Ri(t),
其中Si(t)为合成地震记录,Wi(t)为地震子波,Ri(t)为反射系数序列,*为褶积符号;
c.导入地震解释层位和断层,依据地震解释层位和断层建立正演模型;
d.依据建立的正演模型,对每个有效参数分别进行正演模拟,并对应进行特征响应分析,当正演模拟剖面具有与井对比剖面相似储层特征,并且相对上下围岩,能够有效的区分储层时,选择对应的响应参数为区分储层的最优参数,如图4所示为最终确定的最优参数的正演结果。
(4)利用测井数据和地震数据,对经过正演模拟得到的最优参数进行储层参数反演,将得到的反演结果与所述正演模拟结果进行对比分析,从而确定合理的储层预测结果,包含以下步骤:
a.提取地震子波或制作理论地震子波;
b.根据地震波传播的褶积模型,利用测井曲线制作合成地震记录,并将合成记录与过井地震剖面进行对比,确定地震剖面与测井曲线之间的时深关系,并精确标定各岩性界面在地震剖面上的反射位置;
c.根据研究区地震解释层位、断层和主要目的层的沉积相分布特征,建立地质框架模型,根据框架内的地层接触关系是整合、上超、下超或是削截进一步完善地质框架模型,在此地质框架模型的基础上,将测井的阻抗信息按照反距离乘方的方式在三维空间内进行内插和外推,建立三维的初始波阻抗约束模型;
d.对原始地震数据进行波阻抗反演;
e.在波阻抗反演结果基础上进行最优参数的属性反演,得到最优参数反演结果,如图5所示;
f.对所述反演结果进行切片处理,并进行地质解释;
g.将经步骤f处理过的反演结果与所述正演模拟结果进行对比分析,最终确定合理的储层预测结果。例如,从图4的正演模拟剖面上可以看出,目的层表现为上倾尖灭的储层特征,图5的反演剖面上储层特征也表现为上倾尖灭,两者结果比较一致,那么将此预测结果确定为最终储层预测结果。这说明本发明实施例提供的一种储层预测方法提高了储层预测的精度,对油气有利位置的指示更为明显。
以上给出了一种优选的实施方式,但本发明不局限于所描述的实施方式。本发明的基本思路在于上述基本方案,对本领域普通技术人员而言,根据本发明的教导,设计出各种变形的模型、公式、参数并不需要花费创造性劳动。在不脱离本发明的原理和精神的情况下对实施方式进行的变化、修改、替换和变型仍落入本发明的保护范围内。
Claims (3)
1.一种储层预测方法,其特征在于,包含如下步骤:
1)收集研究区内的地震数据以及相应的测井数据;
2)对测井数据进行交汇分析,确定能够区分储层的有效参数;
3)利用测井数据和地震数据,对每一个有效参数分别进行正演模拟,并对应进行特征响应分析,确定能够区分储层的最优参数;
4)利用测井数据和地震数据,对所述最优参数进行储层参数反演,根据反演结果,最终确定储层预测结果;
步骤1)中所述的测井数据包括测井曲线,测井曲线包括声波曲线、电阻率曲线、自然电位曲线、孔隙度曲线、饱和度曲线、密度曲线、伽马曲线和泥质含量曲线;
步骤2)包含如下步骤:
b1、对所述测井曲线进行横波预测,得到横波曲线;
b2、利用纵波曲线、横波曲线、密度曲线,进行弹性参数计算;
b3、对所述弹性参数进行交汇分析,确定能够区分储层的有效参数;
步骤3)包含如下步骤:
c1、提取地震子波或制作理论地震子波;
c2、利用所述测井曲线制作合成地震记录,并将合成地震记录与过井地震剖面进行对比,确定地震剖面与所述测井曲线之间的时深关系;
c3、导入地震解释层位和断层,依据地震解释层位和断层建立正演模型;
c4、依据建立的正演模型,对每一个有效参数分别进行正演模拟,并对应进行特征响应分析,确定能够区分储层的最优参数。
2.根据权利要求1所述的一种储层预测方法,其特征在于,对正演模拟结果和反演结果进行对比分析,最终确定储层预测结果。
3.根据权利要求1所述的一种储层预测方法,其特征在于,步骤4)包含如下步骤:
d1、提取地震子波或制作理论地震子波;
d2、利用所述测井曲线制作合成地震记录,并将合成地震记录与过井地震剖面进行对比,确定地震剖面与所述测井曲线之间的时深关系,并精确标定各岩性界面在地震剖面上的反射位置;
d3、根据研究区地震解释层位、断层和主要目的层的沉积相分布特征,建立地质框架模型,根据框架内的地层接触关系是整合、上超、下超或是削截进一步完善地质框架模型,在此地质框架模型基础上,将测井的阻抗信息按照反距离乘方的方式在三维空间内进行内插和外推,建立三维的初始波阻抗约束模型;
d4、对原始地震数据进行波阻抗反演;
d5、在波阻抗反演结果的基础上进行所述最优参数的属性反演;
d6、对属性反演结果进行切片处理及地质解释;
d7、依据经步骤d6处理过的反演结果,最终确定合理的储层预测结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510490804.9A CN105652323B (zh) | 2015-08-11 | 2015-08-11 | 一种储层预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510490804.9A CN105652323B (zh) | 2015-08-11 | 2015-08-11 | 一种储层预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105652323A CN105652323A (zh) | 2016-06-08 |
CN105652323B true CN105652323B (zh) | 2017-12-29 |
Family
ID=56481689
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510490804.9A Active CN105652323B (zh) | 2015-08-11 | 2015-08-11 | 一种储层预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105652323B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105954804B (zh) * | 2016-07-15 | 2017-12-01 | 中国石油大学(北京) | 页岩气储层脆性地震预测方法及装置 |
CN106443779B (zh) * | 2016-08-17 | 2018-08-07 | 中国石油化工股份有限公司 | 一种中小型砂砾岩体的储层预测方法 |
CN106291713B (zh) * | 2016-08-24 | 2019-08-06 | 中国石油天然气股份有限公司 | 一种确定地层岩性的处理方法及装置 |
CN106569280B (zh) * | 2016-10-28 | 2021-08-20 | 中国石油化工股份有限公司 | 一种提高断阶构造地震解释精度的标定方法 |
CN107831537B (zh) * | 2017-01-05 | 2020-07-14 | 中国石油化工股份有限公司 | 高砂地比浊积岩储层描述方法 |
CN108375804A (zh) * | 2018-01-22 | 2018-08-07 | 中国石油大港油田勘探开发研究院 | 一种陆相碎屑岩“四相合一”沉积微相判别方法 |
CN108646294B (zh) * | 2018-06-29 | 2019-09-06 | 中海石油(中国)有限公司 | 一种复杂断裂系统下深水浊积砂岩储层高精度反演方法 |
CN110967742B (zh) * | 2018-09-28 | 2022-05-13 | 中国石油化工股份有限公司 | 一种孔隙度反演方法及系统 |
CN114063157B (zh) * | 2020-07-30 | 2024-05-17 | 中国石油化工股份有限公司 | 薄互层地震特征分析的方法、系统、存储介质及电子设备 |
CN112363243A (zh) * | 2020-10-26 | 2021-02-12 | 中国石油天然气集团有限公司 | 预测页岩油优质储层的方法、装置、计算机设备及介质 |
CN112363226A (zh) * | 2020-11-06 | 2021-02-12 | 大庆油田有限责任公司 | 一种非常规油气有利区地球物理预测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7657413B2 (en) * | 2002-07-11 | 2010-02-02 | Institut Francais Du Petrole | Method of constraining a heterogeneous permeability field representing an underground reservoir by dynamic data |
CN101980053A (zh) * | 2010-10-18 | 2011-02-23 | 中国石油化工股份有限公司 | 一种复杂礁滩储层预测方法 |
CN102253415A (zh) * | 2011-04-19 | 2011-11-23 | 中国石油大学(华东) | 基于裂缝等效介质模型的地震响应模式建立方法 |
CN104570128A (zh) * | 2013-10-16 | 2015-04-29 | 中国石油化工股份有限公司 | 一种基于岩石脆性因子的非常规油气储层预测方法 |
CN104749617A (zh) * | 2013-12-26 | 2015-07-01 | 中国石油化工股份有限公司 | 一种多尺度裂缝储层正演模型建立方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10422922B2 (en) * | 2012-05-24 | 2019-09-24 | Exxonmobil Upstream Research Company | Method for predicting rock strength by inverting petrophysical properties |
-
2015
- 2015-08-11 CN CN201510490804.9A patent/CN105652323B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7657413B2 (en) * | 2002-07-11 | 2010-02-02 | Institut Francais Du Petrole | Method of constraining a heterogeneous permeability field representing an underground reservoir by dynamic data |
CN101980053A (zh) * | 2010-10-18 | 2011-02-23 | 中国石油化工股份有限公司 | 一种复杂礁滩储层预测方法 |
CN102253415A (zh) * | 2011-04-19 | 2011-11-23 | 中国石油大学(华东) | 基于裂缝等效介质模型的地震响应模式建立方法 |
CN104570128A (zh) * | 2013-10-16 | 2015-04-29 | 中国石油化工股份有限公司 | 一种基于岩石脆性因子的非常规油气储层预测方法 |
CN104749617A (zh) * | 2013-12-26 | 2015-07-01 | 中国石油化工股份有限公司 | 一种多尺度裂缝储层正演模型建立方法 |
Non-Patent Citations (1)
Title |
---|
利用多元地震属性预测测井特性;石双虎 等;《西北地震学报》;20061231;第28卷(第4期);第309-313页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105652323A (zh) | 2016-06-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105652323B (zh) | 一种储层预测方法 | |
US10572611B2 (en) | Method and system for characterizing fractures in a subsurface region | |
Stephen et al. | Multiple-model seismic and production history matching: a case study | |
US10712472B2 (en) | Method and system for forming and using a subsurface model in hydrocarbon operations | |
US9026416B2 (en) | System and method of grid generation for discrete fracture modeling | |
CN111596978A (zh) | 用人工智能进行岩相分类的网页显示方法、模块和系统 | |
AU2015280679B2 (en) | Instantaneous isochron attribute-based geobody identification for reservoir modeling | |
Soleimani et al. | Integrated petrophysical modeling for a strongly heterogeneous and fractured reservoir, Sarvak Formation, SW Iran | |
CN110632657B (zh) | 泥岩涂抹型断层封闭性分析方法及装置 | |
CN108397130A (zh) | 钻井方法和装置 | |
CN105005079A (zh) | 一种井曲线反演方法 | |
Graves et al. | Observed and simulated ground motions in the San Bernardino basin region for the Hector Mine, California, earthquake | |
US20160377752A1 (en) | Method of Digitally Identifying Structural Traps | |
Zhang et al. | Architecture characteristics and characterization methods of fault-controlled karst reservoirs: A case study of the Shunbei 5 fault zone in the Tarim Basin, China | |
Hollis et al. | Uncertainty management in a giant fractured carbonate field, Oman, using experimental design | |
CN115880455A (zh) | 基于深度学习的三维智能插值方法 | |
Gómez et al. | Comparative analysis of the K-nearest-neighbour method and K-means cluster analysis for lithological interpretation of well logs of the Shushufindi Oilfield, Ecuador | |
CN112379435A (zh) | 相控岩溶型缝洞集合体刻画方法及装置 | |
Li et al. | Three-dimensional reservoir architecture modeling by geostatistical techniques in BD block, Jinhu depression, northern Jiangsu Basin, China | |
CN111815769A (zh) | 逆冲推覆构造带构造的建模方法、计算设备及存储介质 | |
Kondratyev et al. | Submarine fan reservoir architecture and heterogeneity influence on hard-to-recover reserves. Achimov Fm | |
Correia | Integration of reservoir characterization with history matching guided by pilot wells: Application to the Norne Field | |
Stephen et al. | A multiple model approach to history matching and uncertainty analysis using time-lapse seismic | |
Liu et al. | Application of Facies-controlled Modeling in Reservoir Prediction of Sublacustrine Fan | |
Singh et al. | Hybrid Strategy for Porosity Distribution Mapping on a Heterogeneous Reservoir Using Artificial Intelligence Method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |