发明名称 一种覆冰输电导线倾角预测计算方法
摘要 本发明提出了一种覆冰输电导线倾角的预测计算方法。首先基于覆冰输电导线倾角的历史数据构建其差分时间序列,然后对差分时间序列建立自回归计量模型并运用最终预测误差准则估算模型阶p,运用岭回归方法估算模型参数,最后应用建立的模型预测计算覆冰输电线倾角在未来若干小时内的数值。此方法可动态预测覆冰输电导线的倾角,具有实时性强和精确度高的特点,有效解决了覆冰输电导线倾角预测的技术难题,为输电塔线体系的冰致失效分析提供了数据支持。
申请公布号 CN104751003A 申请公布日期 2015.07.01
申请号 CN201510175647.2 申请日期 2015.04.15
申请人 华北电力大学 发明人 武鑫;苏睿;芮晓明
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 代理人
主权项 一种覆冰输电导线倾角的预测计算方法,其具体步骤如下:1)设一个n期时间序列的覆冰输电导线的倾角数据为X={x<sub>1</sub> x<sub>2</sub> … x<sub>n</sub>};2)X的一阶差分序列为Y={y<sub>1</sub> y<sub>2</sub> … y<sub>n‑1</sub>},y<sub>i</sub>=x<sub>i+1</sub>‑x<sub>i</sub> 1≤i≤n‑1;3)差分序列Y的自回归模型为:y<sub>t</sub>=θ<sub>1</sub>y<sub>t‑1</sub>+θ<sub>2</sub>y<sub>t‑2</sub>+…+θ<sub>p</sub>y<sub>t‑p</sub>+ξ<sub>t</sub>式中p为模型的滞后期;t为时间序列Y中的时间刻度,可取大于等于p且小于等于n‑1的任意自然数;y<sub>t</sub>,y<sub>t‑1</sub>,…,y<sub>t‑p</sub>为时间序列Y中任意连续的p个元素;θ<sub>1</sub>,θ<sub>2</sub>,…,θ<sub>p</sub>为自回归参数;ξ<sub>t</sub>为误差项;4)对于模型的滞后期或模型阶数p,可以采用最终预测误差(FPE)准则判断;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>FPE</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msup><mover><mi>&sigma;</mi><mo>^</mo></mover><mn>2</mn></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mrow><mo>(</mo><mfrac><mrow><mi>n</mi><mo>-</mo><mn>1</mn><mo>+</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn><mo>-</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></mfrac><mo>)</mo></mrow><mo>,</mo><mi>k</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mi>p</mi></mrow>]]></math><img file="FSA0000115672560000011.GIF" wi="1129" he="159" /></maths>式中<img file="FSA0000115672560000014.GIF" wi="142" he="85" />为取k=1,2,…,p时的样本误差方差σ<sup>2</sup>的估计;<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>FPE</mi><mrow><mo>(</mo><mover><mi>p</mi><mo>^</mo></mover><mo>)</mo></mrow><munder><mi>min</mi><mrow><mn>1</mn><mo>&le;</mo><mi>k</mi><mo>&le;</mo><mi>p</mi></mrow></munder><mi>FPE</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000115672560000012.GIF" wi="575" he="120" /></maths>即为所求模型阶数p的最终预测误差准则估计;5)自回归模型中的自回归参数,θ<sub>1</sub>,θ<sub>2</sub>,…,θ<sub>p</sub>,可采用岭回归方法获得其估计值<img file="FSA0000115672560000013.GIF" wi="313" he="114" />采用的岭回归方法如下所示:差分序列Y的自回归模型可改写为:y<sub>t</sub>=Φ·θ+ξ<sub>t</sub>式中θ=[θ<sub>1</sub> θ<sub>2</sub> … θ<sub>p</sub>]<sup>T</sup>Φ=[y<sub>t‑1</sub> y<sub>t‑2</sub> … y<sub>t‑p</sub>]通过岭回归方法对自回归参数的预估公式为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mover><mi>&theta;</mi><mo>^</mo></mover><mo>=</mo><msup><mi>&Phi;</mi><mo>+</mo></msup><mo>[</mo><mi>Q</mi><mo>]</mo></mrow>]]></math><img file="FSA0000115672560000015.GIF" wi="275" he="127" /></maths>式中Φ<sup>+</sup>=Φ<sup>T</sup>[ΦΦ<sup>T</sup>+αI]<sup>‑1</sup>,Φ<sup>T</sup>为Φ的转置矩阵,I为单位矩阵,α为偏常数(0≤α<∞);通过取不同的α值,作出y<sub>t</sub>(α)的岭迹图,进行岭迹分析,求得优化的α值;6)第n+1,…,n+j期输电导线覆冰倾角的预测值<img file="FSA0000115672560000016.GIF" wi="338" he="84" />可用下述公式计算:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mover><mi>y</mi><mo>^</mo></mover><mi>n</mi></msub><mo>=</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>x</mi><mi>n</mi></msub><mo>=</mo><msub><mover><mi>&theta;</mi><mo>^</mo></mover><mn>1</mn></msub><msub><mi>y</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mover><mi>&theta;</mi><mo>^</mo></mover><mn>2</mn></msub><msub><mi>y</mi><mrow><mi>n</mi><mo>-</mo><mn>2</mn></mrow></msub><mo>+</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>+</mo><msub><mover><mi>&theta;</mi><mo>^</mo></mover><mi>p</mi></msub><msub><mi>y</mi><mrow><mi>n</mi><mo>-</mo><mi>p</mi></mrow></msub><mo>&DoubleRightArrow;</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>x</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mi>x</mi><mi>n</mi></msub><mo>+</mo><msub><mover><mi>&theta;</mi><mo>^</mo></mover><mn>1</mn></msub><msub><mi>y</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mover><mi>&theta;</mi><mo>^</mo></mover><mn>2</mn></msub><msub><mi>y</mi><mrow><mi>n</mi><mo>-</mo><mn>2</mn></mrow></msub><mo>+</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>+</mo><msub><mover><mi>&theta;</mi><mo>^</mo></mover><mi>p</mi></msub><msub><mi>y</mi><mrow><mi>n</mi><mo>-</mo><mi>p</mi></mrow></msub></mrow>]]></math><img file="FSA0000115672560000021.GIF" wi="2056" he="95" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>n</mi><mo>+</mo><mi>j</mi></mrow></msub><mo>=</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>n</mi><mo>+</mo><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mover><mi>&theta;</mi><mo>^</mo></mover><mn>1</mn></msub><msub><mover><mi>y</mi><mo>^</mo></mover><mrow><mi>n</mi><mo>+</mo><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mover><mi>&theta;</mi><mo>^</mo></mover><mn>2</mn></msub><msub><mover><mi>y</mi><mo>^</mo></mover><mrow><mi>n</mi><mo>+</mo><mi>j</mi><mo>-</mo><mn>2</mn></mrow></msub><mo>+</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>+</mo><msub><mover><mi>&theta;</mi><mo>^</mo></mover><mi>p</mi></msub><msub><mover><mi>y</mi><mo>^</mo></mover><mrow><mi>n</mi><mo>+</mo><mi>j</mi><mo>-</mo><mi>p</mi></mrow></msub></mrow>]]></math><img file="FSA0000115672560000022.GIF" wi="1071" he="112" /></maths>式中j>n;“^”代表预测值;<img file="FSA0000115672560000023.GIF" wi="468" he="101" />
地址 102206 北京市昌平区回龙观北农路2号华北电力大学