天然氣高速非達西滲流動態(tài)產(chǎn)能計算

崔傳智，劉慧卿，耿正玲，郭 奇，呂恒宇

(1.中國石油大學，北京 102249;2.中國石油大學，山東 青島 266555)

天然氣高速非達西滲流動態(tài)產(chǎn)能計算

崔傳智1，劉慧卿1，耿正玲2，郭 奇2，呂恒宇2

(1.中國石油大學，北京 102249;2.中國石油大學，山東 青島 266555)

真實氣體的地下滲流是個不穩(wěn)定的非線性滲流過程，氣井產(chǎn)能受流動規(guī)律的影響，并隨地層壓力及PVT參數(shù)的變化而變化。綜合考慮天然氣高速非達西滲流特征和PVT參數(shù)動態(tài)特征，建立了外邊界封閉氣藏中平面徑向不穩(wěn)定滲流數(shù)學模型。計算結果表明，考慮高速非達西流和氣體PVT參數(shù)變化的產(chǎn)能高于達西線性流和PVT參數(shù)不變時的產(chǎn)能;計算得到了不同地層壓力下不同投產(chǎn)時間的動態(tài)IPR曲線，避免了靜態(tài)IPR曲線中許多不合理的假設，更加符合生產(chǎn)實際。

高速非達西流;動態(tài)IPR曲線;數(shù)學模型;真實氣體PVT參數(shù);氣井產(chǎn)能

引 言

真實氣體的地下滲流是個不穩(wěn)定的非線性滲流過程，在氣藏開發(fā)中，地層壓力隨時間發(fā)生變化，氣體的高壓物性參數(shù)隨壓力發(fā)生變化，氣體按非達西滲流規(guī)律流動，滲流過程較為復雜。目前在國內(nèi)外學者對氣井產(chǎn)能的研究中，部分學者是基于氣體流動滿足達西線性定律［1－5］，從穩(wěn)定滲流和不穩(wěn)定滲流入手，建立氣井的產(chǎn)能預測模型，部分學者考慮了氣體高速非達西滲流特征［6－9］。在以上研究中都忽略了真實氣體的PVT變化，假定氣體的高壓物性參數(shù)是常數(shù)。本文建立了考慮真實氣體PVT參數(shù)隨壓力變化的高速非達西滲流單井模型，采用數(shù)值方法進行求解，分析了真實氣體PVT參數(shù)和高速非達西流對產(chǎn)能的影響，計算了外邊界封閉條件下的動態(tài)IPR曲線，對氣藏產(chǎn)能計算具有指導意義。

1 真實氣體高速非達西滲流數(shù)學模型

基本假設:氣體滲流為等溫過程;天然氣為真實氣體，黏度、壓縮因子、壓縮系數(shù)等隨壓力變化;氣體流動為高速非達西滲流;考慮巖石彈性。

1.1 運動方程

氣體的高速非達西滲流用Forcheimer的二次方程來描述［10］:

式中:p為地層壓力，MPa;v為滲流速度，cm/s;K為滲透率，μm2;ρg為天然氣密度，g/cm3;μ(p)為天然氣黏度，是壓力的函數(shù)，mPa·s;β=7.644×1010/K3/2，為高速速度系數(shù)，cm－1［4，10－11］。

1.2 狀態(tài)方程

式中:γg為天然氣相對密度;Cg(p)為天然氣壓縮系數(shù)，是壓力的函數(shù)，MPa－1;Cf為巖石壓縮系數(shù)，MPa－1;z為天然氣壓縮因子;T為地層溫度，K;Ma為空氣分子量，g/mol;R為通用氣體常數(shù);φ為孔隙度;p0為大氣壓力，MPa;φ0為壓力p0對應的孔隙度。

1.3 天然氣連續(xù)性方程

式(7)是真實氣體不穩(wěn)定等溫滲流的綜合微分方程，式中考慮了非達西流動以及天然氣黏度、壓縮因子、壓縮系數(shù)等隨壓力的變化，是高度非線性的不穩(wěn)定滲流方程。

2 天然氣PVT參數(shù)計算模型

2.1 壓縮因子計算模型

天然氣壓縮因子計算采用Dranchuk等人通過擬合SK圖版提出的方程［5］:

式中:ρpr為中間變量，ρpr=0.27ppr/(ZTpr);ppr為對應壓力;Tpr為對應溫度;A1～A8為常數(shù)。

2.2 天然氣黏度計算模型

3 壓力和擬壓力的換算

壓力函數(shù)的計算有解析法、半解析法和數(shù)值積分法。本文采用了數(shù)值積分方法計算擬壓力，其公式為:

4 模型求解

考慮到井底附近壓力梯度大的特點，區(qū)域離散采用等對數(shù)步長，并將一維問題的柱坐標轉換成直角坐標。

式(12)為三對角矩陣方程，根據(jù)內(nèi)外邊界條件，采用追趕法解此方程組。

5 高速非達西流和氣體PVT參數(shù)變化對產(chǎn)能影響

已知一天然氣藏，外邊界封閉，γg=0.6，ppc=4.62 MPa，Tpc=200K，T=395K，K=0.5 μm2，巖石壓縮系數(shù) Cf=3.0×10－4MPa－1。計算出壓縮因子、黏度、壓縮系數(shù)和擬壓力，研究投產(chǎn)時地層壓力為30 MPa，投產(chǎn)1 d時氣體產(chǎn)能變化特征。

5.1 真實氣體PVT參數(shù)對產(chǎn)能影響

在達西線性流條件下計算了真實氣體PVT參數(shù)變化對產(chǎn)能的影響(圖1)。在井底流壓接近原始地層壓力時，各參數(shù)變化對產(chǎn)能影響較小;隨井底流壓的降低，產(chǎn)能差異增加。相同井底流壓下，PVT參數(shù)不變時的產(chǎn)能最低，PVT參數(shù)綜合變化的產(chǎn)能最高;黏度變化對產(chǎn)能的影響高于壓縮系數(shù)和壓縮因子變化對產(chǎn)能的影響。因此假定氣體PVT參數(shù)為常數(shù)時的產(chǎn)能要小于實際產(chǎn)能。

5.2 高速非達西流對產(chǎn)能影響

圖2為PVT參數(shù)變化和PVT參數(shù)不變2種情況下達西流和非達西流對應的IPR曲線，可以看出PVT參數(shù)變化和PVT參數(shù)不變2種情況下高速非達西流的產(chǎn)能高于達西線性流時的產(chǎn)能;隨井底流壓降低，產(chǎn)量差異變大。

圖2 達西流與高速非達西流情況下的產(chǎn)能曲線

6 真實氣體高速非達西流下的動態(tài)產(chǎn)能

利用編制的計算機程序可以計算真實氣體高速非達西滲流情況下不同地層壓力、不同井底流壓的產(chǎn)量隨時間變化曲線。圖3中給出了氣井在地層壓力分別30、20、10 MPa下開始生產(chǎn)，時間分別為1、10、100 d時的產(chǎn)量與流壓的關系曲線。可以看出，相同地層壓力下開始生產(chǎn)，時間不同，IPR曲線不同，隨時間增加，產(chǎn)量降低;不同地層壓力下投產(chǎn)，氣井的IPR曲線不同，隨投產(chǎn)時地層壓力的降低，產(chǎn)量降低。因此保持地層能量可以相對保持氣井獲得較高的產(chǎn)量。

圖3 不同地層壓力投產(chǎn)氣井的動態(tài)IPR曲線

從不同地層壓力投產(chǎn)氣井的動態(tài)無阻流量可以看出(圖4)，隨氣井投產(chǎn)時地層壓力降低，氣井的無阻流量降低;在一定地層壓力投產(chǎn)后隨生產(chǎn)時間的增加，氣井無阻流量逐漸降低。

圖4 不同地層壓力投產(chǎn)的動態(tài)無阻流量曲線

7 結論

(1)建立了考慮真實氣體PVT參數(shù)隨壓力變化的高速非達西不穩(wěn)定滲流數(shù)學模型，采用數(shù)值差分技術對圓形封閉地層中心一口井定井底壓力生產(chǎn)的數(shù)學模型進行了求解，編制了計算程序。可計算不同地層壓力下投產(chǎn)氣井在不同井底流壓不同時間的產(chǎn)量。

(2)考慮氣體PVT參數(shù)隨壓力變化和高速非達西流的產(chǎn)能高于PVT不變和達西線性流時的產(chǎn)能，并隨井底流壓降低，差異變大。

(3)氣井產(chǎn)能受投產(chǎn)時刻的地層壓力影響，地層壓力降增大時，氣井的產(chǎn)能和無阻流量降低;在一定地層壓力投產(chǎn)后，隨生產(chǎn)時間的增加，氣井的產(chǎn)能和無阻流量逐漸降低。

(4)本文方法所建立的動態(tài)IPR曲線，考慮了氣體PVT參數(shù)變化和高速非達西流特征，避免了靜態(tài)IPR曲線中許多不合理的假設，更加符合生產(chǎn)實際。

［1］郭建春，羅天雨，趙金洲，等.氣井壓后穩(wěn)態(tài)產(chǎn)能模型及求取無阻流量的變系數(shù)一點法［J］.特種油氣藏，2005，12(2):52 －54.

［2］黎洪，彭蘇萍.高壓氣井產(chǎn)能評價方法研究［J］.石油勘探與開發(fā)，2001，28(6):77 －79.

［3］李曉平，胡勇.氣水同產(chǎn)井瞬態(tài)流入動態(tài)關系曲線探討［J］.天然氣工業(yè)，2001，21(3):65－67.

［4］楊景海，何風珍，鄒艷華.利用不穩(wěn)定試井資料確定氣井動態(tài)產(chǎn)能方法研究［J］.天然氣工業(yè)，2007，27(3):95－96.

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

［6］康曉東，李相方，郝偉.氣井高速非達西流動附加壓降計算公式的修正［J］.油氣井測試，2004，13(5):4－5.

［7］張烈輝，朱水橋，王坤 等.高速氣體非達西滲流數(shù)學模型［J］.新疆石油地質，2004，25(2):165－167.

［8］李傳亮，姚淑影.氣井試井分析中氣體物性參數(shù)使用原始物性參數(shù)之探討［J］.特種油氣藏，2010，17(5):123－124.

［9］ Zeng Fanhua，Zhao Gang.Gas well production analysis with non － Darcy flow and real－ gas PVT behavior［J］．Journal of Petroleum Science and Engineering，2007，59(3):169－182.

［10］Huang H，Ayoub J.Applicability of the forchheimer equation for Non － Darcy Flow in porous media［J］.SPE Journal，2008，13(1):112 －122.

［11］生如巖，李相方.一種考慮紊流影響的產(chǎn)水氣井開采動態(tài)預測模型［J］.中國海上油氣，2004，16(5):328－330.

Calculation of dynamic productivity of natural gas in high velocity non－Darcy flow

CUI Chuan– zhi1，LIU Hui－qing1，GENG Zheng– ling2，GUO Qi2，LV Heng– yu2
(1．China University of Petroleum，Beijing 102249，China;
2．China University of Petroleum，Qingdao，Shandong 266555，China)

The actual flow of gas underground is an unsteady nonlinear process and the productivity of gas wells is affected by the flow patterns and it changes with the change of formation pressure and PVT parameters．In this paper，a mathematical model for the two dimensional radial unsteady flow in external boundary closed gas reservoirs was established considering the high velocity non－Darcy flow characteristics of natural gas and the dynamic features of PVT parameters．The calculation results indicated that the productivity considering the high velocity non－Darcy flow and the changes of PVT was higher than that of Darcy linear flow with constant PVT parameters．Dynamic IPR curves at different times under different formation pressures were also obtained from the calculations，eliminating unreasonable assumptions in the static IPR curves and conforming better to the actual production．

high velocity non－Darcy flow;dynamic IPR curve;mathematical model;real gas PVT parameter;gas well productivity

TE312

A

1006－6535(2011)06－0080－03

20110524;改回日期20110702

國家重點基礎研究發(fā)展計劃“973”計劃“溫室氣體提高石油采收率的資源化利用及地下埋存”(2006CB705804)

崔傳智(1970－)，男，1996年畢業(yè)于石油大學(華東)油氣田開發(fā)工程專業(yè)，獲碩士學位，2005年畢業(yè)于中國地質大學(北京)礦產(chǎn)普查與勘探專業(yè)，獲博士學位，現(xiàn)為中國石油大學(北京)在站博士后，從事油氣田開發(fā)工程科研工作。

編輯 孟凡勤