CN102979517A - 一种复杂油气储层饱和度定量评价方法 - Google Patents
一种复杂油气储层饱和度定量评价方法 Download PDFInfo
- Publication number
- CN102979517A CN102979517A CN2012105141577A CN201210514157A CN102979517A CN 102979517 A CN102979517 A CN 102979517A CN 2012105141577 A CN2012105141577 A CN 2012105141577A CN 201210514157 A CN201210514157 A CN 201210514157A CN 102979517 A CN102979517 A CN 102979517A
- Authority
- CN
- China
- Prior art keywords
- porosity
- formula
- value
- quantitatively evaluating
- gas reservoir
- 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.)
- Granted
Links
Images
Landscapes
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Abstract
本发明涉及一种复杂油气储层饱和度定量评价方法,它基于麦克斯韦导电方程,推导出了岩石导电孔隙度的计算公式,并结合十二种岩性的岩电实验结果,推导出胶结指数与其相应的有效孔隙度和导电孔隙度差值之间关系,从而准确地建立了阿尔奇公式中胶结指数的高精度计算模型。本发明还建立了有效孔隙度与导电孔隙度之间关系,从而胶结指数的计算模型能够以导电孔隙度为桥梁,推导出阿尔奇公式中胶结指数与有效孔隙度之间的普适性表达式,因此本发明具有较强的通用性,能够更全面地反映复杂储层的导电机理。本发明为利用阿尔奇公式定量评价复杂孔隙型储层的饱和度,提供了一种精度高、通用性强的有效方法,它可以广泛用于各种复杂孔隙型油气储层饱和度的定量评价中。
Description
技术领域
本发明涉及一种定量评价方法,特别是关于一种复杂油气储层饱和度定量评价方法。
背景技术
常用的油气储层饱和度定量评价方法是根据阿尔奇(Archie)公式建立起来的。阿尔奇公式中包括两个公式,其中,
地层因素F的计算公式为:
电阻增大系数I的计算公式为:
式中:Rt为含油气纯岩石电阻率,单位为Ω·m;b为与岩性有关的比例系数,无量纲;Sw为含水饱和度,无量纲;n为饱和指数,无量纲。
联立公式(1)和(2),便可以得到含水饱和度Sw的公式:
根据上述公式,当阿尔奇参数a、b、m和n取固定值时(例如a=b=1,m=n=2时),常规砂岩储层的含油气饱和度是可以被准确计算出来的。但是,对于孔隙结构十分复杂的火山岩和碳酸盐岩储层,“一成不变”的阿尔奇参数所计算的油气饱和度难以满足当前勘探开发的生产需求。由于阿尔奇参数受诸多因素影响,因此在使用时不应为一固定不变的定值,而应随孔隙结构的变化而变化。准确地建立可变阿尔奇参数的计算方法,是解决复杂油气储层饱和度定量评价的方法之一。
发明内容
针对上述问题,本发明的目的是提供一种建立高精度胶结指数m值模型来确定复杂油气储层饱和度定量评价方法。
计算出被评价地区复杂油气储层的含水饱和度Sw。
③当有效孔隙度全部导电时,根据式(5)能够推导出:
④联立式(6)和式(7),求解出系数c1和c2:
⑥根据雨果·弗瑞克修改的麦克斯韦导电方程的通用形式:
所述步骤3)中,公式(2)通过雨果·弗瑞克修改的麦克斯韦导电方程的通用形式,
本发明由于采取以上技术方案,其具有以下优点:1、本发明基于麦克斯韦导电方程理论上推导出了岩石导电孔隙度的计算公式,并结合十二种岩性的岩电实验结果,推导出胶结指数与其相应的有效孔隙度和导电孔隙度差值之间关系,从而准确地建立了阿尔奇公式中胶结指数的计算模型。2、本发明由于还建立了有效孔隙度与导电孔隙度之间关系,从而胶结指数的计算模型能够以导电孔隙度为桥梁,推导出阿尔奇公式中胶结指数与有效孔隙度之间的普适性表达式,因此本发明具有较强的通用性,能够更全面地反映复杂储层的导电机理。本发明为利用阿尔奇公式定量评价复杂孔隙型储层的饱和度,提供了一种精度高、通用性强的有效方法,它可以广泛用于各种复杂孔隙型油气储层饱和度的定量评价中。
附图说明
图1是本发明白云岩储层胶结指数与有效孔隙度和导电孔隙度差值的关系图
图2是本发明白云岩储层有效孔隙度与导电孔隙度的关系图
图3是本发明白云岩储层的结果示意图
图4是本发明石灰岩储层胶结指数与有效孔隙度和导电孔隙度差值的关系图
图5是本发明石灰岩储层有效孔隙度与导电孔隙度的关系图
图6是本发明石灰岩储层的结果示意图
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
本发明方法是基于以下原理建立:由于胶结指数m值的物理意义是岩石导电路径的长度和岩石长度之比,因此岩石的导电路径与胶结指数m值之间存在较为密切的关系。而岩石的导电孔隙度在一定程度上能够刻画出岩石导电路径的长度,那么可以采用导电孔隙度来建立高精度的胶结指数m值的计算模型,又由于岩石导电孔隙度可以通过麦克斯韦(Maxwell)导电方程来进行计算。因此,本发明方法包括以下步骤:
3)基于胶结指数m值模型,建立油气储层饱和度定量评价方法。
下面分别加以详细描述:
雨果·弗瑞克(Hugo·Fricke)将麦克斯韦导电方程改写为更通用的形式:
式中:导电孔隙度是关于G值和地层因素F的函数,地层因素F可以通过岩电实验求出,因此,只要能求得G值即可以求出导电孔隙度而G值的获得与孔隙几何形状参数x有关,而对于复杂孔隙结构的岩石,x则是很难准确被计算的一个参数,从而导致导电孔隙度的准确计算一直是个难点。
式中:c1、c2均为系数。
联立式(9)和式(10),可以求解出系数c1和c2:
①通过岩电实验测量十二种岩性的地层因素F、胶结指数m值和有效孔隙度
②根据公式(13)
③根据公式(7)
式中:d1、d2均为系数。
3)基于胶结指数m值模型,建立油气储层饱和度定量评价方法将步骤2)建立的胶结指数m值模型带入阿尔奇公式:
便可以建立出复杂油气储层含水饱和度Sw的定量评价方法。
下面以碳酸盐岩中的白云岩、石灰岩储层作为实施例,阐述应用高精度胶结指数m值模型定量评价白云岩、石灰岩储层饱和度的过程。
实施例1:以某油田A井的白云岩储层为例
2)根据公式(13)
3)根据公式(7)
4)根据公式(14)
将式(16)代入式(15),可以得出某油田A井白云岩储层可变m值的计算公式:
5)将公式(17)中白云岩储层胶结指数m值的计算模型,代入阿尔奇公式,计算含水饱和度Sw:
本发明可以通过GeoFrame、Forward、Lead、Geolog等软件来计算该油田A井白云岩储层的含水饱和度Sw。
如图3所示,是采用GeoFrame软件绘制出的含水饱和度Sw,从图中可以看到:
在第七道,含水饱和度Sw的普通方法计算值是基于GeoFrame软件中ELAN-Plus模块最优化所得的含水饱和度,是目前计算复杂储层饱和度相对比较准确的一种方法;含水饱和度Sw的实验值是基于压汞实验分析所得的含水饱和度,其能最准确地反映储层的饱和度。从第七道可以清楚地看出,含水饱和度Sw的普通方法计算值和实验值之间存在一定的误差。
在第八道,胶结指数m值的实验值是基于岩电实验分析得出的;胶结指数m值的计算值是基于第五道有效孔隙度的准确计算,然后根据公式(17)计算的可变胶结指数m值。从第八道可以看出,胶结指数m值的计算值与实验值的符合程度较好,从而说明公式(17)计算的可变m值比较准确。
在第九道,含水饱和度Sw的计算值是基于第八道计算的可变m值,利用阿尔奇公式计算的含水饱和度;含水饱和度Sw的实验值是基于压汞实验分析所得出的,它能够最准确地反映储层的饱和度。从第九道可以看出,含水饱和度Sw的计算值和实验值之间的误差更小,说明基于该油田A井白云岩储层高精度m值模型所计算的饱和度比第七道使用其他的普通方法更准确,因此,改进了白云岩储层含水饱和度Sw的计算精度。
实施例2:以某油田B井的石灰岩储层为例
2)根据公式(13)
3)根据公式(7)
4)根据公式(14)
将式(19)代入式(18),可以得出石灰岩储层可变胶结指数m值的计算公式:
5)将公式(20)中石灰岩储层胶结指数m值的计算模型,代入阿尔奇公式,计算含水饱和度Sw:
如图6所示,是采用GeoFrame软件绘制出的含水饱和度Sw,从图中可以看到:
在第七道,含水饱和度Sw的普通方法计算值也是基于GeoFrame软件中ELAN-Plus模块最优化得出的;含水饱和度Sw的实验值也是基于压汞实验分析得出的。从第七道可以清楚地看出,含水饱和度Sw的普通方法计算值和实验值之间存在一定的误差。
在第八道,胶结指数m值的实验值是基于岩电实验分析得出的;胶结指数m值的计算值是基于第七道有效孔隙度的准确计算,然后根据公式(20)计算的可变胶结指数m值。从第八道可以看出,胶结指数m值的计算值与实验值的符合程度较好,从而说明公式(20)计算的可变m值比较准确。
在第九道,含水饱和度Sw的计算值是基于第九道计算的可变m值,利用阿尔奇公式计算的;含水饱和度Sw的实验值是基于压汞实验分析得出的。从第九道可以看出,含水饱和度Sw的计算值和实验值之间的误差更小,说明基于该油田B井石灰岩储层高精度m值模型所计算的饱和度比一般方法更准确,因此,改进了石灰岩储层含水饱和度Sw的计算精度。
上述各实例仅用于说明本发明,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。
Claims (6)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210514157.7A CN102979517B (zh) | 2012-12-04 | 2012-12-04 | 一种复杂油气储层饱和度定量评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210514157.7A CN102979517B (zh) | 2012-12-04 | 2012-12-04 | 一种复杂油气储层饱和度定量评价方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102979517A true CN102979517A (zh) | 2013-03-20 |
CN102979517B CN102979517B (zh) | 2015-06-17 |
Family
ID=47853834
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210514157.7A Active CN102979517B (zh) | 2012-12-04 | 2012-12-04 | 一种复杂油气储层饱和度定量评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102979517B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103527172A (zh) * | 2013-10-16 | 2014-01-22 | 中国石油集团川庆钻探工程有限公司 | 可变岩电耦合指数含水饱和度计算方法 |
CN103225506B (zh) * | 2013-04-19 | 2015-06-03 | 中国石油大学(华东) | 三组分自动混联导电饱和度模型的建立方法 |
CN104879126A (zh) * | 2015-05-14 | 2015-09-02 | 中国石油天然气集团公司 | 一种基于阵列感应测井的碎屑岩储层流体识别方法 |
CN105556345A (zh) * | 2013-07-18 | 2016-05-04 | 雪佛龙美国公司 | 用于估计地下储层中的孔隙度分布的系统和方法 |
CN106321087A (zh) * | 2015-07-06 | 2017-01-11 | 中国石油化工股份有限公司 | 一种获取岩石地层因素的方法 |
CN109386281A (zh) * | 2017-08-02 | 2019-02-26 | 中国石油化工股份有限公司 | 一种获取裂缝性低孔低渗储层测井饱和度的方法 |
CN112668814A (zh) * | 2019-09-30 | 2021-04-16 | 中国石油天然气股份有限公司 | 胶结等级确定方法及装置 |
CN113009115A (zh) * | 2021-03-11 | 2021-06-22 | 中海石油(中国)有限公司 | 一种多孔缝介质储层连通孔缝度的定量计算方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4751646A (en) * | 1986-05-15 | 1988-06-14 | Restech, Inc. | Method for determining original saturations in a produced field |
CN101042048A (zh) * | 2006-03-24 | 2007-09-26 | 中国石油天然气股份有限公司 | 复杂断块河流相储层油水井动用状况劈分系统 |
CN102565858A (zh) * | 2011-12-21 | 2012-07-11 | 西南石油大学 | 一种多孔介质含水饱和度的计算方法 |
-
2012
- 2012-12-04 CN CN201210514157.7A patent/CN102979517B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4751646A (en) * | 1986-05-15 | 1988-06-14 | Restech, Inc. | Method for determining original saturations in a produced field |
CN101042048A (zh) * | 2006-03-24 | 2007-09-26 | 中国石油天然气股份有限公司 | 复杂断块河流相储层油水井动用状况劈分系统 |
CN102565858A (zh) * | 2011-12-21 | 2012-07-11 | 西南石油大学 | 一种多孔介质含水饱和度的计算方法 |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103225506B (zh) * | 2013-04-19 | 2015-06-03 | 中国石油大学(华东) | 三组分自动混联导电饱和度模型的建立方法 |
CN105556345B (zh) * | 2013-07-18 | 2019-05-07 | 雪佛龙美国公司 | 用于估计地下储层中的孔隙度分布的系统和方法 |
CN105556345A (zh) * | 2013-07-18 | 2016-05-04 | 雪佛龙美国公司 | 用于估计地下储层中的孔隙度分布的系统和方法 |
CN103527172B (zh) * | 2013-10-16 | 2016-07-06 | 中国石油集团川庆钻探工程有限公司 | 可变岩电耦合指数含水饱和度计算方法 |
CN103527172A (zh) * | 2013-10-16 | 2014-01-22 | 中国石油集团川庆钻探工程有限公司 | 可变岩电耦合指数含水饱和度计算方法 |
CN104879126A (zh) * | 2015-05-14 | 2015-09-02 | 中国石油天然气集团公司 | 一种基于阵列感应测井的碎屑岩储层流体识别方法 |
CN104879126B (zh) * | 2015-05-14 | 2017-09-26 | 中国石油天然气集团公司 | 一种基于阵列感应测井的碎屑岩储层流体识别方法 |
CN106321087A (zh) * | 2015-07-06 | 2017-01-11 | 中国石油化工股份有限公司 | 一种获取岩石地层因素的方法 |
CN106321087B (zh) * | 2015-07-06 | 2019-05-07 | 中国石油化工股份有限公司 | 一种获取岩石地层因素的方法 |
CN109386281A (zh) * | 2017-08-02 | 2019-02-26 | 中国石油化工股份有限公司 | 一种获取裂缝性低孔低渗储层测井饱和度的方法 |
CN112668814A (zh) * | 2019-09-30 | 2021-04-16 | 中国石油天然气股份有限公司 | 胶结等级确定方法及装置 |
CN112668814B (zh) * | 2019-09-30 | 2024-03-08 | 中国石油天然气股份有限公司 | 胶结等级确定方法及装置 |
CN113009115A (zh) * | 2021-03-11 | 2021-06-22 | 中海石油(中国)有限公司 | 一种多孔缝介质储层连通孔缝度的定量计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102979517B (zh) | 2015-06-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102979517B (zh) | 一种复杂油气储层饱和度定量评价方法 | |
CN104007482B (zh) | 一种基于各向异性有效场的泥页岩岩石物理模型方法 | |
CN105804734A (zh) | 一种利用核磁共振测井识别稠油储层的方法 | |
CN105626009A (zh) | 一种缝洞型碳酸盐岩油藏单井注水替油效果定量评价方法 | |
CN103422853A (zh) | 一种确定地层含水饱和度的方法及装置 | |
CN103225506B (zh) | 三组分自动混联导电饱和度模型的建立方法 | |
CN105386753A (zh) | 一种利用核磁共振测井构造伪毛管压力曲线方法 | |
CN102540252B (zh) | 基于互相关的高精度中值叠加方法 | |
CN103778298B (zh) | 改进的模拟多孔介质中二维水流运动的多尺度有限元方法 | |
CN101135732A (zh) | 地层流体饱和度评价方法 | |
CN106019207A (zh) | 一种电能计量校准方法 | |
CN108019206B (zh) | 一种高介电常数下随钻电磁波电阻率仪器量程扩展方法 | |
CN103675460A (zh) | 基于相敏解调的电容耦合式非接触电导的测量装置及方法 | |
CN104912551A (zh) | 一种油水相渗曲线和驱油效率的标定方法 | |
CN103335996B (zh) | 一种镀锌板钝化膜膜重的测量方法 | |
CN104463345B (zh) | 一种预测地层油水界面位置的方法及装置 | |
CN101216564A (zh) | 基于小信号提取技术的早期试井分析方法 | |
王业飞 et al. | Study on sensitivities of numerical simulation parameters in polymer flooding | |
CN103527172A (zh) | 可变岩电耦合指数含水饱和度计算方法 | |
CN104374679A (zh) | 一种确定储层渗透率贡献值的方法 | |
CN109582923A (zh) | 一种基于双参数计算水合物饱和度的方法及处理终端 | |
CN107301255B (zh) | 一种计算储层裂缝孔隙度的方法 | |
CN103422852B (zh) | 一种不同井间气测值转换对比方法 | |
RU2015134392A (ru) | Способ моделирования подземного объема | |
CN104181090A (zh) | 一种基于岩电结构系数的含油气饱和度评价方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CP01 | Change in the name or title of a patent holder |
Address after: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Co-patentee after: CNOOC research institute limited liability company Patentee after: China Offshore Oil Group Co., Ltd. Address before: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Co-patentee before: CNOOC Research Institute Patentee before: China National Offshore Oil Corporation |
|
CP01 | Change in the name or title of a patent holder |