【顶级EI复现】基于混合整数线性规划的电力系统蓄意攻击多阶段级联故障脆弱性评估模型（Matlab代码实现）

原创于 2026-06-25 01:01:33 发布 · 47 阅读

本内容遵循CC 4.0 BY-SA版权协议

标签

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

🎁完整资源、论文复现、期刊合作、论文辅导及科研仿真定制事宜点击：

👉👉👉本文完整资源下载

⛳️赠与读者

👨‍💻做科研，涉及到一个深在的思想系统，需要科研者逻辑缜密，踏实认真，但是不能只是努力，很多时候借力比努力更重要，然后还要有仰望星空的创新点和启发点。建议粉丝按目录次序逐一浏览，免得骤然跌入幽暗的迷宫找不到来时的路，它不足为你揭示全部问题的答案，但若能让人胸中升起一朵朵疑云，也未尝不会酿成晚霞斑斓的别一番景致，万一它给你带来了一场精神世界的苦雨，那就借机洗刷一下原来存放在那儿的“躺平”上的尘埃吧。

或许，雨过云收，神驰的天地更清朗.......🔎🔎🔎

💥1 概述

为实现文献所提上层模型的完整复现，本文选取IEEE39节点电力系统作为仿真算例，以系统负荷损失最小化为核心优化目标。该仿真模型主要用于模拟电力调度中心的调控机制，通过优化调整发电机组出力与负荷接入规模，实现对电力系统故障级联扩散过程的有效遏制。在模型构建过程中，本文引入线路运行状态、线路跳闸情况、节点相角等多级级联故障特征变量，完成故障场景的精细化建模。同时，将攻击线路数目、初始故障阶段线路运行状态、线路潮流分布及系统功率平衡条件纳入模型约束体系，保障仿真过程贴合电力系统实际运行规则。本文依托Matlab与Cplex联合仿真求解平台开展数值运算，精准筛选出k重初始故障线路及后续级联扩散引发的故障线路，完整求解故障传播全流程的系统负荷损失数据，所得结果可为该上层模型的有效性验证与性能分析提供核心数据支撑。

摘要

现代电力系统网架结构日趋复杂，输电线路一旦遭受蓄意协同攻击，极易引发潮流大范围转移、支路连锁过载跳闸，进而形成电网孤岛，造成大规模持续性负荷损失，给电网安全稳定运行带来严峻挑战。针对现有静态故障评估模型无法复现级联故障时序演化过程、难以精准刻画孤岛切机切负荷效应的局限性，本文以 IEEE 39 节点新英格兰标准测试系统为研究载体，构建一套面向最坏蓄意攻击场景的多阶段时序级联故障混合整数线性规划评估模型。模型依托直流潮流理论描述电网有功功率传输规律，引入二进制状态变量统一表征初始攻击线路、各阶段支路运行工况、过载保护跳闸行为、发电机组并网状态与负荷供电连通状态；借助大 M 线性化方法将拓扑连通判定、支路过载跳闸等非线性逻辑关系转化为标准线性约束，规避模型非凸性问题，保障求解全局最优。本文将故障演化过程划分为初始源发攻击阶段与多级级联传播阶段，完整还原 “线路蓄意切除 — 全网潮流重分配 — 支路过载跳闸 — 电网孤岛形成 — 机组脱网、负荷切除” 的连锁故障演化全流程。以限定攻击线路数量下系统最终永久负荷损失最大化作为优化目标，求解得到同等攻击强度下危害性最强的关键输电通道组合，同步输出故障演化全过程各阶段连锁停运支路与时序负荷损失数据。所提模型可精准辨识电网高脆弱性输电支路、量化电网抵御蓄意线路攻击的韧性水平，为电网关键线路差异化防护、连锁故障阻断策略制定、极端停电事故应急预案推演提供完备的理论工具与量化分析支撑。

关键词：电力系统；蓄意线路攻击；级联故障；混合整数线性规划；脆弱性评估；电网韧性；孤岛效应

1 引言

1.1 研究背景与意义

电力系统作为国家关键基础设施的核心组成部分，承担着工业生产、居民生活、公共服务等多领域能源供给职能，其安全稳定运行直接关系社会经济平稳运转与公共安全。随着新型电力系统建设持续推进，电网网架规模不断扩张，交直流混联、新能源高占比、多区域互联等特征日益凸显，网络耦合关系愈发紧密，系统运行不确定性显著提升。与此同时，网络空间与物理电网深度融合，针对输电网络的人为蓄意协同攻击已成为电网安全威胁的重要来源。相较于随机自然故障，蓄意攻击具备明确目的性，攻击者会优先选择多条关键输电线路实施同步切除，造成电网拓扑结构性破坏；潮流在剩余完好支路上重新分配后，极易出现多条线路传输功率突破额定容量的现象，触发继电保护连锁动作，故障沿网架持续扩散，最终形成大面积孤岛与持续性停电事故，造成巨大经济损失与社会影响。

传统电网安全分析手段多采用单故障、N-1 静态安全校核模式，仅能评估单一线路停运后的系统运行状态，无法模拟多条线路协同攻击引发的多级连锁跳闸过程，同时忽略拓扑割裂后机组、负荷孤岛停运带来的切机切负荷损失，评估结果偏于乐观，难以反映极端蓄意攻击下系统真实失效程度。现有级联故障仿真方法多采用时序迭代仿真思路，依靠潮流计算逐次模拟故障传播，仅能得到单次随机故障场景，无法主动搜寻对系统危害最大的最坏攻击组合，存在场景遍历不完备、最优脆弱通道难以识别等缺陷。在此背景下，构建能够主动搜寻极端攻击场景、完整复现级联故障时序演化、精准量化孤岛停电损失的优化评估模型，对电网脆弱性分析、韧性提升与主动防御体系构建具备重要理论价值与工程实用意义。

1.2 国内外研究现状

当前国内外围绕电网蓄意攻击与级联故障开展了大量研究，现有成果可分为故障仿真模拟、优化模型脆弱性分析两大方向。

在故障时序仿真领域，主流研究基于直流潮流或交流潮流搭建逐阶段迭代仿真框架，依次模拟线路停运、潮流重分配、过载跳闸全过程，直观复现连锁故障传播路径。该类方法逻辑简单、贴合电网实际故障处置流程，但存在明显短板：仿真依赖随机故障触发，仅能复现单一故障演化场景，无法主动搜寻危害性最大的协同攻击组合；若需遍历全部线路攻击组合，计算量随网架规模呈指数级增长，大规模电网下不具备工程可行性。部分研究引入图论、复杂网络指标辨识电网薄弱支路，依靠介数、度中心性等拓扑指标筛选关键线路，但此类方法仅考虑网络几何结构，忽略电网有功功率传输特性与支路容量约束，筛选得到的脆弱线路与实际潮流过载诱发的连锁故障场景存在显著偏差。

在优化建模脆弱性分析方向，众多学者采用混合整数规划方法构建最坏攻击模型，通过优化算法主动求解最大化系统损失的攻击线路组合。早期相关研究多构建单阶段静态故障模型，仅评估初始线路切除后的系统运行状态，未考虑过载引发的次级、三级连锁跳闸，低估极端攻击下系统停电规模；部分多阶段优化模型简化孤岛效应判定逻辑，直接设定固定比例负荷切除，未依据网架连通关系区分机组与负荷供电状态，损失量化精度不足。此外，部分模型引入交流潮流约束，模型内部存在非线性等式，求解难度大幅提升，易陷入局部最优，难以保证最坏攻击场景的全局最优解。

综合现有研究不足，本文针对性开展改进优化：一是构建多阶段时序统一优化框架，将初始攻击、多级连锁跳闸纳入同一套线性约束体系，无需迭代仿真即可一次性求解完整故障演化路径；二是精细化刻画孤岛连通逻辑，依据支路完好状态判定机组、负荷是否与主网连通，精准量化拓扑割裂带来的切机、切负荷规模；三是基于直流潮流搭建全线性约束体系，依托大 M 法完成全部逻辑关系线性转换，模型为标准混合整数线性规划，可通过商业求解器快速获取全局最优攻击场景。

1.3 论文主要工作与章节安排

本文主要研究工作如下：

完成电力系统元件、故障阶段、网络拓扑关联关系的标准化定义，搭建适配 MATPOWER 标准算例的多阶段级联故障建模基础框架；
引入多类型二进制状态变量区分攻击行为、支路运行、跳闸保护、机组并网、负荷供电五类工况，建立分层时序约束体系刻画故障动态演化过程；
结合直流潮流理论构建全网分阶段功率平衡约束，依托连通性判定逻辑精细化表征孤岛切机、切负荷物理过程；
基于大 M 线性化方法转化支路过载跳闸、拓扑连通判定等耦合逻辑，消除模型非线性项，形成标准 MILP 评估模型；
以全网永久负荷损失最大化为优化目标，实现同等攻击规模下最坏协同攻击线路组合的全局寻优；
基于 IEEE 39 节点新英格兰系统开展模型框架验证，明确模型输出指标与工程应用场景，总结模型创新优势与适用范围。

论文章节安排如下：第一部分为引言，阐述研究背景、国内外研究现状与本文核心工作；第二部分介绍模型基础定义与系统参数划分，明确故障阶段分层机制；第三部分详细阐述模型变量设计思路，区分二进制逻辑变量与连续物理变量的物理含义；第四部分系统梳理模型全部约束条件，逐条阐释约束对应的电网物理机理与时序演化逻辑；第五部分说明优化目标函数构建思路，解释最大化负荷损失的工程意义；第六部分整体剖析模型架构、级联故障完整演化机理，阐述模型在电网脆弱性评估、韧性测算、防御策略优化等场景的应用价值；第七部分介绍模型求解工具、求解流程与量化评估输出指标；第八部分归纳模型核心创新点与相较于现有方法的优势；最后为全文总结与后续研究展望。

2 模型基础定义与故障阶段划分

2.1 电力系统元件集合基础定义

本文模型依托 MATPOWER 标准电力系统算例搭建，统一划分四类核心电网元件集合，分别为全网母线节点、输电支路、同步发电机组、有功负荷节点，同时区分平衡松弛节点与普通非平衡节点。平衡节点作为全网功率兜底单元，承担故障演化全过程中所有功率盈余与缺额，其余普通节点遵循基尔霍夫功率守恒定律，支路、机组、负荷均依托母线节点完成拓扑耦合。

为量化描述元件之间的拓扑连接关系，构建三类标准化关联矩阵，分别表征支路与母线、发电机组与母线、负荷与母线的接入对应关系。关联矩阵仅通过 0、1 元素标记元件接入位置，可快速实现全网功率注入、支路潮流的矩阵化运算，简化大规模电网功率平衡约束的书写与求解流程。针对输电支路电气参数，仅保留电抗参数搭建直流潮流运算体系，忽略电压幅值、无功功率变化，在保证有功潮流计算精度的前提下，大幅降低模型变量维度与求解复杂度。

模型预先设定固定基础常数，包含蓄意攻击线路总数量、支路潮流额定传输上限、大 M 线性化常数、极小松弛阈值四类参数。其中大 M 常数取极大数值，用于实现 0-1 二进制变量与连续潮流变量耦合逻辑的线性转换；极小松弛阈值用于消除二进制约束求解过程中的数值误差，保证连通性判定逻辑稳定收敛。

2.2 级联故障时序阶段分层机制

为完整复现蓄意攻击诱发连锁故障的动态传播过程，本文将故障演化全过程划分为两层时序阶段，分别为初始源发攻击阶段与多级级联传播阶段，全部阶段纳入同一套优化模型同步求解，无需分步迭代计算。

第一层为初始源发攻击阶段，作为故障演化的起点。该阶段模拟攻击者同步切除指定数量输电线路的攻击行为，电网拓扑发生结构性突变，未遭受攻击的线路维持完好运行状态，被攻击线路直接退出输电服务。拓扑改变后全网有功潮流重新分配，基于支路完好状态判定机组、负荷孤岛工况，完成首轮切机、切负荷规模计算；同时对比各支路实际传输功率与额定容量，识别本阶段潮流越限、即将触发保护跳闸的输电支路，为下一阶段拓扑更新提供依据。

第二层为多级级联传播阶段，覆盖故障持续扩散的全过程，可根据研究需求自定义设置传播阶段总数。每一级传播阶段均遵循统一演化逻辑：基于上一阶段最终支路运行状态更新当前电网拓扑，永久移除上一阶段识别出的过载跳闸线路；重新求解全网直流潮流，再次判定孤岛机组与负荷，更新全网负荷损失总量；遍历全部完好支路识别新增过载跳闸线路，标记为下一阶段停运支路。支路一旦在任意阶段触发跳闸，后续所有故障阶段均保持永久停运状态，无法恢复供电，以此模拟输电线路永久性故障保护动作逻辑。当完成全部预设传播阶段计算后，取最后一个稳定阶段的负荷损失总量作为评估蓄意攻击危害程度的核心指标。

3 模型优化变量设计与物理含义

本文模型变量按照功能属性划分为二进制逻辑状态变量与连续物理变量两大类，全部变量均绑定故障阶段时序属性，能够动态反映每一级故障演化过程中电网拓扑、潮流、功率的实时运行状态。

3.1 二进制 0-1 逻辑状态变量

二进制变量是实现故障时序、拓扑连通、保护跳闸逻辑建模的核心载体，全部变量取值仅为 0 或 1，分别对应两种互斥的电网运行工况，共设计五类二进制变量。

第一类为初始攻击线路标识变量，仅用于区分线路是否遭受蓄意协同攻击。变量取值为 0 时代表该线路被攻击者选定切除，在初始阶段直接退出运行；取值为 1 代表线路未遭受初始攻击，初始拓扑中保持完好投运。通过对该类变量求和约束，精准控制攻击者单次能够切断的线路总条数，实现攻击强度的定量限定。

第二类为分阶段支路运行状态变量，对应每一条输电线路在每一个故障阶段的投运工况。变量取值为 1 代表当前阶段线路完好，参与全网有功功率传输；取值为 0 代表线路因初始攻击或前期过载跳闸永久停运，无功率流通。该变量具备时序递推属性，串联起全部故障阶段的电网拓扑变化，是复现连锁故障拓扑演化的核心变量。

第三类为分阶段支路过载跳闸标识变量，用于标记单阶段内潮流越限、触发继电保护动作的输电线路。变量取值为 1 代表当前阶段线路传输功率突破额定容量，执行永久跳闸操作；取值为 0 代表支路潮流处于安全区间，持续正常运行。该变量与支路运行状态变量联动，实现故障阶段之间拓扑更新的递推逻辑。

第四类为分阶段发电机组并网连通状态变量，用于判定机组是否通过完好输电支路与主网形成连通。若机组全部接入联络支路均处于停运状态，该变量取值为 0，判定机组形成独立孤岛，机组有功出力直接归零；若存在至少一条完好联络支路，变量取值为 1，机组保持并网运行，输出额定初始有功功率。该变量精准刻画拓扑割裂带来的机组脱网效应，避免模型高估系统供电能力。

第五类为分阶段负荷供电连通状态变量，逻辑与机组连通变量保持一致，用于判定负荷节点是否具备完整供电通道。负荷全部接入支路停运时变量取 0，负荷完全切除，供电功率归零；存在完好联络支路时变量取 1，负荷全额正常供电。依托该变量可精准量化每一故障阶段全网切负荷规模，为损失评估提供直接计算依据。

3.2 连续物理变量

连续变量用于表征电网运行过程中的模拟物理量，无取值离散限制，包含节点相角、支路有功潮流、机组实时出力、负荷实际供电功率、单负荷永久停电损失五类变量。

节点相角变量对应每一个故障阶段全部母线的电压相角，是直流潮流方程的核心中间变量，通过相角差值直接计算各支路有功传输功率，搭建节点电压状态与支路潮流之间的线性映射关系。支路有功潮流变量记录单阶段每条输电线路的实际传输有功功率，结合支路额定容量判定过载跳闸行为，同时满足全网功率守恒约束。

机组实时出力变量由机组并网二进制状态变量与初始额定出力相乘得到，机组脱网时出力归零，并网时维持基准出力，无需额外设置机组出力上下限约束，简化模型复杂度。负荷实际供电功率变量同理，由负荷连通状态与基准负荷数值耦合计算，直观反映单负荷节点当前供电水平。永久负荷损失变量为初始基准负荷与最终稳定阶段实际供电负荷的差值，代表蓄意攻击引发级联故障后，该负荷节点永久性停电规模，全部负荷损失求和即为系统整体失效程度量化指标。

4 模型约束体系与电网物理机理阐释

本文模型全部约束均为线性等式或线性不等式，无任何非线性项，依托分层设计思路分为攻击规模约束、拓扑时序递推约束、直流潮流约束、机组与负荷连通判定约束、全网功率平衡约束、支路过载跳闸逻辑约束六大模块，各模块约束相互耦合，完整覆盖蓄意攻击、潮流传输、孤岛形成、连锁跳闸全流程物理规律。

4.1 蓄意攻击规模定量约束

该组约束针对初始攻击行为进行定量管控，通过对标记攻击线路的二进制变量求和，限定攻击者同步切除的输电线路总数量。约束保证模型寻优过程中，攻击线路条数严格等于预设数值，避免出现攻击强度偏离设定值的求解结果，实现不同攻击强度场景下电网脆弱性的对比分析。通过调整约束内攻击线路数量参数，可分别模拟少量关键线路协同攻击、大规模多线路同步攻击等不同极端场景，适配多样化电网安全评估需求。

4.2 多阶段拓扑时序递推约束

拓扑递推约束分为初始阶段拓扑约束与级联传播阶段拓扑更新约束两部分，搭建起故障时序演化的逻辑骨架。初始阶段约束直接将支路运行状态与初始攻击标识变量绑定，被攻击线路第一阶段直接停运，完好线路维持投运状态，确定故障演化起点的基础网架结构。

级联传播阶段拓扑更新约束体现故障因果传递关系，当前阶段支路完好状态等于上一阶段运行状态减去本阶段新增跳闸线路。物理层面代表线路一旦因过载触发跳闸，后续所有故障阶段永久退出运行，网架仅能持续简化、无法恢复，完全贴合输电线路永久性故障的保护动作逻辑。依靠该组约束，模型无需分步迭代，仅通过一组线性等式即可统一表征多阶段拓扑连续变化过程。

4.3 分阶段直流潮流线性约束

直流潮流等式约束建立节点相角与支路有功潮流的线性对应关系，适用于全部故障阶段，依托支路电抗参数构建标准化映射矩阵，矩阵运算替代逐支路潮流计算公式，大幅简化大规模电网约束书写流程。

配套增设大 M 不等式约束实现停运支路潮流强制归零：当线路处于停运状态时，不等式约束限制支路潮流取值区间，结合潮流等式唯一推导出潮流为零，模拟停运支路无功率传输的物理现象；当线路完好投运时，不等式约束无实质限制，支路潮流可随节点相角自由变化，完整还原正常输电工况。整套潮流约束无平方、绝对值等非线性算子，全程保持线性特征，保障模型求解稳定性。

4.4 机组、负荷孤岛连通判定约束

连通判定约束是精细化刻画孤岛效应的核心线性不等式组，分别针对发电机组、负荷节点设计两套结构一致的上下界约束，依托支路运行状态求和结果判定元件是否与主网连通。约束引入极小松弛阈值消除数值求解误差，分两种工况实现二进制状态变量的唯一确定。

当机组或负荷节点全部接入联络支路均处于停运状态时，支路完好状态求和结果为零，上下界不等式共同强制并网、供电二进制变量取值为零，判定元件孤岛脱离主网；当存在任意一条完好联络支路时，求和结果大于零，约束强制二进制变量取值为 1，元件维持正常并网、供电状态。相较于传统模型固定比例切负荷的简化处理方式，该组约束严格依托网架拓扑连通关系判定停电规模，评估结果更贴合电网实际故障运行特征。

4.5 全网分阶段有功功率平衡约束

功率平衡约束完整复现基尔霍夫功率守恒定律，分层实现机组负荷功率映射、非平衡节点净注入功率计算、全网潮流守恒、平衡节点功率兜底四项物理规律。首先通过二进制连通状态变量耦合初始基准功率，得到当前阶段机组实际出力与负荷实际供电功率；其次通过关联矩阵矩阵运算，快速计算全部非平衡母线的净有功注入功率，即并网机组总出力减去正常供电负荷总量。

针对普通非平衡母线，构建潮流守恒等式，要求母线净注入功率等于所有相连支路的净潮流代数和；针对全网唯一平衡松弛节点，构建兜底功率平衡约束，平衡节点承担全部非平衡节点的功率盈余或缺额，消除全网功率不平衡量，匹配直流潮流标准计算规则。整套功率平衡约束采用矩阵化表达，逻辑清晰、扩展性强，可直接适配任意节点规模的 MATPOWER 标准算例。

4.6 支路过载跳闸大 M 线性化约束

支路过载跳闸约束为两组耦合大 M 不等式，用于将 “潮流绝对值超过额定容量则线路跳闸” 的非线性逻辑转化为线性约束，适配 MILP 求解框架。约束同时考虑支路运行状态、实际传输功率、额定容量、跳闸标识四类变量，分三种工况实现逻辑自洽判定。

当线路完好投运且潮流超出安全上限时，不等式组强制跳闸标识变量取 1，标记线路在下一阶段永久停运；当线路完好投运、潮流处于安全区间时，约束强制跳闸标识变量取 0，支路持续正常输电；当线路已处于停运状态时，支路潮流天然归零，不会触发过载跳闸，跳闸标识自动置零。整套约束无需引入分段函数、绝对值算子等非线性结构，仅依靠线性不等式即可完整模拟输电线路过载保护动作逻辑。

5 优化目标函数构建与工程内涵

本文模型定位为最坏蓄意攻击场景搜寻模型，优化目标为在限定攻击线路数量的前提下，最大化多阶段级联故障演化结束后全网永久性负荷损失总量。负荷损失定义为正常工况下基准负荷与故障全部传播阶段结束后稳定工况实际供电负荷的差值，代表攻击引发连锁故障后无法恢复的持续性停电规模，是衡量电网遭受蓄意攻击失效程度最直观、最贴合工程实际的量化指标。

最大化负荷损失的优化导向对应电网安全分析中的极限最坏场景：模型通过全局寻优遍历全部符合攻击数量约束的线路组合，筛选出能够诱发最大规模停电事故的协同攻击方案，该组合即为网架内危害最高的脆弱输电通道集合。相较于随机选取攻击线路、依靠仿真得到故障场景的传统方法，本文优化目标可主动定位电网防护重点，避免遗漏高风险线路组合。

现有商业整数规划求解器仅支持极小化标准形式优化问题，因此在输入求解器时对原始最大化目标函数取负值转换为极小化形式；模型求解收敛后，将求解得到的目标函数数值取相反数，即可得到最优攻击场景下全网最大永久负荷损失总量，用于后续电网脆弱性分析与韧性量化评估。

6 模型整体架构、故障演化机理与工程应用价值

6.1 模型整体分层架构

本文所提多阶段级联故障 MILP 模型分为三层递进架构，三层架构变量、约束相互耦合，统一纳入单次优化求解流程，不存在分步迭代计算，计算效率与全局最优性显著优于传统时序仿真模型。

第一层为上层攻击决策层，以初始攻击线路二进制变量为核心，依托攻击规模约束管控攻击强度，通过目标函数寻优实现最坏攻击组合全局搜索，完成电网脆弱线路识别的核心功能；第二层为中层时序演化层，依托多阶段支路运行、跳闸状态变量与拓扑递推约束，串联初始攻击与多级连锁跳闸全过程，完整刻画故障随时间扩散的动态过程；第三层为底层电网物理层，包含直流潮流、功率平衡、孤岛连通、过载跳闸全部物理约束，精准还原电网有功传输、拓扑割裂、保护动作的客观运行规律。三层架构全部约束保持线性，模型为标准混合整数线性规划，不存在非凸、非线性求解难点，可直接调用成熟商业求解器快速收敛。

6.2 级联故障完整时序演化机理

模型通过统一优化框架完整复现蓄意攻击诱发连锁故障的因果演化链条，全过程遵循严格时序因果逻辑，无逻辑断层与简化假设。故障演化起点为初始源发攻击阶段，选定预设数量线路同步切除，网架拓扑发生突变，潮流重新分配后判定孤岛机组与负荷，计算首轮停电损失；同步遍历全部完好支路，识别潮流越限、即将跳闸的输电线路。

进入多级级联传播阶段后，以上一阶段支路停运状态为基础更新当前网架，永久移除前期过载跳闸线路；再次求解全网直流潮流，重新判定孤岛工况并更新负荷损失总量；识别当前阶段新增过载支路，标记为下一阶段停运线路，循环完成全部预设传播阶段计算。当所有级联阶段运算完成后，取最终稳定工况下全网负荷损失作为攻击危害评价指标，同步输出最优攻击线路、各阶段连锁跳闸支路、逐阶段负荷损失时序数据，完整还原故障扩散路径与停电规模变化趋势。

6.3 模型多场景工程应用价值

电网高脆弱输电支路辨识：通过求解最大化负荷损失的最优攻击线路组合，定位同等攻击强度下引发最大停电事故的关键输电通道，将该类线路划定为电网核心防护对象，差异化配置在线监测、后备保护、故障限流装置，降低线路遭受攻击后的连锁扩散风险。
电网韧性量化评估：通过调整攻击线路数量、故障传播阶段数两类核心参数，对比不同攻击强度、不同故障扩散时长下系统最大负荷损失，构建电网韧性量化评价指标体系，直观衡量网架抵御蓄意协同攻击的能力，为网架补强、线路增容规划提供数据支撑。
连锁故障阻断防御策略优化：依托模型输出的故障时序传播路径与连锁跳闸支路序列，针对性设计故障阻断方案，例如在高频连锁跳闸支路上配置柔性直流装置、储能快速功率调节设备、可控负荷快速切除装置，从潮流源头抑制支路过载，阻断级联故障持续扩散。
极端停电事故应急预案推演：模型可完整输出极端蓄意攻击场景下各阶段停运线路、孤岛范围、负荷损失时序曲线，电网调度机构可依托该类极限故障数据制定分级故障处置预案，明确不同停电规模下的机组调控、负荷恢复、抢修作业流程，提升大停电事故应急处置能力。
多类型电网算例通用分析：模型框架完全适配 MATPOWER 标准电力系统算例，仅替换母线、支路、机组、负荷原始数据即可应用于 IEEE 14 节点、30 节点、57 节点、118 节点等不同规模测试系统，同时可拓展应用于区域互联实际电网，具备极强通用性与工程迁移能力。

7 模型求解框架与量化评估输出指标

7.1 模型求解工具与流程

本文模型基于 MATLAB 平台搭载 YALMIP 优化建模工具箱完成变量、约束、目标函数标准化搭建，模型构建完成后调用 CPLEX 或 Gurobi 商业整数规划求解器完成全局最优求解。两类求解器针对大规模混合整数线性规划问题具备成熟分支定界求解算法，求解速度快、收敛稳定性强，可精准输出全局最优最坏攻击场景，不存在局部最优解偏差。

完整求解流程分为四步：第一步导入 MATPOWER 标准电网算例，提取母线、支路、机组、负荷全部电气参数，构建拓扑关联矩阵与直流潮流映射矩阵；第二步完成各类二进制逻辑变量、连续物理变量初始化，按照分层约束体系逐条写入攻击规模、拓扑递推、直流潮流、连通判定、功率平衡、过载跳闸约束；第三步构建最大化负荷损失的优化目标函数并转换为求解器标准极小化形式，配置求解器计算精度、输出日志等参数；第四步调用求解器完成模型运算，求解收敛后提取全部变量最优解，整理、计算、可视化各类评估指标。

7.2 核心量化评估输出指标

模型求解收敛后可输出四类核心量化分析指标，支撑电网脆弱性定量评估：

最优协同攻击线路集合：同等攻击数量约束下，造成全网最大永久负荷损失的输电线路编号组合，代表网架内风险最高、防护优先级最强的关键脆弱通道；
全网最大永久负荷损失总量：全部级联故障阶段结束后系统无法恢复的总停电负荷数值，用于定量衡量对应攻击强度下电网失效严重程度；
分阶段连锁跳闸支路序列：每一级级联传播阶段新触发过载跳闸的输电线路清单，清晰展示故障沿网架扩散的完整路径，定位连锁故障传播关键节点；
时序负荷损失变化曲线：记录每一个故障阶段全网累计切负荷总量，通过曲线直观反映负荷损失随故障演化持续增长的变化规律，区分初始攻击阶段与级联传播阶段的停电贡献占比。

8 模型核心创新点与相较于现有方法的优势

完整多阶段时序一体化建模，无需迭代仿真 区别于现有单阶段静态故障模型、分步时序迭代仿真模型，本文将初始攻击、多级连锁跳闸、孤岛演化全过程纳入同一套 MILP 优化框架，一次性求解完整故障时序信息，不存在迭代误差、场景遗漏问题，可主动搜寻全局最坏攻击场景，解决传统仿真方法无法遍历全部高风险攻击组合的缺陷。
全线性化建模，保证全局最优解 模型全程依托直流潮流理论与大 M 线性化方法，将拓扑连通判定、支路过载跳闸、孤岛切机切负荷等全部非线性耦合逻辑转化为线性等式与不等式，模型为标准混合整数线性规划，无任何非凸约束，商业求解器可稳定输出全局最优攻击方案，规避交流潮流优化模型易陷入局部最优的短板。
精细化孤岛效应拓扑判定，负荷损失量化精度更高 摒弃传统模型固定比例切负荷的简化假设，依托支路完好拓扑状态逐台、逐节点判定机组与负荷孤岛工况，严格依据实际供电通道连通状态计算切机、切负荷规模，精准还原拓扑割裂后电网孤岛停电的真实物理过程，负荷损失评估结果更贴合工程实际故障工况。
框架通用性强，适配多尺度电网分析 模型完全兼容 MATPOWER 通用电力系统数据格式，无需大幅修改代码与约束体系，仅替换算例原始数据即可应用于小规模标准测试系统、区域互联实际电网；攻击线路数量、故障传播阶段数可自由调整，灵活适配不同攻击强度、不同故障扩散时长的电网安全评估需求，拓展应用场景丰富。
输出指标体系完备，兼顾机理分析与工程落地 模型不仅可输出最优攻击线路、总负荷损失等宏观评估指标，还能够完整输出故障时序传播路径、逐阶段停电规模等微观演化数据，既能够支撑电网脆弱性机理理论分析，也可直接为电网规划、调度防御、应急预案编制提供量化数据支撑，理论价值与工程实用价值兼备。

9 总结与后续研究展望

9.1 全文总结

针对蓄意线路攻击诱发多级连锁故障、大规模孤岛停电的电网安全问题，本文构建一套基于混合整数线性规划的多阶段时序级联故障脆弱性评估模型。模型依托直流潮流简化电网有功传输计算，引入多类型二进制状态变量完整表征攻击行为、支路运行、过载跳闸、机组并网、负荷供电五类工况；通过大 M 线性化方法实现全部逻辑约束线性转换，搭建分层时序约束体系完整复现 “线路协同攻击 — 潮流重分配 — 连锁过载跳闸 — 电网孤岛切负荷” 故障演化全流程。以限定攻击线路数量下永久负荷损失最大化为优化目标，全局搜寻危害性最强的协同攻击线路组合，同步输出故障时序传播路径与逐阶段停电损失数据。

基于 IEEE 39 节点新英格兰标准测试系统完成模型框架搭建与逻辑验证，所提模型具备全局最优求解、孤岛效应精细化刻画、多阶段时序一体化建模、算例通用性强等多重优势，可精准辨识电网高脆弱输电通道、量化电网抵御蓄意攻击的韧性水平，为电网关键线路差异化防护、连锁故障阻断策略优化、极端大停电应急预案推演提供完善的量化分析工具。

9.2 后续研究展望

本文当前模型基于直流潮流理论搭建，仅考虑有功功率传输与有功负荷损失，未纳入无功电压约束、交流潮流非线性特性，后续可进一步拓展交流潮流线性化建模方法，兼顾电压越限诱发的保护跳闸工况，提升故障仿真完整度；模型当前仅考虑输电线路单一攻击目标，后续可拓展发电机组、变电站、柔性直流换流站等多元关键设备协同攻击场景，完善多目标物理设备蓄意攻击评估框架；现有模型仅针对静态网架开展脆弱性分析，后续可接入储能、可调节负荷、新能源机组等柔性资源调控约束，研究柔性资源对级联故障的阻断效果，构建主动防御优化模型；此外，可结合图论、机器学习算法对大规模电网算例进行脆弱线路预筛选，降低高维 MILP 模型求解维度，进一步提升大规模实际电网的计算效率。

📚2 运行结果

仿真分析结果表明，在初始故障发生后，伴随连锁故障的持续蔓延，电力系统的负荷损失规模将呈现逐步增大的趋势。因此，如何有效管控连锁故障引发的负荷损耗，是本次研究的核心关键。现有相关研究多依托底层建模方式开展调度优化分析，并借助对偶变换方法，将复杂的双层优化模型转化为单层模型求解，有效降低了模型运算难度。该研究思路可为同领域科研人员提供扎实的研究基础，学者可在此之上创新优化，构建全新的负荷损失抑制策略与模型。以此作为研究创新点开展课题研究，能够有效提升论文内容的创新性与完整性，同时助力学术论文的高效、规范撰写。

本程序不仅可以自动识别初始故障所引发的线路停运事件，还能够完整复现后续连锁故障逐层蔓延所造成的多条支路跳闸结果。为保证事故仿真结果的有效性，本节从连通性校验与负荷数值比对两个维度完成结果验证。

首先验证单一初始故障场景，事故造成 37、39、46 三条线路停运。在电力网络模型中将上述支路切除后进行连通性分析，结果显示全网负荷节点始终保持电气互联，电网未形成孤立子系统，不存在任何负荷功率损失，和程序初始故障的输出结论保持一致。

随后研究初始故障与首轮级联故障共同作用的极端工况，本次事故总计切除 8 条故障线路。潮流仿真观测到 21、23、24 节点发生失电孤岛。依托程序生成的 mpc 数据文件，可查询得到这三个节点的负荷数据分别为 274、247.5 和 308.6。对失电负荷进行加总计算后可以发现，总负荷损失量与负荷损失曲线里 s=2 工况下的数值高度吻合。该项定量对比充分证明，程序能够准确刻画连锁故障传播路径，并且可靠计算事故带来的停电负荷规模。

部分代码：

%% 求解结果

if result.problem == 0

% 提取最优解

attack_lines = find(value(ul) <= 0.1);

total_load_loss = value(Objective);

fprintf('需要攻击的线路序号为：%s\n', mat2str(attack_lines'));

fprintf('对应的最终负荷损失总和为：%.2f MW\n', total_load_loss);

% 每个阶段的负荷损失

loss_per_stage = zeros(m,1);

for s = 1:m

loss_per_stage(s) = sum(Dd0) - sum(value(Ds(:,s)));