李文軍 劉婷 讀剛 王偉

【摘 要】利用Matlab的GUI功能，開發(fā)了流化床煤氣化計(jì)算系統(tǒng)。該圖形用戶界面系統(tǒng)可以根據(jù)建立的煤炭氣化過(guò)程產(chǎn)品的平衡組成計(jì)算模型，方便、高效地通過(guò)改變相關(guān)參數(shù)模擬出煤炭氣化平衡組成的重要數(shù)據(jù)。實(shí)例表明：該流化床煤氣化計(jì)算系統(tǒng)可以輔助研究氣化操作參數(shù)對(duì)氣化組成的影響規(guī)律，對(duì)生產(chǎn)中有針對(duì)性的控制操作參數(shù)，調(diào)節(jié)煤炭氣化產(chǎn)品組成和相對(duì)含量具有實(shí)際指導(dǎo)意義。

【關(guān)鍵詞】Matlab；流化床；煤氣化；模擬仿真

中圖分類號(hào)： TP391 文獻(xiàn)標(biāo)識(shí)碼： A 文章編號(hào)： 2095-2457（2018）28-0260-003

DOI：10.19694/j.cnki.issn2095-2457.2018.28.122

【Abstract】The calculation system of fluidized bed coal gasification is developed with the GUI function of Matlab. The graphical user interface system can calculate the balance composition of coal gasification process products according to the established model， conveniently and efficiently simulate the important data of coal gasification balance composition by changing the relevant parameters. The example shows that the fluidized bed gasification calculation system can assist in studying the influence law of gasification operation parameters on gasification composition， and has practical guiding significance for controlling operation parameters in production and regulating the composition and relative content of coal gasification products.

【Key words】Matlab； Fluidized bed； Gasification； Simulation

0 引言

MATLAB是美國(guó)The MathWorks公司出品的一款商業(yè)數(shù)學(xué)軟件。它具有出色的數(shù)值計(jì)算能力，并在世界數(shù)值計(jì)算軟件中占據(jù)主導(dǎo)地位。數(shù)值計(jì)算的功能包括變量和數(shù)值、數(shù)值微積分、求解線性方程組、矩陣和數(shù)組、多項(xiàng)式計(jì)算、數(shù)據(jù)分析等。除此以外，MATLAB可以調(diào)用其它語(yǔ)言編寫的程序，比如有C語(yǔ)言程序，C++和FORTRAN等。我們可以用MATLAB創(chuàng)建用戶界面，執(zhí)行矩陣運(yùn)算、繪制各種函數(shù)并根據(jù)需要分析數(shù)據(jù)，還可以用來(lái)連接其他編程語(yǔ)言的程序。MATLAB是一種功能強(qiáng)大的軟件，在廣泛應(yīng)用在工程軟件中。

流化床具有類似流體特性，較輕的物體可以浮在流化床上面；當(dāng)容器傾斜時(shí)，床層頂部保持水平；當(dāng)兩個(gè)流化床層連通時(shí)，其床面自動(dòng)找平；固體顆粒可以從一個(gè)容器流道另一個(gè)容器中；床層中任意兩點(diǎn)的壓力差等于兩點(diǎn)間床層的靜壓頭。因此，流化床技術(shù)有許多優(yōu)點(diǎn)：

（1）床層溫度均勻，顆粒流動(dòng)平穩(wěn)，操作可連續(xù)自控；

（2）顆粒混合迅速，反應(yīng)器處于等溫狀態(tài)，操作簡(jiǎn)單，控制可靠；

（3）顆粒能在兩個(gè)反應(yīng)器之間循環(huán)流動(dòng)，傳熱方便；

（4）適宜于大規(guī)模操作；

（5）氣體與固體之間的高傳熱和傳質(zhì)速率；

（6）流化床與其中浸沒(méi)的構(gòu)建之間的傳熱速率高，所需換熱面積小。

由于流化床中的固體顆粒具有上述性質(zhì)，使得固體能像流體一樣進(jìn)行各種處理，流體化技術(shù)在工業(yè)上有許多應(yīng)用。目前，流化床技術(shù)在煤炭轉(zhuǎn)化過(guò)程已有大量應(yīng)用，如煤的氣化與干餾、費(fèi)-托合成、乙烯氧化、石油的催化裂化等。簡(jiǎn)而言之，流化技術(shù)已經(jīng)是一種具有廣泛應(yīng)用和成熟工業(yè)應(yīng)用的技術(shù)。

本文主要用MATLAB對(duì)煤氣化過(guò)程平衡進(jìn)行編程，分析反應(yīng)溫度、氧氣濃度、煤焦含碳量、水蒸氣量和氧氣預(yù)熱溫度對(duì)氣化體系熱力學(xué)平衡組成的影響，以及干燃?xì)鉄嶂怠袢細(xì)鉄嶂岛托实淖兓闆r，并進(jìn)行數(shù)值解模擬，得到數(shù)值解數(shù)據(jù)，為實(shí)際生產(chǎn)提供初步的理論指導(dǎo)和參考。

1 流化床煤氣化計(jì)算系統(tǒng)

1.1 設(shè)計(jì)原則

本文利用MATLAB的強(qiáng)大的數(shù)值計(jì)算功能，對(duì)煤焦炭氣化過(guò)程組分含量計(jì)算進(jìn)行了編程模擬，軟件系統(tǒng)采用適宜的模型和先進(jìn)的算法，使用模塊化程序結(jié)構(gòu)便于系統(tǒng)的擴(kuò)充，通過(guò)GUI設(shè)計(jì)的可視化圖形界面，讓用戶使用人機(jī)會(huì)話的輸入方式，通過(guò)界面上的按鈕和菜單，根據(jù)顯示屏的提示，從鍵盤輸入必要的信息，即可完成一系列復(fù)雜的計(jì)算并按指定方式輸出結(jié)果，不僅可以完成預(yù)定的任務(wù)，也可以根據(jù)實(shí)際需要自行設(shè)計(jì)實(shí)驗(yàn)內(nèi)容、修改實(shí)驗(yàn)參數(shù)來(lái)獲取滿足自己需要的實(shí)驗(yàn)結(jié)果，充分體現(xiàn)了系統(tǒng)的可操作性和交互性原則。

1.1.1 可靠性

軟件系統(tǒng)的可靠性直接關(guān)系到設(shè)計(jì)自身的聲譽(yù)，生存和發(fā)展競(jìng)爭(zhēng)力。在運(yùn)行過(guò)程中，該系統(tǒng)可以抵抗異常情況的干擾，保證系統(tǒng)正常運(yùn)行，功能平穩(wěn)。

1.1.2 健壯性

軟件健壯性是一個(gè)比較籠統(tǒng)的概念，但它是軟件質(zhì)量的一個(gè)非常重要的衡量標(biāo)準(zhǔn)。該系統(tǒng)能夠判定出不符合規(guī)范要求的輸入，給出相應(yīng)的處理方式。

1.1.3 可修改性

該系統(tǒng)有良好的結(jié)構(gòu)和完備的文檔，便于系統(tǒng)之后的調(diào)整。

1.1.4 可讀性

軟件的可讀性是其可靠性和可修改性的前提。該系統(tǒng)文檔清晰可讀，系統(tǒng)本身結(jié)構(gòu)簡(jiǎn)單明了。

1.1.5 實(shí)用易用性

該系統(tǒng)條理清晰，各個(gè)功能操作簡(jiǎn)單、快捷，易于使用，具有良好的可操作性，能夠解決實(shí)際問(wèn)題。

1.2 系統(tǒng)結(jié)構(gòu)框架

整個(gè)系統(tǒng)可以分為三個(gè)階段，預(yù)處理、數(shù)據(jù)輸入和計(jì)算輸出。各階段的相互關(guān)系如圖1所示。

預(yù)處理部分根據(jù)具體的任務(wù)完成反應(yīng)方程的編輯、編譯和執(zhí)行程序的生成；輸入部分負(fù)責(zé)數(shù)據(jù)輸入和存儲(chǔ)；計(jì)算輸出部分負(fù)責(zé)整個(gè)系統(tǒng)的計(jì)算和輸出。系統(tǒng)主界面如下圖所示。

2 應(yīng)用案例

為了展示如何使用Matlab GUI構(gòu)建流化床煤氣化計(jì)算系統(tǒng)，本文通過(guò)具體實(shí)例說(shuō)明了系統(tǒng)的設(shè)計(jì)過(guò)程、GUI圖形界面的效果及實(shí)驗(yàn)結(jié)果。

煤的氣化反應(yīng)是由熱解產(chǎn)生的揮發(fā)分、殘余焦炭顆粒和氣化劑的復(fù)雜反應(yīng)。本文以煤炭氣化過(guò)程為示例進(jìn)行模擬，生成的氣體主要是煤焦炭和氣化劑反應(yīng)得到，其主要產(chǎn)物是：可燃性氣體CO和H2，小部分碳氧化生成CO2，另外還會(huì)有少量的水生成和體系中的氮?dú)饨M分，本文以這五種組分為計(jì)算的主要對(duì)象，研究不同因素對(duì)這些組分的含量的影響。

在氣化過(guò)程以100kg焦炭和2kmol水蒸氣為計(jì)算基準(zhǔn)的條件下，假設(shè)煤炭焦中含碳量為62.25%，其余為灰分，碳轉(zhuǎn)化率為96%，氧氣摩爾分?jǐn)?shù)為21%，氧氣預(yù)熱溫度0℃，反應(yīng)壓力為0.1MPa，散熱損失為5%。

在實(shí)驗(yàn)中，分別改變氣化反應(yīng)溫度、氧氣濃度、煤炭焦含碳量、水蒸氣量及空氣預(yù)熱溫度，模擬得出5組氧氣消耗量、燃?xì)饨M分的摩爾質(zhì)量，總?cè)細(xì)饬?，濕燃?xì)鉄嶂岛透扇細(xì)鉄嶂狄约癊x效率的計(jì)算結(jié)果。

根據(jù)計(jì)算結(jié)果，分別以反應(yīng)溫度、氧氣濃度、煤炭焦含碳量、水蒸氣量及空氣預(yù)熱溫度為橫坐標(biāo)，CO、CO2、H2、H2O、N2、氧氣消耗量、濕燃?xì)饪倸饬?、濕燃?xì)鉄嶂?、干燃?xì)鉄嶂?、效率為縱坐標(biāo)，得到5組各組分摩爾含量隨著橫坐標(biāo)的曲線圖，如圖3-7所示。

從圖3（a）看，氣化平衡過(guò)程在一定溫度范圍內(nèi)隨著溫度的提高，平衡組分中N2、氧氣消耗量和濕燃?xì)饪倸饬侩S著增多，組分中CO和H2的摩爾含量略有下降，CO2和H2O摩爾含量略有上升，從圖3（b）看，升高溫度，濕燃?xì)鉄嶂岛透扇細(xì)鉄嶂到档洼^快，效率降低；

從圖4（a）看，氣化平衡過(guò)程在一定溫度范圍內(nèi)隨著氧氣濃度的提高，平衡組分中N2、氧氣消耗量和濕燃?xì)饪倸饬恐饾u減少，CO含量增加，組分中H2的摩爾含量略有增加，CO2和H2O摩爾含量略有基本相對(duì)穩(wěn)定，從圖4（b）看，升高溫度，濕燃?xì)鉄嶂岛透扇細(xì)鉄嶂翟黾虞^快，效率有所提高，在生產(chǎn)過(guò)程中提高氧氣濃度對(duì)增加燃?xì)鉄嶂涤兄苯佑行У囊嫣?，但是在反?yīng)體系中氧氣濃度高，工藝條件要求就比較苛刻，且操作危險(xiǎn)性較大，較難實(shí)現(xiàn)高氧氣濃度進(jìn)料。因此，可以參考現(xiàn)有的工藝技術(shù)選取氧氣濃度操作參數(shù)；

由圖5（a）看，可以了解到氣化過(guò)程隨著煤焦含碳量的增加，氧氣耗量、氮?dú)?、一氧化碳摩爾含量和濕燃?xì)饪倸饬吭黾用黠@，其他組分相對(duì)變化較??；有圖5（b）看，溫度升高，濕燃?xì)夂透扇細(xì)鉄嶂翟黾樱视兴岣撸?/p>

由圖6（a）看，可以了解到氣化過(guò)程隨著水蒸氣含量的增加，氧氣耗量變化不大，組分中氮?dú)饬咳匀惠^大、濕燃?xì)饪倸饬可杂性黾?，其他組分相對(duì)變化較??；由圖6（b）看，水蒸氣含量增多，濕燃?xì)夂透扇細(xì)鉄嶂刀加兴陆?，效率降低?/p>

由圖7（a）看，隨著氧氣預(yù)熱溫度的逐漸升高，氧氣耗量、組分氮?dú)夂?，濕燃?xì)饪偭恐饾u減少，一氧化碳的含量逐漸略有增多，二氧化碳的有所減少，從圖7（b）分析預(yù)熱溫度升高，濕燃?xì)鉄嶂岛透扇細(xì)鉄嶂抵饾u提升，效率略有提高。

通過(guò)以上模擬研究氣化溫度、氧氣濃度、煤炭焦含碳量、水蒸氣量及、氧氣預(yù)熱溫度的影響規(guī)律，可在一定程度上了解氣化組分含量隨著各個(gè)因素改變的變化趨勢(shì)，可以為生產(chǎn)工藝條件的改進(jìn)作參考。上述因素的影響，具有一定的條件限制，只有結(jié)合實(shí)際的生產(chǎn)條件，通過(guò)經(jīng)濟(jì)學(xué)的比較，才能確定各種參數(shù)的取舍。

3 結(jié)論

流化床煤氣化計(jì)算系統(tǒng)可以輔助研究氣化操作參數(shù)對(duì)氣化組成的影響規(guī)律，對(duì)生產(chǎn)中有針對(duì)性的控制操作參數(shù)，調(diào)節(jié)煤炭氣化產(chǎn)品組成和相對(duì)含量具有指導(dǎo)性的實(shí)際意義。同時(shí)對(duì)使用Matlab軟件進(jìn)行建模能力的提升起著積極的推動(dòng)作用。

【參考文獻(xiàn)】

[1]朱自強(qiáng)，吳有庭.化工熱力學(xué)[M].北京：化學(xué)工程出版社，2013.9.

[2]薛定宇，陳陽(yáng)泉.基于MATLAB/Simulink的系統(tǒng)仿真技術(shù)與應(yīng)用[M].北京：清華大學(xué)出版社，2002.

[3]楊小剛，趙穎.煤氣化過(guò)程熱力學(xué)平衡模型研究進(jìn)展[J].化學(xué)工程與裝備，2012，12：146-149.