一、 前言
2020年，我国正式提出了碳达峰、碳中和（“双碳”）战略目标，以积极应对气候变化并着力化解资源环境约束。京津冀地区作为我国三大城市群之一，是能源消耗密集区 [1]，碳排放量约占全国的11% [2]，碳排放强度高出全国平均水平约40%。在推进“双碳”目标的进程中，京津冀地区“分量”较重。然而，京津冀三省市所处发展阶段有差异，面临的低碳转型挑战也有区别 [3]：北京市已完成工业化，碳达峰目标容易实现 [4]，但受限于风光资源短缺、绿电供应不足等约束，能源、建筑、交通等部门的电气化转型面临挑战，碳减排进展相对缓慢 [5]；天津市以制造业为支柱（第二产业占比为34.1%且保持较快增长 [6]），能耗和碳排放强度较高，产业低碳转型以及结构性调整面临挑战；河北省仍处于工业化阶段（第一、二产业占比分别为10.1%、37.6% [7]），钢铁等重工业占比较高，在淘汰落后产能、推动产业转型升级、提升能效方面压力较大。
作为国家战略，京津冀协同发展以产业协同、交通一体化等为重点领域，构成推进区域内“双碳”工作的有力依托：产业协同有利于三省市发挥各自比较优势，共同形成京津冀地区的产业链和创新链，协同开展产业转型升级和能效优化；能源协同可同步解决河北省绿电消纳、北京市绿电短缺等问题；交通一体化在优化运输结构、提高交通效率、助力交通部门减排方面潜力突出 [8]。在“双碳”背景下，全面把握京津冀协同发展的新内涵，辨识并提出有利于推动“双碳”工作的协同措施，对京津冀地区的更高质量协同发展具有重要意义。
目前，学术层面的京津冀协同发展以定性研究居多，而定量研究不显充分；现有的京津冀地区协同减碳路径定量研究主要包括三类。① 针对京津冀地区的单一或多个城市进行建模，如北京市 [9,10]、张家口市 [11]、京津冀及周边“2+26”城市 [12]；多针对道路交通、钢铁等单一部门 / 行业，所采用的模型结构精细度不高、适应性一般。② 将京津冀地区当作整体，开展单一部门的建模研究，如交通 [13,14]、钢铁 [15]、电力 [16]、能源 [17,18]等，而未考虑区域协同效应；将三省市的数据直接加总并整体性构建模型，未考虑三省市在经济结构、能源结构、发展政策等方面的不同，没有反映京津冀协同中的差异性定位。③ 定量分析单一协同措施对京津冀地区的影响作用，如产业转移 [19]、能源协同 [20,21]、污染治理协同 [22]。
梳理现有文献可见，有关京津冀协同减碳研究存在一些不足之处：① 政策背景考虑不充分，如京津冀“双碳”路径研究的情景设定没有及时全面地引入“十四五”规划内容，导致情景设置与最新规划脱节（“十四五”时期是低碳转型、碳达峰的窗口期 [23]，对京津冀“双碳”发展路径有着重要影响）；② 未能整体性考虑协同措施，有关定量研究基本没有将交通一体化、产业转移与协同升级、能源协同发展等跨领域协同措施纳入统一分析框架；③ 模型结构不够精细，应用于钢铁等重要部门时难以充分体现行业能效提高、产量控制等对减碳工作的影响。
针对于此，本文尝试在以下方面开展探索：尽量全面考虑京津冀协同发展、“十四五”规划等政策，基于情景分析法评估现有政策对京津冀地区的碳减排效应；自下而上地构建考虑京津冀地区特征的长期能源替代规划系统（LEAP-BTH）模型，将交通一体化、产业协同升级、能源协同发展等跨领域协同措施纳入统一框架，力求系统分析京津冀协同发展的碳减排潜力；合理提高模型精细度，体现高能耗部门的能效提高、产量控制等对碳减排的影响，改善模型结果的准确性。
二、 研究方法和基础数据
（一） 基于LEAP模型的能源需求与碳排放核算方法
1. 模型框架
LEAP模型 [24]支持基于情景分析的能源需求预测，加总各部门的能源需求，比较并分析本地资源产出能否满足需求；若不满足需求则计算由此引起的资源进出口量，实现资源、能源需求量、产量、进出口的平衡。在LEAP模型中，根据各经济部门的历史数据和政策措施，输入各个情景的具体设置，据此计算给定情景下未来年的能源需求，进而计算能源结构和温室气体排放量；具有较高的建模灵活性、较强的情景分析能力，在区域经济、能源、碳排放的情景分析研究方面获得较多应用 [9~15]。本文基于LEAP模型框架，结合京津冀地区的能源结构、产业结构等特征，构建了LEAP-BTH模型（见图1）；包含了关键假设、终端能源需求、能源加工转换等主要模块，将各部门的活动简化为能源、物质的消耗与产出。LEAP-BTH模型以2020年为基准年，预测期为2021—2060年，时间步长为1年，区域范围是北京市、天津市、河北省。
图1 LEAP–BTH模型结构
2. 关键假设
常住人口及城镇化率。2020年北京市、天津市的城镇化率约为85%，处于较高水平；北京市人口趋于稳定，天津市、河北省近5年平均变化幅度均小于0.3%，没有出现明显的增 / 减趋势。因此，假设京津冀地区在预测期内的常住人口保持不变。受数据可得性的限制，将北京市、天津市的城市和农村居民生活部门的能源消耗合并计算。2020年，河北省城镇化率为60.1%，据历史趋势预测2030年城镇化率上升至75%，之后保持不变。
地区生产总值及三次产业结构。根据京津冀协同的整体布局、三省市各自的功能定位及其“十四五”规划，参照实际发展趋势及相关文献预测结果，结合对相关管理部门的访谈内容，开展京津冀地区年均经济增速、三次产业结构等参数的评估和假定（见表1、表2）。例如，预测天津市第二产业占比提高的主要原因是政策层面高度重视工业发展，近5年的占比有逐渐提高趋势 [6]；预测北京市经济增速的主要依据是有关文献测算的中长期经济增长率（0.2% [25]至4% [10]）以及当前发展政策、管理部门调研内容等；其余参数设置的考虑类似，为凝练表达而不赘述。研究期内北京市、天津市、河北省的第三产业占比分别于2035年、2025年、2035年之后保持不变。
表1 京津冀地区年均经济增速预测
表2 京津冀地区三次产业结构现状及预测
3. 能源加工转化
能源加工转化模块包括天然气、煤炭、石油等一次能源开采，石油加工、焦化等一次能源加工，电力、热力生产与供应等过程；主要应用LEAP模型的缺省峰值负荷曲线、各发电类型的最大可利用率。
4. 能源需求
能源需求模块涵盖居民生活、交通运输、工业、服务业、农业、建筑业等主要用能部门。工业部门下各产业的能耗计算方法有两种：针对钢铁、水泥、平板玻璃、造纸等高能耗部门，采用产品产量代表活动水平、单位产品能耗代表能源强度；对于其他工业部门，采用增加值代表活动水平、单位增加值能耗代表能源强度。
交通运输分为客运和货运，其中客运包括私人汽车、公交车、出租车、轨道交通，活动水平以行驶里程衡量。2020年北京市交通部门中各类交通方式的活动水平与能源结构如表3所示，天津市、河北省的数据结构与北京市相似，而行驶里程、能耗类型占比等参数依据地区特征而设置。
表3 北京市交通部门活动水平与能源结构（2020年）
5. 碳排放
碳排放的核算范围是与能源生产、转换、消耗相关的各种经济活动及其产生的直接碳排放，不考虑电力等二次能源调入产生的间接碳排放以及农业、畜牧业、林业、土地利用等方面的碳排放变化，甲烷等非CO 2温室气体排放等。采用联合国政府间气候变化专门委员会（IPCC，2006）方法计算碳排放量。
（二） 基础数据与情景设置
表4 研究采用的核心指标及数据来源
表5 3类情景的主要内涵及子情景
1. 基准情景
2021—2060年京津冀三省市各部门的活动水平、能源强度等参数，主要依据国家和区域“十四五”规划、现行的区域协同发展战略等政策文件及相关文献等进行设置。
对于交通部门，客运交通出行结构不断优化，绿色出行占比持续增加；轨道交通出行量年增长率在“十四五”时期为10%，随后降为3%；主要参数有公共交通、轨道交通的出行占比，机动车保有量，新能源私人汽车、公交车、出租车、货车的保有量 [34]等。
对于工业部门，用能结构调整、生产技术提高对工业的影响等是考虑因素，主要参数有“十四五”时期能源强度、终端用能单位电气化率等；有关具体工业产品的产量预测，主要参照当前产量的变化趋势以及有关产能控制的政策。需要说明的是，针对主要工业部门，采用产品产量（如钢铁产量）表示活动水平，由单位产品能耗（如吨钢生产能耗）表示能源强度；课题组还调研访谈了相关管理部门，以更好把握相关数据的后续变化趋势。例如，在推行钢铁行业“去产能”政策后，河北省的钢铁产量在2021—2023年每年下降8%，之后每年下降1%；2025年之前吨钢生产能耗每年下降0.4%，之后每年下降1%。
对于居民生活、建筑业与服务业部门，推广装配式建筑等低碳绿色的建筑方式、对现有建筑进行节能改造等措施是考虑因素，主要参数有建筑业、服务业的能源强度，居民生活的能源消耗强度，居民终端用能电气化率等。“十四五”时期，建筑业、服务业的能源强度每年下降3%，之后每年下降1%；2035年居民终端用能电气化率达到90%，完全淘汰煤炭消耗。
对于能源部门，煤电、气电、光伏发电、风电的装机容量和发电量，供暖能源消耗总量、供暖方式、供暖能源种类 [35]等是考虑因素。
2. 低碳情景
挖掘现行政策的碳减排潜力，结合文献资料，预测各子情景下相关经济部门的变化情况。以北京市为例，低碳情景的4个子情景及其包含的关键措施如表6所示。
表6 北京市低碳情景的子情景及关键措施
低碳交通包括交通电气化、运输方式转变：对于前者，在“十四五”规划的基础上加快电动汽车推广并在2050年完全取代汽油车，公交车、出租车分别在2030年、2035年完成电气化替代 [10]；对于后者，提高公共交通对私人交通的替代水平，实现轨道交通客运量、公交车数量的稳定增长。
对于可再生能源替代，在基准情景的基础上，加快煤电退役进度，推广可再生能源发电和供暖模式 [36]；北京市挖掘分布式光伏和地热供暖的潜力，河北省发展风电和光伏发电，天津市提早煤电退役时间。
对于工业绿色升级，在基准情景的基础上，推广应用低碳节能技术，促进能源强度持续下降，提高用能终端的电气化替代水平；吨钢生产能耗在“十四五”时期每年下降0.6%，之后每年下降1.5%。
对于低碳建筑，加快用能终端的电气化替代， 2030年居民生活部门的煤炭被完全取代，电气化率为90%；提高建筑的能源效率，居民生活的能源强度在2035年之后不再增长，建筑业、服务业在2035年之后能源强度每年下降1.5%。
3. 协同情景
在低碳情景的基础上，强化京津冀协同发展的政策力度，进一步提出有利于区域低碳发展的协同措施，同时体现一定的政策前瞻性。协同情景的4个子情景及其包含的关键措施如表7所示。
表7 协同情景的子情景及关键措施
对于交通一体化，结合“京津冀交通一体化”远景规划，2035年基本实现“轨道上的京津冀”，轨道交通成为出行的主要方式，而私人汽车的年行驶里程稳步减少；发展氢能产业链，“十四五”时期推广1×10 4辆氢能源车，2035年氢能源车占私人汽车的比例为5%；共建共享全域路网的电动车配套基础设施，电动车占比进一步提高；铁路网基础设施建设稳步实施，加速货运部门“公转铁”并发展“港铁联运”，提高铁路货运占比及其电气化率。
对于产业协同发展，北京市转移出一般性的工业和制造业，发展智能和高端制造业，优化工业结构，发挥科技引领作用；天津市、河北省承接产业转移和对口产业帮扶，整合产业资源，推动产业链协同升级，进一步提高综合能效 [37]，2060年吨钢生产能耗下降53% [38]。
对于低碳技术共享，将新建城区（如北京城市副中心、河北雄安新区）的低碳建筑技术扩散到整个区域，推广低能耗建筑 [38]；北京市建筑业能源强度进一步降低，2060年天津市和河北省的建筑业、服务业、城市居民生活的能源强度下降至北京市低碳情景水平。
对于能源协同，京津冀地区统筹构建高效、智能电力系统，支撑更高比例的可再生能源上网，提高可再生能源消纳率，促进风光资源的充分开发利用 [39]；2050年煤电、气电机组基本退役并转为调峰与储备能源，电力完全由可再生能源提供。
三、 京津冀地区能源需求与碳排放模拟结果
（一） 基准情景的能源需求与碳排放预测
图2 基准情景下的京津冀地区能源需求与碳排放情况
1. 北京市的预测结果
在预测期内，北京市的能源需求量持续上升，未能达到峰值。2035年前，能源需求仍以较快速度增长，“十四五”时期的终端能源消费量从6.29×10 7 tce增长到7.12×10 7 tce，对应增速为13.2%，低于“十四五”规划中的增长率（19%）。到2035年，终端能源消费量继续增加到8.97×10 7 tce。由于未考虑航空和管道运输，基准年的能源需求估算略小于实际值，因而2035年北京市的实际能源需求会超过“十四五”规划提出的9.0×10 7 tce控制目标，也不能实现“十四五”规划要求的2035年能源需求达峰目标。在2035年之后，能源需求的增长速度放缓，但增长量依然较大。到2060年，能源需求为1.085×10 8 tce，较基准年增长72.4%。
根据预测，北京市2060年碳排放量下降为2.255×10 7 t，较2020年下降59.3%，不能实现碳中和目标（相关研究表明 [42,43]，2060年我国将通过碳捕集利用与封存、碳移除等技术捕集约2.3×10 9 t CO 2，约占2020年CO 2排放量的20%；本研究判断区域能否实现碳中和目标，标准定为2060年该区域的碳排放量下降为2020年的20%以下）。受新型冠状病毒肺炎（COVID-19）疫情对行业开工率、交通活动水平的影响，基准年的碳排放低于正常情况；2021—2023年的碳排放量比基准年略有上升（约为7.34×10 8 t CO 2），合理推测北京市可以实现“十四五”规划要求的碳排放“稳中有降”目标。2025年之后，碳排放的下降速度有所加快，能效提高、结构转型带来的减排效果逐步显现。2030年、2035年碳排放量分别下降为6.98×10 7 t CO 2、6.29×10 7 t CO 2，较2020年降低4.8%、14.2%。尽管2060年碳排放降幅接近60%，但相较碳中和目标还有一定差距。
2. 天津市的预测结果
在预测期内，天津市的终端能源需求快速增长，未能达到峰值（2060年为2.834×10 8 tce，约为基准年的3.4倍），相应增长主要由工业部门驱动。“十四五”时期，天津市制定的经济增长目标为6%，主要依赖制造业等工业部门拉动经济增长，第二产业增加值的占比进一步提高，导致能源需求保持高位。
根据预测，天津市2060年的碳排放量下降为7.26×10 7 t CO 2，较2020年下降59.9%，不能实现碳中和目标。究其原因，煤电存在碳锁定效应，而天津市的煤电机组难以在2060年完全退役；天津市作为我国北方的制造业中心和航运中心，工业和服务业在较长时期内都将保持一定的增速，而基准情景下的能效提高不足以抵消经济增长对能源的增量需求。
3. 河北省的预测结果
河北省的终端能源需求在2020—2023年随着工业压减产能而有所下降，2024年之后将恢复增长，在预测期内未能达到峰值。河北省的能耗集中在工业部门，受“十四五”前期的产能控制政策影响，钢铁等高能耗工业产品的产能出现下降，导致短期内能耗降低；在长期视角下，随着政策压力减轻、经济增长恢复，相应的能源需求重新出现增长趋势。
（二） 碳排放总量方面的情景模拟分析
2060年，预计基准情景、低碳情景、协同情景下的京津冀地区碳排放量分别为5.87×10 8 t CO 2、4.88×10 8 t CO 2、2.32×10 8 t CO 2（见图3）。在基准情景、低碳情景下，京津冀地区的碳排放量都将保持下降趋势；北京市、河北省、天津市先后达峰，但都不能实现碳中和目标。仅在协同情景下，2060年京津冀地区的碳排放量比2020年下降80.4%，能够实现碳中和目标。
图3 3类情景下的京津冀地区碳排放量
对于北京市，在低碳情景、协同情景下，2060年的碳排放量分别为1.46×10 7 t CO 2、9.3×10 6 t CO 2，分别为基准年的19.9%、12.7%，即在低碳情景下2060年可以实现碳中和目标，在协同情景下可以提前实现碳中和目标。
对于天津市，在低碳情景下，2060年碳排放量为4.76×10 7 t CO 2，为基准年的73.7%；在协同情景下，2060年碳排放量为下2.76×10 7 t CO 2，比基准年下降84.7%，即协同情景下可以提前实现碳中和目标。
对于河北省，3类情景下2020年后的碳排放量均呈下降趋势。在基准情景下，碳排放量下降幅度不足以实现2060年碳中和目标；在低碳情景下，2060年碳排放量为4.259×10 8 t CO 2；在协同情景下，2060年碳排放量为1.954×10 8 t CO 2，为基准年的79%，能够实现碳中和目标。
四、 京津冀协同推进碳达峰碳中和路径分析
（一） 重点减排部门
获得了预测期内、协同情景下的京津冀地区各部门累计碳减排潜力占比（见图4），对比各经济部门的累计减排量可以发现，北京市的重点减排部门是交通运输，河北省、天津市的重点减排部门是工业和发电。具体而言，北京市碳减排的显著贡献部门是交通部门（占31.2%），其次是服务业（占19.4%）、供热（占19%）、发电（占16.1%）；天津市碳减排的显著贡献部门是工业（占31.5%）和发电（占30%），其次是供热（占13.5%）、服务业（占8%）；河北省碳减排的显著贡献部门是工业（占41.3%）和发电(占25.8%)，其次是供热（占9.6%）、交通（占6.8%）。
图4 协同情景下的各部门累计碳减排潜力占比
由此可见，京津冀三省市的重点减排领域有所不同。北京市的能源结构以天然气为主、依赖外调电，发电部门的碳减排潜力不大；而在终端排放中交通部门碳排放量的占比最高，加之服务业发达，因而减少碳排放依赖交通、服务业等部门。天津市、河北省的能源结构都是以煤炭为主且产业结构中的工业占比较高，因而重点减排部门都是工业和发电。此外，河北省的碳排放量相当于京津两市的数倍，因而京津冀地区的整体碳排放情况主要由河北省决定，故京津冀地区的重点减排部门主要是工业、发电、热力生产。
（二） 关键减排路径
分析各项措施累计减排的潜力占比（见图5），据此辨识京津冀三省市及整个地区的关键性减排路径。
图5 各项措施的累计碳减排潜力占比
北京市的关键减排措施是低碳交通、低碳建筑，减排潜力占比分别为28.6%、26.4%；随后是交通一体化、能源协同发展，减排潜力占比分别为18.8%、10.7%。具体而言，加快推广新能源汽车，力争2035年电动车占私家车的比例为60%，氢能源车占比为5%；推动轨道交通等绿色交通基础设施的共建共享，使轨道交通成为出行的主要方式。也要注意到，低碳情景下北京市各项措施的累计减排潜力占比为61.7%，较协同情景下减排措施的效果更明显，这是因为北京市的交通部门电气化、运输模式转变、低碳建筑推广等措施依赖当地政策的推动，而受京津冀协同的影响相对不大。
天津市的关键减排措施是可再生能源替代、工业绿色升级，减排潜力占比分别为33%、19.2%；随后是低碳技术共享、能源协同发展，减排潜力占比分别为14.1%、10.4%。具体而言，煤电和气电机组到2050年基本完成退役，转为调峰和储备能源，电力完全由可再生能源提供；推动工业绿色升级，更好维持天津市的制造业中心地位。低碳情景下天津市各项措施的累计减排潜力占比为59.1%，也较协同情景下减排措施的效果更明显。究其原因，天津市的发电结构以煤炭为主且能源需求部门耗电量大，故煤电退役、可再生能源替代对当地碳减排的影响更为直接，而受京津冀协同的影响相对不大。
河北省的关键减排措施是产业协同、能源协同，减排潜力占比分别为34.1%、22.2%；随后是工业绿色升级、可再生能源替代。河北省钢铁产量极高，超过40%的碳排放源于钢铁生产的能源消耗过程，因而压减钢铁产能、应用节能技术进行产业升级在碳减排工作中起着关键作用。河北省拥有较为丰富的风、光自然资源，通过京津冀能源协同发展，可以运用区域内的资金和技术优势，更好将自然资源优势转化为低碳发展优势；京津冀地区包括包括“风光荷储”一体化建设在内的能源基础设施协同，有利于促进河北省的可再生能源就近消纳，对河北省的碳减排工作也有积极作用。
对于整个京津冀地区而言，协同措施的碳减排潜力占比超过50%，又以产业协同、能源协为主要贡献类别。具体而言，河北省在承接北京市制造业转移的基础上，推动产业链协同创新，整合优质产业资源，进一步提高能效；2060年工业能源强度相比基准年下降75%，吨钢生产能耗下降53%。
五、 研究结论与发展建议
（一） 研究结论
按照现有政策实施后的发展趋势（即基准情景），北京市、河北省、天津市依次实现碳达峰，但三省市在2060年均不能实现碳中和目标。在低碳情景下，北京市可在2060年实现碳中和目标，天津市、河北省的碳排放量将进一步降低，但仍不能实现碳中和目标。对于经济体量大、增长速度较快、工业和城市密集的京津冀地区，碳中和对经济社会转型提出了更高的要求。
仅在协同情景下，京津冀三省市才能都在2060年（基本）实现碳中和目标。这就说明京津冀协同发展具有极大的碳减排潜力，是三省市完成碳中和目标的有效途径。
京津冀三省市的重点减排部门和关键路径存在明显差异：北京市重在低碳交通、低碳建筑、交通一体化、能源协同发展，天津市重在可再生能源替代、工业绿色升级、低碳技术共享、能源协同发展，河北省重在产业协同、能源协同。
着眼碳中和目标，京津冀协同发展对河北省的碳减排贡献将最为突出，也高于河北省加强自身减碳措施执行力度（低碳情景）的效果。河北省实现碳中和目标最为重要的产业协同和能源协同措施，将依赖北京市、天津市溢出的资金和技术优势。对于北京市，交通部门电气化、运输模式转变、低碳建筑推广等措施更多依赖于本地政策的推动，故京津冀协同的收益作用相对较小。对于天津市，可再生能源替代和工业绿色升级等关键减排措施对碳减排的影响更为直接，而受京津冀协同的影响相对不大。
（二） 发展建议
北京市增强交通、建筑部门的碳减排政策力度。对于交通部门，倾向公共交通、自行车和行人基础设施建设，将城市副中心的步行友好示范社区建设经验逐步推广到全市，改善公共交通与自行车、步行等的连接便利性；优化市郊轨道交通布局，使城区与周边县区及外部城市的连接更为高效便捷。对于建筑部门，参照国际先进的建筑标准，提高被动式建筑、绿色屋顶等低碳建筑技术的应用水平，进一步降低建筑的能耗；加快实施现存建筑的节能化改造，因地制宜应用分布式光伏资源。北京市加强与河北省的可再生能源合作开发，支持绿电入京，合理降低碳中和技术成本；运用科技优势，与河北省承德市、张家口市等风光资源丰富的地区协同开发风光资源，为外地绿电入京提供优先保障；在构建京津冀氢能产业链过程中发挥引领作用，发挥技术和市场优势，建设并完善氢能基础设施，推广应用氢能汽车。
天津市侧重煤电和工业部门碳减排、京津冀能源协同及交通一体化。严控新增煤电，尽快明确现役煤电的退役进度，发展分布式光伏等可再生能源，积极引入绿电。对于工业部门，以科技创新引领制造业低碳发展，尽量控制制造业带来的碳排放增量。对于在交通部门，深度参与京津冀交通一体化建设，发挥区域“海上门户”的交通枢纽作用；增强货运“公转铁”水平，提高交通电气化、智能化程度。
河北省把握京津冀能源协同、产业协同的发展机遇。实施京津冀地区的能源基础设施协同建设，发挥西部和北部的风光资源优势。与北京市、天津市合作，解决可再生能源开发、储存、消纳相关的资金与技术难题，增强可再生能源的外送能力，扩大可再生能源当地消纳范围。推动京津冀产业协同发展，发挥北京市、天津市技术转化基地作用，促进工业部门的能源结构转型与低碳发展，降低化石能源占比并提高能源效率。
利益冲突声明
本文作者在此声明彼此之间不存在任何利益冲突或财务冲突。
Received date: November 21, 2022; Revised date: January 5, 2023
Corresponding author: Lu Xi is a professor from the School of Environment of Tsinghua University. His major research fields include carbon neutral energy systems, greenhouse gas and air pollutantemissions. E-mail: xilu@tsinghua.edu.cn
Funding project: Chinese Academy of Engineering project “Research on the Development Strategy of Engineering Science and Technology in the Field of Environment in the Next 20 Years” (2021-XBZD-13); National Social Science Fund of China (22AZD094)