-
全球能源系统正在往脱碳化方向发展并向净零情景转型[1],各个国家按照《巴黎协定》要求持续深度脱碳以应对气候变化,但自2016年各国签定《巴黎协定》以来,除2020年由于新冠疫情封控和防疫措施的影响碳排放量下降外,其余年份碳排放量持续增加。根据2022年第71版《BP世界能源年鉴》[2],2021年能源产生的CO2排放量比上年增长5.9%,达到338.8亿t;2021年因能源、过程排放、甲烷和放空燃烧产生的CO2当量排放比上年增长5.7%,达到389.8亿tCO2当量,基本接近疫情前的水平。2023年3月国际能源署发布的《2022年CO2排放》[3],2022年全球能源相关CO2排放量超过368亿t,比2021年增加3.21亿t,增长0.9%。中国2022年的能源相关CO2排放量达到114.8亿t,占世界碳排放的31%左右,与2021年相比基本持平,下降了0.2%,即2 300万t。同时,这也是自我国推动降碳减排以来的首次年度CO2排放下降[4]。可见,减少CO2的排放对于我国和世界各国仍是巨大的挑战。
碳捕集、利用和封存(Carbon Capture, Utilization and Storage, CCUS)技术是一种实现快速和深度脱碳所需的低碳能源与技术,其中利用CO2加氢催化合成甲醇[5-7]具有可行性,通过电解水得到的氢气,与捕集的CO2合成甲醇,不仅可实现碳资源的大规模有效利用,还可以将产物甲醇作为氢气的储存载体[8-9]。
甲醇可以作为替代燃料用于交通运输、发电和工业生产等领域,其能量密度高、易于储存和运输[10],可与传统石油燃料混合使用或直接替代燃料,具有较好的适应性,甲醇是一种理想的替代能源,能够大规模高效储能,也是一种重要的化工原料,用于生产甲醛和乙酸等有机化工产品[11-12]。常压下甲醇的沸点和熔点分别是64.8 ℃、−97.8 ℃,便于以液体形式储存、运输及配送,可充分利用现有石油、柴油等液体能源基础设施。2021年,我国的甲醇行业CO2排放量在化工全产业碳排放总量中占比超过17%,达1.96亿t[13],因此甲醇生产的绿色发展已成为必然趋势。
利用新能源电解水生产的绿氢,与全球排放的CO2催化合成甲醇和水,不仅可以实现CO2的减排利用,还可将可再生能源储存于液体燃料,用绿氢将CO2转化为甲醇对于解决CO2排放和能源紧缺问题极具经济价值和工业化应用前景。
目前,国内外有关CO2加氢催化合成甲醇的研究在催化剂[14-16]方面报道最多,主要包括Cu基催化剂、贵金属以及过渡元素催化剂、非贵金属氧化物催化剂3大类。另有一部分研究进行了相关工艺流程模拟,但未深入探究反应合成参数的选取,且更多仅关注甲醇合成过程,未与电解水结合。文章聚焦甲醇合成工艺模拟,考虑氢源来自电解水制氢,详细分析了温度、压力、CO2/H2进料比对甲醇产量的影响,为CO2催化加氢制甲醇工艺、电解水与甲醇合成结合提供了理论依据。
-
本文利用Aspen Plus模拟软件,基于Kiss[9]和Nyari[11]等人的模型构建了甲醇合成工艺,工艺流程如图1所示。
考虑H2气源直接来自于电解水制氢[17-19]得到的稳定5 MPa H2,简化考虑CO2气源同为5 MPa,两组气源和循环气于混合器内充分混合,然后加热至250 ℃。根据甲醇合成反应方程式的化学计量比,H2气源和CO2气源的进料量分别为7 500 kg/h和54 975 kg/h。
加热后的混合气体进入合成反应器进行甲醇合成反应,主要发生公式(1)~公式(3)的可逆反应[20],也就是H2和CO2可以直接反应生成甲醇,还可以发生逆水煤气变换反应,再由H2和CO2反应生成甲醇。
$$ \mathrm{CO+2H}_{\mathrm{2}}\leftrightharpoons\mathrm{CH}_{\mathrm{3}}\mathrm{OH}\;\Delta\mathit{H}_{\mathrm{298K}}\mathrm{=-90.8\;kJ/mol} $$ (1) $$ \mathrm{CO}_{\mathrm{2}}\mathrm{+H}_{\mathrm{2}}\leftrightharpoons\mathrm{CO+H}_{\mathrm{2}}\mathrm{O}\;\Delta\mathit{H}_{\mathrm{298K}}\mathrm{=+41.2\;kJ/mol} $$ (2) $$ \mathrm{CO}_{\mathrm{2}}\mathrm{+3H}_{\mathrm{2}}\leftrightharpoons\mathrm{CH}_{\mathrm{3}}\mathrm{OH+H}_{\mathrm{2}}\mathrm{O}\;\Delta\mathit{H}_{\mathrm{298K}}\mathrm{=-49.1\;kJ/mol} $$ (3) 以上公式(1)~公式(3)的可逆反应分别遵循速率方程公式(4)~公式(6)。
$$ \begin{split} &{r}_{A}={r}_{\mathrm{C}\mathrm{H}_3\mathrm{O}\cdot \mathrm{C}\mathrm{O}}=\\&{k}_{{A}}\frac{{K}_{\mathrm{C}\mathrm{O}}\left[{f}_{\mathrm{C}\mathrm{O}}\right.{f}_{{\mathrm{H}}_2}^{\tfrac{3}{2}}-{f}_{\mathrm{C}\mathrm{H}_3\mathrm{O}\mathrm{H}}/\left({K}_{A}\right.\sqrt{{f}_{\mathrm{H}_2}}\left.)\right]}{\left(1\right.+{K}_{\mathrm{C}\mathrm{O}}{f}_{\mathrm{C}\mathrm{O}}+{K}_{\mathrm{C}\mathrm{O}_2}{f}_{\mathrm{C}\mathrm{O}_2}\left)\right[\sqrt{{f}_{\mathrm{H}_2}}+\left({{K}_{\mathrm{H}_2\mathrm{O}}/\sqrt{{K}_{\mathrm{H}}})f}_{\mathrm{H}_2{\mathrm{O}}}\right]} \end{split} $$ (4) $$ \begin{split} &{r}_{B}={r}_{\mathrm{C}\mathrm{O}}={r}_{\mathrm{H}_2\mathrm{O}}=\\&{k}_{B}\frac{{K}_{\mathrm{C}\mathrm{O}2}\left[{f}_{\mathrm{C}\mathrm{O}_2}\right.{f}_{\mathrm{H}_2}-{f}_{\mathrm{H}_2\mathrm{O}}{f}_{\mathrm{C}\mathrm{O}}/{K}_{B}]}{\left(1\right.+{K}_{\mathrm{C}\mathrm{O}}{f}_{\mathrm{C}\mathrm{O}}+{K}_{\mathrm{C}\mathrm{O}_2}{f}_{\mathrm{C}\mathrm{O}_2}\left)\right[\sqrt{{f}_{\mathrm{H}_2}}+\left({{K}_{\mathrm{H}_2\mathrm{O}}/\sqrt{{K}_{\mathrm{H}}})f}_{\mathrm{H}_2{\mathrm{O}}}\right]} \end{split} $$ (5) $$ \begin{split} &{r}_{C}={r}_{\mathrm{C}\mathrm{H}_3\mathrm{O}\cdot \mathrm{C}\mathrm{O}_2}=\\&{k}_{C}\frac{{K}_{\mathrm{C}\mathrm{O}_2}\left[{f}_{\mathrm{C}\mathrm{O}_2}\right.{f}_{{\mathrm{H}}_2}^{\tfrac{3}{2}}-{{f}_{\mathrm{H}_2\mathrm{O}}f}_{\mathrm{C}\mathrm{H}_3\mathrm{O}\mathrm{H}}/({f}_{{\mathrm{H}}_2}^{\tfrac{3}{2}}{K}_{C}\left.)\right]}{\left(1\right.+{K}_{\mathrm{C}\mathrm{O}}{f}_{\mathrm{C}\mathrm{O}}+{K}_{\mathrm{C}\mathrm{O}_2}{f}_{\mathrm{C}\mathrm{O}_2}\left)\right[\sqrt{{f}_{\mathrm{H}_2}}+\left({{K}_{\mathrm{H}_2\mathrm{O}}/\sqrt{{K}_{\mathrm{H}}})f}_{\mathrm{H}_2{\mathrm{O}}}\right]} \end{split} $$ (6) 模拟中合成反应器选用了等温管式反应器——RPlug型反应器,该反应器遵从Langmuir-Hinshelwood-Hougen-Watson(LHHW)动力学模型。LHHW模型的内置表达式如公式(7)所示。
$$ r=\frac{k_{\mathrm{f}}\cdot D}{A} $$ (7) 式中:
$ {k}_\mathrm{f} $ ——动力学因子;
$ {D} $ ——驱动力表达式;
$ {A} $ ——吸附项。
动力学因子可以通过阿累尼乌斯公式(8)来表示:
$$ k_{\mathrm{f}}=K{{\mathrm{e}}}^{-{E}_{{\mathrm{a}}}/{\mathrm{R}}T} $$ (8) 式中:
$ K $ ——前因子;
$ {E}_{{\mathrm{a}}} $ ——活化能;
$ T $ ——反应温度;
$ {\mathrm{R}} $ ——摩尔气体常数,为8.31 J/(mol·K)。
甲醇合成反应的动力学因子关键参数如表1所示。
反应 K Ea/[J·(kmol)−1] 1 4.063 8×10−6 1.169 5×107 2 9.042 1×108 1.128 6×108 3 1.518 8×10−33 2.660 1×108 Table 1. Parameter choice of kinetic factors
驱动力表达式为对应速率方程的分子,可以写成公式(9)~公式(11):
$$ 反应A: {K}_{\mathrm{C}\mathrm{O}}{f}_{\mathrm{C}\mathrm{O}}{f}_{{\mathrm{H}}_2}^{\tfrac{3}{2}}-\frac{{K}_{\mathrm{C}\mathrm{O}}}{{K}_{A}}{f}_{\mathrm{C}\mathrm{H}_3\mathrm{O}\mathrm{H}}{f}_{{\mathrm{H}}_2}^{-\tfrac{1}{2}}\left[{Pa}^{3/2}\right] $$ (9) $$ 反应B: {K}_{\mathrm{C}\mathrm{O}_2}{f}_{\mathrm{C}\mathrm{O}_2}{f}_{\mathrm{H}_2}-\frac{{K}_{\mathrm{C}\mathrm{O}_2}}{{K}_{B}}{f}_{\mathrm{H}_2\mathrm{O}}{f}_{\mathrm{C}\mathrm{O}}\left[Pa\right] $$ (10) $$ 反应C: {K}_{\mathrm{C}\mathrm{O}_2}{f}_{\mathrm{C}\mathrm{O}_2}{f}_{{\mathrm{H}}_2}^{\frac{3}{2}}-{\frac{{K}_{C\mathrm{O}_2}}{{K}_{C}}{f}_{\mathrm{H}_2\mathrm{O}}}f_{\mathrm{C}\mathrm{H}_3\mathrm{O}\mathrm{H}}\left[{Pa}^{3/2}\right] $$ (11) 驱动力表达式中的驱动力常数项$ K $,可以由对数形式的公式(12)来表示。
$$ \mathrm{ln}\left(K\right)=A + \frac{B}{T} $$ (12) 同一反应的正反应和逆反应的驱动力常数项$ K $中A、B参数的取值如表2所示。
正反应 逆反应 A B A B −23.20 14 225 28.895 2 385 −22.48 9 777 −28.120 15 062 −22.48 9 777 23.974 3 222 Table 2. Parameter choice of driving force constant terms and coefficients
吸附项是速率方程公式(4)~ 公式(6)的分母,3个反应的吸附项是相同的,均可表示为公式(13)。
$$ \left(1+{K}_{\mathrm{C}\mathrm{O}}{f}_{\mathrm{C}\mathrm{O}}+{K}_{\mathrm{C}\mathrm{O}_2}{f}_{\mathrm{C}\mathrm{O}_2}\right) \left[\sqrt{{f}_{\mathrm{H}_2}}+\left({K}_{\mathrm{H}_2\mathrm{O}}/\sqrt{{K}_{\mathrm{H}}}\right)f_{\mathrm{H}_2\mathrm{O}}\right] $$ (13) 另外,根据进料量确定反应器管束为2 600个,催化剂中重量为2 776.5 kg。
由于甲醇合成反应单程转化率低,反应器出口产物中除甲醇外还有大量未反应气体和其他气体产物。因此,反应器出口产物首先经过冷凝器冷却至30 ºC,进入高压闪蒸罐内分离成气相和含甲醇液相,气相产物除少量作为驰放气排出,大部分压缩至5 MPa后成为循环气,液相进入低压闪蒸罐内再次进行气液分离,进一步除去未反应气体和其他气体产物。低压闪蒸罐出口气相与高压闪蒸罐的气相合为一股,压缩至5 MPa后成为循环气,进入合成反应器前的混合器与H2气源、CO2气源合成一股,重新进入反应器反应。低压闪蒸罐出口液相即为粗甲醇产物,粗甲醇组成简单,基本为甲醇和水。
-
甲醇合成工艺关键物流模拟结果如表3所示,表中数字1~11与图1中数字1~11位置对应。驰放气按照鲁奇工艺选取,驰放比例为0.75%。在H2气源和CO2气源的进料量分别为7 500 kg/h和54 975 kg/h时,CH3OH产量(数字10对应位置)为38 226.94 kg/h,质量分数为62.85%,碳的总转化率达到了89.5%。考虑年运行8 000 h,甲醇生产规模达30万t/a。
参数 1 2 3 4 5 6 7 8 9 10 11 温度/℃ 25.00 25.00 250.00 30.00 30.00 30.00 30.00 42.32 30.00 30.00 30.00 压力/MPa 5.00 5.00 5.00 5.00 4.50 4.50 4.50 5.00 1.00 1.00 4.50 摩尔流量/(kmol·h−1) 3 720.46 1 249.15 27 135.70 24 748.65 22 332.72 2 415.93 22 165.23 22 166.09 1.80 2 414.13 167.50 H2 3 720.46 0.00 21 467.83 17 882.76 17 882.19 0.58 17 748.07 17 747.37 0.53 0.05 134.12 CO 0.00 0.00 599.17 603.66 603.63 0.04 599.10 599.17 0.03 0.01 4.53 CO2 0.00 1 249.15 4 989.52 3 791.50 3 767.15 24.35 3 738.90 3 740.36 1.22 23.13 28.25 H2O 0.00 0.00 12.31 1 210.33 12.40 1 197.93 12.31 12.31 0.00 1 197.92 0.09 CH3OH 0.00 0.00 66.87 1 260.39 67.35 1 193.04 66.85 66.87 0.02 1 193.02 0.51 摩尔分数H2 1.000 0 0.000 0 0.791 1 0.722 8 0.800 7 0.000 2 0.800 7 0.800 7 0.293 4 0.000 0 0.800 7 摩尔分数CO 0.000 0 0.000 0 0.022 1 0.024 4 0.027 0 0.000 0 0.027 0 0.027 0 0.017 0 0.000 0 0.027 0 摩尔分数CO2 0.000 0 1.000 0 0.183 9 0.153 2 0.168 7 0.010 1 0.168 7 0.168 7 0.674 9 0.009 6 0.168 7 摩尔分数H2O 0.000 0 0.000 0 0.000 5 0.048 9 0.000 6 0.495 9 0.000 6 0.000 6 0.002 3 0.496 2 0.000 6 摩尔分数CH3OH 0.000 0 0.000 0 0.002 5 0.050 9 0.003 0 0.493 8 0.003 0 0.003 0 0.012 4 0.494 2 0.003 0 质量流量/(kg·h−1) 7 500.00 54 975.00 282 011.80 282 011.80 221 129.40 60 882.34 219 470.90 219 536.75 56.27 60 826.08 1 658.47 H2 7 500.00 0.00 43 276.57 36 049.51 36 048.34 1.16 35 777.98 35 776.57 1.07 0.10 270.36 CO 0.00 0.00 16 783.07 16 908.86 16 907.86 1.01 16 781.05 16 783.07 0.86 0.15 126.81 CO2 0.00 54 975.00 219 587.60 166 863.20 165 791.70 1 071.53 164 548.30 164 612.64 53.55 1 017.98 1 243.44 H2O 0.00 0.00 221.85 21 804.43 223.45 21 580.98 221.77 221.85 0.08 21 580.91 1.68 CH3OH 0.00 0.00 2 142.62 40 385.73 2 158.07 38 227.66 2 141.89 2 142.62 0.72 38 226.94 16.19 质量分数H2 1.000 0 0.000 0 0.153 5 0.127 8 0.163 0 0.000 0 0.163 0 0.163 0 0.019 0 0.000 0 0.163 0 质量分数CO 0.000 0 0.000 0 0.059 5 0.060 0 0.076 5 0.000 0 0.076 5 0.076 4 0.015 3 0.000 0 0.076 5 质量分数CO2 0.000 0 1.000 0 0.778 6 0.591 7 0.749 8 0.017 6 0.749 8 0.749 8 0.951 7 0.016 7 0.749 8 质量分数H2O 0.000 0 0.000 0 0.000 8 0.077 3 0.001 0 0.354 5 0.001 0 0.001 0 0.001 3 0.354 8 0.001 0 质量分数CH3OH 0.000 0 0.000 0 0.007 6 0.143 2 0.009 8 0.627 9 0.009 8 0.009 8 0.012 7 0.628 5 0.009 8 Table 3. Key stream simulation results of methanol synthesis process
-
甲醇产量在反应器条件为5 MPa、220~300 ℃温度变化范围的模拟结果如图2所示。可以看出,在反应器温度为220 ℃时,甲醇质量分数约62.73%,甲醇产量为37 394 kg/h;随着温度升高,甲醇的纯度和产量升高,温度升高到250 ℃左右,甲醇的纯度和产量到达最高,质量分数为62.85%,产量为38 226 kg/h;温度再升高后,甲醇的纯度和产量均开始降低,在反应器温度为300 ℃时,甲醇质量分数约62.43%,甲醇产量为36 776 kg/h。这主要是由于甲醇合成反应,从热力学角度来说,正反应是一个放热反应,温度降低利于正反应进行;但从动力学角度来说,温度升高,分子运动加快,反应速率增加。因此,存在一个最佳反应温度,在250 ℃左右,也是目前甲醇合成工艺中通常选取的温度。
甲醇产量在反应器条件为250 ℃、 4~6.5 MPa压力变化范围的模拟结果如图3所示。由于甲醇合成反应为分子数减小的反应,增加压力会使反应向正向进行。从图3也可以看出,随着压力升高,反应器出口的甲醇的质量分数和产量也随之升高,但增长速率逐渐减缓,比如4.5~5 MPa,甲醇产量增加133.7 kg/h;5~5.5 MPa,甲醇产量增加113.2 kg/h;5.5~6 MPa,甲醇产量增加93.1 kg/h。而压力升高,压缩机功耗随之势必会增加,另外考虑甲醇合成与电解水制氢直接结合,5 MPa的电解水制氢产物压力是目前国产质子膜电解水制氢工艺里可达到的最大压力。因此,甲醇合成器选取5 MPa为较合理的参数选择。
-
固定H2进料量7 500 kg/h,在CO2进料量变化,进料比为6~8.7时,甲醇产量的变化结果如表4所示。
H2/(kg·h−1) CO2/(kg·h−1) 进料比 甲醇产量/(kg·h−1) 7 500 45 000 6.0 32 525.78 7 500 47 250 6.3 34 130.90 7 500 49 500 6.6 35 686.52 7 500 51 750 6.9 37 078.29 7 500 54 000 7.2 38 025.43 7 500 54 975 7.33 38 226.94 7 500 56 250 7.5 38 399.76 7 500 58 500 7.8 38 453.77 7 500 60 750 8.1 38 364.57 7 500 63 000 8.4 38 236.43 7 500 65 250 8.7 38 094.35 Table 4. Methanol yield under different H2/CO2 feed ratios
可以看出,甲醇产量随进料比增加,呈现先增后减的趋势,但均可保持较高的甲醇产量。进料比在7.8左右,甲醇产量达到最高,约38 453.8 kg/h,进料比再增加后,产量随之降低。根据CO2与H2合成甲醇的化学反应式计算,CO2与H2的进料比为7.33左右,反应正好完成,与模拟得到的最佳进料比7.8接近。因此,甲醇反应器在启停或波动期间,进料比可略有波动,当反应器稳定运行时,进料比控制在7.8左右,甲醇产量最优。
-
随着“双碳”政策的不断推进和深化,可再生能源电力成本降低,氢能应用规模化和市场逐渐成熟,通过电解水技术制取绿氢将是未来低碳经济的主流方向,成为H2的主要来源。电解水制氢是通过电能使电解槽内水分子在电极上发生电化学反应,目前商业化的水电解制氢主要是碱性电解水制氢和质子交换膜制氢技术。对于碱性制氢装置来说,经过电解槽和氢氧分离器后,H2纯度可达到99.8%;对于质子膜制氢装置来说,经过电解槽和氢氧分离器后,H2纯度可达到99.9%。2种制氢装置均要经过H2干燥装置,H2干燥装置主要是除去H2中的水分,将H2纯度处理到99.99%以上。
电解水制氢通常采用变温吸附工艺干燥提纯,分子筛对水分的吸附达到平衡后,再利用H2对分子筛进行再生,再生温度约170 ℃。对于7 500 kg/h的H2而言,99.8%~99.9%的H2相当于含水量75~150 kg/h。
根据表5不同含水率对应的甲醇纯度和产量表可以看出,H2的含水量对甲醇的产量有一定的影响,每增加100 kg/h的含水量,甲醇质量分数降低约0.1%,产量变化仅为几kg/h数量级,基本可忽略不计。因此,电解水氢气产物的含水量对甲醇的纯度和产量影响很小,而取消H2干燥装置,可以减少约8 000万元的设备费用,同时还可减少再生用的氢气损耗,使得制氢和新能源的跟随性更好,启动速度加快。那么,电解水得到的H2不经过干燥处理,直接与甲醇合成工艺结合是可行的,且对甲醇产量基本无影响。
H2/
(kg·h−1)CO2/
(kg·h−1)H2O/
(kg·h−1)甲醇质量分数 甲醇产量/
(kg·h−1)7 500 54 975 600 0.622 402 38 212.37 7 500 54 975 500 0.623 413 38 213.63 7 500 54 975 400 0.624 419 38 216.62 7 500 54 975 300 0.625 419 38 218.61 7 500 54 975 200 0.626 404 38 217.91 7 500 54 975 100 0.627 411 38 221.97 Table 5. Methanol purity and yield under different water contents
-
本文通过流程模拟对CO2加氢催化制甲醇工艺系统进行了分析,确定了关键工艺参数的选取依据,探究了甲醇合成与电解水结合的可行性。在合成反应器温度、压力分别为250 ℃、5 MPa的条件下,模拟流程H2和CO2的进料量分别为7 500 kg/h和54 975 kg/h,CO2/H2进料比为7.33,CH3OH产量可达38 226.94 kg/h,碳的总转化率为89.5%。目前甲醇合成工艺中通常选取的反应器温度为250 ℃,是最佳合成温度。考虑与电解水工艺的结合,5 MPa的反应器压力是合理的压力选择值。甲醇合成系统稳定运行时,CO2/H2进料比控制在7.8左右,甲醇产量最优,与化学计量数的进料比相差不大。电解水得到的H2不经过干燥处理,直接与甲醇合成工艺结合是可行的,可简化工艺,且对甲醇产量基本无影响。
根据全球碳中和战略的发展要求,新能源电力电解水制氢结合碳捕集技术生产甲醇是重要发展趋势。该技术目前面临的主要挑战在于波动性的新能源电力导致的氢气原料供应波动性如何匹配下游稳定的合成甲醇工艺。此外,高压电解水制氢、甲醇合成催化剂等也是新能源制氢耦合CO2加氢合成甲醇工艺研究的重要方向。随着风电、光电等可再生能源发电成本的下降及甲醇合成适应波动性氢气原料波动技术的发展,新能源制氢耦合CO2加氢合成甲醇工艺是具有应用前景的甲醇合成技术。
System Analysis of Methanol Synthesis by Carbon Dioxide Hydrogenation
doi: 10.16516/j.ceec.2024.5.05
- Received Date: 2023-08-14
- Rev Recd Date: 2023-09-25
- Available Online: 2024-09-30
- Publish Date: 2024-09-10
-
Key words:
- methanol /
- methanol synthesis /
- carbon dioxide /
- CO2 to methanol /
- hydrogen
Abstract:
Citation: | ZHAI Yunchu, JIN Liyan, XIE Jiajia, et al. System analysis of methanol synthesis by carbon dioxide hydrogenation [J]. Southern energy construction, 2024, 11(5): 50-56 doi: 10.16516/j.ceec.2024.5.05 |