专利名称：一种区域电网污闪指数预测方法
技术领域：
本发明属于输配电监测技术领域，特别涉及一种区域电网污闪指数预测方法。
背景技术：
电网在运行过程中，由于绝缘子表面发生污秽闪络对供电系统有很不利影响，目 前，针对绝缘子表面污秽闪络的预测方法有很多，主要集中在对于绝缘子表面等值附盐密 度的预测、绝缘子污秽闪络电压的预测、及基于泄漏电流的污秽闪络预测方面。其中，等值盐密的神经网络预测、污闪电压的支持向量机预测都是将描述污闪的 这两个参数与环境气象参数之间的多维非线性关系用非线性建模方法来进行拟合，这种预 测方法在国内外均处于领先地位，但这两种方法仍存在建模对象单一，建立的预测模型及 其预测功能与电网防污工作要求之间存在一定差异，尚无法在电网的污闪防治运行管理中 得到实际应用，基于泄漏电流的污闪预测方法的基础数据丰富，在此基础上应用人工神经 网络等方法能够取得与工程实际比较接近的拟合效果，但泄漏电流的监测点分布范围有 限，无法对电网的污闪状态进行较为全面的实时监测和预测。
发明内容
针对现有技术的不足，本发明提供一种区域电网污闪指数预测方法。本发明一种区域电网污闪指数预测方法，步骤如下步骤1、应用绝缘子表面等值盐密预测模型，实时预测等值盐密当前值；将电网现场采集的环境参数和等值盐密历史值输入绝缘子表面等值盐密预测模 型，绝缘子表面等值盐密预测模型输出即为实时预测等值盐密当前值；步骤2、应用绝缘子污闪临界电压预测模型，预测污闪临界电压；将实时预测等值盐密当前值和当前采集环境参数输入绝缘子的污闪临界电压预 测模型，绝缘子的污闪临界电压预测模型的输出即为污闪临界电压预测值；步骤3、应用污闪分级预测预警模型，预测电网绝缘子污闪指数；将污闪临界电压预测值输入污闪分级预测预警模型，污闪分级预测预警模型的输 出即为预测的电网污闪指数；步骤4、当电网污闪指数为0和5 %时，不进行污闪预警；当电网污闪指数为20 % 时，发布污闪III级预警；当电网污闪指数为50%和85%时，污闪发生概率已经大于50%， 发布污闪II级预警；当电网污闪指数为100%时，污闪发生概率已经相当大，该区域电网随 时有可能发生污闪，发布污闪I级预警。所述的绝缘子表面等值盐密预测模型的应用，按如下步骤进行1)建立多变量等值盐密时间序列；在固定时间间隔对等值盐密数据进行测量，在一系列时刻、，t2，. . .，&得到的离 散有序集合{Xl，x2, . . . , xj称为离散等值盐密时间序列，简称为等值盐密时间序列；多变量等值盐密时间序列是由等值盐密时间序列及在相同时刻内的气象参数时间序列构成的多维等值盐密时间序列，是由包括等值盐密时间序列在内的多维等值盐密非 线性动力学系统的表现形式； M维等值盐密时间序列X1;X2，...，Xn，N代表时刻数，其中Xi = (xlji; X2ji,... ,χΜ, D，即 式中，i = 1，2，...，N，xm,n表示第M个变量在N时刻的数值，Xfti表示第M个变量 在i时刻的数值；2)重构多变量等值盐密时间序列的相空间多变量等值盐密时间序列相空间重构的相点为 XM,N-(n,M-i)rM表示第M个变量在N时刻在延迟时间为τ Μ嵌入维数为mM的重构相空 间中的值；其中！！表示第！！时亥！！，”^；^恐丨叫-!)^+!, τ i*!!^为第i个时间序列的延迟时
间和嵌入维数，重构相空间的嵌入维数m = HI^m2+. . . +mM, M为时间序列的维数；多变量等值盐密时间序列的相空间重构参数延迟时间τ的选择采用互信 息法，互信息法是以互信息第一次达到最小时的延时作为相空间重构的延迟时间，由
是等值盐密时间序列时间跨度为(i + Ι) τ的自
相关函数，τ为相空间重构参数延迟时间；嵌入维数m由 决定，其中 Xi (m+1)是(m+1)维重构的等值盐密系统相空间中的第i个相点，n(i，m)是在m 维等值盐密系统相空间中使相点Xn(i,m)(m)是相点XiOii)的最邻近点的整数，I I · I I是等值 盐密系统相空间上的欧式距离；3)等值盐密时间序列确定性检验本发明中采用李雅谱诺夫指数法进行等值盐密时间序列的确定性检验，该指数是 相空间中邻近轨道的平均指数发散率的数值表征，用以刻画混沌运动的初态敏感性，该指数作为沿轨道长期平均的结果，是一种整体特征，其值总是实数；判断等值盐密时间序列的非线性特性通过计算最大李雅普诺夫指数得到，该方法 计算由yGO曲线的回归直线斜率为最大指数，其中，
li(k)表示对重构的等值盐密相空间中每一对最邻近点，计算k个离散时间后的欧式距离， M为时间序列的维数；4)全局预测多变量等值盐密时间序列根据泰肯斯延时嵌入定理，只要嵌入维数m和延迟时间T选择合理，重构相空 间在嵌入空间的轨迹就与微分同胚意义下的等值盐密动力学系统等价，且存在光滑映射 fmr’使得
表示重构相空间中第i+i个相点，应用非线性逼近方
法构造映射~来近似逼近f，并使~满足
最小，其中
表示表示第M个变量在n时刻在延迟时间为T M嵌入维数
为mM的重构相空间中的值，、表示第M个变量的延迟时间，mM表示第M个变量的嵌入维 数；5)利用支持向量机模型求解多变量等值盐密时间序列预测模型通过求解预测模型中的非线性映射p，确定等值盐密时间序列预测模型，并使预测 模型在求得的非线性映射7■下的预测误差满足要求，支持向量机理论可以有效地解决等值 盐密数据样本容量偏小的情况下的等值盐密非线性时间序列预测模型中非线性映射y 的求 解问题，用于逼近等值盐密时间序列预测模型中非线性映射关系的支持向量机方法是支持 向量回归；设等值盐密系统相空间相点构成的样本集为S = {(Xi,yi)，i = 1，2，. . .，M}，(Xi， Yi)表示重构相空间中的任一相点，若存在一个超平面g(x) = <wx>+b,w G Rn，b G R，w、 b表示向量参数，为了构造超平面g(x)，使得|yi-g(Xi)|彡e成立，其中，〈 表示向量内 积，i = 1，2，...^^为等值盐密时间序列的维数，则样本集3= {(Xi，yi)，i = 1,2,...,
|〈W ；C〉+ 6— 只 | ￡
M}为 ￡ 的近似集，有|<W.x>+b-yi| ( ￡，即 其中
为S的点到超平面f (x)的距离屯，则有
，即集合S中的点到超平面的距离最大值为
通过最大化S中的点到
超平面距离的上界可得到集合S的最优近似超平面，则最优近似超平面可通过最大化式
得到，因此求解I !wl I2的最小化问题即可得到集合S的最优近似超平面，由于 等值盐密系统是非线性系统，必须用一个非线性映射识把等值盐密系统相空间中的相点Xi映射到一个高维空间，然后在高维空间里进行线性回归，由于优化过程中涉及到高维空间 的内积运算，为了避免内积运算，用核函数W (Xi，Xi+1)代替内积〈炉xw)〉来实现等值 盐密系统相空间中非线性回归，此时，等值盐密系统相空间上的支持向量回归问题可转化 为如下的I kl I2优化问题 其中，i = 1，2，. . .，M，上式为二次规划问题，其Lagrange函数为 其中，a i和被称为拉格朗日乘子，对任何i = 1，2，. . .，M，都有等式a,_ xa/ = 0， a i >0，a,20成立；在进行等值盐密系统相空间中非线性映射函数逼近时，由于求得的回归函数与实 际函数之间不可避免的存在误差，因此引入松弛变量I i 彡 0,^；>0, i = 1，2， ，M，I 丄表示松弛变量；此时的优化为 式中c为惩罚参数，且c > 0 ；可得Lagrange对偶问题为 上式作为绝缘子表面等值盐密的预测；求解上式即可得等值盐密系统相空间中非线性映射的回归函数g(x)，核函 数与作为输入量，即可得到下一时刻的等值盐密值，核函数中的参数和支持向量机模型 参数中的惩罚因子是决定支持向量机方法求得的非线性映射^ 预测性能的最主要因素，采 用交叉验证法进行支持向量机模型参数的选择。
所述的绝缘子的污闪临界电压预测模型的应用，按如下步骤进行1)污闪临界电压的确定污闪临界电压值选用50%污闪电压试验获得的数据U5(l%，对等值盐密及气象参数 与污闪临界电压值U5(l%之间的多维非线性关系进行人工神经网络建模，实现对已知等值盐 密及气象参数的绝缘子污闪临界电压队⑽进行预测；2)人工神经网络模型的建立本发明中采用基于BP方法的污闪临界电压预测人工神经网络模型；确定污闪临界电压预测BP人工神经网络模型的输入层、隐层和输出层三层；污闪临界电压预测BP人工神经网络模型包括输入层、隐层和输出层三层结构，根 据不同的人工污秽实验数据，隐层分别为两层或三层节点组成；X1, X2, ... , Xn为输入层节点，包括等值盐密值、温度、湿度、气压、风速和雨量，hi; h2，. . .，hm为隐层节点，ο为输出层节点，是绝缘子污闪临界电压预测值W1，V2,..., Vffl为输 入层至隐层的权值，WnW2,...，WmS隐层至输出层权值，污闪临界电压人工神经网络模型的 训练使用人工污秽耐压试验数据作为训练数据，整个网络的结构和参数由训练进行优化， 网络的输入层维数为n，输出维数为一维即相关预测值的输出，隐层单元维数m在网络训练 学习中优化确定，预测模型的预测值与实际人工污秽耐压试验的结果进行比较，进而对网 络结构和权值进行修正和完善；3)在人工神经网络模型的基础上增设动量项并对对BP方法的学习率进行自适应 调整为提高污闪临界电压预测人工神经网络的训练速度，在权值调整公式中增加一个 动量项，若用W代表污闪临界电压预测人工神经网络中某层的权值矩阵，X代表某层输入向 量，则含有动量项的污闪临界电压预测人工神经网络权值调整向量表达式为Aff(t) = η δ χ+α Δ W(t_l)，公式中字母表示意思α是动量系数，设α e (0， 1)，η是神经网络的学习率，AW(t-l)是前一次的权值调整量，AW(t)是本次的权值调整 量，动量项反映了以前的调整经验，对于t时刻的调整起阻尼作用，当误差曲面出现骤然起 伏时，可减小振荡趋势，提高训练速度；在污闪临界电压预测人工神经网络建模中对BP方法的学习率进行自适应调整， 学习率η e (0,1)表示比例系数，设一个初始学习率，若经过一次权值调整后使总误差增 加，则本次调整无效；若经过一次权值调整后使总误差下降，则本次调整有效；4)在人工神经网络模型的基础上引入陡度因子，使权值调整脱离平坦区在污闪临界电压预测人工神经网络模型训练过程中引入陡度因子，误差曲面上存 在着平坦区域，权值调整进入平坦区代表是污闪临界电压预测人工神经网络的神经元输出 进入了转移函数的饱和区，如果在进入平坦区后，设法压缩神经元的净输入，使其输出并退 出转移函数的饱和区，就可以改变误差函数的形状，从而使调整脱离平坦区，具体做法是， 在原转移函数中引入一个陡度因子ζ，使输出为
net
γ πO --
net
l + e ζ其中，net为各层节点的输出值，当发现ΔE接近零而d_o值仍较大时，即认为进入平坦区，此时令(> 1 ；当退出平坦区后，再令4 = 1，当4 > 1时，net坐标压缩了倍， 污闪临界电压预测人工神经网络的神经元的转移函数曲线的敏感段变长，从而使net值退 出饱和值，当(=1时，转移函数恢复原状，对较小的net值具有较高的灵敏度；建立改进的BP人工神经网络预测模型，求解公式0(j) = f(net(j))即可得到绝 缘子的污闪临界电压预测值，其中j = 1，2，. . .，1，1为输出层的节点数，net(j)表示污闪 临界电压预测神经元j的输入总和，f表示一种非线性映射关系。所述的污闪分级预测预警模型的应用，按如下步骤进行确定绝缘子的污闪指数，污闪指数预测模型通过公式 「0 UFF>l.UJOP 5% 20% 50% 85% 100%
值；Um为运行电压，污闪临界电压的预测值为该模型的输入，K为该模型的输出，即为预测 的污闪指数。本发明的优点提出了基于相空间重构的多变量等值盐密时间序列的预测模型， 并且用支持向量机模型对其求解，解决了在等值盐密数据小样本且存在噪声情况下的预测 问题，提高了预测精度。
K -
1.05Uop<Uff <\.\Uop 1.0Uop<Uff<1.05Uop 0.9WOP < UFF <l.0UOP 0.95UOP < UFF < 0.98UOP 0.S5UOP < UFF < 0.95UOP
决定，其中K为污闪指数；UFF为污闪临界电压的预测
图1 (a)本发明区域电网污闪指数预测方法流程图；图1 (b)本发明应用绝缘子表面等值盐密的预测模型流程图；图1 (c)本发明应用绝缘子的污闪临界电压预测模型流程图；图2本发明人工神经网络模型图；图3本发明等值盐密监测点等值盐密测量值与计算值比较；图4本发明污闪临界电压测量值与计算值比较图；图5本发明污闪指数预测预警模型图。
具体实施例方式本发明一种区域电网污闪指数预测方法结合实施例和附图加以说明。该方法，步骤如下，如图1 (a)所示步骤1、应用绝缘子表面等值盐密预测模型，实时预测等值盐密当前值；将电网现场采集的环境参数和等值盐密历史值输入绝缘子表面等值盐密预测模 型，绝缘子表面等值盐密预测模型输出即为实时预测等值盐密当前值；其中现场采集的环境参数包括风速值、温度值、气压值、雨量值、湿度值和等值盐 密历史值；风速、温度、气压、降雨量、湿度的测量数据和等值盐密历史值构成一个六维的等 值盐密多变量时间序列，选取这六个变量的同一时间段的20个不同时刻值构成时间序列具体参数见下式
12
上式中第一行为风速值，单位为m/s ；第二行为温度值，单位为。C ；第三行为气压 值，单位为hPa ；第四行为雨量值，单位为mm ；第五行为相对湿度值，单位为％ ；第六行为等值盐密值，单位为mg/cm2。步骤2、应用污闪临界电压预测模型，预测污闪临界电压；将实时预测等值盐密当前值和当前采集环境参数输入绝缘子的污闪临界电压预 测模型，绝缘子的污闪临界电压预测模型的输出即为污闪临界电压预测值；步骤3、应用污闪分级预测预警模型，预测电网绝缘子污闪指数；将污闪临界电压预测值输入污闪分级预测预警模型，污闪分级预测预警模型的输 出即为预测的电网污闪指数；步骤4、当电网污闪指数为0和5 %时，不进行污闪预警；当电网污闪指数为20 % 时，发布污闪III级预警；当电网污闪指数为50%和85%时，污闪发生概率已经大于50%， 发布污闪II级预警；当电网污闪指数为100%时，污闪发生概率已经相当大，该区域电网随 时有可能发生污闪，发布污闪I级预警，如图5所示。污闪预测模型总的输入条件为气象数据，包括历史记录数据和预报，输出为污秽 闪络电压等级，根据污秽闪络电压等级即可对污秽闪络发生风险做出判断。所述的绝缘子表面等值盐密预测模型的应用，按如下步骤进行，如图1(b)所示1)建立多变量等值盐密时间序列；其中现场采集的环境参数包括风速值、温度值、气压值、雨量值、湿度值和等值盐 密历史值；风速、温度、气压、降雨量、湿度的测量数据和等值盐密历史值构成一个六维的等 值盐密多变量时间序列，选取这六个变量的同一时间段的20个不同时刻值构成时间序列具体参数见下式
上式中第一行为风速值，单位为m/s ；第二行为温度值，单位为。C ；第三行为气压 值，单位为hPa ；第四行为雨量值，单位为mm ；第五行为相对湿度值，单位为％ ；第六行为等值盐密值，单位为mg/cm2。2)重构多变量等值盐密时间序列的相空间以延迟时间t = 3和嵌入维数m = 6对等值盐密多变量时间序列进行相空间重 构，在重构的等值盐密相空间中，以相空间中的所有相点构成等值盐密时间序列支持向量 机模型的训练样本，建立支持向量机模型，对等值盐密时间序列全局预测模型中的非线性 映射7 进行拟合；相空间中的全部相点构成的训练样本为 其中相点总数为N = 6(n-(m-l) τ )；3)等值盐密时间序列确定性检验通过改进方法计算最大Lyapimov指数，结果约等于0. 047，由最大Lyapunov指数 大于零可判断等值盐密时间序列为非线性混沌时间序列；4)全局预测多变量等值盐密时间序列基于以上的三个步骤，建立了多变量等值盐密是时间序列全局预测模型；5)利用支持向量机模型求解多变量等值盐密时间序列预测模型根据等值盐密时间序列重构相空间相点构成的支持向量机模型训练样本集容量 及等值盐密非线性系统特性，应用交叉检验法，选择支持向量机模型核函数及模型各参数 为核函数选择高斯核函数；核函数参数Y = 0. 6 ；惩罚c = 50 ；不敏感损失函数参数ε =0. 29 ；部分等值盐密监测点等值盐密时间序列支持向量机预测模型的预测结果如图3 所示，所建立的等值盐密预测模型中，等值盐密预测值的误差基本上控制在12%之内。所述的绝缘子的污闪临界电压预测模型的应用，按如下步骤进行，如图1(c)所 示1)污闪临界电压的确定污闪临界电压值选用50%污闪电压试验获得的数据U5(l%，对等值盐密及气象参数 与污闪临界电压值U5(l%之间的多维非线性关系进行人工神经网络建模，实现对已知等值盐 密及气象参数的绝缘子污闪临界电压队⑽进行预测；其中污闪电压试验获得的数据如下
气压的标准化值湿度的标准化值温度的标准化值ESDD的标准化值U5。％的标准化值0. 39130.87800. 53540.49000. 30700. 86690.87850. 62830.15800. 65800. 43320.75600. 44860.08000. 65880. 61000.69180. 65140.63600. 2343标准化值是由以下公式确定的a, =i为获得的U5。%数据的个数；
max min
a,为某一参数的标准化值，h为该参数的任一值，bmin为该参数所有值中的最小 值，b_为该参数所有值中的最大值；选择一组数据作为检验样本，其它数据作为人工神经 网络的训练样本；2)针对处理后的人工污秽试验数据，建立如附图2的人工神经网络模型；3)在人工神经网络模型的基础上增设动量项并对对BP方法的学习率进行自适应 调整以步骤1)中的标准化数据为训练样本，设定网络参数初始值和网络输出误差允 许值，对网络进行训练，当网络输出误差小于允许值后，保持相应的网络结构参数；4)在人工神经网络模型的基础上引入陡度因子，使权值调整脱离平坦区，得到污 闪临界电压预测模型部分污闪临界电压预测模型的预测结果与实测值比较如图4所示，所建立的污闪 临界电压预测模型中，污闪临界电压值的误差基本上控制在士6%之内；所述的污闪分级预测预警模型的应用，按如下步骤进行 污闪指数预测模型通过公式K =
污闪指数；UFF为污闪临界电压的预测值；
'0uFF>5%\.Q5Uop<Uff20%\.QUOP<UFF<-‘50%0.9WCp < U FF85%0.95UoP < U ff100%0.85UoP < UFFUm为运行电压。
<1.1^
OP
‘OP
决定，其中K为
OP
< 0.98[/,
OP
< 0.95U,
OP
1权利要求
一种区域电网污闪指数预测方法，其特征在于步骤如下步骤1、应用绝缘子表面等值盐密预测模型，实时预测等值盐密当前值；将电网现场采集的环境参数和等值盐密历史值输入绝缘子表面等值盐密预测模型，绝缘子表面等值盐密预测模型输出即为实时预测等值盐密当前值；步骤2、应用绝缘子污闪临界电压预测模型，预测污闪临界电压；将实时预测等值盐密当前值和当前采集环境参数输入绝缘子的污闪临界电压预测模型，绝缘子的污闪临界电压预测模型的输出即为污闪临界电压预测值；步骤3、应用污闪分级预测预警模型，预测电网绝缘子污闪指数；将污闪临界电压预测值输入污闪分级预测预警模型，污闪分级预测预警模型的输出即为预测的电网污闪指数；步骤4、当电网污闪指数为0和5％时，不进行污闪预警；当电网污闪指数为20％时，发布污闪III级预警；当电网污闪指数为50％和85％时，污闪发生概率已经大于50％，发布污闪II级预警；当电网污闪指数为100％时，污闪发生概率已经相当大，该区域电网随时有可能发生污闪，发布污闪I级预警。
2.按权利要求1所述的区域电网污闪指数预测方法，其特征在于所述的绝缘子表面 等值盐密预测模型的应用，按如下步骤进行1)建立多变量等值盐密时间序列；在固定时间间隔对等值盐密数据进行测量，在一系列时刻t1; t2，. . .，tn得到的离散有 序集合Ix1, x2, ... , xj称为离散等值盐密时间序列，简称为等值盐密时间序列；多变量等值盐密时间序列是由等值盐密时间序列及在相同时刻内的气象参数时间序 列构成的多维等值盐密时间序列，是由包括等值盐密时间序列在内的多维等值盐密非线性 动力学系统的表现形式；M维等值盐密时间序列式，&，…，Xn，N代表时刻数，其中Xi = (xiji,x2ji,... ,xMji),即Xl,l Xi,2'" Xl,i'" Xl,N X2t\ 太2’2·.. X2广· X2,NxMiItrt i …xMiN j式中，1 = 1，2，...，队xM,N表示第M个变量在N时刻的数值，xM, i表示第M个变量在i 时刻的数值；2)重构多变量等值盐密时间序列的相空间； 多变量等值盐密时间序列相空间重构的相点为^n = (Xl, , xIiM-T1'-"'xIjn-(Znl-I)T1 ,…,ΧΜ,η，ΧΜ，η-τΜ ,…，ΧΜ,η-(/ηΜ-1)τΜ ) ‘Vi —- (X’,-, xIJ-Ti ,…，Xl,j-(ml-l)ri '—'XM,i' xMj-Tu , — 'XM,i-(mu -1)τΜ ) ^N = (xI1JV，xIN-Ti，·.” Xl,W-(ml-l)r,,…，XM，N , ΧΜ,Ν-τμ，···,ΧΚ ,Ν-[ριΜ- )τΜ )( -表示第M个变量在N时刻在延迟时间为T M嵌入维数为mM的重构相空间中的值；其中n表示第n时刻， 和叫为第i个时间序列的延迟时间和 嵌入维数，重构相空间的嵌入维数m = mi+m2+. . . +mM, M为时间序列的维数；多变量等值盐密时间序列的相空间重构参数延迟时间t的选择采用互信息 法，互信息法是以互信息第一次达到最小时的延时作为相空间重构的延迟时间，由 是等值盐密时间序列时间跨度为(i+1) X的自相关函数，T为相空间重构参数延迟时间；嵌入维数m由 决定，其中 (m+1)是(m+1)维重构的等值盐密系统相空间中的第i个相点，n(i，m)是在m维等 值盐密系统相空间中使相点X&^Oii)是相点XiOii)的最邻近点的整数，| | | |是等值盐密 系统相空间上的欧式距离；3)等值盐密时间序列确定性检验；本发明中采用李雅谱诺夫指数法进行等值盐密时间序列的确定性检验，该指数是相空 间中邻近轨道的平均指数发散率的数值表征，用以刻画混沌运动的初态敏感性，该指数作 为沿轨道长期平均的结果，是一种整体特征，其值总是实数；判断等值盐密时间序列的非线性特性通过计算最大李雅普诺夫指数得到，该方法计算由y(k)曲线的回归直线斜率火为最大指数，其中 表示对重构的等值盐密相空间中每一对最邻近点，计算k个离散时间后的欧式距离，M为时 间序列的维数；4)全局预测多变量等值盐密时间序列；根据泰肯斯延时嵌入定理，只要嵌入维数m和延迟时间T选择合理，重构相空间在 嵌入空间的轨迹就与微分同胚意义下的等值盐密动力学系统等价，且存在光滑映射f: JT。JT，使得:Vi+1 = f ，Vi+1表示重构相空间中第i+1个相点，应用非线性逼近方法构造映射__来近似逼近f，并使_满足 ’最小，其中 表示表示第M个变量在n时刻在延迟时间为T m嵌入维数为mM的重构相空间中的值，、表示第M个变量的延迟时间，mM表示第M个变量的嵌入维 数；5)利用支持向量机模型求解多变量等值盐密时间序列预测模型；通过求解预测模型中的非线性映射，确定等值盐密时间序列预测模型，并使预测模型 在求得的非线性映射7 下的预测误差满足要求，支持向量机理论可以有效地解决等值盐密 数据样本容量偏小的情况下的等值盐密非线性时间序列预测模型中非线性映射 > 的求解问 题，用于逼近等值盐密时间序列预测模型中非线性映射关系的支持向量机方法是支持向量 回归；设等值盐密系统相空间相点构成的样本集为 Xi，Yi) 表示重构相空间中的任一相点，若存在一个超平面g(x) = <wx>+b,w G Rn，b G R，w、b表 示向量参数，为了构造超平面g(x)，使得|yi-g(Xi)| ( e成立，其中，〈 表示向量内积， i = 1，2，…，M，M为等值盐密时间序列的维数，则样本集S = {(Xi，yi)，i = 1，2，. . .，M} 为 ￡ 的近似集，有|<w. x>+b-yi| 彡 e，即一~I M ||2 ^ J , > i = 1，2，...，M; 其中， 为s的点到超平面f(x)的距离di，则有 M，即集合S中的点到超平面的距离最大值为+ ’通过最大化S中的点到超平面距离的上界可得到集合S的最优近似超平面，则最优近似超平面可通过最大化式^ + 得到，因此求解| k| |2的最小化问题即可得到集合S的最优近似超平面，由于等值盐密系统是非线性系统，必须用一个非线性映射P把等值盐密系统相空间中的相点Xi 映射到一个高维空间，然后在高维空间里进行线性回归，由于优化过程中涉及到高维空间 的内积运算，为了避免内积运算，用核函数w (Xi，Xi+1)代替内积〈识{xi+1 )〉来实现等值 盐密系统相空间中非线性回归，此时，等值盐密系统相空间上的支持向量回归问题可转化 为如下的I kl I2优化问题 其中，i = 1，2，...，M，上式为二次规划问题，其Lagrange函数为 其中，、和％'被称为拉格朗日乘子，对任何i = 1，2，...，M，都有等式a,xa/=0， a i >0，a,20成立；在进行等值盐密系统相空间中非线性映射函数逼近时，由于求得的回归函数与实际函 数之间不可避免的存在误差，因此引入松弛变量ξi≥ 0；ξi*0≥,i = 1，2，...，Μ，ξ i 表示松弛变量；此时的优化为 式中c为惩罚参数，且c>0;可得Lagrange对偶问题为 上式作为绝缘子表面等值盐密的预测；求解上式即可得等值盐密系统相空间中非线性映射_ (1)的回归函数g(x)，核函数与71 作为输入量，即可得到下一时刻的等值盐密值，核函数中的参数和支持向量机模型参数中 的惩罚因子是决定支持向量机方法求得的非线性映射7预测性能的最主要因素，采用交叉 验证法进行支持向量机模型参数的选择。
3.按权利要求1所述的区域电网污闪指数预测方法，其特征在于所述的绝缘子的污 闪临界电压预测模型的应用，按如下步骤进行1)污闪临界电压的确定；污闪临界电压值选用50%污闪电压试验获得的数据U5(l%，对等值盐密及气象参数与污 闪临界电压值队⑽之间的多维非线性关系进行人工神经网络建模，实现对已知等值盐密及 气象参数的绝缘子污闪临界电压队⑽进行预测；2)人工神经网络模型的建立；本发明中采用基于BP方法的污闪临界电压预测人工神经网络模型；确定污闪临界电压预测BP人工神经网络模型的输入层、隐层和输出层三层；污闪临界电压预测BP人工神经网络模型包括输入层、隐层和输出层三层结构，根据不 同的人工污秽实验数据，隐层分别为两层或三层节点组成；X1, X2, ...，Xn为输入层节点，包括等值盐密值、温度、湿度、气压、风速和雨量，hi; h2，. . .，hm为隐层节点，ο为输出层节点，是绝缘子污闪临界电压预测值W1，V2,..., Vffl为输 入层至隐层的权值，WnW2,...，WmS隐层至输出层权值，污闪临界电压人工神经网络模型的 训练使用人工污秽耐压试验数据作为训练数据，整个网络的结构和参数由训练进行优化， 网络的输入层维数为n，输出维数为一维即相关预测值的输出，隐层单元维数m在网络训练 学习中优化确定，预测模型的预测值与实际人工污秽耐压试验的结果进行比较，进而对网 络结构和权值进行修正和完善；3)在人工神经网络模型的基础上增设动量项并对对BP方法的学习率进行自适应调整；为提高污闪临界电压预测人工神经网络的训练速度，在权值调整公式中增加一个动量 项，若用W代表污闪临界电压预测人工神经网络中某层的权值矩阵，X代表某层输入向量， 则含有动量项的污闪临界电压预测人工神经网络权值调整向量表达式为Aff(t) = n 6X+a AW(t-l)，公式中字母表示意思a是动量系数，设a G (0,1), n 是神经网络的学习率，AW(t-l)是前一次的权值调整量，AW(t)是本次的权值调整量，动 量项反映了以前的调整经验，对于t时刻的调整起阻尼作用，当误差曲面出现骤然起伏时， 可减小振荡趋势，提高训练速度；在污闪临界电压预测人工神经网络建模中对BP方法的学习率进行自适应调整，学习 率n E (0,1)表示比例系数，设一个初始学习率，若经过一次权值调整后使总误差增加，则 本次调整无效；若经过一次权值调整后使总误差下降，则本次调整有效；4)在人工神经网络模型的基础上引入陡度因子，使权值调整脱离平坦区； 在污闪临界电压预测人工神经网络模型训练过程中引入陡度因子，误差曲面上存在着 平坦区域，权值调整进入平坦区代表是污闪临界电压预测人工神经网络的神经元输出进入 了转移函数的饱和区，如果在进入平坦区后，设法压缩神经元的净输入，使其输出并退出转 移函数的饱和区，就可以改变误差函数的形状，从而使调整脱离平坦区，具体做法是，在原 转移函数中引入一个陡度因子(，使输出为 其中，net为各层节点的输出值，当发现AE接近零而d_o值仍较大时，即认为进入平 坦区，此时令(> 1 ；当退出平坦区后，再令4 =1，当4 > 1时，net坐标压缩了倍，污闪 临界电压预测人工神经网络的神经元的转移函数曲线的敏感段变长，从而使net值退出饱 和值，当(=1时，转移函数恢复原状，对较小的net值具有较高的灵敏度；建立改进的BP人工神经网络预测模型，求解公式0(j) = f(net(j))即可得到绝缘子 的污闪临界电压预测值，其中j = 1，2，. . .，1，1为输出层的节点数，net(j)表示污闪临界 电压预测神经元j的输入总和，f表示一种非线性映射关系。
4.按权利要求1所述的区域电网污闪指数预测方法，其特征在于所述的污闪分级预 测预警模型的应用，按如下步骤进行确定绝缘子的污闪指数，污闪指数预测模型通过公式'0uFF >uuOP5%\.05UoP<UFF<l.lUOP20% 50%\.0UOP 0.9SUo<UFF<l.05UOP Ftt ^ ”决定，其中K为污闪指 P <uFF <l.0UOP数；UFF为污闪临界电压的预测85%0.95UoP <UFF<0.98UOP100%0.85UOio<UFF <0.95UOP4为运行电压，污闪临界电压的预测值为该模型的输入，K为该模型的输出，即为预测的污闪指数。
全文摘要
一种区域电网污闪指数预测方法，属于输配电监测技术领域。方法如下1)应用绝缘子表面等值盐密预测模型，实时预测等值盐密当前值；2)应用绝缘子污闪临界电压预测模型，预测污闪临界电压；3)应用污闪分级预测预警模型，预测电网绝缘子污闪指数；4)判断状态。本发明的优点提出基于相空间重构的多变量等值盐密时间序列预测模型，用支持向量机模型对其求解，解决在等值盐密数据小样本且存在噪声情况下的预测问题，提高了预测精度。
文档编号G06N3/08GK101893674SQ201010224050
公开日2010年11月24日 申请日期2010年7月12日 优先权日2010年7月12日
发明者徐建源, 李永祥, 林莘, 滕云, 苏蔚 申请人:沈阳工业大学