油氣藏流線模擬方法與工程應用研究進展

張 楠，孫 健，任 龍，曹 杰，陳明強

(1.西安石油大學 石油工程學院，陜西 西安 710065；2.油氣藏地質及開發(fā)工程國家重點實驗室(成都理工大學)，四川 成都 610059)

引 言

流線、跡線和脈線是流體力學中常用的基本概念。流線是某一確定瞬時流場中的一組空間曲線，每條曲線上的每一點都與速度矢量相切；跡線是指特定的流體質點在不同時刻經過的路徑；脈線是相繼通過空間一個固定點的流體質點連接而成的曲線，也稱為煙線。在定常流中，這三者在空間的軌跡上重合。流線模擬與基于網(wǎng)格的數(shù)值模擬方法不同，其流動方程的求解是沿著每條一維流線的，而非基于三維離散網(wǎng)格。由于其在高效計算和流場識別方面的特色，流線模擬在空氣動力學、地下水科學與工程、環(huán)境工程、石油與天然氣工程等領域都有廣泛的應用。

作為油藏數(shù)值模擬的一種核心算法，流線能夠根據(jù)儲層地質結構、油氣藏和流體物性、井網(wǎng)布置、完井方式、開發(fā)方式，說明流體的來向、流向與去向，直觀反映瞬時流場和儲層動用情況，動態(tài)展現(xiàn)注采井組的受效關系。這些數(shù)據(jù)很難通過其他數(shù)值模擬方法得到，因此，研究流線模擬方法和技術對科學高效開發(fā)油氣藏具有重要意義。

在石油與天然氣工程領域，流線模擬運算主要包括3個步驟：首先，在瞬時流場中追蹤流線；其次，沿著流線求解流體運移方程；最后，應用流線模擬結果解決工程問題。自1988年Pollock[1]提出了有效的3D流線追蹤方法以來，流線模擬技術在這3個方面都取得了快速發(fā)展：一是研究不同流動介質之內和之間復雜的流體交換和運移機制，進一步提高流線追蹤方法的精確性和計算效率；二是克服傳統(tǒng)流線模擬模型的局限性，將其拓展至可壓縮流體、非牛頓流體等運移問題的求解；三是充分利用流線模擬在高效計算和流場識別方面的特色，將其與其他數(shù)學模型或方法相結合發(fā)展新技術，應用于不同的工程問題。本文主要總結近年來流線模擬研究在這3方面取得的最新進展，分析其中的科學意義和創(chuàng)新性。

1 流線追蹤方法

流線追蹤是指在給定油藏地質模型、井網(wǎng)類型、完井方式、開發(fā)方案、油藏和流體物性的情況下，追蹤流線的軌跡并計算飛行時間[2-4]。流線追蹤是流線模擬的第一步，流線追蹤中出現(xiàn)的錯誤會在模擬、應用中被進一步放大，導致結果謬誤，因此在流線模擬中占據(jù)著重要地位。流線追蹤是在離散網(wǎng)格壓力已知的條件下進行的。依據(jù)Pollock[1]提出的理論，流線追蹤的一般步驟為：依據(jù)網(wǎng)格壓力或網(wǎng)格界面流量求解流線在本地網(wǎng)格中的出入口位置和飛行時間，并定位此流線的上、下游網(wǎng)格，依次重復以上2個步驟直至注采井，得到完整的流線?？梢?，流線追蹤的關鍵在于離散網(wǎng)格的劃分和速度場的逼近，速度場逼近的精確度決定了流線追蹤的精確度。

常用的油藏數(shù)值模擬商業(yè)軟件中的流線模塊是基于傳統(tǒng)的Pollock追蹤方法，在應用于復雜的地質模型時，常會出現(xiàn)流線追蹤結果不穩(wěn)定、甚至出錯的現(xiàn)象[5]。這主要是因為Pollock流線追蹤方法對速度場逼近的簡化處理：它假設同一網(wǎng)格邊界上的法向速度都相同，速度矢量沿著自身方向呈線性變化[1]。但實際上，網(wǎng)格界面上各點的速度不同，尤其是當?shù)刭|模型復雜或流動介質改變時，傳統(tǒng)的Pollock流線追蹤方法將無法滿足運算精度上的要求[6-7]。為保證流線模擬能夠成功應用于各類油氣藏，國內外學者依據(jù)流動介質和流體交換、運移機制的不同，從多個方面著手改進流線追蹤方法：(1)提出新的多項式來逼近速度場；(2)利用網(wǎng)格細化及流量重構描述快速變化的速度場；(3)對流動介質界面及附近區(qū)域進行局部處理。

1.1 速度場的多項式逼近

網(wǎng)格內部速度場逼近的前提是已知網(wǎng)格壓力或網(wǎng)格界面流量，因此，可以依據(jù)有限元等求解壓力或流量的數(shù)值方法提出速度場的逼近多項式。為保證計算精度，速度場的逼近多項式必須滿足物質平衡方程。Matringe 等[8]延用混合有限元方法中的有限元空間(最低階的Raviart-Thomas空間RT0和一階的Brezzi-Douglas-Marini空間BDM1)推導出了不規(guī)則三角形和四邊形網(wǎng)格內的速度場逼近多項式及相應流函數(shù)的半解析解。Juanes和Matringe[9]將使用高階多項式逼近速度場的流線追蹤方法正式命名為高階流線追蹤方法，并證明了使用高階方法比低階方法得到的速度場和流線更加精確，計算效率更高，受網(wǎng)格效應影響更小。

雖然速度場逼近多項式的提出是基于某種特定的網(wǎng)格類型和數(shù)值方法，但其應用并不局限于此。例如，Matringe等[10]通過將BDM1速度逼近多項式與多點通量逼近[11]方法中的分流量聯(lián)系在一起，將BDM1流線追蹤方法成功應用于有限體積方法。

然而，以上研究所使用的速度逼近多項式都是基于某種數(shù)值方法提出的，并未從實際物理規(guī)律的角度考慮其提出的依據(jù)。考慮到實際油藏中壓力分布是連續(xù)的，壓力梯度是對流流動的主要驅動力，且速度場與壓力梯度成正比，Zhang等[12]在2018年依據(jù)壓力分布規(guī)律推導速度場逼近多項式，提出了Bilinear(低階)和Cubic(高階)方法。通過聯(lián)系速度場逼近多項式和網(wǎng)格內部壓力之間的關系，使這些流線追蹤方法適用于所有能夠求解網(wǎng)格內部壓力的數(shù)值方法。在各向同性儲層中，使用Analytical、Pollock、BDM1、Cubic方法解釋x方向速度場的結果如圖1所示。通過對比以上方法發(fā)現(xiàn)，Cubic方法在CPU耗時相同的條件下，得到的速度場結果與解析解(Analytical)最接近，誤差最小，計算效率最高。

圖1 各向同性儲層中，使用Analytical、Pollock、BDM1、Cubic方法解釋x方向速度場Fig.1 Interpretation results of x-direction velocity field using Analytical,Pollock,BDM1,and Cubic methods in heterogeneous reservoir

1.2 網(wǎng)格細化及流量重構

網(wǎng)格細化和速度場重建是另外一種常見的提高流線追蹤精度的方法，這類方法將原始網(wǎng)格細化分為多個子網(wǎng)格，將原網(wǎng)格界面流量分配至各子網(wǎng)格界面，在子網(wǎng)格內追蹤流線。

Zhang等[7]系統(tǒng)地論述了有限差分、有限體積、有限元、混合有限元方法的共性和特性，提出了多邊形網(wǎng)格中逼近速度場的方法：將多邊形分解成多個子四邊形或三角形，對子網(wǎng)格界面分配流量，再選用速度場逼近多項式在子網(wǎng)格中逼近速度場。Jimenez 等[6]提出通過網(wǎng)格細化和界面層流量重建在斷層附近追蹤流線。如圖2所示，整個過程可以簡單概括為：在斷層界面細化出若干子網(wǎng)格，將子網(wǎng)格轉化為正六面體；重構流量并保證網(wǎng)格界面流量連續(xù)和原網(wǎng)格物質守恒；在正六面體內追蹤流線，再將正六面體內的流線變換至實際空間中。該方法保證了流量在斷層界面的連續(xù)性，使得流線追蹤較為穩(wěn)定。2018年，Chen等[13]對Jimenez等[6]的流線追蹤方法進行了拓展，將嵌入式模型中的裂縫-基質界面看作是一個類似于斷層-基質的界面，使其適用于嵌入式離散裂縫模型。

1.3 流動介質界面及附近區(qū)域處理辦法

速度場在流動介質界面及附近區(qū)域內的變化較為特殊，需要針對這些變化規(guī)律對其進行特殊處理，以提高流線追蹤的精確度。例如，井筒內的流動是符合Navier-Stokes公式的管流，與油藏基質內符合Darcy公式的滲流規(guī)律不同，流體在從儲層流入井筒的過程中，流動介質和流動控制方程發(fā)生很大變化，壓降和速度場的變化相比油藏中部更加劇烈。當速度場變化規(guī)律符合某種解析解時，則可以不使用離散網(wǎng)格而直接求解流線。據(jù)此，Wang等[14]依據(jù)利用邊界元方法的基本解在近井周圍求得壓力和速度的解析解后，通過積分直接得到流線。此方法避免了局部網(wǎng)格細化，求解速度快且精度高。但由于其計算所依賴的邊界元方法假設直井裸眼完井，這種方法的應用具有一定局限性。Tang 等[15]考慮到近井周圍壓力分布呈指數(shù)變化形式，使用柱坐標在扇形網(wǎng)格中提出了針對近井周圍區(qū)域的流線追蹤方法。此方法確保了模擬區(qū)域內的壓力連續(xù)，提高了近井周圍流線追蹤的精確度。在近井儲層具有非均質性的情況下，此優(yōu)勢更為顯著。

圖2 復雜斷層流線追蹤方法[6]Fig.2 Streamline tracing method for complex fault systems[6]

此外，如果油藏中部的流動介質發(fā)生變化，流動規(guī)律有顯著差異，則同樣需要進行特殊處理。對于縫洞型油藏來說，流動分為兩種，多孔介質基質和微裂縫內的達西滲流，與大裂縫和溶洞內符合Navier-stokes方程的自由流動，屬于典型的滲流-自由流耦合流動[16]。Peng 等[17]使用了Darcy-stokes模型模擬縫洞油藏中流體的流動，通過Pollock方法追蹤達西滲流中的流線，用流體力學中廣義的流線定義求解自由流動區(qū)域內的流線，再將兩者結合起來得到完整的流線。研究發(fā)現(xiàn)流線方法相比有限體積方法的計算速度更高，收斂速度更快。這項研究對3D溶洞中的流線追蹤具有重要的啟發(fā)意義和開拓作用，并證明了流線模擬在這類油藏的模擬中具有優(yōu)越性。

1.4 發(fā)展方向

基于以上討論可知，近10年來流線追蹤的基本方法有了進一步的發(fā)展，從速度場逼近多項式、網(wǎng)格細化及流量重組、局部區(qū)域特殊處理三方面著手，流線追蹤方法已經從低階流線發(fā)展至高階流線；應用范圍從普通砂巖油藏逐步向斷層、天然裂縫油藏、縫洞型油藏進行拓展；計算精確度、穩(wěn)定度和速度不斷提高。但對于現(xiàn)代油氣藏工程來說，流線追蹤方法仍面臨著許多的問題，制約著流線模擬的可靠程度和適用范圍，需要進一步發(fā)展：

(1)復雜結構井近井周圍的流線追蹤。復雜結構井，如水平井、多分支井等是非常規(guī)油氣藏開發(fā)的重要技術手段。復雜結構井的完井方式優(yōu)化、智能完井中的井下節(jié)流閥操控等等都與油氣藏的采收率有密切的聯(lián)系[18-19]。如何通過復雜結構井近井周圍的流線追蹤，研究復雜結構井完井方式與油氣在油藏中的流動路徑之間的關系，是下一步的重點方向之一。

(2)多場多尺度耦合流動規(guī)律的流線追蹤。碳酸鹽巖油氣藏、頁巖/致密油氣藏的孔隙結構都具有明顯的多尺度性，存在納米和微米尺度孔隙喉道以及微米到米尺度的天然裂縫和人工裂縫；大規(guī)模壓裂、化學驅、熱力采油等工程措施使得油氣滲流與應力場、化學場、熱場存在著耦合作用特征[4]。流動模擬的多尺度升級與多場全耦合的大規(guī)模計算是未來現(xiàn)代滲流理論的發(fā)展方向，考慮多尺度升級和多場耦合作用的流線追蹤方法研究尚處于起步階段。

2 基于流線的流體運移問題求解方法

流線是3D流管的中軸線，通過對每條流線分配流量、求解飛行時間能夠表示相應流管的孔隙體積。當流體沿著流線流動時，沒有流體進入或離開流管，利用流線的飛行時間重構3D笛卡爾坐標中的物質平衡方程，即可沿著1D流線求解運移問題[20]。如此，將非均質網(wǎng)格中的3D運移問題轉化為一系列沿著流線的1D運移問題，大幅度提高了計算速度。

在使用飛行時間重構物質平衡方程時，通常假設油藏及流體不可壓縮，重力和毛管力可忽略。在此條件下，流線變化緩慢，可選用較大的時間步長更新流線，計算效率高。這些假設在一定程度上限制了流線模擬方法的應用，因此，傳統(tǒng)觀念認為流線模擬適用于對流主導流動，將其應用在大型非均質油藏，而非彈性驅動、毛管自吸現(xiàn)象、重力驅替等其他因素對采收率影響較大的流動中[2-3]。經過多年發(fā)展，如今的流線模擬方法已經全面突破了這些限制，使用全三維模型模擬可壓縮流體[21]、毛細管壓力[22-24]、組分模型[25-26]、熱水驅油[27]等問題，已經成為了油藏數(shù)值模擬和油藏管理的重要方法[28]。

2.1 可壓縮流體的流線模擬方法

傳統(tǒng)流線模擬的假設之一是流體和油藏的不可壓縮性。不可壓縮的流動系統(tǒng)意味著油藏中的源、匯都是注采井，流線的兩邊連接著一組注采井；每條流線攜帶的流量恒定，同一流線上各組分的運移速度與各相流體流速線性相關，求解流體運移問題較為簡單。如果流體可壓縮，這樣的情況將會改變。由于油藏壓力分布差異，流線可能在任意網(wǎng)格中起始或結束；流體性質與壓力變化相關，沿著流線的流量發(fā)生了改變。

Cheng等[21]參考不可壓縮流體物質平衡方程的形式提出了有效密度的概念，得出可壓縮流體的物質平衡方程，沿著流線量化描述流體體積隨壓力變化的規(guī)律。在追蹤流線時，沿著流線計算飛行時間和有效密度；在計算飽和度時，加入與有效密度相關的源、匯來描述流體體積隨壓力的變化；最后，通過定期更新壓力、重新分配飽和度、更新流線來處理流場變化的情況。此方法保留了流線模擬中沿著1D流線求解飽和度方程的優(yōu)勢，只需沿著流線加入有效密度便可計算可壓縮流體的運移問題。Beraldo 等[29]進一步發(fā)展了此方法，通過實時核算沿著流線的累積體積和流量變化規(guī)律，確保了物質平衡條件，提高了飽和度計算的精確度。

2.2 對流-擴散類問題的流線模擬方法

按照流體力學的定義，流線上的每一點都與流體的合速度相切。當重力、毛管力等因素不可忽略時，各相流體的速度并不平行，組分不再完全沿著流線運移。此時需要通過算子分裂算法(operator splitting)[22,30-31]求解對流-擴散方程，計算毛管力和重力的影響。當流線追蹤完成后，先沿著流線求解對流方程，將沿流線求得的飽和度通過加權平均轉換至3D離散網(wǎng)格中，再用有限差分方法求解考慮毛管力和重力影響的擴散方程。

算子分裂方法需要通過原3D網(wǎng)格來計算流線之間的流體交換，在流線和網(wǎng)格之間進行飽和度轉換時會引入誤差，同時也縮短了流線更新的最大時間步長，降低了計算精度和效率[22]。為此，Tanaka 等[23-24]提出將毛管力和重力作用下的流量分解為沿著流線方向和垂直于流線方向兩類，將沿著流線方向的流量并入對流方程進行求解，再將垂直于流線方向的流量通過算子分裂在原3D網(wǎng)格中進行求解。

雖然流線模擬方法在計算對流-擴散類問題時的計算效率比對流主導流動時要低，但在模擬大型油藏模型時，流線模擬的計算速度仍然明顯快于有限差分等基于網(wǎng)格的數(shù)值模擬方法[3,24]。

2.3 考慮物理化學作用的流線模擬方法

我國老油田已普遍進入高含水、高采出程度、快速遞減的開發(fā)階段，提高采收率技術在我國老油田穩(wěn)產中具有核心地位。在化學復合驅等驅油過程中存在擴散、彌散、吸附、滯留、化學反應、熱傳導等一系列物理化學作用。流動過程對這些物理化學變化有影響，物理化學因素對流體流動也有影響[32]。考慮注入劑的物理化學及流體力學特性，使流線模擬技術在熱水驅、聚合物驅等問題上取得了一定突破。

Zhu等[27]將熱傳導方程引入流線方法中，使用算子分裂方法模擬熱水驅油過程。Al-sofi和Blunt[33]提出了基于流線模擬的聚合物非牛頓流動模擬方法。文章通過迭代方法求解非牛頓流動中壓力梯度、剪切應力、流動速度之間的相關關系。Thiele等[34]利用有限差分法計算聚合物驅中涉及到的聚合物運移與流體物性的相互關系，使用流線計算由注采速度、地質結構、油藏非均質性等因素決定的運移路徑，顯著提高了計算效率，降低了因數(shù)值彌散導致的誤差。

2.4 發(fā)展方向

通過使用有效密度和算子分裂等方法，流線模擬方法可以突破自身的限制求解可壓縮流體運移、擴散、熱傳導等問題。沿著流線求解流體運移問題的方法正在逐漸趨于成熟，能夠通過流線方法模擬的流體運移機制越來越多。在保證精度的條件下，提高流線模擬在計算復雜流體運移問題時的效率，是突破流線模擬技術瓶頸的重要課題。

3 現(xiàn)代流線模擬方法的應用

從油氣藏工程應用的角度看，現(xiàn)代流線模擬方法最重要的特征在于計算高效和流場識別[35]。發(fā)揮流線模擬方法在計算速度方面的優(yōu)勢，能降低數(shù)值模擬在計算機硬件和計算時間上的要求，為歷史擬合、地質模型篩選、降低油藏不確定參數(shù)等大型計算問題提供解決方案，為實時輔助決策等智能應用奠定基礎。發(fā)揮流場識別特色能夠建立注采井組與油氣藏、注采井組之間、油氣藏內部滲流通道之間的對應關系，對油藏經營管理中的開發(fā)井網(wǎng)部署、注采參數(shù)優(yōu)化、重建滲流場、提高經濟效益等具有理論指導意義和實際參考價值。這兩方面特色決定了現(xiàn)代流線模擬技術在解決油藏工程問題中的相關性和有效性，并推動流線模擬在應用層面的研究不斷發(fā)展。

3.1 發(fā)揮流線模擬計算性能的應用

計算速度是評價某一個油藏模擬方法實用性最重要的指標之一。如今百萬網(wǎng)格級別的地質模型已經十分常見，這使得數(shù)值模擬速度和油藏模型規(guī)模之間的矛盾更加突出，很大程度上限制了傳統(tǒng)的數(shù)值模擬技術在大型油氣藏開發(fā)中的應用[2]。另外，受限于計算速度，傳統(tǒng)的油藏模擬流程是一個確定的過程，對一個給定的油藏模型進行一次模擬計算，較少評價模型中的不確定參數(shù)對預測的影響。得益于計算能力的提高，現(xiàn)在的流線模擬技術已經可以動態(tài)模擬多個地質模型來篩選最接近實際的地質模型，流線模擬已經成為地質模型篩選[36-37]、大型油藏模擬[38-40]、注水開發(fā)優(yōu)化[35,40-41]、歷史擬合[42-43]等方面的核心技術，并在不斷發(fā)展完善之中。

Aarnes等[38]將多尺度模擬與流線模擬相結合，進一步提高了大型油藏模型的運算速度。此方法先使用混合多尺度模擬求解粗網(wǎng)格內的壓力和其子網(wǎng)格內的速度場，再根據(jù)子網(wǎng)格內的速度場求解全域流線，最后沿著流線解飽和度方程。這種方法大幅度提高了壓力方程的求解速度，避免了傳統(tǒng)有限差分模擬中必要的網(wǎng)格粗化，可以用來直接模擬大型或精細的油藏模型，有效保留了精細網(wǎng)格中重要的流動特征。此后，Natvig等[44]提出了一種新的多尺度與流線模擬相結合的方法用于模擬裂縫性油藏，更真實地反映出細小裂縫對流動的影響，飽和度變化與裂縫走向的關系。

近年來發(fā)展迅速的并行算法也可以與流線模擬方法相結合實現(xiàn)計算加速。并行算法是指在并行機上同時用很多個處理器聯(lián)合求解問題的方法和步驟[45]。流線模擬沿著1D流線求解運移問題，在一個時間步長內各流線之間的計算相互獨立，與并行算法的運算邏輯相一致。L?f等[46]首次建立了以流線為獨立單位的并行計算方法，從計算理念上提高了流線模擬的計算速度。

除與其他算法相結合實現(xiàn)計算提速之外，也可對流線模擬本身進行求解加速。Andrianov和Bratvedt[47]提出通過自適應網(wǎng)格加密算法加速求解沿著給定流線的1D運移問題。相比使用原始網(wǎng)格，自適應網(wǎng)格的計算速度提高了5倍。

另外，基于流線模擬技術發(fā)展得到了歷史擬合、計算剩余油分布、數(shù)值試井等問題的快速求解方法。Kam等[42]在歷史擬合過程中，使用了多尺度模擬和流線模擬相結合的方法，依據(jù)實際含水率、油氣比以及井底壓力的歷史數(shù)據(jù)調整地質模型參數(shù)。Batycky等[48]將流線與物質平衡方法相融合，提出了快速計算剩余油分布的方法。姚軍等[49-51]充分發(fā)揮流線模擬計算速度快的優(yōu)勢，依據(jù)油藏非均質性、多井、油水兩相流、堿-聚合物復合驅油等對流動規(guī)律的影響，建立了各自對應的流線數(shù)值試井解釋模型。

3.2 發(fā)揮流場識別特色的應用

流線模擬的流場識別功能能夠反映注采井組之間的受效關系，完井方式對油氣藏動態(tài)的影響，及油氣藏開發(fā)的層間和層內矛盾。以下將從這3個方面分別闡述流線模擬在油氣藏經營管理方面的典型應用。

發(fā)揮流線模擬流場識別特色的應用十分多樣，最常見的例子是評價每對注采井組之間的受效關系。這類應用通常以注采井組為研究對象調整注采方式，以達到擴大波及體積、提高驅油效率、提高油氣藏開發(fā)效益的目的[52]。基于流線優(yōu)化注采參數(shù)的經典方法有兩大類：以Batycky和Thiele[52]為代表的優(yōu)化注采分配系數(shù)；以Alhuthali等[53]、Park和Datta-Gupta[54]為代表的平衡各注采井組之間的飛行時間。

這類理念同樣可以應用在聚合物驅的優(yōu)化問題中。Clemens等[55]將油藏整體的優(yōu)化問題通過流線分解為各個注采井組之間的優(yōu)化問題，通過比較累積產油量與聚合物注入量之間隨時間的變化關系，篩選出受效較差的注采井組，再針對此井組進行注入?yún)?shù)調整，以達到提高聚合物驅驅替效率的目標。Choudhuri等[56]成功將流線模擬應用于一個實際大型油藏的聚合物驅優(yōu)化問題，評價并提高了注采井組的聚合物驅效率。Al-Sofi和Blunt[57]利用流線模擬并行計算的方法優(yōu)化聚合物驅參數(shù)，提高經濟效益。類似的流線模擬技術在優(yōu)化注氣[58]、我國低滲透油田注水[59]、老油田注水[60]、稠油注水[61]、交替注水注氣[62]等開發(fā)過程的優(yōu)化問題中也發(fā)揮了重要作用，取得了巨大的經濟效益。

在通過流線調整注采井組受效關系理念的基礎上，Tanaka等[63]沿著流線計算靜態(tài)和動態(tài)經濟參數(shù)，預測配注量對每對注采井組凈現(xiàn)值的影響，優(yōu)化以提高注采井組經濟效益為直接目標函數(shù)的方法，且計算效率及經濟效益更高。

從井與油藏動態(tài)的關系來看，井模型與流線模型的耦合能夠直觀展現(xiàn)完井方式和工作參數(shù)與油藏流體流動之間的對應關系，能夠通過調整井下設備或注采參數(shù)，動用或防止特定區(qū)域流體流入井內。Cao等[19,64]提出了水平井-油藏全耦合模型，對油藏和井筒分別建立模型求解，再通過2個模型共同計算的區(qū)域——近井油藏區(qū)域的流動條件進行迭代計算，實現(xiàn)了傳統(tǒng)油藏模型與水平井-近井模型的快速精確耦合求解。該方法能夠精確模擬射孔位置、井控設備等結構對流動方向的影響。如圖3所示，通過將流線模型和井模型相結合，能夠快速且精確地描述水平井完井方式、井下流量控制閥等對油藏中流場分布的動態(tài)影響。此技術有望在智能完井領域為完井方式優(yōu)化、流體控制設備調整等提供重要依據(jù)[63]。

從油藏開發(fā)的層間和層內矛盾來看，流線能夠直觀反映并量化分析油藏流場分布不均的情況，為擴大波及體積和提高驅油效率尋找最佳目標潛力區(qū)域。流線模擬在此方面最具代表性的應用是流動診斷技術。流動診斷是指對流線分布特征及其驅油效果進行量化評價的快速數(shù)值方法。其在實際油藏工程中最普遍的應用是研究開發(fā)方式對流場的影響，并以降低流場分布的非均質性為目的對開發(fā)方式進行優(yōu)化[65]。目前，對普通砂巖的水驅油的流動診斷方法已經較為成熟[65]，產生了巨大的經濟價值；近年來，隨著提高采收率技術的高速發(fā)展，聚合物驅替的流動診斷機制已被建立[66]，并已被成功應用于實際油藏開發(fā)之中。

圖3 流線模擬與水平井模型相耦合模擬水驅油過程[19]Fig.3 Simulation of water flooding process by coupling of streamline simulation and a horizontal well model[19]

肖康等[67]研究了優(yōu)勢通道下剩余油改變流線動用機制，提出了改變流線的動用策略，對剩余油進行挖潛。賈虎和鄧力琿[68]根據(jù)流線概念和水驅特征，選取流線空間位置、驅替效率、波及效率、水相驅動能力等方面的參數(shù)來量化表征流場分布特征，并通過流線聚類方法區(qū)分不同類型流場，據(jù)此準確找出潛在優(yōu)勢流場分布，能夠有效識別油藏中無效注水循環(huán)通道以及具有開發(fā)潛力的區(qū)域。

3.3 發(fā)展方向

根據(jù)以上討論可知，流線模擬技術的應用正處在快速發(fā)展之中。在油藏經營管理方面，綜合發(fā)揮流線模擬在計算高效和流場識別方面的特色，能夠直觀反映注采井組的受效關系、井對油藏動態(tài)的影響、油藏開發(fā)的層間和層內矛盾，提供來自井組、單井、整體油藏的綜合信息，建立配注方案、井網(wǎng)部署、流量控制閥開度等工程參數(shù)和地下油藏動態(tài)流場之間的相關關系，實現(xiàn)提高油藏開發(fā)經濟效益的目標。在油氣田智能化發(fā)展的大背景下，流線模擬技術擁有巨大的發(fā)展?jié)摿Α?/p>

4 結 論

流線追蹤的實質在于離散網(wǎng)格內速度場的逼近，追蹤方法由低階方法向高階方法發(fā)展，適用網(wǎng)格由規(guī)則網(wǎng)格向多邊形不規(guī)則網(wǎng)格轉變，模擬對象由普通砂巖油藏向復雜斷層、縫洞型油藏、非常規(guī)油氣藏進行拓展，模擬區(qū)域由油藏向近井地帶進行擴展，這樣的發(fā)展變化正處在不斷加速之中。此外，通過使用有效密度和算子分裂等方法，流線模擬技術已經能夠應用于可壓縮流體，解決對流-擴散類問題，模擬驅替機理較為復雜的流體運移問題。

在數(shù)值模擬應用中，流線方法相對于傳統(tǒng)數(shù)值模擬方法具有計算高效和流場識別兩方面的特色。通過發(fā)展基于流線的多尺度模擬、并行算法、適應網(wǎng)格加密算法等，計算速度得到了進一步的提升；以流線為基礎的歷史擬合、剩余油分布模擬、試井解釋等問題的求解速度也得到了提高。

在油藏經營管理中，流線像一座橋梁，通過計算流線路徑和飛行時間，建立注采井組之間的受效關系，反映單井對油氣藏動態(tài)的影響、油氣藏開發(fā)的層間和層內矛盾，將地下油藏流場動態(tài)與人工可控的工程參數(shù)緊密相連。調整井組時，可沿著流線計算靜態(tài)和動態(tài)經濟參數(shù)，預測配注量對每對注采井組凈現(xiàn)值的影響，實現(xiàn)井組效益最大化。優(yōu)化單井時，流線模擬和井模型相結合能夠描述單井完井方式、工作制度等對油藏中流場的動態(tài)影響，從泄油面積和驅油效率的角度出發(fā)，提高單井效益。在對油藏流場進行整體調整時，流線能夠直觀反映并量化表征油藏流場分布不均的情況，有效識別油藏中無效注水循環(huán)通道，為提高驅油效率尋找最佳目標潛力區(qū)域，提高油藏整體的動用程度和經濟開發(fā)效益。