[发明专利]一种基于临滑区变形信息的滑坡预报方法

摘要

本发明涉及露天矿山边坡、道路边坡及山体边坡的一种滑坡灾变时间预报预测方法。该方法根据现场监测数据和滑坡灾害的蠕变规律和机理，提出了临近滑坡但未滑落部位（临滑区）边坡蠕变曲线和变形发展四个阶段，推导了临滑区边坡变形、速率、加速度曲线方程及基于临滑区变形信息的滑坡预报时间模型，并建立了基于临滑区变形信息的Verhulst灰色模型以确定该滑坡预报时间模型中的参数。本发明的有益效果是：较现有其它滑坡预报方法工作量大大减少，同时提高了精度，简化了滑坡灾害时间预报的判断方法，更加方便和可靠。

主权项

一种基于临滑区变形信息的滑坡预报方法，其特征在于，具体包括以下步骤：(1)选定临滑区边坡，布设监测点进行数据采集，获得原始监测数据列；(2)由原始数据列得到累加生成数据列与紧邻均值生成序列，带入以下公式(1)：[m，n]T＝(BTB)‑1BTY    (1)其中，式中 B = - z ( 1 ) ( 2 ) ( z ( 1 ) ( 2 ) ) 2 - z ( 1 ) ( 3 ) ( z ( 1 ) ( 3 ) ) 2 . . . . . . - z ( 1 ) ( n ) ( z ( 1 ) ( n ) ) 2 , Y = x ( 0 ) ( 2 ) x ( 0 ) ( 3 ) . . . x ( 0 ) ( n ) ; x(0)(i)为原始数据列；z(1)(i)为x(1)(i)紧邻均值生成序列；x(1)(i)‑x(0)(i)为累加生成数据列，求解公式1获得系数m、n；(3)将上述步骤的系数m、n带入坡变形、速率、加速度曲线方程，如公式(2)所示： x = m x 1 n x 1 + ( m - n x 1 ) e m ( t - t 1 ) v = - m 2 x 1 ( m - n x 1 ) e m ( t - t 1 ) [ n x 1 + ( m - n x 1 ) e m ( t - t 1 ) ] 2 a = m 3 x 1 ( m - n x 1 ) e m ( t - t 1 ) [ ( m - n x 1 ) e m ( t - t 1 ) - n x 1 ] [ n x 1 + ( m - n x 1 ) e m ( t - t 1 ) ] 3 其中，x1为初始位移值t1为初始时间；得到边坡变形位移x、边坡变形速度v和为边坡变形加速度a；(4)根据步骤3得到的边坡变形位移x、边坡变形速度v和为边坡变形加速度a，再将步骤2得到系数m、n带入公式(3)： t r = 1 m ln ( n x 1 m - n x 1 ) + t 1 - - - ( 3 ) 得到滑坡发生的时间tr。