CN107563092B - 一种矿山动力灾害的全息预警方法 - Google Patents
一种矿山动力灾害的全息预警方法 Download PDFInfo
- Publication number
- CN107563092B CN107563092B CN201710851433.1A CN201710851433A CN107563092B CN 107563092 B CN107563092 B CN 107563092B CN 201710851433 A CN201710851433 A CN 201710851433A CN 107563092 B CN107563092 B CN 107563092B
- Authority
- CN
- China
- Prior art keywords
- coal
- gas
- mine
- early warning
- holographic
- 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.)
- Active
Links
Images
Landscapes
- Alarm Systems (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种矿山动力灾害的全息预警方法,以直接影响矿山动力灾害发生与否的模态化参数集为输入变量集,以冲击地压发生的可能性和煤与瓦斯突出发生的可能性作为输出变量,建立人工神经网络模型;向学习样本库中添加样本,得到一个矿山动力灾害的全息模态化预警器;对矿山任意一个区域进行预警时,利用该全息数据挖掘转换器把矿山实验测试和安全生产中获得的相关信息转化成预警模型的输入参数,再利用全息模态化预警器给出预警结果。全息预警方法可以实现矿山动力灾害的全息模态化在线预测预报和预警,预警时系统可根据冲击地压事故发生的可能性和煤与瓦斯突出事故发生的可能性分别实现预警。
Description
技术领域
本发明涉及矿山和地下工程的灾害预警和防治技术领域,尤其涉及一种矿山动力灾害的全息预警方法。
背景技术
矿山动力灾害,包括冲击地压(岩爆)、煤与瓦斯突出、突水透水等都属于矿山重大灾害,一旦发生这类灾害对矿山的人员生命和财产都将造成重大损失。因此,对矿山动力灾害的超前预测、预报和预警就显得致关重要。由于矿山动力灾害涉及的因素众多,至今对这类灾害的发生机理、演化过程和诱发方式没有形成一套完整的理论体系,也没有形成可用于可靠监测的、编程实现的和准确率较高的预警方法。现有的大多数预测和预警方法大都是采用层次分析法的思想,通过各层次预警参数权值的定义,实现综合加权法进行灾害预警。但是,这些综合加权法存在一些不足之处,比如所选的预警参数大部分是间接参数,而不是直接参数,例如采深、地质构造、生产工艺、通风方式等;而且很多间接参数都和某一种直接参数有关,而且不只是线性关系,导致权值难以确定。现有方法即使考虑了一些直接参数,例如,冲击倾向性、突出倾向性等,也没有给出这些指标的合理计算方法,例如国家标准GB/T 25217.2-2010,利用加权综合评价方法,有8种情况无法判别都无法给出冲击倾向性结果。把直接指标和间接指标放到一个表达式中进行加权,无法体现各种指标的层次关系,以及对预警结果的支持度。综合加权法无法体现迹象越多,发生事故越大的特点,一旦确定了预警参数集,新的迹象和显现将无法加入,一旦加入,就要重新分配权值,导致预警结果既不确定,也不科学。
发明内容
为了克服上述现有技术中的不足,本发明提供一种矿山动力灾害的全息预警方法,方法包括:
步骤一,以直接影响矿山动力灾害发生与否的模态化参数集为输入变量集,以冲击地压发生的可能性和煤与瓦斯突出发生的可能性作为输出变量,建立人工神经网络模型;
步骤二,向学习样本库中添加样本,当学习样本库中添加的样本达到预设值时,运用预设的学习算法对学习样本库中的样本进行学习,以确定人工神经网络模型中的权值和阈值,得到一个矿山动力灾害的全息模态化预警器;
步骤三,根据全息模态化预警器设置一包括空间检索、插补、降噪、反飘移和信号重建功能的全息数据挖掘转换器;
对矿山任意一个区域进行预警时,利用该的全息数据挖掘转换器把矿山实验测试和安全生产中获得的相关信息转化成预警模型的输入参数,再利用全息模态化预警器给出预警结果。
优选地,矿山实验测试和安全生产中获得的相关信息包括:矿山地理信息、地质信息、生产信息和相关属性信息;
相关属性信息包括:地形、岩层、断层、陷落柱、褶皱、岩浆侵入、火烧等构造信息,地温、水压、岩移、沉陷、煤柱、井巷、硐室、工作面、采空区空间分布信息,支护、卸压、充填、瓦斯抽放灾害防治工艺等参数,以及冲击倾向性测试数据,防突测试数据,微震信号、地音信号、爆破信号、电磁辐射信号、应力信号、矿压信号、锚杆阻力、离层和位移信号、瓦斯浓度、风速、瓦斯抽放流量、水位、涌水量监测数据。
优选地,步骤三中,利用该的全息数据挖掘转换器把矿山实验测试和安全生产中获得的相关信息转化成预警模型的输入参数还包括:
利用岩性测试和煤体测试得到的单轴抗压强度Rc、煤层的冲击能量指数KE、煤层的单性能量指数WET,煤层的动态破坏时间DT,用下列方法计算出围岩和煤层冲击倾向性X24;
P1=PDT+PRC-PDT×PRC
P2=PWET+PKE-PWET×PKE
X24=0.7×P1+0.5×P2-0.2×P1×P2
(2)利用煤层的瓦斯试验数据,例如瓦斯放散初速度Δp、瓦斯涌出初速度q1、解吸指标Δh2、煤层的硬度系数f、综合指标K=Δp/f、k1钻屑瓦斯解吸指标、每个钻孔沿孔长的最大钻屑量Smax、每个钻孔的最大钻孔涌出初速度qmax、钻孔平均综合指标R=(Smax-1.8)×(qmax-4),用下列方法计算出煤层瓦斯突出倾向性X14;
p1=pΔp+pq-pΔp×pq
p2=pΔh2+pk1-pΔh2×pk1
p3=pK+pR-pK×pR
p4=p1+p2-p1×p2
X14=0.7×p3+0.5×p4-0.2×p3×p4
(3)利用微震定位结果,通过反演岩层或煤层的地震波波速vp和横波波速vs,利用式即可解出岩层或煤层的当前状态对应的围岩或煤层弹性模量E(即X26)和泊松比μ;
其中,ρ为岩层或煤层的密度,是已知量;
(4)利用岩层或煤层的当前状态对应的弹性模量E和原始弹性模量E0,即可求出损伤变量D;
(5)利用损伤变量D和原始抗压强度S0,即可求当前的围岩或煤层强度S(即X4);
S=S0(1-D)
(6)利用三轴应力监测系统直接监测测点的绝对地应力(σx,σy,σz),即(X1,X2,X3),也可以先根据三维空间模型,利用解析法计算出初步的应力分布,在用包括但不限于应力监测、矿压监测的实时数据进行修正,得到(X1,X2,X3);
(7)工作面顶板压力Q(即X6)和围岩变形量U(X7),可以用支架阻力传感器和位移传感器直接测量;
(8)根据地应力(X1,X2,X3)和井巷空间模型,就可计算出围岩或煤体抵抗线长W(即X5);
(9)可用下列公式计算围岩体积变形能UV(X8)、围岩形状变形能Uf(X9)和工作面顶板弯曲能Uw(X10);
其中的l,h和b1分别为工作面顶板的长度,厚度和宽度,q2顶板受均匀载荷,E为弹性模量;
G为剪切模量,可由下式计;
(10)利用下列公式计算煤层瓦斯含量Vt(即X11);
Vt=Qf+Qc+Vc
其中Qf为流量涌出量,可用巷道断面、风速监测值和瓦斯浓度监测值计算;
Qc为瓦斯抽放流量,可用瓦斯抽放系统的纯瓦斯流量监测值得到;
Vc为煤体瓦斯残存含量,可用地面试验获得;
(11)利用下列公式计算煤层瓦斯压力P(即X12);
其中:η为煤层孔隙度、a和b2为瓦斯吸附常数,可通过试验获得;
p0为标准大气压;
(12)区域自由空间中绝对气压(即X25)可直接监测;
(13)对应X1的地应力增量X15、对应X2的地应力增量X16、对应X3的地应力增量X17、围岩或煤体强度增量X18、围岩或煤体抵抗线长增量X19、工作面顶板压力增量X20、围岩变形量增量X21、煤层瓦斯含量增量X22、煤层瓦斯压力增量X23可由当前时段的状态量减去前一时段的相应状态量得到。
优选地,步骤一还包括:
选取预警区域和地点的动力灾害公用参数X方向地应力X1,Y方向地应力X2,Z方向地应力X3,围岩或煤层强度X4,围岩或煤体抵抗线长X5,工作面顶板压力X6,围岩变形量X7,围岩体积变形能X8,围岩形状变形能X9以及工作面顶板弯曲能X10,与瓦斯突出强相关的煤层瓦斯含量X11、煤层瓦斯压力X12,当前状态的持续时间X13,煤层瓦斯突出倾向性X14,对应X1的地应力增量X15、对应X2的地应力增量X16、对应X3的地应力增量X17、围岩或煤体强度增量X18、围岩或煤体抵抗线长增量X19、工作面顶板压力增量X20、围岩变形量增量X21、煤层瓦斯含量增量X22、煤层瓦斯压力增量X23、围岩和煤层冲击倾向性X24、区域自由空间绝对气压X25、围岩或煤层弹性模量X26,为26个输入变量;选取冲击地压发生的可能性Y1和煤与瓦斯突出发生的可能性Y2,为2个为输出变量;设置一个单隐含层,隐含层的节点个数L,建立矿山动力灾害预警的神经网络模型。
从以上技术方案可以看出,本发明具有以下优点:
本发明包括实验测试数据、地理空间数据、安全生产数据、监测监控数据在内的矿山全息数据库,可以实现矿山动力灾害的全息模态化在线预测预报和预警,预警时系统可根据冲击地压事故发生的可能性和煤与瓦斯突出事故发生的可能性分别实现预警,并可把预警信息发布到各种显示终端。
本发明可以实现矿山动力灾害的全息模态化在线预测预报和预警,以其达到矿山动力灾害的有效预警和提前防治。
附图说明
为了更清楚地说明本发明的技术方案,下面将对描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为矿山动力灾害的全息预警方法流程图;
图2为矿山动力灾害预警的神经网络模型图;
图3为矿山动力灾害的全息模态化预警方法及系统架构图。
具体实施方式
为使得本发明的发明目的、特征、优点能够更加的明显和易懂,下面将运用具体的实施例及附图,对本发明保护的技术方案进行清楚、完整地描述,显然,下面所描述的实施例仅仅是本发明一部分实施例,而非全部的实施例。基于本专利中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本专利保护的范围。
需要说明的是:本方案中X1-X26可以用X1-X26,它们在本方案中代表的含义相同,其他存在上述情况的也是一样,本领域技术人员根据所具有的知识能够清楚的理解,本方案中出现此种情况是为了便于描述记载,不存在异议问题。
本实施例提供一种矿山动力灾害的全息预警方法,如图1和图3所示,方法包括:
步骤一,以直接影响矿山动力灾害发生与否的模态化参数集为输入变量集,以冲击地压发生的可能性和煤与瓦斯突出发生的可能性作为输出变量,建立人工神经网络模型;
步骤二,向学习样本库中添加样本,当学习样本库中添加的样本达到预设值时,运用预设的学习算法对学习样本库中的样本进行学习,以确定人工神经网络模型中的权值和阈值,得到一个矿山动力灾害的全息模态化预警器;
步骤三,根据全息模态化预警器设置一包括空间检索、插补、降噪、反飘移和信号重建功能的全息数据挖掘转换器;
对矿山任意一个区域进行预警时,利用该的全息数据挖掘转换器把矿山实验测试和安全生产中获得的相关信息转化成预警模型的输入参数,再利用全息模态化预警器给出预警结果。
其中,矿山实验测试和安全生产中获得的相关信息包括:矿山地理信息、地质信息、生产信息和相关属性信息;相关属性信息包括:地形、岩层、断层、陷落柱、褶皱、岩浆侵入、火烧等构造信息,地温、水压、岩移、沉陷、煤柱、井巷、硐室、工作面、采空区空间分布信息,支护、卸压、充填、瓦斯抽放灾害防治工艺等参数,以及冲击倾向性测试数据,防突测试数据,微震信号、地音信号、爆破信号、电磁辐射信号、应力信号、矿压信号、锚杆阻力、离层和位移信号、瓦斯浓度、风速、瓦斯抽放流量、水位、涌水量监测数据。
本实施例中,步骤一,选取直接影响矿山动力灾害发生与否的模态化参数集为输入变量集,以冲击地压发生的可能性和煤与瓦斯突出发生的可能性作为输出变量,建立人工神经网络模型,建立人工神经网络模型。
以冲击地压和煤与瓦斯突出灾害为例,可选取预警区域和地点的动力灾害公用参数X方向地应力X1、Y方向地应力X2、Z方向地应力X3、围岩或煤层强度X4、围岩或煤体抵抗线长X5、工作面顶板压力X6、围岩变形量X7、围岩体积变形能X8、围岩形状变形能X9以及顶板弯曲能X10,与瓦斯突出强相关的煤层瓦斯含量X11、煤层瓦斯压力X12,当前状态(允许有小波动)的持续时间X13,煤层瓦斯突出倾向性X14,地应力增量X15(ΔX1)、地应力增量X16(ΔX2)、地应力增量X17(ΔX3)、围岩或煤体强度X18(ΔX4)、围岩或煤体抵抗线长X19(ΔX5)、工作面顶板压力X20(ΔX6)、围岩变形量X21(ΔX7)、煤层瓦斯含量X22(ΔX11)、煤层瓦斯压力X23(ΔX12)、围岩和煤层冲击倾向性X24、区域自由空间绝对气压X25、围岩或煤层弹性模量X26等26个输入变量,即M=26。选取冲击地压发生的可能性Y1和煤与瓦斯突出发生的可能性Y2等2个为输出变量,即J=2。设置一个单隐含层,隐含层的节点个数L=84,建立矿山动力灾害预警的神经网络模型如图2所示。
步骤二,向学习样本库中添加样本,当学习样本库中添加的样本达到预设值时,运用预设的学习算法对学习样本库中的样本进行学习,以确定人工神经网络模型中的权值和阈值,得到一个矿山动力灾害的全息模态化预警器。
这里可以分别采用理论计算法、专家评估法和案例追溯法不断添加学习样本库中的样本。
所述理论计算法就是利用一组理论计算公式,对任意给定的一组输入值(X1,X2,…,X26),计算出一组(Y1,Y2),添加该组数据到样本库;所述专家评估法就是对任意给定的一组输入值(X1,X2,…,X26),通过专家评价给出一组(Y1,Y2),添加该组数据到样本库;所述案例追溯法就是每当获取到矿井动力灾害事故或异常征兆,就追溯该时段及之前的多组输入参数(X1,X2,…,X26)的值及相应的显现程度(Y1,Y2),并将其添加到样本库。一旦样本库中有较多的样本(比如超过200条)就可以利用学习算法进行学习,以确定人工神经网络模型中的权值和阈值。这个过程可以新样本的加入或无用样本的淘汰循环进行,不断优化神经网络模型中的权值和阈值,不断提高预警的准确率。我们称确定了权值和阈值的矿山动力灾害预警人工神经网络模型为预警器。
由于用专家评估法和案例追溯法在短时间内不易获取足够多的学习样本,所以本发明首先给出了一个快捷计算方法,对任意给定的一组输入值(X1,X2,…,X26),可以快速计算出显现程度(Y1,Y2),计算方法如下:
P1=min{1,py}
P2=min{1,pc}
P5=x24
P:=P1
P:=P+P2-P×P2
P:=P+P3-P×P3
P:=P+P5-P×P5
P:=P+P6-P×P6
Y1=P
P5=x14
P:=P2
P:=P+P3-P×P3
P:=P+P6-P×P6
P4:=P4+P5-P4×P5
P:=0.2×P+0.9×P4-0.1×P×P4
Y2=P
通过随机产生500组(X1,X2,…,X26)的值,并利用上述计算方法,计算出对应的500组(Y1,Y2),得到500个学习样本,如表1所示。
表1用快捷计算法产生的学习样本
经过4200次学习,其结果列于表2。
表2神经网络学习结果
由表2可以看出学习精度是很高的,说明当取隐含层的节点个数L=84时,该神经网络结构可很好的用于矿山动力灾害的全息学习和预警,由此,可建立一个基于神经网络的矿山动力灾害预警器,简称神经网络预警器。
步骤三,开发一个包括空间检索、插补、降噪、反飘移和信号重建等功能的全息数据挖掘转换器,对于矿山任意一个预警区域,在预警时,首先利用该数据挖掘转换器把矿山实验测试和安全生产中能够获得的相关信息(这些信息包括矿山地理信息、地质信息、生产信息和相关属性信息,例如地形、岩层、断层、陷落柱、褶皱、岩浆侵入、火烧等构造信息,地温、水压、岩移、沉陷、煤柱、井巷、硐室、工作面、采空区等空间分布信息,支护、卸压、充填、瓦斯抽放等灾害防治工艺等参数,以及冲击倾向性测试和相关防突测试数据和相关的微震信号、地音信号、爆破信号、电磁辐射信号、应力信号、矿压信号、锚杆阻力、离层和位移信号、瓦斯浓度、风速、瓦斯抽放流量、水位、涌水量等监测数据),转化成预警模型的输入参数(X1,X2,…,X26)。再利用预警器给出预警结果(Y1,Y2)。
该步骤的关键是如何把能够获得的全息信息,例如地形、岩层、断层、陷落柱、褶皱、岩浆侵入、火烧等构造信息,地温、水压、岩移、沉陷、煤柱、井巷、硐室、工作面、采空区等空间分布信息,支护、卸压、充填、瓦斯抽放等灾害防治工艺等参数,以及冲击倾向性测试和相关防突测试数据和相关的微震信号、地音信号、爆破信号、电磁辐射信号、应力信号、矿压信号、锚杆阻力、离层和位移信号、瓦斯浓度、风速、瓦斯抽放流量、水位、涌水量等监测数据,通过转化成预警模型的输入参数(X1,X2,…,X26)。下面给出一个这样的转化方案:
(1)利用岩性测试和煤体测试得到的单轴抗压强度Rc、煤层的冲击能量指数KE、煤层的单性能量指数WET,煤层的动态破坏时间DT,用下列方法计算出围岩和煤层冲击倾向性X24;
P1=PDT+PRC-PDT×PRC
P2=PWET+PKE-PWET×PKE
X24=0.7×P1+0.5×P2-0.2×P1×P2
(2)利用煤层的瓦斯试验数据,例如瓦斯放散初速度Δp、瓦斯涌出初速度q1、解吸指标Δh2、煤层的硬度系数f、综合指标K=Δp/f、k1钻屑瓦斯解吸指标、每个钻孔沿孔长的最大钻屑量Smax、每个钻孔的最大钻孔涌出初速度qmax、钻孔平均综合指标R=(Smax-1.8)×(qmax-4),用下列方法计算出煤层瓦斯突出倾向性X14;
p1=pΔp+pq-pΔp×pq
p2=pΔh2+pk1-pΔh2×pk1
p3=pK+pR-pK×pR
p4=p1+p2-p1×p2
X14=0.7×p3+0.5×p4-0.2×p3×p4
(3)利用微震定位结果,通过反演岩层或煤层的地震波波速vp和横波波速vs,利用式即可解出岩层或煤层的当前状态对应的围岩或煤层弹性模量E(即X26)和泊松比μ;
其中,ρ为岩层或煤层的密度,是已知量;
(4)利用岩层或煤层的当前状态对应的弹性模量E和原始弹性模量E0,即可求出损伤变量D;
(5)利用损伤变量D和原始抗压强度S0,即可求当前的围岩或煤层强度S(即X4);
S=S0(1-D)
(6)利用三轴应力监测系统直接监测测点的绝对地应力(σx,σy,σz),即(X1,X2,X3),也可以先根据三维空间模型,利用解析法计算出初步的应力分布,在用包括但不限于应力监测、矿压监测的实时数据进行修正,得到(X1,X2,X3);
(7)工作面顶板压力Q(即X6)和围岩变形量U(X7),可以用支架阻力传感器和位移传感器直接测量;
(8)根据地应力(X1,X2,X3)和井巷空间模型,就可计算出围岩或煤体抵抗线长W(即X5);
(9)可用下列公式计算围岩体积变形能UV(X8)、围岩形状变形能Uf(X9)和工作面顶板弯曲能Uw(X10);
其中的l,h和b1分别为工作面顶板的长度,厚度和宽度,q2顶板受均匀载荷,E为弹性模量;
G为剪切模量,可由下式计;
(10)利用下列公式计算煤层瓦斯含量Vt(即X11);
Vt=Qf+Qc+Vc
其中Qf为流量涌出量,可用巷道断面、风速监测值和瓦斯浓度监测值计算;
Qc为瓦斯抽放流量,可用瓦斯抽放系统的纯瓦斯流量监测值得到;
Vc为煤体瓦斯残存含量,可用地面试验获得;
(11)利用下列公式计算煤层瓦斯压力P(即X12);
其中:η为煤层孔隙度、a和b2为瓦斯吸附常数,可通过试验获得;
p0为标准大气压;
(12)区域自由空间中绝对气压(即X25)可直接监测;
(13)地应力增量X15(ΔX1)、地应力增量X16(ΔX2)、地应力增量X17(ΔX3)、围岩或煤体强度增量X18(ΔX4)、围岩或煤体抵抗线长增量X19(ΔX5)、工作面顶板压力增量X20(ΔX6)、围岩变形量增量X21(ΔX7)、煤层瓦斯含量增量X22(ΔX11)、煤层瓦斯压力增量X23(ΔX12)可由当前时段的状态量减去前一时段的相应状态量得到。
不妨假设监测时段的状态不变,即X15=X16=…=X23=0
(14)状态持续时间(即X13)可根据状态变化情况直接计算。不妨假设X13=10。
至此,就完成了矿山全息参数到预警神经网络输入参数的转换,
X=(x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x14,x15,x16,x17,x18,x19,x20,x21,x22,x23,x24,x25,x26)
=(8,7.2,6.5,10.345,8,6,0.2,7.848,0.0001405,2.5,12.5,0.761,10,0.593,0,0,0,0,0,0,0,0,0,0.933,0.1024,5)。
把X输入到第一不建立的神经网络预警器得到:(冲击地压发生的可能性,煤与瓦斯突出发生的可能性)=(Y1,Y2)=(0.9735,0.7788)。如果把预警级别定义为:
{红色,橙色,黄色,蓝色,绿色}={[0.9,1],[0.8,0.9),[0.7,0.8),[0.4,0.6),[0,0.4)}
对该例来说冲击地压达到红色预警级别,煤与瓦斯突出达到黄色预警级别,系统可以按这些级别发布相应预警信息。
值得说明的是:除了本方案给出的计算公式外,也可以通过其它更科学的方法实现矿山全息参数到输入参数(X1,X2,…,X26)的转换。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (4)
1.一种矿山动力灾害的全息预警方法,其特征在于,方法包括:
步骤一,以直接影响矿山动力灾害发生与否的模态化参数集为输入变量集,以冲击地压发生的可能性和煤与瓦斯突出发生的可能性作为输出变量,建立人工神经网络模型;
步骤二,向学习样本库中添加样本,当学习样本库中添加的样本达到预设值时,运用预设的学习算法对学习样本库中的样本进行学习,以确定人工神经网络模型中的权值和阈值,得到一个矿山动力灾害的全息模态化预警器;该过程包括新样本的加入或无用样本的淘汰,不断优化神经网络模型中的权值和阈值;
步骤三,根据全息模态化预警器设置一包括空间检索、插补、降噪、反飘移和信号重建功能的全息数据挖掘转换器;
对矿山任意一个区域进行预警时,利用该全息数据挖掘转换器把矿山实验测试和安全生产中获得的相关信息转化成预警模型的输入参数,再利用全息模态化预警器给出预警结果。
2.根据权利要求1所述的矿山动力灾害的全息预警方法,其特征在于,
矿山实验测试和安全生产中获得的相关信息包括:矿山地理信息、地质信息、生产信息和相关属性信息;
相关属性信息包括:地形、岩层、断层、陷落柱、褶皱、岩浆侵入、火烧的构造信息,地温、水压、岩移、沉陷、煤柱、井巷、硐室、工作面、采空区空间分布信息,支护、卸压、充填、瓦斯抽放灾害防治工艺的参数,以及冲击倾向性测试数据,防突测试数据,微震信号、地音信号、爆破信号、电磁辐射信号、应力信号、矿压信号、锚杆阻力、离层和位移信号、瓦斯浓度、风速、瓦斯抽放流量、水位、涌水量监测数据。
3.根据权利要求1或2所述的矿山动力灾害的全息预警方法,其特征在于,
步骤三中,利用该全息数据挖掘转换器把矿山实验测试和安全生产中获得的相关信息转化成预警模型的输入参数还包括:
(1)利用岩性测试和煤体测试得到的单轴抗压强度Rc、煤层的冲击能量指数KE、煤层的单性能量指数WET,煤层的动态破坏时间DT,用下列方法计算出围岩和煤层冲击倾向性X24;
P1=PDT+PRC-PDT×PRC
P2=PWET+PKE-PWET×PKE
X24=0.7×P1+0.5×P2-0.2×P1×P2
(2)利用煤层的瓦斯试验数据,瓦斯放散初速度Δp、瓦斯涌出初速度q1、解吸指标Δh2、煤层的硬度系数f、综合指标K=Δp/f、k1钻屑瓦斯解吸指标、每个钻孔沿孔长的最大钻屑量Smax、每个钻孔的最大钻孔涌出初速度qmax、钻孔平均综合指标R=(Smax-1.8)×(qmax-4),用下列方法计算出煤层瓦斯突出倾向性X14;
p1=pΔp+pq-pΔp×pq
p2=pΔh2+pk1-pΔh2×pk1
p3=pK+pR-pK×pR
p4=p1+p2-p1×p2
X14=0.7×p3+0.5×p4-0.2×p3×p4
(3)利用微震定位结果,通过反演岩层或煤层的地震波波速vp和横波波速vs,利用式即可解出岩层或煤层的当前状态对应的围岩或煤层弹性模量E(X26)和泊松比μ;
其中,ρ为岩层或煤层的密度,是已知量;
(4)利用岩层或煤层的当前状态对应的弹性模量E和原始弹性模量E0,即可求出损伤变量D;
(5)利用损伤变量D和原始抗压强度S0,即可求当前的围岩或煤层强度S(X4);
S=S0(1-D)
(6)利用三轴应力监测系统直接监测测点的绝对地应力(σx,σy,σz),即(X1,X2,X3),也可以先根据三维空间模型,利用解析法计算出初步的应力分布,再 用包括但不限于应力监测、矿压监测的实时数据进行修正,得到(X1,X2,X3);
(7)工作面顶板压力Q(X6)和围岩变形量U(X7),可以用支架阻力传感器和位移传感器直接测量;
(8)根据地应力(X1,X2,X3)和井巷空间模型,就可计算出围岩或煤体抵抗线长W(X5);
(9)可用下列公式计算围岩体积变形能UV(X8)、围岩形状变形能Uf(X9)和工作面顶板弯曲能Uw(X10);
其中的l,h和b1分别为工作面顶板的长度,厚度和宽度,q2顶板受均匀载荷,E为弹性模量;G为剪切模量,可由下式计;
(10)利用下列公式计算煤层瓦斯含量Vt(X11);
Vt=Qf+Qc+Vc
其中Qf为流量涌出量,可用巷道断面、风速监测值和瓦斯浓度监测值计算;
Qc为瓦斯抽放流量,可用瓦斯抽放系统的纯瓦斯流量监测值得到;
Vc为煤体瓦斯残存含量,可用地面试验获得;
(11)利用下列公式计算煤层瓦斯压力P(X12);
其中:η为煤层孔隙度、a和b2为瓦斯吸附常数,可通过试验获得;
p0为标准大气压;
(12)区域自由空间中绝对气压(X25)可直接监测;
(13)对应X1的地应力增量X15、对应X2的地应力增量X16、对应X3的地应力增量X17、围岩或煤体强度增量X18、围岩或煤体抵抗线长增量X19、工作面顶板压力增量X20、围岩变形量增量X21、煤层瓦斯含量增量X22、煤层瓦斯压力增量X23可由当前时段的状态量减去前一时段的相应状态量得到。
4.根据权利要求1或2所述的矿山动力灾害的全息预警方法,其特征在于,
步骤一还包括:
选取预警区域和地点的动力灾害公用参数X方向地应力X1,Y方向地应力X2,Z方向地应力X3,围岩或煤层强度X4,围岩或煤体抵抗线长X5,工作面顶板压力X6,围岩变形量X7,围岩体积变形能X8,围岩形状变形能X9以及工作面顶板弯曲能X10,与瓦斯突出强相关的煤层瓦斯含量X11、煤层瓦斯压力X12,当前状态的持续时间X13,煤层瓦斯突出倾向性X14,对应X1的地应力增量X15、对应X2的地应力增量X16、对应X3的地应力增量X17、围岩或煤体强度增量X18、围岩或煤体抵抗线长增量X19、工作面顶板压力增量X20、围岩变形量增量X21、煤层瓦斯含量增量X22、煤层瓦斯压力增量X23、围岩和煤层冲击倾向性X24、区域自由空间绝对气压X25、围岩或煤层弹性模量X26,为26个输入变量;选取冲击地压发生的可能性Y1和煤与瓦斯突出发生的可能性Y2,为2个为输出变量;设置一个单隐含层,隐含层的节点个数L,建立矿山动力灾害预警的神经网络模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710851433.1A CN107563092B (zh) | 2017-09-19 | 2017-09-19 | 一种矿山动力灾害的全息预警方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710851433.1A CN107563092B (zh) | 2017-09-19 | 2017-09-19 | 一种矿山动力灾害的全息预警方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107563092A CN107563092A (zh) | 2018-01-09 |
CN107563092B true CN107563092B (zh) | 2020-08-04 |
Family
ID=60981811
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710851433.1A Active CN107563092B (zh) | 2017-09-19 | 2017-09-19 | 一种矿山动力灾害的全息预警方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107563092B (zh) |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108506041B (zh) * | 2018-01-31 | 2019-07-19 | 山东蓝光软件有限公司 | 一种基于实时监测数据的动力灾害模态预警方法 |
CN108470222B (zh) * | 2018-02-08 | 2022-11-11 | 山东蓝光软件有限公司 | 一种用于典型事件预测预报的多征兆增信方法 |
CN108919358B (zh) * | 2018-04-23 | 2020-07-24 | 中国矿业大学 | 一种矿震灾害判别及信号重构方法 |
CN108876009A (zh) * | 2018-05-14 | 2018-11-23 | 中国神华能源股份有限公司 | 煤矿事故预测模型确定及监测方法、存储介质和电子设备 |
CN109471378B (zh) * | 2018-12-07 | 2021-06-25 | 重庆大学 | 突出煤层多物理场参数可视化采集方法及系统 |
CN109854303A (zh) * | 2019-02-12 | 2019-06-07 | 青岛理工大学 | 地下开采矿山岩爆监测预警方法 |
CN109632016A (zh) * | 2019-02-20 | 2019-04-16 | 湖北理工学院 | 岩土体硐室开挖与围岩应力、应变监测试验装置及其方法 |
CN110130985B (zh) * | 2019-05-01 | 2020-07-17 | 山东科技大学 | 一种煤矿岩爆危险性预报方法 |
CN110159347B (zh) * | 2019-05-05 | 2020-05-08 | 北京科技大学 | 一种深部高应力坚硬顶板采场动力灾害监测预警方法 |
CN110533887B (zh) * | 2019-08-07 | 2022-08-26 | 山东蓝光软件有限公司 | 一种基于实时监测数据的煤与瓦斯突出灾害离散模态预警方法,装置及存储介质 |
CN111027687B (zh) * | 2019-12-29 | 2023-01-10 | 鄂尔多斯市伊化矿业资源有限责任公司 | 煤矿工作面回采期间矿压预测治理方法及系统 |
CN111177841B (zh) * | 2020-01-02 | 2021-09-14 | 同济大学 | 一种数物空间实时融合的隧道结构火灾试验系统及方法 |
CN111472840B (zh) * | 2020-04-07 | 2020-11-13 | 北京科技大学 | 一种采动围岩地压灾害智能预测方法及系统 |
CN112031872B (zh) * | 2020-07-15 | 2022-12-27 | 中煤科工开采研究院有限公司 | 一种冲击地压全息预警方法及装置、存储介质、计算设备 |
CN114647221B (zh) * | 2022-05-23 | 2022-08-02 | 承德京联电子科技有限公司 | 一种基于物联网的矿山生产链智能管控系统 |
CN117607973B (zh) * | 2024-01-23 | 2024-04-12 | 山东科技大学 | 一种可判别陷落柱时空动态活化程度的方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070230902A1 (en) * | 2006-03-31 | 2007-10-04 | Masstech Group Inc. | Dynamic disaster recovery |
CN103912310B (zh) * | 2014-04-15 | 2016-05-25 | 辽宁工程技术大学 | 一种矿井动力灾害预测方法 |
CN105673075A (zh) * | 2016-01-13 | 2016-06-15 | 中国矿业大学(北京) | 一种煤岩动力灾害多参量无线监测综合预警技术与方法 |
CN106437853B (zh) * | 2016-09-27 | 2019-02-15 | 西安科技大学 | 一种煤岩冲击地压动力灾害预警方法 |
-
2017
- 2017-09-19 CN CN201710851433.1A patent/CN107563092B/zh active Active
Non-Patent Citations (2)
Title |
---|
基于BP神经网络的煤与瓦斯突出危险性预测的研究;李春辉;《工程科技Ⅰ辑》;20120215;43-47 * |
煤与瓦斯突出、冲击地压复合动力灾害一体化研究;潘一山;《煤炭学报》;20160130;105、107 * |
Also Published As
Publication number | Publication date |
---|---|
CN107563092A (zh) | 2018-01-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107563092B (zh) | 一种矿山动力灾害的全息预警方法 | |
Li et al. | Microseismic monitoring and numerical simulation on the stability of high-steep rock slopes in hydropower engineering | |
Ding et al. | Stability of large underground caverns excavated in layered rock masses with steep dip angles: a case study | |
CN103559383B (zh) | 核电站防波堤抗强震性能预测与评价方法 | |
Havaej et al. | The possible role of brittle rock fracture in the 1963 Vajont Slide, Italy | |
CN108842821B (zh) | 一种钻爆法修建海底隧道合理埋深的计算方法 | |
CN110489926B (zh) | 一种基于地层变形的海底隧道施工过程控制方法及系统 | |
Dong et al. | Potential failure patterns of a large landslide complex in the Three Gorges Reservoir area | |
CN105182412A (zh) | 煤矿地下水库煤柱坝体地震安全性的检测方法 | |
CN105701712A (zh) | 判断煤与瓦斯是否会突出的方法 | |
CN112016201B (zh) | 基于dfos应变重构深部采场超前支承压力演化模型的方法 | |
Wang et al. | A novel method for analyzing the factors influencing ground settlement during shield tunnel construction in upper-soft and lower-hard fissured rock strata considering the coupled hydromechanical properties | |
CN113591338A (zh) | 一种荷载和地下水开采引发地面沉降三维变参数全耦合模拟计算方法 | |
Jia et al. | Stability analysis of shallow goaf based on field monitoring and numerical simulation: A case study at an Open-Pit Iron Mine, China | |
Sun et al. | Visualization and early warning analysis of damage degree of surrounding rock mass in underground powerhouse | |
CN116006264A (zh) | 基于声波探测的矿井冲击地压与水害协同预警方法 | |
Zheng et al. | Evolutionary mechanisms of deep-seated toppling deformation of rock slope under different free-face conditions | |
Wu et al. | Assessment of water inrush vulnerability from overlying aquifer using GIS–AHP-based ‘three maps–two predictions’ method: a case study in Hulusu coal mine, China | |
Yang et al. | Evolution of Mining-Induced Stress in Downward Mining of Short-Distance Multiseam | |
Zhang et al. | Research on the underground water inrush mechanism based on the influence of fault | |
Zhang et al. | Influence of Underground Mining Direction Based on Particle Flow on Deformation and Failure of Loess Gully Area | |
CN116910886B (zh) | 一种岩溶地区的冲击能量预警方法及系统 | |
Cabangon et al. | Seismic vulnerability of shallow tunnels subjected to far-field long-period ground motions | |
He et al. | Study on Seismic Dynamic Response of Deep-Buried Tunnel Based on Tectonic Stress | |
CN116842789B (zh) | 一种考虑生产时间窗的加密井井壁稳定分析方法 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |