CN115563802A - Landslide and critical-sliding-time forecasting method based on self-adaptive time window and multiple criteria - Google Patents

Landslide and critical-sliding-time forecasting method based on self-adaptive time window and multiple criteria Download PDF

Info

Publication number
CN115563802A
CN115563802A CN202211304893.XA CN202211304893A CN115563802A CN 115563802 A CN115563802 A CN 115563802A CN 202211304893 A CN202211304893 A CN 202211304893A CN 115563802 A CN115563802 A CN 115563802A
Authority
CN
China
Prior art keywords
time
deformation
time window
landslide
tangent angle
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.)
Pending
Application number
CN202211304893.XA
Other languages
Chinese (zh)
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.)
Beijing Institute of Technology BIT
Chongqing Innovation Center of Beijing University of Technology
Original Assignee
Beijing Institute of Technology BIT
Chongqing Innovation Center of Beijing University of Technology
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 Beijing Institute of Technology BIT, Chongqing Innovation Center of Beijing University of Technology filed Critical Beijing Institute of Technology BIT
Priority to CN202211304893.XA priority Critical patent/CN115563802A/en
Publication of CN115563802A publication Critical patent/CN115563802A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Abstract

The invention discloses a landslide impending slip time forecasting method based on a self-adaptive time window and multiple criteria. According to the method, a plurality of time windows are set, and the optimal time window is selected based on the tangent angle time curve, so that the tangent angle calculation error is avoided, and the intelligence and the accuracy of landslide time prediction are improved; and then, the Pearson correlation coefficient is used as a second criterion for identifying the critical slip stage, namely, when the tangent angle value is large enough, the landslide time is predicted by adopting a speed inverse method, and when the Pearson correlation coefficient reaches a certain range, a landslide time prediction result is input. The method adaptively determines the optimal time window, and simultaneously adopts double criteria to ensure the accuracy of the identification in the critical slip stage, thereby improving the reliability of the forecast.

Description

Landslide and critical-sliding-time forecasting method based on self-adaptive time window and multiple criteria
Technical Field
The invention relates to the technical field of alarms responding to disaster events, in particular to a landslide impending landslide time forecasting method based on an adaptive time window and multiple criteria.
Background
Based on the landslide three-stage deformation theory of M.Saito, T.A.Fukuzono proves that the reciprocal speed of the landslide time-to-slip has a linear relation with time through analytical experiments. When the landslide deformation speed is increased, the expected landslide occurrence time can be calculated according to a speed reciprocal method, but the time forecast has reliability only when the landslide is in a near-slip stage. At present, the mainstream forecasting method is as follows: and judging the critical sliding stage of the landslide by adopting the tangent angle value, and forecasting the time by utilizing a speed reciprocal method. This method has several drawbacks: 1) The time window for calculating the tangent angle and the tangent speed is not definite, some researchers adopt the sampling time interval of monitoring data points, some researchers adopt a fixed time window (such as 3h,7 h), the tangent angle value has higher sensitivity to the time window, and an improper time window easily causes the calculation error of the tangent angle; 2) Only depending on a criterion of the tangent angle value, the accuracy of the identification in the critical sliding stage is not enough.
Disclosure of Invention
In view of this, the invention provides a landslide imminent-slip time forecasting method based on a self-adaptive time window and multiple criteria, which can effectively improve the reliability of landslide time forecasting.
The invention discloses a landslide impending slide time forecasting method based on a self-adaptive time window and multiple criteria, which comprises the following steps:
s1, monitoring a scene to acquire a deformation time curve; self-adaptively obtaining an optimal time window; the method for acquiring the optimal time window specifically comprises the following steps:
s11, continuously monitoring a scene in real time to obtain a deformation time curve; setting a plurality of fixed time windows, respectively calculating tangent angles of the deformation time curves under the fixed time windows, and moving the fixed time windows to obtain tangent angle curves corresponding to the fixed time windows;
s12, sequencing the fixed time windows from small to large, and sequentially calculating the coincidence rate of tangent angle curves under the adjacent fixed time windows; when the coincidence rate exceeds a set threshold value B for the first time, the corresponding fixed time window is the optimal time window;
s2, calculating tangent angles of all moments based on the optimal time window obtained in the step S1;
s3, calculating the reciprocal of the speed at each moment based on the optimal time window obtained in the step S1, and further obtaining the Pearson correlation coefficient at each moment;
and S4, extracting the moments when the tangent angle is larger than the threshold theta and the Pearson correlation coefficient is smaller than the threshold beta, and outputting the landslide critical slip time predicted by adopting an inverse speed method at the moments.
Preferably, the tangent angle under the fixed time window Δ t is calculated as follows:
first, the accumulated deformation S is calculated i Average rate of (B) i
Figure BDA0003905419020000021
Wherein S is 1 、S i Respectively, t in the deformation time curve 1 Time t and i accumulated deformation corresponding to the moment;
then, using the accumulated deformation S i Divided by the average rate B i Obtaining T consistent with the dimension of time T i
Figure BDA0003905419020000022
Then t is i Time of day, deformation time curve at fixed time window Δ t (Δ t = t) i -t i-n ) Lower tangent angle alpha i Comprises the following steps:
Figure BDA0003905419020000023
preferably, in step S12, the method for calculating the overlapping ratio of the tangent angle curves in the adjacent fixed time windows is as follows:
calculating tangent angle difference values on tangent angle curves corresponding to adjacent time windows at the same moment, and when the difference value is less than or equal to a set threshold value A, considering that the information of the two tangent angle curves at the moment is overlapped; and the ratio of the information coincidence duration of the two tangent angle curves to the total duration is the coincidence rate of the two curves.
Preferably, in S3, the pearson correlation coefficient calculation method is as follows:
calculating the inverse velocity IV in a fixed time window deltat i
Figure BDA0003905419020000031
In the formula (7), Δ t is the optimal time window obtained by S1; s i 、S i-n Respectively, t in the deformation time curve i Time t and i-n accumulated deformation corresponding to the moment;
then t i Pearson's correlation coefficient r at time i Comprises the following steps:
Figure BDA0003905419020000032
wherein IV i-n ~IV i Are each t i-n ~t i The inverse of the velocity at that moment.
Preferably, in S3, the speed reciprocal at each time is filtered to remove the maximum value and the negative value.
Preferably, the threshold θ depends on the type of landslide, θ ∈ [60 °,85 ° ].
Preferably, the threshold β depends on the type of landslide, β ∈ [ -0.90, -0.70].
Preferably, in S1, a deformation time curve is obtained by using a ground GNSS device or a ground-based interferometric radar.
Preferably, in the step S1, a ground-based interference radar is used to obtain an accumulated deformation amount of a surface of the monitoring scene; determining a tangent angle based on the surface accumulated deformation, specifically comprising:
s101, designing a plurality of deformation thresholds according to the maximum value in the surface accumulated deformation quantity;
s102, acquiring an area-time curve corresponding to each deformation threshold:
aiming at each deformation threshold, calculating the total area of the pixel points of which the accumulated deformation exceeds the deformation threshold at the current moment in real time; each moment corresponds to a total area value of the pixel points exceeding the deformation threshold, and an area-time curve is formed along with the lapse of time;
and S103, calculating tangent angles of the area-time curves corresponding to the deformation thresholds at the current moment aiming at the fixed time windows based on the self-adaptive time window method, wherein the largest tangent angle of the area-time curves corresponding to the deformation thresholds is the tangent angle at the current moment corresponding to the fixed time window.
Preferably, in S101, the N deformation thresholds are designed in the following manner:
maximum accumulated deformation value S in monitoring scene of current measurement period max When the thickness is less than or equal to 100mm,
Figure BDA0003905419020000041
deformation threshold S k =10·k,k=1~N,
Figure BDA0003905419020000042
Represents rounding down;
maximum accumulated deformation value S in monitoring scene of current measurement period max When the thickness is larger than 100mm,
Figure BDA0003905419020000043
deformation threshold
Figure BDA0003905419020000044
Has the beneficial effects that:
(1) According to the method, a plurality of time windows are set, the optimal time window is selected based on the tangent angle time curve, the tangent angle calculation error is avoided, and the intelligence and the accuracy of landslide time prediction are improved; and then, the Pearson correlation coefficient is used as a second criterion for identifying the critical slip stage, namely, when the tangent angle value is large enough, the landslide time is predicted by adopting an inverse speed method, and when the Pearson correlation coefficient reaches a certain range, a landslide time prediction result is input. The method adaptively determines the optimal time window, and simultaneously adopts double criteria to ensure the identification accuracy of the critical slip stage, thereby improving the reliability of forecasting.
(2) And the different tangent angle curves judge whether the tangent angles are overlapped or not according to the distances of the corresponding tangent angles at the same moment, so that the method is simple and the operability is strong.
(3) When the ground interference radar is adopted to obtain the surface accumulated deformation amount of the monitoring scene, the tangent angle can be determined based on the surface accumulated deformation amount, compared with the traditional early warning method based on single-point deformation information, the method adopting the area time sequence curve can identify the whole deformation trend of the landslide surface, the early warning information is more comprehensive and accurate, and the accuracy of judgment in the landslide early warning stage can be effectively improved.
Drawings
FIG. 1 is a flow chart of the present invention.
Fig. 2 is a plot of tangent angle α for multiple time windows.
Fig. 3 shows the coincidence of the tangent angle alpha curves for multiple time windows.
Fig. 4 is a plot of tangent angle α under an adaptive time window.
Fig. 5 shows the reciprocal velocity values at various times.
Fig. 6 shows the reciprocal velocity values (after filtering) at each time.
Fig. 7 shows the pearson correlation coefficient at each time.
FIG. 8 is t i Time forecast at time instant =51.33 h.
FIG. 9 is a landslide time forecast based on an adaptive time window and multiple criteria.
FIG. 10 is a flow chart of adaptive time window acquisition,
Fig. 11 is a tangent angle acquisition method based on a plane time curve.
Detailed Description
The invention is described in detail below by way of example with reference to the accompanying drawings.
The invention provides a landslide impending slip time forecasting method based on a self-adaptive time window and multiple criteria, which comprises the following steps of:
1. adaptive time window acquisition
When the tangent angle is calculated based on the deformation time curve, the calculated value of the tangent angle has high sensitivity to the time window, if the undersize time window easily causes the jump of the tangent angle, the oversize time window can cause the delay of the early warning information, and the fixed time window suitable for all conditions does not exist for deformation measuring equipment with different monitoring frequencies and deformation characteristics of different types of disasters. Therefore, this step first adaptively selects an optimal time window from a plurality of time windows.
First, the accumulated deformation S is calculated i Average rate of (B) i
Figure BDA0003905419020000051
Wherein S is i 、S 1 Are each t i Time and t 1 The accumulated deformation amount corresponding to the moment.
Then, using the accumulated deformation S i Divided by the average rate B i Obtaining T consistent with the dimension of time T i
Figure BDA0003905419020000052
Calculating T i Angle of tangency (slope of the curve) alpha in time window deltat time on the time curve i ,α i Is a dimensionless natural number.
Figure BDA0003905419020000061
In the formula (3), Δ t is the tangent angle time window, α i Is t i Angle of tangency at time, S i 、t i Respectively corresponding accumulated deformation and time S of the ith monitoring data i-1 、t i-1 Respectively the accumulated deformation amount and the time corresponding to the i-1 th monitoring data.
Now, Δ t is set to a plurality of fixed values, and the tangent angle calculation method is changed from equation (3) to equation (4).
Figure BDA0003905419020000062
In equation (4), Δ t is a fixed value, and the value of n depends on Δ t.
Then m fixed time windows are selected:
Δt j =j,j=1,2,3,…,m (5)
in the formula (5), Δ t j The unit is h, and m is set according to the deformation frequency of the monitored object.
Calculating a time window Δ t j Corresponding tangential angle alpha i,j And obtaining m tangent angle time curves.
Figure BDA0003905419020000063
The fixed time windows are sequenced from small to large, and each time t is calculated i Lower, tangent angle difference C corresponding to adjacent time window i,j
C i,j =|α i,ji,j+1 | (7)
C i,j When t is less than or equal to 3, t is considered to be i Time alpha i,j Curve and alpha i,j+1 Curve information is superposed, and the ratio of the information superposition time length to the total time length is alpha i,j Curve and alpha i,j+1 Coincidence ratio C of curves j (equation (8)). The curves with a coincidence rate of 80% have high similarity, so when C is j When the time is more than or equal to 80 percent of the first time, the model adaptively determines j hours as the optimal time window.
Figure BDA0003905419020000064
2. Preliminary determination of tangent angle
The optimum time window is substituted into equation (3) to calculate the tangent angle α value at each time. If the tangent angle alpha value is larger than or equal to the threshold theta, time forecasting calculation is needed, otherwise, time forecasting is not needed. The value of the threshold θ depends on the type of landslide, typically θ ∈ [60 °,85 ° ].
3. Pearson correlation coefficient quadratic decision
The inverse speed of the landslide is in linear relation with time, and the linear correlation between the inverse speed and the time is measured by calculating the Pearson correlation coefficient of the inverse speed, so that whether the landslide is in a temporary slip state or not can be judged.
Inverse speed IV i The same time window as the tangent angle is used for the calculation of (c).
Figure BDA0003905419020000071
In the above formula, Δ t is the optimal time window obtained by the previous calculation; s i 、S i-n Respectively, t in the deformation time curve i Time t and i-n and (4) accumulated deformation corresponding to the moment.
Then t i Pearson's correlation coefficient r at time i Comprises the following steps:
Figure BDA0003905419020000072
wherein IV i-n ~IV i Are each t i-n ~t i The inverse of the velocity at that moment.
If the Pearson correlation coefficient is less than or equal to the threshold value beta, time forecasting calculation is needed, otherwise, time forecasting is not needed. The value of the threshold value beta depends on the type of landslide, usually beta e-0.90, -0.70.
4. Output time forecast
When the tangent angle and the pearson correlation coefficient satisfy the above conditions at the same time, linear fitting is performed on the inverse velocity scatter (time on the abscissa and inverse velocity value on the ordinate) in the time window, and the intersection of the fitted straight line and the abscissa is the predicted disaster occurrence time. And finally, outputting the landslide forecast time.
The deformation time curve of the invention can be obtained by adopting ground GNSS equipment or ground-based interferometric radar (GB-InSAR).
Particularly, for a mode of monitoring by using GB-InSAR, considering the technical advantage that the GB-InSAR can realize surface scene deformation monitoring, the tangent angle calculation method based on the area-time curve is provided based on the characteristic that a deformation area can show regular expansion in the landslide development process. And by setting a plurality of deformation threshold values, fusing tangent angle calculation results of a plurality of deformation time curves to obtain an optimal tangent angle for judging the landslide early warning stage. Compared with the traditional early warning method based on single-point deformation information, the method of the area time sequence curve can identify the whole change trend of the deformation of the landslide surface, the early warning information is more comprehensive and accurate, and the accuracy of judgment in the landslide early warning stage can be effectively improved. The specific method flow is shown in fig. 6, and specifically includes the following steps:
s101, designing a deformation threshold value
And monitoring the scene by adopting GB-InSAR to obtain the accumulated deformation of the scene under the current time window delta t. According to the maximum value S of the accumulated deformation quantity of the scene under the current time window delta t max Designing N deformation thresholds:
Figure BDA0003905419020000081
the number and the numerical value of the deformation threshold can be determined by experience or numerical simulation. The design of this embodiment is as follows:
when S is max When the thickness is less than or equal to 100mm,
Figure BDA0003905419020000082
deformation threshold S k =10·k,k=1~N,
Figure BDA0003905419020000083
Meaning rounding down.
When S is max When the thickness is larger than 100mm,
Figure BDA0003905419020000084
deformation threshold
Figure BDA0003905419020000085
S102, acquiring an area-time curve
The deformation time curve based on single-point monitoring can better identify the landslide danger state of the point in the side slope deformation process. However, after a significant slope failure, the deformation rate of the single point may be reduced, but the deformation region may continue to expand. The method of the area-time curve can be used for more accurately identifying the whole deformation trend of the sliding slope surface.
For each deformation threshold S for the accumulated amount of deformation over a time window Δ t period k Calculating each time t in the time window delta t i Lower accumulated deformation amount exceeding S k Total area A of pixel points k,i All A within the time window Δ t k,i The set of (A) is the area-time curve A k (t) of (d). During the development of landslide, A k (t) is an increasing curve, A k The acceleration rate of (t) is the expansion rate of the strain region.
S103, calculating the tangent angle of the area-time curve
For each deformation threshold S k Corresponding area time curve A k (t) calculating an area-time curve A k (t) initial time t 1 To each time t i Average speed B of k,i (equation (11)).
Figure BDA0003905419020000091
Dividing the area by the average speed to obtain the ordinate A of the area-time curve k,i Transformed into T in accordance with the dimension of the abscissa time T k,i (equation (12)).
Figure BDA0003905419020000092
Calculating t i Tangential angle (slope) alpha at time k,i (equation (13)), α k,i Is a dimensionless natural number.
Figure BDA0003905419020000093
S104, obtaining the tangent angle of the early warning judgment
t i At the moment, each deformation threshold S k Corresponding to a tangent angle alpha k,i N tangent angles can be obtained by N deformation thresholds, and the maximum tangent angle alpha is obtained i As the tangent angle of the time window (equation (14)).
Figure BDA0003905419020000094
Case(s)
Taking the earth surface deformation data of the mountain and west sharp mountain iron ore side slope monitored by the foundation interference radar as an example, the forecasting and calculating process of the landslide critical-sliding time under the self-adaptive time window and multiple criteria is as follows:
1) Adaptive time window calculation: the multiple time windows, Δ t =1h,2h,3h,4h,5h,6h, \8230, 10h, are set, and the tangent angle α curves under the multiple time windows are shown in fig. 2. The coincidence rate of the tangent angle alpha curves under multiple time windows is shown in figure 3, and the coincidence rate of the tangent angle alpha curves under 4h and 5h time windows reaches 80 percent first, so that the self-adaptive time window is 4h.
2) Criterion 1 tangent angle: substituting the 4h time window into equation (3) calculates the tangent angle value at each time, see fig. 4. The threshold value theta in criterion 1 is set to 65 deg., the red part in fig. 4 corresponding to the moment of the tangent angle > 65 deg..
3) Calculating the speed reciprocal: substituting the 4h time window into equation (9) calculates the inverse velocity value at each time, see fig. 5. The inverse velocities in fig. 5 are filtered to remove maxima and negative values. In this example, values greater than 2000h/m and negative values are rejected to obtain the inverse effective velocity (FIG. 6).
4) Criterion 2 Pearson correlation coefficient: substituting the 4h time window and the filtered inverse velocity into a formula (10) to obtain the Pearson correlation coefficient r at each moment i See fig. 7. The threshold value beta in criterion 2 is set to-0.75 and the red part in fig. 7 corresponds to the moment when the correlation coefficient is < -0.75.
5) And (3) time forecast calculation: with t i The time instant of 51.33h is taken as an example, the tangent angle and the correlation coefficient corresponding to the time instant are 67.79 degrees and 0.98 degrees respectively, and the two criteria of the critical sliding stage are met, so the time pre-prediction needs to be carried outReporting and calculating. Linear fitting is carried out on the reciprocal velocity under the 4h time window, and the intersection point of the fitted straight line and the horizontal axis is t i =53.41h (fig. 8). Namely: t is t i The time of =51.33h, t is calculated i When the deformation is large in the case of =53.41h, the prediction is sent 2.98h in advance.
6) Time forecasting under double criteria:
and when the tangent angle is more than or equal to 65 degrees and the correlation coefficient is less than or equal to-0.75, performing time prediction calculation and outputting a prediction result, otherwise, no prediction is performed. Fig. 9 shows the forecast at each time point during monitoring, where the horizontal axis is 51.33h, and the vertical axis is 2.98h.
In summary, the above description is only a preferred embodiment of the present invention, and is not intended to limit the scope of the present invention. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims (10)

1. A landslide critical-slip time forecasting method based on an adaptive time window and multiple criteria is characterized by comprising the following steps:
s1, monitoring a scene to obtain a deformation time curve; self-adaptively obtaining an optimal time window; the method for acquiring the optimal time window specifically comprises the following steps:
s11, continuously monitoring a scene in real time to obtain a deformation time curve; setting a plurality of fixed time windows, respectively calculating tangent angles of the deformation time curves under the fixed time windows, and moving the fixed time windows to obtain tangent angle curves corresponding to the fixed time windows;
s12, sequencing the fixed time windows from small to large, and sequentially calculating the coincidence rate of tangent angle curves under the adjacent fixed time windows; when the coincidence rate exceeds a set threshold value B for the first time, the corresponding fixed time window is the optimal time window;
s2, calculating tangent angles of all moments based on the optimal time window obtained in the step S1;
s3, calculating the reciprocal of the speed at each moment based on the optimal time window obtained in the step S1, and further obtaining the Pearson correlation coefficient at each moment;
and S4, extracting the moments when the tangent angle is larger than the threshold theta and the Pearson correlation coefficient is smaller than the threshold beta, and outputting the landslide critical slip time predicted by adopting an inverse speed method at the moments.
2. The adaptive time window and multi-criterion based landslide critical-slip time forecasting method according to claim 1, characterized in that the tangent angle at a fixed time window Δ t is calculated as follows:
first, the accumulated deformation S is calculated i Average rate of (B) i
Figure FDA0003905419010000011
Wherein S is 1 、S i Respectively, t in the deformation time curve 1 Time t and i accumulated deformation corresponding to the moment;
then, using the accumulated deformation S i Divided by the average rate B i Obtaining T consistent with the dimension of time T i
Figure FDA0003905419010000012
Then t i Time of day, deformation time curve at fixed time window Δ t (Δ t = t) i -t i-n ) Lower tangent angle alpha i Comprises the following steps:
Figure FDA0003905419010000021
3. the adaptive time window and multi-criterion based landslide critical-slip time forecasting method according to claim 1, wherein in step S12, the coincidence rate of tangent angle curves in adjacent fixed time windows is calculated as follows:
calculating tangent angle difference values on tangent angle curves corresponding to adjacent time windows at the same moment, and when the difference value is less than or equal to a set threshold value A, considering that the information of the two tangent angle curves at the moment is overlapped; and the ratio of the information coincidence duration of the two tangent angle curves to the total duration is the coincidence rate of the two curves.
4. The adaptive time window and multi-criterion based landslide critical-slip time forecasting method according to claim 1, 2 or 3, wherein in S3, the Pearson' S correlation coefficient is calculated as follows:
calculating the reciprocal velocity IV in a fixed time window delta t i
Figure FDA0003905419010000022
In the formula (7), Δ t is the optimal time window obtained by S1; s i 、S i-n Respectively, t in the deformation time curve i Time and t i-n Accumulated deformation corresponding to the moment;
then t i Pearson's correlation coefficient r at time i Comprises the following steps:
Figure FDA0003905419010000023
wherein IV i-n ~IV i Are each t i-n ~t i The inverse of the velocity at that moment.
5. The adaptive time window and multi-criterion based landslide critical-slip time forecasting method according to claim 1, wherein in S3, the inverse speed at each time is filtered to eliminate maximum values and negative values.
6. The adaptive time window and multi-criterion based landslide time-to-slide forecasting method according to claim 1, characterized in that the threshold value θ depends on the type of landslide, θ e [60 °,85 ° ].
7. The adaptive time window and multi-criterion based landslide critical-slip time forecasting method according to claim 1, wherein the threshold β depends on the type of landslide, β e [ -0.90, -0.70].
8. The adaptive time window and multi-criterion based landslide critical-slip time forecasting method according to claim 1, wherein in S1, a deformation time curve is obtained by using a ground GNSS device or a ground-based interferometric radar.
9. The method for forecasting the landslide impending landslide time based on adaptive time window and multi-criteria as claimed in claim 8, wherein in S1, ground based interferometric radar is used to obtain the surface accumulated deformation amount of the monitored scene; determining a tangent angle based on the surface accumulated deformation, specifically comprising:
s101, designing a plurality of deformation thresholds according to the maximum value in the surface accumulated deformation quantity;
s102, acquiring an area-time curve corresponding to each deformation threshold:
aiming at each deformation threshold, calculating the total area of the pixel points of which the accumulated deformation exceeds the deformation threshold at the current moment in real time; each moment corresponds to a pixel point total area value exceeding the deformation threshold, and an area time curve is formed along with the time;
and S103, calculating tangent angles of the area-time curves corresponding to the deformation thresholds at the current moment aiming at the fixed time windows based on the self-adaptive time window method, wherein the largest tangent angle of the area-time curves corresponding to the deformation thresholds is the tangent angle at the current moment corresponding to the fixed time window.
10. The adaptive time window and multi-criterion-based landslide critical-slip time forecasting method according to claim 9, wherein in S101, N deformation thresholds are designed as follows:
maximum accumulated deformation value S in monitoring scene of current measurement period max When the thickness is less than or equal to 100mm,
Figure FDA0003905419010000031
deformation threshold S k =10·k,k=1~N,
Figure FDA0003905419010000032
Represents rounding down;
maximum accumulated deformation value S in monitoring scene of current measurement period max When the thickness is more than 100mm,
Figure FDA0003905419010000033
deformation threshold
Figure FDA0003905419010000034
CN202211304893.XA 2022-10-24 2022-10-24 Landslide and critical-sliding-time forecasting method based on self-adaptive time window and multiple criteria Pending CN115563802A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211304893.XA CN115563802A (en) 2022-10-24 2022-10-24 Landslide and critical-sliding-time forecasting method based on self-adaptive time window and multiple criteria

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211304893.XA CN115563802A (en) 2022-10-24 2022-10-24 Landslide and critical-sliding-time forecasting method based on self-adaptive time window and multiple criteria

Publications (1)

Publication Number Publication Date
CN115563802A true CN115563802A (en) 2023-01-03

Family

ID=84746057

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211304893.XA Pending CN115563802A (en) 2022-10-24 2022-10-24 Landslide and critical-sliding-time forecasting method based on self-adaptive time window and multiple criteria

Country Status (1)

Country Link
CN (1) CN115563802A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117192582A (en) * 2023-11-08 2023-12-08 航天宏图信息技术股份有限公司 Improved tangent angle real-time monitoring method, device and equipment for GNSS displacement data

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117192582A (en) * 2023-11-08 2023-12-08 航天宏图信息技术股份有限公司 Improved tangent angle real-time monitoring method, device and equipment for GNSS displacement data
CN117192582B (en) * 2023-11-08 2024-01-30 航天宏图信息技术股份有限公司 Improved tangent angle real-time monitoring method, device and equipment for GNSS displacement data

Similar Documents

Publication Publication Date Title
CN107255803B (en) Ship track logic initiation algorithm
CN102789545B (en) Based on the Forecasting Methodology of the turbine engine residual life of degradation model coupling
CN112348237B (en) Abnormal trend detection method for dynamic drilling data
CN109396375B (en) A kind of crystallizer bleedout prediction electric thermo method based on feature vector and hierarchical clustering
CN104392610A (en) Distributed video based highway traffic incident influence range dynamic monitoring method
CN111582188B (en) Concrete pouring side station supervision method based on artificial intelligence
CN107424441B (en) Based on Hotelling&#39;s T2Aircraft track inflection point detection and estimation method of statistics
CN115563802A (en) Landslide and critical-sliding-time forecasting method based on self-adaptive time window and multiple criteria
CN113284143B (en) Die casting deburring precision detection system based on image data processing
CN107895014B (en) Time series bridge monitoring data analysis method based on MapReduce framework
CN114360214B (en) Extra-large scale glacier debris flow early warning method
CN110569890A (en) Hydrological data abnormal mode detection method based on similarity measurement
CN113159374B (en) Data-driven urban traffic flow rate mode identification and real-time prediction early warning method
CN110196962A (en) A kind of air speed abnormality recognition method based on Density Estimator
CN111582380A (en) Ship track density clustering method and device based on space-time characteristics
CN115050214B (en) AIS data-based ship collision risk prediction method
CN116308305B (en) Bridge health monitoring data management system
CN115685189A (en) Method for improving landslide early warning stage discrimination based on self-adaptive time window
CN104180801A (en) Method and system for predicting track points based on ADS-B system
CN115457300A (en) Ship abnormal behavior detection method based on distance measurement and isolation mechanism
CN107730882B (en) Road congestion prediction system and method based on artificial intelligence
Wang et al. Classification of ship trajectories by using naive Bayesian algorithm
CN116681722A (en) Traffic accident detection method based on isolated forest algorithm and target tracking
CN113687357B (en) Multi-radar cross-regional networking multi-target tracking method
CN114623950A (en) Self-adaptive alarm constant value setting method for optical fiber temperature measurement system

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