裂縫性油藏稠油蒸汽注采數值模擬自適應網格法的研究.pdf

1、蒸汽注采是在裂縫性稠油油藏開發(fā)中應用非常廣泛且驅油效率較高的開采方法。對于裂縫性油藏稠油蒸汽注采的數值模擬，其困難在于在流場相變區(qū)域附近存在一個非常陡峭的溫度和飽和度鋒面。由于物理量在跨越鋒面界面時變化非常劇烈，在鋒面附近必須采用很小的網格尺寸才能滿足精度的要求。如果對全場都采用均勻細網格，計算量將非常巨大。因此，對于這類問題非常有必要采用能夠自動追蹤物理量變化劇烈的鋒面的自適應網格法，在物理量梯度大的位置采用細網格，而在物理量變化緩慢

2、的區(qū)域采用較粗的網格。 本文工作的目標是將自適應網格法應用于裂縫性油藏蒸汽注采的數值模擬。對于均勻裂縫性油藏，建議對裂縫和基質采用完全相同的自適應網格結構，這種結構下的自適應網格算法既簡便又有效。數值算例的結果表明該自適應網格法相比全精細網格計算能夠在保證計算精度的前提下大幅提高計算速度。同時，裂縫性油藏蒸汽輔助重力驅油(SAGD)過程算例的結果顯示了在SAGD過程中由于受到裂縫滲透能力強的影響，蒸汽腔容易較早地接觸油藏頂部，對

3、蒸汽熱能的利用率有著不利的影響。 然而，在非均勻裂縫性情況下實施自適應網格法，需要面臨一些復雜的數學物理問題。 首先，本文發(fā)現，由于裂縫與基質之間的流體交換受基質滲透率大小的控制，而且相鄰基質塊之間的聯(lián)系很微弱，因此基質滲透率在空間的急劇變化會導致相鄰基質飽和度的空間不連續(xù)，直接對基質粗網格的變量進行插值或平均將帶來問題的復雜性。為了解決這一問題，本文對裂縫和基質方程分別使用不同的網格系統(tǒng)——僅對裂縫使用自適應網格結構，

4、而對基質使用全精細網格。由于在DP模型下不考慮基質之間的影響，本文在Newton—Raphson迭代過程中使用一種特殊的裂縫．基質方程解耦方法，將基質方程和變量完全消去，僅剩下自適應網格下的裂縫方程和變量需要參與大型稀疏矩陣的求解。此外，這一解耦方法也同樣適合于DK模型，只需將相鄰基質網格的變量暫時“顯式化”，即在進行Newton—Raphson迭代過程中取上一迭代步的值。數值算例的結果表明，本文提出的針對非均勻裂縫性油藏的自適應網格法

5、相比于全精細網格計算也擁有很大的計算速度優(yōu)勢，并且保持了較高的計算精度。數值計算效率的提高來源于對裂縫采用自適應網格結構及方程求解時采用了這一特殊的裂縫—基質變量解耦方法。 其次，在非均勻情況下對粗網格等效滲透率的計算一直是國際上關注的難題。特別是對于多層網格結構下的等效滲透率求解方法，問題變得更加復雜。這是由于大多數等效滲透率方法蘊含有統(tǒng)計屬性(需要粗網格包含有大量的精細網格)或者需要使用周期性邊界條件。對于多層網格結構中的較

6、細層次，它們包含的精細網格數量太少，可能不足以使這些等效滲透率算法有效。為了解決這一問題，本文建議在估算一個粗網格的等效滲透率時，不應只關注其自身所攜帶的有限信息，而應考慮其他區(qū)域對他的影響。于是，本文從大尺度平均理論出發(fā)，利用大尺度平均理論在全場求得“輔助變量”，它描述了壓力和速度的脈動與大尺度平均的壓力梯度之間的關系。利用“輔助變量”攜帶的信息，可以求解任意一個局部區(qū)域的等效滲透率。這一方法提供了一個恰當的方式來計算多層網格結構下粗

7、網格的等效滲透率，而且無需對它引入任何邊界條件。數值算例的結果表明這種方法相對其他常用的方法對等效滲透率的估算在多層網格下有著更好的近似。 自適應網格法的計算速度很快，可以使人們在不同參數情況下快速進行大量的數值模擬，有利于發(fā)現油藏開采中的問題及優(yōu)化開采策略。在對一維裂縫性油藏蒸汽驅過程進行大量數值實驗后，發(fā)現毛管力滲吸驅油機制只在某些條件下起作用。當改變油藏參數導致進入油藏的液態(tài)水量減少時，裂縫中可能出現過熱蒸汽。在這些情況下