删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

耦合变压吸附简化模型的提纯回用氢网络协调优化

本站小编 Free考研考试/2020-04-15

邓春 1 , 周业扬 1 , 江苇 1 , 冯霄 2
1.中国石油大学(北京) 化学工程学院, 重质油国家重点实验室, 北京 102249;
2.西安交通大学 化学工程与技术学院, 西安 710049

收稿日期: 2015-08-25
基金项目: 国家“九七三”重点基础研究资助项目(2012CB720500);国家自然科学基金资助项目(21576287,U1162121);中国博士后科学基金资助项目(2015M570215);中国石油大学(北京)科研基金资助项目(2462015YQ0305,2462015BJB02)
作者简介: 邓春(1984-),男,副教授。E-mail:chundeng@cup.edu.cn


摘要:采用变压吸附(PSA)装置提纯含氢流股并回用至加氢过程,可以缓解炼油厂氢气亏缺的现状。该文采用了变压吸附简化模型,构建了提纯回用氢网络模型,以氢公用工程用量为目标函数,建立优化的数学模型,使用商业优化软件GAMS(general algebraic modeling system)平台建模,用DICOPT(discrete and continuous optimizer)作为求解器。案例研究结果表明:变压吸附装置存在最优入口氢气纯度。当吸附和解吸压力比增大、吸附选择性减小时,变压吸附的氢气回收率增大,氢公用工程用量减小。从氢网络优化的角度来说,一味地增加变压吸附装置入口流股氢气纯度以提高回收率的手段并不可取。
关键词: 过程系统 氢网络 优化 数学规划法 变压吸附(PSA) 提纯回用
Coordinated optimization of hydrogen networks with purification reuse coupled with a shortcut model of the pressure swing adsorption
DENG Chun1, ZHOU Yeyang1, JIANG Wei1, FENG Xiao2
1.State Key Laboratory of Heavy Oil Processing, College of Chemical Engineering, China University of Petroleum-Beijing, Beijing 102249, China;
2.School of Chemical Engineering and Technology, Xi'an Jiaotong University, Xi'an 710049, China


Abstract:Pressure swing adsorption (PSA) is used to purify hydrogen streams with the upgraded streams reused by hydroprocessing units to increase the hydrogen supply in refineries. This study used a PSA shortcut model to construct a hydrogen network design with purification reuse. The mathematical model used the minimum hydrogen utility as the objective function with the commercial optimization software GAMS (general algebraic modeling system) to solve the problem with DICOPT (discrete and continuous optimizer) as the solver. The results of a case study show that the PSA device has the optimal inlet hydrogen purity. As the adsorption and desorption pressure ratio increases and the adsorbent selectivity decreases, the PSA hydrogen recovery increases and the hydrogen utility decreases. The hydrogen network optimization shows that designers should not blindly increase the PSA inlet hydrogen purity to enhance the hydrogen recovery.
Key words: process systemhydrogen networkoptimizationmathematical programmingpressure swing adsorption (PSA)purification reuse
目前,重质劣质原油的加工量在日益增加,而环保法规对清洁油品的硫含量的要求却越来越低,这必然导致炼油厂增加加氢过程的比例,使得炼油厂对氢气的亏缺正成为普遍问题,氢气在炼油厂成为了一种宝贵资源。工业中,天然气重整制氢应用最为广泛,然而该过程是典型的高耗能过程。为了节约成本,炼油厂愈发重视炼厂气的回收利用,广泛采用变压吸附装置提纯氢气流股并回用至加氢过程,从而降低制氢装置的负荷。此外,氢系统集成是过程集成的分支,能够有效地优化氢气网络、提高氢气利用率、节约氢气消耗,实现节能降耗的目的。
通常,氢系统集成与优化的方法主要有2类: 夹点法和数学规划法。代表性的夹点法有氢剩余图[1]、 物料回收夹点图[2]、 气体级联图[3]、 复合表算法和改进极限复合曲线[4]以及改进问题表法[5]等。数学规划法能够解决多约束大规模复杂氢网络集成优化的问题,一直是研究的热点。Hallale和Liu[6]首次提出了超结构优化氢气分配网络的方法。随后,很多****提出了新的优化模型和求解方法,例如多阶段优化模型[7]、 状态-空间超结构模型[8]、 集成闪蒸计算模型[9]、 不同情境的对比分析模型[10]、 中间浓度氢气管网模型[11]、 集成加氢过程简化模型[12]等。
变压吸附(pressure swing adsorption,PSA)提纯氢气过程的性能参数包括吸附材料和物系的性质、操作温度和压力、装置的处理能力、入口氢气流股的浓度、提纯产品气的纯度和尾气的纯度、氢气回收率等。目前,变压吸附提纯氢气过程的模型一般简化为固定提纯产品氢气的纯度[6]和固定氢气回收率[6]。显然,目前的研究缺乏这些性能参数之间关联性的考量,例如提纯产品氢气的纯度与氢气回收率之间的关系、吸附温度和压力、吸附材料和物系的性质等也会影响吸附选择性和回收率等。Liu和Zhang[13]引入了变压吸附的简化模型,并探讨了回收率与吸附选择性、吸附操作压力和入口氢气流股的浓度之间的关系。然而,文中并未深入分析变压吸附过程性能参数(例如,入口氢气流股的浓度、氢气回收率、吸附材料和物系的性质、操作温度和压力等)对氢网络优化目标(例如,最小的氢气公用工程用量、最优的提纯产品气流率等)的影响。
本文旨在构建耦合变压吸附简化模型的提纯回用氢网络优化模型,考察变压吸附氢气提纯过程性能参数与氢网络目标值之间的关联性,实现两者的协调优化。
1 问题描述给定一个炼油厂氢网络,它包括氢公用工程、过程氢源、氢阱、燃料系统、压缩机和提纯器(图 1)。氢公用工程(u∈NHU)主要来自制氢厂。每个过程氢源(s∈NSR)都有规定的最大流率(FSRs),c组分的浓度为(ys,cout(c∈NC)),压力为PSRs。每个氢阱(k∈NSK)都有它自身的流率要求,最小流率要求为FKkin。每个氢阱都对入口的气体流股有最小的氢纯度要求以及压力要求(PKkin)。需要设置压缩机 (i∈NI)以提升氢源的压力满足氢阱的需求。引入提纯器(p∈NP)变压吸附装置(PSA),提升部分过程氢源的纯度,从而减少氢公用工程的用量。本文旨在实现变压吸附过程和氢网络的协调优化。
图 1 提纯回用氢网络模型
图选项





2 变压吸附模型对于变压吸附提纯氢气的装置(PSA),氢气回收率并不是常数,它与入口气体流股的纯度、产品气体的纯度、吸附材料和物系的性质、操作温度和压力等因素有关。Lacava等[14]提出了PSA系统的设计原则。本文采用Knaebel[15]提出的简化模型,文中提出氢气回收率可以用式(1)表示,即
${{R}_{\text{PSA}}}=\left( 1-\theta \right)1-\frac{1}{\frac{{{P}_{\text{H}}}}{{{P}_{\text{L}}}}y_{p,c}^{\text{in}}}.$ (1)
其中: RPSA代表PSA的氢气回收率; yp,cin代表PSA的入口氢气纯度; PH是指变压吸附循环的吸附压力(高压); PL是指变压吸附循环的解吸压力(低压),均为绝压; θ是吸附剂选择性,其取值范围在0到1之间,与吸附材料和物系的性质以及操作温度有关。
由式(1)可知,当吸附压力与解吸压力的相对比值增大时,氢气回收率增大。当PSA的入口气体流股的氢气纯度增大时,氢气回收率也增大。吸附选择性θ值越小,氢气回收率越大。
相对于严格模型而言,变压吸附简化模型不需要严格的数学模型和大量的数据(例如吸附循环的次数、操作温度等),也能够预测回收率和产品量。这些正是与氢网络优化相关联的变量。因此,本文利用简化模型代替严格模型,用于提纯回用氢网络优化。
3 提纯回用氢网络模型提纯回用氢网络模型基于邓春等提出的氢气网络超结构模型[16]
1) 与第u个氢气公用工程相关的公式
第u个氢气公用工程分离节点的流量平衡为:
$\begin{align} & \text{FH}{{\text{U}}_{u}}=\sum\limits_{k\in \text{NSK}}{\text{FU}{{\text{K}}_{u,k}}}+\sum\limits_{i\in \text{NI}}{\text{FU}{{\text{I}}_{u,i}}}, \\ & \forall u\in \text{NHU}\text{.} \\ \end{align}$ (2)
其中: FHUu是氢气公用工程的流量; FUKu,k是氢气公用工程分配到氢阱的流量; FUIu,i是氢气公用工程分配到压缩机的流量。
第u个氢气公用工程可以提供的最大流量约束为:
$\text{FH}{{\text{U}}_{u}}\le \text{FHU}_{u}^{\text{UB}}.$ (3)
2) 与第s个过程氢源相关的公式
第s个过程氢源分离节点的流量平衡为:
$\begin{align} & \text{FS}{{\text{R}}_{s}}=\sum\limits_{k\in \text{NSK}}{\text{FS}{{\text{K}}_{s,k}}}+\sum\limits_{i\in \text{NI}}{\text{FS}{{\text{I}}_{s,i}}+} \\ & \sum\limits_{p\in \text{NP}}{\text{FS}{{\text{P}}_{s,p}}}+\sum\limits_{f\in \text{NF}}{\text{FS}{{\text{F}}_{s,f}},}\forall s\in \text{NSR}\text{.} \\ \end{align}$ (4)
第s个过程氢源可以提供的最大流量约束为:
$\text{FS}{{\text{R}}_{s}}\le \text{FSR}_{s}^{\text{UB}}$ (5)
3) 与第i个压缩机相关的公式
第i个压缩机混合节点前的流量平衡及组分质量平衡为:
$\begin{align} & \text{FI}_{i}^{\text{in}}=\sum\limits_{u\in \text{NHU}}{\text{FU}{{\text{I}}_{u,i}}}+\sum\limits_{s\in \text{NSR}}{\text{FS}{{\text{I}}_{s,i}}}+ \\ & \sum\limits_{p\in \text{NP}}{\text{FPI}_{p,i}^{\text{prod}}}+\sum\limits_{\begin{smallmatrix} {i}'e i \\ {i}'\in \text{NI} \end{smallmatrix}}{\text{FI}{{\text{I}}_{{i}',i}},}\forall i\in \text{NI}\text{.} \\ & \text{FI}_{i}^{\text{in}}\cdot y_{i,c}^{\text{in}}=\sum\limits_{u\in \text{NHU}}{\text{FU}{{\text{I}}_{u,i}}\cdot y_{u,c}^{\text{out}}+\sum\limits_{s\in \text{NSR}}{\text{FS}{{\text{I}}_{s,i}}\cdot y_{s,c}^{\text{out}}}}+ \\ \end{align}$ (6)
$\sum\limits_{\begin{smallmatrix} {i}'e i \\ {i}'\in \text{NI} \end{smallmatrix}}{\text{FI}{{\text{I}}_{{i}',i}}\cdot y_{{i}',c}^{\text{out}}},\forall i\in \text{NI,}\forall c\in \text{NC}\text{.}$ (7)
其中:FIiin是压缩机入口流量; yi,cin是压缩机入口c组分的纯度; yu,cout是氢气公用工程出口c组分的纯度; ys,cout是过程氢源c组分的纯度。
针对第i个压缩机的出口和入口,分别作出流量平衡方程和c组分的质量平衡方程为:
$\text{FI}_{i}^{\text{out}}=\text{F}{{\text{I}}^{\text{in}}},\forall i\in \text{NI}\text{.}$ (8)
$y_{i,c}^{\text{out}}=y_{i,c}^{\text{in}},\forall i\in \text{NI,}\forall c\in \text{NC}\text{.}$ (9)
第i个压缩机之后分离节点的流量平衡为:
$\begin{align} & \text{FI}_{i}^{\text{out}}=\sum\limits_{k\in \text{NSK}}{\text{FI}{{\text{K}}_{i,k}}}+\sum\limits_{p\in \text{NP}}{\text{FI}{{\text{P}}_{i,p}}}+ \\ & \sum\limits_{\begin{smallmatrix} {i}'e i \\ {i}'\in \text{NI} \end{smallmatrix}}{\text{FI}{{\text{I}}_{i,{i}'}}},\forall i\in \text{NI}\text{.} \\ \end{align}$ (10)
压缩机最大入口流量约束为:
$\text{FI}_{i}^{\text{in}}\le \text{FI}_{i}^{\text{in,UB}},\forall i\in \text{N}{{\text{I}}_{\text{exist}}}.$ (11)
其中NIexist表示现有的压缩机的集合。
4) 与第p个提纯器相关的公式
p个提纯器之前混合节点的流量平衡和组分质量平衡为:
$\begin{align} & \text{FP}_{p}^{\text{in}}=\sum\limits_{s\in \text{NSR}}{\text{FS}{{\text{P}}_{s,p}}}+\sum\limits_{i\in \text{NI}}{\text{FI}{{\text{P}}_{i,p}}}, \\ & \forall p\in \text{NP}\text{.} \\ \end{align}$ (12)
$\begin{align} & \text{FP}_{p}^{\text{in}}\cdot y_{p,c}^{\text{in}}=\sum\limits_{x\in \text{NSR}}{\text{FS}{{\text{P}}_{s,p}}\cdot y_{s,c}^{\text{out}}+\sum\limits_{i\in \text{NI}}{\text{FI}{{\text{P}}_{i,p}}\cdot y_{i,c}^{\text{out}},}} \\ & \forall p\in \text{NP,}\forall c\in \text{NC}\text{.} \\ \end{align}$ (13)
其中: FPpin是提纯器的入口流量; yp,cin是提纯器入口c组分的纯度。
第p个提纯器的流量平衡和组分质量平衡为:
$\text{FP}_{p}^{\text{in}}=\text{FP}_{p}^{\text{prod}}+\text{FP}_{p}^{\text{resd}},\forall p\in \text{NP;}$ (14)
$\begin{align} & \text{FP}_{p}^{\text{in}}\cdot y_{p.c}^{\text{in}}=\text{FP}_{p}^{\text{prod}}\cdot y_{p,c}^{\text{prod}}+\text{FP}_{p}^{\text{read}}\cdot y_{p,c}^{\text{resd}}, \\ & \forall p\in \text{NP,}\forall c\in \text{NC}\text{.} \\ \end{align}$ (15)
其中: FPpprod是提纯器产品气的流量; yp,cprod是提纯器产品气c组分的纯度; FPpresd是提纯器解析气的流量; yp,cresd是提纯器解吸气c组分的纯度。
氢气回收率(RRp,c)可以定义为: 产品部分中c组分的流量与第p个提纯器原料部分c组分流量之比。
氢气回收率定义式如式(16)所示,即
$\text{FP}_{p}^{\text{prod}}\cdot y_{p,c}^{\text{prod}}=\text{R}{{\text{R}}_{p,c}}\cdot \text{FP}_{p}^{\text{in}}\cdot y_{p,c}^{\text{in}}.$ (16)
值得注意的是,与以往的优化模型不同,这里的氢气回收率不是固定值,而是由式(1)计算得出。
第p个提纯器之后分离节点的流量平衡为:
$\begin{align} & \text{FP}_{p}^{\text{prod}}=\sum\limits_{k\in \text{NSK}}{\text{FPK}_{p.k}^{\text{prod}}}+\sum\limits_{k\in \text{NI}}{\text{FPI}_{p.i}^{\text{prod}}}, \\ & \forall p\in \text{NP}\text{.} \\ \end{align}$ (17)
$\text{FP}_{p}^{\text{resd}}=\sum\limits_{f\in \text{NF}}{\text{FPF}_{p,f}^{\text{resd}}},\forall p\in \text{NP}\text{.}$ (18)
其中FPFp,fresd是提纯器的剩余气送往燃料系统的流量。
5) 与第k个氢阱相关的公式
k个氢阱之前混合节点的流量平衡和组分质量平衡为
$\begin{align} & \text{FK}_{k}^{\text{in}}=\sum\limits_{u\in \text{NHU}}{\text{FU}{{\text{K}}_{u,k}}+}\sum\limits_{s\in \text{NSR}}{\text{FS}{{\text{K}}_{s,k}}}+ \\ & \sum\limits_{p\in \text{NP}}{\text{FPK}_{p,k}^{\text{prod}}}+\sum\limits_{i\in \text{NI}}{\text{FI}{{\text{K}}_{i,k}}},\forall k\in \text{NSK}\text{.} \\ \end{align}$ (19)
$\begin{align} & \text{FK}_{k}^{\text{in}}\cdot y_{k,c}^{\text{in}}=\sum\limits_{u\in \text{NHU}}{\text{FU}{{\text{K}}_{u,k}}\cdot y_{u,c}^{\text{out}}}+ \\ & \sum\limits_{s\in \text{NSR}}{\text{FS}{{\text{K}}_{s,k}}\cdot y_{s,c}^{\text{out}}}+\sum\limits_{p\in \text{NP}}{\text{FPK}_{p,k}^{\text{prod}}\cdot y_{p,c}^{\text{prod}}}+ \\ & \sum\limits_{i\in \text{NI}}{\text{FI}{{\text{K}}_{i,k}}\cdot y_{i,c}^{\text{out}}},\forall k\in \text{NSK}. \\ \end{align}$ (20)
氢阱入口氢气纯度要比氢阱规定的下界高,即
$y_{k,c}^{\text{in}}\ge y_{k,c}^{\text{in,LB}},\forall c=\left\{ {{\text{H}}_{2}} \right\}.$ (21)
6) 与燃料系统相关的公式
燃料系统混合节点的流量平衡和组分平衡为:
$\text{F}{{\text{F}}^{\text{in}}}=\sum\limits_{s\in \text{NSR}}{\text{FS}{{\text{F}}_{s}}+}\sum\limits_{p\in \text{NP}}{\text{FPF}_{p}^{\text{resd}}.}$ (22)
$\begin{align} & \text{F}{{\text{F}}^{\text{in}}}\cdot y_{f,c}^{\text{in}}=\sum\limits_{s\in \text{NSR}}{\text{FS}{{\text{F}}_{s}}\cdot y_{s,c}^{\text{out}}}+ \\ & \sum\limits_{p\in \text{NP}}{\text{FPF}_{p}^{\text{resd}}\cdot y_{p,c}^{\text{resd}}},\forall c\in \text{NC}\text{.} \\ \end{align}$ (23)
其中:FFin是送往燃料系统的总流率; yinf,c是送往燃料系统流股的c组分的平均纯度。
7) 连接和压力约束
连接存在的充分必要条件是流率不为零。式(23)将作为连续变量的流率与二元变量关联起来,则
$\begin{align} & {{F}_{a,b}}-{{z}_{a,b}}\cdot F_{a,b}^{\text{UB}}\le 0,{{F}_{a,b}}-{{z}_{a,b}}\cdot F_{a,b}^{\text{LB}}\ge 0. \\ & z\text{U}{{\text{K}}_{u,k}},z\text{U}{{\text{P}}_{u,p}},z\text{U}{{\text{I}}_{u,i}},z\text{S}{{\text{K}}_{s,k}},z\text{S}{{\text{P}}_{s,p}}, \\ & {{z}_{a,b}}\in z\text{S}{{\text{I}}_{s,i}},z\text{S}{{\text{F}}_{s}},z\text{PK}_{p,k}^{\text{prod}},z\text{PI}_{p,i}^{\text{prod}}, \\ & z\text{PF}_{p}^{\text{resd}},z\text{I}{{\text{K}}_{i,k}},z\text{I}{{\text{P}}_{i,p}},z\text{I}{{\text{I}}_{i,{i}'}} \\ & \text{FU}{{\text{K}}_{u,k}},\text{FU}{{\text{P}}_{u,p}},\text{FU}{{\text{I}}_{u,i}},\text{FS}{{\text{K}}_{s,k}},\text{FS}{{\text{P}}_{s,p}}, \\ & {{F}_{a,b}}\in \text{FS}{{\text{I}}_{s,i}},\text{FS}{{\text{F}}_{s}},\text{FPK}_{p,k}^{\text{prod}},\text{FPI}_{p,i}^{\text{prod}}, \\ & \text{FPF}_{p}^{\text{resd}},\text{FI}{{\text{K}}_{i,k}},\text{FI}{{\text{P}}_{i,p}},\text{FI}{{\text{I}}_{i,{i}'}} \\ \end{align}$ (24)
其中: Fa,b是流率变量; za,b是二元变量,取值为0时表示连接不存在,取值为1时表示连接存在; Fa,bUB是流率变量Fa,b的上界; Fa,bLB是流率变量Fa,b的下界。
压力约束描述如式(25)所示,即
$\begin{align} & \left( {{P}_{b}}-{{P}_{a}} \right)\cdot {{F}_{a,b}}\le 0, \\ & {{P}_{a}}\in \text{PHU}_{u}^{\text{out}},\text{PSR}_{s}^{\text{out}},\text{P}_{p}^{\text{prod}},\text{P}_{p}^{\text{pred}},\text{PI}_{i}^{\text{out}}, \\ & {{P}_{b}}\in \text{PK}_{k}^{\text{in}},\text{PP}_{p}^{\text{in}},\text{PI}_{i}^{\text{in}}. \\ \end{align}$ (25)
8) 目标函数
本文将目标函数设为氢气公用工程用量,目的是使氢气公用工程用量最小,即
$\min \text{FHU=}\sum\limits_{u\in \text{NHU}}{\text{FH}{{\text{U}}_{u}}}.$ (26)
约束条件为式(1)—式(25),其中式(1)为变压吸附模型方程,式(24)为连接约束条件,式(25)为压力约束条件。因式(7)、 式(13)、 式(15)、 式(20)和式(23)存在非线性项,式(24)包含二元变量,该模型为混合整数非线性规划问题(MINLP),利用商业优化软件GAMS(general algebraic modeling system)平台建模,选用DICOPT(discrete and continuous optimizer)作为求解器。
4 案例分析本文采用的案例数据来源于Elkamel等[17]发表的文献。A厂有5个加氢过程,它们分别是加氢裂化单元(HC)、 汽油加氢处理单元(GOHT)、 渣油加氢处理单元(RHT)、 柴油加氢处理单元(DHT)和石脑油加氢处理单元(NHT)。这些耗氢单元的入口相当于过程氢阱,而它们的出口相当于过程氢源。此外,催化重整单元(CR)可以副产,氢气纯度为80%。制氢工厂(HP)的氢气纯度为95%,出口压力为2.069 MPa,最大负荷为89 280 Nm3/h。具体氢源数据如表 1所示,氢阱数据如表 2所示,压缩机数据如表 3所示。本文中PSA的数据是假定的。PSA的入口压力和出口产品部分压力均设定为2.069 MPa,剩余部分的压力被设定为0.103 MPa,出口产品物流的氢气纯度为95%。氢气回收率按式(1)计算。
表 1 氢源数据
氢源流率
Nm3·h-1
氢气纯度
体积百分比
压力
MPa
HP89 280952.069
CR16 182802.069
HC尾气9 099808.276
GOHT尾气11 597752.414
RHT尾气6 466752.759
DHT尾气1 600652.414
NHT尾气2 495601.379


表选项






表 2 氢阱数据
氢阱流率
Nm3·h-1
氢气纯度
体积百分比
压力
MPa
HC41 7189513.79
GOHT38 96593.23.448
RHT19 757904.138
DHT6 068803.448
NHT4 380752.069


表选项






表 3 压缩机数据
压缩机入口压力
MPa
出口压力
MPa
最大负荷
Nm3·h-1
K12.06913.7935 154
K22.06913.7935 154
K32.0694.13835 154


表选项






该提纯回用氢网络优化模型的规模为: 变量数目为218个,等式约束方程数为37个,不等式约束方程数为267个。决策变量有氢源、氢阱、压缩机、 PSA和燃料系统之间的流率、 PSA的入口氢气纯度和氢气回收率; 非决策变量有燃料系统的氢气纯度、吸附和解吸压力比以及吸附选择性。其中,PSA的入口氢气纯度和氢气回收率这2个决策变量对于氢公用工程用量的影响较大,在后文的分析中将得到验证。
分为3种情景进行分析。
情景一?固定θ=0.05,PH/PL=20。
当不设定PSA入口流股氢气纯度,可以得到最小氢公用工程用量为72 677 Nm3/h,此时提纯器入口氢气纯度为76.02%,提纯器的入口原料流率为23 235 Nm3/h,产品气体物流流率为16 502 Nm3/h,氢气回收率为88.75%。优化氢网络如图 2所示。在图 2中,上方的数字表示气体流股的流率,下方的数字表示气体流股的氢气纯度。
图 2 PSA入口氢气纯度为76.02%的优化氢网络
图选项





研究PSA入口流股氢气纯度对回收率和氢公用工程用量的影响。将提纯器入口流股氢气纯度由70%增加至80%,分别得到不同情况下的氢气回收率和最小氢公用工程用量。图 3显示PSA的入口氢气纯度对氢公用工程用量的影响结果。
图 3 PSA的入口氢气纯度对氢公用工程用量的影响
图选项





显然,随着入口流股氢气纯度的增加,氢气回收率逐渐增加,从提高氢气回收率的角度来说,提高PSA的入口流股氢气纯度是有利的。然而,随着入口流股氢气纯度的增加,氢气公用工程用量先降低至最低值,然后缓慢增加。当入口氢气纯度为76.02%时,此时氢公用工程用量最小,为 72 677 Nm3/h。当继续增加入口氢气纯度,PSA的入口气体流率增加,出口产品气体流率增加,氢公用工程用量并未减少,反而略有增加。当PSA入口氢气纯度为80%时,优化后得到的网络结构如图 4所示。在图 4中,上方的数字表示气体流股的流率,下方的数字表示气体流股的氢气纯度。
图 4 PSA入口氢气纯度为80%的优化氢网络
图选项





对比图 2图 4,虽然PSA的产品气流率由 16 502 Nm3/h 增加到22 047 Nm3/h,意味着有更多的产品气被回用至氢网络,直观判断应该认为氢公用工程用量减少,而结果显示略有增加。原因在于为了达到PSA入口氢气纯度为80%,大量的高纯度氢气流股被送往PSA,这样原来被送往氢阱的高纯度氢气流股的流率就会减少,需要增加氢公用工程用量来弥补,从而导致氢公用工程用量略有增加。可见,从氢网络优化的角度来说,一味地增加PSA入口流股氢气纯度以提高回收率的手段并不可取。
情景二?固定θ=0.05,yp,cin=76.02%。
为了分析PSA吸附和解吸压力比对氢气回收率和氢公用工程用量的影响,求解优化模型,结果如图 5所示。由图 5可以看出,当吸附和解吸压力比由5增加到30,氢气回收率由70.0%增加到90.8%,因而会有更多的产品气回用至氢网络,氢公用工程用量由76 162 Nm3/h减少到72 290 Nm3/h,减少了5.1%。然而,吸附和解吸压力比增大,意味着PSA的操作费用会增加。
图 5 PSA的吸附和解吸压力比对氢公用工程用量的影响
图选项





情景三?固定PH/PL=20,yp,cin=76.02%。
为了分析PSA吸附选择性对氢气回收率和氢公用工程用量的影响,求解模型,结果如图 6所示。由图 6可以看出,当吸附选择性由0.8减少到 0.02时,氢气回收率由18.7%增加到91.6%,因而会有更多的产品气回用至氢网络,氢公用工程量由85 705 Nm3/h减少到72 156 Nm3/h,减少了15.8%。一般保持较小吸附选择性,以提高氢气回收率。
图 6 PSA的吸附选择性对氢公用工程用量的影响
图选项





5 结论本文提出了耦合变压吸附简化模型的提纯回用氢网络优化的数学模型,并实现两者的协调优化。案例分析表明,随着PSA入口流股氢气纯度由70%增加到80%,氢公用工程用量先减少至最小值,然后缓慢增加。本案例中PSA最优入口氢气纯度为76.02%,此时氢公用工程用量最小,为 72 677 Nm3/h,比原始氢公用工程用量(89 280 Nm3/h)减少了18.6%。当PSA吸附和解吸压力比PH/PL由5增加到30,氢公用工程用量减少了 5.1%。当PSA吸附选择性由0.8减少到0.02,氢公用工程用量减少了15.8%。

参考文献
[1] Journal of Central South University(Science and Technology), 41(2):649-654.-->Alves J J, Towler G P. Analysis of refinery hydrogen distribution systems[J]. Industrial & Engineering Chemistry Research, 2002, 41(23) : 5759–5769.
[2] Journal of Central South University(Science and Technology), 41(2):649-654.-->El-Halwagi M M, Gabriel F, Harell D. Rigorous graphical targeting for resource conservation via material recycle/reuse networks[J]. Industrial & Engineering Chemistry Research, 2003, 42(19) : 4319–4328.
[3] Journal of Central South University(Science and Technology), 41(2):649-654.-->Foo D C Y, Manan Z A. Setting the minimum utility gas flowrate targets using cascade analysis technique[J]. Industrial & Engineering Chemistry Research, 2006, 45(17) : 5986–5995.
[4] Journal of Central South University(Science and Technology), 41(2):649-654.-->Agrawal V, Shenoy U V. Unified conceptual approach to targeting and design of water and hydrogen networks[J]. American Institute of Chemical Engineers Journal, 2006, 52(3) : 1071–1082.DOI:10.1002/(ISSN)1547-5905
[5] Journal of Central South University(Science and Technology), 41(2):649-654.-->DENG Chun, ZHOU Yuhang, CHEN Chengliang, et al. Systematic approach for targeting interplant hydrogen networks[J]. Energy, 2015, 90(1) : 68–88.
[6] Journal of Central South University(Science and Technology), 41(2):649-654.-->Hallale N, LIU Fang. Refinery hydrogen management for clean fuels production[J]. Advances in Environmental Research, 2001, 6(1) : 81–98.DOI:10.1016/S1093-0191(01)00112-5
[7] Journal of Central South University(Science and Technology), 41(2):649-654.-->Ahmad M I, ZHANG Nan, Jobson M. Modelling and optimisation for design of hydrogen networks for multi-period operation[J]. Journal of Cleaner Production, 2010, 18(9) : 889–899.DOI:10.1016/j.jclepro.2010.01.003
[8] Journal of Central South University(Science and Technology), 41(2):649-654.-->LIAO Zuwei, WANG Jingdai, YANG Yongrong, et al. Integrating purifiers in refinery hydrogen networks:A retrofit case study[J]. Journal of Cleaner Production, 2010, 18(3) : 233–241.DOI:10.1016/j.jclepro.2009.10.011
[9] Journal of Central South University(Science and Technology), 41(2):649-654.-->JIA Nan, ZHANG Nan. Multi-component optimisation for refinery hydrogen networks[J]. Energy, 2011, 36(8) : 4663–4670.DOI:10.1016/j.energy.2011.03.040
[10] Journal of Central South University(Science and Technology), 41(2):649-654.-->DENG Chun, PAN Huaimin, LI Yantao, et al. Comparative analysis of different scenarios for the synthesis of refinery hydrogen network[J]. Applied Thermal Engineering, 2014, 70(2) : 1162–1179.DOI:10.1016/j.applthermaleng.2014.04.036
[11] Journal of Central South University(Science and Technology), 41(2):649-654.-->DENG Chun, PAN Huaimin, LEE Juiyuan, et al. Synthesis of hydrogen network with hydrogen header of intermediate purity[J]. International Journal of Hydrogen Energy, 2014, 39(25) : 13049–13062.DOI:10.1016/j.ijhydene.2014.06.129
[12] Journal of Central South University(Science and Technology), 41(2):649-654.-->Umana B, Shoaib A, ZHANG Nan, et al. Integrating hydroprocessors in refinery hydrogen network optimisation[J]. Applied Energy, 2014, 133(10) : 169–182.
[13] Journal of Central South University(Science and Technology), 41(2):649-654.-->LIU Fang, ZHANG Nan. Strategy of purifier selection and integration in hydrogen networks[J]. Chemical Engineering Research & Design, 2004, 82(10) : 1315–1330.
[14] Journal of Central South University(Science and Technology), 41(2):649-654.-->Lacava A I, Shirley A I, Ramachandran R. How to specify pressure-swing adsorption units[J]. Chemical Engineering, 1998, 105(6) : 110–118.
[15] Journal of Central South University(Science and Technology), 41(2):649-654.-->Knaebel, K S. The basics of adsorber design[J]. Chemical Engineering, 1999, 106(4) : 92–101.
[16] Journal of Central South University(Science and Technology), 41(2):649-654.-->邓春, 周业扬, 陈杰, 等. 石化园区厂际提纯回用氢气系统优化[J]. 化工学报, 2014, 65(12) : 4914–4920.DENG Chun, ZHOU Yeyang, CHEN Jie, et al. Optimization of inter-plant hydrogen system with purification reuse in petrochemical complex[J]. Chemical Industry and Engineering Society of China Journal, 2014, 65(12) : 4914–4920.(in Chinese)
[17] Journal of Central South University(Science and Technology), 41(2):649-654.-->Elkamel A, Alhajri I, Almansoori A, et al. Integration of hydrogen management in refinery planning with rigorous process models and product quality specifications[J]. International Journal of Process Systems Engineering, 2011, 1(3) : 302–330.

相关话题/工程 网络

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 面向复杂网络的威胁度量及聚合方法
    邓辉,刘晖,张宝峰,毛军捷,郭颖,熊琦,谢仕华中国信息安全测评中心,北京100085收稿日期:2016-01-25基金项目:国家自然科学基金资助项目(61472448)作者简介:邓辉(1985-),女,助理研究员。E-mail:gcdh2014@126.com摘要:在复杂网络中,威胁模型结构庞大、行 ...
    本站小编 Free考研考试 2020-04-15
  • 高校师生群体间人员接触网络研究
    马勋,申世飞,倪顺江,雍诺清华大学工程物理系,公共安全研究院,北京100084收稿日期:2015-12-27基金项目:国家自然科学基金资助项目(71203118)作者简介:马勋(1991-),男,博士研究生。通讯作者:倪顺江,讲师,E-mail:sjni@tsinghua.edu.cn摘要:为了研究 ...
    本站小编 Free考研考试 2020-04-15
  • 跪求老师解答,18年毕业的 专科学习的是电子信息工程技术,可以报考贵校的小学教育专硕专业吗?
    提问问题:跪求老师解答,18年毕业的专科学习的是电子信息工程技术,可以报考贵校的小学教育专硕专业吗?学院:教育学部提问人:15***09时间:2019-09-2309:06提问内容:跪求老师解答,18年毕业的专科学习的是电子信息工程技术,可以报考贵校的小学教育专硕专业吗?请问一下,同等学力有什么具体 ...
    本站小编 天津师范大学 2019-11-27
  • 软件工程专业问题
    提问问题:软件工程专业问题学院:计算机与信息工程学院提问人:15***84时间:2019-09-1910:38提问内容:老师您好!新的专业目录发布,请问老师能否提供下软工专业命题数据结构与C++程序设计的参考书,以及题型和分数分布;软件工程专硕培育方向是什么?回复内容:推荐书目请浏览计算机与信息工程 ...
    本站小编 天津师范大学 2019-11-27
  • 软件工程专业相关问题
    提问问题:软件工程专业相关问题学院:计算机与信息工程学院提问人:17***35时间:2018-09-2115:38提问内容:老师您好!贵校2019年招生简章中将计算机与信息学院的软件工程专业的专业课自命题由软件基础综合改为数据结构与C++程序设计,1.请问老师能否提供修改后的考试大纲,参考书目,题型 ...
    本站小编 天津师范大学 2019-11-27
  • 网络教育本科
    提问问题:网络教育本科学院:教育科学学院提问人:15***99时间:2018-09-2011:32提问内容:老师,你好请问网络教育本科明年毕业但现在学历校验未通过报名时能通过吗回复内容:2019年9月1日前可取得国家承认本科毕业证书的网络教育本科生,须凭学生证方可办理网上报名现场确认手续。 ...
    本站小编 天津师范大学 2019-11-27
  • 网络教育
    提问问题:网络教育学院:提问人:15***63时间:2016-09-1912:03提问内容:老师您好,网络教育本科的毕业生,是按普通本科毕业生的身份对待,还是按同等学力的身份对待?回复内容:按普通本科报考即可 ...
    本站小编 天津师范大学 2019-11-27
  • 光学工程
    提问问题:光学工程学院:激光与红外系统集成技术教育部重点实验室提问人:15***58时间:2019-09-2015:48提问内容:您好!想了解下激光与红外系统集成技术实验室的光工与信息工程学院难度有没有区别?院线如何?查到的报录比只有信息工程学院的回复内容:建议查看“山东大学研究生招生信息网”中的报 ...
    本站小编 山东大学 2019-11-26
  • 动力工程
    提问问题:动力工程学院:能源与动力工程学院提问人:17***21时间:2019-09-2011:53提问内容:老师您好,请问动力工程专硕今年缩招了一半?后期调整会变动很大吗?回复内容:以官网专业目录中的招生人数为准,学校后期可能会根据推免人数有所调整(预计是10月下旬)。 ...
    本站小编 山东大学 2019-11-26
  • 125603工业工程管理
    提问问题:125603工业工程管理学院:管理学院提问人:we***om时间:2019-09-2011:33提问内容:请问全日制专科2012年毕业,可以报考125603工业工程管理专业吗?是否全日制,谢谢回复内容:可以报考,是非全日制。 ...
    本站小编 山东大学 2019-11-26