李建 崔久濤 何曉昶 夏春暉

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

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

中圖分類號：TQ051.1+9 ??????文獻(xiàn)標(biā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 ?概 況

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

2 ?存在問題

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

3 ?問題分析

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

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

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

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

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

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

4 ?改進(jìn)措施

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

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

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

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

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

觸媒筐結(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è)計

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

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

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

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

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

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

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

湍流動能方程k：

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

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

本仿真模擬未優(yōu)化情況下反應(yīng)器在設(shè)計工況下的流場模型，各參數(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)過于集中，而反應(yīng)器上部反應(yīng)量較少，導(dǎo)致溫度分布不均勻[5]，反應(yīng)器底部溫度溫度過高。CFD 最終計算結(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)器殼體間隙時速度增大，如圖6所示，導(dǎo)致混合氣更容易進(jìn)入催化劑床層反應(yīng)。

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

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

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

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

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

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

由于對觸媒筐的優(yōu)化設(shè)計，減小了反應(yīng)器底部孔隙率來限制底部進(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)過反應(yīng)氣進(jìn)氣口、觸媒筐組件進(jìn)入催化劑床層。由于反應(yīng)器結(jié)構(gòu)參數(shù)的優(yōu)化設(shè)計，反應(yīng)氣相較于優(yōu)化前分布均勻。

6 ?結(jié)束語

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

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

參考文獻(xiàn)：

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

[2]吳楊.列管式固定床反應(yīng)器內(nèi)部流場模擬[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]鄭丹星.流體與過程熱力學(xué)[M].第二版.北京：化學(xué)工業(yè)出版社，2010.

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

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

[7]韓占中. FLUENT——流體工程仿真計算實(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.