Grid-based emission inventory compilation method based on pollution zones

Abstract

The application provides a grid-based emission inventory compiling method based on pollution zones, relates to the geographic information system technology and the field of spatial analysis, and comprises the following steps: determining the geographic information grid of a research area; generating a pollution zone vector diagram of an emission source based on the emission amount of a target pollutant of each emission source in the research area and the migration and transformation characteristic information of the target pollutant, and determining the pollutant emission amount per unit area of the pollution zone; determining the projection area of the pollution zone on the grid by superimposing the pollution zone vector diagram and the geographic information grid, thereby determining the target pollutant emission amount of each pollution zone on each grid; and establishing the emission inventory of the pollutant by aggregating and statistically analyzing the target pollutant emission amount from different pollution zones in each grid. The application introduces the pollution zone into the grid-based inventory compiling method, so that the grid-based inventory is more in line with the actual situation and is convenient for pollution control research.

Description

Technical Field

[0001] This application relates to geographic information system (GIS) technology and spatial analysis, and in particular to a method for compiling a gridded emission inventory based on polluted areas. Background Technology

[0002] With the rapid development of science and technology, human activities have gradually become a significant factor affecting the ecological environment. Emissions from industrial production, transportation and shipping, agriculture, and residential life all contribute to environmental pressure. Due to the spatial heterogeneity of emission sources, the environmental impacts of emitted pollutants also vary. High spatial resolution pollutant emission inventories are essential for in-depth and accurate target pollutant control, treatment, and policy formulation.

[0003] Most existing pollutant emission inventory compilation methods are "top-down," which primarily rely on macroeconomic or demographic statistics to estimate the total stock of pollutant emissions. For example, they might acquire gridded population data, calculate the percentage of the population in each grid relative to the total population of the region, and then allocate the total pollutant emissions to each grid according to this percentage, ultimately obtaining the region's pollutant emission inventory. Some studies also allocate pollutant emissions based on road network density, land use type, or building area within the grid. All of these methods require the assumption that pollutant emissions are correlated with a certain indicator before allocation, and a proxy value must be selected based on the correlation coefficient. This increases the uncertainty of the results and introduces errors. Furthermore, commonly used population density data often comes from land use data and nighttime light indices, which suffer from data lag, low accuracy, and poor timeliness. Summary of the Invention

[0004] This application provides a grid-based emission inventory compilation method based on polluted areas, which aims to incorporate the pollution radius of pollutant emission sources into the inventory compilation method to improve the accuracy and reliability of the emission inventory.

[0005] The method for compiling a grid-based emission inventory based on contaminated areas provided in this application includes the following steps:

[0006] The study area is gridded to obtain a geographic information grid of the study area.

[0007] Obtain information on the emission amounts of target pollutants from each emission source within the study area, as well as the migration and transformation characteristics of the target pollutants;

[0008] Based on the emission amount of the target pollutant and the migration and transformation characteristics information, a vector map of the pollution area of ​​the emission source is generated, and the area of ​​each pollution area is determined in the vector map of the pollution area of ​​the emission source.

[0009] Based on the emission amount of the target pollutant and the area of ​​the polluted area, determine the pollutant emission amount per unit area of ​​each polluted area;

[0010] The vector map of the polluted area and the geographic information grid are overlaid to determine the projected area of ​​each polluted area on the corresponding grid.

[0011] The target pollutant emission amount for each polluted area on the corresponding grid is determined based on the pollutant emission amount per unit area of ​​each polluted area and the projected area of ​​each polluted area on the corresponding grid.

[0012] The emission amounts of target pollutants from different polluted areas within each grid are summarized and statistically analyzed to establish an emission inventory of pollutants.

[0013] The gridded emission inventory compilation method based on polluted areas provided in this application involves: first, gridding the study area to obtain a geographic information grid; second, acquiring the emission amounts and migration / transformation characteristics of target pollutants from each emission source within the study area; third, generating a vector map of the polluted areas from the emission sources based on the emission amounts and migration / transformation characteristics, and determining the area of ​​each polluted area within the vector map; fourth, determining the pollutant emission amount per unit area of ​​each polluted area based on the emission amounts of the target pollutants and the area of ​​the polluted area; fifth, overlaying the polluted area vector map and the geographic information grid to determine the projected area of ​​each polluted area on the corresponding grid; sixth, determining the target pollutant emission amount of each polluted area on the corresponding grid based on the pollutant emission amount per unit area of ​​each polluted area and the projected area of ​​each polluted area on the corresponding grid; and finally, summarizing and statistically analyzing the target pollutant emission amounts of each grid to establish an emission inventory of pollutants. In the process of compiling the emission inventory, the emission amount and migration and transformation characteristics of the target pollutants are used as the only input data to determine the vector map of the pollution area of ​​the emission source. There is no need to select any proxy values ​​such as gridded population data or road network density, so as to improve the accuracy and reliability of the high-resolution emission inventory and reduce data errors. The high-resolution emission inventory can provide an important reference for formulating emission reduction strategies and precise prevention and control of regional pollution. Attached Figure Description

[0014] To more clearly illustrate the technical solutions of the embodiments of this application, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are some embodiments of this application.

[0015] Figure 1 A flowchart illustrating a grid-based emission inventory compilation method based on polluted areas, provided as an embodiment of this application;

[0016] Figure 2 A schematic diagram of a geographic information grid provided in an embodiment of this application;

[0017] Figure 3 A schematic diagram of a vector map of a contaminated area provided in an embodiment of this application;

[0018] Figure 4 This is a schematic diagram of an emission inventory provided in one embodiment of this application. Detailed Implementation

[0019] The technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, not all embodiments. Based on the embodiments of this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.

[0020] The flowchart shown in the attached diagram is for illustrative purposes only and does not necessarily include all content and operations / steps, nor does it necessarily have to be performed in the described order. For example, some operations / steps can be broken down, combined, or partially merged, so the actual execution order may change depending on the actual situation.

[0021] This application provides a method for compiling a grid-based emission inventory based on polluted areas. This method can be applied to terminal devices, such as mobile phones, tablets, laptops, desktop computers, personal digital assistants, and wearable devices.

[0022] The following detailed description of some embodiments of this application is provided in conjunction with the accompanying drawings. Unless otherwise specified, the following embodiments and features can be combined with each other.

[0023] Please refer to Figure 1 , Figure 1 This is a flowchart illustrating a gridded emission inventory compilation method based on polluted areas, provided as an embodiment of this application.

[0024] like Figure 1 As shown, the gridded emission inventory compilation method based on polluted areas includes steps S101 to S107.

[0025] Step S101: Grid the study area to obtain the geographic information grid of the study area.

[0026] For example, the study area corresponding to the emission inventory to be compiled is determined. It is understood that the study area includes, but is not limited to, countries, cities, regions, etc.

[0027] For example, the study area is gridded to create a geographic information grid covering the study area. It should be understood that the geographic information grid is a spatial grid that contains multiple grids and completely covers the study area.

[0028] In some embodiments, the study area is gridded to obtain a geographic information grid of the study area, including: importing the extent data of the study area into a geographic information system (GIS) to obtain a geographic information grid of the study area.

[0029] For example, vector data of the study area is acquired; the vector data is then meshed based on a preset mesh precision to create a spatial mesh covering the study area.

[0030] Understandably, this involves acquiring vector data corresponding to the study area and inputting it into a Geographic Information System (GIS). The GIS's segmentation tools then be used to grid the study area. Specifically, the vector data corresponding to the study area stores the geometric location and attribute information of geographic features; the acquired vector data includes the area, boundary range, and coordinates of the study area.

[0031] In the specific implementation process, an appropriate grid precision can be selected and input into the Geographic Information System (GIS) according to the actual user needs. This allows for gridding processing based on the user-defined grid precision, improving the applicability of the method. Grid precision includes, but is not limited to, 1°×1° latitude and longitude grid precision, or 1km×1km kilometer grid precision.

[0032] By gridding the study area, the distribution of target pollutant emission sources in each grid can be determined. Geographic information grids can provide a reference for formulating localized emission reduction strategies and precise pollution prevention and control.

[0033] Step S102: Obtain the emission amount of the target pollutant from each emission source within the study area and the migration and transformation characteristics of the target pollutant.

[0034] For example, after determining the study area, at least one emission source of pollutants located within the study area is obtained. It is understood that the emission source includes, but is not limited to, sewage discharge sources, harmful gas emission sources, etc.

[0035] For example, the target pollutant emitted by the pollutant emission source and the emission amount of the target pollutant can be obtained through data crawling or by obtaining it on the pollution source monitoring platform.

[0036] For example, information on the migration and transformation characteristics of the target pollutant is obtained, wherein the migration characteristics of the target pollutant are obtained from literature or national standards. It should be noted that after determining the emission volume of the target pollutant, the pollution radius of the target pollutant can be determined by combining the obtained migration and transformation characteristics information. For instance, when the target pollutant is industrial wastewater, the pollution radius of the industrial wastewater is determined based on the location of the wastewater source, the emission volume, and the migration and transformation characteristics of the industrial wastewater in the current aquatic environment, thereby identifying the pollution area of ​​the industrial wastewater.

[0037] In a specific implementation, the coordinates of the emission source are determined by the latitude and longitude of the emission source's location.

[0038] In some embodiments, the migration and transformation characteristics of the target pollutant are determined based on the emission time and emission space of the target pollutant.

[0039] It is understandable that the migration and transformation characteristics of the same target pollutant will differ at different emission times and in different emission spaces. Emission time indicates the duration of emission of the target pollutant, while emission space indicates the medium to which the target pollutant is emitted, such as air or water.

[0040] It should be understood that, in this application, migration and transformation characteristics information can be used to indicate changes that occur during the migration and transformation process of a target pollutant, such as the pollution caused to a river when wastewater is discharged into it, in order to determine the pollution area of ​​the target pollutant.

[0041] Step S103: Based on the emission amount of the target pollutant and the migration and transformation characteristics information, generate a vector map of the emission source pollution area, and determine the area of ​​each pollution area in the vector map of the emission source pollution area.

[0042] For example, a vector map of the pollution area of ​​the emission source is generated based on the emission amount of the target pollutant and the migration and transformation characteristics of the target pollutant. It should be noted that the vector map of the pollution area of ​​the emission source includes the pollution area of ​​the target pollutant emitted by multiple emission sources. Therefore, the area of ​​each pollution area included in the vector map of the emission source pollution area can be determined.

[0043] For example, target pollutants include industrial wastewater and exhaust gases, traffic air pollutants, and takeout packaging waste, which, through atmospheric transport, surface runoff, and transportation, create an environmental pollution zone around them. In other words, the pollution zone is the area actually polluted by the target pollutants emitted from the pollutant emission sources.

[0044] In some embodiments, generating a polluted area vector map layer based on the emission amount of the target pollutant from the pollutant emission source and the migration and transformation characteristics of the target pollutant, and determining the area of ​​each polluted area in the polluted area vector map layer, includes: inputting the polluted area vector map layer into a geographic information system (GIS) to obtain the area of ​​each polluted area.

[0045] For example, each pollution zone is identified in the vector map of the pollution zone of the emission source using a geographic information system (GIS), and the area corresponding to each pollution zone is determined.

[0046] In the specific implementation process, the pollutant emission source is a point source, and the polluted area can be defined as a circular area. Understandably, the area of ​​this circular area is the area corresponding to the polluted area.

[0047] Please see Figure 2 , Figure 3 , Figure 2 This is a schematic diagram of a gridded layer provided in an embodiment of this application. Figure 3 This is a schematic diagram of the contaminated area provided in an embodiment of this application.

[0048] Step S104: Based on the emission amount of the target pollutant and the area of ​​the polluted area, determine the pollutant emission amount per unit area of ​​each polluted area.

[0049] For example, after determining the area of ​​the contaminated zone, the pollutant emission per unit area of ​​the contaminated zone is determined based on the emission of the target pollutant within the contaminated zone.

[0050] In the specific implementation process, determining the pollutant emission per unit area of ​​each polluted area based on the emission amount of the target pollutant and the area of ​​the polluted area includes: calculating the pollutant emission per unit area of ​​each polluted area according to the following formula:

[0051]

[0052] For example, It is a positive integer used to indicate the first... One contaminated area; The emission amount of the target pollutant; The area of ​​the contaminated zone; This represents the amount of pollutants emitted per unit area of ​​the polluted zone.

[0053] Determine the pollutant emission per unit area of ​​the polluted area, and then determine the target pollutant emission of the polluted area on the grid by using the pollutant emission per unit area of ​​the polluted area and the projected area of ​​the polluted area on the grid.

[0054] Step S105: Overlay the polluted area vector map and the geographic information grid to determine the projected area of ​​each polluted area on the corresponding grid.

[0055] For example, a vector map of the contaminated area and a geographic information grid are overlaid, wherein the geographic information grid can be as follows: Figure 2 As shown, the vector diagram of the contaminated area can be seen as follows: Figure 3 As shown, after overlaying the polluted area vector map and the geographic information grid, the polluted areas included in the polluted area vector map are projected onto the geographic information grid, thereby enabling the determination of the grid associated with each polluted area and the projected area of ​​the polluted area on the associated grid.

[0056] It should be noted that a geographic information grid contains multiple grids. In some cases, a single contaminated area may be projected onto multiple grids, and the projected area of ​​a single contaminated area may differ from that of each grid. By determining the grid corresponding to the contaminated area, i.e., the grid where the contaminated area's projection resides, the projected area of ​​a single contaminated area and its respective grid can be calculated.

[0057] In some embodiments, the step of overlaying the polluted area vector map and the geographic information grid to determine the projected area of ​​each polluted area on the corresponding grid includes: inputting the polluted area vector map and the geographic information grid into a geographic information system so that the polluted area vector map layer and the geographic information grid are overlaid in the same coordinate system; identifying each polluted area based on a GIS identification tool, determining the grid where each polluted area is located based on the identification, and calculating the projected area of ​​each polluted area on its respective grid.

[0058] In the specific implementation process, the geographic information grid is input into the geographic information system so that the geographic information grid is superimposed on the already input polluted area vector map in the same coordinate system. Based on the GIS identification tool, each polluted area is identified, so that the grid where each polluted area is located can be determined according to the identification, so as to realize the calculation of the projected area of ​​each polluted area on its grid.

[0059] Understandably, since the geographic information grid and the polluted area vector map are layers of two geographic objects in the same study area, and are generated under the same spatial reference system, the geographic information grid and the polluted area vector map can be overlaid in the same coordinate system.

[0060] Step S106: Determine the target pollutant emission amount of each polluted area on the corresponding grid based on the pollutant emission amount per unit area of ​​each polluted area and the projected area of ​​each polluted area on the corresponding grid.

[0061] For example, the projected area of ​​each contaminated area on the corresponding grid is determined, thereby determining the target pollutant emission amount of the contaminated area on the corresponding grid. This allows for the statistical calculation of the target pollutant emission amounts corresponding to different contaminated areas within each grid, and based on the target pollutant emission amounts corresponding to different contaminated areas on the same grid, the total pollutant emission amount for that grid can be determined. It is understood that at least one of the following situations may exist during the calculation process: one contaminated area is projected onto one grid; one contaminated area is projected onto multiple grids; multiple contaminated areas are projected onto one grid; or there is no projection of a contaminated area into the grid.

[0062] For example, when there is no projection of a polluted area on the grid, the target pollutant emission of the grid is determined to be zero.

[0063] For example, the emission amount of the target pollutant corresponding to the polluted area on the corresponding grid is determined based on the emission amount of the target pollutant corresponding to the polluted area and the projected area of ​​the polluted area on the grid.

[0064] In some embodiments, each point in the contaminated area is contaminated to the same degree by the target pollutant.

[0065] It should be understood that in this application, the pollution level of the target pollutant is the same within the contaminated area. Since the pollution level of each point within the contaminated area is the same, the target pollutant of the grid can be determined based on the projected area of ​​the contaminated area on the corresponding grid and the pollutant emission per unit area of ​​the contaminated area.

[0066] In the specific implementation process, the study area contains multiple emission sources. When a single polluted area is projected onto multiple grids, the projected area of ​​a single polluted area and its corresponding grid is determined. For example, the polluted area of ​​emission source a corresponds to grids A and B. The target pollutant emission amount of the polluted area in grids A and B is determined based on the pollutant emission per unit area of ​​the polluted area and the projected area of ​​the polluted area in grids A and B. For example, if the projected area of ​​the polluted area of ​​emission source a in grid A is 15π square meters and the projected area of ​​the polluted area of ​​emission source a in target grid B is 20π square meters, then the pollutant emission per unit area of ​​the polluted area of ​​emission source a is multiplied by 15π square meters and 20π square meters respectively to obtain the target pollutant emission amount of emission source a in grid A and the target pollutant emission amount of emission source a in grid B.

[0067] In the specific implementation process, determining the target pollutant emission amount of each polluted area on the corresponding grid includes: calculating the target pollutant emission amount of each polluted area on the grid corresponding to each polluted area according to the following formula:

[0068]

[0069] in, It is a positive integer used to indicate the first... One contaminated area; For the first The target pollutant emissions for each polluted area within its respective grid; For the first The projected area of ​​each contaminated area on its respective grid; For the first Pollutant emissions per unit area of ​​a polluted zone.

[0070] Step S107: Summarize and statistically analyze the emission amounts of target pollutants from different polluted areas within each grid, and establish an emission inventory of pollutants.

[0071] For example, the emissions of target pollutants from different polluted areas within each grid are aggregated to establish an emissions inventory of pollutants.

[0072] In some embodiments, the step of summarizing and statistically analyzing the target pollutant emissions from different polluted areas within each grid to establish a pollutant emission inventory includes: summarizing and statistically analyzing the target pollutant emissions of the grid according to the following formula:

[0073]

[0074] in, It is a positive integer used to indicate the first... One contaminated area; For the first The target pollutant emissions for each polluted area on the grid.

[0075] In the specific implementation process, if a grid contains projections of pollution areas from emission source a and emission source b, the total pollutant emission of the grid is determined based on the target pollutant emission amounts corresponding to the pollution areas of emission source a and emission source b on that grid. For example, if the projected area of ​​the pollution area of ​​emission source a on the grid is 0.2 m², the projected area of ​​the pollution area of ​​emission source b on the grid is 0.3 m², the pollutant emission per unit area of ​​the pollution area of ​​emission source a is 15 kg / m², and the pollutant emission per unit area of ​​the pollution area of ​​emission source b is 20 kg / m², then the total pollutant emission of the grid is 9 kg.

[0076] It should be noted that the pollutant emission per unit area of ​​the polluted area, as well as the number of grids and projected area corresponding to the polluted area, are illustrative examples. They will be determined based on the actual data obtained and the processing procedures provided above during application. This application does not limit the above data.

[0077] Please see Figure 4 , Figure 4 This is a schematic diagram of a pollutant emission inventory provided in one embodiment of this application.

[0078] The gridded emission inventory compilation method based on polluted areas provided in the above embodiments obtains a geographic information grid of the study area by performing gridding processing on the study area; acquires the emission amount of target pollutants from each emission source within the study area and the migration and transformation characteristics of the target pollutants; generates a pollution area vector map of emission sources based on the emission amount of the target pollutants and the migration and transformation characteristics, and determines the area of ​​each pollution area in the pollution area vector map; determines the pollutant emission amount per unit area of ​​each pollution area based on the emission amount of the target pollutants and the area of ​​the pollution area; and performs a mapping between the pollution area vector map and the geographic information grid. The method involves overlay processing to determine the projected area of ​​each polluted area on its corresponding grid. Based on the pollutant emission per unit area of ​​each polluted area and its projected area on the corresponding grid, the target pollutant emission amount for each polluted area on its corresponding grid is determined. The target pollutant emissions from different polluted areas within each grid are then aggregated to establish a pollutant emission inventory. This method eliminates the need for surrogate value calculations, reducing the time required for such calculations and minimizing errors. This improves the resolution, accuracy, and reliability of the emission inventory. A high-resolution emission inventory provides crucial reference for developing emission reduction strategies and implementing refined regional pollution control. This application innovatively considers the differences and dynamic changes in the pollution range of each pollutant and incorporates the pollution radius into the grid inventory compilation method. A bottom-up approach is used to achieve high-resolution emission inventory compilation, aiming to analyze the pollutant emissions of a single inventory component and then extrapolate the results to the inventory level. Overall, the bottom-up approach can assess the impact of environmental changes on pollutant emissions in more detail than the top-down approach, thus providing more specific guidance for developing source reduction policies.

[0079] It should be understood that the terminology used in this specification is for the purpose of describing particular embodiments only and is not intended to limit the scope of the application. As used in this specification and the appended claims, the singular forms “a,” “an,” and “the” are intended to include the plural forms unless the context clearly indicates otherwise.

[0080] It should also be understood that the term "and / or" as used in this specification and the appended claims refers to any combination and all possible combinations of one or more of the associated listed items, and includes such combinations. It should be noted that, herein, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or system that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or system. Without further limitation, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or system that includes that element.

[0081] The sequence numbers of the embodiments in this application are for descriptive purposes only and do not represent the superiority or inferiority of the embodiments. The above descriptions are merely specific implementations of this application, but the scope of protection of this application is not limited thereto. Any person skilled in the art can easily conceive of various equivalent modifications or substitutions within the technical scope disclosed in this application, and these modifications or substitutions should all be covered within the scope of protection of this application. Therefore, the scope of protection of this application should be determined by the scope of the claims.

Claims

1. A method for grid-based emissions inventory development based on pollution zones, characterized by, The method comprises the following steps: griding the research area to obtain a geographic information grid of the research area; obtaining the emission amount of a target pollutant of each emission source in the research area and migration and transformation characteristic information of the target pollutant, the migration and transformation characteristic information of the target pollutant being determined according to the emission time and emission space of the target pollutant; generating an emission source pollution area vector diagram according to the emission amount of the target pollutant and the migration and transformation characteristic information, and determining the area of each pollution area in the emission source pollution area vector diagram, comprising: determining the pollution radius of the target pollutant according to the emission amount of the target pollutant and the migration and transformation characteristic information, determining the pollution area of the target pollutant according to the pollution radius, thereby obtaining the emission source pollution area vector diagram, and identifying each pollution area in the emission source pollution area vector diagram based on a geographic information system, wherein the pollution area is a circular area, and the area corresponding to the circular area is the area of the pollution area; determining the emission amount of the target pollutant per unit area of each pollution area based on the emission amount of the target pollutant and the area of the pollution area; superimposing the pollution area vector diagram and the geographic information grid to determine the projection area of each pollution area on the corresponding grid, comprising: inputting the pollution area vector diagram and the geographic information grid into the geographic information system to superimpose the pollution area vector diagram layer and the geographic information grid in the same coordinate system; identifying each pollution area based on an identification tool, and determining the grid where each pollution area is located and calculating the projection area of each pollution area on the grid based on the identification; determining the emission amount of the target pollutant of each pollution area on the corresponding grid according to the emission amount of the target pollutant per unit area of each pollution area and the projection area of each pollution area on the corresponding grid; summarizing and statistically the emission amount of the target pollutant from different pollution areas in each grid to establish an emission inventory of the pollutant.

2. The gridded emissions inventory compilation method based on pollution zones of claim 1, wherein, The degree of pollution of each point in the pollution area by the target pollutant is the same.

3. The method of grid-based emissions inventory development based on pollution zones of claim 1, wherein, The griding of the research area to obtain a geographic information grid of the research area comprises: importing shp vector data of the research area into a geographic information system to obtain a geographic information grid of the research area.

4. The gridded emissions inventory compilation method based on pollution zones according to any one of claims 1 to 3, characterized in that, The determination of the emission amount of the target pollutant per unit area of each pollution area based on the emission amount of the target pollutant and the area of the pollution area comprises: the emission amount of the target pollutant per unit area of each pollution area is calculated according to the following formula: wherein, is a positive integer, used to indicate the th pollution area; is the emission amount of the target pollutant; is the area of the pollution area; is the pollutant emission amount per unit area of the pollution area.

5. The method of grid-based emissions inventory development based on pollution zones of claim 1, wherein, The determination of the emission amount of the target pollutant of each pollution area on the corresponding grid comprises: the emission amount of the target pollutant of each pollution area on the corresponding grid of each pollution area is calculated according to the following formula: wherein, is a positive integer, indicating the th pollution zone; is the target pollutant emission of the th pollution zone on the grid where it is located; is the projected area of the th pollution zone on the grid where it is located; is the pollutant emission per unit area of the th pollution zone.

6. The method of grid-based emissions inventory development based on pollution zones of claim 1, wherein, The summarizing and statistical of the emission amount of the target pollutant from different pollution areas in each grid to establish an emission inventory of the pollutant comprises: the emission amount of the target pollutant of the grid is summarized and statistically according to the following formula: wherein, is a positive integer indicating the th pollution zone; is the th pollution zone target pollutant emission on the grid.

Images