CN115390093A - 一种用臭氧雷达探测大气边界层的方法 - Google Patents

一种用臭氧雷达探测大气边界层的方法 Download PDF

Info

Publication number
CN115390093A
CN115390093A CN202211155179.9A CN202211155179A CN115390093A CN 115390093 A CN115390093 A CN 115390093A CN 202211155179 A CN202211155179 A CN 202211155179A CN 115390093 A CN115390093 A CN 115390093A
Authority
CN
China
Prior art keywords
ozone
signal
height
atmospheric
aerosol
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
CN202211155179.9A
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.)
Beijing Huantuo Technology Co ltd
Original Assignee
Beijing Huantuo Technology Co ltd
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 Beijing Huantuo Technology Co ltd filed Critical Beijing Huantuo Technology Co ltd
Priority to CN202211155179.9A priority Critical patent/CN115390093A/zh
Publication of CN115390093A publication Critical patent/CN115390093A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S17/00Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/88Lidar systems specially adapted for specific applications
    • G01S17/95Lidar systems specially adapted for specific applications for meteorological use
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/48Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
    • G01S7/4802Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Electromagnetism (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

一种用臭氧雷达探测大气边界层的方法,属于大气污染监测技术领域,包括以下步骤:S1臭氧雷达产生激光信号发射到大气中,受到气溶胶的反射产生回波信号,采集回波信号;S2选择回波信号中的一组,记作A;S3设定一定高度以上的大气层为洁净大气层,将其回波信号作为背景噪声,记作Am,将信号A依次减去Am,获得去噪后的信号数据;S4使用Fernald反演算法处理去除噪声后的信号数据,计算出大气气溶胶消光系数,记作E;S5使用小波协方差函数计算出大气边界层高度H。本发明计算结果满足要求,实现臭氧雷达一物多用;同时为分析大气逆温层,评估臭氧污染的扩散势能,提供了数据支撑。

Description

一种用臭氧雷达探测大气边界层的方法
技术领域
本发明属于大气污染监测技术领域,涉及一种用臭氧雷达探测大气边界层的方法。
背景技术
近年来,臭氧污染逐步成为环境污染的主要参考因素,当其浓度超过一定限值后,会刺激人的呼吸系统,破坏中枢神经,严重影响人们的正常生产与生活。大气边界层高度是逆温层的一个主要参考指标,而逆温层则直接影响了城市臭氧污染的扩散。
公布号为CN103135113A的中国发明申请,公开了一种测量大气边界层高度的方法。该方法采用基于重力波理论的算法——重力波梯度法获取边界层高度,直接将气溶胶在重力波作用下的垂直分布引入到算法理论中。
公布号为CN112146623A的中国发明申请,公开了一种大气边界层顶高度估算方法,利用探空资料统计某区域逆温层的季节分布特征;在近地面逆温占主导地位的季节,统计近地面逆温顶部所在高度的数分布特征,取其上四分之一位数作为该季节的数阈值;在对流层低层逆温占主导地位的季节,统计对流层低层逆温底部所在高度的数分布特征,取其下四分之一位数作为该季节的数阈值。
公布号为CN114282170A的中国发明申请,公开了一种准确计算大气边界层高度的方法,首先通过三次样条差值法,把探空资料的原始数据转换成垂直分辨率为10m的m行×12列的数据,从转化后的数据文件中读取温度值、海拔高度值和气压值;然后计算绝对温度和位温值;再计算两层位温差,计算时设定从距地面600m以上开始,最早出现该层位温与100m以下层的位温差超过0.3℃,同时设定该层以上100m处位温与该层位温差超过0.3℃,表明该层所在高度为大气稳定出现逆温层的高度,则将该层底部以下范围确定为边界层高度。
公布号为CN110058258A的中国发明申请,公开了一种基于混合型激光雷达的大气边界层探测方法。其基本原理为:利用一台混合型激光雷达,集成直接探测气溶胶激光雷达系统和相干测风激光雷达系统,同时探测气溶胶光学特性和湍流强度的垂直廓线,分别反演出大气边界层高度,利用两者差值与选定阈值的比较,判定出不同类型的大气边界层。
上述现有技术在探测大气边界层时,都需要专门的设备。如果在臭氧雷达探测臭氧时空分布的同时,利用臭氧雷达探测出边界层高度,一方面可以使臭氧雷达一物多用、降低城市污染监测治理成本;另一方面可以将臭氧污染分布与边界层高度快速结合,为分析城市臭氧污染的扩散特性,预测未来时间内臭氧的变化趋势提供数据参考。
发明内容
本发明的目的是提供一种用臭氧雷达探测大气边界层的方法,在臭氧雷达探测臭氧时空分布的同时,利用臭氧雷达探测出边界层高度,实现臭氧雷达一物多用,降低城市污染监测治理成本;同时将臭氧污染分布与边界层高度快速结合,为分析城市臭氧污染的扩散特性,预测未来时间内臭氧的变化趋势提供数据参考。本发明的目的通过以下技术方案实现。
一种用臭氧雷达探测大气边界层的方法,其特征在于,包括以下步骤:
S1臭氧雷达产生三组不同波长的激光信号发射到大气中,这三组激光信号在大气中受到气溶胶的反射产生回波信号,计算机采集回波信号。
S2选择计算机所采集的回波信号中的一组,记作A。所选择的回波信号具备高度特性,若高度分辨率为30米,监测总高度为30千米,则共计有1000个信号梯度点。
S3根据经验一定高度以上为洁净的大气层,这个范围内的回波信号基本是由大气分子所产生,故将其作为背景噪声,记作Am,将信号A依次减去Am去除背景噪声,获得去噪后的信号数据。
S4使用Fernald反演算法对去除噪声后的信号数据进行处理,计算出大气气溶胶消光系数,记作E。
S5根据气溶胶消光系数E在大气边界层位置附近有突变的特性,使用小波协方差函数计算出大气边界层高度H。
上述方法中,步骤S1同时是臭氧雷达探测臭氧时空分布的步骤,臭氧雷达探测臭氧时空分布的过程为:通过高能紫外激光器发射多组波长接近的脉冲激光,激光经过扩束器发射到大气中与臭氧、气溶胶等发生相互作用,后向散射光被望远镜接收,得到各波长的回波信号,由差分吸收激光雷达算法反演出大气中臭氧的浓度。
进一步地,步骤S2选择三组激光信号中波长最长的那一组的回波信号。
进一步地,步骤S3所述的一定高度为25千米。
进一步地,步骤S4所述Fernald反演算法的计算公式如下:
Figure BDA0003858063060000021
式中:
αa(R)为R处气溶胶消光系数,无因次;
αm(R)为R处大气分子消光系数,无因次,根据标准大气分子模型确定;
R为探测高度,单位千米;
Rc为不含气溶胶的洁净大气层所在高度,单位千米,可设定为25千米;
P(R)为R处的回波信号,无因次;
r为高度,单位千米;
Sm为大气分子后向散射比,无因次,取值8π/3;
Sa为气溶胶消光后向散射比,无因次。
进一步地,步骤S5所述小波协方差函数的计算公式如下:
Figure BDA0003858063060000031
式中:
zt为信号h距离上限,单位米;
zb为信号h距离下限,单位米;
a为相邻信号距离间隔,单位米;
z为高度,单位米;
E(z)为Z处信号对应的气溶胶消光系数;
Wf(a,b)取得最小值时,对应的位置b即为边界层高度。
本发明具有以下有益技术效果:
在臭氧雷达监测臭氧时空分布的同时,可计算出大气边界层高度,计算结果与专门的气溶胶雷达监测结果相比,满足误差要求,实现臭氧雷达一物多用、降低城市污染监测治理成本;同时为分析大气逆温层,评估臭氧污染的扩散势能,提供了有效的数据支撑。
附图说明
图1是实施例1在时间点T1、T2分别采用臭氧雷达和气溶胶雷达获得的反演消光数据。
图2是实施例1在时间点T3、T4分别采用臭氧雷达和气溶胶雷达获得的反演消光数据。
图3是实施例1分别采用臭氧雷达和气溶胶雷达获得的边界层高度。
图4是实施例2利用臭氧雷达计算的边界层高度曲线。
具体实施方式
下面结合说明书附图,对本发明的技术方案进行清楚、完整的描述。显然,所描述的实施方式仅仅是本发明一部分实施方式,而不是全部的实施方式。基于本发明的实施方式,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施方式,都属于本发明的保护范围。
实施例1
臭氧雷达激光器产生266nm激光,通过拉曼管激发出另外的289nm和316nm两种不同波段的激光,三个波段激光通过扩束器放大光束直径后进入发射接收装置。三个波段的激光经过反射镜反射后进入大气环境中,与大气中的气溶胶颗粒等反射产生回波信号,回波信号进入发射接收装置后汇集并输出至光谱仪。光谱仪接收到回波信号后,经过内部处理分离出316nm回波信号,并输出至采集卡进行光子计数。
选取某地区的臭氧雷达四天的监测数据,为了确定本方法计算的准确性,选择臭氧雷达附近的一台具备认证资质的气溶胶激光雷达同时段的监测数据作为标准比较。
因为大气气溶胶消光是边界层高度的计算输入,为了验证本方法的正确性,首先需要确定大气气溶胶消光系数计算的准确性。
从大气污染较重的一天中选择两个时间点T1、T2,使用臭氧雷达T1、T2时间点监测数据计算气溶胶消光系数,同时选择对应时间气溶胶激光雷达计算消光系数,绘制廓线图。使用3千米内的数据做比对,如图1所示。从图1可以看出,臭氧雷达和气溶胶雷达计算得到的消光系数非常相似,通过相似性计算,T1时间点两者相似度为93.7%,T2时间点两者相似度为90.5%,满足国标中关于激光雷达测量消光系数偏差不大于20%的标准。
从大气污染较轻的一天中选择两个时间点T3、T4,使用臭氧雷达T3、T4时间点监测数据计算气溶胶消光系数,同时选择对应时间气溶胶激光雷达计算消光系数,绘制廓线图。使用8千米内的数据做比对,如图2所示。从图2可以看出,臭氧雷达和气溶胶雷达计算得到的消光系数非常相似,通过相似性计算,T3时间点两者相似度为93.7%,T4时间点两者相似度为91.5%,满足国标中关于激光雷达测量消光系数偏差不大于20%的标准。
上述两组比对说明,臭氧雷达反演的消光系数,在不同天气条件下与气溶胶雷达反演的结果相似性很高。由于气溶胶消光系数在大气边界层位置附近有突变的特性,故而二者得到的边界层高度也非常一致,结果如图3所示。通过相似性计算,两者相似度为91.3%。表明通过臭氧雷达计算出的大气边界层高度,结果可靠。
实施例2
通过本发明提供的方法,使用某一段时间臭氧雷达监测数据,计算边界层高度数据并叠加至臭氧时空分布图谱,如图4所示。其中黑色线条为大气边界层高度。
根据图4直观获知当边界层高度降低时,臭氧在近地面聚集,污染加重,反之当边界层高度升高时,臭氧在大气中扩散能力增强,近地面臭氧污染降低。这有利于研究臭氧污染变化的特征,首先可以根据边界层高度变化趋势,分析预测未来一段时间内臭氧污染变化情况;其次当监测到边界层高度降低时,可以及时采取手段,管控参与臭氧光化学反应前体物氮氧化物或挥发性有机物的排放,减少臭氧的产生,从而避免臭氧在近地面的堆积,形成臭氧污染爆表。
尽管上面已经示出和描述了本发明的实施方式,可以理解的是,上述实施方式是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施方式进行变化、修改、替换和变型。本发明的保护范围由权利要求书及其等同技术方案限定。

Claims (5)

1.一种用臭氧雷达探测大气边界层的方法,其特征在于,包括以下步骤:
S1臭氧雷达产生三组不同波长的激光信号发射到大气中,这三组激光信号在大气中受到气溶胶的反射产生回波信号,计算机采集回波信号;
S2选择计算机所采集的回波信号中的一组,记作A;
S3设定一定高度以上的大气层为洁净大气层,将其回波信号作为背景噪声,记作Am,将信号A依次减去Am去除背景噪声,获得去噪后的信号数据;
S4使用Fernald反演算法对去除噪声后的信号数据进行处理,计算出大气气溶胶消光系数,记作E;
S5根据气溶胶消光系数E在大气边界层位置附近有突变的特性,使用小波协方差函数计算出大气边界层高度H。
2.根据权利要求1所述的方法,其特征在于,步骤S2选择三组激光信号中波长最长的那一组的回波信号。
3.根据权利要求1所述的方法,其特征在于,步骤S3所述的一定高度为25千米。
4.根据权利要求1所述的方法,其特征在于,步骤S4所述Fernald反演算法的计算公式如下:
Figure FDA0003858063050000011
式中:
αa(R)为R处气溶胶消光系数,无因次;
αm(R)为R处大气分子消光系数,无因次,根据标准大气分子模型确定;
R为探测高度,单位千米;
Rc为不含气溶胶的洁净大气层所在高度,单位千米,设定为25千米;
P(R)为R处的回波信号,无因次;
r为高度,单位千米;
Sm为大气分子后向散射比,无因次,取值8π/3;
Sa为气溶胶消光后向散射比,无因次。
5.根据权利要求1所述的方法,其特征在于,步骤S5所述小波协方差函数的计算公式如下:
Figure FDA0003858063050000021
式中:
zt为信号h距离上限,单位米;
zb为信号h距离下限,单位米;
a为相邻信号距离间隔,单位米;
z为高度,单位米;
E(z)为Z处信号对应的气溶胶消光系数;
Wf(a,b)取得最小值时,对应的位置b即为边界层高度。
CN202211155179.9A 2022-09-22 2022-09-22 一种用臭氧雷达探测大气边界层的方法 Pending CN115390093A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211155179.9A CN115390093A (zh) 2022-09-22 2022-09-22 一种用臭氧雷达探测大气边界层的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211155179.9A CN115390093A (zh) 2022-09-22 2022-09-22 一种用臭氧雷达探测大气边界层的方法

Publications (1)

Publication Number Publication Date
CN115390093A true CN115390093A (zh) 2022-11-25

Family

ID=84126158

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211155179.9A Pending CN115390093A (zh) 2022-09-22 2022-09-22 一种用臭氧雷达探测大气边界层的方法

Country Status (1)

Country Link
CN (1) CN115390093A (zh)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105158230A (zh) * 2015-09-08 2015-12-16 中国科学院合肥物质科学研究院 基于ccd成像激光雷达测量大气边界层污染气体的装置
KR20170116805A (ko) * 2016-04-12 2017-10-20 부경대학교 산학협력단 다중 각도에서 태양광 관측기반의 대기 분자의 라만 산란특성을 이용한 에어로솔 고도 산출방법
CN110031868A (zh) * 2019-05-15 2019-07-19 国耀量子雷达科技有限公司 一种基于相干测风激光雷达载噪比反演边界层高度的方法
CN110058258A (zh) * 2019-05-15 2019-07-26 国耀量子雷达科技有限公司 一种基于混合型激光雷达的大气边界层探测方法
WO2020025984A1 (en) * 2018-08-01 2020-02-06 Pantazis Alexandros Method of use of a lidar device and operatively associated lidar data processing unit for providing real-time monitoring of meteorological parameters
WO2020071562A1 (ja) * 2018-10-05 2020-04-09 英弘精機株式会社 気象観測用ライダー
CN111427064A (zh) * 2020-04-27 2020-07-17 安徽皖仪科技股份有限公司 一种基于微脉冲激光雷达的气溶胶消光系数反演方法
CN113687380A (zh) * 2021-07-07 2021-11-23 中科经纬检测科技(苏州)有限公司 一种基于水平扫描激光雷达的臭氧浓度反演方法
CN215813353U (zh) * 2021-06-03 2022-02-11 北京环拓科技有限公司 一种光化学与臭氧激光雷达协同分析的臭氧监测系统

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105158230A (zh) * 2015-09-08 2015-12-16 中国科学院合肥物质科学研究院 基于ccd成像激光雷达测量大气边界层污染气体的装置
KR20170116805A (ko) * 2016-04-12 2017-10-20 부경대학교 산학협력단 다중 각도에서 태양광 관측기반의 대기 분자의 라만 산란특성을 이용한 에어로솔 고도 산출방법
WO2020025984A1 (en) * 2018-08-01 2020-02-06 Pantazis Alexandros Method of use of a lidar device and operatively associated lidar data processing unit for providing real-time monitoring of meteorological parameters
WO2020071562A1 (ja) * 2018-10-05 2020-04-09 英弘精機株式会社 気象観測用ライダー
US20210333363A1 (en) * 2018-10-05 2021-10-28 Eko Instruments Co., Ltd. Meteorological lidar
CN110031868A (zh) * 2019-05-15 2019-07-19 国耀量子雷达科技有限公司 一种基于相干测风激光雷达载噪比反演边界层高度的方法
CN110058258A (zh) * 2019-05-15 2019-07-26 国耀量子雷达科技有限公司 一种基于混合型激光雷达的大气边界层探测方法
CN111427064A (zh) * 2020-04-27 2020-07-17 安徽皖仪科技股份有限公司 一种基于微脉冲激光雷达的气溶胶消光系数反演方法
CN215813353U (zh) * 2021-06-03 2022-02-11 北京环拓科技有限公司 一种光化学与臭氧激光雷达协同分析的臭氧监测系统
CN113687380A (zh) * 2021-07-07 2021-11-23 中科经纬检测科技(苏州)有限公司 一种基于水平扫描激光雷达的臭氧浓度反演方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
应传友: "闽江口臭氧传输污染观测研究", 海峡科学, no. 6, 30 June 2018 (2018-06-30), pages 37 - 39 *

Similar Documents

Publication Publication Date Title
CN102944503B (zh) 基于太阳光度计和激光雷达的pm2.5质量浓度值自动反演方法
CN110441777A (zh) 一种基于激光雷达的气溶胶垂直廓线的反演方法
CN103278479B (zh) 大气辐射传输修正系统及修正方法
CN103175759A (zh) 基于多种地基遥感技术获取城市气溶胶复折射指数的方法
CN101520515B (zh) 基于激光雷达的混合层高度的自动反演方法
Hong et al. Ship-based MAX-DOAS measurements of tropospheric NO 2, SO 2, and HCHO distribution along the Yangtze River
KR101012279B1 (ko) 라이다를 이용한 황사 판독방법
CN103293117A (zh) 一种微脉冲差分吸收激光雷达水汽时空分布反演方法
CN107015243A (zh) 一种基于布里渊激光雷达系统的大气温度测量方法
CN112684471B (zh) 一种基于激光雷达的气溶胶微物理特性反演方法
Shi et al. Measurement of CO2 rectifier effect during summer and winter using ground-based differential absorption LiDAR
CN103728612A (zh) 基于目标红外辐射光谱和带模式的被动测距方法
CN104076345B (zh) 一种针对纯转动拉曼激光雷达测量温度的饱和修正方法
Poggio et al. Scintillometer wind measurements over complex terrain
CN109856021A (zh) 一种pm2.5反演方法及监测区域分割方法
CN113295626A (zh) 基于阵列式微脉冲激光雷达的气溶胶消光光谱测量方法和装置
CN114460602A (zh) 近地层无盲区大气温度探测拉曼激光雷达系统及方法
Yufeng et al. A UV multifunctional Raman lidar system for the observation and analysis of atmospheric temperature, humidity, aerosols and their conveying characteristics over Xi'an
Zhu et al. Retrieving the vertical distribution of PM 2.5 mass concentration from lidar via a random forest model
Nan et al. Observations of nitrous acid and its relative humidity dependence in Shanghai
CN116956749B (zh) 一种基于大数据多维度监测碳排放的方法和系统
CN115524264B (zh) 一种基于激光雷达测量的气溶胶分类方法
CN115390093A (zh) 一种用臭氧雷达探测大气边界层的方法
CN113156452B (zh) 一种星载差分吸收激光雷达co2廓线探测最优波段确定方法
Devara et al. Remote sounding of aerosols in the lower atmosphere using a bistatic CW helium-neon lidar

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