CN111242377A - 一种集成深度学习和数据去噪的短期风速预测方法 - Google Patents
一种集成深度学习和数据去噪的短期风速预测方法 Download PDFInfo
- Publication number
- CN111242377A CN111242377A CN202010042211.7A CN202010042211A CN111242377A CN 111242377 A CN111242377 A CN 111242377A CN 202010042211 A CN202010042211 A CN 202010042211A CN 111242377 A CN111242377 A CN 111242377A
- Authority
- CN
- China
- Prior art keywords
- wind speed
- prediction
- sequence
- wavelet
- wavelet transform
- 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 53
- 238000013135 deep learning Methods 0.000 title claims abstract description 16
- 238000013528 artificial neural network Methods 0.000 claims abstract description 10
- 230000009466 transformation Effects 0.000 claims abstract description 9
- 238000012549 training Methods 0.000 claims abstract description 7
- 230000006870 function Effects 0.000 claims description 24
- 210000002569 neuron Anatomy 0.000 claims description 14
- 239000011159 matrix material Substances 0.000 claims description 13
- 238000000354 decomposition reaction Methods 0.000 claims description 12
- 238000007781 pre-processing Methods 0.000 claims description 11
- 238000010606 normalization Methods 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 3
- 230000021615 conjugation Effects 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims 1
- 230000000694 effects Effects 0.000 description 8
- 230000009467 reduction Effects 0.000 description 7
- 238000002474 experimental method Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 230000006872 improvement Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000000306 recurrent effect Effects 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000010801 machine learning Methods 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 241000288105 Grus Species 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000000739 chaotic effect Effects 0.000 description 1
- 238000013527 convolutional neural network Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000002803 fossil fuel Substances 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000000053 physical method Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/21—Design, administration or maintenance of databases
- G06F16/215—Improving data quality; Data cleansing, e.g. de-duplication, removing invalid entries or correcting typographical errors
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/24—Querying
- G06F16/245—Query processing
- G06F16/2458—Special types of queries, e.g. statistical queries, fuzzy queries or distributed queries
- G06F16/2474—Sequence data queries, e.g. querying versioned data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/048—Activation functions
-
- 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
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Economics (AREA)
- Data Mining & Analysis (AREA)
- Health & Medical Sciences (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Computational Linguistics (AREA)
- Mathematical Physics (AREA)
- Databases & Information Systems (AREA)
- General Health & Medical Sciences (AREA)
- Software Systems (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Quality & Reliability (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- General Business, Economics & Management (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Fuzzy Systems (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Public Health (AREA)
- Entrepreneurship & Innovation (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种集成深度学习和数据去噪的短期风速预测方法,步骤为:1)利用离散小波变换对风速序列A进行分解。2)利用小波软阈值去噪方法区分离散小波变换后的噪声信号和有效信号;3)将噪声信号的小波变换系数设为零,并重建风速时间序列,得到去噪后的风速序列B;4)利用风速序列B训练门控循环单元神经网络,得到风速预测模型;6)将实时风速输入到风速预测模型中,完成未来多步风速预测。本发明可以实现短期风速高质量的预测,以确保电力系统的经济调度和安全运行。
Description
技术领域
本发明涉及风电功率预测领域,具体是一种集成深度学习和数据去噪的短期风速预测方法。
背景技术
在当今能源危机的背景下,可再生能源的开发和利用已成为政府和研究机构关注的焦点。截至2017年底,全球可再生能源装机容量已连续三年保持年均8%的增长率,装机容量已达2179GW。其中,风电装机容量是可再生能源的重要组成部分,占据第二大份额(23.59%)。同时,2017年风电装机容量增加94GW,同比增长10%。与水力、化石燃料、核能等其他能源相比,风力的大小主要取决于风速、风向、大气压等自然条件。而风速具有随机性、波动性和不可预测性,因此,风速的预测值难以准确获取。风速预测不佳,导致我国等一些国家采用了风力限电政策,这是对可再生能源的巨大浪费,大大降低了风电场业主的利润,从而极大地抑制了风电产业的发展。为了灵活地控制风力输出,为电力调度员提供方便,需要开发出有效的风速预测方法。
风速预测方法一般分为两大类:统计学习方法和物理方法。物理预报模型通过建立风的气动模型(例如方程组)来计算未来的风速值,从而能够在长期预报中得到准确的风速预报。统计学习方法通过风速的历史数据来预测风速值,在短期风速预测中效果较好。统计学习方法可分为两类:单一预测方法和整体预测方法(即预测模型与数据预处理技术相结合)。单一预测方法本身存在局限性,而整体预测方法可以更好地从数据中提取特征,因此,单一预测方法普遍不如整体预测方法。所以,整体预测方法被广泛地用于风速预测中,这也是本专利研究的重点。
但是,传统的机器学习方法缺乏深度提取能力,无法从数据中有效地提取时间序列信息。一种可行的解决方案便是采用深度学习方法。近年来,随着计算能力的提高和大数据技术的进步,深度学习方法得到了快速发展。深度学习的优势之一在于,与其他机器学习方法相比,它可以在每个计算步骤上自动构造更复杂的特征。在风速预测中,深度学习是提高风速预测精度的有效方法。但是,CNN或LSTM的模型参数通常很大,使得模型训练困难,容易被过度拟合。一种新型循环神经网络(RNN)方法门控循环单元(GRU),专门用于解决上述问题。在LSTM的基础上,对GRU的结构进行优化、集成和升级,有效地提高了GRU的预测速度。而有关于用于风速预测的GRU的研究报道很少。而且,以往的许多研究只关注了预测精度,忽略了预测效率和结果波动性的重要性,这是需要改进的。因此,便选择GRU作为预测方法。结果的波动性表明,如果重复执行预测过程,则每个预测过程获得的结果的大小都会变化。
一方面,模型的预测速度很重要。模型的预测速度越快,风电企业可以更快地得到预测风速值,并根据最新获得的数据更新预测结果,从而制定出更好的风力发电调度方案。另一方面,结果的小波动性也很重要,它反映了一个算法在一个预测时间内获得的预测结果的可信度和可靠性。从实践的角度看,一个好的预测模型必须具备对不同地点良好的适应性,而关于风场预测模型适应性的研究报道较少。
除了预测算法之外,集成模型的另一个关键技术在于数据预处理。目前,广泛应用的数据预处理方法是基于分解的方法(如经验模式分解,EMD)。然而,新获得的数据将极大地影响初始分解子序列的值,因此基于EMD的预测方法在实际预测中很难使用,这导致在通过逆变换将结果的子序列转换回实际值之后,最终预测的精度降低。
发明内容
本发明的目的是解决现有技术中存在的问题。
为实现本发明目的而采用的技术方案是这样的,一种集成深度学习和数据去噪的短期风速预测方法,主要包括以下步骤:
1)获取原始风速序列,并对所述原始风速序列进行预处理,得到风速序列A:{α1,α2,…,αl}。l为数据序列的长度。所述预处理包括线性插值和归一化。
2)利用离散小波变换对风速序列A进行分解,得到包含低频信息的比例系数cm,n和包含高频信息的小波系数dm,n。
利用离散小波变换对风速序列A进行分解的主要步骤如下:
2.1)设置小波函数为二进小波,即令a=2m,b=n2m,m,n∈Z。a为尺度因子,b为时移因子,m为尺度,n是与时移有关的参数。
2.2)利用尺度因子a、时移因子b对母小波ψ(·)进行伸缩和时移变换,得到小波序列ψm,n(t),即:
2.3)基于公式(1)和小波序列定义公式,计算离散小波变换系数Wx,即:
2.4)利用Mallat算法对输入x(t)进行多级离散小波变换分解,得到包含低频信息的比例系数cm,n和包含高频信息的小波系数dm,n。
分解公式如下:
式中,φ(·)为小波函数。
比例系数cm,n和小波系数dm,n分别如下所示:
3)设置阈值λ,并利用小波软阈值去噪方法区分离散小波变换后的噪声信号和有效信号。
区分离散小波变换后的噪声信号和有效信号的方法为:判断信号的小波变换系数绝对值|Wx|≥λ是否成立,若成立,则对应信号为有效信号,若不成立,则对应信号为噪声信号。区分后,将噪声信号的小波变换系数设为0,即:
4)将噪声信号的小波变换系数设为零,并重建风速时间序列,得到去噪后的风速序列B:{b1,b2,...,bl}。
5)利用风速序列B训练门控循环单元神经网络,得到风速预测模型,主要步骤如下:
5.1)建立风速序列矩阵X:{x1,x2,...,xl},即:
式中,n为预测次数,m为预测器的输入次数。l=n+m-1。xi=bi,i=1,2,…,l。
5.2)提取风速序列矩阵X的第i行作为输入向量Xi,即:
Xi={xi,xi+1,…,xi+m-1}1×m。 (8)
式中,i初始值为1。
5.3)基于门控循环单元神经网络,建立预测函数f(·)。预测函数f(·)的更新门zt、复位门(即重置门)rt如下所示:
zt=σ(Wxzxt+Whzht-1+bz)。 (9)
rt=σ(Wxrxt+Whrht-1+br)。 (10)
式中,当前神经元的输入xt为当前t时刻的输入向量Xi。Wxz为当前神经元输入与更新门之间的权值。Whz为隐藏状态与更新门之间的权值。ht-1表示前一个神经元的输出。bz是更新门的偏置。σ为sigmoid函数。Wxr为当前神经元输入与复位门之间的权值。Whr为隐藏状态与复位门之间的权值。br是复位门的偏置。
当前神经元的输出ht如下所示:
5.4)将输入向量Xi输入到预测函数f(·)中,得到输出状态yt,即:
yt=σ(Whoht+bh)。 (12)
式中,Who为隐藏状态与最终输出状态之间的权值。bh是隐态的偏置。
5.5)将输出状态yt赋值给第i个预测时刻的预测输出yi+m,得到:
yi+m=f(Xi)=f({xi,...,xi+m-1}),i=1,2,…,n。 (13)
5.6)令i=i+1,重复步骤5.2至5.4,直至风速序列矩阵每行向量均预测完毕,输出单步风速预测结果Youtput:{y1+m,y2+m,…,yn+m}1×n。
5.7)基于单步时间序列预测,进行多步时间序列的预测。
预测器f(Xi,c)的单轮输出,即第i次预测时第c步预测结果yi+m+(c-1)如下所示:
yi+m+(c-1)=f(Xi,c)=f({xi+(c-1),xi+c,…,xm,yi+m,…yi+m+(c-2)})。(14)
式中,c为预测超前步骤。Xi,c为第i次预测时第c步时预测器的输入向量。
5.8)利用公式(10)对风速序列矩阵每行向量进行预测,整合得到风速预测结果Youtput:{y1+m+(c-1),y2+m+(c-1)…,yn+m+(c-1)}1×n。
6)将实时风速输入到风速预测模型中,完成未来多步风速预测。
本发明的技术效果是毋庸置疑的。本发明在预测精度、预测效率、结果波动性和适应性等方面的综合研究具有创新性和时效性。
本发明采用小波阈值降噪可以滤除冗余的高频信号,为预测提供高质量的数据。本发明利用WSTD方法进行风速预测,解决了针对电力系统的短期风速预测的基于分解的预处理方法会收到最新获取的数据会极大地影响从而在实际的预测中仍很难使用的问题,实现短期风速高质量的预测,以确保电力系统的经济调度和安全运行。
附图说明
图1算法流程图;
图2基于Mallat算法的多层次离散小波分解;
图3循环神经网络的结构。
具体实施方式
下面结合实施例对本发明作进一步说明,但不应该理解为本发明上述主题范围仅限于下述实施例。在不脱离本发明上述技术思想的情况下,根据本领域普通技术知识和惯用手段,做出各种替换和变更,均应包括在本发明的保护范围内。
实施例1:
参见图1至图3,一种集成深度学习和数据去噪的短期风速预测方法,主要包括以下步骤:
1)获取原始风速序列,并对所述原始风速序列进行预处理,得到风速序列A:{α1,α2,…,αl}。l为数据序列的长度。所述预处理包括线性插值和归一化。αl为风速。
2)利用离散小波变换对风速序列A进行分解,得到包含低频信息的比例系数cm,n和包含高频信息的小波系数dm,n。
利用离散小波变换对风速序列A进行分解的主要步骤如下:
2.1)设置小波函数为二进小波,即令a=2m,b=n2m,m,n∈Z。a为尺度因子,b为时移因子,m为尺度,n是与时移有关的参数。
2.2)利用尺度因子a、时移因子b对母小波ψ(·)进行伸缩和时移变换,得到小波序列ψm,n(t),即:
2.3)基于公式(1)和小波序列定义公式,计算离散小波变换系数Wx,即:
2.4)利用Mallat算法对输入x(t)进行多级离散小波变换分解,得到包含低频信息的比例系数cm,n和包含高频信息的小波系数dm,n。
分解公式如下:
式中,φ(·)为小波函数。
比例系数cm,n和小波系数dm,n分别如下所示:
3)设置阈值λ,并利用小波软阈值去噪方法区分离散小波变换后的噪声信号和有效信号。
区分离散小波变换后的噪声信号和有效信号的方法为:判断信号的小波变换系数绝对值|Wx|≥λ是否成立,若成立,则对应信号为有效信号,若不成立,则对应信号为噪声信号。区分后,将噪声信号的小波变换系数设为0,即:
4)将噪声信号的小波变换系数设为零,并重建风速时间序列,得到去噪后的风速序列B:{b1,b2,…,bl}。
5)利用风速序列B训练门控循环单元神经网络,得到风速预测模型,主要步骤如下:
5.1)建立风速序列矩阵X:{x1,x2,...,xl},即:
式中,n为预测次数,m为预测器的输入次数。l=n+m-1。xi=bi,i=1,2,…,l。
5.2)提取风速序列矩阵X的第i行作为输入向量Xi,即:
Xi={xi,xi+1,…,xi+m-1}1×m。 (8)
式中,i初始值为1。
5.3)基于门控循环单元神经网络,建立预测函数f(·),实现原理如下:
zt=σ(Wxzxt+Whzht-1+bz) (9)
式中,当前神经元的输入xt为当前t时刻的输入向量Xi。Wxz为当前神经元输入与更新门之间的权值。Whz为隐藏状态与更新门之间的权值。ht-1表示前一个神经元的输出。bz是更新门的偏置。σ为sigmoid函数。Wxr为当前神经元输入与复位门之间的权值。Whr为隐藏状态与复位门之间的权值。br是复位门的偏置。
首先,当前t时刻的输入向量Xi作为当前神经元的输入xt。xt先通过更新门zt的计算。对于更新门,当它的值较大时,当前神经元需要保留更多的信息,前一个神经元需要保留较少的信息。
此外,复位门rt代表先前隐藏状态的信息被遗忘的程度,其中Wxr为当前神经元输入与复位门之间的权值;Whr为隐藏状态与复位门之间的权值;br是复位门的偏置。
σ表示sigmoid函数可计算为:
之后就可以计算该神经元的隐藏状态ht。
最后确定最终输出状态yt为:
yt=σ(Whoht+bh) (14)
其中Who为隐藏状态与最终输出状态之间的权值;bh是隐态的偏置。
5.5)输出状态yt即可赋值给下一步骤中的预测输出yi+m,得到:
yi+m=f(Xi)=f({xi,...,xi+m-1}),i=1,2,…,n。 (15)
5.6)令i=i+1,重复步骤5.2至5.4,直至风速序列矩阵每行向量均预测完毕,输出单步风速预测结果Youtput:{y1+m,y2+m,…,yn+m}1×n。
5.7)基于单步时间序列预测,进行多步时间序列的预测。
预测器f(Xi,c)的单轮输出,即第i次预测时第c步预测结果yi+m+(c-1)如下所示:
yi+m+(c-1)=f(Xi,c)=f({xi+(c-1),xi+c,…,xm,yi+m,…yi+m+(c-2)})。(16)
式中,c为预测超前步骤。Xi,c为第i次预测时第c步时预测器的输入向量。
5.8)利用公式(10)对风速序列矩阵每行向量进行预测,整合得到风速预测结果Youtput:{y1+m+(c-1),y2+m+(c-1)…,yn+m+(c-1)}1×n。
6)将实时风速输入到风速预测模型中,完成未来多步风速预测。
实施例2:
一种验证集成深度学习和数据去噪的短期风速预测方法的实验,主要包括以下步骤:
1)获取原始风速序列;通过数据预处理对收集的原始风速序列进行缺失值的线性插值和数据的归一化,得到整理后的风速序列A:{α1,α2,…,αl}(l为数据序列的长度)。
由于风速原始数据存在许多丢失和错误的值,因此使用线性插值技术对数据进行预处理:首先,将缺失的点和错误的点替换为名为“非数字(NaN)”的标志。然后,采用线性插值方法对NaN位置进行修补。同样,将分辨率为几分钟的原始数据序列转换为分辨率为小时的数据。此外,这些时间序列被归一化至-1到1的区间内。
2)利用离散小波变换对风速原始时间序列进行分解;用二进小波作为小波函数进行离散小波变换,则有a=2m,b=n2m,m,n∈Z;其中,a为尺度因子,b为时移因子,m为尺度,n是与时移有关的参数;
5)通过Mallat算法进行多级离散小波变换分解,分解后得到相应的包含低频信息的比例系数和包含高频信息的小波系数,如图2所示,即,
6)通过计算阈值λ,利用小波软阈值去噪(WSTD),从而区分噪声信号和有效信号;即,
7)将噪声的小波变换系数设为零,重建风速时间序列得到去噪后的风速序列B:{b1,b2,…,bl};
8)衡量WSTD的降噪效果:采用Duffing混沌系统来获得包含1000个数据点的时间序列。该系统的公式为:
采用Runge-Kutta法求解,分别设定参数σ,a,F,δ和时间步长h。变量x作为实验值。为了消除暂态点对系统的影响,放弃了前30000个暂态点。从而得到了以下1000个数据点的时间序列。
信号噪声比(SNR)用于评估数字降噪的效果。SNR值越大,意味着信号中会有更多有用的信息(即更少的无用噪声)。SNR值越大,也意味着相应的降噪方法更有可能在风速预测中有更大的精度改进。SNR的定义可以表示为:
其中,H1和H2分别为:
其中,sp为无噪声的纯信号,sn为有高斯噪声的信号。
在Duffing系统的x值中加入一定的高斯噪声,分别采用低通滤波技术、WHTD技术和WSTD技术对高斯噪声信号进行降噪。降噪效果对比如表1所示。数值结果表明,该方法具有较好的降噪效果。
表1去噪效果的比较
9)利用降噪时间序列训练门控循环单元(GRU)神经网络,如图3所示,得到一步预测未来多步预测模型;图3中a部分是RNN结构的三维示意图;b部分是单隐层RNN的二维示意图;c部分是RNN在时域上的扩展示意图;d部分是GRU块(作为RNN的一种类型)的结构示意图。
单步时间序列预测:已知现有输入风速系列X:{x1,x2,…,xl}且xi=bi,i=1,2,…,l,则:
其中,n为预测次数,m为预测器的输入次数,满足l=n+m-1;风速系列的测试数据并不包含在训练数据集中。
提出矩阵的一行,得到输入向量:
Xi={xi,xi+1,…,xi+m-1}1×m;
假设预测器是“f(·)”,则:
yi+m=f(Xi)=f({xi,…,xi+m-1}),i=1,2,…,n
其中,yi+m为第i个预测时刻的预测输出;
得到风速预测的最终输出向量:
Youtput:{y1+m,y2+m,…,yn+m}1×n
则单步时间序列的预测完成。
基于单步时间序列预测,进行多步时间序列的预测;采用递推公式可知预测器的输出表示为:
yi+m+(c-1)=f(Xi,c)=f({xi+(c-1),xi+c,…,xm,yi+m,…yi+m+(c-2)})
其中,c为预测超前步骤(如本文实验中c={1,2,3});yi+m+(c-1)为第i次预测实验第c步的预测结果,为预测器的单轮输出;Xi,c为第i个预测实验第c步时预测器的输入向量;
最后汇总所有的yi+m+(c-1),得到风速预测的最终结果向量为:
Youtput:{y1+m+(c-1),y2+m+(c-1)…,yn+m+(c-1)}1×n。
10)此外,为了评价结果的波动性,将所有预测方法的过程重复20次,以获得标准的平均值和标准的标准偏差。预测实验精度的量化有三个最常用的标准:平均绝对百分比误差(MAPE)、平均绝对误差(MAE)和均方根误差(RMSE)。
结果如表2所示。因此,与其他八种常用的预测模型相比,WSTD-GRU方法具有较高的多步预测精度,并且能够在美国从不同经度获得的所有数据集中实现出良好的效果。尤其是通过与WSTD结合,GRU具有更好的预测性能。
表2多步预测结果精度的比较
9)采用训练后的WSTD-GRU模型(即风速预测模型),进行未来多步风速预测。
Claims (5)
1.一种集成深度学习和数据去噪的短期风速预测方法,其特征在于,主要包括以下步骤:
1)获取所述原始风速序列,并对所述原始风速序列进行预处理,得到风速序列A:{α1,α2,…,αl};l为数据序列的长度。
2)利用离散小波变换对风速序列A进行分解,得到包含低频信息的比例系数cm,n和包含高频信息的小波系数dm,n。
3)设置阈值λ,并利用小波软阈值去噪方法区分离散小波变换后的噪声信号和有效信号;
4)将噪声信号的小波变换系数设为零,并重建风速时间序列,得到去噪后的风速序列B:{b1,b2,…,bl};
5)利用风速序列B训练门控循环单元神经网络,得到风速预测模型;
6)将实时风速输入到风速预测模型中,完成未来多步风速预测。
2.根据权利要求1或2所述的一种集成深度学习和数据去噪的短期风速预测方法,其特征在于:所述预处理包括线性插值和归一化。
3.根据权利要求1或2所述的一种集成深度学习和数据去噪的短期风速预测方法,其特征在于,利用离散小波变换对风速序列A进行分解的主要步骤如下:
1)设置小波函数为二进小波,即令a=2m,b=n2m,m,n∈Z;a为尺度因子,b为时移因子,m为尺度,n是与时移有关的参数;
2)利用尺度因子a、时移因子b对母小波ψ(·)进行伸缩和时移变换,得到小波序列ψm,n(t),即:
3)基于公式(1)和小波序列定义公式,计算离散小波变换系数Wx,即:
4)利用Mallat算法对输入x(t)进行多级离散小波变换分解,得到包含低频信息的比例系数cm,n和包含高频信息的小波系数dm,n;
分解公式如下:
式中,φ(·)为小波函数;
比例系数cm,n和小波系数dm,n分别如下所示:
5.根据权利要求1所述的一种集成深度学习和数据去噪的短期风速预测方法,其特征在于,建立风速预测模型的主要步骤如下:
1)建立风速序列矩阵X:{x1,x2,…,xl},即:
式中,n为预测次数,m为预测器的输入次数;l=n+m-1;xi=bi,i=1,2,…,l;
2)提取风速序列矩阵X的第i行作为输入向量Xi,即:
Xi={xi,xi+1,…,xi+m-1}1×m; (8)
式中,i初始值为1;
3)基于门控循环单元神经网络,建立预测函数f(·);预测函数f(·)的更新门zt、复位门rt如下所示:
zt=σ(Wxzxt+Whzht-1+bz); (9)
rt=σ(Wxrxt+Whrht-1+br); (10)
式中,当前神经元的输入xt为当前t时刻的输入向量Xi;Wxz为当前神经元输入与更新门之间的权值;Whz为隐藏状态与更新门之间的权值;ht-1表示前一个神经元的输出;bz是更新门的偏置;σ为sigmoid函数;Wxr为当前神经元输入与复位门之间的权值;Whr为隐藏状态与复位门之间的权值;br是复位门的偏置;
当前神经元的输出ht如下所示:
4)将输入向量Xi输入到预测函数f(·)中,得到输出状态yt,即:
yt=σ(Whoht+bh); (12)
式中,Who为隐藏状态与最终输出状态之间的权值;bh是隐态的偏置;
5)将输出状态yt赋值给第i个预测时刻的预测输出yi+m,得到:
yi+m=f(Xi)=f({xi,…,xi+m-1}),i=1,2,…,n; (13)
6)令i=i+1,重复步骤2至4,直至风速序列矩阵每行向量均预测完毕,输出单步风速预测结果Youtput:{y1+m,y2+m,…,yn+m}1×n;
7)基于单步时间序列预测,进行多步时间序列的预测;
预测器f(Xi,c)的单轮输出,即第i次预测时第c步预测结果yi+m+(c-1)如下所示:
yi+m+(c-1)=f(Xi,c)=f({xi+(c-1),xi+c,…,xm,yi+m,...yi+m+(c-2)}); (14)
式中,c为预测超前步骤;Xi,c为第i次预测时第c步时预测器的输入向量;
8)利用公式(10)对风速序列矩阵每行向量进行预测,整合得到风速预测结果Youtput:{y1+m+(c-1),y2+m+(c-1)...,yn+m+(c-1)}1×n。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010042211.7A CN111242377B (zh) | 2020-01-15 | 2020-01-15 | 一种集成深度学习和数据去噪的短期风速预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010042211.7A CN111242377B (zh) | 2020-01-15 | 2020-01-15 | 一种集成深度学习和数据去噪的短期风速预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111242377A true CN111242377A (zh) | 2020-06-05 |
CN111242377B CN111242377B (zh) | 2023-08-25 |
Family
ID=70877826
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010042211.7A Active CN111242377B (zh) | 2020-01-15 | 2020-01-15 | 一种集成深度学习和数据去噪的短期风速预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111242377B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111709588A (zh) * | 2020-06-23 | 2020-09-25 | 苏州联电能源发展有限公司 | 一种用电量预测方法和系统 |
CN112379274A (zh) * | 2020-11-16 | 2021-02-19 | 河南科技大学 | 一种动力电池剩余寿命预测方法 |
CN112509602A (zh) * | 2020-12-18 | 2021-03-16 | 平安科技(深圳)有限公司 | 一种居家监护方法、装置、设备及存储介质 |
CN112559827A (zh) * | 2020-12-08 | 2021-03-26 | 上海上实龙创智能科技股份有限公司 | 一种基于深度学习的测量参数预测及污水处理控制方法 |
CN112649623A (zh) * | 2020-11-30 | 2021-04-13 | 中国科学院声学研究所南海研究站 | 一种环境自适应的超声测风系统 |
CN112801332A (zh) * | 2020-11-18 | 2021-05-14 | 国网江苏省电力有限公司江阴市供电分公司 | 一种基于灰度共生矩阵的短期风速预测方法 |
CN113282990A (zh) * | 2021-05-31 | 2021-08-20 | 中铁十六局集团北京轨道交通工程建设有限公司 | 一种盾构运动轨迹的智能化实时确定方法、终端和介质 |
CN113537586A (zh) * | 2021-07-12 | 2021-10-22 | 浙江大学 | 集成小波分解与深度神经网络的水质指标多步预测方法 |
WO2022033258A1 (zh) * | 2020-08-14 | 2022-02-17 | 中国科学院分子细胞科学卓越创新中心 | 基于自动储层神经网络的多步预测未来风速的方法和系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1991018362A1 (en) * | 1990-05-18 | 1991-11-28 | Yale University | Method and apparatus for machine processing using wavelet/wavelet-packet bases |
KR20140018497A (ko) * | 2012-08-01 | 2014-02-13 | 한국전력공사 | 단기 풍속 예측을 통한 풍력 발전량 예측 방법과, 그 기능을 이용한 배전선로 전압 예측 방법 |
CN109146162A (zh) * | 2018-08-07 | 2019-01-04 | 河海大学 | 一种基于集成循环神经网络的概率风速预测方法 |
CN110414045A (zh) * | 2019-06-18 | 2019-11-05 | 东华大学 | 基于vmd-gru的短期风速预测方法 |
-
2020
- 2020-01-15 CN CN202010042211.7A patent/CN111242377B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1991018362A1 (en) * | 1990-05-18 | 1991-11-28 | Yale University | Method and apparatus for machine processing using wavelet/wavelet-packet bases |
KR20140018497A (ko) * | 2012-08-01 | 2014-02-13 | 한국전력공사 | 단기 풍속 예측을 통한 풍력 발전량 예측 방법과, 그 기능을 이용한 배전선로 전압 예측 방법 |
CN109146162A (zh) * | 2018-08-07 | 2019-01-04 | 河海大学 | 一种基于集成循环神经网络的概率风速预测方法 |
CN110414045A (zh) * | 2019-06-18 | 2019-11-05 | 东华大学 | 基于vmd-gru的短期风速预测方法 |
Non-Patent Citations (1)
Title |
---|
薛阳等: "一种结合CNN 和GRU 网络的超短期风电预测模型", 《可再生能源》, vol. 37, no. 3, pages 456 - 462 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111709588A (zh) * | 2020-06-23 | 2020-09-25 | 苏州联电能源发展有限公司 | 一种用电量预测方法和系统 |
CN111709588B (zh) * | 2020-06-23 | 2023-08-15 | 苏州联电能源发展有限公司 | 一种用电量预测方法和系统 |
WO2022033258A1 (zh) * | 2020-08-14 | 2022-02-17 | 中国科学院分子细胞科学卓越创新中心 | 基于自动储层神经网络的多步预测未来风速的方法和系统 |
CN112379274A (zh) * | 2020-11-16 | 2021-02-19 | 河南科技大学 | 一种动力电池剩余寿命预测方法 |
CN112801332A (zh) * | 2020-11-18 | 2021-05-14 | 国网江苏省电力有限公司江阴市供电分公司 | 一种基于灰度共生矩阵的短期风速预测方法 |
CN112801332B (zh) * | 2020-11-18 | 2024-03-26 | 国网江苏省电力有限公司江阴市供电分公司 | 一种基于灰度共生矩阵的短期风速预测方法 |
CN112649623A (zh) * | 2020-11-30 | 2021-04-13 | 中国科学院声学研究所南海研究站 | 一种环境自适应的超声测风系统 |
CN112559827A (zh) * | 2020-12-08 | 2021-03-26 | 上海上实龙创智能科技股份有限公司 | 一种基于深度学习的测量参数预测及污水处理控制方法 |
CN112509602A (zh) * | 2020-12-18 | 2021-03-16 | 平安科技(深圳)有限公司 | 一种居家监护方法、装置、设备及存储介质 |
CN113282990A (zh) * | 2021-05-31 | 2021-08-20 | 中铁十六局集团北京轨道交通工程建设有限公司 | 一种盾构运动轨迹的智能化实时确定方法、终端和介质 |
CN113537586A (zh) * | 2021-07-12 | 2021-10-22 | 浙江大学 | 集成小波分解与深度神经网络的水质指标多步预测方法 |
CN113537586B (zh) * | 2021-07-12 | 2024-04-05 | 浙江大学 | 集成小波分解与深度神经网络的水质指标多步预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111242377B (zh) | 2023-08-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111242377B (zh) | 一种集成深度学习和数据去噪的短期风速预测方法 | |
CN111860982B (zh) | 一种基于vmd-fcm-gru的风电场短期风电功率预测方法 | |
Zhang et al. | Short-term offshore wind power forecasting-A hybrid model based on Discrete Wavelet Transform (DWT), Seasonal Autoregressive Integrated Moving Average (SARIMA), and deep-learning-based Long Short-Term Memory (LSTM) | |
CN111091233B (zh) | 一种风电场短期风电预测建模方法 | |
CN111079683B (zh) | 基于卷积神经网络的遥感图像云雪检测方法 | |
CN114549925B (zh) | 一种基于深度学习的海浪有效波高时间序列预测方法 | |
CN115293415A (zh) | 计及时间演变和空间相关的多风电场短期功率预测方法 | |
CN112330951B (zh) | 一种基于生成对抗网络实现路网交通数据修复的方法 | |
CN112434891A (zh) | 基于wcnn-alstm的太阳辐照度时间序列的预测方法 | |
CN112288193A (zh) | 基于注意力机制的gru深度学习的海洋站表层盐度预测方法 | |
CN115688579A (zh) | 一种基于生成对抗网络的流域多点水位预测预警方法 | |
CN111242351A (zh) | 基于自编码器和gru神经网络的热带气旋轨迹预测方法 | |
CN117114192B (zh) | 基于多尺度时空图变换网络的海上风电功率预测方法及装置 | |
Yao et al. | Wave height forecast method with multi-step training set extension LSTM neural network | |
CN115204035A (zh) | 基于多尺度时序数据融合模型的发电机组运行参数预测方法、装置及存储介质 | |
CN116152611B (zh) | 一种多级多尺度点云补全方法、系统、设备及存储介质 | |
Cao et al. | UTrans-Net: A model for short-term precipitation prediction | |
CN117154680A (zh) | 一种基于非平稳Transformer模型的风电功率预测方法 | |
CN117033987A (zh) | 一种基于小波的风电场发电功率效率预测方法 | |
CN116702831A (zh) | 一种考虑数据大量丢失的混合短期风电功率预测方法 | |
Wang et al. | Filling gaps in significant wave height time series records using bidirectional gated recurrent unit and cressman analysis | |
Zhu et al. | Wind Speed Short-Term Prediction Based on Empirical Wavelet Transform, Recurrent Neural Network and Error Correction | |
CN116777066B (zh) | 基于地基云图引导深度分解的光伏功率概率预测方法及装置 | |
CN117094431A (zh) | 一种多尺度熵门控的DWTformer气象数据时序预测方法及设备 | |
Ding et al. | A hybrid CEEMDAN-VMD-TimesNet model for significant wave height prediction in the South Sea of China |
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 |