摘要:地下儲氣庫多周期運行盤庫,是研究氣庫運行規(guī)律、分析漏失、進一步提高氣庫運行效率和降低氣庫運行成本的關(guān)鍵環(huán)節(jié)。已建立的未水淹氣藏型地下儲氣庫盤庫數(shù)學(xué)模型卻并不適用于水淹枯竭氣藏型地下儲氣庫,評價結(jié)果明顯高于實際注采氣量。為此,通過深入分析水淹枯竭氣藏型地下儲氣庫氣驅(qū)排水?dāng)U容機理、多周期注采運行過程中氣水兩相宏觀與微觀運動規(guī)律及可能的分布狀態(tài),創(chuàng)新性提出了可動用庫存量的概念,以注采氣量與視地層壓力在一個注采周期內(nèi)滿足定容壓升降方程為評價準(zhǔn)則,建立了水淹枯竭氣藏型地下儲氣庫盤庫模型。實例驗證結(jié)果表明,該模型計算的結(jié)果與實際吻合,可用于水淹枯竭氣藏型地下儲氣庫多周期運行盤庫計算和分析。為儲氣庫運行、管理及調(diào)整提供了依據(jù)。
關(guān)鍵詞:地下儲氣庫;水淹;枯竭氣藏;盤庫;庫存量;模型
我國已建成的用于城市調(diào)峰的地下儲氣庫位于大港板橋油氣區(qū),以水淹枯竭氣藏型為主[1~4]。針對未水淹氣藏型儲氣庫自身特點,以真實庫存量為盤庫出發(fā)點,建立了未水淹氣藏型地下儲氣庫盤庫數(shù)學(xué)模型[5],但該模型用于庫群水淹枯竭氣藏型地下儲氣庫多周期運行動態(tài)評價,結(jié)果卻高于實際2~4倍。未水淹氣藏型地下儲氣庫的氣水關(guān)系簡單,多周期運行過程中,儲存在地下的天然氣能及時有效動用;而庫群水淹枯竭氣藏型地下儲氣庫的儲層物性中等,中孔中滲,非均質(zhì)性強,建庫前經(jīng)歷近30a衰竭式開發(fā),水驅(qū)強度為弱-中等,構(gòu)造頂部氣井水淹停噴,氣水關(guān)系極為復(fù)雜,含氣孔隙中存在封閉氣,侵入水占據(jù)大量孔隙,多周期運行時井間孔隙未能有效波及,導(dǎo)致可動含氣孔隙體積大大減少。因此,非常有必要建立適合水淹枯竭氣藏型地下儲氣庫盤庫的評價模型。
1 水淹枯竭氣藏型儲氣庫的運行機理
受儲層物性及潤濕性影響,水淹枯竭氣藏內(nèi)部存在較多的封閉氣。由于儲層巖石具有親水性,在氣水兩相流動過程中,當(dāng)驅(qū)替壓差不大時,無論是孔隙還是喉道,氣水分布及流動方式主要表現(xiàn)為水包氣,水沿管壁流動,氣體在孔道中央流動,水驅(qū)氣過程形成不同形式的封閉氣,表現(xiàn)為繞流、卡斷、孔隙盲端、角隅以及“H”型孔道形成的封閉氣[6~7]。在氣藏開發(fā)過程中,由于地層水或邊水侵入后占據(jù)了一定的孔隙空間,凝析氣反凝析損失也會占據(jù)一部分孔隙空間,從而減少了可動含氣孔隙體積。建庫后多周期運行過程中氣驅(qū)水淹區(qū)主要對象仍然是以大孔道為主,微細(xì)孔道難以有效驅(qū)替,有效供氣半徑減小,從而降低了注采井網(wǎng)對砂體控制程度,使得部分氣體不能及時動用。為此,提出了可動用庫存量概念(其模型如圖1所示),即注采氣階段壓力波及范圍內(nèi)能有效動用的庫存量。庫存分為可動用庫存和不可動用庫存兩部分,彼此間相互作用和轉(zhuǎn)化,隨著注采井網(wǎng)及注采方式進一步完善,可動用庫存量增加,動用程度也隨之提高。
2 盤庫數(shù)學(xué)模型的建立
在應(yīng)用物質(zhì)平衡方法建立可動用庫存量及可動含氣孔隙體積模型基礎(chǔ)上,聯(lián)立真實氣體狀態(tài)方程,建立了水淹枯竭氣藏型儲氣庫盤庫模型,為盤庫分析奠定堅實的理論基礎(chǔ),同時為國內(nèi)同類型地下儲氣庫建設(shè)和生產(chǎn)運行管理提供科學(xué)依據(jù)。
2.1 可動用庫存量
對于由弱-中等水驅(qū)氣藏改建的地下儲氣庫,忽略巖石和束縛水彈性膨脹作用,假設(shè)注采氣量與視地層壓力在一個注采周期內(nèi)滿足定容壓升降方程,建立了可動用庫存量數(shù)學(xué)模型:
當(dāng)缺乏高壓物性資料時,采用摩爾體積加權(quán)得到混合流體密度,根據(jù)經(jīng)驗公式計算偏差系數(shù),然后迭代求解可動用庫存量[8~10]。
2.2 地下可動含氣孔隙體積
將式(1)可動用庫存量反算到地下儲氣庫地層條件下,得到地下可動含氣孔隙體積,數(shù)學(xué)表達(dá)式為:
2.3 可動庫容量
當(dāng)?shù)叵聝鈳爝\行到上限壓力時,庫內(nèi)可動用的天然氣在地面標(biāo)準(zhǔn)條件下的體積為可動庫容量,數(shù)學(xué)表達(dá)式為:
2.4 可動墊氣量
當(dāng)?shù)叵聝鈳爝\行壓力降低到下限壓力時,庫內(nèi)可動用的天然氣在地面標(biāo)準(zhǔn)條件下的體積為可動墊氣量,數(shù)學(xué)表達(dá)式為:
2.5 工作氣量
當(dāng)?shù)叵聝鈳鞆纳舷迚毫\行到下限壓力時采出的天然氣量,即地下儲氣庫可動庫容量與可動墊氣量之差為工作氣量,數(shù)學(xué)表達(dá)式為:
2.6 總墊氣量
當(dāng)?shù)叵聝鈳爝\行壓力降低到下限壓力時,庫內(nèi)天然氣在地面標(biāo)準(zhǔn)條件下的體積為總墊氣量,即為不可動用庫存量與可動墊氣量之和,數(shù)學(xué)表達(dá)式為:
2.7 總庫容
當(dāng)?shù)叵聝鈳爝\行到上限壓力時,庫內(nèi)天然氣在地面標(biāo)準(zhǔn)條件下的體積為總庫容,即為工作氣量與總墊氣量之和,數(shù)學(xué)表達(dá)式為:
2.8 墊氣損耗量
從下限壓力升至上限壓力時地下儲氣庫的注氣量與從上限壓力降至下限壓力時地下儲氣庫的采氣量之差為墊氣損耗量[11~13],其大小為本周期與上一周期總墊氣量之差,數(shù)學(xué)表達(dá)式為:
3 應(yīng)用實例
3.1 盤庫計算基礎(chǔ)數(shù)據(jù)
國內(nèi)某弱邊水凝析氣藏型地下儲氣庫設(shè)計運行壓力區(qū)間為30.5~13.0MPa,建庫前剩余地質(zhì)儲量為9.88×108m3,地層壓力為11.48MPa,地層溫度為102℃,剩余凝析氣的相對密度為0.7570,凝析油相對密度為0.7500。地下儲氣庫投入運行后,完成了6個完整的注氣和采氣周期,累計注氣24.59×108m3,累計采氣19.87×108m3,凈注氣量為4.75×108m3,回采率為81%,運行動態(tài)數(shù)據(jù)詳見表1。
3.2 盤庫計算結(jié)果
根據(jù)新建的盤庫數(shù)學(xué)模型,采用迭代算法,利用某地下儲氣庫多周期運行注采動態(tài)系列數(shù)據(jù),連續(xù)計算了各周期末的盤庫技術(shù)指標(biāo),見表2。
3.3 盤庫計算結(jié)果分析
根據(jù)盤庫結(jié)果,分析可動用庫存量、可動含氣孔隙體積、工作氣量、墊氣損耗量及損耗率等主要技術(shù)指標(biāo)的變化規(guī)律。
3.3.1可動用庫存量
可動用庫存量逐漸增加,不可動用庫存量總體略有降低,庫存的動用程度隨之增加,其變化示意圖見圖2。經(jīng)過6個周期后,可動用庫存量由2.85×108m3增至11.63×108m3,增加了3.1倍;不可動用庫存量由8.62×108m3降至8.32×108m3,降幅為3.5%;庫存動用程度從25%增至58%,提高了1.3倍,但目前仍有42%庫存不可動用。
3.3.2可動含氣孔隙體積
可動含氣孔隙體積總體上具有明顯上升趨勢,運行初期較快,之后減緩趨穩(wěn),其變化示意圖見圖3。經(jīng)過6個周期后,可動含氣孔隙體積由149×104m3增至456×104m3,增加了2倍。第6周期末庫內(nèi)可動含氣孔隙體積僅為總孔隙體積的43%,未能動用的孔隙體積比例高達(dá)57%。
3.3.3工作氣量
工作氣能力快速提高,工作氣比例總體增加,其變化示意圖見圖4。經(jīng)過6個周期后,工作氣量由1.97×108m3增至6.10×108m3,增加2.1倍。工作氣比例由15.9%增至30.8%,增幅為90%。
3.3.4墊氣損耗量及損耗率
運行初期氣庫快速擴容,墊氣損耗量增加,損耗率較高,之后進入穩(wěn)定擴容階段,墊氣損耗量開始減少,損耗率隨之大幅降低,其變化示意圖見圖5。前3個周期損耗量由0.572×108m3增至1.101×108m3,增幅為92%;隨后3個周期快速降至0.221×108m3,降幅達(dá)80%。墊氣損耗率由35.9%降低到3.6%,降幅達(dá)90%。
3.4 小結(jié)
1) 由于地下儲氣庫在建庫前大面積水淹,受儲層物性及潤濕性影響,孔隙中存在大量封閉氣,建庫后注入氣可能被再次捕集。因此,水淹枯竭氣藏型地下儲氣庫多周期運行過程中部分含氣孔隙不能及時動用,即存在不可動用庫存。
2) 含氣孔隙體積總體呈現(xiàn)上升趨勢。但由于建庫前邊水侵入占據(jù)孔隙空間,多周期運行注采井網(wǎng)控制范圍有限,有效供氣半徑小,導(dǎo)致大量孔隙仍不能有效動用。
3) 由于地下可動含氣孔隙體積上升,使氣庫總體運行效率得以提高,工作氣量及工作氣比例大幅度增加,墊氣損耗量及損耗率總體降低。
4 結(jié)論
1) 針對水淹枯竭氣藏型地下儲氣庫的自身特點,首次提出了可動用庫存量概念,建立了適合于弱 中等水驅(qū)的水淹枯竭氣藏型地下儲氣庫盤庫的數(shù)學(xué)模型,進一步豐富和發(fā)展了儲氣庫盤庫理論。通過國內(nèi)某水淹枯竭氣藏型地下儲氣庫6個周期實例計算驗證了新建盤庫數(shù)學(xué)模型的正確性。
2) 受儲層物性及潤濕性影響,水淹枯竭氣藏內(nèi)部存在較多封閉氣,建庫后氣驅(qū)水淹區(qū)主要對象為儲層相對發(fā)育帶,其他區(qū)域難以有效驅(qū)替,因此,多周期運行過程中部分庫存不可動用。
3) 國內(nèi)某水淹枯竭氣藏型地下儲氣庫多周期運行盤庫分析表明:盤庫技術(shù)指標(biāo)的可靠性在很大程度上取決于注采氣末平均地層壓力,因此,建議在注采轉(zhuǎn)換時留夠平衡期,加強壓力監(jiān)測。
4) 該方法已應(yīng)用于國內(nèi)某庫群多周期運行動態(tài)跟蹤評價中,取得較好效果,為水淹枯竭氣藏型地下儲氣庫調(diào)整方案編制與實施、多周期運行配注配產(chǎn)優(yōu)化提供了科學(xué)依據(jù),同時對提高氣庫運行效率、降低成本具有重要指導(dǎo)意義。
符號說明
Gr(i)為某周期末庫存量,108m3;Grm(i-1)為某周期初可動用庫存量,108m3;Grm max(i)為某周期可動庫容量,108m3;Grm min(i)為某周期可動墊氣量,108m3;Gr work(i)為某周期工作氣量,108m3;Gr min(i)為某周期總墊氣量,108m3;Gr max(i)某周期總庫容量,108m3;QSH(i)為某周期墊氣損耗量,108m3;Vm(i)為某周期可動用含氣孔隙體積,108m3;Q(i)為某周期注采氣量,108m3;p(i)為某周期末地層壓力,MPa;T(i)為某周期末地層溫度,℃;Z(i)為某周期天然氣偏差系數(shù);psc為標(biāo)準(zhǔn)狀態(tài)下的壓力,MPa;Tsc為標(biāo)準(zhǔn)狀態(tài)下的溫度,K;i表示注采周期數(shù),其中奇(偶)數(shù)分別表示注(采)氣周期。
參考文獻
[1] 王起京,張余,劉旭.大張坨地下儲氣庫地質(zhì)動態(tài)及運行效果分析[J].天然氣工業(yè),2003,23(2):89-92.
[2] 楊樹合,何書梅,楊波,等.大張坨地下儲氣庫運行實踐及評價[J].天然氣地球科學(xué),2003,14(5):425-428.
[3] 張幸福,謝廣祿,曾杰,等.大張坨地下儲氣庫運行模式分析[J].天然氣地球科學(xué),2003,14(4):240-244.
[4] 王皆明,姜風(fēng)光.地下儲氣庫注采動態(tài)預(yù)測模型[J].天然氣工業(yè),2009,29(2):108-110.
[5] 王皆明,胡旭健.凝析氣藏型地下儲氣庫多周期運行盤庫方法[J].天然氣工業(yè),2009,29(9):100-102.
[6] WU Haojiang,ZHOU Fangde,WU Yuyuan.Intelligent identification system of flow regime of oil-gas-water multiphase flow [J].International Journal of Multiphase Flow,2001,27(3):459-475.
[7] CIESLINSKI J T,MOSDORF R.Gas bubble dynamics-experiment and fractal analysis[J].International Journal of Heat and Mass Transfer,2005,48(9):1808-1818.
[8] 李登偉,張烈輝,周克明,等.可視化微觀孔隙模型中氣水兩相滲流機理[J].中周石油大學(xué)學(xué)報:自然科學(xué)版,2008,32(3).80-83.
[9] 周克明,李寧,張清秀,等.氣水兩相滲流及封閉氣的形成機理實驗研究[J].天然氣工業(yè),2002,22(增刊):122-125.
[10] 楊繼盛,劉建儀.采氣實用計算[M].北京:石油工業(yè)出版社,1994.
[11] 劉旭,王鳳田,王起京.Q/SY 195.1—2007地下儲氣庫天然氣損耗計算方法[S].北京:石油工業(yè)出版社,2007.
[12] 白雪峰,房愛兵.地下儲氣庫天然氣回采的節(jié)流脫水分析[J].天然氣工業(yè),2006,26(6):121-123.
[13] 譚羽非,林濤.利用地下含水層儲存天然氣應(yīng)考慮的問題[J].天然氣工業(yè),2006,26(6):114-117.
(本文作者:胥洪成 王皆明 李春 中國石油勘探開發(fā)研究院廊坊分院)
贊 賞 分享
您可以選擇一種方式贊助本站
支付寶轉(zhuǎn)賬贊助
微信轉(zhuǎn)賬贊助
- 注解:本資料由會員及群友提供僅供閱讀交流學(xué)習(xí),不得用于商業(yè)用途!