买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

龙图腾网恭喜中国科学院合肥物质科学研究院申请的专利基于地面遥感测量的点源排放强度预测方法获国家发明授权专利权，本发明授权专利权由国家知识产权局授予，授权公告号为：CN115480032B 。

龙图腾网通过国家知识产权局官网在2025-04-18发布的发明授权授权公告中获悉：该发明授权的专利申请号/专利号为：202211110049.3，技术领域涉及：G01N33/00；该发明授权基于地面遥感测量的点源排放强度预测方法是由叶函函;王先华;施海亮;熊伟;吴时超设计研发完成，并于2022-09-13向国家知识产权局提交的专利申请。

本基于地面遥感测量的点源排放强度预测方法在说明书摘要公布了：本发明提供一种基于地面遥感测量的点源排放强度预测方法。在烟羽测量点位和背景测量点位测量到达地表的太阳光谱，反演获得点源排放烟羽中特定污染气体在大气中的测量浓度值。利用高斯扩散模型，仿真计算特定污染气体的排放烟羽3维结构，获得烟羽测量点位特定污染气体的仿真浓度值。利用测量浓度值和仿真浓度值进行反演，计算出点源的排放强度和风速。根据反演的排放强度和风速，重新计算高斯扩散模型，重构出烟羽3维结构。本发明采用便携的地面遥感设备在点源周边开展烟羽的柱浓度测量，利用多角度测量数据开展烟羽3维结构的反演重构，不仅能够同步计算出排放强度，还能够降低传统排放强度计算方法中风场误差的影响。

本发明授权基于地面遥感测量的点源排放强度预测方法在权利要求书中公布了：1.一种基于地面遥感测量的点源排放强度预测方法，其特征在于，包括：步骤S1，利用气象仪测量获得排放污染气体的点源的周边区域的气象数据，利用测距仪测量获得所述点源的高度，利用定位仪测量获得所述点源的经纬度，根据所述气象数据中的风向和风速来估算所述点源的排放烟羽的范围，在太阳光线穿过所述点源的排放烟羽在地面上形成的烟羽投影区域中选择多个烟羽测量点位，在所述烟羽投影区域以外的地面上选择至少1个背景测量点位，对所述多个烟羽测量点位和所述至少1个背景测量点位中的各个测量点位，利用光谱仪测量获得到达地表的太阳光谱的光谱数据并记录对应的测量时刻，利用定位仪测量获得所述测量点位的经纬度，利用测距仪测量获得所述测量点位与所述点源的距离；步骤S2，对步骤S1中的每个所述测量点位，利用该测量点位的所述经纬度和所述测量时刻，计算出该测量点位的在所述测量时刻的太阳高度角和太阳方位角；步骤S3，利用步骤S1中获得的所述各个测量点位的光谱数据、同步获得的气象数据、测量时刻以及步骤S2中计算出的太阳高度角和太阳方位角，根据光谱反演算法进行反演计算，计算出在所述各个测量点位测得的所述点源排放的污染气体中的特定污染气体在大气中的测量浓度值；步骤S4，以从所述点源沿铅垂方向向地面引的垂线与地面的交点为原点、以所述气象数据中的风向为x轴方向、以所述点源的高度方向即铅垂方向为z轴方向、以与所述x轴方向和所述z轴方向垂直的方向为y轴方向，建立3维坐标系，根据步骤S1中获得的所述气象数据中的风速与作为先验排放数据的所述点源的排放时间、高度和排放强度，在所述3维坐标系中，利用高斯扩散模型进行仿真获得所述排放烟羽的3维仿真结构；步骤S5，在步骤S4中建立的所述3维坐标系中，根据步骤S1中获得的所述多个烟羽测量点位的经纬度相对于所述点源的经纬度的差值以及所述多个测量点位与所述点源的距离，定位出步骤S1中的所述多个烟羽测量点位，对每个所述烟羽测量点位，根据其测量时刻的太阳高度角和太阳方位角来构建测量光路，计算出该测量光路的穿烟羽部分的所述特定污染气体的浓度积分值，获得所述多个烟羽测量点位的所述特定污染气体的仿真浓度值；和步骤S6，利用步骤S3中反演计算出的所述特定污染气体的所述测量浓度值和步骤S5中仿真计算出的与所述测量浓度值对应的所述仿真浓度值进行反演，计算出作为反演结果的所述点源的反演排放强度和反演风速，在步骤S2中，太阳高度角h0的计算公式如下：h0＝sin-1[sinψsinσ+cosψcosσcos15t0+λ-300]式中ψ为纬度，单位是度，λ为经度，单位是度，t0为观测时刻，单位是时，σ为太阳方位角，单位是度，h0的单位是度，太阳方位角的计算公式如下： 式中θ0＝360dn365，单位是弧度，其中dn为所选日期在本年度中的次序，取值范围为0至364，在步骤S3中，对步骤S1中获得的所述多个烟羽测量点位和所述至少1个背景测量点位的所述光谱数据，采用最优估算法，利用Tikhonov-Philips约束条件，结合非线性逐次迭代，构建出反演计算公式，将步骤S1中获得的所述各个测量点位的所述光谱数据转换成对数光谱数据，与同步获得的气象数据、测量时刻以及步骤S2中计算出的太阳高度角和太阳方位角一并输入所述反演计算公式，计算出所述特定污染气体的柱浓度，将与各个所述烟羽测量点位对应的各个所述特定污染气体的柱浓度分别减去与所述至少1个背景测量点位对应的所述特定污染气体的柱浓度的平均值，获得与各个所述烟羽测量点位对应的所述特定污染气体的所述测量浓度值，在步骤S4中，当设所述风速为μ、所述排放时间为t、所述高度为H、所述排放强度为Q时，计算所述排放烟羽的3维仿真结构的所述高斯扩散模型由下述的1式表示： 其中，x,y,z是所述3维坐标系中的烟羽位置的坐标，Cx,y,z为所述3维坐标系中x,y,z位置处的所述特定污染气体的浓度，σx、σy、σz分别为x轴方向、y轴方向、z轴方向上的扩散系数，所述风速μ的单位是米秒，所述排放时间t的单位是秒，所述高度H的单位是米，所述排放强度Q的单位是gs，所述浓度Cx,y,z的单位是gm2，所述扩散系数σx、σy、σz是根据与大气稳定度等级相应的经验公式计算出的值，单位是米，在步骤S6中，在进行反演时，利用步骤S3中反演计算出的所述特定污染气体的所述测量浓度值和步骤S5中仿真计算出的与所述测量浓度值对应的所述仿真浓度值，采用最优估算法，利用Gauss-Newton迭代反演出所述特定污染气体的排放强度和风速，求解公式如下：ui+1＝ua+KiTSε-1Ki+Sa-1-1KiTSε-1[v-Fui-Kiui-ua]2式中u表示所述特定污染气体的排放强度和风速，下标i表示第i次迭代，ui+1为第i+1次迭代结果，ua为所述特定污染气体的排放强度和风速的初始猜值，为权重函数，T表示矩阵转置，Sε为所述特定污染气体测量误差协方差，Sa为所述特定污染气体的排放强度和风速的先验误差协方差矩阵，v表示所述多个烟羽测量点位的所述特定污染气体的测量浓度值，Fu表示所述多个烟羽测量点位的所述特定污染气体的仿真浓度值，在迭代计算过程中，当u的更新值小于初始猜值的万分之一时，认为迭代收敛，此时ui+1的值即为作为所述反演结果的所述点源的反演排放强度和反演风速。

如需购买、转让、实施、许可或投资类似专利技术，可联系本专利的申请人或专利权人中国科学院合肥物质科学研究院，其通讯地址为：230031 安徽省合肥市蜀山湖路350号；或者联系龙图腾网官方客服，联系龙图腾网可拨打电话0551-65771310或微信搜索“龙图腾网”。