第42卷第5期发电技术Vol.42No.52021年10月PowerGenerationTechnologyOct.2021DOI：10.12096/j.2096-4528.pgt.21010中图分类号：TK91;TQ116.2碳中和背景下高温固体氧化物电解制氢的过程建模与热力学分析王丹丹，李亚楼，李芳，孙璐（电网安全与节能国家重点实验室(中国电力科学研究院有限公司)，北京市海淀区100192）ProcessModellingandThermodynamicAnalysisofHydrogenProductionbyHighTemperatureSolidOxideElectrolysisUndertheBackgroundofCarbonNeutralityWANGDandan,LIYalou,LIFang,SUNLu(StateKeyLaboratoryofPowerGridSafetyandEnergyConservation(ChinaElectricPowerResearchInstitute),HaidianDistrict,Beijing100192,China)摘要：电解制氢能够将不稳定的可再生能源转化为氢能，实现大规模、季节性储能，是助力实现我国碳中和目标的关键手段之一。对高温固体氧化物电解制氢技术进行了探讨与研究，介绍了高温固体氧化物电解的工作原理；针对其实际运行特性，建立了电解过程电化学模型与热力学模型，并通过与实验数据对比验证了模型的有效性；进一步地，分析了电流密度和操作温度对高温固体氧化物电解极化损失的影响，研究了高温电解过程热力学性能的变化情况。结果表明：在高温下对固体氧化物电解进行整体热集成与热管理，有利于减少电能消耗；若将高温电解过程与其他工业过程相结合，利用低品位的工业余热或废热，则能够进一步提升高温固体氧化物电解的整体能效。关键词：碳中和；可再生能源；电解制氢；高温固体氧化物；电化学模型；热力学分析ABSTRACT:Waterelectrolysiscantransformunstablerenewableenergyintothehydrogenenergyandrealizelarge-scaleandseasonalenergystorage,whichisoneofthekeymeanstoachievethegoalof“carbonneutrality”inChina.Thetechnologyofhightemperaturesolidoxideelectrolysiscell(SOEC)wasdiscussedandstudied,andtheworkingprincipleofthehightemperaturesolidoxideelectrolysiswasintroduced.Accordingtoitsactualoperationcharacteristics,theelectrochemicalmodeland基金项目：国家重点研发计划项目(2018YFE0208400)；中国电力科学研究院自筹基金项目(5242002000RF)。ProjectSupportedbyNationalKeyResearchandDevelopmentProgramofChina(2018YFE0208400);ProjectofChinaElectricPowerResearchInstitute(5242002000RF).thermodynamicmodelofelectrolysisprocesswereestablished,andthemodelswereverifiedbycomparingwiththeexperimentaldata.Furthermore,theeffectsofcurrentdensityandoperatingtemperatureonthehightemperaturesolidoxideelectrolysiswereanalyzed,andthethermodynamicanalysiswasconducted.TheresultsshowthatthethermalintegrationandmanagementofSOECathightemperaturewillbeconducivetoreducethepowerconsumption.Ifthehigh-temperatureelectrolysisprocessiscombinedwithotherindustrialprocessesandutilizedthelow-gradeindustrialwasteheat,theoverallefficiencyofhigh-temperaturesolidoxideelectrolysisisexpectedtobefurtherimproved.KEYWORDS:carbonneutrality;renewableenergy;hydrogenproductionbyelectrolysis;hightemperaturesolidoxide;electrochemicalmodel;thermodynamicanalysis0引言2020年9月，中国正式提出“二氧化碳排放力争于2030年前达到峰值，努力争取2060年前实现碳中和”的目标。积极发展清洁能源有助于促进我国碳达峰、碳中和工作的加速进行，加快产业结构的优化。氢能作为一种低碳、高效、环境友好、应用场景丰富的绿色能源，被誉为21世纪的“终极能源”。将可再生能源弃电用于电解水制备氢能，能够为可再生能源季节性大规模消纳提供一种有效的解决方式，同时也是助力实现碳中和的关键手段[1-4]。第42卷第5期发电技术555高渗透率可再生能源为电解制氢提供了契机，有效推动了低碳氢能的规模化发展。当前我国对新能源的消纳能力不足，制约了可再生能源的规模化发展[5-6]。2018年中国可再生能源弃电量高达1020亿kWh，超过了三峡电站的总发电量[7]。电解水制氢可有效提高电力系统对可再生能源的消纳能力，从供给侧增强对化石能源的替代作用，并实现全生命周期碳减排。目前，电解水制氢技术可分为碱性电解水技术、固体聚合物电解水技术和高温固体氧化物电解水(solidoxideelectrolysiscell，SOEC)技术[8-9]。其中，SOEC技术的电解效率高，在高温(600~1000℃)条件下将原料水转化为氢气和氧气，能够实现电能和热能向化学能的高效转化，是极具前景的高效大规模制氢方法[10-11]。由于高温固体氧化物电解较传统低温电解更高效，发展潜力巨大，近十年来SOEC的研究以指数形式增长[12]。为了进一步优化SOEC的生产过程，提高其整体性能，对SOEC建立数值模型进行研究十分有必要。Spacil等[13]首次建立了SOEC的简单模型，分析了电解池性能与欧姆电势、开路电压等参数的关系。Ni等[14]对固体氧化物电解水制氢过程建立了零维模型，基于此模型进行了能量分析。Stoots等[15]开展了合成气的高温共电解过程的研究，并将模拟数据与实验数据进行了对比验证。Hawkes等[16]采用FLUENT自定义子程序，建立了高温固体氧化物电解模型，通过几何结构研究中性电势、气体产率以及电流密度在空间中的分布情况。本文对高温固体氧化物电解工作原理进行介绍，针对SOEC的实际运行特性，考虑由于各类极化现象导致的过电压，建立SOEC的电化学模型以及热力学模型；在此基础上，进一步分析电解过程的关键参数对高温固体氧化物电解极化损失的影响，并研究不同操作条件下高温电解过程热力学性能的变化情况。1高温固体氧化物电解模型建立1.1高温固体氧化物电解工作原理高温固体氧化物电解池由阴极、电解质、阳极和连接体构成，如图1所示。从反应过程上来讲，SOEC可以看作高温固体燃料电池(solidoxidefuelcell，SOFC)的逆运行。在高温下，将一定的直流电压施加在两侧电极上，阴极通入水蒸气和部分氢气(提供还原性氛围)，原料水在阴极被还原分解，产生氢气和氧离子，在电动势的驱动作用下，氧离子穿过致密的固体氧化物电解质层到达阳极，在那里失去电子并析出氧气。阴极电解质层阳极电输入热输入600~1000℃氧气蒸汽氢气蒸汽氢气图1高温固体氧化物电解过程Fig.1Hightemperaturesolidoxideelectrolysisprocess从热力学角度看，高温固体氧化物电解过程是水分解的过程，对于一个理想的电解过程(处于平衡状态下)，其理论所需能量为所需电能和热能之和[17]：HGTS(1)Nernst/()EGnF(2)式中：H为反应焓变，即水分解所需要的理论最小能量；G为吉布斯自由能，为电解过程所需要的最小电能；电解所需热能TS可由外部热源或电能提供，其中T为电解温度，S为反应熵变；n为电子转移数，取值为2；F为法拉第常数，其值为96485C/mol；ENernst为能斯特电压(电解水理论分电压)。标准状况下SOEC所需能量、能斯特电压与温度的关系如图2所示。可以看出，在100℃时，电解所需电能G占比较大，约占全部所需能量H的93%；随着温度进一步升高，H略有增加，G逐渐降低，所需热能TS逐渐增加，当温度升高到1000℃时，电能G占比降低为73%左右，更多电能能够被热能所替代；随着G的降低，能斯特电压ENernst也随之降低。1.2高温固体氧化物电解模型在实际电解过程中，反应往往偏离平衡状态。556王丹丹等：碳中和背景下高温固体氧化物电解制氢的过程建模与热力学分析Vol.42No.5各种不可逆损失会导致电极电位偏离平衡电位，这种现象被称为极化现象。极化现象导致实际工作电压E比电解水理论分电压Er高，从而使部分电能转化成热能。所需能量/(kJ·mol1)热能TS总焓变H电能G能斯特电压Enernst3002502001005001502.01.61.20.80.40.0能斯特电压/V温度/℃020040060080010001200图2SOEC所需能量、能斯特电压与温度的关系Fig.2Relationshipbetweenenergy,NernstvoltageandtemperatureforSOEC极化损失包括欧姆极化Eohm、活化极化Eact以及浓差极化Econ，其实际工作电压E可表示为rconactohmEEEEE(3)式中Er为SOEC的平衡电压，根据能斯特方程，可表示为2220.5HOrHO(,)ln()PPGRTETPnFnFP(4)式中：2HP、2HOP、2OP分别为氢气、蒸汽和氧气分压；R为理想气体常数，其值为8.3145J/(mol·K)。吉布斯自由能ΔG计算式为32()24480049.182.7210GTTT(5)1）浓差过电压浓差过电压Econ的产生是由各种气体分子在多孔电极中的传输、电解质内反应物和产物的溶解与分解、反应物在反应界面的扩散等所引起的[18]。Econ可以用电极表面与电极–电解质界面之间的气体浓度差表示，包括阴极浓差过电压Econ,c和阳极浓差过电压Econ,a，可分别表示为：22222222222222effHOHHHOcon,ceffHHOHHOeffHHOHeffHHOHO[/(2)]ln2[/(2)]1/(2)ln21/(2)PPJRTdFDRTEFPPJRTdFDJRTdFDPRTFJRTdFDP(6)2222OOgcon,aO+/(2)ln4PJRTdFBRTEFP（）(7)式中：2Hd、2Od分别为阴极和阳极的电极厚度；2effHOD为组分H2O的有效扩散系数；J为电流密度；为O2的动力学黏度；Bg为气体渗透率。2）活化过电压化学反应中，反应物需要克服能量势垒，即活化能。在固体氧化物电解过程中会发生电荷迁移的电化学反应，反应物需要克服热能势垒和电势，克服反应活化能垒的过电压为活化过电压。对于活化过电压Eact,i(ic,a分别表示阴极、阳极)，主要根据Butler-Volmer公式[19]计算：2act,0,0,ln[+()1]22iiiRTJJEFJJ(8)式中J0,i代表交换电流密度。对于阴极和阳极，其交换电流密度计算公式分别表示为：act,c0,ccexpEJRT（）(9)act,a0,aaexpEJRT（）(10)式中：c和a分别为阴极、阳极的指前因子；Eact,c和Eact,a分别为阴极、阳极的活化能。3）欧姆过电压电流在流经连接板、电极以及电解质时不可避免地产生欧姆极化，由于电解质电导远小于连接板与电极，因此欧姆过电压主要受电解质影响，其大小由电流密度J、电解质厚度dE和运行温度T决定，可表示为5ohmE103002.9910exp()EJdT(11)4）SOEC热力学模型在SOEC实际运行过程中，上述SOEC的过电压(不可逆)运行导致热能的产生。电解过程的熵产Sgen可表示为ohmconactgen()nFEEEST(12)电解过程吸/放热量Qst与熵产Sgen直接相关：stoutingen()QTSSTS(13)式中Sout、Sin分别为输出熵和输入熵。第42卷第5期发电技术557此外，电解过程实际消耗电能Pst以及氢气产率2Hn分别表示为：stSOECPEJAN(14)2SOECH2JANnF(15)式中：ASOEC为单个电解池面积；N为SOEC电解池数量。可以看出，通过电解的电化学模型，能够定量获得不同条件(电流密度、温度、入口蒸汽分压)下电解槽的电流密度–电压特性曲线；通过进一步计算，可得到电解槽工作在不同操作条件下所需的热量和所消耗电能。2结果分析与讨论2.1模型验证本文借助MATLAB/Simulink软件，对SOEC制氢建立了电化学模型和热力学模型。电解模型的操作条件和结构参数见表1。表1电解模型的操作条件和结构参数Tab.1Operatingconditionsandstructuralparametersofelectrolysismodel操作条件和结构参数数值操作压力P/kPa100操作温度T/K973阴极入口成分及其体积分数100%O2阳极入口成分及其体积分数90%H2O，10%H2电子转移数n2阴极厚度dH2/μm312.5阳极厚度dO2/μm17.5电解质厚度dE/μm12.5阴极指前因子c/(GA/m2)0.39阳极指前因子a/(GA/m2)1.39阴极活化能Eact,c/(MJ/mol)0.1阳极活化能Eact,a/(MJ/mol)0.12孔隙曲度5孔隙率0.3孔隙平均半径r/m0.5瑞士洛桑联邦理工学院的SOEC模型[20]电压实验值与模拟值对比如表2所示。通过再现电解池的电流密度–运行电压特性，验证了SOEC建模的有效性。从表2可以看出，在SOEC运行温度923K和973K条件下，当电流密度从0A/m2升高至8000A/m2时，模型电压模拟值与实验测量值较为吻合，证明该SOEC模型具有可行性。表2SOEC模型电压模拟值与实验值对比Tab.2ComparisonofvoltagesimulationvaluesandexperimentalvaluesofSOECmodel电流密度/(A/m2)电压模拟值/V电压实验值/V923K973K923K973K00.9340.9150.8950.88010001.1311.0331.0640.97220001.2491.1221.2361.07730001.3371.1931.3851.17840001.4101.2521.5101.27750001.4741.3031.5801.36060001.5321.3501.6201.42070001.5861.3921.6551.46680001.6361.4321.6801.5102.2关键参数对极化损失的影响极化过电压Eover包括活化极化Eact、欧姆极化Eohm和浓差极化Econ。其中，Eact是主要的不可逆极化损失，其次是Eohm和Econ。各类极化过电压分布情况如图3所示。Eover/VT1073K电流密度/(A/m2)0200040006000800010000EconEohmEact0.300.250.200.150.100.050.00图3各类极化过电压分布情况Fig.3DistributionofvariouspolarizationovervoltageEact由电极动力学决定，在化学物质和电极之间传输的电荷需要克服此电势。电流密度、温度对活化过电压Eact的影响如图4所示，可以看出，随着电流密度从0A/m2升高至10000A/m2，Eact逐渐变大。在973K运行温度下，Eact从0V增长至0.4V，这主要是由阴极活化极化的影响所致；当温度高于973K时，随着电流密度的增加，Eact呈线性增加。在相同电流密度下(5000A/m2)，随着SOEC的运行温度从873K增长至1373K，Eact大幅度降低，由0.49V降低至0.01V。558王丹丹等：碳中和背景下高温固体氧化物电解制氢的过程建模与热力学分析Vol.42No.5Eact/V电流密度/(A/m2)0200040006000800010000873K973K1073K1173K1273K1373K0.60.50.40.30.20.10.0图4电流密度、温度对活化过电压的影响Fig.4Influenceofcurrentdensityandtemperatureonactivationovervoltage电流密度、温度对欧姆过电压Eohm的影响如图5所示，可以看出，Eohm随电流密度的增大而呈线性增加。在973K运行温度下，当电流密度从0A/m2升高至10000A/m2时，Eohm从0V增长至0.148V。Eohm受温度影响较大，在相同的电流密度下(5000A/m2)，随着SOEC运行温度从873K升高至1373K，Eohm由0.249V降低至0.030V，这主要是由于离子电导率和电极反应随着温度升高而增加。Eohm/V电流密度/(A/m2)0200040006000800010000873K973K1073K1173K1273K1373K0.50.40.30.20.10.0图5电流密度、温度对欧姆过电压的影响Fig.5Influenceofcurrentdensityandtemperatureonohmicovervoltage浓差极化分为阴极和阳极两侧，阴极的浓差极化占据主导地位。电流密度、温度对浓差过电压Econ的影响如图6所示，可以看出，在973K运行温度下，随着电流密度增加，Econ从0V升高至0.037V。这是由于阴极支撑层较厚，随着电流密度增加，反应物消耗增加，蒸汽供给不足，扩散速度缓慢。在相同电流密度下(5000A/m2)，相比于活化极化和欧姆极化，温度对浓差极化的影响较小。与上述2种极化方式不同，浓差极化随运行温度升高而增加，这是由于摩尔扩散率随着温度的升高而降低，导致了更高的浓差过电压。Econ/V电流密度/(A/m2)0200040006000800010000873K973K1073K1173K1273K1373K0.070.060.050.040.030.020.010.00图6电流密度、温度对浓差过电压的影响Fig.6Influenceofcurrentdensityandtemperatureonconcentrationovervoltage随着运行温度上升，活化极化和欧姆极化的下降幅度显著大于浓差极化的升高幅度，因此，总的过电压随着温度上升而降低。进一步地，SOEC的总电解电压会随着温度上升而下降。2.3SOEC热力性能分析为了明确SOEC运行时所需的热能和电能，计算因极化导致的不可逆损失产热以及实际电能消耗是非常重要的。图7、8分别为电流密度、温度对SOEC运行时实际所需热能和电能的影响。在相同运行温度下(1023K)，随着SOEC电流密度从1000A/m2增加至5000A/m2，制备1mol氢气所需热能由74.9kJ降低至30.1kJ，所需电能由176.6kJ提高至217.7kJ。这主要是由于极化损所需热能/kJ电流密度/(A/m2)973K1023K1073K0100020003000400060005000806040200图7电流密度、温度对SOEC所需热能的影响Fig.7InfluenceofcurrentdensityandtemperatureonheatenergyrequiredbySOEC第42卷第5期发电技术559所需电能/kJ电流密度/(A/m2)973K1023K1073K0100020003000400060005000260240220200180160图8电流密度、温度对SOEC所需电能的影响Fig.8InfluenceofcurrentdensityandtemperatureonpowerrequiredbySOEC失造成的熵产增大，使得部分电能转化为热能，从而导致SOEC热量需求降低、电量需求增加。在相同电流密度下(3000A/m2)，随着SOEC的运行温度由973K增加至1073K，制备1mol氢气所需热能从31.0kJ增加至62.8kJ，所需电能从216.5kJ降至185.4kJ。可以看出，如果在高温下合理采取系统集成方式，将过程热能应用于电解，则能有效降低电能的消耗。3结论电解水制氢将电能转化为氢气进行存储，能够极大地拓宽弃电的应用范围。从这个角度出发，对高温固体氧化物电解制氢技术进行了探讨与研究，介绍了高温固体氧化物电解的工作原理，建立了电解过程电化学模型与热力学模型，得到以下结论：1）各类极化过电压都随电流密度增大而升高；电解过程的总体过电压随着温度上升呈下降趋势。2）随着电流密度增大，部分电能转化为热能，导致电解热量需求显著降低；随着运行温度的增加，SOEC热量需求增加，电量需求则逐渐降低。3）在高温下对SOEC进行整体热集成与热管理，有利于减少电能消耗，若再利用低品位的工业余热或废热，则能够进一步提升高温固体氧化物电解的整体能效。为了实现我国“2060年碳中和”目标，大力发展氢能产业是优选路线。高温固体氧化物电解制氢技术制氢效率较高，随着技术的成熟，有望在化工、分布式能源等领域获得较为广泛的应用。研究结果为进一步探索高温固体氧化物电解过程与其他高温工业过程相结合提供了参考，为后续的高温制氢系统集成与优化提供了理论依据。参考文献[1]张文强，于波．高温固体氧化物电解制氢技术发展现状与展望[J]．电化学，2020，26(2)：212-229．ZHANGWQ，YUB．Developmentstatusandprospectsofhydrogenproductionbyhightemperaturesolidoxideelectrolysis[J]．JournalofElectrochemistry，2020，26(2)：212-229．[2]牟树君，林今，邢学韬，等．高温固体氧化物电解水制氢储能技术及应用展望[J]．电网技术，2017，41(10)：3385-3391．MUSJ，LINJ，XINGXT，etal．Technologyandapplicationprospectofhigh-temperaturesolidoxideelectrolysiscell[J]．PowerSystemTechnology，2017，41(10)：3385-3391．[3]雷超，李韬．碳中和背景下氢能利用关键技术及发展现状[J]．发电技术，2021，42(2)：207-217．LEIC，LIT．Technologiesanddevelopmentstatusofhydrogenenergyutilizationunderthebackgroundofcarbonneutrality[J]．PowerGenerationTechnology，2021，42(2)：207-217．[4]严思韵，王晨，周登极．含氢能气网掺混输运的综合能源系统优化研究[J]．电力工程技术，2021，40(1)：10-17．YANSY，WANGC，ZHOUDJ．Optimizationofintegratedelectricityandgassystemconsideringhydrogen-natural-gasmixturetransportation[J]．ElectricPowerEngineeringTechnology，2021，40(1)：10-17．[5]周强，汪宁渤，冉亮，等．中国新能源弃风弃光原因分析及前景探究[J]．中国电力，2016，49(9)：7-12．ZHOUQ，WANGNB，RANL，etal．CauseanalysisonwindandphotovoltaicenergycurtailmentandprospectresearchinChina[J]．ElectricPower，2016，49(9)：7-12．[6]陈建明，肖佳璇．基于弃风弃光问题的氢储能可行性研究[J]．技术与市场，2019，26(11)：23-25．CHENJ，XIAOJX．Feasibilitystudyofhydrogenstoragebasedontheproblemofabandonedwindandlight[J]．TechnologyandMarket，2019，26(11)：23-25．[7]中国电力企业联合会．2018年全国电力工业统计快560王丹丹等：碳中和背景下高温固体氧化物电解制氢的过程建模与热力学分析Vol.42No.5报一览表[EB/OL]．(2018-01-31)[2021-01-02]．https://max.book118.com/html/2019/0520/8056127046002024.shtm．ChinaFederationofElectricPowerEnterprises．Listofnationalelectricpowerindustrystatisticsexpressin2018[EB/OL]．(2018-01-31)[2021-01-02]．https://max.book118.com/html/2019/0520/8056127046002024.shtm．[8]CHIJ，YUH．Waterelectrolysisbasedonrenewableenergyforhydrogenproduction[J]．ChineseJournalofCatalysis，2018，39(3)：390-394．[9]俞红梅，衣宝廉．电解制氢与氢储能[J]．中国工程科学，2018(3)：58-65．YUHM，YIBL．Hydrogenforenergystorageandhydrogenproductionfromelectrolysis[J]．StrategicStudyofCAE，2018(3)：58-65．[10]霍现旭，王靖，蒋菱，等．氢储能系统关键技术及应用综述[J]．储能科学与技术，2016，5(2)：197-203．HUOXX，WANGJ，JIANGL，etal．Reviewonkeytechnologiesandapplicationsofhydrogenenergystoragesystem[J]．EnergyStorageScienceandTechnology，2016，5(2)：197-203．[11]SCHIEBAHNS，GRUBET，ROBINIUSM，etal．Powertogas：Technologicaloverview，systemsanalysisandeconomicassessmentforacasestudyinGermany[J]．InternationalJournalofHydrogenEnergy，2015，40(12)：4285-4294．[12]O’BRIENJE，STOOTSCM，HERRINGJS，etal．Performancemeasurementsofsolid-oxideelectrolysiscellsforhydrogenproduction[J]．JournalofFuelCellScienceandTechnology，2005，2(3)：156-163．[13]SPACILHS，TEDMONCS．Electrochemicaldissociationofwatervaporinsolidoxideelectrolytecells[J]．JournaloftheElectrochemicalSociety，1969，116(12)：1618．[14]NIM，LEUNGM，LEUNGD．Energyandexergyanalysisofhydrogenproductionbysolidoxidesteamelectrolyzerplant[J]．InternationalJournalofHydrogenEnergy，2007，32(18)：4648-4660．[15]STOOTSCM，OBRIENJE，HERRINGJS，etal．Syngasproductionviahigh-temperaturecoelectrolysisofsteamandcarbondioxide[J]．JournalofFuelCellScience&Technology，2009，6(1)：1-12．[16]HAWKESGL，O'BRIENJE，STOOTSCM，etal．CFDmodelofaplanarsolidoxideelectrolysiscellforhydrogenproductionfromnuclearenergy[J]．NuclearTechnology，2005，158(2)：132-144．[17]ATKINSPW，PAULAJCD．PhysicalChemistry[M]．7thed．Oxford：OxfordUniversityPress，2002．[18]陆玉正．太阳能与固体氧化物电解池联合制氢关键技术的研究[D]．南京：东南大学，2017．LUYZ．Researchonkeyissueonhydrogenproductioncombingsolarenergyandlowtemperatureofsolidoxideelectrolysiscells[D]．Nanjing：SoutheastUniversity，2017．[19]EDUARDOHP，SINGHD，HUTTONPN，etal．Amacro-levelmodelfordeterminingtheperformancecharacteristicsofsolidoxidefuelcells[J]．JournalofPowerSources，2004，138(1/2)：174-186．[20]DIETHELMS，HERLEJV，MONTINAROD，etal．Electrolysisandco-electrolysisperformanceofSOEshortstacks[J]．FuelCells，2013，13(4)：631-637．收稿日期：2021-02-18。作者简介：王丹丹(1991)，女，博士，中级工程师，研究方向为综合能源与绿色氢能，wangdandan_iet@163.com；王丹丹李亚楼(1974)，男，博士，教授级高工，研究方向为电力系统分析及控制，综合能源系统基础理论和关键技术等，本文通信作者，liyalou@epri.sgcc.com.cn；李亚楼李芳(1979)，女，博士，教授级高工，长期从事电力系统仿真、综合能源优化研究，lifang@epri.sgcc.com.cn；李芳孙璐(1982)，女，硕士，中级工程师，研究方向为综合能源与电力系统仿真，sunlu1@epri.sgcc.com.cn。孙璐(责任编辑尚彩娟)