曾 靜，高德利(中國科學(xué)院院土)，王宴濱，房 軍

(中國石油大學(xué)石油工程教育部重點實驗室)

在頁巖氣開采中，橋塞分段壓裂工藝需通過壓裂泵向套管內(nèi)輸送壓裂液，此過程中套管和水泥環(huán)受到壓裂液溫度和不斷變化內(nèi)壓力的耦合作用，此交變載荷易導(dǎo)致水泥環(huán)產(chǎn)生微裂縫或微環(huán)隙，為儲層氣體運移至井口提供了通道，最終導(dǎo)致井口套壓不斷升高。國內(nèi)外學(xué)者在井筒完整性方面已開展了相關(guān)研究，比較典型的包括：房軍[1-2]研究了套管-水泥-地層系統(tǒng)在非均勻地應(yīng)力下的應(yīng)力分布規(guī)律；殷有泉[3]給出了非均勻地應(yīng)力下套管-水泥環(huán)力學(xué)特性的理論解；王宴濱[4]對非均勻地應(yīng)力下多層套管柱組合的受載特性進行了研究。許紅林[5]和張智[6]對溫度和內(nèi)壓作用下水泥環(huán)的力學(xué)特性進行了分析；劉奎[7]研究了水泥環(huán)的封隔性能，認(rèn)為第一界面和第二界面是固井失效的危險點。在水泥環(huán)的破壞形式方面，初緯[8]考慮水泥環(huán)塑性變形的影響，研究了變內(nèi)壓條件下微環(huán)隙的形成機理；Goodwin[9]認(rèn)為水泥環(huán)在過高的套管內(nèi)壓和溫變條件下容易形成徑向裂縫；Shadravan[10]進行了高溫和變內(nèi)壓條件下水泥環(huán)破壞實驗，在循環(huán)內(nèi)壓作用下，水泥環(huán)的抗拉強度降低。

綜上所述，現(xiàn)有研究對頁巖氣壓裂過程水泥環(huán)在溫度場和非均勻地應(yīng)力場耦合作用下的力學(xué)特性和失效形式研究不足。一般而言，水泥環(huán)的抗拉強度較低，例如威遠示范區(qū)使用的固井水泥抗拉強度為3.2～4.3 MPa[11]，而水泥環(huán)在體積壓裂過程中承受著較高的內(nèi)壓載荷和溫度變化，容易出現(xiàn)拉伸斷裂失效，是導(dǎo)致環(huán)空帶壓的一個重要因素。

本文重點研究了非均勻地應(yīng)力與熱應(yīng)力耦合作用下由于水泥環(huán)拉伸失效問題，考慮壓裂井筒溫度場的影響，分析了頁巖氣井壓裂過程中水泥環(huán)的環(huán)向應(yīng)力分布規(guī)律，并探討了井筒溫度和地應(yīng)力非均勻程度對井筒完整性的影響規(guī)律。

一、頁巖氣井壓裂井筒傳熱模型

1.假設(shè)條件

(1)忽略井筒套管柱和地層的縱向傳熱。

(2)壓裂液與井壁之間為穩(wěn)態(tài)傳熱，地層內(nèi)傳熱為非穩(wěn)態(tài)傳熱。

(3)套管柱、地層巖石為均質(zhì)各向同性的線彈性材料。

(4)套管柱、地層巖石及壓裂液物性隨溫度和壓力保持不變。

2. 壓裂工況井筒傳熱模型

在井深z處，取體積流量為Q、長度為dz的壓裂液微元體為研究對象，其傳熱模型如圖1所示。

圖1 套管內(nèi)壓裂液微元體熱量傳遞分析模型

如圖1所示，以對流方式帶入和帶出的熱量為Φz和Φz+dz、側(cè)面套管壁的傳熱量為Φk，流體摩擦產(chǎn)生的熱量為Φf。根據(jù)能量守恒定律可得[12]：

(1)

式中：T—壓裂液溫度，℃；Tw—套管內(nèi)壁溫度，℃；r1—套管壁內(nèi)半徑，m；ρm—密度，kg/m3；cm—比熱容，J(kg·℃)；f—摩阻系數(shù)；v—壓裂液流速，m/s；U—井筒總傳熱系數(shù)，W/(m2·℃)。

套管與微元體之間的傳熱量和套管與地層之間的傳熱量相等：

(2)

式中：ke—水泥環(huán)傳熱系數(shù)，W/(m·℃)；Te—地層溫度，℃；f(t)—無因次時間。

壓裂作業(yè)的初始條件和邊界條件為：

(3)

式中：Tin—壓裂液入口溫度，℃。

二、井筒完整性力學(xué)分析模型

套管-水泥環(huán)-地層組合系統(tǒng)力學(xué)模型如圖2所示。組合系統(tǒng)可應(yīng)用平面應(yīng)變理論進行分析,根據(jù)彈性力學(xué)疊加原理對其進行受力分析[13]。

圖 2 套管-水泥環(huán)-地層力學(xué)特性分析模型

1.組合系統(tǒng)各接觸面在熱應(yīng)力下的力學(xué)特性

根據(jù)熱彈性力學(xué)理論，套管-水泥環(huán)-地層組合系統(tǒng)在熱應(yīng)力作用下的應(yīng)力分布為[13]：

(4)

式中：u1,i—第i層管柱的徑向位移，mm；σ1,r,i和σ1,θ,i—第i層管柱的徑向應(yīng)力和環(huán)向應(yīng)力，MPa；Ti—第i個接觸面的溫度，℃；αi—第i層管柱的熱膨脹系數(shù)，1/℃；Ei—第i層管柱的彈性模量，MPa；vi—第i層管柱的泊松比；C1,i和C2,i—待定系數(shù)。

套管-水泥環(huán)-地層組合系統(tǒng)接觸面的溫度為：

(5)

式中：ki—第i層管柱的導(dǎo)熱系數(shù)，W/(m·℃)；ri—第i層管柱內(nèi)半徑，m。

如圖2所示的5層管柱系統(tǒng)的平面應(yīng)變模型，根據(jù)邊界條件和連續(xù)條件可得10個獨立方程。

(1)邊界條件：

(6)

(2)連續(xù)條件：

根據(jù)式(6)和式(7)可求得10個待定系數(shù)，將其帶入式(4)即可求得組合體各接觸面的熱應(yīng)力。

2.組合系統(tǒng)接觸面在非均勻地應(yīng)力下的力學(xué)特性

關(guān)于組合系統(tǒng)各接觸面在非均勻地應(yīng)力下的力學(xué)特性問題已有大量學(xué)者進行研究，組合系統(tǒng)在非均勻地應(yīng)力下應(yīng)力分布的計算方法參見文獻[4]。

3. 組合體應(yīng)力狀態(tài)分析

疊加原理得到組合系統(tǒng)接觸面的總應(yīng)力分布：

(8)

式中：σr,i、σθ,i和σz,i—分別為組合系統(tǒng)第i個接觸面上的徑向應(yīng)力、環(huán)向應(yīng)力和軸向應(yīng)力，MPa；σ2,r,i和σ2,θ,i—分別為組合體在非均勻地應(yīng)力下的徑向應(yīng)力和環(huán)向應(yīng)力，MPa。

三、拉伸失效判據(jù)

采用摩爾-庫侖準(zhǔn)則預(yù)測水泥環(huán)的失效方式[5]。當(dāng)三向應(yīng)力(徑向、環(huán)向、軸向)均處于拉伸狀態(tài)時，其失效準(zhǔn)則為：

σ1≥σt

(9)

式中：σ1—水泥環(huán)的最大主應(yīng)力，MPa；σt—水泥環(huán)的抗拉強度，MPa。

當(dāng)三向應(yīng)力處于拉伸-壓縮-壓縮應(yīng)力狀態(tài)或拉伸-拉伸-壓縮應(yīng)力狀態(tài)時，其失效準(zhǔn)則為：

(10)

式中：σ3—水泥環(huán)的最小主應(yīng)力，MPa；σc—水泥環(huán)的抗壓強度，MPa。

四、套管-水泥環(huán)-地層系統(tǒng)應(yīng)力分析

以某頁巖氣井為例，該井井深為4 120 m，水平段長度為1 364 m，井身結(jié)構(gòu)及計算參數(shù)分別如表1和表2所示。地應(yīng)力梯度為0.023 MPa/m，最大水平地應(yīng)力為36.8 MPa，最小水平地應(yīng)力為27.02 MPa，井口泵壓為100 MPa，水泥環(huán)抗拉強度4.3 MPa。根據(jù)本文介紹的計算方法，計算得到排量為14 m3/min、注入溫度為20℃、注入時間4 h時，井深1 600 m處內(nèi)層水泥環(huán)的應(yīng)力分布如圖3所示。

表1 井身結(jié)構(gòu)

表 2 套管-水泥環(huán)-地層系統(tǒng)物性參數(shù)

圖3 內(nèi)層水泥環(huán)的應(yīng)力分布

由圖3可以看出，水泥環(huán)受非均勻地應(yīng)力場的影響較為顯著，內(nèi)層水泥環(huán)內(nèi)壁的環(huán)向應(yīng)力位于最大水平地應(yīng)力方位，最大值為4.475 MPa，大于其抗拉強度(4.3 MPa)，因此，水泥環(huán)將從內(nèi)壁處最大水平地應(yīng)力方位出現(xiàn)拉伸斷裂并形成徑向裂縫。

1.壓裂液溫度的影響

根據(jù)前述傳熱模型，注入溫度為2℃和25℃時井筒壓裂液溫度分布如圖4所示。從圖4可以看出，當(dāng)注入溫度為2℃，排量由8 m3/min增加至16 m3/min時，井底溫度由46.42℃降低至24.89℃；當(dāng)注入溫度為25℃時，隨著排量增加，井底溫度由54.42℃降低至39.77℃。

不同注入溫度和不同壓裂液排量對水泥環(huán)環(huán)向應(yīng)力的影響如圖5所示。由圖5可知，當(dāng)注入溫度不變時，水泥環(huán)內(nèi)壁的最大環(huán)向應(yīng)力隨排量增加而減小。當(dāng)注入溫度為25℃時，排量應(yīng)大于16 m3/min；當(dāng)注入溫度為2℃時，排量應(yīng)大于8 m3/min。由圖3和圖5可知，當(dāng)排量為14 m3/min，注入溫度為25℃、20℃和2℃時，水泥環(huán)最大環(huán)向應(yīng)力分別為4.771 MPa、4.475 MPa和3.396 MPa，降低注入溫度能降低水泥環(huán)最大環(huán)向應(yīng)力。

圖4 井筒壓裂液溫度分布

圖5 壓裂液排量對內(nèi)層水泥環(huán)環(huán)向應(yīng)力的影響

2. 水平地應(yīng)力非均勻程度的影響

假設(shè)最大水平主應(yīng)力不變，最小水平主應(yīng)力對水泥環(huán)內(nèi)壁的環(huán)向應(yīng)力分布規(guī)律的影響如圖6所示。由圖6可知，隨著水平地應(yīng)力非均勻程度的減小，水泥環(huán)內(nèi)壁的最大、最小環(huán)向應(yīng)力的差值逐漸減小。當(dāng)最小水平地應(yīng)力與最大水平地應(yīng)力相同時，水泥環(huán)內(nèi)壁環(huán)向應(yīng)力小于0。在地應(yīng)力不均勻程度較大的區(qū)塊，水泥環(huán)容易形成徑向裂縫，非均勻地應(yīng)力加劇了水泥環(huán)破壞的風(fēng)險。

圖6 最小水平主應(yīng)力對內(nèi)層水泥環(huán)環(huán)向應(yīng)力分布的影響

五、結(jié)論

(1)根據(jù)熱應(yīng)力、套管內(nèi)壓、非均勻地應(yīng)力耦合條件下套管-水泥環(huán)-地層系統(tǒng)組合體的力學(xué)模型，對水泥環(huán)的環(huán)向應(yīng)力分布和失效位置進行了參數(shù)敏感性分析?？估瓘姸仁怯绊懰嗝芊庑阅艿谋∪跻蛩?，最大水平地應(yīng)力方位的內(nèi)層水泥環(huán)內(nèi)壁最易出現(xiàn)拉伸破壞，是決定井筒完整性的關(guān)鍵位置。

(2)壓裂前后井筒溫差越大和注入溫度越低，水泥環(huán)環(huán)向應(yīng)力越??；地應(yīng)力的非均勻程度越大，水泥環(huán)拉伸破壞的風(fēng)險越大。注入溫度為25℃，排量應(yīng)高于16 m3/min，防止水泥環(huán)拉伸失效形成徑向裂縫。