Wetland credit enhancement quantification method fusing macro remote sensing and micro engineering CAD data

Abstract

The application discloses a wetland credit promotion quantification method fusing macro remote sensing and micro engineering CAD data, belongs to the field of ecological environment, and comprises the following steps: determining the project range and measure type of wetland ecological restoration; identifying the ecological patches of forest land, grassland and water body in the engineering range before and after the implementation of ecological restoration and determining the area size of the ecological patches; identifying and determining the ecosystem service type for different ecological patch types, and calculating the ecosystem service; summing the ecosystem services; determining the total sum of the ecosystem services of various types of ecological patches in the engineering range of wetland ecological restoration; calculating the change of the ecosystem services before and after the implementation of the wetland ecological restoration engineering; fusing the wetland credit grid calculation method of macro remote sensing data and engineering CAD data; and finally determining the promotion amount of the wetland credit. The application adopts the above method, which is helpful to enhancing the implementation of the wetland ecological restoration project and the accuracy of the evaluation decision, and promoting the input intensity of the wetland ecological restoration.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of ecological environment, in particular to a wetland credit promotion quantification method fusing macro remote sensing and micro engineering CAD data. BACKGROUND

[0002] Wetland ecosystems provide important ecosystem services such as water purification and climate regulation. It is estimated that global wetlands provide at least 47.4 trillion dollars in ecosystem service value annually. However, with the continuous development and utilization of human activities on wetland ecosystems, global wetland ecosystem services are facing a dramatic decrease. In response to this challenge, countries have explored credit trading mechanisms based on wetland ecological restoration, aiming to provide incentives and support for wetland ecological restoration through market-based means.

[0003] Wetland credit refers to a unified measure of the services provided by wetland ecosystems. After carrying out ecological restoration activities or measures such as wetland conservation, enhancement, restoration, and new construction, the services provided by wetland ecosystems will be improved. Therefore, the subject of credit trading based on wetland ecological restoration is "wetland credit", and the essence of the transaction is the services of wetland ecosystems.

[0004] Establishing a comprehensive quantitative method for assessing wetland credit is the first step in the credit trading mechanism based on wetland ecological restoration. However, existing technologies mainly rely on qualitative or semi-quantitative evaluation methods, resulting in three systematic defects: first, the non-structured characteristics of the evaluation index system cause the credit value accounting dimension to be incomplete; second, the discrete evaluation method leads to the problem of credit distortion caused by overestimation or underestimation of credit measurement results; third, the lack of uniformity in the measurement benchmark leads to the problem of fuzzy credit rights, which may cause the risk of repeated trading of credit. Moreover, these methods consider fewer types of wetland ecosystem services. Currently, there is no wetland credit quantification method that combines macro remote sensing data and micro engineering data and considers multiple types of ecosystem services. This greatly limits the promotion and implementation of market-based ecological compensation mechanisms based on wetland ecological restoration. Therefore, it is necessary to establish a fully quantitative method to more accurately quantify the credit promotion of wetland ecosystems after ecological restoration to determine the amount of "credit". SUMMARY

[0005] The purpose of the present application is to provide a wetland credit promotion quantification method fusing macro remote sensing and micro engineering CAD data to solve the problems in the background art.

[0006] To achieve the above-mentioned purpose, the present application provides a wetland credit promotion quantification method fusing macro remote sensing and micro engineering CAD data, comprising the following steps:

[0007] S1, determining the project scope and measure type of wetland ecological restoration, determining the boundary and scope of the project according to the planning and design drawing of wetland ecological restoration; according to the specific restoration measures and implementation schedule of the project, the time factor of ecological restoration is determined;

[0008] S2, establishing an ecological patch recognition system, recognizing the ecological patches of forest land, grassland and water body within the project scope before and after the implementation of ecological restoration and determining the area size of the ecological patches;

[0009] S3, for different types of ecological patches, the type of ecosystem service is identified and determined, and the ecosystem service of each type of ecological patch is calculated respectively;

[0010] S4, the ecosystem services of each type of ecological patch are added;

[0011] S5, determining the total ecosystem service of each type of ecological patch within the scope of the wetland ecological restoration project;

[0012] S6, calculating the change of the ecosystem service before and after the implementation of the wetland ecological restoration project;

[0013] S7, using CAD engineering drawings and remote sensing image maps, multi-angle and multi-level calculation is performed on the ecological patches, and grid processing and weighted superposition operation are performed combined with the characteristic attributes of each ecological patch;

[0014] S8, finally determining the amount of wetland credit improvement.

[0015] Preferably, the ecosystem services calculated in S3 are as follows:

[0016] a. Carbon sequestration and oxygen release BC:

[0017] ;

[0018] In the formula, is the biological carbon sequestration (sej / yr); is the carbon sequestration of the first ecological patch type per unit area (g C / m 2 / yr); is the average turnover time of the first ecological patch type (yr); is the area of the ecological patch type (m 2 ); is the energy value conversion rate of the ecological patch type (sej / g);

[0019] b. Constructed soil (FS):

[0020] ;

[0021] In the formula, It is about building soil (sej / yr), and building soil organic matter ( ) and minerals ( The sum of the two parts;

[0022] c. Add bottom sediment (IS):

[0023] ;

[0024] ;

[0025] in, It is the first Increased sediment (sej / yr) in different ecological patch types (such as rivers, lakes and swamps); For the first Annual organic matter deposition (g C / m³) of each ecological patch type 2 / yr); This represents the percentage of organic matter absorbed by wetland plants, with a value of 0.78. It is the first The conversion factor (kcal / g) of g to kcal in each ecological patch type. It is the conversion factor from kcal to J, with a value of 4186 J / kcal; No. Area of ​​each ecological patch type (m²) 2 ); It is the first Energy conversion rate of deposited organic matter in each ecological patch type (sej / J); Organic sedimentary material accounts for the majority. The proportion of NPP for each ecological patch type was 30.37%. It is the first NPP levels (g C / m³) for each ecological patch type 2 / yr).

[0026] d. Groundwater recharge (RW):

[0027] ;

[0028] in, It is the first Recharged groundwater for each ecological patch type (sej / yr); For the first Average annual precipitation (m / yr) within the range of each ecological patch type. For the first Area of ​​each ecological patch type (m²) 2 ); The density of water (kg / m³)3 is the recharge coefficient of precipitation; 1000 is the conversion coefficient from kg to g; is the energy conversion rate of precipitation infiltrating into groundwater, because it is groundwater recharge, so the energy conversion rate of rainwater (sej / g) is used here.

[0029] e, Purification of the atmosphere (PA):

[0030] ;

[0031] In the formula, is the purification of the atmosphere (sej / yr) of the first ecological patch type, which is the sum of the reduction in human health loss ( ) and the reduction in ecosystem loss ( );

[0032] f, Purification of water quality (PW):

[0033] ;

[0034] In the formula, is the purification of water quality (sej / yr) of the first ecological patch type, which is the sum of the reduction in human health loss ( ) and the reduction in ecological resource loss ( );

[0035] g, Soil retention (RS):

[0036] ;

[0037] ;

[0038] In the formula, is the soil retention (sej / yr) of the first ecological patch type; is the total annual soil retention of the first i ecological patch type (t / yr); 106 is the conversion coefficient from t to g; is the percentage of organic matter content in the soil of the first ecological patch type; is the conversion coefficient from g to kcal, which is 5.4 (kcal / g); is the conversion coefficient from mg to kg, which is ; is the conversion coefficient from kcal to J, which is 4186 J / kcal; ​Energy conversion rate of soil (sej / j); Potential soil erosion modulus (t / km 2 / yr) of the th ecological patch type; Real soil erosion modulus (t / km 2 / yr) of the th ecological patch type; Area (m 2 ) of the th ecological patch type;

[0039] h, regulating local temperature and humidity (RM):

[0040] ;

[0041] Wherein, Regulating local temperature and humidity (sej / yr) of the th ecological patch type; Annual evapotranspiration (g / m 2 / yr) of the th ecological patch type; Area (m 2 ) of the th ecological patch type; Energy conversion rate of evapotranspiration (sej / g);

[0042] i, regulating climate (RC):

[0043] ;

[0044] Wherein, Regulating climate of the th ecological patch type, which is the sum of the two parts of carbon sink reducing the impact of climate change on human health ( ) and carbon sink reducing the impact of ecosystem quality ( ).

[0045] Preferably, the formula for calculating organic matter ( ) and mineral matter ( ) is as follows:

[0046] Organic matter ( )

[0047] ;

[0048] Wherein Building soil organic matter (sej / yr); Renewable resources corresponding to the energy of the ecological patch type (sej / yr); is the proportion of vegetation litter in biomass (g / g, %); is the carbon content in vegetation litter (g / g, %);

[0049] mineral substance (s ) in soil:

[0050] ;

[0051] In the formula, is the construction of soil mineral substance (sej / yr); is the first ecological patch type soil in the first mineral substance in the total mineral substance percentage (%); is the soil bulk density of the first ecological patch type (g / cm3); is the soil thickness of the first ecological patch type (cm); is the area of the first ecological patch type (m 2 ); is the percentage of soil mineral substance in total soil weight (%), taking 95%; 10000 is the conversion coefficient from m 2 to cm 2 ; is the turnover time of the first mineral substance in the first ecological patch type soil (years); is the energy conversion rate of the first mineral substance in the first ecological patch type soil (sej / g).

[0052] The preferred amount of reduction of human health loss (Y ) and the amount of reduction of ecosystem loss (Y ) are calculated as follows:

[0053] The amount of reduction of human health loss (Y ):

[0054] ;

[0055] ;

[0056] In the formula, is the amount of reduction of human health loss after purifying atmospheric pollutants in the first ecological patch type (sej / yr); is the first ecological patch type purifying the first mineral substance (sej / yr);Capacity for various air pollutants (kg / ha / yr); For the first Area of ​​each ecological patch type (m²) 2 ); 0.0001 is m 2 The coefficient converted to ha; It is the first The impact factor of each air pollutant in the Eco-indicator 99 assessment framework, i.e., per kg of the first pollutant. Disability life-adjusted years (person-years / kg) caused by various air pollutants in humans. The energy value (sej / person) corresponding to the per capita total regional health expenditure; Total regional health expenditure (yuan); The regional energy value-to-currency ratio (SEJ / Yuan); The number of permanent residents in the region (in people);

[0057] Ecosystem loss reduction ( ):

[0058]

[0059] ;

[0060] In the formula, Indicates the first The energy value (sej / yr) of the reduction in ecosystem loss after the purification of air pollutants by each ecological patch type; For the first The ability of each ecological patch type to purify the jth type of air pollutant (kg / ha / yr); Indicates being subject to the The impact of various air pollutants on the potential extinction rate of species (PDF×m) 2 ×yr); 0.0001 is m 2 The conversion coefficient to ha; Indicates the first Energy required (sej / yr) for species in the region where each ecological patch type is located; This represents the total number of ecological patch types.

[0061] Preferred reduction in human health loss ( ) and reduction in ecological resource loss ( The calculation formula is as follows:

[0062] The reduction in human health loss ( ):

[0063] ;

[0064] in, It is about reducing harm to the human body (sej / yr); It is the first Purification of each type of ecological patch The ability of water pollutants to be converted (mg / kg); 0.001 is the conversion of kg into g The conversion coefficient; It is the first NPP (g C / m³) for each ecological patch type 2 / yr); It is the first Ecological patch types (m) 2 ); It is the first Disability life adjustment years (CAP) caused by water pollutants yr / kg) The coefficient for converting mg to kg is given, and its value is [value missing]. ; The energy value (sej / person) corresponding to the per capita total regional health expenditure; It is the first Turnover time (yr) for Class II water pollutants, with a value of 1000yr;

[0065] Reduction in ecological resource loss ( ):

[0066] ;

[0067] in, It is about reducing ecosystem damage (sej / yr); It is the first Purification of each type of ecological patch The ability of water pollutants to be converted (mg / kg); 0.001 is the conversion of kg into g The conversion coefficient; It is the first NPP (gC / m³) for each ecological patch type 2 / yr); The coefficient for converting mg to kg is given, and its value is [value missing]. ; It is the first Proportion of potential species extinctions caused by Class II water pollutants (%×m) 2 ×yr×kg-1); It is the first The energy required by species in the region where each ecological patch type is located; It is the first The turnover time (yr) of pollutants in water bodies is taken as 1000yr.

[0068] Preferably, carbon sequestration reduces the impact of climate change on human health. ) and the impact of carbon sinks on ecosystem quality ( The calculation formula for ) is as follows:

[0069] Carbon sequestration reduces the impact of climate change on human health. ):

[0070] ;

[0071] In the formula, It is the first The impact of climate change on human health due to reduced carbon sequestration in several ecological patch types (sej / yr); It is the first on a global scale The impact of each ecological patch type on greenhouse gases per unit area Average fixed amount (fixed amount (g / m) 2 / yr); 0.001 is g Converted to kg; It is per kg of the first Disability life-adjusted years due to greenhouse gases (person-years / kg); For the first Greenhouse gas life cycle (yr); It is the per capita total regional health expenditure (sej / cap); It is the first after VF correction Area (ha) of each ecological patch type;

[0072] Carbon sinks reduce the impact on ecosystem quality. ):

[0073] ;

[0074] In the formula, It is the first The impact of carbon sinks of individual ecological patch types on ecosystem quality (sej / yr); It is the first on a global scale The impact of each ecological patch type on greenhouse gases per unit area Average fixed amount (fixed amount (g / m) 2 / yr); 0.001 is g Converted to kg; For the first Greenhouse gas life cycle (yr); It is the first The proportion of species potential extinction caused by greenhouse gases; It is the first The energy required by species in the region where each ecological patch type is located.

[0075] Preferably, the content of S4 is as follows:

[0076] Total ecosystem services of aquatic patches:

[0077] ;

[0078] in, This represents the sum of ecosystem services provided by a water patch. Indicates biological carbon sequestration in water patches. This indicates an increase in organic matter in the bottom sediment, indicating an increase in water patch formation. This indicates the groundwater supply to the water patches. This indicates the regulation of local temperature and humidity by water patches. Indicates the purified water quality of water patches. This indicates the climate moderating effect of water patches. The unit for each ecosystem service is sej / yr.

[0079] Total ecosystem services of grassland patches:

[0080] ;

[0081] in, This represents the sum of ecosystem services provided by grassland patches. Indicates biological carbon sequestration in grassland patches. Indicates the soil that forms grassland patches. This indicates the groundwater supply to the grassland patches. This indicates the regulation of local temperature and humidity by grassland patches. This indicates that the grassland patches purify the atmosphere. This indicates a reduction in soil erosion in grassland patches. This indicates the climate moderating effect of grassland patches. The unit for each ecosystem service is sej / yr.

[0082] Total ecosystem services of forest patches:

[0083] ;

[0084] in, This represents the sum of ecosystem services provided by forest patches. Indicates biological carbon sequestration in forest patches. Indicates the soil that forms woodland patches. This indicates the groundwater replenishment of woodland patches. This indicates the role of forest patches in regulating local temperature and humidity. This indicates that woodland patches purify the atmosphere. This indicates a reduction in soil erosion in forest patches. This indicates the climate moderating effect of forest patches. The unit for each ecosystem service is sej / yr.

[0085] Preferably, the content of S5 is as follows:

[0086] = ;

[0087] In the formula, Representing the Ecosystem services of wetland ecological restoration projects Representing the Ecosystem services of water patches in wetland ecological restoration projects Representing the Ecosystem services of grassland patches in wetland ecological restoration projects Representing the Ecosystem services of forest patches in wetland ecological restoration projects.

[0088] Preferably, the content of S6 is as follows:

[0089] ;

[0090] In the formula, Representing the Changes in ecosystem services of a wetland ecological restoration project. Representing the Ecosystem services following the implementation of wetland ecological restoration projects Representing the Ecosystem services prior to the implementation of a wetland ecological restoration project.

[0091] Preferably, the wetland credit enhancement limit in S8 depends on the improvement of ecosystem services before and after the implementation of the wetland ecological restoration project, while also considering the time risk factor and weight of the wetland ecological restoration project. The calculation formula is as follows:

[0092] ;

[0093] In the formula, It refers to the first The impact of wetland ecological restoration projects on wetland creditworthiness. Representing the Changes in ecosystem services of a wetland ecological restoration project. Representing the The time risk factors of wetland ecological restoration projects Indicates the first The weights of each factor in a wetland ecological restoration project.

[0094] Therefore, the wetland credit enhancement quantification method that integrates macroscopic remote sensing and microscopic engineering CAD data, as described above, has the following beneficial effects:

[0095] (1) Take into account the types of wetland ecological restoration measures, the types of wetland ecological restoration patches, the various services of the wetland ecosystem, and the changes in services before and after ecological restoration.

[0096] (2) It helps to enhance the accuracy of wetland ecological restoration project implementation and evaluation decisions, and at the same time helps to provide a basis for credit transactions for wetland ecological restoration, thereby promoting investment in wetland ecological restoration.

[0097] The technical solution of the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. Attached Figure Description

[0098] Figure 1 This is a flowchart illustrating the wetland credit enhancement quantification method that integrates macroscopic remote sensing and microscopic engineering CAD data according to the present invention.

[0099] Figure 2 This is a schematic diagram illustrating the multi-level calculation of the wetland credit enhancement quantification method that integrates macroscopic remote sensing and microscopic engineering CAD data according to the present invention. Detailed Implementation

[0100] The following detailed description of embodiments of the invention provided in the accompanying drawings is not intended to limit the scope of the claimed invention, but merely to illustrate selected embodiments of the invention. All other embodiments obtained by those skilled in the art based on the embodiments of the invention without inventive effort are within the scope of protection of the invention.

[0101] Please see Figure 1 A quantitative method for improving wetland credit by integrating macroscopic remote sensing and microscopic engineering CAD data includes the following steps:

[0102] S1. Determine the project scope and types of measures for wetland ecological restoration. Based on the planning and design drawings for wetland ecological restoration, determine the project boundaries and scope. Based on the specific restoration measures and implementation schedule of the project, determine the corresponding time factors for ecological restoration.

[0103] Given the complexity of wetland ecosystems, at a macro scale, comprehensive consideration is given to macroeconomic factors such as administrative boundaries, built-up areas, land use, water planning, and landscape planning, as well as natural environmental characteristics such as topography, climate, and vegetation. Based on the wetland ecological restoration planning and design drawings, the project's boundaries and scope are determined; and based on the specific restoration measures and implementation schedule, the corresponding time factors for ecological restoration are determined.

[0104] S2. Establish an ecological patch identification system to identify ecological patches in forest land, grassland and water bodies within the project area before and after the implementation of ecological restoration, and determine the area size of the ecological patches.

[0105] This research developed rapid cropping and color balancing techniques for wetland remote sensing images based on multi-source monitoring data, as well as rapid production techniques for small ecological patch monitoring base maps. It integrated various ecological and environmental monitoring data. By integrating these data and combining macroscopic remote sensing and microscopic CAD data, it determined forest patches, grassland patches, and water body patches, along with their area sizes, considering factors such as the spatial distribution of ecological restoration projects and vegetation characteristics. Based on the characteristics of different regions and engineering sites, it identified and classified small ecological engineering patches in forest land, grassland, and water bodies.

[0106] For different ecological patches, the types of ecosystem services were identified and determined, as shown in Table 1.

[0107]

[0108] S3. For different ecological patch types, identify and determine the types of ecosystem services, and calculate the ecosystem services for each ecological patch type.

[0109] (1) Carbon fixation and oxygen release (BC)

[0110] ;

[0111] In the formula, It is biological carbon sequestration (sej / yr); It is the first Carbon sequestration per unit area of ​​each ecological patch type (gC / m²) 2 / yr); It is the first Average turnover time (yr) for each ecological patch type; It is the area of ​​the ecological patch type (m) 2 ); It is the energy conversion rate (sej / g) of the ecological patch type.

[0112] (2) Soil Construction (FS)

[0113] ;

[0114] In the formula It is about building soil (sej / yr), which is about building soil organic matter ( ) and minerals ( The sum of the two parts.

[0115] ① Organic matter ( )

[0116] ;

[0117] In the formula It is to build soil organic matter (sej / yr); It is the energy value (sej / yr) corresponding to the renewable resources of the ecological patch type; It is the proportion of vegetation litter to biomass (g / g, %). It is the carbon content (g / g, %) in vegetation litter.

[0118] ②Minerals ( )

[0119] ;

[0120] In the formula, It is the building block of soil minerals (sej / yr); It is the first The first ecological patch type of soil Percentage of each mineral in the total minerals (%) It is the first Soil bulk density (g / cm³) of each ecological patch type 3 ); For the first Soil thickness (cm) for each ecological patch type; For the first Ecosystem area (m²) of each ecological patch type 2 ); The percentage of soil minerals by weight (%) is taken as 95%; 10000 is m 2 to cm 2 The conversion coefficient; For the first The first ecological patch type of soil Turnover time (in years) for a type of mineral; It is the first The first ecological patch type of soil Energy conversion rate (sej / g) of a mineral.

[0121] (3) Increase bottom sediment (IS)

[0122] ;

[0123] ;

[0124] in, It is the first Increased sediment (sej / yr) in different ecological patch types (such as rivers, lakes and swamps); For the first Annual organic matter deposition (g C / m³) of each ecological patch type 2 / yr); This represents the percentage of organic matter absorbed by wetland plants, with a value of 0.78. It is the first The conversion factor (kcal / g) of g to kcal in each ecological patch type. It is the conversion factor from kcal to J, with a value of 4186 J / kcal; No. Area of ​​each ecological patch type (m²) 2 ); It is the first Energy conversion rate of deposited organic matter in each ecological patch type (sej / J); Organic sedimentary material accounts for the majority. The proportion of NPP for each ecological patch type was 30.37%. It is the first NPP levels (g C / m³) for each ecological patch type 2 / yr).

[0125] (4) Groundwater recharge (RW)

[0126] ;

[0127] in, It is the first Recharged groundwater for each ecological patch type (sej / yr); For the first Average annual precipitation (m / yr) within the range of each ecological patch type. For the first Area of ​​each ecological patch type (m²) 2 ); The density of water (kg / m³) 3 ); 1000 is the precipitation infiltration recharge coefficient; 1000 is the coefficient for converting kg to g. The energy conversion rate is the energy conversion rate of precipitation infiltrating into groundwater. Since the groundwater is replenished by precipitation, the energy conversion rate of rainwater (sej / g) is used here.

[0128] (5) Air purification (PA)

[0129] ;

[0130] In the formula, It is the first The air purification (SEJ / YR) of each ecological patch type is the amount of reduction in human health loss ( ) and reduction in ecosystem loss ( The sum of the two parts.

[0131] ① The reduction in population health loss ( )

[0132] ;

[0133] ;

[0134] In the formula, For the first The reduction in human health loss after the purification of air pollutants by each ecological patch type (sej / yr). For the first Purification of each type of ecological patch Capacity for various air pollutants (kg / ha / yr); For the first Area of ​​each ecological patch type (m²) 2 ); 0.0001 is m 2 The coefficient converted to ha; It is the first The impact factor of each air pollutant in the Eco-indicator 99 assessment framework, i.e., per kg of the first pollutant. Disability life-adjusted years (person-years / kg) caused by various air pollutants in humans. The energy value (sej / person) corresponding to the per capita total regional health expenditure; Total regional health expenditure (yuan); The regional energy value-to-currency ratio (SEJ / Yuan); This refers to the number of permanent residents (in people) in the region.

[0135] ② Reduction in ecosystem loss ( )

[0136]

[0137] ;

[0138] In the formula, Indicates the first The energy value (sej / yr) of the reduction in ecosystem loss after the purification of air pollutants by each ecological patch type; For the first Purification of each type of ecological patch Capacity for various air pollutants (kg / ha / yr); Indicates being subject to the The impact of various air pollutants on the potential extinction rate of species (PDF×m) 2×yr); 0.0001 is m 2 The conversion coefficient to ha; Indicates the first Energy required for species in each ecological patch type region (sej / yr); This represents the total number of ecological patch types.

[0139] (6) Purification of water quality (PW)

[0140] ;

[0141] In the formula, It is the first The purified water quality (sej / yr) of each ecological patch type is the amount of reduction in human health loss ( ) and reduction in ecological resource loss ( The sum of the two parts.

[0142] ① The reduction in human health loss ( )

[0143] ;

[0144] in, It is about reducing harm to the human body (sej / yr); It is the first Purification of each type of ecological patch The capacity of pollutants in water bodies (mg / kg); 0.001 is the conversion coefficient from kg to g; It is the first NPP (g C / m³) for each ecological patch type 2 / yr); It is the first Ecological patch types (m) 2 ); It is the first Disability life adjustment years (CAP) caused by water pollutants yr / kg) The coefficient for converting mg to kg is given, and its value is [value missing]. ; This represents the energy value (sej / person) corresponding to the per capita total regional health expenditure. It is the first The turnover time (yr) of water pollutants is set to 1000yr.

[0145] ② Reduction in ecological resource loss ( )

[0146] ;

[0147] in, It is about reducing ecosystem damage (sej / yr); It is the first Purification of each type of ecological patch The capacity of pollutants in water bodies (mg / kg); 0.001 is the conversion coefficient from kg to g; It is the first NPP (gC / m³) for each ecological patch type 2 / yr);

[0148] The coefficient for converting mg to kg is given, and its value is [value missing]. ; It is the first Proportion of potential species extinctions caused by Class II water pollutants (%×m) 2 ×yr×kg-1); It is the first The energy required by species in the region where each ecological patch type is located; It is the first The turnover time (yr) of pollutants in water bodies is taken as 1000yr.

[0149] (7) Soil retention (RS)

[0150] ;

[0151] ;

[0152] In the formula, It is the first Persistent soils of each ecological patch type (sej / yr); For the first Annual soil retention volume (t / yr) for each ecological patch type; 106 is the conversion coefficient from t to g; For the first Percentage of organic matter content in soil of each ecological patch type; The conversion factor for g to kcal is 5.4 (kcal / g). The coefficient for converting mg to kg is given, and its value is [value missing]. ; The conversion factor from kcal to J is 4186 J / kcal; The energy conversion rate of the soil (sej / j); For the first Potential soil erosion modulus (t / km²) for each ecological patch type 2 / yr); For the first Realistic erosion modulus (t / km²) for each ecological patch type2 / yr); For the first Area of ​​each ecological patch type (m²) 2 ).

[0153] (8) Regulate local temperature and humidity (RM)

[0154] ;

[0155] in, For the first The regulation of local temperature and humidity by ecological patch type (sej / yr); For the first Annual evapotranspiration of each ecological patch type (g / m³) 2 / yr); For the first Area of ​​each ecological patch type (m²) 2 ); The energy conversion rate of evaporation (sej / g).

[0156] (9) Climate Regulation (RC)

[0157] ;

[0158] In the formula, It is the first Each type of ecological patch regulates climate; it is a carbon sink that reduces the impact of climate change on human health. ) and the impact of carbon sinks on ecosystem quality ( The sum of the two parts.

[0159] ①Reduce carbon sequestration and the impact of climate change on human health ( )

[0160] ;

[0161] In the formula, It is the first The impact of climate change on human health due to reduced carbon sequestration in several ecological patch types (sej / yr); It is the first on a global scale The impact of each ecological patch type on greenhouse gases per unit area Average fixed amount (fixed amount (g / m) 2 / yr); 0.001 is the conversion of g to kg; It is per kg of the first Disability life-adjusted years due to greenhouse gases (person-years / kg); For the first Greenhouse gas life cycle (yr); It is the per capita total regional health expenditure (sej / cap); It is the first after VF correction The area (ha) of each ecological patch type.

[0162] ② The impact of carbon sink reduction on ecosystem quality ( )

[0163] ;

[0164] In the formula, It is the first The impact of carbon sinks of individual ecological patch types on ecosystem quality (sej / yr); It is the first on a global scale The impact of each ecological patch type on greenhouse gases per unit area Average fixed amount (fixed amount (g / m) 2 / yr); 0.001 is the conversion of g to kg; For the first Greenhouse gas life cycle (yr); It is the first The proportion of species potential extinction caused by greenhouse gases; It is the energy required by the species in the region where the i-th ecological patch type is located.

[0165] S4. Sum the ecosystem services of each ecological patch type.

[0166] Wetland ecological restoration projects are self-organizing, complex, and integrated wetland ecological spaces composed of interconnected natural elements such as water bodies, herbs, and trees. By calculating ecosystem services in woodland patches, grassland patches, and water patches within wetland ecological spaces, more accurate calculations of wetland credit can be achieved.

[0167] (1) Total ecosystem services of aquatic patches

[0168] ;

[0169] in, This represents the sum of ecosystem services provided by a water patch. Indicates biological carbon sequestration in water patches. This indicates an increase in organic matter in the bottom sediment, indicating an increase in water patch formation. This indicates the groundwater supply to the water patches. This indicates the regulation of local temperature and humidity by water patches. Indicates the purified water quality of water patches. This indicates the climate moderating effect of water patches. The unit for each ecosystem service is sej / yr.

[0170] (2) Total ecosystem services of grassland patches

[0171] ;

[0172] in, This represents the sum of ecosystem services provided by grassland patches. Indicates biological carbon sequestration in grassland patches. Indicates the soil that forms grassland patches. This indicates the groundwater supply to the grassland patches. This indicates the regulation of local temperature and humidity by grassland patches. This indicates that the grassland patches purify the atmosphere. This indicates a reduction in soil erosion in grassland patches. This indicates the climate moderating effect of grassland patches. The unit for each ecosystem service is sej / yr.

[0173] (3) Total ecosystem services of forest patches

[0174] ;

[0175] in, This represents the sum of ecosystem services provided by forest patches. Indicates biological carbon sequestration in forest patches. Indicates the soil that forms woodland patches. This indicates the groundwater replenishment of woodland patches. This indicates the role of forest patches in regulating local temperature and humidity. This indicates that woodland patches purify the atmosphere. This indicates a reduction in soil erosion in forest patches. This indicates the climate moderating effect of forest patches. The unit for each ecosystem service is sej / yr.

[0176] S5. Determine the total ecosystem services of each type of ecological patch within the scope of the wetland ecological restoration project.

[0177] = ;

[0178] In the formula, Representing the Ecosystem services of wetland ecological restoration projects Representing the Ecosystem services of water patches in wetland ecological restoration projects Representing the Ecosystem services of grassland patches in wetland ecological restoration projects Representing the Ecosystem services of forest patches in wetland ecological restoration projects.

[0179] S6. Calculate the changes in ecosystem services before and after the implementation of wetland ecological restoration projects.

[0180] ;

[0181] In the formula, Representing the Changes in ecosystem services of a wetland ecological restoration project. Representing the Ecosystem services following the implementation of wetland ecological restoration projects Representing the Ecosystem services prior to the implementation of a wetland ecological restoration project.

[0182] S7. Utilizing CAD engineering drawings and remote sensing imagery, perform multi-view, multi-level calculations on ecological patches, and combine the characteristic attributes of each ecological patch for gridding and weighted overlay calculations. Develop an initial state credit accounting method, and based on land use characteristics, develop corresponding "ecological benefit" accounting methods for different ecological patch types such as forest land, wetlands, and grasslands. Combine CAD engineering drawings for gridding calculations at scales of 1:100 and 1:10000, and remote sensing imagery for summing and statistical analysis according to different boundaries of ecological patches in ecological engineering, achieving effective measurement from micro-scale ecological engineering patches to a macro-scale 30m×30m range and correcting the measurement results at different scales; compare and integrate CAD engineering data and experimental data from micro-scale ecological restoration projects with data from large-scale major ecological projects; and classify ecological restoration projects by analyzing coastal zone ecological characteristics. Based on patch identification, perform gridded data correction, feature identification, and weighted overlay calculations based on the original data and characteristic attributes of ecological patches and ecological restoration projects, thereby determining the time risk factors and weights of wetland ecological projects. like Figure 2 As shown.

[0183] S8. Finalize the increase in wetland credit limit.

[0184] The credit limit for wetlands depends on the improvement in ecosystem services before and after the implementation of wetland ecological restoration projects, while also taking into account the time risk factor and weight of the wetland ecological restoration projects. The calculation formula is as follows:

[0185] ;

[0186] In the formula, It refers to the first The impact of wetland ecological restoration projects on wetland creditworthiness. Representing the Changes in ecosystem services of a wetland ecological restoration project. Representing the The time risk factors of wetland ecological restoration projects Indicates the first The weights of each factor in a wetland ecological restoration project.

[0187] Therefore, this invention employs the aforementioned wetland credit enhancement quantification method that integrates macroscopic remote sensing and microscopic engineering CAD data. Combining macroscopic remote sensing data and microscopic engineering CAD data, and utilizing ArcGIS, AutoCAD, and SketchUp software, it considers the spatial distribution characteristics, vegetation characteristics, remote sensing data, Normalized Difference Vegetation Index (NDVI), and Digital Elevation Model (DEM) of coastal ecological restoration projects to extract ecological patches for wetland ecological restoration projects. Based on the identification of ecological patches, and according to the combination of vegetation, water bodies, and other landscape elements, these patches are further divided into four major categories and nine types of micro-patterns: waterside woodland patches (including trees and shrubs), waterside grassland patches (including meadows, grasslands, and clumps), wetland patches (including estuaries, mudflats, and oceans), and auxiliary engineering patches (artificial non-ecological parts). Combining this with a vegetation and tree species database, the vegetation information of the ecological patches is extracted and organized to form a CAD-SU ecological engineering patch feature identification technology. The amount of wetland credit enhancement is determined based on the changes in ecosystem services before and after project implementation, as well as their temporal risk factors and weights.

[0188] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and not to limit them. Although the present invention has been described in detail with reference to preferred embodiments, those skilled in the art should understand that modifications or equivalent substitutions can still be made to the technical solutions of the present invention, and these modifications or equivalent substitutions cannot cause the modified technical solutions to deviate from the spirit and scope of the technical solutions of the present invention.

Claims

1. A wetland credit enhancement quantification method that fuses macro remote sensing and micro engineering CAD data, characterized in that, Comprising the following steps: S1, determining the project scope and measure type of wetland ecological restoration, determining the boundary and scope of the project according to the planning and design drawing of wetland ecological restoration; according to the specific restoration measures and implementation schedule of the project, the corresponding time factor of ecological restoration is determined; S2, establish an ecological patch identification system, identify the ecological patches of forest land, grassland and water body within the project scope before and after the implementation of ecological restoration and determine the area size of the ecological patches; S3, for different types of ecological patches, identify and determine the types of ecosystem services, and calculate the ecosystem services of each type of ecological patch respectively; The calculated ecosystem services are as follows: a. Carbon fixation and oxygen release BC: ； wherein, is the biological carbon sequestration, sej / yr; is the carbon sequestration of the th ecological patch type, g C / m 2 / yr; is the average turnover time of the th ecological patch type, yr; is the area of the ecological patch type, m 2 ; is the emergy conversion rate of the ecological patch type, sej / g; b. Soil construction FS: ； wherein is constructed soil, sej / yr, is constructed soil organic matter and minerals the sum of the two parts; c. Increase in sediment IS: ； ； in, It is the first Increased sediment in individual ecological patch types, sej / yr; For the first Annual organic matter deposition per ecological patch type, g C / m³ 2 / yr; This represents the percentage of organic matter absorbed by wetland plants, with a value of 0.

78. It is the first The conversion coefficient of g to kcal in each ecological patch type, kcal / g; It is the conversion factor from kcal to J, with a value of 4186 J / kcal; No. The area of ​​each ecological patch type, m 2 ; It is the first Energy conversion rate of deposited organic matter in each ecological patch type, sej / J; Organic sedimentary material accounts for the majority. The proportion of NPP for each ecological patch type was 30.37%. It is the first NPP amount for each ecological patch type, g C / m 2 / yr; d. Recharge of groundwater RW: ； wherein, is the recharge of groundwater for the th ecological patch type, sej / yr; is the average annual precipitation within the th ecological patch type, m / yr; is the area of the th ecological patch type, m 2 ; is the density of water, kg / m 3 ; is the precipitation infiltration recharge coefficient for the th ecological patch type; 1000 is the conversion factor from kg to g; is the emergy conversion rate of precipitation infiltration groundwater, here the emergy conversion rate of rainwater is used, sej / g; e. Purification of atmosphere PA: ； wherein is the amount of atmospheric purification, sej / yr, of the ecological patch type, which is the sum of the reduction in human health losses and the reduction in ecosystem losses both of which are given by the following equations: f. Purification of water quality PW: ； In the formula, is the net water quality purification of the jth ecological patch type, sej / yr, is the sum of the reduction in human health loss and the reduction in ecological resource loss of the jth ecological patch type, sej / yr, is the sum of the reduction in human health loss and the reduction in ecological resource loss g. Soil retention RS: ； ； In the formula, It is the first Persistent soil of one ecological patch type, sej / yr; For the first Annual soil retention volume for each ecological patch type, t / yr; 106 is the conversion coefficient from t to g; For the first Percentage of organic matter content in soil of each ecological patch type; is the conversion factor from g to kcal, kcal / g; is the coefficient for converting mg to kg; The conversion factor from kcal to J is 4186 J / kcal; Sej represents the energy conversion rate of the soil. For the first Potential soil erosion modulus for each ecological patch type, t / km 2 / yr; For the first Realistic erosion modulus for each ecological patch type, t / km 2 / yr; For the first The area of ​​each ecological patch type, m 2 ; h. Adjustment of local temperature and humidity RM: ； in, For the first The regulation of local temperature and humidity by ecological patch type, sej / yr; For the first Annual evapotranspiration of each ecological patch type, g / m³ 2 / yr; For the first The area of ​​each ecological patch type, m 2 ; Sej / g represents the energy conversion rate of evaporation. i. Climate regulation RC: ； In the formula, is the first regulating climate of the ecological patch type, is a carbon sink to reduce the impact of climate change on human health and carbon sink to reduce the impact of climate change on ecosystem quality The sum of the two parts; S4, add the ecosystem services of each type of ecological patch; S5, determine the total ecosystem services of each type of ecological patch in the wetland ecological restoration project scope; S6, calculate the change of ecosystem services before and after the implementation of wetland ecological restoration project; S7, use CAD engineering drawing and remote sensing image map to calculate the ecological patches from multiple perspectives and multiple levels, and combine the characteristics of each ecological patch to perform grid processing and weighted superposition operation; S8, finally determine the promotion amount of wetland credit; The promotion amount of wetland credit depends on the improvement of ecosystem services before and after the implementation of wetland ecological restoration project, and also considers the time risk factor and weight of wetland ecological restoration project, and the calculation formula is as follows: ； In the formula, It refers to the first The impact of wetland ecological restoration projects on wetland creditworthiness. Representing the Changes in ecosystem services of a wetland ecological restoration project. Representing the The time risk factors of wetland ecological restoration projects Indicates the first The weights of each factor in a wetland ecological restoration project.

2. The wetland credit enhancement quantification method that fuses macro remote sensing and micro engineering CAD data according to claim 1, characterized in that, organic matter and minerals The calculation formula is as follows: Organic matter : ； wherein is the constructed soil organic matter, sej / yr; is the first is the renewable resource corresponding to the energy of the ecotope type, sej / yr; is the proportion of vegetation litter to biomass; is the carbon content in the vegetation litter; Minerals : ； In the formula, is the soil mineral matter, sej / yr; is the soil mineral matter of the th ecological patch type; is the percentage of the th mineral matter in the total mineral matter in the soil; is the soil bulk density of the 3 th ecological patch type, g / cm is the soil thickness of the th ecological patch type, cm; is the area of the th ecological patch type, m 2 ; is the percentage of the soil mineral matter in the total soil weight, 95%; 10000 is the conversion coefficient from m 2 to cm 2 ; is the turnover time of the th mineral matter in the soil of the th ecological patch type; is the energy conversion rate of the th mineral matter in the soil of the th ecological patch type, sej / g.

3. The wetland credit enhancement quantification method that fuses macro remote sensing and micro engineering CAD data according to claim 2, characterized in that, Amount of reduction in population health loss Amount of reduction in ecosystem loss The calculation formula is as follows: Reduction in population health loss : ； ； wherein, is the amount of health loss reduction of the i-th ecological patch type, sej / yr; is the amount of health loss reduction of the i-th ecological patch type, sej / yr; is the amount of health loss reduction of the i-th ecological patch type, sej / yr; is the amount of health loss reduction of the i-th ecological patch type, sej / yr; is the amount of health loss reduction of the i-th ecological patch type, sej / yr; is the amount of health loss reduction of the i-th ecological patch type, sej / yr; is the amount of health loss reduction of the i-th ecological patch type, sej / yr; 2 ; 0.0001 is the coefficient for converting m 2 to ha; is the amount of health loss reduction of the i-th ecological patch type, sej / yr; is the amount of health loss reduction of the i-th ecological patch type, sej / yr; is the amount of health loss reduction of the i-th ecological patch type, sej / yr; is the amount of health loss reduction of the i-th ecological patch type, sej / yr; is the amount of health loss reduction of the i-th ecological patch type, sej / yr; is the amount of health loss reduction of the i-th ecological patch type, sej / yr; is the amount of health loss reduction of the i-th ecological patch type, sej / yr; Ecosystem loss reduction amount : ； In the formula, Indicates the first The energy value of the reduction in ecosystem loss after the purification of air pollutants by each ecological patch type, sej / yr; For the first Purification of each type of ecological patch Capacity to detect various air pollutants, kg / ha / yr; Indicates the subject of the first The impact of various air pollutants on the potential extinction rate of species, PDF×m 2 ×yr; 0.0001 is m 2 The conversion coefficient to ha; Indicates the first Energy required by species in the region of each ecological patch type, sej / yr; This represents the total number of ecological patch types.

4. The wetland credit enhancement quantification method that fuses macro remote sensing and micro engineering CAD data according to claim 3, characterized in that, Amount of human health loss reduction and amount of ecological resource loss reduction The calculation formula is as follows: Amount of human health loss reduction : ； wherein, is the reduction of human injury, sej / yr; is the ability of the th type of ecological patch to purify the th type of water pollutant, mg / kg; 0.001 is the conversion factor from kg to g ; is the NPP of the th type of ecological patch, g C / m 2 / yr; is the th type of ecological patch, m 2 ; is the disability-adjusted life year caused by the th type of water pollutant, cap yr / kg ; 0.001 is the conversion factor from mg to kg; is the emery of the total regional health cost per capita is the turnover time of the th type of water pollutant, with a value of 1000 yr; ecological resource loss reduction : ； wherein, is the reduction of ecosystem loss, sej / yr; is the ability of the th type of ecological patch to purify the th type of water body pollutant, mg / kg; 0.001 is the conversion factor of kg to g ; is the NPP of the th type of ecological patch, gC / m 2 / yr; is the conversion factor of mg to kg; is the proportion of potential species extinction caused by the th type of water body pollutant, % x m 2 x yr x kg -1 ; is the energy required by the species in the area where the th type of ecological patch is located; is the turnover time of the th type of water body pollutant, with a value of 1000 yr.

5. The wetland credit enhancement quantification method that fuses macro remote sensing and micro engineering CAD data according to claim 4, characterized in that, Carbon sinks reduce the effects of climate change on human health and the effects of carbon sinks on ecosystem quality The calculation formula is as follows: Carbon sinks reduce climate change on human health : ； In the formula, It is the first The impact of climate change on human health due to reduced carbon sequestration in several ecological patch types, sej / yr; It is the first on a global scale The impact of each ecological patch type on greenhouse gases per unit area The average fixed amount; 0.001 is g Convert to kg; It is per kg of the first Disability life adjustment years caused by greenhouse gases; For the first The life cycle of greenhouse gases; It is the per capita total regional health expenditure, SEJ / cap; It is the first after VF correction Area of ​​each ecological patch type; Carbon sink reduces the quality of the ecosystem : ； wherein, is the impact of the thecarbon sink ecosystem type on the quality of the ecosystem, sej / yr; is the average amount of greenhouse gas fixed per unit area of the thecarbon sink ecosystem type at a global scale; 0.001 is converted to kg; is the life cycle of the g thgreenhouse gas; is the proportion of potential species extinction caused by the thgreenhouse gas; is the proportion of potential species extinction caused by the thgreenhouse gas; is the amount of energy required by the species in the area where the thecarbon sink ecosystem type is located.

6. The wetland credit enhancement quantification method of fusing macro remote sensing and micro engineering CAD data according to claim 5, characterized in that, The content of S4 is as follows: The total ecosystem services of water body patches: ； wherein, represents the sum of ecosystem services of the water body patch, represents the carbon sequestration of the water body patch, represents the increase of organic matter in the sediment of the water body patch, represents the recharge of groundwater of the water body patch, represents the regulation of local temperature and humidity of the water body patch, represents the water quality purification of the water body patch, represents the climate regulation of the water body patch; the unit of each ecosystem service is sej / yr; The total ecosystem services of grassland patches: ； wherein, Ecosystem service total of grassland patch, Carbon fixation of grassland patch, Soil formation of grassland patch, Groundwater recharge of grassland patch, Local climate regulation of grassland patch, Air purification of grassland patch, Soil erosion reduction of grassland patch, Climate regulation of grassland patch; the unit of each ecosystem service is sej / yr; The total ecosystem services of forest land patches: ； wherein, represents the total of ecosystem services of the forest land patch, represents the carbon sequestration of the forest land patch, represents the soil formation of the forest land patch, represents the groundwater recharge of the forest land patch, represents the regulation of local temperature and humidity of the forest land patch, represents the air purification of the forest land patch, represents the reduction of water and soil loss of the forest land patch, represents the climate regulation of the forest land patch; the unit of each ecosystem service is sej / yr.

7. The wetland credit enhancement quantification method of fusing macro remote sensing and micro engineering CAD data according to claim 6, characterized in that, The content of S5 is as follows: = ； wherein represent the ecosystem services of the first wetland ecological restoration project, represent the ecosystem services of the water body patch in the first wetland ecological restoration project, represent the ecosystem services of the grassland patch in the first wetland ecological restoration project, represent the ecosystem services of the forest land patch in the first wetland ecological restoration project.

8. The wetland credit enhancement quantification method of fusing macro remote sensing and micro engineering CAD data according to claim 7, characterized in that, The content of S6 is as follows: ； In the formula, Representing the Changes in ecosystem services of a wetland ecological restoration project. Representing the Ecosystem services following the implementation of wetland ecological restoration projects Representing the Ecosystem services prior to the implementation of a wetland ecological restoration project.

Images