Spatial optimization method and device for carbon sink and carbon-related ecosystem service protection, electronic equipment and storage medium

Abstract

The embodiment of the invention provides a space optimization method and device for carbon sink and carbon-related ecological system service protection, electronic equipment and a storage medium, and the method comprises the steps: determining a carbon-related ecological system service in a to-be-planned space according to a carbon sink value in the to-be-planned space and the correlation between at least one kind of ecological system service; determining carbon sensitive risk distribution in the to-be-planned space; obtaining ecological characteristics of each cell in the to-be-planned space, wherein the ecological characteristics comprise a carbon sink value, a value of a carbon-related ecological system service and a value of a non-carbon-related ecological system service; and performing aggregation processing on the cells according to the ecological characteristics of each cell and a protection target of a preset planning characteristic, determining a carbon sensitive risk of the aggregated space partition according to the carbon sensitive risk distribution, and determining a target space partition of the space to be planned after space planning according to the carbon sensitive risk. Therefore, according to the scheme, the spatial planning effect can be improved by combining ecological system service and carbon sink benefits.

Description

Technical Field

[0001] The embodiments of the present application relate to the field of land space planning technology, and in particular to a spatial optimization method, device, electronic device and storage medium for protecting carbon sinks and carbon-related ecosystem services. Background Art

[0002] As global climate change intensifies, improving carbon sink benefits has become an important task in national land space planning. However, existing technologies often ignore the coordinated protection of ecosystem services, and are unable to improve carbon sink benefits and maximize multiple ecological benefits. Therefore, how to achieve dual protection of ecosystem services and carbon sink benefits in the process of spatial planning is an urgent problem to be solved in this application. Summary of the invention

[0003] In view of this, embodiments of the present application provide a method, device, electronic device and storage medium for spatial optimization of carbon sinks and carbon-related ecosystem service protection to at least partially solve the above-mentioned problems.

[0004] According to a first aspect of an embodiment of the present application, a spatial optimization method for the protection of carbon sinks and carbon-related ecosystem services is provided, including: determining carbon-related ecosystem services in the space to be planned based on the correlation between the carbon sink value and at least one ecosystem service in the space to be planned; determining the carbon-sensitive risk distribution in the space to be planned; obtaining the ecological characteristics of each cell in the space to be planned, the ecological characteristics including the carbon sink value, the value of the carbon-related ecosystem service and the value of the non-carbon-related ecosystem service; aggregating the cells according to the ecological characteristics of each cell and the protection targets of preset planning characteristics, determining the carbon-sensitive risk of the aggregated spatial partition according to the carbon-sensitive risk distribution, and determining the target spatial partition of the space to be planned after spatial planning according to the carbon-sensitive risk.

[0005] According to a second aspect of an embodiment of the present application, a spatial optimization device for protecting carbon sinks and carbon-related ecosystem services is provided, including: a first determination module, used to determine the carbon-related ecosystem services in the space to be planned based on the correlation between the carbon sink value in the space to be planned and at least one ecosystem service; a second determination module, used to determine the carbon-sensitive risk distribution in the space to be planned; an acquisition module, used to acquire the ecological characteristics of each cell in the space to be planned, the ecological characteristics including the carbon sink value, the value of the carbon-related ecosystem service and the value of the non-carbon-related ecosystem service; a zoning planning module, used to aggregate the cells according to the ecological characteristics of each cell and the preset planning characteristic protection targets, determine the carbon-sensitive risk of the aggregated spatial zoning according to the carbon-sensitive risk distribution, and determine the target spatial zoning after spatial planning according to the carbon-sensitive risk.

[0006] According to the third aspect of an embodiment of the present application, there is provided an electronic device comprising: a processor, a memory, a communication interface and a communication bus, wherein the processor, the memory and the communication interface communicate with each other through the communication bus; the memory is used to store at least one executable instruction, and the executable instruction enables the processor to perform operations corresponding to the spatial optimization method for protecting carbon sinks and carbon-related ecosystem services described in any one of the above items.

[0007] According to a fourth aspect of an embodiment of the present application, a computer storage medium is provided, on which a computer program is stored, and when the program is executed by a processor, any of the above-mentioned spatial optimization methods for protecting carbon sinks and carbon-related ecosystem services is implemented.

[0008] According to the spatial optimization method for protecting carbon sinks and carbon-related ecosystem services provided in each embodiment of the present application, carbon-related ecosystem services are determined based on the correlation between carbon sink values ​​and ecosystem services in the space to be planned, so that ecosystem services are introduced into spatial planning, and the ecological characteristics of each cell in the space to be planned, such as carbon sink values, values ​​of carbon-related ecosystem services, and values ​​of non-carbon-related ecosystems, can be further determined; by determining the distribution of carbon-sensitive risks in the space to be planned, the cells are aggregated by comprehensively considering carbon sink benefits and ecological characteristics in the subsequent spatial planning process, and conditional constraints are implemented on the aggregated cells based on carbon-sensitive risks, so that the planned target space partition can achieve multiple ecological benefits with the lowest carbon-sensitive risk, and then the planned target space can perform dual protection of ecosystem services and carbon sink benefits, thereby improving the effect of spatial planning. BRIEF DESCRIPTION OF THE DRAWINGS

[0009] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments recorded in the embodiments of the present application. For ordinary technicians in this field, other drawings can also be obtained based on these drawings.

[0010] Figure 1 is a flow chart of a spatial optimization method for protecting carbon sinks and carbon-related ecosystem services according to an embodiment of the present application;

[0011] Figure 2 is a schematic diagram of a spatial optimization device for carbon sequestration and carbon-related ecosystem service protection according to an embodiment of the present application;

[0012] Figure 3 It is a schematic diagram of an electronic device according to an embodiment of the present application. DETAILED DESCRIPTION

[0013] Spatial optimization methods for protecting carbon sinks and carbon-related ecosystem services

[0014] The embodiment of the present application provides a spatial optimization method for the protection of carbon sinks and carbon-related ecosystem services. Through this method, dual protection of ecosystem services and carbon sink benefits can be achieved, thereby improving the effect of spatial planning.

[0015] Figure 1 1 is a flow chart of a spatial optimization method for protecting carbon sinks and carbon-related ecosystem services according to an embodiment of the present application. Figure 1 As shown, the spatial optimization method for protecting carbon sinks and carbon-related ecosystem services in this embodiment includes the following steps:

[0016] Step 101: Determine carbon-related ecosystem services in the space to be planned based on the correlation between the carbon sink value in the space to be planned and at least one ecosystem service.

[0017] In some embodiments, ecological data in the space to be planned may be obtained, wherein the ecological data may include solar radiation distribution data, precipitation distribution data, temperature distribution data, vegetation type distribution data, and vegetation index distribution data, and these data may be normalized.

[0018] In some embodiments, the carbon sink distribution of the area to be planned can be determined based on the normalized ecological data in combination with a light energy utilization model, such as a CASA (Carnegie-Ames-Stanford Approach, a remote sensing method based on light energy utilization) model. Among them, the carbon sink can be used to indicate the carbon absorption capacity of vegetation at different locations in the area to be planned. As an example, the difference between the total organic carbon fixed by plants through photosynthesis in the area to be planned and the carbon consumed by plant respiration can be determined as the carbon sink value, or the carbon sink value of the area to be planned can be determined based on information such as vegetation biomass, soil carbon storage, and decomposition rate of dead vegetation in the area to be planned.

[0019] In some embodiments, the spatial distribution of ecosystem services in the area to be planned can be calculated by an ecological model, and the ecosystem services include at least one of regulating services, supporting services, provisioning services, and cultural services.

[0020] In some embodiments, the space to be planned can be divided into several cells, and the area corresponding to each cell corresponds to the area in the space to be planned. Based on the carbon sink distribution and the spatial distribution of ecosystem services in the area to be planned, the carbon sink value and the ecosystem service value in each cell can be obtained through further extraction and analysis. Among them, the ecosystem service value of the area to be planned can be determined by the total amount of the value corresponding to each type of ecological service in the area to be planned, and the value corresponding to each type of ecological service can be determined by the existing ecosystem service evaluation model.

[0021] The correlation between each ecosystem service value and carbon sink value in each cell can be calculated by the Pearson correlation coefficient or the Spearman correlation coefficient to obtain a correlation coefficient, which can measure the dependence between the two variables. The ecosystem service type corresponding to the ecosystem service value with a correlation coefficient greater than the preset coefficient threshold can be determined as a carbon-related ecosystem service. As an example, if the correlation coefficient r ranges from -1 to 1, the larger the correlation coefficient, the stronger the positive correlation synergy between the ecosystem service value and the carbon sink value; the smaller the correlation coefficient, the stronger the negative correlation synergy between the ecosystem service value and the carbon sink value, and 0 can indicate that there is no correlation between the ecosystem service value and the carbon sink value. The preset coefficient threshold can be determined as 0.7, and the ecosystem service corresponding to the ecosystem service value with a correlation coefficient greater than 0.7 can be determined as a carbon-related ecosystem service, that is, the relationship between the ecosystem service value of the carbon-related ecosystem service and the carbon sink value is a positive correlation synergy.

[0022] Ecosystem services other than carbon-related ecosystem services among all ecosystem services can be identified as non-carbon-related ecosystem services.

[0023] Step 102: Determine the carbon sensitivity risk distribution in the space to be planned.

[0024] In some embodiments, the space to be planned may be divided into administrative regions, and the carbon emissions of each administrative region of the space to be planned in at least two years may be obtained.

[0025] In some embodiments, if the number of years obtained is two, the carbon sensitivity risk in each administrative region can be determined based on formula (1).

[0026]

[0027] In formula (1), RI represents the carbon sensitivity risk of the administrative region, CE A represents the carbon emissions of the administrative region in year A, CE B It represents the carbon emission of the administrative region in the year B, where the year B is earlier than the year A.

[0028] In some embodiments, if the number of years obtained is greater than or equal to three, the method for determining carbon sensitivity risk can be based on the following formula (2).

[0029]

[0030] In formula (2), RI represents the carbon sensitivity risk of the administrative region, n represents the number of years, i represents the i-th year, and CE i It represents the carbon emissions of the administrative region in the i-th year. Among them, the year of i-1 is earlier than the year of i.

[0031] In some embodiments, the carbon sensitivity risk distribution in the space to be planned may be determined based on the carbon sensitivity risk corresponding to each administrative region.

[0032] In some embodiments, the method for obtaining the carbon emissions of each administrative area of ​​the planned space in at least two years may be as follows:

[0033] The planning space can be divided based on administrative regions to obtain the construction land area and non-construction land area in each administrative region in each year, wherein the non-construction land area can include at least one of forest land, grassland, cultivated land and water area.

[0034] The direct carbon emissions corresponding to non-construction land in the administrative area can be calculated based on formula (3).

[0035]

[0036] In formula (3), q represents the total direct carbon emissions in the administrative region; i represents each land use type; Q i represents the carbon emission of the i-th land use type; V i represents the area occupied by the i-th land use type, E i It represents the carbon emission factor corresponding to the i-th land use type. The direct carbon emission of the administrative region can be directly determined by summing up the product of the occupied area of ​​each land use type in the non-construction land area of ​​the administrative region and the carbon emission factor of the land use type.

[0037] The indirect carbon emissions of the built-up land area in each administrative region in each year can be calculated based on formula (3).

[0038]

[0039] In formula (3), Q b represents the total indirect carbon emissions in the administrative region, i represents the type of energy consumed in the administrative region, Q iCarbon emissions of the i-th energy type, B i represents the energy consumption of the i-th energy type, E i represents the carbon emission factor of the i-th energy type, F i The standard coal energy conversion factor represents the i-th energy type. The indirect carbon emissions of the administrative region can be determined as the cumulative sum of the carbon emission factor corresponding to each energy type in the construction land area of ​​the administrative region, the energy consumption of each energy type and the standard coal energy conversion factor corresponding to each energy type. Among them, the standard coal energy conversion factor can be determined according to the energy type, and the standard coal energy conversion factor is the conversion factor of each type of energy into standard coal.

[0040] The sum of the direct carbon emissions and indirect carbon emissions of an administrative region can be determined as the total carbon emissions of the administrative region, and the carbon emission distribution of the administrative region in the planned space can be determined based on the total carbon emissions of each administrative region.

[0041] In some embodiments, the method for determining the carbon sensitivity risk distribution of the space to be planned may also be as follows:

[0042] At least one of the land cover change rate, population growth rate, building land ratio, road network density and population density in the planned space can be determined, and the pressure index of the planned space can be obtained after normalizing these data.

[0043] At least one of the proportion of landscape area in the space to be planned, the maximum landscape patch index, the carbon footprint intensity and the per capita carbon footprint can be determined, and these data are normalized to obtain the sensitive indicators of the space to be planned.

[0044] The ecological support index, forest land proportion and vegetation normalization in the planned space can be determined, and these data can be determined as the resilience indicators of the planned space;

[0045] Each pressure indicator, each sensitivity indicator and each elasticity indicator are weighted and superimposed to obtain the carbon sensitivity risk distribution in the planned space. Among them, the pressure indicator, sensitivity indicator and elasticity indicator can also be determined according to actual conditions, and this application does not limit it.

[0046] Step 103: Obtain the ecological characteristics of each cell in the space to be planned.

[0047] In some embodiments, the planned space division cells can be treated as, for example, the planned space division cells can be treated as fishnet grids, or as geometric units or any basic units based on the natural or social characteristics of the planned space, each cell corresponding to an area in the planned space.

[0048] In some embodiments, ecological characteristics corresponding to each cell of the space to be planned may be obtained, and the ecological characteristics may include carbon sink values, values ​​of carbon-related ecosystem services, and values ​​of non-carbon-related ecosystem services.

[0049] Among them, the value of the carbon-related ecosystem service in the cell may be the sum of the ecosystem service values ​​corresponding to the carbon-related ecosystem services in the cell, and the value of the non-carbon-related ecosystem service in the cell may be the sum of the ecosystem service values ​​corresponding to the non-carbon-related ecosystem services in the cell; or, the value of the carbon-related ecosystem service in the cell may be the set of ecosystem service values ​​corresponding to each type of carbon-related ecosystem service in the cell, and the value of the non-carbon-related ecosystem service in the cell may be the set of ecosystem service values ​​corresponding to each type of non-carbon-related ecosystem service in the cell; or, the value of the carbon-related ecosystem service in the cell may be the percentage of the ecosystem service value corresponding to each type of carbon-related ecosystem service in the cell as a percentage of the ecosystem service value of that type of ecosystem service. The value of non-carbon-related ecosystem services in the cell may be a set of percentages of the sum of the ecosystem service values ​​corresponding to each type of non-carbon-related ecosystem services in the cell to the sum of the ecosystem service values ​​corresponding to that type of ecosystem services; or the value of carbon-related ecosystem services in the cell may be the percentage of the sum of the ecosystem service values ​​corresponding to each type of carbon-related ecosystem services in the cell to the sum of the ecosystem service values ​​corresponding to all ecosystem services in the space to be planned, and the value of non-carbon-related ecosystem services in the cell may be the percentage of the sum of the ecosystem service values ​​corresponding to each type of non-carbon-related ecosystem services in the cell to the sum of the ecosystem service values ​​corresponding to all ecosystem services in the space to be planned.

[0050] Step 104: Aggregate the cells according to the ecological characteristics of each cell and the protection targets of the preset planning characteristics, determine the carbon sensitivity risk of the aggregated spatial partitions according to the carbon sensitivity risk distribution, and determine the target spatial partition of the space to be planned after spatial planning according to the carbon sensitivity risk.

[0051] In some embodiments, the protection target of the preset ecological feature can be obtained. The protection target of the preset ecological feature may include a target carbon sink value corresponding to a preset planning partition in the space to be planned, a first target value of a carbon-related ecosystem service, and a second target value of a non-carbon-related ecosystem service. Among them, the target carbon sink value may be the sum of the carbon sink values ​​corresponding to all cells in the preset planning partition in the space to be planned; the first target value may be the first target value of the sum of the ecosystem service values ​​corresponding to the carbon-related ecosystem services corresponding to all cells in the preset planning partition in the space to be planned, and the second target value may be the second target value of the sum of the ecosystem service values ​​corresponding to the non-carbon-related ecosystem services corresponding to all cells in the preset planning partition in the space to be planned; or, the first target value may be the percentage of the sum of the ecosystem service values ​​corresponding to the carbon-related ecosystem services in all cells in the preset planning partition in the space to be planned to the sum of the ecosystem service values ​​corresponding to the carbon-related ecosystem services in the space to be planned, and the second target value may be the percentage of the sum of the ecosystem service values ​​corresponding to the non-carbon-related ecosystem services in all cells in the preset planning partition in the space to be planned to the sum of the ecosystem service values ​​corresponding to the non-carbon-related ecosystem services in the space to be planned.

[0052] In some embodiments, the cells can be aggregated based on the protection target. The aggregation method can be implemented based on the simulated annealing algorithm. Based on the carbon-sensitive risk as a constraint, the space to be planned is divided into several spatial partitions. The first carbon sink value, the first ecosystem service value, and the second ecosystem service value in the spatial partition obtained after aggregation are obtained. Among them, the first carbon sink value is used to indicate the sum of the carbon sink values ​​corresponding to all cells in the spatial partition, the first ecosystem service value is used to indicate the sum of the values ​​of carbon-related ecosystem services corresponding to all cells in the spatial partition, and the second ecosystem service value is used to indicate the sum of the values ​​of non-carbon-related ecosystem services corresponding to all cells in the spatial partition.

[0053] If the difference between the first carbon sink value and the target carbon sink value is within the first threshold, the difference between the first ecosystem service value and the first target value is within the second threshold, and the difference between the second ecosystem service value and the second target value is within the third threshold, the spatial partition is determined as a candidate spatial partition. The first threshold, the second threshold, and the third threshold can be determined according to actual conditions.

[0054] Among them, the method for aggregating cells can be: the number of target space partitioning schemes generated by the simulated degradation algorithm and the number of iterations of each partitioning scheme can be set, spatial iterative operations can be performed, and finally the space partitioning scheme with the smallest objective function is selected from all partitioning schemes as the partitioning scheme for the selected space partition.

[0055] As an example, the objective function for aggregating cells according to the ecological characteristics of each cell and the preset planning characteristic target quantity can be shown as formula (4).

[0056]

[0057] In formula (4), ∑ PUs Cost represents the cost factor of spatial partition, that is, the carbon sensitivity risk of the spatial partition corresponding to the aggregated cell, which can be determined by the carbon sensitivity risk distribution; the edge length adjustment coefficient BLM represents the boundary factor, ∑ PUs Boundary represents the total boundary length corresponding to the spatial partition, FeaturePenalty represents the ecological feature, and FPF represents the penalty value. Among them, BLM and FPF can be determined through pre-training. f represents the objective function value corresponding to the planned spatial partition. The smaller the f value, the better the plan corresponding to the spatial partition. A preset threshold can be set for f, and the spatial partition corresponding to the plan that satisfies the formula f is determined as the candidate spatial partition.

[0058] In some embodiments, the zoning scheme of the candidate spatial partition obtained after each iteration can be obtained, and the value of the carbon sensitive risk of the candidate spatial partition in each zoning scheme is determined according to the carbon sensitive risk distribution in the planned space, and the candidate spatial partition with the smallest carbon sensitive risk value in the zoning scheme is determined as the target spatial partition. The optimal zoning scheme corresponding to the target spatial partition is spatially visualized in the geographic information system to complete the spatial planning based on carbon sink and carbon-related ecosystem service protection.

[0059] The preset planning spatial partitions can include three, namely ecological protection zone, ecological buffer zone and sustainable development zone. As an example, if the protection target of the preset planning feature is as shown in Table 1, the protection target of the preset planning feature of the ecological protection zone includes a target carbon sink value of 75%, a first target value of 75%, and a second target value of 0%. The three spatial partitions obtained after the aggregation of the cells must all meet the protection target of the preset planning feature, and the carbon sensitivity risk corresponding to each partition is small enough.

[0060] Table 1 Planning characteristics and protection targets of pre-set spatial partitions

[0061] Preset planning features Ecological protection area Ecological buffer zone Sustainable Development Zone Target carbon sink value 75％ 25％ 0％ First target value 75％ 25％ 0％ Second target value 0％ 75％ 25％

[0062] In an embodiment of the present application, through the correlation between the carbon sink value and at least one ecosystem service in the space to be planned, the carbon-related ecosystem services and non-carbon-related ecosystem services in the space to be planned can be determined, so that the ecosystem service value and carbon sink value of each cell in the space to be planned can be accurately evaluated; in determining the carbon-sensitive risk distribution in the space to be planned, the carbon-sensitive risk distribution of the space to be planned is determined by changing the sum of indirect carbon emissions from construction land areas and direct carbon emissions from non-construction land areas in each administrative region for at least two years, or the carbon-sensitive risk distribution of the space to be planned is determined by weighted superposition of pressure indicators, sensitivity indicators and elasticity indicators in the space to be planned, so that the carbon-sensitive risk degree of different areas in the space to be planned can be determined; by aggregating the cells according to the ecological characteristics of each cell and the protection targets of the preset planning characteristics, and iteratively determining the optimal target spatial partition from multiple schemes of the aggregated spatial partition based on the carbon-sensitive risk distribution, the spatial planning can improve the dual protection of carbon sink benefits and ecosystems, and improve the spatial planning effect.

[0063] Spatial optimization device for carbon sinks and protection of carbon-related ecosystem services

[0064] Corresponding to the above-mentioned spatial optimization method embodiment for protecting carbon sinks and carbon-related ecosystem services, Figure 2 A schematic diagram of a spatial optimization device for protecting carbon sinks and carbon-related ecosystem services according to an embodiment of the present application is shown. Figure 2 As shown, the spatial optimization device 200 for protecting carbon sinks and carbon-related ecosystem services includes:

[0065] A first determination module 201 is used to determine the carbon-related ecosystem services in the space to be planned according to the correlation between the carbon sink value in the space to be planned and at least one ecosystem service;

[0066] The second determination module 202 is used to determine the carbon sensitivity risk distribution in the space to be planned;

[0067] An acquisition module 203 is used to acquire ecological characteristics of each cell in the planned space, where the ecological characteristics include carbon sink values, values ​​of carbon-related ecosystem services, and values ​​of non-carbon-related ecosystem services;

[0068] The zoning planning module 204 is used to aggregate the cells according to the ecological characteristics of each cell and the protection targets of the preset planning characteristics, determine the carbon sensitivity risks of the aggregated spatial partitions according to the carbon sensitivity risk distribution, and determine the target spatial partitions of the space to be planned after spatial planning according to the carbon sensitivity risks.

[0069] It should be noted that the spatial optimization device for protecting carbon sinks and carbon-related ecosystem services in this embodiment is used to implement the corresponding spatial optimization method for protecting carbon sinks and carbon-related ecosystem services in the aforementioned method embodiments, and has the beneficial effects of the corresponding method embodiments, which will not be repeated here.

[0070] Electronic devices

[0071] Figure 3 A schematic diagram of the structure of an electronic device according to an embodiment of the present application is shown. The specific embodiment of the present application does not limit the specific implementation of the electronic device.

[0072] like Figure 3 As shown, the electronic device may include: a processor (processor) 302 , a communication interface (Communications Interface) 304 , a memory (memory) 306 , and a communication bus 308 .

[0073] in:

[0074] The processor 302 , the communication interface 304 , and the memory 306 communicate with each other via the communication bus 308 .

[0075] The communication interface 304 is used to communicate with other electronic devices or servers.

[0076] The processor 302 is used to execute the program 310, and specifically can execute the relevant steps in the above-mentioned embodiment of the spatial optimization method for protecting carbon sinks and carbon-related ecosystem services.

[0077] Specifically, the program 310 may include program codes, which include computer operation instructions.

[0078] The processor 302 may be a central processing unit (CPU), or an application-specific integrated circuit (ASIC), or one or more integrated circuits configured to implement the embodiments of the present application. The one or more processors included in the smart device may be processors of the same type, such as one or more CPUs; or may be processors of different types, such as one or more CPUs and one or more ASICs.

[0079] The memory 306 is used to store the program 310. The memory 306 may include a high-speed RAM memory, and may also include a non-volatile memory (non-volatile memory), such as at least one disk memory.

[0080] The program 310 may be specifically configured to enable the processor 302 to perform the following operations:

[0081] In an optional implementation, the program 310 is also used to enable the processor 302 to implement each step in the program 310. The specific implementation can refer to the corresponding steps and corresponding descriptions of the apparatus in the above-mentioned spatial optimization method embodiment for protecting carbon sinks and carbon-related ecosystem services, which will not be repeated here. Those skilled in the art can clearly understand that for the convenience and simplicity of description, the specific working process of the above-described devices and modules can refer to the corresponding process description in the above-mentioned method embodiment, which will not be repeated here.

[0082] An embodiment of the present application also provides a computer program product, including computer instructions, which instruct a computing device to perform operations corresponding to any of the spatial optimization methods for protecting carbon sinks and carbon-related ecosystem services in the above-mentioned multiple method embodiments.

[0083] It should be pointed out that, according to the needs of implementation, the various components / steps described in the embodiments of the present application can be split into more components / steps, or two or more components / steps or partial operations of components / steps can be combined into new components / steps to achieve the purpose of the embodiments of the present application.

[0084] The above-mentioned method according to the embodiment of the present application can be implemented in hardware, firmware, or implemented as software or computer code that can be stored in a recording medium (such as a CD ROM, RAM, floppy disk, hard disk or magneto-optical disk), or implemented as a computer code originally stored in a remote recording medium or a non-temporary machine-readable medium downloaded through a network and to be stored in a local recording medium, so that the method described herein can be stored in such software processing on a recording medium using a general-purpose computer, a dedicated processor or programmable or dedicated hardware (such as an ASIC or FPGA). It can be understood that a computer, a processor, a microprocessor controller or programmable hardware includes a storage component (e.g., RAM, ROM, flash memory, etc.) that can store or receive software or computer code, and when the software or computer code is accessed and executed by a computer, a processor or hardware, the spatial optimization method for protecting carbon sinks and carbon-related ecosystem services described herein is implemented. In addition, when a general-purpose computer accesses code for implementing the spatial optimization method for protecting carbon sinks and carbon-related ecosystem services shown herein, the execution of the code converts the general-purpose computer into a dedicated computer for executing the spatial optimization method for protecting carbon sinks and carbon-related ecosystem services shown herein.

[0085] Those of ordinary skill in the art will appreciate that the units and method steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of the embodiments of the present application.

[0086] The above implementation methods are only used to illustrate the embodiments of the present application, and are not limitations on the embodiments of the present application. Ordinary technicians in the relevant technical field can make various changes and modifications without departing from the spirit and scope of the embodiments of the present application. Therefore, all equivalent technical solutions also belong to the scope of the embodiments of the present application. The scope of patent protection of the embodiments of the present application should be limited by the claims.

Claims

1. A spatial optimization method for protecting carbon sinks and carbon-related ecosystem services, characterized in that: The method comprises: Determining carbon-related ecosystem services in the space to be planned based on the correlation between the carbon sink value and at least one ecosystem service in the space to be planned; Determining the distribution of carbon-sensitive risks in the space to be planned; Obtaining ecological characteristics of each cell in the space to be planned, wherein the ecological characteristics include a carbon sink value, a value of the carbon-related ecosystem service, and a value of a non-carbon-related ecosystem service; The cells are aggregated according to the ecological characteristics of each cell and the protection targets of the preset planning characteristics, the carbon sensitivity risks of the aggregated spatial partitions are determined according to the carbon sensitivity risk distribution, and the target spatial partitions of the space to be planned after spatial planning are determined according to the carbon sensitivity risks.

2. The method according to claim 1, characterized in that The step of determining the carbon-related ecosystem service in the space to be planned according to the correlation between the carbon sink value in the space to be planned and at least one ecosystem service comprises: Dividing the space to be planned into cells, where areas corresponding to the cells correspond one-to-one to areas in the space to be planned; determining a correlation coefficient corresponding to a correlation between a carbon sink value in the cell and an ecosystem service value corresponding to the ecosystem service; The ecosystem service corresponding to the ecosystem service value whose correlation coefficient is greater than a preset coefficient threshold is determined as the carbon-related ecosystem service.

3. The method according to claim 1, characterized in that The determining of the carbon sensitive risk distribution in the space to be planned includes: Obtain the carbon emissions of the administrative area of ​​the planned space in at least two years; If the number of years is two, the carbon sensitivity risk of the administrative area is determined based on the following formula; In the above formula, RI represents carbon sensitivity risk, CE A represents the carbon emissions of the administrative region in year A, CE B It represents the carbon emission of the administrative region in the year B, where the year B is earlier than the year A.

4. The method according to claim 3, characterized in that If the number of years is greater than or equal to three, the carbon sensitivity risk of the administrative region is determined based on the following formula: In the above formula, n represents the total number of years, i represents the i-th year, and CE i represents the carbon emissions of the administrative region in the i-th year, and s represents the carbon sensitivity risk of the administrative region; wherein the year i-1 is earlier than the year i; The carbon sensitivity risk distribution in the space to be planned is determined according to the carbon sensitivity risk of the administrative area.

5. The method according to claim 3, characterized in that: The obtaining of carbon emissions of the administrative area of ​​the space to be planned in at least two years includes: Divide the space to be planned based on the administrative region, and obtain the construction land area and non-construction land area in the administrative region in each year; Determine the indirect carbon emissions of the administrative area by summing up the product of the carbon emission factor corresponding to each energy type in the construction land area, the energy consumption of the energy type and the standard coal energy conversion factor corresponding to the energy type; The cumulative sum of the products of the occupied area of ​​each land use type in the non-construction land area and the carbon emission factor of the land use type is determined as the direct carbon emission of the administrative area; wherein the land type in the non-construction land area includes at least one of forest land, grassland, cultivated land and water area; The sum of the indirect carbon emissions and the direct carbon emissions of the administrative region is determined as the carbon emissions of the administrative region in each year.

6. The method according to claim 1, characterized in that Determine the carbon sensitivity risk distribution in the planned space, including: Determine the pressure index of the space to be planned according to at least one of the land cover change rate, population growth rate, building land ratio, road network density and population density in the planned space; Determine a sensitive index of the space to be planned according to at least one of the proportion of landscape area, the maximum landscape patch index, the carbon footprint intensity and the per capita carbon footprint in the space to be planned; Determining the resilience index of the administrative region according to at least one of the ecological support index, forest land proportion and vegetation normalization in the planned space; The pressure index, the sensitivity index and the elasticity index are weighted and superimposed to obtain the carbon sensitivity risk distribution corresponding to the space to be planned.

7. The method according to claim 1, characterized in that The unit cells are aggregated according to the ecological characteristics of each unit cell and the protection target of the preset planning characteristics, the carbon sensitivity risk of the aggregated spatial partition is determined according to the carbon sensitivity risk distribution, and the target spatial partition of the space to be planned after spatial planning is determined according to the carbon sensitivity risk, including: Obtaining a protection target for the preset planning feature, the protection target for the preset planning feature comprising a target carbon sink value corresponding to all cells in a preset planning partition in the to-be-planned space, a first target value for the sum of ecosystem service values ​​corresponding to carbon-related ecosystem services, and a second target value for the sum of ecosystem service values ​​corresponding to non-carbon-related ecosystem services; Aggregating the cells to obtain a first carbon sink value, a first ecosystem service value, and a second ecosystem service value within the spatial partition obtained after aggregation; wherein the first carbon sink value is used to indicate the sum of the carbon sink values ​​corresponding to all cells in the spatial partition, the first ecosystem service value is used to indicate the sum of the values ​​of carbon-related ecosystem services corresponding to all cells in the spatial partition, and the second ecosystem service value is used to indicate the sum of the values ​​of non-carbon-related ecosystem services corresponding to all cells in the spatial partition; If the difference between the first carbon sink value and the target carbon sink value is within a first threshold, the difference between the first ecosystem service value and the first target value is within a second threshold, and the difference between the second ecosystem service value and the second target value is within a third threshold, the spatial partition is determined as a candidate spatial partition; Determine the total amount of carbon sensitive risk of the selected spatial partition according to the carbon sensitive risk distribution; The candidate spatial partition with the smallest total amount of carbon sensitive risk is determined as the target spatial partition.

8. A spatial optimization device for protecting carbon sinks and carbon-related ecosystem services, characterized in that: include: A first determination module is used to determine the carbon-related ecosystem services in the space to be planned according to the correlation between the carbon sink value in the space to be planned and at least one ecosystem service; A second determination module determines the carbon sensitivity risk distribution in the space to be planned; An acquisition module, used for acquiring ecological characteristics of each cell in the space to be planned, wherein the ecological characteristics include carbon sink value, value of carbon-related ecosystem services and value of non-carbon-related ecosystem services; A zoning planning module is used to aggregate the cells according to the ecological characteristics of each cell and the preset planning characteristic protection targets, determine the carbon sensitivity risk of the aggregated spatial partition according to the carbon sensitivity risk distribution, and determine the target spatial partition of the to-be-planned space after spatial planning according to the carbon sensitivity risk.

9. An electronic device, comprising: A processor, a memory, a communication interface and a communication bus, wherein the processor, the memory and the communication interface communicate with each other via the communication bus; The memory is used to store at least one executable instruction, and the executable instruction enables the processor to perform operations corresponding to the spatial optimization method for protecting carbon sinks and carbon-related ecosystem services according to any one of claims 1 to 7.

10. A computer storage medium having a computer program stored thereon, wherein when the program is executed by a processor, the spatial optimization method for protecting carbon sinks and carbon-related ecosystem services as described in any one of claims 1 to 7 is implemented.

Images