发明名称 智能配电网高效性评估的属性区间辨识方法
摘要 本发明是一种智能配电网高效性评估的属性区间辨识方法,该方法建立智能配电网基础数据到高效性之间的三层属性测度矩阵,在指标数据标准化基础上,基于最大熵原理计算性能指标、准则指标和综合指标的属性测度区间与综合属性测度值,自动辨识智能配电网高效性的属性区间,具体步骤为:第一步:读取评估所需信息,第二步:建立综合指标、准则指标和性能指标的属性测度矩阵,第三步:指标数据标准化处理,第四步:计算性能指标、准则指标、综合指标的属性测度区间,第五步:计算高效性指标的综合属性测度,辨识出智能配电网高效性的属性区间,从而实现对智能配电网的高效性评估。
申请公布号 CN103473450A 申请公布日期 2013.12.25
申请号 CN201310393981.6 申请日期 2013.09.02
申请人 河海大学;南京河海科技有限公司;国家电网公司;江苏省电力公司南京供电公司;中国电力科学研究院;江苏省电力公司;国网重庆市电力公司 发明人 陈星莺;余昆;陈楷;李子韵;徐石明;黄建勇;王晓晶;廖迎晨;王平
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 南京苏高专利商标事务所(普通合伙) 32204 代理人 柏尚春
主权项 1.一种智能配电网高效性评估的属性区间辨识方法,其特征在于建立智能配电网基础数据到高效性之间的三层属性测度矩阵,在指标数据标准化基础上,基于最大熵原理计算性能指标、准则指标和综合指标的属性测度区间与综合属性测度值,自动辨识智能配电网高效性的属性区间,具体步骤如下:第一步:读取评估所需信息,包括电网建设和结构信息、电网运行信息、用户可靠性需求信息、分布式电源容量及运行信息、电磁-噪声监测信息、信息系统数据、配电网业务流程信息、用户互动信息;第二步:建立综合指标、准则指标和性能指标的属性测度矩阵,具体如下:评估指标的属性测度空间有序分割集用C<sub>K</sub>,m个评估指标的属性测度矩阵如下所示:<img file="FDA00003759992700011.GIF" wi="860" he="418" />其中,C<sub>1</sub>~C<sub>K</sub>属性测度空间的元素;I<sub>j</sub>为评估指标,j=1,2,…m;[a<sub>jk</sub>,b<sub>jk</sub>]为第j个指标在属性区间F上的第k个分割区间,满足a<sub>jk</sub>≤b<sub>jk</sub>,k=1,2,…,K,记A=[a<sub>jk</sub>]<sub>m×K</sub>为下界标准矩阵,B=[b<sub>jk</sub>]<sub>m×K</sub>为上界标准矩阵;第三步:指标数据标准化处理,具体如下:规定C<sub>1</sub>级评价标准中评价指标j的下界标准a<sub>j1</sub>规格化为相对隶属度<u>s</u><sub>j1</sub>=0,C<sub>K</sub>级评价标准中评价指标j下界标准a<sub>jk</sub>规格化为相对隶属度<u>s</u><sub>jk</sub>=1,通过公式<img file="FDA00003759992700012.GIF" wi="318" he="165" />确定C<sub>K</sub>级评价标准中指标j的下界标准a<sub>jk</sub>的相对隶属度,因此,下界标准矩阵A=[a<sub>jk</sub>]<sub>m×K</sub>变换为下界标准相对隶属度矩阵S=[s<sub>jk</sub>]<sub>m×K</sub>,其中:指标值越大其等级越大的指标为正指标,反之,指标值越大其等级越小的指标为负指标,同时采用下式将评价样本矩阵X=[x<sub>ij</sub>]<sub>n×m</sub>变换为对A的相对隶属度矩阵<u>F</u>=[<u>f</u><sub>ij</sub>]<sub>n×m</sub>;正指标:<maths num="0001"><![CDATA[<math><mrow><msub><munder><mi>f</mi><mo>&OverBar;</mo></munder><mi>ij</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>0</mn><mo>,</mo><msub><mi>x</mi><mi>ij</mi></msub><mo>&lt;</mo><msub><mi>a</mi><mrow><mi>j</mi><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><mfrac><mrow><msub><mi>x</mi><mi>ij</mi></msub><mo>-</mo><msub><mi>a</mi><mrow><mi>j</mi><mn>1</mn></mrow></msub></mrow><mrow><msub><mi>a</mi><mi>jk</mi></msub><mo>-</mo><msub><mi>a</mi><mrow><mi>j</mi><mn>1</mn></mrow></msub></mrow></mfrac><mo>,</mo><msub><mi>a</mi><mrow><mi>j</mi><mn>1</mn></mrow></msub><mo>&le;</mo><msub><mi>x</mi><mi>ij</mi></msub><mo>&le;</mo><msub><mi>a</mi><mi>jk</mi></msub></mtd></mtr><mtr><mtd><mn>1</mn><mo>,</mo><msub><mi>x</mi><mi>ij</mi></msub><mo>></mo><msub><mi>a</mi><mi>jk</mi></msub></mtd></mtr></mtable></mfenced><mo>;</mo></mrow></math>]]></maths>或负指标:<maths num="0002"><![CDATA[<math><mrow><msub><munder><mi>f</mi><mo>&OverBar;</mo></munder><mi>ij</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>0</mn><mo>,</mo><msub><mi>x</mi><mi>ij</mi></msub><mo>></mo><msub><mi>a</mi><mrow><mi>j</mi><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><mfrac><mrow><msub><mi>x</mi><mi>ij</mi></msub><mo>-</mo><msub><mi>a</mi><mrow><mi>j</mi><mn>1</mn></mrow></msub></mrow><mrow><msub><mi>a</mi><mi>jk</mi></msub><mo>-</mo><msub><mi>a</mi><mrow><mi>j</mi><mn>1</mn></mrow></msub></mrow></mfrac><mo>,</mo><msub><mi>a</mi><mi>jk</mi></msub><mo>&le;</mo><msub><mi>x</mi><mi>ij</mi></msub><mo>&le;</mo><msub><mi>a</mi><mrow><mi>j</mi><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><mn>1</mn><mo>,</mo><msub><mi>x</mi><mi>ij</mi></msub><mo>&lt;</mo><msub><mi>a</mi><mi>jk</mi></msub></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>同理,将B变换为上界标准相对隶属度矩阵<img file="FDA00003759992700023.GIF" wi="266" he="106" />X变换为对B的相对隶属度矩阵<maths num="0003"><![CDATA[<math><mrow><mover><mi>F</mi><mo>&OverBar;</mo></mover><mo>=</mo><msub><mrow><mo>[</mo><msub><mover><mi>f</mi><mo>&OverBar;</mo></mover><mi>ij</mi></msub><mo>]</mo></mrow><mrow><mi>n</mi><mo>&times;</mo><mi>m</mi></mrow></msub><mo>,</mo></mrow></math>]]></maths>第四步:计算性能指标、准则指标、综合指标的属性测度区间,具体如下:n个数据样本对k个评价等级的属性测度区间矩阵如下:<img file="FDA00003759992700025.GIF" wi="894" he="480" />其中:<img file="FDA00003759992700026.GIF" wi="192" he="91" />为第i个数据样本x<sub>i</sub>相对于第C<sub>K</sub>评价等级的属性测度区间;<u>μ</u><sub>ik</sub>为下界属性测度,由样本x<sub>i</sub>与A计算得到,满足<img file="FDA00003759992700027.GIF" wi="542" he="142" />为上界属性测度,由样本x<sub>i</sub>与B计算得到,满足<img file="FDA00003759992700028.GIF" wi="453" he="142" />定义下界广义加权距离<u>D</u><sub>ik</sub>表示第i个评价样本与第k级的“下界差异”,即<maths num="0004"><![CDATA[<math><mrow><msub><munder><mi>D</mi><mo>&OverBar;</mo></munder><mi>ik</mi></msub><mo>=</mo><msub><munder><mi>&mu;</mi><mo>&OverBar;</mo></munder><mi>ik</mi></msub><msub><mi>d</mi><mi>ik</mi></msub><mo>=</mo><msub><munder><mi>&mu;</mi><mo>&OverBar;</mo></munder><mi>ik</mi></msub><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>j</mi></msub><mo>|</mo><msub><munder><mi>f</mi><mo>&OverBar;</mo></munder><mi>ij</mi></msub><mo>-</mo><msub><munder><mi>s</mi><mo>&OverBar;</mo></munder><mi>jk</mi></msub><mo>|</mo><mo>)</mo></mrow></mrow></math>]]></maths>其中<img file="FDA000037599927000210.GIF" wi="458" he="142" />是评价样本的广义距离,<u>μ</u><sub>ik</sub>需要使全体待评价样本与评价标准之间的加权广义距离之和最小,即<img file="FDA000037599927000211.GIF" wi="672" he="149" />同时,将<u>μ</u><sub>ik</sub>视为第i个评价样本属于第C<sub>K</sub>等级的“下界概率”;根据最大熵原理,得到下界属性测度<u>μ</u><sub>ik</sub>的计算式如下:<maths num="0005"><![CDATA[<math><mrow><msub><munder><mi>&mu;</mi><mo>&OverBar;</mo></munder><mi>ik</mi></msub><mo>=</mo><mfrac><mrow><mi>exp</mi><mo>[</mo><mo>-</mo><mi>B</mi><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>j</mi></msub><mo>|</mo><msub><munder><mi>f</mi><mo>&OverBar;</mo></munder><mi>ij</mi></msub><mo>-</mo><msub><munder><mi>s</mi><mo>&OverBar;</mo></munder><mi>jk</mi></msub><mo>|</mo><mo>)</mo></mrow><mo>]</mo></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><mi>exp</mi><mo>[</mo><mo>-</mo><mi>B</mi><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>j</mi></msub><mo>|</mo><msub><munder><mi>f</mi><mo>&OverBar;</mo></munder><mi>ij</mi></msub><mo>-</mo><msub><munder><mi>s</mi><mo>&OverBar;</mo></munder><mi>jk</mi></msub><mo>|</mo><mo>)</mo></mrow><mo>]</mo></mrow></mfrac></mrow></math>]]></maths>其中:ω<sub>j</sub>为各项评价指标的权重,且满足<img file="FDA00003759992700032.GIF" wi="340" he="142" />;B为正的常数;同理,上界属性测度<img file="FDA00003759992700033.GIF" wi="72" he="88" />的计算式为<maths num="0006"><![CDATA[<math><mrow><msub><mover><mi>&mu;</mi><mo>&OverBar;</mo></mover><mi>jk</mi></msub><mo>=</mo><mfrac><mrow><mi>exp</mi><mo>[</mo><mo>-</mo><mi>B</mi><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>j</mi></msub><mo>|</mo><msub><mover><mi>f</mi><mo>&OverBar;</mo></mover><mi>ij</mi></msub><mo>-</mo><msub><mover><mi>s</mi><mo>&OverBar;</mo></mover><mi>jk</mi></msub><mo>|</mo><mo>)</mo></mrow><mo>]</mo></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><mi>exp</mi><mo>[</mo><mo>-</mo><mi>B</mi><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>j</mi></msub><mo>|</mo><msub><mover><mi>f</mi><mo>&OverBar;</mo></mover><mi>ij</mi></msub><mo>-</mo><msub><mover><mi>s</mi><mo>&OverBar;</mo></mover><mi>jk</mi></msub><mo>|</mo><mo>)</mo></mrow><mo>]</mo></mrow></mfrac><mo>,</mo></mrow></math>]]></maths>从而得到评估指标的属性测度区间<img file="FDA00003759992700035.GIF" wi="322" he="94" />第五步:计算高效性指标的综合属性测度,辨识出智能配电网高效性的属性区间,具体如下:首先采用<img file="FDA00003759992700036.GIF" wi="452" he="98" />计算高效性指标的均化综合属性测度,得到样本x<sub>i</sub>属于第k类的综合属性测度μ<sub>ik</sub>,然后采用加强属性识别模型对智能配电网高效性的综合属性测度属性区间识别,计算公式如下,<maths num="0007"><![CDATA[<math><mrow><msub><mi>q</mi><mi>xi</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></munderover><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>&mu;</mi><mi>im</mi></msub><mo>+</mo><mfrac><mi>a</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></mfrac><mo>)</mo></mrow></mrow></math>]]></maths>其中,a=min(μ<sub>1k</sub>,μ<sub>2k</sub>…,μ<sub>mk</sub>),由此得到第i个样本x<sub>i</sub>的分数q<sub>xi</sub>。
地址 211100 江苏省南京市江宁区佛城西路8号