A method and system for optimizing an urban agglomeration ecological network from the perspective of toughness improvement

Abstract

The present application relates to the technical field of urban ecological planning and ecological network optimization, and discloses a kind of urban agglomeration ecological network optimization method and system under the perspective of resilience promotion.The method obtains the ecological environment, terrain, land use and social economic data of research area, identifies ecological source and constructs ecological resistance surface;Ecological network is constructed based on ecological source and resistance surface, and ecological corridor and its obstruction area are extracted and abstracted into topological network;According to the area of obstruction area and land use type, it is converted into edge weight;Under the constraint of preset repair quantity, traverse repair combination, respectively calculate the repair cost and the network resilience of each combination after updating;With the minimum cost and the maximum resilience as the goal, pareto frontier screening is carried out, and the repair strategy solution set under different input conditions is obtained.The present application can unify the quantification of corridor obstruction degree, repair cost and network resilience promotion, and provide decision basis for urban agglomeration ecological network resilience promotion and corridor repair priority determination.

Description

Technical Field

[0001] This invention relates to the field of urban ecological planning and ecological network optimization technology, specifically to a method and system for optimizing urban agglomeration ecological networks from the perspective of resilience enhancement. Background Technology

[0002] Rapid urbanization has exacerbated the problems of natural habitat shrinkage, ecological patch fragmentation, and ecological corridor blockage. Ecological networks, by identifying ecological source areas and constructing ecological corridors, connect isolated ecological patches into a holistic ecological security pattern, serving as a crucial technical means to maintain regional ecological security, protect biodiversity, and enhance landscape connectivity. Existing ecological network construction typically involves steps such as ecological source area extraction, ecological resistance surface generation, and ecological corridor identification. However, the initially identified ecological corridors are often fragmented by high-resistance land types such as roads, farmland, and construction land, resulting in limitations on species migration, energy flow, and the transfer of ecological processes.

[0003] To ensure the effective operation of ecological networks, it is necessary to identify and restore obstructive areas within ecological corridors. Restoring obstructive areas typically involves land use adjustments or spatial controls, resulting in real costs such as agricultural output losses, opportunity costs for construction land, and restrictions on secondary and tertiary industry activities. At the urban agglomeration scale, the large number of ecological corridors, the complexity of obstructive area types, and significant constraints on restoration investment mean that simply identifying key corridors or key obstructive points is insufficient to meet practical decision-making needs. Therefore, how to select restoration combinations under constraints of limited funds or a limited number of restoration projects to maximize the improvement of ecological network functions is a pressing technical problem to be solved in the optimization of urban agglomeration ecological networks.

[0004] The existing technology CN112883471B discloses a dynamic simulation method for the structural resilience of urban ecological networks. It abstracts ecological source areas as nodes and ecological corridors as edges to construct a topological network. It uses network efficiency and maximum connectivity as resilience indicators, and simulates attacks by removing important nodes proportionally and calculates the resilience threshold. The current scheme has the following shortcomings: (1) It uses building height as the single source of edge resistance value. The degree of different land use types that hinder ecological flow is mixed in the edge weight of the topological network and cannot be distinguished. When multiple corridors share the same obstacle area, the area is repeatedly included in each corridor. (2) It calculates resilience using unweighted network efficiency and maximum connectivity. The size of the obstacle area and the difference in land use type do not participate in the resilience calculation. The loss of farmland output and the opportunity cost of secondary and tertiary industries involved in the restoration of the obstacle area are not within the evaluation framework. There is no quantitative correlation between the restoration cost and the resilience improvement effect. (3) It simulates the output of the resilience threshold by removing nodes, which only reflects the overall resilience level of the network. It cannot quantitatively compare the resilience improvement effect of different restoration combinations under the given restoration quantity constraint. It is difficult to quantitatively determine the priority of corridor restoration.

[0005] Therefore, there is an urgent need for an urban agglomeration ecological network optimization method that integrates the identification and area allocation of hindering areas by land use type, the dual-track accounting of restoration costs and the weighted multidimensional resilience quantification, and the traversal and screening of restoration combinations, so as to provide a quantitative decision-making basis for determining the priority of corridor restoration under the constraint of limited restoration resources. Summary of the Invention

[0006] The purpose of this invention is to overcome the shortcomings of existing technologies and provide a method and system for optimizing urban agglomeration ecological networks from the perspective of resilience enhancement. This invention can solve the problems of insufficient coupling between restoration costs and network resilience enhancement, and the difficulty in quantitatively determining corridor restoration priorities in existing ecological network optimization methods.

[0007] To achieve the above objectives, in a first aspect, the present invention provides a method for optimizing urban agglomeration ecological networks from the perspective of resilience enhancement, comprising the following steps: S1. Obtain ecological environment data, topographic data, land use data, and socio-economic data of the study area, and identify the ecological source area based on the ecological environment data and land use data, and construct the ecological resistance surface based on the topographic data and land use data; S2. Using the ecological source areas as network nodes and the ecological resistance surfaces as the basis for ecological flow resistance, construct an ecological network for the study area; extract the ecological corridors connecting different ecological source areas and the obstruction areas located on the ecological corridors in the ecological network, and abstract the ecological network into a topological network; S3. Calculate the area and land use type of the obstructing area according to the ecological corridor, and determine the edge weight of the corresponding ecological corridor based on the area of ​​the obstructing area; S4. For the ecological corridor, under the preset constraint of the number of ecological corridors to be restored, generate multiple ecological corridor restoration combinations; for each ecological corridor restoration combination, calculate the restoration cost based on the obstruction area area and land use type corresponding to the ecological corridor to be restored in the combination. S5. For each ecological corridor restoration combination, update the topology network according to the edge weight changes of the restored ecological corridors in the combination, and calculate the resilience of the updated ecological network. S6. Based on the restoration cost and the ecological network resilience, the multiple ecological corridor restoration combinations are screened to obtain a solution set of ecological network resilience optimization restoration strategies.

[0008] Furthermore, the ecological environment data includes habitat quality assessment data; the topographic data includes elevation data and slope data calculated from the elevation data; the land use data includes at least one of forest, shrub, grassland, water body, farmland, construction land and impermeable surface; and the socio-economic data includes at least one of population density, output value of secondary industry in districts and counties, output value of tertiary industry in districts and counties, grain crop price and grain crop yield.

[0009] Furthermore, identifying the ecological source area in step S1 includes: S11. Perform binarization processing on the land use data, setting land types with ecological connectivity functions as the foreground and other land types as the background to obtain the binarized land use data; S12. Perform morphological spatial pattern analysis on the binarized land use data to obtain the ecological core area; S13. Calculate the habitat quality of the study area to obtain the habitat quality results; S14. Take the spatial intersection of the ecological core area and the habitat quality results as the candidate ecological source area; remove patches with an area smaller than a preset threshold from the candidate ecological source areas to obtain the ecological source area.

[0010] Furthermore, the construction of the ecological resistance surface in step S1 includes: selecting slope, elevation, and land use type as ecological resistance factors; classifying each ecological resistance factor into levels and assigning resistance values; and weighting and superimposing the resistance values ​​according to the weight of each ecological resistance factor to obtain a comprehensive ecological resistance surface.

[0011] Furthermore, the construction of the ecological network in step S2 includes: S21. Determine the ecological corridors between the ecological source areas based on circuit theory or the minimum cumulative resistance theory; S22. Identify obstructions affecting ecological flow along the ecological corridor; S23. Abstract the ecological source area as a node of the topological network, and the ecological corridor as an edge of the topological network to construct an ecological network.

[0012] Furthermore, in step S3, the obstruction areas are divided into first-class obstruction areas and second-class obstruction areas according to land use type. The first-class obstruction areas are farmland obstruction areas, and the second-class obstruction areas are obstruction areas of other land use types besides farmland. The edge weight of each ecological corridor is determined based on the area of ​​the first-class obstruction area and the area of ​​the second-class obstruction area corresponding to the ecological corridor.

[0013] Furthermore, the edge weight of each ecological corridor is positively correlated with the total area of ​​the obstruction region corresponding to that ecological corridor, or it is obtained by normalizing the weighted sum of the areas of the first type of obstruction region and the areas of the second type of obstruction region corresponding to that ecological corridor.

[0014] Furthermore, the repair cost in step S4 includes the farmland output cost corresponding to the first type of obstruction area and the secondary and tertiary industry output costs corresponding to the second type of obstruction area; The farmland output cost corresponding to the first type of obstruction area is calculated according to the following formula: , , in, For correction factor, Prices of major local agricultural products, This accounts for the local main crop yield. To hinder the area area, The equivalent value of ecosystem services per unit area. The cost of farmland output corresponding to the first type of obstacle area; The output costs of the secondary and tertiary industries corresponding to the second type of barrier area are calculated as follows: Cost models are constructed using data on the GDP of the secondary and tertiary industries, the proportion of construction land, and population density of all districts and counties within the urban agglomeration. , Obtained through the above model , , Based on the corresponding parameters of this urban agglomeration, a production cost model for the secondary and tertiary industries corresponding to the second type of barrier region is established: , in, For districts and counties The sum of the gross domestic product of the secondary and tertiary industries For districts and counties Construction land ratio For districts and counties Population density Category II Obstruction Area The sum of the output costs of the secondary and tertiary industries. Category II Obstruction Area The proportion of construction land, Category II Obstruction Area population density For total factor productivity, The output elasticity of the proportion of construction land, The output elasticity is the ratio of population density.

[0015] Furthermore, the ecological network resilience in step S5 is calculated using the following formula: , in, As a network size indicator, As a weighted network density metric, As a weighted network centrality metric, This is a weighted global network efficiency metric.

[0016] Furthermore, let's consider a topology network. , For the number of nodes, Let be the number of edges. Let be the accessibility weight of the ecological corridor between node i and node j. Let the sum of the shortest path weights from node i to node j be . Let be the node strength of node i. If the maximum value among all access weights is given, then: , , , , Furthermore, the access weight It is obtained by reverse conversion of the edge weight of the ecological corridor, so that the larger the area of ​​the obstruction area, the smaller the access weight, and the corresponding access weight increases after the ecological corridor is restored.

[0017] Furthermore, the generation of multiple ecological corridor restoration combinations in step S4 specifically includes: when the given number of ecological corridor restorations is K, selecting K ecological corridors from all ecological corridors in the ecological network to form an ecological corridor restoration combination, and traversing all ecological corridor restoration combinations that satisfy the constraint of the number of ecological corridor restorations.

[0018] Furthermore, the screening method in step S6 is Pareto front screening, and the specific steps include: S61. The restoration cost and ecological network resilience corresponding to each ecological corridor restoration combination are combined into a two-dimensional evaluation point; S62. Remove other evaluation points from the evaluation points, such that the repair cost is no higher than that evaluation point and the ecological network resilience is no lower than that evaluation point, to obtain the remaining evaluation points; S63. The ecological corridor restoration combination corresponding to the retained evaluation points is used as the solution set of ecological network resilience optimization and restoration strategy.

[0019] Secondly, the present invention provides an urban agglomeration ecological network optimization system from the perspective of resilience enhancement, the system being based on any one of the above methods, comprising: The data acquisition module is used to acquire ecological environment data, topographic data, land use data, and socio-economic data of the study area, and to identify the ecological source areas based on the ecological environment data and land use data; and to construct ecological resistance surfaces based on the topographic data and land use data. The network construction module is used to construct an ecological network of the study area using the ecological source areas as network nodes and the ecological resistance surfaces as the basis for ecological flow resistance; extract the ecological corridors connecting different ecological source areas and the obstruction areas located on the ecological corridors in the ecological network, and abstract the ecological network into a topological network. The edge weight determination module is used to calculate the area and land use type of the obstructing area according to the ecological corridor, and determine the edge weight of the corresponding ecological corridor based on the area of ​​the obstructing area. The restoration combination and cost module is used to generate multiple ecological corridor restoration combinations under the preset ecological corridor restoration quantity constraint; for each ecological corridor restoration combination, the restoration cost is calculated based on the obstruction area area and land use type corresponding to the restored ecological corridor in the ecological corridor restoration combination. The resilience calculation module is used to update the topology network based on the edge weight changes of the restored ecological corridors in each ecological corridor restoration combination, and to calculate the resilience of the updated ecological network. The screening module filters the multiple ecological corridor restoration combinations based on the restoration cost and the ecological network resilience to obtain a solution set of optimized restoration strategies for ecological network resilience.

[0020] Thirdly, the present invention provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, it implements the steps of the method for optimizing urban agglomeration ecological networks from the perspective of resilience enhancement.

[0021] Compared with the prior art, the present invention has at least the following beneficial effects: 1. The ecological corridor obstruction area is classified into farmland obstruction area and other land use obstruction area according to land use type, and the area is statistically analyzed and converted into topological network edge weights. When multiple corridors share the same obstruction area, the area is allocated according to the number of corridors, so that the different obstruction degree of different land use types can be mapped to the topological network edge.

[0022] 2. Establish a dual-track repair cost model for farmland output costs and secondary and tertiary industry output costs. After converting edge weights into accessibility weights, incorporate four weighted indicators—network size, weighted network density, network centrality, and global network efficiency—to jointly calculate resilience, thereby establishing a quantitative correlation between repair cost input and resilience improvement effect.

[0023] 3. Under a given constraint on the number of repairs, traverse the repair combinations and perform Pareto front screening to output the solution set of repair strategies under different cost input conditions, so that the priority of corridor repair can be quantitatively determined. Attached Figure Description

[0024] Figure 1 This is a flowchart of an urban agglomeration ecological network optimization method from the perspective of resilience enhancement, according to Embodiment 1 of the present invention. Figure 2 This is a schematic diagram of the ecological network topology of Embodiment 2 of this aspect; Figure 3 This is a schematic diagram of the cost-toughness Pareto front screening results for Embodiment 2 of the present invention; Figure 4 This is a schematic diagram of an urban agglomeration ecological network optimization system from the perspective of resilience enhancement, as described in Embodiment 3 of the present invention. Detailed Implementation

[0025] To make the objectives, technical solutions, and advantages of this application clearer, the following detailed description is provided in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the scope of this application.

[0026] Example 1 like Figure 1 As shown in the preferred embodiment of the present invention, a method for optimizing urban agglomeration ecological networks from the perspective of resilience enhancement includes the following steps: S1. Obtain ecological environment data, topographic data, land use data, and socio-economic data of the study area, and identify the ecological source area based on the ecological environment data and land use data, and construct the ecological resistance surface based on the topographic data and land use data; In this embodiment, the ecological environment data includes habitat quality assessment data; the topographic data includes elevation data and slope data calculated from the elevation data; the land use data includes forests, shrubs, grasslands, water bodies, farmland, construction land, and impermeable surfaces; and the socio-economic data includes population density, output value of the secondary industry in the district / county, output value of the tertiary industry in the district / county, grain crop prices, and grain crop yields. In this embodiment, identifying the ecological source area specifically involves: S11. Perform binarization processing on the land use data, setting forests, shrubs, grasslands, and water bodies as the foreground and other land types as the background to obtain the binarized land use data; S12. Perform morphological spatial pattern analysis on the binarized land use data to obtain the ecological core area; S13. Calculate the habitat quality of the study area and obtain the habitat quality results; S14. Take the spatial intersection of the ecological core area and the habitat quality results as the candidate ecological source area; remove patches with an area smaller than a preset threshold from the candidate ecological source areas to obtain the ecological source area.

[0027] In this embodiment, the construction of the ecological resistance surface is specifically as follows: slope, elevation, and land use type are selected as ecological resistance factors; each ecological resistance factor is classified into levels and assigned a resistance value; the resistance values ​​are weighted and superimposed according to the weight of each ecological resistance factor to obtain a comprehensive ecological resistance surface; each resistance factor is classified into levels using the natural breakpoint method, the equidistant grading method, or the expert assignment method.

[0028] S2. Using the ecological source areas as network nodes and the ecological resistance surfaces as the basis for ecological flow resistance, construct an ecological network for the study area; extract the ecological corridors connecting different ecological source areas and the obstruction areas located on the ecological corridors in the ecological network, and abstract the ecological network into a topological network; In this embodiment, ecological corridors between the ecological source areas are determined based on circuit theory or the minimum cumulative resistance theory, and obstructive regions affecting ecological flow are identified on the ecological corridors. The ecological source areas are abstracted as nodes in a topological network, and the ecological corridors are abstracted as edges in the topological network.

[0029] S3. Calculate the area and land use type of the obstructing area according to the ecological corridor, and determine the edge weight of the corresponding ecological corridor based on the area of ​​the obstructing area; In this embodiment, the area of ​​obstructing regions and their land use types are statistically analyzed according to ecological corridors, dividing the obstructing regions into two categories: Category I obstructing regions and Category II obstructing regions. Category I obstructing regions are farmland obstructing regions, while Category II obstructing regions are obstructing regions with other land use types besides farmland. The edge weight of each ecological corridor is determined based on the areas of the corresponding Category I and Category II obstructing regions. The edge weight can be positively correlated with the total area of ​​the obstructing regions, or it can be obtained by weighted summation of the areas of the two categories of obstructing regions followed by normalization.

[0030] When the same obstruction area is traversed by two or more ecological corridors, the area of ​​the obstruction area is evenly distributed according to the number of ecological corridors passing through it. The distributed area is then included in the obstruction area area of ​​the corresponding ecological corridor. This process avoids shared obstruction areas being counted repeatedly by multiple corridors, improving the accuracy of edge weight calculation.

[0031] S4. For the ecological corridor, under the preset constraint of the number of ecological corridors to be restored, generate multiple ecological corridor restoration combinations; for each ecological corridor restoration combination, calculate the restoration cost based on the obstruction area area and land use type corresponding to the ecological corridor to be restored in the combination. In this embodiment, under the constraint of a preset number of ecological corridors to be restored, a preset number of ecological corridors are selected from all ecological corridors in the ecological network to form a restoration combination, and all restoration combinations that meet the constraint are traversed. For each restoration combination, the farmland output cost corresponding to the first type of obstruction area and the secondary and tertiary industry output costs corresponding to the second type of obstruction area are calculated respectively, and the two are used as the restoration cost of the restoration combination.

[0032] The agricultural output cost for Category I obstructed areas is calculated based on the area of ​​the obstructed area, the local price of the main crop, the local yield of the main crop, and the equivalent value of ecosystem services per unit area. The calculation formula is as follows: , , in, For correction factor, Prices of major local agricultural products, This accounts for the local main crop yield. To hinder the area area, The equivalent value of ecosystem services per unit area. The cost of farmland output corresponding to the first type of obstacle area; Cost models are constructed using data on the GDP of the secondary and tertiary industries, the proportion of construction land, and population density of all districts and counties within the urban agglomeration. , Obtained through the above model , , Based on the corresponding parameters of this urban agglomeration, a production cost model for the secondary and tertiary industries corresponding to the second type of barrier region is established: , in, For districts and counties The sum of the gross domestic product of the secondary and tertiary industries For districts and counties Construction land ratio For districts and counties Population density Category II Obstruction Area The sum of the output costs of the secondary and tertiary industries. Category II Obstruction Area The proportion of construction land, Category II Obstruction Area population density For total factor productivity, The output elasticity of the proportion of construction land, The output elasticity is the ratio of population density.

[0033] S5. For each ecological corridor restoration combination, update the topology network according to the edge weight changes of the restored ecological corridors in the combination, and calculate the resilience of the updated ecological network. In this embodiment, since the edge weight in step S5 is used to characterize the degree of obstruction of the ecological corridor, the larger the edge weight, the higher the degree of obstruction; when calculating the resilience of the ecological network, the edge weight is converted into the corridor access weight in reverse, so that the access weight is smaller when the area of ​​the obstructed area is larger, and the access weight increases after the ecological corridor is repaired.

[0034] Let the topology network be defined. , For the number of nodes, Let be the number of edges. Let be the accessibility weight of the ecological corridor between node i and node j. Let the sum of the shortest path weights from node i to node j be . Let be the node strength of node i. The ecological network resilience is the maximum value among all access weights. Network size indicators Weighted network density index Network centrality index and global network efficiency metrics Jointly determined: , , , , , in, Used to characterize the overall accessibility efficiency between any two nodes in a network. Used to characterize the scale of network nodes, edges, and their access weights. Used to characterize the connection density of weighted networks This is used to characterize the centralization level of node strength relative to the maximum accessibility weight. Through the above model, changes in edge weights caused by the same repair combination can be transformed into changes in ecological network resilience.

[0035] S6. Based on the restoration cost and the ecological network resilience, the multiple ecological corridor restoration combinations are screened to obtain a solution set of ecological network resilience optimization restoration strategies.

[0036] In this embodiment, the repair cost and ecological network resilience corresponding to each repair combination are used to form a two-dimensional evaluation point. Other evaluation points with repair costs no higher than the evaluation point and ecological network resilience no lower than the evaluation point are eliminated. The repair combinations corresponding to the evaluation points that have not been eliminated are used as the solution set of ecological network resilience optimization repair strategies.

[0037] Example 2 This embodiment takes a certain urban agglomeration as the research area. Based on the method described in Embodiment 1, it further illustrates the optional methods in specific implementation, and the steps are as follows: Step S1. Obtain ecological environment data, topographic data, land use data, and socio-economic data of the study area, and identify ecological source areas based on the ecological environment data and land use data, and construct ecological resistance surfaces based on the topographic data and land use data; In this embodiment, a specific urban agglomeration is used as the study area. Land use data, elevation data, population density data, output value of secondary and tertiary industries in districts and counties, grain crop prices, and grain crop yields are obtained for this urban agglomeration. Land use data can be CLCD land cover data, LUCC land use data, or other data with equivalent spatial representation capabilities. Elevation data is used to calculate slope, and population density and industrial output value are used to estimate the socioeconomic output costs of the second type of barrier area. Land use data was reclassified, with ecological land types such as forests, shrubs, grasslands, and water bodies set as the foreground, and construction land, bare land, and other non-ecological land types set as the background. Morphological spatial pattern analysis was used to obtain the ecological core area. Habitat quality was calculated using a habitat quality assessment model. Core areas with high habitat quality and spatial continuity were selected as candidate ecological source areas, while scattered patches with areas smaller than a preset threshold were removed to form the final ecological source areas. An ecological resistance surface was constructed based on elevation, slope, and land use type. Each resistance factor was divided into multiple levels using the natural discontinuity method, the equidistant grading method, or the expert assignment method, and resistance values ​​were assigned according to the ease of ecological mobility. After weighting and superimposing the resistance factors, a comprehensive ecological resistance surface was obtained.

[0038] Step S2. Using the ecological source areas as network nodes and the ecological resistance surfaces as the basis for ecological flow resistance, construct an ecological network for the study area; extract the ecological corridors connecting different ecological source areas and the obstruction areas located on the ecological corridors in the ecological network, and abstract the ecological network into a topological network; In this embodiment, as Figure 2 As shown, ecological source areas are used as nodes, and ecological corridors are constructed using comprehensive ecological resistance surfaces to obtain the spatial pattern of the ecological network. Obstructive areas along the ecological corridors are identified and statistically analyzed, recording the area and land use type of each obstructive area within each corridor. The ecological network is abstracted as a topological network, with each ecological source area corresponding to a node and each ecological corridor corresponding to an edge.

[0039] Step S3. Calculate the area and land use type of the obstructing area according to the ecological corridor, and determine the edge weight of the corresponding ecological corridor based on the area of ​​the obstructing area; In this embodiment, the edge weight is calculated based on the area of ​​the obstruction region corresponding to each edge of the topology network. The larger the area of ​​the obstruction region, the stronger the ecological flow obstruction that the ecological corridor was subjected to before restoration. If an obstruction region is traversed by multiple corridors at the same time, the area of ​​the obstruction region is divided by the number of corridors traversed, and the area after apportionment is included in each corridor.

[0040] Step S4. For the ecological corridor, under the preset constraint of the number of ecological corridors to be restored, generate multiple ecological corridor restoration combinations; for each ecological corridor restoration combination, calculate the restoration cost based on the obstruction area area and land use type corresponding to the ecological corridor to be restored in the combination. In this embodiment, the number of ecological corridors to be restored is set to K, and K corridors are selected from all ecological corridors to form a restoration combination. The cost of farmland obstruction areas is calculated using the farmland output cost model, while the cost of obstruction areas of other land use types is calculated using the secondary and tertiary industry output cost models. For obstruction areas spanning multiple districts and counties, the costs are calculated separately according to the area ratio of the obstruction area in each district and county, and then summed.

[0041] Step S5. For each ecological corridor restoration combination, update the topology network according to the edge weight changes of the restored ecological corridors in the combination, and calculate the resilience of the updated ecological network. In this embodiment, when repairing a corridor, the corresponding edge weights are updated based on the state of the obstructed area after repair, and the updated edge weights are converted into accessibility weights before recalculating network resilience. For each repair combination, its total repair cost and the repaired ecological network resilience are calculated.

[0042] S6. Based on the restoration cost and the ecological network resilience, the multiple ecological corridor restoration combinations are screened to obtain a solution set of ecological network resilience optimization restoration strategies.

[0043] In this embodiment, as Figure 3 As shown, all restoration combinations are mapped to a cost-resilience two-dimensional space. The Pareto dominance relation is used to screen non-dominated solutions to obtain the Pareto front. Each solution on the Pareto front represents a restoration strategy with better resilience enhancement effect at the corresponding cost input level. The output results include the ecological corridor number to be restored, the area of ​​the obstruction area corresponding to each ecological corridor to be restored, the land use type of the obstruction area, the restoration cost, the ecological network resilience value, and the corresponding Pareto front position.

[0044] Example 3 like Figure 4 As shown, a preferred embodiment of the present invention provides an urban agglomeration ecological network optimization system from the perspective of resilience enhancement, comprising: The data acquisition module is used to acquire ecological environment data, topographic data, land use data, and socio-economic data of the study area, and to identify ecological source areas based on the ecological environment data and land use data, and to construct ecological resistance surfaces based on the topographic data and land use data. The network construction module is used to construct an ecological network of the study area with the ecological source area as the network node and the ecological resistance surface as the basis of ecological flow resistance, extract the ecological corridors connecting different ecological source areas and the obstruction areas located on the ecological corridors, and abstract the ecological network into a topological network including nodes and edges. The edge weight determination module is used to calculate the area and land use type of the obstructing area according to the ecological corridor, and determine the edge weight of the corresponding ecological corridor based on the area of ​​the obstructing area. The restoration combination and cost module is used to generate multiple ecological corridor restoration combinations under a given constraint on the number of ecological corridors to be restored, and to calculate the restoration cost for each ecological corridor restoration combination based on the area of ​​the obstruction area and the land use type corresponding to the restored ecological corridor. The resilience calculation module is used to update the topology network for each ecological corridor restoration combination based on the edge weight changes of the restored ecological corridors, and to calculate the resilience of the updated ecological network. The screening module is used to screen multiple ecological corridor restoration combinations with the optimization objectives of minimizing restoration costs and maximizing ecological network resilience, and obtain the solution set of ecological network resilience optimization restoration strategies under different cost input conditions.

[0045] This embodiment correlates the area of ​​the obstructing region of the ecological corridor with the land use type and the edge weight of the topological network, so that the degree of obstruction of the corridor can be included in the network resilience calculation. When multiple corridors share the obstructing area, the area is allocated according to the number of corridors. The output costs of farmland and the output costs of secondary and tertiary industries are included in the restoration cost model, so that the restoration plan takes into account both ecological benefits and practical input. Under the given restoration quantity constraint, the restoration combination is traversed and Pareto front screening is performed to obtain the solution set of restoration strategies under different cost input conditions, providing a quantitative decision-making basis for improving the ecological network resilience of urban agglomerations and determining the priority of corridor restoration.

[0046] Example 4 This embodiment also provides a preferred embodiment of a computer-readable storage medium storing a computer program for a pedestrian identification method. When the computer program is processed, it implements the steps of the urban agglomeration ecological network optimization method from the perspective of resilience enhancement.

[0047] In summary, this invention identifies obstructive areas in ecological corridors and classifies their areas by land use type, transforming these obstructive areas into topological network edge weights. When multiple corridors share the same obstructive area, the area is allocated according to the number of corridors, thus mapping the degree of obstruction to the network topological edges in a differentiated manner based on land use type. By establishing a dual-track restoration cost model that combines farmland output costs and secondary and tertiary industry output costs, edge weights are converted inversely into accessibility weights and incorporated into four weighted indicators—network size, weighted network density, network centrality, and global network efficiency—to jointly calculate resilience, establishing a quantitative correlation between restoration cost input and resilience improvement effects. By traversing restoration combinations under given restoration quantity constraints and performing Pareto front screening, the invention outputs a set of restoration strategy solutions under different cost input conditions, providing a quantitative decision-making basis for improving the resilience of urban agglomeration ecological networks and determining corridor restoration priorities.

[0048] The above description is only a preferred embodiment of the present invention. It should be noted that for those skilled in the art, several improvements and substitutions can be made without departing from the technical principles of the present invention, and these improvements and substitutions should also be considered within the scope of protection of the present invention.

Claims

1. A method for optimizing urban agglomeration ecological networks from the perspective of resilience enhancement, characterized in that, Includes the following steps: S1. Obtain ecological environment data, topographic data, land use data, and socio-economic data of the study area, and identify the ecological source area based on the ecological environment data and land use data, and construct the ecological resistance surface based on the topographic data and land use data; S2. Using the ecological source areas as network nodes and the ecological resistance surfaces as the basis for ecological flow resistance, construct an ecological network for the study area; extract the ecological corridors connecting different ecological source areas and the obstruction areas located on the ecological corridors in the ecological network, and abstract the ecological network into a topological network; S3. Calculate the area and land use type of the obstructing area according to the ecological corridor, and determine the edge weight of the corresponding ecological corridor based on the area of ​​the obstructing area; S4. For the ecological corridor, under the preset constraint of the number of ecological corridors to be restored, generate multiple ecological corridor restoration combinations; for each ecological corridor restoration combination, calculate the restoration cost based on the obstruction area area and land use type corresponding to the ecological corridor to be restored in the combination. S5. For each ecological corridor restoration combination, update the topology network according to the edge weight changes of the restored ecological corridors in the combination, and calculate the resilience of the updated ecological network. S6. Based on the restoration cost and the ecological network resilience, the multiple ecological corridor restoration combinations are screened to obtain a solution set of ecological network resilience optimization restoration strategies.

2. The method for optimizing urban agglomeration ecological networks according to claim 1, characterized in that, The ecological environment data mentioned in step S1 includes habitat quality assessment data; the topographic data includes elevation data and slope data calculated from the elevation data; the land use data includes at least one of forest, shrub, grassland, water body, farmland, construction land and impermeable surface; the socio-economic data includes at least one of population density, output value of secondary industry in districts and counties, output value of tertiary industry in districts and counties, grain crop price and grain crop yield.

3. The method for optimizing urban agglomeration ecological networks according to claim 1, characterized in that, The steps in step S1 to identify the ecological source area include: S11. Perform binarization processing on the land use data, setting land types with ecological connectivity functions as the foreground and other land types as the background to obtain the binarized land use data; S12. Perform morphological spatial pattern analysis on the binarized land use data to obtain the ecological core area; S13. Calculate the habitat quality of the study area to obtain the habitat quality results; S14. Take the spatial intersection of the ecological core area and the habitat quality results as the candidate ecological source area; remove patches with an area smaller than a preset threshold from the candidate ecological source areas to obtain the ecological source area.

4. The method for optimizing urban agglomeration ecological networks according to claim 3, characterized in that, Step S2, which involves constructing the ecological network, includes: S21. Determine the ecological corridors between the ecological source areas based on circuit theory or the minimum cumulative resistance theory; S22. Identify obstructions affecting ecological flow along the ecological corridor; S23. Abstract the ecological source area as a node of the topological network, and the ecological corridor as an edge of the topological network to construct an ecological network.

5. The method for optimizing urban agglomeration ecological networks according to claim 1, characterized in that, The obstruction areas mentioned in step S3 are divided into Class I obstruction areas and Class II obstruction areas according to land use type. Class I obstruction areas are farmland obstruction areas, and Class II obstruction areas are land use types other than farmland obstruction areas. The edge weight of each ecological corridor is determined based on the area of ​​the Class I obstruction area and the area of ​​the Class II obstruction area corresponding to the ecological corridor.

6. The method for optimizing urban agglomeration ecological networks according to claim 5, characterized in that, The repair costs mentioned in step S4 include the farmland output costs corresponding to the first type of obstruction area and the secondary and tertiary industry output costs corresponding to the second type of obstruction area; The farmland output cost corresponding to the first type of obstruction area is calculated according to the following formula: ， ， in, For correction factor, Prices of major local agricultural products, This accounts for the local main crop yield. To hinder the area area, The equivalent value of ecosystem services per unit area. The cost of farmland output corresponding to the first type of obstacle area; The output costs of the secondary and tertiary industries corresponding to the second type of barrier area are calculated as follows: Cost models are constructed using data on the GDP of the secondary and tertiary industries, the proportion of construction land, and population density of all districts and counties within the urban agglomeration. ， Obtained through the above model , , Based on the corresponding parameters of this urban agglomeration, a production cost model for the secondary and tertiary industries corresponding to the second type of barrier region is established: ， in, For districts and counties The sum of the gross domestic product of the secondary and tertiary industries, For districts and counties Construction land ratio For districts and counties Population density Category II Obstruction Area The sum of the output costs of the secondary and tertiary industries. Category II Obstruction Area The proportion of construction land, Category II Obstruction Area population density For total factor productivity, The output elasticity of the proportion of construction land, The output elasticity is the ratio of population density.

7. The method for optimizing urban agglomeration ecological networks according to claim 1, characterized in that, The ecological network resilience described in step S5 is calculated using the following formula: ， Wherein, NS is the network size index, ND is the weighted network density index, NC is the weighted network centrality index, and NV is the weighted global network efficiency index.

8. The method for optimizing urban agglomeration ecological networks according to claim 1, characterized in that, The screening method in step S6 is Pareto front screening, and the specific steps include: S61. The restoration cost and ecological network resilience corresponding to each ecological corridor restoration combination are combined into a two-dimensional evaluation point; S62. Remove other evaluation points from the evaluation points, such that the repair cost is no higher than that evaluation point and the ecological network resilience is no lower than that evaluation point, to obtain the remaining evaluation points; S63. The ecological corridor restoration combination corresponding to the retained evaluation points is used as the solution set of ecological network resilience optimization and restoration strategy.

9. A system for optimizing urban agglomeration ecological networks from the perspective of resilience enhancement, characterized in that, The system, based on the method according to any one of claims 1-8, includes: The data acquisition module is used to acquire ecological environment data, topographic data, land use data, and socio-economic data of the study area, and to identify the ecological source areas based on the ecological environment data and land use data; and to construct ecological resistance surfaces based on the topographic data and land use data. The network construction module is used to construct an ecological network of the study area using the ecological source areas as network nodes and the ecological resistance surfaces as the basis for ecological flow resistance; extract the ecological corridors connecting different ecological source areas and the obstruction areas located on the ecological corridors in the ecological network, and abstract the ecological network into a topological network. The edge weight determination module is used to calculate the area and land use type of the obstructing area according to the ecological corridor, and determine the edge weight of the corresponding ecological corridor based on the area of ​​the obstructing area. The restoration combination and cost module is used to generate multiple ecological corridor restoration combinations under the preset ecological corridor restoration quantity constraint; for each ecological corridor restoration combination, the restoration cost is calculated based on the obstruction area area and land use type corresponding to the restored ecological corridor in the ecological corridor restoration combination. The resilience calculation module is used to update the topology network based on the edge weight changes of the restored ecological corridors in each ecological corridor restoration combination, and to calculate the resilience of the updated ecological network. The screening module filters the multiple ecological corridor restoration combinations based on the restoration cost and the ecological network resilience to obtain a solution set of optimized restoration strategies for ecological network resilience.

10. A computer device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the steps of the method according to any one of claims 1 to 8.

Images