Estimation of hydraulic conductivity of landslides based on support vector machine method optimized with genetic algorithm
-
摘要: 求解库岸边坡岩土体的渗透系数是研究滑坡渗流场及多场演化的基础，一般通过原位试验和室内试验求得，但试验成本较高且试验位置具有一定的随机性。本文以三峡库区马家沟滑坡为例，提出一种利用地下水位动态观测资料反演滑坡岩土层渗透系数的方法。具体步骤为：（1）依据滑坡的勘察资料和水位观测数据，构建滑坡数值模型；（2）利用SPSS生成不同渗透系数正交试验组合，并将渗透系数代入数值模型中计算监测井的水位，得到不同渗透系数及其对应的模拟水位数据；（3）应用遗传算法优化的支持向量机构建坡体模拟水位与渗透系数的非线性映射关系，再通过代入实际动态监测水位值求得滑坡岩土层的渗透系数；（4）将求得的渗透系数代入数值模型，用计算的模拟水位与实际观测水位进行对比验证。研究结果表明：遗传算法优化的支持向量机具有良好的学习预测效果，能准确预测渗透系数与水位的关系。该反演方法具有高效、准确的优点，反演结果的精度满足实际应用需要。Abstract: Estimation of hydraulic conductivities (K) of the rock media in a landslide is the basis for the study of the seepage field and multi-dimensional evolution of the reservoir bank slope. Traditionally, in-situ tests and indoor tests are used to determine the hydraulic conductivity of landslide rock and soil, but this method is costly and the test location has a certain randomness. In this study, the Majiagou landslide in the Three Gorges Reservoir area is taken as an example, and a method for inverting the K values of the deformed rock and soil mass using the groundwater level dynamic monitoring data is proposed. The basic idea is as follows. First, build a numerical model of the landslide based on the landslide survey data and water level observation data. Afterwards, SPSS is used to generate different orthogonal test combinations of hydraulic conductivity, substitute the hydraulic conductivity into the numerical model to calculate the water levels of the monitoring wells, and obtain the data of hydraulic conductivity and corresponding simulated water levels. Finally, the support vector machine (SVM) optimized with the genetic algorithm (GA) is used to construct a nonlinear mapping relationship between slope water level and hydraulic conductivities (K). The results obtained are then replaced for the monitored water levels to obtain the hydraulic conductivities of the landslide rock and soil which is used to develop the finite element model. The model is then verified by comparing the simulated water levels with the observed water levels. The inversion of the Majiagou landslide hydraulic conductivity shows that the SVM optimized with GA yields a good agreement between the simulated and real data and has a very efficient and accurate search results. The inversion accuracy of K based on the GA-SVM method meets the needs of practical applications.
-
库水位的涨落改变坡体力学平衡和软化岩层，决定坡体内渗流场的变化，对库岸滑坡稳定性有着显著影响[1]。因此，研究库岸边坡渗流场变化对研究库岸边坡变形演化过程具有重要意义[2-4]。滑坡岩土体渗透系数的求解是研究滑坡渗流场及渗流场变化引起的多场响应的前提，常见的渗透系数获取方法有原位试验和室内试验。然而，大型滑坡岩性组成较为复杂，试验结果随取样位置不同而存在较大随机性，且难以对滑坡深部的岩土体进行原位试验，因此测试结果往往仅具有参考性，不能反映岩层整体的渗透性质，无法直接用于数值模拟应用。随着机器学习方法和优化算法的发展和应用，利用反演法求解渗透系数受到了许多学者的青睐。
构建变量之间的映射关系是参数反演的前提，支持向量机（support vector machine, SVM）模型已被证明在构建非线性函数映射关系上具有很高的精度[5-7]。Prasad等[8]提出了一种遗传算法和数值模型结合的全局优化方法，并应用于马西右岸运河（MRBC）项目区渗透系数和补给参数的估算。魏进兵等[9]建立了泄滩滑坡饱和－非饱和渗流有限元计算模型，结合滑坡体地下水位监测资料，采用BP（back propagation）神经网络反演泄滩滑坡各地层渗透系数。崔皓东等[10]利用三维渗流场有限元结合遗传算法，建立二滩高拱坝坝基渗流场反演模型，通过测点水位反演获得了坝基渗透系数。姜谙男等[11]利用支持向量机建立水头与渗透系数之间的映射关系，以识别误差目标函数为适应值，通过粒子群优化算法反馈搜索得到了渗透系数。倪沙沙等[12]针对糯扎渡高心墙堆石坝，通过支持向量机建立渗透系数与水头之间的映射关系，参照大坝测压点水头，利用粒子群优化算法对大坝渗透系数进行了反演。向家松等[13]对比分析了不同地下水监测频率和时长对滑坡渗透系数反演结果的影响。
此外，部分学者还将机器学习与优化算法结合的方法应用于反演岩土体力学参数和变形参数。李端有等[14]以清江杨家槽滑坡为例进行研究，提出了将均匀设计、遗传算法与BP神经网络结合起来应用于滑坡反分析的方法，依据实测位移对滑坡岩土体力学参数进行了反演。Su等[15]利用支持向量机建立了大坝的变形预测模型，通过粒子群算法实现大坝变形预测模型的参数优化。Sun等[16]在三峡库区三门洞滑坡力学参数试验和位移观测资料的基础上，通过BP神经网络构建滑坡岩土体力学参数与位移量之间的关系，在此基础上，利用遗传算法反演得到了滑坡岩土体的力学参数。
但现有研究大多是利用某一时刻的水位或变形来反演相关参数，而库岸边坡水位变化实际上是一个动态的过程，因此，本文利用动态监测水位建立了一个新的反演方法，其基本思想是：根据边坡岩性和水文地质试验，确定待求渗透系数的取值区间；利用有限元软件生成渗透系数与对应水位的学习样本，在此基础上，通过机器学习方法构建水位与渗透系数的映射关系；最后，代入实测动态水位求解相应的岩土介质渗透系数。本文利用该方法反演了马家沟滑坡岩土体的渗透系数，为马家沟滑坡渗流场研究提供参数支持。
1. 学习方法简介
1.1 支持向量机
支持向量机是机器学习算法的一种，1964年由Vapnik等[17]在统计学习理论的基础上提出，该算法基于结构风险最小化原理，对高度复杂的非线性问题具备良好的适应性，可以用于解决分类和回归问题。
对于线性回归问题，假设n维空间训练样本集{Xi，Yi}（Xi∈Rn，Yi∈R，i=1，2，···，k，k为样本个数），构造最优线性函数：y=w·x+b，使得所有训练样本在精度
$\varepsilon $ 内都能用该线性函数拟合，且使得所有样本点产生的总偏差最小，引入允许拟合误差松弛因子${\xi _i}$ 、$\xi _i^*$ ，目标函数如下[15]：$$\begin{array}{l} \left\{ \begin{aligned} & \min \left( {\frac{1}{2}{{\left\| w \right\|}^2} + C\sum\limits_{i = 1}^k {\left( {{\xi _i} + \xi _i^*} \right)} } \right) \\ & {y_i} - w \cdot {x_i} - b \leqslant \varepsilon + {\xi _i} \\ & w \cdot {x_i} + b - {y_i} \leqslant \varepsilon + \xi _i^* \\ & {\xi _i} \geqslant 0 \\ & \xi _i^* \geqslant 0 \\ & i = 1,2, \cdot \cdot \cdot ,k \\ \end{aligned} \right\} \end{array} $$ (1) 式中：
$C$ ——惩罚因子，控制对超出误差$\varepsilon $ 的样本的惩 罚程度；w ——超平面函数法向量；
b ——超平面函数截距。
通过构建拉格朗日函数，得到对偶形式：
$$ \begin{split} &\max Q\left( {\alpha ,{\alpha ^*}} \right) = - \frac{1}{2}\sum\limits_{i = 1}^k {\sum\limits_{j = 1}^k {\left( {{\alpha _i} - \alpha _i^*} \right)} } \left( {{\alpha _j} - \alpha _j^*} \right)\left( {{x_i} \cdot {x_j}} \right)\\ &\quad\quad\quad\quad\quad\quad\;\;{ - \varepsilon \sum\limits_{i = 1}^k {\left( {{\alpha _i} + \alpha _i^*} \right)} + \sum\limits_{i = 1}^k {\left( {{\alpha _i} + \alpha _i^*} \right)} {y_i}}\\ &{\sum\limits_{i = 1}^k {\left( {{\alpha _i} - \alpha _i^*} \right)} = 0(0 \leqslant {\alpha _i},\alpha _i^* \leqslant C,i = 1,2, \cdots ,k}) \end{split} $$ (2) 得到拟合样本集的估计函数公式：
$$f\left( x \right) = \sum\limits_{i = 1}^k {\left( {{\alpha _i} - \alpha _i^*} \right)} \left( {{x_i} \cdot {x_j}} \right) + b$$ (3) 对于非线性回归问题，通过映射的方法将原始样本数据变换到高维特征空间（H），在高维特征空间中设计线性回归支持向量机，再通过原空间的函数内积运算，得到在原空间的非线性回归的估计函数。
不同的核函数能构造输入空间不同类型的非线性决策面，常见的满足Mercer条件的核函数主要有线性核函数、多项式核函数、多层感知机核函数（Sigmoid）、径向基核函数（Radial Basis Function, RBF）。其中，径向基核函数能实现非线性映射，参数较少且具有较宽的收敛域，适合作为通用模式识别的核函数。因此，本文使用径向基核函数。
1.2 遗传算法
遗传算法（genetic algorithm，GA）是一种基于自然选择机理的全局优化求解方法，20世纪70年代由美国学者Holland[17]提出，它模拟自然界中生物自然选择和遗传过程中发生的遗传、杂交和变异现象。遗传算法是一种基于“适者生存”的高度并行、随机和自适应的优化算法，通过复制、交叉、变异将解编码表示的“染色体”群一代代不断进化，最终收敛到最适应的群体，从而求得问题的最优解。其优点是原理和操作简单、通用性强、不受限制条件的约束，且具有隐含并行性和全局解搜索能力，在组合优化问题中得到广泛应用。
本文中，利用遗传算法优化寻找支持向量机学习效果最好的参数，即利用遗传算法优化支持向量机模型参数（遗传算法终止代数设置为50，初始种群数量设置为20，以向量机预测的渗透系数和样本渗透系数值差值最小为适应度函数，计算寻找最优的核函数参数），使得支持向量机模型达到最好的数据样本训练效果。
2. 基于遗传算法-支持向量机（GA-SVM）的马家沟滑坡渗透系数反演
2.1 马家沟滑坡概况
三峡库区沿岸地质条件较复杂，滑坡、崩塌众多，对周边居民的生命财产安全有重大影响[18]。马家沟滑坡位于湖北省秭归县，在长江支流吒溪河的左岸，是三峡库区典型的库岸边坡之一（图1）。马家沟滑坡自2003年三峡库区初次蓄水至135 m后发生明显变形，坡体后缘出现大规模的拉裂缝；在2007年做了抗滑桩工程措施，但由于深层滑带的存在以及库水位的周期涨落作用，滑坡的变形并未停止[19-20]。中国地质大学工程学院胡新丽课题组对马家沟滑坡进行了充分的工程勘察和长期的野外综合观测，掌握了马家沟滑坡的地质情况和多年综合观测资料，具备良好的研究基础。马家沟滑坡东西向展布，整体呈缓坡型，平均坡度为15°；滑坡滑动方向为291°，基本垂直于长江支流吒溪河。滑坡纵向长度为540 m，前后缘宽度分别为150，210 m（图2）。
滑坡前缘滑面位置在最低水位145 m以下，后缘高程为280 m。根据工程地质勘察，确定滑坡边界范围，计算出滑坡面积约为9.68×104 m2，滑坡深层滑带深度约30 m，初步估计失稳体积约为2.52×106 m3。
由图3可见，马家沟滑坡岩土体由第四系堆积物和下部的基岩构成。第四系松散堆积层广泛分布于坡体浅中层，一般厚度为8～15 m，总体上由滑坡后缘向前缘厚度递增；堆积物根据成因类型可以分为冲洪积物、崩坡积物、残坡积物，主要岩性为粉质黏土和碎石土；滑坡下部基岩地层主要是侏罗系上统遂宁组J3s，岩性以灰白色长石石英砂岩、细砂岩为主，夹有紫红色砂岩和粉砂质泥岩，力学强度较低，且泥岩夹层遇水易软化，是三峡地区易滑地层之一；岩层层面产状为270°～290°∠25°～30°，节理较发育，以微风化为主，上部岩石风化较严重，根据破坏模式和滑体成分，该滑坡属于顺层岩质滑坡。图中，G01、G02、G03、G04、G05为地表位移监测点，OFS1、OFS2、OFS3、JC1、JC3、JC6、JC8为深部位移监测点，JC1、JC3、JC8布设了水位计和渗压计，同时可以监测地下水位。
2.2 马家沟滑坡模型的建立
选取马家沟滑坡主剖面的前半部分，利用有限元软件Geo-Studio建立马家沟滑坡二维模型（图4），模型水平长度为385 m，垂直高度为230 m。第四系堆积层主要由含块石粉质黏土和碎石土组成，由于黏土仅分布于坡面前缘小部分和坡体后缘潜水面以上部分，因此在数值模型中将该层岩性简化为碎石土。下部基岩依据风化程度，在数值模型中以软弱夹层为界将基岩地层分为两层，上层为风化较严重的砂岩夹粉砂质泥岩，下层为稳定基岩。因此，滑坡数值模型中岩土介质根据岩性从坡体表面向内分别概化为碎石土、砂岩夹粉砂质泥岩（强风化基岩）和稳定基岩三种不同的岩土材料[20]。
各层岩土材料的非饱和渗流特性主要与渗透系数函数和土水特征曲线有关，模型中土水特征曲线选择软件中内置相应岩性的样本函数，输入饱和体积含水量（约为岩土层的孔隙率，见表2）得到，具体如图5所示。随后选择Fredlund和Xing方法估计渗透系数函数，结合土水特征曲线及饱和渗透系数得到渗透系数函数。数值模型左边界中库水位以下设置为定水头边界，水头大小等于随时间变化的库水位，右侧边界设置为定水头边界，水头大小等于随时间变化的监测孔JC8的水位，其他边界设置为无流量边界。在整个滑坡模型上有3个水位监测孔，从河岸往滑坡后缘依次为JC1、JC3、JC8（图3）。模型的初始水位根据库水位和三个观测水位相连而成，初始水位选择2013年1月16日马家沟滑坡初次监测的水位。
表 2 马家沟滑坡岩土体渗透系数取值范围表Table 2. Range of K of rock and soil for the Majiagou landslide滑坡地层 岩土体岩性组成 孔隙率 渗透系数范围/（cm·s−1） 第四系松散堆积层 碎石土 0.4 1.0×10−2～1.0×10−1 砂岩夹粉砂质泥岩 含裂隙岩体 0.3 1.0×10−3～1.0×10−2 基岩 稳定基岩 0.2 1.0×10−4～1.0×10−3 本文中地下水渗流问题的数学模型如下：
$$ \left\{\begin{split} &\frac{\partial }{\partial x}\left({k}_{x}\frac{\partial H}{\partial x}\right)+\frac{\partial }{\partial y}\left({k}_{y}\frac{\partial H}{\partial y}\right)+Q={m}_{{\rm{w}}}{\gamma }_{{\rm{w}}}\frac{\partial H}{\partial t}\\ &\quad\quad(0<{{x}},{{y}}<\infty ;t>0)\\ &{H\left(x,y\right)|}_{t=0}={H}_{0}\left(x,y\right)(0<{{x}},{{y}}<\infty) \\ &{H\left(x,y,t\right)|}_{\Gamma 1}={H}_{1}\left(x,y,t\right)(({{x}},{{y}})\in \Gamma 1;t>0)\\ &{T\frac{\partial H}{\partial n}\Big|}_{\Gamma 2}={q}_{2}\left(x,y,t\right)(({{x}},{{y}})\in \Gamma 2;t>0)\end{split}\right.$$ (4) 式中：kx、ky——x、y方向的渗透系数，土体非饱和时 是体积含水量和基质吸力的函数/ （m·s−1）；
H——总水头/m；
Q——施加的边界流量/（m3·s−1）；
mw——土水特征曲线的斜率；
γw——水的容重/（N·m−3）；
t——时间/s；
$\Gamma 1$ 、$\Gamma 2$ ——第一类、二类边界；n——边界的外法向量方向；
T——导水系数/（m2·s−1）；
q——单宽流量/（m·s−1）。
根据马家沟滑坡的岩性剖面，待反演的岩土体岩性组成从地表向下依次为第四系松散堆积层（碎石土）、强风化砂岩夹粉砂质泥岩（含裂隙基岩）以及稳定基岩。图6为现场入渗试验布点情况，参考现场入渗试验结果（表1）中各试验段岩性对应渗透系数大小，结合文[20]中所取渗透系数范围和模型试算结果，给定三种岩土体的渗透系数数量级范围，见表2。
表 1 马家沟滑坡入渗试验结果Table 1. Infiltration test results for the Majiagou landslide试验编号 试验深度/m 试验段岩性 渗透系数 /（cm·s−1） ZK4 0.50～0.67 粉质黏土（含块石） 1.38×10−5 ZK10 4.3～4.5 块石土 6.4×10−2 ZK11 6.1～6.3 块石土 0.5 ZK1 6.4～6.6 块石土 1.5 ZK2 10.0～10.3 砂岩块石（强风化） 2.5×10−2 ZK6 10.4～10.6 泥岩（强风化） 6×10−3 依据滑坡不同岩土体渗透系数取值范围，设计碎石土K1（第四系堆积物）、含裂隙岩体K2（强风化砂岩夹粉砂质泥岩）、基岩K3（稳定基岩）3个因子，每个因子在取值范围内取5个值，利用SPSS设置正交试验25组（表3）。将构造的渗透系数组合代入有限元软件进行模拟计算，获取相应监测点JC1、JC3的水头值，得到25组渗透系数及其对应模拟水位值H1（JC1）、H2（JC3）的样本。
表 3 数值模型计算方案表Table 3. Calculation schemes with the numerical model样本编号 K1/（cm·s−1） K2/（cm·s−1） K3/（cm·s−1） 样本编号 K1/（cm·s−1） K2/（cm·s−1） K3/（cm·s−1） 1 1.00×10−1 7.50×10−3 1.00×10−4 14 5.00×10−2 7.50×10−3 5.00×10−4 2 2.50×10−2 1.00×10−3 1.00×10−3 15 1.00×10−1 1.00×10−2 7.50×10−4 3 7.50×10−2 1.00×10−2 1.00×10−3 16 5.00×10−2 5.00×10−3 1.00×10−3 4 1.00×10−2 1.00×10−3 1.00×10−4 17 1.00×10−2 1.00×10−2 5.00×10−4 5 2.50×10−2 7.50×10−3 7.50×10−4 18 1.00×10−1 2.50×10−3 1.00×10−3 6 7.50×10−2 2.50×10−3 1.00×10−4 19 1.00×10−2 5.00×10−3 2.50×10−4 7 1.00×10−1 5.00×10−3 5.00×10−4 20 2.50×10−2 5.00×10−3 1.00×10−4 8 5.00×10−2 1.00×10−2 1.00×10−4 21 1.00×10−2 2.50×10−3 7.50×10−4 9 1.00×10−1 1.00×10−3 2.50×10−4 22 5.00×10−2 1.00×10−3 7.50×10−4 10 2.50×10−2 2.50×10−3 5.00×10−4 23 2.50×10−2 1.00×10−2 2.50×10−4 11 7.50×10−2 7.50×10−3 2.50×10−4 24 7.50×10−2 5.00×10−3 7.50×10−4 12 5.00×10−2 2.50×10−3 2.50×10−4 25 7.50×10−2 1.00×10−3 5.00×10−4 13 1.00×10−2 7.50×10−3 1.00×10−3 2.3 基于GA-SVM方法的滑坡渗透系数反演
将表3中的渗透系数-水位数据作为机器学习的样本，应用台湾大学Lin Chih-Jen开发的Libsvm（A Library for Support Vector Machines）工具箱[21]和李洋改进的LIBSVM faruto Version版本[22]，通过Matlab调用支持向量机模型进行计算，模型的训练效果见图7。图7（a1）（b1）（c1）分别表示K1、K2、K3训练过程中适应性曲线的变化，可以看到进化代数在20代之后，适应度值已基本接近于0，说明寻优的速度快，效果好；图7（a2）（b2）（c2）分别表示K1、K2、K3的训练效果，图中均方差（MSE）表示各数据误差平方和的平均数，其值越小表示预测值越接近实际值；线性回归决定系数（R2）表示自变量和应变量的相关关系，越接近于1，表示相关关系越强，预测性越好。可以看到，GA-SVM支持向量机反演渗透系数的方法具有较好的学习预测效果。
将实际观测水位输入到训练后的支持向量机，得到滑坡岩土体的渗透系数（表4）。将基于支持向量机反演得到的滑坡岩土体渗透系数代入数值模型中，得到不同时刻的地下水水位。
表 4 滑坡岩土体渗透系数反演值Table 4. Inversion values of K of landslide rock and soil mass岩土体材料 K1/（cm·s−1） K2/（cm·s−1） K3/（cm·s−1） 渗透系数反演值 1.31×10−1 1.11×10−2 9.95×10−4 图8为马家沟滑坡地下水水位监测情况，监测周期为2013年1月16日至2017年8月7日，监测频率约为2月/次。本文为了减小模型计算时间和数据量，选取2013年1月16日到2014年12月10日的监测数据来与模拟水位进行验证对比。
![]()
图 8 马家沟滑坡地下水水位监测结果Figure 8. Monitoring results of groundwater levels of the Majiagou landslide图9为马家沟滑坡钻孔监测水位与模拟水位对比结果，可以看出，滑坡体监测孔的实测水位值与模拟水位值很接近，说明反演结果比较可靠，基于GA-SVM支持向量机的反演方法具有良好的效果，准确可行。最终得到参数反演的结果为：K1=1.31×10−1 cm/s，K2=1.11×10−2 cm/s，K3=9.95×10−4 cm/s。
![]()
图 9 监测孔地下水水位模拟值与实测值的对比Figure 9. Comparison of the measured values with the simulated values of groundwater levels3. 结论
（1）基于遗传算法优化的支持向量机（GA-SVM）方法具有收敛快、学习预测效果好的优点，构建的水位与渗透系数映射关系非常精准，且训练后的支持向量机具有很好的泛化预测能力，能够预测动态水位对应的岩土体渗透系数，代入实测水位，得到马家沟滑坡岩土介质的渗透系数：碎石土（第四系松散堆积物）K1=1.31×10−1cm/s；含裂隙岩体（强风化砂岩夹粉砂质泥岩）K2=1.11×10−2 cm/s；稳定基岩K3=9.95×10−4 cm/s。
（2）通过代入反演所得岩土介质的渗透系数到数值模型中，计算模拟水位，并将其与实测水位对比验证，结果表明：模拟值与实测值拟合较好，说明反演结果可信度高，方法可行性强；同时，该方法可以给滑坡岩土体力学参数、强度参数的反演提供参考。
-
![]()
图 8 马家沟滑坡地下水水位监测结果
Figure 8. Monitoring results of groundwater levels of the Majiagou landslide
![]()
图 9 监测孔地下水水位模拟值与实测值的对比
Figure 9. Comparison of the measured values with the simulated values of groundwater levels
表 2 马家沟滑坡岩土体渗透系数取值范围表
Table 2 Range of K of rock and soil for the Majiagou landslide
滑坡地层 岩土体岩性组成 孔隙率 渗透系数范围/（cm·s−1） 第四系松散堆积层 碎石土 0.4 1.0×10−2～1.0×10−1 砂岩夹粉砂质泥岩 含裂隙岩体 0.3 1.0×10−3～1.0×10−2 基岩 稳定基岩 0.2 1.0×10−4～1.0×10−3 
下载: 导出CSV
表 1 马家沟滑坡入渗试验结果
Table 1 Infiltration test results for the Majiagou landslide
试验编号 试验深度/m 试验段岩性 渗透系数 /（cm·s−1） ZK4 0.50～0.67 粉质黏土（含块石） 1.38×10−5 ZK10 4.3～4.5 块石土 6.4×10−2 ZK11 6.1～6.3 块石土 0.5 ZK1 6.4～6.6 块石土 1.5 ZK2 10.0～10.3 砂岩块石（强风化） 2.5×10−2 ZK6 10.4～10.6 泥岩（强风化） 6×10−3
下载: 导出CSV
表 3 数值模型计算方案表
Table 3 Calculation schemes with the numerical model
样本编号 K1/（cm·s−1） K2/（cm·s−1） K3/（cm·s−1） 样本编号 K1/（cm·s−1） K2/（cm·s−1） K3/（cm·s−1） 1 1.00×10−1 7.50×10−3 1.00×10−4 14 5.00×10−2 7.50×10−3 5.00×10−4 2 2.50×10−2 1.00×10−3 1.00×10−3 15 1.00×10−1 1.00×10−2 7.50×10−4 3 7.50×10−2 1.00×10−2 1.00×10−3 16 5.00×10−2 5.00×10−3 1.00×10−3 4 1.00×10−2 1.00×10−3 1.00×10−4 17 1.00×10−2 1.00×10−2 5.00×10−4 5 2.50×10−2 7.50×10−3 7.50×10−4 18 1.00×10−1 2.50×10−3 1.00×10−3 6 7.50×10−2 2.50×10−3 1.00×10−4 19 1.00×10−2 5.00×10−3 2.50×10−4 7 1.00×10−1 5.00×10−3 5.00×10−4 20 2.50×10−2 5.00×10−3 1.00×10−4 8 5.00×10−2 1.00×10−2 1.00×10−4 21 1.00×10−2 2.50×10−3 7.50×10−4 9 1.00×10−1 1.00×10−3 2.50×10−4 22 5.00×10−2 1.00×10−3 7.50×10−4 10 2.50×10−2 2.50×10−3 5.00×10−4 23 2.50×10−2 1.00×10−2 2.50×10−4 11 7.50×10−2 7.50×10−3 2.50×10−4 24 7.50×10−2 5.00×10−3 7.50×10−4 12 5.00×10−2 2.50×10−3 2.50×10−4 25 7.50×10−2 1.00×10−3 5.00×10−4 13 1.00×10−2 7.50×10−3 1.00×10−3
下载: 导出CSV
表 4 滑坡岩土体渗透系数反演值
Table 4 Inversion values of K of landslide rock and soil mass
岩土体材料 K1/（cm·s−1） K2/（cm·s−1） K3/（cm·s−1） 渗透系数反演值 1.31×10−1 1.11×10−2 9.95×10−4
下载: 导出CSV
-
[1] 向家松, 文宝萍, 陈明, 等. 结构复杂滑坡活动对库水位变化的响应特征—以三峡库区柴湾滑坡为例[J]. 水文地质工程地质,2017,44(4):71 − 77. [XIANG Jiasong, WEN Baoping, CHEN Ming, et al. Activity response of a landslide with complex structure to fluctuation of reservoir water level: a case study of the Chaiwan landslide in the Three Gorges Reservoir[J]. Hydrogeology & Engineering Geology,2017,44(4):71 − 77. (in Chinese with English abstract) [2] SUN G H, ZHENG H, TANG H M, et al. Huangtupo landslide stability under water level fluctuations of the Three Gorges reservoir[J]. Landslides,2016,13(5):1167 − 1179. DOI: 10.1007/s10346-015-0637-7
[3] HU X L, ZHANG M, SUN M J, et al. Deformation characteristics and failure mode of the Zhujiadian landslide in the Three Gorges Reservoir, China[J]. Bulletin of Engineering Geology and the Environment,2015,74(1):1 − 12. DOI: 10.1007/s10064-013-0552-x
[4] 周剑, 邓茂林, 李卓骏, 等. 三峡库区浮托减重型滑坡对库水升降的响应规律[J]. 水文地质工程地质,2019,46(5):136 − 143. [ZHOU Jian, DENG Maolin, LI Zhuojun, et al. Response patterns of buoyancy weight loss landslides under reservoir water level fluctuation in the Three Gorges Reservoir area[J]. Hydrogeology & Engineering Geology,2019,46(5):136 − 143. (in Chinese with English abstract) [5] 黄发明, 殷坤龙, 何涛, 等. 库岸滑坡地下水位时间序列混沌特征识别与PSO-LSSVM模型预测[J]. 地质科技情报,2015,34(6):186 − 192. [HUANG Faming, YIN Kunlong, HE Tao, et al. Chaotic characteristics identification and prediction using PSO-LSSVM model of reservoir landslide groundwater level time series[J]. Geological Science and Technology Information,2015,34(6):186 − 192. (in Chinese with English abstract) [6] SURYANARAYANA C, SUDHEER C, MAHAMMOOD V, et al. An integrated wavelet-support vector machine for groundwater level prediction in Visakhapatnam, India[J]. Neurocomputing,2014,145:324 − 335. DOI: 10.1016/j.neucom.2014.05.026
[7] HE Z B, WEN X H, LIU H, et al. A comparative study of artificial neural network, adaptive neuro fuzzy inference system and support vector machine for forecasting river flow in the semiarid mountain region[J]. Journal of Hydrology,2014,509:379 − 386. DOI: 10.1016/j.jhydrol.2013.11.054
[8] LAKSHMI PRASAD K, RASTOGI A K. Estimating net aquifer recharge and zonal hydraulic conductivity values for Mahi Right Bank Canal project area, India by genetic algorithm[J]. Journal of Hydrology,2001,243(3/4):149 − 161.
[9] 魏进兵, 邓建辉, 高春玉, 等. 三峡库区泄滩滑坡非饱和渗流分析及渗透系数反演[J]. 岩土力学,2008,29(8):2262 − 2266. [WEI Jinbing, DENG Jianhui, GAO Chunyu, et al. Unsaturated seepage analysis and back analysis of permeability coefficient for Xietan landslide in Three Gorges Reservoir area[J]. Rock and Soil Mechanics,2008,29(8):2262 − 2266. (in Chinese with English abstract) DOI: 10.3969/j.issn.1000-7598.2008.08.045 [10] 崔皓东, 朱岳明. 二滩高拱坝坝基渗流场的反演分析[J]. 岩土力学,2009,30(10):3194 − 3199. [CUI Haodong, ZHU Yueming. Back analysis of seepage field of Ertan high arch dam foundation[J]. Rock and Soil Mechanics,2009,30(10):3194 − 3199. (in Chinese with English abstract) DOI: 10.3969/j.issn.1000-7598.2009.10.052 [11] 姜谙男, 梁冰. 基于粒子群支持向量机的三维含水层渗流参数反馈识别[J]. 岩土力学,2009,30(5):1527 − 1531. [JIANG Annan, LIANG Bing. Feedback identifying seepage parameters of 3D aquifer based on particle swarm optimization and support vector machine[J]. Rock and Soil Mechanics,2009,30(5):1527 − 1531. (in Chinese with English abstract) DOI: 10.3969/j.issn.1000-7598.2009.05.059 [12] 倪沙沙, 迟世春. 基于粒子群支持向量机的高心墙堆石坝渗透系数反演[J]. 岩土工程学报,2017,39(4):727 − 734. [NI Shasha, CHI Shichun. Back analysis of permeability coefficient of high core rockfill dam based on particle swarm optimization and support vector machine[J]. Chinese Journal of Geotechnical Engineering,2017,39(4):727 − 734. (in Chinese with English abstract) DOI: 10.11779/CJGE201704019 [13] 向家松, 文宝萍, 高幼龙, 等. 地下水位监测频率和时长对滑体渗透系数反演结果的影响[J]. 水文地质工程地质,2018,45(5):86 − 92. [XIANG Jiasong, WEN Baoping, GAO Youlong, et al. Effects of frequency and interval of groundwater monitoring on the inversion coefficients of permeability of materials of a landslide[J]. Hydrogeology & Engineering Geology,2018,45(5):86 − 92. (in Chinese with English abstract) [14] 李端有, 甘孝清. 滑坡体力学参数反分析研究[J]. 长江科学院院报,2005,22(6):44 − 48. [LI Duanyou, GAN Xiaoqing. Mechanical parameter back analysis of landslide[J]. Journal of Yangtze River Scientific Research Institute,2005,22(6):44 − 48. (in Chinese with English abstract) DOI: 10.3969/j.issn.1001-5485.2005.06.014 [15] SU H Z, LI X, YANG B B, et al. Wavelet support vector machine-based prediction model of dam deformation[J]. Mechanical Systems and Signal Processing,2018,110:412 − 427. DOI: 10.1016/j.ymssp.2018.03.022
[16] SUN G H, ZHENG H, HUANG Y Y, et al. Parameter inversion and deformation mechanism of Sanmendong landslide in the Three Gorges Reservoir region under the combined effect of reservoir water level fluctuation and rainfall[J]. Engineering Geology,2016,205:133 − 145. DOI: 10.1016/j.enggeo.2015.10.014
[17] 陈海洋, 滕彦国, 王金生. 基于GA-SVR的渗透系数参数反演方法[J]. 水文地质工程地质,2011,38(2):14 − 18. [CHEN Haiyang, TENG Yanguo, WANG Jinsheng. Methods of estimation of hydraulic conductivity with genetic algorithm-support vector regression machine[J]. Hydrogeology & Engineering Geology,2011,38(2):14 − 18. (in Chinese with English abstract) DOI: 10.3969/j.issn.1000-3665.2011.02.003 [18] 刘广润, 徐开祥. 三峡水库沿岸移民区地质灾害防治研究[J]. 中国地质灾害与防治学报,2003,14(4):1 − 4. [LIU Guangrun, XU Kaixiang. Investigation on prevention and control of geologic hazards in the migration area along bank of the Three Gorges reservoir[J]. The Chinese Journal of Geological Hazard and Control,2003,14(4):1 − 4. (in Chinese with English abstract) DOI: 10.3969/j.issn.1003-8035.2003.04.001 [19] ZHANG Y M, HU X L, TANNANT D D, et al. Field monitoring and deformation characteristics of a landslide with piles in the Three Gorges Reservoir area[J]. Landslides,2018,15(3):581 − 592. DOI: 10.1007/s10346-018-0945-9
[20] 张玉明. 水库运行条件下马家沟滑坡—抗滑桩体系多场特征与演化机理研究[D]. 武汉: 中国地质大学（武汉）, 2018. ZHANG Yuming. Multi-field characteristics and evolution mechanism of Majiagou landslide-stablizing piles system under reservoir operations[D]. Wuhan: China University of Geosciences（Wuhan）, 2018. (in Chinese with English abstract)
[21] CHANG C C, LIN C J. Libsvm:a library for support vector machinea library for support vector machine[J]. ACM Transactions on Intelligent Systems and Technology,2011,2(3):1 − 27.
[22] 史峰, 王小川, 郁磊, 等. MATLAB神经网络30个案例分析[M]. 北京: 北京航空航天大学出版社, 2010: 102−135. SHI Feng, WANG Xiaochuan, YU Lei, et al. Analysis of 30 cases of MATLAB neural network[M]. Beijing: Beihang University Press, 2010: 102−135. (in Chinese)
-
期刊类型引用(6)
1. 冯谕，曾怀恩，涂鹏飞. 遗传算法下的滑坡蠕滑位移预测模型研究. 中国地质灾害与防治学报. 2024(01): 82-91 . 
百度学术
2. 李瑞晨，侯木舟，孔梦麟，谢昊含. 基于VMD-SegSigmoid-XGBoost-ClusterLSTM算法的山体滑坡表面位移预测. 科技通报. 2024(09): 111-115 . 百度学术
3. 牟家琦，庄建琦，王世宝，孔嘉旭，杜晨辉. 基于深度神经网络模型的雅安市滑坡易发性评价. 中国地质灾害与防治学报. 2023(03): 157-168 . 百度学术
4. 吴延浩，江思珉，吴自军. 地下水污染强度及渗透系数场的反演识别研究. 水文地质工程地质. 2023(04): 193-203 . 本站查看
5. 熊超，孙红月. 基于多因素-多尺度分析的阶跃型滑坡位移预测. 吉林大学学报(地球科学版). 2023(04): 1175-1184 . 百度学术
6. 杨康，薛喜成，李识博. 信息量融入GA优化SVM模型下的地质灾害易发性评价. 安全与环境工程. 2022(03): 109-118 . 百度学术
其他类型引用(3)
计量
- 文章访问数: 372
- HTML全文浏览量: 170
- PDF下载量: 321
- 被引次数: 9
邮件订阅
RSS