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

Abstract

The present application relates to the safety factor range determination method, device, equipment and medium of karst cave slope, method includes: obtaining the maximum and minimum of the slope in X direction, Y direction;Based on the maximum and minimum of the slope in X direction, according to the vertical division interval, the slope is divided along X direction, and a plurality of vertical units are obtained;Based on the maximum and minimum of the slope in Y direction, the unit body of random height is generated in each vertical unit;According to the unit body, the distribution of rock-soil mass and karst cave in each vertical unit is determined, and the simulation model of karst cave slope is constructed;The simulation model of karst cave slope is calculated, and the safety factor of the slope is obtained;A plurality of groups of safety factors of the slope are obtained by iterative calculation, and the range of the safety factor of the slope is determined based on each group of safety factors of the slope.The present application can provide reliable data support for the engineering design of karst cave slope, and safety and engineering economy are taken into account.

Description

Technical Field

[0001] This invention relates to the field of slope stability analysis technology, and in particular to a method, apparatus, equipment and medium for determining the safety factor range of a karst cave slope. Background Technology

[0002] A karst slope refers to a slope structure formed in a karst region. Due to the unique characteristics of karst landforms, its slope stability is usually low. Karst caves are formed under the long-term action of groundwater, and have irregular cavities and complex internal structures, making the mechanical properties of the slope difficult to predict.

[0003] Traditional methods for determining the safety factor of karst cave slopes are typically based on limit equilibrium theory and finite element analysis. However, these methods struggle to fully account for the randomness and irregularity of karst caves, especially when the caves are large or located close to the sliding surface. In such cases, the impact of karst caves on the slope is particularly significant, and neglecting this factor can lead to an overestimation of the slope's safety factor, thereby posing a potential safety hazard to the project.

[0004] Therefore, to improve the safety of karst cave slopes, engineering designs often employ over-reinforcement measures, such as significantly increasing the number of anti-slide piles, increasing anchor density, or conducting extensive grouting reinforcement to enhance the slope's safety factor. However, while this over-reinforcement strategy can improve slope stability to some extent, it is often accompanied by high engineering costs and resource waste. Furthermore, because the reinforcement measures are not tailored to the specific morphology and distribution characteristics of the karst caves, they may not effectively address the potential hazards posed by the caves to the slope structure, resulting in limited improvement in the safety factor. In addition, over-reinforcement may alter the original mechanical equilibrium of the slope, increasing local stress concentration or causing unnecessary environmental damage. Summary of the Invention

[0005] Therefore, it is necessary to provide a method, apparatus, equipment, and medium for determining the safety factor range of karst cave slopes, which can calculate the safety factor range of karst cave slopes, provide reliable data support for engineering design, and balance safety and engineering economy, in order to address the above-mentioned technical problems.

[0006] A method for determining the safety factor range of a karst cave slope, the method comprising:

[0007] Obtain the maximum and minimum values ​​of the slope in the X and Y directions;

[0008] Set vertical division intervals. 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 intervals to obtain several vertical units.

[0009] Based on the maximum and minimum values ​​of the slope in the Y direction, units of random height are generated within each vertical unit; the distribution of soil and rock masses and karst caves in each vertical unit is determined according to the units.

[0010] A slope karst cave simulation model is constructed based on the distribution of rock and soil and karst caves in each vertical unit;

[0011] The slope safety factor is obtained by assigning values ​​to the slope karst cave simulation model.

[0012] Several sets of slope safety factors are obtained through iterative calculations, and the range of slope safety factors is determined based on each set of slope safety factors.

[0013] A device for determining the safety factor range of a karst cave slope, the device comprising:

[0014] The initial data acquisition module is used to obtain the maximum and minimum values ​​of the slope in the X and Y directions;

[0015] The vertical unit division module is used to set the vertical division interval. Based on the maximum and minimum values ​​of the slope in the X direction, the slope is divided along the X direction according to the vertical division interval to obtain several vertical units.

[0016] The geological distribution random generation module is used to generate unit cells of random height within each vertical unit based on the maximum and minimum values ​​of the slope in the Y direction; and to determine the distribution of soil and rock masses and karst caves in each vertical unit based on the unit cells.

[0017] The simulation modeling module is used to construct a slope karst cave simulation model based on the distribution of rock and soil and karst caves in each vertical unit;

[0018] The slope safety factor calculation module is used to assign values ​​to the slope karst cave simulation model and obtain the slope safety factor.

[0019] The slope safety factor range determination module is used to obtain several sets of slope safety factors through iterative calculations, and to determine the range of the slope safety factor based on each set of slope safety factors.

[0020] A computer device includes a memory and a processor, the memory storing a computer program, and the processor executing the computer program to implement the steps of the method for determining the safety factor range of a karst cave slope.

[0021] A computer-readable storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor, implements the steps of the method for determining the safety factor range of a karst cave slope.

[0022] The aforementioned method, apparatus, equipment, and medium for determining the safety factor range of karst slopes involve: obtaining the maximum and minimum values ​​of the slope in the X and Y directions; setting vertical division intervals; dividing the slope along the X direction based on the maximum and minimum values ​​in the X direction to obtain several vertical units; generating units of random height within each vertical unit based on the maximum and minimum values ​​in the Y direction; determining the distribution of soil and rock masses and karst caves in each vertical unit based on the unit units; constructing a slope-karst cave simulation model based on the distribution of soil and rock masses and karst caves in each vertical unit; assigning values ​​to the slope-karst cave simulation model to obtain the slope safety factor; obtaining several sets of slope safety factors through iterative calculations; and determining the range of the slope safety factor based on each set of slope safety factors.

[0023] This invention discretizes the slope into several vertical units through reasonable vertical partitioning intervals, making it closer to the geological distribution during vertical drilling in exploration. By generating random heights for the soil and rock masses and karst caves within the vertical units, the complex geological structure of intersecting soil and rock masses and karst caves can be effectively simulated, resulting in simulation results that are closer to reality and thus yielding a more accurate slope safety factor. Furthermore, due to the randomness of the generated karst cave heights, multiple iterative calculations avoid the bias of a single result and reduce deviations caused by uncertainties in the calculation, thereby solving for a more reliable range of slope safety factors. This provides reliable data support for subsequent assessment of slope stability and optimization of engineering design, and also avoids over-design, saving engineering costs. Attached Figure Description

[0024] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on the structures shown in these drawings without creative effort.

[0025] Figure 1 This is a flowchart illustrating a method for determining the safety factor range of a karst cave slope in one embodiment.

[0026] Figure 2 This is a schematic diagram of the outer contour of a slope in one embodiment;

[0027] Figure 3 This is a schematic diagram illustrating the division of a slope into several vertical units along the X direction in one embodiment;

[0028] Figure 4 This is a schematic diagram showing the distribution of soil and rock masses and karst caves in a vertical unit generated along the Y direction in one embodiment;

[0029] Figure 5 This is an initial cave model generated in one embodiment;

[0030] Figure 6 This is a schematic diagram of the outer contour model of a slope constructed using a FLAC3D numerical model in one embodiment.

[0031] Figure 7 This is a schematic diagram of a slope karst cave simulation model constructed using the FLAC3D numerical model in one embodiment.

[0032] Figure 8 A structural block diagram of a device for determining the safety factor range of a karst cave slope in one embodiment;

[0033] Figure 9 This is an internal structural diagram of a computer device in one embodiment.

[0034] The realization of the objective, functional features and advantages of the present invention will be further explained in conjunction with the embodiments and with reference to the accompanying drawings. Detailed Implementation

[0035] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative effort are within the scope of protection of the present invention.

[0036] 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 positional relationship and movement of each component in a certain specific posture (as shown in the figure). If the specific posture changes, the directional indication will also change accordingly.

[0037] Furthermore, in this invention, descriptions involving "first," "second," etc., are for descriptive purposes only and should not be construed as indicating or implying their relative importance or implicitly specifying the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include at least one of that feature. In the description of this invention, "a plurality of" means at least two, such as two, three, etc., unless otherwise explicitly specified.

[0038] It is understood that the technical solutions of the various embodiments of the present invention can be combined with each other, but only if they are based on the ability of those skilled in the art to implement them. When the combination of technical solutions is contradictory or cannot be implemented, it should be considered that such combination of technical solutions does not exist and is not within the scope of protection claimed by the present invention.

[0039] The embodiments of the present invention will now be described in detail with reference to the accompanying drawings.

[0040] Example 1

[0041] This embodiment discloses a method for determining the safety factor range of a karst slope. By rationally dividing the slope into vertical units, it more closely approximates the geological distribution during vertical drilling in geological exploration. However, since geological exploration involves vertical point drilling, information such as the height and distribution of karst caves in the exploration data is obtained statistically, resulting in some uncertainty in the cave height. This invention generates random heights for the soil and rock masses and karst caves within the vertical units, effectively simulating the complex geological structure of the intersecting soil and rock masses and karst caves. This makes the simulation results closer to reality, leading to a more accurate slope safety factor. Furthermore, because the generated cave heights have a certain degree of randomness, multiple iterative calculations avoid the one-sidedness of a single result and reduce deviations caused by uncertainties in the calculation. This allows for the determination of a more reliable slope safety factor range, providing reliable data support for subsequent slope stability assessment and engineering design optimization of karst slopes. It also avoids over-design and saves engineering costs.

[0042] like Figure 1 As shown, a method for determining the safety factor range of a karst cave slope is provided, including the following steps:

[0043] Step 201: Obtain the maximum and minimum values ​​of the slope in the X and Y directions.

[0044] Step 202: Set the vertical division interval. 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 several vertical units.

[0045] Step 203: Based on the maximum and minimum values ​​of the slope in the Y direction, generate unit cells of random height within each vertical unit; determine the distribution of soil and rock masses and karst caves in each vertical unit based on the unit cells.

[0046] Step 204: Construct a slope karst cave simulation model based on the distribution of rock and soil and karst caves in each vertical unit.

[0047] Step 205: Assign values ​​to the slope karst cave simulation model to obtain the slope safety factor.

[0048] Step 206: Several sets of slope safety factors are obtained through iterative calculation, and the range of slope safety factors is determined based on each set of slope safety factors.

[0049] In the specific implementation of step 201, the outer contour information of the slope is first obtained through geological exploration data, and the coordinate values ​​of each corner point of the slope to be calculated are determined. (See reference...) Figure 2 This embodiment shows the outer contour shape of the slope and the coordinate values ​​of each corner point of the slope in one of the scenarios provided. Based on 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 X max And determine the minimum value Y of the slope in the Y direction. min With the maximum value Y max .

[0050] In the specific implementation of step 202, based on information such as the location of drilling points in the exploration survey, vertical division intervals are set, and then based on the minimum value X of the slope in the X direction... min With the maximum value X max The slope is divided into multiple long, narrow vertical units along the X-direction, such as... Figure 3 The image shows the vertical elements created when the vertical interval is set to 1m. The number of vertical elements, n, is denoted by Ceil(X). max -X min Ceil indicates rounding up.

[0051] It is understandable that the distribution of karst caves in exploration data is obtained through statistical analysis of local geological information from each borehole, thus introducing uncertainty. Therefore, a reasonable vertical partitioning interval is set based on the distribution of drilling points. This interval discretizes the slope into several vertical units. In subsequent calculations, it is equivalent to considering only the height distribution of karst caves within each vertical unit, simplifying the calculation process. Furthermore, this partitioning confines uncertainty to each vertical unit, preventing it from accumulating throughout the entire model, thus improving the accuracy of karst cave estimation and obtaining a more accurate slope safety factor.

[0052] In the specific implementation process of step 203, firstly, based on the exploration data, the distribution of the rock and soil mass and karst caves of the entire slope is obtained; based on the distribution of the rock and soil mass and karst caves of the entire slope, the initial distribution pattern of the rock and soil mass and 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 pattern of the rock and soil mass and the normal distribution pattern of the karst caves of the slope, the initial distribution information of the rock and soil mass and karst caves in each vertical unit is determined.

[0053] Then, the vertical units are uniformly divided into zones. Any region is randomly selected, and the geological body type is chosen based on the initial distribution information of the soil and rock mass and karst caves within the vertical units. If the geological body type is soil and rock mass, a first unit representing the height of the soil and rock mass is randomly generated based on the initial distribution information, and the relevant data is saved. When saving the data for the first unit, the saving format is: (X... temp X temp +1, Y _initial Y _initial +H rock Rock).

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

[0055] Among them, X temp Indicates the current vertical cell coordinate x-value; X temp +1 indicates that the x-value of the vertical unit coordinate is increased by 1; Y _initial This represents the initial value along the Y direction; Y _initial +H rock This represents the y-coordinate value of the first unit cell; Y _initial +H karst This represents the y-coordinate value of the second unit body; Rock and Karst represent the geological body types.

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

[0057] In addition, when generating units with random heights, a maximum height and a minimum height are preset; units with random heights are generated based on the maximum height and the minimum height.

[0058] Specifically, assume that the preset minimum height is 1m and the maximum height is 4m.

[0059] First, initialize the initial value Y in the Y direction. _initial =Y min .

[0060] Then, the vertical units are uniformly divided, any area is randomly selected, and the geological body type (rock and soil body or karst cave) is selected according to the initial distribution information of the rock and soil body and karst cave in the vertical unit, where the probability of rock and soil body and karst cave being selected is 0.5.

[0061] If the selected area's geological body type is rock and soil, a random height of H will be generated between the minimum height of 1m and the maximum height of 4m. rock The first unit is recorded, and the soil and rock mass model data is saved in the format (X). temp X temp +1, Y _initial Y _initial +H rock (Rock). If the selected area's geological type is a karst cave, a random height of H will be generated between the minimum height of 1m and the maximum height of 4m. karst The second unit body, and at the same time save the data of the cave, the save format is (X temp X temp +1, Y _initial Y _initial +H karst (Karst).

[0062] When Y _initial +H rock Or Y _initial +H karst Greater than or equal to Y max When this time, it indicates that the analysis of the distribution of soil and rock masses and karst caves in the current vertical unit has been completed, thus generating a result such as Figure 4 The diagram shows an initial, elongated cave model. When a randomly selected area is shown, its Y-direction value is Y. max When saving data modified 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 .

[0063] Repeat the above steps until the distribution of soil and rock masses and karst caves in all vertical units has been analyzed, and the results are as follows: Figure 5 The initial cave model shown.

[0064] In the specific implementation of step 204, firstly according to Figure 2 The outer contour of the slope is generated using the FLAC3D numerical model to create a complete outer contour model of the slope, such as... Figure 6 As shown, the model consists entirely of rock and soil, with no karst caves.

[0065] Then, based on the second unit data stored in step 203, in Figure 6 In the slope outer contour model, obtain the position coordinates (x, y) of the corresponding element; when the corresponding element in the slope outer contour model is the second element, delete that element, and finally generate the model as shown below. Figure 7 The simulated model of the karst cave on the slope is shown.

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

[0067] Table 1. Mechanical parameters of soil and rock

[0068]

[0069] In the specific implementation of step 206, since the height of the karst caves in the vertical unit has a certain degree of randomness, the slope safety factor calculated through a slope karst cave simulation model has a certain degree of bias. Therefore, several sets of slope safety factors are obtained through iterative calculation, and then the average value and standard deviation of each set 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 as follows:

[0070] ;

[0071] In the formula, This indicates the minimum range of the slope safety factor; This indicates the maximum range of the slope safety factor. This represents the average value; It represents the standard deviation.

[0072] Although this embodiment Figure 1 The steps are shown sequentially as indicated by the arrows, but they are not necessarily executed in the order indicated by the arrows. Unless otherwise specified in this document, there is no strict order in which these steps are performed; they can be executed in other orders. Figure 1 At least some of the steps in the process may include multiple sub-steps or multiple stages. These sub-steps or stages are not necessarily completed 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 some of the sub-steps or stages of other steps.

[0073] Example 2

[0074] Based on the method for determining the safety factor range of a karst cave slope in Embodiment 1, this embodiment discloses a device for determining the safety factor range of a karst cave slope, as shown in Figure 8. The device includes: 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:

[0075] The initial data acquisition module 401 is used to acquire the maximum and minimum values ​​of the slope in the X and Y directions.

[0076] The vertical unit division module 402 is used to set the vertical division interval. Based on the maximum and minimum values ​​of the slope in the X direction, the slope is divided along the X direction according to the vertical division interval to obtain several vertical units.

[0077] The random geological distribution generation module 403 is used to generate units of random height within each vertical unit based on the maximum and minimum values ​​of the slope in the Y direction; and to determine the distribution of soil and rock masses and karst caves in each vertical unit based on the units.

[0078] The simulation modeling module 404 is used to construct a slope karst cave simulation model based on the distribution of rock and soil and karst caves in each vertical unit.

[0079] The slope safety factor calculation module 405 is used to assign values ​​to the slope karst cave simulation model to obtain the slope safety factor.

[0080] The slope safety factor range determination module 406 is used to obtain several sets of slope safety factors through iterative calculation, and to determine the range of slope safety factors based on each set of slope safety factors.

[0081] 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 Embodiment 1, and therefore will not be described again in this embodiment. Each unit module can be implemented entirely or partially through software, hardware, or a combination thereof. Each unit module can be embedded in or independent of the processor in a computer device in hardware form, or it can be stored in the memory of a computer device in software form, so that the processor can call and execute the operations corresponding to each of the above unit modules.

[0082] Example 3

[0083] like Figure 9The diagram illustrates a terminal device disclosed in this embodiment, comprising a transmitter, a receiver, a memory, and a processor. The transmitter transmits instructions and data, the receiver receives instructions and data, the memory stores computer-executed instructions, and the processor executes the computer-executed instructions stored in the memory to implement the method described in Embodiment 1 above.

[0084] It is important to note that the aforementioned memory can be either standalone or integrated with the processor. When the memory is set up independently, the terminal device also includes a bus for connecting the memory and the processor.

[0085] Example 4

[0086] This embodiment discloses a computer-readable storage medium storing computer-executable instructions. When a processor executes the computer-executable instructions, it implements the method in Embodiment 1 above.

[0087] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium. When executed, the computer program can include the processes of the embodiments of the above methods. Any references to memory, storage, databases, or other media used in the embodiments provided in this application can 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 a variety of forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), dual data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link DRAM (SLDRAM), RAMbus direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.

[0088] The technical features of the above embodiments can be combined in any way. For the sake of brevity, 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.

[0089] The embodiments described above are merely illustrative of several implementations of the present invention, and while the descriptions are specific and detailed, they should not be construed as limiting the scope of the invention. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of the present invention, and these modifications and improvements all fall within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be determined by the appended claims.

Claims

1. A method for determining the safety factor range of a karst cave slope, characterized in that, The method includes: Obtain the maximum and minimum values ​​of the slope in the X and Y directions; Set vertical division intervals. 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 intervals to obtain several vertical units. Based on the maximum and minimum values ​​of the slope in the Y direction, units of random height are generated within each vertical unit; the distribution of soil and rock masses and karst caves in each vertical unit is determined according to the units. A slope karst cave simulation model is constructed based on the distribution of rock and soil and karst caves in each vertical unit; The slope safety factor is obtained by assigning values ​​to the slope karst cave simulation model. Several sets of slope safety factors are obtained through iterative calculations, and the range of slope safety factors is determined based on each set of slope safety factors. Based on the maximum and minimum values ​​of the slope in the Y direction, elements of random height are generated within each vertical unit. The distribution of soil and rock masses and karst caves within each vertical unit is determined according to these elements, including: Based on the exploration data, the initial distribution information of the rock and soil mass and karst caves in the vertical unit was obtained; The vertical unit is uniformly divided into zones, any area is randomly selected, and the geological body type is selected based on the initial distribution information of the rock and soil and karst caves in the vertical unit. If the geological body type is rock and soil, then the first unit representing the height of the rock and soil body is randomly generated based on the initial distribution information of the rock and soil body, and the relevant data is saved; If the geological body type is a karst cave, a second unit representing the height of the karst cave is randomly generated based on the initial distribution information of the karst cave, and the relevant data is saved; Based on the first unit and the second unit, the distribution of rock and soil and karst caves in each vertical unit is determined.

2. The method for determining the safety factor range of a karst cave slope according to claim 1, characterized in that, When dividing the slope 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 karst cave slope according to claim 1 or 2, characterized in that, Obtain the initial distribution information of the soil and rock mass and karst caves in the vertical unit, including: Based on the exploration data, the initial distribution pattern of the rock and soil mass and karst caves on the slope was obtained, wherein the rock and soil mass is uniformly distributed and the karst caves are normally distributed; Based on the uniform distribution pattern of the soil and rock mass of the slope and the normal distribution pattern of the karst caves, the initial distribution information of the soil and rock mass and the karst caves in each vertical unit is determined.

4. The method for determining the safety factor range of a karst cave slope according to claim 1 or 2, characterized in that, Before generating units of random height, the following is also included: Preset maximum and minimum height values; generate unit cells with random heights based on the maximum and minimum height values.

5. The method for determining the safety factor range of a karst cave slope according to claim 1 or 2, characterized in that, When saving data to the first unit body, the saving format is: (X temp X temp +1, Y _initial Y _initial +H rock Rock); When saving data to 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 cell coordinate x-value; X temp +1 indicates that the x-value of the vertical unit coordinate is increased by 1; Y _initial This represents the initial value along the Y direction; Y _initial +H rock This represents the y-coordinate value of the first unit cell; Y _initial +H karst This represents the y-coordinate value of the second unit body; Rock and Karst represent the geological body types.

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

7. A device for determining the safety factor range of a karst cave slope, characterized in that, The device includes: The initial data acquisition module is used to obtain the maximum and minimum values ​​of the slope in the X and Y directions; The vertical unit division module is used to set the vertical division interval. Based on the maximum and minimum values ​​of the slope in the X direction, the slope is divided along the X direction according to the vertical division interval to obtain several vertical units. The geological distribution random generation module is used to generate unit cells of random height within each vertical unit based on the maximum and minimum values ​​of the slope in the Y direction; and to determine the distribution of soil and rock masses and karst caves in each vertical unit based on the unit cells. The simulation modeling module is used to construct a slope karst cave simulation model based on the distribution of rock and soil and karst caves in each vertical unit; The slope safety factor calculation module is used to assign values ​​to the slope karst cave simulation model and obtain the slope safety factor. The slope safety factor range determination module is used to obtain several sets of slope safety factors through iterative calculations, and to determine the range of the slope safety factor based on each set of slope safety factors; In the geological distribution random generation module, based on the maximum and minimum values ​​of the slope in the Y direction, units of random height are generated within each vertical unit; the distribution of soil and rock masses and karst caves in each vertical unit is determined according to the units, including: Based on the exploration data, the initial distribution information of the rock and soil mass and karst caves in the vertical unit was obtained; The vertical unit is uniformly divided into zones, any area is randomly selected, and the geological body type is selected based on the initial distribution information of the rock and soil and karst caves in the vertical unit. If the geological body type is rock and soil, then the first unit representing the height of the rock and soil body is randomly generated based on the initial distribution information of the rock and soil body, and the relevant data is saved; If the geological body type is a karst cave, a second unit representing the height of the karst cave is randomly generated based on the initial distribution information of the karst cave, and the relevant data is saved; Based on the first unit and the second unit, the distribution of rock and soil and karst caves in each vertical unit is determined.

8. A computer device comprising a memory and a processor, wherein the memory stores a computer program, characterized in that, When the processor executes the computer program, it implements the steps of the method for determining the safety factor range of a karst cave slope according to any one of claims 1 to 6.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it implements the steps of the method for determining the safety factor range of the karst cave slope according to any one of claims 1 to 6.

Images