CN109063247A - 基于深度信念网络的滑坡灾害预报方法 - Google Patents
基于深度信念网络的滑坡灾害预报方法 Download PDFInfo
- Publication number
- CN109063247A CN109063247A CN201810667868.5A CN201810667868A CN109063247A CN 109063247 A CN109063247 A CN 109063247A CN 201810667868 A CN201810667868 A CN 201810667868A CN 109063247 A CN109063247 A CN 109063247A
- Authority
- CN
- China
- Prior art keywords
- layer
- landslide
- network
- disaster
- sample
- 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
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000012549 training Methods 0.000 claims abstract description 40
- 238000012544 monitoring process Methods 0.000 claims abstract description 13
- 238000012360 testing method Methods 0.000 claims abstract description 13
- 230000002068 genetic effect Effects 0.000 claims abstract description 7
- 230000000007 visual effect Effects 0.000 claims description 32
- 230000001939 inductive effect Effects 0.000 claims description 30
- 230000006870 function Effects 0.000 claims description 24
- 210000004027 cell Anatomy 0.000 claims description 21
- 230000004913 activation Effects 0.000 claims description 15
- 238000013528 artificial neural network Methods 0.000 claims description 15
- 238000013277 forecasting method Methods 0.000 claims description 14
- 238000012216 screening Methods 0.000 claims description 11
- 230000035772 mutation Effects 0.000 claims description 6
- 210000000349 chromosome Anatomy 0.000 claims description 5
- 238000004891 communication Methods 0.000 claims description 4
- 238000006073 displacement reaction Methods 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 239000011148 porous material Substances 0.000 claims description 3
- 239000002689 soil Substances 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 239000010410 layer Substances 0.000 claims 31
- 239000011229 interlayer Substances 0.000 claims 1
- 238000005457 optimization Methods 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 101000742844 Homo sapiens RNA-binding motif protein, Y chromosome, family 1 member A1 Proteins 0.000 description 2
- 102100038040 RNA-binding motif protein, Y chromosome, family 1 member A1 Human genes 0.000 description 2
- 230000006378 damage Effects 0.000 description 2
- 238000013135 deep learning Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 241000345998 Calamus manan Species 0.000 description 1
- 235000002568 Capsicum frutescens Nutrition 0.000 description 1
- 101001046999 Homo sapiens Kynurenine-oxoglutarate transaminase 3 Proteins 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000001537 neural effect Effects 0.000 description 1
- 239000002243 precursor Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 235000012950 rattan cane Nutrition 0.000 description 1
- 230000002087 whitening effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/086—Learning methods using evolutionary algorithms, e.g. genetic algorithms or genetic programming
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Computing Systems (AREA)
- Software Systems (AREA)
- Mathematical Physics (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Computer Hardware Design (AREA)
- Physiology (AREA)
- Evolutionary Biology (AREA)
- Geometry (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
本发明公开的基于深度信念网络的滑坡灾害预报方法,首先建立滑坡体监测预警系统,采集大量灾害诱发因子,运用MIV算法筛选,筛选出主要灾害诱发因子;将筛选出的灾害诱发因子进行标准化过程,按特定比例分为测试样本和训练样本;然后,构建基于深度信念网络的滑坡灾害预报模型,设置结构为两层RBM和三层BP网络组成;采用CD算法对RBM预训练,更新网络参数;采用遗传算法监督训练学习,保证DBN网络整体最优;最后,重构已优化的滑坡灾害预报模型,对输出结果进行滑坡等级划分,预报出滑坡发生可能性。本发明公开的方法通过提取特征灾害诱发因子,加快收敛速度,防止陷入局部最优,提升滑坡灾害预报精度。
Description
技术领域
本发明属于地质灾害预报方法技术领域,涉及一种基于深度信念网络的滑坡灾害预报方法。
背景技术
滑坡是重要的地质灾害类型之一,它威胁着人类的生命及财产安全,并对灾害发生地区的基础设施和生态环境具有极大破坏作用。因此,如何运用相应技术手段对滑坡灾害进行实时性监测预报,最大程度降低损失,成为我们关注的主要内容。
现有的滑坡灾害预报方法有很多,研究阶段分为多个时期。第一阶段处于20世纪60~70年代,以现象预报和经验预报为主,即著名的“斋藤法”。专家凭借滑坡破坏现象对滑坡失稳进行推断,但此方法仅适用于已有滑坡前兆的情况下,实时性差且预报精度低。第二阶段,20世纪80年代,国内外许多学者大量引入各种数学方法与理论模型,例如,Hoek据智利Chuqicamata矿滑坡监测时间-位移曲线提出了外延法,根据所建的模型作外推预报,但其仅能对滑坡近期行为进行推测,存在一定的限制。第三阶段,Haruyama M等基于灰色系统理论开始被应用于滑坡预报研究,但因其拟合预测距实测值偏差较大,由于灰导数等概念引入白化微分方程,使得预测精度不高。第四阶段,BP人工神经网络算法应用于滑坡灾害预报领域中,此种算法具有非线性映射、自适应性学习、大规模并行处理、较强容错性等优点,对滑坡预测预报方面有很高的研究价值。但同时也存在弊端,如神经网络容易陷入局部最优问题、收敛速度慢等。而深度学习作为第三代神经网络,能够从繁多无标签样本中提取出特征量,凭借其多层交互结构及逐层训练的稳定性,将样本进行训练,更适合滑坡灾害预测模型建立,增加预报精确度。
鉴于此,本发明提出一种基于深度信念网络(Deep BeliefNetwork,DBN)的滑坡体灾害预报方法,首先建立滑坡体监测预警系统,运用平均影响值(Mean Impact Value,MIV)算法筛选,筛选出主要灾害诱发因子,然后对两层RBM(Restricted Boltzmann Machine,RBM)进行预训练并提取特征信息,利用遗传算法对模型调优,采用三层BP网络(Back-Propagation Network,BP)模拟滑坡灾害诱发因子与灾害发生等级的关系进行滑坡预报。
发明内容
本发明目的是提供一种基于深度信念网络的滑坡灾害预报方法,解决现有灾害预报中所采用的算法收敛速度慢,预报精确度不高的问题,通过提取特征灾害诱发因子,加快收敛速度,防止陷入局部最优,提升滑坡灾害预报精度。
本发明所采用的技术方案是,基于深度信念网络的滑坡灾害预报方法,具体的操作步骤如下:
步骤1.建立滑坡体监测预警系统,采集大量灾害诱发因子,运用MIV算法筛选,筛选出主要灾害诱发因子;
步骤2.将筛选出的灾害诱发因子进行标准化过程,按特定比例分为测试样本、训练样本和调优样本;
步骤3.构建基于深度信念网络的滑坡灾害预报模型,设置结构为两层RBM和三层BP网络组成;
步骤4.采用CD算法对RBM预训练,更新网络参数;
步骤5.采用遗传算法监督训练学习,使DBN网络整体最优;
步骤6.重构已优化的滑坡灾害预报模型,对输出结果进行滑坡等级划分,预报出滑坡发生可能性。
本发明的其他特点还在于,
步骤1中的滑坡体监测预报系统包括现场预警终端和若干个无线智能传感器模块,现场预警终端与PC端控制中心信号连接,现场预警终端与若干个无线智能传感器模块通过Zigbee无线传输方式通讯;现场预警终端包括主控芯片STM32,主控芯片STM32分别连接有触摸屏模块、USB接口、储存模块、GPRS模块和Zigbee模块;无线智能传感器模块包括主控芯片MSP430,主控芯片MSP430分别连接有存储模块、RS232模块、RS485模块、Zigbee无线通讯模块以及模拟量、数字量以及脉冲量的数据接口。
步骤1中MIV算法筛选灾害诱发因子的具体过程为:
步骤1.1获取智能无线传感器模块采集到的初步完整数据,设样本矩阵为B,并使样本中的每个变量分别加、减10%,构成两个新的训练样本矩阵B1,B2:
步骤1.2如公式1所示,运用拉格朗日插值法将样本矩阵B1和B2进行运算得到对应的结果L1和L2,并将L1和L2用MATLAB进行仿真得出MIV的值,根据滑坡重要影响因素进行排序,从而筛选出MIV主要灾害诱发因子;
其中,L1,L2表示每个变量的差值,Xk表示均值,Xi表示变量。
步骤2中标准化的过程如下:
步骤2.1令输入量X1=土壤含水率,X2=降雨量,X3=裂缝位移,X4=土压力,X5=孔隙水压力,X6=滑体滑面倾角,X7=次声频率;
步骤2.2对采集到的数据进行归一化处理,如公式2所示:
步骤2.3依据特定比例将处理后的样本数据分为预训样本、调优样本、测试样本;其中,在滑坡预报应用中,若需要采集数据特别多的时候,特定比例取98:1:1,反之则取6:2:2;预训练样本包含大量无标签样本数据,调优样本和测试样本包含少量带有类标签的样本数据。
步骤3的具体过程如下:
步骤3.1构建DBN网络模型,它的关键组成部分是RBM,一层RBM包括一个隐含层H和一个可视层V,定义相同层间节点两两互联,相同层内节点不相接,V层有m个可视单元v=(v1,v2,v3,L vm),H层有n个隐含单元h=(h1,h2,h3,L hn);
一层RBM的能量函数可表示为,
其中,vi为可视单元i的取值,hj为隐含单元j的取值;当i,j=0时处于未激活状态,当i,j=1时处于激活状态,ai表示可视偏执层的偏执量,bj为隐含层的偏执量,ωij为可视单元vi与隐含单元hj之间的连接权重,n是隐含层节点数量,m为可视层节点数量,θ=(a,b,ω)为RBM的模型参数;
基于RBM的能量函数,(v,h)的联合概率分布可以表示为:
其中,为归一化因子,即,剖分函数;
当给定可视单元的状态时,隐含单元的激活概率是条件独立的,记第j个隐单元的激活概率可以表示为:
同理,当给定隐含单元的状态时,可视单元的激活概率是条件独立的,记第i个可视单元的激活概率可以表示为:
其中,σ(x)表示激活函数,如公式7所示,当x>0时,σ(x)=1,当x<0时,σ(x)=0;
对于一层RBM所包含的m个可见单元和n个隐单元,根据独立性,(v,h)的条件概率可表示为:
公式8中给定可视单元节点,得到隐含单元节点;
公式9中给定隐含单元节点,得到可视单元节点;
结合公式3,两层RBM的能量函数如公式10所示:
其中,hj0表示第一层隐单元的第j个节点的值,hk2表示第二个隐单元的第k个节点的值;
其联合概率分布由公式4和公式10计算得到:
根据公式11求解出条件分布概率,如公式12所示:
步骤3.2三层BP网络用于输出滑坡等级分类预报结果,将第二层RBM2的隐蔽层h2节点作为BP神经网络的可见层输入端,输出层由5个结点单元组成,数据输入层为v1,由此输入采集到的滑坡灾害诱发因子样本数据,设定其模型输入为7维。
步骤4的具体过程如下:
步骤4.1确定网络参数:样本x={x1,x2,x3,x4,x5,x6,x7}T,学习速率ε,隐含层单元数n;参数初始化:网络参数θ,可视层单元初始值v1=x,并设置最大迭代次数E,极小值化ω,对应偏置值a与b;
步骤4.2开始预训练:
对隐含层所有隐单元根据公式6计算,从P(h1j|v1)抽取h1j∈{0,1};对可视层所有可视单元根据公式5,从P(v2i|h1)中提取v2i∈{0,1};根据公式6,对隐含层所有隐单元进行计算;
步骤4.3参数更新:
ω←ω+ε[(P(h1=1|v1)v1 T-P(h2=1|v2)v2 T]
a←a+ε(v1-v2)
b←b+ε(P(h1=1|v1)-P(h2=1|v2))
步骤4.4重复执行步骤4.2和步骤4.3,完成E次迭代后结束,完成参数更新,最后一次的迭代,将上一层RBM训练完成的状态h2作为下一层RBM的输入v1,训练下一个RBM。
步骤5的具体计算过程如下:
步骤5.1将整个网络等价为BP神经网络,其网络初始权值与偏置在预训练中已设定,GA和BP参数设置如下:种群规模Gs=100,交叉概率Pc0=0.65,变异概率Pm0=0.01,终止代数T=500,动量因子α=0.6;
步骤5.2采用染色体来记录神经网络各层之间的阈值和连接权重,以编码方式生成初始种群;
步骤5.3将步骤1中获得的主要灾害诱发因子经归一化后,输入到BP神经网络的输入层,即v1层,其评价函数表示为:
其中,e(i)=yi d-yi,N为学习样本,yi d为第i组样本的理想输出值,为第i组样本的实际输出值,e(i)为两者之间的误差;
步骤5.4通过适应度函数对网络权值和阈值进行评价,舍弃适应值较低的权值和阈值,保留适应值较高的权值和阈值并进行自我复制;
步骤5.5在高性能的权值和阈值群体中进行交叉和变异两种操作,再按适应度函数执行对权值和阈值群体评价,直至权值和阈值达到期望或达到最大迭代次数停止,至此基于深度信念网络的滑坡灾害预报模型训练完毕;否则继续重复步骤5.3-步骤5.5。
本发明的有益效果是,基于深度信念网络的滑坡灾害预报方法,解决现有灾害预报中所采用的算法收敛速度慢,预报精确度不高的问题;具有以下的有益效果:(1)建立滑坡监测预警系统,能够通过传感器采集到大量的灾害诱发因子,并通过MIV算法筛选出主要灾害诱发因子,加快提取速率,提升预报精确性;(2)DBN对诱发因子样本数据进行特征学习,提升了数据可用性,进一步提升预报精确度;(3)利用遗传算法对整个DBN模型进行调优,降低收敛时间,防止陷入局部最优,保证网络整体最优。
附图说明
图1是本发明的基于深度信念网络的滑坡灾害预报方法的流程图;
图2是本发明的基于深度信念网络的滑坡灾害预报方法所使用的滑坡体监测预警系统结构示意图;
图3是本发明的基于深度信念网络的滑坡灾害预报方法所使用的滑坡体监测预警系统中RBM结构示意图;
图4是深度信念网络模型结构图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明的基于深度信念网络的滑坡灾害预报方法,如图1所示具体的操作步骤如下:
步骤1.建立滑坡体监测预警系统,采集大量灾害诱发因子,运用MIV算法筛选,筛选出主要灾害诱发因子;
步骤2.将筛选出的灾害诱发因子进行标准化过程,按特定比例分为测试样本、训练样本和调优样本;
步骤3.构建基于深度信念网络的滑坡灾害预报模型,设置结构为两层RBM和三层BP网络组成;
步骤4.采用CD算法对RBM预训练,更新网络参数;
步骤5.采用遗传算法监督训练学习,使DBN网络整体最优;
步骤6.重构已优化的滑坡灾害预报模型,对输出结果进行滑坡等级划分,预报出滑坡发生可能性。
步骤1中的滑坡体监测预报系统,如图2所示,包括现场预警终端和若干个无线智能传感器模块,所述现场预警终端与PC端控制中心信号连接,所述现场预警终端与若干个无线智能传感器模块通过Zigbee无线传输方式通讯;所述现场预警终端包括主控芯片STM32,所述主控芯片STM32分别连接有触摸屏模块、USB接口、储存模块、GPRS模块和Zigbee模块;所述无线智能传感器模块包括主控芯片MSP430,所述主控芯片MSP430分别连接有存储模块、RS232模块、RS485模块、Zigbee无线通讯模块以及模拟量、数字量以及脉冲量数据接口。
步骤1中MIV算法筛选灾害诱发因子的具体过程为:
步骤1.1获取智能无线传感器模块采集到的初步完整数据,设样本矩阵为B,并使样本中的每个变量分别加、减10%,构成两个新的训练样本矩阵B1,B2:
步骤1.2如公式1所示,运用拉格朗日插值法将样本矩阵B1和B2进行运算得到对应的结果L1和L2,,并将L1和L2用MATLAB进行仿真得出MIV的值,根据滑坡重要影响因素进行排序,从而筛选出MIV主要灾害诱发因子。
其中,L1,L2表示每个变量的差值,Xk表示均值,Xi表示变量。
步骤2中标准化的过程如下:
步骤2.1令输入量X1=土壤含水率,X2=降雨量,X3=裂缝位移,X4=土压力,X5=孔隙水压力,X6=滑体滑面倾角,X7=次声频率;其输出量为:Y1=成灾概率0~20%,Y2=成灾概率20~40%,Y3=成灾概率40~60%,Y4=成灾概率60~80%,Y5=成灾概率80~100%,表1为成灾概率及其对应的编码
表1成灾概率及其编码
步骤2.2对采集到的数据进行归一化处理,如公式2所示:
步骤2.3依据特定比例将处理后的样本数据分为预训样本、调优样本、测试样本;其中,预训练样本包含大量无标签样本数据,调优样本和测试样本包含少量带有类标签的样本数据。其中的特定比例指的是:人们运用机器学习传统方法的时候,一般将训练集、测试集和验证集,划为6:2:2,若处理少量数据,这一比例划分相对来说具有科学性;在深度学习中若是数据很大,则可以将训练集、验证集、测试集比例调整为98:1:1。因此,在滑坡预报应用中,若需要采集数据特别多的时候,特定比例取98:1:1,反之则6:2:2。
步骤3的具体过程如下:
步骤3.1构建DBN网络模型,它的关键组成部分是RBM,一层RBM包括一个隐含层H和一个可视层V,定义相同层间节点两两互联,相同层内节点不相接,结构如图3所示,V层有m个可视单元v=(v1,v2,v3,L vm),H层有n个隐含单元h=(h1,h2,h3,L hn);
类似于热动力学原理,一层RBM的能量函数表示为,
其中,vi为可视单元i的取值,hj为隐含单元j的取值;当i,j=0时处于未激活状态,当i,j=1时处于激活状态,ai表示可视偏执层的偏执量,bj为隐含层的偏执量,ωij为可视单元vi与隐含单元hj之间的连接权重,n是隐含层节点数量,m为可视层节点数量,θ=(a,b,ω)为RBM的模型参数;
基于RBM的能量函数,(v,h)的联合概率分布可以表示为:
其中,为归一化因子,即,剖分函数;
由RBM的性质,当给定可视单元的状态时,隐含单元的激活概率是条件独立的,记第j个隐单元的激活概率可以表示为:
同理,当给定隐含单元的状态时,可视单元的激活概率是条件独立的,记第i个可视单元的激活概率可以表示为:
其中,σ(x)表示激活函数,如公式7所示,当x>0时,σ(x)=1,当x<0时,σ(x)=0;
对于一层RBM所包含的m个可见单元和n个隐单元,根据独立性,(v,h)的条件概率可表示为:
公式8中给定可视单元节点,得到隐含单元节点;
公式9中给定隐含单元节点,得到可视单元节点;
本发明将由两层RBM和三层BP网络构成深度信念网络模型,如图4所示:结合公式3,两层RBM的能量函数如公式10所示;
其中,hj0表示第一层隐单元的第j个节点的值,hk2表示第二个隐单元的第k个节点的值;
其联合概率分布由公式4和公式10计算得到:
根据公式11求解出条件分布概率,如公式12所示:
步骤3.2三层BP网络用于输出滑坡等级分类预报结果,将RBM2隐蔽层h2节点作为BP神经网络的可见层输入端,输出层由5个结点单元组成,数据输入层为v1,由此输入采集到的滑坡灾害诱发因子样本数据,设定其模型输入为7维(即m=7)。
步骤4的具体过程如下:
步骤4.1确定网络参数:样本x={x1,x2,x3,x4,x5,x6,x7}T,学习速率ε,隐含层单元数n;参数初始化:网络参数θ,可视层单元初始值v1=x,并设置最大迭代次数E,极小值化ω,对应偏置值a与b,其中,在RBM1中,可视单元数为m=7,隐单元数目提前设定为n=30,训练周期T=50,学习速率ε=0.01;
步骤4.2开始预训练:
对隐含层所有隐单元根据公式6计算,从P(h1j|v1)抽取h1j∈{0,1};对可视层所有可视单元根据公式5,从P(v2i|h1)中提取v2i∈{0,1};根据公式6,对隐含层所有隐单元进行计算;
步骤4.3参数更新:
ω←ω+ε[(P(h1=1|v1)v1 T-P(h2=1|v2)v2 T]
a←a+e(v1-v2)
b←b+e(P(h1=1|v1)-P(h2=1|v2))
步骤4.4重复执行步骤4.2和步骤4.3,完成E次迭代后结束,完成参数更新,最后一次的迭代,将上一层RBM训练完成的状态h2作为下一层RBM的输入v1,训练下一个RBM;
步骤5的具体计算过程如下:
步骤5.1将整个网络等价为BP神经网络,其网络初始权值与偏置在预训练中已设定,GA和BP参数设置如下:种群规模Gs=100,交叉概率Pc0=0.65,变异概率Pm0=0.01,终止代数T=500,动量因子α=0.6;
步骤5.2采用染色体来记录神经网络各层之间的阈值和连接权重,以编码方式生成初始种群;
步骤5.3为了从染色体种群中选择最优染色体来初始化神经网络的阈值和连接权值,因此,将步骤1中获得的主要灾害诱发因子经归一化后,输入到BP神经网络的输入层,即v1层,其评价函数表示为:
其中,e(i)=yi d-yi,N为学习样本,yi d为第i组样本的理想输出值,为第i组样本的实际输出值,e(i)为两者之间的误差;
步骤5.4通过适应度函数对网络权值和阈值进行评价,舍弃适应值较低的权值和阈值,保留适应值较高的权值和阈值并进行自我复制;
步骤5.5在高性能的权值和阈值群体中进行交叉和变异两种操作,再按适应度函数执行对权值和阈值群体评价,直至权值和阈值达到期望或达到最大迭代次数停止,至此基于深度信念网络的滑坡灾害预报模型训练完毕;否则继续重复步骤5.3-步骤5.5。
步骤6中预警等级划分如表2所示,
表2预警对应等级划分
经过试验发现,利用MIV算法筛选灾害诱发因子,能够及时筛选主要灾害因子,提升预报时效性;并且构建的网络结构层数和样本训练周期与预测准确率存在一定关系,本发明采用2层RBM网络结构,训练周期定为200时,预报精确度达到较好的效果;采用遗传算法调优,防止BP网络陷入局部最优,从而进一步提升预报准确度。
Claims (7)
1.基于深度信念网络的滑坡灾害预报方法,其特征在于,具体的操作步骤如下:
步骤1.建立滑坡体监测预警系统,采集大量灾害诱发因子,运用MIV算法筛选,筛选出主要灾害诱发因子;
步骤2.将筛选出的灾害诱发因子进行标准化过程,按特定比例分为测试样本、训练样本和调优样本;
步骤3.构建基于深度信念网络的滑坡灾害预报模型,设置结构为两层RBM和三层BP网络组成;
步骤4.采用CD算法对RBM预训练,更新网络参数;
步骤5.采用遗传算法监督训练学习,使DBN网络整体最优;
步骤6.重构已优化的滑坡灾害预报模型,对输出结果进行滑坡等级划分,预报出滑坡发生可能性。
2.如权利要求1所述的基于深度信念网络的滑坡灾害预报方法,其特征在于,所述步骤1中的滑坡体监测预报系统包括现场预警终端和若干个无线智能传感器模块,所述现场预警终端与PC端控制中心信号连接,所述现场预警终端与若干个无线智能传感器模块通过Zigbee无线传输方式通讯;所述现场预警终端包括主控芯片STM32,所述主控芯片STM32分别连接有触摸屏模块、USB接口、储存模块、GPRS模块和Zigbee模块;所述无线智能传感器模块包括主控芯片MSP430,所述主控芯片MSP430分别连接有存储模块、RS232模块、RS485模块、Zigbee无线通讯模块以及模拟量、数字量以及脉冲量的数据接口。
3.如权利要求1所述的基于深度信念网络的滑坡灾害预报方法,其特征在于,所述步骤1中MIV算法筛选灾害诱发因子的具体过程为:
步骤1.1获取智能无线传感器模块采集到的初步完整数据,设样本矩阵为B,并使样本中的每个变量分别加、减10%,构成两个新的训练样本矩阵B1,B2:
步骤1.2如公式1所示,运用拉格朗日插值法将样本矩阵B1和B2进行运算得到对应的结果L1和L2,并将L1和L2用MATLAB进行仿真得出MIV的值,根据滑坡重要影响因素进行排序,从而筛选出MIV主要灾害诱发因子;
其中,L1,L2表示每个变量的差值,Xk表示均值,Xi表示变量。
4.如权利要求1所述的基于深度信念网络的滑坡灾害预报方法,其特征在于,所述步骤2中标准化的过程如下:
步骤2.1令输入量X1=土壤含水率,X2=降雨量,X3=裂缝位移,X4=土压力,X5=孔隙水压力,X6=滑体滑面倾角,X7=次声频率;
步骤2.2对采集到的数据进行归一化处理,如公式2所示:
步骤2.3依据特定比例将处理后的样本数据分为预训样本、调优样本、测试样本;其中,在滑坡预报应用中,若需要采集数据特别多的时候,特定比例取98:1:1,反之则取6:2:2;预训练样本包含大量无标签样本数据,调优样本和测试样本包含少量带有类标签的样本数据。
5.如权利要求1所述的基于深度信念网络的滑坡灾害预报方法,其特征在于,所述步骤3的具体过程如下:
步骤3.1构建DBN网络模型,它的关键组成部分是RBM,一层RBM包括一个隐含层H和一个可视层V,定义相同层间节点两两互联,相同层内节点不相接,V层有m个可视单元v=(v1,v2,v3,L vm),H层有n个隐含单元h=(h1,h2,h3,L hn);
一层RBM的能量函数可表示为:
其中,vi为可视单元i的取值,hj为隐含单元j的取值;当i,j=0时处于未激活状态,当i,j=1时处于激活状态,ai表示可视偏执层的偏执量,bj为隐含层的偏执量,ωij为可视单元vi与隐含单元hj之间的连接权重,n是隐含层节点数量,m为可视层节点数量,θ=(a,b,ω)为RBM的模型参数;
基于RBM的能量函数,(v,h)的联合概率分布可以表示为:
其中,为归一化因子,即,剖分函数;
当给定可视单元的状态时,隐含单元的激活概率是条件独立的,记第j个隐单元的激活概率可以表示为:
同理,当给定隐含单元的状态时,可视单元的激活概率是条件独立的,记第i个可视单元的激活概率可以表示为:
其中,σ(x)表示激活函数,如公式7所示,当x>0时,σ(x)=1,当x<0时,σ(x)=0;
对于一层RBM所包含的m个可见单元和n个隐单元,根据独立性,(v,h)的条件概率可表示为:
公式8中给定可视单元节点,得到隐含单元节点;
公式9中给定隐含单元节点,得到可视单元节点;
结合公式3,两层RBM的能量函数如公式10所示:
其中,hj0表示第一层隐单元的第j个节点的值,hk2表示第二个隐单元的第k个节点的值;
其联合概率分布由公式4和公式10计算得到:
根据公式11求解出条件分布概率,如公式12所示:
步骤3.2三层BP网络用于输出滑坡等级分类预报结果,将第二层RBM2的隐蔽层h2节点作为BP神经网络的可见层输入端,输出层由5个结点单元组成,数据输入层为v1,由此输入采集到的滑坡灾害诱发因子样本数据,设定其模型输入为7维。
6.如权利要求5所述的基于深度信念网络的滑坡灾害预报方法,其特征在于,所述步骤4的具体过程如下:
步骤4.1确定网络参数:样本x={x1,x2,x3,x4,x5,x6,x7}T,学习速率ε,隐含层单元数n;参数初始化:网络参数θ,可视层单元初始值v1=x,并设置最大迭代次数E,极小
步骤4.2开始预训练:值化ω,对应偏置值a与b;
对隐含层所有隐单元根据公式6计算,从P(h1j|v1)抽取h1j∈{0,1};对可视层所有可视单元根据公式5,从P(v2i|h1)中提取v2i∈{0,1};根据公式6,对隐含层所有隐单元进行计算;
步骤4.3参数更新:
ω←ω+ε[(P(h1=1|v1)v1 T-P(h2=1|v2)v2 T]
a←a+ε(v1-v2)
b←b+ε(P(h1=1|v1)-P(h2=1|v2))
步骤4.4重复执行步骤4.2和步骤4.3,完成E次迭代后结束,完成参数更新,最后一次的迭代,将上一层RBM训练完成的状态h2作为下一层RBM的输入v1,训练下一个RBM。
7.如权利要求1所述的基于深度信念网络的滑坡灾害预报方法,其特征在于,所述步骤5的具体计算过程如下:
步骤5.1将整个网络等价为BP神经网络,其网络初始权值与偏置在预训练中已设定,GA和BP参数设置如下:种群规模Gs=100,交叉概率Pc0=0.65,变异概率Pm0=0.01,终止代数T=500,动量因子α=0.6;
步骤5.2采用染色体来记录神经网络各层之间的阈值和连接权重,以编码方式生成初始种群;
步骤5.3将步骤1中获得的主要灾害诱发因子经归一化后,输入到BP神经网络的输入层,即v1层,其评价函数表示为:
其中,e(i)=yi d-yi,N为学习样本,yi d为第i组样本的理想输出值,为第i组样本的实际输出值,e(i)为两者之间的误差;
步骤5.4通过适应度函数对网络权值和阈值进行评价,舍弃适应值较低的权值和阈值,保留适应值较高的权值和阈值并进行自我复制;
步骤5.5在高性能的权值和阈值群体中进行交叉和变异两种操作,再按适应度函数执行对权值和阈值群体评价,直至权值和阈值达到期望或达到最大迭代次数停止,至此基于深度信念网络的滑坡灾害预报模型训练完毕;否则继续重复步骤5.3-步骤5.5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810667868.5A CN109063247B (zh) | 2018-06-26 | 2018-06-26 | 基于深度信念网络的滑坡灾害预报方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810667868.5A CN109063247B (zh) | 2018-06-26 | 2018-06-26 | 基于深度信念网络的滑坡灾害预报方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109063247A true CN109063247A (zh) | 2018-12-21 |
CN109063247B CN109063247B (zh) | 2023-04-18 |
Family
ID=64821584
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810667868.5A Active CN109063247B (zh) | 2018-06-26 | 2018-06-26 | 基于深度信念网络的滑坡灾害预报方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109063247B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109637090A (zh) * | 2019-01-30 | 2019-04-16 | 深圳市地质局 | 一种基于soa架构的灾害监测预警平台建设方法 |
CN111639813A (zh) * | 2020-06-01 | 2020-09-08 | 福州大学 | 基于深度学习的弃渣场危险性预警方法及系统 |
CN111784070A (zh) * | 2020-07-09 | 2020-10-16 | 中国地质环境监测院 | 一种基于XGBoost算法的滑坡短临智能预警方法 |
CN112232486A (zh) * | 2020-10-19 | 2021-01-15 | 南京宁麒智能计算芯片研究院有限公司 | 一种yolo脉冲神经网络的优化方法 |
CN112308214A (zh) * | 2020-11-09 | 2021-02-02 | 海南核电有限公司 | 一种用于冷源致灾物灾害预测的深度学习改进算法 |
CN112735094A (zh) * | 2020-12-17 | 2021-04-30 | 中国地质环境监测院 | 基于机器学习的地质灾害预测方法、装置和电子设备 |
CN112949189A (zh) * | 2021-03-07 | 2021-06-11 | 湖南安全技术职业学院(长沙煤矿安全技术培训中心) | 一种基于深度学习的多因素诱发滑坡预测的建模方法 |
CN112966722A (zh) * | 2021-02-07 | 2021-06-15 | 南昌大学 | 基于半监督随机森林模型的区域滑坡易发性预测方法 |
CN113065455A (zh) * | 2021-03-30 | 2021-07-02 | 中国水利水电科学研究院 | 一种基于深度学习的滑坡风险巡检方法及系统 |
CN113470316A (zh) * | 2021-05-24 | 2021-10-01 | 电子科技大学 | 一种基于自编码单分类模型的泥石流监测预警方法 |
CN113822522A (zh) * | 2021-06-23 | 2021-12-21 | 中国科学院空天信息创新研究院 | 一种滑坡易发性评估方法、装置、设备及可读存储介质 |
CN114821975A (zh) * | 2022-06-02 | 2022-07-29 | 重庆蜀通岩土工程有限公司 | 一种高边坡安全在线监测方法 |
CN116704708A (zh) * | 2023-08-08 | 2023-09-05 | 山东省减灾中心 | 一种基于大数据的智能自然灾害预警系统及方法 |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH06325270A (ja) * | 1993-05-18 | 1994-11-25 | Hochiki Corp | 火災判定装置及びニューラルネットワークの学習方法 |
US20060020992A1 (en) * | 2002-11-12 | 2006-01-26 | Pugel Michael A | Weather/disaster alert system using a data network |
KR100880715B1 (ko) * | 2008-03-26 | 2009-02-02 | (주)성경기술공사 | 수치지도를 이용한 산사태 발생지 예측 방법 |
US20090205015A1 (en) * | 2006-06-02 | 2009-08-13 | Javier Baliosian | Method for Forecasting Unstable Policy Enforcements |
CN104951425A (zh) * | 2015-07-20 | 2015-09-30 | 东北大学 | 一种基于深度学习的云服务性能自适应动作类型选择方法 |
CN105869100A (zh) * | 2016-04-29 | 2016-08-17 | 中国地质大学(武汉) | 一种基于大数据思维的滑坡多场监测数据的融合及预测方法 |
CN106815971A (zh) * | 2016-12-21 | 2017-06-09 | 西安工程大学 | 一种基于rbf神经网络的滑坡灾害预报方法 |
WO2017129030A1 (zh) * | 2016-01-29 | 2017-08-03 | 阿里巴巴集团控股有限公司 | 磁盘的故障预测方法和装置 |
CN107203810A (zh) * | 2017-05-22 | 2017-09-26 | 河海大学 | 一种基于深度网络的降水量预测方法 |
CN107247259A (zh) * | 2017-07-10 | 2017-10-13 | 西安电子科技大学 | 基于神经网络的k分布海杂波形状参数估计方法 |
CN107423670A (zh) * | 2017-04-20 | 2017-12-01 | 上海交通大学 | 基于深度置信网络的modis雾监测方法 |
CN107506590A (zh) * | 2017-08-26 | 2017-12-22 | 郑州大学 | 一种基于改进深度信念网络的心血管疾病预测模型 |
CN107656152A (zh) * | 2017-09-05 | 2018-02-02 | 西安工程大学 | 一种基于ga‑svm‑bp变压器故障诊断方法 |
US20180157977A1 (en) * | 2016-12-05 | 2018-06-07 | Tata Consultancy Services Limited | Training inductive logic programming enhanced deep belief network models for discrete optimization |
-
2018
- 2018-06-26 CN CN201810667868.5A patent/CN109063247B/zh active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH06325270A (ja) * | 1993-05-18 | 1994-11-25 | Hochiki Corp | 火災判定装置及びニューラルネットワークの学習方法 |
US20060020992A1 (en) * | 2002-11-12 | 2006-01-26 | Pugel Michael A | Weather/disaster alert system using a data network |
US20090205015A1 (en) * | 2006-06-02 | 2009-08-13 | Javier Baliosian | Method for Forecasting Unstable Policy Enforcements |
KR100880715B1 (ko) * | 2008-03-26 | 2009-02-02 | (주)성경기술공사 | 수치지도를 이용한 산사태 발생지 예측 방법 |
CN104951425A (zh) * | 2015-07-20 | 2015-09-30 | 东北大学 | 一种基于深度学习的云服务性能自适应动作类型选择方法 |
WO2017129030A1 (zh) * | 2016-01-29 | 2017-08-03 | 阿里巴巴集团控股有限公司 | 磁盘的故障预测方法和装置 |
CN105869100A (zh) * | 2016-04-29 | 2016-08-17 | 中国地质大学(武汉) | 一种基于大数据思维的滑坡多场监测数据的融合及预测方法 |
US20180157977A1 (en) * | 2016-12-05 | 2018-06-07 | Tata Consultancy Services Limited | Training inductive logic programming enhanced deep belief network models for discrete optimization |
CN106815971A (zh) * | 2016-12-21 | 2017-06-09 | 西安工程大学 | 一种基于rbf神经网络的滑坡灾害预报方法 |
CN107423670A (zh) * | 2017-04-20 | 2017-12-01 | 上海交通大学 | 基于深度置信网络的modis雾监测方法 |
CN107203810A (zh) * | 2017-05-22 | 2017-09-26 | 河海大学 | 一种基于深度网络的降水量预测方法 |
CN107247259A (zh) * | 2017-07-10 | 2017-10-13 | 西安电子科技大学 | 基于神经网络的k分布海杂波形状参数估计方法 |
CN107506590A (zh) * | 2017-08-26 | 2017-12-22 | 郑州大学 | 一种基于改进深度信念网络的心血管疾病预测模型 |
CN107656152A (zh) * | 2017-09-05 | 2018-02-02 | 西安工程大学 | 一种基于ga‑svm‑bp变压器故障诊断方法 |
Non-Patent Citations (3)
Title |
---|
CHIRAG I PATEL: "《Handwritten Character Recognition using Neural Network》" * |
周世超: "《基于深度信念网络的语音服务文本分类》", 《计算机工程与应用》 * |
王鲁昆: "《带高斯噪声的受限玻尔兹曼机在变压器故障诊断中的应用》", 《中国海洋大学学报》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109637090A (zh) * | 2019-01-30 | 2019-04-16 | 深圳市地质局 | 一种基于soa架构的灾害监测预警平台建设方法 |
CN111639813A (zh) * | 2020-06-01 | 2020-09-08 | 福州大学 | 基于深度学习的弃渣场危险性预警方法及系统 |
CN111784070A (zh) * | 2020-07-09 | 2020-10-16 | 中国地质环境监测院 | 一种基于XGBoost算法的滑坡短临智能预警方法 |
CN112232486A (zh) * | 2020-10-19 | 2021-01-15 | 南京宁麒智能计算芯片研究院有限公司 | 一种yolo脉冲神经网络的优化方法 |
CN112308214B (zh) * | 2020-11-09 | 2024-07-16 | 海南核电有限公司 | 一种用于冷源致灾物灾害预测的深度学习改进算法 |
CN112308214A (zh) * | 2020-11-09 | 2021-02-02 | 海南核电有限公司 | 一种用于冷源致灾物灾害预测的深度学习改进算法 |
CN112735094A (zh) * | 2020-12-17 | 2021-04-30 | 中国地质环境监测院 | 基于机器学习的地质灾害预测方法、装置和电子设备 |
CN112966722A (zh) * | 2021-02-07 | 2021-06-15 | 南昌大学 | 基于半监督随机森林模型的区域滑坡易发性预测方法 |
CN112949189A (zh) * | 2021-03-07 | 2021-06-11 | 湖南安全技术职业学院(长沙煤矿安全技术培训中心) | 一种基于深度学习的多因素诱发滑坡预测的建模方法 |
CN113065455A (zh) * | 2021-03-30 | 2021-07-02 | 中国水利水电科学研究院 | 一种基于深度学习的滑坡风险巡检方法及系统 |
CN113065455B (zh) * | 2021-03-30 | 2021-09-17 | 中国水利水电科学研究院 | 一种基于深度学习的滑坡风险巡检方法及系统 |
CN113470316A (zh) * | 2021-05-24 | 2021-10-01 | 电子科技大学 | 一种基于自编码单分类模型的泥石流监测预警方法 |
CN113470316B (zh) * | 2021-05-24 | 2022-03-25 | 电子科技大学 | 一种基于自编码单分类模型的泥石流监测预警方法 |
CN113822522A (zh) * | 2021-06-23 | 2021-12-21 | 中国科学院空天信息创新研究院 | 一种滑坡易发性评估方法、装置、设备及可读存储介质 |
CN113822522B (zh) * | 2021-06-23 | 2024-09-20 | 中国科学院空天信息创新研究院 | 一种滑坡易发性评估方法、装置、设备及可读存储介质 |
CN114821975A (zh) * | 2022-06-02 | 2022-07-29 | 重庆蜀通岩土工程有限公司 | 一种高边坡安全在线监测方法 |
CN116704708A (zh) * | 2023-08-08 | 2023-09-05 | 山东省减灾中心 | 一种基于大数据的智能自然灾害预警系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109063247B (zh) | 2023-04-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109063247B (zh) | 基于深度信念网络的滑坡灾害预报方法 | |
CN109214575B (zh) | 一种基于小波长短期记忆网络的超短期风电功率预测方法 | |
CN103077267B (zh) | 一种基于改进bp神经网络的参量声源建模方法 | |
CN102622515B (zh) | 一种天气预测方法 | |
CN107993012B (zh) | 一种时间自适应的电力系统在线暂态稳定评估方法 | |
CN106897826A (zh) | 一种交通事故风险评估方法及系统 | |
CN105931116A (zh) | 基于深度学习机制的自动化信用评分系统及方法 | |
CN108399470B (zh) | 一种基于多示例遗传神经网络的室内pm2.5预测方法 | |
CN102622418A (zh) | 一种基于bp神经网络的预测装置及设备 | |
CN110991690A (zh) | 一种基于深度卷积神经网络的多时次风速预测方法 | |
CN113408190B (zh) | 基于Bayes-LSTM模型的公路隧道施工期围岩变形预测方法 | |
CN105243259A (zh) | 基于极限学习机的脉动风速快速预测方法 | |
CN111652425A (zh) | 一种基于粗糙集和长短期记忆网络的河流水质预测方法 | |
CN103617147A (zh) | 一种矿井突水水源层识别方法 | |
CN106355540A (zh) | 一种基于gra‑bp神经网络的中小型水库大坝安全评价方法 | |
CN105259331A (zh) | 一种节理岩体单轴强度预测方法 | |
CN103530818A (zh) | 一种基于brb系统的供水管网建模方法 | |
CN113219871B (zh) | 一种养护室环境参数检测系统 | |
CN112215410B (zh) | 基于改进深度学习的电力负荷预测方法 | |
CN105044775A (zh) | 一种地震流体反演识别方法和装置 | |
CN114548482A (zh) | 一种基于蠕变型滑坡动能变化率临滑预警方法 | |
CN111310331B (zh) | 基于条件变分自编码的地壳模型构造方法 | |
CN115409369A (zh) | 基于机理和数据混合驱动的综合能源系统可靠性评估方法 | |
CN115018331A (zh) | 一种传统村落保护利用价值的智能评价方法及系统 | |
CN105568732A (zh) | 盘磨机控制方法 |
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 |