基于数理统计的短期风速预测修正方法
2023-12-12
来源:乌哈旅游
研究与开发 基于数理统计的短期风速预测修正方法 王世谦 田春筝黄景慧 450052) (国网河南省电力公司经济技术研究院,郑州摘要准确的短期风速预测,可以有效地减轻风电场对电力系统的不利影响,提高风电场在 电力市场中的竞争力。由于风速具有很强的随机性和波动性,常规的风速预测方法在风速突变点 均存在较大的预测误差,且通过改善预测方法本身无法得到有效地解决。本文从预测数据后期修 正的角度出发,结合数理统计思想,提出了一种基于历史风速数据波动特性及一定置信水平的短 期风速预测修正方法,将原始风速预测值进行修正处理,得到优化风速预测值。该方法具有通用 性,可以应用于不同方法风速预测结果的修正。对基于灰色GM(I,1)模型的风速预测结果进行 修正,验证了该方法的有效性和可行性。 关键词:风速预测;波动特性;数理统计;数据修正;风力发电 Correction Method for Short Term Wind Speed Forecasting Results Based on Mathematical Statistics Wang Shiqian Tian Chunzheng Huang Jinghui (He’nan Electric Power Corporation Economic Research Institute,Zhengzhou 450052) Abstract To accurately forecast short—term wind speed can effectively reduce the adverse effects of wind power in power systems and improve the competition of wind power in power markets.Because of the strong randomicity of wind,the forecasting method at wind mutation points are greatly improved, and through improve the forecasting method itself cannot obtain more favorable effect.From the angle of revision,a revision method was proposed and applied to revision the forecast wind speed in this paper, which based on the historical data wave characteristics and confidence leve1.Under this method,we can turned the original forecast wind speed into the optimized forecast wind speed.The revision method applies to all short—term wind prediction methods.The validity and feasibility of the method was verified through the real example of a wind speed forecasting method based on the Grey mode1. Key words:wind speed forecast;wave characteristics:mathematics statistics;data revision; wind power 近年来,随着化石能源危机、环境污染和温室 效应等问题的日益加剧,风能已成为当今世界增长 最快的可再生能源【l】。但由于风电具有很强的随机 性和波动性,风电事业的发展受到了很大的限制。 对风电场和风速进行准确的预测,将有利于电力系 Analysis1[4-5]、人工神经网络法(ANN)[6-7]、模糊逻 辑法(Fuzzy Logic1[8-9】、空间相关性法(Spatial Correlation1[8-9]、小波分析法 以及上述方法的各种 组合[11-1 3]等,利用这些方法进行的风电场风速预测 的误差在25% ̄40%左右。预测误差不仅与预测方 统调度部门调整调度计划,减轻风电对电网的不利 影响,同时有效地减少电力系统的运行成本以及旋 转备用,并且有利于在电力市场环境下制定正确的 电能交换计划【2J。 法有关,还与预测周期以及预测地点的风速特性有 关。一般来说,预测周期越短,预测地点的风速变 化越缓和,预测误差就会越小;反之,预测误差就 会越大[I4]。 由于预测的滞后性,现有的风速预测方法均存 在着在风速突变点预测误差较大,进而降低整体预 目前,风速预测方法主要有,卡尔曼滤波法 (Kalman filter)[31、时间序列分析法(Time Series 1 0 l电囊|{技术 2013年第1 1期 I 测效果的情况,且通过改善预测方法本身并不能得 到有效地解决。同时,当前对风速预测的研究更多 的是从预测方法改进的角度出发,但预测效果并没 有随着预测方法的日益复杂而得到有效的提高,尤 其是风速突变点预测误差较大的情况迄今未得到有 效的改善。关于风速预测修正方面的研究尚存在较 大改进空间。 鉴于上述分析,本文从短期风速预测修正的角 度出发,提出了一种基于历史风速数据波动特性及 置信水平的风速预测修正方法。该方法先对风速历 史数据进行统计分析,根据不同的置信水平确定的 预测时刻的风速波动范围,之后对现有成熟的预测 方法获取的预测结果进行修正,从而提高整体预测 精度。 l 风速突变点预测问题分析 风速预测方法,通常是指一套系统的处理过去 与未来风速的数学方法,即在满足一定精度要求的 前提下,利用历史风速数据,尤其是与预测时刻相 关性最强的几个数据的内在规律特性确定未来某特 定时刻的风速值或某特定时段的平均风速值。 由持续预测法预测风速的原理可知,相邻两组 数据的相关性最强,且随着相隔时间的增加相关性 减弱,这就表明预测时刻的前两组数据的内在规律 变化对特定时刻的风速预测值影响最大,这也就意 味着当预测时刻的前两组风速数据呈现递增的规律 时,特定时刻的风速预测值呈现继续增加的趋势, 反之亦然,这也与实际情况相符。 风速突变点,通常是指风速变化规律在某一固定 监测时刻发生改变时对应的风速监测值,即风速突变 点的风速值与上一监测时刻的风速值的增减规律和 上一监测时刻风速值与上上监测时刻的风速值的增 减规律不同。如图1所示,图中第4个时刻的风速值为 风速突变点,且可以发现在风速突变点预测误差较大 的原因,即原有预测继续呈现增加的趋势,而实际风 速值出现减少,这就意味着风速突变点的实际值与预 测值之间的相对差值比一般时刻的相对差值要大。选 例中预测风速值由基于灰色模型GM(1,1)的风速预 测方法得到,原始风速序列长度为5。 对于这些突变点,通过进一步完善风速预测方 法,通常很难提高预测精度,这是由预测方法的滞后 性所决定的。因此,本文引入历史风速数据的波动特 性来确定未来特定时刻的风速波动范围,以此来弥补 研究与开发 预测方法的滞后性所引发的较大的预测误差。 11,q4 。。 = 篇麟 篓 ll 10.02 10.07 l0.19 —一预测德一实测值 2 3 4 预测时刻 图1风速突变点预测情况示意图 2 基于波动特性及置信水平的短期风速预 测修正模型 风电功率随着风速的随机变化而呈现很强的非 规律性的波动是风力发电最显著的特点。从长时间 来看,风速的变化幅度很大,且周期性不强,但风 速是逐渐变化的,即在秒级数据记录中,大部分时 间内风速的变化率是很平缓的。 由基于持续预测法的风速预测基本原理,可知 相邻监测的风速值v 与v 1相关性最大,这也与实际 经验相符。具体到风电机组而言,可以通过单位时 间风速变化率和风速波动率两个指标来表示其对应 的风速波动特性。定义如下: △ 二 : △ △ f= 二 : l00% (2) 式中,Av 、 分别为固定监测时间间隔对应的单 位时问风速变化率(单位:m/s )和风速波动率(单位: 1),f为监测点编号, 、Vf(f_f)为相邻监测风速值, △ 为监测时问间隔(单位:s)。 对历史风速数据进行上述处理,并对Av 、Ap 进行分段统计,可以得知不同单位时间风速变化率 和风速波动率出现的频率,进而可以得到基于一定 置信水平 的单位时间风速变化率Ay 和风速波动 率△ 最大值△ ^、△ ,即单位时间风速变 化率小于△ 一, 和风速波动率小于 , 的概率 值为 。 由单位时间风速变化率△V 、风速波动率 的最大值△ ^、△ i 州和监测风速值 (f_】),可 确定下一时刻风速Vi,的波动范围,即有 2013. ̄11期电毒囊l技徽I 1 1 研究与开发 [rI n(( (f—1)一△ n ×△ ,( o 小时共计144组的风速数据。 修正对象为基于灰色GM(1,1)模型的风速预测 1)二+_gi (t 1)× ))l 3 r ̄x((vi(t 1)-{-AVimax,^×△ ,(Vi(t 简记为 [ ffmin, ,V itm“。^] (4) 此时可建立基于一定置信水平 的风速预测初 步修正公式 f “m ≤ …, :{【 ,V ii一, m≥ in,2≤V一,t it≤ 。一. (5) 式中, ,,为初步修正预测风速, 为原始预测风 速, , 、 i 。 分别为基于历史风速数据及一 定置信水平 确定的f时刻的风速波动范围的上、 下限。 由基于一定置信水平 确定的初步修正预测风 速仍需满足由历史风速数据确定的整体风速的波动 范围[vimin,vim ],其中vim 可设定为0m/s,Vim 可 设定为风电机组的切出风速。 最后,可建立完善的风速预测修正公式 f ≤ f ,:I ,, Vfmi ≤ ≤V max (6) lVimax ≥Vimax 式中,V一 为修正预测风速, 为初步修正预测风 速,Vi 、Vi i 分别为基于历史风速数据确定的整 体风速波动范围的上、下限。 分析可知,当 =0时,原有预测方法将变更为 持续预测法的一种情况,即把最近一点的风速观测 值作为下一点的风速预测值,此时预测误差较大, 这也说明置信水平不能设置过低;当 =1时,修正 预测风速等于原始预测风速,修正方法将不起作用, 修正效果得不到体现。 3算例分析 通过上节内容可知,基于历史风速数据的波动 特性及一定置信水平建立的修正模型对风速预测方 法并没有特殊要求,限于篇幅,本文仅介绍基于GM (1,1)灰色模型的短期风速预测方法,并用实际算 例验证该方法的有效性与可行性。 3.1研究数据及修正对象 历史数据采用我国北方某风电场在2009年3月7 日至17日以10分钟为采样间隔的实际监测风速数据 (共计l584组)。预测日期及时段为3月18日持续24 12 l电 技】lt 2ol3 ̄g 11期 值,鉴于篇幅有限,GM(1,1)模型可参见文献[151, 数据窗长度为5,对3月18日共计144组的风速的预测 数据进行修正。 3.2修正步骤 第一步:由历史风速数据得到单位时间风速变 化率和风速波动率,并完成分段统计。结果见表1、 表2。 第二步:设定修正置信水平 ,进而确定与之 对应的单位时间风速变化率、风速波动率的最大值 Aw ril'ax. 、Ap max. 。 如设定置信水平 =0.85,由历史风速数据统计 结果可知对应的单位时间风速变化率△ 嗽 等于 2.13 10。m/s ,以及对应不同风速段的风速波动率 的最大值△p “,如Vi( )位于9~10m/s时, / ̄pimax, =13.5%,具体结果见表3(a)、(b)。 表1 基于10rain监测间隔的单位时间风速变化率统计数据 单位时间风速变化率(1O。m/s 1 概率/% 0.O~0.4 28.17 0.4~0.8 20.66 0.8~1.2 15.4l 1.2~1.6 12.00 1.6~2.0 7.O1 2.O~3.O 9.41 3.O~6.0 6_32 >6.0 1.O1 表2基于10rain监测间隔的风速波动率统计数据 风速 风速波动率概率/% /(m/s) 0—10 l0—20 20 40 40-8O 80—15O >l50 O~l 1-39 2.78 5.56 15.98 63.20 11.11 1~2 17.78 12.22 27.78 25.55 12.22 4.44 2~3 27.68 17.86 25.90 17.85 lO.72 0.00 3~4 28.7l 26.73 31.68 10.89 1.99 0.O0 4~5 40.16 27.05 22.96 8_2O 1 64 0.00 5~6 49.11 25.0O l8.75 6-25 0 9 0.O0 6~7 63.O3 25 21 l0.O8 1.68 O.00 0.O0 7~8 75.5O 17-22 5.96 l-32 O.OO 0.O0 8~9 73.64 2O.16 5.43 O.78 O.00 O.O0 9~10 73.11 23.53 3 36 0.00 O.00 0.O0 l0~11 76.19 19.05 3.57 1.19 0.O0 0.0O 11~12 83 l3 15.66 1.2O O.O0 O.O0 0.OO l2~13 78_33 20.O0 1.67 0.00 0.00 0.00 l3~14 88.52 11.48 O.00 0.O0 O.OO 0.OO I4~l5 80.00 20.00 0.O0 0.OO 0.O0 0.00 15~l6 89.47 5.26 5.2632 O.O0 O.O0 0.0O >l6 1O0 0.O0 0.0O 0.00 0.00 0.0O 第三步:基于修正公式(4)、 (5)完成预测风 速数据修正,得到优化预测风速值。 如置信水平 =0.85时,修正前、后预测风速 与实际监测风速对比结果如图2所示。由图2可知, 通过修正处理,优化预测风速值较原始预测风速值 更加接近于实际值,有效地降低了预测风速值与实 际风速值的偏离程度,弥补了预测方法的滞后性引 起的较大误差,尤其是在风速突变点修正效果得到 明显的体现,提高了整体预测效果。 鬟 戆 矗 l 逢 曼最 一 多厂 i 一 — 实际监测风速 一修正前预测风速 修正后预测风速 ” ∞ 61时 0 85 如0挖 坞 图2预测风速修正前后与实际监测风速对比图 3.3修正误差分析 修正误差的确定采用绝对百分比误差 PE和平 均绝对百分比误差MAPE。定义如下: P =L—— : 二 ×1O0% × % (7)( ) MAPE:一1 APE (8) N 式中,Ⅳ是参与计算的样本点数。最大、最小绝对 百分比误差分别记为MAx APE、MIN APE。 为了判断不同置信水平对修正效果的影响,本 文基于不同的置信水平对修正对象进行了多组修 正,其中不同置信水平对应的单位时间风速变化率、 风速波动率的最大值Avimax, 、Ap,一, 见表3、表4, 修正结果对比见表5、图3。修正对象原始预测风速 的MAXAPE(%)、MINAPE(%)、MAPE(%)分 别为60.35%、0.1l%、14.44%。 由图2、表5可以看出,经过本文方法的数据修正, 预测精度较修正前有了明显的提高。与修正对象的原始 预测数据相比,其平均绝对百分比误差、最大绝对百分 比误差均有了较大幅度的下降,整体预测效果得到明显 提高。如置信度水平 =0.85时,其平均绝对百分比误 差降低为12.29%,预测精度提高了2.15%,最大绝对 百分比误差也下从60.35%降低到了49.02%,有效的缓 解了风速突变点预测误差较大的情况。 研究与开发 表3不同置信水平及其对应Av 一, 置信度 O.65 0.7 0.75 0.8 0.85 0.9 0.95 Av/max,2 1.23 l_37 1.55 1.8 2.13 2.60 3.5O 表4不同置信水平及其对应 , 风速段 置信水平 m∞ O.65 0.70 O.75 0.80 O.85 0.9 0.95 0~1 lO0 100 10O 1O0 1O0 1 74 298. 1~2 48.2 54.3 63.9 68.0 80.2 95.3 112 2~3 32.8 37.1 44-3 51.4 57.9 84.0 96.6 3~4 22-3 24.3 27.1 28.9 35.6 43.6 60.5 4~5 l9.0 21.O 22.O 24.2 28.5 37.8 5O.5 5~6 15.1 16.5 20.O 23.1 25.0 31.9 42.7 6~7 l0.2 11.2 12.8 l3.6 16.8 20.6 33.6 7~8 7.73 8.76 9.94 1O.9 13.2 15.9 24.2 8~9 7.58 8.57 l0_3 11_3 l3.1 15.7 21.4 9~10 7.22 8.65 l0.1 12.2 13.5 15-3 l7.6 10~l1 6.83 7.55 9.52 10.9 13.5 14.6 19.5 11~l2 5-36 6.17 6.54 7.90 11.2 13.7 15.3 12~l3 6.14 7.13 9.10 10.O l1_3 13.4 19.5 l3~l4 4.35 4.77 6_35 7.15 7.84 10.1 l2.9 l4~15 5_32 5-39 6.64 6.88 10.1 l1.4 17.9 l5~l6 5.83 6.76 7.0l 7.12 8_36 8.69 l6.4 >16 4.84 4.84 6.46 6.46 7.42 7.42 12.1 十平均绝对百分比误差 啪最_丈绝对百分比误差 一/ 一一 ~一 一0.5 0.,55 0.6 0.65 0 7 0 75 0.8 0 8,5 0.9 0 95 置信度水平 图3不同置信水平修正结果比较 表5不同置信水平修正结果 置信度2 MAPE|% M X4PE| MlNAPE|% 0.5O 8.13 38.94 O.O3 O.55 8.61 40-25 0.04 0.6O 8.95 41.03 O.O3 0.65 9.52 41.79 O.05 0.70 l0.03 43.15 O.15 O.75 10.70 44.71 0.00 0.80 l1.27 45.93 0.42 0.85 1’,0 49.O2 O.01 0.90 13.18 52.5O 0_24 0.95 14-20 60-35 O.11 同时,由图3、表5可以看出,随着置信水平 的不断降低,修正对象的平均绝对百分比误差、最 大绝对百分比误差呈现逐步下降的趋势,尤其是当 置信度水平 =0.50时,其平均绝对百分比误差、 最大绝对百分比误差两项指标下降幅度分别高达 I2013年第1 1期电 技 I 13 研究与开发 42.75%、35.45%,预测效果得到根本性的提升, 但 此时置信水平较低,可信度较差。 总体而言,通过历史风速数据的统计处理,可 以较好的发现风电机组对应的特定地区的风速的变 化规律,有效的改善风速预测模型在风速突变点预 测误差较大的情形,提高整体风速预测效果,为准 确的风电功率预测奠定了良好的基础,对减轻风电 场对电力系统的不利影响,提高风电场在电力市场 中的竞争能力具有重要的意义。但是,仍然存在个 别监测点预测误差偏大的现象,仍需进一步的研究, 如 =0.85时,其最大绝对百分比误差仍然达到 49.O2%。 最后,基于上述分析可知,获取准确的基于一 定置信水平的单位时间风速变化率、风速波动率的 最大值Avi , 、Api , 对修正效果起着至关重 要的作用,此时也可以获知不断增加历史风速数据 的总量将逐步提高风速预测修正效果,进而提高风 速预测模型的整体预测效果。 4结论 本文从预测风速预测结果修正的角度出发,提 出了一种基于波动特性及置信水平的风速预测修正 方法,即通过分析大量历史风速数据的波动特性建 立具有一定置信水平的预测时刻的风速波动范围, 利用风速波动范围将原始风速预测值进行修正处 理,得到了最终的优化风速预测值。对我国某风电 场基于灰色模型GM(1,1)的风速预测结果进行了修 正处理,并讨论了不同置信度对预测结果的影响。 结果表明,经过修正处理可以得到更接近于实 际值的优化风速预测值,有效地缓解了因预测方法 的滞后性引起的风速突变点预测误差较大的情况, 且随着置信水平的不断减小,整体的预测效果呈现 逐步提高的趋势。该修正模型可有效地提高整体风 速预测效果,风速突变点的预测精度提高得尤为明 显,对提高风电在电力市场中的竞争力具有重要意 义。该方法与预测方法无关,理论上可以应用于不 同方法风速预测结果的修正。 参考文献 [1]杨茂,熊吴,严干贵,等.基于数据挖掘和模糊聚类的风 电功率实时预测研究[J].电力系统保护与控制,2013, 41(1):1-6. 14 I电 技琳 2013年第11期 [2】谷国利,王维庆,张新燕,等.风电场风速预测方法的研 究【J】.农业工程学报:新能源产业,2009(6):22—24. [3] 卿湘运,杨富文,王行愚.采用贝叶斯一克里金一卡尔曼 模型的多风电场风速短期预测[J].中国电机工程学 报,2012,32(35):107—1 14. [4】TORRES J L,GARCIA A,BLAS M D,et a 1.Forecast of hourly average wind speed with ARMA models in navarre(Spain)[J].Solar Energy,2005,79:65—77. [5】张宏宇,印永华,申洪,等.基于概率测度变换的风速时间 序列建模方法[J].电力系统自动化,2013,37(2):7—10. 【6]BARBOUNIS T GTHEOCHARIS J B,ALEXIADIS M C.et a1.Long—term wind speed and power forecasting using local recurrent neural network Models[J].Transactions on energy conversion,2006, 21(1):273—284. [7】黄小华,李德源,吕文阁.基于人工神经网络模型的风 速预测[J]_太阳能学报,2011,32(2):193-197. [8]ALEXIADIS M,DOKOPOULOS SAHSAMANOGLOU H.et a1.Short term forecasting of wind speed and related electrical power[J].Solar Energy,1998,63(1): 61—68. [9】 BARBOUNIS T G,THEOCHARIS J B.A locally recurrent fuzzy neural network with application to the wind speed prediction using spatial correlation[J]. Neurocomputing,2007,70(1 5):25-42. [10]郭虎全,刘吉臻,柳玉,等.基于小波包分析的风速预测 研究[JJ.华东电力,2011,39(12):2077—2079 [1 1]罗文,王莉娜.基于小波分解与遗传算法和支持向量机 的短期风速预NIJ[j].太阳能学报,2012,33(8):1327—1333. [12]刘辉,田红旗,李燕飞.基于小波分析法与滚动式时间 序列法的风电场风速短期预测优化算法[J].中南大 学学报(自然科学版),2010,41(1):370—375. [13】蔡凯,谭伦农,等.时间序列与神经网络法相结合的短 期风速预测[J].电网技术,2008,32(8):82—85. [14】杨秀嫒,肖洋,陈树勇.风电场风速和发电功率预测研 究[J].中国电机工程学报,2005,25(11):1-5. [15】杜松怀,温步瀛,蒋传文.电力市场[M].3版.北京:中国 电力出版社.2008:107.108. 作者简介 王世谦(1988.),男,河南濮阳人,硕士,主要研究方向为电力系 统分析、电网规划。