CN107302208B - 一种煤矿高压电网速断设置并行检验方法 - Google Patents

Abstract

本发明公开了一种煤矿高压电网速断设置并行检验方法，该方法首先以煤矿高压供电系统电源支路节点为起点，将煤矿高压供电系统支路节点分为n个集合，如果不同集合中存在相同的联络开关节点，则将这些集合合并；最后以此为基础基于关联矩阵进行并行网络拓扑分析，依据获得的所有供电关联矩阵Gi生成整个高压供电系统的最终供电关联矩阵FG，基于关联矩阵FG完成煤矿高压电网速断设置并行检验。本发明基于并行计算技术构造煤矿高压电网的拓扑分析调度模型，实现整个煤矿高压供电系统的速断设置并行检验，能够有效提高速断设置检验速度。

Description

一种煤矿高压电网速断设置并行检验方法

技术领域

本发明公开了一种煤矿高压电网速断设置并行检验方法，属于煤矿高压供电网络继电保护整定计算领域。

背景技术

35kV以上的电力网中存在多个电源，属于复杂闭式电网，而矿井高压供电系统为6kV或10kV等级，两个电源应采用分列运行方式，或者是一路使用一路备用，属于单电源开式电网。当前，故障计算和速断设置的仿真检验是针对复杂闭式电网进行的，在复杂闭式电网中，网络结构复杂，需要针对整个供电网络故障情况构造相应的数学模型完成故障计算，再对供电网络速断设置的可靠性和选择进行仿真检验。已有的基于关联矩阵的矿井高压电网速断设置仿真检验方法则主要是针对矿井高压电网(单电源开式电网)的网络结构特点，通过获得支路节点和支路节点的关联矩阵构造矿井高压电网的网络拓扑模型，完成矿井高压电网供电网络的拓扑辨识；并在此基础上依据每个高压开关的瞬时速断设置、捕获到的故障电流，以及支路节点和支路节点的关联矩阵得到最终受该故障影响需要跳闸的所有开关，执行跳闸操作；最终，依据仿真得到的开关跳闸结果判断在此故障情况下，当前的速断设置是否满足可靠性和选择性要求，以达到速断设置仿真检验的目的。

但在矿井高压供电系统中节点数量较多，基于关联矩阵完成整个矿井高压电网速断设置检验，时间开销大，为了能够以较少的时间开销基于关联矩阵完成矿井高压电网的速断设置检验。本发明提出了一种煤矿高压电网速断设置并行检验方法，该方法基于并行计算技术构造煤矿高压电网的拓扑分析调度模型，并以此为基础基于并行计算技术实现整个煤矿高压供电系统的速断设置并行检验。

发明内容

针对某节点集合中存在的节点基于关联矩阵完成网络拓扑分析，其使用的网络拓扑分析方法主要是采用文献“单母线分多段运行的矿井高压电网自适应整定计算方法”中的网络拓扑分析模型，该模型依据支路节点开关状态矩阵、联络节点开关状态矩阵、母线节点和联络节点的关联矩阵、母线节点和支路节点的关联矩阵生成支路节点与支路节点之间的最终供电关系矩阵，构建单母线分多段运行的矿井高压电网网络拓扑模型，为自适应整定计算提供基础；该模型构建网络拓扑分析模型的计算时间开销由该节点集合中存在的支路节点数量、供电系统级数、母线最大分段数、母线节点数量和联络开关节点数量共同决定。

设置煤矿高压供电系统的电源支路节点，电源支路节点是指由上级供电部门直接供电的支路节点。

在某高压供电系统图中，电源支路节点所在的变电所为第1级变电所，由电源支路节点供电的变电所是第2级变电所，由第2级变电所供电的变电所为第3级变电所，由第(N₁-1)级变电所供电的变电所为第N₁级变电所，如果不存在由第N₁级变电所供电的变电所，则该高压供电系统图的供电系统级数为N₁。

在某高压供电系统图中，如果某段母线上的联络开关数量为N₂，则该母线分段数为(N₂+1)，如果该高压供电系统图中所有母线分段数的最大值用N₃表示，则该高压供电系统图的母线最大分段数为N₃。

假定煤矿高压供电系统图中包含Z个支路节点，则依据支路节点上出线开关的开闭状态，生成支路节点开关状态矩阵SK，SK包含Z个元素，SK＝[SK₁ SK₂ ... SK_i ... SK_Z]，1≤i≤Z；在SK中，第i个元素对应的开关状态闭合，则SK_i＝1；反之，则SK_i＝0。

在矿井高压电网中，针对某支路节点直接供电的线路设置故障点，并对该故障点设置相应的故障类型，故障类型为二相短路或三相短路，然后依据该故障点的故障类型计算其对应的最小两相短路电流或最大三相短路电流；针对所有支路节点对应的出线开关设置相应的瞬时速断定值，保存在瞬时速断定值矢量SD中，依据支路节点编号可以在SD中找到相应支路节点对应的出线开关设置的瞬时速断定值。

假定开关状态为合闸的电源支路节点的数量有n个，以煤矿高压供电系统电源支路节点为起点，将煤矿高压供电系统支路节点分为n个集合，每个节点集合用S_i表示，1≤i≤n；具体步骤如下：

步骤1)、将煤矿高压供电系统中所有开关状态为合闸的电源支路节点加入到集合PS中，假定开关状态为合闸的电源支路节点的数量有n个；从PS中取出一个电源支路节点，i的数值设置为1，执行步骤2)；

步骤2)、该电源支路节点用A_i表示，将A_i加入到节点集合S_i中，在煤矿高压供电系统中查找由该电源支路节点A_i直接供电的母线加入对母线集合B中，从母线集合B中取出一个元素，执行步骤3)；

步骤3)取出的母线用B_i表示，将母线B_i直接连接的开关状态为合闸的联络开关节点加入到集合TS和集合S_i中，将母线B_i直接连接的支路节点加入到集合TS和集合S_i中，执行步骤4)；

步骤4)如果母线集合B中仍然有元素未被遍历，则取出下一个元素，执行步骤3)；如果母线集合B为空，执行步骤5)；

步骤5)针对集合TS中存在的支路节点查找其直接供电的母线，如果查找的母线节点存在，执行步骤6)；如果不存在，执行步骤7)；

步骤6)将集合TS中每一个支路节点直接供电的母线加入集合B中，删除集合TS中的所有元素，从母线集合B中取出一个元素，执行步骤3)；

步骤7)如果集合PS不为空，从集合PS中取出下一个电源支路节点，将i的数值加1，执行步骤2)；如果集合PS为空，则可获得n个支路节点集合，每个节点集合用S_i表示。

针对n个节点集合S_i，其中1≤i≤n；如果不同集合中存在相同的联络开关节点，则将这些集合合并；具体步骤如下：

步骤1)、将n个节点集合S_i加入到集合AS中，UN数值设置为n，从集合AS中取出一个节点集合，执行步骤2)；

步骤2)、取出的节点集合用AS_j表示，集合AS_j与集合AS中保存的每一个支路节点集合AS_k进行比较，如果集合AS_j和集合AS_k中存在相同的联络开关节点，则从集合AS中取出支路节点集合AS_k，将集合AS_j和集合AS_k合并后，加入到集合NS中；如果集合AS_j与集合AS中保存的每一个支路节点集合AS_k不存在相同的联络开关节点，则将集合AS_j加入到集合NS中；

步骤3)、如果集合AS不为空，从集合AS中取出一个节点集合，执行步骤2)；如果集合AS为空，执行步骤4)；

步骤4)、如果集合NS中元素的数量小于UN，则将UN的数值设置为集合NS中元素的数量，将集合NS中的所有元素加入到集合AS中，从集合AS中取出一个支路节点集合，执行步骤2)；如果集合NS中元素的数量等于UN，则将集合NS中所有元素加入到集合AS中，完成了节点集合分组。

构造节点时间权重表PTS，节点时间权重表中包含六个字段：支路节点数量、供电系统级数、母线最大分段数、母线节点数量、联络开关节点数量和时间开销权重；支路节点数量用PTSZN表示，供电系统级数用PTSLN表示，母线最大分段数用PTSSN表示，母线节点数量用PTSBN表示，联络开关节点数量用PTSCN表示，时间开销用PTST表示；依据支路节点数量PTSZN、供电系统级数PTSLN、母线最大分段数PTSSN、母线节点数量PTSBN和联络开关节点数量PTSCN，能够在节点时间权重表PTS中查找到其对应的时间开销。

计算所有集合AS_i完成拓扑分析需要的总的时间开销，具体步骤如下：

步骤1)假定集合AS中包含Q个元素，从集合AS中取出一个节点集合，执行步骤2)；

步骤2)该节点集合表示为AS_i，如果节点集合AS_i中存在重复的节点，则将集合AS_i中重复的节点删除，同一个节点只保留一个；依据集合AS_i中支路节点构成的煤矿高压供电系统图统计其包含的支路节点数量PTSZN_i、供电系统级数PTSLN_i，母线最大分段数PTSSN_i、母线节点数量PTSBN_i和联络开关节点数量PTSCN_i，在时间权重表中查找PTSZN_i、PTSLN_i、、PTSSN_i、PTSBN_i和PTSCN_i对应的时间开销K_i，如果能够找到相应的记录，则节点集合AS_i对应的时间开销等于K_i；如果在时间权重表中找不到相应的记录，则依据集合AS_i中节点构成的供电系统图，基于关联矩阵计算获得支路节点与支路节点之间的最终供电关系矩阵所需要的时间K_i，在时间权重表中增加一条记录，该记录支路节点数量PTSZN等于PTSZN_i，供电系统级数PTSLN等于PTSLN_i，母线最大分段数PTSSN等于PTSSN_i，母线节点数量PTSBN等于PTSBN_i，联络开关节点数量PTSCN等于PTSCN_i，时间开销字段等于K_i；执行步骤3)；

步骤3)如果集合AS不为空，从集合AS中取出一个节点集合，执行步骤2)；如何集合AS为空，则所有集合AS_i完成拓扑分析需要的总的时间开销

获取系统当前活动线程数A1及系统允许建立的最大线程设置数B1，则允许建立线程数W＝B1－A1，基于负载均衡技术将Q个节点集合AS_i分别加入到h个队列Q_j中，其中1≤j≤h，h≤W；具体执行步骤如下：

步骤1)、获取系统当前活动线程数A1及系统允许建立的最大线程设置数B1，则允许建立线程数W＝B1－A1；

步骤2)、每个线程平均时间开销按照每一个节点集合AS_i对应的时间开销K_i(1≤i≤Q)以从大到小的顺序对所有节点集合重新进行排序，按照从大到小的顺序将排序后的节点集合依次加入到队列QS中，j的初始值设置为1，从队列QS中按照先进先出的原则取出一个节点集合，将该节点集合加入到队列Q_j中；执行步骤3)；

步骤3)、从队列QS中按照先进先出的原则取出一个节点集合，该集合用G表示，该节点集合对应的时间开销用D表示，假定在队列Q_j中保存的所有节点集合对应的总的时间开销用F_j表示，如果(F_j+D)>AVG，则将节点集合G加入到队列QT中，执行步骤4)；如果(F_j+D)≤AVG，将节点集合G加入到队列Q_j中，执行步骤4)；

步骤4)、如果队列QS中仍然有节点集合未被取出，执行步骤3)；如果队列QS中所有节点集合都已被取出，执行步骤5)；

步骤5)、如果队列QT不为空，则将j的数值加1，将队列QT中的所有节点集合按照先进先出的顺序加入到队列QS中，从队列QS中按照先进先出的原则取出一个节点集合，将该节点集合加入到队列Q_j中，执行步骤3)；如果队列QT为空，则将j的值赋予h，执行步骤6)；

步骤6)、如果h<m，将h的值赋予m，删除队列Q₁，…，Q_h中包含的所有节点集合，1≤j≤h，执行步骤2)；如果h和m相等，则Q个节点集合AS_i被分别加入到h个队列Q_j中，1≤j≤h。

创建h个新的空闲线程，将h个新的空闲线程加入到空闲线程队列QC中，针对每个队列Q_j中保存的每一个节点集合，基于关联矩阵进行并行网络拓扑分析，获得每个节点集合AS_i对应的最终供电关联矩阵G_i(1≤i≤h)；矩阵G_i中的第a行第b列的元素用G_iab表示，如果某支路节点a由支路节点b供电，则G_iab＝1，反之则G_iab＝0；在完成矩阵运算时，将矩阵元素与矩阵元素乘法运算定义为二进制与运算，将矩阵元素与矩阵元素的加法运算定义为二进制或运算；具体步骤如下：

步骤1)、创建h个新的空闲线程，将h个新的空闲线程加入到空闲线程队列QC中，j的初始值设置为1；

步骤2)、从空闲线程队列QC中取出一个空闲线程，将队列Q_j绑定在新建的空闲线程中，将此线程设置为繁忙线程，加入到繁忙线程队列B2中，执行步骤3)；

步骤3)、如果j<h，将j的数值加1，执行步骤2)；如果j≥h，执行步骤4)；

步骤4)、针对繁忙线程队列B2中的每一个繁忙线程，执行步骤5)；

步骤5)、取出所述绑定队列Q_j的繁忙线程，执行该繁忙线程，该线程针对队列Q_j中包含的每一个节点集合执行步骤6)；

步骤6)、该节点集合用Ps表示，对Ps中存在的节点基于关联矩阵进行并行网络拓扑分析，获得该节点集合Ps对应的最终供电关联矩阵G_i(1≤i≤h)；线程执行完毕后，释放该繁忙线程；

在步骤6)中，基于关联矩阵进行并行网络拓扑分析的方法主要是采用文献“单母线分多段运行的矿井高压电网自适应整定计算方法”中的网络拓扑分析模型，该模型依据支路节点开关状态矩阵、联络节点开关状态矩阵、母线节点和联络节点的关联矩阵、母线节点和支路节点的关联矩阵生成支路节点与支路节点之间的最终供电关系矩阵G_i。

依据获得的所有供电关联矩阵G_i(1≤i≤Q)生成整个高压供电系统的最终供电关联矩阵FG；在完成矩阵运算时，将矩阵元素与矩阵元素乘法运算定义为二进制与运算，将矩阵元素与矩阵元素的加法运算定义为二进制或运算；具体执行步骤如下：

步骤1)、依据关联矩阵G_i(1≤i≤Q)生成其对应的左侧变换矩阵LG_i，假定关联矩阵G_i为m_i行m_i列，变换矩阵LG_i为Z行m_i列；将变换矩阵LG_i中所有元素的数值设置为0；针对每个变换矩阵LG_i(1≤i≤Q)执行步骤2)；

步骤2)、如果i＝1，则将变换矩阵LG_i中第q行第q列元素的数值设置为1，其中1≤q≤m_i；如果i>1，则将变换矩阵LG_i中第行第q列元素的数值设置为1，其中₁≤q≤m_i；

步骤3)、依据关联矩阵G_i(1≤i≤Q)生成其对应的右侧变换矩阵RG_i，变换矩阵RG_i为m_i行Z列；将变换矩阵RG_i中所有元素的数值设置为0；针对每个变换矩阵RG_i执行步骤4)；

步骤4)、如果i＝1，则将变换矩阵RG_i中第q行第q列元素的数值设置为1，其中1≤q≤m_i；如果i>1，则将变换矩阵RG_i中第q行第列元素的数值设置为1，其中1≤q≤m_i；

步骤5)、依据供电关联矩阵G_i、左侧变换矩阵LG_i和右侧变换矩阵RG_i(1≤i≤Q)生成整个高压供电系统的最终供电关联矩阵FG；最终供电关联矩阵

依据支路节点与支路节点供电关联矩阵FG获得瞬时速断定值小于故障电流的支路节点集合TT，具体步骤如下：

步骤1)、假定所设置的某线路故障点是由支路节点k直接供电的，依据支路节点k在供电关联矩阵FG中对应的行号查找到相应的行，然后找到该行中数值为1的所有元素对应的列号，再依据获得的所有列号找到对应的支路节点集合SQ；SQ是受短路故障影响的支路节点序列，在支路节点序列SQ中所有元素的故障电流均为所设置故障点当前的短路电流；

步骤2)、针对支路节点序列SQ中的每一个支路节点，在瞬时速断定值矢量SD中找到其相应的速断定值，然后将支路节点序列SQ中每个元素的故障电流和该元素所设置的瞬时速断定值进行比较，如果故障电流大于瞬时速断定值，则将该元素加入到支路节点集合TT中，TT表示获得的需要跳闸的部分支路节点开关集合。

当TT中某支路节点开关跳闸时，为了能够准确地反映现场环境，由该开关供电的所有下级开关也要全部跳闸；因此需要依据新的开关状态矢量NSK和TT获得因短路故障需要跳闸的所有开关集合WW，依据获得的集合WW将相应开关的状态修改为分闸，执行跳闸操作；具体步骤如下：

步骤1)、针对TT中存在的每一个支路节点，在开关状态矢量SK中将对其相应的开关状态数值修改为0，获得一个新的开关状态矢量NSK；

步骤2)、基于新的开关状态矢量NSK重新计算支路节点与支路节点供电关系的关联矩阵FG；FG即是TT中的支路节点跳闸后，支路节点与支路节点的最终供电关联矩阵FF；

步骤3)、针对每一个支路节点重复如下过程：支路节点i直接给线路i供电，依据支路节点i在最终供电关联矩阵FF中对应的行号找到相应的行，然后找到该行中数值为1的所有元素对应的列号，再依据获得的每一个列号找到其对应的支路节点集合P，P是所有给线路i供电的支路节点集合；如果集合P中不存在电源支路节点，则将支路节点i加入到WW中；WW是因短路故障需要跳闸的所有支路节点开关集合；

步骤4)、集合WW中所有元素对应的支路节点，如果在SK中的状态为1，则需要将该支路节点的开关状态修改为分闸，执行跳闸操作；最终，依据仿真得到的开关跳闸结果判断在此故障情况下，当前的速断设置是否满足可靠性和选择性要求，以达到速断设置仿真检验的目的。

附图说明

图1是煤矿高压供电系统图。

具体实施方式

在煤矿高压供电系统图中，如果开关用黑色填充，则表示其开关状态为分闸；反之，开关状态为合闸。

在附图1所示的煤矿高压供电系统图中，开关状态为合闸的电源支路节点的数量有8个，以煤矿高压供电系统电源支路节点为起点，将煤矿高压供电系统支路节点分为n个集合，每个节点集合用S_i表示，1≤i≤n；则在附图1所示的煤矿高压供电系统图中，S₁＝{X1，X9，X10，X25，X26，X27，<4>，X44，X45}，S₂＝{X2，X11，X12，X28，X29，X30，<4>，X46，X47，X48，X49}，S₃＝{X3，X13，X14，<3>，X31，X32，X33，X50，X51}，S₄＝{X4，X15，X16，<3>，X34，X35，X36，X52，X53}，S₅＝{X5，X17，X18，X37，X38，X39}，S₆＝{X6，X19，X20}，S₇＝{X7，X21，X22，X40，X41}，S₈＝{X8，X23，X24，X42，X43}。

在附图1所示的煤矿高压供电系统图中，针对n个节点集合S_i，其中1≤i≤n，n等于8；如果不同集合中存在相同的联络开关节点，则将这些集合合并；则集合AS＝{{X1，X9，X10，X25，X26，X27，<4>，X44，X45，X2，X11，X12，X28，X29，X30，<4>，X46，X47，X48，X49}，{X3，X13，X14，<3>，X31，X32，X33，X50，X51，X4，X15，X16，<3>，X34，X35，X36，X52，X53}，{X5，X17，X18，X37，X38，X39}，{X6，X19，X20}，{X7，X21，X22，X40，X41}，{X8，X23，X24，X42，X43}}。

构造节点时间权重表PTS，节点时间权重表中包含六个字段：支路节点数量、供电系统级数、母线最大分段数、母线节点数量、联络开关节点数量和时间开销权重；支路节点数量用PTSZN表示，供电系统级数用PTSLN表示，母线最大分段数用PTSSN表示，母线节点数量用PTSBN表示，联络开关节点数量用PTSCN表示，时间开销用PTST表示；依据支路节点数量PTSZN、供电系统级数PTSLN、母线最大分段数PTSSN、母线节点数量PTSBN和联络开关节点数量PTSCN，能够在节点时间权重表PTS中查找到其对应的时间开销。

计算所有集合AS_i完成拓扑分析需要的总的时间开销，具体步骤如下：

步骤1)假定集合AS中包含Q个元素，从集合AS中取出一个节点集合，执行步骤2)；

步骤2)该节点集合表示为AS_i，如果节点集合AS_i中存在重复的节点，则将集合AS_i中重复的节点删除，同一个节点只保留一个；依据集合AS_i中支路节点构成的煤矿高压供电系统图统计其包含的支路节点数量PTSZN_i、供电系统级数PTSLN_i，母线最大分段数PTSSN_i、母线节点数量PTSBN_i和联络开关节点数量PTSCN_i，在时间权重表中查找PTSZN_i、PTSLN_i、、PTSSN_i、PTSBN_i和PTSCN_i对应的时间开销K_i，如果能够找到相应的记录，则节点集合AS_i对应的时间开销等于K_i；如果在时间权重表中找不到相应的记录，则依据集合AS_i中节点构成的供电系统图，基于关联矩阵计算获得支路节点与支路节点之间的最终供电关系矩阵所需要的时间K_i，在时间权重表中增加一条记录，该记录支路节点数量PTSZN等于PTSZN_i，供电系统级数PTSLN等于PTSLN_i，母线最大分段数PTSSN等于PTSSN_i，母线节点数量PTSBN等于PTSBN_i，联络开关节点数量PTSCN等于PTSCN_i，时间开销字段等于K_i；执行步骤3)；

步骤3)如果集合AS不为空，从集合AS中取出一个节点集合，执行步骤2)；如何集合AS为空，则所有集合AS_i完成拓扑分析需要的总的时间开销

获取系统当前活动线程数A1及系统允许建立的最大线程设置数B1，则允许建立线程数W＝B1－A1，基于负载均衡技术将Q个节点集合AS_i分别加入到h个队列Q_j中，其中1≤j≤h，h≤w。

创建h个新的空闲线程，将h个新的空闲线程加入到空闲线程队列QC中，针对每个队列Q_j中保存的每一个节点集合，基于关联矩阵进行并行网络拓扑分析，获得每个节点集合AS_i对应的最终供电关联矩阵G_i(1≤i≤h)。

依据获得的所有供电关联矩阵G_i(1≤i≤Q)生成整个高压供电系统的最终供电关联矩阵FG。

依据支路节点与支路节点供电关联矩阵FG获得瞬时速断定值小于故障电流的支路节点集合TT。

当TT中某支路节点开关跳闸时，为了能够准确地反映现场环境，由该开关供电的所有下级开关也要全部跳闸；因此需要依据新的开关状态矢量NSK和TT获得因短路故障需要跳闸的所有开关集合WW，依据获得的集合WW将相应开关的状态修改为分闸，执行跳闸操作。

Claims (1)

1.一种煤矿高压电网速断设置并行检验方法，其特征在于，所描述的速断设置并行检验方法包括如下步骤：

步骤11，假定开关状态为合闸的电源支路节点的数量有n个，以煤矿高压供电系统电源支路节点为起点，将煤矿高压供电系统支路节点分为n个集合，每个节点集合用S_i表示，1≤i≤n；

步骤12，针对n个节点集合S_i，其中1≤i≤n；如果不同集合中存在相同的联络开关节点，则将这些集合合并；

步骤13，构造节点时间权重表PTS，节点时间权重表PTS中包含六个字段：支路节点数量、供电系统级数、母线最大分段数、母线节点数量、联络开关节点数量和时间开销权重；支路节点数量用PTSZN表示，供电系统级数用PTSLN表示，母线最大分段数用PTSSN表示，母线节点数量用PTSBN表示，联络开关节点数量用PTSCN表示，时间开销权重用PTST表示；依据支路节点数量PTSZN、供电系统级数PTSLN、母线最大分段数PTSSN、母线节点数量PTSBN和联络开关节点数量PTSCN，能够在节点时间权重表PTS中查找到其对应的时间开销权重PTST；

步骤14，计算所有节点集合AS_i完成拓扑分析需要的总的时间开销；

步骤15，获取系统当前活动线程数A1及系统允许建立的最大线程设置数B1，则允许建立线程数W＝B1－A1，基于负载均衡技术将Q个节点集合AS_i分别加入到h个队列Q_j中，其中1≤j≤h，h≤W；

步骤16，创建h个新的空闲线程，将h个新的空闲线程加入到空闲线程队列QC中，针对每个队列Q_j中保存的每一个节点集合，基于关联矩阵进行并行网络拓扑分析，获得每个节点集合AS_i对应的最终供电关联矩阵G_i(1≤i≤Q)；

步骤17，依据获得的最终供电关联矩阵G_i(1≤i≤Q)生成整个高压供电系统的支路节点与支路节点供电关联矩阵FG；

步骤18，依据支路节点与支路节点供电关联矩阵FG获得瞬时速断定值小于故障电流的支路节点集合TT；

步骤19，当支路节点集合TT中某支路节点开关跳闸时，为了能够准确地反映现场环境，由该开关供电的所有下级开关也要全部跳闸；因此需要依据新的开关状态矢量NSK和支路节点集合TT获得因短路故障需要跳闸的所有开关集合WW，依据获得的开关集合WW将相应开关的状态修改为分闸，执行跳闸操作；

在步骤11中，主要进行如下步骤；

步骤111、将煤矿高压供电系统中所有开关状态为合闸的电源支路节点加入到集合PS中，假定开关状态为合闸的电源支路节点的数量有n个；从集合PS中取出一个电源支路节点，i的数值设置为1，执行步骤112；

步骤112、该电源支路节点用A_i表示，将该电源支路节点A_i加入到节点集合S_i中，在煤矿高压供电系统中查找由该电源支路节点A_i直接供电的母线加入对母线集合B中，从母线集合B中取出一个元素，执行步骤113；

步骤113、取出的母线用B_i表示，将母线B_i直接连接的开关状态为合闸的联络开关节点加入到集合TS和集合S_i中，将母线B_i直接连接的支路节点加入到集合TS和集合S_i中，执行步骤114；

步骤114、如果母线集合B中仍然有元素未被遍历，则取出下一个元素，执行步骤113；如果母线集合B为空，执行步骤115；

步骤115、针对集合TS中存在的支路节点查找其直接供电的母线，如果查找的母线节点存在，执行步骤116；如果不存在，执行步骤117；

步骤116、将集合TS中每一个支路节点直接供电的母线加入集合B中，删除集合TS中的所有元素，从母线集合B中取出一个元素，执行步骤113；

步骤117、如果集合PS不为空，从集合PS中取出下一个电源支路节点，将i的数值加1，执行步骤112；如果集合PS为空，则可获得n个支路节点集合，每个节点集合用S_i表示；

在步骤12中，主要进行如下步骤：

步骤121、将n个节点集合S_i加入到集合AS中，UN的数值设置为n，从集合AS中取出一个节点集合，执行步骤122；

步骤122、取出的节点集合用AS_j表示，集合AS_j与集合AS中保存的每一个支路节点集合AS_k进行比较，如果集合AS_j和集合AS_k中存在相同的联络开关节点，则从集合AS中取出支路节点集合AS_k，将集合AS_j和集合AS_k合并后，加入到集合NS中；如果集合AS_j与集合AS中保存的每一个支路节点集合AS_k不存在相同的联络开关节点，则将集合AS_j加入到集合NS中；

步骤123、如果集合AS不为空，从集合AS中取出一个节点集合，执行步骤122；如果集合AS为空，执行步骤124；

步骤124、如果集合NS中元素的数量小于UN，则将UN的数值设置为集合NS中元素的数量，将集合NS中的所有元素加入到集合AS中，从集合AS中取出一个支路节点集合，执行步骤122；如果集合NS中元素的数量等于UN，则将集合NS中所有元素加入到集合AS中，完成了节点集合分组；

在步骤14中，主要进行如下步骤；

步骤141、假定集合AS中包含Q个元素，从集合AS中取出一个节点集合，执行步骤142；

步骤142、该节点集合表示为AS_i，如果节点集合AS_i中存在重复的节点，则将节点集合AS_i中重复的节点删除，同一个节点只保留一个；依据节点集合AS_i中支路节点构成的煤矿高压供电系统图统计其包含的支路节点数量PTSZN_i、供电系统级数PTSLN_i，母线最大分段数PTSSN_i、母线节点数量PTSBN_i和联络开关节点数量PTSCN_i，在节点时间权重表PTS查找支路节点数量PTSZN_i、供电系统级数PTSLN_i、母线最大分段数PTSSN_i、母线节点数量PTSBN_i和联络开关节点数量PTSCN_i对应的时间开销K_i，如果能够找到相应的记录，则节点集合AS_i对应的时间开销等于K_i；如果在节点时间权重表PTS找不到相应的记录，则依据节点集合AS_i中节点构成的供电系统图，基于关联矩阵计算获得支路节点与支路节点之间的最终供电关系矩阵所需要的时间K_i，在节点时间权重表PTS中增加一条记录，该记录支路节点数量PTSZN等于支路节点数量PTSZN_i，供电系统级数PTSLN等于供电系统级数PTSLN_i，母线最大分段数PTSSN等于母线最大分段数PTSSN_i，母线节点数量PTSBN等于母线节点数量PTSBN_i，联络开关节点数量PTSCN等于联络开关节点数量PTSCN_i，时间开销字段等于K_i；执行步骤143；

步骤143、如果集合AS不为空，从集合AS中取出一个节点集合，执行步骤142；如何集合AS为空，则所有节点集合AS_i完成拓扑分析需要的总的时间开销