張利杰 賀新紅

（浙江省電力建設(shè)有限公司，浙江寧波 315000）

0 引言

天然氣是現(xiàn)今使用較為廣泛的清潔能源之一，在中國(guó)能源結(jié)構(gòu)調(diào)整中占據(jù)著關(guān)鍵位置，也是改善生態(tài)環(huán)境的主要途徑之一。調(diào)查研究數(shù)據(jù)顯示，截至2021年底，中國(guó)天然氣消費(fèi)量達(dá)到了3 726億m3，在能源消費(fèi)中占比較大，并呈現(xiàn)出逐年遞增的趨勢(shì)。由此可見(jiàn)，天然氣發(fā)展情況的好壞直接關(guān)系著中國(guó)社會(huì)經(jīng)濟(jì)發(fā)展?fàn)顩r。管道是天然氣運(yùn)輸?shù)年P(guān)鍵載體，目前，中國(guó)已經(jīng)完成全國(guó)管網(wǎng)的建設(shè)，主要包含主干線(xiàn)、聯(lián)絡(luò)線(xiàn)、省際線(xiàn)、配氣線(xiàn)等。隨著經(jīng)濟(jì)水平的不斷提升，天然氣供應(yīng)需求量急劇增加，對(duì)天然氣管道建設(shè)提出了更高的要求，管道能否穩(wěn)定、安全運(yùn)行成為影響天然氣供應(yīng)的關(guān)鍵所在[1]。

近年來(lái)，天然氣管道斜接管安全事故頻發(fā)，給相關(guān)部門(mén)留下了慘痛的教訓(xùn)。以某省管道爆炸事故為例，該事故造成了多人受傷與死亡，經(jīng)濟(jì)損失高達(dá)600萬(wàn)元，并致使事故區(qū)域天然氣供應(yīng)停滯約1個(gè)月時(shí)間，這表明中國(guó)在天然氣管道監(jiān)控方面還存在著很多的不足。對(duì)天然氣管道事故進(jìn)行統(tǒng)計(jì)分析，發(fā)現(xiàn)由斜接管強(qiáng)度不足引發(fā)的事故在頻率上占比第一，而且在提出與實(shí)施多條監(jiān)管措施的背景下，此類(lèi)事故發(fā)生頻率依舊居高不下，故提出了載荷作用下天然氣管道斜接管強(qiáng)度評(píng)估方法的研究。

1 天然氣管道斜接管強(qiáng)度評(píng)估方法研究

1.1 斜接管有限元模型構(gòu)建

為提升斜接管強(qiáng)度評(píng)估的精度，首要任務(wù)就是構(gòu)建斜接管有限元模型，并對(duì)其進(jìn)行有限元分析，為后續(xù)極限載荷的計(jì)算奠定堅(jiān)實(shí)的基礎(chǔ)[2]。

斜接管有限元分析結(jié)果與其單元形式之間存在著緊密的聯(lián)系，單元形式選取是否正確會(huì)影響強(qiáng)度評(píng)估結(jié)果的精度，由此可見(jiàn)，正確選取單元形式是有限元分析的基礎(chǔ)與前提。單元形式主要包括面、線(xiàn)、板殼和三維實(shí)體[3]，在斜接管有限元模型構(gòu)建過(guò)程中，通過(guò)每種單元形式優(yōu)劣勢(shì)的分析與探究，結(jié)合斜接管的實(shí)際情況，選取三維實(shí)體單元形式作為有限元模型構(gòu)建單元，其能夠清晰地模擬斜接管完整應(yīng)力分布情況，并且該單元形式適應(yīng)性較強(qiáng)，能與多種形狀相協(xié)調(diào)，更加能夠體現(xiàn)斜接管與管道連接位置的實(shí)際情況，為有限元分析提供支撐。

斜接管有限元模型構(gòu)建過(guò)程中，網(wǎng)格劃分也是一個(gè)至關(guān)重要的環(huán)節(jié)，網(wǎng)格質(zhì)量與數(shù)量是影響后續(xù)分析與計(jì)算的關(guān)鍵因素。為了提升構(gòu)建模型的精密性，將其分割為三個(gè)部分，分別為管道、斜接管及二者連接區(qū)域。需要注意的是，對(duì)管道與斜接管連接區(qū)域的網(wǎng)格劃分應(yīng)該更加精細(xì)。雖然精細(xì)網(wǎng)格劃分會(huì)延長(zhǎng)有限元模型構(gòu)建的時(shí)間，但卻能夠獲得精度較高的分析結(jié)果，對(duì)研究來(lái)說(shuō)，利大于弊，故應(yīng)最大限度地精細(xì)劃分斜接管有限元模型網(wǎng)格[4]。

為簡(jiǎn)化研究過(guò)程，利用斜接管結(jié)構(gòu)的對(duì)稱(chēng)性，依據(jù)上述選取的單元形式、確定的網(wǎng)格劃分規(guī)則，截取斜接管結(jié)構(gòu)的1/2進(jìn)行有限元建模，具體如圖1所示。

圖1中，L與D表示管道的長(zhǎng)度與內(nèi)徑；l與d表示斜接管的長(zhǎng)度與內(nèi)徑；θ表示斜接管與管道的夾角。

圖1 斜接管有限元模型示意圖

以上述構(gòu)建的斜接管有限元模型為基礎(chǔ)，對(duì)其載荷與邊界條件進(jìn)行深入探究，為后續(xù)極限載荷的計(jì)算提供支撐。在斜接管實(shí)際應(yīng)用過(guò)程中，其載荷來(lái)自多個(gè)方面的共同作用，但以?xún)?nèi)壓為主[5]。常規(guī)情況下，隨著內(nèi)壓的不斷加大，斜接管與管道連接區(qū)域受到的壓力也會(huì)增加，從而產(chǎn)生變形現(xiàn)象，若變形超過(guò)一定的閾值，則會(huì)導(dǎo)致管道破裂的事故。

斜接管載荷計(jì)算公式為：

由于構(gòu)建模型具有對(duì)稱(chēng)性，將斜接管斷面全部節(jié)點(diǎn)進(jìn)行旋轉(zhuǎn)變換，并對(duì)其進(jìn)行位移約束，即可獲得斜接管邊界信息，為后續(xù)研究的進(jìn)行提供便利。

上述過(guò)程完成了斜接管有限元模型構(gòu)建，奠定了研究進(jìn)行的基礎(chǔ)與前提。

1.2 斜接管極限載荷計(jì)算

以上述構(gòu)建的斜接管有限元模型為基礎(chǔ)，衡量斜接管應(yīng)力集中系數(shù)，以此確定斜接管的極限載荷，為最終斜接管強(qiáng)度評(píng)估的實(shí)現(xiàn)提供依據(jù)。

應(yīng)力集中系數(shù)反映的是在載荷作用下斜接管局部范圍內(nèi)應(yīng)力明顯增大的程度，計(jì)算公式為：

式中：β為應(yīng)力集中系數(shù)；υθ為實(shí)時(shí)應(yīng)力值；αθ′為最大基本應(yīng)力；P為極限載荷；k為斜接管管道材料厚度。

在斜接管應(yīng)用過(guò)程中，隨著載荷數(shù)值的不斷增加，斜接管與管道連接區(qū)域最先開(kāi)始屈服，應(yīng)力由此向外擴(kuò)展。但是，管道具有一定的彈性，當(dāng)載荷數(shù)值在閾值以下時(shí)，斜接管處于彈性階段，應(yīng)力分布呈現(xiàn)一定的規(guī)律性，并不會(huì)對(duì)管道安全產(chǎn)生較大的影響[6]。當(dāng)載荷數(shù)值超過(guò)閾值時(shí)，斜接管處于塑性階段，應(yīng)力分布呈現(xiàn)出無(wú)規(guī)律性，管道開(kāi)始出現(xiàn)形變。形變?nèi)羰沁^(guò)大，管道就會(huì)發(fā)生破裂，導(dǎo)致天然氣泄漏，如在此時(shí)遇到明火，就會(huì)發(fā)生爆炸事故。需要注意的是，應(yīng)力變化情況與斜接管實(shí)際結(jié)構(gòu)沒(méi)有關(guān)系，但其與斜接管、管道的夾角θ存在著較為密切的關(guān)系。常規(guī)情況下，應(yīng)力值υθ隨著夾角θ的增大而減小，兩者呈現(xiàn)明顯的負(fù)相關(guān)關(guān)系。因此，在管道建設(shè)過(guò)程中，應(yīng)該盡可能采用夾角大于30°的斜接管結(jié)構(gòu)。

與此同時(shí)，隨著載荷數(shù)值的增加，斜接管變形程度加大。獲取斜接管應(yīng)力分布數(shù)值，將其代入公式（2），經(jīng)過(guò)方程式求解，即可獲得斜接管極限載荷P，為后續(xù)研究的進(jìn)行提供助力。

1.3 斜接管內(nèi)腐蝕預(yù)測(cè)

隨著應(yīng)用時(shí)間的增加，斜接管內(nèi)部會(huì)產(chǎn)生腐蝕現(xiàn)象，這也是影響其強(qiáng)度的關(guān)鍵因素，故此節(jié)對(duì)斜接管內(nèi)部腐蝕情況進(jìn)行預(yù)測(cè)，為最終強(qiáng)度評(píng)估的實(shí)現(xiàn)提供依據(jù)[7]。

根據(jù)已有文獻(xiàn)研究可知，造成管道內(nèi)部腐蝕現(xiàn)象的主要參數(shù)為溫降、壓降與平均腐蝕速率。其中，溫降計(jì)算公式為：

式中：T為溫降參數(shù)數(shù)值；T0為管道內(nèi)部初始溫度；TQ為斜接管內(nèi)部實(shí)時(shí)溫度；M為輔助參數(shù)；CP為斜接管應(yīng)力變化因子；P0為斜接管初始應(yīng)力數(shù)值。

壓降參數(shù)無(wú)法直接測(cè)量獲得，需要通過(guò)下述方程計(jì)算而得：

式中：Z為壓降參數(shù)數(shù)值；R為斜接管內(nèi)部半徑數(shù)值；g為重力加速度；χ為輔助運(yùn)算因子，取值范圍為0～1。

平均腐蝕速率需要依據(jù)實(shí)際天然氣運(yùn)輸情況進(jìn)行預(yù)測(cè)，表達(dá)式為：

式中：Vτ為平均腐蝕速率數(shù)值；p與q為隨機(jī)整數(shù)；F（v）為天然氣運(yùn)輸過(guò)程中管道的損失系數(shù)；G（t）為腐蝕時(shí)間系數(shù)；Vc為標(biāo)準(zhǔn)腐蝕速率。

依據(jù)公式（3）（4）（5）計(jì)算結(jié)果，預(yù)測(cè)斜接管內(nèi)腐蝕情況，為了方便研究，將其轉(zhuǎn)換為斜接管內(nèi)壁面損失，表達(dá)式為：

式中：ζ為斜接管內(nèi)壁面損失；σ2為斜接管內(nèi)腐蝕程度參量。

一般情況下，壁面損失越大，則斜接管厚度越薄，其強(qiáng)度也就越差。

上述過(guò)程完成了斜接管內(nèi)部腐蝕情況的預(yù)測(cè)，為研究目標(biāo)的實(shí)現(xiàn)提供了支撐。

1.4 斜接管強(qiáng)度評(píng)估實(shí)現(xiàn)

以上述獲得的斜接管內(nèi)壁面損失數(shù)值為基礎(chǔ)，制定斜接管強(qiáng)度評(píng)估程序與規(guī)則，獲取精準(zhǔn)的斜接管強(qiáng)度結(jié)果，為天然氣管道的安全使用提供保障。

斜接管強(qiáng)度評(píng)估是天然氣管道運(yùn)行過(guò)程中的關(guān)鍵環(huán)節(jié)之一，是管道檢查與維修的主要依據(jù)。依據(jù)斜接管實(shí)際作業(yè)情況，制定強(qiáng)度評(píng)估程序，具體如圖2所示。

圖2 斜接管強(qiáng)度評(píng)估程序圖

如圖2程序所示，斜接管安全系數(shù)的計(jì)算是強(qiáng)度評(píng)估的主要數(shù)據(jù)依據(jù)，計(jì)算公式為：

式中：Г為斜接管安全系數(shù)；ξm為腐蝕預(yù)測(cè)的分項(xiàng)安全系數(shù)；t為時(shí)間；η為固定的最小拉伸強(qiáng)度；ξd為腐蝕深度的分項(xiàng)安全系數(shù)；Q為長(zhǎng)度校正系數(shù)。

以公式（7）計(jì)算結(jié)果為基礎(chǔ)，制定強(qiáng)度評(píng)估規(guī)則，具體如下：

（1）當(dāng)Г處于[0.0，0.3]范圍內(nèi)時(shí)，認(rèn)定斜接管強(qiáng)度較差；

（2）當(dāng)Г處于[0.4，0.5]范圍內(nèi)時(shí)，認(rèn)定斜接管強(qiáng)度一般；

（3）當(dāng)Г處于[0.6，0.8]范圍內(nèi)時(shí)，認(rèn)定斜接管強(qiáng)度尚可；

（4）當(dāng)Г處于[0.9，1.0]范圍內(nèi)時(shí)，認(rèn)定斜接管強(qiáng)度較好。

依據(jù)上述制定的評(píng)估規(guī)則，即可實(shí)現(xiàn)天然氣管道斜接管強(qiáng)度的精準(zhǔn)評(píng)估，為天然氣穩(wěn)定傳輸提供助力。

2 實(shí)驗(yàn)與結(jié)果分析

2.1 實(shí)驗(yàn)對(duì)象概況

選取某段天然氣管道作為實(shí)驗(yàn)對(duì)象，其長(zhǎng)度約為8 700 km，內(nèi)部包含25個(gè)斜接管，并且斜接管應(yīng)用的時(shí)間長(zhǎng)度存在著明顯的差異。管道采用的是國(guó)產(chǎn)X80管線(xiàn)鋼，內(nèi)部包含Mo、Cr、Nb、Ni、P、C等元素，具有較好的耐腐蝕性能，可以保障天然氣管道安全、經(jīng)濟(jì)、高效運(yùn)行。

隨著天然氣供給需求的增加，實(shí)驗(yàn)對(duì)象——天然氣管道面臨著較大的挑戰(zhàn)，對(duì)斜接管性能的要求更高，因此需要對(duì)斜接管強(qiáng)度進(jìn)行實(shí)時(shí)監(jiān)測(cè)，符合所提方法應(yīng)用性能的測(cè)試需求。

2.2 實(shí)驗(yàn)結(jié)果分析

以上述選取的實(shí)驗(yàn)對(duì)象為依據(jù)，進(jìn)行天然氣管道斜接管強(qiáng)度評(píng)估實(shí)驗(yàn)。為凸顯本文所提方法的應(yīng)用性能，選取含等壁厚體積型缺陷油氣管道的剩余強(qiáng)度評(píng)價(jià)[8]與艦船內(nèi)腐蝕海水管路剩余強(qiáng)度預(yù)測(cè)模型及試驗(yàn)驗(yàn)證[9]作為對(duì)比方法1與2，并選取斜接管極限載荷計(jì)算誤差與斜接管強(qiáng)度評(píng)估結(jié)果作為評(píng)價(jià)指標(biāo)。

通過(guò)實(shí)驗(yàn)獲得斜接管極限載荷計(jì)算誤差數(shù)據(jù)如圖3所示。

圖3 斜接管極限載荷計(jì)算誤差數(shù)據(jù)圖

如圖3數(shù)據(jù)所示，相較于兩種對(duì)比方法，應(yīng)用本文方法獲得的斜接管極限載荷計(jì)算誤差數(shù)值更小，最小值為0.2%。

通過(guò)實(shí)驗(yàn)獲得斜接管強(qiáng)度評(píng)估結(jié)果如表1所示。

表1 斜接管強(qiáng)度評(píng)估結(jié)果表

如表1內(nèi)容所示，應(yīng)用本文方法獲得的斜接管強(qiáng)度評(píng)估結(jié)果與實(shí)際情況保持一致，而對(duì)比方法1與2獲得的斜接管強(qiáng)度評(píng)估結(jié)果與實(shí)際情況存在較大偏差，說(shuō)明應(yīng)用本文方法獲得的斜接管強(qiáng)度評(píng)估精度較高。

上述實(shí)驗(yàn)數(shù)據(jù)顯示，與兩種對(duì)比方法相比較，應(yīng)用本文方法獲得的斜接管極限載荷計(jì)算誤差較小，斜接管強(qiáng)度評(píng)估結(jié)果與實(shí)際情況趨于一致，充分證實(shí)了本文方法應(yīng)用性能更好。

3 結(jié)語(yǔ)

天然氣是現(xiàn)今社會(huì)應(yīng)用極為廣泛的清潔能源之一，其主要通過(guò)管道進(jìn)行運(yùn)輸，管道是否安全直接關(guān)系天然氣企業(yè)的發(fā)展，故提出載荷作用下的天然氣管道斜接管強(qiáng)度評(píng)估方法研究。本文所提方法有效降低了斜接管極限載荷計(jì)算誤差，提升了斜接管強(qiáng)度評(píng)估精度，為天然氣管道管控提供了更加精準(zhǔn)的數(shù)據(jù)支撐，也為相關(guān)研究提供了一定的借鑒與參考。