wiw

油藏自流注水开发机理及影响因素分析

苏海洋;穆龙新;韩海英;刘永革;李波

【摘要】针对自流注水开发现有理论还不够成熟的问题,研究了自流注水需要满足

的油藏条件,定义了“自流注水门限压力”的概念,并根据流体力学和油藏工程的基

本原理建立了其表达式,提出了自流注水过程中注水量、产油量、水层压力、油层

压力、累计注水量和累计产油量等参数的计算方法.根据提出的计算方法编制了相

应的计算程序,并运用程序进行了自流注水实例计算和影响因素分析,结果表明,提出

的自流注水计算方法原理简单,计算结果合理,能够用于自流注水参数计算;在水油储

量比小于60时,自流注水稳产期随水油储量比的增大而延长,采收率随水油储量比

的增大而增大,但在水油储量比超过60后,水油储量比的影响不明显;自流注水时机

对油田开发稳产期和采收率有重要影响,自流注水时机越晚,地层能量的利用越充分,

稳产期越长,采收率也越高.

【期刊名称】《石油勘探与开发》

【年(卷),期】2015(042)005

【总页数】6页(P632-637)

【关键词】自流注水;开发机理;门限压力;稳产期;采收率;注水时机

【作者】苏海洋;穆龙新;韩海英;刘永革;李波

【作者单位】中国石油勘探开发研究院;中国石油勘探开发研究院;中国石油勘探开

发研究院;中国石油大学(华东);中国石油勘探开发研究院

【正文语种】中文

【中图分类】TE341

自流注水是指高压水层的水在压差作用下通过套管自然流入低压油层以保持油藏压

力并驱替原油的过程[1]。目前自流注水多用于人工注水前地层能量的适当补充，

对于地面水源缺乏的地区如沙漠等也比较适用[2-3]。自流注水技术在国内应用较

少，仅在海上平湖油田有过应用先例[4]。国外从20世纪70年代起即开始这方面

的研究，其应用多见于中东地区[5-8]。

关于自流注水的研究多见于应用实践方面的报道，而对自流注水的开发机理缺乏系

统的研究。DaviesCA等[1]对自流注水注水量的计算方法进行了探讨，但未考虑

油井见水后的情况。本文对自流注水机理进行研究，运用油藏工程方法，建立一种

描述自流注水过程中油层注水量、产量、油层压力以及水层压力等参数的方法，并

对该方法的应用效果进行分析。

1.1自流注水门限压力

自流注水需要在水层和油层之间有足够的压差克服自流注水过程中水层到井筒的产

水压差、井筒内的摩擦损失、井筒内的水柱压力以及井筒到油层的注水

压差[9-10]。为了量化这一过程，参考文献[11]的研究成果，笔者定义了“自流注

水门限压力”，即通过自流注水能够使油层在设计产量下达到注采平衡所需的水层

压力。而油层在设计产量下达到注采平衡所需的注水量称为自流注水门限注水量。

根据压力关系，有如下关系式：

水层到井筒的产水压差为：

井筒到油层的注水压差为：

由达西-威斯巴哈公式[12]，井筒内的摩擦损失为：

其中，水力摩擦系数为：

井筒水柱压力为：

将（2）—（6）式带入（1）式可得：

由（7）式即可计算自流注水门限压力。只有水层压力大于自流注水门限压力，油

层才可以在满足产液量要求的条件下采用自流注水技术开发。

1.2自流注水参数计算

油层进行自流注水开发时，自流注水的注水量受水层压力和油层压力的变化所控制，

反过来注水量的变化又影响水层压力和油层压力的变化，进而影响油藏的稳产期和

采收率，因而注水量的求取很重要。假设在自流注水某一时刻t，油层注水量为

qiw，产油量为qo，产水量为qw，水层压力为pw，油层压力为po，由压力关

系可得：

同（2）—（7）式的推导过程，可得：

t时刻油层累计注水量、累计产油量、累计产水量分别为：

对于油层，由物质平衡方程可得[13]：

对于水层，由物质平衡方程可得：

油层在定液量条件下生产：

随着自流注水的进行，水层压力不断下降，水层与油层之间的压差随之不断下降，

由（9）式可知自流注水的注水量不断下降。当油层注水量低于产液量时，油层开

始亏空，油层压力由上升转为下降，当油层压力下降到低于维持油层定液产量的临

界压力时，油层转为定井底流压生产，此时油层产液量可以表示为：

临界压力的计算方法为：

油井见水前，有qw=0，Wp=0，qL=qo，油层产液指数为常数。由（9）—

（13）、（15）或（16）6个方程可依次求解得到qiw、Wiw、qo、Np、po、

pw6个未知量，具体计算程序见图1。

见水时，由Buckley-Leverett水驱油理论可得[14]：

则见水时刻有

由油水相渗曲线数据可求得fw-Sw关系与fw′-Sw关系，进而可求出fw′（Swf）

值，再由（19）式可求得见水时刻的累计注水量Wiw(T)，从而可以判断见水时间。

见水后，同理可得：

由油水相渗曲线数据可求得fw-Sw关系与fw′-Sw关系，然后由（21）式求得的

fw′（Swo）可进一步求得油井出口端含水率fw（Swo），则对应见水后某一时刻

的产水量与产油量为：

见水后，油层产液指数JL随含水率变化而变化。由文献[15]可知，见水后某一时

刻油井的产液指数与含水率为0时的产液指数比值为：

由（24）式可求得见水后不同含水率下产液指数的变化。

油井见水后，由（9）式—（15）式（或（16）式）、（21）式—（24）式共11

个方程以及相渗曲线数据可依次求得qiw、Wiw、JL、qL、fw、qo、Np、qw、

Wp、pw、po11个未知量，具体计算程序见图1。

由自流注水参数的求解过程可以看出，自流注水与人工注水的区别体现在（9）式，

即自流注水的注水量由水层压力、油层压力的变化控制。

（9）式是关于qiw的非线性方程，可用牛顿迭代法进行求解[16]。令F（qiw）

为qiw的函数：

给F（qiw）=0设定一个初值解qiw(n=0)，则有：

设定误差限为ε，若|qiw(n+1)−qiw(n)|≤ε，则qiw(n+1)是方程F(qiw）=0的解，

即（9）式的解，于是得到任意时刻注水量大小。

根据以上原理，编制了自流注水计算程序，可计算自流注水过程中任意时刻的注水

量、产油量、产水量、水层压力、油层压力、累计注水量、累计产油量、累计产水

量等参数。

2.1实例计算

中东N油藏地处沙漠地区，地面水资源缺乏。在N油藏以上1478.9m处的M

层发育有储量丰富的水层（见图2）。N油藏岩性为石英砂岩，非均质性较弱；孔

隙度18%～21%，渗透率（500～700）×10−3μm2；原始油藏压力约39.34

MPa，饱和压力约19.24MPa，储集层厚度15m。水层储量约为油层原油储量

的15倍；孔隙度约22.4%，渗透率约600×10−3μm2；原始水层压力25.48

MPa。油层吸水指数116.74m3/（d·MPa），产液指数104.87m3/（d·MPa），

水层产水指数137.16m3/（d·MPa）。原始油藏压力下原油体积系数为1.49，地

层水黏度0.68mPa·s，原油黏度约为0.90mPa·s，地层原油压缩系数

24.5×10−4MPa−1，地层水压缩系数5.66×10−4MPa−1，岩石压缩系数

8.73×10−4MPa−1。注采井网为排状注水，井距500m，排距500m，井筒直

径0.1778m。油藏开发要求最低压力高于饱和压力，因此生产井井底最低流压

pwf为20MPa。

由油藏条件可知，油层和水层高孔高渗，非均质性较弱，水层产水能力及油层吸水

能力较强，且原油黏度较低，适合进行自流注水开发。将N油藏与M水层连通，

当水层与油层之间有足够的压差时，则M水层中的水通过井筒流入N油藏进行能

量补充（见图3）。

2.1.1自流注水门限压力计算

N油藏原始油层压力为39.34MPa，远高于饱和压力19.24MPa和原始水层压力

25.48MPa，因此首先对油藏进行弹性驱动开发，使油层压力降低，水层和油层

之间建立较大的压差，然后进行自流注水开发。本算例中，当油层压力降低到不能

保持稳产的临界压力时开始转为自流注水开发，由（17）式可得临界压力为

23.79MPa。此压力即为开始自流注水时的油层压力poi。

计算时，稳产期内油层产量维持在397.5m3/d。由自流注水门限压力定义，开始

自流注水时，若要达到定产液量要求并维持注采平衡，注水量和产量应为592

m3/d，其他参数为：油层吸水指数116.74m3/（d·MPa），井筒与油层之间的

注水压差5.07MPa，水层产水指数137.16m3/（d·MPa），水层与井筒直接的

产水压差4.31MPa，水层到油层的高度差1478.9m，井筒内的摩擦损失为0.01

MPa，井筒水柱压力14.48MPa。将上述参数代入（7）式求自流注水门限压力。

N油藏求得的自流注水门限压力18.69MPa，远远低于M水层压力25.48MPa，

因此N油藏可以进行自流注水开发。

2.1.2自流注水参数计算

利用编制的自流注水计算程序，对N油藏进行自流注水参数计算，评价N油藏进

行自流注水的开发效果。

计算得到5年内注水量、产油量、油层压力、水层压力等参数随时间的变化，并

与Eclip数值模拟模型的计算结果进行对比（见图4—图5）。Eclip模型中

对于井筒内的摩擦损失采用多段井模型进行计算：

（28）式与本文计算摩擦损失的（4）式略有不同，因自流注水过程中井筒内摩擦

损失非常小，所以对计算结果的影响可忽略不计。由计算结果对比可知，本文的计

算模型与Eclip计算结果相近，表明本文建立的自流注水计算方法可靠。

N油藏弹性驱动待油藏压力降至23.79MPa时，开始实施自流注水。初始阶段油

层与水层压差较大，自流注水量较大，因而油层压力先小幅上升，然后由于注水量

不断降低，产液量保持不变，油层逐渐产生亏空，油层压力缓慢衰减。油层压力继

续保持在原始油层压力60%以上2.92a，水层压力始终不断下降（见图4）。

自流注水刚开始时注水量最大，为1005.77m3，然后由于油层压力与水层压力

之间的压力差逐渐减小，自流注水量逐渐递减，油层维持定产油量1.78a后开始

见水，见水后含水率迅速上升，产油量下降（见图5）。5a末自流注水累计产油

量438891m3，累计产水量167759m3。

如不实施自流注水，油层仅靠弹性驱动能量生产，则仅能维持稳产0.6a，自流注

水在此基础上可延长稳产期1.78a（见图6），提高采收率6.37%。

自流注水计算表明，从保持地层压力、维持稳产期及提高采收率角度考虑，N油

藏自流注水开发效果较好。

2.2自流注水影响因素

2.2.1水油储量比

本算例中，水层储量为油层储量的15倍，假设水油储量比分别为1、15、30、

45、60、75、90、105、120时，利用自流注水程序进行计算，研究不同水油储

量比对自流注水开发效果的影响（见图7）。

为考察不同水油储量比对自流注水稳产期的影响，将程序设置为定产油量生产。计

算结果表明，随水油储量比的增大，自流注水稳产期先是迅速增加，当水油储量比

达到60以上，水油储量比增大对稳产期影响不大，此时，自流注水稳产期大约

5.46a。

计算结果表明，随水油储量比增大，自流注水采收率先是迅速增加，同样，当水油

储量比达到60以上，水油储量比增大对自流注水采收率影响不大，此时，自流注

水采收率约为19.08%。

2.2.2自流注水时机

本例中，对N油藏先进行弹性驱动开发，当油藏压力降低到不能保持稳产的临界

压力（初始油藏压力的60%）时，开始转为自流注水开发。在水油储量比为15的

条件下，假设弹性驱动开发至油藏压力降到初始油藏压力的100%、90%、80%、

70%、60%时开始转为自流注水开发，对编制的自流注水程序略做改进，同时考

虑初期弹性驱动与后期转自流注水开发两个阶段，研究自流注水时机对开发效果的

影响（见图8）。

由图8可见，弹性驱动转自流注水时机越晚，对地层弹性能量的利用越充分，弹

性驱动+自流注水稳产期也就越长，弹性驱动+自流注水采收率也越高。因此在满

足其他开发要求的前提下，应尽量推迟弹性驱动转自流注水的时机。但如果地层压

力下降太多，溶解气析出会对开发产生不利影响，因此，转自流注水时机最晚应该

在地层压力降到饱和压力之前。

自流注水要求水层和油层之间有足够的压力差，用于克服自流注水过程中水层与井

筒之间的产水压差、井筒内的摩擦损失、井筒内的水柱压差以及井筒与油层之间的

注水压差，即水层压力应该高于自流注水门限压力。

根据流体力学原理和油藏工程基本原理，建立了自流注水门限压力的计算表达式；

并给出了自流注水过程中注水量、产油量、产水量、水层压力、油层压力、累计注

水量、累计产油量等参数的计算方法，通过VB语言编制了自流注水计算程序。实

例计算表明，该方法原理简单，计算结果合理可靠，能够用于自流注水参数计算。

运用自流注水计算程序进行了自流注水影响因素分析。结果表明，在水油储量比小

于60时，自流注水稳产期随水油储量比的增大而延长，采收率随水油储量比的增

大而增大，但水油储量比超过60以后，水油

储量比的影响不明显；自流注水时机对油田开发稳产期和采收率有重要影响，弹性

驱动转自流注水时机越晚，稳产期越长，采收率也越高。

符号注释：

A——流动区域截面积，m2；A0——井筒截面积，m2；a，b——单位换算系数，

a取2.29668×10−15，b取±1×10−6，水层在油层上部时b取负值，水层在油

层下部时b取正值；Bo——原油体积系数，m3/m3；Boi——原始原油体积系数，

m3/m3；Bw——水体积系数，m3/m3；Cto——油层综合压缩系数，MPa−1；

Ctw——水层综合压缩系数，MPa−1；d——井筒内径，m；f——Fanning系数；

fw——含水率，f；fw′——含水率的导数；g——重力加速度，m/s2；h——水

层到油层的高度差，m；Iw——油层吸水指数，m3/（d·MPa）；JL——油层产

液指数，m3/（d·MPa）；JL0——含水率为0时的产液指数，m3/（d·MPa）；

Jw——水层产水指数，m3/（d·MPa）；Kro——油相相对渗透率；Krw——水

相相对渗透率；L——注采井距，m；n——迭代步数；N——迭代截止步数；

No——油层储量，m3；Np——油层累计产油量，m3；Nw——水层储量，m3；

pc——油层临界压力，MPa；pf——井筒内的摩擦损失，MPa；ph——井筒水

柱压力，MPa；po——油层压力，MPa；pob——井筒与油层之间的注水压差，

MPa；poi——开始自流注水时油层的压力，MPa；pth——自流注水门限压力，

MPa；pw——水层压力，MPa；pwb——水层与井筒之间的产水压差，MPa；

pwf——最低井底流压，MPa；pwi——原始水层压力，MPa；qc——油层定液

产量，m3/d；qiw——自流注水注水量，m3/d；qL——油层产液量，m3/d；

qo——油层产油量，m3/d；qth——自流注水门限注水量，即水层门限产水量，

m3/d；qw——油层产水量，m3/d；Re——雷诺数，无因次；Sw——含水饱和

度，f；Swf——水驱前缘含水饱和度，f；Swi——原始含水饱和度，f；Swo——

采油井出口端含水饱和度，f；t——注水时间，d；T——见水时间，d；v——水

在井筒中的流动速度，m/s；Wiw——油层累计注水量，m3；Wiw(T)——油井

见水时累计注水量，m3；Wp——油层累计产水量，m3；λ——水力摩擦系数，

无因次；μo——原油黏度，mPa·s；μw——水黏度，mPa·s；ρw——地层水的

密度，kg/m3；φ——油层孔隙度，%。

【相关文献】

[1]oryandpracticeofmonitoringandcontrollingdumpfloods[R].SPE

3733,1972.

[2]ShizawiW,SubhiH,RashidiA,ementofoilrecoverythrough“dump-flood”

waterinjectionconceptinsatellitefield[R].SPE142361,2011.

[3]remaintenancebyformationwaterdumpingfortheRatawilimestone

oilrervoir,offshoreKhafji[J].JournalofPetroleumTechnology,1982,34(4):738-754.

[4]周俊昌,罗勇,严维锋.国内第一口自流注水井钻井实践[J].中国海上油气,2011,23(1):

Junchang,LuoYong,ngpracticeofthedomesticfirstwaterdumping

well[J].ChinaOffshoreOilandGas,2011,23(1):43-45.

[5]Al-GamberAA,Al-TowailibAA,licationofstand-aloneinjection

systemsinremoteand/orhighlypopulatedareasreducesconstructioncosts[R].SPE

63168,2000.

[6]QuttainahR,airdumpfloodpilotproject,theapplicabilityof

dumpfloodtoenhancesweep&maintainrervoirpressure[R].SPE68721,2001.

[7]QuttainahR,airproductionplateauextension:Theapplicability

offullfielddumpfloodinjectiontomaintainrervoirpressureandextendproduction

plateau[R].SPE97624,2005.

[8]MamdouhMI,JamesWS,HeshamLS,nmentallyfriendlyandeconomic

waterfloodsystemforOctoberfieldatgulfofSuez,Egypt[R].SPE112311,2008.

[9]IkawaH,MercadoG,licationinacarbonateoffshoreoil

field,U.A.E.[R].SPE117918,2008.

[10]ChangM,CullenR,UtomoB,zingwaterfloodingconsideringdipinthe

Wafrafield[R].SPE125916,2010.

[11]王良善,朱光亚,刘雄志,等.哈法亚油田注水开发适应性研究[R].北京:中国石油勘探开发研究院中

东研究所,angshan,ZhuGuangya,LiuXiongzhi,bilitystudyofwater

injectioninHalfayaoilfield[R].Beijing:MiddleEastRearchInstitute,RIPED,2011.

[12]袁恩熙.工程流体力学[M].北京:石油工业出版社,2001:eringfluid

mechanics[M].Beijing:PetroleumIndustryPress,2001:124-166.

[13]姜汉桥,姚军,姜瑞忠.油藏工程原理与方法[M].东营:中国石油大学出版社,2006:

Hanqiao,YaoJun,oirengineeringprincipleand

method[M].Dongying:ChinaUniversityofPetroleumPress,2006:200-225.

[14]张建国,杜殿发,侯健,等.油气层渗流力学[M].东营:中国石油大学出版社,2009:

Jianguo,DuDianfa,HouJian,gasflowmechanicsinporous

medium[M].Dongying:ChinaUniversityofPetroleumPress,2009:214-221.

[15]赵静,刘义坤,赵泉.低渗透油藏采液采油指数计算方法及影响因素[J].新疆石油地

质,2007,28(5):ng,LiuYikun,ationandinfluencingfactors

ofliquidproductivityindexinlowpermeabilityrervoir[J].XinjiangPetroleum

Geology,2007,28(5):601-603.

[16]张晓丹.应用计算方法教程[M].北京:机械工业出版社,2008:d

computationalmethods[M].Beijing:ChinaMachinePress,2008:40-60.