CN111220517B - PM based on satellite remote sensing2.5Inversion method for components with different particle sizes - Google Patents
PM based on satellite remote sensing2.5Inversion method for components with different particle sizes Download PDFInfo
- Publication number
- CN111220517B CN111220517B CN201911240283.6A CN201911240283A CN111220517B CN 111220517 B CN111220517 B CN 111220517B CN 201911240283 A CN201911240283 A CN 201911240283A CN 111220517 B CN111220517 B CN 111220517B
- Authority
- CN
- China
- Prior art keywords
- aod
- components
- component
- particle size
- misr
- 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
- 239000002245 particle Substances 0.000 title claims abstract description 82
- 238000000034 method Methods 0.000 title claims abstract description 21
- 238000009826 distribution Methods 0.000 claims abstract description 20
- 230000000704 physical effect Effects 0.000 claims abstract description 7
- 239000000443 aerosol Substances 0.000 claims description 23
- 239000000203 mixture Substances 0.000 claims description 19
- 238000012360 testing method Methods 0.000 claims description 6
- 230000003287 optical effect Effects 0.000 claims description 5
- 239000013618 particulate matter Substances 0.000 claims description 5
- 238000010586 diagram Methods 0.000 claims description 4
- 230000002159 abnormal effect Effects 0.000 claims description 2
- 238000000605 extraction Methods 0.000 abstract 1
- 230000007613 environmental effect Effects 0.000 description 5
- 238000012544 monitoring process Methods 0.000 description 3
- 230000002265 prevention Effects 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000003915 air pollution Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000008821 health effect Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 210000002345 respiratory system Anatomy 0.000 description 1
- 239000000725 suspension Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/06—Investigating concentration of particle suspensions
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/02—Investigating particle size or size distribution
- G01N15/0205—Investigating particle size or size distribution by optical means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/06—Investigating concentration of particle suspensions
- G01N15/075—Investigating concentration of particle suspensions by optical means
Landscapes
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
The invention provides a PM based on satellite remote sensing2.5An inversion method of components with different particle sizes. The invention utilizes the physical property, concentration and other information of MISR AOD data and PM2.5Surface concentration data, extraction of different PM by MISR AOD2.5The proportion of the particle size components, and constructing a PM concentration inversion model with specific particle size to finish different PMs2.5The concentration of the particle size components is predicted, the data of the PM product with the specific particle size generated according to the method has good model precision, and the PM of a certain region can be truly reflected2.5Particle size distribution.
Description
Technical Field
The invention belongs to the technical field of air pollution prevention and control, and particularly relates to PM based on satellite remote sensing2.5An inversion method of components with different particle sizes.
Background
For forty years, along with the rapid increase of the economy of China, the pollution of dust haze in China gradually develops into a prominent environmental problem in China. In recent years, PM has been increasing with the increase of environmental consciousness of our country2.5Pollution has become one of the environmental issues of government and public concern, and PM2.5Effective acquisition of concentration information is to develop effective PM2.5The basis of environmental management.
PM2.5The pollution prevention and control of the method needs a fine management mode, and researches show that the suspension residence time of particles with smaller particle sizes in the air is longer, the particles are easier to be absorbed by human respiratory tracts, and the method possibly has stronger harmfulness. However, the current environmental air quality standard (GB3095-2012) in China is about PM2.5Only the standard value of the total mass concentration is considered, and the pollution influence caused by the particles with different particle size components is not considered. Due to PM2.5The health effects of particles with different particle sizes are different, and the existing management mode can cause the deviation of regional atmospheric pollution evaluation to a certain extent, so that effective information support cannot be provided for human health evaluation. On the other hand, the existing atmospheric environment monitoring network in China only has PM2.5Monitoring of total mass concentration, without monitoring of concentrations of components with different particle sizes, PM cannot be supported2.5Fine prevention and control.
Estimation of ground PM by Aerosol Optical Depth (AOD) using satellite remote sensing inversion2.5Is a new technology which is rapidly developed in recent years. However, at present, PM at home and abroad2.5PM (particulate matter) with multiple concerns in satellite remote sensing inversion research2.5For PM2.5The study of different particle size distributions is lacking.
Disclosure of Invention
In order to solve the technical problem, the invention provides a PM based on satellite remote sensing2.5An inversion method of components with different particle sizes. The following presents a simplified summary in order to provide a basic understanding of some aspects of the disclosed embodiments. This summary is not an extensive overview and is intended to neither identify key/critical elements nor delineate the scope of such embodiments. Its sole purpose is to present some concepts in a simplified form as a prelude to the more detailed description that is presented later.
The invention adopts the following technical scheme:
in some optional embodiments, a PM based on satellite remote sensing is provided2.5The inversion method of different particle size components comprises the following steps:
extracting the original MISR AOD mixture concentration data by using IDL and removing abnormal values;
obtaining the values of the AOD of eight specific components at characteristic time points and space points according to the extracted concentration data of the MISR AOD mixture and the known classification of the MISR AOD mixture which is successfully tested;
calculating the particle diameter relative mass distribution of the MISR eight AOD components;
drawing an AOD component particle size distribution diagram by utilizing the particle diameter relative mass distribution information of the MISR eight AOD components, and comparing the area occupied by the components with the particle size of less than x microns in each component with the area occupied by the components with the particle size of less than 2.5 microns;
calculating the component PM with particle size less than x microns at specific time points and space pointsxWith component PM having particle size of less than 2.5 μm2.5The comparison result of (1);
based on ground PM2.5The data of the concentration is stored in a memory,the particle size of the components PM with the specific time point and space point is less than x micronsxWith component PM having particle size of less than 2.5 μm2.5According to the specific time point, space point and ground PM2.5The concentration data are matched and multiplied to obtain PM of specific time points and space pointsxConcentration data.
In some optional embodiments, the method previously comprises: obtaining MISR aerosol product data; wherein the MISR aerosol product data comprises: concentration data of the aerosol and physical and photochemical properties data of the aerosol.
In some alternative embodiments, the process of calculating the particle diameter relative mass distribution of each AOD component comprises: and (3) calculating the particle diameter relative mass distribution of each AOD component according to parameter information provided by the obtained MISR aerosol product data by using the aerosol physical and optical property data of eight AODs in the MISR aerosol data and assuming that each AOD component conforms to the log-normal distribution and has constant particle density.
In some alternative embodiments, the values of the eight specific component AODs for the specific time point and spatial point are obtained by the following formula:
wherein, AODi,kIs the AOD value of AOD component i at MISR observation point k, AODjIs the total AOD value of AOD mixture j, Fi,jIndicates the contribution of AOD component i to AOD mixture j, and α ═ 1 indicates that mixture j passed the significance test, otherwise indicates that the test failed.
In some optional embodiments, the area occupied by the component with the particle size less than x microns is compared with the area occupied by the component with the particle size less than 2.5 microns in each component to obtain the comparison result
This result reflects the ratio of the concentrations of particulate matter of the respective particle sizes in AOD component i.
In some alternative embodiments, the component PM having a particle size less than x microns at a particular point in time and space is calculated from the following equationxWith component PM having particle size of less than 2.5 μm2.5The comparison result of (1):
wherein, AODiIs the AOD value of AOD component i;
the comparison results were obtained for the area of the components with a particle size of less than x microns compared to the area of the components with a particle size of less than 2.5 microns.
The invention has the following beneficial effects: the invention utilizes the physical property, concentration and other information of MISR AOD data and PM2.5Prediction of concentration by MISR AOD with different PMs2.5The proportion of the particle size components, and constructing a PM concentration inversion model with specific particle size to finish different PMs2.5The concentration of the particle size components is predicted, the data of the PM product with the specific particle size generated according to the method has good model precision, and the PM of a certain region can be truly reflected2.5Particle size distribution.
Drawings
FIG. 1 is a schematic flow diagram of the present invention.
Detailed Description
The following description sufficiently presents specific embodiments of the invention to enable those skilled in the art to practice them. Other embodiments may incorporate structural, logical, electrical, process, and other changes. The examples merely typify possible variations. Individual components and functions are optional unless explicitly required, and the sequence of operations may vary. Portions and features of some embodiments may be included in or substituted for those of others.
As shown in FIG. 1, in some illustrative embodiments, a PM based on satellite remote sensing is provided2.5Is differentParticle size component inversion method, physical property and concentration information of MISR AOD data and PM2.5Prediction of concentration by MISR AOD and PM2.5Constructing an inversion model by using the scale factors of different particle size components, which specifically comprises the following steps:
s1: MISR aerosol product data were acquired.
Wherein the MISR aerosol product data comprises: concentration data of the aerosol and physical and photochemical properties data of the aerosol.
The MISR aerosol product data can be downloaded by the atmospheric science data center of the NASA Lanli research center.
S2: the raw MISR AOD mixture concentration data was extracted using IDL and outliers were rejected.
S3: from the extracted MISR AOD mixture concentration data and the known successful inspection MISR AOD mixture classes, values for eight specific component AODs at specific time and space points were obtained:
wherein, AODi,kIs the AOD value of AOD component i at MISR observation point k, AODjIs the total AOD value of AOD mixture j, Fi,jIndicates the contribution of AOD component i to AOD mixture j, and α ═ 1 indicates that mixture j passed the significance test, otherwise indicates that the test failed.
S4: the particle diameter relative mass distributions of the eight MISR AOD components were calculated based on the physical and optical properties of the MISRAOD components.
Specifically, using the aerosol physical and optical property data of eight AODs in the MISR aerosol data, assuming that each AOD component conforms to a log-normal distribution and has a constant particle density, the particle diameter relative mass distribution of each AOD component is calculated based on the parametric information provided by the obtained MISR aerosol product data.
S5: utilizing the particle diameter relative mass distribution information of the MISR eight AOD components to draw the particle diameter distribution diagram of the AOD components and separating the particles in each componentComparing the area occupied by the component with the diameter less than x microns with the area occupied by the component with the diameter less than 2.5 microns, and recording the comparison result obtained by calculation as
This result reflects the ratio of the concentrations of particulate matter of the respective particle sizes in AOD component i.
S6: calculating the PM of the component having a particle diameter of less than x μm at a specific time point and space point based on the data calculated in step S5xWith component PM having particle size of less than 2.5 μm2.5The comparison result of (1):
wherein, AODiIs the AOD value of AOD component i;
the results are obtained by comparing the area occupied by the component with the particle size smaller than x with the area occupied by the component with the particle size smaller than 2.5 in each component.
S7: based on ground PM2.5Concentration data of the component PM having particle size of less than x μm at specific time point and space point calculated in step S6xWith component PM having particle size of less than 2.5 μm2.5According to the specific time point, space point and ground PM2.5The concentration data are matched and multiplied to obtain PM of specific time points and space pointsxConcentration data.
Those of skill would further appreciate that the various illustrative logical blocks, modules, circuits, and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, computer software, or combinations of both. To clearly illustrate this interchangeability of hardware and software, various illustrative components, blocks, modules, circuits, and steps have been described above generally in terms of their functionality. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the overall system. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present disclosure.
Claims (5)
1. PM based on satellite remote sensing2.5The inversion method of different particle size components comprises the following steps:
extracting the original MISR AOD mixture concentration data by using IDL and removing abnormal values;
obtaining the values of the AOD of eight specific components at characteristic time points and space points according to the extracted concentration data of the MISR AOD mixture and the known classification of the MISR AOD mixture which is successfully tested;
calculating the particle diameter relative mass distribution of the MISR eight AOD components;
drawing an AOD component particle size distribution diagram by utilizing the particle diameter relative mass distribution information of the MISR eight AOD components, and comparing the area occupied by the components with the particle size of less than x microns in each component with the area occupied by the components with the particle size of less than 2.5 microns;
calculating the component PM with particle size less than x microns at specific time points and space pointsxWith component PM having particle size of less than 2.5 μm2.5The comparison result of (1);
calculating the component PM with particle size less than x microns at specific time point and space point according to the following formulaxWith component PM having particle size of less than 2.5 μm2.5The comparison result of (1):
wherein, AODiIs the AOD value of AOD component i;
comparing the area occupied by the components with the particle size of less than x microns in each component with the area occupied by the components with the particle size of less than 2.5 microns to obtain a comparison result;
based on ground PM2.5Concentration data, composition with particle size less than x micron at the specific time point and space pointPMxWith component PM having particle size of less than 2.5 μm2.5According to the specific time point, space point and ground PM2.5The concentration data are matched and multiplied to obtain PM of specific time points and space pointsxConcentration data.
2. The PM based on satellite remote sensing according to claim 12.5The inversion method of components with different particle sizes is characterized in that data are required to be acquired before the method, and the acquired data comprise: MISR aerosol product data; wherein the MISR aerosol product data comprises: concentration data of the aerosol and physical and photochemical properties data of the aerosol.
3. PM based on satellite remote sensing according to claim 22.5The inversion method of components with different particle diameters is characterized in that the process of calculating the particle diameter relative mass distribution of each AOD component comprises the following steps: and (3) calculating the particle diameter relative mass distribution of each AOD component according to parameter information provided by the obtained MISR aerosol product data by using the aerosol physical and optical property data of eight AODs in the MISR aerosol data and assuming that each AOD component conforms to the log-normal distribution and has constant particle density.
4. PM based on satellite remote sensing according to claim 32.5The inversion method of components with different particle sizes is characterized in that the values of the AOD of the eight specific components at the specific time points and the specific space points are obtained by the following formula:
wherein, AODi,kIs the AOD value of AOD component i at MISR observation point k, AODjIs the total AOD value of AOD mixture j, Fi,jIndicates the contribution of AOD component i to AOD mixture j, and α ═ 1 indicates that mixture j passed the significance test, otherwise indicates that the test failed.
5. PM based on satellite remote sensing according to claim 42.5The inversion method of components with different particle sizes is characterized in that the comparison result obtained by comparing the area occupied by the components with the particle sizes of less than x microns in each component with the area occupied by the components with the particle sizes of less than 2.5 microns is
This result reflects the ratio of the concentrations of particulate matter of the respective particle sizes in AOD component i.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911240283.6A CN111220517B (en) | 2019-12-06 | 2019-12-06 | PM based on satellite remote sensing2.5Inversion method for components with different particle sizes |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911240283.6A CN111220517B (en) | 2019-12-06 | 2019-12-06 | PM based on satellite remote sensing2.5Inversion method for components with different particle sizes |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111220517A CN111220517A (en) | 2020-06-02 |
| CN111220517B true CN111220517B (en) | 2020-12-22 |
Family
ID=70829088
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911240283.6A Active CN111220517B (en) | 2019-12-06 | 2019-12-06 | PM based on satellite remote sensing2.5Inversion method for components with different particle sizes |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111220517B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109030301A (en) * | 2018-06-05 | 2018-12-18 | 中南林业科技大学 | Aerosol optical depth inversion method based on remotely-sensed data |
| WO2019065069A1 (en) * | 2017-09-27 | 2019-04-04 | 京セラ株式会社 | Measurement device component |
| CN110411918A (en) * | 2019-08-02 | 2019-11-05 | 中国科学院遥感与数字地球研究所 | A kind of PM2.5 concentration remote-sensing evaluation method based on satellite polarization technology |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107202750A (en) * | 2017-05-17 | 2017-09-26 | 河北中科遥感信息技术有限公司 | A kind of satellite-ground integrated monitoring quantitative remote sensing method for amalgamation processing of Atmospheric particulates |
| CN108426815B (en) * | 2018-04-20 | 2021-04-27 | 中国科学院遥感与数字地球研究所 | Method for estimating concentration of components of near-surface fine particulate matters |
| CN110160924B (en) * | 2019-06-27 | 2023-10-31 | 中国科学院遥感与数字地球研究所 | Particulate matter concentration detection method |
-
2019
- 2019-12-06 CN CN201911240283.6A patent/CN111220517B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2019065069A1 (en) * | 2017-09-27 | 2019-04-04 | 京セラ株式会社 | Measurement device component |
| CN109030301A (en) * | 2018-06-05 | 2018-12-18 | 中南林业科技大学 | Aerosol optical depth inversion method based on remotely-sensed data |
| CN110411918A (en) * | 2019-08-02 | 2019-11-05 | 中国科学院遥感与数字地球研究所 | A kind of PM2.5 concentration remote-sensing evaluation method based on satellite polarization technology |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111220517A (en) | 2020-06-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Liu et al. | The assessment of traffic accident risk based on grey relational analysis and fuzzy comprehensive evaluation method | |
| CN109784383B (en) | Rail crack identification method based on graph domain feature and DS evidence theory fusion | |
| KR20150086297A (en) | Method for characterising particles by image analysis | |
| CN104484410A (en) | Data fusion method and system applied to big data system | |
| CN113642666B (en) | Active enhancement soft measurement method based on sample expansion and screening | |
| CN108711266B (en) | Thunder and lightning based on atmospheric electric field is short to face local method for early warning | |
| CN111220517B (en) | PM based on satellite remote sensing2.5Inversion method for components with different particle sizes | |
| Ding et al. | The method of MEMS gyroscope random error compensation based on ARMA | |
| CN106338651A (en) | Particle filter analysis method applied to lower frequency oscillation mode identification of power system | |
| CN111141472B (en) | Anti-seismic support and hanger detection method and system | |
| CN110413949B (en) | Data processing method with increasing or decreasing trend | |
| Buryak et al. | Measurement of large-and superlarge-scale structures of the universe | |
| CN115964615A (en) | Track smoothness evaluation method based on centerline point cloud data | |
| Cautun et al. | Understanding the cosmic web | |
| CN105046056B (en) | The method of inspection and device of Cloud motion wind data | |
| Saghi et al. | Three-dimensional metrology and fractal analysis of dendritic nanostructures | |
| CN113255510A (en) | Rainbow cloud point charge moving path imaging method based on multi-time scale DBSCAN and sample entropy | |
| Sanders et al. | The size distributions of nanoscale Fe‐Ni‐S droplets in Stardust melted grains from comet 81P/Wild 2 | |
| Mehta et al. | Characterization of traffic and structure in the US airport network | |
| CN116429488B (en) | Method for inverting in-situ micro-regional chronology history of heterogeneous mineral and application | |
| Dai et al. | Intelligent image analysis technology and application for rail track inspection | |
| Fung et al. | Neural network modelling to estimate particle size distribution based on other particle sections and meteorological parameters 2 | |
| Fang et al. | Application of rough sets theory to sources analysis of atmospheric particulates | |
| Lu et al. | A parallel approach on clustering traffic data stream based on the density | |
| Asmi et al. | Ultrafine Aerosol Dynamics Module for Modal Aerosol Dynamic Model for Europe |
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 |