申请/专利权人：西南交通大学

申请日：2017-05-18

公开（公告）日：2020-07-17

公开（公告）号：CN107203842B

主分类号：G06Q10/06(20120101)

分类号：G06Q10/06(20120101);G06Q50/26(20120101)

优先权：

专利状态码：有效-授权

法律状态：2020.07.17#授权;2017.10.27#实质审查的生效;2017.09.26#公开

摘要：本发明公开了一种基于扩展云相似度与逼近理想解的谐波污染水平评估方法，首先以电压总谐波畸变率和谐波电流有效值作为指标，建立了评估指标体系；再采用群决策特征根法整合多名专家的权重信息，得到谐波污染水平评估专家群体的综合权重；接着采用三角模糊数刻画谐波污染水平评估指标的模糊性，采用扩展云相似性测度将三角模糊数进行统一计算；最后在构造统一的基准等级和正负理想解的情况下，用逼近理想解法给出合理的谐波污染水平评估结果。本发明方法可以全面评估谐波污染水平，考虑一段时间谐波指标的波动性和评估决策的群体性，较好地整合多位专家的权重信息，给出合理的谐波污染水平评估结果。

主权项：1.一种基于扩展云相似度与逼近理想解的谐波污染水平评估方法，其特征在于，包括以下步骤：步骤1：构建谐波污染水平评估指标步骤1.1：将谐波电压总谐波畸变率和5次、7次、11次和13次谐波电流纳入评估指标集，根据国标GBT14549-1993指定4个谐波污染标准等级，将国标限值作为第2等级，其中谐波电流指标等级间距为国标值的20％，谐波电压总谐波畸变率等级间距为2％；将谐波污染水平分为4个基准等级F1,F2,F3,F4，同时设置正负理想解F+和F-，4个等级与正负理想解之间形成的区间分别代表谐波水平合格、轻度谐波污染、中度谐波污染、重度谐波污染和严重谐波污染；构造谐波污染水平决策矩阵X，4个基准等级和正负理想解看作6个方案，方案1～6依次为F+,F1,F2,F3,F4,F-；每个等级下有5个指标，指标1～5依次代表电压总谐波畸变率和5次、7次、11次和13次谐波电流，谐波污染水平决策矩阵X为： 谐波污染水平决策矩阵X中的任一元素xij为方案i在第j个指标下的属性值；步骤1.2：谐波污染水平决策矩阵X的属性值xij被刻画为三角模糊数aij,bij,cij；其中，j＝1时，a,b,c分别为电压总谐波畸变率、95％概率值和最大值；1＜j≤5时，a,b,c分别为各次谐波电流有效值的平均值、95％概率值和最大值；至此，谐波污染水平决策矩阵成为每个xij包含aij,bij,cij三个元素的矩阵；步骤2：确定谐波污染水平评估的指标权重设o个参评人员，5个评估指标依次为电压总谐波畸变率和5次、7次、11次和13次谐波电流；设评估参评人员集为E＝{E1,E2,…,Eo}；参评人员Em对各指标赋权为Wm＝[wm1,wm2,…,wm5]，其中，wmj为参评人员m对第j个指标的赋权，wmj∈[0,1],1≤m≤o，且所有参评人员的指标赋权构成了赋权矩阵W，o×5维矩阵；记理想评估参评人员为E*，代表整个参评人员群体，E*给出的各指标权重W*＝[w1,w2,…,w5]，即整合所有参评人员意见后的综合权重；令M＝WTW，因wij＞0，故M＝WTW为非负不可约矩阵；由Perron-Frobenius定理知，M的最大特征根对应的归一化正特征向量即为谐波污染水平评估的综合权重W*，W*采用数值幂法求得；步骤3：谐波污染水平决策矩阵X的规范化对第j个指标对应的第j列元素aij,bij,cij1≤i≤6，根据下式计算 式中，max[aij]，max[bij]，max[cij]为在第j个指标下模糊数各属性值的最大值，分别为规范化后的决策矩阵模糊属性值；步骤4：将模糊属性值转换为三角云三角云表示为YEx,En1,En2,He1,He2，由两个正态半云Ex,En1,He1和Ex,En2,He2组成；根据下式将x＝a,b,c转化为三角云： 式中，Ex为三角云的期望，En1和En2为三角云的熵，He1和He2为三角云的超熵，k为常数；至此，决策矩阵X中各元素xij对应三角云决策矩阵Y中Yij＝{Exij,En1ij,En2ij,He1ij,He2ij}五个元素；步骤5：计算每个属性下的理想解正理想解Y+表示谐波无污染的理想情况，负理想解Y-表示谐波污染最严重的理想情况，正负理想解均由各指标下的正负理想解组成；Y+＝[Y1+,Y2+,...Y5+]Y-＝[Y1-,Y2-,...Y5-]其中，分别为第j个指标的正、负理想解；正负理想解的超熵分别设为0和0.02，各指标下的正、负理想解通过下式计算： 其中，max[Exij],max[En1ij],max[En2ij]分别为三角云决策矩阵第j列元素中Exij,En1ij,En2ij的最大值；min[Exij],min[En1ij],min[En2ij]分别为三角云决策矩阵第j列元素中Exij,En1ij,En2ij的最小值；步骤6：计算各等级与正理想解的贴近度步骤6.1：分别计算三角云决策矩阵Y中各元素Yij{Exij,En1ij,En2ij,He1ij,He2ij}与该指标下正理想解的相似度和负理想解的相似度 所述步骤6.1中计算扩展云相似度计算过程为：对两个三角云YEx,En1,En2,He1,He2和Y′Ex′,En1′,En2′,He1′,He2′，向量和分别表示YEx,En1,En2,He1,He2的上半云和下半云；相应的，和为YEx1′,En1′,En2′,He1′,He2′的上半云和下半云；两个三角云YEx,En1,En2,He1,He2和Y′Ex′,En1′,En2′,He1′,He2′，二者之间的相似度为： 式中，α是调节系数；步骤6.2：带入综合权重W*，依据下式计算各基准等级与正负理想解的距离和 步骤6.3：依据下式计算各基准等级与正理想解的贴近度Ri； 步骤7：给出评估值区间根据评估值贴近度给出表征谐波水平合格、轻度谐波污染、中度谐波污染、重度谐波污染和严重谐波污染的评估值区间；若步骤8：评估母线谐波污染等级输入待评估母线的谐波污染水平评估指标，每条母线均输入5个指标，指标1～5依次代表电压总谐波畸变率和5次、7次、11次和13次谐波电流，构建母线谐波污染水平评估决策矩阵 母线谐波污染水平评估决策矩阵G中任一元素gnj为母线n在第j个指标下的属性值；gnj包含三个元素，j＝1时，p,q,l分别为电压总谐波畸变率、95％概率值和最大值；1＜j≤5时，p,q,l分别为各次谐波电流有效值的平均值、95％概率值和最大值；至此，母线谐波污染水平评估决策矩阵G每个元素对应三个数值；对母线谐波污染水平评估决策矩阵G，采用步骤3进行决策矩阵规范化；对规范化决策矩阵，采用步骤4转化为三角云决策矩阵；对三角云决策矩阵，利用步骤5求得的正负理想解，采用步骤6求各母线与正理想解的贴近度Rn；根据步骤7得到的评估值区间划分各母线的谐波污染水平。

全文数据：基于扩展云相似度与逼近理想解的谐波污染水平评估方法技术领域[0001]本发明属于电力系统领域，特别涉及一种基于扩展云相似度与逼近理想解的谐波污染水平评估方法。背景技术[0002]近年来，随着电力电子设备等非线性设备的大量应用，电网谐波污染加剧，谐波已成为电能质量中的突出问题。为分析母线的谐波污染是否超标，有必要利用谐波量测数据和谐波源定位结果对母线或区域的谐波污染水平进行综合评估。[0003]目前，鲜有谐波污染水平评估的研究报道，但已有大量的电能质量综合评估方法，且较为成熟，为谐波污染水平评估提供一定参考。电能质量评估的重点是指标的定量表达和不同指标的综合评价。矢量代数法将各分项指标归一化、层次分析法计算主观权重、物元分析法计算不同等级关联度、模糊数学法定量表达指标的模糊性，这些方法具有一定片面性。神经网络法和遗传算法则无需主观赋权，但原理复杂、实现困难。同时，上述评估方法仅以谐波电压作为衡量谐波污染水平的指标，未能较全面地反映谐波污染水平。[0004]因此，有必要综合考虑谐波污染因素，构建一个谐波污染水平评估指标集。同时，考虑一段时间内谐波参数的变化，提出一种谐波污染水平综合评估方法，为供电部门和用户提供一个全面的谐波污染评估结果。发明内容[0005]本发明所要解决的技术问题是提供一种基于扩展云相似度与逼近理想解的谐波污染水平评估方法，构建谐波污染水平评估指标集，进行一段时间内谐波污染水平综合评估，得到全面的谐波污染评估结果。[0006]为解决上述技术问题，本发明采用的技术方案是：[0007]—种基于扩展云相似度与逼近理想解的谐波污染水平评估方法，包括以下步骤：[0008]步骤1:构建谐波污染水平评估指标[0009]步骤1.1:将谐波电压总谐波畸变率和5次、7次、11次和13次谐波电流纳入评估指标集，根据国标GBT14549-1993指定4个谐波污染标准等级，将国标限值作为第2等级，其中谐波电流指标等级间距为国标值的20%，电压总谐波畸变率等级间距为2%;将谐波污染水平分为4个基准等级F1J2,F3，F4,同时设置正负理想解F+和F，4个等级与正负理想解之间形成的区间分别代表谐波水平合格、轻度谐波污染、中度谐波污染、重度谐波污染和严重谐波污染；[0010]构造谐波污染水平决策矩阵X，4个基准等级和正负理想解看作6个方案，方案1〜6依次为F+，F1，F2，F3，F4，F;每个等级下有5个指标，指标1〜5依次代表电压总谐波畸变率和5次、7次、11次和13次谐波电流，谐波污染水平决策矩阵X为：[0011][0012]谐波污染水平决策矩阵X中的任一元素X1伪方案i在第j个指标下的属性值；[0013]步骤1.2:谐波污染水平决策矩阵X的属性值Xlj被刻画为三角模糊数au，1¾，ClJ;其中，j=1时，a，b，c分别为电压总谐波畸变率、95%概率值和最大值；I0,故M=WtW为非负不可约矩阵；由Perron-Frobenius定理知，M的最大特征根对应的归一化正特征向量即为谐波污染水平评估的综合权重W'采用数值舉法求得；[0018]步骤3:谐波污染水平决策矩阵X的规范化[0019]对第j个指标对应的第j列元素叫，blj，Clj1纪0，故M=WtW为非负不可约矩阵。由Perron-Frobenius定理可知，M的最大特征根对应的归一化正特征向量即为谐波污染水平评估的综合权重W'可采用数值舉法求得。[0077]步骤3:谐波污染水平决策矩阵X的规范化[0078]根据评估指标的成本型特性，对决策矩阵进行属性值的规范化，对第j个指标对应的第j列元素a„，b„，C„1彡i6，根据下式计算[0079][0080]式中，max[aij]，max[bij]，max[cij]为在第j个指标下模糊数各属性值的最大值，分别为规范化后的决策矩阵模糊属性值。[0081]步骤4:将模糊属性值转换为三角云[0082]为将模糊数进行统一计算，将三角模糊数转换三角云。三角云可表示为YEx，Em，En2，Hei，He2，由两个正态半云Ex，Eni，Hei和Εχ，Επ2，Ηθ2组成。根据下式将X=a，b，c转化为三角云：[0083][0084]式中，Ex为三角云的期望，Em和En2为三角云的熵，HejPHe2为三角云的超熵，k为常数。至此，决策矩阵X中各元素Xij对应三角云决策矩阵Y中Yij={Exij,Eniij,Ermj,Heiij,He2ij}五个元素。[0085]步骤5:计算每个属性下的理想解[0086]正理想解Y+表示谐波无污染的理想情况，负理想解Γ表示谐波污染最严重的理想情况，正负理想解均由各指标下的正负理想解组成。[0087]Y+=[Yi+,Y2+,...Ys+][0088]Y—=[ΥΓ，Υ2—，···Υ5—][0089]其中，Υ，ΥΓ1彡衫5分别为第j个指标的正、负理想解。[0090]正负理想解的超熵分别设为0和0.02，各指标下的正、负理想解通过下式计算：[0093]其中，max[Exij]，max[Eniij]，max[En2ij]分别为三角云决策矩阵第j列元素中的Exij,Eniij，En2ij的最大值;min[Exij]，min[Eniij]，min[En2ij]分别为三角云决策矩阵第j列元素中的Exij,Eniij，En2ij的最小值。[0094]步骤6:计算各等级与正理想解的贴近度[0095]步骤6·1:分别计算三角云决策矩阵Y中各元素Yij{EXij，Enlij，En2ij，Helij，He2lJ与该指标下正理想解的相似度和负理想解*的相似度[0096][0097][0098]步骤6.2:带入综合权重W*，依据下式计算各基准等级与正负理想解的距离di+和1Γ；[0099][0100][0101]步骤6.3:依据下式计算各基准等级与正理想解的贴近度Ri。[0102][0103]步骤7:给出评估值区间[0104]根据评估值贴近度给出表征谐波水平合格、轻度谐波污染、中度谐波污染、重度谐波污染和严重谐波污染的评估值区间。[0105][0106]步骤8:评估母线谐波污染等级[0107]输入待评估母线的谐波污染水平评估指标，每条母线均输入5个指标，指标1〜5依次代表电压总谐波畸变率和5次、7次、11次和13次谐波电流，构建母线谐波污染水平评估决策矩阵[0108][0109]母线谐波污染水平评估决策矩阵G中任一元素gnj为母线η在第j个指标下的属性值;gnj包含三个元素，j=1时，P，q，1分别为电压总谐波畸变率、95%概率值和最大值；Ij5时，p，q，1分别为各次谐波电流有效值的平均值、95%概率值和最大值;至此，母线谐波污染水平评估决策矩阵G每个元素对应三个数值。[0110]对母线谐波污染水平评估决策矩阵G，采用步骤3进行决策矩阵规范化；[0111]对规范化决策矩阵，采用步骤4转化为三角云决策矩阵；[0112]对三角云决策矩阵，利用步骤5求得的正负理想解，采用步骤6求各母线与正理想解的贴近度Rn。[0113]根据步骤7得到的评估值区间划分各母线的谐波污染水平。[0114]其中，所述步骤2中数值幂法求取综合权重W*的具体操作是：[0115]步骤2·1:迭代计数k=0，y。=[In，In，…，In]TeEn;[0116]步骤2.2:令yi=Myo，zi=yi||yi|!2;[0117]步骤2.3:更新让=1^+1，5^+1=]\^，21{+1=5^+1||5^+1||2;[0118]步骤2.4:求，若其满足精度要求εζε，ε为迭代精度，则zk+1为谐波污染水平评估的综合权重W*，否则转向上一步。[0119]所述步骤步骤6.1中计算扩展云相似度的具体为：[0120]对两个三角云YEx,Em,En2,Hei,He2和Y'Ex',Eni',En2',He",He2O，向量丨和丨分别表示YEx，Em,Em,He1,He2的上半云和下半云，相应地，为Y'Ex,En,Er^,He，He，）的上半云和下半云。[0121]两个三角云YEx，Em，En2，He1，He2和Y'Ex'Em人En2'He1'He2二者之间的相似度为：[0122][0123]式中，α是调节系数。[0124]为验证本发明一种基于扩展云相似度与逼近理想解的谐波污染水平评估方法，以某市电网中7条IOkV母线Busl〜Bus7A相的电压总谐波畸变率和谐波电流数据为例进行仿真，测量数据如表1所示，指标属性值表示为三角模糊数。[0125]表1各母线的指标模糊属性值[0126][0127]据统计，该电网IOkV母线的谐波电流最大值可达80A。为适应绝大多数工况，将仿真中的负理想解值统一设置为100。将各基准等级和正负理想解的指标属性值表示为三角模糊数，如表2所示。[0128]表2各基准等级和理想解的模糊属性值[0129][0130]根据本发明方法，得到表征谐波水平合格、轻度谐波污染、中度谐波污染、重度谐波污染和严重谐波污染的评估值区间如表3所示。[0131]表3各基准等级的谐波污染水平评估值[0132][0133]根据表1所示的待评估母线的评估指标属性值，形成母线谐波污染水平评估决策矩阵，采用本发明方法，求得各母线谐波污染水平与正理想解的贴近度，即母线谐波污染水平评估值：[0134]Rn=[0.28710.43810.30030.24120.26400.38060.3917][0135]各母线的谐波污染水平评估值和基准等级的关系如图2所示。由评估值可知，母线2、母线6和母线7的谐波水平合格，母线1、母线3和母线5属于轻度谐波污染，而母线4的谐波水平达到中度谐波污染。采用24小时谐波数据的平均值、95%概率值和最大值进行模糊评估，划定的谐波污染水平等级是合理的，这是因为24小时的谐波参数是变化的，采用三角模糊数、扩展云和逼近理想解法进行阶段性评估有效适应了这些变化。例如，母线5的电压总谐波畸变率，母线1、母线3和母线4的5次谐波电流，以及母线4的7次谐波电流，虽然它们的95%概率值和最大值较大，但其24小时的平均值较小。

权利要求：1.一种基于扩展云相似度与逼近理想解的谐波污染水平评估方法，其特征在于，包括以下步骤：步骤1:构建谐波污染水平评估指标步骤1.1:将谐波电压总谐波畸变率和5次、7次、11次和13次谐波电流纳入评估指标集，根据国标GBT14549-1993指定4个谐波污染标准等级，将国标限值作为第2等级，其中谐波电流指标等级间距为国标值的20%，电压总谐波畸变率等级间距为2%;将谐波污染水平分为4个基准等级F1，F2，F3，F4，同时设置正负理想解F+和Γ，4个等级与正负理想解之间形成的区间分别代表谐波水平合格、轻度谐波污染、中度谐波污染、重度谐波污染和严重谐波污染；构造谐波污染水平决策矩阵X，4个基准等级和正负理想解看作6个方案，方案1〜6依次为F+，F1，F2，F3，F4，Γ;每个等级下有5个指标，指标1〜5依次代表电压总谐波畸变率和5次、7次、11次和13次谐波电流，谐波污染水平决策矩阵X为：谐波污染水平决策矩阵X中的任一元素为方案i在第j个指标下的属性值；步骤1.2:谐波污染水平决策矩阵X的属性值Xij被刻画为三角模糊数aij，bij，Cij;其中，j=1时，a，b，c分别为电压总谐波畸变率、95%概率值和最大值；I0，故M=WtW为非负不可约矩阵；由Perron-Frobenius定理知，M的最大特征根对应的归一化正特征向量即为谐波污染水平评估的综合权重W'采用数值幂法求得；步骤3:谐波污染水平决策矩阵X的规范化对第j个指标对应的第j列元素a^b^cdli6，根据下式计算式中，max[aij]，max[bij]，max[cij]为在第j个指标下模糊数各属性值的最大值，分别为规范化后的决策矩阵模糊属性值；步骤4:将模糊属性值转换为三角云三角云表示为YEx，Em，Em，Hei，He2，由两个正态半云（Ex，Em，Hei和（Ex，Em，He2组成;根据下式将X=a，b，c转化为三角云：式中，Ex为三角云的期望，Em和En2为三角云的熵，He^He2为三角云的超熵，k为常数；至此，决策矩阵X中各元素Xij对应三角云决策矩阵Y中Yij={Exij，Eniij，Ermj，Heiij，He2ij}五个元素；步骤5:计算每个属性下的理想解正理想解Y+表示谐波无污染的理想情况，负理想解Y_表示谐波污染最严重的理想情况，正负理想解均由各指标下的正负理想解组成；其中分别为第j个指标的正、负理想解；正负理想解的超熵分别设为O和0.02，各指标下的正、负理想解通过下式计算：其中，max[Exij]，max[Eniij]，max[En2ij]分别为三角云决策矩阵第j列元素中Exij，Eniij,En2ij的最大值;min[Exij],min[Eniij],min[En2ij]分别为三角云决策矩阵第j列元素中Exij,Eniij，En2ij的最小值；步骤6:计算各等级与正理想解的贴近度步骤6.1:分别计算三角云决策矩阵Y中各元素Yij{Exij,Eniij,En2ij，Heiij，He2ij}与该指标下正理想f•的相似度和负理想解的相似度》^;步骤6.2:带入综合权重W'依据下式计算各基准等级与正负理想解的距离:和;步骤6.3:依据下式计算各基准等级与正理想解的贴近度R1;步骤7:给出评估值区间根据评估值贴近度给出表征谐波水平合格、轻度谐波污染、中度谐波污染、重度谐波污染和严重谐波污染的评估值区间；步骤8:评估母线谐波污染等级输入待评估母线的谐波污染水平评估指标，每条母线均输入5个指标，指标1〜5依次代表电压总谐波畸变率和5次、7次、11次和13次谐波电流，构建母线谐波污染水平评估决策矩阵母线谐波污染水平评估决策矩阵G中任一元素gnj为母线η在第j个指标下的属性值;gnj包含三个元素，j=1时，P，q，1分别为电压总谐波畸变率、95%概率值和最大值;Ij5时，P，q，1分别为各次谐波电流有效值的平均值、95%概率值和最大值;至此，母线谐波污染水平评估决策矩阵G每个元素对应三个数值；对母线谐波污染水平评估决策矩阵G，采用步骤3进行决策矩阵规范化;对规范化决策矩阵，采用步骤4转化为三角云决策矩阵；对三角云决策矩阵，利用步骤5求得的正负理想解，采用步骤6求各母线与正理想解的贴近度Rn;根据步骤7得到的评估值区间划分各母线的谐波污染水平。2.如权利要求1所述的基于扩展云相似度与逼近理想解的谐波污染水平评估方法，其特征在于，所述步骤2中数值幂法求取综合权重f的具体为：步骤2.1:迭代计数k=0，y〇=[1n,1n,…，ln]TGEn;步骤2.2令yi=Myo，zi=yi||yi|!2;步骤2.3:更亲万1^=1^+1，5^+1=]\^1{，21{+1=5^+1||5^+1||2;步骤2.4::，若其满足精度要求εζε，ε为迭代精度，则zk+1为谐波污染水平评估的综合权重W*，否则转向步骤2.3。3.如权利要求1所述的基于扩展云相似度与逼近理想解的谐波污染水平评估方法，其特征在于，所述步骤6.1中计算扩展云相似度的具为：对两个三角云YExJnhEnhHehHe2和Y'Ex',En1',En,He1',He，向量分别表示¥^，£]113112，1161，!162的上半云和下半云；相应的，J为Y'Ex,En,Er^,He，He，）的上半云和下半云；两个三角云YEx，Em，En2，He1，He2和Y'Ex'Em人En2'，He1人He2〇，二者之间的相似度为：式中，α是调节系数。

百度查询： 西南交通大学 基于扩展云相似度与逼近理想解的谐波污染水平评估方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD