CN110288700A - 一种岩质边坡结构面自动分组及位移预测方法 - Google Patents
一种岩质边坡结构面自动分组及位移预测方法 Download PDFInfo
- Publication number
- CN110288700A CN110288700A CN201910560846.3A CN201910560846A CN110288700A CN 110288700 A CN110288700 A CN 110288700A CN 201910560846 A CN201910560846 A CN 201910560846A CN 110288700 A CN110288700 A CN 110288700A
- Authority
- CN
- China
- Prior art keywords
- formula
- displacement
- visual field
- slope
- artificial fish
- 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.)
- Pending
Links
- 238000006073 displacement reaction Methods 0.000 title claims abstract description 53
- 239000011435 rock Substances 0.000 title claims abstract description 48
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000012544 monitoring process Methods 0.000 claims abstract description 15
- 235000019688 fish Nutrition 0.000 claims description 46
- 230000000007 visual effect Effects 0.000 claims description 46
- 241000251468 Actinopterygii Species 0.000 claims description 20
- 210000002569 neuron Anatomy 0.000 claims description 19
- 230000008859 change Effects 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000019637 foraging behavior Effects 0.000 claims description 8
- 238000012937 correction Methods 0.000 claims description 6
- 230000006399 behavior Effects 0.000 claims description 4
- 235000013305 food Nutrition 0.000 claims description 4
- 238000012804 iterative process Methods 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 230000006978 adaptation Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 2
- 244000062793 Sorghum vulgare Species 0.000 claims description 2
- 238000012512 characterization method Methods 0.000 claims description 2
- 235000019713 millet Nutrition 0.000 claims description 2
- 238000012360 testing method Methods 0.000 claims description 2
- 238000002485 combustion reaction Methods 0.000 claims 1
- 238000013528 artificial neural network Methods 0.000 abstract description 7
- 238000005516 engineering process Methods 0.000 abstract description 7
- 238000004458 analytical method Methods 0.000 description 3
- 239000002245 particle Substances 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000002068 genetic effect Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B7/00—Measuring arrangements characterised by the use of electric or magnetic techniques
- G01B7/02—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring length, width or thickness
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
-
- 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/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T19/00—Manipulating 3D models or images for computer graphics
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Geometry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Computer Graphics (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Computational Linguistics (AREA)
- Mathematical Physics (AREA)
- Multimedia (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Hardware Design (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Health & Medical Sciences (AREA)
- Image Analysis (AREA)
Abstract
本发明的一种岩质边坡结构面自动分组及位移预测方法,包括:通过无人机获取岩质边坡的不同角度的多张图像,对不同角度的多张图像进行集成,构建岩质边坡的三维模型;对岩质边坡的结构面分组;对岩质边坡位移进行预测。本发明基于UAV(Unmanned Aerial Vehicle,无人机)技术提出的一种新型的岩质边坡结构面采集方法,实现了对结构面信息高精度、高效率的获取;同时采用一种新型的分类方法对UAV测量所得的结构面进行分组;最后通过UAV技术和神经网络方法实现了对边坡位移的监测与预测。以解决岩质边坡灾害预警问题。
Description
技术领域
本发明属于岩质边坡灾害预警技术领域,涉及一种岩质边坡结构面自动分组及位移预测方法。
背景技术
对边坡进行地质调查以及产状统计对于研究岩质边坡稳定问题具有非常大的意义,传统测量方法测线法、窗口统计法、借助罗盘等工具直接对边坡进行近地测量,但是效率低下,工作量偏大且具有一定的危险性,对于高陡边坡来说,测量人员是无法近距离接近边坡的。三维激光和近景摄影测量作为两种新兴的测量技术,基于三维激光技术对边坡的进行扫描,获取了边坡的三维点云模型,并对结构面信息进行获取,并对边坡结构面进行聚类分组。但也存在一定的局限性,三维激光扫描可以存储海量云点坐标,但是由于设备不便于携带,且花费成本较高,在测量中往往使用的较少,而Shapematrix3D***测量虽可以采取非接触摄影测量,但是对于测量高陡边坡,由于人工拍摄角度的不确定性,而且拍摄过程中需要在在坡体底部进行标定,具有一定的危险性。岩质边坡是由许多节理面、裂隙面、断层面等各种地质构造切割岩体所形成的非常复杂的结构体系,与土质边坡的构造具有较大的差异,其稳定性不但受制于岩石的强度,更多程度上依赖于边坡结构面的产状,因此对结构面进行统计和分组,得到其分布规律,目前对结构面分组的方法存在较大的误差。
岩质边坡破坏具有一定的突发性,需要对边坡进行实时监测,及时反馈边坡的稳定状况。根据边坡位移的变化规律来准确推演位移信息的变化趋势,能大大减少人力物力的损耗,掌握边坡变形的规律,对超前预测边坡的破坏以及判断边坡的稳定状态具有非常重要的影响。目前虽然监测技术发展得比较快,但监测资料的分析研究还处于待进一步发展的状态。边坡的监测数据虽然多,但是难以与边坡的安全稳定性状态建立直观方便的模型,由于二者之间的相关性的建立是十分复杂和闲难的,因而我们难以通过边坡的监测数据对其安全状况作直接的评判。于理论方法的不足,基于边坡所处的环境的复杂性,并且影响边坡安企稳定性的因素比较多,因此必须在边坡的监测数据反馈分析上加强研究。近年来,随着计算机技术的高速发展,神经网络,遗传算法等也被用于边坡稳定的判别总。关于边坡位移的预报预测方法,一般分为以下几种:1)神经网络与各类算法相结合,包括BP网络和遗传算法,PSO和LSSVM算法,SVM和马尔科夫链。2)考虑边坡位移的模糊性,采用时间序列及其衍生方法关于边坡位移的预测,包括时间序列与灰色***,时间序列与神经网络,混沌时间序列。3)小波变换和神经网络,包括小波变换和神经网络,小波分析和灰色,这几种类型的方法,各有优缺点,应用上具有一定的局限性。
发明内容
本发明的目的是提供一种岩质边坡结构面自动分组及位移预测方法,以解决岩质边坡灾害预警问题。
本发明提供一种岩质边坡结构面自动分组及位移预测方法,包括如下步骤:
步骤1:通过无人机获取岩质边坡的不同角度的多张图像,对不同角度的多张图像进行集成,构建岩质边坡三维模型;
步骤2:对无人机采集获得的边坡三维模型的结构面信息进行分组;
步骤3:根据边坡结构面的分组和无人机采集的数据对岩质边坡位移进行预测。
本发明的种岩质边坡结构面自动分组及位移预测方法,基于UAV(UnmannedAerial Vehicle,无人机)技术,实现了对结构面信息高精度、高效率的获取;同时采用一种新型的分类方法对UAV测量所得的结构面进行分组,并研究了经分组后不同概率分布的结构面产状分布形式下边坡关键块体的稳定情况;最后通过UAV技术和神经网络方法实现了对边坡位移的监测与预测。对岩质边坡的灾害预警达到很好的效果。
附图说明
图1是本发明的一种岩质边坡结构面自动分组及位移预测方法的流程图;
图2是无人机采集和处理图片的流程图;
图3结构面空间坐标系;
图4a为边坡的三维网络模型图;
图4b位移变化方向图;
图5a监测点连线图;
图5b等高线位移变化图;
图6为ASFA-RSK算法流程图。
具体实施方式
如图1所示本发明的一种岩质边坡结构面自动分组及位移预测方法,包括如下步骤:
步骤1:通过无人机获取岩质边坡的不同角度的多张图像,对不同角度的多张图像进行集成,构建岩质边坡的三维模型,如图2所示,所述步骤1具体包括如下步骤:
步骤1.1:地面进行控制点选取,通过无人机对待测岩质边坡进行多角度拍摄,包括垂直拍摄和倾斜拍摄;
步骤1.2:对岩质边坡的不同角度的多张图像进行集成,进行超高密度点云图像构建,以生成点云模型;
步骤1.3:对点云模型进行纹理修正、精细化处理和去噪处理,获得岩质边坡的三维模型;
步骤1.4:获得边坡结构面点云坐标和出露迹长对结构面产状进行拟合,具体提取的结构面参数包括:倾向、倾角、空间坐标x、空间坐标y、空间坐标z、粘聚力和内摩擦角。
步骤2:对无人机采集获得的边坡三维模型的结构面信息进行分组,所述步骤2具体包括:
步骤2.1:结构面视作无厚度的无限延伸的平面,并用结构面法向量来表示,所述步骤2.1具体为:
用倾角α(0≤α≤90o)和倾向β(0≤β≤360o)表征结构面的空间位置,建立如图3所示的XYZ坐标系,X轴代表正东方向,Y轴代表正北方向,Z轴代表极射球面的正上方,球面半径为1,(n1、n2、n3)为结构面法向量投影在X、Y、Z轴的坐标大小,(n1、n2、n3)为垂直于结构面的法向量,将(n1、n2、n3)的坐标值等效在球面坐标系中,按照投影原理,倾角为α,倾向为β的结构面可由垂直于该结构面的单位法向量n表示,法向量n垂直于结构面,投影在坐标面上表达式为:
则结构面法向量可表示:
n=(n1,n2,n3)=(sinαsinβ,sinαcosβ,cosα) (2)。
步骤2.2:将人工鱼群算法与K-Means算法结合,建立AFSA-RSK结构面分类算法,利用鱼群算法对结构面产状聚心集进行搜寻,并通过K-Means算法进行聚类,如图6所示,所述步骤2.2具体为:
1)人工鱼在寻优时,每次移动之后通过评估适应度来确定下一步的行为,当适应度越大,结构面样本中各结构面参数离聚心的距离也就越小,聚类误差也会越小,适应度计算公式为:
式中:k为聚心个数,l为结构面样本个数,||d(ci-xj)||2为分组结构面中样本x离聚心c的距离,距离计算公式为:
式中:nc、nj分别为聚心结构面的单位法向量和j结构面的单位法向量。
通过移动方式更改和搜索聚心坐标:
2)初始人工鱼在Xi位置处时在视野Visual中随机移动到Xj处,在假定初始聚心的情况下,人工鱼在结构面样本机视野内随机搜寻,搜寻公式为:
Xj=Xi+Visual×Rand() (5)
按照上式随机移动后,到达新的位置Xj处,若Xj处结果适应度更大,则向前一步,得到前进公式为:
Xt+1 i为多次执行觅食行为过程中新聚心的位置;
3)人工鱼在Xi位置处搜索视野内的伙伴数量,计算当前视野内鱼群数目的均值坐标,即为中心坐标值,若中心位置适应度较大,输出中心坐标值,同时按式(7)进行移动,得到新的聚心坐标,并与中心坐标值进行比较,较优的记录下来,中心位置计算方法如式(8);
式中,Xj为视野内某结构面的坐标,g为伙伴数目,Xcen为鱼群经过聚群得到的结构面新聚心;
4)人工鱼Xi搜索得到适应度最大伙伴Xmax,移动方式如式(9),按照式(7)可以得到新的聚类中心Xnext的坐标值,将Xnext和Xmax中的较优者记录下来:
式中,Xmax为视野中最大适应度的个体。
步骤2.3:聚类完成后,对人工鱼群算法的步长和视野进行修正,动态调整聚类过程,所述步骤2.3具体为:
1)采用Xie-Beni指标SXB、分类熵指标Hm两项指标检验聚类效果,若满足ASFA-RSK算法的两项指标误差均小于K-Means算法的两项指标误差,将结果记录,并继续迭代,超过最大迭代次数输出聚类结果;
2)若不满足就对步长和视野进行修正,步长修正公式和视野修正公式如式(10),式(12)所示:
Stepmax选为0.16,则相应的改进觅食行为策略为:
随着迭代步长的增加,步长会随着式(11)变化,从而避免迭代过程中跳过最优值;
具体的,视野修正公式为:
式中:Vmax选为1,Tmax为最大迭代次数。
对鱼群视野进行修正,建立视野随迭代次数动态的变化过程函数,有利于寻优精度的提升和运算速度的加快。
步骤2.4:修正视野与步长后,每执行一步与前步结果进行比较,并将两者间误差较小的值记录下来,达到最大迭代次数后,输出记录的最优结果,并与K-Means算法计算结果比较,若优于K-Means算法计算结果,则输出结果;否则,再次修改步长,重新迭代,超过最大迭代次数后输出结果。
步骤3:根据边坡结构面的分组和无人机采集的数据对岩质边坡位移进行预测,所述步骤3具体为:
步骤3.1:位移数据获取,采用无人机对边坡位移进行实时监测,获得一段时间内的边坡位移序列,所述步骤3.1具体为:
对位移方向进行预判,将坡体点云坐标大致移动方向进行统计,在如图4a所示的三维网络模型上用箭头标出位移变化方向,如图4b所示。连接监测点位置,对应得到等高线位置上两监测点连线的点云坐标变化情况,继而得到坡体连线上各点位移点云变化值,如图5a和5b所示。位移计算步骤如下:
1)根据测点进行分组,N个测点对应N个分组,监测标定盘的直径为d,计算标定盘附近点的坐标S,若点K与间距S之间间距小于d,则S、K同为一组;
2)对其他组数据采取与1)相同的计算方式,计算完毕得到N组点云数据;
3)对经建模得到的点云数据进行均值计算,得到中心点的坐标,令中心坐标为Scen(X,Y,Z),并作为该分组的最终坐标值,以中心坐标Scen(X,Y,Z)的变化代表该测点的位置变化,中心坐标Scen(X,Y,Z)计算方法如下式所示:
则每一次测量间隔内,位移变化S为:
其中,t代表测量次数,两次测量的位移变化值为S(t)。
步骤3.2:将人工鱼群算法与Elman网络相结合,建立改进AFSA-Elman边坡位移预测网络,利用改进AFSA-Elman边坡位移预测网络对边坡位移序列进行预测,步骤3.2具体为:
1)利用人工鱼群算法对Elman网络进行改造,鱼群算法寻找的是群体中最优的个体,即食物浓度最大的位置,即误差最小的位置,令网络适应度函数为:
式中:G代表权值变化中适应度最大的位置,Q为网络期望输出与实际输出的误差;
2)选取四层Elman网络,令r为输入层神经元个数,m为隐含层神经元个数,t为承接层神经元个数,任意选取两条人工鱼Rm和Rx,考虑权值变化,利用各层神经元来计算各层连接之间输入与输出的欧氏距离,任意两条人工鱼Rm和Rx之间的距离为:
式(16)中wij为i层隐含层神经元和j层输入层神经元权值,Lgi为g层承接层神经元和i层隐含层权值,Vkg为输入层k层神经元和g层承接层神经元权值,wi0为隐含层神经元阈值,wg0为承接层神经元阈值,wk0为输出层神经元阈值;借鉴粒子群算法中粒子运动时粒子更新速度,对鱼群引入鱼群游速ω,使得AFSA更接近真实鱼群活动过程。
3)对网络进行训练,随机指定某条人工鱼状态为Rm,在鱼群视野范围内随机寻找另一个人工鱼,改变权值变化,令其状态为Rx,若该处食物浓度较大,误差计算,即Gx>Gm,则执行觅食行为,反之,向任意方向进行随机移动,网络权值和阈值被改变,移动公式为:
Lgi、Vkg、wi0、wg0、wk0觅食行为与上述wij方式一致,式(17)中step为人工鱼移动步长,Random为(0~1)中随机一个数;Tmax为最大迭代次数,wij(m)、wij(m+1)分别为人工鱼当前状态和进行移动后的下一步状态;
4)状态为Rm的人工鱼视野中,视野内人工鱼集合可表示为:
Sc={Rm|dm,x≤Visual} (18)
式(18)中Visual代表着人工鱼群的视野,即为权值的变换集合,Sc为人工鱼集合,即网络集合;若Sc≠φ,则在集合内部寻找在中心处的人工鱼,搜寻公式为:
式(19)中Lgi、Vkg、wi0、wg0、wk0为鱼群中初始个体移动权值;
5)视野中鱼群扫描完全后,根据下式计算得视野中人工鱼中心参数:
式(20)中f为视野内搜索得人工鱼个数,执行聚群移动策略,移动方式如(21)式所示:
若人工鱼Lgi、Vkg、wi0、wg0、wk0执行前述觅食行为,移动方式与wij一致;
6)人工鱼视野中寻找Xmax,之后执行追尾移动策略,如(22)式所示
式中,Wij(max)为视野浓度最大处个体,其它层移动方式与wij保持一致。
步骤3.3:对AFSA-Elman边坡位移预测网络中的迭代步长进行修正,所述步骤3.3具体为:
防止迭代过程中跳过最优解,将迭代步长进行修改,公式如式(23)所示:
式中Tmax为最大迭代次数、t为迭代次数,设置最大迭代次数为200,Stepmax、Stepmean为最大迭代步长和平均迭代步长。
以上所述仅为本发明的较佳实施例,并不用以限制本发明的思想,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种岩质边坡结构面自动分组及位移预测方法,其特征在于,包括如下步骤:
步骤1:通过无人机获取岩质边坡的不同角度的多张图像,对不同角度的多张图像进行集成,构建岩质边坡三维模型;
步骤2:对无人机采集获得的边坡三维模型的结构面信息进行分组;
步骤3:根据边坡结构面的分组和无人机采集的数据对岩质边坡位移进行预测。
2.如权利要求1所述的岩质边坡结构面自动分组及位移预测方法,其特征在于,所述步骤1具体包括:
步骤1.1:地面进行控制点选取,通过无人机对待测岩质边坡进行多角度拍摄;
步骤1.2:对岩质边坡的不同角度的多张图像进行集成,进行超高密度点云图像构建,以生成点云模型;
步骤1.3:对点云模型进行纹理修正、精细化处理和去噪处理,获得岩质边坡三维模型;
步骤1.4:获得边坡结构面点云坐标和出露迹长,对结构面产状进行拟合。
3.如权利要求1所述的岩质边坡结构面自动分组及位移预测方法,其特征在于,所述步骤2具体包括:
步骤2.1:结构面视作无厚度的无限延伸的平面,并用结构面法向量来表示;
步骤2.2:将人工鱼群算法与K-Means算法结合,建立AFSA-RSK结构面分类算法,利用鱼群算法对结构面产状聚心集进行搜寻,并通过K-Means算法进行聚类;
步骤2.3:聚类完成后,对人工鱼群算法的步长和视野进行修正,动态调整聚类过程;
步骤2.4:修正视野与步长后,每执行一步与前步结果进行比较,并将两者间误差较小的值记录下来,达到最大迭代次数后,输出记录的最优结果,并与K-Means算法计算结果比较,若优于K-Means算法计算结果,则输出结果;否则,再次修改步长,重新迭代,超过最大迭代次数后输出结果。
4.如权利要求3所述的岩质边坡结构面自动分组及位移预测方法,其特征在于,所述步骤2.1具体为:
用倾角α(0≤α≤90°)和倾向β(0≤β≤360°)表征结构面的空间位置,建立XYZ坐标系,X轴代表正东方向,Y轴代表正北方向,Z轴代表极射球面的正上方,球面半径为1,(n1、n2、n3)为结构面法向量投影在X、Y、Z轴的坐标大小,(n1、n2、n3)为垂直于结构面的法向量,将(n1、n2、n3)的坐标值等效在球面坐标系中,按照投影原理,倾角为α,倾向为β的结构面可由垂直于该结构面的单位法向量n表示,法向量n垂直于结构面,投影在坐标面上表达式为:
则结构面法向量可表示:
n=(n1,n2,n3)=(sinαsinβ,sinαcosβ,cosα) (2)。
5.如权利要求3所述的岩质边坡结构面自动分组及位移预测方法,其特征在于,所述步骤2.2具体为:
1)人工鱼在寻优时,每次移动之后通过评估适应度来确定下一步的行为,当适应度越大,结构面样本中各结构面参数离聚心的距离也就越小,聚类误差也会越小,适应度计算公式为:
式中:k为聚心个数,l为结构面样本个数,||d(ci-xj)||2为分组结构面中样本x离聚心c的距离,距离计算公式为:
式中:nc、nj分别为聚心结构面的单位法向量和j结构面的单位法向量;
2)初始人工鱼在Xi位置处时在视野Visual中随机移动到Xj处,在假定初始聚心的情况下,人工鱼在结构面样本机视野内随机搜寻,搜寻公式为:
Xj=Xi+Visual×Rand() (5)
按照上式随机移动后,到达新的位置Xj处,若Xj处结果适应度更大,则向前一步,得到前进公式为:
Xt+1 i为多次执行觅食行为过程中新聚心的位置;
3)人工鱼在Xi位置处搜索视野内的伙伴数量,计算当前视野内鱼群数目的均值坐标,即为中心坐标值,若中心位置适应度较大,输出中心坐标值,同时按式(7)进行移动,得到新的聚心坐标,并与中心坐标值进行比较,较优的记录下来,中心位置计算方法如式(8);
式中,Xj为视野内某结构面的坐标,g为伙伴数目,Xcen为鱼群经过聚群得到的结构面新聚心;
4)人工鱼Xi搜索得到适应度最大伙伴Xmax,移动方式如式(9),按照式(7)可以得到新的聚类中心Xnext的坐标值,将Xnext和Xmax中的较优者记录下来:
式中,Xmax为视野中最大适应度的个体。
6.如权利要求3所述的岩质边坡结构面自动分组及位移预测方法,其特征在于,所述步骤2.3具体为:
1)采用Xie-Beni指标SXB、分类熵指标Hm两项指标检验聚类效果,若满足ASFA-RSK算法的两项指标误差均小于K-Means算法的两项指标误差,将结果记录,并继续迭代,超过最大迭代次数输出聚类结果;
2)若不满足就对步长和视野进行修正,步长修正公式和视野修正公式如式(10),式(12)所示:
Stepmax选为0.16,则相应的改进觅食行为策略为:
随着迭代步长的增加,步长会随着式(11)变化,从而避免迭代过程中跳过最优值;
视野修正公式:
式中:Vmax选为1,Tmax为最大迭代次数。
7.如权利要求1所述的岩质边坡结构面自动分组及位移预测方法,其特征在于,所述步骤3具体为:
步骤3.1:采用无人机对边坡位移进行实时监测,获得一段时间内的边坡位移序列;
步骤3.2:将人工鱼群算法与Elman网络相结合,建立改进的AFSA-Elman边坡位移预测网络,利用改进的AFSA-Elman边坡位移预测网络对边坡位移序列进行预测;
步骤3.3:对AFSA-Elman边坡位移预测网络中的迭代步长进行修正。
8.如权利要求7所述的岩质边坡结构面自动分组及位移预测方法,其特征在于,所述步骤3.1具体为:
对位移方向进行预判,将坡体点云坐标大致移动方向进行统计,用箭头标出;连接测点位置,对应得到等高线位置上两测点连线的点云坐标变化情况,继而得到坡体连线上各点位移点云变化值,位移计算步骤如下:
1)根据测点进行分组,N个测点对应N个分组,监测标定盘的直径为d,计算标定盘附近点的坐标S,若点K与间距S之间间距小于d,则S、K同为一组;
2)对其他组数据采取与1)相同的计算方式,计算完毕得到N组点云数据;
3)对经建模得到的点云数据进行均值计算,得到中心点的坐标,令中心坐标为Scen(X,Y,Z),并作为该分组的最终坐标值,以中心坐标Scen(X,Y,Z)的变化代表该测点的位置变化,中心坐标Scen(X,Y,Z)计算方法如下式所示:
则每一次测量间隔内,位移变化S为:
其中,t代表测量次数,两次测量的位移变化值为S(t)。
9.如权利要求8所述的岩质边坡结构面自动分组及位移预测方法,其特征在于,所述步骤3.2具体为:
1)利用人工鱼群算法对Elman网络进行改造,鱼群算法寻找的是群体中最优的个体,即食物浓度最大的位置,即误差最小的位置,令网络适应度函数为:
式中:G代表权值变化中适应度最大的位置,Q为网络期望输出与实际输出的误差;
2)选取四层Elman网络,令r为输入层神经元个数,m为隐含层神经元个数,t为承接层神经元个数,任意选取两条人工鱼Rm和Rx,考虑权值变化,利用各层神经元来计算各层连接之间输入与输出的欧氏距离,任意两条人工鱼Rm和Rx之间的距离为:
式(16)中wij为i层隐含层神经元和j层输入层神经元权值,Lgi为g层承接层神经元和i层隐含层权值,Vkg为输入层k层神经元和g层承接层神经元权值,wi0为隐含层神经元阈值,wg0为承接层神经元阈值,wk0为输出层神经元阈值;
3)对网络进行训练,随机指定某条人工鱼状态为Rm,在鱼群视野范围内随机寻找另一个人工鱼,改变权值变化,令其状态为Rx,若该处食物浓度较大,误差计算,即Gx>Gm,则执行觅食行为,反之,向任意方向进行随机移动,网络权值和阈值被改变,移动公式为:
Lgi、Vkg、wi0、wg0、wk0觅食行为与上述wij方式一致,式(17)中step为人工鱼移动步长,Random为(0~1)中随机一个数;Tmax为最大迭代次数,wij(m)、wij(m+1)分别为人工鱼当前状态和进行移动后的下一步状态;
4)状态为Rm的人工鱼视野中,视野内人工鱼集合可表示为:
Sc={Rmdm,x≤Visual}(18)
式(18)中Visual代表着人工鱼群的视野,即为权值的变换集合,Sc为人工鱼集合,即网络集合;若Sc≠φ,则在集合内部寻找在中心处的人工鱼,搜寻公式为:
5)视野中鱼群扫描完全后,根据下式计算得视野中人工鱼中心参数:
式(20)中f为视野内搜索得人工鱼个数,执行聚群移动策略,移动方式如(21)式所示:
若人工鱼的Lgi、Vkg、wi0、wg0、wk0执行前述觅食行为,移动方式与wij一致;
6)人工鱼视野中寻找Xmax,之后执行追尾移动策略,如(22)式所示
式中,Wij(max)为视野浓度最大处个体,其它层移动方式与wij保持一致。
10.如权利要求8所述的岩质边坡结构面自动分组及位移预测方法,其特征在于,所述步骤3.3具体为:
防止迭代过程中跳过最优解,将迭代步长进行修改,公式如式(23)所示:
式中Tmax为最大迭代次数、t为迭代次数,设置最大迭代次数为200,Stepmax、Stepmean为最大迭代步长和平均迭代步长。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910560846.3A CN110288700A (zh) | 2019-06-26 | 2019-06-26 | 一种岩质边坡结构面自动分组及位移预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910560846.3A CN110288700A (zh) | 2019-06-26 | 2019-06-26 | 一种岩质边坡结构面自动分组及位移预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110288700A true CN110288700A (zh) | 2019-09-27 |
Family
ID=68006141
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910560846.3A Pending CN110288700A (zh) | 2019-06-26 | 2019-06-26 | 一种岩质边坡结构面自动分组及位移预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110288700A (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110866531A (zh) * | 2019-10-15 | 2020-03-06 | 深圳新视达视讯工程有限公司 | 一种基于三维建模的建筑物特征提取方法、***及存储介质 |
CN110929599A (zh) * | 2019-11-08 | 2020-03-27 | 宁波大学 | 基于小波分析的岩体结构面轮廓曲线形态分离方法 |
CN111539568A (zh) * | 2020-04-22 | 2020-08-14 | 深圳市地质局 | 基于无人机和三维建模技术的安全监测***及方法 |
CN111985507A (zh) * | 2020-08-28 | 2020-11-24 | 东北大学 | 一种岩体三维点云节理迹线提取方法 |
CN113111595A (zh) * | 2021-06-15 | 2021-07-13 | 西南交通大学 | 一种边坡智能预警方法和装置 |
CN113298939A (zh) * | 2021-04-08 | 2021-08-24 | 湖南科技大学 | 露天矿山边坡岩体结构面抗剪强度分级确定***及方法 |
CN113362459A (zh) * | 2021-04-02 | 2021-09-07 | 成都理工大学 | 一种岩质滑坡破坏边界的特征提取方法 |
CN114353876A (zh) * | 2022-01-07 | 2022-04-15 | 兰州大学 | 一种黄土公路边坡健康监测方法 |
CN115169243A (zh) * | 2022-07-28 | 2022-10-11 | 中铁三局集团有限公司 | 基于ga-pso-glssvm算法的土岩复合地层深基坑变形时序预测方法 |
CN115457229A (zh) * | 2022-11-11 | 2022-12-09 | 中国科学院地质与地球物理研究所 | 岩体结构面三维网络模型确定方法、装置、设备及介质 |
CN117113124A (zh) * | 2023-05-09 | 2023-11-24 | 成都理工大学 | 不连续面优势产状分组方法、***、电子设备以及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108489402A (zh) * | 2018-06-08 | 2018-09-04 | 绍兴文理学院 | 基于三维激光扫描的露天矿山边坡岩体节理规模快速精细取值方法 |
CN109357617A (zh) * | 2018-10-25 | 2019-02-19 | 东北大学 | 一种基于无人机的高陡岩质边坡位移变形监测方法 |
-
2019
- 2019-06-26 CN CN201910560846.3A patent/CN110288700A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108489402A (zh) * | 2018-06-08 | 2018-09-04 | 绍兴文理学院 | 基于三维激光扫描的露天矿山边坡岩体节理规模快速精细取值方法 |
CN109357617A (zh) * | 2018-10-25 | 2019-02-19 | 东北大学 | 一种基于无人机的高陡岩质边坡位移变形监测方法 |
Non-Patent Citations (2)
Title |
---|
王述红等: "《一种改进AFSA - Elman 边坡位移预测网络》", 《东北大学学报(自然科学版)》 * |
王述红等: "《一种改进鱼群聚类算法在结构面分组中的应用》", <东北大学学报(自然科学版)> * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110866531A (zh) * | 2019-10-15 | 2020-03-06 | 深圳新视达视讯工程有限公司 | 一种基于三维建模的建筑物特征提取方法、***及存储介质 |
CN110929599B (zh) * | 2019-11-08 | 2023-05-26 | 宁波大学 | 基于小波分析的岩体结构面轮廓曲线形态分离方法 |
CN110929599A (zh) * | 2019-11-08 | 2020-03-27 | 宁波大学 | 基于小波分析的岩体结构面轮廓曲线形态分离方法 |
CN111539568A (zh) * | 2020-04-22 | 2020-08-14 | 深圳市地质局 | 基于无人机和三维建模技术的安全监测***及方法 |
CN111985507A (zh) * | 2020-08-28 | 2020-11-24 | 东北大学 | 一种岩体三维点云节理迹线提取方法 |
CN111985507B (zh) * | 2020-08-28 | 2023-07-28 | 东北大学 | 一种岩体三维点云节理迹线提取方法 |
CN113362459A (zh) * | 2021-04-02 | 2021-09-07 | 成都理工大学 | 一种岩质滑坡破坏边界的特征提取方法 |
CN113298939A (zh) * | 2021-04-08 | 2021-08-24 | 湖南科技大学 | 露天矿山边坡岩体结构面抗剪强度分级确定***及方法 |
CN113111595B (zh) * | 2021-06-15 | 2021-08-31 | 西南交通大学 | 一种边坡智能预警方法和装置 |
CN113111595A (zh) * | 2021-06-15 | 2021-07-13 | 西南交通大学 | 一种边坡智能预警方法和装置 |
CN114353876A (zh) * | 2022-01-07 | 2022-04-15 | 兰州大学 | 一种黄土公路边坡健康监测方法 |
CN115169243A (zh) * | 2022-07-28 | 2022-10-11 | 中铁三局集团有限公司 | 基于ga-pso-glssvm算法的土岩复合地层深基坑变形时序预测方法 |
CN115457229A (zh) * | 2022-11-11 | 2022-12-09 | 中国科学院地质与地球物理研究所 | 岩体结构面三维网络模型确定方法、装置、设备及介质 |
CN117113124A (zh) * | 2023-05-09 | 2023-11-24 | 成都理工大学 | 不连续面优势产状分组方法、***、电子设备以及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110288700A (zh) | 一种岩质边坡结构面自动分组及位移预测方法 | |
Chen et al. | Machine learning-based classification of rock discontinuity trace: SMOTE oversampling integrated with GBT ensemble learning | |
CN106570468B (zh) | 一种重建LiDAR原始点云建筑物轮廓线的方法 | |
WO2023029461A1 (zh) | 海量高维ais轨迹数据聚类方法 | |
Raumonen et al. | Massive-scale tree modelling from TLS data | |
CN106441292B (zh) | 一种基于众包imu惯导数据的楼宇室内平面图建立方法 | |
Zhang et al. | Automated method for extracting and analysing the rock discontinuities from point clouds based on digital surface model of rock mass | |
Yurtseven et al. | Determination and accuracy analysis of individual tree crown parameters using UAV based imagery and OBIA techniques | |
CN108489402A (zh) | 基于三维激光扫描的露天矿山边坡岩体节理规模快速精细取值方法 | |
CN108489403A (zh) | 基于三维激光扫描的露天矿山边坡岩体节理产状快速精细取值方法 | |
CN110363299B (zh) | 面向露头岩层分层的空间案例推理方法 | |
CN108061901A (zh) | 基于机载激光雷达点云数据重建3d电力线模型的方法 | |
CN109299739A (zh) | 基于法向量的曲面拟合对车载激光点云滤波的方法 | |
CN110443810A (zh) | 基于快速邻接体素查询的点云平面分割方法 | |
CN112068153B (zh) | 一种基于地基激光雷达点云的冠层间隙率估算方法 | |
Wang et al. | UAV photogrammetry and AFSA-Elman neural network in slopes displacement monitoring and forecasting | |
CN107247926A (zh) | 一种人体检测方法及装置 | |
CN108830317B (zh) | 基于数字摄影测量的露天矿山边坡岩体节理产状快速精细取值方法 | |
CN109191519A (zh) | 一种树木树干体积估算模型构建方法、体积估算方法及*** | |
CN114387506A (zh) | 输电杆塔监测方法、装置、计算机设备、存储介质 | |
Lalonde et al. | Automatic three-dimensional point cloud processing for forest inventory | |
CN107705344A (zh) | 激光扫描环境点云数据中树冠模型提取方法 | |
Xiao et al. | Filtering method of rock points based on BP neural network and principal component analysis | |
CN117197677A (zh) | 一种基于激光雷达点云数据的热带雨林乔灌分离方法 | |
Nguyen et al. | Grouping discontinuities in representative sets: influence on the stability analysis of slope cuts |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190927 |
|
WD01 | Invention patent application deemed withdrawn after publication |