張立松，張士巖，閆相禎，曹宇光

（1.山東省油氣儲運安全省級重點實驗室青島市環(huán)海油氣儲運技術重點實驗室，山東 青島 266580；2．中國石油大學（華東）儲運與建筑工程學院，山東 青島 266580）

土體力學性質對輸氣管道穿越地震斷層的影響分析

張立松1，2，張士巖2，閆相禎2，曹宇光2

（1.山東省油氣儲運安全省級重點實驗室青島市環(huán)海油氣儲運技術重點實驗室，山東 青島 266580；2．中國石油大學（華東）儲運與建筑工程學院，山東 青島 266580）

為獲取土體力學性質對輸氣管道應變行為的影響，給出等效邊界彈簧設置在管道彈塑性變形段的本構關系，并據(jù)此建立考慮等效邊界彈簧的輸氣管道穿越地震斷層殼-接觸有限元模型。利用該有限元模型，分析地震斷層作用下土體力學性質對輸氣管道受拉側和受壓側最大軸向應變的影響規(guī)律。結果表明：在較小的斷層位移下，土體彈性模量和內摩擦角對最大軸向應變有較大影響；隨著斷層位移的增大，土體彈性模量和內摩擦角對最大軸向應變的影響降低。無論較小斷層位移還是較大斷層位移，土體內聚力對輸氣管道最大軸向應變的影響均較小。

輸氣管道；地震斷層；土體力學性質；等效邊界彈簧；殼-接觸模型

0 引言

國內西氣東輸、川氣東送、新粵浙等大型輸氣管線項目的不斷建設，使得X70、X80鋼級的大口徑輸氣管道大量應用。由于長距離輸氣管道不可避免地穿越地質災害區(qū)，特別是地震斷層區(qū)，使得地震斷層作用下的輸氣管道安全評價成為工程設計的關鍵問題[1-3]。

目前，國內外學者在這方面開展了大量研究，并取得了豐富的成果，國內學者趙林等[4]引入土彈簧模型，建立了輸氣管道穿越地震斷層的有限元模型。劉愛文等[5]引入等效邊界彈簧，建立了輸氣管道穿越地震斷層的殼-彈簧有限元模型。徐龍軍等[6]基于ABAQUS軟件研究了地震斷層作用下的管道力學行為，并獲取了管線的位移、應力、應變結果。李鴻晶[7-8]重點分析了油氣管道穿越地震斷層的有限元建模問題，并認為有限元模型的準確性是解決輸氣管道穿越地震斷層的關鍵點。趙雷等[9]利用ANSYS軟件，分別討論了輸氣管道壓縮屈曲失效和Euler梁失穩(wěn)失效的臨界條件。國外學者NEWMARK[10]、KENNEDY[11]、WANG等[12]對于埋地管道穿越地震斷層的前期研究上起到了重要作用。近幾年，Trifonov等[13]從解析模型入手，建立了考慮服役載荷的輸氣管道穿越斷層平面應力彈塑性模型。Vazouras等[14]基于ABAQUS軟件，討論了徑厚比、斷層位錯量、管道與斷層交角等影響因素對管道壓縮和拉伸失效的影響。

雖然國內外學者取得了豐富的成果，但少有學者專門開展地震斷層作用下土體力學性質（如土體彈性模量、內摩擦角及內聚力）對輸氣管道最大軸向應變的影響研究?；诖?，本文提出了考慮等效邊界彈簧的輸氣管道穿越地震斷層的殼-接觸有限元模型，并據(jù)此有限元模型討論了斷層土體力學性質對輸氣管道力學行為的影響規(guī)律。

1 考慮等效邊界彈簧的管道穿越斷層殼-接觸模型

1.1 力學模型

現(xiàn)行的殼-接觸有限元模型，普遍假設有限元模型長度為60倍管徑，且模型端部采用固定端約束。然而，由于較小的模型尺寸，管道在端部仍然存在軸向位移。為此，提出了考慮等效邊界彈簧的輸氣管道穿越地震斷層殼-接觸力學模型，如圖1所示。

1.2 等效邊界彈簧的本構關系

劉愛文等[5]最早將等效邊界彈簧引入輸氣管道穿越地震斷層的有限元分析中，然而并未給出等效邊界彈簧設置在管道彈塑性段的本構關系?；诖?，利用管土相互作用理論，給出了等效邊界彈簧設置在管線彈塑性段的本構關系式為

圖1 考慮等效邊界彈簧的殼-接觸有限元模型

式中：F——模型端部的管道拉力，N；

Δx——模型端部的管道變形，m；

x——以斷層面為原點的管道截面位置，m；

σs——屈服強度，MPa；

fs——管土靜摩擦力，N；

A——管道橫截面積，m2；

σo——斷層面處的管道應力，MPa；

E——彈性模量，MPa；

u0——土體屈服位移，m；

Es——塑性模量，MPa。

2 有限元模型

2.1 管殼單元

利用殼單元勢函數(shù)變?yōu)榈扔诹愕臈l件，獲取殼單元的剛度矩陣為

式中：[Db]——彎曲變形剛度矩陣；

[Ds]——剪切變形剛度矩陣；

[Bb]——彎曲變形應變矩陣；

[Bs]——剪切變形應變矩陣；

|J|——雅可比矩陣；

dξ、dη——積分變量。

2.2 土體單元

在本次分析中，采用Mohr-Coulomb準則模擬管土接觸面附近土體的彈塑性力學行為。根據(jù)塑性流動法則，總應變增量可表示為

式中：εe、εp——彈性和塑性應變；

D——彈性矩陣；

F——屈服函數(shù)；

H′——硬化系數(shù)。

3 結果分析

3.1 有限元模型及結果

新粵浙輸氣管線為國內新鋪設的重點管網(wǎng)，主干線設計為X80鋼，管徑1 219mm，壁厚26.4mm，設計壓力12MPa，彈性模量206GPa，泊松比0.3，屈服強度555MPa。該管線在陜西渭南穿越了一剪滑地震斷層，斷層傾角為90°，管道穿越斷層交角為42°，斷層位錯量1.5m，管道埋深1.4m，土體容重16.8kN/m3，土體內摩擦角22°，土體內聚力24kPa。

利用上述基礎參數(shù)，建立考慮等效邊界彈簧的輸氣管道穿越地震斷層殼-接觸有限元模型。其中，模型尺寸為73.14m×13.42m×6.10m，等效邊界彈簧的剛度系數(shù)為18.36×106N/m。分別施加0.5，1.0，1.5，2.0，2.5m的斷層位移，分析不同斷層位錯量下的輸氣管道彈塑性力學行為，并與實驗結果對比驗證有限元模型的正確性（見圖2）。

圖2 不同斷層位錯下管道應變數(shù)值結果與實驗結果對比示意圖

從圖中可知，利用本模型獲取的輸氣管道穿越地震斷層的應變數(shù)值結果與實驗結果吻合較好。對于管道受拉側的最大軸向應變，其數(shù)值結果和實驗結果最大相對誤差僅為7.52%；對于管道受拉側的最大軸向應變，其數(shù)值結果和實驗結果最大相對誤差僅為6.10%。

3.2 影響因素分析

3.2.1 土體彈性模量的影響

圖3給出了地震斷層作用下土體彈性模量對輸氣管道最大軸向拉、壓應變的影響。在斷層位移較小時，土體彈性模量對輸氣管道受拉側的最大軸向應變有較大影響。然而，隨著斷層位移的增加，土體彈性模量對輸氣管道受拉側最大軸向應變的影響降低。對于0.5m的斷層位移，當土體彈性模量從50MPa增加到500MPa時，輸氣管道受拉側的最大軸向應變從0.25%增大到0.78%。然而，當斷層位移增加到2.5m時，隨著土體彈性模量從50MPa變化到500MPa，受拉側的最大軸向應變僅從2.62%增加到2.74%。這說明在斷層位移較小時，土體彈性模型對輸氣管道受拉側最大軸向應變有較大影響。

圖3 土體彈性模量對輸氣管道最大軸向應變的影響

3.2.2 土體內摩擦角的影響

圖4土體內摩擦角對輸氣管道最大軸向應變的影響

圖4 給出了地震斷層作用下土體內摩擦角對輸氣管道最大軸向拉、壓應變的影響。由圖可知，在斷層位移為0.5m時，隨著土體內摩擦角的增大，輸氣管道受拉側最大軸向應變明顯增大。然而，當斷層位移為2.5m時，土體內摩擦角對輸氣管道受拉側最大軸向應變的影響降低。其原因主要為土體內摩擦角是控制土體屈服的關鍵參數(shù)。由于土體內摩擦角的增大，土體屈服變得困難。在較小的斷層位移作用下，這一現(xiàn)象尤為明顯，直接增大了輸氣管道受拉側的最大軸向應變。當斷層位移增大到2.5m時，這一作用降低，因此最大軸向應變隨土體內摩擦角的影響變小。

3.2.3 土體內聚力的影響

圖5給出了地震斷層作用下土體內聚力對輸氣管道最大軸向拉、壓應變的影響，土體內聚力對輸氣管道的最大軸向應變影響較小。其主要原因為：雖然土體內聚力是土體屈服的控制參數(shù)之一，但并非關鍵控制參數(shù)，其對土體屈服的影響較小，這點可從Mohr-Coulomb準則的直線包絡線中看出。在Mohr-Coulomb準則直線包絡線中，土體內聚力僅是直線的截距，直線斜率由土體內摩擦角的正切值控制。

圖5 土體內聚力對輸氣管道受拉側最大軸向應變的影響

4 結束語

本文推導了等效邊界彈簧設置在管道彈塑性變形段的本構關系，并據(jù)此建立了輸氣管道穿越管道地震斷層殼-接觸力學模型。在此基礎上，分析了土體力學性質對輸氣管道穿越地震斷層的應變響應規(guī)律。具體結論包括：

1）通過對比有限元結果與實驗結果，得知最大相對誤差僅為7.52%，驗證了所提輸氣管道穿越地震斷層有限元模型的準確性。

2）輸氣管道在0.5m的斷層位移下發(fā)生局部屈曲失效，而在1.0m以上的斷層位移下發(fā)生拉伸失效。

3）在0.5m的斷層位移下，土體彈性模量和內摩擦角對輸氣管道最大軸向應變有較大影響。但當斷層位移增大到2.5m時，土體彈性模量和內摩擦角對最大軸向應變的影響明顯降低。

[1]劉銘剛，楊秀娟，閆相禎，等．基于多重非線性管-土耦合模型的管道應變設計方法[J]．中國測試，2016，42(7)：9-14．

[2]KARAMITROSD K，BOUCKOVALASG D，KOURETZIS G P.Stress analysis of buried steel pipelines at strikeslip fault crossings[J].Soil Dyn Earthquake Eng，2007（27）：200-211.

[3]TRIFONOV O V，CHERNIY V P.A semi-analytical approach to a nonlinear stress-strain analysis of buried steel pipelines crossing active faults[J].Soil Dyn Earthquake Eng，2010（30）：1298-1308.

[4]趙林，馮啟民．埋地管線有限元建模方法研究[J]．地震工程與工程振動，2001，21（2）：53-57．

[5]劉愛文，胡聿賢，趙鳳新，等．地震斷層作用下埋地管線殼有限元分析的等效邊界方法[J]．地震學報，2004，26（增刊）：141-147．

[6]徐龍軍，劉慶陽，謝禮立．海底跨斷層輸氣管道動力特性數(shù)值模擬與分析[J].工程力學，2015，32（12）：99-107．

[7]李鴻晶，金瀏．穿越斷層埋地管線反應數(shù)值模擬——建模中的幾個具體問題[J]．自然災害學報，2011，20（3）：151-156．

[8]金瀏，李鴻晶．穿越逆沖斷層的埋地管道非線性反應分析[J]．防災減災工程學報，2010，30（2）：130-134．

[9]趙雷，唐暉，彭小波，等．埋地鋼管在逆斷層作用下失效模式研究[J]．應用基礎與工程科學學報，2010，18（增刊）：111-118．

[10]NEWMARK N M，HALL W J.Pipeline design to resist large fault displacement[C]∥Proceedings of US National Conference on Earthquake Engineering.Ann Arbor：University of Michigan，1975：416-425.

[11]KENNEDY R P，CHOW A W，WILLIAMSON R A. Fault movement effects on buried oil pipeline[J].Journal of the Transportation Engineering Division ASCE，1977，103（5）：617-633.

[12]WANG L R L，YEH Y A.Refined seismic analysis and design of buried pipeline for fault movement[J]. Journal of Earthquake Engineering and Structural Dynamics，1985，13（1）：75-96.

[13]TRIFONOV O V，CHERNIY V P.Elastoplastic stressstrain analysis of buried steel pipelines subjected to fault displacements with account for service loads[J]. Soil Dyn Earthquake Eng，2012（33）：54-62.

[14]VAZOURAS P，KARAMANOS S A，DAKOULAS P. Finite element analysis of buried steel pipelines under strike-slip fault displacements[J].Soil Dyn Earthquake Eng，2010（30）：1361-1376.

（編輯：李妮）

Effects of soil properties on gas pipelines crossing faults

ZHANG Lisong1，2，ZHANG Shiyan2，YAN Xiangzhen2，CAO Yuguang2
（1.Shandong Provincial Key Laboratory of Oil&Gas Storage and Transportation Safety，Qingdao Key Laboratory of Circle Sea Oil&Gas Storage and Transportation Technology，Qingdao 266580，China；2.College of Pipeline and Civil Engineering，China University of Petroleum，Qingdao 266580，China）

To analyze the influence of the soil properties on the strain response of the pipeline，the constitutive relation of the equivalent boundary spring acting on the segment where the pipeline yielding occurs was derived.Considering the equivalent boundary spring，the shell-contact finite element model of the pipeline crossing faults was established.Using the finite element model，the effect of the soil properties on the maximum axial strain on the tensile and compressive sides of the pipeline was investigated.The results show that the elastic modulus and internal friction angle of the soil have obvious effects on the maximum axial strain under smaller fault displacements.However，with the increase of the fault displacement，the effect of the elastic modulus and internal friction angle of the soil decreases significantly.Despite the degree of the fault displacement，the effect of the cohesion of the soil on the maximum axial strain is very low.

gas pipeline；eathquake fault；soil mechanical property；equivalent boundary spring；shell-contact model

A

1674-5124（2017）05-0120-04

10.11857/j.issn.1674-5124.2017.05.025

2016-08-20；

2016-11-10

國家自然科學基金（11472309）；CNPC重點實驗室課題（2014A-4214）

張立松（1982-），男，山東濰坊市人，副教授，博士，研究方向為油氣管道強度及穩(wěn)定性。