席 巖，李 軍，柳貢慧，2，付永強(qiáng)，李玉梅

(1.中國石油大學(xué)(北京)，北京 102249；2.北京工業(yè)大學(xué)，北京 100124；3.中國石油西南油氣田分公司，四川 成都 610051)

?

頁巖儲層各向異性對套管應(yīng)力影響敏感性研究

席 巖1，李 軍1，柳貢慧1，2，付永強(qiáng)3，李玉梅1

(1.中國石油大學(xué)(北京)，北京 102249；2.北京工業(yè)大學(xué)，北京 100124；3.中國石油西南油氣田分公司，四川 成都 610051)

頁巖儲層各向異性特征明顯，基于巖石橫向各向同性本構(gòu)關(guān)系理論，建立了套管—水泥環(huán)—地層組合體三維有限元模型，對頁巖儲層各向異性及其與非均勻地應(yīng)力、水泥石力學(xué)特性、水泥環(huán)形態(tài)的耦合作用對套管應(yīng)力的影響進(jìn)行研究。研究結(jié)果表明：頁巖儲層各向異性增大了套管應(yīng)力，增加了套管損壞的風(fēng)險；地應(yīng)力非均勻度越高，儲層各向異性對于套管應(yīng)力的影響越大；水泥石力學(xué)特性不變時，考慮儲層各向異性時套管應(yīng)力比不考慮時高；套管偏心時，在套管應(yīng)力較小的區(qū)域，考慮各向異性時比不考慮各向異性時更小，在套管應(yīng)力較大的區(qū)域，考慮各向異性時要比不考慮各向異性時更大。研究結(jié)果對于解決多級壓裂過程中的套管損壞問題，提高頁巖氣開發(fā)產(chǎn)能具有重要的理論和應(yīng)用價值。

頁巖氣；各向異性；地應(yīng)力；有限元；套管應(yīng)力

0 引 言

威遠(yuǎn)—長寧地區(qū)頁巖氣井多級壓裂過程中套管損壞問題突出。部分專家學(xué)者就改造區(qū)域不對稱、施工壓力大、井筒環(huán)空束縛流體收縮等多個影響因素進(jìn)行了分析[1-6]，但是均未考慮頁巖儲層各向異性對套管應(yīng)力的影響。多級壓裂過程中套管處于復(fù)雜、極端的力學(xué)環(huán)境中，套管內(nèi)壓持平或者已經(jīng)超過其抗內(nèi)壓強(qiáng)度。對套管應(yīng)力影響因素考慮是否全面對于明確套管損壞的主控因素具有重要影響。因此，研究頁巖儲層各向異性對套管應(yīng)力影響的敏感性就十分必要。考慮層理性頁巖橫向各向同性，結(jié)合實(shí)驗(yàn)數(shù)據(jù)，研究頁巖彈性各向異性與地應(yīng)力非均勻性、水泥石力學(xué)特性、水泥環(huán)形態(tài)的耦合作用對套管應(yīng)力的影響，可為頁巖氣水平井多級壓裂過程中套管損壞問題的分析提供指導(dǎo)意見。

1 層理性頁巖各向異性本構(gòu)模型

頁巖具有明顯的層理性特征[7-15]，受這種特征影響，儲層平行層理和垂直層理表現(xiàn)出不同的力學(xué)效應(yīng)，在工程上通常被稱為橫觀各向同性，可以用5個獨(dú)立的彈性常數(shù)來表示。

當(dāng)XZ平面為橫觀各向同性平面時，對于各向異性材料，其應(yīng)力應(yīng)變關(guān)系為[16]：

ε=Sσ

(1)

(2)

式中：ε為應(yīng)變張量；S為柔度矩陣；σ為應(yīng)力張量，GPa；Ex、Ez為平行于各向同性面的彈性模量， GPa；Ey為垂直于各向同性面的彈性模量，GPa；vxz為平行于各向同性面的泊松比；vxy、vyz為垂直于各向同性面的泊松比。

當(dāng)XZ平面為橫觀各向同性平面時，Ex=Ez=Eh，Ey=Ev，vxy=vyx=vyz=vzy=vv，vxz=vzx=vh。Eh、Ev為橫觀各向同性面和垂直于橫觀各向同性面的彈性模量，GPa；vh、vv為橫觀各向同性面和垂直于橫觀各向同性面的泊松比。

各向同性面內(nèi)的剪切模量為：

(3)

式中：Gzx為平行于各向同性面的剪切模量，GPa；Gh為橫觀各向同性面的剪切模量，GPa。

Batugin和Nirenburg等人提出了垂直于各向同性面的Gyz=Gxy=Gv的數(shù)學(xué)解法，并且通過大量的實(shí)驗(yàn)數(shù)據(jù)對其進(jìn)行了驗(yàn)證[17]：

(4)

式中：Gyz、Gxy為垂直于各向同性面的剪切模量，GPa；Gv為垂直于橫觀各向同性面的剪切模量，GPa。

式(2)中最終只剩下Eh、Ev、vh、vv、Gv5個獨(dú)立的彈性參數(shù)。

2 有限元數(shù)值模型

頁巖氣水平井井眼軌跡一般都是沿著最小主應(yīng)力方向?？紤]井眼方向以及層理性頁巖橫向各向同性特征，建立套管—水泥環(huán)—地層組合體三維有限元模型，再建立局部坐標(biāo)系，局部坐標(biāo)系X′Y′Z′坐標(biāo)軸設(shè)置采取與全局坐標(biāo)系坐標(biāo)軸XYZ一一對應(yīng)方式，XZ、X′Z′面位于各向同性面內(nèi)。

模型大小為3 m×3 m×3 m，采用結(jié)構(gòu)網(wǎng)格，網(wǎng)格劃分時，為減小計(jì)算干擾，采用變密度網(wǎng)格劃分方法?；A(chǔ)力學(xué)參數(shù)以頁巖氣井W1井實(shí)際工程數(shù)據(jù)為準(zhǔn)，獨(dú)立彈性參數(shù)根據(jù)本構(gòu)模型計(jì)算得出(表1)。套管選用N80剛級，壁厚為9.17 mm的套管為研究對象。地層最大主應(yīng)力為48 MPa，地層最小主應(yīng)力為29 MPa，地層垂向主應(yīng)力為35 MPa。

表1 套管—水泥環(huán)—地層組合體幾何及力學(xué)參數(shù)

3 儲層各向異性對套管應(yīng)力影響敏感性分析

3.1 彈性模量與泊松比各向異性對套管應(yīng)力影響分析

頁巖儲層各向異性包含彈性模量各向異性和泊松比各向異性。定義頁巖層理面彈性模量與法向彈性模量的比值為K，K=Eh/Ev，層理面泊松比與層理面法向泊松比的比值為K′，K′=vh/vv。K和K′均為表征材料各向異性的參考值。圖1、2為套管最大Mise應(yīng)力隨彈性模量各向異性和泊松比各向異性的變化曲線。由圖1可知：K值不斷增大，套管承受應(yīng)力不斷增大，增速逐漸變緩；K>1時，彈性模量各向異性越強(qiáng)，套管應(yīng)力越大；K′值不同時，3條曲線幾乎發(fā)生重合，說明K′值對套管應(yīng)力的影響較小。由圖2可知：K′值不斷增大時，套管承受應(yīng)力不斷減小，但是變化幅度相對較小。

圖1 套管最大Mise應(yīng)力隨彈性模量各向異性變化曲線

采用重慶石柱縣黑色寒武系頁巖對頁巖巖石力學(xué)特性參數(shù)開展實(shí)驗(yàn)研究。實(shí)驗(yàn)結(jié)果證明：隨著圍壓不斷增大，頁巖彈性模量不斷增大，Eh相比Ev增幅更大。其主要原因?yàn)椋喉搸r存在明顯層狀節(jié)理，圍壓作用時，層理面間的孔隙和微裂隙被壓實(shí)，抑制了側(cè)向變形，而在垂直層理方向，層理間孔隙和微裂隙的壓實(shí)作用受圍壓影響較小，導(dǎo)致圍壓越大，K值越大。

圖2 套管最大Mise應(yīng)力隨泊松比各向異性變化曲線

該規(guī)律與衡帥[18]等人得到的規(guī)律較為相似，根據(jù)其實(shí)驗(yàn)數(shù)據(jù)可知：當(dāng)井下圍壓達(dá)到30 MPa時，K值達(dá)到2.348。在該圍壓條件下，K′也達(dá)到了1.65。由于泊松比各向異性對套管應(yīng)力的影響較小，因此，選擇K=2.348，K′=1.00作為多因素作用分析時的基礎(chǔ)計(jì)算條件。

3.2 多因素作用下各向異性對套管應(yīng)力影響敏感性分析

套管在井下受非均勻地應(yīng)力、水泥石力學(xué)特性、水泥環(huán)形態(tài)等多種因素的影響，這些因素與儲層各向異性會發(fā)生耦合作用。明確耦合作用對套管應(yīng)力的影響規(guī)律，對于精確計(jì)算復(fù)雜、極端力學(xué)環(huán)境下的套管應(yīng)力和準(zhǔn)確判斷套管是否發(fā)生屈服破壞極為必要。

3.2.1 非均勻地應(yīng)力與各向異性耦合影響分析

水平井最大水平地應(yīng)力與垂向地應(yīng)力差值越大，井筒受力就越不均勻，套管損壞的風(fēng)險就越大[19]。定義Kf(Kf=FH/Fv)為水平和垂向地應(yīng)力非均勻度，其中，F(xiàn)H為最大水平地應(yīng)力，MPa；Fv為垂向地應(yīng)力，MPa。當(dāng)最大水平地應(yīng)力保持不變時，改變Fv以獲得不同的Kf。

圖3為考慮(不考慮)各向異性時地應(yīng)力非均勻度對于套管最大Mise應(yīng)力影響曲線。由圖3可知：隨著地應(yīng)力非均勻度的不斷增加，套管最大Mise應(yīng)力不斷增加，增加趨勢不斷變緩；Kf=1時，儲層各向異性帶來的影響較小；Kf>1時，考慮地層各向異性時套管最大Mise應(yīng)力比同等條件下不考慮地層各向異性時要大，2種條件下，套管最大Mise應(yīng)力隨著地應(yīng)力非均勻度的不斷增加，差值不斷增大。

圖3 考慮(不考慮)各向異性時非均勻地應(yīng)力對套管最大Mise應(yīng)力的影響

3.2.2 水泥石力學(xué)特性與各向異性耦合影響分析

水泥石的力學(xué)特性對于套管應(yīng)力大小有顯著影響，高強(qiáng)度、低剛度的水泥石力學(xué)特性對于套管保護(hù)最為有利[20]。圖4為考慮(不考慮)各向異性時水泥石彈性模量變化對套管最大Mise應(yīng)力影響曲線。由圖4可知：水泥石彈性模量不斷增大，套管最大Mise應(yīng)力不斷減??；考慮各向異性情況時應(yīng)力比不考慮時要大，同等力學(xué)條件下套管損壞的風(fēng)險也更大；2種條件下套管應(yīng)力差值隨著水泥環(huán)彈性模量的增大不斷減小。圖5為考慮(不考慮)各向異性時水泥石泊松比變化對套管最大Mise應(yīng)力影響曲線。由圖5可知：水泥石泊松比不斷增大，套管應(yīng)力不斷減?。豢紤]頁巖儲層各向異性時，套管應(yīng)力比不考慮時要大；2種條件下套管最大Mise應(yīng)

圖4 考慮(不考慮)各向異性時水泥石彈性模量對套管最大Mise應(yīng)力的影響

圖5 考慮(不考慮)各向異性時水泥石泊松比對套管最大Mise應(yīng)力的影響

力差幾乎不變。

3.2.3 水泥環(huán)形態(tài)與各向異性耦合影響分析

套管偏心和水泥環(huán)缺失是固井過程中常見的水泥環(huán)形態(tài)。圖6為考慮(不考慮)各向異性時偏心角對套管周向應(yīng)力的影響(偏心角以井眼水平方向右側(cè)為0 °，逆時針方向增加)。由圖6可知：偏心導(dǎo)致套管周向應(yīng)力左右不對稱，各向異性則帶來了“放大效應(yīng)”，在套管應(yīng)力較小的區(qū)域，考慮各向異性時比不考慮各向異性更小，如圖6中藍(lán)色標(biāo)識區(qū)域；在套管應(yīng)力較大的區(qū)域，考慮各向異性時要比不考慮各向異性更大。圖7為偏心角一定時，考慮(不考慮)各向異性時偏心距對套管應(yīng)力的影響。由圖7可知：隨著偏心距的增大，套管最大Mise應(yīng)力不斷增大；偏心距一定時，考慮各向異性時要比不考慮各向異性時應(yīng)力值更大，且偏心距越大，2種條件下的套管應(yīng)力差值越大。

圖6 偏心角對套管周向應(yīng)力的影響(偏心距為18mm)

為進(jìn)一步明確各向異性與水泥環(huán)缺失帶來的耦合影響。選定偏心角為0 °、偏心距為18 mm作為基礎(chǔ)計(jì)算條件，計(jì)算缺失角從0 °到180 °時套管最大Mise應(yīng)力的變化(圖8)。由圖8可知，缺失角不斷增大，套管最大Mise應(yīng)力先增大后減小。考慮儲層各向異時要比不考慮儲層各向異性時套管最大Mise應(yīng)力更大，但是兩者差值幾乎不隨缺失角的改變而改變。

圖7 套管偏心距對套管最大Mise應(yīng)力的影響(偏心角為0°)

圖8 水泥環(huán)缺失角對套管最大Mise應(yīng)力影響

4 結(jié) 論

(1) 建立了考慮頁巖儲層各向異性的套管—水泥環(huán)—地層有限元模型，研究發(fā)現(xiàn)在一定井筒力學(xué)條件下，儲層各向異性提高了套管應(yīng)力。彈性模量各向異性對于套管應(yīng)力大小比較敏感，K值越大，套管最大Mise應(yīng)力也越大。

(2) 在考慮非均勻地應(yīng)力、水泥石力學(xué)特性、水泥環(huán)形態(tài)工程實(shí)際情況的基礎(chǔ)上，建立數(shù)值模型，計(jì)算了這些因素與頁巖儲層各向異性的耦合影響。計(jì)算結(jié)果表明：地應(yīng)力非均勻度越高，儲層各向異性對于套管最大Mise應(yīng)力的影響越大；水泥石力學(xué)特性相同時，考慮儲層各向異性時套管最大Mise應(yīng)力要比不考慮時大，同種條件下的套管損壞風(fēng)險也更大；套管偏心時，儲層各向異性使得套管周向上應(yīng)力小處更小、大處更大。

(3) 頁巖氣多級壓裂過程中，井筒處于較為復(fù)雜、極端的力學(xué)環(huán)境中，工程上建議在套管損壞分析過程中對于頁巖儲層各向異性予以考慮，準(zhǔn)確計(jì)算套管所受應(yīng)力，正確判斷套管損壞發(fā)生的原因。

[1] 滕柏路，程林松，黃世軍，等.考慮近井流動的頁巖氣壓裂水平井理論模型[J].大慶石油地質(zhì)與開發(fā)，2014，33(3)：158-163.

[2] 劉新，張玉瑋，鐘顯東，等.美國典型頁巖油氣藏水平井壓裂技術(shù)[J].大慶石油地質(zhì)與開發(fā)，2014，33 (6)：160-164.

[3] 李志強(qiáng)，趙金洲，胡永全，等.頁巖氣壓裂水平井產(chǎn)能模擬與布縫模式[J].大慶石油地質(zhì)與開發(fā)，2015，34 (5)：162-165.

[4] 尹帥，丁文龍，王濡岳，等.陸相致密砂巖及泥頁巖儲層縱橫波波速比與巖石物理參數(shù)的關(guān)系及表征方法[J].油氣地質(zhì)與采收率，2015，22(3)：22-28.

[5] 田中蘭，石林，喬磊.頁巖氣水平井井筒完整性問題及對策[J].天然氣工業(yè)，2015，35(9)：70-77.

[6] LIAN Zhanghua， YU Hao， LIN Tiejun，et al. A study on casing deformation failure during multi-stage hydraulic fracturing for the stimulated reservoir volume of horizontal shale wells [J]. Journal of Natural Gas Science and Engineering，2015，23(1)：538-546.

[7] TALIERCIO A，LANDERIANI G S.A failure condition for layered rock[J]. International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts，1988，25(5)：299-305.

[8] 李進(jìn)步，盧雙舫，陳國輝，等.基于礦物學(xué)和巖石力學(xué)的泥頁巖儲層可壓裂性評價[J].大慶石油地質(zhì)與開發(fā)，2015，34 (6)：159-164.

[9] 鄒孟飛，隋微波，王旭東，等.基于非?？焖倌M退火法的頁巖巖心雙重區(qū)域重構(gòu)方法[J].油氣地質(zhì)與采收率，2015，22(5)：117-122.

[10] 索彧，汪道兵， 左乾坤，等.單軸加載頁巖破壞全過程的聲發(fā)射特征[J].大慶石油地質(zhì)與開發(fā)，2016，35 (2)：159-164.

[11] 康毅力，陳強(qiáng)，游利軍，等.頁巖氣藏水相圈閉損害實(shí)驗(yàn)研究及控制對策——以四川盆地東部龍馬溪組露頭頁巖為例[J].油氣地質(zhì)與采收率，2014，21(6)：87-91.

[12] 侯振坤，楊春和，郭印同，等.單軸壓縮下龍馬溪組頁巖各向異性特征研究[J].巖土力學(xué)，2015，36(9)：2541-2550.

[13] 王倩，王鵬，項(xiàng)德貴，等.頁巖力學(xué)參數(shù)各向異性研究[J].天然氣工業(yè)，2012，32(12)：62-65.

[14] 王小瓊，葛洪魁，申潁浩.富有機(jī)質(zhì)頁巖波速各向異性的研究進(jìn)展[J].特種油氣藏，2013，20(6)：1-5.

[15] 唐杰.各向異性巖石的靜態(tài)模量與動態(tài)模量實(shí)驗(yàn)研究[J].巖石力學(xué)與工程學(xué)報(bào)，2014，33(1)：3185-3191.

[16] 沈觀林，胡更開.復(fù)合材料力學(xué)[M].北京：清華大學(xué)出版社，2006：28-47.

[17] 李玉梅，柳貢慧，李軍，等.頁巖儲層各向異性對裂縫起裂壓力的影響[J].特種油氣藏，2014，21(6)：133-137.

[18] 衡帥，楊春和，張保平，等. 頁巖各向異性特征的試驗(yàn)研究[J].巖土力學(xué)，2015，36(3)：609-616.

[19] 殷有泉，蔡永恩，陳朝偉，等.非均勻地應(yīng)力場中套管載荷的理論解[J].石油學(xué)報(bào)，2006，27(4)：133-138.

[20] 王耀鋒，李軍強(qiáng)，楊小輝，等.水泥環(huán)彈性模量和泊松比與地層性質(zhì)匹配關(guān)系研究[J].石油鉆探技術(shù)，2008，36(6)：25-29.

編輯 孟凡勤

20160620；改回日期：20160822

中國石油天然氣股份有限公司西南油氣田分公司項(xiàng)目“威遠(yuǎn)長寧頁巖氣水平井固井質(zhì)量對井筒完整性的影響”(XNS21JS2014-04)、“大興場構(gòu)造鉆完井工藝技術(shù)研究(實(shí)驗(yàn)檢驗(yàn))”(XNS05JS2015-37)

席巖(1985-)，男，2008年畢業(yè)于中國石油大學(xué)(北京)石油工程專業(yè)，現(xiàn)為該校油氣井工程專業(yè)在讀博士研究生，從事巖石力學(xué)、井筒完整性方面的研究。

10.3969/j.issn.1006-6535.2016.06.029

TE21

A

1006-6535(2016)06-0128-05