Ecological infrastructure prioritization method based on stormwater regulation service supply and demand assessment

Abstract

The present application relates to a method for prioritizing ecological infrastructure based on supply and demand evaluation of rainwater regulation services, comprising: dividing the study area into several urban functional area units and processing multi-source geographic and hydrological data; calculating the supply and demand of rainwater regulation services for each unit; identifying the spatial aggregation characteristics of rainwater regulation service supply and demand; dividing each urban functional area unit into different supply and demand matching types; analyzing the coupling relationship between the imbalance degree of rainwater regulation service supply and demand and the waterlogging depth under different supply and demand matching types; constructing a priority index and dividing the priority configuration level according to the priority index to develop an ecological infrastructure configuration strategy. The present application has the beneficial effect of breaking through the limitations of traditional macro-scale evaluation, constructing a rainwater regulation service supply and demand quantification framework for urban functional areas, and realizing the accurate identification of the spatial pattern of rainwater regulation service supply and demand at the meso-micro scale.

Description

Technical Field

[0001] This invention relates to the field of ecological infrastructure planning technology, and more specifically, to a method for prioritizing the allocation of ecological infrastructure based on an assessment of the supply and demand of stormwater regulation services. Background Technology

[0002] Against the backdrop of global climate change and rapid urbanization, extreme rainfall events are becoming more frequent and intense. This, combined with the continuous expansion of impervious urban surfaces, has a significant cumulative effect, leading to a substantial decline in urban stormwater storage capacity and a severe mismatch between the supply and demand of stormwater regulation services. Stormwater disasters have become a major bottleneck threatening urban safety and hindering high-quality and sustainable urban development. In recent years, my country has issued a series of policies and standards related to sponge city construction and systematic urban flood control, explicitly defining "zoning-based governance and precise policy implementation" as the core principle of stormwater management, and promoting the transformation of the governance system from macro-level, comprehensive control to meso- and micro-level refined and differentiated governance.

[0003] Currently, academic and engineering communities both domestically and internationally have conducted extensive theoretical research and engineering practice on stormwater regulation service supply and demand assessment and ecological infrastructure layout optimization, forming a series of basic assessment models and layout technologies. These provide basic technical references for urban stormwater management. However, three core technical bottlenecks still exist in practical engineering applications, resulting in a significant gap between current technologies and the needs of refined management: First, existing technologies mostly focus on macro-scales such as national and watershed levels, neglecting the micro-scale of urban functional zones. This hinders the development of stormwater regulation service supply and demand assessment methods tailored to different urban functional zones. Second, most technologies have not established a quantitative coupling relationship between the degree of imbalance between stormwater regulation service supply and demand and the actual stormwater runoff depth and the priority allocation of ecological infrastructure, leading to a disconnect between stormwater service supply and demand assessment results and ecological infrastructure engineering practices. Third, existing ecological stormwater facility configuration technologies often rely on subjective experience for uniform layout across the entire area or prioritize configuration based solely on a single runoff index. This fails to achieve differentiated and precise configuration based on the characteristics of different urban functional zones, resulting in low prevention and control benefits and unreasonable resource allocation, making it difficult to maximize the benefits of stormwater management. Summary of the Invention

[0004] The purpose of this invention is to overcome the shortcomings of the prior art and provide a method for prioritizing the allocation of ecological infrastructure based on the supply and demand assessment of stormwater regulation services.

[0005] Firstly, it provides a method for prioritizing the allocation of ecological infrastructure based on supply and demand assessments of stormwater regulation services, including:

[0006] S1. Divide the study area into several urban functional zones and collect and standardize multi-source geographic and hydrological data;

[0007] S2. Calculate the supply of stormwater regulation services for each unit based on the rainfall interception index, and calculate the demand for stormwater regulation services for each unit based on the natural demand index.

[0008] S3. Using spatial hotspot analysis, identify the spatial clustering characteristics of stormwater regulation service supply and demand;

[0009] S4. Based on the spatial agglomeration characteristics of supply and demand, each urban functional area unit is divided into different supply and demand matching types;

[0010] S5. Simulate the water depth under the target rainstorm return period and analyze the coupling relationship between the degree of supply and demand imbalance of rainstorm regulation services and water depth under different supply and demand matching types;

[0011] S6. Based on the degree of supply and demand imbalance and water depth, construct a priority index, and classify priority allocation levels according to the priority index to formulate differentiated ecological infrastructure allocation strategies.

[0012] Preferably, S1 includes:

[0013] S101. Based on the urban scale of the study area, the area range of urban functional zone units is delineated using the adaptability rule to ensure the uniformity of functional attributes within the units.

[0014] S102. Collect natural geographic data and hydrological data, and complete any missing data.

[0015] S103. Unify the integrated data to the preset resolution and coordinate system.

[0016] Preferably, S2 includes:

[0017] S201. The rainfall interception index is used as a quantitative indicator of stormwater regulation service supply. Surface runoff is calculated using the SCS-CN method, and the rainfall interception index is obtained by combining it with total rainfall. The calculation formula is as follows:

[0018]

[0019] in, The rainfall interception index for the i-th urban functional area unit; Provide a quantitative value for the stormwater regulation service of the i-th urban functional area unit; This represents the total rainfall for the i-th urban functional zone unit; Let be the surface runoff of the i-th urban functional zone unit.

[0020] S202. The Natural Demand Index is adopted as a quantitative indicator of stormwater regulation service demand. This index is constructed based on a weighted average of natural attribute indicators affecting stormwater collection and risk. The calculation formula is as follows:

[0021]

[0022] in, Let be the natural demand index of the i-th urban functional area unit; The quantitative value of the stormwater regulation service demand for the i-th urban functional area unit; Let be the impermeable surface ratio of the i-th urban functional zone unit; Let be the topographic slope of the i-th urban functional area unit; The historical rainfall and flood frequency of the i-th urban functional area unit; Let be the soil permeability coefficient of the i-th urban functional zone unit; Let be the relative elevation difference of the i-th urban functional zone unit; , , , , The weights are respectively: impermeable surface ratio, topographic slope, historical rainfall frequency, soil permeability coefficient, and relative elevation difference.

[0023] Preferably, in S3, the spatial hot and cold spot analysis method is Getis-Ord. The analysis method involves calculating the Z(t) of each unit. ) index, and based on Z( The significance test of the Z-index is used to identify hot and cold areas; The index is derived from Getis-Ord The standardized Z-value of the statistic calculated by the analytical method is the score.

[0024] As a preferred option, in S4, the demand-side Z ( The index is represented by the X-axis and the supply-side Z-axis. The index is used as the Y-axis to classify supply and demand matching types into multiple combination types to cover all urban functional area scenarios.

[0025] As a preferred embodiment, S5 includes:

[0026] S501. The hydrological and hydrodynamic model is used to simulate the water depth in the study area under the target rainstorm return period, and the model is validated.

[0027] S502. Calculate the degree of imbalance between supply and demand for stormwater regulation services and analyze its linear relationship with simulated water depth under different supply and demand matching types.

[0028] As a preferred embodiment, S6 includes:

[0029] S601. Construct an ecological infrastructure allocation priority index obtained by weighted summation of the degree of imbalance between supply and demand of stormwater regulation services and water depth, with the following formula:

[0030]

[0031] in, Let i be the priority index of the i-th urban functional area unit; The weights are respectively the degree of imbalance between supply and demand for stormwater regulation services and the depth of water accumulation; Let be the average water depth of the i-th urban functional area unit; For the demand side of the i-th urban functional area index; Supply side of the i-th urban functional area index.

[0032] S602. Based on the numerical distribution of the priority index, divide it into different priority configuration levels;

[0033] S603. For regions with different priority levels, formulate and implement ecological infrastructure configuration strategies with corresponding density and type.

[0034] Secondly, a priority allocation system for ecological infrastructure based on stormwater regulation service supply and demand assessment is provided for use in any of the methods described in the first aspect, including:

[0035] The first partitioning module is used to divide the study area into several urban functional zone units and to collect and standardize multi-source geographic and hydrological data.

[0036] The calculation module is used to calculate the supply of stormwater regulation services for each unit based on the rainfall interception index and the demand for stormwater regulation services for each unit based on the natural demand index.

[0037] The identification module is used to identify the spatial clustering characteristics of the supply and demand of stormwater regulation services by employing spatial hotspot analysis methods.

[0038] The second division module is used to divide each urban functional area unit into different supply and demand matching types based on the spatial agglomeration characteristics of supply and demand.

[0039] The simulation module is used to simulate the water depth under the target rainfall return period and analyze the coupling relationship between the degree of supply and demand imbalance of rainfall regulation services and water depth under different supply and demand matching types.

[0040] The module is used to construct a priority index based on the degree of supply and demand imbalance and water depth, and to classify priority configuration levels according to the priority index to formulate differentiated ecological infrastructure configuration strategies.

[0041] Thirdly, a computer storage medium is provided, wherein a computer program is stored therein; when the computer program is run on a computer, the computer causes the computer to perform any of the methods described in the first aspect.

[0042] Fourthly, an electronic device is provided, comprising:

[0043] Memory, used to store computer programs;

[0044] A processor for executing the computer program to implement the method as described in any of the first aspects.

[0045] The beneficial effects of this invention are:

[0046] 1. This invention provides a quantitative assessment method for the supply and demand of stormwater regulation services at the urban functional zone scale, breaking through the limitations of traditional macro-scale assessments and constructing a quantitative framework for the supply and demand of stormwater regulation services in urban functional zones. On the supply side, a slope-corrected SCS-CN method is used based on the rainfall interception index. On the demand side, a natural demand index is constructed based on the multi-index entropy weight method, combined with Getis-Ord... Spatial hotspot and colds analysis enables accurate identification of the spatial pattern of supply and demand for stormwater regulation services at the meso- and micro-scale.

[0047] 2. This invention achieves a quantitative coupling between the supply and demand imbalance of stormwater regulation services and stormwater risk. It quantitatively correlates the degree of supply and demand imbalance of stormwater regulation services with water depth. Through MIKE FLOOD model simulation and linear fitting, it establishes a response relationship of "degree of supply and demand imbalance of stormwater regulation services - water depth", so that the supply and demand assessment of stormwater regulation services is directly linked to the actual stormwater disaster risk, providing empirical evidence for allocation.

[0048] 3. This invention adopts an ecological infrastructure configuration priority index and differentiated ecological infrastructure configuration. Specifically, the priority index is constructed based on the degree of imbalance between supply and demand of stormwater regulation services and water depth. Through natural fault point classification and differentiated facility configuration, a closed loop of "assessment-identification-configuration-optimization" is formed to maximize the benefits of stormwater management and optimize the allocation of ecological rainwater facility resources. Attached Figure Description

[0049] Figure 1 A flowchart of the ecological infrastructure prioritization method based on stormwater regulation service supply and demand assessment provided for this application;

[0050] Figure 2 Schematic diagram of the spatial distribution of stormwater regulation services provided in this application;

[0051] Figure 3 A schematic diagram illustrating the spatial distribution of stormwater regulation service demand provided in this application;

[0052] Figure 4 A schematic diagram illustrating the supply and demand matching types of stormwater regulation services provided in this application;

[0053] Figure 5 Spatial distribution map of supply and demand for stormwater regulation services provided in this application;

[0054] Figure 6 The simulation diagram of water depth provided for this application;

[0055] Figure 7a A diagram illustrating the supply and demand relationship of stormwater regulation under the "high demand - medium supply" scenario provided in this application;

[0056] Figure 7b A diagram illustrating the supply and demand relationship of stormwater regulation under the "high demand - low supply" scenario provided in this application;

[0057] Figure 8 Priority zoning maps are provided for the ecological infrastructure in this application;

[0058] Figure 9 A schematic diagram of the ecological infrastructure priority configuration system provided in this application. Detailed Implementation

[0059] The present invention will be further described below with reference to embodiments. The description of the embodiments below is only for the purpose of helping to understand the present invention. It should be noted that those skilled in the art can make several modifications to the present invention without departing from the principle of the present invention, and these improvements and modifications also fall within the protection scope of the claims of the present invention.

[0060] Example 1:

[0061] To address the problems of existing technologies, Embodiment 1 of this application provides a method for prioritizing the allocation of ecological infrastructure based on the supply and demand assessment of stormwater regulation services. By quantifying the supply and demand of stormwater regulation services, deeply coupling stormwater regulation services with waterlogging, and prioritizing the allocation of ecological infrastructure, the method achieves precise layout and maximizes the benefits of stormwater management.

[0062] like Figure 1 As shown, the ecological infrastructure prioritization method based on stormwater regulation service supply and demand assessment provided in this application includes:

[0063] S1. Divide the study area into several urban functional zones and collect and standardize multi-source geographic and hydrological data.

[0064] S1 includes:

[0065] S101. Based on the urban scale of the study area, the area range of urban functional zone units is delineated using the adaptability rule to ensure the uniformity of functional attributes within the units.

[0066] Specifically, a hierarchical urban functional zone unit delineation is implemented: the study area is divided into several urban functional zone units, covering typical functional types such as residential living, comprehensive services, commercial business, industrial development, green space and leisure, and transportation hubs; the "city size adaptation rule" is adopted: the unit area of ​​large cities (urban population ≥ 5 million) is controlled at 0.5-1 km², medium-sized cities (1 million ≤ urban population < 5 million) are controlled at 1-2 km², and small cities (urban population < 1 million) are controlled at 2-3 km², ensuring the uniformity of functional attributes within the unit and solving the problem of poor adaptability of traditional fixed unit division to cities of different sizes.

[0067] For example, taking the main urban area of ​​City A as the research object, it is divided into 1,610 urban functional area units, covering six functional types: residential living, comprehensive services, commercial business, industrial development, green space and leisure, and transportation hub, with a total area of ​​880 square kilometers.

[0068] S102. Collect natural geographic data and hydrological data, and complete any missing data.

[0069] Specifically, this involves integrating and replacing missing data from multiple sources: collecting natural geographic data (land cover, DEM, building density, soil type) and hydrological data (historical rainfall points, rainfall at different return periods); and adopting a "two-dimensional missing data replacement scheme": when historical rainfall frequency data is unavailable, a composite index of "impermeable surface area × 0.6 + terrain slope × 0.4" is used as a replacement; when soil permeability coefficient data is unavailable, the default values ​​for corresponding soil types in the "Classification and Grading Standard for Soil Erosion" (SL190-2007) are referenced and combined with regional climate correction, thus breaking through the dependence of traditional methods on complete data.

[0070] For example, the natural geographic data includes 2024 land cover data at 30m resolution from the Resource and Environmental Science Data Center of the Chinese Academy of Sciences and 2024 DEM data at 30m resolution from the Geospatial Data Cloud; the hydrological data includes historical stormwater maps (2014-2024) from the Water Resources Bureau of City A and rainfall data with different return periods derived from the engineering design standard "Calculation Standard for Rainstorm Intensity" (DB33 / T 1191-2020) of the province where City A is located.

[0071] S103. Unify the integrated data to the preset resolution and coordinate system.

[0072] Specifically, data standardization is performed: the "minimum-maximum" standardization method is used to process the data, unifying it to a 10m resolution and the CGCS2000 coordinate system; a cross-coordinate system adaptation solution is added: if the data in the application area is in other coordinate systems such as WGS84, it is converted using ArcGIS "projection transformation" tools, with the error controlled within ±0.5m, improving cross-regional universality.

[0073] For example, the "minimum-maximum" normalization method was used to process all data. The 30m resolution DEM data, 1km resolution population density and GDP data were uniformly converted to 10m resolution using the "nearest neighbor sampling method" of ArcGIS 10.8 software. All data adopted the CGCS2000 coordinate system.

[0074] S2. Calculate the supply of stormwater regulation services for each unit based on the rainfall interception index, and calculate the demand for stormwater regulation services for each unit based on the natural demand index.

[0075] S2 includes:

[0076] S201. The Rainfall Retention Index (RII) is used as a supply quantification indicator. Surface runoff is calculated using the SCS-CN method, and the Rainfall Retention Index is obtained by combining it with total rainfall. The calculation formula is as follows:

[0077]

[0078] in, For the rainfall interception index of the i-th urban functional area unit, this invention will use the rainfall interception index... As a quantitative indicator of the supply of stormwater regulation services . The total rainfall (in mm) for the i-th urban functional zone unit. Let be the surface runoff of the i-th urban functional zone unit (unit: mm).

[0079] In S201, rainfall Based on the rainstorm intensity formula derived from the Zhejiang Provincial Engineering Design Standard "Calculation Standard for Rainstorm Intensity" (DB33 / T 1191-2020), the surface runoff Qi is calculated using the slope-corrected runoff curve number (CN). i ).

[0080] Soil classification and CN iThe basis for the value correction is as follows: referring to the "US Soil Conservation Service SCS Manual" and "Chinese Soil Classification and Codes" (GB / T 17296-2009), the soil type is determined by soil particle composition analysis (hydrometer method), and the soil is divided into four categories: A (sandy soil), B (sandy loam), C (loam), and D (clay). Among them, the permeability coefficient of soil in category A is > 7.6 cm / h, that of category B is 2.5-7.6 cm / h, that of category C is 0.76-2.5 cm / h, and that of category D is < 0.76 cm / h. The basis for the derivation of the CN value formula (2) after slope correction is as follows: based on the measured data of 12 cities in the eastern plain area of ​​my country (sample size n=360), the slope and CN values ​​are obtained through linear regression analysis. i The correlation model for the values ​​(R²=0.78, P<0.01) is as follows:

[0081]

[0082] In the formula, CN is the slope-corrected value for the i-th urban functional area unit; The average slope of the i-th urban functional area unit is obtained from DEM data using the ArcGIS 10.8 "Slope Calculation" tool; The base CN value for the i-th urban functional zone unit, without considering slope correction, is determined based on land cover type (e.g., green space). =60, Construction Land =90).

[0083] For example, the formula for the intensity of heavy rainfall in areas such as District A and District B of City A is:

[0084]

[0085] The formula for area C is:

[0086]

[0087] The formula for zone D is:

[0088]

[0089] Finally, the RII value of City A is calculated using formulas (1) and (2).

[0090] S202. The Natural Demand Index for Rainwater Regulation Services (NDSI) is adopted as the quantitative indicator of demand. This index is constructed based on a weighted average of natural attribute indicators affecting rainwater collection and risk, and the calculation formula is as follows:

[0091]

[0092] in, The natural demand index of the i-th urban functional area unit This invention will use the natural demand index As a quantitative indicator of the supply of stormwater regulation services. , , , , The weights for impermeable surface ratio, topographic slope, historical rainfall frequency, soil permeability coefficient, and relative elevation difference are respectively, satisfying the following conditions. + + + + The weights are determined using the entropy weight method.

[0093] Let be the impermeable surface ratio (%) of the i-th urban functional area unit; Let be the average slope (°) of the i-th urban functional zone unit; is the historical rainfall frequency (times / year) of the i-th urban functional area unit; Let be the soil permeability coefficient (cm / h) of the i-th urban functional zone unit. Let be the relative elevation difference (m) of the i-th urban functional zone unit.

[0094] For example, in city A, =0.32、 =0.18、 =0.25、 =0.15、 =0.10, calculate the NDSI value of City A using formula (6).

[0095] S3. Using spatial hotspot analysis, identify the spatial clustering characteristics of stormwater regulation service supply and demand.

[0096] Specifically, the spatial hot and cold spot analysis method is Getis-Ord. The analytical method calculates the supply-side separately. index and demand side index Further calculations , Index, combined , The significance test of the index identifies hot and cold areas of supply and demand for stormwater regulation services, and clarifies the spatial distribution characteristics of supply and demand.

[0097] Based on Getis-Ord Spatial hotspot and cold spot analysis methods to construct supply-side structural analysis. The calculation formula is as follows:

[0098]

[0099] in, This reflects the spatial clustering characteristics of stormwater regulation services. A positive difference indicates that the supply of stormwater regulation services in the i-th urban functional area and its surrounding areas is higher than the overall regional average. Therefore, This indicates that the area surrounding the unit is a spatially concentrated area where the supply of stormwater regulation services is higher than the regional average, i.e., a supply hotspot. Conversely, This indicates that the supply value of stormwater regulation services in the surrounding area is lower than the regional average, meaning that the unit and its surrounding area are spatial clusters with low supply value of stormwater regulation services, i.e., supply cold spots. To determine whether the i-th urban functional area unit and the j-th urban functional area unit are spatially adjacent, the Queen adjacency rule is used. If the i-th urban functional area unit and the j-th urban functional area unit share an edge or a vertex, then... Otherwise, non-adjacent units . The normalized stormwater regulation service supply for the j-th urban functional area unit is calculated by formula (1) and the value is normalized. The average normalized stormwater regulation service supply across all urban functional zones; The standard deviation of the normalized stormwater regulation service supply for all urban functional zones in the region; n is the total number of urban functional zones.

[0100] Similarly, on the demand side The formula for calculating the index is as follows:

[0101]

[0102] To reflect the calculation results , The extent to which the exponential rate deviates from the average level (expected value) under random conditions requires further calculation. , An index transforms a raw statistic into a deviation expressed in standard deviations. Because it converts the raw statistic into a deviation... , Transform into , The methods are consistent, and will be used uniformly below. replace , ,use replace , .

[0103] The index is calculated as follows:

[0104]

[0105] In the formula, The expected value under the assumption of random distribution, i.e. The average value of the index; for The variance of the exponent. (Through...) Significance tests are used to determine the reliability of hot or cold spot clustering.

[0106] An index ≥ 1.65 indicates a strong and significant hot spot, while 1.28 ≤ <1.65 indicates a significant hotspot, -1.65< ≤-1.28 is a significant cold spot. ≤-1.65 indicates a strongly significant cold spot.

[0107] For example, areas in City A with strong significant hot spots, significant hot spots, significant cold spots, and strong significant cold spots in the supply and demand of stormwater regulation services, such as... Figure 2 and Figure 3 As shown.

[0108] S4. Based on the spatial agglomeration characteristics of supply and demand, each urban functional area unit is divided into different supply and demand matching types.

[0109] S5. Simulate the water depth under the target rainwater return period and analyze the coupling relationship between the degree of supply and demand imbalance of rainwater regulation services and water depth under different supply and demand matching types.

[0110] S6. Based on the degree of supply and demand imbalance and water depth, construct a priority index, and classify priority allocation levels according to the priority index to formulate differentiated ecological infrastructure allocation strategies.

[0111] Example 2:

[0112] Building upon Example 1, Example 2 of this application provides a more specific method for prioritizing the allocation of ecological infrastructure based on an assessment of the supply and demand of stormwater regulation services, including:

[0113] S1. Divide the study area into several urban functional zones and collect and standardize multi-source geographic and hydrological data.

[0114] S2. Calculate the supply of stormwater regulation services for each unit based on the rainfall interception index, and calculate the demand for stormwater regulation services for each unit based on the natural demand index.

[0115] S3. Using spatial hotspot analysis, identify the spatial clustering characteristics of stormwater regulation service supply and demand.

[0116] S4. Based on the spatial agglomeration characteristics of supply and demand, each urban functional area unit is divided into different supply and demand matching types.

[0117] In S4, demand side index For the X-axis, supply side index Using the Y-axis as the basis, supply and demand matching types are divided into multiple combination types to cover all urban functional area scenarios.

[0118] For example, such as Figure 4 As shown, the supply and demand matching types are divided into 9 categories, including "high demand - high supply", "high demand - medium supply", "high demand - low supply", "medium demand - high supply", "medium demand - medium supply", "medium demand - low supply", "low demand - high supply", "low demand - medium supply", and "low demand - low supply".

[0119] Among them, the "high demand - low supply" type is mainly distributed in commercial and business areas, industrial development zones, and transportation hubs, while the "low demand - high supply" type is mainly in green spaces and leisure areas (such as...). Figure 5 (As shown).

[0120] S5. Simulate the water depth under the target rainwater return period and analyze the coupling relationship between the degree of supply and demand imbalance of rainwater regulation services and water depth under different supply and demand matching types.

[0121] S5 includes:

[0122] S501. A hydrological and hydrodynamic model was used to simulate the water depth in the study area under the target rainstorm return period, and the model was validated.

[0123] Specifically, the MIKE FLOOD model was used to simulate the waterlogging depth at flood points in the study area under the target rain-flood return period. The key parameters of the model were set as follows: the roughness coefficient of the river channel was 0.025-0.035 (0.025 for concrete channels and 0.030-0.035 for natural channels), the Manning coefficient of the surface was 0.015-0.040 (0.015 for asphalt pavement and 0.030-0.040 for green space), the grid resolution was consistent with the resolution after data preprocessing (10m), and the time step was set to 60s. The model validation criteria were: the maximum error with the measured waterlogging depth was ≤10%, and the Nash efficiency coefficient was ≥0.65 to ensure the reliability of the simulation results.

[0124] For example, the MIKE FLOOD model was used to simulate the water depth at flood-prone areas under a 50-year return period for rainstorms. Field measured data, including river network data, river cross-section data, and elevation data, were input. Parameters such as the roughness coefficient were adjusted until the maximum error between the model output value and the measured water depth on July 19, 2023, was 8.1%, and the Nash efficiency coefficient was 0.70. Using the validated MIKE FLOOD model, the water depth in the main urban area of ​​City A was simulated under rainfall conditions under a 50-year return period for rainstorms. Figure 6 ).

[0125] S502. Calculate the degree of imbalance between supply and demand for stormwater regulation services and analyze its linear relationship with simulated water depth under different supply and demand matching types.

[0126] Specifically, construct coupling relationships for all types: through formulas Calculate the degree of imbalance between supply and demand of stormwater regulation services, and analyze the linear relationship between water depth and the degree of imbalance between supply and demand of stormwater regulation services under different supply and demand matching types.

[0127] For example, calculating the degree of imbalance between supply and demand in stormwater regulation revealed that the average water depth at flood-prone areas (68.7 cm) was greater than that in areas with high demand and low supply (24.1 cm), and the two showed a linear positive correlation (e.g., Figure 7a and Figure 7b (As shown).

[0128] S6. Based on the degree of supply and demand imbalance and water depth, construct a priority index, and classify priority allocation levels according to the priority index to formulate differentiated ecological infrastructure allocation strategies.

[0129] S6 includes:

[0130] S601. Construct a priority index obtained by weighted summation of the degree of imbalance between supply and demand of stormwater regulation services and water depth. The formula is:

[0131]

[0132] in, Let i be the priority index of the i-th urban functional area unit; The weights are respectively the degree of imbalance between supply and demand for stormwater regulation services and the depth of water accumulation; Let be the average water depth of the i-th urban functional area unit; For the demand side of the i-th urban functional area index; Supply side of the i-th urban functional area index.

[0133] S602. Based on the numerical distribution of the priority index, it is divided into different priority configuration levels.

[0134] For example, based on The specific values ​​were determined using Jenkins' natural breakpoint method, dividing the priority allocation area into first-level, second-level, and third-level priority allocation areas. Figure 8 ).

[0135] S603. For regions with different priority levels, formulate and implement ecological infrastructure configuration strategies with corresponding density and type.

[0136] Specifically, for the first-priority allocation areas, high-permeability and high-capacity ecological infrastructure will be densely deployed; for the second-priority allocation areas, low-intervention methods will be used to configure ecological infrastructure, which can be further optimized in the long term; for the third-priority allocation areas, the status quo will be maintained or the ecological infrastructure will be appropriately strengthened. Among them, ecological infrastructure includes permeable pavement, rain gardens, sunken green spaces, bioretention ponds, green roofs, grassed swales, and water storage bodies.

[0137] For example, in the first-priority areas (mainly located in old streets in the city core and newly built important functional areas), facilities such as bioretention ponds, rain gardens, and sunken green spaces are densely arranged; in the second-priority areas, existing building roofs and streets are equipped with green roofs, grassed swales, and other facilities; in the third-priority areas (mainly located in high-supply areas), the status quo is maintained or ecological infrastructure is moderately strengthened.

[0138] It should be noted that the parts in this embodiment that are the same as or similar to those in Embodiment 1 can be referred to each other, and will not be repeated in this application.

[0139] Example 3:

[0140] Building upon Example 2, Example 3 of this application provides an ecological infrastructure priority allocation system based on stormwater regulation service supply and demand assessment, such as... Figure 9 As shown, it includes:

[0141] The first partitioning module is used to divide the study area into several urban functional zone units and to collect and standardize multi-source geographic and hydrological data.

[0142] The calculation module is used to calculate the supply of stormwater regulation services for each unit based on the rainfall interception index and the demand for stormwater regulation services for each unit based on the natural demand index.

[0143] The identification module is used to identify the spatial clustering characteristics of the supply and demand of stormwater regulation services by employing spatial hotspot analysis methods.

[0144] The second division module is used to divide each urban functional area unit into different supply and demand matching types based on the spatial agglomeration characteristics of supply and demand.

[0145] The simulation module is used to simulate the water depth under the target rainfall return period and analyze the coupling relationship between the degree of supply and demand imbalance of rainfall regulation services and water depth under different supply and demand matching types.

[0146] The module is used to construct a priority index based on the degree of supply and demand imbalance and water depth, and to classify priority configuration levels according to the priority index to formulate differentiated ecological infrastructure configuration strategies.

[0147] It should be noted that the system provided in this embodiment is the system corresponding to the method provided in embodiment 2. Therefore, the parts that are the same as or similar to those in embodiment 2 in this embodiment can be referred to each other, and will not be described again in this application.

Claims

1. A method for prioritizing the allocation of ecological infrastructure based on supply and demand assessment of stormwater regulation services, characterized in that, include: S1. Divide the study area into several urban functional zones and collect and standardize multi-source geographic and hydrological data; S2. Calculate the supply of stormwater regulation services for each unit based on the rainfall interception index, and calculate the demand for stormwater regulation services for each unit based on the natural demand index. S3. Using spatial hotspot analysis, identify the spatial clustering characteristics of stormwater regulation service supply and demand; S4. Based on the spatial agglomeration characteristics of supply and demand, each urban functional area unit is divided into different supply and demand matching types; S5. Simulate the water depth under the target rainstorm return period and analyze the coupling relationship between the degree of supply and demand imbalance of rainstorm regulation services and water depth under different supply and demand matching types; S6. Based on the degree of supply and demand imbalance and water depth, construct a priority index, and classify priority allocation levels according to the priority index to formulate differentiated ecological infrastructure allocation strategies.

2. The method for prioritizing the allocation of ecological infrastructure based on the supply and demand assessment of stormwater regulation services as described in claim 1, characterized in that, S1 includes: S101. Based on the urban scale of the study area, the area range of urban functional zone units is delineated using the adaptability rule to ensure the uniformity of functional attributes within the units. S102. Collect natural geographic data and hydrological data, and complete any missing data. S103. Unify the integrated data to the preset resolution and coordinate system.

3. The method for prioritizing the allocation of ecological infrastructure based on the supply and demand assessment of stormwater regulation services according to claim 2, characterized in that, S2 include: S201. The rainfall interception index is used as a quantitative indicator of stormwater regulation service supply. Surface runoff is calculated using the SCS-CN method, and the rainfall interception index is obtained by combining it with total rainfall. The calculation formula is as follows: ； in, The rainfall interception index for the i-th urban functional area unit; Provide a quantitative value for the stormwater regulation service of the i-th urban functional area unit; This represents the total rainfall for the i-th urban functional zone unit; Let be the surface runoff of the i-th urban functional zone unit; S202. The Natural Demand Index is adopted as a quantitative indicator of stormwater regulation service demand. This index is constructed based on a weighted average of natural attribute indicators affecting stormwater collection and risk. The calculation formula is as follows: ； in, Let be the natural demand index of the i-th urban functional area unit; The quantitative value of the stormwater regulation service demand for the i-th urban functional area unit; Let be the impermeable surface ratio of the i-th urban functional zone unit; Let be the topographic slope of the i-th urban functional area unit; The historical rainfall and flood frequency of the i-th urban functional area unit; Let be the soil permeability coefficient of the i-th urban functional zone unit; Let be the relative elevation difference of the i-th urban functional zone unit; , , , , The weights are respectively: impermeable surface ratio, topographic slope, historical rainfall frequency, soil permeability coefficient, and relative elevation difference.

4. The method for prioritizing the allocation of ecological infrastructure based on the supply and demand assessment of stormwater regulation services according to claim 3, characterized in that, In S3, the spatial hot and cold spot analysis method is Gittis-Odd. Spatial hot and cold spot analysis method, by calculating the hot and cold spots of each unit Index, and based on Significance tests of the index are used to identify hot and cold areas; The index is composed of Gittis-Odd The standardized Z-value score of the statistic calculated by the spatial hot and cold spot analysis method.

5. The method for prioritizing the allocation of ecological infrastructure based on the supply and demand assessment of stormwater regulation services according to claim 4, characterized in that, In S4, demand side index For the X-axis, supply side index Using the Y-axis as the basis, supply and demand matching types are divided into multiple combination types to cover all urban functional area scenarios.

6. The method for prioritizing the allocation of ecological infrastructure based on the supply and demand assessment of stormwater regulation services according to claim 5, characterized in that, S5 include: S501. The hydrological and hydrodynamic model is used to simulate the water depth in the study area under the target rainstorm return period, and the model is validated. S502. Calculate the degree of imbalance between supply and demand for stormwater regulation services and analyze its linear relationship with simulated water depth under different supply and demand matching types.

7. The method for prioritizing the allocation of ecological infrastructure based on the supply and demand assessment of stormwater regulation services according to claim 6, characterized in that, S6 include: S601. Construct a priority index obtained by weighted summation of the degree of imbalance between supply and demand of stormwater regulation services and water depth, with the following formula: ； in, Let i be the priority index of the i-th urban functional area unit; The weights are respectively the degree of imbalance between supply and demand for stormwater regulation services and the depth of water accumulation; Let be the average water depth of the i-th urban functional area unit; For the demand side of the i-th urban functional area index; Supply side of the i-th urban functional area index; S602. Based on the numerical distribution of the priority index, divide it into different priority configuration levels; S603. For regions with different priority levels, formulate and implement ecological infrastructure configuration strategies with corresponding density and type.

8. A priority allocation system for ecological infrastructure based on supply and demand assessment of stormwater regulation services, characterized in that: For performing the method according to any one of claims 1 to 7, comprising: The first partitioning module is used to divide the study area into several urban functional zone units and to collect and standardize multi-source geographic and hydrological data. The calculation module is used to calculate the supply of stormwater regulation services for each unit based on the rainfall interception index and the demand for stormwater regulation services for each unit based on the natural demand index. The identification module is used to identify the spatial clustering characteristics of the supply and demand of stormwater regulation services by employing spatial hotspot analysis methods. The second division module is used to divide each urban functional area unit into different supply and demand matching types based on the spatial agglomeration characteristics of supply and demand. The simulation module is used to simulate the water depth under the target rainfall return period and analyze the coupling relationship between the degree of supply and demand imbalance of rainfall regulation services and water depth under different supply and demand matching types. The module is used to construct a priority index based on the degree of supply and demand imbalance and water depth, and to classify priority configuration levels according to the priority index to formulate differentiated ecological infrastructure configuration strategies.

9. A computer storage medium, characterized in that, The computer storage medium stores a computer program; when the computer program is run on the computer, it causes the computer to perform the method described in any one of claims 1 to 7.

10. An electronic device, characterized in that, include: Memory, used to store computer programs; A processor for executing the computer program to implement the method as described in any one of claims 1 to 7.

Images