CN105548736A - 一种谐波责任分摊量化评估方法 - Google Patents

Abstract

本发明提供了一种谐波责任分摊量化评估方法，包括：S1)对需要量化分析谐波责任的公共耦合节点电气数据进行测量；S2)根据测得的电气数据，分析和计算等效电路系统侧与用户侧的谐波责任分摊；S3)分析和计算等效电路用户侧各用户之间的谐波责任分摊；S4)重复步骤S1-S3复数次，将步骤S2和S3得到的结果取平均值作为最终的谐波责任分摊量化评估结果。该方法可较为全面的量化计算节点处各用户和系统侧的谐波责任，为电力部门和用户进行科学、公正的谐波责任分摊提供技术支撑。

Description

一种谐波责任分摊量化评估方法

技术领域

本发明属于电能质量领域，具体涉及一种谐波责任分摊量化评估方法。

背景技术

随着工业的快速发展，大量非线性负荷接入电网，并向电网注入谐波电流，且在各母线上形成谐波电压，从而影响电网中其它电气设备的正常工作。如何准确识别谐波源并将各个谐波源所产生的谐波电流分离，确定其各自责任，仍然是一个值得探讨的问题。工程上一般根据单个污染用户的监测数据进行谐波发射水平评估，用电企业依此确定谐波治理方案，电力部门依此进行经济惩奖；但公共连接点一般可能存在多个谐波源用户，用户之间如何精确地进行谐波责任分摊是进行谐波治理和经济惩奖的科学根据；而对于采用谐波治理措施的用户，当其投入滤波器组后，外界谐波电流除注入系统外，还有部分电流会流入滤波支路。外界谐波电流的注入一方面会造成滤波器过载运行，甚至会导致电容器爆裂事故发生；另一方面也会造成对各用户谐波源谐波责任分摊的判别不准确。

发明内容

本发明所要解决的技术问题是：针对现有多用户谐波源系统中存在的谐波责任分摊问题，提出了一种谐波责任分摊量化评估方法。该方法可较为全面的量化计算节点处各用户和系统侧的谐波责任，为电力部门和用户进行科学、公正的谐波责任分摊提供技术支撑。

本发明所采用的技术方案是：一种谐波责任分摊量化评估方法，包括：

S1)对需要量化分析谐波责任的公共耦合节点电气数据进行测量；

S2)根据测得的电气数据，分析和计算等效电路系统侧与用户侧的谐波责任分摊；

S3)分析和计算等效电路用户侧各用户之间的谐波责任分摊；

S4)重复步骤S1-S3复数次，将步骤S2和S3得到的结果取平均值作为最终的谐波责任分摊量化评估结果。

所述的方法，步骤S1的测量方法包括：通过传感器测量需评估谐波责任的公共耦合节点电气数据，电气数据包括基波和谐波电压、电流、相角数据。

所述的方法，步骤S2的计算方法包括：用诺顿等效电路，将系统侧和用户侧均等效为谐波阻抗和谐波电流源的并联，然后分别计算系统侧和用户侧的各次谐波电流责任分摊和系统侧和用户侧的各次谐波电压责任分摊。

所述的方法，系统侧和用户侧的h次谐波电流责任分摊的计算方法包括：

计算I_sh-pcc-pf和I_ech-pcc-pf的比值即得到所述的电流责任分摊，计算公式为

其中，为系统侧h次等效谐波电流源，为用户侧h次等效谐波电流源，为和共同作用时系统侧流向用户侧的h次谐波电流，I_sh-pcc-pf为在上的投影值，为单独作用时用户侧流向系统侧的h次谐波电流，I_ech-pcc-pf为在上的投影值，为单独作用时系统侧流向用户侧的h次谐波电流；α_ih、α_ish-pcc、α_iech-pcc分别为电流的相角；

并且，和的计算方法为

${\stackrel{\·}{I}}_{sh-pcc}=({\stackrel{\·}{I}}_h+{\stackrel{\·}{U}}_h/Z_{sh})*Z_{sh}/(Z_{sh}+Z_{eq})$

${\stackrel{\·}{I}}_{ech-pcc}=({\stackrel{\·}{U}}_h/Z_{eq}-{\stackrel{\·}{I}}_h)*Z_{eq}/(Z_{sh}+Z_{eq})$

其中，为和共同作用时公共耦合节点处的h次谐波电压，Z_eq为的谐波阻抗，Z_sh为的谐波阻抗。

所述的方法，系统侧和用户侧的h次谐波电压责任分摊的计算方法包括：

计算U_sh-pcc-pf和U_ech-pcc-pf的比值即得到所述的电压责任分摊，计算公式为

$U_{sh-pcc-pf}:U_{ech-pcc-pf}=\left|{\stackrel{\·}{U}}_{sh-pcc}\right|cos({\mathrm{\α}}_{uh}-{\mathrm{\α}}_{ush-pcc}):\left|{\stackrel{\·}{U}}_{ech-pcc}\right|cos({\mathrm{\α}}_{uh}-{\mathrm{\α}}_{uech-pcc})$

其中，U_sh-pcc-pf为在上的投影值，为单独作用时公共耦合节点处的h次谐波电压，U_ech-pcc-pf为在上的投影，为单独作用时公共耦合节点处的h次谐波电压，α_uh、α_ush-pcc、α_uech-pcc分别为电压的相角；

并且，和的计算方法为

${\stackrel{\·}{U}}_{sh-pcc}={\stackrel{\·}{I}}_{sh-pcc}{Z}_{eq}$

${\stackrel{\·}{U}}_{ech-pcc}={\stackrel{\·}{I}}_{ech-pcc}{Z}_{sh}.$

所述的方法，步骤S3的计算方法包括：根据诺顿定理，将用户侧每个用户的负荷支路等效成用户谐波电流源、谐波阻抗和滤波阻抗的并联，然后分别计算用户侧的各用户、各次谐波电流责任分摊和各用户、各次谐波电压责任分摊。

所述的方法，设有n个用户负荷支路，各用户谐波电流责任分摊的计算方法是计算两个不同支路i、j的I_hci-pcc-pf和I_hcj-pcc-pf的比值，其计算公式为

$I_{hci-pcc-pf}:I_{hcf-pcc-pf}=|{\stackrel{\·}{I}}_{hi}+\frac{Z_{chi}+Z_{fhi}}{Z_{chi}{Z}_{fhi}}{\stackrel{\·}{U}}_h|\cos ({\mathrm{\α}}_{iech}-{\mathrm{\α}}_{ihci-pcc}):|{\stackrel{\·}{I}}_{hj}+\frac{Z_{chj}+Z_{fhj}}{Z_{chj}{Z}_{fhj}}{\stackrel{\·}{U}}_h|\cos ({\mathrm{\α}}_{iech}-{\mathrm{\α}}_{jhci-pcc})$

其中，I_hci-pcc-pf为在方向上的投影值，为第i个支路h次用户谐波电流源单独作用时流入系统的谐波电流，为h次用户总体等效谐波电流源，为所有h次用户谐波电流源共同作用时在第i个支路上的谐波电流，Z_chi为第i个支路h次谐波阻抗，Z_fhi为第i个支路h次滤波阻抗，为所有h次用户谐波电流源共同作用时公共耦合节点处产生的谐波电压，α_iech、α_ihci-pcc分别为电流的相角；所有带下标i的参数符号换成j即为第j个支路相对应的参数；

然后重复上述计算步骤，直至计算出所有不同支路之间的电流责任分摊比例，即得到各用户谐波电流责任分摊。

所述的方法，各用户谐波电压责任分摊的计算方法是计算两个不同支路i、j的U_hci-pcc-pf和U_hcj-pcc-pf的比值，其计算公式为

$U_{hci-pcc-pf}:U_{hcj-pcc-pj}=|{\stackrel{\·}{I}}_{hi}+\frac{Z_{chi}+Z_{fhi}}{Z_{chi}{Z}_{fhi}}{\stackrel{\·}{U}}_h|\cos ({\mathrm{\α}}_{uech}-{\mathrm{\α}}_{uhci-pcc}):|{\stackrel{\·}{I}}_{hj}+\frac{Z_{chj}+Z_{fhj}}{Z_{chj}{Z}_{fhj}}{\stackrel{\·}{U}}_h|\cos ({\mathrm{\α}}_{uech}-{\mathrm{\α}}_{uhcj-pcc})$

其中，U_hci-pcc-pf为在方向上的投影值，为第i个支路h次用户谐波电流源单独作用时公共耦合节点处产生的谐波电压，为公共耦合节点处产生的h次等效谐波电压，其值与相等；α_uech、α_uhci-pcc分别为电压的相角；所有带下标i的参数符号换成j即为第j个支路相对应的参数；

然后重复上述计算步骤，直至计算出所有不同支路之间的电压责任分摊比例，即得到各用户谐波电压责任分摊。

本发明的优点：该方法针对现有多用户谐波源系统中存在的谐波责任分摊问题，提出了一种逐步分摊计算电压、电流的谐波责任分摊方法，该方法综合考虑了系统中的谐波治理装置、谐波电压和电流的初始相角和谐波源用户数等影响因素，可较为全面的量化计算节点处各用户和系统侧的谐波责任，为电力部门和用户进行科学、公正的谐波责任分摊提供技术支撑。

附图说明

附图1为谐波责任分摊量化评估计算流程图。

附图2为系统和用户等效电路。

附图3为单谐波源的等效电路。

附图4为多谐波源等效电路。

具体实施方式

本发明提供了一种谐波责任分摊量化评估方法，包括：S1)对需要量化分析谐波责任的pcc(公共联接点)节点电压、电流以及谐波数据进行测量；S2)分析计算用户侧与系统侧的谐波责任，应用诺顿原理将用户等效为一个谐波源，量化评估pcc节点两侧的谐波责任；S3)分析等效用户侧各用户之间的谐波责任，分别进行谐波电压和谐波电流责任分摊；S4)按上述步骤多次测量计算，取算术平均值，得到各用户和系统的量化谐波责任值。该方法针对现有多用户谐波源系统中存在的谐波责任分摊问题，提出了一种逐步分摊计算电压、电流的谐波责任分摊方法，该方法综合考虑了系统中的谐波治理装置、谐波电压和电流的初始相角和谐波源用户数等影响因素，可较为全面的量化计算节点处各用户和系统侧的谐波责任，为电力部门和用户进行科学、公正的谐波责任分摊提供理论依据和技术支撑。

下面结合附图进一步详述本发明。

一种谐波责任分摊量化评估方法，包括：

S1)对需要量化分析谐波责任的pcc(公共联接点)节点电压、电流以及谐波数据进行测量；

S2)分析计算用户侧与系统侧的谐波责任，应用诺顿原理将用户等效为一个谐波源，量化评估pcc节点两侧的谐波责任；

S3)分析等效用户侧各用户之间的谐波责任，分别进行谐波电压和谐波电流责任分摊；

S4)按上述步骤多次测量计算，取算术平均值，得到各用户和系统的量化谐波责任值。

上述步骤S1中，测量需评估谐波责任的pcc节点数据，包括基波和谐波电压、电流数据等。

上述步骤S2中，采用诺顿等效电路，将系统和用户均等效为谐波阻抗和恒定电流源的并联，并假设PCC处电压、电流向量及两侧谐波阻抗已知，图2所示的等效电路对PCC处的某次谐波电压和电流进行分析，图3为分解后所示的h次谐波对应的等效电路。

图2中，分别为等效用户侧谐波电流源和系统侧谐波电流源；Z_eq、Z_sh分别为对应的谐波阻抗；分别为单独作用时系统侧流向用户侧的h次谐波电流和PCC处的h次谐波电压；分别为单独作用时用户侧流向系统侧的h次谐波电流和PCC处的h次谐波电压；分别为和共同作用时系统侧流向用户侧的h次谐波电流和PCC处的h次谐波电压。

根据图2，可得到系统侧和等效用户侧谐波电流源的电流模型：

$\left\{\begin{array}{c}\stackrel{\·}{I_{sh}}=\stackrel{\·}{I_h}+\stackrel{\·}{U_h}/Z_{sh}\\ \stackrel{\·}{I_{ech}}=\stackrel{\·}{U_h}/Z_{eq}-\stackrel{\·}{I_h}\end{array}\right.$

对各用户单谐波源电路进行分析，等效电路如图3所示。

分别单独作用时流向系统侧的h次谐波电流和在PCC处产生的h次谐波电压可表示为：

$\left\{\begin{array}{c}{\stackrel{\·}{I}}_{sh-pcc}={\stackrel{\·}{I}}_{sh}*z_{sh}/(Z_{sh}+Z_{eq})\\ {\stackrel{\·}{I}}_{ech-pcc}={\stackrel{\·}{I}}_{ech}*Z_{eq}/(Z_{sh}+Z_{eq})\\ {\stackrel{\·}{U}}_{sh-pcc}={\stackrel{\·}{I}}_{sh-pcc}{Z}_{eq}\\ {\stackrel{\·}{U}}_{ech-pcc}={\stackrel{\·}{I}}_{ech-pcc}{Z}_{sh}\end{array}\right.$

进一步可得：

${\stackrel{\·}{I}}_{sh-pcc}=({\stackrel{\·}{I}}_h+{\stackrel{\·}{U}}_h/Z_{sh})*Z_{sh}/(Z_{sh}+Z_{eq})$

${\stackrel{\·}{I}}_{ech-pcc}=({\stackrel{\·}{U}}_h/Z_{eq}-{\stackrel{\·}{I}}_h)*Z_{eq}/(Z_{sh}+Z_{eq})$

根据上式，可求出在上的投影值I_sh-pcc-pf，表示如下：

$I_{sh-pcc-pf}=\left|{\stackrel{\·}{I}}_{sh-pcc}\right|cos({\mathrm{\α}}_{ih}-{\mathrm{\α}}_{ish-pcc})$

式中，α_ih、α_ish-pcc分别为电流的相角。

同理，可求出在上的投影值I_ech-pcc-pf，表示如下：

$I_{ech-pcc-pf}=\left|{\stackrel{\·}{I}}_{ech-pcc}\right|cos({\mathrm{\α}}_{ih}-{\mathrm{\α}}_{iech-pcc})$

式中，α_iech-pcc为电流的相角。

进而可得到等效系统侧和用户侧的谐波电流责任分摊：

$I_{sh-pcc-pf}:I_{ech-pcc-pf}=\left|{\stackrel{\·}{I}}_{sh-pcc}\right|cos({\mathrm{\α}}_{ih}-{\mathrm{\α}}_{ish-pcc}):\left|{\stackrel{\·}{I}}_{ech-pcc}\right|cos({\mathrm{\α}}_{ih}-{\mathrm{\α}}_{iech-pcc})$

同理，可求出等效系统侧和用户侧的谐波电压责任分摊：

$U_{sh-pcc-pf}:U_{ech-pcc-pf}=\left|{\stackrel{\·}{U}}_{sh-pcc}\right|cos({\mathrm{\α}}_{uh}-{\mathrm{\α}}_{ush-pcc}):\left|{\stackrel{\·}{U}}_{ech-pcc}\right|cos({\mathrm{\α}}_{uh}-{\mathrm{\α}}_{uech-pcc})$

上述步骤S3中，对用户侧多谐波源谐波责任分别进行计算，需要分析多个谐波源接入时的谐波电压、谐波电流，即分析的多个负荷工况数据必须是同时测量的。图4为多用户谐波源作用时谐波等效电路，Z_ch1、Z_ch2…Z_chn分别是各支路等效谐波阻抗，Z_fh1、Z_fh2…Z_fhn分别是各用户滤波支路等效阻抗；各用户谐波源单独作用时流入系统的谐波电流分别是各用户谐波源单独作用时在PCC处产生的谐波电压分别是分别是所有用户谐波源共同作用时流入系统的谐波电流和在PCC处产生的谐波电压，这两个参数值分别等于等效用户谐波电流和等效用户谐波电压分别是所有用户谐波源共同作用时在各谐波支路上的谐波电流。

根据图4可知，

$\left\{\begin{array}{c}{\stackrel{\·}{I}}_h={\stackrel{\·}{I}}_{h1}+{\stackrel{\·}{I}}_{h2}+\mathrm{...}+{\stackrel{\·}{I}}_{hn}\\ {\stackrel{\·}{U}}_h={\stackrel{\·}{U}}_{h1}={\stackrel{\·}{U}}_{h2}=\mathrm{...}={\stackrel{\·}{U}}_{hn}\end{array}\right.$

图4所有开关S1～Sn均闭合，电流源单独作用时，可得：

$Z_{eq}=\frac{Z_{ch1}{Z}_{fh1}\mathrm{...}{Z}_{chn}{Z}_{fhn}}{\underset{i\≠fhn}{\underset{i=ch1-fhn}{\Π}}{Z}_i+\underset{i\≠chn}{\underset{i=ch1-fhn}{\Π}}{Z}_i+\underset{i\≠fh1}{\underset{i=ch1-fhn}{\Π}}{Z}_i+\underset{i\≠ch1}{\underset{i=ch1-fhn}{\Π}}{Z}_i}$

其中，(连乘的项目包括下标ch1…chn-1、fh1…fhn，只是不包括chn)表示从下标ch1的阻抗连乘到fhn的阻抗，但不连乘chn。

$\left\{\begin{array}{c}{\stackrel{\·}{I}}_{hci-pcc}=\frac{Z_{eq}}{Z_{eq}+Z_{sh}}{\stackrel{\·}{I}}_{chi}\\ {\stackrel{\·}{I}}_{chi}={\stackrel{\·}{I}}_{hi}+\frac{Z_{chi}+Z_{fhi}}{Z_{chi}{Z}_{fhi}}{\stackrel{\·}{U}}_h\end{array}\right.$

由上式可知，在等效用户谐波电流(数值与相等)方向上的谐波投影值I_hci-pcc-pf，可表示为：

$I_{hci-pcc-pf}=\left|{\stackrel{\·}{I}}_{hci-pcc}\right|cos({\mathrm{\α}}_{iech}-{\mathrm{\α}}_{ihci-pcc})$

其中，α_iech、α_ihci-pcc分别为电流的相角。

由上式，即可得到所有用户谐波源之间谐波电流的责任分摊比例：

$I_{hci-pcc-pf}:I_{hcj-pcc-pf}=|{\stackrel{\·}{I}}_{hi}+\frac{Z_{chi}+Z_{fhi}}{Z_{chi}{Z}_{fhi}}{\stackrel{\·}{U}}_h|\cos ({\mathrm{\α}}_{iech}-{\mathrm{\α}}_{ihci-pcc}):|{\stackrel{\·}{I}}_{hj}+\frac{Z_{chj}+Z_{fhj}}{Z_{chj}{Z}_{fhj}}{\stackrel{\·}{U}}_h|\cos ({\mathrm{\α}}_{iech}-{\mathrm{\α}}_{jhci-pcc})$

当开关S1、S2、S3…Sn同时闭合，电流源单独作用时，可得：

${\stackrel{\·}{U}}_{hci-pcc}={\stackrel{\·}{I}}_{hci}{Z}_{sh}$

可得到在等效用户谐波电压(数值与相等)方向上的投影为U_ech-pcc-pf，可表示为：

$U_{ech-pcc-pf}=\left|{\stackrel{\·}{U}}_{hci-pcc}\right|cos({\mathrm{\α}}_{uech}-{\mathrm{\α}}_{uhci-pcc})$

式中，α_uech、α_uhci-pcc分别为U_ech、U_hci-pcc的相角。

由上式，即可得到所用用户谐波源之间谐波责任分摊：

$U_{hci-pcc-pf}:U_{hcj-pcc-pj}=|{\stackrel{\·}{I}}_{hi}+\frac{Z_{chi}+Z_{fhi}}{Z_{chi}{Z}_{fhi}}{\stackrel{\·}{U}}_h|\cos ({\mathrm{\α}}_{uech}-{\mathrm{\α}}_{uhci-pcc}):|{\stackrel{\·}{I}}_{hj}+\frac{Z_{chj}+Z_{fhj}}{Z_{chj}{Z}_{fhj}}{\stackrel{\·}{U}}_h|\cos ({\mathrm{\α}}_{uech}-{\mathrm{\α}}_{uhcj-pcc})$

其中，i,j＝1,2,...n(i≠j)。

上述步骤S4中，考虑系统部分阻抗具有时变性，需多次测量计算责任分摊比例，然后取平均值，得到最终评估结果。

Claims (8)

1.一种谐波责任分摊量化评估方法，其特征在于包括：

S1)对需要量化分析谐波责任的公共耦合节点电气数据进行测量；

S2)根据测得的电气数据，分析和计算等效电路系统侧与用户侧的谐波责任分摊；

S3)分析和计算等效电路用户侧各用户之间的谐波责任分摊；

S4)重复步骤S1-S3复数次，将步骤S2和S3得到的结果取平均值作为最终的谐波责任分摊量化评估结果。

2.根据权利要求1所述的方法，其特征在于，步骤S1的测量方法包括：通过传感器测量需评估谐波责任的公共耦合节点电气数据，电气数据包括基波和谐波电压、电流、相角数据。

3.根据权利要求2所述的方法，其特征在于，步骤S2的计算方法包括：用诺顿等效电路，将系统侧和用户侧均等效为谐波阻抗和谐波电流源的并联，然后分别计算系统侧和用户侧的各次谐波电流责任分摊和系统侧和用户侧的各次谐波电压责任分摊。

4.根据权利要求3所述的方法，其特征在于，系统侧和用户侧的h次谐波电流责任分摊的计算方法包括：

计算I_sh-pcc-pf和I_ech-pcc-pf的比值即得到所述的电流责任分摊，计算公式为

$I_{sh-pcc-pf}:I_{ech-pcc-pf}=\left|{\stackrel{\·}{I}}_{sh-pcc}\right|cos({\mathrm{\α}}_{ih}-{\mathrm{\α}}_{ish-pcc}):\left|{\stackrel{\·}{I}}_{ech-pcc}\right|cos({\mathrm{\α}}_{ih}-{\mathrm{\α}}_{iech-pcc})$

其中，为系统侧h次等效谐波电流源，为用户侧h次等效谐波电流源，为和共同作用时系统侧流向用户侧的h次谐波电流，I_sh-pcc-pf为在上的投影值，为单独作用时用户侧流向系统侧的h次谐波电流，I_ech-pcc-pf为在上的投影值，为单独作用时系统侧流向用户侧的h次谐波电流；α_ih、α_ish-pcc、α_iech-pcc分别为电流的相角；

并且，和的计算方法为

${\stackrel{\·}{I}}_{sh-pcc}=({\stackrel{\·}{I}}_h+{\stackrel{\·}{U}}_h/Z_{sh})*Z_{sh}/(Z_{sh}+Z_{eq})$

${\stackrel{\·}{I}}_{ech-pcc}=({\stackrel{\·}{U}}_h/Z_{eq}-{\stackrel{\·}{I}}_h)*Z_{eq}/(Z_{sh}+Z_{eq})$

其中，为和共同作用时公共耦合节点处的h次谐波电压，Z_eq为的谐波阻抗，Z_sh为的谐波阻抗。

5.根据权利要求4所述的方法，其特征在于，系统侧和用户侧的h次谐波电压责任分摊的计算方法包括：

计算U_sh-pcc-pf和U_ech-pcc-pf的比值即得到所述的电压责任分摊，计算公式为

$U_{sh-pcc-pf}:U_{ech-pcc-pf}=\left|{\stackrel{\·}{U}}_{sh-pcc}\right|cos({\mathrm{\α}}_{uh}-{\mathrm{\α}}_{ush-pcc}):\left|{\stackrel{\·}{U}}_{ech-pcc}\right|cos({\mathrm{\α}}_{uh}-{\mathrm{\α}}_{uech-pcc})$

其中，U_sh-pcc-pf为在上的投影值，为单独作用时公共耦合节点处的h次谐波电压，U_ech-pcc-pf为在上的投影，为单独作用时公共耦合节点处的h次谐波电压，α_uh、α_ush-pcc、α_uech-pcc分别为电压的相角；

并且，和的计算方法为

${\stackrel{\·}{U}}_{sh-pcc}={\stackrel{\·}{I}}_{sh-pcc}{Z}_{eq}$

${\stackrel{\·}{U}}_{ech-pcc}={\stackrel{\·}{I}}_{ech-pcc}{Z}_{sh}.$

6.根据权利要求2所述的方法，其特征在于，步骤S3的计算方法包括：根据诺顿定理，将用户侧每个用户的负荷支路等效成用户谐波电流源、谐波阻抗和滤波阻抗的并联，然后分别计算用户侧的各用户、各次谐波电流责任分摊和各用户、各次谐波电压责任分摊。

7.根据权利要求6所述的方法，其特征在于：

设有n个用户负荷支路，各用户谐波电流责任分摊的计算方法是计算两个不同支路i、j的I_hci-pcc-pf和I_hcj-pcc-pf的比值，其计算公式为

$I_{hci-pcc-pf}:I_{hcj-pcc-pj}=|{\stackrel{\·}{I}}_{hi}+\frac{Z_{chi}+Z_{fhi}}{Z_{chi}{Z}_{fhi}}{\stackrel{\·}{U}}_h|\cos ({\mathrm{\α}}_{iech}-{\mathrm{\α}}_{ihci-pcc}):|{\stackrel{\·}{I}}_{hj}+\frac{Z_{chj}+Z_{fhj}}{Z_{chj}{Z}_{fhj}}{\stackrel{\·}{U}}_h|\cos ({\mathrm{\α}}_{iech}-{\mathrm{\α}}_{jhch-pcc})$

其中，I_hci-pcc-pf为在方向上的投影值，为第i个支路h次用户谐波电流源单独作用时流入系统的谐波电流，为h次用户总体等效谐波电流源，为所有h次用户谐波电流源共同作用时在第i个支路上的谐波电流，Z_chi为第i个支路h次谐波阻抗，Z_fhi为第i个支路h次滤波阻抗，为所有h次用户谐波电流源共同作用时公共耦合节点处产生的谐波电压，α_iech、α_ihci-pcc分别为电流的相角；所有带下标i的参数符号换成j即为第j个支路相对应的参数；

然后重复上述计算步骤，直至计算出所有不同支路之间的电流责任分摊比例，即得到各用户谐波电流责任分摊。

8.根据权利要求7所述的方法，其特征在于：

各用户谐波电压责任分摊的计算方法是计算两个不同支路i、j的U_hci-pcc-pf和U_hcj-pcc-pf的比值，其计算公式为

$U_{hci-pcc-pf}:U_{hcj-pcc-pj}=|{\stackrel{\·}{I}}_{hi}+\frac{Z_{chi}+Z_{fhi}}{Z_{chi}{Z}_{fhi}}{\stackrel{\·}{U}}_h|\cos ({\mathrm{\α}}_{uech}-{\mathrm{\α}}_{uhci-pcc}):|{\stackrel{\·}{I}}_{hj}+\frac{Z_{chj}+Z_{fhj}}{Z_{chj}{Z}_{fhj}}{\stackrel{\·}{U}}_h|\cos ({\mathrm{\α}}_{uech}-{\mathrm{\α}}_{uhcj-pcc})$

其中，U_hci-pcc-pf为在方向上的投影值，为第i个支路h次用户谐波电流源单独作用时公共耦合节点处产生的谐波电压；为公共耦合节点处产生的h次等效谐波电压，其值与相等；α_uech、α_uhci-pcc分别为电压的相角；所有带下标i的参数符号换成j即为第j个支路相对应的参数；

然后重复上述计算步骤，直至计算出所有不同支路之间的电压责任分摊比例，即得到各用户谐波电压责任分摊。