CN106706566B - A kind of calculation method of laser radar detection SEQUENCING VERTICAL visibility - Google Patents

A kind of calculation method of laser radar detection SEQUENCING VERTICAL visibility Download PDF

Info

Publication number
CN106706566B
CN106706566B CN201710155843.2A CN201710155843A CN106706566B CN 106706566 B CN106706566 B CN 106706566B CN 201710155843 A CN201710155843 A CN 201710155843A CN 106706566 B CN106706566 B CN 106706566B
Authority
CN
China
Prior art keywords
extinction coefficient
point
visibility
sequencing vertical
vertical
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
Application number
CN201710155843.2A
Other languages
Chinese (zh)
Other versions
CN106706566A (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.)
Hefei Institutes of Physical Science of CAS
Original Assignee
Hefei Institutes of Physical Science of CAS
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 Hefei Institutes of Physical Science of CAS filed Critical Hefei Institutes of Physical Science of CAS
Priority to CN201710155843.2A priority Critical patent/CN106706566B/en
Publication of CN106706566A publication Critical patent/CN106706566A/en
Application granted granted Critical
Publication of CN106706566B publication Critical patent/CN106706566B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • G01N21/49—Scattering, i.e. diffuse reflection within a body or fluid
    • G01N21/53—Scattering, i.e. diffuse reflection within a body or fluid within a flowing fluid, e.g. smoke
    • G01N21/538—Scattering, i.e. diffuse reflection within a body or fluid within a flowing fluid, e.g. smoke for determining atmospheric attenuation and visibility
    • 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
    • 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

Abstract

The present invention discloses a kind of calculation method of laser radar detection SEQUENCING VERTICAL visibility, comprising the following steps: obtains the original echoed signals of laser radar vertical measurement and judges whether there is cloud in original echoed signals；In cloudless signal, SEQUENCING VERTICAL extinction coefficient is calculated using Fernald method；The biggish point of error on SEQUENCING VERTICAL delustring profile is rejected, and interpolation processing is done to the point；Judge a certain height point on SEQUENCING VERTICAL extinction coefficient profile and its at least continuously have extinction coefficient≤0.1956 of 5 points later, taking the point is atmospheric extinction coefficient boundary height point n*L；Measurement origin is calculated to the atmospheric transmittance T between boundary height point n*L1；According to atmospheric transmittance T1, inverting measurement origin to the Zenith Distance extinction coefficient α between boundary height point n*L；SEQUENCING VERTICAL visibility V is calculated, inventive algorithm is simple, when not having cloud right above laser radar apparatus, can in real time, automatically provide SEQUENCING VERTICAL visibility calculated result.

Description

A kind of calculation method of laser radar detection SEQUENCING VERTICAL visibility
Technical field
The present invention relates to the application field of laser radar, the meter of specifically a kind of laser radar detection SEQUENCING VERTICAL visibility Calculation method.
Background technique
Visibility is one of main contents of meteorological element, and the visibility on SEQUENCING VERTICAL direction both can reflect atmosphere Degree of stability, also can be used as determine property of air mass and Pollution Study important factor.The visibility in SEQUENCING VERTICAL direction Measurement can be applied to the fields such as traffic, military affairs, science and technology, provide necessary safety guarantee for takeoff and landing and people's daily life. Currently, conventional visibility meter product is measured primarily directed to horizontal meteorological optical range, do not have vertical visibility measurement Function.Nanjing Information engineering Univ Bu Lingbing, Qiu Zujing, Hou Yuyun et al. are proposed using the forward direction for measuring horizontal meteorological optical range Scattering formula visibility meter and the Laser-ceilometer of vertical measurement measure horizontal visibility value and atmosphere under the conditions of different weather respectively Aerosol vertical extinction coefficient α establishes the corresponding relationship of different visibility values and Aerosol Extinction Coefficients, eventually by anti- The transmitance at different distance height is drilled to solve the vertical visibility of whole atmosphere, this method, which exists, to be needed visibility meter and swash Light ceilometer repeatedly measures calibration simultaneously, and measurement data amount is big, is readily incorporated error, while algorithm is more complicated asks Topic.Laser radar is to pass through quantitative analysis atmospheric aerosol to physical effects such as the scatterings, absorption, delustring of laser according to atmosphere The backscattering echo that (floating dust, particulate matter, flue dust and sand and dust) generate laser come remote sensing Aerosol Extinction Coefficients, The effective tool of mass concentration, cloud layer and visibility information.But its application in terms of visibility measurement is also to be directed to water at present Flat meteorological optical range, the measurement for vertical visibility, there is presently no a kind of particularly effective calculation methods.Therefore, it is necessary to Serious analysis is done from algorithm, is found a kind of calculation method method of effective laser radar detection SEQUENCING VERTICAL visibility, is Laser radar is applied to vertical visibility detection and provides a kind of new approaches.
Summary of the invention
The present invention provides a kind of calculation method of laser radar detection SEQUENCING VERTICAL visibility, substantially former based on laser radar Reason, the backscattering echo that laser is generated by atmospheric aerosol particle in quantitative analysis vertical direction, and then obtain atmosphere Aerosol vertical extinction coefficient profile passes through measurement point in calculating vertical direction to the average extinction coefficient between boundary height point α realizes the purpose for calculating SEQUENCING VERTICAL visibility according to the calculation formula of SEQUENCING VERTICAL visibility V.
In order to solve the above-mentioned technical problem, the technical solution adopted by the present invention is that:
A kind of calculation method of laser radar detection SEQUENCING VERTICAL visibility, comprising the following steps:
(1), the original echoed signals of laser radar vertical measurement are obtained and judge whether there is cloud in original echoed signals；
(2), in cloudless signal, SEQUENCING VERTICAL extinction coefficient profile α is calculated using Fernald method1*L、α2*L、 α3*L..., wherein L is the distance resolution of laser radar；
(3), the biggish point of error on SEQUENCING VERTICAL extinction coefficient profile is rejected, and interpolation processing is done to the point, wherein tool The method of body interpolation is that the SEQUENCING VERTICAL extinction coefficient before and after the point is taken to be weighted and averaged；
(4), judge a certain height point on SEQUENCING VERTICAL extinction coefficient profile and its at least continuously have the delustring of 5 points later Coefficient≤0.1956, i.e. the point and its later horizontal visibility V >=20000m of 5 points, taking the point is atmospheric extinction coefficient side Boundary height point n*L；
(5), measurement origin is calculated to the atmospheric transmittance T between boundary height point n*L1；
(6), according to atmospheric transmittance T1, inverting measurement origin to the Zenith Distance delustring system between boundary height point n*L Number α；
(7), according to calculation formula, SEQUENCING VERTICAL visibility V is calculated.
Preferably, atmospheric transmittance T in the step (5)1Calculation formula it is as follows:
Wherein, α (r) is SEQUENCING VERTICAL extinction coefficient at distance r in formula.
Preferably, in the step (6) origin to the Zenith Distance extinction coefficient α between boundary height point n*L calculating Formula is as follows:
Preferably, the calculation formula of SEQUENCING VERTICAL visibility V is as follows in the step (7):
V=3.912 α-1。
The beneficial effects of the present invention are:
The present invention provides a kind of calculation method of laser radar detection SEQUENCING VERTICAL visibility, without setting up visibility meter, When avoiding the atmospheric extinction coefficient that the visibility meter equipment that visibility meter under different visibility conditions is surveyed is surveyed and being demarcated The error of introducing；Algorithm is simple, when there is no cloud right above laser radar apparatus, can in real time, automatically provide SEQUENCING VERTICAL Visibility calculated result；The application range of laser radar has been widened, is provided for the observation of SEQUENCING VERTICAL visibility businessization a kind of new Thinking.
Detailed description of the invention
The present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments.
Fig. 1 is rudimentary algorithm flow diagram of the invention.
Specific embodiment
It, below will be to specific in order to illustrate more clearly of the specific embodiment of the invention or technical solution in the prior art Embodiment or attached drawing needed to be used in the description of the prior art be briefly described, it should be apparent that, it is described below Attached drawing is some embodiments of the present invention, for those of ordinary skill in the art, before not making the creative labor It puts, is also possible to obtain other drawings based on these drawings.
As shown in Figure 1, the present invention provides a kind of calculation method of laser radar detection SEQUENCING VERTICAL visibility, including following Step:
(1), the original echoed signals of laser radar vertical measurement are obtained and judge whether there is cloud in original echoed signals, by It is much larger than the backscattering cross of aerosol in the backscattering cross of cloud, so when laser transmits encounter cloud in an atmosphere, Its backscattering echo the jump signal increased rapidly will occur, at this time then be to have cloud signal.Since atmospheric aerosol echo is full Sufficient Mie scattering lidar equation is the basis for carrying out SEQUENCING VERTICAL visibility remote sensing survey.When there is cloud, cloud dissipates laser It is emitted back towards wave and is unsatisfactory for Mie scattering lidar equation, the present invention terminates to calculate；
(2), in cloudless signal, SEQUENCING VERTICAL extinction coefficient profile α is calculated using Fernald method1*L、α2*L、 α3*L..., wherein L is the distance resolution of laser radar；
(3), the biggish point of error on SEQUENCING VERTICAL extinction coefficient profile is rejected, and interpolation processing is done to the point, wherein tool The method of body interpolation is that the SEQUENCING VERTICAL extinction coefficient before and after the point is taken to be weighted and averaged；
(4), judge a certain height point on SEQUENCING VERTICAL extinction coefficient profile and its at least continuously have the delustring of 5 points later Coefficient≤0.1956, i.e. the point and its later horizontal visibility V >=20000m of 5 points, taking the point is atmospheric extinction coefficient side Boundary height point n*L；
(5), measurement origin is calculated to the atmospheric transmittance T between boundary height point n*L1；
(6), according to atmospheric transmittance T1, inverting measurement origin to the Zenith Distance delustring system between boundary height point n*L Number α；
(7), according to calculation formula, SEQUENCING VERTICAL visibility V is calculated.
Atmospheric transmittance T in the step (5)1Calculation formula it is as follows:
Wherein, α (r) is SEQUENCING VERTICAL extinction coefficient at distance r in formula.
In the step (6) origin to the Zenith Distance extinction coefficient α between boundary height point n*L calculation formula such as Under:
The calculation formula of SEQUENCING VERTICAL visibility V is as follows in the step (7):
V=3.912 α-1。
The embodiments of the present invention described above are not intended to limit the scope of the present invention, any in the present invention Spirit and principle within made modifications, equivalent substitutions and improvements etc., should be included in claim protection model of the invention Within enclosing.

Claims (2)

1. a kind of calculation method of laser radar detection SEQUENCING VERTICAL visibility, which comprises the following steps:
(1), the original echoed signals of laser radar vertical measurement are obtained and judge whether there is cloud in original echoed signals；
(2), in cloudless signal, SEQUENCING VERTICAL extinction coefficient profile α is calculated using Fernald method1*L、α2*L、α3*L..., Wherein L is the distance resolution of laser radar；
(3), the biggish point of error on SEQUENCING VERTICAL extinction coefficient profile is rejected, and interpolation processing is done to the point, wherein is specific to insert The method of value is that the SEQUENCING VERTICAL extinction coefficient before and after the point is taken to be weighted and averaged；
(4), judge a certain height point on SEQUENCING VERTICAL extinction coefficient profile and its at least continuously have the extinction coefficient of 5 points later ≤ 0.1956, i.e. the point and its later horizontal visibility V >=20000m of 5 points, taking the point is that atmospheric extinction coefficient boundary is high Spend point n*L；
(5), measurement origin is calculated to the atmospheric transmittance T between boundary height point n*L1；
(6), according to atmospheric transmittance T1, inverting measurement origin to the Zenith Distance extinction coefficient α between boundary height point n*L；
(7), according to calculation formula, SEQUENCING VERTICAL visibility V is calculated, the calculation formula of SEQUENCING VERTICAL visibility V is as follows:
V=3.912 α-1。
2. the calculation method of laser radar detection SEQUENCING VERTICAL visibility according to claim 1, which is characterized in that described Atmospheric transmittance T in step (5)1Calculation formula it is as follows:
Wherein, α (r) is SEQUENCING VERTICAL extinction coefficient at distance r in formula；
The calculation formula of origin to the Zenith Distance extinction coefficient α between boundary height point n*L are as follows in the step (6):
CN201710155843.2A 2017-03-16 2017-03-16 A kind of calculation method of laser radar detection SEQUENCING VERTICAL visibility Active CN106706566B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710155843.2A CN106706566B (en) 2017-03-16 2017-03-16 A kind of calculation method of laser radar detection SEQUENCING VERTICAL visibility

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710155843.2A CN106706566B (en) 2017-03-16 2017-03-16 A kind of calculation method of laser radar detection SEQUENCING VERTICAL visibility

Publications (2)

Publication Number Publication Date
CN106706566A CN106706566A (en) 2017-05-24
CN106706566B true CN106706566B (en) 2019-05-10

Family

ID=58887063

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710155843.2A Active CN106706566B (en) 2017-03-16 2017-03-16 A kind of calculation method of laser radar detection SEQUENCING VERTICAL visibility

Country Status (1)

Country Link
CN (1) CN106706566B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108279221B (en) * 2017-12-07 2021-04-13 中国科学院国家天文台 Method for acquiring atmospheric transparency of local sky area
CN108490451B (en) * 2018-03-29 2022-03-25 中国民航大学 Method for inverting slope visibility by utilizing atmospheric extinction coefficient
CN108627812A (en) * 2018-05-28 2018-10-09 成都信息工程大学 A kind of laser radar atmospheric visibility measurement method and device
CN112505651B (en) * 2020-12-23 2022-06-21 北京遥测技术研究所 Automatic processing method for atmospheric detection laser radar

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4419731A (en) * 1979-03-19 1983-12-06 The Secretary Of State For Defence In Her Britannic Majesty's Government Of The United Kingdom Of Great Britain And Northern Ireland Apparatus for estimating slant visibility in fog
CN1948953A (en) * 2006-10-19 2007-04-18 中国人民解放军海军工程大学 Based on laser rediffusion type atmospheric visibility testing method and testing instrument thereof
CN102590145A (en) * 2012-01-08 2012-07-18 杨少辰 Method for measuring atmospheric visibility by using laser and laser radar visibility meter
CN102621102A (en) * 2012-03-31 2012-08-01 中国科学院安徽光学精密机械研究所 Method for measuring horizontal visibility based on CCD (Charge Coupled Device) laser radar
CN104819963A (en) * 2015-05-07 2015-08-05 南京信息工程大学 Measurement method and monitoring system of atmosphere vertical visibility

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4419731A (en) * 1979-03-19 1983-12-06 The Secretary Of State For Defence In Her Britannic Majesty's Government Of The United Kingdom Of Great Britain And Northern Ireland Apparatus for estimating slant visibility in fog
CN1948953A (en) * 2006-10-19 2007-04-18 中国人民解放军海军工程大学 Based on laser rediffusion type atmospheric visibility testing method and testing instrument thereof
CN102590145A (en) * 2012-01-08 2012-07-18 杨少辰 Method for measuring atmospheric visibility by using laser and laser radar visibility meter
CN102621102A (en) * 2012-03-31 2012-08-01 中国科学院安徽光学精密机械研究所 Method for measuring horizontal visibility based on CCD (Charge Coupled Device) laser radar
CN104819963A (en) * 2015-05-07 2015-08-05 南京信息工程大学 Measurement method and monitoring system of atmosphere vertical visibility

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
激光雷达在气象和大气环境监测中的应用;尹青 等;《气象与环境学报》;20091031;第25卷(第5期);第48-56页 *
用于反演半导体激光云高仪垂直能见度的算法;杨成武 等;《强激光与粒子束》;20120229;第24卷(第2期);第307-311页 *

Also Published As

Publication number Publication date
CN106706566A (en) 2017-05-24

Similar Documents

Publication Publication Date Title
CN106706566B (en) A kind of calculation method of laser radar detection SEQUENCING VERTICAL visibility
KR101751642B1 (en) Method for correction of extinction coefficient obtained from atmospheric Light Detection And Ranging(LIDAR)
CN110006848B (en) Method and device for obtaining extinction coefficient of aerosol
CN110095784B (en) Ocean-low atmosphere laser transmission modeling method under influence of complex environment
CN109827906B (en) Inversion method of laser radar slope visibility
CN108490451B (en) Method for inverting slope visibility by utilizing atmospheric extinction coefficient
CN106291590A (en) The method calculating whole atmosphere aerosol optical depth based on lidar measurement data
CN104316443A (en) PM2.5 concentration monitoring method based on CCD back scattering
CN106383207A (en) Computing method applied to atmospheric aerosol mass concentration horizontal route distribution
Chubb et al. Observations of high droplet number concentrations in Southern Ocean boundary layer clouds
CN105044039B (en) A kind of method according to laser radar data automatic inversion horizontal visibility
CN104819963A (en) Measurement method and monitoring system of atmosphere vertical visibility
Cauteruccio et al. Wind tunnel validation of a particle tracking model to evaluate the wind‐induced bias of precipitation measurements
CN110006849B (en) Method and device for obtaining extinction coefficient of aerosol
CN105891078B (en) Sand devil sand dust mass concentration quantitative inversion evaluation method based on wind profile radar
CN109507072B (en) Fine particle turbulent flux measurement method
CN116879899A (en) Method based on aerial precipitation particle spectrum inversion
Mayor et al. Evaluating large-eddy simulations using volume imaging lidar data
CN105891064B (en) A kind of detection method and its device of aerosol particle composite rate
CN110057731A (en) Based on laser beam Intensity Analysis oceanic turbulence and particle cognitive method and device
CN115541461A (en) Particulate pollution tracing method
Voermans et al. Wave anomaly detection in wave measurements
CN105928881B (en) A kind of the window mirror delustring detection method and detection device of transmission-type visual range visibility meter
JP3571268B2 (en) Fog observation radar device
Chan Application of LIDAR backscattered power to visibility monitoring at the Hong Kong International Airport: some initial results

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant