褐煤直接燃烧会造成严重环境污染, 必须提质后才能加以利用。褐煤提质技术包括LFC[1-2]、Toscoal[3-4]、COED[5]、Lurgi[6]、Garrett[7]和ECOPRO[8]等工艺, 其中, 美国的LFC技术流程简洁, 是目前商业化程度最为成熟的褐煤提质技术。运用流程模拟软件Aspen Plus对LFC工艺进行模拟, 可以准确地描述整个过程的运行情况, 对工程设计和装置的运行具有一定的指导意义。㶲分析方法是基于热力学第二定律的分析方法, 可以用来计算系统各操作单元的㶲损失及㶲效率, 直观反映整个系统的用能情况[9-11]。李超[12]等采用Aspen Plus对煤制天然气的气化-热解耦合系统进行了模拟并进行㶲分析, 考察了操作条件对耦合系统气化性能的影响, 提出了优化的操作条件。邱朋华等[13]通过对准东煤气化过程的模拟和㶲分析, 比较了相同气化条件下准东煤与烟煤经济性的差异。Chen等[14]采用Aspen Plus对一种低阶煤低温热解系统进行了流程模拟并分析了系统的㶲效率, 提出了改进方法。张子昂[15]等对煤气化合成甲醇进行了㶲分析, 在流程模拟的基础上计算了煤制甲醇的㶲效率为53.98%, 供热㶲效率为2.96%, 㶲计算结果表明燃烧和化学反应㶲损很大, 并提出了改进建议。笔者通过对LFC实验装置得到的初始数据在流程模拟的基础上进行㶲分析, 计算系统中各操作单元的㶲损失, 分析㶲损失产生的原因和降低㶲损失的方法。
1 LFC工艺LFC工艺是一种褐煤轻度热解工艺。LFC工艺流程见图 1。褐煤经破碎和筛分后送入干燥炉, 被来自干燥燃烧炉的热烟气加热脱除水分, 离开干燥炉的气体经旋风除尘后分为两部分, 一部分与干燥燃烧炉进口气体混合作为循环气体, 另一部分送入烟气脱硫装置, 经脱硫后排入大气。干煤进入热解器, 与来自燃烧炉的高温烟气接触, 发生热解反应, 热解气经过旋风分离器分离出煤粉后, 经过冷却、焦油静电捕集器得到焦油, 剩余的不凝气分成三部分, 一部分作为循环气体, 一部分进入干燥燃烧炉, 一部分进入热解燃烧炉, 为干燥、热解过程提供热量, 煤气燃烧提供的热量并不够, 须补充燃料甲烷。固体产物经过冷却、钝化后形成产品半焦。
为了得到流程模拟需要的初始数据, 进行褐煤低温热解提质实验。褐煤热解提质装置示意图见图 2(处理量为10.00 kg/h)。氮气作为热载气, 干燥器和热解器温度由电炉子来控制, 煤在其中的停留时间由螺旋给料器的转速控制。
内蒙古芒来褐煤储量巨大, 实验中采用这一典型褐煤作为原料煤, 控制干燥器内为200 ℃, 热解器内为520 ℃, 停留时间均为1 h, 原煤、半焦和焦油的性质分析见表 1。
化工流程模拟软件Aspen Plus模拟系统广泛用于化工和石油化学工业过程模拟[16-18]。以1 000.00 kg/h原煤处理量工业化装置为对象, 运用Aspen plus流程模拟平台对LFC工艺进行流程模拟。
模拟过程作以下假设:
(1) 过程处于稳定运行状态, 不考虑开车和停车状态。
(2) 原煤的干燥过程作为化学反应对待, 以方便实现模拟干燥过程。
(3) 灰分在热解过程中作为惰性物质, 不与其他化合物发生反应, 热解结束后, 灰分都富集在半焦中。
(4) 过程中的气体都按照理想气体对待, 方便计算。
(5) 在相同模块内, 气相和固相的温度相同。
(6) 煤气包括H2、CH4、CO、CO2、C2H4、C2H6、C3H8。
原料煤、半焦等固体设定为非常规固体, 它们的属性可以通过工业分析和元素分析确定, 而对于焦油这种极端复杂的液体, Aspen中没有与之对应的物质类型, 只能以某些常规的化合物替代焦油。由于苯环是焦油中主要成分, 在Aspen Plus软件中可用苯代替焦油[19]。Ji等[20]用C6H6O模拟焦油, Yan等[21]用一组化合物C12H26S、C15H33N、C14H12O2和C10H24O4代表焦油, Yi等[22]用另外一组化合物C7H8(<170 ℃)、C6H6O(172~210 ℃)、C10H8(210~230 ℃)、C12H10(230~300 ℃)、C14H10(300~360 ℃)和C16H10(>360 ℃)代表褐煤的煤焦油进行模拟, 用化合物和所在馏分段的质量分数模拟煤焦油, 本文中采用此方法模拟芒来褐煤生产煤焦油。
LFC工艺模拟流程图如图 3所示, 整个工艺主要分为4个操作单元:干燥单元, 热解单元, 焦油处理单元和燃烧供热单元。由于流程中含有常规组分、非常规组分以及固体物料, 选择Aspen Plus中固体类型为MIXNCPSD, 流量基准为质量流量, 系统压力为101 kPa。采用HCOALGEN模型和DCOALIGT模型计算非常规组分的焓和密度, PR-BM方法为流程模型的全局物性方法。干燥器选用Rstoic(化学计量反应器)模块, 利用此反应器不考虑反应动力学数据, 只规定反应器内所发生的化学反应的反应程度或者转化率来模拟反应过程。煤干燥过程的化学反应方程式可以表示为:coal(湿)→coal(干)+0.055H2O。热解反应非常复杂, 热解器选用Ryield(收率反应器)模块, 利用此反应器既不须考虑反应动力学数据, 也不须考虑所发生的化学反应, 仅规定反应器内褐煤热解生成各产物的收率。焦油处理单元主要包括冷却塔和电捕焦油器, 分别用RadFac(严格精馏)模块和SEP型分离模块来模拟。燃烧供热单元由RGibbs反应器来模拟, RGibbs反应模块基于吉布斯自由能最小化的原理, 计算同时达到化学平衡和相平衡时的系统组成和相分布, 最适合模拟复杂混合物的燃烧过程。
为了保证模拟结果的可靠性, 对比了热解产品产率的实验值和模拟值, 结果见表 2。可以看出, 实验值和模拟值非常接近, 平均相对误差仅为1.45%。
从干燥器出来的气体中有58.55%作为循环气体, 41.45%作为烟气排入大气。从电捕焦油器中出来的气体中32.50%的进入了干燥用燃烧炉, 57.29%进入热解用燃烧炉, 10.21%的气体作为循环气。由于煤气产生的热量不能满足褐煤的干燥和热解, 因此向两个燃烧炉中分别补充甲烷, 热解用燃烧炉的补充的甲烷为24.50 kg/h, 干燥用燃烧炉补充的甲烷为13.30 kg/h。当原料煤处理量为1 000.00 kg/h时, 需要空气为1 708.00 kg/h, 须补充甲烷燃料37.80 kg/h。生产焦油66.31 kg/h, 半焦475.63 kg/h。
通过模拟得到各个操作单元流入、流出物流股的温度、流量、组成等参数, 结果见表 3。
㶲不仅能反映热力系统发展过程中能量在数量上的变化, 还能反映其在质量上的变化, 它的表达式为热力学参数焓和熵的函数。稳定流动系统的㶲主要由动能㶲、位能㶲、物理㶲和化学㶲4部分组成。LFC工艺中只考虑物理㶲和化学㶲, 动能㶲和位能㶲忽略不计。
4.1 㶲值计算物理㶲(
$\begin{align} Ex^{\text{PH}}=(H-H_{0})-T_{0}(S-S_{0}). \end{align}$ | (1) |
式中, H和S为流体处于某状态的焓和熵; H0和S0为流体在基准态下的焓和熵。
物流的物理㶲值由下式计算:
$\begin{align} \dot{E}x^{\text{PH}}=∑_i\dot n_iEx^{\text{PH}}_i. \end{align}$ | (2) |
化学㶲是由于系统与环境发生物质交换或化学反应, 达到与环境平衡时所具有的有效能值。气相物流的化学㶲由下式计算:
$\begin{align} \dot{E}x^{\text{CH}}=∑_i\dot n_i\left(Ex^0_i+RT_0\text{ln}\frac{\dot n_i}{∑\dot n_i}\right). \end{align}$ |
式中,
煤、半焦和焦油所含化合物非常复杂, 同时又缺乏热力学参数, 不能采取常规方法计算化学㶲。Eisermann[23]介绍了一种碳基液体和固体燃料化学㶲的计算方法, 即
$\begin{align} Ex^{CH}=ξQ_{LHV} , \end{align}$ | (4) |
$\begin{align} ξ=1.040 1+0.172 8w_{\text{H/C}}+0.043 2w_{\text{O/C}}+\ 0.216 9w_{\text{N/C}}(1-2.062 8w_{\text{S/C}}). \end{align}$ | (5) |
式中, QLHV为燃料的低位燃烧热, kJ/kg; w为质量分数。
物流的化学㶲用下式计算:
$\begin{align} \dot Ex^{\text{CH}}=\dot m_{i}Ex^{\text{CH}}, \end{align}$ | (6) |
稳流系统的㶲平衡方程式为
$\begin{align} \dot Ex_{\text{in}}=\dot Ex_{\text{loss}}+\dot Ex_{\text{out}}. \end{align}$ | (7) |
㶲效率定义为:在系统的能量传递和转换过程中, 能被有效利用的㶲与总消耗的比值, 即
$\begin{align} η=\frac{\dot Ex_{\text{out}}}{\dot Ex_{\text{in}}}. \end{align}$ | (8) |
LFC各个操作单元物流股的㶲值见表 4, 各个单元㶲损失和进出物流的㶲见图 4, 图 5为各个单元的㶲损失所占系统总㶲损失比例。燃烧供热单元的㶲损失最大达到2.48×106 kJ/h, 㶲效率为44.76%, 是由于烟气、甲烷和空气的混合和燃烧的不可逆性引起的㶲损失, 㶲损失占系统总㶲损失的64.00%。干燥单元㶲效率为96.03%, 㶲损失排在第二位。在干燥器中, 当煤的温度达到200 ℃, 水变成蒸汽离开煤, 另外一部分羧基会分解成CO2, 因此既有物理变化又有化学变化, 物理变化是主要的, 水的气化和分离的不可逆性是造成㶲损失的主要原因。热解单元㶲效率为98.55%, 㶲损失为2.84×105 kJ/h, 主要是由于热解反应的不可逆性产生的㶲损失。热解单元和干燥单元拥有比较高的㶲效率, 说明干燥热解是比较好的褐煤提质方法。焦油处理单元㶲效率为58.25%, 㶲损失为4.04×105 kJ/h, 主要是要由于热传递的不可逆性带来的㶲损失。LFC工艺总的㶲损失为3.97×106 kJ/h, 总㶲效率为74.39%。
燃烧供热单元中, 减小烟气、甲烷和空气燃烧前后的温度差可以减少由于温差、混合产生的㶲损失。在热解气进入焦油处理单元之前增加换热装置, 把一部分热量转移到甲烷、空气中, 提高其温度, 可以提高燃烧供热单元的㶲效率, 同时减轻焦油处理单元的热负荷。改善干燥器和热解器的传质传热效果能够降低热载气的使用量, 可以提高干燥、热解单元的㶲效率。半焦离开热解单元时约为460 ℃, 具有较高的物理㶲, LFC工艺中并没有充分利用这部分物理㶲。向半焦中喷淋适量的水, 半焦会迅速降温终止热解反应, 水被加热气化后输送至燃烧炉中, 热量得到重复利用, 从而提高整个系统的㶲效率。
5 结论(1) 当原料煤的处理量为1 000.00 kg/h时, 需要空气为1 708.00 kg/h, 须补充甲烷燃料37.80 kg/h, 生产焦油66.31 kg/h, 半焦475.63 kg/h。
(2) 主要㶲损失是由于传热不可逆性和化学反应的不可逆性引起的, LFC工艺的㶲效率为74.39%, 燃烧供热单元的㶲损失最大, 㶲效率最低。
(3) 通过改进工艺、回收利用半焦和热解气的物理㶲提高系统的㶲效率。
[1] |
ENCOAL Corporation. ENCOAL mild coal gasification project: project performance summary[R]. USA: USDOE, 1999.
|
[2] |
ENCOAL Corporation. ENCOAL mild coal gasification project: ENCOAL project final report[R]. USA: USDOE, 1997.
|
[3] |
MARK T A, BERNARD L S. The toscoal process-pyrolysis of western coals and lignites for char and oil production[J]. Preprint Papers American Chemical Society, Division of Petroleum Chemistry, 1977, 22(2): 233-252. |
[4] |
CORTEZ D H, LADELFA C J. Coproduction of syncrude and power using the toscoal process[C]//Pan-Pacific Synfuels Conference, Japan Petroleum Institute, 1982:331-348.
|
[5] |
STROM A H, EDDINGER R T. COED plant for coal conversion[J]. Chemical Engineering Progress, 1971, 67(1): 75-80. |
[6] |
RAMMLER R W, LUIGI K. Synthetic fuels from Lurgi coal pyrolysis[J]. Energy Program, 1982, 2(2): 121-129. |
[7] |
SASS A. Garrett coal pyrolysis process[J]. Chemical Engineering Progress, 1974, 70(1): 72-73. |
[8] |
YABE H, KAWAMURA T, KOZURU H, et al. Development of coal partial hydro-pyrolysis process[J]. Technology of Nippon Steel Corporation, 2005, 92(3): 8-15. |
[9] |
赵冬, 冯霄, 王东亮. 煤制天然气过程模拟与㶲分析[J]. 化工进展, 2015, 34(4): 990-996. ZHAO Dong, FENG Xiao, WANG Dongliang. Simulation and exergy analysis of coal to SNG process[J]. Chemical Industry and Engineering Progress, 2015, 34(4): 990-996. |
[10] |
孙丽萍, 刘月娥, 朱文坚, 等. 煤制气甲烷化过程能量利用评价方法比较[J]. 能源化工, 2015, 36(3): 27-31. SUN Liping, LIU Yue'e, ZHU Wenjian, et al. Comparison of evaluation methods of energy utilization in the syngas methanation process[J]. Energy Chemical Industry, 2015, 36(3): 27-31. |
[11] |
程文伟, 于戈文, 张文兴, 等. GSP干粉煤气化系统的热力学分析[J]. 山西煤炭, 2014, 34(5): 21-25. CHENG Wenwei, YU Gewen, ZHANG Wenxing, et al. Comparison of evaluation methods of energy utilization in the syngas methanation process[J]. Shanxi Coal, 2014, 34(5): 21-25. |
[12] |
李超, 代正华, 杨骥, 等. 用于煤制天然气的气化热解耦合系统模拟研究及能量分析[J]. 燃料化学学报, 2015, 43(7): 779-789. LI Chao, DAI Zhenghua, YANG Ji, et al. Modelling and energy analysis of an integrated coal gasification and pyrolysis system for synthetic natural gas[J]. Journal of Fuel Chemistry and Technology, 2015, 43(7): 779-789. |
[13] |
邱朋华, 陈希叶, 李丹丹, 等. 准东煤气化过程的热力学分析[J]. 煤炭转化, 2015, 38(1): 33-38. QIU Penghua, CHEN Xiye, LI Dandan, et al. Thermodynamic analysis for gasification process of Zhundong coal[J]. Coal Convertion, 2015, 38(1): 33-38. |
[14] |
CHEN Xiaohui, ZHENG Danxing, GUO Jing, et al. Energy analysis for low-rank coal based process system to co-produce semicoke, syngas and light oil[J]. Energy, 2013, 52(52): 279-288. |
[15] |
张子昂, 马素霞. 分级给氧膜式壁气化炉煤基联产系统能耗分析[J]. 煤炭转化, 2014, 37(2): 57-62. ZHANG Zi'ang, MA Suxia. Energy analysis of polygeneration on systems of the stage-feed oxygen entrained-flow gasifier with water-cooled wall[J]. Coal Conversion, 2014, 37(2): 57-62. |
[16] |
李军, 王汝军, 王丁丁, 等. 差压热耦合反应精馏应用于乙酸甲酯水解研究[J]. 中国石油大学学报(自然科学版), 2014, 38(1): 155-160. LI Jun, WANG Rujun, WANG Dingding, et al. Simulation study of different pressure thermally coupled reactive distillation column for hydrolysis of methyl acetate[J]. Journal of China University of Petroleum(Edition of Natural Science), 2014, 38(1): 155-160. |
[17] |
ADEYEMI I, JANAJREH I. Modeling of the entrained flow gasification: kinetics-based ASPEN plus model[J]. Renewable Energy, 2015, 82: 77-84. DOI:10.1016/j.renene.2014.10.073 |
[18] |
MA Youfu, ZHANG Hua, YUAN Yichao, et al. Optimization of a lignite-fired open pulverizing system boiler process based on variations in the drying agent composition[J]. Energy, 2015, 81(73): 304-316. |
[19] |
ZHANG R Z, LIU C Y, YIN R H, et al. Experimental and kinetic study of the NO-reduction by tar formed from biomass gasification, using benzene as a tar model component[J]. Fuel Process Technology, 2011, 92(1): 132-138. DOI:10.1016/j.fuproc.2010.09.016 |
[20] |
JI P, FENG W, CHEN B. Production of ultrapure hydrogen from biomass gasification with air[J]. Chemical Engineering Science, 2009, 64(3): 582-592. DOI:10.1016/j.ces.2008.10.015 |
[21] |
YAN H M, ZHANG D K. Modeling of a low temperature pyrolysis process using Aspen Plus[J]. Developments in Chemical Engineering and Mineral Processing, 1999, 7(5): 577-591. |
[22] |
YI Qun, FENG Jie, LU Bingchuan, et al. Energy evaluation for lignite pyrolysis by solid heat carrier coupled with gasification[J]. Energy Fuels, 2013, 27(8): 4523-4533. DOI:10.1021/ef400865h |
[23] |
EISERMANN W, JOHNSON P, CONGER W L. Estimating thermodynamic properties of coal, char, tar and ash[J]. Fuel Process Technology, 1980, 3(1): 39-53. DOI:10.1016/0378-3820(80)90022-3 |