CN104049175A - 具有多个多相母线的电网中的故障分析 - Google Patents

Abstract

本发明提出一种用于具有多个多相母线的电网中的故障分析。通过针对感兴趣的多相母线计算等效电路来执行用于具有多个多相母线的电网的故障分析。每个等效电路包括N×N阻抗矩阵，N对应于针对其计算该等效电路的多相母线的相位的数目。基于多个实验电流注入矢量所产生的电压以及因式分解的故障前导纳矩阵确定阻抗矩阵的元素。故障前导纳矩阵代表无故障的多相母线的节点导纳。针对感兴趣的每个多相母线和故障类型，基于针对该多相母线确定的等效电路，确定故障电流注入矢量。针对感兴趣的每个多相母线和故障类型，基于因式分解的故障前导纳矩阵以及针对该多相母线和故障类型确定的故障电流注入矢量，确定故障电压矢量。

Description

具有多个多相母线的电网中的故障分析

技术领域

本申请涉及电网，并且更具体地涉及用于电网的故障分析。

背景技术

用于电力配电系统的配电管理系统(DMS)使用配电网络的系统模型以及假设的故障位置、类型和故障阻抗通过在测量的故障电流与预测的故障电流之间的模式匹配来识别故障位置。常规的基于对称分量的方法不适合于在构建和操作两者中都不平衡的配电系统。需要基于相位的方法用于准确分析。然而由于个别相位的显式表示所产生的大量节点，即使对于中等馈线，所得系统仍然可能很大。基于阻抗的方法既存储密集又耗费时间。基于导纳矩阵的方法在存储效率上高得多，但是耗费时间。基于补偿和基于矩阵求逆引理的方法据称减少与反复地对导纳矩阵求逆相关联的计算负担。然而这些常规方法对于在具有各种故障类型的许多母线位置处的反复故障分析而言计算密集并且涉及到构建和因式分解导纳矩阵，或者甚至为每个故障位置和故障类型对导纳矩阵求逆。使用对称分量的传统方法也可以处理不平衡故障，但是不能处理不平衡网络建模。

发明内容

根据本文描述的示例实施例，可以在电网中、比如在电力输电或者配电系统中的许多位置处在很高效的过程中执行反复故障分析，而仅有故障前导纳矩阵的单次因式分解而没有因式分解的后续更新。所研究的系统和所分析的故障两者都可以是平衡的或者不平衡的。运用故障前系统的稀疏导纳建模以执行故障分析，无论分析的故障类型和位置的数目如何，该故障分析都使用相同的因式分解的故障前导纳矩阵。

根据用于具有多个多相母线的电网的故障分析的方法的一个实施例，该方法包括：针对感兴趣的多相母线计算等效电路，每个等效电路包括N×N阻抗矩阵，其中N对应于针对其计算等效电路的多相母线的相位的数目；基于多个实验电流注入矢量所产生的电压以及因式分解的故障前导纳矩阵，确定阻抗矩阵的元素，故障前导纳矩阵代表无故障的多相母线的节点导纳；针对感兴趣的每个多相母线和故障类型，基于针对该多相母线确定的等效电路，确定故障电流注入矢量；并且针对感兴趣的每个多相母线和故障类型，基于因式分解的故障前导纳矩阵以及针对该多相母线和故障类型确定的故障电流注入矢量，确定故障电压矢量。

根据存储计算机程序的非瞬态计算机可读介质的一个实施例，该计算机程序可操作为执行用于具有多个多相母线的电网的故障分析，该计算机程序包括针对感兴趣的多相母线计算等效电路的程序指令，每个等效电路包括N×N阻抗矩阵，其中N对应于针对其计算该等效电路的多相母线的相位的数目。该计算机程序进一步包括用于基于多个实验电流注入矢量所产生的电压以及因式分解的故障前导纳矩阵确定阻抗矩阵的元素的程序指令，故障前导纳矩阵代表无故障的多相母线的节点导纳。该计算机程序也包括：用以针对感兴趣的每个多相母线和故障类型，基于针对该多相母线确定的等效电路确定故障电流注入矢量的程序指令；以及用以针对感兴趣的每个多相母线和故障类型，基于因式分解的故障前导纳矩阵以及针对该多相母线和故障类型所确定的故障电流注入矢量来确定故障电压矢量的程序指令。

根据与具有多个多相母线的电网通信的计算机系统的一个实施例，该计算机系统包括处理电路，该处理电路可操作用于针对感兴趣的多相母线计算等效电路，每个等效电路包括N×N阻抗矩阵，其中N对应于针对其计算该等效电路的多相母线的相位的数目。该处理电路进一步可操作用于基于多个实验电流注入矢量所产生的电压以及因式分解的故障前导纳矩阵确定阻抗矩阵的元素，故障前导纳矩阵代表无故障的多相母线的节点导纳。该处理电路也可操作用于：针对感兴趣的每个多相母线和故障类型，基于针对该多相母线确定的等效电路确定故障电流注入矢量；并且针对感兴趣的每个多相母线和故障类型，基于因式分解的故障前导纳矩阵以及针对该多相母线和故障类型确定的故障电流注入矢量来确定故障电压矢量。

本领域技术人员将在阅读以下具体描述并且在查看附图时认识附加特征和优点。

附图说明

图中的部件不必按比例，代之重点放在图示本发明的原理。另外，在图中，相似附图标记标示对应部分。在附图中：

图1图示包括多个多相母线的电网和用于实现针对网络的故障分析的控制系统的示意图。

图2图示在电网的不同母线之间的各种相位连接的示意图。

图3图示在电网的示例性三相母线的节点k与1之间的线模型的示意图。

图4A至4D图示用于电网的各种故障类型的示意图。

图5图示用于包括多个多相母线的电网的故障分析方法的一个实施例的流程图。

图6图示建模为多相戴维南等效(TE)电路的感兴趣的故障母线的示意图。

图7图示用于包括多个多相母线的电网的故障分析方法的另一实施例的流程图。

具体实施方式

图1图示可以具有网状或者放射状架构的电网，比如电力输电或者配电系统的一个非限制示例性实施例。电网具有一个或者多个变电站100和连接到变电站100用于提供功率的多个母线102。母线102中的一些母线是多相，即三相或者两相。母线120中的其它母线可以是单相。图1中的具有单端箭头的线代表的支线可以连接到母线102用于向各种服务变压器(为了易于图示而未示出)分配功率。在电网中的许多位置处执行反复故障分析。运用故障前系统的稀疏导纳建模以执行故障分析，无论分析的故障类型和位置数目如何，该故障分析都使用相同因式分解的故障前导纳矩阵。

可以使用本文描述的故障分析方法作为用于在线故障位置分析和保护协调分析的支持功能。分析中的网络和故障二者可以是平衡或者不平衡的。该方法也可以应用于输电系统，其中系统在除了故障点之外的相位之间平衡。虽然本文就单个位置短路故障(在指定的母线处的各种类型的短路故障)方面描述故障分析过程，但是该过程也可以在简单的扩展下应用于在多个位置处的同时故障以及系列故障。

本文描述的方法减少在每个位置处的故障分析的计算工作。该方法使用另一应用已经可用的，或者计算一次并再使用的因式分解的故障前导纳矩阵，来执行最少增加的操作。因而，在电网内的许多不同位置处执行的故障分析被降低至非常流线型过程的针对每个故障位置一次的反复操作，该流线型过程涉及数目非常有限的简单计算。该过程因此可以在许多系统位置处很高效地执行故障分析。该过程对于在多个位置处的同时故障而言类似地有效。

更具体而言，提供控制系统104用于监视和控制电网的操作，包括故障分析。控制系统104可以经由如图1中所示虚线连接指示的有线或者无线连接来连接到电网，或者包含于变电站100中的一个中。控制系统104包括处理电路106，该处理电路可以包括诸如一个或者多个控制器、处理器、ASIC(专用集成电路)等用于运行执行诸如故障分析的网络控制功能的程序代码的数字和/或模拟电路。为此，控制系统104包括在处理电路106中包括的或者与处理电路106相关联的用于执行本文描述的故障分析过程和对应计算的故障分析器108。控制系统104也具有用于存储程序代码以及在程序代码运行期间通过处理电路106和故障分析器108处理和访问有关数据的存储器110，比如DRAM(动态随机存取存储器)和HDD(硬盘驱动)112。控制系统104也具有用于发送和接收信息，包括从电网接收可以在故障分析过程中用来识别网络中的可能故障的电压和电流测量信息的I/O(输入/输出)电路114。

由控制系统104实施的故障分析方法可以例如由于系统构建、操作或者故障中的相位不平衡，对母线102的个别相位建模。图2示出在电网的不同母线(母线i、母线j、母线k、母线l等)与用于每个单独的线路区段的对应自阻抗(等)模型之间的各种相连接(相位a、相位相b、相位c)以及在相同多相线路区段(等)的相位之间的互耦合。图3示出在示例性三相母线的节点k与l之间的对应线路模型。每相具有在节点k的电压和电流以及在节点l的电压和电流可以在母线的节点k与l之间将线路阻抗建模为和而对应相电流为和将在相位a与b之间的互阻抗建模为将在相位b与c之间的互阻抗建模为并且将在相位a与c之间的互阻抗建模为也可以针对每个节点针对相位与地之间的等效导纳(Y)和电流建模。

由控制系统104实施的故障分析方法可以针对电网的母线102的单独的相位对各种类型的短路故障建模。例如图4A至4D图示可以在三相母线上出现的各种类型的短路故障，其中Z_f1、Z_f2、Z_f3和Z_e代表故障阻抗。图4A示出三相接地短路，图4B示出单线接地故障，图4C示出双线接地故障，并且图4D示出线间故障。

图5图示由控制系统104执行的故障分析的一个实施例。该方法包括针对感兴趣的多相母线102计算等效电路，每个等效电路包括N×1电压源矢量和N×N阻抗矩阵Z_TE，其中N对应于针对其计算该等效电路的多相母线102的相位数目(块200)。在一个实施例中，等效电路是戴维南等效电路。利用简单的变换，可以构造多相诺顿等效，如果这样做是有益的或者另外更有利则在基于诺顿等效电路在后续故障计算中使用多相诺顿等效。

图6示出在感兴趣的故障母线处的多相戴维南等效(TE)电路表征的系统行为。由下式给出的三相电压源矢量：

$E=\left(\begin{array}{c}{E}_a\\ E_b\\ E_c\end{array}\right)---\left(1\right)$

以及下式给出的N×N阻抗矩阵：

$Z_{\mathrm{TE}}=\left(\begin{array}{ccc}{z}_{\mathrm{aa}}& z_{\mathrm{ab}}& z_{\mathrm{ac}}\\ z_{\mathrm{ba}}& z_{\mathrm{bb}}& z_{\mathrm{bc}}\\ z_{\mathrm{ca}}& z_{\mathrm{cb}}& z_{\mathrm{cc}}\end{array}\right)---\left(2\right)$

代表TE电路，其中N对应于相位的数目。Z_TE的每个对角线元素代表自阻抗，并且每个对角线外元素代表用于感兴趣的故障母线的戴维南等效电路的相间互阻抗。例如Zaa代表相位a的自阻抗，Zbb代表代表相位b的自阻抗，Zcc代表相位c的自阻抗，Zab代表相位a与b之间的互阻抗，Zac代表相位a与c之间的互阻抗等等。矢量E是戴维南等效电压源并且从故障前条件获得，并且不随故障位置而改变。

回到图5，故障分析方法进一步包括基于针对感兴趣的每个故障母线102的多个实验电流注入矢量所产生的电压以及因式分解的故障前导纳矩阵Y_pre，针对该感兴趣的每个故障母线102确定阻抗矩阵(Z_TE)的元素(块200)。这些电流注入矢量是实验性的，因为未向电网中注入实际电流。取而代之，这些电流注入矢量是用于在某些节点电流条件之下对网络行为建模(为等效电路)的假设的或者仿真的电流注入矢量。在电力工程设计中的导纳矩阵是描述具有n个母线的电力系统的N×N矩阵，并且代表无故障的母线的节点导纳。在包含许多母线的实际电网中，导纳矩阵相当稀疏。

更具体而言，I_pre是复数节点源电流矢量，并且V_pre表示复数故障前节点电压矢量。源电流矢量I_pre代表电网的电流源节点，例如变电站节点。故障前复数导纳矩阵Y_pre可以使用任何适当常规技术由控制系统104计算或者被提供至控制系统104，用于在实施故障分析方法时使用。针对电网的导纳矩阵的计算为电力工程设计领域技术人员所知，因此这里未提供进一步说明。

故障前导纳矩阵Y_pre从另一应用已经可用或者由控制系统104计算一次并且再次使用，以执行用于电网的故障分析。为了符号表示方便，可以针对相同母线102首先对所有源节点(例如变电站节点)编号并且按照相位(a、b等)川页序。对于仅有单个三相源的网络，I_pre和V_pre由下式给出：

I_pre＝[i₁，i₂，i₃，0，0，0，...，0，0]^T                                 (3)

V_pre＝[v₁，v₂，v₃，v₄，v₅，...，v_n]^T                                  (4)

这些等式满足：

Y_preV_pre＝I_pre                                 (5)

$E=\left(\begin{array}{c}{E}_a\\ E_b\\ E_c\end{array}\right)---\left(5a\right)$

以及

V_pre＝[v₁，v₂，v₃，v₄，v₅，...，v_n]^T                                 (5b)

控制系统104可以对故障前导纳矩阵Y_pre进行因式分解，或者这样的因式分解可以已经完成。在一个实施例中，Y_pre的LU因式分解产生下三角矩阵L、上三角矩阵U、置换矩阵P和列重新排序矩阵Q，每个都从故障前导纳矩阵Y_pre因式分解。然后针对感兴趣的故障母线102可以获得阻抗矩阵Z_TE的第一列，一般通过求解下式：

YV＝I_a                                 (6)

并且在一个实施例中更具体地通过求解下式：

V＝[v₁，v₂，v₃，v₄，...，，v_a，v_b，v_c]^T＝Solve(L，U，P，Q，I_a)                           (7)其中Solve(L，U，P，Q，I_a)代表针对指定的实验电流注入矢量I_a，使用已经因式分解的导纳矩阵Y_pre来求解电压矢量V。电压矢量V是故障前导纳矩阵Y_pre的因式分解的分量L、U、P、Q的函数，并且I_a表示在分析之下的针对故障母线102的节点的相位的实验电流注入矢量。

对于Z_TE的第一列，实验电流注入矢量I_a具有如下式给出的分析中的故障母线102的节点的在相位的单位注入电流和在所有其它节点处的零注入：

Ia＝[0，0，0，0，0，0，...，0，0，１，0，0]^T                           (8)

使用这样具有单位电流注入的实验电流注入矢量产生关于Z_TE的第一列如下：

$Z_{\mathrm{TE},a}=\left[\begin{array}{c}{z}_{\mathrm{aa}}\\ z_{\mathrm{ba}}\\ z_{\mathrm{ca}}\end{array}\right]=\left[\begin{array}{c}{v}_a\\ v_b\\ v_c\end{array}\right]---\left(9\right)$

可以通过节点b和节点c处的单位电流注入用相似方式针对分析中的每条故障母线102获得Z_TE的其它两列，分别如下式给出：

I_b＝[0，0，0，0，0，0，...，0，0，0，1，0]^T                           (10)

V＝[v₁，v₂，v₃，v₄，...，，v_a，v_b，v_c]^T＝Solve(L，U，P，Q，I_b)                 (11)

$Z_{\mathrm{TE},b}=\left[\begin{array}{c}{z}_{\mathrm{ab}}\\ z_{\mathrm{bb}}\\ z_{\mathrm{cb}}\end{array}\right]=\left[\begin{array}{c}{v}_a\\ v_b\\ v_c\end{array}\right]---\left(12\right)$

I_c＝[0，0，0，0，0，0，...，0，0，0，0，1]^T                               (13)

V＝[v₁，v₂，v₃，v₄，...，，v_a，v_b，v_c]^T＝Solve(L，U，P，Q，I_c)           (14)

$Z_{\mathrm{TE},c}=\left[\begin{array}{c}{z}_{\mathrm{ac}}\\ z_{\mathrm{bc}}\\ z_{\mathrm{cc}}\end{array}\right]=\left[\begin{array}{c}{v}_a\\ v_b\\ v_c\end{array}\right]---\left(15\right)$

其中Solve(L，U，P，Q，I_x)代表针对指定的实验电流注入矢量，使用已经因式分解的导纳矩阵来求解电压矢量V。电压矢量V是故障前导纳矩阵Y_pre的因式分解的分量L、U、P、Q以及I_x的函数，该I_x具有在分析之下的故障母线的相位x节点处的单位电流注入和在所有其它节点处的零注入。也就是说，针对感兴趣的每个多相母线102确定的实验电流注入矢量I_x具有与该母线102的每个相位对应的单位值元素以及其他的零值元素。

再次回到图5，一旦针对感兴趣的多相母线102获得等效电路(Z_TE)，可以针对感兴趣的每个多相母线102和故障类型，基于针对该多相母线102确定的等效电路，来确定故障电流注入矢量I_f(块210)。可以使用任何已知方法针对感兴趣的不同故障类型计算故障电流。对于故障类型中的许多故障类型，可以容易推导分析表达式。本文省略针对节点故障电流的计算过程的进一步细节，因为可以使用许多方式，这些方式中的每种方式对于具有基本工程背景的技术人员所熟知并且由先前确定的等效电路给出。基于等效电路所计算的故障电流可以如下式给出的3×1矢量中表示：

$I_{f,\mathrm{eq}.}=\left[\begin{array}{c}{I}_{f,\mathrm{eq}.-a}\\ I_{f,\mathrm{eq}.-b}\\ I_{f,\mathrm{eq}.-c}\end{array}\right]---\left(16\right)$

来自等式(16)的故障电流矢量用于构造应用于故障前导纳矩阵的故障后电流注入矢量I_f。I_f是针对感兴趣的每个多相母线102和故障类型计算的并且是等式(16)代表的(向右维度展开的)故障后电流矢量I_f，eq.与等式(3)代表的故障前源电流注入矢量I_pre之和，如下式给出：

I_f＝[i₁，i₂，i₃，0，0，0，...，0，0]^T+[0，0，0，，...，0，0，i_f，eq.-a，_if，eq.-b，i_f，eq.-c]^T   (17)

针对感兴趣的每个多相母线102和故障类型确定的故障电流注入矢量I_f具有与在故障前条件中的源母线以及与该故障母线的相位对应的非零值元素以及其他的零值元素。

故障分析方法继续用控制系统104针对感兴趣的每个多相母线102和故障类型基于针对该多相母线102和故障类型确定的故障电流注入矢量I_f以及因式分解后的故障前导纳矩阵Y_pre确定故障电压矢量V_f(块220)。可以如下式给出的那样针对感兴趣的每个多相母线102和故障类型确定故障电压矢量V_f：

V_f＝[v₁，v₂，v₃，v₄，v₅，...，，v_n]^T＝Solve(L，U，P，Q，I_f)              (18)

其中Solve(L，U，P，Q，I_f)代表求解V_f作为故障前导纳矩阵Y_pre的因式分解后的分量L、U、P、Q以及针对分析中的故障节点母线102确定的故障电流注入矢量I_f的函数。

利用在感兴趣的所有母线102处已知的故障后电压，可以针对在电网内的感兴趣的任何支路使用基本等式来计算支路上的故障电流。可以通过控制系统104基于在假设的故障位置处并且针对不同故障类型在测量的故障电流与计算的故障电流之间的模式匹配来执行故障位置分析。在一个实施例中，在电网中测量并且经由I/O电路114向控制系统104报告一个或者多个故障电流。处理电路106对一个或者多个测量的故障电流与故障分析器108确定的计算的故障电流进行比较，以识别电网中的一个或者多个可能故障条件。这样的方式在速度重要的联机应用中特别有利。

图7图示本文描述的故障分析方法的更详细的流程图。在块300计算或者查明故障前导纳矩阵Y_pre并且在块302因式分解故障前导纳矩阵Y_pre。在块304根据等式(5)基于故障前源电流矢量I_pre和故障前导纳矩阵Y_pre计算故障前电压矢量V_pre。控制系统104然后在感兴趣的所有母线102上迭代故障分析过程(块306)。在感兴趣的所有母线上的迭代涉及到在块308处根据等式(6)和(7)针对感兴趣的每个母线102计算故障前电压，然后针对感兴趣的每个故障母线102在阻抗矩阵Z_TE的所有列上迭代(块310)。针对感兴趣的每个故障母线102在阻抗矩阵Z_TE的所有列上迭代包括在块312处根据等式(8)针对故障母线102的每个相位节点设置电流注入矢量I_x并且在块314处根据等式(8)求解对应节点电压。此时，针对感兴趣的多相母线102获得等效电路(E，Z_TE)。

故障分析过程然后在感兴趣的所有故障类型上迭代以完成该过程(块316)。在感兴趣的所有故障类型上迭代包括根据等式(16)求解在每个故障节点处的故障电流(块318)，根据等式(17)设置故障后电流注入矢量I_f(块320)并且根据等式(18)求解对应的故障后电压(块322)。利用在感兴趣的所有母线102处已知的故障后电压，可以使用电网内的感兴趣的任何支路的基本等式来计算支路上的故障电流。

术语如“第一”、“第二”等用来描述各种单元、区域、区段等并不旨在于限制。贯穿说明书中相似术语指代相似单元。

如这里所用，术语“具有”、“包含”、“包括”等是开放式术语，这些指示陈述的单元或者特征的存在，但是不排除附加单元或者特征。除非上下文另有明示，冠词“一”、“一个”和“该”旨在于包括复数以及单数。

谨记上述变化和应用的范围，应当理解本发明既不受前文描述限制，也不受附图限制。取而代之，本发明仅受所附权利要求及其法律等同限制。

Claims (20)

1.一种用于具有多个多相母线的电网的故障分析的方法，所述方法包括：

针对感兴趣的所述多相母线计算等效电路，每个等效电路包括N×N阻抗矩阵，其中N对应于针对其计算所述等效电路的所述多相母线的相位的数目；

基于多个实验电流注入矢量所产生的电压以及因式分解的故障前导纳矩阵，确定所述阻抗矩阵的元素，所述故障前导纳矩阵代表无故障的所述多相母线的节点导纳；

针对感兴趣的每个多相母线和故障类型基于针对所述多相母线确定的所述等效电路确定故障电流注入矢量；并且

针对感兴趣的每个多相母线和故障类型基于所述因式分解的故障前导纳矩阵以及针对所述多相母线和故障类型确定的所述故障电流注入矢量确定故障电压矢量。

2.根据权利要求1所述的方法，其中所述等效电路是戴维南等效电路，每个戴维南等效电路进一步包括N×1电压源矢量，所述N×1电压源矢量具有被确定为在分析中的所述对应母线处的故障前电压的元素。

3.根据权利要求1所述的方法，其中针对感兴趣的每个多相母线确定的所述实验电流注入矢量具有与该母线的每相对应的单位值元素和其他的零值元素。

4.根据权利要求3所述的方法，其中确定所述阻抗矩阵的所述元素包括：针对每个阻抗矩阵，求解所述阻抗矩阵的所述元素为从所述故障前导纳矩阵因式分解的下三角矩阵、从所述故障前导纳矩阵因式分解的上三角矩阵、从所述故障前导纳矩阵因式分解的置换矩阵、从所述故障前导纳矩阵因式分解的列重新排序矩阵以及针对对于其计算所述阻抗矩阵的所述多相母线确定的所述实验电流注入矢量的函数。

5.根据权利要求1所述的方法，其中针对感兴趣的每个多相母线和故障类型确定的所述故障电流注入矢量具有与所述母线的所述相位对应的非零值元素和其他的零值元素。

6.根据权利要求5所述的方法，其中针对感兴趣的每个多相母线和故障类型确定所述故障电压矢量包括，针对感兴趣的每个多相母线和故障类型求解所述故障电压矢量的每个元素为从所述故障前导纳矩阵因式分解的下三角矩阵、从所述故障前导纳矩阵因式分解的上三角矩阵、从所述故障前导纳矩阵因式分解的置换矩阵、从所述故障前导纳矩阵因式分解的列重新排序矩阵以及针对所述多相母线和故障类型确定的所述故障电流注入矢量的函数。

7.根据权利要求1所述的方法，其中针对感兴趣的每个多相母线和故障类型确定所述故障电流注入矢量包括：

针对感兴趣的每个多相母线和故障类型，基于针对所述多相母线确定的所述等效电路计算故障后电流矢量；以及

将每个故障后电流矢量与代表所述电网的电流源节点的复数源电流矢量相加，以针对感兴趣的每个多相母线和故障类型确定所述故障电流矢量。

8.根据权利要求1所述的方法，进一步包括基于所述故障前导纳矩阵和针对感兴趣的每个多相母线确定的所述故障电压矢量，计算所述多相母线上的故障电流。

9.根据权利要求8所述的方法，进一步包括：

测量所述电网中的一个或者多个故障电流；以及

将所述一个或者多个测量的故障电流与所述计算的故障电流比较，以识别所述电网中的一个或者多个可能故障条件。

10.一种存储计算机程序的非瞬态计算机可读介质，所述计算机程序可操作用于针对具有多个多相母线的电网执行故障分析，所述计算机程序包括：

用于针对感兴趣的所述多相母线计算等效电路的程序指令，每个等效电路包括N×N阻抗矩阵，其中N对应于针对其计算所述等效电路的所述多相母线的相位的数目；

用于确定所述电压源矢量的元素为在故障分析之下的每个母线处的故障前电压的程序指令；

用于基于多个实验电流注入矢量所产生的电压以及因式分解的故障前导纳矩阵，确定所述阻抗矩阵的元素的程序指令，所述故障前导纳矩阵代表无故障的所述多相母线的节点导纳；

用于针对感兴趣的每个多相母线和故障类型基于针对所述多相母线确定的所述等效电路来确定故障电流注入矢量的程序指令；以及

用于针对感兴趣的每个多相母线和故障类型基于所述因式分解的故障前导纳矩阵以及针对所述多相母线和故障类型确定的所述故障电流注入矢量来确定故障电压矢量的程序指令。

11.根据权利要求10所述的非瞬态计算机可读介质，进一步包括用以基于所述故障前导纳矩阵和针对感兴趣的每个多相母线确定的所述故障电压矢量计算所述多相母线上的故障电流，并且将所述电网中的一个或者多个测量的故障电流与所述计算的故障电流比较，以识别所述电网中的一个或者多个可能故障条件的程序指令。

12.一种与具有多个多相母线的电网通信的计算机系统，所述计算机系统包括处理电路，所述处理电路可操作用于：

针对感兴趣的所述多相母线计算等效电路，每个等效电路包括N×N阻抗矩阵，其中N对应于针对其计算所述等效电路的所述多相母线的相位的数目；

基于多个实验电流注入矢量所产生的电压以及因式分解的故障前导纳矩阵，确定所述阻抗矩阵的元素，所述故障前导纳矩阵代表无故障的所述多相母线的节点导纳；

针对感兴趣的每个多相母线和故障类型，基于针对所述多相母线确定的所述等效电路来确定故障电流注入矢量；以及

针对感兴趣的每个多相母线和故障类型，基于所述因式分解的故障前导纳矩阵以及针对所述多相母线和故障类型确定的所述故障电流注入矢量，来确定故障电压矢量。

13.根据权利要求12所述的计算机系统，其中所述处理电路可操作用于将所述感兴趣的多相母线建模为戴维南等效电路，每个戴维南等效电路进一步包括N×1电压源矢量，所述N×1电压源矢量具有由所述处理电路确定的，作为在分析之下的所述对应母线处的故障前电压的元素。

14.根据权利要求12所述的计算机系统，其中所述处理电路可操作用以针对感兴趣的每个多相母线确定所述实验电流注入矢量为对应于所述母线的每相的单位值元素和其他的零值元素。

15.根据权利要求14所述的计算机系统，其中对于每个阻抗矩阵，所述处理电路可操作用于求解所述阻抗矩阵的每个元素为从所述故障前导纳矩阵因式分解的下三角矩阵、从所述故障前导纳矩阵因式分解的上三角矩阵、从所述故障前导纳矩阵因式分解的置换矩阵、从所述故障前导纳矩阵因式分解的列重新排序矩阵以及对于针对其计算所述阻抗矩阵的所述多相母线确定的所述实验电流注入矢量的函数。

16.根据权利要求12所述的计算机系统，其中所述处理电路可操作用于针对感兴趣的每个多相母线和故障类型确定所述故障电流注入矢量为与所述母线的相对应的非零值元素和其他的零值元素。

17.根据权利要求16所述的计算机系统，其中对于感兴趣的每个多相母线和故障类型，所述处理电路可操作用于求解所述故障电压矢量的每个元素为从所述故障前导纳矩阵因式分解的下三角矩阵、从所述故障前导纳矩阵因式分解的上三角矩阵、从所述故障前导纳矩阵因式分解的置换矩阵、从所述故障前导纳矩阵因式分解的列重新排序矩阵以及针对该多相母线和故障类型确定的所述故障电流注入矢量的函数。

18.根据权利要求12所述的计算机系统，其中所述处理电路可操作用于为针对感兴趣的每个多相母线和故障类型，基于针对该多相母线确定的所述等效电路，计算故障后电流矢量，并且将每个故障后电流矢量与代表所述电网的电流源节点的复数源电流矢量相加，以针对感兴趣的每个多相母线和故障类型，确定所述故障电流矢量。

19.根据权利要求12所述的计算机系统，其中所述处理电路进一步可操作用于基于所述故障前导纳矩阵和针对感兴趣的每个多相母线确定的所述故障电压矢量，计算所述多相母线处的故障电流。

20.根据权利要求19所述的计算机系统，其中所述处理电路进一步可操作用于将所述电网中的一个或者多个测量的故障电流与所述计算的故障电流比较以识别所述电网中的一个或者多个可能故障条件。