CN110674583A - Method for constructing inter-city atmospheric pollution transmission matrix - Google Patents
Method for constructing inter-city atmospheric pollution transmission matrix Download PDFInfo
- Publication number
- CN110674583A CN110674583A CN201910930648.1A CN201910930648A CN110674583A CN 110674583 A CN110674583 A CN 110674583A CN 201910930648 A CN201910930648 A CN 201910930648A CN 110674583 A CN110674583 A CN 110674583A
- Authority
- CN
- China
- Prior art keywords
- atmospheric pollution
- city
- pollution
- transmission matrix
- inter
- 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
- 230000005540 biological transmission Effects 0.000 title claims abstract description 27
- 239000011159 matrix material Substances 0.000 title claims abstract description 21
- 238000000034 method Methods 0.000 title claims description 5
- 239000003344 environmental pollutant Substances 0.000 claims abstract description 19
- 231100000719 pollutant Toxicity 0.000 claims abstract description 19
- 238000004088 simulation Methods 0.000 claims abstract description 12
- 239000000779 smoke Substances 0.000 claims description 5
- 238000010276 construction Methods 0.000 claims description 4
- 102000016979 Other receptors Human genes 0.000 claims description 2
- 238000004519 manufacturing process Methods 0.000 claims description 2
- 238000012067 mathematical method Methods 0.000 claims description 2
- 238000003915 air pollution Methods 0.000 abstract description 4
- 238000009792 diffusion process Methods 0.000 abstract description 4
- 238000011160 research Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 5
- 239000000809 air pollutant Substances 0.000 description 1
- 231100001243 air pollutant Toxicity 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 102000005962 receptors Human genes 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
Landscapes
- Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Strategic Management (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention researches the mutual influence of air pollution between cities under two conditions of north wind and south wind, and mainly comprises the following steps: a WRF mode and a CALPUFF mode are combined to simulate the diffusion process of atmospheric pollution of different cities in a certain area under two meteorological conditions of a northern wind and a southern wind, the contribution rate of pollution sources of each city to the atmospheric pollution of other receptors is calculated, a pollutant transmission matrix between cities is obtained, and then the rule of mutual influence of the atmospheric pollution between the cities is analyzed. And (3) establishing a transmission matrix of the inter-city atmospheric pollution by analyzing the inter-city pollutant transmission condition obtained by diffusion mode simulation, and analyzing the rule of the mutual influence of the inter-city atmospheric pollution.
Description
Technical Field
The invention belongs to the field of atmospheric pollution control, and particularly relates to a construction method of an inter-city atmospheric pollution transmission matrix.
Background
Due to the change of meteorological conditions, the area range of the atmospheric pollution is changed, and the atmospheric pollutants discharged at a certain position not only affect the area where the emission source is located, but also affect the atmospheric environment of most areas around the emission source. As the condition of multi-source atmospheric pollution emission is more common, the problem of urban atmospheric pollution is also concerned widely, and pollutants emitted by a plurality of pollution sources are diffused in the atmosphere and mutually influenced. The regional air pollution characteristics of China are undergoing important transformation, the high-intensity and multi-emission-source air pollutants are difficult to self-clean in the emission area, and then are transmitted and transformed among high-density aggregated urban groups, so that the regional and composite air pollution problem is increasingly prominent. Due to the diffusion of pollutants from multiple sources, the areas will interact with each other. Therefore, it is necessary to study the influence of air pollution among cities in a city group and evaluate the degree of influence of pollutants other than self-discharge in the city.
Disclosure of Invention
The invention aims to provide a construction method of an inter-city atmospheric pollution transmission matrix. By software simulation, the transmission mechanism of pollutants under different meteorological conditions is analyzed, and the transmission rule of atmospheric pollution among cities is analyzed.
In order to achieve the purpose, the invention provides the following technical scheme: a construction method of an inter-city atmospheric pollution transmission matrix comprises the following steps:
1) counting the emission conditions of atmospheric pollutants in each city;
2) under different meteorological conditions, calculating the contribution rate of pollution sources of each city to the atmospheric pollution of other receptor cities, and manufacturing an inter-city atmospheric pollution transmission matrix according to the contribution rate;
3) and analyzing the inter-city atmospheric pollution transmission rule based on the inter-city atmospheric pollution transmission matrix.
The invention takes PM2.5 as a pollutant, and the method comprises the following specific steps:
1) setting the range of a CALPUFF model simulation area according to the researched area;
2) converting and extracting a pollution emission source by using ArcGIS, and introducing pollution source data and data comprising chimney height, chimney inner diameter, smoke emission temperature and smoke emission rate into CALPUFF software;
3) simulating meteorological conditions of a region by using a WRF model, and importing simulation data into CALPUFF software;
4) constructing an atmospheric pollution transmission matrix by a mathematical method according to the model output data;
5) and analyzing the rule of mutual influence of the atmospheric pollution among cities according to the transmission matrix calculated in the last step.
The invention has the following advantages:
1. and the CALPUFF and WRF models are used for simulating the regional atmospheric pollution transmission process, so that the speed is high and the precision is accurate. The atmospheric pollution time-space characteristics of the region can be accurately revealed, key factors influencing pollutant diffusion such as meteorological conditions can be changed, and the transmission rule of pollutants among regions can be quantitatively analyzed.
2. On the basis of theoretical results, a regional atmosphere pollution control and regulation strategy can be constructed from the perspective of implementing regional joint defense joint control and optimizing urban spatial patterns.
Drawings
FIG. 1 is a schematic view of a study area grid setup;
FIG. 2 is a schematic view of a pollution source treatment setup;
FIG. 3 is a diagram of a WRF weather scenario simulation;
FIG. 4 is a schematic diagram of a simulation scenario setup scenario;
FIG. 5 is a diagram of emission source scenario simulation values;
fig. 6 is a schematic diagram of an atmospheric pollution transmission matrix in a specific situation.
Detailed Description
The invention is further illustrated by the following figures and specific examples. The examples are given to better understand the present invention for those skilled in the art, and do not limit the present invention in any way. All other embodiments, which can be derived by a person skilled in the art from the examples given herein without making any creative effort, shall fall within the protection scope of the present invention.
First, the study area grid is set up. As shown in fig. 1, CALPUFF covers the range of the simulation area in shandong province. The universal Lambert orthoconic (LCC) projection is adopted, the geodetic datum plane is WGS-84, the longitude and latitude of the center point are 118.12704 degrees E and 36.340947 degrees N, and the standard latitude lines are 35 degrees N and 45 degrees N respectively. The distance of the simulated area in the east-west direction is 740km, the distance of the simulated area in the south-north direction is 148 km, the distance of the simulated area in the south-north direction is 455km, the distance of the simulated area in the south-north direction is 91 km, and the grid distance is 5 km. . In the vertical direction, 11 layers are arranged, and the top heights of the layers are respectively 0m, 20m, 40m, 80m, 160m, 320m, 640m, 1200m, 2000m, 3000m and 4000 m.
And secondly, pollution source treatment. And briefly converting and extracting the gridded emission data by utilizing ArcGIS, taking the center of each grid as a point source, and taking the emission of the point source as the total emission of the grid. 247 emission sources are obtained, and the spatial distribution is shown in fig. 2. The altitude of the emission source is provided by a terrain elevation file, and other input parameters of the emission source, such as the height and the inner diameter of a chimney, the temperature and the speed of smoke emission, and the like, are determined according to relevant standards (such as integrated emission standard of atmospheric pollutants GB16297-1996, emission standard of atmospheric pollutants GB13271-2001 of a boiler, and the like).
Then, WRF weather scene simulation. As shown in fig. 3, in order to more accurately obtain meteorological features of the Shandong province, the WRF mode adopts a double nesting scheme, wherein the outer layer grid is domain01, the inner layer grid is domain02, and domain01 provides boundary conditions for domain 02. domain02 covers the CALPUFF simulation area with a center point of latitude and longitude (118.12704 ° E, 36.340947 ° N). The spatial resolution of the inner and outer grids was 5km and 30km, respectively, and the number of horizontal grids was 91 × 148 and 100 × 127, respectively. Parameters such as meteorological conditions, wind direction, emissions source, and emissions type may then be altered to simulate different meteorological scenarios, as shown in FIG. 4.
And fourthly, constructing an inter-city atmospheric pollution transmission matrix. Fig. 5 is a diagram of emission source scenario simulation values.
The research is carried out by taking 17 cities in the area as the receptors of the pollutants, and the concentration contribution rate of the sources of k city to i city is assumed to be RikThen, then
Wherein, CikThe average concentration of pollutants in the city i under the condition of only considering the emission of the source in the city k; a. theiThe average concentration of pollutants in city when all sources in the city group are discharged simultaneously. With RikAs matrix elements, one can constructThe urban atmospheric pollution transmission matrix, as shown in fig. 6, clearly reflects the source and destination of pollutants in each city.
Although the preferred embodiments of the present invention have been described above with reference to the accompanying drawings, the present invention is not limited to the above-described embodiments, which are merely illustrative and not restrictive, and those skilled in the art can make various changes and modifications within the spirit and scope of the present invention without departing from the spirit and scope of the appended claims.
Claims (2)
1. The construction method of the urban atmospheric pollution transmission matrix is characterized by comprising the following steps:
1) counting the emission conditions of atmospheric pollutants in each city;
2) under different meteorological conditions, calculating the contribution rate of pollution sources of each city to the atmospheric pollution of other receptor cities, and manufacturing an inter-city atmospheric pollution transmission matrix according to the contribution rate;
3) and analyzing the inter-city atmospheric pollution transmission rule based on the inter-city atmospheric pollution transmission matrix.
2. The method for constructing the urban atmospheric pollution transmission matrix according to claim 1, wherein PM2.5 is used as a pollutant, and the method comprises the following specific steps:
1) setting the range of a CALPUFF model simulation area according to the researched area;
2) converting and extracting a pollution emission source by using ArcGIS, and introducing pollution source data and data comprising chimney height, chimney inner diameter, smoke emission temperature and smoke emission rate into CALPUFF software;
3) simulating meteorological conditions of a region by using a WRF model, and importing simulation data into CALPUFF software;
4) constructing an atmospheric pollution transmission matrix by a mathematical method according to the model output data;
5) and analyzing the rule of mutual influence of the atmospheric pollution among cities according to the transmission matrix calculated in the last step.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910930648.1A CN110674583A (en) | 2019-09-29 | 2019-09-29 | Method for constructing inter-city atmospheric pollution transmission matrix |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910930648.1A CN110674583A (en) | 2019-09-29 | 2019-09-29 | Method for constructing inter-city atmospheric pollution transmission matrix |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110674583A true CN110674583A (en) | 2020-01-10 |
Family
ID=69080053
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910930648.1A Pending CN110674583A (en) | 2019-09-29 | 2019-09-29 | Method for constructing inter-city atmospheric pollution transmission matrix |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110674583A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113191096A (en) * | 2021-04-13 | 2021-07-30 | 中南建筑设计院股份有限公司 | WRF and XLow coupling-based multi-fineness fusion pollutant diffusion analysis method |
| CN113674135A (en) * | 2021-08-03 | 2021-11-19 | 北京三易思创科技有限公司 | Calculation method for realizing CALPUFF high performance based on workflow |
| CN113838539A (en) * | 2021-09-24 | 2021-12-24 | 国网福建省电力有限公司 | Detection and analysis system and analysis method for fire residues of electrical equipment of transformer substation |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2797041A1 (en) * | 2013-04-26 | 2014-10-29 | Università Del Salento - Dipartimento Di Ingegneria Dell'Innovazione | Method for estimating pollutant concentration of urban environments using fluid-dynamics clustering |
| CN107563562A (en) * | 2017-09-08 | 2018-01-09 | 新奥泛能网络科技股份有限公司 | The appraisal procedure and device of gridding urban atmospheric pollution thing discharge |
| CN109902913A (en) * | 2019-01-10 | 2019-06-18 | 天津大学 | A kind of identification of atmosphere pollution superposed risk area and evaluation method |
| CN109932988A (en) * | 2019-03-27 | 2019-06-25 | 四川瞭望工业自动化控制技术有限公司 | A kind of city raised dust contamination forecasting system and method |
-
2019
- 2019-09-29 CN CN201910930648.1A patent/CN110674583A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2797041A1 (en) * | 2013-04-26 | 2014-10-29 | Università Del Salento - Dipartimento Di Ingegneria Dell'Innovazione | Method for estimating pollutant concentration of urban environments using fluid-dynamics clustering |
| CN107563562A (en) * | 2017-09-08 | 2018-01-09 | 新奥泛能网络科技股份有限公司 | The appraisal procedure and device of gridding urban atmospheric pollution thing discharge |
| CN109902913A (en) * | 2019-01-10 | 2019-06-18 | 天津大学 | A kind of identification of atmosphere pollution superposed risk area and evaluation method |
| CN109932988A (en) * | 2019-03-27 | 2019-06-25 | 四川瞭望工业自动化控制技术有限公司 | A kind of city raised dust contamination forecasting system and method |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113191096A (en) * | 2021-04-13 | 2021-07-30 | 中南建筑设计院股份有限公司 | WRF and XLow coupling-based multi-fineness fusion pollutant diffusion analysis method |
| CN113191096B (en) * | 2021-04-13 | 2022-04-29 | 中南建筑设计院股份有限公司 | WRF and XLow coupling-based multi-fineness fusion pollutant diffusion analysis method |
| CN113674135A (en) * | 2021-08-03 | 2021-11-19 | 北京三易思创科技有限公司 | Calculation method for realizing CALPUFF high performance based on workflow |
| CN113838539A (en) * | 2021-09-24 | 2021-12-24 | 国网福建省电力有限公司 | Detection and analysis system and analysis method for fire residues of electrical equipment of transformer substation |
| CN113838539B (en) * | 2021-09-24 | 2023-05-26 | 国网福建省电力有限公司 | System and method for detecting and analyzing fire residues of electrical equipment of transformer substation |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112749478B (en) | Atmospheric pollution traceable diffusion analysis system and method based on Gaussian diffusion model | |
| CN106649960B (en) | Three-dimensional iterative calculation method for atmospheric multi-pollutant environment capacity | |
| CN110674583A (en) | Method for constructing inter-city atmospheric pollution transmission matrix | |
| CN111652975B (en) | Urban building group available solar energy resource evaluation method and system | |
| CN113624921A (en) | Multi-mode comprehensive pollution tracing method | |
| CN109902913A (en) | A kind of identification of atmosphere pollution superposed risk area and evaluation method | |
| CN114580310A (en) | Method for realizing scale reduction processing of WRF (hand-wrenching simulation) wind field based on PALM (PALM fiber laser) | |
| CN105675244A (en) | Experimental device for research on flow field and pollutant dispersion in urban continuous street canyons | |
| CN111709644B (en) | Wind power plant wind resource calculation method utilizing unit SCADA data | |
| CN115420854A (en) | Atmospheric pollutant tracing method based on forward and backward model combination | |
| Song et al. | Spatiotemporal distribution of air pollution characteristics in Jiangsu Province, China | |
| Wei | Research on Reducing Carbon Consumption in Residential Community Spaces as Influenced by Microclimate Environments | |
| Wang et al. | Integrating multiple models into computational fluid dynamics for fine three-dimensional simulation of urban waterfront wind environments: A case study in Hangzhou, China | |
| CN114417467A (en) | Method, device, medium and equipment for establishing ventilation corridor of urban dense area | |
| Wang et al. | Green space layout optimization based on microclimate environment features | |
| CN113777223A (en) | Atmospheric pollutant tracing method and system | |
| Huang et al. | Sensitivity analysis of WRF-CFD-based downscaling methods for evaluation of urban pedestrian-level wind | |
| CN109635317A (en) | A kind of high altitude localities CFD artificially generated terrain topological method | |
| CN110727984A (en) | Method for researching mechanism between land utilization change and urban thermal environment | |
| CN112767194A (en) | Analysis method for solar energy utilization potential of building group | |
| Mortensen et al. | Comparison of resource and energy yield assessment procedures | |
| CN103823976A (en) | Photo-thermal environment calculating method of solar greenhouse | |
| CN111123409A (en) | Atmospheric turbulence structural constant forecasting method based on mesoscale numerical model WRF | |
| Gastón et al. | Wind resources map of Spain at mesoscale. Methodology and validation | |
| Graabak et al. | Developing a wind and solar power data model for Europe with high spatial-temporal resolution |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20200110 |