CN111126872A - Atmospheric pollution remediation effect evaluation method and system - Google Patents

Atmospheric pollution remediation effect evaluation method and system Download PDF

Info

Publication number
CN111126872A
CN111126872A CN201911398832.2A CN201911398832A CN111126872A CN 111126872 A CN111126872 A CN 111126872A CN 201911398832 A CN201911398832 A CN 201911398832A CN 111126872 A CN111126872 A CN 111126872A
Authority
CN
China
Prior art keywords
pollutant
grid
concentration
enterprise
evaluated
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
Application number
CN201911398832.2A
Other languages
Chinese (zh)
Inventor
马鹏飞
胡奎伟
陈辉
张连华
周春艳
张丽娟
张玉环
陈翠红
翁国庆
王中挺
厉青
毛慧琴
王延龙
赵爱梅
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Satellite Application Center for Ecology and Environment of MEE
Original Assignee
Satellite Application Center for Ecology and Environment of MEE
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Satellite Application Center for Ecology and Environment of MEE filed Critical Satellite Application Center for Ecology and Environment of MEE
Priority to CN201911398832.2A priority Critical patent/CN111126872A/en
Publication of CN111126872A publication Critical patent/CN111126872A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/0004Gaseous mixtures, e.g. polluted air
    • G01N33/0009General constructional details of gas analysers, e.g. portable test equipment
    • G01N33/0062General constructional details of gas analysers, e.g. portable test equipment concerning the measuring method or the display, e.g. intermittent measurement or digital display
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/0004Gaseous mixtures, e.g. polluted air
    • G01N33/0009General constructional details of gas analysers, e.g. portable test equipment
    • G01N33/0062General constructional details of gas analysers, e.g. portable test equipment concerning the measuring method or the display, e.g. intermittent measurement or digital display
    • G01N33/0067General constructional details of gas analysers, e.g. portable test equipment concerning the measuring method or the display, e.g. intermittent measurement or digital display by measuring the rate of variation of the concentration
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10Services
    • G06Q50/26Government or public services

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Tourism & Hospitality (AREA)
  • General Health & Medical Sciences (AREA)
  • Educational Administration (AREA)
  • Development Economics (AREA)
  • Immunology (AREA)
  • Theoretical Computer Science (AREA)
  • Medicinal Chemistry (AREA)
  • Biochemistry (AREA)
  • Analytical Chemistry (AREA)
  • Combustion & Propulsion (AREA)
  • Marketing (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Food Science & Technology (AREA)
  • General Business, Economics & Management (AREA)
  • Pathology (AREA)
  • Primary Health Care (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The embodiment of the invention provides a method and a system for evaluating the effect of atmospheric pollution remediation, which comprises the following steps: dividing a region to be evaluated into grids, and acquiring the proportional reduction of the average concentration of pollutants in the region to be evaluated and the proportional reduction of the average concentration of each pollutant in each grid; acquiring the emission amount of each pollutant of each enterprise in all preset key industries in an area to be evaluated; acquiring the contribution rate of each enterprise to the concentration of each pollutant of the grid; calculating the pollutant concentration reduction caused by the renovation of each enterprise; and acquiring the contribution rate of remediation of each preset key industry to the reduction of the pollutant concentration according to the preset key industry corresponding to each enterprise, the pollutant concentration reduction amount caused by remediation of each enterprise and the proportional reduction amount of the average concentration of all pollutants in the area to be assessed. According to the method, the effect of pollutant concentration reduction caused by enterprise renovation is calculated, and the reduction contribution rate of various preset key industry enterprises to the pollutant concentration of the whole area to be evaluated is obtained.

Description

Atmospheric pollution remediation effect evaluation method and system
Technical Field
The invention relates to the technical field of computers, in particular to an assessment method and an assessment system for an atmospheric pollution remediation effect.
Background
Existing evaluation techniques focus primarily on macroscopic evaluation as follows: 1) the concentration change condition of 6 items of conventional monitoring elements on the ground of a city is utilized, and the 6 items of conventional monitoring elements on the ground comprise PM2.5、PM10、SO2、NO2、O3And CO; 2) simulating the concentration and proportion change of atmospheric particulate component sulfate, nitrate and ammonium salt by using an atmospheric chemical mode; 3) and simulating the influence of meteorological factors and non-meteorological factors on the quality to evaluate the overall air quality improvement condition of each city.
The three technologies only analyze the control effect of each pollutant through the concentration change conditions of a plurality of pollutants, the influence conditions of meteorological factors and non-meteorological factors on the air quality and the change conditions of each component of atmospheric particulate matters, macroscopically evaluate the overall air quality improvement condition of the city, and provide pollution treatment opinions in a large scale range of the region.
The methods cannot go deep into each type of key industry, cannot refine to a single enterprise, and estimate the contribution rate of the gas-involved enterprise reformation which does not go out of the key industry to the improvement of the overall air quality. Therefore, the proposed proposal is relatively general and has poor operability.
Disclosure of Invention
In order to solve the above problems, embodiments of the present invention provide a method and a system for evaluating an atmospheric pollution remediation effect.
In a first aspect, an embodiment of the present invention provides a method for evaluating an effect of remediation of atmospheric pollution, including: dividing a region to be evaluated into a plurality of grids, obtaining the average concentration of each pollutant in the region to be evaluated and the average concentration of each pollutant in each grid, and obtaining the proportional reduction of the average concentration of all pollutants in the region to be evaluated and the proportional reduction of the average concentration of each pollutant in each grid;
acquiring the pollutant emission amount of each enterprise in all preset key industries in the area to be evaluated;
acquiring the contribution rate of each enterprise to the concentration of each pollutant of the grid according to the emission amount of each pollutant of each enterprise and the average concentration of each pollutant of each grid;
calculating pollutant concentration reduction caused by the renovation of each enterprise according to the proportional reduction of each pollutant concentration contribution rate of each enterprise to the grid and each pollutant average concentration of each grid;
and acquiring the contribution rate of the remediation of each preset key industry to the reduction of the pollutant concentration according to the preset key industry corresponding to each enterprise, the pollutant concentration reduction amount caused by the remediation of each enterprise and the proportional reduction amount of the average concentration of all pollutants in the area to be assessed.
Preferably, the dividing the region to be evaluated into a plurality of grids, obtaining the average concentration of each pollutant in the region to be evaluated and the average concentration of each pollutant in each grid, and obtaining the proportional decrease amount of the average concentration of all pollutants in the region to be evaluated and the proportional decrease amount of the average concentration of each pollutant in each grid specifically includes:
acquiring the concentration of each pollutant at each pixel point in the area to be evaluated in the current year and the concentration of each pollutant at each pixel point in the area to be evaluated in the previous year by utilizing a satellite remote sensing inversion technology, and further acquiring the average concentration of each pollutant in the area to be evaluated in the current year and the concentration of each pollutant at each pixel point in the area to be evaluated in the previous year;
dividing the area to be evaluated into a plurality of grids, obtaining the average concentration of each pollutant at each grid in the current year according to the concentration of each pollutant at each pixel point in the area to be evaluated in the current year, and obtaining the average concentration of each pollutant at each grid in the previous year according to the concentration of each pollutant at each pixel point in the area to be evaluated in the previous year;
obtaining the proportional reduction of the average concentration of all pollutants in the area to be evaluated according to the average concentration of each pollutant in the area to be evaluated in the current year and the average concentration of each pollutant in the area to be evaluated in the previous year;
and obtaining the proportional reduction of the average concentration of each pollutant of each grid according to the average concentration of each pollutant of each grid in the current year and the average concentration of each pollutant of each grid in the previous year.
Preferably, the obtaining of the contribution rate of each enterprise to each pollutant concentration of the grid in which each enterprise is located according to each pollutant emission amount of each enterprise and each pollutant average concentration of each grid specifically includes:
for any grid, acquiring each pollutant emission of each enterprise in the grid according to the following formula:
Figure BDA0002347011260000031
wherein E isjRepresenting the average concentration of the jth pollutant, m representing the number of enterprises in all preset key industries of any grid, en,jRepresenting the discharge amount of j pollutant of n enterprise in any grid;
and acquiring the contribution rate of each enterprise to each pollutant concentration of any grid according to the emission rate of each pollutant in any grid and the emission rate of each pollutant of each enterprise in any grid.
Preferably, the method includes obtaining each pollutant concentration contribution rate of each enterprise to any one grid according to the emission amount of each pollutant in any one grid and each pollutant emission amount of each enterprise in any one grid, and accordingly, a specific calculation formula of each pollutant concentration contribution rate of each enterprise in any one grid to any one grid is as follows:
Figure BDA0002347011260000032
wherein n isn,jRepresenting the contribution rate of the nth enterprise in any grid to the j pollutant emission amount of any grid, en,jRepresents the j pollutant discharge amount of the nth enterprise, EjRepresents the total emission amount of the j pollutant of any grid.
Preferably, the calculating the pollutant concentration decrease amount caused by the remediation of each enterprise according to the proportional decrease amount of each pollutant concentration contribution rate of each enterprise to the grid where each enterprise is located and each pollutant average concentration of each grid specifically includes:
Figure BDA0002347011260000033
wherein, DeltasumIndicating the amount of decrease in contaminant concentration caused by each enterprise remediation,Di,j,1representing the concentration of the j-th pollutant in the ith grid of the previous year, Di,j,2Represents the concentration of the j-th pollutant in the ith grid of the current year, ni,jAnd representing the contribution rate of the ith enterprise in any grid to the j pollutant emission amount of any grid.
Preferably, the proportional decrease of the average concentration of all the contaminants in the region to be evaluated is obtained by the following formula:
Figure BDA0002347011260000041
wherein, delta'sumRepresenting the commensurately decreasing amount of the average concentration of all contaminants in said area to be assessed, Di,j,1Represents the concentration of the j-th pollutant in the ith grid of the previous year, Di,j,2Indicating the concentration of the jth pollutant in the ith grid of the current year.
Preferably, each pollutant comprises PM2.5、NO2And SO2
In a second aspect, an embodiment of the present invention provides an evaluation system for an effect of remediating atmospheric pollution, including: the system comprises a comparing module, a calculating module and a judging module, wherein the comparing module is used for dividing a region to be evaluated into a plurality of grids, acquiring the average concentration of each pollutant in the region to be evaluated and the average concentration of each pollutant in each grid, and acquiring the comparing descending amount of the average concentration of all pollutants in the region to be evaluated and the comparing descending amount of the average concentration of each pollutant in each grid;
the enterprise emission module is used for acquiring the emission amount of each pollutant of each enterprise in all preset key industries in the area to be evaluated;
the contribution rate module is used for acquiring the contribution rate of each enterprise to the concentration of each pollutant of the grid according to the emission amount of each pollutant of each enterprise and the average concentration of each pollutant of each grid;
the regulation module is used for calculating the pollutant concentration reduction amount caused by regulation of each enterprise according to the contribution rate of each enterprise to the pollutant concentration of the grid where each enterprise is located and the proportional reduction amount of each pollutant average concentration of each grid;
and the evaluation module is used for acquiring the reduction contribution rate of each enterprise remediation to the pollutant concentration according to the pollutant concentration reduction amount caused by each enterprise remediation and the proportional reduction amount of the average concentration of all pollutants in the area to be evaluated.
In a third aspect, an embodiment of the present invention provides an electronic device, which includes a memory, a processor, and a computer program stored in the memory and capable of running on the processor, where the processor executes the computer program to implement the steps of the method for evaluating an effect of remediating atmospheric pollution provided by the first aspect of the present invention.
In a fourth aspect, an embodiment of the present invention provides a non-transitory computer-readable storage medium, on which a computer program is stored, where the computer program, when executed by a processor, implements the steps of the method for evaluating the effect of treating air pollution provided by the first aspect of the present invention.
According to the method and the system for evaluating the atmospheric pollution remediation effect, provided by the embodiment of the invention, the area to be evaluated is divided into a plurality of grids, the average concentration of each pollutant in each grid is calculated, the contribution rate of an enterprise to the pollutant concentration of each grid is calculated according to the pollutant emission amount of all enterprises in each grid, and the effect of pollutant concentration reduction caused by enterprise remediation is calculated according to the geometric variation condition of the pollutant concentration in each grid, so that the reduction contribution rate of various preset key industry enterprises to the pollutant concentration of the whole area to be evaluated is obtained.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and those skilled in the art can also obtain other drawings according to the drawings without creative efforts.
Fig. 1 is a flowchart of a method for evaluating an effect of remediation of atmospheric pollution according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of grid and enterprise information partitioning in an embodiment of the present invention;
fig. 3 is a schematic structural diagram of an atmospheric pollution remediation effect evaluation system according to an embodiment of the present invention;
fig. 4 is a schematic physical structure diagram of an electronic device according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Fig. 1 is a flowchart of an atmospheric pollution remediation effect evaluation method according to an embodiment of the present invention, and as shown in fig. 1, an atmospheric pollution remediation effect evaluation method according to an embodiment of the present invention includes:
s1, dividing a region to be evaluated into a plurality of grids, obtaining the average concentration of each pollutant in the region to be evaluated and the average concentration of each pollutant in each grid, and obtaining the proportional reduction of the average concentration of all pollutants in the region to be evaluated and the proportional reduction of the average concentration of each pollutant in each grid;
s2, acquiring the pollutant emission amount of each enterprise in all preset key industries in the area to be evaluated;
s3, acquiring the contribution rate of each enterprise to each pollutant concentration of the grid according to the emission amount of each pollutant of each enterprise and the average concentration of each pollutant of each grid;
s4, calculating the pollutant concentration reduction amount caused by the renovation of each enterprise according to the proportional reduction amount of each pollutant concentration contribution rate of each enterprise to the grid and each pollutant average concentration of each grid;
and S5, acquiring the reduction contribution rate of the remediation of each preset key industry to the pollutant concentration according to the preset key industry corresponding to each enterprise, the pollutant concentration reduction amount caused by the remediation of each enterprise and the proportional reduction amount of the average concentration of all pollutants in the area to be assessed.
Firstly, the average concentration of each pollutant in an area to be evaluated is obtained, in the embodiment of the invention, the change of the average concentration of the atmospheric pollutants in the area to be evaluated is mainly detected, and the atmospheric pollutants comprise PM2.5、NO2And SO2It is also the specific meaning of each pollutant referred to in the embodiments of the present invention, which can be adjusted according to actual needs, and therefore, the PM in the area to be evaluated is obtained2.5Average concentration of (3), NO2Average concentration of (3) and SO2Is measured.
FIG. 2 is a schematic diagram of grid and enterprise information division in an embodiment of the present invention, and as shown in FIG. 2, an area to be evaluated is divided into a plurality of regular grids, and a PM in each grid is obtained2.5Average concentration of (3), NO2Average concentration of (3) and SO2The average concentration of the PM in the area to be evaluated needs to be obtained, and the proportional decrease of the average concentration of all the pollutants in the area to be evaluated is the decrease of the average concentration of the three pollutants in the current year and the average concentration of the three pollutants in the previous year2.5Decrease of mean concentration of (1), NO2Is proportionally decreased amount of SO and2and then obtaining the same-ratio reduction of the average concentration of all the pollutants.
It is also necessary to obtain a proportional reduction of the average concentration of each pollutant for each grid, where the proportional reduction of the average concentration of each pollutant for each grid refers to the current year PM in each grid2.5To the average concentration of PM in the last year2.5The amount of decrease in average concentration of, the current year NO in each grid2To the last year NO2The decrease in average concentration of, the current annual SO in each grid2Is compared with the last year SO2Average of (2)The amount of decrease in concentration.
And then acquiring each pollutant emission amount of each enterprise in all preset key industries in the area to be evaluated, wherein the preset key industries in the embodiment of the invention comprise an electric power heating industry, a metallurgical industry, a building material industry and a petrochemical industry, and the preset key industries can be specifically adjusted according to actual conditions to acquire all enterprises in the industries and acquire PM of each enterprise2.5Emission amount of NO2Emission and SO2And (4) discharging the amount.
In the embodiment of the invention, the PM of each enterprise is obtained by acquiring the inventory data of the enterprise2.5Emission amount of NO2Emission and SO2The discharge amount and the list data comprise the geographical position and discharge amount information of the enterprise, including longitude and latitude coordinates and SO2Emission amount of NOxEmission amount of primary particulate matter PM2.5Emissions, BC and OC emissions, wherein enterprise PM2.5The emission of primary particulate matter PM2.5The emission, BC and OC emissions add up, i.e. ei,j=ei,j(PM2.5)+ei,j(BC)+ei,j(OC),NO2The discharge amount is NOxAnd (4) discharging the amount.
Then according to PM of each enterprise2.5Emission amount of NO2Emission and SO2Emission and PM per grid2.5Average concentration of (3), NO2Average concentration of (3) and SO2The average concentration of the network nodes is obtained, the contribution rate of each enterprise to the concentration of each pollutant of the grid where the enterprise is located is obtained, specifically, any enterprise is taken as an example for explanation, the grid where the enterprise is located is obtained according to the longitude and latitude coordinates of the enterprise, and then the PM of the enterprise is used2.5Emission divided by PM in the grid2.5The average concentration of the PM in the grid of the enterprise can be obtained2.5Contribution rate of concentration, NO of the enterprise2The emission divided by NO in the grid2The average concentration of the NO in the grid of the enterprise can be obtained2Contribution rate of concentration, SO of the enterprise2The emission is divided by the SO in the grid2The average concentration of the SO in the grid of the enterprise can be obtained2Concentration of tributeThe contribution rate.
Obtaining the pollutant concentration reduction amount caused by the renovation of each enterprise according to the contribution rate of each pollutant concentration of the grid where each enterprise is located and the proportional reduction amount of each pollutant average concentration in each grid, taking any enterprise as an example for explanation, and enabling the enterprise to carry out PM (particulate matter) reduction on the grid2.5Concentration contribution multiplied by PM in the grid2.5The same proportion reduction amount of the average concentration can obtain PM brought by the enterprise remediation2.5The average concentration is decreased, and the enterprise is applied to NO in the grid2The contribution of the concentration multiplied by the NO in the grid2The same proportion of the average concentration is decreased, and NO brought by the enterprise remediation can be obtained2The average concentration is decreased, and the enterprise is applied to the SO in the grid2The contribution of the concentration multiplied by the SO in the grid2The same proportion of the average concentration is decreased, SO brought by the enterprise remediation can be obtained2The amount of decrease in average concentration.
And finally, obtaining the reduction contribution rate of each pollutant concentration caused by the renovation of each industry according to the preset key industry corresponding to each enterprise, the pollutant concentration reduction amount caused by the renovation of each enterprise and the proportional reduction amount of the average concentration of all pollutants in the area to be evaluated. Specifically, taking the electric heating industry as an example for explanation, all enterprises in the industry are acquired, and then the PM brought by each enterprise is regulated2.5The sum of the decrease amounts of the average concentration divided by the PM in the region to be evaluated2.5The PM can be regulated in the power heating industry according to the proportional reduction of the average concentration2.5The effect of the decrease in average concentration. NO brought about by every enterprise2The sum of the decrease in the average concentration divided by the NO in the region to be evaluated2The proportional decrease of the average concentration can obtain the NO regulation of the power heat supply industry2The effect of the decrease in average concentration. SO brought by each enterprise renovation2The sum of the drops in the average concentration divided by the SO in the area to be evaluated2The proportional decrease of the average concentration can obtain the SO regulation in the power heat supply industry2The effect of the decrease in average concentration.
Table 1 is an assessment table of the remediation effect of enterprises in the preset key industry in the area to be assessed in the embodiment of the present invention, as shown in table 1:
TABLE 1
Figure BDA0002347011260000091
According to the method for evaluating the atmospheric pollution remediation effect, the area to be evaluated is divided into a plurality of grids, the average concentration of each pollutant in each grid is calculated, the contribution rate of an enterprise to the pollutant concentration of each grid is calculated according to the pollutant emission amount of all enterprises in each grid, the effect of pollutant concentration reduction caused by enterprise remediation is calculated according to the same-proportion change condition of the pollutant concentration in each grid, and therefore the reduction contribution rate of various preset key industry enterprises to the pollutant concentration of the whole area to be evaluated is obtained.
On the basis of the foregoing embodiment, preferably, the dividing the region to be evaluated into a plurality of grids, obtaining an average concentration of each pollutant in the region to be evaluated and an average concentration of each pollutant in each grid, and obtaining a proportional decrease amount of the average concentration of all pollutants in the region to be evaluated and a proportional decrease amount of the average concentration of each pollutant in each grid specifically includes:
acquiring the concentration of each pollutant at each pixel point in the area to be evaluated in the current year and the concentration of each pollutant at each pixel point in the area to be evaluated in the previous year by utilizing a satellite remote sensing inversion technology, and further acquiring the average concentration of each pollutant in the area to be evaluated in the current year and the concentration of each pollutant at each pixel point in the area to be evaluated in the previous year;
dividing the area to be evaluated into a plurality of grids, obtaining the average concentration of each pollutant at each grid in the current year according to the concentration of each pollutant at each pixel point in the area to be evaluated in the current year, and obtaining the average concentration of each pollutant at each grid in the previous year according to the concentration of each pollutant at each pixel point in the area to be evaluated in the previous year;
obtaining the proportional reduction of the average concentration of all pollutants in the area to be evaluated according to the average concentration of each pollutant in the area to be evaluated in the current year and the average concentration of each pollutant in the area to be evaluated in the previous year;
and obtaining the proportional reduction of the average concentration of each pollutant of each grid according to the average concentration of each pollutant of each grid in the current year and the average concentration of each pollutant of each grid in the previous year.
Firstly, a satellite remote sensing technology is utilized to obtain PM at each pixel point in an area to be evaluated in the current year2.5Concentration, NO2Concentration and SO2Concentration, PM at each pixel point in the region to be evaluated in the previous year2.5Concentration, NO2Concentration and SO2Concentration according to PM at each pixel point in the current year area to be evaluated2.5Concentration, NO2Concentration and SO2Concentration, obtaining PM in the current year area to be evaluated2.5Average concentration, NO2Average concentration and SO2Average concentration according to PM at each pixel point in the region to be evaluated in the previous year2.5Concentration, NO2Concentration and SO2Concentration, obtaining PM in the area to be evaluated in the last year2.5Average concentration, NO2Average concentration and SO2Average concentration.
The PM in the area to be evaluated in the last year is processed2.5Average concentration, NO2Average concentration and SO2Average concentration minus PM in the current year area to be evaluated2.5Average concentration, NO2Average concentration and SO2The average concentration, the PM in the area to be evaluated can be obtained2.5Decrease of mean concentration by same ratio, NO2Decrease in mean concentration by same ratio and SO2The decrease in average concentration was proportional.
Then, the area to be evaluated is divided into a plurality of regular grids, and the PM of each pixel point in the grids is calculated according to the current year2.5Concentration, NO2Concentration and SO2Concentration, obtaining the current yearPM in grid2.5Average concentration, NO2Average concentration and SO2Average concentration according to PM of each pixel point in the grid in the last year2.5Concentration, NO2Concentration and SO2Concentration, obtaining PM in the grid in last year2.5Average concentration, NO2Average concentration and SO2Average concentration.
For any grid, PM in the grid in the last year2.5Average concentration, NO2Average concentration and SO2Average concentration minus PM in the grid at the current year2.5Average concentration, NO2Average concentration and SO2Average concentration, PM in the grid can be obtained2.5Decrease of mean concentration by same ratio, NO2Decrease in mean concentration by same ratio and SO2The decrease in average concentration was proportional.
Specifically, based on MODIS data, a geographical weighting model is utilized to obtain near-ground PM of each pixel point in the current year to-be-evaluated area2.5Concentration and near-ground PM of each pixel point in last year region to be evaluated2.5And (4) concentration.
Based on TropomI data, a difference absorption method is utilized to obtain troposphere NO of each pixel point in the region to be evaluated in the current year2The column concentration and troposphere NO of each pixel point in the last year region to be evaluated2Column concentration.
Based on TroPOMI data, the SO of each pixel point in the current year to-be-evaluated area is obtained by using an optimal estimation method2Concentration and SO of each pixel point in last year region to be evaluated2And (4) concentration.
On the basis of the foregoing embodiment, preferably, the obtaining a contribution rate of each enterprise to each pollutant concentration of the grid where each enterprise is located according to each pollutant emission amount of each enterprise and each pollutant average concentration of each grid specifically includes:
for any grid, acquiring each pollutant emission of each enterprise in the grid according to the following formula:
Figure BDA0002347011260000111
wherein E isjRepresenting the average concentration of the jth pollutant, m representing the number of enterprises in all preset key industries of any grid, en,jRepresenting the discharge amount of j pollutant of n enterprise in any grid;
obtaining the contribution rate of each enterprise to each pollutant concentration of any grid according to the emission amount of each pollutant in any grid and the emission amount of each pollutant of each enterprise in any grid
Taking any grid as an example for illustration, the PM of each enterprise in the grid is used2.5Adding the emission to obtain the PM of all enterprises in the grid2.5Emission of NO from each enterprise in the grid2Adding the discharge amount to obtain NO of all enterprises in the grid2Emission of SO for each enterprise in the grid2Adding the emission to obtain SO of all enterprises in the grid2And (4) discharging the amount.
Then PM of all enterprises in the grid2.5Emission divided by PM of the grid2.5The emission can obtain the PM of the enterprises in the grid2.5Contribution rate of concentration, NO of all enterprises in the grid2The emission divided by the NO of the grid2The emission can obtain the NO of the enterprise in the grid2Contribution rate of concentration, SO of all enterprises in the grid2The emission divided by the SO of the grid2The emission can obtain the SO of the enterprise in the grid2Contribution rate of concentration.
On the basis of the above embodiment, preferably, the method obtains each pollutant concentration contribution rate of each enterprise to any grid according to the emission amount of each pollutant in any grid and each pollutant emission amount of each enterprise in any grid, and for any grid, the concrete calculation formula of each pollutant concentration contribution rate of enterprises in any grid to any grid is as follows:
Figure BDA0002347011260000121
wherein n isn,jRepresenting the contribution rate of the nth enterprise in any grid to the j pollutant emission amount of any grid, en,jRepresents the j pollutant discharge amount of the nth enterprise, EjRepresents the total emission amount of the j pollutant of any grid.
On the basis of the foregoing embodiment, preferably, the calculating, according to a percentage contribution of each enterprise to each pollutant concentration of the grid where the enterprise is located and a proportional decrease of each pollutant average concentration of each grid, a pollutant concentration decrease caused by remediation of each enterprise specifically includes:
Figure BDA0002347011260000122
wherein, DeltasumIndicates the amount of decrease in contaminant concentration, D, caused by each enterprise remediationi,j,1Representing the concentration of the j-th pollutant in the ith grid of the previous year, Di,j,2Represents the concentration of the j-th pollutant in the ith grid of the current year, ni,jAnd representing the contribution rate of the ith enterprise in any grid to the j pollutant emission amount of any grid.
In the above formula, the first and second carbon atoms are,
Figure BDA0002347011260000123
can be regarded as the proportional drop of the j pollutant.
On the basis of the above embodiment, preferably, the proportional decrease of the average concentration of all the pollutants in the area to be evaluated is obtained by the following formula:
Figure BDA0002347011260000124
wherein, delta'sumRepresenting the same proportional decrease, D ', of the average concentration of all contaminants in the area to be assessed'i,j,1Representing the j < th > pollutant in the ith grid of the previous yearConcentration of (D)'i,j,2Indicating the concentration of the jth pollutant in the ith grid of the current year.
Fig. 3 is a schematic structural diagram of an atmospheric pollution abatement effect evaluation system according to an embodiment of the present invention, and as shown in fig. 3, the atmospheric pollution abatement effect evaluation system includes: a year-on-year module 301, an enterprise discharge module 302, a contribution rate module 303, an remediation module 304, and an assessment module 305. Wherein:
the same-proportion module 301 is configured to divide the area to be evaluated into a plurality of grids, obtain an average concentration of each pollutant in the area to be evaluated and an average concentration of each pollutant in each grid, and obtain a same-proportion decrease amount of the average concentrations of all pollutants in the area to be evaluated and a same-proportion decrease amount of the average concentrations of each pollutant in each grid;
the enterprise emission module 302 is configured to obtain an emission amount of each pollutant of each enterprise in all preset key industries in the area to be evaluated;
the contribution rate module 303 is configured to obtain a contribution rate of each enterprise to each pollutant concentration of the grid in which each enterprise is located according to each pollutant emission amount of each enterprise and each pollutant average concentration of each grid;
the remediation module 304 is configured to calculate a pollutant concentration reduction amount caused by remediation of each enterprise according to a proportional reduction amount of each pollutant concentration contribution rate of each enterprise to the grid where the enterprise is located and each pollutant average concentration of each grid;
the evaluation module 305 is configured to obtain a decreasing contribution rate of each enterprise remediation to the pollutant concentration according to the pollutant concentration decreasing amount caused by each enterprise remediation and the proportional decreasing amount of the average concentration of all pollutants in the area to be assessed.
The system embodiment provided in the embodiments of the present invention is for implementing the above method embodiments, and for details of the process and the details, reference is made to the above method embodiments, which are not described herein again.
Fig. 4 is a schematic entity structure diagram of an electronic device according to an embodiment of the present invention, and as shown in fig. 4, the electronic device may include: a processor (processor)401, a communication Interface (communication Interface)402, a memory (memory)403 and a bus 404, wherein the processor 401, the communication Interface 402 and the memory 403 complete communication with each other through the bus 404. The communication interface 402 may be used for information transfer of an electronic device. Processor 401 may call logic instructions in memory 403 to perform a method comprising:
dividing a region to be evaluated into a plurality of grids, obtaining the average concentration of each pollutant in the region to be evaluated and the average concentration of each pollutant in each grid, and obtaining the proportional reduction of the average concentration of all pollutants in the region to be evaluated and the proportional reduction of the average concentration of each pollutant in each grid;
acquiring the pollutant emission amount of each enterprise in all preset key industries in the area to be evaluated;
acquiring the contribution rate of each enterprise to the concentration of each pollutant of the grid according to the emission amount of each pollutant of each enterprise and the average concentration of each pollutant of each grid;
calculating pollutant concentration reduction caused by the renovation of each enterprise according to the proportional reduction of each pollutant concentration contribution rate of each enterprise to the grid and each pollutant average concentration of each grid;
and acquiring the contribution rate of the remediation of each preset key industry to the reduction of the pollutant concentration according to the preset key industry corresponding to each enterprise, the pollutant concentration reduction amount caused by the remediation of each enterprise and the proportional reduction amount of the average concentration of all pollutants in the area to be assessed.
In addition, the logic instructions in the memory 403 may be implemented in the form of software functional units and stored in a computer readable storage medium when the software functional units are sold or used as independent products. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the above-described method embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
In another aspect, an embodiment of the present invention further provides a non-transitory computer-readable storage medium, on which a computer program is stored, where the computer program is implemented to perform the transmission method provided in the foregoing embodiments when executed by a processor, and for example, the method includes:
dividing a region to be evaluated into a plurality of grids, obtaining the average concentration of each pollutant in the region to be evaluated and the average concentration of each pollutant in each grid, and obtaining the proportional reduction of the average concentration of all pollutants in the region to be evaluated and the proportional reduction of the average concentration of each pollutant in each grid;
acquiring the pollutant emission amount of each enterprise in all preset key industries in the area to be evaluated;
acquiring the contribution rate of each enterprise to the concentration of each pollutant of the grid according to the emission amount of each pollutant of each enterprise and the average concentration of each pollutant of each grid;
calculating pollutant concentration reduction caused by the renovation of each enterprise according to the proportional reduction of each pollutant concentration contribution rate of each enterprise to the grid and each pollutant average concentration of each grid;
and acquiring the contribution rate of the remediation of each preset key industry to the reduction of the pollutant concentration according to the preset key industry corresponding to each enterprise, the pollutant concentration reduction amount caused by the remediation of each enterprise and the proportional reduction amount of the average concentration of all pollutants in the area to be assessed.
The above-described embodiments of the apparatus are merely illustrative, and the units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. One of ordinary skill in the art can understand and implement it without inventive effort.
Through the above description of the embodiments, those skilled in the art will clearly understand that each embodiment can be implemented by software plus a necessary general hardware platform, and certainly can also be implemented by hardware. With this understanding in mind, the above-described technical solutions may be embodied in the form of a software product, which can be stored in a computer-readable storage medium such as ROM/RAM, magnetic disk, optical disk, etc., and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the methods described in the embodiments or some parts of the embodiments.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.

Claims (10)

1. An evaluation method for the remediation effect of atmospheric pollution, comprising:
dividing a region to be evaluated into a plurality of grids, obtaining the average concentration of each pollutant in the region to be evaluated and the average concentration of each pollutant in each grid, and obtaining the proportional reduction of the average concentration of all pollutants in the region to be evaluated and the proportional reduction of the average concentration of each pollutant in each grid;
acquiring the pollutant emission amount of each enterprise in all preset key industries in the area to be evaluated;
acquiring the contribution rate of each enterprise to the concentration of each pollutant of the grid according to the emission amount of each pollutant of each enterprise and the average concentration of each pollutant of each grid;
calculating pollutant concentration reduction caused by the renovation of each enterprise according to the proportional reduction of each pollutant concentration contribution rate of each enterprise to the grid and each pollutant average concentration of each grid;
and acquiring the contribution rate of the remediation of each preset key industry to the reduction of the pollutant concentration according to the preset key industry corresponding to each enterprise, the pollutant concentration reduction amount caused by the remediation of each enterprise and the proportional reduction amount of the average concentration of all pollutants in the area to be assessed.
2. The method according to claim 1, wherein the step of dividing the area to be evaluated into a plurality of grids, obtaining the average concentration of each pollutant in the area to be evaluated and the average concentration of each pollutant in each grid, and obtaining the proportional reduction of the average concentration of all pollutants in the area to be evaluated and the proportional reduction of the average concentration of each pollutant in each grid specifically comprises:
acquiring the concentration of each pollutant at each pixel point in the area to be evaluated in the current year and the concentration of each pollutant at each pixel point in the area to be evaluated in the previous year by utilizing a satellite remote sensing inversion technology, and further acquiring the average concentration of each pollutant in the area to be evaluated in the current year and the concentration of each pollutant at each pixel point in the area to be evaluated in the previous year;
dividing the area to be evaluated into a plurality of grids, obtaining the average concentration of each pollutant at each grid in the current year according to the concentration of each pollutant at each pixel point in the area to be evaluated in the current year, and obtaining the average concentration of each pollutant at each grid in the previous year according to the concentration of each pollutant at each pixel point in the area to be evaluated in the previous year;
obtaining the proportional reduction of the average concentration of all pollutants in the area to be evaluated according to the average concentration of each pollutant in the area to be evaluated in the current year and the average concentration of each pollutant in the area to be evaluated in the previous year;
and obtaining the proportional reduction of the average concentration of each pollutant of each grid according to the average concentration of each pollutant of each grid in the current year and the average concentration of each pollutant of each grid in the previous year.
3. The method according to claim 2, wherein the step of obtaining the contribution rate of each enterprise to each pollutant concentration of the grid on which the enterprise is located according to the emission amount of each pollutant of each enterprise and the average concentration of each pollutant of each grid specifically comprises:
for any grid, acquiring each pollutant emission of each enterprise in the grid according to the following formula:
Figure FDA0002347011250000021
wherein E isjRepresenting the average concentration of the jth pollutant, m representing the number of enterprises in all preset key industries of any grid, en,jRepresenting the discharge amount of j pollutant of n enterprise in any grid;
and acquiring the contribution rate of each enterprise to each pollutant concentration of any grid according to the emission rate of each pollutant in any grid and the emission rate of each pollutant of each enterprise in any grid.
4. The method according to claim 3, wherein the specific calculation formula of each pollutant concentration contribution rate of each enterprise on any grid is obtained according to the emission amount of each pollutant in any grid and the emission amount of each pollutant of each enterprise in any grid, and for any grid, the specific calculation formula of each pollutant concentration contribution rate of each enterprise in any grid on any grid is as follows:
Figure FDA0002347011250000022
wherein n isn,jRepresenting the contribution rate of the nth enterprise in any grid to the j pollutant emission amount of any grid, en,jRepresents the j pollutant discharge amount of the nth enterprise, EjRepresents the total emission amount of the j pollutant of any grid.
5. The method according to claim 1, wherein the step of calculating the pollutant concentration reduction amount caused by the remediation of each enterprise according to the proportional reduction amount of each pollutant concentration contribution rate of each enterprise to the grid and each pollutant average concentration of each grid comprises:
Figure FDA0002347011250000031
wherein, DeltasumIndicates the amount of decrease in contaminant concentration, D, caused by each enterprise remediationi,j,1Representing the concentration of the j-th pollutant in the ith grid of the previous year, Di,j,2Represents the concentration of the j-th pollutant in the ith grid of the current year, ni,jAnd representing the contribution rate of the ith enterprise in any grid to the j pollutant emission amount of any grid.
6. The method for evaluating the atmospheric pollution remediation effect of claim 1, wherein the proportional decrease in the average concentration of all pollutants in the area to be evaluated is obtained by the following formula:
Figure FDA0002347011250000032
wherein, delta'sumRepresenting the area to be evaluatedProportionally reduced amount of average concentration of all contaminants in D'i,j,1Representing the concentration, D ', of the j-th pollutant in the ith grid of the last year'i,j,2Indicating the concentration of the jth pollutant in the ith grid of the current year.
7. The method of claim 1, wherein each pollutant comprises PM2.5、NO2And SO2
8. An atmospheric pollution remediation effect evaluation system, comprising:
the system comprises a comparing module, a calculating module and a judging module, wherein the comparing module is used for dividing a region to be evaluated into a plurality of grids, acquiring the average concentration of each pollutant in the region to be evaluated and the average concentration of each pollutant in each grid, and acquiring the comparing descending amount of the average concentration of all pollutants in the region to be evaluated and the comparing descending amount of the average concentration of each pollutant in each grid;
the enterprise emission module is used for acquiring the emission amount of each pollutant of each enterprise in all preset key industries in the area to be evaluated;
the contribution rate module is used for acquiring the contribution rate of each enterprise to the concentration of each pollutant of the grid according to the emission amount of each pollutant of each enterprise and the average concentration of each pollutant of each grid;
the regulation module is used for calculating the pollutant concentration reduction amount caused by regulation of each enterprise according to the contribution rate of each enterprise to the pollutant concentration of the grid where each enterprise is located and the proportional reduction amount of each pollutant average concentration of each grid;
and the evaluation module is used for acquiring the reduction contribution rate of each enterprise remediation to the pollutant concentration according to the pollutant concentration reduction amount caused by each enterprise remediation and the proportional reduction amount of the average concentration of all pollutants in the area to be evaluated.
9. An electronic device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the steps of the method for assessing the effectiveness of atmospheric pollution remediation according to any one of claims 1 to 7.
10. A non-transitory computer-readable storage medium having stored thereon a computer program, wherein the computer program when executed by a processor implements the steps of the method for assessing the effectiveness of atmospheric pollution remediation of any one of claims 1 to 7.
CN201911398832.2A 2019-12-30 2019-12-30 Atmospheric pollution remediation effect evaluation method and system Pending CN111126872A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911398832.2A CN111126872A (en) 2019-12-30 2019-12-30 Atmospheric pollution remediation effect evaluation method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911398832.2A CN111126872A (en) 2019-12-30 2019-12-30 Atmospheric pollution remediation effect evaluation method and system

Publications (1)

Publication Number Publication Date
CN111126872A true CN111126872A (en) 2020-05-08

Family

ID=70505539

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911398832.2A Pending CN111126872A (en) 2019-12-30 2019-12-30 Atmospheric pollution remediation effect evaluation method and system

Country Status (1)

Country Link
CN (1) CN111126872A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111626624A (en) * 2020-05-29 2020-09-04 四川省环境政策研究与规划院 Air quality improvement evaluation method
CN111695832A (en) * 2020-06-19 2020-09-22 邵峥业 Lead-related industry comprehensive improvement effect evaluation system and method
CN111929263A (en) * 2020-08-11 2020-11-13 生态环境部卫星环境应用中心 Ozone profile and sulfur dioxide column concentration cooperative inversion method

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111626624A (en) * 2020-05-29 2020-09-04 四川省环境政策研究与规划院 Air quality improvement evaluation method
CN111695832A (en) * 2020-06-19 2020-09-22 邵峥业 Lead-related industry comprehensive improvement effect evaluation system and method
CN111695832B (en) * 2020-06-19 2022-07-22 苏州喜全软件科技有限公司 Lead-related industry comprehensive improvement effect evaluation system and method
CN111929263A (en) * 2020-08-11 2020-11-13 生态环境部卫星环境应用中心 Ozone profile and sulfur dioxide column concentration cooperative inversion method
CN111929263B (en) * 2020-08-11 2021-05-11 生态环境部卫星环境应用中心 Ozone profile and sulfur dioxide column concentration cooperative inversion method

Similar Documents

Publication Publication Date Title
CN111126872A (en) Atmospheric pollution remediation effect evaluation method and system
Paatero et al. Methods for estimating uncertainty in factor analytic solutions
Szpiro et al. Efficient measurement error correction with spatially misaligned data
Kumar et al. Forecasting of air quality in Delhi using principal component regression technique
Petrişor Using CORINE data to look at deforestation in Romania: Distribution & possible consequences
Benedetti et al. Tools to support a model-based methodology for emission/immission and benefit/cost/risk analysis of wastewater systems that considers uncertainty
Ashok et al. Development of a response surface model of aviation's air quality impacts in the United States
US20090024407A1 (en) Indicating which of forecasting models at different aggregation levels has a better forecast quality
Jia et al. Estimation and analysis of emissions from on-road vehicles in Mainland China for the period 2011–2015
Brunelle-Yeung et al. Assessing the impact of aviation environmental policies on public health
Hantush et al. Bayesian framework for water quality model uncertainty estimation and risk management
CN110988269A (en) Deviation correction method and device for atmospheric pollution source emission list and storage medium
Sharma et al. Application of extreme value theory for predicting violations of air quality standards for an urban road intersection
Rehana et al. Modeling of extreme risk in river water quality under climate change
Dennis et al. Observable indicators of the sensitivity of PM2. 5 nitrate to emission reductions—Part II: Sensitivity to errors in total ammonia and total nitrate of the CMAQ-predicted non-linear effect of SO2 emission reductions
CN116776073B (en) Pollutant concentration evaluation method and device
Vinod et al. New approach for risk based inspection of H2S based process plants
Zheng et al. Quantification of variability and uncertainty using mixture distributions: evaluation of sample size, mixing weights, and separation between components
Nguyen et al. The influence of model structure uncertainty on water quality assessment
Tokuno et al. Codesign-oriented performability modeling for hardware-software systems
JP5738140B2 (en) Substance release amount estimation apparatus, method and program
CN114705249B (en) Pollutant emission monitoring method based on artificial intelligence and related equipment
Song et al. Development of a multimedia model (POPsLTEA) to assess the influence of climate change on the fate and transport of polycyclic aromatic hydrocarbons in East Asia
Chettouh et al. Industrial fire simulation and uncertainty associated with the emission dispersion model
Wikle et al. A mechanistic model of annual sulfate concentrations in the United States

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20200508