1、基于时间序列的风速和风功率预测 孙睿男 王飞 江苏师范大学 摘 要： 短期的风速和风功率预测对风场的运行及并网发电是非常重要的。采用时间序列方法, 通过 Matlab 编程, 建立了 ARIMA 和 SARIMA 预测模型对 24 小时后的平均风速和风功率进行预测。最后, 对两个预测模型的结果作了对比分析。结果表明 SARIMA 模型的预测结果是好于 ARIMA 模型。关键词： 时间序列; 风速; 风功率; 预测; ARIMA 模型; SARIMA 模型; 作者简介：孙睿男 (1975) , 男, 江苏徐州人, 高级工, 主要从事电工技术实验教学;作者简介：王飞 (1976) , 男, 内蒙古
2、赤峰人, 实验师, 硕士, 研究方向为可再生能源发电和经济负荷调度。收稿日期：2017-09-15基金：江苏师范大学实验室建设与管理研究项目 (L2017Y11) Forecast of Wind Speed and Wing Power Based on Time SeriesSUN Rui-nan WANG Fei Jiangsu Normal University; Abstract： Short-term forecasting of wind speed and wind power is great importance to wind farm operation and gri
3、d-generation. Based on the time series method and MATLAB software, in this article we have used the ARIMA (Autoregressive Integrated Moving Average) and SARIMA (Seasonal Autoregressive Integrated Moving Average) models to forecast the hourly average wind speed and wind power up to 24 h in advance. F
4、inally, a pair of comparison analysis Based on ARIMA and SARIMA model is run with average relative error. Results indicate that the SARIMA models significantly better in forecast than the ARIMA model.Keyword： time series; wind speed; wind power; forecast; ARIMA model; SARIMA model; Received： 2017-09
5、-151 概述近些年来, 风力发电一直都是世界上发展最快的新能源之一。随着风力发电的大力发展, 风能的随机性、季节性、波动性以及不确定性使得当风力发电大规模的并入电网时, 对于电力系统的稳定性波动影响较大, 给电网的调度以及规划带来了困难1, 为了减少风电的功率储备, 方便电力调度, 降低成本提高电网的风能质量, 便于风电场安排发电机组的检修、保养, 便要求对风速及发电功率有比较准确的预测研究。目前文献中, 相关的研究者使用不同的技术方法对风速和风功率预测。通常这些技术能被分成以下四种:数值天气预报, 统计方法, 人工智能和机器学习方法, 时空模型。这些预测法方法各有特点和适用的场合, 数值天
6、气预报方法对于短期的风速和风功率预测不能得到一个可靠的结果, 但是它在长期预测上有着很好地预测结果。统计方法, 其中最著名的就是自回归滑动平均 (ARMA) 模型, 这种方法能有效地考虑时间序列的有效性及随机波动的干扰性, 从而能准确地预测事件的短期变化趋势。人工智能和机器学习方法预测效果很好, 但是不容易确定网络结构, 学习速度慢。时空模型要求原始数据量大, 计算不便。时间序列模型最为简单, 也是花费最小的预测模型2-5。考虑到风速和风功率的时间序列通常具有趋势性、季节性以及非平稳性, 本文采用差分自回归滑动平均模型 (ARIMA) 和带有季节因子的差分自回归滑动平均模型 (SARIMA)
7、对风电场的风速和风功率进行短期预测。2 时间序列模型概述时间序列是指按照时间顺序来排列的、随时间而变化并且相互关联的数据序列。分析时间序列的方法即时间序列分析, 时间序列分析最早是起源于 1927 年, 数学家耶尔 (Yule) 为了研究预测市场变化规律而建立的自回归 (AR) 模型, 之后在 1931 年, 瓦尔格 (Walker) 在 AR 模型启发下又建立了滑动平均 (MA) 模型以及自回归滑动平均 (ARMA) 模型6, 从而奠定了时间序列分析的基础。按照时间序列统计的特性来分类, 可以分为平稳时间序列和非平稳时间序列。2.1 平稳时间序列模型最常见的平稳时间序列模型是自回归移动平均模
8、型 (AR-MA) , 设X t, t=0, 1, 2, 是零均值平稳序列, 满足下列模型:其中 t是零均值、方差是 的平稳白噪声。则称 Xt是阶数为 p, q 的自回归滑动平均序列, 简记为 ARMA (p, q) 序列, 当 q=0 时, 它是自回归模型 AR (p) 序列;当 p=0 时, 它变成滑动平均模型 MA (q) 序列。= ( 1, 2, , p) 称为自回归参数向量, 其分量 j, j=1, 2, , p称为自回归系数, 它决定前一时刻时间序列的值多大程度上影响当前时刻的值。= ( 1, 2, , Q) 称为滑动平均参数向量, 其分量 j, j=1, 2, , q称为滑动平均
9、系数, 它决定前一时刻高斯随机变量的值影响现在值的程度, 这些参数都需要识别。整理的 ARMA 模型的一般形式如下:2.2 非平稳时间序列模型ARMA 是对平稳时间序列分析的方法, 但在现实生活中的时间序列很多都是非平稳的序列, 它们往往带有趋势性、季节性以及非平稳性, 在研究这类时间序列时要先通过差分变换将其转化为平稳时间序列。2.2.1 自回归差分移动平均 (ARIMA) 模型引入差分算子=1-B, 对原有非平稳时间序列 Xt进行一阶差分变换得到原时间序列进行 d 阶差分后得到:原时间序列用 ARMA 模型来描述可表示为:由上式可以看出, W tWt是一个平稳的 ARMA 序列, 且它是非
10、平稳序列 XtXt的 d 阶差分。由此, 一个非平稳时间序列就转换成了一个平稳的时间序列了, 并且把这个 Xt序列称为 ARI-MA (p, d, q) 序列。2.2.2 季节性自回归差分移动平均 (SARIMA) 模型如果遇到的非平稳时间序列有很明显的季节性或由季节性因素引起的周期性的变化, 我们分析此类模型时可采用季节性自回归差分移动平均 (SARIMA) 模型。一个时间序列X t, t=0, 1, 2, , s 为一个正整数, 满足下列模型:其中 S为季节差分算子 , 它定义为 s=1-B, S表示间隔为 S 步的 D 阶差分。Et一般不必是白噪声, 而可设它是另一个 ARIMA (p,
11、 d, q) 序列综合式 (6) 和 (7) 可以得到季节性自回归差分移动平均 (SARIMA) 模型表达形式:3 基于时序分析的风速和风功率预测模型的建立3.1 AMIMA 模型的建立ARIMA 的建模过程通常分为四步:时间序列的平稳性检验及平稳化、模型识别、参数估计以及建立预测模型。本文使用我国某风电场的一组实测风速和风功率数据进行分析, 时间序列为间隔一小时, 取此序列前 260 点数据为原始数据建立预测模型, 并且用所建立的模型预测后 24 个小时的数据。3.1.1 原始时间序列的平稳性及平稳化判断原风速和风功率序列是否为平稳的序列时, 可以通过观察它的自相关函数 (AC) 是否呈指数
12、迅速衰减到零, 如果是则是平稳的序列, 否则是非平稳。图 1为前 260 个样本值及其前 20 个自相关函数以及偏相关函数 (PAC) 值。图 1 风速序列 下载原图从图 1 可以直观地看出自相关函数不能够迅速的衰减至零, 原始风速和风功率序列是非平稳的, 需要进行平稳化处理。对于非平稳的时间序列, 本文进行差分运算以消除它的趋势性, 进行一阶差分运算之后, 得到两个新的差分序列如图 2, 计算这两个一阶差分序列的自相关函数值和偏相关函数值。图 2 一阶差分运算后风速序列 下载原图从图可以看出原始的风速和风功率序列经一阶差分后序列变得基本平稳, 可以确定相应的模型和阶数。3.1.2 确定 AR
13、IMA 模型及阶数本文利用最小信息准则 (AIC) 准则7来进行参数拟合, 选择最适合的阶数, 由于运算较大, 文中采用 Matlab 进行编程做参数拟合。式中, L 为对数似然值;n 为观测值的数目。ARIMA 模型定阶的过程就是选择适合的 p、q 使得 AIC 的值最小, 这样模型就越接近实际模型, 误差越小。根据AIC 准则, 从程序的运行结果来看, 当 R=1, M=3 时, AIC 取得最小值, 最符合实际模型。所以选用 ARIMA (1, 1, 3) 模型, 进行风速预测。选择 R=1, M=0时, AIC 的值最小, 则建立 ARIMA (1, 1, 0) 的预测模型对风电功率进
14、行预测。3.2 SARIMA 模型的建立在 ARIMA 建模的数据分析中已经知道, 原始风速和风功率时间序列是一个非平稳的时间序列, 那么首先要进行平稳化的处理, 由于采样点是时间间隔为一小时, 所以选取周期为一天 24 个小时, 即 S=24。首先对原风速序列进行季节差分处理, 然后进行取消趋势性的一阶差分运算。3.2.1 原始时间序列的预处理在确定了序列的周期 S 后, 先对风速和风功率时间序列进行趋势差分和季节性差分, 图 3 得到一个平稳的新的时间序列。图 3 原始风速序列 下载原图风速和风功率经过季节性差分和趋势性差分的计算后得到的新序列的自相关函数以及偏相关函数值如图 4 所示。从
15、图形 4 中可以看出, 经过季节差分和趋势性差分后, 风速和风功率得到的新序列已经基本平稳了。3.2.2 确定 SARIMA 模型及阶数依旧根据 AIC 准则并采用 Matlab 进行编程对模型进行定阶运算。从程序的运行结果可以看出, 选择 R=3, M=25 时, AIC 的值最小, 所以风速序列选取 SARIMA预测模型阶数为 R=3, M=25。从 Matlab 的运行结果可以看出取 R=1, M=12 时, AIC 的取值最小, 建立的风功率模型能够最接近实际模型。图 4 风速新序列 下载原图3.3 结果对比分析3.3.1 风速预测结果分析使用前文确定的模型进行风速和风功率预测, 将
16、ARIMA 模型的风速预测结果以及 SARIMA 模型的风速预测结果进行对比分析, 图 5 为实际风速以及 ARIMA 预测风速和 SARIMA 预测风速的趋势对比图。由图 5 的预测结果对比可以看出, ARIMA 模型的预测结果符合实际序列的走势, 但转折没有, 所以预测效果不太好, 而 SARIMA 模型的预测结果存在转折性, 但和实际转折点有误差, 但总体来说预测效果较好。从 24 个预测值误差分析可以看出, SARIMA 模型的预测值的平均相对误差为 8.55%, 而 ARI-MA 模型的预测值平均误差为 9.94%, 所以 SARIMA 模型预测的精确度更高, 更加符合实际的模型,
17、这和图 5 直观的预测结果对比图分析出的结果是一致的。3.3.2 风功率预测结果分析分析比较 ARIMA 模型和 SARIMA 模型的预测结果, 连续 24 小时的实际风电功率数据和 ARIMA 模型风电功率预测数据及 SARIMA 模型风电功率预测数据的趋势对比如图 6 所示。从图 6 已经可以直观地看出 SARIMA 模型的预测效果要比 ARIMA 模型的预测效果好, 从 24 个预测值的对比 SARIMA 模型的预测值的平均相对误差为 17.3%, 而ARIMA 模型的预测值平均误差为 74.8%, 结果可以很明显看出 SARIMA 模型的预测结果精度要大大超过 ARIMA 模型的预测结
18、果。4 结束语1) 通过对实际风电场数据序列的分析, 风速和风功率序列具有时序性, 自相关性及非平稳性, 适合建立相应的时间序列模型。2) 本文采用差分的方法, 根据风电场实际数据建立 ARI-MA 和 SARIMA 模型, 并对两个模型的预测结果进行对比分析, SARIMA 模型比 ARIMA 模型预测的结果好, 更能反映序列的变化, 具有一定的指导意义。图 5 风速预测结果对比分析图 下载原图图 6 风功率预测结果对比图 下载原图3) 为了更好地预测风速和风功率这种非平稳的时间序列, 可以从数据的预处理方法和根据模型的优缺点改进现有模型两个方面研究。参考文献1薛禹胜, 郁琛, 赵俊华.关于
19、短期及超短期风电功率预测的评述J.电力系统自动化, 2015, 39 (6) :141-150. 2潘迪夫, 刘辉, 李燕飞.基于时间序列分析和卡尔曼滤波算法的风电场风速预测优化模型J.电网技术, 2008, 32 (7) :82-86. 3高山, 梁伟平.基于时间序列的风电场短期风速预测J, 电气时代, 2014 (11) :76-78. 4杨秀媛, 肖洋, 陈树勇.风电场风速和发电功率预测研究J, 中国电机工程学报, 2005, 25 (11) :1-5. 5白春艳, 王飞.基于三段式变步长电导增量法 MPPT 技术研究J.电脑知识与技术, 2013, 13, 3155-3158. 6 (英) 恩德斯等著, 杜江, 谢志超译.时间序列分析M.北京:高等教育出版社, 2006.6 7常太华, 王璐, 马巍.基于 AR、ARIMA 模型的风速预测J.华东电力, 2010, 38 (1) :59-62.
