基于吸水剖面资料的油藏层间平均剩余油饱和度计算方法

张世明，吴晓东，李 坤，崔传智

（1.中国石油大学（北京）石油工程学院，北京102249；2.中国石化胜利油田分公司地质科学研究院，山东东营257015；3.中国石油大学（华东）石油工程学院，山东青岛266580）

基于吸水剖面资料的油藏层间平均剩余油饱和度计算方法

张世明1，2，吴晓东1，李 坤3，崔传智3

（1.中国石油大学（北京）石油工程学院，北京102249；2.中国石化胜利油田分公司地质科学研究院，山东东营257015；3.中国石油大学（华东）石油工程学院，山东青岛266580）

层间采出程度差异大是多层合采油藏在高含水期面临的主要问题之一，正确认识层间采出程度状况对开发措施的制定具有重要的指导意义。吸水剖面资料是反映层间吸水状况差异的重要信息，分层吸水量的差异主要受各小层的储层物性和剩余油饱和度的差异影响。利用吸水剖面资料，以水电相似原理和非活塞式水驱油理论为基础，建立了多层合采油藏各小层平均剩余油饱和度的计算方法，并将计算结果与油藏数值模拟结果进行对比，两者的相对误差小于5%。

多层油藏层间矛盾吸水剖面渗流阻力剩余油饱和度

陆相沉积油藏纵向层系多，储层物性差异大，通常采用多个小层组合成一套层系进行开发，层间采出程度差异较大。油田开发措施制定的重点之一就是正确认识层间采出程度差异状况，主要包括室内物理实验、油藏数值模拟和矿物监测等方法。近年来中国学者先后开展了基于室内物理实验或矿场监测方法获取层间动用状况的研究，马圣贤等通过室内实验对桥口油田多层水驱油藏进行研究，认识到储层物性差异致使层间具有不同的动用程度［1］；崔英琢等利用监测资料，从平面、层间及储层特征3个方面分析水驱油藏动用状况变化规律及其原因［2］；冯其红等以单层注采平衡为依据，综合考虑地质和开发条件，分阶段预测单层剩余油饱和度［3］；李淑霞等研究了利用井间示踪剂确定剩余油饱和度的方法［4］。

油藏数值模拟是非均质油藏剩余油定量化研究的重要技术手段，中外学者多采用油藏数值模拟技术对剩余油分布特征进行研究［5-8］。对比分析表明，利用饱和度测井、产液剖面和取心井资料等矿场测试资料分析纵向剩余油分布特征直观准确可靠，但取心井资料少，仅能分析井点纵向剩余油分布；采用油藏数值模拟求取剩余油饱和度能够反映多种因素的影响，但由于该方法需要大量的基础数据准备和专业的油藏动态历史拟合技术，且研究周期长，矿场普及应用难度较大。

吸水剖面测试工艺简单易行，是反映层间吸水状况以及调剖堵水措施效果而采取的常规监测手段。分层的吸水状况直接反映了层间的吸水能力，也间接反映了储层物性及采出程度。目前利用吸水剖面资料确定纵向剩余油分布的研究较少，为此，笔者利用吸水剖面资料，建立了求取各小层平均剩余油饱和度的计算方法，以期为深化认识水驱油藏剩余油分布及制定合理的开发调整措施提供依据。

1 分层平均剩余油饱和度计算方法

1.1 基本思路

以注水井为中心，并将注水井及其所对应的油井视为一个注采井组，计算中将该井组作为一个整体来考虑，利用注水井的吸水剖面资料求取井组各层平均剩余油饱和度。为简单起见，将注采井组等效为一注一采的注采单元，其注采井距为井组中各注采井距的平均值，按照注采单元的控制面积与井组的控制面积相等原则，得到该单元的宽度；根据井组中各油井产水量和产油量计算得到油井的含水率。

根据各小层储层物性和吸水剖面求取其平均剩余油饱和度的主要步骤包括：①将n个合采小层作为一个整体，根据油井含水率求得地层中的平均剩余油饱和度，结合相对渗透率曲线求得地层的总渗流阻力；②应用水电相似原理和水井的总注水量得到单元的等效生产压差；③再根据等效生产压差和每个小层的吸水量求得每个小层的渗流阻力，继而求得各小层油水相的视粘度；④根据建立的视粘度与平均含水饱和度的关系曲线，即可求得各小层的平均剩余油饱和度。

1.2 注采单元相关数据的求取

1.2.1 储层参数平均值

多层合采时的单元总流量等于各小层流量之和，按照等效渗流的原理，其注采单元的平均有效厚度、平均渗透率和平均孔隙度计算式分别为

式中：-h为注采单元的平均有效厚度，cm；i为小层序号；n为总层数；hi为第i小层的有效厚度，cm；-K为注采单元的平均渗透率，μm2；Ki为第i小层的渗透率，μm2；-φ为注采单元的平均孔隙度；φi为第i小层的孔隙度。

1.2.2 平均含水饱和度

根据已知相对渗透率曲线建立注采单元储层含水率与含水饱和度的关系曲线，并由油井含水率可得到油井处的含水饱和度。根据水驱油时岩心平均含水饱和度和岩心出口端含水饱和度之间的关系［9-11］可求得注采单元中平均含水饱和度为

中华民族“站起来”“富起来”“强起来”的奋斗目标和民族伟大复兴的共同理想，把爱国主义和拥护中国共产党的领导、坚持和发展中国特色社会主义内在地统一起来。正如习近平总书记所指出的那样：“我国爱国主义始终围绕着实现民族富强、人民幸福而发展，最终汇流于中国特色社会主义。祖国的命运和党的命运、社会主义的命运是密不可分的。”[注]习近平：《习近平在中共中央政治局第二十九次集体学习时强调：大力弘扬爱国主义精神 为实现中国梦提供精神支柱》，《人民日报》2015年12月31日，第1版。

在已知油井含水饱和度的条件下，利用注采单元储层含水率与含水饱和度的关系曲线，通过作图法或者迭代法求得该饱和度下的含水率。根据式（4）可得到地层的平均含水饱和度。

1.2.3 等效生产压差

注采单元等效生产压差并不是井组的实际生产压差，而是依据注采单元的渗流阻力和注水量计算得到的生产压差。

根据油水相对渗透率曲线，可得到平均含水饱和度下相应的相对渗透率，进而求取注采单元的总渗流阻力为

式中：Rt为注采单元总渗流阻力，MPa/（cm3· s-1）；L为注采单元的长度，cm；B为注采单元的宽度，cm；Kro为油相相对渗透率；μo为地层原油粘度，mPa·s；Krw为水相相对渗透率；μw为地层水粘度，mPa·s。

根据水电相似原理，由注水井总注水量可得到注采单元的等效生产压差为

式中：Δp为等效生产压差，10-1MPa；Qt为总注水量，cm3/s。

1.2.4 分层平均含水饱和度

根据水电相似原理，多层合采时各小层的生产压差相同，皆等于注采单元的等效生产压差，由此求出各小层的渗流阻力为

式中：Ri为第i小层的渗流阻力，MPa/（cm3· s-1）；Qi为第i小层的吸水量，cm3/s。

根据达西定律，各小层的渗流阻力表达式为

式中：μa为油水相的视粘度，mPa·s。

式（8）和式（9）整理后得到各小层视粘度为

式中：μai为注采单元第i小层的视粘度，mPa·s。

由式（10）可以看出，当油水粘度一定时，各小层的视粘度是其含水饱和度的函数。为此，可先建立注采单元油水相对渗透率曲线，且各小层采用相同的相对渗透率曲线（图1），再根据油水相对渗透率曲线建立各小层在油、水粘度分别为40和0.45 mPa·s时的视粘度与含水饱和度的关系曲线（图2）。

图2 注采单元小层含水饱和度与视粘度的关系

当小层的视粘度一定时，其对应的平均含水饱和度是一定的，所以根据式（10）确定的各小层的视粘度，结合小层视粘度与含水饱和度关系曲线可以得到各小层的平均剩余油饱和度。

2 应用实例

在实际矿场应用过程中，分层平均剩余油饱和度无法直接通过检测手段获取，需要采用油藏数值模拟技术对所建方法进行验证。建立了一注四采五点法井组模型，模型规模为300m×300m。纵向上分为5个小层，小层厚度分别为1，2，3，4和5m；渗透率分别为100×10-3，200×10-3，500×10-3，800×10-3和1000×10-3μm2；孔隙度均为0.2。各小层之间不连通，没有油水串流。5个小层合采合注，油井采用定液量方式生产。当油藏含水率达到90%时，根据油藏模拟中各小层的吸水量，采用新建方法计算得到各小层的平均剩余油饱和度，并与油藏数值模拟计算得到的各小层的平均剩余油饱和度进行对比（表1）。从表1中可以看出，两者的相对误差在5%以内，表明新建方法计算结果准确，可以快速准确地对层间剩余油饱和度进行预测。

表1 各小层平均剩余油饱和度计算结果对比

3 结束语

分层吸水量的差异主要是由各小层渗流阻力的差异决定的，而各小层的渗流阻力受小层储层物性和剩余油饱和度的综合影响，在储层物性一定的情况下渗流阻力随剩余油饱和度的变化而变化。因此基于小层吸水量得到的渗流阻力，建立小层平均剩余油饱和度的预测方法是可行的。吸水剖面测试是评价油藏层间采出程度差异的重要方法。利用吸水剖面资料建立了层间平均剩余油饱和度的快速计算方法，经油藏数值模拟验证，方法准确可靠，为油田制定合理的开发措施提供了依据。

［1］ 马圣贤，呼舜兴，耿师江，等.多层水驱油实验在桥口油田开发中的应用［J］.西安石油学院学报：自然科学版，2003，18（6）：48-50.

［2］ 崔英琢，冯佩真，陈代秀，等.胡12块水驱油藏储层动用状况变化分析［J］.江汉石油学院学报，2000，22（4）：66-68.

［3］ 冯其红，隋园园，陈先超，等.单层剩余油饱和度定量预测研究［J］.西南石油大学学报：自然科学版，2011，33（6）：117-120.

［4］ 李淑霞，陈月明.利用井间示踪剂确定剩余油饱和度的方法［J］.石油勘探与开发，2001，28（2）：73-75.

［5］ 俞启泰.关于剩余油研究的探讨［J］.石油勘探与开发，1997，24 （2）：46-50.

［6］ 陈岑，胡望水，徐博，等.高集油田高6块阜宁组剩余油分布规律［J］.油气地质与采收率，2013，20（4）：88-90，93.

［7］ 张少波.特高含水期剩余油分布的油藏数值模拟研究［J］.长江大学学报：自然科学版，2010，7（3）：218-220.

［8］ 冯其红，王守磊，白军伟，等.层间非均质油藏提液效果数值模拟［J］.油气地质与采收率，2013，20（3）：49-52.

［9］ BuckleySE，LeverettMC.Mechanismoffluiddisplacementin sands［J］.TransAME，1942，14（46）：107-116.

［10］徐庆岩，杨正明，何英，等.特低渗透多层油藏水驱前缘研究［J］.油气地质与采收率，2013，20（2）：74-76.

［11］张建国，杜殿发，侯建，等.油气层渗流力学［M］.2版.东营：中国石油大学出版社，2010.

编辑王 星

TE319

：A

：1009-9603（2014）05-0098-03

2014-07-08。

张世明，男，高级工程师，博士，从事油田开发研究及油藏数值模拟工作。联系电话：（0546）8715085，E-mail:zhang⁃shm855.slyt@sinopec.com。

国家科技重大专项“整装油田特高含水期提高水驱采收率技术”（2011ZX05011-002）。