CN106706566A - Calculation method of atmosphere vertical visibility via laser radar detection - Google Patents

Calculation method of atmosphere vertical visibility via laser radar detection Download PDF

Info

Publication number
CN106706566A
CN106706566A CN201710155843.2A CN201710155843A CN106706566A CN 106706566 A CN106706566 A CN 106706566A CN 201710155843 A CN201710155843 A CN 201710155843A CN 106706566 A CN106706566 A CN 106706566A
Authority
CN
China
Prior art keywords
point
laser radar
vertical
extinction coefficient
visibility
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.)
Granted
Application number
CN201710155843.2A
Other languages
Chinese (zh)
Other versions
CN106706566B (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

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Immunology (AREA)
  • Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Pathology (AREA)
  • Analytical Chemistry (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Electromagnetism (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Optical Radar Systems And Details Thereof (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

The invention discloses a calculation method of atmosphere vertical visibility via laser radar detection. The method comprises the following steps: obtaining an original echo signal for laser radar vertical measurement and judging whether cloud exists in the original echo signal; if no cloud signal exists, calculating the atmosphere vertical extinction coefficient by a Fernald method; removing the point with relatively great error on the atmosphere vertical extinction outline, and performing interpolation processing of the point; judging whether the extinction coefficient of an altitude point and at least 5 following consecutive points on the atmosphere vertical extinction outline is less than or equal to 0.1956, wherein the point is marked as the boundary altitude point n*L of atmosphere extinction coefficient; calculating the atmospheric transmissivity T1 between the measurement origin and the boundary altitude point n*L; inverting the average atmosphere extinction coefficient alpha between the measurement origin and the boundary altitude point n*L; and calculating the atmosphere vertical visibility V. In the invention, the algorithm is simple; and when no cloud exists above the laser radar equipment, the calculation result of atmosphere vertical visibility can be automatically provided in real time.

Description

A kind of computational methods of laser radar detection SEQUENCING VERTICAL visibility
Technical field
The present invention relates to the application field of laser radar, specifically a kind of meter of laser radar detection SEQUENCING VERTICAL visibility Calculation method.
Background technology
Visibility is one of main contents of meteorological element, and the visibility on SEQUENCING VERTICAL direction can both reflect atmosphere Degree of stability, also can be used as the important factor for judging property of air mass and Pollution Study.The visibility in SEQUENCING VERTICAL direction Measurement can be applied to the fields such as traffic, military affairs, science and technology, be that takeoff and landing and people's daily life provide necessary safety guarantee. At present, conventional visibility meter product is measured primarily directed to horizontal meteorological optical range, does not possess vertical visibility measurement Function.Nanjing Information engineering Univ foretells makes soldier, Qiu Zujing, Hou Yuyun et al. propose the forward direction using the horizontal meteorological optical range of measurement Scattering formula visibility meter and the Laser-ceilometer of vertical survey determine horizontal visibility value and air under the conditions of different weather respectively Aerosol vertical extinction coefficient α, sets up the corresponding relation of different visibility values and Aerosol Extinction Coefficients, eventually through anti- Drill transmitance that different distance highly locates to solve the vertical visibility of whole atmosphere, the method is present to be needed visibility meter and swash Light ceilometer repeatedly measures calibration simultaneously, and measurement data amount is big, error is readily incorporated, while algorithm is complex asking Topic.Laser radar is to physical effects such as scattering, absorption, the delustrings of laser, by quantitative analysis atmospheric aerosol according to air (floating dust, particulate matter, flue dust and sand and dust) backscattering echo that laser is produced is 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 directed to water at present Flat meteorological optical range, for the measurement of vertical visibility, there is presently no a kind of particularly effective computational methods.Therefore, it is necessary to Serious analysis are done from algorithm, a kind of computational methods method of effective laser radar detection SEQUENCING VERTICAL visibility is found, is Laser radar is applied to vertical visibility detection and provides a kind of new approaches.
The content of the invention
The present invention provides a kind of computational methods of laser radar detection SEQUENCING VERTICAL visibility, substantially former based on laser radar Reason, the backscattering echo produced to laser by atmospheric aerosol particle in quantitative analysis vertical direction, and then obtain air Aerosol vertical extinction coefficient profile, by calculating in vertical direction measurement point to the average extinction coefficient between boundary height point α, computing formula according to SEQUENCING VERTICAL visibility V realizes calculating the purpose of SEQUENCING VERTICAL visibility.
In order to solve the above-mentioned technical problem, the technical solution adopted by the present invention is：
A kind of computational methods of laser radar detection SEQUENCING VERTICAL visibility, comprise the following steps：
(1) original echoed signals of laser radar vertical survey, are obtained and judges whether there is cloud in original echoed signals；
(2), in cloudless signal, SEQUENCING VERTICAL extinction coefficient profile α is calculated using Fernald methods1*L、α2*L、 α3*L..., wherein L is the range resolution ratio of laser radar；
(3) the larger point of error on SEQUENCING VERTICAL delustring profile, is rejected, and interpolation processing is done to the point, wherein, it is specific to insert The SEQUENCING VERTICAL extinction coefficient that the method for value can take before and after the point does weighted average；
(4), judge on SEQUENCING VERTICAL extinction coefficient profile a certain height point and its at least continuously there are 5 delustrings of point afterwards Coefficient≤0.1956, the i.e. point and its afterwards 5 horizontal visibility V >=20000m of point, take the point for atmospheric extinction coefficient side Boundary highly point n*L；
(5), computation and measurement origin is 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 computing formula, SEQUENCING VERTICAL visibility V is calculated.
Preferably, atmospheric transmittance T in the step (5)1Computing formula it is as follows：
Wherein, α (r) is the SEQUENCING VERTICAL extinction coefficient at 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 computing formula of SEQUENCING VERTICAL visibility V is as follows in the step (7)：
V=3.912 α-1。
The beneficial effects of the invention are as follows：
The present invention provides a kind of computational methods 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 surveyed surveyed and being demarcated The error of introducing；Algorithm is simple, when there is no cloud directly over laser radar apparatus, can in real time, automatically provide SEQUENCING VERTICAL Visibility result of calculation；The range of application of laser radar is widened, for the observation of SEQUENCING VERTICAL visibility businessization provides a kind of new Thinking.
Brief description of the drawings
The present invention is further detailed explanation with reference to the accompanying drawings and detailed description.
Fig. 1 is rudimentary algorithm schematic flow sheet of the invention.
Specific embodiment
In order to illustrate more clearly of the specific embodiment of the invention or technical scheme of the prior art, below will be to specific The accompanying drawing to be used needed for implementation method or description of the prior art is briefly described, it should be apparent that, in describing below Accompanying drawing is some embodiments of the present invention, for those of ordinary skill in the art, before creative work is not paid Put, other accompanying drawings can also be obtained according to these accompanying drawings.
As shown in figure 1, the present invention provides a kind of computational methods of laser radar detection SEQUENCING VERTICAL visibility, including it is following Step：
(1) original echoed signals of laser radar vertical survey, are obtained and judges whether there is cloud in original echoed signals, by In the backscattering cross of the backscattering cross much larger than aerosol of cloud, so when laser transmits run into cloud in an atmosphere, Its backscattering echo will appear from the jump signal of rapid increase, be then now have cloud signal.Because 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 to laser It is emitted back towards ripple 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 methods1*L、α2*L、 α3*L..., wherein L is the range resolution ratio of laser radar；
(3) the larger point of error on SEQUENCING VERTICAL delustring profile, is rejected, and interpolation processing is done to the point, wherein, it is specific to insert The SEQUENCING VERTICAL extinction coefficient that the method for value can take before and after the point does weighted average；
(4), judge on SEQUENCING VERTICAL extinction coefficient profile a certain height point and its at least continuously there are 5 delustrings of point afterwards Coefficient≤0.1956, the i.e. point and its afterwards 5 horizontal visibility V >=20000m of point, take the point for atmospheric extinction coefficient side Boundary highly point n*L；
(5), computation and measurement origin is 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 computing formula, SEQUENCING VERTICAL visibility V is calculated.
Atmospheric transmittance T in the step (5)1Computing formula it is as follows：
Wherein, α (r) is the SEQUENCING VERTICAL extinction coefficient at r in formula.
In the step (6) origin to the Zenith Distance extinction coefficient α between boundary height point n*L computing formula such as Under：
The computing formula of SEQUENCING VERTICAL visibility V is as follows in the step (7)：
V=3.912 α-1。
Invention described above implementation method, is not intended to limit the scope of the present invention., any in the present invention Spirit and principle within modification, equivalent and the improvement made etc., should be included in claim protection model of the invention Within enclosing.

Claims (3)

1. a kind of computational methods of laser radar detection SEQUENCING VERTICAL visibility, it is characterised in that comprise the following steps：
(1) original echoed signals of laser radar vertical survey, are obtained and judges whether there is cloud in original echoed signals；
(2), in cloudless signal, SEQUENCING VERTICAL extinction coefficient profile α is calculated using Fernald methods1*L、α2*L、α3*L..., Wherein L is the range resolution ratio of laser radar；
(3) the larger point of error on SEQUENCING VERTICAL delustring profile, is rejected, and interpolation processing is done to the point, wherein, specific interpolation The SEQUENCING VERTICAL extinction coefficient that method can take before and after the point does weighted average；
(4), judge on SEQUENCING VERTICAL extinction coefficient profile a certain height point and its at least continuously there are 5 extinction coefficients of point afterwards ≤ 0.1956, the i.e. point and its afterwards 5 horizontal visibility V >=20000m of point, take the point for atmospheric extinction coefficient border is high Degree point n*L；
(5), computation and measurement origin is 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 computing formula, SEQUENCING VERTICAL visibility V is calculated.
2. computational methods of laser radar detection SEQUENCING VERTICAL visibility according to claim 1, it is characterised in that described Atmospheric transmittance T in step (5)1Computing formula it is as follows：
T 1 = exp ( - ∫ 0 n * l α ( r ) d r ) = exp ( - Σ i = 1 n ( α i * L * L ) )
Wherein, α (r) is the SEQUENCING VERTICAL extinction coefficient at r in formula.
The computational methods of laser radar detection SEQUENCING VERTICAL visibility according to claim 1, it is characterised in that the step Suddenly origin is as follows to the computing formula of the Zenith Distance extinction coefficient α between boundary height point n*L in (6)：
a = - L n ( T 1 ) n * L .
3. computational methods of laser radar detection SEQUENCING VERTICAL visibility according to claim 1, it is characterised in that described The computing formula of SEQUENCING VERTICAL visibility V is as follows in step (7)：
V=3.912 α-1。
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 true CN106706566A (en) 2017-05-24
CN106706566B 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)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108279221A (en) * 2017-12-07 2018-07-13 中国科学院国家天文台 A kind of part day area atmospheric transparency acquisition methods
CN108490451A (en) * 2018-03-29 2018-09-04 中国民航大学 A method of utilizing atmospheric extinction coefficient inverting slant visibility
CN108627812A (en) * 2018-05-28 2018-10-09 成都信息工程大学 A kind of laser radar atmospheric visibility measurement method and device
CN112505651A (en) * 2020-12-23 2021-03-16 北京遥测技术研究所 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
尹青 等: "激光雷达在气象和大气环境监测中的应用", 《气象与环境学报》 *
杨成武 等: "用于反演半导体激光云高仪垂直能见度的算法", 《强激光与粒子束》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108279221A (en) * 2017-12-07 2018-07-13 中国科学院国家天文台 A kind of part day area atmospheric transparency acquisition methods
CN108279221B (en) * 2017-12-07 2021-04-13 中国科学院国家天文台 Method for acquiring atmospheric transparency of local sky area
CN108490451A (en) * 2018-03-29 2018-09-04 中国民航大学 A method of utilizing atmospheric extinction coefficient inverting slant visibility
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
CN112505651A (en) * 2020-12-23 2021-03-16 北京遥测技术研究所 Automatic processing method for atmospheric detection laser radar
CN112505651B (en) * 2020-12-23 2022-06-21 北京遥测技术研究所 Automatic processing method for atmospheric detection laser radar

Also Published As

Publication number Publication date
CN106706566B (en) 2019-05-10

Similar Documents

Publication Publication Date Title
CN102183442B (en) Visibility measuring method based on atmospheric fine particle spectrometer
CN106706566A (en) Calculation method of atmosphere vertical visibility via laser radar detection
Pal et al. Elastic-backscatter-lidar-based characterization of the convective boundary layer and investigation of related statistics
CN110095784B (en) Ocean-low atmosphere laser transmission modeling method under influence of complex environment
CN109827906B (en) Inversion method of laser radar slope visibility
CN106446307B (en) Aerosol foundation data-based AOD (automated optical inspection) vertical correction effect evaluation method and system
Wang et al. Transport characteristics of Chinese haze over Northern Taiwan in winter, 2005–2014
CN108490451B (en) Method for inverting slope visibility by utilizing atmospheric extinction coefficient
Hon et al. Application of LIDAR‐derived eddy dissipation rate profiles in low‐level wind shear and turbulence alerts at H ong K ong I nternational A irport
CN110411927B (en) Atmospheric fine particle AOD and earth surface polarization reflectivity collaborative inversion method
CN108106979A (en) A kind of PM2.5 inversion methods merged based on MODIS and machine learning model
CN110910448B (en) Marine oil spill drift track tracing method
Chubb et al. Observations of high droplet number concentrations in Southern Ocean boundary layer clouds
Shi et al. Multiple technical observations of the atmospheric boundary layer structure of a red-alert haze episode in Beijing
Finn et al. Analysis of urban atmosphere plume concentration fluctuations
CN104819963A (en) Measurement method and monitoring system of atmosphere vertical visibility
Pueschel et al. Visibility and aerosol size frequency distribution
Ohneiser et al. Self-lofting of wildfire smoke in the troposphere and stratosphere: simulations and space lidar observations
CN108170927A (en) A kind of PM2.5 remote sensing inversion methods based on MODIS
CN106596355A (en) Correction method of extinction coefficient below low clouds in laser radar retrieval
Watkins et al. A review of the wind conditions experienced by a moving vehicle
CN105891078B (en) Sand devil sand dust mass concentration quantitative inversion evaluation method based on wind profile radar
CN107784261B (en) Road surface state detection method and device
Höpner et al. Investigation of two optical methods for aerosol‐type classification extended to a Northern Indian Ocean site
Wu et al. Development of a spectrum-based ship fuel sulfur content real-time evaluation method

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