尹　虎 張韻洋

西南石油大學(xué)石油與天然氣工程學(xué)院

?

溫度作用影響套管抗擠強(qiáng)度的定量評價方法
——以頁巖氣水平井大型壓裂施工為例

尹虎 張韻洋

西南石油大學(xué)石油與天然氣工程學(xué)院

尹虎等. 溫度作用影響套管抗擠強(qiáng)度的定量評價方法——以頁巖氣水平井大型壓裂施工為例.天然氣工業(yè), 2016,36(4):73-77.

摘 要目前通常采用大規(guī)模水力壓裂的方法對頁巖氣儲層進(jìn)行改造，由于施工排量大、施工時間長，造成井筒溫度下降速度快且幅度大，由溫度作用而產(chǎn)生的拉應(yīng)力會降低套管抗擠強(qiáng)度，不利于安全生產(chǎn)。針對頁巖氣水平井壓裂施工工藝及工況條件，建立了水力壓裂時的熱傳導(dǎo)控制方程，并利用有限差分法對其進(jìn)行求解，分別探討了不同井口溫度下常規(guī)水力壓裂和大規(guī)模水力壓裂時井底溫度的變化情況；然后再根據(jù)API套管擠毀壓力計算公式求得了不同鋼級套管在不同溫度變化時套管抗擠強(qiáng)度的變化值。結(jié)果表明：①較之于常規(guī)水力壓裂，大型水力壓裂時，井底溫度變化對套管抗擠強(qiáng)度影響較大且對不同鋼級的套管抗擠強(qiáng)度的影響程度也有所不同，其中高鋼級套管的抗擠強(qiáng)度受影響較小，而低鋼級套管抗擠強(qiáng)度受影響則較大；②當(dāng)在冬季進(jìn)行施工及水平井垂深更深時，這種由于溫度作用造成套管抗擠強(qiáng)度的降低更為明顯。結(jié)論認(rèn)為，在對頁巖氣水平井大規(guī)模水力壓裂時，必須考慮因為井底溫度變化所造成的套管強(qiáng)度降低的問題，并進(jìn)行合理的套管設(shè)計。

關(guān)鍵詞頁巖氣 水平井 大規(guī)模水力壓裂 井筒溫度 熱應(yīng)力 套管 抗擠強(qiáng)度 定量評價 套管設(shè)計

中國頁巖氣分布非常廣泛，開發(fā)潛力巨大，具有極高的經(jīng)濟(jì)價值[1]。但頁巖氣儲層具有低孔隙度和極低基質(zhì)滲透率的特征，頁巖氣水平井單井產(chǎn)量低，必須進(jìn)行大規(guī)模水力壓裂從而形成復(fù)雜縫網(wǎng)以提高產(chǎn)量[2]，與常規(guī)水力壓裂相比，套管直接與壓裂液接觸，承受高內(nèi)壓和變溫作用，套管力學(xué)環(huán)境更為復(fù)雜。同時，這種壓裂方式需要的排量大且作業(yè)時間長，造成在施工周期內(nèi)井底溫度下降幅度大，下降速度較快[3-4]。而井底溫度大幅下降，套管在熱脹冷縮作用下，會產(chǎn)生軸向拉力，導(dǎo)致套管抗擠強(qiáng)度降低[5]。因此，定量分析頁巖氣套管大規(guī)模壓裂工況下，溫度變化對套管抗擠強(qiáng)度的影響程度，對合理設(shè)計頁巖氣水平井套管強(qiáng)度具有重要意義。

1 套管壓裂井筒溫度場模型

1.1 基本方程

近年來，國內(nèi)外學(xué)者建立了很多模型來模擬井筒溫度場的分布規(guī)律[6-9]。Ramey Jr[10]引入了總導(dǎo)熱系數(shù)，計算了長時間處于穩(wěn)定狀態(tài)下的井筒溫度變化。Eickmeimer等[11]建立了計算井筒傳熱的數(shù)值方法，該方法可以模擬在較短時間內(nèi)注入液體井筒溫度的變化。王鴻勛等[12]改進(jìn)了Eickmeimer的方法，建立了井筒不穩(wěn)態(tài)傳熱數(shù)值計算方法。苗和平等[13]針對由于在水平井中斜井段和水平段的存在造成整個井筒溫度分布不一致的特點，分別建立了斜井段和水平段溫度場的數(shù)值計算模型。這些模型建立是主要考慮采用油管壓裂時候井筒溫度場分布情況，與目前頁巖氣水平井采用套管壓裂存在一定的差異。

對于頁巖氣水平井大規(guī)模水力壓裂，壓裂液直接作用在套管上，其在井筒中的循環(huán)過程為：①由地面進(jìn)入套管，經(jīng)套管向下流動；②壓裂液在水平段進(jìn)入地層。

進(jìn)行下列假設(shè)以建立合適的熱傳導(dǎo)方程[14]：

1）壓裂液熱容、密度、熱傳導(dǎo)系數(shù)等參數(shù)隨溫度的變化程度很小。

2）流體中的軸向熱傳導(dǎo)與軸向?qū)α飨啾群雎圆挥嫛?/p>

3）由于套管內(nèi)的流體基本處于紊流狀態(tài)，因此假設(shè)套管內(nèi)的流體沒有徑向穩(wěn)定梯度。

4）由流體黏性耗散產(chǎn)生的熱量忽略不計。

5）地層溫度與深度呈線性關(guān)系，即

式中Tz表示地層某一點的溫度，℃；Tb表示地層恒溫點處的溫度，℃；α表示地溫梯度，℃/1 000 m； z表示地層的某一點的深度，m；b表示地層恒溫點處的深度，m。

在套管內(nèi)[15]有：

式中Ac表示套管橫截面積，m2；ρ表示壓裂液密度，g/cm3；Vc表示套管內(nèi)流體的流速，m/s；Cp表示壓裂液比熱，J/（kg·℃）；Tc表示套管內(nèi)溫度，℃；ri表示套管內(nèi)徑，m；U表示流固表面對流換熱系數(shù)，w/（m2·℃）；Tca表示套管壁溫度，℃； t表示壓裂時間，s。

對于套管壁，水泥環(huán)和地層的熱傳導(dǎo)方程可以簡寫為：

式中K表示熱傳導(dǎo)系數(shù)，w/（m·℃）；r表示徑向長度，m。當(dāng)ri＜r＜rca時，K取套管熱傳導(dǎo)系數(shù)Kca，ρ取套管密度ρca，Cp取套管比熱Cpca，此時為套管熱傳導(dǎo)方程；同理當(dāng)rca≤r＜rB時，K取水泥環(huán)熱傳導(dǎo)系數(shù)Kce，ρ取水泥環(huán)密度ρce，Cp取水泥環(huán)比熱Cpce，此時為水泥環(huán)熱傳導(dǎo)方程；當(dāng)rB≤r時，K取地層熱傳導(dǎo)系數(shù)Kf，ρ取地層密度ρf，Cp取地層比熱Cpf，此時為地層熱傳導(dǎo)方程。其中，rca表示套管外徑；rB表示井眼半徑。

由于存在斜井段和水平段，因此整個井筒溫度場變化與垂直段井筒溫度場變化不同。需要分別討論垂直段，斜井段和水平段的邊界條件及初始條件。

1.1.1 垂直段

壓裂液入口溫度已知：

遠(yuǎn)場地層溫度等于初始同一深度地層溫度：

垂直段初始條件：假設(shè)在初始時刻地層，水泥環(huán)和套管的溫度為原始地層溫度，即

式中Tc表示壓裂液溫度；Tf表示地層溫度；Tca表示套管溫度；Tce表示水泥環(huán)溫度。

1.1.2 斜井段

針對斜井段，需要將井深轉(zhuǎn)化為垂深，此時邊界條件：

斜井段遠(yuǎn)場地層溫度等于同一垂深的地層溫度：

斜井段初始條件：在初始時刻地層，水泥環(huán)和套管的溫度等于同一垂深時原始地層溫度，即

1.1.3 水平段

水平段遠(yuǎn)場地層溫度等于同一垂深地層溫度（當(dāng)r=a時）：

水平段初始條件：在初始時刻地層，水泥環(huán)和套管的溫度等于油層中部的溫度，即

采用有限差分方法對控制方程進(jìn)行數(shù)值求解。

1.2 計算實例

以某水平井為例進(jìn)行分析，水平井垂深2 400 m，地溫梯度25 ℃/km。分別計算在常規(guī)壓裂和大規(guī)模壓裂情況下井底溫度變化。一般情況下，常規(guī)水力壓裂排量介于2.5～3.5 m3/min，壓裂施工時間介于30～60 min；而頁巖氣水平井大規(guī)模水力壓裂排量介于10～15 m3/min[16-17]，壓裂施工時間通常大于60 min。根據(jù)以上熱傳導(dǎo)方程，對不同井口溫度下（20℃和0 ℃，分別對應(yīng)在夏季和冬季進(jìn)行水力壓裂施工），大規(guī)模和常規(guī)水力壓裂時的井底溫度變化進(jìn)行了模擬，計算參數(shù)如表1所示，計算結(jié)果如圖1和圖2所示。

表1　井底溫度計算參數(shù)表

圖1　不同井口溫度在大型水力壓裂時井底溫度變化圖

圖2　不同井口溫度在普通水力壓裂時井底溫度變化圖

從模擬結(jié)果可以看出：大排量水力壓裂施工時，當(dāng)壓裂液溫度為20 ℃時，井底溫度從87 ℃下降到34 ℃；壓裂液溫度為0 ℃時，井底溫度從87 ℃下降到20 ℃。大排量水力壓裂施工時，井底溫度能在30 min以內(nèi)迅速達(dá)到熱平衡。常規(guī)水力壓裂施工時，井底溫度下降較慢，通常在30 min之后，仍未達(dá)到熱平衡狀態(tài)。因此大規(guī)模壓裂時，井底溫度下降速度較快，在施工周期內(nèi)井底壓力下降幅度大。

2 溫度作用下套管抗擠強(qiáng)度的變化

當(dāng)套管溫度迅速下降時，由于熱脹冷縮作用，在套管的兩端會產(chǎn)生軸向的內(nèi)應(yīng)力。如果產(chǎn)生的軸向內(nèi)應(yīng)力較大，則會明顯降低套管的強(qiáng)度。

2.1 溫度對軸向應(yīng)力的影響

當(dāng)溫度變化時，鋼的長度變化為：

式中?L表示長度變化值，m；β表示鋼的熱膨脹系數(shù)，12.45×10-6，1/℃；L表示長度，m；?t表示溫度變化值，℃。

當(dāng)鋼材固定，不讓其收縮，在套管柱內(nèi)就會產(chǎn)生內(nèi)應(yīng)力[18]。由胡克定律得：

式中E表示鋼的彈性模量，MPa。

由式（11）可知，由溫度變化產(chǎn)生的軸向應(yīng)力與套管的長度沒有關(guān)系，僅僅取決于溫度的變化值。

2.2 雙軸應(yīng)力下的套管抗擠強(qiáng)度變化

對于彈性材料，最常用的套管破壞強(qiáng)度理論是米塞斯最大應(yīng)變能理論。因此API公式就是基于這一破壞準(zhǔn)則推導(dǎo)出來的。當(dāng)套管發(fā)生擠毀時，根據(jù)米塞斯方程：

當(dāng)徑向應(yīng)力S3=0時等式（12）兩端同時除以Yp2，則有：

用(pca/pco)替代上式中的(S2/Yp)，(Sa/Yp)替代S1/Yp)[19]，則有：

解出式（14），得到：

式中pca表示軸向載荷下的抗擠毀壓力（強(qiáng)度），MPa；pco表示無軸向載荷下的抗擠毀壓力（強(qiáng)度），MPa；Sa表示軸向拉伸載荷的應(yīng)力，MPa；Yp表示管體的屈服強(qiáng)度，MPa。

由此根據(jù)式（15）計算出不同溫度變化下，不同鋼級的套管抗擠強(qiáng)度變化如圖3所示。

圖3　不同溫度變化情況下套管抗擠強(qiáng)度變化趨勢圖

從計算結(jié)果可以看出：不同鋼級的套管在不同的溫度變化時，抗擠強(qiáng)度降低值不一樣，具體表現(xiàn)為較高鋼級套管抗擠強(qiáng)度受溫度變化影響較小，而較低鋼級套管抗擠強(qiáng)度受溫度變化影響較大。當(dāng)溫度變化達(dá)到70 ℃時，TP140套管抗擠強(qiáng)度為原始強(qiáng)度的90％，而C90和N80套管抗擠強(qiáng)度僅為原始強(qiáng)度的80％，溫度變化對套管抗擠強(qiáng)度的影響不可忽視。

仍以某頁巖氣水平井為例進(jìn)行分析，套管為P110鋼級，在冬季進(jìn)行了大規(guī)模水力壓裂時，施工作業(yè)排量為10 m3/min，其井底溫度降低值70 ℃，套管抗擠強(qiáng)度降低幅度14％。

3 結(jié)論

與傳統(tǒng)水力壓裂相比，大規(guī)模水力壓裂時井底溫度迅速下降，在較短的時間內(nèi)溫度下降可超過50℃。井底溫度下降時，由于溫度效應(yīng)在套管內(nèi)產(chǎn)生軸向應(yīng)力會降低套管的抗擠強(qiáng)度。溫度下降越大，對套管抗擠強(qiáng)度影響越顯著：溫度下降超過70 ℃時，高鋼級套管（TP140）的抗擠強(qiáng)度只有初始強(qiáng)度的90％；低鋼級的套管（N80）的抗擠強(qiáng)度下降到初始強(qiáng)度的80％。由此可知，如水平井垂深更深，且在冬季施工時，套管抗擠強(qiáng)度下降更為明顯。在頁巖氣水平井大規(guī)模水力壓裂時，必須考慮因為井底溫度變化所造成的套管抗擠強(qiáng)度降低問題，進(jìn)行合理的套管強(qiáng)度設(shè)計。

參 考 文 獻(xiàn)

[1] 劉洪林, 王紅巖, 劉人和, 趙群, 林英姬. 中國頁巖氣資源及其勘探潛力分析[J]. 地質(zhì)學(xué)報, 2010, 84(9): 1374-1378. Liu Honglin, Wang Hongyan, Liu Renhe, Zhao Qun, Lin Yingji. China shale gas resources and prospect potential[J]. Acta Geologica Sinica, 2010, 84(9): 1374-1378.

[2] 董大忠, 鄒才能, 楊樺, 王玉滿, 李新景, 陳更生, 等. 中國頁巖氣勘探開發(fā)進(jìn)展與發(fā)展前景[J]. 石油學(xué)報, 2012, 33(增刊1): 107-114. Dong Dazhong, Zou Caineng, Yang Hua, Wang Yuman, Li Xinjing, Chen Gengsheng, et al. Progress and prospects of shale gas exploration and development in China[J]. Acta Petrolei Sinica, 2012, 33(S1): 107-114.

[3] 葉靜, 胡永全, 葉生林, 王丹. 頁巖氣藏水力壓裂技術(shù)進(jìn)展[J].天然氣勘探與開發(fā), 2012, 35(4): 64-68. Ye Jing, Hu Yongquan, Ye Shenglin, Wang Dan. Technical progress of hydraulic fracturing in shale-gas reservoirs[J]. Natural Gas Exploration and Development, 2012, 35(4): 64-68.

[4] Sugden C, Ring G A, Chambers M R, Suryanarayana PV. Special considerations in the design optimization of high rate, multistage fractured shale wells[C]//IADC/SPE Drilling Conference and Exhibition, 6-8 March 2012, San Diego, California, USA. DOI:http// dx.doi.org/10.2118/151470-MS.

[5] Clark HC. Mechanical design considerations for fracture-treating down casing strings (includes associated papers 17086 and 17124) [J]. SPE Drilling Engineering, 1987, 2(2): 116-126.

[6] Barelli A, Palamà A. A new method for evaluating formation equilibrium temperature in holes during drilling[J]. Geothermics, 1981, 10(2): 95-102.

[7] Wooley GR. Computing downhole temperatures in circulation, injection, and production wells[J]. Journal of Petroleum Technology, 1980, 32(9): 1509-1522.

[8] 盧德唐, 曾億山, 郭永存. 多層地層中的井筒及地層溫度解析解[J]. 水動力學(xué)研究與進(jìn)展: A輯, 2002, 17(3): 382-390. Lu Detang, Zeng Yishan, Guo Yongcun. Analytical solution of temperature in well-bore and formation in multi-layer[J]. Journal of Hydrodynamics, 2002, 17(3): 382-390.

[9] 周趙川, 王輝, 代向輝, 李登, 徐明海. 海上采油井筒溫度計算及隔熱管柱優(yōu)化設(shè)計[J]. 石油機(jī)械, 2014, 42(4): 43-48. Zhou Zhaochuan, Wang Hui, Dai Xianghui, Li Deng, Xu Minghai. Calculation of offshore oil production wellbore temperature and optimal design of heat-proof string[J]. China Petroleum Machinery, 2014, 42(4): 43-48.

[10] Ramey Jr JJ. Wellbore heat transmission[J]. Journal of Petroleum Technology, 1962, 14(4): 427-435.

[11] Eickmeier JR, Ersoy D, Ramey HJ. Wellbore temperatures and heat losses during production or injection operations[J]. Journal of Canadian Petroleum Technology, 1970, 9(2): 115-121.

[12] 王鴻勛, 李平. 水力壓裂過程中井筒溫度的數(shù)值計算方法[J].石油學(xué)報, 1987, 8(2): 91-99. Wang Hongxun, Li Ping. Numerical method for calculating wellbore temperature during hydraulic fracture[J]. Acta Petrolei Sinica, 1987, 8(2): 91-99.

[13] 苗和平, 王鴻勛. 水平井壓裂施工中溫度場的計算方法[J]. 石油鉆采工藝, 1994, 16(3): 45-49. Miao Heping, Wang Hongxun. A method for calculating the wellbore temperature feld in the fracture treatment of horizontal well[J]. Oil Drilling & Production Technology, 1994, 16(3): 45-49.

[14] Gu H. A study of propagation of hydraulically fractures[D]. Austin: The University of Texas, 1987.

[15] Osisanya SO, Harris OO. Evaluation of equivalent circulating density of drilling fluids under high pressure/high temperature conditions[C]//SPE Annual Technical Conference and Exhibition, 9-12 October 2005, Dallas,Texas, USA. DOI:http://dx.doi. org/10.2118/97018-MS.

[16] 陳作, 薛承瑾, 蔣廷學(xué), 秦鈺銘. 頁巖氣井體積壓裂技術(shù)在我國的應(yīng)用建議[J]. 天然氣工業(yè), 2010, 30(10): 30-32. Chen Zuo, Xue Chengjin, Jiang Tingxue, Qin Yuming. Proposals for the application of fracturing by stimulated reservoir volume (SRV) in shale gas wells in China[J]. Natural Gas Industry, 2010, 30(10): 30-32.

[17] 范鑫, 李根生, 黃中偉, 陳龍, 張啟龍, 李浩. 水力噴射多級壓裂中水力封隔效果的數(shù)值模擬[J].石油機(jī)械, 2015, 43(4): 82-88. Fan Xin, Li Gensheng, Huang Zhongwei, Chen Long, Zhang Qilong, Li Hao. Numerical simulation of effect of hydraulic isolation in hydraulic jet multi-stage fracturing[J]. China Petroleum Machinery, 2015, 43(4): 82-88.

[18] 《鉆井手冊》編寫組. 鉆井手冊[M]. 2版. 北京：石油工業(yè)出版社，2013. Drilling Handbook Compile Group. Drilling handbook[M]. 2nded. Beijing: Petroleum Industry Press, 2013.

[19] American Petroleum Institute. API TR 5C3-2008 Technical report on equations and calculations for casing, tubing, and line pipe used as casing or tubing, and performance properties tables for casing and tubing[S]. Washington DC: API, 2008.

(修改回稿日期 2016-01-11 編 輯 凌 忠)

A quantitative evaluation method for the effect of temperature on casing collapsing strength: A case study of large-scale hydraulic fracturing in shale gas horizontal wells

Yin Hu, Zhang Yunyang
(School of Oil and Gas Engineering, Southwest Petroleum University, Chengdu, Sichuan 610500, China)

NATUR. GAS IND. VOLUME 36, ISSUE 4, pp.73-77, 4/25/2016. (ISSN 1000-0976; In Chinese)

Abstract:At present, a large-scale hydraulic fracturing technology is generally used in the stimulation of shale gas reservoirs. Due to high operation discharge and long operation time, however, the wellbore temperature drops quickly and drastically, and consequently tensile stress is generated, resulting in the reduction of casing collapsing strength, which is unfavorable for safe production. In this paper, a heat conduction equation of hydraulic fracturing was established on the basis of the operation conditions and fracturing technologies for shale gas horizontal wells, and then solution was figured out by the finite differential method. Comparative analysis was performed on the bottomhole temperature changes between conventional hydraulic fracturing and large-scale hydraulic fracturing at different wellhead temperatures. A reasonable casing design was made taking into account the impact of bottomhole temperature change on casing strength in the large-scale fracturing shale-gas well. And finally, the collapsing strength changes of casings with different steel grades at different temperature changes were calculated by using the API casing collapsing stress formula. It is shown that the casing collapsing strength is more affected by bottomhole temperature change in large-scale hydraulic fracturing than in conventional hydraulic fracturing, and the impact varies with the steel grades. The casings of high steel grades are less affected, while the effect on those of low steel grades is more. If the operation is carried out in winter, the casing collapsing strength related with the temperature effect decreases significantly when a horizontal well is vertically deeper. It is concluded that when large-scale hydraulic fracturing is performed in shale gas horizontal wells, it is necessary to take into consideration the reduction of casing collapsing strength caused by the change of bottom hole temperature if a proper casing design is made.

Keywords:Shale gas; Horizontal well; Large-scale hydraulic fracturing; Wellbore temperature; Thermal stress; Casing; Collapsing strength; Quantitative evaluation; Casing design

DOI:10.3787/j.issn.1000-0976.2016.04.011

基金項目：國家科技重大專項“鉆井工程設(shè)計和工藝軟件”（編號：2011ZX05021-006）。

作者簡介：尹虎，1978年生，副教授，博士；從事油氣井工程方面的教學(xué)和科研工作。地址：（610500）四川省成都市新都區(qū)新都大道8號西南石油大學(xué)。電話：13808170963。ORCID:0000-0001-5297-0218。E-mail:huyinswpu@outlook.com