CN102879668B - 含逆变型分布式电源配电网的非对称故障分析方法 - Google Patents

Abstract

本发明公开了含逆变型分布式电源配电网的非对称故障分析方法，本发明结合分布式电源的故障穿越控制策略与输出特性，首先计算配电网正常运行时的公共联结点PCC的电压幅值U_PCC，接着建立配电网故障时的公共联结点PCC的正序电压幅值的求解方程组，求解配电网故障时公共联结点的正序电压，从而实现对含逆变型分布式电源配电网的故障分析。与现有技术相比，本发明提高了故障分析的准确性。

Description

含逆变型分布式电源配电网的非对称故障分析方法

技术领域

本发明涉及电力***故障分析方法，特别涉及含逆变型分布式电源配电网的非对称故障分析方法。

背景技术

对称分量法是非对称故障分析最主要的方法，对传统的同步电机而言，在正序网络中可采用电压源等值，而在负序网络中可采用阻抗等值。然而，逆变型分布式电源与传统的同步电机类型电源完全不同，其输出特性在非对称故障条件下完全由其控制策略决定，在正、负序网络中都不能采用传统的同步电机等值模型。因此，现有的故障分析方法必然会产生较大的误差。为实现含逆变型分布式电源配电网的非对称故障精确分析，必须从逆变型电源的等值模型和配电网故障分析模型上予以改进。

发明内容

为了克服现有技术的上述缺点与不足，本发明针对正序分量控制的逆变型分布式电源，提出了计及故障穿越控制策略与输出特性的配电网非对称故障分析方法，实现了含逆变型分布式电源配电网的非对称故障的精确分析。

本发明的目的通过以下技术方案实现：含逆变型分布式电源配电网的非对称故障分析方法，包括以下步骤(下述步骤中的各变量都表示标幺值)：

S1计算配电网正常运行时的公共联结点（PCC）的电压幅值U_PCC；

S2建立配电网故障时的复合序网，分布式电源只包含在正序网络中，建立配电网故障时的PCC的正序电压幅值的求解方程组，具体包括以下步骤：

S21根据配电网故障时的复合序网得到故障时的节点电压方程：

$\left[Y^{\′}\right]{\stackrel{\·}{U}}_f={\stackrel{\·}{I}}_f---\left(1\right)$

其中， $\left[Y^{\′}\right]=\left[\begin{array}{cccc}{Y^{\′}}_{11}& {Y^{\′}}_{12}& \·\·\·& {Y^{\′}}_{1m}\\ {Y^{\′}}_{21}& {Y^{\′}}_{22}& \·\·\·& {Y^{\′}}_{2m}\\ \·& \·& & \·\\ \·& \·& & \·\\ \·& \·& & \·\\ {Y^{\′}}_{m1}& {Y^{\′}}_{m2}& \·\·\·& {Y^{\′}}_{\mathrm{mm}}\end{array}\right]$ 表示故障时的节点导纳矩阵；对角元素Y_ii为节点i故障时的自导纳，其值等于故障时接于节点i的所有支路导纳之和；非对角元素Y′_ij为节点i、j间故障时的互导纳，当节点i、j间存在支路时，Y′_ij等于直接联接于节点i、j间的支路导纳的负值；当节点i、j间不存在支路时，Y′_ij=0；

${\stackrel{\·}{U}}_f=\left[\begin{array}{c}{\stackrel{\·}{U}}_{1\·f}\\ {\stackrel{\·}{U}}_{2\·f}\\ \·\\ \·\\ \·\\ {\stackrel{\·}{U}}_{\mathrm{PCC}\·f}^{+}\\ \·\\ \·\\ \·\\ {\stackrel{\·}{U}}_{\mathrm{PCC}\·f}^{-}\\ \·\\ \·\\ \·\\ {\stackrel{\·}{U}}_{m\·f}\end{array}\right]$ 表示故障时节点电压，其中分别为故障时PCC的正、负序电压；

${\stackrel{\·}{I}}_f=\left[\begin{array}{c}{\stackrel{\·}{I}}_{1\·f}\\ {\stackrel{\·}{I}}_{2\·f}\\ \·\\ \·\\ \·\\ {\stackrel{\·}{I}}_{\mathrm{DG}\·f}\\ \·\\ \·\\ \·\\ {\stackrel{\·}{I}}_{m\·f}\end{array}\right]$ 表示故障时节点的注入电流，其中为故障时分布式电源(DG)注入PCC的电流；

表示为： ${\stackrel{\·}{I}}_{\mathrm{DG}\·f}=I_{d\·f}-{\mathrm{jI}}_{q\·f}---\left(2\right)$

其中， $\left\{\begin{array}{c}{I}_{d\·f}=P_{\left(0\right)}/U_{\mathrm{PCC}\·f}^{+}\\ I_{q\·f}=I_{q\left(0\right)}+k_q(U_{\mathrm{PCC}}-U_{\mathrm{PCC}\·f}^{+})\end{array}\right.$

式中，I_d·f、I_q·f分别表示故障时分布式电源输出有功电流与无功电流；P₍₀₎、I_q(0)分别表示正常运行时分布式电源输出的有功功率与无功电流，Q₍₀₎表示正常运行时分布式电源输出的无功功率，k_q表示系数；

S22对式（1）进行线性变换，得到的求解方程组；

S3求解

S4求解I_d·f、I_q·f；

S5定义I_ad＝I_d·f+I_q·f，并判断I_ad是否超过逆变器额定电流I_VSC·n，若未超过则进行步骤S6，否则，令I_d·f=I_VSC·n-I_q·f，将I_d·f、I_q·f表达式代入式（1），并重新计算后进行步骤S6；

S6利用式（1）计算中除外的节点电压；

S7根据下式计算节点j与节点k之间的支路电流：

${\stackrel{\·}{I}}_{\mathrm{jk}\·f}=\frac{{\stackrel{\·}{U}}_{i\·f}-{\stackrel{\·}{U}}_{k\·f}}{Z_{\mathrm{jk}}}$

其中，Z_jk分别表示节点j与k之间的支路电流与支路阻抗，分别表示配电网故障时节点j、k的电压。

步骤S1所述计算配电网正常运行时PCC的电压幅值U_PCC，具体包括以下步骤：

S11根据正常运行时配电网等值图得到正常运行时的节点电压方程：

$\left[Y\right]\stackrel{\·}{U}=\stackrel{\·}{I}$

其中， $\left[Y\right]=\left[\begin{array}{cccc}{Y}_{11}& Y_{12}& \·\·\·& Y_{1n}\\ Y_{21}& Y_{22}& \·\·\·& Y_{2n}\\ \·& \·& & \·\\ \·& \·& & \·\\ \·& \·& & \·\\ Y_{n1}& Y_{n2}& \·\·\·& Y_{\mathrm{nn}}\end{array}\right]$ 表示正常运行时的节点导纳矩阵；对角元素Y_ij为节点i正常运行时的自导纳，其值等于正常运行时接于节点i的所有支路导纳之和；非对角元素Y_ij为节点i、j间正常运行时的互导纳，当节点i、j间存在支路时，Y_ij等于直接联接于节点i、j间的支路导纳的负值；当节点i、j间不存在支路时，Y_ij=0；

$\stackrel{\·}{U}=\left[\begin{array}{c}{\stackrel{\·}{U}}_1\\ {\stackrel{\·}{U}}_2\\ \·\\ \·\\ \·\\ {\stackrel{\·}{U}}_{\mathrm{PCC}}\\ \·\\ \·\\ \·\\ {\stackrel{\·}{U}}_n\end{array}\right]$ 表示正常运行时的节点电压，为正常运行时PCC的电压；

$\stackrel{\·}{I}=\left[\begin{array}{c}{\stackrel{\·}{I}}_1\\ {\stackrel{\·}{I}}_2\\ \·\\ \·\\ \·\\ {\stackrel{\·}{I}}_{\mathrm{DG}}\\ \·\\ \·\\ \·\\ {\stackrel{\·}{I}}_n\end{array}\right]$ 表示正常运行时节点的注入电流，为正常运行时DG注入PCC节点的电流；

S12根据节点电压方程求得然后取绝对值得到U_PCC。

与现有技术相比，本发明具有以下优点和有益效果：本发明通过计及逆变型分布式电源控制策略与输出特性，建立新的分布式电源暂态等值模型，能更加真实地反映分布式电源的故障电流特性；在此基础上建立新的配电网故障分析模型，提高了故障分析的准确性。本方法为含逆变型分布式电源配电网的设备选型与保护整定等方面提供了科学的依据，在工程实践中具有很强的实用性。

附图说明

图1为本发明的实施例的配电网单线图。

图2为本发明的实施例的含逆变型分布式电源配电网的非对称故障分析方法的流程图。

图3为本发明的实施例中图1所示的配电网故障时的复合序网图。

具体实施方式

下面结合实施例及附图，对本发明作进一步地详细说明，但本发明的实施方式不限于此。

实施例

本实施例以图1所示的配电网为例，采用本发明含逆变型分布式电源配电网的非对称故障分析方法进行故障分析，如图2所示，包括以下步骤：

S1计算配电网正常运行时PCC的电压幅值U_PCC，具体包括以下步骤：

S11根据正常运行时配电网等值图得到正常运行时的节点电压方程：

$\left[Y\right]\stackrel{\·}{U}=\stackrel{\·}{I}$

本实施例中正常运行时的节点电压方程为：

$\left[\begin{array}{cc}\frac{1}{Z_s}+\frac{1}{Z_{L1}}& -\frac{1}{L_1}\\ -\frac{1}{Z_{L1}}& \frac{1}{Z_{L1}}+\frac{1}{Z_{L2}}\end{array}\right]\left[\begin{array}{c}{\stackrel{\·}{U}}_M\\ {\stackrel{\·}{U}}_{\mathrm{PCC}}\end{array}\right]=\left[\begin{array}{c}\frac{{\stackrel{\·}{E}}_s}{Z_s}\\ {\stackrel{\·}{I}}_{\mathrm{DG}}\end{array}\right]$

其中分别为正常运行时节点M、PCC的自导纳，为节点M、PCC之间的互导纳；Z_s、Z_L1、Z_L2分别表示***的等值阻抗、PCC上游线路L1阻抗及PCC下游线路L2阻抗；表示***等值电势，取值为分布式电源额定电流；为正常运行时节点M的电压；为正常运行时PCC的电压；为注入M点的电流。

S12根据节点电压方程求得然后取绝对值得到U_PCC；

S2建立配电网故障时的复合序网（如图3所示），分布式电源只包含在正序网络中，建立电网故障时的PCC的正序电压的求解方程组，具体步骤如下：

S21根据配电网故障时的复合序网得到故障时的节点电压方程：

$\left[Y^{\′}\right]{\stackrel{\·}{U}}_f={\stackrel{\·}{I}}_f---\left(1\right)$

其中， $\left[Y^{\′}\right]=\left[\begin{array}{cccc}{Y^{\′}}_{11}& {Y^{\′}}_{12}& \·\·\·& {Y^{\′}}_{1m}\\ {Y^{\′}}_{21}& {Y^{\′}}_{22}& \·\·\·& {Y^{\′}}_{2m}\\ \·& \·& & \·\\ \·& \·& & \·\\ \·& \·& & \·\\ {Y^{\′}}_{m1}& {Y^{\′}}_{m2}& \·\·\·& {Y^{\′}}_{\mathrm{mm}}\end{array}\right]$ 表示故障时节点的导纳矩阵；对角元素Y_ii为节点i故障时的自导纳，其值等于故障时接于节点i的所有支路导纳之和；非对角元素Y′_ij为节点i、j间故障时的互导纳，当节点i、j间存在支路时，Y′_ij等于直接联接于节点i、j间的支路导纳的负值，当节点i、j间不存在支路时，Y′_ij=0；

即 $\left[\begin{array}{ccc}\frac{1}{Z_s}+\frac{1}{Z_{L1}}& -\frac{1}{Z_{L1}}& 0\\ -\frac{1}{Z_{L1}}& \frac{1}{Z_{L1}}+\frac{1}{2\mathrm{\β}{Z}_{L2}}& -\frac{1}{2\mathrm{\β}{Z}_{L2}}\\ 0& -\frac{1}{2\mathrm{\β}{Z}_{L2}}& \frac{1}{2\mathrm{\β}{Z}_{L2}}+\frac{1}{Z_s}+\frac{1}{Z_{L1}}\end{array}\right]\left[\begin{array}{c}{\stackrel{\·}{U}}_{M\·f}\\ {\stackrel{\·}{U}}_{\mathrm{PCC}\·f}^{+}\\ {\stackrel{\·}{U}}_{\mathrm{PCC}\·f}^{-}\end{array}\right]=\left[\begin{array}{c}\frac{{\stackrel{\·}{E}}_s}{Z_s}\\ {\stackrel{\·}{I}}_{\mathrm{DG}\·f}\\ 0\end{array}\right]$

其中分别为故障时节点M、PCC正序电压节点、PCC负序电压节点的自导纳；为故障时节点M、PCC正序电压节点之间的互导纳；，表示PCC正序电压节点与PCC负序电压节点之间的互导纳，β为表示在线路L2中的故障位置，取值为0~100%；分别为故障时PCC的正、负序电压；为故障时DG注入PCC的电流；根据分布式电源控制策略，表示为：

${\stackrel{\·}{I}}_{\mathrm{DG}\·f}=I_{d\·f}-{\mathrm{jI}}_{q\·f}---\left(2\right)$

其中， $\left\{\begin{array}{c}{I}_{d\·f}=P_{\left(0\right)}/U_{\mathrm{PCC}\·f}^{+}\\ I_{q\·f}=I_{q\left(0\right)}+k_q(U_{\mathrm{PCC}}-U_{\mathrm{PCC}\·f}^{+})\end{array}\right.$

式中，I_d·f、I_q·f分别表示故障时分布式电源输出的故障电流、有功电流与无功电流；P₍₀₎、I_q(0)分别表示正常运行时分布式电源输出的有功功率与无功电流，Q₍₀₎表示正常运行时分布式电源输出的无功功率，k_q表示系数；

S22对式（1）进行线性变换，得到的求解方程组，具体为：

以电压相位为基准，则 ${\stackrel{\·}{E}}_s=E_s(\cos \mathrm{\α}+j\sin \mathrm{\α}),$ 得到的求解方程组：

$\left\{\begin{array}{c}a{\left(U_{\mathrm{PCC}\·f}^{+}\right)}^2+e{U}_{\mathrm{PCC}\·f}^{+}-c=E_s{U}_{\mathrm{PCC}\·f}^{+}\cos \mathrm{\α}\\ b{\left(U_{\mathrm{PCC}\·f}^{+}\right)}^2+f{U}_{\mathrm{PCC}\·f}^{+}-d=E_s{U}_{\mathrm{PCC}\·f}^{+}\sin \mathrm{\α}\end{array}\right.$

式中，a+jb=[1/(2βZ_L2+Z_s+Z_L1)+1/(Z_s+_L1)-jk_q/U_n](Z_s+Z_L1)，c+jd=U_PCCI_d(0)(Z_s+Z_L1)，e+jf=j(I_q(0)+k_qU_PCC/U_n)(Z_s+Z_L1)；

S3求解

S4求解I_d·f、I_q·f；

S5定义I_ad＝I_d·f+I_q·f，并判断I_ad是否超过逆变器额定电流I_VAC·n，若未超过则进行步骤S6，否则，令I_d·f=I_VSC·n-I_q·f，将I_d·f、I_q·f表达式代入式（1）及式（2），并重新计算后进行步骤S6；

其中，重新计算的过程如下：

将I_d·f、I_q·f表达式代入式（1）及式（2），得到的求解方程组：

$\left\{\begin{array}{c}{U}_{\mathrm{PCC}\·f}^{+}(a-K_q/U_n)c{E}_s\cos \mathrm{\α}+d{E}_s\sin \mathrm{\α}=I_{\max }-(I_{q\left(0\right)}+K_q{U}_{\mathrm{PCC}\·f}^{+}/U_n)\\ U_{\mathrm{PCC}\·f}^{+}(b-K_q/U_n)-d{E}_s\cos \mathrm{\α}-c{E}_s\sin \mathrm{\α}=-(I_{q\left(0\right)}{K}_q{U}_{\mathrm{PCC}}/U_n)\end{array}\right.$

式中，a+jb=1/(2βZ_L2+Z_s+Z_L1)+1/(Z_s+Z_L1)，c+jd=1/(Z_s+Z_L1)；S6利用式（1）计算中除外的节点电压；

S7根据下式计算节点j与k之间的支路电流：

${\stackrel{\·}{I}}_{\mathrm{jk}\·f}=\frac{{\stackrel{\·}{U}}_{i\·f}-{\stackrel{\·}{U}}_{k\·f}}{Z_{\mathrm{jk}}}$

其中，Z_jk分别表示节点j与k之间的支路电流与支路阻抗，分别表示配电网故障时节点j、k的电压。

本实施例的图1中，配电网***的等值阻抗Z_s为1.3j(Ω)，线路L1的等值阻抗Z_L1、L2等值阻抗Z_L2分别为1.18+3.56j(Ω)、0.59+1.78j(Ω)，分布式电源额定容量与逆变器接口额定容量分别为4MW、8MVA，k_q为2。由上述条件可得到I_VSC·n＝0.46kA。故障前分布式电源有功输出为额定值，无功输出为零。

下面列举两种不同短路条件予以说明：

情况1：

两相短路故障发生在L2线路末端，即β=100%，进行步骤S1~S4，计算得到的值为7.26kV，U_PCC为10.1kV。由U_PCC、可得I_d·f、I_q·f、I_DG·f分别为0.32kA、0.13kA、0.35kA。由于I_d·f+I_q·f<I_VSC·n，将I_d·f、I_q·f代入配电网故障时节点电压方程可求得为4.10kV，分别为0.365kA、0.48kA， 都为0.48kA。

情况2：

两相短路故障发生在L2线路70%处，即β=70%，进行步骤S1~S4，并得到求解方程组，计算的值为6.92kV，U_PCC为10.1kV。由U_PCC、可得I_d·f、I_q·f分别为0.33kA、0.14kA。由于I_d·f+I_q·f>I_VSC·n，重新计算I_d·f、I_DG·f为0.32kA、0.35kA。最后求得为4.10kV，分别为0.39kA、0.52kA，都为0.52Ka。

上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受所述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。

Claims (1)

1.含逆变型分布式电源配电网的非对称故障分析方法，其特征在于，包括以下步骤：

S1计算配电网正常运行时的公共联结点PCC的电压幅值U_PCC,具体为：

S11根据正常运行时配电网等值图得到正常运行时的节点电压方程：

$\left[Y\right]\stackrel{.}{U}=\stackrel{.}{I}$

其中， $\left[Y\right]=\left[\begin{array}{cccc}{Y}_{11}& Y_{12}& ...& Y_{1n}\\ Y_{21}& Y_{22}& ...& Y_{2n}\\ \&CenterDot;& \&CenterDot;& & \&CenterDot;\\ \&CenterDot;& \&CenterDot;& & \&CenterDot;\\ \&CenterDot;& \&CenterDot;& & \&CenterDot;\\ Y_{n1}& Y_{n2}& ...& Y_{\mathrm{nn}}\end{array}\right],$ 表示正常运行时的节点导纳矩阵；对角元素Y_ii为节点i正常运行时的自导纳，其值等于正常运行时接于节点i的所有支路导纳之和；非对角元素Y_ij为节点i、j间正常运行时的互导纳，当节点i、j间存在支路时，Y_ij等于直接联接于节点i、j间的支路导纳的负值；当节点i、j间不存在支路时，Y_ij＝0；

$\stackrel{\&CenterDot;}{U}=\left[\begin{array}{c}{\stackrel{\&CenterDot;}{U}}_1\\ {\stackrel{\&CenterDot;}{U}}_2\\ \&CenterDot;\\ \&CenterDot;\\ \&CenterDot;\\ {\stackrel{\&CenterDot;}{U}}_{\mathrm{PCC}}\\ \&CenterDot;\\ \&CenterDot;\\ \&CenterDot;\\ {\stackrel{\&CenterDot;}{U}}_n\end{array}\right]$ 表示正常运行时的节点电压，为正常运行时公共联结点PCC的电压；

$\stackrel{\&CenterDot;}{I}=\left[\begin{array}{c}{\stackrel{\&CenterDot;}{I}}_1\\ {\stackrel{\&CenterDot;}{I}}_2\\ \&CenterDot;\\ \&CenterDot;\\ \&CenterDot;\\ {\stackrel{\&CenterDot;}{I}}_{\mathrm{DG}}\\ \&CenterDot;\\ \&CenterDot;\\ \&CenterDot;\\ {\stackrel{\&CenterDot;}{I}}_n\end{array}\right]$ 表示正常运行时节点的注入电流，为正常运行时分布式电源注入公共联结点PCC节点的电流；

S12根据节点电压方程求得然后取绝对值得到U_PCC；

S2建立配电网故障时的复合序网，分布式电源只包含在正序网络中，建立配电网故障时的公共联结点PCC的正序电压幅值的求解方程组，具体包括以下步骤：

S21根据配电网故障时的复合序网得到故障时的节点电压方程：

$\left[Y^{\&prime;}\right]{\stackrel{\&CenterDot;}{U}}_f={\stackrel{\&CenterDot;}{I}}_f---\left(1\right)$

其中， $\left[Y^{\&prime;}\right]=\left[\begin{array}{cccc}{Y^{\&prime;}}_{11}& {Y^{\&prime;}}_{12}& ...& {Y^{\&prime;}}_{1m}\\ {Y^{\&prime;}}_{21}& {Y^{\&prime;}}_{22}& ...& {Y^{\&prime;}}_{2m}\\ \&CenterDot;& \&CenterDot;& & \&CenterDot;\\ \&CenterDot;& \&CenterDot;& & \&CenterDot;\\ \&CenterDot;& \&CenterDot;& & \&CenterDot;\\ {Y^{\&prime;}}_{m1}& {Y^{\&prime;}}_{m2}& ...& {Y^{\&prime;}}_{\mathrm{mm}}\end{array}\right]$ 表示故障时的节点导纳矩阵；对角元素Y'_ii为节点i故障时的自导纳，其值等于故障时接于节点i的所有支路导纳之和；非对角元素Y'_ij为节点i、j间故障时的互导纳，当节点i、j间存在支路时，Y'_ij等于直接联接于节点i、j间的支路导纳的负值，当节点i、j间不存在支路时，Y'_ij＝0；

${\stackrel{\&CenterDot;}{U}}_f=\left[\begin{array}{c}{\stackrel{\&CenterDot;}{U}}_{1\&CenterDot;f}\\ {\stackrel{\&CenterDot;}{U}}_{2\&CenterDot;f}\\ \&CenterDot;\\ \&CenterDot;\\ \&CenterDot;\\ {\stackrel{\&CenterDot;}{U}}_{\mathrm{PCC}\&CenterDot;f}^{+}\\ .\\ \&CenterDot;\\ \&CenterDot;\\ {\stackrel{\&CenterDot;}{U}}_{\mathrm{PCC}\&CenterDot;f}^{-}\\ \&CenterDot;\\ \&CenterDot;\\ \&CenterDot;\\ {\stackrel{\&CenterDot;}{U}}_{m\&CenterDot;f}\end{array}\right]$ 表示故障时节点电压，其中分别为故障时公共联结点PCC的正、负序电压；

${\stackrel{\&CenterDot;}{I}}_f=\left[\begin{array}{c}{\stackrel{\&CenterDot;}{I}}_{1\&CenterDot;f}\\ {\stackrel{\&CenterDot;}{I}}_{2\&CenterDot;f}\\ \&CenterDot;\\ \&CenterDot;\\ \&CenterDot;\\ {\stackrel{\&CenterDot;}{I}}_{\mathrm{DG}\&CenterDot;f}\\ \&CenterDot;\\ \&CenterDot;\\ \&CenterDot;\\ {\stackrel{\&CenterDot;}{I}}_{m\&CenterDot;f}\end{array}\right]$ 表示故障时节点的注入电流，其中为故障时分布式电源注入公共联结点PCC的电流；

表示为： ${\stackrel{.}{I}}_{\mathrm{DG}\&CenterDot;f}=I_{d\&CenterDot;f}-{\mathrm{jI}}_{q\&CenterDot;f}---\left(2\right)$

其中， $\left\{\begin{array}{c}{I}_{d\&CenterDot;f}=P_{\left(0\right)}/U_{\mathrm{PCC}\&CenterDot;f}^{+}\\ I_{q\&CenterDot;f}=I_{q\left(0\right)}+k_q(U_{\mathrm{PCC}}-U_{\mathrm{PCC}\&CenterDot;f}^{+})\end{array}\right.$

式中，I_d·f、I_q·f分别表示故障时分布式电源输出的有功电流与无功电流；P₍₀₎、I_q(0)分别表示正常运行时分布式电源输出的有功功率与无功电流，Q₍₀₎表示正常运行时分布式电源输出的无功功率，k_q表示系数；

S22对式(1)进行线性变换，得到的求解方程组；

S3求解

S4求解I_d·f、I_q·f；

S5定义I_ad＝I_d·f+I_q·f，并判断I_ad是否超过逆变器额定电流I_VSC·n，若未超过则进行步骤S6，否则，令I_d·f＝I_VSC·n-I_q·f，将I_d·f、I_q·f表达式代入式(1)及式(2)，并重新计算后进行步骤S6；

S6利用式(1)计算中除外的节点电压；

S7根据下式计算节点j与节点k之间的支路电流：

${\stackrel{\&CenterDot;}{I}}_{\mathrm{jk}\&CenterDot;f}=\frac{{\stackrel{\&CenterDot;}{U}}_{j\&CenterDot;f}-{\stackrel{\&CenterDot;}{U}}_{k\&CenterDot;f}}{Z_{\mathrm{jk}}}$

其中，分别表示节点j与k之间的支路电流与支路阻抗，分别表示配电网故障时节点j、k的电压。