CN109272721A - 基于kpca-foa-lssvm的滑坡灾害预报方法 - Google Patents
基于kpca-foa-lssvm的滑坡灾害预报方法 Download PDFInfo
- Publication number
- CN109272721A CN109272721A CN201811294315.6A CN201811294315A CN109272721A CN 109272721 A CN109272721 A CN 109272721A CN 201811294315 A CN201811294315 A CN 201811294315A CN 109272721 A CN109272721 A CN 109272721A
- Authority
- CN
- China
- Prior art keywords
- landslide
- formula
- value
- disaster
- equation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G08—SIGNALLING
- G08B—SIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
- G08B31/00—Predictive alarm systems characterised by extrapolation or other computation using updated historic data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2135—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on approximation criteria, e.g. principal component analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
-
- G—PHYSICS
- G08—SIGNALLING
- G08B—SIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
- G08B21/00—Alarms responsive to a single specified undesired or abnormal condition and not otherwise provided for
- G08B21/02—Alarms for ensuring the safety of persons
- G08B21/10—Alarms for ensuring the safety of persons responsive to calamitous events, e.g. tornados or earthquakes
Landscapes
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Biology (AREA)
- Business, Economics & Management (AREA)
- Emergency Management (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- Computing Systems (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
本发明公开的基于KPCA‑FOA‑LSSVM的滑坡灾害预报方法,首先建立滑坡体实时监测预警系统,获取监测区的实时数据,将其标准化处理后,运用核主成分分析法筛选出滑坡发生主要影响因素作为输入变量;构建基于LSSVM滑坡灾害预报模型;然后运用果蝇算法进行参数寻优,更新网络参数;最后重构已优化的滑坡灾害预报模型,输出滑坡发生概率对应发生等级,完成预报。本发明公开的方法通过建立滑坡体实时监测预警系统获取监测数据,经过核主成分分析筛选出主要影响因素后利用基于果蝇算法优化的最小二乘向量机模型中训练输出滑坡发生概率,提升预报效率,增加精确度。
Description
技术领域
本发明属于地质灾害监测技术领域,涉及一种基于KPCA-FOA-LSSVM的滑坡灾害预报方法。
背景技术
滑坡是一种常见的自然灾害现象,它的出现能够使农田与房舍被摧毁、破坏道路与水利水电设施,从而导致停电、停水、停工等后果,严重时还可能威胁到人类的生命安全。因此,如何采取及时而准确的滑坡灾害预报方法成为人们所关注的热点。
针对滑坡发生的随机性、高频率、影响范围广等特点,国内外学者对其进行了深入的研究,但使用的方法各有利弊。例如灰色模型,无需大量训练样本且计算工作量小,但由于其内在机理复杂,难以分析,导致准确率不够稳定;神经网络模型具有较好的自主学习与容错能力,但其收敛速度较慢,时间过长容易陷入局部极值;朴素贝叶斯判别模型,对于多分类问题处理效果好,但在分析较为复杂的样本时误判率会上升;支持向量机(SVM)模型,具有良好的泛化能力,在训练样本时易于寻求全局最优解,但缺点是其训练机理为一个有约束的二次规划问题,在样本容量大时,可能导致训练时间过长,效率低下。而最小二乘支持向量机(LSSVM)的出现,改变了标准SVM的约束条件与风险函数,在继承支持向量机优点的同时,还能够很好地解决数据量少、过学习以及高维数等实际问题,提升了训练效率与准确性,在预测预报领域中体现出较大优势。
鉴于此,本发明提出一种基于KPCA-FOA-LSSVM的滑坡灾害预报方法,首先运用核主成分分析法(Kernel principal components analysis,KPCA)筛选出研究区滑坡发生主要影响因素,将标准化处理后的数据输入于最小二乘支持向量机(Least square supportvector machines,LSSVM)滑坡预测模型中进行运算,然后利用果蝇优化算法(Fruit FlyOptimization Algorithm,FOA)参数优化,更新网络参数,最后输出滑坡发生概率对应发生等级完成滑坡预报。
发明内容
本发明的目的是提供一种基于KPCA-FOA-LSSVM的滑坡灾害预报方法,解决现有灾害预报方法采用的算法过学习、效率低、精确度不高的问题。
本发明所采用的技术方案是,基于KPCA-FOA-LSSVM的滑坡灾害预报方法,具体过程包括如下步骤:
步骤1.建立滑坡体实时监测预警系统,获取监测区滑坡的实时数据,将其标准化处理后,运用核主成分分析法筛选出滑坡发生主要影响因素作为输入变量;
步骤2.构建基于最小二乘支持向量机滑坡灾害预报模型;
步骤3.运用果蝇算法进行参数寻优,更新网络参数;
步骤4.根据步骤3优化的结果,重构滑坡灾害预报模型,输出滑坡发生概率对应的发生等级,实现滑坡灾害的实时监测。
本发明的其他特点还在于,
步骤1中获取监测区滑坡的实时数据,将其标准化处理后,运用核主成分分析法筛选出滑坡发生主要影响因素作为输入变量的具体过程如下:
步骤1.1分析不同的监测区存在不同的滑坡灾害影响因素,根据地质环境背景与历史灾害发生记载情况,确定初始影响因素为:降雨量、土壤含水率、裂缝位移值、孔隙水压力、次声频率、坡向、坡度、高程、平面曲率、地层振动指数、归一化植被系数和岩组工12个影响因素;将滑坡体监测预警系统采集的实时数据作为原始数据集;
步骤1.2将步骤1.1获取的原始数据标准化处理,如公式1所示:
其中,x为标准化后的值,xOriginal为输入的实际值,xmin和xmax分别为所有输入值中的最小值和所有输入值中的最大值;
步骤1.3采用核主成分分析法筛选出滑坡发生的主要影响因素,具体过程如下:
设训练样本集为X={x1,x2,…xn},其中xi∈RP,yi=Rp(RP为输入空间,P为数据维数,i=1,…n,n为训练样本总数),输入空间非线性映射X→F,则对应的协方差矩阵C如式2所示:
其中,C为协方差矩阵,n为训练样本数,为输入空间非线性映射量;
对协方差矩阵C进行特征分解,如式3所示:
λνi=Cλ (3)
其中,λ≥0,且特征向量如式4所示:
由于特征向量ν是由非线性映射空间组成,因此,式3与式5等价:
将式2、式4代入式5中,且令其中K为核矩阵,得到式6:
Kαn=nλiαn (6)
其中,核矩阵的特征向量为α1,α2…αn,特征值为nλi;
选取前m(m<n)个特征值对应的归一化特征向量α1,α2…αm,其中对在νr投影如式7所示:
其中,j=1,2…,n;r=1,2,…,m;gr(xj)为对应于的第r个非线性主元分量;
令所有投影值g(xj)=(g1(xj),g2(xj)…gm(xj))作为样本特征值,利用核函数代替空间点积计算,则式7转换为式8所示:
根据μj贡献率选取主要成分Tj,如公式9所示:
其中,i为主元成分个数,μi为主成分贡献率,λi为第i个特征值,为总特征值;
式9中取累计方差贡献率μj≥85%作为主要滑坡灾害影响因子数据,因此,筛选出的主要灾害影响因素为:x1=降雨量,x2=土壤含水率,x3=裂缝位移量,x4=孔隙水压力,x5=次声频率,x6=岩组,x7=归一化植被系数。
步骤2中构建基于最小二乘支持向量机的具体过程如下:
步骤2.1将步骤1中筛选出的7个主要影响因素作为模型输入部分,设其中共有K个样本,将样本设置为8:1的分配比例,分别为训练样本和测试样本;且将模型输出设置为一个5维向量,其输出量为:Y1=成灾概率0~20%,Y2=成灾概率20~40%,Y3=成灾概率40~60%,Y4=成灾概率60~80%,Y5=成灾概率80~100%,且各分量的值分别为1或-1;
步骤2.2构建基于最小二乘支持向量机的滑坡灾害预报模型:
对于步骤2.1中给定的样本作非线性映射φ:Rn→F,设定训练样本集(xi,yi),i=1,…,k,k为样本总数;xi∈Rn,yi∈R,Rn为输入空间,n为数据维数;
因此,被估计函数表达式如式10所示:
y=f(x)=wTφ(x)+b (10)
其中,w和b均为结构风险化最小化模型参数,且w代表空间F的权向量,b∈R代表为偏差量;f(x)为被估计函数;
在确定决策参数w、b时,求解问题如式11所示:
其中,ξk∈R为误差变量,γ为可调超参数;
由式11定义拉格朗日函数形式如式12所示:
其中,αk∈R为拉格朗日乘子;
对式12中的各变量求偏导并整理为线性方程组,得到式13:
其中,Il=[1,1,…,1];Ωij=K(xi,xj)=φ(xi)Tφ(xj),K(xi,xj)为核函数;α=[α1,α2…αl],i,j=1,2,…l;
在此,核函数选择RBF高斯径向基核函数,其表达式如式14所示:
由式13中得出最小二乘支持向量机的函数估计如式15所示:
其中,不为零的元素αi所对应的样本(xi,yi)为支持向量;
由式13和式15得到第m等级的分类功能如式16所示:
其中,如果分类结果属于第m类(m=1,2,3,4,5),则ymi=1,否则ymi=-1;fm(x)为分类出的第m个等级;
如果对应的各类的分类函数,即公式16的输出结果是1,则属于该等级;如果输出结果是-1,则不属于该等级。
步骤3中运用果蝇算法进行参数寻优,更新网络参数的具体过程如下:
步骤3.1初始化果蝇算法参数,包括种群规模、迭代次数以及果蝇初始位置,针对最小二乘支持向量机模型中需要优化的两个参数:可调超参数γ和核参数σ,选取样本中随机的两个数,令初始坐标分别为
步骤3.2赋予每只果蝇随机距离和随机方向,定义坐标为 且计算果蝇个体与坐标原点的倒数值,计算出味道浓度判定值
步骤3.3确定需要优化的参数范围即γi∈(0,1000],σ∈(0,100];
步骤3.4将步骤1.3中已经归一化理过的训练样本作为训练样本,令分类准确率作为适应度函数如式17所示:
Smelli=Fitness(Ci,σi)=accuracy(Ci,σi) (17)
其中,Smelli代表果蝇个体位置的味道浓度;
寻找Fitness的最大值对应分类精确值最大的果蝇个体,进入迭代寻优,判断最高分类精确值是否由于前一代,若高于,则保留精确最大值以及对应的坐标值,并更新初始坐标的值;反之,返回步骤3.2继续运算,直至寻求最优为止。
本发明的有益效果是,基于KPCA-FOA-LSSVM的滑坡灾害预报方法,解决现有灾害预报方法采用的算法过学习、效率低、精确度不高的问题。相比于现有方法具有以下优势:
(1)建立滑坡体监测预警系统,实时采集数据,保证时效性;通过KPCA方法筛选滑坡发生主要影响因素,削弱次要因素,并且避免了特征向量线性不可分问题,防止维数灾难。
(2)利用LSSVM模型对灾害影响因素样本特征学习,构造最优决策函数,简化计算方法,提升预报效率;
(3)利用果蝇算法进行全局优化,保证网络整体最优,提升预报精确度。
附图说明
图1是本发明的基于KPCA-FOA-LSSVM的滑坡灾害预报方法所使用的滑坡灾害实时监测预警系统结构示意图;
图2是本发明的基于KPCA-FOA-LSSVM的滑坡灾害预报方法的流程图;
图3是最小二乘支持向量机体系图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明的滑坡灾害实时监测装置如图1所示,包括设置在监测区域的环境因子采集模块,环境因子采集模块将采集的数据传输给滑坡预警遥测终端(RTU),滑坡预警遥测终端设置在监测区的坡体上,滑坡遥测终端的内部设置有电源模块、无线通信模块和报警模块,滑坡遥测终端的外部安装有门禁锁装置,门禁锁为非接触式IC一卡通,采用射频识别技术使得管理员通过一卡通出入滑坡遥测终端;滑坡预警遥测终端通过内部的无线通信模块将环境因素采集模块采集的数据传送到中心站,中心站将收集的数据导入相应预报模型中进行处理并判断预警结果,通过无线通信的方式上传至客户端。
环境因子采集模块包括设置在监测区域坡体上的雨量计、测缝计、渗压计、倾斜计和次声计。雨量计用于获取监测坡体区域内一段时间的累计降雨量;测缝计用于测量坡体的伸缩缝或裂缝变化的位移量;渗压计用于实时监测地质区域内土体所承受的压力值;倾斜计用于实时测量坡体发生倾斜角度的变化量;次声计用于获取坡体固体颗粒之间碰撞产生次声波的变化量。
中心站包括服务器、预警系统、数据库系统、数据分析系统和防火墙,数据库系统用于储存滑坡预警遥测终端传送的实时数据,数据库系统将存储的数据传送到数据分析系统,数据分析系统将数据进行判断分析后发送到服务器进行计算,预警系统的作用是将预警信息上报至客户端,防火墙的是实现数据保密与完整性的重要方式,对访问控制数据的安全性管理提供保证。
电源模块采用锂电池和太阳能光伏发电板供电。当光照充足时,太阳能光伏发电板为整个系统供电;当光照不足或光伏发电板提供的直流电压小于12V时,选取备用的12V锂电池为整个系统供电。电源模块为滑坡预警遥测终端提供电源能力,保证其能够实时监测。
无线通信模块采用无线wifi网络,负责将滑坡预警遥测终端的监测的实时数据上传到中心站,和接收中心站下发的命令进行无线传输。
报警模块采用爆闪灯,当中心站根据得到的滑坡实时数据处理后的结果,若结果表示为警情,中心站自动下发报警信号给滑坡预警遥测终端的报警模块和客户端,现场爆闪灯闪烁进行报警提示。
滑坡灾害实时监测装置还包括设置在坡体现场的视频监控模块,负责监控坡体现场实时的状况以及安全防盗,视频监控模块通过无线wifi将视频数据上传至中心站。
客户端是PC机客户端或者手机客户端,管理者通过在PCJ机或手机上登录其帐号,能够观测滑坡预警遥测终端采集到的数据量以及中心站处理发出的预警结果。同时,管理者能够通过PC机客户端或手机客户端下达指令控制滑坡预警遥测终端的开启与关闭。
本发明的基于KPCA-FOA-LSSVM的滑坡灾害预报方法,如图2所示,具体过程包括如下步骤:
步骤1.建立滑坡体实时监测预警系统,获取监测区滑坡的实时数据,将其标准化处理后,运用核主成分分析法筛选出滑坡发生主要影响因素作为输入变量;
步骤2.构建基于最小二乘支持向量机滑坡灾害预报模型;
步骤3.运用果蝇算法进行参数寻优,更新网络参数;
步骤4.根据步骤3优化的结果,重构滑坡灾害预报模型,输出滑坡发生概率对应的发生等级,实现滑坡灾害的实时监测。
步骤1中获取监测区滑坡的实时数据,将其标准化处理后,运用核主成分分析法筛选出滑坡发生主要影响因素作为输入变量的具体过程如下:
步骤1.1分析不同的监测区存在不同的滑坡灾害影响因素,根据地质环境背景与历史灾害发生记载情况,确定初始影响因素为:降雨量、土壤含水率、裂缝位移值、孔隙水压力、次声频率、坡向、坡度、高程、平面曲率、地层振动指数、归一化植被系数和岩组工12个影响因素;将滑坡体监测预警系统采集的实时数据作为原始数据集;
步骤1.2将步骤1.1获取的原始数据标准化处理,如公式1所示:
其中,x为标准化后的值,xOriginal为输入的实际值,xmin和xmax分别为所有输入值中的最小值和所有输入值中的最大值;
步骤1.3采用核主成分分析法筛选出滑坡发生的主要影响因素,具体过程如下:
设训练样本集为X={x1,x2,…xn},其中xi∈RP,yi=Rp(RP为输入空间,P为数据维数,i=1,…n,n为训练样本总数),输入空间非线性映射X→F,则对应的协方差矩阵C如式2所示:
其中,C为协方差矩阵,n为训练样本数,为输入空间非线性映射量;
对协方差矩阵C进行特征分解,如式3所示:
λνi=Cλ (3)
其中,λ≥0,且特征向量如式4所示:
由于特征向量ν是由非线性映射空间组成,因此,式3与式5等价:
将式2、式4代入式5中,且令其中K为核矩阵,得到式6:
Kαn=nλiαn (6)
其中,核矩阵的特征向量为α1,α2…αn,特征值为nλi;
选取前m(m<n)个特征值对应的归一化特征向量α1,α2…αm,其中对在νr投影如式7所示:
其中,j=1,2…,n;r=1,2,…,m;gr(xj)为对应于的第r个非线性主元分量;
令所有投影值g(xj)=(g1(xj),g2(xj)…gm(xj))作为样本特征值,利用核函数代替空间点积计算,则式7转换为式8所示:
根据μj贡献率选取主要成分Tj,如公式9所示:
其中,i为主元成分个数,μi为主成分贡献率,λi为第i个特征值,为总特征值;
式9中取累计方差贡献率μj≥85%作为主要滑坡灾害影响因子数据,因此,筛选出的主要灾害影响因素为:x1=降雨量,x2=土壤含水率,x3=裂缝位移量,x4=孔隙水压力,x5=次声频率,x6=岩组,x7=归一化植被系数。主要影响因素输入模式如表1所示。
如图3所示,步骤2中构建基于最小二乘支持向量机的具体过程如下:
步骤2.1将步骤1中筛选出的7个主要影响因素作为模型输入部分,设其中共有K个样本,将样本设置为8:1的分配比例,分别为训练样本和测试样本;且将模型输出设置为一个5维向量,其输出量为:Y1=成灾概率0~20%,Y2=成灾概率20~40%,Y3=成灾概率40~60%,Y4=成灾概率60~80%,Y5=成灾概率80~100%,且各分量的值分别为1或-1;成灾概率与对应编码如表2所示:
步骤2.2构建基于最小二乘支持向量机的滑坡灾害预报模型:
表1主要影响输入模式
表2成灾概率及其编码
对于步骤2.1中给定的样本作非线性映射φ:Rn→F,设定训练样本集(xi,yi),i=1,…,k,k为样本总数;xi∈Rn,yi∈R,Rn为输入空间,n为数据维数;
因此,被估计函数表达式如式10所示:
y=f(x)=wTφ(x)+b (10)
其中,w和b均为结构风险化最小化模型参数,且w代表空间F的权向量,b∈R代表为偏差量;f(x)为被估计函数;
在确定决策参数w、b时,求解问题如式11所示:
其中,ξk∈R为误差变量,γ为可调超参数;
由式11定义拉格朗日函数形式如式12所示:
其中,αk∈R为拉格朗日乘子;
对式12中的各变量求偏导并整理为线性方程组,得到式13:
其中,Il=[1,1,…,1];Ωij=K(xi,xj)=φ(xi)Tφ(xj),K(xi,xj)为核函数;α=[α1,α2…αl],i,j=1,2,…l;
在此,核函数选择RBF高斯径向基核函数,这种核函数能够较好处理线性不可分情形,而且需要确定的参数少,无论在大样本还是小样本中均具有较好的性能,其表达式如式14所示:
由式13中得出最小二乘支持向量机的函数估计如式15所示:
其中,不为零的元素αi所对应的样本(xi,yi)为支持向量;
由式13和式15得到第m等级的分类功能如式16所示:
其中,如果分类结果属于第m类(m=1,2,3,4,5),则ymi=1,否则ymi=-1;fm(x)为分类出的第m个等级;
如果对应的各类的分类函数,即公式16的输出结果是1,则属于该等级;如果输出结果是-1,则不属于该等级。
步骤3中运用果蝇算法进行参数寻优,更新网络参数的具体过程如下:
步骤3.1初始化果蝇算法参数,包括种群规模、迭代次数以及果蝇初始位置,针对最小二乘支持向量机模型中需要优化的两个参数:可调超参数γ和核参数σ,选取样本中随机的两个数,令初始坐标分别为
步骤3.2赋予每只果蝇随机距离和随机方向,定义坐标为 且计算果蝇个体与坐标原点的倒数值,计算出味道浓度判定值
步骤3.3确定需要优化的参数范围即γi∈(0,1000],σ∈(0,100];
步骤3.4将步骤1.3中已经归一化理过的训练样本作为训练样本,令分类准确率作为适应度函数如式17所示:
Smelli=Fitness(Ci,σi)=accuracy(Ci,σi) (17)
其中,Smelli代表果蝇个体位置的味道浓度;
寻找Fitness的最大值对应分类精确值最大的果蝇个体,进入迭代寻优,判断最高分类精确值是否由于前一代,若高于,则保留精确最大值以及对应的坐标值,并更新初始坐标的值;反之,返回步骤3.2继续运算,直至寻求最优为止。
步骤4中将相对应的滑坡发生概率分为五个等级,分级情况如表3所示:
表3预警等级划分
Claims (4)
1.基于KPCA-FOA-LSSVM的滑坡灾害预报方法,其特征在于,具体过程包括如下步骤:
步骤1.建立滑坡体实时监测预警系统,获取监测区滑坡的实时数据,将其标准化处理后,运用核主成分分析法筛选出滑坡发生主要影响因素作为输入变量;
步骤2.构建基于最小二乘支持向量机滑坡灾害预报模型;
步骤3.运用果蝇算法进行参数寻优,更新网络参数;
步骤4.根据步骤3优化的结果,重构滑坡灾害预报模型,输出滑坡发生概率对应的发生等级,实现滑坡灾害的实时监测。
2.如权利要求1所述的基于KPCA-FOA-LSSVM的滑坡灾害预报方法,其特征在于,所述步骤1中获取监测区滑坡的实时数据,将其标准化处理后,运用核主成分分析法筛选出滑坡发生主要影响因素作为输入变量的具体过程如下:
步骤1.1分析不同的监测区存在不同的滑坡灾害影响因素,根据地质环境背景与历史灾害发生记载情况,确定初始影响因素为:降雨量、土壤含水率、裂缝位移值、孔隙水压力、次声频率、坡向、坡度、高程、平面曲率、地层振动指数、归一化植被系数和岩组工12个影响因素;将滑坡体监测预警系统采集的实时数据作为原始数据集;
步骤1.2将步骤1.1获取的原始数据标准化处理,如公式1所示:
其中,x为标准化后的值,xOriginal为输入的实际值,xmin和xmax分别为所有输入值中的最小值和所有输入值中的最大值;
步骤1.3采用核主成分分析法筛选出滑坡发生的主要影响因素,具体过程如下:
设训练样本集为X={x1,x2,…xn},其中xi∈RP,yi=Rp(RP为输入空间,P为数据维数,i=1,…n,n为训练样本总数),输入空间非线性映射X→F,则对应的协方差矩阵C如式2所示:
其中,C为协方差矩阵,n为训练样本数,为输入空间非线性映射量;
对协方差矩阵C进行特征分解,如式3所示:
λνi=Cλ (3)
其中,λ≥0,且特征向量如式4所示:
由于特征向量ν是由非线性映射空间组成,因此,式3与式5等价:
将式2、式4代入式5中,且令其中K为核矩阵,得到式6:
Kαn=nλiαn (6)
其中,核矩阵的特征向量为α1,α2…αn,特征值为nλi;
选取前m(m<n)个特征值对应的归一化特征向量α1,α2…αm,其中对在νr投影如式7所示:
其中,j=1,2…,n;r=1,2,…,m;gr(xj)为对应于的第r个非线性主元分量;
令所有投影值g(xj)=(g1(xj),g2(xj)…gm(xj))作为样本特征值,利用核函数代替空间点积计算,则式7转换为式8所示:
根据μj贡献率选取主要成分Tj,如公式9所示:
其中,i为主元成分个数,μi为主成分贡献率,λi为第i个特征值,为总特征值;
式9中取累计方差贡献率μj≥85%作为主要滑坡灾害影响因子数据,因此,筛选出的主要灾害影响因素为:x1=降雨量,x2=土壤含水率,x3=裂缝位移量,x4=孔隙水压力,x5=次声频率,x6=岩组,x7=归一化植被系数。
3.如权利要求2所述的基于KPCA-FOA-LSSVM的滑坡灾害预报方法,其特征在于,所述步骤2中构建基于最小二乘支持向量机的具体过程如下:
步骤2.1将步骤1中筛选出的7个主要影响因素作为模型输入部分,设其中共有K个样本,将样本设置为8:1的分配比例,分别为训练样本和测试样本;且将模型输出设置为一个5维向量,其输出量为:Y1=成灾概率0~20%,Y2=成灾概率20~40%,Y3=成灾概率40~60%,Y4=成灾概率60~80%,Y5=成灾概率80~100%,且各分量的值分别为1或-1;
步骤2.2构建基于最小二乘支持向量机的滑坡灾害预报模型:
对于步骤2.1中给定的样本作非线性映射φ:Rn→F,设定训练样本集(xi,yi),i=1,…,k,k为样本总数;xi∈Rn,yi∈R,Rn为输入空间,n为数据维数;
因此,被估计函数表达式如式10所示:
y=f(x)=wTφ(x)+b (10)
其中,w和b均为结构风险化最小化模型参数,且w代表空间F的权向量,b∈R代表为偏差量;f(x)为被估计函数;
在确定决策参数w、b时,求解问题如式11所示:
其中,ξk∈R为误差变量,γ为可调超参数;
由式11定义拉格朗日函数形式如式12所示:
其中,αk∈R为拉格朗日乘子;
对式12中的各变量求偏导并整理为线性方程组,得到式13:
其中,Il=[1,1,…,1];Ωij=K(xi,xj)=φ(xi)Tφ(xj),K(xi,xj)为核函数;α=[α1,α2…αl],i,j=1,2,…l;
在此,核函数选择RBF高斯径向基核函数,其表达式如式14所示:
由式13中得出最小二乘支持向量机的函数估计如式15所示:
其中,不为零的元素αi所对应的样本(xi,yi)为支持向量;
由式13和式15得到第m等级的分类功能如式16所示:
其中,如果分类结果属于第m类(m=1,2,3,4,5),则ymi=1,否则ymi=-1;fm(x)为分类出的第m个等级;
如果对应的各类的分类函数,即公式16的输出结果是1,则属于该等级;如果输出结果是-1,则不属于该等级。
4.如权利要求2所述的基于KPCA-FOA-LSSVM的滑坡灾害预报方法,其特征在于,所述步骤3中运用果蝇算法进行参数寻优,更新网络参数的具体过程如下:
步骤3.1初始化果蝇算法参数,包括种群规模、迭代次数以及果蝇初始位置,针对最小二乘支持向量机模型中需要优化的两个参数:可调超参数γ和核参数σ,选取样本中随机的两个数,令初始坐标分别为
步骤3.2赋予每只果蝇随机距离和随机方向,定义坐标为 且计算果蝇个体与坐标原点的倒数值,计算出味道浓度判定值
步骤3.3确定需要优化的参数范围即γi∈(0,1000],σ∈(0,100];
步骤3.4将步骤1.3中已经归一化理过的训练样本作为训练样本,令分类准确率作为适应度函数如式17所示:
Smelli=Fitness(Ci,σi)=accuracy(Ci,σi) (17)
寻找Fitness的最大值对应分类精确值最大的果蝇个体,进入迭代寻优,判断最高分类精确值是否由于前一代,若高于,则保留精确最大值以及对应的坐标值,并更新初始坐标的值;反之,返回步骤3.2继续运算,直至寻求最优为止。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811294315.6A CN109272721B (zh) | 2018-11-01 | 2018-11-01 | 基于kpca-foa-lssvm的滑坡灾害预报方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811294315.6A CN109272721B (zh) | 2018-11-01 | 2018-11-01 | 基于kpca-foa-lssvm的滑坡灾害预报方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109272721A true CN109272721A (zh) | 2019-01-25 |
CN109272721B CN109272721B (zh) | 2020-10-23 |
Family
ID=65191047
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811294315.6A Active CN109272721B (zh) | 2018-11-01 | 2018-11-01 | 基于kpca-foa-lssvm的滑坡灾害预报方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109272721B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109887240A (zh) * | 2019-03-22 | 2019-06-14 | 福州大学 | 一种基于人工智能的滑坡灾害安全监测与预警方法 |
CN110196448A (zh) * | 2019-05-21 | 2019-09-03 | 成都信息工程大学 | 一种滑坡次声信号识别方法 |
CN110456217A (zh) * | 2019-08-08 | 2019-11-15 | 辽宁工程技术大学 | 一种基于wpd-foa-lssvm双模型的mmc故障定位方法 |
CN110598914A (zh) * | 2019-08-21 | 2019-12-20 | 武汉理工大学 | 一种多因素影响下矿井灾害气体浓度区间预测方法及系统 |
CN111783866A (zh) * | 2020-06-28 | 2020-10-16 | 中国电子科技集团公司第二十研究所 | 一种基于改进foa-svm的生产物流预警信息多分类方法 |
CN111898249A (zh) * | 2020-07-02 | 2020-11-06 | 华能澜沧江水电股份有限公司 | 滑坡位移非参数概率密度预测方法、设备及存储介质 |
CN112132308A (zh) * | 2020-07-27 | 2020-12-25 | 国网湖南省电力有限公司 | 电网暴雨长期预测关键因子选取方法、装置及预测系统 |
CN112381115A (zh) * | 2020-10-21 | 2021-02-19 | 西安工程大学 | 一种Bagging滑坡预报方法 |
CN112507780A (zh) * | 2020-10-21 | 2021-03-16 | 西安工程大学 | 一种基于tcn深度学习的滑坡敏感性评价方法 |
CN112597647A (zh) * | 2020-12-17 | 2021-04-02 | 河海大学 | 一种快速收敛的超高频微波降雨数据离散化方法 |
CN113723541A (zh) * | 2021-09-02 | 2021-11-30 | 国网福建省电力有限公司电力科学研究院 | 一种基于混合智能算法的边坡位移预测方法 |
CN115019478A (zh) * | 2022-08-09 | 2022-09-06 | 中地宝联(北京)国土资源勘查技术开发集团有限公司 | 一种应用于临界滑坡的监测预警方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010266318A (ja) * | 2009-05-14 | 2010-11-25 | Shizuoka Prefecture | 湿球温度とwbgtの予測方法、wbgt計、および熱中症危険度判定装置 |
CN104091216A (zh) * | 2014-07-29 | 2014-10-08 | 吉林大学 | 基于果蝇优化最小二乘支持向量机的交通信息预测方法 |
CN108154651A (zh) * | 2018-01-29 | 2018-06-12 | 西安工程大学 | 一种用于滑坡灾害实时监测预警的装置和方法 |
-
2018
- 2018-11-01 CN CN201811294315.6A patent/CN109272721B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010266318A (ja) * | 2009-05-14 | 2010-11-25 | Shizuoka Prefecture | 湿球温度とwbgtの予測方法、wbgt計、および熱中症危険度判定装置 |
CN104091216A (zh) * | 2014-07-29 | 2014-10-08 | 吉林大学 | 基于果蝇优化最小二乘支持向量机的交通信息预测方法 |
CN108154651A (zh) * | 2018-01-29 | 2018-06-12 | 西安工程大学 | 一种用于滑坡灾害实时监测预警的装置和方法 |
Non-Patent Citations (2)
Title |
---|
彭令等: "基于核主成成分分析和粒子群优化支持向量机的滑坡位移预测", 《武汉大学学报•信息科学版》 * |
谢洋洋等: "基于相空间重构的FOA-GLSSVM深基坑变形预测模型研究", 《大地测量与地球动力学》 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109887240A (zh) * | 2019-03-22 | 2019-06-14 | 福州大学 | 一种基于人工智能的滑坡灾害安全监测与预警方法 |
CN110196448B (zh) * | 2019-05-21 | 2020-06-16 | 成都信息工程大学 | 一种滑坡次声信号识别方法 |
CN110196448A (zh) * | 2019-05-21 | 2019-09-03 | 成都信息工程大学 | 一种滑坡次声信号识别方法 |
CN110456217B (zh) * | 2019-08-08 | 2021-06-01 | 辽宁工程技术大学 | 一种基于wpd-foa-lssvm双模型的mmc故障定位方法 |
CN110456217A (zh) * | 2019-08-08 | 2019-11-15 | 辽宁工程技术大学 | 一种基于wpd-foa-lssvm双模型的mmc故障定位方法 |
CN110598914B (zh) * | 2019-08-21 | 2021-12-10 | 武汉理工大学 | 一种多因素影响下矿井灾害气体浓度区间预测方法及系统 |
CN110598914A (zh) * | 2019-08-21 | 2019-12-20 | 武汉理工大学 | 一种多因素影响下矿井灾害气体浓度区间预测方法及系统 |
CN111783866A (zh) * | 2020-06-28 | 2020-10-16 | 中国电子科技集团公司第二十研究所 | 一种基于改进foa-svm的生产物流预警信息多分类方法 |
CN111898249A (zh) * | 2020-07-02 | 2020-11-06 | 华能澜沧江水电股份有限公司 | 滑坡位移非参数概率密度预测方法、设备及存储介质 |
CN112132308A (zh) * | 2020-07-27 | 2020-12-25 | 国网湖南省电力有限公司 | 电网暴雨长期预测关键因子选取方法、装置及预测系统 |
CN112132308B (zh) * | 2020-07-27 | 2023-07-25 | 国网湖南省电力有限公司 | 电网暴雨长期预测关键因子选取方法、装置及预测系统 |
CN112381115A (zh) * | 2020-10-21 | 2021-02-19 | 西安工程大学 | 一种Bagging滑坡预报方法 |
CN112507780A (zh) * | 2020-10-21 | 2021-03-16 | 西安工程大学 | 一种基于tcn深度学习的滑坡敏感性评价方法 |
CN112597647A (zh) * | 2020-12-17 | 2021-04-02 | 河海大学 | 一种快速收敛的超高频微波降雨数据离散化方法 |
CN112597647B (zh) * | 2020-12-17 | 2021-06-01 | 河海大学 | 一种快速收敛的超高频微波降雨数据离散化方法 |
CN113723541A (zh) * | 2021-09-02 | 2021-11-30 | 国网福建省电力有限公司电力科学研究院 | 一种基于混合智能算法的边坡位移预测方法 |
CN113723541B (zh) * | 2021-09-02 | 2024-02-09 | 国网福建省电力有限公司电力科学研究院 | 一种基于混合智能算法的边坡位移预测方法 |
CN115019478A (zh) * | 2022-08-09 | 2022-09-06 | 中地宝联(北京)国土资源勘查技术开发集团有限公司 | 一种应用于临界滑坡的监测预警方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109272721B (zh) | 2020-10-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109272721A (zh) | 基于kpca-foa-lssvm的滑坡灾害预报方法 | |
WO2023029461A1 (zh) | 海量高维ais轨迹数据聚类方法 | |
CN101354757B (zh) | 一种精细尺度下的动态风险及易损性预测方法 | |
CN209168369U (zh) | 一种滑坡灾害实时监测装置 | |
CN106778841A (zh) | 异常用电检测模型的建立方法 | |
CN109886456B (zh) | 基于pca与混合核函数lssvr的泥石流预测方法 | |
CN112966722A (zh) | 基于半监督随机森林模型的区域滑坡易发性预测方法 | |
Ismail et al. | A hybrid model of self organizing maps and least square support vector machine for river flow forecasting | |
CN105241524A (zh) | 一种基于径向基函数神经网络模型的城市洪水预警系统及其方法 | |
CN104240507A (zh) | 一种基于多视角融合的交通小区划分方法 | |
Iqbal et al. | Groundwater level prediction model using correlation and difference mechanisms based on boreholes data for sustainable hydraulic resource management | |
CN115270632A (zh) | 饮用水水源地水华风险预测系统和方法 | |
CN115482138A (zh) | 基于特征筛选与差分进化算法优化的滑坡危险性评估方法 | |
CN107194605A (zh) | 能源消耗估计系统及方法 | |
CN117238126A (zh) | 一种连续流道路场景下的交通事故风险评估方法 | |
CN115600513A (zh) | 一种岩溶塌陷监测预警及防控一体化信息化模拟研判系统 | |
CN117788218B (zh) | 一种碳排放评估方法及系统 | |
Momeni et al. | Pattern‐based calibration of cellular automata by genetic algorithm and Shannon relative entropy | |
CN117933701A (zh) | 轨道交通工程建设安全风险监控方法及系统 | |
Zanella et al. | Internet of things for hydrology: Potential and challenges | |
CN117808214A (zh) | 水利工程数据分析系统 | |
CN117370785A (zh) | 智慧城市的环境监测方法、终端及环境监测系统 | |
CN209433536U (zh) | 一种泥石流监测预警装置 | |
CN115829150A (zh) | 一种积水预测系统 | |
CN116307755A (zh) | 基于gis的国土空间规划管控层级预警系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |