一种基于相关向量机的输电线路覆冰预测方法技术

本发明专利技术属于电力系统灾害预警技术领域，具体涉及一种基于相关向量机的输电线路覆冰预测方法。本预测方法根据覆冰现象的特点，具有针对性的选择和处理覆冰预测模型的输入量及其权重指标；采用相关向量机方法建立输电线路覆冰预测模型；采用样本数据对模型进行训练，并利用量子粒子群算法和K-fold交叉验证法对模型进行优化；根据测试数据对输电线路的覆冰厚度及概率分布进行预测，并进一步通过重复训练对模型进行修正以提高预测精度。该预测方法综合考虑了多种因素对输电线路覆冰的影响，能够准确预测输电线路的覆冰厚度，具有很强的预测精度和泛化能力。

【技术实现步骤摘要】

本专利技术属于电力系统灾害预警
，具体涉及一种基于相关向量机的输电线 路覆冰预测方法。
技术介绍
输电线路覆冰厚度超过其设计标准，会造成闪络、跳闸甚至断线、倒塔等严重事故 的发生。近年来，因覆冰导致电网事故的次数越来越多，给社会和人民财产带来巨大损失。 因此，亟需对输电线路的覆冰预测方法进行研究，为电力部门的抗冰减灾工作提供预警等 决策支持。 目前，国内外关于导线覆冰预测主要有物理模型和经验模型两种。由于存在测量 技术限制等原因，物理模型中所需的部分信息在实际线路中难以获取或精度不高，很难直 接应用于实际线路的覆冰预测；经验模型基于模糊逻辑理论、支持向量机等方法实现覆冰 预测，但存在着泛化能力较差、预测精度较低等缺点。 相关向量机（Relevance Vector Machine, RVM)是一种基于贝叶斯框架下的稀疏 概率模型。作为机器学习领域的研究热点之一，RVM较之支持向量机等其它智能类算法具 有更强的泛化能力和更高的预测精度。因此，本专利技术的一种基于相关向量机的输电线路覆 冰预测方法，其建模方便、预测结果精确、泛化能力强等特点，具有重要推广应用价值。
技术实现思路
本专利技术针对现有技术的不足，提出了一种基于相关向量机的输电线路覆冰预测方 法。 本预测方法根据覆冰现象的特点，具有针对性的选择和处理覆冰预测模型的输入 量及其权重指标；采用相关向量机方法建立输电线路覆冰预测模型；采用样本数据对模型 进行训练，并利用量子粒子群算法和K-fold交叉验证法对模型进行优化；根据测试数据对 输电线路的覆冰厚度及概率分布进行预测，并进一步通过重复训练对模型进行修正以提高 预测精度。该预测方法综合考虑了多种因素对输电线路覆冰的影响，能够准确预测输电线 路的覆冰厚度，具有很强的预测精度和泛化能力。 本专利技术的，包括以下步骤： -种基于相关向量机的输电线路覆冰预测方法，其特征在于，包括以下步骤： 步骤1，根据输电线路覆冰的物理规律，选择考虑权重的覆冰历史数据和微气象数 据作为输入数据，覆冰厚度作为输出数据，并进行归一化处理； 其中，考虑权重的覆冰历史数据定义为Yh，基于以下公式； [0011 ] Yk-m = [gi (Wh) X yk-m，· · ·，gn (Wh) X yk_ J 其中：yk_m为k-m时刻的覆冰厚度的历史数据；m为覆冰历史数据和覆冰预测值之 间的时间尺度，由采样周期决定，单位为15min ;n为选择覆冰历史数据的个数，为用户自定 义，其值是大于〇的正整数;wh为权重变化率；gi (Wh)为第i个覆冰历史数据的权重； 其中，第i个覆冰历史数据的权重表征为： gj(wh) = l~iXwh 微气象数据包括温度、温度变化率和相对湿度，用uk_m表示，具体指代k-m时刻的 微气象数据； 覆冰厚度yk用表示，具体指代k时刻的覆冰厚度预测值，单位为厘米； 数据的归一化处理，基于以下公式： X* = (X-Xmin) / (Xmax-Xmin) 其中：X为上述原始的输入、输出数据；X#为经过归一化处理的数据；X min表示X的 最小值；Xniax表示X的最大值。 步骤2,基于相关向量机的预测原理，建立输电线路的覆冰预测模型，基于以下公 式：本文档来自技高网...

【技术保护点】
一种基于相关向量机的输电线路覆冰预测方法，其特征在于，包括以下步骤：步骤1，根据输电线路覆冰的物理规律，选择考虑权重的覆冰历史数据和微气象数据作为输入数据，覆冰厚度作为输出数据，并进行归一化处理；其中，考虑权重的覆冰历史数据定义为Yk‑m，基于以下公式；Yk‑m＝[g1(wh)×yk‑m,...,gn(wh)×yk‑m‑n+1]其中：yk‑m为k‑m时刻的覆冰厚度的历史数据；m为覆冰历史数据和覆冰预测值之间的时间尺度，由采样周期决定，单位为15min；n为选择覆冰历史数据的个数，为用户自定义，其值是大于0的正整数；wh为权重变化率；gi(wh)为第i个覆冰历史数据的权重；其中，第i个覆冰历史数据的权重表征为：gi(wh)＝1‑i×wh微气象数据包括温度、温度变化率和相对湿度，用uk‑m表示，具体指代k‑m时刻的微气象数据；覆冰厚度yk用表示,具体指代k时刻的覆冰厚度预测值，单位为厘米；数据的归一化处理，基于以下公式：x*＝(x‑xmin)/(xmax‑xmin)其中：x为上述原始的输入、输出数据；x*为经过归一化处理的数据；xmin表示x的最小值；xmax表示x的最大值；步骤2，基于相关向量机的预测原理，建立输电线路的覆冰预测模型，基于以下公式：yk=Σi=1NwiK(xk,xi)+w0]]>其中：yk为k时刻的覆冰厚度预测值；wi为模型的权重值；xk为k时刻对应的输入量，并且有xk＝(Yk‑m,uk‑m)；K(xk,xi)为核函数；其中，核函数K(xk,xi)表征为：K(xk,xi)=exp(-||xk-xi||2δ2)]]>其中：δ2为核函数的宽度；步骤3，针对覆冰预测模型及输入变量中的待定参数，基于量子粒子群算法和交叉检验法进行辨识，实现模型优化；覆冰预测模型及输入变量中的待定参数包括权重变化率wh和核函数宽度δ2；wh的取值范围为[0，0.2]，δ2的取值范围为(0，1]；基于量子粒子群算法和交叉检验法的参数辨识，其适应度函数f基于以下公式：f=Σi=1mei]]>其中：m为交叉验证法的折数；ei为每折的均方根误差；其中，每折的均方根误差表征为：ei=Σj=1nΔyj2/n]]>其中：n为每折的样本数量；Δyj为样本j的覆冰厚度预测误差；步骤4，根据步骤1中的输入输出数据和步骤3得到的优化参数对步骤2建立的输电线路的覆冰预测模型进行训练，得到模型的后验概率分布均值μ，后验协方差矩阵Σ以及超参数αMP和σMP2；步骤5，对输电线路的覆冰厚度及概率分布进行预测，基于以下公式：p(t*|t,αMP,σMP2)=∫p(t*|w,σMP2)p(w|t,αMP,σMP2)dw=N(t*|y*,σ*2)]]>其中：t*为需要预测的覆冰厚度；αMP和为步骤4中得到的超参数α和σ2；y*为预测均值，可以作为t*的预测结果；σ*2为预测方差；其中，预测均值可以表征为：其中：为核函数向量；其中，预测方差可以表征为：步骤6，重复步骤4，通过重复训练的方法对覆冰预测模型的权重向量进行修正具体方法是：将最新实测得到的微气象数据和对应的覆冰厚度历史数据作为训练样本代入模型进行重复训练，对模型原有的权重向量w进行修正，并利用修正的权重向量w*对覆冰厚度进行预测。...

【技术特征摘要】
1. 一种基于相关向量机的输电线路覆冰预测方法，其特征在于，包括以下步骤： 步骤1，根据输电线路覆冰的物理规律，选择考虑权重的覆冰历史数据和微气象数据作 为输入数据，覆冰厚度作为输出数据，并进行归一化处理； 其中，考虑权重的覆冰历史数据定义为Yk_m，基于以下公式； Yk_m = Lg1 (wh) Xyk_m, . . . , gn(wh) Xyk_m_n+1] 其中：yk-m为k-m时刻的覆冰厚度的历史数据；m为覆冰历史数据和覆冰预测值之间的 时间尺度，由采样周期决定，单位为15min ;n为选择覆冰历史数据的个数，为用户自定义， 其值是大于〇的正整数;wh为权重变化率；gi (Wh)为第i个覆冰历史数据的权重； 其中，第i个覆冰历史数据的权重表征为： gi (wh) = l_i Xwh 微气象数据包括温度、温度变化率和相对湿度，用uk_m表示，具体指代k-m时刻的微气 象数据； 覆冰厚度yk用表示，具体指代k时刻的覆冰厚度预测值，单位为厘米； 数据的归一化处理，基于以下公式： X* = (X-Xmin) / (Xmax-Xmin) 其中：x为上述原始的输入、输出数据；x#为经过归一化处理的数据 ；Xmin表示X的最小 值；Xnlax表示X的最大值； 步骤2,基于相关向量机的预测原理，建立输电线路的覆冰预测模型，基于以下公式：其中：yk为k时刻的覆冰厚度预测值;Wi为模型的权重值；Xk为k时刻对应的输入量， 并且有 xk = (YkI, uk_J ;K (xk, Xi)为核函数； 其中，核函数K (xk，Xi)表征为：其中：S 2为核函数的宽度； 步骤3,针对覆冰预测模型及输入变量中的待定参数，基于量子粒子群算法和交叉检验 法进行辨识，实现模型优化； 覆冰预测模型及输入变量中的待定参数包括权重变化率Wh和核函数宽度S 2;Wh的取 值范围为[〇,〇.2]，S2的取值范围为（〇，1]; 基于量子粒子群算法和交叉检验法的参数辨识，其适应度函数f基于以下公式：其中：m为交叉验证法的折数；ei为每折的均方根误差； 其中，每折的均方根误差表征为：其中：n为每折的样本数量；△ h为样本j的覆冰厚度预测误差； 步骤4,根据步骤1中的输入输出数据和步骤3得...

【专利技术属性】
技术研发人员：赵洁，王骏，刘涤尘，王力，刘田，赵语，贾骏，唐飞，
申请(专利权)人：武汉大学，
类型：发明
国别省市：湖北;42