买专利,只认龙图腾
首页 专利交易 科技果 科技人才 科技服务 商标交易 会员权益 IP管家助手 需求市场 关于龙图腾
 /  免费注册
到顶部 到底部
清空 搜索

【发明授权】基于改进贝叶斯网络的滑坡灾害预测方法_西安工程大学_201810836499.8 

申请/专利权人:西安工程大学

申请日:2018-07-26

公开(公告)日:2021-11-19

公开(公告)号:CN109146144B

主分类号:G06Q10/04(20120101)

分类号:G06Q10/04(20120101);G06K9/62(20060101)

优先权:

专利状态码:有效-授权

法律状态:2021.11.19#授权;2019.01.29#实质审查的生效;2019.01.04#公开

摘要:本发明公开的基于改进贝叶斯网络的滑坡灾害预测方法,首先建立滑坡体监测预警系统,获取初始滑坡灾害影响因子数据,经标准化处理后运用PCA算法提取并筛选出主要影响因子;将筛选出的主要影响因子样本数据按比例分为训练样本集与测试样本集,并划分出滑坡发生等级;然后构建基于改进贝叶斯的滑坡预测模型,通过加权处理计算出引入特征量后的滑坡发生概率,并添加闭环学习过程,完成滑坡灾害预测;本发明公开的方法通过筛选滑坡灾害主要影响因子,在基于改进贝叶斯模型中对特征量加权处理并添加闭环学习环节输出滑坡发生概率,增加容错率,提升预报精度。

主权项:1.基于改进贝叶斯网络的滑坡灾害预测方法,其特征在于,具体的操作步骤如下:步骤1.建立滑坡体监测预警系统,获取初始滑坡灾害影响因子数据,经标准化处理后运用主成分分析法提取并筛选出主要影响因子;所述步骤1中经标准化处理后运用主成分分析法提取并筛选出主要影响因子过程如下:步骤1.1将数据标准化处理,如公式1所示, 式中,Xnor为经标准化的滑坡灾害发生影响因子数据,XOriginal为原始数据,Xmin,Xmax分别为原始数据最小值和最大值;步骤1.2主成分分析法提取并筛选主要灾害影响因子的过程如下:设式1中的Xnor含M个变量、N个样本的数据集,第k个主成分向量如公式2所示:Pk=p1k,p2k,…pMkT2其中,k=1,2,…,N;Pk为第k个主成分向量;在第k个主成分情况下,每个数据分别表示为P1K,P2K…PMk;Xnor在主成分向量矩阵P的投影为样本的主成分矩阵T=XnorP,主成分的矩阵元素如公式3所示: 其中,j=1,2,…,N,Tj为主成分矩阵元素,Pk为主成分向量,Xnor为经标准化的数据集;主成分矩阵相互正交如公式4所示:TjTTj=XnorPkTXnorPk=PkTXnorTXnorPk=PkTDPk=Λ4式中,Tj为主成分矩阵,TjT主成分矩阵的转置,Xnor为经标准化的数据集,D为协方差矩阵,Λ为对角矩阵;根据μj贡献率选取主要成分Tj,如公式5所示: 其中,j为主元成分个数,μj为累计方差贡献率,λj为第j个特征值,为总特征值;在公式5中取累计方差贡献率μj≥85%作为主要滑坡灾害影响因子数据,运用Python所具有的矩阵计算功能优势实现对滑坡灾害诱发数据的主成分分析,经过主成分变换,得到主要影响因子的矩阵特征值与特征向量,忽略无关且错误的信息;步骤2.将筛选出的主要影响因子样本数据按比例分为训练样本集与测试样本集,并划分出滑坡灾害发生等级;所述步骤2的具体过程如下:步骤2.1令步骤1筛选出的主要灾害影响因子特征量分别为:X1=降雨量,X2=土壤含水率,X3=裂缝位移量,X4=滑坡坡向,X5=坡型,X6=岩性;步骤2.2选取滑坡影响因子共K个样本,按7:3的比例分为训练样本和测试样本;步骤2.3将滑坡可能发生概率等级分为五个等级,分别为1,2,3,4,5个等级;步骤3.将步骤2得到的滑坡灾害发生等级构建基于改进贝叶斯网络的滑坡预测模型,通过加权处理计算出引入特征量后的滑坡发生概率,并添加闭环学习过程,完成滑坡灾害预测;所述步骤3的具体过程如下:步骤3.1构建基于改进贝叶斯的滑坡预测模型,首先引入相关特征量:令XP=X1,X2…,XK为某滑坡发生概率等级所对应影响因子所对应的属性的特征向量,X1,X2,…XK为特征分量;令M={Mi|i=1,2,3,4,5}为已知滑坡发生概率等级;根据筛选出的主要灾害影响因子,构建预测指标体系;步骤3.2针对样本集中Mi类滑坡发生等级及其对应特征,计算Mi类滑坡发生等级次数及相应灾害因子Xj出现次数,同时计算非Mi类滑坡发生等级次数及相应灾害因子Xj出现次数;根据贝叶斯概率公式,如公式6所示: 其中,PMiX为在Mi滑坡发生等级时,该等级的特征向量为X的概率;PMi为先验概率,即,Mi滑坡等级的发生概率;PXP为某滑坡发生等级的特征向量等于X发生的概率;公式6中PXP为: 其中,K为样本总数,KM为第M类等级的样本数;同时,公式6中PXPMi为: 其中,在特征分量相互独立的条件下,PXjMi为滑坡灾害发生等级为Mi的特征分量Xj出现的概率;Xj为X的第j个特征分量值;由于滑坡灾害影响要素是互为关联的,因此将分量之间增加权重,并按照每个成灾因子的所占比重进行排序;对公式8进行加权处理,改进后的贝叶斯算法为: 其中,为加权之后,滑坡发生等级为Mi的特征分量Xj出现的概率;步骤3.3计算灾害因子在不同滑坡灾害发生等级下的后验概率:当滑坡灾害等级为Mi时,先验概率GMi为: 最大似然比表示为: 滑坡等级为Mi时的后验概率由公式10和公式11得到:GMi|Xj=LXj|Mi·GMi12选取GMi|Xj中最大的Mn作为目前最可能出现的滑坡等级类别,则表示为: 将训练样本代入滑坡灾害模型中进行运算,最终由公式13中计算出的滑坡等级类别,与专家知识和存储在地区历史数据库中的知识进行对比,若一致或相差较小,则完成滑坡预测;否则进入闭环学习阶段,将错误结果进行修正后再返回滑坡模型中并重新计算滑坡可能发生概率。

全文数据:基于改进贝叶斯网络的滑坡灾害预测方法技术领域本发明属于地质灾害预报方法技术领域,涉及一种基于改进贝叶斯网络的滑坡灾害预测方法。背景技术滑坡是一种重要的地质灾害类型,对地区交通运输、电站、厂矿、城乡建设等造成很大程度的影响,并且严重威胁着人民的生命和财产安全。滑坡的发生受多重因素的影响,如何通过滑坡的影响因素来提前预测其发生的可能性成为人们所关注的焦点。针对滑坡灾害发生频率高、影响范围广等特点,近年来国内学者基于不同方法对其进行预测预报,但仍然存在一定的局限性。如神经网络模型,具有良好的自学习和容错能力,但易于陷入局部最优且收敛速度较慢;模糊综合评价模型,能对蕴藏信息呈现模糊性的资料作出比较切合实际的量化评价,但计算过程复杂且对指标权重矢量的确定主观性较强;SVM支持向量机模型,具有较好的“鲁棒性”,但对大规模训练样本难以实施;传统贝叶斯判别模型,具有稳定的分类效率、能够处理多分类任务,是处理不确定问题的有效工具,但在判断较为复杂的样本时容易出现误判。而改进的贝叶斯算法,通过引入加权处理与闭环学习系统,可对特征量进行权重排序,并能够修正错误信息,达到提升容错率的效果。鉴于此,本发明提出一种基于改进贝叶斯网络ImprovedBayesnet,IBN的滑坡灾害预测方法,首先运用主成分分析法principalcomponentsanalysis,PCA对滑坡特征的样本进行处理并筛选主要灾害影响因子,将主要灾害影响因子按比例构成训练样本集和测试样本集并划分滑坡发生等级,然后基于改进贝叶斯的预测模型中,建立预测评估体系,并计算出引入加权后的特征量存在滑坡发生概率,添加闭环学习过程,完成滑坡灾害预测。发明内容本发明目的是提供一种基于改进贝叶斯网络的滑坡灾害预测方法,解决现有灾害预测中所采用的算法容错率低,精确度不高的问题。本发明所采用的技术方案是,基于改进贝叶斯网络的滑坡灾害预测方法,具体的操作步骤如下:步骤1.建立滑坡体监测预警系统,获取初始滑坡灾害影响因子数据,经标准化处理后运用主成分分析法提取并筛选出主要影响因子;步骤2.将筛选出的主要影响因子样本数据按比例分为训练样本集与测试样本集,并划分出滑坡灾害发生等级;步骤3.将步骤2得到的滑坡灾害发生等级构建基于改进贝叶斯网络的滑坡预测模型,通过加权处理计算出引入特征量后的滑坡发生概率,并添加闭环学习过程,完成滑坡灾害预测。本发明的其他特点还在于,步骤1中的滑坡体监测预报系统包括现场预警终端和若干个无线智能传感器模块,现场预警终端与PC端控制中心信号连接,现场预警终端与若干个无线智能传感器模块通过Zigbee无线传输方式通讯;现场预警终端包括主控芯片STM32,主控芯片STM32分别连接人机友好界面、USB接口、储存模块、GPRS模块和Zigbee模块;无线智能传感器模块包括主控芯片MSP430,主控芯片MSP430分别连接有4-20mA模块、RS232模块、RS485模块、Zigbee无线通讯模块。步骤1中经标准化处理后运用主成分分析法提取并筛选出主要影响因子过程如下:步骤1.1将数据标准化处理,如公式1所示,式中,X为经标准化的滑坡灾害发生影响因子数据,XOriginal为原始数据,Xmin,Xmax分别为原始数据最小值和最大值;步骤1.2主成分分析法提取并筛选主要灾害影响因子的过程如下:设X是滑坡影响因子数据经标准化后含M个变量、N个样本的数据集,第k个主成分向量如公式2所示:Pk=p1k,p2k,…pMkT2其中,k=1,2,…,N;Pk为第k个主成分向量;在第k个主成分情况下,每个数据分别表示为P1K,P2K…PMk;X在主成分向量矩阵P的投影为样本的主成分矩阵T=XP,主成分的矩阵元素如公式3所示:其中,j=1,2,…,N,Tj为主成分矩阵元素,Pk为主成分向量,X为经标准化的数据集。主成分矩阵相互正交如公式4所示:TjΤTj=XPkΤXPk=PkΤXΤXPk=PkΤDPk=Λ4式中,Tj为主成分矩阵,TjT主成分矩阵的转置,X为经标准化的数据集,D为协方差矩阵,Λ为对角矩阵;根据μj贡献率选取主要成分Tj,如公式5所示:其中,j为主元成分个数,μj为累计方差贡献率,λj为第j个特征值,为总特征值;在公式5中取累计方差贡献率μj≥85%作为主要滑坡灾害影响因子数据,运用Python所具有的矩阵计算功能优势实现对滑坡灾害诱发数据的主成分分析,经过主成分变换,得到主要影响因子的矩阵特征值与特征向量,忽略无关且错误的信息。步骤2的具体过程如下:步骤2.1令步骤1筛选出的主要灾害影响因子特征量分别为:X1=降雨量,X2=土壤含水率,X3=裂缝位移量,X4=滑坡坡向,X5=坡型,X6=岩性;步骤2.2选取滑坡影响因子共K个样本,按7:3的比例分为训练样本和测试样本;步骤2.3将滑坡可能发生概率等级分为五个等级,分别为1,2,3,4,5个等级。步骤3的具体过程如下:步骤3.1构建基于改进贝叶斯的滑坡预测模型,首先引入相关特征量:令X=X1,X2…,XK为某滑坡发生概率等级所对应影响因子所对应的属性的特征向量,X1,X2,…XL为特征分量;令M={Mi|i=1,2,3,4,5}为已知滑坡发生概率等级;根据筛选出的主要灾害影响因子,构建预测指标体系;步骤3.2针对样本集中Mi类滑坡发生等级及其对应特征,计算Mi类滑坡发生等级次数及相应灾害因子Xj出现次数,同时计算非Mi类滑坡发生等级次数及相应灾害因子Xj出现次数;根据贝叶斯概率公式,如公式6所示:其中,PMiX为在Mi滑坡发生等级时,该等级的特征向量为X的概率;PMi为先验概率,即,Mi滑坡等级的发生概率;PX为某滑坡发生等级的特征向量等于X发生的概率;公式6中PX为:其中,K为样本总数,KM为第M类等级的样本数;同时,公式6中PXMi为:其中,在特征分量相互独立的条件下,PXjMi为滑坡灾害发生等级为Mi的特征分量Xj出现的概率;Xj为X的第j个特征分量值;由于滑坡灾害影响要素是互为关联的,因此将分量之间增加权重,并按照每个成灾因子的所占比重进行排序;对公式8进行加权处理,改进后的贝叶斯算法为:其中,为加权之后,滑坡发生等级为Mi的特征分量Xj出现的概率;步骤3.3计算灾害因子在不同滑坡灾害发生等级下的后验概率:当滑坡灾害等级为Mi时,先验概率GMi为:最大似然比表示为:滑坡等级为Mi时的后验概率由公式10和公式11得到:GMi|Xj=LXj|Mi·GMi12选取GMi|Xj中最大的Mn作为目前最可能出现的滑坡等级类别,则表示为:将步骤2中得到的训练样本代入滑坡灾害模型中进行运算,最终由公式13中计算出的滑坡等级类别,与专家知识和存储在地区历史数据库中的知识进行对比,若一致或相差较小,则完成滑坡预测;否则进入闭环学习阶段,将错误结果进行修正后再返回滑坡模型中并重新计算滑坡可能发生概率。本发明的有益效果是,基于改进贝叶斯网络的滑坡灾害预测方法,解决了现有灾害预测中所采用的算法容错率低,预测精确度不高的问题;本发明的方法具有以下的有益效果:1建立滑坡监测预警系统,通过智能传感器采集到大量灾害影响因子,并通过PCA算法筛选出主要灾害影响因子,减少指标选择工作量;2通过IBN对灾害影响因子样本数据进行特征学习,通过特征量加权处理改进独立性假设,解决影响因子之间互为关联问题,进一步提升预报精确度;3添加闭环学习环节,对错误数据纠错,提升网络容错率。附图说明图1是本发明的基于改进贝叶斯网络的滑坡灾害预测方法的流程图;图2是本发明的基于改进贝叶斯网络的滑坡灾害预测方法所使用的滑坡体监测预警系统结构示意图;图3是改进贝叶斯网络与预测评估关系图;图4是闭环学习环节图。具体实施方式下面结合附图和具体实施方式对本发明进行详细说明。本发明的基于改进贝叶斯网络的滑坡灾害预测方法,如图1所示,具体的操作步骤如下:步骤1.建立滑坡体监测预警系统,获取初始滑坡灾害影响因子数据,经标准化处理后运用主成分分析法PCA提取并筛选出主要影响因子;步骤2.将筛选出的主要影响因子样本数据按比例分为训练样本集与测试样本集,并划分出滑坡灾害发生等级;步骤3.将步骤2得到的滑坡灾害发生等级,建立基于改进贝叶斯的滑坡预测模型,通过加权处理计算出引入特征量后的滑坡发生概率,并添加闭环学习过程,完成滑坡灾害预测;步骤1中的滑坡体监测预报系统,如图2所示,包括现场预警终端和若干个无线智能传感器模块,现场预警终端与PC端控制中心信号连接,现场预警终端与若干个无线智能传感器模块通过Zigbee无线传输方式通讯;现场预警终端包括主控芯片STM32,主控芯片STM32分别连接人机友好界面、USB接口、储存模块、GPRS模块和Zigbee模块;无线智能传感器模块包括主控芯片MSP430,主控芯片MSP430分别连接有4-20mA模块、RS232模块、RS485模块、Zigbee无线通讯模块。步骤1中具体过程如下:步骤1.1建立滑坡体监测预警系统,获取初始滑坡灾害影响因子数据;步骤1.2将数据标准化处理,如公式1所示,式中,X为经标准化的滑坡灾害发生影响因子数据,XOriginal为原始数据,Xmin,Xmax分别为原始数据最小值和最大值;步骤1.3PCA算法提取并筛选主要灾害影响因子的过程如下,设X是滑坡影响因子数据经标准化后含M个变量、N个样本的数据集,第k个主成分向量如公式2所示:Pk=p1k,p2k,…pMkT2其中,k=1,2,…,N;Pk为第k个主成分向量;在第k个主成分情况下,每个数据分别表示为P1K,P2K…PMk;X在主成分向量矩阵P的投影为样本的主成分矩阵T=XP,本文主成分的矩阵元素如公式3所示:其中,j=1,2,…,N,Tj为主成分矩阵元素,Pk为主成分向量,X为经标准化的数据集;主成分矩阵相互正交如公式4所示:TjΤTj=XPkΤXPk=PkΤXΤXPk=PkΤDPk=Λ4式中,Tj为主成分矩阵,TjT主成分矩阵的转置,X为经标准化的数据集,D为协方差矩阵,Λ为对角矩阵;根据μj贡献率选取主要成分Tj,如公式5所示:其中,j为主元成分个数,μj为累计方差贡献率,λj为第j个特征值,为总特征值;在公式5中取累计方差贡献率μj≥85%作为主要滑坡灾害影响因子数据,运用Python所具有的矩阵计算功能优势实现对滑坡灾害诱发数据的主成分分析,经过主成分变换,得到主要影响因子的矩阵特征值与特征向量,忽略掉无关且错误的信息。步骤2中具体过程如下:步骤2.1令步骤1筛选出的主要灾害影响因子特征量分别为:X1=降雨量,X2=土壤含水率,X3=裂缝位移量,X4=滑坡坡向,X5=坡型,X6=岩性;步骤2.2选取滑坡影响因子共K个样本,按7:3的比例分为训练样本和测试样本;步骤2.3将滑坡可能发生概率等级分为五个等级,分别为1,2,3,4,5个等级,分别对应为发生可能性很小,发生可能性小,发生可能性中度,发生可能性大,发生可能性很大,具体如表1所示,表1步骤3的具体过程如下,步骤3.1构建基于改进贝叶斯的滑坡预测模型,首先引入相关特征量;令X=X1,X2…,XK为某滑坡发生概率等级所对应影响因子所对应的属性的特征向量,X1,X2,…XL为特征分量;令M={Mi|i=1,2,3,4,5}为已知滑坡发生概率等级;根据筛选出的主要灾害影响因子,构建预测评估指标体系,如表2所示;且改进贝叶斯网络与预测评估关系图如图3所示,从两者的实现流程可以看出,改进贝叶斯网络建模与滑坡预测流程是环环对应的。由开始的构建预测指标体系对应IBN结构学习;再到计算指标权重对应IBN参数学习;然后建立指标融合的模型对应IBN推理计算;最后预测评估结果对应滑坡发生的概率值。表2步骤3.2针对样本集中Mi类滑坡发生等级及其对应特征,计算Mi类滑坡发生等级次数及相应灾害因子Xj出现次数,同时计算非Mi类滑坡发生等级次数及相应灾害因子Xj出现次数;根据贝叶斯概率公式,如公式6所示:其中,PMiX为在Mi滑坡发生等级时,该等级的特征向量为X的概率;PMi为先验概率,即,Mi滑坡等级的发生概率;PX为某滑坡发生等级的特征向量等于X发生的概率;公式6中PX为:其中,K为样本总数,KM为第M类等级的样本数;同时,公式6中PXMi为:其中,在特征分量相互独立的条件下,PXjMi为滑坡灾害发生等级为Mi的特征分量Xj出现的概率;Xj为X的第j个特征分量值;当在某一滑坡发生等级中,主要灾害影响因子可能呈现出不同的现象,不同滑坡发生等级其特征向量也不完全相同,因此向量的特征分量也要最大程度地反映该滑坡发生可能性的特征。但传统贝叶斯算法运算时的特征分量是相互独立的,由于滑坡灾害影响要素是互为关联的,因此将分量之间增加权重,并按照每个成灾因子的所占比重进行排序。对公式8进行加权处理,改进后的贝叶斯算法为:其中,为加权之后,滑坡发生等级为Mi,特征分量Xj出现的概率;步骤3.3计算灾害因子在不同滑坡发生等级下的后验概率;当滑坡等级为Mi时,先验概率GMi可以表示为:最大似然比可以表示为:滑坡等级为Mi时的后验概率可由公式10和公式11得到:GMi|Xj=LXj|Mi·GMi12选取GMi|Xj中最大的Mn作为目前最可能出现的滑坡等级类别,如公式13所示:将步骤2得到的训练样本代入滑坡灾害模型中进行运算,最终由公式13中计算出的滑坡等级类别,与专家知识和存储在地区历史数据库中的知识进行对比,若一致或相差较小,则完成滑坡预测;否则进入闭环学习阶段,其闭环学习系统图如图4所示,将待分类的数据输入到分类系统中,得到分类结果,然后再次将分类结果经过学习后输入到分类系统,即就是其将错误结果进行修正后再返回滑坡模型中并重新计算滑坡可能发生概率。传统贝叶斯网络经过一次学习分类后将不再从实际训练过程中获取新的信息,缺乏实时学习能力。然而,本发明改进后的贝叶斯网络将对传统贝叶斯模型进行修正,使用了闭环学习系统,使其性能更为优化。

权利要求:1.基于改进贝叶斯网络的滑坡灾害预测方法,其特征在于,具体的操作步骤如下:步骤1.建立滑坡体监测预警系统,获取初始滑坡灾害影响因子数据,经标准化处理后运用主成分分析法提取并筛选出主要影响因子;步骤2.将筛选出的主要影响因子样本数据按比例分为训练样本集与测试样本集,并划分出滑坡灾害发生等级;步骤3.将步骤2得到的滑坡灾害发生等级构建基于改进贝叶斯网络的滑坡预测模型,通过加权处理计算出引入特征量后的滑坡发生概率,并添加闭环学习过程,完成滑坡灾害预测。2.如权利要求1所述的基于改进贝叶斯网络的滑坡灾害预测方法,其特征在于,所述步骤1中的滑坡体监测预报系统包括现场预警终端和若干个无线智能传感器模块,所述现场预警终端与PC端控制中心信号连接,所述现场预警终端与若干个无线智能传感器模块通过Zigbee无线传输方式通讯;所述现场预警终端包括主控芯片STM32,所述主控芯片STM32分别连接人机友好界面、USB接口、储存模块、GPRS模块和Zigbee模块;所述无线智能传感器模块包括主控芯片MSP430,所述主控芯片MSP430分别连接有4-20mA模块、RS232模块、RS485模块、Zigbee无线通讯模块。3.如权利要求1所述的基于改进贝叶斯网络的滑坡灾害预测方法,其特征在于,所述步骤1中经标准化处理后运用主成分分析法提取并筛选出主要影响因子过程如下:步骤1.1将数据标准化处理,如公式1所示,式中,X为经标准化的滑坡灾害发生影响因子数据,XOriginal为原始数据,Xmin,Xmax分别为原始数据最小值和最大值;步骤1.2主成分分析法提取并筛选主要灾害影响因子的过程如下:设X是滑坡影响因子数据经标准化后含M个变量、N个样本的数据集,第k个主成分向量如公式2所示:Pk=p1k,p2k,…pMkT2其中,k=1,2,…,N;Pk为第k个主成分向量;在第k个主成分情况下,每个数据分别表示为P1K,P2K…PMk;X在主成分向量矩阵P的投影为样本的主成分矩阵T=XP,主成分的矩阵元素如公式3所示:其中,j=1,2,…,N,Tj为主成分矩阵元素,Pk为主成分向量,X为经标准化的数据集。主成分矩阵相互正交如公式4所示:TjΤTj=XPkΤXPk=PkΤXΤXPk=PkΤDPk=Λ4式中,Tj为主成分矩阵,TjT主成分矩阵的转置,X为经标准化的数据集,D为协方差矩阵,Λ为对角矩阵;根据μj贡献率选取主要成分Tj,如公式5所示:其中,j为主元成分个数,μj为累计方差贡献率,λj为第j个特征值,为总特征值;在公式5中取累计方差贡献率μj≥85%作为主要滑坡灾害影响因子数据,运用Python所具有的矩阵计算功能优势实现对滑坡灾害诱发数据的主成分分析,经过主成分变换,得到主要影响因子的矩阵特征值与特征向量,忽略无关且错误的信息。4.如权利要求1所述的基于改进贝叶斯网络的滑坡灾害预测方法,其特征在于,所述步骤2的具体过程如下:步骤2.1令步骤1筛选出的主要灾害影响因子特征量分别为:X1=降雨量,X2=土壤含水率,X3=裂缝位移量,X4=滑坡坡向,X5=坡型,X6=岩性;步骤2.2选取滑坡影响因子共K个样本,按7:3的比例分为训练样本和测试样本;步骤2.3将滑坡可能发生概率等级分为五个等级,分别为1,2,3,4,5个等级。5.如权利要求4所述的基于改进贝叶斯网络的滑坡灾害预测方法,其特征在于,所述步骤3的具体过程如下:步骤3.1构建基于改进贝叶斯的滑坡预测模型,首先引入相关特征量:令X=X1,X2…,XK为某滑坡发生概率等级所对应影响因子所对应的属性的特征向量,X1,X2,…XL为特征分量;令M={Mi|i=1,2,3,4,5}为已知滑坡发生概率等级;根据筛选出的主要灾害影响因子,构建预测指标体系;步骤3.2针对样本集中Mi类滑坡发生等级及其对应特征,计算Mi类滑坡发生等级次数及相应灾害因子Xj出现次数,同时计算非Mi类滑坡发生等级次数及相应灾害因子Xj出现次数;根据贝叶斯概率公式,如公式6所示:其中,PMiX为在Mi滑坡发生等级时,该等级的特征向量为X的概率;PMi为先验概率,即,Mi滑坡等级的发生概率;PX为某滑坡发生等级的特征向量等于X发生的概率;公式6中PX为:其中,K为样本总数,KM为第M类等级的样本数;同时,公式6中PXMi为:其中,在特征分量相互独立的条件下,PXjMi为滑坡灾害发生等级为Mi的特征分量Xj出现的概率;Xj为X的第j个特征分量值;由于滑坡灾害影响要素是互为关联的,因此将分量之间增加权重,并按照每个成灾因子的所占比重进行排序;对公式8进行加权处理,改进后的贝叶斯算法为:其中,为加权之后,滑坡发生等级为Mi的特征分量Xj出现的概率;步骤3.3计算灾害因子在不同滑坡灾害发生等级下的后验概率:当滑坡灾害等级为Mi时,先验概率GMi为:最大似然比表示为:滑坡等级为Mi时的后验概率由公式10和公式11得到:GMi|Xj=LXj|Mi·GMi12选取GMi|Xj中最大的Mn作为目前最可能出现的滑坡等级类别,则表示为:将训练样本代入滑坡灾害模型中进行运算,最终由公式13中计算出的滑坡等级类别,与专家知识和存储在地区历史数据库中的知识进行对比,若一致或相差较小,则完成滑坡预测;否则进入闭环学习阶段,将错误结果进行修正后再返回滑坡模型中并重新计算滑坡可能发生概率。

百度查询: 西安工程大学 基于改进贝叶斯网络的滑坡灾害预测方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息,力求客观、公正,但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解,仅供参考使用,不能作为本公司承担任何法律责任的依据或者凭证。