(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202210952117.4 (22)申请日 2022.08.09 (71)申请人 贵州电网有限责任公司 地址 550002 贵州省贵阳市南明区滨河路 17号 (72)发明人 曾华荣　杨旗　许逵　文屹　 黄晓旭　陈沛龙　马覃峰　黄欢　 殷蔚翎　 (74)专利代理 机构 南京禹为知识产权代理事务 所(特殊普通 合伙) 32272 专利代理师 褚晓英 (51)Int.Cl. G06Q 10/06(2012.01) G06Q 50/06(2012.01) (54)发明名称 一种基于导纳阵快速三角分解的静态N-1快 速评估方法 (57)摘要 本发明公开了一种基于导纳 阵快速三角分 解的静态N ‑1快速评估 方法， 包括： 读 入电网基础 数据， 形成电纳阵和导纳阵； 选择1条断开的线 路， 进行关联矩阵元素分析， 并基于关联元素对 电纳阵和导纳阵相应的行进行三角分解； 基于修 改后并重新三角分解的电纳阵， 进行二次直流潮 流的过载范围分析； 基于重新三角分解的导纳阵 和初始电流向量进行静态N ‑1评估。 本发明解决 了基于导纳阵进行静态N ‑1计算分析的效率问 题、 并提高准确性， 基于直流潮流电导阵进行分 析范围的界定， 然后基于导纳阵三角分解的非零 元影响范围分析减少三角分析的计算量， 综合提 升静态N‑1计算的效率和准确性， 提高电网覆冰 系统安全分析的效率。 权利要求书2页 说明书6页 附图3页 CN 115375116 A 2022.11.22 CN 115375116 A 1.一种基于导纳阵快速三角分解的静态N ‑1快速评估方法， 其特 征在于， 包括： 读入电网基础数据， 形成电纳阵和导纳阵； 选择1条断开的线路， 进行关联矩阵元素分析， 并基于关联元素对电纳阵和导纳阵相应 的行进行三角分解； 基于修改后并重新 三角分解的电纳阵， 进行二次直 流潮流的过 载范围分析； 基于重新三角分解的导纳阵和初始电流向量进行静态N ‑1评估。 2.如权利要求1所述的基于导纳阵快速三角分解的静态N ‑1快速评估方法， 其特征在 于， 所述电网基础数据包括， 电网潮流数据和电网潮流的计算结果数据。 3.如权利要求2所述的基于导纳阵快速三角分解的静态N ‑1快速评估方法， 其特征在 于， 还包括， 在形成电纳阵和导纳阵后进行初始计算， 所述初始计算包括， 基于电网基础数 据所形成的电纳阵和导纳阵， 完成初次三角分解； 进 行初始直流潮流计算， 得到所有节点的 初始角度 θ0； 根据初始潮流计算结果， 反推出 所有节点的注入电流向量 Iequ0。 4.如权利要求2或3所述的基于导纳阵快速三角分解的静态N ‑1快速评估方法， 其特征 在于， 进行关联矩阵元素分析之前， 对选择的1条断开线路进行修改， 所述修改为删除电纳 阵和导纳阵上断开线路对应的元 素。 5.如权利要求4所述的基于导纳阵快速三角分解的静态N ‑1快速评估方法， 其特征在 于， 所述关联矩阵元 素分析为任选电纳阵和导纳阵其中一个矩阵进行分析， 包括， A1： 建立一个存储节点号的数组， 将断开支路两侧节点 号存入； A2： 提取出所述数组中的第一个节点 号ibus； A3： 在矩阵中寻找 节点号ibus对应行的对角元左侧的第一个非零元 素； A4： 存储所述非零元 素对应的列号到节点 号数组中； A5： 在矩阵中继续 寻找节点号ibus对应行的对角元左侧的下一个非零元 素； A6： 如果存在下一个非零元 素， 则返回A4； 否则， 进行 下一步； A7： 选择节点 号数组中下一个节点 号ibus； A8： 若还有下一个节点 号， 则返回A3， 否则， 结束 整个搜索过程； A9： 保留节点号数组中的所有节点， 按照节点号从小到大的顺序排序并存储， 所述节点 号即为行号。 6.如权利要求5所述的基于导纳阵快速三角分解的静态N ‑1快速评估方法， 其特征在 于， 所述基于关联 元素对电纳阵和导纳阵相应的行进行三角分解， 包括： 读入存储关联节点 号的数组； 选择第一个节点 ibus[i]， i ＝1； 对ibus[i]行进行 标幺化， 即右侧所有非零元 素除对角元； 消除下三角矩阵中第ib us[i]列的所有非零元， 即判断后面所有关联行的ib us[i]列元 素是否是非零元， 若是， 则采用ibus[i]行元素消去； 否则， 跳过。 直到ibus[i]列所有的非零 元变为0为止； 选择下一个节点， 即i＝i+1； 若已经处理完所有关联行， 则计算完成， 退出； 否则返回到 标幺化操作， 继续处 理ibus[i+1]行。 7.如权利要求6所述的基于导纳阵快速三角分解的静态N ‑1快速评估方法， 其特征在 于， 二次直 流潮流的过 载范围分析 得到新的各节点的角度 θ；权　利　要　求　书 1/2 页 2 CN 115375116 A 2根据线路两侧角度差变化估计线路功率变化 率； 存储功率变化 率大于一定限制的线路， 作为检查过 载的线路， 并存 储。 8.如权利要求7所述的基于导纳阵快速三角分解的静态N ‑1快速评估方法， 其特征在 于， 所述功率变化 率表示为： 所述限制表示 为： λ> λlim 其中， θi0、 θj0是计算支路两侧的初始角度， θi、 θj是断开后两侧的角度， λlim为判断变化 率标准。 9.如权利要求8所述的基于导纳阵快速三角分解的静态N ‑1快速评估方法， 其特征在 于， 所述基于重新三角分解的导纳阵和初始电流向量进行静态N ‑1评估， 包括， 计算新的电 压V； 针对直流潮流分析中确定的过载线路， 计算线路电流， 判断是否过载并存储所有过载 的线路。 10.如权利要求9所述的基于导纳阵快速三角分解的静态N ‑1快速评估方法， 其特征在 于， 所述新的电压V， 表示 为： V＝Y‑1Iequ0 其中， Iequ0为注入电流向量， 根据初始状态确定； Y是系统导纳阵。权　利　要　求　书 2/2 页 3 CN 115375116 A 3