李建 崔久濤 何曉昶 夏春暉

摘 ?????要：天然氣凈化脫硫脫碳裝置胺再生產(chǎn)的酸性氣，需要進(jìn)行達(dá)標(biāo)處理并回收硫磺，硫磺回收可以很好地處理含H2S酸氣。對(duì)于本裝置產(chǎn)生的低濃度酸性氣，使用常規(guī)的克勞斯工藝難以達(dá)到硫回收的技術(shù)要求，因此采用直接氧化法。采用流體仿真軟件Fluent對(duì)反應(yīng)器進(jìn)行了流場(chǎng)仿真模擬分析，并對(duì)反應(yīng)時(shí)反應(yīng)器內(nèi)的溫度、壓力、流線等參數(shù)進(jìn)行了計(jì)算，通過(guò)對(duì)比分析反應(yīng)器不同結(jié)構(gòu)下的仿真計(jì)算結(jié)果，最終確定使用該反應(yīng)器的結(jié)構(gòu)和參數(shù)。針對(duì)該廠反應(yīng)器存在的問(wèn)題，對(duì)該廠的反應(yīng)器進(jìn)行了流場(chǎng)仿真，仿真模型涵蓋了進(jìn)氣口、觸媒反應(yīng)區(qū)域、出氣口;根據(jù)催化劑參數(shù)以及反應(yīng)數(shù)據(jù)仿真分析，計(jì)算初步算例、反應(yīng)發(fā)生程度以及放熱量。利用CFD進(jìn)行整個(gè)流場(chǎng)的仿真分析，尋找理論超溫點(diǎn)。通過(guò)優(yōu)化流場(chǎng)以及反應(yīng)器結(jié)構(gòu)設(shè)計(jì)，并對(duì)整流裝置以及觸媒筐組件設(shè)計(jì)，根據(jù)流體壓強(qiáng)改變反應(yīng)器局部滲透率，使反應(yīng)氣能平穩(wěn)的通過(guò)反應(yīng)器間隙及催化劑床層，保證反應(yīng)氣速均勻分布，保證流場(chǎng)條件滿足催化劑技術(shù)指標(biāo)。同時(shí)通過(guò)加裝混合格柵，實(shí)現(xiàn)反應(yīng)氣濃度在催化劑中均勻分布，保證硫化氫的催化氧化反應(yīng)均勻進(jìn)行。并通過(guò)混合格柵的優(yōu)化設(shè)計(jì)，使反應(yīng)氣溫度分布均勻。

關(guān) ?鍵 ?詞：硫磺回收;流場(chǎng)仿真;CFD;氧化

中圖分類(lèi)號(hào)：TQ051.1+9 ??????文獻(xiàn)標(biāo)識(shí)碼： ??????文章編號(hào)： 1671-0460（2019）08-1888-06

Abstract： Acid gas from natural gas purification desulphurization and decarbonization plant needs to be treated and sulfur needs to be recovered. Sulphur recovery is beneficial to treating H2S acid gas well. For the low concentration of acid gas produced by this device， it is difficult to meet the technical requirements for sulfur recovery using the conventional Klaus process. Therefore， direct oxidation method is always used. Fluid simulation software Fluent was used to simulate the flow field in the reactor， and the parameters of temperature， pressure and streamline in the reactor were calculated. By comparing and analyzing the simulation results of different structures of the reactor， the structure and parameters of the reactor were finally determined. In view of the problems existing in the reactor， the flow field simulation of the reactor was carried out， including air intake， catalytic reaction area and air outlet. According to the parameters of catalyst and the simulation analysis of reaction data， the preliminary examples， the degree of reaction and the heat discharge were calculated. CFD was used to simulate and analyze the whole flow field， the theoretical hyperthermia point was found. By optimizing the flow field and reactor structure design， and designing the rectifier and the catalyst basket assembly， the local permeability of the reactor was changed according to the fluid pressure， so that the reaction gas passed smoothly through the reactor gap and the catalyst bed to ensure the uniform distribution of the reaction gas speed， and ensure that flow field conditions met the technical specifications of catalyst. At the same time， the reaction gas concentration was evenly distributed in the catalyst by adding mixed grille to ensure that the catalytic oxidation reaction of hydrogen sulfide was carried out evenly. And through the optimization design of the hybrid grid， the temperature distribution of the reaction gas was uniform.

Key words： Sulphur recovery;Flow field simulation;CFD;oxidization

1 ?概 況

國(guó)內(nèi)某天然氣凈化廠脫硫脫碳裝置胺再生產(chǎn)生的酸性氣的脫硫設(shè)備開(kāi)工生產(chǎn)后，出現(xiàn)了反應(yīng)器局部溫度過(guò)高、脫硫效率遠(yuǎn)遠(yuǎn)低于設(shè)計(jì)水平，而且尾氣中二氧化硫的含量超過(guò)國(guó)家標(biāo)準(zhǔn)，造成空氣污染和硫資源的浪費(fèi)。這些影響了設(shè)備的安全生產(chǎn)和平穩(wěn)運(yùn)行，是亟需解決的問(wèn)題。

2 ?存在問(wèn)題

自硫磺回收設(shè)備投入生產(chǎn)以后，發(fā)現(xiàn)硫磺回收裝置溫度控制存在一定問(wèn)題，在局部會(huì)出現(xiàn)一定的超溫現(xiàn)象，致使反應(yīng)無(wú)法正常進(jìn)行，導(dǎo)致了恒溫反應(yīng)器出口的硫化氫含量高于預(yù)設(shè)值，超過(guò)了下游絕熱反應(yīng)器的允許進(jìn)口溫度，進(jìn)而影響整個(gè)硫磺回收系統(tǒng)的進(jìn)程，導(dǎo)致硫回收率下降嚴(yán)重。

3 ?問(wèn)題分析

3.1 ?反應(yīng)器結(jié)構(gòu)影響

混合反應(yīng)氣由裝置下方進(jìn)氣口加壓進(jìn)入反應(yīng)器，導(dǎo)致裝置下部的反應(yīng)氣速度快、壓強(qiáng)大，通過(guò)催化劑床層的速度過(guò)快，造成底部局部反應(yīng)量過(guò)大，現(xiàn)在設(shè)備的換熱效率難以滿足要求，出現(xiàn)超溫現(xiàn)象[1]。所以要在一定程度上對(duì)結(jié)構(gòu)進(jìn)行改進(jìn)，減小裝置下部反應(yīng)量，保證整個(gè)反應(yīng)器中均勻反應(yīng)。

3.2 ?反應(yīng)器預(yù)熱溫度影響

通常情況下，隨著反應(yīng)氣預(yù)熱溫度增加，硫化氫氣體催化氧化反應(yīng)的效率先增加后降低。溫度過(guò)高會(huì)造成選擇氧化催化劑失活，導(dǎo)致H2S氣體被直接氧化成SO2，放出大量熱，影響反應(yīng)使用效率。所以要在一定程度上調(diào)整反應(yīng)氣預(yù)熱溫度，尋找最佳溫度點(diǎn)。

3.3 ?反應(yīng)器組分比例影響

O2/H2S值的大小會(huì)影響生成SO2的副反應(yīng)的發(fā)生情況。通常情況下，隨著O2/H2S值的減小，主反應(yīng)發(fā)生的比例呈先上升再下降的趨勢(shì)。副反應(yīng)相較主反應(yīng)來(lái)說(shuō)會(huì)放出更多的熱量，導(dǎo)致超溫現(xiàn)象，與此同時(shí)還降低了反應(yīng)效率。所以要在一定程度上調(diào)整空氣和酸性氣的配比，保證主反應(yīng)的正常發(fā)生。

4 ?改進(jìn)措施

4.1 ?觸媒筐結(jié)構(gòu)改進(jìn)

觸媒筐是用于裝盛觸媒并且合成反應(yīng)的發(fā)生地，其設(shè)計(jì)需要考慮有效去反應(yīng)熱，保證氣體均勻的通過(guò)觸媒層，盡可能多接觸媒，結(jié)構(gòu)簡(jiǎn)單可靠。

目前觸媒筐基本都是多床層、徑向或軸徑向流的觸媒床層帶內(nèi)換熱器的內(nèi)件結(jié)構(gòu)型式[2]。徑向或軸徑向流的內(nèi)件結(jié)構(gòu)設(shè)計(jì)，即觸媒床層填充在2 個(gè)筒體之間，氣體通過(guò)觸媒床層的路程短，可以使得觸媒床的阻力降最小，更適用于大中型塔。為在一定程度上利用觸媒筐控制反應(yīng)量，從而控制反應(yīng)區(qū)域溫度，多層媒觸筐的開(kāi)孔以及孔的位置分布是關(guān)鍵。

觸媒筐壁面的進(jìn)氣口直徑D（mm）在保證總體孔隙率不變的情況下根據(jù)流場(chǎng)壓強(qiáng)大小進(jìn)行優(yōu)化設(shè)計(jì)。由于恒溫反應(yīng)器底部流體速度快壓強(qiáng)大，故減小觸媒筐底部孔徑，并隨高度增加逐步擴(kuò)大孔徑。

同時(shí)考慮流體自旋因素，觸媒筐殼體的孔延與豎直方向夾角為5°的直線進(jìn)行排布，來(lái)保證整個(gè)反應(yīng)器各部分反應(yīng)均勻發(fā)生，并在一定程度上減少超溫區(qū)域的反應(yīng)量，來(lái)保證安全裕度。

觸媒筐結(jié)構(gòu)改進(jìn)目的是在一定程度上限制裝置下部反應(yīng)的發(fā)生，是反應(yīng)器內(nèi)反應(yīng)均勻發(fā)生。觸媒筐的改進(jìn)重點(diǎn)主要集中在孔徑大小以及孔的位置排布上（圖1-2）。

4.2 ?散熱管排布優(yōu)化設(shè)計(jì)

熱管置于反應(yīng)器筒體內(nèi)，是反應(yīng)器主要換熱元件，其排布會(huì)對(duì)局部的換熱效率產(chǎn)生極大的影響，其設(shè)計(jì)需要考慮反應(yīng)器內(nèi)溫度的分布情況，并根據(jù)換熱效率的要求設(shè)計(jì)換熱管的排布情況。

換熱管在管板上的正三角形排列形式可以實(shí)現(xiàn)在同樣的管板面積上排列最多的管數(shù)，故保留原有的三角形排布設(shè)計(jì)。由于反應(yīng)混合氣進(jìn)入觸媒筐就開(kāi)始反應(yīng)，導(dǎo)致觸媒筐外圍反應(yīng)量較大，溫度較高。所以在原有的內(nèi)外圈排布設(shè)計(jì)的基礎(chǔ)上，重新進(jìn)行內(nèi)、中、外圈三層換熱排布設(shè)計(jì)。

由于反應(yīng)器強(qiáng)度要求，外圈的換熱管的分布設(shè)計(jì)保持原有設(shè)計(jì)，孔間間距62 mm，共282個(gè);中圈孔間間距80 mm，共234個(gè);內(nèi)圈空間間距100 mm，共60個(gè)。由于換熱管的分布設(shè)計(jì)未考慮加劑管、拉桿管等功能管的位置分布，故對(duì)換熱管進(jìn)行余量設(shè)計(jì)（圖3-4）。

4.3 ?爐體長(zhǎng)徑比優(yōu)化

原則上在保證酸性氣在反應(yīng)器停留時(shí)間的基礎(chǔ)上，對(duì)爐體直徑和長(zhǎng)度進(jìn)行優(yōu)化?？s短長(zhǎng)度，增加直徑，不僅提高了反應(yīng)爐的剛度，降低了整體撓度，而且對(duì)于反應(yīng)器內(nèi)氣體的流通和對(duì)內(nèi)插管取走反應(yīng)熱均有幫助[3]。尤其是長(zhǎng)度縮短，它可以減少熱應(yīng)力對(duì)于附屬固定設(shè)備的工作應(yīng)力，熱位移量也比改造前大大降低，不僅如此，因撓度的減小，對(duì)設(shè)備的吊運(yùn)安裝都有好處。

湍流出現(xiàn)在速度變動(dòng)的地方。這種波動(dòng)使得流體介質(zhì)之間相互交換動(dòng)量、能量和濃度變化，而且引起了數(shù)量的波動(dòng)。由于這種波動(dòng)是小尺度且是高頻率的，所以在實(shí)際工程計(jì)算中直接模擬的話對(duì)計(jì)算機(jī)的要求會(huì)很高。實(shí)際上瞬時(shí)控制方程可能在時(shí)間上、空間上是均勻的，或者可以人為的改變尺度，這樣修改后的方程耗費(fèi)較少的計(jì)算機(jī)。但是，修改后的方程可能包含有未知變量，湍流模型需要用已知變量來(lái)確定這些變量[6]。

在本論文中采用的是標(biāo)準(zhǔn)k-ε模型。它是從實(shí)驗(yàn)現(xiàn)象中總結(jié)出來(lái)的半經(jīng)驗(yàn)的公式。

湍流動(dòng)能方程k：

方程中是由層流速度梯度引起的湍流動(dòng)能，是由浮力引起的湍流動(dòng)能，是由在可壓縮湍流中，擴(kuò)散引起的波動(dòng)，是常量。

（1）未優(yōu)化設(shè)計(jì)工況CFD仿真

本仿真模擬未優(yōu)化情況下反應(yīng)器在設(shè)計(jì)工況下的流場(chǎng)模型，各參數(shù)設(shè)置為酸性氣入口流量13 500 Nm3/h，配風(fēng)量1 152 Nm3/h，酸性氣溫度160 ℃，H2S濃度為2.32%（體積分?jǐn)?shù)），壓力60 kPa。

① ?溫度標(biāo)量分布

由于反應(yīng)器下部混合氣速度快、壓強(qiáng)大，更容易進(jìn)入催化劑床層發(fā)生反應(yīng)，造成反應(yīng)器中下部反應(yīng)過(guò)于集中，而反應(yīng)器上部反應(yīng)量較少，導(dǎo)致溫度分布不均勻[5]，反應(yīng)器底部溫度溫度過(guò)高。CFD 最終計(jì)算結(jié)果：硫磺回收恒溫反應(yīng)器溫度標(biāo)量分布云圖，如圖5所示。其最高溫度達(dá)到283 ℃，而頂部最低溫度僅有218 ℃，溫度差達(dá)到68 ℃（圖5）。

② ?反應(yīng)器內(nèi)速度標(biāo)量云圖

由于反應(yīng)器整體結(jié)構(gòu)和催化劑床層的影響，混合氣在進(jìn)入觸媒筐組件和反應(yīng)器殼體間隙時(shí)速度增大，如圖6所示，導(dǎo)致混合氣更容易進(jìn)入催化劑床層反應(yīng)。

混合氣經(jīng)過(guò)反應(yīng)氣進(jìn)氣口、觸媒筐組件進(jìn)入催化劑床層。由于反應(yīng)器入口壓強(qiáng)和裝置壓損的影響，反應(yīng)氣集中分布在反應(yīng)器中下部。反應(yīng)器內(nèi)流線如圖8所示，可視化了從反應(yīng)器入口開(kāi)始的無(wú)質(zhì)量粒子的路徑，其大小由該點(diǎn)的流場(chǎng)速度、壓力、湍流強(qiáng)度等函數(shù)確定。流線的路徑反映出無(wú)質(zhì)量粒子在不同區(qū)域內(nèi)受到局部湍流情況的影響。煙氣速度梯度的大小以流線的顏色來(lái)區(qū)分[7]。

（2）優(yōu)化后的設(shè)計(jì)工況CFD仿真

本仿真模擬經(jīng)設(shè)計(jì)優(yōu)化后反應(yīng)器在設(shè)計(jì)工況下的流場(chǎng)模型，各參數(shù)設(shè)置為酸性氣入口流量13 500 Nm3/h ，酸性氣溫度160 ℃，H2S濃度2.32%（mol），配風(fēng)量1 152 Nm3/h，壓力60 kPa。

① 溫度標(biāo)量分布

由于觸媒筐組件的優(yōu)化設(shè)計(jì)，改變了反應(yīng)器不同位置的滲透率，避免了反應(yīng)器底部淤積的現(xiàn)象[8]，使反應(yīng)器整體的反應(yīng)分布趨于均勻。CFD 最終計(jì)算結(jié)果：硫磺回收恒溫反應(yīng)器溫度標(biāo)量分布云圖，如圖9所示。其最高溫度為253 ℃，而最低溫度為243 ℃，溫度為10 ℃，符合設(shè)計(jì)要求。

② 反應(yīng)器內(nèi)速度標(biāo)量云圖

由于對(duì)觸媒筐的優(yōu)化設(shè)計(jì)，減小了反應(yīng)器底部孔隙率來(lái)限制底部進(jìn)氣量，使底部催化劑床層氣體流速降低，如圖10所示，反應(yīng)區(qū)域較未優(yōu)化前上移，底部反應(yīng)量減少，在反應(yīng)器內(nèi)分布更加平均。

③ ?反應(yīng)器內(nèi)速度流線分布云圖

反應(yīng)器中截面速度標(biāo)量整體分布云圖，如圖11所示。混合氣經(jīng)過(guò)反應(yīng)氣進(jìn)氣口、觸媒筐組件進(jìn)入催化劑床層。由于反應(yīng)器結(jié)構(gòu)參數(shù)的優(yōu)化設(shè)計(jì)，反應(yīng)氣相較于優(yōu)化前分布均勻。

6 ?結(jié)束語(yǔ)

SOC催化氧化工藝可以完成含低濃度硫化氫酸氣的硫處理，實(shí)現(xiàn)硫磺回收，但是因其為放熱反應(yīng)，所以對(duì)于溫度控制和反應(yīng)器結(jié)構(gòu)設(shè)計(jì)要求較高[9]。

同時(shí)采用增加混合格柵，改進(jìn)觸媒框結(jié)構(gòu)和對(duì)散熱管的排布重新設(shè)計(jì)這三種改進(jìn)方案，可以解決對(duì)該廠脫硫恒溫反應(yīng)器局部溫度過(guò)熱的問(wèn)題。

參考文獻(xiàn)：

[1]傅杰.列管式固定床反應(yīng)器流體流動(dòng)、傳熱及其數(shù)值特征[D].華東理工大學(xué)，2013.

[2]吳楊.列管式固定床反應(yīng)器內(nèi)部流場(chǎng)模擬[D].天津大學(xué)，2010.

[3]Stephen Santo a，Mahin Rameshni.the challenges of designing grass root Sulphur recovery units with a wide range of H2S concentration from natural gas[J]. Journal of Naturals Gas Science and Engineering， 2014，18：137-148.

[4]鄭丹星.流體與過(guò)程熱力學(xué)[M].第二版.北京：化學(xué)工業(yè)出版社，2010.

[5]張曉華，于艷秋.特大型克勞斯反應(yīng)爐運(yùn)行問(wèn)題分析及改造[J].煉油技術(shù)與工程.2016.46（3）：40-44.

[6]于福義.現(xiàn)代汽車(chē)空調(diào)系統(tǒng)數(shù)值模擬仿真[D].重慶大學(xué)，2005.

[7]韓占中. FLUENT——流體工程仿真計(jì)算實(shí)例與應(yīng)用[M].北京：北京理工大學(xué)出版社，2010.

[8]P.He，M.liu，J.L.Luo.Stabilization of Platinum Anode Catalyst in a H2S-O2 Oxide Fuel Cell with an Intermediate TiO2 Layer[J].Journal of The Electrochemical Socity，2002，149（7）：808-814.

[9]侯倩倩.硫磺回收裝置尾氣二氧化硫達(dá)標(biāo)排放工藝優(yōu)化[J].技術(shù)研究，2018（1）：101.