煤炭气化气流床气化炉的数学模拟

第24卷第4期煤炭转化2001年10月COAL CONVERSIONOct.2001煤炭气化气流床气化炉的数学模拟步学朋)彭万旺1)徐振刚2)摘要简要介绍了煤炭气流床气化的原理,总结了到目前为止煤炭气化气流床气化炉数学模拟情况,包括简单平衡模型和动力茡模型(一维或多维),给岀了这些数学模型模拟的主要內容(对气化过程流体力学、热力学、化学反应和质量、能量及动量平衡考虑情况)和模型的主要结论,以及典型气流床气化炉的模拟煤气组成和煤炭转化率数值与实验值或实际操作值旳比较情况,结果显示主要组分模拟误差较小关键词煤炭气化,气流床气化,数学模拟中图分类号TQ529造,而我国第一座IGCC示范电厂也已立项,其气化0引言岛也将采用气流床气化工艺.因此应重视对气流床气化工艺的研究及开发工作煤炭气化是煤洁净利用的关键技术之一,它以对气化炉气化过程数学模拟的研究不仅有利于提高碳转化率、冷煤气效率,降低气化过程旳氧耗及深入了解气化过程规律,而且可以用于指导气化炉煤耗为目标,并向加压及液态排渣、大型化等方向发设计及生产过程优化控制.对气流床气化炉的数学展,从而达到改善环境及降低产品成本的目的.加压模拟主要始于20世纪70年代,国内外对此进行了气流床气化技术是国内外优先发展方向之,气流床大量研究,取得了许多成果气化技术有如下特点:①能够气化任何变质程度的煤及煤的加氫残渣、石油焦等;②气化强度很高,碳1气流床气化原理转化率高;③产品煤气中不含焦油和酚类物质,环境友好;①其缺点是氧耗高、需设置磨粉、显热回收1.1气流床气化过程描述及除尘等较庞大的辅助装置.已工业化或正在示范从气流床气化工艺分析看,煤粉(或水煤浆)与的加压气流床气化技术包括湿法加料的 Texaco,气化剂(O2/HO)经喷嘴喷入气化炉的燃烧区,由于Detec和干法加料的 Shell,, Prenflo及GSP等,并气化炉该处温度高达1500C~2000C,因此煤粉在世界范围内得到广泛应用受热升温速度很快(>104C/s),可认为煤粉中的我国20世纪⑦0年代先后进行过熔渣炉气化及残余水分快速(瞬间)蒸发,同时由于热分解反应速K-T气化研究,也建立过工业化装置,后因耐火材度大大高于煤粉的燃烧及气化反应速度,所以细小料等问题而停止.20世纪80年代后在鲁南化肥厂、的煤粉颗粒开始发生快速热分解,即脱挥发分,生成上海焦化有限公司、渭河化肥厂、淮南化肥厂等引进半焦和气体产物.在富含氧气及高温条件下,挥发分了 Texaco气化炉用于生产合成氨及甲醇,通过多中的活性可燃成分如CO,H2,CH1及焦油与O2发年的实践积累了丰富的经验,已达到了长周期、高负生气相燃烧下应生成CO和HO,并放出大量热荷、安全和稳定运行的要求,且在喷嘴、耐火砖及水量中国煤化工应的进行.由于气相燃煤浆泵等设备国产化方面取得了突破性进展目前,烧CNMHG认为在氧气存在的情况利用干法加料的 Shell气化工艺正用于某化肥厂改下,上述气相燃烧反应达到完全,亦即在氧气存在国家“九五”科技攻关计划项目(97-A26-02-02-05)1)高级工程师;2)研究员,煤炭科学研究总院北京煤化学研究所,100013北京收稿日期:2001-01-16;修回日期:2001-05-20煤炭转化2001年时,气相中不含CO,H2,CH4及焦油.煤中的挥发分反应平衡时的煤气组成及平衡温度.项友谦用能析岀后,发生半焦燃烧及气化反应,与水蒸气及CO2量最小原理,建立了加压气化平衡模型,并用4种方反应,如此时仍有氧气存在,则在气相中仍发生CO法对微分方程求解模拟结果显示平衡模型对气流和H燃烧反应.气化炉中的氧气反应完后,半焦与床的模拟效果要好于固定床. Watkinson提出平水蒸气、CO2和H2等继续发生气化反应,同时气相衡模型,通过质量和能量平衡及反应平衡方程式关中还有水煤气变换反应和甲烷裂解反应等对气流联,可以得到产品煤气组成、产率和最佳适宜温度床气化,一般将变换反应和甲烷化反应视为平衡反并对9种气化炉型的工业化和半工业化气化数据进行了比较.对产品煤气中的CO和H2含量误差在士从如上分析可以认为在对气流床模拟时,考虑%之内,H2S和COS浓度可以准确的预测,但气化炉中气固相的组成分别为气相:O,H1O,CO,CO2预测值的准确性要差一些.结果还表明该模型CO,H2, CH,, N2(E NH3 ), H,S(+COS), Tar对气流床气化炉模拟效果最好,流化床次之,而固定相:C,H,O,N,S,Ash, Moisture.另外考虑气固相由于一些不确定因素如挥发分含量组成等的存温度ˆ,T∷数学模拟的目的就是求出气化炉内每在,模拟误差较大 Ruprecht等建立了平衡模型质点处的上述参数的值并用于评价试验数据分析.平衡模型可用于对出囗组成及温度进行简单预测,由于其假定的条件较理1.2气化反应动力学想,如假定所有反应都达到平衡,在实际过程是不可煤粉热解反应速度很快,主要根据经验式对热能的,因此其用途受到一定限制,对过程控制及气化解速度和产物进行估算,大致有wen等给出的炉设计参考价值较少几种表达形式.气固非均相反应,由于气化炉内反应2.2考虑气化过程动力学的模型温度很高,最高温度可达2000C以上,所以半焦与为了更准确的用数学模型对气化及燃烧过程进O2,H2O,CO2的反应速度由气膜扩散及灰层扩散控行模拟,需考虑气化过程三传一反及动力学行为.许制,而半焦与H的反应较慢,它由化学反应控制.多研究者对此进行了研究,建立了一维、二维及多维从非均相反应动力学模型看,分为未反应核收缩模模型,Fild等对1967年以前的工作进行了归纳,模型及灰层剥离模型,文献[5认为灰层在反应过程中型属活塞流 Ubhayakar等忽略了表面反应,但考剥离,因此反应受气膜扩散和化学反应动力学控制虑了轴向混合、挥发分的气相反应及热解反应而文献[1,6认为,由于气流床气化炉中煤粉颗粒所 Smith等先后建立了一维及二维煤粉气流床燃占体积很小(<1%),且停留时间只有2s~3s,因烧及气化模型,模型中考虑了煤粉颗粒大小分布对此煤粉颗粒碰撞几率较小,在煤粉颗粒表面因反应反应速度的影响.Wen等给出一维模型来描述形成的灰层可认为仍保留而不剥离,所以反应速度 naco气化炉中的混合情况,每一微元被处理成气可由未反应核收缩模型得到反应速度其它学者也相是完全混合固相以活塞流通过整个反应器,固相研究了气化反应动力学问题,基本一致的观点是要粒子的速度用 Stokes方程近似描述模拟了中试装同时考虑化学反应控制及扩散控制.使用比较多的置使用煤及加氫残渣气化情况操作参数的变化对如对加氬气化过程模拟时,采用了 Johnosn的加煤气组成的影响等. Govind等对上述模型进行了氢热解动力学模型;wen等提出的未反应核收缩修正,加入了动量平衡.考查了煤气组成与进料速模型也使用较多, Smith等也采用了独特的动力度、氧/煤比及汽/煤比的关系. Sprouse等模拟学表达式.加氢气化反应,前者使用了经验关联式,但不适合次2气流床气化数学模型综述烟r中国煤化工,但因其复杂性导致其应CNMHG了4种煤沿气化炉的轴向冋气柜组成分巾,订论了氧碳比及蒸汽量等2.1平衡模型对模拟结果的影响. Vamyuka建立了一维稳态模型,发现气固相最大温度位于氧气耗尽处,考察了关平衡模型认为气化过程所有化学反应达到平键参数对操作的影响.我国科技部“九五”期间也安衡,通过对气化过程质量及能量平衡方程求解得到排了对流化床、干法及湿法加料气流床气化炉进行第4期步学朋等煤炭气化气流床气化炉的数学模拟数学模拟的攻关课题,取得了丰硕成果.现将最近的仼何模型模拟的目的是能反映气化炉内物料研究成果及文献[6,12]对以前的煤粉气流床气化及动及化学反应情况,并获得出口处煤气组成及温度、燃烧模拟情况的总结一并列于表1中碳转化率、气化效率等参数.气流床气化炉数学模拟与实际(或实验)结果的对比情况见表2,可见对几3模拟煤气组成与实际(或实验)结果种不同炉型煤气中的主要组分,模拟误差小,说明模对比拟基本上是成功的.表1气流床气化数学模拟汇总Table 1 Reviews of models for coal entrained-bed gasifiersAuther(s) ReferencesMain Contents/ DescriptieMajor ResultsMehta1976 Phenomenological model combination of plug-flow Describes the outlet conditions for aand perfectly-stirred calculations, gas phase chemi- three-stage entrained flow gasifieral equilibrium, overall char gasification with Arrh1979quilibrium chemistry; one-step devolatilization and around the particles of Rock well iOne-dimensional hydrogasification; free-stream ehydrogasification modelFHP gasifier, compared well with outmeasurements for coalFinson et al1978 Steady, one-dimensional model. with gaseous ki- Predictions compared with laboratornetics: kinetics from data for coal pyrolysis het- gasification measurements of axial tem-erogeneous oxidation: considered the char struc- perature and gas compositionUbhayakar1977 One-dimensional steady-state; two-step de- Applied to gasification combustion andStickler & Cannonbined diffusion & surface reac- hydropyrolysis with good agreementBlake et al1977-1979 Three-dimensional, transient, mixed finite-ele- Study the effect of particle sizement, finite-difference scheme with separate Eule- conservation of mass, momenturm andparticle equations; two equation gaseous tur- energy, studied the loctbulence: separate turbulent kinetics energy for par- vectorsticlesBarnhart et alCombination of plug-flow and perfectly-stirred ele- Predicts carbon burnout product gasments; gas phase in equilibrium with two-equation composition, temperature; reportediment&. theory for cyclone gasifier er.kinetics for pyrolysis, char oxidation and CO oxi- fair to good comparison between exper1979-1981 One-dimensional steady-state, two-step de- Good agreement with one-dimensionaltion rates, one-dimensional zonal radiation, equi- gasifier measurementslibrium gas phase reactionsSmith. Fleche*,12 1979-1981 Two-dimensional, axisymmetric, with recircula- Gas-phase components reported andion, k-E turbulence model, Lagrangian partdevolatilization, diffusicerogeneous rates, four-flux radiation with scatterCl1978fier into three zones: one-step de- Predicted temperature and gas compogas phase combustion completely or sition data for Texaco gasifierequilibrium gas phase reactions, consider diffusion pared with pilot plant gasified coal andand kinetic rates for char reactionH-coal residue, with good agreement.1980 One-dimensional steady-state model, Johnsons Modeling bench-scale hydrogasifier ofhydropyrolysis model, continuity equations for Citys Service R& D Coboth phases, mixture momentum balance and ther中国煤化工Govind and shan1984efined the Wen & Chuang mo-composition and carbon conmomentum balanceTHCN MH Gpends on three essential pameters: the fuel rate, the oxygen tofuel ratio and the steam-fuel ratio: theoptimum steam-fuel ratio is betweension: the steam-fuel ratio significantlyaffects the gas product composition.煤炭转化续表1Auther(s)eferences Yeain Contents/DescriptionMajor ResultsGong Suning et al1987 Image the gasifier being composed of several paral- Axial gas/solid temperature and conlel reactors of which each has respective homoge- version profiles were got, the initial reneous particle size and also its own material and action rate very fast, and then slowheat balances; one-dimensional equations similar to down significantly the wall tempera-Smith's: the finite element method was used to cal- ture has influence on reaction. andculate the wall temperature distribution and heat should be taken as an important param5 1988 Refined the Smith two-dimensional model, gave Modeling the radial and axial gas phaseof full/ partheterogection rates and heatloss on prediction results; the effect ofO/C and steam on gas compoRuprecht et al1988 Considered only theaction, the gas composi- The agreement between prediction andtion include CO, H2. COz, no CHa: equilibrium gas measured data is very good, the dis-phase reactions, carbon conversion controlled by crepancy normally amounts to 1% orO/C: mass and energy balance.less: it has been used for the evaluationka et al1995 One-dimensional steady-state, plug-flow, assh didhe boundary between combustionhe reacting particle surfadistinctbined diffusion and chemical reaction kinetics. usedepend on the gasTGa data, equilibrium gas phase reactions, the sure: peak temperature of gas and solidmomentum balance neglectednear the boundary coal conversion isnot sensitive to the pyrolysis rate, bisensitive to pressure: the effect of opNi & willia17 1995 A multivariable model was set up on the basis of e- The oxygen to coal ratio is the mostquilibrium, mass balance and energy balance, the important control variable for theShell gasifieral model the gas composieffhasIngture of inlet flow, steam to coal ratiohas influence on the efficiency.Liu g s et al1999 One-dimensional steady-state, plug-flow. the gas The predicted carbon conversion agreeturbulence is neglected, kinetic data from Ptga with those measured in the gasifier of 8ests: considered the pore structure and reactivity coals; the better prediction results canof char: the heat balance between the wall, gas, be got by considered char structure;modeled the CSIRO atmospheBu XuepengOne-dimensional steady-state, divided gasifier into The 48 t/d and 2 600 t/d Prenflo gaPeng Wanwangthree zones: the devolatilization composition by fiers were modeled, axial gas-solidXu Zhengangcalculated, gas phase combustion completely, oth- temperature and composition were goter gas reactions equilibrium: considered diffusion main gas composition agree with reaand kinetic rates for gas-char reactionsdata; the effect of oxygen to coal ratiogastrication results wasYu Zunhong2000 Based on cold model experiment, divided gasifier Based on Texaco gasifier realistic opernto three zones, the reactions different in each ation data, predicted the feature tem-ong Xin et alne the mixed time and reaction time are consid- perature and gas composition at threeered as main variables, chemical reactions equilibri- different load, they are agreement withLi zhenDivided gasifier into lots of zones along axial diPredicted axial gas-solid temperaturetion, in order to consider thrV凵中国煤化工sition of Texaco gasifier;Wang Tianjiao et altructure and size; gas phasetion kinetics: retention timeCNMHGn, coal kinds on operationGao Juzhong* 2000 Equilibrium model, only thermodynamics was con- Predicted the gas composition and gasiwang ningboidered took the relation between carbon conver- fication indexes of two coals gasified inZhang yaruion and O/C ratio, reaction extent as boundary Shell gasifiers in pilot plant and demonstration plant respectivelyFrom reference 12,** See95"item reports第4期步学朋等煤炭气化气流床气化炉的数学模拟表2模拟结果与实际结果的对比Table 2 Comparison of computational results from the models with experimental/ operation resultsAuthor(s)/Gas composition/%I)CO CO2 CH. H,S+COS N2+Ar H20 Carbon conversion/%37.3047.8614,450.050.26atkinsonxperiment24.3047.1013.200.090.4012.7Model7047.2013,300,332.210.3013,0WatkinsonExperiment30.6061.501.600.001.30ShellModel30.3061,501.400.085.38WatkinsonExperiment34.6055.407.000.001.951.01K-T34,9055.406.700.001.931.03Qizhi NiL17.8065.401.700.000.940.001.9029.7165,651.360.01Yu Zunhong7.4745.8416.0.1038.3144,8415.010.110.822)Li Zheng29.8041.0511.100.8Texcao30.5540.8010.270.13Gao JuzhongOperationShell63,790.800.0357.981.19Bu XuepengPrenfloModel25.9360.212.660.020.7710.41Note: 1) Percent of volume: 2) Include NH3; 3) Others 0. 21础上,建立二维及多维动态模型,并充分考虑颗粒大4结束语小分布、气化炉内气固相流体力学行为等.另外为适应气化炉使用中国煤种的操作,需要研究适宜中国囯外对气流床气化炉的数学模拟研究起步较煤粉的髙温髙压快速热解反应动力学以确定其热解早,经过最近几年的努力,我国对气流床气化炉的数速度及热解煤气组成,研究适宜煤种的气化反应动学模拟也取得了初步成功,但正如专家们指岀的那力学,为数学模拟提供基础数据.最后对建立的数学样,目前离真正应用还有一定距离.根据在模型建立模型还要到生产实际中进行检验,对某些参数进行及求解过程的经验,笔者认为下一步应继续进行气修正,以求能较准确反应气化炉内气化行为并用于流床气化炉的数学模拟工作.在目前一维稳态的基指导气化炉设计及生产操作[] Wen C Y, Chuang T Z. Entrained-bed Coal Gasification Modeling Interim Report. US Department of Er2274-T1,1979.1-41[2] Anthony D B. Howard J B Coal Devolatilization and Hydrogasification, AIChE J. 1976.22(4):625-656[3 Ubhayakar S K, Stichler D B, Ganon R E Modelling of Entrained Bed Pulverized Coal Gasifiers. Fuel, 1977.56(2):281-286[4] Fu Weibiao, Zhang Yanping, Han Hongqiao et al. A General Model of Pulverized Coal Devolatilization. Fuel, 1989, 68(5)505-509[53 Brown B W, Smoot L D, Smith P J et al. Measurement and Prediction of Entrained-flow Gasification Processes. AIChE J1988,34(3):435-446中国煤化工[6 Govind R, Shah J. Modeling and Simulation of an Entrained FlowCNMH30(1):79-92[7 Johnson J L. Kinetics of Coal Gasification. New York: John Willy8 Smith P J. Smoot L D Prediction and Measurement of Optimum Operating Conditions for Entrained Coal Gasification Pro-esses.DE82015610,1981.1-40[9]项友谦.煤气化过程热力学平衡组成的理论计算.煤气与热力,1986(1):4-9[10 Watkinson A P. Lucas J P. Lim C J. A Prediction of Performance of Commercial Coal Gasifiers. Fuel. 1991.70(4):519煤炭转化2001年[11 Ruprecht P, Schafer W, Wallace P. A Computer Model of Entrained Coal Gasification Fuel,1988,67(6):739-743[12 Smith P J. Smoot L, D. Prediction and Measurement of Optimum Operating Conditions for Entrained Coal Gasification Pro-cesses,DE82015611,1981.1-30[13] Spouse K M. Modeling Pulverized Coal Conversion in Erd flows. AIChE J,1980,126(6):964-975[14] Goyal A. Mathematical Modeling of Entrained Flow Coal Gasification Reactors: [Ph D Thesis]. Chicago: lllinool,1980.1-10[15 Vamvuka D, Woodburn E T, Senior P R Modeling of An Entrained Flow Coal Gasifier. Fuel,1995,74(10):1452-1465u Zhihong. Improved Simulation System of Pulverized Coal Gasifier. J of Chemical Industry and Engine[17] Ni Qizi, Williams A, A Simulation Study on the Performance of An Entrained-flow Coal Gasifier. Fuel,1995,74(1):102[18] Liu G S, Rezaei H R, Tate A G et al. Modeling of a Pressurised Entrained Flow Coal Gasifier. In: Prospects for Coal Sciencein the 21st Century. Taiyuan: Shanxi Science & Technology Press,1999.579-582MATHEMATICAL MODELING OF COAL ENTRAINEDBED GASIFIERSBu Xuepeng Peng Wanwang and Xu Zhengang( Beijing Research Institute of Coal Chemistry, China Coal Research Institute, 100013 BeijingABSTRACt The coal entrained-bed gasification principle was briefly introduced in this paper. Based on relative documents, the mathematical modeling of coal entrained-bed gasifiers werereviewed. These include simple equilibrium models and kinetic models (one-dimensional or two/three-dimensional). The main contents-about hydromechanics, thermodynamics, chemical reac-ion kinetics, and mass balance, energy balance, momentum balance, of the models were givenout. The major results of models were also given. Finally, the comparison of computational results from the models with experimental or operational results were given It can be seen that theerrors of the main gas composition are very loKEY WORDS coal gasification, entrained-bed gasification, mathematical modeling(上接第6页)DEVELOPMENT OF COAL STRUCTURECheng Jun Zhou Anning and Li Jianwei( Department of Material Engineering, Xi'an University ofience and Technology, 710054, Xi'an)aBSTRaCt The coal aggregative structure is discussed, including the origin, structure andswelling behavior of several macerals. Coal chemicalH中国煤化工 I structure modeland synthesized models are analyzed and discussedCN MH Gructure research israised. At last. the utilization of coal structure research in new material field is summarizedwhich reveals the importance of coal structure research in ftKEY WORDS coal structure, physical structure, comprehensive structure, computer aidemolecular design