第 42 卷第 3 期
2017 年 3 月
煤
炭
学
报
JOURNAL OF CHINA COAL SOCIETY
Vol.42 No.3
Mar.
2017
温廷新,陈晓宇,邵良杉,等.参数优化 GA-ELM 模型在露天煤矿抛掷爆破的预测[J]. 煤炭学报,2017,42( 3) :630- 638. doi:10.
13225 / j.cnki.jccs.2016.0572
Wen Tingxin,Chen Xiaoyu,Shao Liangshan,et al. Prediction on parameters optimized GA-ELM model for cast blasting in open-pit mine
[J].Journal of China Coal Society,2017,42(3) :630-638.doi:10.13225 / j.cnki.jccs.2016.0572
参数优化 GA-ELM 模型在露天煤矿抛掷爆破的预测
温廷新1
2,陈晓宇1,邵良杉1,窦 融2,魏 鹏2
,
(1. 辽宁工程技术大学 系统工程研究所,辽宁 葫芦岛 125105;2. 辽宁工程技术大学 工商管理学院,辽宁 葫芦岛 125105)
摘 要: 为有效指导露天煤矿制定正确生产计划,提高露天煤矿抛掷爆破预测的准确率,在分析露
天煤矿抛掷爆破影响因素的基础上,通过“试错法”确定模型最优隐含层节点参数,进而提出一种
参数优化后遗传算法( GA) 和极限学习机( ELM) 相结合的抛掷爆破预测模型。选取网络的输入输
出相关参数,针对现有 ELM 输入权值矩阵和隐含层偏差,采用遗传算法对其进行优化选择; 利用某
露天煤矿抛掷爆破监测数据对该模型进行实例分析,并将 RBF,BP,SVM,GA-BP 模型预测结果与
该模型进行对比分析; 并引入 Weibull 模型,通过预测控制参数 ɑ,β 模拟爆堆形态。研究结果表
明: ① 通过“试错法”确定 GA-ELM 模型最优隐含层节点数为 39,有效降低系统的仿真误差,该参
数下仿真误差值为 0. 137 7; ② 相较于传统 ELM 预测模型,通过遗传算法优化后,有效抛掷率,松
散系数以及抛掷距离均得出更小的均方误差 MSE 值( 0. 258 0,1. 748 5×10-4,3. 618 4) 和更高的决
定系数 R2值( 0. 986 4,0. 995 3,0. 970 6) ,改进后的 GA-ELM 具有更好的拟合效果和泛化能力; ③
通过与其他智能算法如 BP,RBF,SVM,GA-BP 相比,改进后的 GA-ELM 测试结果( 均方误差,决定
系数,仿真误差) 明显优于其他预测模型,有效提高预测精度; ④ 利用训练完成的 GA-ELM 网络预
测爆堆形态时,控制参数 a,β 的预测误差均未超过 5%,预测爆堆曲线接近真实爆堆曲线。
关键词: 露天煤矿; 抛掷爆破; 预测; GA-ELM 模型; 试错法
中图分类号: TD824. 3
文章编号: 0253-9993( 2017) 03-0630-09
文献标志码: A
Prediction on parameters optimized GA-ELM model for cast
blasting in open-pit mine
WEN Ting-xin1,2,CHEN Xiao-yu1,SHAO Liang-shan1,DOU Rong2,WEI Peng2
(1.System Engineering Institute,Liaoning Technical University,Huludao
125105,China;2.School of Business Administration,Liaoning Technical University,
125105,China)
Huludao
Abstract:In order to make a correct production plan effectively and improve the accuracy of cast blasting in open-pit
mine,based on the analysis of cast blasting factors,the optimized hidden layer node parameters are determined by“tri-
al-and-error”and then one kind prediction model is proposed with the combination of Genetic Algorithm ( GA) and
Extreme Learning Machine ( ELM) .On the basis of network input and output parameters’selection,the optimization
and selection are made using GA for ELM model’s input weight matrix and hidden layer deviation.Using a cast blas-
ting monitoring data in an open-pit mine,the model analysis was conducted,then the prediction results of RBF,BP,
SVM,and GA-BP models were compared with the model results.In addition,Weibull model was introduced to simulate
the coal blasting form by predicting control parameters a,β. The results show that:① The optimal hidden layer node
收稿日期: 2016-05-04
基金项目: 国家自然科学基金资助项目(71371091) ;辽宁省社科基金资助项目( L14BTJ004)
作者简介: 温廷新(1974—) ,男,山西太谷人,教授,硕士生导师,博士研究生。E-mail:wen_tx@ 163. com。通讯作者: 陈晓宇( 1993—) ,男,硕
修回日期: 2016-09-14
责任编辑: 许书阁
士研究生。E-mail:chenxy93@ qq. com
中国煤炭期刊网 www.chinacaj.net
第 3 期
温廷新等:参数优化 GA-ELM 模型在露天煤矿抛掷爆破的预测
136
number is 39,which reduces the system simulation error effectively,and the simulation error is 0.137 7;② Compared
with traditional ELM prediction model,after optimized by GA,Effective Casting Blasting Ratio,Loose Coefficient and
Cast Distance all have lower Mean Square Error ( MSE:0.258 0,1.748 5×10 -4,3.618 4) and higher Decision Coeffi-
cients ( R2:0.986 4,0.995 3,0.970 6) .Improved GA-ELM model has a better fitting effect and generalization ability;
③ Compared with other intelligent algorithm,such as BP,RBF,SVM,GA-BP,the test set results ( MSE,R2,SE) of
improved GA-ELM model are much better and improve the prediction accuracy effectively;④ Using trained GA-ELM
network to predict coal blasting form,the prediction error of control parameters a and β are not more than 5% and the
prediction of blasting heap curve is close to the real blasting heap.
Key words:open-pit mine;cast blasting;prediction;GA-ELM model;trial-and-error
在露天煤矿开采中应用高台阶抛掷爆破技术,可
以降低开采过程中的剥离成本———在无需进行采,
运,排作业的情况下,直接将 30% ~ 60% 的剥离物抛
掷到采空区,减少拉斗铲剥离量。因此,深度挖掘露
天煤矿抛掷爆破监测数据的时变规律,通过智能方法
准确预测露天煤矿抛掷爆破变化趋势,对于露天煤矿
抛掷爆破具有十分现实的经济意义。
目前,国内有关论述露天煤矿抛掷爆破抛掷的
文献仍为数不多。近年来,大多数学者关于露天煤
矿抛掷 爆 破 预 测 与 模 拟 的 研 究 主 要 集 中 为 两 类。
一种是通过 力 学 模 型 进 行 理 论 分 析 和 现 场 试 验 探
讨露天煤矿抛掷爆破效果:自 2006 年开始,李克民
等采用理论结合经验公式的方法,并结合 D’Appo-
lonia“图解 法 ”最 终 对 抛 掷 爆 破 参 数 作 出 优 化 分
析[1];周伟等利用弹道理论以及露天采矿学的基本
原理,建立了抛掷速度反演模型和有效抛掷率离散
模型[2];李祥龙等通过理论分析和露天煤矿现场实
面的预测分析。此外,由于 BP 神经网络具有易陷入
局部最优、收敛速率慢,训练样本需求大等局限性,用
其对抛掷爆破进行预测所得结果误差较大,并不理
想;而 SVM 模型核函数参数以及参数惩罚项的选择、
优化均对 预 测 结 果 有 较 大 影 响; 鉴 于 此,笔 者 结 合
ELM 算法 训 练 速 度 快,泛 化 能 力 好 等 优 点,采 用
GA( 遗传算法) 优化 ELM( 极限学习机) 输入权值矩
阵 和 隐 含 层 偏 差,对 露 天 煤 矿 抛 掷 爆 破 进 行 预
。
测[9-12]
1 GA-ELM 基本理论
1. 1 ELM 基本原理
极限学习机 ( ELM) 算法是 2004 年南洋理工大
学黄广斌副教授提出来的一种网络学习算法,其本质
上属于典型的的单隐层前馈神经网络。但相较于传
统的单隐层前馈神经网络,ELM 克服了传统梯度算
法易陷入局部极值的局限性,具有学习速度快,分类
验,研究台 阶 高 度 对 抛 掷 爆 破 效 果 的 影 响[3];马 力
准确 率 准 确 率 高,泛 化 能 力 好,调 节 参 数 少 等 优
等在弹道理论和平面药包法的理论基础上对( 炸药
点[13-16]
。
单耗,炮 孔 倾 角) 与 抛 掷 距 离 的 关 系 展 开 深 入 研
究[4];另一种 是 以 抛 掷 爆 破 相 关 数 据 为 基 础,结 合
智能化数学模型展开分析预测:如刘希亮等提出遗
传算法优化 的 支 持 向 量 机 模 型 对 抛 掷 爆 破 有 效 抛
掷率进行 预 测[5];孙 文 希 等 以 抛 掷 爆 破 台 阶 高 度、
孔距等参数的 MIV 值作为评价依据,研究输入参数
对抛掷距离等输出参数的影响权重[6];孙文希等通
过设置不同参数,研究 BP 神经网络的收敛速度和
其预测结果的准确性,并与 RBF,SVM 预 测 模 型 的
结果进行对比[7];刘干等在分析了抛掷爆破有效抛
掷率影响因素的基础上,通过广义回归神经网络建
立有效抛率预测模型[8]
。
在上述学者的研究中,大多以单评价因子建立预
测模型。本文综合考虑有效抛掷率,松散系数,抛掷
距离 3 项重要指标,并通过预测 Weibull 模型控制参
数模拟爆破堆形态,对露天煤矿抛掷爆破进行更为全
ELM 实现的过程,对于 N 个任意的各不相同的
样本( xi
xi = [xi1
,yi
,xi2
) ,其中
,…,xin
]T "Rm
则前馈神经网络激励函数为 g( x) ,隐含层节点
]T "Rn;yi = [yi1
,…,yin
,yi2
数为 L 时,输出可以表示:
L
i = 1
fL
) ,
β i G( ai·xi + bi
( x) = ∑
x ∈ Rn,ai ∈ Rn,β i ∈ Rm
,ai2
(1)
]T 为输入层到隐含层节点 i
,
,βi2
表
]T 为连接隐含层节点 i 的输出权值;ai ·xi
为第 i 个隐含层节点的阈值;β i =[βi1
,…,ain
式中,ai =[ai1
的输入;bi
…,βim
示向量 ai
和 xi
的内积。
G( x ) 为 激 励 函 数,可 以 采 用 ”RBF ”,”Sine ”
或”Sigmoid”等。
如果这个隐含层节点数为 L 的前馈神经网络能
以误差 0 逼近 N 个样本,则 ai
,bi
,β i
存在关系如下:
中国煤炭期刊网 www.chinacaj.net
煤
炭
学
报
2017 年第 42 卷
236
fL
L
( x) = ∑
i = 1
β i G( ai * xi + bi
i = 1,2,…,L
) = yi
,
(2)
式(2) 可简化为:Hβ = Y',式中,
H( a1
,…,aL
g( a1* x1 + b1
,…,bL
,b1
,…,xN
) =
)
) … g( a1* x1 + bL
,x1
…
…
…
g( a1* xN + b1
β =
)
β T
1
) … g( a1* xN + bL
,T =
yT
β T
yT
1
L
L
N ×L
H 为隐含层的输出矩阵。一旦 ELM 算法中输出
权值和阈值被随机确定,隐层的输出矩阵 H 就变成
惟一确定的。而训练前馈神经网络就可以转化成一
个求解线性系统最小二乘解的问题。输出权重 β 也
可以被确定:β^ = H'Y,其中,H'是矩阵 H 的 Moore -
Penrose 广义逆。
传统 ELM 算法的流程如下:
Step 1:随机给定网络输入隐含层权值 αi
,i = 1,…,L:其中隐含层节点个数为 L;
Step 2:选择隐含层神经元的激活函数为一个无
以及阈
值 bi
限可微的函数,并对隐含层输出矩阵 H 计算;
Step 3:计算网络隐含层节点与输出节点之间的
的连接权值 β^ :β^ = H'Y。
1. 2 GA 优化 ELM 预测模型
随机给定 ELM 的输入权值矩阵和隐含层偏差,
可能部分数值为 0,使得一些隐含层节点失效,为达
到预计精度只有提高隐含层节点的数目。但随之带
来模型 对 未 知 样 本 适 应 能 力差,导 致 泛 化 能 力 低
下[17-20]
。针对上述问题,本文采用遗传算法对极限
学习机中输入层权值和隐含层偏差进行优化选择 。
训练步骤如下:
(1) 产生种群。种群个数 N 一般设置为 20 ~ 40,
种群中的个体是由输入权值矩阵和隐含层偏差构成,
个体长度为 L = ( n + 1) k,式中,k 为隐含层节点数;n
为输入层神经元数,即输入向量维数。
(2) 对种群中的任一个体,采用 ELM 算法( 隐含
层激活函数为 Sig 函数) 计算输出权值矩阵,通过代
入 P _train 计算,得出初始化种群每个个体的分类准
确率,将该参数作为遗传优化算法的适应度。计算种
群第 1 个个体 θ1 的适应度值 Fitness_1,将 Fitness_1
赋值给最优适应度 Fitness_best;并将 θ1 作为 pBest;然
后执行以下伪代码:
Step 1:for i = 1:N { 计算 θi的适应度值 Fitness_1
If Fitness _ 1 > Fitness _ best
then
pBest( i)= θi
Else pBest( i)= θi-1}
令 Fitness_i 对应的粒子 θ 为全局极值 gBest
Step 2:for i = 1:P %P 为最大迭代次数
{ 利用交叉,变异操作对当前群体进化;
If 满足约束条件 or
i>P
Break;
Else loop Step 1 }
(3) 执行上述 伪 代 码 流 程,计 算 出 最 优 适 应 度
所对应个体,即 ELM 算法所需的最优输入权值矩阵
和隐含层偏差,利用 ELM 算法继续求得输出权值矩
阵。
2 GA-ELM 预测模型因子选取
2. 1 影响因子选取
露天煤矿抛掷爆破是一种非线性的复杂动力系
统[21],它时刻受地质条件、岩层性质、炸药特性等多
种因素影响而不断发生着变化。目前研究文献中,露
天煤矿抛掷爆破影响因素主要分为岩石性质、炸药特
性、爆破实施 3 类[1,5,22-25]
。
由于露天煤矿地质构造及岩石特性较为单一,且
为提高开采效率大多选择岩性比较均匀的矿区,故在
本次研究中,笔者将岩性条件视为同质。此外,抛掷
爆破实施过程中通常采用相同的炮孔直径、炸药类
型、装药结构、炮孔堵塞方式等,因此,抛掷爆破优化
设计更多的是通过调整台阶高度、孔距、排距等地形
和设计参数来实现。
因此,笔者选取以下 10 个参数作为 GA-ELM 预
测模型的影响因子:炸药单耗、台阶高度、煤层厚度、
采空区上口宽、采空区下口宽、最小抵抗线、坡面角、
孔距、排距、剖面宽,如图 1 所示。
图 1 露天煤矿抛掷爆破影响因子
Fig. 1 Effective factors of cast blasting
2. 2 评价因子选取
有效抛掷率,松散系数,抛掷距离是反映露天矿
抛掷爆破效果的重要评价因子[5,7,24],各参数计算示
意[4]如图 2 所示。
中国煤炭期刊网 www.chinacaj.net
第 3 期
温廷新等:参数优化 GA-ELM 模型在露天煤矿抛掷爆破的预测
336
计算松散系数。
取 H( X) 为 Weibull 分布的概率密度函数,有
图 2 计算示意
Fig. 2 Calculation schemes
有效抛掷率 ηa
指 被 爆 岩 土 体 积 中 通 过 爆 破 技
术不需要 2 次转运直接抛掷到采空区或者坡下废
系数越大,被抛掷
石体积所占的百分比系数。该 ηa
到采空区的岩石量也就越大,拉斗铲的作业量也就
越小,降低了剥离成本。按图 2 所 示,即 为 A 部 分
岩石所占 A,B,C 部 分 岩 石 总 体 积 的 百 分 比,表 示
如下:
ηa = VA / ( VA + VB + VC
) × 100%
松散系数 ηb
指爆破前后松散体积与岩石体积之
系数值影响着爆堆沉降高度、拉斗铲铲装
比。该 ηb
效率及辅助设备作业量等相关变量。按图 2 所示,即
为 A,B,C 部分作业后破碎岩石总体积与 C,D 部分
原石台阶岩体体积的比值,表示为
ηb = ( VA + VB + VC
) / ( VC + VD
)
抛掷距离 ηc
指从爆前自由面到爆后倒堆面之间
的水平距离,实际作业过程中利用爆后自由面数据与
倒堆面数据相交后的叠加面进行采集。该指标是保
证抛掷爆破安全性的重要指标之一。
鉴于此,并综合考虑前人相关研究,笔者选取抛
掷距离、松散系数及抛掷率作为本次 GA -ELM 露天
煤矿抛掷爆破预测模型的评价因子。
2. 3 Weibull 模拟爆堆形态预测参数选取
用 Weibull 分布模拟爆堆形态其难点在于 a,β
值的选取,a,β 是控制 Weibull 分布曲线形状的重要
参数。Weibull 分 布 模 型 理 论 是 1939 年 Weibull 提
出,以质量守恒定律为力学基础,被广泛应用于处理
各种随机现象的数量性特征[26]
。
根据质量守恒定律,爆破前后岩体质量不变,即
lm
∫
0
ρhh( x) dx = S0ρq
(3)
分别为爆破前后岩石密度,kg / m3;h( x) 为
为待爆岩体的台阶剖面
,ρh
式中,ρq
在 x 轴上的爆堆高度,m;S0
面积,m2;lm
为岩块最远抛距,m。
将式(3) 无量纲化,得到
lm
∫
0
H( X) dX = 1
(4)
式中,X =
x
,H =
h
ξS槡 0
ξS槡 0
,Lm =
lm
ξS槡 0
,ξ' =
ρh
ρq
,其中 ξ' 为
H( X) =
0
(
β
α
(
) X
α
) β - 1
(0 ≤ X ≤ Lm
)
( X > Lmor X < 0)
[
exp -
) β
]
(5)
(
X
α
其中,ɑ,β 为控制曲线形状的参数( β>1) ,如果 ɑ,β
选择合理的话,H( X) 在 lm
出变化很小,则式( 4) 可以
H( X) dX = 1。
写成: ∫
+#
-#
多排微差起爆的台阶爆破由于最后一排炮孔的
爆破漏斗作用,在所形成的爆堆中会出现一个漏斗状
的凹陷[27-28]
。真实爆堆形态轮廓线与 Weibull 分布
曲线仅表现为主要部分符合。在利用 Weibull 模型
模拟爆堆形态需要通过无量钢化因子对真实的爆堆
形态曲线无量纲化,2 者位于同一坐标系下;式( 1) 无
量纲化的推导过程中,由 ξ = ρh / ρq
,可得无量纲化因
子为 1
。 根据现场统计资料和经验分析,确定爆
S0槡 ξ
破漏斗深度为(0. 17 ~ 0. 38) H',H'为台阶高度。
3 实证分析
经 整 理,本 次 研 究 所 选 取 的 某 露 天 煤 矿 在
2009—2010 年共进行 16 次有效爆破作业成果,通过
对该矿区近年来的现场抛掷爆破实验进行数据搜集,
共获得了 188 组适用数据,选取其中 178 组观测数据
作为 GA-ELM 预测模型的训练数据,其余 10 组观测
数据作为测试样本数据。经整理计算,表 1 显示部分
输入参数。
3. 1 露天煤矿抛掷爆破 GA -ELM 预测模型参数优
化
本文应用遗传算法工具箱 Sheffield
[29],相关参数
设置如下:个体编码方式为二进制编码并设置二进制
位数为 10,适应度函数采用排序的适应度分配函数:
FitnV = ranking( obj) ,其中 obj 为预测样本的预测值
与期望值误差矩阵的范数,选择算子采用随机遍历抽
样,交叉算子采用单点交叉算子。遗传算法其他参数
分别设置为:种群个体数目 NIND 为 20,最大遗传代
数 MAXGEM 为 100,交叉概率 px
为 0. 7,变异概率 pm
为 0. 01,代沟 GGAP 为 0. 95。
ELM 算法,作为单隐层前馈神经网络的一种,如
何确定隐含层的神经元数目十分关键。隐含层神经
元数太少会导致 ELM 网络“欠适配”,而隐含层神经
元数太多会导致 ELM 网络“过适配”。为解决这个
中国煤炭期刊网 www.chinacaj.net
436
剖面
编号
09e1-2
09e1-3
09e1-4
09e1-5
…
煤
炭
学
报
2017 年第 42 卷
表 1 GA-ELM 模型输入数据集( 部分)
Table 1
Input data set of GA-ELM model ( part)
输入层参数
输出层参数
炸药单耗 /
( kg·m-3 )
台阶高
度 / m
煤层
采空区
采空区
厚度 /
上口
下口
m
款 / m
宽 / m
抵抗
坡面
孔距 /
排距 /
剖面
有效抛
松散
线 / m
角 / ( °)
m
m
宽 / m
掷率 / %
系数
抛掷
距离
Weibull 分
布曲线控
/ m
制参数 ɑ
Weibull 分
布曲线控
制参数 β
计算松
散系数
ξ'
0. 72
0. 72
0. 72
0. 72
…
34. 61
28. 50 131. 38 95. 91
31. 40
30. 55 128. 37 92. 96
30. 60
35. 11 132. 99 99. 50
30. 36
30. 88 133. 50 98. 70
…
…
…
…
8. 0
7. 0
7. 0
7. 0
…
56
59
63
66
…
11. 20 7. 11 66. 50
28. 2
1. 02
220. 99 1. 171 4
1. 076 9 1. 296 3
11. 00 8. 18 75. 11
35. 3
1. 07
202. 08 1. 110 2
1. 034 5 1. 511 1
11. 00 8. 25 74. 56
36. 7
1. 18
205. 89 1. 107 4
1. 091 3 1. 561 4
11. 01 8. 01 74. 19
33. 0
1. 19
197. 89 1. 116 9
0. 090 4 1. 511 5
…
…
…
…
…
…
…
…
问题,往往需要设计者反复多次实验和根据自身经验
来确定。本文借鉴“试错法”的思想,为优化测试样
本的仿真误差,通过不断试验和连续地改变系统的隐
含节点参数,试验系统所作出的应答,进而寻求较优
的 GA-ELM 网络拓补结构。
综合以上情况,本文选取 GA -ELM 网络隐含层
单元数区间为[1,100],测试该预测模型隐含层节点
数从 1 到 100,100 种不同情况下 GA-ELM 网络的仿
真误差,比较其拟合情况。输出模型仿真误差参数如
图 3 所示。
由于图 3 中试验结果无法精确选出最优参数,进
一步通过仿真误差值从上述试验结果中选取隐含层
节点参数较优的几种情况,再次进行 3 次重复测试,
见表 2。
图 3 抛掷爆破预测模型———“试错法”过程
Fig. 3 Prediction model of cast blasting-“trial-and-
error”procedure
表 2
3 次重复测试结果比较
Table 2 Triplicate comparing test results
隐含层节点
13
18
23
34
39
47
57
63
83
测试第 1 次仿真误差 0. 334 1
测试第 2 次仿真误差 0. 149 9
测试第 3 次仿真误差 0. 160 5
0. 196 4
0. 334 1
0. 162 4
0. 096 3
0. 470 8
0. 613 4
0. 844 8
0. 364 9
0. 322 8
0. 621 4
0. 292 9
4. 93×10 -4
0. 907 0
0. 862 8
0. 315 7
0. 144 9
0. 189 9
0. 087 5
0. 209 2
0. 137 7
0. 367 2
0. 348 7
0. 668 7
0. 200 4
依 据 表 2 中 相 关 数 据 ,当 参 数 设 置 为 39 时 ,
系 统 仿 真 误 差 在 3 次 重 复 测 试 中 均 明 显 优 于 其
他 情 况 ,最 终 确 定 本 次预 测 模 型 隐 含 节 点 参 数 为
39 。
3. 2 GA-ELM 与 ELM 模型优化预测效果比较
将表 2 所得最优参数,作为 GA -ELM 网络的隐
含层节点数,并选择 Sig 为激活函数,代入训练样本
数据,通过多次试验,GA -ELM 预测模型经过 100 次
进化计算得到一个最佳适应度的稳定迭代值[30]如图
4 所示。
在确定最大遗传变异代数为 100 的基础上,进行
试验,获得最优的 GA -ELM 模型。训练样本的回代
准确率为 100%,决定系数 R2 值为 1。继续使用优化
图 4 进化过程
Fig. 4 Evolutionary process
后的 GA-ELM 模型对测试样本进行仿真预测,并与
未经遗传算法优化的 ELM 模型进行比较,测试样本
预测结果如图 5 所示。
中国煤炭期刊网 www.chinacaj.net
第 3 期
温廷新等:参数优化 GA-ELM 模型在露天煤矿抛掷爆破的预测
536
由图 5 可知,采用 GA -ELM 模型对该露天煤矿
抛 掷 爆 破 进 行 预 测 的 结 果 分 别 为: 有 效 抛 掷 率
的 MSE 值为 0. 258 0,决定系数 R2
值为 0. 986 4;松散
系数的 MSE 值为 1. 748 5 × 10 -4,决定系数 R2 值为
0. 970 6;抛掷距离的 MSE 值为 3. 618 4,决定系数 R2
值为 0. 995 3;而采用 ELM 模型对露天煤矿抛掷爆破
进行预测 的 结 果 分 别 为: 有 效 抛 掷 率 的 MSE 值 为
3. 658 0,决定系数 R2 值为 0. 888 5;松散系数的 MSE
值为 0. 002 7,决定系数 R2 值为 0. 717 9;抛掷距离
的 MSE 值为 99. 199 7,决定系数 R2 值为 0. 830 8;通
过对比可知,采用遗传算法优化后的 ELM 模型相较
于 ELM 模型预测效果具有更高的准确率。
3. 3 不同模型预测效果对比分析
为进一步测试 GA-ELM 预测模型是否较其他智
能算法模型在露天煤矿抛掷爆破预测方面具有优越
性,本 文 选 择 常 用 的 BP 神 经 网 络,RBF 神 经 网
络,SVM 模型以及 GA-BP 神经网络模型与 GA-ELM
模型进行预测精度对比,从均方误差 MSE,决定系数
R2 以及预测精度,3 个角度对模型的精度分别进行描
述对比情况,见表 3。
由表 3 可以看出,经综合比较,采用单一地 ELM
网络,BP 网络,RBF 或者 SVM 进行模型预测,预测结
果的均方误差 MSE 均高于通过遗传算法改进的 2 个
模型,且预测精度和决定系数 R2 较小。由此反映出
图 5 测试集( 有效抛掷率,松散系数,抛掷距
离) -预测结果对比
Fig. 5 Test set( Effective stripping ratio,Loose coefficient,
Thrown distance) comparison of prediction results
表 3 不同模型预测效果测试样本结果比较
Table 3 Comparison of different model effects
均方误差 MSE
有效抛掷率
松散系数
抛掷距离
有效抛掷率
决定系数 R2
松散系数
抛掷距离
仿真误差
0. 258 0
1. 748 50×10 -4
3. 618 4
0. 986 4
0. 970 6
0. 995 3
0. 137 7
3. 658 0
0. 002 70
99. 199 7
0. 888 5
0. 717 9
0. 830 8
8. 283 7
7. 521 0
0. 002 30
266. 488 3
0. 589 9
0. 717 1
0. 014 3
25. 861 1
4. 105 0
0. 001 20
64. 673 7
0. 793 2
0. 869 8
0. 109 5
14. 287 4
4. 859 4
0. 000 93
10. 734 6
0. 725 9
0. 879 7
0. 421 9
15. 043 4
模型
GA-ELM
ELM
BP 神经网络
RBF 神经网络
SVM 模型
GA-BP 神经网络 ( 相同参数)
0. 822 8
0. 000 81
4. 326 2
0. 920 5
0. 890 8
0. 920 4
0. 453 6
单一预测模型针对露天煤矿抛掷爆破自适应较弱;而
经过遗传算法优化后,模型预测精度和拟合效果明显
提高,表明遗传算法能有效降低预测模型的误差,提
高模型的泛化能力和预测效果。并且相较于 GA-BP
模型,本文所选用 GA -ELM 模型的各项 MSE,R2 和
仿真误差,在相同参数设置下均优于该模型,具有更
好的预测能力。
3. 4 基于 WeiBull 模型爆堆形态预测
将抛掷爆破实验试验中所得样本作为检验源,限
于篇幅,本 文 仅 列 出 部 分 样 本 ( 表 4 ) ,用 于 已 训
练 GA-ELM 网络进行仿真测试,进行模拟爆堆形态
对比分析。
通过 GA-ELM 网络预测所得 Weibull 控制参数
ɑ,β 预测值与其真实值的误差对比见表 5。
在上述分析基础上,进一步对预测的 Weibull 概
率密度函数反无量纲化,得到预测模拟爆堆的曲线函
数。将该函数预测爆堆形态与真实爆堆形态以及模
拟爆堆形态进行对比,如图 6 所示。
中国煤炭期刊网 www.chinacaj.net
636
煤
炭
学
报
2017 年第 42 卷
表 4 仿真测试输入输出层样本参数( 部分)
Table 4 Parameter of input and output samples in simulation testing ( Part)
样本
编号
炸药单
耗 /
(kg·m-3 )
台阶高
度 / m
…
e23-11
e23-12
…
0. 67
0. 67
e23-13
0. 67
e23-14
0. 67
…
…
…
36. 3
37. 9
37. 5
37. 9
…
煤层
采空区
采空区
厚度 /
上口
下口
m
…
款 / m
宽 / m
…
…
28. 5
110. 3
89. 32
28. 6
107. 8
88. 54
25. 7
114. 5
91. 36
25. 9
113. 9
89. 81
…
…
…
抛掷
距离 /
Weibull 曲
线控制
参数 ɑ
Weibull 曲
线控制
参数 β
计算松
散系数
ξ'
m
…
抵抗
坡面
孔距 /
排距 /
剖面
有效抛
松散
线 / m
角 / (°)
m
m
宽 / m
掷率 / %
系数
…
7. 0
7. 0
7. 0
7. 0
…
…
63
54
66
71
…
…
…
…
…
…
…
…
11. 00
7. 0
72. 25 35. 46
1. 22
228. 00
1. 068 0
1. 097 9
1. 311 9
11. 00
7. 0
75. 92 36. 48
1. 16
221. 13
1. 088 0
1. 050 0
1. 530 0
11. 01
7. 0
70. 33 36. 07
1. 24
217. 12
1. 138 0
1. 100 0
1. 600 0
11. 01
7. 0
70. 29 32. 06
1. 14
217. 00
1. 170 0
1. 526 0
1. 548 1
…
…
…
…
…
…
…
…
…
表 5 Weibull 模型控制参数 ɑ,
β
误差比较
Table 5 Error comparison of Weibull model’s control parameters ɑ,
β
控制参数 ɑ
控制参数 β
计算松散系数 ξ'
预测值
真实值
1. 113 7
1. 114 8
1. 120 0
1. 123 2
1. 068 0
1. 088 0
1. 138 0
1. 170 0
误差 / %
4. 275 9
2. 647 5
1. 581 3
3. 998 7
预测值
真实值
1. 106 7
1. 081 2
1. 105 3
1. 496 5
1. 097 9
1. 050 0
1. 100 0
1. 526 0
误差 / %
0. 798 5
2. 967 2
0. 481 2
1. 935 2
预测值
真实值
1. 356 6
1. 546 7
1. 524 2
1. 506 8
1. 311 9
1. 530 0
1. 600 0
1. 548 1
误差 / %
3. 409 3
1. 091 9
4. 739 9
2. 665 2
样本编号
e23-11
e23-12
e23-13
e23-14
4 结
论
(1) 在采用遗传算法对 ELM 输入权值矩阵和隐
含层偏差进行优化,建立 GA -ELM 预测模型的基础
上,根据“试错法”思想,以预测模型仿真误差为评判
依据,通过两次筛选确定最优隐含层节点数。预测模
型在保证模型自身泛化能力的同时,避免输入权值矩
阵、隐含层偏差和隐含层节点数随机性对 ELM 预测
精度的影响,一定程度上减少网络过拟合所造成的隐
患,提高模型预测的准确率。
(2) 通过对露天煤矿抛掷爆破有效抛掷率,松散
系数,抛掷距离的预测可以保证拉斗铲倒堆工艺系统
的可靠性,为优化抛掷爆破设计、改善抛掷爆破效果
提供有力的依据;进一步预测模拟爆堆形态,为实现
计算机技术模拟爆破过程和预测爆破效果提供理论
基础;对实际生产具有一定的指导意义。
(3) 在关于露天煤矿 GA-ELM 模型预测中,实际
观测样本不可避免受噪声影响,进而对预测效果产生
一定的偏差;此外遗传算法对 ELM 的改进,尽管一定
程度上提升模型预测精度,但预测效率相对较低,耗
费大量的时间成本,占用运载空间较大;为降低噪声
影响,笔者应当进一步还原检测数据自身准确性,另
外如何寻找到预测精度和预测效率的平衡点也是模
型下一步需要改进的地方。
图 6 预测爆堆形态与真实爆堆形态以及模拟爆堆形态对比
Fig. 6 Comparison among prediction blasting form,reality of
blasting form and blasting simulation of blasting form
中国煤炭期刊网 www.chinacaj.net
第 3 期
温廷新等:参数优化 GA-ELM 模型在露天煤矿抛掷爆破的预测
736
参考文献( References) :
[1] 李克民,张幼蒂,傅洪贤.露天煤矿抛掷爆破参数分析[J]. 采矿
与安全工程学报,2006,23(4) :423-426.
Li Kemin,Zhang Youdi,Fu Hongxian. Analysis of casting blast pa-
rameters in surface coal mines[J].Journal of Mining & Safety Engi-
neering,2006,23(4) :423-426.
[2] 周伟,才庆祥,李克民. 露天煤矿抛掷爆破有效抛掷率预测模型
[J].采矿安全与工程学报,2011,28(4) :614-617.
Zhou Wei,Cai Qingxiang,Li Kemin. Prediction model of effective
stripping ratio of casting blast in open pit[J]. Journal of Mining
& Safety Engineering,2011,28(4) :614-617.
[3] 李祥龙,刘殿书,何丽华,等.露天煤矿的台阶高度对抛掷率的影
响[J].爆炸与冲击,2012,32(2) :211-215.
Li Xianglong,Liu Dianshu,He Lihua,et
of
bench height of an open-pit coal mine on cast percentage[J].Explo-
sion and Shock Waves,2013,32(2) :211-215.
al. Influences
[4] 马力,李克明,丁小华,等.抛掷爆破岩体炮制距离影响因素研究
[J].工程爆破,2013,19(1-2) :50-53.
Ma Li,Li Keming,Ding Xiaohua,et al.Study on the influencing fac-
tors on casting distance of rock by casting blast[J]. Engineering
Blasting,2013,19(1-2) :50-53.
[5] 孙文希,刘希亮,王洪斌,等. 基于 MIV 的抛掷爆破影响因子权
重分析[J].中国矿业大学学报,2012,41(6) :993-998.
Sun Wenxi,Liu Xiliang,Wang Hongbin,et al.Weight analysis of cast
blasting effective factors based on MIV method[J].Journal of China
University of Mining & Technology,2012,41(6) :993-998.
[6] 刘希亮,赵学胜,陆锋,等. 基于 GA - SVM 的露天煤矿抛掷爆破
抛掷率预测[J].煤炭学报,2012,37(12) :1999-2005.
Liu Xiliang,Zhao Xuesheng,Lu Feng,et al.A GA-SVM based model
for throwing rate prediction in the open-pit cast blasting[J]. Journal
of China Coal Society,2012,37(12) :1999-2005.
对民房破坏的预测分析[J].中国安全生产技术,2015,11( 8) :
119-124.
Wen Tingxin,Zhu Chengwei,Kong Xiangbo. Predicting analysis
on damage to residential house by blasting vibration in open
pit mining based on PCA-ELM model[J].Journal of Safety Science
and Technology,2015,11(8) :119-124.
[12] 韩亮,刘殿书,李洪江,等.基于 Weibull 模型的高台阶抛掷爆破
爆堆形态 BP 神 经 网 络 预 测[J]. 煤 炭 学 报,2013,38 ( 11 ) :
1947-1952.
Han Liang,Liu Dianshu,Li Hongjiang,et al. BP neural network
forecast of blasting muck pile form of high bench cast blasting
based on Weibull model[J]. Journal of China Coal Society,2013,
38(11) :1947-1952.
[13] Huang Guangbin,Zhu Qinyu,Siew C K.Extreme learning machine:
Theory and application[J].Neurocomputing,2006,70( 1 /3) :489-
501.
[14] Huang Guangbin,Zhu Qinyu,Siew C K.Real-time learning capabil-
ity of neural networks[J].Neurocomputing,2006,70:863-878.
[15] Lan Yuan,Soh Y C,Huang Guangbin. Ensemble of online sequen-
tial extreme learning machine[J]. Neurocomputing,2009,72 ( 13 /
15) :3391-3395.
[16] Huang Guangbin,Ding Xiaojian,Zhou Hongming. Optimiza-
tion method based extreme learning machine for classification
[J].Neurocomputing,2010,74(1 /3) :155-163.
[17] 张玲,张钹.遗传算法机理的研究[J]. 软件学报,2000,11( 7) :
945-952.
Zhang Ling,Zhang Bo.Research on the mechanism of genetic algo-
rithms[J].Journal of Software,2000,11(7) :945-952.
[18] 戴晓辉,李敏强,寇纪淞.遗传算法理论研究综述[J].控制与决
策,2000,15(3) :263-268.
Dai Xiaohui,Li Minqiang,Kou Jisong.Survey on the theory of ge-
netic algorithms[J]. Control and Decision,2001,15 ( 3 ) : 263 -
[7] 孙文希,刘希亮,谭正龙,等.基于抛掷爆破预测的 BP 神经网络
268.
参数优化方法[J].煤炭学报,2012,37(1) :59-64.
Sun Wenxi,Liu Xiliang,Tan Zhenglong,et al.Parameter optimization
of BP-neural network based on the forecast of cast blasting[J].Jour-
nal of China Coal Society,2012,37(1) :59-64.
[8] 刘干,李克民,肖双双,等. 露天煤矿抛掷爆破有效抛掷率预测
[J].金属矿山,2014(4) :65-69.
Liu Gan,Li Kemin,Xiao Shuangshuang,et al. Prediction of effective
stripping ratio of casting blast in surface coal mines[J].Metal Mine,
2014(4) :65-69.
[9] 周志坚,毛宗源. 一种最优模糊神经网络控制器[J]. 控制与决
策,2000,15(3) :358-360.
Zhou Zhijian,Mao Zongyuan. An optimal fuzzy neural network con-
troller[J].Control and Decision,2000,15(3) :358-360.
[10] 温廷新,戚磊,邵良杉. 露天采矿爆破振动特征参量的 Logistic-
ELM 预测[J].计算机工程与设计,2015,36(10) :2791-2795.
Wen Tingxin,Qi Lei,Shao Liangshan. Prediction blasting vibra-
tion characteristic parameters in open-pitmining based on Logistic-
ELM[J].Computer Engineering and Design,2015,36(10) :2791-
2795.
[11] 温廷新,朱成伟,孔祥博.基于 PCA-ELM 模型的露采爆破振动
[19] 玄光男,程润伟,于歆杰,等. 遗传算法与工程优化[M]. 北京:
清华大学出版社,2004:24-30.
Xuan Guangnan,Cheng Runwei,Yu Xinjie,et al.Genetic algorithms
and engineering optimization [M]. Beijing: Tsinghua University
Press,2004:24-30.
[20] 丁小华,李克民,狐为民,等.基于非线性理论的抛掷爆破爆堆
形态预测[J].中国矿业大学学报,2012,41(5) :764-769.
Ding Xiaohua,Li Kemin,Hu Weimin,et al. Prediction of opti-
mum muck pile casting shape during blasting: A nonlinear theory
[J]. Journal of China University of Mining & Technology,2012,
41(5) :764-769.
[21] 傅洪贤,李克民.露天煤矿高台阶抛掷爆破参数分析[J]. 煤炭
学报,2006,31(4) :442-445.
Fu Hongxian,Li Kemin.Analysis of high bench cast blasting param-
eters in surface coal mines[J]. Journal of China Coal Society,
2006,31(4) :442-445.
[22] 缪伟,赵庆,王斯日古楞.哈尔乌素露天煤矿抛掷爆破参数的确
定[J].煤矿安全,2013(11) :211-213.
Miao Wei,Zhao Qing,Wang Siriguleng. Determination of casting
blast parameters for Haerwusu open-pit coal mine[J]. Safety in
中国煤炭期刊网 www.chinacaj.net