CN109752302B - 一种预测致密砂岩储层上倾方向孔隙度的方法 - Google Patents

一种预测致密砂岩储层上倾方向孔隙度的方法 Download PDF

Info

Publication number
CN109752302B
CN109752302B CN201711061971.7A CN201711061971A CN109752302B CN 109752302 B CN109752302 B CN 109752302B CN 201711061971 A CN201711061971 A CN 201711061971A CN 109752302 B CN109752302 B CN 109752302B
Authority
CN
China
Prior art keywords
sandstone reservoir
porosity
upward
value
inclination direction
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
Application number
CN201711061971.7A
Other languages
English (en)
Other versions
CN109752302A (zh
Inventor
崔宇
付立新
吴雪松
李宏军
王辉
段润梅
冯建园
赵勇刚
楼达
张津宁
刁丽颖
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201711061971.7A priority Critical patent/CN109752302B/zh
Publication of CN109752302A publication Critical patent/CN109752302A/zh
Application granted granted Critical
Publication of CN109752302B publication Critical patent/CN109752302B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明公开了一种预测致密砂岩储层上倾方向孔隙度的方法,属于油田地质勘探技术领域。所述预测致密砂岩储层上倾方向孔隙度的方法通过获取致密砂岩储层的探井区域的孔隙度取值、石英次生加大边的体积百分比、高岭石的体积百分比、正长石的体积百分比和斜长石的体积百分比,计算得到溶蚀系数,再根据孔隙度取值、石英次生加大边的体积百分比、高岭石的体积百分比和溶蚀系数,计算得到致密砂岩储层上倾方向孔隙度的预测取值,实现根据已有的致密砂岩储层的物性资料推测量化得到其上倾方向孔隙度,方法简单易行,经济成本低,预测结果可靠性高。

Description

一种预测致密砂岩储层上倾方向孔隙度的方法
技术领域
本发明涉及石油地质勘探技术领域,特别涉及一种预测致密砂岩储层上倾方向孔隙度的方法。
背景技术
在对致密砂岩储层进行分析中,需要搜集钻井资料中致密砂岩储层的孔隙度、渗透率、含油饱和度、含水饱和度等储层物性资料,用以预测油藏的储量,并为钻井决策提供依据。
现有技术中致密砂岩储层上倾方向的孔隙度是通过钻井取芯,对取出的岩心进行分析得到的,而不能根据已测得的致密砂岩储层的物性资料预测得到其上倾方向的孔隙度。
在实现本发明的过程中,本发明人发现现有技术中至少存在以下问题:
现有技术在获取致密砂岩储层上倾方向的孔隙度时,需要钻井队对致密砂岩储层上倾方向的岩层进行钻探,获取岩心,并对岩心进行分析,过程复杂,经济成本高,施工时间长,工作效率低。
发明内容
鉴于此,本发明提供一种预测致密砂岩储层上倾方向孔隙度的方法,可以根据已有的致密砂岩储层的物性资料推测量化得到其上倾方向孔隙度,且简单易行。
具体而言,包括以下的技术方案:
一种预测致密砂岩储层上倾方向孔隙度的方法,所述方法包括:
获取致密砂岩储层的探井区域的孔隙度取值;
获取致密砂岩储层的探井区域的石英次生加大边的体积百分比和高岭石的体积百分比;
获取致密砂岩储层的探井区域中正长石的体积百分比和斜长石的体积百分比;
根据所述正长石的体积百分比和斜长石的体积百分比,计算得到溶蚀系数;
根据致密砂岩储层的探井区域的所述孔隙度取值、所述石英次生加大边的体积百分比、所述高岭石的体积百分比和所述溶蚀系数,计算得到致密砂岩储层上倾方向孔隙度的预测取值,且所述致密砂岩储层上倾方向孔隙度的预测取值的计算公式为:
φ=φ1+V1+V2+a(V1+V2)
式中:φ为致密砂岩储层上倾方向孔隙度的预测取值;φ1为致密砂岩储层的探井区域的孔隙度取值;V1为致密砂岩储层的探井区域的石英次生加大边的体积百分比;V2为致密砂岩储层的探井区域的高岭石的体积百分比;a为溶蚀系数。
进一步地,所述溶蚀系数的计算公式为:
a=0.88v1+0.95v2
式中:v1为所述正长石的体积百分比;v2为所述斜长石的体积百分比。
进一步地,所述探井区域的石英次生加大边的体积百分比和高岭石的体积百分比是利用显微镜对探井得到的致密砂岩储层的岩石薄片进行观察得到的。
进一步地,所述计算得到致密砂岩储层上倾方向孔隙度的预测取值之后,所述方法还包括:对致密砂岩储层上倾方向进行钻井,得到所述致密砂岩储层上倾方向的岩心。
进一步地,所述对致密砂岩储层上倾方向进行钻井,得到所述致密砂岩储层上倾方向的岩心之后,所述方法还包括:根据所述岩心,得到致密砂岩储层上倾方向孔隙度的实际取值。
进一步地,所述根据所述岩心,得到致密砂岩储层上倾方向孔隙度的实际取值之后,所述方法还包括:根据所述致密砂岩储层上倾方向孔隙度的实际取值和所述致密砂岩储层上倾方向孔隙度的预测取值,计算得到取值误差。
进一步地,若所述取值误差小于等于预设的误差阈值,则确定所述致密砂岩储层上倾方向孔隙度的预测取值有效;若所述取值误差大于所述预设的误差阈值,则确定所述致密砂岩储层上倾方向孔隙度的预测取值无效。
进一步地,所述取值误差的计算公式为:
Figure BDA0001454964810000021
式中:δ为取值误差;φ'为致密砂岩储层上倾方向孔隙度的实际取值。
本发明实施例提供的技术方案的有益效果:
本发明实施例提供的一种预测致密砂岩储层上倾方向孔隙度的方法,通过获取致密砂岩储层的探井区域的孔隙度取值、石英次生加大边的体积百分比、高岭石的体积百分比、正长石的体积百分比和斜长石的体积百分比,计算得到溶蚀系数,再根据孔隙度取值、石英次生加大边的体积百分比、高岭石的体积百分比和溶蚀系数,计算得到致密砂岩储层上倾方向孔隙度的预测取值,实现根据已有的致密砂岩储层的物性资料推测量化得到其上倾方向孔隙度,方法简单易行,经济成本低,预测结果可靠性高。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种预测致密砂岩储层上倾方向孔隙度的方法的方法结构示意图。
具体实施方式
为使本发明的技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
本发明实施例提供了一种预测致密砂岩储层上倾方向孔隙度的方法,以大港油田扣村油田庄63井为例,其方法结构示意图如图1所示,包括如下的步骤:
步骤101:获取致密砂岩储层的探井区域的孔隙度取值。
具体地,施工人员选定探井区域,通过钻井施工队对致密砂岩储层进行钻井,获取致密砂岩储层处的岩心;对岩心进行分析,得到探井区域的物性参数资料,包括孔隙度、渗透率、含水饱和度和含油饱和度等;从探井区域的物性参数资料中得到致密砂岩储层的探井区域的孔隙度取值。
在本发明实施例的预测致密砂岩储层上倾方向孔隙度的方法中,庄63井所在的区域的孔隙度取值为6.3%。
步骤102:获取致密砂岩储层的探井区域的石英次生加大边的体积百分比和高岭石的体积百分比。
需要说明的是,致密砂岩储层岩性致密,其中石英颗粒及长石颗粒占到总体积的95%以上,长石颗粒受到大气淡水林滤作用的溶蚀,形成副产物次生高岭石及石英,两者顺水流而下,沿储层下倾方向堆积,在下倾方向的储层中以石英次生加大边及自生高岭石的状态赋存于孔隙中,造成下倾方向储层物性的变差,同时,由于参加反应的长石及生成物高岭石和石英摩尔体积不同,反应过程中有一定的体积损耗,使得致密砂岩储层中不同高度储层的孔隙度取值不同,且孔隙度的取值与石英次生加大边的体积和高岭石的体积相关。因此,在已知探井区域的致密砂岩储层的孔隙度取值的情况下,求取探井区域致密砂岩储层上倾方向的孔隙度取值的大小时,就需要利用到石英次生加大边的体积百分比与高岭石的体积百分比关系。
在本步骤中,探井区域的石英次生加大边的体积百分比和高岭石的体积百分比是利用显微镜对探井得到的致密砂岩储层的岩石薄片进行观察,并记录得到的。
在本发明实施例的预测致密砂岩储层上倾方向孔隙度的方法中,通过显微镜对从庄63井得到的致密砂岩储层的岩石薄片观察记录,确定次生石英加大边的体积百分比为1%,高岭石的体积百分比为3%。
步骤103:获取致密砂岩储层的探井区域中正长石的体积百分比和斜长石的体积百分比。
具体地,可以通过显微镜对从探井区域得到的致密砂岩储层的岩石薄片观察记录,得到致密砂岩储层的探井区域中正长石的体积百分比和斜长石的体积百分比。
步骤104:根据正长石的体积百分比和斜长石的体积百分比,计算得到溶蚀系数。
具体地,溶蚀系数的计算公式为:
a=0.88v1+0.95v2
式中:v1为正长石的体积百分比;v2为斜长石的体积百分比。
需要说明的是,根据正长石溶蚀体积变化量表,得到系数取值0.88,如下表表1所示;根据斜长石溶蚀体积变化量表,得到系数取值0.95,如下表表2所示。
表1正长石溶蚀体积变化量表
Figure BDA0001454964810000051
注:产物、反应物体积之比为0.88(长石溶蚀之后总体积损耗12%)
表2斜长石溶蚀体积变化量表
Figure BDA0001454964810000052
注:产物、反应物体积之比为0.95(长石溶蚀之后总体积损耗5%)
在本发明实施例的预测致密砂储层上倾方向孔隙度的方法中,通过显微镜对庄63井得到的致密砂岩储层的岩石薄片观察记录,得到正长石的体积百分比为50%,斜长石的体积百分比为50%,那么,根据溶蚀系数的计算公式,得到溶蚀系数a为0.92。
步骤105:根据致密砂岩储层的探井区域的孔隙度取值、石英次生加大边的体积百分比、高岭石的体积百分比和溶蚀系数,计算得到致密砂岩储层上倾方向孔隙度的预测取值。
具体地,致密砂岩储层上倾方向孔隙度的预测取值的计算公式为:
φ=φ1+V1+V2+a(V1+V2)
式中:φ为致密砂岩储层上倾方向孔隙度的预测取值;φ1为致密砂岩储层的探井区域的孔隙度取值;V1为致密砂岩储层的探井区域的石英次生加大边的体积百分比;V2为致密砂岩储层的探井区域的高岭石的体积百分比;a为溶蚀系数。
需要说明的是,在致密砂岩储层上倾方向孔隙度的预测取值计算中,由于致密砂岩储层已经固结成岩石,因此,在计算孔隙度的预测取值时不考虑孔隙度与垂向高度的关系,并将计算模型理想化,限定在本发明的预测致密砂岩储层上倾方向孔隙度的方法中孔隙度的取值只和溶蚀作用有关。
根据致密砂岩储层上倾方向孔隙度的预测取值的计算公式,带入庄63井的孔隙度取值、石英次生加大边的体积百分比和高岭石的体积百分比,得到致密砂岩储层上倾方向孔隙度的预测取值,计算得到φ的取值为13.98%。
进一步地,为了验证预测结果的准确性,在计算得到致密砂岩储层上倾方向孔隙度的预测取值之后,对致密砂岩储层上倾方向进行钻井,得到致密砂岩储层上倾方向的岩心;根据岩心,得到致密砂岩储层上倾方向孔隙度的实际取值;根据致密砂岩储层上倾方向孔隙度的实际取值和致密砂岩储层上倾方向孔隙度的预测取值,计算得到取值误差,取值误差的计算公式为:
Figure BDA0001454964810000061
式中:δ为取值误差;φ'为致密砂岩储层上倾方向孔隙度的实际取值。
若取值误差小于等于预设的误差阈值,则确定致密砂岩储层上倾方向孔隙度的预测取值有效;若取值误差大于预设的误差阈值,则确定致密砂岩储层上倾方向孔隙度的预测取值无效。需要说明的是,在本发明实施例的预测致密砂岩储层上倾方向孔隙度的方法中,预设的误差阈值为10%。
在本发明实施例的预测致密砂岩储层上倾方向孔隙度的方法中,通过对庄63井上倾方向进行钻取扣24井,得到致密砂岩储层上倾方向的岩心,根据岩心,得到致密砂岩储层上倾方向孔隙度的实际取值,即为15.1%。将孔隙度的预测取值和孔隙度的实际取值带入到取值误差公式中,计算得到取值误差为7.4%。
因此,计算得到的取值误差小于预设的误差阈值,那么确定致密砂层储层上倾方向孔隙度的预测取值有效,证明预测致密砂岩储层上倾方向孔隙度的方法有效。
本发明实施例通过获取致密砂岩储层的探井区域的孔隙度取值、石英次生加大边的体积百分比、高岭石的体积百分比、正长石的体积百分比和斜长石的体积百分比,计算得到溶蚀系数,再根据孔隙度取值、石英次生加大边的体积百分比、高岭石的体积百分比和溶蚀系数,计算得到致密砂岩储层上倾方向孔隙度的预测取值,实现根据已有的致密砂岩储层的物性资料推测量化得到其上倾方向孔隙度,方法简单易行,经济成本低,预测结果可靠性高。
以上所述仅是为了便于本领域的技术人员理解本发明的技术方案,并不用以限制本发明。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (6)

1.一种预测致密砂岩储层上倾方向孔隙度的方法,其特征在于,所述方法包括:
获取致密砂岩储层的探井区域的孔隙度取值;
获取致密砂岩储层的探井区域的石英次生加大边的体积百分比和高岭石的体积百分比;
获取致密砂岩储层的探井区域中正长石的体积百分比和斜长石的体积百分比;
根据所述正长石的体积百分比和斜长石的体积百分比,计算得到溶蚀系数;
根据致密砂岩储层的探井区域的所述孔隙度取值、所述石英次生加大边的体积百分比、所述高岭石的体积百分比和所述溶蚀系数,计算得到致密砂岩储层上倾方向孔隙度的预测取值,且所述致密砂岩储层上倾方向孔隙度的预测取值的计算公式为:
φ=φ1+V1+V2+a(V1+V2)
式中:φ为致密砂岩储层上倾方向孔隙度的预测取值;φ1为致密砂岩储层的探井区域的孔隙度取值;V1为致密砂岩储层的探井区域的石英次生加大边的体积百分比;V2为致密砂岩储层的探井区域的高岭石的体积百分比;a为溶蚀系数。
2.根据权利要求1所述的预测致密砂岩储层上倾方向孔隙度的方法,其特征在于,所述溶蚀系数的计算公式为:
a=0.88v11+0.95v22
式中:v11为所述正长石的体积百分比;v22为所述斜长石的体积百分比。
3.根据权利要求1所述的预测致密砂岩储层上倾方向孔隙度的方法,其特征在于,所述探井区域的石英次生加大边的体积百分比和高岭石的体积百分比是利用显微镜对探井得到的致密砂岩储层的岩石薄片进行观察得到的。
4.根据权利要求1所述的预测致密砂岩储层上倾方向孔隙度的方法,其特征在于,所述计算得到致密砂岩储层上倾方向孔隙度的预测取值之后,所述方法还包括:对致密砂岩储层上倾方向进行钻井,得到所述致密砂岩储层上倾方向的岩心;根据所述岩心,得到致密砂岩储层上倾方向孔隙度的实际取值;根据所述致密砂岩储层上倾方向孔隙度的实际取值和所述致密砂岩储层上倾方向孔隙度的预测取值,计算得到取值误差。
5.根据权利要求4所述的预测致密砂岩储层上倾方向孔隙度的方法,其特征在于,若所述取值误差小于等于预设的误差阈值,则确定所述致密砂岩储层上倾方向孔隙度的预测取值有效;若所述取值误差大于所述预设的误差阈值,则确定所述致密砂岩储层上倾方向孔隙度的预测取值无效。
6.根据权利要求4所述的预测致密砂岩储层上倾方向孔隙度的方法,其特征在于,所述取值误差的计算公式为:
Figure FDA0003132694490000021
式中:δ为取值误差;φ'为致密砂岩储层上倾方向孔隙度的实际取值。
CN201711061971.7A 2017-11-02 2017-11-02 一种预测致密砂岩储层上倾方向孔隙度的方法 Active CN109752302B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711061971.7A CN109752302B (zh) 2017-11-02 2017-11-02 一种预测致密砂岩储层上倾方向孔隙度的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711061971.7A CN109752302B (zh) 2017-11-02 2017-11-02 一种预测致密砂岩储层上倾方向孔隙度的方法

Publications (2)

Publication Number Publication Date
CN109752302A CN109752302A (zh) 2019-05-14
CN109752302B true CN109752302B (zh) 2021-09-28

Family

ID=66397773

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711061971.7A Active CN109752302B (zh) 2017-11-02 2017-11-02 一种预测致密砂岩储层上倾方向孔隙度的方法

Country Status (1)

Country Link
CN (1) CN109752302B (zh)

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2716536B1 (fr) * 1994-02-22 1996-04-26 Geophysique Cie Gle Procédé et dispositif pour mesurer la perméabilité d'un milieu rocheux .
US6860147B2 (en) * 2002-09-30 2005-03-01 Alberta Research Council Inc. Process for predicting porosity and permeability of a coal bed
CN104655540B (zh) * 2013-11-19 2017-09-12 中国石油大学(北京) 一种确定成藏期致密砂岩古孔隙度方法
CN103760081B (zh) * 2013-12-31 2016-01-06 中国石油天然气股份有限公司 基于孔隙结构特征的碳酸盐岩储层的气藏预测方法及系统
EP3074913B1 (en) * 2014-01-02 2023-09-27 Landmark Graphics Corporation History matching multi-porosity solutions
CN104360412B (zh) * 2014-11-14 2017-12-26 中国石油大学(北京) 致密深盆气成藏预测方法和装置
CN105510207B (zh) * 2016-01-29 2018-02-06 河海大学 确定裂隙岩体不同方向渗透率对孔隙度敏感程度的方法
CN106290105B (zh) * 2016-07-20 2017-08-11 中国石油大学(华东) 一种碳酸盐岩储层溶蚀孔隙体积含量预测方法
CN106597548B (zh) * 2016-12-02 2018-12-28 中国石油大学(华东) 一种多因素定量评价地质时期3d孔隙度的方法
CN106841001B (zh) * 2017-01-17 2020-03-10 重庆科技学院 一种基于储层质量主控因素分析的致密砂岩孔隙度、渗透率预测方法

Also Published As

Publication number Publication date
CN109752302A (zh) 2019-05-14

Similar Documents

Publication Publication Date Title
CN104806232B (zh) 一种确定碳酸盐岩储层孔隙度下限的方法
US20240151139A1 (en) Systems and Methods for Identifying Two or More Charges into Reservoir Using Downhole Fluid Analysis
CN106223941B (zh) 基于测井资料的页岩气储层有机孔隙度确定方法
CN104833789B (zh) 利用致密砂岩微观孔隙结构确定气水关系的方法
Male et al. Forecasting production from Bakken and three forks wells using a segregated flow model
CN104428492A (zh) 使用矿物学确定压裂刺激点的布置的系统、方法、和计算机程序产品
CN107725044B (zh) 基于阵列感应、侧向测井的砂岩含气储层产水率预测的方法
US9535187B2 (en) Device and method to determine conductivity for high pressure-high temperature service
CN111963159A (zh) 一种砾岩致密油储层流体性质的识别方法
CN112145165A (zh) 一种微裂缝-孔隙型储层动静态渗透率转换方法
CN111764895B (zh) 一种适用于页岩气储层地质模型的测井评价方法
CN109752302B (zh) 一种预测致密砂岩储层上倾方向孔隙度的方法
CN110244358B (zh) 一种构造破坏导致的油气逸散区的判识方法
CN110244356B (zh) 一种构造破坏形成的后期油气藏的判识方法
CN112282743A (zh) 用于对钻井泥岩地层压力进行预测的方法
Chunhua et al. Depression and apophysis structure type and its controls on deposition in the Jiyang Depression, Bohai Bay Basin
La Ode Ahdyar et al. Integrated carbonate non-matrix characterization in Banyu Urip field
CN104463362A (zh) 一种预测储层流体产量的方法及装置
CN115099014A (zh) 一种基于随钻录井的天然气井地质探明储量计算方法
CN114060015B (zh) 一种致密砂岩含气性的评价方法及装置
RU2515629C1 (ru) Способ определения хрупких зон коллекторов
US20200399987A1 (en) System and Method for Predicting Fluid Behavior in an Unconventional Shale Play
KR101818098B1 (ko) 물리검층 자료를 이용한 지층 내 점토 함량 산출방법
CN110930020A (zh) 非常规油气资源的经济可采资源量确定方法
CN110346416A (zh) 基于声波和电阻率的特征参数曲线识别火山岩岩性的方法

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