A method, device, storage medium and product for rapid quantification of local and transport contributions of pollutants

Abstract

The present invention relates to a method for quantifying the local and transmission contributions of six air quality pollutants, in particular to a method, device, storage medium and product for rapid quantification of the local and transmission contributions of pollutants. A method for rapid quantification of the local and transmission contributions of pollutants, comprising the following steps: using Meteoinfo software to calculate according to a certain pollutant concentration threshold value to obtain a grid data result; using ArcGIS software to perform grid point assignment calculation according to the obtained grid data result to obtain the absolute proportion of each grid point to the evaluation site; correcting the concentration according to the obtained absolute proportion of each grid point to the evaluation site and outputting the result. A method for rapid quantification of the local and transmission contributions of pollutants provided by the present invention combines Meteoinfo software with ArcGIS software, re-assigns and calculates the grid data of the shp file of PSCF in ArcGIS, can rapidly and quantitatively evaluate the external transmission and local contributions of the point, and has the advantages of small calculation amount, flexibility, simple and easy operation, etc.

Description

Technical Field

[0001] The present invention relates to a method for quantifying the local and transmission contributions of six air quality pollutants, in particular to a method, device, storage medium and product for quickly quantifying the local and transmission contributions of pollutants. Background Art

[0002] Ambient air quality is affected by the emission intensity of local sources on the one hand, and by the impact of meteorology on the transmission, diffusion and deposition of pollutants on the other hand. Therefore, it is particularly important to be able to quickly and quantitatively assess the impact of local and external transmission of air quality pollutants.

[0003] In recent years, air quality models have become one of the important means to study the spatiotemporal distribution characteristics of atmospheric pollutants and predict the quality of the atmospheric environment. Common models include the CAMx model, CMAQ model, NAQPMS model, and WRF-Chem model, which are typical comprehensive regional scale models. They have played a huge role in simulating various complex air quality problems and regional complex air pollution problems such as O3, particulate matter, visibility, acid rain, and even climate change, and can evaluate the impact of local and transmission contributions. However, due to the complexity of atmospheric pollution at the comprehensive regional scale, there are complex physical and chemical reactions and phase transitions between various pollutants, which cannot be simulated by simple small and medium scale models. At the same time, air quality models have certain limitations on the use of models due to problems such as incomplete clarity in the atmospheric chemical reaction mechanism, high hardware requirements, large amount of calculations, and high costs. Summary of the invention

[0004] To solve the above problems, the present invention provides a method for rapid quantification of local and transport contributions of pollutants. The specific technical solution is as follows:

[0005] A method for rapid quantification of local and transport contributions of pollutants, characterized by comprising the following steps:

[0006] Meteoinfo software is used to calculate the concentration threshold of a certain pollutant to obtain gridded data results;

[0007] According to the grid data results, ArcGIS software is used to calculate the grid point assignment to obtain the absolute proportion of each grid point to the evaluation site;

[0008] The concentration is corrected based on the absolute proportion of each grid point to the assessment site and the results are output.

[0009] Preferably, the calculation using Meteoinfo software includes the following steps:

[0010] Determine the concentration threshold of a pollutant to be assessed. After determining the assessment point, select the concentration threshold of a pollutant to be assessed;

[0011] The gridded data results were obtained using Meteoinfo software.

[0012] Furthermore, the determination of the concentration threshold of a pollutant to be evaluated includes the following methods: using a random forest method, PM 2.5 The 10th percentile of pollutant concentration was selected as the threshold, CO, PM 10 , SO2, NO2 and O3 take the average of excellent and low-value good as the threshold, and the low-value good AQI is <75.

[0013] Preferably, the ArcGIS software includes the following steps when performing grid point value calculation:

[0014] Add city or county basemap;

[0015] Add the gridded data results to obtain the PSCF.shp file;

[0016] Change the coordinate system of the PSCF.shp file to a projected coordinate system;

[0017] Set parameters for the PSCF.shp file and calculate the main grid area;

[0018] Identify the elements of the PSCF.shp file;

[0019] Perform sub-grid area calculation on the PSCF.shp file with completed element identification;

[0020] Calculate the proportion of the secondary grid area to the primary grid area;

[0021] Calculate the absolute proportion of each grid Mij;

[0022] Output grid calculation results.

[0023] Further, the adding of the gridded data results includes: adding the shp file of the potential source area obtained by the Meteoinfo software to obtain a PSCF layer;

[0024] The setting of parameters for the PSCF.shp file and the main grid calculation includes: right-clicking the PSCF layer, selecting "Open Attribute Table" - "Add Field" - naming it as "MJ" - selecting double precision in "Type", right-clicking "Calculate Geometry", and changing the unit to "Square Kilometers";

[0025] The element identification of the PSCF.shp file includes: "Analysis Tools" - "Identification", inputting relevant elements, clicking "OK", and a new layer appears, named "PSCFBS";

[0026] The sub-grid area calculation includes: right-clicking the PSCFBS layer, selecting "Open Attribute Table" - "Add Field" - naming it as "BSMJ" - selecting double precision in "Type", and clicking "OK"; right-clicking the "BSMJ" field - "Calculate Geometry", changing the unit to "Square Kilometers", and clicking "OK";

[0027] The calculation of the proportion of the secondary grid area in the primary grid area includes: right-clicking the PSCFBS layer, selecting "Open Attribute Table" - "Add Field" - naming it "ZB" - selecting double precision in "Type", and clicking "OK"; right-clicking the "ZB" field - "Field Calculator", clicking "MJBS / MJ", and clicking "OK";

[0028] The calculation of the absolute proportion of Mij in each grid includes: right-clicking the PSCFBS layer, selecting "Open Attribute Table" - "Add Field" - naming it "Finally" - selecting double precision in "Type", and clicking "OK"; right-clicking the "Finally" field - "Field Calculator", clicking "Mij*ZB", and clicking "OK";

[0029] The outputting of the grid calculation results includes: selecting "Conversion Tools" - "Table to Excel", inputting relevant elements and output location and name, and outputting Excel of the grid calculation results.

[0030] Preferably, the modified concentration includes:

[0031] The concentration threshold of a pollutant under evaluation is used as the local concentration;

[0032] The mean concentration of pollutants exceeding the threshold is used as the base of transmission concentration for other areas;

[0033] Multiply the transmission concentration base by the absolute proportion obtained to obtain the contribution concentration of other regions;

[0034] The contribution ratio of different regions and local areas is calculated by the sum and ratio of local concentration and other regions’ contribution concentration:

[0035] Contribution ratio of a certain region = transmission concentration in this region / (local concentration + transmission concentration in this region)

[0036] A device for a method for rapidly quantifying local and transport contributions of pollutants comprises a memory, a processor and a computer program stored in the memory, wherein the processor executes the computer program to implement the steps of the method for rapidly quantifying local and transport contributions of pollutants.

[0037] A computer-readable storage medium stores a computer program / instruction thereon, which, when executed by a processor, implements the steps of the above-mentioned method for rapid quantification of local and transport contributions of pollutants.

[0038] A computer program product comprises a computer program / instruction, which, when executed by a processor, implements the steps of the above-mentioned method for rapid quantification of local and transport contributions of pollutants.

[0039] Compared with the prior art, the present invention has the following beneficial effects:

[0040] The invention provides a method for rapid quantification of local and transport contributions of pollutants. The method combines Meteoinfo software with ArcGIS software, re-assigns and calculates the grid data of the shp file of PSCF in ArcGIS, and can rapidly and quantitatively evaluate the external transmission and local contribution of a point, thus overcoming the shortcomings of the air quality model, such as long time consumption and high cost, and the semi-quantitative results of Meteoinfo software. BRIEF DESCRIPTION OF THE DRAWINGS

[0041] Figure 1 is a flow chart of this application;

[0042] Figure 2 It is a schematic diagram of the PSCF calculation of the pollutant trajectory number in Meteoinfo software;

[0043] Figure 3 It is a schematic diagram of sub-grid area calculation. DETAILED DESCRIPTION

[0044] The present invention will be further described with reference to the accompanying drawings.

[0045] The purpose of the rapid quantification method for the transport contribution of six air quality pollutants is to evaluate the transport contribution of six pollutants (SO2, NO2, CO, O3-8h, PM 10 、PM 2.5 The realization of the rapid quantification method of the transmission contribution of six air quality pollutants can not only qualitatively evaluate the local and external transmission impacts, but also quickly and easily evaluate the quantitative transmission contribution ratios of local and external cities in the pollution process, month, quarter and year, and quickly realize the demand for evaluation results and timely understand and grasp the external transmission impact ratio.

[0046] This method is based on the combination of Meteoinfo software and ArcGIS software to quickly and quantitatively evaluate the transmission contribution of air quality parameters.

[0047] After using Meteoinfo software to obtain the potential source area results, this method loads the gridded data results into ArcGIS, and uses ArcGIS software to perform grid value calculations to obtain the absolute proportion of each grid point to the assessment site. On the one hand, this method can reduce labor and time costs. On the other hand, this method has a small amount of calculation, is simple to learn and easy to operate, and can quickly obtain quantitative assessment results. It solves the problems of existing air quality model assessments such as CAMx model, CMAQ model, WRF-Chem model, etc., which require high hardware, long time, and high cost due to the consideration of meteorological fields, source emission models, and photochemical models in the entire assessment process.

[0048] Meteoinfo software is based on the meteorological data of the Global Data Assimilation System (GDAS) of the National Oceanic and Atmospheric Administration (NOAA) of the United States and the grid statistical analysis method based on the backward trajectory model. It uses the potential source contribution factor method (PSCF) and the concentration weighted trajectory analysis method (CWT) to only semi-quantitatively obtain the distribution of pollution sources. In addition, Meteoinfo software can quickly and conveniently obtain the distribution of pollution sources semi-quantitatively through cluster analysis, potential source areas, etc., but it cannot evaluate the absolute proportion of pollution source contribution points to the evaluation points. However, by combining with ArcGIS software, the grid can be assigned values ​​and subsequent calculations can be performed to achieve quantitative result evaluation, which has the advantages of small calculation amount, flexibility, and simple and easy operation.

[0049] PSCF was calculated by Meteoinfo software in the early stage. According to the software principle, PSCF needs to set a threshold for the pollution factor. When the element value corresponding to the point on the backward trajectory is higher than this threshold, the trajectory point is considered to be a pollution trajectory. Mij is the number of points on the pollution trajectory passing through the grid (i, j) in the selected study area, that is, the path and area where the trajectory exceeding the threshold is mainly concentrated. The larger the value, the greater the contribution of the pollution source in the area. The shp file of PSCF calculated by Meteoinfo software is imported into GIS, the Mij value in the PSCF file is assigned to each grid, and finally the absolute proportion of the area is calculated.

[0050] like Figure 1 A method for rapid quantification of local and transport contributions of pollutants is shown, comprising the following steps:

[0051] According to the estimated concentration threshold of a pollutant, the Meteoinfo software is used to calculate and obtain the grid data results;

[0052] According to the grid data results obtained, ArcGIS software is used to calculate the grid point assignment to obtain the absolute proportion of each grid point to the assessment site; the absolute proportion refers to the ratio of the grid point value to the sum of all grid points, which is a one-to-one grid calculation, and the proportion of the sub-grids 1.1 and 1.2 to the main grid 1 is calculated. For example: after the identification is completed, if the area of ​​the main grid 1 is divided into two parts, the main grid 1 will have two rows, one corresponding to 1.1 and the other corresponding to 1.2. Both rows 1.1 and 1.2 correspond to the area of ​​the main grid 1, and only one-to-one grid calculation is required.

[0053] The concentration is corrected based on the absolute proportion of each grid point to the assessment site and the results are output.

[0054] Specifically, a method for rapid quantification of local and transport contributions of pollutants is divided into three parts to achieve quantitative transport assessment:

[0055] 1. Get the PSCF result shp file through Meteoinfo software:

[0056] After determining the assessment point, select the concentration threshold of a pollutant to be assessed, and use Meteoinfo software to calculate the shp file of PSCF.

[0057] like Figure 2 As shown, Meteoinfo software mainly calculates Mij, which represents the number of pollutant trajectories (Mij) passing through a certain grid (i, j) in the selected study area. This result is one of the contents of the PSCF calculation results performed by Meteoinfo software.

[0058] Determination of the concentration threshold of a certain pollutant to be evaluated: PM 2.5 The 10th percentile of pollutant concentration was selected as the threshold, CO, PM 10 , SO2, NO2 and O3 use the average of excellent and low-quality (AQI < 75) as the threshold.

[0059] The average of excellent days (AQI less than 50) and low-value good days (AQI greater than 50 and less than or equal to 75) in the evaluation period is selected as the threshold, which is to calculate the average of the data of excellent days and low-value good days in the corresponding period.

[0060] 2. Grid computing through ArcGIS software:

[0061] (1) Add a city or county base map;

[0062] (2) Add the shp file of the potential source area obtained by the meteoinfo software, taking the "PSCF.shp" layer as an example;

[0063] (3) Calculation of main grid area: First, change the coordinate system to the projected coordinate system, right-click the PSCF layer, select "Open Attribute Table" - "Add Field" - Name it "MJ" - Select "Type" as double precision. Right-click "Calculate Geometry" and change the unit to "Square Kilometers";

[0064] The main grid area calculation is the grid calculation of the PSCF.shp file, which can also be understood as the Mij grid area (because the PSCF.shp file includes Mij results).

[0065] (4) Element identification: "Analysis Tools" - "Identification", enter the relevant elements, click "OK", a new layer will appear, named "PSCFBS";

[0066] Among them, step (1) and step (2) are two independent layers loaded and do not affect each other. Step (1) contains city information, and step (2) contains grid data results. These two layers are combined into one for subsequent calculations. Step (4) assigns the content of the city or county base map added in step (1) to the PSCF.shp layer of step (2), thereby forming a new layer. The boundary line of the city or county base map in the new layer coincides with the grid of Mij, and the grid of Mij is re-divided. The grid of Mij is divided into smaller grids by the boundary line to obtain a number of sub-grids.

[0067] (5) Subgrid area calculation: Right-click the PSCFBS layer, select "Open Attribute Table" - "Add Field" - Name it "BSMJ" - Select "Type" as double precision, and click "OK". Right-click the "BSMJ" field - "Calculate Geometry", change the unit to "Square Kilometers", and click "OK";

[0068] like Figure 3 As shown, step (3) is to calculate each grid as a whole, step (4) is to divide each grid in step (3) into multiple smaller grids according to the boundary of the base map after assignment, and step (5) is to calculate the area of ​​the smaller grid after the division;

[0069] (6) Calculate the ratio of the marked grid area to the total area: right-click the PSCFBS layer, select "Open Attribute Table" - "Add Field" - Name it "ZB" - Select "Type" as double precision, and click "OK". Right-click the "ZB" field - "Field Calculator", click "BSMJ / MJ", and click "OK";

[0070] The ratio of the grid area after the calculation mark to the total area refers to the ratio of 1.1 and 1.2 to 1 after the calculation is divided: ZB = MJBS / MJ;

[0071] (7) Calculate the absolute proportion of Mij in each grid. Right-click the PSCFBS layer, select "Open Attribute Table" - "Add Field" - Name it "Finally" - Select "Type" as double precision, and click "OK". Right-click the "Finally" field - "Field Calculator", click "Mij*ZB", and click "OK".

[0072] Among them, absolute proportion = Mij×ZB.

[0073] (8) Convert table to Excel. Select "Conversion Tools" - "Convert Table to Excel", enter the relevant elements and output location and name, and output the Excel of the grid calculation results.

[0074] 3. Result correction:

[0075] The concentration threshold of a certain pollutant under assessment is taken as the local concentration; the average concentration of pollutants exceeding the threshold is taken as the transmission concentration base for other regions. The absolute proportion (Finally) result is multiplied by the transmission concentration base to obtain the contribution concentration of other regions. The contribution proportion of different regions and the local region is calculated by adding the local concentration and the contribution concentration of other regions and calculating the proportion.

[0076] The concentration threshold of a pollutant under evaluation is used as the local concentration, i.e. the SO2, NO2, CO, O3-8h, PM 10 、PM 2.5 The threshold value of any one of the six parameters.

[0077] The contribution ratio of a certain region = transmission concentration in this region / (local concentration + transmission concentration in this region).

[0078] Meteoinfo software (model) is used to cluster the trajectories. The principle is to use the Euclidean distance to represent the distance between two similar trajectories, and then classify and group the backward trajectories according to the similarity of the trajectory space. For example, if there are n trajectories, at the beginning, all the trajectories are divided into n categories, the two closest trajectories are selected as the same category, and the average latitude of the two trajectory nodes within one hour is calculated to obtain the average trajectory track. Then, cluster diagrams are established based on these two results. Finally, the optimal path is determined by analyzing the relationship between the coordinates of the cluster points and the maximum acceleration at each moment.

[0079] The potential source contribution function (PSCF) is also known as the residence time analysis method. This method is a method for identifying source areas based on conditional probability functions and airflow trajectory analysis. The PSCF value is determined mainly based on the ratio of the number of pollutant trajectories (Mij) passing through a certain grid (i, j) in the selected study area to the number of all trajectories (Nij) passing through this grid. The value indicates the probability that the air mass passing through this grid affects the air quality of the study point. A certain pollution factor (such as the concentration of pollutants, etc.) is added to the trajectory obtained during the study observation period, and a corresponding threshold is set for the pollution factor. When the factor value corresponding to a point on the backward trajectory is higher than the set threshold, the point is regarded as a pollution trajectory, thereby determining the possible emission source area of ​​the pollutant. The corresponding formula is as follows:

[0080] PSCF ij =M ij / N ij ;

[0081] Since PSCF is a method of analysis using conditional probability, a weight factor is introduced to reduce the error. The PSCF value is multiplied by a weight function (W ij ),Right now:

[0082] WPSCF=W ij ×PSCF;

[0083] Although the PSCF method can give the location of the potential source area of ​​pollutants, it cannot quantitatively give the extent to which the potential source area contributes to the actual pollutants at the receiving point. Therefore, the pollution weight index passing through each grid trajectory can be calculated by the concentration weight trajectory analysis method (CWT) to reflect the contribution of the grid to the pollutant concentration at the receiving point. The calculation of CWT is shown in the following formula. When calculating the CWT value, the same weight function Wij as the PSCF method is used for weighted calculation, that is, WCWT=CWT×W ij value.

[0084] ;

[0085] In the formula, CWT ij is the average weighted concentration of the contaminant in the grid ij, C l It's a track l The pollutant concentration corresponding to the grid ij, τ ijl It's a track l The residence time on the grid ij.

[0086] Based on the calculation results of Meteoinfo software, ArcGIS software is used to perform grid point assignment calculations, thereby achieving absolute quantification of the transmission contribution results of the assessment points.

[0087] A device for a method for rapidly quantifying local and transport contributions of pollutants comprises a memory, a processor and a computer program stored in the memory, wherein the processor executes the computer program to implement the steps of the method for rapidly quantifying local and transport contributions of pollutants.

[0088] A computer-readable storage medium stores a computer program / instruction thereon, which, when executed by a processor, implements the steps of the above-mentioned method for rapid quantification of local and transport contributions of pollutants.

[0089] A computer program product comprises a computer program / instruction, which, when executed by a processor, implements the steps of the above-mentioned method for rapid quantification of local and transport contributions of pollutants.

[0090] The memory and the processor are electrically connected directly or indirectly to achieve data transmission or interaction. For example, these elements can be electrically connected to each other through one or more communication buses or signal lines, such as through a bus connection. The memory stores computer-executable instructions for implementing the data access control method, including at least one software function module that can be stored in the memory in the form of software or firmware. The processor executes various functional applications and data processing by running the software programs and modules stored in the memory.

[0091] The memory may be, but is not limited to, a random access memory (RAM), a read only memory (ROM), a programmable read-only memory (PROM), an erasable programmable read-only memory (EPROM), an electrically erasable read-only memory (EEPROM), etc. The memory is used to store programs, and the processor executes the programs after receiving the execution instruction.

[0092] The processor may be an integrated circuit chip having the ability to process signals. The above-mentioned processor may be a general-purpose processor, including a central processing unit (CPU), a network processor (NP), etc. The methods, steps, and logic block diagrams disclosed in the embodiments of the present application may be implemented or executed. The general-purpose processor may be a microprocessor or the processor may be any conventional processor, etc.

[0093] The embodiments of the present invention are described with reference to flowcharts of methods, terminal devices (systems), and computer program products according to the embodiments of the present invention. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing terminal device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing terminal device generate means for implementing the functions specified in the flowchart and / or.

[0094] These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing terminal device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce a product including an instruction device that implements the functions specified in the flowchart.

[0095] These computer program instructions can also be loaded onto a computer or other programmable data processing terminal device so that a series of operation steps are executed on the computer or other programmable terminal device to produce computer-implemented processing, whereby the instructions executed on the computer or other programmable terminal device provide steps for implementing the functions specified in the flowchart.

[0096] The technical principle of the present invention is described above in conjunction with specific embodiments. These descriptions are only for explaining the principle of the present invention and cannot be interpreted as limiting the protection scope of the present invention in any way. Based on the explanations herein, those skilled in the art can associate other specific implementations of the present invention without paying creative labor, and these methods will fall within the protection scope of the claims of the present invention.

Claims

1. A method for rapid quantification of local and transport contributions of pollutants, characterized in that: The following steps are involved: Meteoinfo software is used to calculate the concentration threshold of a certain pollutant to obtain gridded data results; According to the grid data results, ArcGIS software is used to calculate the grid point assignment to obtain the absolute proportion of each grid point to the assessment site; Correct the concentration according to the absolute proportion of each grid point to the assessment site and output the result; The calculation using Meteoinfo software includes the following steps: Determine the concentration threshold of a pollutant to be assessed. After determining the assessment point, select the concentration threshold of a pollutant to be assessed; Use Meteoinfo software to obtain gridded data results; The ArcGIS software includes the following steps when performing grid point value calculation: (1) Add a city or county base map; (2) Add the gridded data results to obtain the PSCF.shp file; (3) Change the coordinate system of the PSCF.shp file to a projected coordinate system; Set parameters for the PSCF.shp file and calculate the main grid area; (4) Identify the elements of the PSCF.shp file; (5) Calculate the sub-grid area of ​​the PSCF.shp file after element identification; (6) Calculate the proportion of the secondary grid area to the primary grid area; (7) Calculate the absolute proportion of each grid Mij; (8) Output grid calculation results; Among them, step (1) contains city information, and step (2) contains grid data results. These two layers are combined into one for subsequent calculations. Step (4) assigns the content of the city or county base map added in step (1) to the PSCF.shp layer of step (2), thereby forming a new layer. The boundary line of the city or county base map in the new layer coincides with the grid of Mij, and the grid of Mij is re-divided. The grid of Mij is divided into smaller grids by the boundary line to obtain a number of sub-grids; The adding of gridded data results includes: adding a shp file of the potential source area obtained by Meteoinfo software to obtain a PSCF layer; The setting of parameters for the PSCF.shp file and the main grid calculation includes: right-clicking the PSCF layer, selecting "Open Attribute Table" - "Add Field" - naming it as "MJ" - selecting double precision in "Type", right-clicking "Calculate Geometry", and changing the unit to "Square Kilometers"; The element identification of the PSCF.shp file includes: "Analysis Tools" - "Identification", input relevant elements, click "OK", a new layer appears, named "PSCFBS"; The sub-grid area calculation includes: right-clicking the PSCFBS layer, selecting "Open Attribute Table" - "Add Field" - naming it "BSMJ" - selecting double precision in "Type", and clicking "OK"; right-clicking the "BSMJ" field - "Calculate Geometry", changing the unit to "Square Kilometers", and clicking "OK"; The calculation of the proportion of the secondary grid area in the primary grid area includes: right-clicking the PSCFBS layer, selecting "Open Attribute Table" - "Add Field" - naming it "ZB" - selecting double precision in "Type", and clicking "OK"; right-clicking the "ZB" field - "Field Calculator", clicking "MJBS / MJ", and clicking "OK"; The calculation of the absolute proportion of Mij in each grid includes: right-clicking the PSCFBS layer, selecting "Open Attribute Table" - "Add Field" - naming it "Finally" - selecting double precision in "Type", and clicking "OK"; right-clicking the "Finally" field - "Field Calculator", clicking "Mij*ZB", and clicking "OK"; Outputting the grid calculation results includes: selecting "Conversion Tools" - "Table to Excel", inputting relevant elements and output location and name, and outputting Excel of the grid calculation results; The corrected concentrations include: The concentration threshold of a pollutant under evaluation is used as the local concentration; The mean concentration of pollutants exceeding the threshold is used as the base of transmission concentration for other areas; Multiply the transmission concentration base by the absolute proportion obtained to obtain the contribution concentration of other regions; The contribution ratio of different regions and local areas is calculated by the sum and ratio of local concentration and other regions’ contribution concentration: The contribution ratio of a certain region = transmission concentration in this region / (local concentration + transmission concentration in this region).

2. A method for rapid quantification of local and transport contributions of pollutants according to claim 1, characterized in that: The method for determining the concentration threshold of a pollutant to be evaluated includes the following methods: using the random forest method, PM 2.5 The 10th percentile of pollutant concentration was selected as the threshold, CO, PM 10 , SO2, NO2 and O3 take the average of excellent and low-value good as the threshold, and the low-value good AQI is <75.

3. A device for rapid quantification of local and transport contributions of pollutants, comprising a memory, a processor, and a computer program stored in the memory, characterized in that: The processor executes the computer program to implement the steps of a method for rapid quantification of local and transport contributions of pollutants as claimed in claim 1 or 2.

4. A computer-readable storage medium having a computer program / instruction stored thereon, characterized in that: When the computer program / instructions are executed by a processor, the steps of a method for rapid quantification of local and transport contributions of pollutants as claimed in claim 1 or 2 are implemented.

5. A computer program product comprising a computer program / instructions, characterized in that When the computer program / instructions are executed by a processor, the steps of a method for rapid quantification of local and transport contributions of pollutants as claimed in claim 1 or 2 are implemented.

Images