基于J函数的含油饱和度建模方法—以陆丰A油田为例

刘萝卜锅

品质源于技术 服务源于态度

这是阿什卡微信公众号的第471篇原创文章

首发于2018年1月14日

▽

作者

戴宗1，朱义东1，李威1，居字龙1，徐辰2

1. 中海石油（中国）有限公司深圳分公司，广东深圳，518000

2. 阿什卡能源科学咨询服务联合公司，北京，100101

Dai Zong1，Zhu Yidong1，Li Wei1，Ju Zilong1，Xu Chen2

1.Institute of CNOOC Shenzhen Branch； Shenzhen, Guangdong 518000, China

2.Essca EnergyScience & Consulting Associate，Beijing 100101，China

作者简介：戴宗，高级工程师，主要从事油气田开发和油田地质研究工作。地址：广东深深圳市后海滨路3168号中海油大厦（邮编518000），Email：daizong@cnooc.com.cn.项目来源：中海石油（中国）有限公司深圳分公司自立科研项目。项目名称：南海东部在生产油田挖潜技术研究。项目编号：YXKY-ZL-SZ-2016-04

摘要：随着注水开发的进行，油层各项参数随之发生变化。由于水淹的影响，测井解释的含油饱和度与油藏原始含油饱和度在数值和分布上都有显著差异，不能代表油藏的原始分布状态。依据平均化后的毛管压力曲线，分流动单元对各类型储层进行原始含油饱和度的计算，建立基于J函数的含油饱和度模型，减小不同钻井时期测井解释含油饱和度的误差，对储量复算、油藏数值模拟等工作具有重要意义。

关键词：陆丰A油田  J函数  饱和度模型测井毛管压力曲线

1.  引言

随着油田注水开发的不断进行，油层的孔隙结构发生改变，含油饱和度发生变化。一般表现为泥质含量下降，孔隙度、渗透率增加以及电阻率的下降。

因此，对开发后期的新钻井来说，测井解释含油饱和度不能代表油层原始含油饱和度，建立在测井解释基础上的含油饱和度模型不能代表油藏原始的分布状态。

图1 第7类样品毛管压力曲线

传统饱和度建模一般基于饱和度解释曲线采用空间插值算法完成，该方法可靠性很大程度上依赖于控制点数量，因此对于原始饱和度建模来说只能用少量的早期未水淹井做为控制[1]。

而且插值方法无法考虑油水界面及毛管力的影响，无法准确表征油水过渡带。利用毛管压力计算的原始含油饱和度只与储层孔隙结构、油柱高度、油水密度差有关，与电阻率无关[2-4]。

基于Ｊ函数的饱和度建模方法主要使用受流体影响较小的孔隙度和渗透率属性场，因而控制点数更多。

另外，该方法最大的特点是考虑毛细管力所引起的流体过渡带，能够为数值模拟提供近似平衡的饱和度场，更有利于后续的数值模拟研究。

2. 基本原理

2.1 利用J函数计算原始含油饱和度

油藏形成过程实际上是石油浮力不断克服毛管压力而发生垂向或侧向运移，在遇到遮挡层而在储集岩聚集而成[8]。石油浮力与油柱高度及油水密度差有关，毛管压力则与界面角、接触张力和孔喉结构有关(式1，2)[5-6]。两者的相对大小决定了油气是否能进入储层孔隙，也就决定了油藏的原始分布状态(式3)[7]。

当油气充注结束，毛管压力与浮力达到平衡时，有以下关系：

式中，F浮为石油浮力，N；△ρ为油水密度差，g/cm3；Pc为毛细管压力，MPa；σ为流体两相的界面张力，mN/m；θ为流体与固体的接触角，°；r为毛细管半径，μm；Ｈ为油柱高度，ｍ；Pc为毛细管压力，Mpa。

利用这个基本原理可以根据毛管压力曲线计算油藏的原始含油饱和度。岩块样品的毛管压力曲线只能代表该点上的油藏特征，在实际研究过程中，为了得到整个油藏范围内各类储层条件下的油藏特征曲线，需要对多块岩心样品进行平均化处理。

Leverett于1941年提出无因次J函数(式4，5)[8]，拟合Swn与J(Swn)以得到代表各类储层条件的毛管压力曲线，从而实现对原始含油饱和度的计算。

式中Sw，岩心含水饱和度，%；Swi，岩心束缚水饱和度，%；K，岩心水平渗透率，mD；φ，岩心有效孔隙度，%。

实际储层中，由于非均质性的影响，一条毛管压力曲线不能代表储层每个点的性质，因此分类前需要根据物性和储集层质量系数指数(Irq)(式6)[9]对样品进行分类[10]。Irq可以定性反映储层孔隙结构，Irq值越大，储层孔隙结构越好，孔喉半径越大。利用Irq指数分类后，每类样品具有类似的孔隙结构，因此也可以同一毛管压力曲线代表。

式中Irq，储集层质量系数指数，无量纲；K和φ同上。

获得平均毛管压力曲线后，将实验室条件下的毛管压力曲线转变为储层条件下的毛管压力曲线(式7)[11]。

式中Pc，毛管压力，Mpa；σ，界面张力，dyn/cm；θ，润湿接触角，°；Pc(res)表示油藏条件下的毛管压力值，Mpa；Pc(lab)表示实验室条件下的毛管压力测量值，Mpa。

由于平衡时毛管压力与浮力相等，对于已知的任意高度，可以计算其原始含油饱度。

2.2 基于J函数的含油饱和度建模

根据上述理论，在孔隙度、渗透率模型的基础上，结合流体界面深度认识及界面张力、润湿角等参数，可以完成基于J函数的含油饱和度建模。

对于开发成熟的油气田，利用Ｊ函数建立的含油饱和度模型，更符合油气田开发实际，可以更真实地反映出地层的原始含油(气)状态，为储量复算和油藏数值模型提供准确依据。

刘萝卜锅

3. 实例运用

陆丰A油田位于珠江口盆地北部坳陷带，为一平缓、完整的背斜构造。主力油层2500层储层岩性主要为块状砂岩夹泥岩。储层物性好，产量高。岩心分析孔隙度3.6%-26.17%，平均19.16%；渗透率0.2-8843.9mD，平均134mD。经过20多年的投产开发，处于特高含水的挖潜后期。

本次研究收集到岩心毛管压力资料共3口井47块。Irq值为0.11-8.32，按照聚类分析结果将其划分为11类，个别偏离较大的样品会造成J函数拟合相关系数不高，删除后对各类样品进行拟合(图1，图2)，得到11类J函数(表1)。

将J函数代入式(式4)[8]，得到实验室条件下11类平均毛管压力曲线[12]。

此次收集到的毛管压力所用实验方法为半渗透隔板法，实验介质为空气-盐水。半渗透隔板法与压汞法相比，实验系统更接近地下油藏实际情况，测试精度更高，数据更可靠。将得到的毛管压力曲线从实验室条件转换为地下油藏条件(式7)[11]，σ和θ的取值参考了标准值以及邻区陆丰B油田的实验数据，取σlabcosθlab为72，σrescosθres为28。得到11类平均毛管压力曲线(图3)。从第1类到第11类毛管压力曲线的启动压力、饱和度中值压力逐渐增高，平缓段变短，非润湿相饱和度增加。这说明随着Irq减小，样品物性逐渐变差，孔喉半径逐渐减小。

图2 第7类样品J函数拟合

据高压物性分析、地层水分析资料，取ρw=1.02g/cm3，ρo=0.87g/cm3，H为所在深度距自由水面的距离(m)，将得到的毛管压力曲线代入式3，计算得到任意高度的原始含水饱和度。

实际计算中，依据测井解释的孔隙度、渗透率对数据点进行分类，按照所属类别计算该数据点的J函数、毛管压力曲线，并最终得到原始含水饱和度。

图3 平均毛管压力曲线

图4为计算结果，从左到右，第一道为自然伽马曲线、钻头尺寸、井径尺寸；第二道为层位、第三道为中子、密度曲线、第四道为测井解释孔隙度、第五道为J函数计算含水饱和度及测井解释饱和度对比。蓝色曲线为印尼公式计算的含水饱和度，枚红色曲线为J函数计算原始含水饱和度。早期钻井(左侧)J含水计算饱和度与测井解释饱和度吻合较好，随着开发的进行，开发中晚期钻井(右侧)测井饱和度较原始含水饱和度明显增加，油水界面随之上升。

图4 不同时期钻井J函数计算含水饱和度与测井解释含水饱和度对比

本次研究，基于J函数的含油饱和度建模结果见图5，主要过程如下：

(1)以相关反演结果为约束条件，采用续贯高斯方法模拟孔隙度分布。

(2)以孔隙度模型为约束条件，采用续贯高斯方法模拟渗透率分布。

(3)在孔隙度、渗透率模型的基础上，根据上述11类储层分类标准，计算储层分类模型。

(4)结合流体界面深度认识及界面张力、润湿角等参数，分不同储层类型完成基于J函数的饱和度建模。

图5 陆丰A区块主要参数地质模型

从所建含水饱和度模型可以看出：在油水界面之上，高孔渗区对应高含油饱和度，与不同储层类型所对应的相渗端点值一致，体现了物性好的储层束缚水饱和度相对较低的规律；在油水界面上方，不同储层类型的具有不同的油水过渡带高度。上述模型特点更符合油藏实际情况，而传统饱和度插值方法是难以实现的。

4. 主要结论

1) 通过J函数计算原始含油饱和度的方法，计算的结果只与储层孔隙结构、油水密度、油柱高度有关，与电阻率无关，从而解决了水淹层测井解释含油饱和度不能代表原始含油饱和度的问题；利用储层质量系数分类的方法，充分考虑了储层的孔隙结构差异，准确地计算各类储层的原始水饱和度，在非均质性强的储层中得到了很好的运用。

2) 基于Ｊ函数方法建立的含油饱和度模型，与不同岩石类型相渗一致，能够体现不同岩石类型毛细管力所引起的流体过渡带，与传统插值方法所建立的含油饱和度模型相比，能够更加真实地反映油气藏原始含油饱和度分布。

参考文献：

[1] 程超，蔡华，廖恒杰等. 海上油气田J函数约束饱和度地质建模方法[J]. 石油地质与工程,2016,30(1):103-105.

[2] 陈康，赵新军. 确定储层含油饱和度的测井~毛管资料综合法[J].江汉石油学院学报,1995,17(2):51-55.

[3] 赵国欣，朱家俊，关丽. 用毛管压力资料求取原始含油饱和度的方法[J]. 中国石油大学学报(自然科学版), 2008,32(4):38-41.

[4] 唐衔，侯加根，许凡，等.利用压汞曲线求取油藏原始含油饱和度方法的研究~以克拉玛依油田为例[J]. 吐哈油气, 2008,13(4):366-399.

[5] 邵才瑞, 张鹏飞, 张福明, 等. 用J函数提高致密砂岩气层饱和度测井评价精度[J]. 中国石油大学学报(自然科学版), 2015, 40(4): 57-65.

[6] 廖敬, 彭彩珍, 吕文均, 等. 毛管压力曲线平均化及J函数处理[J]. 特种油气藏, 2008, 15(6): 73-75.

[7] 卢艳，水淹层原始含水饱和度评价方法研究[J]. 国外测井技术,2015,(1):16-19.

[8] Leverett M C. Capillary behavior in porous solids[J].Transactions of the AIME, 1941, 142(1): 151-169.

[9] 司马立强, 李清, 杨毅, 等. 用 J 函数法求取碳酸盐岩储层饱和度方法探讨[J]. 岩性油气藏, 2014, 26(6): 106-110.

[10] 胡勇, 于兴河, 陈恭洋, 等. 平均毛管压力函数分类及其在流体饱和度计算中的应用[J].石油勘探与开发, 2012,39(6): 733-738.

[11] 章成广，秦瑞宝. 用毛管压力曲线解释原始含水饱和度[J]. 江汉石油大学学报,1999,21(4):8-10.

[12] 管耀. 利用毛管压力资料求原始含油饱和度方法探讨[J]. 勘探地球物理进展,2009,32(5):365-369.

相关文章

• 国际权威期刊论文：这样解决缝洞型碳酸盐岩油气藏水力压裂评价难题
• 国际权威期刊论文｜Dionisos解决渤海湾凹陷砂体描述难题
• 国际权威期刊论文｜关于沉积正演数值模拟研究
• 国际期刊论文｜渤海湾三维精细储层地质模型实例
• 用户论文 | 微注压降测试在页岩气井的应用
  
  

-END-

（关注微信号 每期早看到↘）

--------------版权声明--------------

阿什卡原创文章

转发请注明出处

--------------版权声明--------------

郑漂一数模

刘萝卜锅

郑漂一数模 发表于 2018-2-11 16:04

等了那么久，您终于同意了，谢谢