A method, medium and device for evaluating the resilience of a road traffic disaster reduction route selection in a strong earthquake mountainous area

Abstract

The present application relates to railway and highway route selection technical field, disclose a kind of strong earthquake mountainous area land transportation disaster mitigation route selection flexibility evaluation method, medium and equipment, specific method includes: according to the comprehensive geographic information in research area, establish geographic information model in research area;And based on the construction of geographic information model seismic risk field model;The selected line scheme is scattered into multiple section units, and the flexibility of each section unit of the selected line scheme is calculated by combining the geographic information model in research area and the seismic risk field model, and the flexibility integral corresponding to each section unit of the selected line scheme is obtained;Based on flexibility integral, determine the comprehensive flexibility level index of the selected line scheme;According to the minimum flexibility integral in the selected line scheme and the comprehensive flexibility level index, the disaster mitigation route selection flexibility of the selected line scheme is evaluated to obtain evaluation result.The present application can realize the comprehensive evaluation of the safety and post-disaster recovery capacity of different route selection schemes under the action of strong earthquake disaster.

Description

Technical Field

[0001] This invention relates to the field of railway and highway route selection technology, and specifically discloses a method, medium and equipment for assessing the resilience of land transportation in mountainous areas during strong earthquakes to mitigate disasters. Background Technology

[0002] As the construction of new-generation high-speed railways and highways gradually advances into the complex mountainous areas of western China, the external environment for line construction is becoming increasingly severe. On the one hand, these complex mountainous areas are mostly located in earthquake-prone zones, with extremely complex and uncertain geological conditions, and frequent natural disasters such as landslides, mudslides, and flash floods. On the other hand, the complex terrain imposes strict constraints on the route alignment, coupled with multiple environmental requirements such as ecological protection and land use control, resulting in a significant reduction in the feasible space for the route and inevitably requiring it to traverse some areas prone to strong earthquakes.

[0003] Existing research on line disaster prevention and safety has the following limitations: 1. Primarily post-evaluation: This approach focuses on post-construction (post-evaluation) analyses after the project is completed or operational, making it difficult to provide timely guidance for route selection and optimization during the route selection phase. 2. Single evaluation object: Most methods focus on the safety or vulnerability analysis of individual structures or local sections such as bridges, tunnels, and roadbeds, and lack a systematic assessment of the entire line in the whole process of pre-disaster, during-disaster and post-disaster. 3. Qualitative analysis of indicators: Existing studies mostly use qualitative methods such as expert scoring or experience grading, which are insufficient in terms of quantitative analysis and make it difficult to reflect the differences in the overall resilience level and post-disaster recovery capacity of the line schemes under strong earthquakes.

[0004] Especially in mountainous areas prone to strong earthquakes, existing technologies struggle to conduct a comprehensive, systematic, and quantitative evaluation of the disaster mitigation capabilities and resilience of different route options when traversing high-risk disaster areas during the route selection phase. This makes it difficult to effectively shift from "post-event evaluation" to "pre-event assessment" and to provide scientific support for the optimal selection of route corridors and the comparison of different options. Summary of the Invention

[0005] The purpose of this invention is to provide a method, medium, and equipment for assessing the resilience of mountainous land transportation disaster reduction routes under strong earthquakes, enabling a comprehensive evaluation of the safety and post-disaster recovery capabilities of different route selection schemes. The specific technical solution is as follows: This invention provides a method, medium, and equipment for assessing the resilience of land transportation routes in mountainous areas prone to strong earthquakes, comprising the following steps: S1: Establish a geographic information model of the study area based on comprehensive geographic information within the study area; and construct an earthquake risk field model based on the geographic information model. S2: Discretize the candidate route scheme into multiple segment-level units, and combine the geographic information model and seismic risk field model of the study area to perform resilience calculations on each segment-level unit of the candidate route scheme to obtain the resilience integral corresponding to each segment-level unit of the candidate route scheme. S3: Based on the resilience integrals corresponding to each segment-level unit of the candidate route scheme obtained in step S2, determine the overall resilience level index of the candidate route scheme; evaluate the disaster mitigation and route selection resilience of the candidate route scheme according to the minimum resilience integral and the overall resilience level index of the candidate route scheme to obtain the evaluation result.

[0006] Optionally, in S1, the geographic information model includes topographic information, geological information, distribution of rescue stations, and existing road network; the seismic risk field model includes at least the spatial distribution of ground peak ground acceleration (PGA) within the study area, used to characterize the intensity level of ground motion at different locations.

[0007] Optionally, S2 specifically includes: S2.1. Based on the existing three-dimensional route space model, construct a segment-level discrete space model for the candidate route schemes, specifically: S2.1.1 Obtain the spatial location of the candidate route scheme based on the existing three-dimensional route space model; S2.1.2 Discretize the candidate route schemes along the centerline of the route at preset intervals to obtain multiple segment-level units, and assign corresponding segment-level attributes to each segment-level unit, including: structure type, bridge height, bridge span, tunnel depth, whether the tunnel crosses a fault, embankment filling height, road cut excavation depth, and distance to the nearest rescue station. S2.2 Based on the segment-level discrete spatial model of the candidate route scheme, and combined with the geographic information model and seismic risk field model of the study area, the resilience of each segment-level unit of the candidate route scheme is calculated to obtain the resilience integral corresponding to each segment-level unit of the candidate route scheme.

[0008] Optionally, S2.2 includes: S2.2.1. Based on the seismic risk field model and the segment-level attributes of each segment-level unit, the structural vulnerability of each segment-level unit is calculated using a unified log-normal vulnerability curve paradigm; where: structural vulnerability includes the probability of occurrence of different structural damage states, specifically: no damage The probability of occurrence, minor injury The probability of occurrence, moderate injury The probability of occurrence and the severity of injury The probability of occurrence; S2.2.2 Calculate the remaining function of each segment-level unit based on the structural vulnerability of each segment-level unit; S2.2.3 Calculate the toughness integral of each segment-level element based on the residual function of each segment-level element, specifically: The equivalent recovery curves of the candidate route schemes are calculated based on the geographic information model, the seismic risk field model, and the residual function of each segment-level unit. ; For equivalent recovery curve Within the control time interval Integrating over time, we obtain the first... The integral of the segment's toughness is calculated using the following formula: ; in: For the first The toughness integral of the segment-level unit, To control time, For a moment.

[0009] Optionally, in S2.2.1, a specific seismic motion level is given. No damage probability of occurrence The specific calculation formula is as follows: ; in: The condition is damaged. Minor injury The probability of exceeding; It is the standard normal distribution function; PAG is the ground motion intensity index in the earthquake risk field model; This indicates that the damage to the structure exceeds the level of minor damage. hour The middle threshold; Indicates confirmation The logarithmic standard deviation at that time; minor injury probability of occurrence The specific calculation formula is as follows: ; in: Moderate injury The probability of exceeding; This indicates that the damage to the structure exceeds the level of moderate damage. hour The middle threshold; Indicates confirmation The logarithmic standard deviation at that time; Moderate injury probability of occurrence The specific calculation formula is as follows: ; in: For severe injury The probability of exceeding; The structural damage condition exceeds severe damage. hour The middle threshold; To indicate certainty The logarithmic standard deviation at that time; Severe injury probability of occurrence The specific calculation formula is as follows: ; In S2.2.1, for tunnel sections in the candidate route that do not cross fault zones and have a burial depth exceeding 400 meters, there is no damage. probability of occurrence The probability of occurrence of other damage states is 0, with a value of 1.

[0010] Optionally, S2.2.2 includes: ① Based on the structural vulnerability of each segment-level unit, the structural damage probability is mapped to the corresponding operational state, specifically: No damage Supports unlimited speed operation Minor injury Dual-line operation, slight speed limit Moderate injury Single-line operation, significantly reduced speed Severe injury Corresponding closed line For tunnel sections crossing fault zones, when their damage condition is moderate... At that time, its corresponding operational status is directly determined to be closed. ; ② Set the speed factor according to the operating status corresponding to each segment level element. and line usage factor Specifically: velocity factor Set as the ratio of the post-earthquake restricted speed to the design speed; Set to the ratio of the throughput capacity of the number of lines currently open in the operating state to the throughput capacity when both lines are open; ③ Based on the velocity factor and line usage factor The remaining functional coefficient for each operating state is calculated using the following formula: ; in: Damaged state The residual function coefficient at that time, for , , , Any damage state; ④ Based on residual function coefficient The remaining functionality of each segment-level unit is calculated using the following formula: ; in: For remaining functions, For the first Segment-level unit damage status The probability of occurrence.

[0011] Optionally, in step S2.2.3, the equivalent recovery curve of the candidate route scheme is calculated based on the geographic information model, the seismic risk field model, and the residual function of each segment-level unit. Specifically, it includes: ① The rescue accessibility index is calculated based on the geographic information model and the earthquake risk field model. The specific calculation formula is as follows: ; in, This indicates the actual distance of the segment-level unit from the nearest rescue station; Indicates the base reachable distance; ② Utilize the rescue accessibility index The typical repair time is corrected to obtain the corrected typical repair time. The specific formula is as follows: ; in: This is the typical repair time, which is the time required to restore 80% of the functionality; ③ An exponential recovery curve is used to characterize each damage state, and the recovery rate is calibrated using the corrected typical repair time to obtain the recovery curve for each damage state. The specific expression is as follows: ; in: Damaged state Down The function of time; It is a natural constant; Damaged state Recovery rate at the following levels; ④ Based on the principle of equivalent area, the recovery curves under the four injury states are weighted and fitted into a single equivalent recovery curve. Where: Let the first The segment is in a damaged state after the earthquake. The probability is The fitting process is as follows: ; in: The fitted recovery rate; Single equivalent recovery curve The specific expression is as follows: .

[0012] Optionally, S3 includes: S3.1. The equivalent downtime is defined based on a single equivalent recovery curve of segment-level units, and the specific formula is as follows: ; in, The downtime equivalent to functional degradation within the controlled time period; S3.2. Take the segment-level unit corresponding to the minimum resilience integral in the candidate route scheme as the most unfavorable segment, and then calculate the resilience integral of the most unfavorable segment. To determine the resilience level of the proposed route options, specifically: At that time, the toughness level was excellent, and it entered S3.3; At this point, the toughness level is good, and it enters S3.3; When the resilience level is low, the evaluation result of the proposed route is that it is not recommended as the preferred option. in: Set a threshold for the first one. Set a threshold for the second one; S3.3 Calculate the overall resilience level index of the candidate route scheme. The specific calculation formula is as follows: ; in: To satisfy the toughness integral of the entire line segment The cumulative length of the line segment, Set a threshold for the third; The total length of the line; S3.4, Based on the overall resilience level index of the entire line The overall resilience of the proposed route options is assessed, specifically as follows: when At that time, the overall resilience level of the proposed route was low, which did not meet the expected disaster mitigation route selection resilience requirements. The assessment result was that it was not recommended as the preferred route. when At that time, the overall resilience level of the proposed route was high, meeting the expected disaster mitigation route selection resilience requirements, and the assessment result recommended it as the preferred option.

[0013] The present invention also provides a readable storage medium storing computer program instructions, which, when executed by a processor, implement the above-described method for assessing the resilience of land transportation in mountainous areas during strong earthquakes.

[0014] The present invention also provides an electronic device, comprising: at least one processor, at least one memory, and computer program instructions stored in the memory, wherein the computer program instructions are executed by the processor to perform the above-described method for assessing the resilience of land transportation disaster reduction in mountainous areas during strong earthquakes.

[0015] The application of the technical solution of the present invention has the following beneficial effects: (1) This invention proposes a method for assessing the resilience of railway route selection in mountainous areas prone to strong earthquakes. It establishes a geographic information model by integrating geographic information parameters and constructs a seismic risk field model based on the peak ground acceleration (PGA) distribution. On this basis, the proposed route scheme is discretized at the segment level, and a segment-level resilience calculation model is established to obtain the resilience integral of each route segment. Based on this, the overall resilience of the route scheme is evaluated. This invention enables quantitative analysis and comprehensive evaluation of the resilience of the entire line during the route selection stage, filling a gap in the research on "pre-resilience assessment" of route schemes. It can provide strong technical support for the comparison and optimization of railway route selection in mountainous areas prone to strong earthquakes.

[0016] (2) When establishing the geographic information model, this invention comprehensively considers multiple influencing factors, including key geographic information parameters such as peak ground acceleration (PGA), distribution of rescue stations, distribution of existing road networks, topographic information, and geological information. This provides a complete data foundation for subsequent segment-level discretization and risk characterization of the route plan, which is conducive to more realistically and accurately reflecting the external environmental conditions of the route in the mountainous area affected by strong earthquakes, and realizing a systematic assessment of the entire route in the whole process of pre-disaster, during-disaster, and post-disaster.

[0017] (3) This invention constructs a segment-level discrete spatial model of the route scheme. Based on the geographic information model and the route scheme, it extracts and assigns complete attribute information to each segment, including structure type (bridge, tunnel, roadbed), bridge height, bridge span, tunnel burial depth, whether the tunnel crosses a fault, embankment filling height, road cut excavation depth, distance to the nearest rescue station, etc., forming a segment-level discrete spatial model of the route scheme and a corresponding segment-level attribute table, providing an efficient and unified data carrier for carrying out structural vulnerability analysis and toughness calculation, and is highly operable.

[0018] (4) This invention constructs a segment-level resilience calculation model, comprehensively considering seismic vulnerability, post-earthquake residual function, and recovery process, and converts the functional level of each segment over time into a segment-level resilience integral, ultimately achieving a quantitative assessment of the overall resilience of the proposed route during the route selection stage. This method can intuitively reflect the safety and recovery capability of different route schemes under strong earthquakes, and has good engineering applicability and promotion value.

[0019] In addition to the objectives, features, and advantages described above, the present invention has other objectives, features, and advantages. The invention will now be described in further detail with reference to the figures. Attached Figure Description

[0020] The accompanying drawings, which form part of this application, are used to provide a further understanding of the invention. The illustrative embodiments of the invention and their descriptions are used to explain the invention and do not constitute an undue limitation of the invention. In the drawings: Figure 1 This is a schematic diagram of the geographic information model in an embodiment of the present invention; Figure 2 This is a schematic diagram of the longitudinal profile geometry of the segment-level discrete space model of the candidate route scheme in this embodiment of the invention. Figure 3 This is a schematic diagram of the planar geometric shape of the segment-level discrete space model of the candidate route scheme in the embodiment of the present invention; Figure 4 This is a schematic diagram of a single recovery curve in an embodiment of the present invention. Detailed Implementation

[0021] To make the objectives, features, and advantages of the present invention more apparent and understandable, specific embodiments of the present invention will be described in detail below with reference to the accompanying drawings. Several embodiments of the present invention are shown in the drawings. However, the present invention can be implemented in many different forms and is not limited to the embodiments described herein.

[0022] In this invention, unless otherwise explicitly specified and limited, the terms "installation," "connection," "linking," and "fixing," etc., should be interpreted broadly. For example, they can refer to a fixed connection, a detachable connection, or an integral connection; they can refer to a mechanical connection or an electrical connection; they can refer to a direct connection or an indirect connection through an intermediate medium; and they can refer to the internal connection of two components. Those skilled in the art can understand the specific meaning of the above terms in this invention according to the specific circumstances. The terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of indicated technical features. Thus, a feature defined with "first" or "second" may explicitly or implicitly include one or more of that feature.

[0023] In this invention, unless otherwise expressly specified and limited, "above" or "below" a second feature can include direct contact between the first and second features, or contact between the first and second features through another feature between them. Furthermore, "above," "over," and "on top" of a second feature includes the first feature being directly above or diagonally above the second feature, or simply indicates that the first feature is at a higher horizontal level than the second feature. "Below," "below," and "under" of a second feature includes the first feature being directly below or diagonally below the second feature, or simply indicates that the first feature is at a lower horizontal level than the second feature. The terms "vertical," "horizontal," "left," "right," "above," "below," and similar expressions are for illustrative purposes only and do not indicate or imply that the device or element referred to must have a specific orientation, be constructed or operated in a specific orientation, and therefore should not be construed as limiting the invention.

[0024] The embodiments of the present invention will be described in detail below with reference to the accompanying drawings. However, the present invention can be implemented in many different ways, which may be limited to or covered by the present invention.

[0025] In one embodiment, a method for assessing the resilience of land transportation routes in mountainous areas during strong earthquakes includes the following steps: S1: Establish a geographic information model of the study area based on comprehensive geographic information within the study area; and construct an earthquake risk field model based on the geographic information model. In S1, the geographic information model includes topographic information, geological information, distribution of rescue stations, and existing road network; the seismic risk field model includes at least the spatial distribution of ground peak ground acceleration (PGA) within the study area, used to characterize the intensity level of ground motion at different locations.

[0026] See Figure 1 The study area was divided into × Divide the data into square cells and add the corresponding peak ground acceleration (PGA) value to each square cell to characterize the seismic intensity level at different locations, where: ; in, and These represent the length and width of the study area, respectively. The width of the square cell is represented; the center coordinates of each square represent its spatial location, and each square cell is assigned comprehensive geographic information parameters within the study area; S2: Discretize the candidate route scheme into multiple segment-level units, and combine the geographic information model and seismic risk field model of the study area to perform resilience calculations on each segment-level unit of the candidate route scheme to obtain the resilience integral corresponding to each segment-level unit of the candidate route scheme. See Figure 2 , The first longitudinal section A change of slope point, The first longitudinal section +1 slope change point, wherein S2 specifically includes: S2.1. Based on the existing three-dimensional route space model, construct a segment-level discrete space model for the candidate route schemes, specifically: S2.1.1 Obtain the spatial location of the candidate route scheme based on the existing three-dimensional route space model; S2.1.2 Discretize the candidate route schemes along the centerline of the route at preset intervals to obtain multiple segment-level units, and assign corresponding segment-level attributes to each segment-level unit. The segment-level attributes include: structure type, bridge height, bridge span, tunnel depth, whether the tunnel crosses a fault, embankment embankment height, cut excavation depth, and distance to the nearest rescue station. In this embodiment, the preset interval is 20m. Preferably, the structure types include bridges, tunnels, and roadbeds. The structure type can be determined by the difference between the route scheme and the ground elevation: according to the "Railway Line Design Code," when the embankment height is greater than or equal to 8m, a bridge is used instead of the roadbed; when the cut depth is greater than or equal to 15m, a tunnel is used instead of the roadbed. The embankment and cut heights can be obtained by subtracting the ground elevation from the route design elevation. See Figure 3 , For the plane of the first One intersection point, For the plane of the first +1 intersection point. For the calculation of the distance to the nearest rescue station, the geographic information model constructed in step S1 is used to construct a traffic network map using the existing road network distribution. The location of the rescue station and the center point of the route segment are projected onto the network. The shortest path search method is used to calculate the shortest path distance from the representative point of each route segment to the nearest rescue station along the road network. This distance is used as the distance to the nearest rescue station for each segment. S2.2 Based on the segment-level discrete spatial model of the candidate route scheme, and combined with the geographic information model and seismic risk field model of the study area, the resilience of each segment-level unit of the candidate route scheme is calculated to obtain the resilience integral corresponding to each segment-level unit of the candidate route scheme.

[0027] S2.2 includes: S2.2.1. Based on the seismic risk field model and the segment-level attributes of each segment-level unit, the structural vulnerability of each segment-level unit is calculated using a unified log-normal vulnerability curve paradigm; where: structural vulnerability includes the probability of occurrence of different structural damage states, specifically: no damage The probability of occurrence, minor injury The probability of occurrence, moderate injury The probability of occurrence and the severity of injury The probability of occurrence; In S2.2.1, a specific seismic motion level is given. No damage probability of occurrence The specific calculation formula is as follows: ; in: The condition is damaged. Minor injury The probability of exceeding; It is the standard normal distribution function; PAG is the ground motion intensity index in the earthquake risk field model; This indicates that the damage to the structure exceeds the level of minor damage. hour The middle threshold; Indicates confirmation The logarithmic standard deviation at that time; minor injury probability of occurrence The specific calculation formula is as follows: ; in: Moderate injury The probability of exceeding; This indicates that the damage to the structure exceeds the level of moderate damage. hour The middle threshold; Indicates confirmation The logarithmic standard deviation at that time; Moderate injury probability of occurrence The specific calculation formula is as follows: ; in: For severe injury The probability of exceeding; The structural damage condition exceeds severe damage. hour The middle threshold; To indicate certainty The logarithmic standard deviation at that time; Severe injury probability of occurrence The specific calculation formula is as follows: .

[0028] In this embodiment, the structure is in a damage state no lower than a certain level. The transcendence probability is expressed as: ; in, It is the standard normal distribution function; PAG is the seismic intensity index. Indicates that the structural damage condition exceeds hour The middle threshold; Indicates confirmation The logarithmic standard deviation at that time; In S2.2.1, for tunnel sections in the candidate route that do not cross fault zones and have a burial depth exceeding 400 meters, there is no damage. probability of occurrence The probability of occurrence of other damage states is 0, with a value of 1.

[0029] In this embodiment, the segment-level units of different structural segments are distinguished only by parameters. , The values ​​are distinguished as follows: Bridge Section: Based on characteristics such as bridge height, bridge length, and maximum span, a matching bridge vulnerability curve is selected to determine the corresponding damage state. , ; Tunnel Section: Based on factors such as tunnel depth and whether it crosses active fault zones, appropriate tunnel vulnerability curves are selected to determine the corresponding damage states. , ; Subgrade Section: Based on parameters such as embankment fill height and cut excavation depth, select the corresponding subgrade vulnerability curve to determine the corresponding damage state. , .

[0030] S2.2.2 Calculate the remaining function of each segment-level unit based on the structural vulnerability of each segment-level unit; S2.2.2 includes: ① Based on the structural vulnerability of each segment-level unit, the structural damage probability is mapped to the corresponding operational state, specifically: No damage Supports unlimited speed operation Minor injury Dual-line operation, slight speed limit Moderate injury Single-line operation, significantly reduced speed Severe injury Corresponding closed line For tunnel sections crossing fault zones, when their damage condition is moderate... If the conditions for safe operation are not met, the corresponding operational status will be directly determined as closed. ; ② Set the speed factor according to the operating status corresponding to each segment level element. and line usage factor Specifically: velocity factor Set post-earthquake speed limit With design speed The ratio, that is: In this embodiment, Take 350 km / h; The throughput capacity of the number of lines currently open in the operating state is set to the ratio of the throughput capacity when both lines are open. In this embodiment, the throughput capacity of single-line operation is 0.3, and the throughput capacity of dual-line operation is 1. This information can be obtained from the "Interim Provisions on Post-Earthquake Emergency Response for High-Speed ​​Railway Infrastructure". In this implementation, the values ​​of speed factors and line usage factors for each operating state are shown in Table 1. Table 1. Values ​​of speed factor and line usage factor for each operating state.

[0031] ③ Based on the velocity factor and line usage factor The remaining functional coefficient for each operating state is calculated using the following formula: ; in: Damaged state The residual function coefficient at that time, for , , , Any damage state; ④ Based on residual function coefficient The remaining functionality of each segment-level unit is calculated using the following formula: ; in: For remaining functions, For the first Segment-level unit damage status The probability of occurrence.

[0032] S2.2.3 Calculate the toughness integral of each segment-level element based on the residual function of each segment-level element, specifically: The equivalent recovery curves of the candidate route schemes are calculated based on the geographic information model, the seismic risk field model, and the residual function of each segment-level unit. Specifically, it includes: ① The rescue accessibility index is calculated based on the geographic information model and the earthquake risk field model. The specific calculation formula is as follows: ; in, This indicates the actual distance of the segment-level unit from the nearest rescue station; Indicates the base reachable distance; in this embodiment, ; ② Utilize the rescue accessibility index The typical repair time is corrected to obtain the corrected typical repair time. The specific formula is as follows: ; in: This represents the typical repair time, i.e., the time required to restore 80% of the functionality. In this embodiment, to avoid unreasonable amplification or reduction of the results due to extreme cases during correction, Take 2 at the time. Take 0.5 at a time; ③ An exponential recovery curve is used to characterize each damage state, and the recovery rate is calibrated using the corrected typical repair time to obtain the recovery curve for each damage state. The specific expression is as follows: ; in: Damaged state Down The function of time; It is a natural constant; For recovery rate; ④ Based on the principle of equivalent area, the recovery curves under the four injury states are weighted and fitted into a single equivalent recovery curve. In this embodiment, let the first... The segment is in a damaged state after the earthquake. The probability is The fitting process is as follows: ; in: The fitted recovery rate; See Figure 4 The specific expression for a single equivalent recovery curve is as follows: .

[0033] For equivalent recovery curve Within the control time interval Integrating over time, we obtain the first... The integral of the segment's toughness is calculated using the following formula: ; in: For the first The toughness integral of the segment-level unit, To control time, For a moment.

[0034] S3: Based on the resilience integrals corresponding to each segment-level unit of the candidate route scheme obtained in step S2, determine the overall resilience level index of the candidate route scheme; evaluate the disaster mitigation and route selection resilience of the candidate route scheme according to the minimum resilience integral and the overall resilience level index of the candidate route scheme to obtain the evaluation result.

[0035] S3 includes: S3.1. The equivalent downtime is defined based on a single equivalent recovery curve of segment-level units, and the specific formula is as follows: ; in, The downtime equivalent to functional degradation within the controlled time period; S3.2. Take the segment-level unit corresponding to the minimum resilience integral in the candidate route scheme as the most unfavorable segment, and then calculate the resilience integral of the most unfavorable segment. In this embodiment, the resilience level of the proposed route scheme is determined. Set a threshold for the first one. Set a threshold for the second one; Take 0.9, Take 0.8, specifically: That is, the proportion of equivalent downtime At that time, it maintained a high functional level and excellent resilience level for most of the time after the earthquake, entering S3.3; That is, the proportion of equivalent downtime At that time, there were certain speed limits and interruptions after the earthquake, but the overall repair was relatively fast, the resilience level was good, and it entered S3.3; That is, the proportion of equivalent downtime At that time, the post-earthquake operation and recovery capabilities were clearly insufficient, and the resilience level was low. The evaluation result of this candidate route was that it was not recommended as the preferred option. S3.3 Calculate the overall resilience level index of the candidate route scheme. The specific calculation formula is as follows: ; in: To satisfy the toughness integral of the entire line segment The cumulative length of the line segment, In this embodiment, a threshold is set for the third step. Take 0.8; The total length of the line; S3.4, Based on the overall resilience level index of the entire line The overall resilience of the proposed route options is assessed, specifically as follows: when At that time, the overall resilience level of the proposed route was low, which did not meet the expected disaster mitigation route selection resilience requirements. The assessment result was that it was not recommended as the preferred route. when At that time, the overall resilience level of the proposed route was high, meeting the expected disaster mitigation route selection resilience requirements, and the assessment result recommended it as the preferred option.

[0036] The threshold in this embodiment can be adjusted based on actual engineering experience, actual geographical conditions, and relevant specifications.

[0037] This embodiment also includes a readable storage medium storing computer program instructions, which, when executed by a processor, implement the above-described method for assessing the resilience of land transportation disaster reduction in mountainous areas during strong earthquakes.

[0038] It should be noted that the device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate, and the components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the modules can be selected to achieve the purpose of this embodiment according to actual needs. Furthermore, in the accompanying drawings of the device embodiments provided by this invention, the connection relationships between modules indicate that they have communication connections, which can be specifically implemented as one or more communication buses or signal lines. Those skilled in the art can understand and implement this without any creative effort.

[0039] This embodiment also includes an electronic device, comprising: at least one processor, at least one memory, and computer program instructions stored in the memory, wherein the computer program instructions are executed by the processor to perform the mining slope stability analysis method as described above.

[0040] For example, the computer program may be divided into one or more modules / units, which are stored in the memory and executed by the processor to complete the present invention. The one or more modules / units may be a series of computer program instruction segments capable of performing a specific function, which describe the execution process of the computer program in the electronic device.

[0041] The electronic device can be a mobile phone, desktop computer, laptop, handheld computer, cloud server, or other computing device. The electronic device may include, but is not limited to, processors and memory. For example, the electronic device may also include input / output devices, network access devices, buses, etc.

[0042] The processor can be a Central Processing Unit (CPU), or other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor can be a microprocessor or any conventional processor. The processor is the control center of the electronic device, connecting all parts of the electronic device via various interfaces and lines.

[0043] The memory can be used to store the computer program and / or modules. The processor implements the computer program by running or executing the computer program and / or modules stored in the memory, and by calling data stored in the memory. The memory may mainly include a program storage area and a data storage area. The program storage area may store the operating system, at least one application program required for a function (such as sound playback function, image playback function, etc.), etc.; the data storage area may store data created according to the use of the mobile phone (such as audio data, phonebook, etc.). In addition, the memory may include high-speed random access memory, and may also include non-volatile memory, such as hard disk, memory, plug-in hard disk, smart media card (SMC), secure digital card (SD) card, flash card, at least one disk storage device, flash memory device, or other volatile solid-state storage device.

[0044] Wherein, if the modules / units integrated in the electronic device are implemented as software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, all or part of the processes in the methods of the above embodiments of the present invention can also be implemented by a computer program instructing related hardware. The computer program can be stored in a readable storage medium, and when executed by a processor, it can implement the steps of the various method embodiments described above. The computer program includes computer program code, which can be in the form of source code, object code, executable files, or certain intermediate forms. The computer-readable medium can include: any entity or device capable of carrying the computer program code, recording media, USB flash drives, portable hard drives, magnetic disks, optical disks, computer memory, read-only memory (ROM), random access memory (RAM), electrical carrier signals, telecommunication signals, and software distribution media, etc.

[0045] The above description is merely a preferred embodiment of the present invention and is not intended to limit the invention. Various modifications and variations can be made to the present invention by those skilled in the art. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.

Claims

1. A method for assessing the resilience of land transportation routes in mountainous areas prone to strong earthquakes, characterized in that, Includes the following steps: S1: Establish a geographic information model of the study area based on comprehensive geographic information within the study area; and construct an earthquake risk field model based on the geographic information model. S2: Discretize the candidate route scheme into multiple segment-level units, and combine the geographic information model and seismic risk field model of the study area to perform resilience calculations on each segment-level unit of the candidate route scheme to obtain the resilience integral corresponding to each segment-level unit of the candidate route scheme. S2 specifically includes: S2.

1. Based on the existing three-dimensional route space model, construct a segment-level discrete space model for the candidate route schemes, specifically: S2.1.1 Obtain the spatial location of the candidate route scheme based on the existing three-dimensional route space model; S2.1.2 Discretize the candidate route schemes along the centerline of the route at preset intervals to obtain multiple segment-level units, and assign corresponding segment-level attributes to each segment-level unit, including: structure type, bridge height, bridge span, tunnel depth, whether the tunnel crosses a fault, embankment filling height, road cut excavation depth, and distance to the nearest rescue station. S2.2 Based on the segment-level discrete spatial model of the candidate route scheme, the resilience of each segment-level unit of the candidate route scheme is calculated by combining the geographic information model and the seismic risk field model of the study area to obtain the resilience integral corresponding to each segment-level unit of the candidate route scheme. S2.2 includes: S2.2.

1. Based on the seismic risk field model and the segment-level attributes of each segment-level unit, the structural vulnerability of each segment-level unit is calculated using a unified log-normal vulnerability curve paradigm; where: structural vulnerability includes the probability of occurrence of different structural damage states, specifically: no damage The probability of occurrence, minor injury The probability of occurrence, moderate injury The probability of occurrence and the severity of injury The probability of occurrence; S2.2.2 Calculate the remaining function of each segment-level unit based on the structural vulnerability of each segment-level unit; S2.2.2 includes: ① Based on the structural vulnerability of each segment-level unit, the structural damage probability is mapped to the corresponding operational state, specifically: No damage Supports unlimited speed operation Minor injury Dual-line operation, slight speed limit Moderate injury Single-line operation, significantly reduced speed Severe injury Corresponding closed line For tunnel sections crossing fault zones, when their damage condition is moderate... At that time, its corresponding operational status is directly determined to be closed. ; ② Set the speed factor according to the operating status corresponding to each segment level element. and line usage factor Specifically: velocity factor Set as the ratio of the post-earthquake restricted speed to the design speed; Set to the ratio of the throughput capacity of the number of lines currently open in the operating state to the throughput capacity when both lines are open; ③ Based on the velocity factor and line usage factor The remaining functional coefficient for each operating state is calculated using the following formula: ； in: Damaged state The residual function coefficient at that time for , , , Any damage state; ④ Based on residual function coefficient The remaining functionality of each segment-level unit is calculated using the following formula: ； in: For the remaining functions, For the first Segment-level unit damage status The probability of occurrence; S2.2.3 Calculate the toughness integral of each segment-level element based on the residual function of each segment-level element, specifically: The equivalent recovery curves of the candidate route schemes are calculated based on the geographic information model, the seismic risk field model, and the residual function of each segment-level unit. ; For equivalent recovery curve Within the control time interval Integrating over time, we obtain the first... The integral of the toughness of the segment is calculated using the following formula: ； in: For the first The toughness integral of the segment-level unit, To control time, For a specific moment; In step S2.2.3, the equivalent recovery curve of the candidate route scheme is calculated based on the geographic information model, the seismic risk field model, and the residual function of each segment-level unit. Specifically, it includes: ① The rescue accessibility index is calculated based on the geographic information model and the earthquake risk field model. The specific calculation formula is as follows: ； in, This indicates the actual distance of the segment-level unit from the nearest rescue station; Indicates the base reachable distance; ② Utilize the rescue accessibility index The typical repair time is corrected to obtain the corrected typical repair time. The specific formula is as follows: ； in: This is the typical repair time, which is the time required to restore 80% of the functionality; ③ An exponential recovery curve is used to characterize each damage state, and the recovery rate is calibrated using the corrected typical repair time to obtain the recovery curve for each damage state. The specific expression is as follows: ； in: Damaged state Down The function of time; It is a natural constant; Damaged state Recovery rate at the following levels; ④ Based on the principle of equivalent area, the recovery curves under the four injury states are weighted and fitted into a single equivalent recovery curve. Where: Let the first The segment is in a damaged state after the earthquake. The probability is The fitting process is as follows: ； in: The fitted recovery rate; Single equivalent recovery curve The specific expression is as follows: ； S3: Based on the resilience integrals corresponding to each segment-level unit of the candidate route scheme obtained in step S2, determine the overall resilience level index of the candidate route scheme; evaluate the disaster mitigation and route selection resilience of the candidate route scheme according to the minimum resilience integral and the overall resilience level index of the candidate route scheme to obtain the evaluation result.

2. The method for assessing the resilience of earthquake-prone mountainous land transportation routes for disaster reduction as described in claim 1, characterized in that, In S1, the geographic information model includes topographic information, geological information, distribution of rescue stations, and existing road network; the seismic risk field model includes at least the spatial distribution of ground peak ground acceleration (PGA) within the study area, used to characterize the intensity level of ground motion at different locations.

3. The method for assessing the resilience of earthquake-prone mountainous land transportation routes for disaster reduction as described in claim 1, characterized in that, In S2.2.1, a specific seismic motion level is given. No damage probability of occurrence The specific calculation formula is as follows: ； in: The condition is damaged. Minor injury The probability of exceeding; It is the standard normal distribution function; PAG is the ground motion intensity index in the earthquake risk field model; This indicates that the damage to the structure exceeds the level of minor damage. hour The median threshold; Indicates confirmation The logarithmic standard deviation at that time; minor injury probability of occurrence The specific calculation formula is as follows: ； in: Moderate injury The probability of exceeding; This indicates that the damage to the structure exceeds the level of moderate damage. hour The median threshold; Indicates confirmation The logarithmic standard deviation at that time; Moderate injury probability of occurrence The specific calculation formula is as follows: ； in: For severe injury The probability of exceeding; The structural damage condition exceeds severe damage. hour The median threshold; To indicate certainty The logarithmic standard deviation at that time; Severe injury probability of occurrence The specific calculation formula is as follows: ； In S2.2.1, for tunnel sections in the candidate route that do not cross fault zones and have a burial depth exceeding 400 meters, there is no damage. probability of occurrence The probability of occurrence of other damage states is 0, with a value of 1.

4. The method for assessing the resilience of earthquake-prone mountainous land transportation routes for disaster reduction as described in claim 1, characterized in that, S3 includes: S3.

1. The equivalent downtime is defined based on a single equivalent recovery curve of segment-level units, and the specific formula is as follows: ； in, The downtime equivalent to functional degradation within the controlled time period; S3.

2. Take the segment-level unit corresponding to the minimum resilience integral in the candidate route scheme as the most unfavorable segment, and then calculate the resilience integral of the most unfavorable segment. To determine the resilience level of the proposed route options, specifically: At that time, the toughness level was excellent, and it entered S3.3; At this point, the toughness level is good, and it enters S3.3; When the resilience level is low, the evaluation result of the proposed route is that it is not recommended as the preferred option. in: Set a threshold for the first one. Set a threshold for the second one; S3.3 Calculate the overall resilience level index of the candidate route scheme. The specific calculation formula is as follows: ； in: To satisfy the toughness integral of the entire line segment The cumulative length of the line segment, Set a threshold for the third; The total length of the line; S3.4, Based on the overall resilience level index of the entire line The overall resilience of the proposed route options is assessed, specifically as follows: when At that time, the overall resilience level of the proposed route was low, which did not meet the expected disaster mitigation route selection resilience requirements. The assessment result was that it was not recommended as the preferred route. when At that time, the overall resilience level of the proposed route was high, meeting the expected disaster mitigation route selection resilience requirements, and the assessment result recommended it as the preferred option.

5. A readable storage medium, characterized in that, It stores computer program instructions, which, when executed by a processor, implement the earthquake-prone mountain land transportation disaster reduction route selection resilience assessment method as described in any one of claims 1 to 4.

6. An electronic device, characterized in that, include: The method for assessing the resilience of mountainous land transportation disaster reduction routes as described in any one of claims 1 to 4 includes at least one processor, at least one memory, and computer program instructions stored in the memory, which are executed by the processor when the computer program instructions are executed.

Images