CN112100859A - Multi-inversion method for shear strength parameters of slope rock-soil mass - Google Patents
Multi-inversion method for shear strength parameters of slope rock-soil mass Download PDFInfo
- Publication number
- CN112100859A CN112100859A CN202010986182.XA CN202010986182A CN112100859A CN 112100859 A CN112100859 A CN 112100859A CN 202010986182 A CN202010986182 A CN 202010986182A CN 112100859 A CN112100859 A CN 112100859A
- Authority
- CN
- China
- Prior art keywords
- internal friction
- slope
- rock stratum
- friction angle
- rock
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 239000002689 soil Substances 0.000 title claims abstract description 8
- 239000011435 rock Substances 0.000 claims abstract description 31
- 230000015572 biosynthetic process Effects 0.000 claims description 11
- 238000005065 mining Methods 0.000 abstract description 2
- 230000001360 synchronised effect Effects 0.000 abstract 1
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009877 rendering Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/23—Dune restoration or creation; Cliff stabilisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
Abstract
A multi-inversion method for shear strength parameters of slope rock-soil mass belongs to the field of open-pit mining. The method comprises the following steps: firstly, determining the cohesive force c and the volume weight gamma of each rock stratum of the slope; drawing up an internal friction angle of the rock stratum A, and continuously changing the internal friction angle of the rock stratum B until the integral stability coefficient of the slope is between 0.98 and 0.99, and marking the internal friction angle of the rock stratum B at the moment; thirdly, obtaining an internal friction angle of the rock stratum B with the integral stability coefficient of the side slope between 1.01 and 1.02 in a synchronous second step mode; fourthly, repeating the second step to obtain a plurality of groups of internal friction angle data of two rock strata which correspond one to one, and drawing a relation curve; fifthly, repeating the step three, similarly obtaining a group of internal friction angle data of two rock stratums which correspond one to one, and drawing a relation curve; and sixthly, the intersection points of the two relation curves are the finally determined internal friction angles of the rock stratum A and the rock stratum B respectively.
Description
Technical Field
The invention belongs to the field of surface mining, and particularly relates to a side slope rock-soil body shear strength parameter multiple inversion method.
Background
The shear strength parameter is data which is essential in the research of the slope problem, but the determination of the shear strength parameter is difficult, and the shear strength parameter is always widely concerned by geotechnical engineering researchers at home and abroad, and the reasonable shear strength parameter plays a key role in the evaluation of the slope stability and the landslide prevention and control engineering. The shear strength parameter inversion method is a common method for obtaining the data at present, the method is an inverse process of landslide stability calculation, the obtained parameters are more in line with the deformation condition of the landslide, and meanwhile, the method can also be used as a reference for experimental data selection, if no experimental data exists, the method can be directly applied to stability calculation, engineering design and the like, many successful examples at home and abroad are worth reference, but the traditional parameter inversion method is only suitable for the conditions of homogeneity, single sliding surface, landslide model and clear boundary conditions, and when some composite slope stability problems are analyzed, in order to ensure the reliability of selecting the shear strength parameters of the rock-soil mass, the shear strength parameters of the multilayer rock-soil mass need to be inverted at the same time, so that a method capable of inverting a plurality of shear strength parameters at the same time is urgently needed, and the method has great significance for analyzing the slope stability problems.
Disclosure of Invention
Aiming at the defects in the prior art, a residual thrust method in the rigid body limit balance theory is improved based on the rigid body limit balance theory, a novel inversion method for determining a plurality of shear strength parameters is provided, and a novel scientific means is provided for solving the problem of simultaneously determining a plurality of shear strength parameters.
In order to achieve the purpose, the technical scheme adopted by the invention is as follows:
in the residual thrust method, the slope is divided into a plurality of strips (figure 1), the nth strip is subjected to stress analysis (figure 2), and the stability coefficient of the slope is calculated according to the stress analysis result of the nth strip. However, the traditional inversion method can only establish an equation according to the existing parameters to solve one unknown shear strength parameter, but cannot simultaneously solve a plurality of unknown parameters. According to the method, the integral stability coefficient of the side slope and the local stability coefficient of the side slope in the limit balance state before the side slope is subjected to landslide are determined by analyzing the mechanical cause of the landslide, and two equations containing unknown shear strength parameters can be established according to the two stability coefficients and the existing parameters, so that the two unknown parameters can be inverted simultaneously.
In the shear strength parameters, the internal friction angle is difficult to measure in the experiment of rock mechanical parameters and plays a great role in the analysis of slope stability. The method can simultaneously invert two internal friction angles.
The invention discloses a slope rock-soil mass shear strength parameter multiple inversion method, which comprises the following steps:
the method comprises the following steps: determining the cohesive force c and the volume weight gamma of each rock stratum of the slope;
step two: the internal friction angle of the rock stratum AIs intended to be defined asContinuously changing the internal friction angle of the rock stratum B according to an inversion formula of the shear strength parameterUntil the internal friction angle of the rock stratum B is marked as the integral stability coefficient Fs of the slope is between 0.98 and 0.99The inversion formula of the shear strength parameter is as follows:
in the formula: wnThe weight (Kg) of the nth block is calculated by the volume weight gamma;
Fsis a stability factor;
anthe included angle (DEG) between the bottom edge of the nth block and the horizontal plane;
an-1the included angle (DEG) between the bottom edge of the (n-1) th block and the horizontal plane;
cnthe cohesion (kPa) of the nth block is equal to the rock formation penetrated by the sliding surfaceCohesion;
cn-1the cohesion (kPa) of the (n-1) th block is equal to the cohesion of the formation penetrated by the sliding surface;
lnthe length (m) of the bottom surface of the nth block;
ln-1the length (m) of the bottom surface of the (n-1) th block;
A=Wn-1sinan-1+Fn-2cos(an-2-an-1)
B=Wn-1cosan-1+Fn-2sin(an-2-an-1)
Wn-1the weight (Kg) of the (n-1) th block is calculated by the volume weight gamma;
Fn-2residual thrust (N) for the upper bar;
an-2the included angle (DEG) between the bottom edge of the (n-2) th block and the horizontal plane;
step three: using internal friction angle of rock stratum AContinuously changing the internal friction angle of the rock stratum B according to the inversion formula of the shear strength parameterUntil the integral stability factor Fs of the slope is between 1.01 and 1.02, the internal friction angle of the rock stratum B is marked
Step four: repeating the second step to obtain a plurality of groups which are in one-to-one correspondence and can meet the requirement that the integral stability coefficient Fs of the side slope is between 0.98 and 0.99Data is then plottedAndcurve of relationship between
Step five: repeating the third step to obtain a plurality of groups which can meet the one-to-one correspondence of the slope local stability coefficients Fs between 1.01 and 1.02Data is then plottedAndcurve of relationship between
Step six: curve lineAnd curveCorresponding to an intersection ofAndi.e., the finally determined internal friction angles of formation a and formation B, respectively (fig. 3).
From FIG. 3, the curvesAnd curveThe parameters at the intersection point respectively satisfy the integral stability coefficient and the local stability coefficient of the slope, so that the parameters obtained by inversion have higher reliability and can be used forThe slope stability analysis.
The invention has the beneficial effects that:
the invention is based on the rigid body limit balance theory, improves the residual reasoning method in the rigid body limit balance theory, provides a new inversion method for determining the shear strength parameters, can simultaneously determine a plurality of shear strength parameters by inversion under the condition of determining the landslide mechanical cause, provides a solving process, and makes up the defects of the traditional inversion method.
Drawings
FIG. 1 is a schematic diagram of slope calculation according to the present invention; and F is the horizontal direction load.
FIG. 2 is a force analysis chart of the nth block in the method of the present invention.
FIG. 3 is a graph of multiple inversions of internal friction angles of formation A and formation B in the method of the present invention.
FIG. 4 is a diagram illustrating an inversion result of an internal friction angle when an overall stability coefficient of a slope is between 0.98 and 0.99 according to an embodiment of the present invention;
FIG. 5 shows the inversion result of the internal friction angle when the local stability factor of the slope is between 1.01 and 1.02 in the embodiment of the present invention;
FIG. 6 is a multiple inversion curve of friction angles in mudstone and weak zones according to an embodiment of the present invention.
Detailed Description
Examples
And analyzing the slope stability and the landslide mechanism aiming at the slope where landslide occurs. The formation lithology includes: and in a third system, the mudstone, the coal and the weak layer, rock mass mechanical parameters of main rock stratums of the side slope can be determined by referring to the previous side slope stability research and geological results, but the mudstone and the weak layer rock mass mechanical parameters are not accurate.
According to the position of the collapse zone and the form of the weak layer, the inclination angle of the weak layer at the rear edge of the sliding body is larger, the front edge is smaller, and the rear edge rock mass extrudes the lower rock mass to slide along the weak layer under the action of self weight, so that the mechanical cause of the landslide is the pushing type. Therefore, the overall stability of the slope before landslide is smaller than the local stability of the slope at the front edge, and the basic principle of the parameter multiple inversion method is met. Therefore, the shear strength parameter multiple inversion method can be adopted to calculate the rock mass mechanical parameters of the mudstone and the weak layer, and the specific method comprises the following steps:
the method comprises the following steps: determination of mudstone cohesion C Mudstone26 kPa; weak layer cracking and weak layer cohesion CWeak layer=0kPa;
Step two: drawing up the internal friction angle of the weak layerContinuously changing the internal friction angle of the mudstone according to an inversion formula of the shear strength parameterUntil the overall stability factor F of the slopeSBetween 0.98 and 0.99, equal to 0.98, the internal friction angle of the mudstone is marked
Step three: using internal friction angle of weak layerContinuously changing the internal friction angle of the mudstone according to an inversion formula of the shear strength parameterUntil the local stability factor F of the slopeSBetween 1.01 and 1.02, equal to 1.02, the internal friction angle of the mudstone is marked
Step four: repeating the operation of the second step, and continuously drawing up a new internal friction angle of the weak layerObtaining a plurality of groups of one-to-one correspondence which can meet the condition that the integral stability factor Fs of the side slope is between 0.98 and 0.99The data, in particular,time of flight FS=0.983,Time of flightFS=0.984,Time of flightFS0.99; using multiple groups in step two and step fourData renderingAndcurve of relationship betweenAs shown in fig. 4;
step five: repeating the operation of the third step, and utilizing the new internal friction angle of the weak layerObtaining a plurality of groups of one-to-one correspondence which can meet the condition that the local stability factor Fs of the side slope is between 1.01 and 1.02The data, in particular,time of flightFS=1.19,Time of flightFS=1.016,Time of flightFS1.018; using the multiple groups in the third step and the fifth stepData renderingCurve of relationship betweenAs shown in fig. 5;
Claims (1)
1. A slope rock-soil mass shear strength parameter multiple inversion method is characterized by comprising the following steps:
the method comprises the following steps: determining the cohesive force c and the volume weight gamma of each rock stratum of the slope;
step two: the internal friction angle of the rock stratum AIs intended to be defined asContinuously changing the internal friction angle of the rock stratum B according to an inversion formula of the shear strength parameterUntil the internal friction angle mark of the rock stratum B is marked when the integral stability factor Fs of the side slope is between 0.98 and 0.99Is composed ofThe inversion formula of the shear strength parameter is as follows:
in the formula: wnThe weight (Kg) of the nth block is calculated by the volume weight gamma;
Fsis a stability factor;
anthe included angle (DEG) between the bottom edge of the nth block and the horizontal plane;
an-1the included angle (DEG) between the bottom edge of the (n-1) th block and the horizontal plane;
cnthe cohesion (kPa) of the nth piece is equal to the cohesion of the formation penetrated by the sliding surface;
cn-1the cohesion (kPa) of the (n-1) th block is equal to the cohesion of the formation penetrated by the sliding surface;
lnthe length (m) of the bottom surface of the nth block;
ln-1the length (m) of the bottom surface of the (n-1) th block;
A=Wn-1sinan-1+Fn-2cos(an-2-an-1)
B=Wn-1cosan-1+Fn-2sin(an-2-an-1)
Wn-1the weight (Kg) of the (n-1) th block is calculated by the volume weight gamma;
Fn-2residual thrust (N) for the upper bar;
an-2the included angle (DEG) between the bottom edge of the (n-2) th block and the horizontal plane;
step three: using internal friction angle of rock stratum AContinuously changing the internal friction angle of the rock stratum B according to the inversion formula of the shear strength parameterUntil the integral stability factor Fs of the slope is between 1.01 and 1.02, the internal friction angle of the rock stratum B is marked
Step four: repeating the second step to obtain a plurality of groups which are in one-to-one correspondence and can meet the requirement that the integral stability coefficient Fs of the side slope is between 0.98 and 0.99Data is then plottedAndcurve of relationship between
Step five: repeating the third step to obtain a plurality of groups which can meet the one-to-one correspondence of the slope local stability coefficients Fs between 1.01 and 1.02Data is then plottedAndcurve of relationship between
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010986182.XA CN112100859B (en) | 2020-09-18 | 2020-09-18 | Multiple inversion method for shear strength parameters of slope rock and soil mass |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010986182.XA CN112100859B (en) | 2020-09-18 | 2020-09-18 | Multiple inversion method for shear strength parameters of slope rock and soil mass |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112100859A true CN112100859A (en) | 2020-12-18 |
CN112100859B CN112100859B (en) | 2024-03-12 |
Family
ID=73759419
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010986182.XA Active CN112100859B (en) | 2020-09-18 | 2020-09-18 | Multiple inversion method for shear strength parameters of slope rock and soil mass |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112100859B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113408150A (en) * | 2021-07-16 | 2021-09-17 | 中南大学 | Parameter inversion method and system for geotechnical engineering random process |
CN116244813A (en) * | 2023-05-11 | 2023-06-09 | 中国铁路设计集团有限公司 | Soil body strength parameter inversion method based on sounding technology |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2034147C1 (en) * | 1992-05-15 | 1995-04-30 | Александр Яковлевич Егоров | Method for field determination of shear strength parameters of rocks in rock mass |
US20170284911A1 (en) * | 2016-03-31 | 2017-10-05 | Powerchina Huadong Engineering Corporation Limited | Integrated style shear apparatus for rock structural plane and a shear experimental method for rock structural plane |
CN109299562A (en) * | 2018-10-12 | 2019-02-01 | 中国电建集团西北勘测设计研究院有限公司 | A kind of ancient slide sliding surface Comprehensive Shear Strength metric history inversion method |
CN109726444A (en) * | 2018-12-07 | 2019-05-07 | 青岛理工大学 | Inversion determination method for rainfall type landslide shear strength parameter |
CN109800473A (en) * | 2018-12-26 | 2019-05-24 | 武汉大学 | Rock mechanical parameters inversion method based on differential evolution method |
-
2020
- 2020-09-18 CN CN202010986182.XA patent/CN112100859B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2034147C1 (en) * | 1992-05-15 | 1995-04-30 | Александр Яковлевич Егоров | Method for field determination of shear strength parameters of rocks in rock mass |
US20170284911A1 (en) * | 2016-03-31 | 2017-10-05 | Powerchina Huadong Engineering Corporation Limited | Integrated style shear apparatus for rock structural plane and a shear experimental method for rock structural plane |
CN109299562A (en) * | 2018-10-12 | 2019-02-01 | 中国电建集团西北勘测设计研究院有限公司 | A kind of ancient slide sliding surface Comprehensive Shear Strength metric history inversion method |
CN109726444A (en) * | 2018-12-07 | 2019-05-07 | 青岛理工大学 | Inversion determination method for rainfall type landslide shear strength parameter |
CN109800473A (en) * | 2018-12-26 | 2019-05-24 | 武汉大学 | Rock mechanical parameters inversion method based on differential evolution method |
Non-Patent Citations (1)
Title |
---|
曹兰柱;王珍;王东;宋子岭;白继元;: "扎哈淖尔露天煤矿南帮压煤回采边坡稳定性分析", 安全与环境学报, no. 01 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113408150A (en) * | 2021-07-16 | 2021-09-17 | 中南大学 | Parameter inversion method and system for geotechnical engineering random process |
CN116244813A (en) * | 2023-05-11 | 2023-06-09 | 中国铁路设计集团有限公司 | Soil body strength parameter inversion method based on sounding technology |
CN116244813B (en) * | 2023-05-11 | 2023-07-18 | 中国铁路设计集团有限公司 | Soil body strength parameter inversion method based on sounding technology |
Also Published As
Publication number | Publication date |
---|---|
CN112100859B (en) | 2024-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Faccenna et al. | Initiation of subduction at Atlantic‐type margins: Insights from laboratory experiments | |
CN103982179B (en) | A kind of paleopressure quantitative inversion detection method of oil reservoir | |
Cooke et al. | Bedding-plane slip in initial stages of fault-related folding | |
CN112100859A (en) | Multi-inversion method for shear strength parameters of slope rock-soil mass | |
Oldow et al. | Large-scale longitudinal extension in the southern Apennines contractional belt, Italy | |
Nunn et al. | Gravity anomalies and flexure of the lithosphere at the Middle Amazon Basin, Brazil | |
Gudmundsson | Fluid overpressure and flow in fault zones: field measurements and models | |
Nemčok | Rifts and passive margins: structural architecture, thermal regimes, and petroleum systems | |
CN103225509B (en) | Large-mining-height fully-mechanized mining stope roof classification and support resistance determination method | |
CN104141496A (en) | Rectangular roadway surrounding rock deformation and failure control method | |
CN102678118A (en) | Method for determining height of caving zone of stope | |
Duerto et al. | The role of syntectonic sedimentation in the evolution of doubly vergent thrust wedges and foreland folds | |
Bull et al. | Buckling of the oceanic lithosphere from geophysical data and experiments | |
CN113804617B (en) | Slope stability evaluation method under intermediate bridge retaining effect | |
Ansari et al. | Investigating Effect of Induced Stresses due to Coal Panel Extraction on Next Panel Strata behavior during Mechanized Longwall Mining: a Case Study | |
Hynes et al. | Subduction of continental margins and the uplift of high-pressure metamorphic rocks | |
Morley | Lateral and vertical changes of deformation style in the Osen-Røa thrust sheet, Oslo Region | |
Konon | Successive episodes of normal faulting and fracturing resulting from progressive extension during the uplift of the Holy Cross Mountains, Poland | |
Sapigni et al. | Engineering geological characterization and comparison of predicted and measured deformations of a cavern in the Italian Alps | |
Xu et al. | Seismic interpretation of the Kelasu triangle zone in the southern Tian Shan foothills, northwestern China | |
CN115201917A (en) | Method for predicting transverse wave velocity of sea-land transition phase shale formation | |
Wu et al. | Study on dynamics of tectonic evolution in the Fushun Basin, Northeast China | |
Haneberg | Drape folding of compressible elastic layers—I. Analytical solutions for vertical uplift | |
Vendeville | Scaled experiments on the interaction between salt flow and overburden faulting during syndepositional extension | |
CN115807664B (en) | Abnormal pressure prediction model construction method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |