氧脫木素動力學模型的建立

氧脫木素動力學模型的建立

氧氣脫木素過程是以naoh溶液為介質(zhì)和氧化氧化的過程。在制備砂漿泡沫行業(yè)中，提出了將殘余木材氧化脫離出來的方法。相比傳統(tǒng)的蒸煮深度脫木素及含氯漂白而言,它能顯著提高紙漿的得率,減少后續(xù)漂白過程化學品的消耗量,有效降低漂白廢液中的COD、BOD、有機氯化物含量及色值,從而提高漂白過程中廢水的回用率。為了進一步發(fā)揮現(xiàn)行氧脫木素工藝節(jié)能減排的功效,并為工藝的改進提供指導,有必要對該過程木素脫除的規(guī)律進行深入的研究。目前,用于描述氧脫木素過程中木素脫除動態(tài)過程規(guī)律的數(shù)學模型主要有兩種,即:將紙漿中的殘余木素根據(jù)其反應活性的不同進而建立的雙一級脫木素模型和假設脫木素反應為非一級反應、用單一動力學方程描述氧脫木素過程的數(shù)學模型。在關于建立接近工業(yè)實踐氧脫木素體系的脫木素模型的報道中,由于均假設工藝的主要參數(shù)之一——NaOH濃度在反應過程中保持不變,因此這種經(jīng)驗模型并沒有反映出氧脫木素體系的真實情況。其結果不僅導致擬合出的各種模型參數(shù)具有一定的偏差,而且會增加模型參數(shù)的變異性,進而影響了模型的準確性和可靠性。很顯然,在氧脫木素工藝中紙漿上殘余木素的脫除與堿的加入量有關,在實際氧脫木素過程中NaOH濃度也是隨著反應的進行發(fā)生顯著的變化,Vincent等人對木素脫除反應中化學計量的堿的消耗計算結果表明,當木素脫除50%時,堿濃度約降低3%~8%。然而本實驗對有機酸含量較多的桉木漿氧脫木素過程中殘余有效堿的實際測定結果表明,當木素脫除30%時,由于大部分有效堿參與中和體系中的有機酸,堿濃度約降低43%~95%,將直接影響到后續(xù)脫木素反應的速率。盡管有研究采用超低漿濃以達到NaOH濃度恒定的目的,進而建立了脫木素動力學模型;但是,這種模型適用的體系與工業(yè)生產(chǎn)的實際情況相差甚遠。因此,只有將NaOH濃度的變化引入經(jīng)驗模型中進而建立氧脫木素過程的數(shù)學模型,才能更客觀地描述氧脫木素過程中木素脫除的動態(tài)規(guī)律。桉木屬于闊葉木種,是近年來在我國大力推廣種植的一種制漿原料樹種。本實驗將以桉木的硫酸鹽漿為研究對象,通過改變氧脫木素過程中的3個主要工藝參數(shù)(溫度、氧壓和用堿量),考察該過程中紙漿上殘余木素(用卡伯值表征)的脫除規(guī)律及殘堿的變化規(guī)律;將過程中殘堿的變化作為變量之一,結合氧脫木素的動力學方程導出了木素脫除的數(shù)學模型,從而建立能更客觀描述氧脫木素規(guī)律的數(shù)學模型。1實驗1.1反應轅/蒸煮器桉木(由廣東省雷州林業(yè)局提供)硫酸鹽漿為實驗室自制,經(jīng)洗滌、篩選、風干后備用,其初始卡伯值為32.78。反應釜(型號為Parr845);蒸煮器(型號為Greenwood2110-2);NaOH溶液(質(zhì)量濃度51.74g/L);硫酸鎂溶液(質(zhì)量分數(shù)為1%)。1.2聚乙烯目標液的制備氧脫木素工藝條件如表1所示。每次氧脫木素實驗前,先將相當20.00g絕干漿的未漂風干漿用一定量的水(保證漿濃10%的總加水量-藥品溶液加入量-清洗封口袋用水量)在聚乙烯封口袋中先浸泡約10h,然后將用量0.2%的MgSO4(對絕干漿)溶液及NaOH溶液加入封口袋中并分別搓揉,再將混合均勻的漿料全部倒入反應器內(nèi),用一定量(約40mL)水清洗封口袋并一起移入反應器內(nèi)。密封反應器,同時啟動轉子(300r/min),開始升溫,待O2升到目標溫度時,通入O2開始計時,在反應進行至指定時間(0min、15min、30min、60min、90min)時,迅速放下加熱套,開始排氣,放鍋后濾出部分廢液,放冰箱保存待測。然后將漿料洗滌干凈,風干備用。1.3檢測和分析紙漿卡伯值及廢液殘堿的檢測方法見參考文獻。2結果與討論2.1氧壓和用量對脫木素速率的影響圖1為氧脫木素過程中,不同溫度、氧壓及用堿量下,紙漿卡伯值(正比與脫木素的程度)隨反應時間的變化規(guī)律。由圖1可知,在不同反應條件下紙漿卡伯值隨反應時間的增加均呈現(xiàn)先較快降低而后慢速下降的趨勢;而且在影響脫木素速率的3個工藝參數(shù)中,溫度對脫木素速率影響尤為顯著,表現(xiàn)為一定時間內(nèi)木素脫除率隨溫度的升高幾乎成比例增大(氧脫木素90min時,3種溫度的木素脫除率分別為85℃(32%)、100℃(45%)、115℃(56%)),而氧壓和用堿量增大至一定值(0.6MPa,3.0%)時,木素脫除率增大不再明顯。在相同反應時間下脫木素程度的上升量隨著氧壓和NaOH用量的增加逐漸下降。上述這些實驗數(shù)據(jù)為下面的建模提供了依據(jù)。2.2不同初始堿量氧脫木素含量的堿樣的擬合氧脫木素是利用堿性條件下的O2對紙漿中木素進行氧化脫除的反應,因此,反應途徑可由式(1)表示。L+O2+OH-→Ld(1)式中,L為紙漿中的殘余木素;Ld為斷裂的木素小分子。Ji等人在連續(xù)攪拌微分反應器(CSTR)中的理想狀態(tài)下氧脫木素的結果表明,氧脫木素反應在溫度、氧壓、堿濃不變的情況下,紙漿殘余木素的脫除為一級反應。根據(jù)化學反應原理及前人的研究結果,氧脫木素反應速率常數(shù)與溫度的關系可用阿倫尼烏斯方程表示、與氧壓、[OH-]的關系為冪指數(shù)關系。據(jù)此,氧脫木素的動力學方程可用式(2)表達。式中,CL為紙漿殘余木素含量,g/100g漿;t為氧脫木素時間,min;A為阿倫民烏斯方程指前因子;Ea為脫木素反應的活化能,kJ/mol;R為摩爾氣體常量,J/(mol·K);T為絕對溫度,K,P為氧壓,MPa;[OH-]為堿濃度,g/L;n為氧壓的反應級數(shù);m為堿濃度的反應級數(shù)。由于卡伯值與紙漿殘余木素之間呈正比關系,即:CL=r·K,則可將式(2)轉化為式(3)或式(4)。式中,K為紙漿卡伯值;r為紙漿卡伯值與木素含量的指標系數(shù),g/100g漿。其中,;在給定的溫度及氧壓條件下,k不隨反應時間而變。圖2為氧脫木素實驗中所測定過程溶液中的殘余NaOH濃度隨卡伯值自然對數(shù)(lnK)計算值的變化規(guī)律。由圖2可知,在給定不同的初始用堿量下,殘余NaOH濃度與ln(K)之間均呈現(xiàn)表觀截距不同而斜率相同的線性關系,而反應溫度和氧壓不同時,殘余NaOH濃度與ln(K)呈現(xiàn)近乎一致的線性關系,即溫度和氧壓的影響極不明顯。據(jù)此,殘余NaOH濃度與ln(K)之間可用式(5)表示。[OH-]=a+bln(K)(5)式中,a為堿濃度方程的截距;b為堿濃度方程的斜率。令不同用堿量下氧脫木素過程溶液的殘余NaOH濃度與ln(K)之間線性關系的斜率一致,而截距a則不同,并利用殘余NaOH濃度的測定值與ln(K)的計算值進行線性擬合得到不同NaOH用量下的a值,結果見表2;將表2中的a值對初始堿濃做線性擬合,得到式(6)。a=-8.73+0.865[OH-]0(R2=0.99)(6)合并式(5)和式(6)得到式(7)。[OH-]=-8.73+0.865×[OH-]0+2.342×ln(K)(7)將式(7)代入式(5),積分可得到紙漿卡伯值的數(shù)學模型見式(8)。2.3模型的評估和比較2.3.1堿濃、氧壓的反應動力學模型將本實驗不同反應條件及其對應的卡伯值數(shù)據(jù)代入式(8),利用ExcelSolver擬合出該模型的各個參數(shù),結果如表3所示。由表3可知,本實驗中的氧脫木素動力學模型對堿濃、氧壓的反應級數(shù)及氧脫木素活化能的擬合值分別為0.372、0.316及54.50kJ/mol,這與Ji等人用CSTR反應器理想體系下氧脫木素動力學模型的研究結果(分別為0.412,0.305和54.5kJ/mol)相近。圖3是本實驗數(shù)據(jù)對模型的擬合的相關結果。由圖3可知,該模型能很好地預測氧脫木素過程中不同工藝條件下的紙漿卡伯值(R2=0.980),這將為桉木硫酸鹽漿氧脫木素工藝條件的設定及控制提供理論的指導。2.3.2本實驗模型參數(shù)設置的合理性分析表4為兩種經(jīng)典的氧脫木素模型(多級脫木素模型、兩相脫木素模型)與本實驗建立的動力學模型參數(shù)及擬合效果的對比。由表4可知,盡管這兩種經(jīng)典模型也能對本實驗數(shù)據(jù)進行很好的擬合;但是,本實驗建立的動力學模型的參數(shù)量顯然比上述兩種模型的參數(shù)要少。另外,在利用上述兩種模型進行擬合過程中發(fā)現(xiàn),卡伯值的反應級數(shù)在0.5~5的范圍內(nèi)確定任意數(shù)值均能對該模型進行很好的擬合,這說明卡伯值的反應級數(shù)并不能通過多級脫木素模型確定,從而降低了該模型的可靠性;而兩相脫木素模型的多參數(shù)使得對本實驗數(shù)據(jù)進行擬合的過程中,模型的每個參數(shù)的可變性都很大。因此,從模型參數(shù)設置的合理性及參數(shù)量的控制而言,本實