CN113947197A - 基于cnn考虑冲击地压前兆信息的微震事件危险性预测方法 - Google Patents

基于cnn考虑冲击地压前兆信息的微震事件危险性预测方法 Download PDF

Info

Publication number
CN113947197A
CN113947197A CN202111558332.8A CN202111558332A CN113947197A CN 113947197 A CN113947197 A CN 113947197A CN 202111558332 A CN202111558332 A CN 202111558332A CN 113947197 A CN113947197 A CN 113947197A
Authority
CN
China
Prior art keywords
microseismic event
microseismic
event
energy
data
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202111558332.8A
Other languages
English (en)
Other versions
CN113947197B (zh
Inventor
李海涛
张海宽
杨冠宇
郑伟钰
杜伟升
薛珊珊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Coal Research Institute CCRI
Original Assignee
China Coal Research Institute CCRI
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Coal Research Institute CCRI filed Critical China Coal Research Institute CCRI
Priority to CN202111558332.8A priority Critical patent/CN113947197B/zh
Publication of CN113947197A publication Critical patent/CN113947197A/zh
Application granted granted Critical
Publication of CN113947197B publication Critical patent/CN113947197B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/04—Architecture, e.g. interconnection topology
    • G06N3/045—Combinations of networks
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/282—Application of seismic models, synthetic seismograms
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/288—Event detection in seismic signals, e.g. microseismics
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/23—Clustering techniques
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/08—Learning methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Remote Sensing (AREA)
  • Business, Economics & Management (AREA)
  • General Engineering & Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Environmental & Geological Engineering (AREA)
  • Evolutionary Computation (AREA)
  • Computational Linguistics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Molecular Biology (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • General Health & Medical Sciences (AREA)
  • Geophysics (AREA)
  • Computing Systems (AREA)
  • Biophysics (AREA)
  • Geology (AREA)
  • Biomedical Technology (AREA)
  • Acoustics & Sound (AREA)
  • Economics (AREA)
  • Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Strategic Management (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Marketing (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Emergency Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Game Theory and Decision Science (AREA)
  • Development Economics (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)

Abstract

本发明涉及一种基于CNN考虑冲击地压前兆信息的微震事件危险性预测方法，首先以时间间隔、各微震事件至中心点距离和能量为特征建立微震事件时空和能量信息预测模型，以时序集中度、时空扩散度、地点集中度和能量为特征指标建立微震事件危险性预测模型，然后利用微震事件时空和能量信息预测模型对未来一次微震事件发生时间、地点和能量数据进行预测，基于预测结果，利用微震事件危险性预测模型预测未来一次微震事件的危险性。解决了目前对于矿井中微震事件危险性预测研究中需要大量不同类型微震监测数据和未考虑冲击地压前兆信息的问题，为冲击地压等矿井中危险事件的防治提供帮助，同时可以为冲击地压等危险性动力灾害发生机理的研究提供帮助。

Description

基于CNN考虑冲击地压前兆信息的微震事件危险性预测方法
技术领域
本发明涉及矿井中冲击地压等危险灾害领域，具体为一种基于CNN（卷积神经网络）考虑冲击地压前兆信息的微震事件危险性预测方法。
背景技术
矿井中冲击地压等危险动力灾害的发生往往造成巷道的毁坏和工作人员的伤亡，因此准确预测矿井中的危险灾害对防治具有十分重要的意义。已有基于微震监测数据对冲击地压等危险灾害预测的方法取得了一定成果，为冲击地压等危险灾害的分析和预测提供了一定依据，但仍存在部分不足。在微震监测数据使用方面，以往研究往往需要微震事件的时序、地点、能量和应力等复杂的监测数据，但实际工程中多数矿井仅有微震的时序、地点、能量数据。在模型设计方面，多数预测或预警模型使用传统机器学习或经验模型，这些模型的相关理论虽然已经成熟且在其他领域取得了一定成果，但难以应对研究复杂、庞大、多变的矿井微震数据的任务。以冲击地压为例，目前关于其发生机理研究中，很多学者总结了煤矿冲击地压发生的很多前兆信息，但大部分预测模型却往往不考虑利用冲击地压发生机理和相关前兆信息构建模型而直接使用简单处理的微震监测信息。
发明内容
针对上述现有技术中存在的不足，本发明提出一种基于CNN考虑冲击地压前兆信息的微震事件危险性预测方法，用于预测未来一次微震事件是否对应冲击地压或地面震动等危险事件，包括如下步骤：
S1，建立微震事件时空和能量信息预测模型
S11，确定微震事件时空和能量信息预测模型的特征指标
以时间间隔、各微震事件发生地点至中心距离和能量作为微震事件时空和能量信息预测模型的特征指标；所述时间间隔为连续多次微震事件中相邻微震事件的时间间隔，其中第一个时间间隔赋值0；所述各微震事件发生地点至中心距离为连续多次微震事件中每次微震发生地点与这多次微震事件中心点之间距离；所述能量为连续多次微震事件中每次微震事件的能量；
S12，建立微震事件时空和能量信息预测模型数据集
在连续N次微震事件中，以第a次至第a+b次微震事件的监测数据作为数据集中一个样本的数据源，1≤a≤N-b，最多可以建立N-b个样本，每个样本包括特征数据和标签数据，所述特征数据包括第a次至第a+b-1次微震事件的时间间隔、第a次至第a+b-1次微震事件中各微震事件发生地点到中心距离和第a次至第a+b-1次微震事件的能量值；所述标签数据为第a+b次微震事件的时间间隔、地点和能量，标签数据需进行One-Hot编码处理；
优选的，步骤S12中，标签数据的第a+b次微震事件的时间间隔为第a+b次微震事件与第a+b-1次微震事件的时间间隔，第a+b次微震事件的地点用三维坐标x、y、z表示。
优选的，步骤S12中，预测之前首先利用Kmeans聚类算法或分段方法对所有微震事件的时间间隔、地点和能量进行分级；标签为One-Hot编码后的分级标签，时间间隔、地点和能量对应级别的标签为1，其它为0。
S13，建立微震事件时空和能量信息预测模型
基于一维卷积神经网络建立一个三层的预测模型，每一层由卷积层、批正则化和激活层构成，以Relu函数为激活函数，预测模型的最后设置两层全连接层用于分类，使用Softmax函数进行标签的输出；
优选的，步骤S13中，三个卷积层的卷积核个数分别为12、24、12，卷积核大小分别为1×1、3×3、3×3，两层全连接层的神经元个数分别为20和10。
优选的，步骤S13中，训练过程中选择Adam优化器来反向传播梯度信息从而更新网络的参数，使用MSELoss损失函数来表示模型计算结果与实际数据之间的差距程度。
优选的，步骤S13中，训练过程中对标签分别为时间间隔、地点坐标和能量的模型进行单独训练。
S2，建立微震事件危险性预测模型
S21，确定微震事件危险性预测模型的特征指标
以时序集中度Q、时空扩散度d、地点集中度D和能量E作为预测模型的特征指标，其中
Figure 492672DEST_PATH_IMAGE002
Figure 940971DEST_PATH_IMAGE004
Figure 501134DEST_PATH_IMAGE006
式中：
Figure 630764DEST_PATH_IMAGE008
分别表示连续n次微震事件中相邻微震事件时间间隔的方差与平均值；
Figure 241874DEST_PATH_IMAGE010
分别表示连续n次微震事件中各微震事件发生地点与连续n次微震事件发生地点中心之间距离的方差与平均值；
S22，建立数据集并进行样本不平衡性处理
在连续M次微震事件中，以第m次至第m+n-1次微震事件的数据作为数据集中一个样本的数据源，1≤m≤M-n+1，最多可以建立M-n+1个样本，每个样本包括特征数据和标签数据，所述特征数据包括第m次至第m+n-1次微震事件的时序集中度Q、时空扩散度d、地点集中度D和第m+n-1次微震事件的能量值；所述标签数据为第m+n-1次微震事件的危险类别，标签数据需进行One-Hot编码处理；样本不平衡性处理，使用上采样法，在采样过程中对数据集中样本标签为危险事件的特征数据加入轻微扰动后复制，以使不同标签的样本数量平衡。
优选的，步骤S22中，One-Hot编码处理后危险事件对应标签为1，不危险事件对应标签为0。
优选的，步骤S22中，轻微扰动的方法为加上随机数。
S23，建立微震事件危险性预测模型
基于一维卷积神经网络建立一个三层的预测模型，每一层由卷积层、批正则化和激活层构成，以Relu函数为激活函数，预测模型的最后设置两层全连接层用于分类，使用Softmax函数进行分类标签的输出；
优选的，步骤S23中，三个卷积层的卷积核个数分别为12、24、12，卷积核大小均为1×1，两层全连接层的神经元个数分别为4和2。
优选的，步骤S23中，训练过程中选择Adam优化器来反向传播梯度信息从而更新网络的参数，使用MSELoss损失函数来表示模型计算结果与实际数据之间的差距程度；
优选的，还包括S24，微震事件危险性预测模型性能验证，通过建立含混合标签样本和单独含危险事件标签样本两个测试集评价微震事件危险性预测模型的性能。
S3，对未来一次微震事件时空和能量信息进行预测
假设未来一次将要发生的微震事件为第d次，之前的d-1次数据为现场监测数据，以包括第d-b次至第d-1次微震事件的时间间隔、第d-b次至第d-1次微震事件的分别发生地点到中心距离和第d-b次至第d-1次微震事件的能量值的特征数据作为输入值，代入步骤S1确定的微震事件时空和能量信息预测模型，输出标签数据，即未来的第d次微震事件的时间间隔、地点坐标和能量；
优选的，每个需要预测的时间间隔、地点坐标和能量特征指标进行单独预测。
S4，预测未来一次将要发生的微震事件的危险性
以包括第d-n+1次至第d次微震事件的时序集中度Q、时空扩散度d、地点集中度D和第d次微震事件的能量值作为输入值，其中第d次微震事件的时间、地点和能量数据为步骤S3的预测值，第d-n+1次至第d-1次微震事件的时间、地点和能量数据为实测值，代入步骤S2确定的微震事件危险性预测模型，输出标签数据，即未来的第d次微震事件是否对应冲击地压等危险事件。
有益效果：本发明以时间间隔、各微震事件至中心点距离和能量为特征建立微震事件时空和能量信息预测模型，以时序集中度、时空扩散度、地点集中度和能量为特征指标建立微震事件危险性预测模型，然后利用微震事件时空和能量信息预测模型对未来一次微震事件发生时间、地点和能量数据进行预测，基于预测结果，利用微震事件危险性预测模型预测未来一次微震事件的危险性。解决了目前对于矿井中微震事件危险性预测研究中需要大量不同类型微震监测数据和未考虑冲击地压前兆信息的问题，为冲击地压等矿井中危险事件的防治提供帮助，同时可以为冲击地压等危险性动力灾害发生机理的研究提供帮助。
附图说明
图1 微震事件时空和能量信息预测模型框架图；
图2 微震事件时空和能量信息预测模型的卷积神经网络模型；
图3（a） 步骤S21确定微震事件危险性预测模型的特征指标；
图3（b） 步骤S22建立数据集并进行样本不平衡性处理；
图3（c） 步骤S23建立微震事件危险性预测模型；
图4 微震事件危险性预测模型的卷积神经网络模型；
图5 利用某矿井微震监测数据评价微震事件危险性预测模型性能的实例。
具体实施方式
下面结合附图，对本发明的技术方案进行更为详细的描述。
如图1-5所示，基于CNN考虑冲击地压前兆信息的微震事件危险性预测方法，用于预测将要发生的一次（未来一次）微震事件是否对应冲击地压或地面震动等危险事件，包括如下步骤：
S1，建立微震事件时空和能量信息预测模型
S11，确定微震事件时空和能量信息预测模型的特征指标
以时间间隔（连续多次微震事件中相邻微震事件的时间间隔，第一个时间间隔赋值0）、各微震事件发生地点至中心距离（连续多次微震事件中每次微震发生地点与这多次微震事件中心点之间距离）和能量（连续多次微震事件中每次微震事件的能量）作为微震事件时空和能量信息预测模型的特征指标；
S12，建立微震事件时空和能量信息预测模型数据集
在连续N次微震事件中，以第a次至第a+b次（连续b+1次）微震事件的监测数据作为数据集中一个样本的数据源，1≤a≤N-b，则最多可以建立N-b个样本，每个样本包括特征数据和标签数据，其中特征数据包括第a次至第a+b-1次微震事件的时间间隔（两两之间时间间隔、第一个时间间隔赋值0）、第a次至第a+b-1次微震事件中各微震事件发生地点到中心距离（连续b次微震事件的分别发生地点与这b次微震事件发生地点中心之间距离）和第a次至第a+b-1次微震事件的能量值；所述标签数据为第a+b次微震事件的时间间隔（第a+b次微震事件与第a+b-1次微震事件的时间间隔）、地点（第a+b次微震事件的地点，可以用三维坐标x、y、z表示）和能量（第a+b次微震事件的能量）；标签数据需进行One-Hot编码处理。
例如，在连续N=7次微震事件中，以第a次至第a+5次（N=7，b=5）微震事件的监测数据作为数据集中一个样本的数据源，1≤a≤7-5=2，则最多可以建立2条数据；第1个样本的特征数据包括第1次至第5次微震事件的两两时间间隔（第一个时间间隔赋值0）、第1次至第5次微震事件中各微震事件发生地点到中心距离和第1次至第5次微震事件的能量值，所述标签数据为第6次微震事件与第5次微震事件的时间间隔、第6次微震事件地点和第6次微震事件的能量值；第2个样本的特征数据包括第2次至第6次微震事件的两两时间间隔（第一个时间间隔赋值0）、第2次至第6次微震事件分别至中心距离和第2次至第6次微震事件的能量值，所述标签数据为第7次微震事件与第6次微震事件的时间间隔、第7次微震事件地点和第7次微震事件的能量值。
优选的，预测之前首先利用Kmeans聚类算法或分段方法对所有微震事件（包括特征数据和标签数据）的时间间隔、各微震事件发生地点至中心距离和能量进行分类或分段（即分级），以简化运算；标签为One-Hot编码后的分级标签，即时间间隔、地点坐标和能量对应级别的标签为1，其它为0。
以能量分级为例，假设能量值的范围为0-10000，可以分为10个级别，每个级别范围为1000，即可以以0-1000（含1000）为第1级，1000-2000（不含1000，含2000）为第2级，以此类推，第10级为9000-10000（不含9000，含10000）；某个样本的能量为7200时，其对应的能量级别为7，则在One-Hot编码后其标签采用如表1所示；
表1 能量分级及One-Hot编码示例
能量级别 1 2 3 4 5 6 7 8 9 10
One-Hot编码 0 0 0 0 0 0 1 0 0 0
S13，建立微震事件时空和能量信息预测模型
如图1-2所示，基于一维卷积神经网络建立一个三层的预测模型，每一层由卷积层、批正则化（Batch Normalization）和激活层构成，三个卷积层的卷积核个数分别为12、24、12，卷积核大小分别为1×1、3×3、3×3，以Relu函数为激活函数，预测模型的最后设置两层全连接层用于分类，神经元个数分别为20和10，最后使用Softmax函数进行分类标签的输出；训练过程中选择Adam优化器来反向传播梯度信息从而更新网络的参数，使用MSELoss损失函数来表示模型计算结果与实际数据之间的差距程度；
优选的，步骤S13中，训练过程中对标签分别为时间间隔、地点坐标和能量的模型进行单独训练。
S2，建立微震事件危险性预测模型
S21，确定微震事件危险性预测模型的特征指标
如图3（a）所示，建立指标过程中主要参考冲击地压的前兆信息，冲击地压发生前由于围岩系统与外界能量的交换，围岩结构处于非稳定调整期，因此在冲击地压发生前往往出现一段微震事件的活跃期，该微震事件活跃期的特征可以通过时空指标（时序集中度、时空扩散度、地点集中度）体现，此外能量也是判断地压等危险事件是否发生的重要指标。
据此确定，以时序集中度Q、时空扩散度d、地点集中度D和能量E作为预测模型的特征指标，其中
Figure 595495DEST_PATH_IMAGE002
Figure 596949DEST_PATH_IMAGE004
Figure 327008DEST_PATH_IMAGE006
式中：
Figure 543357DEST_PATH_IMAGE008
分别表示连续n次微震事件中相邻微震事件时间间隔的方差与平均值；
Figure 271141DEST_PATH_IMAGE010
分别表示连续n次微震事件中各微震事件发生地点与连续n次微震事件发生地点中心之间距离的方差与平均值；
S22，建立数据集并进行样本不平衡性处理
如图3（b）所示，在连续M次微震事件中，以第m次至第m+n-1次（连续n次）微震事件的数据作为数据集中一个样本的数据源（在模型训练过程中使用第m+n次微震事件的监测数据，在模型预测过程中使用经过S1中微震事件时空和能量信息模型预测得到的第m+n次微震事件的数据），1≤m≤M-n+1，则最多可以建立M-n+1个样本，每个样本包括特征数据和标签数据，其中特征数据包括第m次至第m+n-1次微震事件的时序集中度Q、时空扩散度d、地点集中度D和第m+n-1次微震事件的能量值；所述标签数据为第m+n-1次微震事件的危险类别，即标签数据为是否发生了危险事件（冲击地压或地面震动），标签数据需进行One-Hot编码处理，处理后危险事件对应标签为1，不危险事件对应标签为0；
例如，在连续7次微震事件中，以第m次至第m+5-1次微震事件的监测数据作为数据集中一个样本的数据源（M=7，n=5），1≤m≤7-5+1=3，则最多可以建立3条数据；第1个样本的特征数据包括第1次至第5次微震事件的时序集中度Q、时空扩散度d、地点集中度D和第5次微震事件的能量值，所述标签数据为第5次微震事件的危险类别；第2个样本的特征数据包括第2次至第6次微震事件的时序集中度Q、时空扩散度d、地点集中度D和第6次微震事件的能量值，所述标签数据为第6次微震事件的危险类别；第3个样本的特征数据包括第3次至第7次微震事件的时序集中度Q、时空扩散度d、地点集中度D和第7次微震事件的能量值，所述标签数据为第7次微震事件的危险类别。
样本不平衡性处理，矿井微震事件中危险事件数量要远远小于不危险事件数量，故使用上采样法，在采样过程中对数据集中样本标签为危险事件的特征数据加入轻微扰动（如加上随机数）后复制，以使不同标签（危险类别，包括危险事件和不危险事件）的样本数量平衡；
S23，建立微震事件危险性预测模型
如图3（c）和图4所示，基于一维卷积神经网络建立一个三层的预测模型，每一层由卷积层、批正则化（Batch Normalization）和激活层构成，卷积层的作用是信息融合和提取，三个卷积层的卷积核个数分别为12、24、12，卷积核大小均为1×1，以Relu函数为激活函数，预测模型的最后设置两层全连接层用于分类，神经元个数分别为4和2，最后使用Softmax函数进行分类标签的输出；训练过程中选择Adam优化器来反向传播梯度信息从而更新网络的参数，使用MSELoss损失函数来表示模型计算结果与实际数据之间的差距程度；
S24，微震事件危险性预测模型性能验证
考虑到对冲击地压等危险事件预测正确的重要性，通过建立含混合标签样本和单独含危险事件标签样本两个测试集评价微震事件危险性预测模型的性能；
如图5所示，基于某矿井微震事件监测数据对微震事件危险性预测模型的性能进行评价，随机选取训练样本进行了20次实验，每次实验中微震危险性预测模型在含混合标签样本测试集和单独含危险事件标签样本测试集上同时达到的最高准确率，最高准确率分别为0.913和0.838。
S3，对未来一次微震事件时空和能量信息进行预测
假设未来一次将要发生的微震事件为第d次（之前的d-1次数据为现场监测数据），以包括第d-b次至第d-1次微震事件（连续b次）的时间间隔（两两之间时间间隔、第一个时间间隔赋值0）、第d-b次至第d-1次微震事件的分别发生地点到中心距离（连续b次微震事件的分别发生地点与这b次微震事件发生地点中心之间距离）和第d-b次至第d-1次微震事件的能量值的特征数据作为输入值，代入步骤S1确定的微震事件时空和能量信息预测模型，输出标签数据（每种信息单独预测），即未来一次将要发生（第d次）的微震事件的时间间隔（第d次微震事件与第d-1次微震事件的时间间隔）、地点（可以用三维坐标x、y、z表示）和能量；
S4，预测未来一次将要发生的微震事件的危险性
以包括第d-n+1次至第d次微震事件（连续n次）的时序集中度Q、时空扩散度d、地点集中度D和第d次微震事件的能量值作为输入值，其中第d次微震事件的时间、地点和能量数据为步骤S3的预测值，第d-n+1次至第d-1次微震事件的时间、地点和能量数据为实测值，代入步骤S2确定的微震事件危险性预测模型，输出标签数据，即未来一次将要发生的（第d次）微震事件是否对应冲击地压等危险事件。

Claims (9)

1.基于CNN考虑冲击地压前兆信息的微震事件危险性预测方法，其特征在于，包括如下步骤：
S1，建立微震事件时空和能量信息预测模型
S11，确定微震事件时空和能量信息预测模型的特征指标
以时间间隔、各微震事件发生地点至中心距离和能量作为特征指标；所述时间间隔为连续多次微震事件中相邻微震事件的时间间隔，其中第一个时间间隔赋值0；所述各微震事件发生地点至中心距离为连续多次微震事件中每次微震发生地点与这多次微震事件中心点之间距离；所述能量为连续多次微震事件中每次微震事件的能量；
S12，建立微震事件时空和能量信息预测模型数据集
在连续N次微震事件中，以第a次至第a+b次微震事件的监测数据作为数据集中一个样本的数据源，1≤a≤N-b，最多可以建立N-b个样本，每个样本包括特征数据和标签数据，所述特征数据包括第a次至第a+b-1次微震事件的时间间隔、第a次至第a+b-1次微震事件中各微震事件发生地点到中心距离和第a次至第a+b-1次微震事件的能量值；所述标签数据为第a+b次微震事件的时间间隔、地点和能量，标签数据需进行One-Hot编码处理；
S13，建立微震事件时空和能量信息预测模型
基于一维卷积神经网络建立一个三层的预测模型，每一层由卷积层、批正则化和激活层构成，以Relu函数为激活函数，预测模型的最后设置两层全连接层用于分类，使用Softmax函数进行标签的输出；
S2，建立微震事件危险性预测模型
S21，确定微震事件危险性预测模型的特征指标
以时序集中度Q、时空扩散度d、地点集中度D和能量E作为特征指标，其中
Figure 347016DEST_PATH_IMAGE001
Figure 401560DEST_PATH_IMAGE002
Figure 611830DEST_PATH_IMAGE003
式中：
Figure 384614DEST_PATH_IMAGE004
分别表示连续n次微震事件中相邻微震事件时间间隔的方差与平均值；
Figure 190896DEST_PATH_IMAGE005
分别表示连续n次微震事件中各微震事件发生地点与连续n次微震事件发生地点中心之间距离的方差与平均值；
S22，建立数据集并进行样本不平衡性处理
在连续M次微震事件中，以第m次至第m+n-1次微震事件的数据作为数据集中一个样本的数据源，1≤m≤M-n+1，最多可以建立M-n+1个样本，每个样本包括特征数据和标签数据，所述特征数据包括第m次至第m+n-1次微震事件的时序集中度Q、时空扩散度d、地点集中度D和第m+n-1次微震事件的能量值；所述标签数据为第m+n-1次微震事件的危险类别，标签数据需进行One-Hot编码处理；样本不平衡性处理，在采样过程中对数据集中样本标签为危险事件的特征数据加入轻微扰动后复制，以使不同标签的样本数量平衡；
S23，建立微震事件危险性预测模型
基于一维卷积神经网络建立一个三层的预测模型，每一层由卷积层、批正则化和激活层构成，以Relu函数为激活函数，预测模型的最后设置两层全连接层用于分类，使用Softmax函数进行分类标签的输出；
S3，对未来一次微震事件时空和能量信息进行预测
假设未来一次微震事件为第d次，之前的d-1次数据为现场监测数据，以包括第d-b次至第d-1次微震事件的时间间隔、第d-b次至第d-1次微震事件的分别发生地点到中心距离和第d-b次至第d-1次微震事件的能量值的特征数据作为输入值，代入步骤S1确定的微震事件时空和能量信息预测模型，输出标签数据，即未来的第d次微震事件的时间间隔、地点坐标和能量；
S4，预测未来一次将要发生的微震事件的危险性
以包括第d-n+1次至第d次微震事件的时序集中度Q、时空扩散度d、地点集中度D和第d次微震事件的能量值作为输入值，其中第d次微震事件的时间、地点和能量数据为步骤S3的预测值，第d-n+1次至第d-1次微震事件的时间、地点和能量数据为实测值，代入步骤S2确定的微震事件危险性预测模型，输出标签数据，即未来的第d次微震事件是否对应冲击地压危险事件。
2.根据权利要求1所述的微震事件危险性预测方法，其特征在于，步骤S12中，标签数据的第a+b次微震事件的时间间隔为第a+b次微震事件与第a+b-1次微震事件的时间间隔，第a+b次微震事件的地点用三维坐标x、y、z表示。
3.根据权利要求1或2所述的微震事件危险性预测方法，其特征在于，步骤S12中，预测之前首先利用Kmeans聚类算法或分段方法对所有微震事件的时间间隔、各微震事件发生地点至中心距离和能量进行分级；标签为One-Hot编码后的分级标签，时间间隔、地点和能量对应级别的标签为1，其它为0。
4.根据权利要求1所述的微震事件危险性预测方法，其特征在于，步骤S13中，三个卷积层的卷积核个数分别为12、24、12，卷积核大小分别为1×1、3×3、3×3，两层全连接层的神经元个数分别为20和10；和/或步骤S23中，三个卷积层的卷积核个数分别为12、24、12，卷积核大小均为1×1，两层全连接层的神经元个数分别为4和2。
5.根据权利要求4所述的微震事件危险性预测方法，其特征在于，步骤S13和S23中，训练过程中选择Adam优化器来反向传播梯度信息从而更新网络的参数，使用MSELoss损失函数来表示模型计算结果与实际数据之间的差距程度。
6.根据权利要求1或4所述的微震事件危险性预测方法，其特征在于，步骤S13中，训练过程中对标签分别为时间间隔、地点坐标和能量的模型进行单独训练；和/或步骤S3中，每个需要预测的时间间隔、地点坐标和能量特征指标进行单独预测。
7.根据权利要求1所述的微震事件危险性预测方法，其特征在于，步骤S22中，One-Hot编码处理后危险事件对应标签为1，不危险事件对应标签为0。
8.根据权利要求1或7所述的微震事件危险性预测方法，其特征在于，步骤S22中，轻微扰动的方法为加上随机数。
9.根据权利要求1所述的微震事件危险性预测方法，其特征在于，还包括S24，微震事件危险性预测模型性能验证，通过建立含混合标签样本和单独含危险事件标签样本两个测试集评价微震事件危险性预测模型的性能。
CN202111558332.8A 2021-12-20 2021-12-20 基于cnn考虑冲击地压前兆信息的微震事件危险性预测方法 Active CN113947197B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111558332.8A CN113947197B (zh) 2021-12-20 2021-12-20 基于cnn考虑冲击地压前兆信息的微震事件危险性预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111558332.8A CN113947197B (zh) 2021-12-20 2021-12-20 基于cnn考虑冲击地压前兆信息的微震事件危险性预测方法

Publications (2)

Publication Number Publication Date
CN113947197A true CN113947197A (zh) 2022-01-18
CN113947197B CN113947197B (zh) 2022-04-26

Family

ID=79339353

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111558332.8A Active CN113947197B (zh) 2021-12-20 2021-12-20 基于cnn考虑冲击地压前兆信息的微震事件危险性预测方法

Country Status (1)

Country Link
CN (1) CN113947197B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114519392A (zh) * 2022-02-14 2022-05-20 煤炭科学研究总院有限公司 冲击地压的预警方法及装置、电子设备和存储介质
CN114757266A (zh) * 2022-03-27 2022-07-15 中国矿业大学 专家知识与数据融合驱动的冲击地压预测模型构建方法
CN114994760A (zh) * 2022-08-08 2022-09-02 煤炭科学研究总院有限公司 基于Transformer和自监督算法的微震事件危险性分析方法
CN114997742A (zh) * 2022-07-19 2022-09-02 煤炭科学研究总院有限公司 基于注意力机制和聚类算法的微震事件危险性分析方法
CN116500678A (zh) * 2023-06-26 2023-07-28 煤炭科学研究总院有限公司 一种多源数据因果融合驱动的微震时空预测方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106680867A (zh) * 2016-11-17 2017-05-17 大连理工大学 一种微震事件精确定位的动态参数方法
CN107203003A (zh) * 2017-07-27 2017-09-26 河北煤炭科学研究院 一种矿井水害微震监测时空簇分析方法
EP3619560A1 (en) * 2017-05-05 2020-03-11 ConocoPhillips Company Stimulated rock volume analysis
CN111123355A (zh) * 2020-01-07 2020-05-08 山东大学 一种基于微震监测数据的岩爆预测方法及系统
CN113723595A (zh) * 2021-09-02 2021-11-30 重庆大学 基于微震事件定量预测的煤矿冲击地压智能预警方法
CN113756870A (zh) * 2021-09-09 2021-12-07 中煤科工开采研究院有限公司 冲击地压微震监测系统布置方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106680867A (zh) * 2016-11-17 2017-05-17 大连理工大学 一种微震事件精确定位的动态参数方法
EP3619560A1 (en) * 2017-05-05 2020-03-11 ConocoPhillips Company Stimulated rock volume analysis
CN107203003A (zh) * 2017-07-27 2017-09-26 河北煤炭科学研究院 一种矿井水害微震监测时空簇分析方法
CN111123355A (zh) * 2020-01-07 2020-05-08 山东大学 一种基于微震监测数据的岩爆预测方法及系统
CN113723595A (zh) * 2021-09-02 2021-11-30 重庆大学 基于微震事件定量预测的煤矿冲击地压智能预警方法
CN113756870A (zh) * 2021-09-09 2021-12-07 中煤科工开采研究院有限公司 冲击地压微震监测系统布置方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
JOZEF KABIESZ等: "神经网络预测矿震危险性", 《煤炭科技》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114519392A (zh) * 2022-02-14 2022-05-20 煤炭科学研究总院有限公司 冲击地压的预警方法及装置、电子设备和存储介质
CN114519392B (zh) * 2022-02-14 2023-05-23 煤炭科学研究总院有限公司 一种基于微震时序数据的冲击地压危险预测算法
CN114757266A (zh) * 2022-03-27 2022-07-15 中国矿业大学 专家知识与数据融合驱动的冲击地压预测模型构建方法
CN114997742A (zh) * 2022-07-19 2022-09-02 煤炭科学研究总院有限公司 基于注意力机制和聚类算法的微震事件危险性分析方法
CN114994760A (zh) * 2022-08-08 2022-09-02 煤炭科学研究总院有限公司 基于Transformer和自监督算法的微震事件危险性分析方法
CN114994760B (zh) * 2022-08-08 2022-10-14 煤炭科学研究总院有限公司 基于Transformer和自监督算法的微震事件危险性分析方法
CN116500678A (zh) * 2023-06-26 2023-07-28 煤炭科学研究总院有限公司 一种多源数据因果融合驱动的微震时空预测方法
CN116500678B (zh) * 2023-06-26 2023-09-26 煤炭科学研究总院有限公司 一种多源数据因果融合驱动的微震时空预测方法

Also Published As

Publication number Publication date
CN113947197B (zh) 2022-04-26

Similar Documents

Publication Publication Date Title
CN113947197B (zh) 基于cnn考虑冲击地压前兆信息的微震事件危险性预测方法
CN113255815B (zh) 用户行为异常分析方法、装置、设备及存储介质
US11709979B1 (en) Bridge damage identification method considering uncertainty
CN112039903B (zh) 基于深度自编码神经网络模型的网络安全态势评估方法
CN104636449A (zh) 基于lsa-gcc的分布式大数据系统风险识别方法
CN110119853B (zh) 基于时间序列监测数据分析的供水管网漏损报警阈值选取方法
CN114037022A (zh) 一种考虑监测数据时空特性的微震事件危险性分类方法
CN114266289A (zh) 一种复杂装备健康状态评估方法
CN116341901B (zh) 一种滑坡面域-单体危险性预警的一体化评价方法
CN116307103A (zh) 一种基于硬参数共享多任务学习的交通事故预测方法
Qin et al. Evaluation of goaf stability based on transfer learning theory of artificial intelligence
CN114662405A (zh) 基于少样本度量和集成学习的岩爆预测方法
Kalogirou et al. Package ‘lctools’
CN111737993B (zh) 一种配电网设备的故障缺陷文本提取设备健康状态方法
CN109656904A (zh) 一种案件风险检测方法及系统
Nanda et al. A tri-stage cluster identification model for accurate analysis of seismic catalogs
Fuentes et al. Spatio-temporal seismicity prediction in Chile using a multi-column ConvLSTM
CN111582313A (zh) 样本数据生成方法、装置及电子设备
CN114553468A (zh) 一种基于特征交叉与集成学习的三级网络入侵检测方法
CN114997742B (zh) 基于注意力机制和聚类算法的微震事件危险性分析方法
Senkevich Search for anomalies in pulse flows of acoustic and electromagnetic emissions
CN115982968B (zh) 一种揭示地面时空间演化差异的算法
CN110288140B (zh) 一种基于地理相关性模型的阿片类药物空间传播预测方法
Piccarreta et al. Assessing and comparing models for sequence data by microsimulation (with Supplementary Material)
Sadeghian The effects of zoning methods on forecasting of the next earthquake occurrences, according to semi-Markov models

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information

Address after: 100010 No. 5, qingqinggou Road, Chaoyang District, Beijing

Applicant after: General Coal Research Institute Co.,Ltd.

Address before: 100010 No. 5, qingqinggou Road, Chaoyang District, Beijing

Applicant before: CHINA COAL Research Institute

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant