鹽穴儲油（氣）庫水溶造腔管柱動力特性初探

李銀平，楊春和,，屈丹安，楊長來，施錫林

（1. 中國科學院武漢巖土力學研究所 巖土力學與工程國家重點實驗室，武漢 430071；2. 重慶大學 西南資源開發及環境災害控制工程教育部重點實驗室，重慶 400044；3. 中國石油化工股份有限公司天然氣分公司，北京 100120）

1 引 言

深部鹽巖構造是石油、天然氣及高放核廢料地下存儲的理想場所[1]。目前建造鹽巖儲庫普遍采用單井油墊對流法水溶開采[2]。如圖1所示，地面的淡水或淡鹵水經造腔內管（或者造腔內管和造腔外管之間環空）注入腔體，腔內鹵水從造腔內管與造腔外管之間環空（或者造腔內管）流出到地面，流動的非飽和鹵水不斷溶蝕腔壁鹽巖，腔體的體積逐步擴大；造腔外管和生產套管之間充滿柴油或其他液體，下行形成油墊層隔開淡鹵水和腔頂鹽巖，達到抑制上溶并促使腔體橫向溶漓，或者上行撤除油墊層，讓淡鹵水充分接觸腔頂鹽巖，腔體向上發展，從而達到建造較為理想的儲庫腔體形狀的目的。

圖1 對流法水溶造腔示意圖Fig.1 Sketch diagram of the solution mining process

建腔目的層一般埋深在500～2000 m之間，造腔管柱（內管和外管）從地面井口自由懸伸至地下溶腔中，考慮到修井和腔體測試需要，造腔內管和外管內壁之間以及造腔外管和生產套管內壁之間一般不設任何支撐。

造腔工程實際中，由于造腔內管破損導致的事故時有發生。如圖2所示，某儲庫造腔內管發生嚴重彎曲變形。修井時需要提升至地面，由于彎曲致使其在有限圓筒空間上升受阻，強行提升就會造成外管及套管的破損，工程上有時不得不割斷彎曲嚴重的造腔內管，永久棄之于溶腔之中。

圖2 造腔現場提出的S型彎曲造腔內管Fig.2 The crooked inner casing in-situ

造腔外管破損現象較少，但在某埋深較大儲庫建設過程中，發現下端最后1根φ177.8 mm造腔外管在距母扣1/3處有明顯斷裂現象（見圖3）。

另外，造腔管柱接箍損壞現象也時有發生（見圖4），這會導致出水口深度發生改變和腔體形狀控制失敗。

造腔管柱的破損，不僅會提高造腔成本，延長工期，嚴重的還會造成工程事故，造腔失敗。但國內外關于水溶造腔管柱破損的研究并不多見。近年來施錫林等[3－4]針對我國鹽礦一般含有眾多難溶夾層，如硬石膏層、鈣芒硝層和泥巖層等特點，系統研究了造腔過程中夾層垮塌導致造腔內管被砸彎、砸壞以及套管被卡等事故的機制。認為自下而上水溶造腔過程中，腔頂以及腔體中懸空的難溶夾層垮塌，下行不規則巖塊在溶漓流程中曲線運動，側向撞擊造腔內管，導致管柱的彎曲變形及托箍等（見圖5）。文獻[5]還結合現場研究，提出了相應的夾層垮塌控制技術。

圖3 造腔現場提出的斷裂造腔外管Fig.3 The fractured casing in-situ

圖4 套管接箍損壞Fig.4 Casing coupling damage

圖5 難溶夾層的垮塌及其造成的管柱損壞Fig.5 Collapse of interbeds and tubing string damage

本文基于對造腔現場勘察以及對工況的分析認為，管柱破損除了與造腔過程中管柱所受外力撞擊和擠壓等有關外，更重要的是與其自身運動的動力特性緊密相關，管柱在受限空間的自激振動和動力失穩是導致管柱發生嚴重彎曲或者破損的主要原因。本文對管柱動力特性開展了初步研究，考慮液-固耦聯作用分析造腔管柱的臨界流速，還對管柱水擊等現象開展的初步分析，最后對下一步結合實際工況更深入研究管柱動力特性提出了建議。

2 輸流管道液-固耦聯振動及臨界流速

由于輸流管道被廣泛應用于石油、化工、能源、宇航及海洋工程等領域，因此，一直受到學術界的重視[6]。由于管道與流體的耦合作用，系統的運動方程中含有陀螺力項（又稱哥氏加速度項），這使得輸流管道的振動現象千變萬化，所蘊含的振動機制也極其復雜。被關注較多的是兩端支撐管及懸壁管，此兩類管道在力學上有較大的不同，懸壁管屬于非保守系統，其管道-流體耦合系統的能量隨時間而發生變化，相對較為復雜[7]。

2.1 懸臂輸流管振動方程

圖1中，若不考慮各管道之間及與套管的相互干擾，即振動不受空間限制，則造腔內管和外管都可以看成是懸臂輸流管道（見圖6）。這里先以造腔內管為例開展分析，當處于正循環工況時，淡水從地面注入，通過內管輸送射入溶腔，管道為排液管（cantilevered pipe discharging fluid）；反循環時，鹵水從溶腔進入內管，并排出地面，造腔內管稱為吸流管（cantilevered pipe aspirating fluid）。

圖6 懸臂輸流管道Fig.6 Fluid conveying cantilevered pipe

當管內無液體時，若不計重力、管道剪切變形、截面轉動慣量的影響，管道沿橫向彎曲振動的微分方程可寫為[8]

式中：EI為管道的彎曲剛度；m為管道單位長度的質量；T為下端軸向力（拉力為正，暫不計重力）；管道的橫向變形 w=w(x, t)，管道外激發力q=q(x, t)；x和t分別為軸向坐標和時間。

當管內有液體時，管道的振動將會有很大不同，假設流體為線流、不可壓縮且無黏性，液體單位長度內作用在管道上的慣性力為

式中：ρ為液體密度；A為液體橫截面積；u為液體流速。

考慮液體壓力P和上述慣性力F后，管道彎曲振動微分方程為

對于懸臂管道，其邊界條件為

式中：l為管道總長度。

2.2 懸臂輸流管道臨界流速分析

令式（4）中q(x, t)＝0即可得到管道自由振動運動方程為

運算后寫成矩陣的形式為

其中ω0j為充液管(u=0)的第j階振型固有頻率。

通過以上離散，問題化解為一個二階常微分方程組，即式（8）。它的解有如下形式：

將式（9）代入式（8）得

于是頻率方程為

式（11）的求解屬于一個代數復特征值問題，可求出臨界流速ucr。一般情況下，s為復數，其虛部就是輸液管的固有頻率，若s具有正實部，這時振動幅值將隨時間增長而迅速最大，狀態變成了不穩定，這種現象被稱為顫振（flutter）。

3 造腔管柱臨界流速計算及討論

3.1 造腔管柱臨界流速計算

對于造腔管柱的臨界流速，這里先考慮最簡化的情形來進行分析。造腔內管先不考慮造腔外管的空間限制，簡化為從井口直至溶腔內的懸臂輸流管；而對于造腔外管，則先不考慮外部生產套管及裸眼段的空間限制，也不考慮其內部造腔內管的影響，當然參與振動的液體只是內外管環空中液體。這樣就都簡化為如圖6所示的懸壁輸流管。

對于懸臂管，當u＝0時，其前兩階固有頻率分別為

對應的前兩階振型函數分別為

將振型函數式（13）代入頻率方程式（6），可求得懸臂輸流管臨界流速為

水溶造腔工程中常用的114.3 mm造腔內管及177.8 mm造腔外管的橫斷面如圖7所示，空白部分為所考慮液體的橫斷面。材料的彈性模量 E＝200 GPa。由式（14）可得到

造腔內管臨界流速（m/s）為

造腔外管臨界流速（m/s）為

圖7 造腔管柱橫截面Fig.7 Cross-sections of tubing string

3.2 造腔管柱臨界流速討論

水溶造腔施工過程中，一般注水量為 40～100 m3/h，因此，造腔內管中流速為1.37～3.43 m/s；造腔外管中環空液體流速為 1.08～2.71 m/s。由式（15）和（16）可以看出，造腔內管臨界流速較低，而其內液體流速較高，因此，造腔內管一般會先于外管達到臨界流速而失穩。由于一般溶腔埋深都大于 500 m，由式（15）可知，造腔內管臨界流速ucr<0.83 m/s，相應的臨界流量Qcr<24.19 m3/h，這也就是說對于長度為500 m的造腔內管，當流量達到24.19 m3/h時就失穩了。

以上計算表明，若不考慮造腔外管空間約束，水溶造腔工程中造腔內管一般會失穩甚至發生破壞。

實際上，造腔管柱之間以及外部套管的相互約束是必須考慮的因素。由圖1可知，造腔內管一般比外管長（工程上一般長為30～60 m），懸伸于外管之外的部分是自由的，而嵌套于外管內的部分要受到外管的約束，如圖8所示。很顯然，這種空間約束將限制管道的橫向振動，無疑會提高其臨界流速，但內外管及與生產套管之間的碰撞、擠壓會帶來管柱的磨損和破壞。

圖8 受約束的懸臂輸流管道Fig.8 Cantilevered fluid conveying pipe with space constraint

因此，考慮造腔管柱之間相互干擾，受限空間內外管柱振動耦合分析，這將是一個十分復雜的問題，是下一步研究的重點之一。高德利[9]針對受井眼約束的受壓管柱，假設管柱和井眼均為圓形且二者保持連續接觸，采用細長梁理論開展了管柱屈曲分析。但對于主要受拉力的造腔管柱，如何給出合理的管柱之間及與井眼的接觸條件，需要進一步研究分析。

4 造腔管柱臨界流速計算力學模型初步探討

上一節討論了外部約束（受限空間）對輸流管動力特性影響，分析表明，管柱抗彎剛度及管柱懸伸長度是決定管柱動力特性的兩個關鍵因素，其實對于實際工程問題還有諸多因素需要進一步探討，由于問題的復雜性，這里僅對各因素的影響做初步探討和分析。

4.1 非均勻軸力及液體壓力的影響

造腔管柱一般是垂直懸伸至溶腔中，軸向要受到重力、浮力以及管內外流體摩擦力作用，沿軸向的軸力應該是近似線性非均勻分布，軸力對輸流管動力特性有明顯的影響，拉力將會提高管道的振動頻率，進而將提高其臨界流速。此外，溶腔內的液壓一般較高，1000 m埋深溶腔，其液體壓力一般為10 MPa以上，而且液體壓力沿管柱軸向近似呈線性變化，液體壓力對管柱動力特性的影響有多大也需要進一步研究。當然附著在管道外壁的液體對于管道的橫向振動也會有一定影響。當然，決定管柱動力特性的關鍵因素還是管柱抗彎剛度和懸伸長度。

4.2 正反循環造腔工藝的影響

水溶造腔有正循環和反循環兩種方式，如圖 1所示，從內管注入淡水，從內外管環空排出鹵水為正循環方式，反之為反循環方式。正循環方式時，內管為排液管，外管則為吸流管。已有研究表明，排液管和吸流管的動力特性是不相同的，懸臂輸流管道屬于非保守系統, 這意味著管道-流體耦合系統的能量隨時間而發生變化，能量的增減取決于管內流體的輸送方向和流速大小。關于吸流管，學術界還發生了一次著名的爭論[10]：吸流管道在無窮小流速下是否會發生失穩以及管道自由端口處的流-固耦合機制如何。針對這一疑問, 學術界相繼提出一些理論（或假說）和闡釋[11－12]。事實上，對這一問題的探討不僅是學術上的需要, 而且具有重要的工程實際意義。

4.3 水擊和管道橫向振動耦合影響分析

水等液體在壓力管道中流動時，當遇到閥門突然關閉（或開啟）、水泵突然停機（或啟動）時，水的流速會發生突然的急劇變化，由于流體的慣性和壓縮性，可引起管道中流動的液體壓力發生反復的、急劇的周期性變化，此時壓力將大大超過正常工作壓力，并伴隨著管壁的擴張和收縮，發出強烈的振動和噪音，有如管道受到錘擊的聲音，這種現象稱為水擊或水錘[13]。水擊現象發生時，會產生一個很大的壓力躍升，使管壁材料及管道上的設備及附件承受很大的壓力，而可能產生嚴重的變形以致破壞。如果此時管道系統存在缺陷，則有可能對管道系統或設備造成破壞，導致事故的發生。

對于造腔管柱，水擊現象導致的壓力脈動可能導致套管接箍的損壞。可以想象，當軸向的水擊現象和橫向的彎曲振動耦合時，本來就承受較大軸向拉力套管接箍就極有可能發生松動，甚至脫扣破壞（見圖 4）帶來工程事故。因此，水擊和管道橫向振動耦合影響分析也將是下一步需要研究的課題。

4.4 造腔管柱中間支撐影響分析及可行性

從以上分析可以看出，造腔管柱由地面懸伸至溶腔，長度500 m以上，盡管水平空間受限，但工程上在內外管之間以及外管和套管之間普遍不設置水平支撐，這導致其臨界流速較低。結合工程，探討是否可能在不影響修井過程上提管柱和溶腔聲納測試的前提下，在外管和生產套管之間以及內外管之間設置易于拆卸的水平支撐，并開展支撐對輸液管柱振動特性影響分析，對支撐設計時數量及位置進行分析，并結合工程實踐開展應用研究。

5 結 論

針對鹽穴儲氣庫單井水溶造腔過程中造腔管柱動力特性開展了初步探討，主要結論如下：

（1）造腔工程中彎曲和管柱破損除了與造腔過程中管柱所受外力撞擊和擠壓等有關外，管柱在受限空間的自激振動和動力失穩是導致管柱發生嚴重彎曲或者破損的主要原因。

（2）分析了無空間約束條件下內外管柱的振動特性，將其簡化為懸臂輸流管，得到了管柱臨界流速的計算式，算例表明，無約束條件下造腔管柱一般會由于液-固耦聯振動而發生動力失穩，而且造腔內管一般會先于造腔外管發生失穩破壞。

（3）對影響造腔管柱臨界流速的其他因素進行了探討，初步分析了管道軸力、液體壓力、造腔循環方式、管道水平支撐，以及軸向水擊和橫向振動耦合的影響，為進一步結合工程實際開展空間約束條件下造腔管柱動力特性全面分析提出了建議。

[1] 楊春和, 李銀平, 陳鋒. 層狀巖鹽力學理論與工程[M].北京: 科學出版社, 2009.

[2] 王清明. 鹽類礦床水溶開采[M]. 北京: 化學工業出版社, 2003.

[3] 施錫林, 李銀平, 楊春和, 等. 鹵水浸泡對泥質夾層抗拉強度影響的試驗研究[J]. 巖石力學與工程學報, 2009,28(11): 2301－2308.SHI Xi-lin, LI Yin-ping, YANG Chun-he, et al.Experimental study of the tensile strength weakening of mudstone interlayer samples soaked in brine[J]. Chinese Journal of Rock Mechanics and Engineering, 2009,28(11): 2301－2308.

[4] 施錫林, 李銀平, 楊春和, 等. 鹽穴儲氣庫水溶造腔夾層垮塌力學機制研究[J]. 巖土力學, 2009, 30(12): 3615－3620, 3626.SHI Xi-lin, LI Yin-ping, YANG Chun-he, et al. Research on mechanical mechanism of interlayer collapse in solution mining for salt cavern gas storage[J]. Rock and Soil Mechanics, 2009, 30(12): 3615－3620, 3626.

[5] 施錫林, 李銀平, 楊春和, 等. 多夾層鹽礦油氣儲庫水溶建腔夾層垮塌控制技術[J]. 巖土工程學報, 2011,33(12): 1957－1963.SHI Xi-lin, LI Yin-ping, YANG Chun-he, et al.Technology of controlling interbeds collapse in solution mining for oil/gas storage in multi-interbedded salt formation[J]. Chinese Journal of Geotechnical Engineering, 2011, 33(12): 1957－1963.

[6] 王琳, 匡友弟, 黃玉盈, 等. 輸液管振動與穩定性研究的新進展: 從宏觀尺度到微納米尺度[J]. 固體力學學報, 2010, 31(5): 481－492.WANG Lin, KUANG You-di, HUANG Yu-ying, et al.Recent development on vibration and stability of pipes conveying fluid: From macro-scale to micro-and-nanoscales[J]. Chinese Journal of Solid Mechanics, 2010,31(5): 481－492.

[7] 黃玉盈, 鄒時智, 錢勤, 等. 輸液管的非線性振動、分叉與混沌——現狀與展望[J]. 力學進展, 1998, 28(1):30－42.HUANG Yu-ying, ZOU Shi-zhi, QIAN Qin, et al.Advances and trends of nonlinear dynamics of pipes conveying fluid[J]. Advance in Mechanics, 1998, 28(1):30－42.

[8] 黃玉盈. 結構振動分析基礎[M]. 武漢: 華中工學院出版社, 1988.

[9] 高德利. 油氣井管柱力學與工程[M]. 北京: 中國石油大學出版社, 2006.

[10] 倪樵, 王琳, 黃玉盈. 吸流管道動力學模型的研究現狀與展望[J] . 應用力學學報, 2008, 25(3): 450－454.NI Qiao, WANG Lin, HUANG Yu-ying. Advances and trends of a dynamical model: Pipes aspirating fluid[J].Chinese Journal of Applied Mechanics, 2008, 25(3):450－454.

[11] PAIDOUSSIS M P, LI G X. Pipes conveying fluid: A model dynamical problem[J]. Journal of Fluids and Structures, 1993, 8: 137－204.

[12] PAIDOUSSIS M P. The canonical problem of the fluid conveying pipe and radiation of the knowledge gained to other dynamics problems across applied mechanics[J].Journal of Sound and Vibration, 2008, 310(3): 462－492.

[13] 陳鑫, 魯傳敬, 李長俊, 等. 輸油管道系統水擊分析與計算[J]. 油氣儲運, 2002, 21(12): 23－26.CHEN Xin, LU Chuan-jing, LI Chang-jun, et al. Water hammer analysis and its application to oil pipeline system[J]. Oil/Gas Storage and Transportation, 2002,21(12): 23－26.