非線性數(shù)值模擬進展評價

1、非線性數(shù)值模擬進展評價通過對程林松、姜瑞忠、黃延章、何順利、姚軍等人的相關文章調研發(fā)現(xiàn)，目前對于低速非達西模型的描述主要為啟動壓力梯度和毛管力動態(tài)效應兩種表述方式。1. 單相啟動壓力梯度大量的實驗表明，關于低滲透性油層中滲流特征和滲流規(guī)律的表述可歸納為：具有啟動壓力梯度的非線性律，如圖1所示。其基本滲流特征為：(1) 當壓力梯度小于某個值時，流體不流動。這個壓力梯度稱為最小啟動壓力梯度(真實啟動壓力梯度)，以a 表示.(2) 當壓力梯度在比較低的范圍時，滲流速度的增加呈下凹型曲線。(3) 當壓力梯度比較高時，滲流速度呈直線性增加。(4) 該直線段的延伸與壓力梯度軸交于某點而不通過坐標原點. 此

2、點稱為平均啟動壓力梯度 (擬啟動壓力梯度)，以b表示, c為最大啟動壓力梯度。圖1 滲流速度與壓力梯度的關系目前常用的單相啟動壓力梯度描述模型主要分為擬啟動壓力梯度模型、分段模型和連續(xù)模型。1.1 擬啟動壓力梯度模型一些研究人員基于擬啟動壓力梯度模型，針對擬啟動壓力梯度的求解方法進行了實驗和理論研究，實驗流體從單一流體到多種流體，考慮的影響因素也從單一的滲透率變?yōu)闈B透率、粘度等多種因素的綜合作用。李忠興，程林松1-3等人（2004）使用抽真空煤油進行室內驅替實驗，依據(jù)考慮擬啟動壓力梯度的達西定律：Q=KAPL-對兩邊取常用對數(shù)得到：lgPL-=lgQ-lgKA該式中流量為零時驅替壓力梯度即為啟

3、動壓力梯度，通過回歸實驗數(shù)據(jù)，得出冪律關系式：=K-n最后研究人員繪制了確定低滲透油田啟動壓力梯度的理論圖版，如圖2所示。該求解方法僅考慮了滲透率和啟動壓力梯度的關系。郝斐，程林松4等人（2006）使用模擬油、地層水、注入水和蒸餾水進行物理模擬，實驗結果表明：雙對數(shù)坐標中啟動壓力梯度與滲透率呈線性關系且斜率為-1(圖3所示)，這與理論研究相符合：lgG=-lgk+lgb×驗證了壓力梯度與滲透率的冪函數(shù)關系，并且研究了不同粘度對啟動壓力梯度的影響，結果表明粘度越大，啟動壓力梯度越大。許建紅，程林松5-6等人（2007）使用抽真空煤油進行驅替實驗，擬合得到驅動壓力梯度與滲流速度的關系式：

4、v=ap2+bp+c pG再與考慮擬啟動壓力梯度的達西定律進行聯(lián)立求解，可以得到下式：Gi=-aP2+b-KiiP+cKii該求解方法不僅考慮了滲透率與啟動壓力梯度的關系，還考慮了流體粘度、驅動壓力梯度對啟動壓力梯度的影響。圖2 確定低滲透地層啟動壓力梯度大小理論圖版圖3 不同流體啟動壓力梯度與滲透率的關系也有一些研究人員討論了考慮擬啟動壓力梯度的數(shù)學模型，并進行了數(shù)值模擬或者尋找簡化方法求得解析解。韓洪寶，程林松1-3,6-9等人（2004）引入了邊界層理論5，將油相粘度用表觀粘度來代替：B=A1+1-A2同時考慮了擬啟動壓力梯度（如圖4點C），將結果應用到黑油模型中，得到了修正后的三維三相

5、黑油模型：油方程：kKroBoBpo-ogD-Go+qoosc=tSoBo 水方程 ：kKrwBwwpw-wgD+qwwsc=tSwBw 氣方程：kKroRsoBoBpo-ogD-Go+kKrwRswBwwpw-wgD+kKrgBggpg-ggD+qggsc=tSoRsoBo+SwRswBw+SgBg圖4 典型的非達西滲流關系曲線李莉，董平川等人（2006）基于擬啟動壓力梯度模型，建立了三維油水兩相非達西滲流數(shù)學模型，利用差分法得到數(shù)值模型，并設計了正方形注采單元進行數(shù)值模擬，通過改變注采井距，的餓到了見水時地層壓力、壓力梯度等參數(shù)的空間分布，最后得到了有效驅動注采井距。一些研究人員也以擬啟動

6、壓力梯度為基礎建立了低滲透油藏非達西滲流模型，使用不同的方法探索模型的精確解。郝斐，程林松10-12等人（2006）建立了考慮擬啟動壓力梯度低滲透油藏不穩(wěn)定滲流模型，并在內邊界條件中考慮了動邊界：pr-Gr=Rt=0p=pe rRt利用油藏壓力近似表達式：p=a0lnrRt+a1+a2rRt,rwrRt求得近似解和半解析解，分析了啟動壓力梯度對地層能量傳播和流體運移的影響。兩種求解結果有較好的一致性，解決了數(shù)值解法計算量大的弊端。李松泉，程林松13等人（2008）在前人研究的基礎上，建立了考慮擬啟動壓力梯度和介質變形的穩(wěn)定、不穩(wěn)定滲流數(shù)學模型，利用近似求解方法，給出了定產、變產和定壓條件下的解

7、。同時也低滲透油藏產能遞減規(guī)律、井距確定和含水率變化進行了分析。含有動邊界條件的低滲透油藏非達西不穩(wěn)定滲流數(shù)學模型具有較強的非線性，難以求出精確解析解。劉文超，姚軍14-17等人（2013）建立了能反映動邊界移動規(guī)律的低滲透油藏不穩(wěn)定非達西滲流數(shù)學模型，模型中不含動邊界條件，經過無因次化，模型如下：pDtD=1rDvDpDrD-rDvD'pDrD×2pDrD2pDtD=0=0rDpDrDrD=1=-1+DpDrD=reD=0vDpDrD=0 0-pDrDDpDrD+D -pDrD>DvD'pDrD=0 0-pDrDD-1 -pDrD>D通過采用Dougla

8、s-Jones預估校正有限差分法，求得了該模型瞬時壓力的精確數(shù)值解。計算結果表明井底瞬時壓力的雙對數(shù)曲線存在拐點，拐點前的曲線反映動邊界到達定壓邊界前的狀態(tài)，拐點后的曲線反映動邊界到達定壓邊界后的狀態(tài)，如圖5所示。圖5 瞬時井底壓力的雙對數(shù)曲線上述關于擬啟動壓力梯度問題的探討存在一定的局限性：井-網格方程未考慮擬啟動壓力梯度或因處理太簡單而不合理，流體在網格、網格之間的流動與井、網格之間的流動會產生矛盾，前者可使用模型直接離散化實現(xiàn)，后者需要對井附近的流動進行解析研究，建立井-網格方程，從而修正了內邊界條件。趙國忠18,19（2006）建立了考慮擬啟動壓力梯度的三維三相滲流數(shù)值模型，通過等效半

9、徑（如圖6所示）re=xy/12將井與其所在網格關聯(lián)起來。圖6 單位時間步內的井-網格穩(wěn)定流動示意圖王建忠，姚軍等人（2013）建立了三維兩相數(shù)值模型，利用Peaceman等效半徑將網格壓力和井底壓力聯(lián)系起來：pi,j,k-pwf=Q2KxKy12hlnrorw+ro-rwG進而推導了考慮啟動壓力梯度影響的定井底壓力、定產油量和定產液量條件下的內邊界條件。現(xiàn)有的低滲透油藏非線性流動數(shù)值模擬，大都采用傳統(tǒng)的有限差分法或控制體積有限差分，此類方法基于相鄰單元中心點的物理量來構造數(shù)值計算個事，求取壓力梯度，進而進行求解。而非均質低滲油藏中不同區(qū)域甚至不同網格單元的滲流模式不盡相同，必須獲取各個網格單

10、元的壓力梯度，并選擇相應滲流模式，因此傳統(tǒng)的數(shù)值模擬方法無法滿足要求。姚軍，黃濤等人（2014）基于擬啟動壓力梯度模型，結合動態(tài)滲透率概念，建立了單相數(shù)值模型。在控制體積有限差分法的基礎上引入界面中心點壓力和速度，提出一種新的數(shù)值算法，該算法基于單個網格來構造數(shù)值計算格式，課精確獲得每個網格的壓力梯度，單元間通過流量連續(xù)條件聯(lián)系起來。總結：擬啟動壓力梯度模型較為成熟，但不能體現(xiàn)出非線性滲流的特點，且當壓力梯度小于擬啟動壓力梯度時油藏流體無法動用，大大縮小了特低滲油藏的流動范圍。1.2 分段模型分段模型起步較早，黃延章22（1998）總結出目前主要的三種分段模型為：直線逼近法、冪律函數(shù)逼近法和二

11、次函數(shù)逼近法，如下式所示。直線逼近法 v=0              PLa v=K1PL-a  a<PL<b v=K2PL-c bPL冪律函數(shù)逼近法v=0                       PLav=a1

12、KPL-Gia2    a<PL<b v=K2PL-c         bPL二次函數(shù)逼近法v=0              PLav=APL2+BPL+c a<PL<b v=K2PL-c            &#

13、160;      bPL尹芝靈，孫文靜等人（2011）以分段模型為基礎，結合了姚軍教授提出的動態(tài)滲透率概念：K=f1pKe GA<pm<Gcf1p=2ap+b2aGc+b建立了涵蓋非線性段和擬線性段的油藏數(shù)值模擬模型，其運動方程為：v1=0, pl<GA f1plKmaxKrllpl-GA, GA<pl<GcKmaxKrllpl-GB, pl>Gc通過全隱式差分格式建立了數(shù)值模型，采用預處理共軛梯度法進行求解。結果表明：啟動壓力梯度對壓力場有顯著地影響，且用此非線性滲流模型計算的壓力變化比擬線性情

14、況平緩，更符合油藏開發(fā)實際（如圖7所示）圖7 不同滲流模型條件下注采井間的壓力分布鑒于現(xiàn)有分段模型的連續(xù)性和確定性問題，難以直接應用于建模，劉文超，姚軍24等人（2012）從低滲透多孔介質滲透率漸變理論出發(fā)，將冪律分段模型進行了公式化，得到了非線性運動方程實數(shù)域內的表示：v=0, dPdrCv=-dPdrdPdrkA-BA-CndPdr-Cn,C<dPdrAv=-dPdrdPdrkdPdr-B, dPdrA該表達式中滲流速度關于壓力梯度的一階導數(shù)是連續(xù)的，得到的是滲透率和擬啟動壓力梯度也是連續(xù)變化的，符合低滲介質滲透率連續(xù)變化的客觀規(guī)律：kb=-v'Gb=dPdr-vv'

15、之后，通過建立單相非線性徑向流動的數(shù)值模型，采用Douglas-Jones預估校正發(fā)求得了精確數(shù)值解，結果表明非線性滲流模型世界語擬線性和達西模型之間的一種理想模型?？偨Y：分段模型比較精確，描述了不同階段的滲流規(guī)律，對于精確的科學研究可以使用。但是應用較為麻煩，并且對于各滲流階段的臨界點判斷較難。1.3 連續(xù)模型鄧英爾，劉慈群25（2001）首次提出了連續(xù)函數(shù)模型，經典的三參數(shù)模型，為非線性滲流定量分析奠定了基礎：va1+a21+bv=-P該模型的缺點是無法反應滲流存在最小啟動壓力梯度的現(xiàn)象。楊清立，楊正明26等人（2007）提出了一個新的模型，既反映了最小啟動壓力梯度現(xiàn)象，又很好地描述了非線

16、性段滲流特征，該模型如下：v=KP1-1a+bP時宇，楊正明27（2008）根據(jù)大慶油田不同區(qū)塊低滲透巖心的實驗結果，特出了反應真實啟動壓力梯度和非線性滲流段的兩參數(shù)模型：v=Kgradp1-b1+agradp以上模型為早期的連續(xù)函數(shù)模型，這些模型主要利用唯像學方法，根據(jù)非線性滲流物理現(xiàn)象擬合出數(shù)學方程，難以反映產生非線性的內在原因和影響因素。時宇，楊正明28（2009）依據(jù)毛細管模型和邊界層理論，建立了低滲透非線性滲流模型，定量描述了產生非線性滲流的內在因素，模型如下：v=ackP8P-cpPckPPminckP=rmax4-rmin44 PPmaxrmax4-r44Pmax>PPmi

17、ncpP=rmax3-rmin33 PPmaxrmax3-r33Pmax>PPmin楊仁峰，姜瑞忠等人29-32（2011）提出了一個新模型，并建立了相應的數(shù)學模型，編制了三維三相非線性滲流數(shù)值模擬軟件，為非線性滲流模型的工程化應用奠定了基礎，該模型如下：Q=KA1-1P-12PP-2P該模型引入的兩個新參數(shù)具有豐富的內涵，體現(xiàn)了屈服應力和邊界層對滲流的影響，較好地體現(xiàn)出低滲透油藏最小啟動壓力梯度 、壓力梯度隨驅替壓力梯度增加而增加、非線性滲流等滲流特征。作者同時考慮了壓敏效應，建立了三維三相數(shù)學模型，構造了響應的有限差分離散化格式，確定了考慮非線性滲流的井-網格流動方程。姜瑞忠，李林凱

18、等人33（2012）基于毛細管模型，結合邊界層理論，通過引入描述低滲透儲層滲流的特征參數(shù)c1和c2，建立了新的數(shù)學模型，從根本上解釋了啟動壓力梯度和非線性滲流產生的原因，該模型如下：v=K1-c1P-c2P該模型及包含了達西模型，又包含了擬啟動壓力梯度模型，且反映了最小啟動壓力梯度的現(xiàn)象。黃延章，楊正明等人34（2013）提出了一個新的非線性滲流模型，該方程各項參數(shù)由實驗直接得出，且含義明確，是一個通用方程，模型如下：v=KodPdx1-cdPdx+c-a 該模型是一個通用方程，當a=c0時，該方程為賓漢流體的滲流方程；當a=0時，該方程變成無真實啟動壓力梯度的非線性滲流方程；當a=c=0時，

19、該方程為達西線性滲流方程。 以上模型為近期的連續(xù)函數(shù)模型，此類模型適用性較強，具有連續(xù)性易于離散后建立數(shù)值模型，發(fā)展前景大；但是模型中的參數(shù)均需要實驗或者擬合得到，且有些參數(shù)具體含義不明確。2. 兩相啟動壓力梯度蘇法卿，孫志剛35等人（2004）實驗研究了油水兩相最小啟動壓力梯度，實驗發(fā)現(xiàn)：（1） 最小啟動壓力梯度與孔隙滲透率呈冪函數(shù)關系；（2） 最小啟動壓力梯度隨含水飽和度升高先增大后減?。唬?） 兩相啟動壓力梯度比單相的值增大510倍。郝斐，程林松4等人（2006）比較了單相油、水及束縛水、殘余油狀態(tài)下的啟動壓力梯度。單相滲流時，啟動壓力梯度較?。粌上囹屘鏁r，啟動壓力梯度值較大，即隨著含水

20、飽和度增大，兩相啟動壓力梯度逐漸增大，在某一含水飽和度時達到最大值，隨著含水飽和度繼續(xù)增大，啟動壓力梯度反而減小，如圖8所示。分析認為，兩相啟動壓力梯度大很多的原因主要是毛管力影響的結果。圖8 典型巖樣啟動壓力梯度與含水飽和度關系鄧玉珍，劉慧卿36（2006）通過室內試驗測定了油水兩相啟動壓力梯度，提出了適合兩相滲流的啟動壓力梯度模型，認為兩相啟動壓力梯度與流度呈冪指數(shù)關系，其中引入了油水兩相綜合流度的概念，兩相流體粘度用油水兩相粘度的飽和度加權形式。G=k-n隨著流度增加，啟動壓力梯度逐漸減小，如圖9所示。圖9 油水兩相啟動壓力梯度與流體流度的關系李松泉，程林松13等人（2008）實驗測量了

21、束縛水和殘余油飽和度狀態(tài)下的啟動壓力梯度值，結果與郝斐（2006）的結果一直，分析原因是賈敏效應使得連續(xù)的液流變?yōu)榉稚⒌囊旱危蟠笤黾恿怂層退璧膲毫?，反映在啟動壓力梯度的增大。同時，作者在油相和水相運動方程中加入擬啟動壓力梯度，給出了油水兩相非活塞驅替的分流量、壓力、產量方程。時宇，楊正明27等人（2008）根據(jù)巖心實驗提出了單相啟動壓力梯度的兩參數(shù)模型，并在B-L方程的基礎上建立了非線性非混相驅替方程進行了求解，得出的含水率公式如下：fw=krwswwkrwsww+kroswo1+badpdx-1楊仁峰，姜瑞忠等人37（2011）深入分析了相對滲透率引入的基本假設，將建立的單相模型推廣到

22、兩相體系，得出了多相非線性滲流新模型：Ql=KKrllA1-1lPl-1l2lPlPl-2lPl式中1l、2l應取單相的值。同時，作者分析了現(xiàn)有兩相非線性滲流理論研究方面存在的誤區(qū)，指出現(xiàn)有連續(xù)模型中的參數(shù)均應取單相滲流時的測量值?？偨Y：現(xiàn)有關于兩相啟動壓力梯度的研究多為實驗測量并定性地分析規(guī)律，而關于非線性滲流數(shù)學模型也僅僅是將單相模型推廣到兩相模型的油相運動方程中，并沒有專門描述兩相啟動壓力梯度的數(shù)學模型。3. 毛管力動態(tài)效應多數(shù)國內學者只研究了油水界面平衡狀態(tài)下的靜態(tài)毛細管壓力及其對低滲透油藏水驅油效果的影響，認為毛細管壓力只是含水飽和度的函數(shù)。以Hassanizadeh 38-44為代

23、表的國外學者研究了動態(tài)毛細管壓力的效應，認為在非穩(wěn)態(tài)運動過程中，毛細管壓力不斷變化，其不僅是濕相流體飽和度的函數(shù)，還受到濕相流體飽和度變化率的影響，即：Pc=Po-Pw=fSw,Swt Hassanizadeh44(2002)提出動態(tài)毛管力與含水飽和度變化率有以下關系：Pc,dyn=Po-Pw=Pc,statSw+10-6Swt田樹寶，何順利45等人（2012）首次實驗研究了低滲透油藏中毛細管理的動態(tài)效應，實驗發(fā)現(xiàn)巖心滲透率和毛細管動態(tài)系數(shù)呈反比，如圖10所示；同時，作者建立了數(shù)值模型，考慮了動態(tài)毛細管壓力，沒有單獨引入啟動壓力梯度，即將Pc,dyn=Po-Pw=Pc,statSw+

24、10-6Swt引入了輔助方程，對模型進行了求解，計算結果與滲流實驗規(guī)律基本吻合，表明低滲透油藏油水兩相動態(tài)毛細管壓力是導致油水兩相產生非達西流動的原因之一。 （a） （b）圖10 不同進汞速度條件下測定的低滲透巖心和高滲透巖心毛細管壓力曲線對比舒衛(wèi)兵，許鶴華46等人（2014）建立了二次函數(shù)非線性滲流模型，自定義了三種不同變化率的飽和度時變曲線，如圖11所示。Sw1=0.1T Sw2=0.4167lnT+1 Sw3=0.01T2圖11 數(shù)值模擬采用的三種飽和度時變曲線通過有限元數(shù)值模擬一維驅油過程，研究了毛管力動態(tài)效應。結果表明：三種不同飽和度時變曲線得到的毛細管壓力存在明顯的差異（圖12所示

25、）圖12 模擬所得到的毛細管壓力曲線同時，計算結果表明動態(tài)系數(shù)與滲透率呈負相關，與潤濕相流體粘滯系數(shù)呈正相關。參考文獻：1 李忠興，程林松.特低滲油藏啟動壓力梯度新的求解方法及應用.2 馮文光，葛家理.單一介質、雙重介質中非定常非達西低速滲流問題.3 馮文光.非達西低滲滲流的研究現(xiàn)狀與展望.4 郝斐，程林松.特低滲透油藏啟動壓力梯度研究.5 許建紅，程林松.低滲透油藏啟動壓力梯度新算法及應用.6 韓洪寶，程林松.特低滲透油藏考慮啟動壓力梯度的物理模擬及數(shù)值模擬方法.7 黃延章，嚴衡文，楊正明.低滲透油層滲流機理.8 鄧英爾，閆慶來.界面分子力作用于滲透率的關系及其對滲流的影響.9 黃延章.低滲

26、透油層分線性滲流特征.10 郝斐，程林松等.考慮啟動壓力梯度的低滲透油藏不穩(wěn)定滲流模型.11 李凡華，劉慈群.含啟動壓力梯度的不定常滲流的壓力動態(tài)分析.12 程時清，徐論勛等.低速非達西滲流試井典型曲線擬合法.13 李松泉，程林松.特地滲透油藏非線性滲流模型.14 劉文超，姚軍等.低滲透油藏非達西滲流瞬時壓力計算方法.15 程時清，黃延章等.低速非達西滲流流動邊界問題的積分解.16 郭永存，曾清紅等.動邊界低滲透油藏的無網格方法數(shù)值模擬.17 Wu J S.Flow and displacement of Bingham non Newtonian fluids in porous media

27、18 趙國忠.變啟動壓力梯度三維三相滲流數(shù)值模擬方法.19 郭永存，曾清紅燈.有啟動壓力梯度滲流的數(shù)學模型.20 王建忠，姚軍等.考慮井網格影響的油水兩相低滲油藏數(shù)值模擬.21 姚軍，黃濤.非均質低滲透油藏非線性流動數(shù)值模擬.22 尹芝靈，孫文靜.動態(tài)滲透率三維油水兩相低滲透油藏數(shù)值模擬.23 姚軍，劉順.基于動態(tài)滲透率效應的低滲透油藏試井解釋模型.24 劉文超，姚軍等.基于滲透率連續(xù)變化的低滲透多孔介質非線性滲流模型.25 鄧英爾，劉慈群.低滲油藏非線性滲流規(guī)律數(shù)學模型及其應用.26 楊清立，楊正明等.特低滲透油藏滲流理論研究.27 時宇，楊正明.低滲透油藏非線性兩相滲流研究.28 時宇，楊正明.低滲透儲層非線性滲流模型研究.29 楊仁峰，姜瑞忠等.特地滲透油藏非線性滲流數(shù)值模擬.30 楊清立，楊正明等.特地滲透油藏滲流理論研究.31 徐運亭，徐啟等.低滲透油藏滲流機理研究及應用.32 楊正明，于榮澤等.特地滲透油藏非線性滲流數(shù)值模擬.33 姜瑞忠，李林凱等.低滲透油藏非線性滲流新模型及試井分析.34 黃延章，楊正明等.低滲透多孔介質中的非線性滲流理論.35 蘇法卿，孫志剛等.低滲透砂巖油藏油水兩相最小啟動壓力梯度實驗研究.36 鄧玉珍，劉慧卿.低滲透巖心中油