• <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ù)值模擬
    19禁男女啪啪无遮挡网站| 757午夜福利合集在线观看| 91九色精品人成在线观看| 最新在线观看一区二区三区| 欧美日韩一级在线毛片| e午夜精品久久久久久久| 亚洲欧美一区二区三区黑人| 国产精品98久久久久久宅男小说| 久久草成人影院| 国内精品久久久久久久电影| 变态另类丝袜制服| 午夜影院日韩av| 欧美乱色亚洲激情| 欧美乱色亚洲激情| 两性午夜刺激爽爽歪歪视频在线观看 | АⅤ资源中文在线天堂| 亚洲国产精品sss在线观看| 999久久久精品免费观看国产| 校园春色视频在线观看| 色尼玛亚洲综合影院| 狂野欧美激情性xxxx| 欧美成人性av电影在线观看| 变态另类丝袜制服| 免费电影在线观看免费观看| 亚洲一区中文字幕在线| 亚洲一区二区三区不卡视频| 精品国产国语对白av| 日日夜夜操网爽| 欧美国产日韩亚洲一区| 午夜精品在线福利| 国产麻豆成人av免费视频| 成人三级做爰电影| 老司机福利观看| 精品国产乱子伦一区二区三区| 视频在线观看一区二区三区| 大型av网站在线播放| 午夜精品在线福利| 午夜久久久在线观看| 很黄的视频免费| 最近最新中文字幕大全免费视频| 男女之事视频高清在线观看| 一区二区三区精品91| 国产色视频综合| 中文字幕人妻熟女乱码| 欧美人与性动交α欧美精品济南到| 午夜精品在线福利| 给我免费播放毛片高清在线观看| 在线观看午夜福利视频| 欧美色视频一区免费| 久久久国产成人精品二区| 久久人妻av系列| 精品免费久久久久久久清纯| 欧美日韩乱码在线| 美女午夜性视频免费| 日本精品一区二区三区蜜桃| 亚洲av电影在线进入| 亚洲激情在线av| 欧美日本视频| 亚洲专区国产一区二区| 国产精品综合久久久久久久免费| 久久国产乱子伦精品免费另类| 久久午夜综合久久蜜桃| 久9热在线精品视频| 亚洲欧美激情综合另类| 中出人妻视频一区二区| 亚洲第一电影网av| 免费看日本二区| 久久午夜综合久久蜜桃| 亚洲专区字幕在线| 一边摸一边做爽爽视频免费| 亚洲 欧美 日韩 在线 免费| 午夜福利高清视频| 欧美在线黄色| 少妇被粗大的猛进出69影院| 久久精品国产亚洲av香蕉五月| 亚洲激情在线av| 欧美黑人精品巨大| 日韩欧美三级三区| 波多野结衣av一区二区av| 国产精品久久久人人做人人爽| www日本黄色视频网| 国内少妇人妻偷人精品xxx网站 | av福利片在线| 久久久久九九精品影院| 女人被狂操c到高潮| 亚洲一卡2卡3卡4卡5卡精品中文| 十八禁人妻一区二区| 亚洲国产精品成人综合色| 男女下面进入的视频免费午夜 | 99在线视频只有这里精品首页| 别揉我奶头~嗯~啊~动态视频| 1024手机看黄色片| 99精品久久久久人妻精品| 99国产精品一区二区蜜桃av| 免费女性裸体啪啪无遮挡网站| 久久久久久久久久黄片| 2021天堂中文幕一二区在线观 | 制服诱惑二区| 亚洲精品美女久久av网站| 欧美激情 高清一区二区三区| 两个人视频免费观看高清| 一个人观看的视频www高清免费观看 | 精品国产美女av久久久久小说| 亚洲国产欧洲综合997久久, | 韩国av一区二区三区四区| 久久香蕉精品热| 午夜福利18| 日韩三级视频一区二区三区| 日本a在线网址| 欧美日韩一级在线毛片| 伦理电影免费视频| 国产精品乱码一区二三区的特点| 身体一侧抽搐| 搡老岳熟女国产| 成人亚洲精品一区在线观看| www.www免费av| 国产视频内射| 可以在线观看毛片的网站| 亚洲成人免费电影在线观看| 久久久精品欧美日韩精品| 少妇的丰满在线观看| 岛国在线观看网站| 99在线人妻在线中文字幕| 免费在线观看亚洲国产| 国产精品一区二区精品视频观看| av有码第一页| 日韩欧美在线二视频| av免费在线观看网站| 欧美激情 高清一区二区三区| av福利片在线| 国产伦在线观看视频一区| 日韩精品免费视频一区二区三区| 午夜精品在线福利| 不卡av一区二区三区| 亚洲片人在线观看| 国产精品美女特级片免费视频播放器 | 日本一区二区免费在线视频| 欧美绝顶高潮抽搐喷水| 99热只有精品国产| 久久精品国产99精品国产亚洲性色| 757午夜福利合集在线观看| 午夜免费观看网址| 亚洲欧美激情综合另类| 婷婷亚洲欧美| 国产一区二区激情短视频| 国产区一区二久久| 久久精品91蜜桃| 女性生殖器流出的白浆| 中出人妻视频一区二区| 在线观看66精品国产| 哪里可以看免费的av片| 欧美色欧美亚洲另类二区| 激情在线观看视频在线高清| 女生性感内裤真人,穿戴方法视频| 首页视频小说图片口味搜索| 不卡一级毛片| 亚洲av熟女| 久久伊人香网站| 黄色 视频免费看| 极品教师在线免费播放| 亚洲第一电影网av| 最好的美女福利视频网| 91国产中文字幕| 9191精品国产免费久久| 亚洲成人国产一区在线观看| 国产精品亚洲av一区麻豆| 一级a爱视频在线免费观看| www.999成人在线观看| 精品无人区乱码1区二区| 少妇熟女aⅴ在线视频| 久久精品国产清高在天天线| 最新在线观看一区二区三区| 亚洲精品中文字幕在线视频| 啪啪无遮挡十八禁网站| 一区二区三区国产精品乱码| 精品国产一区二区三区四区第35| 别揉我奶头~嗯~啊~动态视频| 久久精品亚洲精品国产色婷小说| 十八禁网站免费在线| 欧美色欧美亚洲另类二区| 国产精品1区2区在线观看.| 久久精品影院6| 十分钟在线观看高清视频www| 国产精品免费视频内射| 免费观看人在逋| 亚洲欧洲精品一区二区精品久久久| 国产爱豆传媒在线观看 | 欧美中文日本在线观看视频| 亚洲国产欧美日韩在线播放| 国产精品 国内视频| 久久青草综合色| 国产成人系列免费观看| 在线观看免费午夜福利视频| 变态另类成人亚洲欧美熟女| 俄罗斯特黄特色一大片| 成人一区二区视频在线观看| 少妇被粗大的猛进出69影院| 免费看日本二区| 国产99白浆流出| 黄频高清免费视频| 麻豆一二三区av精品| 99re在线观看精品视频| 国产国语露脸激情在线看| 日韩精品中文字幕看吧| 久久久久久大精品| xxx96com| 搡老岳熟女国产| 久久午夜亚洲精品久久| 精品福利观看| 久久中文字幕人妻熟女| 国产99久久九九免费精品| www日本黄色视频网| 午夜精品久久久久久毛片777| 欧美 亚洲 国产 日韩一| 亚洲第一青青草原| 美女扒开内裤让男人捅视频| 国产精品久久久久久亚洲av鲁大| 午夜精品在线福利| 国产黄片美女视频| 中文字幕另类日韩欧美亚洲嫩草| 中文字幕av电影在线播放| 欧美黑人精品巨大| bbb黄色大片| 在线十欧美十亚洲十日本专区| www.999成人在线观看| 国产v大片淫在线免费观看| 999久久久精品免费观看国产| а√天堂www在线а√下载| 精品国产乱子伦一区二区三区| 精华霜和精华液先用哪个| 国产欧美日韩一区二区精品| 国产精品野战在线观看| 久久精品人妻少妇| 色综合站精品国产| 无人区码免费观看不卡| 狂野欧美激情性xxxx| 一级毛片女人18水好多| 亚洲精品在线美女| 亚洲第一青青草原| bbb黄色大片| 老司机福利观看| 在线十欧美十亚洲十日本专区| 97人妻精品一区二区三区麻豆 | 亚洲一区高清亚洲精品| 波多野结衣巨乳人妻| 日韩免费av在线播放| 男女下面进入的视频免费午夜 | 夜夜爽天天搞| 成人三级做爰电影| 亚洲国产精品合色在线| 日韩欧美一区视频在线观看| 欧美zozozo另类| 高清在线国产一区| 精品少妇一区二区三区视频日本电影| 午夜福利一区二区在线看| 一边摸一边抽搐一进一小说| 精品第一国产精品| 国产精品1区2区在线观看.| 精品乱码久久久久久99久播| 男人操女人黄网站| 日日爽夜夜爽网站| 欧美一区二区精品小视频在线| 99国产精品一区二区蜜桃av| 日韩国内少妇激情av| 在线天堂中文资源库| www.熟女人妻精品国产| 男女下面进入的视频免费午夜 | 黄色 视频免费看| 日本 欧美在线| 国产伦一二天堂av在线观看| 成人亚洲精品一区在线观看| 黄片播放在线免费| 国产一卡二卡三卡精品| 欧美国产日韩亚洲一区| 日韩精品免费视频一区二区三区| 欧美成人免费av一区二区三区| 男人舔奶头视频| 成人手机av| 欧美激情久久久久久爽电影| 国产黄色小视频在线观看| 一级毛片高清免费大全| 国产精品免费一区二区三区在线| 人妻久久中文字幕网| 日本五十路高清| 99在线视频只有这里精品首页| 麻豆成人av在线观看| 中文字幕另类日韩欧美亚洲嫩草| 国产熟女xx| 亚洲精品国产精品久久久不卡| 国产人伦9x9x在线观看| 啦啦啦免费观看视频1| 亚洲欧美激情综合另类| 午夜福利在线在线| bbb黄色大片| 亚洲中文字幕一区二区三区有码在线看 | 亚洲欧美日韩无卡精品| 国产日本99.免费观看| 国产一区二区三区视频了| 男女视频在线观看网站免费 | 18禁观看日本| 久久香蕉激情| 国产一级毛片七仙女欲春2 | 99国产极品粉嫩在线观看| tocl精华| 麻豆久久精品国产亚洲av| 丝袜人妻中文字幕| 日本成人三级电影网站| 在线观看午夜福利视频| 国产亚洲av高清不卡| 免费人成视频x8x8入口观看| 老熟妇乱子伦视频在线观看| 色播亚洲综合网| 国产精品亚洲美女久久久| 亚洲av日韩精品久久久久久密| 黄色丝袜av网址大全| av天堂在线播放| 久久精品国产亚洲av高清一级| 美女 人体艺术 gogo| 在线国产一区二区在线| 久久国产乱子伦精品免费另类| 久久精品91无色码中文字幕| 欧美最黄视频在线播放免费| 午夜福利在线在线| 俄罗斯特黄特色一大片| 中文字幕精品亚洲无线码一区 | av免费在线观看网站| 亚洲一区二区三区色噜噜| 18禁观看日本| 午夜福利免费观看在线| 99久久无色码亚洲精品果冻| 成人18禁在线播放| 国产亚洲精品久久久久久毛片| АⅤ资源中文在线天堂| 久久久久国产一级毛片高清牌| 一边摸一边做爽爽视频免费| 亚洲狠狠婷婷综合久久图片| 国产亚洲欧美精品永久| 少妇粗大呻吟视频| 91av网站免费观看| 日本一本二区三区精品| 精品人妻1区二区| 日韩高清综合在线| 亚洲av成人一区二区三| 国产主播在线观看一区二区| 操出白浆在线播放| 国产精品国产高清国产av| 99精品欧美一区二区三区四区| 亚洲成a人片在线一区二区| 婷婷亚洲欧美| 日韩三级视频一区二区三区| 精品第一国产精品| 亚洲成av片中文字幕在线观看| 国产aⅴ精品一区二区三区波| 草草在线视频免费看| 精品日产1卡2卡| 波多野结衣av一区二区av| 国产区一区二久久| 18禁黄网站禁片免费观看直播| 欧美中文日本在线观看视频| 99re在线观看精品视频| 国产久久久一区二区三区| 一区福利在线观看| 黄色丝袜av网址大全| 午夜福利一区二区在线看| 亚洲国产精品久久男人天堂| 久久精品国产清高在天天线| 91老司机精品| 啦啦啦 在线观看视频| 国产精品电影一区二区三区| 欧美日韩中文字幕国产精品一区二区三区| 色综合婷婷激情| a级毛片a级免费在线| 色综合婷婷激情| 亚洲av电影在线进入| 人人妻人人看人人澡| 国产精品爽爽va在线观看网站 | 亚洲成人国产一区在线观看| 我的亚洲天堂| 久久久久久国产a免费观看| 看片在线看免费视频| 女性生殖器流出的白浆| 欧美乱码精品一区二区三区| 国产黄色小视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美zozozo另类| 曰老女人黄片| 18禁裸乳无遮挡免费网站照片 | 亚洲一码二码三码区别大吗| 欧美不卡视频在线免费观看 | 黄色丝袜av网址大全| 香蕉av资源在线| 观看免费一级毛片| 啦啦啦 在线观看视频| 久久亚洲精品不卡| 丰满人妻熟妇乱又伦精品不卡| 欧美午夜高清在线| 我的亚洲天堂| 久久久久久大精品| 欧美绝顶高潮抽搐喷水| 男女视频在线观看网站免费 | 国产精品,欧美在线| 精品第一国产精品| 亚洲精品国产精品久久久不卡| 成人精品一区二区免费| 少妇裸体淫交视频免费看高清 | 欧美三级亚洲精品| 亚洲成人久久爱视频| 国产精品久久久久久亚洲av鲁大| 成年女人毛片免费观看观看9| 精华霜和精华液先用哪个| 亚洲男人的天堂狠狠| 成人三级黄色视频| 久久久精品国产亚洲av高清涩受| 亚洲一区二区三区色噜噜| 欧美黑人巨大hd| 国产黄片美女视频| 国产亚洲av嫩草精品影院| 亚洲国产欧美日韩在线播放| 老司机午夜十八禁免费视频| 国产真人三级小视频在线观看| 亚洲国产精品久久男人天堂| 男女床上黄色一级片免费看| 国内精品久久久久精免费| av有码第一页| 欧美绝顶高潮抽搐喷水| 中国美女看黄片| 在线观看www视频免费| 在线观看日韩欧美| 男人舔女人的私密视频| 十八禁人妻一区二区| 亚洲av五月六月丁香网| 午夜福利免费观看在线| 最近最新中文字幕大全电影3 | 久久这里只有精品19| 日本五十路高清| 99久久综合精品五月天人人| 日韩欧美一区二区三区在线观看| 男女做爰动态图高潮gif福利片| 99热这里只有精品一区 | 精品国产亚洲在线| 欧美日韩福利视频一区二区| 国产亚洲精品第一综合不卡| 美女扒开内裤让男人捅视频| 亚洲男人天堂网一区| 久久 成人 亚洲| 极品教师在线免费播放| 国产精品99久久99久久久不卡| 亚洲五月婷婷丁香| 波多野结衣高清无吗| www国产在线视频色| 午夜福利一区二区在线看| 淫妇啪啪啪对白视频| 视频区欧美日本亚洲| 老司机深夜福利视频在线观看| 欧美zozozo另类| 熟女电影av网| 日韩有码中文字幕| 欧美av亚洲av综合av国产av| 欧美国产日韩亚洲一区| 国产黄a三级三级三级人| 黑人操中国人逼视频| 亚洲一码二码三码区别大吗| 精品国产亚洲在线| 国产日本99.免费观看| 国产一区在线观看成人免费| 日本一本二区三区精品| 黑丝袜美女国产一区| 少妇熟女aⅴ在线视频| 十八禁人妻一区二区| 18禁观看日本| 成熟少妇高潮喷水视频| 岛国视频午夜一区免费看| 中文字幕高清在线视频| 亚洲精品中文字幕一二三四区| 精品少妇一区二区三区视频日本电影| 免费在线观看黄色视频的| 久久国产乱子伦精品免费另类| 日本在线视频免费播放| 久久香蕉国产精品| 亚洲va日本ⅴa欧美va伊人久久| 久久精品国产亚洲av高清一级| 亚洲一区中文字幕在线| а√天堂www在线а√下载| 在线十欧美十亚洲十日本专区| 男人舔女人下体高潮全视频| 免费看美女性在线毛片视频| av有码第一页| www.熟女人妻精品国产| 黄片小视频在线播放| 一本精品99久久精品77| 久久天堂一区二区三区四区| 国产亚洲欧美精品永久| 国内揄拍国产精品人妻在线 | 少妇熟女aⅴ在线视频| 久久精品国产亚洲av高清一级| 欧美国产精品va在线观看不卡| 午夜亚洲福利在线播放| netflix在线观看网站| 久久久久久大精品| 国产成人精品久久二区二区免费| av有码第一页| 成人一区二区视频在线观看| 亚洲精品美女久久久久99蜜臀| 午夜久久久久精精品| 香蕉av资源在线| 色尼玛亚洲综合影院| 国产精品久久久人人做人人爽| 亚洲国产精品久久男人天堂| 日本 av在线| 久久伊人香网站| 两个人免费观看高清视频| 熟女少妇亚洲综合色aaa.| 亚洲国产看品久久| 精品久久久久久成人av| 久热这里只有精品99| 美女扒开内裤让男人捅视频| 一区二区日韩欧美中文字幕| 午夜久久久在线观看| 国产av一区二区精品久久| 欧美中文综合在线视频| 日韩欧美一区视频在线观看| svipshipincom国产片| 欧美午夜高清在线| 欧美久久黑人一区二区| 久久精品影院6| 国产精品电影一区二区三区| 精品免费久久久久久久清纯| 久久婷婷成人综合色麻豆| 中出人妻视频一区二区| 久久精品影院6| 少妇熟女aⅴ在线视频| 好男人在线观看高清免费视频 | 最近最新中文字幕大全电影3 | 久久国产亚洲av麻豆专区| 亚洲国产精品成人综合色| 欧美一区二区精品小视频在线| 一卡2卡三卡四卡精品乱码亚洲| 日韩三级视频一区二区三区| 激情在线观看视频在线高清| 国产极品粉嫩免费观看在线| 草草在线视频免费看| 人人妻人人看人人澡| 亚洲精品久久国产高清桃花| 午夜福利在线在线| 欧美大码av| 亚洲九九香蕉| 成人特级黄色片久久久久久久| 亚洲成a人片在线一区二区| 国产精品,欧美在线| 久久精品亚洲精品国产色婷小说| 成年人黄色毛片网站| 精华霜和精华液先用哪个| 国产99久久九九免费精品| 国产成人精品无人区| 12—13女人毛片做爰片一| 成人三级做爰电影| av有码第一页| 一级黄色大片毛片| 免费在线观看亚洲国产| 欧美激情 高清一区二区三区| 熟妇人妻久久中文字幕3abv| tocl精华| 亚洲精品久久成人aⅴ小说| 看黄色毛片网站| 1024视频免费在线观看| 亚洲真实伦在线观看| 日本一本二区三区精品| av视频在线观看入口| 久久久久久九九精品二区国产 | 在线观看66精品国产| 99久久精品国产亚洲精品| 88av欧美| 美女国产高潮福利片在线看| 黄色片一级片一级黄色片| 午夜免费观看网址| 亚洲自偷自拍图片 自拍| 国产真实乱freesex| av在线播放免费不卡| 看黄色毛片网站| 高清毛片免费观看视频网站| 97人妻精品一区二区三区麻豆 | 岛国视频午夜一区免费看| 免费观看精品视频网站| 亚洲国产精品成人综合色| 国产精品久久久久久亚洲av鲁大| 两个人看的免费小视频| 成人午夜高清在线视频 | 亚洲成av人片免费观看| 欧美又色又爽又黄视频| 国产精品永久免费网站| 日韩欧美免费精品| 免费在线观看影片大全网站| svipshipincom国产片| 免费看美女性在线毛片视频| 少妇粗大呻吟视频| 国产片内射在线| 搡老妇女老女人老熟妇| 大型黄色视频在线免费观看| 亚洲天堂国产精品一区在线| 搡老妇女老女人老熟妇| 久久精品人妻少妇| 欧美久久黑人一区二区| 日韩av在线大香蕉| 日韩精品青青久久久久久| xxx96com| 国产精品久久久久久人妻精品电影| 久久久久亚洲av毛片大全| 国产在线精品亚洲第一网站| 亚洲av熟女| 日本五十路高清| 久久青草综合色| www.www免费av| 日韩 欧美 亚洲 中文字幕| netflix在线观看网站| 最近在线观看免费完整版|