发明名称 一种考虑工频电压的高压输电线路雷击跳闸率分析方法
摘要 一种考虑工频电压的高压输电线路雷击跳闸率分析方法:S1构建概率函数;S2建立联合密度函数;S3考虑工频电压的反击跳闸率;S4考虑工频电压的绕击跳闸率。本发明主要是针对工程上常忽略工频电压值所计算的雷击跳闸率过于保守的问题,提出考虑工频电压下求解雷击跳闸率的新方法,该方法是结合雷电概率分布的特性,并利用概率学中的二维联合密度函数的方法,提出了考虑工频电压下雷击跳闸率的计算公式,利用该公式可以更精确的计算出雷击跳闸率,以便更加准确地评估输电线路的安全性能,且本发明实现方法简便,具有很好的实用性。
申请公布号 CN102854414A 申请公布日期 2013.01.02
申请号 CN201210279298.5 申请日期 2012.08.06
申请人 广东电网公司东莞供电局;华南理工大学 发明人 陈锡阳;许彬;汪逍旻;王伟然;刘刚;杨挺;孔慧超;林建华;尹创荣;王凯;张弦
分类号 G01R31/00(2006.01)I 主分类号 G01R31/00(2006.01)I
代理机构 广州知友专利商标代理有限公司 44104 代理人 周克佑
主权项 1.一种考虑工频电压的高压输电线路雷击跳闸率分析方法,包括以下步骤:S1构建概率函数将雷电流分布函数P<sub>I</sub>求导,得出雷电流密度函数p<sub>I</sub>,然后根据耐雷水平计算公式得出不同大小雷电流幅值对应的绝缘子串两端产生的雷击过电压U<sub>g</sub>,推导出绝缘子串两端雷击过电压的密度函数如下式所示:设I<sub>m</sub>=U<sub>50%</sub>/R,则有:U<sub>g</sub>=I/R    (1)计算雷击过电压时仅考虑绝缘子串闪络之前的瞬间状态;按规程中给出的雷电流密度函数计算,将I=U<sub>g</sub>/R带入公式<img file="FDA00001980335400011.GIF" wi="270" he="115" />得出<maths num="0001"><![CDATA[<math><mrow><msub><mi>P</mi><mi>g</mi></msub><mo>=</mo><msup><mn>10</mn><mrow><mo>-</mo><mfrac><mi>Ug</mi><mrow><mn>88</mn><mi>R</mi></mrow></mfrac></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>设m=1/(88R),则有p<sub>g</sub>=(mln10)0.1<sup>mUg</sup>    (3)其中,I<sub>m</sub>为耐雷水平幅值;I为任意幅值雷电流;P<sub>g</sub>为雷击过电压不大于U<sub>g</sub>的概率分布函数;p<sub>g</sub>为绝缘子串雷击后产生不同电压值的概率分布函数;同理根据工频电压公式U<sub>f</sub>=U<sub>m</sub>sinωt推导出工频电压的概率密度函数p<sub>f</sub>和工频电压值的概率分布函数P<sub>f</sub>;S2建立联合密度函数因为雷电流在绝缘子串两端产生电压值与工频电压值之间互为独立,所以根据概率学可知:p<sub>U</sub>=p<sub>g</sub>p<sub>f</sub>                  (4)P<sub>U</sub>=∫∫p<sub>g</sub>p<sub>f</sub>dU<sub>g</sub>dU<sub>f</sub>    (5)已知U<sub>g</sub>=U-U<sub>f</sub>而此处需计算U<sub>g</sub>+U<sub>f</sub>&gt;U<sub>50%</sub>概率,则<maths num="0002"><![CDATA[<math><mrow><msub><mi>P</mi><mi>U</mi></msub><mrow><mo>(</mo><mi>U</mi><mo>></mo><msub><mi>U</mi><mrow><mn>50</mn><mo>%</mo></mrow></msub><mo>)</mo></mrow><mo>=</mo><msubsup><mo>&Integral;</mo><mrow><mi>U</mi><mn>50</mn><mo>%</mo><mo>+</mo><mi>Um</mi></mrow><mrow><mi>U</mi><mn>50</mn><mo>%</mo><mo>+</mo><mi>Um</mi></mrow></msubsup><msub><mi>p</mi><mi>g</mi></msub><mi>d</mi><msub><mi>U</mi><mi>g</mi></msub><msubsup><mo>&Integral;</mo><mrow><mo>-</mo><mi>Um</mi></mrow><mi>Um</mi></msubsup><msub><mi>p</mi><mi>f</mi></msub><mi>d</mi><msub><mi>U</mi><mi>f</mi></msub><mo>+</mo><msub><mi>P</mi><mi>g</mi></msub><mrow><mo>(</mo><msub><mi>U</mi><mi>g</mi></msub><mo>></mo><msub><mi>U</mi><mrow><mn>50</mn><mo>%</mo></mrow></msub><mo>+</mo><msub><mi>U</mi><mi>m</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,P<sub>U</sub>(U&gt;U<sub>50%</sub>)为叠加上工频电压后绝缘子串两端电压超过U<sub>50%</sub>的概率;S3考虑工频电压的反击跳闸率先求三相叠加后的概率分布函数,如式(7):P<sub>f3</sub>=|2(U<sub>f</sub>/U<sub>m</sub>-0.5)|                (7)其中,P<sub>f3</sub>为工频电压不大于U<sub>f</sub>的概率分布,而U<sub>f</sub>取值为[U<sub>m</sub>/2,U<sub>m</sub>];分析输电线路耐雷水平时,采用集中参数模型代替多波阻抗模型,其耐雷水平如公式(8)所示:I<sub>m</sub>(1-K<sub>c</sub>)[β<sub>g</sub>(R<sub>ch</sub>+L<sub>gt</sub>/τ)+h<sub>d</sub>/2.6]=U<sub>50%</sub>(8)式中,β<sub>g</sub>为杆塔分流系数;K<sub>c</sub>为导线与避雷线耦合系数;R<sub>ch</sub>为接地电阻;L<sub>gt</sub>为杆塔等效电抗;τ为雷电流波头时间;h<sub>d</sub>为导线对地高度;根据直击耐雷水平计算公式,知R=(1-K<sub>c</sub>)[β<sub>g</sub>(R<sub>ch</sub>+L<sub>gt</sub>/τ)+h<sub>d</sub>/2.6](9)则,雷击过电压公式为U<sub>g</sub>=IR                (10)将雷击过电压概率密度函数和三相工频电压概率分布函数联合,最后可以得出雷击塔顶时雷击跳闸率的计算公式如式(11)所示:<maths num="0003"><![CDATA[<math><mrow><msub><mi>P</mi><mi>n</mi></msub><mo>=</mo><msubsup><mo>&Integral;</mo><mrow><mi>U</mi><mn>0</mn></mrow><mrow><mi>U</mi><mn>1</mn></mrow></msubsup><msub><mi>p</mi><mi>g</mi></msub><mo>|</mo><mn>2</mn><mrow><mo>(</mo><mrow><mo>(</mo><msub><mi>U</mi><mi>g</mi></msub><mo>-</mo><msub><mi>U</mi><mrow><mn>50</mn><mo>%</mo></mrow></msub><mo>)</mo></mrow><mo>/</mo><msub><mi>U</mi><mi>m</mi></msub><mo>-</mo><mn>0.5</mn><mo>)</mo></mrow><mo>|</mo><mi>d</mi><msub><mi>U</mi><mi>g</mi></msub><mo>+</mo><msubsup><mo>&Integral;</mo><mrow><mi>U</mi><mn>1</mn></mrow><mrow><mo>+</mo><mo>&infin;</mo></mrow></msubsup><msub><mi>p</mi><mi>g</mi></msub><mi>d</mi><msub><mi>U</mi><mi>g</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,P<sub>n</sub>为考虑工频电压影响后雷击塔顶的雷击跳闸率;p<sub>g</sub>为雷击过电压概率密度函数;U<sub>0</sub>=U<sub>50%</sub>-U<sub>m</sub>、U<sub>1</sub>=U<sub>50%</sub>-U<sub>m</sub>/2;S4考虑工频电压的绕击跳闸率包括以下子步骤;S4-1根据EGM推算出输电线路的最小绕击跳闸电流I<sub>min</sub>和最大绕击跳闸电流I<sub>max</sub>;S4-2根据EGM计算出不同雷电流I对应的绕击率P<sub>α</sub>,将绕击率带进雷电流概率密度函数中得到绕击雷电流概率密度函数p<sub>Ir</sub>=P<sub>I</sub>P<sub>α</sub>;S4-3根据绕击耐雷水平计算公式I<sub>r</sub>=U<sub>50%</sub>/Z,可以得到雷电流转换成绝缘子串两端电压的系数R=Z,利用转换系数将绕击雷电流密度函数转换成雷击过电压值密度函数p<sub>gr</sub>,Z为导线并联波阻抗值;S4-4雷电流绕击导线时仅需要考虑被绕击导线的工频电压对于跳闸率的影响,故工频电压的概率密度函数可利用单项概率密度函数分析,如式(12)(13)所示:<maths num="0004"><![CDATA[<math><mrow><msub><mi>p</mi><mrow><mi>f</mi><mn>1</mn></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>U</mi><mi>f</mi></msub><mo>/</mo><msub><mi>U</mi><mi>m</mi></msub></mrow><mn>2</mn></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msub><mi>P</mi><mrow><mi>f</mi><mn>1</mn></mrow></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mfrac><mrow><mo>-</mo><msub><mi>U</mi><mi>f</mi></msub><mo>/</mo><msub><mi>U</mi><mi>m</mi></msub></mrow><mn>2</mn></mfrac><mo>,</mo><mo>-</mo><msub><mi>U</mi><mi>m</mi></msub><mo>&lt;</mo><mi>U</mi><mo>&lt;</mo><mn>0</mn></mtd></mtr><mtr><mtd><mfrac><mrow><mn>1</mn><mo>+</mo><msub><mi>U</mi><mi>f</mi></msub><mo>/</mo><msub><mi>U</mi><mi>m</mi></msub></mrow><mn>2</mn></mfrac><mo>,</mo><mn>0</mn><mo>&lt;</mo><mi>U</mi><mo>&lt;</mo><msub><mi>U</mi><mi>m</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,P<sub>f1</sub>为单项工频电压值不大于U的概率分布函数;p<sub>f1</sub>为单项工频电压值的概率密度函数;将工频电压的概率密度函数和绕击雷击过电压概率密度函数联合,可以得出绕击时考虑工频电压影响下的雷击跳闸率的计算式(14):<maths num="0006"><![CDATA[<math><mrow><msub><mi>P</mi><mi>r</mi></msub><mo>=</mo><msubsup><mo>&Integral;</mo><mrow><mi>U</mi><mn>0</mn></mrow><mrow><mi>U</mi><mn>1</mn></mrow></msubsup><msub><mi>p</mi><mi>gr</mi></msub><msub><mi>P</mi><mi>fr</mi></msub><mo>|</mo><mfrac><mrow><mo>(</mo><msub><mi>U</mi><mi>gr</mi></msub><mo>-</mo><msub><mi>U</mi><mrow><mn>50</mn><mo>%</mo></mrow></msub><mo>)</mo></mrow><msub><mi>U</mi><mi>m</mi></msub></mfrac><mo>|</mo><mi>d</mi><msub><mi>U</mi><mi>gr</mi></msub><mo>+</mo><msubsup><mo>&Integral;</mo><mrow><mi>U</mi><mn>1</mn></mrow><mrow><mo>+</mo><mo>&infin;</mo></mrow></msubsup><msub><mi>p</mi><mi>gr</mi></msub><msub><mi>P</mi><mi>fr</mi></msub><mi>d</mi><msub><mi>U</mi><mi>gr</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>14</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,P<sub>r</sub>为考虑工频电压影响后雷电流绕击导线的雷击跳闸率;U<sub>0</sub>=U<sub>50%</sub>-U<sub>m</sub>;U<sub>1</sub>=U<sub>50%</sub>+U<sub>m</sub>。
地址 523008 广东省东莞市东城区东城楼239号