謝 姍,焦 揚(yáng),艾慶琳,伍 勇,袁繼明
動(dòng)儲(chǔ)量是指氣藏中實(shí)際參與滲流的那部分地質(zhì)儲(chǔ)量,是確定氣井合理產(chǎn)能和井網(wǎng)密度的重要依據(jù),是編制、調(diào)整開(kāi)發(fā)技術(shù)政策的物質(zhì)基礎(chǔ)[1]。最常用的動(dòng)儲(chǔ)量評(píng)價(jià)方法是利用封閉氣藏的物質(zhì)平衡方程,即視壓力差與累計(jì)產(chǎn)氣量存在線性關(guān)系,當(dāng)視壓力差為0時(shí),其對(duì)應(yīng)的累計(jì)產(chǎn)氣量為動(dòng)儲(chǔ)量(圖1藍(lán)色線)。與常規(guī)氣井動(dòng)儲(chǔ)量評(píng)價(jià)相比,產(chǎn)水氣井會(huì)受到水侵的影響,井底獲得外來(lái)壓力補(bǔ)給,使得視壓力差與累計(jì)產(chǎn)氣量曲線上翹偏離直線段(圖 1紅色曲線),導(dǎo)致不可直接利用早期直線段的延長(zhǎng)段評(píng)價(jià)氣井動(dòng)儲(chǔ)量[2]。
隨著水驅(qū)氣藏相關(guān)理論的不斷發(fā)展,地下水體的影響即水侵量的準(zhǔn)確計(jì)算被認(rèn)為是評(píng)價(jià)此類氣藏動(dòng)儲(chǔ)量的關(guān)鍵。目前水侵計(jì)算模型已相對(duì)成熟(表1),主要有穩(wěn)態(tài)模型、修正的穩(wěn)態(tài)模型、不穩(wěn)態(tài)模型等,但這些模型由于需要不確定性的地下水體參數(shù)、計(jì)算繁復(fù)等原因?qū)е缕鋺?yīng)用受到很大限制[3-5]。需要說(shuō)明的是,其中相對(duì)應(yīng)用廣泛的Fetkovich水侵模型雖不依賴于地下水體形狀、性質(zhì),但忽略了水侵開(kāi)始的瞬間時(shí)間段,導(dǎo)致計(jì)算的水侵量偏小[6-8]。
長(zhǎng)慶油田M區(qū)氣井普遍面臨產(chǎn)水情形,研究一種不依賴于水體形狀、性質(zhì),且簡(jiǎn)單、準(zhǔn)確、便于現(xiàn)場(chǎng)應(yīng)用的動(dòng)儲(chǔ)量求取方法十分必要。
以水侵封閉氣藏物質(zhì)平衡原理為基礎(chǔ),結(jié)合M區(qū)的產(chǎn)水特征,研究水侵氣藏的動(dòng)態(tài)特征[9,10]:當(dāng)氣藏流體被開(kāi)采出來(lái)之后,氣藏壓力下降,在氣藏和周圍的水體之間產(chǎn)生壓差,水體遵循孔隙介質(zhì)中流體流動(dòng)的基本規(guī)律,通過(guò)原始?xì)馑缑媲秩耄划?dāng)氣藏受到水侵后,流體相由原來(lái)的氣相變?yōu)闅狻⒁簝上唷?/p>
圖 2表現(xiàn)的就是氣藏前后條件發(fā)生改變的情形,反映在壓降曲線圖1中則是:紅色曲線反映氣藏受到水侵的生產(chǎn)特征,藍(lán)色直線是封閉氣藏(無(wú)水侵)的生產(chǎn)特征,而紅色曲線向上偏轉(zhuǎn)出現(xiàn)在紅藍(lán)曲線兩者之間的垂向距離(視壓力差pPΔ),則反映了水侵的大小。
圖1 受到水侵影響氣井視壓力差和累計(jì)產(chǎn)氣量關(guān)系
表1 水侵計(jì)算模型適應(yīng)性分析
圖2 水侵封閉氣藏物質(zhì)平衡示意圖
若已知不同生產(chǎn)時(shí)間下壓降曲線的視壓力差(ppΔ),便可計(jì)算得到該時(shí)間下的水侵量,再利用這個(gè)水侵量便可進(jìn)一步校正氣井動(dòng)儲(chǔ)量。
由圖2水驅(qū)氣藏物質(zhì)平衡原理可定義水侵體積系數(shù)為:
可知水侵量為:
對(duì)于產(chǎn)水氣藏,若考慮其為正常壓力系統(tǒng),并忽略束縛水的膨脹作用和巖石壓縮引起的膨脹,其物質(zhì)平衡方程可寫(xiě)為[11]:
已知封閉氣藏物質(zhì)平衡方程為:
將水侵體積系數(shù)引入式(3),并與式(4)相減,可得:
計(jì)算時(shí),利用式(5),每計(jì)算一個(gè)視壓力差,就可獲得一個(gè)水侵系數(shù)ω,將封閉氣藏物質(zhì)平衡得到的動(dòng)儲(chǔ)量代入式(2)得到該時(shí)間下的水侵量eW ,再將計(jì)算得到的水侵量eW 代入式(3)可校正水侵后的動(dòng)儲(chǔ)量。該方法僅需氣藏的生產(chǎn)動(dòng)態(tài)數(shù)據(jù)和少量測(cè)壓數(shù)據(jù),可達(dá)到預(yù)測(cè)早期水侵下動(dòng)儲(chǔ)量的目的。
X1井是M氣區(qū)的一口典型產(chǎn)水氣井,目前日產(chǎn)氣0.5×104m3,日產(chǎn)水13.6 m3。該井生產(chǎn)時(shí)間較長(zhǎng),具備相對(duì)豐富的測(cè)壓數(shù)據(jù)。追蹤評(píng)價(jià)壓降曲線發(fā)現(xiàn),生產(chǎn)933天壓降曲線發(fā)生上翹偏轉(zhuǎn)(圖3),結(jié)合生產(chǎn)實(shí)際表明X1井發(fā)生水侵。
通過(guò)視壓力差和 Fetkovich水侵模型兩種方法對(duì)水侵量進(jìn)行計(jì)算,并校正氣井水侵后動(dòng)儲(chǔ)量(表2),結(jié)果表明:①隨著水侵量的增加,曲線偏轉(zhuǎn)越加嚴(yán)重,動(dòng)儲(chǔ)量呈減小趨勢(shì);②Fetkovich水侵模型至少需要3組壓力測(cè)試數(shù)據(jù)才可進(jìn)行擬合計(jì)算,而視壓力差則可單獨(dú)計(jì)算每一個(gè)水侵時(shí)間下的水侵量,且計(jì)算結(jié)果在水侵早期與利用全部測(cè)壓數(shù)據(jù)的Fetkovich方法結(jié)果接近,表明新方法可在早期預(yù)測(cè)水侵量;③視壓力差計(jì)算得到的目前動(dòng)儲(chǔ)量為 1.4×108m3,是未發(fā)生水侵情形動(dòng)儲(chǔ)量的82%左右,利用Fetkovich計(jì)算水侵量校正動(dòng)儲(chǔ)量為1.9×108m3,大于未發(fā)生偏轉(zhuǎn)的測(cè)壓數(shù)據(jù)計(jì)算得到的動(dòng)儲(chǔ)量,認(rèn)為結(jié)果不可靠。出現(xiàn)該錯(cuò)誤的原因一是Fetkovich模型忽略了開(kāi)始瞬間時(shí)間段,計(jì)算的水侵量偏小,二是Fetkovich模型需要大量測(cè)壓數(shù)據(jù)擬合水體形狀、性質(zhì)來(lái)避免多解性,而第3 848天水侵量準(zhǔn)確需要后續(xù)更多的測(cè)壓數(shù)據(jù)才可準(zhǔn)確校正。
表2 視壓力差和Fetkovich水侵模型計(jì)算X1井水侵量與動(dòng)儲(chǔ)量結(jié)果
圖3 X1井壓降曲線
(1)利用視壓力差法計(jì)算動(dòng)儲(chǔ)量不依賴水體形狀、性質(zhì),僅需氣藏的生產(chǎn)動(dòng)態(tài)數(shù)據(jù)和少量測(cè)壓數(shù)據(jù),即達(dá)到預(yù)測(cè)早期水侵下的動(dòng)儲(chǔ)量的目的,具有快速簡(jiǎn)便、相對(duì)可靠的特點(diǎn)。
(2)新方法的提出有力支撐了產(chǎn)水氣井動(dòng)儲(chǔ)量評(píng)價(jià)技術(shù),為加密井部署、工作制度優(yōu)化等方面夯實(shí)了基礎(chǔ)。
符號(hào)注釋
Gp為累計(jì)產(chǎn)氣量,108m3;ppΔ為壓降差,MPa; ω為水侵體積系數(shù);eW 為水侵量,104m3;pW 為累計(jì)產(chǎn)水量,104m3;wB 為地層水體積系數(shù);G為氣井動(dòng)儲(chǔ)量,108m3;giB為原始?xì)獠貕毫ο碌臍怏w體積系數(shù);P,iP為分別為氣藏和原始?xì)獠貕毫Γ琈Pa;Z,iZ為分別為氣藏和原始?xì)獠貕毫ο碌臍怏w偏差因子;T,iT為分別為氣藏和原始條件下地層溫度,K;scP 為地面標(biāo)況下的壓力,MPa。
參考文獻(xiàn)
[1] 王東旭,王鴻章,李躍剛.長(zhǎng)慶氣田難采儲(chǔ)量動(dòng)用程度評(píng)價(jià)[J].天然氣工業(yè),2000,20(5):64–66.
[2] 李傳亮.氣藏水侵量的計(jì)算方法研究[J].新疆石油地質(zhì),2003,24(5):430–431.
[3] SCHITHUIS R J.Active oil and reservoir energy[R].SPE 936033,1936:33–52.
[4] VAN EVERDINGEN A F,HURST W.The application of the Laplace trans-formation to flow problems in reservoirs[R].SPE 949305,1949:305–324.
[5] HURST W.Water influx into a reservoir and its application to the equation of volumetric balance[R].SPE 943057,1943:57–72.
[6] FETKOVICH M J.A simplified approach to water influx calculations finite aquifer systems[J].Journal of Petroleum Technology,1971,23(7):814–828.
[7] SARFRAZ A,okhio,Djebbar Tiab.Establishing inflow performance relationship (IPR) for gas condensate wells[C].Paper SPE 75503 presented at the 2002 SPE Gas Technology Symposium,Alberta,April:20–26.
[8] 塔雷克·艾哈邁德,保羅·麥金尼.現(xiàn)代油藏工程[M].北京:石油工業(yè)出版社,2012:307–310.
[9] 葛家理.油氣層滲流力學(xué)[M].北京:石油工業(yè)出版社,1982:29–31.
[10] 陳元千.氣田天然水侵的判斷方法[J].石油勘探與開(kāi)發(fā),1978,5(3):51–57.
[11] 王怒濤,黃炳光,張崇軍,等.水驅(qū)氣藏動(dòng)態(tài)儲(chǔ)量及水侵量計(jì)算新方法[J].西南石油學(xué)院學(xué)報(bào),2000,22(4):26–32.