張 旭， 張哲平， 楊尚諭， 王雪剛， 宋 琳

1天津鋼管制造有限公司 2中國石油集團(tuán)工程材料研究院有限公司 3石油管材及裝備材料服役行為與結(jié)構(gòu)安全國家重點(diǎn)實(shí)驗(yàn)室 4中國石油新疆油田公司工程技術(shù)研究院

0 引言

隨著頁巖氣、煤層氣開采規(guī)模逐漸加大[1- 4]，套損問題變得尤為突出[5- 9]，頁巖氣及煤層氣在壓裂過程中經(jīng)常發(fā)生套管變形現(xiàn)象，損壞了井筒結(jié)構(gòu)的完整性，影響后續(xù)壓裂及生產(chǎn)作業(yè)。因此，摸清套管的真實(shí)服役狀態(tài)，提供安全可靠的管材性能，避免安全系數(shù)過剩，有助于進(jìn)行科學(xué)的固井設(shè)計(jì)選材。

國內(nèi)諸多學(xué)者對(duì)套管的抗外擠額定值做了研究，韓建增等[10]給出了一種考慮套管幾何尺寸的抗擠毀計(jì)算公式，但該公式依據(jù)彈性力學(xué)推導(dǎo)無法計(jì)算后屈曲行為；孫坤忠[11]引用了K-T公式，該公式同樣不能計(jì)算套管的后屈曲行為；曾義金等[12]研究了鹽膏層蠕變產(chǎn)生的擠壓應(yīng)力隨時(shí)間的變化，沒有給出套管抗擠毀性能的評(píng)價(jià)方法；申照熙等[13]研究了影響套管抗擠毀性能的因素，并給出了計(jì)算公式，但沒有考慮到塑形套管擠毀過程的后屈曲行為；覃成錦等[14]分析了不同應(yīng)力狀態(tài)對(duì)套管抗擠毀性能的影響，仇偉德等[15]深入分析了影響抗擠毀性能的各項(xiàng)因素，但沒有給出計(jì)算公式；張?jiān)旅舻萚16]介紹了TP130TT高抗擠毀套管的成功應(yīng)用，在使用過程中存在下入磨阻大問題。因此，準(zhǔn)確快速計(jì)算套管的抗擠毀值，可以加強(qiáng)對(duì)套管性能認(rèn)識(shí)，避免過度設(shè)計(jì)，降低材料消耗，所以如何有效評(píng)價(jià)套管的抗擠毀性能成為一個(gè)關(guān)鍵問題。當(dāng)外壓逐漸增大，套管發(fā)生屈服變形，隨著變形量的增加，套管抗擠毀能力提高(見圖1)，第一次變形稱為線性特征值屈曲或前屈曲；第二次變形稱為非線性屈曲或后屈曲[17- 18]。

圖1 載荷位移響應(yīng)示意圖

實(shí)際生產(chǎn)、運(yùn)輸儲(chǔ)存、使用過程不可避免地存在缺欠，或?yàn)椴牧先鼻坊驗(yàn)閹缀稳鼻罚@種缺欠按照API 5CT也是允許存在的，而API 5C3對(duì)此處理較為粗糙，通過線性及非線性分析方法可以較好的解決該問題。

1 套管的線性屈曲分析

套管的線性特征值屈曲，其線彈性屈曲載荷{P0}與位移{U0}關(guān)系：

[Ke]{U0}={P0}

(1)

式中:[Ke]—彈性剛度矩陣;

U0—由{P0}引起的位移;

{σ}—與{U0}對(duì)應(yīng)的應(yīng)力。

在任意狀態(tài)下的增量平衡方程：

{ΔP}=([Ke]+[KG(σ)]){ΔU}

(2)

式中:[KG(σ)]—初始應(yīng)力矩陣。

假設(shè)線彈性屈曲是一個(gè)外加載荷{P0}的線性函數(shù)，即:

{P}=λ{(lán)P0}

{U}=λ{(lán)U0}

{σ}=λ{(lán)σ0}

則，可得：

[KG(σ)]=λ[KG(σ0)]

(3)

因此，其增量平衡方程變?yōu)?

{ΔP}=([Ke]+λ[KG(σ0)]){ΔU}

(4)

當(dāng)發(fā)生線性屈曲時(shí){ΔP}=0，于是式(4)變?yōu)?

([Ke]+λ[KG(σ0)]){ΔU}={0}

(5)

那么必須有det{[Ke]+λ[KG(σ0)]}=0，由該方程解出關(guān)于λ的n階多項(xiàng)式，解出最小特征值λ，λ代表臨界抗擠毀載荷比例因子，最小特征值乘以初始載荷得到臨界載荷，即抗擠毀值。線性特征值分析對(duì)于解決線彈性變形問題可以提供較好的解決方案，但是對(duì)于屈曲及后屈曲問題計(jì)算結(jié)果偏大，為此需要進(jìn)行弧長法分析。實(shí)際分析當(dāng)中，對(duì)于徑壁比較大套管，使用線性特征值分析即可得到準(zhǔn)確解；對(duì)于徑壁比較小的套管需要采用兩種方法進(jìn)行計(jì)算分析。

2 弧長法計(jì)算套管的抗擠毀強(qiáng)度

2.1 弧長法計(jì)算理論分析

通過引用弧長因子(ΔL)，將載荷因子λ和位移增量{ΔU}相聯(lián)系，彌補(bǔ)了Newton-Raphson方法極值點(diǎn)發(fā)散的缺點(diǎn)，可以精確求解后屈曲載荷，求解示意圖如圖1所示。

弧長法增量平衡迭代方程：

(6)

(7)

(8)

(9)

(10)

曲面控制弧長法約束方程可表述為：

(11)

令α=1，當(dāng)β=1時(shí)，為球面弧長法；當(dāng)β=0時(shí)，為柱面弧長法。當(dāng)β=當(dāng)前剛度數(shù)值時(shí)，為橢球面參數(shù)法。

下面以柱面弧長法為例，進(jìn)行闡述說明，其約束方程為

(12)

第i+1載荷步的位移增量關(guān)系為：

(13)

將式(10)帶入式(13)，得:

(14)

將式(14)帶入式(12)，得:

(15)

(16)

(17)

(18)

2.2 模型及其參數(shù)

選取油氣井常用的典型規(guī)格套管進(jìn)行建模計(jì)算,套管規(guī)格及材料參數(shù)見表1，模型參數(shù)長度2 000 mm。模型的分析單位類型為C3D8R,即8結(jié)點(diǎn)線性六面體單元，網(wǎng)格劃分在壁厚方向上分三層，對(duì)管體兩端進(jìn)行全約束，如圖2所示。

圖2 約束載荷施加情況

2.3 特征值計(jì)算

首先采用特征值分析方法對(duì)模型進(jìn)行線性分析，加載方式為施加均布外壓載荷，分析采用六階特征值方法，取第一階結(jié)果作為后續(xù)分析的模型如表1所示。一階結(jié)果如表2所示。

表1 套管規(guī)格及材料參數(shù)

表2 線性屈曲計(jì)算結(jié)果

2.4 非線性屈曲分析

取線性分析結(jié)果的一階模態(tài)進(jìn)行后屈曲分析，考慮壁厚的5%作為缺欠進(jìn)行缺欠敏感性分析。計(jì)算結(jié)果如表3所示，各個(gè)不同規(guī)格的套管計(jì)算結(jié)果與實(shí)物試驗(yàn)擠毀值比較可見，計(jì)算誤差在3%左右。通過與表2的線性屈曲計(jì)算結(jié)果相對(duì)比，特征值計(jì)算結(jié)果對(duì)于彈性擠毀套管計(jì)算較為準(zhǔn)確，但是對(duì)于塑形擠毀套管結(jié)果明顯偏大，說明線性分析存在局限性。雖然線性分析計(jì)算的結(jié)果偏大，但其可作為非線性分析的基礎(chǔ)模態(tài)，使非線性分析計(jì)算更為準(zhǔn)確。高度接近實(shí)物試驗(yàn)結(jié)果說明弧長法計(jì)算套管抗擠強(qiáng)度具有可行性。

表3 非線性分析計(jì)算結(jié)果

3 材料缺欠和彎曲度對(duì)抗擠性能的影響

3.1 彎曲度對(duì)抗擠毀性能的影響

以?139.7 mm×12.70 mm進(jìn)行建模。管體兩端全約束，管體中間通過剛體施加位移，模型及推彎過程如圖3所示。

圖3 推彎狀態(tài)

3.2 材料缺欠對(duì)抗擠性能的影響

套管自身存在的缺欠一定程度上影響套管的抗擠毀性能，為此對(duì)模型施加不同的缺欠模態(tài)，分別為壁厚的5%、10%、15%，初始載荷保持不變，計(jì)算結(jié)果如表4，不同規(guī)格的套管隨著缺欠模態(tài)的增加，lpf載荷比例因子逐漸減小，即套管抗擠毀性能逐漸降低，表明缺欠模態(tài)的增加降低套管的抗擠毀性能。

表4 有限元計(jì)算峰值因子

套管串在井下隨井眼軌跡增、降斜作業(yè)存在彎曲狀態(tài)，建模時(shí)通過在套管中間部位施加不同的位移，使其產(chǎn)生不同的彎曲度，分別為8°/30 m、16°/30 m、25°/30 m、41°/30 m。采用弧長法計(jì)算套管的抗擠毀性能，彎曲應(yīng)力疊加局部應(yīng)力集中導(dǎo)致套管抗擠毀性能降低，8°/30 m的彎曲度導(dǎo)致抗擠毀性能降低5%左右，到41°/30 m的抗擠毀性能降低35%，如圖4所示。通過弧長法可以有效表征套管在彎曲度存在的條件下抗擠毀強(qiáng)度損失情況,如圖5所示。

圖4 彎曲度與擠毀值的關(guān)系

圖5 擠毀狀態(tài)

4 實(shí)物試驗(yàn)與有限元計(jì)算對(duì)比

該分析方法可以準(zhǔn)確的計(jì)算鋼管的后屈曲行為及其結(jié)果，對(duì)比圖6～圖9發(fā)現(xiàn)鋼管壓潰的FEA模擬形貌和實(shí)物試驗(yàn)形貌一致，具有高度吻合性，即彈性擠毀套管發(fā)生彈性失穩(wěn)變形，變形部位呈現(xiàn)微橢圓狀態(tài)；塑形擠毀套管發(fā)生后屈曲變形，變形部位呈現(xiàn)擠毀形貌。結(jié)合數(shù)值比較和形貌比較，說明該方法可用于套管的計(jì)算分析。

圖6 FEA模擬失效形貌圖及實(shí)物擠毀圖(?139 mm×12.70 mm)

圖7 FEA模擬失效形貌圖及實(shí)物擠毀圖(?177 mm×11.51 mm)

圖8 FEA模擬失效形貌圖及實(shí)物擠毀圖(?244 mm×11.99 mm)

圖9 FEA模擬失效形貌圖及實(shí)物擠毀圖(?399 mm×13.06 mm)

結(jié)合表3數(shù)據(jù)說明，對(duì)于大尺寸套管以及處于彈性擠毀范圍的中度尺寸的套管可采用線性特征值分析直接進(jìn)行計(jì)算，其計(jì)算結(jié)果與Riks弧長法計(jì)算結(jié)果很接近，可以忽略兩者的差別。

5 結(jié)論

(1)特征值和弧長法計(jì)算可以有效計(jì)算套管的抗擠毀值。對(duì)于彈性擠毀套管，可直接用線性屈曲分析計(jì)算;對(duì)于塑形擠毀套管需要在線性屈曲分析的基礎(chǔ)上采用弧長法進(jìn)行計(jì)算。

(2)通過有限元計(jì)算與實(shí)物試驗(yàn)對(duì)比分析，單軸載荷條件下，套管擠毀值與實(shí)物試驗(yàn)結(jié)果一致，且FEA模擬失效形貌與實(shí)物試驗(yàn)失效形貌一致。

(3)隨著缺欠模態(tài)的增加，lpf載荷比例因子逐漸減小，套管擠毀值逐漸降低;彎曲度顯著影響擠毀值，隨著彎曲度增大，套管抗擠毀能力逐漸降低。鉆井過程中應(yīng)加強(qiáng)井眼軌跡控制，防止局部出現(xiàn)狗腿度過大現(xiàn)象。