鄭 偉， 張 宏，2， 劉嘯奔， 陳嚴(yán)飛， 梁樂才

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

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

鄭 偉1， 張 宏1，2， 劉嘯奔1， 陳嚴(yán)飛1， 梁樂才1

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

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

新粵浙管道；活動斷層；基于應(yīng)變設(shè)計；有限元；抗震設(shè)計

0 前言

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

1 管道基于應(yīng)變設(shè)計準(zhǔn)則

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

基于應(yīng)變的管道強(qiáng)度設(shè)計方法的數(shù)學(xué)表達(dá)式為

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

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

F—設(shè)計系數(shù)；

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

極限應(yīng)變求取方法按照CSA－Z662－03：2006標(biāo)準(zhǔn)附錄中方法求取，見公式（2）和公式（3），

δ—表觀CTOD韌性；

λ—區(qū)相比；

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

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

t—管道壁厚;

D—管道外徑;

P—設(shè)計壓力;

Fy—有效屈服強(qiáng)度。

2 有限元模型

2.1 管材模型

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

圖1 管材拉伸應(yīng)力－應(yīng)變曲線

2.2 管土相互作用

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

2.3 單元選擇

圖2 土彈簧作用示意圖

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

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

2.4 模型建立

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

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

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

圖3 管殼耦合模型示意圖

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

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

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

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

3 工程應(yīng)用

3.1 基于應(yīng)變設(shè)計流程

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

3.2 案例背景

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

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

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

表1 土彈簧參數(shù)

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

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

3.3 設(shè)計過程

3.3.1 初步設(shè)計

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

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

3.3.2 改變穿越角度

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

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

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

3.3.3 設(shè)計建議

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

4 結(jié) 語

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

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

（2）通過新粵浙斷層區(qū)管道設(shè)計實例介紹了管道基于應(yīng)變設(shè)計流程；當(dāng)管道設(shè)計應(yīng)變超出許用應(yīng)變時，需采取抗震措施，包括松砂回填、增加管道鋼級、加大壁厚等；改變管道與斷層交角可以有效降低設(shè)計應(yīng)變，對于空間斜滑斷層，最優(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，油氣輸送管道線路工程抗震技術(shù)規(guī)范[S].

［4］張宏，崔紅升.基于應(yīng)變的管道強(qiáng)度設(shè)計方法的適用性[J].油氣儲運(yùn)， 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］余志峰，史航，佟雷，等.基于應(yīng)變設(shè)計方法在西氣東輸二線的應(yīng)用[J].油氣儲運(yùn),2010,29(02):143－147.

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

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

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

［11］劉嘯奔，陳嚴(yán)飛，張宏，等.受壓時跨斷層X80管道設(shè)計應(yīng)變研究[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，西氣東輸二線管道工程強(qiáng)震區(qū)和活動斷層區(qū)段埋地管道基于應(yīng)變設(shè)計導(dǎo)則[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

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

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

中石化石油工程設(shè)計有限公司科學(xué)研究項目：基于應(yīng)變設(shè)計技術(shù)研究（項目編號：2013406）。

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

2014－12－11

羅 剛