鄭臘年，安紅燕，達麗娜，霍明會

（1.中國石油長慶油田分公司蘇里格氣田研究中心氣藏動態(tài)研究，陜西西安　710018；2.低滲透油氣田勘探開發(fā)國家工程實驗室，陜西西安　710018；3.中國石油長慶油田分公司，陜西西安　710018）

?

蘇里格氣田含水氣藏廢棄地層壓力確定的新方法

鄭臘年1，2，安紅燕1，2，達麗娜3，霍明會1，2

（1.中國石油長慶油田分公司蘇里格氣田研究中心氣藏動態(tài)研究，陜西西安710018；
2.低滲透油氣田勘探開發(fā)國家工程實驗室，陜西西安710018；3.中國石油長慶油田分公司，陜西西安710018）

摘要：蘇里格氣田是典型的致密砂巖氣藏，儲層不但致密，非均質(zhì)性強，氣井產(chǎn)量低，且多數(shù)產(chǎn)水，氣藏采收率低。廢棄地層壓力是氣藏采收率的關(guān)鍵影響因素之一，對于不產(chǎn)水氣藏可以通過單相二項式產(chǎn)能方程結(jié)合經(jīng)濟極限產(chǎn)量綜合確定，但對于產(chǎn)水氣藏廢棄地層壓力的確定國內(nèi)外研究不多。本文針對蘇里格氣田產(chǎn)水氣藏，以最小攜液速度下的產(chǎn)氣量與產(chǎn)水時經(jīng)濟極限產(chǎn)氣量相結(jié)合共同確定了廢棄產(chǎn)量區(qū)間范圍，應(yīng)用質(zhì)量守恒原理推導(dǎo)得出氣水兩相產(chǎn)能公式，最終求取含水氣藏廢棄地層壓力的一種新方法，并對影響因素進行分析。本文利用VC編程采用多次迭代計算，提高了計算精度，為產(chǎn)水氣藏廢棄地層壓力的求取及采收率的確定提供了一種新思路。

關(guān)鍵詞：含水氣藏；二項式產(chǎn)能方程；廢棄產(chǎn)量；廢棄地層壓力

廢棄地層壓力是計算氣藏采收率或可采儲量的重要參數(shù)，也是有關(guān)地面工程論證和設(shè)計的重要指標和依據(jù)。對于不產(chǎn)水氣藏可以通過單相二項式產(chǎn)能方程結(jié)合經(jīng)濟極限產(chǎn)量綜合確定［1-3］，但對于產(chǎn)水氣藏，目前有余貝貝、唐海等［4］提出的以氣井最小攜液速度為約束條件，通過復(fù)合半徑模型確定地層產(chǎn)水時氣井的穩(wěn)定二項式產(chǎn)能方程進行求取，但它存在的問題是復(fù)合半徑參數(shù)求取的不確定性，以及復(fù)合半徑模型不能得出不同水氣比下的傷害程度，給這種方法的應(yīng)用帶來一定的局限性。

本文根據(jù)質(zhì)量守恒原理推導(dǎo)出來的產(chǎn)水氣井二項式產(chǎn)能方程，以多種集輸條件下的氣井最小攜液流量和產(chǎn)水時經(jīng)濟極限產(chǎn)量為約束，最終獲得了產(chǎn)水氣井廢棄地層壓力，并對影響廢棄地層壓力因素進行分析。

1　含水氣藏廢棄地層壓力確定方法

1.1廢棄井口壓力

自噴開采井廢棄井口壓力等于最低外輸氣壓力，蘇里格上古井流程外輸壓力為3MPa，下古井流程外輸壓力為6.4MPa。隨著氣田深入開發(fā)，在一級增壓以及多級增壓的情況下，氣井的最小集輸井口壓力可低至0.5MPa。

1.2廢棄井口產(chǎn)量

對于產(chǎn)水氣井，隨著地層壓力的下降，氣井最大產(chǎn)氣量不斷降低，當最大產(chǎn)氣量低于最小攜液速度下的產(chǎn)氣量時，此時氣井將不能正常生產(chǎn)，對應(yīng)產(chǎn)氣量即為廢棄井口產(chǎn)量。

目前，比較常見的臨界攜液流量模型主要有Turner模型、李閩模型和王毅忠模型。對于蘇里格氣田，王毅忠模型與實際生產(chǎn)情況比較接近，故本文采用王毅忠模型來計算廢棄井口壓力下對應(yīng)的臨界攜液流量即廢棄井口產(chǎn)量：

其中：

式中：qSC-最小攜液產(chǎn)氣量，m3/d；pwf-井底流動壓力，MPa；vg-氣體流速；A=πd2/4-油管內(nèi)截面積，m2；d-油管內(nèi)徑，m；T-氣流溫度，K；Z-天然氣偏差系數(shù)；σ-界面張力，對水取60 mN/m；ρl-液體密度，kg/m3，水取1 074 kg/m3；ρg-天然氣密度，kg/m3；γg-天然氣相對密度。

1.3廢棄井底流壓

已知廢棄井口壓力和井口產(chǎn)量，代入到井底流壓計算公式中，便可算出廢棄井底流壓，其公式如下：

式中：pwh－井口油壓，MPa；H－氣層中深，m；e1－管壁粗糙度，m；μg－一定溫度壓力下的天然氣黏度，mPa·s。

由于垂直段隨著井深的變化井筒的壓力和溫度都會發(fā)生變化，從而引起了偏差系數(shù)和黏度的變化，為了計算更加準確，本文將垂直井段以每米為單位劃分單元，每個單元計算出相應(yīng)的平均偏差系數(shù)和黏度進行迭代計算。步驟如下：

（1）劃分單元，確定基本參數(shù)。

（2）假設(shè)下一單元的初始壓力，讓它等于上一單元的壓力。取它們的平均值求出本段相應(yīng)的壓力。

（3）根據(jù)平均壓力和平均溫度計算出擬壓力和擬溫度，進而計算出下一段的偏差系數(shù)和黏度等參數(shù)。

（4）采用式（4）求出下一段壓力。

（5）用計算出來的壓力跟假設(shè)的壓力值作比較，若兩者比較接近時，則計算出來的壓力即為所求，否則將計算出的壓力作為迭代值，重復(fù)步驟（2）～（4）繼續(xù)計算直到滿足精度為止。

（6）重復(fù)上面的步驟，直到計算到井底。

根據(jù)上述步驟，利用VC編制相應(yīng)程序，從而得出氣井垂直井筒壓力分布曲線（見圖1），從圖1中可以看出井筒壓力隨井深呈線性增加。

圖1　井筒壓力隨井深變化的關(guān)系曲線

1.4廢棄地層壓力確定方法

1.4.1產(chǎn)水氣井產(chǎn)能方程的確定對于產(chǎn)水氣井，引入氣水兩相擬壓力函數(shù)。

同時考慮表皮系數(shù)和高速非達西滲流的影響，推導(dǎo)出產(chǎn)水氣井二項式產(chǎn)能方程，如下式所示：

其中：

式中：ψ（p）-擬壓力，kg·MPa/（m3·mPa·s）；krg-氣相相對滲透率；krw-水相相對滲透率；μw-一定溫度壓力下的水相黏度，mPa·s；ρgsc-地面標準狀況下氣體密度，kg/m3；ρwsc-地面標準狀況下水密度，kg/m3；WGR-水氣體積流量之比，m3/m3；re-供給半徑，m；rw-井口半徑，m；S-表皮系數(shù)；k-儲層滲透率，10-3μm2；h-氣層有效厚度，m。

該產(chǎn)能方程中引入兩相擬壓力公式比用壓力平方計算精度高，此外引入水氣比這個參數(shù)，進而可以計算不同水氣比對廢棄地層壓力影響。

1.4.2兩相擬壓力與地層壓力的關(guān)系根據(jù)儲層巖心氣水相對滲透率曲線，可以得出兩相擬壓力與地層壓力的關(guān)系曲線（見圖2），從圖2中可以看出，隨著壓力的增大，氣水兩相擬壓力隨之增大。

圖2　兩相擬壓力與壓力的關(guān)系曲線

1.4.3二項式產(chǎn)能系數(shù)與地層壓力的關(guān)系從式（10）、（11）中可以看出，在氣井開采過程中假設(shè)表皮系數(shù)S不變，當?shù)貙訅毫ο陆禃r，二項式系數(shù)A值不變，式（11）中的黏度將會發(fā)生變化從而影響B(tài)值的變化。假設(shè)氣井在不同開采階段的二項式產(chǎn)能系數(shù)分別為B1和B2且對應(yīng)的天然氣黏度分別為μg1和μg2，那么根據(jù)式（11）有：

根據(jù)上式作出二項式系數(shù)B隨地層壓力變化的關(guān)系曲線（見圖3）。

圖3　B值隨地層壓力變化曲線

由圖3知，B值隨地層壓力的降低逐漸增加，因此若始終采用初期建立的穩(wěn)定產(chǎn)能方程來計算廢棄地層壓力必將帶來誤差。

1.4.4廢棄地層壓力計算由1.2和1.3節(jié)確定出廢棄產(chǎn)量和廢棄井底流壓后，代入到二項式產(chǎn)能公式（9）中，從而可以求出對應(yīng)的廢棄地層壓力。但是具體計算中存在的難點是：在求解二項式系數(shù)B時需要先得到廢棄條件下的天然氣黏度，而求解又需要知道廢棄地層壓力，然而氣井的廢棄地層壓力本身就是需要求解的參數(shù)，是未知的，因此需要采用迭代計算方法。利用VC編制程序，求解廢棄地層壓力，其具體步驟為：

（1）令廢棄地層壓力p的初始迭代值p0為原始地層壓力；

（3）按式（9）計算ψ（p），根據(jù)ψ（p）與p的關(guān)系，求出對應(yīng)的p值；

（4）若計算出來的值p與假定的值p0相差很小，則計算出來的壓力即為所求的廢棄地層壓力；否則，取p0=p，重復(fù)（2）～（3）步，直到滿足要求。

2　實例應(yīng)用

本文以SX區(qū)塊為例，該區(qū)塊滲透率介于0.06 mD～2.457 mD，平均值為0.510 mD。有效厚度介于2m～22.7m，平均值為8.6m，氣層平均中部埋藏深度3 556m，原始地層壓力31MPa，水氣比平均為0.8m3/104m3。

考慮多種集輸條件下的情況，計算出對應(yīng)的臨界攜液流量，進而求出廢棄井底流壓，再利用產(chǎn)水氣井二項式產(chǎn)能公式，得到產(chǎn)水氣藏廢棄地層壓力，同時計算出產(chǎn)水經(jīng)濟極限產(chǎn)量下對應(yīng)的廢棄地層壓力，其結(jié)果（見表1）。

表1　不同集輸條件下廢棄地層壓力情況

從表1中可以看出，不同集輸條件下，廢棄地層壓力差異比較大。在自然生產(chǎn)條件下，其臨界攜液流量即為其極限產(chǎn)量，計算出來的廢棄地層壓力要大于采取措施下的廢棄地層壓力，因此對于產(chǎn)水氣藏廢棄地層壓力最小值為采取措施情況下計算的廢棄地層壓力值，最大值是自然生產(chǎn)遞減到臨界攜液流量時對應(yīng)的廢棄地層壓力，此壓力范圍對產(chǎn)水氣藏的開發(fā)具有實際指導(dǎo)意義。

此外，考慮產(chǎn)水井需采取排水采氣措施，其經(jīng)營成本有所增加，計算出產(chǎn)水和不產(chǎn)水情況下對應(yīng)的廢棄地層壓力，其結(jié)果（見表2）。

表2　產(chǎn)水和不產(chǎn)水情況下廢棄地層壓力對比表

從表2中可以看出，水的存在會使廢棄地層壓力增加，同時井口集輸壓力越小，水的影響越大，廢棄地層壓力增加的幅度也越大。因此對于產(chǎn)水氣井廢棄地層壓力的計算需要采用產(chǎn)水氣井二項式產(chǎn)能方程，不能用單相產(chǎn)能公式代替。

3　廢棄地層壓力影響因素分析

根據(jù)上述理論公式，結(jié)合SX區(qū)塊各項基本參數(shù)，分別作出氣層滲透率、有效厚度、井口最低集輸壓力、水氣比對廢棄地層壓力的影響曲線（見圖4～圖7）。

3.1滲透率的影響

從圖4可以看出，在滲透率較低時，廢棄地層壓力隨滲透率的增加明顯減小，而在滲透率較高時，廢棄地層壓力隨滲透率的增加變化很小，有一明顯分界點，即在滲透率高于0.2 mD后廢棄地層壓力趨于平緩。說明在低滲條件下，廢棄地層壓力對滲透率非常敏感，而在滲透率較高時則影響很小。

圖4　滲透率對廢棄地層壓力的影響

3.2氣層有效厚度的影響

考慮不同地層滲透率條件下有效厚度變化對廢棄地層壓力的影響，從圖5中可以看出，隨著有效厚度的增加，廢棄地層壓力逐漸降低。此外，當儲層滲透率較低時，隨有效厚度的增加，廢棄地層壓力降低幅度比較大；而在滲透率較高時，廢棄地層壓力降低的幅度較小。說明儲層滲透率較低時，廢棄地層壓力對有效厚度比儲層滲透率高時更敏感。

圖5　有效厚度對廢棄地層壓力的影響

3.3井口最低集輸壓力的影響

不同井口最低集輸壓力下廢棄地層壓力隨有效厚度變化的關(guān)系曲線，從圖6中可以看出，廢棄地層壓力隨井口最低集輸壓力的增加而增大。

圖6　井口最低集輸壓力對廢棄地層壓力的影響

圖7　水氣比對廢棄地層壓力的影響

3.4水氣比的影響

不同水氣比條件下廢棄地層壓力隨有效厚度變化的關(guān)系曲線，從圖7中可以看出，廢棄地層壓力隨水氣比的增加而增大。由此可知，水的存在對廢棄地層壓力的影響不容忽視。

4　結(jié)論

（1）針對產(chǎn)水區(qū)塊，在自然生產(chǎn)情況下，提出了以臨界攜液流量為極限產(chǎn)量，利用產(chǎn)水氣井二項式產(chǎn)能方程評價廢棄地層壓力。該方法計算的結(jié)果大于采取措施情況下的廢棄地層壓力，為廢棄地層壓力的上限值，對產(chǎn)水氣藏的開發(fā)有實際指導(dǎo)意義。

（2）隨著氣井地層壓力的下降，二項式系數(shù)A值不變，B值不斷增加，因此必需用廢棄時的產(chǎn)能方程計算廢棄地層壓力。

（3）廢棄地層壓力隨儲層滲透率的增加逐漸降低，當滲透率小于0.2 mD，滲透率對廢棄地層壓力的影響十分明顯；當滲透率大于0.2 mD，滲透率對廢棄地層壓力的影響不明顯。

（4）廢棄地層壓力隨有效厚度的增加而減小，且儲層滲透率較低時，廢棄地層壓力對有效厚度比儲層滲透率高時更敏感。

（5）隨著井口最低集輸壓力和水氣比的增加，廢棄地層壓力逐漸增大，水的存在對廢棄地層壓力的影響不容忽視。

參考文獻：

［1］郝玉鴻，閻紀輝.確定氣井廢棄地層壓力的新方法［J］．油氣采收率技術(shù)，1999，6（4）：77-82．

［2］鄒建波，李閩，代平，等．氣藏廢棄壓力和采收率影響因素［J］.新疆石油地質(zhì)，2006，27（6）：736-739．

［3］盧濤，蘭義飛，何磊．一種計算氣藏廢棄壓力的方法-二項式產(chǎn)能方程法［J］．天然氣勘探與開發(fā)，2009，32（2）：18-20．

［4］余貝貝，唐海，李東林，等．產(chǎn)水氣藏廢棄壓力的確定方法［J］．新疆石油地質(zhì)，2008，29（5）：641-643.

A new model of determination of abandonment pressure in water bearing gas reservoir in Sulige gasfield

ZHENG Lanian1，2，AN Hongyan1，2，DA Lina3，HUO Minghui1，2
（1.Sulige Gasfield Research Center of PetroChina Changqing Oilfield Company，Xi'an Shanxi 710018，China；2.National Engineering Laboratory for Exploration & Development of Low-Permeability Oil & Gasfield，Xi'an Shanxi 710018，China 3.PetroChina Changqing Oilfield Company，Xi'an Shanxi 710018，China）

Abstract:Sulige gasfield is typical of tight sandstone gas reservoirs. The reservoir not only compact，strong heterogeneity，low gas production，and high water production. The gas reservoir recovery has many factors. Abandonment formation pressure is one of the key factors that influence gas reservoir recovery. For no water gas reservoir can be integrated to determine by single-phase binomial productivity equation combined economic output limit，but for producing water gas reservoir，research of abandonment formation pressure is not much. For Sulige gasfield，this paper proposed a new method for calculating abandonment formation pressure，which used gas water two phase productivity formula derived from principle of mass conservation and combine with abandonment gas production rate range by determined by gas production at the minimum liquid carring rate and economic limit production of waterproduction. To improve accuracy，we use VC programming using multiple iterative calculations. This paper provides a new idea for caculating water gas reservoir abandonment formation pressure and recovery.

Key words：water bearing gas reservoir；binomial deliverability equation；abandonment production rate；abandonment formation pressure

中圖分類號：TE377

文獻標識碼：A

文章編號：1673-5285（2016）06-0075-05

DOI:10.3969/j.issn.1673-5285.2016.06.019

*收稿日期：2016－05-16

作者簡介：鄭臘年（1986-），2013年畢業(yè)于西南石油大學油氣田開發(fā)工程專業(yè)，現(xiàn)在長慶油田蘇里格氣田研究中心工作，從事氣藏工程方面的研究工作，郵箱：zln.0404@163.com。