关键词:
滑坡
机器学习
变形预测
预测区间
稳定性分析
参数反演
摘要:
滑坡灾害预测预报和稳定性分析研究具有重大的工程意义和应用价值。机器学习作为目前炙手可热的人工智能方法,在各个行业均有所渗透,其也被逐步引入滑坡防治领域,大大促进了该领域的发展。然而,机器学习算法主要在滑坡位移确定性预测和区域滑坡易发性评价中应用较为广泛,在滑坡位移区间预测和滑坡稳定性分析中应用有限。因此,本文以机器学习为主要研究方法,借助机器学习算法在回归和聚类方面的优势,探索其在滑坡变形区间预测和稳定性评价上新的可能。首先根据滑坡变形的非线性特点,提出了三种基于机器学习算法的滑坡变形区间预测模型,随后又应用机器学习算法分别研究了库岸滑坡水力学参数全局敏感性分析、边坡临界滑面搜索、考虑滑面不确定性的滑坡抗剪强度参数概率反演等稳定性分析问题,主要的研究原因、过程和结果叙述如下:(1)准确可靠的滑坡位移预测滑坡的早期预警非常重要,机器学习方法因其强大的非线性处理能力而被广泛用于滑坡位移的确定性预测。然而,由于滑坡系统存在不确定性,传统的确定性预测方法不可避免地会产生预测误差。为了量化确定性预测结果的不确定性大小,提出了三种区间预测方法,即B-LSSVM、DS-LSSVM、DES-PSO-ELM,对滑坡位移进行区间预测,而不是确定性预测。BLSSVM方法计算过程是首先利用Bootstrap和LSSVM算法计算预测值的均值和模型方差,随后建立遗传算法优化的LSSVM模型对预测模型的偏差进行训练和预测,最后将预测模型的回归均值作为点预测值,位移的预测区间为预测均值、模型方差、噪声方差三者的和。DS-LSSVM方法是在原始LUBE方法框架下所提出的一种滑坡位移直接区间预测方法,该方法直接使用两个LSSVM模型来生成未来位移的上下边界,并通过DS算法最小化基于预测区间的适应度函数得到最佳模型参数。DES-PSO-ELM算法计算步骤为首先运用DES来预测累积滑坡位移数据的趋势项(线性分量),DES的预测残差表示非线性位移,然后采用集成的PSO-ELM方法构造非线性位移的预测区间,累积位移的最终预测区间通过将趋势位移预测值和非线性位移的预测区间相加得出。将这三种算法用于三峡库区白水河滑坡、树坪滑坡和谭家河滑坡的位移区间预测,并与已有算法进行对比,验证了所提方法的可靠性。(2)水力参数是库岸滑坡地下水位计算的关键数据,对水库滑坡稳定性评价至关重要,水力参数的敏感性分析可为库岸滑坡地下水位计算模型参数选取提供量化依据。引入了一种新的全局敏感性方法PAWN,并首次对库岸堆积层滑坡水力参数进行了全局敏感性分析。以三峡库区典型的库岸堆积层滑坡树坪滑坡为例进行研究,首先采用超拉丁立方采样法对饱和渗透系数、饱和体积含水量、剩余体积含水量和van Genuchten模型拟合参数进行采样,然后采用非侵入有限元法将其赋值到有限元模型中进行批量的渗流计算。根据有限元计算结果,采用PAWN方法计算了各水力参数的敏感度指标,划分了主要影响参数和非主要影响参数。全局敏感性结果对库岸滑坡现场调查和模型矫正有一定的参考价值,所展示的全局敏感性框架为类似滑坡的全局敏感性分析提供了参考。(3)基于强度折减的数值方法在边坡稳定性分析只能提供一条剪切带而无法提供一条清晰的临界滑动面。广泛使用的可视化技术来定义临界滑动面容易受到主观判断,并且该方法在批处理分析和三维分析中效率低下。当斜坡破坏时,滑面两侧的位移将明显不同,据此现象提出了一种自动识别临界滑动面的方法。首先采用k均值聚类算法将节点位移自动分为两类,分别代表滑体和稳定体,其次基于Alpha-Shapes算法来提取滑体的边界点,随后通过点云集合之间的布尔运算得到分界面附近的散点,最后采用分段三次样条拟合方法对提取的散点进行拟合,得到光滑的二维临界滑动面。对于三维边坡,采用薄板光滑样条拟合技术对滑面进行拟合。以均质边坡、含软弱夹层边坡和树坪滑坡为例,验证了所提方法的有效性和准确性。(4)先前关于滑坡抗剪强度参数的概率反分析的研究通常采用固定的滑动面,由于根据有限的观察结果推测出的滑动面位置存在不确定性,因此采用固定滑面进行参数概率反演可能会导致计算结果不可靠。基于贝叶斯定理,提出了考虑滑面不确定性的滑坡滑面抗剪强度参数概率反分析的框架,并用马尔可夫链蒙特卡罗模拟方法求解贝叶斯推理中抗剪强度参数的后验分布。为了提高计算效率,建立了基于极限学习机的响应面函数来逼近抗剪强度参数与相应的安全系数和临界滑动面位置之间的关系。通过一个已知实际抗剪强度参数和滑面的假想滑坡来验证所提方法以及对比所提出的方法和传统方法。研究了滑面测量误差和抗剪强度参数先验分布对概率反分析结果的影响。结果表明,忽略滑面不确定性的传统概率反分析所得的抗剪强度参数与实际值有较大偏差,且受抗剪强度参数先验均值的影响较大,所提方法反分析结果优于传统方法,且受抗剪强度参数先验