CN105467466B - 一种基于多尺度信息的致密储层成岩相预测方法 - Google Patents
一种基于多尺度信息的致密储层成岩相预测方法 Download PDFInfo
- Publication number
- CN105467466B CN105467466B CN201510895321.7A CN201510895321A CN105467466B CN 105467466 B CN105467466 B CN 105467466B CN 201510895321 A CN201510895321 A CN 201510895321A CN 105467466 B CN105467466 B CN 105467466B
- Authority
- CN
- China
- Prior art keywords
- diagenesis
- quantitative
- sample point
- log
- parameter
- 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
- Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
Abstract
本发明涉及一种基于多尺度信息的致密储层成岩相预测方法,包括以下步骤:分析每一样本岩心薄片的成岩相,确定其成岩定量三参数(视压实率、视胶结率和视溶蚀率);读取每一薄片的成岩定量三参数和该深度测井值;每个样本的测井值和成岩定量三参数组成学习样本,利用BP神经网络预测全井段的成岩定量三参数;利用成岩定量三参数确定全井段测井成岩相;利用岩心数据分别建立地震弹性参数与成岩定量三参数的定量关系;利用建立的上述定量关系,用地震弹性参数反演结果分别预测成岩定量三参数平面分布;根据成岩定量参数平面分布确定成岩相平面分布。本发明实现了致密储层井间成岩相定量预测。
Description
技术领域
本发明涉及致密储层的地质技术领域,尤其涉及一种基于多尺度信息的致密储层成岩相预测方法。
背景技术
致密油气藏是近年来油气勘探开发的新领域,也是今后油气开发新的增长点,但是对于致密油气储层来说,如何寻找相对优质储层是最大的难题。对于裂缝不发育的致密储层来说,其物性主要受控于两方面因素:沉积和成岩作用。其中,沉积控制储层的原始物性,而成岩作用则控制了储层后期的物性变化,同时,这两方面因素也是相互联系的,沉积是成岩作用的基础。成岩相是沉积物在特定的沉积和成岩环境中,受到流体变化、构造运动等作用,经历复杂的成岩作用和成岩演化的产物,包含岩石组分、颗粒接触关系、胶结物、孔隙特征、裂缝等综合特征,能较宏观的反映不同成岩作用及成岩作用组合对储层物性的影响。目前对于油气储层的成岩作用研究和成岩相划分主要是根据能够反映岩心样品微观特征的扫描电镜、铸体薄片、阴极发光资料的分析来完成。然而,基于取样困难和节约成本的考虑,岩心薄片资料有限,利用薄片分析资料只能确定某个样本点的成岩相,而不能连续地反映整段地层的成岩相。
长期以来,对储层沉积作用和沉积特征的研究开展较多,但是对成岩作用的研究主要还是局限在取心段薄片分析的基础上,对于非取心段或者井间的成岩相,缺乏有效的研究方法。基于岩心样品分析的成岩相研究已经比较成熟,最直接的方法就是薄片观察。通过岩石薄片观察,确定样品中的主要成岩作用类型,进而确定成岩相类型。但是这种研究方法是依赖于岩石样品的薄片分析,对于非取心段来说,只有测井资料,没有岩心薄片,不能确定其成岩相,而对于井间地层更是如此。因此,目前成岩相研究的最大难点是井点的非取心段和井间地层。
前人研究中也有部分学者探索过利用测井资料解释成岩相的方法,这些方法主要有以下几类:一是建立成岩相与特定测井曲线之间的关系,利用单一测井曲线解释储层成岩相,但是实践效果均不理想,究其原因,单一测井曲线都是对储层某一种物理性质的反映,没有一种测井信息是与成岩相有直接的、良好的对应关系的,因此单一测井曲线解释成岩相的效果不好;二是利用测井相定性解释成岩相,这一方法的缺点是定量化程度不够,不便于不同储层之间成岩强度的比较;三是利用蜘蛛网图等分析方法,采用多测井参数解释成岩相,这类方法的缺点在于其是一种定性化的成岩相解释方法,判别结果不便于不同储层的比较,而且主观性较强,不便于大规模的工业化解释。对于井间成岩相的预测来说,难度就更大了。
发明内容
本发明的目的在于解决上述现有技术存在的缺陷,提供一种基于多尺度信息的致密储层成岩相预测方法。
一种基于多尺度信息的致密储层成岩相预测方法,包括如下步骤:
(a)使用偏光显微镜分析目标地层的岩心薄片,观察每片薄片中储层成岩作用类型,确定该薄片所对应样品点的储层成岩相类型;
(b)通过薄片分析,确定每一片薄片所对应样品点处储层的成岩定量参数,所述成岩定量参数包括视压实率、视胶结率和视溶蚀率;
(c)对岩心深度进行测井深度归位,读取深度归位后岩心薄片所在深度位置的测井曲线值和成岩定量参数,第一个样本点的岩心薄片对应的数据为
[dep1,log11,log21,…,logi1,Comp1,Res1,Cem1],
其中dep1为第一个样本点的深度;
log11为第一个样本点的第一条测井曲线的测井值;
log21为第一个样本点的第二条测井曲线的测井值;
logi1为第一个样本点的第i条测井曲线的测井值,i=1,2,…;
Comp1为第一个样本点薄片观察确定的视压实率;
Res1为第一个样本点薄片观察确定的视溶蚀率;
Cem1为第一个样本点薄片观察确定的视胶结率;
每个样本点薄片得到一组样本数据,依次编号,第j(j=1,2,3…)个样本点的薄片对应的数据为
[depj,log1j,log2j,…,logij,Compj,Resj,Cemj],
其中depj为第j个样本点的深度;
log1j为第j个样本点的第一条测井曲线的测井值;
log2j为第j个样本点的第二条测井曲线的测井值;
logij为第j个样本点的第i条测井曲线的测井值;
Compj为第j个样本点薄片观察确定的视压实率;
Resj为第j个样本点薄片观察确定的视溶蚀率;
Cemj为第j个样本点薄片观察确定的视胶结率;
(d)以步骤(c)中建立的岩心薄片样本点数据为学习样本,利用BP神经网络预测全井段的成岩定量参数,分别得到全井段的成岩定量参数曲线;
(e)按照不同成岩相类型的定量标准,利用测井解释的成岩定量参数曲线解释全井段成岩相;所述定量标准为:视压实率、视胶结率和视溶蚀率标准;
(f)在井点位置的地震弹性参数反演结果上,统计每个样本点的地震弹性参数值,利用井点数据分别建立地震弹性参数与每个成岩定量参数之间的定量关系;
(g)利用地震弹性参数反演结果和步骤(f)建立的定量关系,将地震弹性参数反演数据体转化为成岩定量参数数据,得到井间的成岩定量参数分布;
(h)根据不同成岩相的定量标准,利用井间的成岩定量参数预测结果,解释井间成岩相分布。
本发明综合岩心、测井和三维地震资料三类不同尺度的信息,联合预测不同尺度上的储层成岩相,实现成岩相空间分布的预测,解决了致密储层研究中的井间成岩相预测难题。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面对本发明中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的具体实施步骤如下:
(a)使用偏光显微镜分析目标地层的岩心薄片,观察每片薄片中储层成岩作用类型,确定该薄片所对应样品点的储层成岩相类型;
(b)通过薄片分析,确定每一片薄片所对应样品点处储层的成岩定量参数,包括视压实率、视胶结率和视溶蚀率;
(c)对岩心深度进行测井深度归位,读取深度归位后岩心薄片所在深度位置的测井曲线值和成岩定量参数,第一个样本点的岩心薄片对应的数据为
[dep1,log11,log21,…,logi1,Comp1,Res1,Cem1],
其中dep1为第一个样本点的深度;
log11为第一个样本点的第一条测井曲线的测井值;
log21为第一个样本点的第二条测井曲线的测井值;
logi1为第一个样本点的第i条测井曲线的测井值,i=1,2,…;
Comp1为第一个样本点薄片观察确定的视压实率;
Res1为第一个样本点薄片观察确定的视溶蚀率;
Cem1为第一个样本点薄片观察确定的视胶结率;
每个样本点薄片得到一组样本数据,依次编号,第j(j=1,2,3…)个样本点的薄片对应的数据为
[depj,log1j,log2j,…,logij,Compj,Resj,Cemj],
其中depj为第j个样本点的深度;
log1j为第j个样本点的第一条测井曲线的测井值;
log2j为第j个样本点的第二条测井曲线的测井值;
logij为第j个样本点的第i条测井曲线的测井值;
Compj为第j个样本点薄片观察确定的视压实率;
Resj为第j个样本点薄片观察确定的视溶蚀率;
Cemj为第j个样本点薄片观察确定的视胶结率;
(d)以步骤(c)中建立的岩心薄片样本点数据为学习样本,利用BP神经网络预测全井段的成岩定量参数,分别得到全井段的成岩定量参数曲线(视压实率曲线、视胶结率曲线和视溶蚀率曲线);
(e)按照不同成岩相类型的定量标准(视压实率、视胶结率和视溶蚀率标准),利用测井解释的成岩定量参数曲线解释全井段成岩相;
(f)在井点位置的地震弹性参数反演结果上,统计每个样本点的地震弹性参数值,利用井点数据分别建立地震弹性参数与每个成岩定量参数之间的定量关系;
(g)利用地震弹性参数反演结果和步骤(f)建立的定量关系,将地震弹性参数反演数据体转化为成岩定量参数数据,得到井间的成岩定量参数分布;
(h)根据不同成岩相的定量标准,利用井间的成岩定量参数预测结果,解释井间成岩相分布。
本发明综合岩心、测井和三维地震资料三类不同尺度的信息,联合预测不同尺度上的储层成岩相,实现成岩相空间分布的预测,解决了致密储层研究中的井间成岩相预测难题。
应用实例
我国东海盆地A气田面积140km2,仅有钻井3口,主要的天然气产层属于深层致密砂岩储层,裂缝不发育,其中的D5段气层中相对优质储层(“甜点”储层)的发育主要受成岩作用控制,有利的成岩相类型为中压实中溶蚀成岩相,需要开展成岩相分布预测,得到中压实中溶蚀成岩相的平面分布,从而预测“甜点”储层的分布,指导下一步气藏的钻探开发。在该实例中按照以下步骤开展工作:
(a)利用偏光显微镜分析目标地层的岩心薄片,对3口井D5段地层的210个岩石薄片样品进行观察分析,确定每片薄片中储层成岩作用类型,确定了研究区发育4种成岩相类型,包括中溶蚀-中压实成岩相、中溶蚀-中强压实成岩相、强溶蚀-强压实成岩相、中强溶蚀-强压实成岩相;
(b)通过薄片分析,确定每一片薄片所对应样品点处储层的成岩定量参数,包括视压实率、视胶结率和视溶蚀率,得到研究区210个样品点的成岩定量参数数据;
(c)对岩心深度进行测井深度归位,读取深度归位后岩心薄片所在深度位置的测井曲线值和成岩定量参数。本研究区选取的测井曲线包括声波时差、补偿中子、自然伽马、电阻率和密度测井;
(d)以步骤(c)中建立的岩心薄片样本点数据为学习样本,利用三层BP神经网络预测全井段的成岩定量三参数,分别得到全井段的成岩定量参数曲线(视压实率曲线、视胶结率曲线和视溶蚀率曲线);
(e)按照下表(表1)的不同成岩相定量标准(视压实率、视胶结率和视溶蚀率标准),利用测井解释的成岩定量参数曲线解释全井段成岩相;
表1成岩相参数的定量标准
(f)在研究区井点位置的地震弹性参数(纵横波速度比、纵波阻抗、横波阻抗)反演结果上,统计目的层段每个样本点的地震弹性参数值(纵横波速度比、纵波阻抗、横波阻抗),利用步骤(e)解释的井点成岩定量参数数据分别建立地震弹性参数值与每个成岩定量参数之间的定量关系,研究区的成岩相类型中只涉及到溶蚀作用和压实作用,所以本研究实例中建立的成岩定量参数与地震反演参数间的关系式,没有对视胶结率进行预测:
Comp=-1.68126e-6Ip+3.83944e-5Is-204.3903488Vp/Vs+103.14788
Res=4.07672e-6Ip+3.84935e-5Is-244.484517Vp/Vs+80.16936878
其中,Comp为视压实率,Res为视溶蚀率,Ip为纵波阻抗,Is为横波阻抗,Vp/Vs为纵横波速度比
(g)根据步骤(f)建立的定量关系,利用地震反演的纵横波速度比、横波阻抗、纵波阻抗三个反演数据体,转换得到视溶蚀率和视压实率参数数据体,得到井间的视溶蚀率和视压实率参数分布;
(h)按照表1中不同成岩相的定量标准,根据步骤(g)得到的井间成岩定量参数预测结果,解释井间成岩相分布。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (1)
1.一种基于多尺度信息的致密储层成岩相预测方法,其特征在于,包括如下步骤:
(a)使用偏光显微镜分析目标地层的岩心薄片,观察每片薄片中储层成岩作用类型,确定该薄片所对应样品点的储层成岩相类型;
(b)通过薄片分析,确定每一片薄片所对应样品点处储层的成岩定量参数,所述成岩定量参数包括视压实率、视胶结率和视溶蚀率;
(c)对岩心深度进行测井深度归位,读取深度归位后岩心薄片所在深度位置的测井曲线值和成岩定量参数,第一个样本点的岩心薄片对应的数据为
[dep1,log11,log21,…,logi1,Comp1,Res1,Cem1],
其中dep1为第一个样本点的深度;
log11为第一个样本点的第一条测井曲线的测井值;
log21为第一个样本点的第二条测井曲线的测井值;
logi1为第一个样本点的第i条测井曲线的测井值,i=1,2,…;
Comp1为第一个样本点薄片观察确定的视压实率;
Res1为第一个样本点薄片观察确定的视溶蚀率;
Cem1为第一个样本点薄片观察确定的视胶结率;
每个样本点薄片得到一组样本数据,依次编号,第j(j=1,2,3…)个样本点的薄片对应的数据为
[depj,log1j,log2j,…,logij,Compj,Resj,Cemj],
其中depj为第j个样本点的深度;
log1j为第j个样本点的第一条测井曲线的测井值;
log2j为第j个样本点的第二条测井曲线的测井值;
logij为第j个样本点的第i条测井曲线的测井值;
Compj为第j个样本点薄片观察确定的视压实率;
Resj为第j个样本点薄片观察确定的视溶蚀率;
Cemj为第j个样本点薄片观察确定的视胶结率;
(d)以步骤(c)中建立的岩心薄片样本点数据为学习样本,利用BP神经网络预测全井段的成岩定量参数,分别得到全井段的成岩定量参数曲线;
(e)按照不同成岩相类型的定量标准,利用测井解释的成岩定量参数曲线解释全井段成岩相;所述定量标准为:视压实率、视胶结率和视溶蚀率标准;
(f)在井点位置的地震弹性参数反演结果上,统计每个样本点的地震弹性参数值,利用井点数据分别建立地震弹性参数与每个成岩定量参数之间的定量关系;
(g)利用地震弹性参数反演结果和步骤(f)建立的定量关系,将地震弹性参数反演数据体转化为成岩定量参数数据,得到井间的成岩定量参数分布;
(h)根据不同成岩相的定量标准,利用井间的成岩定量参数预测结果,解释井间成岩相分布。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510895321.7A CN105467466B (zh) | 2015-12-07 | 2015-12-07 | 一种基于多尺度信息的致密储层成岩相预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510895321.7A CN105467466B (zh) | 2015-12-07 | 2015-12-07 | 一种基于多尺度信息的致密储层成岩相预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105467466A CN105467466A (zh) | 2016-04-06 |
CN105467466B true CN105467466B (zh) | 2018-01-09 |
Family
ID=55605362
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510895321.7A Active CN105467466B (zh) | 2015-12-07 | 2015-12-07 | 一种基于多尺度信息的致密储层成岩相预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105467466B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108090278B (zh) * | 2017-12-15 | 2021-06-01 | 长江大学 | 碎屑岩储层成岩相的划分方法 |
CN108375804A (zh) * | 2018-01-22 | 2018-08-07 | 中国石油大港油田勘探开发研究院 | 一种陆相碎屑岩“四相合一”沉积微相判别方法 |
CN108345046B (zh) * | 2018-02-06 | 2020-07-03 | 长江大学 | 碎屑岩储层埋藏过程中成岩相划分方法 |
CN114113537B (zh) * | 2020-09-01 | 2024-03-26 | 中国石油天然气股份有限公司 | 定量识别致密砂岩成岩作用的方法和装置 |
CN115203970B (zh) * | 2022-08-03 | 2023-04-07 | 中国石油大学(北京) | 基于人工智能算法的成岩参数预测模型训练方法和预测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103913556A (zh) * | 2014-03-10 | 2014-07-09 | 中国石油天然气股份有限公司吉林油田分公司勘探开发研究院 | 确定致密油气储层优势成岩相的岩心分析技术 |
CN104360039A (zh) * | 2014-10-31 | 2015-02-18 | 中国石油化工股份有限公司 | 一种致密砂岩储层成岩相定量评价方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB0017227D0 (en) * | 2000-07-14 | 2000-08-30 | Schlumberger Ind Ltd | Fully coupled geomechanics in a commerical reservoir simulator |
-
2015
- 2015-12-07 CN CN201510895321.7A patent/CN105467466B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103913556A (zh) * | 2014-03-10 | 2014-07-09 | 中国石油天然气股份有限公司吉林油田分公司勘探开发研究院 | 确定致密油气储层优势成岩相的岩心分析技术 |
CN104360039A (zh) * | 2014-10-31 | 2015-02-18 | 中国石油化工股份有限公司 | 一种致密砂岩储层成岩相定量评价方法 |
Non-Patent Citations (4)
Title |
---|
"Micro-fracture characteristics of tight sandstone reservoirs and its";Peng Zhu 等;《Journal of Natural Gas Science and Engineering》;20150822;第71卷(第3期);第90-97页 * |
"川中低孔渗砂岩成岩相定量评价与快速预测--以遂宁—蓬溪须二段为例";蒋裕强 等;《中国地质》;20140430;第41卷(第2期);第437-449页 * |
"成岩相的形成_分类与定量评价方法";邹才能 等;《石油勘探与开发》;20081031;第35卷(第5期);第526-540页 * |
"致密砂岩储集层成岩相定量划分及其测井识别方法";赖锦 等;《矿物岩石地球化学通报》;20150131;第34卷(第1期);第128-138页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105467466A (zh) | 2016-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sun et al. | Shale gas exploration and development in China: Current status, geological challenges, and future directions | |
Xu et al. | Multiscale faults and fractures characterization and their effects on shale gas accumulation in the Jiaoshiba area, Sichuan Basin, China | |
CN105467466B (zh) | 一种基于多尺度信息的致密储层成岩相预测方法 | |
Zhu et al. | Key factors of marine shale conductivity in southern China—Part II: The influence of pore system and the development direction of shale gas saturation models | |
CN103869052B (zh) | 一种砂岩储层成岩测井相定量表征的方法 | |
CN103852787B (zh) | 一种砂岩储层成岩地震相表征方法 | |
CN1275048C (zh) | 用于浅流探测的方法 | |
CN103529474A (zh) | 采用岩性细分实现岩相精细描述的方法 | |
Wang et al. | Development characteristics and dominant factors of fractures and their significance for shale reservoirs: a case study from Є1b2 in the Cen'gong block, southern China | |
CN112577875A (zh) | 一种碳酸盐岩多种原生孔隙高效探测技术 | |
Wang et al. | Multi-scale natural fracture prediction in continental shale oil reservoirs: a case study of the Fengcheng Formation in the Mahu Sag, Junggar Basin, China | |
Suyun et al. | Migration and accumulation mechanisms and main controlling factors of tight oil enrichment in a continental lake basin | |
Yu et al. | Evaluation of low permeability conglomerate reservoirs based on petrophysical facies: A case study from the Triassic Baikouquan Formation, northern Mahu Sag, Junggar Basin, China | |
Hussain et al. | Non-destructive techniques for linking methodology of geochemical and mechanical properties of rock samples | |
Zhang et al. | Effects of diagenesis on natural fractures in tight oil reservoirs: a case study of the Permian Lucaogou Formation in Jimusar Sag, Junggar Basin, NW China | |
Galford et al. | Quantifying organic porosity from logs | |
Han et al. | Formation and accumulation of lower Jurassic tight gas sands field in Kekeya area of Tuha Basin, northwestern China | |
Fan et al. | Semi-supervised learning–based petrophysical facies division and “Sweet Spot” identification of low-permeability sandstone reservoir | |
Zhang et al. | Multi-parameters logging identifying method for sand body architectures of tight sandstones: A case from the Triassic Chang 9 Member, Longdong area, Ordos Basin, NW China | |
Deng et al. | Factors controlling reservoir quality of a retreating delta-front in shallow-water lacustrine in the Songliao Basin, Northeast China | |
Li et al. | Study on Anisotropic Petrophysical Modeling of Shale: A Case Study of Shale Oil in Qingshankou Formation in Sanzhao Sag, Songliao Basin, China | |
CN110133722B (zh) | 一种应用声速定性识别泥质分布类型的方法 | |
CN110568149B (zh) | 沉积盆地烃源岩生排烃史精细快速定量模拟方法 | |
Zhang et al. | Improved stochastic simulation using stratigraphic forward modeling: a case study of the lithological distribution of tide-dominated estuary in JE-AW Oil Field, Ecuador | |
Grier et al. | Reservoir quality: Part 6. Geological methods |
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 |