CN103927596A - 基于复合数据源自回归模型的风电功率超短期预测方法 - Google Patents

Abstract

本发明公开了一种基于复合数据源自回归模型的风电功率超短期预测方法，包括输入数据得到自回归模型参数；以及输入风电功率预测所需输入数据到根据自回归模型的参数确定的自回归模型中得到预测结果；输入数据得到自回归模型参数具体包括，输入模型训练基础数据，采用残差方差图法对自回归模型AR(p)定阶，采用矩估计方法对定阶的AR(p)模型参数进行估计。对风力发电过程中的风电功率进行预测，为新能源发电实时调度、新能源发电日前计划、新能源发电月度计划、新能源发电能力评估和弃风电量估计提供关键信息。引入复合数据源有效提高风电功率超短期预测精度，从而实现在保障电网安全稳定经济运行的前提下有效提高新能源上网电量目的。

Description

基于复合数据源自回归模型的风电功率超短期预测方法

技术领域

本发明涉及新能源发电过程中风电功率预测技术领域，具体地，涉及一种基于复合数据源自回归模型的风电功率超短期预测方法。

背景技术

我国风电进入规模化发展阶段以后所产生的大型新能源基地多数位于“三北地区”(西北、东北、华北)，大型新能源基地一般远离负荷中心，其电力需要经过长距离、高电压输送到负荷中心进行消纳。由于风、光资源的间歇性、随机性和波动性，导致大规模新能源基地的风电、光伏发电出力会随之发生较大范围的波动，进一步导致输电网络充电功率的波动，给电网运行安全带来一系列问题。

截至2014年4月，甘肃电网并网风电装机容量已达707万千瓦，约占甘肃电网总装机容量的22％，成为仅次于火电的第二大主力电源。目前，甘肃电网风电、光伏发电装机超过甘肃电网总装机容量的1/3。随着新能源并网规模的不断提高，风电、光伏发电不确定性和不可控性给电网的安全稳定经济运行带来诸多问题。准确预估可利用的发电风资源是对大规模风电优化调度的基础。对风力发电过程中的风电功率进行预测，可为新能源发电实时调度、新能源发电日前计划、新能源发电月度计划、新能源发电能力评估和弃风电量估计提供关键信息。

发明内容

本发明的目的在于，针对上述问题，提出一种基于复合数据源自回归模型的风电功率超短期预测方法，以实现在保障电网安全稳定经济运行的前提下有效提高新能源上网电量的优点。

为实现上述目的，本发明采用的技术方案是：

一种基于复合数据源自回归模型的风电功率超短期预测方法，包括输入数据得到自回归模型参数，以及输入风电功率预测所需输入数据到根据上述自回归模型的参数确定的自回归模型中得到预测结果；

所述输入数据得到自回归模型参数具体包括步骤101、输入模型训练基础数据，

步骤102、采用残差方差图法对自回归模型AR(p)定阶，

步骤103、采用矩估计方法对定阶的AR(p)模型参数进行估计。

根据本发明的优选实施例，所述步骤101输入模型训练基础数据，输入数据包括，风电场基础信息、历史风速数据、历史功率数据和地理信息系统数据。

根据本发明的优选实施例，所述步骤102采用残差方差图法对自回归模型AR(p)定阶：

具体为设x_t为需要估计的项，x_t-1,x_t-2,...,x_t-n为已知历史功率序列，自回归模型AR(p)，模型定阶就是确定模型中参数p的值；

用系列阶数逐渐递增的模型拟合原始序列，每次都计算残差平方和然后画出阶数和的图形，当阶数由小增大时，会显著下降，达到真实阶数后的值会逐渐趋于平缓，甚至反而增大，

拟合误差的平方和/(实际观测值个数-模型参数个数)，

实际观测值个数指拟合模型时实际使用的观察值项数，对于具有N个观察值的序列，拟合AR(p)模型，则实际使用的观察值最多为N-p，模型参数个数指所建立的模型中实际包含的参数个数，对于含有均值的模型，模型参数个数为模型阶数加1，对于N个观测值的序列，AR模型的残差估计式为：

根据本发明的优选实施例，所述步骤103采用矩估计方法对定阶的AR(p)模型参数进行估计具体步骤为：

将风电场历史功率数据利用数据序列x₁,x₂,...,x_t表示，其样本自协方差定义为

${\widehat{\mathrm{\γ}}}_k=\frac{1}{n}\underset{t=k+1}{\overset{n}{\mathrm{\Σ}}}{x}_t{x}_{t-k},$

其中，k＝0,1,2,...,n-1，x_t和x_t-k均为数据序列x₁,x₂,...,x_t中的数值；

则 ${\widehat{\mathrm{\γ}}}_0=\frac{1}{n}\underset{t=1}{\overset{n}{\mathrm{\Σ}}}{x}_t^2$

则历史功率数据样本自相关函数为：

${\widehat{\mathrm{\ρ}}}_k=\frac{{\widehat{\mathrm{\γ}}}_k}{{\widehat{\mathrm{\γ}}}_0}=\frac{\frac{1}{n}\underset{t=k+1}{\overset{n}{\mathrm{\Σ}}}{x}_t{x}_{t-k}}{\frac{1}{n}\underset{t=1}{\overset{n}{\mathrm{\Σ}}}{x}_t^2}=\frac{\underset{t=k+1}{\overset{n}{\mathrm{\Σ}}}{x}_t{x}_{t-k}}{\underset{t=1}{\overset{n}{\mathrm{\Σ}}}{x}_t^2},$

其中，k＝0,1,2,...,n-1。

AR部分的矩估计为，

令

则协方差函数为

用的估计代替γ_k，

可得参数

根据本发明的优选实施例，所述输入风电功率预测所需输入数据到根据上述自回归模型的参数确定的自回归模型中得到预测结果的步骤包括，

步骤201、输入功率预测基础数据；

步骤202、对输入的基础数据进行噪声滤波及数据预处理；

步骤203、根据确定的参数建立自回归模型，并将处理后的数据输入从而得到预测结果；

步骤204、将预测结果输出，并通过图表及曲线展示预测结果。

根据本发明的优选实施例，所述输入功率预测基础数据包括资源监测系统数据和运行监测系统数据，所述资源监测系统数据包含风资源监测数据；所述运行监测系统数据包括风机监测数据、升压站监测数据和数据采集与监视控制系统数据。

根据本发明的优选实施例，所述噪声滤波及数据预处理具体为：噪声滤波模块对监测系统实时采集得到的带有噪声的数据进行滤波处理，去除坏数据和奇异值；数据预处理模块对数据进行对齐、归一化处理和分类筛选处理。

根据本发明的优选实施例，所述自回归模型为：

其中，是系数，α_t是白噪声序列。

本发明的技术方案具有以下有益效果：

本发明的技术方案通过对风力发电过程中的风电功率进行预测，为新能源发电实时调度、新能源发电日前计划、新能源发电月度计划、新能源发电能力评估和弃风电量估计提供关键信息。通过引入复合数据源有效提高风电功率超短期预测精度，从而实现在保障电网安全稳定经济运行的前提下有效提高新能源上网电量目的。

下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。

附图说明

图1为本发明实施例所述的基于复合数据源自回归模型的风电功率超短期预测方法的原理框图。

具体实施方式

以下结合附图对本发明的优选实施例进行说明，应当理解，此处所描述的优选实施例仅用于说明和解释本发明，并不用于限定本发明。

一种基于复合数据源自回归模型的风电功率超短期预测方法，包括输入数据得到自回归模型参数，

以及输入风电功率预测所需输入数据到根据上述自回归模型的参数确定的自回归模型中得到预测结果；其中输入数据得到自回归模型参数具体包括步骤101、输入模型训练基础数据，

步骤102、采用残差方差图法对自回归模型AR(p)定阶，

步骤103、采用矩估计方法对定阶的AR(p)模型参数进行估计。

含大规模风电的电力系统运行依赖庞大的、准确的数据集，而风电功率预测若能将这些数据有效融合利用则可有效提高预测精度。与常规电力系统SCADA监测不同，在各类电气、机械和热力等数据之外，风电监测数据还包含大量的资源监测、运行监测及地理信息等。

如图1所示，本发明技术方案提出的风电功率超短期预测可分为两个阶段：模型训练阶段和功率预测阶段。

阶段1：模型训练

步骤1.1：模型训练基础数据输入

风功率预报系统模型训练所需输入数据包括，风电场基础信息、历史风速数据、历史功率数据，地理信息系统(GIS)数据(风电场/风机坐标、测风塔坐标、升压站坐标等)。将基础数据输入到预测模型中进行模型训练。

步骤1.2：模型定阶

由于事先无法确定需要使用多少已知时间序列的项来建立估计函数，所以需要对模型进行定阶判断。

设x_t为需要估计的项，x_t-1,x_t-2,...,x_t-n为已知历史功率序列，对于自回归模型AR(p)，模型定阶就是确定模型中参数p的值。

采用残差方差图法进行模型定阶。假定模型是有限阶自回归模型，如果设置的阶数小于真实阶数，则是一种不足拟合，因而拟合残差平方和必定偏大，此时通过提高阶数可以显著降低残差平方和。反之，如果阶数已经达到真实值，那么再增加阶数，就是过度拟合，此时增加阶数不会令残差平方和显著减小，甚至会略有增加。

这样用系列阶数逐渐递增的模型来拟合原始序列，每次都计算残差平方和然后画出阶数和的图形。当阶数由小增大时，会显著下降，达到真实阶数后的值会逐渐趋于平缓，有时甚至反而增大。残差方差的估计式为：

拟合误差的平方和/(实际观测值个数-模型参数个数)

“实际观测值个数”是指拟合模型时实际使用的观察值项数，对于具有N个观察值的序列，拟合AR(p)模型，则实际使用的观察值最多为N-p。“模型参数个数”是指所建立的模型中实际包含的参数个数，对于含有均值的模型，模型参数个数为模型阶数加1。对于N个观测值的序列，相应AR模型的残差估计式为：

      (式1)

其中，公式中，Q为拟合误差的平方和函数，是模型系数，N是观测序列长度，是模型参数中的常数项，的常识值是根据不同的进行变化的常数项，不同的对照不同的值。

步骤1.3：模型参数估计

采用矩估计方法对ARMA(p)的模型参数进行估计。首先，将风电场历史功率数据利用数据序列x₁,x₂,...,x_t表示，其样本自协方差定义为

${\widehat{\mathrm{\γ}}}_k=\frac{1}{n}\underset{t=k+1}{\overset{n}{\mathrm{\Σ}}}{x}_t{x}_{t-k}$         (式2)

其中，k＝0,1,2,...,n-1，x_t和x_t-k均为数据序列x₁,x₂,...,x_t中的数值。

特别的，

${\widehat{\mathrm{\γ}}}_0=\frac{1}{n}\underset{t=1}{\overset{n}{\mathrm{\Σ}}}{x}_t^2$           (式3)

则历史功率数据样本自相关函数为：

${\widehat{\mathrm{\ρ}}}_k=\frac{{\widehat{\mathrm{\γ}}}_k}{{\widehat{\mathrm{\γ}}}_0}=\frac{\frac{1}{n}\underset{t=k+1}{\overset{n}{\mathrm{\Σ}}}{x}_t{x}_{t-k}}{\frac{1}{n}\underset{t=1}{\overset{n}{\mathrm{\Σ}}}{x}_t^2}=\frac{\underset{t=k+1}{\overset{n}{\mathrm{\Σ}}}{x}_t{x}_{t-k}}{\underset{t=1}{\overset{n}{\mathrm{\Σ}}}{x}_t^2}$        (式4)其中，k＝0,1,2,...,n-1。

AR部分的矩估计为

        (式5)

令

           (式6)

则协方差函数为

      (式7)

用的估计代替γk，有

            (式8)

可得参数

通过上述求解过程发现，要求解时间序列模型的阶数，就要得到时间序列的预测值；要得到时间序列的预测值，必须先建立具体的预测函数；要建立具体的预测函数，必须知道模型的阶数。

根据实践验证，时间序列模型阶数一般不超过5阶。所以在该算法具体实现时，可以首先假设模型为1阶，利用步骤1.3中的参数估计方法得到一阶模型的参数，进而建立估计函数便可以求得一阶模型时间序列模型估计得到各个项的预测值，从而求得一阶模型的残差方差；之后，假设模型为二阶，用上述方法求得二阶模型的残差；以此类推，可以得到1到5阶模型的残差，选残差最小的模型的阶数作为最终模型的阶数。确定模型阶数后，便可计算得到参数的值。

阶段2：功率预测

步骤2.1：功率预测基础数据输入

风电功率预测所需输入数据包括资源监测系统数据和运行监测系统数据两部分，其中，资源监测系统数据包含风资源监测数据；运行监测系统数据包括风机监测数据、升压站监测数据和数据采集与监视控制系统(SCADA)数据等。

步骤2.2：噪声滤波及数据预处理

噪声滤波模块对监测系统实时采集得到的带有噪声的数据进行滤波处理，去除坏数据和奇异值；数据预处理模块对数据进行对齐、归一化处理和分类筛选等操作，以便使得输入的数据可以为模型所用。

步骤2.3：超短期功率预测

将模型参数估计出来之后，结合已估计的模型阶数，便可得到用于风电功率超短期预测的时间序列方程。根据上述步骤2和步骤3得出的p值，以及的值建立自回归模型；

自回归模型如下：

          (式9)

其中，是系数，α_t是白噪声序列。

步骤2.4：预测结果输出及展示

该步骤首先对预测结果进行输出，并通过图形和表格等形式对预测结果进行展示。

最后应说明的是：以上所述仅为本发明的优选实施例而已，并不用于限制本发明，尽管参照前述实施例对本发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (9)

1.一种基于复合数据源自回归模型的风电功率超短期预测方法，其特征在于，包括输入数据得到自回归模型参数；

以及输入风电功率预测所需输入数据到根据上述自回归模型的参数确定的自回归模型中得到预测结果；

所述输入数据得到自回归模型参数具体包括步骤101、输入模型训练基础数据，

步骤102、采用残差方差图法对自回归模型AR(p)定阶，

步骤103、采用矩估计方法对定阶的AR(p)模型参数进行估计。

2.根据权利要求1所述的基于复合数据源自回归模型的风电功率超短期预测方法，其特征在于，所述步骤101输入模型训练基础数据，输入数据包括，风电场基础信息、历史风速数据、历史功率数据和地理信息系统数据。

3.根据权利要求2所述的基于复合数据源自回归模型的风电功率超短期预测方法，其特征在于，所述步骤102采用残差方差图法对自回归模型AR(p)定阶：

具体为设x_t为需要估计的项，x_t-1,x_t-2,...,x_t-n为已知历史功率序列，自回归模型AR(p)，模型定阶就是确定模型中参数p的值；

用系列阶数逐渐递增的模型拟合原始序列，每次都计算残差平方和然后画出阶数和的图形，当阶数由小增大时，会显著下降，达到真实阶数后的值会逐渐趋于平缓，甚至反而增大，

拟合误差的平方和/(实际观测值个数-模型参数个数)，

实际观测值个数指拟合模型时实际使用的观察值项数，对于具有N个观察值的序列，拟合AR(p)模型，则实际使用的观察值最多为N-p，模型参数个数指所建立的模型中实际包含的参数个数，对于含有均值的模型，模型参数个数为模型阶数加1，对于N个观测值的序列，AR模型的残差估计式为：

其中，Q为拟合误差的平方和函数，是模型系数，N是观测序列长度，是模型参数中的常数项。

4.根据权利要求3所述的基于复合数据源自回归模型的风电功率超短期预测方法，其特征在于，所述步骤103采用矩估计方法对定阶的AR(p)模型参数进行估计具体步骤为：

将风电场历史功率数据利用数据序列x₁,x₂,...,x_t表示，其样本自协方差定义为

${\widehat{\mathrm{\γ}}}_k=\frac{1}{n}\underset{t=k+1}{\overset{n}{\mathrm{\Σ}}}{x}_t{x}_{t-k},$

其中，k＝0,1,2,...,n-1，x_t和x_t-k均为数据序列x₁,x₂,...,x_t中的数值；

则 ${\widehat{\mathrm{\γ}}}_0=\frac{1}{n}\underset{t=1}{\overset{n}{\mathrm{\Σ}}}{x}_t^2$

则历史功率数据样本自相关函数为：

${\widehat{\mathrm{\ρ}}}_k=\frac{{\widehat{\mathrm{\γ}}}_k}{{\widehat{\mathrm{\γ}}}_0}=\frac{\frac{1}{n}\underset{t=k+1}{\overset{n}{\mathrm{\Σ}}}{x}_t{x}_{t-k}}{\frac{1}{n}\underset{t=1}{\overset{n}{\mathrm{\Σ}}}{x}_t^2}=\frac{\underset{t=k+1}{\overset{n}{\mathrm{\Σ}}}{x}_t{x}_{t-k}}{\underset{t=1}{\overset{n}{\mathrm{\Σ}}}{x}_t^2},$

其中，k＝0,1,2,...,n-1；

AR部分的矩估计为，

令

则协方差函数为

用的估计代替γ_k，

可得参数

5.根据权利要求4所述的基于复合数据源自回归模型的风电功率超短期预测方法，其特征在于，所述输入风电功率预测所需输入数据到根据上述自回归模型的参数确定的自回归模型中得到预测结果的步骤包括，

步骤201、输入功率预测基础数据；

步骤202、对输入的基础数据进行噪声滤波及数据预处理；

步骤203、根据确定的参数建立自回归模型，并将处理后的数据输入从而得到预测结果。

6.根据权利要求5所述的基于复合数据源自回归模型的风电功率超短期预测方法，其特征在于，还包括，

步骤204、将预测结果输出，并通过图表及曲线展示预测结果。

7.根据权利要求6所述的基于复合数据源自回归模型的风电功率超短期预测方法，其特征在于，所述输入功率预测基础数据包括资源监测系统数据和运行监测系统数据，所述资源监测系统数据包含风资源监测数据；所述运行监测系统数据包括风机监测数据、升压站监测数据和数据采集与监视控制系统数据。

8.根据权利要求6所述的基于复合数据源自回归模型的风电功率超短期预测方法，其特征在于，所述噪声滤波及数据预处理具体为：噪声滤波模块对监测系统实时采集得到的带有噪声的数据进行滤波处理，去除坏数据和奇异值；数据预处理模块对数据进行对齐、归一化处理和分类筛选处理。

9.根据权利要求6所述的基于复合数据源自回归模型的风电功率超短期预测方法，其特征在于，所述自回归模型为：

其中，是系数，α_t是白噪声序列。