(19)中华 人民共和国 国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202111486621.1 (22)申请日 2021.12.07 (71)申请人 浙江中控技 术股份有限公司 地址 310053 浙江省杭州市滨江区六和路 309号 (72)发明人 蒋益波　王宽心　江泉达　韩文巧　 黄小婷　吴哲峰　郭廉洁　陈镭　 (74)专利代理 机构 北京易捷胜知识产权代理事 务所(普通 合伙) 11613 代理人 李会娟 (51)Int.Cl. G06F 30/20(2020.01) C01B 3/48(2006.01) C01B 3/38(2006.01) C01B 3/56(2006.01)G06F 119/08(2020.01) (54)发明名称 炼厂多制氢装置原料气供给方案优化方法 和系统 (57)摘要 本申请涉及一种炼厂多制氢装置原料气供 给方案优化方法， 该方法包括： S10、 将炼厂原料 气数据输入到预先建立的生产过程模 型中， 得到 每种原料气在每个制氢装置的产氢比和单位产 氢成本； S2 0、 根据单位产氢成本的大小对原料气 与制氢装置组合进行排序， 得到组合队列； S30、 按照单位产氢成本从低到高的顺序， 基于原料气 数据产氢比和每个制氢装置的产氢上下限约束， 计算各组合的原料气实际供应量和实际产氢量； S40、 按照单位产氢成本从低到高的顺序选取组 合， 直至选取组合的总产氢量满足预设产氢量要 求； S50、 将选 取的组合和相应的原料气实际供应 量作为方案优化结果。 本优化方法不增加装置和 成本， 可对运行参数进行定量 化调节。 权利要求书4页 说明书15页 附图4页 CN 114297829 A 2022.04.08 CN 114297829 A 1.一种炼厂多制氢装置原料气供给方案优化方法， 其特征在于， 各制氢装置以天然气 和/或炼厂干气为原料气， 采用烃类水蒸汽转化法造气， 并经PSA法净化提纯制取氢气， 该方 法包括： S10、 获取炼厂原料气数据， 将所述炼厂原料气数据输入到预先建立的生产过程模型 中， 得到每种原料气在每 个制氢装置的产氢比和单位产氢成本； S20、 根据所述单位产氢成本， 对原料气与制氢装置组合进行排序， 得到组合队列； S30、 按照单位产氢成本从低到高的顺序， 基于所述原料气数据、 所述产氢比和每个制 氢装置的产氢上下限约束， 获取所述组合队列中每个组合的原料气实际供应量和实际产氢 量； S40、 根据每个组合的实际产 氢量， 按照单位产氢成本从低到高的顺序从所述组合队列 中选取组合， 直至 选取组合的总产氢量满足预设的产氢量要求； S50、 将选取的组合和相应的原料气实际供应量作为炼厂多制氢装置原料气供给方案 优化结果。 2.根据权利要求1所述的炼厂多制氢装置原料气供给方案优化方法， 其特征在于， 所述 炼厂原料气数据包括潜在原料气种类、 流体压缩前压力、 原料气流量的可供给上限、 潜在原 料气组成、 原料气流 量约束。 3.根据权利要求2所述的炼厂多制氢装置原料气供给方案优化方法， 其特征在于， 所述 S10之前， 还 包括， 建立所述 生产过程模型； 具体包括： S01、 针对每种原料气， 获取各制氢装置在每种原料气下的相关运行历史数据， 所述相 关运行历史数据包括 LIMS分析数据、 原料组成热力学 数据、 历史运行 数据； 所述LIMS分析数据包括原料气组成分析 数据、 中变气组成分析 数据； 所述原料组成热力学 数据包括 生成焓、 吉布斯自由能； 所述历史运行数据包括原料气流量、 装置产氢量、 除盐水流量、 反应消耗水量、 除盐水 温度、 配汽量、 配汽冷凝温度、 转化炉温度、 中变器温度、 流体压缩后压力、 流体压缩前的温 度、 能源价格； 所述能源价格包括电价、 天然气价格、 除盐水价格； S02、 基于所述转化炉温度、 所述中变器温度、 所述生成焓及所述吉布斯自由能计算转 化炉中烃类转化为一氧化碳与氢气的平衡常数K1、 烃类转化为二氧化碳与氢气的平衡常数 K2、 中变器中一氧化 碳转化为二氧化碳的平衡常数 K3； α1CnHm+β1H2O＝＝γ1CO+δ1H2 α2CnHm+β2H2O＝＝γ2CO2+δ2H2 α3CO+β3H2O＝＝γ3CO2+δ3H2 其中： α 、 β 、 γ、 δ 分别为 化学方程式的化学计量数； 每个平衡常数在温度为T时的计算式如下： 权　利　要　求　书 1/4 页 2 CN 114297829 A 2其中： KT、 Kθ分别表示温度为T与Tθ下的平衡常数； G表示 吉布斯自由能， 表示反应吉 布斯自由能变； H表示焓， 表示反应焓变； R表示气体体积常数， T表示反应器中的温度， Tθ代表标准状况 下的温度； S03、 基于所述原料气流量、 所述原料气组成分析数据、 所述配汽量、 所述平衡常数K1、 K2、 K3及反应过程物质守恒定律， 计算得到中变气组成理论数据； S04、 将所述中变气组成理论数据与中变气组成分析数据比对， 采用偏最小二乘回归拟 合制氢装置转化反应与中变反应的催化剂性能系数； 基于所述催化剂性能系数确定最终的 装置产氢比； S05、 基于所述原料气组成分析数据、 所述中变气组成分析数据计算反应焓变， 计算式 如下： ΔH＝ΔH1+ΔH2 其中： 分别代表反应过程中消耗的烃类与蒸汽的物质的量， nCO分别代 表反应过程中生成的二氧化碳、 一氧化碳的物质的量； 分别代表烃类、 蒸汽、 二氧化碳与一氧化碳的标准生成焓； C代表恒压热容； T、 Tθ分别代表组 成温度与标 况温度； Δ H1、 ΔH2、 ΔH分别代表标况下反应焓变、 反应的升温吸热与反应焓变； S06、 基于所述除盐水温度、 所述除盐水流量、 所述反应消耗水量、 所述配汽冷凝温度计 算反应热损耗； Q＝[(T2‑T1)×CH+H]×(F2‑F1)×A 其中： Q代表热损耗， T2和T1分别代表配汽冷凝温度与除盐水温度， CH代表水的热容， H代 表汽化热， F2和F1分别代表除盐水流 量、 反应消耗水量， A 代表配汽冷凝率； S07、 依据 所述中变气组成分析数据和所述装置产氢量计算得到解析气热值， 依据 所述 反应焓变、 所述反应热损耗、 所述解析气热值以及能源价格， 计算得到各制氢装置在每种原 料气下的单位产氢的热能成本， 依据所述装置产氢量与所述原料气流量计算得到各制氢装 置在每种原料气下的产氢比： S08、 基于模型的输入数据以及所述单位产氢的热能成本和所述产氢比， 通过以下公式 计算得到单位产氢成本， 将所述单位产氢成本和所述产氢比作为模型输出；权　利　要　求　书 2/4 页 3 CN 114297829 A 3