CN109506606A - A kind of prediction side slope faces the linear speed counting backward technique of sliding time - Google Patents

A kind of prediction side slope faces the linear speed counting backward technique of sliding time Download PDF

Info

Publication number
CN109506606A
CN109506606A CN201811507753.6A CN201811507753A CN109506606A CN 109506606 A CN109506606 A CN 109506606A CN 201811507753 A CN201811507753 A CN 201811507753A CN 109506606 A CN109506606 A CN 109506606A
Authority
CN
China
Prior art keywords
time
speed
sliding
side slope
formula
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
Application number
CN201811507753.6A
Other languages
Chinese (zh)
Other versions
CN109506606B (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

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 the linear speed counting backward techniques that a kind of prediction side slope faces the sliding time, it calculates the velocity of displacement of side slope, make and accelerate sliding phase speed inverse-time scatter plot the following steps are included: 1, selected specific slip mass to be predicted；2, speed inverse-time scatter plot is fitted using formula of the invention, show that the side slope of prediction faces the sliding time.The invention has the advantages that making the calculating for facing the sliding time that there is specific physical significance, and the present invention can provide more accurate prediction.

Description

A kind of prediction side slope faces the linear speed counting backward technique of sliding time
Technical field
The invention belongs to geological disaster Control Technology fields, and in particular to a kind of prediction side slope faces the linear speed of sliding time Counting backward technique.
Background technique
China is the country on mountain more than one, and mountain area area accounts for the 2/3 of national land area, and mountain area geology construction is complicated, Ecological environment frailty, mountain region disaster are frequently broken out, and the loss caused by mountain region disaster is annual accounts for total caused by all kinds of natural calamities The a quarter of loss.Landslide is used as common geological disaster, and monitoring and warning is the core missions of diaster prevention and control work.So logical Corresponding monitoring technology is crossed, the time that look-ahead landslide occurs has great realistic meaning.
Currently, common landslide monitoring method is mostly based on the displacement monitoring of side slope surface, by face accordingly sliding criterion come Prediction side slope faces the sliding time.Speed counting backward technique is to be widely used at present and fairly simple prediction technique.Speed counting backward technique Calculation formula are as follows:
In formula, A is empirical, tfFor landslide time, t*For the monitoring moment in the stage that accelerates the failure, v*For with t*It is corresponding Velocity of displacement.
In practical applications, the speed for accelerating sliding phase is calculated by displacement monitoring data, draws speed and falls Number-time scatter plot；It is then assumed that speed inverse when facing sliding is 0, the matched curve of speed inverse-time scatter plot is utilized Predict to face sliding time t with the intersection point of time shaftf。
The characteristics of according to sliding phase speed inverse-time graph is accelerated, landslide can be mainly divided into two classes: 1, straight line Type is as shown in Figure 1, speed inverse-time scatterplot curve is in apparent straight line；2, asymptotic type as shown in Fig. 2, i.e. speed it is reciprocal- Time graph scatterplot is in the asymptote for approaching time shaft.
This method has the following problems: 1, the intersection point physical significance of speed inverse-time graph and time shaft is indefinite； It 2, is that the landslide of asymptotic type makes because matched curve and time shaft do not have intersection point for speed inverse-time scatter plot With this method, cannot directly give side slope faces the sliding time, needs to specify by experience, and prediction error is big.
Summary of the invention
The problem to be solved in the present invention is to provide the linear speed counting backward techniques that a kind of prediction side slope faces the sliding time, it can make The calculating for facing the sliding time has specific physical significance, and can improve the prediction accuracy on lienar for landslide.
The technical problem to be solved by the present invention is to technical solutions in this way to realize, it the following steps are included:
Step 1, selected specific slip mass to be predicted, are accelerating sliding phase, to its original displacement monitoring data into Row calculates, and is accordingly monitored the speed at moment, makes speed inverse-time scatter plot of boost phase；
If scatterplot is the asymptotic distribuion for tending to time shaft in speed inverse-time plane, it is long that dynamic is carried out to speed Average or dynamic is short average:
It is short average for dynamic when taking n=3, it is long average for dynamic when n=7；
Step 2 is fitted speed inverse-time diagram, the formula of the fitting are as follows:
In formula, t*For the monitoring moment in the stage that accelerates the failure, v*For with t*Corresponding velocity of displacement, t0Accelerate to slide for side slope The a certain moment of dynamic initial stage, v0For t0The gliding cable structure at moment, tfTo face the sliding time；
Step 3 is fitted and is faced to step 1 resulting speed inverse-time scatter plot using the formula of step 2 Sliding time tf。
The solution have the advantages that:
Compared with current velocity Huo existing speed counting backward technique, the present invention predict side slope face the sliding time when, can by speed inverse-when Between the fitting formula of relationship immediately arrive at and face sliding time tf, definitely, and precision of prediction is high for physical significance.
Detailed description of the invention
Detailed description of the invention of the invention is as follows:
Fig. 1 is the Velocity-time figure on linear type landslide；
Fig. 2 is the Velocity-time figure on asymptotic type landslide；
Fig. 3 is speed inverse-time scatter plot of embodiment；
Fig. 4 is speed inverse-time diagram after the long average treatment of raw velocity dynamic.
Specific embodiment
Present invention will be further explained below with reference to the attached drawings and examples:
The present embodiment selects Stromboli to come down, according to the following steps:
The selected specific slip mass to be predicted of step 1, is accelerating sliding phase, is carrying out to its original displacement monitoring data It calculates, is accordingly monitored the speed at moment, according to the calculation formula of speed counting backward technique, make the speed inverse-of boost phase Time scatter plot；Fig. 3 is the speed inverse-time scatter plot of the present embodiment.
If scatterplot is the asymptotic distribuion for tending to time shaft in speed inverse-time diagram, the data of early period can be passed through Processing, such as dynamic length averagely (LMA) or dynamic short average (SMA) are carried out to speed:
It is short average for dynamic when this patent takes n=3, it is long average for dynamic when n=7；Make linearly to be divided after processing Cloth speed inverse-time diagram；Speed inverse-time diagram such as Fig. 4 institute after carrying out the long average treatment of dynamic to the speed of Fig. 3 Show.
Step 2 is fitted speed inverse-time scatterplot after the long average treatment of dynamic, is specifically expressed Formula finds out and faces sliding time tf；
According to " Fukuzono T (1985) A new method for predicting the failure time of slopes.Proceedings,4th International Conference&Field Workshop on Landslides, A kind of Tokyo, pp145-150 " (" new method of Fukuzono T (1985) slope failure time prediction.4th international landslide Scientific seminar and procceedings, Tokyo, the page number 145-150 ") propose landslide tertiary creep stage velocity and acceleration warp Test formula:
v-α* (2) a=A
In formula (2), v is the rate of side slope, and a is the acceleration of side slope, and α, A are empirical.
Formula (2) is integrated and is expressed as follows (α > 1) with rate:
In formula, vfFor the rate of landslide failure, tfTo face the sliding time.
Assuming that in any time t=t in the stage that accelerates the failure*, ν=v*, transform (3):
Assuming that speed v when landslide faces slidingf=∞, as α=2, formula (4) simplifies are as follows:
Enable t0Accelerate a certain moment of sliding initial stage, v for side slope0For gliding cable structure at this time, then have:
By formula (5) compared with formula (6) both ends of the equation, obtain:
Step 1 resulting speed inverse-time scatter plot is fitted using formula (7), obtains facing sliding time tf；
It, can be directly using tools such as MATALAB when linearly distribution is α=2 to scatterplot in speed inverse-time diagram Scatter plot is fitted according to formula (7), so that is predicted faces sliding time tf。
Δ t=t*-t0To monitor moment t*With t0Between time interval.
Face sliding time t by what the present invention was calculatedf3.5 were only differed with practical landslide time 205 hours for 201.5 hours Hour.
For the precision of prediction for further verifying present method invention, calculate eight landslides through the above steps faces the sliding time And compared with true landslide time, specific calculated result is shown in Table 1:
Table 1
Come down title True landslide time/day / day is faced the sliding time in prediction
Its waste level ground switchyard 527 525
Temple of Recumbent Dragon 66 67
Crow temple 466 463
White mountain is not when paddy (Japan) 113 112
Yellow thatch 544 540
La Chenaula 884 886
Benny mountain 258 258
Selborne 600 599
As can be seen from Table 1: the present invention predicts that facing sliding time and the difference of true landslide time is very small, Er Qieqi In the predictions of some cases to face the sliding time identical as true landslide time.Therefore this method can be used in actual engineering Landslide time is predicted in invention, to reduce life and property loss caused by landslide to the maximum extent.

Claims (1)

1. the linear speed counting backward technique that a kind of prediction side slope faces the sliding time, characterized in that the following steps are included:
Step 1, selected specific slip mass to be predicted, are accelerating sliding phase, original displacement monitoring data are counted to its It calculates, is accordingly monitored the speed at moment, make speed inverse-time scatter plot of boost phase；If scatterplot is fallen in speed It is the distribution asymptotic to time shaft in number-time plane, then dynamic length is carried out averagely to speed or dynamic is short average:
It is short average for dynamic when taking n=3, it is long average for dynamic when n=7；
Step 2 is fitted speed inverse-time diagram, the formula of the fitting are as follows:
In formula, t*For the monitoring moment in the stage that accelerates the failure, v*For with t*Corresponding velocity of displacement, t0Accelerate sliding initial for side slope The a certain moment in stage, v0For t0The gliding cable structure at moment, tfTo face the sliding time；
Step 3 is fitted to obtain when facing sliding using the formula of step 2 to step 1 resulting speed inverse-time scatter plot Between 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 true CN109506606A (en) 2019-03-22
CN109506606B 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)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111796113A (en) * 2020-06-19 2020-10-20 西南交通建设集团股份有限公司 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
李德营等: "速度倒数方法在不同运动形式滑坡中的应用", 《工程地质学报》 *
王延平等: "速度倒数法滑坡预警模型流变试验研究", 《岩土力学》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111796113A (en) * 2020-06-19 2020-10-20 西南交通建设集团股份有限公司 Slope damage time determination method and system based on angular velocity reciprocal method
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

Also Published As

Publication number Publication date
CN109506606B (en) 2020-10-09

Similar Documents

Publication Publication Date Title
Tighe et al. Model for the scaling of stresses and fluctuations in flows near jamming
CN108538026A (en) A kind of multiple attribute synthetical evaluation and method for early warning of slope stability
CN111861107B (en) Landslide dynamic stability evaluation method based on intensity reduction method
CN109506606A (en) A kind of prediction side slope faces the linear speed counting backward technique of sliding time
CN109543341A (en) A kind of prediction side slope faces the power function speed counting backward technique of sliding time
CN109680669B (en) Foundation pit stability monitoring and early warning method based on underground water dynamic loading
CN103646181A (en) Method for determining wriggle slide type artificial side slope stability coefficient and early warning criteria
CN103605841A (en) Method for building numerical simulation model of hydraulic engineering asphalt concrete
CN111784536B (en) Method for estimating water level overrun time in open channel according to actually measured water level change condition
Duncan et al. Evaluating reliability in geotechnical engineering
True Penetration of projectiles into seafloor soils
CN116227941B (en) Risk simulation calculation evaluation method and system for water diversion project
CN111931429B (en) Simulation method for change of effective density of flocs along with particle size
CN103870656A (en) Method for determining downburst crosswind profile
CN109861201A (en) A kind of new energy power output state computation method and system
CN112763181B (en) Method for determining sampling parameters of pulsating pressure wind tunnel test signals
CN109783760B (en) Method for calculating initial speed of explosion and collapse fragments of reinforced concrete target plate
CN108874735B (en) Method and device for determining pressure of ancient fluid in sedimentary basin
Bechteler et al. Stochastic model for particle movement in turbulent open channel flow
CN110427655A (en) A kind of extracting method for the sensitiveness that comes down
Hou et al. Vessel traffic flow distribution model of bridge area waterway in the middle stream of Yangtze River
CN113379138A (en) Inverse square root method for predicting slope damage time
CN109826626A (en) A kind of intelligent coalcutter cut mode identifying system
KR102102957B1 (en) Data collecting system for evaluating deteriorated level of facility, and method for the same
CN112133062B (en) Landslide overall stability early warning method based on multi-monitoring-point synergistic effect

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