三维非均质油藏剩余油分布数值模拟预测方法.docxVIP

三维非均质油藏剩余油分布数值模拟预测方法 数值模拟法是预测油层剩余油的饱和分布最常用的方法之一。在数学模型的构建过程中，基本原则是物质稳定性、动态稳定性和西部定律的规则性。经过严格的数学推导,建立了描述油层内压力、温度、饱和度以及不同驱油条件下各种组分在油层中渗流规律的数学模型,对这种数学模型用离散数学方法(有限差分法、有限元法和边界元等)组成非线性方程组,然后经过一系列近似处理,将非线性方程组化成线性方程组,然后综合利用迭代法和直接求解方法,求得这种线性方程组的解,也就得到了油层内各种场的分布以及井点动态指标变化规律。 到目前为止,油藏数值模拟理论模型、求解方法和计算机软件已经达到了近似完美的程度。在理论模型方面,以多重介质组分模型为基础,可以将人们能够想象到的任何一项油层中的条件加入到模型中而得到更完美的模型,同时可根据某种油藏的特殊性,减掉所有在这种油藏中不曾出现的因素,而得到特殊用途的简化模型。在解法方面,全隐式解法和牛顿迭代法被认为是最稳定和最可靠的方法,用这种方法组合,可以消除任何解的不稳定、不收敛问题。但是,油藏数值模拟存在两个难题:①解方程组必须经过多次迭代直接求解或解大的线性方程组;②要进行反复的历史拟合,须耗费大量的人工和机时。在油田应用中,需要不断地更换容量更大、计算速度更快的计算机,而且要有专门的技术人员操作数值模拟软件,即使这样,也难以做到对油田所有开发区块完成数值模拟计算。近年来,国内外学者在非线性渗流、数模简化算法等方面取得了大量研究成果,但仍须求解复杂的渗流方程。 笔者提出一种基于注入孔隙体积倍数与采出量关系的剩余油饱和度分布预测方法,此方法避免了数值模拟中多次迭代和解巨大线性方程组的运算,可大大提高预测剩余油分布及动态指标的时间。 1 在非稳定流条件下，剩余油饱和与注入孔体积的体积关系为 1.1 函数与sw的关系 根据Buckley-Leverett方程,对于图1所示的一维单向流,在开始注入的t时刻,若累积注入量为Wi(t),则含水饱和度为Sw的平面所在的位置x可表示为 x-x0=Wi(t)?Aφ(Sw)(1)x?x0=Wi(t)?Aφ(Sw)(1) 式中 x0为注入井点起始坐标;?为孔隙度;A为油层横截面积;φ(Sw)是水的分流量fw关于含水饱和度Sw的导数。 在稳定流条件下,fw可表示为 fw=Κrwμw/(Κrwμw+Κroμo)(2)fw=Krwμw/(Krwμw+Kroμo)(2) 式中 μw和Krw分别为水相的粘度和相对渗透率;μo和Kro分别为油相的粘度和相对渗透率。则有 φ(Sw)=dfwdSwφ(Sw)=dfwdSw 在t时刻,图1中x1和x2位置处的Sw1和Sw2与累积注入量的关系分别由式(1)得 x1-x0=Wi(t)?Aφ(Sw1)(3)x2-x0=Wi(t)?Aφ(Sw2)(4)x1?x0=Wi(t)?Aφ(Sw1)(3)x2?x0=Wi(t)?Aφ(Sw2)(4) 由式(3)和式(4)得 x2-x1=Wi(t)?A[φ(Sw2)-φ(Sw1)](5)x2?x1=Wi(t)?A[φ(Sw2)?φ(Sw1)](5) 令 Ιpv=Wi(t)?A(x2-x1)Ipv=Wi(t)?A(x2?x1) 由于模型中流体假设为不可压缩的稳定流,所以累积注入量与累积产液量以及流过任何一个断面的累积液量均相等,则Ipv恰好是流入(x2-x1)单元的累积孔隙体积倍数,而从式(5)可得 φ(Sw2)=φ(Sw1)+1Ιpv(6)φ(Sw2)=φ(Sw1)+1Ipv(6) 式(6)表明:油层中某个单元的φ函数是上游单元φ函数加上流过本单元累积液量孔隙体积倍数的倒数。另外,从式(1)可得 φ(Sw)=?A(x-x0)Wi(t)φ(Sw)=?A(x?x0)Wi(t) 由此可见,φ(Sw)是从注入端算起到x处孔隙内累积注水的体积倍数。而在产出端,φ函数则是注入到这一地层中累积水量孔隙体积倍数的倒数,也就是累积产液孔隙体积倍数的倒数。 若图1中x1和x2之间网格中Sw的平均值为ˉSwSˉˉw,则有 φ(ˉSw)=φ(Sw1)+12Ιpvφ(Sˉˉw)=φ(Sw1)+12Ipv 在实际应用中,φ(Sw1)可作为上游网格的值。φ函数与Sw的关系一般如图2所示。 φ函数随Sw的变化规律是先升后降的过程,而用以上诸式表示的φ函数则始终是随着Ipv的增加而降低,也就对应着Sw不断升高。而在油层实际渗流过程中,只有油水前缘推进到的位置,含水饱和度才会不断上升;否则,无论流过去的流体量有多少,对应的含水饱和度总是束缚水饱和度,而此时计算出的φ函数的值也将大于前缘含水饱和度对应的值。据此可以得到利用φ函数求Sw时更严格的表达式为 Sw={φ-1(Sw)?φ(Sw)≤φ(Swf)Swc?φ(Sw)＞φ(Swf)(7)