CN111178786B - Emission source position determining method and system for guaranteeing regional air quality - Google Patents
Emission source position determining method and system for guaranteeing regional air quality Download PDFInfo
- Publication number
- CN111178786B CN111178786B CN202010016852.5A CN202010016852A CN111178786B CN 111178786 B CN111178786 B CN 111178786B CN 202010016852 A CN202010016852 A CN 202010016852A CN 111178786 B CN111178786 B CN 111178786B
- Authority
- CN
- China
- Prior art keywords
- grid
- pollutant
- area
- research
- sensitive
- 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
- 238000000034 method Methods 0.000 title claims abstract description 60
- 239000003344 environmental pollutant Substances 0.000 claims abstract description 221
- 231100000719 pollutant Toxicity 0.000 claims abstract description 221
- 239000000356 contaminant Substances 0.000 claims abstract description 30
- 230000008569 process Effects 0.000 claims abstract description 12
- 238000009825 accumulation Methods 0.000 claims abstract description 11
- 238000009792 diffusion process Methods 0.000 claims abstract description 10
- 230000002776 aggregation Effects 0.000 claims description 23
- 238000004220 aggregation Methods 0.000 claims description 23
- 238000012876 topography Methods 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 description 3
- CBENFWSGALASAD-UHFFFAOYSA-N Ozone Chemical compound [O-][O+]=O CBENFWSGALASAD-UHFFFAOYSA-N 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification 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
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06395—Quality analysis or management
-
- 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)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Educational Administration (AREA)
- Economics (AREA)
- Development Economics (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Primary Health Care (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a method for determining the position of an emission source for guaranteeing the air quality of an area, which comprises the following steps: the air quality model simulates the pollutant discharge process of each grid in the research area and the diffusion process of pollutants under the influence of meteorological data of a corresponding research time period to form a pollutant concentration distribution field; identifying a contaminant gathering area and a sensitive area; generating a pollutant gathering area layer and a sensitive area layer; and removing the pollutant accumulation zone layer and the sensitive zone layer in the research zone layer, wherein the residual zone is a final determined discharge source setting position for guaranteeing the air quality of the zone. The advantages are that: only from the influence of topography and weather, two main influencing factors of the atmospheric environment quality of the influence area are determined: a pollutant collection zone and a receptor sensitive zone, and thereby determine the location of the pollutant discharge source. The complexity of pollutant emission source position determination is simplified, the time for pollutant emission source position determination is shortened, and the efficiency of pollutant emission source position determination is improved.
Description
Technical Field
The invention belongs to the technical field of discharge source site selection, and particularly relates to a method and a system for determining a discharge source position for guaranteeing regional air quality.
Background
The emission source position determining method for guaranteeing the air quality of the area is used for determining the position of a pollutant emission source on the premise of guaranteeing the air quality of the area, wherein the pollutant emission source comprises but is not limited to an industrial park, a holiday concentrated firework setting-off place and the like.
In the prior art, a very complex method is needed to determine the position of the pollutant emission source, and the problem of low efficiency of the pollutant emission source position determination process exists.
Disclosure of Invention
Aiming at the defects existing in the prior art, the invention provides a method and a system for determining the position of an emission source for guaranteeing the air quality of a region, which can effectively solve the problems.
The technical scheme adopted by the invention is as follows:
the invention provides a method for determining the position of an emission source for guaranteeing the air quality of a region, which comprises the following steps:
step 1, gridding a research area to form a plurality of grids; setting a grid number of each grid; obtaining a grid attribute of each grid;
step 2, setting the same pollutant emission parameters for each grid based on the pollutants studied at the time; storing the corresponding relation among the grid number, the grid attribute and the pollutant emission parameter into a pollutant emission source list to form a final pollutant emission source list;
step 3, determining a research time period of the current research; then acquiring meteorological data corresponding to a research time period and a research area;
step 4, an air quality model is established, and the pollutant emission source list obtained in the step 2 and the meteorological data obtained in the step 3 are input into the air quality model;
the air quality model simulates the pollutant discharge process of each grid in the research area and the diffusion process of pollutants under the influence of meteorological data of the corresponding research time period, tracks the movement track of the pollutants discharged by each grid, and obtains the final concentration of the pollutants corresponding to each grid to form a pollutant concentration distribution field;
step 5, setting a first pollutant concentration threshold value; sorting the final concentration of the pollutants of each grid from small to large, identifying all grids with final concentration of the pollutants larger than a first pollutant concentration threshold value, and forming a pollutant aggregation area by the identified grids;
step 6, determining sensitive points of a research area, wherein the number of the sensitive points is more than or equal to 1; for each sensitive point, the corresponding sensitive area is determined by adopting the following method:
step 6.1, determining a grid to which the sensitive point belongs, namely grid A; the grid A corresponds to the final concentration of the pollutant;
step 6.2, tracing the final concentration of the pollutants of the grid A determined in the step 6.1 by adopting the air quality model to obtain a grid set contributing to the final concentration of the pollutants of the grid A;
step 6.3, sorting each grid in the grid set according to the concentration contribution to the final concentration of the pollutants of the grid A from small to large;
step 6.4, setting a second pollutant concentration threshold value; identifying grids with concentration contribution larger than a second pollutant concentration threshold value from the grid set, wherein the identified grids form a sensitive area;
step 7, generating a pollutant aggregation area layer from the pollutant aggregation area in the research area by using a geographic information system; generating a sensitive area layer from a sensitive area in a research area;
step 8, generating a research area layer from the research area; and removing the pollutant accumulation zone layer and the sensitive zone layer in the research zone layer, wherein the residual zone is a final determined discharge source setting position for guaranteeing the air quality of the zone.
Preferably, in step 1, the grid attribute includes longitude and latitude coordinates of a grid center point and altitude of the grid center point.
Preferably, the step 5 specifically comprises:
step 5.1, sorting the final concentration of pollutants of each grid from small to large;
step 5.2, determining a first percentage L 1 ；
Step 5.3, calculating (n-1) and L 1 Is a product value of (2); wherein n is the total number of grids in the research area;
step 5.4, the integer in the product value is recorded as A 1 finishing The method comprises the steps of carrying out a first treatment on the surface of the The decimal in the product value is noted as B 1 decimal ；
Step 5.5, in the ranking in step 5.1, locating to rank in A 1 finishing The final concentration of contaminants at the site, designated C A1 The method comprises the steps of carrying out a first treatment on the surface of the Positioned to be ordered in A 1 finishing The final concentration of contaminant at +1 is designated C A1+1 ；
Step 5.6, calculate (1-B 1 decimal )C A1 +B 1 decimal C A1+1 Is denoted as R 1L ；
Step 5.7, R 1L Namely a first pollutant concentration threshold value; thus, it is identified that the final concentration of the contaminant is greater than the first contaminant concentration threshold R 1L The identified grid forms a contaminant gathering area.
Preferably, step 6.4 specifically comprises:
step 6.4.1, sorting the concentration contributions of each grid in the grid set from small to large;
step 6.4.2 determining a second percentage L 2 ；
Step 6.4.3, calculate (m-1) and L 2 Is a product value of (2); wherein m is the total number of grids contained in the grid set;
step 6.4.4, the integer in the product value is denoted A 2 finishing The method comprises the steps of carrying out a first treatment on the surface of the The decimal in the product value is noted as B 2 decimal ；
Step 6.4.5, in the ordering of step 6.3, locating to order at A 2 finishing The concentration contribution of bits, noted C A2 The method comprises the steps of carrying out a first treatment on the surface of the Positioned to be ordered in A 2 finishing The final concentration of contaminant at +1 is designated C A2+1 ；
Step 6.4.6, calculate (1-B 2 decimal )C A2 +B 2 decimal C A2+1 Is denoted as R 2L ；
Step 6.4.7, R 2L Namely a second pollutant concentration threshold value; thus, it is identified that the concentration contribution is greater than the second contaminant concentration threshold R 2L Is aware of all grids of (1)The other grids form sensitive areas.
The invention also provides an emission source position determining system of an emission source position determining method for guaranteeing regional air quality, which comprises the following steps:
the gridding module is used for gridding the research area to form a plurality of grids; setting a grid number of each grid; obtaining a grid attribute of each grid;
the pollutant emission source list generation module is used for setting the same pollutant emission parameters for each grid based on the pollutants studied at the time; storing the corresponding relation among the grid number, the grid attribute and the pollutant emission parameter into a pollutant emission source list to form a final pollutant emission source list;
the acquisition module is used for determining a research time period of the current research; then acquiring meteorological data corresponding to a research time period and a research area;
the air quality model building module is used for building an air quality model;
the air quality model operation module is used for inputting a pollutant emission source list and meteorological data into the air quality model; the air quality model simulates the pollutant discharge process of each grid in the research area and the diffusion process of pollutants under the influence of meteorological data of the corresponding research time period, tracks the movement track of the pollutants discharged by each grid, and obtains the final concentration of the pollutants corresponding to each grid to form a pollutant concentration distribution field;
the pollutant accumulation zone identification module is used for setting a first pollutant concentration threshold value; sorting the final concentration of the pollutants of each grid from small to large, identifying all grids with final concentration of the pollutants larger than a first pollutant concentration threshold value, and forming a pollutant aggregation area by the identified grids;
the sensitive area identification module is used for determining sensitive points of a research area, wherein the number of the sensitive points is more than or equal to 1; for each sensitive point, the corresponding sensitive area is determined by adopting the following method:
1) Determining a grid to which the sensitive points belong, namely a grid A; the grid A corresponds to the final concentration of the pollutant;
2) Tracing the final concentration of the pollutants of the grid A by adopting the air quality model to obtain a grid set contributing to the final concentration of the pollutants of the grid A;
3) Ordering the grids in the grid set according to the concentration contribution to the final concentration of the pollutants in the grid A from small to large;
4) Setting a second contaminant concentration threshold; identifying grids with concentration contribution larger than a second pollutant concentration threshold value from the grid set, wherein the identified grids form a sensitive area;
the pollutant aggregation area layer generation module is used for generating a pollutant aggregation area layer from the pollutant aggregation area in the research area by using the geographic information system;
the sensitive area layer generation module is used for generating a sensitive area layer from a sensitive area in a research area by using a geographic information system;
the emission source setting position determining module is used for generating a research area image layer from a research area; and removing the pollutant accumulation zone layer and the sensitive zone layer in the research zone layer, wherein the residual zone is a final determined discharge source setting position for guaranteeing the air quality of the zone.
The emission source position determining method and system for guaranteeing regional air quality provided by the invention have the following advantages:
according to the invention, by avoiding the influence of complex emission source arrangement, two main influencing factors of the atmospheric environment quality of an influence area are determined only under the influence of terrain and weather: a pollutant collection zone and a receptor sensitive zone, and thereby determine the location of the pollutant discharge source. The method simplifies the complexity of determining the position of the pollutant discharge source, shortens the time for determining the position of the pollutant discharge source and improves the efficiency of determining the position of the pollutant discharge source.
Drawings
Fig. 1 is a schematic flow chart of an emission source position determining method for guaranteeing regional air quality.
Detailed Description
In order to make the technical problems, technical schemes and beneficial effects solved by the invention more clear, the invention is further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the invention.
The invention provides a method for determining the position of an emission source for guaranteeing regional air quality, which relates to the technical fields of an air quality model, a geographic information system and an atmosphere planning and evaluating, and is applied to the fields of regional industry layout, industrial park site selection, holiday centralized firework setting site selection, incineration site selection, atmospheric environment regional planning site selection and the like. The method can provide data support for industrial park layout site selection, holiday centralized setting-off, incineration site selection and atmospheric environment area planning, avoids complex emission source setting, simplifies site selection flow and shortens site selection time in the process of guaranteeing regional air quality site selection.
The location process of the pollutant emission source is an optimal zone selection process taking into account some kind or several kinds of limiting factors. The method determines two factors with great influence on regional air quality, namely a pollutant gathering area and a sensitive area; wherein the pollutant aggregation zone refers to: besides the influence of pollutant discharge amount, the method is only influenced by topography, meteorological factors and diffusion characteristics of pollutant species, and is easy to accumulate in areas. Sensitive area: based on the sensitive points in the research area, the air quality model is used for tracing, and the area with larger pollutant concentration contribution value is the sensitive area.
Generally, the location of the pollutant emission source needs to be input with pollutant emission parameters after the position of the pollutant emission source is primarily determined, so that pollutant diffusion simulation is performed, the method is complex, and the results of different models are greatly different. The method avoids the influence of different models, and only finds out the area where the pollutants are easy to gather, namely the pollutant gathering area, under the influence of the terrain and the weather; and removing the pollutant accumulation area and other areas of the sensitive area, namely recommending sites for the finally determined pollutant emission source.
Therefore, the emission source position determining method for guaranteeing the regional air quality disclosed by the invention only determines two main influencing factors influencing the regional atmospheric environment quality under the influence of terrain and weather by avoiding the influence of complex emission source setting: a pollutant collection zone and a receptor sensitive zone, and thereby determine the location of the pollutant discharge source. The method simplifies the complexity of determining the position of the pollutant discharge source, shortens the time for determining the position of the pollutant discharge source and improves the efficiency of determining the position of the pollutant discharge source.
Referring to fig. 1, an emission source position determining method for guaranteeing regional air quality includes the steps of:
step 1, gridding a research area to form a plurality of grids; setting a grid number of each grid; obtaining a grid attribute of each grid; wherein the grid attribute comprises longitude and latitude coordinates of the grid center point and altitude of the grid center point.
Step 2, setting the same pollutant emission parameters for each grid based on the pollutants, such as ozone pollutants, of the study; the pollutant emission parameter may be emission quantity and the like. In this application, each grid has the same pollutant emission parameters. That is, each mesh continues to release contaminants outwardly at the set contaminant discharge parameters as the contaminant discharge process simulation is performed using the air quality model.
Storing the corresponding relation among the grid number, the grid attribute and the pollutant emission parameter into a pollutant emission source list to form a final pollutant emission source list;
and gridding a research area, extracting the longitude and latitude and the altitude of each grid in ArcGis, and then carrying out batch parameter setting on the corresponding longitude and latitude and altitude in excel according to the input format of a point source in a calpuff.inp file and replacing the content in the corresponding inp file.
The pollutant emission parameters of different air quality models are input in different forms, and all parameters except geographic information and terrains need to be set uniformly in the invention. And specifically completing a pollutant emission source list according to the actual air quality model.
The pollutant emission source list is a readable list generated based on pollutant emission parameters.
Storing the corresponding relation among the grid number, the grid attribute and the pollutant emission parameter into a pollutant emission source list to form a final pollutant emission source list;
step 3, determining a research time period of the current research; then acquiring meteorological data corresponding to a research time period and a research area;
meteorological data may be obtained through WRF model output.
Step 4, an air quality model is established, and the pollutant emission source list obtained in the step 2 and the meteorological data obtained in the step 3 are input into the air quality model;
the air quality model can be calpuff, cmaq, wrf-chem, etc.
The air quality model simulates the pollutant discharge process of each grid in the research area and the diffusion process of pollutants under the influence of meteorological data of the corresponding research time period, tracks the movement track of the pollutants discharged by each grid, and obtains the final concentration of the pollutants corresponding to each grid to form a pollutant concentration distribution field;
specifically, the research area simulates an aerial image field by using WRF, wrout is processed into 3D.dat which can be identified by calwrf, and the starting point and the grid number of the research area are configured in calmet to generate met.dat. The inp file is loaded into the calpuff model and the met. Dat is loaded and then the pollutant diffusion simulation is performed.
Step 5, setting a first pollutant concentration threshold value; sorting the final concentration of the pollutants of each grid from small to large, identifying all grids with final concentration of the pollutants larger than a first pollutant concentration threshold value, and forming a pollutant aggregation area by the identified grids;
the step 5 is specifically as follows:
step 5.1, sorting the final concentration of pollutants of each grid from small to large;
step 5.2, determining a first percentage L 1 The method comprises the steps of carrying out a first treatment on the surface of the For example, a first percentage L 1 0.9, 0.8, etcSpecifically, the device is arranged according to actual needs.
Step 5.3, calculating (n-1) and L 1 Is a product value of (2); wherein n is the total number of grids in the research area;
step 5.4, the integer in the product value is recorded as A 1 finishing The method comprises the steps of carrying out a first treatment on the surface of the The decimal in the product value is noted as B 1 decimal ；
Step 5.5, in the ranking in step 5.1, locating to rank in A 1 finishing The final concentration of contaminants at the site, designated C A1 The method comprises the steps of carrying out a first treatment on the surface of the Positioned to be ordered in A 1 finishing The final concentration of contaminant at +1 is designated C A1+1 ；
Step 5.6, calculate (1-B 1 decimal )C A1 +B 1 decimal C A1+1 Is denoted as R 1L ；
For example, assume that the total number of study area grids n is 1000; l (L) 1 0.9; then (n-1) and L 1 The product value of (2) is 899.1. Thus, A 1 finishing ＝899；B 1 decimal =0.1. In the sorting in step 5.1, the final concentration of the contaminant positioned at position 899 in the sorting is C A1 The method comprises the steps of carrying out a first treatment on the surface of the Positioned to be ordered in A 1 finishing The final concentration of contaminant at +1=900, designated C A1+1 The method comprises the steps of carrying out a first treatment on the surface of the From this, R can be calculated 1L 。
Step 5.7, R 1L Namely a first pollutant concentration threshold value; thus, it is identified that the final concentration of the contaminant is greater than the first contaminant concentration threshold R 1L The identified grid forms a contaminant gathering area.
Step 6, determining sensitive points of a research area, wherein the number of the sensitive points is more than or equal to 1; for each sensitive point, the corresponding sensitive area is determined by adopting the following method:
the sensitive points are determined according to the land planning of the research area, and can be places such as schools, hospitals and the like.
Step 6.1, determining a grid to which the sensitive point belongs, namely grid A; the grid A corresponds to the final concentration of the pollutant;
step 6.2, tracing the final concentration of the pollutants of the grid A determined in the step 6.1 by adopting the air quality model to obtain a grid set contributing to the final concentration of the pollutants of the grid A;
step 6.3, sorting each grid in the grid set according to the concentration contribution to the final concentration of the pollutants of the grid A from small to large;
step 6.4, setting a second pollutant concentration threshold value; identifying grids with concentration contribution larger than a second pollutant concentration threshold value from the grid set, wherein the identified grids form a sensitive area;
the step 6.4 specifically comprises the following steps:
step 6.4.1, sorting the concentration contributions of each grid in the grid set from small to large;
step 6.4.2 determining a second percentage L 2 ；
Step 6.4.3, calculate (m-1) and L 2 Is a product value of (2); wherein m is the total number of grids contained in the grid set;
step 6.4.4, the integer in the product value is denoted A 2 finishing The method comprises the steps of carrying out a first treatment on the surface of the The decimal in the product value is noted as B 2 decimal ；
Step 6.4.5, in the ordering of step 6.3, locating to order at A 2 finishing The concentration contribution of bits, noted C A2 The method comprises the steps of carrying out a first treatment on the surface of the Positioned to be ordered in A 2 finishing The final concentration of contaminant at +1 is designated C A2+1 ；
Step 6.4.6, calculate (1-B 2 decimal )C A2 +B 2 decimal C A2+1 Is denoted as R 2L ；
Step 6.4.7, R 2L Namely a second pollutant concentration threshold value; thus, it is identified that the concentration contribution is greater than the second contaminant concentration threshold R 2L The identified grid forms the sensitive area.
Specifically, DAT data are converted into lattice point data and loaded into ArcGis, the lattice point data and the fishing net data are subjected to spatial join operation, the lattice file data are linked to the fishing net to obtain A1 and A2, and the value of A1 and A2 is greater than or equal to '90 percentile' by using spatial analyst tools-extract by attributes functions of ArcGis and matching with simple SQL sentences, so that a pollutant aggregation area and a sensitive area are respectively identified.
In practical application, based on the definition of the sensitive area, the Arcgis is used for selecting the sensitive area, the grid contained in the sensitive area is extracted, the trace back function of the calpost is used for tracing the sensitive points of the output data of the calpuff, and particularly, attention is required to be paid to the trace back that only one sensitive point can be traced at a time, a batch script is required to be written by the shell under linux, and the tracing result is also subjected to 90 percentiles.
Step 7, generating a pollutant aggregation area layer from the pollutant aggregation area in the research area by using a geographic information system; generating a sensitive area layer from a sensitive area in a research area;
step 8, generating a research area layer from the research area; and removing the pollutant accumulation zone layer and the sensitive zone layer in the research zone layer, wherein the residual zone is a final determined discharge source setting position for guaranteeing the air quality of the zone.
Specifically, reintroducing the pollutant gathering area and the sensitive area into ArcGi, and performing spat i al j oin operation on the point data and the fishing net data to obtain a layer file which should be avoided by the setting position of the emission source
The invention also provides an emission source position determining system of an emission source position determining method for guaranteeing regional air quality, which comprises the following steps:
the gridding module is used for gridding the research area to form a plurality of grids; setting a grid number of each grid; obtaining a grid attribute of each grid;
the pollutant emission source list generation module is used for setting the same pollutant emission parameters for each grid based on the pollutants studied at the time; storing the corresponding relation among the grid number, the grid attribute and the pollutant emission parameter into a pollutant emission source list to form a final pollutant emission source list;
the acquisition module is used for determining a research time period of the current research; then acquiring meteorological data corresponding to a research time period and a research area;
the air quality model building module is used for building an air quality model;
the air quality model operation module is used for inputting a pollutant emission source list and meteorological data into the air quality model; the air quality model simulates the pollutant discharge process of each grid in the research area and the diffusion process of pollutants under the influence of meteorological data of the corresponding research time period, tracks the movement track of the pollutants discharged by each grid, and obtains the final concentration of the pollutants corresponding to each grid to form a pollutant concentration distribution field;
the pollutant accumulation zone identification module is used for setting a first pollutant concentration threshold value; sorting the final concentration of the pollutants of each grid from small to large, identifying all grids with final concentration of the pollutants larger than a first pollutant concentration threshold value, and forming a pollutant aggregation area by the identified grids;
the sensitive area identification module is used for determining sensitive points of a research area, wherein the number of the sensitive points is more than or equal to 1; for each sensitive point, the corresponding sensitive area is determined by adopting the following method:
1) Determining a grid to which the sensitive points belong, namely a grid A; the grid A corresponds to the final concentration of the pollutant;
2) Tracing the final concentration of the pollutants of the grid A by adopting the air quality model to obtain a grid set contributing to the final concentration of the pollutants of the grid A;
3) Ordering the grids in the grid set according to the concentration contribution to the final concentration of the pollutants in the grid A from small to large;
4) Setting a second contaminant concentration threshold; identifying grids with concentration contribution larger than a second pollutant concentration threshold value from the grid set, wherein the identified grids form a sensitive area;
the pollutant aggregation area layer generation module is used for generating a pollutant aggregation area layer from the pollutant aggregation area in the research area by using the geographic information system;
the sensitive area layer generation module is used for generating a sensitive area layer from a sensitive area in a research area by using a geographic information system;
the emission source setting position determining module is used for generating a research area image layer from a research area; and removing the pollutant accumulation zone layer and the sensitive zone layer in the research zone layer, wherein the residual zone is a final determined discharge source setting position for guaranteeing the air quality of the zone.
The emission source position determining method and system for guaranteeing regional air quality provided by the invention have the following advantages:
according to the invention, by avoiding the influence of complex emission source arrangement, two main influencing factors of the atmospheric environment quality of an influence area are determined only under the influence of terrain and weather: the pollutant collecting area and the receptor sensitive area realize the determination of the area with unfavorable air quality only under the influence of terrain and climate, and the determination of the position of pollutant emission source. The method simplifies the complexity of determining the position of the pollutant discharge source, shortens the time for determining the position of the pollutant discharge source and improves the efficiency of determining the position of the pollutant discharge source.
The foregoing is merely a preferred embodiment of the present invention and it should be noted that modifications and adaptations to those skilled in the art may be made without departing from the principles of the present invention, which is also intended to be covered by the present invention.
Claims (3)
1. An emission source position determining method for guaranteeing regional air quality is characterized by comprising the following steps:
step 1, gridding a research area to form a plurality of grids; setting a grid number of each grid; obtaining a grid attribute of each grid;
step 2, setting the same pollutant emission parameters for each grid based on the pollutants studied at the time; storing the corresponding relation among the grid number, the grid attribute and the pollutant emission parameter into a pollutant emission source list to form a final pollutant emission source list;
step 3, determining a research time period of the current research; then acquiring meteorological data corresponding to a research time period and a research area;
step 4, an air quality model is established, and the pollutant emission source list obtained in the step 2 and the meteorological data obtained in the step 3 are input into the air quality model;
the air quality model simulates the pollutant discharge process of each grid in the research area and the diffusion process of pollutants under the influence of meteorological data of the corresponding research time period, tracks the movement track of the pollutants discharged by each grid, and obtains the final concentration of the pollutants corresponding to each grid to form a pollutant concentration distribution field;
step 5, setting a first pollutant concentration threshold value; sorting the final concentration of the pollutants of each grid from small to large, identifying all grids with final concentration of the pollutants larger than a first pollutant concentration threshold value, and forming a pollutant aggregation area by the identified grids;
the step 5 is specifically as follows:
step 5.1, sorting the final concentration of pollutants of each grid from small to large;
step 5.2, determining a first percentage L 1 ；
Step 5.3, calculating (n-1) and L 1 Is a product value of (2); wherein n is the total number of grids in the research area;
step 5.4, the integer part of the product value is denoted as A 1 finishing The method comprises the steps of carrying out a first treatment on the surface of the The decimal part in the product value is denoted as B 1 decimal ；
Step 5.5, in the ranking in step 5.1, locating to rank in A 1 finishing The final concentration of contaminants at the site, designated C A1 The method comprises the steps of carrying out a first treatment on the surface of the Positioned to be ordered in A 1 finishing The final concentration of contaminant at +1 is designated C A1+1 ；
Step 5.6, calculate (1-B 1 decimal ）C A1 + B 1 decimal C A1+1 Is denoted as R 1L ；
Step 5.7, R 1L Namely a first pollutant concentration threshold value; thus, it is identified that the final concentration of the contaminant is greater than the first contaminant concentration threshold R 1L The identified grids form a contaminant gathering area;
step 6, determining sensitive points of a research area, wherein the number of the sensitive points is more than or equal to 1; for each sensitive point, the corresponding sensitive area is determined by adopting the following method:
step 6.1, determining a grid to which the sensitive point belongs, namely grid A; the grid A corresponds to the final concentration of the pollutant;
step 6.2, tracing the final concentration of the pollutants of the grid A determined in the step 6.1 by adopting the air quality model to obtain a grid set contributing to the final concentration of the pollutants of the grid A;
step 6.3, sorting each grid in the grid set according to the concentration contribution to the final concentration of the pollutants of the grid A from small to large;
step 6.4, setting a second pollutant concentration threshold value; identifying grids with concentration contribution larger than a second pollutant concentration threshold value from the grid set, wherein the identified grids form a sensitive area;
the step 6.4 specifically comprises the following steps:
step 6.4.1, sorting the concentration contributions of each grid in the grid set from small to large;
step 6.4.2 determining a second percentage L 2 ；
Step 6.4.3, calculate (m-1) and L 2 Is a product value of (2); wherein m is the total number of grids contained in the grid set;
step 6.4.4, the integer part of the product value is denoted A 2 finishing The method comprises the steps of carrying out a first treatment on the surface of the The decimal part in the product value is denoted as B 2 decimal ；
Step 6.4.5, in the ordering of step 6.3, locating to order at A 2 finishing The concentration contribution of bits, noted C A2 The method comprises the steps of carrying out a first treatment on the surface of the Positioned to be ordered in A 2 finishing The final concentration of contaminant at +1 is designated C A2+1 ；
Step 6.4.6, calculate (1-B 2 decimal ）C A2 + B 2 decimal C A2+1 Is denoted as R 2L ；
Step 6.4.7, R 2L Namely a second pollutant concentration threshold value; thus, it is identified that the concentration contribution is greater than the second contaminant concentration threshold R 2L The identified grids form a sensitive area;
step 7, generating a pollutant aggregation area layer from the pollutant aggregation area in the research area by using a geographic information system; generating a sensitive area layer from a sensitive area in a research area;
step 8, generating a research area layer from the research area; and removing the pollutant accumulation zone layer and the sensitive zone layer in the research zone layer, wherein the residual zone is a final determined discharge source setting position for guaranteeing the air quality of the zone.
2. The method for determining the location of an emission source for guaranteeing regional air quality according to claim 1, wherein in step 1, the grid attribute includes longitude and latitude coordinates of a grid center point and an altitude of the grid center point.
3. An emissions source position determining system of an emissions source position determining method for guaranteeing regional air quality as claimed in any one of claims 1-2, comprising:
the gridding module is used for gridding the research area to form a plurality of grids; setting a grid number of each grid; obtaining a grid attribute of each grid;
the pollutant emission source list generation module is used for setting the same pollutant emission parameters for each grid based on the pollutants studied at the time; storing the corresponding relation among the grid number, the grid attribute and the pollutant emission parameter into a pollutant emission source list to form a final pollutant emission source list;
the acquisition module is used for determining a research time period of the current research; then acquiring meteorological data corresponding to a research time period and a research area;
the air quality model building module is used for building an air quality model;
the air quality model operation module is used for inputting a pollutant emission source list and meteorological data into the air quality model; the air quality model simulates the pollutant discharge process of each grid in the research area and the diffusion process of pollutants under the influence of meteorological data of the corresponding research time period, tracks the movement track of the pollutants discharged by each grid, and obtains the final concentration of the pollutants corresponding to each grid to form a pollutant concentration distribution field;
the pollutant accumulation zone identification module is used for setting a first pollutant concentration threshold value; sorting the final concentration of the pollutants of each grid from small to large, identifying all grids with final concentration of the pollutants larger than a first pollutant concentration threshold value, and forming a pollutant aggregation area by the identified grids;
the sensitive area identification module is used for determining sensitive points of a research area, wherein the number of the sensitive points is more than or equal to 1; for each sensitive point, the corresponding sensitive area is determined by adopting the following method:
1) Determining a grid to which the sensitive points belong, namely a grid A; the grid A corresponds to the final concentration of the pollutant;
2) Tracing the final concentration of the pollutants of the grid A by adopting the air quality model to obtain a grid set contributing to the final concentration of the pollutants of the grid A;
3) Ordering the grids in the grid set according to the concentration contribution to the final concentration of the pollutants in the grid A from small to large;
4) Setting a second contaminant concentration threshold; identifying grids with concentration contribution larger than a second pollutant concentration threshold value from the grid set, wherein the identified grids form a sensitive area;
the pollutant aggregation area layer generation module is used for generating a pollutant aggregation area layer from the pollutant aggregation area in the research area by using the geographic information system;
the sensitive area layer generation module is used for generating a sensitive area layer from a sensitive area in a research area by using a geographic information system;
the emission source setting position determining module is used for generating a research area image layer from a research area; and removing the pollutant accumulation zone layer and the sensitive zone layer in the research zone layer, wherein the residual zone is a final determined discharge source setting position for guaranteeing the air quality of the zone.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010016852.5A CN111178786B (en) | 2020-01-08 | 2020-01-08 | Emission source position determining method and system for guaranteeing regional air quality |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010016852.5A CN111178786B (en) | 2020-01-08 | 2020-01-08 | Emission source position determining method and system for guaranteeing regional air quality |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111178786A CN111178786A (en) | 2020-05-19 |
| CN111178786B true CN111178786B (en) | 2023-04-25 |
Family
ID=70654526
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010016852.5A Active CN111178786B (en) | 2020-01-08 | 2020-01-08 | Emission source position determining method and system for guaranteeing regional air quality |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111178786B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111612064B (en) * | 2020-05-20 | 2021-02-19 | 中科三清科技有限公司 | PM2.5Method and device for tracing pollution air mass, electronic equipment and storage medium |
| CN111860996B (en) * | 2020-07-15 | 2021-09-03 | 中科三清科技有限公司 | Emission source data generation method and device of air quality model and medium |
| CN112214913A (en) * | 2020-11-16 | 2021-01-12 | 中科三清科技有限公司 | Method and device for identifying dominant precursor of ozone, electronic equipment and storage medium |
| CN116757366B (en) * | 2023-08-14 | 2023-11-10 | 中科三清科技有限公司 | Method, device, medium and electronic equipment for locating high-pollution emission pollution source |
| CN117408856B (en) * | 2023-12-13 | 2024-03-26 | 浙江省生态环境监测中心(浙江省生态环境信息中心) | Pollutant tracing method and device, storage medium and electronic equipment |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN201697870U (en) * | 2010-06-24 | 2011-01-05 | 中科宇图天下科技有限公司 | Automatic flue gas blackness monitoring and analyzing system |
| CN102628852A (en) * | 2012-03-13 | 2012-08-08 | 北京工业大学 | Atmospheric pollution source grading method based on pollutant source identification technology |
| CN105403664A (en) * | 2015-10-19 | 2016-03-16 | 电力规划设计总院 | WRF-CHEM-based large point pollution source atmosphere environment influence evaluating method |
| CN106503447A (en) * | 2016-10-31 | 2017-03-15 | 中科宇图科技股份有限公司 | A kind of non-point pollution emission inventories immediate processing method |
| KR20180000619A (en) * | 2016-06-23 | 2018-01-03 | (주)해동기술개발공사 | Soil loss evaluation method based GIS |
| CN107563646A (en) * | 2017-09-05 | 2018-01-09 | 环境保护部华南环境科学研究所 | A kind of management system of burning city domestic garbage treatment facility to population health risk |
| CN107832910A (en) * | 2017-09-30 | 2018-03-23 | 东南大学 | A kind of method that evaluation path traffic characteristic influences on carbon monoxide emission concentration |
| CN207623771U (en) * | 2017-12-27 | 2018-07-17 | 国网冀北电力有限公司电力科学研究院 | A kind of coal-burning power plant's atmosphere pollution diffusion monitoring system |
| CN109711681A (en) * | 2018-12-12 | 2019-05-03 | 天津大学 | A kind of region-wide Environment Risk Assessment Methods based on ArcGIS |
| CN109902913A (en) * | 2019-01-10 | 2019-06-18 | 天津大学 | A kind of identification of atmosphere pollution superposed risk area and evaluation method |
| CN110009226A (en) * | 2019-04-03 | 2019-07-12 | 中国科学院地理科学与资源研究所 | A kind of evaluation and test of contaminated sites risk class and sensitive receptors influence identifying system |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100121680A1 (en) * | 2008-09-09 | 2010-05-13 | Propel Biofuels, Inc. | Method and system for identifying installation sites of alternative fuel stations |
| US20180144353A1 (en) * | 2016-11-21 | 2018-05-24 | Propel Biofuels, Inc. | Identifying installation sites for alternative fuel stations |
-
2020
- 2020-01-08 CN CN202010016852.5A patent/CN111178786B/en active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN201697870U (en) * | 2010-06-24 | 2011-01-05 | 中科宇图天下科技有限公司 | Automatic flue gas blackness monitoring and analyzing system |
| CN102628852A (en) * | 2012-03-13 | 2012-08-08 | 北京工业大学 | Atmospheric pollution source grading method based on pollutant source identification technology |
| CN105403664A (en) * | 2015-10-19 | 2016-03-16 | 电力规划设计总院 | WRF-CHEM-based large point pollution source atmosphere environment influence evaluating method |
| KR20180000619A (en) * | 2016-06-23 | 2018-01-03 | (주)해동기술개발공사 | Soil loss evaluation method based GIS |
| CN106503447A (en) * | 2016-10-31 | 2017-03-15 | 中科宇图科技股份有限公司 | A kind of non-point pollution emission inventories immediate processing method |
| CN107563646A (en) * | 2017-09-05 | 2018-01-09 | 环境保护部华南环境科学研究所 | A kind of management system of burning city domestic garbage treatment facility to population health risk |
| CN107832910A (en) * | 2017-09-30 | 2018-03-23 | 东南大学 | A kind of method that evaluation path traffic characteristic influences on carbon monoxide emission concentration |
| CN207623771U (en) * | 2017-12-27 | 2018-07-17 | 国网冀北电力有限公司电力科学研究院 | A kind of coal-burning power plant's atmosphere pollution diffusion monitoring system |
| CN109711681A (en) * | 2018-12-12 | 2019-05-03 | 天津大学 | A kind of region-wide Environment Risk Assessment Methods based on ArcGIS |
| CN109902913A (en) * | 2019-01-10 | 2019-06-18 | 天津大学 | A kind of identification of atmosphere pollution superposed risk area and evaluation method |
| CN110009226A (en) * | 2019-04-03 | 2019-07-12 | 中国科学院地理科学与资源研究所 | A kind of evaluation and test of contaminated sites risk class and sensitive receptors influence identifying system |
Non-Patent Citations (3)
| Title |
|---|
| 刘永红 ; 王成 ; 余志 ; .佛山市禅城区火电企业减排效果研究.科技导报.2011,(第11期),全文. * |
| 吕楠 ; 赵敬源 ; 张鹏 ; .基于GIS技术的城市大气污染物扩散模型研究.西安建筑科技大学学报(自然科学版).2019,(第05期),全文. * |
| 梁津 ; 刘志红 ; 姚琳 ; .基于WRF/CMAQ的成都空气质量模拟与预报.高原山地气象研究.2016,(第03期),全文. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111178786A (en) | 2020-05-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111178786B (en) | Emission source position determining method and system for guaranteeing regional air quality | |
| Sohl et al. | Spatially explicit modeling of 1992–2100 land cover and forest stand age for the conterminous United States | |
| Falloon et al. | Regional estimates of carbon sequestration potential: linking the Rothamsted Carbon Model to GIS databases | |
| CN110598953A (en) | Space-time correlation air quality prediction method | |
| CN109060617B (en) | Simulation and analysis method for influences of different agricultural ammonia emission situations on PM2.5 concentration | |
| CN103336894B (en) | A kind of city plot ratio automatic Mesh Partition Method | |
| CN116011317B (en) | Small-scale near-real-time atmospheric pollution tracing method based on multi-method fusion | |
| CN107016095A (en) | Climate change comprehensive estimation method based on multi-source carbon number evidence | |
| CN112884088A (en) | Forest carbon reserve calculation method based on neural network model | |
| CN117171223A (en) | Microorganism culture scheme recommendation method and system in microorganism repair process | |
| CN116739224B (en) | Method, apparatus and storage medium for determining contribution concentration of heavy point pollution source | |
| CN117036112B (en) | Geographic information system and method for land planning | |
| Issiako et al. | Prospective mapping of land cover and land use in the classified forest of the upper alibori based on satellite imagery | |
| CN115239027B (en) | Method and device for forecasting air quality check set | |
| CN116739232A (en) | Method, device, equipment and storage medium for distinguishing rural sewage treatment mode | |
| CN115860189A (en) | Method and system for optimizing land utilization spatial pattern under low-carbon target | |
| CN113408895B (en) | Ecological quality index construction method and system based on pixel scale | |
| CN115759810A (en) | Enterprise evaluation model iterative optimization method, enterprise evaluation method and assembly line system | |
| Badr et al. | Toward large-scale crop production forecasts for global food security | |
| CN114299238A (en) | Quantitative prediction method for changes of rare animal and plant habitat suitable areas along railway | |
| Gianquintieri et al. | Land use influence on ambient PM2. 5 and ammonia concentrations: Correlation analyses in the Lombardy region, Italy | |
| Haripersaud | Collecting biodiversity | |
| CN113592664B (en) | Crop production space prediction simulation method, device, model and storage medium | |
| CN116049342B (en) | Habitat quality monitoring method and device | |
| CN117972397B (en) | Atmospheric dry-wet sedimentation model simulation optimization monitoring method and medium based on big data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: A Method and System for Determining the Location of Emission Sources to Ensure Regional Air Quality Effective date of registration: 20230904 Granted publication date: 20230425 Pledgee: Zhongguancun Beijing technology financing Company limited by guarantee Pledgor: CHINA SCIENCE MAPUNIVERSE TCHNDOGY Co.,Ltd. Registration number: Y2023990000439 |