图1 BP神经网络的拓扑构架流程图
纸质出版日期:2020-09-25,
收稿日期:2019-07-27
扫 描 看 全 文
引用本文
阅读全文PDF
以乐昌峡鹅公带滑坡体作为研究对象,考虑日降雨量、渗压对边坡变形的影响,建立了BP、SVM、PSO-BP、PSO-SVM四种滑坡体变形预测模型。从乐昌峡安全检测系统导出近4年研究数据,筛选使用其中410组数据进行训练,取30组变形位移作为输出,分析后发现PSO-SVM模型为最佳模型。以PSO-SVM模型为研究对象,对粒子群算法迭代次数、种群规模、速度位置相关系数(k)等因素进行研究,得知三者分别为100、30、0.5时得到最优的PSO-SVM模型,此时的RMSE、MAPE、R2分别为0.202 mm、0.589%、0.985。相对于大型有限元仿真软件、多元线性回归模型等传统方法,文章所提出的预测模型可以减少计算成本;在面对非线性问题时也能够获得更好的处理效果。
Taking Lechangxia Egongdai landslide as the research object, the influence of daily rainfall and osmotic pressure on slope deformation is considered. By establishing BP, SVM, PSO-BP, PSO-SVM four landslide body deformation prediction models, the research data of the last 4 years is derived from the Lechangxia safety inspection system, and 410 sets of data are used for training through screening, and 30 sets of deformation displacements are taken as an output,after analysis, the PSO-SVM model is found to be the accurate model. Taking the PSO-SVM model as the basic model, the factors such as the number of iterations of the particle swarm algorithm, the population size, and velocity position correlation coefficient (k) are studied, and the best PSO-SVM is obtained when the three are 100, 30, and 0.5, respectively. In this model, the RMSE, MAPE, and R2 are 0.202 mm, 0.589%, and 0.985, respectively. Compared with traditional methods such as large-scale finite element simulation software and multiple linear regression models, the prediction model proposed in this article can reduce the computational cost and obtain better processing results in the face of nonlinear problems. At the same time, it can reduce the lack of fitting accuracy caused by incomplete factor analysis.
在水电站大坝的蓄水和施工建设阶段,受所在区域地质条件等因素的影响,加上工人施工过程中的操作不当,往往会对坝区周边稳定性欠佳的地质环境造严重成破坏,进而引发滑坡现象[
目前,对滑坡体位移的预测主要是利用各种边坡稳定性有限元软件[
李桥等[
虽然,许多学者对PSO优化算法进行了改进,并将这些算法应用于相关模型参数优化中 [
因此,本文考虑降雨、渗压对乐昌峡鹅公带滑坡体变形的影响,从乐昌峡安全检测系统导出近4年研究数据,筛选使用其中410组数据进行训练,选取30组变形位移作为预测对象,并在此基础上构建了四种神经网络模型。并基于PSO算法对SVM进行了优化,利用PSO算法在解空间中寻找到最佳的粒子,从而得到最佳惩罚参数c及核函数参数g。最后,通过对PSO迭代步数、种群规模、k等因素进行研究,选出最优迭代步数、种群规模、k数值,进一步优化了PSO-SVM模型。结果表明,这一模型的计算精度更高,是一种可靠的滑坡体变形预测方法,能够为大坝安全运行监测提供有价值的参考。
1.1.1 BP神经网络简述
BP神经网络可以被视为一种典型的前馈神经网络模型,由输入层、隐藏层以及输出层等三部分组成(
图1 BP神经网络的拓扑构架流程图
Fig. 1 Topology framework flow chart of BP neural network
1.1.2 粒子群算法简述
粒子群优化算法(particle swarm optimization, PSO)是许多学者经过大量的研究后提出的,它所对应的模型基本原理相对简单,并且也未提出目标函数可微的要求[
当采用PSO算法时,关于潜在解的问题都可以被视为搜索空间的粒子。对于任何一个粒子而言,在优化函数的作用下它都会形成一个适应度值,同时还对应有一个额外的速度,它决定了粒子们飞的方向以及距离。粒子位置更新方法如
图2 每代粒子位置的更新方式
Fig. 2 Updating method of particle location for each generation
我们可以假设在D维的目标搜索空间内部,里面有N个粒子组成的群体,此时第i个离子则可以被视为D维向量,有:
Xi=(xi1,xi2,…,xiD), i=1,2,…,N | (1) |
第i个v也可以表示为D维向量,表示为:
Vi=(vi1,vi2,…,viD), i=1,2,…,N | (2) |
第i个粒子寻找到的最优解作为单体的极值,为:
pi=(pi1,pi2,…,piD), i=1,2,…,N | (3) |
所有粒子一起找到的整体最优解被称为全局的极值,表示为:
gbest=(pg1,pg2,…,pgD) | (4) |
粒子速度和位置更新公式为:
Vj+1i=ω×Vji+c1×rand( )×(pbest-xji)+c2×rand( )×(gbest-xji) | (5) |
xj+1i=xji+vji | (6) |
其中,N指的便是种群规模;ω指的是惯性权重; Vji指的是粒子i经过多次迭代且在第j次迭代中所对应的速度;rand () 指的是介于0至1区间内的随机数; xji指的是粒子i经过多次迭代且在第j次迭代中所对应的位置;c1和c2均指的是学习因子; pbest指的是粒子i所对应的个体最优值;最后,gbest指的是群体最优值。
为了防止粒子的盲目搜索,一般建议将其位置和速度限制在一定的区间[-Xmax,Xmax]、[-Vmax,Vmax]。则粒子的飞行速度的取值范围为[-Vmax,Vmax]= k×[-Vmax,Vmax],若k取值太小,粒子飞行速度过慢,导致全区搜索能力过低,算法计算效率低;k取值范围过大,算法则存在着粒子容易飞出边界或飞过最优解的问题。借鉴刘欣蔚等[
min R (ω,ξi,ξ*i)=12ω×ω+cn∑i=1(ξi+ξ*i) | (7) |
其约束条件为:
{y-(ω×Φ(x)+b)<ε+ξi(ω×Φ(x)+b)-yi<ε+ξ*i ξi, ξ*i≥0, i=1,2,…,n | (8) |
(7)式中等号右侧的首项指的是泛化能力,紧接着的第二项则为惩罚项; c 为惩罚参数;ξi,ξi *为损失函数(当预测存在误差时都大于 0,否则等于0)。
在拉格朗日乘子的作用下,可以实现对原问题的处理,进一步转换为对偶优化问题,经相关求解后可以得到回归方程,具体有:
f(x)=W×F(x)+b=∑ni=1(αi-α*i)K(x,xi)+b | (9) |
K(x,xi)=exp (-‖x-z‖2g2) | (10) |
αi和α*i为Language 乘子,对应点即为支持向量。
本文考虑降雨、渗压对边坡变形的影响,以日降雨量、渗压作为输入,位移变形作为输出。通过对PSO迭代步数、种群规模、k等参数寻优,然后通过最优PSO算法来搜索获取最优的c和g以提高SVM预测结果的准确性。基于PSO-SVM模型的滑坡体位移预测流程见
图3 滑坡体位移预测流程图
Fig. 3 Flow chart of landslide displacement prediction
(1)展开有关于初始化粒子群的参数设置工作,主要涉及粒子位置、速度以及种群规模。
(2)计算每个粒子的适应度值Fit[i]。
(3)使用所有粒子的适应度值Fit[i]和个体的极端值pbest(i)进行比较。若Fit[i]>pbest(i),则用Fit[i]替换掉pbest(i)。
(4)围绕粒子适应度值Fit[i]和全局极值gbest(i)展开对比分析。如果出现Fit[i]>gbest(i)的情况,则意味着满足Fit[i]替换掉gbest(i)的条件。
(5)依据式(
(6)如满足,获得优化的c和g;否则返回(2)步骤。
(7)对结果进行误差分析。对迭代步数、种群规模、k参数进行寻优。
(8)寻找到优化参数,即建立满足条件的滑坡体位移预测模型。
乐昌峡水利枢纽工程是当前乐昌市境内最为重要的水利项目之一,它与乐昌市区的距离仅为14 km,且与韶关市相连,它在整个北江上游防控工作中发挥出极为重要的作用。该水利工程以防洪控制为主,还兼具发电与航运等多种应用功能。
参考相关工程资料得知,坝址上游库区多年来均存在岸坡体稳定性不足的问题,并且鹅公带滑坡体最为严重。此部分的坡体总方量达到240万m3,对平面布置图进行分析得知其主要呈椭圆形,长度达到400 m,平均宽度达到240 m;区域周边的山脊边坡陡度较大,但地形总体来说较为完整,上部后缘陡度普遍达到45°。对滑坡体的地质进行分析,主要以全风化粉砂质板岩为主,且含有大量的强风化岩块,在近几年多次发生小型滑坡现象。乐昌峡地区部分滑坡见
图 4 乐昌峡地区部分滑坡图
Fig. 4 Some landslide maps in Lechangxia
图5 鹅公带滑坡体平面示意
Fig. 5 Plane indication of landslide body in Egongdai
本文以乐昌峡鹅公带滑坡体作为研究对象,考虑1-1,剖面滑坡体日降雨量及P3点渗压对P3点边坡变形的影响。建立了BP、SVM、PSO-BP、PSO-SVM四种滑坡体变形预测模型,迭代次数为100,此外种群规模与k值分别对应为20与0.6。从乐昌峡安全检测系统导出近4年研究数据,筛选使用其中410组数据进行训练,对30组变形位移作为输出,相关数据见
样本 | 日降雨/mm | 渗压/m | 位移变形/mm | 样本 | 日降雨/mm | 渗压/m | 位移变形/mm |
---|---|---|---|---|---|---|---|
1 | 0.000 | 233.000 | 0.660 | 16 | 0.000 | 234.047 | 0.172 |
2 | 10.500 | 233.908 | 0.667 | 17 | 0.000 | 234.307 | 0.245 |
3 | 32.500 | 233.941 | 0.572 | 18 | 11.000 | 234.430 | 0.339 |
4 | 11.000 | 233.891 | 0.169 | 19 | 0.000 | 234.386 | 0.110 |
5 | 24.500 | 233.833 | 0.433 | 20 | 0.000 | 234.321 | 0.110 |
6 | 0.000 | 233.861 | 0.393 | … | … | … | … |
7 | 2.000 | 233.853 | 0.350 | 432 | 1.000 | 237.212 | 0.683 |
8 | 1.000 | 233.838 | 0.397 | 433 | 0.000 | 235.288 | 1.154 |
9 | 1.500 | 233.826 | 0.401 | 434 | 0.000 | 234.611 | 1.677 |
10 | 5.500 | 233.846 | 0.361 | 435 | 0.000 | 234.576 | 1.122 |
11 | 16.000 | 233.866 | 0.386 | 436 | 0.000 | 235.416 | 1.599 |
12 | 11.500 | 233.863 | 0.266 | 437 | 0.500 | 234.403 | -4.113 |
13 | 23.000 | 233.843 | 0.276 | 438 | 1.500 | 233.605 | -4.044 |
14 | 7.500 | 233.835 | 0.470 | 439 | 0.500 | 234.497 | -3.408 |
15 | 4.000 | 233.869 | 0.198 | 440 | 0.000 | 234.455 | -3.527 |
为了进一步量化各个模型之间的可靠性和准确性,使用三种指标展开评估工作,具体如下:
RMSE=√1nn∑i=1(Di-di)2 | (11) |
MAPE=1nn∑i=1|Di-didi|×100% | (12) |
R2=1-∑ni=1(Di-di)2∑ni=1(¯di-di)2 | (13) |
式中, di指的是实际值,¯di指的是实际值所对应的算术平均值;Di为预测值;n为预测样本的个数。
本文通过建立BP、PSO-BP、SVM、PSO-SVM等四种模型来进行P3点位移预测,迭代次数为100,此外种群规模与k值分别对应为20与0.6。
从
图 6 神经网络预测结果
Fig. 6 Neural network prediction result
(a)BP prediction model、(b)PSO-BP prediction model(c)SVM prediction model、(d)PSO-SVM prediction model
由
在本文的分析中,首先对PSO算法展开了迭代步数寻优操作,使其种群规模为20,k为0.6,迭代步数从10 -1 000之间进行寻找最佳迭代步数。然后,选择最佳迭代步数,令种群规模在10-200之间进行寻优。最后,选择最佳迭代步数和种群规模,探索最佳k值,选出最优PSO-SVM模型。
在对迭代次数进行研究后,我们从
图 7 计算时间和结果随迭代次数的变化
Fig. 7 Calculation time and results change with the number of iterations
当迭代次数与k值分别为100与0.6时,此时要求种群规模介于10~200之间,在此基础上展开寻优工作。
如
图 8 计算时间和结果随种群规模的变化
Fig. 8 Calculation time and results RMSE, MAPE, R2 as a function of population size
基于对
图 9 计算时间和结果的RMSE、MAPE、R2随k变化图
Fig. 9 RMSE, MAPE, R2 under different calculation time and results change with k
从图(b) 曲线可以看出RMSE、MAPE在0~0.5之间呈减小趋势,在0.5~1之间呈增长趋势。R2曲线变化则与之相反。不难看出,当k为0.5时,PSO-SVM模型预测准确率最高。此时,RMSE、MAPE、R2分别为0.202 mm、0.589%、0.985。
通过上述研究,发现迭代次数、种群规模、k分别为100、30、0.5时得到最优PSO-SVM模型,此时RMSE、MAPE、R2分别为0.202 mm、0.589%、0.985,如
图 10 PSO‒SVM模型参数优化预测对比图
Fig. 10 PSO-SVM model parameter optimization prediction comparison chart
针对滑坡体的变行预测问题,本文建立了BP、PSO-BP、SVM、PSO-SVM四种模型,在此基础上开展了参数寻优工作,由此得到了PSO-SVM模型。通过参数优化分析,明确了三大参数与模型变化之间的相关关系。进行了PSO算法的迭代次数、种群规模、k参数分析,由此确定出最佳值,加强了PSO 对SVM模型c、g参数的寻优功能,进一步提高了预测结果的准确性。最终的计算结果表明,参数寻优后的PSO-SVM模型,优于BP、SVM、PSO-BP模型,是一种有效的滑坡体位移预测新方法,为目前滑坡体监测困难提供了一条简便的监测方法。
李明,曾金志. 基于ARIMA模型的大坝滑坡体表面位移监测分析预报[J]. 东华理工大学学报(自然科学版), 2016(S1): 17-20. [百度学术]
LI M, ZENG J Z. Monitoring, analysis and forecast of surface displacement of dam landslide based on ARIMA model [J]. Journal of East China University of Technology: Natural Science, 2016(S1): 17-20. [百度学术]
DAI F C, LEE CF, NGAI Y Y. Landslide risk assessment and management: an overview[J]. Engineering Geology, 2002,64(1): 65-87. [百度学术]
JIANG J W, EHRET D,XIANG W,et al. Numerical simulation of Qiaotou Landslide deformation caused by drawdown of the three gorges reservoir, China[J]. Environmental Earth Sciences, 2011, 62(2): 411-419. [百度学术]
LEE K T, HO J Y. Prediction of landslide occurrence based on slope-instability analysis and hydrological model simulation[J]. Journal of Hydrology, 2009,375(3): 489-497. [百度学术]
CHENG L,ZENG Z G,YAO W, et al. Multiple neural networks switched prediction for landslide displacement[J]. Engineering Geology, 2015,186: 91-99. [百度学术]
ZENGY,YAN E C,LI C F, et al. Application of multivariable time series based on RBF neural network in prediction of landslide displacement[C]// Proceedings of 2008 International Workshop on Chaos-Fractals Theories and Applications & the 9th International Conference for Young Computer Scientists,2008:15-20. [百度学术]
CAI Z L,XU W Y, MENG Y D, et al. Prediction of landslide displacement based on GA-LSSVM with multiple factors[J]. Bulletin of Engineering Geology & the Environment, 2016,75(2): 637-646. [百度学术]
白航, 王伟竹,姚清河. 某水下机器人阻力特性的数值模拟[J].中山大学学报(自然科学版), 2018,57(1): 150-159. [百度学术]
BAI H,WANG W Z,YAO Q H. Numerical simulation of the resistance character of an underwater vehicle[J]. Acta Scientiarum Naturalium Universitatis Sunyatseni, 2018,57(1): 150-159. [百度学术]
冯立文,姚清河,苏炜等人. 洪湾渔港工程潮流泥沙的输运特性[J]. 中山大学学报(自然科学版), 2019, 58(1): 56-64. [百度学术]
FENG L W,YAO Q H,SU W,et al.Transportation characteristic analysis of tidal current sediment in Hongwan fishing port[J]. Acta Scientiarum Naturalium Universitatis Sunyatseni, 2019,58(1): 56-64. [百度学术]
李桥, 巨能攀, 黄健,等.基于EEMD与SE的IPSO-LSSVM耦合模型在坝肩边坡变形预测中的应用[J]. 长江科学院院报, 2019,36(12):47-53. [百度学术]
LI Q, JU N P, HUANG J, et al.A prediction method for abutment slope deformation based on IPSO-LSSVM model integrating ensemble empirical mode decomposition and sample entropy[J]. Journal of Yangtze River Scientific Research Institute, 2019,36(12):47-53. [百度学术]
魏博文,熊威,李火坤,等.融合混沌残差的大坝位移蛙跳式组合预报模型[J].武汉大学学报(信息科学版),2016,41(9):1272-1278. [百度学术]
WEI B W, XIONG W, LI H K, et al., Dam Deformation forecasting of leapfrog combined model merging residual errors of chaos[J]. Geomatics & Information Science of Wuhan University, 2016,41(9):1272-1278. [百度学术]
LIAN C,ZENG Z G,YAO W, et al. Displacement prediction model of landslide based on a modified ensemble;empirical mode decomposition and extreme learning machine[J]. Natural Hazards, 2013,66(2): 759-771. [百度学术]
彭令, 牛瑞卿, 吴婷. 时间序列分析与支持向量机的滑坡位移预测[J]. 浙江大学学报(工学版), 2013,47(9): 1672-1679. [百度学术]
PENG L,NIU R Q,WU T. Time series analysis and support vector machine for landslide displacement prediction[J]. Journal of Zhejiang University(Engineering Science),2013,47(9): 1672-1679. [百度学术]
ZHOU C,YIN K L,CAO Y, et al.A novel method for landslide displacement prediction by integrating advanced computational intelligence algorithms[J]. Scientific Reports, 2018,8(1):47-53. [百度学术]
易庆林, 曾怀恩,黄海峰. 利用BP神经网络进行水库滑坡变形预测[J]. 水文地质工程地质, 2013,40(1): 124-128. [百度学术]
YI Q L,ZENG H E,HUANG H F. Reservoir landslide deformation forecast using BP neural network[J].Hydrogeology and Engineering Geology, 2013,40(1): 124-128. [百度学术]
MIAO F S , WU Y , XIE Y P , et al. Research on progressive failure process of Baishuihe landslide based on Monte Carlo model[J]. Stochastic Environmental Research and Risk Assessment, 2017,31:1683-1696. [百度学术]
KANGA F, LIA J J, ZHAOA S Z, et al. Structural health monitoring of concrete dams using long-term air temperature for thermal effect simulation[J]. Engineering Structures , 2019,180: 642–653. [百度学术]
张文明, 董增川, 朱成涛,等.基于粒子群算法的水文模型参数多目标优化研究[J]. 水利学报, 2008,39(5): 20-26. [百度学术]
ZHANG W M,DONG Z C,ZHU C T, et al. Automatic calibration of hydrologic model based on multi-objective particle swarm optimization method[J].Journal of Hydraulic Engineering, 2008,39(5): 20-26. [百度学术]
王问宇,袁鹏,邵骏,等. 改进粒子群算法在新安江模型参数优选中的应用[J]. 水力发电, 2008,34(8): 14-16. [百度学术]
WANG W Y,YUAN P,SHAO J, et al. Application of improved particle swarm algorithm to the parameter selection optimization of Xin'anjiang Model[J]. Journal of Hydropower, 2008, 34(8): 14-16. [百度学术]
齐红,任洪娥,贾鹤鸣,等.基于改进PSO神经网络的微米木纤维切削参数优化[J]. 中山大学学报(自然科学版), 2018,57(2): 116-122. [百度学术]
QI H,REN H E,JIA H M, et al. Optimization study of cutting parameters of micron wood fiber based on improved PSO neural network[J]. Acta Scientiarum Naturalium Universitatis Sunyatseni, 2018,57(2): 116-122. [百度学术]
SHI Y , EBERHART R C . Empirical study of particle swarm optimization[C]// Congress on Evolutionary Computation. IEEE, 2002. [百度学术]
王维博, 林川,郑永康. 粒子群算法中参数的实验与分析[J]. 西华大学学报(自然科学版), 2008,27(1): 76-80. [百度学术]
WANG W B,LING C, ZHENG Y K. Experiment and analysis of parameters in particle swarm optimization[J].Journal of Xihua University(Natural Science Edition), 2008,27(1): 76-80. [百度学术]
丁根宏,曹文秀. 改进粒子群算法在水库优化调度中的应用[J]. 南水北调与水利科技, 2014,12(1): 118-121. [百度学术]
DING G H,CAO W X. Application of improved particle swarm optimization algorithm in optimal operation of reservoir[J]. South-to-North Water Transfers and Water Science & Technology, 2014,12(1): 118-121. [百度学术]
KENNEDY J , EBERHART R . Particle swarm optimization[C]// Proceedings of ICNN'95 - International Conference on Neural Networks. IEEE, 1995. [百度学术]
43
浏览量
59
下载量
3
CSCD
相关文章
相关作者
相关机构