甲醇制烯烴反應動力學及反應器模型研究進展

文堯順，南海明，吳秀章，關豐忠，孫保全

（1中國神華煤制油化工有限公司，北京 100011；2中國石油大學（北京）重質油國家重點實驗室，北京 102249）

特約評述

甲醇制烯烴反應動力學及反應器模型研究進展

文堯順1，2，南海明1，吳秀章1，關豐忠1，孫保全1

（1中國神華煤制油化工有限公司，北京 100011；2中國石油大學（北京）重質油國家重點實驗室，北京 102249）

甲醇制烯烴（MTO）工藝是現代煤化工領域的研究熱點，MTO反應動力學及其反應器模型研究是高效反應器開發和工業裝置操作優化的基礎。本文綜述了甲醇制烯烴反應動力學研究進展，詳細論述了機理型動力學模型、八集總動力學模型、五集總動力學模型，指出集總動力學模型適用于描述MTO反應過程，如何考慮水、積炭等因素的影響是MTO動力學研究的難點和關鍵。結合現有動力學模型，評述了MTO反應過程在工業規模的固定床反應器、提升管反應器、循環流化床反應器、湍動流化床反應器中產物分布和轉化率模擬情況，結果表明：循環流化床反應器和湍動流化床反應器適合MTO工業過程。最后指出，甲醇制烯烴反應動力學下一步研究方應集中于工業規模流化床反應器氣固兩相流動模擬，以及與動力學模型結合獲得準確預測工業反應結果的MTO反應器模型。

甲醇；烯烴；反應動力學；反應器；模型

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

目前，科學界對甲醇制烯烴催化劑[8-9]、反應特性[10]、工藝技術[11]、催化劑積炭[12-13]等方面研究總結較多，但是關于反應動力學研究進展的評述鮮有報道。動力學研究不僅是探究MTO過程反應機理的重要手段，對高效催化劑研發的具有指導意義，而且也是新型反應器和工藝開發的重要基礎。本文系統總結甲醇制烯烴過程反應動力學研究的最新進展，并結合動力學研究對不同類型MTO反應器模擬情況進行了介紹和評述。

1 MTO反應機理

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

根據烴池反應機理，類似積炭的烴池活性物種隨著反應進行，活性逐漸降低，最終形成積炭。積炭一方面影響小孔SAPO-34分子篩形狀擇型性，導致產物選擇性發生變化；另一方面積炭會堵塞部分催化劑孔道，造成催化劑活性降低[19]，MTO總反應速率隨著積炭含量而變化。因此，如何考慮積炭影響對MTO動力學研究十分重要。

2 MTO反應動力學

根據目前甲醇制烯烴動力學的研究進展，可以將動力學研究分成兩大類[20]：一類是從反應機理出發，列出MTO反應中可能出現的所有基元反應，通過實驗、量子理論計算得到的機理型反應動力學；另一類是采用集總方法，忽略基元反應步驟，獲得的集總反應動力學模型。機理型反應動力學能夠描述反應歷程和細節，而集總動力學模型則在反應器設計方面實用性更強。表1匯總了具有代表性甲醇制烯烴反應動力學研究結果，下面分別進行論述。

表1 甲醇制烯烴反應動力學研究匯總

2.1 機理型反應動力學

1983年，Mihail等[21]報道了甲醇轉化為碳氫化合物的反應動力學模型，基于碳烯反應機理，建立了包括53個基元反應的甲醇烴類化合物的反應動力學方程組，將基元反應分成12種主要的反應步驟，獲得了動力學參數，并在固定床反應器上進行數學模擬。Park 和Froment[22-23]根據正碳離子機理，建立了包括氫轉移、甲基化、β消去等726個基元反應和225種反應中間體的動力學模型。采用Hougen-Watson模型描述從二甲醚形成到各種產物生成的步驟，烯烴的生成反應采用正碳離子機理描述。反應網絡利用計算算法生成，各個基元反應的速率常數采用單事件方法求取。基元反應的活化能通過Evans-Polanyi關聯正碳離子和烯烴異構體的能級獲得。Forment等最終獲得33個參數的機理型動力學模型，并對絕熱固定床、等溫固定床、以及流化床反應器進行模擬計算[24]。Forment的甲醇制烯烴動力學參數基于HZSM-5分子篩催化劑實驗結果和正碳離子機理獲得，與業界廣泛接受的烴池反應機理不同，影響了模型的實用性。Alwahabi隨后也報道了基于SAPO-34分子篩的MTO單事件反應動力學模型。

反應機理型動力學模型涉及反應較多，對反應過程的單個反應進行分析，對揭示反應機理和反應過程規律具有重要意義。但是模型預測產物分布的可靠性和準確性還有待于進一步證實和完善。

2.2 集總動力學

集總動力學模型是根據分子的動力學特性，將反應體系劃分為若干個集總組分，建立集總動力學模型。集總動力學模型對于處理復雜反應體系具有廣泛應用，例如煉油工業中的催化裂化和催化重整等采用集總動力學模型取得了成功。與催化裂化相比，MTO反應過程反應原料為單一甲醇組分，相對比較簡單，但是反應產物包含C1～C6組分，反應網絡比較復雜，也適用于采用集總方法來描述。

2.2.1 八集總動力學模型

圖1 Bos甲醇制烯烴反應八集總動力學網絡[25]

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

Bos的動力學模型包含組分較多，研究系統全面，也是后續研究者采用較多的動力學模型，具有重要的借鑒意義。但是該模型未考慮反應溫度和水含量等重要因素對反應結果的影響。模型僅給出了單個反應溫度下（450℃）的速率常數，未給出各個反應的活化能，難以外推其他反應溫度條件下的反應速率。

圖2 Chen甲醇制烯烴反應八集總動力學網絡[26]

2007年，挪威科技大學的Chen研究小組等[26]在Hydro公司資助下，基于電子震蕩天平（TEOM）反應器實驗數據建立了MTO八集總反應動力模型，同樣根據烴池反應機理，但是反應網絡集總劃分與Bos的模型不同。如圖2所示，研究認為甲烷的生成機理與其他產物不同，甲烷的產生來源于催化劑表面的甲氧基分解產生，單獨作為集總。認為所有的烯烴都由二甲醚直接產生，甲醇迅速轉化為二甲醚且其濃度難以測定，將甲醇和二甲醚作為一個集總。甲醇和二甲醚作為含氧化合物集總的劃分方法得到很多研究者的采用。八集總包括甲烷、乙烯、丙烯、丁烯、烷烴（乙烷+丙烷）、C5、C6、和焦炭。所有的烯烴生成反應都是二甲醚的一級平行反應，烯烴發生二次反應生成烷烴（乙烷+丙烷）。積炭對各個反應速率常數的影響，擬合得到失活函數與焦炭含量呈線性關系，失活速率常數隨反應產物分子碳數增加而減小。模型求解獲得乙烯、丙烯、丁烯、C5、C6、乙烷+丙烷反應動力學速率常數、失活速率常數、指前因子和活化能等動力學參數。Chen動力學模型得到乙烯、丙烯和丁烯的反應活化能分別為38.4kJ/mol、27.0kJ/mol和26.9kJ/mol?；罨芟鄬ζ渌芯空邎蟮澜Y果低，Chen認為是由于受晶內擴散的限制[27]。

需要指出的是Chen的動力學反應網絡中給出了甲烷和焦炭兩種反應產物，由于甲烷含量很低而且重復性較差，未能給出兩種反應產物的動力學參數。研究者對Chen實驗數據和參數進行重新回歸擬合后，發現數據的相關性差，模型的實用性還有待于進一步驗證。

2.2.2 五集總動力學

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

圖3 Gayubo MTO反應動力學模型反應網絡[28]

國內研究者也在甲醇制烯烴反應動力學方面開展了大量的研究，取得了不錯的進展。中國石化的齊國禎對MTO反應動力學進行了系統和全面的研究，提供了豐富的動力學反應數據，主要研究結果見文獻[31-32]。齊國禎[32]對五集總動力學模型進行了修正，將丁烯集總變為C4集總，將C5+集總調整為烷烴集總（包括乙烷、丙烷和C5+），如圖4所示。同時選擇反應為1.5min的反應結果進行動力學數據擬合，更能準確的反映初始階段反應特性。模型認為反應過程中生成的水對主反應和副反應均有貢獻，在反應中引入水的吸附阻力項，定量分析水對MTO各個反應的影響程度。

圖4 齊國禎MTO五集總動力學模型反應網絡[32]

圖5 齊國禎MTO五集總動力學模型速率方程[32]

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

齊國禎的MTO動力學模型考慮了水的影響、反應溫度的影響以及催化劑積炭失活的影響。

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

表2 胡浩五集總動力學模型參數[34]

胡浩動力學模型的數據基礎為乙烯加丙烯選擇性最高的時間段，工業MTO裝置操作基本在上述階段，因而采用胡浩的五集總動力學模型應用于工業MTO裝置的產物分布預測，相對誤差可能較小。

3 MTO動力學在反應器模擬方面的研究

在MTO反應動力學研究基礎上，結合反應器內催化劑流動特性、氣固傳質、傳熱等進行MTO反應器模擬研究。一方面給出不同種類反應器MTO反應轉化率和產物分布，優選適合于甲醇轉化生產低碳烯烴工業過程的理想反應器；另一方面，也為工業MTO反應器的設計和MTO裝置的優化操作提供基礎數據。業界在MTO反應器模型研究方面開展了大量的工作。

3.1 固定床反應器

Bos[25]和Alwahabi等[35]模擬研究表明固定床不宜作為MTO反應器，主要有以下4個方面的缺陷：①MTO反應放熱量大，固定床反應器取熱困難；②催化劑失活較快，需要多個反應器進行切換；③固定床反應器再生困難，容易出現飛溫；④反應器壓降大。此外，工業規模固定床反應器若采用等溫操作，采用列管式反應器造價高，裝置投資大；當采用多層絕熱反應器操作時，設置層間取熱和多反應器切換，工藝較為復雜。

3.2 流化床反應器

流化床反應器能夠實現催化劑的連續反應和再生，反應器和再生器取熱方便。目前，工業化的MTO技術均采用流化床工藝。研究者對不同流化床反應器如提升管反應器、循環流化床反應器、湍動流化床反應器進行了MTO反應模擬。

3.2.1 提升管反應器

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

3.2.2 循環流化床反應器

循環流化床反應器包括密相提升管反應器部分和用于換熱的循環部分，Bos 模擬結果表明[25]，當催化劑停留時間為7min時，催化劑平均積炭含量為8.2%（質量分數），乙烯丙烯比例接近1。當停留時間進一步增大時，乙烯/丙烯比增大，反應器體積也明顯增大。Soundararajan等[36]也研究了MTO在循環流化床反應器的產物分布，由于催化劑停留時間短，單程積炭較少，經旋風分離的催化劑大部分返回提升管底部，僅有少部分再生，從而提高提升管內催化劑的平均積炭含量。研究發現提升管出口結構形式和積炭含量對產物分布有較大的影響，當積炭含量為5%時，反應產物低碳烯烴選擇性最高。該研究結果對于循環流化床反應器設計和操作條件選取具有重要指導意義。Schoenfelder等[37]也在MTO循環流化床反應器模擬開展了研究工作。

3.2.3 湍動流化床

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

Bos等[25]湍動流化床模擬結果證實，對于30萬噸/年乙烯生產規模，湍動流化床反應器直徑為8.5m即可滿足要求，反應氣速為1m/s，低于循環流化床的3m/s。隨著停留時間延長，丙烯選擇性在初始階段略有下降，然后基本保持不變，而乙烯選擇性則逐漸增大。

鼓泡流化床操作氣速低于湍動流化床，學者也對模擬了MTO鼓泡流化床反應性能，比較有代表性的是Alwahabi等[35]和齊國禎[32]。前者模擬了年產55000噸乙烯和丙烯的生產規模的MTO裝置流化床反應結果，進料溫度為430℃，出口溫升21℃，甲醇轉化率為90%，乙烯丙烯收率相當于進料甲醇的12%和13%（質量分數）。需要指出的是，前述所有的MTO動力學基礎數據均是基于固定床或者微反試驗取得的，Kaarsholm等[42]報道了基于流化床反應器試驗數據的動力學MTO動力學研究結果，對于指導工業MTO反應器的設計具有意義。但是研究采用的催化劑為磷改性的ZSM-5催化劑，影響了其實用性。

4 結語和展望

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

符 號 說 明

Ea,i——反應活化能，kJ/mol

FMO——甲醇進料速率，g/h

Kd——失活速率常數

Kw——水吸附阻力因子

k1——甲醇轉化為C1的速率常數

k2——甲醇轉化為乙烯的速率常數

k3——甲醇轉化為丙烯的速率常數

k4——甲醇轉化為碳四的速率常數

k5——甲醇轉化為烷烴的速率常數

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

it——反應時間，s

W——催化劑裝填量，g

X1——C1的摩爾分率

XA——甲醇轉化摩爾分率

Xb——碳四的摩爾分率

Xe——乙烯的摩爾分率

Xp——丙烯的摩爾分率

Xr——烷烴的摩爾分率

Xw——反應過程中水的摩爾分率

Xwo——甲醇進料中水的摩爾分率

[1]吳秀章. 煤制低碳烯烴工藝與工程[M]. 北京：化學工業出版社. 2014：2-5.

[2]陳香生，劉昱，陳俊武. 煤基甲醇制烯烴（MTO）工藝生產低碳烯烴的工程技術及投資分析[J]. 煤化工，2005，5：6-11.

[3]李宏圖. 煤制低碳烯烴的技術路線及現狀分析[J]. 中國煤炭，2006，32（10）：47-49.

[4]張玉卓. 煤潔凈轉化工程（神華煤制燃料和合成材料技術探索與工程實踐）[M]. 北京：煤炭工業出版社， 2011.

[5]項東，彭麗娟，楊思宇，等. 石油與煤路線制烯烴過程技術評述[J]. 化工進展，2013，32（5）：959-920.

[6]吳秀章. 煤制低碳烯烴工業示范工程最新進展[J]. 化工進展，2014，33（4）：787-794.

[7]杜善明. 現代煤化工產業發展現狀及發展建議[J]. 科技和產業，2012，12（12）：41-44.

[8]黃曉昌，方奕文，喬曉輝. 甲醇制烯烴催化劑及其反應機理研究進展[J]. 工業催化，2008，16（1）：23-27.

[9]吳磊，劉子玉，夏林，等. SAPO-34晶粒形貌對甲醇轉化制低碳烯烴反應的影響[J]. 催化學報，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]. 化工學報，2010，61（7）：1674-1684.

[12]邢愛華，朱偉平，岳國，等. 甲醇制烯烴反應催化劑積炭問題研究進展[J]. 化工進展，2011，30（8）：1717-1725.

[13]南海明，文堯順，吳秀章，等. 甲醇制烯烴工業裝置催化劑積炭特性模擬研究[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]吳秀章. 基于工業裝置探究甲醇制烯烴反應機理[J]. 石化技術與應用，2013，31（5）：363-369.

[19]胡浩，葉麗萍，應衛勇，等. 甲醇制烯烴反應機理和動力學研究進展[J]. 工業催化，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，Froment 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，Froment 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，Froment 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，Fuglerud 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]齊國禎，馬濤，劉紅星，等. 甲醇制烯烴反應動力學[J]. 化工學報，2005，56（12）：2326-2331

[32]齊國禎. 甲醇制烯烴（MTO）反應過程研究[D]. 上海：華東理工大學，2006：31-51.

[33]胡浩，葉麗萍，應衛勇，等. SAPO-34分子篩催化劑上甲醇制烯烴反應的本征動力學[J].華東理工大學學報：自然科學版，2009，35（5）：655-660.

[34]胡浩. 甲醇制烯烴（MTO）催化反應過程的研究[D]. 上海：華東理工大學，2010：38-51.

[35]Alwahabi S M，Forment 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]吳秀章. 流態化技術在甲醇制烯烴（MTO）工業裝置中的應用[C]//第七屆全國流態化會議，北京，2013.

[39]鄭康，成有為，李希. 甲醇制烯烴流化床反應器的模擬與分析[J].高?；瘜W工程學報，2012，21（1）：69-76.

[40]鄭康，成有為，李希. 甲醇制烯烴流化床反應器的模擬與分析[D].杭州：浙江大學，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。

神華集團科技創新項目（SHJT-13-24）。

文堯順（1982—），男，博士，工程師，從事甲醇制烯烴領域技術研究。E-mail wenyaoshun@163.com。