Rock slope stability evaluation and early warning method, system, equipment and medium

Abstract

The invention relates to the technical field of slope risk early warning based on slope data processing, in particular to a rock slope stability evaluation and early warning method, system, equipment and medium. The method comprises the following steps: firstly, obtaining an initial instability score and a slope stability type of a rock slope based on a gradient, rock-soil body shear strength and a structural surface development degree; then, considering the influence of potential slope instability in a rock slope group on the stability of a neighborhood slope, updating the initial instability score of the slope, and obtaining an initial comprehensive instability score of the slope; then, calculating an environmental instability score increment of the rock slope caused by rainfall, and updating the initial comprehensive instability score to obtain an updated comprehensive instability score of the rock slope; and finally, performing rock slope risk early warning according to a comparison result of the updated comprehensive instability score of the slope and an instability score threshold value. When the method is applied to the safety early warning of the rock slope, the accuracy and timeliness of the safety early warning of the rock slope can be improved.

Description

Technical Field

[0001] The present invention relates to the technical field of slope risk early warning based on slope data processing, and in particular to a rock slope stability assessment and early warning method, system, equipment and medium. Background Art

[0002] In the construction of transportation, water conservancy, mining and other projects, slope stability directly affects the safety and durability of the project. By evaluating slope stability, we can provide a safety design basis for the project design, discover potential dangers in advance, take preventive measures, and protect people's lives and property.

[0003] Traditional slope stability assessment methods are mostly based on static geological parameters and fail to fully consider the impact of dynamic environmental factors, such as rainfall, on slopes. In practical engineering, rainfall intensity and duration can significantly alter the seepage field and rock and soil strength of a slope. The current lack of an early warning mechanism to dynamically adjust the comprehensive stability value based on these environmental factors results in poor slope safety risk assessments. Summary of the Invention

[0004] The purpose of the present invention is to provide a rock slope stability assessment and early warning method, system, equipment and medium.

[0005] The technical solutions of the present invention are as follows:

[0006] A rock slope stability assessment and early warning method includes the following operations:

[0007] S1. Based on the slope, shear strength, and structural surface development of the rock slope, several instability contribution values ​​are obtained and constructed into a matrix to obtain an instability matrix. The weight matrix formed by the slope weight, shear strength weight, and structural surface development weight of the rock slope is multiplied by the instability matrix and normalized to obtain a slope classification assessment matrix. Based on the slope classification assessment matrix, the initial instability score and slope stability type of the rock slope are obtained.

[0008] S2. Based on the number of potential slopes within the neighborhood, and / or the sliding surface area, and / or the distance to the potential slope, obtain the forced instability value of each rock slope in the slope group, and update the initial instability score of each rock slope to obtain an initial comprehensive instability score for each rock slope;

[0009] S3. Obtain the fissure water filling degree and fissure water pressure increment of the rock slope within a preset time range based on the rainfall intensity and rainfall duration; obtain the environmental instability score increment of the rock slope based on the fissure water filling degree and the fissure water filling degree critical value; and obtain an updated comprehensive instability score of the slope based on the environmental instability score increment and the initial comprehensive instability score;

[0010] S4. When the updated comprehensive instability score of the slope exceeds the instability score threshold, a rock slope risk warning is issued.

[0011] The specific operation of obtaining the initial instability score in S1 is as follows: multiplying the slope classification assessment matrix with the instability weight matrix formed by the stable slope instability weight, the general stable slope instability weight and the unstable slope instability weight to obtain the initial instability value; and taking the product of the initial instability value and the total score as the initial instability score.

[0012] The specific method for obtaining the environmental instability score increment in S3 is as follows: determine whether the fissure water filling degree of the rock slope is less than the critical value of the fissure water filling degree; if it is less than, perform linear processing on the slope anti-sliding force attenuation value obtained based on the fissure water pressure increment to obtain the environmental instability score increment; if it is not less than, perform pressure superposition analysis on the fissure water filling degree, rainfall intensity and rainfall duration to obtain the environmental instability score increment.

[0013] The fracture water filling degree is calculated by the following formula:

[0014]

[0015] η is the water filling degree of the crack, Q n,f is the water filling capacity of the nth crack, Q n,f =I f ·t·A n,f , I f =min(q,C r ·r f ), I f is the effective rainfall infiltration rate, t is the rainfall duration within the preset range, A n,f is the crack area of ​​the nth crack, q is the rainfall intensity, C r is the fracture permeability, r f is the fracture permeability coefficient, V n,f is the crack volume of the nth crack, and N is the total number of cracks.

[0016] In S1, the method for obtaining the slope weight, rock and soil shear strength weight, and structural surface development weight of the rock slope is as follows: normalize the slope sample data, rock and soil sample data, and structural surface development sample data respectively to obtain the slope sample standard data, rock and soil sample standard data, and structural surface sample standard data; obtain the sample proportion of each standard data in the slope sample standard data, rock and soil sample standard data, and structural surface sample standard data respectively; based on the sample proportion of all standard data and the total number of sample data of the corresponding type, obtain the information entropy of each sample data to obtain the slope sample information entropy, rock and soil sample information entropy, and structural surface sample information entropy; based on the slope sample information entropy, rock and soil sample information entropy, and structural surface sample information entropy, obtain the slope weight, rock and soil shear strength weight, and structural surface development weight.

[0017] In S2, the forced instability value of a stable slope is obtained based on the number of potential danger slopes in the neighborhood; the forced instability value of a generally stable slope or unstable slope is obtained based on the number and sliding surface area of ​​potential danger slopes in the neighborhood, as well as the distance from the potential danger slopes.

[0018] The slope with hidden danger is an unstable slope whose sliding direction is toward a stable slope, a generally stable slope or an unstable slope.

[0019] A rock slope stability assessment and early warning system, used to implement the above-mentioned rock slope stability assessment and early warning method, comprises:

[0020] The module for generating an initial instability score and slope stability type is used to obtain a number of instability contribution values ​​based on the slope, shear strength of the rock soil mass, and the degree of structural surface development of the rock slope, and construct a matrix to obtain an instability matrix. The weight matrix formed by the slope weight of the rock slope, the shear strength weight of the rock soil mass, and the degree of structural surface development is multiplied by the instability matrix and normalized to obtain a slope classification assessment matrix. Based on the slope classification assessment matrix, the initial instability score and slope stability type of the rock slope are obtained.

[0021] An initial comprehensive instability score generation module is used to obtain the forced instability value of each rock slope in the slope group based on the number of potential slopes in the neighborhood, and / or the sliding surface area, and / or the distance to the potential slope, and to update the initial instability score of each rock slope to obtain the initial comprehensive instability score of each rock slope;

[0022] An updated comprehensive instability score generation module is used to obtain the fissure water filling degree and fissure water pressure increment of the rock slope within a preset time range based on rainfall intensity and rainfall duration; obtain the environmental instability score increment of the rock slope based on the fissure water filling degree and the fissure water filling degree critical value; and obtain an updated comprehensive instability score of the slope based on the environmental instability score increment and the initial comprehensive instability score;

[0023] The rock slope risk warning generation module is used to issue a rock slope risk warning when the slope's updated comprehensive instability score exceeds the instability score threshold.

[0024] A rock slope stability assessment and early warning device comprises a processor and a memory, wherein the processor implements the above-mentioned rock slope stability assessment and early warning method when executing a computer program stored in the memory.

[0025] A computer-readable storage medium is used to store a computer program, wherein when the computer program is executed by a processor, it implements the above-mentioned rock slope stability assessment and early warning method.

[0026] The beneficial effects of the present invention are:

[0027] The present invention provides a rock slope stability assessment and early warning method. First, based on the slope gradient, rock and soil shear strength, and structural surface development of the rock slope, an initial instability score of each rock slope is calculated; and based on the mapping relationship between the initial instability score of the slope and the slope stability type, the slope stability type of each slope is accurately assessed and obtained; then, considering the impact of the instability of the hidden danger slope within the rock slope group on the stability of the neighboring slopes, the forced instability values ​​of the slopes with different slope stability types are obtained, and the initial instability score of the slope is updated to obtain a more accurate initial comprehensive instability score of the slope; then, from the perspective that rainfall can cause an increase in the fissure water pressure of the rock slope, the environmental instability score increment of the rock slope caused by rainfall is calculated, and the initial comprehensive instability score is updated to obtain a more accurate updated comprehensive instability score of the rock slope; finally, based on the comparison result of the updated comprehensive instability score of the slope with the instability score threshold, a rock slope risk early warning is performed; the method is applied in rock slope safety early warning to improve the accuracy and timeliness of rock slope safety early warning. DETAILED DESCRIPTION

[0028] This embodiment provides a rock slope stability assessment and early warning method, including the following operations:

[0029] S1. Based on the slope, shear strength, and structural surface development of the rock slope, several instability contribution values ​​are obtained and constructed into a matrix to obtain an instability matrix. The weight matrix formed by the slope weight, shear strength weight, and structural surface development weight of the rock slope is multiplied by the instability matrix and normalized to obtain a slope classification assessment matrix. Based on the slope classification assessment matrix, the initial instability score and slope stability type of the rock slope are obtained.

[0030] S2. Based on the number of potential slopes within the neighborhood, and / or the sliding surface area, and / or the distance to the potential slope, obtain the forced instability value of each rock slope in the slope group, and update the initial instability score of each rock slope to obtain an initial comprehensive instability score for each rock slope;

[0031] S3. Obtain the fissure water filling degree and fissure water pressure increment of the rock slope within a preset time range based on the rainfall intensity and rainfall duration; obtain the environmental instability score increment of the rock slope based on the fissure water filling degree and the fissure water filling degree critical value; and obtain an updated comprehensive instability score of the slope based on the environmental instability score increment and the initial comprehensive instability score;

[0032] S4. When the updated comprehensive instability score of the slope exceeds the instability score threshold, a rock slope risk warning is issued. The specific operation details are as follows.

[0033] S1. Based on the slope of the rock slope, the shear strength of the rock soil, and the degree of structural surface development, the respective instability contribution values ​​are obtained and constructed into a matrix to obtain the instability matrix; the weight matrix formed by the slope weight of the rock slope, the shear strength weight of the rock soil, and the degree of structural surface development is multiplied by the instability matrix, and after normalization processing, a slope classification assessment matrix is ​​obtained; based on the slope classification assessment matrix, the initial instability score and slope stability type of the rock slope are obtained.

[0034] This embodiment calculates the initial instability score of each rock slope based on the slope gradient, rock and soil shear strength, and structural surface development; and accurately evaluates the slope stability type of each rock slope based on the mapping relationship between the initial instability score of the slope and the slope stability type.

[0035] Firstly, based on the slope of the rock slope, the shear strength of the rock and soil mass, and the degree of structural surface development, several instability contribution values ​​are obtained and constructed into a matrix to obtain the instability matrix.

[0036] The specific method for obtaining the instability contribution value is to substitute the slope, shear strength of the rock and soil mass, and the degree of structural surface development (structural surfaces include joints, fissures, and faults) of the rock slope into their respective instability segmentation evaluation functions to obtain the slope instability contribution value, the rock and soil mass instability contribution value, and the structural surface instability contribution value. The structural surface development degree is calculated based on the average number of structural surfaces within the region of interest of all slopes. A greater degree of structural surface development indicates a denser development of structural surfaces in the rock slope, a greater susceptibility to the presence of stress-weak zones, and a lower stability of the rock slope.

[0037] Since rock slopes are mainly composed of rocks, their stability is greatly affected by the slope (reflecting the distribution of gravitational potential energy and the magnitude of the sliding force of the rock slope), the development of structural surfaces (such as joints, cracks, and faults), and the type of rock and soil mass (reflecting the shear strength of the rock mass). Therefore, the unstable segmented evaluation function designed in this embodiment is as follows.

[0038] The unstable segmented evaluation function of slope or structural surface development is as follows:

[0039] When α≤α1, α α,1 =1,q α,2 =0,q α,3 =0;

[0040] When α1<α<α2, q α,3 =1;

[0041] When α≥α2, q α,1 =0,q α,2 =0,q α,3 =1;

[0042] In the unstable segmented evaluation function, there are multiple levels of stability evaluation values ​​in one segment of the evaluation function. Each level of stability evaluation value reflects the possibility that the rock slope belongs to each type of stable slope. θ is the slope or structural surface development degree, θ1 is the first threshold of the slope or structural surface development degree, θ2 is the second threshold of the slope or structural surface development degree, and q θ,1 is the first-level stability assessment value of the slope or structural surface development (the larger the first-level stability assessment value, the greater the possibility that the slope stability type belongs to the stable type), q θ,2 is the secondary stability assessment value of the slope or structural surface development (the larger the secondary stability assessment value, the greater the possibility that the slope stability type belongs to the general stability type), q θ,3 It is the third-level stability assessment value of the slope or structural surface development degree (the larger the third-level stability assessment value is, the greater the possibility that the slope stability type is unstable).

[0043] The unstable segmented evaluation function of the shear strength of rock and soil is as follows:

[0044] When τ≤τ1, qτ,1 =0,q τ,2 =0,q τ,3 =1;

[0045] When τ1<τ<τ2, q τ,1 =1,

[0046] When τ≥τ2, q τ,1 =1,q τ,2 =0,q τ,3 =0;

[0047] τ is the shear strength of rock mass, which reflects the anti-sliding ability of rock mass in natural state. It is based on the cohesion c and internal friction angle of rock slope. Got, σ is the normal stress of the structural surface. The greater the shear strength of the rock and soil, the better the slope stability. τ1 and τ2 are the first and second thresholds of the shear strength of the rock and soil, respectively. q τ,1 is the first-level stability assessment value of the shear strength of the rock and soil mass (the larger the first-level stability assessment value, the greater the possibility that the slope stability type belongs to the stable type), q τ,2 is the secondary stability assessment value of the shear strength of the rock and soil mass (the larger the secondary stability assessment value, the greater the possibility that the slope stability type belongs to the general stability type), q τ,3 It is the third-level stability assessment value of the shear strength of the rock and soil mass (the larger the third-level stability assessment value is, the greater the possibility that the slope stability type is unstable).

[0048] Therefore, for rock slopes, a matrix is ​​constructed from several instability contribution values, and the instability matrix is ​​obtained as follows:

[0049]

[0050] q 1,i is the stability assessment value of level i corresponding to the slope, q 2,i is the stability assessment value of level i corresponding to the shear strength of the rock and soil mass, q 3,i is the stability assessment value of level i corresponding to the degree of structural surface development, i = 1, 2, 3.

[0051] Then, the weight matrix formed by the rock slope's slope weight, the rock and soil shear strength weight, and the structural surface development weight is multiplied by the instability matrix and normalized to obtain a one-dimensional slope classification assessment matrix. In the slope classification assessment matrix, the first, second, and third data represent the rock slope's tendency to be stable, moderately stable, or unstable, respectively.

[0052] In the actual processing process, due to the irregular shape of the rock slope, the complexity of the geological composition and structural distribution, the slope, rock and soil mass, and structural surface development degree of the rock slope are usually taken as the average value of several sample data. At the same time, there are also differences in the quality of different types of sample data (for example, if the difference between samples in the same type of sample data is small, the sample data is evenly distributed and the sample quality is good). Taking this problem into consideration, this embodiment improves the objectivity of weight allocation and facilitates understanding the key risk factors of the rock slope through sample data. The slope weight, rock and soil mass shear strength weight, and structural surface development weight of the rock slope are designed. The weights are obtained by performing weight allocation processing on the slope sample data, rock and soil mass sample data, and structural surface development sample data based on data discretization analysis.

[0053] The specific steps of weight allocation processing based on data discretization analysis are as follows:

[0054] Step a: normalize the slope sample data, rock and soil sample data, and structural surface development sample data respectively to achieve sample data standardization, and obtain slope sample standard data, rock and soil sample standard data, and structural surface sample standard data.

[0055] Step b: Obtain the sample ratio of each standard data in the slope sample standard data, the rock and soil sample standard data, and the structural surface sample standard data respectively.

[0056] The sample proportion of standard data is calculated by the following formula:

[0057]

[0058] C' i,j is the sample proportion of the jth standard data in the i-th sample data, C i,j is the jth standard data in the i-th sample data, and J is the total number of standard data in the i-th sample data.

[0059] Step c: Based on the sample proportion of all standard data and the total number of sample data of the corresponding type, the information entropy of each sample data is obtained to obtain the slope sample information entropy, the rock and soil sample information entropy, and the structural surface sample information entropy.

[0060] The information entropy of each sample data is calculated using the following formula:

[0061]

[0062] E i is the information entropy of the i-th sample data, E1, E2, and E3 correspond to the information entropy of slope samples, the information entropy of rock and soil samples, and the information entropy of structural surface samples, respectively.

[0063] Step d: Based on the slope sample information entropy, the rock and soil sample information entropy, and the structural surface sample information entropy, the slope weight, the rock and soil shear strength weight, and the structural surface development weight are obtained.

[0064] The weight of each sample data is calculated using the following formula:

[0065]

[0066] w is the weight of the i-th sample data, I is the total number of sample data, I = 3, w1, w2, and w3 correspond to the slope weight, the shear strength weight of the rock and soil mass, and the weight of the structural surface development, respectively.

[0067] Finally, based on the slope classification assessment matrix, the initial instability score of the rock slope is obtained. Based on the mapping relationship between the initial instability score and the slope stability type, the slope stability type is obtained. The slope stability types include stable, generally stable, and unstable.

[0068] The above-mentioned operation to obtain the initial instability score is as follows: multiply the slope classification assessment matrix with the instability weight matrix (the instability weight matrix can be determined based on historical experience) formed by the instability weights of stable slopes, generally stable slopes, and unstable slopes. This quantifies the relationship between the inclination and instability of each type of slope and obtains the initial instability value. The product of the initial instability value and the total score is used as the initial instability score. For example, if the slope classification assessment matrix is ​​[0.64, 0.2, 0.1] and the instability weight matrix is ​​[0, 0.5, 1], the initial instability value is 0.2.

[0069] S2. Based on the number of potential slopes in the neighborhood, and / or the sliding surface area, and / or the distance to the potential slope, obtain the forced instability value of each rock slope in the slope group, update the initial instability score of each rock slope, and obtain the initial comprehensive instability score of each rock slope.

[0070] Considering the impact of the instability of hidden danger slopes within a rock slope group on the stability of neighboring slopes, the forced instability values ​​of slopes with different slope stability types are obtained, and then the initial instability scores of the slopes are updated to obtain a more accurate initial comprehensive instability score of the slopes.

[0071] First, considering that the distances between slopes in a slope group are relatively close, the instability of one slope may cause changes in the stresses of the surrounding slopes, affecting the stability of the surrounding slopes. Taking into account the mutual influence of such slope instability, this embodiment designs different methods based on the number of potential slopes in the neighborhood, and / or the sliding surface area, and / or the distance from the potential slope, according to the slope stability type, to obtain the forced instability value of each rock slope in the slope group. The specific methods are as follows.

[0072] When the rock slope is a stable slope, it is relatively stable and less sensitive to the instability disturbance of nearby slopes. If a large number of unstable slopes become unstable within the neighborhood, the stable slope will be affected and become unstable. Therefore, in this embodiment, the forced instability value of the stable slope is set based on the number of potential slopes with sliding surfaces sliding towards the stable slope within the neighborhood, which can be obtained by the following formula:

[0073] R1=1-e -βM ,

[0074] R1 is the forced instability value of the stable slope, M is the total number of potential danger slopes, and the potential danger slope is an unstable slope whose sliding direction is toward the stable slope. β is the quantity influence coefficient, which can be determined by statistical analysis of a large amount of slope data and using methods such as maximum likelihood estimation.

[0075] Since generally stable and unstable slopes are highly sensitive to the instability disturbance of nearby slopes and are easily affected by the external environment, the forced instability value of generally stable or unstable slopes is designed in this embodiment. It is obtained based on the number and sliding surface area of ​​potential danger slopes in the neighborhood and the distance from the potential danger slopes, which can be obtained by the following formula:

[0076]

[0077] R i is the forced instability value of the generally stable slope (i=2) or the unstable slope (i=3), M is the total number of potential danger slopes, and the potential danger slope is an unstable slope with the sliding direction of the sliding surface toward the generally stable slope or the unstable slope within the neighborhood of the generally stable slope or the unstable slope, γ is the area influence coefficient, S m is the sliding surface area of ​​the mth hidden danger slope, G(l m ) is the distance attenuation value of the mth hidden danger slope, which is the attenuation effect of the slope spacing on the degree of influence of the general stable slope or unstable slope neighborhood, ρ is the distance influence coefficient, l m is the distance between the mth hidden danger slope and the generally stable slope or unstable slope.

[0078] Then, according to the forced instability value of each slope in the slope group, the initial instability score of the corresponding slope is updated to obtain the initial comprehensive instability score of each slope. The initial comprehensive instability score can be calculated by the following formula: F'=W0F0+W B R i F max , F' is the initial comprehensive instability score, F0 is the initial instability score, F max is the total score, W0, WB are the initial scoring weight and the forced instability scoring weight, which respectively reflect the importance of the slope's own characteristics and the correction strength of the neighborhood influence, and are preferably 0.7 and 0.3 respectively. i is the forced instability value, i=1, 2, 3, corresponding to stable slope, generally stable slope, and unstable slope, respectively.

[0079] S3. Obtain the fissure water filling degree and fissure water pressure increment of the rock slope within a preset time range based on rainfall intensity and rainfall duration; obtain the environmental instability score increment of the rock slope based on the fissure water filling degree and the critical value of the fissure water filling degree; and obtain an updated comprehensive instability score of the slope based on the environmental instability score increment and the initial comprehensive instability score.

[0080] Rainfall not only directly scours the slope surface but also enters the cracks in the structural surfaces of the rock slope, changing the physical and mechanical properties of the rock and soil mass, increasing the slope's seepage pressure, and reducing its stability. Therefore, this embodiment calculates the incremental environmental instability score of the rock slope caused by rainfall from the perspective that rainfall can increase the water pressure in the rock slope's cracks, and updates the initial comprehensive instability score to obtain a more accurate updated comprehensive instability score for the rock slope.

[0081] First, based on the rainfall intensity and duration within a preset time range, the fracture water filling degree of the rock slope and the fracture water pressure increment reflecting the slope water pressure increment under the influence of rainfall are obtained.

[0082] The fracture water filling degree (the average value of the water filling degree of all fractures) is calculated by the following formula:

[0083]

[0084] η is the water filling degree of the crack, Q n,f In order to consider the effective crack filling capacity of the nth crack caused by rainfall intensity and rainfall duration, Q n,f =I f ·t·A n,f , I f =min(q,C r ·r f ), I f is the effective rainfall infiltration rate (unit: mm / h), t is the rainfall duration within the preset range, A n,f is the crack area of ​​the nth crack, is the average cross-sectional area of ​​the structural surface (including faults) whose crack width is greater than the width threshold, q is the rainfall intensity (the average rainfall intensity within the preset time range), C r is the fracture permeability (unit: mm / h), r f is the fracture permeability coefficient, the empirical value is 3600, Vn,f is the crack volume of the nth crack, and N is the total number of cracks.

[0085] The fracture water pressure increment (the sum of the water pressure increments in all fractures) is calculated using the following formula:

[0086]

[0087] ΔP is the increment of fracture water pressure, γ w is the bulk density of water, h n,f is the crack depth of the nth crack, θ n,f is the crack inclination angle of the nth crack.

[0088] Then, based on the fissure water filling degree and the critical value of the fissure water filling degree, the incremental environmental instability score of the rock slope is obtained. Specifically, the following steps are performed: determine whether the fissure water filling degree of the rock slope is less than the critical value of the fissure water filling degree; if it is less than the critical value, the internal pressure of the rock slope increases linearly and slowly. The slope anti-sliding force attenuation value obtained based on the fissure water pressure increment is linearly processed to obtain the incremental environmental instability score of the rock slope; if it is not less than the critical value, the internal pressure of the rock slope is in a state of rapid superposition, and the slope pressure effect increases suddenly. The fissure water filling degree, rainfall intensity, and rainfall duration are subjected to a pressure superposition analysis to obtain the incremental environmental instability score of the rock slope.

[0089] Under rainfall conditions, the environmental instability score increment of rock slopes is calculated using the following formula:

[0090] When η<η l hour,

[0091] When η≥η l When ΔF=min(F η,2 ,k η ·q·t·η);

[0092] η is the critical value of fracture water filling, ΔF is the increment of environmental instability score, is the attenuation value of the slope anti-sliding force obtained based on the Mohr-Coulomb principle, c is the cohesion, is the internal friction angle, σ n is the initial normal stress, which can be determined by the stress relief method of the rock and soil on site, F η,1 、F η,2 are the total score of the first increase in the fissure and the total score of the second increase in the fissure, respectively. η,1 =30, F η,2 =40, t is the duration of rainfall within the preset range, and q is the rainfall intensity (the average rainfall intensity within the preset time range).

[0093] Finally, based on the incremental environmental instability score and the initial comprehensive instability score, the updated comprehensive instability score of the rock slope is obtained. The updated comprehensive instability score is the sum of the incremental environmental instability score and the initial comprehensive instability score.

[0094] S4. When the updated comprehensive instability score of the slope exceeds the instability score threshold, a rock slope risk warning is issued.

[0095] This embodiment further provides a rock slope stability assessment and early warning system, which is used to implement the above-mentioned rock slope stability assessment and early warning method, including:

[0096] The module for generating an initial instability score and slope stability type is used to obtain a number of instability contribution values ​​based on the slope, shear strength of the rock soil mass, and the degree of structural surface development of the rock slope, and construct a matrix to obtain an instability matrix. The weight matrix formed by the slope weight of the rock slope, the shear strength weight of the rock soil mass, and the degree of structural surface development is multiplied by the instability matrix and normalized to obtain a slope classification assessment matrix. Based on the slope classification assessment matrix, the initial instability score and slope stability type of the rock slope are obtained.

[0097] An initial comprehensive instability score generation module is used to obtain the forced instability value of each rock slope in the slope group based on the number of potential slopes in the neighborhood, and / or the sliding surface area, and / or the distance to the potential slope, and to update the initial instability score of each rock slope to obtain the initial comprehensive instability score of each rock slope;

[0098] An updated comprehensive instability score generation module is used to obtain the fissure water filling degree and fissure water pressure increment of the rock slope within a preset time range based on rainfall intensity and rainfall duration; obtain the environmental instability score increment of the rock slope based on the fissure water filling degree and the fissure water filling degree critical value; and obtain an updated comprehensive instability score of the slope based on the environmental instability score increment and the initial comprehensive instability score;

[0099] The rock slope risk warning generation module is used to issue a rock slope risk warning when the slope updated comprehensive instability score exceeds the instability score threshold.

[0100] This embodiment also provides a rock slope stability assessment and early warning device, including a processor and a memory, wherein the processor implements the above-mentioned rock slope stability assessment and early warning method when executing a computer program stored in the memory.

[0101] This embodiment further provides a computer-readable storage medium for storing a computer program, wherein the computer program, when executed by a processor, implements the above-mentioned rock slope stability assessment and early warning method.

[0102] This embodiment provides a rock slope stability assessment and early warning method. First, based on the slope gradient, rock and soil shear strength, and structural surface development of the rock slope, an initial instability score is calculated for each rock slope. Then, based on the mapping relationship between the initial instability score and the slope stability type, the slope stability type of each slope is accurately assessed. Then, considering the impact of potential slope instability within a rock slope group on the stability of neighboring slopes, forced instability values ​​are obtained for slopes of different slope stability types. The initial instability score of the slope is then updated to obtain a more accurate initial comprehensive instability score. Next, based on the perspective that rainfall can cause an increase in fissure water pressure on the rock slope, the environmental instability score increment of the rock slope caused by rainfall is calculated, and the initial instability score is updated to obtain a more accurate updated comprehensive instability score of the rock slope. Finally, based on the comparison result of the updated comprehensive instability score with the instability score threshold, a rock slope risk early warning is issued. This method, when applied in rock slope safety early warning, can improve the accuracy and timeliness of rock slope safety early warning.

Claims

1. A rock slope stability assessment and early warning method, characterized in that: The following operations are included: S1. Based on the slope, shear strength, and structural surface development of the rock slope, several instability contribution values ​​are obtained and constructed into a matrix to obtain an instability matrix. The weight matrix formed by the slope weight, shear strength weight, and structural surface development weight of the rock slope is multiplied by the instability matrix and normalized to obtain a slope classification assessment matrix. Based on the slope classification assessment matrix, the initial instability score and slope stability type of the rock slope are obtained. S2. Based on the number of potential slopes within the neighborhood, and / or the sliding surface area, and / or the distance to the potential slope, obtain the forced instability value of each rock slope in the slope group, and update the initial instability score of each rock slope to obtain an initial comprehensive instability score for each rock slope; S3. Obtaining the fissure water filling degree and fissure water pressure increment of the rock slope within a preset time range based on the rainfall intensity and rainfall duration; and obtaining the environmental instability score increment of the rock slope based on the fissure water filling degree and the fissure water filling degree critical value; Based on the incremental environmental instability score and the initial comprehensive instability score, an updated comprehensive instability score of the slope is obtained; S4. When the updated comprehensive instability score of the slope exceeds the instability score threshold, a rock slope risk warning is issued.

2. A rock slope stability assessment and early warning method according to claim 1, characterized in that: The specific operations for obtaining the initial instability score in S1 are: The slope classification assessment matrix is ​​multiplied with the instability weight matrix formed by the instability weight of the stable slope, the instability weight of the generally stable slope, and the instability weight of the unstable slope to obtain the initial instability value; the product of the initial instability value and the total score is used as the initial instability score.

3. A rock slope stability assessment and early warning method according to claim 1, characterized in that: The specific method for obtaining the incremental environmental instability score in S3 is: Determine whether the fissure water filling degree of the rock slope is less than the critical value of the fissure water filling degree; if so, perform linear processing on the slope anti-sliding force attenuation value obtained based on the fissure water pressure increment to obtain the environmental instability score increment; If it is not less than, the fracture water filling degree, rainfall intensity and rainfall duration are subjected to pressure superposition analysis to obtain the incremental environmental instability score.

4. A rock slope stability assessment and early warning method according to claim 1 or 3, characterized in that: The fracture water filling degree is calculated by the following formula: η is the water filling degree of the crack, Q n,f is the water filling capacity of the nth crack, Q n,f =I f ·t·A n,f , I f =min(q,C r ·r f ), I f is the effective rainfall infiltration rate, t is the rainfall duration within the preset range, A n,f is the crack area of ​​the nth crack, q is the rainfall intensity, C r is the fracture permeability, r f is the fracture permeability coefficient, V n,f is the crack volume of the nth crack, and N is the total number of cracks.

5. A rock slope stability assessment and early warning method according to claim 1, characterized in that: The method for obtaining the slope weight of the rock slope, the shear strength weight of the rock and soil mass, and the weight of the structural surface development degree in S1 is as follows: Normalize the slope sample data, rock and soil sample data, and structural surface development sample data to obtain slope sample standard data, rock and soil sample standard data, and structural surface sample standard data; Obtain the sample ratio of each standard data in the slope sample standard data, rock and soil sample standard data, and structural surface sample standard data respectively; Based on the sample proportion of all standard data and the total number of sample data of the corresponding type, the information entropy of each sample data is obtained, and the information entropy of slope sample, rock and soil sample, and structural surface sample is obtained; Based on the information entropy of slope samples, rock and soil sample information entropy, and structural surface sample information entropy, the slope weight, rock and soil shear strength weight, and structural surface development weight are obtained.

6. A rock slope stability assessment and early warning method according to claim 1, characterized in that: In S2, the forced instability value of a stable slope is obtained based on the number of potential danger slopes in the neighborhood; the forced instability value of a generally stable slope or unstable slope is obtained based on the number and sliding surface area of ​​potential danger slopes in the neighborhood, as well as the distance from the potential danger slopes.

7. A rock slope stability assessment and early warning method according to claim 6, characterized in that: The slope with hidden danger is an unstable slope whose sliding direction is toward a stable slope, a generally stable slope or an unstable slope.

8. A rock slope stability assessment and early warning system, used to implement the rock slope stability assessment and early warning method according to claim 1, characterized in that: include: The module for generating an initial instability score and slope stability type is used to obtain a number of instability contribution values ​​based on the slope, shear strength of the rock soil mass, and the degree of structural surface development of the rock slope, and construct a matrix to obtain an instability matrix. The weight matrix formed by the slope weight of the rock slope, the shear strength weight of the rock soil mass, and the degree of structural surface development is multiplied by the instability matrix and normalized to obtain a slope classification assessment matrix. Based on the slope classification assessment matrix, the initial instability score and slope stability type of the rock slope are obtained. An initial comprehensive instability score generation module is used to obtain the forced instability value of each rock slope in the slope group based on the number of potential slopes in the neighborhood, and / or the sliding surface area, and / or the distance to the potential slope, and to update the initial instability score of each rock slope to obtain the initial comprehensive instability score of each rock slope; Updated the comprehensive instability score generation module to obtain the fracture water filling degree and fracture water pressure increment of the rock slope within a preset time range based on rainfall intensity and duration; and based on the fracture water filling degree and fracture water filling critical value, obtain the environmental instability score increment of the rock slope; Based on the incremental environmental instability score and the initial comprehensive instability score, an updated comprehensive instability score of the slope is obtained; The rock slope risk warning generation module is used to issue a rock slope risk warning when the slope updated comprehensive instability score exceeds the instability score threshold.

9. A rock slope stability assessment and early warning device, characterized in that: The method comprises a processor and a memory, wherein when the processor executes the computer program stored in the memory, the method for evaluating and warning the stability of a rock slope according to any one of claims 1 to 7 is implemented.

10. A computer-readable storage medium, characterized in that Used to store a computer program, wherein when the computer program is executed by a processor, a rock slope stability assessment and early warning method according to any one of claims 1 to 7 is implemented.

Images