长岩芯驱替中不同位置的相渗曲线计算*

2014-06-01孙雷闫成海潘毅陈亮李士伦

西南石油大学学报(自然科学版) 2014年2期

关键词：调和油水渗流

孙雷，闫成海，潘毅，陈亮，李士伦

1.“油气藏地质及开发工程”国家重点实验室·西南石油大学，四川成都610500 2.中国石油新疆油田分公司陆梁油田作业区，新疆克拉玛依834000

长岩芯驱替中不同位置的相渗曲线计算*

孙雷1，闫成海1，潘毅1，陈亮2，李士伦1

1.“油气藏地质及开发工程”国家重点实验室·西南石油大学，四川成都610500 2.中国石油新疆油田分公司陆梁油田作业区，新疆克拉玛依834000

相对渗透率是一个解释多相不可压缩流体在多孔介质中流动的复杂影响的重要参数，以便校正达西单相流动方程。但是影响相对渗透率的因素很多，如：流体饱和度、岩石物性、润湿性、历史饱和度（滞后效应）上覆压力、黏土及颗粒含量、黏度、温度、表面张力、束缚水饱和度、驱替速度和末端效应等。利用桓冠仁、沈平平的均质岩芯一维渗流相渗推导结果，考虑长岩芯驱替室内实验中一维岩芯排列的非均质性和长岩芯逐级调和排序的特点，利用Matlab软件拟合出逐级调和平均渗透率与岩芯累计长度的关系，并反求出绝对渗透率倒数与累计岩芯长度的关系，将关系式由离散模型转化成连续模型，并代入一维均质相渗模型，利用WC油田长岩芯驱替实验数据计算出长岩芯驱替过程中不同位置的相对渗透率曲线，计算结果是可靠的。

相渗曲线；非均质；长岩芯驱替；拟合

孙雷，闫成海，潘毅，等．长岩芯驱替中不同位置的相渗曲线计算[J]．西南石油大学学报：自然科学版，2014，36（2）：139–144．

Sun Lei,Yan Chenghai,Pan Yi,et al．Calculation of Relative Permeability Curve at Different Position in Long Core Displacement Experiment[J]．Journal of Southwest Petroleum University：Science&Technology Edition，2014，36（2）：139–144．

在考虑流体饱和度、岩石物性、润湿性、驱替速度和末端效应等对相对渗透率影响的基础上[1-4]，利用桓冠仁[5]等的推导结果，以更符合油藏实际的实验室非稳态相渗求解方法确定相渗曲线。在非稳态法的求解中较早的有JBN法，但是由于以差分代替求导等因素的影响，相渗曲线形态并不规整。此外还有历史拟合法求解相渗曲线，其中有指数拟合法、幂函数拟合法、三次样条拟合法和B样条拟合法等，在特定范围内适应性较好。

Krause M H[1]考虑了毛管非均质性（微尺度上的非均质性）对相对渗透率的影响，结果认为即便是低非均质性的岩石，精细尺度上非均质性的影响依然相当明显。Li Heng等[6]用EnKF法以幂函数为相对渗透率函数历史拟合出了相对渗透率曲线，由于不需要对目标函数求导而很容易应用。Sheng Q等[2]用CT成像法、核磁共振法等数值求解相渗曲线。此外还有两相渗流的非线性模型[7-11]，应用逾渗理论求解相渗曲线，数值试井等，给相渗曲线的求解提供了新的方向。

本文用历史拟合法在岩芯离散调和排序的基础上，拟合逐级调和平均渗透率与累计长度的关系，并通过一定的数学方法使离散渗透率在拟合程度较好的情况下转换为随长度变化的连续渗透率，并求解出相对渗透率曲线，求解简单，结果在一定程度上是可靠的。

1 模型推导

（1）假定流体在模型中服从达西渗流[5，12-13]

式中：vo，vw—油、水相线性流速，cm/s；K—绝对渗透率，mD；Kro，Krw—油、水相相对渗透率，无因次；µo，µw—油、水相黏度，mPa·s；p—油水体系压力，MPa；x—流动空间维度。

（2）假定流体在模型中服从Buckley-Leverett一维渗流规律

式中：qo，qw—油、水相体积流量，cm3/s；ϕ—孔隙度，%；A—流动截面积，cm2；So，Sw—油、水相流动饱和度，%；t—时间，s。

由于qi=A×vi(i=o,w)，qt=qw+qo，则

式中：qt—油水相体积流量之和，cm3/s；

K（x）—与长度有关的绝对渗透率，mD。

对p关于x在区间[0，L]上积分，得

式中：Δp—流体流动出入口压差，MPa；L—流体流动距离，cm。

式中：α—两相阻力系数，mD−1；Θ—两相总阻力与单相阻力的比值，cm−2；fw—含水率，%；—无因次流动长度，无因次。

代入连续性方程，并用特征函数法求解，得

式中：x(S,t)—等饱和度面的位移函数，cm；

Ψ(S)—含水率关于饱和度的导数；

Vc(t)—累计产液量，cm3。

其中，

将式（10）、式（11）代入式（6），并进行微分，得

从而由式（5）、式（7）可得

式中：Ψe—含水率关于饱和度的导数在出口端的值；Vw—累计产水量，cm3；Swc—束缚水饱和度，%；Se—出口端含水饱和度，%。

对于长岩芯驱替，如果要采用1 m左右的完整天然岩芯作驱替实验，从取芯技术来讲是不可行的[14]，因此，国内外普遍采用常规短岩芯按一定的排列方式拼成长岩芯。为了消除岩石的末端效应，每块短岩芯之间用滤纸连接。经加拿大Hycal公司的Tomas等人论证，当岩芯足够长（1 m左右）通过在每块小岩芯之间加滤纸可将末端效应降低到一定程度。每块岩芯的排列顺序按下列调和平均方式排列。由式（18）调和平均法算出值，然后将值与所有岩芯的渗透率作比较，取渗透率与最接近的那块岩芯放在出口站第一位；然后再由剩余岩芯求出K，将新求出的值与所有剩下的（n−1）岩芯作比较，取渗透率与新的值最接近的那块岩芯放在出口端第二位；依次类推便可得出岩芯排列顺序。

利用上述岩芯排序计算方法，得到所选岩芯从出口端到入口端的排序结果（表1）。

表1 岩芯逐级调和平均值与岩芯长度数据Tab.1 The data of harmonious average permeability and length of cores

图1和图2分别给出高渗组和低渗组长岩芯实验数据的拟合结果。表2给出了逐级调和平均值与岩芯长度之间关系的拟合结果。可以看出，两组拟合效果较好，由于低渗组拟合渗透率级差较大（＞100），除个别点外使得RMSE（误差均方根）偏大外，拟合效果较好，可以满足预测要求。

图1 低渗组拟合曲线与残差图Fig.1 The fitted curve and residuals of lower permeable cores

图2 高温高压岩芯滤失仪简易图Fig.1 Diagram of high temperature and pressure core filtration instrumentation

表2 逐级调和平均值与岩芯长度拟合结果Tab.2 The fitted results of harmonious average permeability and accumulated length of core_s

由于拟合得到的值是调和渗透率，需将调和平均值转换为单块岩芯的渗透率值，由调和平均式得

式中：

Kk—第k块岩芯的渗透率反演值，mD；

Kk+1—前k+1块岩芯的渗透率调和平均值，mD；

lk—第k块岩芯的长度，cm；

lj—第j块岩芯的长度，cm；

n—岩芯总块数；

k—从入口端到出口端岩芯序号。

将单个岩芯简化为质点，即lk趋于零时，并将回归所得公式代入式（19），可获得K的连续表达式为

2 相渗曲线计算

（1）作Δp，V0，Vw与t的关系曲线，把实验数据做一次光滑处理

（2）计算Ψe，Θ，Swa

式中：Q—总产液量，cm3/s；

Swa—平均含水饱和度，%。

（3）用公式逼近Θ～Ψe，Swa～Ψe曲线

拟合得[15]a1=0.5343,b1=0.04398,c1=−0.4338,d1=−0.5202

由上式可知

结合式（13）、式（21），可求得水驱油过程相渗曲线。

2.1 高渗组长岩芯不同位置相对渗透率曲线特征

图3给出高渗组长岩芯距注入端1.0，35.0，65.0 cm拟合的相对渗透率曲线。由图3可以看出：

（1）不同位置等渗点相对渗透率几乎相等，在0.760 6～0.763 9。不同位置的点几乎在相同的含水饱和度上两相渗滤能力达到最弱，这主要是由于岩芯同一润湿性、在统一的温度、压力系统下两相的渗流特性有一定的统一性。油水相相对渗透率随距入口距离的增加而渗滤能力不断减小。

（2）Kro+Krw的值从开始的1.000 0下降为0.102 7。下降幅度较大，说明在驱替末段尽管驱替相地层水的流动能力在压差放大的情况下有所增大，但是油水两相的渗滤能力在不断下降，有将近50%的下降幅度，地层的两相渗滤条件在减弱。这主要是由于岩芯胶结疏松，且泥质含量较高，流体冲刷造成的颗粒运移堵塞渗滤空间所造成的。

（3）驱替后段，驱替相地层水流动能力增加，形成一定的冲洗效果，再加上两相渗滤能力的减弱，有更多的剩余油被驱替出来，虽然油量较少但下降较缓，驱替效果在一定程度上是可观的。

图3 高渗组岩芯1.0，35.0，65.0 cm处相对渗透率曲线Fig.3 The relative permeability curve of more permeable cores in 1，35，65 cm

2.2 低渗组长岩芯不同位置相对渗透率曲线特征

图4给出低渗组长岩芯距注入端10.8，45.0，83.0 cm处拟合的相对渗透率曲线。由图4可以看出：

（1）不同位置等渗点相对渗透率几乎相等，但其值在0.811 4～0.816 4（大于高渗组长岩芯）。不同位置的点几乎在相同的含水饱和度上两相渗滤能力达到最弱，这主要是由于岩芯同一润湿性、在统一的温度、压力系统下两相的渗流特性有一定的统一性。

（2）同样，油水相相对渗透率随距入口距离的增加而渗滤能力不断减小。这也是由于从入口到出口的压力梯度分布所造成的，距入口越远的点，其与出口间的压差就越小，从出口端冲刷过来的泥质颗粒便越会在出口附近被阻滞、沉降下来堵塞孔道，引起两相渗滤能力下降。

图4 低渗组岩芯10.8，45.0，83.0 cm处相对渗透率曲线Fig.4 The relative permeability curve of low permeable cores in 10.8，45.0，83.0 cm

3 结论

（1）在约翰逊公式的基础上，根据桓冠仁、沈平平计算公式，考虑了实验室长岩芯驱替中，岩芯渗透率随岩芯长度的变化，并拟合出岩芯渗透率随调和平均渗透率的变化趋势。

（2）利用渗透率与岩芯长度的变化关系，可求解出长岩芯驱替过程中，长岩芯不同位置的相对渗透率变化曲线。

（3）拟合结果显示，在调和渗透率与岩芯长度拟合差的点，相渗曲线可能会产生一定的变形，但是对调和平均值与累计长度拟合好的点而言，相渗规律可靠。所得规律能够用于分析长岩芯驱替对渗流特征的影响。

[1]Krause M H.Modeling and investigation of the influence of capillary heterogeneity on relative permeability[C].SPE 160909，2012.

[2]Sheng Q，Thompson K E，Fredrich J T.Numerical prediction of relative permeability from micro CT images：comparison of steady-state versus displacement methods[C]. SPE 147431，2011.

[3]Bennion D B，Thomas F B.Recent improvements in experimental and analytical techniques for the determination of relative permeability data from unsteady state flow experiments[C]//10thTechncialConferenceandExposition，Port of Spain，Trinidad，SPE 1991.

[4]赵春鹏，岳湘安.特低渗透油藏超前注水长岩芯实验研究[J].西南石油大学学报：自然科学版，2011，33（3）：105–108.Zhao Chunpeng，Yue Xiang’an.Long-core experimental study on advance water injection of extralow permeability reservoir[J].Journal of Southwest Petroleum University：Science&Technology Edition，2011，33（3）：105–108.

[5]桓冠仁，沈平平.一种非稳态油水相对渗透率曲线计算方法[J].石油勘探与开发，1982（2）：52–58. Huan Guanren，Shen Pingping.Calculating relative permeability[J].Petroleum Exploration and Development，1982（2）：52–58.

[6]LiHeng，ChenShengnan，YangDaoyong，etal. Ensemble-based relative permeability estimation using B-spline model[J].Transport in Porous Media，2010，85（3）：703–721.

[7]邓英尔，刘慈群.两相流体非线性渗流模型及其应用[J].应用数学和力学，2003，24（10）：1049–1056. Deng Ying’er，Liu Ciqun.Mathematical model of twophasefluidnonlinearflowinlow-permeabilityporousmedia with applications[J].Applied Mathematics and Mechanics，2003，24（10）：1049–1056.

[8]王坤，张烈辉.考虑应力敏感超低渗油藏油水相对渗透率的计算[J].石油天然气学报，2011，33（11）：117–119. Wang Kun，Zhang Liehui.Calculation of relative permeability of extra-low permeability reservoirs considering stress sensitivity[J].Journal of Oil and Gas Technology，2011，33（11）：117–119.

[9]王金勋，刘庆杰，杨普华，等.应用逾渗理论计算非稳态法油水相渗曲线[J].石油勘探与开发，2001，28（2）：79–82. Wang Jinxun，Liu Qingjie，Yang Puhua，et al.Calculating the relative permeability from displacement experiments basedonpercolationtheory[J].PetorleumExplorationand Development，2001，28（2）：79–82.

[10]向祖平，张烈辉，陈辉，等.相渗曲线对油水两相渗流数值试井曲线的影响[J].西南石油大学学报，2007，29（4）：74–78. Xiang Zuping，Zhang Liehui，Chen Hui，et al.Influnce of the relative permeability curve on numerical well test curve of oil-water-phase flow[J].Journal of Southwest Petroleum Univesity，2007，29（4）：74–78.

[11]杨小平.精确计算相对渗透率的方法[J].石油勘探与开发，1998，25（6）：63–66. Yang Xiaoping.A method to calculate the relative permeability accurately[J].Petorleum Exploration and Development，1998，25（6）：63–66.

[12]Johnson E F，Bossler D P，Bossler V O.Calculation of relative permeability from displacement experiments[C]. SPE 1023，1959.

[13]王学武，杨正明，时宇，等.核磁共振研究低渗透砂岩油水两相渗流规律[J].科技导报，2009，27（15）：56–58. Wang Xuewu，Yang Zhengming，Shi Yu，et al.Experimental study of water-oil two-phase fluid flow in low permeabilityreservoirbynuclearmagneticresonance[J].Science&Technology Review，2009，27（15）：56–58.

[14]李士伦，张正卿，冉新权.注气提高石油采收率技术[M].成都：四川科学技术出版社，2001.

[15]张志涌.精通Matlab R2011a[M].北京：北京航空航天大学出版社，2011.

编辑：牛静静

编辑部网址：http：//zk.swpuxb.com

Calculation of Relative Permeability Curve at Different Position in Long Core Displacement Experiment

Sun Lei1,Yan Chenghai1,Pan Yi1,Chen Liang2,Li Shilun1
1.State of Key Laboratory of Oil&Gas Reservoir Geology and Exploitation，Southwest Petroleum University，Chengdu，Sichuan 610500，China 2.The Luliang Operation Section of Xinjiang Oil Field，Petrol China，Karamay，Xinjiang 834000，China

Relative permeability is an empirical parameter，which interprets the complex effects of the flow of the multi-phase and incompressible fluid in porous media，so as to calibrate Darcy’s single phase flow equation.The relative permeability can be influenced by many factors，such as fluid saturation，physical rock properties，wettability，saturation history，overburden stress，clay and fines content，temperature，interfacial tension，viscosity，magnitude of initial phase saturations，and displacement rates and capillary outlet phenomena.Based on the results of Huan Guanren and Shen Pingping’s derivation of relative permeability curve and by considering the character in the displacement experiment whose core arrangement is heterogeneous in one dimension and harmonious arrangement，we got the relationship between cumulative core length and harmonious average permeability by Matlab.And then we reversely solved the relationship between absolute permeability and cumulative core length.We also converted the discrete model into the continuous model.By substituting the above results into the single dimension homogenous model，we obtained the relative permeability curve in different position of the displaced cores in WC oilfield.The result is reliable except for few points that fitted poorly.

relative permeability curve；heterogeneity；long core displacement；fit

http：//www.cnki.net/kcms/doi/10.11885/j.issn.1674-5086.2013.03.20.01.html