发明名称 基于多元物理量数学模型的输电线路覆冰预测方法
摘要 一种基于多元物理量数学模型的输电线路覆冰预测方法,包括如下步骤:(1)获取当前时刻的温湿度及风力气象参数;(2)分别判断当前各气象数据是否满足覆冰条件,温度-8~0℃,湿度90%RH以上,风速2~7m/s,若满足条件则转到(3),否则进入下一次采样周期,即回到(1);(3)通过张力传感器获得的数据值与初始值做比较,若<img file="DDA00002111904900011.GIF" wi="377" he="92" />说明导线有覆冰,转到步骤(4),否则返回到(1);(4)通过地面终端获取当前环境气压值,计算覆冰厚度b;(5)通过覆冰厚度,温湿度等参数计算任意时刻t<sub>i</sub>导线的综合比载;(6)通过综合比载得到档距端点处的轴向应力σ<sub>x</sub>,令σ<sub>x</sub>=σ<sub>max</sub>,计算断裂剩余时间t,通过终端发出预警。本发明实现有效监测、准确性良好。
申请公布号 CN102927949B 申请公布日期 2015.08.05
申请号 CN201210330175.X 申请日期 2012.09.07
申请人 浙江工业大学 发明人 陈庆章;李兴华;王凯;倪云峰;陈宇铮;范聪玲;吴荣杰;杨帆;宦若红;赵小敏;毛科技;雷艳静;任佳萍;丁舒羽
分类号 G01B21/08(2006.01)I 主分类号 G01B21/08(2006.01)I
代理机构 杭州天正专利事务所有限公司 33201 代理人 王兵;王利强
主权项 一种基于多元物理量数学模型的输电线路覆冰预测方法,其特征在于:所述预测方法包括如下步骤:(1)采用温湿度传感器、风力传感器获取当前时刻的温湿度及风力气象参数;(2)分别判断当前各气象数据是否满足覆冰条件,即:温度‑8~0℃,湿度90%RH以上,风速2~7m/s,若满足条件则转到(3),否则进入下一次采样周期,即返回到(1);(3)通过张力传感器获得的数据值与初始值做比较,若<img file="FDA0000701785000000015.GIF" wi="421" he="126" />说明导线有覆冰,转到步骤(4),否则返回到(1);(4)通过地面终端获取当前环境气压值,计算覆冰厚度r,计算公式如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>r</mi><mo>=</mo><msubsup><mo>&Integral;</mo><mn>0</mn><msub><mi>t</mi><mi>i</mi></msub></msubsup><mi>f</mi><mrow><mo>(</mo><mi>t</mi><mo>,</mo><mi>T</mi><mo>,</mo><mi>V</mi><mo>,</mo><mi>p</mi><mo>)</mo></mrow><mi>dt</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000701785000000011.GIF" wi="1066" he="116" /></maths>单位覆冰厚度Δr表示为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>&Delta;r</mi><mo>=</mo><mfrac><mrow><mi>t</mi><msqrt><msup><mrow><mo>(</mo><mi>VW</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>P&rho;</mi><mn>0</mn></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mrow><mrow><msub><mi>&rho;</mi><mi>i</mi></msub><mi>&pi;</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000701785000000012.GIF" wi="1133" he="189" /></maths>其中,V是风速,单位为m/s,W是空气中的液态水含量,单位为g/m<sup>3</sup>,P是降雨强度,单位为mm/h,ρ<sub>0</sub>是水的密度,单位为g/cm<sup>3</sup>,ρ<sub>i</sub>是该处覆冰的密度,单位为g/cm<sup>3</sup>,t是降雨时间,单位为s;在具体运算时需要将各单位换算成国际标准单位,所以计算结果Δr的单位为m;W=aP<sup>b</sup>,其中a,b是常数,将降雨强度P与湿度转换,P满足如下关系:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>P</mi><mo>=</mo><mfrac><mn>1</mn><mrow><msub><mi>&rho;</mi><mn>0</mn></msub><mi>g</mi></mrow></mfrac><msubsup><mo>&Integral;</mo><mn>0</mn><msub><mi>p</mi><msub><mi>z</mi><mn>0</mn></msub></msub></msubsup><mi>qdp</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000701785000000013.GIF" wi="1045" he="163" /></maths>其中,q为比湿,p为气压,单位是MPa,ρ<sub>0</sub>是水的密度,单位是g/cm<sup>3</sup>,g为重力加速度,单位是m/s<sup>2</sup>,<img file="FDA0000701785000000014.GIF" wi="87" he="91" />为地面气压,单位是MPa;比湿q根据水汽压e,单位是hPa计算出来,它和水汽压的关系为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>q</mi><mo>=</mo><mfrac><mi>&epsiv;e</mi><mrow><mi>p</mi><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&epsiv;</mi><mo>)</mo></mrow><mi>e</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000701785000000021.GIF" wi="1022" he="147" /></maths>式中,<img file="FDA0000701785000000027.GIF" wi="229" he="153" />μ<sub>v</sub>,μ<sub>d</sub>分别为水汽和干空气的平均摩尔质量;采用湿度传感器得到的湿度值是一个相对湿度,其表达式为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>f</mi><mo>=</mo><mfrac><mi>e</mi><mi>E</mi></mfrac><mo>&times;</mo><mn>100</mn><mo>%</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000701785000000022.GIF" wi="1032" he="132" /></maths>其中,E为同一温度下的饱和水汽压,单位是hPa;结合公式(7)‑(10),得到单位覆冰厚度Δr关于时间t、温度T、风速V和气压p的表达式:Δr=f(t,T,V,p)   (11)将公式(11)转化为微分形式,两边积分可得到公式(13);(5)通过覆冰厚度和温湿度计算任意时刻t<sub>i</sub>导线的综合比载g<sub>5</sub>,计算公式为:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>g</mi><mn>5</mn></msub><mo>=</mo><msqrt><msup><msub><mi>g</mi><mn>3</mn></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>g</mi><mn>4</mn></msub><mn>2</mn></msup></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>20</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000701785000000023.GIF" wi="1068" he="111" /></maths>其中,g<sub>4</sub>为水平比载,g<sub>3</sub>为垂直比载;<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>g</mi><mn>3</mn></msub><mo>=</mo><msub><mi>g</mi><mn>1</mn></msub><mo>+</mo><msub><mi>g</mi><mn>2</mn></msub><mo>=</mo><mfrac><mi>gw</mi><mi>A</mi></mfrac><mo>+</mo><mfrac><mrow><msub><mi>g&rho;</mi><mi>i</mi></msub><mi>&pi;b</mi><mrow><mo>(</mo><msub><mrow><mn>2</mn><mi>r</mi></mrow><mn>0</mn></msub><mo>+</mo><mi>b</mi><mo>)</mo></mrow><mo>&times;</mo><msup><mn>10</mn><mrow><mo>-</mo><mn>3</mn></mrow></msup></mrow><mi>A</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>17</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000701785000000024.GIF" wi="1358" he="156" /></maths>其中,g<sub>1</sub>为导线比载,g<sub>2</sub>为覆冰比载,g为重力加速度,w为导线覆冰前的单位长度质量,单位是kg/m,ρ<sub>i</sub>是该处覆冰的密度,单位是g/cm<sup>3</sup>,r<sub>0</sub>为导线的半径,单位是mm,A为导线横截面积,单位是mm<sup>2</sup>;<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mfenced open='' close=''><mtable><mtr><mtd><msub><mi>g</mi><mn>4</mn></msub><mo>=</mo><mn>0.6128</mn><msub><mi>K</mi><mi>z</mi></msub><mi>&alpha;C</mi><mrow><mo>(</mo><msub><mrow><mn>2</mn><mi>r</mi></mrow><mn>0</mn></msub><mo>+</mo><mn>2</mn><mi>b</mi><mo>)</mo></mrow><mfrac><msub><mi>V</mi><mn>2</mn></msub><mi>A</mi></mfrac><mo>&times;</mo><msup><mn>10</mn><mrow><mo>-</mo><mn>3</mn></mrow></msup></mtd></mtr><mtr><mtd><mo>=</mo><mn>1.2256</mn><msub><mi>K</mi><mi>z</mi></msub><mi>&alpha;C</mi><mrow><mo>(</mo><msub><mi>r</mi><mn>0</mn></msub><mo>+</mo><mi>b</mi><mo>)</mo></mrow><mfrac><msup><mi>V</mi><mn>2</mn></msup><mi>A</mi></mfrac><mo>&times;</mo><msup><mn>10</mn><mrow><mo>-</mo><mn>3</mn></mrow></msup></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>19</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000701785000000025.GIF" wi="1321" he="344" /></maths>其中,K<sub>z</sub>是风压高度变化系数,α是风速不均匀系数,C是风载体型系数;(6)通过综合比载g<sub>5</sub>得到电力线档距端点处的轴向应力σ<sub>x</sub>,令σ<sub>x</sub>=σ<sub>max</sub>,计算断裂剩余时间t,并发出预警信号。
地址 310014 浙江省杭州市下城区朝晖六区
您可能感兴趣的专利