(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202211119913.6 (22)申请日 2022.09.15 (71)申请人 四川大学 地址 610064 四川省成 都市一环路南 一段 24号 (72)发明人 周子琛　罗梁骁　李振霖　晁苗苗　 唐文楚　王杨　 (74)专利代理 机构 西安正华恒远知识产权代理 事务所(普通 合伙) 61271 专利代理师 黄鑫 (51)Int.Cl. G06Q 10/06(2012.01) G06Q 50/06(2012.01) (54)发明名称 一种谐波污染分区方法 (57)摘要 本发明公开了一种谐波污染分区方法， 包括 以下步骤： S1、 采集配电网各节点的谐波电压数 据， 并对谐波电压数据进行填充， 得到谐波电压 序列； S2、 根据谐波电压序列， 计算配电网两节点 间的关联度矩阵； S3、 根据配电网两节点间的关 联度矩阵， 计算平均值和节点数期望值； S4、 根据 平均值和节点数期望值， 对谐波污染进行分区， 并计算谐波污染分区评价值， 得到最佳谐波污染 分区结果； 本发 明解决了 现有谐波污染治理方法 存在谐波污染分区精度不高的问题。 权利要求书2页 说明书5页 附图1页 CN 115204757 A 2022.10.18 CN 115204757 A 1.一种谐波污染分区方法， 其特 征在于， 包括以下步骤： S1、 采集配电网各节点的谐波电压数据， 并对谐波电压数据进行填充， 得到谐波电压序 列； S2、 根据谐波电压序列， 计算配电网两节点间的关联度矩阵； S3、 根据配电网两节点间的关联度矩阵， 计算平均值和节点数期望值； S4、 根据平均值和节点数期望值， 对谐波污染进行分区， 并计算谐波污染分区评价值， 得到最佳谐波污染分区结果。 2.根据权利要求1所述的谐波污染分区方法， 其特征在于， 所述S1中对谐波电压数据进 行填充的方法为： 根据采样时间去查找谐波电压数据中的丢失点， 取丢失点左右相邻采样 时间点的谐波电压值的平均值作为丢失点的谐波电压值。 3.根据权利要求2所述的谐波污染分区方法， 其特征在于， 所述S1中进行填充的公式 为： 其中， 为第 个采样时间点填充的谐波电压值， 为第 个采样时间点的谐波 电压值， 为第 个采样时间点的谐波电压值。 4.根据权利要求1所述的谐波污染分区方法， 其特征在于， 所述S2中计算配电网两节点 间的关联度矩阵的公式为： 其中， 为配电网两节点间的关联度矩阵， 为第 个节点的谐波电压序列， 为第 个节点的谐波电压序列， 为第 个节点与第 个节点间的距离， 为配电网中节点的数量， 为二范数运算。 5.根据权利要求1所述的谐波污染分区方法， 其特征在于， 所述S3中计算平均值的公式 为： 其中， 为将配电网两节点间的关联度矩阵 的每一行元素升 序排列， 为升序排列后 的矩阵， 为升序排列后 的矩阵中第 行 第 列的元素， 为第 列所有元 素的平均值， 为配电网中节点的数量。 6.根据权利要求5所述的谐波污染分区方法， 其特征在于， 所述S3中计算节点数期望的 公式为：权　利　要　求　书 1/2 页 2 CN 115204757 A 2其中， 为第 个节点以平均值 为半径的邻域中节点数量， 为配电网中节点的 数量， 为节点数期望 。 7.根据权利要求6所述的谐波污染分区方法， 其特 征在于， 所述S4包括以下分步骤： S41、 令 =1； S42、 从配电网中选取一节点， 作为 参考节点； S43、 判断将参考节点作为中心， 以平均值 为半径的圆内是否存在不少于节点 数期望 的节点， 若是， 则该参考节点和其邻域内节点为一个谐波污染分区， 并跳转至步骤S44， 若否， 则该参 考点为噪声点， 并跳转至步骤S4 4； S44、 从配电网中选取下一节点， 作为参考节点， 并跳转至步骤S43， 直到配电网中所有 节点遍历完； S45、 自加1， 并跳转至步骤S42， 直到所有平均值 和节点数期望 遍历完； S46、 计算谐波污染分区评价值， 从经过步骤S41~S45得到的所有谐波污染分区选取最 佳谐波污染分区结果。 8.根据权利要求7所述的谐波污染分区方法， 其特征在于， 所述S46中计算谐波污染分 区评价值公式为： 其中， 为第 个节点的谐波污染分区评价值， 为第 个节点与其所在谐波 污染分区内其他节点的平均值， 为第 个节点与其他谐波污染分 区内节点的平均值 中的最小值， 为选取 和 的较大值。权　利　要　求　书 2/2 页 3 CN 115204757 A 3