发明名称 一种基于变权层次分析法的绝缘子污闪预警方法
摘要 本发明公开了一种基于变权层次分析法的绝缘子污闪预警方法,包括以下步骤:(1)对绝缘子的固有因素进行量化评价,所述的绝缘子的固有因素包括有爬电比距、绝缘子憎水性、绝缘子老化程度;(2)采用变权层次分析法计算得到绝缘子固有因素指标<i>m</i>;根据大气环境参数和气象条件计算绝缘子表面污秽程度指标<i>P<sub>f</sub></i>;(3)根据绝缘子固有因素指标<i>m</i>和绝缘子表面污秽程度指标<i>P<sub>f</sub></i>的乘积判断绝缘子污闪综合预警指标<i>R</i>,确定预警等级。本方法能够给出所在区域绝缘子的综合预警指标并发布预警等级,实现科学的绝缘子污秽闪络预警,从而有效地提高供电可靠性,降低电网的经济损失,提高电网的可靠运行能力。
申请公布号 CN104361216A 申请公布日期 2015.02.18
申请号 CN201410594395.2 申请日期 2014.10.29
申请人 国家电网公司;国网河南省电力公司电力科学研究院;重庆大学 发明人 李哲;熊小伏;梁允;王建;刘松;张小斐
分类号 G06F19/00(2011.01)I;G01R31/12(2006.01)I 主分类号 G06F19/00(2011.01)I
代理机构 郑州联科专利事务所(普通合伙) 41104 代理人 刘建芳;崔卫琴
主权项 一种基于变权层次分析法的绝缘子污闪预警方法,其特征在于:包括以下步骤:(1)、对绝缘子的固有因素进行量化评价,所述的绝缘子的固有因素包括有爬电比距、绝缘子憎水性、绝缘子老化程度,计算绝缘子爬电比距指标k<sub>1</sub>、绝缘子憎水性指标k<sub>2</sub>、绝缘子老化指标k<sub>3</sub>,得到绝缘子固有因素指标向量K=(k<sub>1</sub>,k<sub>2</sub>,k<sub>3</sub>):其中绝缘子爬电比距指标k<sub>1</sub>计算方法为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>k</mi><mn>1</mn></msub><mo>=</mo><mfrac><mrow><msub><mi>&lambda;</mi><mi>max</mi></msub><mo>-</mo><mi>&lambda;</mi></mrow><mrow><msub><mi>&lambda;</mi><mi>max</mi></msub><mo>-</mo><msub><mi>&lambda;</mi><mi>min</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000596732940000011.GIF" wi="1630" he="209" /></maths>式中,λ为绝缘子爬电比距,λ<sub>max</sub>为绝缘子最大爬电比距,λ<sub>min</sub>为绝缘子最小爬电比距;其中绝缘子憎水性指标k<sub>2</sub>是通过评分的方式确定,并对评分等级D采取归一化处理,归一化处理后的绝缘子憎水性指标k<sub>2</sub>为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>k</mi><mn>2</mn></msub><mo>=</mo><mfrac><mrow><mi>D</mi><mo>-</mo><mn>1</mn></mrow><mrow><mn>7</mn><mo>-</mo><mn>1</mn></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000596732940000012.GIF" wi="1200" he="149" /></maths>式中,D为绝缘子憎水性评分等级;其中绝缘子老化指标k<sub>3</sub>是采用威布尔分布描述的,威布尔分布的故障密度函数定义为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>f</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><msup><mrow><mi>&beta;</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>&gamma;</mi><mo>)</mo></mrow></mrow><mrow><mi>&beta;</mi><mo>-</mo><mn>1</mn></mrow></msup><mi>&alpha;</mi></mfrac><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><msup><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>&gamma;</mi><mo>)</mo></mrow><mi>&beta;</mi></msup><mi>&alpha;</mi></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000596732940000013.GIF" wi="1394" he="172" /></maths>式中,t为绝缘子运行时间,β为形状参数,α为尺度参数,γ为位置参数;绝缘子老化指标k<sub>3</sub>直接用绝缘子累计故障函数表示,即:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>k</mi><mn>3</mn></msub><mo>=</mo><mi>F</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mn>1</mn><mo>-</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><msup><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>&gamma;</mi><mo>)</mo></mrow><mi>&beta;</mi></msup><mi>&alpha;</mi></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000596732940000014.GIF" wi="1273" he="171" /></maths>式中,t为绝缘子运行时间,β为形状参数,α为尺度参数;(2)、采用变权层次分析法计算得到绝缘子固有因素指标m;具体包括以下步骤:2.1)、根据变权层次分析法对绝缘子固有因素相对重要性的两两比较,得到绝缘子固有因素指标比较矩阵A如式(5)所示,比较矩阵A的元素取值为:对角线上的元素a<sub>ii</sub>=1;非对角线上的第i行j列元素a<sub>ij</sub>表示i因素对j因素的重要程度,同等重要取1,稍微重要取3,较强重要取5,强烈重要取7,极端重要9,两相邻判断的中间值可取2、4、6、8,此时,a<sub>ji</sub>=1/a<sub>ij</sub>;<img file="FDA0000596732940000021.GIF" wi="1372" he="350" />2.2)、通过规范列平均法计算初始权重向量<img file="FDA0000596732940000022.GIF" wi="444" he="111" />即:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mi>&omega;</mi><mi>i</mi><mn>0</mn></msubsup><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mn>3</mn></munderover><msub><mi>a</mi><mi>ij</mi></msub></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>3</mn></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mn>3</mn></munderover><msub><mi>a</mi><mi>ij</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000596732940000023.GIF" wi="1313" he="380" /></maths>得到初始权重:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msubsup><mi>&omega;</mi><msub><mi>k</mi><mi>i</mi></msub><mn>0</mn></msubsup><mo>=</mo><msup><mrow><mo>(</mo><msubsup><mi>&omega;</mi><msub><mi>k</mi><mn>1</mn></msub><mn>0</mn></msubsup><mo>,</mo><msubsup><mi>&omega;</mi><msub><mi>k</mi><mn>2</mn></msub><mn>0</mn></msubsup><mo>,</mo><msubsup><mi>&omega;</mi><msub><mi>k</mi><mn>3</mn></msub><mn>0</mn></msubsup><mo>)</mo></mrow><mi>T</mi></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000596732940000024.GIF" wi="1467" he="167" /></maths>2.3)、利用变权分析法计算变权重向量<img file="FDA0000596732940000025.GIF" wi="450" he="116" />绝缘子固有因素指标任何一个值过大都会使绝缘子处于危险水平,即它们的权重应随指标值增大而增大,因此应该选择激励型变权,变权系数a>1,其变权公式为:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>&omega;</mi><msub><mi>k</mi><mi>i</mi></msub></msub><mrow><mo>(</mo><msub><mi>k</mi><mn>1</mn></msub><mo>,</mo><msub><mi>k</mi><mn>2</mn></msub><mo>,</mo><msub><mi>k</mi><mn>3</mn></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msubsup><mi>&omega;</mi><msub><mi>k</mi><mi>i</mi></msub><mn>0</mn></msubsup><msubsup><mi>k</mi><mi>i</mi><mrow><mi>a</mi><mo>-</mo><mn>1</mn></mrow></msubsup></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>3</mn></munderover><msubsup><mi>&omega;</mi><msub><mi>k</mi><mi>i</mi></msub><mn>0</mn></msubsup><msubsup><mi>k</mi><mi>i</mi><mrow><mi>a</mi><mo>-</mo><mn>1</mn></mrow></msubsup></mrow></mfrac><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2,3</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000596732940000026.GIF" wi="1686" he="304" /></maths>2.4)、计算绝缘子固有因素指标m:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mi>m</mi><mo>=</mo><mi>K</mi><mo>&times;</mo><mi>&omega;</mi><mo>=</mo><msub><mi>&omega;</mi><msub><mi>k</mi><mn>1</mn></msub></msub><msub><mi>k</mi><mn>1</mn></msub><mo>+</mo><msub><mi>&omega;</mi><msub><mi>k</mi><mn>2</mn></msub></msub><msub><mi>k</mi><mn>2</mn></msub><mo>+</mo><msub><mi>&omega;</mi><msub><mi>k</mi><mn>3</mn></msub></msub><msub><mi>k</mi><mn>3</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000596732940000027.GIF" wi="1446" he="94" /></maths>式中,m为绝缘子固有因素指标;k<sub>1</sub>为绝缘子爬电比距指标;k<sub>2</sub>为绝缘子憎水性指标;k<sub>3</sub>为绝缘子老化指标;(3)、根据环保部门提供的SO<sub>2</sub>、NO<sub>2</sub>、PM10浓度测量数据和气象局发布的气象预报信息,计算绝缘子表面污秽程度指标P<sub>f</sub>,具体包括以下步骤:3.1)、计算大气质量指数P和绝缘子等值盐密值ρ:将环保部门提供的SO<sub>2</sub>、NO<sub>2</sub>、PM10浓度测量数据代入式(10)计算大气质量指数P,然后将计算所得的P值代入等值盐密值与大气质量指数经验式(11)计算绝缘子等值盐密值ρ:<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><mi>P</mi><mo>=</mo><msqrt><mfrac><msub><mi>U</mi><mrow><mi>PM</mi><mn>10</mn></mrow></msub><msub><mi>B</mi><mrow><mi>PM</mi><mn>10</mn></mrow></msub></mfrac><mo>&times;</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><mo>|</mo><mfrac><msub><mi>U</mi><msub><mi>SO</mi><mn>2</mn></msub></msub><msub><mi>B</mi><msub><mi>SO</mi><mn>2</mn></msub></msub></mfrac><mo>+</mo><mfrac><msub><mi>U</mi><msub><mi>NO</mi><mn>2</mn></msub></msub><msub><mi>B</mi><msub><mi>NO</mi><mn>2</mn></msub></msub></mfrac><mo>+</mo><mfrac><msub><mi>U</mi><mrow><mi>PM</mi><mn>10</mn></mrow></msub><msub><mi>B</mi><mrow><mi>PM</mi><mn>10</mn></mrow></msub></mfrac><mo>|</mo></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000596732940000031.GIF" wi="1525" he="216" /></maths>ρ=‑0.015+0.0806P   (11)式中:U为污染物月度实测浓度,mg/m<sup>3</sup>;B为污染物评价标准值,mg/m<sup>3</sup>;SO<sub>2</sub>、NO<sub>2</sub>、PM10的评价标准值分别取0.06mg/m<sup>3</sup>,0.05mg/m<sup>3</sup>,0.1mg/m<sup>3</sup>;公式中的U<sub>PM10</sub>/B<sub>PM10</sub>独立项可以变换,选取区域影响最突出的因素;应用大气质量指数P预测绝缘子积污的发展趋势,控制绝缘子表面等值盐密值,合理安排清扫计划;当绝缘子的等值盐密值与允许值有较大裕度时,可以延长清扫周期;当清扫工作完成后,系统绝缘子表面等值盐密值清零,若清零时间在当月第t日,则清零后月末的等值盐密值为<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><msubsup><mi>&rho;</mi><mi>m</mi><mo>&prime;</mo></msubsup><mo>=</mo><mfrac><msub><mi>&rho;</mi><mi>m</mi></msub><mn>30</mn></mfrac><mo>&times;</mo><mrow><mo>(</mo><mn>30</mn><mo>-</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000596732940000032.GIF" wi="1348" he="141" /></maths>式中,ρ<sub>m</sub>为按大气质量指数计算的月累积等值盐密值;3.2)、确定绝缘子表面污秽程度指标P<sub>f</sub>:绝缘子表面污秽程度指标P<sub>f</sub>使用绝缘子在某等值盐密值ρ下发生闪络的概率为P<sub>f</sub>表示,则单串绝缘子在运行电压U<sub>c</sub>下,P<sub>f</sub>与等值盐密值ρ的关系可表示为:<maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><msub><mi>P</mi><mi>f</mi></msub><mo>=</mo><mn>1</mn><mo>-</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><msup><mrow><mo>(</mo><mfrac><mn>1</mn><mi>&epsiv;</mi></mfrac><mrow><mo>(</mo><mfrac><msup><mi>&rho;</mi><mi>&sigma;</mi></msup><mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>nc</mi><mo>)</mo></mrow><mi>AL</mi><mo>/</mo><msub><mi>U</mi><mi>c</mi></msub></mrow></mfrac><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>)</mo></mrow><mi>h</mi></msup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000596732940000033.GIF" wi="1533" he="216" /></maths><maths num="0012" id="cmaths0012"><math><![CDATA[<mrow><mi>&epsiv;</mi><mo>=</mo><mfrac><mi>nc</mi><mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>nc</mi><mo>)</mo></mrow><msup><mrow><mo>(</mo><mi>ln</mi><mn>2</mn><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mi>h</mi></mrow></msup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>14</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000596732940000041.GIF" wi="1288" he="155" /></maths><maths num="0013" id="cmaths0013"><math><![CDATA[<mrow><mi>h</mi><mo>=</mo><mfrac><mn>1.38</mn><mrow><mi>ln</mi><mi>n</mi><mo>-</mo><mi>ln</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>15</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000596732940000042.GIF" wi="1243" he="169" /></maths>式中:n为切断参数;c为绝缘子的标准方差与其50%的耐受电压的比值;A为与绝缘子形状和污秽度有关的系数;L为绝缘子长度,m;σ为等值盐密值影响污闪电压的特征指数;ε和h为中间计算量;(4)、根据绝缘子固有因素指标m和绝缘子表面污秽程度指标P<sub>f</sub>的乘积判断绝缘子污闪综合预警指标R,确定预警等级:具体包括以下步骤:4.1)、计算绝缘子污闪综合预警指标R:R=m×P<sub>f</sub>   (16)式中:R为绝缘子污闪综合预警指标;m为绝缘子固有因素指标;P<sub>f</sub>为绝缘子表面污秽程度指标;4.2)、确定预警等级:根据绝缘子污闪综合预警指标R确定预警等级,预警等级共分黄色、橙色、红色三个等级,具体如表1所示:表1  预警等级表<img file="FDA0000596732940000043.GIF" wi="2081" he="339" />对存在危险的绝缘子,根据其预警等级的大小,发布相应的预警信息:预警信息内容包括:预警信息编号、预警发布时间、预警等级、线路名称、杆塔ID。
地址 100031 北京市西城区西长安街86号