一种基于主成分分析的Lorenz扰动风速预测方法技术

本发明专利技术公开了一种基于主成分分析的Lorenz扰动风速预测方法，属于发电技术领域。本发明专利技术是根据原始数据进行主成分分析，再利用最小二乘支持向量机(LS‑SVM)模型对主成分和风速数据进行风速预测，最后利用Lorenz大气扰动序列对风速的初步预测值进行修正，以提高风速的预测水平。仿真结果显示，在风速预测过程中充分考虑大气系统的扰动影响，显著提高了风速的预测精度。本发明专利技术弥补了风速预测领域对大气系统的忽视作用，有助于风电并网的稳定性和风资源的大规模开发利用。

A method of Lorenz disturbance wind speed prediction based on principal component analysis

The invention discloses a Lorenz wind speed forecasting method based on principal component analysis, which belongs to the field of power generation technology. The present invention is a principal component analysis based on the original data, and then using the least squares support vector machine (LS SVM) model of wind speed on the principal components and wind speed data prediction, finally using the Lorenz atmospheric disturbance was corrected preliminary prediction on wind speed sequence, in order to improve the level of wind speed forecast. The simulation results show that in the process of wind speed prediction, the disturbance effect of the atmospheric system is taken into account, and the prediction accuracy of wind speed is significantly improved. The invention compensates the neglect effect of the wind speed prediction field on the atmosphere system, and is helpful for the stability of the wind power grid connection and the large-scale exploitation and utilization of the wind resources.

【技术实现步骤摘要】
一种基于主成分分析的Lorenz扰动风速预测方法
本专利技术涉及一种能够准确预测风电场风速的方法，属于发电

技术介绍
当今世界，化石燃料的供给紧张带来了环境、生态和全球气候变化等一系列问题，各国都自觉行动起来，主动破解困局，加快能源转型，发展可再生能源。风能作为一种环保、可再生的清洁能源，在低碳能源技术中占有举足轻重的地位，受到了全球的广泛关注。同时，随着风力发电在电力生产中的份额越来越大，电网运营商在电力平衡、电能质量、并网的稳定性、负荷调度计划等方面面临着挑战。因此，在风力发电的并网过程中，必须考虑可靠的风电功率预测方法。国内外研究学者提出了很多关于风电功率预测的研究方法。目前风电功率预测按照预测周期的长短分为长期预测、中期预测、短期预测和超短期预测；按照建立的数学模型可分为物理模型、统计模型、人工智能技术和组合模型。鉴于预测方法的多样性，国内外出现了许多成熟稳定的预报系统。国外从事风力发电功率预测研究工作较早，最具代表性的预报工作包括Prediktor、WPPT、Previento、eWind等等。国内开发的风电功率预报系统主要有WINPOP系统，WPPS系统，WPFS系统、T213L31系统等。以上的风电功率预测方法大多是基于算法模型以及预报系统的改进研究，而忽视了对于大气动力系统的扰动影响研究，因而预测效果不是很理想。
技术实现思路
本专利技术所要解决的技术问题是提供一种基于主成分分析的Lorenz扰动风速预测方法，以达到提高风速预测的目的。为解决上述技术问题，本专利技术所采取的技术方案是：一种基于主成分分析的Lorenz扰动风速预测方法，其特征在于：(1)对待测风电场每隔一段时长t'采集1组原始数据，在时间t-1时刻采集到n组原始数据，所述原始数据包括风速数据(v)＝(vt-1-(n-1)t',vt-1-(n-2)t',......,vt-1)以及对风速的影响因素原始数据X1,X2,......,Xm分别表示对风速的影响因素集，m表示对风速的影响因素个数；(2)将t-1时刻所采集到的n组原始数据中的对风速的影响因素原始数据(X)进行主成分分析，根据累计贡献率确定主成分数据；(3)将主成分数据与将t-1时刻所所采集到的n组原始数据中的风速数据(v)结合作为最小二乘支持向量机(LS-SVM)模型的输入，将t时刻的风速vt作为模型的输出对数据进行训练，最后得到风速的初步预测序列v'(t)；(4)给定初值条件和参数，求解Lorenz方程，得到三维的大气扰动序列，并对扰动序列做标准化处理，得到无数量级和量纲的标准化扰动序列；(5)计算标准化数据与初始值之间的二阶闵可夫斯基(Minkowski)距离，得到Lorenz综合扰动流；(6)利用Lorenz综合扰动流对风速的初步预测序列v'(t)进行扰动修正，并以实现风速实际值和预测值之间的平均绝对误差最小为目标得到扰动修正的最优扰动强度和最优扰动系数；(7)利用Lorenz综合扰动流得到的最优扰动强度和最优扰动系数对风速的初步预测序列进行扰动修正，得到风速的扰动修正序列v"(t)。进一步的技术方案在于，所述对风速的影响因素包括风向D、气温T、气压P、比容a、比湿H和地表粗糙度R；进一步的技术方案在于，所述主成分分析的过程如下：a、以风向D的正弦值、风向D的余弦值、气温T、气压P、比容a、比湿H和地表粗糙度R为对风速的影响因素；将采集的数据构建主成分分析的原始数据，如下所示：X1,X2,......,Xm分别表示风向D的正弦值、风向D的余弦值、气温T、气压P、比容a、比湿H和地表粗糙度R的因素集，这里m＝7；b、对7个因素集的n×7个原始数据(X)标准化处理，以消除变量间在数量级上的差异产生的影响，使每个变量的均值为0，方差为1，得到标准化矩阵Y＝{ynm}；其中，xm和Sm分别表示Xm的均值和方差。c、根据标准化矩阵得出相关系数矩阵A；d、计算相关系数矩阵A的特征值及特征向量，并计算特征值贡献率和累计贡献率；其中，λm表示A的特征值。e、根据累计贡献率求主成分Z＝[Z1,Z2,...,Zm]Zm＝Yam其中，am表示λm的特征向量；m≤7。进一步的技术方案在于，所述步骤(3)将主成分数据与将t-1时刻所所采集到的n组原始数据中的风速数据(v)结合作为最小二乘支持向量机模型的输入矩阵M，将t时刻的风速vt作为模型的输出矩阵N对数据进行训练，最后得到风速的初步预测序列v'(t)；其中，式中，dmn表示主成分得到的数据；进一步的技术方案在于，所述Lorenz方程如下：其中，式中x,y,z是状态变量，分别表示对流运动的强度、对流运动流体的水平方向温差和垂直温差对无对流时的偏离程度。σ,r,b都是无量纲的正实数，分别表示普朗特数、瑞利数和与气候区域大小相关的量。对于经典的Lorenz系统的参数值σ＝10,给定初值h＝(1.1,1,1)，令r＝23,r＝27，我们会得到两组不同的数值解，进而得到两组不同的Lorenz吸引子。由于Lorenz方程的数值解是一组三维的扰动序列，因而先对方程的解进行标准化处理，得到无数量级和量纲的标准化扰动序列具体标准化处理公式为：式中xn,yn,zn表示Lorenz方程的数值解，xmin,ymin,zmin和xmax,ymax,zmax分别表示x,y,z的最小值和最大值。进一步的技术方案在于，所述步骤(5)中所述二阶闵可夫斯基(Minkowski)距离如下：计算标准化数据与初始值(0,0,0)之间的二阶闵可夫斯基距离：其中，Cn表示C0表示C0(x0,y0,z0)＝C0(0,0,0)。进一步的技术方案在于，采用平均绝对误差MAE和均方根误差RMSE来定量的评价模型的优劣，误差值越小，表示模型的预测能力越好。其公式为其中，y(t)和v(t)分别表示风速的实际值和预测值，k为预测样本个数。采用上述技术方案所产生的有益效果在于：在风速预测过程中引入描述大气运动状态的Lorenz系统对风速预测结果进行扰动修正，并利用风电场的实际数据进行仿真模拟。实验结果表明，经过Lorenz扰动修正后的模型显著提高了风电场的预测水平，有助于风能资源的大规模开发和利用。附图说明下面结合附图和具体实施方式对本专利技术作进一步详细的说明。图1(a)和(b)分别为风电场每隔十分钟的风速分布情况和其风向玫瑰图；图2(a)和(b)分布为瑞利数时的Lorenz吸引子形态；图3是在欧式距离下Lorenz系统生成的一组Lorenz综合扰动流；图4为时的Lorenz扰动强度，横坐标表示扰动的数据个数，即预测样本数据，纵坐标表示扰动强度的大小；图5是各个模型的风速预测曲线；图6是本专利技术的流程图。具体实施方式下面结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本专利技术的一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。在下面的描述中阐述了很多具体细节以便于充分理解本专利技术，但是本专利技术还可以采用其他不同于在此描述的其它方式来实施，本领域技术人员可以在不违背本专利技术内涵的情况下做类似推广，因此本专利技术不受下面公开的具体实施例的限制。本专利技术提供了一种本文档来自技高网...

【技术保护点】
一种基于主成分分析的Lorenz扰动风速预测方法，其特征在于：(1)对待测风电场每隔一段时长t'采集1组原始数据，在时间t‑1时刻采集到n组原始数据，所述原始数据包括风速数据(v)＝(vt‑1‑(n‑1)t',vt‑1‑(n‑2)t',......,vt‑1)以及对风速的影响因素原始数据

【技术特征摘要】
1.一种基于主成分分析的Lorenz扰动风速预测方法，其特征在于：(1)对待测风电场每隔一段时长t'采集1组原始数据，在时间t-1时刻采集到n组原始数据，所述原始数据包括风速数据(v)＝(vt-1-(n-1)t',vt-1-(n-2)t',......,vt-1)以及对风速的影响因素原始数据X1,X2,......,Xm分别表示对风速的影响因素集，m表示对风速的影响因素个数；(2)将t-1时刻所所采集到的n组原始数据中的对风速的影响因素原始数据(X)作为主成分分析，根据累计贡献率确定主成分数据；(3)将主成分数据与将t-1时刻所所采集到的n组原始数据中的风速数据(v)结合作为最小二乘支持向量机模型的输入，将t时刻的风速vt作为模型的输出对数据进行训练，最后得到风速的初步预测序列v'(t)；(4)给定初值条件和参数，求解Lorenz方程，得到三维的大气扰动序列，并对扰动序列做标准化处理，得到无数量级和量纲的标准化扰动序列；(5)计算标准化数据与初始值之间的二阶闵可夫斯基距离，得到Lorenz综合扰动流；(6)利用Lorenz综合扰动流对风速的初步预测序列v'(t)进行扰动修正，并以实现风速实际值和预测值之间的平均绝对误差最小为目标得到扰动修正的最优扰动强度和最优扰动系数；(7)利用Lorenz综合扰动流得到的最优扰动强度和最优扰动系数对风速的初步预测序列进行扰动修正，得到风速的扰动修正序列v"(t)。2.根据权利要求1所述得一种基于主成分分析的Lorenz扰动风速预测方法，其特征在于：所述对风速的影响因素包括风向D、气温T、气压P、比容a、比湿H和地表粗糙度R。3.根据权利要求1所述得一种基于主成分分析的Lorenz扰动风速预测方法，其特征在于：所述主成分分析的过程如下：a、以风向D的正弦值、风向D的余弦值、气温T、气压P、比容a、比湿H和地表粗糙度R为对风速的影响因素；将采集的数据构建主成分分析的原始数据，如下所示：X1,X2,......,Xm分别表示风向D的正弦值、风向D的余弦值、气温T、气压P、比容a、比湿H和地表粗糙度R的因素集，这里m＝7；b、对7个因素集的n×7个原始数据(X)标准化处理，以消除变量间在数量级上的差异产生的影响，使每个变量的均值为0，方差为1，得到标准化矩阵Y＝{yn...

【专利技术属性】
技术研发人员：张亚刚，王鹏卉，王增平，
申请(专利权)人：华北电力大学保定，
类型：发明
国别省市：河北,13