考慮裂縫干擾的氣藏壓裂水平井產(chǎn)能預(yù)測(cè)模型

王樹平，史云清，嚴(yán) 謹(jǐn)，鄭榮臣

（中石化石油勘探開發(fā)研究院，北京 100083）

考慮裂縫干擾的氣藏壓裂水平井產(chǎn)能預(yù)測(cè)模型

王樹平，史云清，嚴(yán) 謹(jǐn)，鄭榮臣

（中石化石油勘探開發(fā)研究院，北京 100083）

基于壓裂水平井井身結(jié)構(gòu)復(fù)雜，需要建立考慮裂縫參數(shù)變化且符合水平井滲流特征的產(chǎn)能預(yù)測(cè)模型.為準(zhǔn)確預(yù)測(cè)壓裂水平井產(chǎn)能，通過保角變換和勢(shì)的疊加原理，建立新模型，不僅考慮裂縫長度、裂縫間距、裂縫夾角的變化對(duì)壓裂水平井產(chǎn)能的影響，還考慮裂縫間的相互干擾.結(jié)果表明：壓裂水平井產(chǎn)能隨著裂縫條數(shù)、長度和裂縫夾角的增加而增大.模型通過現(xiàn)場(chǎng)驗(yàn)證，計(jì)算單井無阻流量與試采結(jié)果對(duì)比，相對(duì)誤差為8.3%，計(jì)算結(jié)果較為準(zhǔn)確.

壓裂水平井；產(chǎn)能預(yù)測(cè)；計(jì)算模型；裂縫干擾；低滲氣藏

0 引言

近年來，國內(nèi)外許多學(xué)者對(duì)壓裂水平井產(chǎn)能預(yù)測(cè)開展相關(guān)研究.Giger F M等［1］首次利用數(shù)學(xué)模型分析壓裂水平井產(chǎn)能，但其模型不能很好地耦合流體在地層和裂縫中的流動(dòng)；Karcher R J等［2］和Soliman M Y等［3］修正了Giger模型，利用數(shù)值模擬方法，考慮橫向裂縫和縱向裂縫時(shí)，提出了定壓條件下致密氣藏壓裂水平井產(chǎn)能模型，但沒有考慮裂縫參數(shù)的變化；Joshi S D［4］建立了一個(gè)壓裂水平井解析模型，但裂縫假設(shè)為無限導(dǎo)流能力；Raghavan R和Joshi S D［5］建立了預(yù)測(cè)多條垂直裂縫的壓裂水平井產(chǎn)能模型，模型中利用有效井徑的概念模擬流體向裂縫井的流動(dòng)，忽略了流體在裂縫內(nèi)的流動(dòng)階段；Guo B和Schechter D S［6］、Wei Y和Economides M J［7］及Guo Boyun，Yu Xiance［8］建立的壓裂水平井模型不僅可以考慮油藏的徑向流和線性流，還考慮裂縫線性流和裂縫內(nèi)的徑向流，在水平井近井地帶滲流特征的刻畫上比較完善.20世紀(jì)90年代以來，我國許多學(xué)者對(duì)水平井壓裂產(chǎn)能預(yù)測(cè)及其影響因素進(jìn)行多方面、不同層次的分析研究，早期郎兆新、范子菲等［9－12］建立油藏壓裂水平井產(chǎn)能預(yù)測(cè)模型.在此基礎(chǔ)上，韓樹剛等［13］提出氣藏壓裂水平井產(chǎn)能預(yù)測(cè)模型，模型考慮地層、裂縫和井筒的三者耦合；孫福街等［14］建立考慮加速和摩擦損失的低滲氣藏壓裂水平井產(chǎn)能方程；徐嚴(yán)波等［15］建立油藏壓裂水平井，考慮各條裂縫的長度、縫寬、導(dǎo)流能力、縫間距以及裂縫平面與水平井筒之間的夾角對(duì)各條裂縫產(chǎn)量的影響.目前已有模型在滲流特征和裂縫參數(shù)方面考慮比較全面，但還沒有在考慮流體在地層和裂縫中流動(dòng)的同時(shí)，考慮裂縫參數(shù)的變化，尤其是裂縫間干擾對(duì)產(chǎn)能的影響.筆者利用保角變換和勢(shì)的疊加原理，建立一個(gè)新的壓裂水平井產(chǎn)能預(yù)測(cè)模型，模型不僅慮裂縫半長、裂縫間距和裂縫夾角的變化，還考慮裂縫間相互干擾對(duì)壓裂水平井產(chǎn)能的影響，從而達(dá)到準(zhǔn)確預(yù)測(cè)氣藏壓裂水平井產(chǎn)能的目的，進(jìn)而為水平井優(yōu)化設(shè)計(jì)、氣井合理配產(chǎn)和開發(fā)方案經(jīng)濟(jì)評(píng)價(jià)提供依據(jù).

1 預(yù)測(cè)模型

低滲氣藏中（川西中淺層和華北）沉積環(huán)境以河流相為主，縱向上主力層隔層穩(wěn)定.河道中一口壓裂水平井見圖1.圖1（a）中，Re為泄氣半徑；L1，L2，L3均為裂縫間距；a為河道寬度；c為保角變換后的河道寬度；d為保角變換后的裂縫距泄氣邊界距離；qfsc1，qfsc2，qfsc3，qfsc4均為單條裂縫產(chǎn)量.假設(shè)裂縫長度、裂縫夾角和裂縫間距各不相同，儲(chǔ)層泄氣面積為條帶狀；平行水平井兩側(cè)邊界及儲(chǔ)層上下邊界封閉；氣體先流入裂縫再流入井筒；穩(wěn)態(tài)等溫；裂縫在縱向上穿透儲(chǔ)層，高度為氣藏厚度h，在橫向上不完全穿透儲(chǔ)層.這種壓裂分布特征為建立產(chǎn)能預(yù)測(cè)模型帶來一定難度，因此需要利用保角變化把不規(guī)則的裂縫變?yōu)橐?guī)則形狀，見圖1（b）.

1.1 保角變換方法

利用保角變換方法［16－17］把不規(guī)則裂縫變?yōu)榈乳L、與儲(chǔ)層寬度相同的裂縫，以單條縫為例，變化過程見圖2.圖2中垂直裂縫的長度為Lf，高度為hf，這里裂縫完全穿透氣藏，即hf＝h，垂直裂縫位于氣藏任意位置，與垂直裂縫平行的兩側(cè)為恒壓邊界，pe為恒壓邊界壓力，pf為裂縫壓力，氣藏頂?shù)走吔?、與垂直裂縫相垂直的兩側(cè)均為封閉邊界，即

圖1 壓裂水平井多條裂縫示意

p＝pe，與垂直裂縫平行的兩側(cè)；

p＝pf，裂縫中.

當(dāng)裂縫與水平井夾角為α?xí)r，裂縫保角變換見圖3.圖2俯視平面圖的1／2見圖3（a）；把裂縫投射到軸線上，可得圖3（b）；再通過保角變換到圖3（c）的半橢圓平面，該半橢圓平面又可以進(jìn)一步變換到圖3（d）線性系統(tǒng)；因此，垂直裂縫橫向上不完全穿透氣藏的壓裂水平井壓力分布可通過求解線性系統(tǒng)的拉普拉斯方程獲得.

圖2 垂直裂縫示意

圖3 裂縫保角變換示意

式中：φe為地層擬壓力；φf為裂縫處擬壓力；K為滲透率；T為氣藏溫度；Lf為裂縫半長；qfsc為裂縫產(chǎn)量.

1.2 數(shù)學(xué)描述

對(duì)于多條裂縫，相互之間存在壓力干擾，設(shè)裂縫條數(shù)為n條，每條裂縫的產(chǎn)氣量為qfsc1，qfsc2，…，qfscn，裂縫半長分別為Lf1，Lf2，…，Lfn，裂縫寬分別為hf1，hf2，…，hfn，裂縫滲透率分別為Kf1，Kf2，…，Kfn，裂縫與水平段夾角為α1，α2，…，αn，裂縫1和裂縫2的間距為L1，裂縫2和裂縫3的間距為L2，裂縫n－1和裂縫n的間距為Ln－1，見圖1（a）.

根據(jù)勢(shì)的疊加原理，考慮裂縫間的相互干擾，則裂縫1總的干擾擬壓力降為

從裂縫流到水平井筒壓降為

式中：rw為井筒半徑.

裂縫1在條形儲(chǔ)層內(nèi)的總的壓降計(jì)算公式為

式中：φwf1為井筒中第i條裂縫擬壓力；Li為裂縫間距.

水平井產(chǎn)量不高的情況下井筒壓降對(duì)產(chǎn)能的影響可忽略.對(duì)于有n條裂縫壓裂水平井，在定壓生產(chǎn)的情況下，給出井底流動(dòng)擬壓力φwf，則各裂縫處的壓力可知.建立n個(gè)方程為

每條裂縫處壓力差都是n條裂縫產(chǎn)量的函數(shù)，即

其中自變量qfsc有n個(gè)，而方程也有n個(gè)，從而可求每條裂縫的產(chǎn)氣量，井的產(chǎn)量為各條裂縫產(chǎn)量之和.

2 模型應(yīng)用

2.1 實(shí)例分析

DP35－1井是大牛地氣田一口尾管射孔壓裂水平井，其基本參數(shù)見表1，計(jì)算結(jié)果見表2.

表1 DP35－1井基本數(shù)據(jù)

表2 DP35－1井試采計(jì)算結(jié)果

與試采計(jì)算無阻流量（試采穩(wěn)定后一點(diǎn)法計(jì)算無阻）對(duì)比，文中所建模型計(jì)算無阻流量相對(duì)誤差為8.3%，小于10%.

2.2 影響因素

影響因素分析主要考慮裂縫條數(shù)、裂縫長度、裂縫夾角和裂縫位置對(duì)氣井產(chǎn)量的影響，氣藏及水平井參數(shù)見表1.

（1）裂縫條數(shù).裂縫條數(shù)與產(chǎn)量關(guān)系曲線見圖4.由圖4可知，在泄氣面積、裂縫間距和泄氣半徑不變的情況下，水平井產(chǎn)量隨裂縫條數(shù)的增加而增加，但裂縫條數(shù)增到一定數(shù)量后，氣井產(chǎn)量增加逐漸變緩.考慮裂縫間的干擾性，確定大牛地氣藏合理的裂縫間距.

（2）裂縫長度.裂縫長度與產(chǎn)量關(guān)系曲線見圖5.由圖5可知，隨著裂縫半長的增加，產(chǎn)量逐漸增加.通過結(jié)合工程工藝水平和經(jīng)濟(jì)評(píng)價(jià)指標(biāo)，指導(dǎo)大牛地氣藏壓裂水平井裂縫半長的優(yōu)化.

圖4 裂縫條數(shù)與產(chǎn)量關(guān)系曲線

圖5 裂縫長度與產(chǎn)量關(guān)系曲線

（3）裂縫夾角.水平井產(chǎn)量與裂縫夾角的關(guān)系曲線見圖6.由圖6可知，裂縫產(chǎn)氣量隨著裂縫與水平段夾角的增加而增大，這是因?yàn)殡S著夾角的增大，各條裂縫之間的垂直距離變大，也使得各條裂縫之間的相互干擾減小.說明壓裂裂縫與水平井井筒垂直是最好的，因此在水平井方向設(shè)計(jì)時(shí)要考慮儲(chǔ)層地應(yīng)力狀況，從而確定合理的裂縫夾角.

（4）裂縫位置.3種裂縫位置分布見圖7，裂縫位置分別為均勻分布、裂縫在水平段兩端密集和裂縫在水平段中間密集3種情況，分析其產(chǎn)量大小.裂縫位置對(duì)產(chǎn)量的影響見表3.

由表3可知，對(duì)于不同的裂縫布局，不同位置的裂縫應(yīng)該盡量錯(cuò)開排列.為求得最大產(chǎn)量，根部和端部裂縫的間距應(yīng)小于內(nèi)部裂縫的間距.

圖7 裂縫位置分布示意

表3 裂縫位置對(duì)產(chǎn)氣量的影響

3 結(jié)論

（1）所建氣藏壓裂水平井產(chǎn)能預(yù)測(cè)模型不僅考慮了裂縫參數(shù)（裂縫長度、半長、夾角等）的變化，同時(shí)能結(jié)合儲(chǔ)層特征考慮裂縫間的干擾，在理論上更加完善.

（2）結(jié)合現(xiàn)場(chǎng)實(shí)際數(shù)據(jù)，利用所建模型計(jì)算單井產(chǎn)能與試采結(jié)果對(duì)比，相對(duì)誤差為8.3%.

（3）通過分析裂縫參數(shù)對(duì)產(chǎn)能的影響，表明壓裂井水平井產(chǎn)能隨著裂縫條數(shù)、長度和夾角的增加而增加，且裂縫分布方式對(duì)產(chǎn)能影響也較大.

［1］ Giger F M，Resis L H，Jourdan A P.The reservoir engineering aspects of horizontal drilling［C］.SPE 13024，1984.

［2］ Karcher R J，Giger F M，Combe J.Some practical formulas to predicthorizontal well behavior［C］.SPE 15430，1986.

［3］ Soliman M Y，HuntJ L，Raba A M.Fracturing aspects of horizontal well［J］.JPT，1990（8）：66－73.

［4］ Joshi S D.Augmentation of well production using slantand horizontal wells［C］.SPE 15375，1986.

［5］ Raghavan R，Joshi S D.Productivity of multiple drainholes or fractured horizontal wells［C］.SPE Formation Evaluation，1993（3）：11－16.

［6］ Guo B，Schechter D S.A simple and rigorous mathematical model for estimating inflow performance of wells intersecting long fractures［C］.SPE 38104，1997.

［7］ Wei Y，Economides M J.Transverse hydraulic fractures from a horizontal well［C］.SPE 94671，2005.

［8］ Guo Boyun，Yu Xiance.A simple and accurate mathematical model for predicting productivity of multifractured horizontal wells［C］.SPE 114452，2008.

［9］ 郎兆新，張麗華，程林松.壓裂水平井產(chǎn)能研究［J］.石油大學(xué)學(xué)報(bào)：自然科學(xué)版，1994，18（2）：43－46.

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

［11］ 范子菲，方宏長.裂縫性油藏水平井穩(wěn)態(tài)解產(chǎn)能公式研究［J］.石油勘探與開發(fā)，1996，23（3）：52－57.

［12］ 劉振宇，劉洋，賀麗艷，等.人工壓裂水平井研究綜述［J］.大慶石油學(xué)院學(xué)報(bào)，2002，26（4）：96－99.

［13］ 韓樹剛，程林松，寧正福.氣藏壓裂水平井產(chǎn)能預(yù)測(cè)新方法［J］.石油大學(xué)學(xué)報(bào)：自然科學(xué)版，2002，26（4）：37－40.

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

［15］ 徐嚴(yán)波，齊桃，楊鳳波，等.壓裂后水平井產(chǎn)能預(yù)測(cè)新模型［J］.石油學(xué)報(bào)，2006，27（1）：89－92.

［16］ 曲占慶，溫慶志.水平井壓裂技術(shù)［M］.北京：石油工業(yè)出版社，2009.

［17］ 李春蘭，張士誠.魚骨型分支井穩(wěn)態(tài)產(chǎn)能公式［J］.大慶石油學(xué)院學(xué)報(bào)，2010，34（1）：56－59.

A new productivity prediction model of fracturing horizontal well in gas reservoir taking into accountthe mutual interference between the fractures／2012，36（3）：67－71

WANG Shu－ping，SHI Yun－qing，YAN Jin，ZHENG Rong－chen
（Petroleum Exploration ＆Production Research Institute Sinopec，Beijing100083，China）

Because the horizontal well fracturing technique can increase drainage area of the single well，and effectively connectthe reservoir，ithas become an importantmeans to improve gas well flow rate in low－permeability gas reservoirs.During its application，itis difficultto build a productivity predicting model of fracturing horizontal well which takes the changes of reservoir characteristics and fracture parameters into account.Aiming atthis problem，the paper established a new fracturing horizontal well productivity prediction model，ittakes into accountnotonly the impactof the fracture half－length，fracture spacing and fracture angle changes on horizontal well productivity，butalso the mutual interference between the fractures.Through on－site verification，the results of the model are accurate.

fracturing horizontal well；productivity prediction；calculation model；mutual interference of fractures；low－permeability gas reservoir

book=3,ebook=69

TE313.8

A

1000－1891（2012）03－0067－05

2011－11－21；編輯：關(guān)開澄

中石化科技攻關(guān)項(xiàng)目（G5800－10－ZS－YK004－2）

王樹平（1978－），男，博士，高級(jí)工程師，主要從事油氣藏工程方面的研究.