Multi-scale seismic geology disaster assessment method considering deformation field and ground motion characteristics

Abstract

The application discloses a kind of multi-scale seismic geology disaster assessment method considering deformation field and ground motion characteristics, comprising the following steps: based on deformation field and seismic ground motion parameter, the risk assessment of regional scale is carried out;With the background of the risk assessment result of regional scale, the dynamic risk score of corridor scale is carried out;Based on the dynamic risk score of corridor scale, the disaster process simulation is carried out, the disaster prediction of site scale is carried out, and the final seismic disaster visualization map and report are generated.The application is coupled with multi-scale dynamic coupling and multi-physical field simulation, combines deformation field, seismic ground motion parameter and numerical simulation technology, through the hierarchical evaluation method of different scales (region, corridor, site), can more comprehensively identify geological disaster risk under the action of earthquake, and adapt to different application scenarios, improve the accuracy and engineering applicability of evaluation.

Description

TECHNICAL FIELD

[0001] The present application belongs to the technical field of earthquake engineering, and particularly relates to a multi-scale seismic geological disaster evaluation method considering deformation field and ground motion characteristics. BACKGROUND

[0002] Traditional geological disaster evaluation methods are mainly based on environmental factors such as topography, geotechnical mechanical properties, and rainfall, and perform hazard evaluation through empirical models or statistical analysis methods. These methods often lack comprehensive consideration of deformation field and ground motion parameters, and are difficult to accurately reflect the geological disaster risk at different spatial scales under the action of earthquakes. In addition, current evaluation methods mostly focus on a single scale, such as regional-scale seismic landslide susceptibility evaluation or local geological disaster risk evaluation of engineering corridors, lacking a multi-scale unified evaluation framework that adapts to different scenarios. Secondly, some studies use historical landslide data and topographic factors to train machine learning models to realize regional landslide susceptibility mapping. The limitations are: insufficient consideration of deformation field information, making it difficult to dynamically evaluate the impact of earthquakes on landslide stability. Limited by historical data, it is difficult to predict the potential seismic triggering landslide risk in the future. Furthermore, some studies use InSAR deformation field information combined with geological conditions to evaluate landslide susceptibility and identify potential high-risk areas, but this method still has limitations: mainly focusing on deformation itself, without combining ground motion parameters, making it difficult to evaluate the potential mechanism of earthquake-triggered landslides. Most studies are limited to regional scale, lacking fine-grained analysis of engineering corridors or specific site scales.

[0003] Therefore, it is an important research direction to establish an effective seismic geological disaster response evaluation technology that comprehensively considers deformation field and ground motion characteristics and performs risk evaluation for different spatial scales, in order to ensure the safety of regional infrastructure and improve disaster response capabilities. SUMMARY

[0004] In view of the above deficiencies in the prior art, the multi-scale seismic geological disaster evaluation method considering deformation field and ground motion characteristics provided by the present application solves the problems of low evaluation accuracy, single scale, insufficient precision, and lack of dynamic evaluation capability in the prior art.

[0005] In order to achieve the above-mentioned application purposes, the technical scheme adopted by the present application is as follows: a multi-scale seismic geological disaster evaluation method considering deformation field and ground motion characteristics, comprising the following steps:

[0006] Performing regional-scale hazard evaluation based on deformation field and ground motion parameters;

[0007] Performing dynamic hazard scoring at the corridor scale based on the results of regional-scale hazard evaluation;

[0008] Based on corridor-scale dynamic hazard scoring, we conduct disaster process simulation, perform site-scale disaster prediction, and generate the final earthquake disaster visualization map and report.

[0009] Furthermore, the specific method for regional-scale hazard assessment based on deformation field and seismic motion parameters is as follows:

[0010] Obtain millimeter-level surface deformation time series, annual average deformation rate, spatial gradient, and distribution of fault-locked segments in the study area; construct a regional deformation field based on the above data.

[0011] A three-dimensional ground motion parameter field, including peak ground acceleration, spectral period, and Arias intensity, was constructed by combining regional ground motion attenuation models, historical strong earthquake records, and fault rupture simulation results.

[0012] The regional deformation field and the three-dimensional seismic motion parameter field are superimposed using the evidence weight method, and a multi-level buffer zone along the fault zone is established based on the superimposed data using the kernel density estimation method.

[0013] Establish a coupled weighted model of regional deformation field and three-dimensional seismic motion parameter field;

[0014] A coupled weighted model of regional deformation field and three-dimensional seismic motion parameter field is used to assess the hazard of multi-level buffer zones along the fault zone and output a regional-scale hazard classification map.

[0015] Furthermore, the coupling weight model of the regional deformation field and the three-dimensional seismic motion parameter field is as follows:

[0016] H R = α·deformation rate + β·peak ground acceleration + γ·topographic slope + δ·lithological vulnerability index

[0017] Among them, H R The risk score represents the comprehensive assessment value of the seismic geological hazard risk in the study area; α is the weighting coefficient of deformation rate; β is the weighting coefficient of peak ground acceleration; γ is the weighting coefficient of topographic slope; δ is the weighting coefficient of lithological vulnerability index; deformation rate is obtained from the regional deformation field, and peak ground acceleration is obtained from the three-dimensional seismic motion parameter field.

[0018] Furthermore, specific methods for conducting dynamic hazard scoring at the corridor scale based on regional-scale hazard assessment results include:

[0019] Screening for high-risk segments with deformation-seismic coupling;

[0020] The parameters of the high-risk deformation-seismic coupling section are refined to obtain the assessment range at the corridor scale;

[0021] The evaluation range of the corridor scale is dynamically scored, i.e., a corridor danger index CSI is calculated;

[0022] According to the corridor danger index CSI, automatic grading is performed in ArcGIS by using natural breakpoints, and a corridor scale danger grading map is output.

[0023] Further, the specific method for screening the deformation-seismic coupling high-risk section is:

[0024] Based on the regional deformation field, a sliding window method is used in combination with spatial clustering analysis to extract a continuous section along the engineering corridor where the deformation rate exceeds a threshold, i.e., a deformation intense section;

[0025] Based on the three-dimensional seismic motion parameter field, a region where the peak ground acceleration is greater than 0.2g is identified, and the spatial intersection of this region and the deformation intense section is calculated to obtain a region that is severely affected by seismic motion;

[0026] A double-index threshold method is used to divide the region that is severely affected by seismic motion into a deformation-seismic coupling high-risk section, wherein the double-index includes the deformation rate and the peak ground acceleration.

[0027] Further, the parameter fine-tuning is a fine-tuning of related parameters in the deformation-seismic coupling high-risk section, and the specific method includes:

[0028] Increasing the time sampling rate to obtain a local large-gradient deformation field of the corridor;

[0029] Taking the terrain slope, elevation, and lithology as inputs of a terrain effect amplification model, the peak ground acceleration is corrected;

[0030] The seismic motion directionality coefficient is quantified by the different angles between the fault strike and the corridor;

[0031] The duration of seismic motion is corrected by using a radiation pattern function.

[0032] Further, the dynamic danger scoring calculation formula is:

[0033] CSI = deformation rate anomaly coefficient × seismic motion amplification coefficient × engineering importance weight

[0034] Wherein, CSI is the corridor danger index; the deformation rate anomaly coefficient is Wherein v i is the deformation rate of the i-th unit, is the average deformation rate of all units in the region, and σ vThe standard deviation of the deformation rate is the standard deviation of the deformation rate; the ground motion amplification factor is the terrain effect amplification model multiplied by the ground motion directionality coefficient; the engineering importance weight is specifically: the engineering importance weight of the high-speed railway is 1.5, the engineering importance weight of the highway is 1.0, and the engineering importance weight of the lifeline engineering is 1.3.

[0035] Further, the specific method for simulating the catastrophic process is:

[0036] The paragraph with CSI>0.6 is taken as the object of site scale analysis, and a three-dimensional finite element model is established;

[0037] The non-saturated seepage model under the rainstorm working condition is simulated by using the Richards equation, and the seismic duration obtained after the parameter fine correction is taken as the input of the three-dimensional finite element model;

[0038] The limit equilibrium method and the strength reduction method are used to calculate the stability coefficient of the slope under the action of earthquake-rainfall coupling, and the state of the slope is divided according to the critical threshold value; the state of the slope includes unstable state, potential unstable state and stable state;

[0039] The Newmark slider calculation model is used to calculate the displacement of the earthquake-triggered landslide for the slope in the unstable state, and if the displacement is greater than the displacement threshold value, it indicates that the slope has obvious deformation;

[0040] The PFC3D is used to simulate the damage process of the slope with obvious deformation, analyze the motion trajectory, impact energy, accumulation range and threat to the project of the landslide, and generate the final seismic disaster visualization map and report.

[0041] The beneficial effects of the present application are:

[0042] 1. Through multi-level data dimension reduction and parameter transmission of regional scale-corridor scale-site scale, the whole chain of seismic geological disaster analysis from macroscopic law to mesoscopic focus to microscopic mechanism is realized for the research area, the one-sidedness problem of traditional single scale evaluation is solved, and the seismic geological disaster risk is more comprehensively identified, so that the analysis result is more accurate and reliable.

[0043] 2. By constructing the coupling weight model of the regional deformation field and the three-dimensional seismic parameter field, the joint analysis of seismic parameters and topographic and geological factors is realized, the correlation between topography and earthquake is revealed, and the seismic geological disaster in the complex topographic region can be more accurately evaluated.

[0044] 3. By combining the seismic duration and the rainfall model for numerical simulation, the real-time influence of the change of pore water pressure on the stability of the slope is considered, and the dynamic response evaluation of the geological disaster risk under different working conditions is realized.

[0045] 4. By outputting visualized hazard classification diagrams at various scales, the earthquake geological hazard assessment results of the study area can be displayed more intuitively, providing a clear, intuitive and reliable basis for engineering protection, which has high practical value. Attached Figure Description

[0046] Figure 1 This is a flowchart illustrating the method proposed in this invention;

[0047] Figure 2 A regional-scale hazard classification map;

[0048] Figure 3 A hazard classification map at the corridor scale;

[0049] Figure 4 This is a schematic diagram simulating a disaster process at the site scale. Detailed Implementation

[0050] The specific embodiments of the present invention are described below to enable those skilled in the art to understand the present invention. However, it should be understood that the present invention is not limited to the scope of the specific embodiments. For those skilled in the art, various changes are obvious as long as they are within the spirit and scope of the present invention as defined and determined by the appended claims. All inventions utilizing the concept of the present invention are protected.

[0051] like Figure 1 As shown, in one embodiment of the present invention, a multi-scale seismic geological hazard assessment method considering deformation field and seismic motion characteristics includes the following steps:

[0052] S1. Regional-scale hazard assessment based on deformation field and seismic motion parameters:

[0053] Millimeter-level surface deformation time series, annual average deformation rate, spatial gradient, and fault zone locking segment distribution on the southeastern edge of the Tibetan Plateau were acquired using InSAR technology, GNSS (Global Navigation Satellite System), or ground monitoring stations. A regional deformation field was constructed based on this data. Millimeter-level surface deformation monitoring can promptly detect minute displacement changes in geological features such as mountains and slopes, identify potential landslide hazard areas, and help reveal the crustal movement status of the study area, providing reliable data support for subsequent geological hazard assessments.

[0054] Combine regional seismic ground motion attenuation model, historical strong earthquake record and fault rupture simulation results to construct three-dimensional seismic ground motion parameter field including peak ground acceleration (PGA), spectral period (T) and Arias intensity (Ia). The seismic ground motion attenuation model can well reflect the variation law of seismic ground motion with distance and magnitude, and the fault rupture simulation results can reflect the influence of seismic rupture process on seismic ground motion. By combining the historical strong earthquake record with the above two, the prediction of seismic ground motion parameters is more refined, and the spatial distribution characteristics of seismic ground motion can be more truly reflected. At the same time, the three-dimensional seismic ground motion parameter field constructed can more comprehensively show the potential influence of earthquake on different structures and facilities through multi-parameter comprehensive evaluation.

[0055] The evidence weight method (Weights of Evidence) is used to stack the regional deformation field and the three-dimensional seismic ground motion parameter field in a grid, and the kernel density estimation method is used to establish multi-level buffer zones (such as 0-5km, 5-10km, >10km) along the fault zone according to the stacked grid data.

[0056] A coupling weight model of regional deformation field and three-dimensional seismic ground motion parameter field is established:

[0057] H R = α · deformation rate + β · peak ground acceleration + γ · topographic slope + δ · lithology brittleness index

[0058] Wherein, H R is the hazard score, representing the comprehensive evaluation value of the seismic geological disaster hazard of the research area; α is the weight coefficient of deformation rate, representing the relative importance of deformation field in the overall evaluation; β is the weight coefficient of peak ground acceleration, representing the relative importance of seismic intensity in the hazard evaluation; γ is the weight coefficient of topographic slope, representing the influence of topographic steepness on geological disaster susceptibility; δ is the weight coefficient of lithology brittleness index, representing the relative influence of different lithology types in the disaster susceptibility evaluation, which is obtained by machine training (such as random forest); the deformation rate is obtained from the regional deformation field, and the peak ground acceleration is obtained from the three-dimensional seismic ground motion parameter field.

[0059] The coupling weight model of regional deformation field and three-dimensional seismic ground motion parameter field is innovatively proposed, which combines deformation field and seismic ground motion parameters, can not only provide the change of crust before and after the earthquake, but also can reflect the energy release and propagation when the earthquake occurs, which lays a solid foundation for subsequent seismic geological disaster evaluation, and further, provides more comprehensive data support for seismic science research and more scientific basis for engineering seismic design.

[0060] The coupling weight model of regional deformation field and three-dimensional ground motion parameter field is used to evaluate the risk of multi-stage buffer zones along the fault zone. Based on the risk score, the natural breakpoints are automatically classified and the regional scale risk classification map (including low, medium and high risk) is output, as shown in FIG. 1. Figure 2

[0061] S2, the risk assessment results of the regional scale are used as the background to conduct dynamic risk scoring of the corridor scale:

[0062] S2-1, screening of deformation-seismic coupling high-risk section:

[0063] Based on the regional deformation field, the sliding window method (for example, the window size is 5x5 pixels) is used in combination with spatial clustering analysis to extract the continuous section of the engineering corridor (such as highway, railway) along which the deformation rate exceeds the threshold value (in this embodiment, the threshold value is 5mm / year), that is, the deformation intense section;

[0064] Based on the three-dimensional ground motion parameter field, the region with peak ground acceleration greater than 0.2g is identified, and the spatial intersection of the region and the deformation intense section is calculated to obtain the region affected by intense ground motion;

[0065] The double-index threshold method is used to divide the region affected by intense ground motion into deformation-seismic coupling high-risk sections, wherein the double-index includes deformation rate and peak ground acceleration. Specifically, the division is as follows: high-risk section—deformation rate>5mm / year and PGA>0.3g; medium-risk section—deformation rate>3mm / year and 0.2g≤PGA≤0.3g; low-risk section—deformation rate<3mm / year or PGA<0.2g.

[0066] S2-2, parameter fine correction of deformation-seismic coupling high-risk section to obtain the evaluation range of corridor scale:

[0067] The small baseline set (SBAS-InSAR) method is used to improve the time sampling rate to quarter, and the local large gradient deformation field of the corridor is obtained;

[0068] The terrain slope, elevation and lithology are used as the input of the topographic amplification factor (TAF) to correct the peak ground acceleration;

[0069] The seismic motion directionality coefficient is quantified by the angle between the fault strike and the corridor;

[0070] The duration of ground motion is corrected by using the radiation pattern function.

[0071] ​By fine-tuning the above parameters, the interaction between seismic activity and crustal deformation can be more accurately reflected, resulting in more accurate corridor evaluation results. In addition, many factors need to be considered for engineering site selection in high-risk earthquake areas. Fine-tuned parameters can also help identify relatively safe areas and provide more comprehensive references for major engineering site selection.

[0072] S2-3, dynamically score the evaluation range of the corridor scale, i.e. calculate the corridor danger index CSI;

[0073] CSI = deformation rate anomaly coefficient x ground motion amplification coefficient x engineering importance weight

[0074] Wherein, CSI is the corridor danger index; the deformation rate anomaly coefficient is Where v i is the deformation rate of the i-th unit, is the average deformation rate of all units in the region, and σ v is the standard deviation of the deformation rate; the ground motion amplification coefficient = topographic effect amplification model x ground motion directionality coefficient; the engineering importance weight is specifically: the engineering importance weight of high-speed railway is 1.5, the engineering importance weight of highway is 1.0 (the engineering importance weight of national highway is 1.1), and the engineering importance weight of lifeline engineering is 1.3.

[0075] S2-4, according to the corridor danger index CSI, automatically classify in ArcGIS using natural breakpoints, and output the corridor scale danger classification map, as shown in Figure 3 .

[0076] S3, based on the dynamic danger score of the corridor scale, carry out disaster process simulation and site scale disaster prediction, and generate the final seismic disaster visualization map and report:

[0077] The paragraph with CSI>0.6 is taken as the object of site scale analysis, and a three-dimensional finite element model is established;

[0078] The unsaturated seepage model under storm conditions simulated by Richards equation and the ground motion duration after parameter fine-tuning are taken as the input of the three-dimensional finite element model; thereby the pore water pressure distribution and change of key horizons (such as near the sliding surface) are monitored. In order to make the prediction result more fine-tuned, Mohr-Coulomb criterion (peak strength parameter and residual strength parameter) can be used simultaneously to input the three-dimensional finite element model combined with strain softening model, so as to realize the consideration of cumulative damage under seismic cyclic load and thus carry out disaster prediction.

[0079] The limit equilibrium method and the strength reduction method are used to calculate the stability coefficient (Fs) of the slope under the coupling effect of earthquake and rainfall, and the state of the slope is divided according to the critical threshold (the critical threshold of the stability coefficient of the slope is generally determined according to engineering geology and research experience); the state of the slope includes unstable state, potential unstable state and stable state; in this embodiment, the specific state division is: if Fs≥1.20, it is a stable state; if 1.00≤Fs<1.20, it is a potential unstable state; if Fs<1.00, it is an unstable state; if Fs<1.2, the landslide deformation trend needs to be further evaluated.

[0080] The Newmark slider calculation model is used to calculate the landslide displacement triggered by the earthquake for the slope in the unstable state, to determine whether the landslide will have obvious deformation, and if the landslide displacement is greater than a displacement threshold (for example, 18 cm), it indicates that the slope has obvious deformation and may develop into a large-scale damage, and damage mode simulation needs to be performed;

[0081] The PFC3D is used to simulate the damage process of the slope with obvious deformation, to analyze the motion trajectory, impact energy, accumulation range and threat to the project of the landslide, as shown in Figure 4 Figures (a), (b) and (c) represent seismic duration diagrams, and (a-1), (b-1) and (c-1) represent the deformation process of the slope under different seismic durations. Figure (a-1) shows that after the application of the earthquake, the landslide deforms at 2s-8s, the green material represents the area where the landslide first produces a crack, and the crack rapidly expands to form a tensile stress concentration zone. Figure (b-1) shows that at 9s-16s, the rear edge of the landslide rapidly cracks, shear failure occurs at the slope toe, the sliding zone is penetrated at this time, and the landslide body begins to accelerate and disintegrate. Figure (c-1) shows that after 20s, the front edge of the landslide particles dive towards the valley, and then accumulate at the valley. Through the damage process simulation of the slope with obvious deformation, the whole process from the initial deformation to the final destruction of the slope can be clearly seen by combining theoretical analysis and actual observation, the deformation characteristics and damage mode at different stages are understood, and more comprehensive basis is provided for the evaluation of earthquake geological disasters.

[0082] In summary, the present application can more comprehensively identify the geological disaster risk under the action of the earthquake through multi-scale dynamic coupling and multi-physical field simulation, combined with deformation field, seismic parameters and numerical simulation technology, through the hierarchical evaluation method of different scales (region, corridor, site), adapt to different application scenarios, and improve the accuracy and engineering applicability of the evaluation.

Claims

1. A multi-scale seismic geology disaster assessment method considering deformation field and ground motion characteristics, characterized in that, The method comprises the following steps: Performing regional-scale risk assessment based on the deformation field and ground motion parameters; Performing corridor-scale dynamic risk scoring in the background of the regional-scale risk assessment results; Performing disaster process simulation based on the corridor-scale dynamic risk scoring, performing site-scale disaster prediction, and generating final seismic disaster visualization maps and reports; The specific method for performing regional-scale risk assessment based on the deformation field and ground motion parameters is: Obtain the millimeter-level surface deformation time series, annual average deformation rate, spatial gradient, and distribution of locked segments of the fault zone in the study area; construct the regional deformation field according to the millimeter-level surface deformation time series, annual average deformation rate, spatial gradient, and distribution of locked segments of the fault zone in the study area; Combine the regional ground motion attenuation model, historical strong earthquake records, and fault rupture simulation results to construct a three-dimensional ground motion parameter field including peak ground acceleration, spectral period, and Arias intensity; Use the evidence weight method to perform grid stacking of the regional deformation field and the three-dimensional ground motion parameter field, and use the kernel density estimation method to establish a multi-level buffer zone along the fault zone according to the grid-stacked data; Establish a coupling weight model of the regional deformation field and the three-dimensional ground motion parameter field; Use the coupling weight model of the regional deformation field and the three-dimensional ground motion parameter field to perform risk assessment on the multi-level buffer zone along the fault zone, and output a regional-scale risk classification map; The coupling weight model of the regional deformation field and the three-dimensional ground motion parameter field is: H R = a • deformation rate + b • peak ground acceleration + g • terrain slope + d • lithology vulnerability index wherein H R is a hazard score, representing a comprehensive evaluation value of seismic geological disaster hazard of the study area; a is a weight coefficient of the deformation rate; β is a weight coefficient of the peak ground acceleration; γ is a weight coefficient of the topographic slope; δ is a weight coefficient of the lithology vulnerability index; the deformation rate is obtained from the regional deformation field, and the peak ground acceleration is obtained from the three-dimensional ground motion parameter field; The specific method for performing corridor-scale dynamic risk scoring in the background of the regional-scale risk assessment results includes: Screening deformation-seismic coupling high-risk segments; Performing parameter fine-tuning on the deformation-seismic coupling high-risk segments to obtain an evaluation range at the corridor scale; Performing dynamic risk scoring on the evaluation range at the corridor scale, i.e., calculating the corridor risk index CSI; Using natural breakpoints in ArcGIS to automatically classify the corridor risk index CSI, and outputting a corridor-scale risk classification map.

2. The multi-scale seismic geological disaster assessment method considering deformation field and ground motion characteristics according to claim 1, wherein the specific method for screening deformation-seismic coupling high-risk segments is: Based on the regional deformation field, use the sliding window method combined with spatial clustering analysis to extract continuous segments along the engineering corridor where the deformation rate exceeds the threshold value, i.e., the deformation intense segment; Based on the three-dimensional ground motion parameter field, identify the region where the peak ground acceleration is greater than 0.2g, and calculate the spatial intersection of the region and the deformation intense segment to obtain the region affected by intense seismic motion; Use the double-index threshold method to divide the region affected by intense seismic motion into deformation-seismic coupling high-risk segments, wherein the double-index includes the deformation rate and the peak ground acceleration.

3. The multi-scale seismic geological disaster assessment method considering deformation field and ground motion characteristics according to claim 2, wherein the parameter fine-tuning is fine-tuning of related parameters in the deformation-seismic coupling high-risk segments, and the specific method includes: Increase the time sampling rate to obtain a local large-gradient deformation field of the corridor; Use the terrain slope, elevation, and lithology as inputs of the terrain effect amplification model to correct the peak ground acceleration; The directionality coefficient of ground motion is quantified by the angle between the fault strike and the corridor; The duration of ground motion is modified by the radiation pattern function.

4. The multi-scale seismic geological disaster assessment method considering deformation field and ground motion characteristics according to claim 3, wherein the dynamic risk score calculation formula is: CSI = deformation rate anomaly coefficient × ground motion amplification coefficient × engineering importance weight.

5. The multi-scale seismic geological disaster assessment method considering deformation field and ground motion characteristics according to claim 4, wherein the specific method for carrying out disaster process simulation is: taking the paragraph with CSI > 0.6 as the analysis object of site scale, and establishing a three-dimensional finite element model; taking the duration of ground motion obtained by simulating the unsaturated seepage model under storm conditions by using the Richards equation and after parameter fine correction as the input of the three-dimensional finite element model; calculating the stability coefficient of the slope under the action of earthquake-rainfall coupling by using the limit equilibrium method and the strength reduction method, and dividing the state of the slope according to the critical threshold value; the state of the slope includes unstable state, potential unstable state and stable state; calculating the landslide displacement triggered by the earthquake by using the Newmark slider calculation model for the slope in the unstable state, and indicating that the slope has obvious deformation if the landslide displacement is greater than the displacement threshold value; simulating the failure process of the slope with obvious deformation by using PFC3D, analyzing the motion trajectory, impact energy, accumulation range of the landslide and the threat to the project, and generating the final seismic disaster visualization map and report. wherein, CSI is the corridor risk index; the abnormal coefficient of deformation rate is wherein is the deformation rate of the i unit, is the average deformation rate of all units in the region, is the standard deviation of the deformation rate; the ground motion amplification coefficient = topographic effect amplification model × ground motion directionality coefficient; the engineering importance weight is specifically: the engineering importance weight of high-speed railway is 1.5, the engineering importance weight of highway is 1.0, and the engineering importance weight of lifeline engineering is 1.

3. ​ ​ ​ ​ ​ ​

Images