本发明专利技术涉及一种克劳斯硫磺回收过程中反应炉动态建模方法。本发明专利技术在原子物料守恒、Gibbs自由能最小原则及能量守恒计算的基础上，首次提出通过将热量损失计入能量守恒方程进行迭代计算判断反应炉内的温度是否合理，进而计算反应炉内的组分，为克劳斯硫磺回收过程中反应炉动态建模提出了一种新的方法。在计算过程中，首先假定克劳斯硫磺回收过程中反应炉内的温度及各组分的量，通过物料平衡方程、Gibbs自由能最小原则计算反应炉内组分含量，然后通过能量平衡方程判断假定的温度是否合理。本发明专利技术初次提出一种克劳斯硫磺回收过程中反应炉动态建模方法，为克劳斯硫磺回收过程的动态建模计算提供了一条新的途径。
【技术实现步骤摘要】
本专利技术属于化工
，涉及一种克劳斯硫磺回收过程中反应炉动态建模方法。
技术介绍
随着全球含硫原油和天然气资源的大量开发，工厂和生活过程中产生的酸性气体已经对生态环境产生了恶劣的影响，因此，硫磺回收装置已成为大型天然气净化厂、炼油厂、石油化工厂加工含硫天然气、含硫原油时不可缺少的配套装置。目前，通过使用克劳斯法回收硫磺的工艺已成为天然气或炼厂气加工的一个重要组成部分。近年来，克劳斯硫磺回收工艺的研究主要集中在催化过程中催化剂的研究、吸收剂的选择、工艺参数的简化计算及稳态计算模型的研究（包含反应炉数学模型的建立、转化器数学模型的建立等）等。然而，实际生产过程是一个动态过程，反应器内的温度、压力、组成等均受到进料、出料及操作等因素的影响，稳态模型不足以解决实际生产问题，因此开发一种克劳斯硫磺回收过程中动态模型显得尤为必要。
技术实现思路
本专利技术针对现有技术的不足，提出了一种新的克劳斯硫磺回收过程中反应炉动态建模方法。克劳斯硫磺回收过程中参与的酸性气体中除含有H2S,还含有CO2，N2和H2O等，来自炼厂气、焦炉气或水煤气脱硫装置中的酸性气体还可能含有NH3和HCN等。因此，反应炉内实际发生的反应非常复杂，但主要反应有以下几种：2H2S + 3O2 →2SO2 + 2H2O2H2S + SO2 → SX + 2H2OCOS + H2S → CS2 + H2OH2S + CO2 → COS + H2OH2 + CO2 → CO + H2O只需确定反应炉内SX的形态，就可以确定反应炉内的主要热反应，而在气相中的硫蒸汽主要以S2，S6，S8的形式存在，并有如下平衡关系：3S2 → S64S2 → S84S6 → 3S8一般而言，大于1000℃时，气相中主要是S2，小于450℃时，气相中主要是S6，S8。由于反应炉的出炉气温度一般在1000℃以上，所以在反应炉出炉过程气中，硫蒸汽中主要是S2。在克劳斯硫磺回收过程的反应炉动态建模过程中，主要用到的方程包含：物料平衡计算，Gibbs自由能最小的原则及能量平衡计算，通过以下四个具体步骤完成反应炉的动态建模。步骤（1），建立从某一个初始时刻t到下一个时刻t+Δt之间的动态建模过程中的物料平衡方程，具体是：已知t时刻反应器内各组分的物质的量、进料组分及流量，在反应炉中，先将t时刻反应器内的物料与Δt时间段内的进料物质混合，视为此轮计算的初始反应物；建立物料平衡方程，其模型表达式为： 式中，M为此轮计算初始反应物的量，mol；为反应炉内原有物料的累积量，单位为mol；为进料组分i的摩尔流量，单位为mol/h；为每轮计算所需的时间，单位为h。反应炉内Δt内进入的物料、t内积累的物料与t+Δt时刻反应开始后物料尊循原子守恒，其模型如下：硫：BS =+ 氢：BH =+氧：BO =+碳 ：BC =+式中，为t+Δt时刻，反应开始后j元素的原子总数；为t时刻反应炉内累积的元素j的原子总数；为Δt时刻内进入反应炉内的物料中元素j的原子总数；j分别表示硫、氢、氧、碳元素中的一种；步骤（2），设克劳斯硫磺回收过程中反应炉温度的初值为T，由于系统Gibbs自由能为温度和组成的函数，其总Gibbs自由能可表示如下：其中，G为混合物系统的总Gibbs自由能，kJ；为（）时间段内反应器中所有组分i累积的总量，mol；为组分i的Gibbs自由能，KJ/mol，其计算表达式为： 式中，是标态下组分i的Gibbs自由能，kJ/mol；fi为组分i在当前温度压力下对应的逸度，kPa；为组分i在标准状态下的逸度，kPa；R为热力学常数，J/(mol·K) ；T表示温度，K；系统的总Gibbs自由能可以表示为：由于在反应达到平衡时，整个系统的Gibbs自由能会取得其最小值；结合步骤（1）的原子数守恒公式，可以将该时刻反应达到平衡时整个系统的组成和相分率的求解问题转化为求解如下最优化问题: 采用拉格朗日因子法解决此最优化问题：首先，将原子数守恒条件乘以拉格朗日因子，然后加和到Gibbs自由能表达式上；在设定的P、T下，对系统Gibbs自由能表达式求的偏导，其偏导为0，则有： + ( ) = 0 + ( ) = 0·· + = 0联立原子守恒方程与Gibbs方程，即可通过牛顿迭代求解出当前温度压力下对应的物料组成。步骤（3），根据能量平衡方程,建立反应炉内的能量平衡模型，计算在时间段反应器内物料的总焓值，其表达式如下： 式中，(MH)t 为反应器内累积物料的总焓值，为t时刻反应器内累积物料的摩尔焓，KJ/mol；为反应器内每轮计算所需要的时间，h；为进料物流的摩尔焓（可以是多股进料），KJ/mol；为外界提供反应炉的热量，kJ/h；为热损失的量，kJ/h。忽略热辐射的影响，在反应炉内的热损失包含三个部分，分别为：反应炉内部气体对流传热、保温层及器壁传热、空气侧对流。假定保温层内壁温度为，计算热损失的量。在反应炉内部气体与保温层内壁的对流换热属于内强制对流换热，其对流传热系数的计算公式为： 式中，为炉内气体对流传热系数，；为气体在平均温度下的导热系数，W/(m·℃)；为气体在平均温度下的运动粘度，/s；为气体在平均温度下的普朗特准数；为当量直径，m；为气体温度与保温层内壁温度的修正系数；为受热面的相对长度修正系数；为辐射换热修正系数；为气体流速，m/s。平均温度为反应炉内部气体温度与保温层内壁温度的平均值。在内外介质温度保持不变的情况下，反应炉通过壁面向外传热损失的热量恒定不变，因此，单位面积内通过对流传热的热损失q为： 假设保温层与器壁在同一时刻始终在热平衡状态，稳定传热时，保温层壁内外之间的热通量相等，由传热边界条件得： 式中，Tw为保温层外壁温度，℃；h为平均对流传热系数，；可通过上式根据已知的保温层外壁温度Tw ，求出器壁外侧温度Two；在确定空气温度为t的情况下，通过上述传热边界方程可推测计算空气在自然对流状态下的对流传热系数：根据值是否在1～10之间，确定假设的保温层内壁温度是否合理。若的值超过这个取值范围，则重新假定保温层内壁温度，重新开始计算；如果合理，则继续计算反应炉内总的传热系数： 式中，分别为内外侧壁表面与流体的对流传热系数，；为构成反应炉壁各层材料的导热系数，；反应炉壁各种材料的厚度，m。通过总的传热系数，有效传热面积A，计算反应炉内热量的总损失： 步骤（4），由于物料的摩尔焓为温度、压力及组成的函数，其表达式为： 通过已知的压力P、步骤（2）中本文档来自技高网...
【技术保护点】
一种克劳斯硫磺回收过程中反应炉动态建模方法，其特征在于该方法包括以下步骤：步骤（1），建立从某一个初始时刻t到下一个时刻t+Δt之间的动态建模过程中的物料平衡方程，具体是：已知t时刻反应器内各组分的物质的量、进料组分及流量，在反应炉中，先将t时刻反应器内的物料与Δt时间段内的进料物质混合，视为此轮计算的初始反应物；建立物料平衡方程，其模型表达式为：式中，M为此轮计算初始反应物的量，mol；为反应炉内原有物料的累积量，单位为mol；为进料组分i的摩尔流量，单位为mol/h；为每轮计算所需的时间，单位为h；反应炉内Δt内进入的物料、t内积累的物料与t+Δt时刻反应开始后物料尊循原子守恒，其模型如下：硫：BS =+氢：BH =+氧：BO =+碳 ：BC =+式中，为t+Δt时刻，反应开始后j元素的原子总数；为t时刻反应炉内累积的元素j的原子总数；为Δt时刻内进入反应炉内的物料中元素j的原子总数；j分别表示硫、氢、氧、碳元素中的一种；步骤（2），设克劳斯硫磺回收过程中反应炉温度的初值为T，由于系统Gibbs自由能为温度和组成的函数，其总Gibbs自由能可表示如下：其中，G为混合物系统的总Gibbs自由能，kJ；为（）时间段内反应器中所有组分i累积的总量，mol；为组分i的Gibbs自由能，KJ/mol，其计算表达式为：式中，是标态下组分i的Gibbs自由能，kJ/mol；fi为组分i在当前温度压力下对应的逸度，kPa；为组分i在标准状态下的逸度，kPa；R为热力学常数，J/(mol·K) ；T表示温度，K；系统的总Gibbs自由能可以表示为：由于在反应达到平衡时，整个系统的Gibbs自由能会取得其最小值；结合步骤（1）的原子数守恒公式，可以将该时刻反应达到平衡时整个系统的组成和相分率的求解问题转化为求解如下最优化问题: 采用拉格朗日因子法解决此最优化问题：首先，将原子数守恒条件乘以拉格朗日因子，然后加和到Gibbs自由能表达式上；在设定的P、T下，对系统Gibbs自由能表达式求的偏导，其偏导为0，则有：+ () = 0+ () = 0··+= 0联立原子守恒方程与Gibbs方程，即可通过牛顿迭代求解出当前温度压力下对应的物料组成；步骤（3），根据能量平衡方程,建立反应炉内的能量平衡模型，计算在时间段反应器内物料的总焓值，其表达式如下：式中，(MH)t 为反应器内累积物料的总焓值，为t时刻反应器内累积物料的摩尔焓，KJ/mol；为反应器内每轮计算所需要的时间，h；为进料物流的摩尔焓（可以是多股进料），KJ/mol；为外界提供反应炉的热量，kJ/h；为热损失的量，kJ/h；忽略热辐射的影响，在反应炉内的热损失包含三个部分，分别为：反应炉内部气体对流传热、保温层及器壁传热、空气侧对流；假定保温层内壁温度为，计算热损失的量；在反应炉内部气体与保温层内壁的对流换热属于内强制对流换热，其对流传热系数的计算公式为：式中，为炉内气体对流传热系数，；为气体在平均温度下的导热系数，W/(m·℃)；为气体在平均温度下的运动粘度，/s；为气体在平均温度下的普朗特准数；为当量直径，m；为气体温度与保温层内壁温度的修正系数；为受热面的相对长度修正系数；为辐射换热修正系数；为气体流速，m/s；平均温度为反应炉内部气体温度与保温层内壁温度的平均值；在内外介质温度保持不变的情况下，反应炉通过壁面向外传热损失的热量恒定不变，因此，单位面积通过对流传热的热损失q为：假设保温层与器壁在同一时刻始终在热平衡状态，稳定传热时，保温层壁内外之间的热通量相等，由传热边界方程得：式中，Tw为保温层外壁温度，℃；h为平均对流传热系数，；可通过上式根据已知的保温层外壁温度Tw求出器壁外侧温度Two；在确定空气温度为t的情况下，通过上述传热边界方程可推测计算空气在自然对流状态下的对流传热系数：根据值是否在1～10之间，确定假设的保温层内壁温度是否合理；若的值超过这个取值范围，则重新假定保温层内壁温度，重新开始计算；如果合理，则继续计算反应炉内总的传热系数：式中，分别为内外侧壁表面与流体的对流传热系数，；为构成反应炉壁各层材料的导热系数，；反应炉壁各种材料的厚度，m；通过总的传热系数，有效传热面积A，计算反应炉内热量的总损失：步骤（4），由于物料的摩尔焓为温度、压力及组成的函数，其表达式为：通过已知的压力P、步骤（2）中假定的温度T和求解的组成Zi，计算反应后反应器内物料的摩尔焓，那物料的总焓值表达式为：H根据能量平衡方程,反应前后的总能量应该相等，建立反应炉内的能量平衡模型，其表达式如下：= |‑|若，则返回步骤（2），估算新的温度T，重新开始计算；若，则输出计算结果。...
【技术特征摘要】
1. 一种克劳斯硫磺回收过程中反应炉动态建模方法，其特征在于该方法包括以下步骤：
步骤（1），建立从某一个初始时刻t到下一个时刻t+Δt之间的动态建模过程中的物料平衡方程，具体是：
已知t时刻反应器内各组分的物质的量、进料组分及流量，在反应炉中，先将t时刻反应器内的物料与Δt时间段内的进料物质混合，视为此轮计算的初始反应物；建立物料平衡方程，其模型表达式为：
式中，M为此轮计算初始反应物的量，mol；为反应炉内原有物料的累积量，单位为mol；为进料组分i的摩尔流量，单位为mol/h；为每轮计算所需的时间，单位为h；
反应炉内Δt内进入的物料、t内积累的物料与t+Δt时刻反应开始后物料尊循原子守恒，其模型如下：
硫：BS =+ 氢：BH =+氧：BO =+碳 ：BC =+式中，为t+Δt时刻，反应开始后j元素的原子总数；为t时刻反应炉内累积的元素j的原子总数；为Δt时刻内进入反应炉内的物料中元素j的原子总数；j分别表示硫、氢、氧、碳元素中的一种；
步骤（2），设克劳斯硫磺回收过程中反应炉温度的初值为T，由于系统Gibbs自由能为温度和组成的函数，其总Gibbs自由能可表示如下：
其中，G为混合物系统的总Gibbs自由能，kJ；为（）时间段内反应器中所有组分i累积的总量，mol；为组分i的Gibbs自由能，KJ/mol，其计算表达式为：
式中，是标态下组分i的Gibbs自由能，kJ/mol；fi为组分i在当前温度压力下对应的逸度，kPa；为组分i在标准状态下的逸度，kPa；R为热力学常数，J/(mol·K) ；T表示温度，K；系统的总Gibbs自由能可以表示为：
由于在反应达到平衡时，整个系统的Gibbs自由能会取得其最小值；结合步骤（1）的原子数守恒公式，可以将该时刻反应达到平衡时整个系统的组成和相分率的求解问题转化为求解如下最优化问题:
采用拉格朗日因子法解决此最优化问题：首先，将原子数守恒条件乘以拉格朗日因子，然后加和到Gibbs自由能表达式上；在设定的P、T下，对系统Gibbs自由能表达式求的偏导，其偏导为0，则有：
+ ( ) = 0
+ ( ) = 0
·
·
+ = 0
联立原子守恒方程与Gibbs方程，即可通过牛顿迭代求解出当前温度压力下对应的物料组成；
步骤（3），根据能量平衡方程,建立反应炉内的能量平衡模型，计算在时间段反应器内物料的总焓值，其表达式如下：
式中，(MH)t 为反应器内累积物料的总焓值，为t...
【专利技术属性】
技术研发人员：葛铭，魏江，葛文锋，郑小青，郑松，李小爽，
申请(专利权)人：杭州电子科技大学，
类型：发明
国别省市：浙江;33
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。