CN109506606B - Linear speed reciprocal method for predicting slope critical slip time - Google Patents
Linear speed reciprocal method for predicting slope critical slip time Download PDFInfo
- Publication number
- CN109506606B CN109506606B CN201811507753.6A CN201811507753A CN109506606B CN 109506606 B CN109506606 B CN 109506606B CN 201811507753 A CN201811507753 A CN 201811507753A CN 109506606 B CN109506606 B CN 109506606B
- Authority
- CN
- China
- Prior art keywords
- time
- speed
- reciprocal
- formula
- velocity
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
- G01B21/02—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness
Abstract
The invention discloses a linear speed reciprocal method for predicting the critical slip time of a side slope, which comprises the following steps: 1. selecting a concrete landslide body to be predicted, calculating the displacement speed of the side slope, and making a speed reciprocal-time scatter diagram in an accelerated sliding stage; 2. the formula of the invention is used for fitting the speed reciprocal-time scatter diagram to obtain the predicted slope critical slip time. The invention has the advantages that: the calculation of the critical sliding time has definite physical significance, and the invention can provide more accurate prediction.
Description
Technical Field
The invention belongs to the technical field of geological disaster prevention and control, and particularly relates to a linear velocity reciprocal method for predicting slope critical slip time.
Background
China is a multi-mountain country, the mountain area accounts for 2/3 of the land area of the country, the mountain area is complex in geological structure and fragile in ecological environment, mountain disasters frequently burst, and the loss caused by the mountain disasters accounts for one fourth of the total loss caused by various natural disasters. Landslide is a common geological disaster, and monitoring and early warning of the landslide are core tasks of disaster prevention and control work. Therefore, the method has great practical significance for predicting the landslide occurrence time in advance through the corresponding monitoring technology.
At present, the common landslide monitoring method is based on the monitoring of the surface displacement of the side slope, and the critical slide time of the side slope is predicted through a corresponding critical slide criterion. The speed reciprocal method is a prediction method which is widely applied and is relatively simple at present. The calculation formula of the speed reciprocal method is as follows:
wherein A is an empirical constant, tfFor landslide time, t*To accelerate the monitoring of the destruction phase, v*Is given as*The corresponding displacement speed.
In practical application, the speed of an accelerated sliding stage is obtained through displacement monitoring data calculation, and a speed reciprocal-time scatter diagram is drawn; then, assuming that the inverse velocity of the temporary slip is 0, the temporary slip time t is predicted by using the intersection point of the fitted curve of the inverse velocity-time scatter diagram and the time axisf。
According to the characteristics of the speed reciprocal-time curve in the accelerated sliding stage, landslides can be mainly divided into two types: 1. the straight line type is shown in FIG. 1, and the curve of the reciprocal velocity-time scatter point is a clear straight line; 2. the asymptotic profile is shown in FIG. 2, in which the divergence point of the velocity-time curve is an asymptote approaching the time axis.
This method has the following problems: 1. the physical meaning of the intersection point of the speed reciprocal-time curve and the time axis is not clear; 2. for a slope with a speed reciprocal-time scatter diagram of a gradual type, because a fitting curve does not have an intersection point with a time axis, the method cannot directly give the slide time of the slope, and needs to be specified through experience, so that the prediction error is large.
Disclosure of Invention
The invention aims to provide a linear speed reciprocal method for predicting the critical slip time of a side slope, which can enable the calculation of the critical slip time to have clear physical significance and can improve the prediction accuracy of the linear type slope.
The technical problem to be solved by the invention is realized by the technical scheme, which comprises the following steps:
step 1, selecting a concrete landslide body to be predicted, calculating original displacement monitoring data of the landslide body in an accelerated sliding stage to obtain the speed of a corresponding monitoring moment, and making a speed reciprocal-time scatter diagram in the accelerated stage;
if the scatter points are asymptotic towards the time axis in the reciprocal velocity-time plane, then the velocity is dynamically long averaged, or dynamically short averaged:
taking dynamic short average when n is 3 and dynamic long average when n is 7;
and 2, fitting the speed reciprocal-time diagram, wherein the fitting formula is as follows:
in the formula, t*To speed up the monitoring moment of the destruction phase, v*Is given as*Corresponding displacement speed, t0Is a certain moment v of the initial stage of the slope accelerated sliding0Is t0Sliding rate of time, tfThe critical slip time;
step 3, fitting the speed reciprocal-time scatter diagram obtained in the step 1 by using the formula in the step 2 to obtain the temporary slip time tf。
The invention has the technical effects that:
compared with the conventional speed reciprocal method, the method can directly obtain the critical slip time t through a fitting formula of speed reciprocal-time relation when predicting the critical slip time of the side slopefThe physical significance is more definite, and the prediction precision is high.
Drawings
The drawings of the invention are illustrated as follows:
FIG. 1 is a velocity-time diagram for a straight-line landslide;
FIG. 2 is a speed-time diagram of an asymptotic landslide;
FIG. 3 is a velocity reciprocal-time scatter plot of an embodiment;
FIG. 4 is a plot of inverse velocity versus time after the raw velocity dynamic long average process.
Detailed Description
The invention is further illustrated by the following examples in conjunction with the accompanying drawings:
in this embodiment, a Stromboli landslide is selected, and the following steps are performed:
step 1, selecting a concrete landslide body to be predicted, calculating original displacement monitoring data of the landslide body in an accelerated sliding stage to obtain the speed of a corresponding monitoring moment, and making a speed inverse-time scatter diagram in the accelerated stage according to a calculation formula of a speed inverse method; FIG. 3 is a velocity reciprocal-time scattergram of the present example.
If the scatter points are distributed asymptotically towards the time axis in the velocity reciprocal-time diagram, the velocity can be processed by a previous data processing, such as dynamic long average (LMA) or dynamic short average (SMA):
the patent takes the dynamic short average when n is 3 and the dynamic long average when n is 7; obtaining a linear distribution velocity reciprocal-time graph after the treatment; the velocity inverse-time graph after the velocity of fig. 3 is subjected to the dynamic long average processing is shown in fig. 4.
Step 2, fitting the speed reciprocal-time scatter point after the dynamic long average treatment to obtain a specific expression, and solving the critical slip time tf;
The empirical formula of speed and acceleration in accelerated creep phase of a landslide proposed by "Fukuzono T (1985) A new method for predicting the failure time of slopes, proceedings,4th International Conference & Field works on landslopes, Tokyo, pp 145-150" ("Fukuzono T (1985) slope failure time prediction, fourth International society for landslide Conference, Tokyo, Page number 145-150"):
v-α*a=A (2)
in the formula (2), v is the velocity of the slope, a is the acceleration of the slope, and α and a are empirical constants.
The rate at which the formula (2) is integrated is expressed as follows (α > 1):
in the formula, vfRate of landslide instability, tfThe critical slip time.
Let t be t at any time during the accelerated destruction phase*,ν=v*Converting the formula (3) to obtain:
assuming the speed v at which a landslide is temporarily slipperyfWhen α is 2, formula (4) reduces to:
let t0Is a certain moment v of the initial stage of the slope accelerated sliding0For the sliding rate at this time, there are:
from a comparison of equation (5) with equation (6) at both ends, we obtain:
fitting the velocity reciprocal-time scatter diagram obtained in the step 1 by using the formula (7) to obtain the critical slip time tf;
When the scatter points are linearly distributed in the reciprocal velocity-time diagram, that is, α is 2, the scatter points can be directly fitted to the scatter diagram according to the formula (7) by using a tool such as MATALAB, so as to obtain the predicted critical slip time tf。
Δt=t*-t0For monitoring the time t*And t0The time interval in between.
The temporary slip time t calculated by the inventionf201.5 hours differs from the actual landslide time of 205 hours by only 3.5 hours.
In order to further verify the prediction accuracy of the method, the temporary slip time of eight landslides is calculated through the steps and is compared with the real landslide time, and the specific calculation result is shown in a table 1:
TABLE 1
Name of landslide | Real landslide time/day | Predicting the critical slip time/day |
Switch station for terrace | 527 | 525 |
Temple with sleeping dragon | 66 | 67 |
Chicken temple | 466 | 463 |
White mountain Bierdanggu (Japan) | 113 | 112 |
Huang Zz | 544 | 540 |
La Chenaula | 884 | 886 |
Benishan | 258 | 258 |
Selburn | 600 | 599 |
As can be seen from table 1: the difference between the predicted critical sliding time and the real sliding time is very small, and the predicted critical sliding time of some cases is the same as the real sliding time. Therefore, the method can be used for predicting the landslide time in actual engineering, so that the life and property loss caused by landslide is reduced to the maximum extent.
Claims (1)
1. A linear speed reciprocal method for predicting the critical slip time of a side slope is characterized by comprising the following steps:
step 1, selecting a concrete landslide body to be predicted, calculating original displacement monitoring data of the landslide body in an accelerated sliding stage to obtain the speed of a corresponding monitoring moment, and making a speed reciprocal-time scatter diagram in the accelerated stage; if the scatter points are distributed in the reciprocal velocity-time plane in a manner of gradually approaching to the time axis, the velocity is subjected to dynamic long average or dynamic short average:
taking dynamic short average when n is 3 and dynamic long average when n is 7;
and 2, fitting the speed reciprocal-time diagram, wherein the fitting formula is as follows:
in the formula, t*To speed up the monitoring moment of the destruction phase, v*Is given as*Corresponding displacement speed, t0Is a certain moment v of the initial stage of the slope accelerated sliding0Is t0Sliding rate of time, tfThe critical slip time;
step 3, fitting the speed reciprocal-time scatter diagram obtained in the step 1 by using the formula in the step 2 to obtain the temporary slip time tf。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811507753.6A CN109506606B (en) | 2018-12-11 | 2018-12-11 | Linear speed reciprocal method for predicting slope critical slip time |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811507753.6A CN109506606B (en) | 2018-12-11 | 2018-12-11 | Linear speed reciprocal method for predicting slope critical slip time |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109506606A CN109506606A (en) | 2019-03-22 |
CN109506606B true CN109506606B (en) | 2020-10-09 |
Family
ID=65752048
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811507753.6A Active CN109506606B (en) | 2018-12-11 | 2018-12-11 | Linear speed reciprocal method for predicting slope critical slip time |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109506606B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111796113B (en) * | 2020-06-19 | 2022-03-04 | 西南交通建设集团股份有限公司 | Slope damage time determination method and system based on angular velocity reciprocal method |
CN113379138A (en) * | 2021-06-21 | 2021-09-10 | 重庆大学 | Inverse square root method for predicting slope damage time |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20080016057A (en) * | 2006-08-17 | 2008-02-21 | 정용호 | The coordinate displacement system which uses the length displacement sensor |
CN103424099A (en) * | 2013-08-26 | 2013-12-04 | 招商局重庆交通科研设计院有限公司 | Slope monitoring and early warning system and method based on deformation data |
CN103615962A (en) * | 2013-12-09 | 2014-03-05 | 三峡大学 | Landslide mass surface displacement measuring method |
CN103886176A (en) * | 2014-01-20 | 2014-06-25 | 青岛理工大学 | Method for determining dynamic stability coefficient of creep-type side slope |
-
2018
- 2018-12-11 CN CN201811507753.6A patent/CN109506606B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20080016057A (en) * | 2006-08-17 | 2008-02-21 | 정용호 | The coordinate displacement system which uses the length displacement sensor |
CN103424099A (en) * | 2013-08-26 | 2013-12-04 | 招商局重庆交通科研设计院有限公司 | Slope monitoring and early warning system and method based on deformation data |
CN103615962A (en) * | 2013-12-09 | 2014-03-05 | 三峡大学 | Landslide mass surface displacement measuring method |
CN103886176A (en) * | 2014-01-20 | 2014-06-25 | 青岛理工大学 | Method for determining dynamic stability coefficient of creep-type side slope |
Non-Patent Citations (2)
Title |
---|
速度倒数方法在不同运动形式滑坡中的应用;李德营等;《工程地质学报》;20141231;第299-303页第1-3节及图1-7 * |
速度倒数法滑坡预警模型流变试验研究;王延平等;《岩土力学》;20150630;第36卷(第6期);第1608页第3节及图2 * |
Also Published As
Publication number | Publication date |
---|---|
CN109506606A (en) | 2019-03-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109506606B (en) | Linear speed reciprocal method for predicting slope critical slip time | |
EP4300965A3 (en) | Interpolation for inter prediction with refinement | |
CN106406219B (en) | It is a kind of to exempt from programming electronic cam curve generation method for crosscutting | |
CN109543341B (en) | Power function speed reciprocal method for predicting side slope critical slip time | |
CN104614995A (en) | General design method for second-order system limited time slip form controller | |
CN108470101A (en) | Mechatronic Systems Y type sealing structure reliability estimation methods based on agent model | |
CN105574318A (en) | Oil well production automatic splitting device and method | |
AU2015416311A1 (en) | Automated upscaling of relative permeability using fractional flow in systems comprising disparate rock types | |
KR101632940B1 (en) | Method for identifying friction parameter for linear module | |
CN107704428B (en) | Bayesian resampling method for solving structural failure probability function | |
CN104359772B (en) | The normal stress threshold value determination method of rock mass discontinuity peak value angle of friction | |
CN103870656A (en) | Method for determining downburst crosswind profile | |
CN111159957A (en) | Method for rapidly predicting aerodynamic characteristics of bluff body shape transition flow region | |
CN108804829A (en) | Power hardening elastoplasticity bending crack opening displacement computational methods under a kind of dynamic load effect | |
CN112419493B (en) | Shale reservoir three-dimensional attribute model building method and device | |
CN104359771B (en) | The fractal dimension of a kind of rock mass discontinuity peak value angle of friction determines method | |
CN105094114B (en) | A kind of method of stability determining LPV controller | |
MX2021007046A (en) | Point cloud encoding method and encoder. | |
CN103575278B (en) | A kind of robust filtering method of losing for parameter uncertainty and observation | |
CN102322826A (en) | A kind of improved measuring method of object dimensional surface data | |
CN104973250B (en) | A kind of be applicable to glide vehicle press down trajectory angle of attack section defining method | |
Smith | Example of validation of statistical extrapolation | |
CN109919327A (en) | A kind of bearing maintenance opportunity acquisition methods | |
Ayas et al. | Energy release rate during the cracking of composite materials | |
CN111963143B (en) | Method for determining displacement boundary condition and load boundary condition of underground service casing |
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 |