石灰石脫硫的反應活性分析論文

導語：石灰石脫硫的反應活性分析論文一文來源于網友上傳，不代表本站觀點，若需要原創文章可咨詢客服老師，歡迎參考。

向流化床鍋爐的燃燒區加入石灰石，首先發生的是石灰石的高溫分解，分解產物為CaO。CaO顆粒在O2過量條件下，與SO2發生硫化反應，生成CaSO4，即：CaO+SO2+1/2O2=CaSO4+486kJ/mol(1)

石灰石的反應活性對反應式(1)的反應程度影響很大。因此，國內外研究人員對此進行了大量工作〔1～3〕，實驗主要采用熱重分析(TGA)法，測定的對象一般為CaO的硫酸鹽化程度，并以此為基礎研究石灰石的脫硫反應活性。由于該方法的測定對象為固態，故簡稱之為“固測法”。但是，用固測法研究石灰石的活性，有時存在較大的偏差。TGA法的吸硫曲線的增重趨勢總是被認為是按反應式(1)中CaO吸收SO2和O2生成CaSO4所造成。而實際上，石灰石中除了主要成分CaCO3外，還含有許多其它雜質成分。一些雜質成分經高溫分解后產生的一些堿性氧化物同樣也能與SO2和O2反應生成硫酸鹽，另外一些雜質成分還會生成一些復雜的復合物，從而間接地影響反應(1)的進行，這些情況會給熱重分析帶來不可避免的誤差，而影響石灰石活性數據的準確性。

本文的研究方法是將測定對象由TGA法的固測改為對SO2的氣測。在流化床脫硫模擬試驗臺上，通過監測SO2濃度的變化研究石灰石的反應活性。在特定的工況條件下，SO2在通過吸硫劑石灰石時濃度變化可以認為是由于石灰石固有特性所引起的，它體現了石灰石總體吸硫效果，這樣就避免了由于僅僅考慮CaO的轉化而忽略了其它雜質成分影響帶來的活性數值的偏差。本文將以流化床典型運行溫度850℃下的反應速度常數作為石灰石脫硫活性的指標。

1材料與方法

試驗采用流化床反應器模擬法。全部試驗在石灰石脫硫反應活性試驗臺上進行，圖1為試驗流程簡圖。N2、CO2和SO2經流量計A進入氣體混合器B，在a點用MIS-2000煙道氣體分析儀E對混合氣體的SO2初始濃度進行測定。當SO2流經a點進入反應器C時，脫硫反應開始進行。每間隔10min在b點用氣體分析儀E對SO2的濃度進行一次測定，直到SO2的濃度不再改變時為止，此時，試驗結束。為了模擬流化床的實際工況，石灰石的熱解與脫硫在反應器中同時進行(這與TGA法不同)，反應溫度設置為流化床典型運行溫度850℃。

A.浮子式氣體流量計B.氣體混合器

C.流化床反應器D.隔膜泵E.MSI-2000氣體分析儀

Fig.1Diagramofexperimentunit

A.flaotgasmeterB.gasmixer

C.fluidizedbedreactorD.meteringpumpMSI-2000gasanalyzer

2動力學模型的建立

反應式(1)的速度公式為：

r=-dCA/dt=KCxACyB(2)

式中：r為脫硫反應速度；t為反應時間，CA為t時刻的SO2濃度；CB為t時刻的O2濃度；x和y分別為SO2和O2的表觀反應級數。

Lee〔4〕等人研究石灰石在流化床中與SO2的反應活性時指出：因為反應是在氧過量的條件下進行的，所以固硫化反應對于SO2是一級反應，而對O2是零級反應。于是，根據這一結論可將(2)式簡化為：

r=-dCA/dt=KCA(3)

采用積分法，將(3)式移項積分

∫cc°dCA/CA=∫°tKdt(4)

當t=0時，SO2的濃度為C0；當t=t時，SO2的濃度為C。對(4)式積分后即得到：

C=C0exp(-Kt)(5)

3石灰石脫硫反應活性的求取

將(5)式兩邊取對數，得

lnC=lnC0-Kt(6)

(6)式可轉換為：Y=b+aX(7)

由式(7)可見，變量X與Y之間成線性關系。因此，可以線性回歸確定參數a、b及線性相關系數R。使用ppp-60程序上機運算，結果列于表1。

由表1可見，線性相關系數在0.95以上，查t分布的ta值表，t0.0005=-5.405。表1中所有的t檢驗值均小于-5.405，說明有1-α=1-0.0005=99.95%的把握拒絕μ=0的零假設〔5〕。

4石灰石化學成分與活性關系

以反應活性K作為因變量，化學成分為自變量，建立多元線性回歸方程。為了保證方程的穩定性，采用逐步回歸分析法〔6〕將一些次要成分剔除，最后得到回歸方程為：

=0.01382+9.471×10-5〔CaO〕-1.432×10-4〔MgCO3〕-3.084×10-5〔SiO2〕-1.054×10-3〔Al2O3〕+1.080×10-3〔Fe2O3〕(8)

復相關系數R=0.9722，殘差平方和Q=1.275×10-5，剩余標準差S=9.544×10-4，F檢驗值F=80.48?？梢?，方程(8)的回歸擬合優度較高。

表1運算結果

Table1Calculatingresults

序回歸式=b+aX反應速度常數相關系數Rt檢驗植

號850℃

1=8.000-1.794×10-2X1.794×10-20.9712-9.990

2=7.886-1.374×10-2X1.374×10-20.9506-7.505

3=7.968-8.843×10-3X8.843×10-30.9960-27.42

4=7.807-1.794×10-2X1.794×10-20.9764-11.08

5=7.912-1.409×10-2X1.409×10-20.9936-21.63

6=7.884-1.708×10-2X1.708×10-20.9885-16.03

7=7.849-7.160×10-3X7.160×10-30.9644-8.935

8=7.792-1.795×10-2X1.795×10-20.9794-11.89

9=7.879-1.519×10-2X1.519×10-20.9891-16.48

10=7.778-1.665×10-2X1.665×10-20.9837-13.38

11=7.752-1.808×10-2X1.808×10-20.9864-14.70

12=7.857-1.022×10-2X1.022×10-20.9708-9.905

13=7.857-1.431×10-2X1.431×10-20.9831-13.15

14=7.794-2.001×10-2X2.001×10-20.9816-15.43

15=7.777-2.048×10-2X2.048×10-20.9833-13.23

16=7.783-1.677×10-2X1.677×10-20.9760-10.96

17=7.868-1.510×10-2X1.510×10-20.9788-11.71

18=7.788-1.788×10-2X1.788×10-20.9714-10.02

19=7.900-1.811×10-2X1.811×10-20.9761-11.01

20=7.844-1.749×10-2X1.749×10-20.9790-11.76

5結論

5.1采用氣測法獲得的石灰石脫硫反應速率常數K，可以全面地反映石灰石中脫硫能力，避免了固測法由于石灰石雜質造成的誤差。

5.2由式(8)可知，石灰石的反應活性與雜質成分有一定的數量關系。值得注意的是〔MgCO3〕在方程中對K具有負效應，這似乎有些不可思議，這是因為〔MgCO3〕的分解產物MgO的反應活性與壓力有關。MgO由于在常壓下與SO2的反應非常之慢，可視為惰性，所以MgCO3一般在增壓流化床內才具有良好的脫硫作用。

5.3利用式(8)可對石灰石的反應活性進行預測，例如：當〔CaO〕=50%,〔MgCO3〕=6.2%,〔SiO2〕=3.42%,〔Fe2O3〕=0.08%,〔Al2O3〕=0.2%時，可求得K=1.763×10-2。

參考文獻

1DennisJS，HayhurstAN.ChemicalEngineeringScience,1986,41:25～36

2JamesRH,etal.TestingofLimestoneSamplefromtheTVARegionasSulfurDioxideSorbentinAFBC,InstituteforMiningangMineralsandResearchKentuckyCenterforEnergy,ResearchLaboratoryandTennessee,ValleyAuthority,19843國家環保局.大氣防治技術研究.北京：科學出版社，1993.370～373

4LeeDC,etal.AICHEJ,1981,27:472～480

5中國科學院數學研究所.回歸分析方法.北京：科學出版社，1974.29～32

6邢英等.經濟統計分析及預測.北京：清華大學出版社，1982.109～119