CN106772422A - 高空大气密度在线探测及计算方法 - Google Patents
高空大气密度在线探测及计算方法 Download PDFInfo
- Publication number
- CN106772422A CN106772422A CN201610985247.2A CN201610985247A CN106772422A CN 106772422 A CN106772422 A CN 106772422A CN 201610985247 A CN201610985247 A CN 201610985247A CN 106772422 A CN106772422 A CN 106772422A
- Authority
- CN
- China
- Prior art keywords
- laser
- density
- atmospheric
- laser beam
- video camera
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/88—Lidar systems specially adapted for specific applications
- G01S17/89—Lidar systems specially adapted for specific applications for mapping or imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/47—Scattering, i.e. diffuse reflection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/47—Scattering, i.e. diffuse reflection
- G01N2021/4704—Angular selective
- G01N2021/4709—Backscatter
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Immunology (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Pathology (AREA)
- Health & Medical Sciences (AREA)
- Computer Networks & Wireless Communication (AREA)
- Electromagnetism (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Optical Radar Systems And Details Thereof (AREA)
Abstract
本发明公开了一种高空大气密度在线探测及计算方法,以连续激光为发射光源,利用摄像机接收来自大气的激光回波信号,通过对激光回波图像进行分析处理,得到大气密度信息。本发明采用的设备轻巧、具有高时空分辨率,并具有快速探测、可昼夜连续运行、自动化程度高、无盲区和易于移动等优点可用于大多数军用、民用飞行器对大气密度的测量。
Description
技术领域
本发明属于一种激光探测技术领域,具体是一种高空大气密度在线探测及计算方法。
背景技术
高空大气密度是影响高空飞行器速度和方向的重要因素。地基所采用的现场测量方法难以得到高层大气特别是高层大气实时背景条件下的大气密度参数。本方法使用激光遥感手段,安装在高空飞行器上,实时得到运动路径周围的大气密度参量,为飞行器提供第一手的数据。
发明内容
本发明提供一种高空大气密度在线探测及计算方法,采用激光发射和接收分离设计,使用二极管泵浦连续激光器发射激光、摄像机接收信号;通过几何结构设计,利用中焦距镜头将测量范围全部纳入摄像机视场,通过计算机进行计算分析,得到大气密度分布。
本发明采用的技术方案是:
高空大气密度在线探测及计算方法,其特征在于,以连续激光为发射光源,利用摄像机接收来自大气的激光回波信号,通过对激光回波图像进行分析处理,得到大气密度信息;具体包括以下步骤:
(1)采用二极管泵浦连续激光器发射激光,激光器的前方光路上设置分光镜,分光镜的反射光路上设置功率监测仪,激光器、功率监测仪与计算机连接;
(2)激光器输出的激光束,一部分激光束透过分光镜射向大气中,一部分激光束镜反光镜反射后进入功率监测仪;
(3)射向大气中的激光束,受到大气中大气分子的作用,产生瑞利散射,瑞利散射的一部分散射光,被光学接收部分收集,光学接收部分包括有滤光片、CCD摄像机,CCD摄像机通过A/D转换器与计算机连接;
(4)通过CCD摄像机采集激光回波图像,每幅图像曝光时间一定,激光回波图像由几列光点和背景组成,几列长条光点由密到疏对应高度分布,计算机对采集得到的激光回波图像进行图像处理、平均和去噪,再进行计算分析,反演得出大气密度;
(5)数据反演方程可表示为:
Er=KElATatmzTatmRn(z)σRaydθ/D (1)
其中,Er为发射距离z处视场角dθ内单位象素接收能量,在光斑图像区域内假设dθ不变;
σRay为所用波长的瑞利分子后项散射截面;
D为CCD与激光器的距离;
EL曝光时间内激光脉冲总能量;
A为接收有效孔径面积;
dθ为光斑图像单位角度对应距离分辨率;
K是设备定标常数,取决于系统光学效率;
n(z)为距离z处大气分子数密度;
Tatmz为从激光器到z处的大气透过率;
TatmR为测量点到CCD摄像机间的大气透过率;对于532nm波长从30公里到80公里大气透过率变化小于0.5%(根据中分辨率大气辐射传输软件Modtran),这里测量距离也较短,大气透过率可以当作常数;
大气密度与大气分子数密度n有如下关系:
ρ(z)=n(z)k/R (2)
ρ(z)为z高度上的大气密度;k为波尔兹曼常数,R气体常数;
(6)最后,根据公式(1)和(2)计算得到大气密度。
所述的高空大气密度在线探测及计算方法,其特征在于,所述激光器采用激光二极管泵浦连续激光器,波长532nm,能量2W。
所述的高空大气密度在线探测及计算方法,其特征在于,所述CCD摄像机的像素单元数为1004x 1002,使用中焦距镜头,镜头口径为60mm、焦距为50mm。
本发明的大气密度探测原理:
主要通过测量大气中气溶胶和分子等的后向散射信号,反演得到测量区域的大气参数。高度约从30公里至100公里的高层大气中,气溶胶含量极微,除特殊情况外,一般可把这一高度范围的大气当作无气溶胶的清洁大气。所以激光雷达接收的大气回波主要是大气分子的的Rayleigh散射信号。在设备测量范围内高层大气的衰减很小,大气消光影响很小,大气散射信号的强度仅与大气分子密度有关,可以由散射信号强度反演出大气密度廓线。
本发明的优点是:
本发明采用的设备轻巧、具有高时空分辨率,并具有快速探测、可昼夜连续运行、自动化程度高、无盲区和易于移动等优点可用于大多数军用、民用飞行器对大气密度的测量。
附图说明
图1是激光发射和光学接收部分的结构框图。
图2是激光雷达结构原理图。
具体实施方式
如图1、2所示,高空大气密度在线探测及计算方法,以连续激光为发射光源,利用摄像机接收来自大气的激光回波信号,通过对激光回波图像进行分析处理,得到大气密度信息;具体包括以下步骤:
(1)采用二极管泵浦连续激光器发射激光,激光器的前方光路上设置分光镜,分光镜的反射光路上设置功率监测仪,激光器、功率监测仪与计算机连接;
(2)激光器输出的激光束,一部分激光束透过分光镜射向大气中,一部分激光束镜反光镜反射后进入功率监测仪;
(3)射向大气中的激光束,受到大气中大气分子的作用,产生瑞利散射,瑞利散射的一部分散射光,被光学接收部分收集,光学接收部分包括有滤光片、CCD摄像机,CCD摄像机通过A/D转换器与计算机连接;
(4)通过CCD摄像机采集激光回波图像,每幅图像曝光时间一定,激光回波图像由几列光点和背景组成,几列长条光点由密到疏对应高度分布,计算机对采集得到的激光回波图像进行图像处理、平均和去噪,再进行计算分析,反演得出大气密度;
(5)数据反演方程可表示为:
Er=KElATatmzTatmRn(z)σRaydθ/D (1)
其中,Er为发射距离z处视场角dθ内单位象素接收能量,在光斑图像区域内假设dθ不变;
σRay为所用波长的瑞利分子后项散射截面;
D为CCD与激光器的距离;
EL曝光时间内激光脉冲总能量;
A为接收有效孔径面积;
dθ为光斑图像单位角度对应距离分辨率;
K是设备定标常数,取决于系统光学效率;
n(z)为距离z处大气分子数密度;
Tatmz为从激光器到z处的大气透过率;
TatmR为测量点到CCD摄像机间的大气透过率;对于532nm波长从30公里到80公里大气透过率变化小于0.5%(根据中分辨率大气辐射传输软件Modtran),这里测量距离也较短,大气透过率可以当作常数;
大气密度与大气分子数密度n有如下关系:
ρ(z)=n(z)k/R (2)
ρ(z)为z高度上的大气密度;k为波尔兹曼常数,R气体常数;
(6)最后,根据公式(1)和(2)计算得到大气密度。
激光器采用激光二极管泵浦连续激光器,波长532nm,能量2W。
激光发射部分的技术参数如表1所示:
表1激光发射部分技术参数
名称 | 参数 |
泵浦方式 | 激光二极管泵浦 |
激光波长 | 532nm |
能量 | 2W |
M2factor | <2.0 |
Power stability | ≤1% |
出光孔径 | ≤3mm |
CCD摄像机的像素单元数为1004x 1002,使用中焦距镜头,镜头口径为60mm、焦距为50mm。为了消除杂散光干扰,接收系统安装窄带滤光片。根据测量的远近设计发射和接收夹角,1msec一组数据。信号接收部分技术参数如表2所示。
表2信号接收部分技术参数
Claims (3)
1.高空大气密度在线探测及计算方法,其特征在于,以连续激光为发射光源,利用摄像机接收来自大气的激光回波信号,通过对激光回波图像进行分析处理,得到大气密度信息;具体包括以下步骤:
(1)采用二极管泵浦连续激光器发射激光,激光器的前方光路上设置分光镜,分光镜的反射光路上设置功率监测仪,激光器、功率监测仪与计算机连接;
(2)激光器输出的激光束,一部分激光束透过分光镜射向大气中,一部分激光束镜反光镜反射后进入功率监测仪;
(3)射向大气中的激光束,受到大气中大气分子的作用,产生瑞利散射,瑞利散射的一部分散射光,被光学接收部分收集,光学接收部分包括有滤光片、CCD摄像机,CCD摄像机通过A/D转换器与计算机连接;
(4)通过CCD摄像机采集激光回波图像,每幅图像曝光时间一定,激光回波图像由几列光点和背景组成,几列长条光点由密到疏对应高度分布,计算机对采集得到的激光回波图像进行图像处理、平均和去噪,再进行计算分析,反演得出大气密度;
(5)数据反演方程可表示为:
Er=KElATatmzTatmRn(z)σRaydθ/D (1)
其中,Er为发射距离z处视场角dθ内单位象素接收能量,在光斑图像区域内假设dθ不变;
σRay为所用波长的瑞利分子后项散射截面;
D为CCD与激光器的距离;
EL曝光时间内激光脉冲总能量;
A为接收有效孔径面积;
dθ为光斑图像单位角度对应距离分辨率;
K是设备定标常数,取决于系统光学效率;
n(z)为距离z处大气分子数密度;
Tatmz为从激光器到z处的大气透过率;
TatmR为测量点到CCD摄像机间的大气透过率;
大气密度与大气分子数密度n有如下关系:
ρ(z)=n(z)k/R (2)
ρ(z)为z高度上的大气密度;k为波尔兹曼常数,R气体常数;
(6)最后,根据公式(1)和(2)计算得到大气密度。
2.根据权利要求1所述的高空大气密度在线探测及计算方法,其特征在于,所述激光器采用激光二极管泵浦连续激光器,波长532nm,能量2W。
3.根据权利要求1或2所述的高空大气密度在线探测及计算方法,其特征在于,所述CCD摄像机的像素单元数为1004x 1002,使用中焦距镜头,镜头口径为60mm、焦距为50mm。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610985247.2A CN106772422A (zh) | 2016-10-25 | 2016-10-25 | 高空大气密度在线探测及计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610985247.2A CN106772422A (zh) | 2016-10-25 | 2016-10-25 | 高空大气密度在线探测及计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106772422A true CN106772422A (zh) | 2017-05-31 |
Family
ID=58972854
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610985247.2A Pending CN106772422A (zh) | 2016-10-25 | 2016-10-25 | 高空大气密度在线探测及计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106772422A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109253797A (zh) * | 2018-10-15 | 2019-01-22 | 湖北航天技术研究院总体设计所 | 一种激光辐照效应试验系统及方法 |
CN110849769A (zh) * | 2019-10-28 | 2020-02-28 | 北京空间机电研究所 | 一种基于可调谐激光的掩星大气密度廓线测量系统及方法 |
CN111398993A (zh) * | 2020-05-20 | 2020-07-10 | 青岛镭测创芯科技有限公司 | 一种无盲区脉冲相干测风激光雷达系统 |
CN113624640A (zh) * | 2021-06-30 | 2021-11-09 | 北京空间机电研究所 | 用于探测大气温度和密度廓线的临边散射探测装置及方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020118352A1 (en) * | 2001-02-23 | 2002-08-29 | Japan Atomic Energy Research Institute | Fast gate scanning three-dimensional laser radar apparatus |
CN101952709A (zh) * | 2007-11-15 | 2011-01-19 | 爱克斯崔里斯科技有限公司 | 颗粒探测 |
CN102891063A (zh) * | 2012-09-20 | 2013-01-23 | 中国科学院空间科学与应用研究中心 | 空间用四极滤质器及同时获得大气密度和大气成分的方法 |
CN103267969A (zh) * | 2013-04-25 | 2013-08-28 | 中国科学院安徽光学精密机械研究所 | 基于激光光柱的成像激光雷达测量大气光学湍流廓线的方法 |
-
2016
- 2016-10-25 CN CN201610985247.2A patent/CN106772422A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020118352A1 (en) * | 2001-02-23 | 2002-08-29 | Japan Atomic Energy Research Institute | Fast gate scanning three-dimensional laser radar apparatus |
CN101952709A (zh) * | 2007-11-15 | 2011-01-19 | 爱克斯崔里斯科技有限公司 | 颗粒探测 |
CN102891063A (zh) * | 2012-09-20 | 2013-01-23 | 中国科学院空间科学与应用研究中心 | 空间用四极滤质器及同时获得大气密度和大气成分的方法 |
CN103267969A (zh) * | 2013-04-25 | 2013-08-28 | 中国科学院安徽光学精密机械研究所 | 基于激光光柱的成像激光雷达测量大气光学湍流廓线的方法 |
Non-Patent Citations (1)
Title |
---|
王超杰等: "Research on Atmospheric Density Online Measurement Technology Based on Spacecraft Platform", 《2015 INTERNATIONAL CONFERENCE ON MODELLING,SIMULATION AND APPLIED MATHEMATICS》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109253797A (zh) * | 2018-10-15 | 2019-01-22 | 湖北航天技术研究院总体设计所 | 一种激光辐照效应试验系统及方法 |
CN110849769A (zh) * | 2019-10-28 | 2020-02-28 | 北京空间机电研究所 | 一种基于可调谐激光的掩星大气密度廓线测量系统及方法 |
CN110849769B (zh) * | 2019-10-28 | 2022-07-29 | 北京空间机电研究所 | 一种基于可调谐激光的掩星大气密度廓线测量系统及方法 |
CN111398993A (zh) * | 2020-05-20 | 2020-07-10 | 青岛镭测创芯科技有限公司 | 一种无盲区脉冲相干测风激光雷达系统 |
CN111398993B (zh) * | 2020-05-20 | 2023-03-14 | 青岛镭测创芯科技有限公司 | 一种无盲区脉冲相干测风激光雷达系统 |
CN113624640A (zh) * | 2021-06-30 | 2021-11-09 | 北京空间机电研究所 | 用于探测大气温度和密度廓线的临边散射探测装置及方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11480509B2 (en) | System and method for individual particle sizing using light scattering techniques | |
US9927356B2 (en) | Systems and methods for detecting gases, airborne compounds, and other particulates | |
Mehendale et al. | Review on lidar technology | |
US6836285B1 (en) | Lidar with streak-tube imaging,including hazard detection in marine applications; related optics | |
CN106772422A (zh) | 高空大气密度在线探测及计算方法 | |
CN102621102B (zh) | 基于ccd激光雷达的水平能见度的测量方法 | |
Henriksson et al. | Continuously scanning time-correlated single-photon-counting single-pixel 3-D lidar | |
CN110045391B (zh) | 一种用于气溶胶尺度谱测量的高光谱激光雷达系统 | |
CN105911562B (zh) | 一种用于环境监测的三维扫描式多参量激光雷达系统 | |
CN110095423B (zh) | 一种多通道污染排放应急遥测成像系统及方法 | |
WO2019152787A1 (en) | Apparatuses and methods for gas flux measurements | |
CN110264570A (zh) | 一种基于无人机的林地自主巡检系统 | |
CN102759753A (zh) | 隐藏危险品检测方法及设备 | |
CN106646429B (zh) | 一种用于激光雷达的自标定几何因子的装置及方法 | |
CN104316443A (zh) | 一种基于ccd后向散射的pm 2.5浓度监测方法 | |
CN108717718A (zh) | 基于层析成像的测量瞬态燃烧场三维结构空间分布的方法 | |
CN106597464A (zh) | 基于沙姆定律的扫描光谱激光雷达 | |
CN110161280A (zh) | 混合探测多普勒激光雷达风速测量系统及其测量方法 | |
Kapustin et al. | Active pulse television measuring systems for ensuring navigation of transport means in heavy weather conditions | |
Kong et al. | Three-wavelength polarization Scheimpflug lidar system developed for remote sensing of atmospheric aerosols | |
CN108562887A (zh) | 一种基于波分复用的多波长激光雷达 | |
Steinvall et al. | Active and passive short-wave infrared and near-infrared imaging for horizontal and slant paths close to ground | |
CN102435562A (zh) | 机载大气痕量气体二维分布快速监测系统及方法 | |
CN102788678A (zh) | 一种海洋湍流参数的非接触测量方法及系统 | |
CN202693258U (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170531 |