基于中心复合响应面的覆盖层岩块粒径与密度实验分析_杨贺.pdf

收稿日期2019-12-20 基金项目国家 “十二五” 科技支撑计划项目 （编号 2013BAB02B08） 。 作者简介杨贺 （1994） ， 男， 硕士研究生。通讯作者路增祥 （1965） ， 男， 教授， 博士。 总第 528 期 2020 年第 6 期 金属矿山 METAL MINE 基于中心复合响应面的覆盖层岩块粒径与 密度实验分析 杨贺路增祥张国建 1 （辽宁科技大学矿业工程学院， 辽宁 鞍山 114051） 摘要无底柱分段崩落法崩落的矿石与覆盖层岩石直接接触， 使矿石损失贫化加大。覆盖层是保证井下正 常生产的关键， 覆盖层形成不合理会加大矿石的损失贫化。针对覆盖层岩石颗粒的粒度与密度， 基于中心复合设 计进行了实验， 建立了粒度与密度对岩石混入率影响的二阶响应预测模型， 对实验结果进行方差分析与显著性检 验， 检验结果表明预测模型的显著性较高， 实验值与预测值的拟合程度良好。通过绘制三维响应面图， 分析了粒度 与密度的交互作用对混入率的影响， 结果表明覆盖层岩石颗粒的粒度与密度对岩石混入率的影响情况呈现出非线 性， 当岩石颗粒的密度值在3.05～3.46103kg/m3时， 增加岩石颗粒粒度有助于减少岩石混入。 关键词无底柱分段崩落法覆盖层粒度密度中心复合设计响应面岩石混入率 中图分类号TD853文献标志码A文章编号1001-1250 （2020） -06-150-06 DOI10.19614/ki.jsks.202006023 Experimental Analysis of Particle Size and Density of Cover Layer Based on Central Composite Response Surface Yang HeLu ZengxiangZhang Guojian2 （School of Mining Engineering， University of Science and Technology Liaoning， Anshan 114051， China） AbstractThe caved ores by non-pillar sublevel caving have directly contact with the overburden rock，which makes the ore loss and dilution． The cover layer is extremely important for ensuring the normal operation of the underground． The unreasonable ation of cover layer will increase the loss of ore． According to the grain size and density of the cover rocks，the experiment is carried out based on the central composite design，and the second-order response prediction model for the influence of grain size and density on the rock mixing rate is established． The experimental results are analyzed by variance analysis and significance test． The test results show that the significance of the prediction model is high，and the ex⁃ perimental value fit well with the prediction value． The influence of the interaction both particle size and density on the mix⁃ ing rate was analyzed by drawing the three-dimensional response surface． The experimental results showed that the influence of size and density of the cover layer rocks on the rock mixing rate was nonlinear． When the density of rock particles ranged from（3．05～3．46） 103kg/m3， increasing the size of rock particles would help to reduce rock mixing． KeywordsNon-pillar sublevel caving ， Cover layer， Grain size， Density， Central composite design， Response surface， Rock mixing rate Series No． 528 June 2020 无底柱分段崩落法因其开采方式简单、 安全、 机 械化程度高， 在我国得到了较快的发展， 我国的和睦 山铁矿 ［1］、 小官庄铁矿［2］、 眼前山铁矿［3］、 北洺河铁 矿 ［4］、 大红山铁矿［5］、 梅山铁矿［6］等矿山采用无底柱分 段崩落法开采， 取得了很好的应用效果。但是， 该采 矿方法的的特点是覆岩下放矿， 崩落的矿石与废石 直接接触， 导致了矿石损失率增加， 贫化率加大， 损 失贫化大成为影响该采矿方法使用的一大技术难 题。 为尽可能地减少矿石的损失和贫化， 专家和学 者对覆盖层进行了大量的研究工作。如朱志根等 ［7］ 通过物理模型实验发现， 矿岩块度不均匀性是引起 150 ChaoXing 矿石贫化的原因之一， 并且废石层细小颗粒含量增 加时矿石贫化率明显增加。张国建等 ［8-9］发现了覆盖 岩层的自然分级现象， 分级的程度与大小颗粒的分 布有关。自然分级现象的发现， 为降低矿石损失贫 化提供了新的方法。张志军等 ［10］通过物理模拟实验 证明了自然分级现象的存在， 并且， 覆盖岩层下移过 程中会不断出现平衡拱的形成与破坏的现象， 其中 废石颗粒越小移动速度越大， 在截止品位放矿条件 下， 覆盖层岩石块度大有利于放矿。徐国元等 ［11］运 用数值模拟的方法研究覆盖岩层下矿岩流动特性， 发现离散元数值模拟能够较准确、 直观地反映矿石 与废石的运动规律以及矿岩接触面变化情况， 增大 覆盖层岩石块度有利于减少矿石贫化。张杨等 ［12］从 矿山选取不同容重的岩石作为覆盖层进行物理模型 实验， 发现岩块的容重对岩石混入率有一定影响， 在 放矿的其他条件不变情况下， 岩石混入率随着覆盖 层岩块容重的增大而增大。 无底柱分段崩落法应用过程中， 组成覆盖层岩 石的粒径与密度不同， 对矿石损失贫化的影响也不 同。本项目通过中心复合设计响应面进行物理模拟 实验， 分析不同的岩石颗粒粒度和密度对岩石混入 率的影响规律， 建立岩石混入率的回归预测模型， 最 后根据所得到的回归预测模型进行优化分析， 以期 为矿山生产提供理论指导。 1实验设计 1. 1中心复合响应面实验设计 响应面法是一种基于数学和统计学的方法， 可 以对多因素、 多变量进行目标优化， 除了分析自变量 的影响外， 该法还可生成数学模型， 分析因素间的相 互作用， 根据分析结果可进行新的产品设计 ［13］。中 心复合设计 （CCD） 是比较常用的响应面设计， 可以量 化一个或多个测量的输出响应与输入响应之间的关 系 ［14］， CCD在二水平因子设计点的基础上， 增加了轴 向点和中心点来完善响应面的模拟 ［15］。中心点是用 来进行非线性测试与误差估计； 轴向点是用来预测 回归方程的。中心点和轴向点设计示意图如图1所 示。 在进行实验前要设计CCD的编码点， CCD编码 点是按因素数与水平数进行选取的。本次实验考 虑的因素为覆盖层岩石颗粒的密度与粒度， 选取的 覆盖层岩石材料为石英、 角闪石、 石榴子石， 密度值 分别为 2.64103kg/m3、 3.06103kg/m3和 4.07103 kg/m3， 并且每种材料选取 3个粒度 5 mm、 10 mm、 15 mm， 综合考虑为二因素三水平， 依据标准的二因素三 水平CCD编码点选取原则， 选取编码值如表1所示。 1. 2实验方案 根据选取的因素数与水平数， 进行了中心复合 实验设计， 共进行13次实验， 且每次实验的密度与粒 度编码值点如表2所示。 在进行实验前要设计放矿箱尺寸和放矿参数， 依据矿山实际的设计参数， 按1 ∶50比例设计了实验 室物理实验模型， 模型的长、 宽、 高分别为 400 mm、 300 mm、 800 mm， 进路口尺寸为80 mm70 mm， 崩矿 步距大小为60 mm， 铺设的矿石层高度为240 mm， 覆 盖岩层厚度为500 mm。实验过程如图2所示。 首先进行放矿箱的装填， 如图 2 中 （a） 所示， 先 在放矿箱中装入粒度为 10 mm的磁铁矿， 当矿石层 2020年第6期杨贺等 基于中心复合响应面的覆盖层岩块粒径与密度实验分析 151 ChaoXing 高度达到240 mm刻度线时停止装矿， 进行覆盖层的 装填， 当覆盖层高度达到 740 mm 刻度线时停止装 填， 此时放矿箱装填完毕。然后开始放矿， 如图2中 （b） 所示。当放出的矿石品位达到截止品位时， 停止 放矿， 如图 2中 （c） 所示， 最后进行实验数据的统计 与处理。 为确定各因素之间的关系， 将密度、 粒度2个独 立变量转化为无量纲变量x1和x2， 变量安排及实验次 数如表 2， 因素变量与输出结果间的函数关系式如 下 y β0∑ j 1 p β0 xj∑ i ＜ j βijxixj∑ j 1 p βijx2 j ε．（1） 式 （1） 表示了预测的反映值 （y） 与回归系数 （β0） 、 自变 量 （x） 、 因素个数 （p） 、 误差项 （ε） 间的关系。 2实验结果及分析 实验主要考虑覆盖层岩石颗粒的密度与粒度对 岩石混入率的影响， 混入率依据当次混入率的变化 进行选取。混入率指的是总的放出量中岩石所占的 比例， 衡量的是总体的岩石混入情况， 而当次混入率 指的是当次放出量中废石所占的比例， 衡量的是每 一次出矿岩石的混入情况。图3为3种覆盖层岩石 材料的混入率与当次混入率变化曲线。 从图3可以看出， 随着出矿次数的增加， 混入率 与当次混入率都在增加， 但是在当次混入率达到图 中虚线时， 当次混入率变化明显趋于平稳。若假定 每一铲子的铲取量相同， 随着出矿次数的增加， 每一 铲子中矿石所占比例明显减小， 在当次混入率达到 图中虚线时， 矿石所占比例小于10， 不利于回收， 而 此时混入率依然在升高。根据矿石资源最大程度地 回收以及减少采矿成本的原则， 选取在当次混入率 达到90时的混入率值作为实验结果， 研究密度与 粒度及它们的相互作用对混入率的影响， 实验结果 如表3。 3中心复合设计和拟合回归模型 采用CCD研究了放矿效果评价指标 （混入率） 和 2个可以控制的因素 （覆盖层岩石粒度和密度） 之间 的关系， 这些关系可以利用式 （1） 进行表述， 该关系 式不仅包含2个可控因素粒度与密度对放矿评价指 标的一阶影响， 还有2个因素的相互作用的影响。将 所得实验数据结果进行方差分析， 得到拟合回归方 程如下 金属矿山2020年第6期总第528期 152 ChaoXing y1 224．717 - 118．085x1 1．843x2- 0．578x1x2 18．681x2 1 0．006x 2 2． 2 从方差分析中获得的统计参数如表4所示。其 中， R2代表拟合回归方程的拟合度， 拟合度越高， 表示 方程可靠性越大； Prob ＞ F 代表显著性水平， 小于 0.05表示显著； Std. Dev.表示标准差； C.V. 表示变异 系数， 用来度量数据的分散程度； PRESS表示预测残 差平方和， 提供模型与观察样本拟合的总度量； Mean 代表样本均值。 从表4可以看出， 方程拟合效果较好， 各因素对 实验结果影响效果显著， 可以用于分析密度与粒度 对混入率的影响。 为提高模型预测精度， 需要对各次实验结果的 异常值进行分析。残差可以用来判断线性回归中的 异常点， 学生化残差是残差与它的标准差的比值， 即 使误差的标准差相同， 但是样本中的标准差也会因 为数据点的不同而不同。学生化外残差是学生化残 差的一种， 是异常点检验的t统计量， 可以有效地预 测样本差异点。 图4中在刻度4、 -4之间的横线表示置信区间， 实验运行点超出置信区间表明该点为差异点。图中 各实验点均在置信区间内， 差异效果良好， 因此进行 了各实验结果点的预测， 表5列出了实验值与预测值 的结果比较。 为了清晰直观地对比实验值与预测值之间的差 异， 按照表5绘制了预测值与实验值比较图， 见图5。 预测值与实验值曲线相互重叠， 预测曲线变化 趋势与实验值曲线相一致。第4次、 第5次、 第6次、 第8次、 第11次、 第12次实验预测值与实验值曲线变 化值相差不多， 第7次预测值与实验值相差较大， 但 仅为3， 预测值误差在理想范围内， 预测准确。 由式 （2） 可知， 覆盖层岩石的颗粒密度与粒度的 相互性对岩石混入率也存在影响， 图6反映了2因素 的相关性分析曲线。2条曲线中A表示密度， B表示 粒度。岩石密度在3.44103kg/m3时， A、 B2条曲线相 交， 在交点外A曲线或高于或低于B曲线， 表明2因素 的相互性对岩石混入率有影响。 针对2因素的相互性对岩石混入率有影响的情 2020年第6期杨贺等 基于中心复合响应面的覆盖层岩块粒径与密度实验分析 153 ChaoXing ［1］ ［2］ ［3］ ［4］ ［5］ ［6］ ［7］ 况， 绘制了2因素的扰动图。扰动图有助于比较设计 中各因素对实验结果的敏感程度， 曲线陡峭， 表明该 因素对实验结果变化敏感； 曲线相对平坦， 表明该因 素对实验结果变化不敏感。如图7所示， A曲线陡峭， B曲线平坦， 由此可知， 密度对混入率的影响更为敏 感。 通过方差分析与回归模型的建立， 得知， 粒度与 密度的交互作用对混入率有明显的影响， 对于2因素 的交互作用， 使用三维响应面图进行分析是一个很 好的解释手段 ［16-18］。由于回归模型的预测值整体变 化趋势与实验值相符， 且误差很小， 根据此模型， 建 立覆盖层岩块粒径与密度对混入率影响的响应面， 如图8、 图9所示。 图9中响应面表现为中间低、 两端高且前后倾斜 的U字形状。当密度值在 （3.05～3.46） 103kg/m3区 间内， 混入率随着粒度的增大逐渐减小， 在粒度值达 到15 mm时， 混入率达到最低值35.78； 当密度值小 于3.05103kg/m3， 混入率随着粒度的增加而增加； 当密度值大于3.46103kg/m3， 混入率随着粒度的增 加逐渐减小。混入率或低或高的变化表明粒度与密 度对混入率的影响是非线性的。 4结论 （1） 方差分析和拟合度结果表明建立的二阶响 应面预测模型可以用于覆盖岩层岩石颗粒粒度与密 度对岩石混入率的预测， 预测曲线变化趋势与实验 值曲线基本一致。 （2） 覆盖岩层岩石颗粒粒度与密度对岩石混入 率影响效果显著， 显著性水平为0.006， 远小于0.05。 （3） 覆盖层岩石颗粒粒度与密度对岩石混入率 的影响呈现出非线性， 当密度值在 （3.05～3.46） 103 kg/m3区间内时， 增加岩石颗粒粒度有助于减少岩石 混入率。 （4） 调整覆盖层岩石颗粒粒度与密度大小， 有助 于减小岩石混入率和降低矿山开采成本。 参 考 文 献 黄小忠， 刘发平． 和睦山铁矿的新型高效采矿方法 ［J］ ． 采矿技 术， 2008， 8 （6） 1-2． Hang Xiaozhong， Liu Faping． The new efficient mining in Hemushan iron mine ［J］ ． Mining Technology， 2008， 8 （6） 1-2． 辛国帅． 小官庄铁矿北区地下开采对地表沉陷及井筒变形影响 数值模拟研究 ［D］ ． 沈阳 东北大学， 2014． Xin Guoshuai． Numerical Simulation Study on the Influence of Un⁃ derground Mining on Surface Subsidence and Shaft Deation in North Region of Xiaoguanzhuang Iron Mine［D］ ． Shenyang Norteasten University， 2014． 李鹏． 无底柱分段崩落法铲装效率实时监测研究 ［D］ ． 鞍山 辽宁科技大学， 2018． Li Peng． Study on Remote Real-Time Lifting Effect Monitoring of Underground Mine ［D］ ． AnshanUniversity of Science and Technol⁃ ogy Liaoning， 2018． 何荣兴， 任凤玉， 宋德林， 等． 无底柱分段崩落法崩矿步距优化 方法及应用 ［J］ ．矿业研究与开发， 2015， 35 （8） 1-4． He Rongxing， Ren Fengyu， Song Delin， et al． Optimization and ap⁃ plication on caving space of non-pillar sublevel caving ［J］ ． Mining Research and Development， 2015， 35 （8） 1-4． 陈发兴， 张志雄． 大参数无底柱分段崩落法在大红山铁矿的运 用 ［J］ ． 有色金属设计， 2009， 36 （3） 18-21． Chen Faxing， Zhang Zhixiong． Application of big parameter of no- pillar sublevel caving in Dahongshan iron mine ［J］ ． Nonferrou Met⁃ als Design， 2009， 36 （3） 18-21． 金闯， 董振民， 贡锁国， 等． 梅山铁矿无底柱分段崩落法加大 结构参数研究 ［J］ ． 金属矿山， 2000 （4） 16-19． Jin Chuang， Dong Zhenmin， Gong Suoguo， et al． Study on the en⁃ larged-structure parameters of sublevel caving in Meishan iron mine ［J］ ． Metal Mine， 2000 （4） 16-19． 朱志根， 吴爱祥， 习泳． 自然崩落法覆盖岩层下放矿贫化实验 金属矿山2020年第6期总第528期 154 ChaoXing ［8］ ［9］ ［10］ ［11］ ［12］ ［13］ ［14］ ［15］ ［16］ ［17］ ［18］ 研究 ［J］ ． 江西有色金属， 2005， 19 （4） 7-9． Zhu Zhigen， Wu Aixiang， Xi Yong． Study of drawing dilution under overburden in block caving based on experiments ［J］ ． Jiangxi Non⁃ ferrous Metals， 2005， 19 （4） 7-9． 张国建， 徐国权． 崩落法松散覆盖岩层自然分级对矿石损失贫 化的影响 ［J］ ． 金属矿山， 2010 （5） 21-24． Zhang Guojian， Xu Guoquan． Effect of the natural grade classifica⁃ tion of loose-covered rock with caving on ore loss and dilution ［J］ ． Metal Mine， 2010 （5） 21-24． 张国建， 蔡美峰．无底柱分段崩落法应进一步研究的几个问题 ［J］ ． 中国矿业， 2003， 12 （11） 41-43． Zhang Guojian， Cai Meifeng． Some critical issues in application of pillarless sublevel caving ［J］ ． China Mining Magazine， 2003， 12 （11） 41-43． 张志军， 明世祥， 宋洪勇． 崩落法采场松散覆盖岩层移动规律研 究 ［J］ ． 采矿技术， 2008， 8 （2） 11-13． Zhang Zhijun， Ming Shixiang， Song Hongyong． Study on the move⁃ ment law of loose cover laywer in caving stope ［J］ ． Mining Technol⁃ ogy， 2008， 8 （2） 11-13． 徐国元， 赵建平， 黄仁东， 等．覆盖岩层下矿岩流动过程数值模 拟 ［J］ ．金属矿山， 2004 （12） 1-4． Xu Guoyuan， Zhao Jianping， Huang Rendong， et al． Numerical sim⁃ ulation of ore and rock flow under overburden ［J］ ． Metal Mine， 2004 （12） 1-4． 张杨， 张国建．覆盖层岩块容重对岩石混入的影响 ［J］ ．现代矿 业， 2019（3） 50-52． Zhang Yang， Zhang Guojian． Influence of rock weight density of overburden stratum on waste in-ore ［J］ ． Modern Mining， 2019 （3） 50-52． Deniz Baş， İsmail H Boyacı． Modeling and optimization I Usability of response surface ology ［J］ ． Journal of Food Engineering， 2007， 78 （3） 836-845． 张家雨， 唐德文， 邹树梁， 等 ． 基于中心复合响应曲面法的 42CrMo高速铣削实验分析 ［J］ ．制造技术与机床， 2019 （2） 94- 100． Zhang Jiayu， Tang Dewen， Zou Shuliang， et al． Experimental analy⁃ sis of 42CrMo high speed milling based on central composite re⁃ sponse surface ology ［J］ ． Manufacturing Technology Ma⁃ chine Tool， 2019 （2） 94-100． 孙开畅， 尹志伟， 李权． 基于中心复合设计的RSM边坡可靠度 研究 ［J］ ． 水力发电， 2018， 44 （4） 19-23． Sun Kaichang， Yin Zhiwei， Li Quan． Slope reliability study of re⁃ sponse surface based on central composite design ［J］ ． Wa⁃ ter Power， 2018， 44 （4） 19-23． Mason R L， Gunst R F， Hess J J． Statistical Design and Analysis of Experiments with Applications to Engineering and Science［M］ ． Hoboken， NJ John Wiley and Sons Inc， 2003． Montgomery D C， Runger G C， Hubele N F． Engineering Statistics ［M］ ． Hoboken， NJJohn Wiley and Sons Inc， 2001． Vining G G．Statistical s for Engineers ［M］ ． LondonDuxburg Press， 2003． 2020年第6期杨贺等 基于中心复合响应面的覆盖层岩块粒径与密度实验分析 155 ChaoXing