[发明专利]基于幅值波动量峰度检测原理的谐波责任量化方法

权利要求书

1.一种基于幅值波动量峰度检测原理的谐波责任量化方法，其特征在于，该方法包括以下步骤：

步骤10)采集谐波实测数据，形成谐波幅值数据序列，所述谐波实测数据包括采集谐波电压幅值和谐波电流幅值数据；

步骤20)设置滑动窗口参数，所述滑动窗口W的参数包括滑动窗口长度N和滑动距离D；

步骤30)计算峰度值并筛选子序列：对所述步骤10)形成的谐波幅值数据序列按滑动窗口顺序计算谐波电压幅值与谐波电流幅值的峰度值，筛选出峰度值不等于期望值的子序列，具体包括：

步骤301)滑动窗口W从谐波幅值数据序列的第1个元素开始，针对滑动窗口W覆盖的谐波电压幅值数据序列x_U＝{x_U(n):n＝1,2,……,N}和谐波电流幅值数据序列x_I＝{x_I(n):n＝1,2,……,N}，利用下式分别计算该谐波电压幅值数据序列和谐波电流幅值数据序列的无偏估计G_U与G_I：

式中，G_U与G_I为x_U和x_I的无偏估计，x_U(n)和x_I(n)分别为x_U和x_I中的第n个元素，和分别是x_U的平均值和标准差；和分别是x_I的平均值和标准差；

步骤302)根据下式确定该滑动窗口覆盖的x_U和x_I的峰度值K_u-U和K_u-I：

式中，K_u-U和K_u-I分别为x_U和x_I的峰度值，E(G_U)和E(G_I)为G_U和G_I的期望值；var(G_U)和var(G_I)为G_U和G_I的方差；q_U和q_I为G_U和G_I的置信区间，上述参数取值为：

步骤303)判断滑动窗口W是否覆盖了整个监测时段，若不是则该滑动窗口W向后移动滑动距离D，返回步骤301)，若是则筛选出K_u-U≠E(G_U)且K_u-I≠E(G_I)所对应W的子序列；

步骤40)通过波动量法估算筛选出的子序列的系统侧等效谐波阻抗，并计算测量时段内的谐波电压责任。

2.根据权利要求1所述的基于幅值波动量峰度检测原理的谐波责任量化方法，其特征在于，所述的步骤10)具体包括：利用谐波监测装置，采集公共连接点在监测时段内的谐波电压幅值和谐波电流幅值两个谐波数据，并将两个谐波实测数据形成如下式所示的谐波幅值数据序列：

式中，U_PCC(1)和I_PCC(1)分别表示序列中第1个谐波电压幅值和谐波电流幅值数据；U_PCC(2)和I_PCC(2)分别表示序列中第2个谐波电压幅值和谐波电流幅值数据；U_PCC(m)和I_PCC(m)分别表示序列中第m个谐波电压幅值和谐波电流幅值数据；m表示谐波幅值数据序列长度。

3.根据权利要求2所述的基于幅值波动量峰度检测原理的谐波责任量化方法，其特征在于，所述步骤40)中根据下式确定系统侧等效谐波阻抗：

式中，|Z_S|cosθ为系统侧等效谐波阻抗；△U_PCC和△I_PCC分别为公共连接点谐波电压幅值和谐波电流幅值的波动量。

4.根据权利要求3所述的基于幅值波动量峰度检测原理的谐波责任量化方法，其特征在于，所述步骤40)中根据下式计算测量时段内的谐波电压责任：

式中，T_C为谐波电压责任。