氣藏動態(tài)儲量計算中的幾個關鍵參數(shù)探討

摘 要

摘要：傳統(tǒng)的計算方法以物質平衡公式為基礎，利用關井測壓數(shù)據(jù)來計算動態(tài)儲量。現(xiàn)今以“氣井生產動態(tài)曲線特征圖版擬合法”為主的動態(tài)計算儲量方法，特點是利用流動壓力
摘要：傳統(tǒng)的計算方法以物質平衡公式為基礎，利用關井測壓數(shù)據(jù)來計算動態(tài)儲量?，F(xiàn)今以“氣井生產動態(tài)曲線特征圖版擬合法”為主的動態(tài)計算儲量方法，特點是利用流動壓力計算氣井動態(tài)儲量。在對比分析這兩種方法的基礎上，結合國內外已開發(fā)氣田動態(tài)儲量計算實例，研究了關井測壓數(shù)據(jù)和巖石有效壓縮系數(shù)對動態(tài)儲量計算的影響。結果表明，在沒有關井測壓資料的情況下，利用現(xiàn)今的方法可以計算動態(tài)儲量，并能夠準確擬合生產歷史，但會使動態(tài)儲量計算結果存在很大的不確定性。因此關井測壓資料對動態(tài)儲量的準確計算十分重要。對于異常高壓氣藏，儲層巖石有效壓縮系數(shù)對動態(tài)儲量計算結果的影響不可忽視，巖石有效壓縮系數(shù)取值越小，動態(tài)儲量計算結果越大，在開發(fā)早期計算動態(tài)儲量時應該確定可靠的巖石有效壓縮系數(shù)值。
關鍵詞：異常高壓氣藏；儲量；動態(tài)；分析；計算；測試；巖石；壓縮系數(shù)
1 氣藏動態(tài)儲量計算方法概述
    針對不同的氣藏類型，以及不同的開發(fā)階段，都有與之相對應的動態(tài)儲量計算方法。其中最常用的也是最傳統(tǒng)的方法就是以物質平衡為基礎的計算方法。近期又出現(xiàn)了以“氣井生產動態(tài)曲線特征圖版擬合法”為主的動態(tài)計算儲量方法。
1.1 物質平衡法
根據(jù)氣藏物質平衡通式：
 
上式經(jīng)過變化后得到[1]：
 
Ep為變容系數(shù)，表示巖石和束縛水彈性膨脹所提供的能量，表達式為：
 
   定容封閉氣藏的壓降曲線法、水驅氣藏的視地質儲量法(Havlena-Odeh法)、曲線擬合法以及異常高壓氣藏的視地層壓力校正法[2]、Hammmerlindl修正系數(shù)法[3]、求解圖法(Roach、Poston等)、Becerra-Arteaga壓力圖法等，都是基于上述物質平衡基礎上得來的。物質平衡法計算儲量的關鍵是能夠獲取氣藏(井)的平均地層壓力。
1.2 氣井生產動態(tài)曲線特征圖版擬合法[4～6]
    “氣井生產動態(tài)曲線特征圖版擬合法”又稱“產量不穩(wěn)定分析(rate transient analysis，RTA)”、“高級遞減分析(advanced decline analysis)”，是以氣井生產動態(tài)(產量、流動壓力等)為基礎，引入不穩(wěn)定試井分析的基本思想，對傳統(tǒng)的產量遞減法進行了改進，建立了氣井生產曲線特征圖版，把氣井早期的不穩(wěn)定流動段和后期的邊界流動段結合起來，通過特征圖版擬合計算儲層物性、表皮系數(shù)、井控半徑、動態(tài)儲量。主要分析特征圖版包括Arps圖版、Fetkovich圖版、Balasingame圖版、Agwarl-Gardner圖版、NPI(normalized pressure integrate)圖版、Transient圖版、流動物質平衡(flowing material balance)等。
    氣井生產動態(tài)曲線特征圖版擬合法以氣井流壓和產量為基礎，除原始地層壓力外，不需要關井測壓數(shù)據(jù)。
2 關于動態(tài)儲量計算的幾個問題探討
2.1 準確求取地層壓力是動態(tài)法儲量計算的關鍵
    壓力是氣藏能量大小的直接體現(xiàn)。最準確的也是最傳統(tǒng)的確定方法是進行全氣藏關井測壓，然后計算氣藏平均壓力。氣田開發(fā)歷史比較長的川渝氣區(qū)，基本上每個氣田在生產過程中都要安排數(shù)次全氣藏關井測壓求動態(tài)儲量。但隨著投入開發(fā)氣藏類型日趨復雜化，這一監(jiān)測方式也受到挑戰(zhàn)，比如對于主力氣藏由于供氣需要，全氣藏關井受到限制；對于異常高壓氣藏以及高含硫氣藏，由于井下監(jiān)測風險大，可能無法像常規(guī)氣井那樣進行井下測靜壓；還有就是低滲透氣藏，需要長時間的關井地層壓力才能恢復，這些都會給氣藏壓力分析帶來影響。
    近期出現(xiàn)的“氣井生產動態(tài)曲線特征圖版擬合法”可以在不需要關井測壓的情況下計算氣井動態(tài)儲量，為復雜類型氣藏動態(tài)分析提供了有效手段，該方法的應用也逐漸受到技術人員的關注。值得強調的是，通過關井進行靜壓監(jiān)測是不可替代的。下面列舉墨西哥灣-氣井(G-111井，原始地層壓力122.7MPa，儲層溫度177℃)的實例來說明靜壓資料獲取對動態(tài)儲量計算的重要性。利用“氣井生產動態(tài)曲線特征圖版擬合法”中的“流動物質平衡”計算井控儲量5.35×108m3(圖1)。根據(jù)這一儲量計算結果，建立單井模型擬合生產歷史，井底流壓擬合程度高(圖2)，如果沒有關井壓力作參照，可以認為流動物質平衡計算的井控儲量是可靠的。但如果參照關井壓力數(shù)據(jù)，就會發(fā)現(xiàn)模型計值靜壓值偏低，說明井控儲量計算存在誤差，根據(jù)這一情況，對模型中井控儲量進行了調整，并考慮了滲透率的應力敏感性，完全擬合后，得到井控儲量為8.5×108m3。由此可見靜壓數(shù)據(jù)的獲取對動態(tài)儲量計算的重要性。
 

2.2 異常高壓氣藏巖石有效壓縮系數(shù)對動態(tài)儲量計算影響
    對異常高壓氣藏，由于儲層的欠壓實作用，使得巖石有效壓縮系數(shù)(Cf)比常壓氣藏高一個數(shù)量級，在4×10-6×10-3MPa-1，而表征氣體彈性膨脹能力的壓縮系數(shù)Cg值遠低于常壓情況，其數(shù)量級也在10-3MPa-1，也就是說在高壓階段Cg值與Cf值數(shù)量級相當，因此在高壓階段不能忽略巖石孔隙壓縮所提供的彈性膨脹能量，這一特點使得高壓氣藏的壓降曲線出現(xiàn)兩段式(見圖3)，直接用高壓階段的壓降曲線計算動態(tài)儲量會導致結果偏高。針對這一情況，異常高壓氣藏物質平衡法考慮到了巖石有效壓縮系數(shù)Cf影響，對壓降曲線進行了修正。圖4就是A氣藏壓降曲線圖，根據(jù)壓降曲線回歸，傳統(tǒng)的p/Z-Gp曲線計算動態(tài)儲量為3616×108m3，是容積法儲量的1.3倍，考慮到巖石有效壓縮系數(shù)的情況下，經(jīng)過修正后的壓降曲線計算動態(tài)儲量為2412×108m3，與容積法儲量基本符合。由此可見對異常高壓氣藏，巖石彈性膨脹為氣體的采出提供了不可忽視的能量。
 

    實驗分析表明，對于同一個氣藏，Cf值與其他物性參數(shù)一樣，也同樣具有非均質性。圖5為A氣藏巖心分析結果，從中可以看出，Cf與儲層滲透率有關，高滲儲層Cf值低，低滲儲層Cf值高，對于該氣田不同滲透率巖樣，Cf的變化范圍為2×10-3～10×10-3MPa-1。針對不同的Cf值，動態(tài)儲量計算結果也不同，Cf值越大，動態(tài)儲量計算結果越?。籆f值越小，動態(tài)儲量計算結果越大。
 

   在數(shù)值模擬中，巖石壓縮系數(shù)作為已知參數(shù)輸入到模型中，其取值會給歷史擬合過程中儲量的擬合帶來影響，在其他參數(shù)相向、流壓和關井壓力完全擬合的情況下，不同的Cf值對應不同擬合儲量(見圖6)，從另一個角度來看，對于一個未投入開發(fā)異常高壓氣藏，以容積法儲量為基礎進行生產動態(tài)預測時，不同的Cf值會得到不同的早期壓力變化趨勢。由此可見，對于異常高壓氣藏，獲得可靠的巖石壓縮系數(shù)值對早期動態(tài)儲量計算和生產特征認識非常重要。
 

2.3 動態(tài)儲量發(fā)生變化的情況
   由于動態(tài)儲量代表的是參與流動的地質儲量。在有些情況下，氣藏中的流動供給發(fā)生了變化，也會導致動態(tài)法計算的儲量發(fā)生變化，比較常見的是非均質性很強的低滲氣藏和具有邊底水的氣藏。
    對于非均質性比較強的低滲氣藏，儲層中的流動存在著啟動壓力。在開采初期，參與流動的主要是滲透率相對較高的或縫洞比較發(fā)育的儲層中的氣體，隨著天然氣的采出，儲層中的壓降漏斗進一步增大，低滲儲層中的天然氣也開始參與流動。這類氣藏的壓降曲線(p/Z-Gp曲線)具有多段特征[7～8](見圖7)。如果利用早期的壓降曲線，會導致動態(tài)儲量計算結果偏低。
 

    對于具有邊底水氣藏，地層水沿裂縫或高滲層侵入氣藏(氣井)后，占據(jù)了流動通道，封存了部分天然氣儲量，使參與流動的天然氣儲量降低，從而導致動態(tài)儲量減少[9]。
3 結論
    1) 動態(tài)儲量計算方法在傳統(tǒng)的物質平衡方法的基礎上，近期又出現(xiàn)了以單井生產歷史為基礎的“流動物質平衡法”和“氣井生產動態(tài)曲線特征圖版擬合法”動態(tài)儲量計算方法，針對日趨復雜的氣藏類型，新方法的出現(xiàn)提供了有利的分析手段。
    2) 要獲得可靠的動態(tài)儲量計算結果，關井測壓資料是必不可少的。
    3) 針對異常高壓氣藏進行動態(tài)法儲量計算時，必須獲得準確的巖石有效壓縮系數(shù)值。
    4) 氣藏的流動供給發(fā)生變化會導致動態(tài)儲量發(fā)生變化。
符號說明
    pi為原始地層壓力，MPa；p為地層壓力，MPa；△p為地層壓降，MPa；pwf為井底流運壓力，MPa；G為原始地質儲量，104m3；Gp為累產氣量，104m3；Bgi、Bi分別為對應壓力pi和p時氣體體積系數(shù)；Bw為地層水體積系數(shù)；We為累積水侵量，104m3；Wp為累產水量，104m3；Zi、Z分別為對應壓力pi和p時氣體壓縮因子；Cf為巖石有效壓縮系數(shù)，MPa-1；Cg、Cw分別為天然氣、地層水壓縮系數(shù)，MPa-1；Swi為原始含水飽和度，小數(shù)。
參考文獻
[1] 李士倫.天然氣工程[M].北京：石油工業(yè)出版社，2000.
[2] 陳元千，李璗.現(xiàn)代油藏工程[M].北京：石油工業(yè)出版社，2008.
[3] POSTON S W，B ERG R R. Overpressured gas reservoirs[M].Texas，U.S.A.：Society of Petroleum Engineers，Inc. Richardson，1997.
[4] MATTAR L，ANDERSON D.M.A Systematic and comprehensive methodology for advanced analysis of production data[C∥SPE Annual Technical Conference and Exhibition. Denver，Colorado：SPE，2003.
[5] ANDERSON D，MATTAR L. Practical diagnostics using production data and flowing pressures[C]∥SPE Annual Technical Conference and Exhibition. Houston，Texas：SPE，2004.
[6] ANDERSON D M，MATTARET L. Practical diagnostics using production data and flowing pressures[C]∥SPE Annual Technical Conference and Exhibition，San Antonio，Texas，USA：SPE，2006.
[7] 郝玉鴻，許敏，徐小蓉.正確計算低滲透氣藏的動態(tài)儲量[J].石油勘探與開發(fā)，2002，29(5)：66-68.
[8] 鄧紹強，黃全華，肖莉，等.低滲透氣藏儲量早期預測[J].西南石油學院學報，2005，27(6)：33-35.
[9] 王玉文.中壩氣田須二氣藏排水采氣開發(fā)效果分析及開發(fā)前景展望[J].天然氣工業(yè)，1995，15(5)：28-31.
 
(本文作者：劉曉華1，2 1.中國石油勘探開發(fā)研究院廊坊分院；2.中國地質大學(北京)能源學院)