CN109902913A - A kind of identification of atmosphere pollution superposed risk area and evaluation method - Google Patents
A kind of identification of atmosphere pollution superposed risk area and evaluation method Download PDFInfo
- Publication number
- CN109902913A CN109902913A CN201910024455.XA CN201910024455A CN109902913A CN 109902913 A CN109902913 A CN 109902913A CN 201910024455 A CN201910024455 A CN 201910024455A CN 109902913 A CN109902913 A CN 109902913A
- Authority
- CN
- China
- Prior art keywords
- pollution
- atmosphere
- atmosphere pollution
- risk
- region
- 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
- 238000011156 evaluation Methods 0.000 title claims abstract description 10
- 239000003344 environmental pollutant Substances 0.000 claims abstract description 22
- 231100000719 pollutant Toxicity 0.000 claims abstract description 22
- 230000005540 biological transmission Effects 0.000 claims abstract description 13
- 239000011159 matrix material Substances 0.000 claims abstract description 10
- 230000009467 reduction Effects 0.000 claims abstract description 3
- 238000000034 method Methods 0.000 claims description 6
- 239000000809 air pollutant Substances 0.000 claims description 4
- 231100001243 air pollutant Toxicity 0.000 claims description 4
- 238000011160 research Methods 0.000 claims description 3
- 230000008676 import Effects 0.000 claims description 2
- 238000012067 mathematical method Methods 0.000 claims description 2
- 238000004088 simulation Methods 0.000 abstract description 5
- 238000009792 diffusion process Methods 0.000 abstract description 3
- 239000000356 contaminant Substances 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 7
- 238000003915 air pollution Methods 0.000 description 2
- 238000007599 discharging Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention relates to a kind of identification of atmosphere pollution superposed risk area and evaluation methods, the Polluted area influenced by multiple sources of atmospheric pollution is considered as pollution overlap-add region, risk zones identification and grading evaluation are carried out to the severity of overlying pollution suffered by it, specifically include that the diffusion process for combining certain Regional And Multi-source atmosphere pollution under WRF and CALPUFF pattern simulation DIFFERENT METEOROLOGICAL CONDITIONS, certain pollution sources is calculated to the contribution rate of other region of acceptance atmosphere pollution, the contaminant transmission matrix of each pollution sources is obtained, and then analyzes the interactional rule of atmosphere overlying pollution.The group of cities pollutant spatial distribution state simulated by analyzing dispersal pattern identifies atmosphere pollution superposed risk region and carries out grading evaluation according to certain index, provides guidance for each department emission reduction work.
Description
Technical field
The invention belongs to Air Pollution Control fields, the in particular to a kind of identification of atmosphere pollution superposed risk area and evaluation side
Method.
Background technique
Atmosphere pollution is influenced by meteorological condition, the migration between easy generation area, in the atmosphere pollution of certain discharge
Substance not only has an impact discharge source region, but will be to the atmospheric environment in region most of around influencing.Now, for more
The case where source air pollution emission, is more common, and multiple pollution sources, which discharge pollutants, to spread in an atmosphere and influence each other.I
Great shift is occurring for state's regionality atmosphere pollution character, and high-intensitive, multiple pollutant sources discharge atmosphere pollution is in its discharge
It is difficult to self-cleaning inside region, and then transmits, convert between the group of cities that high density is gathered, leads to regional, compound atmosphere
Pollution problem is increasingly prominent.Numerous studies discovery, atmosphere pollution show certain region homogeneity feature.High-intensitive source row
Putting no doubt is the main reason for causing Regional Atmospheric Pollution, and multiple rows of diffusion for putting source pollutants, can be to each interregional generation phase
It mutually influences, further initiation area atmosphere combined pollution aggravates the exposure of crowd.Therefore, it is necessary to study multi-source atmosphere
The combined pollution effect generated to region is polluted, the pollutant effect degree except itself suffered discharge of evaluation region plans as a whole area
Domain discharge, mitigates atmosphere pollution.
Summary of the invention
It is an object of the invention to overcome the deficiencies in the prior art, provide a kind of identification of atmosphere pollution superposed risk area and comment
Valence method, is simulated by software, analyzes the mechanism of transmission of the pollutant under DIFFERENT METEOROLOGICAL CONDITIONS, quantization areas is by overlying pollution
Degree.
To achieve the above object, the invention provides the following technical scheme: a kind of atmosphere pollution superposed risk area identifies and comments
Valence method, includes the following steps:
1) quantify various Air Pollutants status transmissions under DIFFERENT METEOROLOGICAL CONDITIONS, construct the atmosphere pollution between pollution sources
Object transmission matrix；
2) other pollution sources are calculated to the pollution condition contribution rate in certain region, and refine degree of risk on this basis, are drawn
Determine risk class；
3) it is based on above-mentioned group of cities atmosphere pollution transportation law, atmosphere overlying pollution risk area risk class index, is calculated
Each pollution sources emission reduction scale degree.
The present invention is using PM2.5 as polluter, and specific step is as follows for this method:
1) according to institute's survey region, CALPUFF modeling regional scope is set；
2) disposal of pollutants source is converted and is extracted with ArcGIS, by pollution sources point data and discharge amount and including
Chimney height and the data of internal diameter, smog discharge temperature and rate import CALPUFF software；
3) meteorological condition in WRF modeling region is used, and analogue data is imported into CALPUFF software；
4) according to model output data, atmosphere pollution transmission matrix is constructed with mathematical method；
5) according to other regions in atmosphere pollution transmission matrix calculating research range to some region of pollutant contribution rate；
6) risk class is divided according to the calculated contribution rate of previous step, completes identification and evaluation to atmosphere pollution superposed risk area.
The present invention has the advantage that
1. utilizing CALPUFF and WRF modeling area atmosphere pollution transmission process, the fast precision of speed is quasi-.It can be accurate
This area's atmosphere pollution space-time characteristic is disclosed, atmosphere overlying pollution risk area Distribution Pattern is disclosed.And meteorological condition can be changed etc.
The key factor of pollutant diffusion is influenced, quantitative analysis pollutant is in interregional transportation law, to predict atmosphere combined pollution
Situation provides guidance.
2. can divide from the angle implemented zone defence joint control, optimize regional space pattern on the basis of theoretical result
With each region energy-saving and emission-reduction task, Regional Atmospheric Pollution prevention and treatment and regulating strategy are constructed.
Detailed description of the invention
Fig. 1 is that schematic diagram is arranged in survey region grid；
Fig. 2 is pollution sources processing setting schematic diagram；
Fig. 3 is WRF meteorology Scene Simulation schematic diagram；
Fig. 4 is simulated scenario plan of establishment schematic diagram；
Fig. 5 is emission source Scene Simulation numerical value schematic diagram；
Fig. 6 is the atmosphere pollution transmission matrix schematic diagram under specified context；
Fig. 7 is delimitation and the grade scale of overlying pollution risk area.
Specific embodiment
Below by the drawings and specific embodiments, the present invention is further illustrated.The embodiment of the present invention is in order to more
So that those skilled in the art is more fully understood the present invention well, any limitation is not made to the present invention.Based in the present invention
Embodiment, every other embodiment obtained by those of ordinary skill in the art without making creative efforts,
It shall fall within the protection scope of the present invention.
First, research area's grid setting.As shown in Figure 1, the simulated domain range of CALPUFF covers Jing-jin-ji region group of cities.
General Lambert conformality circular cone (LCC) is taken to project, geodetic level WGS-84, central point longitude and latitude (116.21342 ° of E,
39.584055 ° of N), standard parallel is respectively 35 ° of N and 45 ° of N.Simulated domain east-west direction distance 576km, totally 96 lattice points, southern
The north is to distance 738km, totally 123 lattice points, grid spacing 6km.Vertical direction is arranged 11 layers altogether, and each layer heights of roofs is respectively
0m, 20m, 40m, 80m, 160m, 320m, 640m, 1200m, 2000m, 3000m and 4000m.
Followed by pollution sources are handled.Briefly conversion is done using emissions data of the ArcGIS to gridding and is extracted, with each net
As point source, the discharge amount of the point source is the total release of the grid at lattice center.Obtain 362 emission sources, spatial distribution state
As shown in Figure 2.Emission source substrate height above sea level is provided by landform altitude file, remaining emission source input parameter for example chimney height with
Internal diameter, smog discharge temperature and rate etc., referring to relevant criterion (such as " discharge standard of air pollutants GB16297-
1996 ", " emission standard of air pollutants for boilers GB13271-2001 " etc.) and determine.
Then, WRF meteorology Scene Simulation.As shown in figure 3, special in order to more accurately obtain the meteorology of Jing-jin-ji region group of cities
Sign, WRF mode use double nested scheme, and outer layer grid is domain01, and internal layer grid is domain02, and domain01 is
Domain02 provides boundary condition.Domain02 cover CALPUFF simulated domain, central point longitude and latitude be (116.21342 ° of E,
39.584055° N).The spatial resolution of inside and outside grid is respectively 6km and 18km, and horizontal grid number is respectively 96 × 123 and 60
×63.The parameters such as such as meteorological condition, wind direction, emission source and discharge species can be then modified, so as to different meteorologies
Scene is simulated, as shown in Figure 4.
4th step is intercity atmosphere pollution transmission matrix building.Fig. 5 is emission source Scene Simulation numerical value schematic diagram.
The receptor of 13 cities as pollutants is studied using in region, it is assumed that pollutant concentration of the source in the city k to the city i
Contribution rate is Rik, then
Wherein, CikUnder scene only to consider the city k source emission, the mean concentration of the city i pollutant；AiIt is active for group of cities
When discharging simultaneously, the mean concentration of the city i pollutant.With RikAs matrix element, group of cities atmosphere pollution transmission square can be constructed
Battle array, as shown in fig. 6, clearly reflecting the source of pollutant and whereabouts in each city.
It is finally the identification and classification of overlying pollution risk area.If the pollutant concentration in a certain area is by the dirt in multiple cities
Dye source dominates, then this area can suffer from overlying pollution；If the pollutant in a certain area is almost entirely from city
Pollution sources, then this area does not have the risk of overlying pollution.The receptor of each grid as pollutant is studied using in region, false
If the source in the city k is r to the pollutant concentration contribution rate of j-th of gridjk, then
Wherein, SjkUnder scene only to consider the city k source emission, the total quantity control on emission of j-th of grid；GjFor group of cities
Institute is active while when discharging, the total quantity control on emission of j-th of grid.
If there are two or more than two city, to the contribution rate r of the pollutant concentration of a certain grid in 20% to 80%,
Then the grid faces the risk of overlying pollution.In order to refine degree of risk, according to rjkThe distribution of value delimit risk class,
As shown in Fig. 7.For example, if a certain city is more than 80% to the contribution rate of grid j, which is not present superposed risk；
If certain two or three city is to the contribution rate of grid j 30%~40% or 60%~Between 70%, then the grid faces superposition
Pollute second level risk.
Although the preferred embodiment of the present invention is described above in conjunction with attached drawing, the invention is not limited to upper
The specific embodiment stated, the above mentioned embodiment is only schematical, be not it is restrictive, this field it is common
Technical staff under the inspiration of the present invention, without breaking away from the scope protected by the purposes and claims of the present invention, may be used also
By make it is many in the form of specific transformation, within these are all belonged to the scope of protection of the present invention.
Claims (2)
1. a kind of identification of atmosphere pollution superposed risk area and evaluation method, which comprises the steps of:
1) quantify various Air Pollutants status transmissions under DIFFERENT METEOROLOGICAL CONDITIONS, the atmosphere pollution constructed between pollution sources passes
Defeated matrix；
2) other pollution sources are calculated to the pollution condition contribution rate in certain region, and refine degree of risk on this basis, delimit wind
Dangerous grade；
3) it is based on above-mentioned group of cities atmosphere pollution transportation law, atmosphere overlying pollution risk area risk class index, calculates each dirt
Dye source emission reduction scale degree.
2. a kind of atmosphere pollution superposed risk area identification according to claim 1 and evaluation method, which is characterized in that with
PM2.5 is polluter, and specific step is as follows for this method:
1) according to institute's survey region, CALPUFF modeling regional scope is set；
2) disposal of pollutants source is converted and is extracted with ArcGIS, by pollution sources point data and discharge amount and including chimney
The data of height and internal diameter, smog discharge temperature and rate import CALPUFF software；
3) meteorological condition in WRF modeling region is used, and analogue data is imported into CALPUFF software；
4) according to model output data, atmosphere pollution transmission matrix is constructed with mathematical method；
5) according to other regions in atmosphere pollution transmission matrix calculating research range to some region of pollutant contribution rate；
6) risk class is divided according to the calculated contribution rate of previous step, complete the identification to atmosphere pollution superposed risk area and commented
Valence.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910024455.XA CN109902913A (en) | 2019-01-10 | 2019-01-10 | A kind of identification of atmosphere pollution superposed risk area and evaluation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910024455.XA CN109902913A (en) | 2019-01-10 | 2019-01-10 | A kind of identification of atmosphere pollution superposed risk area and evaluation method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109902913A true CN109902913A (en) | 2019-06-18 |
Family
ID=66943588
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910024455.XA Pending CN109902913A (en) | 2019-01-10 | 2019-01-10 | A kind of identification of atmosphere pollution superposed risk area and evaluation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109902913A (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110428104A (en) * | 2019-08-01 | 2019-11-08 | 软通动力信息技术有限公司 | A kind of genes' contamination ratio determines method, apparatus, electronic equipment and storage medium |
| CN110674583A (en) * | 2019-09-29 | 2020-01-10 | 天津大学 | Method for constructing inter-city atmospheric pollution transmission matrix |
| CN110727984A (en) * | 2019-09-29 | 2020-01-24 | 天津大学 | Method for researching mechanism between land utilization change and urban thermal environment |
| CN110849837A (en) * | 2019-11-26 | 2020-02-28 | 南昌工程学院 | Atmospheric pollution superposition risk area identification and evaluation method |
| CN111178786A (en) * | 2020-01-08 | 2020-05-19 | 中科宇图科技股份有限公司 | Emission source position determining method and system for guaranteeing regional air quality |
| CN111257250A (en) * | 2020-01-21 | 2020-06-09 | 河北省农林科学院农业资源环境研究所 | Identification method and differentiation processing method for high-risk area of land agricultural product producing area |
| CN114965883A (en) * | 2022-05-17 | 2022-08-30 | 南京信息职业技术学院 | Method for judging main atmospheric pollutants and main atmospheric pollution sources in area |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102628852A (en) * | 2012-03-13 | 2012-08-08 | 北京工业大学 | Atmospheric pollution source grading method based on pollutant source identification technology |
| CN106407714A (en) * | 2016-10-14 | 2017-02-15 | 珠海富鸿科技有限公司 | Air pollution assessment method and device based on CALPUFF system |
| CN107194139A (en) * | 2016-03-14 | 2017-09-22 | 日电（中国）有限公司 | Source of atmospheric pollution stage division and computing device |
| US20180082199A1 (en) * | 2016-09-21 | 2018-03-22 | International Business Machines Corporation | System, method and computer program product for pollution source attribution |
| CN107917987A (en) * | 2017-11-09 | 2018-04-17 | 北京伟瑞迪科技有限公司 | Urban air pollutant source tracing analysis method |
| CN107941994A (en) * | 2017-11-09 | 2018-04-20 | 北京伟瑞迪科技有限公司 | It is a kind of towards industrial park pollutant Source Tracing method |
-
2019
- 2019-01-10 CN CN201910024455.XA patent/CN109902913A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102628852A (en) * | 2012-03-13 | 2012-08-08 | 北京工业大学 | Atmospheric pollution source grading method based on pollutant source identification technology |
| CN107194139A (en) * | 2016-03-14 | 2017-09-22 | 日电（中国）有限公司 | Source of atmospheric pollution stage division and computing device |
| US20180082199A1 (en) * | 2016-09-21 | 2018-03-22 | International Business Machines Corporation | System, method and computer program product for pollution source attribution |
| CN106407714A (en) * | 2016-10-14 | 2017-02-15 | 珠海富鸿科技有限公司 | Air pollution assessment method and device based on CALPUFF system |
| CN107917987A (en) * | 2017-11-09 | 2018-04-17 | 北京伟瑞迪科技有限公司 | Urban air pollutant source tracing analysis method |
| CN107941994A (en) * | 2017-11-09 | 2018-04-20 | 北京伟瑞迪科技有限公司 | It is a kind of towards industrial park pollutant Source Tracing method |
Non-Patent Citations (2)
| Title |
|---|
| 孟伟 等: "北京及周边地区大气污染数值模拟研究", 《环境科学研究》 * |
| 李文慧 等: "关中二氧化硫污染特征及城市间相互影响分析", 《环境科学与技术》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110428104A (en) * | 2019-08-01 | 2019-11-08 | 软通动力信息技术有限公司 | A kind of genes' contamination ratio determines method, apparatus, electronic equipment and storage medium |
| CN110674583A (en) * | 2019-09-29 | 2020-01-10 | 天津大学 | Method for constructing inter-city atmospheric pollution transmission matrix |
| CN110727984A (en) * | 2019-09-29 | 2020-01-24 | 天津大学 | Method for researching mechanism between land utilization change and urban thermal environment |
| CN110849837A (en) * | 2019-11-26 | 2020-02-28 | 南昌工程学院 | Atmospheric pollution superposition risk area identification and evaluation method |
| CN111178786A (en) * | 2020-01-08 | 2020-05-19 | 中科宇图科技股份有限公司 | Emission source position determining method and system for guaranteeing regional air quality |
| CN111178786B (en) * | 2020-01-08 | 2023-04-25 | 中科宇图科技股份有限公司 | Emission source position determining method and system for guaranteeing regional air quality |
| CN111257250A (en) * | 2020-01-21 | 2020-06-09 | 河北省农林科学院农业资源环境研究所 | Identification method and differentiation processing method for high-risk area of land agricultural product producing area |
| CN114965883A (en) * | 2022-05-17 | 2022-08-30 | 南京信息职业技术学院 | Method for judging main atmospheric pollutants and main atmospheric pollution sources in area |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109902913A (en) | A kind of identification of atmosphere pollution superposed risk area and evaluation method | |
| Yang et al. | Air pollution dispersal in high density urban areas: Research on the triadic relation of wind, air pollution, and urban form | |
| CN106649960B (en) | Three-dimensional iterative calculation method for atmospheric multi-pollutant environment capacity | |
| CN104865353B (en) | Industrial park atmosphere pollution collecting method based on unmanned plane | |
| CN102628852B (en) | Atmospheric pollution source grading method based on pollutant source identification technology | |
| Shafiev et al. | Studies of the influence of vegetation cover on the process of transfer and diffusion of harmful substances in the atmosphere | |
| Wang et al. | Strategies for sustainable urban development and morphological optimization of street canyons: Measurement and simulation of PM2. 5 at different points and heights | |
| CN115420854A (en) | Atmospheric pollutant tracing method based on forward and backward model combination | |
| Kwok et al. | To what extent can urban ventilation features cool a compact built-up environment during a prolonged heatwave? A mesoscale numerical modelling study for Hong Kong | |
| KR100661595B1 (en) | Decision Support System for Air pollution source planning: ApspDSS | |
| Peng et al. | An experimental work to investigate the capabilities of plants to remove particulate matters in an enclosed greenhouse | |
| Jiang et al. | Revealing the spatiotemporal characteristics and drivers of the block-scale thermal environment near a large river: Evidences from Shanghai, China | |
| Zhang et al. | Natural and human factors influencing urban particulate matter concentrations in central heating areas with long-term wearable monitoring devices | |
| Zhou et al. | The impact of urban morphology on multiple ecological effects: Coupling relationships and collaborative optimization strategies | |
| CN110674583A (en) | Method for constructing inter-city atmospheric pollution transmission matrix | |
| Harkat et al. | The impact of industrial air pollution on the urban environment of setif: Modeling and mapping of total suspended particles | |
| Hsieh et al. | Improving the local wind environment through urban design strategies in an urban renewal process to mitigate urban heat island effects | |
| Wang et al. | Influence of grass lawns on the summer thermal environment and microclimate of heritage sites: a case study of Fuling mausoleum, China | |
| CN105844073A (en) | Source/sink landscape recognition method and apparatus of non-point source pollution | |
| Kang et al. | The Influence of Meteorological Conditions and Complex Topography on Ozone Concentrations in a Valley Area near Coastal Metropolitan Cities. | |
| Musco et al. | Mitigation of and adaptation to UHI phenomena: The Padua case study | |
| CN103617452A (en) | Photometry network layout method in large-scale photovoltaic base area | |
| Biao et al. | The dynamic coupling model and its application of urbanization and eco-environment in Hexi Corridor | |
| Chiang et al. | A Study on the cooling effects of the park dispersions applied in PUD | |
| EA | OPTIMUM FORM AND PLACEMENT OF URBAN BLOCKS TO MAXIMIZE THE USE OF SOLAR ENERG–A CASE STUDY |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190618 |
|
| WD01 | Invention patent application deemed withdrawn after publication |