CN103972886A - 一种主动配电网谐波潮流分析方法 - Google Patents

Abstract

本发明公开了一种主动配电网谐波潮流分析方法，包含以下内容，首先，根据主动配电网中非线性负荷和分布式电源(DG)的特点，分析非线性负荷和DG并网接口电力电子装置的工作特性及其注入系统各次谐波电流的含量；然后，通过计算主动配电网基波潮流得出非线性负荷和DG并网的基波电压，由非线性负荷和DG的功率，可计算出其基波电流；最后，通过非线性负荷和DG的基波电流求出其注入主动配电网的各次谐波电流，并计算出整个主动配电网谐波潮流分布情况。本发明无需迭代计算即可得出所有节点的各次谐波电压分布情况，从根本上改善了谐波潮流分析方法的计算速度和效率，克服了谐波潮流迭代收敛性的问题，具有计算过程清晰，易于编程等特点。

Description

一种主动配电网谐波潮流分析方法

技术领域

本发明属于电力系统分析与计算领域，具体涉及一种主动配电网谐波潮流分析方法。

背景技术

谐波潮流计算是谐波研究领域的重要分支，是了解主动配电网谐波特性和进行谐波分析的重要手段，可计算出关键节点的谐波指标，进而提出抑制谐波的有效措施。随着主动配电网(ADN)技术迅速发展，越来越多的分布式电源(DG)接入主动配电网，以及随着科学技术水平和生活水平不断提高，主动配电网中接入了大量的非线性负荷，而这些非线性负荷和带有大量电力电子器件的分布式源必将给主动配电网的电能质量带来新的问题和挑战。毫无疑问，在电力公司和工业生产中，从规划到设计阶段，主动配电网谐波分析是在分布式电源、非线性负荷不断增加的前提下，保证主动配电网和设备可靠运行并预测潜在问题的有效手段。研究DG和非线性负荷对主动配电网谐波分布的不利影响，有助于找到解决措施；分析DG对谐波分布的有利影响，有利于充分发挥其优势来改善电能质量，因此，研究DG和非线性负荷对主动配电网谐波和电能质量的影响，对主动配电网的发展具有重要意义。本发明从谐波潮流分析方法着手，提出一种主动配电网谐波潮流分析方法，给研究主动配电网谐波影响提供一种有效的技术方法。

发明内容

发明目的：针对上述现有技术存在的问题和不足，本发明的目的是提供一种主动配电网谐波潮流分析方法，该方法采用非迭代式谐波潮流分析方法，无需迭代计算即可得出主动配电网中所有节点的各次谐波电压分布情况，从根本上改善了谐波潮流分析方法的计算速度和效率，克服了谐波潮流迭代收敛性的问题，较大的提高了谐波潮流的计算速度和效率，具有较高的计算精度，计算过程清晰，易于编程等特点。

技术方案：为实现上述发明目的，本发明采用的技术方案为一种主动配电网谐波潮流分析方法，该方法包括以下步骤：

步骤A、获取主动配电网的网络参数，包括系统总节点数，总支路数，节点序号，支路序号，独立节点数，独立支路数，关联支路数，环网个数，参考节点，平衡节点，支路阻抗，线性负荷功率，非线性负荷功率，线性负荷节点序号，非线性负荷节点序号，网络拓扑结构，DG并网个数，DG输出功率，以及DG并网的节点序号，如针对具有N个节点、K个环网、W个非线性负荷和M个DG的主动配电网，假设首节点作为参考节点和平衡节点，则主动配电网独立节点数为n＝N-1，独立支路数即树支数为c＝n，关联支路数即连支数为d＝K，总支路数为b＝c+d，其中DG为分布式电源的英文缩写且代表分布式电源；

步骤B、分析非线性负荷运行特性，通过傅氏分析公式计算出其各次谐波电流分量与其基波电流的比值C(h)，其中，C(h)为非线性负荷各谐波电流分量与其基波电流的比值，h为谐波次数；

步骤C、分析DG并网接口电力电子装置的工作特性，通过傅氏分析公式计算出其注入主动配电网各次谐波电流分量与其基波电流的比值ρ(h)，其中，ρ(h)为DG注入主动配电网第h次谐波电流分量与其基波电流的比值，h为谐波次数；

步骤D、通过常规潮流算法计算主动配电网基波潮流，得出各非线性负荷节点的基波电压和各DG并网节点的基波电压其中，为非线性负荷节点的基波电压，为DG并网节点的基波电压，下标“NL”表示非线性负荷，下标“DG”表示分布式电源，下标“i”表示非线性负荷节点和DG并网节点各自对应的节点序号，上标“(1)”表示基波分量；

步骤E、根据主动配电网中非线性负荷的有功功率P_NL,i和无功功率Q_NL,i，计算出非线性负荷基波电流为其中，为非线性负荷基波电流，为非线性负荷节点的基波电压，P_NL,i和Q_NL,i分别为非线性负荷的有功功率和无功功率，上标“*”表示取共轭复数，j为复数的虚部单位；

步骤F、根据DG输出的有功功率P_DG,i和无功功率Q_DG,i，计算出DG注入主动配电网的基波电流为其中，为DG注入主动配电网的基波电流，为DG并网节点的基波电压，P_DG,i和Q_DG,i分别为DG输出的有功功率和无功功率，上标“*”表示取共轭复数，j为复数的虚部单位；

步骤G、根据非线性负荷基波电流计算出该非线性负荷的第h次谐波电流为其中，为非线性负荷基波电流，为非线性负荷的第h次谐波电流，h为谐波次数；

步骤H、根据DG注入主动配电网的基波电流计算出该DG注入主动配电网的第h次谐波电流为其中，为DG注入主动配电网的基波电流，为DG注入主动配电网的第h次谐波电流，h为谐波次数；

步骤I、根据公式计算主动配电网中所有谐波源节点向主动配电网输入的各次谐波电流其中，为谐波源节点i向主动配电网输入的第h次谐波电流，为非线性负荷的第h次谐波电流，为DG注入主动配电网的第h次谐波电流，h为谐波次数；

步骤J、根据公式计算主动配电网第h次谐波潮流，得出各独立节点i的第h次谐波电压其中，为独立节点i的第h次谐波电压，V^(h)为主动配电网中各独立节点i的第h次谐波电压组成的n×1阶相量矩阵，即i＝1,2,…,n，I^(h)为主动配电网中各谐波源节点向主动配电网输入的第h次谐波电流组成的n×1阶相量矩阵，即i＝1,2,…,n，T为主动配电网的道路矩阵，T为n×b阶矩阵，T^T为矩阵T的转置矩阵，为支路l的第h次谐波阻抗，为基于对角元素分别为各支路l的第h次谐波阻抗组成的b×b阶对角矩阵，即l＝1,2,…,b，B为主动配电网的回-支关联矩阵，B为d×b阶矩阵，BT为矩阵B的转置矩阵，为矩阵的逆矩阵，即为主动配电网的第h次谐波回路阻抗矩阵，为d×d阶矩阵，且有n为主动配电网独立节点数，b为主动配电网总支路数，d为主动配电网连支数，h为谐波次数；

步骤K、根据不同谐波次数h，重复步骤G、步骤H、步骤I和步骤J即可求出主动配电网所有节点的各次谐波电压分布情况，以及计算出各节点的电压总谐波畸变率和电压单个谐波畸变率。

有益效果：本发明针对现有谐波潮流分析方法存在的不足，提出了一种主动配电网谐波潮流分析方法，该方法根据主动配电网中非线性负荷和分布式电源(DG)的特点，分析非线性负荷和DG并网接口电力电子装置的工作特性及其注入系统各次谐波电流的含量；然后，通过计算主动配电网基波潮流得出非线性负荷和DG并网的基波电压，由非线性负荷和DG的功率，可计算出其基波电流；最后，通过非线性负荷和DG的基波电流求出其注入主动配电网的各次谐波电流，并计算出整个主动配电网谐波潮流分布情况。由于采用非迭代式谐波潮流分析方法，无需迭代计算即可得出主动配电网中所有节点的各次谐波电压分布情况，从根本上改善了谐波潮流分析方法的计算速度和效率，克服了多谐波源系统谐波潮流迭代收敛困难的问题，极大提高了谐波潮流的计算速度和效率，整个方法的计算过程清晰，编程简单，计算速度快，具有较高的计算精度和效率。通过测试算例验证了本发明的有效性和正确性。由此可见，本发明具有很好的工程应用价值和借鉴意义。

附图说明

图1为本发明的总体流程图；

图2为33母线主动配电网示意图。

具体实施方式

下面结合附图和具体实施例，进一步阐明本发明，应理解这些实施例仅用于说明本发明而不用于限制本发明的范围，在阅读了本发明之后，本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。

图1为本发明的总体流程图，具体方法如下：

针对具有N个节点、K个环网、W个非线性负荷和M个分布式电源(DG)的主动配电网，假设首节点作为参考节点和平衡节点，则主动配电网独立节点数为n＝N-1，独立支路数即树支数为c＝n，关联支路数即连支数为d＝K，总支路数为b＝c+d。

对于连通图中一颗选定的树，由于基本回路中仅包含一条连支，基本回路数等于连支数，基本回路-支路(下面简称“回-支”)关联矩阵B描述基本回路、树支、连支之间的联系。其中回-支关联矩阵B是一个d×b阶矩阵，假定连支支路的正方向都是从大号节点指向小号节点，基本回路的正方向与连支支路的正方向相同，如果支路r在回路j上，且二者方向相同，则B(j,r)＝1，如果支路r在回路j上，且二者方向相反，则B(j,r)＝-1，如果支路r不在回路j上，则B(j,r)＝0。

一个节点的道路是指节点沿树到根所经过的路径上的支路集合，节点的道路强调的是路径上的支路，对于一个给定的树，节点的道路是唯一的、只由树支组成，可用道路矩阵T描述。其中道路矩阵T是一个n×b阶矩阵，假定道路的正方向都是从电源点(即根节点)指向各节点，各树支支路正方向与道路正方向相同，如果支路r在道路i上，则T(i,r)＝1，反之T(i,r)＝0。道路矩阵T是一个稀疏矩阵，利用稀疏技术可以降低内存需求。

1)获取主动配电网的网络参数，包括系统总节点数，总支路数，节点序号，支路序号，独立节点数，独立支路数，关联支路数，环网个数，参考节点，平衡节点，支路阻抗，线性负荷功率，非线性负荷功率，线性负荷节点序号，非线性负荷节点序号，网络拓扑结构，DG并网个数，DG输出功率，以及DG并网的节点序号。

2)分析非线性负荷运行特性，通过傅氏分析公式计算出其各次谐波电流分量与其基波电流的比值C(h)，其中，C(h)为非线性负荷各谐波电流分量与其基波电流的比值，h为谐波次数。

3)分析DG并网接口电力电子装置的工作特性，通过傅氏分析公式计算出其注入主动配电网各次谐波电流分量与其基波电流的比值ρ(h)，其中，ρ(h)为DG注入主动配电网第h次谐波电流分量与其基波电流的比值，h为谐波次数。

4)通过常规潮流算法计算主动配电网基波潮流，得出各非线性负荷节点的基波电压和各DG并网节点的基波电压其中，为非线性负荷节点的基波电压，为DG并网节点的基波电压，下标“NL”表示非线性负荷，下标“DG”表示分布式电源，下标“i”表示非线性负荷节点和DG并网节点各自对应的节点序号，上标“(1)”表示基波分量。

步骤4)中的主动配电网基波潮流和各非线性负荷节点和DG并网节点的基波电压计算过程如下：

在主动配电网中，设k为迭代次数，为第k次迭代时独立节点i注入基波电流，为第k次迭代时回路f基波电流(也即连支f基波电流)，为第k次迭代时支路l的基波电流，为第k次迭代时支路l的基波电压，为第k次迭代时各独立节点i注入基波电流组成的n×1阶相量矩阵，即i＝1,2,…,n，n为主动配电网独立节点数，为第k次迭代时各连支f基波电流(也即各回路f基波电流)组成的d×1阶相量矩阵，即f＝1,2,…,d，d为主动配电网连支数，为第k次迭代时各支路l基波电流组成的b×1阶相量矩阵，即l＝1,2,…,b，为第k次迭代时各支路l基波电压组成的b×1阶相量矩阵，即l＝1,2,…,b，b为主动配电网总支路数。

根据KCL电流定理，支路基波电流相量矩阵与独立节点注入基波电流相量矩阵回路基波电流相量矩阵满足如下等式：

$I_b^{\left(1\right)}=T^T{I}_{n_k}^{\left(1\right)}+B^T{I}_{d_k}^{\left(1\right)}---\left(1\right)$

其中，T为分布式发电系统的n×b阶道路矩阵，T^T为矩阵T的转置矩阵。

基于欧姆定理，支路基波电流相量矩阵和支路基波电压相量矩阵满足如下等式：

$V_{b_k}^{\left(1\right)}=Z_b^{\left(1\right)}{I}_{b_k}^{\left(1\right)}---\left(2\right)$

其中，为支路l的基波阻抗，为基于对角元素分别为各支路l的基波阻抗形成的b×b阶对角矩阵，即l＝1,2,…,b，b为主动配电网总支路数。

根据KVL电压定律，基本回路约束为

${\mathrm{Bv}}_{b_k}^{\left(1\right)}=0---\left(3\right)$

由式(1)和式(2)代入式(3)可得：

${\mathrm{BZ}}_b^{\left(1\right)}{B}^T{I}_{d_k}^{\left(1\right)}+{\mathrm{BZ}}_b^{\left(1\right)}{T}^T{I}_{n_k}^{\left(1\right)}=0---\left(4\right)$

令为回路基波阻抗矩阵，其逆矩阵为则有

$I_{d_k}^{\left(1\right)}=-Y_d^{\left(1\right)}{\mathrm{BZ}}_b^{\left(1\right)}{T}^T{I}_{n_k}^{\left(1\right)}---\left(5\right)$

$I_{b_k}^{\left(1\right)}=T^T{I}_{n_k}^{\left(1\right)}-B^T{Y}_d^{\left(1\right)}{\mathrm{BZ}}_b^{\left(1\right)}{T}^T{I}_{n_k}^{\left(1\right)}---\left(6\right)$

$V_{b_k}^{\left(1\right)}=Z_b^{\left(1\right)}{I}_{b_k}^{\left(1\right)}=(Z_b^{\left(1\right)}-Z_b^{\left(1\right)}{B}^T{Y}_b^{\left(1\right)})T^T{I}_{n_k}^{\left(1\right)}---\left(7\right)$

在主动配电网的模型中，可知任一独立节点与首节点的基波电压差等于从此独立节点开始沿着该独立节点的道路到达首节点所经独立支路的支路基波电压之和，即：

$\mathrm{\Δ}{V}_{n_k}^{\left(1\right)}={\mathrm{TV}}_{b_k}^{\left(1\right)}=T(Z_b^{\left(1\right)}-Z_b^{\left(1\right)}{B}^T{Y}_d^{\left(1\right)}{\mathrm{BZ}}_b^{\left(1\right)})T^T{I}_{n_k}^{\left(1\right)}---\left(8\right)$

其中为各独立节点与首节点的基波电压差矩阵，为n×1阶矩阵。

由每一独立节点与首节点的基波电压差，可根据下式(9)求出主动配电网各独立节点基波电压

$V_{n_k}^{\left(1\right)}={\mathrm{\λ}}_n{\stackrel{\·}{V}}_0^{\left(1\right)}-\mathrm{\Δ}{V}_{n_k}^{\left(1\right)}---\left(9\right)$

其中，为首节点的基波电压，λ_n＝[1,1，…,1]^T为n维列向量，为第k次迭代时独立节点i的基波电压，为第k次迭代时各独立节点基波电压组成的n×1阶相量矩阵，即i＝1,2,…,n，n为主动配电网独立节点数。

给主动配电网各独立节点基波电压赋初始值为并令迭代次数k＝1，由下式(10)计算出第k次迭代时各独立节点注入基波电流

${\stackrel{\·}{I}}_{n,i_k}^{\left(1\right)}={({\stackrel{\‾}{S}}_i/{\stackrel{\·}{V}}_{n,i_{(k-1)}}^{\left(1\right)})}^{*}-Y_i{\stackrel{\·}{V}}_{n,i_{(k-1)}}^{\left(1\right)}-{[(P_{\mathrm{DG},i}+{\mathrm{jQ}}_{\mathrm{DG},i})/{\stackrel{\·}{V}}_{n,i_{(k-1)}}^{\left(1\right)}]}^{*}---\left(10\right)$

其中，为独立节点i注入复功率，Y_i为独立节点i并联导纳之和，为第k-1次迭代后求出的独立节点i基波电压，上标“*”表示取共轭复数，P_DG,i和Q_DG,i分别为DG输出的有功功率和无功功率，下标“DG”表示分布式电源，下标“i”表示该DG并网的节点序号。

由式(10)、(8)和(9)重复迭代计算，直到满足收敛精度要求，则停止迭代，输出基波潮流结果，得出主动配电网各独立节点基波电压接着，可求出各非线性负荷节点的基波电压以及各DG并网节点的基波电压其中ε一般取10^-4～10^-6，i＝1,2,…,n，n为主动配电网独立节点数。

5)根据主动配电网中非线性负荷的有功功率P_NL,i和无功功率Q_NL,i，计算出非线性负荷基波电流为其中，为非线性负荷基波电流，为非线性负荷节点的基波电压，P_NL,i和Q_NL,i分别为非线性负荷的有功功率和无功功率，上标“*”表示取共轭复数，j为复数的虚部单位。

6)根据DG输出的有功功率P_DG,i和无功功率Q_DG,i，计算出DG注入主动配电网的基波电流为其中，为DG注入主动配电网的基波电流，为DG并网节点的基波电压，P_DG,i和Q_DG,i分别为DG输出的有功功率和无功功率，上标“*”表示取共轭复数，j为复数的虚部单位。

7)根据非线性负荷基波电流计算出该非线性负荷的第h次谐波电流为其中，为非线性负荷基波电流，为非线性负荷的第h次谐波电流，h为谐波次数。

8)根据DG注入主动配电网的基波电流计算出该DG注入主动配电网的第h次谐波电流为其中，为DG注入主动配电网的基波电流，为DG注入主动配电网的第h次谐波电流，h为谐波次数。

9)根据公式计算主动配电网中所有谐波源节点向主动配电网输入的各次谐波电流其中，为谐波源节点i向主动配电网输入的第h次谐波电流，为非线性负荷的第h次谐波电流，为DG注入主动配电网的第h次谐波电流，h为谐波次数。

10)根据公式计算主动配电网第h次谐波潮流，得出各独立节点i的第h次谐波电压其中，为独立节点i的第h次谐波电压，V^(h)为主动配电网中各独立节点i的第h次谐波电压组成的n×1阶相量矩阵，即i＝1,2,…,n，I(h)为主动配电网中各谐波源节点向主动配电网输入的第h次谐波电流组成的n×1阶相量矩阵，即i＝1,2,…,n，T为主动配电网的道路矩阵，T为n×b阶矩阵，T^T为矩阵T的转置矩阵，为支路l的第h次谐波阻抗，为基于对角元素分别为各支路l的第h次谐波阻抗组成的b×b阶对角矩阵，即l＝1,2,…,b，B为主动配电网的回-支关联矩阵，B为d×b阶矩阵，BT为矩阵B的转置矩阵，为矩阵的逆矩阵，即为主动配电网的第h次谐波回路阻抗矩阵，为d×d阶矩阵，且有n为主动配电网独立节点数，b为主动配电网总支路数，d为主动配电网连支数，h为谐波次数。

11)根据不同谐波次数h，重复步骤7)、步骤8)、步骤9)和步骤10)即可求出主动配电网所有节点的各次谐波电压分布情况，以及计算出各节点的电压总谐波畸变率和电压单个谐波畸变率。

步骤7)中各节点的电压总谐波畸变率THD的计算公式为：

$\mathrm{TH}{D}_i(\%)=\frac{\sqrt{\underset{h=2}{\overset{H}{\mathrm{\Σ}}}{\left|{\stackrel{\·}{V}}_i^{\left(h\right)}\right|}^2}}{\left|{\stackrel{\·}{V}}_i^{\left(1\right)}\right|}\×100\%---\left(11\right)$

电压单个谐波畸变率IHD的计算公式为：

${\mathrm{IhD}}_i^{\left(h\right)}(\%)=\frac{\left|{\stackrel{\·}{V}}_i^{\left(h\right)}\right|}{\left|{\stackrel{\·}{V}}_i^{\left(1\right)}\right|}\×100\%---\left(12\right)$

其中H为谐波的最高次数，h为谐波次数，i为节点序号，为节点i的第h次谐波电压，为节点i的基波电压，i＝1,2,…,n，n为主动配电网独立节点数。

算例仿真

如图2为33母线主动配电网测试系统，基于本发明算法的主动配电网谐波潮流计算结果如表1所示(算法的收敛精度为10^-6)。

表133母线主动配电网谐波潮流结果

从表1中可看出，仿真测试验证了本发明一种主动配电网谐波潮流分析方法的有效性和正确性。

Claims (2)

1.一种主动配电网谐波潮流分析方法，其特征是该方法包括以下步骤：

步骤A、获取主动配电网的网络参数，如针对具有N个节点、K个环网、W个非线性负荷和M个DG的主动配电网，假设首节点作为参考节点和平衡节点，则主动配电网独立节点数为n＝N-1，独立支路数即树支数为c＝n，关联支路数即连支数为d＝K，总支路数为b＝c+d，其中DG为分布式电源的英文缩写且代表分布式电源；

步骤B、分析非线性负荷运行特性，通过傅氏分析公式计算出其各次谐波电流分量与其基波电流的比值C(h)，其中，C(h)为非线性负荷各谐波电流分量与其基波电流的比值，h为谐波次数；

步骤C、分析DG并网接口电力电子装置的工作特性，通过傅氏分析公式计算出其注入主动配电网各次谐波电流分量与其基波电流的比值ρ(h)，其中，ρ(h)为DG注入主动配电网第h次谐波电流分量与其基波电流的比值，h为谐波次数；

步骤D、通过常规潮流算法计算主动配电网基波潮流，得出各非线性负荷节点的基波电压和各DG并网节点的基波电压其中，为非线性负荷节点的基波电压，为DG并网节点的基波电压，下标“NL”表示非线性负荷，下标“DG”表示分布式电源，下标“i”表示非线性负荷节点和DG并网节点各自对应的节点序号，上标“(1)”表示基波分量；

步骤E、根据主动配电网中非线性负荷的有功功率P_NL,i和无功功率Q_NL,i，计算出非线性负荷基波电流为其中，为非线性负荷基波电流，为非线性负荷节点的基波电压，P_NL,i和Q_NL,i分别为非线性负荷的有功功率和无功功率，上标“*”表示取共轭复数，j为复数的虚部单位；

步骤F、根据DG输出的有功功率P_DG,i和无功功率Q_DG,i，计算出DG注入主动配电网的基波电流为其中，为DG注入主动配电网的基波电流，为DG并网节点的基波电压，P_DG,i和Q_DG,i分别为DG输出的有功功率和无功功率，上标“*”表示取共轭复数，j为复数的虚部单位；

步骤G、根据非线性负荷基波电流计算出该非线性负荷的第h次谐波电流为其中，为非线性负荷基波电流，为非线性负荷的第h次谐波电流，h为谐波次数；

步骤H、根据DG注入主动配电网的基波电流计算出该DG注入主动配电网的第h次谐波电流为其中，为DG注入主动配电网的基波电流，为DG注入主动配电网的第h次谐波电流，h为谐波次数；

步骤I、根据公式计算主动配电网中所有谐波源节点向主动配电网输入的各次谐波电流其中，为谐波源节点i向主动配电网输入的第h次谐波电流，为非线性负荷的第h次谐波电流，为DG注入主动配电网的第h次谐波电流，h为谐波次数；

步骤J、根据公式计算主动配电网第h次谐波潮流，得出各独立节点i的第h次谐波电压其中，为独立节点i的第h次谐波电压，V^(h)为主动配电网中各独立节点i的第h次谐波电压组成的n×1阶相量矩阵，即i＝1,2,…,n，I^(h)为主动配电网中各谐波源节点向主动配电网输入的第h次谐波电流组成的n×1阶相量矩阵，即i＝1,2,,,n，T为主动配电网的道路矩阵，T为n×b阶矩阵，T^T为矩阵T的转置矩阵，为支路l的第h次谐波阻抗，为基于对角元素分别为各支路l的第h次谐波阻抗组成的b×b阶对角矩阵，即l＝1,2,…,b，B为主动配电网的回-支关联矩阵，B为d×b阶矩阵，B^T为矩阵B的转置矩阵，为矩阵的逆矩阵，即为主动配电网的第h次谐波回路阻抗矩阵，为d×d阶矩阵，且有n为主动配电网独立节点数，b为主动配电网总支路数，d为主动配电网连支数，h为谐波次数；

步骤K、根据不同谐波次数h，重复步骤G、步骤H、步骤I和步骤J即可求出主动配电网所有节点的各次谐波电压分布情况，以及计算出各节点的电压总谐波畸变率和电压单个谐波畸变率。

2.如权利要求1所述的主动配电网谐波潮流分析方法，其特征是，主动配电网的网络参数包括系统总节点数，总支路数，节点序号，支路序号，独立节点数，独立支路数，关联支路数，环网个数，参考节点，平衡节点，支路阻抗，线性负荷功率，非线性负荷功率，线性负荷节点序号，非线性负荷节点序号，网络拓扑结构，DG并网个数，DG输出功率，以及DG并网的节点序号。