CN109506606B - Linear speed reciprocal method for predicting slope critical slip time - Google Patents

Linear speed reciprocal method for predicting slope critical slip time Download PDF

Info

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
Application number
CN201811507753.6A
Other languages
Chinese (zh)
Other versions
CN109506606A (en
Inventor
周小平
刘林江
付梁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing University
Original Assignee
Chongqing University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chongqing University filed Critical Chongqing University
Priority to CN201811507753.6A priority Critical patent/CN109506606B/en
Publication of CN109506606A publication Critical patent/CN109506606A/en
Application granted granted Critical
Publication of CN109506606B publication Critical patent/CN109506606B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/02Measuring 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

Linear speed reciprocal method for predicting slope critical slip time
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:
Figure BDA0001899895990000011
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:
Figure BDA0001899895990000021
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:
Figure BDA0001899895990000022
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):
Figure BDA0001899895990000031
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):
Figure BDA0001899895990000032
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:
Figure BDA0001899895990000033
assuming the speed v at which a landslide is temporarily slipperyfWhen α is 2, formula (4) reduces to:
Figure BDA0001899895990000034
let t0Is a certain moment v of the initial stage of the slope accelerated sliding0For the sliding rate at this time, there are:
Figure BDA0001899895990000035
from a comparison of equation (5) with equation (6) at both ends, we obtain:
Figure BDA0001899895990000041
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:
Figure FDA0001899895980000011
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:
Figure FDA0001899895980000012
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
CN201811507753.6A 2018-12-11 2018-12-11 Linear speed reciprocal method for predicting slope critical slip time Active CN109506606B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (4)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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