Method, device and equipment for determining safety coefficient range of karst cave slope and medium

Abstract

The invention relates to a method, device and equipment for determining the safety coefficient range of a karst cave slope and a medium. The method comprises the steps that the maximum value and the minimum value of the slope in the X direction and the Y direction are obtained; based on the maximum value and the minimum value of the slope in the X direction, the slope is divided in the X direction according to the vertical division interval, and a plurality of vertical units are obtained; based on the maximum value and the minimum value of the slope in the Y direction, unit bodies with random heights are generated in all the vertical units; according to the unit bodies, determining the distribution condition of the rock-soil body and the karst cave in each vertical unit, and constructing a slope karst cave simulation model; performing assignment calculation on the slope karst cave simulation model to obtain a slope safety coefficient; and a plurality of groups of slope safety coefficients are obtained through iterative calculation, and the range of the slope safety coefficients is determined based on each group of slope safety coefficients. Reliable data support can be provided for engineering design of the karst cave slope, and safety and engineering economy are both considered.

Description

Technical Field

[0001] The invention relates to the technical field of slope stability analysis, and in particular to a method, device, equipment and medium for determining the safety factor range of a karst cave slope. Background Art

[0002] The karst slope refers to the slope structure formed in the karst area. Usually, due to the particularity of karst landform, the stability of the slope is low. The karst cave is formed under the long-term action of groundwater, with irregular cavities and complex internal structures, which makes the mechanical properties of the slope difficult to predict.

[0003] Traditional methods for determining the safety factor of cave slopes are usually based on calculation methods such as limit equilibrium theory and finite element analysis. However, these methods are difficult to fully consider the randomness and irregularity of caves. Especially when the caves are large in scale or close to the sliding surface, the impact of caves on the slopes is particularly significant. Ignoring this will lead to an overestimation of the slope safety factor, thus posing a hidden danger to engineering safety.

[0004] Therefore, in order to improve the safety of cave slopes, excessive reinforcement measures are usually taken in engineering design, such as significantly increasing the number of anti-slide piles, increasing the density of anchor rods, or conducting extensive grouting reinforcement to improve the safety factor of the slope. However, although this excessive reinforcement strategy can improve the stability of the slope to a certain extent, it is often accompanied by high engineering costs and waste of resources. At the same time, since the reinforcement measures are not targeted at the specific morphology and distribution characteristics of the caves, they may not be able to effectively solve the hidden dangers of the caves to the slope structure, resulting in limited effect on improving the safety factor. In addition, excessive reinforcement may also change the original mechanical balance of the slope, increase local stress concentration on the slope, or cause unnecessary environmental damage. Summary of the invention

[0005] Based on this, it is necessary to provide a method, device, equipment and medium for determining the safety factor range of the cave slope that can calculate the safety factor range of the cave slope, provide reliable data support for engineering design, and take into account both safety and engineering economy in response to the above technical problems.

[0006] A method for determining a safety factor range of a karst cave slope, the method comprising: Get the maximum and minimum values ​​of the slope in the X and Y directions; Setting a vertical division interval, based on the maximum and minimum values ​​of the slope in the X direction, dividing the slope along the X direction according to the vertical division interval to obtain a plurality of vertical units; Based on the maximum and minimum values ​​of the slope in the Y direction, a unit body with a random height is generated in each vertical unit; and the distribution of the rock and soil mass and the karst cave in each vertical unit is determined according to the unit body; A slope cave simulation model is constructed based on the distribution of rock and soil bodies and caves in each vertical unit; Performing value calculation on the slope cave simulation model to obtain the slope safety factor; Several groups of slope safety factors are obtained through iterative calculation, and the range of slope safety factors is determined based on each group of slope safety factors.

[0007] A device for determining a safety factor range of a karst cave slope, the device comprising: The initial data acquisition module is used to obtain the maximum and minimum values ​​of the slope in the X and Y directions; A vertical unit division module is used to set a vertical division interval, and based on the maximum and minimum values ​​of the slope in the X direction, divide the slope along the X direction according to the vertical division interval to obtain a number of vertical units; The geological distribution random generation module is used to generate a unit body of random height in each vertical unit based on the maximum and minimum values ​​of the slope in the Y direction; and determine the distribution of rock and soil bodies and karst caves in each vertical unit according to the unit body; A simulation modeling module is used to construct a slope cave simulation model based on the distribution of rock and soil bodies and caves in each vertical unit; A slope safety factor calculation module is used to perform value calculation on the slope cave simulation model to obtain the slope safety factor; The slope safety factor range determination module is used to obtain several groups of slope safety factors through iterative calculation, and determine the range of the slope safety factor based on each group of slope safety factors.

[0008] A computer device comprises a memory and a processor, wherein the memory stores a computer program, and the processor implements the steps of the method for determining the range of the safety factor of a karst cave slope when executing the computer program.

[0009] A computer-readable storage medium stores a computer program, which, when executed by a processor, implements the steps of the method for determining the range of the safety factor of a karst cave slope.

[0010] The method, device, equipment and medium for determining the range of the safety factor of the karst cave slope obtain the maximum and minimum values ​​of the slope in the X direction and the Y direction; set a vertical division interval, and divide the slope along the X direction according to the vertical division interval based on the maximum and minimum values ​​of the slope in the X direction to obtain a number of vertical units; generate a unit body of random height in each vertical unit based on the maximum and minimum values ​​of the slope in the Y direction; determine the distribution of rock and soil bodies and karst caves in each vertical unit according to the unit body; construct a slope karst cave simulation model according to the distribution of rock and soil bodies and karst caves in each vertical unit; assign values ​​to the slope karst cave simulation model to obtain the slope safety factor; obtain a number of groups of slope safety factors through iterative calculation, and determine the range of the slope safety factor based on each group of slope safety factors.

[0011] The present invention discretizes the slope into a number of vertical units through reasonable vertical division intervals, so that it is closer to the geological distribution during vertical drilling during exploration. By generating random heights of rock and soil bodies and caves in vertical units, the complex geological structure of rock and soil bodies and caves can be effectively simulated, so that the obtained simulation results are closer to the actual situation, thereby obtaining a more accurate slope safety factor. In addition, since the height of the generated cave has a certain randomness, through multiple iterative calculations, the one-sidedness of a single result is avoided, and the deviation caused by uncertainty in the calculation is reduced, thereby solving a more reliable range of slope safety factors, providing reliable data support for subsequent evaluation of the stability of the edge skin and engineering design optimization, and avoiding over-design and saving engineering costs. BRIEF DESCRIPTION OF THE DRAWINGS

[0012] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on the structures shown in these drawings without paying creative work.

[0013] Figure 1 A schematic flow chart of a method for determining a safety factor range of a karst cave slope in one embodiment; Figure 2 is a schematic diagram of the outer contour of a slope in one embodiment; Figure 3 is a schematic diagram of dividing a slope into a plurality of vertical units along the X direction in one embodiment; Figure 4 It is a schematic diagram of the distribution of rock and soil bodies and caves in a vertical unit generated along the Y direction in one embodiment; Figure 5 An initial cave model generated in one embodiment; Figure 6 A schematic diagram of a slope outer contour model constructed by a FLAC3D numerical model in one embodiment; Figure 7 A schematic diagram of a slope cave simulation model constructed by a FLAC3D numerical model in one embodiment; Figure 8 It is a structural block diagram of a device for determining the safety factor range of a karst cave slope in one embodiment; Fig. 9 FIG. 4 is a diagram showing the internal structure of a computer device in one embodiment.

[0014] The realization of the purpose, functional features and advantages of the present invention will be further explained in conjunction with embodiments and with reference to the accompanying drawings. DETAILED DESCRIPTION

[0015] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0016] It should be noted that all directional indications (such as up, down, left, right, front, back, etc.) in the embodiments of the present invention are only used to explain the relative position relationship, movement status, etc. between the components under a certain specific posture (as shown in the accompanying drawings). If the specific posture changes, the directional indication will also change accordingly.

[0017] In addition, in the present invention, descriptions such as "first", "second", etc. are only used for descriptive purposes and cannot be understood as indicating or implying their relative importance or implicitly indicating the number of the indicated technical features. Therefore, the features defined as "first" or "second" may explicitly or implicitly include at least one of the features. In the description of the present invention, the meaning of "plurality" is at least two, such as two, three, etc., unless otherwise clearly and specifically defined.

[0018] It can be understood that the technical solutions between the various embodiments of the present invention can be combined with each other, but it must be based on the fact that ordinary technicians in the field can implement it. When the combination of technical solutions is contradictory or cannot be implemented, it should be deemed that such combination of technical solutions does not exist and is not within the scope of protection required by the present invention.

[0019] The following will describe the implementation of the present invention in detail with reference to the accompanying drawings in the embodiment diagram of the present invention.

[0020] Example 1 The present embodiment discloses a method for determining the safety factor range of a cave slope. By reasonable vertical division intervals, the slope is discretized into a number of vertical units, so that it is closer to the geological distribution when vertical drilling is performed during exploration. However, since the geological exploration is a vertical point drilling method, the information about the height and distribution of the cave in the exploration data is obtained by statistics, so there is also a certain uncertainty in the height of the cave. The present invention can effectively simulate the complex geological structure of the rock and soil body and the cave by generating random heights of the rock and soil body and the cave in the vertical unit, so that the obtained simulation results are closer to the actual situation, and then a more accurate slope safety factor is obtained. In addition, since the generated cave height has a certain randomness, through multiple iterative calculations, the one-sidedness of a single result is avoided, and the deviation caused by uncertainty in the calculation is reduced, so as to solve a more reliable slope safety factor range, for the subsequent evaluation of slope stability, to provide reliable data support for the engineering design optimization of the cave slope, can also avoid over-design, save engineering costs.

[0021] like Figure 1 As shown, a method for determining the safety factor range of a cave slope is provided, comprising the following steps: Step 201, obtaining the maximum and minimum values ​​of the slope in the X direction and the Y direction.

[0022] Step 202, setting a vertical division interval, based on the maximum and minimum values ​​of the slope in the X direction, dividing the slope along the X direction according to the vertical division interval to obtain a plurality of vertical units.

[0023] Step 203, based on the maximum and minimum values ​​of the slope in the Y direction, generate a unit body with a random height in each vertical unit; and determine the distribution of the rock and soil mass and the cave in each vertical unit according to the unit body.

[0024] Step 204, constructing a slope cave simulation model according to the distribution of rock and soil bodies and caves in each vertical unit.

[0025] Step 205, performing value assignment calculation on the slope cave simulation model to obtain the slope safety factor.

[0026] Step 206, obtaining several groups of slope safety factors through iterative calculation, and determining the range of the slope safety factor based on each group of slope safety factors.

[0027] In the specific implementation process of step 201, firstly, the outer contour information of the slope is obtained through geological exploration data, and the coordinate values ​​of each corner point of the slope to be calculated are determined. Figure 2 , is the outer contour shape of the slope and the coordinate values ​​of each corner point of the slope in one case provided in this embodiment. According to the coordinate values ​​of each corner point of the slope, the minimum value X of the slope in the X direction is determined. min With the maximum value Xmax , and determine the minimum slope value Y in the Y direction min With the maximum value Y max .

[0028] In the specific implementation process of step 202, according to the information such as the location of the drilling point for survey and measurement, the vertical division interval is set, and then according to the minimum value X of the slope in the X direction min With the maximum value X max , divide the slope into multiple long vertical units along the X direction according to the vertical division interval, such as Figure 3 As shown in FIG. 1 , the vertical units are divided when the vertical division interval is set to 1 m. When the vertical units are divided, the number of vertical units n is Ceil(X max -X min ), Ceil means round up.

[0029] It can be understood that the distribution of caves in the exploration data is obtained through statistics after obtaining local geological information from each borehole, so there is uncertainty. Based on this, a reasonable vertical division interval is set according to the distribution position of the drilling points, and the slope is discretized into several vertical units according to the vertical division interval. In the subsequent calculation, it is equivalent to only considering the height distribution of the caves in the vertical unit, which simplifies the calculation process. In addition, through such division, the uncertainty is confined to each vertical unit without having to be accumulated in the entire model, which improves the accuracy of the cave estimation and obtains a more accurate slope safety factor.

[0030] In the specific implementation process of step 203, first, the distribution of the rock and soil mass and the karst caves of the entire slope is obtained based on the exploration data; based on the distribution of the rock and soil mass and the karst caves of the entire slope, the initial distribution form of the rock and soil mass and the karst caves of the slope is determined, wherein the rock and soil mass is uniformly distributed and the karst caves are normally distributed; based on the uniform distribution form of the rock and soil mass of the slope and the normal distribution form of the karst caves, the initial distribution information of the rock and soil mass and the karst caves in each vertical unit is determined.

[0031] Then, the vertical unit is evenly partitioned, any area is randomly selected, and the geological body type is selected according to the initial distribution information of the rock and soil body and the cave of the vertical unit; if the geological body type is rock and soil body, the first unit body representing the height of the rock and soil body is randomly generated according to the initial distribution information of the rock and soil body, and the relevant data is saved; when the data of the first unit body is saved, the saving format is: (X temp , X temp +1, Y _initial , Y _initial +H rock , Rock).

[0032] If the geological body type is a cave, a second unit body representing the height of the cave is randomly generated according to the initial distribution information of the cave, and the relevant data is saved; when the data of the second unit body is saved, the saving format is: (X temp , X temp +1, Y _initial , Y _initial +H karst , Karst).

[0033] Among them, X temp Indicates the current vertical unit coordinate x value; X temp +1 means the vertical unit coordinate x value is increased by 1; Y _initial Indicates the initial value along the Y direction; Y _initial +H rock Represents the y coordinate value of the first unit cell; Y _initial +H karst Represents the y coordinate value of the second unit body; Rock and Karst represent the geological body type.

[0034] Finally, the distribution of rock and soil bodies and caves in each vertical unit is determined based on the first unit body and the second unit body.

[0035] In addition, when generating a unit body of random height, a maximum height value and a minimum height value are preset; and a unit body of random height is generated based on the maximum height value and the minimum height value.

[0036] Specifically, it is assumed that the preset minimum height is 1 m and the maximum height is 4 m. First, initialize the initial value Y in the Y direction _initial =Y min .

[0037] Then the vertical units are evenly partitioned, any area is randomly selected, and the geological body type (rock body or cave) is selected according to the initial distribution information of the rock body and cave of the vertical unit, where the probability of rock body and cave being selected is 0.5.

[0038] If the selected regional geological body type is rock and soil body, a random height H is generated between the minimum height of 1m and the maximum height of 4m. rock The first unit body of the rock and soil model data is saved at the same time, and the saving format is (X temp , X temp +1, Y _initial , Y _initial +H rock , Rock). If the selected regional geological body type is a cave, a random height of H is generated between the minimum height of 1m and the maximum height of 4m. karst The second unit of the cave data is saved at the same time, and the saving format is (Xtemp , X temp +1, Y _initial , Y _initial +H karst , Karst).

[0039] When Y _initial +H rock or Y _initial +H karst Greater than or equal to Y max When , it means that the distribution of rock and soil and karst caves in the current vertical unit has been analyzed, resulting in the following Figure 4 The initial cave model of a long strip is shown. When the randomly selected area has a Y value max When saving data, modify it to (X temp , X temp +1, Y _initial , Y max , Rock) or (X temp , X temp +1, Y _initial , Y max , Karst), followed by Y _initial The update is: Y _initial =Y _initial +H rock or Y _initial =Y _initial +H karst .

[0040] Repeat the above operation until the distribution of rock and soil and caves in all vertical units is analyzed and the following is obtained: Figure 5 The initial cave model is shown.

[0041] In the specific implementation process of step 204, firstly according to Figure 2 The slope outer contour in the FLAC3D numerical model is used to generate a complete slope outer contour model, such as Figure 6 As shown in the figure, the model consists of rock and soil without any cave areas.

[0042] Then, based on the second unit volume data stored in step 203, Figure 6 Obtain the position coordinates (x, y) of the corresponding unit body from the slope outer contour model in the figure; when the corresponding unit body in the slope outer contour model is the second unit body, delete the unit body, and finally generate the following Figure 7 The simulation model of the slope cave is shown.

[0043] In the specific implementation process of step 205, the slope cave simulation model is assigned with values ​​and calculated to obtain the slope safety factor. The mechanical parameters of the rock and soil mass are shown in Table 1.

[0044] Table 1 Mechanical parameters of rock and soil

[0045] In the specific implementation process of step 206, since the height of the cave in the vertical unit is generated randomly, the slope safety factor calculated by a slope cave simulation model has a certain one-sidedness. Therefore, several groups of slope safety factors are obtained through iterative calculation, and then the average value and standard deviation of each group of slope safety factors are solved. The range of the slope safety factor is calculated based on the average value and standard deviation. The calculation expression is: ; In the formula, Indicates the minimum value of the slope safety factor range; Indicates the maximum value of the slope safety factor range; represents the average value; Represents standard deviation.

[0046] Although this embodiment Figure 1 The steps in the process are shown in sequence as indicated by the arrows, but the steps are not necessarily executed in the order indicated by the arrows. Unless otherwise specified in this document, there is no strict order restriction for the execution of the steps, and the steps can be executed in other orders. Figure 1 At least part of the steps may include multiple sub-steps or multiple stages. These sub-steps or stages are not necessarily executed at the same time, but can be executed at different times. The execution order of these sub-steps or stages is not necessarily sequential, but can be executed in turn or alternately with other steps or at least part of the sub-steps or stages of other steps.

[0047] Example 2 Based on the method for determining the safety factor range of a karst cave slope in Example 1, this embodiment discloses a device for determining the safety factor range of a karst cave slope. As shown in FIG8 , the device for determining the safety factor range of a karst cave slope comprises: an initial data acquisition module 401, a vertical unit division module 402, a geological distribution random generation module 403, a simulation modeling module 404, a slope safety factor calculation module 405, and a slope safety factor range determination module 406, wherein: The initial data acquisition module 401 is used to obtain the maximum and minimum values ​​of the slope in the X direction and the Y direction.

[0048] The vertical unit division module 402 is used to set the vertical division interval, and based on the maximum and minimum values ​​of the slope in the X direction, divide the slope along the X direction according to the vertical division interval to obtain a number of vertical units.

[0049] The geological distribution random generation module 403 is used to generate a unit body of random height in each vertical unit based on the maximum and minimum values ​​of the slope in the Y direction; and determine the distribution of rock and soil bodies and caves in each vertical unit according to the unit body.

[0050] The simulation modeling module 404 is used to construct a slope cave simulation model according to the distribution of rock and soil bodies and caves in each vertical unit.

[0051] The slope safety factor calculation module 405 is used to perform value assignment calculation on the slope cave simulation model to obtain the slope safety factor.

[0052] The slope safety factor range determination module 406 is used to obtain a plurality of groups of slope safety factors through iterative calculation, and determine the range of the slope safety factor based on each group of slope safety factors.

[0053] In this embodiment, the specific working process and working principle of the initial data acquisition module 401, the vertical unit division module 402, the geological distribution random generation module 403, the simulation modeling module 404, the slope safety factor calculation module 405 and the slope safety factor range determination module 406 are the same as those in the method of embodiment 1, so they are not repeated in this embodiment. Each unit module can be implemented in whole or in part by software, hardware and a combination thereof, and each unit module can be embedded in or independent of a processor in a computer device in the form of hardware, or can be stored in a memory in a computer device in the form of software, so that the processor can call and execute the operations corresponding to each of the above unit modules.

[0054] Example 3 like Fig. 9 The terminal device disclosed in this embodiment includes a transmitter, a receiver, a memory, and a processor. The transmitter is used to send instructions and data, the receiver is used to receive instructions and data, the memory is used to store computer-executable instructions, and the processor is used to execute the computer-executable instructions stored in the memory to implement the method in the above-mentioned embodiment 1.

[0055] It should be noted that the above memory can be independent or integrated with the processor. When the memory is independently provided, the terminal device further includes a bus for connecting the memory and the processor.

[0056] Example 4 This embodiment discloses a computer-readable storage medium, in which computer-executable instructions are stored. When a processor executes the computer-executable instructions, the method in the above-mentioned embodiment 1 is implemented.

[0057] Those of ordinary skill in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be implemented by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, storage, database or other media used in the embodiments provided in this application may include non-volatile and / or volatile memory. Non-volatile memory may include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM) or flash memory. Volatile memory may include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in many forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.

[0058] The technical features of the above embodiments may be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0059] The above-mentioned embodiments only express several implementation modes of the present invention, and the description thereof is relatively specific and detailed, but it cannot be understood as limiting the scope of the invention. It should be pointed out that, for a person of ordinary skill in the art, several modifications and improvements can be made without departing from the concept of the present invention, and these all belong to the protection scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the attached claims.

Claims

1. A method for determining the safety factor range of a cave slope, characterized in that: The method comprises: Get the maximum and minimum values ​​of the slope in the X and Y directions; Setting a vertical division interval, based on the maximum and minimum values ​​of the slope in the X direction, dividing the slope along the X direction according to the vertical division interval to obtain a plurality of vertical units; Based on the maximum and minimum values ​​of the slope in the Y direction, a unit body with a random height is generated in each vertical unit; and the distribution of rock and soil bodies and karst caves in each vertical unit is determined according to the unit body; A slope cave simulation model is constructed based on the distribution of rock and soil bodies and caves in each vertical unit; Performing value calculation on the slope cave simulation model to obtain the slope safety factor; Several groups of slope safety factors are obtained through iterative calculation, and the range of slope safety factors is determined based on each group of slope safety factors.

2. The method for determining the safety factor range of a cave slope according to claim 1, characterized in that: When the slope is divided along the X direction according to the vertical division interval, the number of the vertical units is rounded up.

3. The method for determining the safety factor range of a cave slope according to claim 1, characterized in that: Based on the maximum and minimum values ​​of the slope in the Y direction, a unit cell with random height is generated in each vertical unit; The distribution of rock and soil bodies and karst caves in each vertical unit is determined according to the unit body, including: Based on the exploration data, the initial distribution information of the rock and soil bodies and karst caves of the vertical units is obtained; Uniformly partitioning the vertical unit, randomly selecting any area, and selecting the type of geological body according to the initial distribution information of the rock and soil body and the karst cave of the vertical unit; If the geological body type is rock and soil body, a first unit body representing the height of the rock and soil body is randomly generated according to the initial distribution information of the rock and soil body, and the relevant data is saved; If the geological body type is a cave, a second unit body representing the height of the cave is randomly generated according to the initial distribution information of the cave, and the relevant data is saved; The distribution of rock and soil bodies and karst caves in each vertical unit is determined according to the first unit body and the second unit body.

4. The method for determining the safety factor range of a cave slope according to claim 3, characterized in that: Obtain the initial distribution information of the rock and soil mass and karst caves of the vertical unit, including: According to the exploration data, the initial distribution forms of the rock and soil mass and the karst caves of the slope are obtained, wherein the rock and soil mass is evenly distributed and the karst caves are normally distributed; According to the uniform distribution form of the rock and soil mass of the slope and the normal distribution form of the karst cave, the initial distribution information of the rock and soil mass and the karst cave in each vertical unit is determined.

5. The method for determining the safety factor range of a karst cave slope according to claim 3, characterized in that: Before generating units of random height, it also includes: Preset the maximum and minimum height values; generate units with random heights based on the maximum and minimum height values.

6. The method for determining the safety factor range of a cave slope according to claim 3, characterized in that: When data is saved for the first unit, the saving format is: (X temp , X temp +1, Y _initial , Y _initial +H rock , Rock); When data is saved for the second unit, the saving format is: (X temp , X temp +1, Y _initial , Y _initial +H karst , Karst); Among them, X temp Indicates the current vertical unit coordinate x value; X temp +1 means the vertical unit coordinate x value is increased by 1; Y _initial Indicates the initial value along the Y direction; Y _initial +H rock Represents the y coordinate value of the first unit cell; Y _initial +H karst Represents the y coordinate value of the second unit body; Rock and Karst represent the geological body type.

7. The method for determining the safety factor range of a karst cave slope according to claim 1, characterized in that: The range of slope safety factors is determined based on each group of slope safety factors, including: Several groups of slope safety factors are obtained through iterative calculation, the average value and standard deviation of each group of slope safety factors are solved, and the range of slope safety factors is calculated based on the average value and standard deviation.

8. A device for determining the safety factor range of a cave slope, characterized in that: The device comprises: The initial data acquisition module is used to obtain the maximum and minimum values ​​of the slope in the X and Y directions; A vertical unit division module is used to set a vertical division interval, and based on the maximum and minimum values ​​of the slope in the X direction, divide the slope along the X direction according to the vertical division interval to obtain a number of vertical units; The geological distribution random generation module is used to generate a unit body of random height in each vertical unit based on the maximum and minimum values ​​of the slope in the Y direction; and determine the distribution of rock and soil bodies and karst caves in each vertical unit according to the unit body; A simulation modeling module is used to construct a slope cave simulation model based on the distribution of rock and soil bodies and caves in each vertical unit; A slope safety factor calculation module is used to perform value calculation on the slope cave simulation model to obtain the slope safety factor; The slope safety factor range determination module is used to obtain several groups of slope safety factors through iterative calculation, and determine the range of the slope safety factor based on each group of slope safety factors.

9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, wherein: When the processor executes the computer program, the steps of the method for determining the range of the safety factor of the cave slope described in any one of claims 1 to 7 are implemented.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method for determining the range of the safety factor of a cave slope as described in any one of claims 1 to 7 are implemented.

Images