Lagrange插值對相滲曲線歸一化方法的改進
——以M油田相滲曲線為例

路 強,任宜偉,尹彥君,張 震,王中華,陳秋月

(中海油能源發(fā)展股份有限公司工程技術(shù)分公司,天津 300457)

相對滲透率指巖石孔隙中飽和多相流體時,巖石對每一相流體的有效滲透率與巖石絕對滲透率的比值。油水相對滲透率與含水飽和度的關(guān)系曲線稱為相對滲透率曲線[1],即相滲曲線。該曲線是研究儲層多相流動的重要參數(shù),是油田開發(fā)動態(tài)分析、開發(fā)方案設(shè)計、油藏數(shù)值模擬等工作的基礎(chǔ)資料[2]。油藏開采初期,實驗室需要測定多條相對滲透率曲線,為了獲得一條能夠準(zhǔn)確反映儲層流體滲流特征的相對滲透率曲線,需要對相對滲透率曲線進行必要的插值,以得到一條歸一化的曲線。對于不同類型的相滲曲線應(yīng)該如何進行插值計算,針對這一問題,本文剖析了Lagrange插值對相滲曲線歸一化的改進。

1 相滲曲線的分類

不同儲層的油相相滲曲線變化趨勢基本一致,表現(xiàn)為一條光滑的單調(diào)遞減曲線,而水相相滲曲線則差異較大,王國先[3]按照水相相滲曲線形態(tài)將相滲曲線分為五類,中高滲儲層常見的是水相上凹型(圖1a),低滲儲層常見的是水相上凹型和水相直線型(圖1b),低滲敏感性儲層常見的是水相下凹型(圖1c)和水相上凸型(圖1d),低滲裂縫型儲層主要是水相靠椅型(圖1e)。

圖1 不同類型相滲曲線

水相上凹型:實驗室測得的相滲曲線絕大多數(shù)都是水相上凹型,這類曲線為標(biāo)準(zhǔn)型油水相對滲透率曲線。隨含水飽和度增加,水相曲線光滑、無拐點、單調(diào)遞增,曲線近三次多項式。

水相直線型:曲線形態(tài)與水相上凹型曲線相似,也較為常見,受儲層孔隙填充物膨脹影響,曲線不上凹,隨含水飽和度增加,水相曲線光滑、無拐點、單調(diào)遞增,曲線近一次多項式。

水相下凹型:曲線反映出儲層黏土礦物含量高且已發(fā)生鹽(水)敏,黏土礦物遇水膨脹堵塞喉道,水相流動阻力增大,隨含水飽和度增加,水相曲線光滑、無拐點、單調(diào)遞增,曲線近二次多項式。

水相上凸型:水相上凸型的水相曲線與水相下凹型曲線相似,但尾部明顯下掉,曲線形態(tài)近似駝峰,隨含水飽和度增加,水相曲線光滑、一個拐點,曲線近二次多項式。

水相靠椅型:裂縫較發(fā)育儲層的相滲曲線多表現(xiàn)為水相靠椅型,隨含水飽和度增加,水相曲線光滑、兩個拐點,曲線近三次多項式。

由于一次多項式和二次多項式是三次多項式的特殊形式,因此,五類相滲曲線的水相曲線形態(tài)均與三次多項式曲線近似。

2 相滲歸一化處理方法及改進

相滲歸一化處理有三種方法:平均法、經(jīng)驗公式法[4-8]、相對滲透率特征曲線法[9]。平均法是相滲曲線歸一化處理的一種常用方法,該方法不需要對實驗數(shù)據(jù)進行擬合,能最大程度地保存原始相滲曲線的形態(tài),但是潘婷婷[10]認(rèn)為該方法精度低且曲線不光滑。

2.1 平均法

平均法包括平均相對滲透率法和平均飽和度法[11],兩種處理方法類似,處理過程如下:

1)選取具有代表性的油水相對滲透率曲線數(shù)據(jù);

2)對各巖心樣品的實驗數(shù)據(jù)進行無因次化處理(式1),并繪制無因次后的油水相對滲透率曲線;

(1)

(2)

4)將各樣品的相滲曲線特征值分別進行算術(shù)平均(式3),并將平均值作為平均相滲曲線的特征值;

(3)

5)將求得的平均無因次相對滲透率曲線換算成無因次化之前的相滲曲線,其公式為:

(4)

2.2 平均法的改進

Lagrange插值法原理是:在插值點的附近選取若干合適的結(jié)點,構(gòu)造一個簡單的插值函數(shù)y=L(x),要求L(x)穿過選取的結(jié)點。在所選數(shù)據(jù)區(qū)間用L(x)作為原函數(shù)f(x)的值,使得f(xi)=L(xi)(i=0,1,…,n-1)成立[16]。

1)構(gòu)建Lagrange插值基函數(shù)。給定n個結(jié)點xi(i=0,1,…,n-1)上的函數(shù)值yi=f(x),則插值基函數(shù)的表達式為:

(5)

基函數(shù)滿足lj(xj)=1,lj(xi)=0(i≠j)。

2)n個結(jié)點構(gòu)建的Lagrange多項式L(x),其表達式為:

(6)

其特點是滿足L(xi)=yi(i=0,1,…,n-1)。

2.3 Lagrange插值法的改進及實現(xiàn)

鄭洪印[17]1993年提出利用Lagrange插值法進行插值,插值時選取距離插值點最近的三個相鄰結(jié)點進行Lagrange插值,但該方法存在以下問題:三個插值結(jié)點構(gòu)造成二次多項式,如果多項式對稱軸介于三個插值結(jié)點之間,曲線會出現(xiàn)拐點。

一般情況下,插值結(jié)點越多,預(yù)測結(jié)果越準(zhǔn)確,但插值次數(shù)越多,越容易出現(xiàn)龍格現(xiàn)象。五類相滲曲線的形態(tài)均與三次多項式曲線近似,4個結(jié)點可以確定唯一的三次多項式,因此在n個結(jié)點中選取距離插值點最近的4個結(jié)點進行插值,且使指定插值點t位于4個結(jié)點之間。

3 實例分析

以M油田相滲曲線為例,實驗室共測得四塊巖心的相滲曲線數(shù)據(jù)(表1),每塊巖心選取三個實驗數(shù)據(jù)點作為盲點,對比線性插值法、四點Lagrange插值法、三點Lagrange插值法的精度。

表1 H油田巖心相滲曲線數(shù)據(jù)

結(jié)果表明:線性插值法盲點最高相對誤差可達40%,油相相對滲透率平均誤差為6.13%～17.95%,水相相對滲透率平均誤差為0.88%～4.83%;三點Lagrange插值法盲點最高相對誤差為13.33%,油相相對滲透率平均誤差為0.83%～4.99%,水相相對滲透率平均誤差為0～0.1%;四點Lagrange插值法盲點最高相對誤差為0.41%,油相相對滲透率平均誤差為0～0.14%,水相相對滲透率平均誤差為0～0.09%(表2)。因此,Lagrange插值法的精度明顯高于線性插值法,而四點Lagrange插值法精度高于三點Lagrange插值法。

表2 不同插值法相對誤差對比

4 結(jié)論

1)五類典型相滲曲線形態(tài)與三次多項式曲線形態(tài)類似,Lagrange插值法可以構(gòu)建一個多項式,保證實驗數(shù)據(jù)點誤差為零,并且能夠最大程度保持每一類相滲曲線的形態(tài)。

2)采用線性插值法的相對誤差最高,而Lagrange插值法的相對誤差較小。相較于三點Lagrange插值法,四點Lagrange插值法構(gòu)造的曲線為三次多項式,曲線形態(tài)更符合五類相滲曲線的形態(tài),計算精度更高。