Greenhouse gas management and control partition division method and system based on spatial difference distribution

Abstract

The embodiment of the invention discloses a greenhouse gas management and control partition division method and system based on spatial difference distribution, and belongs to the technical field of greenhouse gas emission, and the method comprises the steps: evaluating the accessibility of parameters according to a target region, unifying the spatial precision of evaluation work, combining the accessibility of the evaluation parameters, and selecting an evaluation calculation method; according to the selected evaluation calculation method, collecting data required by evaluation, calculating a land ecosystem carbon sink function importance index and an atmospheric pollutant and greenhouse gas emission influence importance index, and evaluating the land ecosystem carbon sink function importance index and the atmospheric pollutant and greenhouse gas emission influence importance index. According to the embodiment of the invention, a complex traceability operation process can be omitted, the guidance and operability on practice are improved, the delimited threshold value of the greenhouse gas priority protection area is determined, and the problem that the threshold value is difficult to determine in a traditional grading evaluation method is effectively solved.

Description

Technical Field

[0001] The present application belongs to the technical field of greenhouse gas emissions, and specifically, relates to a method and system for zoning greenhouse gas control based on spatial differential distribution. Background Art

[0002] We currently attach great importance to the issue of climate change, and actively responding to climate change will be regarded as a major strategy and policy orientation for economic and social development. Monitoring the concentration and changing trends of greenhouse gases in the atmosphere is an important basic work.

[0003] Among them, spatial zoning is an important means to achieve refined and differentiated control of pollutant and greenhouse gas emissions. Existing greenhouse gas spatialization technologies are mainly based on the preparation and generation of gridded emission inventories, which allocate the total regional greenhouse gas emissions to grid units of specific spatial resolution according to certain rules to achieve spatialization of total emissions; some use remote sensing monitoring methods to identify areas with high greenhouse gas concentrations based on differences in spatial distribution of concentrations.

[0004] However, existing technologies have defects. The existing greenhouse gas control zoning demarcation methods lack systematicity and only consider carbon emissions or greenhouse gas concentrations, but lack consideration of carbon absorption. At the same time, it is difficult to locate the real emission source based on the identification of high-concentration areas based on spatial differences in concentration distribution, resulting in poor guidance for practical applications. There is a lack of systematic quantitative assessment methods for spatial area division and threshold determination, and accurate evaluation results cannot be obtained. Summary of the invention

[0005] In order to solve the above problems and technical defects, the embodiment of the present application adopts the following technical solution, which is a method for dividing greenhouse gas control zones based on spatial differential distribution, including the following steps:

[0006] Step 1: Unify the spatial accuracy of the assessment work based on the availability of assessment parameters in the target area;

[0007] Step 2: Based on the importance of the carbon sink function of the terrestrial ecosystem in the target area and the emission status of atmospheric pollutants and greenhouse gases, combined with the availability of assessment parameters, select the assessment calculation method;

[0008] Step 3: Collect the data required for the evaluation according to the selected evaluation calculation method;

[0009] Step 4: Based on the evaluation calculation method and the data required for the evaluation, calculate the importance index of the carbon sink function of the terrestrial ecosystem and the importance index of the impact of atmospheric pollutants and greenhouse gas emissions;

[0010] Step 5: Based on the cumulative contribution calculation principle and the calculation results, evaluate the importance index of terrestrial ecosystem carbon sink function and the importance index of atmospheric pollutants and greenhouse gas emissions;

[0011] Step 6: Based on the assessment results, zoning is carried out to divide the corresponding areas into greenhouse gas priority protection areas, greenhouse gas key control areas and greenhouse gas general control areas.

[0012] Preferably, in step 1, the evaluation parameters of the target area are screened and identified, the availability of the evaluation parameters after screening and identification is obtained, and the spatial accuracy of the target area for carrying out the evaluation work is determined and unified according to the availability of the evaluation parameters, and the spatial accuracy is used as the basic evaluation unit;

[0013] The basic assessment units include: the spatial resolution of evaluating the importance of carbon sink function of terrestrial ecosystems, and the spatial resolution of evaluating the importance of the impact of atmospheric pollutants and greenhouse gas emissions.

[0014] Preferably, after collecting the data required for evaluation in step 3, the data required for evaluation is preprocessed to convert the data required for evaluation into gridded raster data;

[0015] The data required for the assessment include basic geographic information data, land use data, vegetation cover data, remote sensing images, industrial enterprise location data, road network data, population data, and fused emission inventory data of atmospheric pollutants and greenhouse gases.

[0016] Preferably, the terrestrial ecosystem carbon sink function importance index in step 4 is calculated by characterizing the carbon sink capacity of regional forest ecosystems, and the forest ecosystem carbon sink accounting method is as follows:

[0017] Obtain vegetation coverage and its change types, area, vegetation and soil activity data within the accounting boundary of the target area through investigation and statistics;

[0018] Combined with the emission factor, the carbon stock change within the boundary is calculated and converted into carbon dioxide equivalent. The carbon dioxide equivalent calculation formula is as follows:

[0019]

[0020] Among them, Sink LAND,t is the carbon sink of terrestrial ecosystems within the accounting boundary in year t, is the carbon storage of carbon pool p in vegetation cover and its change type i within the boundary in the initial year t1, is the carbon storage of carbon pool p in vegetation cover and its change type i within the boundary of the accounting period in year t2, is the activity data of carbon pool p in vegetation cover and its change type i within the accounting boundary in year t, is the emission factor of carbon pool p in vegetation cover and its change type i within the accounting boundary in year t, i is vegetation cover and its change type, p is carbon pool type, T is the accounting time interval, t is the accounting year, is a constant for converting carbon into carbon dioxide equivalent.

[0021] Furthermore, the importance index of air pollutant and greenhouse gas emissions in step 4 is calculated by characterizing the amount of air pollutant emissions and greenhouse gas emissions in the spatial grid unit. First, the total amount of emissions is calculated, the emission space is allocated, and the emission equivalent is converted. Then, the air pollutants are calculated to obtain the pollutant emission equivalent. Finally, the pollutant emission equivalent is counted to the corresponding grid unit. The calculation formula is as follows:

[0022]

[0023] Among them, E AP is the main pollutant equivalent, Q i is the emission of the i-th pollutant, W i is the equivalent value of the ith pollutant.

[0024] Furthermore, the emission spatial allocation includes point source processing, mobile source processing and surface source processing. Point source processing is directly positioned on the grid unit according to the longitude and latitude coordinates. Mobile source processing uses a substitute parameter weight method to allocate road network data emissions to the grid. Surface source processing uses a substitute parameter weight method to allocate population data emissions to the grid.

[0025] Furthermore, in step 5, the importance index has three levels, including general importance, important and extremely important.

[0026] Furthermore, the process of evaluating the importance index in step 5 is as follows:

[0027] The calculation results are calculated using a grid calculator to obtain normalized importance evaluation grid maps, including importance evaluation grid maps of terrestrial ecosystem carbon sink function, and importance evaluation grid maps of atmospheric pollutants and greenhouse gas emissions;

[0028] Derive the raster data attribute table from the importance evaluation raster map. The attribute table records the importance value of each raster pixel. The importance values ​​are arranged in descending order and the accumulated importance value is calculated.

[0029] Two ratio values ​​are preset to obtain the grid values ​​corresponding to the two ratio values ​​of the cumulative importance value in the total value, and the two grid values ​​are used as the dividing points for the evaluation and classification of the importance index of terrestrial ecosystem carbon sink function, atmospheric pollutants and greenhouse gas emissions;

[0030] The reclassification tool was used in combination with the cut-off points to evaluate and grade the importance index of terrestrial ecosystem carbon sink function, atmospheric pollutants and greenhouse gas emissions.

[0031] Preferably, the specific process of partitioning in step 6 is as follows:

[0032] The extremely important areas and important areas in the grid map of the importance of carbon sink function of terrestrial ecosystems are designated as priority protection areas for greenhouse gases;

[0033] Overlay analysis is performed on the extremely important areas and important areas in the greenhouse gas and air pollutant emission importance assessment grid map, and the areas with spatial overlap are regarded as key greenhouse gas control areas;

[0034] Other areas are designated as general greenhouse gas control areas.

[0035] The greenhouse gas control zoning system based on spatial differential distribution includes:

[0036] The unit accuracy identification module unifies the spatial accuracy of the assessment work according to the availability of assessment parameters in the target area;

[0037] The evaluation calculation selection module selects the evaluation calculation method based on the importance of the carbon sink function of the terrestrial ecosystem in the target area and the emission status of atmospheric pollutants and greenhouse gases, combined with the availability of evaluation parameters;

[0038] The data collection and processing module collects the data required for the evaluation according to the selected evaluation calculation method;

[0039] Importance index calculation module, which calculates the importance index of terrestrial ecosystem carbon sink function and the importance index of atmospheric pollutants and greenhouse gas emissions based on the evaluation calculation method and the data required for the evaluation;

[0040] The importance assessment and grading module evaluates the importance index of terrestrial ecosystem carbon sink function, atmospheric pollutants and greenhouse gas emissions impact based on the cumulative contribution calculation principle and combined with the calculation results;

[0041] The zoning demarcation module carries out zoning demarcation based on the assessment results, dividing the corresponding areas into greenhouse gas priority protection areas, greenhouse gas key control areas and greenhouse gas general control areas.

[0042] Compared with the prior art, the beneficial effects of the embodiments of the present application are:

[0043] (1) This application considers carbon absorption, carbon emissions, and the coordination of pollution reduction and carbon reduction to obtain the spatial resolution of the importance evaluation of the carbon sink function of terrestrial ecosystems and the spatial resolution of the importance evaluation of the impact of atmospheric pollutants and greenhouse gas emissions. It takes into account a comprehensive and systematic technical method for the delineation of greenhouse gas control zones covering the entire region, which can effectively improve the efficiency of the coordinated control of carbon emissions and pollution emissions. At the same time, compared with the traditional concentration-based division method, it not only saves the complicated tracing operation process, but also improves the guidance and operability of practice;

[0044] (2) This application uses a method for delineating key greenhouse gas control areas based on the importance index of atmospheric pollutants and greenhouse gas emissions, introduces the carbon sink function importance index as the basis for grading the importance of regional carbon sink functions, introduces an atmospheric pollutant and greenhouse gas emission equivalent algorithm and an emission impact importance index as the basis for grading the importance of regional atmospheric pollutants and greenhouse gas emissions, and determines the threshold for delineating greenhouse gas priority protection areas as the basis for importance evaluation and grading, effectively solving the problem of difficulty in determining thresholds in traditional grading assessment methods. BRIEF DESCRIPTION OF THE DRAWINGS

[0045] In the attached picture:

[0046] Figure 1 A schematic diagram of the method steps of an embodiment of the present application;

[0047] Figure 2 A schematic diagram of surface vegetation coverage according to an embodiment of the present application;

[0048] Figure 3 This is a schematic diagram of the distribution of point sources of industrial enterprise locations according to an embodiment of the present application;

[0049] Figure 4 A schematic diagram of road distribution network data of an embodiment of the present application;

[0050] Figure 5 A schematic diagram of population distribution data of an embodiment of the present application;

[0051] Figure 6 A schematic diagram of the carbon sink function importance index of an embodiment of the present application;

[0052] Figure 7 A schematic diagram of the importance index of air pollutant emissions according to an embodiment of the present application;

[0053] Figure 8 A schematic diagram of the greenhouse gas emission impact importance index according to an embodiment of the present application;

[0054] Fig. 9 A schematic diagram of the importance classification of carbon sink functions in an embodiment of the present application;

[0055] Fig.10 This is a schematic diagram of the importance classification of the impact of air pollutant emissions in an embodiment of the present application;

[0056] Fig.11 A schematic diagram of the importance classification of the impact of greenhouse gas emissions according to an embodiment of the present application;

[0057] Fig.12 A schematic diagram of greenhouse gas control zones according to an embodiment of the present application;

[0058] Fig.13 Schematic diagram of the system structure of an embodiment of the present application. DETAILED DESCRIPTION

[0059] In order to make the purpose, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are part of the embodiments of the present application, not all of the embodiments. Generally, the components of the embodiments of the present application described and shown in the drawings here can be arranged and designed in various different configurations.

[0060] Example 1

[0061] like Figure 1 As shown, the greenhouse gas control zoning method based on spatial differential distribution includes the following steps:

[0062] Unify the spatial accuracy of assessment work based on the availability of assessment parameters in the target area;

[0063] Screen and identify the evaluation parameters of the target area, obtain the availability of the evaluation parameters after screening and identification, determine and unify the spatial accuracy of the target area for evaluation work based on the availability of the evaluation parameters, and use spatial accuracy as the basic evaluation unit;

[0064] The basic assessment units include: the spatial resolution of evaluating the importance of carbon sink function of terrestrial ecosystems, and the spatial resolution of evaluating the importance of the impact of atmospheric pollutants and greenhouse gas emissions.

[0065] Select the assessment calculation method based on the importance of the carbon sink function of the terrestrial ecosystem and the emission status of atmospheric pollutants and greenhouse gases in the target area, combined with the availability of assessment parameters;

[0066] Collect the data required for the evaluation according to the selected evaluation calculation method;

[0067] After collecting the data required for the evaluation, the data required for the evaluation are pre-processed and converted into gridded raster data;

[0068] The data required for the assessment include basic geographic information data, land use data, vegetation cover data, remote sensing images, industrial enterprise location data, road network data, population data, and fused emission inventory data of atmospheric pollutants and greenhouse gases.

[0069] Based on the evaluation calculation method and the data required for the evaluation, calculate the importance index of the carbon sink function of terrestrial ecosystems and the importance index of the impact of atmospheric pollutants and greenhouse gas emissions;

[0070] The importance index of terrestrial ecosystem carbon sink function is calculated by characterizing the carbon sink capacity of regional forest ecosystems. The forest ecosystem carbon sink accounting method is as follows:

[0071] Obtain vegetation coverage and its change types, area, vegetation and soil activity data within the accounting boundary of the target area through investigation and statistics;

[0072] Combined with the emission factor, the carbon stock change within the boundary is calculated and converted into carbon dioxide equivalent. The carbon dioxide equivalent calculation formula is as follows:

[0073]

[0074] Among them, Sink LAND,t is the carbon sink of terrestrial ecosystems within the accounting boundary in year t, is the carbon storage of carbon pool p in vegetation cover and its change type i within the boundary in the initial year t1, is the carbon storage of carbon pool p in vegetation cover and its change type i within the boundary of the accounting period in year t2, is the activity data of carbon pool p in vegetation cover and its change type i within the accounting boundary in year t, is the emission factor of carbon pool p in vegetation cover and its change type i within the accounting boundary in year t, i is vegetation cover and its change type, p is carbon pool type, T is the accounting time interval, t is the accounting year, is a constant for converting carbon into carbon dioxide equivalent.

[0075] The importance index of atmospheric pollutant and greenhouse gas emissions is characterized by the amount of atmospheric pollutant emissions and greenhouse gas emissions in the spatial grid unit. First, the total amount of emissions is calculated, the emission space is allocated, and the emission equivalent is converted. Then, the atmospheric pollutants are calculated to obtain the pollutant emission equivalent. Finally, the pollutant emission equivalent is counted to the corresponding grid unit. The calculation formula is as follows:

[0076]

[0077] Among them, E AP is the main pollutant equivalent, Q i is the emission of the i-th pollutant, W i is the equivalent value of the ith pollutant.

[0078] The spatial allocation of emissions includes point source processing, mobile source processing and area source processing. Point source processing is directly located on the grid unit according to the longitude and latitude coordinates. Mobile source processing uses the substitute parameter weight method to allocate road network data emissions to the grid. Area source processing uses the substitute parameter weight method to allocate population data emissions to the grid.

[0079] Based on the calculation principle of cumulative contribution and combined with the calculation results, the importance index of terrestrial ecosystem carbon sink function and the importance index of atmospheric pollutants and greenhouse gas emissions are evaluated;

[0080] There are three levels of importance index, including general importance, important and extremely important.

[0081] The process of evaluating the importance index is as follows:

[0082] The calculation results are calculated using a grid calculator to obtain normalized importance evaluation grid maps, including importance evaluation grid maps of terrestrial ecosystem carbon sink function, and importance evaluation grid maps of atmospheric pollutants and greenhouse gas emissions;

[0083] Derive the raster data attribute table from the importance evaluation raster map. The attribute table records the importance value of each raster pixel. The importance values ​​are arranged in descending order and the accumulated importance value is calculated.

[0084] Two ratio values ​​are preset to obtain the grid values ​​corresponding to the two ratio values ​​of the cumulative importance value in the total value, and the two grid values ​​are used as the dividing points for the evaluation and classification of the importance index of terrestrial ecosystem carbon sink function, atmospheric pollutants and greenhouse gas emissions;

[0085] The reclassification tool was used in combination with the cut-off points to evaluate and grade the importance index of terrestrial ecosystem carbon sink function, atmospheric pollutants and greenhouse gas emissions.

[0086] Based on the assessment results, zoning is carried out and the corresponding areas are divided into greenhouse gas priority protection areas, greenhouse gas key control areas and greenhouse gas general control areas.

[0087] The specific process of zoning is as follows:

[0088] The extremely important areas and important areas in the grid map of the importance of carbon sink function of terrestrial ecosystems are designated as priority protection areas for greenhouse gases;

[0089] Overlay analysis is performed on the extremely important areas and important areas in the greenhouse gas and air pollutant emission importance assessment grid map, and the areas with spatial overlap are regarded as key greenhouse gas control areas;

[0090] Other areas are designated as general greenhouse gas control areas.

[0091] Example 2

[0092] Hunan Province was selected as an implementation case of this system. The spatial resolution of the basic assessment unit was first determined. According to the availability of regional assessment parameters in Hunan Province, the spatial resolution for evaluating the importance of ecosystem carbon sink function was determined to be 30 meters, and the spatial resolution for evaluating the importance of atmospheric pollutants and greenhouse gas emissions was determined to be 3 km.

[0093] Reconstruct the calculation equations of the importance index of ecosystem carbon sink function, and the importance index of atmospheric pollutants and greenhouse gas emissions;

[0094] The carbon sink function importance index is characterized by the carbon sink capacity of regional forest ecosystems. The forest ecosystem carbon sink accounting method is as follows:

[0095] Through surveys, statistics and other means, we obtain the vegetation coverage and its change types, areas, vegetation, soil and other activity data within the accounting boundary, calculate the carbon storage changes within the boundary in combination with emission factors, and convert them into carbon dioxide equivalents:

[0096]

[0097] Among them, Sink LAND,t is the carbon sink of the terrestrial ecosystem within the accounting boundary in year t. A negative value indicates carbon emissions, in tons of carbon dioxide equivalent / year (t CO2e / yr);

[0098] is the carbon storage of carbon pool p in vegetation cover and its change type i within the boundary in the initial year t1 of the accounting period, in tons of carbon (t C);

[0099] is the carbon storage of carbon pool p in vegetation cover and its change type i within the boundary of year t2 at the end of accounting, in tons of carbon (t C);

[0100] The activity data of carbon pool p in vegetation cover and its change type i within the accounting boundary in year t. Its functional unit is related to the type of activity data, for example, the area unit is hectare (hm 2 ), volume unit is cubic meter (m 3 ), the unit of mass is ton (t);

[0101] is the emission factor of carbon pool p in vegetation cover and its change type i within the accounting boundary in year t, in tons of carbon / activity data unit;

[0102] i is the vegetation cover and its change type, dimensionless;

[0103] p is the carbon pool type, dimensionless;

[0104] T is the time interval for calculation, T = t2-t1, unit: year (yr);

[0105] t is the accounting year, t1≤t≤t2, the unit is: year (yr);

[0106] is a constant for converting carbon into carbon dioxide equivalent.

[0107] The importance index of air pollutant and greenhouse gas emissions is represented by the amount of air pollutant emissions and greenhouse gas emissions in the spatial grid units, respectively. It specifically includes three steps: total emission accounting, spatial allocation, and emission equivalent conversion.

[0108] The total emission accounting method refers to the relevant requirements in the Technical Guidelines for the Preparation of Integrated Emission Inventories of Air Pollutants and Greenhouse Gases (Trial Implementation);

[0109] The spatial distribution of emissions is processed according to point sources, mobile sources and surface sources. Point sources are directly located on grid cells according to their longitude and latitude coordinates. Mobile sources and surface sources are distributed to grids based on road network data and population data using the substitute parameter weight method.

[0110] If multiple air pollutants are involved, they should be converted into pollutant emission equivalents according to the following formula, and then added and counted to the corresponding grid unit;

[0111] The emission equivalents of major pollutants are calculated according to the equivalent values ​​of various pollutants stipulated in the Environmental Protection Tax Law. The calculation formula is as follows:

[0112]

[0113] Among them, E AP is the main pollutant equivalent, dimensionless, expressed as AP;

[0114] Q i is the emission of the ith pollutant, in kilograms (kg);

[0115] W i is the equivalent value of the ith pollutant, in kilograms (kg).

[0116] Then collect the basic data required for the above model calculation, such as basic geographic information data, vegetation cover data, industrial enterprise location, road network data, population data, air pollutant and greenhouse gas fusion emission inventory data, etc.;

[0117] Vegetation cover data such as Figure 2 As shown, Figure 2 The resolution is 30 meters for Hunan, and the locations of industrial enterprises are as follows: Figure 3 As shown, Figure 3 Taking Hunan industrial enterprises as an example, the road network data is as follows Figure 4 As shown, population data Figure 5 As shown;

[0118] The basic data type for evaluation is raster data. Non-raster data should be preprocessed and uniformly converted into gridded raster data that is convenient for spatial calculation.

[0119] Then calculate the importance index;

[0120] The calculation of the importance index of the carbon sink function of the terrestrial ecosystem selected 2001 and 2020 as the start and end time nodes of the model calculation, and calculated the carbon sink amount with a spatial resolution of 30 meters according to the carbon sink function calculation equation as the evaluation index of the importance of the carbon sink function of the terrestrial ecosystem in Hunan Province, such as Figure 6 As shown;

[0121] The importance index of atmospheric pollutants and greenhouse gas emissions was calculated by selecting 2020 as the base year, and the importance index of atmospheric pollutants and greenhouse gas emissions at 3 km spatial resolution was calculated according to the calculation equations, as follows: Figure 7 and Figure 8 shown.

[0122] Then, the evaluation and grading were carried out, and through calculation, the grid map of the importance evaluation of the carbon sink function of the Hunan ecosystem, and the grid map of the importance evaluation of the emission of atmospheric pollutants and greenhouse gases were obtained respectively;

[0123] In the geographic information system software, use the raster calculator and input the formula "Int([raster data of a certain function] / [the maximum value of raster data of a certain function]×100)" to obtain the normalized importance evaluation raster map;

[0124] Export the raster data attribute table, which records the importance value of each raster pixel, arranges the importance values ​​in descending order, and calculates the cumulative importance value;

[0125] The grid values ​​corresponding to the proportion of the cumulative importance value to the total value of 50% and 80% are used as the dividing points for the evaluation and classification of the importance of ecosystem carbon sink function, atmospheric pollutants and greenhouse gas emissions;

[0126] Using the reclassification tool of geographic information system software, the importance of ecosystem carbon sink function, atmospheric pollutants and greenhouse gas emissions were divided into three levels, namely, extremely important, important and moderately important;

[0127] The importance of ecosystem carbon sink function is as follows Fig. 9 As shown in Figure 2, the importance of atmospheric pollutants and greenhouse gas emissions is as follows: Fig.10 and Fig.11 shown.

[0128] Finally, the greenhouse gas control zones are delineated, and the extremely important and important areas with carbon sink functions in the evaluation are extracted as the priority greenhouse gas protection areas;

[0129] On the basis of demarcating priority protection areas, the important and extremely important areas affected by greenhouse gas emissions and air pollutant emissions are spatially superimposed, and the intersection of the two is taken as the key greenhouse gas control area. The areas outside the priority protection and key control areas are taken as general greenhouse gas control areas, such as Fig.12 shown.

[0130] Example 3

[0131] The greenhouse gas control zoning system based on spatial differential distribution includes:

[0132] The unit accuracy identification module unifies the spatial accuracy of the evaluation work according to the availability of the evaluation parameters of the target area. It has an evaluation parameter and unit accuracy identification program to screen and identify the evaluation index parameters and determine the spatial accuracy of the evaluation work in the target area, providing a basis for determining the subsequent importance evaluation method. According to the availability of the evaluation parameters of the target area, the spatial resolution of the evaluation work is unified as the basic evaluation unit;

[0133] The evaluation and calculation selection module consists of quantitative evaluation methods for the importance of ecosystem carbon sink functions, and the importance of atmospheric pollutants and greenhouse gas emissions. Based on the importance of terrestrial ecosystem carbon sink functions and the status of atmospheric pollutants and greenhouse gas emissions in the target area, combined with the availability of evaluation parameters, quantitative evaluation and calculation methods for the importance of carbon sink functions, and the importance of atmospheric pollutants and greenhouse gas emissions are selected;

[0134] The data collection and processing module collects the data required for the evaluation according to the selected evaluation calculation method and the actual operation requirements of the importance evaluation method, and uniformly converts it into gridded raster data that is convenient for spatial calculation, providing input data for the next step of model operation, such as basic geographic information data, land use data, vegetation cover data, remote sensing images, atmospheric pollutants and greenhouse gas fusion emission inventory data, etc.;

[0135] The importance index calculation module, based on the evaluation calculation method and the data required for the evaluation, inputs the data into the calculation equation of the importance quantitative evaluation method, calculates the importance index of the carbon sink function of the terrestrial ecosystem and the importance index of the impact of atmospheric pollutants and greenhouse gas emissions, as the basis for importance classification;

[0136] The importance assessment grading module evaluates the importance index of terrestrial ecosystem carbon sink function, atmospheric pollutants and greenhouse gas emission impacts based on the cumulative contribution calculation principle and the calculation results, and divides the importance index into three levels: general importance, important and extremely important.

[0137] The zoning demarcation module conducts zoning demarcation based on the assessment results, and divides the corresponding areas into greenhouse gas priority protection areas, greenhouse gas key control areas and greenhouse gas general control areas. The extremely important and important areas with ecosystem carbon sink function in the importance assessment classification results are designated as greenhouse gas priority protection areas, and the extremely important and important areas affected by greenhouse gas and air pollutant emissions are superimposed for analysis. The areas with spatial overlap are identified as greenhouse gas key control areas, and other areas are identified as greenhouse gas general control areas.

[0138] The above-mentioned embodiments only express the preferred implementation methods of the present application, and the descriptions thereof are relatively specific and detailed, but they cannot be understood as limiting the scope of the present application. It should be pointed out that, for ordinary technicians in this field, several modifications, improvements and substitutions can be made without departing from the concept of the present application, and these all belong to the protection scope of the present application.

Claims

1. A greenhouse gas control zoning method based on spatial differential distribution, characterized in that: The following steps are involved: Step 1: Unify the spatial accuracy of the assessment work based on the availability of assessment parameters in the target area; Step 2: Based on the importance of the carbon sink function of the terrestrial ecosystem in the target area and the emission status of atmospheric pollutants and greenhouse gases, combined with the availability of assessment parameters, select the assessment calculation method; Step 3: Collect the data required for the evaluation according to the selected evaluation calculation method; Step 4: Based on the evaluation calculation method and the data required for the evaluation, calculate the importance index of the carbon sink function of the terrestrial ecosystem and the importance index of the impact of atmospheric pollutants and greenhouse gas emissions; Step 5: Based on the cumulative contribution calculation principle and the calculation results, evaluate the importance index of terrestrial ecosystem carbon sink function and the importance index of atmospheric pollutants and greenhouse gas emissions; Step 6: Based on the assessment results, zoning is carried out to divide the corresponding areas into greenhouse gas priority protection areas, greenhouse gas key control areas and greenhouse gas general control areas.

2. The greenhouse gas control zoning method based on spatial differential distribution according to claim 1 is characterized in that: In the step 1, the evaluation parameters of the target area are screened and identified, and the availability of the evaluation parameters after screening and identification is obtained. According to the availability of the evaluation parameters, the spatial accuracy of the target area for carrying out the evaluation work is determined and unified, and the spatial accuracy is used as the basic evaluation unit; The basic assessment units include: the spatial resolution of evaluating the importance of carbon sink function of terrestrial ecosystems, and the spatial resolution of evaluating the importance of the impact of atmospheric pollutants and greenhouse gas emissions.

3. The greenhouse gas control zoning method based on spatial differential distribution according to claim 1 is characterized in that: After collecting the data required for the evaluation in step 3, preprocessing the data required for the evaluation is performed to convert the data required for the evaluation into gridded raster data; The data required for the assessment include basic geographic information data, land use data, vegetation cover data, remote sensing images, industrial enterprise location data, road network data, population data, and fused emission inventory data of atmospheric pollutants and greenhouse gases.

4. The greenhouse gas control zoning method based on spatial differential distribution according to claim 1 is characterized in that: The importance index of the terrestrial ecosystem carbon sink function in step 4 is calculated by characterizing the carbon sink capacity of regional forest ecosystems. The forest ecosystem carbon sink accounting method is as follows: Obtain vegetation coverage and its change types, area, vegetation and soil activity data within the accounting boundary of the target area through investigation and statistics; Combined with the emission factor, the carbon stock change within the boundary is calculated and converted into carbon dioxide equivalent. The carbon dioxide equivalent calculation formula is as follows: Among them, Sink LAND,t is the carbon sink of terrestrial ecosystems within the accounting boundary in year t, is the carbon storage of carbon pool p in vegetation cover and its change type i within the boundary in the initial year t1, is the carbon storage of carbon pool p in vegetation cover and its change type i within the boundary of the accounting period in year t2, is the activity data of carbon pool p in vegetation cover and its change type i within the accounting boundary in year t, is the emission factor of carbon pool p in vegetation cover and its change type i within the accounting boundary in year t, i is vegetation cover and its change type, p is carbon pool type, T is the accounting time interval, t is the accounting year, is a constant for converting carbon into carbon dioxide equivalent.

5. The greenhouse gas control zoning method based on spatial differential distribution according to claim 4 is characterized in that: The importance index of air pollutant and greenhouse gas emissions in step 4 is calculated by characterizing the amount of air pollutant emissions and greenhouse gas emissions in the spatial grid unit. First, the total amount of emissions is calculated, the emission space is allocated, and the emission equivalent is converted. Then, the air pollutants are calculated to obtain the pollutant emission equivalent. Finally, the pollutant emission equivalent is counted to the corresponding grid unit. The calculation formula is as follows: Among them, E AP is the main pollutant equivalent, Q i is the emission of the i-th pollutant, W i is the equivalent value of the ith pollutant.

6. The greenhouse gas control zoning method based on spatial differential distribution according to claim 5 is characterized in that: The emission spatial allocation includes point source processing, mobile source processing and surface source processing. Point source processing is directly positioned on the grid unit according to the longitude and latitude coordinates. Mobile source processing uses the substitute parameter weight method to allocate road network data emissions to the grid. Surface source processing uses the substitute parameter weight method to allocate population data emissions to the grid.

7. The greenhouse gas control zoning method based on spatial differential distribution according to claim 1 is characterized in that: There are three levels of importance index in step 5, including general importance, important and extremely important.

8. The greenhouse gas control zoning method based on spatial differential distribution according to claim 7 is characterized in that: The process of evaluating the importance index in step 5 is as follows: The calculation results are calculated using a grid calculator to obtain normalized importance evaluation grid maps, including importance evaluation grid maps of terrestrial ecosystem carbon sink function, and importance evaluation grid maps of atmospheric pollutants and greenhouse gas emissions; Derive the raster data attribute table from the importance evaluation raster map. The attribute table records the importance value of each raster pixel. The importance values ​​are arranged in descending order and the accumulated importance value is calculated. Two ratio values ​​are preset to obtain the grid values ​​corresponding to the two ratio values ​​of the cumulative importance value in the total value, and the two grid values ​​are used as the dividing points for the evaluation and classification of the importance index of terrestrial ecosystem carbon sink function, atmospheric pollutants and greenhouse gas emissions; The reclassification tool was used in combination with the cut-off points to evaluate and grade the importance index of terrestrial ecosystem carbon sink function, atmospheric pollutants and greenhouse gas emissions.

9. The greenhouse gas control zoning method based on spatial differential distribution according to claim 8 is characterized in that: The specific process of zoning in step 6 is as follows: The extremely important areas and important areas in the grid map of the importance of carbon sink function of terrestrial ecosystems are designated as priority protection areas for greenhouse gases; Overlay analysis is performed on the extremely important areas and important areas in the greenhouse gas and air pollutant emission importance assessment grid map, and the areas with spatial overlap are regarded as key greenhouse gas control areas; Other areas are designated as general greenhouse gas control areas.

10. A greenhouse gas control zoning system based on spatial differential distribution, characterized in that: include: The unit accuracy identification module unifies the spatial accuracy of the assessment work according to the availability of assessment parameters in the target area; The evaluation calculation selection module selects the evaluation calculation method based on the importance of the carbon sink function of the terrestrial ecosystem in the target area and the emission status of atmospheric pollutants and greenhouse gases, combined with the availability of evaluation parameters; The data collection and processing module collects the data required for the evaluation according to the selected evaluation calculation method; Importance index calculation module, which calculates the importance index of terrestrial ecosystem carbon sink function and the importance index of atmospheric pollutants and greenhouse gas emissions based on the evaluation calculation method and the data required for the evaluation; The importance assessment and grading module evaluates the importance index of terrestrial ecosystem carbon sink function, atmospheric pollutants and greenhouse gas emissions impact based on the cumulative contribution calculation principle and combined with the calculation results; The zoning demarcation module carries out zoning demarcation based on the assessment results, dividing the corresponding areas into greenhouse gas priority protection areas, greenhouse gas key control areas and greenhouse gas general control areas.

Images