基于多滑面B型深孔测斜曲线的滑动面位置精准确定方法与流程

本发明属于滑坡监测，具体涉及一种基于多滑面b型深孔测斜曲线的滑动面位置精准确定方法。

背景技术：

1、在滑坡稳定性研究中，滑坡位移监测是重要的研究方向。通过滑坡位移监测，不仅可以得到坡体变形位移速率、方向等数据，且通过分析监测资料，可深化滑坡变形机制及变形特征的认识，为失稳斜坡整治提供重要依据。滑坡位移监测中，确定滑坡滑动面的准确位置是评价斜坡稳定性并对其进行有效整治的前提。

2、滑动面辨识方法分为非确定性方法和确定性方法。非确定性方法主要有简易力学判识法﹑物探法﹑数值模拟法等；非确定性方法辨识的滑面为推测性滑面，需确定性方法进行验证；此外，非确定性方法一般为辅助手段或进行初步研究的手段。而滑坡治理工程勘查、设计和施工阶段则采用确定性方法进行滑动面识别，确定性方法主要有野外地质判识法、勘探识别法和深层位移观测法。

3、目前，深层位移观测法是主要的滑动面判识方法，该方法根据钻孔在监测期内土体的深部累计位移曲线变形特征直观判断滑动面所在位置。常见的滑坡深孔测斜曲线（即深部累积位移-时间曲线）主要有“v”型、“b”型、“r”型、“钟摆”型及“复合”型等。

4、在“b”型深孔测斜曲线的滑动面位置识别上，传统的直观判断法存在一些明显的缺陷：1) 缺乏科学可靠的辨识依据，并无一套完整的数学计算理论方法将滑动面位置的确定可计算化、精确化；2) 受测斜仪中测点布置间隔的影响该方法只能粗略地判断滑动面位置，并不能精准地获取滑动面位置信息；3)深孔测斜曲线整体呈鼓包状近似“b”型，存在多个滑动面且滑动区间范围大，各滑动区间内滑面所在位置为鼓包凸起最大值处，但通常难以直接准确获取最大凸起处对应的深度；4)深部累计位移-时间曲线图的横纵坐标观测尺度也会极大地影响滑动面位置的识别，观测尺度过小曲线中的无关变形特征被放大，不利于滑面特征捕捉；相反，观测尺度过大，累计位移曲线变化特征不明显，难以确定滑动面的位置。

5、目前，在现有的深层位移观测法中，相关的研究主要集中在理论分析与实际变形曲线的对比分析上，鲜有测斜曲线监测数据与滑面位置辨识相结合的深入研究，基于不同类型的深孔测斜曲线都面临着一个问题，就是缺乏科学可靠的数学处理手段和数据挖掘方法。

技术实现思路

1、本发明提供了一种基于多滑面b型深孔测斜曲线的滑动面位置精准确定方法，目的在于提供一套完整的数学计算理论方法将滑动面位置的确定可计算化、精确化。

2、为此，本发明采用如下技术方案：

3、基于多滑面b型深孔测斜曲线的滑动面位置精准确定方法，具体步骤如下：

4、1）导入深孔测斜仪传感器的位移监测数据，依据监测数据绘制各传感器对应的钻孔深部累计位移-时间曲线图（即深孔测斜曲线），判断曲线图的变形特征是否为“b”型（图1），若为“b”型则继续执行以下步骤。

5、2）取步骤1）的深孔测斜仪的监测数据，采用式①分别计算该深孔测斜仪不同深度处的土体在监测期内的位移速率；先取深孔测斜仪某一深度的位移数据计算位移速率，计算完成后再取下一深度的位移数据进行计算，依次计算得到传感器所在各个深度的位移速率。计算模型如图2所示。

6、vi=(si+1-si)/(ti+1-ti) ①

7、式①中：计算参数的对象均为同一深度处的土体，ti为深部位移采集时间；ti+1为与ti相邻的下一次深部位移采集时间；si为ti当日的累计位移；si+1为ti+1当日的累计位移；vi为ti至ti+1内的位移速率。

8、a.选择步骤2）的数据，绘制钻孔不同深度处的位移速率-时间曲线图，在位移速率-时间曲线图上确定位移速率曲线分簇特征显著的日期；以时间为横坐标速率值为纵坐标，则理想状态下土体位移速率沿深度方向的速率曲线分布特征如图3所示。

9、b.提取步骤a分簇特征显著日期内不同深度处的土体位移速率值，以深度为横坐标，土体位移速率为纵坐标，绘制位移速率-深度散点图，根据散点图的分布特征，初步确定滑坡体各个滑动区间的深度范围；如图4所示。

10、3）由式①计算得到不同深度处土体在监测期内的位移速率值，基于式①的位移速率计算结果，运用式②计算不同深度处土体在监测期内的土体加速度；先取深孔测斜仪某一深度的位移速率数据计算加速度，计算完成后再取下一深度的位移速率数据进行计算，依次计算得到传感器所在各个深度的加速度；

11、 ai=(vi+1-vi)/(ti+1-ti)      ②

12、式②中：计算参数的对象均为同一深度处的土体，ti为深部位移采集时间；ti+1为与ti相邻的下一次深部位移采集时间；vi为ti当日的位移速率；vi+1为ti+1当日的位移速率； ai为ti至ti+1内的位移加速度。

13、4）选择步骤3）的数据，绘制钻孔不同深度处的土体加速度-时间曲线图，掌握坡体的整体运动状态。其中，若将图3中的位移速率-时间曲线转变为加速度-时间曲线，则两者的曲线分簇特征相似，且由于不同位置处土体彼此间的位移速率比例明显小于加速度比例，使得曲线簇之间的间隔区间更大。因此，在加速度-时间曲线图上确定加速度曲线分簇特征显著的日期更加科学合理。

14、5）提取步骤4）分簇特征显著日期内不同深度处的土体加速度值，以深度为横坐标，土体加速度为纵坐标绘制加速度-深度散点图；根据散点图的分布特征，并在步骤b的基础上，进一步缩小滑动区间深度范围。本发明采用加速度-时间曲线图可划定更小的滑动区间深度范围，便于后续计算并精确确定滑动面的深度，降低计算误差。

15、6）提取任一滑动区间内离散点的加速度和土体深度数据，以深度为横坐标，土体加速度为纵坐标绘制加速度-深度散点图，如图5所示。从中可以发现，各滑面土体加速度分布区间内的散点分布更加密集，且各区间内的加速度散点簇具有更加显著的游离特征；不同的土体加速度分布区间彼此之间的散点簇间隔更大，即散点簇间隔区间1、2明显，相较于速率散点图能够更加直接准确的辨识坡体滑动区间及其对应的深度范围。

16、由步骤6）可以确定各个滑动区间深度范围，提取任一滑动区间内离散点的加速度和土体深度数据，选择合适的数学方法计算该区间内最大加速度值及其对应的土体深度。

17、为了能够穿越区间内的所有数据点，并且保持曲线连续光滑的特性，使其符合“b”型曲线的变形特征，本发明选择采用“三次样条插值法”。相比于其他数学方法，该方法具有更高的稳定性、贴合度、以及精度。

18、三次样条插值法的本质是在相邻的数据节点间构成一段三次多项式，再把每一段三次多项式按顺序拼接成光滑的曲线，具有更好的收敛性，能够有效克服牛顿插值法、拉格朗日插值法存在的尖点问题以及多项式阶数过高时存在的龙格现象，且不需要求解过多的导数信息、计算更加简便，适合用于加速度最大值点的计算。

19、因此，选择三次样条插值法更加符合“b”型曲线滑动区间内离散点的数学运算，其数学计算模型如下

20、运用三次样条插值法计算该区间内最大加速度值及其对应的土体深度。

21、其中，三次样条插值法计算原理如下：

22、假设一系列互不相等的数据节点[ x i, y i]( i=0,1...n)存在某种关系 y i =f ( x i)，其中 a=x 0 <x 1 <...<x n =b。根据三次样条插值理论，以区间[ x j, x j+1]为例， y j =f( x j) =s( x j)， s( x)的二阶导数 s ，，( x) =m  j ( j=0,1...n)，则该区间上 s( x)的三次样条插值表达式如公式③所示：

23、③

24、由此可以得到各滑动区间内的最大加速度及其对应的土体深度，由于相同时间内滑动面所在位置处的土体位移速率最快，故该位置处的土体加速度也是最大的。因此，可将加速度最大值对应的土体深度作为滑动面的深度，则由步骤8的计算结果可准确获取滑动面的深度位置。

25、本发明的有益效果在于：

26、1）具有更加科学可靠的理论依据。研究基于滑坡体概化模型结合“b”型曲线滑坡体变形特征，考虑不同位置处的土体变形规律，深入分析了模型变形与位移速率及加速度之间的潜在联系，并在大量的实际案例分析中得到验证。

27、2）充分挖掘“b”型深孔测斜曲线的关键特征，具有广泛的适用性。研究表明不同滑动区间的土体位移速率差异显著，且当滑动面的发展变形达到一定程度时，土体的位移速率沿深度方向会形成明显曲线簇。进一步地，将位移速率转化为加速度，则土体加速度曲线的分簇特征更加显著，能够明显区分土体的滑动区间，可为滑动面位置的准确辨识提供依据。

28、3）具有合理的数学依据，有着更高的可靠度。研究提出的方法通过提取监测期内曲线分簇特征显著的某日土体加速度数据，运用三次样条插值法计算各滑动面所在滑动区间内的加速度最大值，并以该值对应的土体深度作为滑面的计算深度，从而使滑动面位置的确定可计算化。

29、4）相比于传统的直观判别法和理论推导法，本发明提出的计算方法能够有效克服深孔测斜曲线传统辨识方法中存在的滑动面定义不清晰、数值不确定的问题，避免了因观测尺度选取不当导致曲线特征误判以及人为主观判别导致滑面辨识误差较大的情况，有助于降低滑坡施工治理成本。