CN104598990B - 基于回声状态网络的架空输电线路最高温度预测方法 - Google Patents
基于回声状态网络的架空输电线路最高温度预测方法 Download PDFInfo
- Publication number
- CN104598990B CN104598990B CN201410827951.6A CN201410827951A CN104598990B CN 104598990 B CN104598990 B CN 104598990B CN 201410827951 A CN201410827951 A CN 201410827951A CN 104598990 B CN104598990 B CN 104598990B
- Authority
- CN
- China
- Prior art keywords
- temperature
- max
- sample
- local
- wire
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/06—Electricity, gas or water supply
Abstract
本发明涉及基于回声状态网络的架空输电线路最高温度预测方法,属于电力系统输电线路继电保护领域,该方法包括:基于回声状态网络对输电沿线最大温差ΔTmax(n)进行预测;基于回声状态网络对温升变化量A(n)进行预测;基于回声状态网络对温升变化量B(n)进行预测;根据A(n)和B(n),结合当前时刻的采集的环境温度数据Ta,进行迭代,预测变电站本地导线温度Tlocal(n);将预测的沿线最大温差ΔTmax(n)与预测的本地导线温度Tlocal(n)相加,得到预测的沿线最高温度Tmax(n);本发明尽可能挖掘了线路潜在载流能力,对导线温度进行预估并自适应调整保护动作时间,为调度和紧急控制系统采取措施消除过负荷争取更多时间。
Description
技术领域
本发明属于电力系统保护和控制技术领域,特别涉及一种应用于架空输电线路的基于回声状态网络的线路最高温度预测方法。
背景技术
架空输电线温度是导线焦耳吸热、从外界吸热及向外界散热的共同结果。目前采用的导线的暂态热平衡方程如式(1)所示:
式(1)中,m为导线质量(也可由密度表示),Cp为导线比热容,TC为导线温度,t为时间;Pj为焦耳吸热、Ps为日照吸热;Pc为对流散热、Pr为辐射散热。
为提高求解效率,通常将暂态热平衡公式近似视为常系数一阶微分方程,求得通解得到导线暂态温升近似表达式(2):
式(2)中,t为当前时刻,t0为初始时刻;Tc(t0)为导线初始时刻温度, Ta(t)为t时刻环境温度;参数A包含了焦耳吸热温度不相关分量和日照吸热分量;参数B包含了焦耳吸热温度相关分量、对流散热分量和辐射散热分量。
上述方法的主要缺点是将参数A、参数B视为常数,没有考虑实际运行中二者的时变特性,这会增大温度预测误差。
发明内容
本发明的目的是为解决传统导线暂态温升计算方法精度不高的问题,提出一种基于回声状态网络的架空输电线路最高温度预测方法,本发明应用于架空输电线路的自适应过负荷保护的导线暂态温升预测,可尽可能挖掘线路潜在载流能力,对导线温度进行预估并自适应调整保护动作时间,为调度和紧急控制系统采取措施消除过负荷争取更多时间。
本发明提出的基于回声状态网络的架空输电线路最高温度预测方法,其特征在于,该方法包括以下步骤:
1)基于回声状态网络对输电沿线最大温差ΔTmax(n)进行预测;
11)采集变电站本地一段时间内各时刻的导线温度数据,构成本地导线温度样本Tlocal(1),Tlocal(2),……,Tlocal(R1),R1为本地导线温度的样本总数;
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)计算第k时刻第一温升变化量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为导线直径,单位为米;Ta(k)是环境温度第k时刻的采样数据,单位为℃;I(k)是电流第 k时刻采样数据,单位为安培;QS(k)是日照强度第k时刻采样数据,单位为W/ 平方米;k为当前采样时刻编号;
23)用步骤22)中得到的第k时刻第一温升变化量A(k)构建参数A的样本A(1),A(2),……,A(R2),R2为参数A的样本总数,一般R2的取值不小于200;
24)用步骤23)得到的参数A的样本构建输入样本uA(1),uA(2),……, uA(S3)和输出样本yA(1),yA(2),……,yA(S4);S3为输入样本总数,S4为输出样本总数,S3=S4=R2/2;其中输入样本第i2个数据uA(i2)=[A(i2), A(i2+1),……,A(i2+S3-1)]T,i2为1至S3间任一整数,输出样本第j2个数据 yA(j2)=[A(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)计算第k时刻第二温升变化量B(k):
B(k)=[mCp(Tc(k)-Tc(k-1))/Δt-A(k)]/(Tc(k)-Ta(k)) (2)
式(2)中,m为导线单位质量,即密度,单位为kg/千米;Cp是导线比热,单位为J/(kg·K);Δt为采样间隔,单位为s;Tc(k)是导线温度第k时刻采样数据,单位为℃;Tc(k-1)是导线温度第k-1时刻采样数据,单位为℃;Ta(k)是环境温度第k时刻采样数据,单位为℃;A(k)是第k时刻第一温升变化量(在步骤22) 根据式(1)计算得到);k为当前采样时刻编号;k-1是上一个采样时刻编号。
33)用步骤32)中得到的第k时刻第二温升变化量B(k)构建参数B的样本B(1),B(2),……,B(R3),R3为参数B的样本总数,一般R3的取值不小于200;
34)用步骤33)得到的参数B的样本构建输入样本uB(1),uB(2),……, uB(S5)和输出样本yB(1),yB(2),……,yB(S6);S5为输入样本总数,S6为输出样本总数,S5=S6=R3/2;其中输入样本第i3个数据uB(i3)=[B(i3),B(i3+1),……, B(i3+S5-1)]T,i3为1至S5间任一整数,输出样本第j3个数据yB(j3)=[B(j3+S6)]T, j3为1至S6间任一整数;
35)用输入样本、输出样本对回声状态网络进行初始化、训练,利用已完成训练的回声状态网络预测未来时刻n的第二温升变化量B(n);
4)根据步骤2)中得到的第一温升变化量A(n)和步骤3)中得到的第二温升变化量B(n),结合当前时刻的采集的环境温度数据Ta,按照式(3)进行迭代,可预测变电站本地导线温度Tlocal(n);
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%。即形成矩阵其中第k个列向量s(k)=[u(k);x(k);y(k -1)]T。
h)利用y(t0)、y(t0+1)……y(t1)构造期望输出矩阵E:
其中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)
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)计算第k时刻第一温升变化量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为导线直径,单位为米;Ta(k)是环境温度第k点采样数据,单位为℃;I(k)是电流第k点采样数据,单位为安培;QS(k)是日照强度第k点采样数据,单位为W/平方米; k为当前采样点编号;
23)用步骤22)中得到的第k时刻第一温升变化量A(k)构建参数A的样本A(1),A(2),……,A(R2),R2为参数A的样本总数,本实施例取值R2=200;即参数A的样本为A(1),A(2),……,A(200);
24)用步骤23)得到的参数A的样本构建输入样本uA(1),uA(2),……, uA(S3)和输出样本yA(1),yA(2),……,yA(S4)。S3为输入样本总数,S4为输出样本总数,S3=S4=R2/2=100。其中输入样本第i2个数据uA(i2)=[A(i2), A(i2+1),……,A(i2+S3-1)]T,i2为1至S3间任一整数,输出样本第j2个数据 yA(j2)=[A(j2+S4)]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)及导线参数,按式(4)计算第k时刻第二温升变化量B(k):
B(k)=[mCp(Tc(k)-Tc(k-1))/Δt-A(k)]/(Tc(k)-Ta(k)) (4)
式(4)中,m为导线单位质量,即密度,单位为kg/千米;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)中得到的第k时刻第二温升变化量B(k)构建参数B的样本B(1),B(2),……,B(R3),R3为参数B的样本总数,本实施例取值R3=200;即参数B的样本为B(1),B(2),……,B(200);
34)用步骤33)得到的参数B的样本构建输入样本uB(1),uB(2),……, uB(S5)和输出样本yB(1),yB(2),……,yB(S6);S5为输入样本总数,S6为输出样本总数,S5=S6=R3/2;其中输入样本第i3个数据uB(i3)=[B(i3),B(i3+1),……, B(i3+S5-1)]T,i3为1至S5间任一整数,输出样本第j3个数据yB(j3)=[B(j3+S6)]T, j3为1至S6间任一整数;
35)用输入样本、输出样本对回声状态网络进行初始化、训练,利用已完成训练的回声状态网络预测未来时刻n的第二温升变化量B(n);参见步骤 15);
4)根据步骤25)中得到的第一温升变化量A(n)和步骤35)中得到的第二温升变化量B(n),结合当前时刻的环境温度数据Ta,按照式(5)进行迭代,可以预测变电站本地导线温度Tlocal(n);
5)将步骤17)中预测的沿线最大温差ΔTmax(n)与步骤4)中预测的本地导线温度Tlocal(n)相加,得到预测的沿线最高温度Tmax(n)。
Tmax(n)=Tlocal(n)+△Tmax(n) (6)。
Claims (2)
1.一种基于回声状态网络的架空输电线路最高温度预测方法,其特征在于,该方法包括以下步骤:
1)基于回声状态网络对输电沿线最大温差ΔTmax(n)进行预测;
2)基于回声状态网络对第一温升变化量A(n)进行预测;
3)基于回声状态网络对第二温升变化量B(n)进行预测;
4)根据步骤2)中得到的第一温升变化量A(n)和步骤3)中得到的第二温升变化量B(n),结合当前时刻的采集的环境温度数据Ta,按照式(1)进行迭代,可预测变电站本地导线温度Tlocal(n);
<mrow>
<msub>
<mi>T</mi>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>c</mi>
<mi>a</mi>
<mi>l</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>T</mi>
<mi>a</mi>
</msub>
<mo>-</mo>
<mfrac>
<mrow>
<mi>A</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>B</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>+</mo>
<mo>&lsqb;</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>c</mi>
<mi>a</mi>
<mi>l</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>a</mi>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<mi>A</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>B</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
<msup>
<mi>e</mi>
<mfrac>
<mrow>
<mi>B</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>mC</mi>
<mi>P</mi>
</msub>
</mrow>
</mfrac>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
式(1)中:m为导线单位质量,即密度,单位为kg/千米;Cp是导线比热,单位为J/(kg·K);n为未来某时刻;
5)将步骤1)中预测的沿线最大温差ΔTmax(n)与步骤4)中预测的本地导线温度Tlocal(n)相加,得到预测的沿线最高温度Tmax(n);
Tmax(n)=Tlocal(n)+ΔTmax(n) (2);
所述步骤2)具体包括以下步骤:
21)采集变电站本地一段时间内各时刻的环境温度数据Ta(k)、电流数据I(k)、日照强度数据Qs(k);
22)根据采集的Ta(k)、I(k)、Qs(k)及导线参数Rdc20、αj、kj、αs、D,按式(3)计算第k时刻第一温升变化量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;D为导线直径,单位为米;Ta(k)是环境温度第k时刻的采样数据,单位为℃;I(k)是电流第k时刻采样数据,单位为安培;QS(k)是日照强度第k时刻采样数据,单位为W/平方米;k为当前采样时刻编号;
23)用步骤22)中得到的第k时刻第一温升变化量A(k)构建参数A的样本A(1),A(2),……,A(R2),R2为参数A的样本总数;
24)用步骤23)得到的参数A的样本构建输入样本uA(1),uA(2),……,uA(S3)和输出样本yA(1),yA(2),……,yA(S4);S3为输入样本总数,S4为输出样本总数,S3=S4=R2/2;其中输入样本第i2个数据uA(i2)=[A(i2),A(i2+1),……,A(i2+S3-1)]T,i2为1至S3间任一整数,输出样本第j2个数据yA(j2)=[A(j2+S4)]T,j2为1至S4间任一整数;
25)用步骤24)的输入样本、输出样本对回声状态网络进行初始化、训练,利用已完成训练的回声状态网络预测未来时刻n的第一温升变化量A(n);
所述步骤3)具体包括以下步骤:
31)采集变电站本地一段时间内各时刻的导线温度数据Tc(k);
32)根据步骤31)中的第k时刻的Tc(k),步骤21)中的Ta(k)、A(k)及导线参数m、Cp,按式(4)计算第k时刻第二温升变化量B(k):
B(k)=[mCp(Tc(k)-Tc(k-1))/Δt-A(k)]/(Tc(k)-Ta(k)) (4);
式(4)中,m为导线单位质量,即密度,单位为kg/千米;Cp是导线比热,单位为J/(kg·K);Δt为采样间隔,单位为s;Tc(k)是导线温度第k时刻采样数据,单位为℃;Tc(k-1)是导线温度第k-1时刻采样数据,单位为℃;Ta(k)是环境温度第k时刻采样数据,单位为℃;A(k)是第k时刻第一温升变化量;k为当前采样时刻编号;k-1是上一个采样时刻编号;
33)用步骤32)中得到的第k时刻第二温升变化量B(k)构建参数B的样本B(1),B(2),……,B(R3),R3为参数B的样本总数;
34)用步骤33)得到的参数B的样本构建输入样本uB(1),uB(2),……,uB(S5)和输出样本yB(1),yB(2),……,yB(S6);S5为输入样本总数,S6为输出样本总数,S5=S6=R3/2;其中输入样本第i3个数据uB(i3)=[B(i3),B(i3+1),……,B(i3+S5-1)]T,i3为1至S5间任一整数,输出样本第j3个数据yB(j3)=[B(j3+S6)]T,j3为1至S6间任一整数;
35)用输入样本、输出样本对回声状态网络进行初始化、训练,利用已完成训练的回声状态网络预测未来时刻n的第二温升变化量B(n)。
2.如权利要求1所述的方法,其特征在于,所述步骤1)具体包括以下步骤
11)采集变电站本地一段时间内各时刻的导线温度数据,构成本地导线温度样本Tlocal(1),Tlocal(2),……,Tlocal(R1),R1为本地导线温度的样本总数;
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)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410827951.6A CN104598990B (zh) | 2014-12-25 | 2014-12-25 | 基于回声状态网络的架空输电线路最高温度预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410827951.6A CN104598990B (zh) | 2014-12-25 | 2014-12-25 | 基于回声状态网络的架空输电线路最高温度预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104598990A CN104598990A (zh) | 2015-05-06 |
CN104598990B true CN104598990B (zh) | 2017-12-26 |
Family
ID=53124758
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410827951.6A Active CN104598990B (zh) | 2014-12-25 | 2014-12-25 | 基于回声状态网络的架空输电线路最高温度预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104598990B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113161983B (zh) * | 2021-03-29 | 2022-06-03 | 重庆大学 | 一种计及动态热特性的输电线路自适应过负荷保护方法 |
CN113241762B (zh) * | 2021-05-31 | 2023-09-26 | 东北大学 | 基于事件触发的回声状态网络自适应负荷频率控制方法 |
CN115642567B (zh) * | 2022-11-07 | 2024-04-05 | 黑龙江科技大学 | 一种电路电器保护方法、装置、设备及存储介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101614602A (zh) * | 2009-08-03 | 2009-12-30 | 电子科技大学 | 输电线路监测方法及装置 |
CN101672699A (zh) * | 2009-09-21 | 2010-03-17 | 华南理工大学 | 电缆导体暂态温度确定方法及确定装置 |
CN102830314A (zh) * | 2012-08-30 | 2012-12-19 | 广东电网公司佛山供电局 | 架空导线载流量检测方法 |
CN102982393A (zh) * | 2012-11-09 | 2013-03-20 | 山东电力集团公司聊城供电公司 | 一种输电线路动态容量的在线预测方法 |
CN103336908A (zh) * | 2013-07-23 | 2013-10-02 | 哈尔滨工业大学 | 单调回声状态网络的剩余寿命预测方法 |
CN103808426A (zh) * | 2014-03-07 | 2014-05-21 | 中国科学院半导体研究所 | 电缆线芯温度的间接测量方法 |
CN103886374A (zh) * | 2014-04-22 | 2014-06-25 | 武汉大学 | 一种基于rbf神经网络的电缆接头导线温度预测方法 |
-
2014
- 2014-12-25 CN CN201410827951.6A patent/CN104598990B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101614602A (zh) * | 2009-08-03 | 2009-12-30 | 电子科技大学 | 输电线路监测方法及装置 |
CN101672699A (zh) * | 2009-09-21 | 2010-03-17 | 华南理工大学 | 电缆导体暂态温度确定方法及确定装置 |
CN102830314A (zh) * | 2012-08-30 | 2012-12-19 | 广东电网公司佛山供电局 | 架空导线载流量检测方法 |
CN102982393A (zh) * | 2012-11-09 | 2013-03-20 | 山东电力集团公司聊城供电公司 | 一种输电线路动态容量的在线预测方法 |
CN103336908A (zh) * | 2013-07-23 | 2013-10-02 | 哈尔滨工业大学 | 单调回声状态网络的剩余寿命预测方法 |
CN103808426A (zh) * | 2014-03-07 | 2014-05-21 | 中国科学院半导体研究所 | 电缆线芯温度的间接测量方法 |
CN103886374A (zh) * | 2014-04-22 | 2014-06-25 | 武汉大学 | 一种基于rbf神经网络的电缆接头导线温度预测方法 |
Non-Patent Citations (2)
Title |
---|
"基于遗传算法优化BP神经网络的电缆接头导线温度预测";周辉 等;《电气应用》;20131231(第S2期);第581-585页 * |
"架空输电线路输电能力的研究";梁立凯;《中国博士学位论文全文数据库工程科技Ⅱ辑》;20131015(第10期);第C042-57页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104598990A (zh) | 2015-05-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Premalatha et al. | Prediction of solar radiation for solar systems by using ANN models with different back propagation algorithms | |
Hasan et al. | An improved dynamic thermal current rating model for PMU-based wide area measurement framework for reliability analysis utilizing sensor cloud system | |
Ishaque et al. | The performance of perturb and observe and incremental conductance maximum power point tracking method under dynamic weather conditions | |
Mellit et al. | Performance prediction of 20 kWp grid-connected photovoltaic plant at Trieste (Italy) using artificial neural network | |
CN110321601B (zh) | 一种架空线路动态载流能力超前预测方法及系统 | |
Voyant et al. | Twenty four hours ahead global irradiation forecasting using multi‐layer perceptron | |
JP6104235B2 (ja) | 建物のエネルギー消費の予測方法 | |
CN104598990B (zh) | 基于回声状态网络的架空输电线路最高温度预测方法 | |
Rao et al. | Estimation of daily global solar radiation using temperature, relative humidity and seasons with ANN for Indian stations | |
CN109598061A (zh) | 一种变压器群平均寿命损失的监测方法 | |
CN102982393A (zh) | 一种输电线路动态容量的在线预测方法 | |
CN105678439A (zh) | 基于bp神经网络的输电线路动态增容运行风险评估方法 | |
CN107110995A (zh) | 日射量预测方法 | |
CN112152206B (zh) | 一种输电线路动态增容方法及系统 | |
Mohanty | ANFIS based prediction of monthly average global solar radiation over Bhubaneswar (State of Odisha) | |
Yadav et al. | ANN based prediction of daily global solar radiation for photovoltaics applications | |
Nomiyama et al. | A study on global solar radiation forecasting using weather forecast data | |
JP5735283B2 (ja) | 使用電力量予測システム、及び使用電力量予測システムの制御方法 | |
Kumar et al. | New regression model to estimate global solar radiation using artificial neural network | |
CN107194507A (zh) | 一种基于组合支持向量机的风电场短期风速预测方法 | |
CN104253586B (zh) | 一种太阳能电池板电气参数在线测量评价装置及方法 | |
Prusty et al. | A detailed formulation of sensitivity matrices for probabilistic load flow assessment considering electro-thermal coupling effect | |
Yang et al. | Overhead conductor thermal dynamics identification by using echo state networks | |
Androvitsaneas et al. | Artificial neural network methodology for the estimation of ground enhancing compounds resistance | |
Mohanty et al. | Comparison and prediction of monthly average solar radiation data using soft computing approach for Eastern India |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant |