劉忠華 祁俊武 崔艷敏 李小剛

(1.重慶科技學院復雜油氣田勘探開發(fā)重慶市重點實驗室，重慶 401331;2.長慶井下技術作業(yè)公司，西安 710018;3.重慶科技學院組織人事部，重慶 401331)

目前，國內(nèi)外對穩(wěn)定試井的解釋方法常用的有2種:二項式法和指數(shù)式法。前者是理論方法，后者是經(jīng)驗方法。氣井產(chǎn)能指數(shù)式法是分析氣井動態(tài)和確定氣井無阻流量的重要方法之一，其提出已有近百年的時間。由于它使用簡便，在國外尤其是在美國得到了廣泛的應用［1］。指數(shù)式屬于達西穩(wěn)定流產(chǎn)能方程的修正式，因此又稱為經(jīng)驗方法。國內(nèi)的學者主要是用該方法預測氣井的無阻流量等［1－5］。

本文依托第二屆全國石油工程設計大賽的原始數(shù)據(jù)，從氣井指數(shù)式產(chǎn)能方程出發(fā)，推導出了符合指數(shù)式產(chǎn)能方程的氣井的無因次IPR方程，并繪制了氣井無因次IPR曲線圖版。進一步應用無因次IPR方程推導出了氣井的預測IPR方程，并繪制了4口氣井的IPR曲線。通過實例驗證該無因次IPR曲線圖版能夠快速準確的擬合氣井動態(tài)指數(shù)，氣井預測IPR曲線能夠很好的預測氣井的井底流入動態(tài)。

1 方法推導

Rawlins等人提出修正的氣井指數(shù)式產(chǎn)能方程為:

式中:qg—氣井產(chǎn)量，104m3d;pR—地層壓力，MPa;pwf—井底流動壓力，MPa;c—氣井產(chǎn)能系數(shù);n —氣井的動態(tài)指數(shù)(0.5≤n≤1.0)，當 n=0.5時，氣井的生產(chǎn)動態(tài)完全為湍流控制，當n=1.0時，氣井的生產(chǎn)動態(tài)完全為達西層流所控制，實際中 n取 0.6 ～0.9［1］。

為了確定c和n的值，將式(1)兩邊取常用對數(shù)，有

當 pwf=0.101 MPa時=qAOF，得到利用指數(shù)式產(chǎn)能方程確定氣井絕對無阻流量的公式為:

將式(3)兩邊取常用對數(shù)，有

將式(8)變形得到氣井指數(shù)式無因次IPR方程為:

將式(6)改寫成下式:

同時，將式(7)改寫成下式:

將式(10)和式(11)代入式(9)，得到氣井指數(shù)式IPR方程為:

2 無因次IPR曲線圖版的建立與應用

2.1 無因次IPR曲線圖版的建立

根據(jù)式(9)，氣井動態(tài)指數(shù)n分別取0.50、0.55、0.60、0.65、0.70、0.75、0.80、0.85、0.90、0.95、1.00，使用Matlab軟件繪制氣井指數(shù)式無因次IPR曲線圖版如圖1所示。從圖中可以看出，指數(shù)式無因次IPR曲線在實數(shù)坐標系中為一曲線簇。氣井動態(tài)指數(shù)n越大，無因次IPR曲線遞減越快;反之，越慢。

圖1 氣井指數(shù)式無因次IPR曲線圖

2.2 氣井動態(tài)指數(shù)n的擬合

收集整理了 A3.1 井［2］、克 006 井［3］、50056井［3］、克77 井［4］的一點法測試數(shù)據(jù)以及文中計算無阻流量數(shù)據(jù)(見表1)。然后計算各井的無因次壓力、無因次產(chǎn)量和動態(tài)指數(shù)n(見表1)。使用計算的無因次壓力和無因次產(chǎn)量結果，在氣井指數(shù)式無因次IPR曲線圖版上描點(如圖2所示)，分別得到了4 口井的氣井動態(tài)指數(shù) n 為 0.63、0.80、0.80、0.94，擬合值與計算值完全一致。說明該圖版擬合效果好，精度高。另外，從數(shù)據(jù)處理過程和擬合過程來看，使用氣井指數(shù)式無因次IPR曲線圖版擬合求取氣井動態(tài)指數(shù)n方便快捷。

圖2 氣井動態(tài)指數(shù)n的擬合

2.3 預測IPR曲線

整理M1井、M4井、M5井和M6井的系統(tǒng)試井測試數(shù)據(jù)如表2所示(資料來源于第二屆全國石油工程設計大賽)，根據(jù)式(2)分別繪制了4口井的指數(shù)式關系圖如圖3至圖6所示，可以看出擬合精度較高，擬合效果較好。4口井的氣井動態(tài)指數(shù)n和氣井產(chǎn)能系數(shù)c見表2。

表1 氣井一點法測試數(shù)據(jù)表

表2 氣井系統(tǒng)試井數(shù)據(jù)及指數(shù)式產(chǎn)能方程計算結果

根據(jù)式(12)，利用計算得到的參數(shù)值繪制了4口井的指數(shù)式關系(見圖3—圖6)和理論IPR曲線(見圖7)，再將實測得到的氣井產(chǎn)量和井底流壓值在相同坐標系中描點，結果實測值完全在理論IPR曲線上，說明理論計算結果準確。使用本文推導關系式預測氣井井底流入動態(tài)準確可行。

圖3 M1井指數(shù)式關系圖

圖4 M4井指數(shù)式關系圖

圖5 M5井指數(shù)式關系圖

圖6 M6井指數(shù)式關系圖

圖7 氣井預測IPR曲線圖

3 結論

基于氣井的指數(shù)式產(chǎn)能方程，繪制了氣井無因次IPR曲線圖，提出了圖版擬合確定氣井動態(tài)指數(shù)的方法，并利用無因次IPR方程推導出預測氣井井底流入動態(tài)的方程。經(jīng)實例驗證該理論正確，計算和擬合結果準確可靠。

［1］陳元千.油氣藏工程實踐［M］.北京:石油工業(yè)出版社，2011:161-166.

［2］李治平，鄔云龍，青永固，等.氣藏動態(tài)分析與預測方法［M］.北京:石油工業(yè)出版社，2002:49.

［3］劉英，胡新平，柳海，等.一種新的一點法穩(wěn)定試井在克拉瑪依氣藏中的應用［J］.新疆石油地質(zhì)，1999，20(4):343-347.

［4］張興權.大牛地氣田一點法試井對計算無阻流量的影響［J］.油氣井測試，2007，16(6):19-20.

［5］楊洪志，馮文光，陳友蓮.二項式產(chǎn)能試井擬合圖版［J］.礦物巖石，2000，20(2):79-81.