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 PDFInfo
- 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
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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
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
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。
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)
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)
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 creeping 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 creeping slope |
Non-Patent Citations (2)
Title |
---|
李德营等: "速度倒数方法在不同运动形式滑坡中的应用", 《工程地质学报》 * |
王延平等: "速度倒数法滑坡预警模型流变试验研究", 《岩土力学》 * |
Cited By (3)
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 | |
Jakob et al. | Runout prediction methods | |
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 | |
CN110210621A (en) | A kind of object detection method based on residual error network improvement | |
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) | Determination method for stability coefficient and early warning criterion of creeping type artificial slope | |
CN111784536B (en) | Method for estimating water level overrun time in open channel according to actually measured water level change condition | |
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 | |
JP3821437B2 (en) | River condition simulation method | |
CN103870656A (en) | Method for determining downburst crosswind profile | |
CN109861201A (en) | A kind of new energy power output state computation method and system | |
Li et al. | On-line diagnosis method of crack behavior abnormality in concrete dams based on fluctuation of sequential parameter estimates | |
CN108897718A (en) | Quantitative description method for rock mass structure | |
CN109783760B (en) | Method for calculating initial speed of explosion and collapse fragments of reinforced concrete target plate | |
Haque et al. | Analytical study on steepness of ripples and dunes | |
Hou et al. | Vessel traffic flow distribution model of bridge area waterway in the middle stream of Yangtze River | |
Cruden et al. | A technique for estimating the complete creep curve of a sub-bituminous coal under uniaxial compression | |
CN113379138A (en) | Inverse square root method for predicting slope damage time | |
CN109826626A (en) | A kind of intelligent coalcutter cut mode identifying system | |
CN112133062B (en) | Landslide overall stability early warning method based on multi-monitoring-point synergistic effect | |
CN115983529B (en) | River channel overcurrent capacity analysis method, system and equipment based on geographic space recognition | |
Burkham et al. | General study of the modified Einstein method of computing total sediment discharge |
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 |