发明名称 基于多元物理量数学模型的输电线路覆冰预测方法
摘要 一种基于多元物理量数学模型的输电线路覆冰预测方法,包括如下步骤:(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,通过终端发出预警。本发明实现有效监测、准确性良好。
申请公布号 CN102927949A 申请公布日期 2013.02.13
申请号 CN201210330175.X 申请日期 2012.09.07
申请人 浙江工业大学 发明人 陈庆章;李兴华;王凯;倪云峰;陈宇铮;范聪玲;吴荣杰;杨帆;宦若红;赵小敏;毛科技;雷艳静;任佳萍;丁舒羽
分类号 G01B21/08(2006.01)I 主分类号 G01B21/08(2006.01)I
代理机构 杭州天正专利事务所有限公司 33201 代理人 王兵;王利强
主权项 1.一种基于多元物理量数学模型的输电线路覆冰预测方法,其特征在于:所述预测方法包括如下步骤:(1)采用温湿度传感器、风力传感器获取当前时刻的温湿度及风力气象参数;(2)分别判断当前各气象数据是否满足覆冰条件,即:温度-8~0℃,湿度90%RH以上,风速2~7m/s,若满足条件则转到(3),否则进入下一次采样周期,即返回到(1);(3)通过张力传感器获得的数据值与初始值做比较,若<img file="FDA00002111904600011.GIF" wi="377" he="92" />说明导线有覆冰,转到步骤(4),否则返回到(1);(4)通过地面终端获取当前环境气压值,计算覆冰厚度b,计算公式如下:<maths num="0001"><![CDATA[<math><mrow><mi>b</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>]]></maths>单位覆冰厚度Δr表示为:<maths num="0002"><![CDATA[<math><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><mi>P</mi><msub><mi>&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>]]></maths>其中,V是风速(m/s),W是空气中的液态水含量(g/m<sup>3</sup>),t是降雨时间(s);P是降雨强度(mm/h),ρ<sub>0</sub>是水的密度(g/cm<sup>3</sup>),t是降雨时间(s);W=aP<sup>b</sup>,其中a,b是常数,将降雨强度P与湿度转换,P满足如下关系:<maths num="0003"><![CDATA[<math><mrow><mi>P</mi><mo>=</mo><mfrac><mn>1</mn><mi>&rho;g</mi></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>]]></maths>其中,q为比湿,p为气压(MPa),ρ为液态水密度(g/cm<sup>3</sup>),g为重力加速度(m/s<sup>2</sup>),<img file="FDA00002111904600015.GIF" wi="63" he="53" />为地面气压(MPa);比湿q根据水汽压e(hPa)计算出来,它和水汽压的关系为:<maths num="0004"><![CDATA[<math><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>]]></maths>式中,<img file="FDA00002111904600017.GIF" wi="180" he="125" />μ<sub>v</sub>,μ<sub>d</sub>分别为水汽和干空气的平均摩尔质量;采用湿度传感器得到的湿度值是一个相对湿度,其表达式为:<maths num="0005"><![CDATA[<math><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>]]></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"><![CDATA[<math><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>]]></maths>其中,g<sub>4</sub>为水平比载,g<sub>3</sub>为垂直比载;<maths num="0007"><![CDATA[<math><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>gq</mi><mi>A</mi></mfrac><mo>+</mo><mfrac><mrow><mi>g&rho;&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>]]></maths>其中,g<sub>1</sub>为导线比载,g<sub>2</sub>为覆冰比载,g为重力加速度,q为单位长度质量(kg/m),ρ是覆冰的密度(g/cm<sup>3</sup>),r<sub>0</sub>为导线的半径(mm),A为导线横截面积(mm<sup>2</sup>);<maths num="0008"><![CDATA[<math><mrow><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><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></mrow></math>]]></maths><maths num="0009"><![CDATA[<math><mrow><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><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>19</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,K<sub>z</sub>是风压高度变化系数,α是风速不均匀系数,C是风载体型系数;(6)通过综合比载g<sub>5</sub>得到电力线档距端点处的轴向应力σ<sub>x</sub>,令σ<sub>x</sub>=σ<sub>max</sub>,计算断裂剩余时间t,并发出预警信号。
地址 310014 浙江省杭州市下城区朝晖六区