专栏名称: 科学温故社
知幽渺之理而收显著之效:从实验到理论,探索催化反应隐藏的信息与机理,打开催化过程“黑匣子”
目录
相关文章推荐
51好读  ›  专栏  ›  科学温故社

【解读】化工进展:基于Aspen Plus的甲烷联合重整制合成气过程热力学计算

科学温故社  · 公众号  ·  · 2024-07-29 08:15

正文




文章信息




基于Aspen Plus的甲烷联合重整制合成气过程热力学计算

庞淑馨 1 ,王昊 2 ,王健宇 3 ,朱卡克 1 , 刘志成 4

1 华东理工大学化学工程联合国家重点实验室,上海 200237; 2 中石化齐鲁分公司研究院,山东 淄博 255400; 3 国防科技大学气象海洋学院,湖南 长沙 410005; 4 中石化(上海)石油化工研究院,绿色化工与工业催化国家重点实验室,上海 201208

引用本文

庞淑馨, 王昊, 王健宇, 等. 基于Aspen Plus的甲烷联合重整制合成气过程热力学计算[J]. 化工进展, 2024, 43(5): 2890-2900.

DOI: 10.16085/j.issn.1000-6613.2023-2118


摘要

甲烷联合重整产物合成气的组成受热力学平衡控制,反应中积炭反应也受操作条件的影响,因此,计算反应条件(原料组成、温度、压力)对产物组成和积炭形成的影响,进行定量热力学分析是工艺设计所需要的。本文使用Aspen Plus软件中的RGibbs模块,建立甲烷联合重整热力学反应模型,采用吉布斯自由能最小化法,考 虑所有相关反应的进料组成、温度和压力对平衡反应混合物组成和积炭的影响。首先,设定化学计量比 CH 4 ∶(CO 2 +H 2 O)=1∶1进料,调控CO 2 、H 2 O相对进料比例,结果发现:CO 2 /CH 4 比增加时,CH 4 和CO 2 转化率升高,但积炭量也相应增加,H 2 /CO比减小;CH 4 和CO 2 的平衡转化率随温度升高而升高,且随着温度升高趋于极限,可通过调节温度从而调控H 2 /CO比达到所需生产要求。然后,研究了压力对反应平衡的影响,结果表明CH 4 和CO 2 的热力学平衡转化率随压力升高而降低,积炭量增加。最后,结合进料比、温度和压力的影响,针对合成气制下游化学产品的不同需求以及特定氢碳比精准调控,模拟计算出合适的反应进料比例,为将来实际生产调控提供了理论计算依据。


甲烷与二氧化碳重整催化反应制合成气,又称甲烷干重整(DRM),是氨、甲醇和费-托合成工艺等的通用原料,被认为是CH 4 和CO 2 这两种温室效应气体转化为高附加值化学品最关键的技术之一。DRM与甲烷水蒸气重整(SRM)相结合的联合重整(BRM)可使用富CO 2 的天然气、沼气进料,还可用于回收来自化石燃料(煤、石油、天然气等)燃烧的废气和各种工业排放的CO 2 ,降低分离操作成本,比传统SRM更具经济性。与DRM相比,BRM具有几方面优势:①富氧化反应体系中的C*更易氧化从而减少了积炭的形成;②通过调节进料中H 2 O和CO 2 的组成,可以生产不同比例的H 2 /CO合成气,应用于不同的下游生产。此外,Choudhary等在甲烷重整反应中同时加入H 2 O和CO 2 ,发现甲烷几乎可以完全转化为合成气,对CO和H 2 的选择性接近100%。工业应用最广泛的甲烷重整技术是SRM,但由于SRM高吸热的特性,目前在寻找可替代工艺。BRM工艺被认为是商用水蒸气重整技术的替代方案,可将能源成本和温室气体排放降至最低。本研究旨在通过将DRM与SRM技术相结合,克服DRM的局限性。BRM反应出口浓度接近平衡组成,通过热力学平衡计算,针对不同的下游工艺,分析在满足工业合成气要求的前提下,寻找较优进料组成及对应的反应温度与压力,为实际生产应用提供理论指导与数据支持

BRM工艺的研究主要集中在催化剂性能的研究上,由于其反应条件 (如反应温度、压力、进料组成等)差异较大,难以对反应条件进行系统分析。另一方面,BRM在接近热力学平衡的操作条件下运行,产物具有热力学平衡特点。因此考察温度、压力和进料组成等对平衡转化率。产物分布和积炭过程的影响,对相关工艺和催化过程开发具有重要指导作用。近年来,利用吉布斯自由能最小化技术计算甲烷重整制合成气的热力学平衡结果已被报道。Amin和Yaw利用拉格朗日待定乘数法对DRM和甲烷部分氧化联合重整制合成气的过程进行热力学平衡分析,计算得到进料参数对产物分布的影响,认为最佳进料条件为CH 4 ∶CO 2 ∶O 2 为1∶(0.8~1.0)∶(0.1~0.2),但计算过程没有考虑积炭的影响。Li等分别研究了甲烷重整过程CO 2 、H 2 O和O 2 组成变化和温度对反应的影响,发现在重整过程中提高反应温度,增加合适的进水量,可以达到减少积炭生成的目的。Song和Pan研究了甲烷三重整,在单个反应器中协同组合DRM、SRM和甲烷部分氧化重整,探究了不同进料组成对反应的影响,计算过程未涉及压力影响。Yang等在研究BRM的过程中,设计开发了联合重整串并联反应器,为后续的技术改进、降低成本和环保评估等提供了参考。目前针对联合重整制合成气将温度、压力和进料组成系统考察的热力学分析较少,旨在提供用于下游产品设计合适的合成气生产工艺开发的热力学数据

本文利用Aspen Plus软件,首先建立CH 4 与CO 2 和H 2 O联合重整反应的反应模型,并采用吉布斯自由能最小化原则,建立热力学计算方法。然后,在CH 4 ∶(CO 2 +H 2 O)=1∶1的进料条件下,改变CO 2 与H 2 O的进料组成,计算不同反应条件下的产物组分。接着,结合重整体系中反应压力、积炭极限的计算分析,增加H 2 O进料量,建立最优反应条件,以满足工业生产的CH 4 干气出口浓度及催化剂抗积炭要求。最后,对所需实际应用合成气衍生品的生产要求模拟计算对应上游工艺反应条件,为未来BRM反应工艺设计提供参考


1

实验模型和计算方法

1.1

联合重整工艺实验模型

本模型建立采用Aspen Plus软件中的 RGibbs模块,RGibbs反应器根据系统的吉布斯自由能趋于最小值的原则,计算同时达到化学平衡和相平衡时的产物组成和相分布,该方法不需要知道具体的反应方程式和反应动力学。RGibbs反应器是唯一能处理气液固三相平衡的反应器设计模块

以RGibbs模块作为反应器建立BRM反应模型,CH 4 、CO 2 和H 2 O作为进料共同注入混合器中,接着送入反应器,设置初始反应参数,通过模型分析工具建立灵敏度分析,选取反应温度、反应压力与进料中各组分比例分别作为变量,定义变量后计算,将所得计算结果进行整合并作图。Aspen Plus流程图如图1所示,S1为水进料,经过预热段B1转化为水蒸气进入混合器B2中,S2、S3分别作为CH 4 与CO 2 的进料段,进入混合器B2后与水蒸气共同进入反应器B3中,反应结束后,合成气由S6段排出

图1 Aspen Plus计算流程图

S1—水进料;S2—CH 4 进料;S3—CO 2 进料;S4—水蒸气;S5—混合气进料;S6—气体排出;B1—水蒸发器;B2—混合器;B3—Gibbs反应器

1.2

联合重整反应网络

BRM反应的主反应有两个,分别是DRM反应CH 4 -CO 2 生成H 2 /CO=1的合成气[式(1)],SRM反应CH 4 -H 2 O反应生成H 2 /CO=3的合成气[式(2)],将两式相结合得到BRM反应式[式(3)]
副反应包括水煤气变换[WGS,式(4)]和积炭反应,沉积在催化剂活性位点上的石墨碳会导致严重的反应器堵塞和造成催化剂失活。这种碳质物可能来源于CH 4 裂解反应[式(5)]、CO歧化(Boudouard)反应[式(6)]、CO加氢反应[式(7)]和CO 2 加氢反应[式(8)]

1.3

计算方法

总吉布斯自由能最小化被认为是定义任何反应系统平衡状态的合适工具。多组分混合物在任何温度和压力条件下的热力学平衡分析都涉及系统吉布斯自由能的计算和最小化。结果是在特定的工艺条件下获得产物的平衡分布,这提供了不同工艺变量之间的直接关系。这反过来又有助于找到最佳的操作条件。通过将所有组分的化学势相加,计算出系统的总吉布斯自由能[式(9) ]。
式中, n i 是组分 i 的物质的量; μ i 是组分 i 的化学势; N 是系统中组分的总数。物质 i 的化学势可定义为式(10)
式中, 是物质 i 的标准生成吉布斯自由能; R 是摩尔气体常数; T 是热力学温度; y i 是组分 i 的摩尔分数; 是使用相关状态方程计算的组分 i 的逸度系数; p 是操作压力; p 0 是标准压力。将 式(10)代入式(9)中,得到吉布斯自由能的表达式为式(11)。

考虑到固体碳(C s )物种生成的吉布斯能表达式,式(11)可以修改为式(12)。

式(12)的最小化可以通过至少两种方式来实现。一般常用的方法是利用拉格朗日待定乘数概念方法。该方法要求在系统化学势的方程[式(10)]中引入拉格朗日变量。具有拉格朗日乘子的吉布斯自由能的总表达式可以等于零。这产生了一个由 N 个非线性方程组成的系统,该系统可以使用非线性方程来求解,如式(13)。
式中, λ k 是元素 k 的拉格朗日待定乘数。
物质约束受到系统中使用的不同元素组分的质量平衡的影响,对于BRM的情况将CH 4 、H 2 O、CO 2 、CO、H 2 和C作为反应混合物。质量平衡约束为式(14)。

式中, a ik 是第 k 个元素的原子数; A k k 元素的总质量。

利用Aspen Plus过程模拟对稳态和非稳态过程系统进行建模,Aspen Plus中提供的RGibbs反应器用作基于吉布斯自由能最小化方法预测生产气体的平衡组成。给定温度和压力的平衡态是给定组分的吉布斯自由能通过最小化来确定的,而不需要对系统中最终可能发生的反应作任何指定。在RGibbs反应器中,计算是根据化学计量约束、反应压力和温度来改变每个相中每个物种的物质的量,并找到一个使系统的总吉布斯能最小化的组成。该模块能够模拟单相或多相的多组分平衡,在吉布斯自由能最小条件下,计算反应体系组成。无论反应网络如何,该反应器通过最小化反应体系中每种现有物质的吉布斯自由能来计算反应产物的分布。Aspen RGibbs主要推荐Peng-Robinson状态方程,因为它在温度和压力方面的适用范围最大,Peng-Robinson状态方程适用于处理含有碳氢化合物、水、空气和燃烧气体的系统。鉴于上述事实,在逸度计算中使用Peng-Robinson状态方程。为了进行计算,在特定的温度和压力下,确定了适当数量的反应物和生成物种类,然后进行最小化以获得平衡生成物和反应物组成。

平衡组成CH 4 转化率、CO 2 转化率和H 2 /CO摩尔比根据式(15)~式(17)定义。

平衡积炭量为以进料条件计算所得反应达到平衡时生成碳物质的量。

1.4

模型验证

为了验证模型的准确性,进行了充分的尝试,通过其他研究者的BRM实验及计算结果来验证计算模型。本研究的模型结果与Jabbour和Cunha等发表的实验BRM结果在CH 4 转化率方面的比较,如图2(a)所示。模拟计算条件为1bar(1bar=0.1MPa)压力、600~850℃下进料CH 4 ∶CO 2 ∶H 2 O=3∶1∶2,他们的结果与本文的计算结果变化趋势非常一致,特别是在高温状态下。由于BRM是一种高吸热反应过程,需要高温才能使该反应过程达到可行的状态。可以看出,反应温度是导致动力学限制的控制因素,在较低的温度下,达到平衡的转化率要低得多。Noureldin等在SRM的研究中报告了炭沉积趋势,如图2(b)所示,也与本文的模 型结果非常一致。因此,使用本文的计算模型进行的平衡分析具有可信度,从而能更好地了解BRM系统。

图2 模型验证结果


2

联合重整反应热力学模拟计算

2.1

对主、副反应吉布斯自由能计算与分析

BRM反应中涉及的主、副反应的吉布斯自由能随反应温度的变化如图3所示,各反应吉布斯自由能在常压条件下计算。图3(a)为主反应过程,分别为DRM[式(1)]、SRM[式(2)]以及式(1)和式(2)结合得到的BRM的吉布斯自由能随温度的变化曲线。它们均为气体体积增大的强吸热反应,Δ G 随着温度的升高而降低,高温有利于反应朝正方向进行;从图3中可以看出,至少640℃以上,DRM、SRM以及BRM的吉布斯自由能为负。因此,总体来说,BRM主反应在高温下热力学是有利的。






请到「今天看啥」查看全文