一种预驱动中长期空气质量预报系统和方法与流程

本发明涉及空气质量预报领域，特别涉及一种预驱动中长期空气质量预报系统和方法。

背景技术：

目前，随着人民群众对空气质量的关注程度日益提高，国家对空气质量预报提出了新的要求，需要空气质量预报的时长能够覆盖到10天乃至更长，而长时间的空气质量预报难以采用传统的统计预报方法开展，基于gfs的传统数值预报方法体系则难以预报月尺度甚至季节尺度的空气质量。

目前cfs(ncepcoupledforecastsystemmodel)产品广泛用于国内外1-3月中长期气象条件预报，并基于气象条件预报开展了水文、生态等多方面的应用，也包括污染形势分析等大气环境相关的工作，但少见其在空气质量数值预报领域的应用。

数值模拟(含预报)技术通常需要进行spin-up，即初始化模拟，由于模式驱动数据的空间分辨率通常与模拟网格不一致，通过初始模拟可以让模式输入数据在空间尺度上更为接近模拟网格分辨率，达到稳定状态。通常情况气象模式需要通过6-24小时的初始模拟；对于空气质量模式而言，由于空气质量模式的初试场通常为均一分布的统计浓度，实际排放则具备更为明确的空间分布属性，污染物排入大气后达到稳定的混合状态需要的时间更长，初始化模拟需要的时间更多，因此通常需要7-15天的初始化模拟时间以期获得更为合理的空气质量模拟结果。

在基于gfs产品的短期空气质量预报中，由于ncepfnl资料是gfs系统整合全球观测数据后的再分析产品，可以使用fnl资料对历史气象条件进行模拟，并通过restartrun方式提供初始场以改善预报效果，但在基于cfs产品的中长期空气质量预报中，由于cfs系统与gfs系统由不同模式产出，垂直层次等属性存在差异，无法使用restartrun的方法为模式提供已完成spin-uip的高时效性初始场数据，从而无法使用fnl预驱动cfs中长期数值预报。

在进行数值预报的时候，由于预报场数据的质量通常不如结合了实测数据订正后的再分析资料，且预报场数据的准确性随着预报时间的增加逐渐降低，故对于数值预报而言初始化模拟更为重要，可一定程度上提高早期空气质量数值预报的准确性。对于基于cfs产品的中长期数值预报而言，因其预报时长通常在一个月到三个月之间，早期预报效果相较后期效果更佳，若使用cfs的早期产品进行初始化模拟，由于空气质量模拟所需的初始化模拟时间较长，污染物浓度并未混合到稳定状态，从而导致该阶段的空气质量预报效果不佳，浪费了cfs产品气象条件预报效果相对较好的阶段。

技术实现要素：

本发明所要解决的技术问题是如何基于再分析数据进行预驱动从而提高cfs中长期预报精度。

本发明解决上述技术问题的技术方案如下：一种预驱动中长期空气质量预报系统，包括：气象数据下载模块、数值模型驱动模块和后处理模块；

所述气象数据下载模块用于下载起报时间前10天的再分析资料和起报时间的cfs预报场资料；

所述数值模型驱动模块用于读取用户提供的排放清单，并利用所述再分析资料和所述排放清单得到起报时间前10天的气象场回顾模拟结果和空气质量回顾模拟结果；

所述数值模型驱动模块还包括：气象初始数据提取子模块和污染物浓度初始数据提取子模块；所述气象初始数据提取子模块根据所述气象场回顾模拟结果得到预报开始第一时间的气象初始数据；所述污染物浓度数据提取子模块根据所述空气质量回顾模拟结果得到预报开始第一时间的污染物浓度初始数据，所述数值模型驱动模块用于使用所述气象初始数据和所述污染物浓度初始数据；替换原始wrfinput数据中基于插值的粗分辨率气象场初值和icon数据中基于统计的平均污染物浓度初值，从而利用回顾模拟结果改善cfs预报前期预报效果，实现预驱动；

所述cfs预报场资料包括cfs地面预报场数据和cfs高空预报场数据；

所述数值模型驱动模块还用于根据所述cfs预报场数据及所述排放清单得到中长期气象条件预报结果及中长期空气质量预报结果；

所述后处理模块用于根据所述中长期气象条件预报结果及所述空气质量预报结果提取点位预报产品，并生成预报报告。

本发明的有益效果是：本方案从回顾性模拟结果中直接获取模式初始化需要的变量，经处理后写入预报数据的初始场中，利用回顾模拟结果对预报进行初始化，实现预驱动。相比现有技术中采用的restartrun热启动方式，本方案所需要的数据不需要具有相同的数据结构，而是采用数值模型驱动模块这个中间程序处理模拟结果，获取初始变量用于预报的初始化工作，故该方法对前后两次模拟的驱动数据无要求，具备较高的灵活性，也解决了cfs系统无法使用fnl等成熟再分析资料初始化cfs中长期数值预报的技术问题，从而达到了提高中长期空气质量预报精度的技术效果；结合再分析资料得到的回顾模拟结果对未来空气质量变化趋势相比使用默认初始值对未来空气质量变化趋势的预报更有参考意义。

在上述技术方案的基础上，本发明还可以做如下改进。

进一步，所述再分析资料包括fnl再分析资料、era-interim或era5资料。

进一步，所述后处理模块包括产品绘图子模块和点位预报数据提取子模块；

所述产品绘图子模块用于根据所述中长期气象条件预报结果及所述中长期空气质量预报结果绘制风速、风向、温度、相对湿度、降水、太阳辐射、no2、so2、pm10、pm2.5、o3、co、ox的空间分布图件；

所述点位预报数据提取子模块用于根据空气质量监测站点经纬度信息从所述中长期气象条件预报结果中提取点位对应的气象要素时间序列预报数据，还用于根据空气质量监测站点经纬度和所述中长期空气质量预报结果提取点位对应的空气质量时间序列预报数据，并根据所述气象要素时间序列预报数据、所述空气质量时间序列预报数据和我国aqi计算规范计算点位aqi数据；

所述后处理模块还用于结合所述空间分布图件、气象要素时间序列预报数据、空气质量时间序列预报数据及aqi数据生成预报报表。

本发明解决上述技术问题的另一技术方案如下：

一种预驱动中长期空气质量预报方法，包括以下步骤：

s1：下载起报时间前10天的再分析资料和起报时间的cfs预报场资料；

s2：读取用户排放清单，利用再分析资料和所述排放清单得到起报时间前10天的气象场回顾模拟结果和空气质量回顾模拟结果；

s3：根据所述气象场回顾模拟结果得到预报开始第一时间的气象初始数据，根据所述空气质量回顾模拟结果得到预报开始第一时间的污染物浓度初始数据；

s4：利用所述预报开始第一时间的气象初始数据替换cfs数据生成的wrfinput中的初始气象数据，并利用所述预报开始第一时间的污染物浓度初始数据替换cmaq模型icon中的初始污染物浓度数据，并分别完成中长期气象数值预报和空气质量数值预报；

s5：根据所述中长期气象条件预报结果及中长期空气质量预报结果得到气象要素及污染物的空间分布图件，及点位气象要素时间序列预报数据和空气质量时间序列预报数据；

s6：根据所述空间分布图件、点位气象要素时间序列预报数据和空气质量时间序列预报数据生成预报报告。

进一步，步骤s2具体为：步骤s2中，所述再分析资料包括fnl再分析资料、era-interim或era5资料。

进一步，步骤s5具体为，根据所述中长期气象条件预报结果及所述中长期空气质量预报结果绘制风速、风向、温度、相对湿度、降水、太阳辐射、no2、so2、pm10、pm2.5、o3、co、ox的空间分布图件；根据空气质量监测站点经纬度信息从所述中长期气象条件预报结果中提取点位对应的气象要素时间序列预报数据，还根据空气质量监测站点经纬度和所述中长期空气质量预报结果提取点位对应的空气质量时间序列预报数据，并根据所述气象要素时间序列预报数据、所述空气质量时间序列预报数据和我国aqi计算规范计算点位aqi数据。

本发明附加的方面的优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明实践了解到。

附图说明

图1为本发明预驱动中长期空气质量预报系统的实施例的系统结构示意图；

图2为本发明预驱动中长期空气质量预报方法的实施例的方法流程图。

具体实施方式

以下结合附图对本发明的原理和特征进行描述，所举实例只用于解释本发明，并非用于限定本发明的范围。

实施例基本如附图1所示：

本实施例中预驱动中长期空气质量预报方法，包括：气象数据下载模块1、数值模型驱动模块2和后处理模块3；

气象数据下载模块1用于下载起报时间前10天的再分析资料4和起报时间的cfs预报场资料5，所述数值模型驱动模块2用于读取用户提供的排放清单13，并利用所述再分析资料和所述排放清单13得到起报时间前10天的气象场回顾模拟结果和空气质量回顾模拟结果；

本实施例中的排放清单13可以包括预报城市及周边地区vocs、nox、co、so2、pm10、pm2.5、nh3等大气污染物排放数据，具备时间、空间和化学物种属性，用于描述不同时间不同位置的大气污染排放情况，可为城市相关科研机构编制，或按区域截取自全国排放清单，本实施例中的气象数据下载子模块用于下载起报时间前10天的再分析资料4和起报时间的cfs预报场资料5；

本实施例中的数值模型驱动模块2包括气象模拟子模块和cmaq空气质量模拟子模块，其中气象模拟子模块用于得到气象场回顾模拟结果6，cmaq空气质量模拟子模块用于得到空气质量回顾模拟结果；数值模型驱动模块2还包括：气象初始数据提取子模块和污染物浓度初始数据提取子模块；气象初始数据提取子模块根据气象场回顾模拟结果6得到预报开始第一时间的气象初始数据7；污染物浓度数据提取子模块根据空气质量回顾模拟结果得到预报开始第一时间的污染物浓度初始数据8，数值模型驱动模块用于使用气象初始数据7和污染物浓度初始数据8替换原始wrfinput数据中基于插值的粗分辨率气象场初值和icon数据中基于统计的平均污染物浓度初值，从而利用回顾模拟结果6改善cfs预报前期预报效果，实现预驱动；

cfs预报场资料5包括cfs地面预报场数据和cfs高空预报场数据；

所述数值模型驱动模块2还用于根据所述cfs预报场资料5及所述排放清单13得到中长期气象条件预报结果12及中长期空气质量预报结果11；

所述后处理模块3用于根据所述中长期气象条件预报结果12及所述空气质量预报结果提取点位预报产品，并生成预报报告10，本实施例中的点位预报产品为空气质量监测站点的中长期气象预报及空气质量预报，包括数据和图标等内容；中长期气象条件预报结果12是指未来1-3个月的气象条件数值预报，中长期空气质量预报结果11是指未来1-3个月的空气质量数值预报，本实施例中的点位预报产品包括气象预报子模块，气象预报子模块用于获取中长期气象预报。

本发明的有益效果是：本方案从回顾性模拟结果中直接获取模拟变量，并提取模拟结果中模式初始化需要的变量，经处理后写入cfs预报的初始场中，利用回顾模拟结果对预报进行初始化，实现预驱动。相比现有技术中采用的restartrun热启动方式，本方案所需要的数据不需要具有相同的数据结构，而是采用数值模型驱动模块2这个中间程序处理模拟结果，获取初始变量用于预报的初始化工作，故该方法对前后两次模拟的数据结构无要求，具备较高的灵活性，也解决了cfs系统无法使用fnl等成熟再分析资料初始化cfs中长期数值预报的技术问题，从而达到了提高中长期空气质量预报精度的技术效果；结合再分析资料4得到的回顾模拟结果6对未来空气质量变化趋势的预报相比使用默认初始值对未来空气质量变化趋势更有参考意义。

在上述技术方案的基础上，本发明还可以做如下改进。

可选的，在一些其它实施例中，再分析资料4包括fnl再分析资料、era-interim或era5资料。

可选的，在一些其它实施例中，后处理模块3包括产品绘图子模块和点位预报数据9提取子模块；

产品绘图子模块用于根据中长期气象条件预报结果12及中长期空气质量预报结果11绘制风速、风向、温度、相对湿度、降水、太阳辐射、no2、so2、pm10、pm2.5、o3、co、ox的空间分布图件；

点位预报数据9提取子模块用于根据空气质量监测站点经纬度信息从中长期气象条件预报结果12中提取点位对应的气象要素时间序列预报数据，还用于根据空气质量监测站点经纬度和中长期空气质量预报结果11提取点位对应的空气质量时间序列预报数据，并根据气象要素时间序列预报数据、空气质量时间序列预报数据和我国aqi计算规范计算点位aqi数据；

后处理模块3还用于结合空间分布图件、气象要素时间序列预报数据、空气质量时间序列预报数据及aqi数据生成预报报表。

本发明解决上述技术问题的另一技术方案如下：

如附图2所示，一种预驱动中长期空气质量预报方法，包括以下步骤：

s1：下载起报时间前10天的再分析资料4和起报时间的cfs预报场资料5；

s2：读取用户排放清单13，利用再分析资料4和排放清单13得到起报时间前10天的气象场回顾模拟结果6和空气质量回顾模拟结果；

s3：根据气象场回顾模拟结果6得到预报开始第一时间的气象初始数据7，根据空气质量回顾模拟结果得到预报开始第一时间的污染物浓度初始数据8；

s4：利用预报开始第一时间的气象初始数据7替换cfs数据生成的wrfinput中的初始气象数据，并利用预报开始第一时间的污染物浓度初始数据8替换cmaq模型icon中的初始污染物浓度数据，并分别完成中长期气象数值预报和空气质量数值预报；

s5：根据中长期气象条件预报结果12及中长期空气质量预报结果11得到气象要素及污染物的空间分布图件，及点位气象要素时间序列预报数据和空气质量时间序列预报数据；

s6：根据空间分布图件、点位气象要素时间序列预报数据和空气质量时间序列预报数据生成预报报告10。

可选的，在一些其它实施例中，步骤s2具体为：步骤s2中，再分析资料4包括fnl再分析资料、era-interim或era5资料。

可选的，在一些其它实施例中，步骤s5具体为，根据中长期气象条件预报结果12及中长期空气质量预报结果11绘制风速、风向、温度、相对湿度、降水、太阳辐射、no2、so2、pm10、pm2.5、o3、co、ox的空间分布图件；根据空气质量监测站点经纬度信息从中长期气象条件预报结果12中提取点位对应的气象要素时间序列预报数据，还根据空气质量监测站点经纬度和中长期空气质量预报结果11提取点位对应的空气质量时间序列预报数据，并根据气象要素时间序列预报数据、空气质量时间序列预报数据和我国aqi计算规范计算点位aqi数据。

需要说明的是，上述各实施例是与上述各方法实施例对应的产品实施例，对于本实施例中各结构装置及可选实施方式的说明可以参考上述各方法实施例中的对应说明，在此不再赘述。

读者应理解，在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。

所属领域的技术人员可以清楚地了解到，为了描述的方便和简洁，上述描述的装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。

以上，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到各种等效的修改或替换，这些修改或替换都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以权利要求的保护范围为准。