呂 晶，劉顯太，孫業(yè)恒，陳德坡，黃迎松

（1.中國(guó)石化勝利石油管理局有限公司博士后科研工作站，山東東營(yíng) 257001；2.中國(guó)石化勝利油田分公司，山東東營(yíng) 257001；3.中國(guó)石化勝利油田分公司勘探開發(fā)研究院，山東東營(yíng) 257015）

水相分流量曲線是研究油藏開發(fā)動(dòng)態(tài)的重要基礎(chǔ)資料之一，利用水相分流量曲線能夠分析油藏含水率上升規(guī)律、水驅(qū)油藏含水飽和度分布、地層平均含水飽和度及采出程度、驅(qū)油效率等重要開發(fā)指標(biāo)，對(duì)分析油藏開發(fā)動(dòng)態(tài)意義重大［1-3］。目前水相分流量曲線的獲取方法是利用室內(nèi)相對(duì)滲透率測(cè)試數(shù)據(jù)，通過水相分流量定義式計(jì)算得到［4-5］。由于室內(nèi)實(shí)驗(yàn)驅(qū)替倍數(shù)有限，不能直接獲得殘余油飽和度，水相分流量方程的末端值是估算所得。理論上水相分流量的數(shù)值為0～1，但現(xiàn)階段針對(duì)水相分流量方程的近似擬合模型，如經(jīng)典邏輯斯蒂方程，其方程形式?jīng)Q定了水相分流量計(jì)算結(jié)果在初始時(shí)刻無法等于0，同時(shí)曲線末端不能等于1，只能無限逼近1，無法計(jì)算殘余油飽和度，水相分流量曲線的末端值仍無法確定。

針對(duì)水相分流量擬合以及曲線末端含水飽和度取值的問題，以室內(nèi)一維小巖心驅(qū)替實(shí)驗(yàn)為研究對(duì)象，根據(jù)水相分流量曲線的實(shí)際特征，結(jié)合水相分流量的實(shí)際物理意義及水相分流量在理論上的初始條件和邊界條件，通過對(duì)比多種S型函數(shù)，確定并建立滿足水相分流量的邊界條件的全過程數(shù)學(xué)模型。該模型通過2 個(gè)調(diào)整指數(shù)，保障曲線S 型特征與實(shí)測(cè)水相分流量分布特征一致。同時(shí)，圍繞水相分流量末端含水飽和度預(yù)測(cè)的問題，通過設(shè)置極限含水飽和度參數(shù)，實(shí)現(xiàn)利用實(shí)測(cè)資料擬合直接獲得殘余油飽和度。該方法利用相對(duì)滲透率測(cè)試報(bào)告即可實(shí)現(xiàn)殘余油飽和度的預(yù)測(cè)，具有較為廣泛的適用性?；谠撃Ｐ涂色@得相對(duì)滲透率比值全過程曲線，能夠明確驅(qū)替從初始時(shí)刻到極限時(shí)刻的油水運(yùn)動(dòng)規(guī)律變化特征，為開展水驅(qū)油藏開發(fā)動(dòng)態(tài)研究，特別是高含水老油田注水開發(fā)動(dòng)態(tài)特征研究和長(zhǎng)期注水導(dǎo)致的油藏局部高耗水等問題的研究提供理論基礎(chǔ)和指導(dǎo)。

1 水相分流量曲線

1.1 形態(tài)特征及物理意義

根據(jù)油藏工程理論，水相分流量的物理意義為任意位置的水相流量與油水兩相流量之和的比值，反映油水兩相流動(dòng)中水相的流動(dòng)能力［6］。水相分流量通常用fw表示，數(shù)值范圍介于0 到1，函數(shù)表達(dá)式為：

（1）式表明，水相分流量與油水黏度比和油水相相對(duì)滲透率比值有關(guān)，對(duì)于特定油藏，油水黏度比可近似認(rèn)為是定值，因此，水相分流量的大小和變化過程實(shí)質(zhì)上反映了油水相相對(duì)滲透率比值隨驅(qū)替過程的變化。由于油水相相對(duì)滲透率比值是含水飽和度的函數(shù)，水相分流量也是含水飽和度的函數(shù)（圖1）。

圖1 水相分流量理論曲線Fig.1 Theoretical curve of water fractional flow

大量室內(nèi)實(shí)驗(yàn)及礦場(chǎng)開發(fā)實(shí)踐表明，水相分流量曲線呈S 型，具有單調(diào)遞增性質(zhì)［7-8］。曲線形態(tài)反映了油水兩相流量在不同驅(qū)替階段的變化過程，初始階段水相分流量較低，表明低含水期水相流速低，地層中流體以油相流動(dòng)為主；進(jìn)入中高含水期后，曲線大幅上升，此時(shí)水相流速逐漸上升，油相流速下降迅速；到了高含水階段，曲線趨近平緩且逐漸接近于1，表明此時(shí)兩相流體以水相流動(dòng)為主，油相流動(dòng)能力極低，驅(qū)替效果越來越差［9-11］。當(dāng)fw=1 時(shí)，地層流體完全為水相單相流，油相不再發(fā)生流動(dòng)，此時(shí)曲線對(duì)應(yīng)的含水飽和度為殘余油飽和度［12-13］。

1.2 邊界條件

殘余油是指被注入水或驅(qū)油劑波及仍驅(qū)替不出來的原油，這部分原油體積占地層孔隙體積的百分?jǐn)?shù)即為殘余油飽和度。對(duì)于室內(nèi)一維小巖心驅(qū)替實(shí)驗(yàn)而言，注入水波及系數(shù)為1，因此，小巖心驅(qū)替的殘余油飽和度是在注入量足夠大甚至達(dá)到無窮大時(shí)出現(xiàn)。由于室內(nèi)實(shí)驗(yàn)注入量有限，利用實(shí)驗(yàn)獲取分流量末端含水飽和度是難以實(shí)現(xiàn)的。

從數(shù)學(xué)角度上，注入量趨近于無限大是存在的，水相分流量理論上可等于1，此時(shí)對(duì)應(yīng)的含水飽和度為水驅(qū)油含水飽和度的極限值。同時(shí)，在注入水即將注入的時(shí)刻，巖心內(nèi)尚未但即將發(fā)生流動(dòng)，此時(shí)水相分流量為0，含水飽和度為原始含水飽和度，由此確定了水相分流量曲線的邊界條件。

初始條件為：

極限條件為：

由于Sorw不確定，直接擬合fw與Sw關(guān)系具有多解性，為提高數(shù)學(xué)模型的精度，定義參數(shù)可采儲(chǔ)量采出程度的物理意義為目前累積采出油量與可動(dòng)油區(qū)間內(nèi)地質(zhì)儲(chǔ)量的比值。其表達(dá)式為：

由于Sw∈[Swi,Sorw]，由（2）式和（3）式可知ηo∈[0,1 ]。將fw與Sw曲線轉(zhuǎn)化成fw與ηo曲線，等同于對(duì)水相分流量曲線進(jìn)行橫向拉伸，縱軸不變，仍可以保持曲線形態(tài)及特征（圖2）。

圖2 fw與ηo關(guān)系曲線Fig.2 Relationship curve between fw and ηo

2 水相分流量曲線模型的建立

2.1 模型基礎(chǔ)結(jié)構(gòu)的確定

通過對(duì)水相分流量曲線邊界條件的分析，將fw∈[0,1 ]，Sw∈[Swi,Sorw]的擬合問題轉(zhuǎn)化成fw∈[0,1 ]，ηo∈[0,1 ]的問題，由于ηo取值范圍明確，對(duì)數(shù)學(xué)模型的選擇也更加嚴(yán)格?？紤]待求曲線的S型特征，結(jié)合曲線定義域和值域均為[0,1 ]的特點(diǎn)，對(duì)目前較為常見的S 型曲線進(jìn)行分析，目前較為常用的S型曲線模型［14-17］如表1所示。

目前大多數(shù)S型曲線的函數(shù)模型由于存在指數(shù)e，導(dǎo)致函數(shù)初始值不能取0，且極限值不能取1，與水相分流量曲線的初始條件和理論極限不符。這些函數(shù)大多可以用于擬合水相分流量曲線的中間區(qū)間，無法用于擬合包含初始值和極限值的全過程分流量曲線。多項(xiàng)式理論上只要階數(shù)足夠高，就可以表征任意復(fù)雜函數(shù)曲線，但考慮多項(xiàng)式的復(fù)雜程度，多項(xiàng)式次冪越高，求解越難，不利于實(shí)際應(yīng)用。有理式具有端值可取、曲線可調(diào)的特征，能夠滿足水相分流量曲線的初始條件和極限條件。同時(shí)，針對(duì)不同類型的水相分流量曲線，有理式能夠通過調(diào)整擬合指數(shù)，滿足中間過程的擬合需要，因此，水相分流量的擬合模型采用有理式形式。

2.2 數(shù)學(xué)模型的建立

參考表1 中S 型曲線有理式的基本形式，令自變量x為可采儲(chǔ)量采出程度，y為水相分流量。由于驅(qū)替極限條件fw=1，令Vmax=kn+1，建立水相分流量曲線的初始模型：

表1 常用S型曲線的數(shù)學(xué)模型及取值范圍統(tǒng)計(jì)結(jié)果Table1 Statistics of common S-curve mathematical model and value range

根據(jù)（5）式的函數(shù)形式，函數(shù)的初始值為ηo=0，fw=0，而極限值滿足條件ηo=1，fw=1，符合水相分流量曲線的邊界條件。數(shù)學(xué)結(jié)構(gòu)上，（5）式包含kn和ηon2個(gè)指數(shù)項(xiàng)，其中kn主要影響曲線初始階段的形態(tài)，kn越小，曲線上升段向左偏移且上升速度越快（圖3）。根據(jù)油藏工程理論，水油黏度比越低（即油水黏度比越高），水相分流量曲線的凸性越突出，表現(xiàn)為油相黏度大，水驅(qū)油見水后含水率上升快，水相分流量曲線上升迅速［18-20］。水油黏度比對(duì)分流量曲線的影響規(guī)律與數(shù)學(xué)模型中kn具有一致性。設(shè)油水黏度比為μR，令，將（5）式改寫為：

實(shí)際油藏開發(fā)中，即使油水黏度比相同，油水運(yùn)動(dòng)規(guī)律也會(huì)因儲(chǔ)層物性及驅(qū)替條件的不同而存在較大差異，水相分流量曲線形態(tài)也存在不同。為保障模型的通用性和擬合精度，對(duì)模型指數(shù)項(xiàng)進(jìn)行區(qū)分，將μR-n改寫為μR-m，進(jìn)一步提高模型在各種條件下的適用性。最終構(gòu)建水相分流量基本函數(shù)式為：

由（7）式可知，改寫后的數(shù)學(xué)模型仍滿足水驅(qū)油過程的邊界條件，主要表現(xiàn)為開始驅(qū)替時(shí)，ηo=0，fw=0；達(dá)到極限驅(qū)替條件時(shí)，巖石孔隙內(nèi)無油相流動(dòng)，采出端只出水，此時(shí)ηo=1，fw=1。對(duì)于特定油藏流體，油水黏度比可認(rèn)為是常數(shù)；而擬合系數(shù)m和n，其主要作用是調(diào)整曲線的S 形態(tài)，使模型在處理不同類型水相分流量曲線均可達(dá)到較好的擬合效果。

將（4）式代入（7）式中，得：

（8）式反映了水相分流量與含水飽和度的關(guān)系，式中未知量包括m，n和Sorw，可通過擬合相對(duì)滲透率測(cè)試中水相分流量測(cè)試結(jié)果，從而得到理論極限驅(qū)替下的殘余油飽和度。

根據(jù)（1）式，油水相相對(duì)滲透率比值的表達(dá)式為：

將（8）式代入（9）式，即可得到油水相對(duì)滲透率比值的全過程曲線。

3 應(yīng)用實(shí)例

圖3 相同n不同k的水相分流量有理式模型計(jì)算結(jié)果Fig.3 Calculation results of rational formula model under same n and different k values

勝利油田中高滲透油藏某砂巖樣品孔隙度為0.34，滲透率為4.67 D，束縛水飽和度為0.18。在50 ℃條件下采用恒壓驅(qū)替，注入水黏度為0.576 3 mPa·s，模擬油黏度為19.75 mPa·s。實(shí)驗(yàn)結(jié)束時(shí)注入量為1 000 PV，含水飽和度為0.728，剩余油飽和度為0.272。油水相相對(duì)滲透率及水相分流量測(cè)試結(jié)果如表2所示。

表2 勝利油田中高滲透油藏某砂巖樣品高注入量相對(duì)滲透率測(cè)試結(jié)果Table2 Test results of relative permeability with high injection volume of sandstone sample from high permeability reservoir of Shengli Oilfield

利用（8）式對(duì)表2中的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行擬合，計(jì)算相對(duì)誤差平均為0.4%，擬合精度極高，計(jì)算結(jié)果具有正確性和可靠性。模型的擬合系數(shù)分別為：m=2.29，n=4.61，Sorw=0.904，進(jìn)而得到理論極限驅(qū)替的殘余油飽和度為0.096，極限采收率為88.26%。

通過（8）式計(jì)算相對(duì)滲透率比值，對(duì)比實(shí)驗(yàn)測(cè)試結(jié)果（圖4）可知，相對(duì)滲透率比值的擬合曲線與實(shí)測(cè)值具有較高的吻合度，表明建立的數(shù)學(xué)模型能夠達(dá)到較好的擬合效果。由于模型考慮了極限條件，得到的相對(duì)滲透率比值曲線具有初始快速下降、中間近似直線段、末端下彎的特征，這一特征與現(xiàn)有油水運(yùn)動(dòng)規(guī)律的認(rèn)識(shí)一致，表明模型具有可靠性。同時(shí)，相對(duì)滲透率比值曲線采用連續(xù)函數(shù)計(jì)算得到，使計(jì)算結(jié)果在數(shù)學(xué)和滲流2 方面的連續(xù)性均得到保障。

圖4 水相分流量及相對(duì)滲透率比值擬合結(jié)果Fig.4 Fitting results between water fractional flow and relative permeability ratio

4 結(jié)論

基于油藏工程理論，利用數(shù)學(xué)方法建立了有理式形式的水相分流量全過程模型。相比以往的指數(shù)e 形式的模型，該模型考慮了水相分流量曲線的端點(diǎn)值，為水驅(qū)油殘余油飽和度的確定以及采收率理論極限的預(yù)測(cè)提供了新的思路和方法。模型所需參數(shù)為相對(duì)滲透率測(cè)試報(bào)告，資料易獲取，應(yīng)用范圍廣泛，具有較高的普及性與實(shí)用性。但目前對(duì)該模型的認(rèn)識(shí)尚淺，根據(jù)油藏工程理論，儲(chǔ)層物性條件、流體性質(zhì)條件及注采條件等因素發(fā)生變化，水相分流量曲線的形態(tài)也會(huì)隨之改變。模型中調(diào)整曲線形態(tài)的參數(shù)為m和n，這2 個(gè)參數(shù)與實(shí)際物性、流體、注采等參數(shù)的關(guān)系及參數(shù)具體物理意義，還需要進(jìn)一步開展研究。

符號(hào)解釋

a,b,c,d,k,r,m,n——擬合系數(shù)；

fw——水相分流量，無量綱；

Kro——油相相對(duì)滲透率，無量綱；

Krw——水相相對(duì)滲透率，無量綱；

qo——油相流速，mL/min；

qw——水相流速，mL/min；

Sorw——?dú)堄嘤惋柡投认聦?duì)應(yīng)的含水飽和度，無量綱；

Sw——含水飽和度，無量綱；

Swi——原始含水飽和度，無量綱；

Vmax——應(yīng)變量的最大值，無量綱；

ηo——可采儲(chǔ)量采出程度，無量綱；

μo——油相黏度，mPa·s；

μR——油水黏度比，無量綱；

μw——水相黏度，mPa·s。