一种基于时序组织型P系统的电网故障诊断方法

一种基于时序组织型p系统的电网故障诊断方法
技术领域
1.本发明属于电网故障诊断技术领域，具体涉及一种基于时序组织型p系统的电网故障诊断方法。


背景技术：

2.电力系统故障诊断旨在通过控制中心获得的警报信号，识别发生故障的元件以及评价保护装置的动作行为。故障警报信息是电力系统故障诊断准确、有效进行的重要前提。电力系统故障警报信号通常从监控和数据采集(supervisory control and data acquisition，scada)系统检索操作信号，但随着电力系统规模的逐步扩大和结构的复杂性日益增加，大量的警报信号涌入调度中心，当这些警报信号过多时，容易让调度员出现误判、漏判等操作。因此，对电力系统快速准确的故障诊断的要求越来越高。
3.到目前为止，国内外已经提出了多种故障诊断方法，例如：专家系统、petri网、人工神经网络等。从现有的大部分故障诊断方法来来，也会面临以下问题：(1)目前大部分故障诊断方法仍然不能很好地处理警报信号的误报、漏报等情况；(2)基于传统组织型p系统的故障诊断模型对保护装置之间动作顺序与逻辑的解释性不强，且其模型复杂度仍然有降低的空间；(3)如今，随着电网规模逐步扩大，电力系统发生故障时需应对多种故障场景，如何在多种故障场景下高效准确的识别出故障元件成为故障诊断方法的关键。
4.电力系统故障诊断的准确性和快速性对于发生故障后快速恢复供电以及维护电网系统安全稳定运行具有重要意义。因此，如何解决好电网实际运行时的故障报警信号的误报、漏报等情况以及描述电力系统故障诊断时各个保护装置之间的动作顺序与逻辑，相关工作亟待研究。


技术实现要素：

5.针对现有技术中的上述不足，本发明提供的基于时序组织型p系统的电网故障诊断方法解决了现有的电网故障诊断中故障信号的误报、漏报及诊断效率较低的问题。
6.为了达到上述发明目的，本发明采用的技术方案为：一种基于时序组织型p系统的电网故障诊断方法，包括以下步骤：
7.s1、通过结线分析法确定待诊断电网中的故障区域，进而确定可疑故障元件；
8.s2、基于scada系统中保护装置的动作信息和时序信息构造并修正初始配置向量w0；
9.s3、根据保护装置之间的动作顺序建立各可疑故障元件的ttps故障诊断模型；
10.s4、根据ttps故障诊断模型建立连接矩阵c和初始配置向量w0，构建转移向量
11.s5、根据ttps故障诊断模型中各可疑故障元件的节点分布，生成故障节点向量s；
12.s6、根据转移向量w1和故障节点向量s，构建判断向量w
end
，并根据判断向量w
end
中的非零元素，确定发出故障的元件，实现电网故障诊断。
13.进一步地，所述步骤s1具体为：
14.s11、设置初始迭代次数k＝1；
15.s12、依次对待诊断电网中的每个元件进行编号，构建元件编号集合p
k
；
16.s13、从集合p
k
中任意取一个元件编号的元件放入元件编号子集合q
k
中，判断待诊断电网中是否存在闭合的断路器与当前选择的元件相连；
17.若是，则进入步骤s14；
18.若否，则进入步骤s15；
19.s14、在待诊断电网中找出所有与该断路器连接的元件，并将其编号加入到集合q
k
中，并返回到步骤s13；
20.s15、使搜索迭代次数k的值增加1；
21.s16、将集合q
k
‑1中的元件编号从集合p
k
‑1中移除，获得新的元件编号集合p
k
；
22.s17、判断当前新的元件编号集合p
k
是否为空；
23.若是，则进入步骤s18；
24.若否，则返回步骤s13；
25.s18、列出元件编号子集合q1,q2,...,q
k
中的所有无源网络，即为故障区域，进而确定可以故障元件；
26.其中，下标k为迭代过程中所获得的元件编号子集合的个数。
27.进一步地，所述步骤s2具体为：
28.s21、基于scada系统中保护装置的动作信息构造初始配置向量w0；
29.s22、判断待诊断电网中的警报信号是否满足时序特性一致性；
30.若是，则进入步骤s23；
31.若否，则进入步骤s24；
32.s23、将当前构造的初始配置向量作为最终的初始配置向量w0，进入步骤s3；
33.s24、通过待诊断电网中警报信号的时序信息对初始配置向量进行修正，得到修正后的初始配置向量w0，进入步骤s3。
34.进一步地，所述步骤s21中的初始配置向量w0表示为w0＝(w
10
,w
20
,...,w
m0
)
t
，w
i0
为组织细胞的初始对象，即代表警报信号的状态值，w
i0
的值为0或者1，1≤i≤m；当初始对象w
i0
所对应的事件发生，即收到该节点i所对应的警报信号，则w
i0
＝1将该节点状态值设为1，当初始对象w
i0
所对应的事件未发生，即未收到该节点i所对应的警报信号，则w
i0
＝0将该节点状态值设为0。
35.进一步地，所述步骤s22中，基于确定警报信号的时序约束和时序推理规则来判断其是否时序特性一致性；其中，警报信号的时序约束包括时间点约束和时间距离约束；
36.所述时间点约束的表达式为：
[0037][0038]
式中，t(i)表示第i个事件在区间内发生，和分别为t(i)的上界和下界，当时，则第i个事件为确定事件，在特定时间点发生；
[0039]
所述时间距离约束的表达式为：
[0040][0041]
式中，d(i,j)表示第i个事件和第j个事件发生时间的间隔必定在时间区间内，和分别表示d(i,j)的下界和上界。
[0042]
进一步地，所述时序推理规则包括前向时序推理和逆向时序推理；
[0043]
其中，前向时序推理为：
[0044]
已知第i个事件的时间点约束，以及其与第j个事件之间的时间距离约束，则推理第j个事件的时间点约束为：
[0045][0046]
逆向时序推理为：
[0047]
已知第j个事件的时间点约束，以及第i个事件与第j个事件的时间距离约束，则推理第j个事件的时间点约束为：
[0048][0049]
进一步地，所述步骤s22中，基于时序约束和时序推理规则确定的时序特性一致性的表达式为：
[0050][0051]
式中，m
i
和m
j
分别为警报信号所对应的第i个事件和第j个事件，t(m
i
)为事件m
i
发生的时间点约束，d(m
i
,m
j
)为事件m
i
和m
j
的时间距离约束。
[0052]
进一步地，所述步骤s3中，一个度为m的ttps故障诊断模型π的表达式为：
[0053]
π＝(o,δ
i
,e,t,i0,syn)
[0054]
式中，o为非空子目标；
[0055]
δ
i
为第i个组织细胞，1≤i≤m，m为组织细胞总数，ttps故障诊断模型π中的组织细胞包括真实细胞和虚拟细胞，其形式分别为δ
i
＝(w
i0
,r
i1
)和δ
i
＝(w
i0
,r
i2
)；
[0056]
w
i0
∈o表示组织细胞中的初始对象多重集；r
i1
表示真实细胞的转运规则，其形式为(i,x/λ,j)，x和λ分别代表细胞i和细胞j中的对象，且细胞i和细胞j分别对应第i个事件和第j个事件，细胞j中的对象λ代表空字符串，执行该规则后，细胞i中的对象x将传递到细胞j中的空字符串λ中；r
i2
表示虚拟细胞的转运规则，其形式为(e,e/λ,vc)，其中vc代表虚拟细胞，表示主保护拒动或断路器拒动，e和λ分别代表组织液环境e和虚拟细胞vc中的对象，虚拟细胞vc中的对象λ代表空字符串；执行该规则后，组织液环境中的对象e将传递到虚拟细胞vc中的空字符串λ中；
[0057]
e＝{e1,...,e
n
}表示组织细胞的组织液环境，其中e
d
∈o,1≤d≤n,代表组织液环境中的对象，n为组织液环境中对象的总个数，当两个细胞之间没有直接的交流通道，则可以通过此环境间接的进行信息交换，当真实组织细胞不满足转运规则r
i1
时，则执行转运规则(i,x/λ,e)，其中x和λ分别代表细胞i和组织液环境e中的对象；执行该规则后，细胞i中的对象x将传递到组织液环境e中的空字符串λ中；
[0058]
t＝{t
δ1
,...,t
δm
}表示各个细胞所对应的时间点集合，其维度为1
×
m，m为细胞总数，t＝{t(t
δi
)|δ
i
∈δ}表示细胞δ
i
的时间点约束集合，表示细胞δ
i
所对应事
件的发生时间应该在时间区间内；d＝{d(δ
i
,δ
j
)|δ
i
,δ
j
∈δ}表示时间距离约束集合，其中表示细胞δ
i
与细胞δ
j
所对应事件的发生时间间隔应该在时间区间内；
[0059]
i0∈{1,...,m}表示输出细胞的标号，每个输出细胞标记为一个可疑故障元件；
[0060]
表示细胞间的连接状态集合，相互连接的细胞可以进行信息交换，对所有细胞(i,j)∈syn,1≤i,j≤m有i≠j。
[0061]
进一步地，所述步骤s4中，建立的连接矩阵c的表达式为：
[0062][0063]
其中，c
ij
表示细胞δ
i
与δ
j
之间存在正向关联，即δ
j
所对应事件是触发δ
i
对应事件的原因，δ
i
是δ
j
发生后产生的结果，c
ij
＝0则表示两个细胞节点之间没有任何关联。
[0064]
进一步地，所述步骤s6中的判断向量w
end
将转移向量w1和故障节点向量s进行逻辑与运算构建出判断向量w
end
，其表达式为：
[0065]
w
end
＝w1∧s
[0066]
式中，w1＝(w
11
,w
21
,...,w
m1
)
t
，w
i1
＝(c
i1
∧w
11
)∨(c
i2
∧w
21
)∨...∨(c
im
∧w
m1
)，i＝1,2,...,m；s＝(s1,s2,...s
n
)
t
，s
i
∈s，s
i
＝0或1，1≤i≤n，当s
i
＝1时代表δ
i
为可疑故障元件，当s
i
＝0时代表为其他设备元件。
[0067]
本发明的有益效果为：
[0068]
(1)本发明在传统组织型p系统基础上那引入了时序因素，提出了一种基于时序组织型p系统(temporal tissue
‑
like p systems,ttps)，并在ttps的诊断模型中，将时间点约束加入到细胞中的对象上，将时间距离约束加入到细胞之间用于信息交流的连接通道上，以此有效利用时序信息进行各个保护装置之间的时序特性一致性约束检查，解决了故障警报信号的误报、漏报等情况；
[0069]
(2)在ttps中引入虚拟细胞和组织液环境，并在此基础上提出可以描述各类保护装置之间动作顺序与逻辑的四种关联关系，进而提出一种新的图形化建模方式，进而提高了模型的解释性，同时降低了模型的复杂度；
[0070]
(3)通过初始配置向量、连接矩阵和故障节点向量之间的逻辑运算推算出故障元件，以应对多种故障场景下的故障诊断，提高了故障诊断效率。因此，利用警报信号的时序特性与组织型p系统相结合，不仅能够较好的处理警报信号中误报、漏报等问题，还可以推理出故障发生时刻所处的时间区间，准确地描述出故障演变过程。
附图说明
[0071]
图1为本发明提供的基于时序组织型p系统的电网故障诊断方法流程图。
[0072]
图2为本发明提供的简单输电系统及ttps故障诊断模型示意图。
[0073]
图3为本发明提供的ttps故障诊断模型中的四种关联关系示意图。
[0074]
图4为本发明提供的标准ieee14节点示意图。
[0075]
图5为本发明提供的母线b13的ttps故障诊断模型示意图。
具体实施方式
[0076]
下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
[0077]
实施例1：
[0078]
本发明实施例提供的基于时序组织型p系统的电网故障诊断方法，包括以下步骤：
[0079]
s1、通过结线分析法确定待诊断电网中的故障区域，进而确定可疑故障元件；
[0080]
s2、基于scada系统中保护装置的动作信息和时序信息构造并修正初始配置向量w0；
[0081]
s3、根据保护装置之间的动作顺序建立各可疑故障元件的ttps故障诊断模型；
[0082]
s4、根据ttps故障诊断模型建立连接矩阵c和初始配置向量w0，构建转移向量
[0083]
s5、根据ttps故障诊断模型中各可疑故障元件的节点分布，生成故障节点向量s；
[0084]
s6、根据转移向量w1和故障节点向量s，构建判断向量w
end
，并根据判断向量w
end
中的非零元素，确定发出故障的元件，实现电网故障诊断。
[0085]
需要说明的是，在执行上述步骤的过程中，步骤s2和s3均是基于步骤s1获得的结果进行的，且步骤s2和步骤s3是并列关系，本实施例中将其表示为s2～s3的顺序仅是为了形成连贯的执行逻辑及便于理解本技术方案，实际实现过程中，可根据需要同时执行步骤s2和s3。
[0086]
本发明实施例的步骤s1具体为：
[0087]
s11、设置初始迭代次数k＝1；
[0088]
s12、依次对待诊断电网中的每个元件进行编号，构建元件编号集合p
k
；
[0089]
s13、从集合p
k
中任意取一个元件编号的元件放入元件编号子集合q
k
中，判断待诊断电网中是否存在闭合的断路器与当前选择的元件相连；
[0090]
若是，则进入步骤s14；
[0091]
若否，则进入步骤s15；
[0092]
s14、在待诊断电网中找出所有与该断路器连接的元件，并将其编号加入到集合q
k
中，并返回到步骤s13；
[0093]
s15、使搜索迭代次数k的值增加1；
[0094]
s16、将集合q
k
‑1中的元件编号从集合p
k
‑1中移除，获得新的元件编号集合p
k
；
[0095]
s17、判断当前新的元件编号集合p
k
是否为空；
[0096]
若是，则进入步骤s18；
[0097]
若否，则返回步骤s13；
[0098]
s18、列出元件编号子集合q1,q2,...,q
k
中的所有无源网络，即为故障区域，进而确定可以故障元件；
[0099]
其中，下标k为迭代过程中所获得的元件编号子集合的个数。
[0100]
在本发明实施例的步骤s2具体为：
[0101]
s21、基于scada系统中保护装置的动作信息构造初始配置向量w0；
[0102]
s22、判断待诊断电网中的警报信号是否满足时序特性一致性；
[0103]
若是，则进入步骤s23；
[0104]
若否，则进入步骤s24；
[0105]
s23、将当前构造的初始配置向量作为最终的初始配置向量w0，进入步骤s3；
[0106]
s24、通过待诊断电网中警报信号的时序信息对初始配置向量进行修正，得到修正后的初始配置向量w0，进入步骤s3。
[0107]
在上述步骤s21中，由scada系统提供初始状态信息，以此作为基于ttps故障诊断方法的初始出入，初始配置向量w0表示为w0＝(w
10
,w
20
,...,w
m0
)
t
，w
i0
为组织细胞的初始对象，即代表警报信号的状态值，w
i0
的值为0或者1，1≤i≤m；当初始对象w
i0
所对应的事件发生，即收到该节点i所对应的警报信号，则w
i0
＝1将该节点状态值设为1，当初始对象w
i0
所对应的事件未发生，即未收到该节点i所对应的警报信号，则w
i0
＝0将该节点状态值设为0。
[0108]
上述步骤s22中，基于确定警报信号的时序约束和时序推理规则来判断其是否时序特性一致性；其中，警报信号的时序约束包括时间点约束和时间距离约束；
[0109]
所述时间点约束的表达式为：
[0110][0111]
式中，t(i)表示第i个事件在区间内发生，和分别为t(i)的上界和下界，当时，则第i个事件为确定事件，在特定时间点发生；
[0112]
所述时间距离约束的表达式为：
[0113][0114]
式中，d(i,j)表示第i个事件和第j个事件发生时间的间隔必定在时间区间内，和分别表示d(i,j)的下界和上界；
[0115]
由于时间点约束和时间距离约束均采用区间进行定义，因此同样基于区间的计算规则来定义时序特性约束推理规则，本实施例中的时序推理规则包括前向时序推理和逆向时序推理；
[0116]
其中，前向时序推理为：
[0117]
已知第i个事件的时间点约束，以及其与第j个事件之间的时间距离约束，则推理第j个事件的时间点约束为：
[0118][0119]
逆向时序推理为：
[0120]
已知第j个事件的时间点约束，以及第i个事件与第j个事件的时间距离约束，则推理第j个事件的时间点约束为：
[0121][0122]
基于上述约束及推理规则可知，警报信号之间具有时序逻辑关联，这种时序逻辑关系可用于判断保护装置或断路器误动，以及警报信号误报等不确定情况；基于时序约束和时序推理规则确定的时序特性一致性的表达式为：
[0123]
[0124]
式中，m
i
和m
j
分别为警报信号所对应的第i个事件和第j个事件，t(m
i
)为事件m
i
发生的时间点约束，d(m
i
,m
j
)为事件m
i
和m
j
的时间距离约束。
[0125]
上述时序一致性表达式的含义为：若事件m
j
发生于t(m
j
)时刻，则t(m
j
)需满足时序特性一致性约束，即t(m
j
)∈(t(m
i
)+d(m
i
,m
j
))。若满足，则事件m
j
的发生时刻t(m
j
)在规定时间限制内完成；否则，相应的警报信号属于错误警报信息。
[0126]
在本实施例的步骤s3中，一个度为m的ttps故障诊断模型π的表达式为：
[0127]
π＝(o,δ
i
,e,t,i0,syn)
[0128]
式中，o为非空子目标；
[0129]
δ
i
为第i个组织细胞，1≤i≤m，m为组织细胞总数，ttps故障诊断模型π中的组织细胞包括真实细胞和虚拟细胞，其形式分别为δ
i
＝(w
i0
,r
i1
)和δ
i
＝(w
i0
,r
i2
)；
[0130]
w
i0
∈o表示组织细胞中的初始对象多重集；r
i1
表示真实细胞的转运规则，其形式为(i,x/λ,j)，x和λ分别代表细胞i和细胞j中的对象，且细胞i和细胞j分别对应第i个事件和第j个事件，细胞j中的对象λ代表空字符串，执行该规则后，细胞i中的对象x将传递到细胞j中的空字符串λ中；r
i2
表示虚拟细胞的转运规则，其形式为(e,e/λ,vc)，其中vc代表虚拟细胞，表示主保护拒动或断路器拒动，e和λ分别代表组织液环境e和虚拟细胞vc中的对象，虚拟细胞vc中的对象λ代表空字符串；执行该规则后，组织液环境中的对象e将传递到虚拟细胞vc中的空字符串λ中；
[0131]
e＝{e1,...,e
n
}表示组织细胞的组织液环境，其中e
d
∈o,1≤d≤n,代表组织液环境中的对象，n为组织液环境中对象的总个数，当两个细胞之间没有直接的交流通道，则可以通过此环境间接的进行信息交换，当真实组织细胞不满足转运规则r
i1
时，则执行转运规则(i,x/λ,e)，其中x和λ分别代表细胞i和组织液环境e中的对象；执行该规则后，细胞i中的对象x将传递到组织液环境e中的空字符串λ中；
[0132]
t＝{t
δ1
,...,t
δm
}表示各个细胞所对应的时间点集合，其维度为1
×
m，m为细胞总数，t＝{t(t
δi
)|δ
i
∈δ}表示细胞δ
i
的时间点约束集合，表示细胞δ
i
所对应事件的发生时间应该在时间区间内；d＝{d(δ
i
,δ
j
)|δ
i
,δ
j
∈δ}表示时间距离约束集合，其中表示细胞δ
i
与细胞δ
j
所对应事件的发生时间间隔应该在时间区间内；
[0133]
i0∈{1,...,m}表示输出细胞的标号，每个输出细胞标记为一个可疑故障元件；
[0134]
表示细胞间的连接状态集合，相互连接的细胞可以进行信息交换，对所有细胞(i,j)∈syn,1≤i,j≤m有i≠j。
[0135]
上述ttps故障诊断模型如图2所示，根据scada系统所提供的保护装置的动作信息，采用ttps建模展现出可疑故障元件和保护装置之间的因果关系。下面，以如图2(a)所示的简单输电系统为例，对本发明方法的建模思路进行阐述，该系统中包括一条线路l、两个断路器cb1和cb2、两个主保护器mlr1和mlr2和两个后备保护blr1和blr2。线路的保护配置为：若线路l发生故障，则mlr1和mlr2动作，跳开断路器cb1和cb2；若mlr1和mlr2未能正常动作，则blr1和blr2将动作跳开断路器cb1和cb2。
[0136]
基于ttps建立此输电系统的故障诊断模型，如图2(b)所示，其中实线有向弧代表
各个保护装置之间的因果关系，其有向弧上的时间区间表示关联关系不同的连接通道具有不同的时间距离约束区间，蓝色虚线代表进行推理算法的方向。该模型由9个细胞组成，可疑故障元件与其保护装置之间的因果关系逻辑由ttps进行建模。首先利用上传到scada系统的保护和断路器的动作信息，判定停电区域中的线路l为可疑元件，其中线路l所对应的细胞δ1被指定为输出细胞，vc(δ1,δ2)表示细胞δ1所代表的事件发生，而细胞δ2所代表的事件不发生。
[0137]
考虑时序特性之后，ttps中的组织细胞和连接组织细胞的连接通道就具备了时间属性。而关联关系不同的连接通道具有不同的时序属性，因此，定义了ttps的四种关联关系。
[0138]
ttps的四种关联关系定义如图3所示:
[0139]
ttps与tps相比，其引入了时间点约束和时间距离约束来提高诊断模型的容错性，使正向推理和时间推理可以一起进行，此外，ttps还引入了虚拟细胞(virtual cell，vc)，以更好地描述电力系统故障诊断的动态逻辑。
[0140]
为了更好地理解ttps模型，对ttps中的四种关联关系进行详细描述。图3
‑
(a)表示故障设备节点δ
j
引发相应的保护节点δ
i
动作；图3
‑
(b)表示保护节点δ
j
引发相应断路器细胞δ
i
动作；图3
‑
(c)表示断路器细胞δ
j
拒动引发相应失灵保护细胞δ
i
动作，vc细胞表示断路器细胞δ
j
拒动引发的细胞；图3
‑
(d)表示主保护细胞δ
j
拒动引发相应的后备保护细胞δ
i
动作，vc细胞表示主保护细胞δ
j
拒动引发的细胞。
[0141]
在ttps中，关联关系不同的连接通道具有不同的时间区间属性，相应的时间距离约束为：
[0142]
(1)(故障发生，主保护动作)∈[10,40]ms；
[0143]
(2)(保护动作，断路器动作)∈[40,60]ms；
[0144]
(3)(断路器拒动，失灵保护动作)∈[210,240]m；
[0145]
(4)(主保护拒动，后备保护动作)∈[940,1030]ms。在本实施例的步骤s4中，为了更好地描述各个细胞的通信规则和因果关系，引入连接矩阵c，其表达式为：
[0146][0147]
其中，c
ij
表示细胞δ
i
与δ
j
之间存在正向关联，即δ
j
所对应事件是触发δ
i
对应事件的原因，δ
i
是δ
j
发生后产生的结果，c
ij
＝0则表示两个细胞节点之间没有任何关联。
[0148]
在本实施例的步骤s4中，转移向量w1为：
[0149][0150]
其中，w1＝(w
11
,w
21
,...,w
m1
)
t
，w
i1
＝(c
i1
∧w
11
)∨(c
i2
∧w
21
)∨...∨(c
im
∧w
m1
)，i＝1,2,...,m；对于或1，定义a∧b＝min(a,b)和a∨b＝max(a,b)；
[0151]
w0＝(w
10
,w
20
,...,w
m0
)
t
代表由scada系统提供的保护装置的动作信息构造的初始配置向量w0；
[0152]
c＝(c
ij
)
m
×
m
，c
ij
＝0或1，代表各个事件所对应细胞的通信规则和因果关系；
[0153]
本实施例的步骤s5中，获取ttps中各个故障元件的节点分布，根据其分布生成故障节点向量s，其反映了ttps中各可疑故障元件的节点分布，s＝(s1,s2,...s
n
)
t
，s
i
∈s，s
i
＝0或1，1≤i≤n，当s
i
＝1时代表δ
i
为可疑故障元件，当s
i
＝0时代表为其他设备元件；
[0154]
本实施例的步骤s6中的判断向量w
end
将转移向量w1和故障节点向量s进行逻辑与运算构建出判断向量w
end
，其表达式为：
[0155]
w
end
＝w1∧s
[0156]
实施例2：
[0157]
本实施例中，以标准ieee14节点系统为诊断对象，以算例1为例，给出具体计算过程，以促进细节理解，图4为标准ieee14节点网络拓扑图。
[0158]
本实施例中的预设故障情景和故障诊断元件结果如表1所示
[0159]
表1算例1的预设故障情景与故障元件诊断结果
[0160][0161]
首先，在故障发生后，利用结线分析法判断出停电区域，得到可疑故障元件母线b13、线路l1213和线路l1314。建立母线b13的ttps故障诊断模型，如图5所示：
[0162]
基于图5母线b13的ttps故障诊断模型，建立14
×
14的连接矩阵c，表示如下：
[0163][0164]
根据scada系统提供的保护装置和断路器的动作信息，初始配置向量设置为w0＝(0,1,1,0,0,0,1,0,0,0,1,1,1,1)
t
。
[0165]
经时序一致性约束检查，所有保护装置均符合时序一致性约束条件，无需修正初
始配置向量w0＝(0,1,1,0,0,0,1,0,0,0,1,1,1,1)
t
。
[0166]
连接矩阵c和初始配置向量w0计算得到转移向量w1。
[0167][0168]
生成故障节点向量s，s中非0元素为母线b13和虚拟细胞δ
12
。
[0169]
将生成的故障节点向量s与转移向量w1进行逻辑“与”运算得到判断向量w
end
。
[0170]
w
end
＝w1∧s＝(1,0,0,0,0,0,0,0,0,0,0,0,0,0)
[0171]
判断得出母线b13发生了故障。
[0172]
同理判断得出线路l1213，断路器cb1314发生拒动。
[0173]
分析故障演变过程，以接收到的第一个警报的时标为基准点，通过分析时序推理，实际故障演变过程为：母线b13在[
‑
20,
‑
10]期间发生故障，0ms时母线主保护b13m动作并引发断路器cb1312和cb1306分别在48ms和51ms动作，断路器cb1314拒动在1989ms引发远后备保护l1413s，并在2041ms跳闸断路器cb1413。