变压吸附氢气纯化过程瞬态分析

doi:10.11949/0438-1157.20211103

摘要/Abstract

摘要：

变压吸附技术是工业上生产高纯氢气最常用的方法之一。然而，在实际生产过程中无法观察到塔内各组分在不同时刻的分布状态，因此借助模拟的手段来研究从投料至系统达到循环稳态期间各组分在塔内的动态变化规律，进而指导工艺改进是很有必要的。采用活性炭和5A分子筛为吸附剂，设计了八塔变压吸附工艺从蒸汽甲烷重整气中纯化氢气，模拟了变压吸附制氢开车过程，分析了开车过程中塔内各组分在吸附、顺放以及冲洗三个阶段以及循环稳态后吸附阶段瞬态吸附行为和塔内温度变化。结果表明，在吸附以及顺放等过程中重组分会随着循环周期向塔顶移动。这一现象是组分间竞争吸附和冲洗再生方式下重组分在床层底部累积两个作用因素共同导致的。这些因素在一定程度上也会造成CO的吸附前沿在吸附阶段就过多进入5A分子筛上，使得CO含量成为限制工艺性能的主要因素。

关键词: 变压吸附, 活性炭, 5A分子筛, 制氢, 数值模拟, 瞬态分析

Abstract:

Pressure swing adsorption (PSA) is one of the most common technologies to produce high-purity hydrogen in the industry. However, in the actual production process, the distribution state of each component in the bed at different time cannot be observed. For this reason, simulation methods are used to study the dynamic change of each component in the bed from the feeding to the system reaching the steady state, which is essential to guide process improvement. In this study, an eight-bed PSA process was designed to purify hydrogen from steam methane reforming (SMR) gas and the activated carbon and 5A zeolite were used as adsorbent. The start-up procedure of PSA hydrogen production was performed, and the transient adsorption behaviors of each component in the bed during the three stages of adsorption, sequential release and flushing, as well as the transient adsorption behavior during the adsorption stage and the temperature variation in the bed after the system reached the cyclic steady state were analyzed. The results show that the heavy components move towards the top of the bed with the cycle during the adsorption as well as the sequential release. This phenomenon is the result of both inter-component competition for adsorption and the accumulation of fractions at the bottom of the bed in the flush regeneration mode. These factors can also, to some extent, cause the adsorption front of CO to enter too much on the 5A zeolite during the adsorption phase, making the CO content become a major limiting factor for process performance.

Key words: pressure swing adsorption, activated carbon, 5A zeolite, hydrogen production, numerical simulation, transient analysis

中图分类号:

TQ 028.1

张超, 陈健, 殷文华, 沈圆辉, 钮朝阳, 余秀鑫, 张东辉, 唐忠利. 变压吸附氢气纯化过程瞬态分析[J]. 化工学报, 2022, 73(1): 308-321.

Chao ZHANG, Jian CHEN, Wenhua YIN, Yuanhui SHEN, Zhaoyang NIU, Xiuxin YU, Donghui ZHANG, Zhongli TANG. Transient analysis of pressure swing adsorption hydrogen purification process[J]. CIESC Journal, 2022, 73(1): 308-321.

图/表 18

图1 八塔PSA整体工艺结构图

Fig.1 Configuration diagram of eight-bed PSA process

图2 PSA循环工艺示意图

Fig.2 Process schematic diagram of PSA

表1 八塔PSA时序

Table 1 The schedule of eight-bed PSA

时间/s	塔1	塔2	塔3	塔4	塔5	塔6	塔7	塔8
40	AD1	ER1	ER3	RP1	BD	PP1	ED3	ED1
40	AD2	PR	ER2	ER4	RP2	PP2	ED4	ED2
40	ED1	AD1	ER1	ER3	RP1	BD	PP1	ED3
40	ED2	AD2	PR	ER2	ER4	RP2	PP2	ED4
40	ED3	ED1	AD1	ER1	ER3	RP1	BD	PP1
40	ED4	ED2	AD2	PR	ER2	ER4	RP2	PP2
40	PP1	ED3	ED1	AD1	ER1	ER3	RP1	BD
40	PP2	ED4	ED2	AD2	PR	ER2	ER4	RP2
40	BD	PP1	ED3	ED1	AD1	ER1	ER3	RP1
40	RP2	PP2	ED4	ED2	AD2	PR	ER2	ER4
40	RP1	BD	PP1	ED3	ED1	AD1	ER1	ER3
40	ER4	RP2	PP2	ED4	ED2	AD2	PR	ER2
40	ER3	RP1	BD	PP1	ED3	ED1	AD1	ER1
40	ER2	ER4	RP2	PP2	ED4	ED2	AD2	PR
40	ER1	ER3	RP1	BD	PP1	ED3	ED1	AD1
40	PR	ER2	ER4	RP2	PP2	ED4	ED2	AD2

表2 吸附塔数学模型

Table 2 Mathematical model of adsorption bed

	数学模型
质量守恒	$- ε b D a x, i ? 2 c i ? z 2 + ε b ? (v g c i) ? z + ε b + (1 - ε b) ε p ? c i ? t + ρ s (1 - ε b) ? q i ? t = 0$	(1)
	$ε b ? (v g c) ? z + ε b + (1 - ε b) ε p ? c ? t + ρ s (1 - ε b) ∑ i = 1 n ? q i ? t = 0$	(2)
	$D a x, i = 0.73 D m, i + v g r p ε b (1 + 9.49 ε b D m, i 2 v g r p)$	(3)
能量守恒	$- ε b k g ? 2 T g ? z 2 + C v g v g ρ g ? T g ? z + ε b C v g ρ g ? T g ? t + P ? v g ? z + h f a p (T g - T s) + 4 h w D b (T g - T 0) = 0$	(4)
	$- k s ? 2 T s ? z 2 + C p s ρ s ? T s ? t + ρ s ∑ i = 1 n (C p a, i q i) ? T s ? t + ρ s ∑ i = 1 n (Δ H i ? q i ? t) - h f a p (T g - T s) = 0$	(5)
	$k w ? 2 T w ? z 2 + C p w ρ w ? T w ? t - h w 4 D b (D b + W t) 2 - D b 2 (T g - T w) + h a m b 4 (D b + W t) 2 (D b + W t) 2 - D b 2 (T w - T a m b) = 0$	(6)
Ergun方程	$- ? P ? z = 150 μ (1 - ε b) 2 ε b 3 (2 r p ψ) 2 v g + 1.75 (1 - ε b) ρ g 2 r p ψ ε b 3 \| v g \| v g$	(7)
LDF方程	$? q i ? t = k L D F, i (q i * - q i)$	(8)
吸附等温式	$q i * = q m, i b i P i 1 + ∑ i = 1 n b i P i$ , $b i = b 0 e x p - Δ H i R g T$	(9)

表2 吸附塔数学模型

Table 2 Mathematical model of adsorption bed

	数学模型
质量守恒	$- ε b D a x, i ? 2 c i ? z 2 + ε b ? (v g c i) ? z + ε b + (1 - ε b) ε p ? c i ? t + ρ s (1 - ε b) ? q i ? t = 0$	(1)
	$ε b ? (v g c) ? z + ε b + (1 - ε b) ε p ? c ? t + ρ s (1 - ε b) ∑ i = 1 n ? q i ? t = 0$	(2)
	$D a x, i = 0.73 D m, i + v g r p ε b (1 + 9.49 ε b D m, i 2 v g r p)$	(3)
能量守恒	$- ε b k g ? 2 T g ? z 2 + C v g v g ρ g ? T g ? z + ε b C v g ρ g ? T g ? t + P ? v g ? z + h f a p (T g - T s) + 4 h w D b (T g - T 0) = 0$	(4)
	$- k s ? 2 T s ? z 2 + C p s ρ s ? T s ? t + ρ s ∑ i = 1 n (C p a, i q i) ? T s ? t + ρ s ∑ i = 1 n (Δ H i ? q i ? t) - h f a p (T g - T s) = 0$	(5)
	$k w ? 2 T w ? z 2 + C p w ρ w ? T w ? t - h w 4 D b (D b + W t) 2 - D b 2 (T g - T w) + h a m b 4 (D b + W t) 2 (D b + W t) 2 - D b 2 (T w - T a m b) = 0$	(6)
Ergun方程	$- ? P ? z = 150 μ (1 - ε b) 2 ε b 3 (2 r p ψ) 2 v g + 1.75 (1 - ε b) ρ g 2 r p ψ ε b 3 \| v g \| v g$	(7)
LDF方程	$? q i ? t = k L D F, i (q i * - q i)$	(8)
吸附等温式	$q i * = q m, i b i P i 1 + ∑ i = 1 n b i P i$ , $b i = b 0 e x p - Δ H i R g T$	(9)

表3 工艺性能评价指标

Table 3 Process performance indicator

指标	数学模型
纯度	$p u r i t y H 2 = ∫ 0 t a d F o u t y o u t, H 2 d t ∫ 0 t a d F o u t d t$	(10)
回收率	recovery $H 2 = ∫ 0 t a d F o u t y o u t, H 2 d t ∫ 0 t a d F i n y i n, H 2 d t$	(11)
生产量	$p r o d u c t i v i t y H 2 = 60 R T s t a n ∫ 0 t a d F o u t y o u t, H 2 d t P s t a n t c y c l e$	(12)

表3 工艺性能评价指标

Table 3 Process performance indicator

指标	数学模型
纯度	$p u r i t y H 2 = ∫ 0 t a d F o u t y o u t, H 2 d t ∫ 0 t a d F o u t d t$	(10)
回收率	recovery $H 2 = ∫ 0 t a d F o u t y o u t, H 2 d t ∫ 0 t a d F i n y i n, H 2 d t$	(11)
生产量	$p r o d u c t i v i t y H 2 = 60 R T s t a n ∫ 0 t a d F o u t y o u t, H 2 d t P s t a n t c y c l e$	(12)

表4 活性炭和5A分子筛吸附剂参数

Table 4 Parameters of activated carbon and 5A zeolite adsorbent

物性参数	活性炭	5A分子筛
ρ_b /(kg·m^-3)	522.0	698.0
C_p_s/(kJ·kg^-1·K^-1)	1.047	0.92
r_p/m	0.00115	0.00157
ε_b	0.433	0.357
ε_p	0.61	0.65
$ψ$	0.9	1.0
k_LDF,CH4/s^-1	0.195	0.147
k_LDF,CO/s^-1	0.15	0.063
k_LDF,CO2/s^-1	0.0355	0.0135
k_LDF,H2/s^-1	0.7	0.7
k_LDF,N2/s^-1	3.0	2.5

表4 活性炭和5A分子筛吸附剂参数

Table 4 Parameters of activated carbon and 5A zeolite adsorbent

物性参数	活性炭	5A分子筛
ρ_b /(kg·m^-3)	522.0	698.0
C_p_s/(kJ·kg^-1·K^-1)	1.047	0.92
r_p/m	0.00115	0.00157
ε_b	0.433	0.357
ε_p	0.61	0.65
$ψ$	0.9	1.0
k_LDF,CH4/s^-1	0.195	0.147
k_LDF,CO/s^-1	0.15	0.063
k_LDF,CO2/s^-1	0.0355	0.0135
k_LDF,H2/s^-1	0.7	0.7
k_LDF,N2/s^-1	3.0	2.5

表5 吸附床参数

Table 5 Parameters of adsorption bed

参数	数值
H_b/m	1.0
W_t/m	0.002
D_b/m	0.2
$ρ w$ /(kg·m^-3)	7800
C_pw/(kJ·kg^-1·K^-1)	0.5024
h_w/(W·m^-2·K^-1)	94.0
T_amb/K	298.15

表5 吸附床参数

Table 5 Parameters of adsorption bed

参数	数值
H_b/m	1.0
W_t/m	0.002
D_b/m	0.2
$ρ w$ /(kg·m^-3)	7800
C_pw/(kJ·kg^-1·K^-1)	0.5024
h_w/(W·m^-2·K^-1)	94.0
T_amb/K	298.15

表6 扩展型Langmuir吸附模型拟合参数

Table 6 Extended Langmuir adsorption model fitting parameters

组分	T/K	活性炭				5A 分子筛
组分	T/K	q_m/(mmol·g^-1)	b₀/bar^-1	R	ΔH/(kJ·mol^-1)	q_m/(mmol·g^-1)	b₀/bar^-1	R	ΔH/(kJ·mol^-1)
H₂	298.15	3.85	5.48×10^-5	0.9996	-12.84	0.76	7.70×10^-4	0.9995	-9.23
	308.15	4.32	5.12×10^-5	0.9997		0.67	8.57×10^-4	0.9996
	318.15	2.21	1.06×10^-4	0.9994		0.65	8.47×10^-4	0.9997
CO	298.15	6.72	9.10×10^-6	0.9998	-22.58	2.24	6.14×10^-6	0.9991	-29.77
	308.15	6.22	9.13×10^-6	0.9998		2.20	6.58×10^-6	0.9993
	318.15	5.50	9.68×10^-6	0.9997		2.27	6.14×10^-6	0.9997
CO₂	298.15	7.96	2.52×10^-6	1.0000	-29.08	4.35	7.60×10^-6	0.9990	-35.97
	308.15	7.32	3.20×10^-6	0.9998		4.14	9.74×10^-6	0.9982
	318.15	6.59	4.17×10^-6	1.0000		3.96	9.12×10^-6	0.9956
CH₄	298.15	7.21	1.73×10^-5	0.9999	-22.70	3.73	5.91×10^-5	0.9996	-20.64
	308.15	6.75	1.98×10^-5	0.9998		3.54	5.37×10^-5	1.0000
	318.15	6.14	2.38×10^-5	0.9998		3.63	5.12×10^-5	1.0000
N₂	298.15	3.48	1.87×10^-4	1.0000	-15.99	2.75	3.98×10^-5	0.9999	-20.65
	308.15	3.38	1.94×10^-4	0.9997		2.67	4.16×10^-5	0.9997
	318.15	3.30	2.01×10^-4	0.9998		2.59	4.33×10^-5	1.0000

表7 吸附床初始条件

Table 7 Initial conditions of adsorption bed

Initial conditions

$y i z = 0,$ $y i n e r t z = 1,$ $q i z = 0,$ $q i n e r t z = q * i n e r t$ ,

$c i z = P y i R g T$ , $P z = P f e e d$ , $T g z = T s z = T w z = T f e e d$

表7 吸附床初始条件

Table 7 Initial conditions of adsorption bed

Initial conditions

$y i z = 0,$ $y i n e r t z = 1,$ $q i z = 0,$ $q i n e r t z = q * i n e r t$ ,

$c i z = P y i R g T$ , $P z = P f e e d$ , $T g z = T s z = T w z = T f e e d$

表8 PSA过程中各步的边界条件

Table 8 The boundary conditions for each step of PSA process

Step	z=0	z=L
吸附 (AD1, AD2)
	$P z = 0 = P i n l e t$	$P z = L = P o u t l e t$
	$- ε b D a x, i ? c i ? z z = 0 = u z = 0 (c i z = 0 + - c i z = 0 -)$	$? c i ? z z = L = 0$
	$k g C v g ρ g ? T g ? z z = 0 = u z = 0 (T g z = 0 + - T g z = 0 -)$	$? T g ? z z = L = 0$
	$k s C p s ρ s ? T s ? z z = 0 = u z = 0 (T s z = 0 + - T s z = 0 -)$	$? T s ? z z = L = 0$
	$k w C p w ρ w ? T w ? z z = 0 = u z = 0 (T w z = 0 + - T w z = 0 -)$	$? T w ? z z = L = 0$
均压降 (ED1, ED2, ED3, ED4)
	$u z = 0 = 0$	$P z = L = P o u t l e t$
	$? c i ? z z = 0 = 0$	$? c i ? z z = L = 0$
	$? T g ? z z = 0 = 0$	$? T g ? z z = L = 0$
	$? T s ? z z = 0 = 0$	$? T s ? z z = L = 0$
	$? T w ? z z = 0 = 0$	$? T w ? z z = L = 0$
顺放 (PP1, PP2)
	$u z = 0 = 0$	$P z = L = P o u t l e t$
	$? c i ? z z = 0 = 0$	$? c i ? z z = L = 0$
	$? T g ? z z = 0 = 0$	$? T g ? z z = L = 0$
	$? T s ? z z = 0 = 0$	$? T s ? z z = L = 0$
	$? T w ? z z = 0 = 0$	$? T w ? z z = L = 0$

表8 PSA过程中各步的边界条件

Table 8 The boundary conditions for each step of PSA process

Step	z=0	z=L
吸附 (AD1, AD2)
	$P z = 0 = P i n l e t$	$P z = L = P o u t l e t$
	$- ε b D a x, i ? c i ? z z = 0 = u z = 0 (c i z = 0 + - c i z = 0 -)$	$? c i ? z z = L = 0$
	$k g C v g ρ g ? T g ? z z = 0 = u z = 0 (T g z = 0 + - T g z = 0 -)$	$? T g ? z z = L = 0$
	$k s C p s ρ s ? T s ? z z = 0 = u z = 0 (T s z = 0 + - T s z = 0 -)$	$? T s ? z z = L = 0$
	$k w C p w ρ w ? T w ? z z = 0 = u z = 0 (T w z = 0 + - T w z = 0 -)$	$? T w ? z z = L = 0$
均压降 (ED1, ED2, ED3, ED4)
	$u z = 0 = 0$	$P z = L = P o u t l e t$
	$? c i ? z z = 0 = 0$	$? c i ? z z = L = 0$
	$? T g ? z z = 0 = 0$	$? T g ? z z = L = 0$
	$? T s ? z z = 0 = 0$	$? T s ? z z = L = 0$
	$? T w ? z z = 0 = 0$	$? T w ? z z = L = 0$
顺放 (PP1, PP2)
	$u z = 0 = 0$	$P z = L = P o u t l e t$
	$? c i ? z z = 0 = 0$	$? c i ? z z = L = 0$
	$? T g ? z z = 0 = 0$	$? T g ? z z = L = 0$
	$? T s ? z z = 0 = 0$	$? T s ? z z = L = 0$
	$? T w ? z z = 0 = 0$	$? T w ? z z = L = 0$

图3 进料流量、活性炭/5A分子筛填充比和冲洗进料比对PSA性能的影响

Fig.3 Effect of feed flow rate, layer ratio of AC/5A zeolite and purge-to-feed ratio on PSA performance

图4 塔顶出口各组分摩尔分数随时间的变化曲线

Fig.4 Dynamic profile of mole fraction of each component in the top of bed

图5 吸附结束时N2在轴向上的气相浓度分布

Fig.5 Gas phase concentration profile of N2 at the end of adsorption

图6 吸附结束时H2、CO、CO2和CH4在轴向上的气相浓度分布

Fig.6 Gas phase concentration of H2, CO, CO2 and CH4 at the end of adsorption

图7 顺放结束时H2、CO、CO2、CH4在轴向上的气相浓度分布

Fig.7 Gas phase concentration profiles of H2, CO, CO2 and CH4 at the end of sequential release

图8 冲洗结束时H2、CO、CO2、CH4在轴向上的固相浓度分布

Fig.8 Solid phase concentration profiles of H2, CO, CO2 and CH4 at the end of flushing

图9 循环稳态后H2、CO、CO2、CH4在吸附阶段随着轴向距离的气相浓度分布曲线

Fig.9 Gas phase concentration of H2, CO, CO2 and CH4 in the adsorption stage after cycling steady state

图10 循环稳态时吸附塔内气相温度轴向分布和压力分布以及塔底解吸气分布

Fig.10 Temperature distribution, pressure profile within adsorption bed, and desorption gas distribution at the bottom of the bed at cycle steady state

参考文献 30

1	Momirlan M, Veziroglu T N. The properties of hydrogen as fuel tomorrow in sustainable energy system for a cleaner planet[J]. International Journal of Hydrogen Energy, 2005, 30(7): 795-802.
2	Yang S I, Choi D Y, Jang S C, et al. Hydrogen separation by multi-bed pressure swing adsorption of synthesis gas[J]. Adsorption, 2008, 14(4/5): 583-590.
3	Delgado J A, Águeda V I, Uguina M A, et al. Adsorption and diffusion of H₂, CO, CH₄, and CO₂ in BPL activated carbon and 13X zeolite: evaluation of performance in pressure swing adsorption hydrogen purification by simulation[J]. Industrial & Engineering Chemistry Research, 2014, 53(40): 15414-15426.
4	Zhu X C, Shi Y X, Li S, et al. Two-train elevated-temperature pressure swing adsorption for high-purity hydrogen production[J]. Applied Energy, 2018, 229: 1061-1071.
5	Ritter J A, Ebner A D. State-of-the-art adsorption and membrane separation processes for hydrogen production in the chemical and petrochemical industries[J]. Separation Science and Technology, 2007, 42(6): 1123-1193.
6	Sircar S, Kratz W C. Simultaneous production of hydrogen and carbon dioxide from steam reformer off-gas by pressure swing adsorption[J]. Separation Science and Technology, 1988, 23(14/15): 2397-2415.
7	Sircar S, Golden T C. Purification of hydrogen by pressure swing adsorption[J]. Separation Science and Technology, 2000, 35(5): 667-687.
8	Brea P, Delgado J A, Águeda V I, et al. Modeling of breakthrough curves of N₂, CH₄, CO, CO₂ and a SMR type off-gas mixture on a fixed bed of BPL activated carbon[J]. Separation and Purification Technology, 2017, 179: 61-71.
9	Ribeiro A M, Grande C A, Lopes F V S, et al. Four beds pressure swing adsorption for hydrogen purification: case of humid feed and activated carbon beds[J]. AIChE Journal, 2009, 55(9): 2292-2302.
10	Ahn S, You Y W, Lee D G, et al. Layered two- and four-bed PSA processes for H₂ recovery from coal gas[J]. Chemical Engineering Science, 2012, 68(1): 413-423.
11	Luberti M, Friedrich D, Brandani S, et al. Design of a H₂ PSA for cogeneration of ultrapure hydrogen and power at an advanced integrated gasification combined cycle with pre-combustion capture[J]. Adsorption, 2014, 20(2/3): 511-524.
12	Moon D K, Lee D G, Lee C H. H₂ pressure swing adsorption for high pressure syngas from an integrated gasification combined cycle with a carbon capture process[J]. Applied Energy, 2016, 183: 760-774.
13	Golmakani A, Fatemi S, Tamnanloo J. Investigating PSA, VSA, and TSA methods in SMR unit of refineries for hydrogen production with fuel cell specification[J]. Separation and Purification Technology, 2017, 176: 73-91.
14	Moon D K, Park Y, Oh H T, et al. Performance analysis of an eight-layered bed PSA process for H₂ recovery from IGCC with pre-combustion carbon capture[J]. Energy Conversion and Management, 2018, 156: 202-214.
15	Shi W R, Yang H W, Shen Y H, et al. Two-stage PSA/VSA to produce H₂ with CO₂ capture via steam methane reforming (SMR)[J]. International Journal of Hydrogen Energy, 2018, 43(41): 19057-19074.
16	Park Y, Kang J H, Moon D K, et al. Parallel and series multi-bed pressure swing adsorption processes for H₂ recovery from a lean hydrogen mixture[J]. Chemical Engineering Journal, 2021, 408: 127299.
17	Divekar S, Arya A, Hanif A, et al. Recovery of hydrogen and carbon dioxide from hydrogen PSA tail gas by vacuum swing adsorption[J]. Separation and Purification Technology, 2021, 254: 117113.
18	Majlan E H, Daud W R W, Iyuke S E, et al. Hydrogen purification using compact pressure swing adsorption system for fuel cell[J]. International Journal of Hydrogen Energy, 2009, 34(6): 2771-2777.
19	Lopes F V S, Grande C A, Rodrigues A E. Activated carbon for hydrogen purification by pressure swing adsorption: multicomponent breakthrough curves and PSA performance[J]. Chemical Engineering Science, 2011, 66(3): 303-317.
20	Abdeljaoued A, Relvas F, Mendes A, et al. Simulation and experimental results of a PSA process for production of hydrogen used in fuel cells[J]. Journal of Environmental Chemical Engineering, 2018, 6(1): 338-355.
21	Zhang N N, Xiao J S, Bénard P, et al. Single-and double-bed pressure swing adsorption processes for H₂/CO syngas separation[J]. International Journal of Hydrogen Energy, 2019, 44(48): 26405-26418.
22	Li H R, Liao Z W, Sun J Y, et al. Modelling and simulation of two-bed PSA process for separating H₂ from methane steam reforming[J]. Chinese Journal of Chemical Engineering, 2019, 27(8): 1870-1878.
23	Yáñez M, Relvas F, Ortiz A, et al. PSA purification of waste hydrogen from ammonia plants to fuel cell grade[J]. Separation and Purification Technology, 2020, 240: 116334.
24	You Y W, Lee D G, Yoon K Y, et al. H₂ PSA purifier for CO removal from hydrogen mixtures[J]. International Journal of Hydrogen Energy, 2012, 37(23): 18175-18186.
25	Golmakani A, Nabavi S A, Manović V. Effect of impurities on ultra-pure hydrogen production by pressure vacuum swing adsorption[J]. Journal of Industrial and Engineering Chemistry, 2020, 82: 278-289.
26	Liu Z, Grande C A, Li P, et al. Multi-bed vacuum pressure swing adsorption for carbon dioxide capture from flue gas[J]. Separation and Purification Technology, 2011, 81(3): 307-317.
27	Yang H W, Yin C B, Jiang B, et al. Optimization and analysis of a VPSA process for N₂/CH₄ separation[J]. Separation and Purification Technology, 2014, 134: 232-240.
28	Sun W N, Shen Y H, Zhang D H, et al. A systematic simulation and proposed optimization of the pressure swing adsorption process for N₂/CH₄ separation under external disturbances[J]. Industrial & Engineering Chemistry Research, 2015, 54(30): 7489-7501.
29	Tang X, Wang Z F, Ripepi N, et al. Adsorption affinity of different types of coal: mean isosteric heat of adsorption[J]. Energy & Fuels, 2015, 29(6): 3609-3615.
30	Wu T B, Shen Y H, Feng L, et al. Adsorption properties of N₂O on zeolite 5A, 13X, activated carbon, ZSM-5, and silica gel[J]. Journal of Chemical & Engineering Data, 2019, 64(8): 3473-3482.

[1]	叶展羽, 山訸, 徐震原. 用于太阳能蒸发的折纸式蒸发器性能仿真[J]. 化工学报, 2023, 74(S1): 132-140.
[2]	张义飞, 刘舫辰, 张双星, 杜文静. 超临界二氧化碳用印刷电路板式换热器性能分析[J]. 化工学报, 2023, 74(S1): 183-190.
[3]	王志国, 薛孟, 董芋双, 张田震, 秦晓凯, 韩强. 基于裂隙粗糙性表征方法的地热岩体热流耦合数值模拟与分析[J]. 化工学报, 2023, 74(S1): 223-234.
[4]	宋嘉豪, 王文. 斯特林发动机与高温热管耦合运行特性研究[J]. 化工学报, 2023, 74(S1): 287-294.
[5]	张思雨, 殷勇高, 贾鹏琦, 叶威. 双U型地埋管群跨季节蓄热特性研究[J]. 化工学报, 2023, 74(S1): 295-301.
[6]	黄琮琪, 吴一梅, 陈建业, 邵双全. 碱性电解水制氢装置热管理系统仿真研究[J]. 化工学报, 2023, 74(S1): 320-328.
[7]	何松, 刘乔迈, 谢广烁, 王斯民, 肖娟. 高浓度水煤浆管道气膜减阻两相流模拟及代理辅助优化[J]. 化工学报, 2023, 74(9): 3766-3774.
[8]	程小松, 殷勇高, 车春文. 不同工质在溶液除湿真空再生系统中的性能对比[J]. 化工学报, 2023, 74(8): 3494-3501.
[9]	刘文竹, 云和明, 王宝雪, 胡明哲, 仲崇龙. 基于场协同和耗散的微通道拓扑优化研究[J]. 化工学报, 2023, 74(8): 3329-3341.
[10]	洪瑞, 袁宝强, 杜文静. 垂直上升管内超临界二氧化碳传热恶化机理分析[J]. 化工学报, 2023, 74(8): 3309-3319.
[11]	邢雷, 苗春雨, 蒋明虎, 赵立新, 李新亚. 井下微型气液旋流分离器优化设计与性能分析[J]. 化工学报, 2023, 74(8): 3394-3406.
[12]	韩晨, 司徒友珉, 朱斌, 许建良, 郭晓镭, 刘海峰. 协同处理废液的多喷嘴粉煤气化炉内反应流动研究[J]. 化工学报, 2023, 74(8): 3266-3278.
[13]	吕龙义, 及文博, 韩沐达, 李伟光, 高文芳, 刘晓阳, 孙丽, 王鹏飞, 任芝军, 张光明. 铁基导电材料强化厌氧去除卤代有机污染物：研究进展及未来展望[J]. 化工学报, 2023, 74(8): 3193-3202.
[14]	黄可欣, 李彤, 李桉琦, 林梅. 加装旋转叶轮T型通道流场的模态分解[J]. 化工学报, 2023, 74(7): 2848-2857.
[15]	史方哲, 甘云华. 超薄热管启动特性和传热性能数值模拟[J]. 化工学报, 2023, 74(7): 2814-2823.