偶國富，AI-barakani A. Ahmed，許 健，葉浩杰，金浩哲

(浙江理工大學(xué)流動腐蝕研究所，杭州 310018)

?

常壓塔頂換熱器出口管道沖蝕特性的數(shù)值模擬

偶國富，AI-barakani A. Ahmed，許 健，葉浩杰，金浩哲

(浙江理工大學(xué)流動腐蝕研究所，杭州 310018)

針對石化工業(yè)常壓塔頂(常頂)系統(tǒng)中換熱器出口管道的沖蝕失效問題，分別利用工藝仿真軟件Aspen和計算流體力學(xué)軟件，進(jìn)行工藝計算和熱器出口管道的三維流場數(shù)值模擬，得到多相流中腐蝕性介質(zhì)的分布規(guī)律和流場結(jié)果。研究發(fā)現(xiàn)：管道內(nèi)油氣水三相流中均存在腐蝕性介質(zhì)，其溶于水后形成的腐蝕性溶液對管道壁面產(chǎn)生腐蝕作用，生成的腐蝕產(chǎn)物保護(hù)膜在壁面剪切應(yīng)力的作用下快速的脫落、再生，進(jìn)而加速了管道的腐蝕破壞；流場中水相主要集中在管道的外側(cè)，水相分率由外側(cè)壁面至內(nèi)側(cè)壁面逐漸降低，在腐蝕性溶液聚集的外側(cè)壁面，各彎管和直管段剪切應(yīng)力沿流動方向逐漸增大；沿流動方向第四只彎管內(nèi)側(cè)由入口至30°之間區(qū)域、其它四只彎管外側(cè)壁面沿流向30°至出口段和第三、四只彎管間直管外側(cè)壁面，為水相分率和壁面剪切應(yīng)力最大區(qū)域，即在腐蝕和流體剪切作用下失效的高風(fēng)險區(qū)域，仿真結(jié)果與管道測厚結(jié)果基本吻合。

換熱器管道；沖蝕特性；工藝分析；數(shù)值模擬

0 引 言

石化工業(yè)是國民經(jīng)濟(jì)的基礎(chǔ)工業(yè)之一，是保障能源安全、實現(xiàn)本世紀(jì)中葉達(dá)到中等發(fā)達(dá)國家這一戰(zhàn)略目標(biāo)的重要組成部分。受高硫、高酸等劣質(zhì)原油加工比例增加的影響，設(shè)備系統(tǒng)的嚴(yán)重腐蝕已成為煉油工業(yè)安全生產(chǎn)的瓶頸[1-3]。常減壓蒸餾作為石油煉制的首道工序，原油劣質(zhì)化帶來的弊端在常減壓裝置上被直接體現(xiàn)出來[4-7]。近年來，常壓塔頂(常頂)系統(tǒng)已多次出現(xiàn)換熱器管束穿孔、原油內(nèi)漏等事故，不僅嚴(yán)重影響企業(yè)的生產(chǎn)計劃，而且嚴(yán)重威脅節(jié)能、環(huán)保與安全。

常頂系統(tǒng)的主要腐蝕失效機(jī)理有多相流沖蝕、胺/銨鹽沉積垢下腐蝕、HCl露點腐蝕等。其中研究沖蝕失效國外最具代表性的是美國俄亥俄州立大學(xué)，主要采用電化學(xué)阻抗譜法、極化曲線法和電化學(xué)噪聲技術(shù)，研究復(fù)雜多相流環(huán)境下的沖蝕特性[8]；Shargay等[9]研究了銨鹽的垢下腐蝕機(jī)理，提出從原料含量、管道材質(zhì)、優(yōu)化設(shè)計等方面緩解銨鹽腐蝕失效；金浩哲等[10-11]采用熱力學(xué)推導(dǎo)的方法，建立了NH4Cl和NH4HS結(jié)晶平衡曲線；Dewan等[12]研究了胺鹽(DEA·HCl)對碳鋼的腐蝕速率，其所測得的DEA·HCl平衡與動力學(xué)數(shù)據(jù)，可為常減壓裝置加注中和劑DEA優(yōu)化提供指導(dǎo)；Zarenezhad等[13]建立了含硫燃燒氣體露點預(yù)測的三層四元人工神經(jīng)網(wǎng)絡(luò)模型，可針對酸性氣體的露點溫度進(jìn)行預(yù)測。與國外比較而言，我國針對常減壓裝置的腐蝕研究相對較晚，在原油劣質(zhì)化過程中存在的腐蝕問題更加突出。

為了揭示常減壓裝置的腐蝕失效機(jī)理，本文以常頂系統(tǒng)為研究對象，確定常頂系統(tǒng)的腐蝕回路，通過Aspen仿真分析多相流中腐蝕介質(zhì)的物性參數(shù)和分布規(guī)律；在此基礎(chǔ)上，針對常減壓裝置換熱器出口管道進(jìn)行結(jié)構(gòu)離散建模，通過流體動力學(xué)方法分析出口管道系統(tǒng)的沖蝕風(fēng)險。以仿真分析結(jié)果為基礎(chǔ)，對現(xiàn)場管道壁面相應(yīng)危險區(qū)域進(jìn)行壁厚檢測，驗證了數(shù)值模擬的準(zhǔn)確性，為常頂系統(tǒng)易腐蝕管道的測厚布點優(yōu)化提供參考依據(jù)。

1 常頂系統(tǒng)腐蝕回路分析

1.1 常頂系統(tǒng)腐蝕回路簡介

本文以某石化企業(yè)常頂系統(tǒng)換熱器出口管道為研究對象，常頂系統(tǒng)的工藝流程圖如圖1所示。原油在常壓塔T102中進(jìn)行蒸餾，常壓塔頂部分離出油氣混合物，混合物在進(jìn)入換熱器前進(jìn)行注氨、緩蝕劑和中和劑，之后分5路分別進(jìn)入換熱器E101/1-E101/5，5臺換熱器前均設(shè)有注水點?；旌衔锝?jīng)換熱器冷卻后進(jìn)入常頂回流罐V102進(jìn)行油、氣、水三相的分離，回流罐中分離后的油相在回流泵P103的作用下，一部分經(jīng)空冷器AC102和水冷器wc101/1、wc101/2冷卻后作為常頂一級油出裝置，另一部分進(jìn)入常壓塔頂部作為常頂回流油回流。回流罐中的氣相經(jīng)空冷器AC101和水冷器wc101/3、wc101/4進(jìn)行冷卻后，進(jìn)入分離罐V103再次進(jìn)行三相分離，分別得到常頂二級油、常頂瓦斯氣和含硫污水?；亓鞴轛102與V103分離出的含硫污水進(jìn)入污水處理罐V110進(jìn)行處理，處理后一部分常頂含硫污水經(jīng)泵P121流出裝置，另一部分作為冷卻水經(jīng)泵P122后進(jìn)入換熱器與常頂混合物進(jìn)行換熱。

圖1 常頂系統(tǒng)工藝流程示意圖

1.2 常頂系統(tǒng)腐蝕介質(zhì)分析

常頂系統(tǒng)中的腐蝕介質(zhì)主要有HCl與H2S。由于目前我國大量煉制國外進(jìn)口的劣質(zhì)原油，原油中存在一定量的硫化物及氯化物等雜質(zhì)。在加工過程中，隨著溫度逐步升高，原油中的MgCl2、CaCl2在水溶液的環(huán)境下逐步發(fā)生水解反應(yīng)，從而產(chǎn)生HCl，反應(yīng)過程如式(1)和式(2)所示。同時原油中開采中添加的有機(jī)氯(如CHCl3)也會因溫度的升高發(fā)生水解反應(yīng)產(chǎn)生HCl。而常頂系統(tǒng)中的H2S主要來源于原油雜質(zhì)中硫化物的受熱分解。

(1)

(2)

隨著原油蒸餾過程的不斷進(jìn)行，常頂系統(tǒng)逐漸產(chǎn)生HCl及H2S氣體，進(jìn)入塔頂位置及塔頂冷凝區(qū)，當(dāng)HCl與H2S以氣體狀態(tài)存在時不發(fā)生腐蝕，但當(dāng)其進(jìn)入塔頂冷凝區(qū)，系統(tǒng)溫度降低至露點溫度時，開始出現(xiàn)液態(tài)水，HCl及H2S氣體極易溶于液態(tài)水形成濃度極高的腐蝕性溶液，與金屬本體發(fā)生化學(xué)反應(yīng)從而破壞管道內(nèi)壁面。腐蝕性溶液在露點溫度處的pH值最低，腐蝕性最強(qiáng)，隨著溫度逐漸降低，液態(tài)水量不斷增加，腐蝕速率也逐漸減小。

1.3 常頂系統(tǒng)沖蝕風(fēng)險分析

常減壓裝置主要設(shè)備的操作工況如下：常壓塔T102的操作壓力為0.160 MPaG，操作溫度為140.6℃；常頂回流罐V102的操作壓力為0.127 MPaG，操作溫度為90.1℃；分離罐V103的操作壓力為0.076 MPaG，操作溫度為40.3℃；常壓注水總量為22000 kg/h。常頂瓦斯氣、常頂含硫污水、常頂一級油和常頂二級油的流量分別為4610.90、13060、133800 kg/h和74100 kg/h，氣、水、油三種介質(zhì)的工況及組成如表1-表4所示。

表1 常頂瓦斯氣組成

根據(jù)常頂系統(tǒng)的典型工況，采用工藝流程計算軟件Aspen進(jìn)行工藝過程的計算。根據(jù)常頂系統(tǒng)工藝過程分析，構(gòu)建常頂冷卻系統(tǒng)工藝計算模型，如圖2所示。模型包括：常頂瓦斯氣(CD-GAS)，含硫污水中的氨、硫、氯組分轉(zhuǎn)化為NH3、H2S、HCl的物流(NSCL-WAT)，常頂一級油(OIL-YJ)，常頂二級油(OIL-EJ)，以及一級油中的H2S(H2S-YJ)，將上述物流混合后，將該股物流的操作溫度、操作壓力調(diào)整為常壓塔塔頂?shù)牟僮鳒囟?、操作壓力，通過計算可得到常頂塔塔頂出口油氣物料的流量、組成。

表2 常頂含硫污水工況及組成

表3 常頂一級油工況及組成

表4 常頂二級油工況及組成

圖2 常頂冷卻系統(tǒng)工藝計算模型

由于常壓塔T102塔頂油氣為飽和狀態(tài)，在注水前，由于體系中不存在液態(tài)水，無需考慮HCl、H2S、NH3在水中的溶解及解離平衡，此時，可以采用PR物性方法。該方法包含所有烴和非烴組分的二元交互作用參數(shù)，對于非庫存組分和烴類虛擬組分，Aspen可自動給出HC-HC交互作用參數(shù)。PR方程的公式為：

(3)

式中：p為壓力，R為氣體常數(shù)，T為溫度，V為摩爾體積；a、b為與流體特性有關(guān)的參數(shù)，Aspen數(shù)據(jù)庫包含了大部分物質(zhì)的a、b值。

經(jīng)Aspen建模仿真計算，可獲得常頂系統(tǒng)油氣水三相的物性參數(shù)，其中常頂系統(tǒng)換熱器出口總管物性參數(shù)見表5。

表5 換熱器出口總管物性參數(shù)

2 換熱器出口管道沖蝕數(shù)值計算

2.1 出口管道結(jié)構(gòu)

換熱器出口管道由直管和多個彎管連接組成，其結(jié)構(gòu)如圖3所示。流向由左至右，其中在管道的壁面存在沖蝕失效隱患。將彎管沿流向按角度劃分坐標(biāo)，彎管入口為0°，出口為90°，彎管之間的直管段以距離坐標(biāo)為基準(zhǔn)。本文對編號為a—e的彎管及彎管間直管進(jìn)行數(shù)值分析。彎管規(guī)格為700-S40，材質(zhì)為20#碳鋼，標(biāo)號為SH3408-1996，GB/T8163-1999。其中a、b、c和d號彎管為90°彎管，e號彎管為45°彎管，曲率半徑R均為1.5倍公稱直徑D。

圖3 換熱器出口管道結(jié)構(gòu)簡圖(長度單位：mm)

2.2 控制方程

本文研究的計算流體為油氣水三相混合介質(zhì)，多相流模擬采用Mixture混合物模型。氣液相均使用不可壓條件，湍流模擬選用RNGk-ε模型，并考慮重力的影響，控制方程如下。

連續(xù)方程：

(4)

動量方程：

(5)

式中：vm、ρ分別為混合相的速度和密度，P為壓強(qiáng)，n是相數(shù)，F(xiàn)是體積力，λα是第α相的相分率，ρα是第α相的密度，μm是混合物粘度，vdr,α是第α相的漂移速度。

湍流模型：

(6)

(7)

湍流耗散方程中加入Rε項，使其在不同應(yīng)變率區(qū)域?qū)ν牧髡扯冗M(jìn)行改進(jìn)，從而使得該模型在高應(yīng)變率影響下有更高的計算精度。Rε的表達(dá)式為：

(8)

2.3 邊界條件和計算格式

根據(jù)工藝分析結(jié)果確定換熱器出口管道入口介質(zhì)的物性參數(shù)見表6?；贔luent對換熱器出口管道內(nèi)的多相流動進(jìn)行流體動力學(xué)數(shù)值模擬。根據(jù)混合物模型計算，湍流采用RNGk-ε模型，并對壁面使用標(biāo)準(zhǔn)壁面函數(shù)處理。計算模型入口采用速度入口，各相速度通過相應(yīng)的體積流量和管徑換算得到，出口采用自由流動出口。壓力-速度耦合采用PISO格式，該方法包含一個預(yù)測步和兩個修正步，可加快迭代步中的收斂速度；壓力差值采用PRESTO!格式，動量方程采用二階迎風(fēng)格式離散，體積分?jǐn)?shù)通過QUICK格式離散。

表6 換熱器出口管道入口介質(zhì)物性參數(shù)

利用ICEM軟件對出口管道進(jìn)行網(wǎng)格劃分，整體采用結(jié)構(gòu)化網(wǎng)格，且管道徑向圓截面采用O型網(wǎng)格。分別劃分網(wǎng)格數(shù)為22萬、33萬、44萬、55萬和66萬的五種網(wǎng)格，使用表6中介質(zhì)的物性參數(shù)和相應(yīng)的邊界條件計算得到不同網(wǎng)格數(shù)下管道出口流量，如圖4所示。從圖4中可以發(fā)現(xiàn)，隨著網(wǎng)格數(shù)的增加，出口流量數(shù)值逐漸趨于穩(wěn)定。本文選用55萬網(wǎng)格計算結(jié)果進(jìn)行分析，認(rèn)為此時已達(dá)到了網(wǎng)格無關(guān)性要求。

圖4 網(wǎng)格無關(guān)性驗證

2.4 計算結(jié)果分析

2.4.1 水相分率分析

圖5為換熱器出口總管流場內(nèi)的水相分率分布情況。定義與彎管外側(cè)壁面相連的直管同側(cè)壁面為外側(cè)，即圖5中表示流向的虛線箭頭所在側(cè)，與之相對位置則為內(nèi)側(cè)壁面?？梢钥吹皆趶澒躠至彎管e之間這一彎管密集分布區(qū)域內(nèi)，水相沿著流場在內(nèi)外側(cè)呈現(xiàn)較大差異的分布。因為油、氣、水三相流進(jìn)入彎管段后，流道結(jié)構(gòu)的變化導(dǎo)致流向改變，經(jīng)過拐彎處在離心力的作用下，三相混合介質(zhì)都有向管道外側(cè)流動的趨勢。但是由于油相和水相的密度遠(yuǎn)大于氣相，相對受到更強(qiáng)烈的離心作用，更容易向外側(cè)運動并聚集，因此在各彎管段均呈現(xiàn)流道外側(cè)水相分率向內(nèi)側(cè)逐漸遞減的分布趨勢。在圖5(a)中注釋位置中，水相分率的最大值集中分布在彎管b出口到彎管d進(jìn)口之間的管道的外側(cè)壁面，彎管e的彎管外側(cè)壁面。流動過程中含有腐蝕性介質(zhì)的水溶液對管道材料有腐蝕作用，流場中的水相分率較低，形成高濃度的腐蝕溶液附著在管道壁面時，將對管道壁面造成極強(qiáng)的腐蝕破壞。而水相集中分布的區(qū)域更容易受到腐蝕作用，為失效危險區(qū)域。

圖5 換熱器出口總管道水相相分率分布

下面提取水相分率數(shù)據(jù)進(jìn)一步對比分析危險區(qū)域的位置。圖6為彎管沿流動方向(彎管入口為0°)以15°為間隔提取內(nèi)外側(cè)中心壁面水相分率的變化曲線?？梢园l(fā)現(xiàn)除了a彎管在沿流向45°位置的外側(cè)壁面相分率達(dá)到最大值，其他彎管外側(cè)壁面的水相分率均沿著流動方向由彎管入口向出口逐漸增大，相應(yīng)的內(nèi)側(cè)壁面水相分率則會呈現(xiàn)減小的趨勢。各個彎管外側(cè)壁面的水相分率基本都大于內(nèi)側(cè)，水相主要在外側(cè)壁面區(qū)域聚集。但是d彎管在入口處內(nèi)側(cè)水相分率會高于外側(cè)，這是因為水相經(jīng)過前面的c彎管后流動聚集到c—d直管段的外側(cè)，進(jìn)而流動至d彎管入口的內(nèi)側(cè)，導(dǎo)致了入口處內(nèi)側(cè)較高的水相分率。綜合對比圖5中水相分布云圖和圖6中水相分率數(shù)值的變化趨勢，可以發(fā)現(xiàn)，彎管a、b、c、e的外側(cè)壁面和彎管d入口處的內(nèi)側(cè)壁面區(qū)域水相分率整體會高于其他位置，在這些區(qū)域水相集中，更容易對壁面造成腐蝕。

圖6 各彎管內(nèi)外側(cè)水相分率分布圖

圖5中還可以看到，在各彎管之間的直管段上同樣存在水相分率較大的區(qū)域(如c—d段)。分別提取各直管段內(nèi)側(cè)、外側(cè)及管內(nèi)流場中心線的水相分率，分布情況如圖7所示。由于流體經(jīng)過彎管后是沿著彎管的外側(cè)進(jìn)入直管段，因此圖7(a)—(d)中直管段的水相分布同樣呈現(xiàn)外側(cè)較高、由外至內(nèi)逐漸減少的趨勢，水相主要集中在管道的外側(cè)區(qū)域。但是由于各彎管的連接結(jié)構(gòu)特征，a—b直管段的內(nèi)側(cè)壁面與b彎管的外側(cè)相連，水相分率在流動靠近b彎管的過程中逐漸增大，最終高于外側(cè)。對比所有直管的水相分率分布，b—c和c—d直管段的相分率數(shù)值范圍整體都要高于其他位置，所以這兩個區(qū)域存在更大腐蝕風(fēng)險。

2.4.2 剪切應(yīng)力分析

腐蝕介質(zhì)的水溶液對于管道壁面的腐蝕作用，在壁面生成腐蝕產(chǎn)物保護(hù)膜，脆弱的保護(hù)膜在流體的流動剪切作用下破裂脫落，進(jìn)而在不斷的腐蝕作用下產(chǎn)生新的保護(hù)膜并脫落，最終會加速壁面材料的損失破壞。提取彎管壁面的剪切應(yīng)力分布云圖如圖8所示，從圖中可以看出，剪切應(yīng)力較大的區(qū)域分布在各個彎管的內(nèi)側(cè)壁面和各直管段的外側(cè)壁面。

(a) a—b段

(b) b—c段

(c) c—d段

(d) d—e段

圖8 換熱器出口總管道剪切應(yīng)力分布

下面對各彎管內(nèi)外側(cè)壁面的剪切應(yīng)力分布情況進(jìn)一步對比分析。各彎管內(nèi)外側(cè)剪切應(yīng)力變化曲線如圖9所示，從圖中可以看到，所有彎管的內(nèi)側(cè)壁面，相比于外側(cè)，都會受到更大的剪切應(yīng)力。但是從前面的分析可知，流場中的水相主要集中在各彎管的外側(cè)壁面區(qū)域，內(nèi)側(cè)基本被氣相占據(jù)，氣相中的腐蝕介質(zhì)難以對壁面形成腐蝕作用，單純的氣相流動剪切作用很難造成破壞作用。因此，在腐蝕作用較大的情況下，分析各彎管外側(cè)壁面剪切應(yīng)力，可以發(fā)現(xiàn)它們具有相同的變化趨勢，進(jìn)出口分布不均勻，基本在30°左右開始至出口位置逐漸增大并達(dá)到相近的剪切應(yīng)力，對于彎管d(入口內(nèi)側(cè))在入口至30°之間逐漸增大達(dá)到峰值。各彎管水相分率較大的區(qū)域剪切應(yīng)力從大到小依次為：d(入口內(nèi)側(cè))、c、e、a、b，其中d(入口內(nèi)側(cè))明顯高于c、e、a、b各外側(cè)壁面。

圖9 各彎管內(nèi)外側(cè)剪切應(yīng)力分布

圖10為各彎管間直管段上內(nèi)外側(cè)壁面剪切應(yīng)力變化曲線。相比于內(nèi)外側(cè)壁面相分率的變化曲線，剪切應(yīng)力沿著流向的變化幅度更大，各直管段基本呈現(xiàn)外側(cè)壁面剪切應(yīng)力整體大于內(nèi)側(cè)的趨勢，除了在b—c段的進(jìn)出口位置內(nèi)側(cè)剪切應(yīng)力會高于外側(cè)，是由于這一段的進(jìn)出口與b、c兩彎管的內(nèi)側(cè)壁面相連，為高速氣相的剪切作用。故各直管段剪切應(yīng)力較大區(qū)域集中在外側(cè)區(qū)域，其中a—b和c—d段整體會高于其他位置。

圖10 不同直管壁面剪切應(yīng)力

2.4.3 彎管壁面測厚驗證

根據(jù)數(shù)值模擬分析得到的沖蝕失效危險區(qū)域，對實際彎管的相應(yīng)位置進(jìn)行了現(xiàn)場測厚工作，處理后得到壁厚減薄量情況如表7所示。從表7中可以發(fā)現(xiàn)，測厚數(shù)據(jù)與計算結(jié)果有相似的變化規(guī)律，a、b、c、e號彎管外側(cè)壁面減薄嚴(yán)重，都在90°達(dá)到了最大減薄量；d號彎管內(nèi)側(cè)壁面減薄嚴(yán)重，在30°達(dá)到了最大值，與圖9中壁面剪切應(yīng)力的變化趨勢相符。

表7 彎管內(nèi)/外側(cè)壁厚減薄量 mm

3 結(jié) 論

采用Aspen工藝分析確定換熱器出口管道流場中存在腐蝕性介質(zhì)，得到油氣水三相混合介質(zhì)的物性參數(shù)和組成。通過數(shù)值模擬方法得到流場內(nèi)水相分率和壁面剪切應(yīng)力的分布情況。結(jié)果發(fā)現(xiàn)：

a) 流場中的腐蝕介質(zhì)分布于油氣水三相之中，其中溶于水后的腐蝕性溶液會使管道壁面生成腐蝕產(chǎn)物保護(hù)膜；保護(hù)膜在壁面剪切應(yīng)力的作用下快速的脫落、再生，從而加速了管道的腐蝕破壞。b) 在彎管結(jié)構(gòu)中離心力的作用下，流場中水相主要集中在管道的外側(cè)，相分率由外側(cè)壁面至內(nèi)側(cè)壁面逐漸降低。在腐蝕性溶液聚集的外側(cè)壁面，各彎管和直管段剪切應(yīng)力沿著流動方向逐漸增大。

c) 在各彎管中a、b、c、e彎管外側(cè)壁面水相集中，壁面剪切應(yīng)力沿流向30°至出口段逐漸增大至峰值；彎管d入口內(nèi)側(cè)水相集中，壁面剪切應(yīng)力沿流向由入口至30°之間逐漸增大至峰值。各直管段中c—d段外側(cè)壁面為水相分率和壁面剪切應(yīng)力最大區(qū)域。以上區(qū)域是管道壁面在腐蝕和流體剪切作用下失效的高風(fēng)險區(qū)域，并且通過對管道壁面失效危險區(qū)域的現(xiàn)場測厚驗證了數(shù)值模擬結(jié)果的準(zhǔn)確性。

[1] 朱岳麟,周健,熊常健,等.煉油設(shè)備腐蝕與防護(hù)技術(shù)新進(jìn)展[J].石油化工設(shè)備,2002,31(1):14-16.

[2] 曹東學(xué).常減壓裝置腐蝕及應(yīng)對措施[J].石油化工設(shè)備技術(shù),2005,26(4):36-39.

[3] ZHU M, SUN L, OU G, et al. Erosion corrosion failure analysis of the elbow in sour water stripper overhead condensing reflux system[J]. Engineering Failure Analysis,2016,62:93-102.

[4] 趙敏,康強(qiáng)利,馬紅杰,等.煉油廠常減壓蒸餾裝置腐蝕防護(hù)現(xiàn)狀[J].腐蝕科學(xué)與防護(hù)技術(shù),2012,24(5):430-432.

[5] 于艷秋,張景生,劉小輝,等.常減壓蒸餾裝置塔頂腐蝕案例分析與控制[J].石油化工腐蝕與防護(hù),2007,24(5):29-30.

[6] 韓紅亮.加工進(jìn)口劣質(zhì)原油常減壓腐蝕原因分析及防護(hù)措施[J].山東化工,2010,39(7):42-45.

[7] 程光旭,馬貞欽,胡海軍,等.常減壓裝置塔頂?shù)蜏叵到y(tǒng)露點腐蝕及銨鹽沉積研究[J].石油化工設(shè)備,2014,43(1):1-8.

[8] 偶國富,朱祖超,楊健,等.加氫反應(yīng)流出物空冷器系統(tǒng)的腐蝕機(jī)理[J].中國腐蝕與防護(hù)學(xué)報,2009,25(1):61-64.

[9] SHARGAY C A, MARCINIEC S. Analyzing water washing requirements for low severity hydroprocessing units[C]// Corrosion, 2004. New Orleans: NACE International,2004:04646.

[10] 金浩哲,偶國富,王寬心,等.加氫處理系統(tǒng)NH4Cl結(jié)晶沉積預(yù)測及優(yōu)化防控[J].石油學(xué)報(石油加工),2014,30(4):662-667.

[11] 金浩哲,王寬心,偶國富,等.基于Aspen Plus的加氫反應(yīng)流出物銨鹽結(jié)晶速率計算模型與分析[J].石油學(xué)報(石油加工),2016,31(6):1444-1449.

[12] DEWAN A K, VALENZUELA D P, DUBEY S T, et al. Corrosion at metal interfaces a study of corrosion rate and solution properties, including electrical conductance, viscosity, and density[J]. Industrial & Engineering Chemistry Research,2002,41(5):914-921.

[13] ZARENEZHAD B, AMINIAN A. A multi-layer feed forward neural network model for accurate prediction of flue gas sulfuric acid dew points in process industries[J]. Applied Thermal Engineering,2010,30(6):692-696.

(責(zé)任編輯： 康 鋒)

Numerical Simulation of Erosion Characteristics of Atmospheric Distillation Tower Overhead Heat Exchanger Outlet Pipe

OUGuofu,AI-barakaniA.Ahmed,XUJian,YEHaojie,JINHaozhe

(Institute of Flow Induced Corrosion, Zhejiang Sci-Tech University, Hangzhou 310018, China)

In a case study of erosion failure of overhead heat exchanger outlet pipe of atmospheric distillation tower in petrochemical industry, the distribution of corrosive medium in multiphase flow is calculated using the process simulation software Aspen and a detailed three-dimensional flow filed of heat exchanger outlet pipe is obtained by computational fluid dynamics (CFD) software Fluent. The results show that the corrosive medium existed in oil, gas and water, and the pipe wall will be corroded by corrosive aqueous solution. The generated corrosion products protective film fall off and regeneration rapidly with the effect of wall shear stress which accelerated the corrosion damage of pipe. The water phase is mainly concentrated on the outer edge of the pipe, and the water phase volume fraction is gradually reduced from outside to inside. The wall shear stress of elbows and straight pipes are gradually increased along the flow direction on the outside that corrosive solution gathered. The inside wall along the flow direction from entrance to 30° of the fourth elbow, the outside wall from 30° to exit section of the other four elbows, and the outside wall of straight pipe between the third and fourth elbow are failure high-risk areas with maximum water phase volume fraction and wall shear stress. And the results of simulation and wall thickness measurement are basic coincidence.

heat exchanger pipe; erosion characteristics; process analysis; numerical simulation

10.3969/j.issn.1673-3851.2017.07.009

2017-03-04 網(wǎng)絡(luò)出版日期： 2017-05-24

國家自然科學(xué)基金委員會-神華集團(tuán)有限公司煤炭聯(lián)合基金項目(U1361107)；浙江省公益技術(shù)應(yīng)用研究計劃項目(2015C31013)；浙江省自然科學(xué)基金項目(LY17E060008)

偶國富(1965-)，男，江蘇太倉人，教授，博士，主要從事流動腐蝕預(yù)測及特種設(shè)備安全保障技術(shù)方面的研究。

金浩哲，E-mail：haozhe2007@163.com

TQ055.8;TE986

A

1673- 3851 (2017) 04- 0518- 09