负荷重分配攻击下电–气系统损失评估与脆弱性分析

周步祥，闵昕玮，臧天磊*，张远洪，陈 阳，赵雯雯

(1.四川大学 电气工程学院，四川 成都 610065；2.智能电网四川省重点实验室（四川大学），四川 成都 610065)

智能电网作为典型的信息物理系统[1],承担了推动建设新型电力系统的重要任务。然而，信息物理的高度融合也使系统的安全风险急剧增加。2015年，乌克兰电网遭受信息攻击，造成其全国大面积停电，停电持续时间之长，影响范围之大，十分罕见[2]。由此可见，信息攻击引发的安全事件会造成严重的经济损失和恶劣的社会影响。同时，伴随着燃气发电技术的快速发展，电网与天然气网耦合程度不断加深[3]。天然气网和电网具有很强的相互依存性，针对天然气网的信息攻击会影响电网的安全运行[4]。因此，全面研究信息攻击机理与影响对电–气耦合系统（integrated electricity natural gas system，IEGS）安全稳定运行意义重大。

现有研究中，针对电力系统的信息攻击主要有以下两大类。第1类是信息–物理协同攻击，通过物理方法破坏电网中元件的同时发动信息攻击，以掩盖电网的真实运行状态，对电网造成严重的损失。田猛等[5]以攻击者和调度中心的目标冲突为基础建立双层规划模型，并针对电网中的脆弱性节点进行分析。伏坚等[6]提出防御性随机规划模型来应对协同攻击带来的风险。阮振等[7]建立了以电力系统最大化损失为目标的双层模型，并给出了节点在信息攻击下的脆弱性分析。以上研究虽然为保护系统信息安全提供了有效的支撑，但并不能直观反映信息攻击对于系统造成的损失，同时也未对系统内节点的脆弱性进行全面的评估，故不能真正反映系统内的高脆弱性节点。此外，协同攻击需要物理攻击的时效性较高，因此在实际中难以对电力系统实施信息物理协同攻击。第2类是单一信息攻击，在实际电力系统的攻击中较为常见。本文主要对信息攻击中的虚假数据注入攻击展开深入研究。虚假数据注入攻击可分为拓扑攻击和负荷重分配攻击（load redistribution，LR）。拓扑攻击通过篡改线路开断数据，误导调度中心的状态感知，危害系统安全。王超超等[8]建立了电力信息物理系统分层模型，考虑了多层耦合的特性，并基于拓扑分析方法对攻击进行了模拟。Liang等[9]提出线路切换的新拓扑攻击方式，采用元启发式优化算法求解攻击模型，但其提出的拓扑攻击方式会对系统造成严重的损失。相比于拓扑攻击，负荷重分配攻击具有更强的隐蔽性，在一定阈值内实施攻击不易引起调度中心的警觉。因此本文主要针对负荷重分配攻击进行研究。Liu等[10]提出虚假数据注入攻击的概念，其构造的攻击可以绕过当时的不良数据检测技术。Yuan等[11]首次提出一种改进虚假数据注入攻击，即LR攻击。LR攻击是一种特殊类型的虚假数据注入攻击，它以系统支路虚假过载为基础干扰安全经济调度。舒隽[12]、Xiang[13–14]等以最大化系统运行成本为目标进行研究，其研究表明虚假数据注入双层优化攻击可能严重危害电力系统安全和经济运行。然而，上述研究未考虑攻击资源的限制及攻击引起的潮流转移会使线路过载，继而导致严重的连锁故障，给电力系统运行造成更大的损失。Che等[15–16]研究了LR攻击对电网造成级联故障的风险并识别了高风险线路，但未能将攻击给电力系统造成的总损失统一度量。以上关于LR攻击的研究均侧重于在电力系统分析攻击效果，且鲜有针对IEGS的信息攻击。Zhao等[4]首次在IEGS中对虚假数据注入攻击进行建模，但其主要研究单一能源子系统被攻击，并未涉及IEGS被协同攻击的情况。Zhao等[17]定义了天然气网中的虚假数据注入攻击，提出了一种日前IEGS协同运行方案来抵御信息攻击的影响，但其主要关注天然气网侧发生的信息攻击，并未考虑可能发生在电网侧的信息攻击。曹茂森等[18]提出了计及连锁故障的电–气耦合系统信息物理协同攻击策略，但其仅在电网侧进行了LR攻击，且未考虑LR攻击资源的限制，也未对系统的脆弱性节点进行分析，不能清晰地反映LR攻击对于IEGS的影响。

综上，本文详细分析在电力–天然气信息物理深度耦合下，双侧协同负荷重分配攻击的机理及流程；计及连锁故障对电力系统造成的影响，建立了LR攻击下的电–气耦合系统损失评估模型，统一度量不同能源子系统中的经济损失，并直观反映双侧协同LR的攻击效果；定义一种新的节点综合脆弱性评估指标，并在不同攻击场景下分析节点的脆弱性变化，进而提出针对不同攻击下的IEGS综合保护策略。仿真验证了保护策略可以有效提高系统抵御攻击的能力，减少系统经济损失。

1 IEGS下LR攻击机理分析

天然气网与电网通过燃气机组紧密耦合，攻击可以在电–气耦合系统中协同篡改电网和气网的负荷量测数据，诱导不同的调度中心做出错误的调整策略，引起支路或管道过载，从而转移发电机组和气源出力，并且导致系统甩负荷，造成严重的经济损失，危害IEGS安全经济运行。

由图1可知，双侧LR协同攻击过程如下：

1）电网侧 LR 攻击。攻击者在有限攻击资源的约束下先篡改电网侧负荷数据，诱导电网调度中心调整发电机出力和削减负荷，并将调整后的燃气机组出力传导到气网侧参与调度。

2）气网侧 LR 攻击。燃气机组输出功率调整完毕后，攻击者再篡改气网侧负荷量测值，诱导气网侧调度中心做出错误决策，削减气转电负荷和气负荷，使天然气系统处于非经济运行状态。

3） 燃气机组出力更新。由于气网侧的攻击会导致气转电负荷严重削减。确定削减后气转电负荷大小，并根据气电转化效率更新燃气机组输出功率。

4）电网侧连锁故障。气转电负荷严重不足会使燃气机组无法参与电网功率自平衡，而燃气机组出力的降低导致电网侧潮流重新平衡[22]，部分过载线路被切除，通过调整燃煤机组出力和切除电力负荷实现功率自平衡。

为突出双侧协同LR攻击效果，本文基于文献[9–21]中的前提假设对部分条件进行简化处理，该处理充分考虑了信息攻击的实际情况及其对IEGS所产生的实质影响[23]。本文在攻击的损失评估中做如下假设：

1）为了使信息攻击的效果更加明显，凸显IEGS的运行风险和脆弱节点，以大多数文献的共通场景为基础进行建模和测试[5–7,11–17]。假设攻击者知晓IEGS的结构参数（如节点导纳矩阵、节点支路矩阵等）和运行方式，并能成功注入虚假数据，躲过状态估计的检测。

2）天然气管存可能会对燃气机组输出功率下降起到缓解作用，但考虑到天然气调度时间尺度远远大于电网调度时间，且管存余量相较于气转电负荷缺口相差甚大，因此可认为管存对保持燃气机组输出功率效果甚微。

3）假设攻击者掌握足够的攻击资源，协同篡改线路和管道潮流量测值作为辅助配合攻击消耗部分资源后，攻击者仍有足量攻击资源用于篡改节点负荷量测值。

2 LR攻击的损失评估模型

负荷损失是LR攻击和连锁故障的直观体现，而由于负荷损失带来的运行成本升高能更真实地反映攻击的危险程度。因此，本文采用攻击后IEGS的总经济损失，统一度量气网侧与电网侧的负荷损失，表征系统抵御攻击的能力。系统经济损失越小，则系统抵御攻击的能力越强。基于攻击者视角建立的LR攻击损失评估模型如图2所示。

2.1 攻击模型

攻击者在保证LR攻击可以躲过IEGS不良数据检测，知晓IEGS系统调度策略的前提下，以IEGS经济损失最大为目标实施攻击。

2.1.1 攻击模型的目标函数

式（1）计及电网连锁故障对系统负荷的影响，采用经济损失统一度量电网与天然气网的负荷损失，表示攻击模型的目标函数。

式中，c1、c2、c3分别为单位发电机组发电成本、电网单位切负荷成本、天然气单位切负荷成本，NG、NLD、NGD分别为电网中发电机组节点数、电网中负荷节点数、天然气网中负荷节点数，s为电网侧连锁故障的阶段数，Ns为连锁故障传播次数。

2.1.2 攻击约束

LR攻击可在躲过IEGS不良数据检测的同时不引起调度中心的警觉，为此本文建立LR攻击约束的数学模型，如式（2）～（8）所示。式（2）、（3）分别表示电网、气网负荷节点负荷量测值的篡改量代数和为0，保证各系统有功功率和天然气质量平衡。式（4）表示双侧协同攻击下的LR攻击变量。篡改量测值相对于原量测值的比值大小与攻击资源消耗正相关，且近似满足一次相关性[7]，在限定信息攻击资源一定时，LR攻击变量满足约束如式（5）所示。式（6）限制了攻击对单个节点负荷量测值的篡改幅度，以确保攻击可以躲过状态估计的检测。式（7）、（8）分别表示攻击对电网支路、气网管道量测值的篡改量。

式（4）～（8）中， ∆Ld为整个IEGS各节点负荷量测值的篡改量， τ、λ、Rc分别为各负荷节点的攻击难度、单个负荷节点篡改量的阈值上限、LR攻击篡改负荷量测值资源上限， ∆F、∆fmn分别为电网线路和气网管道需要篡改的量测值，S1、S2分别为电网、气网转移因子矩阵，Abl、Gbl分别为电网支路负荷关联矩阵和气网管道负荷关联矩阵。若节点b为支路l的起点，则Abl=Gbl=1； 若母线b为支路l的终点，则Abl=Gbl=−1； 否则，Abl=Gbl=0。

2.2 IEGS调度模型

调度人员在满足IEGS中功率平衡、运行约束及系统安全约束的条件下，对IEGS中电网和气网分别进行安全经济调度。

2.2.1 IEGS调度模型的目标函数

IEGS调度模型目标函数如式（9）所示，以最小化电网和气网运行成本为目标分别进行调度。

2.3 电网连锁故障模型

攻击会使燃气机组输出功率下降或电网线路过载，导致电网中燃煤发电机组的出力有较大的调整，极易引起连锁故障，危害电网安全[18]。

电网中的连锁故障会使电网拓扑分解为孤岛状态运行，通过电网内的自动发电控制和低频减载装置可以自动平衡孤岛内的机组出力和负荷需求。由于天然气系统运行时间尺度远大于电网，当气网侧的攻击导致气转电负荷不足时，无法在较短时间内调整燃气机组出力，因此只有燃煤机组参与孤岛功率平衡[22]。孤岛平衡策略主要依据系统备用容量v决定。

调整孤岛内燃煤发电机出力和负荷，得到各节点的净注入功率，运用直流潮流模型计算各支路潮流后，系统功率自平衡。潮流的转移可能会使部分线路越限，从而导致线路停运。如果线路上的潮流大于其上限时，就将其切除，功率自平衡过程如式（28）～（31）所示。阶段s时的节点功率平衡约束为：

3 求解方法

3.1 线性化

3.1.1 连锁故障中支路潮流的线性化

连锁故障中的支路潮流是一个非线性函数，通过引入连续变量φ1～φ4和0–1变量b1进行线性化，将约束式（30）线性化得到式（32）：

式中， θ1、θ2分别为l线首端和末端的相位角，θmax、θmin分别为最大和最小相位角。

3.1.2 管道流量线性化

天然气网相关节点间的管道流量和压差公式是非线性的。事实上，天然气管道中的流量与气源出力和气负荷相关，其传输损耗较小。为了实现式（15）的线性化，本文利用文献[24]的方法对非线性项进行线性化处理，基于天然气网节点质量平衡计算管道流量。式（33）表示对天然气管道流量上限进行了约束，可以有效确保天然气网安全运行。

式中，fmn,max为管道mn流量上限。采用式（33）替代式（15）和（20）参与模型求解。

根据获取的管道流量值和预设的松弛节点压力（选择恒压气源作为松弛节点），可以逐一获取节点压力，可表示为：

求解损失评估模型时，不直接考虑节点压力约束，而是在输出最优解之前进行检验，本文方法在遵循节点气压约束的前提下，简化了损失评估模型的求解。

3.2 Karush–Kuhn–Tucher（KKT）算法

在天然气网管道潮流线性化后，对IEGS调度模型通过KKT条件进行转化，使电–气耦合系统损失评估模型转化为一个带平衡约束的数学规划问题，并将KKT条件作为一组额外的互补约束添加到攻击模型中。KKT算法中的广义拉格朗日函数如下：

式中，f(x)为 IEGS调度模型的目标函数，ci(x)和hj(x)分别为IEGS调度模型中的不等式和等式约束， αi和 βj分别为拉格朗日不等式和等式乘子，k和l分别为IEGS调度模型中不等式和等式约束的个数，x为IEGS调度模型中的变量。由函数推导的KKT条件如下：

采用KKT算法时所引入的拉格朗日乘子会使式（36）中的约束非线性化，借鉴文献[25]的大M法使其线性化。之后，通过式（33）代替式（15）和（20）、式（36）代替IEGS调度模型的约束条件，对电–气耦合系统损失评估模型进行求解；并在求解后，根据式（34）计算和判断节点压力，以对求解可行性进行检验。若可行，则输出结果；否则，调整式（33）参数重新求解。

4 算例分析

本文利用IEEE 39节点电力系统和改进比利时20节点天然气系统组成的电–气耦合系统验证上述损失评估模型。此外，调度中心与发电机组之间的信息联系十分紧密，发电机组的输出功率数据是无法被攻击者篡改的，但这种紧密的信息联系会随着负荷与发电机组距离的增加而减弱[15]，因此，依据系统负荷距离发电机或气源的远近来确定负荷节点的攻击难度，由参数 τ表示[7]，为了更好对比电–气耦合系统所有负荷节点变化情况，现将电、气负荷节点进行统一编号，电网侧负荷节点的编号依次为1～18，气网侧气负荷节点的编号依次为19～25，具体参数见表1。

各母线的负荷损失代价为100 $/MWh，各管道的气负荷损失代价为150 $/Mm3。本文在合理化假设的前提下为了凸显攻击效果，将单个负荷节点篡改量的阈值 λ设置为0.5[6,7,15]。攻击者若成功发动攻击，则攻击者应掌握足够的攻击资源，篡改线路潮流和管道流量所需的攻击资源不会影响篡改各负荷节点攻击资源的分配，且本文重点研究负荷重分配攻击对IEGS的影响，因此，在算例测试中，暂时不计入篡改线路潮流值所需的攻击资源。本文利用MATLAB实现上述过程，并调用CPLEX求解器进行求解。计算机硬件平台为Inter Core i5–8300H处理器，主频2.3 GHz，内存8 GB。

4.1 不同LR攻击方式对系统运行影响分析

研究3种不同攻击方式对电–气耦合系统运行情况的影响。攻击1，仅从电网侧进行LR攻击；攻击2，仅从气网侧进行LR攻击；攻击3，综合考虑电网侧与气网侧的协同LR攻击。

定义Rc为LR攻击资源上限，选取信息攻击资源为 0 .7Rc时，分析不同攻击方式对于系统的破坏效果。攻击1、2、3分别对IEGS造成的经济损失为21.99×103、29.87×103、33.52×103$。不同攻击形式下系统各节点损失负荷情况如图3所示。

由图3可知，不同攻击方式对于系统运行情况的破坏程度有较大差异。总的来说，采用攻击3对系统的影响最大，攻击2次之，攻击1对系统影响程度最小。说明了在攻击1下，天然气机组具有灵活调节能力，紧急情况下调度中心会优先增加燃气机组的输出，有效解决了电网功率缺额的问题，降低了系统的经济损失；同时，在攻击2、3下，当天然气系统受到攻击或者阻塞时，会优先切断气转电负荷，使燃气机组丧失灵活调节的能力。所以，双侧协同LR攻击会削减增大的气转电负荷，导致燃气机组的输出功率大大降低，不能参与电网功率平衡，因此影响电网侧的安全裕度，从而引发连锁故障，使电网负荷有较大的损失。而相比于攻击2，攻击1并不影响系统的备用容量，则连锁故障后形成的孤岛也更易满足系统功率自平衡的需求，所以攻击1对于系统运行情况的影响最小。

值得注意的是，攻击1中电网侧的部分节点负荷损失比例较高，而在其他两种攻击形式中并未发现节点负荷完全失去的现象。如节点3、4中损失负荷比例达到100%，说明了采用攻击1时，电网侧攻击资源充足，引起更多线路过载，切除过载线路后易形成不含发电机的孤岛，导致节点负荷被全部切除。同时，虽然攻击3会使电网系统的负荷损失大于攻击2的损失，但在气网侧对于气负荷的削减量却小于攻击2，这是因为攻击资源总量的限制导致攻击3在气网侧分配的资源小于攻击2，所以攻击2会使气网系统的经济损失更大。

4.2 不同攻击资源对系统运行影响分析

为研究不同攻击资源对系统运行情况的影响，比较总攻击资源为 0.5Rc、 0 .6Rc、0.7Rc和 0 .8Rc下系统总经济损失情况。不同攻击资源对于系统的影响情况如图4所示。

由图4可以看出，同一攻击方式下不同攻击资源时对系统的经济损失具有一定的递进关系，攻击资源越大，对系统造成的经济损失越大。如在攻击资源分别为 0.5Rc、0.6Rc、0.7Rc和 0.8Rc时，攻击2对系统造成损失分别为23.51×103、26.13×103、29.87×103、34.69×103$。显然攻击资源为 0.8Rc时攻击对系统运行影响最大。对比同一攻击资源下不同的攻击形式可知，当攻击资源较充足时，攻击3对系统造成的损失与攻击2的差值越来越小，如在攻击资源为 0 .8Rc时，攻击2、3对系统造成的经济损失分别为34.69×103$和35.68×103$。表明了攻击2造成气网侧经济损失较多的同时，攻击3中由于总攻击资源的限制，电网侧安全经济调度时燃气机组输出功率增加有限，导致气转电负荷削减量与攻击2接近。因此，攻击2、3在攻击资源较多时对系统造成的经济损失较为相近。

4.3 电–气耦合系统负荷节点脆弱性分析

在信息攻击中，攻击者篡改节点负荷的量测数据，达到破环电–气耦合系统运行的目的。定义指标Lvul表征负荷节点遭受LR攻击的程度，反映节点在LR攻击下的脆弱性，如式（37）所示：

图5（a）、（b）和（c）分别给出了不同攻击下各负荷节点的脆弱程度。

由图5（a）可知，在攻击1中，随着攻击资源的不断提升，部分负荷节点的脆弱性变化明显。如节点14、18的脆弱性在攻击资源为 0.5Rc和 0 .8Rc时较高，但在0.6Rc和 0.7Rc时维持在较低水平。表明了在攻击资源总量和一定区域内线路传输容量限制的情况下，不同区域内攻击资源的分配存在一定的波动性。当攻击资源较少时，对于节点18的攻击更容易使系统引发连锁故障，达到切机切负荷、扩大系统经济损失的目的；随着攻击资源的提升，将攻击资源分配给其他区域的负荷节点后对系统造成的损失明显增大，而在攻击资源总量的限制下，分配给节点18的攻击资源会大大减少，因此造成了系统高脆弱性节点的转移；当攻击资源达到 0 .8Rc时，对其他区域内负荷分配的攻击资源趋于饱和，因此又将剩余的资源重新分配给节点18所在区域的负荷节点，导致节点18重新成为高脆弱性节点。为了使得系统总体经济损失达到最高，攻击资源会优先向满载发电机组的邻接负荷和大负荷节点分配，诱使潮流发生大规模变化和发电机出力转移。如负荷占比较大的节点5和邻接多个发电机组的节点1、12的脆弱性一直保持在较高水平，说明了部分负荷节点会成为攻击的首要目标。

由图5（b）可知，在攻击2中，虽然受到不同攻击资源的影响，系统节点脆弱性并未发生较大改变。如节点19、20、21一直保持较高的脆弱性，表明了攻击若仅发生在气网侧，攻击的主要目标是最大化削减气转电负荷。节点20、21与最大的气转电负荷相邻，给节点20和21分配较多的攻击资源会大幅提高气转电负荷削减量。节点19与负载量最大的气源直接相连，针对节点19的攻击可以使整个气网的气源出力产生有效转移，达到最大化削减气网负荷的效果。

由图5（c）可知，攻击3中系统的高脆弱性节点与攻击1、2相比产生了较大的变化，电网与气网中高脆弱性节点数目的变化体现了攻击资源在电网和气网协同分配的过程。随着攻击资源的增加，在攻击3中系统脆弱性节点转移原理与攻击1、攻击2相同，在此不再赘述。为清晰刻画由于攻击方式的转变而导致系统脆弱性发生的变化，对式（37）的节点脆弱性指标Lvul进行修正，如式（38）所示，修正后的Avul指标可以更明确地表征同一负荷节点在不同攻击资源下的脆弱程度，反映了负荷节点的综合脆弱性。式中：Avul为负荷节点的综合脆弱性；n为攻击资源场景个数，本文中n取值为4， αi为攻击资源i可能出现的概率，满足正态分布规律，其详细参数见表2。

表3给出了采用式（38）计算不同攻击方式下系统各节点的综合脆弱性。

由表3可知，在攻击资源有限的前提下，双侧协同LR攻击中的资源分配给气网侧的比重较高。如节点20的脆弱性达到0.430，高于其他负荷节点，反映了在攻击3中，电网侧的LR攻击主要作用是诱使燃气机组出力提升，真正使系统遭受严重经济损失是依靠LR攻击在气网侧削减气转电负荷，降低燃气机组出力，从而在电网侧引起连锁故障实现的。在攻击1和2中系统高脆弱性节点依次分别为12、1、10、19、20，而攻击3中系统高脆弱性节点排序分别为18、5、17、20、21。可见电网中的高脆弱性节点有较明显的转移，说明不同攻击形式对电网的影响有较大差异。同时，气网中节点的脆弱性有一定的变动，表明了虽然不同攻击在气网中的目标一致，但气转电负荷的变动会改变气网中攻击资源的分配。

以上结果表明，当LR攻击资源改变后，系统中的部分高脆弱性节点会发生转移，但仍然有一定比例的负荷节点的脆弱性相对稳定。同时，若改变LR攻击形式，系统的脆弱性节点也会发生较大的转移。所定义的节点综合脆弱性指标有效地反映了不同攻击场景下各节点遭受攻击的程度和高脆弱性节点转移规律。这可以为IEGS防护策略提供相应思路。

4.4 提升系统抵御攻击能力的策略分析

由于LR攻击会对系统安全运行造成严重影响，为提升系统抵御攻击的能力，提出以下3种针对LR攻击的系统保护策略。策略1，加强保护攻击1、2下系统的高脆弱节点；策略2，加强保护攻击3下系统的高脆弱性节点；策略3，限制气转电负荷的削减比例及加强3种攻击下共同的高脆弱性节点。

采用在负荷节点部署冗余量测装置和增强负荷节点量测冗余度的方式保护脆弱性节点，提升负荷节点被攻击的难度。根据表3中负荷节点在不同攻击形式下的综合脆弱性可知：在策略1中选取节点12、1、19作为保护节点；策略2中节点17、18、20为保护节点；策略3中保护节点为5、18、20，并且将气转电负荷削减比例控制在50%以下。

分别采用3种保护策略后系统经济损失的情况如图6所示。

由图6可知，对比加入3种保护策略前后对系统的影响发现，3种保护策略都起到减小系统经济损失，提高系统抵御攻击能力的作用。但对攻击1、2中的高脆弱节点提供保护后，由双侧协同LR攻击造成的系统经济损失并未明显下降，表明保护策略1效果有限，侧面验证了第4.3节中高脆弱性节点转移的正确性。对比3种保护策略，策略3的保护效果显著高于策略1、2，因此，考虑设置气转电负荷削减上限及对共同高脆弱性节点进行保护，可有效减少因LR攻击造成系统严重的经济损失，保障系统安全稳定运行。

5 结 论

本文以电力–天然气信息物理深度耦合系统为研究对象，定量分析了不同LR攻击形式对IEGS造成的经济损失。基于电–气耦合系统的运行特性和双侧协同信息攻击的特点，以负荷重分配的攻击方式最大化系统经济损失，建立损失评估模型，定义新的节点综合脆弱性指标，并提出应对不同攻击的保护策略。仿真结果表明：

1）随着攻击资源的提升，LR攻击对系统造成的损失会进一步加大。但攻击资源提升到一定比例后，攻击对系统造成的经济损失增长缓慢。

2）在攻击资源有限的情况下，合理分配资源对电力系统和天然气系统进行攻击，会使电–气耦合系统遭受最严重的损失；而对脆弱性节点提供保护后，系统抵御攻击能力有所提高。

3）不同的攻击资源和攻击形式都会导致系统的高脆弱性节点发生变化，若只对电网侧LR攻击下的高脆弱性节点进行保护，则双侧协同信息攻击依然会对IEGS造成严重的经济损失。而采用限制气转电负荷的削减比例及加强保护不同攻击下共同高脆弱性节点的策略可以显著提高系统抵御攻击的能力，可为防御者制定合理的防御策略提供重要的依据。

本文在建模时对实际中复杂信息物理系统进行了简化，未考虑以概率模型来表征线路连锁故障，且未计入篡改线路潮流和节点压力所需的攻击资源。因此，下一步将研究信息攻击在电、热（冷）、气等多能源系统的作用方式，采用分布式算法，在更大规模的综合能源系统中建立精细的模型。