Fuzzy hierarchy evaluation method for metro disaster chain evolution time-space partition

Abstract

The invention relates to the technical field of disaster management, in particular to a metro disaster chain evolution space-time partition fuzzy hierarchy evaluation method, which comprises the following steps: dividing a target area into a plurality of space grid areas, recording disaster events of different time slices in each grid area, and extracting disaster space-time characteristics from the disaster events of different time slices; constructing a disaster chain network, and generating a disaster chain space-time diffusion map and risk partitions based on analysis of disaster space-time characteristics; identifying disaster gathering hotspots in the risk subareas, and predicting time-space information of key disasters through a time prediction model; introducing fuzzy logic, performing comprehensive evaluation on disaster risks according to time-space information of key disasters, and determining a high-risk area; and after the high-risk area is determined, performing secondary evaluation on the disaster risk in combination with the subway structure and the dredging capability.

Description

Technical Field

[0001] The present invention relates to the technical field of disaster management, and particularly to a fuzzy hierarchical evaluation method for spatio-temporal zoning of subway disaster chain evolution. Background Art

[0002] A disaster chain refers to a chain reaction phenomenon of a series of secondary disasters triggered by a primary disaster, forming a spatio-temporally correlated disaster sequence.

[0003] Through the analysis of the disaster chain, disasters can be scientifically managed and warned. Disasters can be controlled by dealing with the disaster-causing points and cutting off the propagation paths.

[0004] The subway system is a special system. When studying the subway disaster chain, it is necessary to consider both the underground system (tunnel structure, track, etc.) and the above-ground system (ground buildings, ground traffic, etc.). Most of the existing research focuses on the study of the underground space. For example, the Chinese patent with the publication number CN 118365209 A uses the analytic hierarchy process (AHP) as the research method. This method has the defects of strong subjectivity (easy to introduce human errors) and sensitivity to the 1-9 scale method (the subjectivity of scale selection will amplify the uncertainty of the results); at the same time, this method has great difficulty in dealing with fuzzy scenarios and cannot adapt to the analysis under complex scenarios; because the disaster accident data often has incompleteness or fuzziness, this incompleteness and fuzziness lead to a decrease in the credibility of the scoring value, resulting in the distortion of the calculation of the transfer probability for optimizing the disaster chain. Summary of the Invention

[0005] The present invention considers the above-ground space and the underground space of the subway system, combines the two to form a three-dimensional space, divides the three-dimensional space into multiple three-dimensional grids, analyzes the spatio-temporal characteristics of disasters in each three-dimensional grid, predicts key disasters using a time prediction model, and combines fuzzy logic and a risk value calculation model to evaluate the disaster risk multiple times.

[0006] The technical solution proposed by the present invention is: a fuzzy hierarchical evaluation method for spatio-temporal zoning of subway disaster chain evolution, the method comprising: Dividing the target area into multiple spatial grid areas, each grid area recording disaster events in different time slices, and extracting the spatio-temporal characteristics of disasters from the disaster events in different time slices; Constructing a disaster chain network, generating a spatio-temporal diffusion map and a risk zoning of the disaster chain based on the analysis of the spatio-temporal characteristics of disasters; Identifying the disaster aggregation hotspots within the risk zoning, and predicting the spatio-temporal information of key disasters through a time prediction model; Introducing fuzzy logic, comprehensively evaluating the disaster risk according to the spatio-temporal information of key disasters, and determining the high-risk areas; After determining the high-risk areas, a secondary evaluation of the disaster risk is carried out in combination with the subway structure and the evacuation capacity.

[0007] Preferably, the target area is divided into multiple spatial grid areas, and each grid area records disaster events in different time slices. The disaster spatio-temporal characteristics are extracted from the disaster events in different time slices, including: Obtain the underground data and ground data of the target area; Unify the coordinate system and time stamp, align the underground data and ground data in space and time, and construct an underground-ground three-dimensional space; Divide the underground-ground space into multiple three-dimensional grids; specifically: segment the target area horizontally and layer it vertically to form multiple three-dimensional grids; Obtain the historical disaster data within each three-dimensional grid, and extract the type of disaster event, the number of occurrences of the disaster event, and the occurrence time point of the disaster event from the historical disaster data; Calculate the disaster density within the three-dimensional grid area , where represents the number of disasters occurring within the time period of the th grid area, represents the time span, represents the number of three-dimensional grids.

[0008] Preferably, the construction of the disaster chain network includes: Obtain the accident data within the target area from the database, and extract the disaster-causing factors from the accident data; the accident data includes leakage, settlement, and fire; the disaster-causing factors include rising groundwater level and structural cracks; Determine the causal relationship between the disaster-causing factors through the expert knowledge base, and construct the disaster propagation path, that is, the disaster chain; Taking the disaster-causing factors as nodes and the causal relationship between the disaster-causing factors as edges, construct a disaster chain topology graph, that is, the disaster chain network.

[0009] Preferably, the construction of the disaster chain network further includes: Model the disaster chain network, including: Construct a mechanical-hydrological coupling model based on the modified Mohr-Coulomb constitutive equation and Darcy's law to simulate the soil deformation and groundwater seepage during disasters; The mechanical-hydrological model is ; where represents the effective cohesive force, represents the effective internal friction angle, represents the permeability coefficient, represents the pore water pressure; represents the water head height; represents the stress tensor; Construct a thermal-structural coupling model based on the heat conduction equation and the thermal stress equation to simulate the heat conduction during a disaster; The thermal-structural coupling model is ; where represents heat; represents the thermal diffusivity, represents the heat gradient; represents the coefficient of thermal expansion, represents the elastic modulus; represents the thermal stress; Establish a disaster propagation model among the nodes of the disaster chain through a Bayesian network or a dynamic system equation, specifically including: Calculate the probability of disaster propagation between nodes based on the Bayesian network ; Based on the differential equation, that is, the dynamic system equation, describe the evolution process of the disaster intensity over time: , where represents the disaster intensity, represent the diffusion rate and the attenuation coefficient respectively.

[0010] Preferably, based on the analysis of the spatio-temporal characteristics of the disaster, generate a spatio-temporal diffusion map and a risk partition of the disaster chain, including: Let the three-dimensional grid of the target area be , where represents the indices of the X, Y, and Z axes of the th three-dimensional grid; The time slice of each three-dimensional grid is ; Obtain the monitoring data of each three-dimensional grid and perform normalization processing; the monitoring data includes the settlement amount and the leakage rate; Extract the spatio-temporal characteristics of the disaster, including: Obtain the disaster intensity , where represents the weight coefficient; represent the probability of disaster occurrence, the disaster diffusion speed, and the disaster duration after normalization respectively; Obtain the cumulative disaster intensity ; where represents the disaster attenuation coefficient; If the disaster intensity of the adjacent three-dimensional grid is , then trigger diffusion: ; where represents the diffusion threshold, represents the diffusion efficiency, Indicates the number of interval grids between two three-dimensional grids; Map the cumulative intensity of each three-dimensional grid to a color gradient, and display the risk level of the cumulative disaster intensity through colors; Use time as the X-axis and disaster intensity as the Y-axis to show the process of the spread of disaster intensity; Perform risk zoning through the clustering algorithm K-means++, including: Obtain the risk values of all three-dimensional grids, specifically: Take the maximum cumulative disaster intensity of each three-dimensional grid as its risk value, that is, the risk value ; Indicates the time period when the disaster occurs Within a preset time span ; Determine the optimal number of clusters through the elbow method ; Establish a grading function ; Among them, Indicates The central value of the cluster; Indicates the three-dimensional space of the target area; Set risk level labels, and the real-time level labels include high risk, medium risk, and low risk; When , divide the three-dimensional grid Into high-risk zones; when , divide the three-dimensional grid Into medium-risk zones; when , divide the three-dimensional grid Into low-risk zones; Perform dilation or erosion operations on the clustering results to smooth the boundaries of each risk zone.

[0011] Preferably, the identification of disaster aggregation hotspots within the risk zone includes: Identify disaster aggregation hotspots within the risk zone through the local spatial autocorrelation analysis algorithm, including: Calculate the local Moran's I index within the risk zone ; Among them, Indicates the average disaster intensity value; Indicates the variance of the global disaster intensity value; Indicates the spatial adjacency weight; Indicates the number of three-dimensional grids, Indicates the Disaster intensity value of the th adjacent three-dimensional grid; Calculate the p-value through permutation test or Z-value calculation. If p < 0.05, it is considered that the spatial autocorrelation is significant, and the significantly aggregated areas of disasters within the risk area are screened out accordingly. Calculate the kernel density of the significantly aggregated areas. ; Where: represents the number of disaster events, represents the Gaussian kernel function; represents the bandwidth; represents the coordinates of the center point of the area, represents the coordinates of the point to be evaluated; If the kernel density value is greater than the density threshold, it is determined that the aggregated area is a hot spot, otherwise it is a cold spot; Using the time prediction model to predict the spatio-temporal information of key disasters, including: Construct a spatio-temporal prediction model for key disasters, including: Construct a spatio-temporal autoregressive moving average model STARMA to predict the disaster intensity at a future time point, that is: ; Where, respectively represent the disaster intensity, spatial weight matrix, autoregressive coefficient, and moving average coefficient at position in the aggregated area at time , represents the time lag order; and respectively represent the random error term and the moving error term; Construct a long short-term memory network LSTM and train the long short-term memory network with historical disaster data; use the trained LSTM to predict the type of disaster that will occur at the hot spot at a future time point, that is, the type of key disaster; Use the pre-trained spatio-temporal graph convolutional network ST-GCN to predict the disaster intensity that will occur at the hot spot within a future time period, that is, the intensity of key disasters; Perform Kriging interpolation on the prediction results to generate a continuous risk surface to reflect the risk distribution range; Overlay the prediction results with the hot spot area to identify potential future hot spots.

[0012] Preferably, introduce fuzzy logic to comprehensively evaluate the disaster risk based on the spatio-temporal information of key disasters and determine the high-risk areas, including: Obtain the spatio-temporal information of key disasters within the hot spot area and construct key disaster characteristic variables; the spatio-temporal information of key disasters includes disaster duration, diffusion speed, hot spot kernel density, and tunnel deformation amount; Construct "IF-THEN" fuzzy rules; Based on fuzzy rules, using the Mamdani fuzzy inference algorithm, after inputting the key disaster feature variables, a fuzzy conclusion is output; Through the MAX operation, all the output fuzzy conclusions are aggregated to form a fuzzy set; Through the center of gravity method COG for defuzzification, the center of gravity position of the fuzzy set is calculated; Select the conclusion corresponding to the value of the center of gravity position as the output conclusion; Calculate the comprehensive risk value , where, represents the value of the center of gravity position output by the th fuzzy rule; represents the rule weight; represents the number of fuzzy rules; Set the high-risk threshold , if then the output hot spot is a high-risk area.

[0013] Preferably, after determining the high-risk area, combined with the subway structure and evacuation capacity, a secondary evaluation of the disaster risk is carried out, including: Obtain a list of high-risk areas, including the longitude and latitude of the high-risk area, the type of disaster, and the comprehensive risk value; Obtain subway structure data, including the length of the subway tunnel, the depth of the subway tunnel, the support material, and the support structure; Obtain subway evacuation capacity data, including the location and number of subway entrances and exits, the location and number of emergency exits, the distance from each emergency exit to the high-risk area, and the maximum passenger flow capacity of the subway; Normalize the subway structure data and subway evacuation capacity data; Obtain the structural safety factor , where, respectively represent the structural strength, the structural cross-sectional area, the safety redundancy factor, and the disaster acting force; Obtain the evacuation capacity coefficient ; The actual evacuation time , where represents the maximum number of people to be evacuated in the area; respectively represent the number of people that can pass through per second per unit width of the emergency exit, the width of a single emergency exit, and the number of available emergency exits; represents the evacuation attenuation coefficient; represents the minimum evacuation time, ; Modify the risk formula to obtain the secondary evaluation score ; where, represents the structure weight and the evacuation weight.

[0014] The present invention also provides an electronic device, including a processor, a memory connected to the processor, and a communication module. The electronic device is configured to execute the fuzzy hierarchical evaluation method for spatio-temporal partitioning of subway disaster chain evolution.

[0015] The present invention also provides a computer-readable storage medium storing a computer program, which is executed by a processor to implement the fuzzy hierarchical evaluation method for spatio-temporal partitioning of subway disaster chain evolution.

[0016] Advantages of the present invention: 1. For the subway scenario which includes both ground and underground structures, the present invention considers physical coupling (mechanical and hydrological interactions) and system coupling (causal relationships between disaster events), constructs a mechanical-hydrological and thermal-structural model, and determines the nodes (disasters) and edges (causal relationships) of the disaster chain network through the output of the model. By analyzing the disaster chain network, key nodes can be found and targeted prevention and control can be carried out for the key nodes.

[0017] 2. The present invention divides the subway space into multiple three-dimensional grids, identifies high-risk areas by analyzing the spatio-temporal characteristics of disasters in each three-dimensional grid; further analyzes the high-risk areas to determine disaster aggregation hotspots, and introduces fuzzy logic to preliminarily evaluate the disaster risks of the disaster aggregation points, which can effectively integrate the uncertainties in spatio-temporal information and generate a more practical high-risk area division. Finally, combined with the different structural characteristics and evacuation capabilities of each subway tunnel, the disaster risks are evaluated again. Brief Description of the Drawings

[0018] Figure 1 It is a flowchart of the fuzzy hierarchical evaluation method for spatio-temporal partitioning of subway disaster chain evolution of the present invention. Detailed Embodiments

[0019] The following description is used to disclose the present invention so that those skilled in the art can implement the present invention. The preferred embodiments in the following description are only examples, and those skilled in the art can think of other obvious deformations. The basic principles defined in the following description can be applied to other implementation schemes, deformation schemes, improvement schemes, equivalent schemes, and other technical schemes that do not depart from the spirit and scope of the present invention.

[0020] It can be understood that the term "one" should be understood as "at least one" or "one or more". That is, in one embodiment, the number of an element can be one, while in other embodiments, the number of the element can be multiple. The term "one" cannot be understood as a limitation on the number.

[0021] Embodiment 1

[0022] Reference Figure 1 , the technical solution provided by the present invention is: a fuzzy hierarchical evaluation method for the spatio-temporal partitioning of the evolution of subway disaster chains, including the following steps: Step 1: Divide the target area into multiple spatial grid areas. Each grid area records disaster events in different time slices, and extract the disaster spatio-temporal characteristics from the disaster events in different time slices. This step includes the following sub-steps: Obtain the underground data and ground data of the target area; Unify the coordinate system and time stamp, align the underground data and ground data spatio-temporally, and construct an underground-ground three-dimensional space; Divide the underground-ground space into multiple three-dimensional grids; specifically: segment the target area horizontally and layer it vertically to form multiple three-dimensional grids; Obtain the historical disaster data within each three-dimensional grid, and extract the type of disaster event, the number of occurrences of the disaster event, and the occurrence time point of the disaster event from the historical disaster data; Calculate the disaster density within the three-dimensional grid area , where represents the th grid area The number of disasters occurring within the time period represents the time span, represents the number of three-dimensional grids.

[0023] Step 2: Construct a disaster chain network, and generate a spatio-temporal diffusion map and risk partition of the disaster chain based on the analysis of the disaster spatio-temporal characteristics.

[0024] The steps for constructing the disaster chain network are as follows: Obtain the accident data within the target area from the database, and extract the disaster-causing factors from the accident data; the accident data includes leakage, settlement, and fire; the disaster-causing factors include rising groundwater level and structural cracks; Determine the causal relationship between the disaster-causing factors through an expert knowledge base, and construct a disaster propagation path, that is, a disaster chain; Use the disaster-causing factors as nodes and the causal relationship between the disaster-causing factors as edges to construct a disaster chain topology diagram, that is, a disaster chain network.

[0025] Determine the high-risk paths through the analysis of the disaster chain topology diagram (such as the PageRank algorithm). For example: "heavy rain - rising groundwater level - soil softening - ground collapse" may be a key path. Identify the key nodes in the network (such as "soil softening"), and strengthening them can significantly reduce the global risk.

[0026] Among them, the steps for generating a spatio-temporal diffusion map and risk partition of the disaster chain based on the analysis of the disaster spatio-temporal characteristics are as follows: Let the three-dimensional grid of the target area be , where represents the indices of the X, Y, and Z axes of the th three-dimensional grid; The time slice of each three-dimensional grid is ; Obtain the monitoring data of each three-dimensional grid and perform normalization processing; the monitoring data includes settlement amount and leakage rate; Extract the spatio-temporal characteristics of disasters, including: Obtain the disaster intensity , where represents the weight coefficient; respectively represent the probability of disaster occurrence, the disaster diffusion speed, and the disaster duration after normalization processing; Obtain the cumulative disaster intensity ; where represents the disaster attenuation coefficient; If the disaster intensity of adjacent three-dimensional grids is , then trigger diffusion: ; where represents the diffusion threshold, represents the diffusion efficiency, represents the number of intermediate grids between two three-dimensional grids; Map the cumulative intensity of each three-dimensional grid to a color gradient, and display the risk level of the cumulative disaster intensity through colors; Use time as the X-axis and disaster intensity as the Y-axis to show the process of disaster intensity diffusion; Perform risk zoning through the K-means++ clustering algorithm, including: Obtain the risk values of all three-dimensional grids, specifically: Use the maximum cumulative disaster intensity of each three-dimensional grid as its risk value, that is, the risk value ; represents the time period of disaster occurrence within the preset time span ; Determine the optimal number of clusters through the elbow method ; Establish a hierarchical classification function ; where represents the central value of the cluster; Set risk level labels, and the real-time level labels include high risk, medium risk, and low risk; When , the three-dimensional grid Divided into high-risk zones; when the three-dimensional grid is divided into medium-risk zones; when the three-dimensional grid is divided into low-risk zones; Perform dilation or erosion operations on the clustering results to smooth the boundaries of each risk zone.

[0027] Step 3: Identify the disaster concentration hotspots within the risk zones, and predict the spatio-temporal information of key disasters through a time prediction model.

[0028] The steps for identifying the disaster concentration hotspots within the risk zones are as follows: Identify the disaster concentration hotspots within the risk zones through the local spatial autocorrelation analysis algorithm, including: Calculate the local Moran's I index within the risk zone ; where represents the average disaster intensity value; represents the variance of the global disaster intensity value; represents the spatial adjacency weight; represents the number of three-dimensional grids, represents the th adjacent three-dimensional grid's disaster intensity value; Calculate the p-value through permutation test or Z-value calculation. If p < 0.05, it is considered that the spatial autocorrelation is significant, and thus the significantly clustered areas within the risk zone are screened out; Calculate the kernel density of the significantly clustered areas; ; Where: represents the number of disaster events, represents the Gaussian kernel function; represents the bandwidth; represents the coordinates of the regional center point, represents the coordinates of the point to be evaluated; If the kernel density value is greater than the density threshold, then judge that the clustered area is a hotspot, otherwise it is a cold spot; Among them, predicting the spatio-temporal information of key disasters through a time prediction model includes the following steps: Construct a spatio-temporal prediction model for key disasters, specifically: Construct a spatio-temporal autoregressive moving average model STARMA to predict the disaster intensity at the next time point, that is: ; where respectively represent the disaster intensity, spatial weight matrix, autoregressive coefficient, and moving average coefficient at the position in the clustered area at time ​ Denote the order of time lag; and respectively denote the random error term and the moving error term; Construct a long short-term memory network (LSTM) and train the LSTM with historical disaster data; Use the trained LSTM to predict the type of disaster that will occur at the hot spot at the next time point, that is, the type of key disaster; Use a pre-trained spatio-temporal graph convolutional network (ST-GCN) to predict the intensity of the disaster that will occur at the hot spot within the next time period, that is, the intensity of the key disaster; Perform Kriging interpolation on the prediction results to generate a continuous risk surface to reflect the risk distribution range; Overlay the prediction results with the hot spot area to identify potential future hot spots.

[0029] Step 4: Introduce fuzzy logic to comprehensively evaluate the disaster risk based on the spatio-temporal information of the key disaster and determine the high-risk areas. It includes the following sub-steps: Obtain the spatio-temporal information of the key disaster in the hot spot area and construct the characteristic variables of the key disaster; The spatio-temporal information of the key disaster includes the disaster duration, diffusion speed, hot spot kernel density, and tunnel deformation; Construct "IF-THEN" fuzzy rules; Based on the fuzzy rules, use the Mamdani fuzzy inference algorithm to output a fuzzy conclusion after inputting the characteristic variables of the key disaster; Through the MAX operation, aggregate all the output fuzzy conclusions to form a fuzzy set; Through the center of gravity method (COG) for defuzzification, calculate the center of gravity position of the fuzzy set; Select the conclusion corresponding to the value of the center of gravity position as the output conclusion; Calculate the comprehensive risk value , where denotes the value of the center of gravity position output by the th fuzzy rule; denotes the rule weight; denotes the number of fuzzy rules; Set the high-risk threshold , if then output that the hot spot is a high-risk area.

[0030] For example, when there is a leakage in the subway tunnel, the three-dimensional grid A(50, 60, 2) of the subway tunnel is the hot spot area. Use the triangular membership function to determine the level of membership. Define the boundary of the fuzzy set with the vertex of the triangle (which can be obtained based on historical data analysis).

[0031] The calculated disaster intensity is 75, with a membership degree of "high" and a weight of 0.8; the diffusion speed is 2.5 m / h, with a membership degree of "high" and a weight of 0.7.

[0032] Fuzzy rule 1: IF hotspot area AND disaster intensity is high THEN high risk (rule weight is 0.9); Fuzzy rule 2: IF diffusion speed is fast THEN high risk (rule weight is 0.7); After defuzzification, calculated by the centroid method , if the high-risk threshold is 0.75 and 0.82 > 0.75, then it is determined that this hotspot is a high-risk area.

[0033] Step 5: After determining the high-risk area, combine the subway structure and evacuation capacity to conduct a secondary evaluation of the disaster risk. It includes the following sub-steps: Obtain the list of high-risk areas, including the longitude and latitude, disaster type, and comprehensive risk value of the high-risk areas; Obtain the subway structure data, including the subway tunnel length, subway tunnel depth, support material, and support structure; Obtain the subway evacuation capacity data, including the location and number of subway entrances and exits, the location and number of emergency exits, the distance from each emergency exit to the high-risk area, and the maximum passenger flow capacity of the subway; Normalize the subway structure data and subway evacuation capacity data; Obtain the structural safety factor , where respectively represent the structural strength, structural cross-sectional area, safety redundancy factor (1.2 - 1.5), and disaster acting force; Obtain the evacuation capacity coefficient ; Actual evacuation time , where represents the maximum number of people to be evacuated in the area; respectively represent the number of people that can pass through per second per unit width of the emergency exit, the width of a single emergency exit, and the number of available emergency exits; represents the evacuation attenuation coefficient; represents the minimum evacuation time, ; Modify the risk formula to obtain the secondary evaluation score ; where represents the structure weight and evacuation weight.

[0034] For example, when there is a leakage in the subway tunnel, calculated by the centroid method .

[0035] The structural parameters are as follows: the thickness of the support structure is 0.4 m, the strength of the concrete support structure is 50 MPa, and the water pressure of the disaster is 200 kN. The drainage parameters are , , , , ; Then , minutes; ; , and this area is of low risk.

[0036] It can be seen that the original risk is relatively high, but the structural safety is strong. Even when the drainage capacity is insufficient (much less than 1), the comprehensive risk assessment after secondary assessment is of low risk.

[0037] Embodiment 2: When constructing a disaster chain and a disaster chain network for the subway system, it is necessary to consider the spatio-temporal correlation of ground and underground space disasters. When studying its spatio-temporal correlation, it is necessary to quantify the abstract causal relationship. For this reason, on the basis of Embodiment 1, we propose the following technical solutions: Model the disaster chain network, specifically: Construct a mechanical-hydrological coupling model based on the modified Mohr-Coulomb constitutive equation and Darcy's law to simulate the soil deformation and groundwater seepage during disasters; The mechanical-hydrological model is ; where represents the effective cohesive force, represents the effective internal friction angle, represents the permeability coefficient, represents the pore water pressure; represents the water head height; represents the stress tensor; Construct a thermal-structural coupling model based on the heat conduction equation and the thermal stress equation to simulate the heat conduction during disasters; The thermal-structural coupling model is ; where represents the heat, represents the thermal diffusivity, represents the heat gradient; represents the coefficient of thermal expansion, represents the elastic modulus; represents the thermal stress; Establish a disaster propagation model between the nodes of the disaster chain through a Bayesian network or a dynamic system equation, specifically including: Calculate the probability of disaster propagation between nodes based on the Bayesian network ; Based on differential equations, i.e., dynamic system equations, to describe the evolution process of disaster intensity over time: , where represents the disaster intensity, respectively represent the diffusion rate and the attenuation coefficient.

[0038] For example, the input of the mechanics - hydrology model: the rainfall intensity is 120 mm / h, and the rainfall duration is 6 hours, resulting in a 2 - m rise in the groundwater level.

[0039] The model output: the pore water pressure is 60 kPa, and the predicted settlement is 12 mm.

[0040] The disaster chain is mapped to the rise of the groundwater level - soil softening - tunnel settlement.

[0041] The input of the thermo - structure model: the fire temperature is 800 °C, and the duration is 30 minutes, resulting in a 40% decrease in the concrete strength. The model output shows that the thermal stress is 50 MPa, exceeding the tensile strength of the lining by 45 MPa. The disaster chain is high fire temperature - lining cracking - tunnel collapse.

[0042] The mechanics - hydrology and thermo - structure models provide a physical mechanism support for the disaster chain, transforming the abstract causal relationship into quantitative parameters (such as probability, intensity). The model output directly defines the nodes and edges of the disaster chain network. Through the topological analysis of the disaster chain network, the mapping from the local physical process to the global risk network is realized.

[0043] Meanwhile, it improves the accuracy of disaster prediction, assists in formulating targeted prevention and control strategies, and takes remedial measures (such as strengthening the drainage system, adding fire isolation layers) for key nodes (for example, soil softening) to reduce the operation risk of the subway.

[0044] The present invention also provides an electronic device, including a processor, a memory connected to the processor, and a communication module. The electronic device is used to execute the fuzzy hierarchical evaluation method for the spatio - temporal partitioning of the subway disaster chain evolution.

[0045] The present invention also provides a computer - readable storage medium. The computer - readable storage medium stores a computer program, and the computer program is executed by a processor to implement the fuzzy hierarchical evaluation method for the spatio - temporal partitioning of the subway disaster chain evolution.

[0046] Embodiments disclosed in the present invention, the processes described above with reference to the flowcharts can be implemented as computer software programs. Embodiments disclosed in the present invention include a computer program product, which includes a computer program carried on a computer-readable medium, and the computer program contains program codes for executing the methods shown in the flowcharts. In such an embodiment, the computer program can be downloaded and installed from a network through a communication part, and / or installed from a removable medium. When the computer program is executed by a central processing unit (CPU), the above-mentioned functions defined in the method of the present invention are executed. It should be noted that the computer-readable medium in the present invention can be a computer-readable signal medium or a computer-readable storage medium or any combination of the two. A computer-readable storage medium can, for example but not limited to, be an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination of the above. More specific examples of a computer-readable storage medium can include but are not limited to: an electrical connection having one or more wire segments, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the above. In the present invention, a computer-readable storage medium can be any tangible medium that contains or stores a program, and this program can be used by or combined with an instruction execution system, apparatus, or device. In the present invention, a computer-readable signal medium can include a data signal propagated in a baseband or as part of a carrier wave, in which computer-readable program codes are carried. Such a propagated data signal can take various forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination of the above. A computer-readable signal medium can also be any computer-readable medium other than a computer-readable storage medium, and this computer-readable medium can send, propagate, or transmit a program for use by or combined with an instruction execution system, apparatus, or device. The program codes contained on the computer-readable medium can be transmitted by any appropriate medium, including but not limited to: a wireless segment, a wire segment, an optical cable, RF, etc., or any suitable combination of the above.

[0047] The flowcharts and block diagrams in the accompanying drawings illustrate the possible architectures, functions, and operations of systems, methods, and computer program products according to various embodiments of the present invention. In this regard, each block in the flowchart or block diagram may represent a module, a segment of a program, or a part of code that contains one or more executable instructions for implementing a specified logical function. It should also be noted that in some alternative implementations, the functions marked in the blocks may occur in a different order than that marked in the accompanying drawings. For example, two consecutive blocks shown may actually be executed substantially in parallel, and they may sometimes be executed in the reverse order, depending on the functions involved. It should also be noted that each block in the block diagram and / or flowchart, and combinations of blocks in the block diagram and / or flowchart, may be implemented by a dedicated hardware-based system that performs the specified functions or operations, or may be implemented by a combination of dedicated hardware and computer instructions.

[0048] Those skilled in the art should understand that the embodiments of the present invention described above and shown in the accompanying drawings are only examples and do not limit the present invention. The objectives of the present invention have been fully and effectively achieved. The functions and structural principles of the present invention have been demonstrated and explained in the embodiments. Without departing from the said principles, any changes or modifications can be made to the embodiments of the present invention.

Claims

1. A fuzzy hierarchical evaluation method for the spatio-temporal partitioning of the evolution of subway disaster chains, characterized in that, The method includes: Dividing the target area into multiple spatial grid areas, each grid area recording disaster events in different time slices, and extracting disaster spatio-temporal features from the disaster events in different time slices; Constructing a disaster chain network, generating a disaster chain spatio-temporal diffusion map and risk zoning based on the analysis of disaster spatio-temporal features; Identifying disaster aggregation hotspots within the risk zoning, and predicting the spatio-temporal information of key disasters through a time prediction model; Introducing fuzzy logic to comprehensively evaluate disaster risks according to the spatio-temporal information of key disasters and determining high-risk areas; After determining the high-risk areas, combining the subway structure and evacuation capacity to conduct a secondary evaluation of disaster risks.

2. The fuzzy hierarchical evaluation method for spatio-temporal zoning of subway disaster chain evolution according to claim 1, characterized in that The dividing the target area into multiple spatial grid areas, each grid area recording disaster events in different time slices, and extracting disaster spatio-temporal features from the disaster events in different time slices includes: Obtaining underground data and ground data of the target area; Unifying the coordinate system and time stamp, aligning the underground data and ground data spatio-temporally, and constructing an underground-ground three-dimensional space; Dividing the underground-ground space into multiple three-dimensional grids; specifically: segmenting the target area horizontally and layering it vertically to form multiple three-dimensional grids; Obtaining historical disaster data within each three-dimensional grid, and extracting the type of disaster events, the occurrence times of disaster events, and the occurrence time points of disaster events from the historical disaster data; Calculate the disaster density within the three-dimensional grid area , where represents the number of disasters that occurred during the time period of the -th grid area, represents the time span, represents the number of three-dimensional grids.

3. The fuzzy hierarchical evaluation method for spatio-temporal zoning of subway disaster chain evolution according to claim 2, characterized in that, The constructing the disaster chain network includes: Obtaining accident data within the target area from the database and extracting disaster-causing factors from the accident data; the accident data includes leakage, settlement, and fire; the disaster-causing factors include rising groundwater level and structural cracks; Determining the causal relationship between disaster-causing factors through an expert knowledge base and constructing a disaster propagation path, i.e., a disaster chain; Constructing a disaster chain topology graph, i.e., a disaster chain network, with disaster-causing factors as nodes and the causal relationship between disaster-causing factors as edges.

4. The fuzzy hierarchical evaluation method for spatio-temporal zoning of subway disaster chain evolution according to claim 3, characterized in that, The constructing the disaster chain network also includes: Modeling the disaster chain network, including: Constructing a mechanical-hydrological coupling model based on the modified Mohr-Coulomb constitutive equation and Darcy's law to simulate soil deformation and groundwater seepage during disasters; The mechanical-hydrological model is ; where represents the effective cohesive force, represents the effective internal friction angle, represents the permeability coefficient, represents the pore water pressure; represents the water head height; represents the stress tensor; Constructing a thermal-structural coupling model based on the heat conduction equation and thermal stress equation to simulate heat conduction during disasters; The thermo-structural coupling model is ; where represents heat represents the thermal diffusivity represents the heat gradient represents the coefficient of thermal expansion represents the elastic modulus represents the thermal stress Establishing a disaster propagation model between each node of the disaster chain through a Bayesian network or dynamic system equation, specifically including: Calculating the probability of disaster propagation between nodes based on the Bayesian network ； Based on differential equations, that is, dynamic system equations, to describe the evolution process of disaster intensity over time: , where represents the disaster intensity, represent the diffusion rate and the attenuation coefficient respectively.

5. The fuzzy hierarchical evaluation method for spatio-temporal zoning of subway disaster chain evolution according to claim 4, characterized in that The generating a disaster chain spatio-temporal diffusion map and risk zoning based on the analysis of disaster spatio-temporal features includes: Let the three-dimensional grid of the target area be , where represents the indices of the X, Y, and Z axes of the th three-dimensional grid; The time slice of each three-dimensional grid is ; Obtaining the monitoring data of each three-dimensional grid and performing normalization processing; the monitoring data includes settlement amount and leakage rate; Extracting disaster spatio-temporal features, including: Obtain the disaster intensity , where represents the weight coefficient; respectively represent the disaster occurrence probability, disaster diffusion speed, and disaster duration after normalization processing; Obtain the cumulative intensity of disasters ; where represents the disaster attenuation coefficient; If the disaster intensity of adjacent three-dimensional grids is , then diffusion is triggered: ; where represents the diffusion threshold, represents the diffusion efficiency, represents the number of spaced grids between two three-dimensional grids; Mapping the cumulative intensity of each three-dimensional grid into a color gradient and displaying the risk degree of the disaster cumulative intensity through colors; Taking time as the X-axis and disaster intensity as the Y-axis to display the diffusion process of disaster intensity; Performing risk zoning through the clustering algorithm K-means++, including: Obtaining the risk values of all three-dimensional grids, specifically: Take the maximum cumulative intensity of disasters for each three-dimensional grid as its risk value, that is, the risk value ; Indicates the time period when the disaster occurs Within the preset time span ; Determine the optimal number of clusters by the elbow method ; Establish a grading function ; among which, represents the central value of the cluster; represents the three-dimensional space of the target area; Setting risk level labels, and the real-time level labels include high risk, medium risk, and low risk; When the three-dimensional grid is divided into high-risk zones; when the three-dimensional grid is divided into medium-risk zones; when the three-dimensional grid is divided into low-risk zones; Perform dilation or erosion operations on the clustering results to smooth the boundaries of each risk zone.

6. The fuzzy hierarchical evaluation method for the spatio-temporal zoning of the subway disaster chain evolution according to claim 5, wherein Identifying disaster concentration hotspots within the identified risk zones includes: Identifying disaster concentration hotspots within the risk area through the local spatial autocorrelation analysis algorithm, including: Calculate the local Moran's I index within the risk area ; where represents the average disaster intensity value; represents the variance of the global disaster intensity value; represents the spatial adjacency weight; represents the number of three-dimensional grids, represents the disaster intensity value of the th adjacent three-dimensional grid; Calculating the p-value through permutation test or Z-value. If p < 0.05, it is considered that the spatial autocorrelation is significant, and thus the areas with significant disaster concentration within the risk area are screened out; Calculating the kernel density of the significantly concentrated areas; ; where: represents the number of disaster events, represents the Gaussian kernel function; represents the bandwidth; represents the coordinates of the center point of the area, represents the coordinates of the point to be evaluated; If the kernel density value is greater than the density threshold, then the concentrated area is judged as a hotspot, otherwise it is a cold spot; Predicting the spatio-temporal information of key disasters through the time prediction model, including: Constructing a spatio-temporal prediction model for key disasters, including: Constructing a spatio-temporal autoregressive moving average model STARMA to predict the disaster intensity at a future time point, that is: ; among which, respectively represent the position within the aggregation area at time the disaster intensity, spatial weight matrix, autoregressive coefficient, moving average coefficient, represents the time lag order; and respectively represent the random error term and the moving error term; Constructing a long short-term memory network LSTM and training the LSTM with historical disaster data; using the trained LSTM to predict the type of disaster occurring at the hotspot at a future time point, that is, the type of key disaster; Using the pre-trained spatio-temporal graph convolutional network ST-GCN to predict the disaster intensity occurring at the hotspot within a future time period, that is, the intensity of key disasters; Performing Kriging interpolation on the prediction results to generate a continuous risk surface to reflect the risk distribution range; Overlaying the prediction results with the hotspot area to identify potential future hotspots.

7. The fuzzy hierarchical evaluation method for spatio-temporal zoning of subway disaster chain evolution according to claim 6, characterized in that Introducing fuzzy logic to comprehensively evaluate the disaster risk based on the spatio-temporal information of key disasters and determining high-risk areas, including: Obtaining the spatio-temporal information of key disasters within the hotspot area and constructing characteristic variables of key disasters; the spatio-temporal information of key disasters includes disaster duration, diffusion speed, hotspot kernel density, and tunnel deformation; Constructing "IF-THEN" fuzzy rules; Based on the fuzzy rules, using the Mamdani fuzzy inference algorithm to output fuzzy conclusions after inputting the characteristic variables of key disasters; Aggregating all the output fuzzy conclusions through the MAX operation to form a fuzzy set; Defuzzifying through the center of gravity method COG to calculate the center of gravity position of the fuzzy set; Selecting the conclusion corresponding to the value of the center of gravity position as the output conclusion; Calculate the comprehensive risk value , where represents the value of the centroid position of the output of the fuzzy rule; represents the rule weight; represents the number of fuzzy rules; Set a high-risk threshold , if then the output hotspot is a high-risk area.

8. The fuzzy hierarchical evaluation method for spatio-temporal zoning of subway disaster chain evolution according to claim 7, characterized in that After determining the high-risk areas, combining the subway structure and evacuation capacity to conduct a secondary evaluation of the disaster risk, including: Obtaining a list of high-risk areas, including the longitude and latitude of high-risk areas, disaster types, and comprehensive risk values; Obtaining subway structure data, including subway tunnel length, subway tunnel depth, support materials, and support structures; Obtaining subway evacuation capacity data, including the location and number of subway entrances and exits, the location and number of emergency exits, the distance from each emergency exit to the high-risk area, and the maximum passenger flow capacity of the subway; Normalizing the subway structure data and subway evacuation capacity data; Obtain the structural safety factor , where respectively represent the structural strength, the structural cross-sectional area, the safety redundancy factor, and the disaster acting force; Obtain the dredging capacity coefficient ; Actual evacuation time , where represents the maximum number of people to be evacuated in the area; respectively represent the number of people that can pass through per second per unit width of the emergency exit, the width of a single emergency exit, and the number of available emergency exits; represents the evacuation attenuation coefficient; represents the minimum evacuation time, ; Obtain the secondary evaluation score by modifying the risk formula ; among them, represents the structure weight and the guidance weight.

9. An electronic device, comprising a processor, a memory connected to the processor, and a communication module, characterized in that, The electronic device is used to execute the fuzzy hierarchical evaluation method for the spatio-temporal partition of the subway disaster chain evolution described in any one of claims 1-8 above.

10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program, and the computer program is executed by a processor to implement the fuzzy hierarchical evaluation method for the spatio-temporal partition of the subway disaster chain evolution described in any one of claims 1-8 above.

Images