多段壓裂水平氣井紊流產(chǎn)能模擬模型

            ——以塔里木盆地克拉蘇氣田大北區(qū)塊為例

摘　要：模擬預(yù)測(cè)壓裂水平氣井產(chǎn)量模型通常未考慮高速非達(dá)西流動(dòng)的影響，特別是進(jìn)行了多段壓裂的水平氣井的產(chǎn)能計(jì)算更未涉及。為此，將裂縫劃分為若干單元，應(yīng)用點(diǎn)源理論、氣體不穩(wěn)定滲流公式及勢(shì)疊加原理，考慮N條裂縫同時(shí)產(chǎn)生相互干擾條件t時(shí)刻地層中任一點(diǎn)產(chǎn)生的總壓降計(jì)算模型，以及高速非達(dá)西流動(dòng)條件對(duì)壓裂水平氣井壓后產(chǎn)能的影響，建立了更為完善的模擬預(yù)測(cè)壓裂水平氣井產(chǎn)能的新模型。實(shí)例模擬計(jì)算結(jié)果表明：壓裂水平氣井的產(chǎn)量隨著生產(chǎn)時(shí)間的增加而逐漸降低且降低的幅度也不斷減??；水平氣井端部裂縫的產(chǎn)量最高，中部裂縫產(chǎn)量次之，中心裂縫產(chǎn)量最低；水平氣井壓后生產(chǎn)動(dòng)態(tài)模擬預(yù)測(cè)新模型為優(yōu)選高產(chǎn)水平氣井橫向裂縫數(shù)目和裂縫參數(shù)、確定合理生產(chǎn)壓差提供了更準(zhǔn)確的手段。

關(guān)鍵詞：水平井 氣井 多段壓裂 高速非達(dá)西流動(dòng)不穩(wěn)定滲流 生產(chǎn)模擬 數(shù)學(xué)模型

A simulation model for the turbulent production of a multi-stage fractured gas well：A case sudy from the Dabei block in Kelasu Gas Field，Tarim Basin

Abstract：Generally，the simulation models of predicting the production of fractured horizontal gas wells never take into account the influence of the high speed non-Darcy flow，especially in the case of the production calculation of a multi stage fractured horizontal gas well.Thus，the fractures were first divided into different groups，then the point source theory，the unstahle gas seepage formula and the potential superposition theory were adopted to establish the calculation model considering the total pressure drop of any point in the formation at the time t and under the mutual interferences simultaneously generated by the fractures with the number of N，and finally a better simulation roodel was built of predicting the production of a fractured horizontai gas well considering the impact of the high speed Darcy flow on a post-frac horizontal well．Case studies show that the production of a fractured horizontal well gradually goes down over time，and downward margin decreases too；the fractures at the top of the horizontal well stand out for their highest gas output，followed by those in the middle and the lowest by those in the center；this new simulation model is a more accurate tool to determine an optimized number of horizontal fractures and more suitable fracture parameters for a more reasonable drawdown pressure．

Key words：horizontal well，gas well，multi stage fracturing，high speed non Darcy flow，unstable seepage，production simulation，mathematical model

水平井是有效開發(fā)低滲透儲(chǔ)層的重要技術(shù)，尤其是完井與壓裂工藝技術(shù)的完美結(jié)合開創(chuàng)了低滲透、極低滲透儲(chǔ)層高效開發(fā)的新紀(jì)元^[1]。壓裂水平井生產(chǎn)動(dòng)態(tài)模擬預(yù)測(cè)一直得到國內(nèi)外工程技術(shù)界人士的高度重視。國外率先對(duì)多段壓裂水平氣井的產(chǎn)量預(yù)測(cè)開展研究，主要是針對(duì)地層和裂縫中的達(dá)西滲流規(guī)律建立水平井多段壓裂生產(chǎn)模擬模型^[2-4]。歸結(jié)起來大致可以分為兩類：一類是基于穩(wěn)定滲流理論，把整條裂縫看成一個(gè)點(diǎn)匯或者線匯，然后利用勢(shì)疊加原理或積分變換建立壓裂水平井壓后產(chǎn)能解析模型^[5-9]；另一類是基于不穩(wěn)定滲流理論，將每條裂縫分成若干微元段，然后利用勢(shì)疊加原理或積分變換建立壓裂水平井壓后產(chǎn)能半解析模型^[10-12]。

我國西部克拉蘇氣田大北區(qū)塊儲(chǔ)層埋藏較深、壓裂施工壓力很高、水平井壓裂施工難度極大。借鑒水平井多段壓裂開發(fā)的核心思想，開發(fā)方案設(shè)計(jì)的壓裂水平氣井產(chǎn)量將達(dá)到甚至超過100×10⁴m³／d，高速流動(dòng)產(chǎn)生的紊流效應(yīng)對(duì)水力壓裂和氣井生產(chǎn)的影響必然不可忽略而必須高度重視。一方面用于優(yōu)化水平井裂縫條數(shù)和各條壓裂裂縫參數(shù)，另一方面用以確定壓裂水平氣井的開采工作制度。但在水平氣井多段壓裂紊流(高速非達(dá)西流動(dòng))下產(chǎn)能模擬方面幾乎仍是空白，前述兩類模型雖然都呵以考慮橫向多段壓裂裂縫之間的相互干擾、通過簡(jiǎn)化考慮非對(duì)稱兩翼人工裂縫影響，但都沒有考慮非達(dá)西因素對(duì)壓裂井產(chǎn)能的影響。筆者對(duì)此進(jìn)行探索以滿足實(shí)際工程的迫切需要。

1　壓裂水平氣井產(chǎn)能預(yù)測(cè)新模型

考慮壓裂水平井有N條縫垂直于水平井筒的橫向人工裂縫，建立圖1所示坐標(biāo)。當(dāng)壓裂形成奇數(shù)條裂縫時(shí)，以中間那條裂縫和水平井井筒相交的點(diǎn)作為原點(diǎn)；當(dāng)壓裂形成偶數(shù)條裂縫時(shí)，以中間兩條裂縫間距的平分線與水平井井筒的交點(diǎn)為原點(diǎn)。再由其他裂縫與該裂縫間的幾何位置間距，可以依次確定每條裂縫的坐標(biāo)。

如圖l所示，將裂縫的兩翼分別分成咒單元，每一單元視為一個(gè)點(diǎn)匯進(jìn)行處理。以第j單元的中心坐標(biāo)來表示左翼裂縫上的第j個(gè)點(diǎn)匯的坐標(biāo)(x1j、y1j)：

同理，其右翼裂縫上的第J個(gè)點(diǎn)匯的坐標(biāo)可以表示為：

式中xf1為裂縫左翼的長(zhǎng)度，m；xfr為裂縫右翼的長(zhǎng)度，m；yf為裂縫Y方向坐標(biāo)。

根據(jù)氣體不穩(wěn)定滲流公式^[12]以及勢(shì)疊加原理，N條裂縫同時(shí)生產(chǎn)時(shí)，在t時(shí)刻地層中任一點(diǎn)(x，y)產(chǎn)生的總壓降為：

式中pi為原始地層壓力，Pa；p(x，y，t)為地層中點(diǎn)(x，y)在t時(shí)刻的壓力，Pa；qflij為第i條裂縫左端第j段的產(chǎn)量，m³／s；xfli、xfri分別為第i條裂縫左、右翼的長(zhǎng)度，m；ps為地面標(biāo)況下的壓力，Pa；Z為偏差因子；m為氣體黏度，Pa·S；K為地層滲透率，m²；h為產(chǎn)層厚度，m；t為滲流時(shí)間，s；T為地層溫度，K；ts為地面標(biāo)準(zhǔn)溫度，K；h為導(dǎo)壓系數(shù)；yfi為第i條裂縫y方向坐標(biāo)。

由于水力裂縫的產(chǎn)生，使得地層壓力重新分布。使得地層壓力重新分布。假設(shè)第i條裂縫的左翼和右翼的縫端壓力分別為pflj和pfri。則t時(shí)刻可得第i條裂縫的左翼尖端產(chǎn)生的壓降為：

同理，則t時(shí)刻第i條裂縫右翼尖端產(chǎn)生的壓降為：

式中pfli、pfri分別為第i條裂縫左翼端部和右翼端部的壓力，Pa；如；xflk、xfrk；分別為第k條裂縫左翼和右翼的長(zhǎng)度，m；yfk為第k條裂縫在Y方向的坐標(biāo)。

如果壓裂裂縫關(guān)于井筒不對(duì)稱，進(jìn)行簡(jiǎn)化處理。取裂縫兩端的平均壓力作為裂縫尖端壓力，即

根據(jù)面積相等(xfli+xfri)h=Ri²，可得：

第i條裂縫就可以看成是一個(gè)以流動(dòng)半徑為Ri，厚度為裂縫寬度wi，邊界壓力為裂縫尖端壓力p(xfi,yfi，t)，滲透率為Kf，流壓pwf的微型氣藏，其流動(dòng)方式可以看成徑向流，所以氣體在裂縫內(nèi)的流動(dòng)過程可表示為：

對(duì)于高產(chǎn)氣井，地層和裂縫均可能產(chǎn)生高速非達(dá)西流動(dòng)，即紊流效應(yīng)。當(dāng)?shù)貙訅毫Σ惶咔疑a(chǎn)壓差不是很大的時(shí)候，忽略地層紊流效應(yīng)，認(rèn)為高速非達(dá)西流動(dòng)只發(fā)生于裂縫內(nèi)。那么，高速非達(dá)兩流動(dòng)產(chǎn)生的附加壓降也主要發(fā)生在裂縫內(nèi)，類似于表皮效應(yīng)的處理思路，可以用流量相關(guān)表皮系數(shù)(Dqfi)來表示，則可將式(6)改為：

其中

對(duì)于非達(dá)西影響系數(shù)的描述有很多種，筆者應(yīng)用Jones^[13]公式來描述非達(dá)西影響系數(shù)：

b＝2.0127×10¹¹/Kf^1.55            (8)

式中b為描述裂縫的紊流影響的系數(shù)，m^-1；rg為氣體相對(duì)密度。

2　實(shí)例計(jì)算

某氣藏壓裂層有效藏厚度8m，地層滲透率為10mD，地層溫度80℃，井筒半徑0.1m，水平井段長(zhǎng)度400m；原始地層壓力27.56MPa，井底流壓為23.56MPa；氣體黏度為0.022mPa·S，氣體壓縮因子為0.957。壓裂形成5條橫向裂縫，各條裂縫長(zhǎng)度200m，裂縫滲透率50D，支撐裂縫寬度0.003m。模擬計(jì)算結(jié)果如圖2、3所示。

由圖2可以看出：壓裂水平氣井的產(chǎn)量隨著生產(chǎn)時(shí)間的增加而逐漸降低，且降低的幅度也不斷減?。槐M管壓裂水平氣井壓后初期產(chǎn)量高于后期產(chǎn)量，但非達(dá)西流動(dòng)對(duì)壓裂水平氣井產(chǎn)量的影響非常顯著。這與單裂縫考慮非達(dá)西因素的計(jì)算結(jié)果一致；即產(chǎn)量越高，非達(dá)西因素對(duì)壓裂水平氣井產(chǎn)能影響越明顯。

由圖3可以看出：水平氣井端部裂縫(1號(hào)和5號(hào)裂縫)的產(chǎn)量最高，中部裂縫(2號(hào)和4號(hào)裂縫)產(chǎn)量其次，中心裂縫(3號(hào)裂縫)產(chǎn)量最低。這是由于多裂縫之間干擾相互疊加影響的結(jié)果，水平氣井端部裂縫受到的干擾最小，而中心裂縫受到的干擾最大。與不考慮高速紊流(非達(dá)西流動(dòng))時(shí)的產(chǎn)量變化規(guī)律相似。

3　結(jié)論與認(rèn)識(shí)

1)基于不穩(wěn)定滲流，考慮了水平氣井橫向多裂縫干擾和高速非達(dá)西流動(dòng)對(duì)壓裂水平氣井產(chǎn)能的影響，建立了模擬預(yù)測(cè)壓裂水平氣井產(chǎn)能的新模型。該模型更準(zhǔn)確地反映壓裂水平氣井的生產(chǎn)特性。

2)壓裂水平氣井的兩端裂縫的產(chǎn)量最大，中間裂縫的產(chǎn)量最低。但同時(shí)非達(dá)西因素對(duì)兩端裂縫比對(duì)中間裂縫產(chǎn)量的影響更顯著。

參考文獻(xiàn)

[1]STANOJCIC M，RISPLER K．How to achieve and control branch fracturing for unconventional reservoirs：Two novel multistage-stimulation processes[C]∥paper l36566 presented at the Canadian Unconventional Resources＆International Petroleum Conference，19-21 October 2010，Calgary，Alberta，Canada．New York：SPE，2010．

[2]ROBERTS B E，ENGEN H V．Productivity of muhiply fractured horizontal wells in tight gas reservoirs[C]//paper 23113 presented at the Offshore Europe Conference，3-6 September l991，Aberdeen，United Kingdom．New York：SPE，l991．

[3]SOLIMAN M Y，AZARI M，HUNT J L，et al．Design and analysis of fractured horizontal wells in gas reservoirs[C]//paper 35343 presented at the International Petroleum Conference and Exhibition，5-7 March l996，Villahermosa，Mexico．New York：SPE，1996．

[4]BAGHERIAN B，SARMADIVALEH M，GHALAMBOR A，et al．Optimization of multiple-fractured horizontal tight gas well[C]//paper 27899 presented at the 2010 SPE International Symposium and Exhibition on Formation Damage Control，10-12 February 2010，Lafayette，Louisiana，USA．New York：SPE，2010．

[5]郎兆新，張麗華，程林松．壓裂水平井產(chǎn)能研究[J]．石油大學(xué)學(xué)報(bào)：自然科學(xué)版，l994，18(2)：43-46．

Lang Zhaoxin，Zhang Lihua，Cheng Linsong．Investigation on pkodectivity of fracturen horizontal well[J]．Journal of the University of Petroleum：Natural Science Edition，1994，18(2)：43-46．

[6]寧正福，韓樹剛，程林松，等．低滲透油氣藏壓裂水平井產(chǎn)能計(jì)算方法[J]．石油學(xué)報(bào)，2002，23(2)：68-71．

Ning Zhengfu,Han Shugang,Cheng Linsong,et al.Productivity calculation method of fractured horizontal wells in low permeability oil or gas field[J].Acta Petrolei Sinica,2002,23(2):68-71.

[7]韓樹剛，程林松，寧正福.氣藏壓裂水平井產(chǎn)能預(yù)測(cè)新方法[J].石油大學(xué)學(xué)報(bào)：自然科學(xué)版，2002，26(4):36-39.

Han Shugang，Cheng Linsong，Ning Zhengfu．A new method for predicting productivity of fractured horizontal wells in gas reservoir[J]．Journal of the University of Petroleum：Natural Science Edition，2002，26(4)：36-39．

[8]岳建偉，段永剛，陳偉，等．含多條垂直裂縫的壓裂氣井產(chǎn)能研究[J]．大慶石油地質(zhì)與開發(fā)，2004，23(3)：46-48．

Yue Jianwei，Duan Yonggang，Chen Wei，et al．Study on productivity performances of horizontal fractured gas well with many vertical fractures[J]．Petroleum Geology＆．Oil-field Development in Daqin9，2004，23(3)：46-48．

[9]孫福街，韓樹剛，程林松，等．低滲氣藏壓裂水平井滲流與井筒管流耦合模型[J]．西南石油學(xué)院學(xué)報(bào)，2005，27(1)：32-36．

Sun Fujie，Han Shugang，Cheng Linsong，et al．Low permeability gas reservoir fracturing in horizontal well with well bore seepage pipe coupling model[J]．Journal of Southwest Petroleum Institute，2005，27(1)：32-36．

[10]徐嚴(yán)波，齊桃，楊風(fēng)波，等．壓裂后水平井產(chǎn)能預(yù)測(cè)新模型[J]．石油學(xué)報(bào)．2006，27(1)：89-91．

Xu Yanbo，Qi Tao，Yang Fengbo，et al．New model for productivity test of horizontal well after hydraulic fracturing[J]．Acta Petrolei Sinica，2006，27(1)：89-91．

[11]雷征東，李相方，鄭紅軍．基于不穩(wěn)定滲流壓裂水平氣井產(chǎn)能研究[J]．天然氣工業(yè)，2006，26(4)：102-104．

Lei Zhengdong，Li Xiangfang，Zheng Hongjun．A study of deliverability of fractured horizontal gas well based on unsteady state flow[J]．Natural Gas Industry，2006，26(4)：102-104．

[12]葛家理，劉月田，姚約東．現(xiàn)代油藏滲流力學(xué)原理(上)[M]．北京：石油工業(yè)出版社，2003：167-175．

Ge Jiali，Liu Yuetian，Yao Yuedong．Modern reservoir seepage mechanics：Part I[M]．Beijing：Petroleum Industry Press，2003：l67-175．

[13]JONES S C．Using the inertial coefficient，b，to characterize heterogeneity in reservoir rock[C]//paper l6949 presented at the 62nd Annual Technical Conference and Exhibition of the Society of Petroleum Engineers，September 27-30．1987，Dallas．Texas，USA．New York：SPE．1987．

本文作者：胡永全  嚴(yán)向陽  趙金洲

作者單位：“油氣藏地質(zhì)及開發(fā)工程”國家重點(diǎn)實(shí)驗(yàn)室·西南石油大學(xué)

西南石油大學(xué)研究生院