发明名称 一种基于时空因素影响的输电线路覆冰风险评估方法
摘要 本发明公开了一种基于时空因素影响的输电线路覆冰风险评估方法,其包括以下步骤:步骤1、综合计算输电线路各部件的部件风险值,获取整条输电线路的线路风险值;步骤2、综合考虑整条输电线路的运行环境系数,得到整条输电线路的整体风险值。本发明基于输电线路部件风险的基础上,综合时空因素,构建多因素输电线路风险分析方法拟通过灰色模糊综合评判得到贴近实际的风险评判结果。
申请公布号 CN105930964A 申请公布日期 2016.09.07
申请号 CN201610239202.0 申请日期 2016.04.15
申请人 中国南方电网有限责任公司超高压输电公司检修试验中心 发明人 宋云海;李晋伟;甘振宁;陈岳;王奇;陈浩宁
分类号 G06Q10/06(2012.01)I;G06Q50/06(2012.01)I 主分类号 G06Q10/06(2012.01)I
代理机构 广州科粤专利商标代理有限公司 44001 代理人 孔德超;黄培智
主权项 一种基于时空因素影响的输电线路覆冰风险评估方法,其特征在于,其包括以下步骤:步骤1、综合计算输电线路各部件的部件风险值,获取整条输电线路的线路风险值<img file="FDA0000966350670000011.GIF" wi="65" he="63" />所述整条输电线路由m段输电线路组成,每段输电线路又由n个部件组成;所述步骤1包括以下步骤:步骤11、计算输电线路每个部件的部件风险值R′:其包括以下步骤:步骤111、计算与每个部件的风险评估相关的第i个状态量的风险度Q<sub>i</sub>;其中,1≤i≤p,p为与每个部件的风险评估相关的状态量总数,所述风险度Q<sub>i</sub>的计算为下列情形之一:第一种情形:当第i个状态量为直接观测状态量时:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>Q</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><munderover><mo>&Sigma;</mo><mrow><mi>l</mi><mo>=</mo><mn>1</mn></mrow><mi>A</mi></munderover><mrow><mo>(</mo><msub><mi>q</mi><mi>l</mi></msub><mo>&times;</mo><msup><mn>1.05</mn><mrow><msub><mi>x</mi><mi>l</mi></msub><mo>-</mo><mn>10</mn></mrow></msup><mo>)</mo></mrow></mrow><mi>A</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000966350670000012.GIF" wi="1166" he="199" /></maths>其中,q<sub>l</sub>为第l个打分人员给出的风险度量化值,x<sub>l</sub>为第l个打分人员的工作年限,A为打分人员的总人数,1≤l≤A;第二种情形:当第i个状态量为正状态参量时:所述正状态参量满足:μ<sub>1</sub>&lt;ξ<sub>1</sub>,其中,μ<sub>1</sub>为正状态参量的正常限值,ξ<sub>1</sub>为正状态参量的故障限值;<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>Q</mi><mi>i</mi></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mrow><msub><mi>x</mi><mn>1</mn></msub><mo>&le;</mo><msub><mi>&mu;</mi><mn>1</mn></msub></mrow></mtd></mtr><mtr><mtd><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>x</mi><mn>1</mn></msub><mo>-</mo><msub><mi>&mu;</mi><mn>1</mn></msub></mrow><mrow><msub><mi>&xi;</mi><mn>1</mn></msub><mo>-</mo><msub><mi>&mu;</mi><mn>1</mn></msub></mrow></mfrac><mo>)</mo></mrow><mrow><mi>k</mi><mn>1</mn></mrow></msup></mtd><mtd><mrow><msub><mi>&mu;</mi><mn>1</mn></msub><mo>&lt;</mo><msub><mi>x</mi><mn>1</mn></msub><mo>&lt;</mo><msub><mi>&xi;</mi><mn>1</mn></msub></mrow></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mrow><msub><mi>x</mi><mn>1</mn></msub><mo>&GreaterEqual;</mo><msub><mi>&xi;</mi><mn>1</mn></msub></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000966350670000013.GIF" wi="1213" he="303" /></maths>其中,x<sub>1</sub>为正状态参量的监测值;k1为正状态参量的风险趋势指数;第三种情形:当第i个状态量为负状态参量时:所述负状态参量满足:μ<sub>2</sub>&gt;ξ<sub>2</sub>,其中,μ<sub>2</sub>为负状态参量的正常限值,ξ<sub>2</sub>为负状态参量的故障限值;<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>Q</mi><mi>i</mi></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mrow><msub><mi>x</mi><mn>2</mn></msub><mo>&le;</mo><msub><mi>&xi;</mi><mn>2</mn></msub></mrow></mtd></mtr><mtr><mtd><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>&mu;</mi><mn>2</mn></msub><mo>-</mo><msub><mi>x</mi><mn>2</mn></msub></mrow><mrow><msub><mi>&mu;</mi><mn>2</mn></msub><mo>-</mo><msub><mi>&xi;</mi><mn>2</mn></msub></mrow></mfrac><mo>)</mo></mrow><mrow><mi>k</mi><mn>2</mn></mrow></msup></mtd><mtd><mrow><msub><mi>&xi;</mi><mn>2</mn></msub><mo>&lt;</mo><msub><mi>x</mi><mn>2</mn></msub><mo>&lt;</mo><msub><mi>&mu;</mi><mn>2</mn></msub></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><msub><mi>x</mi><mn>2</mn></msub><mo>&GreaterEqual;</mo><msub><mi>&mu;</mi><mn>2</mn></msub></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000966350670000014.GIF" wi="1246" he="300" /></maths>其中,x<sub>2</sub>为负状态参量的监测值;k2为负状态参量的风险趋势指数;第四种情形:当第i个状态量为偏差性状态参量时:所述偏差性状态参量满足:ξ<sub>3</sub>&lt;μ<sub>3</sub>&lt;μ<sub>4</sub>&lt;ξ<sub>4</sub>,其中,μ<sub>3</sub>和μ<sub>4</sub>分别为偏差性状态参量的正常下限值和正常上限值,ξ<sub>3</sub>和ξ<sub>4</sub>分别为偏差性状态参量的故障下限值和故障上限值;<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>Q</mi><mi>i</mi></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mrow><msub><mi>x</mi><mn>3</mn></msub><mo>&le;</mo><msub><mi>&xi;</mi><mn>3</mn></msub></mrow></mtd></mtr><mtr><mtd><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>&mu;</mi><mn>3</mn></msub><mo>-</mo><msub><mi>x</mi><mn>3</mn></msub></mrow><mrow><msub><mi>&mu;</mi><mn>3</mn></msub><mo>-</mo><msub><mi>&xi;</mi><mn>3</mn></msub></mrow></mfrac><mo>)</mo></mrow><mrow><mi>k</mi><mn>3</mn></mrow></msup></mtd><mtd><mrow><msub><mi>&xi;</mi><mn>3</mn></msub><mo>&lt;</mo><msub><mi>x</mi><mn>3</mn></msub><mo>&lt;</mo><msub><mi>&mu;</mi><mn>3</mn></msub></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><msub><mi>&mu;</mi><mn>3</mn></msub><mo>&le;</mo><msub><mi>x</mi><mn>3</mn></msub><mo>&le;</mo><msub><mi>&mu;</mi><mn>4</mn></msub></mrow></mtd></mtr><mtr><mtd><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>x</mi><mn>3</mn></msub><mo>-</mo><msub><mi>&mu;</mi><mn>4</mn></msub></mrow><mrow><msub><mi>&xi;</mi><mn>4</mn></msub><mo>-</mo><msub><mi>&mu;</mi><mn>4</mn></msub></mrow></mfrac><mo>)</mo></mrow><mrow><mi>k</mi><mn>3</mn></mrow></msup></mtd><mtd><mrow><msub><mi>&mu;</mi><mn>4</mn></msub><mo>&lt;</mo><msub><mi>x</mi><mn>3</mn></msub><mo>&lt;</mo><msub><mi>&xi;</mi><mn>4</mn></msub></mrow></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mrow><msub><mi>x</mi><mn>3</mn></msub><mo>&GreaterEqual;</mo><msub><mi>&xi;</mi><mn>4</mn></msub></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000966350670000021.GIF" wi="1237" he="517" /></maths>其中,x<sub>3</sub>为偏差性状态参量的监测值;k3为偏差性状态参量的风险趋势指数;步骤112、计算每个部件的第i个状态量的状态量隶属度α<sub>i</sub><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>n</mi><mi>&beta;</mi></msub><msub><mi>&beta;</mi><mi>i</mi></msub><mo>+</mo><msub><mi>n</mi><mi>&gamma;</mi></msub><msub><mi>&gamma;</mi><mi>i</mi></msub></mrow><mrow><msub><mi>n</mi><mi>&beta;</mi></msub><mo>+</mo><msub><mi>n</mi><mi>&gamma;</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000966350670000022.GIF" wi="1254" he="151" /></maths>其中,β<sub>i</sub>和γ<sub>i</sub>分别第i个状态量的历史数据隶属度和权重隶属度,n<sub>β</sub>和n<sub>γ</sub>分别是针对历史数据隶属度和权重隶属度的权重系数;步骤113、计算每个部件的第i个状态量的状态量风险值r<sub>i</sub>:r<sub>i</sub>=Q<sub>i</sub>×α<sub>i</sub>        (6)步骤114、计算每个部件的部件风险值R′:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msup><mi>R</mi><mo>&prime;</mo></msup><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></munderover><msub><mi>r</mi><mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000966350670000023.GIF" wi="1126" he="125" /></maths>步骤12、对所述步骤11计算的每个部件的部件风险值R′进行修正,获得修正后每个部件的部件风险值R:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mi>R</mi><mo>=</mo><mfrac><msup><mi>R</mi><mo>&prime;</mo></msup><mn>1.5</mn></mfrac><mo>&times;</mo><msub><mi>FY</mi><mi>t</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000966350670000024.GIF" wi="1165" he="120" /></maths>其中:FY<sub>t</sub>为所述每个部件的年限系数,FY<sub>t</sub>=0.05AG<sub>t</sub>+1        (9)其中,AG<sub>t</sub>为被求取年份的老化指数,<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>AG</mi><mi>t</mi></msub><mo>=</mo><msub><mi>AG</mi><mn>0</mn></msub><mo>&times;</mo><mi>exp</mi><mo>{</mo><mfrac><mrow><mi>B</mi><mrow><mo>(</mo><mi>T</mi><mo>-</mo><msub><mi>T</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow><mrow><msub><mi>age</mi><mi>f</mi></msub></mrow></mfrac><mo>}</mo><mo>&times;</mo><msub><mi>f</mi><mi>mod</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000966350670000031.GIF" wi="1325" he="166" /></maths>其中,AG<sub>0</sub>为每个部件的初始老化指数,B为老化常数;T为被求取年份的老化指数AG<sub>t</sub>对应的年份;T<sub>0</sub>为设备投运年份;age<sub>f</sub>为老化速率修正系数;f<sub>mod</sub>为修正系数;步骤13、计算每段输电线路的段风险值<img file="FDA0000966350670000032.GIF" wi="70" he="68" /><maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><mover><mi>R</mi><mo>&OverBar;</mo></mover><mo>=</mo><mn>1</mn><mo>-</mo><munderover><mo>&Pi;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>R</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000966350670000033.GIF" wi="1022" he="127" /></maths>其中,R<sub>j</sub>为每段输电线路中第j个部件的部件风险值,1≤j≤n;步骤14、计算整条输电线路的线路风险值<img file="FDA0000966350670000034.GIF" wi="70" he="63" /><maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><mover><mi>R</mi><mo>~</mo></mover><mo>=</mo><mn>1</mn><mo>-</mo><munderover><mo>&Pi;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mover><mi>R</mi><mo>&OverBar;</mo></mover><mi>k</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000966350670000035.GIF" wi="1078" he="126" /></maths>其中,<img file="FDA0000966350670000036.GIF" wi="62" he="61" />为整条输电线路中第k段输电线路的段风险值,1≤k≤m;步骤2、综合考虑整条输电线路的运行环境系数tp,得到整条输电线路的整体风险值<img file="FDA0000966350670000037.GIF" wi="75" he="79" /><maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><msup><mover><mi>R</mi><mo>&OverBar;</mo></mover><mo>&prime;</mo></msup><mo>=</mo><mover><mi>R</mi><mo>~</mo></mover><mo>&CenterDot;</mo><mi>t</mi><mi>p</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000966350670000038.GIF" wi="945" he="87" /></maths>其中,运行环境系数tp的表达式为:<img file="FDA0000966350670000039.GIF" wi="1293" he="151" />其中,t和p分别表示线路覆冰运行环境的时间系数和地理位置系数。
地址 510663 广东省广州市黄埔区科学大道223号2号楼检修试验中心