发明名称 空气击穿电压的预测方法与装置
摘要 本发明公开了一种空气击穿电压的预测方法与装置,采用主成分分析法对多个气象指标进行降维,再通过加权叠加,得到一维的大气放电综合指数,该指数能真实地反映各气象指标的综合特性。计算出未来某时间的大气放电综合指数后,参照往期空气击穿电压与当时大气放电综合指数的关系和规律,可以预测出未来空气击穿电压的变化趋势,且预测过程简单,预测结果准确有效。
申请公布号 CN102722660A 申请公布日期 2012.10.10
申请号 CN201210223501.7 申请日期 2012.06.29
申请人 广州供电局有限公司;华南理工大学 发明人 蚁泽沛;牛海清;张耿斌;张尧;吴倩;刘毅刚
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 广州华进联合专利商标代理有限公司 44224 代理人 王茹;曾旻辉
主权项 1.一种空气击穿电压的预测方法,其特征在于,包括:记录空气击穿电压与气象指标的历史数据,根据气象指标的历史数据计算当时的大气放电综合指数;分析历史空气击穿电压与历史大气放电综合指数的关系;获取待预测时间的气象指标的预测数据,据此计算待预测时间的大气放电综合指数;参照历史空气击穿电压与历史大气放电综合指数的关系,根据待预测时间的大气放电综合指数,预测待预测时间的空气击穿电压,大气放电综合指数的计算方法如下:建立气象指标主成分分析法数学模型:<maths num="0001"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>z</mi><mn>1</mn></msub><mo>=</mo><msub><mi>l</mi><mn>11</mn></msub><msub><mi>x</mi><mn>1</mn></msub><mo>+</mo><msub><mi>l</mi><mn>12</mn></msub><msub><mi>x</mi><mn>2</mn></msub><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><msub><mi>l</mi><mrow><mn>1</mn><mi>p</mi></mrow></msub><msub><mi>x</mi><mi>p</mi></msub></mtd></mtr><mtr><mtd><msub><mi>z</mi><mn>2</mn></msub><mo>=</mo><msub><mi>l</mi><mn>21</mn></msub><msub><mi>x</mi><mn>1</mn></msub><mo>+</mo><msub><mi>l</mi><mn>22</mn></msub><msub><mi>x</mi><mn>2</mn></msub><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><msub><mi>l</mi><mrow><mn>2</mn><mi>p</mi></mrow></msub><msub><mi>x</mi><mi>p</mi></msub></mtd></mtr><mtr><mtd><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>z</mi><mi>m</mi></msub><mo>=</mo><msub><mi>l</mi><mrow><mi>m</mi><mn>1</mn></mrow></msub><msub><mi>x</mi><mn>1</mn></msub><mo>+</mo><msub><mi>l</mi><mrow><mi>m</mi><mn>2</mn></mrow></msub><msub><mi>x</mi><mn>2</mn></msub><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><msub><mi>l</mi><mi>mp</mi></msub><msub><mi>x</mi><mi>p</mi></msub></mtd></mtr></mtable></mfenced></math>]]></maths>x<sub>1</sub>,x<sub>2</sub>,...,x<sub>p</sub>表示p个原气象指标,z<sub>1</sub>,z<sub>2</sub>,...,z<sub>m</sub>(m≤p)表示m个主成分气象指标,l<sub>ij</sub>(i=1,2,...,m;j=1,2,...,p)表示原气象指标x<sub>j</sub>在主成分气象指标z<sub>i</sub>上的载荷;计算所述数学模型的相关参数,具体过程如下:计算原气象指标相关系数矩阵<maths num="0002"><![CDATA[<math><mrow><mi>R</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>r</mi><mn>11</mn></msub></mtd><mtd><msub><mi>r</mi><mn>12</mn></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>r</mi><mrow><mn>1</mn><mi>p</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>r</mi><mn>21</mn></msub></mtd><mtd><msub><mi>r</mi><mn>22</mn></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>r</mi><mrow><mn>2</mn><mi>p</mi></mrow></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd></mtd><mtd></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd></mtd><mtd></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>r</mi><mrow><mi>p</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>r</mi><mrow><mi>p</mi><mn>2</mn></mrow></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>r</mi><mi>pp</mi></msub></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>r<sub>ij</sub>(i,j=1,2,...,p)为原气象指标x<sub>i</sub>与x<sub>j</sub>的相关系数,其计算公式为:<maths num="0003"><![CDATA[<math><mrow><msub><mi>r</mi><mi>ij</mi></msub><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><msub><mi>x</mi><mi>ki</mi></msub><mo>-</mo><mover><msub><mi>x</mi><mi>i</mi></msub><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>x</mi><mi>kj</mi></msub><mo>-</mo><mover><msub><mi>x</mi><mi>j</mi></msub><mo>&OverBar;</mo></mover><mo>)</mo></mrow></mrow><msqrt><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>ki</mi></msub><mo>-</mo><mover><msub><mi>x</mi><mi>i</mi></msub><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>kj</mi></msub><mo>-</mo><mover><msub><mi>x</mi><mi>j</mi></msub><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mfrac></mrow></math>]]></maths>计算所述数学模型中原气象指标的特征值与特征向量:通过解特征方程|λI-R|=0求出特征值λ<sub>i</sub>(i=1,2,...,m),并使其按大小顺序排列,再分别求出对应于特征值λ<sub>i</sub>的特征向量z<sub>i</sub>(i=1,2,...,m);计算主成分气象指标的贡献率及累计贡献率:主成分气象指标z<sub>i</sub>的贡献率:<maths num="0004"><![CDATA[<math><mrow><msub><mi>a</mi><mi>i</mi></msub><mo>=</mo><msub><mi>&lambda;</mi><mi>i</mi></msub><mo>/</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></munderover><msub><mi>&lambda;</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>.</mo><mo>.</mo><mo>.</mo><mi>m</mi><mo>)</mo></mrow></mrow></math>]]></maths>主成分气象指标z<sub>i</sub>的累计贡献率<img file="FDA00001833525700022.GIF" wi="291" he="299" />取累计贡献率达预定值的特征值λ<sub>1</sub>,λ<sub>2</sub>…λ<sub>m</sub>,将其对应的第一,第二,……,第m(m≤p)个主成分作为主成分气象指标;采用主成分回归法求解所述数学模型,得到大气放电综合指数:对所述数学模型,以各个主成分气象指标的贡献率为权进行分项叠加,各个主成分气象指标加权后再将其综合得到大气放电综合指数Y:<maths num="0005"><![CDATA[<math><mrow><mi>Y</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>l</mi><mrow><mi>n</mi><mn>1</mn></mrow></msub><msub><mi>a</mi><mi>n</mi></msub><msub><mi>x</mi><mn>1</mn></msub><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>l</mi><mrow><mi>n</mi><mn>2</mn></mrow></msub><msub><mi>a</mi><mi>n</mi></msub><msub><mi>x</mi><mn>2</mn></msub><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>l</mi><mrow><mi>n</mi><mn>3</mn></mrow></msub><msub><mi>a</mi><mi>n</mi></msub><msub><mi>x</mi><mn>3</mn></msub><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>l</mi><mi>np</mi></msub><msub><mi>a</mi><mi>n</mi></msub><msub><mi>x</mi><mi>p</mi></msub><mo>.</mo></mrow></math>]]></maths>
地址 510620 广东省广州市天河区天河南二路2号