龔迪光,曲占慶,郭天魁,鞏法成,田旭新,黃子桐
(中國石油大學(華東)石油工程學院,山東青島266580)
?
徑向井水力壓裂摩阻影響因素與計算公式
龔迪光,曲占慶,郭天魁,鞏法成,田旭新,黃子桐
(中國石油大學(華東)石油工程學院,山東青島266580)
龔迪光等.徑向井水力壓裂摩阻影響因素與計算公式[J].鉆井液與完井液,2016,33(3):102-106.
摘要為研究井眼內徑及壓裂施工參數(shù)對徑向井壓裂摩阻影響規(guī)律和確定井眼內壓裂液摩阻大小,使用中國石油大學(華東)研制的壓裂液摩阻測試系統(tǒng)(主要由摩阻測試控制中心、壓裂液調配釜、可調速螺桿泵、變徑管道、高靈敏度壓力測量儀、電子流量計組成),模擬徑向井中壓裂液的流動狀態(tài),并對壓裂液摩阻進行準確測量,分析了影響壓裂液摩阻的主要因素。實驗結果表明,對徑向井摩阻影響由大到小的因素依次為:井眼內徑、排量、黏度、支撐劑粒徑和砂比,且影響規(guī)律各有不同。采用降阻比原理,通過對322組實驗數(shù)據(jù)進行回歸擬合,建立了考慮井眼內徑、排量、壓裂液黏度、支撐劑粒徑及砂比的徑向井瓜膠壓裂液摩阻損失計算關系式。利用相關系數(shù)檢驗法計算標準估計誤差為0.140,擬合回歸方程有效。實驗發(fā)現(xiàn),瓜膠壓裂液黏度對摩阻有雙重影響:一方面會增加流體內部以及流體與管壁間的剪切應力,導致摩阻損失增大;另一方面,隨著黏度增加,聚合物溶液產生轉捩延遲效應,同時對支撐劑控制能力增強,促使壓裂液摩阻損失減小。
關鍵詞摩阻測試;摩阻影響因素;徑向井;降阻比;水力壓裂;壓裂液;正交實驗
徑向井壓裂技術不但可以改善油藏滲流條件,而且在提高產能的基礎上可大大節(jié)省壓裂成本。但是,由于徑向井井眼較小,壓裂液在徑向井眼中的流動摩阻要比常規(guī)井大很多,摩阻的大小直接影響壓裂設計以及施工的成敗。研究徑向井井眼中壓裂液摩阻損失可以為徑向井壓裂設計和施工提供重要依據(jù)。目前,攜砂壓裂液摩阻的求取大多采用等效修正的方法,將支撐劑顆粒和液體等效看作一種流體,采用理論分析計算摩阻取值,而實驗研究較少。針對這一難題,中國石油大學(華東)壓裂酸化研究中心結合徑向井井眼特點,研發(fā)了一套測量徑向井攜砂壓裂液摩阻的設備,可模擬徑向井中壓裂液的流動狀態(tài),并對壓裂液摩阻進行準確測量,分析其內在影響因素[1-3]。
1.1實驗裝置
壓裂液摩阻測試系統(tǒng)主要由摩阻測試控制中心、壓裂液調配釜、可調速螺桿泵、變徑管道、高靈敏度壓力測量儀、電子流量計組成。實驗裝置設計原則以相似準則為依據(jù),徑向井管徑管長與現(xiàn)場參數(shù)按照1∶1的比例設計。利用可調速螺桿泵模擬壓裂液施工不同排量,每隔一定距離在不同半徑管柱上安裝高靈敏度壓力測量儀,對管柱沿程壓力進行實時測量,近而可得到沿程摩阻損失量。
1.2壓裂液和支撐劑
實驗采用0.3%瓜膠粉+30 ℃清水+不同量交聯(lián)劑(0.4%、0.5%和0.6%)配制成不同黏度(50、100和150 mPa·s)的瓜膠壓裂液,支撐劑選擇卡博陶粒支撐劑。
1.3實驗方法
根據(jù)前人理論研究[4-7]認為,影響壓裂液摩阻的主要因素有排量、壓裂液黏度、管徑、支撐劑粒徑和砂比。針對5種因素設計相應水平進行正交實驗。受實驗條件限制,各因素的水平數(shù)不同,設計L18(43×21)四因素三水平和一因素兩水平混合正交實驗??紤]不同因素實驗水平不一樣,對實驗結果進行極差分析時,需對極差R進行修正,修正后極差R′= dR。其中,r為不同因素每個水平重復試驗次數(shù),d為折算系數(shù),不同水平數(shù)折算系數(shù)不同,見表1[8]。修正極差值越大,表明對應的因子對摩阻損失影響程度越大,通過修正后的極差R′衡量因素的主次順序,更能真實反映各因素對摩阻的影響程度。
表1 不同水平數(shù)折算系數(shù)表
2.1確定摩阻影響因素的主次順序
設計四因素三水平和一因素兩水平的混合正交實驗,實驗數(shù)據(jù)及分析結果見表2和表3。通過比較修正極差值的大小來判斷各個因素對摩阻損失影響程度,確定各因素對徑向井摩阻的影響程度依次為:管徑(井眼內徑)>排量>黏度>支撐劑粒徑>砂比。
圖5為軌道連續(xù)高低變化仿真。由圖5可知,在4根軌枕范圍,即4.5 m內,軌道連續(xù)高低變化對軌道直線度具有顯著影響。在前兩個監(jiān)測探頭剛進入軌道連續(xù)高低變化區(qū)域和后兩個監(jiān)測探頭臨離開軌道連續(xù)高低變化區(qū)域時,直線度達到峰值。
表2 壓裂液摩阻測試混合正交實驗數(shù)據(jù)
表3 實驗結果分析表
2.2單因素分析
以各影響因素水平值為橫坐標,正交實驗指標平均值(ki)對應的摩阻損失值為縱坐標,繪制單因素與摩阻損失關系曲線圖,見圖1~圖5。
由各因素與摩阻損失關系圖可直觀看出因素水平對瓜膠壓裂液摩阻的影響規(guī)律。
1)管徑從20 mm增加到30 mm,壓裂液摩阻損失大幅度降低,降低幅度為78.5%,驗證了井眼內徑是影響壓裂液摩阻的最主要因素。較大內徑可以減小壓裂液在井眼中的紊流強度以及壓裂液與管壁的相互作用,減小摩阻損失。
2)排量是影響摩阻的重要因素。隨著排量的增加,摩阻顯著增大,當實驗排量從3.0 m3/h增加到8.0 m3/h,摩阻增幅為60.1%。增大排量增加了支撐劑相互碰撞幾率及壓裂液與管壁剪切應力;同時,增加排量會改變壓裂液流動狀態(tài),增加紊流強度,導致壓裂液損失較多的能量,增大壓裂液的摩阻損失。
3)黏度是摩阻不可忽視的影響因素。實驗結果與通常認識有所不同。當壓裂液表觀黏度從50 mPa·s增加到150 mPa·s,壓裂液摩阻呈現(xiàn)先增大再減小的規(guī)律。壓裂液黏度對摩阻有雙重影響:①會增加流體內部以及流體與管壁間的剪切應力,導致摩阻損失增大;②隨著黏度繼續(xù)增加,聚合物溶液產生轉捩延遲效應,對紊流狀態(tài)起到抑制作用,流動形態(tài)有所改變,同時壓裂液對支撐劑控制能力增強,減小了支撐劑運動空間、顆粒間相互碰撞以及與管壁摩擦幾率,促使壓裂液摩阻損失減小。
4)支撐劑顆粒大小與砂比對摩阻影響較小。支撐劑粒徑從20/40目減小到40/70目,壓裂液摩阻損失減小6.7%。
5)砂比從5%增加到15%,摩阻增加1.7%。隨著支撐劑粒徑與砂比的增加,壓裂液摩阻都有一定的增大,但相比其他幾種因素,對摩阻影響較小。
圖1 摩阻損失與管徑擬合曲線
圖2 摩阻損失與排量擬合曲線
圖3 摩阻損失與黏度擬合曲線
圖4 摩阻損失與支撐劑粒徑擬合曲線
圖5 摩阻損失與砂比擬合曲線
2.3摩阻計算公式擬合
瓜膠壓裂液在管柱內多為復雜的紊流流動,很難用具體的公式表示描述[9-11]。因此,利用數(shù)學推導給出較準確的壓裂液摩阻計算公式幾乎不太可能。由于回歸方程表達式未知,經(jīng)過大量嘗試與探索,發(fā)現(xiàn)以降阻比原理為指導,擬合壓裂液摩阻計算公式誤差較小。擬通過實驗方法對影響徑向井壓裂液摩阻損失的5種因素進行多元非線性擬合回歸,給出摩阻計算公式,并利用相關系數(shù)檢驗法對擬合結果進行檢驗。
利用實驗數(shù)據(jù)對Lord和MC Cowen等人提出的降阻比公式和摩阻計算公式進行修正,可得到適合徑向井的摩阻計算公式[12-13]。
降阻比公式為:
式中,ΔP0為清水摩阻損失,Pa。清水摩阻損失采用下式[14]計算。
式中,f為范寧摩阻系數(shù)[15];u為管內清水流速,m/s;μ0為清水黏度,Pa·s;(30 ℃清水黏度為0.800 7 mPa·s);Q0為管內清水排量,m3/s;D為管徑,m;ρ為清水密度,kg/m3;L為實驗管道長度,m。
擬合公式為
式中,ΔPp為壓裂液摩阻損失,Pa;d為管徑,mm;Q為排量,m3/min;G為壓裂液稠化劑濃度,kg/m3(實驗濃度為3 kg/m3);C為支撐劑濃度,kg/m3。
式中,c為壓裂液砂比,無因次;ρS為支撐劑體密度,kg/m3。
測量10目到100目(0.154~2.000 mm)卡博陶粒支撐劑體密度,得到支撐劑粒徑與體密度的回歸關系式:
式中,l為支撐劑粒徑,mm,則
將式(5)和式(8)代入式(4),得
根據(jù)322組實驗數(shù)據(jù)對(9)式進行系數(shù)回歸與化簡,得到摩阻損失與管徑、排量、壓裂液表觀黏度、砂比及支撐劑粒徑的關系式:
式中,μ為壓裂液表觀黏度,Pa·s。
利用相關系數(shù)法對擬合方程進行檢驗,計算相關系數(shù)平方根為0.927,相關系數(shù)為0.860,標準估計誤差為0.140,認為擬合的回歸方程有效,可作為徑向井摩阻的計算公式。這也為徑向井水力壓裂施工參數(shù)的設計和摩阻的計算提供了科學依據(jù)和理論方法。
1.考慮影響壓裂液摩阻損失的5種因素(管徑、排量、壓裂液黏度、支撐劑粒徑、砂比),針對徑向井井眼特點和瓜膠壓裂液流動特征,設計了四因素三水平和一因素兩水平的混合正交實驗。實驗結果顯示,影響徑向井摩阻的因素由大到小依次為:井眼內徑(管徑)>排量>黏度>支撐劑粒徑>砂比。
2.采用降阻比原理,對實驗數(shù)據(jù)進行回歸擬合,建立了考慮井眼內徑、排量、壓裂液黏度、支撐劑粒徑及砂比的徑向井瓜膠壓裂液摩阻損失計算關系式。利用相關系數(shù)檢驗法計算標準估計誤差為0.140,擬合回歸方程有效。
3.實驗發(fā)現(xiàn)瓜膠壓裂液黏度對摩阻有雙重影響:一方面會增加流體內部以及流體與管壁間的剪切應力,導致摩阻損失增大;另一方面,隨著黏度增加,聚合物溶液產生轉捩延遲效應,同時對支撐劑控制能力增強,促使壓裂液摩阻損失減小。
參考文獻
[1]MEDIAN B,SHAH S N.A new approach for predicting frictional pressure losses of non-newtonian fluids in coiled tubing[C].SPE 60319.
[2]張廣明,劉合,張勁.油井水力壓裂流-固耦合非線性有限元數(shù)值模擬[J].石油學報,2009,30(1):113-116. ZHANG Guangming,LIU He,ZHANG Jin.Simulation of hydraulic fracturing of oil well based on fluid-solid coupling equation and non-linear finite element [J].Acta Petrolei Sinica,2009,30(1):113-116.
[3]ADACHI J ,SIEBRITS E,PEIRCE A.Computer simulation of hydraulic fractures [J].International Journal of Rock Mechanics and Mining Sciences,2007,44:739-757.
[4]李勇明,趙金洲,毛虎,等.水平井壓裂攜砂液摩阻定量計算模型研究[J].西南石油大學學報(自然科學版),2011,33(6):81-84. LI Yongming,ZHAO Jinzhou,MAO Hu,et al. Quantitive computational model of the carrier fluid drag of hori-zontal-well fracturing[J].Journal of Southwest Petroleum University(Science & Technology Edition),2011,33(6):81-84.
[5]陳俊,張勁軍.用表觀黏度計算冪律流體摩阻的方法及其應用[J].西安石油學院學報(自然科學版),2003,18(6):58-62. CHEN Jun,ZHANG Jinjun.A method of calculating the pressure-drop of power-law fluid by apparent viscosity and its application[J].Journal of Xi′an Petroleum Institute (Natural Science Edition),2003,18(6):58-62.
[6]RESENDE P R,ESCUDIER M P,PRESTI F,et al.Numerical predictions and measurements of reynolds normal stresses in turbulent pipe flow of polymers[J]. International Jounmal of Heat and Fluid Flow,2006,27:204-219.
[7]MCELFRESH P M,WOOD W R,WILLIAMS C F,et al. A study of the friction pressure and proppant transport behavior of surfactant-based gels[J].SPE 77603
[8]曹睿,余剛,陳劍虹.鋁合金和鋼異種金屬CMT焊接性分析[C].第十六次全國焊接學術會議論文摘要集,2011. CAO Rui,Yu Gang,Chen Jianhong.A CMT welding analysis of Aluminum and steel dissimilar metals[C]. Sixteenth National Welding Conference Abstract Book,2011.
[9]CRUMP J B,CONWAY M W.Effects of perforationentry friction on bottom hole treating analysis[J].Journal of Petroleum Technology,1988,40(8):1041-1048.
[10]EI-RABBA A M ,SHAH S N,LORD D L.New perforation pressure-loss correlations for limited-entry fracturing treatments[J].Spe Production & Facilities,1999,14(1):63-71
[11]劉合,張廣明,張勁,等.油井水力壓裂摩阻計算和井口壓力預測[J].巖石力學與工程學報,2010,29(增1):2833-2835 LIU He,ZHANG Guangming,ZHANG Jin,et al. Friction loss calculation and surface pressure prediction in oil well hydraulic fracturing[J].Chinese Journal of Rock Mechanics and Engineering,2010,29(S1):2833-2835.
[12]LORD L.Turbulent flow of stimulation fluids: an evaluation of friction loss scale-up methods[J].SPE 16889.
[13]王鴻勛,張士誠.水力壓裂設計數(shù)值計算方法[M].北京:石油工業(yè)出版社,1998. WANG Hongxun,ZHANG Shicheng.Numerical methods of hydraulic fracturing design[M].Bei Jing:Petroleum Industry Press,1998.
[14]肖榮鴿,王立洋,鄧志安,等.凝析天然氣管道分層流動相間水力摩阻系數(shù)計算式評價[J].西安石油大學學報,2006,21(3):55-57. XIAO Rongge,WANG Liyang,DENG Zhian,et al. Evaluation of the formulas for calculating the interfacial friction factor of the two-phase stratified flow in horizontal condensate natural gas pipeline[J].Journal of Xi′an Shiyou University,2006,21(3):55-57.
[15]LORD D L, MCGOWEN J M.Real-time treating pressure analysis aided by new correlation[R].New Orleans,Louisiana,1986.
Factors Affecting Friction Loss of Hydraulic Fracturing in Ultra-short Radius Radial Wells and the Calculating Equation Thereof
GONG Diguang, QU Zhanqing, GUO Tiankui, GONG Facheng, TIAN Xuxin, HUANG Zitong
(School of Petroleum Engineering, China University of Petroleum, Qingdao, Shandong 266580)
AbstractThe friction loss testing equipment (composed of friction testing controller,fracturing fuid distributor,speed-controllable screw pump,pipes with varied diameters,high-sensitivity pressure gauge,and electric fow meter) developed by China University of Petroleum (East China) has been used to study the effects of borehole diameters and fracturing parameters on the friction loss of fracturing ultra-short radius radial (URR) wells and to calculate the magnitude of the friction losses of the fracturing fuid. The fow conditions in URR wells were simulated and the friction losses of the fracturing fuids used accurately measured. Analysis of the experiment results shows that factors affecting the friction losses of the fracturing fuids in URR wells are,in the order of decreasing importance,hole diameter,fow rate,viscosity,particle size of proppant,and sand content,and the effects varies in different conditions. An equation for the calculation of friction losses of guar gum fracturing fuid in URR wells has been developed taking into account the aforementioned factors,through the regression fts of 322 sets of data and based on the principle of reducing friction ratios. The standard error of estimate calculated from correlation coeffcient testing method is 0.140,indicating that the regression ft equation is valid. It has been found through the experiments that the viscosity of the guar gum fracturing fuid has dual effects on the friction; the viscosity,by increasing the intrinsic shear stress and the shear stress between fuid and the wall of the pipe,increases the friction losses of the fracturing fuid. On the other hand,with an increase in viscosity,a transition delay is being developed in the polymer solution,and the polymer solution has stronger control over the proppants used,thereby decreases the friction losses of the fracturing fuid.
Key wordsFriction loss testing; Factors affecting friction loss; Ultra-short radius radial well; Reduce friction ratio; Hydraulic fracturing;Fracturing fuid; Orthogonal experiment
中圖分類號:TE357.12
文獻標識碼:A
文章編號:1001-5620(2016)03-0102-05
doi:10.3696/j.issn.1001-5620.2016.03.021
基金項目:國家自然科學基金青年基金項目“徑向鉆孔引導水力壓裂裂縫定向擴展機理研究”(51404288);大學生創(chuàng)新創(chuàng)業(yè)訓練計劃項目“徑向井水力壓裂降摩阻技術研究”(201510425014);中國石油大學(華東)研究生創(chuàng)新工程資助項目(YCX2014010)。
第一作者簡介:龔迪光,1983年生,男,中國石油大學(華東)在讀博士研究生,現(xiàn)主要從事油氣田開發(fā)研究工作。電話 17864229383;E-mail:380752913@qq.com。
收稿日期(2016-1-9;HGF=1603F8;編輯付玥穎)