CN112525787A - Method for inverting PM2.5 all-weather fine grid data based on surface high-density meteorological data - Google Patents
Method for inverting PM2.5 all-weather fine grid data based on surface high-density meteorological data Download PDFInfo
- Publication number
- CN112525787A CN112525787A CN202011362233.8A CN202011362233A CN112525787A CN 112525787 A CN112525787 A CN 112525787A CN 202011362233 A CN202011362233 A CN 202011362233A CN 112525787 A CN112525787 A CN 112525787A
- Authority
- CN
- China
- Prior art keywords
- data
- delta
- time
- density
- lattice point
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 24
- 230000008033 biological extinction Effects 0.000 claims description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000003915 air pollution Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000013441 quality evaluation Methods 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
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W1/10—Devices for predicting weather conditions
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Chemical & Material Sciences (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Life Sciences & Earth Sciences (AREA)
- Databases & Information Systems (AREA)
- Dispersion Chemistry (AREA)
- Algebra (AREA)
- Atmospheric Sciences (AREA)
- Ecology (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Biodiversity & Conservation Biology (AREA)
- Environmental Sciences (AREA)
- Operations Research (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Image Processing (AREA)
Abstract
The invention discloses a method for inverting PM2.5 all-weather fine grid data based on surface high-density meteorological data, which utilizes high-density surface visibility and relative humidity meteorological observation data to construct ground PM2.5 all-weather time-by-time high-resolution grid data.
Description
Technical Field
The invention relates to a method for inverting PM2.5 all-weather fine grid data based on surface high-density meteorological data.
Background
The ground PM2.5 monitoring is limited to equipment and operation and maintenance cost, the number of sites is limited, the space distribution is uneven (the urban area is many and suburbs are few), how to acquire all-weather time-by-time fine grid PM2.5 data makes up the problem that the sites are rare or even blank in the suburbs in remote areas is always a hot scientific problem concerned by the academic community, and the method is also an urgent demand for air quality evaluation, PM2.5 human health influence research, air pollution treatment decision and the like; the general spatial interpolation algorithm is limited in that the error is large in a region with rare sites. The academic world provides a plurality of algorithms for fusing earth surface data and performing satellite remote sensing (AOD) inversion to obtain earth surface refined PM2.5, and the effects are good. The algorithm is limited in that satellite remote sensing AOD data are only 1-2 times per day, even existing high-resolution stationary satellites monitor time by time, the PM2.5 on the earth surface can not be monitored when clouds exist, and the cloud interference problem is particularly serious in southern areas. It is necessary to find a method for inverting ground PM2.5 time-by-time fine grid data, which does not depend on remote sensing data and has available precision quality.
The current automatic meteorological station is fast in construction and development, the ground visibility and relative humidity observation station point density which are closely related to PM2.5 are more than 10 times higher than that of the PM2.5 station point, and therefore a feasible basis is provided for establishing a method which only depends on high-density ground observation data.
Disclosure of Invention
According to the method, the ground PM2.5 all-weather time-by-time high-resolution grid data are constructed by utilizing meteorological observation data of high-density earth surface visibility and relative humidity, the problem that the conventional interpolation method has a large error in a rare area of a site is solved, and the problem of data loss caused by cloud interference based on a remote sensing method is also solved.
The technical scheme provided by the invention specifically comprises the following steps:
a method for inverting PM2.5 all-weather fine grid data based on surface high-density meteorological data comprises the following steps:
(1) calculating surface extinction delta from time-to-time or day-to-day visibility (V/km) data of a ground meteorological station in a target area, and performing spatial interpolation by adopting an inverse distance method to obtain delta lattice point data of a target resolution;
(2) performing spatial interpolation on the relative humidity (RH/%) data of the earth surface weather station time by time or day by adopting an inverse distance method to obtain RH lattice point data of a target resolution;
(3) according to PM2.5 observation stations in the target area, matching delta and RH data of corresponding positions, fitting an equation (1) one by one to obtain K and gamma parameters of each station;
RH in formula (1)0Taking 40% and delta being 3.912/V;
(4) adopting inverse distance spatial interpolation to obtain K and gamma lattice point data of target resolution for K and gamma parameters of PM2.5 observation stations in the region;
(5) the obtained K, gamma, delta and RH lattice point data are substituted for the formula (1) to obtain PM2.5 lattice point data with the target resolution.
As an improvement, the visibility data in the step (1) is collected from national stations and regional stations.
The invention has the following advantages: according to the method, the meteorological observation data of high-density earth surface visibility and relative humidity are utilized, the problem of data loss caused by cloud interference based on a remote sensing AOD method is solved, and all-weather time-by-time fine grid data of ground PM2.5 are constructed.
The beneficial effects include: (1) the method is beneficial to improving the capability of an air quality numerical prediction mode; (2) the PM2.5 evaluation precision level of the environmental air quality is improved, more accurate decision reference is provided for the government department to make an accurate atmospheric pollution control strategy, and the method has remarkable social and ecological benefits. (3) And key data support is provided for developing PM2.5 human health influence research.
Drawings
Fig. 1 is a schematic structural view of the present invention.
Detailed Description
With the combination of the attached drawings, the method for inverting PM2.5 all-weather fine grid data based on surface high-density meteorological data is characterized by comprising the following steps:
(1) calculating surface extinction delta from time-to-time or day-to-day visibility (V/km) data of a ground meteorological station in a target area, and performing spatial interpolation by adopting an inverse distance method to obtain delta lattice point data of a target resolution;
(2) performing spatial interpolation on the relative humidity (RH/%) data of the earth surface weather station time by time or day by adopting an inverse distance method to obtain RH lattice point data of a target resolution;
(3) according to PM2.5 observation stations in the target area, matching delta and RH data of corresponding positions, fitting an equation (1) one by one to obtain K and gamma parameters of each station;
RH in formula (1)0Taking 40% and delta being 3.912/V;
(4) adopting inverse distance spatial interpolation to obtain K and gamma lattice point data of target resolution for K and gamma parameters of PM2.5 observation stations in the region;
(5) the obtained K, gamma, delta and RH lattice point data are substituted for the formula (1) to obtain PM2.5 lattice point data with the target resolution.
The visibility data in the step (1) is collected from national stations and regional stations.
The present invention and its embodiments have been described above, and the description is not intended to be limiting, and the drawings are only one embodiment of the present invention, and the actual structure is not limited thereto. In summary, those skilled in the art should appreciate that they can readily use the disclosed conception and specific embodiments as a basis for designing or modifying other structures for carrying out the same purposes of the present invention without departing from the spirit and scope of the invention as defined by the appended claims.
Claims (2)
1. The method for inverting PM2.5 all-weather fine grid data based on the surface high-density meteorological data is characterized by comprising the following steps:
(1) calculating surface extinction delta from time-to-time or day-to-day visibility (V/km) data of a ground meteorological station in a target area, and performing spatial interpolation by adopting an inverse distance method to obtain delta lattice point data of a target resolution;
(2) performing spatial interpolation on the relative humidity (RH/%) data of the earth surface weather station time by time or day by adopting an inverse distance method to obtain RH lattice point data of a target resolution;
(3) according to PM2.5 observation stations in the target area, matching delta and RH data of corresponding positions, fitting an equation (1) one by one to obtain K and gamma parameters of each station;
RH in formula (1)0Taking 40% and delta being 3.912/V;
(4) adopting inverse distance spatial interpolation to obtain K and gamma lattice point data of target resolution for K and gamma parameters of PM2.5 observation stations in the region;
(5) the obtained K, gamma, delta and RH lattice point data are substituted for the formula (1) to obtain PM2.5 lattice point data with the target resolution.
2. The method for inverting PM2.5 all-weather fine grid data based on surface high-density meteorological data as claimed in claim 1, wherein the visibility data in the step (1) is collected from national stations and regional stations.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011362233.8A CN112525787A (en) | 2020-11-27 | 2020-11-27 | Method for inverting PM2.5 all-weather fine grid data based on surface high-density meteorological data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011362233.8A CN112525787A (en) | 2020-11-27 | 2020-11-27 | Method for inverting PM2.5 all-weather fine grid data based on surface high-density meteorological data |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112525787A true CN112525787A (en) | 2021-03-19 |
Family
ID=74994419
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011362233.8A Pending CN112525787A (en) | 2020-11-27 | 2020-11-27 | Method for inverting PM2.5 all-weather fine grid data based on surface high-density meteorological data |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112525787A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113111309A (en) * | 2021-03-25 | 2021-07-13 | 中国气象局广州热带海洋气象研究所(广东省气象科学研究所) | Method for evaluating influence of artificial and meteorological factors on concentration of atmospheric pollutants |
Citations (4)
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 |
CN103942439A (en) * | 2014-04-24 | 2014-07-23 | 中国科学院遥感与数字地球研究所 | Inhalable particle concentration estimating method based on meteorological observation data |
CN108426815A (en) * | 2018-04-20 | 2018-08-21 | 中国科学院遥感与数字地球研究所 | A kind of fine particle concentration of component evaluation method near the ground |
CN111912754A (en) * | 2020-07-23 | 2020-11-10 | 安徽省气象科学研究所 | Remote sensing inversion method for near-surface particulate matter concentration |
-
2020
- 2020-11-27 CN CN202011362233.8A patent/CN112525787A/en active Pending
Patent Citations (4)
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 |
CN103942439A (en) * | 2014-04-24 | 2014-07-23 | 中国科学院遥感与数字地球研究所 | Inhalable particle concentration estimating method based on meteorological observation data |
CN108426815A (en) * | 2018-04-20 | 2018-08-21 | 中国科学院遥感与数字地球研究所 | A kind of fine particle concentration of component evaluation method near the ground |
CN111912754A (en) * | 2020-07-23 | 2020-11-10 | 安徽省气象科学研究所 | Remote sensing inversion method for near-surface particulate matter concentration |
Non-Patent Citations (2)
Title |
---|
李婷苑等: "卫星遥感AOD反演地面细颗粒物浓度方法与效果", 《中国环境科学》 * |
许娈等: "基于逐时降水站点资料空间插值方法对比研究", 《气象与环境学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113111309A (en) * | 2021-03-25 | 2021-07-13 | 中国气象局广州热带海洋气象研究所(广东省气象科学研究所) | Method for evaluating influence of artificial and meteorological factors on concentration of atmospheric pollutants |
CN113111309B (en) * | 2021-03-25 | 2022-10-14 | 中国气象局广州热带海洋气象研究所(广东省气象科学研究所) | Method for evaluating influence of artificial and meteorological factors on concentration of atmospheric pollutants |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110363327B (en) | ConvLSTM and 3D-CNN-based short rainfall prediction method | |
CN111323352B (en) | Regional PM2.5 remote sensing inversion model fusing fine particulate matter concentration data | |
CN112070286A (en) | Rainfall forecast early warning system for complex terrain watershed | |
Li et al. | Urban morphology in China: dataset development and spatial pattern characterization | |
CN107958312B (en) | Electric transmission line galloping prediction method and system based on inversion algorithm and storage medium | |
CN110058329A (en) | A kind of meteorological element smart grid forecast model products space NO emissions reduction method | |
CN106772696A (en) | For the total method for dodging translation proxy humidity of assimilation of Severe Convective Weather Forecasting | |
CN106971396A (en) | Ice sheet freeze thawing detection method based on super-pixel | |
CN108959347B (en) | Aurora observation candidate region determination method based on terrain and night light data | |
Li et al. | Changes in the glacier extent and surface elevation in Xiongcaigangri region, Southern Karakoram Mountains, China | |
CN112525787A (en) | Method for inverting PM2.5 all-weather fine grid data based on surface high-density meteorological data | |
Sgoff et al. | Assimilation of crowd‐sourced surface observations over Germany in a regional weather prediction system | |
Han et al. | Spatiotemporal changes in the urban heat island intensity of distinct local climate zones: Case study of Zhongshan District, Dalian, China | |
CN110806606B (en) | Electric power weather refined forecast early warning system and method | |
CN115511192B (en) | Precipitation prediction method and system based on lightning data assimilation | |
Liu et al. | Building Chinese wind data for Wind Erosion Prediction System using surrogate US data | |
Fujibe | Evaluation of background and urban warming trends based on centennial temperature data in Japan | |
CN112285808A (en) | Method for reducing scale of APHRODITE precipitation data | |
Lussana et al. | Spatial interpolation of two‐metre temperature over Norway based on the combination of numerical weather prediction ensembles and in situ observations | |
CN116258065A (en) | Cloud high space-time resolution near-ground air temperature estimation method | |
CN107688712B (en) | A kind of temperature NO emissions reduction method based on DEM and NDVI | |
KR100975202B1 (en) | Method of estimating distribution of wind velocities | |
Rosli et al. | Sustainable urban forestry potential based quantitative and qualitative measurement using geospatial technique | |
Yang et al. | Mapping wind erosion hazard in Australia using MODIS-derived ground cover, soil moisture and climate data | |
CN113568008A (en) | Sunflower-8 satellite-based minute-level precipitation real-time inversion estimation method |
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 |
Application publication date: 20210319 |
|
RJ01 | Rejection of invention patent application after publication |