Pull-off position prediction method, device, equipment and readable storage medium
By selecting target rock layers in a layered rock mass wading slope, calculating the total tensile stress and constructing a stress balance equation, and considering the coupling effect of pore water pressure and fracture propagation, the problem of inaccurate prediction of tensile fracture location in existing technologies is solved, achieving more accurate disaster early warning and enhanced safety.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2026-06-09
- Publication Date
- 2026-07-10
AI Technical Summary
Existing technologies for predicting the breakage location of flexural-tensile type reservoir landslides on steep slopes suffer from discrepancies between model predictions and actual locations. This leads to false alarms and missed warnings in disaster early warning, hindering timely and accurate disaster warnings and posing safety hazards.
By selecting target rock layers from the layered rock mass wading slope, calculating the total tensile stress, and combining the self-weight of the rock layer, the force of the upper rock layer, and the pore water pressure, a stress balance equation is constructed. The minimum value between the equivalent tensile strength of crack propagation and the tensile strength of water-rock softening is selected as the critical stress for tensile fracture. The fracture location is solved using the length parameter, taking into account the coupling effect of pore water pressure on the anti-sliding capacity and crack propagation.
It significantly improves the accuracy and timeliness of tensile failure location prediction, enabling earlier identification of potential disasters, reducing false alarms and missed alarms, ensuring safety, and guiding the precise deployment of monitoring equipment or reinforcement measures in engineering practice to avoid cost waste and reinforcement failure.
Smart Images
Figure CN122364608A_ABST
Abstract
Description
Technical Field
[0001] This application relates to the field of geological disaster prevention and control technology, and in particular to a method, device, equipment and readable storage medium for predicting the location of a break. Background Technology
[0002] Current research on reservoir landslides with bending and tensile fractures on steep slopes mainly focuses on theoretical discussions and qualitative analyses of the failure mechanisms, while research on predicting the specific fracture locations of landslide tensile fractures is lacking. Existing prediction methods are mostly based on simplified geological assumptions, leading to discrepancies between model predictions and actual fracture locations.
[0003] Furthermore, relying on model predictions for early warning of landslides in steep slopes with bending and tensile cracking is prone to false alarms and missed alarms, and cannot accurately achieve early warning of landslides in steep slopes with bending and tensile cracking. Once a disaster occurs, it will seriously threaten the lives and property of people around the disaster site and the safety of construction workers and machinery at the rescue site.
[0004] In summary, the current prediction of the breakage location of the bending and tensile-fracturing reservoir on the steep slope is not accurate enough, and timely and accurate disaster warnings cannot be provided, posing a safety hazard. Summary of the Invention
[0005] This application provides a method, apparatus, equipment, and readable storage medium for predicting the breakage location, aiming to solve the technical problem that the current prediction of the breakage location of bent-twist reservoirs on steep slopes is not accurate enough, and the disaster warning cannot be given in a timely and accurate manner, which poses a safety hazard.
[0006] In a first aspect, this application provides a method for predicting the location of tensile fractures, applicable to layered rock mass wading slopes, the method comprising: Select target rock layers from the layered rock mass wading slope; Calculate the total tensile stress of the target rock layer based on its own weight, the force exerted by the upper rock layer and the pore water pressure. The minimum value between the equivalent tensile strength of the target rock layer fracture propagation and the tensile strength of water-rock softening is selected as the critical stress for tensile fracture of the target rock layer. To make the total tensile stress of the target rock layer equal to the critical stress at break, a stress balance equation is constructed. By using the length parameter of the target rock layer, the stress balance equation is solved to obtain the fracture location of the target rock layer.
[0007] Optionally, the calculation of the total tensile stress of the target rock layer based on its own weight, the force exerted by the overlying rock layers, and the pore water pressure includes: The total tensile stress of the target rock layer is obtained by superimposing the tensile stress generated by the bending moment of the target rock layer, the compressive stress generated by the gravity component of the upper rock layer, the compressive stress caused by interlayer friction, and the tensile stress generated by the axial shear force. Among them, the compressive stress caused by interlayer friction is calculated based on the interlayer friction minus the pore water pressure.
[0008] Optionally, before obtaining the total tensile stress of the target rock layer by superimposing the tensile stress generated by the bending moment of the target rock layer, the compressive stress generated by the gravity component of the upper rock layer, the compressive stress caused by interlayer friction, and the tensile stress generated by the axial shear force, the following steps are included: Based on the bending moment at the fracture point of the target rock stratum and the bending strength of the target rock stratum, a cantilever plate beam model is used to calculate the tensile stress generated by the bending moment of the target rock stratum.
[0009] Optionally, before selecting the minimum value between the equivalent tensile strength of the target rock layer's fracture propagation and the tensile strength of water-rock softening as the critical tensile stress for the target rock layer, the following steps are included: Based on the fracture parameters, stress intensity factor, and geometric correction factor of the target rock stratum, the equivalent tensile strength of the fracture propagation in the target rock stratum is calculated using the equivalent tensile strength formula.
[0010] Optionally, before selecting the minimum value between the equivalent tensile strength of the target rock layer fracture propagation and the tensile strength of water-rock softening as the critical tensile stress for the target rock layer, the method further includes: Based on the initial tensile strength, water content, and water-rock softening coefficient of the target rock stratum, the tensile strength of the target rock stratum under water-rock softening is calculated.
[0011] Optionally, the stress balance equation is: ; in, Let x be the total tensile stress at location x on the target rock stratum. To the critical stress at which tensile strength is broken, The equivalent tensile strength for the propagation of fractures in the target rock strata. Let Y be the fracture toughness of the target rock layer, Y be the geometric correction factor, and v be the fracture half-length of the target rock layer. The tensile strength of the target rock strata softened by water. k represents the initial tensile strength of the target rock layer. w ω is the water-rock softening coefficient, and ω is the water content.
[0012] Optionally, the step of using the length parameter of the target rock layer to solve the stress balance equation to obtain the tensile fracture location of the target rock layer includes: Based on the stress balance equation, the position x on the target rock layer is solved using an iterative algorithm within the length L of the target rock layer to obtain the tensile fracture position x of the target rock layer.
[0013] Secondly, this application provides a tensile fracture location prediction device, applied to a layered rock mass wading slope, the tensile fracture location prediction device comprising: The selection module is used to select target rock layers from a layered rock mass wading slope. The calculation module is used to calculate the total tensile stress of the target rock layer based on its own weight, the force exerted by the upper rock layer, and the pore water pressure. The determination module is used to select the minimum value between the equivalent tensile strength of the target rock layer fracture propagation and the tensile strength of water-rock softening as the critical stress for tensile fracture of the target rock layer. A module is constructed to make the total tensile stress of the target rock layer equal to the critical stress at break, and to construct the stress balance equation. The solver module is used to solve the stress balance equation using the length parameter of the target rock layer to obtain the fracture location of the target rock layer.
[0014] Thirdly, this application provides a breakage position prediction device, which includes a processor, a memory, and a breakage position prediction program stored in the memory and executable by the processor, wherein when the breakage position prediction program is executed by the processor, it implements the steps of the breakage position prediction method as described above.
[0015] Fourthly, this application provides a readable storage medium storing a breakage position prediction program, wherein when the breakage position prediction program is executed by a processor, it implements the steps of the breakage position prediction method as described above.
[0016] The beneficial effects of the technical solution provided in this application include: In this application, a target rock layer is selected from a layered rock mass wading slope; the total tensile stress of the target rock layer is calculated based on its own weight, the force exerted by the upper rock layer, and the pore water pressure; the minimum value between the equivalent tensile strength of the target rock layer due to fracture propagation and the tensile strength due to water softening is selected as the critical stress for tensile fracture of the target rock layer; the total tensile stress of the target rock layer is made equal to the critical stress for tensile fracture, and a stress balance equation is constructed; the stress balance equation is solved using the length parameter of the target rock layer to obtain the location of the fracture of the target rock layer. To simplify modeling, this application proposes a representative single target rock layer for wading slopes in layered rock masses. Firstly, when calculating the total tensile stress of the target rock layer, in addition to considering its own weight and the force exerted by the upper rock layers, pore water pressure is also taken into account. Since the tensile strength of the rock layer is much lower than its compressive strength, the tensile strength of the target rock layer is used to determine its fracture. The propagation of fissures in the target rock layer weakens its tensile strength. Simultaneously, the water-softening effect also weakens its tensile strength. Based on the barrel effect, rock mass failure occurs along the path of least resistance. Therefore, the equivalent tensile strength of the target rock layer's fissure propagation and the water-softening effect are selected. The minimum value of the tensile strength is taken as the critical stress for tensile failure. When determining the critical stress for tensile failure, the interaction between the crack propagation and water softening effect of the target rock layer is comprehensively considered. When the total tensile stress of the target rock layer is greater than the critical stress for tensile failure, it means that the target rock layer will fail. The stress balance equation is constructed by making the total tensile stress equal to the critical stress for tensile failure. In the stress balance equation, pore water pressure is considered both as a load (which will increase the total tensile stress by reducing interlayer friction) and as a weakening factor of tensile strength. The synergy of the two will significantly increase the net tensile stress and accelerate crack propagation, causing the rock layer to fail at a stress threshold far below the dry state. By innovatively incorporating the coupling effect of the propagation and evolution of internal rock fractures with pore water pressure into the mechanical model, this study overcomes the key technical challenge of traditional analysis methods failing to accurately capture the potential instability mechanism of layered rock mass water-crossing slopes. It reveals from a mechanical perspective the mechanism by which the water pressure wedging effect and stress concentration at the fracture tip synergistically drive the tensile fracturing of rock strata. This allows for more accurate determination of the tensile fracture location of layered rock mass water-crossing slopes, enabling timely disaster warnings and the elimination of safety hazards. Attached Figure Description
[0017] Figure 1 This is a flowchart illustrating an embodiment of the pull-off location prediction method of this application; Figure 2 This is a schematic diagram of a cross-section of a layered rock mass wading slope according to an embodiment of the tensile fracture location prediction method of this application; Figure 3 This is a schematic diagram of the stress on a layered rock mass wading slope according to an embodiment of the tensile fracture location prediction method of this application. Figure 4This is a schematic diagram of the target rock layer stress analysis according to an embodiment of the tensile fracture location prediction method of this application; Figure 5 This is a schematic diagram of the functional modules of an embodiment of the breakage position prediction device of this application; Figure 6 This is a schematic diagram of the hardware structure of the breakage location prediction device involved in the embodiments of this application. Detailed Implementation
[0018] To enable those skilled in the art to better understand the present application, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present application, and not all embodiments. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without creative effort are within the scope of protection of the present application.
[0019] To make the objectives, technical solutions, and advantages of this application clearer, the embodiments of this application will be described in further detail below with reference to the accompanying drawings.
[0020] In a first aspect, embodiments of this application provide a method for predicting the breakage location.
[0021] In one embodiment, reference is made to Figure 1 , Figure 1 This is a flowchart illustrating an embodiment of the breakage location prediction method of this application, as shown below. Figure 1 As shown, the methods for predicting the location of tensile fractures in layered rock mass wading slopes include: Step S10: Select the target rock layer from the layered rock mass wading slope.
[0022] In this embodiment, refer to Figure 2 , Figure 2 This is a schematic diagram of a cross-section of a layered rock mass wading slope according to an embodiment of the tensile strength location prediction method of this application, as shown below. Figure 2 As shown, representative single rock layers in the slope can be identified as target rock layers using geological exploration data. The dip angle of the layered rock mass wading slope is greater than a threshold, for example, set to 40 degrees. Figure 2 The rock strata at the predicted fault location are selected as target rock strata. Selection principles can be based on the lithological characteristics, degree of fracture development, and deformation indicators of the rock strata; for example, priority can be given to rock strata with the most significant deformation or those located at critical stress points. Simplifying the complex layered slope model into a single target rock stratum model not only significantly reduces the complexity of mechanical modeling but also focuses on the critical strata most likely to fail. This simplification significantly improves prediction efficiency while maintaining computational accuracy, enabling engineers to quickly identify risk sources and lay the foundation for subsequent detailed mechanical analysis.
[0023] Step S20: Calculate the total tensile stress of the target rock layer based on its own weight, the force exerted by the upper rock layer, and the pore water pressure.
[0024] In this embodiment, the calculation of total tensile stress is the core of predicting the fracture location. Traditional methods often ignore pore water pressure or treat it as a single load. However, this step considers pore water pressure as a key variable, taking into account both the thrust it exerts on the rock strata as a load and the weakening of its anti-sliding capacity after infiltrating the rock mass. By comprehensively superimposing the target rock stratum's own gravity load, the lateral pressure and gravity components transmitted from the upper rock strata, and the mechanical effects of pore water pressure, a stress field that better reflects the actual operating conditions of the reservoir can be constructed. This multi-factor coupled calculation method effectively overcomes the problem of large prediction deviations in traditional dry or simplified models under reservoir water level fluctuation scenarios, significantly improving the realism and reliability of stress calculation.
[0025] Step S30: Select the minimum value between the equivalent tensile strength of the target rock layer fracture propagation and the tensile strength of water softening as the critical stress for tensile fracture of the target rock layer.
[0026] In this embodiment, rock mass failure often follows the "barrel effect," meaning it occurs along the path of least resistance. The target rock stratum may simultaneously experience strength reduction due to primary fracture propagation and material softening due to water-rock interaction. By calculating the equivalent tensile strength under fracture propagation (based on fracture mechanics) and the tensile strength under water-rock softening (based on a water content softening model), and taking the minimum of the two as the critical stress for tensile failure, the bearing capacity limit of the target rock stratum can be conservatively and accurately defined. This dual verification mechanism avoids the overestimation risk that may arise from single-index assessment, ensuring timely triggering of early warnings under either fracture development or water softening conditions, greatly enhancing the safety of disaster identification.
[0027] Step S40: Make the total tensile stress of the target rock layer equal to the critical stress at break, and construct the stress balance equation.
[0028] In this embodiment, the critical state of tensile failure is defined as the total tensile stress generated within the target rock stratum reaches its maximum withstandable critical stress. By setting the total tensile stress equal to the critical stress for failure, a stress balance equation is established that includes location variables (such as the distance x from the top of the target rock stratum). This equation directly links the mechanical state with the spatial location, transforming the abstract stability assessment into a concrete mathematical problem. Constructing this balance equation is a crucial step in achieving quantitative prediction. It provides clear mathematical constraints for subsequently solving for the specific failure location, ensuring that the prediction results no longer rely on qualitative judgments based on human experience but are based on rigorous mechanical equilibrium principles.
[0029] Step S50: Using the length parameter of the target rock layer, solve the stress balance equation to obtain the fracture location of the target rock layer.
[0030] In this embodiment, the geometric length parameter of the target rock layer is used as a boundary condition to solve the stress balance equation mentioned above. This allows for the calculation of the specific location coordinate x that satisfies the failure condition, which is the predicted landslide break point. Obtaining specific spatial coordinates through numerical solution enables geological disaster prevention personnel to accurately determine where monitoring equipment or reinforcement measures need to be deployed. This quantitative location output directly guides engineering practice, avoiding cost waste caused by excessively large reinforcement areas due to ambiguous locations, or reinforcement failure due to location deviations, thus achieving an optimal balance between disaster prevention investment and safety benefits.
[0031] In this embodiment, through steps S10 to S50 above, to simplify modeling, for the water-eroded slope of layered rock mass, a representative single target rock layer can be selected for modeling. First, when calculating the total tensile stress of the target rock layer, in addition to considering the self-weight of the target rock layer and the force exerted by the upper rock layer, the pore water pressure is also taken into account. This multi-factor coupled calculation method effectively overcomes the problem of large prediction deviations in traditional dry models or simplified models under reservoir water level fluctuation scenarios, and significantly improves the realism and reliability of stress calculation. The tensile strength of rock strata is much lower than their compressive strength. Using the tensile strength of the target rock strata to determine its fracture strength is problematic. The propagation of fissures in the target rock strata weakens its tensile strength, and the softening effect of water also weakens it. Based on the "barrel effect," rock mass failure occurs along the path of least resistance. Therefore, the minimum of the equivalent tensile strength due to fissure propagation and the tensile strength due to water softening is selected as the critical stress for fracture rupture. In determining the critical stress, the coupling effect between fissure propagation and water softening is comprehensively considered. This dual-verification mechanism avoids the overestimation risk that may arise from a single indicator assessment, ensuring timely triggering of early warnings under either unfavorable conditions of fissure development or water softening, greatly enhancing the safety of disaster identification. When the total tensile stress of the target rock stratum exceeds the critical stress for tensile failure, it signifies tensile failure of the target rock stratum. A stress balance equation is constructed by equalizing the total tensile stress to the critical stress for tensile failure. In this equation, pore water pressure is considered both as a load (increasing the total tensile stress by reducing interlayer friction) and as a weakening factor in tensile strength. The combined effect of these two factors significantly increases the net tensile stress and accelerates crack propagation, causing the rock stratum to fail tensilely at stress levels far below the dry stress threshold. Numerical solutions yield specific spatial coordinates, enabling geological disaster prevention personnel to accurately determine the locations where monitoring equipment or reinforcement measures need to be deployed. This quantitative location output directly guides engineering practice, avoiding cost waste due to excessively large reinforcement areas caused by ambiguous locations, or reinforcement failure due to location deviations, thus achieving an optimal balance between disaster prevention investment and safety benefits. By innovatively incorporating the coupling effect of the propagation and evolution of internal rock fractures with pore water pressure into the mechanical model, this study overcomes the key technical challenge of traditional analysis methods failing to accurately capture the potential instability mechanism of layered rock mass water-crossing slopes. It reveals from a mechanical perspective the mechanism by which the water pressure wedging effect and stress concentration at the fracture tip synergistically drive the tensile fracturing of rock strata. This allows for more accurate determination of the tensile fracture location of layered rock mass water-crossing slopes, enabling timely disaster warnings and the elimination of safety hazards.
[0032] Further, in one embodiment, step S20 includes: The total tensile stress of the target rock layer is obtained by superimposing the tensile stress generated by the bending moment of the target rock layer, the compressive stress generated by the gravity component of the upper rock layer, the compressive stress caused by interlayer friction, and the tensile stress generated by the axial shear force. Among them, the compressive stress caused by interlayer friction is calculated based on the interlayer friction minus the pore water pressure.
[0033] In this embodiment, the total tensile stress is not the result of a single force source, but rather the superposition of multiple stress components. Bending moment primarily originates from the cantilever effect of the target rock stratum and is the main source of tensile stress; the gravity component and interlayer friction typically generate compressive stress, which counteracts the tensile stress; while axial shear force may generate additional tensile stress. Specifically, pore water pressure reduces the effective normal stress between layers, thereby decreasing interlayer friction and indirectly weakening the ability to counteract tensile stress, effectively increasing the net tensile stress. By finely decomposing and superimposing these stress components, the true stress level of the rock stratum under complex stress conditions can be comprehensively reflected, avoiding calculation errors caused by neglecting secondary stress components or coupling effects, and ensuring the comprehensiveness and accuracy of the total tensile stress calculation results.
[0034] Specifically, the formula for calculating the total tensile stress of the target rock stratum is: ,in, The tensile stress generated by the bending moment of the target rock stratum. This is the compressive stress generated by the gravitational component of the upper rock strata. This is the compressive stress caused by interlayer friction. The tensile stress generated by axial shear force, refer to Figure 3 and Figure 4 , Figure 3 This is a schematic diagram of the stress on a layered rock mass wading slope according to an embodiment of the tensile failure location prediction method of this application. Figure 4 This is a schematic diagram of target rock stress analysis according to an embodiment of the tensile fracture location prediction method of this application, as shown. Figure 3 and Figure 4 As shown, the first step is to perform a stress analysis on the target rock layer, including the force exerted by the upper rock layer (the lateral pressure exerted by the upper rock layer on the top of the target rock layer). pore water pressure Gravity load of the upper rock strata ,in, The lateral pressure coefficient, The rock is of high density. The distance from the top of the target rock layer to the top of the slope. Pore water pressure, It is water-weighted. The water head height, The width of the rock stratum.
[0035] The target rock layer length is set as The location of any point on the target rock stratum is ,make ,in, For position The remaining length, measured from the other end of the target rock layer, is used to simplify subsequent formulas. Therefore, any point on the target rock layer... The site was severely flooded Lateral pressure is: The gravitational force acting on it is: Axial shear force applied by the upper rock strata Components of normal force in the upper rock strata ; Interlayer friction generated under action ;in, The angle of inclination of the target rock stratum relative to the horizontal plane. for Axial component of force, for Axial component of force, for Normal component, for Normal component, is the coefficient of friction.
[0036] If we consider the target rock stratum as a cantilever beam, then the weight per unit length of the target rock stratum is: ,in, The compressive stress generated by the gravitational component of the upper rock layer is the thickness of the rock layer. The calculation formula is: Compressive stress caused by interlayer friction The calculation formula is: Tensile stress generated by axial shear force The calculation formula is: ,in, Moisture content, It represents the interlayer friction angle.
[0037] Further, in one embodiment, before obtaining the total tensile stress of the target rock layer by superimposing the tensile stress generated by the bending moment of the target rock layer, the compressive stress generated by the gravity component of the upper rock layer, the compressive stress caused by interlayer friction, and the tensile stress generated by the axial shear force, the following steps are included: Based on the bending moment at the fracture point of the target rock stratum and the bending strength of the target rock stratum, a cantilever plate beam model is used to calculate the tensile stress generated by the bending moment of the target rock stratum.
[0038] In this embodiment, the target rock strata in the anti-dip layered slope are treated as a cantilever beam model, conforming to its mechanical characteristics of being fixed at the root and free at the ends. The tensile stress generated by the bending moment can be calculated using mature material mechanics formulas based on this model. The cantilever model not only has clear physical meaning but also readily available parameters (such as the thickness and length of the target rock strata), making the calculation process both theoretically profound and engineeringally feasible. The tensile stress calculated by this model accurately reflects the stress concentration during the bending deformation process of the rock strata, providing a direct mechanical basis for determining whether the rock strata have undergone bending tensile cracking.
[0039] Continue to refer to Figure 4 Gravity components of the target rock strata and the forces caused by the upper rock strata It is the force that causes the target rock layer to bend and deform, located on the target rock layer. The bending moment at point is: The flexural stiffness of the target rock stratum is: Tensile stress generated by bending moment of target rock strata The calculation formula is: ,in, Gravity of the target rock strata The resulting bending moment components, Normal force component of the upper rock strata The resulting bending moment components, interlayer friction The resulting bending moment component (to resist bending moment).
[0040] Further, in one embodiment, before step S30, the following steps are included: Based on the fracture parameters, stress intensity factor, and geometric correction factor of the target rock stratum, the equivalent tensile strength of the fracture propagation in the target rock stratum is calculated using the equivalent tensile strength formula.
[0041] In this embodiment, primary or secondary fractures often exist within the rock mass, and the stress concentration effect at the fracture tips significantly reduces the actual bearing capacity of the rock mass. By introducing stress intensity factors and geometric correction factors from fracture mechanics, combined with parameters such as fracture half-length, the weakening effect of fractures on rock mass strength can be quantified. Through an equivalent tensile strength formula, the complex failure behavior of fractured rock masses is transformed into an equivalent strength index, enabling the prediction model to adapt to geological environments with fracture development. This approach overcomes the limitation of traditional continuum mechanics models in not considering the influence of discrete fractures, significantly improving the accuracy of tensile fracture location prediction in fracture-developed zones.
[0042] Specifically, the half-length of the fracture can be obtained through on-site laser scanning or core analysis. When the crack depth direction is perpendicular to the tensile stress, the stress intensity factor The calculation formula is: The fracture toughness of the target rock strata can be determined by a laboratory three-point bending test. According to the equivalent tensile strength formula, when When fracture occurs, the formula for calculating the stress intensity factor is as follows: ,make Solve by reverse reasoning The formula for calculating the equivalent tensile strength of the target rock stratum fracture propagation is as follows: ,in, The equivalent tensile strength for the propagation of fractures in the target rock strata. The nominal stress (far-field stress) applied to the rock mass. Geometric correction factor The method for determining the location of the crack is shown in Table 1. Table 1 shows the location of the crack. / Target rock layer length With geometric correction factor A table showing the correspondence between them, including the location of the crack. / Target rock layer length Location of the crack Length of the target rock layer The ratio, based on the location of the crack. Length of the target rock layer The ratio of the two values determines the geometric correction factor. The specific value.
[0043] Table 1. Location of cracks / Target rock layer length With geometric correction factor Correspondence table between
[0044] Furthermore, in one embodiment, before step S30, the method further includes: Based on the initial tensile strength, water content, and water-rock softening coefficient of the target rock stratum, the tensile strength of the target rock stratum under water-rock softening is calculated.
[0045] In this embodiment, fluctuations in reservoir water level lead to changes in rock mass water content, which in turn triggers a water-rock softening effect, causing a decrease in the tensile strength of the rock mass. By establishing a functional relationship between water content and tensile strength (such as an exponential decay model) and introducing water-rock softening coefficients for different lithologies, the rock mass strength under different hydrological conditions can be dynamically calculated. This allows the prediction method to adapt to different operating conditions such as reservoir impoundment and release, achieving dynamic early warning. Considering the water-rock softening effect can prevent the risk of rock mass strength reduction going undetected during the rainy season or high water level period, enhancing the sensitivity of the early warning system to environmental changes. The formula for calculating the tensile strength of the target rock stratum under water-rock softening is, for example: ,in, The initial tensile strength of the target rock stratum. Moisture content, This is the water softening coefficient (e.g., 0.15 for sandstone and 0.30 for mudstone).
[0046] Furthermore, in one embodiment, the stress balance equation is: ; in, Location on the target rock strata Total tensile stress at the point, To the critical stress at which tensile strength is broken, The equivalent tensile strength for the propagation of fractures in the target rock strata. The fracture toughness of the target rock stratum. For geometric correction factor, The fracture half-length of the target rock stratum, The tensile strength of the target rock strata softened by water. The initial tensile strength of the target rock stratum. The water-rock softening coefficient, This refers to the moisture content.
[0047] In this embodiment, the stress balance equation clearly defines the critical conditions for failure. The left side of the equation represents the stress state caused by external loads and geological conditions, while the right side represents the threshold of the rock mass's ability to resist failure. By clarifying the physical meaning and calculation formula of each parameter, the equation becomes solvable and verifiable. In particular, defining the critical stress as the minimum of the equivalent tensile strength and the water-softening tensile strength reflects the conservative design principle. This equation serves as a bridge between mechanical analysis and location prediction; its scientific construction directly determines the reliability of the final prediction results, providing a solid mathematical foundation for subsequent numerical solutions.
[0048] Further, in one embodiment, step S50 includes: Based on the stress balance equation, the position x on the target rock layer is solved using an iterative algorithm within the length L of the target rock layer to obtain the tensile fracture position x of the target rock layer.
[0049] In this embodiment, since the stress balance equation is usually nonlinear (containing polynomial or exponential terms), it is difficult to obtain an exact solution directly through analytical methods. An iterative algorithm (such as Newton's iteration method or the bisection method) is used to gradually approximate the root of the equation within the length L of the target rock layer, efficiently and accurately determining the position x that satisfies the equilibrium condition. The iterative algorithm has good convergence and numerical stability, making it suitable for implementation in engineering calculation software. The tensile fracture position x obtained through numerical solution is a specific distance value that can be directly mapped onto a geological profile, providing precise spatial coordinate guidance for engineering remediation and realizing the transformation from theoretical model to engineering application.
[0050] Secondly, embodiments of this application also provide a device for predicting the location of a breakage.
[0051] In one embodiment, reference is made to Figure 5 , Figure 5 This is a functional module diagram of an embodiment of the breakage position prediction device of this application, as shown below. Figure 5 As shown, the tensile strength location prediction device, applied to layered rock mass water-eroded slopes, includes: Select module 10, used to select target rock layers from layered rock mass wading slopes; Calculation module 20 is used to calculate the total tensile stress of the target rock layer based on its own weight, the force exerted by the upper rock layer, and the pore water pressure. Module 30 is used to select the minimum value between the equivalent tensile strength of the target rock layer fracture propagation and the tensile strength of water-rock softening as the critical stress for tensile fracture of the target rock layer. Module 40 is used to construct stress balance equations so that the total tensile stress of the target rock layer is equal to the critical stress at break. The solver module 50 is used to solve the stress balance equation using the length parameter of the target rock layer to obtain the fracture location of the target rock layer.
[0052] Furthermore, in one embodiment, the calculation module 20 is used for: The total tensile stress of the target rock layer is obtained by superimposing the tensile stress generated by the bending moment of the target rock layer, the compressive stress generated by the gravity component of the upper rock layer, the compressive stress caused by interlayer friction, and the tensile stress generated by the axial shear force. Among them, the compressive stress caused by interlayer friction is calculated based on the interlayer friction minus the pore water pressure.
[0053] Furthermore, in one embodiment, the tensile failure location prediction device further includes a bending moment and tensile stress calculation module, used for: Based on the bending moment at the fracture point of the target rock stratum and the bending strength of the target rock stratum, a cantilever plate beam model is used to calculate the tensile stress generated by the bending moment of the target rock stratum.
[0054] Furthermore, in one embodiment, the tensile fracture location prediction device further includes a crack propagation equivalent tensile strength calculation module, used for: Based on the fracture parameters, stress intensity factor, and geometric correction factor of the target rock stratum, the equivalent tensile strength of the fracture propagation in the target rock stratum is calculated using the equivalent tensile strength formula.
[0055] Furthermore, in one embodiment, the tensile fracture location prediction device further includes a water-rock softening tensile strength calculation module, used for: Based on the initial tensile strength, water content, and water-rock softening coefficient of the target rock stratum, the tensile strength of the target rock stratum under water-rock softening is calculated.
[0056] Furthermore, in one embodiment, the stress balance equation is: ; in, Let x be the total tensile stress at location x on the target rock stratum. To the critical stress at which tensile strength is broken, The equivalent tensile strength for the propagation of fractures in the target rock strata. Let Y be the fracture toughness of the target rock layer, Y be the geometric correction factor, and v be the fracture half-length of the target rock layer. The tensile strength of the target rock strata softened by water. k represents the initial tensile strength of the target rock layer. w ω is the water-rock softening coefficient, and ω is the water content.
[0057] Furthermore, in one embodiment, the solving module 50 is used for: Based on the stress balance equation, the position x on the target rock layer is solved using an iterative algorithm within the length L of the target rock layer to obtain the tensile fracture position x of the target rock layer.
[0058] The functions of each module in the above-mentioned tensile fracture location prediction device correspond to the steps in the above-mentioned tensile fracture location prediction method embodiment, and their functions and implementation processes will not be described in detail here.
[0059] Thirdly, embodiments of this application provide a device for predicting the location of a pull-off.
[0060] Reference Figure 6 , Figure 6 This is a schematic diagram of the hardware structure of the breakage position prediction device involved in the embodiments of this application. In the embodiments of this application, the breakage position prediction device may include a processor, a memory, a communication interface, and a communication bus.
[0061] The communication bus can be of any type and is used to interconnect the processor, memory, and communication interface.
[0062] The communication interface includes input / output (I / O) interfaces, physical interfaces, and logical interfaces used for interconnecting internal components of the pull-off position prediction device, as well as interfaces used for interconnecting the pull-off position prediction device with other devices (such as other computing devices or user equipment). Physical interfaces can be Ethernet interfaces, fiber optic interfaces, ATM interfaces, etc.; user equipment can be displays, keyboards, etc.
[0063] Memory can be various types of storage media, such as random access memory (RAM), read-only memory (ROM), non-volatile RAM (NVRAM), flash memory, optical storage, hard disk, programmable ROM (PROM), erasable PROM (EPROM), electrically erasable PROM (EEPROM), etc.
[0064] The processor can be a general-purpose processor, which can call the breakage position prediction program stored in memory and execute the breakage position prediction method provided in the embodiments of this application. For example, the general-purpose processor can be a central processing unit (CPU). The method executed when the breakage position prediction program is called can be referred to in various embodiments of the breakage position prediction method of this application, and will not be repeated here.
[0065] Those skilled in the art will understand that Figure 6 The hardware structure shown does not constitute a limitation of this application and may include more or fewer components than shown, or combine certain components, or have different component arrangements.
[0066] Fourthly, embodiments of this application also provide a readable storage medium.
[0067] The present application has a readable storage medium storing a breakage position prediction program, wherein when the breakage position prediction program is executed by a processor, it implements the steps of the breakage position prediction method as described above.
[0068] The method implemented when the breakage location prediction program is executed can be referred to in various embodiments of the breakage location prediction method of this application, and will not be repeated here.
[0069] It should be noted that the sequence numbers of the embodiments in this application are for descriptive purposes only and do not represent the superiority or inferiority of the embodiments.
[0070] The terms "comprising" and "having," and any variations thereof, in the specification, claims, and accompanying drawings of this application are intended to cover non-exclusive inclusion. For example, a process, method, system, product, or apparatus that includes a series of steps or units is not limited to the listed steps or units, but may optionally include steps or units not listed, or may optionally include other steps or units inherent to such process, method, product, or apparatus. The terms "first," "second," and "third," etc., are used to distinguish different objects, etc., and do not indicate a sequence, nor do they limit "first," "second," and "third" to different types.
[0071] In the description of the embodiments of this application, terms such as "exemplary," "for example," or "for instance" are used to indicate examples, illustrations, or explanations. Any embodiment or design described as "exemplary," "for example," or "for instance" in the embodiments of this application should not be construed as being more preferred or advantageous than other embodiments or designs. Specifically, the use of terms such as "exemplary," "for example," or "for instance" is intended to present the relevant concepts in a concrete manner.
[0072] In the description of the embodiments of this application, unless otherwise stated, " / " means "or". For example, A / B can mean A or B. The "and / or" in the text is merely a description of the relationship between related objects, indicating that there can be three relationships. For example, A and / or B can mean: A exists alone, A and B exist simultaneously, and B exists alone. In addition, in the description of the embodiments of this application, "multiple" means two or more.
[0073] In some processes described in the embodiments of this application, multiple operations or steps are included in a specific order. However, it should be understood that these operations or steps may not be executed in the order they appear in the embodiments of this application, or they may be executed in parallel. The sequence number of the operation is only used to distinguish different operations, and the sequence number itself does not represent any execution order. In addition, these processes may include more or fewer operations, and these operations or steps may be executed sequentially or in parallel, and these operations or steps may be combined.
[0074] Through the above description of the embodiments, those skilled in the art can clearly understand that the methods of the above embodiments can be implemented by means of software plus necessary general-purpose hardware platforms. Of course, they can also be implemented by hardware, but in many cases the former is a better implementation method. Based on this understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product is stored in a storage medium (such as ROM / RAM, magnetic disk, optical disk) as described above, and includes several instructions to cause a terminal device to execute the methods described in the various embodiments of this application.
[0075] The above are merely preferred embodiments of this application and do not limit the patent scope of this application. Any equivalent structural or procedural transformations made using the content of this application's specification and drawings, or direct or indirect applications in other related technical fields, are similarly included within the patent protection scope of this application.
Claims
1. A method for predicting the location of tensile failure, characterized in that, The method for predicting the location of tensile fractures, applied to layered rock mass wading slopes, includes: Select target rock layers from the layered rock mass wading slope; Calculate the total tensile stress of the target rock layer based on its own weight, the force exerted by the upper rock layer and the pore water pressure. The minimum value between the equivalent tensile strength of the target rock layer fracture propagation and the tensile strength of water-rock softening is selected as the critical stress for tensile fracture of the target rock layer. To make the total tensile stress of the target rock layer equal to the critical stress at break, a stress balance equation is constructed. By using the length parameter of the target rock layer, the stress balance equation is solved to obtain the fracture location of the target rock layer.
2. The method for predicting the breakage location as described in claim 1, characterized in that, The calculation of the total tensile stress of the target rock layer based on its own weight, the force exerted by the overlying rock layer, and the pore water pressure includes: The total tensile stress of the target rock layer is obtained by superimposing the tensile stress generated by the bending moment of the target rock layer, the compressive stress generated by the gravity component of the upper rock layer, the compressive stress caused by interlayer friction, and the tensile stress generated by the axial shear force. Among them, the compressive stress caused by interlayer friction is calculated based on the interlayer friction minus the pore water pressure.
3. The method for predicting the breakage location as described in claim 2, characterized in that, Before obtaining the total tensile stress of the target rock layer by superimposing the tensile stress generated by the bending moment of the target rock layer, the compressive stress generated by the gravity component of the upper rock layer, the compressive stress caused by interlayer friction, and the tensile stress generated by the axial shear force, the process includes: Based on the bending moment at the fracture point of the target rock stratum and the bending strength of the target rock stratum, a cantilever plate beam model is used to calculate the tensile stress generated by the bending moment of the target rock stratum.
4. The method for predicting the breakage location as described in claim 1, characterized in that, Before selecting the minimum value between the equivalent tensile strength of the target rock layer fracture propagation and the tensile strength of water-softening as the critical tensile stress for the target rock layer to break, the following steps are included: Based on the fracture parameters, stress intensity factor, and geometric correction factor of the target rock stratum, the equivalent tensile strength of the fracture propagation in the target rock stratum is calculated using the equivalent tensile strength formula.
5. The method for predicting the breakage location as described in claim 1, characterized in that, Before selecting the minimum value between the equivalent tensile strength of the target rock layer fracture propagation and the tensile strength of water-rock softening as the critical tensile stress for the target rock layer, the method further includes: Based on the initial tensile strength, water content, and water-rock softening coefficient of the target rock stratum, the tensile strength of the target rock stratum under water-rock softening is calculated.
6. The method for predicting the breakage location as described in claim 1, characterized in that, The stress balance equation is: ； in, Let x be the total tensile stress at location x on the target rock stratum. To the critical stress at which tensile strength is broken, The equivalent tensile strength for the propagation of fractures in the target rock strata. Let Y be the fracture toughness of the target rock layer, Y be the geometric correction factor, and v be the fracture half-length of the target rock layer. The tensile strength of the target rock strata softened by water. k represents the initial tensile strength of the target rock layer. w ω is the water-rock softening coefficient, and ω is the water content.
7. The method for predicting the breakage location as described in claim 6, characterized in that, The method of using the length parameter of the target rock layer to solve the stress balance equation to obtain the tensile fracture location of the target rock layer includes: Based on the stress balance equation, the position x on the target rock layer is solved using an iterative algorithm within the length L of the target rock layer to obtain the tensile fracture position x of the target rock layer.
8. A device for predicting the location of a tensile failure, characterized in that, The tensile failure location prediction device, applied to layered rock mass wading slopes, includes: The selection module is used to select target rock layers from a layered rock mass wading slope. The calculation module is used to calculate the total tensile stress of the target rock layer based on its own weight, the force exerted by the upper rock layer, and the pore water pressure. The determination module is used to select the minimum value between the equivalent tensile strength of the target rock layer fracture propagation and the tensile strength of water-rock softening as the critical stress for tensile fracture of the target rock layer. A module is constructed to make the total tensile stress of the target rock layer equal to the critical stress at break, and to construct the stress balance equation. The solver module is used to solve the stress balance equation using the length parameter of the target rock layer to obtain the fracture location of the target rock layer.
9. A device for predicting the location of a tensile failure, characterized in that, The breakage location prediction device includes a processor, a memory, and a breakage location prediction program stored in the memory and executable by the processor, wherein when the breakage location prediction program is executed by the processor, it implements the steps of the breakage location prediction method as described in any one of claims 1 to 7.
10. A readable storage medium, characterized in that, The readable storage medium stores a breakage location prediction program, wherein when the breakage location prediction program is executed by a processor, it implements the steps of the breakage location prediction method as described in any one of claims 1 to 7.