CN110031412A - 基于移动ahsa观测的大气污染物排放通量获取方法 - Google Patents

基于移动ahsa观测的大气污染物排放通量获取方法 Download PDF

Info

Publication number
CN110031412A
CN110031412A CN201910338227.XA CN201910338227A CN110031412A CN 110031412 A CN110031412 A CN 110031412A CN 201910338227 A CN201910338227 A CN 201910338227A CN 110031412 A CN110031412 A CN 110031412A
Authority
CN
China
Prior art keywords
observation
ahsa
mobile
emission flux
atmosphere
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
Application number
CN201910338227.XA
Other languages
English (en)
Inventor
刘诚
谈伟
朱一芝
邢成志
刘浩然
季祥光
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Science and Technology of China USTC
Original Assignee
University of Science and Technology of China USTC
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Science and Technology of China USTC filed Critical University of Science and Technology of China USTC
Priority to CN201910338227.XA priority Critical patent/CN110031412A/zh
Publication of CN110031412A publication Critical patent/CN110031412A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry

Landscapes

  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明公开了一种基于移动AHSA观测的大气污染物排放通量获取方法,利用移动AHSA观测获取空间分布的大气污染物(二氧化氮、二氧化硫、甲醛等)的差分斜柱浓度结果,利用气象化学模型模拟获取的大气信息,输入到大气辐射传输中计算出大气光子路径信息,进一步转化差分斜柱浓度结果为对流层垂直柱浓度结果,结合气象化学模型模拟的风场信息计算获取污染源的排放通量结果。该方法使用遥感方法间接、快速、准确量化大气污染物排放通量结果,为大气环境污染提供了有力的监测手段,同时,大气污染物排放通量计算结果,为大气环境治理提供必要的数据支持。

Description

基于移动AHSA观测的大气污染物排放通量获取方法
技术领域
本发明涉及环境科学、地基遥感技术领域,尤其涉及一种基于移动AHSA观测的大气污染物排放通量获取方法。
背景技术
随着经济、社会的快速发展,环境问题日益加重,如何解决环境问题受到国际社会的普遍关注。环境污染的有力监测,可以为环境治理提供必要的数据支持。
目前,我国在区域大气环境污染监测方面,主要是以地面采样为主,辅以适当的探空观测等对区域大气环境的污染状况进行观测研究,缺少对大气污染物排放通量的获取以及有针对性的对大气污染源排放特征的定量分析研究。
移动AHSA(边界层大气成分高光谱扫描与分析仪)观测作为地基遥感的一种,相对于传统的地面采样分析技术具有其快速、准确、信息量大、针对性强等特点在大气环境监测领域具有独特的优势。同时,我国又迫切需要获取城市、化工厂、电厂等重大污染源的污染物排放通量结果,并且污染源排放数据可以同化在大气化学模型数值模拟研究中,以便更准确的呈现大气混合演化过程,为大气环境治理提供有力的数据支持。
然而,目前还没有较为准确的用于大气污染物排放通量的计算方案。
发明内容
本发明的目的是提供一种基于移动AHSA观测的大气污染物排放通量获取方法,用于解决对大气污染物排放通量的量化难题,拓展地基AHSA的应用范围(搭载在移动平台上),实现对主要大气污染物(二氧化氮、二氧化硫、甲醛等)排放的快速、准确量化,以及对主要污染源的有效监测,满足我国大气污染气体的排放总量监测,为大气污染源排放特征量化研究提供关键的数据支持。
本发明的目的是通过以下技术方案实现的:
一种基于移动AHSA观测的大气污染物排放通量获取方法,包括:
基于差分吸收光谱解析方法,使用光谱分析算法反演分析移动AHSA观测的散射光谱数据,获取大气污染物的差分斜柱浓度;
使用气象化学模型模拟得到大气温度、压力廓线,气溶胶以及相关大气污染物的廓线信息;
将模拟得到大气温度、压力廓线,气溶胶以及相关大气污染物的廓线信息输入到大气辐射传输模型中计算获得大气光子路径结果;
使用大气光子路径结果将大气污染物的差分斜柱浓度转化为对流层垂直柱浓度;
结合对流层垂直柱浓度与预先获得的风场信息,采用环路积分方法计算获得大气污染物的排放通量。
由上述本发明提供的技术方案可以看出,其主要获得如下有益效果:
1)光谱拟合反演过程中采用固定参考谱设置,扣除相同的平流层背景信息,有效避免不同观测位置处不同背景值对污染物浓度的影响。
2)基于移动AHSA观测,获取主要大气污染气体(二氧化氮、二氧化硫、甲醛等)的柱浓度的空间分布结果,能够直接呈现污染物的时空分布信息。
3)采用大气辐射传输模型计算大气光子路径信息,较简单的几何近似方法,充分考虑了实际观测的大气信息(温度、压力廓线,气溶胶信息、气体信息等)、以及观测的几何信息(太阳天顶角,观测方位角等),能够获得更为精确的大气光子路径结果,进一步获得高质量的大气污染物垂直柱浓度结果。
4)采用模型模拟获得的高分辨的风场信息,较有限的观测,能够提供更高的空间分辨率(水平方向最高可实现1km×1km,垂直方向高度可达到大气压为5000pa处对应高度)、时间分辨率(最高可实现10min),弥补实际观测时间空间分辨率的不足(实测数据可以用来验证模型模拟结果,确保模型模拟的准确性、可靠性)。
5)建立基于地基遥感获取大气污染物排放通量的方案,实现快速、准确量化大气污染物排放通量,以及对大气污染源排放的有效监测。满足我国大气污染气体的排放总量监测,为大气污染源排放特征量化研究提供关键的数据支持。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域的普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他附图。
图1为本发明实施例提供的一种基于移动AHSA观测的大气污染物排放通量获取方法的流程图。
具体实施方式
下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明的保护范围。
本发明实施例提供一种基于移动AHSA观测的大气污染物排放通量获取方法,其主要包括如下步骤:
步骤一、基于差分吸收光谱解析方法,使用光谱分析算法反演分析移动AHSA观测的散射光谱数据,获取大气污染物的差分斜柱浓度。
本发明实施例中,上述大气污染物主要包括二氧化氮、二氧化硫、甲醛等。
本发明实施例中,所述移动AHSA观测的散射光谱数据为多个角度的散射光谱数据,示例性的,可以为22°或者30°(测量谱),90°(参考谱)的太阳散射光谱。
在实现过程中,可以利用差分吸收光谱解析方法,考虑不同大气污染物的特征吸收结构定性定量分析获得二氧化氮、二氧化硫、甲醛等大气污染物差分斜柱浓度结果。反演过程基于差分吸收光谱解析方法,使用光谱分析算法对测得的太阳散射光谱进行反演,计算二氧化氮、二氧化硫、甲醛等痕量气体的差分斜柱浓度结果。同时在光谱反演之前要扣除仪器暗电流的影响以及矫正光谱偏移。具体的是使用采集的90°光谱作为夫琅禾费参考谱,结合测量谱反演获得不同气体的差分斜柱浓度结果,反演过程中根据拟合的残差以及拟合误差量级判定拟合的质量,获得高精度的差分斜柱浓度结果。
步骤二、使用气象化学模型模拟得到大气温度、压力廓线,气溶胶、以及相关大气污染物的廓线信息。
示例性的,气象化学模型模拟过程中可使用美国国家环境预测中心通过全球预测模型得到的数据(0.5°的空间分辨率和6小时的时间分辨率)作为气象场模拟的初始和边界条件。在模拟大气化学过程中,综合考虑光化学机制与气溶胶相互作用等化学过程。人为排放输入数据可来源于高分辨率的中国排放清单数据。根据研究需要可以获得不同水平、垂直分辨率的模型数据。
步骤三、将模拟得到大气温度、压力廓线,气溶胶以及相关大气污染物的廓线信息输入到大气辐射传输模型中计算获得大气光子路径结果。
大气辐射传输模型能够计算权重函数、大气光子路径、斜柱浓度、垂直光学厚度等多种参量。大气辐射传输模型计算过程考虑了Rayleigh散射、Mie散射、地表反照率和大气折射,可以计算出大气中任意高度多种观测模式的光谱模拟。大气辐射传输模型的模拟输入参数包括了气溶胶参数、压强温度廓线,地面反照率,云等大气模拟参数。
步骤四、使用大气光子路径结果将大气污染物的差分斜柱浓度转化为对流层垂直柱浓度。
本发明实施例中,可通过如下公式将大气污染物的差分斜柱浓度转化为对流层垂直柱浓度:
其中,VCD为大气污染物的垂直柱浓度,DSCDtrop为对流层大气污染物的差分斜柱浓度,AMF为大气光子路径,α为观测仰角。
步骤五、结合对流层垂直柱浓度与风场信息,采用环路积分方法计算获得大气污染物的排放通量。
本发明实施例中,所述预先获得的风场信息为实测的风场信息,或者通过气象化学模型模拟得到的风场信息,图1所示流程以后者为例。
本发明实施例中,采用环路积分方法计算获得大气污染物的排放通量的公式如下:
Flux=∮VCD·V·ds;
将环路积分转化为有限的观测片段累和,则大气污染物的排放通量具体计算过程为:
Flux=∑VCDi·Vi·sinθi·ΔSi
上式中,∮为环路积分符号,d表示增量;V代表观风矢量在行驶方向垂直方向的分量,VCDi第i次观测对应位置处的对流层垂直柱浓度,Vi为第i次观测对应位置处的风速,sinθi第i次观测对应位置处风矢量和速度矢量之间的夹角正弦值,ΔSi第i次观测对应观测车行驶距离。
实际应用过程中,可以对于模型风场的垂直廓线分布使用气体廓线权重平均方法处理,而对于模型风场水平分布(1km×1km)可以采用空间插值的方式匹配观测位置最近处的风场结果。
为进一步验证结果的准确性,可以和对应位置处排放清单结果,以及卫星监测对比验证,在实验中,本发明实施例上述方案取得了较高准确度。
本发明实施例提供的基于移动AHSA观测的大气污染物排放通量获取方法,利用移动AHSA观测获取空间分布的大气污染物(二氧化氮、二氧化硫、甲醛等)的差分斜柱浓度结果,利用气象化学模型模拟获取的大气信息,输入到大气辐射传输模型中计算出大气光子路径信息,进一步转化差分斜柱浓度结果为对流层垂直柱浓度结果,结合气象化学模型模拟的风场信息计算获取污染源的排放通量结果。该方法使用遥感方法间接、快速、准确量化大气污染物排放通量结果,为大气环境污染提供了有力的监测手段,同时,大气污染物排放通量计算结果,为大气环境治理提供必要的数据支持。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到上述实施例可以通过软件实现,也可以借助软件加必要的通用硬件平台的方式来实现。基于这样的理解,上述实施例的技术方案可以以软件产品的形式体现出来,该软件产品可以存储在一个非易失性存储介质(可以是CD-ROM,U盘,移动硬盘等)中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述的方法。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明披露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求书的保护范围为准。

Claims (5)

1.一种基于移动AHSA观测的大气污染物排放通量获取方法,其特征在于,包括:
基于差分吸收光谱解析方法,使用光谱分析算法反演分析移动AHSA观测的散射光谱数据,获取大气污染物的差分斜柱浓度;
使用气象化学模型模拟得到大气温度、压力廓线,气溶胶以及相关大气污染物的廓线信息;
将模拟得到大气温度、压力廓线,气溶胶以及相关大气污染物的廓线信息输入到大气辐射传输模型中计算获得大气光子路径结果;
使用大气光子路径结果将大气污染物的差分斜柱浓度转化为对流层垂直柱浓度;
结合对流层垂直柱浓度与预先获得的风场信息,采用环路积分方法计算获得大气污染物的排放通量。
2.根据权利要求1所述的一种基于移动AHSA观测的大气污染物排放通量获取方法,其特征在于,所述移动AHSA观测的散射光谱数据为多个角度的散射光谱数据。
3.根据权利要求2所述的一种基于移动AHSA观测的大气污染物排放通量获取方法,其特征在于,所述使用大气光子路径结果将大气污染物的差分斜柱浓度转化为对流层垂直柱浓度的公式为:
其中,VCD为大气污染物的垂直柱浓度,DSCDtrop为大气污染物的差分斜柱浓度,AMF为大气光子路径,α为观测仰角。
4.根据权利要求1所述的一种基于移动AHSA观测的大气污染物排放通量获取方法,其特征在于,所述预先获得的风场信息为实测的风场信息,或者通过气象化学模型模拟得到的风场信息。
5.根据权利要求1所述的一种基于移动AHSA观测的大气污染物排放通量获取方法,其特征在于,所述采用环路积分方法计算获得大气污染物的排放通量的公式如下:
将环路积分转化为有限的观测片段累和,则大气污染物的排放通量具体计算过程为:
Flux=∑VCDi·Vi·sinθi·ΔSi
上式中,V代表观风矢量在行驶方向垂直方向的分量,VCDi第i次观测对应位置处的对流层垂直柱浓度,Vi为第i次观测对应位置处的风速,sinθi第i次观测对应位置处风矢量和速度矢量之间的夹角正弦值,ΔSi第i次观测对应观测车行驶距离。
CN201910338227.XA 2019-04-25 2019-04-25 基于移动ahsa观测的大气污染物排放通量获取方法 Pending CN110031412A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910338227.XA CN110031412A (zh) 2019-04-25 2019-04-25 基于移动ahsa观测的大气污染物排放通量获取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910338227.XA CN110031412A (zh) 2019-04-25 2019-04-25 基于移动ahsa观测的大气污染物排放通量获取方法

Publications (1)

Publication Number Publication Date
CN110031412A true CN110031412A (zh) 2019-07-19

Family

ID=67240254

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910338227.XA Pending CN110031412A (zh) 2019-04-25 2019-04-25 基于移动ahsa观测的大气污染物排放通量获取方法

Country Status (1)

Country Link
CN (1) CN110031412A (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111089846A (zh) * 2019-12-26 2020-05-01 中国科学院合肥物质科学研究院 一种机载与车载doas同步观测的污染源排放通量测量方法
CN111289040A (zh) * 2020-04-03 2020-06-16 中科三清科技有限公司 走航监测***
CN113281264A (zh) * 2021-05-08 2021-08-20 安徽大学 一种基于卫星遥感确定环境大气污染点源一氧化碳排放通量的方法
CN113989667A (zh) * 2021-10-26 2022-01-28 中国科学技术大学先进技术研究院 基于卫星遥感技术的近地面甲醛浓度获取方法
CN114324780A (zh) * 2022-03-03 2022-04-12 阿里巴巴达摩院(杭州)科技有限公司 大气污染物排放通量处理方法、存储介质以及计算机终端
CN114324226A (zh) * 2021-11-30 2022-04-12 中国科学技术大学 一种大气污染物立体分布无人机载超光谱遥测***
CN114357894A (zh) * 2022-03-11 2022-04-15 阿里巴巴达摩院(杭州)科技有限公司 大气污染物的处理方法、存储介质以及计算机终端
CN114414525A (zh) * 2021-11-26 2022-04-29 安徽大学 基于半导体激光no2柱浓度的探测装置及探测方法
CN115238327A (zh) * 2022-05-25 2022-10-25 南京大学 基于卫星监测的三维空间甲烷分布格局模拟方法
CN115372287A (zh) * 2022-09-12 2022-11-22 中国科学技术大学 基于高分五号卫星数据的二氧化硫排放清单构建方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108519340A (zh) * 2018-03-28 2018-09-11 刘诚 一种卫星遥感污染气体柱浓度时空分布反演算法
CN108918436A (zh) * 2018-05-08 2018-11-30 刘诚 基于max-doas对气溶胶及痕量污染气体的垂直廓线反演算法
CN108956491A (zh) * 2018-06-05 2018-12-07 中国科学院遥感与数字地球研究所 一种大气hcho遥感监测方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108519340A (zh) * 2018-03-28 2018-09-11 刘诚 一种卫星遥感污染气体柱浓度时空分布反演算法
CN108918436A (zh) * 2018-05-08 2018-11-30 刘诚 基于max-doas对气溶胶及痕量污染气体的垂直廓线反演算法
CN108956491A (zh) * 2018-06-05 2018-12-07 中国科学院遥感与数字地球研究所 一种大气hcho遥感监测方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
WEI TAN ET AL.: "Estimation of winter time NOx emissions in Hefei, a typical inland city of China, using mobile MAX-DOAS observations", 《ATMOSPHERIC ENVIRONMENT》 *
多克辛: "《大气灰霾追因与防治对策-河南省大气灰霾污染专项研究成果》", 31 March 2016, 中国环境出版社 *
胡肇焜: "面向源清单校验的污染气体立体分布及排放通量遥测技术研究", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 *
陈泽青等: "大气成分的遥感监测方法与应用", 《地球科学进展》 *

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111089846A (zh) * 2019-12-26 2020-05-01 中国科学院合肥物质科学研究院 一种机载与车载doas同步观测的污染源排放通量测量方法
CN111289040A (zh) * 2020-04-03 2020-06-16 中科三清科技有限公司 走航监测***
CN113281264A (zh) * 2021-05-08 2021-08-20 安徽大学 一种基于卫星遥感确定环境大气污染点源一氧化碳排放通量的方法
CN113989667A (zh) * 2021-10-26 2022-01-28 中国科学技术大学先进技术研究院 基于卫星遥感技术的近地面甲醛浓度获取方法
CN113989667B (zh) * 2021-10-26 2024-03-22 中国科学技术大学先进技术研究院 基于卫星遥感技术的近地面甲醛浓度获取方法
CN114414525A (zh) * 2021-11-26 2022-04-29 安徽大学 基于半导体激光no2柱浓度的探测装置及探测方法
CN114414525B (zh) * 2021-11-26 2024-05-14 安徽大学 基于半导体激光no2柱浓度的探测装置及探测方法
CN114324226A (zh) * 2021-11-30 2022-04-12 中国科学技术大学 一种大气污染物立体分布无人机载超光谱遥测***
WO2023165336A1 (zh) * 2022-03-03 2023-09-07 阿里巴巴达摩院(杭州)科技有限公司 大气污染物排放通量处理方法、存储介质以及计算机终端
CN114324780A (zh) * 2022-03-03 2022-04-12 阿里巴巴达摩院(杭州)科技有限公司 大气污染物排放通量处理方法、存储介质以及计算机终端
CN114324780B (zh) * 2022-03-03 2022-08-02 阿里巴巴达摩院(杭州)科技有限公司 大气污染物排放通量处理方法、存储介质以及计算机终端
CN114357894A (zh) * 2022-03-11 2022-04-15 阿里巴巴达摩院(杭州)科技有限公司 大气污染物的处理方法、存储介质以及计算机终端
CN115238327A (zh) * 2022-05-25 2022-10-25 南京大学 基于卫星监测的三维空间甲烷分布格局模拟方法
CN115372287A (zh) * 2022-09-12 2022-11-22 中国科学技术大学 基于高分五号卫星数据的二氧化硫排放清单构建方法
CN115372287B (zh) * 2022-09-12 2024-05-03 中国科学技术大学 基于高分五号卫星数据的二氧化硫排放清单构建方法

Similar Documents

Publication Publication Date Title
CN110031412A (zh) 基于移动ahsa观测的大气污染物排放通量获取方法
Abreu et al. Description of atmospheric conditions at the pierre auger observatory using the global data assimilation system (gdas)
US10802009B2 (en) Networked environmental monitoring system and method
Irie et al. Eight-component retrievals from ground-based MAX-DOAS observations
Eckhardt et al. Estimation of the vertical profile of sulfur dioxide injection into the atmosphere by a volcanic eruption using satellite column measurements and inverse transport modeling
Zhuang et al. Optical properties and radiative forcing of urban aerosols in Nanjing, China
Pan et al. The potential of CO2 satellite monitoring for climate governance: A review
Zhao et al. Recent progress of aerosol light-scattering enhancement factor studies in China
Zhang et al. Review of aerosol optical depth retrieval using visibility data
CN114547553B (zh) 二氧化碳排放量的反演方法、装置、设备及存储介质
Shi et al. An inversion method for estimating strong point carbon dioxide emissions using a differential absorption Lidar
CN111707622A (zh) 基于地基max-doas的大气水汽垂直分布及输送通量的测算方法
Moorthy et al. Evolution of aerosol research in India and the RAWEX–GVAX: an overview
CN114739949A (zh) 基于移动ftir遥感的点源碳排放通量测算方法
Jing et al. Cloud overlapping parameter obtained from CloudSat/CALIPSO dataset and its application in AGCM with McICA scheme
CN115480032A (zh) 基于地面遥感测量的点源排放强度预测方法
CN114077798B (zh) 估算小尺度区域大气污染物网格浓度的方法
Pereira et al. Solar irradiance modelling using an offline coupling procedure for the Weather Research and Forecasting (WRF) model
CN103674904A (zh) 红外特性测量中大气传输快速修正方法
Draxler Demonstration of a global modeling methodology to determine the relative importance of local and long-distance sources
Wang et al. Long-term evolution of clear sky surface solar radiation and its driving factors over East Asia
CN116698691B (zh) 大气细粒子aod反演方法、装置、电子设备及存储介质
Wang et al. Coupled mesoscale-LES modeling of air quality in a polluted city using WRF-LES-Chem
CN115372287B (zh) 基于高分五号卫星数据的二氧化硫排放清单构建方法
CN116822141A (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
RJ01 Rejection of invention patent application after publication

Application publication date: 20190719

RJ01 Rejection of invention patent application after publication