甲醇制丙烯工藝的再生過(guò)程模擬與分析

梁 昊， 王 林， 廖祖維， 雍曉靜， 蔣斌波， 王靖岱， 陽(yáng)永榮， 焦洪橋

(1.浙江大學(xué) 化學(xué)工程與生物工程學(xué)學(xué)院 化學(xué)工程聯(lián)合國(guó)家重點(diǎn)實(shí)驗(yàn)室， 浙江 杭州 310027；2.神華寧夏煤業(yè)集團(tuán)， 寧夏 銀川 750411)

?

甲醇制丙烯工藝的再生過(guò)程模擬與分析

梁 昊1， 王 林2， 廖祖維1， 雍曉靜2， 蔣斌波1， 王靖岱1， 陽(yáng)永榮1， 焦洪橋2

(1.浙江大學(xué) 化學(xué)工程與生物工程學(xué)學(xué)院 化學(xué)工程聯(lián)合國(guó)家重點(diǎn)實(shí)驗(yàn)室， 浙江 杭州 310027；2.神華寧夏煤業(yè)集團(tuán)， 寧夏 銀川 750411)

建立了甲醇制丙烯(MTP)固定床再生器軸向一維動(dòng)態(tài)模型，考察了氧體積分?jǐn)?shù)、入口溫度、再生氣流量對(duì)再生效果的影響，以及再生器內(nèi)溫度及炭含量的分布情況。結(jié)果表明，影響再生效果的諸因素中，作用由大到小依次為氧體積分?jǐn)?shù)、入口溫度、再生氣流量。提出“提高再生氣量”的再生方案，可有效減小新鮮氮?dú)獾挠昧浚s短再生時(shí)間，對(duì)工業(yè)操作有一定的指導(dǎo)意義。

甲醇制丙烯(MTP)； 固定床； 反應(yīng)器； 再生器

近年來(lái)，甲醇制烯烴(MTO)和甲醇制丙烯(MTP)工藝迅速發(fā)展，已成為我國(guó)烯烴生產(chǎn)的重要工藝。但該類(lèi)工藝高能耗的特點(diǎn)也制約了其進(jìn)一步發(fā)展。例如，ZSM-5分子篩催化MTP反應(yīng)及其再生過(guò)程的能耗均高于同類(lèi)型石化裝置，故亟待改進(jìn)、提升。

有大量的研究關(guān)注MTP工藝的反應(yīng)過(guò)程，但關(guān)于催化劑再生的研究還不多見(jiàn)。Bibby等[1]認(rèn)為，甲醇制汽油(MTG)工藝中ZSM-5催化劑的燒炭速率比其他工藝中的ZSM-5催化劑的燒炭速率要慢。Andre等[2]擬合得出ZSM-5在MTG過(guò)程后的燒炭動(dòng)力學(xué)方程及再生催化劑的活性恢復(fù)動(dòng)力學(xué)方程。Li等[3]認(rèn)為，ZSM-5催化劑的積炭類(lèi)型可分為兩種，500℃以下為可溶性結(jié)焦，500℃以上為積炭；溫鵬宇等[4]認(rèn)為，可溶性結(jié)焦主要為三甲基苯，其可遷移到催化劑外表面石墨化形成積炭，從而堵塞分子篩孔道。

MTP工藝中，催化劑在再生器中的再生過(guò)程存在著再生時(shí)間長(zhǎng)、氮?dú)庀拇蟆⒛芎母叩膯?wèn)題。雖然存在著較大的改進(jìn)空間，但限于對(duì)其再生規(guī)律的認(rèn)識(shí)不足，相關(guān)工作尚未展開(kāi)。筆者結(jié)合再生動(dòng)力學(xué)與固定床再生器模型，對(duì)實(shí)際工況進(jìn)行動(dòng)態(tài)模擬，獲得催化劑再生規(guī)律的初步認(rèn)識(shí)，并進(jìn)一步考察再生氣氧體積分?jǐn)?shù)、入口溫度、再生氣流量等操作參數(shù)對(duì)再生效果的影響，優(yōu)化了催化劑再生條件。

1 積炭催化劑再生動(dòng)力學(xué)

ZSM-5分子篩的積炭實(shí)際上是高度縮合的碳?xì)浠衔铮銫/H質(zhì)量比在10～20之間，失活后的催化劑炭含量接近18%[4]。一般認(rèn)為，積炭催化劑再生動(dòng)力學(xué)方程rC符合冪函數(shù)形式，如式(1)所示， ZSM-5分子篩催化劑積炭的再生動(dòng)力學(xué)方程[2]如式(2)所示。

(1)

(2)

假設(shè)碳元素燃燒全部生成CO2，氫元素燃燒全部生成H2O，則催化劑表面的燒炭反應(yīng)可用式(3)表示。

CH12/α+(1+3/α)O2→CO2+(6/α)H2O

(3)

所以，當(dāng)催化劑表面的燒炭速率為rC時(shí)，耗氧速率rO2的表達(dá)式為式(4)。

(4)

2 MTP工藝再生器再生過(guò)程模型的建立

假定再生器內(nèi)氣體的溫度和體積分?jǐn)?shù)沿徑向無(wú)變化，沿軸向有變化，忽略床層壓降并認(rèn)為催化劑和再生氣的熱容保持不變，可用一維擬均相非穩(wěn)態(tài)模型描述。

再生器床層中沿軸向取長(zhǎng)度為dz的一維微元建立模型[7]，如圖1所示。氧平衡、碳平衡和熱量平衡方程分別見(jiàn)式(5)～(7)。

(5)

圖1 再生器微元示意圖Fig.1 The schematic of micro unit in regeneration reactor

(6)

Fcp,g(T+273)-(Fcp,g(T+273)+

rcρBπr2dzΔH+2πrdzλ(Τ-Τw)=

(7)

參考劉偉等[9]的工作，通過(guò)MATLAB軟件進(jìn)行數(shù)值求解，即可獲得再生器內(nèi)溫度分布、炭含量分布及氧體積分?jǐn)?shù)分布等關(guān)鍵信息。

3 MTP工藝再生模型的工業(yè)驗(yàn)證

3.1 工業(yè)再生過(guò)程簡(jiǎn)介

MTP反應(yīng)和再生過(guò)程均在同一反應(yīng)器內(nèi)完成，當(dāng)進(jìn)料甲醇轉(zhuǎn)化率低于90%時(shí)，反應(yīng)器須切換再生[5]。為保護(hù)催化劑，再生過(guò)程中需控制再生溫度在(490±2)℃以下。為控制以上操作條件，實(shí)際操作中采用大氣體流量下逐步提升氧體積分?jǐn)?shù)和入口氣體溫度的策略。再生器床層高度依次為0.253、0.286、0.327、0.384、0.462和0.588 m，其他操作參數(shù)及完整的再生程序如表1、表2所示。

表1 MTP工藝模擬條件的設(shè)定Table 1 The simulated conditions for MTP process

表2 MTP工藝催化劑再生過(guò)程中新鮮N2和空氣流量及溫度控制Table 2 The N2, air flows and temperature control of catalyst regeneration in MTP process

3.2 模擬結(jié)果與工業(yè)數(shù)據(jù)的比較

筆者采用的數(shù)據(jù)來(lái)自工業(yè)裝置，包括再生時(shí)新鮮空氣流量、溫度、再生器出口物料濃度等信息，分別如圖2、圖3、表3所示。

圖2 MTP工藝中新鮮空氣流量隨再生時(shí)間的變化Fig.2 Fresh air flow vs regeneration time in MTP process

圖3 MTP工藝中床層溫度隨再生時(shí)間的變化Fig.3 Bed temperature vs regeneration time in MTP process Bed 1 inlet; Bed 4 inlet; Bed 5 inlet; Bed 6 inlet; Bed 4 exit; Bed 5 exit; Bed 6 exit

從圖2可以看出，提氧過(guò)程主要分為4個(gè)階段，所對(duì)應(yīng)的空氣流量約為5300、5800、6400、7600 m3/h。由圖3可知，第1床層再生起始溫度為445℃，是提氧階段，持續(xù)時(shí)間約3 d；起始溫度為460℃，是第6床層燒炭階段，持續(xù)時(shí)間約1 d；起始溫度為490℃，是恒溫階段，持續(xù)時(shí)間約2 d。同時(shí)，圖3還列出了幾條容易發(fā)生飛溫的床層溫度曲線，分別是第4、第5和第6床層，溫度最高可達(dá)490℃，說(shuō)明實(shí)際操作過(guò)程中保持穩(wěn)定狀態(tài)較為困難。

表3 MTP工藝中再生期間廢氣組成Table 3 The composition of waste gas during the regeneration period in MTP process

模擬結(jié)果主要與工業(yè)結(jié)果對(duì)比了床層溫度及出口氧體積分?jǐn)?shù)，如圖4、表4所示。

圖4 MTP工藝中再生器第1、第6床層出口溫度模擬結(jié)果與工業(yè)結(jié)果對(duì)比Fig.4 The comparison of simulated and industrial results of the temperature at bed 1 and bed 6 in MTP process(a) Bed 1; (b) Bed 6(1)Simulation value; (2) Industrial value表4 MTP工藝中再生器出口O2體積分?jǐn)?shù)模擬結(jié)果與 工業(yè)結(jié)果對(duì)比Table 4 The comparison of simulated with industrial results of O2 volume fraction at reactor exit in MTP process

Stepφ(O2)atbedoutlet/%SimulatedvalueIndustrialvalue10.210.4120.510.7431.250.9541.581.5051.88-62.932.90720.3220.57

由上述結(jié)果可知，溫度的模擬結(jié)果與工業(yè)值較吻合，O2體積分?jǐn)?shù)的模擬結(jié)果與工業(yè)值總體變化趨勢(shì)一致，在O2體積分?jǐn)?shù)高于1.5%時(shí)誤差低于5%。所建立的再生模型可定性分析操作條件對(duì)再生結(jié)果的影響。

4 操作條件對(duì)MTP工藝中再生器再生結(jié)果的影響

按表1、表2給出的條件進(jìn)行模擬，可得到整個(gè)床層溫度、積炭量隨再生時(shí)間的變化情況，結(jié)果示于圖5。

圖5 MTP工藝中不同床層厚度下再生器床層溫度(T)及炭含量(wC)隨再生時(shí)間的變化Fig. 5 T and wC vs regeneration time with different bed thickness in MTP process(a) Bed temperature(T)； (b) Coke mass fraction(wC)

由于再生時(shí)采用逐步提升O2體積分?jǐn)?shù)與入口溫度的操作方式，故圖5(a)中的溫度突變?yōu)樵偕僮髑袚Q到下一個(gè)操作步驟引起。

由圖5(a)可知，(1)同一床層溫度隨再生時(shí)間先增加后減小，這是因?yàn)槠鹗茧A段溫度隨反應(yīng)熱量的升高而增加，隨著反應(yīng)的進(jìn)行，積炭逐漸減少，溫度隨反應(yīng)熱量的降低而減小；(2)隨著床層數(shù)增加，其溫度峰值會(huì)延后出現(xiàn)，這是因?yàn)闅怏w每通過(guò)一個(gè)床層，都會(huì)消耗一部分O2，導(dǎo)致下一床層O2體積分?jǐn)?shù)隨時(shí)間延長(zhǎng)而增加，使溫度上升，峰值延后出現(xiàn)。由圖5(b)可知，前面床層的燒炭速率要大于后面床層，表現(xiàn)為其積炭量由紅色變?yōu)樗{(lán)色所需時(shí)間更短。綜合圖5(a)和(b)可知，一開(kāi)始第5、第6床層由于O2體積分?jǐn)?shù)非常低，不能充分利用前面床層燃燒所釋放的熱量進(jìn)行燒炭反應(yīng)，導(dǎo)致再生時(shí)6段床層并不是均勻燃燒，可見(jiàn)再生改進(jìn)空間巨大。

將再生器出口處催化劑炭質(zhì)量分?jǐn)?shù)wC≤0.01%作為再生結(jié)束的判據(jù)，并要求床層最高溫度低于(490±2)℃，以保證催化劑骨架不會(huì)在水熱環(huán)境下遭到破壞。由圖3可以看出，第6床層對(duì)溫度最為敏感，故以下所有分析均以第6床層為對(duì)象。

4.1 O2體積分?jǐn)?shù)的影響

按照表2的操作條件，將O2體積分?jǐn)?shù)從0.4%提高到1.0%，觀察床層溫度和燒炭轉(zhuǎn)化率隨再生時(shí)間的變化，結(jié)果示于圖6。

圖6 MTP工藝中不同O2體積分?jǐn)?shù)(φ(O2))和 燒炭轉(zhuǎn)化率(x(Coke))下第6床層溫度隨再生時(shí)間的變化Fig.6 Temperature at bed 6 vs regeneration time with different φ(O2) and x(Coke) in MTP process φ(O2)=0.4%, x(Coke)=76.31%; φ(O2)=0.5%, x(Coke)=85.20%; φ(O2)=0.6%, x(Coke)=90.89%; φ(O2)=0.7%, x(Coke)=94.46%; φ(O2)=0.8%, x(Coke)=96.68%; φ(O2)=0.9%, x(Coke)=98.01%; φ(O2)=1.0%, x(Coke)=98.76%

從圖6可以看出，床層溫度隨著入口O2體積分?jǐn)?shù)增加而迅速上升。在開(kāi)始階段，若入口O2體積分?jǐn)?shù)高于0.6%，則床層最高溫度迅速超過(guò)490℃，可見(jiàn)該變量較敏感。此外，雖然提高O2體積分?jǐn)?shù)能使積炭轉(zhuǎn)化率迅速下降，但不可避免會(huì)產(chǎn)生高溫，導(dǎo)致催化劑骨架遭到不可逆轉(zhuǎn)的破壞。因此，不可能通過(guò)單純改變?nèi)肟贠2體積分?jǐn)?shù)來(lái)同時(shí)滿(mǎn)足再生過(guò)程中溫度和再生效果的雙重要求，還需要調(diào)節(jié)其他參數(shù)。

4.2 再生氣入口溫度的影響

按照表2的操作條件，僅將再生氣入口溫度從445℃提高到470℃，觀察第6床層溫度及出口積炭量的變化，結(jié)果示于圖7和表5。

圖7 MTP工藝中不同入口溫度下第6床層 溫度隨再生時(shí)間的變化Fig.7 Temperature at bed 6 vs regeneration time with different inlet temperatures in MTP process Inlet temperatures/℃: 445; 450; 455; 460; 465; 470表5 MTP工藝中不同入口溫度下第6床層 出口處的積炭轉(zhuǎn)化率Table 5 Coke conversion at bed 6 with different inlet temperatures in MTP process

Inlettemperature/℃x(Coke)/%Thehighestbedtemperature/℃ΔT1)/℃44567.794551045072.774631345577.414701546081.614802046585.324902547088.5050030

1) The difference between the highest bed temperature and inlet temperature

從圖7和表5可以看出，對(duì)于較低的入口溫度(445～460℃)，床層最高溫度提升約10～15℃，當(dāng)入口溫度大于460℃時(shí)，床層最高溫度可提高20～30℃。低溫下提高10～15℃，高溫下提高20～30℃，這可能是因?yàn)槌跏荚偕鷾囟容^低，并未達(dá)到再生反應(yīng)最佳溫度所致。為控制溫度低于490℃，較佳的入口溫度為465℃。

4.3 再生氣入口流量的影響

從動(dòng)力學(xué)表達(dá)式看，O2體積分?jǐn)?shù)及溫度都通過(guò)影響燒炭反應(yīng)速率來(lái)影響催化劑再生效果，而再生氣入口流量通過(guò)改變氣體移熱能力來(lái)影響再生效果。圖8為不同再生氣流量下第6床層溫度隨再生時(shí)間的變化。

圖8 MTP工藝中不同再生氣流量下第6床層 溫度隨再生時(shí)間的變化Fig.8 Temperature at bed 6 vs regeneration time with different regeneration gas flow rate in MTP process Regeneration gas flow rate/(m3·h-1): 60000; 70000; 80000; 90000; 100000

從圖8可以看出，提高再生氣流量，可以降低床層溫度，但隨著氣體流量的增加，床層溫度的下降幅度越來(lái)越小。因此，再生氣入口流量不應(yīng)過(guò)小，以防止最高溫度超過(guò)490℃；入口流量也不應(yīng)過(guò)大，以防止床層溫度過(guò)低，達(dá)不到最佳反應(yīng)溫度。在表2 操作條件下，較佳的再生氣入口流量是80000 m3/h。

4.4 燒炭再生方案

針對(duì)固定床MTP 6段再生器，在保證床層最高溫度不超過(guò)(490±2)℃的前提下，改變?cè)偕鷼饪倸饬?其中循環(huán)比為0.8)和新鮮空氣量，以整個(gè)床層積炭量wC≤0.01%時(shí)所消耗的新鮮N2較小為目標(biāo)，提出了如表6所列8種再生方案。

從表6可以看出，隨著再生氣量的增加，所需的新鮮N2量越來(lái)越小，最后基本不變。這是因?yàn)闅怏w移熱能力隨氣量增加而提高，此時(shí)可適當(dāng)提高初始O2體積分?jǐn)?shù)，即新鮮空氣的增加幅度可大于再生氣量的增加幅度，故所需的N2量會(huì)減小。溫度對(duì)O2體積分?jǐn)?shù)極其敏感，需嚴(yán)格控制。在臨界值附近，新鮮空氣與再生氣量的增加幅度需持平，此時(shí)所需新鮮N2量基本不變。此外，從表6可知較優(yōu)的再生氣量為120000 m3/h左右，此時(shí)需要的再生時(shí)間為2 d，新鮮N2用量為591500 m3，該值比目前工業(yè)操作所需的新鮮氮?dú)庥昧抗?jié)省約30%。因此，提高再生氣總氣量可有效減小新鮮N2的消耗量，縮短再生時(shí)間，這對(duì)于工業(yè)操作有一定的指導(dǎo)意義。

表6 8種MTP再生器再生方案的模擬結(jié)果Table 6 The simulation results of several regeneration schemes for MTP regeneration reactors

5 結(jié) 論

(1)在影響甲醇制丙烯工藝再生器再生效果的諸因素中，作用由大到小依次為O2體積分?jǐn)?shù)、入口溫度和再生氣流量。由于增大O2體積分?jǐn)?shù)和提高入口溫度都會(huì)升高燒炭溫度，二者均受到最高燒炭溫度的限制；增大入口流量可有效增大熱容，從而松弛對(duì)O2體積分?jǐn)?shù)、入口溫度的約束，有效加快燒炭速率。

(2)對(duì)于不同積炭量的床層，可通過(guò)改變?nèi)肟跍囟取2體積分?jǐn)?shù)、再生氣流量3個(gè)操作參數(shù)來(lái)優(yōu)化再生過(guò)程。如針對(duì)一開(kāi)始第5、第6床層燒炭慢的問(wèn)題，可通過(guò)側(cè)線引入新鮮空氣，提高O2體積分?jǐn)?shù)、加快其燒炭速率的方法，從而使6段床層同時(shí)充分燃燒。

(3)提高現(xiàn)有過(guò)程的總再生氣量可減小新鮮N2的消耗量，縮短再生時(shí)間，可望將再生能耗減少30%。

符號(hào)說(shuō)明：

cp,g——?dú)怏w熱容，kJ/(kg·K)；

cp,s——固體熱容，kJ/(kg·K)；

E——活化能，J/mol；

F——再生氣摩爾流量，mol/h；

FO2——氧氣摩爾流量，mol/h；

ΔH——燃燒熱，kJ/kg；

k0——指前因子，kg C/(kg cat·min·Pa)；

MC——碳元素摩爾質(zhì)量，kg/kmol；

Mg——氧氣摩爾質(zhì)量，kg/kmol；

nO2——氧氣物質(zhì)的量，mol；

pO2——氧氣分壓，×105Pa；

r——再生器半徑，m；

rC——燒炭速率，kg C/(kg cat·min)；

rO2——耗氧速率，kg O2/(kg cat·min)；

R——理想氣體常數(shù)，J/(mol·K)；

t——反應(yīng)時(shí)間，min；

T——床層溫度，℃；

Tw——壁面溫度，℃；

wC——催化劑炭質(zhì)量分?jǐn)?shù)，kg C/kg cat；

z——床層厚度，m；

ρB——催化劑堆密度，kg/m3；

ρg——?dú)怏w密度，kg/m3；

λ——總傳熱系數(shù)，kJ/(K·min·m2)；

α——積炭物質(zhì)的碳/氫質(zhì)量比；

ε——催化劑孔隙率。

[1] BIBBY D M, HOWE R F, MCLELLAN G D. Coke formation in high-silica zeolites[J].Applied Catalysis A General, 1992, 93(1): 1-34.

[2] ANDRES T, AGUAYO, ANA G, et al. Coke combustion and reactivation kinetics of a ZSM-5 zeolite based catalyst used for the transformation of methanol into hydrocarbons[J].Chemical Engineering Communications, 1999, 176(1): 43-63.

[3] LI J, XIONG G, FENG Z, et al. Coke formation during the methanol conversion to olefins in zeolites studied by UV Raman spectroscopy[J].Microporous & Mesoporous Materials, 2000, 39(1-2): 275-280.

[4] 溫鵬宇, 梅長(zhǎng)松, 劉紅星, 等. 甲醇制丙烯過(guò)程中ZSM-5催化劑的失活行為[J].石油學(xué)報(bào)(石油加工), 2008, 24(4): 446-450. (WEN Pengyu, MEI Changsong, LIU Hongxing, et al. Deactivation of ZSM-5 catalysts during methanol to propylene process[J].Acta Petrolei Sinica(Petroleum Processing Section), 2008, 24(4): 446-450.)

[5] 閆錫軍, 王志斌, 夏春江. 甲醇制丙烯(MTP)工藝過(guò)程中催化劑的再生過(guò)程分析[J].廣東化工, 2013, 40(2): 53-54.(YAN Xijun， WANG Zhibin， XIA Chunjian. Analysis of regeneration process for MTP catalyst[J].Guangdong Chemical Industry, 2013, 40(2): 53-54.)

[6] 朱玉川, 趙啟輝. 神華寧煤烯烴MTP反應(yīng)器再生優(yōu)化[J].中國(guó)高新技術(shù)企業(yè)旬刊, 2014, (9): 38-39. (ZHU Yuchuan， ZHAO Qihui. Optimization of MTP reactor regeneration Shenhua Ningxia coal olefins[J].China High-Tech Enterprises, 2014, (9): 38-39.)

[7] 呂高鋒, 吳慧雄. 環(huán)己醇脫氫制環(huán)己酮工藝的動(dòng)態(tài)模擬[J].北京化工大學(xué)學(xué)報(bào)(自然科學(xué)版), 2012, 39(5): 16-19. (Lü Gaofeng, WU Huixiong. Dynamic simulation of the cyclohexanol dehydrogenation process[J].Journal of Beijing University of Chemical Technology (Natural Science Edition), 2012, 39(5): 16-19.)

[8] 劉偉, 趙修仁. 結(jié)焦FCC催化劑再生過(guò)程動(dòng)力學(xué)[J].遼寧工學(xué)院學(xué)報(bào), 1993， (3): 56-61. (LIU Wei， ZHAO Xiuren. Kinetics of the regeneration progress of coked FCC catalysts[J].Journal of Liaoning Institute of Technology, 1993， (3): 56-61.)

[9] 劉偉, 李志林, 趙修仁. 催化劑在固定床中再生過(guò)程的數(shù)學(xué)模型及其求解[J].遼寧工學(xué)院學(xué)報(bào), 1994， (3): 76-79. (LIU Wei， LI Zhilin， ZHAO Xiuren. A model of the regeneration process of coked catalysts and its solution[J].Journal of Liaoning Institute of Technology, 1994， (3): 76-79.)

Simulation and Analysis of the Regeneration Process of Methanol to Propylene

LIANG Hao1, WANG Lin2, LIAO Zuwei1, YONG Xiaojing2, JIANG Binbo1,WANG Jingdai1, YANG Yongrong1, JIAO Hongqiao2

(1.StateKeyLaboratoryofChemicalEngineering,CollegeofChemicalandBiologicalEngineering,ZhejiangUniversity,Hangzhou310027,China; 2.ShenhuaNingxiaCoalGroup,Yinchuan750411,China)

An axial dynamic one-dimensional model of methanol to propylene (MTP) fixed bed regeneration reactor was established to investigate the effects of oxygen volume fraction, inlet temperature, gas flow rate on the regeneration process, as well as the distribution of temperature and carbon content in the regeneration reactor. Results showed that the factors affecting the regeneration of MTP process from large to small were oxygen volume fraction, inlet temperature and gas flow rate. A regeneration scheme was also put forward, by which the consumption of fresh nitrogen could be reduced effectively, as well as the time shortened, and a certain guiding significance would be obtained for industrial operation.

methanol to propylene (MTP); fixed bed; reactor; regeneration reactor

2016-04-13

國(guó)家自然科學(xué)基金項(xiàng)目(61590925)、“863”國(guó)際合作項(xiàng)目(2015DFA40660)、浙江省自然科學(xué)基金項(xiàng)目(LY14B060007)資助

梁昊，男，碩士研究生，從事化工過(guò)程工程方面的研究

廖祖維，男，副教授，博士，從事化工系統(tǒng)工程方面的研究；E-mail： liaozuwei@zju.edu.cn

1001-8719(2017)01-0108-07

TQ021.8

A

10.3969/j.issn.1001-8719.2017.01.015