石油加氫裝置放空管道的動態(tài)分析

申登巒 魏安安 王 瑤

(常州大學機械工程學院，江蘇 常州 213016)

石油加氫裝置放空管道的動態(tài)分析

申登巒 魏安安 王 瑤

(常州大學機械工程學院，江蘇 常州 213016)

化工裝置在突然放空時，放空氣體會對管道產(chǎn)生很大的推力，從而危及管道安全。結(jié)合某工程實例，采用CAESARII軟件對放空管道進行動態(tài)分析，應用模態(tài)分析計算出了管道系統(tǒng)的固有頻率，應用力響應頻譜法分析了管道的安全閥載荷和氣錘載荷對管道系統(tǒng)的影響，找到了受氣錘載荷和安全閥載荷沖擊影響最大、振動最劇烈的位置，并對此位置進行改進。據(jù)改進后的方案重新建模并分析，管道振動情況明顯減緩并安全地通過了模擬工況下的放空試驗。

放空管道 動態(tài)分析 力響應頻譜法 CAESARII軟件

放空管道用于石化裝置在突然停電、超壓及出現(xiàn)故障等意外情況時，通過安全閥的開啟，釋放系統(tǒng)中的介質(zhì)并降壓，以保護裝置的安全，這也是保護裝置安全的最后一道屏障。

管道應力分析軟件CAESARII可以進行管道的靜力和動力分析，功能齊全，因而在全世界范圍內(nèi)被廣泛應用[1,2]。

某石化廠為了節(jié)省占地，在技術(shù)改造中將現(xiàn)有的兩個石油加氫裝置的放空管道合并到一條原有放空管道中。筆者采用軟件CAESARII對管道進行靜動態(tài)分析，以校核此管線的安全性。

石油加氫裝置在安全閥開啟時會產(chǎn)生巨大的氣體排放反力，使得管道安全存在很大隱患。放空管內(nèi)介質(zhì)和安全閥的基本參數(shù)如下：

容器內(nèi)介質(zhì) 飽和蒸汽

容器內(nèi)最高工作溫度 250℃

安全閥型號 HF-CA48Y-100C

安全閥起跳壓力 4 170kPa

安全閥喉部直徑 80mm

安全閥公稱直徑 100mm

泄放量 8.26kg/s

放空管直徑 200mm

放空管出口處的壓強 497kPa

安全閥出口管的直徑 150mm

現(xiàn)用兩種方法計算氣體排放反力。

API-520中計算排氣反力的公式[1]為：

(1)

式中A——泄放出口面積，mm2；

F——泄放反力，N；

p——泄放管出口靜壓，kPa；

T——出口溫度，K；

w——氣體流量，kg/s。

將放空管道各項數(shù)據(jù)代入式(1)可得F=46033.19N。

另外，采用CAESARII軟件計算排氣反力。將管道內(nèi)介質(zhì)和安全閥的參數(shù)導入CASEARII中的Relief Load Synthesis模塊計算所得安全閥放空時的氣體排放反力F=36108.57N。

兩種方法的計算結(jié)果有一定差異，用API-520中的方法計算的結(jié)果較大。為保守起見，筆者取較大的數(shù)值應用于管道的動態(tài)分析，在進行管道系統(tǒng)的動態(tài)分析之前，必須先進行靜態(tài)分析。

2 管道動態(tài)分析

在項目完成初期，該廠做過一次小型放空試驗，發(fā)現(xiàn)第五個彎頭(即節(jié)點249)處振動劇烈(振幅近200mm)，嚴重超過了國家標準?，F(xiàn)著重介紹此管道系統(tǒng)的動態(tài)分析，在進行管道系統(tǒng)的動態(tài)分析前，必須先進行靜態(tài)分析。

2.1靜態(tài)分析

建立管道系統(tǒng)的模型，管道模型與節(jié)點分布如圖1所示。經(jīng)分析管道的一次應力、二次應力和管道熱應力均未超出應力許用值。

圖1 管道模型與節(jié)點分布

2.2動態(tài)分析

影響放空管道振動的主要因素是氣錘載荷和安全閥載荷，二者特性曲線的特點都是載荷值從零突然躍升到某一定值，然后保持一定時間后再突然降到零[3]。求解這種具有脈沖力/時間特性曲線的最常用方法是時程法和力頻譜法。時程分析采用數(shù)值積分方法模擬載荷的作用，這種方法需要迭代，人工計算部分較為繁瑣，因此工程上常用力響應頻譜法計算氣錘對管道的影響[4]。

2.2.1模態(tài)分析

模態(tài)解僅進行本征解計算，在頻譜分析前也需要進行本征解的計算。模態(tài)分析只是為了求解管道系統(tǒng)的固有頻率和振型，因此不需要定義載荷，只需將法蘭剛度輸入即可運算[5,6]。經(jīng)計算，此管線前五階的固有頻率為1.073、1.680、2.254、2.550、3.112Hz。此排放管系的一階固有頻率小于DL/T 5054-1996《火力發(fā)電廠汽水管道設(shè)計技術(shù)》中規(guī)定的一階固有頻率的最小值3.500Hz，極易受外界干擾發(fā)生強烈振動。

2.2.2安全閥排放反力分析

將上述求得的安全閥排放反力依次施加到管道彎頭處，并在每個彎頭處采用CAESARII軟件建立安全閥放空時氣體排放反力的響應頻譜，找到安全閥排放反力對管道影響最強的力，并用此力校核管道和支架應力[7]。以節(jié)點249為例生成的頻譜曲線如圖2所示?？梢钥闯觯l率為0.586Hz時動態(tài)載荷系數(shù)最大(1.997)，即此時安全閥排放的氣體反力對管道的影響最大。

圖2 節(jié)點249處排放反力與時間響應頻譜

計算安全閥排氣工況下固定支架上的受力，結(jié)果見表1。

表1 力作用在節(jié)點249處固定支架的受力

由以上的力，校核管道在安全閥排放工況下固定支架上的應力。

從CAESARII軟件的輸出文件中分別讀取排放反力施加到每個彎頭上時，各支架上受到的動態(tài)推力載荷，并進行管道與支架的應力校核。經(jīng)校核，管道和支架的應力均未超過許用值。

2.2.3氣錘力分析

用文獻[7]中的公式計算得到這條放空管線的氣錘力的大小、方向和作用位置見表2。

表2 氣錘力大小、方向與作用位置

從廠家提供的氣錘脈動時間歷程波形圖(圖3)中，讀出壓力波從零到達極值的時間為0.09s。載荷的持續(xù)時間計算式為：

t=L/c

(2)

式中c——壓力波傳播速度；

L——相鄰彎頭-彎頭對之間的距離。

圖3 氣錘脈動波形

以節(jié)點249處的彎頭為例，在節(jié)點249處產(chǎn)生最大的變形彎頭-彎頭對是249和429。對于249彎頭-彎頭對在6 000mm單元管長計算的最大截止頻率：

fcutoff=SQRT(305E/ρ)/L=136.9Hz

(3)

式中E——管材的彈性模量；

ρ——管材的密度。

取最大頻率150Hz，用CAESARII軟件生成的頻譜分析如圖4所示，頻譜數(shù)據(jù)見表3。

圖4 節(jié)點處氣錘力與時間響應頻譜

階次頻率/Hz動載荷系數(shù)80．75851．371091．21501．8242101．85191．9549112．71131．9049123．84001．8148

由圖4和表3可以看出，在此管系中管道頻率為1.851 9～2.711 3Hz時氣錘力有明顯加強，氣錘力對管道的影響最為明顯。

對每個彎頭對進行上述分析，并輸出各彎頭的位移量，將每個彎頭的最大位移量列于表4，放空過程中管道波動形態(tài)如圖5所示?？梢钥闯?，在氣錘力的作用下節(jié)點249處的彎頭位移量最大(340.33mm)，這與先前的放空試驗結(jié)果吻合。

表4 各彎頭的最大位移量

圖5 管道波動形態(tài)

3 管道系統(tǒng)的優(yōu)化

分析管系受氣錘載荷的影響，發(fā)現(xiàn)在節(jié)點249處的彎頭振動最為劇烈，同時發(fā)現(xiàn)整個管系的固有頻率較低，為此在節(jié)點250處增加一個固定約束，并在管道與支座間增加一個大于5mm的橡膠墊片。在節(jié)點619處(即圖1中的三通處)增加一個支撐約束。

按照修改后的方案，應用CASEARII軟件重新建模，并對管道進行靜態(tài)和動態(tài)分析，結(jié)果表明管道與各支架應力符合國家標準，各彎頭偏移量均在允許范圍內(nèi)，彎頭的最大偏移量縮小到50mm。最后對該放空管道進行模擬真實工況的泄放試驗，結(jié)果表明管道與支架系統(tǒng)均穩(wěn)定安全，管道振動幅度明顯減小，未見明顯的屈服、偏移和變形。目前該放空管道已經(jīng)投入使用，未發(fā)生異常現(xiàn)象。

4 結(jié)束語

結(jié)合某石化廠裝置放空管道的具體實例，介紹了在裝置放空時氣體對管道沖擊力的計算方法，以及采用CAESARII軟件對放空管道的動態(tài)模擬過程，找出了管道振動劇烈的原因，并提出了改進方案。為優(yōu)化管道設(shè)計與布置，以及優(yōu)化運行方式提供了依據(jù)。

[1] 趙晉云,周興濤,劉冰,等.國內(nèi)外輸氣管道放空系統(tǒng)設(shè)計標準分析[J].油氣儲運,2013,32(3):274～278.

[2] 葉源,吳勇.聚乙烯裝置工業(yè)高壓管線設(shè)計及應力分析[J].石油化工技術(shù)與經(jīng)濟,2012,28(5):31～33.

[3] 吳曉南,盧泓方,黃坤,等.基于頻譜分析的地震帶輸氣管道應力分析方法[J].天然氣工業(yè),2014,34(5):152～157.

[4] 張樂川,胡訓棟,李琳.火力發(fā)電廠主要管道精確動態(tài)分析方法及應用[J].發(fā)電設(shè)備,2011,25(6):402～407.

[5] 喻斌,張世彬,左勇,等.CAESAR II用于輸氣管道法蘭泄漏的治理[J].油氣田地面工程,2013,32(3):32～33.

[6] 李方圓,江濤,李玉星,等.川氣東送管道放空參數(shù)的計算[J].油氣儲運,2010,29(6):427～429.

[7] 吳曉南,鮮燕,黃坤,等.運行工況下隧道內(nèi)輸氣管道的應力分析[J].油氣儲運,2012,31(12):927～930.

DynamicAnalysisofVentPipesinOilHydrogenationUnit

SHEN Deng-luan, WEI An-an, WANG Yao

(SchoolofMechanicalEngineering,ChangzhouUniversity,Changzhou213016,China)

Sudden venting operation can incur great thrust on the pipelines and endanger their security. Taking a project as an example, dynamically analyzing the vent pipe was performed with CAESARII software, including having their natural frequency calculated through modal analysis method and the influence of pipeline relief valve’s load and the steam hammer’s load on the pipeline system analyzed by employing a force response spectrum method so as to determine the location where the largest thrust from both hammer’s load and relieve valve’s load and the most violent vibration exist and can be improved. Having improved scheme based to model and analyze the vibration of pipelines shows that pipeline vibration becomes decreased and the vent pipeline system safely passed the vent test under the simulated conditions.

vent pipe, dynamic analysis, force response spectrum method, CAESARII software

TQ055.8

A

1000-3932(2016)02-0173-04

2015-11-15(修改稿)