鄭 偉， 張 宏，2， 劉嘯奔， 陳嚴飛， 梁樂才

（1.中國石油大學（北京）機械與儲運工程學院，北京 102249；2.克拉瑪依理工學院籌建辦公室，新疆克拉瑪依834000)

基于應變設計方法在新粵浙斷層區(qū)管道上的應用*

鄭 偉1， 張 宏1，2， 劉嘯奔1， 陳嚴飛1， 梁樂才1

（1.中國石油大學（北京）機械與儲運工程學院，北京 102249；2.克拉瑪依理工學院籌建辦公室，新疆克拉瑪依834000)

新粵浙管道是我國即將建設的最長煤制氣管道，針對其穿越活動斷層管段采用基于應變的方法設計。介紹了管道基于應變設計準則，明確了基于應變設計校核標準及應用范圍?；诜蔷€性有限元法，給出了穿越活動斷層管道應變計算數(shù)值模型，模型采用管殼耦合的方法在保證計算精度的同時兼顧了經(jīng)濟性。結(jié)合管道基于應變設計流程，給出新粵浙管道穿越活動斷層實例的設計過程，為穿越斷層區(qū)管道設計建設提供了參考。

新粵浙管道；活動斷層；基于應變設計；有限元；抗震設計

0 前言

我國是地震損失最為嚴重的國家之一，地震引發(fā)的地表永久性變形會導致埋地管道運動甚至破壞，是考驗管道系統(tǒng)能否安全可靠運行的重要因素，也是管道設計者和運營者必須關注的問題。新粵浙管道工程即中國石化新疆煤制天然氣外輸管道工程，線路全長7 927 km，途經(jīng)13個?。ㄗ灾螀^(qū)），全線共有29處與活動斷層相交，如何安全地穿越活動斷層是管道工程建設的一個關鍵問題。本研究應用了基于應變的設計方法，明確基于應變設計的流程與校核標準，建立了管道穿越活動斷層有限元模型，對新粵浙斷層區(qū)管道進行設計，確保管道安全穿越活動斷層，為管道建設提供了參考。

1 管道基于應變設計準則

管道基于應變設計準則在國內(nèi)外現(xiàn)行的各類管道設計標準中都有所體現(xiàn)。管道基于應變設計準則是建立在極限狀態(tài)設計思想和位移載荷作用的基礎上，允許管道的應力超過材料屈服應力，此時管道雖然發(fā)生塑性變形，但在未超過許用應變的情況下，仍可以保證管道安全運營，充分發(fā)揮管道材料性能。針對活動斷層、滑坡、采空等位移控制載荷的不良地質(zhì)災害下進行管道設計，更宜采用基于應變設計準則。

基于應變的管道強度設計方法的數(shù)學表達式為

式中：ε—各種載荷作用時管道可能產(chǎn)生的最大應變；

εcrit—管材所能承受的極限應變，為設計系數(shù)；

F—設計系數(shù)；

[ε]—許用應變。如果管道設計應變的計算值不超過許用應變，則滿足安全要求，否則管道失效。

極限應變求取方法按照CSA－Z662－03：2006標準附錄中方法求取，見公式（2）和公式（3），

δ—表觀CTOD韌性；

λ—區(qū)相比；

ξ—缺欠長度與壁厚比率；

η—缺欠深度與壁厚比率；

t—管道壁厚;

D—管道外徑;

P—設計壓力;

Fy—有效屈服強度。

2 有限元模型

2.1 管材模型

新粵浙管道穿越活動斷層管段選用國產(chǎn)L555HD大變形鋼管，管材應力-應變拉伸曲線屬于“圓屋頂”形，不含屈服平臺，相對普通鋼材管道有更好的變形能力和更高的屈服應變。模型應力-應變數(shù)據(jù)選取管材真實應力-應變曲線，在多條試驗拉伸曲線中取最保守一條作為計算用曲線，如圖1所示。

圖1 管材拉伸應力－應變曲線

2.2 管土相互作用

管土相互作用是復雜的非線性問題，ALA《Guideline for the Design of Buried Steel Pipeline》中提出的使用非線性土彈簧描述管土作用的方法在工程中得到了廣泛應用。其采用3個方向的土彈簧描述土壤在軸向(KT)、側(cè)向(KP)和垂直方向(KQ)對管道的作用，作用方式如圖2(a)所示。土彈簧由極限作用力和屈服位移兩個量確定，如圖2(b)所示，其數(shù)值大小可由ALA規(guī)范求得。

2.3 單元選擇

圖2 土彈簧作用示意圖

殼單元主要用來模擬結(jié)構1個方向尺度遠小于其他方向尺度并忽略沿厚度方向應力時的情況，模型與斷層相交附近管道使用ABAQUS中的S4R殼單元模擬，其能夠很好地描述管道的受力情況。兩側(cè)管道使用PIPE31管單元模擬，管單元屬于梁單元的一種，由管道截面屬性定義單元對載荷的響應，計算成本小于殼單元。

在管道上的所有節(jié)點外建立同樣數(shù)量的土節(jié)點，并采用JOINTC單元連接管道與土節(jié)點來模擬土彈簧的作用，JOINTC單元為ABAQUS中的一個特殊連接單元，可以用來描述兩個節(jié)點之間隨相對位移變化產(chǎn)生的相互作用力。

2.4 模型建立

針對大直徑薄壁管道，使用殼單元可以準確模擬薄壁管道的力學行為。建立管道穿越活動斷層有限元模型，管道應足夠長，用以真實地模擬斷層兩側(cè)的受力情況，如果管道全由殼單元模擬，直接導致計算成本的升高。為了縮減計算代價，常用等效邊界的方法模擬管道遠端受力情況。

筆者采用管殼耦合模型分析管道受力，穿越斷層附近管道采用殼單元建模，殼單元管道模型長100 m，環(huán)向20個單元，軸向0.2 m長度一個單元；兩端采用管單元建模，長度分別取1 000 m。

殼單元管道兩端與管單元管道耦合，將管單元管道端部作為控制點，約束殼單元管道端部區(qū)域的全部6個自由度，此區(qū)域的各節(jié)點之間不會發(fā)生相對位移，只會隨著控制點做剛性運動，即將管端行為傳遞給殼單元管道邊界。通過管殼耦合的方式縮減了單元數(shù)量和計算成本。模型如圖3所示。

圖3 管殼耦合模型示意圖

通過與全殼模型對比驗證管殼耦合模型的正確性，取L555管道，管徑1 219mm，壁厚22mm，工作壓力10 MPa；全殼模型全長2 000 m，管道均由殼單元建模；對比管道90°穿越走滑斷層，位錯由0.5~4.5 m增加。

計算結(jié)果如圖4所示。管道穿越斷層產(chǎn)生變形和軸向拉伸，提取管道最大拉壓應變值，可以看出純殼模型與管殼耦合模型結(jié)果一致；由于管殼耦合模型單元數(shù)少于純殼模型，計算效率顯著提高，如圖5所示。

圖4 兩種數(shù)值模型結(jié)果對比圖

圖5 兩種數(shù)值模型計算成本示意圖

3 工程應用

3.1 基于應變設計流程

開展管道抗震設計，需要先校核其抗拉伸和壓縮能力。根據(jù)斷裂帶參數(shù)、土壤參數(shù)和管道性能參數(shù)等計算設計應變；利用經(jīng)驗公式和試驗等確定對應管道的極限應變，從而得到許用應變；當初始設計應變不符合規(guī)范要求時，采取對應抗震措施，直到設計應變小于許用應變，以保證管道安全運營。

3.2 案例背景

斷層參數(shù)：管道與開墾河斷裂帶相交，活動斷層傾角為50°，性質(zhì)為左旋逆斷層，預測地表最大位錯為垂直方向1.3 m，水平方向1.7 m，管道與活動斷層設計交角為83°。

管道參數(shù)：穿越管道管材為L555HD鋼，計算采用真實應力-應變拉伸曲線；管道初始設計參數(shù)為管徑1 219mm，設計壁厚22mm，設計壓力12 MPa，管道中心線埋深2.1 m。

土壤參數(shù)：穿越斷層管段管溝采用松砂回填，根據(jù)ALA《Guideline for the Design of Buried Steel Pipeline》中公式求得土彈簧參數(shù)見表1。

表1 土彈簧參數(shù)

許用應變：大變形L555HD鋼材管道的許用拉伸和壓縮應變求取參照CSA－Z662－03：2006中給出的經(jīng)驗公式，鋼材特征參數(shù)取值參考西氣東輸二線大變形管材數(shù)據(jù)，結(jié)果見表2。

表2 管道許用應變計算結(jié)果

3.3 設計過程

3.3.1 初步設計

初步設計首先計算設計工況下管道應變，不滿足許用應變時采取對管道淺埋、增大管道壁厚等措施。淺埋直接降低管土之間作用力，而增加壁厚加強了管道抵抗變形能力。計算結(jié)果見表3。

表3 初步設計校核結(jié)果 %

3.3.2 改變穿越角度

初步設計校核不滿足應變設計準則時，應考慮改變管道穿越活動斷層的交角，改變穿越角直接影響管道受力情況。實際管道基于應變設計時，斷層類型一般為組合斷層，即正逆斷層與走滑斷層的組合，不同于穿越單一走滑、正逆斷層，管道穿越組合斷層其最大應變與穿越角沒有統(tǒng)一規(guī)律可循。針對此，設計校核時計算不同穿越角度下的設計應變，擇優(yōu)選擇。

管道穿越開墾河斷裂改變穿越角度，結(jié)果如圖6所示。角度大于80°后由于壓應變過大，計算不收斂；軸向最大拉、壓應變隨著穿越角度增大而增大；當穿越角度小于60°時，管道滿足設計準則；原穿越角度為83°，綜合考慮安全設計及改線成本，管道以60°穿越為宜。

圖6 改變穿越角度應變變化結(jié)果

3.3.3 設計建議

基于應變設計準則，綜合設計過程，給出管道建設建議：管道穿越開墾河斷裂時，建議采用26.4mm壁厚L555HD管道，管道1.8 m淺埋并保證其與斷層交角為60°。

4 結(jié) 語

管道穿越活動斷層作為長輸管道建設及運營中的關鍵問題越來越受到重視，筆者介紹了管道基于應變設計準則，針對穿越活動斷層管道建立有限元模型并進行安全設計，得到以下結(jié)論：

（1）穿越活動斷層管道基于應變設計宜采用有限元方法計算設計應變；管殼耦合模型通過對與斷層相交殼單元管道和兩側(cè)對應管單元管道進行耦合，同時兼顧了計算精度與成本。

（2）通過新粵浙斷層區(qū)管道設計實例介紹了管道基于應變設計流程；當管道設計應變超出許用應變時，需采取抗震措施，包括松砂回填、增加管道鋼級、加大壁厚等；改變管道與斷層交角可以有效降低設計應變，對于空間斜滑斷層，最優(yōu)交角需多次試算得出。

［1］CSA－Z662－03:2006,Oil and Gas Pipeline Systems[S].

［2］American Lifelines Alliance： 2001， Guideline for the Design of Buried Steel Pipeline[S].

［3］GB 50470—2008，油氣輸送管道線路工程抗震技術規(guī)范[S].

［4］張宏，崔紅升.基于應變的管道強度設計方法的適用性[J].油氣儲運， 2012， 31(12)： 952－954.

［5］ABAQUS Inc.Abaqus Analysis User’s Manual Volume IV：Element[M].Providence USA:ABAQUS Inc.,2004.

［6］劉愛文.基于殼模型的埋地管線地震抗震分析[D].北京：中國地震局地球物理研究所，2002.

［7］余志峰，史航，佟雷，等.基于應變設計方法在西氣東輸二線的應用[J].油氣儲運,2010,29(02):143－147.

［8］李鶴林，吉玲康.西氣東輸二線高強韌性焊管及保障管道安全運行的關鍵技術[J].世界鋼鐵,2009,20(01)： 56－64.

［9］劉冰，劉學杰，張宏.基于應變的管道設計準則[J].天然氣工業(yè)， 2008， 28(02)： 129－131.

［10］李璞，陶燕麗，周建.基于應變設計管道局部屈曲應變極限值的計算[J].天然氣工業(yè)，2013，33(07)：101－107.

［11］劉嘯奔，陳嚴飛，張宏，等.受壓時跨斷層X80管道設計應變研究[J].天然氣工業(yè)，2014，34(12)：123－130.

［12］LIU Xiaoben， ZHANG Hong， CHEN Yanfei， et al.Design strain of X80 pipeline crossing active fault[C]//The 5th World Conference of Safety of Oil and Gas Industry.Okayama， Japan：[s.n.]， 2014.

[13]Q/SY GJX 0136—2008，西氣東輸二線管道工程強震區(qū)和活動斷層區(qū)段埋地管道基于應變設計導則[S].

Application of the Strain-based Design Method in Area of fault of XinYueZhe Pipeline

ZHENG Wei1,ZHANG Hong1,2,LIU Xiaoben1,LIANG Lecai1,FANG Maoli3
（1.College of Mechanical and Transportation Engineering,China University of Petroleum,Beijing 102249,China;2.Preparation Office of Karamay Institute of Technology,Karamay 834000,Xinjiang,China)

The XinYueZhe pipeline is the longest coal-to-gas pipeline to be constructed in China.The strain-based design method will be applied in the pipeline through active fault.In this article,the design criteria,check standard and applied range of strain-based design method were introduced.The numerical model of pipeline through active fault was set up based on the nonlinear finite element method.This model adopted pipe shell coupling method,which ensures the calculation accuracy and also with good efficiency.Combined with the strain-based design flow,the design process of the pipeline crossing the active area of fault instance was given,which provides reference for design and construction of pipeline through area of fault.

XinYueZhe pipeline;active fault;strain-based design method;finite element method;aseismatic design

TE832

B

1001－3938（2015）04-0038-05

國家自然科學基金（項目編號：51309236）；

教育部博士點基金（項目編號：20120007120009）；

中石化石油工程設計有限公司科學研究項目：基于應變設計技術研究（項目編號：2013406）。

鄭偉（1990—），男，碩士研究生，主要從事油氣管道失效分析、安全評價及石油石化安全研究。

2014－12－11

羅 剛