说明书基于回声状态网络的架空输电线路最高温度预测方法
技术领域
本发明属于电力系统保护和控制技术领域，特别涉及一种应用于架空输电线路的基于回声状态网络的线路最高温度预测方法。
背景技术
架空输电线温度是导线焦耳吸热、从外界吸热及向外界散热的共同结果。目前采用的导线的暂态热平衡方程如式(1)所示：
m C p d T c dt = P j + P s - P c - P r - - - ( 1 ) ]]>
式(1)中，m为导线质量(也可由密度表示)，Cp为导线比热容，TC为导线温度，t为时间；Pj为焦耳吸热、Ps为日照吸热；Pc为对流散热、Pr为辐射散热。
为提高求解效率，通常将暂态热平衡公式近似视为常系数一阶微分方程，求得通解得到导线暂态温升近似表达式(2)：
T c ( t ) = T a ( t ) - A B + [ T c ( t 0 ) - T a ( t ) + A B ] e B ( t - t 0 ) m C P - - - ( 2 ) ]]>
式(2)中，t为当前时刻，t0为初始时刻；Tc(t0)为导线初始时刻温度，Ta(t)为t时刻环境温度；参数A包含了焦耳吸热温度不相关分量和日照吸热分量；参数B包含了焦耳吸热温度相关分量、对流散热分量和辐射散热分量。
上述方法的主要缺点是将参数A、参数B视为常数，没有考虑实际运行中二者的时变特性，这会增大温度预测误差。
发明内容
本发明的目的是为解决传统导线暂态温升计算方法精度不高的问题，提出一种基于回声状态网络的架空输电线路最高温度预测方法，本发明应用于架空输电线路的自适应过负荷保护的导线暂态温升预测，可尽可能挖掘线路潜在载流能力，对导线温度进行预估并自适应调整保护动作时间，为调度和紧急控制系统采取措施消除过负荷争取更多时间。
本发明提出的基于回声状态网络的架空输电线路最高温度预测方法，其特征在于， 该方法包括以下步骤：
1)基于回声状态网络对输电沿线最大温差ΔTmax(n)进行预测；
11)采集变电站本地一段时间内各时刻(时间长度由预测精度决定，一般可在1至24小时内取值，采样间隔为分钟级)的导线温度数据，构成本地温度样本Tlocal(1)，Tlocal(2)，……，Tlocal(R1)，R1为样本总数，一般R1的取值不小于200；
12)从输电导线沿线测量点(为电网布局设定)所得温度数据中选出该段时间内每一时刻的最高温度，构建该段时间内沿线最高温度样本Tmax(1)，Tmax(2)，……，Tmax(R1)；
13)用步骤12)中最高温度样本减去步骤11)中对应时刻的本地温度样本，得到该段时间内沿线最大温差样本ΔTmax(1)，ΔTmax(2)，……，ΔTmax(R1)；
14)用步骤13)得到的沿线最大温差样本构建输入样本u(1),u(2),……,u(S1)和输出样本y(1),y(2),……,y(S2)；S1为输入样本总数，S2为输出样本总数，S1＝S2＝R1/2；其中输入样本第i1个数据u(i1)＝[ΔTmax(i1),ΔTmax(i1+1),……,ΔTmax(i1+S1-1)]T(i1为1至S1间任一整数)，输出样本第j1个数据y(j1)＝[ΔTmax(j1+S2)]T(j1为1至S2间任一整数)；
15)用步骤14)的输入样本、输出样本对回声状态网络进行初始化、训练，利用已完成训练的回声状态网络预测未来时刻n的沿线最大温差ΔTmax(n)(该方法属于常规技术方法)；
2)基于回声状态网络对温升变化量A(n)进行预测
21)采集变电站本地一段时间内各时刻(时间长度由预测精度决定，一般可在1至24小时内取值，采样间隔为分钟级)的环境温度数据Ta(k)、电流数据I(k)、日照强度数据Qs(k)；
22)根据采集的Ta(k)、I(k)、Qs(k)及导线参数Rdc20、αj、kj、αs、D，按式(1)计算温升变化量A(k)：
A(k)＝I(k)2Rdc20[αjkj(Ta(k)-20)+1]+αsQs(k)D/mCp (1)
式(1)中，Rdc20为20℃时的直流电阻；αj为20℃时导线的温度电阻系数；kj为集肤效应的修正系数(对含铁质材料的导线取kj＝1，对不含铁质材料的导线，一般取kj＝1.0123)；αs为导线对光照的吸收率，根据导线风化程度取值范围为0.23～0.95(新导线取0.23，运行1个月左右为0.5，一年后约0.9)；D为导线直径，单位为m；Ta(k)是环境温度第k时刻的采样数据，单位为℃；I(k)是电流第k时刻采样数据，单位为A；QS(k)是日照强度第k时刻采样数据，单位为W/m2；k为当前采样时刻编号；
23)用步骤22)中得到的温升变化量A(k)构建参数A的样本A(1)，A(2)，……，A(R2)，R2为样本总数，一般R2的取值不小于200；
24)用步骤23)得到的参数A的样本构建输入样本u(1),u(2),……,u(S3)和输出样本y(1),y(2),……,y(S4)。S3为输入样本总数，S4为输出样本总数，S3＝S4＝R2/2。其中输入样本第i2个数据u(i2)＝[ΔTmax(i2),ΔTmax(i2+1),……,ΔTmax(i2+S3-1)]T(i2为1至S3间任一整数)，输出样本第j2个数据y(j2)＝[ΔTmax(j2+S4)]T(j2为1至S4间任一整数)；
25)用步骤24)的输入样本、输出样本对回声状态网络进行初始化、训练，利用已完成训练的回声状态网络预测未来时刻n的沿线最大温差A(n)(该方法属于常规技术方法)；
3)基于回声状态网络对温升变化量B(n)进行预测
31)采集变电站本地一段时间内各时刻(时间长度由预测精度决定，一般可在1至24小时内取值，采样间隔为分钟级)的导线温度数据Tc(k)；
32)根据步骤31)中的第k时刻的Tc(k)，步骤21)中的Ta(k)、A(k)及导线参数m、Cp，按式(2)计算温升变化量B(k)：
B(k)＝[mCp(Tc(k)-Tc(k-1))/Δt-A(k)]/(Tc(k)-Ta(k)) (2)
式(2)中，m为导线单位质量，即密度，单位为kg/km；Cp是导线比热，单位为J/(kg·K)；Δt为采样间隔，单位为s；Tc(k)是导线温度第k时刻采样数据，单位为℃；Tc(k-1)是导线温度第k-1时刻采样数据，单位为℃；Ta(k)是环境温度第k时刻采样数据，单位为℃；A(k)是温升变化量(在步骤22)根据式(1)计算得到)；k为当前采样时刻编号；k-1是上一个采样时刻编号。
33)用步骤32)中得到的B(k)构建参数B的样本B(1)，B(2)，……，B(R3)，R3为样本总数，一般R3的取值不小于200；
34)用步骤33)得到的参数B的样本构建输入样本u(1),u(2),……,u(S5)和输出样本y(1),y(2),……,y(S6)。S5为输入样本总数，S6为输出样本总数，S5＝S6＝R3/2。其中输入样本第i3个数据u(i3)＝[ΔTmax(i3),ΔTmax(i3+1),……,ΔTmax(i3+S5-1)]T(i3为1至S5间任一整数)，输出样本第j3个数据y(j3)＝[ΔTmax(j3+S6)]T(j3为1至S6间任一整数)；
35)用输入样本、输出样本对回声状态网络进行初始化、训练，利用已完成训练 的回声状态网络预测未来时刻n的沿线最大温差B(n)(该方法属于常规技术方法)；
4)根据步骤2)中得到的A(n)和步骤3)中得到的B(n)，结合当前时刻的采集的环境温度数据Ta，按照式(3)进行迭代，可预测变电站本地导线温度Tlocal(n)；
T local ( n ) = T a - A ( n ) B ( n ) + [ T local ( n ) - T a + A ( n ) B ( n ) ] e B ( n ) m C P - - - ( 3 ) ]]>
5)将步骤1)中预测的沿线最大温差ΔTmax(n)与步骤4)中预测的本地导线温度Tlocal(n)相加，得到预测的沿线最高温度Tmax(n)；
Tmax(n)＝Tlocal(n)+ΔTmax(n) (4)。
本发明的特点及有益效果：
该方法利用变电站内导线温度、环境温度监测数据进行实时温度预测，同时利用长期在线监测温度数据对沿线最大温差进行建模并预测，两部分预测值相加得到输电线路沿线最高温度。该方法解决了工程中采用最高允许温度作为线路载流能力界定时，对线路最大载流能力利用不足的问题。本发明尽可能挖掘了线路潜在载流能力，对导线温度进行预估并自适应调整保护动作时间，该方法具有在线学习能力、预测精度高的特点，可为调度和紧急控制系统采取措施消除过负荷争取更多时间。
附图说明
图1为本发明基于回声状态网络的架空输电线路最高温度预测方法流程框图；
具体实施方式
本发明提出的基于回声状态网络-递归最小二乘法的沿线最大温度预测结合附图及实施例详细说明如下：
本发明提出的基于回声状态网络的架空输电线路最高温度预测方法，如图1所示，
1)对沿线最大温差ΔTmax(n)的预测；
11)采集变电站本地一段时间内各时刻(本实施例的时间长度为8小时内，采样间隔为2分钟)的导线温度数据，构成本地温度样本Tlocal(1)，Tlocal(2)，……，Tlocal(R1)，R1为样本总数，R1＝200；即本实施例本地温度样本为：Tlocal(1)，Tlocal(2)，……，Tlocal(200)；
12)从电网布局的导线沿线测量点所得温度数据中选出该段时间内每一时刻的最高温度，构建沿线最高温度样本Tmax(1)，Tmax(2)，……，Tmax(200)；
13)用步骤12)中最高温度样本减去步骤11)中对应时刻的本地温度样本，得到沿线最大温差样本ΔTmax(1)，ΔTmax(2)，……，ΔTmax(200)；
14)用步骤13)得到的沿线最大温差样本构建输入样本u(1),u(2),……,u(S1)和输出样本y(1),y(2),……,y(S2)。S1为输入样本总数，S2为输出样本总数，S1＝S2＝R1/2＝100。即输入样本为u(1),u(2),……,u(100)；输出样本为y(1),y(2),……,y(100)。其中输入样本第i1个数据u(i1)＝[ΔTmax(i1),ΔTmax(i1+1),……,ΔTmax(i1+99)]T(i1为1至S1间任一整数)，输出样本第j1个数据y(j1)＝[ΔTmax(j1+100)]T(j1为1至S2间任一整数)；
15)用步骤14)的输入样本、输出样本对回声状态网络进行初始化、训练，利用已完成训练的回声状态网络预测未来时刻n的沿线最大温差ΔTmax(n)；具体包括以下步骤：
15-1)对回声状态网络中网络结构进行初始化：设定输入样本为u(1),u(2),……,u(100)，输出样本为y(1),y(2),……,y(100)。根据步骤14)中定义可知，输入样本中任一元素u(k)是维数M＝100的列向量，输出样本中任一元素y(k)是维数L＝1的列向量；
具体步骤以下步骤：
a)获得输入样本u(1),u(2),……,u(t)和输出样本y(1),y(2),……,y(t)，t为样本总数。u(k)代表输入样本中某一元素，y(k)代表输出样本中某一元素。u(k)可具体表示为u(k)＝[u1(k),u2(k),......,uM(k)]T，获取u(k)的维数M；y(k)可具体表示为y(k)＝[y1(k),y2(k),......,yL(k)]T，获取y(k)的维数L。
b)随机生成内部连接权矩阵W0∈RN×N，N一般取值为50至500的整数，此处N＝200。计算W0的谱半径λmax，根据式(6)对W0进行归一化处理得到W∈RN×N；随机生成输入权值矩阵Win∈RN×M和反馈连接权矩阵Wfb∈RN×L。
W＝W0/λmax (11)
c)对输入输出样本信号进行缩放和位移，保证他们来自同一紧集。
d)生成输出权值矩阵Wout∈RL×(M+N+L)，矩阵中各元素为0；
15-2)对回声状态网络进行训练。其中输入样本为u(1),u(2),……,u(100)，输出样本为y(1),y(2),……,y(100)，样本总数t＝100；
具体包括以下步骤：
a)获得输入样本u(1),u(2),……,u(t)和输出样本y(1),y(2),……,y(t)，获得样本 总数t；
b)从输入输出样本中选择最开始t1个数据，一般t1取值为t的一半。即从输入输出样本中选出[u(1),y(1)]、[u(2),y(2)]、…[u(t1),y(t1)]；
c)从A)回声状态网络结构初始化具体步骤中步骤a)获取L，步骤b)获取N。生成y(0)＝[0；0；......；0]L*1，x(0)＝[0；0；......；0]N*1。取内部神经元激活函数f＝[f1,f2,......,fN]，f1～fN通常选择非线性作用函数如tanh或sigmoid函数；
d)令k＝0；
e)向回声状态网络的输入层输入u(k+1)，向输出层输入向量y(k)，利用公式(12)可求得x(k+1)。k＝k+1；
x(k+1)＝f(Wx(k)+Winu(k+1)+Wfby(k)) (12)
f)若步骤e)求得x(t1)，则进入步骤g)；否则继续执行步骤e)。
g)利用u(1)、u(2)……u(t1)，y(1)、y(2)……y(t1)和x(1)、x(2)……x(t1)构造状态收集矩阵S。为避免初始因素的干扰，从t0点开始收集。t0一般取值为t1的10％至20％，此处t0＝t1×20％。即形成矩阵S＝[s(t0),s(t0+1),...s(t1)]∈R(t1-t0+1)×(M+N+L)，其中第k个列向量s(k)＝[u(k)；x(k)；y(k-1)]T。
h)利用y(t0)、y(t0+1)……y(t1)构造期望输出矩阵E：
E＝[(fout)-1y(t0)；(fout)-1y(t0+1)；...(fout)-1y(t1)]∈R(t1-t0+1)×L。
其中fout＝[f1out,f2out,......,fLout]为输出函数，一般取线性恒等函数。
i)根据式(8)计算出(Wout)T，将结果转置得到输出权值矩阵Wout。
(Wout)T＝S-1E (13)
j)从输入输出样本中选择[u(t1+1),y(t1+1)]、[u(t1+2),y(t1+2)]、…[u(t),y(t)]。
k)令k＝t1+1。
l)选取[u(k),y(k)]，根据步骤b)中得到的x(k)和式(9)计算出yforecast(k)。
yforecast(k+1)＝fout(Wout[x(k+1),u(k+1),y(k)]) (14)
m)根据式(14)计算输出信号的预测值与样本值的差值；根据式(15)计算累计平方误差，其中0＜λ＜1，此处λ＝0.9995；根据递归最小二乘法修正Wout。
e(k)＝yforecast(k)-y(k)(14)
E ( k ) = Σ i = 1 k λ k - i e ( k ) 2 - - - ( 15 ) ]]>
n)令k＝k+1。若k>t，则训练完成；否则执行步骤l)；
15-3)完成步骤15-2)中对回声状态网络的训练后，依据式(1)和式(2)进行迭代，预测出样本之后的沿线最大温差序列ΔTmax(201),ΔTmax(202)……，其中ΔTmax(n)＝yforcast(n)；
x(n+1)＝f(Wx(n)+Winu(n+1)+Wfby(n)) (1)
yforecast(n+1)＝fout(Wout[x(n+1),u(n+1),y(n)]) (2)
2)对温升变化量A(n)的预测
21)采集变电站本地一段时间内各时刻(本实施例的时间长度为8小时内，采样间隔为2分钟)的环境温度数据Ta(k)、电流数据I(k)、日照强度数据Qs(k)；
22)根据Ta(k)、I(k)、Qs(k)及导线参数，按式(1)计算温升变化量A(k)：
A(k)＝I(k)2Rdc20[αjkj(Ta(k)-20)+1]+αsQs(k)D/mCp (3)
式(3)中，Rdc20为20℃时的直流电阻；αj为20℃时导线的温度电阻系数；kj为集肤效应的修正系数(对含铁质材料的导线有kj＝1，对不含铁质材料的导线，一般取kj＝1.0123)；αs为导线对光照的吸收率，根据导线风化程度取值范围为0.23～0.95(新导线取0.23，运行1个月左右为0.5，一年后约0.9)；D为导线直径，单位为m；Ta(k)是环境温度第k点采样数据，单位为℃；I(k)是电流第k点采样数据，单位为A；QS(k)是日照强度第k点采样数据，单位为W/m2；k为当前采样点编号；
23)用步骤22)中得到的A(k)构建参数A的样本A(1)，A(2)，……，A(R2)，R2为样本总数，本实施例取值R2＝200；即参数A的样本为A(1)，A(2)，……，A(200)；
24)用步骤23)得到的参数A的样本构建输入样本u(1),u(2),……,u(S3)和输出样本y(1),y(2),……,y(S4)。S3为输入样本总数，S4为输出样本总数，S3＝S4＝R2/2＝100。其中输入样本第i2个数据u(i2)＝[A(i2),A(i2+1),……,A(i2+99)]T(i2为1至S3间任一整数)，输出样本第j2个数据y(j2)＝[A(j2+100)]T(j2为1至S4间任一整数)；
25)用输入样本、输出样本对回声状态网络进行初始化、训练，利用已完成训练 的回声状态网络预测未来时刻n的温升变化量A(n)；参见步骤15)；
3)对温升变化量B(n)的预测
31)采集变电站本地一段时间内各时刻(本实施例的时间长度为8小时内，采样间隔为2分钟)的导线温度数据Tc(k)；
32)根据步骤31)中的Tc(k)，步骤21)中的Ta(k)、A(k)及导线参数，按式(2)计算温升变化量B(k)：
B(k)＝[mCp(Tc(k)-Tc(k-1))/Δt-A(k)]/(Tc(k)-Ta(k)) (4)
式(2)中，m为导线单位质量，即密度，单位为kg/km；Cp是导线比热，单位为J/(kg·K)；Δt为采样间隔，单位为s；Tc(k)是导线温度第k点采样数据，单位为℃；Tc(k-1)是导线温度第k-1点采样数据，单位为℃；Ta(k)是环境温度第k点采样数据，单位为℃；A(k)是第k点温升变化量，在步骤22)根据式(3)计算得到；k为当前采样点编号；k-1是上一个采样点编号；
33)用步骤32)中得到的B(k)构建参数B的样本B(1)，B(2)，……，B(R3)，R3为样本总数，本实施例取值R3＝200；即参数B的样本为B(1)，B(2)，……，B(200)；
34)用步骤33)得到的参数B的样本构建输入样本u(1),u(2),……,u(S5)和输出样本y(1),y(2),……,y(S6)。S5为输入样本总数，S6为输出样本总数，S5＝S6＝R3/2＝100。其中输入样本第i3个数据u(i3)＝[ΔTmax(i3),ΔTmax(i3+1),……,ΔTmax(i3+99)]T(i3为1至S5间任一整数)，输出样本第j3个数据y(j3)＝[ΔTmax(j3+100)]T(j3为1至S6间任一整数)；
35)用输入样本、输出样本对回声状态网络进行初始化、训练，利用已完成训练的回声状态网络预测未来时刻n的温升变化量B(n)；参见步骤15)；
4)根据步骤25)中得到的A(n)和步骤35)中得到的B(n)，结合当前时刻的环境温度数据Ta，按照式(5)进行迭代，可以预测变电站本地导线温度Tlocal(n)；
T local ( n ) = T a - A ( n ) B ( n ) + [ T local ( n ) - T a + A ( n ) B ( n ) ] e B ( n ) m C P - - - ( 5 ) ]]>
5)将步骤17)中预测的沿线最大温差ΔTmax(n)与步骤4)中预测的本地导线温度Tlocal(n)相加，得到预测的沿线最高温度Tmax(n)。
Tmax(n)＝Tlocal(n)+ΔTmax(n) (6)。