楊炳秀 何東博 王麗娟 孟德偉
1.中國石油勘探與生產(chǎn)公司, 北京 100007;2.中國石油勘探開發(fā)研究院, 北京 100083
?
某氣田氣井壓力前緣到達(dá)滲流邊界時間的判斷
楊炳秀1何東博2王麗娟2孟德偉2
1.中國石油勘探與生產(chǎn)公司,北京100007;2.中國石油勘探開發(fā)研究院,北京100083
致密砂巖氣藏儲層物性差,流體流動阻力大、滲流速度低,壓力前緣經(jīng)過較長時間才能到達(dá)滲流邊界。在進(jìn)行氣井生產(chǎn)指標(biāo)評價時,采用壓力前緣到達(dá)邊界前和到達(dá)邊界后的生產(chǎn)動態(tài)資料所得到的結(jié)果存在較大差異。準(zhǔn)確判斷氣井壓力前緣傳導(dǎo)到邊界的時間對評價氣井生產(chǎn)指標(biāo)尤為重要。以中國某致密砂巖氣田為例,應(yīng)用流線模擬技術(shù)、影響半徑計算公式、產(chǎn)量不穩(wěn)定分析法和氣井動儲量圖版等四種方法進(jìn)行了壓力前緣到達(dá)邊界時間的判斷,并對判斷結(jié)果進(jìn)行了綜合分析,得到該氣田氣井的壓力前緣到達(dá)邊界的時間為200~300 d。該結(jié)果較好地解釋了隨著生產(chǎn)時間的延長,氣井生產(chǎn)指標(biāo)發(fā)生變化的原因,為科學(xué)合理評價氣井生產(chǎn)指標(biāo)提供了指導(dǎo)。
致密砂巖;氣井;壓力前緣;邊界;時間
致密砂巖氣是目前我國開發(fā)規(guī)模最大的非常規(guī)天然氣,在天然氣產(chǎn)量中所占比例較大[1-2]。中國某致密砂巖氣田近年來得到了規(guī)?;行ч_發(fā),已達(dá)到年產(chǎn)250×108m3的生產(chǎn)規(guī)模,成為中國規(guī)模最大的天然氣田[3-4]。
致密砂巖氣藏儲層物性差、非均質(zhì)性強(qiáng)[5],氣體流動阻力大、滲流速度低,氣井的壓力前緣經(jīng)過較長時間才能到達(dá)邊界,儲集體的整體特征才能通過生產(chǎn)動態(tài)體現(xiàn)出來[6]。壓力前緣到達(dá)邊界的時間判斷,對科學(xué)合理地評價氣井的生產(chǎn)指標(biāo)較為重要。本文以中國某致密砂巖氣田為例,對氣井壓力前緣到達(dá)邊界的時間進(jìn)行判斷。
2002年氣田投產(chǎn),2006年開始規(guī)模建產(chǎn),經(jīng)過8年的持續(xù)開發(fā),完成建產(chǎn)任務(wù),年產(chǎn)能達(dá)到250×108m3,進(jìn)入穩(wěn)產(chǎn)階段。由于儲層物性差,非均質(zhì)性強(qiáng),氣井控制儲量和產(chǎn)能較低[9]。直井控制儲量主要集中在1 000×104~3 500×104m3,單井配產(chǎn)1×104m3/d[10-11]。
致密砂巖氣藏的儲層物性差,氣體流動速度慢,同時儲集層非均質(zhì)性強(qiáng),內(nèi)部發(fā)育細(xì)小的夾層,阻礙氣體流動,導(dǎo)致氣體流動前緣到達(dá)滲流邊界的時間較長。在流動前緣到達(dá)滲流邊界之前,隨著生產(chǎn)時間的延長,氣井控制儲量會逐漸增加;當(dāng)流動前緣到達(dá)滲流邊界以后,氣井控制儲量保持不變。所以在流動前緣到達(dá)滲流邊界前后,氣井生產(chǎn)的評價指標(biāo)會有較大差異。下面以某致密砂巖氣田為例,通過流線模擬技術(shù)、影響半徑計算公式、產(chǎn)量不穩(wěn)定分析法和氣井動儲量變化圖版等四種方法計算壓力前緣到達(dá)滲流邊界的時間。
2.1流線模擬技術(shù)
流線模擬技術(shù)在早期流管方法基礎(chǔ)上發(fā)展而來,它可以真實地反映出流體的實際運(yùn)移路線,而在傳統(tǒng)有限差分方法模擬中,流體沿著網(wǎng)格流動,網(wǎng)格尺度越大,流體流動軌跡與真實運(yùn)移路線偏差越大。流線模擬技術(shù)將三維模型轉(zhuǎn)化為一維流線模型,沿著壓力梯度的方向,形成流體流動的流線。對于注入井向生產(chǎn)井追蹤流線,通過流線分布可以清楚地顯示油水驅(qū)替前緣的變化,同樣對于存在滲透阻擋層和流動邊界的流動單元,也可以清楚地看到滲透阻擋層和流動邊界對氣井生產(chǎn)動態(tài)的影響,并量化壓力前緣到達(dá)滲透阻擋層和流動邊界的時間。
本文選用比較成熟的數(shù)值模擬軟件Eclipse的Fron-tsim流線模塊設(shè)計了概念模型,依據(jù)某致密砂巖氣田的儲層參數(shù),設(shè)定了模型基本參數(shù),具體如下:氣水兩相均質(zhì)模型,網(wǎng)格大小5 m×5 m×10 m,網(wǎng)格數(shù)240×240×1,孔隙度0.09,滲透率0.07×10-3μm2,含氣飽和度0.7,束縛水飽和度0.3。模型中間部署1口氣井,氣井采用定井底流壓生產(chǎn),模擬結(jié)果見圖1和表1。
從圖1可以看出,流線直觀地顯示了壓力前緣的變化過程,也將邊界對流體流動的影響可視化。根據(jù)表1模擬統(tǒng)計結(jié)果,對滲透率為0.07×10-3μm2,半徑為300~400 m的均質(zhì)儲層,壓力前緣完全到達(dá)邊界的時間為50~200 d。由于模型是理想的均質(zhì)模型,模型計算結(jié)果比實際儲層的壓力前緣到邊界的時間短。所以某氣田氣井的壓力前緣到達(dá)邊界的時間大于區(qū)間50~200 d。
a) 75 d
b) 320 d
c) 610 d
表1應(yīng)用流線模擬技術(shù)計算的影響半徑與開井時間的關(guān)系表
開井時間t/d影響半徑r/m2.65010.310035.820056.3300198.7400298.3500395.8600
2.2影響半徑公式計算
對于1口氣井生產(chǎn)后,井底壓力開始下降,壓力降從井底逐漸向周圍擴(kuò)展,形成一個壓降漏斗,隨著時間的推移,壓降漏斗不斷擴(kuò)大。根據(jù)影響半徑的定義,在t時刻,壓降漏斗的邊界擴(kuò)展到半徑r的位置,半徑r以內(nèi)的地層受到生產(chǎn)井的擾動,壓力發(fā)生了變化,半徑r以外的地層壓力未發(fā)生變化[12]。
影響半徑的公式:
依據(jù)氣田的儲層物性,K=0.07×10-3μm2,μ=0.023 mPa·s,φ=0.09,Ct=0.02 MPa-1時,代入公式,得到影響半徑與時間的關(guān)系(表2)。氣田有效砂體的大小主要為600 m×800 m,根據(jù)表2,壓力傳導(dǎo)到邊界的時間大致在85~273 d。
表2應(yīng)用影響半徑公式計算的影響半徑與開井時間的關(guān)系表
開井時間t/d影響半徑r/m4.35017.110068.520085.2300273.8400427.8500616.1600
2.3產(chǎn)量不穩(wěn)定分析法
產(chǎn)量不穩(wěn)定分析是基于試井理論,利用動態(tài)資料評價氣井的控制儲量和泄氣面積。該方法考慮裂縫長度、表皮系數(shù)、滲流邊界等一系列參數(shù)建立解析模型,利用井的生產(chǎn)動態(tài)歷史數(shù)據(jù)和儲集層的基本地質(zhì)參數(shù)進(jìn)行擬合,使模型計算結(jié)果與井的實際生產(chǎn)動態(tài)和動儲量一致,進(jìn)而可以確定氣井的泄氣半徑和控制儲量。同時,依據(jù)擬合模型,可預(yù)測氣井今后的生產(chǎn)動態(tài),另外根據(jù)擬合圖版,可以間接得到氣體流動進(jìn)入邊界控制流的時間。
圖2 氣井日產(chǎn)量Blasingame典型曲線擬合圖
圖3 物質(zhì)平衡時間
氣井開始投產(chǎn)時,壓降漏斗的范圍不斷擴(kuò)大,當(dāng)壓力前緣到達(dá)邊界后,流動狀態(tài)逐漸轉(zhuǎn)化成邊界控制流。故壓力前緣到達(dá)邊界的時間會稍早于邊界控制流的時間。故126口氣井的壓力前緣到達(dá)邊界的時間稍早于時間段150~300 d。
2.4氣井動儲量圖版
為了分析氣井控制儲量隨時間的變化,選擇230口生產(chǎn)較穩(wěn)定的典型氣井,分別應(yīng)用30、60、90、120、150、200、300、400、600、1 000、1 500 d的生產(chǎn)動態(tài),選用RTA軟件進(jìn)行擬合分析,擬合不同生產(chǎn)時間的歷史動態(tài),得到不同生產(chǎn)時間的控制儲量。
氣井投產(chǎn)之后,壓力逐漸傳播到邊界,在壓力前緣到達(dá)邊界之前,氣井的動儲量逐漸增加,當(dāng)壓力傳播到邊界后,流體進(jìn)入邊界控制流階段,氣井的控制儲量保持不變[14]。從圖4可以看出,在310 d以前,隨著生產(chǎn)時間的延長,計算的動儲量逐漸增加;超過310 d之后,動儲量基本保持不變。圖4也反映了壓力傳導(dǎo)過程中氣井動儲量的變化,所以從圖4可以判斷,壓力前緣到達(dá)邊界的時間在310 d左右。
圖4 應(yīng)用不同生產(chǎn)時間的生產(chǎn)動態(tài)計算的動儲量
1)應(yīng)用流線模擬技術(shù)、影響半徑計算公式、產(chǎn)量不穩(wěn)定分析法和氣井動儲量圖版等四種方法計算,結(jié)果表明,某氣田氣井投產(chǎn)后,壓力前緣到達(dá)邊界的時間大致為200~300 d。
2)氣井生產(chǎn)指標(biāo)評價應(yīng)選擇投產(chǎn)時間超過300 d的氣井進(jìn)行,否則生產(chǎn)動態(tài)評價的生產(chǎn)指標(biāo)會偏低。
3)四種評價方法中流線模擬技術(shù)得出的氣體到達(dá)滲流邊界的時間偏??;影響半徑計算公式法得出的氣體到達(dá)滲流邊界的時間稍偏??;產(chǎn)量不穩(wěn)定分析法和氣井動儲量模板分析法,評價結(jié)果更具客觀性和科學(xué)性。
[1] 黃鑫,董秀成,肖春躍,等.非常規(guī)油氣勘探開發(fā)現(xiàn)狀及發(fā)展前景[J].天然氣與石油,2012,30(6):38-41.
Huang Xin, Dong Xiucheng, Xiao Chunyue, et al. Situation and Development Prospects of Unconventional Oil and Gas Exploration and Development [J]. Natural Gas and Oil, 2012, 30 (6): 38-41.
[2] 徐博.2020年前中國多氣源供應(yīng)格局展望[J].天然氣工業(yè),2012,32(8):53-55.
Xu Bo. Prospect of Multiple Sources of Gas Supplies in China Before 2020 [J]. Natural Gas Industry, 2012, 32 (8): 53-55.
[3] 李海平,賈愛林,何東博,等.中國石油的天然氣開發(fā)技術(shù)進(jìn)展及展望[J].天然氣工業(yè),2010,30(1):5-7.
Li Haiping, Jia Ailin, He Dongbo, et al. Technical Progress and Outlook of Natural Gas Development for the PetroChina [J]. Natural Gas Industry, 2010, 30 (1): 5-7.
[4] 張磊,張文強(qiáng),李晶,等.蘇6與蘇36-11區(qū)塊水平井井位優(yōu)選[J].天然氣與石油,2013,31(3):61-65.
Zhang Lei, Zhang Wenqiang, Li Jing, et al. Optimization of Horizontal Well Position in Su 6 and Su 36-11 [J]. Natural Gas and Oil, 2013, 31 (3): 61-65.
[5] 李建忠,郭斌程,鄭民,等.中國致密砂巖氣主要類型、地質(zhì)特征與資源潛力[J].天然氣地球科學(xué),2012,23(6):607-615.
Li Jianzhong, Guo Bincheng, Zheng Min, et al. Main Types, Geological Features and Resource Potential of Tight Sandstone Gas in China [J]. Natural Gas Geoscience, 2012, 23 (6): 607-615.
[6] 王記俊,廖新武,趙秀娟,等.注采井網(wǎng)對裂縫形態(tài)及壓裂規(guī)模的影響[J].天然氣與石油,2014,32(5):49-51.
Wang Jijun, Liao Xinwu, Zhao Xiujuan, et al. Effect of Flooding Patterns on Fracture Morphology and Size of the Fracture [J]. Natural Gas and Oil, 2014, 32 (5): 49-51.
[7] 王少飛,安文宏,陳鵬,等.蘇里格氣田致密氣藏特征與開發(fā)技術(shù)[J].天然氣地球科學(xué),2013,24(1):138-145.
Wang Shaofei, An Wenhong, Chen Peng, et al. Characteristic and Development Techniques of Sulige Tight Gas Pool [J]. Natural Gas Geoscience, 2013, 24 (1): 138-145.
[8] 余淑明,劉艷俠,武力超,等.低滲透氣藏水平井開發(fā)技術(shù)難點及攻關(guān)建議——以鄂爾多斯盆地為例[J].天然氣工業(yè),2013,33(1):65-69.
Yu Shuming, Liu Yanxia, Wu Lichao, et al. Technical Difficulties and Proposed Countermeasures in Drilling Horizontal Wells in Low-permeability Reservoirs: A Case Study from the Ordos Basin [J]. Natural Gas Industry, 2013, 33 (1): 65-69.
[9] 胡俊坤,李曉平,肖強(qiáng),等.利用生產(chǎn)動態(tài)資料確定氣井產(chǎn)能方程新方法[J].天然氣地球科學(xué),2013,24(5):1027-1031.
Hu Junkun, Li Xiaoping, Xiao Qiang, et al. A New Method of Using Dynamic Data to Determine Deliverability Equation of Gas Well [J]. Natural Gas Geoscience, 2013, 24 (5): 1027-1031.
[10] 王軍磊,賈愛林,何東博,等.致密氣藏分段壓裂水平井產(chǎn)量遞減規(guī)律及影響因素[J].天然氣地球科學(xué),2014,25(2):278-285.
Wang Junlei, Jia Ailin, He Dongbo, et al. Rate Decline of Multiple Fractured Horizontal Well and Influence Factors on Productivity in Tight Gas Reservoirs [J]. Natural Gas Geoscience, 2014, 25 (2): 278-285.
[11] 賈成業(yè),姬鵬程,賈愛林,等.低滲透砂巖氣藏開發(fā)指標(biāo)數(shù)值模擬預(yù)測[J].西南石油大學(xué)學(xué)報(自然科學(xué)版),2010,32(5):100-104.
Jia Chengye, Ji Pengcheng, Jia Ailin, et al. Numerical Simulation Prediction of the Developmental Indices of Low Permeable Sandstone Gas Reservoir [J]. Journal of Southwest Petroleum University (Science & Technology Edition), 2010, 32 (5): 100-104.
[12] 莊惠農(nóng).氣藏動態(tài)描述和試井[M].北京:石油工業(yè)出版社,2004:136-141.
Zhuang Huinong. The Dynamic Description and Well Test in Gas Reservoir [M]. Beijing: Petroleum Industry Press, 2004: 136-141.
[13] 孫賀東.油氣井現(xiàn)代產(chǎn)量遞減分析方法及應(yīng)用[M].北京:石油工業(yè)出版社,2013:18-66.
Sun Hedong. Advanced Production Decline Analysis and Application [M]. Beijing:Petroleum Industry Press, 2013: 18-66.
[14] 公言杰,柳少波,姜林,等.致密砂巖氣非達(dá)西滲流規(guī)律與機(jī)制實驗研究——以四川盆地須家河組為例[J].天然氣地球科學(xué),2014,25(6):804-809.
Gong Yanjie, Liu Shaobo, Jiang Lin, et al. Experimental Study of Seepage Characteristic and Mechanism in Tight Gas Sands: A Case from Xujiahe Reservoir of Sichuan Basin [J]. Natural Gas Geoscience, 2014, 25 (6): 804-809.
2015-11-07
國家科技重大專項“大型油氣田及煤層氣開發(fā)”(2011 ZX 0515)
楊炳秀(1966-),女,山東招遠(yuǎn)人,高級工程師,博士,主要從事天然氣開發(fā)生產(chǎn)管理工作。
10.3969/j.issn.1006-5539.2016.02.009