CN111458727B - 基于相干激光雷达谱强度数据的大气能见度分析方法 - Google Patents
基于相干激光雷达谱强度数据的大气能见度分析方法 Download PDFInfo
- Publication number
- CN111458727B CN111458727B CN202010309070.0A CN202010309070A CN111458727B CN 111458727 B CN111458727 B CN 111458727B CN 202010309070 A CN202010309070 A CN 202010309070A CN 111458727 B CN111458727 B CN 111458727B
- Authority
- CN
- China
- Prior art keywords
- atmospheric
- extinction coefficient
- laser radar
- alpha
- intensity data
- 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.)
- Active
Links
Images
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/95—Lidar systems specially adapted for specific applications for meteorological use
-
- 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/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/4802—Details 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
-
- 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
- G01N2021/1738—Optionally different kinds of measurements; Method being valid for different kinds of measurement
- G01N2021/1744—Optionally different kinds of measurements; Method being valid for different kinds of measurement either absorption or scatter
-
- 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
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Health & Medical Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Electromagnetism (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Optical Radar Systems And Details Thereof (AREA)
Abstract
基于相干激光雷达谱强度数据的大气能见度分析方法,根据相干激光雷达的信号探测流程分析,其不同距离门谱数据最大强度包含大气回波信息。提取出对应高度的谱强度数据,可对应包含大气消光系数与后向散射系数大气激光雷达方程。利用大气后向散射系数同大气消光系数之间的计算关系,反演出大气消光系数。结合大气能见度同消光系数之间的关系,最终得出不同高度处大气能见度结果。相干激光雷达进行能见度探测,扩展了相干雷达的应用功能,提高了相干系统的数据利用率。相比传统的气溶胶激光雷达探测,谱强度数据信噪比较大,探测距离及精度有明显优势。并且相干激光雷达结构紧凑,光学部件稳定可靠,使用的环境可靠性相比传统气溶胶雷达更有优势。
Description
技术领域
本发明涉及一种激光雷达探测方法,具体涉及一种基于相干激光雷达谱强度数据测量大气能见度信息。
背景技术
大气信息,如风速,温度,湿度,气溶胶,云等信息具有十分重要的作用。无论对天气预报,机场安全高效运行,大气污染物的监测,还是在军事方面的用途,都有十分重要的作用。
激光雷达的测量信息时空分辨率高,能够连续测量,相比传统的大气探测方式而言优势明显。相干激光雷达系统目前已有成熟的应用,但主要用来进行大气风场测量。相干激光雷达中,将回波强度通过傅里叶变换转换为谱数据。谱数据信息包含两个方面,频率变化信息及谱强度信息。
由于大气回波的多普勒效应,谱数据的频率变化信息可用于反演大气风速;谱强度信息中蕴含大气消光信息,可基于谱强度信息反演大气能见度。
常规的相干激光雷达的测量中,仅使用谱数据的频率波动信息,反演大气;而谱数据的强度则不作考虑,限制了相干雷达的应用范围。而常规的大气测量中,大气风场及大气能见度作为两个大气的基本参数,存在共同测量的需要。因此,通常在配备一台相干激光雷达测量大气风场时,还需同时配备气溶胶雷达或能见度仪,来同步获取大气能见度信息。
目前尚未出现利用相干激光雷达进行大气能见度分析的应用。
发明内容
本发明的目的是提供一种基于相干激光雷达谱强度数据的大气能见度分析方法,该方法利用相干激光雷达探测数据,将不同高度处谱强度信息对应激光雷达方程,进行谱数据强度分析,从而反演大气能见度信息。
一种基于相干激光雷达谱强度数据的大气能见度分析方法,其特征在于包括以下步骤:
步骤1,对相干激光雷达的光电探测器采样信号进行变换,获得对应于不同高度的谱数据;
步骤2,对不同高度的谱数据进行强度检测,得出不同高度处最大谱强度数据F(h),最大谱强度数据包含大气信息,其中h为高度;
步骤3,将谱强度数据等效至激光雷达方程,同时也进行了距离校正:
其中,A为谱强度数据同回波强度转换的比例关系,为定值;K为系统常数;h为高度;β(h)为大气总后向散射系数,α(ξ)是关于高度ξ的大气总消光系数;
步骤4,利用大气消光系数与大气后向散射系数之间的对应关系,将大气总后向散射系数由大气总消光系数代替:β(h)=αk(h)/S1,
S1为消光后向散射比,k与激光发射波长、气溶胶离子特性相关,0.67≤k≤1,一般取k=1;
对步骤3所得公式进行对数与微分处理,得到结果:
其中,步骤3公式的常数项做对数变为单独常数项,再微分计算时常数为0;
步骤5,结合伯努利方程解法,获取大气消光系数信息:
αf(h)为后向大气消光系数,α(hf)为其边界条件;
αb(h)为前向大气消光系数,α(hb)为其边界条件;
hf及hb分别为边界条件对应的高度值;
当高度大于等于hf时,大气消光系数α(h)取αf(h);
当高度小于hb时,大气消光系数α(h)取αb(h);
步骤6,在水平方向大气消光均匀的前提下,利用大气能见度同1550nm对应大气消光系数之间关系计算:
V(h)为不同高度的等效能见度,其中修正因子q为:
激光雷达公式里面有些较真系数是需要根据实际情况选取,而一般测量的也是实际情况,在实际操作时,若要统一到一个公式则可通过增加加跟实际情况相关的校正系数来实现。
对采样的时域信号转换为频域信号,其方法包括但不限于分段FFT变换、自相关分析等。
发明优点
相干激光雷达最近几年经过了比较快的发展,已具有了广泛的数据可做进一步研究。而总体而言,相干激光雷达应用仍处于商业发展初期,相干激光雷达已有研究的着重点,基本聚焦于对相干系统硬件性能的改进、不同测风模式的开发、风场测量结果的整合处理等方面,并未对相干回波中隐藏的大气信息进行提取分析。
本发明将常规相干激光雷达所忽略的信息充分利用,基于相干雷达对风测量中的计算流程,应用其中的中间处理数据,用谱强度数据来等效成大气回波强度,再用大气回波的反演方法计算能见度,从而实现了基于相干激光雷达谱强度数据的大气能见度的测量,不但扩展了相干激光雷达的测量功能,也为后续的研究打下基础。
由于相干激光雷达的探测原理,大气远场弱信号同本振信号拍频获取远场回波数据。相干激光雷达可利用较强的本振信号有效提升远场信号信噪比,相比传统气溶胶激光雷达对远场信号的探测更加精确。
相干雷达系统光学部件多采用成熟的集成化模块,系统本身稳定性和环境适应性相比传统的气溶胶雷达系统更有优势,更适合商业性的大范围推广。
附图说明
图1是相干激光雷达谱强度图,
横轴为点数,每100个点为一组,对应一个距离门的谱数据;纵轴为谱强度,对应距离门范围内最高点为对应距离门的谱强度数值。
图2是通过图1的数据获取不同距离门的等效大气回波强度图。
图3是本发明的流程图。
具体实施方式
相干激光雷达的探测过程为激光本振信号同雷达回波信号在空间混频探测。其中,激光本振信号由激光种子源引出一束有一定强度的窄带宽信号。激光发射脉冲在大气中的回波由相同光学器件接收,然后同本振光一起在空间上完成混频,进而由探测器将光信号转化为电信号。
其中,相干激光雷达系统本振信号为:
雷达所接收回波信号强度表示为:
两路信号在探测器上进行混频,光波矢量为:
因此,探测器所探测到的光强为:
Id(x,y,t)=|ud(x,y,t)|2=|uod|2+|usd|2+2Re[uod(x,y,t)*usd(x,y,t)]=Iod+Isd+In (1.4)
Iod=|uod|2
Isd=|usd|2
In=2Re[uod(x,y,t)*usd(x,y,t)]
其中,uod(x,y,t)*usd(x,y,t)为两个复数相乘,Re表示取其中的实部作为实际强度的有效数值。
Iod及Isd均为直流信号,因此可方便的滤除;而外差信号In为:
Δf=fod-fsd (1.6)
光电探测器输出信号使用高速模拟数字采集卡(AD采集卡),获得数字信号。光的传播速度确定,因此系统每个距离门内的一组数据进行FFT变化,获得不同径向高度对应的频率谱强度信号,如图1所示。FFT谱数据强度中,中心频率对应强度同IodIsd成正比,对应频率数值为Δf。
相干激光雷达仅使用频率中心频率相比发射相干中心频率的差值Δf,用于反演风速;而对应频率的强度信息则含有大气回波的信号Isd与本振光强度Iod。
根据以上分析,提取出的不同距离中心频率强度为本振光信号强度同回波信号强度的乘积。而本振光强度稳定,因此变化后的谱强度信息对应了普通激光探测的回波信号。同时,由于本振光强度较大,对应谱强度信号信噪比较高,利用谱强度数据进行大气能见度探测,系统探测距离及准确度相比普通激光雷达方法,优势明显,也具有充分的可靠性,接下来的主要工作就是如何采用准确的算法利用好这些数据来实现大气能见度的测量。
对于相干激光雷达获得的谱强度数据,首先提取出对应距离门的最大谱强度信号,经足够的时间累积平均,获取大小相合激光雷达方程的信号。
设获得的谱中心强度数据为F(h),则与实际大气参数之间关联性如下:
其中,A为谱强度数据同回波强度转换的比例关系,为定值;K为系统常数;h为对应的高度;βaer(h)及βmol(h)分别为大气分子后向散射系数及气溶胶后向散射系数;αaer(ξ)及αmol(ξ)分别为大气分子消光系数及气溶胶消光系数,ξ表示高度值。
获取的谱强度数据可对应激光雷达方程,参照一般的气溶胶反演方法进行反演。且相比普通的气溶胶激光雷达,谱强度数据信噪比较大,探测距离及精度有明显优势。
对于后向散射系数为β(h),消光系数为α(h),有以下公式:
β(h)=βaer(h)+βmol(h) (2.2)
α(h)=αaer(h)+αmol(h) (2.3)
谱强度信号进行距离校正后为X(h),
对公式2.4取对数得:
对公式2.5取微分:
在大气回波的特性中,后向散射系数同消光系数之间存在一定的对应关系,以如下公式进行表示:
β(h)=αk(h)/S1 (2.7)
S1为消光后向散射比,k与激光发射波长、气溶胶离子特性相关,0.67≤k≤1,一般取k=1;
将公式2.7带入2.6,
利用伯努利方程的求解方法,公式2.8可得:
式2.9为后向积分式,α(hf)为边界条件;式2.10为前向积分式,α(hb)为边界条件;边界条件α(hf)和α(hb)的求取,可以假设某一范围内大气均匀并利用斜率法得到。
在水平大气消光均匀的前提下,剩余能量ε:
ε=exp(-αλ(h)V(h)) (2.11)
V(h)为不同高度等效能见度,αλ(h)为波长λ对应大气消光系数,正常人眼的敏感亮度一般等效为550nm,平均亮度对比感阈值ε=0.02,等效关系为550nm光波能量降至0.02时,系统传输所经过的距离。因此获得水平能见度结果为:
V(h)为不同高度对应的等效能见度,α550(h)为等效550nm波长消光系数,h为对应高度。
一般相干激光雷达所用波长为1550nm,因此需将1550nm波长测量消光系数结果等效至550nm的人眼敏感波长:
其中修正因子q为:
根据以上关系可知,使用1550nm作为能见度测量波长时,需根据实际天气条件,对能见度进行计算。当能见度大于50公里,对应测量消光系数小于2.2×10-4/m时,修正因子取1.6。当能见度在6km~50km之间,对应消光系数为2.2×10-4/m~1.53×10-3/m时,修正因子取值为1.3。当能见度小于6km,消光系数大于1.53×10-3/m时,对应修正因子取值为第三种情况。其中,消光系数的计算是根据能见度的边界条件计算的,两者实际是一个条件。
Claims (2)
1.一种基于相干激光雷达谱强度数据的大气能见度分析方法,其特征在于包括以下步骤:
步骤1,对相干激光雷达的光电探测器采样信号进行变换,获得对应于不同高度的谱数据;
步骤2,对不同高度的谱数据进行强度检测,得出不同高度处最大谱强度数据F(h),最大谱强度数据包含大气信息,其中h为高度;
步骤3,将谱强度数据等效至激光雷达方程,同时也进行了距离校正:
其中,A为谱强度数据同回波强度转换的比例关系,为定值;K为系统常数;h为高度;β(h)为大气总后向散射系数,α(ξ)是关于高度ξ的大气总消光系数;
步骤4,利用大气消光系数与大气后向散射系数之间的对应关系,将大气总后向散射系数由大气总消光系数代替:β(h)=αk(h)/S1,
S1为消光后向散射比,k与激光发射波长、气溶胶离子特性相关,0.67≤k≤1,
对步骤3所得公式进行对数与微分处理,得到结果:
其中,步骤3公式的常数项做对数变为单独常数项,再微分计算时常数为0;
步骤5,结合伯努利方程解法,计算步骤4公式的未知函数α(h),而α(h)可结合边界条件分为下面两种表达公式:
αf(h)为后向大气消光系数,α(hf)为其边界条件;
αb(h)为前向大气消光系数,α(hb)为其边界条件;
hf及hb分别为边界条件对应的高度值;
当高度大于等于hf时,大气消光系数α(h)取αf(h);
当高度小于hb时,大气消光系数α(h)取αb(h);
步骤6,在水平方向大气消光均匀的前提下,利用大气能见度同1550nm对应大气消光系数之间关系计算:
V(h)为不同高度的等效能见度,其中修正因子q为:
2.如权利要求1所述的基于相干激光雷达谱强度数据的大气能见度分析方法,其特征在于所述步骤4中,与激光发射波长、气溶胶离子特性相关的参数k,取k=1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010309070.0A CN111458727B (zh) | 2020-04-19 | 2020-04-19 | 基于相干激光雷达谱强度数据的大气能见度分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010309070.0A CN111458727B (zh) | 2020-04-19 | 2020-04-19 | 基于相干激光雷达谱强度数据的大气能见度分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111458727A CN111458727A (zh) | 2020-07-28 |
CN111458727B true CN111458727B (zh) | 2023-03-28 |
Family
ID=71684611
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010309070.0A Active CN111458727B (zh) | 2020-04-19 | 2020-04-19 | 基于相干激光雷达谱强度数据的大气能见度分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111458727B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117930281B (zh) * | 2024-03-25 | 2024-05-24 | 南京信息工程大学 | 一种基于相干多普勒激光雷达后向散射的能见度反演方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006013062A1 (en) * | 2004-08-06 | 2006-02-09 | Alvise Francisco Baseggio | Method and system of detecting, monitoring and defining parameters for flying or stationary objects in airspace and outer space |
CN105866795A (zh) * | 2016-05-19 | 2016-08-17 | 中国科学技术大学 | 基于全光纤扫描f-p干涉仪的大动态风场探测激光雷达 |
CN107807366A (zh) * | 2017-10-30 | 2018-03-16 | 中国科学技术大学 | 一种大气能见度的计算方法、装置、雷达及系统 |
CN108445509A (zh) * | 2018-04-10 | 2018-08-24 | 中国科学技术大学 | 基于gpu的相干激光雷达信号处理方法 |
CN108490451A (zh) * | 2018-03-29 | 2018-09-04 | 中国民航大学 | 一种利用大气消光系数反演斜程能见度的方法 |
CN109827906A (zh) * | 2019-01-31 | 2019-05-31 | 西安理工大学 | 一种激光雷达斜程能见度的反演方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AUPR301401A0 (en) * | 2001-02-09 | 2001-03-08 | Commonwealth Scientific And Industrial Research Organisation | Lidar system and method |
-
2020
- 2020-04-19 CN CN202010309070.0A patent/CN111458727B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006013062A1 (en) * | 2004-08-06 | 2006-02-09 | Alvise Francisco Baseggio | Method and system of detecting, monitoring and defining parameters for flying or stationary objects in airspace and outer space |
CN105866795A (zh) * | 2016-05-19 | 2016-08-17 | 中国科学技术大学 | 基于全光纤扫描f-p干涉仪的大动态风场探测激光雷达 |
CN107807366A (zh) * | 2017-10-30 | 2018-03-16 | 中国科学技术大学 | 一种大气能见度的计算方法、装置、雷达及系统 |
CN108490451A (zh) * | 2018-03-29 | 2018-09-04 | 中国民航大学 | 一种利用大气消光系数反演斜程能见度的方法 |
CN108445509A (zh) * | 2018-04-10 | 2018-08-24 | 中国科学技术大学 | 基于gpu的相干激光雷达信号处理方法 |
CN109827906A (zh) * | 2019-01-31 | 2019-05-31 | 西安理工大学 | 一种激光雷达斜程能见度的反演方法 |
Non-Patent Citations (2)
Title |
---|
《On the origin,nature and uses of Antarctic ice-sheet radio-echo layering》;SIEGERT M J;《Prog Phys Geog》;19991231;第23卷(第2期);第159-179页 * |
防患于未然 汽车毫米波雷达;张伦维;《当代汽车》;20090501(第05期);第98-101页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111458727A (zh) | 2020-07-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108534811B (zh) | 一种短腔光纤法珀传感器的腔长解调算法 | |
CN110058258B (zh) | 一种基于混合型激光雷达的大气边界层探测方法 | |
Sroga et al. | Lidar measurement of wind velocity profiles in the boundary layer | |
May et al. | The accuracy of RASS temperature measurements | |
CN112285742B (zh) | 一种在频域内估算相干测风激光雷达载噪比的方法 | |
CN111458726B (zh) | 基于相干激光雷达谱强度数据的大气气溶胶分析方法 | |
McGill et al. | Comparison of two direct-detection Doppler lidar techniques | |
CN110132138A (zh) | 基于级联干涉仪的双扫频光源测距系统及方法 | |
CN114019525A (zh) | 一种基于光梳的高精度激光光谱测距方法 | |
CN111458727B (zh) | 基于相干激光雷达谱强度数据的大气能见度分析方法 | |
CN112505719A (zh) | 多普勒频率校正二次补偿的激光测风雷达测风方法及系统 | |
JP2837401B2 (ja) | 光学速度計プローブ | |
AU2021105440A4 (en) | Instrument and method for monitoring the soil moisture change by using GPS ground reflection signal | |
CN113447946B (zh) | 微弱激光回波信号的微多普勒信息测量系统 | |
CN103149180B (zh) | 一种土壤光谱反射率和电导率检测方法 | |
CN115453574B (zh) | 一种大气多参数探测多功能激光雷达 | |
CN115826088B (zh) | 中高层大气垂直风廓线的激光外差光谱测量方法及系统 | |
CN111750972A (zh) | 一种适用于双频连续激光雷达的目标振动测量方法 | |
CN111624614B (zh) | 一种通过激光相干探测进行测距的方法及系统 | |
CN115639572A (zh) | 一种星载风、温、气溶胶同步观测激光雷达 | |
Saklakova et al. | Field programmable gate array-based coherent lidar employing the ordinal statistics method for fast Doppler frequency determination | |
Lin et al. | Inter-comparison of wind measurements between Doppler wind lidar and L-band radiosonde | |
CN110568453A (zh) | 一种基于大气分子吸收的测风激光雷达及测风方法 | |
CN117930281B (zh) | 一种基于相干多普勒激光雷达后向散射的能见度反演方法 | |
Shen et al. | Rayleigh Doppler Lidar Technology Based on a Quadruple Dual-pass Fabry–Perot Interferometer |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |