【技术实现步骤摘要】
本专利技术涉及沙尘预报，具体地讲，是涉及一种基于深度学习的沙尘短临预报方法。
技术介绍
1、沙尘是指特定地理环境和下垫面条件下，由强风将地面大量尘土、沙粒卷入空中，使空气混浊、水平能见度下降的天气现象，具有突发性强、持续时间短、影响范围广的特点，会急剧降低能见度，并严重威胁生态环境安全、城市运行和人体健康。提高沙尘预报的精度是当前气象工作的重点和难点。利用对面观测站网，根据天气现象、能见度等观测要素，可以长时间、全天候的获得沙尘发生的时间、地点、强度等特征，同时基于卫星反演的沙尘实况产品能够时空连续地监测沙尘发展过程，这些沙尘监测数据是开展沙尘研究和预报的基础。
2、现有沙尘预报根据预报方法实现原理的不同可分为数值预报方法和机器学习预报方法。数值预报方法通过动力学方程组模拟沙尘排放、输送到沙尘的干或湿沉降来预报未来的沙尘状态，是目前沙尘预报业务上主要采用的方法。成熟的沙尘数值模式系统包括gocart、dream、itr等全球沙尘数值模式，以及cemsys5、cfors、narcm等区域沙尘数值模式，以及中国的cuace-dust，grapes-cuace/dust等。由于数值模式难以准确表示小规模的初始条件及在垂直和水平方向上都缺乏足够的分辨率，使其更适合长期预报。此外数值模式在生成沙尘形成和演变中涉及的关键参数仍然很困难，特别是在地形复杂或观测有限的地区，这使得预报结果具有一定的不确定性。
3、利用机器学习相关方法进行沙尘预报的主要思路是通过建模气象观测数据与沙尘天气相关观测数据之间的关系，实现对沙尘天
4、沙尘短临预报是沙尘灾害预警预报的重要基础，人们更关注短时临近的沙尘预报效果，预报效率较高的潜势预报方法能帮助业务人员提前确定沙尘发生区域和范围，从而将主要精力放在重点区域内的预警。沙尘带来的低能见度可能导致交通事故、航班取消和道路封闭，从而造成人员伤亡和严重的经济损失。因而迫切需要一种效率更快、效果更好的沙尘短临预报方法，为民航中心、高速管理部门提供更好地气象服务支撑。
技术实现思路
1、针对上述现有技术中的上述问题，本专利技术提供一种基于深度学习的沙尘短临预报方法，将沙尘短临预报问题转换为沙尘实况序列预测问题，通过将多个连续的历史帧沙尘实况产品序列预测图输入网络中来得到几个小时后的短临预测图，通过判断预测图的区域和沙尘浓度来查看临近预报的沙尘发展情况，实现准确快速的沙尘短临预报。
2、为了实现上述目的，本专利技术采用的技术方案如下：
3、一种基于深度学习的沙尘短临预报方法，包括以下步骤：
4、s10、根据地面站点的观测数据筛选出沙尘过程数据，并根据沙尘过程数据的时间和空间匹配已有的基于卫星反演的沙尘实况数据；
5、s20、对步骤s10匹配后的数据进行预处理，形成沙尘时空序列样本集；
6、s30、基于predrnn网络构建时空序列的深度学习模型，并使用步骤s20的沙尘时空序列样本集训练获得沙尘浓度预测模型，表示为：
7、pt+1～pt+n＝argmaxp(xt+1～xt+n|xt-j+1～xt)
8、xt+1～xt+n
9、上式中，xt-j+1～xt表示从时刻t-j+1到时刻t的历史沙尘浓度实况数据序列，xt+1～xt+n表示从时刻t+1到时刻t+n的沙尘浓度真值数据序列，pt+1～pt+n表示从时刻t+1到时刻t+n的沙尘浓度预测结果，即概率p最大的预测沙尘浓度时空数据序列，n<j；
10、s40、在发生沙尘天气时，将已有的基于卫星反演的沙尘实况数据输入步骤s30获得的沙尘浓度预测模型中，输出指定区域未来设定时间段内的沙尘分布结果，进行沙尘短临预报。
11、具体地，所述步骤s10中地面站点的观测数据包括地面国家气象站逐小时全要素数据中的能见度和天气现象数据，以及地面环保站的pm10数据；所述已有的基于卫星反演的沙尘实况数据为按时间序列的包含沙尘区和非沙尘区标记的沙尘格点图像。
12、具体地，所述步骤s10中根据沙尘过程数据和已有的基于卫星反演的沙尘实况数据进行时间和空间匹配的过程包括：
13、s11、收集历史沙尘过程中地面站点的观测数据，筛选出包含沙尘发生的开始时间和结束时间在内的时间信息以及包含沙尘发生站点的经纬度信息在内的空间信息；
14、s12、根据步骤s11获得的时间信息，筛选基于卫星反演的沙尘实况数据中含有相同时间信息的数据，完成数据的时间匹配；
15、s13、在时间匹配的条件下，根据步骤s11获得的空间信息从基于卫星反演的沙尘实况数据中提取反演结果数据，判断在指定区域内反演结果数据与观测数据是否一致，若是，则将该时间和空间下的基于卫星反演的沙尘实况数据纳入初始样本中。
16、具体地，所述步骤s20中对数据进行预处理，形成沙尘时空序列样本集的过程为：
17、s21、对步骤s10匹配后的数据进行区域裁剪，裁剪后构成256×512的网格；
18、s22、采用滑动窗口模式从裁剪后的数据中以每j+n个连续时刻的数据获取一个样本，并循环获取样本，形成样本集；
19、s23、采用最大最小值归一化方法，对每个样本进行标准化处理，并将空值赋值为0，获得沙尘时空序列样本集。
20、具体地，所述步骤s10中已有的基于卫星反演的沙尘实况数据为geotiff格式，所述步骤s20中沙尘时空序列样本集中样本格式存储为npy格式。
21、具体地，所述步骤s30中基于predrnn网络构建的时空序列的深度学习模型为堆叠式网络结构，输入的沙尘图像经过特征嵌入层从单通道图像转换为隐空间通道维度d为w的沙尘特征图，沙尘特征图经过堆叠式网络结构学习和特征降维层的处理，转换为单通道沙尘预测图像输出。
22、具体地，所述步骤s30中训练获得沙尘浓度预测模型的过程包括：
23、s31、模型训练初始参数设置：模型循环单元层数设置为4，基本单元之间的隐特征维度设置为64；损失函数设置为mse损失函数，初始学习率设定为0.0001；计划采样停止轮数设置为5000；滤波器大小设置为2；采用adamw优化器来最小化mse损失函数，通过其自适应梯度机制和动量梯度机制来调整学习率，加快模型收敛速度；最大迭代轮次设置为200；输入的时空序列长度设置为j，预测时空序列长度设置为n本文档来自技高网...
【技术保护点】
1.一种基于深度学习的沙尘短临预报方法，其特征在于，包括以下步骤：
2.根据权利要求1所述的基于深度学习的沙尘短临预报方法，其特征在于，所述步骤S10中地面站点的观测数据包括地面国家气象站逐小时全要素数据中的能见度和天气现象数据，以及地面环保站的PM10数据；所述已有的基于卫星反演的沙尘实况数据为按时间序列的包含沙尘区和非沙尘区标记的沙尘格点图像。
3.根据权利要求2所述的基于深度学习的沙尘短临预报方法，其特征在于，所述步骤S10中根据沙尘过程数据和已有的基于卫星反演的沙尘实况数据进行时间和空间匹配的过程包括：
4.根据权利要求3所述的基于深度学习的沙尘短临预报方法，其特征在于，所述步骤S20中对数据进行预处理，形成沙尘时空序列样本集的过程为：
5.根据权利要求4所述的基于深度学习的沙尘短临预报方法，其特征在于，所述步骤S10中已有的基于卫星反演的沙尘实况数据为geotiff格式，所述步骤S20中沙尘时空序列样本集中样本格式存储为npy格式。
6.根据权利要求5所述的基于深度学习的沙尘短临预报方法，其特征在于，所述步骤S
7.根据权利要求6所述的基于深度学习的沙尘短临预报方法，其特征在于，所述步骤S30中训练获得沙尘浓度预测模型的过程包括：
8.根据权利要求7所述的基于深度学习的沙尘短临预报方法，其特征在于，所述基于PredRNN网络构建的时空序列的深度学习模型以ST-LSTM单元为网络核心，包括多层堆叠的神经单元，通过时空记忆单元在层之间上下流动，通过长期记忆单元和隐藏层在层与层之间传递，通过输入时刻t-1的输入沙尘浓度图像Xt-1，获得输出时刻t的预测沙尘浓度图像表示为：
9.根据权利要求8所述的基于深度学习的沙尘短临预报方法，其特征在于，所述输入的时空序列长度设置为6，预测时空序列长度设置为3。
...【技术特征摘要】
1.一种基于深度学习的沙尘短临预报方法，其特征在于，包括以下步骤：
2.根据权利要求1所述的基于深度学习的沙尘短临预报方法，其特征在于，所述步骤s10中地面站点的观测数据包括地面国家气象站逐小时全要素数据中的能见度和天气现象数据，以及地面环保站的pm10数据；所述已有的基于卫星反演的沙尘实况数据为按时间序列的包含沙尘区和非沙尘区标记的沙尘格点图像。
3.根据权利要求2所述的基于深度学习的沙尘短临预报方法，其特征在于，所述步骤s10中根据沙尘过程数据和已有的基于卫星反演的沙尘实况数据进行时间和空间匹配的过程包括：
4.根据权利要求3所述的基于深度学习的沙尘短临预报方法，其特征在于，所述步骤s20中对数据进行预处理，形成沙尘时空序列样本集的过程为：
5.根据权利要求4所述的基于深度学习的沙尘短临预报方法，其特征在于，所述步骤s10中已有的基于卫星反演的沙尘实况数据为geotiff格式，所述步骤s20中沙尘时空序列样本集中样本格式存储为npy格式。
6.根据权利要...
【专利技术属性】
技术研发人员：徐艳琴，刘辉，张翔，陶睿，徐亮亮，段晓梅，乔淼，温建伟，张立，
申请(专利权)人：内蒙古自治区气象数据中心内蒙古自治区气象探测中心，内蒙古自治区气象档案馆，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。