第36卷第1期 2008年2月 燃料化学学报 Vo1.36 No.1 Journal of Fuel Chemistry and Technology Feb.2008 文章编号:0253-2409(2008)01-0099-05 吸收增强式甲烷水蒸气重整制氢循环反应模拟 李振山,蔡宁生 (清华大学热科学与动力工程教育部重点实验室热能工程系,北京100084) 摘要:吸收增强式甲烷水蒸气重整制氢反应可以生成高浓度的H 和较低浓度的CO、CO 。研究建立了考虑钙基吸收剂活 性下降对吸收增强式甲烷水蒸气重整制氢过程影响的多次循环反应模型,在实验数据验证的基础上,计算了三种吸收剂活性 下降特性对吸收增强式重整制氢过程的影响。结果表明,对于石灰石吸收剂,产生高纯H 的时间随循环次数的增加而急剧 下降;白云石循环反应活性提高,产生高纯H 的时间随循环次数的增加而缓慢下降;CaO/Ca。 AI。 O,,的循环使用次数明显大 于石灰石和白云石。 关键词:吸收增强;甲烷水蒸气重整;制氢;循环反应 中图分类号:TQ116.2 5 文献标识码:A Modeling of hydrogen production by sorption enhanced methane steam reforming reactions LI Zhen—shan.CAI Ning—sheng (Key Laboratoryfor Thermal Science and Power Engineering ofMinistry ofEducation,Department fTheromal Engineering, Tsinghua University,Beijing 100084,China) Abstract:High concen ̄ation of H2 with a low concen ̄ation of CO and CO2 Can be produced by the sorption・ enhanced methane steam reforming with the CO,SOrbent addition.A mathematic model for the sorption enhanced methane steam reforming considering the activi ̄loss of Ca..based sorbents was developed and veriied by fhe tex—. eripmental data.The effect of he acttivi ̄decay of three Ca—based sorentbs on the process of hydrogen produc— tion was explored by model calculation.The results indicate that the breakthrough time for high concen ̄ation of hydrogen production reduces rapidly wih tthe cycle numbers increasing for a limestone sorbent,while the break- through time decreases gradually for a dolomite sorbent.The effective cycle umber of CaO/Cal2All4O33 sorbent is far more than that of dolomite and limestone. Key words:sorption enhanced;methane steam reforming;hydrogen production;cyclic reaction 在传统的天然气水蒸气重整制氢反应过程 (Steam Methane Reforming:SMR)中加入CO2吸收 模型,用于预测反应过程的特性和规律。然而,针对 甲烷水蒸气重整制氢过程加入钙基吸收剂以同步吸 收反应过程所产生的CO 的模型却非常有限。Lee 等¨。。建立了以钙基吸收剂为基础的吸收增强式甲 烷水蒸气重整制氢固定床反应器动态模型,模型中 剂以改善系统的整体特性,被称作“吸收增强式重 整反应(sorption enhanced steam reforming reaction: SE.sMR)”。根据化学反应平衡原理,在CH 与 H o的反应过程中,如果产生的CO 能被就地吸 收,则传统的SMR反应热力学限制将被打破,重整 和水气变换反应将会向着H 增加的方向移动,CH 转化率将会提高¨“ 。 许多研究者对吸附增强式制氢过程进行了模 考虑了CaO与CO 反应过程分别受动力学或传质 过程所控制,并给出了相应的处理方法,经与实验结 果比较,证明模型较为合理。但Lee等在其所发表 的文章¨。。中只考虑了钙基吸收剂与CO 反应的单 次循环。 拟,这些模拟工作主要集中在以K CO,一水滑石为 吸附剂的制氢过程。根据水滑石在高温下与CO 吸附平衡特性和吸附过程动力学特性 J,Ding 等 以及Xiu等 .9 建立了考虑多组元气体质量平 衡、压降方程和能量方程等固定床反应器一维动态 由于钙基吸收剂的反应活性不可避免的要随循 环次数的增加而下降,从而导致吸收增强作用的逐 渐减弱。这一点对于采用固定床作为重整反应器的 制氢过程非常重要。无论是SE—SMR制氢相关方 面的研究,还是反应器的设计与操作,吸收剂循环反 收稿日期:2007-05—15;修回日期:2007-09-23。 基金项目:国家重点基础研究发展规划(973计划,2006CB705807);国家高技术研究发展计划(863计划,2003AA501330)。 作者简介:李振山(1975一),男,辽宁康平人,博士,E-mail:lizs@mail.tsinghua.odu.cn。 维普资讯 http://www.cqvip.com l0o 燃料化学学报 第36卷 应特性对制氢过程的影响都必须予以考虑。本研究 在前人有关实验和理论分析的基础上,建立了考虑 钙基吸收剂活性下降对SE—SMR制氢过程影响的 多次循环反应模型,利用数学模型采用数值模拟方 法研究吸收剂循环反应活性变化对SE—SMR制氢 过程的影响。一方面,可以根据有限次制氢循环的 实验数据对后续循环反应工况进行预测。节省财力 物力;另外一方面,根据模型,可以对反应器的设计 与操作进行指导。 l SE—SMR制氢多次循环反应模型 1.1 CaO循环吸收CO2模型研究表明,CaO与 CO 的反应活性会随循环使用次数的增加而快速降 低¨ 。Wang等¨ 给出了描述CaO的活性降低规 律的公式,该公式只用一个参数就能确定CaO的活 性下降趋势,在循环次数较低时(≤100),实验数值 与模型计算值吻合较好,而当循环次数较大(≥ 500)时产生误差。针对此问题,Grasa等 3J对石灰 石和白云石进行了多达500次循环的反应/煅烧实 验,发现随循环次数多达几百次后,CaO的活性不 再下降,维持为一恒定值,并总结出如下的经验公 式: 1= u. + r) (1) 1一 …一 为每次循环的最终转化率,k为吸收剂活性下 降系数,Ⅳ为循环使用次数, 为吸收剂最低转化 率。本研究采用式(1)来描述CaO活性下降规律。 作为比较,采用了三种CO 吸收剂,石灰石、白云石 和CaO/Ca A1 O33。三种吸收剂活性下降特性及 实验数据来源见表1。 表1 吸收剂活性下降系数及数据来源 Table 1 Reactivity decay coefficients of sorbents andtheir references 通过对吸收剂进行不同温度、不同CO 浓度以 及不同压力下的实验,总结出用表观动力学表达式 (2)描述CaO与CO 反应过程第Ⅳ次循环转化率 随时间t的变化: 訾-8.124×(1一 ( 一 , ) (2) P为总压,cc02为反应气体中CO2浓度,ccq.c02是CO2 与CaO反应的平衡浓度,随温度变化的表达式为: o2=4.137×10 2exp(-20474/T)( )(3) 第Ⅳ次循环过程中单位质量的CaO对CO 的 吸收速率可以表达为: =忐 (4) 在t时刻,时间增加 后,CaO的转化率可以表述 为下式: r N,t+At N.t+Mca0 J rN,cbndt (5) 1.2化学反应动力学模型xu等¨ 建立的甲烷 重整反应动力学模型比较成熟,被多数的研究者所 使用。利用如下三个化学反应式描述甲烷水蒸气重 整反应(SMR)和水气变换反应(WGS)。 甲烷水蒸气重整反应(SMR): CH4+H2O=CO+3H2 (6) CH4+2H2O=CO2+4H2 (7) 水气变换反应(WGS): CO+H2O=CO2+H2 (8) SMR、WGS三个反应的化学反应速率可表示 为: (PcH4PH20一 )/(DEN) (9) k2(pcrt,P2 一 )/(删 l0) (PcoPH20一 )/(DEN) (11) DEN=(1+Kc0pc0+Kx2p.2+Kcn4Pcn4+ 2oP.2o/ P ) (12) 其中,尺 (i=1、2、3)为化学反应速率(kmol・kg— cat~・S ),k 为反应速率常数, (i=1、2、3)为化学 平衡常数, (j=CO、CO 、H 、CH 、H 0)为气固 相催化反应吸附常数,DEN为吸附项,表示各组分 对催化剂活性中心竞争产生的影响,上述参数的具 体计算公式可参见文献[9]。 由(6)一(11)可得重整反应中各种气体生成和 消耗的速率: =一(尺l+尺2) (13) rH 0=一(尺l+2尺2+尺3) (14) 维普资讯 http://www.cqvip.com 第1期 rH,=李振山等:吸收增强式甲烷水蒸气重整制氢循环反应模拟 l0l 3Rl+4R2+R3 =R2+R3 (15) (16) (17) 与入口处进料组分浓度和温度一致,反应总压与进 料的总压一致。 I'CO=Rl—R3 I'CO,1.5求解方法方程(18)、(19)经离散后,变成一 1.3 吸收增强式重整反应模型 吸收增强式重整 组带初值的非线性微分方程租,同时该方程组又是 刚性的,因此,采用Matlab中的“odel5s.m”函数来 反应模型在管状固定床催化反应器拟均相模型基础 上,耦合了CaO吸收CO 模型。模型作出如下假 设:反应气体为理想气体;固定床内为轴向活塞流; 径向上没有温度和浓度梯度分布;催化剂与吸收剂 求解上述方程组。计算中所用到的参数及其数值见 表2。 表2计算中所用到的主要参数及其数值 粒径一致;反应器管壁温保持不变。 质量平衡方程: aCi a(HCi) ,, 、 s_ 一—— +.,7c atl(1一s)pcat i— r/c.o(1一s)Pc.orN-cb .(18) i≠C0 ,rN.cb -0;i=C0 ,rN.cb = 其中Ci(i=CO、CO 、H 、CH4、H 0)为各组分气体 的浓度(kmol/m ),t为时间(s ),z为反应器轴向 距离(m),ri(i_CO、CO 、H 、CH4、H 0)为各组分 气体生成和消耗的速率(kmol・kg—cat~・s ),H为气 体流速(m/s),s为固定床空隙率,叩。 为催化剂效 率因子,叩 帕为吸收剂效率因子,P。 和P 帕分别为催 化剂和氧化钙密度(kg/m )。 能量平衡方程: [(1-8)p。Cps+即gCpg]等= 一 -(1 + 麦一.,7ca0(1_s)pcaOFN,cbn +hw(Tw— )麦(19) 其中P。为反应固体密度(kg/m ),P 为气体密度 (kg/m ),CPS为固体比热(1d/(kmol・K)),CPg为气 体比热(1d/(kmol・K)), 为温度(K), 为壁温 (K),HR:为(11)一(13)的反应热(1d/kmo1),Hcb 是 CaO吸收CO 的反应热(1d/kmo1),h 为与壁面的 换热系数(kJ・m /(s・K)),D 为反应器内径(m)。 反应器内的压降采用Ergun方程: 一=一—————一l——— ————一d :一吐[ +1+1..z dp s dpogH … 75]×10一s(20)× I、一 ZI J l 1.4初始条件与边界条件 Ci=C¨nit,T= it当t=0时, C =C¨,T= P=P 当z=0时 (i_CO、CO2、H2、CH4、H2O) 反应开始t=0时,各组分浓度和温度等于给定 的初始值;在反应器入口z=0处,各组分浓度、温度 Tab1e 2 Parameters and their values in simulation Parameter Value Size 200 m一425 m Temperature 650℃ Pressure 1.5№ Reactor length l2.7 X 10I2 in Reactor diameter 1.9 X 10一 in Sorbent mass 2.6 X10一 kg Catalyst mass 1.4 X10一 kg H2O/CH4 4 CH flow 60 mL/min N,flow 200 mL/min 2模拟结果分析 2.1 以白云石为吸收剂循环制氢反应特性 为了 研究循环次数对SE—SMR制氢过程的影响,分别模 拟计算了温和煅烧条件下(850 oc煅烧)第2、15、25、 50、100、500、1 000和6 000循环次数下反应器出口 各气体含量随时间的变化,见图1。为了和模拟结 果进行比较,第2、15、25次循环的实验数据 也同 时显示在图中。从图1可以看出,第2、15、25次循 环中H 含量随时间的变化特性实验数据基本与模 拟结果吻合,从而验证了所建立的多次循环反应模 型的正确性。 图1不同循环次数下产品气中各组分体积分数 随时间的变化 Figure 1 Evolution of product gas composition with time at different cycles (regeneration at 850℃in N2 for dolomite sorbent,com ̄ison ofthe siInulali0n m panial experiment in 25.cycle test ) ,0维普资讯 http://www.cqvip.com 燃料化学学报 第36卷 CaO的循环反应活性对SE—SMR制氢过程影 响极其重要。由于高温煅烧导致CaO颗粒内部微 观结构发生变化,从而使得CaO吸收能力下降。 CaO与CO 反应生成CaCO,后,必须进行煅烧分 解以使CaO再生,以便重复利用。在CaCO 的高 温煅烧过程中,新产生的CaO具有较大的比表面积 和孔隙率,但CaO颗粒内部的烧结也同时发生, CaO晶粒发生聚合,致使比表面积迅速减小,从而 导致CaO的反应活性明显下降。从图1可以看出, 当白云石的循环反应次数超过100次后,由于其活 性的大幅下降,吸收增强阶段所持续的时间下降较 大。随着循环次数的继续增加,由于钙基吸收剂最 后能维持一个恒定的最终转化率¨ ,所以吸收增强 阶段最终能保持在20 min以内,而基本不再受循环 次数的影响。 图1的结果是在温和煅烧下(850℃煅烧)获得 的,煅烧过程中通人N 气体以稀释CO 。但这样 将会使煅烧所释放出的CO 稀释,从而达不到CO 分离富集的效果。为了在煅烧过程中获得高纯度 CO ,必须升高煅烧温度,因为CaCO,的分解受CO 平衡浓度所限制(见式(3))。当煅烧温度提高时, 钙基吸收剂的活性下降得更快¨n驯。为了研究利 用吸收剂制高纯氢的同时,还达到CO 分离减排的 效果,本研究对高温煅烧下吸收剂活性下降对制氢 过程的影响进行了计算。高温下吸收剂活性下降特 性见表1。图2为以白云石为吸收剂的高温煅烧下 (1 000℃煅烧)不同循环次数下反应器出口各气体 含量随时间的变化。 Q 垄 图2 白云石高温煅烧不同循环次数下产品气中 各组分随时间的变化 Figure 2 Evolution of product gas composition with itme at different cycles at high calcinationtemperature wim dolomite as sorbent 从图2可以看出,吸收增强阶段所持续的时间 随着循环次数的增加而快速下降,吸收增强阶段最 终所能保持在10 min左右。 2.2以CaO/Ca。2AI。4O33为吸收剂循环制氢反应 特性以CaO/Ca 2A1 4O33为CO2吸收剂高温煅烧 下(980 oC煅烧)的SE—SMR制氢多次循环模拟结果 见图3。整个制氢过程仍然可以分为吸收增强阶 段、过渡阶段和传统重整反应阶段,但由于CaO/ ca Al o,,吸收剂反应活性随循环使用次数的增加 下降较缓慢,所以,在较大的循环使用次数下,吸收 增强阶段仍能保持较长的时间。 图3 以CaO/Cal2AI。 O 为吸收剂高温煅烧时 不同循环次数下H:体积分数随时间的变化 Figure 3 Evolution of H2 fraction with time at diferent cycles at high calcination temperature with CaO/Cal2AIl40 as sorbent 2.3以石灰石为吸收剂循环制氢反应特性石灰 石也可以作为CO 吸收剂,并且受到了许多研究者 的关注。与白云石和CaO/Ca 2A1 4O33相比,虽然石 灰石中含有较多的CaO,但其与CO 的反应活性随 循环使用次数的增加下降非常迅速。石灰石也可以 用于SE—SMR重整制氢过程,为了揭示不同吸收剂 对制氢过程的影响,本研究也对以石灰石作为SE— SMR重整制氢中的CO 吸收剂进行了多次循环模 拟计算,结果见图4。以石灰石为CO 吸收剂的 SE—SMR制氢第一次循环基本与白云石和CaO/ Ca A1 O,,吸收剂的反应特性类似,吸收增强段持 续了很长时间。但由于其活性的快速下降¨引,第 15次循环反应时吸收增强段的持续时间已经明显 减小。从第50次循环开始,吸收增强段维持在 5 min左右,此后随着循环次数的增加而变化较小, 见图4。 通过对上述SE—SMR制氢多次循环反应过程 的分析可知,经过多次循环反应后,由于CaO吸收 剂活性的快速下降,吸收增强段所持续的时间随循 环次数的增加都会连续减少。对于石灰石吸收剂而 言,产生高纯H 的时间随循环次数的增加而急剧 维普资讯 http://www.cqvip.com 第1期 李振山等:吸收增强式甲烷水蒸气重整制氢循环反应模拟 103 下降,这主要是因为石灰石的活性随循环次数的增 ,o/§ LIo 3B出 }1 石相比,白云石的循环反应活性有所提高。在三种 吸收剂中,CaO/Ca。2A1。 O33吸收剂的循环反应活性 最高。 3结语 建立了SE・SMR甲烷水蒸气重整制氢的多次 循环模型,利用文献中的实验数据对多次循环模型 进行了验证,两者比较吻合,从而证明了模型的正 确性。 吸收剂的反应活性随循环次数的增加而下降, 从而导致吸收增强作用的逐渐减弱。但随着循环次 图4以石灰石为吸收剂时高温煅烧不同循环次数下 数的继续增加,由于钙基吸剂最后能维持一个恒定 H 的体积分数随时间的变化 的最终转化率,所以吸收增强阶段最终能保持在一 Figure 4 Evolution of H2 fraction、Ⅳim time at different cycles at hJIgh calcination temperature wim limestone as sorbent 个稳定的时间段内,而基本不再受循环次数的影响。 不同循环反应活性的吸收剂对SE.SMR的影 加而快速下降的原因。同时,也可以看出,在SE・ 响也不同,CaO/Ca。2AI。 O33的循环使用次数明显大 SMR制氢过程中,产生高纯氢的时间主要取决于吸 于石灰石和白云石。 收剂对CO 的循环吸收能力,吸收剂循环吸收能力 合成或制备具有高循环反应活性和循环稳定性 增加,则对应的吸收增强段的时间也增加。与石灰 的吸收剂是SE.SMR制氢反应的核心问题。 参考文献: [1]李振山,蔡宁生,黄煜煜.吸收增强式甲烷水蒸气重整制氢实验研究[J].燃料化学学报,2007,35(1):79-84. (LI Zhen—shah,CAI Ning—sheng,HUANG Yu—yu.Experimental research on hydrogen production from sorption enhanced steam methane erfor- ming[J].Journal of Fuel Chemistry and Technology,2007,35(1):79-84.) [2] BALAsuBRAMANIAN B,ORTIZ A L,KAYTAKOGLU S,HARRISON D P.Hydrogen from methane in a single—step process[J].Chem Eng Sci,1999,54(15—16):3543—3552. [3]ORTIZ A L,HARRISON D P.Hydrogen production using sorption—enhanced reaction[J].Ind Eng Chem Res,2001,40(23):5102-5109. [4]LI z—S,CAIN—S,YANG J-B.Continuous production ofhydrogenfromsorption-enhanc ̄l steammethanereformingintwo parallelfixed—bed reactors operated in a cyclic manner[J].Ind Eng Chem Res,2006,45(26):8788-8793. [5]CARVILL B T,HuFroN J R,ANAND M.Sorption enhanced reaction process[J].AIChE J,1996,42(10):2765-2772. [6]DING Y,ALPAY E.Equilibria and kinetics ofCO2 adsorption on hydrotalcite adsorbent[J].‘Chem Eng Sci,2000,55(17):3461-3474. [7]DING Y,ALPAY E.Adsorption—enhanced steam—methane reforming[J].Chem Eng Sci,2000,55(18):3929-3940. [8]XIU G H,SOARES J L,LI P.Simulation offive-step one-bed sorption-enhanced reaction process[J].AIChE J,2002,48(12):2817.2832. [9]XIU G H。LI P,RODRIGUES A E.New generalized strategy forimproving sorption-enhanced reaction process[J].Chem Eng Sc。2003,58 (15):3425-3437. [1O]LEEDK,BAEKBIH,YOONWL.Modelingand simulationforthemethane steamreforming enhanced byin sierCO2removal utilizingthe CaO carbonation for H2 production[J].Chem Eng Sci,2004,59(4):931-942. [11]ABANADES J C,ALVAREZ D.Conversion fimits in the reaction of CO2 wiht lime[J].Energy Fuels,2003,17(2):308-315. [12]WANG J。ANTHONY E J.On the decay behavior ofthe CO2 absorption capacity ofCaO-based osrbents[J].Ind Eng Chem Res,2005。44 (3):627-629. [13]GRASAG S,ABANADES JC.CO2 capture capaciytofCaOinlong series of carbonation/calcination cycles[J].IndEngChemRes。2006, 45(26):8846-8851. [14]LI z S,CAI N S,HUANG Y Y,HAN H J.Synthesis,experimentla studies,and analysis of a new calcium.based c ̄bon dioxide sorbent [J].Energy Fuels,2005,19(4):1447—1452. [15]LI z S,CAI N S,HUANG Y Y.Effect of preparation temperature on cyclic CO2 capture and mulitple carbonation-calcination cycles for a new ca-based CO2 sorbent[J].Ind Eng Chem Res,2006,45(6):1911-1917. [16]XU J,FROMENTG F.Methane steam erforming,methanation and water-gas ̄ifl:L Intrinsic iknetics[J].AIChE J,1989,35(1):88-96.
因篇幅问题不能全部显示,请点此查看更多更全内容