CN103867195A - 一种低渗透储层产水率的定量评价方法 - Google Patents
一种低渗透储层产水率的定量评价方法 Download PDFInfo
- Publication number
- CN103867195A CN103867195A CN201410078941.7A CN201410078941A CN103867195A CN 103867195 A CN103867195 A CN 103867195A CN 201410078941 A CN201410078941 A CN 201410078941A CN 103867195 A CN103867195 A CN 103867195A
- Authority
- CN
- China
- Prior art keywords
- water
- low permeability
- relation
- saturation
- water saturation
- 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.)
- Pending
Links
Images
Landscapes
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明涉及一种低渗透储层产水率的定量评价方法,包括以下步骤:1)建立低渗透储层的产水率Fw与含水饱和度Sw之间的关系,即选用神经元非线性函数Sigmoid表述产水率Fw与含水饱和度Sw之间的关系;2)钻取待开发的低渗透油田的某一深度下的岩心,对该岩心进行岩心实验获得Fw-Sw关系、地层水电阻率Rw、以及阿尔奇参数a、b、m和n值;并且根据实验得到的Fw-Sw关系辨识出步骤1)中产水率Fw与含水饱和度Sw关系式里的系数e、f和g从而确定出具体的Fw-Sw函数表达式;3)利用测井仪器测量得到低渗透油田的整个地层段岩石电阻率Rt和有效孔隙度并进一步算出整个地层段的含水饱和度Sw;4)根据步骤2)得到的Fw-Sw函数表达式以及步骤3)得到的整个地层段的含水饱和度Sw计算出整个地层段的产水率Fw。
Description
技术领域
本发明涉及一种石油开发领域的定量评价方法,特别是关于一种低渗透储层产水率的定量评价方法。
背景技术
储层产水率能准确划分油层的水淹级别,从而对油田储量和产量的准确评估具有至关重要的作用。储层产出流体的动态规律,主要服从多相流体在多孔介质微观孔隙中的分布与渗流特性,即遵从达西定律。对于油、水共渗体系,储层的产水率Fw可近似表示为式:
式中:Fw为产水率,无量纲;Qw为水的分流量,单位为m3;Qo为油的分流量,单位为m3;Krw为水的相对渗透率,无量纲;Kro为油的相对渗透率,无量纲;μw为水的粘度,单位为Pa·s;μo为油的粘度,单位为Pa·s;A为修正值,无量纲。
按照传统方法的思路,利用含水饱和度准确求取油水相对渗透率比值,是精确计算储层产水率的关键。根据Buckley-Leverrett的水驱理论和油藏数值模拟方法,含水饱和度Sw与油水相对渗透率比值Kro/Krw之间存在如下关系式(即它们之间呈现为指数关系):
式中:Sw为含水饱和度,无量纲;b、c为系数,无量纲。
综上,只要先利用含水饱和度Sw计算油水相对渗透率比值Kro/Krw,就可以求得储层的产水率Fw。
对于中、高渗透储层(储层渗透率主要范围为10~500mD),基于式(2)能准确建立含水饱和度Sw和油水相对渗透率比值Kro/Krw之间的函数关系,也就能准确计算储层的产水率Fw。然而,对于孔隙空间狭小、喉道非常细小、孔隙连通性差的低渗透储层(储层渗透率分布的主要范围小于1mD),因启动压力梯度和强应力敏感性的存在,使得流体在低渗透储层中的流动规律与中、高渗透储层存在较大差异。流体在低渗透储层中的流动规律并不遵循达西定律,从而使得含水饱和度Sw与油水相对渗透率比值Kro/Krw之间的指数函数关系并不存在。
发明内容
针对上述问题,本发明的目的是提供一种能够准确计算储层产水率Fw的低渗透储层产水率定量评价方法。
为实现上述目的,本发明采用以下技术方案:一种低渗透储层产水率的定量评价方法,包括以下步骤:
1)建立低渗透储层的产水率Fw与含水饱和度Sw之间的关系,即选用神经元非线性函数Sigmoid表述产水率Fw与含水饱和度Sw之间的关系:
式中:e、f和g为Sigmoid函数的系数;
2)钻取待开发的低渗透油田的某一深度下的岩心,对该岩心进行岩心实验获得Fw-Sw关系、地层水电阻率Rw、以及阿尔奇参数a、b、m和n值;并且根据实验得到的Fw-Sw关系辨识出步骤1)中产水率Fw与含水饱和度Sw关系式里的系数e、f和g从而确定出具体的Fw-Sw函数表达式;
式中:a、b为与岩性有关的比例系数,无量纲;m为胶结指数,无量纲;n为饱和指数,无量纲;Rw为地层水的电阻率,单位为Ω·m;Rt为岩石电阻率,单位为Ω·m;为有效孔隙度,无量纲;其中的a、b、m和n统称为阿尔奇参数;
4)根据步骤2)得到的Fw-Sw函数表达式和步骤3)得到的整个地层段的含水饱和度Sw,计算出整个地层段的产水率Fw。
本发明由于采取以上技术方案,其具有以下优点:1、本发明方法通过某一深度下的岩心实验来建立整个地层段的产水率与含水饱和度的关系,因此不用开展大量的岩心实验,能够有效地节约成本,具有较强的经济性。2、本发明方法选取神经元非线性函数Sigmoid来描述低渗透储层产水率Fw和含水饱和度Sw之间的突变关系,使得本发明方法与低渗透储层的渗流机理更为吻合,也保证了定量评价结果的准确性。
附图说明
图1是某低渗透油田岩心实验得到的产水率与含水饱和度之间的关系;
图2是某低渗透油田辨识参数得到的产水率与含水饱和度的拟合曲线;
图3是某低渗透油田计算产水率与实际产水率的对比图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
本发明一种低渗透储层产水率的定量评价方法,包括以下步骤:
1)建立低渗透储层的产水率Fw与含水饱和度Sw之间的关系
由于低渗透储层特殊的渗流机理,含水饱和度Sw与油水相对渗透率比值Kro/Krw之间理想的指数函数关系并不存在。但是可以肯定的是,油、水的相对渗透率Kro、Krw与含水饱和度Sw之间均存在一定的函数关系,即:
Kro=f(Sw),Krw=f(Sw)(3)
那么由式(1)可知,产水率Fw与含水饱和度Sw之间也一定存在某种函数关系,即:
低渗透储层由于启动压力梯度和强应力敏感性的存在,含水饱和度Sw和产水率Fw之间的关系呈现为“突变”,因此可选用神经元非线性函数Sigmoid来描述这一突变特征,Sigmoid函数的表达式如下:
通过改变Sigmoid函数的系数,可以控制其突变的时间和幅度,建立产水率Fw与含水饱和度Sw的关系如下:
式中:e、f和g为Sigmoid函数的系数;通过改变系数e、f和g,Sigmoid函数便能准确表征含水饱和度Sw和产水率Fw的突变现象。
2)钻取待开发的低渗透油田的某一深度下的岩心,对该岩心进行岩心实验获得Fw-Sw关系(如图1所示)、地层水电阻率Rw、以及阿尔奇参数a、b、m和n值;并且根据实验得到的Fw-Sw关系辨识出步骤1)中产水率Fw与含水饱和度Sw关系式里的系数e、f和g从而确定出具体的Fw-Sw函数表达式。
式中:a、b为与岩性有关的比例系数,无量纲;m为胶结指数,无量纲;n为饱和指数,无量纲;Rw为地层水的电阻率,单位为Ω·m;Rt为岩石电阻率,单位为Ω·m;为有效孔隙度,无量纲。其中的a、b、m和n统称为阿尔奇参数。
4)根据步骤2)得到的Fw-Sw函数表达式和步骤3)得到的整个地层段的含水饱和度Sw,计算出整个地层段的产水率Fw。
下面是一具体实施例:对某个待开发的低渗透油田的产水率进行定量评价。
首先建立低渗透储层的产水率Fw与含水饱和度Sw之间的关系:
通过岩心实验,得到了该油田的阿尔奇参数:a=9.75、m=0.87、b=1.19、n=1.61,地层水电阻率Rw=0.11Ω·m,以及如图1所示的Fw-Sw关系;对系数e、f和g进行辨识,从而确定出具体的Fw-Sw函数表达式(如图2所示):
Fw=1/(0.02856×exp(-0.5319×(Sw-50))+0.01023),R2=0.96
式中:R2为相关系数。
最后利用Fw=1/(0.02856×exp(-0.5319×(Sw-50))+0.01023)计算出该油田的整个地层段的产水率。
在该具体实施例中,还对上述方法计算得到的产水率与实际产水率进行了对比(如图3所示),在该低渗透油田,计算产水率与实际产水率之间的相关系数R2约为0.95,平均绝对误差约为2.08%,产水率的计算精度相对较高,满足该油田低渗透储层的勘探开发生产需求。
上述各实施例仅用于说明本发明,其中方法的实施步骤等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。
Claims (1)
1.一种低渗透储层产水率的定量评价方法,包括以下步骤:
1)建立低渗透储层的产水率Fw与含水饱和度Sw之间的关系,即选用神经元非线性函数Sigmoid表述产水率Fw与含水饱和度Sw之间的关系:
式中:e、f和g为Sigmoid函数的系数;
2)钻取待开发的低渗透油田的某一深度下的岩心,对该岩心进行岩心实验获得Fw-Sw关系、地层水电阻率Rw、以及阿尔奇参数a、b、m和n值;并且根据实验得到的Fw-Sw关系辨识出步骤1)中产水率Fw与含水饱和度Sw关系式里的系数e、f和g从而确定出具体的Fw-Sw函数表达式;
式中:a、b为与岩性有关的比例系数,无量纲;m为胶结指数,无量纲;n为饱和指数,无量纲;Rw为地层水的电阻率,单位为Ω·m;Rt为岩石电阻率,单位为Ω·m;为有效孔隙度,无量纲;其中的a、b、m和n统称为阿尔奇参数;
4)根据步骤2)得到的Fw-Sw函数表达式和步骤3)得到的整个地层段的含水饱和度Sw,计算出整个地层段的产水率Fw。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410078941.7A CN103867195A (zh) | 2014-03-05 | 2014-03-05 | 一种低渗透储层产水率的定量评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410078941.7A CN103867195A (zh) | 2014-03-05 | 2014-03-05 | 一种低渗透储层产水率的定量评价方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103867195A true CN103867195A (zh) | 2014-06-18 |
Family
ID=50906157
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410078941.7A Pending CN103867195A (zh) | 2014-03-05 | 2014-03-05 | 一种低渗透储层产水率的定量评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103867195A (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105587297A (zh) * | 2014-10-23 | 2016-05-18 | 中国石油化工股份有限公司 | 复杂断块油藏仿强边水驱技术适应性定量评价方法 |
CN106897531A (zh) * | 2017-03-06 | 2017-06-27 | 中海石油(中国)有限公司 | 一种低渗透石灰岩储层渗透率的定量评价方法 |
CN106970000A (zh) * | 2017-04-21 | 2017-07-21 | 西南石油大学 | 煤/页岩超高压气体吸附和渗流实验评价页岩气吸附方法 |
CN107631974A (zh) * | 2017-09-07 | 2018-01-26 | 中国海洋石油总公司 | 一种致密砂岩蜡封岩心含水饱和度的获取方法 |
CN107725044A (zh) * | 2017-09-27 | 2018-02-23 | 中国石油天然气集团公司 | 基于阵列感应、侧向测井的砂岩含气储层产水率预测的方法 |
CN111350497A (zh) * | 2020-03-23 | 2020-06-30 | 中国石油天然气集团有限公司 | 一种基于含水率的低饱和度油层动态评价方法 |
CN112392464A (zh) * | 2020-12-11 | 2021-02-23 | 中国石油天然气集团有限公司 | 一种基于常规测井资料计算储层产水率的方法 |
CN112800646A (zh) * | 2021-01-08 | 2021-05-14 | 中海石油深海开发有限公司 | 一种修正阿尔奇模型中岩电参数的方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1111002A (zh) * | 1994-04-22 | 1995-11-01 | 郑金安 | 用普通取心—测试资料求水淹油层饱和度的方法 |
CN101649738A (zh) * | 2008-08-13 | 2010-02-17 | 中国石油天然气集团公司 | 一种确定地层含水饱和度的方法 |
CN101942994A (zh) * | 2010-09-16 | 2011-01-12 | 中国石油天然气股份有限公司 | 水淹层产水率定量预测方法及其系统 |
CN102720485A (zh) * | 2012-07-20 | 2012-10-10 | 中国石油天然气股份有限公司 | 一种中高含水油田含水上升率的预测方法 |
CN103225500A (zh) * | 2013-05-02 | 2013-07-31 | 中国石油大学(华东) | 一种应用三参数自洽迭代算法的新型水淹层测井评价方法 |
CN103233726A (zh) * | 2013-05-02 | 2013-08-07 | 中国石油大学(华东) | 一种水淹层饱和度、产水率测井评价模型的实验刻度方法 |
-
2014
- 2014-03-05 CN CN201410078941.7A patent/CN103867195A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1111002A (zh) * | 1994-04-22 | 1995-11-01 | 郑金安 | 用普通取心—测试资料求水淹油层饱和度的方法 |
CN101649738A (zh) * | 2008-08-13 | 2010-02-17 | 中国石油天然气集团公司 | 一种确定地层含水饱和度的方法 |
CN101942994A (zh) * | 2010-09-16 | 2011-01-12 | 中国石油天然气股份有限公司 | 水淹层产水率定量预测方法及其系统 |
CN102720485A (zh) * | 2012-07-20 | 2012-10-10 | 中国石油天然气股份有限公司 | 一种中高含水油田含水上升率的预测方法 |
CN103225500A (zh) * | 2013-05-02 | 2013-07-31 | 中国石油大学(华东) | 一种应用三参数自洽迭代算法的新型水淹层测井评价方法 |
CN103233726A (zh) * | 2013-05-02 | 2013-08-07 | 中国石油大学(华东) | 一种水淹层饱和度、产水率测井评价模型的实验刻度方法 |
Non-Patent Citations (3)
Title |
---|
梅俊伟 等: "利用Sigmoid函数表征超低渗透储层的非达西渗流特征", 《西安石油大学学报(自然科学版)》, vol. 26, no. 2, 25 March 2011 (2011-03-25) * |
王敏 等: "低闭合油藏油水同层储层含水饱和度新模型", 《西南石油大学学报(自然科学版)》, vol. 32, no. 1, 15 February 2010 (2010-02-15) * |
许长福 等: "岩性油藏特征制约下超低渗透油层快速识别方法与模型", 《中南大学学报(自然科学版)》, vol. 43, no. 1, 26 January 2012 (2012-01-26) * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105587297A (zh) * | 2014-10-23 | 2016-05-18 | 中国石油化工股份有限公司 | 复杂断块油藏仿强边水驱技术适应性定量评价方法 |
CN105587297B (zh) * | 2014-10-23 | 2017-12-15 | 中国石油化工股份有限公司 | 复杂断块油藏仿强边水驱技术适应性定量评价方法 |
CN106897531A (zh) * | 2017-03-06 | 2017-06-27 | 中海石油(中国)有限公司 | 一种低渗透石灰岩储层渗透率的定量评价方法 |
CN106897531B (zh) * | 2017-03-06 | 2020-03-27 | 中海石油(中国)有限公司 | 一种低渗透石灰岩储层渗透率的定量评价方法 |
CN106970000A (zh) * | 2017-04-21 | 2017-07-21 | 西南石油大学 | 煤/页岩超高压气体吸附和渗流实验评价页岩气吸附方法 |
CN106970000B (zh) * | 2017-04-21 | 2024-03-22 | 西南石油大学 | 煤/页岩超高压气体吸附和渗流实验评价页岩气吸附方法 |
CN107631974B (zh) * | 2017-09-07 | 2019-11-12 | 中国海洋石油集团有限公司 | 一种致密砂岩蜡封岩心含水饱和度的获取方法 |
CN107631974A (zh) * | 2017-09-07 | 2018-01-26 | 中国海洋石油总公司 | 一种致密砂岩蜡封岩心含水饱和度的获取方法 |
CN107725044A (zh) * | 2017-09-27 | 2018-02-23 | 中国石油天然气集团公司 | 基于阵列感应、侧向测井的砂岩含气储层产水率预测的方法 |
CN107725044B (zh) * | 2017-09-27 | 2020-09-04 | 中国石油天然气集团公司 | 基于阵列感应、侧向测井的砂岩含气储层产水率预测的方法 |
CN111350497A (zh) * | 2020-03-23 | 2020-06-30 | 中国石油天然气集团有限公司 | 一种基于含水率的低饱和度油层动态评价方法 |
CN112392464A (zh) * | 2020-12-11 | 2021-02-23 | 中国石油天然气集团有限公司 | 一种基于常规测井资料计算储层产水率的方法 |
CN112800646A (zh) * | 2021-01-08 | 2021-05-14 | 中海石油深海开发有限公司 | 一种修正阿尔奇模型中岩电参数的方法 |
CN112800646B (zh) * | 2021-01-08 | 2024-01-23 | 中海石油深海开发有限公司 | 一种修正阿尔奇模型中岩电参数的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103867195A (zh) | 一种低渗透储层产水率的定量评价方法 | |
CN106837297B (zh) | 一种识别井间连通性及油水动态预测的方法 | |
CN103352693B (zh) | 一种基于含油气含量的产能预测系统及其方法 | |
CN104847340B (zh) | 一种水淹层测井定量评价方法 | |
CN104834807B (zh) | 一种基于分形理论的应力敏感储层相对渗透率计算方法 | |
CN102789003B (zh) | 利用阵列感应测井资料进行地层参数反演的方法及装置 | |
CN103485770B (zh) | 一种基于人工神经网络求取含油饱和度的方法及系统 | |
CN104695950A (zh) | 火山岩油藏产能预测方法 | |
CN104899411B (zh) | 一种储层产能预测模型建立方法和系统 | |
CN104100263B (zh) | 一种非均匀水淹层剩余油饱和度的确定方法 | |
CN107038268A (zh) | 一种确定非均质储层五点井网水驱波及系数的方法 | |
CN104453874A (zh) | 一种基于核磁共振的砂砾岩储层含油饱和度的计算方法 | |
CN104343445B (zh) | 识别气层与水层、高气油比储层及含气饱和度的方法 | |
CN103670390A (zh) | 一种水淹层测井评价方法及系统 | |
CN105525917A (zh) | 识别储层的方法以及识别储层中流体类型的方法 | |
CN106154343A (zh) | 计算致密油储层的含油饱和度的方法 | |
CN106897531A (zh) | 一种低渗透石灰岩储层渗透率的定量评价方法 | |
CN102536223B (zh) | 煤层裂隙孔隙度快速计算方法 | |
CN104237107A (zh) | 地层中低渗透率储层的视渗透率解释方法及系统 | |
CN105844011A (zh) | 一种基于毛管模型的渗透率计算方法 | |
CN103615239B (zh) | 一种测井地层成像系统及其方法 | |
CN104122182B (zh) | 矿井储层有效厚度下限的获取方法 | |
CN104407397A (zh) | 一种致密油气储层物性下限的确定方法 | |
CN105003258A (zh) | 一种高温高压气层甲烷流体密度骨架参数的获取方法 | |
CN104675391A (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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20140618 |