• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    甲醇制烯烴反應(yīng)動(dòng)力學(xué)及反應(yīng)器模型研究進(jìn)展

    2014-07-05 16:05:53文堯順南海明吳秀章關(guān)豐忠孫保全
    化工進(jìn)展 2014年10期
    關(guān)鍵詞:積炭流化床丙烯

    文堯順,南海明,吳秀章,關(guān)豐忠,孫保全

    (1中國(guó)神華煤制油化工有限公司,北京 100011;2中國(guó)石油大學(xué)(北京)重質(zhì)油國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 102249)

    特約評(píng)述

    甲醇制烯烴反應(yīng)動(dòng)力學(xué)及反應(yīng)器模型研究進(jìn)展

    文堯順1,2,南海明1,吳秀章1,關(guān)豐忠1,孫保全1

    (1中國(guó)神華煤制油化工有限公司,北京 100011;2中國(guó)石油大學(xué)(北京)重質(zhì)油國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 102249)

    甲醇制烯烴(MTO)工藝是現(xiàn)代煤化工領(lǐng)域的研究熱點(diǎn),MTO反應(yīng)動(dòng)力學(xué)及其反應(yīng)器模型研究是高效反應(yīng)器開發(fā)和工業(yè)裝置操作優(yōu)化的基礎(chǔ)。本文綜述了甲醇制烯烴反應(yīng)動(dòng)力學(xué)研究進(jìn)展,詳細(xì)論述了機(jī)理型動(dòng)力學(xué)模型、八集總動(dòng)力學(xué)模型、五集總動(dòng)力學(xué)模型,指出集總動(dòng)力學(xué)模型適用于描述MTO反應(yīng)過程,如何考慮水、積炭等因素的影響是MTO動(dòng)力學(xué)研究的難點(diǎn)和關(guān)鍵。結(jié)合現(xiàn)有動(dòng)力學(xué)模型,評(píng)述了MTO反應(yīng)過程在工業(yè)規(guī)模的固定床反應(yīng)器、提升管反應(yīng)器、循環(huán)流化床反應(yīng)器、湍動(dòng)流化床反應(yīng)器中產(chǎn)物分布和轉(zhuǎn)化率模擬情況,結(jié)果表明:循環(huán)流化床反應(yīng)器和湍動(dòng)流化床反應(yīng)器適合MTO工業(yè)過程。最后指出,甲醇制烯烴反應(yīng)動(dòng)力學(xué)下一步研究方應(yīng)集中于工業(yè)規(guī)模流化床反應(yīng)器氣固兩相流動(dòng)模擬,以及與動(dòng)力學(xué)模型結(jié)合獲得準(zhǔn)確預(yù)測(cè)工業(yè)反應(yīng)結(jié)果的MTO反應(yīng)器模型。

    甲醇;烯烴;反應(yīng)動(dòng)力學(xué);反應(yīng)器;模型

    隨著我國(guó)經(jīng)濟(jì)的高速發(fā)展,對(duì)乙烯、丙烯等基本化工原料的需求持續(xù)增加。2013年我國(guó)乙烯生產(chǎn)能力達(dá)到1823萬噸/年,產(chǎn)量為1623萬噸/年,成為全球第二大乙烯生產(chǎn)國(guó)。但是我國(guó)乙烯、丙烯市場(chǎng)仍然供不應(yīng)求,目前乙烯和丙烯的對(duì)外依存度分別為45%和38%左右[1]。我國(guó)乙烯、丙烯生產(chǎn)原料主要為石腦油、輕柴油和加氫尾油等石油基原料,由于國(guó)內(nèi)油氣資源嚴(yán)重短缺,開發(fā)非石油基的低碳烯烴生產(chǎn)路線尤為重要。煤基低碳烯烴開辟了一條全新的乙烯、丙烯生產(chǎn)路線,將我國(guó)豐富的煤炭資源潔凈高效轉(zhuǎn)化為高價(jià)值的化工產(chǎn)品,在一定程度實(shí)現(xiàn)石油替代[2-5]。煤基烯烴合成路線為煤氣化、合成氣凈化、合成甲醇、甲醇制烯烴,烯烴分離、烯烴聚合。甲醇制烯烴(MTO)是連接傳統(tǒng)煤化工和石油化工之間的橋梁,世界首套甲醇制烯烴裝置于2010年8月投料成功,商業(yè)化穩(wěn)定運(yùn)行滿3年[6]。截至2013年底,在建、擬建以及處于前期規(guī)劃甲醇制烯烴裝置合計(jì)29套,總產(chǎn)能2800萬噸[7]。目前已經(jīng)具有工業(yè)應(yīng)用業(yè)績(jī)的甲醇制烯烴技術(shù)有中國(guó)科學(xué)院大連化學(xué)物理研究所的DMTO技術(shù)、中國(guó)石化上海石油化工研究院的SMTO技術(shù)以及UOP公司的UOP/hydro MTO技術(shù)。

    目前,科學(xué)界對(duì)甲醇制烯烴催化劑[8-9]、反應(yīng)特性[10]、工藝技術(shù)[11]、催化劑積炭[12-13]等方面研究總結(jié)較多,但是關(guān)于反應(yīng)動(dòng)力學(xué)研究進(jìn)展的評(píng)述鮮有報(bào)道。動(dòng)力學(xué)研究不僅是探究MTO過程反應(yīng)機(jī)理的重要手段,對(duì)高效催化劑研發(fā)的具有指導(dǎo)意義,而且也是新型反應(yīng)器和工藝開發(fā)的重要基礎(chǔ)。本文系統(tǒng)總結(jié)甲醇制烯烴過程反應(yīng)動(dòng)力學(xué)研究的最新進(jìn)展,并結(jié)合動(dòng)力學(xué)研究對(duì)不同類型MTO反應(yīng)器模擬情況進(jìn)行了介紹和評(píng)述。

    1 MTO反應(yīng)機(jī)理

    許多科學(xué)家對(duì)甲醇制烯烴反應(yīng)機(jī)理進(jìn)行了大量的研究和討論,通常認(rèn)為包括三步[10]:①甲醇在SAPO系列酸性分子篩上首先脫水生成二甲醚,并迅速達(dá)到平衡;②初始C—C鍵形成和低碳烯烴生成;③低碳烯烴進(jìn)一步生成其他高碳烴類。Haw等[14]則認(rèn)為MTO反應(yīng)分為5步,增加了反應(yīng)誘導(dǎo)期過程和催化劑積炭失活步驟。針對(duì)第二步反應(yīng)中的C—C鍵的形成,研究者先后提出了氧鎓葉立德機(jī)理、碳烯機(jī)理、碳正離子機(jī)理、自由基機(jī)理和烴池機(jī)理[10]。其中Dahl等[15-16]提出的烴池機(jī)理得到業(yè)界廣泛認(rèn)可,烴池物種具有獨(dú)特的化學(xué)性質(zhì),與普通積炭類似;MTO反應(yīng)誘導(dǎo)期就是烴池物種的形成時(shí)期,烴池物種形成后,低碳烯烴產(chǎn)物才大量生成,烯烴的生成反應(yīng)具有平行反應(yīng)特征。劉中民研究團(tuán)隊(duì)[17]研究發(fā)現(xiàn)烴池活性物種主要為多甲基苯,并在反應(yīng)中直接觀測(cè)到七甲基苯正碳離子。工業(yè)MTO裝置運(yùn)轉(zhuǎn)結(jié)果也證實(shí)烴池反應(yīng)機(jī)理的合理性[18]。

    根據(jù)烴池反應(yīng)機(jī)理,類似積炭的烴池活性物種隨著反應(yīng)進(jìn)行,活性逐漸降低,最終形成積炭。積炭一方面影響小孔SAPO-34分子篩形狀擇型性,導(dǎo)致產(chǎn)物選擇性發(fā)生變化;另一方面積炭會(huì)堵塞部分催化劑孔道,造成催化劑活性降低[19],MTO總反應(yīng)速率隨著積炭含量而變化。因此,如何考慮積炭影響對(duì)MTO動(dòng)力學(xué)研究十分重要。

    2 MTO反應(yīng)動(dòng)力學(xué)

    根據(jù)目前甲醇制烯烴動(dòng)力學(xué)的研究進(jìn)展,可以將動(dòng)力學(xué)研究分成兩大類[20]:一類是從反應(yīng)機(jī)理出發(fā),列出MTO反應(yīng)中可能出現(xiàn)的所有基元反應(yīng),通過實(shí)驗(yàn)、量子理論計(jì)算得到的機(jī)理型反應(yīng)動(dòng)力學(xué);另一類是采用集總方法,忽略基元反應(yīng)步驟,獲得的集總反應(yīng)動(dòng)力學(xué)模型。機(jī)理型反應(yīng)動(dòng)力學(xué)能夠描述反應(yīng)歷程和細(xì)節(jié),而集總動(dòng)力學(xué)模型則在反應(yīng)器設(shè)計(jì)方面實(shí)用性更強(qiáng)。表1匯總了具有代表性甲醇制烯烴反應(yīng)動(dòng)力學(xué)研究結(jié)果,下面分別進(jìn)行論述。

    表1 甲醇制烯烴反應(yīng)動(dòng)力學(xué)研究匯總

    2.1 機(jī)理型反應(yīng)動(dòng)力學(xué)

    1983年,Mihail等[21]報(bào)道了甲醇轉(zhuǎn)化為碳?xì)浠衔锏姆磻?yīng)動(dòng)力學(xué)模型,基于碳烯反應(yīng)機(jī)理,建立了包括53個(gè)基元反應(yīng)的甲醇烴類化合物的反應(yīng)動(dòng)力學(xué)方程組,將基元反應(yīng)分成12種主要的反應(yīng)步驟,獲得了動(dòng)力學(xué)參數(shù),并在固定床反應(yīng)器上進(jìn)行數(shù)學(xué)模擬。Park 和Froment[22-23]根據(jù)正碳離子機(jī)理,建立了包括氫轉(zhuǎn)移、甲基化、β消去等726個(gè)基元反應(yīng)和225種反應(yīng)中間體的動(dòng)力學(xué)模型。采用Hougen-Watson模型描述從二甲醚形成到各種產(chǎn)物生成的步驟,烯烴的生成反應(yīng)采用正碳離子機(jī)理描述。反應(yīng)網(wǎng)絡(luò)利用計(jì)算算法生成,各個(gè)基元反應(yīng)的速率常數(shù)采用單事件方法求取?;磻?yīng)的活化能通過Evans-Polanyi關(guān)聯(lián)正碳離子和烯烴異構(gòu)體的能級(jí)獲得。Forment等最終獲得33個(gè)參數(shù)的機(jī)理型動(dòng)力學(xué)模型,并對(duì)絕熱固定床、等溫固定床、以及流化床反應(yīng)器進(jìn)行模擬計(jì)算[24]。Forment的甲醇制烯烴動(dòng)力學(xué)參數(shù)基于HZSM-5分子篩催化劑實(shí)驗(yàn)結(jié)果和正碳離子機(jī)理獲得,與業(yè)界廣泛接受的烴池反應(yīng)機(jī)理不同,影響了模型的實(shí)用性。Alwahabi隨后也報(bào)道了基于SAPO-34分子篩的MTO單事件反應(yīng)動(dòng)力學(xué)模型。

    反應(yīng)機(jī)理型動(dòng)力學(xué)模型涉及反應(yīng)較多,對(duì)反應(yīng)過程的單個(gè)反應(yīng)進(jìn)行分析,對(duì)揭示反應(yīng)機(jī)理和反應(yīng)過程規(guī)律具有重要意義。但是模型預(yù)測(cè)產(chǎn)物分布的可靠性和準(zhǔn)確性還有待于進(jìn)一步證實(shí)和完善。

    2.2 集總動(dòng)力學(xué)

    集總動(dòng)力學(xué)模型是根據(jù)分子的動(dòng)力學(xué)特性,將反應(yīng)體系劃分為若干個(gè)集總組分,建立集總動(dòng)力學(xué)模型。集總動(dòng)力學(xué)模型對(duì)于處理復(fù)雜反應(yīng)體系具有廣泛應(yīng)用,例如煉油工業(yè)中的催化裂化和催化重整等采用集總動(dòng)力學(xué)模型取得了成功。與催化裂化相比,MTO反應(yīng)過程反應(yīng)原料為單一甲醇組分,相對(duì)比較簡(jiǎn)單,但是反應(yīng)產(chǎn)物包含C1~C6組分,反應(yīng)網(wǎng)絡(luò)比較復(fù)雜,也適用于采用集總方法來描述。

    2.2.1 八集總動(dòng)力學(xué)模型

    圖1 Bos甲醇制烯烴反應(yīng)八集總動(dòng)力學(xué)網(wǎng)絡(luò)[25]

    1995年,荷蘭Shell公司的Bos和Tromp在MTO動(dòng)力學(xué)研究中做出了開創(chuàng)性的工作[25],基于平行反應(yīng)的烴池機(jī)理,建立了八集總12反應(yīng)動(dòng)力學(xué)模型。動(dòng)力學(xué)模型反應(yīng)網(wǎng)絡(luò)如圖1所示,八集總包括甲醇、甲烷、乙烯、丙烯、丙烷、混合碳四、混合碳五及以上、焦炭。Bos在MTO固定床微反動(dòng)力學(xué)試驗(yàn)中未發(fā)現(xiàn)乙烷產(chǎn)物,因此動(dòng)力學(xué)模型中未包含乙烷。MTO反應(yīng)動(dòng)力學(xué)模型包括平行反應(yīng)和連續(xù)反應(yīng),其中甲醇生成甲烷、乙烯、丙烯、丙烷、混合碳四、混合碳五、焦炭的反應(yīng)為一級(jí)平行反應(yīng);丙烯具有中間反應(yīng)產(chǎn)物特性,能夠繼續(xù)反應(yīng)生成乙烯、丙烷、混合碳四以及焦炭等。其中丙烯轉(zhuǎn)化為丙烷、混合碳四以及焦炭的反應(yīng)為一級(jí)反應(yīng);而丙烯轉(zhuǎn)化生成乙烯的反應(yīng)具有二級(jí)反應(yīng)特征,對(duì)丙烯和甲醇均為一級(jí)反應(yīng),混合碳四轉(zhuǎn)化生成乙烯的反應(yīng)也具有二級(jí)反應(yīng)特征,對(duì)碳四和甲醇均為一級(jí)反應(yīng)。由于SAPO-34分子篩積炭對(duì)甲醇轉(zhuǎn)化率和產(chǎn)物選擇性產(chǎn)生重要影響,Bos在動(dòng)力學(xué)研究中采用預(yù)結(jié)焦技術(shù)進(jìn)行動(dòng)力學(xué)研究,獲得不同反應(yīng)對(duì)積炭含量的失活方程滿足冪指數(shù)型函數(shù)形式。文獻(xiàn)[25]報(bào)道了新鮮催化劑反應(yīng)速率常數(shù),甲醇轉(zhuǎn)化生成丙烯的速率常數(shù)最大,達(dá)到297.2bar/h(1bar=105Pa,下同),其次是生成乙烯和混合碳四的速率常數(shù)為149.5bar/h和103.8bar/h,甲醇生成丙烷、C5和焦炭的反應(yīng)速率常數(shù)在30~34.8bar/h。丙烯轉(zhuǎn)化的二次反應(yīng)速率常數(shù)比丙烯生成的反應(yīng)速率常數(shù)小2~3個(gè)數(shù)量級(jí),表明二次反應(yīng)程度很低。

    Bos的動(dòng)力學(xué)模型包含組分較多,研究系統(tǒng)全面,也是后續(xù)研究者采用較多的動(dòng)力學(xué)模型,具有重要的借鑒意義。但是該模型未考慮反應(yīng)溫度和水含量等重要因素對(duì)反應(yīng)結(jié)果的影響。模型僅給出了單個(gè)反應(yīng)溫度下(450℃)的速率常數(shù),未給出各個(gè)反應(yīng)的活化能,難以外推其他反應(yīng)溫度條件下的反應(yīng)速率。

    圖2 Chen甲醇制烯烴反應(yīng)八集總動(dòng)力學(xué)網(wǎng)絡(luò)[26]

    2007年,挪威科技大學(xué)的Chen研究小組等[26]在Hydro公司資助下,基于電子震蕩天平(TEOM)反應(yīng)器實(shí)驗(yàn)數(shù)據(jù)建立了MTO八集總反應(yīng)動(dòng)力模型,同樣根據(jù)烴池反應(yīng)機(jī)理,但是反應(yīng)網(wǎng)絡(luò)集總劃分與Bos的模型不同。如圖2所示,研究認(rèn)為甲烷的生成機(jī)理與其他產(chǎn)物不同,甲烷的產(chǎn)生來源于催化劑表面的甲氧基分解產(chǎn)生,單獨(dú)作為集總。認(rèn)為所有的烯烴都由二甲醚直接產(chǎn)生,甲醇迅速轉(zhuǎn)化為二甲醚且其濃度難以測(cè)定,將甲醇和二甲醚作為一個(gè)集總。甲醇和二甲醚作為含氧化合物集總的劃分方法得到很多研究者的采用。八集總包括甲烷、乙烯、丙烯、丁烯、烷烴(乙烷+丙烷)、C5、C6、和焦炭。所有的烯烴生成反應(yīng)都是二甲醚的一級(jí)平行反應(yīng),烯烴發(fā)生二次反應(yīng)生成烷烴(乙烷+丙烷)。積炭對(duì)各個(gè)反應(yīng)速率常數(shù)的影響,擬合得到失活函數(shù)與焦炭含量呈線性關(guān)系,失活速率常數(shù)隨反應(yīng)產(chǎn)物分子碳數(shù)增加而減小。模型求解獲得乙烯、丙烯、丁烯、C5、C6、乙烷+丙烷反應(yīng)動(dòng)力學(xué)速率常數(shù)、失活速率常數(shù)、指前因子和活化能等動(dòng)力學(xué)參數(shù)。Chen動(dòng)力學(xué)模型得到乙烯、丙烯和丁烯的反應(yīng)活化能分別為38.4kJ/mol、27.0kJ/mol和26.9kJ/mol?;罨芟鄬?duì)其他研究者報(bào)道結(jié)果低,Chen認(rèn)為是由于受晶內(nèi)擴(kuò)散的限制[27]。

    需要指出的是Chen的動(dòng)力學(xué)反應(yīng)網(wǎng)絡(luò)中給出了甲烷和焦炭?jī)煞N反應(yīng)產(chǎn)物,由于甲烷含量很低而且重復(fù)性較差,未能給出兩種反應(yīng)產(chǎn)物的動(dòng)力學(xué)參數(shù)。研究者對(duì)Chen實(shí)驗(yàn)數(shù)據(jù)和參數(shù)進(jìn)行重新回歸擬合后,發(fā)現(xiàn)數(shù)據(jù)的相關(guān)性差,模型的實(shí)用性還有待于進(jìn)一步驗(yàn)證。

    2.2.2 五集總動(dòng)力學(xué)

    八集總動(dòng)力學(xué)研究結(jié)果表明,二次反應(yīng)速率遠(yuǎn)小于一次反應(yīng),研究者對(duì)模型進(jìn)行適當(dāng)簡(jiǎn)化,建立了五集總動(dòng)力學(xué)模型。Gayubo研究小組在甲醇轉(zhuǎn)化制烴類化合物動(dòng)力學(xué)方面研究較多,研究了甲醇在SAPO-34[28]、SAPO-18[29]、HZSM-5[30]分子篩催化劑上的反應(yīng)動(dòng)力學(xué)。Gayubo等[28]根據(jù)MTO反應(yīng)的積炭生成規(guī)律,認(rèn)為在反應(yīng)初始階段積炭反應(yīng)速率較大,隨后積炭含量隨運(yùn)轉(zhuǎn)時(shí)間變化較慢。將動(dòng)力學(xué)研究區(qū)間確定為積炭含量變化較小的平穩(wěn)階段,在該反應(yīng)區(qū)間,忽略積炭含量變化對(duì)反應(yīng)速率的影響,建立如圖3所示動(dòng)力學(xué)模型。Gayubo的動(dòng)力學(xué)模型忽略了慢反應(yīng)步驟,同時(shí)考慮了水含量和反應(yīng)溫度對(duì)反應(yīng)的影響,定量表征水含量的影響程度,其動(dòng)力學(xué)模型在MTO反應(yīng)器設(shè)計(jì)中更具實(shí)用性,主要反映乙烯、丙烯和丁烯的反應(yīng)動(dòng)力學(xué)參數(shù)。研究認(rèn)為反應(yīng)體系中的水與反應(yīng)物、反應(yīng)產(chǎn)物競(jìng)爭(zhēng)催化劑表面的活性中心,能夠抑制積炭的生成,有利于低碳烯烴的生成。Gayubo等采用外推法獲得初始反應(yīng)時(shí)刻的產(chǎn)物分布,不是實(shí)驗(yàn)測(cè)定值,不同的外推曲線會(huì)引起較大的誤差,同時(shí)由于忽略了積炭的影響,在預(yù)測(cè)實(shí)際反應(yīng)產(chǎn)物分布時(shí)會(huì)出現(xiàn)一定的偏差。

    圖3 Gayubo MTO反應(yīng)動(dòng)力學(xué)模型反應(yīng)網(wǎng)絡(luò)[28]

    國(guó)內(nèi)研究者也在甲醇制烯烴反應(yīng)動(dòng)力學(xué)方面開展了大量的研究,取得了不錯(cuò)的進(jìn)展。中國(guó)石化的齊國(guó)禎對(duì)MTO反應(yīng)動(dòng)力學(xué)進(jìn)行了系統(tǒng)和全面的研究,提供了豐富的動(dòng)力學(xué)反應(yīng)數(shù)據(jù),主要研究結(jié)果見文獻(xiàn)[31-32]。齊國(guó)禎[32]對(duì)五集總動(dòng)力學(xué)模型進(jìn)行了修正,將丁烯集總變?yōu)镃4集總,將C5+集總調(diào)整為烷烴集總(包括乙烷、丙烷和C5+),如圖4所示。同時(shí)選擇反應(yīng)為1.5min的反應(yīng)結(jié)果進(jìn)行動(dòng)力學(xué)數(shù)據(jù)擬合,更能準(zhǔn)確的反映初始階段反應(yīng)特性。模型認(rèn)為反應(yīng)過程中生成的水對(duì)主反應(yīng)和副反應(yīng)均有貢獻(xiàn),在反應(yīng)中引入水的吸附阻力項(xiàng),定量分析水對(duì)MTO各個(gè)反應(yīng)的影響程度。

    圖4 齊國(guó)禎MTO五集總動(dòng)力學(xué)模型反應(yīng)網(wǎng)絡(luò)[32]

    圖5 齊國(guó)禎MTO五集總動(dòng)力學(xué)模型速率方程[32]

    通過模型求解,齊國(guó)禎最終得到如圖5所示動(dòng)力學(xué)方程,模型類似于雙曲線型動(dòng)力學(xué)方程,分子為速率常數(shù)項(xiàng)與推動(dòng)力項(xiàng)的乘積,分母為水的吸附阻力項(xiàng)與積炭引起的失活阻力項(xiàng)的乘積。反應(yīng)模型能夠預(yù)測(cè)反應(yīng)溫度、水醇比和重時(shí)空速對(duì)反應(yīng)速率的影響。由圖5可知,C1集總的反應(yīng)活化能最高,達(dá)到100.56kJ/mol,隨著產(chǎn)物碳數(shù)的增加,反應(yīng)活化能逐漸降低。乙烯和丙烯反應(yīng)活化能分別為89.05kJ/mol和77.58kJ/mol,乙烯活化能高于丙烯,表明乙烯反應(yīng)速率隨著反應(yīng)溫度升高增加速率高于丙烯。工業(yè)MTO裝置運(yùn)轉(zhuǎn)過程中,提高反應(yīng)溫度,乙烯/丙烯比增大,也證實(shí)乙烯和丙烯反應(yīng)活化能的差異。通過對(duì)不同反應(yīng)溫度條件下失活速率常數(shù)進(jìn)行擬合,得到失活級(jí)數(shù)為2級(jí)。該動(dòng)力學(xué)模型能夠較好的預(yù)測(cè)各種產(chǎn)物轉(zhuǎn)化率隨著甲醇空時(shí)的變化規(guī)律。

    齊國(guó)禎的MTO動(dòng)力學(xué)模型考慮了水的影響、反應(yīng)溫度的影響以及催化劑積炭失活的影響。

    華東理工大學(xué)的胡浩等[33-34]也建立了MTO五集總反應(yīng)動(dòng)力學(xué)模型,與齊國(guó)禎模型不同之處在于將甲烷單獨(dú)作為集總,忽略了CO和CO2。在動(dòng)力學(xué)數(shù)據(jù)選擇上,胡浩與其他研究者不同,認(rèn)為甲醇在SAPO-34分子篩上進(jìn)行MTO反應(yīng)基本沒有反應(yīng)穩(wěn)定期,即使選擇最初使反應(yīng)為1~2min的反應(yīng)數(shù)據(jù),隨著重時(shí)空速變化,反應(yīng)誘導(dǎo)期和積炭失活對(duì)動(dòng)力學(xué)的影響不可消除,難以獲得本征動(dòng)力學(xué)數(shù)據(jù)。在動(dòng)力學(xué)實(shí)驗(yàn)中,選擇乙烯加丙烯總收率最高時(shí)的產(chǎn)物分布作為動(dòng)力學(xué)研究的基本實(shí)驗(yàn)數(shù)據(jù)。根據(jù)平行反應(yīng)機(jī)理獲得冪級(jí)數(shù)型動(dòng)力學(xué)模型,在乙烯丙烯收率最高的時(shí),甲醇轉(zhuǎn)化總的反應(yīng)活化能為52.4kJ/mol,乙烯、丙烯和丁烯反應(yīng)活化能分別為52.7kJ/mol、50.9kJ/mol和47.5kJ/mol,各反應(yīng)活化能比較相近,比齊國(guó)禎所得活化能小,但高于Chen的動(dòng)力學(xué)模型。詳細(xì)動(dòng)力學(xué)參數(shù)如表2所示。該動(dòng)力模型對(duì)反應(yīng)器出口產(chǎn)物分布預(yù)測(cè)效果較好,對(duì)乙烯、丙烯和丁烯平均誤差小于5%,最大相對(duì)誤差絕對(duì)值小于11%。

    表2 胡浩五集總動(dòng)力學(xué)模型參數(shù)[34]

    胡浩動(dòng)力學(xué)模型的數(shù)據(jù)基礎(chǔ)為乙烯加丙烯選擇性最高的時(shí)間段,工業(yè)MTO裝置操作基本在上述階段,因而采用胡浩的五集總動(dòng)力學(xué)模型應(yīng)用于工業(yè)MTO裝置的產(chǎn)物分布預(yù)測(cè),相對(duì)誤差可能較小。

    3 MTO動(dòng)力學(xué)在反應(yīng)器模擬方面的研究

    在MTO反應(yīng)動(dòng)力學(xué)研究基礎(chǔ)上,結(jié)合反應(yīng)器內(nèi)催化劑流動(dòng)特性、氣固傳質(zhì)、傳熱等進(jìn)行MTO反應(yīng)器模擬研究。一方面給出不同種類反應(yīng)器MTO反應(yīng)轉(zhuǎn)化率和產(chǎn)物分布,優(yōu)選適合于甲醇轉(zhuǎn)化生產(chǎn)低碳烯烴工業(yè)過程的理想反應(yīng)器;另一方面,也為工業(yè)MTO反應(yīng)器的設(shè)計(jì)和MTO裝置的優(yōu)化操作提供基礎(chǔ)數(shù)據(jù)。業(yè)界在MTO反應(yīng)器模型研究方面開展了大量的工作。

    3.1 固定床反應(yīng)器

    Bos[25]和Alwahabi等[35]模擬研究表明固定床不宜作為MTO反應(yīng)器,主要有以下4個(gè)方面的缺陷:①M(fèi)TO反應(yīng)放熱量大,固定床反應(yīng)器取熱困難;②催化劑失活較快,需要多個(gè)反應(yīng)器進(jìn)行切換;③固定床反應(yīng)器再生困難,容易出現(xiàn)飛溫;④反應(yīng)器壓降大。此外,工業(yè)規(guī)模固定床反應(yīng)器若采用等溫操作,采用列管式反應(yīng)器造價(jià)高,裝置投資大;當(dāng)采用多層絕熱反應(yīng)器操作時(shí),設(shè)置層間取熱和多反應(yīng)器切換,工藝較為復(fù)雜。

    3.2 流化床反應(yīng)器

    流化床反應(yīng)器能夠?qū)崿F(xiàn)催化劑的連續(xù)反應(yīng)和再生,反應(yīng)器和再生器取熱方便。目前,工業(yè)化的MTO技術(shù)均采用流化床工藝。研究者對(duì)不同流化床反應(yīng)器如提升管反應(yīng)器、循環(huán)流化床反應(yīng)器、湍動(dòng)流化床反應(yīng)器進(jìn)行了MTO反應(yīng)模擬。

    3.2.1 提升管反應(yīng)器

    Bos等[25]基于前述8集總12反應(yīng)動(dòng)力學(xué)模型,以30萬噸乙烯/年的工業(yè)規(guī)模MTO裝置為例,模擬計(jì)算了提升管反應(yīng)器甲醇轉(zhuǎn)化率、產(chǎn)物分布和反應(yīng)器尺寸。研究發(fā)現(xiàn),催化劑采用不完全再生方式,要達(dá)到理想的低碳烯烴選擇性,提升管長(zhǎng)度高達(dá)150m。齊國(guó)禎等[32]采用五集總動(dòng)力學(xué)模型,模擬了直徑0.85m、高度30m、操作線速8m/s、催化劑循環(huán)量60t/h工況下的反應(yīng)結(jié)果,提升管反應(yīng)器出口轉(zhuǎn)化率僅為53.4%,溫升達(dá)到40℃;當(dāng)增加催化劑循環(huán)量到400t/h時(shí),甲醇轉(zhuǎn)化率達(dá)到99%,溫升降低到10℃左右。但是在上述工況下,裝置處理量?jī)H為6.8萬噸。與催化裂化相比,MTO裝置進(jìn)料甲醇分子遠(yuǎn)小于催化原料分子,如果MTO裝置達(dá)到百萬噸甲醇處理規(guī)模,提升管反應(yīng)器將十分巨大。此外催化劑單程失活小,再生不方便。因此,提升管反應(yīng)器不適宜作為工業(yè)規(guī)模MTO反應(yīng)器。

    3.2.2 循環(huán)流化床反應(yīng)器

    循環(huán)流化床反應(yīng)器包括密相提升管反應(yīng)器部分和用于換熱的循環(huán)部分,Bos 模擬結(jié)果表明[25],當(dāng)催化劑停留時(shí)間為7min時(shí),催化劑平均積炭含量為8.2%(質(zhì)量分?jǐn)?shù)),乙烯丙烯比例接近1。當(dāng)停留時(shí)間進(jìn)一步增大時(shí),乙烯/丙烯比增大,反應(yīng)器體積也明顯增大。Soundararajan等[36]也研究了MTO在循環(huán)流化床反應(yīng)器的產(chǎn)物分布,由于催化劑停留時(shí)間短,單程積炭較少,經(jīng)旋風(fēng)分離的催化劑大部分返回提升管底部,僅有少部分再生,從而提高提升管內(nèi)催化劑的平均積炭含量。研究發(fā)現(xiàn)提升管出口結(jié)構(gòu)形式和積炭含量對(duì)產(chǎn)物分布有較大的影響,當(dāng)積炭含量為5%時(shí),反應(yīng)產(chǎn)物低碳烯烴選擇性最高。該研究結(jié)果對(duì)于循環(huán)流化床反應(yīng)器設(shè)計(jì)和操作條件選取具有重要指導(dǎo)意義。Schoenfelder等[37]也在MTO循環(huán)流化床反應(yīng)器模擬開展了研究工作。

    3.2.3 湍動(dòng)流化床

    工業(yè)MTO裝置采用湍動(dòng)流化床反應(yīng)器[38],對(duì)湍動(dòng)流化床反應(yīng)器模型研究,有利于反應(yīng)器設(shè)計(jì)和工業(yè)裝置操作優(yōu)化。浙江大學(xué)鄭康等[39-40]基于Bos的反應(yīng)動(dòng)力學(xué)數(shù)據(jù)對(duì)MTO湍動(dòng)流化床反應(yīng)器進(jìn)行了模擬,采用流化床兩相模型描述流化床MTO反應(yīng)過程,利用拉格朗日顆粒跟蹤法模擬顆粒的積炭和反應(yīng)過程。研究發(fā)現(xiàn),MTO本征反應(yīng)速率遠(yuǎn)大于氣泡相和乳化相的相間傳質(zhì)速率,湍動(dòng)流化床的相間傳質(zhì)是速率控制步驟;強(qiáng)化相間傳質(zhì)和延長(zhǎng)催化劑在反應(yīng)器中的停留時(shí)間均能夠顯著提高催化劑的積炭量,有利于提高乙烯和丙烯的選擇性。催化劑的積炭含量是影響低碳烯烴選擇性的關(guān)鍵性因素,減少催化劑的返混和采用多級(jí)串聯(lián)操作對(duì)于提高低碳烯烴選擇性影響不大,采用多級(jí)氣固逆流接觸反而降低甲醇轉(zhuǎn)化率。需要指出的是,由于缺少工業(yè)裝置和萬噸級(jí)中試裝置數(shù)據(jù),鄭康的流化床模型僅對(duì)比了模擬結(jié)果與中科院大連化物所專利報(bào)道中試結(jié)果[41]。

    Bos等[25]湍動(dòng)流化床模擬結(jié)果證實(shí),對(duì)于30萬噸/年乙烯生產(chǎn)規(guī)模,湍動(dòng)流化床反應(yīng)器直徑為8.5m即可滿足要求,反應(yīng)氣速為1m/s,低于循環(huán)流化床的3m/s。隨著停留時(shí)間延長(zhǎng),丙烯選擇性在初始階段略有下降,然后基本保持不變,而乙烯選擇性則逐漸增大。

    鼓泡流化床操作氣速低于湍動(dòng)流化床,學(xué)者也對(duì)模擬了MTO鼓泡流化床反應(yīng)性能,比較有代表性的是Alwahabi等[35]和齊國(guó)禎[32]。前者模擬了年產(chǎn)55000噸乙烯和丙烯的生產(chǎn)規(guī)模的MTO裝置流化床反應(yīng)結(jié)果,進(jìn)料溫度為430℃,出口溫升21℃,甲醇轉(zhuǎn)化率為90%,乙烯丙烯收率相當(dāng)于進(jìn)料甲醇的12%和13%(質(zhì)量分?jǐn)?shù))。需要指出的是,前述所有的MTO動(dòng)力學(xué)基礎(chǔ)數(shù)據(jù)均是基于固定床或者微反試驗(yàn)取得的,Kaarsholm等[42]報(bào)道了基于流化床反應(yīng)器試驗(yàn)數(shù)據(jù)的動(dòng)力學(xué)MTO動(dòng)力學(xué)研究結(jié)果,對(duì)于指導(dǎo)工業(yè)MTO反應(yīng)器的設(shè)計(jì)具有意義。但是研究采用的催化劑為磷改性的ZSM-5催化劑,影響了其實(shí)用性。

    4 結(jié)語和展望

    自MTO反應(yīng)發(fā)現(xiàn)以來,業(yè)界對(duì)MTO反應(yīng)機(jī)理和反應(yīng)動(dòng)力學(xué)研究一直持續(xù)MTO整個(gè)研究過程。由于對(duì)反應(yīng)機(jī)理一直爭(zhēng)論不斷,當(dāng)烴池機(jī)理被廣泛接受后,MTO反應(yīng)動(dòng)力學(xué)研究取得了很大的進(jìn)展。集總動(dòng)力學(xué)模型對(duì)描述MTO反應(yīng)性能和產(chǎn)物分布擬合效果較好,為甲醇制烯烴過程工業(yè)反應(yīng)器的設(shè)計(jì)和工業(yè)化進(jìn)程提供了強(qiáng)有力的支持。研究者對(duì)MTO流化床內(nèi)流動(dòng)規(guī)律[43]、數(shù)值模擬[44-45]進(jìn)行了研究,取得不錯(cuò)的研究結(jié)果,但是目前的動(dòng)力學(xué)模型對(duì)工業(yè)MTO流化床反應(yīng)器的模擬研究結(jié)果鮮有報(bào)道,對(duì)于工業(yè)MTO流化床反應(yīng)器,催化劑的流動(dòng)規(guī)律、氣固傳質(zhì)、停留時(shí)間分布對(duì)MTO反應(yīng)產(chǎn)物分布具有重要影響。對(duì)MTO工業(yè)反應(yīng)器兩相流進(jìn)行模擬,準(zhǔn)確描述反應(yīng)器中的氣固流場(chǎng)分布,進(jìn)而結(jié)合MTO反應(yīng)動(dòng)力學(xué),建立流化床反應(yīng)器模型,將是MTO動(dòng)力學(xué)下一步的研究重點(diǎn)。此外,由于MTO過程積炭對(duì)催化劑活性和產(chǎn)物形狀擇形性影響顯著。在動(dòng)力學(xué)模型中不僅考慮積炭含量的影響,同時(shí)考慮操作條件對(duì)積炭生成速率的影響,對(duì)建立更為合理的MTO動(dòng)力學(xué)模型具有重要意義,更好的指導(dǎo)工業(yè)MTO反應(yīng)器的設(shè)計(jì)和裝置的操作優(yōu)化。

    符 號(hào) 說 明

    Ea,i——反應(yīng)活化能,kJ/mol

    FMO——甲醇進(jìn)料速率,g/h

    Kd——失活速率常數(shù)

    Kw——水吸附阻力因子

    k1——甲醇轉(zhuǎn)化為C1的速率常數(shù)

    k2——甲醇轉(zhuǎn)化為乙烯的速率常數(shù)

    k3——甲醇轉(zhuǎn)化為丙烯的速率常數(shù)

    k4——甲醇轉(zhuǎn)化為碳四的速率常數(shù)

    k5——甲醇轉(zhuǎn)化為烷烴的速率常數(shù)

    k0——指前因子,mol/(g·kPa·h)

    it——反應(yīng)時(shí)間,s

    W——催化劑裝填量,g

    X1——C1的摩爾分率

    XA——甲醇轉(zhuǎn)化摩爾分率

    Xb——碳四的摩爾分率

    Xe——乙烯的摩爾分率

    Xp——丙烯的摩爾分率

    Xr——烷烴的摩爾分率

    Xw——反應(yīng)過程中水的摩爾分率

    Xwo——甲醇進(jìn)料中水的摩爾分率

    [1]吳秀章. 煤制低碳烯烴工藝與工程[M]. 北京:化學(xué)工業(yè)出版社. 2014:2-5.

    [2]陳香生,劉昱,陳俊武. 煤基甲醇制烯烴(MTO)工藝生產(chǎn)低碳烯烴的工程技術(shù)及投資分析[J]. 煤化工,2005,5:6-11.

    [3]李宏圖. 煤制低碳烯烴的技術(shù)路線及現(xiàn)狀分析[J]. 中國(guó)煤炭,2006,32(10):47-49.

    [4]張玉卓. 煤潔凈轉(zhuǎn)化工程(神華煤制燃料和合成材料技術(shù)探索與工程實(shí)踐)[M]. 北京:煤炭工業(yè)出版社, 2011.

    [5]項(xiàng)東,彭麗娟,楊思宇,等. 石油與煤路線制烯烴過程技術(shù)評(píng)述[J]. 化工進(jìn)展,2013,32(5):959-920.

    [6]吳秀章. 煤制低碳烯烴工業(yè)示范工程最新進(jìn)展[J]. 化工進(jìn)展,2014,33(4):787-794.

    [7]杜善明. 現(xiàn)代煤化工產(chǎn)業(yè)發(fā)展現(xiàn)狀及發(fā)展建議[J]. 科技和產(chǎn)業(yè),2012,12(12):41-44.

    [8]黃曉昌,方奕文,喬曉輝. 甲醇制烯烴催化劑及其反應(yīng)機(jī)理研究進(jìn)展[J]. 工業(yè)催化,2008,16(1):23-27.

    [9]吳磊,劉子玉,夏林,等. SAPO-34晶粒形貌對(duì)甲醇轉(zhuǎn)化制低碳烯烴反應(yīng)的影響[J]. 催化學(xué)報(bào),2013,34(7):1348-1356.

    [10]St?ecker M. Methanol-to-hydrocarbons:Catalytic materials and their behavior[J].Micro. Meso. Mater.,1999,29:3-48.

    [11]朱杰,崔宇,陳元君,等. 甲醇制烯烴過程研究進(jìn)展[J]. 化工學(xué)報(bào),2010,61(7):1674-1684.

    [12]邢愛華,朱偉平,岳國(guó),等. 甲醇制烯烴反應(yīng)催化劑積炭問題研究進(jìn)展[J]. 化工進(jìn)展,2011,30(8):1717-1725.

    [13]南海明,文堯順,吳秀章,等. 甲醇制烯烴工業(yè)裝置催化劑積炭特性模擬研究[J]. 石油化工,2014,43(1):79-86.

    [14]Haw J F,Song W,Marcus D M. The mechanism of methanol to hydrocarbon catalysis[J].Acc. Chem. Res.,2003,36:317-326.

    [15]Dahl I M,Kolboe S. On the reaction-mechanism forhydrocarbon formation in the MTO reaction over SAPO-34[J].Catal. Lett.,1993,(20):329-336.

    [16]Dahl I M,Kolboe S. On the reaction mechanism for hydrocarbon formation from mechanism studies[J].J. Catal.,1994,149:458-464.

    [17]Li J Z,Wei Y X,Chen J R,et al. Observation of heptamethylbenzeniumcation over SAPO-Type molecular sieve DNL-6 under real MTO conversion conditions[J].J. Am. Chem. Soc.,2012,134:836-839.

    [18]吳秀章. 基于工業(yè)裝置探究甲醇制烯烴反應(yīng)機(jī)理[J]. 石化技術(shù)與應(yīng)用,2013,31(5):363-369.

    [19]胡浩,葉麗萍,應(yīng)衛(wèi)勇,等. 甲醇制烯烴反應(yīng)機(jī)理和動(dòng)力學(xué)研究進(jìn)展[J]. 工業(yè)催化,2008,16(3):18-23.

    [20]Keil F J. Methanol-to-hydrocarbon:Process technology[J].Microporous Mesoporous Mater.,1999,29:49-66.

    [21]Mihail R,Straja S,Maria GH,et al. A kinetic model for methanol conversion to hydrocarbons[J].Chem. Eng. Sci.,1983,38(9):1581-1591.

    [22]Park T Y,F(xiàn)roment G F. Kinetic model of the methanol to olefins process. 1. Model formulation[J].Ind. Eng. Chem. Res.,2001,40:4172-4186.

    [23]Park T Y,F(xiàn)roment G F. Kinetic model of the methanol to olefins process. 2. Experimental results,model discrimination,and parameter estimation[J].Ind. Eng. Chem. Res.,2001,40:4187-4196.

    [24]Park T Y,F(xiàn)roment G F. Analysis of fundamental reaction rates in the methanol-to-olefins process on ZSM-5 as a basis for reactor design and operation[J].Ind. Eng. Chem. Res.,2004,43:682-689.

    [25]Bos A N,Tromp P J,Akse H N. Conversion of methanol to lower olefins kinetic modeling,reactor simulation,and selection[J].Ind. Eng. Chem. Res.,1995,34:3808-3816.

    [26]Chen D,Gronvold A,Moljord K,et al. Methanol conversion to light olefins over SAPO-34:Reaction network and deactivation kinetics[J].Ind. Eng. Chem. Res.,2007,46:4116-4123.

    [27]Chen D,Moljord K,F(xiàn)uglerud T,et al. The effect of crystal size of SAPO-34 on the selectivity and deactivation of the MTO reaction[J].Microporous Mesoporous Mater.,1999,29(1):191-203.

    [28]Gayubo A G,Aguayu A T,Sanchez del Campo A E,et al,Kinetic modeling of methanol transformation into olefins on a SAPO-34 catalyst[J].Ind. Eng. Chem. Res.,2000,39:292-300.

    [29]Gayubo A G,Aguayu A T,Alonso A,et al. Kinetic modeling of the methanol to olefins process on a silicoaluminophosphate (SAPO-18)catalyst by considering deactivation and the formation of individual olefins[J].Ind. Eng. Chem. Res.,2007,46:1981-1989.

    [30]Gayubo A G,Aguayu A T,Castilla M,et al. Catalyst reactivation kinetics for methanol transformation into hydrocarbons. Expressions for designing reaction-regeneration cycles in isothermal and adiabatic fixed bed reactor[J].Chem. Eng. Sci.,2001,56:5059-5071.

    [31]齊國(guó)禎,馬濤,劉紅星,等. 甲醇制烯烴反應(yīng)動(dòng)力學(xué)[J]. 化工學(xué)報(bào),2005,56(12):2326-2331

    [32]齊國(guó)禎. 甲醇制烯烴(MTO)反應(yīng)過程研究[D]. 上海:華東理工大學(xué),2006:31-51.

    [33]胡浩,葉麗萍,應(yīng)衛(wèi)勇,等. SAPO-34分子篩催化劑上甲醇制烯烴反應(yīng)的本征動(dòng)力學(xué)[J].華東理工大學(xué)學(xué)報(bào):自然科學(xué)版,2009,35(5):655-660.

    [34]胡浩. 甲醇制烯烴(MTO)催化反應(yīng)過程的研究[D]. 上海:華東理工大學(xué),2010:38-51.

    [35]Alwahabi S M,F(xiàn)orment G F. Conceptual reactor design for the methanol -to- olefins process on SAPO-34[J].Ind. Eng. Chem. Res.,2004,43:5112-5122.

    [36]Soundararajan S,Dalai A K,Berruti F. Modeling of methanol to olefin (MTO) process in a circulating fluidized bed reactor[J].Fuel,2001,80:1187-1197.

    [37]Schoenfelder H,Hinderer J,Werther J,et al. Methanol to olefins prediction of the performance of a circulating fluidized-bed reactor on the basis of kinetic experiments in a fixed bed reactor[J].Chem. Eng.Sci.,1994,49(24B):5377-5390.

    [38]吳秀章. 流態(tài)化技術(shù)在甲醇制烯烴(MTO)工業(yè)裝置中的應(yīng)用[C]//第七屆全國(guó)流態(tài)化會(huì)議,北京,2013.

    [39]鄭康,成有為,李希. 甲醇制烯烴流化床反應(yīng)器的模擬與分析[J].高校化學(xué)工程學(xué)報(bào),2012,21(1):69-76.

    [40]鄭康,成有為,李希. 甲醇制烯烴流化床反應(yīng)器的模擬與分析[D].杭州:浙江大學(xué),2011.

    [41]Qi Y,Liu Z M. Lv Z H,et al. Methanol for producing light olefins from methanol or/and dimethyl ether:US,20100063336A1[P]. 2010-03-11.

    [42]Kaarsholm M,Refii B,Joensen F,et al. Kinetic modeling of methanol-to-olefin reaction over ZSM-5 in fluid bed[J].Ind. Eng. Chem. Res.,2010,49(1):29-38.

    [43]Wang B,Li T,Sun Q W,et al. Experimental study on flow behavior in a gas-solid fluidized bed for the methanol-to-olefins process[J].Chemical Engineering & Technology,2010,33(10):1591-1600.

    [44]Chen X M,Luo Z H,Zhu Y P,et al. Direct concurrent multi-scale CFD modeling:The effect of intraparticle transfer on the flow field in a MTO FBR[J].Chem. Eng. Sci.,2013,104:690-700.

    [45]Zhuang Y Q,Chen X M,Luo Z H,et al. CFD-DEM modeling of gas-solid flow and catalytic MTO reaction in a fluidized bed reactor[J].Computers & Chemical Engineering,2014,60:1-16.

    Advances in reaction kinetics and reactor models for methanol to olefins reaction

    WEN Yaoshun1,2,NAN Haiming1,WU Xiuzhang1,GUAN Fengzhong1,SUN Baoquan1
    (1China Shenhua Coal to Liquid and Chemical Co.,Ltd.,Beijing 100011,China;2State Key Laboratory of Heavy Oil Processing,China University of Petroleum,Beijing 102249,China)

    The methanol to olefins (MTO) process has attracted considerable attention in modern coal chemical industry. MTO reaction kinetics and reactor models play a key role in designing highly effective commercial reactor and optimizing operation conditions of MTO units. In this paper the advance in MTO reaction kinetics research is reviewed. Mechanistic kinetic model,eight lumps model,and five lumps model are discussed. The lumping models are suitable for describing MTO process,while coking and water content have a significant influence on reaction behavior. According to kinetics models,simulations of product distribution and methanol conversion in commercial scale MTO process in fixed bed reactor,riser reactor,circulating fluidized reactor and turbulent fluidized reactor are summarized. The circulating fluidized reactor and turbulent fluidized reactor are ideal reactors for commercial MTO process. Further researches on kinetics of MTO are also proposed.

    methanol;olefin;reaction kinetics;reactor;model

    TQ 013.2

    A

    1000-6613(2014)10-2521-08

    10.3969/j.issn.1000-6613.2014.10.001

    2014-04-22;修改稿日期:2014-06-04。

    神華集團(tuán)科技創(chuàng)新項(xiàng)目(SHJT-13-24)。

    文堯順(1982—),男,博士,工程師,從事甲醇制烯烴領(lǐng)域技術(shù)研究。E-mail wenyaoshun@163.com。

    猜你喜歡
    積炭流化床丙烯
    錢愛康
    《共生》主題系列作品
    大眾文藝(2022年16期)2022-09-07 03:08:04
    苯丙烯菌酮
    流化床丙烷脫氫反應(yīng)段的模擬及優(yōu)化
    淺談發(fā)動(dòng)機(jī)積炭
    關(guān)于循環(huán)流化床鍋爐集控運(yùn)行研究
    液化氣中的丙烯含有的雜質(zhì)對(duì)丙烯聚合反應(yīng)的影響
    重整催化劑Ni2P/Al2O3-SAPO-11積炭失活規(guī)律研究
    HZSM-5催化劑氣相沉積改性及催化甲苯歧化抗積炭性能
    單沉浸管流化床內(nèi)離散顆粒數(shù)值模擬
    免费人成在线观看视频色| 国产一区二区在线观看av| 一级片'在线观看视频| 男男h啪啪无遮挡| 亚洲av综合色区一区| a级一级毛片免费在线观看| 99视频精品全部免费 在线| 久久久久精品久久久久真实原创| 91久久精品国产一区二区成人| 一级av片app| 国产精品人妻久久久久久| 大话2 男鬼变身卡| 国产亚洲欧美精品永久| 下体分泌物呈黄色| a级一级毛片免费在线观看| 永久免费av网站大全| 女人精品久久久久毛片| 一本一本综合久久| 午夜视频国产福利| 这个男人来自地球电影免费观看 | 六月丁香七月| 三上悠亚av全集在线观看 | 插阴视频在线观看视频| 国产色爽女视频免费观看| 免费av不卡在线播放| av播播在线观看一区| 日韩强制内射视频| 亚洲成人一二三区av| 久久亚洲国产成人精品v| 亚洲精品亚洲一区二区| 久久ye,这里只有精品| 久久精品夜色国产| 欧美人与善性xxx| 免费高清在线观看视频在线观看| 我要看日韩黄色一级片| 少妇 在线观看| 成人午夜精彩视频在线观看| 国产精品成人在线| 中文字幕av电影在线播放| 嘟嘟电影网在线观看| 欧美激情国产日韩精品一区| 又大又黄又爽视频免费| 午夜91福利影院| 亚洲不卡免费看| 久久久久久久久久成人| 五月伊人婷婷丁香| 精品卡一卡二卡四卡免费| 亚洲图色成人| 亚洲欧美精品自产自拍| 欧美少妇被猛烈插入视频| 国产 精品1| 免费人成在线观看视频色| 国产成人精品无人区| 777米奇影视久久| 一本一本综合久久| 久久精品熟女亚洲av麻豆精品| 纵有疾风起免费观看全集完整版| 精品久久久久久久久av| 久久久久精品性色| 欧美成人精品欧美一级黄| 一个人看视频在线观看www免费| 中文欧美无线码| 亚洲情色 制服丝袜| 建设人人有责人人尽责人人享有的| 欧美高清成人免费视频www| 汤姆久久久久久久影院中文字幕| 久久久久久久精品精品| 高清在线视频一区二区三区| av网站免费在线观看视频| 2018国产大陆天天弄谢| 日韩一区二区三区影片| 久久亚洲国产成人精品v| 久久99蜜桃精品久久| 观看美女的网站| 亚洲精品中文字幕在线视频 | 色哟哟·www| 国产 精品1| 久久毛片免费看一区二区三区| 午夜日本视频在线| av国产精品久久久久影院| 中文字幕制服av| 免费大片黄手机在线观看| 日韩中字成人| 国产亚洲精品久久久com| 午夜影院在线不卡| 久久久午夜欧美精品| 99九九线精品视频在线观看视频| 欧美 日韩 精品 国产| 女人久久www免费人成看片| 搡女人真爽免费视频火全软件| av国产久精品久网站免费入址| 亚洲国产精品国产精品| 欧美精品人与动牲交sv欧美| 亚洲av电影在线观看一区二区三区| 欧美丝袜亚洲另类| 大陆偷拍与自拍| 亚洲一级一片aⅴ在线观看| 亚洲av免费高清在线观看| 日本欧美视频一区| 国产精品福利在线免费观看| 校园人妻丝袜中文字幕| 极品教师在线视频| 日本-黄色视频高清免费观看| 自拍欧美九色日韩亚洲蝌蚪91 | 69精品国产乱码久久久| 精品国产国语对白av| 国产日韩欧美亚洲二区| freevideosex欧美| 校园人妻丝袜中文字幕| 在线观看一区二区三区激情| 免费av不卡在线播放| 免费观看a级毛片全部| 国产高清国产精品国产三级| 亚洲精品国产av成人精品| 国产色婷婷99| 国产一区二区三区av在线| 男男h啪啪无遮挡| 一个人免费看片子| 欧美亚洲 丝袜 人妻 在线| 亚洲情色 制服丝袜| 三上悠亚av全集在线观看 | 少妇熟女欧美另类| 国产白丝娇喘喷水9色精品| 国产成人免费观看mmmm| 亚洲精品成人av观看孕妇| 亚洲无线观看免费| 久久久久久久久大av| 国产片特级美女逼逼视频| 国产又色又爽无遮挡免| 草草在线视频免费看| 春色校园在线视频观看| 亚洲,欧美,日韩| 亚洲美女视频黄频| 狂野欧美激情性bbbbbb| 亚洲成色77777| 亚洲精品乱久久久久久| 啦啦啦啦在线视频资源| 制服丝袜香蕉在线| 国产精品女同一区二区软件| 丝袜喷水一区| 亚洲av欧美aⅴ国产| 青春草亚洲视频在线观看| 欧美日韩亚洲高清精品| 天天操日日干夜夜撸| 成人影院久久| 免费黄网站久久成人精品| 亚洲电影在线观看av| 国产亚洲午夜精品一区二区久久| 纯流量卡能插随身wifi吗| 色网站视频免费| 一级a做视频免费观看| 色视频www国产| 伦理电影大哥的女人| 国产精品人妻久久久影院| 99热6这里只有精品| 大香蕉久久网| 啦啦啦在线观看免费高清www| 人妻系列 视频| 久久久久久久久久久免费av| 久久久久人妻精品一区果冻| 久久久国产欧美日韩av| 免费观看a级毛片全部| 亚洲综合精品二区| 国产在线视频一区二区| 国产亚洲91精品色在线| 中文字幕免费在线视频6| 亚洲第一区二区三区不卡| 欧美国产精品一级二级三级 | 国产男人的电影天堂91| 狠狠精品人妻久久久久久综合| 亚洲电影在线观看av| 国产精品一区二区在线不卡| av黄色大香蕉| 欧美日韩视频精品一区| 下体分泌物呈黄色| av福利片在线观看| 一级毛片久久久久久久久女| 国产免费又黄又爽又色| 最新中文字幕久久久久| 欧美xxxx性猛交bbbb| 五月开心婷婷网| 国产精品人妻久久久影院| 天堂俺去俺来也www色官网| 国产白丝娇喘喷水9色精品| 国产国拍精品亚洲av在线观看| 亚洲成人一二三区av| 在线天堂最新版资源| 国产精品欧美亚洲77777| 三上悠亚av全集在线观看 | 亚洲国产欧美日韩在线播放 | 色视频www国产| 精品国产一区二区三区久久久樱花| 丰满乱子伦码专区| av在线老鸭窝| 日韩不卡一区二区三区视频在线| 国产精品一二三区在线看| 最近中文字幕2019免费版| √禁漫天堂资源中文www| 如何舔出高潮| 午夜福利网站1000一区二区三区| 天美传媒精品一区二区| 日韩av在线免费看完整版不卡| 青春草亚洲视频在线观看| 新久久久久国产一级毛片| 日产精品乱码卡一卡2卡三| 欧美日韩国产mv在线观看视频| 久久久欧美国产精品| 免费播放大片免费观看视频在线观看| 免费在线观看成人毛片| 高清黄色对白视频在线免费看 | 欧美日韩一区二区视频在线观看视频在线| 免费观看性生交大片5| 国产高清三级在线| 人人妻人人澡人人看| 黄色视频在线播放观看不卡| 熟女电影av网| 少妇猛男粗大的猛烈进出视频| 精品视频人人做人人爽| 亚洲怡红院男人天堂| 久久人人爽人人片av| 成人亚洲欧美一区二区av| 国产视频首页在线观看| 乱码一卡2卡4卡精品| 欧美日韩在线观看h| 国产精品嫩草影院av在线观看| 成年人免费黄色播放视频 | 91成人精品电影| 国产午夜精品久久久久久一区二区三区| 一本久久精品| 黄色一级大片看看| 9色porny在线观看| 两个人免费观看高清视频 | 国产av一区二区精品久久| 搡女人真爽免费视频火全软件| 久久精品夜色国产| 高清黄色对白视频在线免费看 | 高清黄色对白视频在线免费看 | 纵有疾风起免费观看全集完整版| 一区二区三区乱码不卡18| 色5月婷婷丁香| 国产欧美另类精品又又久久亚洲欧美| 九九爱精品视频在线观看| 亚洲精品视频女| 亚洲精品国产av成人精品| 一个人看视频在线观看www免费| 午夜激情久久久久久久| av黄色大香蕉| 亚洲美女视频黄频| 女的被弄到高潮叫床怎么办| 亚洲激情五月婷婷啪啪| 久久久久精品性色| 亚洲欧美精品自产自拍| 一级毛片我不卡| 亚洲久久久国产精品| 九色成人免费人妻av| 高清毛片免费看| 久久99精品国语久久久| 51国产日韩欧美| 精品国产一区二区三区久久久樱花| 色视频在线一区二区三区| 一级毛片aaaaaa免费看小| 菩萨蛮人人尽说江南好唐韦庄| 日本黄色片子视频| 成人毛片a级毛片在线播放| 美女内射精品一级片tv| 久久鲁丝午夜福利片| 精品一区在线观看国产| 久久狼人影院| 日韩大片免费观看网站| 97在线人人人人妻| 日韩精品有码人妻一区| 久久久久久久大尺度免费视频| 亚洲一区二区三区欧美精品| 男人爽女人下面视频在线观看| 国产高清不卡午夜福利| 在线观看免费高清a一片| 尾随美女入室| 成人影院久久| 久久久欧美国产精品| 国产精品人妻久久久影院| 亚洲不卡免费看| 国产真实伦视频高清在线观看| 国产美女午夜福利| 三级国产精品欧美在线观看| 性色avwww在线观看| 成年人免费黄色播放视频 | 日韩电影二区| 国产精品国产三级国产av玫瑰| 中文字幕免费在线视频6| 亚洲真实伦在线观看| 麻豆成人午夜福利视频| 两个人的视频大全免费| 亚洲精品日本国产第一区| 亚洲成人av在线免费| 免费观看性生交大片5| 蜜桃在线观看..| 国产欧美另类精品又又久久亚洲欧美| 91成人精品电影| 少妇人妻 视频| 亚洲av综合色区一区| 亚洲av中文av极速乱| 国内精品宾馆在线| 你懂的网址亚洲精品在线观看| 国产成人a∨麻豆精品| 麻豆成人av视频| 丰满乱子伦码专区| 另类亚洲欧美激情| 激情五月婷婷亚洲| av专区在线播放| xxx大片免费视频| 中文字幕人妻熟人妻熟丝袜美| 五月伊人婷婷丁香| 在线亚洲精品国产二区图片欧美 | 校园人妻丝袜中文字幕| 不卡视频在线观看欧美| 欧美日韩视频高清一区二区三区二| 另类亚洲欧美激情| 亚洲精品国产色婷婷电影| 日韩精品免费视频一区二区三区 | 亚洲国产精品一区二区三区在线| 春色校园在线视频观看| 美女脱内裤让男人舔精品视频| 日本av手机在线免费观看| 日韩成人伦理影院| 美女xxoo啪啪120秒动态图| 少妇的逼水好多| 大片电影免费在线观看免费| 成年女人在线观看亚洲视频| 国产亚洲午夜精品一区二区久久| 成人午夜精彩视频在线观看| 免费av中文字幕在线| 亚洲三级黄色毛片| 亚洲内射少妇av| 国产一区有黄有色的免费视频| 色婷婷av一区二区三区视频| 国产乱人偷精品视频| 亚洲av成人精品一二三区| 中文字幕制服av| 日本午夜av视频| 国产无遮挡羞羞视频在线观看| a级毛片在线看网站| 国产亚洲精品久久久com| 美女脱内裤让男人舔精品视频| 精品卡一卡二卡四卡免费| 国产白丝娇喘喷水9色精品| 国精品久久久久久国模美| 不卡视频在线观看欧美| 久久 成人 亚洲| 大香蕉久久网| 好男人视频免费观看在线| 亚洲av二区三区四区| 久久韩国三级中文字幕| 黑丝袜美女国产一区| 久久韩国三级中文字幕| 大香蕉久久网| 99九九在线精品视频 | 欧美日韩av久久| 人妻夜夜爽99麻豆av| 五月天丁香电影| 天堂中文最新版在线下载| 美女cb高潮喷水在线观看| 国产有黄有色有爽视频| 国产亚洲精品久久久com| 免费观看在线日韩| 夜夜骑夜夜射夜夜干| 肉色欧美久久久久久久蜜桃| 一区二区三区乱码不卡18| 插逼视频在线观看| 一本一本综合久久| 日韩人妻高清精品专区| 秋霞在线观看毛片| 欧美日韩视频精品一区| 久久久久久人妻| 一区二区av电影网| 久久久a久久爽久久v久久| 爱豆传媒免费全集在线观看| 日韩电影二区| 赤兔流量卡办理| 国产亚洲最大av| 97超碰精品成人国产| 美女xxoo啪啪120秒动态图| 人人妻人人爽人人添夜夜欢视频 | 人妻 亚洲 视频| 亚洲欧洲日产国产| 亚洲欧美成人精品一区二区| 精品少妇久久久久久888优播| 精品人妻熟女毛片av久久网站| 人人妻人人添人人爽欧美一区卜| 成年人免费黄色播放视频 | 伦精品一区二区三区| 在线观看美女被高潮喷水网站| 欧美人与善性xxx| 最近中文字幕高清免费大全6| www.av在线官网国产| 少妇猛男粗大的猛烈进出视频| 亚洲成色77777| 一级毛片我不卡| 精品熟女少妇av免费看| 国产美女午夜福利| 中文欧美无线码| 欧美成人精品欧美一级黄| 日韩精品有码人妻一区| 丝袜喷水一区| 97在线视频观看| 18+在线观看网站| 9色porny在线观看| 亚洲中文av在线| 毛片一级片免费看久久久久| 你懂的网址亚洲精品在线观看| 男人和女人高潮做爰伦理| 亚洲av欧美aⅴ国产| 51国产日韩欧美| 亚洲精品,欧美精品| 欧美成人精品欧美一级黄| 精品国产一区二区久久| 少妇精品久久久久久久| 欧美另类一区| 99国产精品免费福利视频| 日日摸夜夜添夜夜爱| 黑丝袜美女国产一区| 热99国产精品久久久久久7| 国产成人精品无人区| 最近的中文字幕免费完整| 自拍偷自拍亚洲精品老妇| 亚洲精品456在线播放app| 日韩精品有码人妻一区| 午夜影院在线不卡| 成人影院久久| 看免费成人av毛片| 婷婷色综合www| 午夜福利网站1000一区二区三区| 日韩一区二区三区影片| 最新中文字幕久久久久| 在线观看一区二区三区激情| 狂野欧美激情性xxxx在线观看| 欧美+日韩+精品| 久久韩国三级中文字幕| 国产成人精品久久久久久| 水蜜桃什么品种好| 中文欧美无线码| 九九爱精品视频在线观看| 免费不卡的大黄色大毛片视频在线观看| 大又大粗又爽又黄少妇毛片口| 成人特级av手机在线观看| 麻豆成人av视频| 国产精品久久久久久av不卡| 国产免费福利视频在线观看| 人妻系列 视频| 3wmmmm亚洲av在线观看| 久久久久久久亚洲中文字幕| 搡老乐熟女国产| 日本黄色日本黄色录像| 成人无遮挡网站| 寂寞人妻少妇视频99o| 欧美xxxx性猛交bbbb| 亚洲熟女精品中文字幕| 中文在线观看免费www的网站| 亚洲内射少妇av| 高清视频免费观看一区二区| 午夜91福利影院| 永久免费av网站大全| 久久99热6这里只有精品| 日韩精品有码人妻一区| 少妇 在线观看| 亚洲一区二区三区欧美精品| 国产欧美日韩一区二区三区在线 | 高清欧美精品videossex| 美女大奶头黄色视频| 亚洲成人一二三区av| 亚洲伊人久久精品综合| 久久婷婷青草| av在线播放精品| 日本免费在线观看一区| 麻豆乱淫一区二区| 亚洲国产av新网站| a级一级毛片免费在线观看| 亚洲精品久久午夜乱码| 高清午夜精品一区二区三区| 99热全是精品| 日韩不卡一区二区三区视频在线| 大又大粗又爽又黄少妇毛片口| 亚洲情色 制服丝袜| av天堂中文字幕网| 久久精品夜色国产| 亚洲欧洲国产日韩| 精品一品国产午夜福利视频| 国产高清有码在线观看视频| 亚洲欧美一区二区三区黑人 | 亚洲国产精品专区欧美| 少妇的逼水好多| 高清不卡的av网站| 国产免费一级a男人的天堂| 啦啦啦啦在线视频资源| 最近最新中文字幕免费大全7| 少妇的逼水好多| 国产免费视频播放在线视频| 欧美精品一区二区大全| 欧美日韩亚洲高清精品| 中文精品一卡2卡3卡4更新| 亚洲性久久影院| 最后的刺客免费高清国语| 亚洲精品乱码久久久久久按摩| 久久久久久久国产电影| 亚洲内射少妇av| 亚洲国产日韩一区二区| 五月天丁香电影| 国产一区二区三区av在线| 久久久国产欧美日韩av| 偷拍熟女少妇极品色| 人体艺术视频欧美日本| 青春草国产在线视频| 国产伦精品一区二区三区视频9| 亚洲丝袜综合中文字幕| 日日爽夜夜爽网站| 亚洲av欧美aⅴ国产| 亚洲国产成人一精品久久久| 亚洲内射少妇av| www.av在线官网国产| 国产片特级美女逼逼视频| 免费观看a级毛片全部| 午夜影院在线不卡| 精品少妇久久久久久888优播| 日本-黄色视频高清免费观看| 欧美丝袜亚洲另类| 亚洲自偷自拍三级| 亚洲欧美日韩东京热| 五月天丁香电影| 观看av在线不卡| 精品一区二区三卡| 亚洲,一卡二卡三卡| 两个人的视频大全免费| 亚洲精品国产av蜜桃| 午夜久久久在线观看| 欧美激情极品国产一区二区三区 | 亚洲综合精品二区| av播播在线观看一区| 又黄又爽又刺激的免费视频.| 亚洲情色 制服丝袜| 精品熟女少妇av免费看| 91精品国产九色| 国产国拍精品亚洲av在线观看| 黄色配什么色好看| 国产成人aa在线观看| 少妇人妻 视频| 久久毛片免费看一区二区三区| 久久午夜福利片| 久久久午夜欧美精品| 涩涩av久久男人的天堂| 国产综合精华液| 日本黄大片高清| 久久久久久久大尺度免费视频| 一边亲一边摸免费视频| 三级经典国产精品| 免费在线观看成人毛片| 一本色道久久久久久精品综合| 在线精品无人区一区二区三| 热99国产精品久久久久久7| 日本午夜av视频| h视频一区二区三区| 99久久综合免费| 国产极品天堂在线| 乱码一卡2卡4卡精品| 色视频在线一区二区三区| 成人毛片60女人毛片免费| 成年av动漫网址| 欧美少妇被猛烈插入视频| 亚洲成人av在线免费| 免费看av在线观看网站| 2018国产大陆天天弄谢| 亚洲av二区三区四区| 三级国产精品片| 国产精品人妻久久久影院| www.av在线官网国产| 久久精品夜色国产| 丰满人妻一区二区三区视频av| 99国产精品免费福利视频| 国产淫片久久久久久久久| 精品国产一区二区久久| 久久亚洲国产成人精品v| 欧美 亚洲 国产 日韩一| av福利片在线| 另类亚洲欧美激情| 高清欧美精品videossex| av福利片在线| 国产精品国产av在线观看| 亚洲综合精品二区| 3wmmmm亚洲av在线观看| 国产午夜精品一二区理论片| 国语对白做爰xxxⅹ性视频网站| 亚洲熟女精品中文字幕| a级一级毛片免费在线观看| 精品久久久精品久久久| 自拍偷自拍亚洲精品老妇| 精品一区二区三卡| 亚洲婷婷狠狠爱综合网| 日本av免费视频播放| 午夜免费观看性视频| 人妻人人澡人人爽人人| av免费观看日本| 少妇的逼水好多| 欧美 亚洲 国产 日韩一| 国内揄拍国产精品人妻在线| 久久亚洲国产成人精品v| 精品一区二区三区视频在线| 极品人妻少妇av视频| 亚洲久久久国产精品| 国产精品99久久99久久久不卡 | 男女国产视频网站| 性高湖久久久久久久久免费观看| 久久精品国产亚洲网站| 啦啦啦在线观看免费高清www| 国产精品人妻久久久影院| 伦理电影大哥的女人| 久久99热6这里只有精品| 九色成人免费人妻av|