CHEMICAL INDUSTRY AND ENGINEERING PROGRESS 2016年第35卷第8期
·2426·
化 工 进 展
典型气流床煤气化炉气化过程的建模
东赫1，刘金昌1,2，解强1，党钾涛1 ，王新1
（1中国矿业大学(北京) 化学与环境工程学院，北京 100083；2九州大学电子和材料应用科学系，日本 福冈
春日 816-8580 ）
摘要：利用Aspen Plus 、基于热力学平衡模型对GSP 煤粉气化炉、GE 水煤浆气化炉及四喷嘴对置式水煤浆气化炉的气化过程建模。根据煤颗粒热转化的历程，将煤气化过程划分为热解、挥发分燃烧、半焦裂解及气化反应4个阶段，利用David Merrick 模型计算热解过程，采用Beath 模型校正压力对热解过程的影
响，选用化学计量反应器模拟挥发分燃烧反应，编制Fortran 程序计算半焦裂解产物收率，最后基于Gibbs 自由能最小化方法计算气化反应。结果表明，采用建立的气流床气化过程模型模拟工业气化过程的结果与生产数据基本吻合，对GSP 煤粉气化炉、GE 水煤浆气化炉及四喷嘴对置式水煤浆气化炉等3种气化炉有效气成分（CO+H 2）体积分数模拟结果的误差均不超过2%，建立模型的可靠性得到验证。 关键词：气流床气化炉；热力学平衡模型；Aspen Plus
中图分类号：TQ 546 文献标志码：A 文章编号：1000–6613（2016）08–2426–06 DOI ：10.16085/j.issn.1000-6613.2016.08.19
Modeling of coal gasification reaction in typical entrained-flow coal
gasifiers
DONG He 1，LIU Jinchang 1，
2，XIE Qiang 1，DANG Jiatao 1，WANG Xin 1
（1School of Chemical and Environmental Engineering ，China University of Mining and Technology （Beijing ），Beijing 100083，China ；2Department of Applied Science for Electronics and Materials ，Kyushu University ，6-1 Kasuga-Koen ，
Kasuga ，Fukuoka 816-8580，Japan ）reactor软件
Abstract ：This paper presents a modeling method for the coal gasification process proceeding in GSP pulverized coal gasification ，GE coal-water slurry gasification and Opposed Multiple-Burner gasification based on the thermodynamic equilibrium with the aid of Aspen Plus. In the light of thermal conversion procedure of fine coal particles ，the coal gasification was interpreted as consisting of four stages including pyrolysis ，volatile combustion ，char decomposition and gasification reaction. Then ，the pyrolysis stage was calculated by the David Merrick model and the effect of pressure on the coal pyrolysis was corrected by means of Beath model. The volatile combustion stage was simulated by using Rstoic reactor and the yield of char decomposition products was calculated via compiling Fortran program. And finally ，the gasification reaction stage was simulated based on the Gibbs free energy minimization. The results revealed that the simulated values from the developed simulation model of gasification processes were in good consistent with the industrial field data. The deviation of simulated results of volume fraction of the effective gas (CO+H 2) of these three typical entrained-flow gasifiers were all less than 2%，which can validate the reliability of the coal gasification model.
第一作者：东赫（1991—），女，硕士研究生。联系人：解强，教授，博士生导师。*********************** 。 收稿日期：2015-10-30；修改稿日期：2015-12-09。 基金项目：国家重点基
础研究发展计划（2014CB238905）及中央高校
基本科研业务费（2009KH10）项目。
第8期东赫等：典型气流床煤气化炉气化过程的建模·2427·Key words：entrained-flow gasifiers；thermodynamic equilibrium model；Aspen Plus
建立煤气化数学模型的核心在于选择合适的反应模型。迄今为止，煤气化数学模型繁多，主要为动力学模型和热力学平衡模型[1-2]。动力学模型虽能真实地反映炉内的气化反应过程，但该模型的建立依赖于详细的动力学反应机理及反应速率表达式，对于复杂的煤气化反应体系，其应用往往受限于清晰定义的反应和组分个数；另一方面，由于动力学模拟方法及软件的限制，仅能对气化炉单一设备进行模拟研究，难以考察气化前及气化后的整个工艺流程。相较而言[1]，热力学平衡模型不针对特定的气化炉型，以一定的平衡限制条件为假定基础，通过建立各元素组分的反应平衡、质量方程及整个气化炉热量平衡方程，求解煤气的平衡组成。热力学平衡模型设计过程简单，虽不能考察气化炉内的流动传热、传质特性及气化反应的过程，但由于求解速度很快且在选取准确的模型参数时能够获得较为准确的结果，故而在气化炉工程模拟和流程模拟中获得了广泛的应用[3]。
DAI等[4]基于Gibbs自由能最小化原则对煤粉气化过程建立模型，在此基础上考察了操作参数对气化过程的影响，并将模拟结果与实验数据进行比较，两者具有一致性；车德勇等[5]模拟流化床气化时将气化过
程分为热解及气化两个部分，并考虑流化床气化炉密相区和稀相区流体动力学参数不同，利用嵌入Fortran气化动力学子程序的2个搅拌釜式反应器模块分别模拟2个区内焦炭的气化反应；汪洋等[6]借助Aspen Plus软件建立了气流床煤气化炉的热力学平衡模型，其中包含裂解、气化及冷却等3个模块，将固体煤经裂解转化为Aspen Plus可处理的常规组分物质，经Gibbs平衡反应器求解产物组成，后经分离模块进行分离得到产品合成气和副产物灰渣。至今，基于热力学平衡模型对煤气化过程的模拟研究主要将气化过程简化为裂解单元和气化单元两部分，且侧重于操作参数对气化炉性能的分析，而关于热力学平衡模型详细的建立过程仍不充分，且实际的气化反应过程较为复杂，上述模型建立时尚未考虑煤气化过程伴随的热解及挥发分燃烧反应，难以系统和准确地反映实际气化过程。KONG等[7]针对水煤浆气化过程建立了三级平衡模型，将气化过程分为热解和燃烧级、气固相反应级及气相反应级3个等级进行模拟，结果表明该模型具有一定的可靠性；徐越等[8]将复杂的煤气化过程分为热解及挥发分燃烧、半焦燃烧气化反应和气化还原反应三部分进行建模，模型预测较为准确。此外，对于煤气化“拔头”工艺煤热解过程，基于Aspen Plus的模拟方法主要使用相关经验关联式计算热解产物组成，而所拟合的相关经验关联式多基于常压条件下的热解实验结果归纳[9]。但对于气流床气化炉，压力对热解过程的影响不可忽略。BEATH[10]研究发现煤热解析出的挥发分得率随着压力升高而减小，并根据实验数据拟合出压力与煤热解产物收率的关系式；LIU等[11]在对气流床气化过程建模过程中，采用Beath模型校正了热解过程，模拟结果更为贴近实际情况。
结合煤（颗粒）气化反应历程，基于热力学平衡模型对气化反应历程中各环节的模拟在以求解产物组成及相态为目的的情况下准确度高且计算量小。本文以实际生产指标为依据，利用Aspen Plus 软件，对3种典型气流床气化炉气化过程进行模拟，深入明晰气化炉煤热解、挥发分燃烧、半焦裂解及气化反应过程模型建立的机理及过程，并依据Beath 压力校准模型校准热解过程。
1 煤气化过程建模
1.1 煤（颗粒）气化历程解析
煤炭气化过程实际上是煤炭在高温下的热化学反应过程，王辅臣[2]、于遵宏[3]等对煤炭气化过程进行了详细分析，认为煤气化反应虽受气化方式、气化剂及煤种等诸多因素影响，但所有的气化过程存在明显的共性，即煤炭颗粒进入气化炉内先后经历了干燥、热解、挥发分燃烧、焦炭燃烧和气化反应过程。
对于气流床气化过程，在干燥和热解阶段，由于大部分煤颗粒小于200目，且炉内平均温度在1300℃以上，煤粉受热速度极快，因而认为干燥阶段煤中水分瞬间蒸发，同时发生快速的热分解脱除挥发分。目前气化过程模型的建立均将水分蒸发过程与挥发分析出过程一起考虑[12]。而后，热解析出的挥发分和形成的半焦进行燃烧和气化反应，此过程中挥发分的燃烧反应产生的热量以供为后阶段的气化反应，不可忽略。气化反应阶段包括气化剂与煤或焦之间的反应、煤与反应物之间的反应、反应产物与气化剂及反应产物之间的反应。
需说明，实际气化过程中各个阶段可能同时进
化 工 进 展 2016年第35卷
·2428·行，难以清晰分割。为便于气化过程模型的建立，结合煤气化反应历程，拟将复杂的煤气化过程进行分段，即划分为热解、挥发分燃烧、半焦裂解及气化反应4个阶段。 1.2 热解过程
1.2.1 David Merrick 模型
David Merrick 建立了煤热解挥发分及焦炭产量的计算经验模型[13]，该模型基于元素平衡计算热解挥发分的成分及焦炭的产率，计算简便且通用性较好[14]。假设挥发分产物有CH 4、C 2H 6、CO 、CO 2、H 2、H 2O 、NH 3、H 2S 和焦油，其他产物只有焦。具体方程如式(1)所示。式(1)是由10个方程组成的10阶矩阵，矩阵每行代表每个方程中10个回归系数。其中，第1行为碳平衡方程，2～5行分别为H 、O 、N 、S 平衡方程，第6行为质量平衡方程，7～10行表示各个物质的收率。
(C)0.750.80.42860.27270.850
000(H)0.250.2000.08210.11110.17650.0588(O)000.57140.72730.04900.888900(N)00000.0090000(S)00000.01000.83250.94121
00000000001000000000
01000000000010000000000100000⎡⎢⎢
⎢⎢
⎢
⎣
w w w w w 42622232CHAR C CH H C H O CO N CO S TAR 1H 1.31H H O 0.22H NH 0.32O H S 0.15O ⎡⎤⎤⎡⎤⎢⎥⎥⎢⎥⎢⎥⎥⎢⎥⎢⎥⎥⎢⎥
⎢⎥⎥⎢⎥⎢⎥⎥⎢⎥⎢⎥⎢⎥⎢⎥
∙=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥
⎢⎥⎢⎥⎢⎥
⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥
⎢⎥⎢⎥⎢⎥
⎢⎥⎢⎥⎢⎥⎦⎣⎦⎣⎦V - (1)
1.2.2 压力校正
气流床气化炉内压力较高，在使用热解模型对原煤热解产物进行计算后需要通过压力校正才能更加准确地模拟气流床气化过程。论文拟根据BEA TH [10] 实验结果提出的烟煤压力校准模型，如式(2)。
**1atm
0.13
=
P t V V P (2) 1.3 挥发分燃烧过程
使用David Merrick 模型可计算出的挥发分有CH 4、C 2H 6、CO 、CO 2、H 2、H 2O 、NH 3、H 2S 及焦油。其中可以发生燃烧反应的有CH 4、C 2H 6、CO 、
H 2、
H 2S 及焦油。假定气化反应过程中气化剂充足，上述五种可燃挥发分以及焦油发生了完全燃烧反应，采用Aspen Plus 软件现有的操作模块化学计量反应器（RSTOIC ）进行模拟，主要反应方程式如式(3)～式(7)。
2C H (0.5)O ++m n m n ==22CO 0.5H O +m n (3)
42CH +2O ==22CO +2H O (4) 22H +0.5O ==2H O (5) 22H S+1.5O ==22H O+SO (6) 262C H +3.5O ==222CO +3H O (7) 1.4 煤(焦)裂解过程
煤、半焦都是复杂的混合物，在模拟软件中无法对组分不清楚的混合物进行计算，将热解后的半焦依据其工业分析和元素分析结果分解为单元素分子（C 、H 2、O 2、N 2、S ）和灰分，裂解组分的收率
通过Aspen Plus 内嵌的Calculator 计算模块中Fortran 语言编写的程序来控制。该单元本身遵循质量守恒，只需考虑稳态情况下系统能量的守恒，且由于煤气化是化学反应过程，因而可采用生成焓来计算反应中气体和固体的焓。煤裂解单元的焓平衡方程及生成焓[15]如式(8)、式(9)。
00coal f ,coal,298coal coal feed f ,prod,298,1
prod
1
()()==∆+∆=∆+
+NP
i i i NP
i
i
P
i m H m H T n H n H T
Q ∑∑ (8)
f,coal,298ar ar ar ar (327.861418.7992.84158.67–)
∆=+++H HHV C H S M (9)
式(8)表示在一定温度T 下，等号左边的两项分别表示进料物流的生成焓和总焓，等号右边前两项分别表示反应产物的生成焓和总焓。
在Aspen Plus 软件模拟时，半焦裂解过程产物依据焦的工业分析和元素分析，通过Fortran 编程语言实现，如下所示[16]。
FACT=(100–WATER)/100 H 2O=WATER/100
ASH=ULT(1)/100*FACT C=ULT(2)/100*FACT H 2=ULT(3)/100*FACT N 2=ULT(4)/100*FACT Cl 2=ULT(5)/100*FACT S=ULT(6)/100*FACT O 2=ULT(7)/100*FACT
第8期 东赫等：典型气流床煤气化炉气化过程的建模
·2429·
其中，Fortran 语言中ULT 为矢量，代表煤的元素分析（ULTANAL ）。 1.5 煤(焦)气化过程
气化反应单元采用平衡反应器RGIBBS 模块，该模块根据系统Gibbs 自由能趋于最小原则，由于气流床煤气化本身温度较高，约达1500℃，压力约为4MPa ，因此可认为反应体系在瞬间达到化学反应平衡，在考虑热损失的前提下，通过建立下述非线性方程组计算同时达到化学平衡和相平衡时的系统组成和分布[17]。
目标函数如式(10)。
1
11
min ,==+==+s
C
P
c j J jl jl J j s L G G n G n G ∑∑∑ (10)
质量守恒约束条件如式(11)。
1
11
1,2,3,,==+==+= s
C
P
c
k j
jk jl jk
j j s l b n m n m k Z ∑∑∑ (11)
焓平衡约束条件如式(12)。
0f ,in,298,,(in,)P 1
1
f ,out,298,,(out,)L
1
1
====∆+∆+=
∆+∆+C C
i i i i T i i i C
C
i
i i i T i i i n H n H Q n H
n H Q ∑∑∑∑(12)
非负约束条件如式(13)。
0i n ≥ (13)
2 Aspen Plus 模拟流程
2.1 建模假设条件
结合煤气化反应特征及Aspen Plus 模拟的适用性，建模需作如下假设[7]。
（1）整个气化反应可以分为热解、挥发分燃烧、半焦裂解和气化4个步骤进行；
（2）热解过程不受气化炉内气氛条件的影响，
热解产物主要考虑焦、焦油（C m H n ）
、CH 4、CO 、CO 2、H 2、H 2O 、N 2和H 2S ；
（3）挥发分燃烧过程中，由热解产生挥发分的可燃成分在富氧的条件下完全燃烧；
（4）热解产生的半焦根据工业分析和元素分析结果裂解为C 、H 、O 、N 、S 的单质和灰分进行气化反应；
（5）气化炉处于稳定运行状态，所有参数不随时间发生变化，无压力降，且不考虑气化炉内的浓度梯度和温度梯度；
（6）半焦中的C 、H 、O 、N 、S 等元素经过反应之后全部转化为气相，而C 根据反应条件的不同不完全转化；其中ASH 为惰性组分，不参与气化反应；
（7）整个气化过程都瞬间完成，且所有煤粒都参与到气化反应中。 2.2 模型选择
建模过程所选用的Aspen Plus 中的模块如表1所示。
2.3 模拟流程
基于前述假设，选择适于处理高压烃热力学性质的PR-BM 状态方程[18]，采用Aspen Plus 现有的单元操作模块对整个煤气化过程进行建模，整个气化过程模拟流程如图1所示。
（1）热解过程 原料（COAL ）和载气（CARRIERG ）进入混合器（MIX1）中混合，混合之后（FEED ）煤粉在RYIELD 反应器（PYROLY1）中发生常压下的热解反应生成挥发分和半焦的混合物P1。P1进入RYIELD 反应器（PYROLY2）中进行热解压力校正，得到热解产物P2。P2经过分离器（SEP2）进行气固相分离，得到挥发分（GAS ）、半焦（CHAR ）和焦油（TAR ）。由于气流床煤气化反应在高温高压富氧条件下进行，热解产生的焦油瞬间发生了完全燃烧反应，因而将热解产生的焦油作为气相的热解挥发分处理。
（2）挥发分燃烧过程 挥发分（GAS ）氧化剂
表1 Aspen Plus 模块应用汇总
单元
模型
功能
作用
MIX1 MIXER
物流混合器，将多股物流汇合成一股物流 煤粉与载气混合，模
拟载气输送煤粉过程；煤粉与水混合，模拟水
煤浆的制备过程
PYROLY1RYIELD
产率收率反应器，在产品收率已知的情况下模拟化学反应过程
根据David Merrick 模型计算热解产物，模拟煤常压热解过程 PYROLY2RYIELD
产率收率反应器，在产品收率已知的情况下模拟化学反应过程 根据挥发分压力校正方程计算煤热解挥发分产物在实际气化压力下的收率 SEPARATE SEP2
组分分离器，根据规定的流率、分率或纯度，将一股物流分成两股物流
热解产物中气固相分离 COMBUST RSTOIC
化学计量转化反应器，根据反应计量学计算化学反应 根据化学反应方程模拟挥发分的燃烧反应 DECOMP RYIELD
产率收率反应器，
在产品收率已知的
情况下模拟化学反
应过程
根据半焦的工业分
析和元素分析结果模
拟半焦的裂解反应
MIX2 MIXER
物流混合器，将多股物流汇合成一股物流 模拟半焦裂解产物、气化剂和挥发分燃烧产物的混合过程
GASIFER RGIBBS
多相化学平衡反应
器，根据Gibbs 自由能最小化法计算气化过程
根据Gibbs 自由能最小化法模拟计算气化过程
化 工 进 展 2016年第35卷
·2430
·
图1 气流床煤气化过程模拟流程图
（OXYGEN ）混合进入到RSTOIC 反应器（COMBUST ）中发生燃烧反应，得到产物（GASPRO ）。
（3）半焦裂解过程 半焦（CHAR ）进入到RYIELD 反应器（DECOMP ）中进行裂解反应，得
到C 、H 2、O 2、N 2、S 及ASH 混合物（DECOMP ）
。 （4）气化过程 挥发分燃烧产物（GASPRO ）、半焦裂解产物（DECOMPUD ）及水蒸气（STEAM ）首先进入混合器（MIX2）中混合，混合物（MIXTURE ）进入RGIBBS 反应器（GASIFIER ）中进行气化，生成产物（PRODUCT ）。
3 模拟与校核数据
3.1 分析数据
（1）煤质分析数据 选用宁东重化工基地某 煤气化厂中GSP 、GE 及四喷嘴对置式水煤浆气化炉实际运行过程的气化原煤，其主要煤质指标如表2所示。
表2 煤的工业分析和元素分析 单位：%
M ad
A d
V d
FC d
C d
N d
S t,d
H d
O d 7.22 17.04 28.51 54.45 66.02
0.73 0.97 4.51
10.73
（2）半焦分析数据 利用David Merrick 模型估算热解产物需用到半焦的工业分析和元素分析数据，将原煤通过热解实验制得半焦，热解后的半焦工业分析与元素分析数据如表3所示。 3.2 热解模型计算
根据David Merrick 模型计算热解产物的结果如表4所示；根据Beath 模型校准的热解产物结果如表5所示。
表3 半焦的工业分析和元素分析 单位：%
M ad
A d
V d
FC d
C d
N d
S t,d
H d
O d 0.58 23.98 9.85 66.17 70.26
0.58 0.86 1.23
3.09
表4 热解产物组成
产物 产率（质量分数）/%
产物
产率（质量分数）/%
半焦 70 CO 2 1.93 焦油
11.26 H 2 2.13
CH 4 7.07 H 2O 2.17 C 2H 6 1.19 NH 3 0.015 CO 4.12 H 2S 0.114
表5 Beath 模型校准的热解产物组成
产物 产率（质量分数）/%
产物
产率（质量分数）/%
半焦 58.5 CO 2 1.61 焦油
9.4 H 2 1.78
CH 4 5.9 H 2O 1.81
C 2H 6 0.99 NH 3 0.01 CO 3.44 H 2S 0.1
3.3 气化模型验证
应用建立的气化模型进行模拟的结果汇总于表6。表6也给出了气化生产的实际数据。
由表6可知，模拟结果与实际工业生产数据基本吻合，GSP 、GE 及四喷嘴对置水煤浆气化炉合成气有效成分（CO+H 2）体积分数模拟结果误差分别为1.48%、0.53%及0.31%，说明该模型能够较为准
表6 GSP 煤粉气化炉、GE 水煤浆气化炉及四喷嘴对置式
水煤浆气化炉模拟结果验证
组分（体积分数） /%
GSP GE
四喷嘴对置式水
煤浆气化炉
生产数据模拟结果生产 数据 模拟 结果 生产数据模拟结果
H 2 25.5 25.8837.30 36.18 33.7733.96CO 69.25
67.39
39.95 41.60
46.5846.88CO 2 4.42 3.92 21.59 15.33
18.9612.60
CH 4 0.05 0.08 0.25 0.20 0.21 0.15 N 2 0.69
0.38 0.69 0.44 0.31 0.33 H 2S 0.09 0.12 0.16 0.19 0.16 0.12 CO+H 2
94.75
93.27
77.25 77.78
80.3580.04
版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。
发表评论