CN106019026A - 一种基于波形匹配的谐波电压责任的划分方法 - Google Patents

Abstract

本发明公开了一种基于波形匹配的谐波电压责任的划分方法，包括：步骤10）采集谐波测试数据，形成谐波数据序列；步骤20）设置基本匹配参数；步骤30）计算欧氏距离：根据步骤20）设置的基本匹配参数，对步骤10）形成的谐波数据序列按子序列顺序进行归一化处理，得到谐波样本数据，计算谐波样本数据中谐波电压与谐波电流的欧氏距离；步骤40）筛选子序列：计算每个子序列的相似度，筛选出相似度不小于设置值的子序列；步骤50）划分谐波责任：通过最小二乘法估算出筛选出的子序列的系统侧等效谐波阻抗，利用等效谐波阻抗的均值，计算谐波电压责任。该方法能够利用谐波幅值数据对馈线的谐波污染责任进行量化，适用于工程实际。

Description

一种基于波形匹配的谐波电压责任的划分方法

技术领域

本发明属于电能质量监控与分析技术领域，具体来说，涉及一种基于波形匹配的谐波电压责任的划分方法。

背景技术

随着电力系统的发展，各种传统非线性负载及光伏、风电等新能源大量接入电网，使得电网中的谐波污染越来越严重。国家标准对公用电网谐波的允许阈值做出了具体规定，但并不能科学地解决供用电双方在制定经济性评估标准上的争论。为了对谐波污染进行有效合理的经济性评估，需要在公共连接点对系统及用户承担的谐波责任进行定量划分。

近年来，“非干预式”方法成为谐波责任量化研究的一个主要方向。该方法通过测量公共连接点的谐波数据来估算系统侧谐波阻抗，进而实现责任划分。该类方法对系统的运行不构成影响，具体包括了线性回归法和波动量法。上述研究方法针对的是谐波向量数据，而目前变电站的电能质量监测系统仅仅给出谐波电压或电流的幅值数据，若进行额外的谐波相角测试，将大大增加测试的费用及工作量。谐波问题对供电企业安全稳定运行及供用电双方的经济效益带来了巨大影响，为了有效利用目前变电站的谐波实测数据，减少附加相位测量的工作量，降低背景谐波波动带来的计算误差，需要提出更加实用的谐波责任量化方法。

发明内容

技术问题：本发明所要解决的技术问题是：提供一种基于波形匹配的谐波电压责任的划分方法，该方法能够利用谐波幅值数据对馈线的谐波污染责任进行量化，避免对谐波相角的直接量测，适用于工程实际。

技术方案：为解决上述技术问题，本发明实施例采用一种基于波形匹配的谐波电压责任的划分方法，包括以下步骤：

步骤10)采集谐波测试数据，形成谐波数据序列；

步骤20)设置基本匹配参数；

步骤30)计算欧氏距离：根据步骤20)设置的基本匹配参数，对步骤10)形成的谐波数据序列按子序列顺序进行归一化处理，得到谐波样本数据，计算谐波样本数据中谐波电压与谐波电流的欧氏距离；

步骤40)筛选子序列：计算每个子序列的相似度，筛选出相似度不小于设置值的子序列；

步骤50)划分谐波责任：通过最小二乘法估算出筛选出的子序列的系统侧等效谐波阻抗，利用等效谐波阻抗的均值，计算谐波电压责任。

作为优选例，所述的步骤10)具体包括：利用谐波测量装置，采集公共连接点在测量时间内的谐波电压幅值和谐波电流幅值，并将该两种谐波测试数据形成如下式所示的谐波数据序列：

$\begin{array}{c}{U}_{pcc-h}\left(1\right),I_{pcc-h}\left(1\right)\\ U_{pcc-h}\left(2\right),I_{pcc-h}\left(2\right)\\ .\\ .\\ .\\ U_{pcc-h}\left(m\right),I_{pcc-h}\left(m\right)\end{array}$

式中，U_pcc-h表示谐波电压幅值；I_pcc-h表示谐波电流幅值；m表示谐波数据序列长度；U_pcc-h(1)表示第一次测量采集的谐波电压幅值，I_pcc-h(1)表示第一次测量采集的谐波电流幅值，U_pcc-h(2)表示第二次测量采集的谐波电压幅值，I_pcc-h(2)表示第二次测量采集的谐波电流幅值，U_pcc-h(m)表示第m次测量采集的谐波电压幅值，I_pcc-h(m)表示第m次测量采集的谐波电流幅值；该谐波数据序列中的每一行为其一子序列。

作为优选例，所述的步骤20)中，设置基本匹配参数包括设置滑动窗口长度L，滑动距离T及相似度设定值S。

作为优选例，所述的步骤30)具体包括：

步骤301)滑动窗口从谐波数据序列的第1组数据开始，将滑动窗口内U_pcc-h与I_pcc-h子序列分别归算到区间[-1，1]内，归一化函数为：

${\stackrel{\‾}{U}}_{pcc-h}\left(p\right)=\frac{2U_{pcc-h}\left(p\right)-U_{pcc-h-\min }-U_{pcc-h-\max }}{U_{pcc-h-\max }-U_{pcc-h-\min }},(p=1,2,3...,M)$

${\stackrel{\‾}{I}}_{pcc-h}\left(p\right)=\frac{2I_{pcc-h}\left(p\right)-I_{pcc-h-min}-I_{pcc-h-max}}{I_{pcc-h-max}-I_{pcc-h-\min }},(p=1,2,3...,M)$

式中，表示归一化后的谐波电压幅值，U_pcc-h表示归一化之前的谐波电压幅值，U_pcc-h-max表示U_pcc-h中的最大值，U_pcc-h-min表示U_pcc-h中的最小值，表示归一化后的谐波电流幅值，I_pcc-h表示归一化之前的谐波电流幅值，I_pcc-h-max表示I_pcc-h中的最大值，I_pcc-h-min表示I_pcc-h中的最小值，M为一个L长度滑动窗口内谐波电压幅值或谐波电流幅值的数据数量；p表示位于滑动窗口内的谐波电压幅值或谐波电流幅值的序号；

步骤302)按下式计算归一化后子序列中谐波电压幅值与谐波电流幅值的欧氏距离：

$DIS=\sqrt{\underset{p=1}{\overset{M}{\mathrm{\Σ}}}{\[{\stackrel{\‾}{U}}_{pcc-h}\left(p\right)-{\stackrel{\‾}{I}}_{pcc-h}\left(p\right)\]}^2}$

式中，DIS表示归一化后子序列的序列点与的欧氏距离；

步骤303)滑动窗口向后移动T，重复步骤301)，直至滑动窗口覆盖了整个谐波测试数据，得到K个子序列的欧氏距离。

作为优选例，所述的步骤40)具体包括：计算K个子序列的欧氏距离均值为每个子序列的谐波电压幅值与谐波电流幅值的相似度表示为：

$S\left(q\right)=\frac{\stackrel{\‾}{DIS}}{DIS\left(q\right)+\stackrel{\‾}{DIS}},(q=1,2,...K)$

式中，S(q)为第q个子序列的相似度，DIS(q)为第q个子序列的欧氏距离；

从K个子序列中筛选出相似度不小于相似度设定值S的子序列，设筛选出K'组子序列。

作为优选例，所述的步骤50)具体包括：

步骤501)对筛选出的K'组子序列利用最小二乘法估算出K'个系统侧等效谐波阻抗值，最小二乘法估算方法为：

将筛选出的K'组子序列中的每一组子序列代入下式：

${\stackrel{\‾}{U}}_{pcc-h}=Z_h\·{\stackrel{\‾}{I}}_{pcc-h}+U_S$

式中，Z_h表示该组子序列的系统侧等效谐波阻抗值，U_S表示该组子序列的系统侧背景谐波电压值；

通过下式求解出该组子序列的Z_h和U_S：

$f(Z_h,U_S)=\underset{p=1}{\overset{M}{\mathrm{\Σ}}}{\[{\stackrel{\‾}{U}}_{pcc-h}\left(p\right)-(Z_h\·{\stackrel{\‾}{I}}_{pcc-h}(p)+U_S)\]}^2$

$\frac{\∂f(Z_h,U_S)}{\∂Z_h}=0$

$\frac{\∂f(Z_h,U_S)}{\∂U_S}=0$

步骤502)对获得的K'个Z_h值，取均值为通过下式求解出测量时间内的谐波电压责任：

$T_C\left(l\right)\%=\frac{{\stackrel{\‾}{Z}}_h\·I_{pcc-h}\left(l\right)}{U_{pcc-h}\left(l\right)}\×100\%$

式中，T_C表示用户侧谐波电压责任百分数，l＝1,2,…m。

有益效果：与现有技术相比，本发明实施例具有以下有益效果：该方法能够利用实测数据准确计算出系统侧等效谐波阻抗值。传统的测算方法针对谐波向量数据进行分析，需要测量谐波电压和谐波电流的相角，而日常电能质量监测系统往往只测试谐波的幅值数据，传统方法不适用于工程实际。本实施例的方法，利用工程实测的谐波电压和谐波电流幅值数据准确计算出谐波阻抗值，方法简单有效，易于工程应用。

附图说明

图1为本发明实施例的流程图。

图2(a)为本发明实施例中谐波电压幅值采集数据波形图。

图2(b)为本发明实施例中谐波电流幅值采集数据波形图。

图3为本发明实施例中谐波电压责任计算结果图。

具体实施方式

下面结合实例和附图，对本发明实施例的技术方案做进一步的说明。

如图1所示，本发明实施例提供一种基于波形匹配的谐波电压责任的划分方法，包括以下步骤：

步骤10)采集谐波测试数据，形成谐波数据序列。

步骤20)设置基本匹配参数。设置基本匹配参数包括设置滑动窗口长度L，滑动距离T及相似度设定值S。

步骤30)计算欧氏距离：根据步骤20)设置的基本匹配参数，对步骤10)形成的谐波数据序列按子序列顺序进行归一化处理，得到谐波样本数据，计算谐波样本数据中谐波电压与谐波电流的欧氏距离。

步骤40)筛选子序列：计算每个子序列的相似度，筛选出相似度不小于设置值的子序列。

步骤50)划分谐波责任：通过最小二乘法估算出筛选出的子序列的系统侧等效谐波阻抗，利用等效谐波阻抗的均值，计算谐波电压责任。

在上述实施例中，所述的步骤10)具体包括：利用谐波测量装置，采集公共连接点(对应英文全称Point of Common Coupling；文中简称PCC点)在测量时间内的谐波电压幅值和谐波电流幅值，并将该两种谐波测试数据形成如下式所示的谐波数据序列：

$\begin{array}{c}{U}_{pcc-h}\left(1\right),I_{pcc-h}\left(1\right)\\ U_{pcc-h}\left(2\right),I_{pcc-h}\left(2\right)\\ .\\ .\\ .\\ U_{pcc-h}\left(m\right),I_{pcc-h}\left(m\right)\end{array}$

式中，U_pcc-h表示谐波电压幅值；I_pcc-h表示谐波电流幅值；m表示谐波数据序列长度；U_pcc-h(1)表示第一次测量采集的谐波电压幅值，I_pcc-h(1)表示第一次测量采集的谐波电流幅值，U_pcc-h(2)表示第二次测量采集的谐波电压幅值，I_pcc-h(2)表示第二次测量采集的谐波电流幅值，U_pcc-h(m)表示第m次测量采集的谐波电压幅值，I_pcc-h(m)表示第m次测量采集的谐波电流幅值；该谐波数据序列中的每一行为其一子序列。

在上述实施例中，所述的步骤30)具体包括：

步骤301)滑动窗口从谐波数据序列的第1组数据开始，将滑动窗口内U_pcc-h与I_pcc-h子序列分别归算到区间[-1，1]内，归一化函数为：

${\stackrel{\‾}{U}}_{pcc-h}\left(p\right)=\frac{2U_{pcc-h}\left(p\right)-U_{pcc-h-\min }-U_{pcc-h-\max }}{U_{pcc-h-\max }-U_{pcc-h-\min }},(p=1,2,3...,M)$

${\stackrel{\‾}{I}}_{pcc-h}\left(p\right)=\frac{2I_{pcc-h}\left(p\right)-I_{pcc-h-min}-I_{pcc-h-max}}{I_{pcc-h-max}-I_{pcc-h-\min }},(p=1,2,3...,M)$

式中，表示归一化后的谐波电压幅值，U_pcc-h表示归一化之前的谐波电压幅值，U_pcc-h-max表示U_pcc-h中的最大值，U_pcc-h-min表示U_pcc-h中的最小值，表示归一化后的谐波电流幅值，I_pcc-h表示归一化之前的谐波电流幅值，I_pcc-h-max表示I_pcc-h中的最大值，I_pcc-h-min表示I_pcc-h中的最小值，M为一个L长度滑动窗口内谐波电压幅值或谐波电流幅值的数据数量；p表示位于滑动窗口内的谐波电压幅值或谐波电流幅值的序号；

步骤302)按下式计算归一化后子序列中谐波电压幅值与谐波电流幅值的欧氏距离：

$DIS=\sqrt{\underset{p=1}{\overset{M}{\mathrm{\Σ}}}{\[{\stackrel{\‾}{U}}_{pcc-h}\left(p\right)-{\stackrel{\‾}{I}}_{pcc-h}\left(p\right)\]}^2}$

式中，DIS表示归一化后子序列的序列点与的欧氏距离；

步骤303)滑动窗口向后移动T，重复步骤301)，直至滑动窗口覆盖了整个谐波测试数据，得到K个子序列的欧氏距离。

在上述实施例中，所述的步骤40)具体包括：计算K个子序列的欧氏距离均值为每个子序列的谐波电压幅值与谐波电流幅值的相似度表示为：

$S\left(q\right)=\frac{\stackrel{\‾}{DIS}}{DIS\left(q\right)+\stackrel{\‾}{DIS}},(q=1,2,...K)$

式中，S(q)为第q个子序列的相似度，DIS(q)为第q个子序列的欧氏距离；

从K个子序列中筛选出相似度不小于相似度设定值S的子序列，设筛选出K'组子序列。

在上述实施例中，所述的步骤50)具体包括：

步骤501)对筛选出的K'组子序列利用最小二乘法估算出K'个系统侧等效谐波阻抗值，最小二乘法估算方法为：

将筛选出的K'组子序列中的每一组子序列代入下式：

${\stackrel{\‾}{U}}_{pcc-h}=Z_h\·{\stackrel{\‾}{I}}_{pcc-h}+U_S$

式中，Z_h表示该组子序列的系统侧等效谐波阻抗值，U_S表示该组子序列的系统侧背景谐波电压值；

通过下式求解出该组子序列的Z_h和U_S：

$f(Z_h,U_S)=\underset{p=1}{\overset{M}{\mathrm{\Σ}}}{\[{\stackrel{\‾}{U}}_{pcc-h}\left(p\right)-(Z_h\·{\stackrel{\‾}{I}}_{pcc-h}(p)+U_S)\]}^2$

$\frac{\∂f(Z_h,U_S)}{\∂Z_h}=0$

$\frac{\∂f(Z_h,U_S)}{\∂U_S}=0$

步骤502)对获得的K'个Z_h值，取均值为通过下式求解出测量时间内的谐波电压责任：

$T_C\left(l\right)\%=\frac{{\stackrel{\‾}{Z}}_h\·I_{pcc-h}\left(l\right)}{U_{pcc-h}\left(l\right)}\×100\%$

式中，T_C表示用户侧谐波电压责任百分数，l＝1,2,…m。

本发明实施例的方法，通过欧氏距离与滑动窗口寻找到测量时间内谐波电压幅值与谐波电流幅值相似度较高的波形，筛选出背景谐波稳定的子序列，再利用最小二乘法求解线性方程，估算出系统侧等效谐波阻抗，进而实现谐波责任划分。现有技术中，谐波阻抗的计算运用谐波向量数据，而目前电能质量日常监测系统仅仅测量谐波电压或谐波电流的幅值数据。本实施例方法的数据需求满足工程要求，可计算变压器低压侧或某条馈线的系统侧等效谐波阻抗值。这对于谐波分析，如谐波发射水平估算、谐波责任量化、谐波污染溯源等，具有重要指导意义，并且为谐波治理及制定奖惩机制提供帮助。

下面例举一具体实施例。

以某个110kV变电站的低压侧10kV母线为例，在母线下的某条馈线配置谐波测量装置进行谐波数据采集。

执行本发明实施例的测量方法，包括以下步骤：

步骤10)采集PCC点的5次谐波电压幅值和5次谐波电流幅值数据，谐波量测设备每3秒给出1组谐波电压幅值和谐波电流幅值数据，每小时可以采集1200组谐波数据。本实例选择了连续采集的20h的24000组谐波数据构成谐波数据序列，序列数据如图2所示，图中横坐标表示测试时间，纵坐标表示谐波电压幅值和谐波电流幅值。本实例中，谐波电压幅值和谐波电流幅值的单位分别为V和A。

步骤20)设置滑动窗口长度L为1h，即1个滑动窗口包括1200组谐波数据；设置滑动距离T为10min，即滑动一次向后移动240组数据；设置相似度S为0.9。

步骤30)滑动窗口第1组谐波数据序列开始，依次进行归一化并计算子序列的欧氏距离，直至滑动窗口覆盖了整个谐波测试数据的24000组数据，计算得到96个子序列的欧氏距离。欧氏距离按从上到下，从左到右的顺利排列如表1所示。

表1欧氏距离计算结果

70.57	84.59	282.66	34.30	272.88	304.16	83.24	65.52
								113.95	56.71	409.26	47.31	427.17	162.72	61.42	46.88
156.26	57.50	505.54	237.45	592.36	127.07	25.34	66.76
								84.74	65.01	498.10	58.01	490.89	132.95	23.10	114.24
226.55	46.68	358.75	69.33	99.46	201.01	13.34	142.34
								423.11	142.85	79.65	40.53	80.33	78.60	34.54	162.22
406.11	82.06	160.85	48.18	691.93	116.00	19.33	185.27
								213.13	326.69	209.56	63.21	676.70	110.91	91.65	239.54
136.38	282.14	80.97	91.74	670.09	44.33	83.22	315.95
								73.87	212.74	119.00	72.40	521.86	66.35	36.54	211.97
100.43	202.03	38.62	220.39	461.91	96.48	53.93	96.43
								84.83	340.73	15.12	212.46	275.96	79.64	56.22	51.83

步骤40)计算96个子序列的欧氏距离均值为177.14，筛选出3个相似度不小于0.9的子序列。每个子序列的谐波电压幅值与谐波电流幅值的相似度计算结果如表2所示。

表2相似度计算结果

0.72	0.68	0.39	0.84	0.39	0.37	0.68	0.73
								0.61	0.76	0.30	0.79	0.29	0.52	0.74	0.79
0.53	0.75	0.26	0.43	0.23	0.58	0.87	0.73
								0.68	0.73	0.26	0.75	0.27	0.57	0.88	0.61
0.44	0.79	0.33	0.72	0.64	0.47	0.93	0.55
								0.30	0.55	0.69	0.81	0.69	0.69	0.84	0.52
0.30	0.68	0.52	0.79	0.20	0.60	0.90	0.49
								0.45	0.35	0.46	0.74	0.21	0.61	0.66	0.43
0.57	0.39	0.69	0.66	0.21	0.80	0.68	0.36
								0.71	0.45	0.60	0.71	0.25	0.73	0.83	0.46
0.64	0.47	0.82	0.45	0.28	0.65	0.77	0.65
								0.68	0.34	0.92	0.45	0.39	0.69	0.76	0.77

步骤50)利用最小二乘法估算筛选出的3个子序列的系统侧等效谐波阻抗值，计算结果如表3所示。

表3系统侧等效谐波阻抗计算结果

相似度	0.92	0.93	0.90
				Zh(Ω)	1.39	1.37	1.56

3个Z_h取均值求解出的测量时间内的谐波电压责任如图3所示，图中横坐标表示时间，纵坐标表示谐波电压责任百分数。

通过变电站母线的额定电压等级与最小短路容量可以计算出系统侧的5次谐波阻抗为1.50Ω。实际由于系统侧与部分馈线并联，所以系统侧等效谐波阻抗应略小于1.50Ω。本实施中计算的结果为1.44Ω，计算结果准确，符合实际情况。

以上显示和描述了本发明的基本原理、主要特征和优点。本领域的技术人员应该了解，本发明不受上述具体实施例的限制，上述具体实施例和说明书中的描述只是为了进一步说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。本发明要求保护的范围由权利要求书及其等效物界定。

Claims (6)

1.一种基于波形匹配的谐波电压责任的划分方法，其特征在于，该方法包括以下步骤：

步骤10)采集谐波测试数据，形成谐波数据序列；

步骤20)设置基本匹配参数；

步骤30)计算欧氏距离：根据步骤20)设置的基本匹配参数，对步骤10)形成的谐波数据序列按子序列顺序进行归一化处理，得到谐波样本数据，计算谐波样本数据中谐波电压与谐波电流的欧氏距离；

步骤40)筛选子序列：计算每个子序列的相似度，筛选出相似度不小于设置值的子序列；

步骤50)划分谐波责任：通过最小二乘法估算出筛选出的子序列的系统侧等效谐波阻抗，利用等效谐波阻抗的均值，计算谐波电压责任。

2.根据权利要求1所述的基于波形匹配的谐波电压责任的划分方法，其特征在于，所述的步骤10)具体包括：利用谐波测量装置，采集公共连接点在测量时间内的谐波电压幅值和谐波电流幅值，并将该两种谐波测试数据形成如下式所示的谐波数据序列：

U_pcc-h(1),I_pcc-h(1)

U_pcc-h(2),I_pcc-h(2)

U_pcc-h(m),I_pcc-h(m)

式中，U_pcc-h表示谐波电压幅值；I_pcc-h表示谐波电流幅值；m表示谐波数据序列长度；U_pcc-h(1)表示第一次测量采集的谐波电压幅值，I_pcc-h(1)表示第一次测量采集的谐波电流幅值，U_pcc-h(2)表示第二次测量采集的谐波电压幅值，I_pcc-h(2)表示第二次测量采集的谐波电流幅值，U_pcc-h(m)表示第m次测量采集的谐波电压幅值，I_pcc-h(m)表示第m次测量采集的谐波电流幅值；该谐波数据序列中的每一行为其一子序列。

3.根据权利要求2所述的基于波形匹配的谐波污染责任量化方法，其特征在于，所述的步骤20)中，设置基本匹配参数包括设置滑动窗口长度L，滑动距离T及相似度设定值S。

4.根据权利要求3所述的基于波形匹配的谐波污染责任量化方法，其特征在于，所述的步骤30)具体包括：

步骤301)滑动窗口从谐波数据序列的第1组数据开始，将滑动窗口内U_pcc-h与I_pcc-h子序列分别归算到区间[-1，1]内，归一化函数为：

${\stackrel{\‾}{U}}_{pcc-h}\left(p\right)=\frac{2U_{pcc-h}\left(p\right)-U_{pcc-h-min}-U_{pcc-h-max}}{U_{pcc-h-max}-U_{pcc-h-\min }}(p=1,2,3...,M)$

${\stackrel{\‾}{I}}_{pcc-h}\left(p\right)=\frac{2I_{pcc-h}\left(p\right)-I_{pcc-h-min}-I_{pcc-h-max}}{I_{pcc-h-max}-I_{pcc-h-\min }}(p=1,2,3...,M)$

式中，表示归一化后的谐波电压幅值，U_pcc-h表示归一化之前的谐波电压幅值，U_pcc-h-max表示U_pcc-h中的最大值，U_pcc-h-min表示U_pcc-h中的最小值，表示归一化后的谐波电流幅值，I_pcc-h表示归一化之前的谐波电流幅值，I_pcc-h-max表示I_pcc-h中的最大值，I_pcc-h-min表示I_pcc-h中的最小值，M为一个L长度滑动窗口内谐波电压幅值或谐波电流幅值的数据数量；p表示位于滑动窗口内的谐波电压幅值或谐波电流幅值的序号；

步骤302)按下式计算归一化后子序列中谐波电压幅值与谐波电流幅值的欧氏距离：

$DIS=\sqrt{\underset{p=1}{\overset{M}{\Σ}}{\[{\stackrel{\‾}{U}}_{pcc-h}\left(p\right)-{\stackrel{\‾}{I}}_{pcc-h}\left(p\right)\]}^2}$

式中，DIS表示归一化后子序列的序列点与的欧氏距离；

步骤303)滑动窗口向后移动T，重复步骤301)，直至滑动窗口覆盖了整个谐波测试数据，得到K个子序列的欧氏距离。

5.根据权利要求4所述的基于波形匹配的谐波污染责任量化方法，其特征在于，所述的步骤40)具体包括：计算K个子序列的欧氏距离均值为每个子序列的谐波电压幅值与谐波电流幅值的相似度表示为：

$S\left(q\right)=\frac{\stackrel{\‾}{DIS}}{DIS\left(q\right)+\stackrel{\‾}{DIS}}(q=1,2,...K)$

式中，S(q)为第q个子序列的相似度，DIS(q)为第q个子序列的欧氏距离；

从K个子序列中筛选出相似度不小于相似度设定值S的子序列，设筛选出K'组子序列。

6.根据权利要求5所述的基于波形匹配的谐波污染责任量化方法，其特征在于，所述的步骤50)具体包括：

步骤501)对筛选出的K'组子序列利用最小二乘法估算出K'个系统侧等效谐波阻抗值，最小二乘法估算方法为：

将筛选出的K'组子序列中的每一组子序列代入下式：

${\stackrel{\‾}{U}}_{pcc-h}=Z_h\·{\stackrel{\‾}{I}}_{pcc-h}+U_S$

式中，Z_h表示该组子序列的系统侧等效谐波阻抗值，U_S表示该组子序列的系统侧背景谐波电压值；

通过下式求解出该组子序列的Z_h和U_S：

$f(Z_h,U_S)=\underset{p=1}{\overset{M}{\Σ}}{\[{\stackrel{\‾}{U}}_{pcc-h}\left(p\right)-(Z_h\·{\stackrel{\‾}{I}}_{pcc-h}(p)+U_S)\]}^2$

$\frac{\∂f(Z_h,U_S)}{\∂Z_h}=0$

$\frac{\∂f(Z_h,U_S)}{\∂U_S}=0$

步骤502)对获得的K'个Z_h值，取均值为通过下式求解出测量时间内的谐波电压责任：

$T_C\left(l\right)\%=\frac{{\stackrel{\‾}{Z}}_h\·I_{pcc-h}\left(l\right)}{U_{pcc-h}\left(l\right)}\×100\%$

式中，T_C表示用户侧谐波电压责任百分数，l＝1,2,…m。