CN109492823A - 一种对输电线路覆冰厚度的预测方法 - Google Patents
一种对输电线路覆冰厚度的预测方法 Download PDFInfo
- Publication number
- CN109492823A CN109492823A CN201811416279.6A CN201811416279A CN109492823A CN 109492823 A CN109492823 A CN 109492823A CN 201811416279 A CN201811416279 A CN 201811416279A CN 109492823 A CN109492823 A CN 109492823A
- Authority
- CN
- China
- Prior art keywords
- information
- ice
- covering thickness
- att
- weather forecast
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- 230000005540 biological transmission Effects 0.000 claims abstract description 39
- 230000005611 electricity Effects 0.000 claims abstract description 24
- 238000012360 testing method Methods 0.000 claims abstract description 23
- 238000012549 training Methods 0.000 claims abstract description 15
- 230000000694 effects Effects 0.000 claims abstract description 13
- 230000007246 mechanism Effects 0.000 claims description 17
- 239000013598 vector Substances 0.000 claims description 15
- 230000008569 process Effects 0.000 claims description 10
- 239000000203 mixture Substances 0.000 claims description 5
- 238000013528 artificial neural network Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000011478 gradient descent method Methods 0.000 claims description 3
- 230000002123 temporal effect Effects 0.000 claims description 3
- 238000011160 research Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000001556 precipitation Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000006403 short-term memory Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000015654 memory Effects 0.000 description 1
- 238000012806 monitoring device Methods 0.000 description 1
- 238000000465 moulding Methods 0.000 description 1
- 210000004218 nerve net Anatomy 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
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—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种对输电线路覆冰厚度的预测方法,包括:步骤1:采集输电线路的微气象、覆冰厚度信息和天气预报信息;步骤2:将采集到的微气象信息、覆冰厚度信息和天气预报信息划分为训练数据集和测试数据集;步骤3:训练得到输电线路覆冰厚度预测模型Ice_Thick_Predict_Model;步骤4:使用测试集测试模型Ice_Thick_Predict_Model的效果,根据测试效果调整模型的超参数,然后重复跳转到步骤3,直至获得理想的模型。步骤5:实现短期输电线路覆冰厚度的预测。
Description
技术领域
本发明属于机器学习技术领域,尤其涉及一种对输电线路覆冰厚度的预测方法。
背景技术
输电线路冰风灾害一直威胁着供电系统的安全。世界各地都曾出现过不同程度的输电线路冰风灾害,我国电网系统也发生过多次冰风灾害,而且这种灾害的发生频率呈现上升趋势。开展电力系统防灾研究、保障电力系统在冰风灾害时的安全稳定运行是当前刻不容缓的挑战和艰巨的任务。
在电力线路覆冰监测方面,国外做了较多的理论研究和应用实践工作。如加拿大、美国、俄罗斯等国家很早就将远程监测系统用于架空输电线路覆冰的实时监测中。早在20世纪70年代,加拿大在跨越魁北克无人区的输电工程中就使用了可自动记录输电线路覆冰的测冰仪,这种观测仪器可以实现数据的远程传输,而在此后的40多年间,这套仪器对输电线路覆冰进行了连续观测,记录了魁北克地区冬季大量的架空线路覆冰数据,为魁北克地区绘制冻雨频率图提供了重要依据。在20世纪80年代早期,美国陆军寒区工程与研究实验室在华盛顿地区建设了一条安装有各种监测设备的架空输电线路,用以监测各种微气象变化下的线路覆冰情况。2005年底,乌克兰国家电网在架空线路上安装了8套覆冰监控系统,用于其国内覆冰区架空线路覆冰状况的实时监测。俄罗斯开发的冰风负荷传感器主要用于导线覆冰的早期发现,该设备主要架设在中部伏尔加格勒市的架空线路上运行。
近年来,尽管我国针对电网输电线路冰风灾害应急管理的研究取得了较快的进展,但由于研究起步晚,与冰风自然灾害带来的严峻挑战相比,输电线路冰风灾害应急管理研究的内容和深度都不够。因此对于输电线路冰风灾害预警管理有待进一步深入细致的进行专题研究。这是因为输电线路冰风灾害不能避免的情况,只有密切监视和研究风险的变化规律,才能有效地进行输电线路冰风灾害风险预警,在灾害的萌芽及演化阶段进行有效应对。有必要将输电线路冰风灾害防御研究重心从事中或事后的应急管理向事前风险预警管理转变,将电网冰风灾害应急管理框架向冰风自然灾害预警领域延伸。只有从灾害源头和灾害演变过程中开展多层次预警管理,才能有效地将预警的时间尺度进行扩展,从而延长留给决策者的应对时间,变被动应急为主动预防,将输电线路冰风灾害损失降低到最小。
我国在在输电线路冰风灾害特征识别和预警管理方面,一直缺乏比较系统的研究。目前已有的输电线路覆冰厚度预测模型主要是对气象数据和覆冰厚度之间的关系建立模型,然后根据建立的模型对覆冰厚度进行预测。该类模型一般很难拟合气象信息与覆冰厚度之间时间序列关系,无法准确对覆冰厚度进行预测。
发明内容
本发明提出一种基于注意力机制和LSTM编解码(Encoder-Decoder)模型的输电线路覆冰厚度预测方法。该方法主要由LSTM编码模块、时序特征权重提取模块(注意力机制)、LSTM解码模块三部分组成,另外该方法在预测时还结合了天气预报信息,实现对短期输电线路覆冰厚度进行预测。
本方法主要有两点优势,一、结合注意力机制和LSTM模型,解决了传统覆冰预测方法无法学习长序列气象信息的问题;二、结合微气象信息和天气预报信息对输电线路覆冰厚度进行预测,弥补了传统预测方法缺失未来天气信息问题,提高了输电线路覆冰厚度预测精度。
为解决上述技术问题,本发明采用的技术方案如下:
包括如下步骤:
步骤1:采集输电线路的微气象、覆冰厚度和天气预报信息;
步骤2:统计微气象、覆冰厚度信息和天气预报信息组成的数据条目数,并按照7:3的比例将微气象、覆冰厚度信息和天气预报信息条目划分为训练集和测试集;
步骤3:使用训练数据拟合预测模型Ice_Thick_Predict_Model:结合LSTM(Longshort-term memory:长短时记忆)模型与Attention注意力机制将当前n个微气象信息、气象预报信息、输电线路覆冰厚度信息、未来m个时间序列的天气预报信息与未来m个时间序列的输电线路覆冰厚度信息进行拟合,得到输电线路覆冰厚度预测模型Ice_Thick_Predict_Model;
步骤4:使用测试集测试模型Ice_Thick_Predict_Model的效果,根据测试效果调整模型的超参数,直至获得理想的模型。
步骤5:将最近(比如最近三天)的微气象信息序列、输电线路覆冰厚度信息序列与天气预报信息序列作为模型的输入序列,并将该序列输入到输电线路覆冰厚度预测模型Ice_Thick_Predict_Model,实现短期内输电线路覆冰厚度的预测。
步骤1中,每隔一小时采集一次微气象信息和覆冰厚度信息;微气象信息包括:psur_obs,t2m_obs,q2m_obs,w10m_obs,d10m_obs,rh2m_obs,u10m_obs,v10m_obs,RAIN_obs,ice_obs(分别表示:地面气压、地面以上2米高度处温度、地面以上2米高度处比湿、地面以上2米高度处相对湿度、地面以上10米高度处风速、地面以上10米高度处风向、地面以上10米高度处经向风、地面以上10米高度处纬向风、地面1小时累计降水量、观测点输电线路覆冰厚度)共10种信息;在输电线路附近架设观察导线,将观察导线上覆冰厚度测量值作为输电线路覆冰厚度值,观察导线需要和输电线路具有相同的规格,覆冰厚度由人工测量;天气预报信息从气象部门获取。
步骤2包括:
步骤2-1,将采集到的微气象、覆冰厚度信息和天气预报信息按照时间信息建立一一对应关系;
步骤2-2,将整理好的微气象信息、覆冰厚度信息和天气预报信息按照时间序列排序,时间间隔是一个小时,其中,前n(n是超参数)个微气象信息、覆冰厚度信息和天气预报信息序列作为LSTM编码器模型的输入;接下来的m(m是超参数)个天气预报信息作为Attention模块的输入信息,m(m是超参数)表示预测的覆冰厚度序列数量。
步骤2-3,统计微气象、覆冰厚度信息和天气预报信息组成的数据条目数,并按照7:3的比例将微气象、覆冰厚度信息和天气预报信息条目划分为训练集和测试集。
步骤3包括:
步骤3-1,使用LSTM编码器(Encoder)模型对长度序列为n(n为超参数,可以取24~96)的微气象信息、覆冰厚度信息和天气预报信息序列进行时间序列编码,编码结果为C1,C2,C3,......,Ct,......,Cn,Ct表示第t个微气象、覆冰厚度和天气预报序列数据的编码结果。
步骤3-2,将未来m个时间点的天气预报信息W′1,W′2,W′3,......,W′t,......,W′m和LSTM编码器编码结果C1,C2,C3,......,Ct,......,Cn作为输入,利用注意力机制计算LSTM解码器所需的输入序列Att1,Att2,Att3,......,Attt,......,Attm,其中,Attt是注意力机制模块对LSTM编码器(Encoder)模型编码结果的加权和;
步骤3-3,将Att1,Att2,Att3,......,Attt,......,Attm和隐变量序列H′0,H′1,H′2,......,H′t-1,......,H′m-1依次输入至LSTM解码器,利用LSTM解码器(Decoder)预测覆冰厚度信息
步骤3-4,计算预测覆冰厚度与实际覆冰厚度的均方误差MSE(Θ);
步骤3-5,反向传播,使用随机梯度下降法对模型参数Θ进行优化;
步骤3-6,读取下一组序列化的微气象信息、覆冰厚度信息和天气预报信息,跳回到步骤3-1,直至将训练集中的数据读取完。
步骤3-7,重复执行步骤3-1至3-6并持续对模型的参数进行优化,直至使MSE(Θ)下降到合理的数量级,建议使MSE(Θ)下降到1以下。
步骤3-1包括:
通过LSTM编码器(Encoder)对序列化的微气象信息Xt进行编码,输出序列化的隐变量序列Ht和LSTM编码器编码结果Ct,计算公式如(1)所示:
Ht,Ct=f(Ht-1,Xt)(1)
其中,Ht表示LSTM编码器输出的第t个隐变量;Ct是LSTM编码器输出的第t个输出结果;Xt是三层全连接神经网络对第t个微气象、覆冰厚度信息和天气预报信息的编码值;H0是随机生成的隐变量,H0与X1作为f(·)函数首次的输入信息;f(·)是LSTM编码器(Encoder)模块。
步骤3-2包括:
输入LSTM编码器编码结果C1,C2,C3,......,Ct,......,Cn和未来m个时间点的天气预报信息W′1,W′2,W′3,......,W′t,......,W′m,未来第i(取值范围[1,m])个时间点上的天气预报信息为W′i,根据输入信息,计算LSTM编码器编码结果C1,C2,C3,......,Ct,......,Cn与天气预报信息W′1,W′2,W′3,......,W′t,......,W′m之间的相关系数Ai,j,计算过程如公式(2)(3)所示:
Ei,j=tanh(w(W′i)+u(Cj)(2)
其中,tanh(·)是双曲正切函数,w(·)和u(·)分别表示一个三层全连接神经网络(这是两个单独的神经网络,他们的参数是自动学习的,最终两个网络会学习到不一样的参数,来拟合预测模型);W′i是未来第i(取值范围[1,m])个时间点上的天气预报信息;Cj是历史序列上第j个微气象、覆冰厚度信息和天气预报信息序列数据使用LSTM编码器(Encoder)的编码结果;Ei,j是利用注意力机制计算输入序列权重的中间变量,是公式(3)计算权重的依据;
exp(·)是指数函数;n代表预测未来天气所用到的历史微气象信息、覆冰厚度信息和天气预报信息的序列长度;i代表需要预测的覆冰厚度信息序列编号(取值范围[1,m]);j,k代表历史气象信息序列编号(取值范围[1,n])。Ei,j和Ei,k是公式(2)输出的中间变量;
根据权重系数Ai,j和LSTM编码器编码结果计算LSTM解码器输入序列Att1,Att2,Att3,......,Attt,......,Attm:
其中,Cj是第j个微气象、覆冰厚度信息和天气预报信息序列使用LSTM编码器(Encoder)的编码结果;Ai,j是公式(3)输出的权重系数;Atti是最终输出序列化的特征向量,这些特征向量Att1,Att2,Att3,......,Attt,......,Attm作为LSTM解码器(Decoder)的输入信息;Attt表示第t个特征向量,j代表历史气象信息序列编号(取值范围[1,n]);n代表历史气象信息序列的总条目数;i代表需要预测的序列编号(取值范围[1,m]),其中m是需要预测的覆冰厚度总数;
步骤3-3包括:
输入序列化特征向量Attt和隐变量H′t-1,利用LSTM解码器(Decoder)模块预测特征向量Attt对应的输电线路覆冰厚度信息,输出对应的覆冰厚度,计算过程如公式(5)所示:
其中,H′t表示LSTM编码器输出的第t个隐变量;是预测的未来第t个时间点的覆冰厚度;Attt是公式(4)求得的LSTM解码器输入信息;H′0是随机生成的隐变量,H′0和Att1作为h(·)函数首次的输入信息;h(·)是LSTM解码器(Decoder)模块。
步骤3-4包括:通过如下公式计算预测覆冰厚度与实际覆冰厚度的均方误差MSE(Θ):
其中,yi代表覆冰厚度的真实值;代表覆冰厚度的预测值;m代表LSTM解码器(Decoder)模块预测的覆冰厚度个数。通过最小化均方误差,使得预测值与真实值良好拟合。
步骤4包括:将测试集在模型Ice_Thick_Predict_Model上测试,检测模型的预测效果,并根据模型的预测效果调整超参数信息,超参数信息包括:学习率和输入微气象信息、覆冰厚度信息和天气预报信息序列长度n等内容。最终,使模型能在测试集上取得满意的表现,即预测的覆冰厚度与真实值之间均方误差在一个合理的范围内,否则跳到步骤3,继续调整超参数信息,直到得到理想的输电线路覆冰厚度预测模型Ice_Thick_Predict_Model。
与现有技术相比,本发明具有如下优点:
本发明是通过LSTM模型建立微气象、历史覆冰厚度和天气预报信息与输电线路未来覆冰厚度值在时序上的关系模型,结合注意力机制和未来天气预报信息学习不同时间点的微气象、历史覆冰厚度和天气预报信息对未来覆冰厚度的影响程度权重,最终本模型可以结合未来天气预报信息对输电线路覆冰厚度进行预测。
(1)本发明使用LSTM模型提取微气象、历史覆冰厚度和天气预报信息与未来覆冰厚度在时间序列上的关系模型。LSTM可以对任意长度的时间序列建模,可以对时序关系进行良好的拟合,提高了覆冰厚度预测的准确性。
(2)本发明使用注意力(Attention)机制,可以学习到不同时间节点微气象信息对输电线路覆冰厚度的影响权重,提高了输电线路覆冰厚度预测的准确性。
(3)本发明在预测覆冰厚度时结合了历史的微气象信息和天气预报信息,提高了覆冰预测的准确性。
附图说明
下面结合附图和具体实施方式对本发明做更进一步的具体说明,本发明的上述或其他方面的优点将会变得更加清楚。
图1是本发明框架图。
图2使用LSTM模型编码微气象、历史覆冰厚度和天气预报信息序列数据示意图。
图3结合注意力机制和未来天气预报信息学习不同时间点的微气象、历史覆冰厚度和天气预报信息对未来覆冰厚度的影响程度权重的过程。
图4 LSTM解码流程图。
具体实施方式
下面结合附图及实施例对本发明做进一步说明。
本发明具体包括:
步骤1:采集输电线路的微气象和覆冰厚度信息,记录输电线路附近的天气预报信息。
在覆冰高发区域的输电线路附近安装微气象信息采集装置、输电线路覆冰信息实时监测系统。该系统可以实时获取输电线路的微气象信息和输电线路的覆冰信息。采集数据的周期和内容如下所示。
数据采集周期:一小时;
数据采集内容:微气象信息、输电线路覆冰厚度信息和天气预报信息。微气象信息包括:psur_obs,t2m_obs,q2m_obs,w10m_obs,d10m_obs,rh2m_obs,u10m_obs,v10m_obs,RAIN_obs,ice_obs(分别表示:地面气压、地面以上2米高度处温度、地面以上2米高度处比湿、地面以上2米高度处相对湿度、地面以上10米高度处风速、地面以上10米高度处风向、地面以上10米高度处经向风、地面以上10米高度处纬向风、地面1小时累计降水量、观测点输电线路覆冰厚度)共10种信息。天气预报信息从气象部门收集。
步骤2:将采集到的微气象、覆冰厚度信息和天气预报信息划分为训练数据集和测试数据集。
首先,将采集到的微气象、覆冰厚度信息和天气预报信息按照时间建立一一对应关系。
其次,将整理好的微气象、覆冰厚度信息和天气预报信息按照时间序列排序。
最后,按照7:3的比例将数据集划分为训练集和测试集。
步骤3:使用训练数据拟合预测模型Ice_Thick_Predict_Model。
结合LSTM编解码(Encoder-Decoder)模型与注意力(Attention)机制对微气象信息与输电线路覆冰厚度信息进行拟合,得到输电线路覆冰厚度预测模型(Ice_Thick_Predict_Model)。
本发明使用的LSTM模型采用了编解码(Encoder-Decoder)的结构,该结构首先将序列化的微气象、历史覆冰厚度和天气预报信息输入编码器(Encoder)得到一个输出序列;然后利用注意力(Attention)机制模块,从编码器(Encoder)输出序列中提取解码器(Decoder)所需的输入向量序列;解码器(Decoder)根据注意力(Attention)机制模块输出的向量序列预测输电线路覆冰厚度。模型的具体训练步骤如下:
步骤3.1:使用LSTM模型对微气象序列数据编码。
输入:序列化的微气象信息、历史覆冰厚度信息和天气预报信息拼接成的序列:X1,X2,X3,......,Xt,......,Xn和隐变量序列H0,H1,H2,......,Ht-1,......,Hn-1;
处理:LSTM编码器(Encoder)对序列化的微气象等信息进行处理;
输出:序列化的隐变量序列H1,H2,H3,......,Ht,......,Hn和LSTM编码器编码结果C1,C2,C3,......,Ct,......,Cn。
对微气象信息编码的过程如图2所示。其中,Xt代表输入;LSTM编码器A代表LSTM编码器(Encoder)模块;Ht-1是输入隐变量,Ht是输出隐变量;Ct是LSTM编码器编码结果;H0是一个随机生成的隐变量,n代表预测未来覆冰厚度所需的历史信息个数。
隐变量H1,H2,H3,......,Ht,......,Hn、LSTM编码器编码结果C1,C2,C3,......,Ct,......,Cn与序列化的微气象信息X1,X2,X3,......,Xt,......,Xn的关系如公式(1)所示。其中,f(·)代表长短时记忆网络(LSTM);Ht和Ct分别代表输出的隐变量和输出编码值;Xt代表输入序列中第t个微气象、历史覆冰厚度和天气预报信息:
Ht,Ct=f(Ht-1,Xt)(1)
步骤3.2:利用注意力机制计算LSTM解码器(Decoder)模块输入序列Att1,Att2,Att3,......,Attt,......,Attm。
输入:步骤3.1输出LSTM编码器编码结果C1,C2,C3,......,Ct,......,Cn和第i个预测时点的天气预报信息W′i。
处理(如图3所示):首先,根据输入信息,计算LSTM编码器编码结果C1,C2,C3,......,Ct,......,Cn与第i个预测时点的天气预报信息W′i之间的相关系数Ai,j,计算方法如公式(2)(3)所示。其次,根据相关系数Ai,j计算LSTM解码器(Decoder)输入序列Att1,Att2,Att3,......,Attt,......,Attm,计算过程如公式(4)所示。
Ei,j=tanh(w(W′i)+u(Cj))(2)
输出:序列化的特征向量Att1,Att2,Att3,......,Attt,......,Attm,这些特征作为LSTM解码器(Decoder)的输入信息。
步骤3.3:利用LSTM解码器(Decoder)预测覆冰厚度信息,具体过程如图4所示。最终输出结果为预测的覆冰厚度信息。
输入:隐变量H′t-1和步骤3.2输出的序列化特征向量Attt,本项目使用H′t-1作为LSTM解码器(Decoder)解码时第t步的隐变量;H′0是随机生成的隐变量。
处理:利用LSTM解码器(Decoder)预测特征向量Attt和隐变量H′t-1对应的输电线路覆冰厚度信息
输出:与序列化特征向量Attt和隐变量H′t-1对应的覆冰厚度序列与H′t。计算过程如公式(5)所示。
步骤3.4:计算预测覆冰厚度与实际覆冰厚度的均方误差(MSE(Θ)),计算过程如公式(6)所示。
步骤3.5:反向传播,使用随机梯度下降法对模型参数Θ进行优化。
步骤3.6:读取下一组序列化的微气象信息,重复执行步骤3.1至3.5,直至将训练集数据读取完成。
步骤3.7:重复执行步骤3.1至3.6并对模型的参数进行优化,直至使MSE(Θ)下降到合理的数量级,建议使MSE(Θ)下降到1以下。
步骤4:使用测试集测试Ice_Thick_Predict_Model的效果,根据测试效果调整训练模型的超参数,然后跳转到步骤3,直至获得理想的模型。
将测试集部署在模型Ice_Thick_Predict_Model上进行测试,检测模型的效果,并根据模型的效果调整超参数信息,如:学习率,输入序列长度等超参数,然后跳转到步骤3重新训练,直到模型能在测试集上取得满意的效果,通过本步骤可以得到最终的输电线路覆冰厚度预测模型Ice_Thick_Predict_Model。
步骤5:将最近(比如最近三天)的微气象信息、输电线路覆冰厚度信息与天气预报信息序列组成模型的输入序列,并将该序列输入到输电线路覆冰厚度预测模型Ice_Thick_Predict_Model,实现短期内输电线路覆冰厚度的预测。如图1所示其中X1,X2,X3,X4,X5代表微气象信息、历史覆冰厚度信息和天气预报信息拼接成的序列;Att1,Att2,Att3,Att4,Att5是使用注意力机制提取到的中间结果;代表该方法最终预测的输电线路覆冰厚度。具体的预测结果的形式如下:
1月1日0点覆冰厚度10mm,
1月1日1点覆冰厚度12mm,
1月1日2点覆冰厚度15mm,
1月1日3点覆冰厚度15mm,
1月1日4点覆冰厚度14mm,
......。
得到该预测结果后,可以由工程维护单位根据覆冰厚度判断有没有倒塔或者断线的风险。
本发明提供了一种对输电线路覆冰厚度的预测方法,具体实现该技术方案的方法和途径很多,以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有技术加以实现。
Claims (6)
1.一种对输电线路覆冰厚度的预测方法,其特征在于,包括如下步骤:
步骤1:采集输电线路的微气象、覆冰厚度信息和天气预报信息;
步骤2:将采集到的微气象、覆冰厚度信息和天气预报信息划分为训练数据集和测试数据集;
步骤3:使用训练数据拟合预测模型Ice_Thick_Predict_Model:结合LSTM模型与Attention注意力机制对微气象、气象预报信息和输电线路覆冰厚度信息进行拟合,得到输电线路覆冰厚度预测模型Ice_Thick_Predict_Model;
步骤4:使用测试集测试模型Ice_Thick_Predict_Model的效果,根据测试效果调整模型的超参数,然后重复跳转到步骤3,直至获得理想的模型;
步骤5:将最近的微气象信息序列、输电线路覆冰厚度信息序列与天气预报信息序列作为模型的输入序列,并将该序列输入到输电线路覆冰厚度预测模型Ice_Thick_Predict_Model,实现短期输电线路覆冰厚度预测。
2.根据权利要求1所述的方法,其特征在于,步骤2包括:
步骤2-1,将采集到的微气象、覆冰厚度信息和天气预报信息按照时间信息建立一一对应关系;
步骤2-2,将整理好的微气象、覆冰厚度信息和天气预报信息按照时间序列排序,时间间隔是一个小时;
步骤2-3,统计微气象、覆冰厚度信息和天气预报信息组成的数据条目数,并按照7:3的比例将微气象、覆冰厚度信息和天气预报信息条目划分为训练集和测试集。
3.根据权利要求2所述的方法,其特征在于,步骤3包括:
步骤3-1,使用LSTM编码器模型对长度序列为n的微气象信息、覆冰厚度信息和天气预报信息序列进行时间序列编码,编码结果为C1,C2,C3,......,Ct,......,Cn,Ct表示第t个微气象、覆冰厚度信息和天气预报信息序列数据的编码结果;
步骤3-2,将未来m个时间点的天气预报信息W′1,W′2,W′3,......,W′t,......,W′m和LSTM编码器编码结果C1,C2,C3,......,Ct,......,Cn作为输入,利用注意力机制计算LSTM解码器所需的输入序列Att1,Att2,Att3,......,Attt,......,Attm,其中,Attt是注意力机制模块对LSTM编码器模型编码结果的加权和;
步骤3-3,利用LSTM解码器Decoder预测覆冰厚度信息;
步骤3-4,计算预测覆冰厚度与实际覆冰厚度的均方误差MSE(Θ);
步骤3-5,反向传播,使用随机梯度下降法对模型参数Θ进行优化;
步骤3-6,读取下一组序列化的微气象信息、覆冰厚度信息和天气预报信息,跳回到步骤3-1,直至将训练集中的数据读取完;
步骤3-7,重复执行步骤3-1至3-6并持续对模型的参数进行优化,直至使MSE(Θ)下降到合理的数量级。
4.根据权利要求3所述的方法,其特征在于,步骤3-1包括:
通过LSTM编码器对序列化的微气象、覆冰厚度信息和天气预报信息X1,X2,X3,......,Xt,......,Xn进行编码,输出序列化的隐变量序列H1,H2,H3,......,Ht,......,Hn和LSTM编码器编码结果C1,C2,C3,......,Ct,......,Cn,计算公式如(1)所示:
Ht,Ct=f(Ht-1,Xt) (1)
其中,Ht表示LSTM编码器输出的第t个隐变量;Ct是LSTM编码器第t个输出结果;Xt是三层全连接神经网络对第t个微气象、覆冰厚度信息和天气预报信息的编码值;H0是随机生成的隐变量,H0与X1作为f(·)函数首次的输入信息;f(·)是LSTM编码器模块。
5.根据权利要求4所述的方法,其特征在于,步骤3-2包括:
输入LSTM编码器编码结果C1,C2,C3,......,Ct,......,Cn和未来m个时间点的天气预报信息W′1,W′2,W′3,......,W′t,......,W′m,根据输入信息,计算LSTM编码器编码结果C1,C2,C3,......,Ct,......,Cn与天气预报信息W′1,W′2,W′3,......,W′t,......,W′m之间的相关系数Ai,j,计算公式如(2)(3)所示:
Ei,j=tanh(w(W′i)+u(Cj)) (2)
其中,tanh(·)是双曲正切函数;w(·)和u(·)分别表示一个三层的全连接神经网络;W‘i是未来序列上第i个时间点的天气预报信息;Cj是历史序列上第j个微气象、覆冰厚度信息和天气预报信息序列数据使用LSTM编码器的编码结果;Ei,j是利用注意力机制计算序列权重的中间变量,是公式(3)计算权重的依据;
exp(·)是计算指数值的函数;n代表预测未来天气所用到的历史微气象信息、覆冰厚度信息和天气预报信息的序列长度,Ei,j和Ei,k是公式(2)输出的中间变量;
根据权重系数Ai,j和LSTM编码器编码结果计算LSTM解码器输入序列Att1,Att2,Att3,......,Attt,......,Attm:
其中,Cj是历史序列上第j个微气象、覆冰厚度信息和天气预报信息序列数据使用LSTM编码器的编码结果;Ai,j是公式(3)输出的权重系数;最终输出序列化的特征向量Att1,Att2,Att3,......,Attt,......,Attm,Attt表示第t个特征向量,这些特征向量作为LSTM解码器的输入信息。
6.根据权利要求5所述的方法,其特征在于,步骤3-3包括:
输入序列化特征向量Attt和隐变量H′i-1,利用LSTM解码器预测序列化的特征向量Attt对应的输电线路覆冰厚度信息,输出对应的覆冰厚度序列计算过程如公式(5)所示:
其中,H′t表示LSTM编码器输出的第t个隐变量;是预测的第t个时间点的覆冰厚度;Attt是公式(4)求得的LSTM解码器输入信息;H′0是随机生成的隐变量;H′0和Att1作为h(·)函数首次的输入信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811416279.6A CN109492823B (zh) | 2018-11-26 | 2018-11-26 | 一种对输电线路覆冰厚度的预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811416279.6A CN109492823B (zh) | 2018-11-26 | 2018-11-26 | 一种对输电线路覆冰厚度的预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109492823A true CN109492823A (zh) | 2019-03-19 |
CN109492823B CN109492823B (zh) | 2021-04-30 |
Family
ID=65696740
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811416279.6A Active CN109492823B (zh) | 2018-11-26 | 2018-11-26 | 一种对输电线路覆冰厚度的预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109492823B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110070530A (zh) * | 2019-04-19 | 2019-07-30 | 山东大学 | 一种基于深度神经网络的输电线路覆冰检测方法 |
CN110321601A (zh) * | 2019-06-14 | 2019-10-11 | 山东大学 | 一种架空线路动态载流能力超前预测方法及系统 |
CN110852563A (zh) * | 2019-10-09 | 2020-02-28 | 天津大学 | 基于深度学习注意力机制的社会风险多目标优化生成方法 |
CN111458769A (zh) * | 2020-05-26 | 2020-07-28 | 南京大学 | 用于输电线路环境气象数据预测的方法及系统 |
CN111696097A (zh) * | 2020-06-15 | 2020-09-22 | 国网山东省电力公司平邑县供电公司 | 一种电缆在线监测方法和装置 |
CN112183897A (zh) * | 2020-11-02 | 2021-01-05 | 成都卡普数据服务有限责任公司 | 基于深度学习的架空输电线路覆冰厚度的长时间预测方法 |
CN113065253A (zh) * | 2021-04-01 | 2021-07-02 | 贵州电网有限责任公司 | 一种基于深度学习的输电线路覆冰类型预测方法 |
CN113128888A (zh) * | 2021-04-26 | 2021-07-16 | 国网湖北省电力有限公司宜昌供电公司 | 基于覆冰特征变量分箱评分卡的输电线路覆冰预测方法 |
CN113358040A (zh) * | 2021-05-14 | 2021-09-07 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 输电线路覆冰厚度获取方法、装置和计算机设备 |
CN113723669A (zh) * | 2021-08-09 | 2021-11-30 | 贵州电网有限责任公司 | 一种基于Informer模型的输电线路覆冰预测方法 |
CN112183897B (zh) * | 2020-11-02 | 2024-06-28 | 成都卡普数据服务有限责任公司 | 基于深度学习的架空输电线路覆冰厚度的长时间预测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103673960A (zh) * | 2012-08-30 | 2014-03-26 | 国际商业机器公司 | 用于预测输电线路上的覆冰状态的方法和装置 |
CN107092983A (zh) * | 2017-04-11 | 2017-08-25 | 北京国网富达科技发展有限责任公司 | 输电导线覆冰厚度预测方法及装置 |
CN107784395A (zh) * | 2017-10-27 | 2018-03-09 | 云南电网有限责任公司电力科学研究院 | 一种基于lstm人工神经网络的输电线路覆冰厚度预测方法 |
CN108154435A (zh) * | 2017-12-26 | 2018-06-12 | 浙江工业大学 | 一种基于循环神经网络的股指价格预测方法 |
CN108345955A (zh) * | 2017-11-24 | 2018-07-31 | 北京国网富达科技发展有限责任公司 | 一种输电线路覆冰短期预测方法、装置及存储介质 |
-
2018
- 2018-11-26 CN CN201811416279.6A patent/CN109492823B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103673960A (zh) * | 2012-08-30 | 2014-03-26 | 国际商业机器公司 | 用于预测输电线路上的覆冰状态的方法和装置 |
CN107092983A (zh) * | 2017-04-11 | 2017-08-25 | 北京国网富达科技发展有限责任公司 | 输电导线覆冰厚度预测方法及装置 |
CN107784395A (zh) * | 2017-10-27 | 2018-03-09 | 云南电网有限责任公司电力科学研究院 | 一种基于lstm人工神经网络的输电线路覆冰厚度预测方法 |
CN108345955A (zh) * | 2017-11-24 | 2018-07-31 | 北京国网富达科技发展有限责任公司 | 一种输电线路覆冰短期预测方法、装置及存储介质 |
CN108154435A (zh) * | 2017-12-26 | 2018-06-12 | 浙江工业大学 | 一种基于循环神经网络的股指价格预测方法 |
Non-Patent Citations (2)
Title |
---|
AN YUAN 等: "Prediction of Icing Thickness on Transmission Lines using ANFIS model", 《2017 GLOBAL CONFERENCE ON MECHANICS AND GIVIL ENGINEERIGN》 * |
吴琼 等: "基于RBF神经网络的输电线路覆冰短期预测研究", 《贵州电力技术》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110070530A (zh) * | 2019-04-19 | 2019-07-30 | 山东大学 | 一种基于深度神经网络的输电线路覆冰检测方法 |
CN110321601A (zh) * | 2019-06-14 | 2019-10-11 | 山东大学 | 一种架空线路动态载流能力超前预测方法及系统 |
CN110852563A (zh) * | 2019-10-09 | 2020-02-28 | 天津大学 | 基于深度学习注意力机制的社会风险多目标优化生成方法 |
CN111458769A (zh) * | 2020-05-26 | 2020-07-28 | 南京大学 | 用于输电线路环境气象数据预测的方法及系统 |
CN111696097A (zh) * | 2020-06-15 | 2020-09-22 | 国网山东省电力公司平邑县供电公司 | 一种电缆在线监测方法和装置 |
CN112183897A (zh) * | 2020-11-02 | 2021-01-05 | 成都卡普数据服务有限责任公司 | 基于深度学习的架空输电线路覆冰厚度的长时间预测方法 |
CN112183897B (zh) * | 2020-11-02 | 2024-06-28 | 成都卡普数据服务有限责任公司 | 基于深度学习的架空输电线路覆冰厚度的长时间预测方法 |
CN113065253A (zh) * | 2021-04-01 | 2021-07-02 | 贵州电网有限责任公司 | 一种基于深度学习的输电线路覆冰类型预测方法 |
CN113128888A (zh) * | 2021-04-26 | 2021-07-16 | 国网湖北省电力有限公司宜昌供电公司 | 基于覆冰特征变量分箱评分卡的输电线路覆冰预测方法 |
CN113358040A (zh) * | 2021-05-14 | 2021-09-07 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 输电线路覆冰厚度获取方法、装置和计算机设备 |
CN113358040B (zh) * | 2021-05-14 | 2024-01-23 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 输电线路覆冰厚度获取方法、装置和计算机设备 |
CN113723669A (zh) * | 2021-08-09 | 2021-11-30 | 贵州电网有限责任公司 | 一种基于Informer模型的输电线路覆冰预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109492823B (zh) | 2021-04-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109492823A (zh) | 一种对输电线路覆冰厚度的预测方法 | |
CN106548645B (zh) | 基于深度学习的车辆路径寻优方法及系统 | |
CN104408900B (zh) | 基于动态优化的神经网络洪水预警装置及其方法 | |
Thomas et al. | Predictions of urban volumes in single time series | |
CN102938021B (zh) | 一种输电线路覆冰负荷定量估计及预测方法 | |
CN102221714B (zh) | 基于大气变量物理分解的低温雨雪冰冻天气的中期预报系统和方法 | |
CN109658695A (zh) | 一种多因素的短时交通流预测方法 | |
CN107516145A (zh) | 一种基于加权欧氏距离模式分类的多通道光伏出力预测方法 | |
CN113177657B (zh) | 轨道交通客流预测方法及装置 | |
CN102930155B (zh) | 获取电力需求的预警参数的方法及装置 | |
CN110211325A (zh) | 一种基于气象大数据的区域道路结冰高精度监测预警系统 | |
CN111311027A (zh) | 一种用于输电线路的风速预测方法 | |
CN105608326A (zh) | 一种山区复杂地形风场大涡模拟入口边界条件输入方法 | |
CN107330243A (zh) | 一种基于典型气象年的架空输电线路载流量计算方法 | |
CN112925825A (zh) | 用于输电线路的多气象因子预测方法 | |
CN105846425A (zh) | 一种基于通用风电预测误差模型的经济调度方法 | |
CN103793887A (zh) | 基于自适应增强算法的短期电力负荷在线预测方法 | |
CN116070676B (zh) | 基于注意力机制和自编码器的高速公路路面温度预报方法 | |
CN109447322A (zh) | 一种基于相空间重构小波神经网络输电线路覆冰预测的方法 | |
CN103914740B (zh) | 一种基于数据驱动的输电线路覆冰预测与自校正方法 | |
CN103914737B (zh) | 一种输变电线路全线路现地气象信息计算方法 | |
CN114021830A (zh) | 一种基于cnn-lstm的多时间范围风速预测方法 | |
CN114357670A (zh) | 一种基于bls和自编码器的配电网用电数据异常预警方法 | |
Zhou et al. | The development of next NCEP global ensemble forecast system | |
CN107403004A (zh) | 一种基于地形数据的遥测雨量站点可疑数值检验方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |