黃曉媛，莫建松，吳忠標(biāo)

1.浙江大學(xué)環(huán)資學(xué)院，浙江 杭州 310058；2.浙江天藍(lán)環(huán)保技術(shù)股份有限公司，浙江 杭州 311202；3.浙江省工業(yè)鍋爐爐窯煙氣污染控制工程技術(shù)研究中心，浙江 杭州 311202

石灰石-石膏濕法煙氣脫硫系統(tǒng)的數(shù)值模擬

黃曉媛1,2,3，莫建松2,3，吳忠標(biāo)1

1.浙江大學(xué)環(huán)資學(xué)院，浙江 杭州 310058；2.浙江天藍(lán)環(huán)保技術(shù)股份有限公司，浙江 杭州 311202；3.浙江省工業(yè)鍋爐爐窯煙氣污染控制工程技術(shù)研究中心，浙江 杭州 311202

石灰石-石膏濕法脫硫工藝是目前應(yīng)用最廣泛的脫硫工藝。為了對脫硫工藝進(jìn)行優(yōu)化，本文基于雙膜理論，以某電廠2×135MW機(jī)組燃煤鍋爐的脫硫系統(tǒng)為研究對象，建立了石灰石-石膏濕法煙氣脫硫系統(tǒng)的數(shù)學(xué)模型，詳細(xì)描述了脫硫工藝中的氣液傳質(zhì)、化學(xué)反應(yīng)的過程和機(jī)理。通過將模型求解結(jié)果與工業(yè)運(yùn)行數(shù)據(jù)進(jìn)行比較，發(fā)現(xiàn)計(jì)算結(jié)果與工業(yè)數(shù)據(jù)的吻合程度較高，并考察了幾個(gè)關(guān)鍵工藝參數(shù)對脫硫效率的影響。計(jì)算結(jié)果顯示，提高液氣比或漿液pH值，脫硫率隨之提高；增大入口SO2分壓或煙氣流量，會(huì)導(dǎo)致脫硫率下降；循環(huán)漿液中的Cl-離子濃度不宜超過21 300 mg/L。本研究所建立的模型及研究結(jié)果，可以為石灰石-石膏法脫硫系統(tǒng)的設(shè)計(jì)和優(yōu)化提供參考。

石灰石-石膏法 煙氣脫硫 噴淋塔 雙膜理論 數(shù)值模擬

SO2是一種主要來自于燃煤鍋爐煙氣的大氣污染物，世界各國都有嚴(yán)格的排放標(biāo)準(zhǔn)[1]?？刂芐O2排放主要通過采用脫硫措施實(shí)現(xiàn)，常見的煙氣脫硫工藝有濕法脫硫、干法脫硫和半干法脫硫等[2]。石灰石-石膏濕法脫硫工藝具有脫硫效率高（92%～98%）、運(yùn)行成本低等優(yōu)點(diǎn)，是目前應(yīng)用最廣泛的煙氣脫硫技術(shù)[2]。

石灰石-石膏法脫硫工藝是以石灰石漿液與煙氣中的SO2接觸并發(fā)生反應(yīng)，達(dá)到脫硫的目的。這是一個(gè)典型的化學(xué)吸收過程，涉及到氣液固三相反應(yīng)過程，且化學(xué)反應(yīng)和傳質(zhì)過程同時(shí)發(fā)生，導(dǎo)致吸收過程非常復(fù)雜。多年來，關(guān)于石灰石-石膏脫硫系統(tǒng)的操作和設(shè)計(jì)主要依賴工程經(jīng)驗(yàn)，缺乏理論依據(jù)。為此，國內(nèi)外學(xué)者對脫硫過程的傳質(zhì)和反應(yīng)機(jī)理做了大量的理論研究工作[3-13]。Kill等[3,4]基于氣液吸收的雙膜理論對采用填料塔的脫硫系統(tǒng)進(jìn)行了模擬，預(yù)測了不同操作條件下的脫硫率以及HCl對脫硫效果的影響。項(xiàng)光明[5]基于雙膜理論對液柱噴射塔脫硫系統(tǒng)進(jìn)行了研究，通過計(jì)算得到了不同操作條件下脫硫塔內(nèi)SO2和漿液pH值沿吸收塔高度的分布情況。Zhong等[6,7]等用溶質(zhì)滲透理論對采用噴淋塔的脫硫過程進(jìn)行了描述，并考察了工藝參數(shù)和噴淋層排布方式對脫硫效果的影響。張華偉[8,9]用基于雙膜理論的數(shù)學(xué)模型對噴淋塔脫硫系統(tǒng)進(jìn)行靜態(tài)仿真和動(dòng)態(tài)擾動(dòng)實(shí)驗(yàn)，模擬結(jié)果與工程運(yùn)行數(shù)據(jù)和經(jīng)驗(yàn)吻合良好。Dou[10]、Kallinikos[11]、Neveux[12]等建立的脫硫數(shù)學(xué)模型均基于雙膜理論，計(jì)算結(jié)果顯示，對于噴淋吸收塔，漿液液滴的直徑對脫硫效率有重要影響，液滴越小，越有利于脫硫效率的提高。Zhu等[13]同樣基于雙膜理論建立了噴淋塔脫硫數(shù)學(xué)模型，在考察不同工藝參數(shù)對脫硫效率影響的基礎(chǔ)上，計(jì)算了不同條件下達(dá)到脫硫要求所需要的吸收區(qū)最小高度，其計(jì)算結(jié)果能夠?yàn)楣I(yè)上設(shè)計(jì)脫硫塔提供理論指導(dǎo)?？梢钥闯?，通過建立數(shù)學(xué)模型進(jìn)行數(shù)值模擬能夠獲得對石灰石-石膏法脫硫工藝更深入的認(rèn)識。然而，這些模型有的過于復(fù)雜龐大[3,6]，有的數(shù)學(xué)求解方法較難掌握[11,12]，局限了這些模型的應(yīng)用。因此，本文擬在這些研究的基礎(chǔ)上，建立一個(gè)相對簡約、便于求解的石灰石-石膏濕法脫硫數(shù)學(xué)模型對脫硫過程進(jìn)行數(shù)值模擬。從已有的研究可以發(fā)現(xiàn)，雙膜理論在描述脫硫過程時(shí)被廣泛采用，這一理論概念簡明、數(shù)學(xué)處理方便，而且對于采用不同類型吸收塔的脫硫系統(tǒng)均能獲得較好的計(jì)算結(jié)果，因此本文同樣采用雙膜理論，建立石灰石-石膏法脫硫系統(tǒng)的數(shù)學(xué)模型。通過數(shù)值模擬，除了研究幾個(gè)關(guān)鍵參數(shù)與脫硫效率之間的關(guān)系、為脫硫系統(tǒng)的運(yùn)行和脫硫工藝的優(yōu)化與設(shè)計(jì)提供參考之外，還擬對其它因素（如Cl-離子濃度）對脫硫效果的影響進(jìn)行考察。

1 脫硫系統(tǒng)

研究中模擬的石灰石-石膏濕法脫硫工藝流程如圖1所示。脫硫系統(tǒng)主要由吸收塔和氧化槽組成，吸收塔采用目前應(yīng)用最廣泛的噴淋塔。漿液循環(huán)泵將脫硫漿液輸送至噴淋層，噴淋層上分布著許多噴嘴，噴嘴將漿液分散成無數(shù)的小液滴。煙氣從吸收塔底部進(jìn)入并向上流動(dòng)，在與液滴接觸的過程中通過吸收并與石灰石反應(yīng)脫去SO2、HCl、HF等有害氣體，最后經(jīng)過除霧器除去水汽后進(jìn)入煙囪排放。在氧化槽底部通過鼓入空氣將脫硫產(chǎn)物氧化為石膏，石膏結(jié)晶后排出送去脫水得到合格的脫硫石膏。

2 數(shù)學(xué)模型

2.1 基本假設(shè)

由于脫硫塔內(nèi)的傳質(zhì)和反應(yīng)過程很復(fù)雜，很難對每一個(gè)過程都建立精確的數(shù)學(xué)模型，為了能夠相對容易、但是又比較準(zhǔn)確的描述脫硫過程，在建模過程中對脫硫系統(tǒng)做了一些基本假設(shè)。

（1）氣體為理想氣體；

（2）在石灰石漿液下落過程中，液滴被看作是嚴(yán)格的球形，且運(yùn)動(dòng)軌跡為垂直的直線；

（3）SO2和CO2進(jìn)入水中后以水合狀態(tài)存在；

（4）漿液里的正負(fù)電荷和化學(xué)反應(yīng)瞬間實(shí)現(xiàn)平衡；

圖1 石灰石-石膏法脫硫工藝Fig.1 Flow diagram of wet limestone-gypsum desulfurization system

（5）氣相流場在徑向上是均勻的，只考慮在軸向上的濃度變化；

（6）氧化槽內(nèi)混合充分均勻，槽內(nèi)各處的成分一致；

（7）忽略煙氣流經(jīng)氧化槽上方時(shí)的吸收量。

2.2 化學(xué)反應(yīng)速率

2.2.1 CaCO3的溶解

石灰石經(jīng)研磨、加水制成漿液，不斷加入SO2吸收塔內(nèi)。石灰石中的CaCO3不斷溶解，生成的CO32-解離后產(chǎn)生OH-，從而穩(wěn)定住漿液的pH值，反應(yīng)式如式(1)所示：

CaCO3的溶解度很小，且其溶解速度受傳質(zhì)過程控制，是脫硫反應(yīng)的關(guān)鍵步驟。石灰石溶解速率被認(rèn)為是其相對飽和度的函數(shù)，如式(4)所示[5,6,9,11]：

2.2.2 SO2的吸收

氣相SO2被水吸收并經(jīng)以下水合反應(yīng)解離：

SO2溶于水產(chǎn)生的H+與CO32-解離后產(chǎn)生OH-中和生成水，促使反應(yīng)平衡不斷向右移動(dòng)。

SO2的化學(xué)吸收反應(yīng)速度很快，而且反應(yīng)物濃度低，因此反應(yīng)主要發(fā)生在液膜內(nèi)。此時(shí)，SO2的吸收速率可用方程(7)描述[5,9,11]：

aA為吸收塔中每單位體積的氣液接觸面積，為漿液液滴平均直徑的函數(shù)：

SO2吸收的增強(qiáng)因子ESO2可用方程(9)表示[11]：

式中，[ST]為漿液中所有四價(jià)硫組分的總濃度，即SO2(aq)、HSO3-和SO32-的濃度之和。

在SO2吸收的同時(shí)，石灰石與H+反應(yīng)后產(chǎn)生的CO2氣體要從漿液中釋放出來，其解吸速率為：

CO2解吸的增強(qiáng)因子與ESO2的形式一致：

式中，[CT]為漿液中所有碳組分的總濃度，即CO2(aq)、HCO3-和CO32-的濃度之和。

2.2.3 HSO3-的氧化

SO2被吸收后解離產(chǎn)生的HSO3-和SO32-離子在氧化空氣的作用下，被氧化成SO42-：

目前能被廣泛接受的在pH值為4.5~5.5內(nèi)（此時(shí)SO32-的存在可以忽略）、Mn2+催化的HSO3-氧化速率為[11]：

吸收段中氧氣的吸收速率用下面的公式描述[3]：

氧化槽內(nèi)，由于空氣被大量鼓入并且通過攪拌來強(qiáng)制進(jìn)行氧化，因此認(rèn)為漿液中的氧含量與空氣中的氧氣達(dá)到了平衡。

2.2.4 石膏的結(jié)晶

HSO3-氧化生成的SO42-與Ca2+結(jié)合生成石膏CaSO4，硫酸鈣結(jié)晶后形成最后的脫硫產(chǎn)物石膏（CaSO4·2H2O）：

石膏的結(jié)晶速率與晶型無關(guān)，晶體的長大以過飽和度為推動(dòng)力，可用下面的經(jīng)驗(yàn)公式表示[3,9,11]：

2.3 質(zhì)量守恒方程

2.3.1 吸收段

在吸收塔中，沿吸收塔的任意長度微元dz對氣相及液相組分做物料平衡[11]。

氣相中SO2的質(zhì)量守恒：

液相中四價(jià)硫的質(zhì)量守恒：

液相中Ca2+的質(zhì)量守恒：

固體物質(zhì)CaCO3、CaSO4的質(zhì)量守恒：

2.3.2 氧化槽

在氧化槽中，通過攪拌器的攪拌，脫硫后的漿液與鼓入的氧化空氣充分接觸，因此可以把氧化槽近似看做一個(gè)全混流反應(yīng)器。

根據(jù)質(zhì)量守恒定律，氧化槽內(nèi)的各種離子和分子組分i應(yīng)當(dāng)滿足以下關(guān)系[5,8]：

累積量=輸入量-排出量+凈產(chǎn)生量

根據(jù)假設(shè)，不考慮漿液在吸收段的蒸發(fā)，因此循環(huán)漿液的流量沒有損失；補(bǔ)充的新鮮漿液直接加入氧化槽內(nèi)，且穩(wěn)定運(yùn)行時(shí)氧化槽內(nèi)漿液量不變，因此漿液流入量與流出量相等，即：Qfeed+Qfall=Qrecycle+Qout；整個(gè)脫硫系統(tǒng)為連續(xù)操作，供漿量和排出量一致，即有：Qfeed=Qout；前面已經(jīng)提到，氧化槽可以看作是全混流反應(yīng)器，因此各組分在氧化槽內(nèi)各處分布均勻，因此槽內(nèi)的組分濃度與排出漿液的組分濃度[i]out一樣。

2.3.3 解離平衡與電荷守恒

漿液中的四價(jià)硫物種、四價(jià)碳物種以及水各自都存在解離平衡：

漿液中陰陽離子的電荷守恒方程如下：

通過聯(lián)立求解解離平衡方程和電荷守恒方程，可以求得漿液的pH值。

2.4 脫硫效率

脫硫效率（η）為吸收塔出入口煙氣中所含SO2量的差值與進(jìn)入吸收塔的SO2量的比值。由于SO2吸收的同時(shí)有CO2放出，煙氣壓力基本保持不變，因此脫硫效率可以簡化為出入口SO2分壓的差值與入口SO2分壓之比，即：

3 數(shù)值模擬結(jié)果及分析

模型計(jì)算所需要的吸收塔結(jié)構(gòu)設(shè)計(jì)參數(shù)和運(yùn)行參數(shù)參考了某電廠2×135 MW機(jī)組的濕法脫硫系統(tǒng)，涉及到的理化參數(shù)、反應(yīng)速率常數(shù)等采用50 ℃時(shí)的值。脫硫數(shù)學(xué)模型在gPROMS流程模擬軟件中采用聯(lián)立方程法進(jìn)行求解。為了驗(yàn)證所建立模型的準(zhǔn)確性，將計(jì)算得到的結(jié)果與工程項(xiàng)目運(yùn)行工況下的數(shù)據(jù)進(jìn)行對比。

3.1 主要運(yùn)行指標(biāo)的對比

在設(shè)計(jì)工況下，脫硫系統(tǒng)達(dá)到穩(wěn)態(tài)運(yùn)行時(shí)主要指標(biāo)的設(shè)計(jì)值與模型計(jì)算結(jié)果的比較如表1所示。從表中可以看出，兩者吻合程度良好。

表1 主要性能參數(shù)設(shè)計(jì)值與模擬計(jì)算值的比較Table 1 Comparison of design value and simulation calculations of main performance parameters

3.2 關(guān)鍵參數(shù)對脫硫效果的影響

脫硫系統(tǒng)運(yùn)行時(shí)影響脫硫效率的工藝參數(shù)有很多，其中液氣比、漿液pH值、入口SO2分壓、煙氣流量的影響尤為關(guān)鍵。此外，漿液中Cl-濃度也對脫硫效果有重要影響。為了考察這些參數(shù)對脫硫效果的影響，我們考察了設(shè)計(jì)工況下，改變其中一個(gè)工藝參數(shù)之后脫硫效率的變化趨勢。設(shè)計(jì)工況下的關(guān)鍵工藝參數(shù)如表2所示。

表2 設(shè)計(jì)工況下的關(guān)鍵工藝參數(shù)值Table 2 Values of key parameters under design-condition

3.2.1 液氣比的影響

如圖2所示，在設(shè)計(jì)工況下，當(dāng)其它條件保持不變時(shí)，脫硫效率η隨著液氣比L/G的增加而增加，但是增加的幅度逐漸變小，并趨于穩(wěn)定。在入口煙氣中SO2分壓為80 Pa（2256 mg/Nm3）的情況下，采用13.4的液氣比可以使出口的SO2濃度降低至現(xiàn)行火電廠大氣污染物排放標(biāo)準(zhǔn)規(guī)定的100 mg/Nm3以下，脫硫率為96.05%，此時(shí)出口SO2濃度約為89.1 mg/Nm3。如要使出口SO2濃度下降至超低排放所要求的35.0 mg/Nm3，需要將液氣比提高至15以上，此時(shí)脫硫率要達(dá)到98.67%。因此，為了滿足環(huán)保排放要求，液氣比最少應(yīng)取13.4。

圖2 脫硫效率η隨液氣比的變化趨勢Fig.2 The relation between desulfurization efficiency and liquid-gas ratio

圖3 脫硫效率η隨漿液pH值的變化趨勢Fig.3 The relation between desulfurization efficiency and pH of slurry

3.2.2 pH的影響

漿液pH作為脫硫工藝中的一個(gè)重要因素，對脫硫效果有非常重大的影響。本工作模擬了設(shè)計(jì)工況下，脫硫效率隨著漿液pH值的變化趨勢。如圖3所示，脫硫率隨著漿液pH值的升高而不斷提高。根據(jù)脫硫過程的化學(xué)反應(yīng)機(jī)理可以推斷，高的pH值有利于SO2的吸收，而低的pH值有利于CaCO3溶出和HSO3-氧化。為了平衡幾個(gè)關(guān)鍵反應(yīng)的反應(yīng)速率，脫硫漿液的pH不宜過高也不能過低，工程應(yīng)用上一般將脫硫漿液的pH值控制為5.0～5.5。

3.2.3 入口SO2分壓的影響

設(shè)計(jì)工況下，脫硫效率隨著入口煙氣中的SO2分壓的變化趨勢如圖4所示。隨著煙氣中SO2分壓的增加，脫硫率不斷下降。在漿液噴淋量及新鮮石灰石漿液供漿量保持不變的情況下，煙氣中的SO2濃度增加，造成系統(tǒng)Ca/S比降低，漿液pH值降低，吸收SO2的速度和能力都會(huì)下降，從而導(dǎo)致脫硫效果變差。

圖4 脫硫效率η隨入口煙氣中SO2分壓pSO2的變化趨勢Fig.4 The relation between desulfurization efficiency and inlet SO2partial pressure

圖5 脫硫效率η隨煙氣流量F的變化趨勢Fig.5 The relation between desulfurization efficiency and flue gas flow rate efficiency and inlet SO2partial pressure

3.2.4 入口煙氣流量的影響

設(shè)計(jì)工況下，采用13.4的液氣比時(shí)脫硫效率隨著煙氣流量的變化趨勢如圖5所示。如果噴淋量保持不變，隨著煙氣流量增加，脫硫效率會(huì)逐漸降低。這是因?yàn)闊煔饬吭黾訉?dǎo)致液氣比和Ca/S比降低，漿液pH降低，從而導(dǎo)致脫硫率下降。

3.2.5 氯離子濃度的影響

雖然本研究在建立脫硫數(shù)學(xué)模型時(shí)未考慮HCl吸收對脫硫效果的影響，但是仍然認(rèn)為循環(huán)漿液中存在一定量的Cl-離子。脫硫系統(tǒng)在運(yùn)行時(shí)，會(huì)將部分石膏脫水產(chǎn)生的水送回脫硫系統(tǒng)用于制作新鮮石灰石漿液，這會(huì)導(dǎo)致Cl-離子回到脫硫系統(tǒng)內(nèi)并不斷累積，其濃度越來越高。當(dāng)Cl-離子的濃度達(dá)到一定的值，通過排放廢水、補(bǔ)充新鮮水的操作，使進(jìn)出脫硫系統(tǒng)的Cl-離子達(dá)到平衡，從而維持漿液中Cl-離子濃度的穩(wěn)定。過高的Cl-離子濃度對脫硫效果是不利的，目前的研究關(guān)于Cl-離子對脫硫吸收塔液位、石膏品質(zhì)、設(shè)備腐蝕等方面的討論比較多，而且廣泛認(rèn)為氯離子濃度過高會(huì)導(dǎo)致漿液pH值下降從而影響脫硫效果[14,15]，但是究竟影響程度如何，目前相關(guān)研究較少。在設(shè)計(jì)工況下，對漿液中氯離子濃度對脫硫效率的影響進(jìn)行了考察。從圖6中可以看出，隨著氯離子濃度增加，脫硫率先有所增加。這是因?yàn)槁入x子在漿液中是以氯化鈣（CaCl2）的形式存在，因此鈣離子濃度也將隨著氯離子濃度增加而增加。漿液中鈣離子增多，意味著脫硫塔內(nèi)實(shí)際的Ca/S變大，因而促進(jìn)了SO2的吸收。而隨著氯離子濃度繼續(xù)增加，漿液pH值下降，抵消了Ca/S變大的促進(jìn)作用，脫硫效率增加越來越小，并在氯離子濃度為600 mol/m3（21 300 mg/L）時(shí)達(dá)到最大值，然后開始緩慢下降。目前在工業(yè)運(yùn)行中氯離子濃度上限值控制在20 000 mg/L左右，這與模型預(yù)測的結(jié)果基本一致。

圖6 脫硫效率η隨CaCl2中的氯離子濃度[Cl-]的變化趨勢Fig.6 The relation between desulfurization efficiency and concentration of Cl-

4 結(jié) 論

分析了石灰石-石膏濕法煙氣脫硫工藝中的化學(xué)反應(yīng)和傳質(zhì)過程之后，基于雙膜理論對脫硫過程的四個(gè)關(guān)鍵反應(yīng)步驟進(jìn)行了描述，并根據(jù)物質(zhì)守恒原理對吸收段、氧化槽兩個(gè)反應(yīng)器建立了數(shù)學(xué)模型，從而得到了完整的脫硫數(shù)學(xué)模型。通過求解數(shù)學(xué)模型，得到了如下結(jié)論：

a）通過將穩(wěn)定運(yùn)行時(shí)的關(guān)鍵工藝指標(biāo)與設(shè)計(jì)指標(biāo)進(jìn)行對比，發(fā)現(xiàn)計(jì)算得到的脫硫率、漿液pH值、石灰石消耗量和石膏產(chǎn)量四個(gè)關(guān)鍵工藝指標(biāo)與設(shè)計(jì)指標(biāo)吻合良好，因此所建立的數(shù)學(xué)模型能夠合理描述脫硫過程，其計(jì)算結(jié)果是可靠的。

b）關(guān)鍵的工藝參數(shù)（液氣比、漿液pH值、入口SO2分壓、煙氣流量、漿液Cl-離子濃度）對脫硫效率有重要的影響。液氣比增加，脫硫效率逐漸提高并趨于穩(wěn)定；提高漿液pH值，有利于SO2吸收但不利于石灰石溶解，應(yīng)控制在合理范圍；增大入口SO2分壓或者煙氣流量，都會(huì)導(dǎo)致脫硫效率下降；漿液中Cl-離子的濃度不宜超過21 300 mg/L。這些計(jì)算結(jié)果，可以為脫硫系統(tǒng)工況發(fā)生變化時(shí)的操作提供指導(dǎo)。

本工作對石灰石-石膏濕法煙氣脫硫系統(tǒng)實(shí)現(xiàn)了全流程的仿真模擬，除本研究所模擬的對象外，所建立的數(shù)學(xué)模型還可以擴(kuò)展至其他脫硫工藝流程，且根據(jù)計(jì)算所取得的結(jié)果，可以對現(xiàn)有工藝進(jìn)行優(yōu)化，也可以作為工藝設(shè)計(jì)的理論依據(jù)。

[1] 燕中凱, 劉 媛, 岳 濤, 等. 我國煙氣脫硫工藝選擇及技術(shù)發(fā)展展望 [J]. 環(huán)境工程, 2013, 31(6):58-61. Yan Zhongkai, Liu Yuan, Yu Tao, et al. Selection of flue gas desulfurization process and prospects of desulfurization technology development [J]. Environmental Engineering, 2013, 31(6):58-61.

[2] Cordoba P. Status of flue gas desulphurization (FGD) systems from coal-fired power plants: Overview of the physic-chemical control process of wet limestone FGDs [J]. Fuel, 2015, 144:274-286.

[3] Kill S, Michelsen M L, Dam-Johansen K. Experimental investigation and modeling of a wet flue gas desulfurization pilot plant [J]. Industrial & Engineering Chemical research, 1998, 37(7):2792-2806.

[4] Kill S, Nygaard H, Johnsson J E. Simulation studies of the influence of HCl absorption on the performance of a wet flue gas desulphurization pilot plant [J]. Chemical Engineering Science, 2002, 57(3):347-354.

[5] 項(xiàng)光明. 液柱噴射煙氣脫硫研究 [D]. 北京: 清華大學(xué), 2003.

[6] Zhong Y, Gao X, Huo W, et al. A model for performance optimization of wet flue gas desulfurization systems of power plants [J]. Fuel Processing Technology, 2008, 89(11):1025-1032.

[7] 鐘 毅. 基于WFGD系統(tǒng)的硫、氮、汞污染物協(xié)同脫除的理論與實(shí)驗(yàn)研究 [D]. 浙江大學(xué), 2008.

[8] 何祖威, 張華偉, 展錦程, 等. 660MW機(jī)組濕法煙氣脫硫噴淋塔的簡化模型及仿真驗(yàn)證 [J]. 重慶大學(xué)學(xué)報(bào), 2009, 32(10):1169-1173. He Zuwei, Zhang Huawei, Zhan Jincheng, et al. Simplifying and validating the model of wet flue gas desulfurization spray tower for 660MW units [J]. Journal of Chongqing University, 2009, 32(10):1169-1173.

[9] 張華偉. 660MW機(jī)組石灰濕式煙氣脫硫系統(tǒng)仿真研究 [D]. 重慶: 重慶大學(xué), 2009.

[10] Dou B, Pan W, Jin Q, et al. Prediction of SO2removal efficiency foe wet flue gas desulfurization [J]. Energy Conversion and Management, 2009, 50(10): 2547-2553.

[11] Kallinikos L E, Farsari E I, Spartinos D N, et al. Simulation of the operation of an industrial wet flue gas desulfurization system [J]. Fuel Processing Technology, 2010, 91(12):1794-1802.

[12] Neveux T, Le Moullec Y. Wet industrial flue gas desulfurization unit: model development and validation on industrial data [J]. Industrial & Engineering Chemistry Research, 2011, 50(12):7579-7592.

[13] Zhu J, Ye S C, Bai J, et al. A concise algorithm for calculating absorption height in spray tower for wet limestone-gypsum flue gas desulfurization [J]. Fuel Processing Technology, 2015, 129:15-23.

[14] 張昭文, 靖汪建, 申利芬. 氯離子濃度對濕法煙氣脫硫工藝的影響 [J]. 工業(yè)安全與環(huán)保, 2009, 35(11): 28-29. Zhang Zhaowen, Jing Wangjian, Shen Lifen. Effect of chlorine ion concentration on wet limestone gypsum desulfurization technology [J]. Industrial Safety and Environmental Protection, 2009, 35(11):28-29.

[15] 靳曉潔. 石灰石-石膏濕法脫硫吸收塔中氯離子問題的探討 [J]. 電力科技與環(huán)保, 2013, 29(1): 46-47. Jin Xiaojie. Research on excessive Cl-of absorption tower on limestone-gypsum desulfurization [J]. Electric Power Technology and Environmental Protectionn, 2013, 29(1):46-47.

Abstract：Wet limestone-gypsum flue gas desulfurization (FGD) technology is currently the most widely used process. In order to optimize the wet limestone-gypsum FGD system, a coal-fired boiler of 2×135MW units of a thermal power plant was taken as the research object and the mathematical model of wet limestone-gypsum FGD system was established based on the theory of double membrane. The process and the mechanism of the gas-liquid mass transfer and chemical reaction in desulfurization process were described in detail. By comparison, it was found that the calculation results were in good accordance with the industrial operation data, and the effects of several key factors on the desulfurization efficiency were investigated. The calculated results indicated that the desulfurization rate increased with higher liquid/gas ratio or the pH value of the slurry. Increasing the SO2partial pressure or flue gas flow at the inlet would lead to the decrease of the desulfurization rate. The concentration of Cl-in the circulating slurry should not exceed 21 300 mg/L.

Numerical Simulation of Wet Limestone-Gypsum Flue Gas Desulfurization System

Huang Xiaoyuan1,2,3, Mo Jiansong2,3, Wu Zhongbiao1
1. College of Environmental & Resource Science, Zhejiang University, Hangzhou 310058, China; 2. Zhejiang Tianlan Environmental Protection Technology Limited Company, Hangzhou 311202, China; 3. Zhejiang Provincial Engineering Research Center of Industrial Boiler & Furnace Flue Gas Pollution Control, Hangzhou 311202, China

wet limestone-gypsum; flue gas desulfurization; spray tower; dual-film theory; numerical simulation

X701.3; TQ019

A

1001—7631 ( 2017 ) 01—0055—10

10.11730/j.issn.1001-7631.2017.01.0055.10

2016-12-22;

2017-01-31。

黃曉媛（1986—），女，博士研究生；莫建松（1977—），男，高級工程師，通訊聯(lián)系人。E-mail:mojs1977@163.com。