摘要:使用二項式產能解釋模型、指數式氣井產能解釋模型、一般單點法氣井產能解釋模型、指數式單點法氣井產能解釋模型和二次單點法氣井產能解釋模型,分別對澀北氣田氣井產能試井資料進行解釋與對比,篩選了適合于澀北氣田的氣井產能試井解釋模型。利用防砂氣井措施前后的產能試井資料擬合得到IPR曲線,評價了澀北氣田80井次的防砂措施改善井底流動條件的效果,其評價結果與氣井生產動態分析結果比較一致,具有較好的可信度,為澀北疏松砂巖氣藏防砂效果評價找到了一種新的手段。
關鍵詞:澀北氣田;氣井;試井;生產能力;防砂;評價
1 氣井產能試井解釋模型
產能試井測得不同井底流壓下的氣井產量,通過解釋分析可擬合IPR曲線并計算出無阻流量等參數,目前產能試井解釋模型主要有以下5種[1~4]。
1.1 二項式產能模型
氣井產能方程的二項式形式為:
pr2-pwf2=AQsc+BQsc2 (1)
式中:pr為平均氣藏壓力,MPa;pwf為井底流壓,MPa;Qsc為地面標況下產氣量,104m3/d;A、B分別為回歸系數。
在使用產能試井資料A、B系數進而計算IPR曲線時,如果未知地層靜壓,需要首先根據測試資料先求出地層靜壓,然后再擬合系數A、B。
絕對無阻流量(QAOF)為:

1.2 指數式氣井產能解釋模型
Rawlins和Schelhardt根據大量氣井生產數據總結出氣井產能方程的指數彩式為:
Qsc=C(pr2-pwf2)n (3)
當n=1時,氣流入井相當于層流,井底沒發生與流量相關的表皮系數,完全符合達西滲流規律;當0.5<n<1時,表示氣流入井符合非達西流動規律。
實際礦場應用中,根據指數方程直接計算C系數的做法很少。一般是根據系統試氣或生產測試資料反求系數C和n。對指數式產能方程兩端取對數,使用一組Q-p資料的線性回歸關系曲線,利用回歸直線的斜率和截距可計算參數C、n,進而繪制IPR曲線和計算絕對無阻流量。
QAOF=C(pr2-0.1012)n (4)
1.3 一般單點法氣井產能解釋模型
一般單點法氣井IPR模型為:

如果已知地層靜壓為pr,則方程中只有一個待定系數QAOF。分別將測試數據點代入可線性擬合得到絕對無阻流量(QAOF)。
1.4 指數式單點法氣井產能解釋模型
指數式單點法氣井IPR模型為:

1.5 二次單點法氣井產能解釋模型
二次單點法氣井產能IPR方程為:

2 氣井產能試井解釋及模型篩選
目前澀北氣田對于防砂井進行了大量的產能試井,得到了大量試井資料。使用上述5種模型分別對試井資料進行解釋,分析其相關性,進而篩選合適的結束模型。以澀北氣田澀2-4井為例進行試井解釋和分析,該井在2005年5月19日防砂前和2005年7月19日防砂后分別進行了產能試井,試井數據如表1所示。
表1 澀2-4井兩次產能試井數據表
2005年5月19日產能試井
|
2005年7月19日產能試井
|
||||||
氣嘴(mm)
|
氣量(104m3/d)
|
流壓(MPa)
|
靜壓(MPa)
|
氣嘴(mm)
|
氣量(104m3/d)
|
流壓(MPa)
|
靜壓(MPa)
|
3.0
|
1.6176
|
11.3559
|
11.5410
|
3.0
|
2.1480
|
11.3150
|
11.628
|
4.0
|
2.6107
|
11.2049
|
11.5148
|
4.0
|
2.5551
|
11.1670
|
11.704
|
4.5
|
3.2814
|
22.1009
|
11.4889
|
4.5
|
3.1536
|
10.9600
|
11.690
|
5.0
|
3.8434
|
10.9899
|
11.4659
|
5.0
|
3.7623
|
10.7650
|
11.677
|
5.5
|
4.7978
|
10.8210
|
11.4459
|
5.5
|
4.5089
|
10.5177
|
11.655
|
使用這5種方法對該井的試井資料進行解釋分析并繪制氣井的流入動態曲線,二項式模型在計算過程中由于有開方計算,在這兩組數據進行計算過程中出現負值,不能正常進行計算,后4種模型計算得到的IPR曲線如圖1所示。
由圖1可以看出不同的模型計算擬合出的IPR曲線不同,無阻流量的計算結果也不同,IPR曲線與橫坐標的交點處即為無阻流量,不同計算模型中所計算出的回歸系數,無阻流量及模型與實際數據點的相關系數如表2所示。
表2 澀2-4井產能試井解釋計算結果數據表
IPR模型
|
2005年5月19日產能試井
|
2005年7月19日產能試井
|
||||
回歸系數
|
無阻流量(104m3/d)
|
相關系數
|
回歸系數
|
無阻流量(104m3/d)
|
相關系數
|
|
指數式模型
|
0.35,0.92
|
32.31120
|
0.892
|
0.02,1.51
|
29.37121
|
0.889
|
一般單點法
|
18.0069
|
18.00617
|
0.958
|
12.4938
|
12.49327
|
0.955
|
指數式單點法
|
19.2393
|
20.07329
|
0.978
|
13.1756
|
13.74677
|
0.958
|
二次單點法
|
21.5462
|
21.54589
|
0.969
|
12.6218
|
12.62166
|
0.952
|
由計算結果可以看出成功計算的4種模型中指數模型的相關系數較低,一般單點法、指數式單點法和二次單點法的相關系數都比較高,說明對于澀北氣田這3種模型的擬合效果較好,在進行產能試井解釋時應優先選用這3種計算方法。使用其余15井次的產能試井資料解釋并對比相關性,驗證了上述結論。
3 產能試井分析在防砂井效果評價中的應用
澀北氣田目前實施了105井次的防砂措施,對這些防砂措施進行客觀準確的綜合效果評價有利于總結成功的經驗,探索失敗的原因,進一步提高防砂效果。防砂措施對氣井生產動態主要產生了擋砂作用、增產作用和改善井底流動條件作用。對應的防砂效果評價也分為擋砂效果評價、增產效果評價與改善井底流動條件效果評價3個層次。其中改善井底流動條件效果評價可以真正體現防砂措施對氣井生產動態的實質影響[5]。
根據防砂前后的產能試井資料計算出的無阻流量可計算相應的變化比,變化比的大小反應防砂措施改善井底流動條件的好壞。假設防砂前、后的氣井無阻流量為QAOF0和QAOF1,則變化比(PI)為:

氣井無阻流量的變化比越大,防砂效果越好;反之,防砂效果越差。而氣井的無阻流量需要根據產能試井資料解釋氣井流入動態關系曲線而得到[6]。
澀4-9井在2003年9月28日實施高壓一次充填作業,并且分別在防砂前20Q3年8月12日和防砂后2003年11月24日進行了產能試井,對試井資料使用一般單點法、指數式單點法和二次單點法進行分析處理,得出措施前后的IPR對比曲線如圖2所示。
經計算該井措施前的無阻流量為26.7×104m3/a,措施后的無阻流量為16.8×104m3/a,無阻流量的變化比PI=0.37??梢钥闯龇郎按胧┖蟮臒o阻流量降低,說明措施增加了井底表皮系數,防砂效果差。
澀3-7井在2002年10月12日實施高壓一次充填作業,并且分別在防砂前2002年8月1日和防砂后2002年11月1日進行了產能試井,對試井資料使用一般單點法、指數式單點法和二次單點法進行分析處理,得出措施前后的IPR對比曲線(圖3)。
經計算該井措施前的無阻流量為20.2×104m3/a,措施后的無阻流量為31.3×104m3/a,無阻流量的變化比為0.55??梢钥闯龇郎按胧┖蟮臒o阻流量增加,說明措施改善了井底流動條件,降低了表皮系數,防砂效果好。
使用上述防砂措施效果評價方法,以產能試井資料解釋為基礎,評價了80井次的防砂措施改善井底流動條件效果,評價結果與氣井生產動態分析結果比較一致,評價結果具有較好的可信度。
4 結論
通過使用5種主要的氣井產能試井解釋模型,分別對澀北氣田氣井產能試井資料進行解釋與對比,篩選了適合于澀北氣田的氣井產能試井解釋模型。適合于澀北氣田氣井的產能試井解釋模型依次為一般單點法、指數式單點法和二次單點法。防砂效果評價也分為擋砂效果評價、增產效果評價與改善井底流動條件效果評價3個層次。其中改善井底流動條件效果評價可以真正體現防砂措施對氣井生產動態的實質影響。以產能試井資料解釋為基礎,評價了澀北氣田80井次的防砂措施改善井底流動條件效果,評價結果與氣井生產動態分析結果比較一致,評價結果具有較好的可信度。
參考文獻
[1] 張文彪,劉啟國,羅玉合.結合二項式與指數式產能方程求取氣井無阻流量[J].內蒙古石油化工,2008,18(10):91-92.
[2] 李士倫.天然氣工程[M].北京:石油工業出版社,2000.
[3] 王怒濤,黃炳光,梁文斌,等.氣井產能分析方法研究[J].大慶石油地質與開發,2004,23(1):33-34.
[4] 王德山,李兆敏,聶立新.水平分支井試井解釋方法研究[J].西南石油大學學報,2007,29(1):70-74.
[5] 張國東.新場氣田氣井產能和無阻流量的確定[J].油氣井測試,2003,22(1):4-6.
[6] 董長銀,李長印,扈福堂,等.油氣井防砂效果評價方法體系的研究與應用[J].油氣地質與采收率,2009,16(1):103-106.
(本文作者:尉亞民1,2 汪天游2 馮勝利2 董長銀1 李進秀2 1.中國石油大學(華東)石油工程學院;2.中國石油青海油田公司鉆采工藝研究院)
贊 賞 分享
您可以選擇一種方式贊助本站
支付寶轉賬贊助
微信轉賬贊助

- 注解:本資料由會員及群友提供僅供閱讀交流學習,不得用于商業用途!