摘要:为提高地下开采引起地表下沉预测结果的精度,提出融合混沌残差的BP强预测器(BP-Adaboost)的地表下沉预测模型。以顾北矿1312(1)实测值为例,分别用融合混沌残差的BP-Adaboost模型、BP神经网络模型和BP-Adaboost模型对最大下沉值点进行稳定期和活跃期的单步预测和多步预测,结果表明,融合混沌残差的BP-Adaboost模型无论是在单步预测还是在多步预测上的精度均最高,尤其在单步预测上有显著的提高。
加入收藏
构建地表下沉预测模型是矿山开采的难点,但对于由开采沉陷引发的地质环境灾害的防治具有重要意义[1]。目前,针对沉降预测的方法有模糊神经网络[2]、灰色模型[3]、支持向量机[4]、径向基神经网络[5]、组合预报模型[6]、多核相关向量机[7,8]、时间序列分析[9]和智能算法优化的BP神经网络[10]等。但利用上述方法进行建模时,模糊神经网络的隶属函数带有一定的主观性,灰色模型精度较低,支持向量机的核函数的参数较难确定,径向基神经网络容易陷入局部最优解(结果不稳定),组合预报模型权系数难以确定,多核相关向量机的核函数难以选择,时间序列分析精度较低,智能算法优化BP神经网络时间较长(迭代次数不易确定)。针对上述问题,BP-Adaboost在寻找全局最优解方面有一定的优势[11]。因为地下开采引起的地表下沉具有显著的混沌特性和非线性特性,本文针对BP-Adaboost没有考虑残差而导致精度降低并且不稳定的问题,运用混沌理论对下沉残差序列进行分析和预测,提出一种新的融合混沌残差的BP-Adaboost预测模型。
1、融合混沌残差的BP-Adaboost模型预测原理
BP神经网络和BP-Adaboost神经网络原理见文献[10,11]。
1.1混沌残差序列的预测
客观事物的运动除了周期、准周期和定常以外,还存在一种更具普遍意义的运动形式,即混沌。混沌具有如下特征:内部似随机性、整体稳定局部不稳定性、初值敏感性、短期可预测长期不可预测性、奇异吸引子、轨道稳定性及分岔、普适性等。混沌残差序列的预测分为3个步骤,分别是残差序列的相空间重构、混沌特性的识别、混沌残差的预测。
1.1.1残差序列的相空间重构
相空间重构理论是混沌序列的基础,可用坐标延迟法对混沌序列进行相空间重构。对于残差序列{d1,d2,d3,…,dn}可以重构一个m维的相空间:
[X1X2⋯XΚ]=[d1d2⋯dΚd1+τd2+τ⋯dΚ+τd1+2τd2+2τ⋯dΚ+2τ⋮⋮⋯d1+(m-1)τd2+(m-1)τ⋯dΚ+(m-1)τ] (1)
式中,m为嵌入维数,τ为坐标延迟时间,K为相点个数,X为相点,K=n-(m-1)τ。运用G-P算法可以获取混沌序列的关联维(G),而如果嵌入维数m≥2G+1,则系统原始状态变量构成的相空间和一维观测值重构相空间的动力学行为等价,即获得了所有状态变量演变的全部信息,从而可获得下一时刻的预测值,同时嵌入维数可以通过G-P算法获得[12]。残差序列的τ取1。
1.1.2残差序列混沌特性的识别
常见判断残差序列混沌特性的方法有最大Lyapunov指数法和关联维分数判别法,下面主要介绍关联维分数判别法。关联维数为分数则说明该系统的动力行为可能具有混沌特征,只要相空间重构的维数足够大,就能描述出系统的混沌吸引子,表示出传统方法无法展示的系统运动规律。但是在实际测算中也可能出现关联维没有饱和的情况,其原因有两种可能:1)m不够大,所以应继续增大m;2)该序列没有无特征尺度区[13]。
1.1.3一阶加权局域预测法
常用的混沌时间序列的预测方法有零阶局域预测法、加权零阶局域预测法、一阶局域预测法、加权一阶局域预测法、全局预测方法和最大Lyapunov指数法等,下面主要介绍一阶加权局域预测法。局域法是将相空间轨迹的最后一点作为中心点,把离中心点最近的若干轨迹点作为相关点,然后对这些相关点进行拟合,再估计轨迹下一点的走向,最后从预测出的轨迹点的坐标中分离出所需要的预测值,其原理即“寻找历史上情况最相似之处”。对于状态空间中重构的状态矢量,假设中心点为Xm,Xmi(i=1,2,…,q,q是邻近点的个数)到Xm的距离是disi,其中的最小值是dismin,Xmi的权值表示为:
pdi=exp(-c(disi-dismin))∑i=1qexp(-c(disi-dismin)) (2)
式中,c为参数,一般取c=1。一阶加权的拟合表达式为:
Xmi+1=αe+bXmi (3)
式中,a、b是拟合系数,e=(1,1,…,1)T,Xmi+1是通过Xmi预测的值。
通过加权最小二乘求取拟合系数a、b,其表达式为:
∑i=1qpdi(Xmi+1-αe+bXmi)2=min (4)
求解出a、b的值,代入(3)式即可求出混沌序列的预测值。
1.2融合混沌残差的BP-Adaboost模型预测原理
将BP-Adaboost神经网络模型的残差d代入式(1),再由式(4)得出拟合系数a、b,代入式(3)得到残差,最终将每个预测对应的残差和BP-Adaboost神经网络模型的预测值pre累加就得到了融合混沌残差的预测值。融合混沌残差的BP-Adaboost预测算法流程如图1所示,图中Etrain为训练误差矩阵。
图1融合混沌残差的BP-Adaboost预测算法流程
融合混沌残差的BP-Adaboost预测算法的表达式为:
y=pre+ε (5)
式中,ε是残差预测值。
2、工程实例
2.1工程应用
以淮南顾北煤矿1312(1)工作面的最大下沉点为研究对象,经过数据预处理,以下沉速度1.6mm/d为临界速度,将下沉点的下沉进程分为开始期、活跃期和稳定期。开始期和活跃期以4d为一个三次样条插值,稳定期以10d为一个三次样条插值,最终得到开始期的天数是1~7d,活跃期的天数是8~108d,稳定期的天数是109~409d。获得的最大下沉点MS23的图像如图2所示。
图2最大下沉点MS23下沉图
2.2最大下沉点活跃期的预测
以最大相对中误差(error1)、平均相对中误差(error2)、均方根误差(error3)和最大误差绝对值(error4)为指标,对BP神经网络模型、BP-Adaboost神经网络模型和融合混沌残差的BP-Adaboost模型进行精度评价。其中网络的输入为活跃期前7个三次样条插值下沉数据,输出为活跃期输入值下一个三次样条插值的下沉数据,隐含层的节点个数是8,网络的拓扑结构是7-8-1,预测最后5个三次样条插值。重复计算10次后,下沉平均值的模型精度比较如表1所示。
表1最大下沉点的活跃期单步预测精度比较
从表1可以看出,融合混沌残差的BP-Adaboost模型在4种精度评判指标中精度均最好,其次是BP-Adaboost神经网络模型,BP神经网络最差。
BP-Adaboost模型拟合下沉值残差的嵌入维和相关维的值如表2所示。从表3可以看出,在进行多步预测的时候,3种模型的预测精度均有所下降,BP神经网络模型精度最低,然后是BP-Adaboost神经网络模型,融合混沌残差的BP-Adaboost模型精度最高。融合混沌残差的BP-Adaboost模型在BP-Adaboost神经网络模型的基础上略微有所提高,两者精度相差不大的原因在于,混沌序列只能进行短期预测,不能进行长期预测,同时又没有顾及数据的新鲜度。
表2BP-Adaboost模型拟合下沉值残差的嵌入维和相关维
表3最大下沉点的活跃期多步预测精度比较
2.3最大下沉点稳定期的预测
用3种模型进行最大观测点稳定期的预测,同样用上节中4个指标进行精度评价。其中网络的输入为稳定期前7个插值的下沉数据,输出为下一个三次样条插值的下沉数据,隐含层的节点个数是8,网络的拓扑结构是7-8-1,预测后5个三次样条插值的下沉值。重复计算10次后,下沉平均值的模型精度比较如表4所示。
表4最大下沉点的稳定期单步预测精度比较
从表4可以看出,在最大下沉点的稳定期单步预测时,3种模型的预测精度均很高,原因在于稳定期下沉值变化相比于活跃期不是很大。其中BP神经网络模型和BP-Adaboost神经网络模型精度差不多,融合混沌残差的BP-Adaboost模型效果最好,具有明显的优势。
从表5可以看出,在最大下沉点的稳定期多步预测时,3种模型的预测精度均很高,其中BP神经网络模型最差,BP-Adaboost神经网络模型和融合混沌残差的BP-Adaboost模型结果精度相差不大,再次说明混沌序列只适合短期预测,由于没有顾及下沉数据的新鲜度,所以精度只是略有提高。
表5最大下沉点的稳定期多步预测精度比较
最大下沉点的稳定期多步BP-Adaboost预测残差如表6所示。
表6最大下沉点的稳定期多步BP-Adaboost预测残差
3、结语
为提高地下开采引起地表下沉预测结果的精度,本文提出融合混沌残差的BP-Adaboost地表下沉预测模型。通过对顾北矿1312(1)实测最大下沉值的预测发现,融合混沌残差的BP-Adaboost地表下沉预测模型单步预测和多步预测精度均高于BP神经网络模型和BP-Adaboost神经网络模型,尤其在对单步预测的时候,精度提高最为明显。
参考文献:
[1]李培现.深部开采地表沉陷规律及预测方法研究[D].徐州:中国矿业大学,2012.
[3]何君,杨国东.灰色预测理论在建筑物沉降中的应用研究[J].测绘通报,2012(3):63-64.
[5]谭洋,牛雪峰,王明常.径向基神经网络在沉降预测中的应用[J].测绘科学,2016,41(4):33-36.
[6]魏博文,熊威,李火坤,等.融合混沌残差的大坝位移蛙跳式组合预报模型[J].武汉大学学报:信息科学版,2016,41(9):1272-1278.
[8]罗亦泳,姚宜斌,赵庆志,等.利用优化的组合核相关向量机算法构建地表下沉预测模型[J].武汉大学学报:信息科学版,2018,43(9):1295-1301.
[9]王晶晶,尹晖.一种建筑沉降叠加预测方法[J].测绘科学,2019,44(3):107-113.
[10]李苋兰,张顶,黄晞.基于BP-AdaBoost神经网络的多参数掌静脉图像质量评价法[J].计算机系统应用,2020,29(3):20-28.
[11]刘艳丽,陈跃东.基于BP-Adaboost的目标跟踪算法应用研究[J].河北工程大学学报:自然科学版,2012,29(3):99-102.
[12]董春娇,邵春福,张辉,等.基于G-P算法的快速路交通流参数相空间重构[J].吉林大学学报:工学版,2012,42(3):594-599
[13]张安兵.动态变形监测数据混沌特性分析及预测模型研究[D].徐州:中国矿业大学2009.
陈兴达,余学祥,池深深,蒋创,赵祥硕.融合混沌残差的BP强预测器的地表下沉预测模型[J].大地测量与地球动力学,2020,40(09):913-917.
基金:淮南矿业(集团)有限责任公司基金(HZMDGB-JF2013-14);淮浙煤电有限责任公司基金(HZMDGB-JF2019-0501).
分享:
随着社会经济不断发展,社会化进程不断加快,国际油田勘探工程也得到了较快的发展,全球工业生产,社会经济的发展对于石油的需求量越来越大,石油产业发展的相关技术研究也得到了业界的关注。目前,国际油田开采存在受到地质环境限制,石油勘探及开采项目开展会受到地质环境的限制,增加了石油勘探的难度。
2024-08-31目前我国地下金属矿山的开采正处于“由浅至深、由易至难、由富至贫”的关键转型期,在一个矿区内采用单一采矿方法进行回采的传统模式已难以实现矿产资源安全、经济、高效、充分回收的目标。近年来,协同开采理念的提出为化解这一难题提供了一条重要途径该理念的核心在于将不同采矿方法在同一矿区内进行最大程度整合。
2024-06-25随着全球化进程的加快和社会的不断发展,工业、军事等领域对矿产资源需求不断增大,但我国主要资源开采带的矿产储量逐渐减少,亟需对我国复杂地形的矿产资源进行勘探开采,因此需要更高性能和更加灵活的地层勘探仪器来探明地层结构,从而保障资源开采过程中人员和设备的安全。
2024-01-03随着全球化进程的加快和社会的不断发展,工业、军事等领域对矿产资源需求不断增大,但我国主要资源开采带的矿产储量逐渐减少,亟需对我国复杂地形的矿产资源进行勘探开采,因此需要更高性能和更加灵活的地层勘探仪器来探明地层结构,从而保障资源开采过程中人员和设备的安全。
2024-01-03油气管道是目前最重要的能源输送方式之一,我国陆上油气管道总里程长、分布广阔、沿线地质环境复杂多变。由于各方面的考虑,管道可能敷设在滑坡地带,也有可能在管道敷设后由于工程活动或者地址变动等种种原因导致沿线发生滑坡。而滑坡灾害对管道危害巨大,能够挤压管道导致管道变形甚至破坏,进而导致管道输送的油气泄露,
2023-12-29由于我国采矿技术的提升,以及煤矿开采范围的扩大,部分矿业公司在经营管理与生产活动中将面临着很多的安全难题,比如对安全技术的了解不足、对安全方面的资金支持欠缺等,使采矿过程中地质灾害问题频频出现。不论是矿业企业的管理者还是其他管理人员,都应当提升采矿过程中的安全生产能力,积极采取正确可行的安全措施,减少矿井地质灾害的风险,维护矿井职工的安全,确保我国矿业的可持续、安全、平稳的发展。
2023-12-06据中国地震台网正式测定,当地时间2023-02-06,土耳其发生4次地震,其中2次震级高达7.8。研究表明,地震前后大气中的水汽含量会出现异常变化[1,2,3,4,5]。大气可降水总量(total precipitable water, TPW)是表征大气中水汽含量的重要指标,而卫星遥感是反演TPW的重要手段之一[6,7]。
2023-11-29天水市是典型的山间河谷盆地型城市。区内黄土分布广,厚度大,垂直节理和裂隙发育,降雨集中分布,致使沟壑纵横,地貌地形支离破碎,地质环境十分脆弱[1,2,3],在自然因素以及人类工程活动影响下,滑坡、泥石流、崩塌、不稳定斜坡失稳等地质灾害现象频繁发生[4]。
2023-11-08在工程建设和岩土施工过程中,地质灾害尤其是滑坡问题时有发生,这主要是由于边坡的稳定性受到了岩土层变化的影响,因此,提高边坡的稳定性对预防滑坡问题至关重要。随着科学技术的不断发展,我国岩土施工技术也在不断地改进中,岩土施工技术的合理运用,对有效治理地质灾害起决定性作用,要求相关技术人员必须认识到科学合理应用其施工技术的重要意义。
2023-11-08乐矿能源(原乐平矿务局,简称乐矿)地处赣东北景德镇市境内,其下属的煤矿均分布在萍乐坳陷东北端,多在20世纪50、60年代建设。乐平矿区所属煤矿煤炭资源分布不均,开采条件相对较差,矿井规模都比较小,而且煤层瓦斯含量高、灾害严重,矿区矿井水文地质条件也较复杂。
2023-11-08人气:6027
人气:3748
人气:3695
人气:3106
人气:2962
我要评论
期刊名称:地球科学进展
期刊人气:1864
主管单位:中国科学院
主办单位:中国科学院资源环境科学信息中心,国家自然科学基金委员会地球科学部,中国科学院前沿科学与教育局
出版地方:甘肃
专业分类:科学
国际刊号:1001-8166
国内刊号:62-1091/P
邮发代号: 54-86
创刊时间:1986年
发行周期:月刊
期刊开本:大16开
见刊时间:一年半以上
影响因子:2.546
影响因子:0.659
影响因子:0.149
影响因子:1.209
影响因子:2.866
您的论文已提交,我们会尽快联系您,请耐心等待!
你的密码已发送到您的邮箱,请查看!