基于氢能设备多状态模型的电氢能源系统可靠性评估方法

本发明涉及电力系统可靠性评估，具体是一种基于氢能设备多状态模型的电氢区域综合能源系统可靠性评估方法。

背景技术：

1、随着制-储-运-用氢等技术的不断成熟，氢能系统与电、气、热等能源系统的耦合愈发紧密，电氢区域综合能源系统能够有效挖掘异质能源系统间的协调互补潜力，充分发挥氢能系统零碳、灵活运行优势，被认为是实现绿色低碳能源系统的重要组成部分。然而，目前针对电氢区域综合能源系统可靠性评估的相关研究极为匮乏。

2、在氢能设备可靠性建模方面，现有研究主要采用电解槽的两状态模型及发电模式下燃料电池的多状态模型，忽略了由于辅助元件故障导致电解槽制氢效率降低的非正常运行状态，以及热电联产模式下不同组件故障与燃料电池热电转化效率间的耦合关系。同时，尚未有研究建立起加氢设备的可靠性评估模型。因此，现有氢能设备可靠性模型难以准确反映电氢区域综合能源系统的可靠性水平。

3、在电氢区域综合能源系统可靠性评估方面，现有研究仅从光伏-热-燃料电池集成系统的设备故障风险优先级对和氢能系统对微电网供电可靠性的提升作用方面进行相关分析，可见电氢区域综合能源系统的供需失衡机理并不明晰，且系统内异质能源可靠性评估指标体系尚未完善。

4、综上，亟需从氢能设备机理分析出发，建立准确的氢能设备可靠性评估模型，并提出电氢区域综合能源系统可靠性评估指标体系及其评估方法。

技术实现思路

1、本发明的目的是解决现有技术中系统内异质能源可靠性评估指标体系尚未完善的技术问题。

2、本发明中的基于氢能设备多状态模型的电氢区域综合能源系统可靠性评估方法，主要包括以下步骤：

3、步骤1)对电氢区域综合能源系统进行描述，并对系统中电解槽、燃料电池和加氢机建立多状态马尔可夫空间，对其他能源设备建立两状态马尔可夫空间；

4、步骤2)根据所有设备的马尔可夫空间，采用马尔科夫链蒙特卡洛法抽样生成设备运行状态序列并构成系统运行状态；

5、步骤3)以弃风弃光惩罚成本与异质能源负荷的削减惩罚成本之和最小为目标函数，计及系统运行约束条件，建立系统最优异质能源负荷削减模型；

6、步骤4)采用异质能源负荷削减概率、异质能源负荷削减时间期望、异质能源负荷缺供能量期望、异质能源负荷缺供能量综合损失和设备可靠性贡献度构建异质能源可靠性评估指标体系；

7、步骤5)对所有系统运行状态下的最优异质能源负荷削减模型进行求解，计算异质能源可靠性评估指标，评估电氢区域综合能源系统可靠性水平。

8、进一步的，步骤1)中的电氢区域综合能源系统包含电能设备、氢能设备和热能设备；

9、所述的电能设备包括电力线路、风电机组、光伏机组和蓄电池；

10、所述的氢能设备包括电解槽、储氢罐、燃料电池和若干加氢机；

11、所述的热能设备包括燃气锅炉和燃料电池；

12、所述的燃料电池通过消耗来自储氢罐的氢气进行热电联产；

13、电氢区域综合能源系统的电负荷由外部电网、风光机组和蓄电池供给，氢负荷由来自电解槽和储氢罐释放的氢气通过若干加氢机供给；

14、热负荷由燃气锅炉消耗天然气产热和燃料电池消耗氢气产热供给。

15、进一步的，步骤1)中，对所述电解槽建立三状态马尔可夫空间，

16、其状态包括，正常运行状态；

17、由于碱性电解池或气体分离装置、循环水泵发生故障所引起的故障停运状态；

18、由热器管程发生故障引起的降额运行状态；

19、对所述燃料电池建立四状态马尔可夫空间，

20、其状态包括，正常运行状态；

21、由于燃料电池堆发生故障或氧气供给不足造成质子膜故障所导致的燃料电池故障停运状态；

22、由于换热器故障导致燃料电池工作温度升高导致的燃料电池输出电功率降低、输出热功率增加状态；

23、由于当加湿器故障时质子膜变干导致其欧姆电阻增加，进而导致的燃料电池输出电功率降低、输出热功率增加状态；

24、计及每台加氢机独立故障，对加氢机整体建立多状态马尔可夫空间，其状态包括，正常运行状态，

25、以及，由于不同数量加氢机同时发生故障所导致的多种故障运行/停运状态；

26、对所述其他能源设备均建立两状态马尔可夫空间，计及正常运行状态和故障停运状态。

27、进一步的，步骤2)中，生成系统运行状态具体包括以下步骤：

28、步骤2.1)确定设备x在初始时刻的设备状态sx,0；

29、步骤2.2)根据当前时刻的设备状态sx,t，根据其马尔可夫空间内各状态转移路径发生的概率进行随机抽样，确定下一时刻的状态sx,t+1；

30、步骤2.3)重复步骤2.2)，直到确定周期内所有时刻设备x的状态，生成设备x的时序状态序列；

31、步骤2.4)对电氢区域综合能源系统内所有设备重复步骤2.1)～2.3)，生成系统内所有设备的状态序列集合，即系统运行状态。

32、进一步的，步骤3)中，最优异质能源负荷削减模型的目标函数表达式为：

33、

34、式中，cdg为弃风弃光的单位惩罚成本；ωwg、ωpv分别为风光机组安装节点集合；△pwg,it、△ppv,it分别为弃风弃光功率；t为调度周期；n为节点数量；z为供给负荷类型集合，包括电、热、氢三类负荷；cz为第z类负荷的单位削减惩罚成本；lz,it为第z类负荷的削减功率，下标i、t分别表示节点i与t时刻。

35、进一步的，步骤3)中，系统运行约束条件包括功率平衡约束、负荷削减约束、设备运行约束和交互功率约束：

36、所述功率平衡约束包括电、热、氢功率平衡约束；

37、所述电功率平衡约束如下所示：

38、pgrid,it+pwg,it+ppv,it+pdis,it+pf,it＝pel,it+pch,it+pload,it-lp,it

39、

40、式中，pgrid,it为能源耦合节点与电网的交互功率；pwg,it、ppv,it分别为风光机组实际输出功率；pch,it、pdis,it分别为蓄电池的充放电功率；pfc,it为燃料电池的产电功率；pel,it为电解槽的耗电功率；pload,it为电负荷需求功率；lp,it为电负荷削减功率；pij,t为电力线路ij传输的有功功率；k∈j表示节点k以节点j为首端节点；ωeh为能源耦合节点集合；

41、所述热功率平衡约束如下所示：

42、hgb,it+hfc,it＝hload,it-lh,it

43、式中，hgb,it为燃气锅炉产热功率；hfc,it为燃料电池产热功率；hload,it为热负荷需求功率；lh,it为热负荷削减功率；

44、所述氢功率平衡约束如下所示：

45、

46、式中，qel,it为电解槽产氢量；qht,in,it、qht,out,it分别为储氢罐的储氢量和放氢量，qfc,it为燃料电池的耗氢量，qel2dp,it、qht2dp,it分别为电解槽和储氢罐供给加氢机的氢气量，qdp,in,t、qdp,out,t分别为所有加氢机输入、输出的氢气量，qload,it为氢负荷需求量，lq,it为氢负荷削减量；

47、所述负荷削减约束如下所示：

48、

49、所述设备运行约束包括电力线路运行约束、风光机组运行约束、蓄电池运行约束、燃气锅炉运行约束、电解槽运行约束、储氢罐运行约束、燃料电池运行约束和加氢机运行约束；

50、所述电力线路运行约束如下所示：

51、-sij,tpij,max≤pij,t≤sij,tpij,max，

52、式中，sij,t为线路ij的运行状态，1表示正常运行，0表示故障停运；pij,max为线路传输有功功率上限；

53、所述风光机组运行约束如下所示：

54、

55、式中，分别为风光机组预测的输出功率；swg,it、spv,it分别为风光机组的运行状态，1表示正常运行，0表示故障停运；

56、所述蓄电池运行约束如下所示：

57、

58、式中，soci,t、soci,t+1分别为t时刻和t+1时刻蓄电池的荷电状态，pch,it、pdis,it分别为蓄电池的充放电功率，ηch、ηdis分别为充放电效率，ses为蓄电池容量，αch,it、βdis,it分别为蓄电池充电和放电状态，0表示无电能充/放，1表示电能充/放，socmax、socmin分别为蓄电池荷电状态上下限约束，soc0、soct分别为研究周期内初始时刻、末尾时刻蓄电池的荷电状，；ses,it为蓄电池的运行状态，1表示正常运行，0表示故障停运；

59、所述燃气锅炉运行约束如下所示：

60、

61、式中，ggb,it为燃气锅炉耗气功率；ηgb为燃气锅炉效率；sgb,it为燃气锅炉运行状态，1表示正常运行，0表示故障停运；为燃气锅炉耗气功率上限；

62、所述电解槽运行约束如下所示：

63、

64、式中，qel,it为电解槽产氢量；ρh2为氢气的标准密度；k1、k2为单位转换系数；f为法拉第常数；ηf为法拉第效率；为单个电解槽处于状态sel,it下的运行电压；为电解槽处于状态sel,it下总耗电功率；μel,it为电解槽启停状态，1为启动状态，0为停机状态；pel,max和pel,min分别为电解槽总耗电功率上下限；

65、所述储氢罐运行约束如下所示：

66、

67、式中，sohi,t、sohi,t+1分别为t时刻和t+1时刻储氢罐的储氢状态，即该时刻储存氢气与容量的比值；qht,in,it、qht,out,it分别为储氢罐的储氢量和放氢量；ηht,in、ηht,out分别为充放氢效率；sht为储氢罐容量；αin,it、βout,it分别为储氢罐充氢和放氢状态，0表示无氢气充/放，1表示氢气充/放；sohmax、sohmin分别为储氢罐的储氢状态上下限约束；soh0、soht分别为研究周期内初始时刻、末尾时刻储氢罐的储氢状态；sht,it为储氢罐的运行状态，1表示正常运行，0表示故障停运；

68、燃料电池运行约束如下所示：

69、

70、式中，qfc,it为燃料电池的耗氢量，单位kg；mfc,it为燃料电池的耗氢量，单位kw；pfc,it、hfc,it分别为燃料电池的产电功率、产热功率；％为燃料电池处于状态sfc,it下输出电功率的降额系数，其中rp1％＝0、rp4％＝1；ufc,it为单个燃料电池的运行电压；kfc为燃料电池的热回收利用系数；△h为氢气的摩尔热值；μfc,it为燃料电池的启停状态，1为启动状态，0为停机状态；pfc,max、pfc,min分别为燃料电池产电功率上下限值；qfc,max、qfc,min分别为燃料电池消耗氢气量上下限值；

71、所述加氢机运行约束如下所示：

72、

73、式中，qdp,in,it、qdp,out,it分别为所有加氢机输入、输出的氢气量；ηdp为加氢效率；cdp,it分别为所有加氢机的可用容量；sdp,it为加氢设备的运行状态；cdp为单台加氢机的可用容量。

74、所述交互功率约束如下所示：

75、

76、式中，为能源耦合节点与电网交互功率的上限；sggrid,it为天然气网等效元件的运行状态，1表示正常运行，0表示故障停运；为能源耦合节点与天然气网交互功率的上限。

77、进一步的，步骤4)中，异质能源可靠性评估指标体系包括异质能源负荷削减概率、异质能源负荷削减时间期望、异质能源负荷缺供能量期望、异质能源负荷缺供能量综合损失和设备可靠性贡献度；

78、所述异质能源负荷削减概率指标如下所示：

79、

80、式中，xz,t为t时刻第z类负荷的削减状态变量，1表示系统进行负荷削减，否则为0；ny为仿真年限；

81、所述异质能源负荷削减时间期望指标如下所示：

82、lolez＝8760lolpz，

83、所述异质能源负荷缺供能量期望指标如下所示：

84、

85、所述异质能源负荷缺供能量综合损失指标如下所示：

86、

87、所述设备可靠性贡献度指标如下所示：

88、rcoex,z＝(eensz-eensz,x)/eensz，

89、式中，rcoex,z为设备x对第z类负荷缺供能量的贡献度指标；eensz,x为仅不考虑设备x故障场景下第z类负荷的异质能源负荷缺供能量期望指标。

90、进一步的，步骤5)中，评估电氢区域综合能源系统可靠性水平具体包括以下步骤：

91、步骤5.1)输入电氢区域综合能源系统及其设备参数相关数据，初始化所有设备的运行状态；

92、步骤5.2)以天为单位，将通过马尔可夫链蒙特卡洛模拟抽样得到的系统运行状态划分成不同断面；

93、步骤5.3)初始化断面计数d＝1和模拟年数ny＝1；

94、步骤5.4)获取第d个断面下系统状态，求解该断面下系统最优异质能源负荷削减问题，得到该断面每个时段的负荷削减状态和最优负荷削减量；

95、步骤5.5)若d>365，转向步骤6)，否则令d＝d+1，转向步骤5.4)；

96、步骤5.6)计算第ny年的系统可靠性评估指标，并令ny＝ny+1；

97、步骤5.7)若ny>15000或系统eens指标方差ε<0.05，转向步骤8)，否则令d＝1，转向步骤4)；

98、步骤5.8)输出系统异质能源可靠性评估指标。

99、与现有技术相比，本发明具体以下优点：

100、本发明提出的电解槽、热电联产燃料电池和加氢机的多状态模型，有助于准确评估电氢区域综合能源系统的异质能源负荷供给可靠性水平，并通过燃料电池热电联产模式提高系统能源利用率，降低了负荷缺供能量和系统经济损失。

101、本发明从不同角度出发构建异质能源可靠性评估指标体系，并结合最优异质能源负荷削减模型，基于马尔可夫链蒙特卡洛法对电氢区域综合能源系统的可靠性进行量化评估，揭示了故障模式下电氢区域综合能源系统内异质能源的供需失衡机理，有效量化了系统供能可靠性水平、经济损失以及设备的可靠性贡献程度，为系统的规划建设和安全运行提供了科学合理的决策依据。