蕭嘉繁,劉鶴欣,譚厚章,杜勇樂,王毅斌

SO2-NO2-NH3-H2O四元反應(yīng)體系中氣溶膠的生成特性

蕭嘉繁,劉鶴欣,譚厚章*,杜勇樂,王毅斌

(西安交通大學(xué),熱流科學(xué)與工程教育部重點實驗室,陜西 西安 710049)

通過煙霧箱實驗,研究了SO2-NO2-NH3-H2O四元反應(yīng)體系在氣-粒轉(zhuǎn)化過程中新形成顆粒物數(shù)濃度與粒徑分布的變化.研究發(fā)現(xiàn), SO2-NO2-NH3-H2O四元反應(yīng)體系具有顯著的成核能力,且其成核強(qiáng)度大,持續(xù)時間短.當(dāng)SO2、NO2濃度為200mg/m3,NH3濃度為12′10-6時,反應(yīng)體系的氣溶膠總個數(shù)濃度在2min時達(dá)到峰值2.5′106cm-3.缺少任一種氣體均會使氣溶膠成核強(qiáng)度下降.SO2、NO2及NH3濃度分別為0的工況下氣溶膠總個數(shù)濃度峰值分別下降了41.0%、83.6%及98.5%.在電廠污染氣體排放濃度區(qū)間內(nèi),NO2對氣溶膠生成影響大于SO2.燃煤電廠控制NO排放濃度對改善煙囪出口氣溶膠數(shù)濃度更有效.在實驗基礎(chǔ)上,對顆粒物成核特性進(jìn)行擬合,反應(yīng)體系在氣-粒轉(zhuǎn)化過程中產(chǎn)生的新顆粒物總個數(shù)濃度及中值粒徑與氣態(tài)前體物濃度線性相關(guān);采用布朗團(tuán)聚模型對氣溶膠成核后的團(tuán)聚過程總個數(shù)濃度及粒徑變化進(jìn)行模擬計算,根據(jù)給定的燃煤電廠SO2、NO2、NH3排放濃度,給出預(yù)測氣溶膠顆粒成核速率、粒徑分布及總個數(shù)濃度變化的方法.

燃煤電廠；灰霾；氣溶膠；煙霧箱；成核；團(tuán)聚

近年來,頻繁爆發(fā)的灰霾天氣逐漸引起了人們的廣泛關(guān)注.灰霾的本質(zhì)是細(xì)粒子氣溶膠污染.灰霾的產(chǎn)生與氣溶膠的生成及增長密切相關(guān),可歸因于高強(qiáng)度的污染物排放,不利的氣象條件以及大量二次氣溶膠的產(chǎn)生.二次無機(jī)氣溶膠的形成,尤其是硫酸鹽、硝酸鹽和銨鹽(SNA)的爆發(fā)式增長是導(dǎo)致灰霾形成的重要因素之一[1].針對我國京津冀地區(qū)[2-4]、長三角地區(qū)[5-7]及珠三角地區(qū)[8-10]嚴(yán)重的灰霾現(xiàn)象,大量研究都表明了灰霾期SNA是微細(xì)顆粒物(PM2.5)的主要組成成分,其濃度變化與PM2.5質(zhì)量濃度變化顯著相關(guān).

高濃度氣態(tài)前體物的氣-粒轉(zhuǎn)化是SNA的主要生成途徑[11].氣態(tài)SO2經(jīng)過氣相氧化或液相反應(yīng)可轉(zhuǎn)化為硫酸鹽顆粒,NO2可轉(zhuǎn)化為硝酸鹽,NH3經(jīng)過反應(yīng)則轉(zhuǎn)化為銨鹽.因此,控制氣態(tài)前體物SO2、NO2、NH3的排放是減少SNA形成的根本途徑.而燃煤電廠廢氣中大量SO2、NO2是SNA氣態(tài)前體物的重要來源[12].統(tǒng)計結(jié)果表明[13],2010年長三角電廠排放的SO2及NO分別達(dá)到826.8,137.3Gg,分別占長三角地區(qū)人為源總排放量的34%及14%.同時,為了實現(xiàn)NO超低排放(<50mg/Nm3)的要求,大量火電機(jī)組均采用SCR脫硝技術(shù),SCR系統(tǒng)不可避免地存在氨逃逸問題.NH3對SO2向SO42-的轉(zhuǎn)化以及NO2向NO3-的轉(zhuǎn)化具有顯著的促進(jìn)作用[14].在脫硫塔之后的煙道中,由于廢氣尚未與環(huán)境空氣混合、稀釋,燃煤廢氣中SO2、NO2、NH3等污染氣體濃度較高.在高濕度環(huán)境下SO2、NO2、NH3及H2O極易在煙道中發(fā)生物理化學(xué)反應(yīng),經(jīng)過成核、團(tuán)聚產(chǎn)生大量SNA氣溶膠顆粒.SO2、NO2、NH3作為SNA的氣態(tài)前體物在燃煤電廠排放廢氣中大量存在,勢必對當(dāng)?shù)鼗姻蔡鞖獾某霈F(xiàn)產(chǎn)生一定的促進(jìn)作用.研究SO2、NO2、NH3在大氣中的氣-粒轉(zhuǎn)化過程,探究不同濃度氣態(tài)污染物在一定條件下發(fā)生化學(xué)反應(yīng)形成氣溶膠顆粒的過程中,氣溶膠顆粒的成核、團(tuán)聚特性可以明確電廠排放污染物對細(xì)粒子生成和成霾的貢獻(xiàn).

本研究通過自主搭建室內(nèi)煙霧箱,模擬燃煤電廠排放的高濃度SO2、NO2及NH3在大氣中發(fā)生化學(xué)反應(yīng)的氣-粒轉(zhuǎn)化過程,探究了不同氣氛下顆粒物的成核強(qiáng)度變化及之后的粒徑長大規(guī)律.在實驗基礎(chǔ)上,將生成的新顆粒物粒徑分布特征與氣態(tài)前體物濃度的變化關(guān)系進(jìn)行擬合,并采用布朗團(tuán)聚模型對成核后氣溶膠顆粒團(tuán)聚過程的粒徑、數(shù)濃度變化進(jìn)行計算,根據(jù)給定的燃煤電廠SO2、NO2、NH3排放濃度,建立預(yù)測SNA氣溶膠成核強(qiáng)度、粒徑及數(shù)濃度變化的方法,為明晰燃煤煙氣排放的氣相污染物向顆粒相轉(zhuǎn)化行為提供參考.

1 材料與方法

1.1 實驗裝置

為了研究SO2、NO2及NH3在大氣中發(fā)生化學(xué)反應(yīng)的氣-粒轉(zhuǎn)化過程及不同氣氛下氣溶膠顆粒物的成核強(qiáng)度變化和之后的粒徑長大規(guī)律,實驗中搭建了如圖1所示的煙霧箱實驗系統(tǒng).煙霧箱長寬均為60cm,高100cm,采用壁厚3mm的有機(jī)玻璃制成.煙霧箱內(nèi)壁貼附有0.05mm厚的聚四氟乙烯薄膜,聚四氟乙烯材料具有優(yōu)良的化學(xué)穩(wěn)定性、耐腐蝕性、高潤滑不粘性和抗老化能力.煙霧箱上下端蓋與側(cè)壁采用法蘭連接,便于定期拆卸進(jìn)行煙霧箱內(nèi)部清潔.煙霧箱上端蓋設(shè)置進(jìn)氣孔5個,下端蓋中心有聚四氟乙烯取樣管可將煙霧箱內(nèi)氣體抽取至TSI 3938L72型掃描電遷移率粒徑譜儀(SMPS)中進(jìn)行氣溶膠粒徑分布測量.煙霧箱連接CA7201型智能數(shù)顯溫濕度控制器、KELLER通用型數(shù)顯壓力表,用來實時監(jiān)測煙霧箱內(nèi)溫濕度及壓力.煙霧箱外設(shè)置有遮光罩以防止雙向光泄漏,遮光罩內(nèi)壁等距布置美國GE生產(chǎn)的F40T8BL型黑光燈8支,UV-A長波紫外線輻射波長365nm,作為實驗中光化學(xué)反應(yīng)的穩(wěn)定模擬光源.遮光罩內(nèi)壁平整敷設(shè)鋁箔,可反射紫外光,提高光利用率.

圖1 煙霧箱實驗系統(tǒng)

實驗用NO2純度99.99%,SO2為5.07%體積分?jǐn)?shù)的SO2標(biāo)準(zhǔn)氣體(高純N2為平衡氣),NH3純度99.9999%.實驗采用靜態(tài)配氣法,為了便于進(jìn)樣,將鋼瓶內(nèi)的氣體先收集于5L的聚四氟乙烯氣袋中.使用1和10mL氣密性進(jìn)樣針配合進(jìn)樣,將NO2、SO2及NH3定量注射入煙霧箱中.煙霧箱中的實驗氣體濃度通過引入反應(yīng)體系的實驗氣體體積以及煙霧箱體積計算得出.

空氣壓縮機(jī)鼓出的空氣經(jīng)多級除塵除水后,部分凈化空氣經(jīng)旁路進(jìn)入鼓泡式水汽發(fā)生器,將含水汽的空氣帶入煙霧箱中,其余凈化空氣則直接通入煙霧箱進(jìn)行沖洗.采用浮子流量計控制壓縮空氣總流量30L/min,使用SMPS檢測煙霧箱內(nèi)未通入SO2、NO2、NH3時煙霧箱中的顆粒物總個數(shù)濃度N及中值粒徑50如圖2所示.空白組中顆粒物數(shù)濃度穩(wěn)定在300~500cm-3.實驗中由氣-粒轉(zhuǎn)化而來的氣溶膠總個數(shù)濃度高達(dá)106cm-3,遠(yuǎn)高于壓縮空氣中殘余的顆粒物數(shù)濃度,因此可以將空氣壓縮機(jī)產(chǎn)生并過濾、干燥后的空氣視作凈化空氣.實驗過程中,煙霧箱內(nèi)溫度穩(wěn)定在(25±1)℃,相對濕度穩(wěn)定在(20±3)%,體系壓力始終保持為1.0′105Pa.

圖2 空白組中氣溶膠Nt及D50分布

SO2=0mg/m3, NO2=0mg/m3, NH3=0

1.2 實驗方法與工況設(shè)置

實驗前,打開黑光燈,使用凈化空氣以30L/min的流量沖洗煙霧箱4h,連接SMPS在線測量煙霧箱中的顆粒物,待顆粒物數(shù)濃度穩(wěn)定至500cm-3以下即清潔完成,記錄此時煙霧箱內(nèi)顆粒物粒徑分布為=0min時的顆粒物粒徑分布狀態(tài).封閉煙霧箱氣體出口,將凈化空氣流量降至1L/min以作為載氣,采用氣密性進(jìn)樣針將一定體積的SO2、NO2、NH3同時注入載氣中,封閉煙霧箱,開啟SMPS在線監(jiān)測煙霧箱內(nèi)的氣溶膠濃度及粒徑變化情況.

TSI 3938L72型掃描電遷移率粒徑譜儀可以對10~1000nm范圍內(nèi)的顆粒物粒徑分布進(jìn)行測量;它主要由3082型靜電分級器(DMA)和3772型凝聚核粒子計數(shù)器(CPC)組成. DMA用來選擇、測量顆粒物的粒徑, CPC 用來測定顆粒物的數(shù)濃度. SMPS進(jìn)樣流量為0.25L/min,設(shè)置粒徑測量下限19.1nm,粒徑測量上限649.4nm,掃描時間及吹掃時間均為60s.測量模式為連續(xù)測量14組,單工況用時30min.

2011年發(fā)布的《火電廠大氣污染物排放標(biāo)準(zhǔn)》[15]規(guī)定:自2014年7月1日起,非重點地區(qū)燃煤鍋爐SO2及NO(以NO2計)限值為100mg/m3.因此實驗設(shè)置SO2、NO2濃度為0~200mg/m3.電廠運行經(jīng)驗表明[16]采用SCR脫硝技術(shù)的機(jī)組氨逃逸濃度一般不超過3′10-6,因此實驗設(shè)置NH3濃度為0~12′10-6.實驗共計13個工況,如表1所示.

表1 實驗工況設(shè)計

SMPS有64個掃描通道,其粒徑掃描范圍為19.1~649.4nm,將粒徑對數(shù)化后平均分為98段,每一小段的顆粒物個數(shù)濃度為N.對N按式(1)進(jìn)行求和即得該工況下的單通道顆粒物總個數(shù)濃度N.記每一工況下顆粒物N隨時間變化的峰值為max.

顆粒物的累計粒度數(shù)量分布百分?jǐn)?shù)達(dá)到50%時所對應(yīng)的粒徑為中值粒徑50,其對應(yīng)的顆粒物個數(shù)濃度記為50.由于每個工況均進(jìn)行了重復(fù)實驗,因此本文中N、50及50均為算術(shù)平均后的結(jié)果.

2 結(jié)果與討論

2.1 顆粒物生成與團(tuán)聚煙霧箱實驗結(jié)果

2.1.1 SO2濃度對顆粒物生成特性的影響 不同SO2濃度下,煙霧箱內(nèi)顆粒物N及50隨時間變化的結(jié)果如圖3所示.SO2、NO2及NH3注入煙霧箱后迅速發(fā)生反應(yīng),產(chǎn)生了大量氣溶膠顆粒,持續(xù)時間較短.在實驗過程中,顆粒物總個數(shù)濃度首先由凈空氣中的102cm-3急劇上升至106cm-3,而后成核作用減弱,團(tuán)聚作用使得顆粒物總個數(shù)濃度在5~10min內(nèi)迅速衰減,最終在30min時下降至5′105cm-3,而50增大了近一倍.工況1為反應(yīng)物濃度最高的工況,其max高達(dá)2.5′106cm-3.隨著不同工況SO2注入量的減少,化學(xué)反應(yīng)速率下降,煙霧箱內(nèi)產(chǎn)生的氣溶膠顆?？倐€數(shù)濃度有序降低,進(jìn)而導(dǎo)致顆粒碰撞頻率降低,團(tuán)聚作用減弱,因此氣溶膠顆粒的50也降低.工況5為無SO2工況,實驗中仍存在大量新顆粒物生成,工況5的氣溶膠max與工況1相比減少了41.0%.可見無SO2的條件下,NO2、NH3及H2O反應(yīng)體系亦可產(chǎn)生大量氣溶膠顆粒物,但是引入SO2可大大提高氣溶膠成核強(qiáng)度.

NO2=200mg/m3NH3=12×10-6

2.1.2 NO2濃度對顆粒物生成特性的影響 如圖4所示為不同NO2濃度下,氣溶膠N及50隨時間的變化曲線.可見SO2與NO2在反應(yīng)體系中有一致的促進(jìn)成核作用,但是在SO2-NO2-NH3-H2O四元反應(yīng)體系中,NO2不僅參與硝酸鹽的生成,而且對SO2向SO42-的氧化過程具有顯著的促進(jìn)作用[17].同時由于該反應(yīng)生成H+,反應(yīng)體系中的NH3可通過酸堿中和吸收H+,加快正反應(yīng)速率,促進(jìn)SO2氧化反應(yīng)的化學(xué)平衡正向移動,產(chǎn)生更多硫酸鹽與銨鹽.

圖4 NO2濃度對氣溶膠Nt及D50的影響

SO2=200mg/m3NH3=12×10-6

工況9為無NO2工況,在反應(yīng)開始的6min內(nèi),氣溶膠顆?？倐€數(shù)濃度曲線緩慢上升,此時成核作用強(qiáng)于團(tuán)聚作用,煙霧箱中有大量新顆粒物生成.新顆粒物粒徑小,因此在2~4min前表現(xiàn)為整體顆粒物50下降;此后,成核作用不斷減弱,6~10min時氣溶膠數(shù)濃度曲線趨于平穩(wěn);10min后氣溶膠數(shù)濃度曲線開始緩慢下降,團(tuán)聚作用占主導(dǎo)地位,顆粒物50緩慢上升.煙霧箱中無NO2時,顆粒物N變化緩慢,此時max與工況1相比減少了83.6%,可見SO2、NH3與H2O反應(yīng)體系的成核速率較慢.楊衛(wèi)衛(wèi)[18]對礦質(zhì)氧化物表面硫酸鹽、硝酸鹽和銨鹽的生成機(jī)制研究發(fā)現(xiàn)NH3對SO2異相轉(zhuǎn)化有催化作用, Benner等[19]也發(fā)現(xiàn)NH3可極大的加速SO2向SO42-的轉(zhuǎn)化,產(chǎn)生以(NH4)2SO4為主的硫酸鹽.所以工況10雖然沒有NO2參與反應(yīng),但仍可觀察到氣溶膠顆粒生成可能存在NH3對SO2的催化作用,產(chǎn)生的顆粒物可能為硫酸銨、硫酸氫銨及部分亞硫酸鹽.

2.1.3 NH3濃度對顆粒物生成特性的影響 圖5為不同NH3濃度下的氣溶膠顆粒N、50隨時間變化的曲線.不同NH3濃度時,氣溶膠總個數(shù)濃度及粒徑分布總體趨勢變化不大,說明在該反應(yīng)體系下,NH3濃度變化對氣溶膠生成的影響較小.但是當(dāng)NH3濃度降為0時,產(chǎn)生新顆粒物極少,在6分鐘時N達(dá)到峰值,max僅有3.8′104cm-3,與工況1相比減少了98.5%.由于此時煙霧箱內(nèi)氣溶膠顆粒較少,碰撞概率低,團(tuán)聚作用弱,因此新顆粒物粒徑增長較慢.再加上新顆粒物粒徑小,所以該工況下顆粒物50僅為41.6nm,遠(yuǎn)小于其他工況的顆粒物50.在SO2-NO2- NH3-H2O四元反應(yīng)體系中,NH3除了對SO2的氧化具有催化作用,NH4+作為主要陽離子對化學(xué)反應(yīng)的進(jìn)行及氣溶膠的產(chǎn)生具有無可取代的作用.在NH3濃度較高時,反應(yīng)體系生成物可能以NH4NO3及(NH4)2SO4為主;在NH3濃度較低時,反應(yīng)產(chǎn)物可能為NH4NO3、NH4HSO4;在無NH3時,反應(yīng)體系最終可能以HNO3及H2SO4氣溶膠為主.

SO2=200mg/m3NO2=200mg/m3

2.2 氣溶膠生成與團(tuán)聚計算結(jié)果

不同氣氛下的煙霧箱實驗結(jié)果顯示了SO2- NO2-NH3-H2O四元反應(yīng)體系具有明顯的成核、團(tuán)聚現(xiàn)象,其氣溶膠成核及團(tuán)聚作用強(qiáng)度與氣態(tài)前體物濃度呈正相關(guān).為了歸納總結(jié)SO2-NO2-NH3-H2O四元反應(yīng)體系的氣-粒轉(zhuǎn)化特性,本文進(jìn)一步對氣溶膠生成及團(tuán)聚過程進(jìn)行擬合計算與模型驗證.

2.2.1 氣溶膠生成特性擬合計算 新顆粒物的生成與成長過程包括顆粒物的成核過程和之后的粒徑成長過程[20].通常意義上,粒徑小于1μm的顆粒物可以分為凝結(jié)核模態(tài)(小于30nm)、Aitken模態(tài)(30~100nm)及積聚模態(tài)(100~1000nm)[21].不同模態(tài)的顆粒物來源不同,凝結(jié)核模態(tài)粒子是由氣態(tài)前體物成核形成的;Aitken模態(tài)顆粒則是由半揮發(fā)性的氣體在固體顆粒物表面凝聚而成;積聚模態(tài)粒子可通過更小模態(tài)粒子碰并形成[22].

實驗中,由于儀器難以測量粒徑小于15nm的凝結(jié)核模態(tài)顆粒物,因此本文中應(yīng)用不同氣氛下經(jīng)過相同反應(yīng)深度后,反應(yīng)體系產(chǎn)生的顆粒物總個數(shù)濃度凈增長量來定義新顆粒物的成核強(qiáng)度.以工況1為例,使用SMPS檢測煙霧箱內(nèi)顆粒物N最大時的粒徑分布情況,對顆粒物粒徑進(jìn)行對數(shù)運算后對顆粒物數(shù)濃度作正態(tài)擬合曲線如圖6所示.擬合方程為式(3).

式中:2.10′106表示擬合曲線下的面積,其物理意義為氣溶膠顆粒的總個數(shù)濃度;=116.3為擬合曲線的均值,其物理意義為氣溶膠顆粒的中值粒徑50;=0.3810表示擬合曲線的標(biāo)準(zhǔn)差.

方程的調(diào)整確定系數(shù)(Adj.2)為0.9912,擬合優(yōu)度較好,新顆粒物粒徑滿足正態(tài)分布.因此,對于未知氣溶膠初始粒徑分布的反應(yīng)體系,可以通過擬合方程模擬實際粒徑分布曲線.對于給定濃度的SO2、NO2、NH3氣氛,只需確定氣態(tài)前體物的濃度與氣溶膠顆粒的N、50及50的關(guān)系即可通過下式求解出新顆粒物的粒徑分布模擬曲線.

對化學(xué)反應(yīng)所引起的成核問題,此時一般認(rèn)為核體的產(chǎn)生速率與前體物的數(shù)目濃度成線性關(guān)系[23].Wang等[17]對霧霾過程中硫酸鹽持續(xù)性生成研究發(fā)現(xiàn)高相對濕度、NH3及NO2氛圍是SO2液相氧化產(chǎn)生大量硫酸鹽的關(guān)鍵影響因素.當(dāng)反應(yīng)體系中缺少任一種氣體時,由于復(fù)雜的復(fù)合效應(yīng)會使得體系產(chǎn)生的氣溶膠顆粒物大幅度減少.考慮到實際情況下,電廠排放的氣體污染物也不存在工況5、9、13所示的無SO2、無NO2、無NH3的情況.因此,在對實驗數(shù)據(jù)進(jìn)行擬合及模擬計算時,本文僅選擇工況1,2,3,4,6,7,8,10,11,12分組進(jìn)行擬合.對SO2、NO2和NH3濃度與顆粒物N、50及50的關(guān)系分別擬合作圖,如圖7所示.擬合方程的參數(shù)如表2所示.擬合方程的Adj.2均接近1,表明實驗數(shù)據(jù)的線性度較好.擬合曲線的斜率反映了該氣態(tài)前體物濃度變化對新顆粒物生成過程強(qiáng)度的影響程度.在實驗所研究的電廠污染物排放濃度范圍內(nèi),NO2在SO2- NO2-NH3-H2O四元反應(yīng)體系中對氣溶膠顆粒成核強(qiáng)度影響大于SO2.實際應(yīng)用中,燃煤電廠控制NO排放濃度對改善煙囪出口氣溶膠數(shù)濃度更有效.

2.2.2 氣溶膠團(tuán)聚模型與計算 對于恒溫靜置煙霧箱,其內(nèi)部產(chǎn)生的氣溶膠顆粒因布朗運動而產(chǎn)生團(tuán)聚.團(tuán)聚過程表現(xiàn)為氣溶膠顆粒粒徑增大,而總個數(shù)濃度減小.對于粒徑為d,體積為v的氣溶膠顆粒,其數(shù)濃度n隨時間變化可表示為式(5)所示的氣溶膠團(tuán)聚動力學(xué)方程[24].

式中:n、n分別代表粒徑為d、d,體積為v、v的顆粒、個數(shù)濃度,滿足v+v=v,即表示顆粒由顆粒和顆粒團(tuán)聚生成,滿足體積守恒條件.(,)表示顆粒和顆粒的團(tuán)聚核函數(shù),定義為單位時間內(nèi)單位數(shù)目濃度顆粒與顆粒發(fā)生團(tuán)聚的次數(shù).方程右邊第一項表示其他小顆粒團(tuán)聚生成顆粒引起的顆粒數(shù)目濃度增加,第二項表示顆粒與其他顆粒團(tuán)聚為大顆粒引起的顆粒數(shù)目濃度減少.

表2 擬合方程參數(shù)

將上述微分方程離散化,得Smoluchowski方程[25]:

實驗中氣溶膠顆粒粒徑為幾十~幾百nm,顆粒位于過渡區(qū),布朗團(tuán)聚的團(tuán)聚核函數(shù)可按下式(7)~式(12)計算[26]:

本文采用Matlab編程,以每一工況2min時的氣溶膠顆粒粒徑分布為團(tuán)聚起始狀態(tài),基于上述過渡區(qū)顆粒碰撞核函數(shù)及Smoluchowski方程模擬此后的氣溶膠粒徑分布,獲得不同工況下的氣溶膠粒徑分布隨時間變化曲線.以工況1為例,由于每一個工況均進(jìn)行了3次重復(fù)試驗,在進(jìn)行模擬計算時,選擇統(tǒng)計數(shù)據(jù)較好的一次試驗結(jié)果作為模擬計算的起始粒徑分布曲線.其氣溶膠粒徑分布隨時間變化曲線如圖8所示,氣溶膠顆粒N及50變化曲線如圖9所示.可見模擬曲線較好地預(yù)測了氣溶膠粒徑分布隨時間變化的規(guī)律,采用該模型能較好地反應(yīng)SO2-NO2-NH3-H2O四元反應(yīng)體系成核后的團(tuán)聚過程.

圖8 工況1氣溶膠粒徑分布的實驗值與模擬值

通過上述方程擬合及模擬計算,對于燃煤電廠某時刻的SO2、NO2、NH3排放濃度,通過如表2所示的擬合方程可直接計算出該工況下可能產(chǎn)生的氣溶膠顆粒的N、50及50,將獲得的參數(shù)代入式(3)所示擬合方程可獲得該工況下的顆粒物粒徑分布模擬曲線.最后,以顆粒物粒徑分布模擬曲線為團(tuán)聚初始狀態(tài),通過布朗團(tuán)聚模型即可預(yù)測此后的氣溶膠顆粒N及50變化情況.

3 結(jié)論

3.1 SO2-NO2-NH3-H2O四元反應(yīng)體系具有顯著的成核能力,且其成核強(qiáng)度大,持續(xù)時間短.當(dāng)SO2、NO2濃度為200mg/m3,NH3濃度為12′10-6時,反應(yīng)體系的氣溶膠總個數(shù)濃度在2min時即達(dá)到峰值2.5′106cm-3.

3.2 缺少任一種氣體均會使氣溶膠成核強(qiáng)度下降.SO2、NO2及NH3濃度分別為0的工況下氣溶膠總個數(shù)濃度峰值分別下降了41.0%、83.6%及98.5%.在電廠污染氣體排放濃度區(qū)間內(nèi),NO2對氣溶膠生成影響大于SO2.燃煤電廠控制NO排放濃度對改善煙囪出口氣溶膠數(shù)濃度更有效.

3.3 氣-粒轉(zhuǎn)化產(chǎn)生的銨鹽氣溶膠總個數(shù)濃度、中值粒徑與氣態(tài)前體物濃度在電廠污染物排放濃度范圍內(nèi)線性相關(guān).反應(yīng)體系產(chǎn)生的顆粒物數(shù)濃度-粒徑滿足正態(tài)分布.采用布朗團(tuán)聚模型可較好地模擬氣溶膠的團(tuán)聚過程,預(yù)測其后的氣溶膠顆粒數(shù)濃度及粒徑分布變化.

[1] Yan J, Chen L, Lin Q, et al. Chemical characteristics of submicron aerosol particles during a long-lasting haze episode in Xiamen, China [J]. Atmospheric Environment, 2015,113:118-126.

[2] Wang Q, Zhuang G, Huang K, et al. Probing the severe haze pollution in three typical regions of China: Characteristics, sources and regional impacts [J]. Atmospheric Environment, 2015,120:76-88.

[3] 謝丹丹,祁建華,張瑞峰.青島不同強(qiáng)度霾天氣溶膠中二次無機(jī)離子的生成及粒徑分布[J]. 環(huán)境科學(xué), 2017,38(7):2667-2678. Xie D D, Qi J H, Zhang R F. Formation and size distribution of the secondary aerosol inorganic ions in different intensity of haze in qingdao [J]. Environmental Science, 2017,38(7):2667-2678.

[4] 陳 靜,楊 鵬,韓軍彩,等.基于高分辨率MARGA數(shù)據(jù)分析石家莊PM2.5成分譜特征 [J]. 中國環(huán)境科學(xué), 2015,35(9):2594-2604. Chen J, Yang P, Han J C, et al. Analysis of PM2.5spectrum characteristics in Shijiazhuang based on high resolution MARGA data [J]. China Environmental Science, 2015,35(9):2594-2604.

[5] 郭振東,朱 彬,王紅磊,等.長江三角洲霾天氣PM2.5中水溶性離子特征及來源解析 [J]. 中國環(huán)境科學(xué), 2019,39(3):928-938. Guo Z D, Zhu B, Wang H L, et al. Characteristics and source analysis of water-soluble ions in PM2.5in the haze weather over in Yangtze River Delta [J]. China Environmental Science, 2019,39(3):928-938.

[6] 王曼婷,朱 彬,王紅磊,等.長三角冬季一次霾過程氣溶膠及其水溶性離子的區(qū)域分布特征[J]. 環(huán)境科學(xué), 2015,36(7):2337-2345. Wang M T, Zhu B, Wang H L, et al. Composition and regional characteristics of atmosphere aerosol and its water soluble ions over the Yangtze River Delta Region in a winter haze period [J]. Environmental Science, 2015,36(7):2337-2345.

[7] 譚成好,趙天良,崔春光,等.近50年華中地區(qū)霾污染的特征 [J]. 中國環(huán)境科學(xué), 2015,35(8):2272-2280. Tan C H, Zhao T L, Cui C G, et al. Characterization of haze pollution over Central China during the past 50years [J]. China Environmental Science, 2015,35(8):2272-2280.

[8] 郭送軍,譚吉華,段菁春,等.廣州市灰霾期大氣PM10中水溶性離子特征[J]. 環(huán)境科學(xué)與技術(shù), 2012,35(11):83-86. Guo S J, Tan J H, Duan J C, et al. Characteristics of water soluble ions in atmospheric PM10during haze periods in Guangzhou [J]. Environmental Science & Technology, 2012,35(11):83-86.

[9] 王開燕,王蓓蕾,王春林,等.廣州2017年1月一次嚴(yán)重灰霾過程的影響因素分析 [J]. 環(huán)境科技, 2018,31(2):79-84. Wang K Y, Wang B L, Wang C L, et al. Analysis of influencing factors during a serious haze process in Guangzhou in January 2017 [J]. Environmental Science and Technology, 2018,31(2):79-84.

[10] 麥健華,鄧 濤,黃燁琪,等.中山市一次灰霾天氣過程污染物來源數(shù)值模擬分析 [J]. 中國環(huán)境科學(xué), 2017,37(9):3258-3267. Mai J H, Deng T, Huang Y Q, et al. Source analysis of a haze event in Zhongshan by numerical simulation [J]. China Environmental Science, 2017,37(9):3258-3267.

[11] Xue J, Yuan Z, Griffith S M, et al. Sulfate Formation Enhanced by a Cocktail of High NO, SO2, Particulate Matter, and Droplet pH during Haze-Fog Events in Megacities in China: An Observation-Based Modeling Investigation [J]. Environmental Science & Technology, 2016,50(14):7325-7334.

[12] 武衛(wèi)玲,陳良富,陶金花,等.利用遙感數(shù)據(jù)評價燃煤電廠空氣質(zhì)量 [J]. 遙感學(xué)報, 2013,17(5):1235-1245. Wu W L, Chen L F, Tao J H, et al. Assessment of air in coal-fired power plants using satellite observations [J]. Journal of Remote Sensing, 2013,17(5):1235-1245.

[13] 沙 桐,馬曉燕,王健穎,等.長江三角洲冬季電廠排放對大氣污染的影響 [J]. 中國環(huán)境科學(xué), 2018,38(9):90-99. Sha T, Ma X Y, Wang J Y, et al. The impact of power plant emission on air pollution during winter over Yangtze River Delta [J]. China Environmental Science, 2018,38(9):90-99.

[14] 韓力慧,向 欣,張海亮,等.北京市開發(fā)區(qū)PM1污染特征及影響霾形成的因素 [J]. 中國環(huán)境科學(xué), 2018,38(8):48-58. Han L H, Xiang X, Zhang H L, et al. Pollution characteristics of PM1and factors affecting the formation of haze pollution at a developed zone in Beijing [J]. China Environmental Science, 2018, 38(8):48-58.

[15] GB 13223-2011 火電廠大氣污染物排放標(biāo)準(zhǔn) [S]. GB 13223-2011 Emission standard of air pollutants for thermal power plants [S].

[16] 馬雙忱,金 鑫,孫云雪,等. SCR煙氣脫硝過程硫酸氫銨的生成機(jī)理與控制 [J]. 熱力發(fā)電, 2010,39(8):12-17. Ma S C, Jin X, Sun Y X, et al. The formation mechanism of ammonium bisulfate in SCR flue gas denitrification process and control thereof [J]. Thermal power generation, 2010,39(8):12-17.

[17] Wang G, Zhang R, Gomez M E, et al. Persistent sulfate formation from London Fog to Chinese haze [J]. Proceedings of the National Academy of Sciences of the United States of America, 2016,48(113):13630- 13635.

[18] Yang W W, He H, Ma Q X, et al. Synergistic formation of sulfate and ammonium resulting from reaction between SO2and NH3on typical mineral dust [J]. Phys. Chem. Chem. Phys., 2016,18:956-964.

[19] Benner W H, Ogorevc B, Novakov T. Oxidation of SO2in thin water films containing NH3[J]. Atmospheric Environment. Part A. general Topics, 1992,26(9):1713-1723.

[20] Kulmala M, Maso M D, M?kel? J M, et al. On the formation, growth and composition of nucleation mode particles [J]. Tellus Series B-chemical & Physical Meteorology, 2001,53(4):479–490.

[21] Gieré R, Querol X. Solid particulate matter in the atmosphere [J]. Elements, 2010,6(4):215-222.

[22] Vu T V, Delgado-Saborit J M, Harrison R M. Review: Particle number size distributions from seven major sources and implications for source apportionment studies [J]. Atmospheric Environment, 2015, 122:114-132.

[23] Landgrebe J D, Pratsinis S E. Gas-phase manufacture of particulates: interplay of chemical reaction and aerosol coagulation in the free- molecular regime [J]. Industrial & Engineering Chemistry Research, 1989,28(10):1474-1481.

[24] Friedlander S K, Marlow W H. Smoke, Dust and Haze: Fundamentals of Aerosol Behavior [J]. Physics Today, 1977,30(9):58-59.

[25] Sarabia R F D, Elvira-Segura L, I González-Gómez, et al. Investigation of the influence of humidity on the ultrasonic agglomeration of submicron particles in diesel exhausts [J]. Ultrasonics, 2004,41(4):277-281.

[26] Park S H, Lee K W, Otto E, et al. The log-normal size distribution theory of Brownian agglomeration for the entire particle size range [J]. Journal of Aerosol Science, 1999,30(28):S23–S24.

Formation characteristics of aerosols in the SO2-NO2-NH3-H2O quaternary system.

XIAO Jia-fan, LIU He-xin, TAN Hou-zhang*, DU Yong-le, WANG Yi-bin

(Key Laboratory of Thermo-Fluid Science and Engineering, Ministry of Education, Xi’an Jiaotong University, Xi’an 710049, China)., 2019,39(7)：2721~2729

The variations in concentration and diameter distribution of newly formed particulate matter, during the gas-to-solid conversion of the SO2-NO2-NH3-H2O quaternary system, were experimentally investigated with the smog chamber in this article. This quaternary system had a significant capability of nucleation, and a large nucleation intensity in such a short duration.When both of SO2and NO2concentrations were 200mg/m3and NH3concentration was 12×10-6, the total concentration of aerosols in this system reached a peak of 2.5×106cm-3over two minute. Moreover, the nucleation strength of the aerosol would be reduced lacking in any one among three gases.Like this typical condition that SO2, NO2or NH3concentration was zero, the maximum concentration of aerosols decreased by 41.0%, 83.6% and 98.5%, respectively. Within the concentration ranges of SO2, NO2or NH3emitted from coal-fired power plants, the NO2had a greater influence than SO2on the formation and agglomeration of newly formed particles. Based on experimental results, the fitting curves for nucleation characteristics of particulate matters were made. A linear correlation was observed between the total number concentration of newly formed particles and median particle diameter, and gaseous pollutant concentrations. The Brown agglomeration model was adopted to simulate the changes in total number concentration and particle size distribution of aerosols after nucleation in the agglomeration process. Finally, a method was proposed for predicting the change in the particle size distribution and total number concentration of aerosol particles, based on the typical emission concentrations of SO2, NO2and NH3from coal-fired power plants.

coal-fired power plant；haze；aerosol；smog chamber；nucleation；agglomeration

X515

A

1000-6923(2019)07-2721-09

蕭嘉繁(1994-),男,福建晉江人,西安交通大學(xué)碩士研究生,主要從事燃煤電廠大氣污染物控制研究.發(fā)表論文1篇.

2018-12-11

國家重點研發(fā)計劃(2016YFB0601504)

*責(zé)任作者, 教授, tanhz@mail.xjtu.edu.cn