A method for evaluating urban waterlogging resilience based on multi-dimensional indexes and coupling models

Abstract

The application discloses a kind of urban waterlogging resilience evaluation methods based on multidimensional index and coupling model, belong to urban waterlogging prevention and resilience assessment technical field.The method includes: building the resilience evaluation index system covering three dimensions of city space, lifeline engineering and rapid recovery;Establish a two-dimensional coupled hydraulic model, simulate different rainfall scenarios, output dynamic waterlogging parameters;Using comprehensive weighting method and TOPSIS model to quantify the resilience score of each spatial unit;Identify low resilience aggregation area and locate key impact indicators;According to the identification results, targeted improvement measures are developed, and the effect of the measures is verified by the coupling model closed loop.The application solves the problems of one-sided index system, lack of dynamic simulation and lack of governance closed loop in existing evaluation methods, and realizes high-precision, dynamic evaluation and scientific optimization of urban waterlogging resilience.

Description

Technical Field

[0001] This invention belongs to the field of urban waterlogging prevention and resilience assessment technology, specifically involving an urban waterlogging resilience evaluation method based on multi-dimensional indicators and coupled models. Background Technology

[0002] The intensification of global climate change and the acceleration of urbanization have led to frequent extreme rainfall events, making urban flooding a key issue restricting sustainable urban development. Since the concept of "resilient cities" was proposed, urban flooding resilience assessment has become a research hotspot in disaster prevention and mitigation. However, existing technologies have significant limitations: First, the indicator system is one-sided, focusing primarily on physical facilities (such as pipe networks and pumping stations) while neglecting socio-economic and emergency response dimensions, making it difficult to cover the entire process before, during, and after a disaster. Second, dynamic simulation is lacking; traditional methods rely on static statistical data and fail to incorporate hydrological and hydrodynamic models to obtain dynamic flooding parameters, resulting in a disconnect between assessment results and actual flooding processes. Third, weight calculations are simplistic, employing only the analytic hierarchy process (subjective) or entropy weighting (objective), leading to significant biases. Fourth, there is a lack of closed-loop improvement; the assessment process is often a one-way "data input - result output" model, making it difficult to guide actual governance work. To address these issues, this invention proposes an evaluation method that integrates multi-dimensional indicators, coupled models, and combined weights, filling a gap in existing technologies. Summary of the Invention

[0003] In view of this, the purpose of this invention is to provide an urban flooding resilience evaluation method based on multi-dimensional indicators and coupled models. By constructing a full-dimensional indicator system, accurate dynamic simulation, scientific quantification, and closed-loop improvement, a high-precision assessment and optimization of urban flooding resilience can be achieved.

[0004] To achieve the above objectives, the present invention provides the following technical solution:

[0005] A method for evaluating urban flood resilience based on multi-dimensional indicators and coupled models includes the following steps:

[0006] S1. Select indicators from three dimensions—urban space, lifeline engineering, and rapid recovery—to construct a three-level urban flood resilience evaluation index system;

[0007] S2. Construct a one-dimensional coupled urban flood simulation model using IFMS / Urban. Collect basic data of the study area and preprocess it. Perform the following operations using IFMS / Urban software: import DEM data to generate a two-dimensional surface grid, and perform elevation uplift processing on the grid in the building area; digitize the river network and drainage network, generalize river cross-sections, pipeline parameters, and manhole attributes, and construct a one-dimensional model; establish a one-dimensional coupling relationship using orifice connection method and set core modeling parameters; input measured rainfall data to simulate rainfall scenarios with different return periods, and extract key flooding parameters from the coupled simulation results for quantification of evaluation indicators.

[0008] S3. Quantitative assessment of urban waterlogging resilience level: Based on key waterlogging parameters and a three-level urban waterlogging resilience evaluation index system, the weight of each index is determined by a comprehensive weighting method combining the analytic hierarchy process and the entropy weight method. The resilience score of each assessment unit is calculated by the superior-inferior solution distance method, and the resilience level is divided into five levels: first resilience, second resilience, third resilience, fourth resilience and fifth resilience.

[0009] It also includes S4 and S5;

[0010] S4. Identify the fifth resilience cluster region through Moran's I index and LISA cluster analysis, and screen key influencing indicators by combining indicator weights and the mean values ​​of indicators for different resilience levels.

[0011] S5. Based on the identified fifth resilience cluster area and key influencing indicators, formulate resilience enhancement measures and cyclically execute S2~S4 to verify the effectiveness of the measures, forming a closed-loop resilience enhancement process.

[0012] The three-level urban flood resilience evaluation index system described in S1 includes a target layer, a criterion layer, and an indicator layer. The target layer is urban flood resilience, the criterion layer includes three dimensions: urban space, lifeline engineering, and rapid recovery, and the indicator layer includes 20 indicators.

[0013] The urban spatial dimension includes seven indicators: green coverage rate, water coverage rate, topographic slope, ground elevation, distance from river network, distance from flood-prone areas, and ground water depth. Among them, green coverage rate, water coverage rate, topographic slope, ground elevation, distance from river network, and distance from flood-prone areas are positive indicators, while ground water depth is a negative indicator.

[0014] The lifeline project dimension includes eight indicators: distance to medical institutions, road network density, distance to transportation hubs, drainage network density, distance to telecommunications maintenance stations, water supply network density, distance to power maintenance stations, and population density. Among these, road network density, drainage network density, and water supply network density are positive indicators, while distance to medical institutions, distance to transportation hubs, telecommunications maintenance stations, power maintenance stations, and population density are negative indicators.

[0015] The rapid recovery dimension includes five indicators: nighttime lighting, distance to drainage pumping stations, distance to sluice gates, distance to emergency rescue stations, and distance to disaster shelters; among them, nighttime lighting, distance to drainage pumping stations, and distance to sluice gates are positive indicators, while distance to emergency rescue stations and distance to disaster shelters are negative indicators.

[0016] S4 includes the following steps:

[0017] S41. Fifth resilience region identification: Load resilience proximity in ArcGIS software. Given raster data, calculate Moran's I index. and This indicates a significant spatial clustering of toughness levels; cluster maps are generated using LISA cluster analysis to identify the concentrated distribution area of ​​the fifth toughness unit;

[0018] S42. Key indicator identification and calculation of weighted standard deviation for each indicator. Compared to the first toughness, With the fifth resilience, that is Regional indicator mean, filtering Key influencing indicators with large and significant differences in means; the criteria for determining significant differences in means are as follows: ;

[0019] The formula for calculating the weighted standard deviation is as follows:

[0020]

[0021] In the formula The standard deviation of the standardized values ​​of the indicator. The mean, This is the overall weight.

[0022] As a further preferred embodiment of the present invention, S2 includes the following steps:

[0023] S21. Collect basic data for the study area: Use a combination of field surveys and data collection to obtain basic data, including basic geographic and infrastructure data required for model building, boundary condition data required for model operation, and measured data required for model accuracy verification.

[0024] S22. Basic data preprocessing and one-dimensional model construction: The collected basic data is standardized and preprocessed using geographic information system software, including digital elevation model filling and smoothing, slope analysis, land use classification system regularization, and verification of river network and pipeline network topology. A one-dimensional river network model is constructed using the river channel generalization method in IFMS / Urban software. A one-dimensional drainage pipeline network model system is formed by structural generalization of the pipeline system, consisting of nodes, pipe segments, and discharge outlets.

[0025] S23. Sub-catchment delineation and 2D model construction: Primary drainage zones are delineated according to the distribution of river network and pipe network, and then secondary sub-catchments are subdivided based on the Thiessen polygon method to extract the core topography and underlying surface parameters of the sub-catchments; a 2D surface model is constructed using unstructured quadrilateral meshing technology, and the mesh scale is determined according to the principle of denser meshing around the river channel and conventional meshing in other areas. Main roads are used as meshing constraints, and the mesh elevation is assigned in combination with topographic data. The effect of buildings on water flow obstruction is simulated by adjusting the elevation.

[0026] S24. Coupling of one-dimensional and two-dimensional models: In the coupling module of IFMS / Urban software, inspection wells are used as key connection nodes of one-dimensional and two-dimensional models. Through orifice connection technology, a two-way hydraulic coupling relationship between the one-dimensional river network-pipeline network model and the two-dimensional surface model is established to realize the dynamic water volume interaction simulation among the river, pipeline and surface.

[0027] S25. Model Validation: Selecting actual rainfall events as input conditions, a combination of matching degree verification of flood-prone areas and quantitative error analysis of flood depth is adopted. By comparing simulation results with actual data, the model's ability to replicate flood processes and its accuracy and reliability are verified.

[0028] S26. Model parameter settings and calculation scheme determination: Based on the characteristics of the underlying surface and the hydraulic motion law of the study area, the core hydraulic parameters of the permeable zone, impermeable zone, pipe network and river channel are set; the design rainfall intensity of different return periods is calculated using the rainfall intensity formula published for the study area; the rainfall time history is allocated by adopting the Chicago rainfall pattern through the spatiotemporal distribution characteristics of regional rainfall; and the model simulation calculation scheme is determined.

[0029] S27. Extraction of key parameters for urban flooding: Extract key parameters for urban flooding in the study area at the current level of the year from the coupled simulation results, including the surface water depth and the distribution of flood-prone points, for the quantification of evaluation indicators.

[0030] As a further preferred embodiment of the present invention, S3 includes the following steps:

[0031] S31. Standardization of indicator data: The original data of 20 indicators are organized into an n×m evaluation matrix in 100m×100m evaluation units, where n is the number of evaluation units and m=20. The original data includes the key waterlogging parameters extracted in S2 and the distance and density parameters obtained from GIS spatial analysis. The extreme value method is used to eliminate the influence of dimensions. The standardized values ​​of both positive and negative indicators are calculated according to the formula.

[0032] For n evaluation objects and m evaluation indicators, establish an evaluation matrix. The evaluation matrix is ​​shown below:

[0033]

[0034] In the formula, i = 1, 2, ..., n; j = 1, 2, ..., m;

[0035] The formula for calculating the positive indicator is as follows:

[0036]

[0037] In the formula, Let j be the standardized value of the j-th indicator in the i-th evaluation unit, ranging from 0 to 1. The original value of the indicator. , These are the maximum and minimum values ​​of the i-th indicator, respectively; if all evaluation units of the indicator have the same value, its standardized value is set to 0.5.

[0038] The formula for calculating the negative indicator is as follows:

[0039]

[0040] In the formula, Let j be the standardized value of the j-th indicator in the i-th evaluation unit, ranging from 0 to 1. The original value of the indicator. , These are the maximum and minimum values ​​of the i-th indicator, respectively; if all evaluation units of an indicator have the same value, its standardized value is set to 0.5.

[0041] S32. Calculate the comprehensive weight using the comprehensive weighting method: calculate the objective weight using the entropy weighting method. First, calculate the proportion of the i-th evaluation unit under the j-th indicator using the formula. Then, the information entropy is calculated using the formula. Finally, the objective weights are calculated using a formula. ; Analytic Hierarchy Process (AHP) for calculating subjective weights First, invite 5-10 experts in water conservancy, planning, and emergency response to construct a judgment matrix for indicators within the same criterion layer using a 1-9 scaling method; then calculate the maximum eigenvalue of the judgment matrix. The corresponding feature vectors are then normalized through a consistency test, and the global subjective weights are calculated by combining the criterion layer weights. The weights for the criteria layer include urban space, lifeline engineering, and rapid recovery, with weights of 0.35, 0.4, and 0.25 respectively; the comprehensive weight is calculated by integrating subjective and objective weights. This reflects the overall impact of the indicators on resilience;

[0042] The formula for calculating the proportion of the i-th evaluation unit under the j-th indicator is as follows:

[0043]

[0044] The formula for calculating information entropy is as follows:

[0045]

[0046] The formula for calculating objective weights is as follows:

[0047]

[0048] The formula for calculating the consistency test is as follows:

[0049]

[0050] In the formula, , RI represents the number of criteria-level indicators, and RI is the average random consistency index.

[0051] The formula for calculating the overall weight is as follows:

[0052]

[0053] S33, the TOPSIS method quantifies resilience level, first evaluating the parameter matrix in S31. Forward processing yields For the normalized matrix Standardization is performed to obtain the standardized decision matrix. Then, the standardized decision matrix is... After weighting, a weighted standardized decision matrix is ​​obtained. Determine the ideal solution With negative ideal solution ;

[0054] Through formula , Calculate the geometric distance between each evaluation unit and the positive and negative ideal solutions; use the formula... Calculate the similarity in toughness. The range is 0~1, with values ​​closer to 1 indicating stronger toughness; the natural discontinuity grading method is used to classify... Classified as the first, that is Second, that is Third, that is Fourth, that is Fifth, namely Five resilience levels;

[0055] Positive processing is for positive indicators = Regarding negative indicators = ;

[0056] A standardized decision matrix is ​​constructed, and the normalized parameter matrix is ​​standardized using a normalization formula, as shown below:

[0057]

[0058] The weighted standardized decision matrix is ​​shown below:

[0059] =

[0060] In the formula, The overall weights calculated for S32 are... This is a standardized decision matrix.

[0061] As a further preferred embodiment of the present invention, S5 includes the following steps:

[0062] S51. Regional measures are formulated. For old urban areas where the terrain is below the preset threshold and the difficulty of pipeline renovation is greater than the preset threshold, water storage tanks are built to collect surface water in the area. For areas where the upstream flow is greater than the preset threshold and the pipeline drainage pressure is greater than the preset threshold, drainage pumping stations are added to improve the area's pumping capacity. In areas where the pipeline drainage capacity is insufficient and the pipeline density is lower than the preset threshold, the drainage capacity is improved by increasing the pipe size and increasing the pipe slope.

[0063] S52. Verification of the effectiveness of the measures: Adjust the corresponding parameters of the measures in the IFMS / Urban one-dimensional coupled urban flood simulation model, re-simulate 10-50 year return period rainfall and calculate the resilience approximation. The proportion of resilience enhancement evaluation units should be ≥80%.

[0064] S53. Closed-loop optimization: Combined with dynamic adjustment measures of urban development planning, a closed-loop resilience enhancement process of assessment, identification, improvement and verification is formed to continuously optimize urban flood resilience.

[0065] The beneficial effects of this invention are as follows:

[0066] This invention presents a method for assessing urban flood resilience and optimizing the configuration of flood control facilities. It integrates resilience theory, disaster chain theory, and hydrological and hydrodynamic mechanisms, overcoming the limitations of traditional assessments. It constructs an indicator system covering multiple dimensions—natural, social, economic, and emergency—and the entire process from pre-disaster to post-disaster, comprehensively reflecting the core connotations of urban flood resilience. Through a one- or two-dimensional coupled IFMS / Urban model, it completes design storm calculations, data processing, modeling simulations, and coupled verification according to specifications, providing high-precision flood parameters with a ≥85% hit rate for flood-prone areas, thus solving the problem of disconnect between static assessments and actual processes. It employs a comprehensive weighting method to balance subjective experience and objective data patterns, combined with the TOPSIS method to intuitively quantify resilience levels, and uses spatial analysis to accurately locate weak areas, avoiding biases from single methods. Targeted improvement measures are formulated for different regions, and the optimization effect is verified through model iterations, forming a dynamic closed-loop improvement mechanism. This method is applicable to long-term flood resilience management in cities of different sizes, playing a core role in scientifically and quantitatively guiding urban flood resilience assessment, optimizing facility layout, and improving control effectiveness.

[0067] Other advantages, objectives, and features of the invention will be set forth in the following description and will be apparent to those skilled in the art in some respects, or may be learned by practice of the invention. The objectives and other advantages of the invention can be realized and obtained through the following description. Attached Figure Description

[0068] To make the objectives, technical solutions, and beneficial effects of this invention clearer, the following figures are provided for illustration:

[0069] Figure 1 is a flowchart of the urban flooding resilience evaluation method of the present invention;

[0070] Figure 2 is Figure 1 System structure diagram of urban waterlogging resilience evaluation index system in medium-sized cities;

[0071] Figure 3 yes Figure 1 A schematic diagram showing the distribution of urban spatial resilience levels in the study area.

[0072] Figure 4 yes Figure 1 A schematic diagram showing the distribution of resilience levels across different urban spatial dimensions in the study area.

[0073] Figure 5 yes Figure 1 A schematic diagram showing the distribution of lifeline engineering resilience levels in the central study area;

[0074] Figure 6 yes Figure 1 A schematic diagram showing the distribution of resilience levels across different zones within the lifeline engineering dimension of the central study area;

[0075] Figure 7 yes Figure 1 A schematic diagram illustrating the distribution of rapid recovery dimension resilience levels in the study area.

[0076] Figure 8 yes Figure 1 A schematic diagram showing the distribution of resilience levels across different regions of the rapid recovery dimension in the study area.

[0077] Figure 9 yes Figure 1 A schematic diagram showing the distribution of the overall resilience level in the study area;

[0078] Figure 10 yes Figure 1 A schematic diagram showing the overall distribution of resilience levels across different zones in the central study area;

[0079] Figure 11 is Figure 1 Average values ​​of various indicators for different resilience levels in the study area;

[0080] Figure 12 is Figure 1 A schematic diagram of LISA clustering in the fifth resilience region of the study area; Detailed Implementation

[0081] like Figures 1-12 As shown, this invention discloses a method for evaluating urban flood resilience based on multi-dimensional indicators and a coupled model.

[0082] The study area is the southwestern part of Suqian City, covering an area of ​​164.28 km², with a terrain that slopes from north to south (ground elevation 15-58 m). The average annual precipitation is 915 mm, with the flood season (June to September) accounting for 70% of the annual rainfall. The area includes 14 rivers such as the ancient Yellow River and the Ximinbian River. The drainage network in the old city is designed for a 1-2 year return period, which aligns with the application scenario of this invention.

[0083] like Figure 1 As shown in the figure, a method for evaluating urban flood resilience based on multi-dimensional indicators and a coupled model in this specific embodiment includes the following steps: Step 1, obtaining the urban flood resilience evaluation indicator system, includes the following steps:

[0084] Step 1: Based on resilience theory and disaster chain theory, identify the influencing factors of urban flooding resilience, select 20 indicators from three dimensions of urban space, lifeline engineering, and rapid recovery, and construct a three-level urban flooding resilience evaluation index system.

[0085] Step 11, Implementation of Urban Spatial Dimension Indicators. ArcGIS 10.8 software was used to process the basic data, obtaining data for seven indicators: green coverage rate, water coverage rate, topographic slope, ground elevation, distance to river network, distance to flood-prone areas, and surface water depth. Of these indicators, the first six are positive, while surface water depth is a negative indicator. The data was processed using 3... The principle test showed no abnormalities.

[0086] Step 12, Implementation of Lifeline Engineering Dimension Indicators. Eight indicators are obtained through POI data and infrastructure data: distance to medical institutions, road network density, distance to transportation hubs, drainage network density, distance to telecommunications maintenance stations, distance to power maintenance stations, water supply network density, and population density. Among these, road network density, drainage network density, and water supply network density are positive indicators, while the rest are negative indicators.

[0087] Step 13, implement rapid recovery of dimensional indicators. Integrate multi-source data to obtain 5 indicators: nighttime light, distance to drainage pumping station, distance to drainage sluice gate, distance to emergency rescue station, and distance to disaster relief shelter. The first 3 are positive indicators, and the last 2 are negative indicators. Combined with the data from steps 11 and 12, a total of 20 comprehensive indicators (x1~x) are finally formed. 20 ).

[0088] Step 14: A three-level indicator system is formed, including the target layer (urban flood resilience), the criterion layer (3 dimensions), and the indicator layer (20 indicators), as shown in Table 1.

[0089] Table 1. Evaluation Index System for Urban Flood Resilience

[0090]

[0091]

[0092] Step 2: Construct a coupled one-dimensional and two-dimensional urban flood simulation model using IFMS / Urban. First, collect hydrological, meteorological, topographic, river network, and land use data for the study area. Then, use IFMS / Urban software to perform the following operations: import DEM data to generate a two-dimensional surface grid, and perform elevation uplift processing on the grid in the building area; digitize the river network and drainage network, generalize river cross-sections, pipe parameters, and manhole attributes to construct a one-dimensional model; establish a one-dimensional and two-dimensional coupling relationship using orifice connections, and set parameters such as permeable depression storage capacity, Manning coefficient, and soil infiltration rate; input measured rainfall data to validate the model, requiring a simulation hit rate of ≥85% for flood-prone points and an average relative error of ≤20% for the maximum inundation depth; simulate rainfall scenarios with different return periods and output data on surface water depth and distribution of flood-prone points.

[0093] Step 21: Collect basic data for the study area, including model construction data such as 5m resolution DEM, 5m resolution land use, measured river cross sections, pipeline network data of manholes and pipes, and building distribution; boundary condition data such as hourly rainfall at rain gauge stations and river water level-flow relationship; and model validation data such as measured rainfall inundation depth and distribution of flood-prone areas.

[0094] Step 22: Basic data preprocessing and one-dimensional model construction. Data preprocessing is completed using ArcGIS software, including DEM depression smoothing and slope analysis, land use reclassification, and river network topology inspection; rivers are generalized and key cross-sections are densified to construct a one-dimensional river network model; the drainage system is generalized into a "manhole-pipe-outlet" system, and parameters such as pipe cross-sectional shape and pipe diameter are defined to construct a one-dimensional drainage network model.

[0095] Step 23: Sub-catchment division and 2D model construction. First, the area is divided into 21 drainage zones according to the distribution of rivers and pipe networks. Then, 2244 sub-catchments are divided using the Thiessen polygon method (farmland / village areas without pipe networks are not divided). Parameters such as sub-catchment area, slope, and proportion of impermeable areas are extracted. 395469 unstructured quadrilateral grids are divided into the modeling area (excluding river channels). The main roads are used as grid control lines. The grid elevation is assigned based on the 2.5m precision DEM. The grid elevation of the building is raised to construct the 2D surface model.

[0096] Step 24: Coupling of one-dimensional and two-dimensional models. In the "Coupled Modeling Module" of IFMS / Urban software, the hydraulic coupling relationship between the one-dimensional river network-pipeline model and the two-dimensional surface model is established using the manhole as the connection node and the orifice connection method, so as to realize the dynamic interaction between the river, pipeline and surface water volume.

[0097] Step 25, Model Validation: Using two measured rainfall events as input, the flooding situation of 21 flood-prone points was simulated. The flooding rates of the flood-prone points in the two rainfall events were 90.5% and 85.7%, respectively. The average relative errors of the measured and simulated flooding depths were 18.1% and 19.3%, respectively (absolute error < 0.1m), thus validating the model accuracy.

[0098] Step 26: Model parameter setting and calculation scheme determination. Set core modeling parameters, including Manning coefficient for permeable areas, Manning coefficient for impermeable areas, Manning coefficient for pipes, and maximum soil infiltration rate. Use the local official rainstorm intensity formula to calculate the design rainstorm intensity for 10-year, 20-year, 30-year, and 50-year return periods. Use the Chicago rain pattern to distribute the rainfall time history from 120 to 180 minutes and determine the simulation calculation scheme.

[0099] Furthermore, this specific example uses Suqian City as an example, so the intensity of the aforementioned rainstorm is as follows:

[0100]

[0101] In the formula, i represents the intensity of the rainstorm (mm / min); T represents the return period (years); and t represents the duration of rainfall (min).

[0102] Calculate the 3-hour design rainfall for 10a to 50a return periods: 105.1 mm cumulative rainfall for 10a return period (maximum rainfall intensity 3.1 mm / min), 121.3 mm for 20a return period (maximum rainfall intensity 3.6 mm / min), 130.8 mm for 30a return period (maximum rainfall intensity 3.9 mm / min), and 142.8 mm for 50a return period (maximum rainfall intensity 4.2 mm / min); allocate the 180-minute rainfall duration according to the Chicago rainfall pattern (peak rainfall coefficient 0.4), with the peak rainfall occurring at the 72nd minute.

[0103] Step 27: Extraction of key flooding parameters. Key flooding parameters are extracted from the coupled simulation results: Surface water depth: Under a 50-year return period rainfall, the water depth in the eastern old urban area is 0.8~2.24m, while in the western rural area it is <0.3m; Distribution of flood-prone points: Of the 21 measured flood-prone points, 19 were simulated to have water accumulation; Pipeline overflow nodes: Of the 3079 inspection wells, 58 experienced overflow. These parameters are used for subsequent quantification of the "surface water depth" and "distance from flood-prone points" indicators.

[0104] Step 3, Quantitative assessment of urban flood resilience: The weight of each indicator is determined by a comprehensive weighting method combining the analytic hierarchy process (subjective weighting) and the entropy weighting method (objective weighting). The resilience score of each assessment unit is calculated by the top-inferior solution distance method (TOPSIS). The resilience level is divided into five levels: first resilience, second resilience, third resilience, fourth resilience, and fifth resilience. Spatial autocorrelation analysis is used to identify resilience clusters and weak areas.

[0105] Step 31, Standardization of Indicator Data: The raw data of the 20 indicators (including the waterlogging parameters extracted in Step 2 and the distance and density parameters obtained from GIS spatial analysis) are divided into evaluation units of 100m×100m, totaling 17317 units. The raw data of the 20 indicators are then organized into a 17317×20 matrix. The extreme value method is used to eliminate the influence of dimensions. Positive indicators (such as green coverage rate) are standardized using the formula, and negative indicators (such as ground water depth) are also standardized using the formula. The water supply network density has a consistent value in some areas, so its standardized value is set to 0.5. After standardization, the data range is 0~1, with no outliers.

[0106] Furthermore, the evaluation matrix is ​​as follows:

[0107] For n evaluation objects and m evaluation indicators, establish an evaluation matrix.

[0108]

[0109] In the formula, i = 1, 2, ..., n; j = 1, 2, ..., m.

[0110] Furthermore, the formula for calculating the positive indicator is as follows:

[0111]

[0112] In the formula, This is the standardized value (range 0~1) of the j-th indicator for the i-th evaluation unit. The original value of the indicator. , These are the maximum and minimum values ​​of the i-th indicator, respectively; if all evaluation units of an indicator have the same value, its standardized value is set to 0.5.

[0113] Furthermore, the formula for calculating the negative index is as follows:

[0114]

[0115] In the formula, This is the standardized value (range 0~1) of the j-th indicator for the i-th evaluation unit. The original value of the indicator. , These are the maximum and minimum values ​​of the i-th indicator, respectively; if all evaluation units of an indicator have the same value, its standardized value is set to 0.5.

[0116] Step 32: Calculate the comprehensive weight using the comprehensive weighting method. Calculate the objective weight using the entropy weight method. First, calculate the proportion of the i-th evaluation unit under the j-th indicator using the formula. Then, the information entropy is calculated using the formula. Finally, the objective weights are calculated using a formula. Calculating subjective weights using the analytic hierarchy process (AHP). First, invite 5-10 experts in water conservancy, planning, and emergency response to construct a judgment matrix for indicators within the same criterion layer using a 1-9 scaling method; then calculate the maximum eigenvalue of the judgment matrix. The corresponding feature vectors are then normalized through a consistency test. The global subjective weights are calculated by combining the criterion layer weights (0.35 for urban space, 0.4 for lifeline engineering, and 0.25 for rapid recovery). The overall weight is calculated by combining subjective and objective factors according to the formula. The results show: ground elevation =0.081, distance from medical institution =0.074, Distance from disaster shelter =0.072 and terrain slope Dj=0.071 are the top four weighted indicators.

[0117] Furthermore, the formula for calculating the proportion of the i-th evaluation unit under the j-th indicator is as follows:

[0118]

[0119] If in the formula ,set up Avoid logarithms being meaningless.

[0120] Furthermore, the formula for calculating the information entropy is as follows:

[0121]

[0122] In the formula .

[0123] Furthermore, the objective weight calculation formula is as follows:

[0124]

[0125] In the formula Must meet .

[0126] Furthermore, the consistency test calculation formula is as follows:

[0127]

[0128] In the formula , RI represents the number of criteria-level indicators, and RI is the average random consistency index.

[0129] Furthermore, the formula for calculating the comprehensive weight is as follows:

[0130]

[0131] Step 33: Quantify resilience level using the TOPSIS method. First, evaluate the parameter matrix from step S31. Forward processing yields For the normalized matrix Standardization is performed to obtain the standardized decision matrix. Then, the standardized decision matrix is... After weighting, a weighted standardized decision matrix is ​​obtained. Determine the ideal solution. With negative ideal solution ;

[0132] Through formula , Calculate the geometric distance between each evaluation unit and the positive and negative ideal solutions; use the formula... Calculate the similarity in toughness. The range is 0~1, with values ​​closer to 1 indicating stronger toughness; the natural discontinuity grading method is used to classify... Classified as the first, that is Second, that is Third, that is Fourth, that is Fifth, namely Five resilience levels, spatial distribution of each resilience level as follows: Figures 3-10 As shown.

[0133] Furthermore, the positiveization process is applied to positive indicators. = Regarding negative indicators = ;

[0134] Furthermore, the construction of the standardized decision matrix involves standardizing the forward-oriented parameter matrix using a normalization formula, as shown below:

[0135]

[0136] Furthermore, the weighted standardized decision matrix is ​​shown below:

[0137] =

[0138] In the formula, The comprehensive weights calculated in step 32 are... For standardized decision matrices;

[0139] Step 4: Identify the fifth resilience cluster region through Moran's I index and LISA cluster analysis, and combine the index weights and the mean values ​​of different resilience levels to screen key influencing indicators such as drainage network density and distance from flood-prone areas.

[0140] Step 41, Fifth, Resilience Region Identification. Load resilience proximity data in ArcGIS software. Based on raster data, Moran's I index was calculated to be 0.708 (P<0.05), indicating significant spatial clustering of resilience levels. LISA cluster analysis generated a cluster map, identifying a "fifth cluster" region accounting for 18.29%, which is a key area for remediation. Overlaying land use data revealed that the fifth resilience region is mostly located in old urban areas and low-lying areas, such as... Figure 12 As shown.

[0141] Step 42, Key Indicator Identification. Calculate the weighted standard deviation of each indicator. and screening Key indicators; comprehensive and The study identified several key indicators that significantly impact resilience assessment results, including green coverage rate, distance from river network, distance from flood-prone areas, sluice gate scheduling capacity, and drainage network density. These indicators are considered the core entry points for resilience enhancement.

[0142] Furthermore, the formula for calculating the weighted standard deviation is as follows:

[0143]

[0144] In the formula The standard deviation of the standardized values ​​of the indicator. This is the mean.

[0145] Step 5: Based on the different target level years and different catchment area needs of the study area, formulate measures such as adding sponge city facilities, upgrading pipe networks, and raising riverbanks. Execute steps 2-4 in a cycle to verify the effectiveness of the measures, ensuring that more than 80% of the area has improved resilience, thus forming a closed-loop resilience improvement process.

[0146] Step 51: Formulate measures for different areas. For old urban areas where the terrain is below the preset threshold and the difficulty of pipe network renovation exceeds the preset threshold, water storage tanks should be constructed to collect surface water. Furthermore, for areas where upstream flow exceeds the preset threshold and pipe network drainage pressure exceeds the preset threshold, consider adding drainage pumping stations to improve the area's pumping capacity. In areas where pipe network drainage capacity is below the preset threshold and pipe network density is below the preset threshold, drainage capacity should be improved by increasing pipe size and slope.

[0147] Step 52: Verify the effectiveness of the measures. Adjust the corresponding parameters of the measures in the IFMS / Urban model. Resimulate using a 50-year return period for rainfall, and repeat steps 3-4 to calculate the resilience approximation. The results showed that the resilience level of 85.8% of the regions improved after the measures were implemented, and the measures achieved their intended effect.

[0148] Step 53: Closed-loop optimization. Combined with dynamic adjustment measures of urban development planning, a closed-loop resilience enhancement process of "assessment-identification-improvement-verification" is formed to ensure continuous optimization of urban flood resilience.

[0149] Finally, it should be noted that the above preferred embodiments are only used to illustrate the technical solutions of the present invention and are not intended to limit it. Although the present invention has been described in detail through the above preferred embodiments, those skilled in the art should understand that various changes can be made to it in form and detail without departing from the scope defined by the claims of the present invention.

Claims

1. A method for evaluating urban flood resilience based on multi-dimensional indicators and a coupled model, characterized in that, Includes the following steps: S1. Select indicators from three dimensions—urban space, lifeline engineering, and rapid recovery—to construct a three-level urban flood resilience evaluation index system; S2. Construct a one-dimensional coupled urban flood simulation model using IFMS / Urban. Collect basic data of the study area and preprocess it. Perform the following operations using IFMS / Urban software: import DEM data to generate a two-dimensional surface grid, and perform elevation uplift processing on the grid in the building area; digitize the river network and drainage network, generalize river cross-sections, pipeline parameters, and manhole attributes, and construct a one-dimensional model; establish a one-dimensional coupling relationship using orifice connection method and set core modeling parameters; input measured rainfall data to simulate rainfall scenarios with different return periods, and extract key flooding parameters from the coupled simulation results for quantification of evaluation indicators. S3. Quantitative assessment of urban waterlogging resilience level: Based on key waterlogging parameters and a three-level urban waterlogging resilience evaluation index system, the weight of each index is determined by a comprehensive weighting method combining the analytic hierarchy process and the entropy weight method. The resilience score of each assessment unit is calculated by the superior-inferior solution distance method, and the resilience level is divided into five levels: first resilience, second resilience, third resilience, fourth resilience and fifth resilience. It also includes S4 and S5; S4. Identify the fifth resilience cluster region through Moran's I index and LISA cluster analysis, and screen key influencing indicators by combining indicator weights and the mean values ​​of indicators for different resilience levels. S5. Based on the identified fifth resilience cluster area and key influencing indicators, formulate resilience enhancement measures and cyclically execute S2~S4 to verify the effectiveness of the measures, forming a closed-loop resilience enhancement process. The three-level urban flood resilience evaluation index system described in S1 includes a target layer, a criterion layer, and an indicator layer. The target layer is urban flood resilience, the criterion layer includes three dimensions: urban space, lifeline engineering, and rapid recovery, and the indicator layer includes 20 indicators. The urban spatial dimension includes seven indicators: green coverage rate, water coverage rate, topographic slope, ground elevation, distance from river network, distance from flood-prone areas, and ground water depth. Among them, green coverage rate, water coverage rate, topographic slope, ground elevation, distance from river network, and distance from flood-prone areas are positive indicators, while ground water depth is a negative indicator. The lifeline project dimension includes eight indicators: distance to medical institutions, road network density, distance to transportation hubs, drainage network density, distance to telecommunications maintenance stations, water supply network density, distance to power maintenance stations, and population density. Among these, road network density, drainage network density, and water supply network density are positive indicators, while distance to medical institutions, distance to transportation hubs, telecommunications maintenance stations, power maintenance stations, and population density are negative indicators. The rapid recovery dimension includes five indicators: nighttime lighting, distance to drainage pumping stations, distance to sluice gates, distance to emergency rescue stations, and distance to disaster shelters; among them, nighttime lighting, distance to drainage pumping stations, and distance to sluice gates are positive indicators, while distance to emergency rescue stations and distance to disaster shelters are negative indicators. S4 includes the following steps: S41. Fifth resilience region identification: Load resilience proximity in ArcGIS software. Given raster data, calculate Moran's I index. and This indicates a significant spatial clustering of toughness levels; cluster maps are generated using LISA cluster analysis to identify the concentrated distribution area of ​​the fifth toughness unit; S42. Key indicator identification and calculation of weighted standard deviation for each indicator. Compared to the first toughness, With the fifth resilience, that is Regional indicator mean, filtering Key influencing indicators with large and significant differences in means; the criteria for determining significant differences in means are as follows: ; The formula for calculating the weighted standard deviation is as follows: In the formula The standard deviation of the standardized values ​​of the indicator. The mean, This is the overall weight.

2. The urban flood resilience evaluation method based on multi-dimensional indicators and a coupled model according to claim 1, characterized in that: S2 includes the following steps: S21. Collect basic data for the study area: Use a combination of field surveys and data collection to obtain basic data, including basic geographic and infrastructure data required for model building, boundary condition data required for model operation, and measured data required for model accuracy verification. S22. Basic data preprocessing and one-dimensional model construction: The collected basic data is standardized and preprocessed using geographic information system software, including digital elevation model filling and smoothing, slope analysis, land use classification system regularization, and verification of river network and pipeline network topology. A one-dimensional river network model is constructed using the river channel generalization method in IFMS / Urban software. A one-dimensional drainage pipeline network model system is formed by structural generalization of the pipeline system, consisting of nodes, pipe segments, and discharge outlets. S23. Sub-catchment delineation and 2D model construction: Primary drainage zones are delineated according to the distribution of river network and pipe network, and then secondary sub-catchments are subdivided based on the Thiessen polygon method to extract the core topography and underlying surface parameters of the sub-catchments; a 2D surface model is constructed using unstructured quadrilateral meshing technology, and the mesh scale is determined according to the principle of denser meshing around the river channel and conventional meshing in other areas. Main roads are used as meshing constraints, and the mesh elevation is assigned in combination with topographic data. The effect of buildings on water flow obstruction is simulated by adjusting the elevation. S24. Coupling of one-dimensional and two-dimensional models: In the coupling module of IFMS / Urban software, inspection wells are used as key connection nodes of one-dimensional and two-dimensional models. Through orifice connection technology, a two-way hydraulic coupling relationship between the one-dimensional river network-pipeline network model and the two-dimensional surface model is established to realize the dynamic water volume interaction simulation among the river, pipeline and surface. S25. Model Validation: Selecting actual rainfall events as input conditions, a combination of matching degree verification of flood-prone areas and quantitative error analysis of flood depth is adopted. By comparing simulation results with actual data, the model's ability to replicate flood processes and its accuracy and reliability are verified. S26. Model parameter settings and calculation scheme determination: Based on the characteristics of the underlying surface and the hydraulic motion law of the study area, the core hydraulic parameters of the permeable zone, impermeable zone, pipe network and river channel are set; the design rainfall intensity of different return periods is calculated using the rainfall intensity formula published for the study area; the rainfall time history is allocated by adopting the Chicago rainfall pattern through the spatiotemporal distribution characteristics of regional rainfall; and the model simulation calculation scheme is determined. S27. Extraction of key parameters for urban flooding: Extract key parameters for urban flooding in the study area at the current level of the year from the coupled simulation results, including the surface water depth and the distribution of flood-prone points, for the quantification of evaluation indicators.

3. The urban flood resilience evaluation method based on multi-dimensional indicators and a coupled model according to claim 2, characterized in that: S3 includes the following steps: S31. Standardization of indicator data: The original data of 20 indicators are organized into an n×m evaluation matrix in 100m×100m evaluation units, where n is the number of evaluation units and m=20. The original data includes the key waterlogging parameters extracted in S2 and the distance and density parameters obtained from GIS spatial analysis. The extreme value method is used to eliminate the influence of dimensions. The standardized values ​​of both positive and negative indicators are calculated according to the formula. For n evaluation objects and m evaluation indicators, establish an evaluation matrix. The evaluation matrix is ​​shown below: In the formula, i = 1, 2, ..., n; j = 1, 2, ..., m; The formula for calculating the positive indicator is as follows: In the formula, Let j be the standardized value of the j-th indicator in the i-th evaluation unit, ranging from 0 to 1. The original value of the indicator. , These are the maximum and minimum values ​​of the i-th indicator, respectively; if all evaluation units of the indicator have the same value, its standardized value is set to 0.

5. The formula for calculating the negative indicator is as follows: In the formula, Let j be the standardized value of the j-th indicator in the i-th evaluation unit, ranging from 0 to 1. The original value of the indicator. , These are the maximum and minimum values ​​of the i-th indicator, respectively; If all evaluation units of a certain indicator have the same value, its standardized value is set to 0.5; S32. Calculate the comprehensive weight using the comprehensive weighting method: calculate the objective weight using the entropy weighting method. First, calculate the proportion of the i-th evaluation unit under the j-th indicator using the formula. ; Then calculate the information entropy using the formula. Finally, the objective weights are calculated using a formula. ; Analytic Hierarchy Process (AHP) for calculating subjective weights First, invite 5-10 experts in water conservancy, planning, and emergency response to construct a judgment matrix for indicators within the same criterion layer using a 1-9 scaling method; then calculate the maximum eigenvalue of the judgment matrix. The corresponding feature vectors are then normalized through a consistency test, and the global subjective weights are calculated by combining the criterion layer weights. The weights for the criteria layer include urban space, lifeline engineering, and rapid recovery, with weights of 0.35, 0.4, and 0.25 respectively; the comprehensive weight is calculated by integrating subjective and objective weights. This reflects the overall impact of the indicators on resilience; The formula for calculating the proportion of the i-th evaluation unit under the j-th indicator is as follows: The formula for calculating information entropy is as follows: The formula for calculating objective weights is as follows: The formula for calculating the consistency test is as follows: In the formula, , RI represents the number of criteria-level indicators, and RI is the average random consistency index. The formula for calculating the overall weight is as follows: S33, the TOPSIS method quantifies resilience level, first evaluating the parameter matrix in S31. Forward processing yields For the normalized matrix Standardization is performed to obtain the standardized decision matrix. Then, the standardized decision matrix is... After weighting, a weighted standardized decision matrix is ​​obtained. Determine the ideal solution With negative ideal solution ; Through formula , Calculate the geometric distance between each evaluation unit and the positive and negative ideal solutions; use the formula... Calculate the similarity in toughness. The range is 0~1, with values ​​closer to 1 indicating stronger toughness; the natural discontinuity grading method is used to classify... Classified as the first, that is Second, that is Third, that is Fourth, that is Fifth, namely Five resilience levels; Positive processing is for positive indicators = Regarding negative indicators = ; A standardized decision matrix is ​​constructed, and the normalized parameter matrix is ​​standardized using a normalization formula, as shown below: The weighted standardized decision matrix is ​​shown below: = In the formula, The overall weights calculated for S32 are... This is a standardized decision matrix.

4. The urban flood resilience evaluation method based on multi-dimensional indicators and a coupled model according to claim 1, characterized in that: S5 includes the following steps: S51. Regional measures are formulated. For old urban areas where the terrain is below the preset threshold and the difficulty of pipeline renovation is greater than the preset threshold, water storage tanks are built to collect surface water in the area. For areas where the upstream flow is greater than the preset threshold and the pipeline drainage pressure is greater than the preset threshold, drainage pumping stations are added to improve the area's pumping capacity. In areas where the pipeline drainage capacity is insufficient and the pipeline density is lower than the preset threshold, the drainage capacity is improved by increasing the pipe size and increasing the pipe slope. S52. Verification of the effectiveness of the measures: Adjust the corresponding parameters of the measures in the IFMS / Urban one-dimensional coupled urban flood simulation model, re-simulate 10-50 year return period rainfall and calculate the resilience approximation. The proportion of resilience enhancement evaluation units should be ≥80%. S53. Closed-loop optimization: Combined with dynamic adjustment measures of urban development planning, a closed-loop resilience enhancement process of assessment, identification, improvement and verification is formed to continuously optimize urban flood resilience.

Images