王胤淵，陳小凡，劉峰，楊博義

（西南石油大學(xué)油氣藏地質(zhì)及開發(fā)工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，四川 成都610500）

0 引言

水驅(qū)曲線是水驅(qū)開發(fā)油田固有的特征曲線，不僅可預(yù)測(cè)油田的可采儲(chǔ)量和采收率，還可對(duì)油田可動(dòng)儲(chǔ)量作出評(píng)估，因而在我國(guó)得到了廣泛應(yīng)用［1-9］。然而，油田開發(fā)指標(biāo)與開發(fā)時(shí)間密切相關(guān)，水驅(qū)曲線卻無法預(yù)測(cè)油藏開發(fā)指標(biāo)與開發(fā)時(shí)間的關(guān)系。開發(fā)模型法也是一種重要的油藏工程預(yù)測(cè)方法，該方法以數(shù)理統(tǒng)計(jì)為基礎(chǔ)，通過開發(fā)歷史擬合確定油氣藏的數(shù)學(xué)模型，并對(duì)未來的生產(chǎn)指標(biāo)進(jìn)行預(yù)測(cè)。其中，威布爾（Weibull）預(yù)測(cè)模型不僅可預(yù)測(cè)油氣田產(chǎn)量、累計(jì)產(chǎn)量隨時(shí)間的變化規(guī)律，還可預(yù)測(cè)可采儲(chǔ)量、最高年產(chǎn)量及發(fā)生時(shí)間，但無法預(yù)測(cè)含水率等開發(fā)指標(biāo)［10-11］。

在目前已有的20 多種水驅(qū)曲線中，甲型水驅(qū)曲線的應(yīng)用范圍最廣，精確度最高，本文將其與Weibull 模型相結(jié)合，建立了一種新的預(yù)測(cè)模型，用以預(yù)測(cè)含水率隨時(shí)間的變化規(guī)律［12-17］。

1 模型推導(dǎo)

1.1 Weibull 預(yù)測(cè)模型

Weibull 預(yù)測(cè)模型是基于數(shù)理統(tǒng)計(jì)學(xué)中的Weibull分布，經(jīng)推導(dǎo)轉(zhuǎn)換得出的，其基本關(guān)系式包括：

式中：qo為年產(chǎn)油量，104t；t 為油田開發(fā)時(shí)間，a；Np為累計(jì)產(chǎn)油量，104t；NR為可采儲(chǔ)量，104t；a，b，c 為模型的經(jīng)驗(yàn)常數(shù)。

將式（1）兩邊對(duì)時(shí)間t 求導(dǎo)，并令dqo/dt=0，得到最高年產(chǎn)量發(fā)生時(shí)間tm：

將式（4）代入式（1），得到最高年產(chǎn)量qm：

聯(lián)立式（2）與式（5），得到最高年產(chǎn)量對(duì)應(yīng)的累計(jì)產(chǎn)油量Npm：

1.2 甲型水驅(qū)曲線

甲型水驅(qū)曲線最早是由馬克西莫夫提出的，1978年童憲章將其命名為甲型水驅(qū)曲線，其基本關(guān)系式為

式中：Wp為累計(jì)產(chǎn)水量，104t；A，B 為經(jīng)驗(yàn)常數(shù)。將式（7）兩邊對(duì)時(shí)間t 求導(dǎo)，可得：

式中：Rwo為生產(chǎn)水油比；qw為年產(chǎn)水量，104t。

結(jié)合Rwo與含水率fw的關(guān)系式，進(jìn)一步推導(dǎo)得出：

1.3 模型聯(lián)解

將式（2）代入式（9），即可得到不同開發(fā)時(shí)間的含水率預(yù)測(cè)公式：

為了確定模型常數(shù)a，b，c，將式（1）兩邊取自然對(duì)數(shù)，并令，則

可采用單變量尋優(yōu)算法對(duì)式（11）進(jìn)行求解，即假設(shè)一系列b 值，并繪制相應(yīng)的關(guān)系曲線，對(duì)曲線進(jìn)行線性回歸，最佳擬合直線對(duì)應(yīng)的b 值即為常數(shù)b。再根據(jù)擬合直線的截距α 和斜率β，由確定參數(shù)a，c。為確保模型常數(shù)a，b，c 的值正確可靠，可先利用式（1）、（2）、（10）計(jì)算理論產(chǎn)油量、累計(jì)產(chǎn)油量和含水率，而后將其分別與實(shí)際值進(jìn)行對(duì)比。只有當(dāng)各生產(chǎn)指標(biāo)達(dá)到最佳擬合效果時(shí)，參數(shù)a，b，c 才具有使用價(jià)值。將可用的模型參數(shù)a，b，c 代入式（10），即可進(jìn)行含水率隨開發(fā)時(shí)間變化規(guī)律的預(yù)測(cè)。

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

借用文獻(xiàn)［14］中大慶南二三開發(fā)區(qū)葡一組油田的開發(fā)數(shù)據(jù)，利用Weibull 預(yù)測(cè)模型與甲型水驅(qū)曲線的聯(lián)解法進(jìn)行含水率等生產(chǎn)指標(biāo)的預(yù)測(cè)，并與實(shí)際結(jié)果進(jìn)行對(duì)比，以驗(yàn)證聯(lián)解模型的可靠性與實(shí)用性。

根據(jù)式（7），在半對(duì)數(shù)坐標(biāo)系中繪制葡一組油田的Wp-Np關(guān)系曲線（見圖1）。可以看出，lg Wp-Np呈較好的線性關(guān)系，對(duì)直線段進(jìn)行線性回歸，相關(guān)系數(shù)高達(dá)0.999 4。由此可見，甲型水驅(qū)曲線能較好地反映該油田的水驅(qū)開發(fā)特征。

圖1 葡一組油田lg Wp-Np 關(guān)系

假設(shè)系列b 值（從1 開始試?。罁?jù)式（11）繪制不同b 值條件下的關(guān)系曲線（見圖2）。對(duì)曲線進(jìn)行線性回歸，可以看出，當(dāng)b=1.125 時(shí)曲線的相關(guān)系數(shù)最大，此時(shí)直線的截距α 為2.927 3，斜率β 為0.002 8，相關(guān)系數(shù)為0.997 6。由a-α 及c-β 關(guān)系式，求得a=18.632 4，c=357.142 9。

將a，b，c 的值分別代入式（4）—（6），得到最高年產(chǎn)油量為176.05×104t，對(duì)應(yīng)的開發(fā)時(shí)間和累計(jì)產(chǎn)油量分別為11.79 a 和1 287.2×104t，與實(shí)際年產(chǎn)油量及累計(jì)產(chǎn)油量隨時(shí)間變化曲線的最高點(diǎn)基本吻合（見圖3）。

同時(shí)，將a，b，c 的值分別代入式（1）、式（2）和式（10），可以得到葡一組油田產(chǎn)油量、累計(jì)產(chǎn)油量和含水率的理論預(yù)測(cè)值。將預(yù)測(cè)數(shù)據(jù)與實(shí)際生產(chǎn)數(shù)據(jù)進(jìn)行對(duì)比（見圖3），可以看出，聯(lián)解模型的預(yù)測(cè)結(jié)果總體上與各開發(fā)指標(biāo)隨開發(fā)時(shí)間的變化趨勢(shì)相吻合。其中，在產(chǎn)量上升階段，年產(chǎn)油量的預(yù)測(cè)值與實(shí)際值的誤差相對(duì)較大，但是在產(chǎn)量遞減階段，預(yù)測(cè)精度較高，誤差較小。分析認(rèn)為，這主要是由于水驅(qū)曲線一般用于含水率大于50%的油田開發(fā)階段。由此可見，確定的a，b，c 值是比較可靠的。為此，應(yīng)用聯(lián)解模型對(duì)葡一組油田未來20 a 的開發(fā)生產(chǎn)指標(biāo)進(jìn)行了預(yù)測(cè)（見圖3）。

圖2 不同b 值下的ln （qo /tb）-tb+1 關(guān)系曲線

圖3 葡一組油田生產(chǎn)指標(biāo)預(yù)測(cè)值與實(shí)際值對(duì)比

3 結(jié)束語

將Weibull 預(yù)測(cè)模型與甲型水驅(qū)曲線相結(jié)合，推導(dǎo)得出開發(fā)指標(biāo)隨開發(fā)時(shí)間變化規(guī)律的聯(lián)解預(yù)測(cè)模型，克服了2 種單一方法的局限性，實(shí)現(xiàn)了含水率隨開發(fā)時(shí)間變化規(guī)律的預(yù)測(cè)。從應(yīng)用結(jié)果看，預(yù)測(cè)結(jié)果與實(shí)際生產(chǎn)數(shù)據(jù)基本吻合，但仍存在一定誤差，有待進(jìn)一步研究與改進(jìn)。

［1］李傳亮.油藏工程原理［M］.北京：石油工業(yè)出版社，2005：211-278.

［2］吳宏利，孔自超，李蕾，等.一種研究含水率上升規(guī)律的新方法［J］．?dāng)鄩K油氣田，2001，8（3）：37-38.

［3］張虎俊.預(yù)測(cè)油田產(chǎn)量的新模型及其應(yīng)用［J］．大慶石油地質(zhì)與開發(fā)，1995，14（4）：40-45.

［4］俞啟泰.兩種水驅(qū)特征曲線及其應(yīng)用［J］.大慶石油地質(zhì)與開發(fā)，1997，16（1）：31-34.

［5］曹立迎.底水油藏水平井水淹規(guī)律經(jīng)驗(yàn)?zāi)Ｐ停跩］．?dāng)鄩K油氣田，2012，19（3）：323-326.

［6］王怒濤.新型水驅(qū)特征曲線系列［J］.西南石油大學(xué)學(xué)報(bào)：自然科學(xué)版，2007，29（2）：103-105.

［7］王華，邴紹獻(xiàn)，張海燕，等.特高含水期水驅(qū)油效率計(jì)算新方法［J］．?dāng)鄩K油氣田，2013，20（2）：201-203.

［8］陳元千．水驅(qū)曲線關(guān)系式的推導(dǎo)［J］．石油學(xué)報(bào)，1985，6（2）：69-78．

［9］閻靜華，許尋，杜永波.計(jì)算相滲曲線的新方法：甲型水驅(qū)曲線法［J］．?dāng)鄩K油氣田，2001，8（1）：38-40.

［10］俞啟泰.Weibull增長(zhǎng)曲線預(yù)測(cè)油田開發(fā)指標(biāo)的0.632 12準(zhǔn)數(shù)［J］．石油勘探與開發(fā)，1998，25（2）：69-78．

［11］Wang Jixia，Miao Yu.Bivariate generalized Weibull distribution model ［J］．Journal of Mathematics，2012，32（4）：69-73．

［12］孫來喜，付國(guó)民，馮文光，等.基于Logistic模型及水驅(qū)曲線的一種綜合預(yù)測(cè)模型［J］.礦物巖石，2002，22（4）：97-100．

［13］趙鋒，陳祥光.Γ模型與水驅(qū)曲線的一種聯(lián)解法［J］．特種油氣藏，2004，11（1）：44-46.

［14］陳元千，趙慶飛.Weibull預(yù)測(cè)模型與乙型水驅(qū)曲線的聯(lián)解法［J］．新疆石油地質(zhì)，2000，21（5）：405-407．

［15］魏洪濤，王怒濤.兩種新型水驅(qū)特征曲線系列［J］.西南石油大學(xué)學(xué)報(bào)：自然科學(xué)版，2009，22（6）：343-347.

［16］張光強(qiáng)，朱有志，劉忠云.一種簡(jiǎn)易的水驅(qū)特征曲線的校正方法［J］.斷塊油氣田，1997，4（3）：30-32.

［17］羅洪，俞啟泰.用增長(zhǎng)曲線一體化法預(yù)測(cè)水驅(qū)油田開發(fā)指標(biāo)［J］.斷塊油氣田，2001，8（3）：19-24.