本发明公开一种电力系统谐波检测器部署方法，包括如下三大步骤：步骤1：初始化；步骤2：优化寻优；步骤3：选出F(l)最小的θ(l)，如果θ(l)中的元素为1，则相应的母线部署谐波检测器；如果θ(l)中的元素为0，则相应的母线不部署。与现有技术相比，本发明提供的电力系统谐波检测器部署方法，方法简单，不仅可以使部署检测器的数量最少，而且还可以使可检测的冗余度最多。
1.一种电力系统谐波检测器部署方法，其特征在于，包括如下三大步骤：步骤1：初始化；
步骤3：选出F(l)最小的θ(l)，如果θ(l)中的元素为1，则相应的母线部署谐波检测器；
如果θ(l)中的元素为0，则相应的母线不部署；其中：步骤1分为如下步骤完成：
步骤1.2：初始化迭代次数k＝1，最大迭代次数kmax，klimit，没有得到更优的解的记录kcount(l)，蜂群的规模NS和待优化解的维数N；随机产生二元优化解θ(l,d)的初始值，l＝1,
2,…,NS；d＝1,2,…,N；θ(l)＝[θ(l,1),θ(l,2),…,θ(l,N)]；
步骤1.3：判断θ(l)元素中1的个数是否超过n个，如果超过，则重新随机产生θ(l)；
步骤1.4：根据θ(l)的元素判断各条母线是否有部署谐波检测器，在此基础上再根据母线的谐波检测器能否检测到其它母线的谐波信息，构造矩阵X；然后由式(1)计算每个θ(l)的代价函数F(l)；
F(l)＝ω1J1+ω2J2 (1)式(1)中，ω1和ω2为权值系数；J1＝ATA，J2＝(R-XA)T(R-XA)，J1越小表示部署的谐波检测器数量越少，J2越小表示部署的谐波检测器所能获得信息冗余度越接近期望冗余值R；
步骤2.1：工蜂寻优阶段，分为以下具体步骤：步骤2.1.1：将l设置为1；
步骤2.1.2：如果l≤NS,将执行步骤2.1.3，否则跳至步骤2.2；
对θ(l)、θ(i)、θ(j)、θbest和zero两两进行随机配对；
对配对后的θ(l)、θ(i)、θ(j)、θbest和zero两两互换它们的部分元素；
由式(1)计算θ(l)、θ(i)、θ(j)、θbest和zero的代价函数，并从中选出代价函数最小的，作为θEB(l)及相应地赋值予FEB(l)；
步骤2.2：根据式(2)计算每个θ(l)相应代价函数的概率；
步骤2.3：工蜂寻优阶段，分为以下具体步骤：步骤2.3.1：将l设置为1；
步骤2.3.2：如果l≤NS,将执行步骤2.3.3，否则跳至步骤2.4；
θOB(l)＝θ(l)+rand[0 1][θ(i)-θ(j)] (3)如果θOB(l)中的元素等于或大于0.5，则把它置为1；否则，置为0；
由式(1)计算θOB(l)、θ(i)和θ(j)的代价函数，并从中选出代价函数最小的，更新θOB(l)及其FOB(l)；
步骤2.4：侦察蜂寻优阶段，分为以下具体步骤：步骤2.4.1：将l设置为1；
步骤2.4.2：如果l≤NS,将执行步骤2.4.3，否则跳至步骤2.5；
步骤2.4.3：如果kcount(l)≥klimit,则执行步骤2.4.4，否则跳至步骤2.5；
随机产生一个二元的N维列向量，记为θSB；如果θSB元素中1的个数超过n个，则重新随机产生；
步骤2.6：如果k大于kmax，则跳至步骤3；否则，跳至步骤2.1.3。
一种电力系统谐波检测器部署方法
技术领域
[0001] 本发明涉及电力系统领域，确切地说是指一种基于二元人工蜂群优化机理的电力系统谐波检测器部署方法。
背景技术
[0002] 近年来，由于新能源并入大电网以及非线性负载的出现，使得电力系统与自然界的环境一样也会受到“污染”，即谐波干扰，谐波不仅会影响电力系统的正常工作，还会引起电力系统故障的产生。为了保证电力系统有效稳定运行，显而易见需要抑制谐波对电力系统的影响，目前采用的谐波抑制方法均须以谐波源为先知信息，但影响电力系统正常工作的谐波源极少是显性的，更多的谐波来源是未知的。
发明内容
[0003] 针对上述缺陷，本发明解决的技术问题在于提供一种电力系统谐波检测器部署方法，方法简单，不仅可以使部署检测器的数量最少，而且还可以使可检测的冗余度最多。
[0004] 为了解决以上的技术问题，本发明提供的一种电力系统谐波检测器部署方法，包括如下三大步骤：
[0005] 步骤1：初始化；
[0006] 步骤2：优化寻优；
[0007] 步骤3：选出F(l)最小的θ(l)，如果θ(l)中的元素为1，则相应的母线部署谐波检测器；如果θ(l)中的元素为0，则相应的母线不部署；其中：
[0008] 步骤1分为如下步骤完成：
[0009] 步骤1.1：给定部署谐波检测器的期望冗余值R；
[0010] 步骤1.2：初始化迭代次数k＝1，最大迭代次数kmax，klimit，没有得到更优的解的记录kcount(l)，蜂群的规模NS和待优化解的维数N；随机产生二元优化解θ(l,d)的初始值，l＝1,2,…,NS；d＝1,2,…,N；θ(l)＝[θ(l,1),θ(l,2),…,θ(l,N)]；
[0011] 步骤1.3：判断θ(l)元素中1的个数是否超过n个，如果超过，则重新随机产生θ(l)；
[0012] 步骤1.4：根据θ(l)的元素判断各条母线是否有部署谐波检测器，在此基础上再根据母线的谐波检测器能否检测到其它母线的谐波信息，构造矩阵X；然后由式(1)计算每个θ(l)的代价函数F(l)；
[0013] F(l)＝ω1J1+ω2J2 (1)
[0014] 式(1)中，ω1和ω2为权值系数；J1＝ATA，J2＝(R-XA)T(R-XA)，J1越小表示部署的谐波检测器数量越少，J2越小表示部署的谐波检测器所能获得信息冗余度越接近期望冗余值R；
[0015] 其中，步骤2分为如下步骤完成：
[0016] 步骤2.1：工蜂寻优阶段，分为以下具体步骤：
[0017] 步骤2.1.1：将l设置为1；
[0018] 步骤2.1.2：如果l≤NS,将执行步骤2.1.3，否则跳至步骤2.2；
[0019] 步骤2.1.3：
[0020] 随机产生l相邻的序号i和j，i≠j≠l；
[0021] 选出F(l)最小的θ(l)，记为θbest；
[0022] 产生一个N维的零列向量，记为zero；
[0023] 对θ(l)、θ(i)、θ(j)、θbest和zero两两进行随机配对；
[0024] 对配对后的θ(l)、θ(i)、θ(j)、θbest和zero两两互换它们的部分元素；
[0025] 由式(1)计算θ(l)、θ(i)、θ(j)、θbest和zero的代价函数，并从中选出代价函数最小的，作为θEB(l)及相应地赋值予FEB(l)；
[0026] 如果FEB(l)
[0027] 步骤2.2：根据式(2)计算每个θ(l)相应代价函数的概率；
[0028]
[0029] 步骤2.3：工蜂寻优阶段，分为以下具体步骤：
[0030] 步骤2.3.1：将l设置为1；
[0031] 步骤2.3.2：如果l≤NS,将执行步骤2.3.3，否则跳至步骤2.4；
[0032] 步骤2.3.3：如果p(l)
[0033] 步骤2.3.4：
[0034] 随机产生l相邻的序号i和j，i≠j≠l；
[0035] 由式(3)计算θOB(l)；
[0036] θOB(l)＝θ(l)+rand[0 1][θ(i)-θ(j)] (3)
[0037] 如果θOB(l)中的元素等于或大于0.5，则把它置为1；否则，置为0；
[0038] 对θ(i)和θ(j)随机地互换它们的部分元素；
[0039] 由式(1)计算θOB(l)、θ(i)和θ(j)的代价函数，并从中选出代价函数最小的，更新θOB(l)及其FOB(l)；
[0040] 如果FOB(l)
[0041] 步骤2.4：侦察蜂寻优阶段，分为以下具体步骤：
[0042] 步骤2.4.1：将l设置为1；
[0043] 步骤2.4.2：如果l≤NS,将执行步骤2.4.3，否则跳至步骤2.5；
[0044] 步骤2.4.3：如果kcount(l)≥klimit,则执行步骤2.4.4，否则跳至步骤2.5；
[0045] 步骤2.4.4：
[0046] 随机产生一个二元的N维列向量，记为θSB；如果θSB元素中1的个数超过n个，则重新随机产生；
[0047] 由式(1)计算θSB的代价函数，记为FSB；
[0048] 如果FSB
[0049] 步骤2.5：k＝k+1；
[0050] 步骤2.6：如果k大于kmax，则跳至步骤3；否则，跳至步骤2.1.3。
[0051] 在本发明中，A为N维的列向量，N为母线数量，本发明将A设定为二元人工蜂群待优化解θ；当A的第i元素为1时，表示第i条母线上部署有谐波检测器，如果为0时，表示它有部署谐波检测器；A元素中1的个数不能超过n个，n为可用谐波检测器的个数。X为N阶的矩阵，即X＝[xij]N×N，xij表示母线i的谐波检测器是否能检测到母线j的谐波信息，1表示可以，0表示不可以。
[0052] 准确辨识船舶电力系统中的未知谐波源，需要融合来源于谐波检测器的实测数据。合理的谐波检测器部署方法将为谐波源的准确辨识提供良好的前提条件。与现有技术相比，本发明的电力系统谐波检测器部署方法，方法简单，不仅可以使部署检测器的数量最少，而且还可以使可检测的冗余度最多。
具体实施方式
[0053] 为了本领域的技术人员能够更好地理解本发明所提供的技术方案，下面结合具体实施例进行阐述。
[0054] 本发明实施例提供的电力系统谐波检测器部署方法，包括如下三大步骤：
[0055] 步骤1：初始化；
[0056] 步骤2：优化寻优；
[0057] 步骤3：选出F(l)最小的θ(l)，如果θ(l)中的元素为1，则相应的母线部署谐波检测器；如果θ(l)中的元素为0，则相应的母线不部署；其中：
[0058] 步骤1分为如下步骤完成：
[0059] 步骤1.1：给定部署谐波检测器的期望冗余值R；
[0060] 步骤1.2：初始化迭代次数k＝1，最大迭代次数kmax，klimit，没有得到更优的解的记录kcount(l)，蜂群的规模NS和待优化解的维数N；随机产生二元优化解θ(l,d)的初始值，l＝1,2,…,NS；d＝1,2,…,N；θ(l)＝[θ(l,1),θ(l,2),…,θ(l,N)]；
[0061] 步骤1.3：判断θ(l)元素中1的个数是否超过n个，如果超过，则重新随机产生θ(l)；
[0062] 步骤1.4：根据θ(l)的元素判断各条母线是否有部署谐波检测器，在此基础上再根据母线的谐波检测器能否检测到其它母线的谐波信息，构造矩阵X；然后由式(1)计算每个θ(l)的代价函数F(l)；
[0063] F(l)＝ω1J1+ω2J2 (1)
[0064] 式(1)中，ω1和ω2为权值系数；J1＝ATA，J2＝(R-XA)T(R-XA)，J1越小表示部署的谐波检测器数量越少，J2越小表示部署的谐波检测器所能获得信息冗余度越接近期望冗余值R；
[0065] 其中，步骤2分为如下步骤完成：
[0066] 步骤2.1：工蜂寻优阶段，分为以下具体步骤：
[0067] 步骤2.1.1：将l设置为1；
[0068] 步骤2.1.2：如果l≤NS,将执行步骤2.1.3，否则跳至步骤2.2；
[0069] 步骤2.1.3：
[0070] 随机产生l相邻的序号i和j，i≠j≠l；
[0071] 选出F(l)最小的θ(l)，记为θbest；
[0072] 产生一个N维的零列向量，记为zero；
[0073] 对θ(l)、θ(i)、θ(j)、θbest和zero两两进行随机配对；
[0074] 对配对后的θ(l)、θ(i)、θ(j)、θbest和zero两两互换它们的部分元素；
[0075] 由式(1)计算θ(l)、θ(i)、θ(j)、θbest和zero的代价函数，并从中选出代价函数最小的，作为θEB(l)及相应地赋值予FEB(l)；
[0076] 如果FEB(l)
[0077] 步骤2.2：根据式(2)计算每个θ(l)相应代价函数的概率；
[0078]
[0079] 步骤2.3：工蜂寻优阶段，分为以下具体步骤：
[0080] 步骤2.3.1：将l设置为1；
[0081] 步骤2.3.2：如果l≤NS,将执行步骤2.3.3，否则跳至步骤2.4；
[0082] 步骤2.3.3：如果p(l)
[0083] 步骤2.3.4：
[0084] 随机产生l相邻的序号i和j，i≠j≠l；
[0085] 由式(3)计算θOB(l)；
[0086] θOB(l)＝θ(l)+rand[0 1][θ(i)-θ(j)] (3)
[0087] 如果θOB(l)中的元素等于或大于0.5，则把它置为1；否则，置为0；
[0088] 对θ(i)和θ(j)随机地互换它们的部分元素；
[0089] 由式(1)计算θOB(l)、θ(i)和θ(j)的代价函数，并从中选出代价函数最小的，更新θOB(l)及其FOB(l)；
[0090] 如果FOB(l)
[0091] 步骤2.4：侦察蜂寻优阶段，分为以下具体步骤：
[0092] 步骤2.4.1：将l设置为1；
[0093] 步骤2.4.2：如果l≤NS,将执行步骤2.4.3，否则跳至步骤2.5；
[0094] 步骤2.4.3：如果kcount(l)≥klimit,则执行步骤2.4.4，否则跳至步骤2.5；
[0095] 步骤2.4.4：
[0096] 随机产生一个二元的N维列向量，记为θSB；如果θSB元素中1的个数超过n个，则重新随机产生；
[0097] 由式(1)计算θSB的代价函数，记为FSB；
[0098] 如果FSB
[0099] 步骤2.5：k＝k+1；
[0100] 步骤2.6：如果k大于kmax，则跳至步骤3；否则，跳至步骤2.1.3。
[0101] 在本发明中，A为N维的列向量，N为母线数量，本发明将A设定为二元人工蜂群待优化解θ；当A的第i元素为1时，表示第i条母线上部署有谐波检测器，如果为0时，表示它有部署谐波检测器；A元素中1的个数不能超过n个，n为可用谐波检测器的个数。X为N阶的矩阵，即X＝[xij]N×N，xij表示母线i的谐波检测器是否能检测到母线j的谐波信息，1表示可以，0表示不可以。
[0102] 准确辨识船舶电力系统中的未知谐波源，需要融合来源于谐波检测器的实测数据。合理的谐波检测器部署方法将为谐波源的准确辨识提供良好的前提条件。与现有技术相比，本发明的电力系统谐波检测器部署方法，方法简单，不仅可以使部署检测器的数量最少，而且还可以使可检测的冗余度最多。
[0103] 对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
