CN108445507A - Aerosol particle size distribution distribution fitting method and system based on haze measurement data - Google Patents

Aerosol particle size distribution distribution fitting method and system based on haze measurement data Download PDF

Info

Publication number
CN108445507A
CN108445507A CN201810096141.6A CN201810096141A CN108445507A CN 108445507 A CN108445507 A CN 108445507A CN 201810096141 A CN201810096141 A CN 201810096141A CN 108445507 A CN108445507 A CN 108445507A
Authority
CN
China
Prior art keywords
lower atmosphere
layer
atmosphere layer
data
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
CN201810096141.6A
Other languages
Chinese (zh)
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 Institute of Environmental Features
Original Assignee
Beijing Institute of Environmental Features
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 Institute of Environmental Features filed Critical Beijing Institute of Environmental Features
Priority to CN201810096141.6A priority Critical patent/CN108445507A/en
Publication of CN108445507A publication Critical patent/CN108445507A/en
Pending legal-status Critical Current

Links

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
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • 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

Abstract

The present invention relates to a kind of aerosol particle size distribution distribution fitting method and system based on haze measurement data, wherein method includes the following steps:Aerosol optical characteristics measurement is carried out using laser radar and heliograph respectively, laser radar data is obtained and solar luminosity counts;The optical thickness of lower atmosphere layer is obtained according to the heliograph data inversion;The extinction coefficient of lower atmosphere layer is obtained according to the laser radar data inverting;The extinction coefficient of the optical thickness and lower atmosphere layer of the lower atmosphere layer based on inverting acquisition, the particulate Spectral structure of lower atmosphere layer is obtained using Retrieved by Extinction Method.The measurement data of present invention combination laser radar and heliograph, the layering of applied atmospheric optics thickness and Retrieved by Extinction Method Aerosol Size distribution characteristics, obtain the volume size distribution fitting function of lower atmosphere layer.

Description

Aerosol particle size distribution distribution fitting method and system based on haze measurement data
Technical field
The present invention relates to atmospheric optics technical field more particularly to a kind of aerosol particle size distributions based on haze measurement data Distribution fitting method and system.
Background technology
Haze is a kind of common weather phenomenon, it is since the total amount of the particulate matter, nitrogen oxides that contain in air is super Normal level is crossed, causes air muddy, visibility degree reduces.Influence and its particulate distribution of the haze to air have Direct relation.For haze weather, small particles concentration or big particle concentration either in air all compare normal weather Big is more.Compare greasy weather gas, and the aerosol small particles in haze weather are more than aerosol small particles concentration in greasy weather gas, and haze day The big particle concentration of aerosol is much smaller than the big particle concentration of aerosol in greasy weather gas in gas.
Based on heliograph measurement data, when carrying out Aerosol Size distribution inversion using light extinction method, due to solar luminosity The atmospheric extinction that meter obtains is whole atmosphere, so thus inverting obtains the Aerosol Size point that Aerosol Size is distributed as whole atmosphere Cloth, and haze is typically distributed across in low latitude, therefore researcher is more concerned with being aerosol in lower atmosphere layer.In lower atmosphere layer, Delustring is mainly influenced by aerosol, and aerosol is mainly gathered in 0-5km, and the delustring of upper atmosphere is mainly molecular extinction, therefore Atmospheric extinction at different height is caused by different reasons, and the Spectral structure research of low layer aerosol thus can be carried out.
Invention content
The technical problem to be solved in the present invention is, flood can only be obtained for currently based on heliograph measurement data The defect of the Aerosol Size distribution of air, provide a kind of aerosol particle size distribution distribution fitting method based on haze measurement data and System, applied atmospheric optics thickness layered approach and Retrieved by Extinction Method Aerosol Size distribution characteristics, obtain different aerosol particles Spectral structure fitting function when sub- size.
In order to solve the above-mentioned technical problem, it is molten to provide a kind of gas based on haze measurement data for first aspect present invention Glue volume size distribution approximating method, this method include:
Aerosol optical characteristics measurement is carried out using laser radar and heliograph respectively, obtain laser radar data and Solar luminosity counts;
The optical thickness of lower atmosphere layer is obtained according to the heliograph data inversion;
The extinction coefficient of lower atmosphere layer is obtained according to the laser radar data inverting;
The extinction coefficient of the optical thickness and lower atmosphere layer of the lower atmosphere layer based on inverting acquisition, it is anti-using light extinction method Drill the particulate Spectral structure for obtaining lower atmosphere layer.
In the aerosol particle size distribution distribution fitting method according to the present invention based on haze measurement data, preferably Ground, it is described according to the heliograph data inversion obtain lower atmosphere layer optical thickness the step of be specially:Based on uniform Air is divided into two layers by the hypothesis of parallel spherical surfaces air, and first layer is lower atmosphere layer, and the second layer is upper atmosphere, and sunray passes through First layer air path and second layer air path change with solar zenith angle;Using the data under multiple solar zenith angles This two layers of air path length is calculated by multiple linear regression analysis method;It is calculated separately based on the air path length Obtain the optical thickness of lower atmosphere layer and upper atmosphere.
In the aerosol particle size distribution distribution fitting method according to the present invention based on haze measurement data, preferably Ground, it is described according to the laser radar data inverting obtain lower atmosphere layer extinction coefficient the step of be specially:According to measurement Laser radar data obtains the extinction coefficient at different detection ranges using laser radar data Inversion Software;Utilize different detections Extinction coefficient at distance calculates the corresponding extinction coefficient of lower atmosphere layer.
In the aerosol particle size distribution distribution fitting method according to the present invention based on haze measurement data, preferably Ground, the extinction coefficient of the optical thickness and lower atmosphere layer of the lower atmosphere layer obtained based on inverting are anti-using light extinction method Drill obtain lower atmosphere layer particulate Spectral structure the step of be specially:By aerosol it is assumed that sphere, using Michaelis principle Establish the single particle optical property database of different scale;Establish the optical thickness of the lower atmosphere layer of heliograph data inversion With contacting for particulate Spectral structure, and the extinction coefficient and aerosol of the lower atmosphere layer of laser radar data inverting are established The contact of volume size distribution;Judge whether haze weather is caused by sand and dust, be, using normal distribution to Aerosol Size be distributed into Row description, is determined normal function parameter, is otherwise retouched to the particulate Spectral structure of lower atmosphere layer using Junge Spectral structures It states, determines Junge parameters.
Second aspect of the present invention provides a kind of aerosol particle size distribution fitting of distribution system based on haze measurement data, Including:
Data capture unit carries out aerosol optical characteristics measurement for obtaining using laser radar and heliograph Laser radar data and solar luminosity count;
First data processing unit, the optics for obtaining lower atmosphere layer according to the heliograph data inversion are thick Degree;
Second data processing unit, the extinction coefficient for obtaining lower atmosphere layer according to the laser radar data inverting;
Fitting of distribution unit, the delustring of the optical thickness and lower atmosphere layer of the lower atmosphere layer for being obtained based on inverting Coefficient obtains the particulate Spectral structure of lower atmosphere layer using Retrieved by Extinction Method.
In the aerosol particle size distribution fitting of distribution system according to the present invention based on haze measurement data, preferably Ground, the first data processing unit are specially according to the step of heliograph data inversion acquisition atmosphere optical thickness:Base Assume air being divided into two layers in uniform parallel Spherical Atmosphere, first layer is lower atmosphere layer, and the second layer is upper atmosphere, sunlight Line is changed by first layer air path and second layer air path with solar zenith angle, using under multiple solar zenith angles Data this two layers of air path length is calculated by multiple linear regression analysis method;Based on the air path length point The optical thickness of lower atmosphere layer and upper atmosphere Ji Suan not obtained.
In the aerosol particle size distribution fitting of distribution system according to the present invention based on haze measurement data, preferably Ground, the fitting of distribution unit:By aerosol it is assumed that sphere, the single particle optical that different scale is established using Michaelis principle are special Property database;The optical thickness of the lower atmosphere layer of heliograph data inversion and contacting for particulate Spectral structure are established, And establish the extinction coefficient of lower atmosphere layer and the contacting for particulate Spectral structure of laser radar data inverting;Judge haze Whether weather is caused by sand and dust, is, is described to Aerosol Size distribution using normal distribution, is determined normal function parameter, no Then the particulate Spectral structure of lower atmosphere layer is described using Junge Spectral structures, determines Junge parameters.
Implement the aerosol particle size distribution distribution fitting method and system based on haze measurement data of the present invention, has following Advantageous effect:The measurement data of present invention combination laser radar and heliograph, applied atmospheric optics thickness are layered and disappear Light method Retrieving Aerosol Size distribution characteristics, obtains the volume size distribution fitting function of lower atmosphere layer;Further, the present invention works as mist When haze weather is caused by sand and dust, the normal distribution fitting function of aerosol particle size distribution is obtained;When haze weather is that non-sand and dust cause When, obtain the Junge Spectral structure fitting functions of aerosol particle size distribution.
Description of the drawings
Fig. 1 is the aerosol particle size distribution distribution fitting method based on haze measurement data according to first embodiment of the invention Flow chart;
Fig. 2 is the aerosol particle size distribution distribution fitting method based on haze measurement data according to second embodiment of the invention Flow chart;
Fig. 3 is the aerosol particle size distribution fitting of distribution system based on haze measurement data according to the preferred embodiment of the present invention Module frame chart.
Specific implementation mode
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is A part of the embodiment of the present invention, instead of all the embodiments.Based on the embodiments of the present invention, ordinary skill people The every other embodiment that member is obtained without making creative work, shall fall within the protection scope of the present invention.
The present invention proposes a kind of aerosol particle size distribution distribution fitting method based on haze measurement data, this method application Atmosphere optical thickness hierarchical algorithm and Retrieved by Extinction Method Aerosol Size distribution characteristics, when obtaining different particulate sizes Spectral structure fitting function.
Referring to Fig. 1, for according to the particulate Spectral structure based on haze measurement data of first embodiment of the invention The flow chart of approximating method.As shown in Figure 1, the embodiment provides the aerosol particle size distribution fitting of distribution based on haze measurement data Method includes the following steps:
First, in step S101, data acquisition step is executed, when there is haze weather, using laser radar and too Sunlight degree meter carries out aerosol optical characteristics measurement respectively, obtains laser radar data and solar luminosity counts.
Then, in step s 102, the first data processing step is executed, is obtained according to the heliograph data inversion The optical thickness of lower atmosphere layer.
Then, in step s 103, the second data processing step is executed, is obtained according to the laser radar data inverting low The extinction coefficient of atmosphere.
Finally, the optical thickness of the lower atmosphere layer obtained based on inverting in step S104 and disappearing for lower atmosphere layer Backscatter extinction logarithmic ratio obtains the particulate Spectral structure of lower atmosphere layer using Retrieved by Extinction Method.
Fig. 2 is please referred to, for according to the aerosol particle size distribution based on haze measurement data of second embodiment of the invention The flow chart of distribution fitting method.As shown in Fig. 2, the particulate based on haze measurement data that the second embodiment provides Spectral structure approximating method specifically includes following steps:
First, data acquisition step includes step S201 and step S205;The two steps arbitrarily can successively execute or Person is performed simultaneously.
In step s 201, aerosol optical characteristics measurement is carried out using heliograph, obtains solar luminosity and counts;
In step S205, aerosol optical characteristics measurement is carried out using laser radar, obtains laser radar data;
Then, in step S202~S204, the first data processing step is executed, is counted instead according to the solar luminosity Drill the optical thickness for obtaining lower atmosphere layer;It specifically includes:
In step S202, air is layered:Assume air being divided into two layers based on uniform parallel Spherical Atmosphere, the One layer is lower atmosphere layer, and atmospheric extinction is mainly influenced by aerosol, and the second layer is upper atmosphere, and atmospheric extinction is mainly by molecule shadow It rings.Sunray is changed by first layer air path and second layer air path with solar zenith angle.
In step S203, this is calculated by multiple linear regression analysis method using the data under multiple solar zenith angles Two layers of air path length.After determining lower atmosphere layer and upper atmosphere respective air path length, so that it may to divide Lower atmosphere layer and upper atmosphere.The air path length of lower atmosphere layer is usually 5km or so, slightly not according to same day weather condition Together.
In step S204, the air path length based on two atmospheres, so that it may with calculate separately obtain lower atmosphere layer and The optical thickness of upper atmosphere.Since count acquisition based on the solar luminosity that aerosol optical characteristics measure is whole atmosphere Optical thickness, therefore the light of wherein lower atmosphere layer can be obtained in the present invention after lower atmosphere layer and upper atmosphere by dividing Thickness is learned to use for subsequent step.
In step S206, the second data processing step is executed, it is big to obtain low layer according to the laser radar data inverting The extinction coefficient of gas, specifically includes:According to the laser radar data of measurement difference is obtained using laser radar data Inversion Software Extinction coefficient at detection range;The corresponding delustring system of lower atmosphere layer is calculated followed by the extinction coefficient at different detection ranges Number.
In step S207~S210, fitting of distribution step is executed, is specifically included:
In step S207, the particulate Spectral structure of lower atmosphere layer is obtained using Retrieved by Extinction Method;Step S207 Including:
1) aerosol is assumed to sphere, the single particle optical characteristic quantity of different scale is established using Michaelis (Mie) principle According to library;
2) contacting for heliograph and lidar measurement data inversion result and particulate Spectral structure is established:
The formula mainly utilized is:
The calculation formula of extinction coefficient σ (λ) is as follows:
Wherein, r is particulate radius, and λ is wavelength, QeFor extinction efficiency factor, by it is aforementioned 1) in the simple grain established Sub-light property database determines;
The calculation formula of opticalthicknessτ is as follows:
Wherein, z is air height.
It is combined by above formula, opticalthicknessτ and the gas for establishing the lower atmosphere layer of heliograph data inversion are molten The contact of glue volume size distribution n (r), and establish the lower atmosphere layer of laser radar data inverting extinction coefficient σ (λ) and gas it is molten The contact of glue volume size distribution n (r);
Specifically, the opticalthicknessτ of lower atmosphere layer is established by following formula (3) with particulate Spectral structure n's (r) Contact:
Specifically, established by aforementioned formula (1) the extinction coefficient σ (λ) of the lower atmosphere layer of laser radar data inverting with The contact of particulate Spectral structure n (r).
In step S208, judge whether haze weather is caused by sand and dust, be, goes to step S210, otherwise go to step S209;The present invention can judge whether haze weather is caused by sand and dust by the weather forecast information on the day of obtaining.
In step S209, the particulate Spectral structure of lower atmosphere layer is described using Junge Spectral structures, is determined Junge parameters.
By the formula (4) of following Junge Spectral structures to the particulate Spectral structure of lower atmosphere layer in step S209 It is described, and contact and extinction coefficient σ (λ) and aerosol of the combination opticalthicknessτ with particulate Spectral structure n (r) The contact of volume size distribution n (r), it is common to solve Junge parameters:
In formula, r assumes that the radius of the particulate for spherical shape, υ are Junge parameters, and A is dependent on aerosol concentration Constant.
In step S210, Aerosol Size distribution is described using normal distribution, determines normal function parameter.
The particulate Spectral structure of lower atmosphere layer is carried out by the formula (5) of following normal distribution in step S210 Description, and contact and extinction coefficient σ (λ) and particulate of the combination opticalthicknessτ with particulate Spectral structure n (r) The contact of Spectral structure n (r), it is common to solve normal function parameter:
In formula, N0Indicate particulate sum in unit volume,Indicate that the mean radius of particle, σ are joined for normal function Number, characterizes the variance of particulate radius.
Referring to Fig. 3, for according to the particulate Spectral structure based on haze measurement data of the preferred embodiment of the present invention The module frame chart of fitting system.As shown in figure 3, the system 300 that the embodiment provides includes:Data capture unit 301, first counts According to processing unit 302, the second data processing unit 303 and fitting of distribution unit 304.
Wherein, data capture unit 301 carries out aerosol optical spy for obtaining using laser radar and heliograph Property measure laser radar data and solar luminosity count.
First data processing unit 302 is connected with data capture unit 301, for being counted instead according to the solar luminosity Drill the optical thickness for obtaining lower atmosphere layer.First data processing unit 302 is handled at process and aforementioned first data of data Reason step is consistent, and details are not described herein.
Second data processing unit 303 is connected with data capture unit 301, for according to the laser radar data inverting Obtain the extinction coefficient of lower atmosphere layer.Second data processing unit 303 handles the process of data and aforementioned second data processing Step is consistent, and details are not described herein.
The optical thickness of the lower atmosphere layer and disappearing for lower atmosphere layer that fitting of distribution unit 304 is used to obtain based on inverting Backscatter extinction logarithmic ratio obtains the particulate Spectral structure of lower atmosphere layer using Retrieved by Extinction Method.The fitting of distribution unit 304 and above-mentioned side Fitting of distribution step is consistent in method, and details are not described herein.
In conclusion the present invention obtains gas in lower atmosphere layer not merely with the atmospheric aerosol measurement data of haze weather The Spectral structure of colloidal sol;And when haze weather is caused by sand and dust, the normal distribution fitting function of aerosol particle size distribution is obtained;When When haze weather is that non-sand and dust cause, the Junge Spectral structure fitting functions of aerosol particle size distribution are obtained.
It should be appreciated that aerosol particle size distribution distribution fitting method and system based on haze measurement data in the present invention Principle it is identical, therefore elaborating to the embodiment of the aerosol particle size distribution distribution fitting method based on haze measurement data It is also applied for the aerosol particle size distribution fitting of distribution system based on haze measurement data.
Finally it should be noted that:The above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations;Although Present invention has been described in detail with reference to the aforementioned embodiments, it will be understood by those of ordinary skill in the art that:It still may be used With technical scheme described in the above embodiments is modified or equivalent replacement of some of the technical features; And these modifications or replacements, various embodiments of the present invention technical solution that it does not separate the essence of the corresponding technical solution spirit and Range.

Claims (10)

1. a kind of aerosol particle size distribution distribution fitting method based on haze measurement data, which is characterized in that this method includes:
Aerosol optical characteristics measurement is carried out using laser radar and heliograph respectively, obtains laser radar data and the sun Luminometer number evidence;
The optical thickness of lower atmosphere layer is obtained according to the heliograph data inversion;
The extinction coefficient of lower atmosphere layer is obtained according to the laser radar data inverting;
The extinction coefficient of the optical thickness and lower atmosphere layer of the lower atmosphere layer based on inverting acquisition, is obtained using Retrieved by Extinction Method Obtain the particulate Spectral structure of lower atmosphere layer.
2. the aerosol particle size distribution distribution fitting method according to claim 1 based on haze measurement data, feature exist In, it is described according to the heliograph data inversion obtain lower atmosphere layer optical thickness the step of be specially:
Assuming air being divided into two layers based on uniform parallel Spherical Atmosphere, first layer is lower atmosphere layer, and the second layer is upper atmosphere, Sunray is changed by first layer air path and second layer air path with solar zenith angle;Using multiple sun days This two layers of air path length is calculated by multiple linear regression analysis method in data under apex angle;
The optical thickness for obtaining lower atmosphere layer and upper atmosphere is calculated separately based on the air path length.
3. the aerosol particle size distribution distribution fitting method according to claim 1 based on haze measurement data, feature exist In, it is described according to the laser radar data inverting obtain lower atmosphere layer extinction coefficient the step of be specially:
According to the laser radar data of measurement the delustring system at different detection ranges is obtained using laser radar data Inversion Software Number;
The corresponding extinction coefficient of lower atmosphere layer is calculated using the extinction coefficient at different detection ranges.
4. the aerosol particle size distribution fitting of distribution side described in any one of claim 1 to 3 based on haze measurement data Method, which is characterized in that the extinction coefficient of the optical thickness and lower atmosphere layer of the lower atmosphere layer obtained based on inverting is adopted It is specially with the step of particulate Spectral structure of Retrieved by Extinction Method acquisition lower atmosphere layer:
By aerosol it is assumed that sphere, the single particle optical property database of different scale is established using Michaelis principle;
Establish the optical thickness of the lower atmosphere layer of heliograph data inversion and contacting for particulate Spectral structure, Yi Jijian The extinction coefficient of lower atmosphere layer of vertical laser radar data inverting and contacting for particulate Spectral structure;
Judge whether haze weather is caused by sand and dust, be, Aerosol Size distribution is described using normal distribution, is determined just Otherwise state function parameter is described the particulate Spectral structure of lower atmosphere layer using Junge Spectral structures, determines Junge Parameter.
5. the aerosol particle size distribution distribution fitting method according to claim 4 based on haze measurement data, feature exist In by following formula, establishing the opticalthicknessτ and particulate Spectral structure of the lower atmosphere layer of heliograph data inversion The contact of n (r), and establish the extinction coefficient σ (λ) and particulate Spectral structure of the lower atmosphere layer of laser radar data inverting The contact of n (r);
Wherein, σ (λ) is extinction coefficient, and r is particulate radius, and λ is wavelength, QeFor extinction efficiency factor, z is that air is high Degree.
6. the aerosol particle size distribution distribution fitting method according to claim 5 based on haze measurement data, feature exist In, it is described that the particulate Spectral structure of lower atmosphere layer is described using Junge Spectral structures, determine the step of Junge parameters It is rapid to be specially:
The particulate Spectral structure of lower atmosphere layer is described by the formula of following Junge Spectral structures, and combines optics Thicl ness T and the contact of particulate Spectral structure n (r) and the connection of extinction coefficient σ (λ) and particulate Spectral structure n (r) System, it is common to solve Junge parameters:
In formula, r assumes that the radius of the particulate for spherical shape, υ are Junge parameters, and A is dependent on the normal of aerosol concentration Number.
7. the aerosol particle size distribution distribution fitting method according to claim 5 based on haze measurement data, feature exist In the use normal distribution is distributed the step of being described, determining normal function parameter to Aerosol Size and is specially:
The particulate Spectral structure of lower atmosphere layer is described by the formula of following normal distribution, and combines optical thickness τ and the contact of particulate Spectral structure n (r) and the contact of extinction coefficient σ (λ) and particulate Spectral structure n (r), altogether With solution normal function parameter:
In formula, N0Indicate particulate sum in unit volume,Indicate that the mean radius of particle, σ are normal function parameter, table The variance of particulate radius is levied.
8. a kind of aerosol particle size distribution fitting of distribution system based on haze measurement data, which is characterized in that including:
Data capture unit, for obtaining the laser for carrying out aerosol optical characteristics measurement using laser radar and heliograph Radar data and solar luminosity count;
First data processing unit, the optical thickness for obtaining lower atmosphere layer according to the heliograph data inversion;
Second data processing unit, the extinction coefficient for obtaining lower atmosphere layer according to the laser radar data inverting;
Fitting of distribution unit, the delustring system of the optical thickness and lower atmosphere layer of the lower atmosphere layer for being obtained based on inverting Number obtains the particulate Spectral structure of lower atmosphere layer using Retrieved by Extinction Method.
9. the aerosol particle size distribution fitting of distribution system according to claim 8 based on haze measurement data, feature exist In the first data processing unit is specially according to the step of heliograph data inversion acquisition atmosphere optical thickness:Base Assume air being divided into two layers in uniform parallel Spherical Atmosphere, first layer is lower atmosphere layer, and the second layer is upper atmosphere, sunlight Line is changed by first layer air path and second layer air path with solar zenith angle, using under multiple solar zenith angles Data this two layers of air path length is calculated by multiple linear regression analysis method;Based on the air path length point The optical thickness of lower atmosphere layer and upper atmosphere Ji Suan not obtained.
10. the aerosol particle size distribution fitting of distribution system according to claim 8 based on haze measurement data, feature exist In the fitting of distribution unit:
By aerosol it is assumed that sphere, the single particle optical property database of different scale is established using Michaelis principle;
Establish the optical thickness of the lower atmosphere layer of heliograph data inversion and contacting for particulate Spectral structure, Yi Jijian The extinction coefficient of lower atmosphere layer of vertical laser radar data inverting and contacting for particulate Spectral structure;
Judge whether haze weather is caused by sand and dust, be, Aerosol Size distribution is described using normal distribution, is determined just Otherwise state function parameter is described the particulate Spectral structure of lower atmosphere layer using Junge Spectral structures, determines Junge Parameter.
CN201810096141.6A 2018-01-31 2018-01-31 Aerosol particle size distribution distribution fitting method and system based on haze measurement data Pending CN108445507A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810096141.6A CN108445507A (en) 2018-01-31 2018-01-31 Aerosol particle size distribution distribution fitting method and system based on haze measurement data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810096141.6A CN108445507A (en) 2018-01-31 2018-01-31 Aerosol particle size distribution distribution fitting method and system based on haze measurement data

Publications (1)

Publication Number Publication Date
CN108445507A true CN108445507A (en) 2018-08-24

Family

ID=63191476

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810096141.6A Pending CN108445507A (en) 2018-01-31 2018-01-31 Aerosol particle size distribution distribution fitting method and system based on haze measurement data

Country Status (1)

Country Link
CN (1) CN108445507A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111999267A (en) * 2020-08-19 2020-11-27 北京环境特性研究所 Particle average extinction efficiency factor fitting method and device
CN115201074A (en) * 2022-06-30 2022-10-18 中国科学院大气物理研究所 Method, system, device and computer-readable storage medium for remote sensing inversion of aerosol component distribution

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101504353A (en) * 2008-02-06 2009-08-12 香港科技大学 Method and system for providing near ground suspending particulate distribution
CN102539336A (en) * 2011-02-01 2012-07-04 环境保护部卫星环境应用中心 Method and system for estimating inhalable particles based on HJ-1 satellite
CN102944503A (en) * 2012-10-11 2013-02-27 中国科学院安徽光学精密机械研究所 PM2.5 mass concentration value automatic inversion algorithm based on sun photometer and lidar
CN103048653A (en) * 2012-12-13 2013-04-17 中国人民解放军63655部队 Micro pulse lidar system constant calibration method
CN103175759A (en) * 2013-02-25 2013-06-26 中国科学院安徽光学精密机械研究所 Method for acquiring complex refractive index of urban aerosol on basis of various ground-based remote sensing technologies
CN105891078A (en) * 2016-06-08 2016-08-24 中国气象局乌鲁木齐沙漠气象研究所 Quantitative inversion estimation method for sand dust storm sand dust mass concentration on basis of wind profile radar
CN106446307A (en) * 2015-08-05 2017-02-22 中国科学院遥感与数字地球研究所 Aerosol foundation data-based AOD vertical correction effect evaluation method and system
CN107219193A (en) * 2017-05-08 2017-09-29 西安电子科技大学 The inversion method of Atmosphere Refractivity Profile

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101504353A (en) * 2008-02-06 2009-08-12 香港科技大学 Method and system for providing near ground suspending particulate distribution
CN102539336A (en) * 2011-02-01 2012-07-04 环境保护部卫星环境应用中心 Method and system for estimating inhalable particles based on HJ-1 satellite
CN102944503A (en) * 2012-10-11 2013-02-27 中国科学院安徽光学精密机械研究所 PM2.5 mass concentration value automatic inversion algorithm based on sun photometer and lidar
CN103048653A (en) * 2012-12-13 2013-04-17 中国人民解放军63655部队 Micro pulse lidar system constant calibration method
CN103175759A (en) * 2013-02-25 2013-06-26 中国科学院安徽光学精密机械研究所 Method for acquiring complex refractive index of urban aerosol on basis of various ground-based remote sensing technologies
CN106446307A (en) * 2015-08-05 2017-02-22 中国科学院遥感与数字地球研究所 Aerosol foundation data-based AOD vertical correction effect evaluation method and system
CN105891078A (en) * 2016-06-08 2016-08-24 中国气象局乌鲁木齐沙漠气象研究所 Quantitative inversion estimation method for sand dust storm sand dust mass concentration on basis of wind profile radar
CN107219193A (en) * 2017-05-08 2017-09-29 西安电子科技大学 The inversion method of Atmosphere Refractivity Profile

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘吉等: ""兰州城区冬季大气气溶胶粒子谱的反演研究"", 《高原气象》 *
李放等: ""边界层气溶胶垂直分布阶梯遥感"", 《遥感技术与应用》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111999267A (en) * 2020-08-19 2020-11-27 北京环境特性研究所 Particle average extinction efficiency factor fitting method and device
CN115201074A (en) * 2022-06-30 2022-10-18 中国科学院大气物理研究所 Method, system, device and computer-readable storage medium for remote sensing inversion of aerosol component distribution

Similar Documents

Publication Publication Date Title
Bäumer et al. Relationship of visibility, aerosol optical thickness and aerosol size distribution in an ageing air mass over South-West Germany
Dayan et al. The influence of meteorological conditions and atmospheric circulation types on PM10 and visibility in Tel Aviv
Ma et al. Aerosol optical properties in the North China Plain during HaChi campaign: an in-situ optical closure study
CN101281142B (en) Method for measuring atmosphere visibility
Pu et al. Impact of long-range transport on aerosol properties at a regional background station in Northern China
Reddy et al. Seasonal aerosol variations over a coastal city, Zhoushan, China from CALIPSO observations
Liu et al. Observations of aerosol color ratio and depolarization ratio over Wuhan
CN110411927B (en) Atmospheric fine particle AOD and earth surface polarization reflectivity collaborative inversion method
CN107561554A (en) Inversion method with multi-wavelength laser radar data is counted based on solar luminosity
CN102183442A (en) Visibility measuring method based on atmospheric fine particle spectrometer
Li et al. Vertical distribution of particulate matter and its relationship with planetary boundary layer structure in Shenyang, Northeast China
Charlson et al. Background aerosol light scattering characteristics: nephelometric observations at Mauna Loa Observatory compared with results at other remote locations
CN110058258A (en) A kind of atmospheric boundary layer detection method based on mixed type laser radar
CN108445507A (en) Aerosol particle size distribution distribution fitting method and system based on haze measurement data
CN104330388A (en) Haze detector and haze detection method
Zheng et al. Aerosol optical properties observation and its relationship to meteorological conditions and emission during the Chinese National Day and Spring Festival holiday in Beijing
Barrett et al. Wet deposition of black carbon: A synthesis
Wandinger et al. Air mass modification over Europe: EARLINET aerosol observations from Wales to Belarus
Lee et al. Analytical approach to estimating aerosol extinction and visibility from satellite observations
Nakajima et al. Consistency of aerosol size distributions inferred from measurements of solar radiation and aerosols
KR20170025709A (en) Dust and haze prediction verification method and system
Gilardoni et al. Using measurements for evaluation of black carbon modeling
CN110907318B (en) Near-ground atmospheric total suspended particulate matter mass concentration remote sensing physical estimation method
Peng et al. Visibility measurements using two-angle forward scattering by liquid droplets
Xue et al. Long-term aerosol optical depth datasets over China retrieved from satellite data

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: 20180824