特高含水期油水两相渗流特征表征方程.docxVIP

特高含水期油水两相渗流特征表征方程 油田开发存在许多问题：高水位、高开采速度、严重的储水损失和高损失速度；注采孔网不完整，水驱动量低。层之间的不均匀性和小差异。油藏综合含量高，采收过程高，剩余油分布分散。注入水效率低下，浪费循环严重，新储水减少，产水减少，挖井的难度增加。如何提高油田最终采收率是油田开发面临的主要问题,也是采油工程面临的最大挑战。以往关于相渗曲线的研究与应用大多停留在油水相对渗透率比与含水饱和度指数式呈线性关系基础上。 进入特高含水开发阶段,油水相对渗透率比与饱和度指数式不再呈线性关系,基于该线性规律得到的分流量曲线无法准确描述特高含水期油水两相渗流特征,因而需要建立新的相渗关系表征方程,指导特高含水期油藏开发实践。 1 特高含水期的处理工艺 以某油田中14-检10井所在区块水驱油试验为基础,研究特高含水期相渗关系,试验结果见图1、2。 由图2可知,进入特高含水期以后,油水相对渗透率比Kro/Krw与含水饱和度Sw半对数关系曲线发生转折,出现新的直线段。因此,以往的相渗关系公式不能再应用。需要对以往相对渗透率比Kro/Krw与含水饱和度Sw的回归公式Kro/Krw=ae-bSw进行修正,进而指导生产实践。 2 高含水期相渗关系回归分析 2.1 回归分析的应用 以适用于中低含水阶段相渗关系表达式Kro/Krw=ae-bSw为基础, 利用数学分析构造函数的方法,确定特高含水期相渗关系的表达式为: ln(Kro/Krw)=a+bS2w+cln(1?Sw)a,b,c≥0(1)ln(Κro/Κrw)=a+bSw2+cln(1-Sw)a,b,c≥0(1) 式中,a为束缚水饱和度对应的半对数图版下的坐标点与高含水期过度到特高含水期对应饱和度在半对数图版下的坐标点所形成的直线在Y轴上的截距;b为该直线的斜率;c为残余油饱和度对应的坐标点与特高含水期对应饱和度在半对数图版下的坐标点所成直线的斜率。 特高含水阶段, 含水率fw≥90%后,原相渗曲线出现下弯段,且取值趋近于无穷大。即方程(1)左边Kro/Krw→0,ln(Kro/Krw)→-∝;方程(1)右边Sw→1时,a+bS2w+cln(1-Sw)→-∝,因而表征方程满足该性质。 以上述油田水驱油试验数据为基础,利用回归分析数学方法,建立适合于该油田的特高含水期相渗关系表达式。对式(1)引入回归分析,确定误差平方最小:记ln(Kro/Krw)=y,Sw=x,则对应x、y的数据为(xi,yi),求取二乘误差E(a,b,c),即: E(a,b,c)=∑i=1k[yi?a?bS2w?cln(1?Sw)]2E(a,b,c)=∑i=1k[yi-a-bSw2-cln(1-Sw)]2 则最小二乘误差为: minE(a,b,c)=∑i=1k[yi?a?bS2w?cln(1?Sw)]2minE(a,b,c)=∑i=1k[yi-a-bSw2-cln(1-Sw)]2 此处给定数据对k=12,得: ?E(a,b,c)?a=?E(a,b,c)?b=?E(a,b,c)?c=0?E(a,b,c)?a=?E(a,b,c)?b=?E(a,b,c)?c=0 利用Matlab软件求解该线性方程组,确定出α=(a,b,c)=(23.96,192.7226,103.9373),则该区块特高含水期相渗方程为:ln(Kro/Krw)=23.96+192.7226S2w+103.9373ln(1-Sw)。 图3是基于水驱油试验数据,利用公式Kro/Krw=ae-bSw与结合软件Matlab回归出来的折线图。 结合统计学方差分析理论,利用软件计算出利用公式Kro/Krw=ae-bSw的回归误差为47.3%;利用适合于特高含水期的公式ln(Kro/Krw)=a+bS2w+cln(1?Sw),a,b,c≥0ln(Κro/Κrw)=a+bSw2+cln(1-Sw),a,b,c≥0回归出误差仅为8.47%。 从对比图形及回归误差均可以看出,该特高含水区块相渗拟合曲线更贴合理论与矿场实际。 2.2 含水率与含水率的关系 结合Kro/Krw与Sw表征关系式,改进分流量曲线(fw与Sw关系曲线)。在忽略毛细管力、重力作用下的分流方程为:fw=11+μwKoμoKwfw=11+μwΚoμoΚw,于是确定含水率与含水饱和度的关系为: fw=[1+μwμoe23.96+192.7226S2w+103.9373ln(1?Sw)]?1=[1+μwμoe23.96+e192.7226S2w+(1?Sw)103.9373]?1(2)fw=[1+μwμoe23.96+192.7226Sw2+103.9373ln(1-Sw)]-1=[1+μwμoe23.96+e192.7226Sw2+(1-Sw)103.9373]-1(2) 利用方程(2)结合水驱油试验数