A three-level refined accounting system for ecological product gross value

Abstract

The application discloses an ecological product total value three-level refined accounting system, which comprises a multi-source data import module, a localization parameter library module, a core accounting engine module, a result display and export module and a batch data management module. The method comprises the following steps: obtaining spatial basic data and statistical data of a county to be accounted and preprocessing; calling a county-specific localization parameter library; based on the spatial basic data and the localization parameter library, the function amount of adjustment service and cultural service is calculated on a grid scale; based on the statistical data, the material product supply service value is accounted on a county scale; the grid-level function amount is automatically matched into village-level administrative boundary vector data through spatial superposition analysis, and is aggregated into a village-level unit based on a grid area weighting method; the material product supply service value is decomposed into each village-level unit according to a distribution weight, and the accounting results of a township level and a county level are obtained through step-by-step collection; the function amount is converted into a value amount, and the three-level refined accounting results of a village level, a township level and a county level are generated and exported into a report.

Description

Technical Field

[0001] This invention relates to the field of interdisciplinary technology of ecological environment assessment and geographic information big data, and in particular to a three-level refined accounting system for the gross ecological product value. Background Technology

[0002] Gross Ecosystem Product (GEP) accounting is an important technical means to quantify the value of ecosystem services and support the construction of ecological civilization and the asset-based management of ecological resources. Currently, accounting systems based on the UN Millennium Ecosystem Assessment (MA) framework and the Economics of Ecosystems and Biodiversity (TEEB) guidelines have been established both domestically and internationally. my country has also successively issued standards such as the "GEP Accounting Standard," providing methodological support for GEP accounting. Existing technologies typically employ tools such as the InVEST model, the CASA model, and the modified generalized soil loss equation to assess the value of three major categories of ecosystem products—material product supply, regulatory services, and cultural services—at the county or regional scale.

[0003] However, existing GEP (Gross Ecosystem Product) accounting methods generally suffer from technical deficiencies such as coarse accounting scale and insufficient spatial resolution. On the one hand, traditional accounting methods often remain at the county-level as a whole, making it difficult to finely decompose the ecological service functions to township and village-level administrative units, thus failing to meet the spatial precision requirements for ecological resource asset confirmation and precise ecological compensation policies. On the other hand, existing methods lack an automatic spatial matching mechanism between raster-level function calculation results and administrative boundary vector data, making it difficult to effectively integrate remote sensing raster data with administrative divisions. Furthermore, the parameter libraries of existing accounting systems mostly contain fixed and general parameters, lacking the ability to dynamically update localized parameters based on the natural geographical characteristics of the county. The automation level of preprocessing steps such as coordinate unification, format conversion, and integrity verification of multi-source heterogeneous data is low, resulting in low accounting efficiency and limited accuracy. Therefore, there is an urgent need for an integrated GEP accounting system that can realize raster-scale function calculation, fine spatial decomposition of three-level administrative units, dynamic updating of localized parameters, and automatic preprocessing of multi-source data. Summary of the Invention

[0004] The purpose of this invention is to provide a three-level refined accounting method for the gross ecological product value.

[0005] To achieve the above objectives, the present invention is implemented according to the following technical solution: This invention includes the following steps: S1. Obtain the spatial basic data and statistical data of the county to be calculated. Perform coordinate system conversion, format conversion, and integrity verification on the spatial basic data. Perform village-level code matching and data integrity verification on the statistical data. The spatial basic data includes village-level administrative boundary vector data, 30m resolution DEM digital elevation data, 30m resolution NDVI normalized vegetation index remote sensing data, multi-year average rainfall raster data, and soil texture raster data. The statistical data includes agricultural product output data, forestry product output data, livestock product output data, fishery product output data, scenic spot tourist visitor data, meteorological monitoring data, and socio-economic statistical data. S2. Call the county-specific localized parameter library and dynamically update the localized parameter library according to the accounting year; the localized parameter library includes soil physicochemical parameters, soil and water conservation measures factors, land use type water use coefficient, surface runoff coefficient, ecological service unit area function quantity parameters, pollutant treatment cost parameters, reservoir unit capacity construction cost parameters, and price index conversion coefficient. S3. Based on the aforementioned spatial foundation data and the localized parameter library, calculate the functional quantities of regulatory services and cultural services at the raster scale; calculate the value of material product supply services at the county scale based on the aforementioned statistical data; the regulatory services include soil conservation services, water source conservation services, carbon sequestration and oxygen release services, climate regulation services, water purification services, flood control services, windbreak and sand fixation services, air purification services, and biodiversity maintenance services; the functional quantity of soil conservation services is calculated based on the modified general soil loss equation, the functional quantity of water source conservation services is calculated based on the water balance equation, and the functional quantity of carbon sequestration and oxygen release services is calculated based on the net primary productivity model; S4. The grid-level ecosystem service function quantities are automatically matched to the village-level administrative boundary vector data through spatial overlay analysis. Based on the grid area weighting method, the function quantities are aggregated to village-level units to obtain the ecosystem service function quantities of each administrative village. The value of material product supply services is decomposed to each village-level unit according to the allocation weight. The ecosystem service function quantities of the village-level units are spatially summarized at each level to obtain the ecosystem service function quantities of townships and the total function quantities of the county. Based on the value quantity accounting parameters in the localized parameter library, the function quantities are converted into value quantities to generate three-level refined accounting results at the village, township, and county levels. S5. Summarize and statistically analyze the results of the three-level refined accounting, and generate county-level master tables, township-level detailed tables, and village-level detailed reports. Support filtering and querying by accounting year, administrative region, and ecological service type, and export them as standardized reports.

[0006] Further, in step S4, the method for automatically matching the raster-level ecosystem service functional quantities to the village-level administrative boundary vector data through spatial overlay analysis includes: Calculate the intersection area between each grid cell and the village-level administrative boundary vector data; use the proportion of the intersection area to the total area of ​​the grid cells as the spatial allocation weight, and allocate the functional quantity of the grid cells to the corresponding village-level cells; sum the functional quantities allocated to all grid cells within the same village-level cell to obtain the ecological service functional quantity of the village-level cell.

[0007] Furthermore, in step S4, the three-tiered refined allocation of the material product supply service includes: The value of agricultural products is allocated to each village-level unit according to the proportion of village-level cultivated land area to the total cultivated land area of ​​the county; the value of forestry products is allocated to each village-level unit according to the proportion of village-level forest land area to the total forest land area of ​​the county; the value of animal husbandry products is allocated to each village-level unit according to the proportion of village-level population to the total population of the county; and the value of fishery products is allocated to each village-level unit according to the proportion of village-level aquaculture pond area to the total aquaculture pond area of ​​the county.

[0008] Furthermore, in step S3, the calculation of the functional quantity of the soil conservation service based on the modified general soil loss equation includes: The slope length-slope factor LS is calculated based on the DEM digital elevation data; the vegetation cover f is calculated based on the NDVI normalized vegetation index remote sensing data, and the vegetation cover and management factor C is obtained based on the vegetation cover; the rainfall erosion factor R is calculated based on the multi-year average rainfall raster data; and the soil erodibility factor K is calculated based on the soil texture raster data. Soil retention capacity was calculated using the modified general soil loss equation. ;in Factors related to soil and water conservation measures; The total value of soil conservation is calculated by separately calculating the value of reducing non-point source pollution and the value of reducing sediment deposition based on soil retention capacity, and then summing the results. ,in To reduce non-point source pollution, To reduce the value of siltation.

[0009] Furthermore, in step S3, the calculation of the functional quantity of water conservation services based on the water balance equation includes: Rainfall amount is obtained based on the multi-year average rainfall raster data. ; Based on the DEM digital elevation data, the NDVI normalized vegetation index remote sensing data, and the potential evapotranspiration in the localized parameter library. Water use coefficient for land use type Calculate the evapotranspiration of each land use type. ; Based on the surface runoff coefficient in the localized parameter library Calculate surface runoff for each land use type ; Calculation of total water conservation source based on water balance equation ,in Let i be the area of ​​the i-th land use type; The total value of water conservation is calculated based on the reservoir unit capacity construction cost in the localized parameter database. ,in This refers to the construction cost per unit capacity of the reservoir.

[0010] Furthermore, in step S3, the functional quantity of the carbon sequestration and oxygen release service, calculated based on the net primary productivity model, includes: Net primary productivity was obtained based on the NDVI normalized vegetation index remote sensing data. ; Calculate carbon sequestration ,in The ratio of the molar mass of carbon dioxide to that of carbon; calculate the amount of oxygen released. ,in The ratio of the molar mass of oxygen to carbon dioxide; the total value of carbon sequestration and oxygen release is calculated based on the carbon dioxide and oxygen prices in the localized parameter library. ,in For the price of carbon dioxide. The price of oxygen.

[0011] Furthermore, in step S3, the calculation of the functional quantity of the climate regulation service includes: Based on the NDVI normalized vegetation index remote sensing data, the transpiration heat consumption per unit area of ​​different ecosystem types was obtained. ; Calculate the energy consumed by transpiration and evaporation in an ecosystem ,in Let i be the area of ​​the i-th ecosystem type. Evaporation from water surface; climate regulation value calculated based on electricity prices in the localized parameter database. ,in For electricity prices.

[0012] Furthermore, in step S3, the calculation of the cultural service value includes: Tourism and wellness services are accounted for using the travel expense method. ,in For the total number of tourists of type j, + , For travel time, The local average daily wage Direct costs for transportation, accommodation, entrance fees, and related expenses; The accounting results are allocated to the corresponding townships and administrative villages according to the service area of ​​the scenic area; the alternative cost method is used to account for leisure and recreation services. ,in Total time for leisure and recreation This refers to the average wage per unit time in the local area.

[0013] Secondly, a three-tiered refined accounting system for the gross ecological product value includes: The multi-source data import module is used to import spatial basic data and statistical data. The spatial basic data includes village-level administrative boundary vector data, 30m resolution DEM digital elevation data, 30m resolution NDVI normalized vegetation index remote sensing data, multi-year average rainfall raster data, and soil texture raster data. The statistical data includes agricultural product output data, forestry product output data, livestock product output data, fishery product output data, scenic spot tourist visitor data, meteorological monitoring data, and socio-economic statistical data. The multi-source data import module is also used to perform coordinate system unification, format conversion, integrity verification, and village-level code matching on the imported data. Localized parameter library module: used to store and dynamically update the county-specific localized parameter library, which includes soil physicochemical parameters, soil and water conservation measures factors, land use type water use coefficient, surface runoff coefficient, ecological service unit area function quantity parameters, pollutant treatment cost parameters, reservoir unit capacity construction cost parameters, and price index conversion coefficient; The core accounting engine module is used to calculate the functional quantities of three major categories of ecological product services—material product supply services, regulatory services, and cultural services—at the raster scale based on the spatial foundation data and the localized parameter library. The core accounting engine module has built-in GIS spatial analysis capabilities, which automatically match the raster-level ecological service functional quantities to the village-level administrative boundary vector data through spatial overlay analysis. Based on the raster area weighting method, the functional quantities are aggregated to village-level units, and the results are spatially summarized at each level to obtain township-level and county-level accounting results. Results Display and Export Module: Used to display the refined accounting results of the county, township and village levels in a tree table format. It supports filtering and querying by accounting year, administrative region and ecological service type, and can export standardized reports with one click. Batch Data Management Module: This module supports batch import of accounting data, addition of new annual accounting data, historical data tracing and querying, modification and updating of single-item service data for a single year, and background operation management of accounting tasks.

[0014] Furthermore, the core accounting engine module includes: a raster calculation submodule, used to calculate soil retention, water conservation, carbon sequestration and oxygen release, climate regulation, water purification, flood control, windbreak and sand fixation, air purification, and biodiversity maintenance based on the DEM digital elevation data, the NDVI normalized vegetation index remote sensing data, the multi-year average rainfall raster data, and the soil texture raster data; and a spatial overlay analysis submodule, used to calculate the intersection area of ​​each raster unit with the village-level administrative boundary vector data, and to calculate the intersection area as a percentage of the total raster area. The proportion of the total area of ​​the grid cells is used as the spatial allocation weight. The functional quantity of the grid cells is allocated to the corresponding village-level cells, and the functional quantities allocated to all grid cells within the same village-level cell are summed to obtain the ecological service functional quantity of the village-level cell. The three-level aggregation submodule is used to decompose the value of the material product supply service to each village-level cell according to the allocation weight. The ecological service functional quantities of the village-level cells are spatially aggregated level by level to obtain the ecological service functional quantities of each township and the total functional quantity of the county. Based on the value calculation parameters in the localized parameter library, the functional quantities are converted into value quantities.

[0015] The beneficial effects of this invention are: This invention is a three-tiered refined accounting system for the gross ecological product value. Compared with existing technologies, this invention has the following technical advantages: This invention automatically matches and aggregates ecosystem service values ​​to village, township, and county-level administrative units through grid-scale functional quantity calculation and spatial overlay analysis, overcoming the limitations of existing technologies in terms of coarse calculation scale and insufficient spatial resolution. By constructing a county-specific localized parameter library and implementing annual dynamic updates, the localization accuracy and reliability of the calculation results are improved. Through automatic preprocessing of multi-source data and batch data management, the calculation efficiency is significantly improved and historical data traceability is supported, providing high-resolution and reproducible technical support for the asset-based management of ecological resources, precise ecological compensation, and performance evaluation. Attached Figure Description

[0016] Figure 1 This is a flowchart illustrating the steps of a three-level refined accounting method for the total value of ecological products according to the present invention. Detailed Implementation

[0017] The present invention will be further described below through specific embodiments. The illustrative embodiments and descriptions herein are used to explain the present invention, but are not intended to limit the present invention.

[0018] The present invention provides a three-tiered refined accounting system for the gross ecological product value, comprising the following steps: like Figure 1 As shown, this embodiment includes the following steps: S1. Obtain the spatial basic data and statistical data of the county to be calculated. Perform coordinate system conversion, format conversion, and integrity verification on the spatial basic data. Perform village-level code matching and data integrity verification on the statistical data. The spatial basic data includes village-level administrative boundary vector data, 30m resolution DEM digital elevation data, 30m resolution NDVI normalized vegetation index remote sensing data, multi-year average rainfall raster data, and soil texture raster data. The statistical data includes agricultural product output data, forestry product output data, livestock product output data, fishery product output data, scenic spot tourist visitor data, meteorological monitoring data, and socio-economic statistical data. In the actual assessment, County A1 was taken as the research object, and the accounting year was 2025; S2. Call the county-specific localized parameter library and dynamically update the localized parameter library according to the accounting year; the localized parameter library includes soil physicochemical parameters, soil and water conservation measures factors, land use type water use coefficient, surface runoff coefficient, ecological service unit area function quantity parameters, pollutant treatment cost parameters, reservoir unit capacity construction cost parameters, and price index conversion coefficient. In the actual assessment, the following parameters were confirmed and updated in the localized parameter library module: forest water use coefficient w=2.0, shrubland w=1.5, grassland w=0.5, farmland w=0.5, and others w=0; surface runoff coefficients for paddy fields, dry land / irrigated land, orchards, forest land, grassland, construction land, and other land uses. The percentages were 34.7%, 46.96%, 8.735%, 6.119%, 10.155%, 65%, and 32.88%, respectively. The construction cost per unit reservoir capacity in 2021 was Prb = 8.64 yuan / m³. 3 (Already converted according to the Consumer Price Index); Potential Evapotranspiration Data obtained from actual measurements taken by the Wulian County Meteorological Monitoring Station in 2025; S3. Based on the aforementioned spatial foundation data and the localized parameter library, calculate the functional quantities of regulatory services and cultural services at the raster scale; calculate the value of material product supply services at the county scale based on the aforementioned statistical data; the regulatory services include soil conservation services, water source conservation services, carbon sequestration and oxygen release services, climate regulation services, water purification services, flood control services, windbreak and sand fixation services, air purification services, and biodiversity maintenance services; the functional quantity of soil conservation services is calculated based on the modified general soil loss equation, the functional quantity of water source conservation services is calculated based on the water balance equation, and the functional quantity of carbon sequestration and oxygen release services is calculated based on the net primary productivity model; S4. The grid-level ecosystem service function quantities are automatically matched to the village-level administrative boundary vector data through spatial overlay analysis. Based on the grid area weighting method, the function quantities are aggregated to village-level units to obtain the ecosystem service function quantities of each administrative village. The value of material product supply services is decomposed to each village-level unit according to the allocation weight. The ecosystem service function quantities of the village-level units are spatially summarized at each level to obtain the ecosystem service function quantities of townships and the total function quantities of the county. Based on the value quantity accounting parameters in the localized parameter library, the function quantities are converted into value quantities to generate three-level refined accounting results at the village, township, and county levels. S5. Summarize and statistically analyze the results of the three-level refined accounting, and generate county-level master tables, township-level detailed tables, and village-level detailed reports. Support filtering and querying by accounting year, administrative region, and ecological service type, and export them as standardized reports.

[0019] In this embodiment, the method for automatically matching the raster-level ecosystem service functional quantities to the village-level administrative boundary vector data through spatial overlay analysis in step S4 includes: Calculate the intersection area between each grid cell and the village-level administrative boundary vector data; use the proportion of the intersection area to the total area of ​​the grid cells as the spatial allocation weight, and allocate the functional quantity of the grid cells to the corresponding village-level cells; sum the functional quantities allocated to all grid cells within the same village-level cell to obtain the ecological service functional quantity of the village-level cell.

[0020] In this embodiment, step S4, the three-tiered refined allocation of the material product supply service, includes: The value of agricultural products is allocated to each village-level unit according to the proportion of village-level cultivated land area to the total cultivated land area of ​​the county; the value of forestry products is allocated to each village-level unit according to the proportion of village-level forest land area to the total forest land area of ​​the county; the value of animal husbandry products is allocated to each village-level unit according to the proportion of village-level population to the total population of the county; and the value of fishery products is allocated to each village-level unit according to the proportion of village-level aquaculture pond area to the total aquaculture pond area of ​​the county.

[0021] In this embodiment, step S3, where the functional quantity of the soil conservation service is calculated based on the modified general soil loss equation, includes: The slope length-slope factor LS is calculated based on the DEM digital elevation data; vegetation cover is calculated based on the NDVI normalized vegetation index remote sensing data. The vegetation cover and management factor C is obtained based on the vegetation cover; the rainfall erosion factor R is calculated based on the multi-year average rainfall raster data; and the soil erodibility factor is calculated based on the soil texture raster data. ; Soil retention capacity was calculated using the modified general soil loss equation. ;in Factors related to soil and water conservation measures; The total value of soil conservation is calculated by separately calculating the value of reducing non-point source pollution and the value of reducing sediment deposition based on soil retention capacity, and then summing the results. ,in To reduce non-point source pollution, To reduce the value of siltation; In actual assessment, , ,in For soil bulk density (t / m3), take 1.44 t / m3. 3 , The values ​​represent the pure nitrogen and phosphorus content in the soil, with nitrogen content at 0.203% and phosphorus content at 0.447%. The environmental engineering degradation cost for nitrogen is 1750 yuan / t, and for phosphorus it is 2800 yuan / t. The sediment deposition coefficient is 0.24, and the reservoir dredging cost is 18.8 yuan / m³. 3 ; KEPIC=(0.2+0.3exp(0.0256 SAN(1-SIL / 100))) (SIL / (CLA+SIL))0.3 (1-(0.25*SOC) / (SOC+exp(3.72-2.95*SOC))) (1-(0.7 (1-SAN / 100)) / ((1-SAN / 100)+exp(22.9*(1-SAN / 100)-5.51))), where SAN is the percentage of sand, SIL is the percentage of silt, CLA is the percentage of clay, and SOC is the percentage of organic carbon, all extracted from soil texture raster data.

[0022] In this embodiment, the calculation of the functional quantity of water conservation services based on the water balance equation in step S3 includes: Rainfall amount is obtained based on the multi-year average rainfall raster data. ; Based on the DEM digital elevation data, the NDVI normalized vegetation index remote sensing data, and the potential evapotranspiration in the localized parameter library. Water use coefficient for land use type Calculate the evapotranspiration of each land use type. ; Based on the surface runoff coefficient in the localized parameter library Calculate surface runoff for each land use type ; Calculation of total water conservation source based on water balance equation ,in Let i be the area of ​​the i-th land use type. ; The total value of water conservation is calculated based on the reservoir unit capacity construction cost in the localized parameter database. ,in Construction cost per unit capacity of the reservoir; In the actual assessment, based on the "Evaluation Standard for Forest Ecosystem Service Functions" (LYT 1721-2008), the construction cost per unit capacity of the reservoir in 2005 was 6.11 yuan / t. Converted to 2021 using the Consumer Price Index, this amount is 8.64 yuan / m³. 3 The water use coefficient for forests is 2, for shrublands it is 1.5, for grasslands it is 0.5, for farmland it is 0.5, and for artificial surfaces, land and deserts it is 0. The average surface runoff coefficients for paddy fields, dry land / irrigated land, orchards, forest land, grassland, construction land and other land uses are 34.7, 46.96, 8.735, 6.119, 10.155, 65 and 32.88, respectively.

[0023] In this embodiment, step S3, where the functional quantity of the carbon sequestration and oxygen release service is calculated based on the net primary productivity model, includes: Net primary productivity was obtained based on the NDVI normalized vegetation index remote sensing data. ; Calculate carbon sequestration ,in The ratio of the molar mass of carbon dioxide to that of carbon; calculate the amount of oxygen released. ,in The ratio of the molar mass of oxygen to carbon dioxide; the total value of carbon sequestration and oxygen release is calculated based on the carbon dioxide and oxygen prices in the localized parameter library. ,in For the price of carbon dioxide. For the price of oxygen; In the actual assessment, the molar mass ratio of carbon dioxide to carbon was 44 / 12; the molar mass ratio of oxygen to carbon dioxide was 32 / 44. Using the latest carbon emission social cost published by the United States, the price of carbon dioxide was US$51 / ton. After conversion based on exchange rate and price, the price of carbon dioxide in 2021 was RMB357.714 / ton. Referring to the oxygen price recommended in the "Evaluation Specification for Forest Ecosystem Service Functions" (LYT 1721-2008), the price of oxygen in 2021 was determined to be RMB1346.76 / ton through price conversion.

[0024] In this embodiment, the calculation of the functional quantity of the climate regulation service in step S3 includes: Based on the NDVI normalized vegetation index remote sensing data, the transpiration heat consumption per unit area of ​​different ecosystem types was obtained. ; Calculate the energy consumed by transpiration and evaporation in an ecosystem ,in Let i be the area of ​​the i-th ecosystem type. Evaporation from water surface; climate regulation value calculated based on electricity prices in the localized parameter database. ,in For electricity price; In the actual assessment, referring to existing research findings, the GPP of the forest, shrubland, grassland, and farmland ecosystems in AI County was determined to be 7.04. 107 KJ / km²·d, 3.92 107 KJ / km²·d, 2.56 107 KJ / km²·d, 3.24 Based on the actual conditions of AI County, the evaporation intensity of natural water bodies in AI County in 2021 was determined to be 1649.0 mm, and the evaporation intensity of inland tidal flats was taken as 73.85 mm; γ is the power consumption (kW·h) to convert 1 m³ of water into steam, taken as 120 kW·h / m³. 3 .

[0025] In this embodiment, step S3, the calculation of the cultural service value, includes: Tourism and wellness services are accounted for using the travel expense method. ,in For the total number of tourists of type j, + , For travel time, The local average daily wage Direct costs for transportation, accommodation, entrance fees, and related expenses; The accounting results are allocated to the corresponding townships and administrative villages according to the service area of ​​the scenic area; the alternative cost method is used to account for leisure and recreation services. ,in Total time for leisure and recreation This refers to the average wage per unit of time in the local area. In actual assessment, Transportation costs for tourists from location j to the natural scenic area Food and lodging expenses Ticket fees Tourism-driven spending on shopping, entertainment, and other related activities. .

[0026] Secondly, a three-tiered refined accounting system for the gross ecological product value includes: The multi-source data import module is used to import spatial basic data and statistical data. The spatial basic data includes village-level administrative boundary vector data, 30m resolution DEM digital elevation data, 30m resolution NDVI normalized vegetation index remote sensing data, multi-year average rainfall raster data, and soil texture raster data. The statistical data includes agricultural product output data, forestry product output data, livestock product output data, fishery product output data, scenic spot tourist visitor data, meteorological monitoring data, and socio-economic statistical data. The multi-source data import module is also used to perform coordinate system unification, format conversion, integrity verification, and village-level code matching on the imported data. Localized parameter library module: used to store and dynamically update the county-specific localized parameter library, which includes soil physicochemical parameters, soil and water conservation measures factors, land use type water use coefficient, surface runoff coefficient, ecological service unit area function quantity parameters, pollutant treatment cost parameters, reservoir unit capacity construction cost parameters, and price index conversion coefficient; The core accounting engine module is used to calculate the functional quantities of three major categories of ecological product services—material product supply services, regulatory services, and cultural services—at the raster scale based on the spatial foundation data and the localized parameter library. The core accounting engine module has built-in GIS spatial analysis capabilities, which automatically match the raster-level ecological service functional quantities to the village-level administrative boundary vector data through spatial overlay analysis. Based on the raster area weighting method, the functional quantities are aggregated to village-level units, and the results are spatially summarized at each level to obtain township-level and county-level accounting results. Results Display and Export Module: Used to display the refined accounting results of the county, township and village levels in a tree table format. It supports filtering and querying by accounting year, administrative region and ecological service type, and can export standardized reports with one click. Batch Data Management Module: This module supports batch import of accounting data, addition of new annual accounting data, historical data tracing and querying, modification and updating of single-item service data for a single year, and background operation management of accounting tasks.

[0027] In this embodiment, the core calculation engine module includes: a raster calculation submodule, used to calculate soil retention, water conservation, carbon sequestration and oxygen release, climate regulation, water purification, flood control, windbreak and sand fixation, air purification, and biodiversity maintenance based on the DEM digital elevation data, the NDVI normalized vegetation index remote sensing data, the multi-year average rainfall raster data, and the soil texture raster data; and a spatial overlay analysis submodule, used to calculate the intersection area of ​​each raster unit with the village-level administrative boundary vector data, and to calculate the intersection area as a percentage of the total area of ​​the village-level administrative boundary vector data. The proportion of the total area of ​​the grid cells is used as the spatial allocation weight. The functional quantity of the grid cells is allocated to the corresponding village-level units, and the functional quantities allocated to all grid cells within the same village-level unit are summed to obtain the ecological service functional quantity of the village-level unit. The three-level aggregation submodule is used to decompose the value of the material product supply service to each village-level unit according to the allocation weight. The ecological service functional quantities of the village-level units are spatially aggregated level by level to obtain the ecological service functional quantities of each township and the total functional quantity of the county. Based on the value calculation parameters in the localized parameter library, the functional quantities are converted into value quantities.

[0028] In the actual assessment, the calculation results were spatially matched to 651 administrative villages, and the total value at the township and county levels was obtained by summarizing them level by level. After the calculation was completed, the soil conservation module was used to select the year 2025 to view the soil conservation value of 12 townships (subdistricts) including Hongning Subdistrict, Gaoze Subdistrict, and Jietou Town. The system automatically generated an Excel report containing details of 651 administrative villages and downloaded it to the local machine to complete this calculation. This calculation yielded the total soil conservation value of AI County in 2025, with Hongning Subdistrict at 47,903,456 yuan, Gaoze Subdistrict at 20,726,852 yuan, and Jietou Town at 45,262,368 yuan. The deviation from the manual calculation result was less than 1%, verifying the accuracy of the calculation results. The total water conservation value of A1 County in 2025, after system calculation, deviated from the manual calculation result by less than 1.5%. The total carbon sequestration and oxygen release value deviated by less than 1.2%, verifying the accuracy and stability of the system in multi-service joint calculation.

[0029] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

Claims

1. A three-tiered refined accounting method for the gross ecological product value, characterized in that, Includes the following steps: S1. Obtain the spatial basic data and statistical data of the county to be calculated. Perform coordinate system conversion, format conversion, and integrity verification on the spatial basic data. Perform village-level code matching and data integrity verification on the statistical data. The spatial basic data includes village-level administrative boundary vector data, 30m resolution DEM digital elevation data, 30m resolution NDVI normalized vegetation index remote sensing data, multi-year average rainfall raster data, and soil texture raster data. The statistical data includes agricultural product output data, forestry product output data, livestock product output data, fishery product output data, scenic spot tourist visitor data, meteorological monitoring data, and socio-economic statistical data. S2. Call the county-specific localized parameter library and dynamically update the localized parameter library according to the accounting year; the localized parameter library includes soil physicochemical parameters, soil and water conservation measures factors, land use type water use coefficient, surface runoff coefficient, ecological service unit area function quantity parameters, pollutant treatment cost parameters, reservoir unit capacity construction cost parameters, and price index conversion coefficient. S3. Based on the aforementioned spatial foundation data and the localized parameter library, calculate the functional quantities of regulatory services and cultural services at the raster scale; calculate the value of material product supply services at the county scale based on the aforementioned statistical data; the regulatory services include soil conservation services, water source conservation services, carbon sequestration and oxygen release services, climate regulation services, water purification services, flood control services, windbreak and sand fixation services, air purification services, and biodiversity maintenance services; the functional quantity of soil conservation services is calculated based on the modified general soil loss equation, the functional quantity of water source conservation services is calculated based on the water balance equation, and the functional quantity of carbon sequestration and oxygen release services is calculated based on the net primary productivity model; S4. The grid-level ecosystem service function quantities are automatically matched to the village-level administrative boundary vector data through spatial overlay analysis. Based on the grid area weighting method, the function quantities are aggregated to village-level units to obtain the ecosystem service function quantities of each administrative village. The value of material product supply services is decomposed to each village-level unit according to the allocation weight. The ecosystem service function quantities of the village-level units are spatially summarized at each level to obtain the ecosystem service function quantities of townships and the total function quantities of the county. Based on the value quantity accounting parameters in the localized parameter library, the function quantities are converted into value quantities to generate three-level refined accounting results at the village, township, and county levels. S5. Summarize and statistically analyze the results of the three-level refined accounting, and generate county-level master tables, township-level detailed tables, and village-level detailed reports. Support filtering and querying by accounting year, administrative region, and ecological service type, and export them as standardized reports.

2. The method for three-level refined accounting of the total value of ecological products according to claim 1, characterized in that, The method for automatically matching the raster-level ecosystem service functions to the village-level administrative boundary vector data through spatial overlay analysis in step S4 includes: Calculate the intersection area between each grid cell and the village-level administrative boundary vector data; use the proportion of the intersection area to the total area of ​​the grid cells as the spatial allocation weight, and allocate the functional quantity of the grid cells to the corresponding village-level cells; sum the functional quantities allocated to all grid cells within the same village-level cell to obtain the ecological service functional quantity of the village-level cell.

3. The method for three-level refined accounting of the total value of ecological products according to claim 1, characterized in that, In step S4, the three-tiered refined allocation of material product supply services includes: The value of agricultural products is allocated to each village-level unit according to the proportion of village-level cultivated land area to the total cultivated land area of ​​the county; the value of forestry products is allocated to each village-level unit according to the proportion of village-level forest land area to the total forest land area of ​​the county; the value of animal husbandry products is allocated to each village-level unit according to the proportion of village-level population to the total population of the county; and the value of fishery products is allocated to each village-level unit according to the proportion of village-level aquaculture pond area to the total aquaculture pond area of ​​the county.

4. The method for three-level refined accounting of the total value of ecological products according to claim 1, characterized in that, In step S3, the functional quantity of the soil conservation service is calculated based on the modified general soil loss equation, including: The slope length-slope factor LS is calculated based on the DEM digital elevation data; the vegetation cover f is calculated based on the NDVI normalized vegetation index remote sensing data, and the vegetation cover and management factor C is obtained based on the vegetation cover; the rainfall erosion factor R is calculated based on the multi-year average rainfall raster data; and the soil erodibility factor K is calculated based on the soil texture raster data. Soil retention capacity was calculated using the modified general soil loss equation. ;in Factors related to soil and water conservation measures; The total value of soil conservation is calculated by separately calculating the value of reducing non-point source pollution and the value of reducing sediment deposition based on soil retention capacity, and then summing the results. ,in To reduce non-point source pollution, To reduce the value of siltation.

5. The method for three-level refined accounting of the total value of ecological products according to claim 1, characterized in that, In step S3, the functional quantities of water conservation services are calculated based on the water balance equation, including: Rainfall amount is obtained based on the multi-year average rainfall raster data. ; Based on the DEM digital elevation data, the NDVI normalized vegetation index remote sensing data, and the potential evapotranspiration in the localized parameter library. Water use coefficient for land use type Calculate the evapotranspiration of each land use type. ; Based on the surface runoff coefficient in the localized parameter library Calculate surface runoff for each land use type ; Calculation of total water conservation source based on water balance equation ,in Let i be the area of ​​the i-th land use type; The total value of water conservation is calculated based on the reservoir unit capacity construction cost in the localized parameter database. ,in This refers to the construction cost per unit capacity of the reservoir.

6. The method for three-level refined accounting of the total value of ecological products according to claim 1, characterized in that, In step S3, the functional quantity of the carbon sequestration and oxygen release service is calculated based on the net primary productivity model, including: Net primary productivity was obtained based on the NDVI normalized vegetation index remote sensing data. ; Calculate carbon sequestration ,in The ratio of the molar mass of carbon dioxide to that of carbon; calculate the amount of oxygen released. ,in The ratio of the molar mass of oxygen to carbon dioxide; the total value of carbon sequestration and oxygen release is calculated based on the carbon dioxide and oxygen prices in the localized parameter library. ,in For the price of carbon dioxide. The price of oxygen.

7. The method for three-level refined accounting of the total value of ecological products according to claim 1, characterized in that, In step S3, the calculation of the functional quantity of the climate regulation service includes: Based on the NDVI normalized vegetation index remote sensing data, the transpiration heat consumption per unit area of ​​different ecosystem types was obtained. ; Calculate the energy consumed by transpiration and evaporation in an ecosystem ,in Let i be the area of ​​the i-th ecosystem type. Evaporation from water surface; climate regulation value calculated based on electricity prices in the localized parameter database. ,in For electricity prices.

8. The method for three-level refined accounting of the total value of ecological products according to claim 1, characterized in that, In step S3, the calculation of the value of cultural services includes: Tourism and wellness services are accounted for using the travel expense method. ,in For the total number of tourists of type j, + , For travel time, The local average daily wage Direct costs for transportation, accommodation, entrance fees, and related expenses; The accounting results are allocated to the corresponding townships and administrative villages according to the service area of ​​the scenic area; the alternative cost method is used to account for leisure and recreation services. ,in Total time for leisure and recreation This refers to the average wage per unit time in the local area.

9. A three-tiered refined accounting system for the gross ecological product value, used to execute the method described in any one of claims 1-8, characterized in that, include: The multi-source data import module is used to import spatial basic data and statistical data. The spatial basic data includes village-level administrative boundary vector data, 30m resolution DEM digital elevation data, 30m resolution NDVI normalized vegetation index remote sensing data, multi-year average rainfall raster data, and soil texture raster data. The statistical data includes agricultural product output data, forestry product output data, livestock product output data, fishery product output data, scenic spot tourist visitor data, meteorological monitoring data, and socio-economic statistical data. The multi-source data import module is also used to perform coordinate system unification, format conversion, integrity verification, and village-level code matching on the imported data. Localized parameter library module: used to store and dynamically update the county-specific localized parameter library, which includes soil physicochemical parameters, soil and water conservation measures factors, land use type water use coefficient, surface runoff coefficient, ecological service unit area function quantity parameters, pollutant treatment cost parameters, reservoir unit capacity construction cost parameters, and price index conversion coefficient; The core accounting engine module is used to calculate the functional quantities of three major categories of ecological product services—material product supply services, regulatory services, and cultural services—at the raster scale based on the spatial foundation data and the localized parameter library. The core accounting engine module has built-in GIS spatial analysis capabilities, which automatically match the raster-level ecological service functional quantities to the village-level administrative boundary vector data through spatial overlay analysis. Based on the raster area weighting method, the functional quantities are aggregated to village-level units, and the results are spatially summarized at each level to obtain township-level and county-level accounting results. Results Display and Export Module: Used to display the refined accounting results of the county, township and village levels in a tree table format. It supports filtering and querying by accounting year, administrative region and ecological service type, and can export standardized reports with one click. Batch Data Management Module: This module supports batch import of accounting data, addition of new annual accounting data, historical data tracing and querying, modification and updating of single-item service data for a single year, and background operation management of accounting tasks.

Images