ISSN 1003-8035 CN 11-2852/P
  • 中国科技核心期刊
  • CSCD收录期刊
  • Caj-cd规范获奖期刊
  • Scopus 收录期刊
  • DOAJ 收录期刊
  • GeoRef收录期刊
欢迎扫码关注“i环境微平台”

滑坡多发区第四系堆积层厚度智能预测与影响因素分析以重庆万州铁峰乡为例

孟佳佳, 吴益平, 柯超, 苗发盛

孟佳佳,吴益平,柯超,等. 滑坡多发区第四系堆积层厚度智能预测与影响因素分析−以重庆万州铁峰乡为例[J]. 中国地质灾害与防治学报,2023,34(2): 1-10. DOI: 10.16031/j.cnki.issn.1003-8035.202202008
引用本文: 孟佳佳,吴益平,柯超,等. 滑坡多发区第四系堆积层厚度智能预测与影响因素分析−以重庆万州铁峰乡为例[J]. 中国地质灾害与防治学报,2023,34(2): 1-10. DOI: 10.16031/j.cnki.issn.1003-8035.202202008
MENG Jiajia,WU Yiping,KE Chao,et al. Intelligent prediction and analysis of influencing factors of Quaternary accumulation layer thickness in landslide-prone areas: A case study in the Tiefeng area of Wanzhou District, Chongqing City[J]. The Chinese Journal of Geological Hazard and Control,2023,34(2): 1-10. DOI: 10.16031/j.cnki.issn.1003-8035.202202008
Citation: MENG Jiajia,WU Yiping,KE Chao,et al. Intelligent prediction and analysis of influencing factors of Quaternary accumulation layer thickness in landslide-prone areas: A case study in the Tiefeng area of Wanzhou District, Chongqing City[J]. The Chinese Journal of Geological Hazard and Control,2023,34(2): 1-10. DOI: 10.16031/j.cnki.issn.1003-8035.202202008

滑坡多发区第四系堆积层厚度智能预测与影响因素分析——以重庆万州铁峰乡为例

基金项目: 国家自然科学基金项目(41977244;42007267)
详细信息
    作者简介:

    孟佳佳(1999-),女,河南开封人,硕士研究生,主要从事地质灾害预测预报研究。E-mail:jjm@cug.edu.cn

    通讯作者:

    吴益平(1971-),女,浙江杭州人,教授,博士生导师,从事地质灾害预测预报与风险评价研究。E-mail:ypwu@cug.edu.cn

  • 中图分类号: P642.22

Intelligent prediction and analysis of influencing factors of Quaternary accumulation layer thickness in landslide-prone areas: A case study in the Tiefeng area of Wanzhou District, Chongqing City

  • 摘要: 堆积层厚度是区域工程地质调查的基础数据,对识别滑坡具有重要意义,在滑坡的稳定性评价和风险评估中发挥着重要的作用。传统的空间插值方法对外界因素考虑不足,难以满足精度要求。文章以三峡库区万州区铁峰乡为研究对象,在实地调查的基础上,总结提出了堆积层的成因模式,并据此添加区域堆积层厚度控制点;基于Apriori算法挖掘影响因子与堆积层厚度分布之间的关联准则,利用已知样本点采用三种机器学习方法建立模型,最后将模型应用于整个区域得到堆积层厚度分布图,并对结果进行对比分析。结果表明:选取的13种影响因子中6种因子与堆积层厚度表现出较强的相关性;数据挖掘显示坡度和地形起伏是影响堆积层厚度空间分布差异的主要因素。3种机器学习模型中GWO-SVM模型的预测结果与实际最为吻合。研究结果揭示了区域第四系堆积层的成因机制,为机器学习技术在堆积层厚度预测领域奠定了基础。
    Abstract: Accumulation layer thickness is the basic data of regional engineering geological investigation, which is of great significance to identifying landslides and plays an important role in landslide stability evaluation and risk assessment. The traditional interpolation method does not take into account external factors, and it is difficult to meet the accuracy requirements. This study takes the Tiefeng Township of Wanzhou District in the Three Gorges Reservoir as the research object, summarizing the formation model of the accumulation layer and regional accumulation layer thickness control points are added accordingly. The Apriori algorithm is used to mine the association rules between the influence factors and the thickness distribution of accumulation layers. Three machine learning methods are used to build the model by using the known sample points, which are applied to the entire area to obtain the thickness distribution map of the accumulation layer, then compare the results. Results show that 6 of the 13 selected influence factors have a strong correlation with the accumulation layer thickness, and slope and topographic undulation are the main factors that cause the difference in the spatial distribution of accumulation layer thickness. Among the three machine learning methods, the prediction result of the GWO-SVM model is the most consistent with the actual situation. The results reveal the formation mechanism of regional quaternary accumulation layer and lay a foundation for machine learning in the field of accumulation layer thickness prediction.
  • 膨胀土边坡的稳定性一直是岩土界广泛关注的问题。目前,边坡稳定性分析的常用方法主要包括了极限平衡法、极限分析法等,都建立在极限平衡理论基础之上,并不适用于膨胀土边坡的稳定性分析[1]。另一种常用的方法是有限元强度折减法,早在1975年该方法就被Zienkiewice等[2]用来求解边坡稳定问题,随着计算机硬件技术和有限元软件技术的飞速发展,运用有限元强度折减法分析边坡稳定已经成为新的趋势[3-10]。国内很多学者将强度折减法运用到膨胀土边坡稳定分析中,取得了一系列成果。

    周健等[11]利用强度折减法研究膨胀土边坡的稳定性,发现干湿循环会导致膨胀土抗剪强度衰减,且随着干湿循环次数的增加,边坡稳定性降低,安全系数减小。刘明维等[12]研究了强度折减法在膨胀土斜坡地基路堤稳定性分析中的应用,发现强度折减法所得结果与实际情况相符。张硕等[3]基于有限元强度折减法研究了雨季土体增重、强度降低和膨胀作用对膨胀土边坡稳定性的影响,发现强度降低是导致边坡失稳的主要原因,膨胀作用次之,土体增重较小。程灿宇等[13]利用MIDAS/GTS、FLAC和ANSYS三种软件采用强度折减法分别对不同工况进行了稳定性分析,发现弱膨胀土边坡无论采用M-C屈服准则,还是D-P屈服准则所得结果差异不大。谭波等[14]采用强度折减法对不同条件下的膨胀土边坡的安全系数进行了计算,发现次生裂隙面发育是导致膨胀土边坡失稳的主要原因之一。杨才等[15]根据强度折减有限元法对不同条件失稳边坡稳定性分析结果,提出以最大塑性应变以及最小塑性应变的量级指标来判定塑性区贯通时刻。

    然而,干湿循环、降雨入渗等因素会引起浅层膨胀土干密度降低、吸力衰减,从而使抗剪强度大幅度下降。目前,在采用强度折减法分析膨胀土边坡稳定性的同时系统考虑抗剪强度衰减影响的研究尚不多见。为此,本文采用试验与数值模拟相结合的方式,系统地考虑了抗剪强度衰减特性的膨胀土边坡稳定性分析。首先对广西宁明膨胀土开展了室内直剪试验,分析了含水量、干密度对膨胀土抗剪强度衰减的影响;再以此为依据,利用Midas有限元分析软件研究考虑抗剪强度衰减特性对膨胀土边坡稳定性安全系数的影响,获取了边坡安全系数随抗剪强度折减的动态变化规律,以期为工程实践提供参考。

    土样取自广西崇左-夏石镇某高速公路膨胀土边坡路段,其天然含水量、最优含水量和天然干密度分别为32.5%,24%和1.40 g/cm3,其他土性指标,比重(Gs),液限(WL),塑限(WP),塑性指数(IP),自由膨胀率(σf)见表1。自由膨胀率为42.8%,按照《膨胀土地区建筑技术规范》[16]的分类,该膨胀土为弱膨胀性膨胀土。

    表  1  宁明膨胀土基本土体参数
    Table  1.  Basic soil parameters of Ningming expansive soil
    参数Gs/(g.cm−3wL/%wP/%IPσf/%
    取值2.8059.1124.6834.4342.8
    下载: 导出CSV 
    | 显示表格

    首先,将现场取回的扰动土试样碾散过2 mm筛,过筛后放入105℃的烘箱中烘24h,使试样具有相同的初始结构,并将烘干土用收纳箱密封保存备用。接着,按目标含水量(控制干密度为1.6 g/cm3)和目标干密度(控制含水量18%)要求配制成湿土,并装入保鲜袋,经闷料24 h后测得土样的最终含水量与目标含水量之间误差不超过1%;最后,为保证环刀试样均匀一致,采用自制的模具(图1)进行制样,并利用液压千斤顶脱模推出,控制试样的直径为61.8 mm,高度为15 mm,目的是使试样在竖直方向上能够充分膨胀,每组平行土样密度差不超过±0.02 g/cm3,否则废弃重做。试样配制过程如图2,最终制成的每个环刀试样表面均平整无破损,且长度误差不超过0.2 mm,则为满足要求的试样。

    图  1  制样模具
    Figure  1.  Sample preparation mould
    图  2  配土过程示意图
    Figure  2.  Diagram of the soil preparation process

    以初始干密度为1.6 g/cm3,含水量分别为9%、12%、15%、18%、21%、24%和27%制取环刀试样7组,每组4个;并以初始含水量为18%,干密度分别为1.4、1.5、1.6和1.7 g/cm3制取环刀试样4组,每组4个,然后进行常规直剪试验(图3),试验施加的竖向压力分别为100 kPa、200 kPa、300 kPa、400 kPa,剪切速率为0.02 mm/min,初始剪切位移均保持在3.850 mm左右,剪切位移量程13.000 mm。

    图  3  四联直剪仪
    Figure  3.  Quadruple direct shear testing device

    为研究广西宁明膨胀土的抗剪强度随含水量变化的规律,对不同含水量的土样进行直剪试验,试验结果如表2所示。

    表  2  宁明膨胀土抗剪强度试验结果表
    Table  2.  Results of shear strength of Ningming expensive soils
    试验参数w/%φ/(°)c/kPa
    试验结果8.8027.3100.36
    11.724.5693.28
    14.621.8067.34
    17.519.8254.64
    20.817.9241.22
    23.315.2030.86
    26.112.389.90
    下载: 导出CSV 
    | 显示表格

    根据表2可绘制出宁明膨胀土黏聚力和内摩擦角与含水量的关系如图4图5所示,拟合后可得到黏聚力和内摩擦角与含水量的关系式:

    图  4  宁明膨胀土黏聚力随含水量变化规律
    Figure  4.  Variation of cohesive force of Ningming expansive soil with water content
    图  5  宁明膨胀土内摩擦角随含水量变化规律
    Figure  5.  Variation of internal friction angle of Ningming expansive soil with water content
    $$ c = { - 5.192}w + 147.9 $$ (1)
    $$ \varphi = - 0.827w + 34.36 $$ (2)

    由式(1)和(2)可知,cφw都存在近似线性的关系,这与文献[17-18]结果一致,含水量每增大5%,其黏聚力约减小26 kPa,内摩擦角减小4.2°左右;为更好的表示cw的衰减规律,参考吕海波等[19]的研究,可计算出c的衰减率为:

    $$ \eta = \frac{{\left| {{c_0} - {c_1}} \right|}}{{{c_0}}} \times 100\% $$ (3)

    式中:η——黏聚力衰减率;

    c0——初始黏聚力;

    c1——随含水量变化后的黏聚力。

    根据表3可知,随着宁明膨胀土含水量的逐渐增大黏聚力不断衰减,在最低目标含水量9%以3%递增至目标含水量27%的过程中,黏聚力的衰减率变化趋势为增大-减小-增大,说明膨胀土在低含水量和接近饱和含水量时,黏聚力对含水量的变化显得十分敏感。

    表  3  宁明膨胀土黏聚力衰减率计算结果表
    Table  3.  Results of cohesion decay rate of Ningming expansive soil
    试验参数w/%c/kPaη/%
    试验结果8.8100.36
    11.793.287.05
    14.667.3427.81
    17.554.6418.86
    20.841.2224.56
    23.330.8625.13
    26.19.967.92
    下载: 导出CSV 
    | 显示表格

    在试样ρd保持一致的情况下(1.6 g/cm3),可从图6图7中看出在相同垂直应力作用下,抗剪强度随着w的增大呈现减小的趋势。

    图  6  不同含水量试样抗剪强度随垂直压力的变化
    Figure  6.  Change of the shear strength with vertical pressure of samples with different water contents
    图  7  不同荷载下试样抗剪强度随含水量的变化
    Figure  7.  Change of the shear strength with water content of specimens undergoing different vertical loads

    上述试验结果表明,宁明膨胀土的抗剪强度随着含水量的改变发生显著变化;主要表现为在含水量增大时黏聚力和内摩擦角发生衰减,其中黏聚力的衰减较内摩擦角更为明显。

    根据表4数据可拟合出试样黏聚力和内摩擦角随干密度的变化规律,如图8图9所示。

    表  4  不同干密度下试样试验结果记录表
    Table  4.  Record table of test results under different dry densities
    试验参数ρd/(g·cm−3c/(kPa)φ/(°)
    试验结果1.797.2626.5
    1.654.6419.82
    1.540.3417.82
    1.437.5716.87
    下载: 导出CSV 
    | 显示表格
    图  8  宁明膨胀土黏聚力随干密度变化规律
    Figure  8.  Variation of cohesive force of Ningming expansive soil with dry density
    图  9  宁明膨胀土内摩擦角随干密度变化规律
    Figure  9.  Variation of internal friction angle of Ningming expansive soil with dry density

    图8图9可观察出宁明膨胀土的黏聚力和内摩擦角随干密度的变化曲线符合乘幂函数的拟合结果,其中:

    $$ c = 0.126{{\rm{e}}^{3.884{\rho _{\rm{d}}}}} $$ (4)
    $$ \varphi = 1.631{{\rm{e}}^{1.614{\rho _{\rm{d}}}}} $$ (5)

    分析式(4)可知试样c随着ρd的减小而减小,且随着ρd的减小,c的衰减速率由快到慢,并最终趋于稳定;而在接近最大干密度(1.78 g/cm3)时变化较为显著,在干密度由1.4 g/cm3增大至1.6 g/cm3时,c增加了17.07 kPa;在干密度由1.6 g/cm3增大至1.7 g/cm3时,c增加了42.62 kPa。而由式(5)能看出φ亦随着ρd的减小而减小,但其整体的变化幅度并不大,干密度1.4 g/cm3与1.7 g/cm3的试样φ相差约9.6°;图10中各级载荷下的抗剪强度都随着试样ρd的减小而降低,且其变化幅度在高垂直应力条件下更为显著。

    图  10  不同干密度下试样抗剪强度随垂直应力的变化
    Figure  10.  Variation of shear strength with vertical stress of specimens of different dry densities

    干密度对宁明膨胀土抗剪强度的影响主要体现在黏聚力上,试样干密度越小,单位体积土体的土颗粒越少,土粒间水膜越薄,其抗剪强度越小;此外,膨胀土干密度越小,其吸力越大,试样的抗剪强度越低;而干密度对于内摩擦角的整体影响并不显著,其变化在10°以内。

    根据广西崇左-夏石镇某高速公路膨胀土边坡为研究对象,并参考该公路的地质勘察报告,该边坡土质主要由填土(①1和①2)、黏土②、强风化泥岩③和中风化泥岩④组成。同时根据地质调查及钻探、探槽揭示,该边坡滑动带基本位于黏土层,且下部强风化泥岩等土体不透水,大气影响深度为7 m,刚好大致为填土厚度和黏土厚度之和,影响急剧层深度为2.5 m。相关土层天然状态下基本参数指标见表5

    表  5  土层相关参数
    Table  5.  Soil layer related parameters
    地层岩性厚度
    /m
    重度
    /(kN·m−3
    内摩擦角
    /(°)
    黏聚力
    /kPa
    其它
    填土①10.2~118.0524成分黏土
    填土①22.5~3.318.8307上层砾砂,
    下层碎石
    黏土②0.3~418.48.435.6中等膨胀土
    强风化泥岩③0.6~119.32545质量等级Ⅴ级
    中风化泥岩④未钻穿19.63565质量等级Ⅴ级
    下载: 导出CSV 
    | 显示表格

    结合上述实际工程地质勘察报告,将膨胀土边坡考虑为非匀质边坡,同时为提高模型求解时间,取黏土弹性模量12000 kPa,容重18.4 N/m3,泊松比0.3,边坡高20 m,坡比1∶1.5。为避免尺寸效应带来的误差和便于模型求解收敛,坡顶取15 m,坡底取25 m,网格按线性梯度(长度)划分,起始长度1.2 m,结束长度0.5 m。由于填土土层由于土体较松散,易膨胀开裂,在降雨作用下容易引发降雨入渗,易软化下部土体,因此实际工程中对该部分填土进行了挖除。填土挖除后,为充分合理考虑到大气影响层对膨胀土边坡中黏土的影响,同时又不会影响到下部不透水泥岩,取大气影响层为距离坡面4 m范围的土体,正好为黏土厚度,急剧层为距离坡面1.5 m范围的土体(图11)。

    图  11  模型示意图
    Figure  11.  Numerical simulation model

    根据室内直剪试验结果,同时考虑到膨胀土具有浅层性,将测得的7个含水量下(干密度均为1.6 g/cm3)的膨胀土抗剪强度参数指标cφ赋予给受大气影响的风化层土体,即距离坡面4 m范围内的黏土。强、中风化泥岩层土体参数指标取地质勘察报告的值,具体数值见表5。计算得到不同含水量w下膨胀土边坡整体位移和潜在滑移面,如图12图13所示。

    图  12  1.6 g/cm3干密度不同含水量条件下的边坡位移
    Figure  12.  Slope displacement with the 1.6 g /cm3 dry density under different moisture content conditions
    图  13  1.6 g/cm3干密度不同含水量条件下的边坡潜在滑移面
    Figure  13.  Potential slip surface of slope with the dry density of 1.6 g/cm3under different moisture content

    分析图12图13可知,随着含水量w的增大,边坡的整体位移整体呈增大趋势,非饱和膨胀土边坡的浅层破坏由受大气影响层膨胀土强度衰减导致。随着含水量的增加,土体的c不断减小,边坡位移不断增大,滑移面逐渐变浅;破坏形式为浅层滑塌式的破坏。边坡失稳的滑移面位置位于大气影响层和不透水泥岩的交界处,且与黏土的底部相切。

    基于相同干密度,不同含水量下膨胀土的剪切试验和地质勘察报告,利用有限元分析软件对边坡进行稳定性分析,可得到随着膨胀土含水量的变化对边坡稳定性安全系数的影响规律,如图14所示的曲线,表达式为:

    图  14  边坡安全系数随含水量的变化规律
    Figure  14.  The variation of slope safety factor with water content
    $$ y = - {\text{0}}{\text{.008}}{x^2} + {\text{0}}{\text{.1884}}x + {\text{2}}{\text{.025}} $$ (6)

    随着w的增大,膨胀土的强度参数指标不断衰减,含水量较高比低含水量情况下的衰减速度更大。同时,膨胀土边坡在天然状况下处于稳定状态,但当w增大至27%时,其Fs为0.850,稳定性转变为失稳状态,发生滑坡、坍塌等工程现象;在此基础上,若继续增大含水量,膨胀土边坡将可能由浅层失稳进入完全失稳状态,这与实际工程中,在长时间降雨后,曾出现的多次滑坡现象类似。

    根据试验结果,将测得的四个干密度下(含水量均为18%)的膨胀土抗剪强度参数指标cφ赋予给距离坡面4 m范围的黏土。强、中风化泥岩层土体抗剪强度参数指标取地质勘察报告值,具体数值见表5。计算得到不同ρd下膨胀土边坡整体位移和潜在滑移面,如图15图16所示。

    图  15  18%含水量不同干密度条件下的边坡位移
    Figure  15.  Slope displacement under different dry densities with the moisture content of 18%
    图  16  18%含水量不同干密度条件下的边坡潜在滑移面
    Figure  16.  Potential slip surface of slope under different dry densities with the 18% moisture content

    图15图16中可以看出试样的ρd越小,边坡位移越大,潜在滑移面变浅;这是因为土体的c随着ρd的减小而减小,使得其抗剪强度降低;此时,边坡的破坏形式由整体滑动变为浅层滑塌。基于相同含水量,不同干密度下膨胀土的剪切试验和地质勘察报告,利用有限元分析软件对边坡进行稳定性分析,可得到随着膨胀土干密度的变化对边坡稳定性安全系数的影响规律,如图17所示的曲线,其表达式为:

    图  17  边坡安全系数随干密度的变化规律
    Figure  17.  The variation of slope safety factor with dry density
    $$ y = {\text{8}}{\text{.375}}{x^2} - {\text{23}}{\text{.24}}x + {\text{18}}{\text{.41}} $$ (7)

    试样ρd越小,其抗剪强度越低;且在ρd越大时其Fs增大趋势越为显著;1.5 g/cm3干密度下的Fs为2.409,比1.4 g/cm3的高出0.124,而1.7 g/cm3干密度下的Fs与1.6 g/cm3条件下的差值为0.459。

    (1)含水量的增大、干密度的减小都会引起膨胀土的峰值抗剪强度、黏聚力以及内摩擦角发生不同程度的衰减,其中,黏聚力的衰减幅度相较于内摩擦角更大。

    (2)通过多次膨胀土强度折减的方法可以很好地模拟降雨过程中由抗剪强度衰减引起的边坡稳定性的动态变化:风化层土体强度接近未风化层土体强度时,边坡处于稳定状态,潜在滑动面穿过分层界面;随着含水量增大、干密度变小,风化层抗剪强度会不断衰减,引起潜在滑动面逐渐外移,边坡稳定性降低。

    (3)数值模拟结果表明:与干密度减小相比,含水量的增大对边坡稳定更为不利,含水量增加到27%以后,膨胀土边坡由稳定状态变为欠稳定状态,因此在分析膨胀土边坡稳定性时,应着重考虑含水量变化的影响。

  • 图  1   Apriori算法示意图

    Figure  1.   Schematic diagram of Apriori algorithm

    图  2   万州区位置及铁峰乡滑坡分布图

    Figure  2.   Location of Wanzhou District and distribution of landslides in Tiefeng Township

    图  3   残坡积示意图

    Figure  3.   Schematic diagram of residual slope

    图  4   顺向坡远程滑动堆积示意图

    Figure  4.   Schematic diagram of long-distance sliding accumulation along the slope

    图  5   顺向坡近距离滑动堆积示意图

    Figure  5.   Schematic diagram of close range sliding accumulation along the slope

    图  6   顺向坡侧向崩塌堆积示意图

    Figure  6.   Schematic diagram of lateral collapse and accumulation along the slope

    图  7   软弱相间岩层分布模式示意图

    Figure  7.   Schematic diagram of the distribution mode of soft and weak alternating rock formation

    图  8   各因子统计图

    Figure  8.   Statistical chart of each factor

    图  9   各模型测试结果对比(测试样本)

    Figure  9.   Comparison of test results of each model (test samples)

    图  10   各模型的区域堆积层厚度预测结果

    Figure  10.   Prediction results of regional accumulation layer thickness of each model

    表  1   各因子与堆积层厚度相关性指标

    Table  1   Correlation index between each factor and accumulation layer thickness

    因子GRAMIC
    因子1坡度0.69230.6977
    因子2坡向0.68780.5285
    因子3流向0.86730.4696
    因子4高程0.70500.0752
    因子5汇水条件0.63490.0417
    因子6地形起伏0.68560.6933
    因子7地层0.75460.0486
    因子8距水系距离0.63660.0525
    因子9距陡崖距离0.60630.0637
    因子10剖面曲率0.69530.2141
    因子11平面曲率0.67830.0784
    因子12斜坡结构0.63640.3271
    因子13植被覆盖0.64400.0943
    下载: 导出CSV

    表  2   各因子属性指标

    Table  2   Attribute indicators of each factor

    因子属性类别
    因子1坡度0~20°F11
    20°~40°F12
    >40°F13
    因子2坡向337.5°~67.5°F21
    67.5°~157.5°F22
    157.5°~247.5°F23
    247.5°~337.5°F24
    因子3流向1&2F31
    4&8F32
    16&32F33
    64&128F34
    因子6地形起伏0~2F61
    2~4F62
    >4F63
    因子10剖面曲率0~20F101
    20~40F102
    40~60F103
    >60F104
    因子12斜坡结构顺向坡F121
    逆向坡F122
    斜交坡F123
    下载: 导出CSV

    表  3   堆积层厚度关联准则

    Table  3   Association criterion of accumulation layer thickness

    规则 ID规则支持度/% 置信度/%
    1YZ3=F31 & YZ6=F62 & YZ2=F21$\Rightarrow $HD50.08100.00
    2YZ3=F31 & YZ2=F21 & YZ1=F12$\Rightarrow $HD50.08100.00
    3YZ2=F22 & YZ3=F34 & YZ1=F11$\Rightarrow $HD30.08100.00
    4YZ2=F22 & YZ3=F34 & YZ6=F61$\Rightarrow $HD30.08100.00
    5YZ3=F34 & YZ1=F12 & YZ2=F24$\Rightarrow $HD20.17100.00
    6YZ2=F23 & YZ10=F102 & YZ1=F11$\Rightarrow $HD20.08100.00
    7YZ10=F102 & YZ2=F22 & YZ1=F11$\Rightarrow $HD20.08100.00
    8YZ3=F31 & YZ6=F63$\Rightarrow $HD118.58100.00
    9YZ3=F31 & YZ1=F13 & YZ6=F63$\Rightarrow $HD117.18100.00
    10YZ2=F22 & YZ1=F13 & YZ6=F63$\Rightarrow $HD117.01100.00
    下载: 导出CSV

    表  4   各模型预测精度对比

    Table  4   Comparison of prediction accuracy of each model

    预测模型MSERMSEQC平均厚度/m
    PSO-SVM401.221.420.470.791.49
    GA-SVM198.571.000.230.761.76
    GWO-SVM152.240.870.190.831.83
    下载: 导出CSV
  • [1] 王强,吴炳方,朱亮. 土壤厚度研究进展[J]. 安徽农业科学,2012,40(9):5273 − 5277. [WANG Qiang,WU Bingfang,ZHU Liang. Review of soil depth[J]. Journal of Anhui Agricultural Sciences,2012,40(9):5273 − 5277. (in Chinese with English abstract) DOI: 10.3969/j.issn.0517-6611.2012.09.071
    [2] 李岩,尚士友,阿拉塔其其格,等. 典型草原栗钙土层厚度空间插值方法研究[J]. 中国农机化学报,2015,36(1):344 − 348. [LI Yan,SHANG Shiyou,A Lataqiqige,et al. Research of spatial interpolation methods of typical grassland chestnut soil thickness[J]. Journal of Chinese Agricultural Mechanization,2015,36(1):344 − 348. (in Chinese with English abstract) DOI: 10.13733/j.jcam.issn.2095-5553.2015.01.090
    [3] 王桂林,向林川,孙帆. 粒子群优化协同克里金法在确定山地斜坡土层厚度中的应用[J]. 土木建筑与环境工程,2018,40(6):60 − 66. [WANG Guilin,XIANG Linchuan,SUN Fan. Application of cooperative Kriging method based on particle swarm optimization in estimation slope soil thickness[J]. Journal of Civil,Architectural & Environmental Engineering,2018,40(6):60 − 66. (in Chinese with English abstract)
    [4]

    SCARPONE C,SCHMIDT M G,BULMER C E,et al. Modelling soil thickness in the critical zone for Southern British Columbia[J]. Geoderma,2016,282:59 − 69. DOI: 10.1016/j.geoderma.2016.07.012

    [5]

    PRADHAN A M S,KIM Y T. Development and evaluation of relative relief based soil thickness model:A comparative study in hilly terrain,south Korea[J]. KSCE Journal of Civil Engineering,2021,25(6):2186 − 2198. DOI: 10.1007/s12205-021-1379-9

    [6]

    LI X C,LUO J H,JIN X L,et al. Improving soil thickness estimations based on multiple environmental variables with stacking ensemble methods[J]. Remote Sensing,2020,12(21):3609. DOI: 10.3390/rs12213609

    [7] 刘磊,殷坤龙,张俊. 三峡库区万州主城区第四系堆积层厚度的估算方法及应用[J]. 地质科技情报,2016,35(1):177 − 183. [LIU Lei,YIN Kunlong,ZHANG Jun. Estimation method of the quaternary deposits thickness and its application in Wanzhou central district,Three Gorges Reservoir region[J]. Geological Science and Technology Information,2016,35(1):177 − 183. (in Chinese with English abstract)
    [8] 段贵娟,李飒,孙立强,等. 基于水平相关距离的土层厚度区域评价方法研究[J]. 水力发电学报,2021,40(2):177 − 186. [DUAN Guijuan,LI Sa,SUN Liqiang,et al. Study of soil thickness regional evaluation using horizontal scales of fluctuation[J]. Journal of Hydroelectric Engineering,2021,40(2):177 − 186. (in Chinese with English abstract) DOI: 10.11660/slfdxb.20210218
    [9] 邵晨灿,向喜琼. 基于成因分类的土层厚度空间分布研究[J]. 中国水运(下半月),2021,21(3):115 − 117. [SHAO Chencan,XIANG Xiqiong. Study on spatial distribution of soil thickness based on genetic classification[J]. China Water Transport,2021,21(3):115 − 117. (in Chinese)
    [10]

    WANG Q,WU B F,STEIN A,et al. Soil depth spatial prediction by fuzzy soil-landscape model[J]. Journal of Soils and Sediments,2018,18(3):1041 − 1051. DOI: 10.1007/s11368-017-1779-0

    [11]

    AKUMU C E,WOODS M,JOHNSON J A,et al. GIS-fuzzy logic technique in modeling soil depth classes:Using parts of the Clay Belt and Hornepayne region in Ontario,Canada as a case study[J]. Geoderma,2016,283:78 − 87. DOI: 10.1016/j.geoderma.2016.07.028

    [12]

    CHENG W,ZHU A X,QIN C Z,et al. Updating conventional soil maps by mining soil-environment relationships from individual soil polygons[J]. Journal of Integrative Agriculture,2019,18(2):265 − 278. DOI: 10.1016/S2095-3119(18)61938-0

    [13]

    AGRAWAL R, IMIELIŃSKI T, SWAMI A. Mining association rules between sets of items in large databases[C]//Proceedings of the 1993 ACM SIGMOD international conference on Management of data. Washington, D. C. , USA. New York: ACM, 1993: 207-216.

    [14]

    VAPNIK V N. The nature of statistical learning theory[M]. New York: Springer , 2000.

  • 期刊类型引用(5)

    1. 贺伟明,石胜伟,蔡强,梁炯. 基于上下限解的膨胀土边坡首次滑动区域分析. 水文地质工程地质. 2025(01): 104-112 . 百度学术
    2. 孙银磊,余川,廖磊,李志妃. 钢渣粉固化改良膨胀性黏土机理研究进展. 水文地质工程地质. 2025(01): 113-129 . 百度学术
    3. 张锐,周豫,兰天,郑健龙,刘昭京,李彬. 高速铁路土工格栅加筋膨胀土边坡作用机制. 铁道科学与工程学报. 2024(01): 1-12 . 百度学术
    4. 纪佑军,熊军,蒋国斌,王泽根. 考虑应变软化的鸡场镇降雨型滑坡数值分析. 水文地质工程地质. 2024(04): 178-188 . 百度学术
    5. 张再江. 基于改进极限平衡原理的膨胀土边坡稳定性计算分析. 水利科技与经济. 2024(07): 48-51 . 百度学术

    其他类型引用(7)

图(10)  /  表(4)
计量
  • 文章访问数:  454
  • HTML全文浏览量:  493
  • PDF下载量:  180
  • 被引次数: 12
出版历程
  • 收稿日期:  2022-02-09
  • 修回日期:  2022-04-06
  • 网络出版日期:  2023-04-10
  • 刊出日期:  2023-04-24

目录

/

返回文章
返回