Seismic stability analysis method and system for mine rock slopes under freeze-thaw conditions

Abstract

The present invention relates to the technical field of slope stability analysis, and in particular to a method and system for analyzing the seismic stability of a mine rock slope under a freeze-thaw state. The method comprises the following steps: collecting historical seismic stability analysis data of a mine rock slope under a freeze-thaw state, and then establishing a slope stability analysis data set; performing plastic deformation analysis of a mine rock slope under a freeze-thaw state according to the slope stability analysis data set, and then constructing a first stability analysis model; establishing an instability penetration prediction model and a second stability analysis model according to the result of the plastic deformation analysis and the first stability analysis model; acquiring real-time seismic stability analysis data, and then using the instability penetration prediction model and the second stability analysis model to determine the real-time seismic stability of the slope. The present invention can simply, quickly and relatively accurately estimate the seismic stability of a mine rock slope under a freeze-thaw state, and provides more options for stability analysis of mine rock slopes in high-cold areas.

Description

Technical Field

[0001] The present invention relates to the technical field of slope stability analysis, and in particular to a method and system for analyzing seismic stability of a mine rock slope under a freeze-thaw state. Background Art

[0002] In mining engineering practice, especially in cold climates, the stability assessment of rock slopes in mines is a crucial link. These slopes are not only affected by conventional geological conditions, but also face additional challenges brought by freeze-thaw cycles. Freeze-thaw can significantly change the physical and mechanical properties of rocks, such as reducing strength and increasing cracks, which in turn affects the overall stability of the slope. At the same time, earthquakes, as an unpredictable and highly destructive natural disaster, pose another major threat to the stability of slopes. The propagation of seismic waves can stimulate stress redistribution in the slope body, induce geological disasters such as landslides and collapses, and cause great harm to mining production and the surrounding environment. Therefore, the stability assessment of rock slopes in mines under cold climates must comprehensively consider the dual effects of freeze-thaw cycles and earthquakes.

[0003] At present, seismic stability analysis methods for rock slopes in mines and stability analysis methods considering the influence of freeze-thaw conditions have been developed. These methods usually rely on the establishment of complex geological, mechanical and mathematical models to simulate the behavior of slopes under different conditions, which can provide reference for the seismic stability analysis of rock slopes in mines under freeze-thaw conditions. However, these methods generally face the problems of large amount of calculation and long time consumption in practical applications, especially in mines in remote areas. Due to limited monitoring data, insufficient computing resources and hardware equipment, the implementation of these methods becomes particularly difficult.

[0004] In addition, the environmental conditions in remote areas are usually harsh, which places higher demands on the simplicity, rapidity and practicality of the analysis methods. Although the existing complex models and methods can provide relatively accurate analysis results, they are often difficult to meet the needs of rapid response and on-site decision-making in practical applications.

[0005] In summary, although the seismic stability and stability analysis of mine rock slopes under freeze-thaw conditions have achieved certain research results and can be applied to the seismic stability analysis of mine rock slopes under freeze-thaw conditions, in practical applications, there is still a lack of a simple, fast and practical stability analysis method. Therefore, it is of great practical significance and application value to explore a seismic stability analysis method of mine rock slopes under freeze-thaw conditions that can take into account both accuracy and efficiency and adapt to the special environmental conditions of remote areas. Summary of the invention

[0006] In view of the defects in the prior art, the present invention provides a method and system for analyzing the seismic stability of a mine rock slope under freeze-thaw conditions.

[0007] In order to achieve the above-mentioned purpose, in the first aspect, the present invention provides a method for analyzing the seismic stability of a mine rock slope under a freeze-thaw state, the method comprising the following steps: collecting historical seismic stability analysis data of a mine rock slope under a freeze-thaw state, and using the historical seismic stability analysis data to establish a slope stability analysis data set; performing plastic deformation analysis of a mine rock slope under a freeze-thaw state according to the slope stability analysis data set, and then constructing a first stability analysis model; establishing an instability penetration prediction model and a second stability analysis model based on the results of the plastic deformation analysis and the first stability analysis model; acquiring real-time seismic stability analysis data, and then using the instability penetration prediction model and the second stability analysis model to determine the real-time seismic stability of a mine rock slope under a freeze-thaw state. The present invention can simply, quickly and relatively accurately estimate the seismic stability of a mine rock slope under a freeze-thaw state, and provides more options for stability analysis of mine rock slopes in high-cold areas.

[0008] Optionally, the collecting of historical earthquake stability analysis data of a mine rock slope under a freeze-thaw state and establishing a slope stability analysis data set using the historical earthquake stability analysis data comprises the following steps:

[0009] Collect historical earthquake stability analysis data of mine rock slopes under freeze-thaw conditions and pre-process them to obtain accurate historical earthquake stability data;

[0010] The slope stability analysis data set is established using the accurate historical earthquake stability data.

[0011] Furthermore, the present invention provides an accurate and reliable data basis for subsequent seismic stability analysis of mine rock slopes under freeze-thaw conditions by constructing a slope stability analysis data set.

[0012] Optionally, performing plastic deformation analysis of a mine rock slope under a freeze-thaw state according to the slope stability analysis data set to further construct a first stability analysis model comprises the following steps:

[0013] Performing plastic deformation analysis on the mine rock slope according to the slope stability analysis data set, and then determining the sliding surface type of the mine rock slope;

[0014] The rock strips are divided according to the sliding surface type, and then the first stability analysis model of the mine rock slope in the limit equilibrium state is constructed based on the moment balance when the horizontal seismic force acts.

[0015] Furthermore, the present invention provides a theoretical basis for constructing a first stability analysis model by performing plastic deformation analysis on the mine rock slope to determine the sliding surface type, thereby improving the accuracy and reliability of seismic stability analysis of the mine rock slope under freeze-thaw conditions.

[0016] Optionally, the first stability analysis model satisfies the following relationship:

[0017]

[0018] Where N is the number of divided rock blocks, is the height of the i-th rock block, is the angle between the bottom surface of the i-th rock strip and the horizontal plane, H is the slope height of the mine rock slope, R is the radius of the circle where the sliding surface of the mine rock slope is located, is the moment of the earthquake force on the i-th rock block, is the rock bulk density of the mine rock slope, is the width of the i-th rock strip, is the internal friction angle of the rock block, is the effective cohesion of the rock block.

[0019] Furthermore, the first stability analysis model of the present invention is obtained based on moment balance after determining the sliding surface type of the mine rock slope, which makes the model accurate and reliable, providing a solid and reliable foundation for further obtaining the instability penetration prediction model and the second stability analysis model, which is beneficial to improving the accuracy and reliability of the seismic stability analysis results of the mine rock slope under freeze-thaw conditions.

[0020] Optionally, establishing the instability penetration prediction model and the second stability analysis model according to the result of the plastic deformation analysis and the first stability analysis model comprises the following steps:

[0021] Determine key analysis indicators for seismic stability analysis of a mine rock slope under freeze-thaw conditions based on the first stability analysis model;

[0022] Constructing the instability penetration prediction model under different earthquake load amplitudes according to the results of plastic deformation analysis and the key analysis indicators;

[0023] A first mathematical relationship between each of the key analysis indicators at different earthquake load amplitudes is determined, and the first mathematical relationship is used as the second stability analysis model.

[0024] Furthermore, the first stability analysis model is used to determine the key analysis indicators, and then a simpler second stability analysis model is constructed to analyze and predict seismic stability. At the same time, an instability penetration prediction model is constructed to measure the stability prediction error caused by classifying the complex rock slope sliding surface into a specific type of sliding surface, which provides a theoretical basis and data foundation for the subsequent acquisition of accurate and reliable real-time seismic stability.

[0025] Optionally, the key analysis indicators include the internal friction angle of the rock block, the height of the mine rock slope, the rock bulk density of the mine rock slope and the effective cohesion of the rock block;

[0026] The instability penetration prediction model satisfies the following relationship:

[0027]

[0028] in, is the instability penetration factor, H is the slope height of the mine rock slope, The plastic deformation analysis shows that the distance of plastic cracking of the rock slope of a mine under earthquake action is perpendicular to the horizontal direction and starts from the top of the slope. The distance of plastic cracking of the rock slope of a mine under earthquake action, obtained from plastic deformation analysis, is the distance of plastic cracking upward in a direction perpendicular to the horizontal direction starting from the bottom of the slope. and is a constant parameter, is the effective cohesion of the rock block, is the rock bulk density of the mine rock slope, is the internal friction angle of the rock block.

[0029] Furthermore, the present invention calculates the instability penetration factor by using the instability penetration prediction model, which can measure the stability prediction error caused by classifying the complex rock slope sliding surface into a specific sliding surface, and can be used to correct the predicted seismic stability to obtain more accurate real-time seismic stability.

[0030] Optionally, the second stability analysis model satisfies the following relationship:

[0031]

[0032] in, is the internal friction angle of the rock block, H is the slope height of the mine rock slope, is the rock bulk density of the mine rock slope, is the effective cohesion of the rock block, , and For variable parameters.

[0033] Furthermore, the second stability analysis model provided by the present invention involves fewer parameters and is easy to obtain, which can improve the efficiency of seismic stability analysis of mine rock slopes under freeze-thaw conditions.

[0034] Optionally, the acquiring of real-time seismic stability analysis data, and then using the instability penetration prediction model and the second stability analysis model to determine the real-time seismic stability of the mine rock slope in a freeze-thaw state comprises the following steps:

[0035] The first rectangular coordinate system is established with the quotient of the product of effective cohesion, slope height and rock bulk density as the abscissa and the tangent of the internal friction angle as the ordinate.

[0036] Obtain real-time seismic stability analysis data, specifically, obtain real-time seismic load and real-time key analysis index data of the mine rock slope under freeze-thaw state;

[0037] Marking points corresponding to the real-time key analysis indicator data in the first rectangular coordinate system and recording them as real-time data points;

[0038] Obtaining a straight line equation determined by the coordinate origin of the first rectangular coordinate system and the real-time data point, and calculating the coordinates of the intersection of the straight line equation and the second stability analysis model under the real-time seismic load;

[0039] The real-time seismic stability of the mine rock slope under the freeze-thaw state is calculated according to the instability penetration factor, the real-time key analysis indicator data and the intersection coordinates.

[0040] Furthermore, the present invention calculates the real-time seismic stability by drawing the second stability analysis model in the first rectangular coordinate system and calculating the intersection of the second stability analysis model with the straight line equation determined by the coordinate origin and the real-time data points. The whole analysis process is simple and fast and does not require the aid of too much equipment, thereby improving the efficiency of seismic stability analysis of mine rock slopes under freeze-thaw conditions.

[0041] Optionally, the real-time seismic stability satisfies the following relationship:

[0042]

[0043] in, is the real-time seismic stability, is the instability penetration factor, is the real-time internal friction angle in the real-time key analysis indicator data, is the ordinate in the intersection coordinates.

[0044] Furthermore, the present invention uses the instability penetration factor to correct the predicted earthquake stability to obtain more accurate real-time earthquake stability. The calculation method is simple and fast and has high accuracy.

[0045] In the second aspect, the present invention provides a seismic stability analysis system for a mine rock slope in a frozen-thaw state, the seismic stability analysis system for a mine rock slope in a frozen-thaw state comprises: a data acquisition device, a processor, a storage device and a data output device, the storage device comprises a computer-readable storage medium, the computer-readable storage medium stores a computer program, the computer program comprises program instructions, and when the program instructions are executed by the processor, the processor implements the seismic stability analysis method for a mine rock slope in a frozen-thaw state provided by the present invention. The system can stably run the method provided by the present invention, further improve the efficiency of the seismic stability analysis for a mine rock slope in a frozen-thaw state, and improve the practicality of the method provided by the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0046] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings required for use in the embodiments will be briefly introduced below. It should be understood that the following drawings only show certain embodiments of the present application and therefore should not be regarded as limiting the scope. For ordinary technicians in this field, other related drawings can be obtained based on these drawings without paying creative work.

[0047] Figure 1 A schematic flow chart of a method for analyzing seismic stability of a mine rock slope under a freeze-thaw state according to an embodiment of the present invention;

[0048] Figure 2 The plastic distribution cloud diagram of the rock slope of a mine under different earthquake load amplitudes according to an embodiment of the present invention is shown in FIG.

[0049] Figure 2 (a) is a cloud diagram of the plastic distribution of a rock slope in a mine when the earthquake load amplitude is 0.1g according to an embodiment of the present invention;

[0050] Figure 2 (b) is a cloud diagram of the plastic distribution of a rock slope in a mine when the earthquake load amplitude is 0.2g according to an embodiment of the present invention;

[0051] Figure 2 (c) is a cloud diagram of the plastic distribution of the rock slope of a mine when the earthquake load amplitude is 0.3g according to an embodiment of the present invention;

[0052] Figure 2 (d) is a cloud diagram of the plastic distribution of the rock slope of a mine when the earthquake load amplitude is 0.4g according to an embodiment of the present invention;

[0053] Figure 3 It is a schematic diagram of the framework of the seismic stability analysis system of a mine rock slope under freeze-thaw state according to an embodiment of the present invention. DETAILED DESCRIPTION

[0054] The specific embodiments of the present invention will be described in detail below. It should be noted that the embodiments described herein are only for illustration and are not intended to limit the present invention. In the following description, a large number of specific details are set forth in order to provide a thorough understanding of the present invention. However, it is obvious to those of ordinary skill in the art that these specific details do not need to be adopted to implement the present invention. In other examples, in order to avoid confusing the present invention, known circuits, software or methods are not specifically described.

[0055] Throughout the specification, references to "one embodiment," "an embodiment," "an example," or "an example" mean that a particular feature, structure, or characteristic described in conjunction with the embodiment or example is included in at least one embodiment of the present invention. Therefore, the phrases "in one embodiment," "in an embodiment," "an example," or "an example" appearing in various places throughout the specification do not necessarily all refer to the same embodiment or example. In addition, particular features, structures, or characteristics may be combined in one or more embodiments or examples in any suitable combination and / or subcombination. In addition, it should be understood by those of ordinary skill in the art that the figures provided herein are for illustrative purposes and that the figures are not necessarily drawn to scale.

[0056] It should be noted in advance that, in an optional embodiment, except for independent explanations, the same symbols or letters appearing in all formulas have the same meanings and values.

[0057] In an alternative embodiment, see Figure 1 The present invention provides a method for analyzing the seismic stability of a mine rock slope under a freeze-thaw state, the method comprising the following steps:

[0058] S1. Collect historical earthquake stability analysis data of the mine rock slope under freeze-thaw state, and use the historical earthquake stability analysis data to establish a slope stability analysis data set.

[0059] Wherein, step S1 specifically includes the following steps:

[0060] S11. Collect historical earthquake stability analysis data of mine rock slopes under freeze-thaw conditions and perform preprocessing to obtain accurate historical earthquake stability data.

[0061] Specifically, in this embodiment, seismic stability analysis data of the target mine rock slope is collected and preprocessed to obtain accurate historical seismic stability data.

[0062] More specifically, historical seismic stability analysis data include historical internal friction angle, rock bulk density, effective cohesion, pore water pressure, slope height, slope length, slope temperature, slope foot, seismic load amplitude, slope displacement and slope image of the target mine rock slope. The steps of preprocessing seismic stability analysis data include missing value filling, outlier removal, image denoising and image enhancement.

[0063] Furthermore, in other optional embodiments, if the historical seismic stability analysis data of the target mine rock slope is insufficient, a digital twin model of the target mine rock slope can be constructed and the seismic load amplitude and the number of freeze-thaw cycles can be adjusted to simulate the changes of the target mine rock slope under different seismic load amplitudes and freeze-thaw cycles, thereby obtaining sufficient seismic stability analysis data as historical seismic stability analysis data, and the historical seismic stability analysis data can also include other types of data.

[0064] S12. Using the accurate historical earthquake stability data, establish the slope stability analysis data set.

[0065] Specifically, in this embodiment, by constructing a slope stability analysis data set, an accurate and reliable data basis is provided for subsequent seismic stability analysis of mine rock slopes under freeze-thaw conditions.

[0066] S2. Perform plastic deformation analysis of the mine rock slope under freeze-thaw state according to the slope stability analysis data set, and then construct a first stability analysis model.

[0067] Wherein, step S2 specifically includes the following steps:

[0068] S21. Performing plastic deformation analysis on the mine rock slope according to the slope stability analysis data set, and then determining the sliding surface type of the mine rock slope.

[0069] Specifically, in this embodiment, a group of data is randomly selected from the slope stability analysis data set as reference data, and the corresponding plastic distribution cloud map is obtained based on the group of reference data using ABAQUS numerical simulation software.

[0070] More specifically, by adjusting the earthquake load amplitude in the reference data, the ABAQUS numerical simulation software can be used to obtain the plastic distribution cloud map of the target mine rock slope under different earthquake load amplitudes when the number of freeze-thaw cycles is the same, as shown in the following figure: Figure 2 shown.

[0071] Further, in Figure 2In the figure, the lighter the color, the greater the plastic damage of the rock mass. It is not difficult to see that the rock mass of the target mine rock slope suffers plastic damage under the reciprocating cycle of seismic loads, and the damage to the rock mass increases with the increase of the seismic load amplitude. The details are as follows: The rock mass on the upper part of the slope begins to suffer damage first, such as Figure 2 As shown in (a), Figure 2 The seismic load amplitude corresponding to (a) is 0.1g, where g is the gravitational acceleration. As the seismic load amplitude increases, it gradually penetrates the upper rock mass, and the rock mass will crack along the plastic deformation area, forming an arc-shaped sliding surface, such as Figure 2 As shown in (b), (c) and (d) in Figure 2 The seismic load amplitudes corresponding to (b), (c) and (d) are 0.2g, 0.3g and 0.4g respectively.

[0072] This embodiment provides a theoretical basis for constructing a first stability analysis model by performing plastic deformation analysis on the mine rock slope to determine the sliding surface type, thereby improving the accuracy and reliability of seismic stability analysis of the mine rock slope under freeze-thaw conditions.

[0073] S22, dividing the rock strips according to the sliding surface type, and then constructing a first stability analysis model of the mine rock slope under a limit equilibrium state based on moment balance when the horizontal seismic force acts.

[0074] Specifically, in this embodiment, because the sliding surface of the target mine rock slope is an arc-shaped sliding surface under the dual effects of freeze-thaw and earthquake, the Bishop method can be referred to to obtain a calculation relationship for calculating the stability coefficient of the target mine rock slope when considering the horizontal seismic force, and then the stability coefficient of the target mine rock slope is set to 1, and then the obtained calculation relationship is simply transformed to obtain the first stability analysis model, and the first stability analysis model satisfies the following relationship:

[0075]

[0076] Where N is the number of divided rock blocks, is the height of the i-th rock block, is the angle between the bottom of the i-th rock block and the horizontal plane, H is the slope height of the mine rock slope, R is the radius of the circle where the sliding surface of the mine rock slope is located, is the moment of the earthquake force on the ith rock strip, is the rock bulk density of the mine rock slope, is the width of the i-th rock strip, is the internal friction angle of the rock block, is the effective cohesion of the rock strip. In order to simplify the calculation, in this first stability analysis model, the slight differences in rock bulk density, effective cohesion and internal friction angle at different locations of the target mine rock slope are ignored.

[0077] The first stability analysis model is obtained based on moment balance after determining the sliding surface type of the mine rock slope, which makes the model accurate and reliable, providing a solid and reliable foundation for further obtaining the instability penetration prediction model and the second stability analysis model, which is beneficial to improving the accuracy and reliability of the seismic stability analysis results of the mine rock slope under freeze-thaw state.

[0078] S3. Establishing an instability penetration prediction model and a second stability analysis model according to the results of the plastic deformation analysis and the first stability analysis model.

[0079] Wherein, step S3 specifically includes the following steps:

[0080] S31. Determine key analysis indicators for seismic stability analysis of the mine rock slope under freeze-thaw conditions based on the first stability analysis model.

[0081] Specifically, in this embodiment, it is not difficult to see from the first stability analysis model that for the target mine rock slope, if the division of rock strips is fixed, then for the same rock strip, its corresponding , , and will not change if further restrictions are , then when the earthquake force does not change, the only parameters that change in the first stability analysis model are effective cohesion, rock bulk density and internal friction angle, and there is a one-to-one correspondence between effective cohesion, rock bulk density and internal friction angle. According to the first stability analysis model, and There should be a one-to-one correspondence, and the key analysis indicators include the internal friction angle of rock strips, the height of the mine rock slope, the rock density of the mine rock slope and the effective cohesion of rock strips.

[0082] Furthermore, in other optional embodiments, the key analysis indicator may not include the height of the rock slope of the mine.

[0083] S32. Constructing the instability penetration prediction model under different earthquake load amplitudes according to the results of the plastic deformation analysis and the key analysis indicators.

[0084] Specifically, in this embodiment, if Figure 2As shown, although the target mine rock slope will eventually form an arc-shaped sliding surface under the dual effects of freeze-thaw and earthquake, the sliding surface of the target mine rock slope is not arc-shaped before reaching the critical seismic load amplitude. Therefore, if the seismic stability of the target mine rock slope is analyzed completely according to the arc-shaped sliding surface, there will be errors. Therefore, this embodiment calculates the degree of formation of the arc-shaped sliding surface by constructing an instability penetration prediction model, that is, the instability penetration factor, which is used to reflect the error degree of analyzing the seismic stability of the target mine rock slope completely according to the arc-shaped sliding surface, and is used to correct the predicted seismic stability in the future to obtain more accurate real-time seismic stability. The instability penetration prediction model satisfies the following relationship:

[0085]

[0086] in, is the instability penetration factor, The plastic deformation analysis shows that the distance of plastic cracking of the rock slope of a mine under earthquake action is perpendicular to the horizontal direction and starts from the top of the slope. The distance of plastic cracking of the rock slope of a mine under earthquake action, obtained from plastic deformation analysis, is the distance of plastic cracking upward in a direction perpendicular to the horizontal direction starting from the bottom of the slope. and is a constant parameter, which can be obtained by fitting using Matlab, but and The value of will change with the change of the rock slope of the target mine and the amplitude of the seismic load.

[0087] like Figure 2 As shown in Figure 2, the plastic deformation area at the bottom of the slope will extend upward with the increase of the earthquake load amplitude. Figure 2 The distance that the plastic deformation area at the bottom of the mid-slope extends upward is The plastic deformation area at the top of the slope will extend downward with the increase of the earthquake load amplitude. Figure 2 The distance that the plastic deformation area at the top of the mid-slope extends downward is The value of .

[0088] Furthermore, under the same earthquake load amplitude, the values ​​of key analysis indicators are adjusted to change Then, the corresponding plastic distribution cloud diagram is obtained using ABAQUS numerical simulation software to determine and The value of The numerical value is the horizontal axis, The value of is taken as the ordinate, and a series of first-class coordinate points are obtained. Then, these first-class coordinate points are fitted using Matlab to obtain the instability penetration prediction model.

[0089] Furthermore, in this step, the slope height is usually not adjusted when adjusting the values ​​of the key analysis indicators. According to the content described in this step, the instability penetration prediction model under different earthquake load amplitudes can be obtained.

[0090] S33. Determine a first mathematical relationship between each of the key analysis indicators at different earthquake load amplitudes, and use the first mathematical relationship as the second stability analysis model.

[0091] Specifically, in this embodiment, under the same earthquake load amplitude, As an independent variable, As the dependent variable, the independent variable is changed by adjusting the key analysis indicators in the independent variable, and then the adjusted key analysis indicators are brought into the first stability analysis model to calculate the corresponding dependent variable, and then a series of second-type coordinate points with the independent variable as the horizontal coordinate and the dependent variable as the vertical coordinate can be obtained, and then these second-type coordinate points are fitted using Matlab to obtain the second stability analysis model. In order to highlight the relationship between the key analysis indicators, this embodiment writes the second stability analysis model as follows:

[0092]

[0093] in, , and For variable parameters, , and The value of will change with the change of the target mine rock slope and seismic load amplitude. Under the fixed target mine rock slope and seismic load amplitude, , and Yes As an independent variable, As dependent variables, three constants were obtained when using MATLAB fitting to obtain the second stability analysis model.

[0094] Furthermore, in this step, when adjusting the value of the key analysis index, the slope height is usually not adjusted. According to the content described in this step, the second stability analysis model under different earthquake load amplitudes can be obtained.

[0095] The second stability analysis model provided in this embodiment involves fewer parameters and is easy to obtain, which can improve the efficiency of seismic stability analysis of mine rock slopes under freeze-thaw conditions.

[0096] S4. Acquire real-time seismic stability analysis data, and then use the instability penetration prediction model and the second stability analysis model to determine the real-time seismic stability of the mine rock slope under the freeze-thaw state.

[0097] Wherein, step S4 specifically includes the following steps:

[0098] S41. Establish the first rectangular coordinate system with the quotient of the effective cohesion and the product of the slope height and the rock bulk density as the horizontal coordinate and the tangent value of the internal friction angle as the vertical coordinate.

[0099] S42. Obtain real-time seismic stability analysis data, specifically, obtain real-time seismic load and real-time key analysis index data of the mine rock slope under freeze-thaw state.

[0100] S43. Mark the points corresponding to the real-time key analysis indicator data in the first rectangular coordinate system and record them as real-time data points.

[0101] Specifically, in this embodiment, the coordinates of the real-time data points can be expressed as , For real-time effective cohesion in key analytical indicator data, Real-time rock density in real-time key analysis indicator data, Real-time internal friction angle in key analytical indicator data.

[0102] S44, obtaining the equation of a straight line determined by the origin of the first rectangular coordinate system and the real-time data point, and calculating the coordinates of the intersection of the straight line equation and the second stability analysis model under the real-time seismic load.

[0103] Specifically, in this embodiment, the coordinates of the intersection of the origin of the first rectangular coordinate system and the straight line equation determined by the real-time data point and the second stability analysis model under the real-time earthquake load are marked as: , is the horizontal coordinate in the intersection coordinates, is the ordinate of the intersection point.

[0104] S45. Calculate the real-time seismic stability of the mine rock slope under the freeze-thaw state according to the instability penetration factor, the real-time key analysis indicator data and the intersection coordinates.

[0105] Specifically, in this embodiment, according to the proportional strength reduction method, the seismic stability of the slope satisfies:

[0106]

[0107] in, The predicted seismic stability is obtained by further introducing the instability penetration factor under real-time seismic load to correct the predicted seismic stability. Satisfies the following relationship:

[0108]

[0109] Furthermore, the calculations involved in steps S41 to S44 are simple and easy to perform and do not require the use of too many equipment. The final real-time seismic stability also has a high accuracy, thereby improving the efficiency of seismic stability analysis of mine rock slopes under freeze-thaw conditions.

[0110] It should be noted that, in some cases, the actions described in the specification can be performed in a different order and still achieve the desired results. In this embodiment, the order of steps given is only to make the embodiment appear clearer and easier to explain, rather than to limit it.

[0111] In an alternative embodiment, see Figure 3 The present invention also provides a seismic stability analysis system for a mine rock slope in a freeze-thaw state. The seismic stability analysis system for a mine rock slope in a freeze-thaw state comprises: a data acquisition device 1, a processor 2, a storage 3 and a data output device 4. The storage 3 comprises a computer-readable storage medium. A computer program is stored in the computer-readable storage medium. The computer program comprises program instructions. When the program instructions are executed by the processor 3, the processor 3 implements the seismic stability analysis method for a mine rock slope in a freeze-thaw state provided by the present invention, and specifically executes the contents described in steps S1 to S4.

[0112] Specifically, in this embodiment, the data acquisition device 1, the processor 2, the storage 3 and the data output device 4 are interconnected, and the data acquisition device 1 and the data output device 4 respectively include a first display screen and a second display screen. Relevant personnel can input the collected historical earthquake stability analysis data and real-time earthquake stability analysis data into the data acquisition device 1 through the first display screen, and the final real-time earthquake stability will be displayed on the second display screen of the data output device 4.

[0113] The system can stably run the method provided by the present invention, further improve the efficiency of seismic stability analysis of mine rock slopes under freeze-thaw conditions, and improve the practicability of the method provided by the present invention.

[0114] In summary, the method first analyzes the sliding surface type of the mine rock slope under the dual effects of freeze-thaw and earthquake, and then constructs a first stability analysis model to determine the key analysis indicators, and further constructs the instability penetration prediction model and the second stability analysis model under different seismic load amplitudes based on the expression form of the first stability analysis model, and then uses the second stability analysis model to realize the prediction of the seismic stability of the mine rock slope under the freeze-thaw state, and finally uses the instability penetration prediction model to correct the predicted results, and finally obtains a more accurate and reliable real-time seismic stability. After determining the instability penetration prediction model and the second stability analysis model of the mine rock slope, only a few key analysis indicators need to be collected and the real-time seismic stability of the mine rock slope under the freeze-thaw state can be obtained through simple calculation, which improves the efficiency of the real-time seismic stability analysis of the mine rock slope under the freeze-thaw state, and provides more options for the stability analysis of the mine rock slope in the high-cold area. In addition, the system can stably run the method provided by the present invention, further improve the efficiency of the seismic stability analysis of the mine rock slope under the freeze-thaw state, and improve the practicality of the method provided by the present invention.

[0115] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or replace some or all of the technical features therein by equivalents. These modifications or replacements do not make the essence of the corresponding technical solutions deviate from the scope of the technical solutions of the embodiments of the present invention, and they should all be included in the scope of the claims and specification of the present invention.

Claims

1. A method for analyzing the seismic stability of a mine rock slope under freeze-thaw conditions, characterized in that: The steps include: Collect historical earthquake stability analysis data of the mine rock slope under freeze-thaw state, and use the historical earthquake stability analysis data to establish a slope stability analysis data set; Performing plastic deformation analysis of the mine rock slope under freeze-thaw state according to the slope stability analysis data set, and then constructing a first stability analysis model; Establishing an instability penetration prediction model and a second stability analysis model according to the results of the plastic deformation analysis and the first stability analysis model; Acquire real-time seismic stability analysis data, and then use the instability penetration prediction model and the second stability analysis model to determine the real-time seismic stability of the mine rock slope under the freeze-thaw state; The plastic deformation analysis of the mine rock slope under the freeze-thaw state is performed according to the slope stability analysis data set, and then the first stability analysis model is constructed, which includes the following steps: Performing plastic deformation analysis on the mine rock slope according to the slope stability analysis data set, and then determining the sliding surface type of the mine rock slope; Dividing rock strips according to the sliding surface type, and then constructing a first stability analysis model of the mine rock slope in a limit equilibrium state based on moment balance when the horizontal seismic force acts; The first stability analysis model satisfies the following relationship: ， Where N is the number of divided rock blocks, is the height of the i-th rock block, is the angle between the bottom surface of the i-th rock strip and the horizontal plane, H is the slope height of the mine rock slope, R is the radius of the circle where the sliding surface of the mine rock slope is located, is the moment of the earthquake force on the i-th rock block, is the rock bulk density of the mine rock slope, is the width of the i-th rock strip, is the internal friction angle of the rock block, is the effective cohesion of the rock block; The step of establishing the instability penetration prediction model and the second stability analysis model based on the results of the plastic deformation analysis and the first stability analysis model comprises the following steps: Determine key analysis indicators for seismic stability analysis of a mine rock slope under freeze-thaw conditions based on the first stability analysis model; Constructing the instability penetration prediction model under different earthquake load amplitudes according to the results of plastic deformation analysis and the key analysis indicators; Determine a first mathematical relationship between each of the key analysis indicators at different earthquake load amplitudes, and use the first mathematical relationship as the second stability analysis model; The key analysis indicators include the internal friction angle of the rock block, the slope height of the mine rock slope, the rock bulk density of the mine rock slope and the effective cohesion of the rock block; The instability penetration prediction model satisfies the following relationship: ， in, is the instability penetration factor, H is the slope height of the mine rock slope, The plastic deformation analysis shows that the distance of plastic cracking of the rock slope of a mine under earthquake action is perpendicular to the horizontal direction and starts from the top of the slope. The distance of plastic cracking of the rock slope of a mine under earthquake action, obtained from plastic deformation analysis, is the distance of plastic cracking upward in a direction perpendicular to the horizontal direction starting from the bottom of the slope. and is a constant parameter, is the effective cohesion of the rock block, is the rock bulk density of the mine rock slope, is the internal friction angle of the rock block; The second stability analysis model satisfies the following relationship: ， in, is the internal friction angle of the rock block, H is the slope height of the mine rock slope, is the rock bulk density of the mine rock slope, is the effective cohesion of the rock block, , and For variable parameters.

2. The method for analyzing seismic stability of a mine rock slope under freeze-thaw conditions according to claim 1, characterized in that: The steps of collecting historical earthquake stability analysis data of a mine rock slope under a freeze-thaw state and using the historical earthquake stability analysis data to establish a slope stability analysis data set include the following steps: Collect historical earthquake stability analysis data of mine rock slopes under freeze-thaw conditions and pre-process them to obtain accurate historical earthquake stability data; The slope stability analysis data set is established using the accurate historical earthquake stability data.

3. The method for analyzing seismic stability of a mine rock slope under freeze-thaw conditions according to claim 1, characterized in that: The step of obtaining real-time seismic stability analysis data and then using the instability penetration prediction model and the second stability analysis model to determine the real-time seismic stability of the mine rock slope under the freeze-thaw state comprises the following steps: The first rectangular coordinate system is established with the quotient of the product of effective cohesion, slope height and rock bulk density as the abscissa and the tangent of the internal friction angle as the ordinate. Obtain real-time seismic stability analysis data, specifically, obtain real-time seismic load and real-time key analysis index data of the mine rock slope under freeze-thaw state; Marking points corresponding to the real-time key analysis indicator data in the first rectangular coordinate system and recording them as real-time data points; Obtaining a straight line equation determined by the coordinate origin of the first rectangular coordinate system and the real-time data point, and calculating the coordinates of the intersection of the straight line equation and the second stability analysis model under the real-time seismic load; The real-time seismic stability of the mine rock slope under the freeze-thaw state is calculated according to the instability penetration factor, the real-time key analysis indicator data and the intersection coordinates.

4. The method for analyzing seismic stability of a mine rock slope under freeze-thaw conditions according to claim 3, characterized in that: The real-time seismic stability satisfies the following relationship: ， in, is the real-time seismic stability, is the instability penetration factor, is the real-time internal friction angle in the real-time key analysis indicator data, is the ordinate in the intersection coordinates.

5. The seismic stability analysis system of the mine rock slope under the freeze-thaw state is characterized by: The seismic stability analysis system for a mine rock slope in a freeze-thaw state comprises: a data acquisition device, a processor, a storage device and a data output device, wherein the storage device comprises a computer-readable storage medium, wherein a computer program is stored in the computer-readable storage medium, and wherein the computer program comprises program instructions, and when the program instructions are executed by the processor, the processor implements the seismic stability analysis method for a mine rock slope in a freeze-thaw state as described in any one of claims 1 to 4.

Images