Working side slope landslide prevention and control method, system and equipment and storage medium

Abstract

The invention discloses a working side slope landslide prevention and control method, system and device and a storage medium, and relates to the field of geological disaster prevention and control, and the method comprises the steps: extracting geological parameters of a soft layer from a three-dimensional geological model constructed based on remote sensing surveying and mapping data, and obtaining a shear strength parameter range of the soft layer according to the geological parameters, generating a two-dimensional slope displacement vector cloud picture and a two-dimensional slope strain rate cloud picture according to the three-dimensional geological model and the shear strength parameter range, generating a coverage range after the landslide occurs in the predicted landslide area according to the geological parameters, the shear strength parameter range, the two-dimensional slope displacement vector cloud picture and the two-dimensional slope strain rate cloud picture, and calculating the coverage range according to the coverage range. And displaying the preset landslide prevention and control method matched with the maximum safety distance. Based on the characteristics of high detection depth and precision of remote sensing surveying and mapping data, the influence of deep geological parameters on the landslide is fully considered in the obtained coverage range, and the landslide dangerous area delimiting precision and the landslide prevention and control measure configuration precision are improved.

Description

Technical Field

[0001] The present application relates to the technical field of geological disaster prevention and control, and in particular to a method, system, equipment and storage medium for preventing and controlling landslides on a working slope. Background Art

[0002] The stability of the working slope is an important factor affecting the safety of open-pit mining. When the working slope is not stable, it will cause slope landslides, resulting in major accidents, seriously affecting the safety risks and economic benefits of open-pit mines. Existing technologies mostly use a combination of manual visual inspection and sensor detection to evaluate the stability of the working slope, and based on the evaluation results, define the landslide risk area and formulate corresponding landslide prevention and control strategies.

[0003] However, the existing method is limited by the blind spot of vision and the detection depth of the sensor, and cannot accurately detect the strength changes of the internal structure of the working slope, resulting in low accuracy in determining the landslide hazard area, which in turn affects the configuration accuracy of the landslide prevention and control measures on the working slope. Summary of the invention

[0004] In view of the above problems, this application provides a method, system, device and storage medium for preventing and controlling landslides on working slopes, so as to achieve the purpose of improving the accuracy of determining landslide hazard areas and the accuracy of configuring landslide prevention and control measures. The specific scheme is as follows:

[0005] The first aspect of the present application provides a method for preventing and controlling landslides on a working slope, comprising:

[0006] Constructing a three-dimensional geological model of the target working slope based on remote sensing mapping data of the target working slope, and extracting geological parameters of the weak layer in the target working slope from the three-dimensional geological model;

[0007] Obtaining a shear strength parameter range of the weak layer according to the geological parameters using a preset limit equilibrium algorithm, wherein the parameter types of the shear strength parameter range include a cohesion range and an internal friction angle range;

[0008] Using a preset cloud map generation software to generate a two-dimensional slope displacement vector cloud map and a two-dimensional slope strain rate cloud map based on the three-dimensional geological model and the shear strength parameter range, the two-dimensional slope displacement vector cloud map represents the predicted landslide area of ​​the target working slope, and the two-dimensional slope displacement vector cloud map represents the landslide earthwork volume of the predicted landslide area;

[0009] Based on the geological parameters, the shear strength parameter range, the two-dimensional slope displacement vector cloud map and the two-dimensional slope strain rate cloud map, generating a coverage range after the landslide occurs in the predicted landslide area;

[0010] The maximum safety distance is determined based on the coverage range, and the coverage range and the preset landslide prevention and control method corresponding to the applicable distance with the smallest difference from the maximum safety distance are displayed.

[0011] In a possible implementation, extracting geological parameters of the weak layer in the target working slope from the three-dimensional geological model includes:

[0012] Determine the stratum at the junction of the strongly weathered mudstone and the strongly weathered sandstone in the three-dimensional geological model as the weak layer, and extract the boundary length of the weak layer;

[0013] The weak layer is divided based on the boundary length and the preset division unit length to obtain a plurality of measurement areas corresponding to respective equal boundary lengths; parameters are extracted from each of the measurement areas in the three-dimensional geological model to obtain geological parameters of the weak layer.

[0014] In a possible implementation, the preset limit equilibrium algorithm includes: a first preset limit equilibrium algorithm and a second preset limit equilibrium algorithm, wherein the preset limit equilibrium algorithm is used to obtain the shear strength parameter range of the weak layer according to the geological parameters, and the parameter types of the shear strength parameter range include a cohesion range and an internal friction angle range, including:

[0015] For each of the measurement areas: based at least on the initial stability coefficient and the angle of repose in the geological parameters in the measurement area, the cohesion of the measurement area is calculated using the first preset limit equilibrium algorithm; and the internal friction angle of the measurement area is calculated based at least on the cohesion of the measurement area using the second preset limit equilibrium algorithm;

[0016] The cohesion and the internal friction angle of each of the measurement areas are summarized to obtain the shear strength parameter range of the weak layer including the cohesion range and the internal friction angle range.

[0017] In a possible implementation, the calculating the cohesion of the measurement area by using the first preset limit equilibrium algorithm based on at least the initial stability coefficient and the angle of repose in the geological parameters in the measurement area includes:

[0018] According to the initial stability coefficient F in the geological parameters of the measurement area i i s and the repose angle θ, by the formula of the first preset limit equilibrium algorithm:

[0019]

[0020] The cohesion C of the measurement area i is obtained i , wherein the Li is the preset division unit length, the W i It is the weight of the soil strip in the geological parameters of the measurement area i.

[0021] In a possible implementation, the calculating the internal friction angle of the measurement area based on the cohesion of the measurement area by using the second preset limit equilibrium algorithm includes:

[0022]

[0023] Obtain the internal friction angle φ of the measurement area i i , wherein the L i is the preset division unit length, the C i is the cohesive force, the W i is the weight of the soil strip in the geological parameters of the measurement area i, and u i is the pore water pressure in the geological parameters of the measurement area i, and Q i is the horizontal force on the measured area i in the geological parameters, and the e i is the normal force on the measured area i in the geological parameters, R i is the sliding surface radius in the geological parameters of the measurement area i, θ is the repose angle in the geological parameters of the measurement area i, and m θ is the calculation coefficient corresponding to the angle of repose in the geological parameters of the measurement area i, and the F i s It is the initial stability coefficient among the geological parameters of the measurement area i.

[0024] In a possible implementation, generating the coverage range of the predicted landslide area after the landslide occurs based on the geological parameters, the shear strength parameter range, the two-dimensional slope displacement vector cloud map and the two-dimensional slope strain rate cloud map includes:

[0025] The geological parameters, the shear strength parameter range, the two-dimensional slope displacement vector cloud map and the two-dimensional slope strain rate cloud map are imported into a preset particle flow simulation model, so that the preset particle flow simulation model simulates the landslide process of the preset landslide area, and outputs a landslide coverage cloud map after the landslide occurs in the predicted landslide area when the simulation result meets the preset conditions; and extracts the parameters of the coverage range from the landslide coverage cloud map.

[0026] In a possible implementation, before the preset particle flow simulation model simulates the landslide process of the preset landslide area, the method further includes:

[0027] The blasting vibration influence coefficient is introduced into the preset particle flow simulation model.

[0028] The second aspect of the present application provides a working slope landslide prevention and control system, comprising:

[0029] A parameter extraction module is used to construct a three-dimensional geological model of the target working slope based on the remote sensing mapping data of the target working slope, and to extract geological parameters of the weak layer in the target working slope from the three-dimensional geological model;

[0030] A parameter range extraction module, used to obtain the shear strength parameter range of the weak layer according to the geological parameters using a preset limit equilibrium algorithm, wherein the parameter types of the shear strength parameter range include a cohesion range and an internal friction angle range;

[0031] A cloud map generation module, used for generating a two-dimensional slope displacement vector cloud map and a two-dimensional slope strain rate cloud map based on the three-dimensional geological model and the shear strength parameter range using a preset cloud map generation software, wherein the two-dimensional slope displacement vector cloud map represents the predicted landslide area of ​​the target working slope, and the two-dimensional slope displacement vector cloud map represents the landslide earthwork volume of the predicted landslide area;

[0032] A coverage range determination module, used for generating a coverage range after a landslide occurs in the predicted landslide area based on the geological parameters, the shear strength parameter range, the two-dimensional slope displacement vector cloud map and the two-dimensional slope strain rate cloud map;

[0033] The data display module is used to determine the maximum safety distance based on the coverage range, and to display the coverage range and the preset landslide prevention and control method corresponding to the applicable distance with the smallest difference from the maximum safety distance.

[0034] In a possible implementation, the parameter extraction module is set to extract the geological parameters of the weak layer in the target working slope from the three-dimensional geological model as follows:

[0035] Determine the stratum at the junction of the strongly weathered mudstone and the strongly weathered sandstone in the three-dimensional geological model as the weak layer, and extract the boundary length of the weak layer;

[0036] The weak layer is divided based on the boundary length and the preset division unit length to obtain a plurality of measurement areas corresponding to respective equal boundary lengths; parameters are extracted from each of the measurement areas in the three-dimensional geological model to obtain geological parameters of the weak layer.

[0037] In a possible implementation, the parameter range extraction module is set to: when the parameter range of the weak layer is obtained according to the geological parameters using a preset limit equilibrium algorithm, and the parameter types of the shear strength parameter range include a cohesion range and an internal friction angle range:

[0038] For each of the measurement areas: based at least on the initial stability coefficient and the angle of repose in the geological parameters in the measurement area, the cohesion of the measurement area is calculated using a first preset limit equilibrium algorithm; based at least on the cohesion of the measurement area, the internal friction angle of the measurement area is calculated using a second preset limit equilibrium algorithm, the preset limit equilibrium algorithm including the first preset limit equilibrium algorithm and the second preset limit equilibrium algorithm;

[0039] The cohesion and the internal friction angle of each of the measurement areas are summarized to obtain the shear strength parameter range of the weak layer including the cohesion range and the internal friction angle range.

[0040] In a possible implementation, the parameter range extraction module is set to: when calculating the cohesion of the measurement area using the first preset limit equilibrium algorithm based on at least the initial stability coefficient and the angle of repose in the geological parameters in the measurement area:

[0041] According to the initial stability coefficient F in the geological parameters of the measurement area i i s and the repose angle θ, by the formula of the first preset limit equilibrium algorithm:

[0042]

[0043] The cohesion C of the measurement area i is obtained i , wherein the L i is the preset division unit length, the W i It is the weight of the soil strip in the geological parameters of the measurement area i.

[0044] In a possible implementation, when the parameter range extraction module calculates the internal friction angle of the measurement area based on the cohesion of the measurement area using the second preset limit equilibrium algorithm, it is set to:

[0045]

[0046] Obtain the internal friction angle φ of the measurement area i i , wherein the L i is the preset division unit length, the C i is the cohesive force, the W iis the weight of the soil strip in the geological parameters of the measurement area i, and u i is the pore water pressure in the geological parameters of the measurement area i, and Q i is the horizontal force on the measured area i in the geological parameters, and the e i is the normal force on the measured area i in the geological parameters, R i is the sliding surface radius in the geological parameters of the measurement area i, θ is the repose angle in the geological parameters of the measurement area i, and m θ is the calculation coefficient corresponding to the angle of repose in the geological parameters of the measurement area i, and the F i s It is the initial stability coefficient among the geological parameters of the measurement area i.

[0047] In a possible implementation, the coverage range determination module is set to:

[0048] The geological parameters, the shear strength parameter range, the two-dimensional slope displacement vector cloud map and the two-dimensional slope strain rate cloud map are imported into a preset particle flow simulation model, so that the preset particle flow simulation model simulates the landslide process of the preset landslide area, and outputs a landslide coverage cloud map after the landslide occurs in the predicted landslide area when the simulation result meets the preset conditions; and extracts the parameters of the coverage range from the landslide coverage cloud map.

[0049] In a possible implementation, the coverage range determination module is further configured to:

[0050] Before the preset particle flow simulation model simulates the landslide process in the preset landslide area, the blasting vibration influence coefficient is introduced into the preset particle flow simulation model.

[0051] A third aspect of the present application provides an electronic device, comprising at least one processor and a memory connected to the processor, wherein:

[0052] The memory is used to store computer programs;

[0053] The processor is used to execute the computer program so that the electronic device can implement the slope landslide prevention and control method of the first aspect or any implementation method of the first aspect.

[0054] A fourth aspect of the present application provides a computer storage medium, which carries one or more computer programs. When the one or more computer programs are executed by an electronic device, the electronic device can help the slope landslide prevention and control method according to the first aspect or any implementation method of the first aspect.

[0055] Through the above-mentioned technical scheme, the present application provides a working slope landslide prevention and control method, system, equipment and storage medium, which constructs a three-dimensional geological model of the target working slope by configuring remote sensing sensing and mapping data based on the target working slope, and extracts the geological parameters of the weak layer in the target working slope from the three-dimensional geological model. The remote sensing and mapping data can truly reflect the characteristics of the subsurface stratigraphic structure, improve the accuracy of determining the maximum safe distance, and thus improve the accuracy of delineating the landslide danger area. Subsequently, the shear strength parameter range of the weak layer is obtained according to the geological parameters by configuring the preset limit equilibrium algorithm, and the parameter types of the shear strength parameter range include the cohesion range and the internal friction angle range. Based on the three-dimensional geological model and the shear strength parameter range, a two-dimensional slope displacement vector cloud map and a two-dimensional slope strain rate cloud map are generated. The two-dimensional slope displacement vector cloud map represents the predicted landslide area of ​​the target working slope, and the two-dimensional slope displacement vector cloud map represents the landslide earthwork volume of the predicted landslide area, so as to realize the prediction of the area where the target working slope may occur, and also realize the prediction of the possible coverage range after the predicted landslide area occurs, so that the present application improves the accuracy of delineating the landslide hazard area compared with the manual analysis method. Finally, by displaying the coverage range and the preset landslide prevention and control method corresponding to the applicable distance with the smallest difference from the maximum safety distance, by pre-configuring the corresponding preset landslide prevention and control method for different applicable distances, and when the maximum safety distance is obtained, the preset landslide prevention and control method corresponding to the applicable distance with the smallest difference from the obtained maximum safety distance is displayed, so as to shorten the generation time of the landslide prevention and control method while improving the configuration accuracy of the landslide prevention and control measures. It can be seen that this application improves the accuracy of determining landslide risk areas and configuring landslide prevention and control measures. BRIEF DESCRIPTION OF THE DRAWINGS

[0056] The above and other features, advantages and aspects of the embodiments of the present disclosure will become more apparent with reference to the following detailed description in conjunction with the accompanying drawings. Throughout the accompanying drawings, the same or similar reference numerals represent the same or similar elements. It should be understood that the drawings are schematic and the originals and elements are not necessarily drawn to scale.

[0057] Figure 1 A flow chart of a method for preventing and controlling landslides on a working slope provided for this application;

[0058] Figure 2 A mine elevation map provided for this application;

[0059] Figure 3 A two-dimensional slope displacement vector cloud map provided for this application;

[0060] Figure 4 A two-dimensional slope strain rate cloud map provided for this application;

[0061] Figure 5 A schematic diagram of a measurement area provided for this application;

[0062] Figure 6 A landslide coverage cloud provided for this application;

[0063] Figure 7 A block diagram of a working aid slope landslide prevention and control system provided for this application;

[0064] Figure 8 A schematic diagram of the structure of an electronic device provided in this application. DETAILED DESCRIPTION

[0065] The following describes the embodiments of the present application in conjunction with the drawings in the embodiments of the present application. The terms used in the implementation method section of the present application are only used to explain the specific embodiments of the present application, and are not intended to limit the present application.

[0066] The embodiments of the present application are described below in conjunction with the accompanying drawings. Those skilled in the art will appreciate that, with the development of technology and the emergence of new scenarios, the technical solutions provided in the embodiments of the present application are also applicable to similar technical problems.

[0067] The terms "first", "second", etc. in the specification and claims of the present application and the above-mentioned drawings are used to distinguish similar objects, and need not be used to describe a specific order or sequential order. It should be understood that the terms used in this way can be interchangeable under appropriate circumstances, which is only to describe the distinction mode adopted by the objects of the same attributes when describing in the embodiments of the present application. In addition, the terms "including" and "having" and any of their variations are intended to cover non-exclusive inclusions, so that the process, method, system, product or equipment comprising a series of units need not be limited to those units, but may include other units that are not clearly listed or inherent to these processes, methods, products or equipment.

[0068] The first aspect of the present application provides a method for preventing and controlling landslides on a working slope, such as Figure 1 As shown in the figure, the methods for preventing and controlling slope landslides in this work include:

[0069] S101, constructing a three-dimensional geological model of the target working slope based on the remote sensing mapping data of the target working slope, and extracting geological parameters of the weak layer in the target working slope from the three-dimensional geological model.

[0070] It should be noted that in actual application scenarios, the existing method of implementing the stability assessment of the working side slope by combining manual visual inspection with sensor detection is as follows: the inspector conducts a visual inspection of the surface conditions of the target working side slope at the mine observation point, determines the detection points with potential risks, then opens a probe hole at the detection point and places a sensor to collect the stratigraphic information at the detection point, and then evaluates the stability of the working side slope through manual analysis. However, due to the large area of ​​the target working side and the presence of blind spots, the inspector is unable to fully observe all parts of the target working side slope. At the same time, due to construction conditions, the depth of the probe hole is shallow, and accurate stratigraphic information cannot be collected. This makes the accuracy of the existing manual visual inspection combined with sensor detection to assess the stability of the working side slope low, which in turn leads to a reduction in the accuracy of demarcation of landslide hazard areas and configuration of landslide prevention and control measures.

[0071] It should be noted that in actual application scenarios, the above remote sensing mapping data can be obtained by extracting and processing the reflected electromagnetic radiation of the target working slope using drones or satellites based on remote sensing (RS) technology. Address data. The above remote sensing mapping data may include various types of data, including but not limited to: stratum depth, Stratigraphic Quality, pore water pressure, formation elevation, etc. Since electromagnetic waves can penetrate the surface and penetrate deep into the formation, the remote sensing mapping data obtained can reflect the structure inside the formation. Therefore, this application constructs a three-dimensional geological model of the target working slope by configuring the remote sensing mapping data based on the target working slope, so that compared with the existing manual visual observation or sensor detection method, the internal geological parameters of the target working slope can be accurately obtained.

[0072] In a possible implementation, the remote sensing mapping data may be data obtained by aggregating multiple types of mapping data, for example, Figure 2 Shown is an elevation map of the mining area collected using satellite. Figure 2 The altitude of each working slope is marked with a gradient color from blue (indicating that the altitude of the area is 940 meters) to red (indicating that the altitude of the area is 1090 meters). After obtaining the above elevation map, based on the regional coordinates of the target working slope, the altitude of the target working slope is extracted from the corresponding regional coordinate area in the elevation map.

[0073] It should be noted that, in actual application scenarios, the above-mentioned three-dimensional geological model can be a digital model for representing the external and internal geological information of the above-mentioned target working slope, which is generated based on the three-dimensional geological modeling technology and the above-mentioned remote sensing mapping data. The above-mentioned three-dimensional geological model can be a three-dimensional grid body, which represents the structure and geometry of the stratum, and each node in each grid stores multiple types of geological parameters of the stratum position corresponding to the grid, such as porosity, permeability, water saturation, etc.

[0074] It should be noted that in actual application scenarios, the above-mentioned weak layer is also called weak rock stratum, which refers to bedrock with low strength and elastic modulus. For the working slope of the mine, since the working slope is mostly mudstone with a small thickness, the rock layer inside the working slope is reduced when it is soaked in groundwater for a long time, resulting in a weak layer. As the mine peeling progresses, when the foot of the working slope is excavated, the slope releases the internal weak layer due to the release of slope stress. Since the weak layer has a low strength, it slides when it is squeezed by the slope. At this time, the slope will slide with the weak layer to the free surface of the weak layer. In this process, the excavation of the foot will cause the anti-sliding force of the front edge of the slope to decrease, and finally cause the working slope to become unstable and cause a landslide. It can be seen that the strength of the weak layer is an important influencing parameter that determines whether a landslide occurs on the working slope. Therefore, this application extracts the geological parameters of the weak layer in the target working slope from the three-dimensional geological model, and subsequently delineates the coverage area and formulates landslide prevention and control methods based on the geological parameters of the weak layer, thereby improving the delineation accuracy of the landslide danger zone and the configuration accuracy of the landslide prevention and control measures.

[0075] S102. Obtaining a shear strength parameter range of the weak layer according to geological parameters using a preset limit equilibrium algorithm, wherein parameter types of the shear strength parameter range include a cohesion range and an internal friction angle range.

[0076] It should be noted that, in actual application scenarios, the above-mentioned preset limit equilibrium algorithm can be an algorithm based on the reverse construction of the limit equilibrium method. The limit equilibrium method analyzes the stress state of the slope under various failure modes (such as groundwater immersion) according to the principle of static equilibrium, and calculates the stability coefficient of the block according to the ratio of the shear force of the block (working slope) along the failure surface acting on the potential failure surface (weak layer) in the rock and soil body (working slope) and the shear force of the block along the failure surface.

[0077] It should be noted that in practical application scenarios, in the field of geological research, shear strength is an important mechanical parameter used to measure the stability of rock mass, which characterizes the critical value of rock mass landslide. The above shear strength includes two types of parameters: cohesion (also known as adhesion) and internal friction angle.

[0078] It should be noted that the present application configures and utilizes a preset limit equilibrium algorithm to obtain the shear strength parameter range of the weak layer, including two parameter types: cohesion range and internal friction angle range, according to geological parameters, so as to simulate the landslide process of the target working slope based on the shear strength parameter range in subsequent steps.

[0079] S103. Generate a two-dimensional slope displacement vector cloud map and a two-dimensional slope strain rate cloud map based on the three-dimensional geological model and the shear strength parameter range using a preset cloud map generation software. The two-dimensional slope displacement vector cloud map represents the predicted landslide area of ​​the target working slope. The two-dimensional slope displacement vector cloud map represents the landslide earthwork volume in the predicted landslide area.

[0080] Those skilled in the art can understand that, in actual application scenarios, the above-mentioned preset cloud map generation software can be software for simulating geological change processes in the field of geological design or exploration, such as slope radar displacement cloud map automatic generation software, three-dimensional excavation and slope finite element analysis software (RS3), three-dimensional finite difference program (Fast Lagrangian Analysis of Continua 3D, FLAC 3D), etc. This application does not make too many restrictions and redundant descriptions on the specific selection and operation process of the above-mentioned preset cloud map generation software.

[0081] It should be noted that in actual application scenarios, the above two-dimensional slope displacement vector cloud map is as follows: Figure 3 As shown in the figure, the above two-dimensional slope strain rate cloud diagram is as follows: Figure 4 As shown. Among them, for Figure 3 The two-dimensional slope displacement vector cloud map shown in the figure uses a gradient from blue to red to represent the areas with different probabilities of landslides in the weak layer of the target working slope. Among them, the blue area has the lowest probability of landslides, and the red area has the highest probability of landslides. Figure 4 The two-dimensional slope strain rate cloud diagram shown uses a gradient color from blue to green to represent the amount of landslide earthwork when landslide occurs in the weak layer of the target working slope.

[0082] S104, based on the geological parameters, the shear strength parameter range, the two-dimensional slope displacement vector cloud map and the two-dimensional slope strain rate cloud map, generate a coverage range of the predicted landslide area after the landslide occurs.

[0083] It should be noted that in actual application scenarios, since the occurrence of landslides is regional, that is, the landslide area is irregular. Therefore, this application generates a coverage range after the predicted landslide area occurs by configuring geological parameters, shear strength parameter range, two-dimensional slope displacement vector cloud map and two-dimensional slope strain rate cloud map, so that the generated coverage range is close to the irregular landslide coverage range when the landslide actually occurs, thereby improving the delineation accuracy of the landslide danger area.

[0084] S105. Determine the maximum safety distance based on the coverage range, and display the coverage range and the preset landslide prevention and control method corresponding to the applicable distance with the smallest difference from the maximum safety distance.

[0085] It should be noted that in actual application scenarios, since the probability of a landslide directly affects the maximum safety distance, and there are differences in landslide prevention and control methods for landslides of different degrees, for example: when the maximum safety distance is small, the probability of a landslide is small. At this time, it is only necessary to clear and strip the working side slope, and remove the waste earthwork accumulated on the working side slope to reduce the pressure of the working side slope on the soft layer to avoid the probability of a landslide. When the maximum safety distance is large, the probability of a landslide is large. At this time, in addition to clearing and stripping the working side slope, it is also necessary to simultaneously drain water to increase the strength of the soft layer, thereby reducing the probability of a landslide. Therefore, this application improves the configuration accuracy of landslide prevention and control measures for the target working side slope by configuring the corresponding relationship between the applicable distance and the preset landslide prevention and control method, and selecting the corresponding preset landslide prevention and control method for display based on the difference between the applicable distance and the maximum safety distance.

[0086] This application constructs a three-dimensional geological model of the target working side slope by configuring remote sensing sensing mapping data based on the target working side slope, and extracts geological parameters of the weak layer in the target working side slope from the three-dimensional geological model. The remote sensing mapping data can truly reflect the characteristics of the subsurface stratigraphic structure, improve the accuracy of determining the maximum safe distance, and thus improve the accuracy of delineating the landslide hazard area. Subsequently, the shear strength parameter range of the weak layer is obtained according to the geological parameters by configuring the preset limit equilibrium algorithm, and the parameter types of the shear strength parameter range include the cohesion range and the internal friction angle range. Based on the three-dimensional geological model and the shear strength parameter range, a two-dimensional slope displacement vector cloud map and a two-dimensional slope strain rate cloud map are generated. The two-dimensional slope displacement vector cloud map represents the predicted landslide area of ​​the target working side slope, and the two-dimensional slope displacement vector cloud map represents the landslide earthwork volume in the predicted landslide area, so as to realize the prediction of the area where the target working side slope may experience landslide, and also realize the prediction of the possible coverage range after the predicted landslide area occurs, so that this application improves the accuracy of delineating the landslide hazard area compared with the manual analysis method. Finally, by displaying the coverage range and the preset landslide prevention and control method corresponding to the applicable distance with the smallest difference from the maximum safety distance, by pre-configuring the corresponding preset landslide prevention and control method for different applicable distances, and when the maximum safety distance is obtained, the preset landslide prevention and control method corresponding to the applicable distance with the smallest difference from the obtained maximum safety distance is displayed, the generation time of the landslide prevention and control method is shortened while improving the configuration accuracy of the landslide prevention and control measures. It can be seen that this application improves the accuracy of determining the landslide danger area and the configuration accuracy of the landslide prevention and control measures.

[0087] In a possible implementation, the geological parameters of the weak layer in the target working slope are extracted from the three-dimensional geological model, including:

[0088] The strata at the junction of strongly weathered mudstone and strongly weathered sandstone in the three-dimensional geological model are determined as weak layers, and the boundary length of the weak layers is extracted;

[0089] The weak layer is divided based on the boundary length and the preset division unit length to obtain a plurality of measurement areas corresponding to each of the equal boundary lengths; the parameters of each measurement area in the three-dimensional geological model are extracted to obtain the geological parameters of the weak layer.

[0090] It should be noted that, in actual application scenarios, the value of the above-mentioned preset division unit length is set based on the soil information of the target working side slope. Since landslides have regional characteristics, that is, for a section of working side slope, since some areas of the slope are soaked by groundwater and the soil is relatively loose and has low strength, the probability of landslides in this area is higher than that in other areas of the working side slope. Therefore, the present application divides the weak layer based on the boundary length and the preset division unit length, obtains a plurality of measurement areas corresponding to each of the equal boundary lengths, and extracts parameters for each measurement area in the three-dimensional geological model, so that in the subsequent landslide simulation based on the obtained geological parameters of the weak layer, the influence of the soil information of different measurement areas on the probability of landslide occurrence is fully considered, thereby improving the accuracy of determining the coverage range.

[0091] It should be noted that, in actual application scenarios, the above-mentioned implementation method of extracting parameters from each measurement area in the three-dimensional geological model to obtain geological parameters of the weak layer may include the following steps A1 to A2.

[0092] Step A1, extract parameters of each measurement area in the three-dimensional geological model, obtain the initial geological parameter set of each measurement area, and set the identifier of the initial geological parameter set as the identifier of the corresponding measurement area, and trigger step A2.

[0093] Step A2: Summarize the initial geological parameter sets to obtain the geological parameters of the weak layer.

[0094] It should be noted that, in an actual application scenario, in order to facilitate the understanding of the above-mentioned measurement area, a possible implementation of the present application is specifically described here:

[0095] like Figure 5 The following is a schematic diagram of the measurement area. Figure 5 As shown in FIG. 1 , the weak layer is located inside the target working slope. After analyzing the soil information of the target working slope, a length with low soil strength is determined as the initial division length. Then, the initial division length with redundant length is determined as the above-mentioned preset division unit length, and the weak layer is divided to obtain the following: Figure 5 Test Area 1 and Test Area 2 shown.

[0096] In a possible implementation, the preset limit equilibrium algorithm includes: a first preset limit equilibrium algorithm and a second preset limit equilibrium algorithm, and the preset limit equilibrium algorithm is used to obtain the shear strength parameter range of the weak layer according to the geological parameters, and the parameter types of the shear strength parameter range include the cohesion range and the internal friction angle range, including:

[0097] For each measurement area: based on at least the initial stability coefficient and the angle of repose in the geological parameters in the measurement area, the cohesion of the measurement area is calculated using a first preset limit equilibrium algorithm; and based on at least the cohesion of the measurement area, the internal friction angle of the measurement area is calculated using a second preset limit equilibrium algorithm;

[0098] The cohesion and internal friction angle of each measured area are summarized to obtain the shear strength parameter range of the weak layer including the cohesion range and the internal friction angle range.

[0099] It should be noted that, in actual application scenarios, the above initial stability coefficient can be a coefficient determined based on the numerical value of the parameter with a high correlation with the probability of landslide in the geological parameters of the measured area. The larger the initial stability coefficient, the lower the probability of landslide in the measured area. The above initial stability coefficient can be determined by: assuming that the parameters with a high correlation with the probability of landslide in the geological parameters include: porosity, permeability and water saturation, setting the weight corresponding to porosity to 30%, setting the weight corresponding to permeability to 30%, setting the weight corresponding to water saturation to 40%, then calling the preset weighted mean algorithm, performing weighted average calculation based on the numerical values ​​of the above parameters and the corresponding weights, and determining the calculation result as the above initial stability coefficient.

[0100] It should be noted that, in actual application scenarios, the above-mentioned angle of repose refers to the minimum angle between the target working side slope and the horizontal plane when the target working side is in a critical transition state of sliding down the target working side slope in the measurement area.

[0101] In a possible implementation, based at least on the initial stability coefficient and the angle of repose among the geological parameters in the measurement area, the cohesion of the measurement area is calculated using a first preset limit equilibrium algorithm, including:

[0102] According to the initial stability coefficient F in the geological parameters of the measured area i i s and the angle of repose θ, by the formula of the first preset limit equilibrium algorithm:

[0103]

[0104] Calculate the cohesion C of the measured area i i , where L i is the preset division unit length, W i It is the weight of the soil strip in the geological parameters of the measurement area i.

[0105] In a possible implementation, calculating the internal friction angle of the measurement area based on the cohesion of the measurement area using a second preset limit equilibrium algorithm includes:

[0106]

[0107] Calculate the internal friction angle φ of the measurement area i i , where L i is the preset division unit length, C i is the cohesion, W i is the weight of the soil strip in the geological parameters of the measured area i, u i is the pore water pressure in the geological parameters of the measured area i, Q i is the horizontal force on the measured area i in the geological parameters, e i is the normal force on the measured area i in the geological parameters, R i is the sliding surface radius in the geological parameters of the measurement area i, θ is the repose angle in the geological parameters of the measurement area i, m θ is the calculation coefficient corresponding to the angle of repose in the geological parameters of the measurement area i, F i s It is the initial stability coefficient of the geological parameters in the measurement area i.

[0108] In a possible implementation, based on geological parameters, shear strength parameter range, two-dimensional slope displacement vector cloud map and two-dimensional slope strain rate cloud map, a coverage range of a predicted landslide area after a landslide occurs is generated, including:

[0109] The geological parameters, shear strength parameter range, two-dimensional slope displacement vector cloud map and two-dimensional slope strain rate cloud map are imported into the preset particle flow simulation model, so that the preset particle flow simulation model simulates the landslide process of the preset landslide area, and outputs the landslide coverage cloud map after the landslide occurs in the predicted landslide area when the simulation result meets the preset conditions; the coverage parameters are extracted from the landslide coverage cloud map.

[0110] It should be noted that, in actual application scenarios, the above-mentioned preset particle flow simulation model can be a numerical simulation software based on particle flow theory (Particle Flow Code in 3Dimensions, PFC3D). Since the landslide process can be regarded as the mechanical behavior under the interaction between soil particles, the preset particle flow simulation model is configured in this application based on geological parameters, shear strength parameter range, two-dimensional slope displacement vector cloud map and two-dimensional slope strain rate cloud map to simulate the landslide process in the preset landslide area, thereby obtaining a landslide coverage cloud map after the landslide is predicted in the landslide area, realizing the visual calibration of the landslide coverage cloud map, and improving the calibration accuracy.

[0111] It should be noted that in actual application scenarios, the above landslide coverage cloud map is as follows: Figure 6As shown, blue represents the coverage of silt, green represents the coverage of silty clay, and red represents the coverage of silty sand. Since the soil of the working slope usually includes silty clay, silty sand and silt, silty clay is not easy to flow due to its high viscosity, and will be squeezed to the upper part of the landslide body during the landslide, while silt has finer particles and better fluidity, so it will flow to the lower part of the landslide body during the landslide. The fluidity and particle size of silt are between silty clay and silty sand, so it is usually located in the middle of the landslide body. The present application configures a particle flow simulation model to output a cloud map of the landslide coverage after a landslide occurs in the landslide area, thereby realizing the coverage of each component of the landslide body, thereby realizing the deployment of corresponding preset landslide prevention and control methods for the coverage of landslide bodies of different soil types, thereby improving the configuration accuracy of landslide prevention and control measures.

[0112] In a possible implementation, assuming that the coverage parameters extracted from the landslide coverage cloud map include the coverage parameters corresponding to the three landslide types of silty clay, silty sand and silty soil, the above Figure 1 The implementation of step S105 shown may be:

[0113] For each landslide type:

[0114] Step B1, extracting the maximum safety distance from the parameters of the coverage area corresponding to the landslide type, and triggering step B2.

[0115] Step B2, respectively calculating the difference between the maximum safety distance obtained in step B1 and each applicable distance in the preset applicable distance list, and triggering step B3.

[0116] Step B3, querying the preset landslide prevention and control method database according to the applicable distance associated with the minimum difference in step B2, and caching the preset landslide prevention and control methods in the preset landslide prevention and control method database that have a corresponding relationship with the applicable distance.

[0117] Those skilled in the art can understand that, in actual application scenarios, the above-mentioned preset conditions are that when the landslide process simulated by the above-mentioned preset particle flow simulation model reaches a stable state, the landslide body parameters in the simulation process are substituted into the physical equation of the force and displacement law and the motion equation of Newton's second law, respectively, so that the physical equation of the force and displacement law and the motion equation of Newton's second law are identical. This application does not make too many restrictions and redundant descriptions on the substitution process of the above-mentioned preset conditions.

[0118] In a possible implementation, before the preset particle flow simulation model simulates the landslide process in the preset landslide area, the method further includes:

[0119] Import the blasting vibration influence coefficient into the preset particle flow simulation model.

[0120] It should be noted that in actual application scenarios, since mining requires blasting, the vibration generated by blasting will also affect the landslide process. Therefore, this application imports the blasting vibration influence coefficient into the preset particle flow simulation model before simulating the landslide process in the preset landslide area, so that the preset particle flow simulation model fully considers the influence of the blasting vibration wave on the landslide process during the landslide simulation process, thereby improving the accuracy of determining the coverage range.

[0121] The second aspect of the present application provides a working slope landslide prevention and control system, such as Figure 7 As shown in the figure, the working slope landslide prevention and control system includes:

[0122] Parameter extraction module 701, used to construct a three-dimensional geological model of the target working slope based on the remote sensing mapping data of the target working slope, and extract geological parameters of the weak layer in the target working slope from the three-dimensional geological model;

[0123] A parameter range extraction module 702 is used to obtain a shear strength parameter range of a weak layer according to geological parameters using a preset limit equilibrium algorithm, wherein the parameter types of the shear strength parameter range include a cohesion range and an internal friction angle range;

[0124] A cloud map generation module 703 is used to generate a two-dimensional slope displacement vector cloud map and a two-dimensional slope strain rate cloud map based on a three-dimensional geological model and a shear strength parameter range using a preset cloud map generation software, wherein the two-dimensional slope displacement vector cloud map represents a predicted landslide area of ​​a target working slope, and the two-dimensional slope displacement vector cloud map represents a landslide earthwork volume in the predicted landslide area;

[0125] The coverage range determination module 704 is used to generate the coverage range of the predicted landslide area after the landslide occurs based on the geological parameters, the shear strength parameter range, the two-dimensional slope displacement vector cloud map and the two-dimensional slope strain rate cloud map;

[0126] The data display module 705 is used to determine the maximum safety distance based on the coverage range, and to display the coverage range and the preset landslide prevention and control method corresponding to the applicable distance with the smallest difference from the maximum safety distance.

[0127] In a possible implementation, the parameter extraction module 701 is set to extract the geological parameters of the weak layer in the target working slope from the three-dimensional geological model as follows:

[0128] The strata at the junction of strongly weathered mudstone and strongly weathered sandstone in the three-dimensional geological model are determined as weak layers, and the boundary length of the weak layers is extracted;

[0129] The weak layer is divided based on the boundary length and the preset division unit length to obtain a plurality of measurement areas corresponding to each of the equal boundary lengths; the parameters of each measurement area in the three-dimensional geological model are extracted to obtain the geological parameters of the weak layer.

[0130] In a possible implementation, the parameter range extraction module 702 is set to: when the parameter type of the shear strength parameter range includes the cohesion range and the internal friction angle range, the shear strength parameter range of the weak layer is obtained according to the geological parameters using the preset limit equilibrium algorithm:

[0131] For each measurement area: based at least on the initial stability coefficient and the angle of repose in the geological parameters in the measurement area, the cohesion of the measurement area is calculated using a first preset limit equilibrium algorithm; based at least on the cohesion of the measurement area, the internal friction angle of the measurement area is calculated using a second preset limit equilibrium algorithm, the preset limit equilibrium algorithm including a first preset limit equilibrium algorithm and a second preset limit equilibrium algorithm;

[0132] The cohesion and internal friction angle of each measured area are summarized to obtain the shear strength parameter range of the weak layer including the cohesion range and the internal friction angle range.

[0133] In a possible implementation, the parameter range extraction module 702 is set to: when calculating the cohesion of the measurement area using the first preset limit equilibrium algorithm based on at least the initial stability coefficient and the angle of repose among the geological parameters in the measurement area:

[0134] According to the initial stability coefficient F in the geological parameters of the measured area i i s and the angle of repose θ, by the formula of the first preset limit equilibrium algorithm:

[0135]

[0136] Calculate the cohesion C of the measured area i i , where L i is the preset division unit length, W i It is the weight of the soil strip in the geological parameters of the measurement area i.

[0137] In a possible implementation, the parameter range extraction module 702 is set to:

[0138]

[0139] Calculate the internal friction angle φ of the measurement area i i , where L i is the preset division unit length, C iis the cohesion, W i is the weight of the soil strip in the geological parameters of the measured area i, u i is the pore water pressure in the geological parameters of the measured area i, Q i is the horizontal force on the measured area i in the geological parameters, e i is the normal force on the measured area i in the geological parameters, R i is the sliding surface radius in the geological parameters of the measurement area i, θ is the repose angle in the geological parameters of the measurement area i, m θ is the calculation coefficient corresponding to the angle of repose in the geological parameters of the measurement area i, F i s It is the initial stability coefficient of the geological parameters in the measurement area i.

[0140] In a possible implementation, the coverage range determination module is set to:

[0141] The geological parameters, shear strength parameter range, two-dimensional slope displacement vector cloud map and two-dimensional slope strain rate cloud map are imported into the preset particle flow simulation model, so that the preset particle flow simulation model simulates the landslide process of the preset landslide area, and outputs the landslide coverage cloud map after the landslide occurs in the predicted landslide area when the simulation result meets the preset conditions; the coverage parameters are extracted from the landslide coverage cloud map.

[0142] In a possible implementation, the coverage range determining module 704 is further configured to:

[0143] Before the preset particle flow simulation model simulates the landslide process in the preset landslide area, the blasting vibration influence coefficient is imported into the preset particle flow simulation model.

[0144] A third aspect of the present application provides an electronic device, comprising at least one processor and a memory connected to the processor, wherein:

[0145] The memory is used to store computer programs;

[0146] The processor is used to execute a computer program so that the electronic device can implement the slope landslide prevention and control method of the first aspect or any implementation method of the first aspect.

[0147] This application Figure 8 A third aspect provides a structural schematic diagram of an electronic device as shown in Figure 8The electronic devices in the embodiments of the present application may include, but are not limited to, fixed terminals such as mobile phones, notebook computers, PDAs (personal digital assistants), PADs (tablet computers), desktop computers, and the like. Figure 8 The electronic device shown is merely an example and should not bring any limitation to the functions and scope of use of the embodiments of the present application.

[0148] like Figure 8 As shown, the electronic device may include a processing device (e.g., a central processing unit, a graphics processing unit, etc.) 801, which can perform various appropriate actions and processes according to a program stored in a read-only memory (ROM) 802 or a program loaded from a storage device 808 to a random access memory (RAM) 803. When the electronic device is powered on, various programs and data required for the operation of the electronic device are also stored in the RAM 803. The processing device 801, the ROM 802, and the RAM 803 are connected to each other via a bus 804. An input / output (I / O) interface 805 is also connected to the bus 804.

[0149] Typically, the following devices may be connected to the I / O interface 805: input devices 806 including, for example, a touch screen, a touch pad, a keyboard, a mouse, a camera, a microphone, an accelerometer, a gyroscope, etc.; output devices 807 including, for example, a liquid crystal display (LCD), a speaker, a vibrator, etc.; storage devices 808 including, for example, a memory card, a hard disk, etc.; and communication devices 809. The communication device 809 may allow the electronic device to communicate with other devices wirelessly or by wire to exchange data. Although Figure 8 An electronic device having various devices is shown, but it should be understood that it is not required to implement or possess all the devices shown. More or fewer devices may be implemented or possessed instead.

[0150] Also provided in an embodiment of the present application is a computer program product, including computer-readable instructions. When the computer-readable instructions are executed on an electronic device, the electronic device implements any one of the working slope landslide prevention and control methods provided in the embodiments of the present application.

[0151] The fourth aspect of the present application also provides a computer-readable storage medium, which carries one or more computer programs. When one or more computer programs are executed by an electronic device, the electronic device can implement any one of the working slope landslide prevention and control methods provided in the embodiments of the present application.

[0152] It should also be noted that the device embodiments described above are merely schematic, wherein the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed over multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the scheme of this embodiment. In addition, in the drawings of the device embodiments provided by the present application, the connection relationship between the modules indicates that there is a communication connection between them, which may be specifically implemented as one or more communication buses or signal lines.

[0153] Through the description of the above implementation mode, the technicians in the field can clearly understand that the present application can be implemented by means of software plus necessary general hardware, and of course, it can also be implemented by special hardware including special integrated circuits, special CPUs, special memories, special components, etc. In general, all functions completed by computer programs can be easily implemented by corresponding hardware, and the specific hardware structure used to implement the same function can also be various, such as analog circuits, digital circuits or special circuits. However, for the present application, software program implementation is a better implementation mode in more cases. Based on such an understanding, the technical solution of the present application is essentially or the part that contributes to the prior art can be embodied in the form of a software product, which is stored in a readable storage medium, such as a computer floppy disk, a U disk, a mobile hard disk, a ROM, a RAM, a disk or an optical disk, etc., including a number of instructions to enable a computer device (which can be a personal computer, a training device, or a network device, etc.) to execute the methods described in each embodiment of the present application.

[0154] In the above embodiments, all or part of the embodiments may be implemented by software, hardware, firmware or any combination thereof. When implemented by software, all or part of the embodiments may be implemented in the form of a computer program product.

[0155] The computer program product includes one or more computer instructions. When the computer program instructions are loaded and executed on a computer, the process or function described in the embodiment of the present application is generated in whole or in part. The computer may be a general-purpose computer, a special-purpose computer, a computer network, or other programmable devices. The computer instructions may be stored in a computer-readable storage medium, or transmitted from one computer-readable storage medium to another computer-readable storage medium. For example, the computer instructions may be transmitted from a website site, a computer, a training device, or a data center by wired (e.g., coaxial cable, optical fiber, digital subscriber line (DSL)) or wireless (e.g., infrared, wireless, microwave, etc.) mode to another website site, computer, training device, or data center. The computer-readable storage medium may be any available medium that a computer can store or a data storage device such as a training device, a data center, etc. that includes one or more available media integrations. The available medium may be a magnetic medium, (e.g., a floppy disk, a hard disk, a tape), an optical medium (e.g., a DVD), or a semiconductor medium (e.g., a solid-state drive (SSD)), etc.

Claims

1. A method for preventing and controlling landslides on working slopes, characterized in that: include: Constructing a three-dimensional geological model of the target working slope based on remote sensing mapping data of the target working slope, and extracting geological parameters of the weak layer in the target working slope from the three-dimensional geological model; Obtaining a shear strength parameter range of the weak layer according to the geological parameters using a preset limit equilibrium algorithm, wherein the parameter types of the shear strength parameter range include a cohesion range and an internal friction angle range; Using a preset cloud map generation software to generate a two-dimensional slope displacement vector cloud map and a two-dimensional slope strain rate cloud map based on the three-dimensional geological model and the shear strength parameter range, the two-dimensional slope displacement vector cloud map represents the predicted landslide area of ​​the target working slope, and the two-dimensional slope displacement vector cloud map represents the landslide earthwork volume of the predicted landslide area; Based on the geological parameters, the shear strength parameter range, the two-dimensional slope displacement vector cloud map and the two-dimensional slope strain rate cloud map, generating a coverage range after the landslide occurs in the predicted landslide area; The maximum safety distance is determined based on the coverage range, and the coverage range and the preset landslide prevention and control method corresponding to the applicable distance with the smallest difference from the maximum safety distance are displayed.

2. The method for preventing and controlling landslides on working slopes according to claim 1, characterized in that: The step of extracting geological parameters of the weak layer in the target working slope from the three-dimensional geological model includes: Determine the stratum at the junction of the strongly weathered mudstone and the strongly weathered sandstone in the three-dimensional geological model as the weak layer, and extract the boundary length of the weak layer; The weak layer is divided based on the boundary length and the preset division unit length to obtain a plurality of measurement areas corresponding to respective equal boundary lengths; parameters are extracted from each of the measurement areas in the three-dimensional geological model to obtain geological parameters of the weak layer.

3. The method for preventing and controlling landslides on working slopes according to claim 2, characterized in that: The preset limit equilibrium algorithm includes: a first preset limit equilibrium algorithm and a second preset limit equilibrium algorithm, wherein the preset limit equilibrium algorithm is used to obtain the shear strength parameter range of the weak layer according to the geological parameters, and the parameter types of the shear strength parameter range include the cohesion range and the internal friction angle range, including: For each of the measurement areas: based at least on the initial stability coefficient and the angle of repose in the geological parameters in the measurement area, the cohesion of the measurement area is calculated using the first preset limit equilibrium algorithm; and the internal friction angle of the measurement area is calculated based at least on the cohesion of the measurement area using the second preset limit equilibrium algorithm; The cohesion and the internal friction angle of each of the measurement areas are summarized to obtain the shear strength parameter range of the weak layer including the cohesion range and the internal friction angle range.

4. The method for preventing and controlling landslides on working slopes according to claim 3 is characterized in that: The step of calculating the cohesion of the measurement area by using the first preset limit equilibrium algorithm based on at least the initial stability coefficient and the angle of repose among the geological parameters in the measurement area comprises: According to the initial stability coefficient F in the geological parameters of the measurement area i i s and the repose angle θ, by the formula of the first preset limit equilibrium algorithm: The cohesion C of the measurement area i is obtained i , wherein the L i is the preset division unit length, the W i It is the weight of the soil strip in the geological parameters of the measurement area i.

5. The method for preventing and controlling landslides on working slopes according to claim 3 is characterized in that: The calculating the internal friction angle of the measurement area based on the cohesion of the measurement area by using the second preset limit equilibrium algorithm includes: Obtain the internal friction angle φ of the measurement area i i , wherein the L i is the preset division unit length, the C i is the cohesive force, the W i is the weight of the soil strip in the geological parameters of the measurement area i, and u i is the pore water pressure in the geological parameters of the measurement area i, and Q i is the horizontal force on the measured area i in the geological parameters, and the e i is the normal force on the measured area i in the geological parameters, R i is the sliding surface radius in the geological parameters of the measurement area i, θ is the repose angle in the geological parameters of the measurement area i, and m θ is the calculation coefficient corresponding to the angle of repose in the geological parameters of the measurement area i, and the F i s It is the initial stability coefficient among the geological parameters of the measurement area i.

6. The method for preventing and controlling landslides on working slopes according to claim 1, characterized in that: The generating of the coverage range of the predicted landslide area after the landslide occurs based on the geological parameters, the shear strength parameter range, the two-dimensional slope displacement vector cloud map and the two-dimensional slope strain rate cloud map comprises: The geological parameters, the shear strength parameter range, the two-dimensional slope displacement vector cloud map and the two-dimensional slope strain rate cloud map are imported into a preset particle flow simulation model, so that the preset particle flow simulation model simulates the landslide process of the preset landslide area, and outputs a landslide coverage cloud map after the landslide occurs in the predicted landslide area when the simulation result meets the preset conditions; and extracts the parameters of the coverage range from the landslide coverage cloud map.

7. The method for preventing and controlling landslides on working slopes according to claim 6, characterized in that: Before the preset particle flow simulation model simulates the landslide process in the preset landslide area, the method further includes: The blasting vibration influence coefficient is introduced into the preset particle flow simulation model.

8. A working slope landslide prevention and control system, characterized in that: include: A parameter extraction module is used to construct a three-dimensional geological model of the target working slope based on the remote sensing mapping data of the target working slope, and to extract geological parameters of the weak layer in the target working slope from the three-dimensional geological model; A parameter range extraction module, used to obtain the shear strength parameter range of the weak layer according to the geological parameters using a preset limit equilibrium algorithm, wherein the parameter types of the shear strength parameter range include a cohesion range and an internal friction angle range; A cloud map generation module, used for generating a two-dimensional slope displacement vector cloud map and a two-dimensional slope strain rate cloud map based on the three-dimensional geological model and the shear strength parameter range using a preset cloud map generation software, wherein the two-dimensional slope displacement vector cloud map represents the predicted landslide area of ​​the target working slope, and the two-dimensional slope displacement vector cloud map represents the landslide earthwork volume of the predicted landslide area; A coverage range determination module, used for generating a coverage range after a landslide occurs in the predicted landslide area based on the geological parameters, the shear strength parameter range, the two-dimensional slope displacement vector cloud map and the two-dimensional slope strain rate cloud map; The data display module is used to determine the maximum safety distance based on the coverage range, and to display the coverage range and the preset landslide prevention and control method corresponding to the applicable distance with the smallest difference from the maximum safety distance.

9. An electronic device, characterized in that: The method comprises at least one processor and a memory connected to the processor, wherein: The memory is used to store computer programs; The processor is used to execute the computer program so that the electronic device can implement the working slope landslide prevention and control method as described in any one of claims 1 to 7.

10. A computer storage medium, characterized in that: The storage medium carries one or more computer programs, and when the one or more computer programs are executed by an electronic device, the electronic device can implement the working slope landslide prevention and control method as described in any one of claims 1 to 7.

Images