储能式低温空分系统热力学与经济性分析

doi:10.11949/0438-1157.20240255

摘要/Abstract

摘要：

低温空气分离设备是大型化工系统的高耗能环节，若与液态空气储能技术相结合，可有效平衡电网峰谷负荷，显著提高系统的经济效益。提出了一种储能式低温空分系统（ASU-ESG），利用谷电制取液态空气并储存，在峰电期内液态空气膨胀发电并参与低温精馏。此方案可有效降低系统在用电峰时的压缩负荷和电力成本，改善空分装置的负荷调节能力。研究结果表明，空气处理量为60000 m³/h的ASU-ESG系统的日均制氧压缩功比功耗为0.378 kW·h/m³，储能和释能过程氧提取率分别为89.46%和93.71%。相比常规ASU系统，该系统的谷时压缩负荷提高了28%，峰时压缩负荷降低了20%，可节约年度用电成本5.72%~6.88%，动态回收期仅为3.3~4.3年，全生命周期净现值可达（19714.6~28074.7）万元。本系统可平衡电网峰谷波动，并利用电价差显著提高系统的经济效益。

关键词: 空气分离, 液态空气储能, 动态回收期

Abstract:

Low-temperature air separation equipment is a high energy consumption link for large chemical systems. If it is combined with liquid air energy storage technology, it can effectively balance the load of the grid peak valley and significantly improve the economic benefits of the system. This paper introduces an air separation unit with energy storage and generation (ASU-ESG). It uses valley electricity to liquefy air and recovers liquid air for electricity generation and air distillation, thus reducing peak compression loads and improving load regulation ability. The results show that for an ASU-ESG system with an air processing capacity of 60000 m³/h, the power consumption for oxygen production is 0.378 kW·h/m³. The oxygen extraction rates for energy storage and energy release processes reach 89.46% and 93.71%, respectively. Compared with conventional ASU, the ASU-ESG improves compression load by 28% during off-peak electricity periods and reduces it by 20% during on-peak electricity periods, resulting in annual electricity cost savings of 5.72%—6.88%. The dynamic payback period for the ASU-ESG system is 3.3—4.3 a, and the net present value (NPV) reaches 19714.6×10⁴—28074.7×10⁴ CNY. The ASU-ESG system can balance the peak-valley fluctuation in power grids as well as obtain economic benefits in ASU operations.

Key words: air separation, liquid air energy storage, dynamic payback period

中图分类号:

TK 02

秦晓巧, 谭宏博, 温娜. 储能式低温空分系统热力学与经济性分析[J]. 化工学报, 2024, 75(7): 2409-2421.

Xiaoqiao QIN, Hongbo TAN, Na WEN. Thermodynamic and economic analysis of air separation unit with energy storage and generation[J]. CIESC Journal, 2024, 75(7): 2409-2421.

图/表 14

图1 常规ASU系统（a）、ASU-ESG系统储能过程（b）和释能过程（c）流程图

Fig.1 Flow diagram of conventional ASU system (a), ASU-ESG system energy storage process (b) and energy release process (c)

表1 ASU-ESG系统主要参数

Table 1 Main parameters of ASU-ESG system

对象	参数	数值
纯化空气	温度	298.15 K ^[12-13]
	压力	0.1013 MPa ^[12-13]
氧气产品	纯度	≥99.6%（摩尔分数） ^[14]
氮气产品	纯度	≥99.999%（摩尔分数）
CAC	绝热效率	80%^[15]
	压比	1.4
PAC	绝热效率	80%
	压比	1.4
BET	增压端绝热效率	75%
	膨胀端等熵效率	78%
	膨胀端出口带液量	≤3%
AT	等熵效率	80%
	排气压力	0.55 MPa
LAT	容积	120 m³
	压力	0.55 MPa
LAP	绝热效率	80% ^[16]
	泵后压力	6 MPa
CACC/PACC/BETC	进口水温	305.15 K ^[17]
CON	最小换热温差	≥1 K ^[18]
MHX/SC	最小换热温差	≥0.5 K
HPC	压力	0.55 MPa
	塔板数	30
LPC	压力	0.13 MPa ^[19]
	塔板数	50

表2 上海工业用电的峰谷时段［3］

Table 2 Peak and valley electricity periods in Shanghai industrials［3］

项目	低谷时段	平时段	高峰时段
时间范围	23∶00~7∶00	7∶00~10∶00，15∶00~18∶00，21∶00~23∶00	10∶00~15∶00，18∶00~21∶00
ASU-ESG系统运行模式	储能模式	常规模式	释能模式
ASU-ESG系统每日运行时长	8 h	8 h	8 h

图2 常规ASU系统与ASU-ESG系统运行模式比较

Fig.2 Comparison of operation modes between conventional ASU system and ASU-ESG system

表3 常规ASU系统和ASU-ESG系统精馏产品模拟结果

Table 3 Distillation product simulation results of conventional ASU and ASU-ESG systems

参数	常规ASU系统	ASU-ESG 系统储能过程	ASU-ESG 系统释能过程
氧气产品
产量/（m³/h）	9649	9650	9649
纯度/% （摩尔分数）	99.60	99.60	99.60
提取率/%	91.53	89.46	93.71
氮气产品
产量/（m³/h）	21910	21910	21910
纯度/% （摩尔分数）	99.999	99.999	99.999
提取率/%	55.47	52.80	58.42
压缩机功耗/kW	3548	4542	2838
液态空气泵功耗/kW	—	—	24
空气膨胀机输出功率/kW	—	—	190
日均制氧压缩功比功耗/（kW·h/m³）	0.368	0.471	0.294

图3 常规ASU系统(a)、ASU-ESG系统储能过程(b)和释能过程(c)㶲流图

Fig.3 Exergy flow in conventional ASU system (a), ASU-ESG system energy storage process (b) and energy release process (c)

图4 常规ASU系统(a)和ASU-ESG系统(b)初始投资情况

Fig.4 Initial investment costs in conventional ASU system (a) and ASU-ESG system (b)

表4 电价假设条件

Table 4 Calculation of electricity price assumptions

序号	峰电单价 $C ˙ e l e, p$ / （CNY/（kW·h））	谷电单价 $C ˙ e l e, v$ / （CNY/（kW·h））	平电单价 $C ˙ e l e, f$ / （CNY/（kW·h））	每月容量电价 $C ˙ c a p a c i t y$ / (CNY/（kV·A））	峰谷价差/ （CNY/（kW·h））
1	1.177	0.329	0.707	20	0.848
2	1.127	0.315	0.677	20	0.812
3	1.077	0.301	0.647	20	0.776
4	1.027	0.287	0.617	20	0.740
5	0.977	0.273	0.587	20	0.704
6	0.928	0.259	0.557	20	0.668
7	0.878	0.245	0.527	20	0.632
8	0.828	0.231	0.497	20	0.596
9	0.778	0.217	0.467	20	0.560
10	0.728	0.203	0.437	20	0.524

表4 电价假设条件

Table 4 Calculation of electricity price assumptions

序号	峰电单价 $C ˙ e l e, p$ / （CNY/（kW·h））	谷电单价 $C ˙ e l e, v$ / （CNY/（kW·h））	平电单价 $C ˙ e l e, f$ / （CNY/（kW·h））	每月容量电价 $C ˙ c a p a c i t y$ / (CNY/（kV·A））	峰谷价差/ （CNY/（kW·h））
1	1.177	0.329	0.707	20	0.848
2	1.127	0.315	0.677	20	0.812
3	1.077	0.301	0.647	20	0.776
4	1.027	0.287	0.617	20	0.740
5	0.977	0.273	0.587	20	0.704
6	0.928	0.259	0.557	20	0.668
7	0.878	0.245	0.527	20	0.632
8	0.828	0.231	0.497	20	0.596
9	0.778	0.217	0.467	20	0.560
10	0.728	0.203	0.437	20	0.524

图5 峰谷电价差对年度用电成本以及电力成本节约率的影响

Fig.5 Impact of peak and valley electricity price differences on Cele and ξ

图6 峰谷电价差对年度现金流量(a)和累计净现值(b)的影响

Fig.6 Impact of peak and valley electricity price differences on NCF (a) and NPV (b)

图7 峰谷电价差对动态回收期的影响

Fig.7 Impact of peak and valley electricity price differences on DPP

表A1 ASU-ESG系统设备的能量平衡方程和㶲平衡方程

Table A1 Energy balance equations and exergy balance equations for ASU-ESG system

设备	能量平衡方程	㶲平衡方程
CAC	$m ˙ 0 (h 0 - h 1) + P C A C = 0$	$m ˙ 0 (e 0 - e 1) + P C A C = E x d, C A C$
PAC	$m ˙ 33 (h 33 - h 34) + P P A C = 0$	$m ˙ 33 (e 33 - e 34) + P P A C = E x d, P A C$
BET	$m ˙ 3 (h 3 - h 4) + m ˙ 5 (h 5 - h 6) + m ˙ 8 (h 8 - h 9) + m ˙ 11 (h 11 - h 12) = 0$	$m ˙ 3 (e 3 - e 4) + m ˙ 5 (e 5 - e 6) + m ˙ 8 (e 8 - e 9) + m ˙ 11 (e 11 - e 12) = E x d, B E T$
	$m ˙ 3 (h 3 - h 4) + m ˙ 5 (h 5 - h 6) + m ˙ 8 (h 8 - h 48) + m ˙ 11 (h 11 - h 50) = 0$	$m ˙ 3 (e 3 - e 4) + m ˙ 5 (e 5 - e 6) + m ˙ 8 (e 8 - e 48) + m ˙ 11 (e 11 - e 50) = E x d, B E T$
CACC	$m ˙ W 1 (h W 1 - h W 2) + m ˙ 1 (h 1 - h 2) = 0$	$m ˙ 1 (e 1 - e 2) = E x d, C A C C$
PACC	$m ˙ W 7 (h W 7 - h W 8) + m ˙ 34 (h 34 - h 35) = 0$	$m ˙ 34 (e 34 - e 35) = E x d, P A C C$
BETC	$m ˙ W 3 (h W 3 - h W 4) + m ˙ W 5 (h W 5 - h W 6) + m ˙ 4 (h 4 - h 5) + m ˙ 6 (h 6 - h 7) = 0$	$m ˙ 4 (e 4 - e 5) + m ˙ 6 (e 6 - e 7) = E x d, B E T C$
MHX	$m ˙ 7 (h 7 - h 8) + m ˙ 10 (h 10 - h 11) + m ˙ 14 (h 14 - h 15) + m ˙ 25 (h 25 - h 26) + m ˙ 28 (h 28 - h 29) + m ˙ 31 (h 31 - h 32) = 0$	$m ˙ 7 (e 7 - e 8) + m ˙ 10 (e 10 - e 11) + m ˙ 14 (e 14 - e 15) + m ˙ 25 (e 25 - e 26) + m ˙ 28 (e 28 - e 29) + m ˙ 31 (e 31 - e 32) = E x d, M H X$
	$m ˙ 7 (h 7 - h 8) + m ˙ 10 (h 10 - h 11) + m ˙ 14 (h 14 - h 15) + m ˙ 35 (h 35 - h 36) + m ˙ 53 (h 53 - h 26) + m ˙ 57 (h 57 - h 29) + m ˙ 55 (h 55 - h 32) = 0$	$m ˙ 7 (e 7 - e 8) + m ˙ 10 (e 10 - e 11) + m ˙ 14 (e 14 - e 15) + m ˙ 35 (e 35 - e 36) + m ˙ 53 (e 53 - e 26) + m ˙ 57 (e 57 - e 29) + m ˙ 55 (e 55 - e 32) = E x d, M H X$
	$m ˙ 7 (h 7 - h 8) + m ˙ 10 (h 10 - h 11) + m ˙ 14 (h 14 - h 15) + m ˙ 25 (h 25 - h 26) + m ˙ 28 (h 28 - h 29) + m ˙ 31 (h 31 - h 32) + m ˙ 41 (h 41 - h 42) + m ˙ 43 (h 43 - h 44) = 0$	$m ˙ 7 (e 7 - e 8) + m ˙ 10 (e 10 - e 11) + m ˙ 14 (e 14 - e 15) + m ˙ 25 (e 25 - e 26) + m ˙ 28 (e 28 - e 29) + m ˙ 31 (e 31 - e 32) + m ˙ 41 (e 41 - e 42) + m ˙ 43 (e 43 - e 44) = E x d, M H X$
CON	$m ˙ 36 (h 36 - h 37) + m ˙ 46 (h 46 - h 47) + m ˙ 48 (h 48 - h 49) + m ˙ 50 (h 50 - h 51) + m ˙ 52 (h 52 - h 53) + m ˙ 54 (h 54 - h 55) + m ˙ 56 (h 56 - h 57) = 0$	$m ˙ 36 (e 36 - e 37) + m ˙ 46 (e 46 - e 47) + m ˙ 48 (e 48 - e 49) + m ˙ 50 (e 50 - e 51) + m ˙ 52 (e 52 - e 53) + m ˙ 54 (e 54 - e 55) + m ˙ 56 (e 56 - e 57) = E x d, C O N$
SC	$m ˙ 16 (h 16 - h 17) + m ˙ 21 (h 21 - h 22) + m ˙ 24 (h 24 - h 25) + m ˙ 27 (h 27 - h 28) + m ˙ 30 (h 30 - h 31) = 0$	$m ˙ 16 (e 16 - e 17) + m ˙ 21 (e 21 - e 22) + m ˙ 24 (e 24 - e 25) + m ˙ 27 (e 27 - e 28) + m ˙ 30 (e 30 - e 31) = E x d, S C$
	$m ˙ 16 (h 16 - h 17) + m ˙ 21 (h 21 - h 22) + m ˙ 24 (h 24 - h 52) + m ˙ 27 (h 27 - h 56) + m ˙ 30 (h 30 - h 54) = 0$	$m ˙ 16 (e 16 - e 17) + m ˙ 21 (e 21 - e 22) + m ˙ 24 (e 24 - e 52) + m ˙ 27 (e 27 - e 56) + m ˙ 30 (e 30 - e 54) = E x d, S C$
TV	$m ˙ 17 (h 17 - h 18) + m ˙ 22 (h 22 - h 23) + m ˙ 37 (h 37 - h 38) = 0$	$m ˙ 17 (e 17 - e 18) + m ˙ 22 (e 22 - e 23) + m ˙ 37 (e 37 - e 38) = E x d, T V$
SEP	$m ˙ 18 (h 18 - h 19 - h 20) + m ˙ 38 (h 38 - h 39 - h 46) = 0$	$m ˙ 18 (e 18 - e 19 - e 20) + m ˙ 38 (e 38 - e 39 - e 46) = E x d, S E P$
LAT	$m ˙ 39 (h 39 - h 40) = 0$	$m ˙ 39 (e 39 - e 40) = 0$
DC	$m ˙ 13 h 13 + m ˙ 15 h 15 + m ˙ 19 h 19 + m ˙ 20 h 20 + m ˙ 23 h 23 + Q r = m ˙ 16 h 16 + m ˙ 21 h 21 + m ˙ 24 h 24 + m ˙ 27 h 27 + m ˙ 30 h 30 + Q c$	$m ˙ 13 e 13 + m ˙ 15 e 15 + m ˙ 19 e 19 + m ˙ 20 e 20 + m ˙ 23 e 23 = m ˙ 16 e 16 + m ˙ 21 e 21 + m ˙ 24 e 24 + m ˙ 27 e 27 + m ˙ 30 e 30 + E x d, D C$
	$m ˙ 13 h 13 + m ˙ 15 h 15 + m ˙ 19 h 19 + m ˙ 20 h 20 + m ˙ 23 h 23 + m ˙ 47 h 47 + Q r = m ˙ 16 h 16 + m ˙ 21 h 21 + m ˙ 24 h 24 + m ˙ 27 h 27 + m ˙ 30 h 30 + Q c$	$m ˙ 13 e 13 + m ˙ 15 e 15 + m ˙ 19 e 19 + m ˙ 20 e 20 + m ˙ 23 e 23 + m ˙ 45 e 47 = m ˙ 16 e 16 + m ˙ 21 e 21 + m ˙ 24 e 24 + m ˙ 27 e 27 + m ˙ 30 e 30 + E x d, D C$
	$m ˙ 13 h 13 + m ˙ 15 h 15 + m ˙ 19 h 19 + m ˙ 20 h 20 + m ˙ 23 h 23 + m ˙ 45 h 45 + Q r = m ˙ 16 h 16 + m ˙ 21 h 21 + m ˙ 24 h 24 + m ˙ 27 h 27 + m ˙ 30 h 30 + Q c$	$m ˙ 13 e 13 + m ˙ 15 e 15 + m ˙ 19 e 19 + m ˙ 20 e 20 + m ˙ 23 e 23 + m ˙ 45 e 45 = m ˙ 16 e 16 + m ˙ 21 e 21 + m ˙ 24 e 24 + m ˙ 27 e 27 + m ˙ 30 e 30 + E x d, D C$
LAP	$m ˙ 40 (h 40 - h 41) + P L A P = 0$	$m ˙ 40 (e 40 - e 41) + P L A P = E x d, L A P$
AT	$m ˙ 42 (h 42 - h 43) + m ˙ 44 (h 44 - h 45) = P A T 1 + P A T 2$	$m ˙ 42 (e 42 - e 43) + m ˙ 44 (e 44 - e 45) = P A T 1 + P A T 2 + E x d, A T$

表A1 ASU-ESG系统设备的能量平衡方程和㶲平衡方程

Table A1 Energy balance equations and exergy balance equations for ASU-ESG system

设备	能量平衡方程	㶲平衡方程
CAC	$m ˙ 0 (h 0 - h 1) + P C A C = 0$	$m ˙ 0 (e 0 - e 1) + P C A C = E x d, C A C$
PAC	$m ˙ 33 (h 33 - h 34) + P P A C = 0$	$m ˙ 33 (e 33 - e 34) + P P A C = E x d, P A C$
BET	$m ˙ 3 (h 3 - h 4) + m ˙ 5 (h 5 - h 6) + m ˙ 8 (h 8 - h 9) + m ˙ 11 (h 11 - h 12) = 0$	$m ˙ 3 (e 3 - e 4) + m ˙ 5 (e 5 - e 6) + m ˙ 8 (e 8 - e 9) + m ˙ 11 (e 11 - e 12) = E x d, B E T$
	$m ˙ 3 (h 3 - h 4) + m ˙ 5 (h 5 - h 6) + m ˙ 8 (h 8 - h 48) + m ˙ 11 (h 11 - h 50) = 0$	$m ˙ 3 (e 3 - e 4) + m ˙ 5 (e 5 - e 6) + m ˙ 8 (e 8 - e 48) + m ˙ 11 (e 11 - e 50) = E x d, B E T$
CACC	$m ˙ W 1 (h W 1 - h W 2) + m ˙ 1 (h 1 - h 2) = 0$	$m ˙ 1 (e 1 - e 2) = E x d, C A C C$
PACC	$m ˙ W 7 (h W 7 - h W 8) + m ˙ 34 (h 34 - h 35) = 0$	$m ˙ 34 (e 34 - e 35) = E x d, P A C C$
BETC	$m ˙ W 3 (h W 3 - h W 4) + m ˙ W 5 (h W 5 - h W 6) + m ˙ 4 (h 4 - h 5) + m ˙ 6 (h 6 - h 7) = 0$	$m ˙ 4 (e 4 - e 5) + m ˙ 6 (e 6 - e 7) = E x d, B E T C$
MHX	$m ˙ 7 (h 7 - h 8) + m ˙ 10 (h 10 - h 11) + m ˙ 14 (h 14 - h 15) + m ˙ 25 (h 25 - h 26) + m ˙ 28 (h 28 - h 29) + m ˙ 31 (h 31 - h 32) = 0$	$m ˙ 7 (e 7 - e 8) + m ˙ 10 (e 10 - e 11) + m ˙ 14 (e 14 - e 15) + m ˙ 25 (e 25 - e 26) + m ˙ 28 (e 28 - e 29) + m ˙ 31 (e 31 - e 32) = E x d, M H X$
	$m ˙ 7 (h 7 - h 8) + m ˙ 10 (h 10 - h 11) + m ˙ 14 (h 14 - h 15) + m ˙ 35 (h 35 - h 36) + m ˙ 53 (h 53 - h 26) + m ˙ 57 (h 57 - h 29) + m ˙ 55 (h 55 - h 32) = 0$	$m ˙ 7 (e 7 - e 8) + m ˙ 10 (e 10 - e 11) + m ˙ 14 (e 14 - e 15) + m ˙ 35 (e 35 - e 36) + m ˙ 53 (e 53 - e 26) + m ˙ 57 (e 57 - e 29) + m ˙ 55 (e 55 - e 32) = E x d, M H X$
	$m ˙ 7 (h 7 - h 8) + m ˙ 10 (h 10 - h 11) + m ˙ 14 (h 14 - h 15) + m ˙ 25 (h 25 - h 26) + m ˙ 28 (h 28 - h 29) + m ˙ 31 (h 31 - h 32) + m ˙ 41 (h 41 - h 42) + m ˙ 43 (h 43 - h 44) = 0$	$m ˙ 7 (e 7 - e 8) + m ˙ 10 (e 10 - e 11) + m ˙ 14 (e 14 - e 15) + m ˙ 25 (e 25 - e 26) + m ˙ 28 (e 28 - e 29) + m ˙ 31 (e 31 - e 32) + m ˙ 41 (e 41 - e 42) + m ˙ 43 (e 43 - e 44) = E x d, M H X$
CON	$m ˙ 36 (h 36 - h 37) + m ˙ 46 (h 46 - h 47) + m ˙ 48 (h 48 - h 49) + m ˙ 50 (h 50 - h 51) + m ˙ 52 (h 52 - h 53) + m ˙ 54 (h 54 - h 55) + m ˙ 56 (h 56 - h 57) = 0$	$m ˙ 36 (e 36 - e 37) + m ˙ 46 (e 46 - e 47) + m ˙ 48 (e 48 - e 49) + m ˙ 50 (e 50 - e 51) + m ˙ 52 (e 52 - e 53) + m ˙ 54 (e 54 - e 55) + m ˙ 56 (e 56 - e 57) = E x d, C O N$
SC	$m ˙ 16 (h 16 - h 17) + m ˙ 21 (h 21 - h 22) + m ˙ 24 (h 24 - h 25) + m ˙ 27 (h 27 - h 28) + m ˙ 30 (h 30 - h 31) = 0$	$m ˙ 16 (e 16 - e 17) + m ˙ 21 (e 21 - e 22) + m ˙ 24 (e 24 - e 25) + m ˙ 27 (e 27 - e 28) + m ˙ 30 (e 30 - e 31) = E x d, S C$
	$m ˙ 16 (h 16 - h 17) + m ˙ 21 (h 21 - h 22) + m ˙ 24 (h 24 - h 52) + m ˙ 27 (h 27 - h 56) + m ˙ 30 (h 30 - h 54) = 0$	$m ˙ 16 (e 16 - e 17) + m ˙ 21 (e 21 - e 22) + m ˙ 24 (e 24 - e 52) + m ˙ 27 (e 27 - e 56) + m ˙ 30 (e 30 - e 54) = E x d, S C$
TV	$m ˙ 17 (h 17 - h 18) + m ˙ 22 (h 22 - h 23) + m ˙ 37 (h 37 - h 38) = 0$	$m ˙ 17 (e 17 - e 18) + m ˙ 22 (e 22 - e 23) + m ˙ 37 (e 37 - e 38) = E x d, T V$
SEP	$m ˙ 18 (h 18 - h 19 - h 20) + m ˙ 38 (h 38 - h 39 - h 46) = 0$	$m ˙ 18 (e 18 - e 19 - e 20) + m ˙ 38 (e 38 - e 39 - e 46) = E x d, S E P$
LAT	$m ˙ 39 (h 39 - h 40) = 0$	$m ˙ 39 (e 39 - e 40) = 0$
DC	$m ˙ 13 h 13 + m ˙ 15 h 15 + m ˙ 19 h 19 + m ˙ 20 h 20 + m ˙ 23 h 23 + Q r = m ˙ 16 h 16 + m ˙ 21 h 21 + m ˙ 24 h 24 + m ˙ 27 h 27 + m ˙ 30 h 30 + Q c$	$m ˙ 13 e 13 + m ˙ 15 e 15 + m ˙ 19 e 19 + m ˙ 20 e 20 + m ˙ 23 e 23 = m ˙ 16 e 16 + m ˙ 21 e 21 + m ˙ 24 e 24 + m ˙ 27 e 27 + m ˙ 30 e 30 + E x d, D C$
	$m ˙ 13 h 13 + m ˙ 15 h 15 + m ˙ 19 h 19 + m ˙ 20 h 20 + m ˙ 23 h 23 + m ˙ 47 h 47 + Q r = m ˙ 16 h 16 + m ˙ 21 h 21 + m ˙ 24 h 24 + m ˙ 27 h 27 + m ˙ 30 h 30 + Q c$	$m ˙ 13 e 13 + m ˙ 15 e 15 + m ˙ 19 e 19 + m ˙ 20 e 20 + m ˙ 23 e 23 + m ˙ 45 e 47 = m ˙ 16 e 16 + m ˙ 21 e 21 + m ˙ 24 e 24 + m ˙ 27 e 27 + m ˙ 30 e 30 + E x d, D C$
	$m ˙ 13 h 13 + m ˙ 15 h 15 + m ˙ 19 h 19 + m ˙ 20 h 20 + m ˙ 23 h 23 + m ˙ 45 h 45 + Q r = m ˙ 16 h 16 + m ˙ 21 h 21 + m ˙ 24 h 24 + m ˙ 27 h 27 + m ˙ 30 h 30 + Q c$	$m ˙ 13 e 13 + m ˙ 15 e 15 + m ˙ 19 e 19 + m ˙ 20 e 20 + m ˙ 23 e 23 + m ˙ 45 e 45 = m ˙ 16 e 16 + m ˙ 21 e 21 + m ˙ 24 e 24 + m ˙ 27 e 27 + m ˙ 30 e 30 + E x d, D C$
LAP	$m ˙ 40 (h 40 - h 41) + P L A P = 0$	$m ˙ 40 (e 40 - e 41) + P L A P = E x d, L A P$
AT	$m ˙ 42 (h 42 - h 43) + m ˙ 44 (h 44 - h 45) = P A T 1 + P A T 2$	$m ˙ 42 (e 42 - e 43) + m ˙ 44 (e 44 - e 45) = P A T 1 + P A T 2 + E x d, A T$

表A2 ASU-ESG系统设备的成本估算模型

Table A2 Equipment investment cost calculation formulas for the ASU-ESG system

项目	设备	成本估算模型
预处理和纯化单元	PPU	$C P P U = 13 % × C e q, A S U$
涡轮机械	CAC	$C C A C = 1.051 × 174.9943 × p C A C ÷ 0.14$ ^[32]
	PAC	$l g (0.14 × C P A C) = K C E P C I K B M {K 1 + K 2 l g (P P A C / 1000) + K 3 [l g (p P A C / 1000)] 2}$ ^[24] K_CEPCI=CEPCI₂₀₂₃/CEPCI₂₀₀₁=708/397 ^[19] K_BM=2.75，K₁=2.2897，K₂=1.3604，K₃=-0.1027
	BET	$C B E T = C B E T - A C + C B E T - A T$
		$C B E T - A C = 174.9943 × p B E T - A C ÷ 0.14$ ^[13]
		$C B E T - A T = 1100 × (p B E T - A T) 0.81 ÷ 0.14$ ^[33]
换热设备	CACC/PACC/BETC	$C C A C C = 12000 × (U A C A C C / 100) 0.6 ÷ 0.14$ ^[34] $U A C A C C = (Q ˙ / L M T D) C A C C$
精馏单元	SC	$C S C = 12000 × (U A S C / 100) 0.6 ÷ 0.14$ ^[34]
	MHX	$C M H X = 1.218 × F m F p e x p {8.821 - 0.65 × l n (10.7639 × A M H X) + 0.0681 ×$ $[l n (10.7639 × A M H X)] 2} ÷ 0.14$ ^[35] $F m = 2.3$ ， $F p = 1.2$ ^[36] $U M H X = 70 W / m 2$ ^[35]
	TV1/TV2	$C T V = 114.5 × (m ˙ T V) 0.67 ÷ 0.14$ ^[34]
	SEP1	$C S E P = 1218 × (65 + 3.5 × m ˙ S E P × 3600 × 10 - 3) ÷ 0.14$ ^[8]
	HPC/LPC	$C H P C = C v e s s e l + C t r a y + C h e c$ ^[37-38]
		$C v e s s e l = (M & S / 280) F c, v e s s e l D 1.066 H 0.802 ÷ 0.14$ ^[39]
		$C t r a y = (M & S / 280) F c, t r a y D 1.55 H ÷ 0.14$ ^[39]
		$C h e c = (M & S / 280) F c, h e c A h e c 0.65 ÷ 0.14$ ^[39]
		$M & S = 1490.2$ ， $F c, v e s s e l = 3966.2$ ， $F c, t r a y = 437.4$ ， $F c, H E C = 1609.16$ ^[39]
LAES单元	CON	$C C O N = 32800 × (U A C O N / 80) 0.68 ÷ 0.14$ ^[34]
	LAT	$C L A T = 320 × V L A T ÷ 0.14$ ^[14]
	LAP	$C L A P = 483 × p L A P ÷ 0.14$ ^[14]
	AT	$C A T = 1100 × p A T 0.81 ÷ 0.14$ ^[33]
	TV3	$C T V = 114.5 × m ˙ T V 0.67 ÷ 0.14$ ^[34]
	SEP2	$C S E P = 1218 × (65 + 3.5 × m ˙ S E P × 3600 × 10 - 3) ÷ 0.14$ ^[8]

表A2 ASU-ESG系统设备的成本估算模型

Table A2 Equipment investment cost calculation formulas for the ASU-ESG system

项目	设备	成本估算模型
预处理和纯化单元	PPU	$C P P U = 13 % × C e q, A S U$
涡轮机械	CAC	$C C A C = 1.051 × 174.9943 × p C A C ÷ 0.14$ ^[32]
	PAC	$l g (0.14 × C P A C) = K C E P C I K B M {K 1 + K 2 l g (P P A C / 1000) + K 3 [l g (p P A C / 1000)] 2}$ ^[24] K_CEPCI=CEPCI₂₀₂₃/CEPCI₂₀₀₁=708/397 ^[19] K_BM=2.75，K₁=2.2897，K₂=1.3604，K₃=-0.1027
	BET	$C B E T = C B E T - A C + C B E T - A T$
		$C B E T - A C = 174.9943 × p B E T - A C ÷ 0.14$ ^[13]
		$C B E T - A T = 1100 × (p B E T - A T) 0.81 ÷ 0.14$ ^[33]
换热设备	CACC/PACC/BETC	$C C A C C = 12000 × (U A C A C C / 100) 0.6 ÷ 0.14$ ^[34] $U A C A C C = (Q ˙ / L M T D) C A C C$
精馏单元	SC	$C S C = 12000 × (U A S C / 100) 0.6 ÷ 0.14$ ^[34]
	MHX	$C M H X = 1.218 × F m F p e x p {8.821 - 0.65 × l n (10.7639 × A M H X) + 0.0681 ×$ $[l n (10.7639 × A M H X)] 2} ÷ 0.14$ ^[35] $F m = 2.3$ ， $F p = 1.2$ ^[36] $U M H X = 70 W / m 2$ ^[35]
	TV1/TV2	$C T V = 114.5 × (m ˙ T V) 0.67 ÷ 0.14$ ^[34]
	SEP1	$C S E P = 1218 × (65 + 3.5 × m ˙ S E P × 3600 × 10 - 3) ÷ 0.14$ ^[8]
	HPC/LPC	$C H P C = C v e s s e l + C t r a y + C h e c$ ^[37-38]
		$C v e s s e l = (M & S / 280) F c, v e s s e l D 1.066 H 0.802 ÷ 0.14$ ^[39]
		$C t r a y = (M & S / 280) F c, t r a y D 1.55 H ÷ 0.14$ ^[39]
		$C h e c = (M & S / 280) F c, h e c A h e c 0.65 ÷ 0.14$ ^[39]
		$M & S = 1490.2$ ， $F c, v e s s e l = 3966.2$ ， $F c, t r a y = 437.4$ ， $F c, H E C = 1609.16$ ^[39]
LAES单元	CON	$C C O N = 32800 × (U A C O N / 80) 0.68 ÷ 0.14$ ^[34]
	LAT	$C L A T = 320 × V L A T ÷ 0.14$ ^[14]
	LAP	$C L A P = 483 × p L A P ÷ 0.14$ ^[14]
	AT	$C A T = 1100 × p A T 0.81 ÷ 0.14$ ^[33]
	TV3	$C T V = 114.5 × m ˙ T V 0.67 ÷ 0.14$ ^[34]
	SEP2	$C S E P = 1218 × (65 + 3.5 × m ˙ S E P × 3600 × 10 - 3) ÷ 0.14$ ^[8]

表A3 资本支出和年度运营支出的计算模型

Table A3 Capital investment and operating cost models

项目	计算模型
系统设备（ $C e q$ ）	见表A2
安装工程（ $C i n s t a l l$ ）	$C i n s t a l l = 10 % × C e q$ ^[24]
仪控电控系统（ $C e e m$ ）	$C e e m = 18 % × C e q$ ^{[15, 40]}
管道（ $C p i p$ ）	$C p i p = 10 % × C e q$ ^[40]
其他（ $C s f$ ）	$C s f = 30 % × C e q$ ^[40]
直接支出（ $C D P C$ ）	$C D P C = C e q + C i n s t a l l + C e e m + C p i p + C s f$
采购（ $C e p$ ）	$C e p = 9.8 % × C e q$ ^[40]
施工（ $C c c f$ ）	$C c c f = 11.9 % × C D P C$ ^[40]
间接支出（ $C I P C$ ）	$C I P C = C e p + C c c f$
资本支出总额（ $C T C I$ ）	$C T C I = C D P C + C I P C$
运营人工（ $C o l a b$ ）	$C o l a b = 0.3 % × 20 × 8000 × 12$
维修维护（ $C m a i n t a i n$ ）	$C m a i n t a i n = 2.5 % × 95 % × C T C I$ ^[24]
保险（ $C i n s u r$ ）	$C i n s u r = 1 % × C T C I$ ^[15]
年度固定运营支出（ $C F O C$ ）	$C F O C = C o l a b + C m a i n t a i n + C i n s u r$
电力（ $C e l e$ ）	见正文式（6）和式（7）
冷却用水（ $C c w$ ）	$C c w = 340 × 10 - 3 × 8 × 3600 × Z c w (m ˙ c w, p + m ˙ c w, v + m ˙ c w, f)$
耗材（ $C c o n s$ ）	$C c o n s = 1 % × C e q$ ^[35]
年度可变运营支出（ $C V O C$ ）	$C V O C = C e l e + C c w + C c o n s$
年度运营支出总额（ $C O & M$ ）	$C O & M = C F O C + C V O C$

表A3 资本支出和年度运营支出的计算模型

Table A3 Capital investment and operating cost models

项目	计算模型
系统设备（ $C e q$ ）	见表A2
安装工程（ $C i n s t a l l$ ）	$C i n s t a l l = 10 % × C e q$ ^[24]
仪控电控系统（ $C e e m$ ）	$C e e m = 18 % × C e q$ ^{[15, 40]}
管道（ $C p i p$ ）	$C p i p = 10 % × C e q$ ^[40]
其他（ $C s f$ ）	$C s f = 30 % × C e q$ ^[40]
直接支出（ $C D P C$ ）	$C D P C = C e q + C i n s t a l l + C e e m + C p i p + C s f$
采购（ $C e p$ ）	$C e p = 9.8 % × C e q$ ^[40]
施工（ $C c c f$ ）	$C c c f = 11.9 % × C D P C$ ^[40]
间接支出（ $C I P C$ ）	$C I P C = C e p + C c c f$
资本支出总额（ $C T C I$ ）	$C T C I = C D P C + C I P C$
运营人工（ $C o l a b$ ）	$C o l a b = 0.3 % × 20 × 8000 × 12$
维修维护（ $C m a i n t a i n$ ）	$C m a i n t a i n = 2.5 % × 95 % × C T C I$ ^[24]
保险（ $C i n s u r$ ）	$C i n s u r = 1 % × C T C I$ ^[15]
年度固定运营支出（ $C F O C$ ）	$C F O C = C o l a b + C m a i n t a i n + C i n s u r$
电力（ $C e l e$ ）	见正文式（6）和式（7）
冷却用水（ $C c w$ ）	$C c w = 340 × 10 - 3 × 8 × 3600 × Z c w (m ˙ c w, p + m ˙ c w, v + m ˙ c w, f)$
耗材（ $C c o n s$ ）	$C c o n s = 1 % × C e q$ ^[35]
年度可变运营支出（ $C V O C$ ）	$C V O C = C e l e + C c w + C c o n s$
年度运营支出总额（ $C O & M$ ）	$C O & M = C F O C + C V O C$

参考文献 40

1	才正彬. 空分装置高效节能优化与应用[D]. 杭州: 浙江大学, 2019.
	Cai Z B. High efficiency and energy saving optimization and application of air separation unit[D]. Hangzhou: Zhejiang University, 2019.
2	张培昆, 王立. 空分短期停车时间阈值对氧气生产调度的影响[J]. 化工学报, 2017, 68(6): 2423-2433.
	Zhang P K, Wang L. Effects of temporary shutdown time-threshold on oxygen production schedule in air separation unit[J]. CIESC Journal, 2017, 68(6): 2423-2433.
3	Liu Y N, Wang L, He X F. An external-compression air separation unit with energy storage and its thermodynamic and economic analysis[J]. Journal of Energy Storage, 2023, 59: 106513.
4	王晨, 折晓会, 张小松. 含空气净化过程的液态空气储能热力学研究[J]. 化工学报, 2020, 71(S1): 23-30.
	Wang C, She X H, Zhang X S. Thermodynamic study of liquid air energy storage with air purification unit[J]. CIESC Journal, 2020, 71(S1): 23-30.
5	Wang K W, Tong L G, Yin S W, et al. Novel ASU-LAES system with flexible energy release: analysis of cycle performance, economics, and peak shaving advantages[J]. Energy, 2024, 288: 129720.
6	Liu Y X, Kong F L, Tong L G, et al. A novel cryogenic air separation unit with energy storage: recovering waste heat and reusing storage media[J]. Journal of Energy Storage, 2024, 80: 110359.
7	He X F, Liu Y N, Rehman A, et al. A novel air separation unit with energy storage and generation and its energy efficiency and economy analysis[J]. Applied Energy, 2021, 281: 115976.
8	He X F, Liu Y N, Rehman A, et al. Feasibility and performance analysis of a novel air separation unit with energy storage and air recovery[J]. Renewable Energy, 2022, 195: 598-619.
9	Kong F L, Liu Y X, Shen M H, et al. A novel economic scheduling of multi-product deterministic demand for co-production air separation system with liquid air energy storage[J]. Renewable Energy, 2023, 209: 533-545.
10	黄欢. 耦合空分装置生产的液化空气储能系统研究[J]. 当代化工研究, 2022(16): 101-103.
	Huang H. Research on the liquefied air energy storage system coupled with air separation unit production[J]. Modern Chemical Research, 2022(16): 101-103.
11	刘波. 基于HYSYS的空分精馏过程仿真分析[D]. 西安: 西安电子科技大学, 2013.
	Liu B. Distillation process simulation analysis based on HYSYS[D]. Xi’an: Xidian University, 2013.
12	Li D, Duan L Q. Techno-economic analysis of a new thermal storage operation strategy for a solar aided liquid air energy storage system[J]. Journal of Energy Storage, 2024, 78: 110029.
13	Adedeji M, Abid M, Dagbasi M, et al. Improvement of a liquid air energy storage system: investigation of performance analysis for novel ambient air conditioning[J]. Journal of Energy Storage, 2022, 50: 104294.
14	初日召, 李勇, 王庄. KDON-5600/14000空分装置运行总结[J]. 中国氯碱, 2017(8): 28-30, 47.
	Chu R Z, Li Y, Wang Z. Introduction and operation summary of KDON-5600/14000 air separation unit[J]. China Chlor-Alkali, 2017(8): 28-30, 47.
15	Huo C J, Sun J J, Song P. Energy, exergy and economic analyses of an optimal use of cryogenic liquid turbine expander in air separation units[J]. Chemical Engineering Research and Design, 2023, 189: 194-209.
16	Li Y H, Fan X Y, Li J X, et al. Novel liquid air energy storage coupled with liquefied ethylene cold energy: thermodynamic, exergy and economic analysis[J]. Applied Thermal Engineering, 2024, 245: 122909.
17	荣杨一鸣, 吴巧仙, 周霞, 等. 空分系统空气压缩余热自利用性能优化研究[J]. 化工学报, 2021, 72(3): 1654-1666.
	Rong Y Y M, Wu Q X, Zhou X, et al. Research on optimization of self-utilization performance of air compression waste heat in air separation system[J]. CIESC Journal, 2021, 72(3): 1654-1666.
18	Yang Y, Tong L G, Liu Y X, et al. A novel integrated system of hydrogen liquefaction process and liquid air energy storage (LAES): energy, exergy, and economic analysis[J]. Energy Conversion and Management, 2023, 280: 116799.
19	Tesch S, Morosuk T, Tsatsaronis G. Comparative evaluation of cryogenic air separation units from the exergetic and economic points of view[M]//Tatiana M, Muhammad S. Low-temperature Technologies. Rijeka: IntechOpen, 2019.
20	Fan X Y, Ji W, Guo L N, et al. Thermo-economic analysis of the integrated system of thermal power plant and liquid air energy storage[J]. Journal of Energy Storage, 2023, 57: 106233.
21	李燕鹏. 低温空气分离装置的流程选型方法研究[D]. 杭州: 浙江大学, 2019.
	Li Y P. Research on process selection method of cryogenic air separation unit[D]. Hangzhou: Zhejiang University, 2019.
22	盛新江. 空分设备项目投资经济性分析[J]. 深冷技术, 2005(4): 5-7.
	Sheng X J. Economic analysis of investing an ASU project[J]. Cryogenic Technology, 2005(4): 5-7.
23	中华人民共和国个人所得税法实施条例[J]. 中华人民共和国国务院公报, 2019(1): 25-30.
	Regulations for the implementation of the individual income tax law of the People’s Republic of China[J]. Gazette of the State Council of the People’s Republic of China, 2019(1): 25-30.
24	田英男. 先进绝热压缩空气储能系统的高等㶲分析及经济性评价[D]. 济南:山东大学, 2023.
	Tian Y N. Advanced exergy analysis and economic evaluation of advanced adiabatic compressed air energy storage system[D]. Jinan: Shandong University, 2023.
25	周程, 苏礼勇. 40000 m³/h空分装置设备选型及优化的可行性分析[J]. 中氮肥, 2020(3): 1-6.
	Zhou C, Su L Y. Feasibility analysis of equipment selection and optimization of 40000 m³/h air separation unit[J]. Nitrogenous Fertilizer Progress, 2020(3): 1-6.
26	李政辰. 整装冷箱形式特大型空分设备成本管理研究[D]. 杭州: 浙江大学, 2023.
	Li Z C. Study of cost control of large size ASU with modular coldbox[D]. Hangzhou: Zhejiang University, 2023.
27	Li J X, Fan X Y, Li Y H, et al. A novel system of liquid air energy storage with LNG cold energy and industrial waste heat: thermodynamic and economic analysis[J]. Journal of Energy Storage, 2024, 86: 111359.
28	Gao Z Z, Ji W, Guo L N, et al. Thermo-economic analysis of the integrated bidirectional peak shaving system consisted by liquid air energy storage and combined cycle power plant[J]. Energy Conversion and Management, 2021, 234: 113945.
29	Liu Y X, Tong L G, Kong F L, et al. An improved ASU distillation process and DIM-LPB method for variable product ratio demand[J]. Separation and Purification Technology, 2021, 277: 119499.
30	陈仕卿. LNG冷能在空气分离系统中的集成与优化研究[D]. 北京: 中国科学院大学（中国科学院工程热物理研究所）, 2019.
	Chen S Q. Research on the integration of LNG regasification and air separation units[D]. Beijing: Institute of Engineering Thermophysics, Chinese Academy of Sciences, 2019.
31	马国光, 李雅娴, 张晨. 基于改良㶲分析方法的LNG冷能空分工艺优化[J]. 天然气工业, 2018, 38(9): 121-128.
	Ma G G, Li Y X, Zhang C. Optimization of the LNG cold energy air separation process based on the advanced exergy analysis method[J]. Natural Gas Industry, 2018, 38(9): 121-128.
32	李佳佳, 李兴朔, 魏凡超, 等. 耦合火电机组的新型压缩空气储能系统技术经济性评估研究[J]. 中国电机工程学报, 2023, 43(23): 9171-9183.
	Li J J, Li X S, Wei F C, et al. Research on techno-economic evaluation of new type compressed air energy storage coupled with thermal power unit[J]. Proceedings of the CSEE, 2023, 43(23): 9171-9183.
33	Bashiri Mousavi S, Ahmadi P, Hanafizadeh P, et al. Dynamic simulation and techno-economic analysis of liquid air energy storage with cascade phase change materials as a cold storage system[J]. Journal of Energy Storage, 2022, 50: 104179.
34	Saleh Kandezi M, Mousavi Naeenian S M. Investigation of an efficient and green system based on liquid air energy storage (LAES) for district cooling and peak shaving: energy and exergy analyses[J]. Sustainable Energy Technologies and Assessments, 2021, 47:101396.
35	Singla R, Chowdhury K. Comparisons of thermodynamic and economic performances of cryogenic air separation plants designed for external and internal compression of oxygen[J]. Applied Thermal Engineering, 2019, 160: 114025.
36	Hamdy S, Morosuk T, Tsatsaronis G. Exergetic and economic assessment of integrated cryogenic energy storage systems[J]. Cryogenics, 2019, 99: 39-50.
37	Ebrahimi A, Ziabasharhagh M. Optimal design and integration of a cryogenic air separation unit (ASU) with liquefied natural gas (LNG) as heat sink, thermodynamic and economic analyses[J]. Energy, 2017, 126: 868-885.
38	Mehrpooya M, Zonouz M J. Analysis of an integrated cryogenic air separation unit, oxy-combustion carbon dioxide power cycle and liquefied natural gas regasification process by exergoeconomic method[J]. Energy Conversion and Management, 2017, 139: 245-259.
39	王之宇. 热耦合空分塔的节能优化与最优结构设计研究[D]. 杭州: 浙江大学, 2021.
	Wang Z Y. Research on energy saving optimization and optimal structure design of heat-intrgrated air separation columns[D]. Hangzhou: Zhejiang University, 2021.
40	Park J H, Heo J Y, Lee J I. Techno-economic study of nuclear integrated liquid air energy storage system[J]. Energy Conversion and Management, 2022, 251: 114937.

[1]	文兆伦, 李沛睿, 张忠林, 杜晓, 侯起旺, 刘叶刚, 郝晓刚, 官国清. 基于自热再生的隔壁塔深冷空分工艺设计及优化[J]. 化工学报, 2023, 74(7): 2988-2998.
[2]	王晨, 折晓会, 张小松. 含空气净化过程的液态空气储能热力学研究[J]. 化工学报, 2020, 71(S1): 23-30.
[3]	田涛, 刘冰, 石梅生, 安亚雄, 马军, 张彦军, 徐新喜, 张东辉. 双塔微型变压吸附制氧机实验和模拟[J]. 化工学报, 2019, 70(3): 969-978.
[4]	汪亚燕, 田彩霞, 丁兆阳, 唐忠利, 张东辉. 基于双回流变压吸附工艺的空气分离模拟及分析[J]. 化工学报, 2019, 70(10): 4002-4011.
[5]	王浩宇1，2，刘应书1，杨雄1. 空气预纯化的技术研究进展[J]. 化工进展, 2014, 33(03): 542-549.
[6]	常亮, 刘兴高. 内部热耦合空分塔的节能优化分析[J]. 化工学报, 2012, 63(9): 2936-2940.
[7]	李军1，2，王汝军1，张月明1，马占华1，孙兰义1. 隔壁精馏塔应用于空气分离的研究[J]. CIESC Journal, 2011, 30(11): 2393-.
[8]	张培昆, 王立, 李正强. 基于改进型吸附等温线模型的空气预纯化TSA双层床模拟 [J]. 化工学报, 2010, 61(10): 2604-2614.
[9]	脱瀚斐, 厉彦忠, 谭宏博. 空气分离与天然气液化组合循环 [J]. 化工学报, 2008, 59(10): 2498-2504.
[10]	卢洪,李成岳. 变压吸附空分制氮过程的研究 [J]. CIESC Journal, 2000, 51(5): 586-591.