甘笑非 周蓮蓮 趙玄之 淦文杰 周旻昊 高晨軒

中國石油西南油氣田公司川中油氣礦

0 引言

深層海相碳酸鹽巖氣藏深度普遍超過5 000 m，資源潛力大，但由于受到多期巖溶作用的影響，儲集空間以中小溶洞為主，次為粒間（溶）孔，孔洞間連通性差。儲層流體滲流規(guī)律也十分復(fù)雜，氣體滲流需要較長時間才能達(dá)到穩(wěn)態(tài)流動。

安岳特大型氣田高石梯—磨溪地區(qū)（以下簡稱為“高磨地區(qū)”）位于四川盆地中部川中古隆中斜平緩帶威遠(yuǎn)至龍女寺構(gòu)造群，上震旦統(tǒng)燈影組第四段（以下簡稱“燈四段”）為儲層以白云巖為主、受地層剝蝕尖滅、巖性及構(gòu)造等控制的高溫常壓大型巖性—地層復(fù)合圈閉氣藏。該氣藏儲層平均深度為5 200 m，全直徑巖心平均孔隙度為3.97%，平均滲透率為0.6 mD，優(yōu)質(zhì)儲層厚度分布在20～60 m，屬于超深、低孔、低滲海相碳酸鹽巖邊水氣藏。由于受成本、安全性等諸多因素制約，很難對氣藏所有氣井開展正規(guī)產(chǎn)能測試工作，現(xiàn)場基本采用常規(guī)一點法經(jīng)驗公式進(jìn)行產(chǎn)能計算，其結(jié)果與實際產(chǎn)能測試結(jié)果存在較大誤差，因此有必要建立適用于深層碳酸鹽巖氣藏的產(chǎn)能計算新方法。

本文利用修正一點法和回歸一點法分別建立新的產(chǎn)能評價模型，優(yōu)選既能滿足氣井產(chǎn)能的快速評價需求、又能減小產(chǎn)能評價結(jié)果誤差的評價方法。在高磨地區(qū)燈四段氣藏的實際應(yīng)用表明，修正一點法、回歸一點法均能降低常規(guī)一點法的計算誤差，適用于該類氣藏的快速產(chǎn)能評價，其中回歸一點法比修正一點法的計算誤差更小。

同時建立了第三種方法——求取二項式產(chǎn)能方程系數(shù)A、B的新方法，適用于產(chǎn)能測試時有3次以上變產(chǎn)量測試數(shù)據(jù)的氣井。這種方法比前兩種方法的計算誤差更小，產(chǎn)能計算結(jié)果更為準(zhǔn)確。

1 氣井產(chǎn)能方程的確定

計算氣井產(chǎn)能時，常規(guī)一點法和氣井二項式產(chǎn)能方程是最為常用的方法[1]。陳元千根據(jù)14個氣田和2個案例共計16口氣井的二項式產(chǎn)能方程修正系數(shù)α值統(tǒng)計得出：α＝0.25，在此基礎(chǔ)上建立了常規(guī)一點法產(chǎn)能評價方法[2]，該方法在致密砂巖氣藏得到了較好的應(yīng)用推廣。

對產(chǎn)能計算方法進(jìn)行改進(jìn)，有如下幾種方式[3-22]：①修正一點法，這種方法應(yīng)用較多，是在陳元千常規(guī)一點法的理論基礎(chǔ)上，結(jié)合正規(guī)產(chǎn)能測試數(shù)據(jù)確定修正系數(shù)（α）值，建立適用于本氣藏的快速產(chǎn)能評價方法；②回歸一點法，這種方法應(yīng)用相對較少，是建立變量參數(shù)半對數(shù)坐標(biāo)系，通過非線性回歸坐標(biāo)系中的散點分布獲得新的產(chǎn)能計算公式，回歸后的產(chǎn)能計算模型能夠與所有參與回歸的測試數(shù)據(jù)保持最小誤差；③確定二項式產(chǎn)能方程系數(shù)A、B的新算法，詳見后述。

1.1 修正一點法產(chǎn)能方程

修正一點法是在二項式產(chǎn)能方程基礎(chǔ)上建立的，利用每口產(chǎn)能測試井的二項式產(chǎn)能系數(shù)A、B求解每口井的α，再通過統(tǒng)計平均的方法獲得整個區(qū)塊的α值。該方法需要5口以上氣井的正規(guī)產(chǎn)能測試數(shù)據(jù)來降低求解α的誤差，對樣本的數(shù)量和數(shù)據(jù)的質(zhì)量（即測試數(shù)據(jù)的準(zhǔn)確性）有一定要求。

1.2 回歸一點法產(chǎn)能方程

回歸一點法是對常規(guī)一點法的一個改進(jìn)，通過變換陳元千一點法產(chǎn)能公式[式（7）]，將產(chǎn)量項和壓力項分別移動到等式的兩邊，得到新的方程形式為：

從變換后的公式可知，常規(guī)一點法產(chǎn)能公式是qg/qAOF與呈正相關(guān)，那么針對有產(chǎn)能試井資料的氣井，可以建立無因次產(chǎn)量（qg/qAOF）與無因次壓力的關(guān)系，通過回歸相關(guān)性變量的方法獲得回歸一點法的產(chǎn)能公式[11-22]。

回歸一點法基于常規(guī)一點法公式，但是不用求解α，而是將具有產(chǎn)能測試數(shù)據(jù)的井的測試點刻畫到散點坐標(biāo)系上，利用二項式曲線回歸散點的分布趨勢而獲得無因次產(chǎn)量和無因次壓力的關(guān)系式。該方法利用了每口井的變產(chǎn)量測試數(shù)據(jù)，因而在產(chǎn)能測試井?dāng)?shù)較少的情況下，依舊能夠依靠數(shù)據(jù)點的分布趨勢進(jìn)行回歸，具有較強的適用性。

1.3 確定系數(shù)A、B的常規(guī)算法

根據(jù)產(chǎn)能試井（多點回壓法）資料建立二項式產(chǎn)能方程的方法有多種，較為常用的是作圖法和最小二乘法。

1）作圖法

在普通方格紙上作圖，縱軸為（pR2-pwf2）/qg，橫軸為qg，用幾組實測試井?dāng)?shù)據(jù)回歸后應(yīng)為一直線，據(jù)此直線方程可以得到截距（A）、斜率（B）。

2）最小二乘法

依據(jù)試井實測數(shù)據(jù)，也可用最小二乘法來確定二項式產(chǎn)能方程的系數(shù)A、B，即

1.4 確定系數(shù)A、B的新算法

2 實例分析

高磨地區(qū)燈四段氣藏自2012年9月投產(chǎn)以來，先后開展了壓力恢復(fù)試井19井次、產(chǎn)能試井5井次。氣藏5口井的地層壓力（pR）、無阻流量（qAOF）、產(chǎn)能方程，即5口井的產(chǎn)能試井結(jié)果如表1所示。在此基礎(chǔ)上，對常規(guī)一點法、前述三種產(chǎn)能評價改進(jìn)方法進(jìn)行實際應(yīng)用，開展對比研究。

表1 高磨地區(qū)燈四段氣藏產(chǎn)能試井結(jié)果表

2.1 修正一點法實例分析

利用氣藏5口氣井的穩(wěn)定產(chǎn)能試井資料，在各氣井的二項式產(chǎn)能方程、無阻流量（qAOF）已知的情況下，對各測試井的修正系數(shù)（α）進(jìn)行了計算，得到α平均值為0.16（表2）。

表2 高磨地區(qū)燈四段氣藏產(chǎn)能試井結(jié)果及修正系數(shù)表

將5口井的修正系數(shù)（α）與A、B的對應(yīng)關(guān)系分別繪制在直角坐標(biāo)系中（圖1）可以看出，α與A呈線性正相關(guān)，且相關(guān)性較好，隨著A值增加，α值也增加；同時，α與B也呈正相關(guān)，但相關(guān)性較差。

圖1 5口測試井的二項式產(chǎn)能方程系數(shù)A、B與α相關(guān)性圖

式（4）顯示，修正系數(shù)（α）是由二項式產(chǎn)能方程的層流系數(shù)（A）、紊流系數(shù)（B）和絕對無阻流量（qAOF）這3個參數(shù)共同計算得出的。那么可在直角坐標(biāo)系中繪制α與3個參數(shù)變量組合的對應(yīng)關(guān)系（圖2），得到回歸公式，通過這種方式也可快速算得α。

圖2 二項式產(chǎn)能方程組合系數(shù)與α相關(guān)性圖

將氣藏的α平均值0.16代入式（6），得到適用于該氣藏的修正一點法無阻流量公式為：

利用式（21），可對氣藏其他沒有開展穩(wěn)定試井的井，用修正一點法代替常規(guī)一點法，進(jìn)行氣井產(chǎn)能的快速評價。同時，為了比對修正一點法比常規(guī)一點法在產(chǎn)能計算誤差上的提升效果，用這兩種方法對已有穩(wěn)定試井產(chǎn)能結(jié)果的5口井進(jìn)行計算，結(jié)果如表3所示。從表3可見，5口井采用常規(guī)一點法所得無阻流量與穩(wěn)定試井所得無阻流量的誤差平均在27.2%，而修正一點法將平均誤差降低到了23.2%，基本每口井都能夠降低誤差3%～5%。

表3 常規(guī)一點法、修正一點法評價氣井產(chǎn)能對比分析表

2.2 回歸一點法實例分析

同樣基于高磨地區(qū)燈四段氣藏5口氣井的穩(wěn)定試井測試資料，在半對數(shù)坐標(biāo)系中建立無因次產(chǎn)量（qg/qAOF）與無因次生產(chǎn)壓差（1-pwf2/pR2）之間的關(guān)系圖版（圖3），通過對散點的非線性回歸，而得到氣井產(chǎn)能的回歸一點法公式：

圖3 高磨地區(qū)燈四段氣藏產(chǎn)能評價一點法回歸分析圖

利用式（22），可對氣藏其他沒有開展穩(wěn)定試井的井，用回歸一點法代替常規(guī)一點法，進(jìn)行氣井產(chǎn)能的快速評價。兩種方法的計算結(jié)果與誤差對比見表4?；貧w一點法產(chǎn)能計算的平均誤差為13.9%，較之常規(guī)一點法的平均誤差27.2%、修正一點法的平均誤差23.2%，可見回歸一點法進(jìn)一步降低了計算誤差，顯著提升了一點法產(chǎn)能評價的適用性和準(zhǔn)確性。

修正一點法、回歸一點法這兩種方法均能快速評價氣井產(chǎn)能，同時降低常規(guī)一點法的計算誤差，具有較好的普適性。但兩種方法的誤差水平同樣取決于產(chǎn)能測試井的樣本數(shù)量，樣本越多，誤差越?。粯颖驹缴?，誤差越大。

2.3 確定A、B的新方法實例分析

應(yīng)用式（20），采用規(guī)劃求解對氣藏已有穩(wěn)定試井產(chǎn)能評價成果的5口井進(jìn)行了產(chǎn)能計算，并對結(jié)果進(jìn)行誤差分析（表5）。從表中可見，采用這種新方法計算出的無阻流量（qAOF）平均誤差僅為8.60%；計算出的地層壓力（pR）較實測值更為接近，平均誤差僅為2.00%?？梢娺@種新方法是準(zhǔn)確可靠的。該方法適用于求解單井的二項式產(chǎn)能方程系數(shù)A、B及地層壓力，由于引入了地層壓力作為未知量，所以需要3次以上變產(chǎn)量的測試數(shù)據(jù)才能夠求解。

表5 A、B新算法所得氣井產(chǎn)能、地層壓力計算結(jié)果與穩(wěn)定試井成果對比分析表

3 結(jié)論

1）在不能全面開展正規(guī)產(chǎn)能測試（即穩(wěn)定試井）的情況下，修正一點法、回歸一點法均能減小常規(guī)一點法的計算誤差，能夠快速、更準(zhǔn)確地評價氣井產(chǎn)能。

2）回歸一點法比修正一點法的計算結(jié)果誤差更小。

3）建立了確定二項式產(chǎn)能方程系數(shù)A、B的新方法，引入了地層壓力作為未知量，在準(zhǔn)確獲取氣井產(chǎn)能的同時也確定了地層壓力，適用于產(chǎn)能測試時有3次以上變產(chǎn)量測試數(shù)據(jù)的氣井。

4）通過確定二項式產(chǎn)能方程系數(shù)A、B的新方法，計算得到的氣井產(chǎn)能結(jié)果，比修正一點法、回歸一點法的計算結(jié)果誤差更小，準(zhǔn)確性更高。