发明名称 一种输电线路雷电绕击跳闸在线预警方法
摘要 本发明公开了一种输电线路雷电绕击跳闸在线预警方法,该方法通过统计历史雷电绕击跳闸信息,应用二维信息扩散理论和条件概率方法得到输电线路绕击跳闸概率分布。选取回波强度、回波顶高,垂直积累液态水含量等雷达预报数据,建立了基于反向传播神经网络的雷电流幅值预测模型,根据预测的雷电流幅值和落雷侧面距离,应用输电线路跳闸概率分布模型,对输电线路绕击跳闸概率进行实时预警并发布预警等级。本发明根据气象雷达实时预报数据,得出输电线路的预报跳闸概率,发布预警信号,给电网调度的运行人员决策分析提供依据,及时做好输电线路调度策略,提高供电可靠性,降低电网的经济损失,提高电网的可靠运行能力。
申请公布号 CN103207340B 申请公布日期 2015.04.08
申请号 CN201310159521.7 申请日期 2013.05.02
申请人 深圳供电局有限公司;重庆大学 发明人 程韧俐;熊小伏;徐旭辉;包博;沈智健;陈洪云;于洋
分类号 G01R31/00(2006.01)I 主分类号 G01R31/00(2006.01)I
代理机构 重庆博凯知识产权代理有限公司 50212 代理人 李海华
主权项 一种输电线路雷电绕击跳闸在线预警方法,其特征在于:按如下步骤进行:1)根据历史雷电绕击跳闸事件,应用信息扩散方法分析得到绕击条件概率,然后结合条件概率方法,得到特定雷电流幅值和落雷侧面距离下输电线路发生雷电绕击跳闸的概率分布;2)根据历史雷达预报数据中的回波强度、回波顶高和垂直积累液态水含量三个参数信息以及对应的雷电流幅值信息,应用bp神经网络训练得到雷达预报数据中回波强度、回波顶高和垂直积累液态水含量三个参数信息跟雷电流幅值之间的关系;3)对输电线路的经过区域进行网格划分,根据雷电设备得到包括回波强度、回波顶高和垂直积累液态水含量在内的每个网格的实际雷达预报数据,根据第2)步所得结果得出每个网格预测的雷电流幅值;同时依据网格信息计算得到每个网格输电线路的落雷侧面距离;4)根据第3)步预测得到的每个网格雷电流幅值和计算得到的落雷侧面距离,结合第1)步的输电线路雷电绕击跳闸概率分布,得出每个网格在实际雷达预报数据下输电线路发生绕击跳闸的概率,并发布预警信息;其中,第1)步中特定雷电流幅值和落雷侧面距离下输电线路发生雷电绕击跳闸的概率分布按如下步骤获得:1.1)获取历史雷电绕击跳闸事件发生时的信息作为概率分析的样本值,包括雷电流幅值和落雷侧面距离:(X,Y)={(x<sub>1</sub>,y<sub>1</sub>),(x<sub>2</sub>,y<sub>2</sub>),…,(x<sub>M</sub>,y<sub>M</sub>)}式中:X表示雷电流幅值;      Y表示落雷侧面距离;1.2)信息扩散方法分析样本数据选取历史M次故障事件作为样本值,分别进行电流幅值指标和侧面距离指标扩展;其指标论域分别为:U={u<sub>1</sub>,u<sub>2</sub>,…,u<sub>m</sub>}V={v<sub>1</sub>,v<sub>2</sub>,…,v<sub>n</sub>}式中;U表示电流幅值指标论域;      V表示侧面距离指标论域;把观测样本(x<sub>k</sub>,y<sub>k</sub>)依下式将其信息扩散到u<sub>i</sub>,v<sub>j</sub>中的所有点;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>f</mi><mi>k</mi></msub><mrow><mo>(</mo><msub><mi>u</mi><mi>i</mi></msub><mo>,</mo><msub><mi>v</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>&pi;</mi><msub><mi>h</mi><mn>1</mn></msub><msub><mi>h</mi><mn>2</mn></msub></mrow></mfrac><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mrow><mo>(</mo><mfrac><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>k</mi></msub><mo>-</mo><msub><mi>u</mi><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><msup><msub><mi>h</mi><mn>1</mn></msub><mn>2</mn></msup></mfrac><mo>+</mo><mfrac><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>k</mi></msub><mo>-</mo><msub><mi>v</mi><mi>j</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><msup><msub><mi>h</mi><mn>2</mn></msub><mn>2</mn></msup></mfrac><mo>)</mo></mrow><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000651848730000021.GIF" wi="860" he="149" /></maths>式中:h<sub>1</sub>,h<sub>2</sub>为扩散系数,可以根据样本集合中样本的最大值b和最小值a及样本个数M来确定,其计算公式为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>h</mi><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1.4230</mn><mo>&times;</mo><mfrac><mrow><mi>b</mi><mo>-</mo><mi>a</mi></mrow><mrow><mi>M</mi><mo>-</mo><mn>1</mn></mrow></mfrac></mtd><mtd><mi>M</mi><mo>&le;</mo><mn>9</mn></mtd></mtr><mtr><mtd><mn>1.4208</mn><mo>&times;</mo><mfrac><mrow><mi>b</mi><mo>-</mo><mi>a</mi></mrow><mrow><mi>M</mi><mo>-</mo><mn>1</mn></mrow></mfrac></mtd><mtd><mi>M</mi><mo>&GreaterEqual;</mo><mn>10</mn></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000651848730000022.GIF" wi="500" he="249" /></maths>令<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>C</mi><mi>k</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>f</mi><mi>k</mi></msub><mrow><mo>(</mo><msub><mi>u</mi><mi>i</mi></msub><mo>,</mo><msub><mi>v</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000651848730000023.GIF" wi="358" he="118" /></maths>相应的模糊子集的隶属函数为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>&mu;</mi><mi>k</mi></msub><mrow><mo>(</mo><msub><mi>u</mi><mi>i</mi></msub><mo>,</mo><msub><mi>v</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>f</mi><mi>k</mi></msub><mrow><mo>(</mo><msub><mi>u</mi><mi>i</mi></msub><mo>,</mo><msub><mi>v</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow><msub><mi>C</mi><mi>k</mi></msub></mfrac></mrow>]]></math><img file="FDA0000651848730000024.GIF" wi="389" he="140" /></maths>其中μ<sub>k</sub>(u<sub>i</sub>,v<sub>j</sub>)为样本的归一化信息分布;对μ<sub>k</sub>(u<sub>i</sub>,v<sub>j</sub>)进行处理,可得到一种效果较好的风险评估结果;令<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>q</mi><mrow><mo>(</mo><msub><mi>u</mi><mi>i</mi></msub><mo>,</mo><msub><mi>v</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>t</mi></munderover><msub><mi>&mu;</mi><mi>k</mi></msub><mrow><mo>(</mo><msub><mi>u</mi><mi>i</mi></msub><mo>,</mo><msub><mi>v</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000651848730000025.GIF" wi="406" he="125" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>Q</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mi>q</mi><mrow><mo>(</mo><msub><mi>u</mi><mi>i</mi></msub><mo>,</mo><msub><mi>v</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000651848730000026.GIF" wi="326" he="133" /></maths>样本落在(u<sub>i</sub>,v<sub>j</sub>)处的频率值为:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mi>p</mi><mrow><mo>(</mo><msub><mi>u</mi><mi>i</mi></msub><mo>,</mo><msub><mi>v</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>q</mi><mrow><mo>(</mo><msub><mi>u</mi><mi>i</mi></msub><mo>,</mo><msub><mi>v</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow><mi>Q</mi></mfrac></mrow>]]></math><img file="FDA0000651848730000027.GIF" wi="349" he="134" /></maths>1.3)得到输电线路雷电绕击跳闸概率分布p(u<sub>i</sub>,v<sub>j</sub>)的物理意义为输电线路发生跳闸时,雷电流幅值为u<sub>i</sub>,侧面距离为v<sub>j</sub>的概率;其本质是条件概率,可记为p(BC|A),A表示跳闸事件,B表示雷电流幅值,C表示侧面距离;根据条件概率公式:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mi>p</mi><mrow><mo>(</mo><mi>BC</mi><mo>|</mo><mi>A</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>p</mi><mrow><mo>(</mo><mi>ABC</mi><mo>)</mo></mrow></mrow><mrow><mi>p</mi><mrow><mo>(</mo><mi>A</mi><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000651848730000028.GIF" wi="361" he="124" /></maths>可得:<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><mi>p</mi><mrow><mo>(</mo><mi>A</mi><mo>|</mo><mi>BC</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>p</mi><mrow><mo>(</mo><mi>ABC</mi><mo>)</mo></mrow></mrow><mrow><mi>p</mi><mrow><mo>(</mo><mi>BC</mi><mo>)</mo></mrow></mrow></mfrac><mo>=</mo><mfrac><mrow><mi>p</mi><mrow><mo>(</mo><mi>BC</mi><mo>|</mo><mi>A</mi><mo>)</mo></mrow><mo>&times;</mo><mi>p</mi><mrow><mo>(</mo><mi>A</mi><mo>)</mo></mrow></mrow><mrow><mi>p</mi><mrow><mo>(</mo><mi>B</mi><mo>)</mo></mrow><mi>p</mi><mrow><mo>(</mo><mi>C</mi><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000651848730000029.GIF" wi="692" he="132" /></maths>式中:p(A|BC)表示在特定雷电流幅值、侧面距离条件下输电线路发生绕击跳闸的概率;p(A)表示输电线路绕击跳闸的概率,p(B)表示雷电流幅值概率,p(C)表示侧面距离分布概率,p(B)和p(C)为相互独立事件;在输电线路绕击概率模型中,绕击跳闸概率可用输电线路绕击率计算得到;计算公式如下所示:p(A)=1‑exp(‑λt)式中:λ为输电线路绕击跳闸率;t为M次故障发生的时间跨度;雷电流幅值概率的计算公式如下:p(B)=P(i&gt;I)‑P(i&gt;I‑I<sub>m</sub>)<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><mi>P</mi><mrow><mo>(</mo><mo>></mo><mi>I</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><msup><mrow><mo>(</mo><mi>I</mi><mo>/</mo><mn>31</mn><mo>)</mo></mrow><mn>2.6</mn></msup></mrow></mfrac></mrow>]]></math><img file="FDA0000651848730000031.GIF" wi="388" he="121" /></maths><maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><msub><mi>I</mi><mi>m</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>u</mi><mi>m</mi></msub><mo>-</mo><msub><mi>u</mi><mn>1</mn></msub></mrow><mi>m</mi></mfrac></mrow>]]></math><img file="FDA0000651848730000032.GIF" wi="212" he="109" /></maths>式中:u<sub>m</sub>,u<sub>1</sub>分别为电流幅值指标论域的最大值,最小值;m为电流幅值指标个数;侧面距离分布概率采用下式计算:<maths num="0012" id="cmaths0012"><math><![CDATA[<mrow><mi>p</mi><mrow><mo>(</mo><mi>C</mi><mo>)</mo></mrow><mo>=</mo><mfrac><msub><mi>N</mi><mi>L</mi></msub><mi>N</mi></mfrac><mfrac><mrow><mrow><mo>(</mo><msub><mi>v</mi><mi>n</mi></msub><mo>-</mo><msub><mi>v</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>/</mo><mi>n</mi></mrow><mrow><mn>0.5</mn><mi>L</mi></mrow></mfrac></mrow>]]></math><img file="FDA0000651848730000033.GIF" wi="424" he="125" /></maths>式中:N<sub>L</sub>表示线路经过网格数,N表示总网格数,L表示网格宽度;v<sub>n</sub>,v<sub>1</sub>分别为侧面距离指标论域的最大值,最小值;n为侧面距离指标个数。
地址 518001 广东省深圳市罗湖区深南东路4020号