矿山空区诱发的岩移特征及覆盖层冒落效应研究

摘要

采空区是采矿工程产生的最主要的灾害形式之一，当暴露面积增大到一定程度时，将可能会诱发灾难性的事故。研究采空区的稳定性及其岩层移动冒落特征和规律意义重大。本文以作者负责的项目“龙桥铁矿采空区监测与顶板岩石冒落规律研究”为研究背景，采用室内试验、数值模拟、理论分析、现场监测与调查相结合的方法，借鉴统计学、多源信息融合技术的思想，研究揭示龙桥铁矿矿体回采产生的开采效应及关键块体特征参数。
　　 系统地对现有的岩石力学参数取值方法存在的缺陷进行分析，给出了用稳健回归方法确定岩石抗剪强度的模型，并对各种岩体参数取值方法进行评述。针对岩石抗剪强度确定中图解法和最小二乘法(LS)存在的问题：(1)只适用于实验数据相关性较高的情况，不适用于数据离散并有异常值存在的情况；(2)在对试验法获得的试验数据进行处理时，传统的最小二乘法由于采用残差平方和，容易夸大试验数据中异常值的影响。鉴于此，提出了岩石抗剪强度参数的稳健(Robust)回归分析方法。该方法在实验数据相关性差、数据离散并有异常值存在的情况下，具备削弱数据离散和对异常值定位的能力，提高了估计参数的稳健性和可靠性。该方法以残差的绝对值之和代替残差平方和，并通过复形法求得力学参数，避免了异常值的二次项，可以有效地减少异常值的影响。研究表明，在试验数据的相关性较好时，两种方法的计算结果相差不大，但当试验数据的相关性较差，存在异常值时，稳健回归方法的计算结果要优于最小二乘法。另外，为方便及提高数值计算精度，测量了龙桥铁矿的地应力。
　　 在前期现场调查、现场原岩应力测定、室内力学试验及岩体力学参数工程处理的研究基础上，采用FLAC3D方法对龙桥铁矿空区形成过程及采空区稳定性进行模拟计算和预测分析，分析了空区围岩的主应力、正应力、剪应力及变形情况，初步得到了围岩体变形和失稳的规律。采空区稳定性数值模拟结果表明，采空区形成后，空区四周各角隅处将首先达到极限剪切破坏状态。随着采空区的增大，角隅处破坏区域逐步延深扩大，顶板中央拉应力分布逐渐明显，最终可能变为拉应力破坏。根据研究结论，指出该矿山后续开采设计和管理中应注意的问题，并提出了合理的开采设计建议。
　　 传统的有限元分析方法很难分析岩体大变形及块体的跨落情况，系统地将关键块体理论、极限平衡分析与可靠度理论相结合，研究关键块体跨落的可能性。在常规可靠度分析理论一次二阶矩方法和二次二阶矩方法中，需要对基本随机变量进行正态变换或当量正态化，必须已知随机变量的概率分布，有时比较麻烦，而且这种变换也是导致功能函数非线性从而成为误差的一种来源。为消除这些误差，本文引入一次渐近积分方法与二次次渐近积分方法，很好的解决了该问题，不再需将变量正态化，提高了计算精度。
　　 科学准确的理解关键块体的特征参数，如关键块体的位置、体积、滑动方向、下滑力及抗滑力，可以有效对采矿方案进行优化，最大限度的提高开采效益，有效防止灾害事故的发生。建立空区上块状岩体的特征分析方法，该方法结合块体理论、可靠度理论及矢量分析可以确定关键块体的大小、移动方向、滑动模式、块体的可靠性及其失效概率等。通过本文建立的模型，根据模型的思路和数学公式编制了Matlab程序进行计算，计算结果显示，在15个开采区域，可动块体有270多个，可动块体中安全系数小于2的关键块体有48个，其中安全系数小于1.5的关键块体有33个，占68.73％。并给出了每个块体特征参数、可靠度及可靠概率。
　　 将安全系数与可靠度比较指出了安全系数在分析中存在的缺陷，由于岩体的复杂性和不确定性，很难准确确定岩石力学计算参数，然而，岩体力学参数的确定直接影响着计算结果，考虑参数的随机性，引入随机可靠度理论对其进行研究。研究发现，用安全系数分析关键块体的稳定性存在很大的误差，甚至有时候会得出完全错误和相反的结论。
　　 根据2009年龙桥铁矿回采后采空区的实际空间分布情况，结合整体稳定性分析对几何建模工作提出的要求，与龙桥铁矿技术人员多次讨论、分析，在确定本研究中矿体开采状况的时空对应关系后，对矿整个矿山开采显现的规律进行研究，主要从最大主应力、最小主应力、产生的垂直位移及地表变形等几个方面展开。
　　 将意大利产的Detector新型探地雷达设备应用到龙桥铁矿采空区覆盖层厚度探测中，并论述了探测原理及具体方法。该设备天线为IDSTR40MHz屏蔽天线，最大探测深度可以达到50米。探测结果表明，该矿山空区垫层厚度大于25m，满足国家安全生产规程规定的无底柱分段崩落法空区垫层厚度必须超过2个分段高度。
　　 讨论了三种GM(1，1)模型与最小二乘支持向量机结合的方法：并联型、串联型和残差型灰色最小二乘支持向量机预测模型。给出了岩石覆盖层冒落及碴石厚度预测的灰色最小二乘支持向量机模型的建模思路及关键步骤，并应用该模型对覆盖层冒落高度及碴石厚度进行预测，研究表明，3种灰色最小二乘支持向量机的预测结果优于单一的GM(1，1)和单纯的LSSVM，各项指标均有所降低，体现了多信息融合的优势。残差型模型要优于并联型和串联型模型，表明用LSSVM来描述覆盖岩层及碴石厚度序列的随机项是正确可靠的。在串联型模型中，该组合预测结果要优于单个GM(1，1)模型预测，在实际应用GM(1，1)模型进行预测时，不能够事先知道哪个模型最优，因此使用该模型可以减小预测的盲目性，有一定的实用价值。充分显示了GLSSVM这种新的信息处理和预测方法的优越性，可以在工程实践中加以推广。