壓裂管柱力學(xué)分析及其應(yīng)用

王尊策，張 辰，溫后珍，陳 思，柯 林

（東北石油大學(xué)機(jī)械科學(xué)與工程學(xué)院，黑龍江大慶 163318）①

壓裂管柱力學(xué)分析及其應(yīng)用

王尊策，張 辰，溫后珍，陳 思，柯 林

（東北石油大學(xué)機(jī)械科學(xué)與工程學(xué)院，黑龍江大慶 163318）①

根據(jù)壓裂管柱的壓裂施工過(guò)程、管柱結(jié)構(gòu)和井下工具特性，建立了壓裂管柱的綜合力學(xué)模型，研究了在施工過(guò)程中能夠影響管柱位移的4種基本效應(yīng)。運(yùn)用相應(yīng)公式及算法編制了壓裂管柱力學(xué)分析軟件，并結(jié)合大慶油田壓裂管柱實(shí)例進(jìn)行了力學(xué)分析計(jì)算。經(jīng)現(xiàn)場(chǎng)測(cè)試，通過(guò)軟件計(jì)算結(jié)果與現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)進(jìn)行比較，誤差在合理范圍內(nèi)，滿足工程需要。軟件所得出的計(jì)算結(jié)果對(duì)于現(xiàn)場(chǎng)壓裂作業(yè)具有一定的指導(dǎo)意義。

壓裂；管柱；軟件；力學(xué)分析；力學(xué)模型

壓裂技術(shù)是一項(xiàng)重要的油田增產(chǎn)技術(shù)。在壓裂施工過(guò)程中，壓裂工藝管柱在油管內(nèi)外壓力、彎矩、軸向力、轉(zhuǎn)矩、溫度等相關(guān)因素的作用下，形成較為復(fù)雜的應(yīng)力和應(yīng)變，有時(shí)甚至?xí)l(fā)生導(dǎo)致管柱失效的事故。例如，管柱的屈服破壞、斷脫破損或發(fā)生永久性的螺旋屈曲，導(dǎo)致管柱無(wú)法正常工作，造成很大的經(jīng)濟(jì)損失。因此，影響管柱壓裂施工是否成功的一項(xiàng)重要關(guān)鍵因素便是在施工過(guò)程中管柱的受力情況［1-2］。

本文建立了壓裂工藝管柱的力學(xué)綜合模型，運(yùn)用VS 2008和Matlab 2010等軟件編制相應(yīng)的管柱力學(xué)分析軟件，對(duì)壓裂管柱在不同工況下進(jìn)行受力和形變分析，找出影響壓裂管柱形變和位移的各種因素，從而對(duì)壓裂管柱的結(jié)構(gòu)和參數(shù)等進(jìn)行優(yōu)化設(shè)計(jì)，對(duì)提高壓裂施工作業(yè)的成功率、保證壓裂管柱正常施工具有重要的意義。

1 壓裂管柱的綜合力學(xué)模型

本文以大慶地區(qū)州扶49-平45井（如圖1）為例，建立深井壓裂管柱綜合力學(xué)模型。

圖1 州扶49-平45井壓裂管柱結(jié)構(gòu)

根據(jù)壓裂管柱的實(shí)際工作狀態(tài)，選取整個(gè)壓裂管柱為研究對(duì)象。在變形前，假設(shè)壓裂管柱與井眼的軸線重合，套管內(nèi)壁為剛性，壓裂管柱及工具均為彈性變形體；變形后壓裂管柱沿井眼軸線與套管內(nèi)壁產(chǎn)生接觸，在接觸處存在接觸反力及摩擦阻力。依據(jù)平衡方程，建立壓裂管柱的綜合力學(xué)模型，得到管柱基本微分方程［3］。

等效軸力微分方程為式中：φ為井斜角；ψ為方位角；Fτe為等效軸向力；N為接觸正壓力；θ為偏轉(zhuǎn)角。

根據(jù)上述3個(gè)方程和具體邊界條件，可以確定相應(yīng)的3個(gè)未知函數(shù)θ、Fτe和N［4］。

2 影響管柱位移的4種效應(yīng)

在進(jìn)行壓裂作業(yè)中，管柱自重是壓裂管柱所受的最主要載荷；同時(shí)由于注入了壓裂劑，導(dǎo)致油管內(nèi)和油套環(huán)空內(nèi)的溫度和壓力發(fā)生了變化，因此產(chǎn)生了4種基本效應(yīng)，這4種基本效應(yīng)能夠影響封隔器管柱的受力和長(zhǎng)度的變化［5］。

2.1 活塞效應(yīng)

由于油管內(nèi)、外壓力的變化從而導(dǎo)致油管長(zhǎng)度的變化，這種現(xiàn)象稱為活塞效應(yīng)，如圖2所示。活塞力是油管內(nèi)、外壓力對(duì)油管的合力，計(jì)算公式為

F1＝（Ap-A′i）pi-（Ap-A0）p0

式中：p0為油套環(huán)空壓力；pi為油管內(nèi)壓力；A0、Ai分別為油管外徑和內(nèi)徑截面積；Ap為封隔器密封腔的橫截面積。

圖2 活塞效應(yīng)示意

2.2 溫度效應(yīng)

井內(nèi)溫度隨井的深度增加而升高，而且在井中注入壓裂液或其他物質(zhì)等時(shí)，管柱溫度會(huì)相應(yīng)地發(fā)生變化，由于這種溫度的變化，管柱將會(huì)相應(yīng)地伸長(zhǎng)或縮短，這種現(xiàn)象稱為溫度效應(yīng)。由于整個(gè)管柱的每一處都會(huì)有溫度效應(yīng)的存在，因此整個(gè)管柱的平均溫度為

由于管柱平均溫度的變化引起的力ΔF2和長(zhǎng)度ΔL的變化為

式中：tS為井口溫度；tB為井底溫度；β為材料的熱膨脹系數(shù)；L為油管長(zhǎng)度；犠為單位長(zhǎng)度油管質(zhì)量。

2.3 鼓脹效應(yīng)

當(dāng)油管內(nèi)部有壓力存在時(shí)，油管內(nèi)的壓力作用于油管壁內(nèi)，使油管的直徑增大，因此油管將縮短，這種現(xiàn)象叫做正鼓脹效應(yīng)；當(dāng)油套環(huán)空內(nèi)有壓力存在時(shí)，油套環(huán)空的壓力將作用于油管壁外，使油管的直徑減小，管柱將伸長(zhǎng)，這種現(xiàn)象叫反向鼓脹效應(yīng)。鼓脹力計(jì)算公式為

式中：Δpia為管柱內(nèi)平均壓力的變化；Δp0a為管柱外平均壓力的變化。

2.4 螺旋屈曲效應(yīng)

由于壓力不僅沿井眼軸線垂直作用在封隔器的膠筒上，還會(huì)水平作用于整個(gè)壓裂管柱的壁面上。如果緊靠封隔器處的活塞力大于管柱發(fā)生螺旋屈曲的臨界力時(shí)，則油管將會(huì)產(chǎn)生螺旋屈曲現(xiàn)象。

從油管底部到中和點(diǎn)的距離為

式中：F為壓縮力；犠為油管單位長(zhǎng)度的平均質(zhì)量。

3 算法設(shè)計(jì)

將管柱的總長(zhǎng)設(shè)為L(zhǎng)，并將其劃分為犿段，則產(chǎn)生犿＋1個(gè)節(jié)點(diǎn)。從井口到井底依次排列這些節(jié)點(diǎn)，節(jié)點(diǎn)的序號(hào)分別為1，2，3，…，犿，犿＋1。節(jié)點(diǎn)的選取包括壓裂管柱結(jié)構(gòu)的改變、材料特性的改變、工具類型的改變及相應(yīng)計(jì)算參數(shù)不連續(xù)的點(diǎn)。然后應(yīng)用差分方法計(jì)算各節(jié)點(diǎn)參數(shù)。

3.1 管柱單元內(nèi)外壓力

管柱單元的內(nèi)外壓的差分方程為

當(dāng)井口油壓、套壓已知時(shí)，由上述方程可以確定每個(gè)節(jié)點(diǎn)處的內(nèi)、外壓力，但應(yīng)注意開(kāi)關(guān)元件、封隔器等處壓力的不連續(xù)性。其中：公式中的下角標(biāo)i為管柱內(nèi)側(cè)參數(shù)；o為管柱外側(cè)參數(shù)；j-1、j、j＋1均為管柱的節(jié)點(diǎn)編號(hào)；下同。

3.2 軸力、接觸壓力、彎矩的計(jì)算

無(wú)論何種工況，由于管柱最底端所受的等效軸力已知，故可根據(jù)由上述基本微分方程所推導(dǎo)出的公式計(jì)算該點(diǎn)的fnj及βj，即

當(dāng)1≤βj＜1.469時(shí)，管柱發(fā)生正弦屈曲，即

3.3 位移計(jì)算

如果管柱上節(jié)點(diǎn)j處的位移已知，那么可以利用下述公式求得相應(yīng)j＋1、j-1處節(jié)點(diǎn)的位移，即

式中：qej為j單元管柱的單位長(zhǎng)度重力，N／m；fj為j單元處管柱與套管間的摩擦因數(shù)；ρij、ρoj分別為j單元處管柱內(nèi)、外流體的密度，kg／m3；△Sj為j單元處管柱的長(zhǎng)度，m；φj為j單元處井眼軸線的井斜角，rad；ψj為j單元處井眼軸線的方位角，rad；nj為j單元處管柱的強(qiáng)度安全系數(shù)；koj為j單元處井眼軸線的曲率，m-1；Fτj為j單元處的實(shí)際軸向力，N。

4 實(shí)例應(yīng)用

根據(jù)上述理論，在Windows平臺(tái)下，應(yīng)用VS 2008，Matlab 2010及Office Access等軟件編寫(xiě)壓裂管柱力學(xué)分析軟件，運(yùn)用差分法對(duì)壓裂管柱進(jìn)行力學(xué)計(jì)算。該軟件可計(jì)算管柱在不同工況下位移和載荷的分布情況。

以大慶地區(qū)州扶49-平45井為例，該井壓裂管柱采用?70.4 mm和?88.9 mm外加厚2種復(fù)合油管，壓裂井段2 470.0～2 510.5 m，最高施工井口壓力為70 MPa，施工排量為2～6 m3／min，采用K344-110型封隔器。利用該軟件輸入井眼軌跡數(shù)據(jù)（如圖3）及各項(xiàng)參數(shù)，最后針對(duì)該管柱進(jìn)行力學(xué)計(jì)算。

圖3 州扶49-平45井井眼軌跡

根據(jù)軟件計(jì)算結(jié)果可知，管柱發(fā)生伸縮變形，最大變形量為5.85 m，發(fā)生在第2個(gè)封隔器位置處。管柱在壓裂作業(yè)中，管柱最小安全系數(shù)為1.1，最大拉力為291 k N、最大應(yīng)力為365 MPa，小于許用值，滿足強(qiáng)度條件。

為了驗(yàn)證軟件計(jì)算的準(zhǔn)確性，在現(xiàn)場(chǎng)施工中，使用相應(yīng)的測(cè)試工具測(cè)量了壓裂管柱在不同位置處的軸向力，并與軟件計(jì)算結(jié)果進(jìn)行比較，如表1～2。通過(guò)比較得知，理論計(jì)算載荷和實(shí)際測(cè)試載荷基本一致，該軟件計(jì)算結(jié)果完全滿足壓裂施工的要求。由此可見(jiàn)，本軟件的計(jì)算結(jié)果可為壓裂管柱設(shè)計(jì)和施工提供理論依據(jù)。

表1 軸向力計(jì)算結(jié)果與實(shí)際測(cè)試結(jié)果對(duì)比

表2 工況參數(shù)及封隔器位移數(shù)據(jù)

5 結(jié)論

1） 該軟件與其他軟件相比，可輸入三維的井眼軌跡參數(shù)，井眼軌跡圖可進(jìn)行3D旋轉(zhuǎn)，更直觀地揭示了管柱作業(yè)時(shí)的危險(xiǎn)點(diǎn)或危險(xiǎn)段。

2） 可根據(jù)現(xiàn)場(chǎng)不同需求進(jìn)行壓裂工藝管柱的不同配置。軟件能自動(dòng)識(shí)別管柱結(jié)構(gòu)與工況參數(shù)，實(shí)現(xiàn)邊界條件的自動(dòng)載入。

3） 該軟件可針對(duì)壓裂管柱進(jìn)行不同工況下的力學(xué)分析，預(yù)測(cè)每個(gè)封隔器在不同工況下產(chǎn)生的位移，并能計(jì)算出管柱各個(gè)位置所受到的不同載荷。與現(xiàn)場(chǎng)測(cè)試得出的結(jié)果相比較，計(jì)算結(jié)果與測(cè)試結(jié)果的誤差在合理范圍內(nèi)，能夠滿足工程需要。

［1］ 杜現(xiàn)飛，王海文，王帥，等.深井壓裂井下管柱力學(xué)分析及其應(yīng)用［J］.石油礦場(chǎng)機(jī)械，2008，37（8）：28-33.

［2］ 李子豐，孫虎，蘇金柱，等.壓裂管柱力學(xué)分析理論與應(yīng)用［J］.應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào)，2012，（5）：846-861.

［3］ Lubinski A，Althouse W S，Logan J L.Helical Buckling of Tubing Sealed in Packers［J］.JPT，1962，14（3）：655-670.

［4］ Hammerlindl，D J Movement.Forces and Stresses Associated with Combination Tubing Strings Sealed in Packers［J］.JPT，1977，29（1）：195-208.

［5］ 黃云.注水工藝管柱力學(xué)行為研究［D］.南充：西南石油學(xué)院，2005.

Mechanical Analysis and Application of Fracturing String

According to the process of fracturing program and the structure of fracturing strings and characteristics of the dowhole tools，a comprehensive mechanical model of fracturing strings is established.The four basic effects which influenced displacement of the fracturing strings were studied.An appropriate formulas and algorithms were used to establish a set of fracturing strings mechanical analysis software.Daqing oilfield fracturing strings was used as an instance to process the mechanics analysis and calculation with the software.According to the results from the field test，the error between the data from software and the field test meets the needs of the error within a reasonable range.Therefore，the results obtained from this software have important implications for field fracturing operations.

fracturing；strings；software；mechanical analysis；mechanical model.

TE934.101

A

10.3969／j.issn.1001-3842.2014.09.008

1001-3482（2014）09-0029-04

2014-03-11

石油裝備制造創(chuàng)新服務(wù)平臺(tái)建設(shè)（2012BAH28F03）

王尊策（1962-），男，黑龍江同江人，教授，博士生導(dǎo)師，主要從事井下工具、流體機(jī)械工作理論及技術(shù)領(lǐng)域的研究，E-mail：wangzc＠nepu.edu.cn。