背景技术
目前,风力发电是国际上具有广阔应用前景的绿色清洁能源之一,风能的准确预报对风能的合理调配利用、电网稳定、商业运营、以及决策服务等方面有着非常重要的作用,它一般采用统计预报和动力预报。通过对风电场风场的预报,结合风电场发电量的情况,可以得到风电场风电量的预报,即风能预报。风能预报的时间尺度,有短期预报(如日预报)和超短期预报(如小时预报)。
20世纪90年代初期,欧洲一些国家就已经开始研制开发风能预报系统并用于预报服务。预报技术多采用中期天气预报模式嵌套高分辨率有限区域模式对风电场发电量进行预报,如丹麦的Prediktor预报系统,目前已用于丹麦、西班牙、爱尔兰和德国的短期风能预报业务。同时,风电功率预测工具(Wind PowerPrediction Tool,简称WPPT),也用于欧洲一些地区的风能预报业务。
90年代中期以后,美国True Wind Solutions 公司也开始商业化的风能预报服务,他们研发的风能预报软件eWind,是由高分辨率的中尺度气象数值模式和统计学模式构成的用于风场和发电量的预报系统。eWind和Prediktor目前在美国加利福尼亚同时用于两个大型风电场的预报服务。
2002年10月,欧盟委员会资助启动了“为陆地和离岸大规模风电场建设发展下一代风资源预报系统”(ANEMOS)计划,目标是发展优于现有方法的、先进的预报模式,重点强调复杂地形和极端气象条件下的预报,同时也发展近海风能预报。加拿大的风能资源数值评估预报软件WEST是将中尺度气象模式MC2与WASP相结合制作分辨率为100~200 m 的风能图谱并进行预报。另外,目前用于风能业务预报的系统,还有德国的 Previento、西班牙的LocalPred 和 RegioPred、以及爱尔兰与丹麦的HIRPOM等。
因此,改进基于天气预报产品为初值的中尺度模式,通过统计降尺度方法,是改善风电场风速预报的主流方法,需要针对某些复杂的地形(如甘肃河西地区)和下垫面,应用数值预报和统计的方法,开发适合本地区的风速预报方法和流程,开展风电场必须的各层风速短期和临近预报。
在现有技术中,一般是改进基于天气预报产品为初值的中尺度模式,通过统计降尺度方法改善风电场风速预报,如何提高复杂地形和极端气象条件下的预报精度目前还是一个技术难题。针对甘肃河西地区复杂的地形和下垫面,没有有效的风速短期和临近预报方法,另外现有模式对突发性天气的预报精度较差,对于风塔10~20min风速的预报存在一定的难度。
在实现本发明的过程中,发明人发现现有技术中至少存在预测精度低、预测准确性差与适用范围小等缺陷。
发明内容
本发明的目的在于,针对上述问题,提出一种包含上下游效应实时监测的超短期预测方法,以实现预测精度高、预测准确性好与适用范围大的优点。
为实现上述目的,本发明采用的技术方案是:一种包含上下游效应实时监测的超短期预测方法,包括:
a、基于T639全球谱模式场库数据源、CALMET风场诊断模型、以及静态数据,通过模式格点,利用WRF-RUC系统,应用WRF3DVAR变分同化技术,获取超短期模式预报结果;
b、基于目标风电基地风塔资料数据库、并结合上下游参考指标站风向风速实时监测数据,进行数值分析与统计,建立相应的参考指标站与各目标风塔之间的上下游效应统计方程,对目标风塔上下游效应进行运算;
c、基于各目标风塔上下游效应的运算结果,对各目标风塔未来超短期的风速变化进行预测预报,并结合所述超短期模式预报结果进行订正,形成目标风电基地风塔的超短期风速变化的预测预报;
d、经多次循环进行上述操作,获取目标风电基地各风电场在目标区域内各高度层未来超短期风速变化的预测预报。
进一步地,步骤a具体包括:
a1、基于CALMET风场诊断模型,对静态数据进行处理,降低WRF中尺度数值预报模式的尺度,生成模式格点;
a2、基于T639全球谱模式场库数据源,读取T639全球谱模式同化资料,对T639全球谱模式同化资料中GRIB格式的气象场资料进行解析,并插值到相应的模式格点;
a3、基于模式格点上的气象场信息,生成初始场和边界条件;利用WRF-RUC系统,应用WRF3DVAR变分同化技术,通过分析,建立用于循环式进行积分预报运送的模式主程序;
a4、启动模式主程序,进行循环运行,实现超短时预报,获取超短期模式预报结果。
进一步地,在步骤a中还包括:
a5、利用绘图设备输出模式产品,输出超短期模式预报结果。
进一步地,在步骤a2中,将T639全球谱模式同化资料中GRIB格式的气象场资料的解析结果插值到相应的模式格点的操作,具体包括:
a21、将T639全球谱模式同化资料中GRIB格式的气象场资料的解析结果,水平插值到相应的模式格点;
a22、将T639全球谱模式同化资料中GRIB格式的气象场资料的解析结果,垂直插值到相应的模式格点。
进一步地,在步骤a21之前,还包括:
获取T639全球谱模式同化资料中GRIB格式的气象场资料的解析结果后,进行格点化准备操作或模式格点化操作。
进一步地,在步骤a4中,所述循环运行具体包括每天8次循环运行;在该8次循环运行中,1200UTC起始的循环为冷启动,其他时次为热启动。
进一步地,步骤b具体包括:
b1、基于目标风电基地风塔资料数据库,获取参考风塔实况监测数据;
b2、针对每个目标风塔,通过最优子集方法,筛选得到不同风向上下游效应相关性最好的参考指标站;
b3、基于参考风塔实况监测数据,通过数值分析与统计,建立参考指标站与各目标风塔之间的上下游效应统计方程;
b4、根据上下游参考指标站风向风速实时监测数据,通过相应的参考指标站与各目标风塔之间的上下游效应统计方程,对目标风塔上下游效应进行运算;
b5、基于各目标风塔上下游效应的运算结果,对各目标风塔未来超短期的风速变化进行预测预报。
进一步地,步骤b4所述上下游参考指标站风向风速实时监测数据,包括风速的上下游效应与高低空效应。
进一步地,在步骤b3中,所述参考指标站与各目标风塔之间的上下游效应统计方程,包括未来0-3小时风电基地风速预报方程;
在步骤b5中,所述未来超短期的设定时间包括5-10分钟。
优选地,在步骤d中,所述目标区域包括10-120m范围内间隔10m的区域,所述各高度层包括10m高度层、70m高度层与100m高度层;所述目标风电基地各风电场在目标区域内各高度层未来超短期风速变化的预测预报中,预报有效时效60h,预报间隔为15min。
本发明各实施例的包含上下游效应实时监测的超短期预测方法,由于包括:基于T639全球谱模式场库数据源、CALMET风场诊断模型、以及静态数据,通过模式格点,利用WRF-RUC系统,应用WRF3DVAR变分同化技术,获取超短期模式预报结果;基于目标风电基地风塔资料数据库、并结合上下游参考指标站风向风速实时监测数据,进行数值分析与统计,建立相应的参考指标站与各目标风塔之间的上下游效应统计方程,对目标风塔上下游效应进行运算;基于各目标风塔上下游效应的运算结果,对各目标风塔未来超短期的风速变化进行预测预报,并结合超短期模式预报结果进行订正,形成目标风电基地风塔的超短期风速变化的预测预报;经多次循环进行上述操作,获取目标风电基地各风电场在目标区域内各高度层未来超短期风速变化的预测预报;可以预报目标风塔10-20min风速,开展风电场必须的各层风速短期和临近预报;从而可以克服现有技术中预测精度低、预测准确性差与适用范围小的缺陷,以实现预测精度高、预测准确性好与适用范围大的优点。
本发明的其它特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
具体实施方式
以下结合附图对本发明的优选实施例进行说明,应当理解,此处所描述的优选实施例仅用于说明和解释本发明,并不用于限定本发明。
根据本发明实施例,如图1和图2所示,提供了一种包含上下游效应实时监测的超短期预测方法。
如图1所示,本实施例的包含上下游效应实时监测的超短期预测方法,包括:
步骤100:基于T639全球谱模式场库数据源,读取T639全球谱模式同化资料;
在步骤100中,T639模式是TL639L60全球谱模式的简称,是国家气象中心在T213L31全球谱模式的基础上升级发展而来。T639模式对一般性降水的雨区范围、位置及移动趋势都做出了较正确的预报,24 h小雨预报TS评分为56;对主要影响系统如高原槽、低空西南急流及西南涡和亚欧中高纬大尺度环流背景均有较准确的刻画;在各种物理量场检验中,反映暴雨动力结构的涡度场、散度场及全风速的模式预报性能稍差,反映水汽条件的比湿及水汽通量散度场预报效果较好;
步骤101:基于CALMET风场诊断模型,对静态数据进行处理,降低WRF中尺度数值预报模式的尺度,生产模式格点;
步骤102:对步骤100所得T639全球谱模式同化资料中GRIB格式的气象场资料,进行解析,获取解析结果,并将解析结果插值到步骤101所得模式格点的相应格点中;
步骤103:基于步骤102中相应模式格点上的气息场信息,生成初始场和边界条件;
步骤104:基于步骤103所得初始场和边界条件,利用WRF-RUC系统,应用WRF3DVAR变同化技术,通过分析,建立用于循环式进行积分预报运算的模式主程序;
步骤105:启动步骤104所得模式主程序,进行每天多次循环运行,实现超短时预报;
在步骤105中,可以使模式主程序进行每天4次的循环运行,在每天的4次循环运行中,12UTC和00UTC起始的循环为冷启动,其他时次为热启动;
步骤106:利用绘图设备输出模式主程序的模式产品,即输出超短期模式预报结果;
步骤107:基于目标风电基地风塔资料数据库,获取参考风塔实况监测数据;
步骤108:基于步骤107所得参考风塔实况监测数据,针对每个目标风塔,通过最优子集方法,筛选得到不同风向上下游效应相关性最好的参考指标站;
步骤109:基于参考风塔实况监测数据,通过数值分析与统计,建立步骤108所得参考指标站与各目标风塔之间的上下游效应统计方程;
在步骤109中,参考指标站与各目标风塔之间的上下游效应统计方程,包括未来0-3小时风电基地风速预报方程;
步骤110:根据上下游参考指标站风向风速实时监测数据,通过相应的参考指标站与各目标风塔之间的上下游效应统计方程,对目标风塔上下游效应进行运算;
在步骤110中,上下游参考指标站风向风速实时监测数据,包括风速的上下游效应与高低空效应;
步骤111:基于步骤110所得目标风塔上下游效应运算结果,对目标风塔未来超短期(如未来5-10分钟)的风速变化进行预测预报;
步骤112:结合步骤106所得超短期模式预报结果,对步骤111所得目标风塔未来5-10分钟超短期的风速变化预测预报结果,进行订正,形成目标风电基地风塔的超短期风速变化的预测预报;
步骤113:经多次循环上述操作,进而完成目标风电基地(即目标区域的风电基地)各风电场10-120m间隔10米各高度层(可以包括10m高度层、70m高度层与100m高度层)未来超短期风速变化的预测预报;
在步骤113中,目标风电基地各风电场在目标区域内各高度层未来超短期风速变化的预测预报中,预报时效为48h,预报间隔为1h。
如图2所示,在上述实施例中,将T639全球谱模式同化资料中GRIB格式的气象场资料的解析结果插值到相应的模式格点的操作,具体包括:
步骤200:读取T639全球谱模式同化资料中GRIB格式的气象场资料,并进行解析,执行步骤201或步骤205或步骤206;
步骤201:获取T639全球谱模式同化资料中GRIB格式的气象场资料的解析结果后,进行格点化准备操作,执行步骤202;
步骤202:经步骤201格点化准备操作后,将T639全球谱模式同化资料中GRIB格式的气象场资料的解析结果,水平插值到相应的模式格点,执行步骤203;
步骤203:经步骤203水平插值完成后,将T639全球谱模式同化资料中GRIB格式的气象场资料的解析结果,垂直插值到相应的模式格点,执行步骤204;
步骤204:获取模式格点上的气象场信息;
步骤205:获取T639全球谱模式同化资料中GRIB格式的气象场资料的解析结果后,进行模式格点化操作,执行步骤202或步骤206;
步骤206:将T639全球谱模式同化资料中GRIB格式的气象场资料的解析结果,水平插值到相应的模式格点,执行步骤207;
步骤207:经步骤206水平插值完成后,将T639全球谱模式同化资料中GRIB格式的气象场资料的解析结果,垂直插值到相应的模式格点,执行步骤204。
上述实施例的包含上下游效应实时监测的超短期预测方法,基于T639全球谱模式同化资料,应用改进和优化的WRF中尺度数值预报模式,开展酒泉风电基地风电场10-120m间隔10米各高度层风速预报;预报有效时效60h,预报间隔为15min,水平分辨率为3km。利用CALMET风场诊断模型,对WRF中尺度模式降尺度,提高预报能力和精度;利用WRF-RUC系统,应用WRF3DVAR变分同化技术,通过分析和循环,每天共计4次循环运行(具体启动方式为:12UTC和00UTC起始的循环为冷启动,其他时次为热启动,实现超短时预报。该预测方法属于气象风能预测技术领域,可以应用于电力行业风电调度。
例如,可以针对甘肃省河西地区复杂的地形和下垫面,应用数值预报和统计的方法,开发适合本地区的风速预报方法和流程,通过分析风塔风速变化特征,研究上游风塔风速的变化、传播对下游风塔的影响,最终建立上下游预报方程,预报目标风塔10~20min风速,开发短时预报订正技术,开展风电场必须的各层风速短期和临近预报。
通过统计酒泉地区风塔变化特征,利用酒泉风电基地区域内风向盛行偏东风和偏西风,风向有较好的一致性,风速的传播存在一定规律的特性,通过分析20个风塔风速变化特征,研究上游风塔风速的变化、传播对下游风塔的影响,最终建立上下游预报方程,预报目标风塔10~20min风速,开发短时预报订正技术;研究风速的上下游效应和高低空效应,建立未来0-3小时风电基地风速预报方程;利用统计学方法,对酒泉风电基地风塔资料及气象测站、自动站资料进行分析,研究风塔风速时空变化规律及风速随高度变化特征,分析风塔风速与周围测站风速的关系,找出上游指标站。
又如,可以结合甘肃省酒泉地区天气气候特点,基于国内T639同化资料,应用优化的WRF中尺度数值预报模式,开展酒泉风电基地风电场风速预报,预报水平精度9km,预报层次10m、70m、100m,有效预报时效60小时,预报间隔15min。尝试用CALMET风场诊断模型,提高风要素预报的水平精度;利用风塔上下游效应,开发短时预报订正预报技术;建立酒泉风电基地风功率预报系统,开发预报准确率较高业务化风电预报产品,为风电调度提供技术支撑。
具体地,可以统计分析酒泉风电基地20个风塔观测资料,针对每个目标风塔挑选出不同风向上下游效应参考指标站,通过最优子集方法筛选相关性最好的参考站,通过资料分析、统计建立参考站与目标风塔之间的上下游效应统计方程。
根据上游参考指标站风向风速实时监测,通过相对应的各风塔上下游效应统计方程运算,预报目标风塔未来5-10分钟的风速变化。结合模式输出预报结果进行订正,形成风塔的超短期风速预测。
上述实施例的包含上下游效应实时监测的超短期预测方法,至少具有以下特点:
⑴改进优化中尺度数值预报模式、降尺度方法和风塔等气象资料同化技术,提高预报准确率;
⑵继续改进短临预报技术,完成风塔布设格距调整,进行风塔观测资料质量控制;
⑶研究上游风塔风速的变化、传播对下游风塔的影响,开发了上下游预报方法,应用于短时临近预报技术,为超短时风力预测提供一种新的思路;
⑷研究风塔风速时空变化规律及风速随高度变化特征,分析风塔风速与周围测站风速的关系,确立风速预报上游参考指标站;
⑸在短时风力预测方面,经过对15分钟70米高度风速预报和风塔资料进行对比效果检验,结果表明,24h内预报相对误差为22.9-30%,平均相对误差为26.97%;绝对误差为1.6-2.3m/s,平均绝对误差为1.8m/s,可以达到电力调度的要求。
最后应说明的是:以上所述仅为本发明的优选实施例而已,并不用于限制本发明,尽管参照前述实施例对本发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。