[发明专利]用于复杂能源互联系统的韧性评估方法

权利要求书

1.一种用于复杂能源互联系统的韧性评估方法，其特征在于，包括如下步骤：

步骤1：分析能源互联系统在典型极端事件下的动态过程，获得在典型极端事件下能源互联系统及能源互联系统下设子系统的功能函数特征，所述子系统包括故障子系统和非故障子系统；

步骤2：基于所述能源互联系统及子系统的功能函数特征，建立能源互联系统韧性评价指标体系，所述韧性评价指标包括故障子系统对非故障子系统的感染指标、非故障子系统对故障子系统的反馈指标；

步骤3：通过抽样方法确定所述能源互联系统的典型极端事件场景及其概率；

步骤4：通过仿真方法获得各个典型极端事件下，能源互联系统及下设子系统的功能函数曲线；基于功能函数曲线，计算获得能源互联系统韧性评价指标的值；

步骤2中，所述韧性评价指标包括初始感染深度指标，通过如下公式获得，

式中，E_infection,0为系统的初始感染深度；F_inf,0(t)为初始感染后系统的功能函数；F(t_d)为t_d时刻系统的功能函数值；F(t_f)为t_f时刻系统的功能函数值；t_f为出现反馈的时刻；t_r为故障系统开始修复时刻；对于未出现反馈作用的场景可以认为t_f＝+∞；

步骤2中，所述韧性评价指标包括反馈次数与反馈感染深度指标，通过如下公式获得，

式中，N_feedback演化过程中出现的反馈次数；E_infection为系统的反馈感染深度；n_m→n代表子系统m对子系统n产生扰动激励的次数；对于未出现反馈的场景，N_feedback＝0且E_infection＝0；A表示故障子系统，B表示非故障子系统；

步骤1中，系统的功能函数为，系统可满足负荷占总负荷需求的比例，对于含多类负荷的系统，功能函数值按每类负荷经济价值进行加权计算，

式中，σ_k(t)表示t时刻可满足k类负荷占k类负荷总需求的比例；ε_k为k类负荷经济价值因子；F(t)为复杂能源互联系统的功能函数值；Ω为典型极端事件场景集合。

2.根据权利要求1所述的用于复杂能源互联系统的韧性评估方法，其特征在于：步骤2中，所述韧性评价指标包括吸收能力指标，通过如下公式获得，

E_absorption＝t_d-t_D＝Δt

式中，E_absorption为故障系统的吸收能力；t_d为故障系统开始出现功能损失的时刻；t_D为恶意攻击发生的时刻；Δt为故障系统维持功能平稳的时长。

3.根据权利要求1所述的用于复杂能源互联系统的韧性评估方法，其特征在于：步骤2中，所述韧性评价指标包括适应能力指标，通过如下公式获得，

式中，E_adaptation为系统的适应能力；F_A(t_d)为故障系统初次调整后的功能函数值；F_A(t₀)为系统正常时的功能函数值。