考慮變質(zhì)量流影響的煤層氣水平井數(shù)值模擬研究

張先敏，馮其紅，張紀遠

(中國石油大學，山東 青島 266580)

引 言

與常規(guī)壓裂直井相比，采用多分支水平井技術(shù)開采低滲透煤層氣具有極大的優(yōu)越性［1］。目前，國內(nèi)外學者已在煤層氣水平井開采模擬方面開展了大量的研究工作:Sung和Ertekin［2］建立了可模擬多個水平井眼的煤層氣排采模型，Maricic等［3］開展了煤層氣多分支結(jié)構(gòu)與羽狀分支結(jié)構(gòu)的水平井數(shù)值模擬研究，高德利等［4］建立了不同流態(tài)下分支井眼臨界長度與直徑計算模型以及多分支井身結(jié)構(gòu)設計模型，Dikken［5］在1989年提出水平井筒內(nèi)的壓降不能忽略，張冬麗等［6］建立了考慮摩擦壓降和加速度壓降影響的煤層氣羽狀水平井開采數(shù)學模型，Keim等［7］對煤層氣多分支井的分支間距、長度以及井眼軌跡進行了優(yōu)化研究。

本文在此基礎上，充分考慮羽狀水平井筒中分支和主支流體匯合的流體流動特點，提出了煤層氣羽狀水平井開采數(shù)學模型方法，進行了相應的數(shù)值模擬分析研究。

1 數(shù)學模型的建立及求解

假設水平井筒中的流體流動為一維、穩(wěn)態(tài)流動，分支流體匯流入主支井筒的過程近似為等溫流動。將煤層氣羽狀水平井的主支井筒從跟端到指端劃分成N個微元井段，沿井筒延伸方向各微元井段長為Δxi;將第j分支井筒從跟端至指端分成Mj段，沿井筒延伸方向各微元井段長為 Δxj，k。對微元井段在主支上且該井段內(nèi)有分支的情況，在主支井筒上取流動微元控制體見圖1。對此匯合點處可簡化為斜三通，根據(jù)伯努利方程［8－9］可得:式中:pwi+1、pwi分別為主支相鄰的上、下游井段中心的壓力，MPa;pwj，1為微元控制體內(nèi)分支井段的壓力，MPa;D為主支井筒直徑，m;dj為第j分支井筒直徑，m;vRi為匯合處分支流體的入流速度，m/s;vi+1、vi分別為流體流入和流出該微元井段的平均流速，m/s為第j分支流入主支的平均入流速度，m/s;vi為該井段主支流體速度的平均值，m/s;ρi為流體的密度，kg/m3;fi+1為i+1主支井段的摩擦系數(shù)，采用 Ouyang［10］的摩擦系數(shù)公式計算;fj，1為第 j分支的摩擦系數(shù);Δxi為主支微元井段長度，m;Δxj，1為第 j分支第 1 井段長度，m;ξi為主支微元井段局部阻力系數(shù);ξj，1為第j分支第1井段局部阻力系數(shù);αi+1、αi、αj，1依次為上游井段、下游井段、第 j分支第1井段動能修正系數(shù)，在實際應用中大多數(shù)情況下令 αi+1≈αi≈αj，1= α，紊流情況下 α 值近似取為1，層流情況下α值取為2。

圖1 井筒流體匯合流動微元控制體示意圖

根據(jù)質(zhì)量守恒原理可知，流入與流出微元井段的質(zhì)量流量速率之差等于該微元井段內(nèi)質(zhì)量增加速率，即:

式中:A為主支井筒截面積，m2;ARj為第j分支井筒截面積，m2;νVi為煤層徑向流入主支微元井段的平均速度，m/s。

將式(3)代入式(1)整理可得主支上、下游井段之間的壓降方程為:

式中:qi為煤層流體流入主支微元井段的流量，m3;QRi為微元井段內(nèi)分支流入主支的流量，m3;Qi為相鄰上游段流入主支井段的流量，m3;fi為主支井段內(nèi)流體和管壁之間的摩擦系數(shù)。

同理，可得到羽狀水平井分支井段與主支井段之間的壓降方程為:

與裸眼完井水平井筒模型的方程［11］相似，可 得到無匯合點的主支井段或分支井段壓降方程為:

式中:pwj，k+1為第j分支第 k+1井段的壓力，MPa;pwj，k為第 j分支第 k 井段的壓力，MPa;Δxj，k為第 j分支第k井段長度，m;ρj，k為第j分支第k井段內(nèi)混合流體的密度，kg/m3;fj，k為第j分支第k井段內(nèi)流體和壁面之間的摩擦系數(shù);QRj，k為從第j分支第k井段相鄰上游段流入該井段的流量，m3;qj，k為煤層徑向流入第j分支第k井段的流量，m3。

根據(jù)壓力連續(xù)原理，煤層流體與井筒流體在井壁處的壓力應相等，結(jié)合井筒變質(zhì)量流動壓降方程以及煤儲層氣—水兩相流動方程［12］，加上相應的初邊界條件，即可耦合構(gòu)成1個完整的煤層氣羽狀水平井開采數(shù)學模型。在定井底流壓條件下，首先令井筒壓力為已知井底壓力;若在定產(chǎn)條件下，則可設為任意合理的井筒壓力值。采用全隱式方法求解煤儲層流動模型可得到煤層壓力和含水飽和度值，然后代入井筒壓降方程組中求解得到井筒壓力值;再將最新的井筒壓力值帶入煤儲層流動模型中重復計算，如此反復迭代，直到滿足一定的收斂條件，即可得到煤層壓力、含水飽和度和各井段的壓力分布值。

2 計算實例

沁水盆地南部是中國煤層氣勘探開發(fā)最有利的區(qū)域之一，區(qū)內(nèi)廣泛發(fā)育的上石炭統(tǒng)太原組15號煤層和下二疊統(tǒng)山西組3號煤層是主要開發(fā)目的層，其中3號煤層厚度較大，橫向分布穩(wěn)定，煤層頂?shù)装宸馍w能力較強，可作為煤層氣羽狀水平井開發(fā)的目的層。煤層基本參數(shù)如下:地層壓力為2.45 MPa，含氣量為 19.60 m3/t，地層溫度為24.2℃，煤層厚度為6.0m，滲透率為0.61×10－3μm ，孔隙度為 0.06，Langmuir壓力為 3.58 MPa，Langmuir體積為48.69 m3/t。以煤層氣羽狀水平井進行模擬開采，主支長為1200 m，分支長為800 m，同側(cè)分支間距為300 m，相鄰不同側(cè)分支間距為50 m，井段總長5700 m，分支與主支夾角為45°。

主支井筒不同位置單位長度徑向入流量沿井長分布是不均勻的，反映出主支井筒附近煤層氣解吸程度的差異，可按照產(chǎn)量上升和產(chǎn)量遞減2個階段來分別分析。在產(chǎn)量上升階段 (圖2a)，主支井筒的單位長度入流量分布曲線呈上凸型，分支與主支匯合點附近煤層降壓程度較高，氣體徑向入流量相對較大，主支井筒兩端部的入流量較低;在產(chǎn)量遞減階段 (圖2b)，主支井筒入流量分布曲線總體上呈下凹型，反映出主支不同位置處供氣范圍的不同。由于中部井段供氣范圍較小以及各分支之間以及與主支之間的壓力干擾，煤層氣解吸程度較高，因此入流量越來越低，尤其分支節(jié)點處的入流量更低;相反，由于主支井筒兩端部的供氣范圍較大，井段氣體入流量保持相對較高水平。

圖2 不同時刻主支井筒氣體入流量分布曲線

圖3為不同時刻主支井筒壓力分布圖，可以看出，沿主支井筒壓力分布是不均勻的。越靠近主支井筒跟端，井筒質(zhì)量流量越大，摩擦引起的壓力損失以及入流、匯流造成的壓力損失亦越大，主支井筒壓力越低，因此井筒壓降在氣井產(chǎn)能動態(tài)評價中不可忽略。另外，隨著煤層氣不斷開采，主支井筒指端與跟端之間的壓差呈先增加后減少的趨勢，這是由于隨著煤層氣體的不斷解吸，流體(氣和水)沿井壁徑向入流和分支流體匯入的量增大，摩擦、動量變化以及匯合引起的壓降損失也逐漸增大(圖3a);隨著生產(chǎn)時間的延長，主支徑向入流量和分支匯入量不斷減小，且此時主支井筒內(nèi)流體流動處于紊流光滑區(qū)，井筒內(nèi)的壓降損失也不斷減少，井筒壓力越來越接近出口端壓力(圖3b)。

3 結(jié)論

(1)根據(jù)煤層氣羽狀水平井復雜井筒結(jié)構(gòu)以及井筒內(nèi)流體流動特征，建立了考慮井筒變質(zhì)量流影響的煤層氣羽狀水平井開采數(shù)學模型。

(2)沁水盆地煤層氣羽狀分支水平井開采模擬分析表明:由于井筒變質(zhì)量流的影響，不同時刻的主支井筒單位長度入流量分布曲線分別呈上凸型和下凹型，且水平井筒指端與跟端之間的壓差呈先增加后減少的趨勢。

圖3 不同時刻主支井筒壓力分布曲線

［1］Palmer I.Coalbed methane completions:a world view［J］.International journal of coal geology，2010，8(2):184 －195.

［2］Sung W，Ertekin T，Schwerer F C.The Development，testing and application of a comprehensive coal seam degasification model［C］.SPE15247，1986:457 － 463.

［3］Maricic N，Mohaghegh S D，Artun E A parametric study on the benefits of drilling horizontal and multilateral wells in coalbed methane reservoirs［C］.SPE96018，2005:1 －8.

［4］高德利，鮮保安.煤層氣多分支井身結(jié)構(gòu)設計模型研究［J］.石油學報，2007，28(6):113－117.

［5］Dikken B J.Pressure drop in horizontal wells and its effect on production performance［J］.JPT，1990，42(11):1426－1433.

［6］張冬麗，王新海.煤層氣羽狀水平井開采數(shù)值模擬研究［J］.煤田地質(zhì)與勘探，2005，33(4):47－51.

［7］Keim S A，Luxbacher K D，Karmis M.A numerical study on optimization of multilateral horizontal wellbore patterns for coalbed methane production in Southern Shanxi Province，China［J］.International journal of coal geology，2011，86(4):306 －317.

［8］黃世軍，程林松，李秀生，等.多分支水平井壓力系統(tǒng)分析模型［J］.石油學報，2003，24(6):81－86.

［9］步玉環(huán)，王瑞和，劉繼林，等.分支井匯合流動壓降計算［J］.鉆采工藝，2007，30(4):10－12.

［10］Ouyang L B，Aziz K.A mechanistic model for gas－liquid flow in pipes with radial influx or outflux［C］.SPE56525，1999:359 －372.

［11］劉想平，張兆順，劉翔鶚，等.水平井筒內(nèi)與滲流耦合的流動壓降計算模型［J］.西南石油學院學報，2000，22(2):36－39.

［12］張先敏，同登科，胡愛梅.沁水盆地煤層氣井網(wǎng)開采數(shù)值模擬［J］.工程力學，2008，25(12):218－222.