双重介质油藏注水井试井解释模型的建立及应用

双重介质油藏注水井试井解释模型的建立及应用姜永;别旭伟;刘洪洲;王迪;吴浩君【摘要】Based on water flooding characteristics of dual porosity reservoir and Buckley-Leverett equation,the injection pressure drop testing interpretation model is place transform and Stehfest numerical inversion are used to solve the model and a pressure drop typical plate of injection well is obtained with dimensionless pressure drop solution.Results show that oil-water two-phase region has large pressure response characteristics,oil-water viscosity difference causes double logarithmic curve upward.The larger the difference of oil-water viscosity,the more two-phase double logarithmic curve upturns.So the upward curve caused by reservoir change or impermeability boundary should be distinguished.The application in an injection well of JZ oilfield shows that the interpretation results agree with field test data.The proposed model can be used to calculate the position of water flood front,and to analyze the formation around the injection well and boundary condition,thus guiding water injection well evaluation and the design of injection scheme.%针对双重介质油藏注水开发特征,基于Buckley-Leverett饱和度分布方程建立了注水井压降试井解释模型,通过Laplace变换和Stehfest数值反演对模型进行了求解,利用无因次压力降落解获得了注水井的压力降落典型图版.分析图版可知,两相区对注水井压力响应特征较大,油水黏度差异会导致双对数压力导数曲线上翘,并且两者差异越大,两相区双对数压力导数曲线上翘幅度越大,因此在试井解释时应区别储层物性变化或不渗透边界造成的压力导数曲线上翘.选取JZ油田1口注水井进行了实例应用,结果表明本文解释图版拟合结果与现场测试数据吻合较好,并且解释所得的储层参数合理.本文建立的试井解释模型可以计算注水前缘的位置,分析注水井周围地层信息及边界情况,对评价双重介质油藏注水井的动态及注水方案的设计具有重要的指导意义.【期刊名称】《中国海上油气》【年(卷),期】2017(029)004【总页数】6页(P98-103)【关键词】双重介质油藏;注水井;试井;解释模型;压力降落典型图版【作者】姜永;别旭伟;刘洪洲;王迪;吴浩君【作者单位】中海石油(中国)有限公司天津分公司天津300459;中海石油(中国)有限公司天津分公司天津300459;中海石油(中国)有限公司天津分公司天津300459;中海石油(中国)有限公司天津分公司天津300459;中海石油(中国)有限公司天津分公司天津300459【正文语种】中文【中图分类】TE344姜永,别旭伟,刘洪洲，等.双重介质油藏注水井试井解释模型的建立及应用[J].中国海上油气，2017,29(4)：98-103.JIANG Yong,BIE Xuwei,LIU Hongzhou,et al.Establishment and application of a well test interpretation model for injection wells in dual porosity reservoir [J].China Offshore Oil and Gas,2017,29(4):98-103.注水井试井可以获取油相和水相渗透率、注水井周围污染情况、注水前缘距离、油藏边界等信息，为注水井措施提供依据。

关于注水井试井解释方法，国内外许多学者开展了相关研究[1-4]。

于九政等[5]针对低渗透油藏注水开发特点，在考虑启动压力梯度方向性和时效性的情况下，提出了低渗储层注水井井下关井试井压力响应机理，建立了低渗储层井下关井条件下注水井的试井解释模型。

刘佳洁等[6]针对注水井试井未考虑含水率、水驱前缘等参数的问题，根据质量守恒原理，结合注水井岩心归一化相渗曲线，在油水两相区运用油水两相流渗流理论，建立了三区流体复合注水井试井模型。

李乃华[7]建立了多层均质定压边界灌注试井数学模型，针对适于现场应用的折线型灌注函数，获得了层间及井底压降的实空间精确解析表达式。

张艳玉等[8]假设测试时间内地层压力均匀分布，油藏为两区复合储层模型，不考虑每个区中饱和度梯度，建立了单一介质注水井压降试井解释模型。

但目前有关双重介质油藏注水井压降试井解释方法的研究较少。

为了更好地掌握双重介质地层渗流动态信息，为注水井的增注措施提供依据，本文考虑地层中饱和度梯度和注水阶段对压力降落试井的影响，建立了双重介质油藏注水井压降试井解释模型，并进行了实例应用，结果表明本文解释图版拟合结果与现场测试数据吻合较好，解释所得的储层参数合理可靠。

1.1 注水阶段的物理模型注水井物理模型见图1。

为研究方便，对模型作如下基本假设：1) 地层水平等厚，油藏由裂缝系统和基质系统组成；2) 注水井以恒定速率注水；3) 流体和岩石微可压缩，且压缩系数为常数；4) 流体在地层等温流动，渗流满足达西定律；5) 忽略重力和毛细管力影响；6) 考虑井储和表皮系数影响；7) 水驱油为非活塞式的，裂缝系统和基质系统饱和度分布满足Buckley-Leverett 方程。

1.2 注水井试井解释模型注水阶段的物理模型假设考虑了油藏内油水饱和度的分布，为方便求解，将油水两相区剖分为N个圆环(图2)，即将连续的饱和度分布用阶梯型分布的饱和度代替，当圆环个数趋近于无限大时，阶梯型饱和度分布趋于连续的饱和度分布(图3)。

剖分后每一圆环内的饱和度为一定值，即可认为每一圆环内的水驱油是活塞式的，其饱和度分布可由Buckley-Leverett方程确定，即每一圆环内的饱和度为将每个圆环内的Swi代入Buckley-Leverett方程，则可计算出每一圆环无因次半径rDi。

对于每一圆环，其扩散方程为裂缝系统基质系统其中；初始条件内边界条件外边界条件1) 无限大外边界2) 定压外边界3) 封闭外边界衔接条件1.3 数学模型的求解定义拉普拉斯变换则数学模型可化为裂缝系统基质系统整理得其中内边界条件衔接条件外边界条件1) 定压外边界2) 封闭外边界3) 无限大外边界对于上述数学模型，其方程为一组零阶的虚宗量Bessel方程，通解为对于N个圆环，共有2N个未知数(Ai，Bi)，N-1个界面提供了2N-2个方程，另外加上内、外边界条件，共有2N个方程，可确定2N个未知数。

使用Gauss消元法[9-12]解出未知数Ai、Bi。

根据式(19)、(20)计算出Laplace空间下井底处无因次压力的解。

借助于Stehfest拉普拉斯数值反演变换算法，最终求得在实空间无因次井底压力pwD。

Stehfest拉普拉斯数值反演变换算法[13]为zi=(ln2/t)i通过求解数学模型可以获得无因次井底压力pwD与无因次时间tD的关系。

假设关井测试阶段水驱前缘不再发生变化，利用叠加原理得到关井阶段的压力降落解为利用无因次压力降落解可获得注水井的压力降落典型图版，如图4所示，计算参数为CD=50、S=5、λ=0.000 01、ω=0.02。

根据导数曲线特征可划分7个流动段,分别对应井筒储集效应阶段、表皮效应阶段、基质向裂缝窜流阶段、水区径向流阶段、油水两相区响应阶段、总系统径向流段和边界响应阶段。

从图4可以看出，在窜流阶段结束后出现第1条水平直线段，表征的是井周围水区的响应；之后由于油水黏度的差异，压力导数曲线呈逐渐上翘的特征，反映的是油水两相区的响应，并且油水黏度差异越大，两相区双对数曲线上翘幅度越大(图5)。

另外，水驱前缘距离越远，油水两相区时间越长，双对数曲线上到达前缘压力响应时间越晚(图6)，但不影响由于流体性质差异导致的曲线上翘幅度。

因此，在试井解释时应区别储层物性变化或不渗透边界造成的压力导数曲线上翘，在到达水驱前缘后导数曲线上出现一个凸起，接着导数曲线出现第2条水平直线段，反映油区的压力响应特征。

JZ油田一口注水井X井，在进行压力降落测试前平均日注水256 m3，注水2个月累注水1.54×104m3。

该注水井及油藏基本参数见表1，地层流体的相对渗透率曲线及地层分流率曲线见图7。

该注水井实测压力数据与本文典型图版拟合情况见图8，可以看出，两者吻合较好。

根据压力拟合可计算出残余油饱和度下的水相渗透率，根据时间拟合可计算出前缘的位置[8]。

该注水井最终解释参数值见表2，可以看出，储层物性较好，注水井存在一定程度污染，水驱波及距离达到183 m。

1) 基于Buckley-Leverett饱和度分布方程建立了双重介质油藏注水井压降试井解释模型，通过Laplace变换和Stehfest数值反演对模型进行了求解，得到了注水井的压力降落典型图版。

分析图版可知，两相区对注水井压力响应特征较大，油水黏度差异会导致双对数曲线上翘，并且油水黏度差异越大，两相区双对数曲线上翘幅度越大，因此在试井解释时应区别储层物性变化或不渗透边界造成的压力导数曲线上翘。

2) 利用本文建立的试井解释模型可以计算注水前缘的位置，分析注水井周围地层信息及边界情况。

实例应用表明，本文解释图版拟合结果与现场测试数据吻合较好，解释所得的储层参数合理可靠。

Sw—含水饱和度，%；rD(Sw)—Sw所对应的无因次距离；q—注水量，m3/d；tp—注水时间，h；h—油藏厚度，m；rw—井半径，m；(Sw)—分流曲线上Sw对应的斜率；Sor—残余油饱和度，%；SwF—水驱油前缘饱和度，%；pmD—基质系统压力,无因次,；pfD—裂缝系统压力,无因次,；tD—时间，无因次,Kf—裂缝渗透率，D；K—地层平均渗透率，D；λ—窜流系数；φf—裂缝孔隙度,%；φ—地层孔隙度,%；Ctf—裂缝系统压缩系数,1/MPa；Ct—地层综合压缩系数,1/MPa；μ—流体黏度,mPa·s；re—油藏半径，m；B—体积系数,m3/m3；t—时间，h；S—表皮系数,无因次；CD—无因次井筒储存系数；pi—原始地层压力,MPa；z—拉普拉斯变量；I0、K0—第1类、第2类零阶修正贝塞尔函数；Ai、Bi—待定常数；rDf—水驱油前缘无因次距离。