段慕白　李　皋　孟英峰　田　旭

“油氣藏地質(zhì)及開發(fā)工程”國家重點實驗室·西南石油大學

?

氣體鉆井地層動態(tài)出水量預(yù)測計算模型

段慕白李皋孟英峰田旭

“油氣藏地質(zhì)及開發(fā)工程”國家重點實驗室·西南石油大學

段慕白等.氣體鉆井地層動態(tài)出水量預(yù)測計算模型.天然氣工業(yè)，2016,36(6):66-71.

摘 要氣體鉆井地層出水容易導(dǎo)致井眼清潔困難、井壁失穩(wěn)等井下復(fù)雜問題。為了準確預(yù)測氣體鉆井打開水層過程中的地層動態(tài)出水量，綜合考慮了非穩(wěn)定滲流對地層壓力分布的影響和水層被打開的程度，建立了氣體鉆井地層動態(tài)出水量定量預(yù)測計算模型。采用該模型的計算結(jié)果表明：①鉆開水層的厚度小于井徑、鉆頭進入水層但未完全鉆穿水層、鉆頭部分鉆穿水層而未完全脫離水層時，出水量主要受水層打開程度的控制，水層打開程度越高出水量越大，在鉆頭完全鉆穿水層時出水量達到峰值，該時間段較短，出水量大，應(yīng)當嚴密監(jiān)測，以便在鉆井過程中及早發(fā)現(xiàn)水層；②完全鉆穿水層后，出水量主要受到地層非穩(wěn)定滲流的影響，隨著地層壓力的衰減，地層出水量隨之衰減，最終趨于穩(wěn)定，該穩(wěn)定值為氣體鉆井適用性地層篩選和轉(zhuǎn)換鉆井方式時機的參考值?，F(xiàn)場隨鉆監(jiān)測驗證結(jié)果表明,井口返出氣體的濕度變化規(guī)律與計算結(jié)果相吻合。結(jié)論認為，該模型能更真實地反映地層出水量的實際變化情況，為氣體鉆井適用性地層篩選和轉(zhuǎn)換鉆井方式時機提供了理論依據(jù)。

關(guān)鍵詞氣體鉆井非穩(wěn)定滲流地層壓力分布地層出水量出水量定量預(yù)測模型水層打開程度鉆井方式地層篩選

氣體鉆井因其鉆速快、保護儲層、防止惡性井漏等優(yōu)勢在全球范圍內(nèi)都得到了廣泛使用，但地層出水容易導(dǎo)致井眼凈化困難、井壁失穩(wěn)等問題，對地層出水很敏感，出水量過大則只有轉(zhuǎn)換鉆井方式。因此鉆前出水量的準確預(yù)測成為篩選氣體鉆井適應(yīng)性地層和選擇轉(zhuǎn)換鉆井方式時機的關(guān)鍵[1-3]。

氣體鉆井過程中可以根據(jù)鉆井參數(shù)和地面返出巖屑形態(tài)的變化來判斷地層是否出水，粗略估計出水量[4]?；诜€(wěn)定滲流力學，李祖光等結(jié)合氣體鉆井工況，根據(jù)水層的打開程度推導(dǎo)出不同滲流模型的出水量計算公式[5]。鄒靈戰(zhàn)等通過判斷流體類型、地層物性解釋、地層壓力解釋和出水量計算等步驟來預(yù)測地層出水量[6]。劉彪等通過鄰井地層對比及測井與錄井資料評價，從流體物性、地層物性、地層出水量3個方面對侵入井筒的地層水量進行了分析[7]。Jiang Yufang等[8]、楊決算[9]和李皋等[10]建立了圓形封閉性天然水域系統(tǒng)下非穩(wěn)定滲流地層出水量計算模型，確定非穩(wěn)定滲流出水量計算公式。Nguyen等耦合了鉆井參數(shù)、井筒流動參數(shù)以及流體參數(shù)等建立了欠平衡鉆井產(chǎn)能預(yù)測模型，并開發(fā)了UBDDRILL軟件評估UBD的產(chǎn)能[11]，此模型也可以用來預(yù)測地層出水量。趙向陽等總結(jié)分析了現(xiàn)有的氣體鉆井地層出水預(yù)測方法，提出隨鉆監(jiān)測方法，通過觀測壓力和注入?yún)?shù)的變化求取出水量[12]。范翔宇等基于隨鉆測井技術(shù)，對地層出水層位的預(yù)測方法進行研究，建立氣體鉆井出水量定量預(yù)測模型[13]。

實際上，地層水隨著水層被打開而流入井筒，滲流過程為非穩(wěn)定滲流，地層水滲出速度隨壓降漏斗擴大而減小，出水斷面隨著水層被逐漸打開而增大，直至完全鉆穿水層，出水量受兩者同時控制。已有的研究沒有將非穩(wěn)定滲流與水層的打開程度相結(jié)合。本文在非穩(wěn)定滲流力學基礎(chǔ)上，結(jié)合機械鉆速和水層被揭開后地層壓力的分布規(guī)律，根據(jù)水層被打開的不同程度，建立從開始打開水層到完全鉆開水層及之后的地層出水量動態(tài)計算模型，從而獲得出水量隨時間變化的動態(tài)關(guān)系。

1　水層打開程度與滲流模型

水層被打開的程度與機械鉆速相關(guān)，水層被打開的程度越大，過流斷面面積就越大，出水量為水的流速與過流斷面之積。水層的地層壓力變化遵循非穩(wěn)定滲流規(guī)律，將地層壓力變化與水層打開程度耦合，機械鉆速和時間決定水層的打開程度，建立非穩(wěn)定滲流與地層打開程度耦合的出水量動態(tài)預(yù)測模型。

在此將水層的打開程度分為：鉆開水層的厚度小于井徑、鉆頭進入水層但未完全鉆穿水層、鉆頭部分鉆穿水層而未完全脫離水層、完全鉆穿水層[14-15]。

1.1鉆開水層的厚度小于井徑

鉆頭從接觸水層到完全進入水層，即h≤rw，地層水通過井底球面涌入井筒，可將其滲流情況視為球面向心流，如圖1所示，出水面積為：

式中As表示出水斷面面積，m2；rw表示井徑，m；h表示鉆頭鉆入水層厚度，m。

圖1　鉆開水層厚度小于井徑時地層出水示意圖

1.2鉆頭尚未鉆穿水層

鉆頭完全進入水層，但尚未鉆穿水層，滿足hwater≥h＞rw，地層水通過鉆頭之上的筒狀地層和鉆頭底部球面地層進入井筒，可將上部視為平面徑向流，下部視為球面向心流，如圖2所示，出水面積為上部井筒面積與下部球面面積之和：

圖2　鉆頭完全進入水層但尚未鉆穿水層時地層出水示意圖

1.3鉆頭部分鉆穿水層

鉆頭頂部鉆穿出水層，而胎體尚未完全脫離水層，滿足hwater+rw＞h＞hwater，仍將上部視為平面徑向流，下部視為球面向心流，如圖3所示，出水面

圖3　鉆頭未完全脫離水層時地層出水流動示意圖

積為上部井筒面積與下部部分球面面積之和：

式中hwater表示水層厚度，m。h =Vpt，Vp表示機械鉆速，m/h，t表示時間，h。

1.4完全鉆穿水層

鉆頭完全脫離水層，h＞hwater+rw，地層水僅通過上部以平面徑向流的方式進入井筒，如圖4所示，此時， 出水面積為：

圖4　鉆頭完全鉆穿水層的地層水流動示意圖

2　非穩(wěn)定滲流地層壓力動態(tài)分布

球面向心非穩(wěn)定流的偏微分方程[16]：

式中p表示地層孔隙壓力，MPa；pi表示原始地層壓力，MPa；pw表示井底壓力，MPa；η表示導(dǎo)壓系數(shù)，cm2/s；φ表示孔隙度， K表示滲透率，μ表示黏度（mPa·s），Ct表示綜合壓縮系數(shù)（MPa-3）。

平面徑向非穩(wěn)定流的偏微分方程：

式中t表示時間，min；r表示徑向距離，m；

結(jié)合邊界條件與初始條件，采用差分方法對上式離散可得，球面向心流近井壁地帶孔隙壓力分布方程：

平面徑向流近井壁地帶的孔隙壓力分布方程：

3　出水量計算

在井壁和井底取薄層巖石微元，如圖5、6所示，入井地層水經(jīng)過該微元流入井筒，求出微元內(nèi)的壓降梯度，從而獲得流經(jīng)微元體的地層水流速，再根據(jù)過流面積，即可求得入井流量。

圖5　井周平面徑向流微元示意圖

圖6　井底球面向心流微元示意圖

井底巖石不斷被鉆頭揭開，井底位置下降，與此同時，由于非穩(wěn)定滲流，地層壓降漏斗在水層被揭開后不斷擴大，但壓降傳播速度比機械鉆速快得多，因此將井底位置降低對壓力分布的影響忽略。根據(jù)式（7）、（8）采用MathCAD軟件編程，求得的地層壓力隨時間和空間的分布情況，然后求取微元體內(nèi)的壓降梯度，從而根據(jù)計算出水量。

地層水滲流速度偏微分方程[16]：

式中v表示地層水滲流速度，m/s。

流量計算公式：

式中Q表示流量，m3/h。

根據(jù)地層壓力隨時間和空間的分布，對壓降進行差分離散得：

從而獲得地層水流量公式：

地層水黏度受溫度的影響明顯，據(jù)本文參考文獻[17]，地層水黏度與溫度的關(guān)系如下：

式中T表示溫度，K；s表示含鹽量，mg/L。

4　實例分析

某油田ND102井鉆至4 616.99m時因嚴重井漏決定采用氣體鉆井，氣體鉆井鉆至4 792m，立壓上升，起鉆遇阻劃眼，氮氣循環(huán)，排砂管線出口呈霧狀返出，判斷地層出水。由鄰井實鉆資料及測井分析該井水層基礎(chǔ)數(shù)據(jù)：水層深度為4 786m，井筒半徑為0.108m，井底流壓為2.1MPa，滲透率為5.0 mD，孔隙度為10%，孔隙壓力為50.2MPa，機械鉆速為8.4m/h，地層水黏度為0.45 mPa·s，水層厚度1.13m。

氣體鉆井揭開水層后， 地層水通過井底和井壁滲入井筒，地層壓力和地層水流速隨滲流的進行不斷變化，由本文所建立的模型，根據(jù)上述基礎(chǔ)數(shù)據(jù)，采用MathCAD軟件編程計算得出1 h內(nèi)平面徑向流和10min內(nèi)球面向心流在近井壁地帶的地層壓力與地層水流速變化，計算結(jié)果如圖7～10所示。

氣體鉆井速度快，可在短時間內(nèi)鉆穿水層，鉆穿水層后井底不再出水，僅由井壁平面徑向流供水，因此，只表示出短時間內(nèi)井底球面向心流產(chǎn)生的地層壓力變化。10min內(nèi)井底球面向心流地層壓力隨時間和徑向距離變化如圖7所示。

圖7　球面向心流近井壁壓力分布圖

圖8　平面徑向流近井壁壓力分布圖

圖9　球面向心流井底微元內(nèi)地層水流速變化曲線圖

圖10　平面徑向流井壁微元內(nèi)地層水流速變化曲線圖

由圖9可知，鉆頭揭開水層后，井底地層水滲出，井底附近地層壓力急劇下降，越靠近井底壓降越大，初始時刻的壓降速度最大，而后減緩，由地層水流速計算公式可知，在揭開井底初始時刻的地層水流速最大，之后逐漸趨于平緩；鉆頭鉆入水層后，形成井眼，地層水隨之滲出，近井壁地帶地層壓力不斷降低，形成壓降漏斗，并隨時間推移不斷擴大，如圖8所示，在形成井眼后，井壁附近地層壓力急劇降低，越靠近井壁下降幅度越大，初始時刻的壓降速度最大，而后減緩，由地層水流速計算公式可知，井壁微元內(nèi)地層水流速度在初始時刻最大而后減緩，如圖10所示。

鉆頭在短時間內(nèi)全部鉆入水層，地層水入井流量隨地層被打開的程度而變化，出水量隨時間的變化關(guān)系如圖11所示。

圖11　出水量隨時間變化關(guān)系圖

由圖11可知，在鉆頭完全鉆穿水層前，隨著鉆頭不斷進入水層，出水量不斷增大，鉆頭鉆穿水層后，出水量隨時間推移而緩慢減小，趨于穩(wěn)定。在鉆頭完全進入地層前，滲流僅為球面向心流，如圖11中的1區(qū)域所示，出水量隨鉆頭吃入深度而增加；鉆頭全部進入水層后，地層水滲流為上部平面徑向流與井底球面向心流的疊加，出水量為兩者之和，由于氣體鉆井速度快，會在短時間內(nèi)鉆穿水層，在這段時間內(nèi)地層壓力也會隨鉆進而衰減，但主要影響出水量的還是水層被鉆開的深度，隨著鉆頭不斷鉆開地層，出水量不斷增加，如圖11中的2區(qū)域所示，而后鉆頭頂部穿出水層，但未完全脫離水層，此時出水量由上部的平面徑向流與底部部分球面向心流疊加，出水量繼續(xù)增加，鉆頭快速穿出水層，此段時間較短，如圖11中的3區(qū)域所示。鉆頭完全脫離水層后，井底不再供水，僅由井壁供水，由于非穩(wěn)定滲流，地層壓力隨時間的推移而降低，在完全鉆穿水層之前地層壓力已有所降低，在完全鉆穿水層之后壓降速度繼續(xù)減小，出水量亦緩慢減小，如圖11中的4區(qū)域所示，鉆頭在接觸水層后8.7min之內(nèi)便鉆穿水層，出水流動形態(tài)轉(zhuǎn)變?yōu)閱我环欠€(wěn)定平面徑向流，出水量逐漸減小并趨于穩(wěn)定值。針對某油田ND102井，本文模型計算得到的最大值為5.2m3/h，穩(wěn)定值為3.1m3/h，與現(xiàn)場監(jiān)測得到的出水量在3～5m3/h間吻合。本模型能真實地反映地層出水量的實際變化情況，更準確地預(yù)測氣體鉆井過程中的地層出水量。

5　結(jié)論

基于非穩(wěn)定滲流理論，綜合考慮地層壓力分布的影響和水層被打開的程度，建立了氣體鉆井地層動態(tài)出水量預(yù)測計算模型。

基于該模型計算得到，鉆開水層的厚度小于井徑、鉆頭進入水層但未完全鉆穿水層、鉆頭部分鉆穿水層而未完全脫離水層時，出水量主要受水層打開程度的控制，水層打開程度越高出水量越大，在鉆頭完全鉆穿水層時出水量達到峰值，該時間段較短，出水量大，應(yīng)當嚴密監(jiān)測，以便在鉆井過程中及早發(fā)現(xiàn)水層；完全鉆穿水層后，出水量主要受到地層非穩(wěn)定滲流的影響，隨著地層壓力的衰減，地層出水量隨之衰減，最終趨于穩(wěn)定，該穩(wěn)定值為氣體鉆井適用性地層篩選和轉(zhuǎn)換鉆井方式時機的參考值。

現(xiàn)場隨鉆監(jiān)測結(jié)果表明：井口返出氣體的濕度變化規(guī)律與計算結(jié)果相吻合，該模型能更真實地反映地層出水量的實際變化情況，為氣體鉆井適用性地層篩選和轉(zhuǎn)換鉆井方式時機提供理論依據(jù)。

參考文獻

[1] Lyons WC,Guo Boyun,Graham RL,Hawley GD.Chpater 8–Air,gas,and unstable foam drilling[M]//Air & Gas Drilling Manual,3rded.New York: McGraw-Hill,2009:183-212.

[2] Guo Boyun,Gao Deli.New development of theories in gas drilling[J].Petroleum Science,2013,10(4): 507-514.

[3] 吳志均,唐紅君.淺談氣體鉆井需要關(guān)注的問題[J].鉆采工藝,2008,31(3): 28-30.Wu Zhijun,Tang Hongjun.Discussion on key points of gas drilling[J].Drilling & Production Technology,2008,31(3): 28-30.

[4] 韓朝輝,張建斌,范興沃,張漢林,賈軍喜.川東北地區(qū)空氣鉆井中幾種復(fù)雜情況診斷與處理[J].鉆采工藝,2008,31(6): 24-25.Han Chaohui,Zhang Jianbin,Fan Xingwo,Zhang Hanlin,Jia Junxi.Diagnostic and disposal of several complicated problemsfor air drilling in Northeast Sichuan area[J].Drilling & Production Technology,2008,31(6): 24-25.

[5] 李祖光,蔣宏偉,翟應(yīng)虎,張旭.氣體鉆井地層出水量理論模型[J].遼寧工程技術(shù)大學學報: 自然科學版,2009,28(增刊1): 8-10.Li Zuguang,Jiang Hongwei,Zhai Yinghu,Zhang Xu.Water output model in gas drilling[J].Journal of Liaoning Technical University: Natural Science Edition,2009,28(S1): 8-10.

[6] 鄒靈戰(zhàn),鄧金根,汪海閣.氣體鉆井地層出水定量預(yù)測技術(shù)[J].石油鉆探技術(shù),2009,37(3): 30-33.Zou Lingzhan,Deng Jingen,Wang Haige.Quantitative prediction of formation water production during gas drilling[J].Petroleum Drilling Techniques,2009,37(3): 30-33.

[7] 劉彪,劉剛,劉成,孫磊,梁昌健.空氣鉆井預(yù)測地層出水量的探討[J].斷塊油氣田,2009,16(3): 93-95.Liu Biao,Liu Gang,Liu Cheng,Sun Lei,Liang Changjian.Further discussion on prediction of formation water discharge by air drilling[J].Fault-Block Oil & Gas Field,2009,16(3): 93-95.

[8] Jiang Yufang,Liu He,Liu Zhihuan,Liu Yonggui,Li Zengle.Study and application of formation water production and borehole wall stability prediction technologies in Daqing gas/foam drilling[C]// IPTC 2009: International Petroleum Technology Conference.Houten: EAGE,2009.

[9] 楊決算.大慶油田氣體鉆井地層出水與井壁穩(wěn)定技術(shù)研究[D].大慶: 東北石油大學,2010.Yang Juesuan.A study on predicting formation water production and wellbore stability for gas drilling in Daqing oilfield[D].Daqing: Northeast Petroleum University,2010.

[10] 李皋,李誠,孟英峰,李永杰,劉金龍,陳一健,等.氣體鉆井隨鉆安全風險識別與監(jiān)控[J].天然氣工業(yè),2015,35(7): 66-72.Li Gao,Li Cheng,Meng Yingfeng,Li Yongjie,Liu Jinlong,Chen Yijian， et al.While-drilling safety risk identification and monitoring in air drilling[J].Natural Gas Industry,2015,35(7): 66-72.

[11] Nguyen C,Somerville JML,Smart BGD.Predicting the production capacity during underbalanced-drilling operations in Vietnam[C]//IADC/SPE Managed Pressure Drilling and Underbalanced Operations Conference & Exhibition,12-13 February 2009,San Antonio,Texas,USA.DOI:http://dx.doi.org/10.2118/122266-MS.

[12] 趙向陽,孟英峰,楊順輝,李皋,李永杰,周玉良.氣體鉆井地層出水隨鉆監(jiān)測與攜水規(guī)律研究[J].斷塊油氣田,2012,19(1): 120-123.Zhao Xiangyang,Meng Yingfeng,Yang Shunhui,Li Gao,Li Yongjie,Zhou Yuliang.Monitoring while drilling and water carrying law for formation water production in gas drilling[J].Fault-Block Oil & Gas Field,2012,19(1): 120-123.

[13] 范翔宇,蘇靜,王俊瑞,張萍,何傳亮.川東北河壩地區(qū)中淺部地層氣體鉆井出水隨鉆預(yù)測[J].中國石油大學學報: 自然科學版,2013,37(1): 85-90.Fan Xiangyu,Su Jing,Wang Junrui,Zhang Ping,He Chuanliang.Water exit prediction while drilling from mid-shallow formation during gas drilling in Heba area,Northeast Sichuan[J].Journal of China University of Petroleum: Edition of Natural Science,2013,37(1): 85-90.

[14] 翟云芳.滲流力學[M].2版.北京: 石油工業(yè)出版社,2003: 13-19.Zhai Yunfang.Permeation fluid mechanics[M].2nded.Beijing: Petroleum Industry Press,2003: 13-19.

[15] 鄧英爾.高等滲流理論與方法[M].北京: 科學出版社,2004: 35-44.Deng Ying'er.Higher seepage theory and method[M].Beijing: Science Press,2004: 35-44.

[16] 葛家理,寧正福.現(xiàn)代油藏滲流力學原理[M].北京: 石油工業(yè)出版社,2001.Ge Jiali,Ning Zhengfu.Modern reservoir seepage mechanics principle[M].Beijing: Petroleum Industry Press,2001.

[17] Kincaid JF,Eyring H,Stearn AE.The theory of absolute reaction rates and its application to viscosity and diffusion in the liquid state[J].Chemical Reviews,1941,28(2): 301-365.

(修改回稿日期2016-03-15編輯凌忠)

A prediction and calculation model for dynamic formation water yield in gas drilling

Duan Mubai,Li Gao,Meng Yingfeng,Tian Xu
(State Key Laboratory of Oil & Gas Reservoir Geology and Exploration， Southwest Petroleum University， Chengdu，Sichuan 610500,China)

NATUR.GAS IND.VOLUME 36,ISSUE 6,pp.66-71,6/25/2016.(ISSN 1000-0976; In Chinese)

Abstract:Water breakthrough during gas drilling tends to result in complex downhole problems， such as borehole instability and the difficulty in hole cleaning.For accurate prediction of the dynamic formation water yield while the water layer is being drilled by gas drilling,a dynamic formation water yield prediction model was built by analyzing comprehensively the effect of unsteady seepage on the formation pressure distribution and the drilling degree of water layers.It is shown from the calculation results of this model that the water yield is mainly controlled by the drilling degree of a water layer when the drilling thickness of the water layer is less than that of the well diameter,or when the bit drills into the water layer but doesn't run through it completely,or when the bit partially penetrates the water layer but doesn't run out of it completely.The higher the drilling degree of water layer is,the more the water yield is.The peak water yield is reached when the water layer is penetrated completely by the bit.And this period lasts short with a large water yield.Therefore,it is necessary to carry out intensive monitoring so as to distinguish water layers in the process of drilling as soon as possible.After the water layer is completely penetrated,the water yield is mainly controlled by the unstable seepage.Formation water yield decreases with the attenuation of formation pressure， and finally tends to be a stable value， which is a reference for applicable formation screening and drilling pattern change timing in gas drilling.It is indicated from the field monitoring and verifying results that the humidity change law of the wellhead returning gas is in accordance with the calculation results.It is concluded that this model can reflect the actual changes of the formation water yield more accurately， and provide a theoretical basis for applicable formation screening and drilling pattern change timing in gas drilling.

Keywords:Gas drilling; Unsteady seepage; Distribution of formation pressure; Formation water yield; Quantitative water yield prediction; Model; Drilling degree of water layer; Drilling pattern; Formation screening

DOI:10.3787/j.issn.1000-0976.2016.06.010

基金項目：四川省杰出青年學術(shù)技術(shù)帶頭人培育計劃“氣體鉆井井下安全風險及控制方法研究”(編號:2014JQ0045)。

作者簡介：段慕白，1984年生，博士；從事氣體鉆井、欠平衡鉆井、裂縫精細描述、裂縫閉合機理、壓力波動裂縫變形機理等方面的研究工作。地址:(610500)四川省成都市新都區(qū)新都大道8號西南石油大學“油氣藏地質(zhì)及開發(fā)工程”國家重點實驗室。電話:(028) 83034786。ORCID:0000-0002-3012-1773。E-mail:duan_mubai@foxmail.com