一种多源谐波责任划分方法
【技术领域】
[0001] 本发明涉及谐波责任划分方法技术领域，更为具体地说，涉及一种多源谐波责任 划分方法。
【背景技术】
[0002] 由于电力电子技术的发展，整流设备、电弧炉、电气化铁路等大容量的非线性负荷 日益增加，使得电网中的谐波污染情况日趋严重，谐波含量不断增加，并产生谐波污染。而 谐波污染所引起的谐波电流，在输电线路上会产生附加能耗，增加线损;线路的谐波阻抗因 高次谐波、间谐波的存在而增加;谐波电流使变压器的铜损增大，引起变压器震荡和发热； 还会影响电能计量的准确度;谐波电压使变压器的磁滞及涡流损耗增加，增大了绝缘材料 承受的电气应力，必然减少变压器的使用寿命，如果通讯线路和电力线发生耦合，谐波会严 重影响通讯质量。电网电能质量问题不仅会影响到电力系统本身，而且还影响到电力用户 和相关行业等。随着计算机、交换机、网络设备、数控设备、半导体制造设备、医疗设备等对 电能质量敏感的用电设备的大量应用，人们对电能质量的要求也越来越高。
[0003] 电能质量污染危害用户的用电设备，影响电力用户生产的产品质量，给用户造成 重大损失。由于缺少有效的电能质量监测手段和较高的检测精度，使得电网运行存在很大 的安全隐患，电能质量也很难使用户满意。这就需要电力部门实时准确掌握电能质量情况， 并及时采取相应措施确保电能质量符合相关标准的要求。快速、准确的谐波辨识技术对提 高电能质量，保证电网和电气设备及各种用电器具的安全、稳定、经济运行，保障用户的正 常生产和生活具有重要意义。
[0004] 目前，谐波阻抗估算方法有二元线性回归法、稳健回归法、偏最小二乘法。其中，二 元线性回归法能反映出当用户侧与系统侧谐波扰动同时作用时，用户侧和系统侧的谐波发 射水平，缺点是该方法利用最小二乘估计，容易受异常数据的干扰，回归方程缺乏稳健性， 无法计算系统谐波复阻抗中的电阻分量，同时对数据统计特性的一致性要求较高。稳健回 归法在二元线性回归法的基础上通过加权处理，能够有效去除一些奇异值，但是在统计特 性上没有对变量的相关性进行分析，在数据处理上仍然会带来一些误差。偏最小二乘法能 较好地克服变量多重相关性在系统建模过程中的不良影响，但在实际运用中发现，它不是 对于所有的多重共线性问题都有效。因此，在实际运用中，上述谐波阻抗估算方法均会出现 误判的情况，从而不能准确的对谐波阻抗进行估算。
【发明内容】
[0005] 本发明的目的是提供一种多源谐波责任划分方法，以解决【背景技术】所述的现有谐 波阻抗估算方法均会出现误判的问题。
[0006] 为了解决上述技术问题，本发明提供如下技术方案：
[0007] -种多源谐波责任划分方法，所述划分方法包括以下步骤：
[0008] S01:对谐波数据进行不间断的连续采样；
[0009] S02:对每3小时所采集的数据进行PCC处谐波指标因子h计算；
[0010] S03:根据谐波指标因子^设定论域L;
[0011] S04:选择谐波责任评价方法；
[0012] S05:根据论域L中的数据特点，对步骤S04中的谐波责任评价方法进行动态赋权 Wi;
[0013] S06:根据不同谐波责任评价方法对各馈线谐波责任进行评价；
[0014] S07:根据谐波责任评价结果对PCC处各馈线的谐波污染责任进行排序，形成线性 序Ψ?;
[0015] S08:判断当前责任评价是否达到8次；
[0016] S09:若当前责任评价达到8次，则计算论域L中元素的Borda数；
[0017] S10:选定最大的Borda数所对应的馈线为此次污染责任最大的谐波源；
[0018] S11:算法结束。
[0019] 优选地，步骤S02中所述各谐波指标因子匕为： 其中是谐波源 f 负荷i产生谐波电压，是谐波畸变电压，贤是^^与的之间的夹角。
[0020] 优选地，步骤S03中所述论域L为:1^={11，12，...山}，其中，11为谐波指标因子的个 数。
[0021] 优选地，步骤S04中所述谐波责任评价方法包括二元线性回归法、稳健回归法和偏 最小二乘法。
[0022] 优选地，步骤S05中所述动态赋权Wi包括：
' t U' m其 中，t e论域L; m表示为m种谐波责任评价方法，Rtll, t2l表示为任意两种谐波责任评价方法的 相对适应度。
[0023] 优选地，所述两种谐波责任评价方法的相对适应度Rtli,t2i为:R tli,t2i= | RItli-RIt2i ，其中RItll和Rim表示为任意时刻任意两种谐波责任评价方法对某谐波负荷的责任评价 结果。
[0024] 优选地，步骤S08中所述判断当前责任评价是否达到8次为:对每3小时所采集的数 据进行一次谐波责任评价，判断当谐波前责任评价在24小时之内是否达到8次。
[0025] 优选地，步骤S09中所述Borda数为R(t)，
其中Wi动态赋权，Ri (t)为在线性序中排在t后面的计算结果的个数，te论域L。
[0026]优选地，算法结束。
[0027]本发明提供的多源谐波责任划分方法包括以下步骤:S01:对谐波数据进行不间断 的连续采样;S02:对每3小时所采集的数据进行PCC处谐波指标因子h计算;S03:根据谐波 指标因子1:设定论域L; S04:选择谐波责任评价方法;S05:根据论域L中的数据特点，对步骤 S04中的谐波责任评价方法进行动态赋权Wl; S06:根据不同谐波责任评价方法对各馈线谐 波责任进行评价;S07:根据谐波责任评价结果对PCC处各馈线的谐波污染责任进行排序，形 成线性序Ψ?; S08:判断当前责任评价是否达到8次；S09:若当前责任评价达到8次，则计算 论域L中元素的Borda数;S10:选定最大的Borda数所对应的馈线为此次污染责任最大的谐 波源;S11算法结束。本发明提供的多源谐波责任划分方法是基于模糊意见集中决策法的责 任划分方法，并根据谐波源样本数据的特点尽可能的综合二元线性回归法、稳健回归法和 偏最小二乘法，并为之赋予不同的评价权重，不仅能够准确掌握问题主体部分的信息，同时 又不会忽略对非主体信息的充分挖掘，使得评价结果准确合理。
【附图说明】
[0028] 为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使 用的附图作简单地介绍，显而易见地，对于本领域普通技术人员而言，在不付出创造性劳动 的前提下，还可以根据这些附图获得其它的附图。
[0029] 图1是本发明实施例提供的多源谐波责任划分方法的流程示意图。
【具体实施方式】
[0030] 本发明实施例提供的多源谐波责任划分方法，解决了现有谐波阻抗估算方法均会 出现误判的问题。
[0031] 为了使本技术领域的人员更好地理解本发明实施例中的技术方案，并使本发明实 施例的上述目的、特征和优点能够更加明显易懂，下面结合附图对本发明实施例中的技术 方案作进一步详细的说明。
[0032]模糊意见集中决策方法是重要的模糊数学方法，模糊决策方法有意见集中法、二 元对比法、综合评判法等，本发明采用意见集中法中的加权Borda法。加权Borda法基本思想 是:设论域U= {ul，u2，···，un}是参与决策评价的集合，即将U中元素进行排序，专家组M=m 人，发表m种意见，记为V= {vl，v2，···，vm}。其中vi是第i种意见序列，即U中元素的某一排 序。令ueU，Bi(u)表示第i种意见序列vi中排在U之后的元素个数，即若u在第i种意见vi中 排在第k位，则Bi(u)=n-k。
[0033]
为u的波达Borda数。若评价意见的重要性不同，可以赋权重，再 求加权Borda数
，各论域U的所有元素可按加权Borda数的大小排成一个线 性序列，此序列就是集中意见之后的一个比较合理的意见，根据排序结果能分出论域中元 素的优劣、先后。
[0034] 请参考附图1，附图1示出了本发明实施例提供的多源谐波责任划分方法的流程示 意图。
[0035]本发明实施例提供的多源谐波责任划分方法是一种基于上述模糊意见集中决策 法的责任划分方法。本发明实施例提供的多源谐波责任划分方法包括以下步骤：
[0036] S01:对谐波数据进行不间断的连续采样；
[0037] S02:对每3小时所采集的数据进行PCC处谐波指标因子h计算；
[0038] S03:根据谐波指标因子乜设定论域L;
[0039] S04:选择谐波责任评价方法；
[0040] S05:根据论域L中的数据特点，对步骤S04中的谐波责任评价方法进行动态赋权 Wi;
[0041 ] S06:根据不同谐波责任评价方法对各馈线谐波责任进行评价；
[0042] S07:根据谐波责任评价结果对PCC处各馈线的谐波污染责任进行排序，形成线