矿井断层构造预测方法
【专利摘要】本发明涉及一种矿井断层构造预测方法，包括：步骤10，确定并统计已开拓区和未开拓区的沉积介质定量参数、以及已开拓区内各统计单元的煤层小断层定量参数；步骤20，利用灰色关联分析法，对沉积介质定量参数进行筛选，从而得到对煤层小断层发育起决定性作用的主控因素；步骤30，利用主成分分析法对主控因素进行降维处理，以消除主控因素之间的相关性，得到多个主成分；步骤40，将仅包含主成分的数据作为样本数据，利用Elman神经网络建立煤层小断层定量化预测模型；步骤50，根据预测模型，预测未开拓区域的煤层小断层。本发明具有费用低廉、参数易选择、学习速度快、预测精度高的特点。
【专利说明】矿井断层构造预测方法
【技术领域】
[0001]本发明涉及煤矿领域，特别涉及一种矿井断层构造预测方法。
【背景技术】
[0002]煤层小断层是指仅出现于煤层(尤其厚煤层、中厚煤层)、顶板或底板内，其平面分布范围不大、垂向距离有限的小型断层。煤层小断层是影响煤矿生产的重要地质因素之一，它对巷道掘进和工作面回采等具有极大的影响和安全威胁。近年来，由于对地层结构和小断层发育情况了解不清，致使煤矿开采过程中造成了很大的经济损失，甚至引发事故。对煤层小断层进行定量化预测预报具有重要的理论指导意义和实用价值。
[0003]由于煤层小断层的发育具有隐蔽性和有限性，一般不易被探测，预测预报难度较大。自20世纪70年代中期以来，国内外学者对地质构造预测做了大量的研究工作。在定性预测和数理统计预测的基础上，国内学者相继提出了等性块段法、构造力学成因解释法、地质构造规律法、地球物理探测法、几何作图法、灰色模糊法、统计分析法等。这些传统方法对解决矿井地质构造的预测预报难题起到了积极作用，但由于预报的构造太宏观、或多解多译问题难以解决、或探测费用太高、或精度有限等问题，未能很好的预测预报煤层小构造。近年来，相关学者对于小构造的预测高度重视，提出了多种预测方法。
[0004]王生全在2001年第23卷第I期的西安工程学院学报期刊上公开了一种小断层预测的方法，论文名为:煤系变形介质条件研究在煤层小断层预测中的应用。该方法利用灰色关联分析选取了定量化预测参数，建立了基于多元线性回归的小断层定量预测模型。但是由于断层定量参数与沉积介质参数间线性关系较差使得回归模型预测精度较差，且未考虑沉积介质参数之间的相关性，导致信息重叠，加大了预测误差。
[0005]黄乃斌等在2006年第34卷第4期的煤田地质与勘探期刊上公开了一种利用回采巷道的中线水平间距或曲线图预测小构造的方法，论文名为:煤矿开采工作面内小构造预测研究。该方法利用工作面两条回采巷道(风巷、机巷)的实测资料，通过制图预测隐伏构造的出现及其性质，简单实用，但是人工操作误差较大，若工作面内存在小构造的情况预测结果较差，且没有给出小构造的发育程度。
[0006]武强等在2010年第35卷第3期煤炭学报期刊上公开了一种基于环套原理的ANN型矿井小构造预测方法，论文名为:基于环套原理的ANN型矿井小构造预测方法与应用。该文章应用“多重环套理论”和人工神经网络(ANN)技术，提出了预测预报矿井小构造的一套完整理论体系和工作方法，解决了诱发矿井重大灾害事故小构造的预测预报难题。但是该方法只对工作面前方是否存在小构造作出了判断，并未定量预测工作面内小断层的发育程度。
[0007]针对上述问题，本发明提出GRA、PCA和Elman神经网络相结合的小断层定量化预测方法。利用灰色关联分析法选取预测小断层的主控因素，使预测参数易于选择；利用主成分分析法进行主成分建模，消除原始输入层数据的相关性，提取新变量作为Elman神经网络的新输入变量，提高模型的预测速度和精度。
【发明内容】
[0008]为满足煤矿生产和安全的需求，本发明的目的是提供一种费用低廉、参数易选择、学习速度快、预测精度高的矿井断层构造预测方法。
[0009]为解决上述技术问题，作为本发明的一个方面，提供了一种矿井断层构造预测方法，包括:步骤10，确定并统计已开拓区和未开拓区的沉积介质定量参数、以及已开拓区内各统计单元的煤层小断层定量参数；步骤20，利用灰色关联分析法，对沉积介质定量参数进行筛选，从而得到对煤层小断层发育起决定性作用的主控因素；步骤30，利用主成分分析法对主控因素进行降维处理，以消除主控因素之间的相关性，得到多个主成分；步骤40，将仅包含主成分的数据作为样本数据，利用Elman神经网络建立煤层小断层定量化预测模型；步骤50，根据预测模型，预测未开拓区域的煤层小断层。
[0010]本发明选取的沉积介质参数是在勘探或开拓过程中获得的，费用低廉，且将参数定量化，实现了煤层小断层定量化预测；利用灰色关联分析法选取了预测煤层小断层的主控因素，使预测参数易选择；利用主成分分析法进行主成分建模，消除原始输入层数据的相关性，对Elman输入数据进行压缩，减少输入维数，提供了 Elman神经网络的预测精度和速度，具有费用低廉、参数易选择、学习速度快、预测精度高的特点。
【专利附图】
【附图说明】
[0011]图1为本发明方法流程图；
[0012]图2为利用灰色关联分析法对众多沉积介质定量参数进行筛选的流程图；
[0013]图3为利用主成分分析法对主控因素进行降维的流程图；
[0014]图4为本发明实施例已开拓区7煤层小断层密度专题图；
[0015]图5a为本发明实施例7煤层厚度等值线专题图；
[0016]图5b为本发明实施例7煤层直接顶板厚度等值线专题图；
[0017]图5c为本发明实施例7煤层底板厚度等值线专题图；
[0018]图5d为本发明实施例7煤层直接顶板与底板岩性组合类别专题图；
[0019]图5e为本发明实施例7煤层顶板强岩厚度等值线专题图；
[0020]图5f为本发明实施例7煤层底板强岩厚度等值线专题图；
[0021]图5g为本发明实施例7煤层顶底板强岩厚度之和等值线专题图；
[0022]图5h为本发明实施例7煤层顶板强岩层效应值等值线专题图；
[0023]图5i为本发明实施例7煤层底板强岩层效应值等值线专题图；
[0024]图5j为本发明实施例7煤层顶底板强岩层效应总值等值线专题图；
[0025]图6为本发明实施例试验网络的煤层小断层密度预测相对误差图；
[0026]图7为本发明实施例Elman网络训练预测值与实际值对比图；
[0027]图8为本发明实施例煤层小断层定量预测专题图。
【具体实施方式】
[0028]以下对本发明的实施例进行详细说明，但是本发明可以由权利要求限定和覆盖的多种不同方式实施。
[0029]本发明涉及一种矿井断层构造预测方法，特别是一种基于灰色关联(GreyRelat1n Analysis,简称 GRA)-主成分分析(Principal Component Analysis,简称PCA)-Elman神经网络的煤层小断层定量化预测方法。
[0030]请参考图1，本发明提供了一种矿井断层构造预测方法，包括:
[0031]步骤10，确定并统计已开拓区和未开拓区的沉积介质定量参数、以及已开拓区内各统计单元的煤层小断层定量参数；
[0032]步骤20，利用灰色关联分析法，对所述沉积介质定量参数进行筛选，从而得到对煤层小断层发育起决定性作用的主控因素；
[0033]步骤30，利用主成分分析法对所述主控因素进行降维处理，以消除所述主控因素之间的相关性，得到多个主成分；
[0034]步骤40，将仅包含所述主成分的数据作为样本数据，利用Elman神经网络建立煤层小断层定量化预测模型；
[0035]步骤50，根据所述预测模型，预测未开拓区域的煤层小断层。以未开拓区内钻孔统计的介质参数为样本首先进行主成分分析，得到Elman神经网络的输入样本，然后利用前面建立好的模型进行模拟仿真，可得到各钻孔处的煤层小断层定量参数预测值。以开拓区内统计的煤层小断层定量参数值和未开拓区内煤层小断层定量参数预测值，绘制等值线图，即可得到研究区煤层小断层定量预测专题图。
[0036]本发明选取的沉积介质参数是在勘探或开拓过程中获得的，费用低廉，且将参数定量化，实现了煤层小断层定量化预测；利用灰色关联分析法选取了预测煤层小断层的主控因素，使预测参数易选择；利用主成分分析法进行主成分建模，消除原始输入层数据的相关性，对Elman输入数据进行压缩，减少输入维数，提供了 Elman神经网络的预测精度和速度。
[0037]优选地，步骤10中，所述沉积介质定量参数包括:煤层厚度、直接顶板厚度、底板厚度、直接顶板与底板岩性组合类别、顶板强岩厚度、底板强岩厚度、顶底板强岩厚度之和、顶板强岩层效应值、底板强岩层效应值、顶底板强岩层效应总值；其中，将所述直接顶板与底板岩性组合类别按照弱岩层-煤层-弱岩层、弱岩层-煤层-强岩层、强岩层-煤层-强岩层、强岩层-煤层-弱岩层4种类别进行划分，分别定量化赋值依次为1、2、3、4 ;所述煤层小断层定量参数是指煤层小断层密度或煤层小断层强度指数。
[0038]优选地，步骤10还包括统计全区内各钻孔的沉积介质定量参数，并绘制专题图。在采掘工程平面图上，将已开拓区内按照一定长度的网格划分为若干统计单元，统计计算出已采区域内每个单元内的煤层小断层定量参数，并绘制专题图。另外在沉积介质定量参数专题图上，可得到每个单元与煤层小断层定量参数相对应的沉积介质定量参数。
[0039]优选地，所述步骤20包括以下步骤:
[0040]步骤21:按下述方式构建原始数据矩阵:
【权利要求】
1.一种矿井断层构造预测方法，其特征在于，包括: 步骤10，确定并统计已开拓区和未开拓区的沉积介质定量参数、以及已开拓区内各统计单元的煤层小断层定量参数； 步骤20，利用灰色关联分析法，对所述沉积介质定量参数进行筛选，从而得到对煤层小断层发育起决定性作用的主控因素； 步骤30，利用主成分分析法对所述主控因素进行降维处理，以消除所述主控因素之间的相关性，得到多个主成分； 步骤40，将仅包含所述主成分的数据作为样本数据，利用Elman神经网络建立煤层小断层定量化预测模型； 步骤50，根据所述预测模型，预测未开拓区域的煤层小断层。
2.根据权利要求1所述的方法，其特征在于，步骤10中，所述沉积介质定量参数包括:煤层厚度、直接顶板厚度、底板厚度、直接顶板与底板岩性组合类别、顶板强岩厚度、底板强岩厚度、顶底板强岩厚度之和、顶板强岩层效应值、底板强岩层效应值、顶底板强岩层效应总值；其中，将所述直接顶板与底板岩性组合类别按照弱岩层-煤层-弱岩层、弱岩层-煤层-强岩层、强岩层-煤层-强岩层、强岩层-煤层-弱岩层4种类别进行划分，分别定量化赋值依次为1、2、3、4 ; 所述煤层小断层定量参数是指煤层小断层密度或煤层小断层强度指数。
3.根据权利要求1所述的方法，其特征在于，所述步骤20包括以下步骤: 步骤21:按下述方式构建原始数据矩阵:
其中，
为已开拓区内统计的所述煤层小断层定量参数统计数据构成母序列； {χ*01 (/)}为所述沉积介质定量参数统计数据构成子序列；
是对第i个因素在第t个统计单元内得到的观测值，其中，i = O，1，2，…，m，是主因素(i = O)及m个子因素的标号，t = I, 2,…，η为统计单元号； 步骤22:对所述原始数据矩阵进行无量纲化处理，采用如下公式:
其中， η为统计单元个数； 步骤23:计算所述子序列与所述母序列之间的关联度为:
I为分辨系数，一般情况下，取0.1~0.5为宜,通常取0.5 ； 步骤24:将所述关联度大于或等于0.4的沉积介质定量参数作为煤层小断层预测的主控因素。
4.根据权利要求1所述的方法，其特征在于，所述步骤30包括: 步骤31，取已开拓区内的主控因素数据作为样本集； 步骤32，对样本集数据归一化处理，得到样本集矩阵X ； 步骤33，将样本集矩阵X用下式变换为相关矩阵，即主成分矩阵R:
其中:xai为第i个主控因素第a个统计单元的数值； V为第i个主控因素所有统计单元数值的平均值； Xaj为第j个主控因素第a个统计单元的数值； A为第j个主控因素所有统计单元数值的平均值； η为统计单元个数； P为主控因素个数； r.j为第i个主控因素与第j个主控因素的相关系数； 步骤34，确定主成分个数，建立主成分模型； 根据主成分矩阵R求出特征值、主成分贡献率和累计贡献率,确定主成分个数m，并按下式建立主成分模型:
Fi = BliX^a2iX2+---+apiXp (i = I, 2,...，m) 其中，Fi为第i主成分；每个方程中的系数向量(aii，a2i，…，api)是分别是特征值入1;λ2,…，Xm所对应的单位特征向量； 步骤35，求取各单元主成分的取值，作为Elman网络的输入。
5.根据权利要求1所述的方法，其特征在于，所述步骤40包括: 步骤41，将已开拓区内的观测样本划分为训练集和测试集，将确定的主成分作为Elman网络的输入值，将相应的煤层小断层定量参数作为网络的输出值； 步骤42，用训练集样本进行训练，给出训练步长、最大训练步数、均方差最小值、隐层节点数、输入节点数、输出节点数，激励函数，在训练过程中每一步随机设定各层的连接权个阈值，一旦在最大训练步数范围内均方差达到规定值时便结束训练； 步骤43，用测试集样本测试Elman网络，判断煤层小断层预测值与实际值是否满足规定的误差要求，若满足，则终止训练，否则，在重新设定隐层节点，然后再重复以上各步，直到满足规定的误差要求为止。
【文档编号】G06Q10/04GK104200284SQ201410439960
【公开日】2014年12月10日 申请日期:2014年9月1日 优先权日:2014年9月1日
【发明者】邱梅, 施龙青, 牛超, 翟培合 申请人:山东科技大学