本发明属于风电场功率预测预报技术领域，尤其是涉及一种基于集合局地多点不同气候特征的风功率预测方法及装置。
背景技术：
风电是一种清洁可再生能源，到2015年年底，全球风电累计装机容量达到432419MW，其中，中国占据总容量的48.4％。由于风力发电受风速和风向变化的影响，这使得风力发电具有明显的波动性、间歇性和随机性，对电力系统的安全及稳定的运行造成了极大的影响，而通过准确的风功率预测，对风电场的出力进行短期预报，可以有力的调控电场运维，减少电量的损失，因此准确的风功率预测会带给风电这一能源形式带来更高的价值。
风功率预测依赖于对各类气象要素的精确预报，其中主要包括风速、风向、湿度等，目前国内外预报的方式主要分为两种，一种是基于历史数据，通过统计学的方法对未来进行预测；另一种是通过不同尺度下的数值天气预报，得到各类气象要素预测值，通过神经网络、卡尔曼滤波等技术最终得到未来的日前风功率预测。目前数值天气预报正值飞速发展，现已可提供未来1-3天高时空分辨率高精度的气象预报。国内短期15分钟分辨率的风功率预报，由于其严格的精度考核要求，因此第二种方法更加适应现有的需求，与此同时，还需要通过不同的方法来提高风功率预测精度。
目前的数值天气预报，主要使用全球高分辨集合结果，但多数将全球尺度集合预报结果，通过中小尺度数值天气模式进行物理降尺度到几公里，人为的选出风电场的地理坐标，通过插值方法将格点的数据计算到电场的单点坐标上，获得风电场的预报风速，最终通过统计模型，将预报的风速转换成风功率。
大气的混沌特性决定了即使采用全球最好的数值模式也很难预报单点的气象结果，模式自身的系统误差与错误预报会间接导致风速的预报误差。
由于现在功率建模技术是基于全场风机平均风速，表现的是整场总体气候状态，根据电场风机排布仅进行粗略的人为选择单点很难描述全场平均的风速状态，亦不能准确的抓住电场的气候状态，因此，预报风速并不适于风电场。
电场空间跨度较大，不同电场所处海拔高度亦不同，边界层内风速随高度变化十分明显，会导致不同电场气候特点有很大差异，仅使用数值模式中单层的或者轮毂高度的结果，将会十分不准确。
综上所述，根据的风机排布，人为选取一个地理坐标，最终将数值天气预报的格点数据插值到该点，如此无法准确的反应整个风电场的整体气候特点，最终导致风功率预测精度的不理想。
技术实现要素：
有鉴于此，本发明旨在提出一种基于集合局地多点不同气候特征的风功率预测方法，通过选择风电场最优的若干点，进行多点合成，可以很好的反映整场的风速状况，来增强气象预报的稳定性，显著提高风功率预测的精度。
为达到上述目的，本发明的技术方案是这样实现的：
一种基于集合局地多点不同气候特征的风功率预测方法，包括如下步骤：
(1)建立电场模型，获得电场周围气象数值预报结果，并通过数据插值方法获得单点预报风速在三维空间的数值；
(2)将单点预报风速合成两点和三点风速，获得合成风速；
(3)利用神经网络方法将合成风速转化成风功率。
进一步的，所述步骤(1)中具体包括如下步骤：
(a1)用Google Earth进行3维电场建模，根据电场风机的排布，选择30公里x 30公里的矩形试验区域机型；
(a2)利用数值天气模式通过并行计算得到试验区域天气预报结果；
(a3)使用NCO(netCDF Operators)提取模式垂直方向1-4层的模式计算结果，并输出模式的地面高度数据、位势高度等，将4层的σ坐标系转换成得到高度坐标系；
(a4)将垂直方向的模式计算结果进行插值计算，获得三维空间1000-3000个单点风速预报结果。
进一步的，所述步骤(a4)中插值计算方法包括线性插值、三次样条插值。
进一步的，所述步骤(2)具体包括如下步骤：
(b1)计算单点预报风速在三维空间的数值的风速预报与电场风机的平均风速之间的相关性系数r，剔除相关性系数后50％的单点结果；
(b2)对保留的50％的单点，分别给出所有两点、三点的组合方案；
(b3)对所有组合方案分别取两点和三点均值，得到合成风速结果；
(b4)计算所有的合成风速结果与电场风机平均风速之间的相关性系数r、平均绝对误差MAE；
(b5)删选出相关性系数最优的前200组组合方案，从200组方案中选出平均绝对误差最小的50组组合方案；
(b6)根据单电场历史实测的风速给出高低风速段的阈值；
(b7)对高低风速段进行权重的订正。
进一步的，所述步骤(3)具体包括如下步骤：
(c1)通过风电场历史风速和功率，经过神经网络模型拟合风速-功率曲线；
(c2)将合成风速放入到神经网络模型中得到预测的风功率；
(c3)通过计算50组预测的风功率的扣电量和精度；
(c4)比较扣电量和精度选出最优的方案。
相对于现有技术，本发明所述的一种基于集合局地多点不同气候特征的风功率预测方法具有以下优势：
(1)本发明在原有单点预报的基础上，整合了单点的结果，形成多点集合预报，减少风速波动，提高风功率预测的精度；
(2)本发明集合的气象数据所选用水平和垂直三维高分辨率格点，综合了不同水平点和不同高度的气候状态，提高风功率预测的精度；
本发明通过整合多点的气候特征，降低风场波动性，不受电场扩容限电等因素的影响，调整预报风速，提高风功率预测的精度。
本发明的另一目的在于提出一种基于集合局地多点不同气候特征的风功率预测装置，以增强气象预报的稳定性，显著提高风功率预测的精度。
为达到上述目的，本发明的技术方案是这样实现的：
一种基于集合局地多点不同气候特征的风功率预测系统，包括
用于建立电场模型，获得电场周围气象数值预报结果，并通过数据插值方法获得单点预报风速在三维空间的数值的电场模型建立装置；
用于将单点预报风速合成两点和三点风速，获得合成风速的风速合成装置；
用于利用神经网络方法将合成风速转化成风功率的风功率转化装置。
进一步的，电场模型建立装置包括
用Google Earth进行三维电场建模，根据电场风机的排布，选择30公里x 30公里的矩形试验区域机型的三维电场建模装置；
用于利用数值天气模式通过并行计算得到试验区域天气预报结果的天气预报结果计算装置；
用于使用NCO(netCDF Operators)提取模式垂直方向1-4层的模式计算结果，并输出模式的地面高度数据、位势高度等，将4层的σ坐标系转换成得到高度坐标系的模式计算结果提取装置；
用于将垂直方向的模式计算结果进行插值计算，获得三维空间1000-3000个单点风速预报结果的插值计算装置。
进一步的，所述插值计算装置包括线性插值计算装置、三次样条插值计算装置。
进一步的，风速合成装置包括
用于计算单点预报风速在三维空间的数值的风速预报与电场风机的平均风速之间的相关性系数r，剔除相关性系数后50％的单点结果的单点结果提取装置；
用于对保留的50％的单点，分别给出所有两点、三点的组合方案的组合方案获得装置；
用于对所有组合方案分别取两点和三点均值，得到合成风速结果的合成风速获得装置；
用于计算所有的合成风速结果与电场风机平均风速之间的相关性系数r、平均绝对误差MAE的相关性系数获得装置以及平均绝对误差获得装置；
用于删选出相关性系数最优的前200组组合方案，从200组方案中选出平均绝对误差最小的50组组合方案最优方案获得装置；
用于根据单电场历史实测的风速给出高低风速段阈值的阈值获取装置；
用于对高低风速段进行权重订正的订正装置。
进一步的，风功率转化装置包括
用于通过风电场历史风速和功率，经过神经网络模型拟合风速-功率曲线的曲线拟合装置；
用于将成风速放入到神经网络模型中得到预测的风功率的风功率预测装置；
用于通过计算50组预测的风功率的扣电量和精度的扣电量和精度获取；
用于比较扣电量和精度选出最优的方案的最优方案获取装置。
所述基于集合局地多点不同气候特征的风功率预测装置与上述基于集合局地多点不同气候特征的风功率预测方法相对于现有技术所具有的优势相同，在此不再赘述。
附图说明
构成本发明的一部分的附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
图1为本发明实施例所述的一种基于集合局地多点不同气候特征的风功率预测方法示意图。
具体实施方式
需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。
下面将参考附图并结合实施例来详细说明本发明。
本发明加入了集合多点优化的理念，重点考虑风电场所处区域大范围的天气状况，水平与垂直三维多点数值预报结果与实测风速、风向的比较，及通过多种方法进行多点合成，最终通过模型转化成风功率结果，再次检验风功率预测的准确性。具体步骤如图1所示，本发明是采用下述的技术方案实现的：
a.电场周围局地多点气象结果的计算；
a.1用Google Earth进行3维电场建模，根据电场风机的排布，选择30公里x 30公里的矩形试验区域，30公里的范围可以涵盖所有风电场的风机，同时可以考虑到电场上游的风场状况，有利于不同季节情况下对电场风速的调优；
a.2利用数值天气模式通过并行计算得到试验区域天气预报结果，数据水平分辨率为3公里、垂直分辨率约50米；
a.3使用NCO(netCDF Operators)提取模式垂直方向1-4层(最高层约200m)的模式计算结果，并输出模式的地面高度数据、位势高度等，将4层的σ坐标系转换成得到高度坐标系；
a.4由于所预报的风速要代表整场的气候特征，同一风电场不同的风机的海拔有很大差距，不同风电场所处的地形和海拔也不尽相同，仅仅将得到的风速插值到轮毂高度是很难满足整场风速的要求。因此需要细化垂直分层，通过边界层风速变化规律见公式1，将垂直插值为20层，最终获得三维空间2000个单点风速预报结果，如此高的分辨率可以精细的考虑风电场所处区域内气候的空间特征，经过试验可得，2000个单点风速预报结果的效果最佳，公式1如下：
其中，Zb，Z:标准参考高度和任一高度
标准参考高度处的平均风速和任一高度处的平均风速
α：地面粗糙指数。
b.合成风速结果，与实测风速对比，筛选合成方案；
b.1为优化计算速度，首先将2000个单点的风速预报分别于电场风机的平均风速进行比较；
b.2考虑到国家电网对风功率预测的考核要求，比较标准为：二者相关性系数r、平均绝对误差MAE、均方根误差RMS等，剔除相关性系数后50％的单点结果；相关性系数计算公式如公式2所示：
其中，r为相关性系数，xi、yi分别表示实发风速和预测风速的第i个值，分别是实发风速和预测风速的算术平均值，n为样本总量。
b.3b.2中保留的1000个单点，分别给出所有两点、三点的组合方案，如此综合两个或三个单点的风速状态，相比单点，多点组合的预报气象结果日内或逐月的波动性均显著减小；
b.4所有组合方案分别取两点和三点均值，得到合成风速结果；
b.5将所有的合成风速结果与电场风机平均风速比较，比较标准为相关性系数r、平均绝对误差MAE；
b.6删选出相关系数最优的前200组组合方案，从200组方案中选出平均绝对误差最小的50组组合方案；
b.7根据单电场历史实测的风速给出高低风速段的阈值；
b.8对高低风速段进行权重的订正。
c.将风速合成方案通过模型处理，对比功率结果；
c.1风速转化成风功率，主要依赖神经网络的方法；
c.2通过风电场历史风速和功率，经过神经网络模型拟合风速-功率曲线；
c.3将50组合成风速放入到神经网络模型中得到50组预测的风功率；
c.4通过计算50组预测的风功率的扣电量和精度；
c.5比较扣电量和精度选出最优的方案。
本发明还提出了一种基于集合局地多点不同气候特征的风功率预测系统，包括
(1)用于建立电场模型，获得电场周围气象结果，并通过数据插值方法获得单点风速在三维空间的数值的电场模型建立装置；
所述电场模型建立装置包括
用Google Earth进行三维电场建模，根据电场风机的排布，选择30公里x 30公里的矩形试验区域机型的三维电场建模装置；
用于利用数值天气模式通过并行计算得到试验区域天气预报结果的天气预报结果计算装置；
用于使用NCO(netCDF Operators)提取模式垂直方向1-4层的模式计算结果，并输出模式的地面高度数据、位势高度等，将4层的σ坐标系转换成得到高度坐标系的模式计算结果提取装置；
用于将垂直方向的模式计算结果进行插值计算，获得三维空间1000-3000个单点风速预报结果的插值计算装置。
所述插值计算装置包括线性插值计算装置、三次样条插值计算装置。
(2)用于将单点风速合成两点和三点风速，获得合成风速的风速合成装置；
所述风速合成装置包括
用于计算单点风速在三维空间的数值的风速预报与电场风机的平均风速之间的相关性系数r，剔除相关性系数后50％的单点结果的单点结果提取装置；
用于对保留的50％的单点，分别给出所有两点、三点的组合方案的组合方案获得装置；
用于对所有组合方案分别取两点和三点均值，得到合成风速结果的合成风速获得装置；
用于计算所有的合成风速结果与电场风机平均风速之间的相关性系数r、平均绝对误差MAE的相关性系数获得装置以及平均绝对误差获得装置；
用于删选出相关性系数最优的前200组组合方案，从200组方案中选出平均绝对误差最小的50组组合方案最优方案获得装置；
用于根据单电场历史实测的风速给出高低风速段阈值的阈值获取装置；
用于对高低风速段进行权重订正的订正装置。
(3)用于利用神经网络方法将合成风速转化成风功率的风功率转化装置；
所述风功率转化装置包括
用于通过风电场历史风速和功率，经过神经网络模型拟合风速-功率曲线的曲线拟合装置；
用于将成风速放入到神经网络模型中得到预测的风功率的风功率预测装置；
用于通过计算50组预测的风功率的扣电量和精度的扣电量和精度获取；
用于比较扣电量和精度选出最优的方案的最优方案获取装置。
本发明可以十分简便的筛选出最优的组合方案，便于之后的实际应用；经实际业务测试后对于不同气候特征的电场风速的预报均有提高，同时对于复杂的山地、丘陵等地形同样均有显著的效果。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。