CN114279915A - 一种大气颗粒物浓度反演方法及相关组件 - Google Patents
一种大气颗粒物浓度反演方法及相关组件 Download PDFInfo
- Publication number
- CN114279915A CN114279915A CN202111604123.2A CN202111604123A CN114279915A CN 114279915 A CN114279915 A CN 114279915A CN 202111604123 A CN202111604123 A CN 202111604123A CN 114279915 A CN114279915 A CN 114279915A
- Authority
- CN
- China
- Prior art keywords
- data set
- signal
- wind
- noise ratio
- atmospheric particulate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 239000008277 atmospheric particulate matter Substances 0.000 title claims abstract description 102
- 238000000034 method Methods 0.000 title claims abstract description 58
- 238000012544 monitoring process Methods 0.000 claims abstract description 20
- 238000005259 measurement Methods 0.000 claims description 90
- 238000003062 neural network model Methods 0.000 claims description 27
- 238000013528 artificial neural network Methods 0.000 claims description 19
- 230000008569 process Effects 0.000 claims description 18
- 238000012549 training Methods 0.000 claims description 18
- 238000007781 pre-processing Methods 0.000 claims description 17
- 238000012935 Averaging Methods 0.000 claims description 16
- 238000012937 correction Methods 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 6
- 238000013527 convolutional neural network Methods 0.000 claims description 6
- 125000004122 cyclic group Chemical group 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 4
- 238000003915 air pollution Methods 0.000 abstract description 7
- 230000008033 biological extinction Effects 0.000 abstract description 7
- 230000008878 coupling Effects 0.000 abstract description 5
- 238000010168 coupling process Methods 0.000 abstract description 5
- 238000005859 coupling reaction Methods 0.000 abstract description 5
- 230000004044 response Effects 0.000 abstract description 4
- 238000004088 simulation Methods 0.000 abstract description 2
- 239000003344 environmental pollutant Substances 0.000 description 10
- 231100000719 pollutant Toxicity 0.000 description 10
- 239000000356 contaminant Substances 0.000 description 8
- 230000009471 action Effects 0.000 description 5
- 230000008859 change Effects 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000009792 diffusion process Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 239000012080 ambient air Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 230000010349 pulsation Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004062 sedimentation Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000013618 particulate matter Substances 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000005180 public health Effects 0.000 description 1
- 230000000541 pulsatile effect Effects 0.000 description 1
- 230000000306 recurrent effect Effects 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Optical Radar Systems And Details Thereof (AREA)
Abstract
本发明公开了一种大气颗粒物浓度反演方法及相关组件,依据空气污染动力学中的标量守恒方程来确定需要获取的数据,将第一风数据集、第一湿度、第二信噪比数据集、第一能见度及第一湍流强度数据集输入至预先训练好的大气颗粒物浓度反演模型以得到第一预设时长内整层大气的第一大气颗粒物浓度。该方案建立了多普勒测风激光雷达、地面气象站监测数据与整层大气的大气颗粒物浓度之间的强耦合关系,且由于多普勒测风激光雷达可以实时观测测量高度阈范围内的整层大气,因此与现有技术相比,该方案能够获得整层大气的大气颗粒物浓度,响应稳定,准确度高,减少了利用基于激光雷达消光系数的经验模型反演大气颗粒物浓度容易造成的模拟误差,反演精度更高。
Description
技术领域
本发明涉及大气污染监测领域,特别是涉及一种大气颗粒物浓度反演方法及相关组件。
背景技术
随着城市工业化发展,大气污染日益严重,PM10和PM2.5作为影响环境空气质量的主要污染物,越来越受到人们的关注。当大气中颗粒物浓度升高时,环境空气质量会随之恶化,这严重影响了公众健康和城市形象,限制了城市发展。为此,环境监测部门常用地面监测仪器对近地面的大气颗粒物浓度及分布情况进行实时监测,并在城市内部或城市之间组成监测网络。但是地面监测仪器只能获取预设监测区域内的近地面的大气颗粒物浓度情况,却无法反映整层大气中的大气颗粒物浓度,因此,寻找一个能够有效反映整层大气中的大气颗粒物浓度方法是至关重要的。
近年来,随着激光雷达技术的逐渐发展,为大气颗粒物大范围、高频率的探测研究提供了有利工具,该技术可实时观测整层大气,为了解大气颗粒物光学特性在整层大气中的空间分布和时间变化及其环境效应的监测研究提供了可靠的技术手段。
在现有技术中,利用激光雷达无法直接测量大气颗粒物的浓度,通常需要首先利用激光雷达的信噪比计算出消光系数,再利用基于激光雷达消光系数的经验模型反演大气颗粒物浓度。但在实际生活中,大气颗粒物的浓度不仅与激光雷达的消光系数相关,而且受大气中物理和化学因素的影响也很大,但该方法中的基于激光雷达消光系数的经验模型并未考虑到这一点,因此缺少空气污染动力学方面的理论支撑;且目前该经验模型在形式上还停留在线性模型和指数模型上,其中包括较多需要人为设定参数,针对这些参数的确定过程往往存在较大误差,使得利用该经验模型得到的反演结果与地面监测仪器得到的监测结果相差较大,普适性较差,往往无法多次使用。
发明内容
本发明的目的是提供一种大气颗粒物浓度反演方法及相关组件,建立了多普勒测风激光雷达、地面气象站测得的数据与整层大气的大气颗粒物浓度之间的强耦合关系,能够获得整层大气的大气颗粒物浓度,响应稳定,准确度高。
为解决上述技术问题,本发明提供了一种大气颗粒物浓度反演方法,应用于多普勒测风激光雷达中的处理器,包括:
获取所述多普勒测风激光雷达在第一预设时长且在测量高度范围内测得的第一风数据集、第一信噪比数据集及第一湿度;
获取地面气象站在第一预设时长内测得的第一能见度;
基于所述第一风数据集及所述多普勒测风激光雷达的预存的风速标准差数据集确定第一湍流强度数据集;
对所述第一信噪比数据集进行预处理后得到第二信噪比数据集;
将所述第一风数据集、所述第二信噪比数据集、所述第一湿度、所述第一能见度及所述第一湍流强度数据集作为输入项输入预先训练好的大气颗粒物浓度反演模型;
确定所述大气颗粒物浓度反演模型的输出项为所述第一预设时长内整层大气的第一大气颗粒物浓度。
优选的,所述第一风数据集包括所述多普勒测风激光雷达的所述测量高度范围内的与各测量返回高度一一对应的第一风向和第一风速;
所述第一信噪比数据集包括与各所述测量返回高度一一对应的第一信噪比A;
所述第一湍流强度数据集包括与各所述测量返回高度一一对应的第一湍流强度B。
优选的,对所述第一信噪比数据集进行预处理后得到第二信噪比数据集,包括:
基于预设信噪比校正关系式对所述第一信噪比数据集进行校正后得到第二信噪比数据集,所述第二信噪比数据集包括与所述第一信噪比数据集中的各个第一信噪比A一一对应的校正后的各个第二信噪比C;
所述预设信噪比校正关系式为C=A*D*E;
其中,D为校正因子;E为与所述A对应的所述测量返回高度的平方。
优选的,基于所述第一风数据集及所述多普勒测风激光雷达的预存的风速标准差数据集确定第一湍流强度数据集,包括:
根据所述多普勒测风激光雷达的预存的风速标准差数据集及预设湍流强度关系式对所述第一风数据集进行处理后确定第一湍流强度数据集;
所述预设湍流强度关系式为B=F/G,其中F为与各个所述B一一对应的各个所述第一风速,G为所述风速标准差数据集中的与各所述测量返回高度一一对应的风速标准差。
优选的,将所述第一风数据集、所述第二信噪比数据集、所述第一湿度、所述第一能见度及所述第一湍流强度数据集作为输入项输入预先训练好的大气颗粒物浓度反演模型,包括:
对于第i个所述测量返回高度,其中1≤i≤所述测量返回高度的总个数且i为整数,执行如下步骤:
将所述第一能见度、所述第一湿度、与第i个所述测量返回高度对应的所述第一风向、所述第一风速、所述第二信噪比C、所述第一湍流强度B作为输入项输入预先训练好的大气颗粒物浓度反演模型以得到与第i个所述测量返回高度对应的输出项;
确定所述大气颗粒物浓度反演模型的输出项为所述第一预设时长内整层大气的第一大气颗粒物浓度,包括:
确定与第i个所述测量返回高度对应的输出项为所述第一预设时长内与第i个所述测量返回高度对应的第一大气颗粒物浓度。
优选的,所述大气颗粒物浓度反演模型的预先训练过程,包括:
S11:获取所述地面环境监测站在各个第二预设时长内测得的整层大气的第二大气颗粒物浓度,其中,所述第二预设时长的个数为预设训练基准阈值;
S12:获取所述多普勒测风激光雷达测得的与第j个第二预设时长内的各个第三预设时长对应的第二风数据集、第三信噪比数据集及第二湿度,其中,第j个所述第二预设时长包括H个所述第三预设时长,H≥1且H为整数,1≤j≤所述预设训练基准阈值且j为整数,其中,所述第二风数据集包括所述多普勒测风激光雷达的所述测量高度范围内的与各所述测量返回高度一一对应的第二风向和第二风速;所述第三信噪比数据集包括与各所述测量返回高度一一对应的第三信噪比;
S13:获取所述地面气象站测得的与各个所述第二预设时长对应的第二能见度;
S14:判断在所有的所述第二湿度中是否存在不小于预设湿度阈值的第二湿度,若是,进入S15;若否,进入S18;
S15:判断在与各个不小于所述预设湿度阈值的第二湿度分别一一对应的各个所述第三信噪比数据集中,与第一预设高度对应的各个第三信噪比中是否包括所述多普勒测风激光雷达返回的无效值,若是,进入S16;若否,进入S17;
S16:将与各个包括所述无效值分别对应的各个所述第三预设时长内的第二风数据集、第三信噪比数据集、第二湿度及第二能见度删除;
S17:对与各个不包括所述无效值分别对应的各个所述第三信噪比数据集进行预处理后得到与各个所述第三信噪比数据集分别一一对应的各个第四信噪比数据集;对与各个不包括所述无效值分别对应的各个所述第二风数据集,分别基于各个所述第二风数据集及所述风速标准差数据集确定与各个所述第二风数据集一一对应的第二湍流强度数据集,其中,所述第四信噪比数据集包括与各所述测量返回高度一一对应的第四信噪比,所述第二湍流强度数据集包括与各所述测量返回高度一一对应的第二湍流强度;
S18:针对与各个小于所述预设湿度阈值的第二湿度分别对应的各个所述第二风数据集,分别基于各个所述第二风数据集及所述风速标准差数据集确定与各个所述第二风数据集一一对应的第二湍流强度数据集;对与各个小于所述预设湿度阈值的第二湿度分别对应的各个所述第三信噪比数据集,进行预处理后得到与各个所述第三信噪比数据集一一对应的各个第四信噪比数据集,其中,所述第四信噪比数据集包括与各所述测量返回高度一一对应的第四信噪比,所述第二湍流强度数据集包括与各所述测量返回高度一一对应的第二湍流强度;
S19:将各个所述第二预设时长及与各个所述第二预设时长一一对应的各个所述第二大气颗粒物浓度作为输出项输入至神经网络模型;
S20:针对与第j个所述第二预设时长对应的H个第三预设时长内的各所述第二风数据集、各所述第四信噪比数据集、各所述第二湿度、各所述第二湍流强度数据集,将各所述第二风数据集中的对应行的对应项进行加和并求平均值后得到第三风数据集,将各所述第四信噪比数据集的对应行的对应项进行加和并求平均值后得到第五信噪比数据集,将各所述第二湿度的对应行的对应项进行加和并求平均值后得到与第j个所述第二预设时长对应的第三湿度,将各所述第二湍流强度数据集的对应行的对应项进行加和并求平均值后得到第三湍流强度数据集;
S21:取所述第三风数据集中与第二预设高度对应的数据作为与第j个所述第二预设时长对应的第四风数据集;取所述第五信噪比数据集中与所述第二预设高度对应的数据作为与第j个所述第二预设时长对应的第六信噪比;取所述第三湍流强度数据集中与所述第二预设高度对应的数据作为与第j个所述第二预设时长对应的第四湍流强度;其中,所述第二预设高度为与各个所述测量返回高度中的最小值对应的高度;
S22:将各个所述第二预设时长、与各个所述第二预设时长分别对应的所述第四风数据集、所述第六信噪比、所述第三湿度、所述第二能见度及所述第四湍流强度作为输入项输入至所述神经网络模型以建立在各个所述第二预设时长下的输入-输出对应关系;
S23:利用所述神经网络模型进行第z次学习,z为大于1的整数;
S24:判断第z次学习的所述神经网络模型的学习精度是否满足预设精度阈值,若是,进入S26;若否,进入S25;
S25:令z=z+1并返回S22;
S26:停止学习并确定此时的所述神经网络模型为所述大气颗粒物浓度反演模型。
优选的,所述神经网络模型包括BP神经网络、卷积神经网络、循环神经网络及长短期记忆神经网络中的一种或多种的组合。
为解决上述技术问题,本发明还提供了一种大气颗粒物浓度反演***,包括:
第一获取单元,用于获取所述多普勒测风激光雷达在第一预设时长且在测量高度范围内测得的第一风数据集、第一信噪比数据集及第一湿度;
第二获取单元,用于获取地面气象站在第一预设时长内测得的第一能见度;
第一确定单元,用于基于所述第一风数据集及所述多普勒测风激光雷达的预存的风速标准差数据集确定第一湍流强度数据集;
第二确定单元,用于对所述第一信噪比数据集进行预处理后得到第二信噪比数据集;
输入单元,用于将所述第一风数据集、所述第二信噪比数据集、所述第一湿度、所述第一能见度及所述第一湍流强度数据集作为输入项输入预先训练好的大气颗粒物浓度反演模型;
输出确定单元,用于确定所述大气颗粒物浓度反演模型的输出项为所述第一预设时长内整层大气的第一大气颗粒物浓度。
为解决上述技术问题,本发明还提供了一种大气颗粒物浓度反演装置,包括:
存储器,用于存储计算机程序;
处理器,用于执行所述计算机程序时实现如上述所述的大气颗粒物浓度反演方法的步骤。
为解决上述技术问题,本发明还提供了一种多普勒测风激光雷达,包括如上述所述的大气颗粒物浓度反演装置。
本发明提供了一种大气颗粒物浓度反演方法及相关组件,依据空气污染动力学中的标量守恒方程来确定需要获取的数据,即获取多普勒测风激光雷达在第一预设时长且在测量高度范围内测得的第一风数据集及第一湿度及第一信噪比数据集,对第一信噪比数据集进行预处理后得到的第二信噪比数据集,并基于第一风数据集及多普勒测风激光雷达的预存的风速标准差数据集确定第一湍流强度数据集,获取地面气象站在第一预设时长内测得的第一能见度,通过将上述数据作为输入项输入至预先训练好的大气颗粒物浓度反演模型,可以得到第一预设时长内整层大气的第一大气颗粒物浓度。该方案建立了多普勒测风激光雷达、地面气象站监测数据与整层大气的大气颗粒物浓度之间的强耦合关系,且由于多普勒测风激光雷达可以实时观测测量高度范围内的整层大气,因此与现有技术相比,该方案能够获得整层大气的大气颗粒物浓度,响应稳定,准确度高,且减少了利用基于激光雷达消光系数的经验模型反演大气颗粒物浓度容易造成的模拟误差,反演精度更高。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对现有技术和实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明提供的一种大气颗粒物浓度反演方法的流程图;
图2为本发明提供的一种大气颗粒物浓度反演***的结构示意图;
图3为本发明提供的一种大气颗粒物浓度反演装置的结构示意图。
具体实施方式
本发明的核心是提供一种大气颗粒物浓度反演方法及相关组件,建立了多普勒测风激光雷达、地面气象站测得的数据与整层大气的大气颗粒物浓度之间的强耦合关系,能够获得整层大气的大气颗粒物浓度,响应稳定,准确度高。
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
请参照图1,图1为本发明提供的一种大气颗粒物浓度反演方法的流程图。
该大气颗粒物浓度反演方法,应用于多普勒测风激光雷达中的处理器,包括:
S31:获取多普勒测风激光雷达在第一预设时长且在测量高度范围内测得的第一风数据集、第一信噪比数据集及第一湿度;
S32:获取地面气象站在第一预设时长内测得的第一能见度;
S33:基于第一风数据集及多普勒测风激光雷达的预存的风速标准差数据集确定第一湍流强度数据集;
S34:对第一信噪比数据集进行预处理后得到第二信噪比数据集;
S35:将第一风数据集、第二信噪比数据集、第一湿度、第一能见度及第一湍流强度数据集作为输入项输入预先训练好的大气颗粒物浓度反演模型;
S36:确定大气颗粒物浓度反演模型的输出项为第一预设时长内整层大气的第一大气颗粒物浓度。
本实施中,考虑到现有技术中为了监测整层大气中的大气颗粒物浓度,一种方式是利用地面监测仪器进行实时监测,但该种方式只能获取预设监测区域内的近地面的大气颗粒物浓度情况;另一种方式是利用激光雷达检测的数据配合基于激光雷达消光系数的经验模型得到大气颗粒物的浓度,但该种方式缺乏理论基础且普适性较差。为解决上述技术问题,本申请依据空气污染动力学中的标量守恒方程来确定需要获取的数据,并预先训练了大气颗粒物浓度反演模型,将获取的数据输入至该大气颗粒物浓度反演模型可以确定整层大气的大气颗粒物浓度。
首先需要说明的是,空气污染动力学中的标量守恒方程如下所示:
其中,表示某区域内单位体积污染物浓度的局地变化,且其中q表示该区域内的污染物浓度,t表示时间;于是可以看出某区域内单位体积污染物浓度的局地变化由污染物的平流输送项湍流扩散项体源项Sc、沉降项R以及化学反应项∑chemis五大部分共同决定。具体的,考虑到大气湍流运动的特性,在边界层气象中可以将湍流运动表示为:
湍流运动=平均运动+脉动运动
同时在大气中风是三维矢量,将其分别在空间坐标系的x、y、z轴上进行投影,分别对应得到U、V、W,则是风在x方向上的平均速度,u′是风在x方向上的脉动速度;是风在y方向上的平均速度,v′是风在y方向上的脉动速度;是风在z方向的平均速度,w′是在z方向上的脉动速度;为污染物的平均部分,q′是污染物的脉动部分。
于是针对平流输送项其表征平均流场下污染物浓度的变化;针对湍流扩散项其表征脉动流场下污染物浓度的变化,其中,为x方向上风速的脉动量和污染物浓度的脉动量的协方差,为y方向上风速的脉动量和污染物浓度的脉动量的协方差、为z方向上风速的脉动量和污染物浓度的脉动量的协方差;针对体源项SC,其表征大气污染源在该区域的排放水平,由于大气污染源的排放在较短时间内的变化不大,可用大气污染物浓度进行表征;针对沉降项R,表示大气污染物在下沉气流、降水等作用下被地面、建筑、植被等捕捉的部分;针对化学反应项∑chemis,表述大气中各污染物相互发生化学反应的部分。因此可以看出上述五项共同作用决定了该区域内单位体积污染物浓度的局地变化。
大气颗粒物可以看作一种污染物,在上述理论部分的基础上可知,想要对大气颗粒物的浓度进行反演,需要确定能够表征上述标量守恒方程中平流输送项、湍流扩散项、体源项、沉降项及化学反应项这五部分的数据。因此,本申请中,考虑到对于大气边界层,可以将其理解为无数个单位体积堆叠起来的空间,其垂直高度即为边界层高度,地面污染物浓度受边界层各高度气象因子的共同作用影响。于是,根据该标量守恒方程,获取多普勒测风激光雷达在第一预设时长且在测量高度范围内测得的第一风数据集、第一信噪比数据集及第一湿度,获取地面气象站在第一预设时长内测得的第一能见度,其中,第一风数据集用以表征该标量守恒方程中的平流输送项;由于下垫面在单个季节内的变化不大,因此将第一湿度和第一能见度一起用以表征该标量守恒方程中的沉降项。需要说明的是,这里的第一预设时长可以为多普勒测风激光雷达的当前返回数据的时间与上一次返回数据的时间之间的时间间隔,如10分钟,本申请在此不作特别的限定,根据具体的多普勒测风激光雷达的工作性能决定;这里的测量高度范围可以为多普勒测风激光雷达能够测量到数据的最大距离,如3公里,本申请在此对该测量高度范围的具体数值不作特别的限定,根据具体的多普勒测风激光雷达的工作性能决定。
随后,基于该第一风数据集及多普勒测风激光雷达的预存的风速标准差数据集确定第一湍流强度数据集,该第一湍流强度数据集用以表征该标量守恒方程中的湍流扩散项,其中该预存的风速标准差数据集为多普勒测风激光雷达根据其工作性能直接设定并预存,本申请在此不作特别的限定;对第一信噪比数据集进行预处理后得到第二信噪比数据集,该第二信噪比数据集用以表征该标量守恒方程中的体源项和化学反应项;将第一风数据集、第二信噪比数据集、第一湿度、第一能见度及第一湍流强度数据集作为输入项输入预先训练好的大气颗粒物浓度反演模型,则确定大气颗粒物浓度反演模型的输出项为第一预设时长内整层大气的第一大气颗粒物浓度。
需要说明的是,为了更直观地展现通过反演得到的第一大气颗粒物浓度,本申请中可以在确定大气颗粒物浓度反演模型的输出项为第一预设时长内整层大气的第一大气颗粒物浓度之后,由处理器控制显示模块显示该第一大气颗粒物浓度,以更直观地方式将第一大气颗粒物浓度显示给开发人员,便于开发人员后续的记录及处理。
综上,本申请提供了一种大气颗粒物浓度反演方法,依据空气污染动力学中的标量守恒方程来确定需要获取的数据,通过预先训练好的大气颗粒物浓度反演模型建立了多普勒测风激光雷达、地面气象站的监测数据与整层大气的大气颗粒物浓度之间的强耦合关系,且由于多普勒测风激光雷达可以实时观测测量高度范围内的整层大气,因此与现有技术相比,该方案能够获得整层大气的大气颗粒物浓度,响应稳定,准确度高,且减少了利用基于激光雷达消光系数的经验模型反演大气颗粒物浓度容易造成的模拟误差,反演精度更高。
在上述实施例的基础上:
作为一种优选的实施例,第一风数据集包括多普勒测风激光雷达的测量高度范围内的与各测量返回高度一一对应的第一风向和第一风速;
第一信噪比数据集包括与各测量返回高度一一对应的第一信噪比A;
第一湍流强度数据集包括与各测量返回高度一一对应的第一湍流强度B。
本申请中,考虑到多普勒测风激光雷达会在测量高度范围内的各测量返回高度内返回与各测量高度对应的数据,于是第一风数据集可以包括多普勒测风激光雷达的测量高度范围内的与各测量返回高度一一对应的第一风向和第一风速;第一信噪比数据集可以包括与各测量返回高度一一对应的第一信噪比A;第一湍流强度数据集包括与各测量返回高度一一对应的第一湍流强度B。需要说明的是,这里的各测量返回高度为相对于地面来说,多普勒测风激光雷达返回数据的高度相对于地面的高度,其具体数值由多普勒测风激光雷达的工作性能决定,本申请在此不作特别的限定。
可见,通过这种方式可以简单有效地获取用于输入大气颗粒物浓度反演模型的输入项。
作为一种优选的实施例,对第一信噪比数据集进行预处理后得到第二信噪比数据集,包括:
基于预设信噪比校正关系式对第一信噪比数据集进行校正后得到第二信噪比数据集,第二信噪比数据集包括与第一信噪比数据集中的各个第一信噪比A一一对应的校正后的各个第二信噪比C;
预设信噪比校正关系式为C=A*D*E;
其中,D为校正因子;E为与A对应的测量返回高度的平方。
本实施例中,考虑到对于多普勒测风激光雷达测量得到的第一信噪比数据集需要进行预处理,因此本申请中可以基于预设信噪比校正关系式对第一信噪比数据集进行校正以得到第二信噪比数据集。
需要说明的是,这里的校正因子与根据多普勒测风激光雷达的工作性能相关,本申请在此不作特别的限定。
可见,通过这种方式可以简单有效地对第一信噪比数据集进行预处理以得到表征标量守恒方程中的体源项和化学反应项的第二信噪比数据集。
作为一种优选的实施例,基于第一风数据集及多普勒测风激光雷达的预存的风速标准差数据集确定第一湍流强度数据集,包括:
根据多普勒测风激光雷达的预存的风速标准差数据集及预设湍流强度关系式对第一风数据集进行处理后确定第一湍流强度数据集;
预设湍流强度关系式为B=F/G,其中F为与各个B一一对应的各个第一风速,G为风速标准差数据集中的与各测量返回高度一一对应的风速标准差。
本申请中,为了得到第一湍流强度数据集,考虑到多普勒测风激光雷达中预存有风速标准差数据集,其中该风速标准差数据集为多普勒测风激光雷达根据其工作性能直接设定并预存,于是根据该风速标准差数据集及预设湍流强度关系式可以确定与各个第一风速一一对应的第一湍流强度B以构成第一湍流强度数据集,可见,通过这种方式可以简单可靠的得到用以表征标量守恒方程中的湍流扩散项的第一湍流强度数据集。
作为一种优选的实施例,将第一风数据集、第二信噪比数据集、第一湿度、第一能见度及第一湍流强度数据集作为输入项输入预先训练好的大气颗粒物浓度反演模型,包括:
对于第i个测量返回高度,其中1≤i≤测量返回高度的总个数且i为整数,执行如下步骤:
将第一能见度、第一湿度、与第i个测量返回高度对应的第一风向、第一风速、第二信噪比C、第一湍流强度B作为输入项输入预先训练好的大气颗粒物浓度反演模型以得到与第i个测量返回高度对应的输出项;
确定大气颗粒物浓度反演模型的输出项为第一预设时长内整层大气的第一大气颗粒物浓度,包括:
确定与第i个测量返回高度对应的输出项为第一预设时长内与第i个测量返回高度对应的第一大气颗粒物浓度。
本实施例中,考虑到第一风数据集中包括多普勒测风激光雷达的测量高度范围内的与各测量返回高度一一对应的第一风向和第一风速;第二信噪比数据集中包括与各测量返回高度一一对应的第二信噪比C;第一湍流强度数据集中包括与各测量返回高度一一对应的第一湍流强度B,于是在输入该大气颗粒物浓度反演模型时需要按照各测量返回高度进行逐高度层输入。具体来说,从第1个测量返回高度开始,将第一能见度、第一湿度、与第1个测量返回高度对应的第一风向、第一风速、第二信噪比C、第一湍流强度B作为输入项输入预先训练好的大气颗粒物浓度反演模型以得到与第1个测量返回高度对应的输出项;再针对2个测量返回高度,将第一能见度、第一湿度、与第2个测量返回高度对应的第一风向、第一风速、第二信噪比C、第一湍流强度B作为输入项输入预先训练好的大气颗粒物浓度反演模型以得到与第2个测量返回高度对应的输出项,依次类推,直到完成对各个测量返回高度的反演,由于多普勒测风激光雷达可以实时检测整层大气,因此最终得到了整层大气中的第一大气颗粒物浓度。
可见,通过这种方式可以简单可靠地依据观测数据和预先训练好的大气颗粒物浓度反演模型实现对整层大气的第一大气颗粒物浓度的观测。
作为一种优选的实施例,大气颗粒物浓度反演模型的预先训练过程,包括:
S11:获取地面环境监测站在各个第二预设时长内测得的整层大气的第二大气颗粒物浓度,其中,第二预设时长的个数为预设训练基准阈值;
S12:获取多普勒测风激光雷达测得的与第j个第二预设时长内的各个第三预设时长对应的第二风数据集、第三信噪比数据集及第二湿度,其中,第j个第二预设时长包括H个第三预设时长,H≥1且H为整数,1≤j≤预设训练基准阈值且j为整数,其中,第二风数据集包括多普勒测风激光雷达的测量高度范围内的与各测量返回高度一一对应的第二风向和第二风速;第三信噪比数据集包括与各测量返回高度一一对应的第三信噪比;
S13:获取地面气象站测得的与各个第二预设时长对应的第二能见度;
S14:判断在所有的第二湿度中是否存在不小于预设湿度阈值的第二湿度,若是,进入S15;若否,进入S18;
S15:判断在与各个不小于预设湿度阈值的第二湿度分别一一对应的各个第三信噪比数据集中,与第一预设高度对应的各个第三信噪比中是否包括多普勒测风激光雷达返回的无效值,若是,进入S16;若否,进入S17;
S16:将与各个包括无效值分别对应的各个第三预设时长内的第二风数据集、第三信噪比数据集、第二湿度及第二能见度删除;
S17:对与各个不包括无效值分别对应的各个第三信噪比数据集进行预处理后得到与各个第三信噪比数据集分别一一对应的各个第四信噪比数据集;对与各个不包括无效值分别对应的各个第二风数据集,分别基于各个第二风数据集及风速标准差数据集确定与各个第二风数据集一一对应的第二湍流强度数据集,其中,第四信噪比数据集包括与各测量返回高度一一对应的第四信噪比,第二湍流强度数据集包括与各测量返回高度一一对应的第二湍流强度;
S18:针对与各个小于预设湿度阈值的第二湿度分别对应的各个第二风数据集,分别基于各个第二风数据集及风速标准差数据集确定与各个第二风数据集一一对应的第二湍流强度数据集;对与各个小于预设湿度阈值的第二湿度分别对应的各个第三信噪比数据集,进行预处理后得到与各个第三信噪比数据集一一对应的各个第四信噪比数据集,其中,第四信噪比数据集包括与各测量返回高度一一对应的第四信噪比,第二湍流强度数据集包括与各测量返回高度一一对应的第二湍流强度;
S19:将各个第二预设时长及与各个第二预设时长一一对应的各个第二大气颗粒物浓度作为输出项输入至神经网络模型;
S20:针对与第j个第二预设时长对应的H个第三预设时长内的各第二风数据集、各第四信噪比数据集、各第二湿度、各第二湍流强度数据集,将各第二风数据集中的对应行的对应项进行加和并求平均值后得到第三风数据集,将各第四信噪比数据集的对应行的对应项进行加和并求平均值后得到第五信噪比数据集,将各第二湿度的对应行的对应项进行加和并求平均值后得到与第j个第二预设时长对应的第三湿度,将各第二湍流强度数据集的对应行的对应项进行加和并求平均值后得到第三湍流强度数据集;
S21:取第三风数据集中与第二预设高度对应的数据作为与第j个第二预设时长对应的第四风数据集;取第五信噪比数据集中与第二预设高度对应的数据作为与第j个第二预设时长对应的第六信噪比;取第三湍流强度数据集中与第二预设高度对应的数据作为与第j个第二预设时长对应的第四湍流强度;其中,第二预设高度为与各个测量返回高度中的最小值对应的高度;
S22:将各个第二预设时长、与各个第二预设时长分别对应的第四风数据集、第六信噪比、第三湿度、第二能见度及第四湍流强度作为输入项输入至神经网络模型以建立在各个第二预设时长下的输入-输出对应关系;
S23:利用神经网络模型进行第z次学习,z为大于1的整数;
S24:判断第z次学习的神经网络模型的学习精度是否满足预设精度阈值,若是,进入S26;若否,进入S25;
S25:令z=z+1并返回S22;
S26:停止学习并确定此时的神经网络模型为大气颗粒物浓度反演模型。
本申请中,利用神经网络模型来实现对大气颗粒物浓度反演模型的预先训练过程。
首先获取地面环境监测站测得的在各个第二预设时长内测得的整层大气的第二大气颗粒物浓度以作为待训练的大气颗粒物浓度反演模型的输出项;获取多普勒测风激光雷达测得的与第j个第二预设时长内的各个第三预设时长对应的第二风数据集、第三信噪比数据集及第二湿度;获取地面气象站测得的与各个第二预设时长对应的第二能见度。需要说明的是,首先,这里的第三预设时长与上述第一预设时长相同,即为多普勒测风激光雷达的当前返回数据的时间与上一次返回数据的时间之间的时间间隔,如10分钟,本申请在此不作特别的限定,根据具体的多普勒测风激光雷达的工作性能决定;且考虑到现有技术中,地面气象站和地面环境监测站每次测得数据之间的时间间隔是相同的,于是这里的第二预设时长可以为该时间间隔的时长,如1小时,于是第二预设时长包括6个第三预设时长,即H=6;同时,为了保证训练的效果,这里的第二预设时长的总个数需要足够多以保证用于训练的数据量足够多,该总个数即为这里的预设训练基准阈值,如在上述说明的基础上以3个月的数据进行训练则该预设训练基准阈值可以为2160,本申请在此不作特别的限定。
发明人进一步考虑到当发生雨雾等天气时,多普勒测风激光雷达返回的数据会受这些天气的影响而没办法准确反映以表征大气颗粒物浓度,因此,需要对该部分天气下测得的数据进行剔除。于是首先判断在所有的第二湿度中是否存在不小于预设湿度阈值的第二湿度,若否,则说明没有发生雨雾天气,则针对与各个小于预设湿度阈值的第二湿度分别对应的各个第二风数据集,分别基于各个第二风数据集及风速标准差数据集确定与各个第二风数据集一一对应的第二湍流强度数据集;对与各个小于预设湿度阈值的第二湿度分别对应的各个第三信噪比数据集,进行预处理后得到与各个第三信噪比数据集一一对应的各个第四信噪比数据集;若是,则说明可能发生了雨雾天气,但为了保证精确性,进一步考虑到当发生雨雾天气时通常多普勒测风激光雷达会在距离地面1公里左右返回的数据中包含无效值,因此继续判断在与各个不小于预设湿度阈值的第二湿度分别一一对应的各个第三信噪比数据集中,与第一预设高度对应的各个第三信噪比中是否包括多普勒测风激光雷达返回的无效值,这里的第一预设高度可以为上述所述的1公里,若是,说明确实发生了雨雾天气,于是将与各个包括无效值分别对应的各个第三预设时长内的第二风数据集、第三信噪比数据集、第二湿度及第二能见度删除;若否,说明虽然湿度达到了预设湿度阈值但是并没有发生雨雾天气,则对与各个不包括无效值分别对应的各个第三信噪比数据集进行预处理后得到与各个第三信噪比数据集分别一一对应的各个第四信噪比数据集;对与各个不包括无效值分别对应的各个第二风数据集,分别基于各个第二风数据集及风速标准差数据集确定与各个第二风数据集一一对应的第二湍流强度数据集。
随后针对与第j个第二预设时长对应的H个第三预设时长内的各第二风数据集、各第四信噪比数据集、各第二湿度、各第二湍流强度数据集,将各第二风数据集中的对应行的对应项进行加和并求平均值后得到第三风数据集,将各第四信噪比数据集的对应行的对应项进行加和并求平均值后得到第五信噪比数据集,将各第二湿度的对应行的对应项进行加和并求平均值后得到第三湿度,将各第二湍流强度数据集的对应行的对应项进行加和并求平均值后得到第三湍流强度数据集;
于是为了建立在各个第二预设时长下的输入-输出对应关系,取第三风数据集中与第二预设高度对应的数据作为与第j个第二预设时长对应的第四风数据集,其中这里的第二预设高度为与各个测量返回高度中的最小值对应的高度,即与返回的各层数据中最底层(距离地面最近的一层)对应的高度,如50米,本申请在此不作特别的限定,根据多普勒测风激光雷达的工作性能决定;取第五信噪比数据集中与第二预设高度对应的数据作为与第j个第二预设时长对应的第六信噪比;取第三湍流强度数据集中与第二预设高度对应的数据作为与第j个第二预设时长对应的第四湍流强度;随后将各个第二预设时长、与各个第二预设时长分别对应的第四风数据集、第六信噪比、第三湿度、第二能见度及第四湍流强度作为输入项输入至神经网络模型,即建立起了在各个第二预设时长下的输入-输出对应关系,需要说明的是,由于建立了各个第二预设时长下的输入-输出关系,若针对某个第二预设时长只有输出且由于数据的剔除过程而没有与之对应的输入时,则神经网络模型不会进行针对该第二预设时长的学习。
于是利用该神经网络模型进行第z次学习,并判断第z次学习的神经网络模型的学习精度是否满足预设阈值,若是,则说明达到训练的要求,停止学习并确定此时的神经网络模型为大气颗粒物浓度反演模型;若否,则说明还未达到训练的要求,令z=z+1并继续开始学习。需要说明的是,这里的预设阈值根据***经网络模型训练最终需要达到的学习精度设定,本申请在此不作特别的限定。
可见,通过这种方式可以可靠精确地建立大气颗粒物浓度反演模型,实现了从多普勒测风激光雷达及地面气象站测得的数据到大气颗粒物浓度的反演。
作为一种优选的实施例,神经网络模型包括BP神经网络、卷积神经网络、循环神经网络及长短期记忆神经网络中的一种或多种的组合。
本实施例中,用于得到大气颗粒物浓度反演模型的神经网络模型可以包括BP神经网络、卷积神经网络、循环神经网络及长短期记忆神经网络中的一种,也可以根据实际应用需要,该神经网络模型也可以包括BP神经网络、卷积神经网络、循环神经网络及长短期记忆神经网络中的多种的组合。
需要说明的是,这里仅列举了BP神经网络、卷积神经网络、循环神经网络及长短期记忆神经网络,在实际使用时也可以利用其他神经网络,本申请在此不作特别的限定,根据实际需要而定。
可见,通过这种方式可以简单可靠地得到用于获得大气颗粒物浓度的大气颗粒物浓度反演模型,且该模型的精确性和可靠性均较高。
请参照图2,图2为本发明提供的一种大气颗粒物浓度反演***的结构示意图。
该大气颗粒物浓度反演***,包括:
第一获取单元41,用于获取多普勒测风激光雷达在第一预设时长且在测量高度范围内测得的第一风数据集、第一信噪比数据集及第一湿度;
第二获取单元42,用于获取地面气象站在第一预设时长内测得的第一能见度;
第一确定单元43,用于基于第一风数据集及多普勒测风激光雷达的预存的风速标准差数据集确定第一湍流强度数据集;
第二确定单元44,用于对第一信噪比数据集进行预处理后得到第二信噪比数据集;
输入单元45,用于将第一风数据集、第二信噪比数据集、第一湿度、第一能见度及第一湍流强度数据集作为输入项输入预先训练好的大气颗粒物浓度反演模型;
输出确定单元46,用于确定大气颗粒物浓度反演模型的输出项为第一预设时长内整层大气的第一大气颗粒物浓度。
对于本发明中提供的大气颗粒物浓度反演***的介绍请参照上述大气颗粒物浓度反演方法的实施例,此处不再赘述。
请参照图3,图3为本发明提供的一种大气颗粒物浓度反演装置的结构示意图。
该大气颗粒物浓度反演装置,包括:
存储器51,用于存储计算机程序;
处理器52,用于执行计算机程序时实现如上述的大气颗粒物浓度反演方法的步骤。
对于本发明中提供的大气颗粒物浓度反演装置的介绍请参照上述大气颗粒物浓度反演方法的实施例,此处不再赘述。
本发明还提供了一种多普勒测风激光雷达,包括如上述的大气颗粒物浓度反演装置。
对于本发明中提供的多普勒测风激光雷达的介绍请参照上述大气颗粒物浓度反演方法的实施例,此处不再赘述。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
还需要说明的是,在本说明书中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
专业人员还可以进一步意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、计算机软件或者二者的结合来实现,为了清楚地说明硬件和软件的可互换性,在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本发明的范围。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其他实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (10)
1.一种大气颗粒物浓度反演方法,其特征在于,应用于多普勒测风激光雷达中的处理器,包括:
获取所述多普勒测风激光雷达在第一预设时长且在测量高度范围内测得的第一风数据集、第一信噪比数据集及第一湿度;
获取地面气象站在第一预设时长内测得的第一能见度;
基于所述第一风数据集及所述多普勒测风激光雷达的预存的风速标准差数据集确定第一湍流强度数据集;
对所述第一信噪比数据集进行预处理后得到第二信噪比数据集;
将所述第一风数据集、所述第二信噪比数据集、所述第一湿度、所述第一能见度及所述第一湍流强度数据集作为输入项输入预先训练好的大气颗粒物浓度反演模型;
确定所述大气颗粒物浓度反演模型的输出项为所述第一预设时长内整层大气的第一大气颗粒物浓度。
2.如权利要求1所述的大气颗粒物浓度反演方法,其特征在于,所述第一风数据集包括所述多普勒测风激光雷达的所述测量高度范围内的与各测量返回高度一一对应的第一风向和第一风速;
所述第一信噪比数据集包括与各所述测量返回高度一一对应的第一信噪比A;
所述第一湍流强度数据集包括与各所述测量返回高度一一对应的第一湍流强度B。
3.如权利要求2所述的大气颗粒物浓度反演方法,其特征在于,对所述第一信噪比数据集进行预处理后得到第二信噪比数据集,包括:
基于预设信噪比校正关系式对所述第一信噪比数据集进行校正后得到第二信噪比数据集,所述第二信噪比数据集包括与所述第一信噪比数据集中的各个第一信噪比A一一对应的校正后的各个第二信噪比C;
所述预设信噪比校正关系式为C=A*D*E;
其中,D为校正因子;E为与所述A对应的所述测量返回高度的平方。
4.如权利要求2所述的大气颗粒物浓度反演方法,其特征在于,基于所述第一风数据集及所述多普勒测风激光雷达的预存的风速标准差数据集确定第一湍流强度数据集,包括:
根据所述多普勒测风激光雷达的预存的风速标准差数据集及预设湍流强度关系式对所述第一风数据集进行处理后确定第一湍流强度数据集;
所述预设湍流强度关系式为B=F/G,其中F为与各个所述B一一对应的各个所述第一风速,G为所述风速标准差数据集中的与各所述测量返回高度一一对应的风速标准差。
5.如权利要求3所述的大气颗粒物浓度反演方法,其特征在于,将所述第一风数据集、所述第二信噪比数据集、所述第一湿度、所述第一能见度及所述第一湍流强度数据集作为输入项输入预先训练好的大气颗粒物浓度反演模型,包括:
对于第i个所述测量返回高度,其中1≤i≤所述测量返回高度的总个数且i为整数,执行如下步骤:
将所述第一能见度、所述第一湿度、与第i个所述测量返回高度对应的所述第一风向、所述第一风速、所述第二信噪比C、所述第一湍流强度B作为输入项输入预先训练好的大气颗粒物浓度反演模型以得到与第i个所述测量返回高度对应的输出项;
确定所述大气颗粒物浓度反演模型的输出项为所述第一预设时长内整层大气的第一大气颗粒物浓度,包括:
确定与第i个所述测量返回高度对应的输出项为所述第一预设时长内与第i个所述测量返回高度对应的第一大气颗粒物浓度。
6.如权利要求1至5任一项所述的大气颗粒物浓度反演方法,其特征在于,所述大气颗粒物浓度反演模型的预先训练过程,包括:
S11:获取所述地面环境监测站在各个第二预设时长内测得的整层大气的第二大气颗粒物浓度,其中,所述第二预设时长的个数为预设训练基准阈值;
S12:获取所述多普勒测风激光雷达测得的与第j个第二预设时长内的各个第三预设时长对应的第二风数据集、第三信噪比数据集及第二湿度,其中,第j个所述第二预设时长包括H个所述第三预设时长,H≥1且H为整数,1≤j≤所述预设训练基准阈值且j为整数,其中,所述第二风数据集包括所述多普勒测风激光雷达的所述测量高度范围内的与各所述测量返回高度一一对应的第二风向和第二风速;所述第三信噪比数据集包括与各所述测量返回高度一一对应的第三信噪比;
S13:获取所述地面气象站测得的与各个所述第二预设时长对应的第二能见度;
S14:判断在所有的所述第二湿度中是否存在不小于预设湿度阈值的第二湿度,若是,进入S15;若否,进入S18;
S15:判断在与各个不小于所述预设湿度阈值的第二湿度分别一一对应的各个所述第三信噪比数据集中,与第一预设高度对应的各个第三信噪比中是否包括所述多普勒测风激光雷达返回的无效值,若是,进入S16;若否,进入S17;
S16:将与各个包括所述无效值分别对应的各个所述第三预设时长内的第二风数据集、第三信噪比数据集、第二湿度及第二能见度删除;
S17:对与各个不包括所述无效值分别对应的各个所述第三信噪比数据集进行预处理后得到与各个所述第三信噪比数据集分别一一对应的各个第四信噪比数据集;对与各个不包括所述无效值分别对应的各个所述第二风数据集,分别基于各个所述第二风数据集及所述风速标准差数据集确定与各个所述第二风数据集一一对应的第二湍流强度数据集,其中,所述第四信噪比数据集包括与各所述测量返回高度一一对应的第四信噪比,所述第二湍流强度数据集包括与各所述测量返回高度一一对应的第二湍流强度;
S18:针对与各个小于所述预设湿度阈值的第二湿度分别对应的各个所述第二风数据集,分别基于各个所述第二风数据集及所述风速标准差数据集确定与各个所述第二风数据集一一对应的第二湍流强度数据集;对与各个小于所述预设湿度阈值的第二湿度分别对应的各个所述第三信噪比数据集,进行预处理后得到与各个所述第三信噪比数据集一一对应的各个第四信噪比数据集,其中,所述第四信噪比数据集包括与各所述测量返回高度一一对应的第四信噪比,所述第二湍流强度数据集包括与各所述测量返回高度一一对应的第二湍流强度;
S19:将各个所述第二预设时长及与各个所述第二预设时长一一对应的各个所述第二大气颗粒物浓度作为输出项输入至神经网络模型;S20:针对与第j个所述第二预设时长对应的H个第三预设时长内的各所述第二风数据集、各所述第四信噪比数据集、各所述第二湿度、各所述第二湍流强度数据集,将各所述第二风数据集中的对应行的对应项进行加和并求平均值后得到第三风数据集,将各所述第四信噪比数据集的对应行的对应项进行加和并求平均值后得到第五信噪比数据集,将各所述第二湿度的对应行的对应项进行加和并求平均值后得到与第j个所述第二预设时长对应的第三湿度,将各所述第二湍流强度数据集的对应行的对应项进行加和并求平均值后得到第三湍流强度数据集;
S21:取所述第三风数据集中与第二预设高度对应的数据作为与第j个所述第二预设时长对应的第四风数据集;取所述第五信噪比数据集中与所述第二预设高度对应的数据作为与第j个所述第二预设时长对应的第六信噪比;取所述第三湍流强度数据集中与所述第二预设高度对应的数据作为与第j个所述第二预设时长对应的第四湍流强度;其中,所述第二预设高度为与各个所述测量返回高度中的最小值对应的高度;S22:将各个所述第二预设时长、与各个所述第二预设时长分别对应的所述第四风数据集、所述第六信噪比、所述第三湿度、所述第二能见度及所述第四湍流强度作为输入项输入至所述神经网络模型以建立在各个所述第二预设时长下的输入-输出对应关系;
S23:利用所述神经网络模型进行第z次学习,z为大于1的整数;
S24:判断第z次学习的所述神经网络模型的学习精度是否满足预设精度阈值,若是,进入S26;若否,进入S25;
S25:令z=z+1并返回S22;
S26:停止学习并确定此时的所述神经网络模型为所述大气颗粒物浓度反演模型。
7.如权利要求6所述的大气颗粒物浓度反演方法,其特征在于,所述神经网络模型包括BP神经网络、卷积神经网络、循环神经网络及长短期记忆神经网络中的一种或多种的组合。
8.一种大气颗粒物浓度反演***,其特征在于,包括:
第一获取单元,用于获取所述多普勒测风激光雷达在第一预设时长且在测量高度范围内测得的第一风数据集、第一信噪比数据集及第一湿度;
第二获取单元,用于获取地面气象站在第一预设时长内测得的第一能见度;
第一确定单元,用于基于所述第一风数据集及所述多普勒测风激光雷达的预存的风速标准差数据集确定第一湍流强度数据集;
第二确定单元,用于对所述第一信噪比数据集进行预处理后得到第二信噪比数据集;
输入单元,用于将所述第一风数据集、所述第二信噪比数据集、所述第一湿度、所述第一能见度及所述第一湍流强度数据集作为输入项输入预先训练好的大气颗粒物浓度反演模型;
输出确定单元,用于确定所述大气颗粒物浓度反演模型的输出项为所述第一预设时长内整层大气的第一大气颗粒物浓度。
9.一种大气颗粒物浓度反演装置,其特征在于,包括:
存储器,用于存储计算机程序;
处理器,用于执行所述计算机程序时实现如权利要求1至7任一项所述的大气颗粒物浓度反演方法的步骤。
10.一种多普勒测风激光雷达,其特征在于,包括如权利要求9所述的大气颗粒物浓度反演装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111604123.2A CN114279915A (zh) | 2021-12-24 | 2021-12-24 | 一种大气颗粒物浓度反演方法及相关组件 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111604123.2A CN114279915A (zh) | 2021-12-24 | 2021-12-24 | 一种大气颗粒物浓度反演方法及相关组件 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114279915A true CN114279915A (zh) | 2022-04-05 |
Family
ID=80875440
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111604123.2A Pending CN114279915A (zh) | 2021-12-24 | 2021-12-24 | 一种大气颗粒物浓度反演方法及相关组件 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114279915A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114859321A (zh) * | 2022-07-11 | 2022-08-05 | 青岛镭测创芯科技有限公司 | 区域大气颗粒物来源比例计算方法、装置、设备、介质 |
CN116359877A (zh) * | 2023-05-30 | 2023-06-30 | 航天宏图信息技术股份有限公司 | 一种基于风廓线雷达的大雾天气监测方法和装置 |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4940327A (en) * | 1988-10-25 | 1990-07-10 | Trc Companies Inc. | Method and apparatus for real time asbestos aerosol monitoring |
US5922976A (en) * | 1995-10-12 | 1999-07-13 | California Institute Of Technology | Method of measuring aerosol particles using automated mobility-classified aerosol detector |
US20120050750A1 (en) * | 2009-04-21 | 2012-03-01 | Michigan Aerospace Corporation | Atmospheric measurement system |
US20140157872A1 (en) * | 2011-08-04 | 2014-06-12 | Cambidge Enterptise Limited | Sensing systems |
CN106596362A (zh) * | 2016-12-14 | 2017-04-26 | 北京市环境保护监测中心 | 一种激光雷达输送通量计算方法和装置 |
US20170241893A1 (en) * | 2016-02-19 | 2017-08-24 | Research Triangle Institute | Devices, systems and methods for detecting particles |
US20180292826A1 (en) * | 2017-04-11 | 2018-10-11 | Defelice Thomas Peter | Intelligent systems for weather modification programs |
CN108763756A (zh) * | 2018-05-28 | 2018-11-06 | 河南工业大学 | 一种气溶胶光学厚度与pm2.5反演订正方法及其*** |
CN109145962A (zh) * | 2018-07-31 | 2019-01-04 | 南京信息工程大学 | 一种基于数字图像的大气参数反演观测方法 |
US20190293539A1 (en) * | 2015-03-06 | 2019-09-26 | Scanit Technologies, Inc. | Airborne particle monitor |
CN111398993A (zh) * | 2020-05-20 | 2020-07-10 | 青岛镭测创芯科技有限公司 | 一种无盲区脉冲相干测风激光雷达*** |
CN113138398A (zh) * | 2020-01-17 | 2021-07-20 | 中国海洋大学 | 一种基于无人机大气激光雷达的气溶胶消光系数反演方法 |
CN113422643A (zh) * | 2021-02-10 | 2021-09-21 | 淮南师范学院 | 基于气象大数据预测的混合fso/rf自适应切换通信方法 |
CN113657028A (zh) * | 2021-08-05 | 2021-11-16 | 长春理工大学 | 一种基于多源信息的气溶胶光学厚度在线预测方法 |
-
2021
- 2021-12-24 CN CN202111604123.2A patent/CN114279915A/zh active Pending
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4940327A (en) * | 1988-10-25 | 1990-07-10 | Trc Companies Inc. | Method and apparatus for real time asbestos aerosol monitoring |
US5922976A (en) * | 1995-10-12 | 1999-07-13 | California Institute Of Technology | Method of measuring aerosol particles using automated mobility-classified aerosol detector |
US20120050750A1 (en) * | 2009-04-21 | 2012-03-01 | Michigan Aerospace Corporation | Atmospheric measurement system |
US20140157872A1 (en) * | 2011-08-04 | 2014-06-12 | Cambidge Enterptise Limited | Sensing systems |
US20190293539A1 (en) * | 2015-03-06 | 2019-09-26 | Scanit Technologies, Inc. | Airborne particle monitor |
US20170241893A1 (en) * | 2016-02-19 | 2017-08-24 | Research Triangle Institute | Devices, systems and methods for detecting particles |
CN106596362A (zh) * | 2016-12-14 | 2017-04-26 | 北京市环境保护监测中心 | 一种激光雷达输送通量计算方法和装置 |
US20180292826A1 (en) * | 2017-04-11 | 2018-10-11 | Defelice Thomas Peter | Intelligent systems for weather modification programs |
CN108763756A (zh) * | 2018-05-28 | 2018-11-06 | 河南工业大学 | 一种气溶胶光学厚度与pm2.5反演订正方法及其*** |
CN109145962A (zh) * | 2018-07-31 | 2019-01-04 | 南京信息工程大学 | 一种基于数字图像的大气参数反演观测方法 |
CN113138398A (zh) * | 2020-01-17 | 2021-07-20 | 中国海洋大学 | 一种基于无人机大气激光雷达的气溶胶消光系数反演方法 |
CN111398993A (zh) * | 2020-05-20 | 2020-07-10 | 青岛镭测创芯科技有限公司 | 一种无盲区脉冲相干测风激光雷达*** |
CN113422643A (zh) * | 2021-02-10 | 2021-09-21 | 淮南师范学院 | 基于气象大数据预测的混合fso/rf自适应切换通信方法 |
CN113657028A (zh) * | 2021-08-05 | 2021-11-16 | 长春理工大学 | 一种基于多源信息的气溶胶光学厚度在线预测方法 |
Non-Patent Citations (7)
Title |
---|
IWONA S. STACHLEWSKA 等: "Effect of Heat Wave Conditions on Aerosol Optical Properties Derived from Satellite and Ground-Based Remote Sensing over Poland", 《REMOTE SENS.》, vol. 9, no. 11, 31 December 2017 (2017-12-31), pages 1199 * |
THOMAS FAHEY 等: "Laser Beam Atmospheric Propagation Modelling for Aerospace LIDAR Applications", 《ATMOSPHERE》, vol. 12, no. 7, 17 July 2021 (2021-07-17), pages 918 * |
XIAOYE WANG 等: "Retrieval and Calculation of Vertical Aerosol Mass Fluxes by a Coherent Doppler Lidar and a Sun Photometer", 《REMOTE SENS.》, vol. 13, no. 16, 18 August 2021 (2021-08-18), pages 3259 * |
付松琳: "星载气溶胶激光雷达的模拟仿真和反演算法应用研究", 《中国博士学位论文全文数据库信息科技辑》, no. 9, 15 September 2021 (2021-09-15), pages 136 - 39 * |
刘松: "基于Mie散射激光雷达的南京仙林地区气溶胶消光特性及PM2.5质量浓度估算模型研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》, no. 1, 15 January 2019 (2019-01-15), pages 027 - 1056 * |
唐新姿 等: "《风电技术》", vol. 1, 31 August 2017, 湘潭大学出版社, pages: 19 * |
靳翔 等: "基于多普勒激光雷达的边界层内湍流参数估算", 《中国激光》, vol. 48, no. 11, 30 June 2021 (2021-06-30), pages 170 - 179 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114859321A (zh) * | 2022-07-11 | 2022-08-05 | 青岛镭测创芯科技有限公司 | 区域大气颗粒物来源比例计算方法、装置、设备、介质 |
CN116359877A (zh) * | 2023-05-30 | 2023-06-30 | 航天宏图信息技术股份有限公司 | 一种基于风廓线雷达的大雾天气监测方法和装置 |
CN116359877B (zh) * | 2023-05-30 | 2023-08-15 | 航天宏图信息技术股份有限公司 | 一种基于风廓线雷达的大雾天气监测方法和装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10830743B2 (en) | Determining the net emissions of air pollutants | |
CN114279915A (zh) | 一种大气颗粒物浓度反演方法及相关组件 | |
US20180197499A1 (en) | Display status adjustment method, display status adjustment device and display device | |
CN109709577A (zh) | 一种基于WRF-Chem模式的气溶胶激光雷达反演PM2.5的三维变分同化方法 | |
CN110058258B (zh) | 一种基于混合型激光雷达的大气边界层探测方法 | |
Rodrigo et al. | Multi-site testing and evaluation of remote sensing instruments for wind energy applications | |
CN111473755B (zh) | 一种远程测距方法及装置 | |
CN116297068B (zh) | 基于地表反射率优化的气溶胶光学厚度反演方法及*** | |
CN110738354A (zh) | 预测颗粒物浓度的方法、装置、存储介质及电子设备 | |
CN113189014A (zh) | 一种融合卫星遥感和地面监测数据的臭氧浓度估算方法 | |
CN112711002B (zh) | 一种基于co2-dial模拟测量的点源co2排放的新型估算方法 | |
CN110595960A (zh) | 一种基于机器学习的pm2.5浓度遥感估算方法 | |
CN112364940A (zh) | 基于多源数据的大气污染物源解析方法、存储介质及设备 | |
CN114662344A (zh) | 基于连续在线观测数据的大气污染溯源预测方法及*** | |
CN116205541A (zh) | 评估本地污染源对环境空气质量影响的方法及装置 | |
CN113360587A (zh) | 一种基于gis技术的土地测绘设备及测绘方法 | |
Zhou et al. | Mobile sensing of point-source gas emissions using Bayesian inference: An empirical examination of the likelihood function | |
CN116626238A (zh) | 一种气流和气体组合感知的双通道检测罗盘***、气体泄漏检测方法、数据融合及追踪方法 | |
Weber | Estimators for the standard deviation of horizontal wind direction | |
CN114492984A (zh) | 粉尘浓度的时空分布预测方法、装置、设备和存储介质 | |
CN116258101B (zh) | 一种基于气体浓度监测数据快速计算排放源强的方法 | |
CN117131441A (zh) | 夜间光污染监测方法、装置、计算机设备和存储介质 | |
CN115639625B (zh) | 一种高空间分辨率的精细化大气水平能见度自动预报*** | |
CN115372287B (zh) | 基于高分五号卫星数据的二氧化硫排放清单构建方法 | |
CN116466368A (zh) | 基于激光雷达和卫星资料的沙尘消光系数廓线估算方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |