CN105891874B - 一种采动煤岩体突水微震监测方法 - Google Patents
一种采动煤岩体突水微震监测方法 Download PDFInfo
- Publication number
- CN105891874B CN105891874B CN201610513797.4A CN201610513797A CN105891874B CN 105891874 B CN105891874 B CN 105891874B CN 201610513797 A CN201610513797 A CN 201610513797A CN 105891874 B CN105891874 B CN 105891874B
- Authority
- CN
- China
- Prior art keywords
- rock
- coal
- microseism
- gushing water
- trigger event
- 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
- 239000011435 rock Substances 0.000 title claims abstract description 181
- 239000003245 coal Substances 0.000 title claims abstract description 152
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 119
- 238000012544 monitoring process Methods 0.000 title claims abstract description 98
- 238000000034 method Methods 0.000 title claims abstract description 60
- 238000004458 analytical method Methods 0.000 claims abstract description 28
- 238000009826 distribution Methods 0.000 claims abstract description 23
- 230000009172 bursting Effects 0.000 claims abstract description 21
- 230000002776 aggregation Effects 0.000 claims abstract description 6
- 238000004220 aggregation Methods 0.000 claims abstract description 6
- 238000009434 installation Methods 0.000 claims description 10
- 230000035939 shock Effects 0.000 claims description 10
- 230000004913 activation Effects 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 9
- 230000005284 excitation Effects 0.000 claims description 8
- 230000001133 acceleration Effects 0.000 claims description 7
- 230000005855 radiation Effects 0.000 claims description 7
- 238000012360 testing method Methods 0.000 claims description 7
- 238000002474 experimental method Methods 0.000 claims description 6
- 238000004880 explosion Methods 0.000 claims description 6
- 238000004519 manufacturing process Methods 0.000 claims description 6
- 238000005457 optimization Methods 0.000 claims description 6
- 230000003245 working effect Effects 0.000 claims description 6
- 238000001514 detection method Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims description 3
- 239000002245 particle Substances 0.000 claims description 3
- 230000021715 photosynthesis, light harvesting Effects 0.000 claims description 3
- 238000004321 preservation Methods 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 230000004907 flux Effects 0.000 claims 1
- 238000009472 formulation Methods 0.000 claims 1
- 230000008569 process Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000005065 mining Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000035945 sensitivity Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 230000006641 stabilisation Effects 0.000 description 2
- 238000011105 stabilization Methods 0.000 description 2
- 238000009412 basement excavation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Business, Economics & Management (AREA)
- Emergency Management (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明具体涉及一种采动煤岩体突水微震监测方法，包括：根据采动煤岩体突水微震监测系统的工程要求及实施条件，确定采动煤岩体突水微震监测的系统架构；对检波器阵列拾取的微震信号，采用阈值算法和STA/LTA算法判断是否存在采动煤岩体突水微震触发事件，若是，保存触发事件波形；采用STA/LTA算法拾取p波到时t，并通过单纯形法迭代求解出触发事件的发生位置；计算触发事件发生时的微震辐射能量；分析采动煤岩体微震辐射能量的时空分布特征和聚集程度，设定微震辐射能量释放阈值，以识别出采动煤岩体潜在突水灾害区域。本发明的技术方案，能实现对采动煤岩体突水的全局、实时和动态监测，监测方法简单，监测效率高，便于推广实施。
Description
技术领域
本发明涉及矿井水文地质技术领域，具体涉及一种采动煤岩体突水微震监测方法。
背景技术
在煤矿开采过程中，煤岩体突水、淹井事故时有发生，严重威胁煤矿的安全生产。特别是近年来，矿井开采深度已经达千米，采动煤岩体突水问题就显得更加突出。煤矿开采过程中不同的地质赋存环境、岩石性质、地应力分布和断裂构造以及破碎带的发育情况等都有可能诱发不同机理的矿井突水。因此，采动煤岩体突水是由多种因素共同作用造成的，对于突水灾害的监测需要综合考虑多种因素的影响，圈定开采过程中可能的突水通道和突水危险区域。
目前，煤矿突水的监测方法主要有钻孔注水法、岩层移动钻孔探测法、水文地质钻孔探测法、三维地球物理探测技术等。近几年，以地质雷达技术、电阻率法、瞬变电磁法和红外探测方法为代表的地球物理探测技术方法得到了飞速发展，这些监测技术均能较好地获得煤层底板潜在的突水点和突水通道，但这些监测结果往往只获得采动煤岩体突水全过程中的某一静态时刻的局部结果，未能实现完整的突水全过程中突水通道形成的连续动态监测，并且受监测环境中电磁信号干扰的影响较大，监测结果可靠性不高，从而制约了采动煤岩体突水的监测和预警技术的发展。
近年来，高精度微震监测技术的发展为岩体微破裂的全局、实时和动态监测提供了一种新方法。微震监测技术通过对岩体在变形破坏过程释放的微地震波形的拾取、反演和分析，获取岩体微破裂的发生的时间、位置和强度等震源信息，进而对岩体破坏的活动范围、稳定性及其发展趋势做出定性、定量评价。采用微震监测技术对采动煤岩体突水进行监测，突破了以流量、水温及水压等表象信息监测为依据进行采动煤岩体突水灾害预警的传统思路，引入微震能量与煤岩体劣化表征关系，提出从采动诱发的煤岩体微破裂导致煤岩体力学性质劣化并形成突水通道的本质出发，以微震监测技术和水文地质资料分析为手段，进行采动煤岩体突水的监测和预警。
发明内容
有鉴于此，本发明提供一种采动煤岩体突水微震监测方法，实现实时监测采动煤岩体可能出现的突水灾害。
为实现以上目的，本发明采用如下技术方案：
一种采动煤岩体突水微震监测方法，包括：
步骤S1、根据采动煤岩体突水微震监测系统的工程要求及实施条件，确定采动煤岩体突水微震监测的系统架构；其中，所述系统架构包括检波器阵列在煤矿回采工作面的回采巷道中的安装方式及安装位置，所述检波器阵列用于拾取被监测采动煤岩体产生的微震信号；
步骤S2、对所述检波器阵列拾取的微震信号，采用阈值算法和STA/LTA算法判断是否存在采动煤岩体突水微震触发事件，若是，保存触发事件波形；
步骤S3、滤除所述触发事件波形中的噪声信号，对滤噪后的触发事件波形，采用STA/LTA算法拾取p波到时t，并根据采动煤岩体整体最优波速vopt通过单纯形法迭代求解出所述触发事件的发生位置；
步骤S4、根据公式(1)计算所述触发事件发生时的微震辐射能量；其中，ρ为煤岩体密度；v为煤岩体弹性波速；R为震源到检波器阵列的距离；Jc为辐射能通量，通过质点速度谱在频域中积分获得；Fc为地震波辐射类型经验系数；
步骤S5、分析所述触发事件发生时微震辐射能量的时空分布特征和聚集程度，设定微震辐射能量释放阈值，以识别出采动煤岩体潜在突水灾害区域。
优选地，所述采动煤岩体突水微震监测方法，还包括：
步骤S6、根据所述触发事件发生时微震辐射能量的时空分布特征及采动煤岩体潜在突水灾害区域，确定采动煤岩体潜在突水通道赋存的主控地质因素，建立煤岩体突水灾害分析的三维地质力学模型，为采动煤岩体突水灾害的监测预警提供模型支持；
步骤S7、根据所述触发事件的时空分布规律和所述触发事件发生时微震辐射能量的时空分布特征，利用能量耗散原理，将煤矿围岩微震辐射能量转化为围岩损伤弱化的力学参数，为采动煤岩体的数值分析提供力学修正参数。
优选地，所述步骤S1具体包括：
步骤S10、根据采动煤岩体突水微震监测系统的工程要求及实施条件，确定采动煤岩体突水微震监测系统的系统架构；
步骤S11、通过现场波速测试试验获取可供选择的岩体弹性波波速范围(vmin,vmax)，其中，vmin和vmax分别表示岩体弹性波波速的最小值和最大值；
步骤S12、对采动煤岩体现场进行定点爆破或敲击试验，人工激发弹性波，滤除所述人工激发弹性波中的噪声信号，并对滤除噪声信号后的人工激发弹性波，采用STA/LTA算法拾取p波到时t；
步骤S13、遍历岩体弹性波波速vi∈(vmin,vmax)内的所有弹性波波速，采用单纯形法迭代求解微震震源位置Xi(xi,yi,zi,vi)，并利用公式minE(vi)＝||Xi-X0||，i＝1,2L n，n≥2(2)，计算出采用单纯形法迭代求解的震源位置Xi(xi,yi,zi,vi)与人工震源位置X0(x0,y0,z0)之间的误差最小值minE(vi)，其对应的弹性波波速为采动煤岩体整体最优波速vopt；
步骤S14、若minE(vi)的值不在预设误差范围内，返回步骤S10，调整采动煤岩体突水微震监测系统的系统架构，直至minE(vi)的值在预设误差范围；其中，调整采动煤岩体突水微震监测系统的系统架构包括调整采动煤岩体突水微震监测系统中检波器阵列的安装方式及安装位置。
需要说明的是，所述预设误差范围根据用户实际需要和微震监测工程要求进行设置。
优选地，所述步骤S2具体包括：
步骤S21、采集所述检波器阵列输出的微震信号，判断所述微震信号初至时的幅值As是否大于第一触发阈值At；其中，At根据采动煤岩体突水微震监测系统工程现场的背景噪声值确定；
步骤S22、若As大于At，判断所述微震信号是否满足值，若是，则判定所述微震信号存在采动煤岩体突水微震触发事件；其中，STA为微震信号短时能量平均值，LTA为微震信号长时能量平均值，T为第二触发阈值。
优选地，所述步骤S21中At根据如下步骤进行确定：
步骤S210、计算预设时段内采动煤岩体突水微震监测系统工程现场的背景噪声的平均振幅值；
步骤S211、若微震信号初至时的幅值As大于所述平均振幅值，则判定所述微震信号中可能存在采动煤岩体突水微震触发事件；
步骤S212、若被监测采动煤岩体实际状态稳定，利用步骤S211判定的可能存在的触发事件每分钟低于30个，且微震信号的信噪比在预设比值范围内，则将所述平均振幅值设为第一触发阈值At；否则，返回步骤S210，调整预设时段的长度重新计算所述平均振幅。
需要说明的是，所述预设比值根据用户实际需要和微震监测工程要求进行设置。
优选地，所述步骤S3中滤除所述触发事件波形中的噪声信号具体为：
对工程现场大于预设幅值的噪声信号进行快速傅里叶变换，以获取噪声信号的时域和频域特征，建立噪声信号特征库；根据噪声信号特征库中噪声信号的时域和频域特征滤除所述触发事件波形中的噪声信号。
需要说明的是，所述预设幅值根据工程经验和现场监测试验确定。
优选地，所述步骤S7具体包括：
根据公式(3)计算微震触发事件发生后岩体的损伤变量D；根据损伤变量D，对采动煤岩体的力学参数进行修正。
其中，η地震效率，E表示震源尺寸范围内的岩体单元损伤后被微震监测系统拾取的辐射能，Ue为岩体单元可释放总应变能，Ue由公式
(4)计算，其中E0为岩体初始的弹性模量，σ1、σ2和σ3分别为岩体第一、第二和第三主应力，ν为岩体泊松比。
优选地，所述步骤S1中采动煤岩体突水微震监测的系统架构包括硬件架构和软件架构，其硬件架构包括：
检波器阵列，包括预设数量的检波器，所述检波器阵列设置在回采工作面的回采巷道中；
信号采集单元，包括：集线器、模数转换器和存储单元，其中，所述模数转换器的输入端通过所述集线器与所述检波器阵列连接，输出端与所述存储单元连接；
信号分析单元，包括：处理子单元和分析子单元，其中，所述处理子单元的输入端与所述存储单元连接，输出端与所述分析子单元连接；所述处理子单元用于根据存储单元中的微震信号，计算震源参数，所述震源参数包括微震发生时刻、微震辐射能量和震源位置；所述分析子单元用于根据处理子单元输出的震源参数，分析采动煤岩体微震事件和微震能量时空分布规律；
预警决策单元，与所述分析子单元连接，用于根据所述采动煤岩体微震事件和微震能量时空分布规律进行突水灾害预警。
优选地，所述检波器为单分量和/或三分量检波器，所述单分量检波器包括单分量加速度检波器和/或单分量速度检波器，所述三分量检波器包括三分量加速度检波器和/或三分量速度检波器。
优选地，所述信号采集单元为多个，当所述检波器为单分量检波器时，每六个检波器为一组，连接一个信号采集单元；当所述检波器为三分量检波器时，每两个检波器为一组，连接一个信号采集单元。
本发明采用以上技术方案，至少具备以下有益效果：
本发明提供的这种采动煤岩体突水微震监测方法，通过对检波器阵列输出的微震信号，采用阈值算法和STA/LTA算法判断是否存在采动煤岩体突水微震触发事件，并采用STA/LTA算法拾取p波到时t，并通过单纯形法迭代求解出触发事件的发生位置，根据触发事件的辐射能量，分析触发事件发生时微震辐射能量的时空分布特征和聚集程度，通过设定微震辐射能量释放阈值，从而识别出采动煤岩体潜在突水灾害区域，相比现有技术，能实现对采动煤岩体突水的全局、实时和动态监测，监测方法简单，监测效率高，便于推广实施。
附图说明
图1为本发明一实施例提供的一种采动煤岩体突水微震监测方法的流程示意图；
图2为本发明另一实施例提供的一种采动煤岩体突水微震监测方法的流程示意图；
图3为本发明一实施例提供的采动煤岩体突水微震监测系统的硬件架构结构示意图。
具体实施方式
下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。
参见图1，本发明一实施例提供的一种采动煤岩体突水微震监测方法，包括：
步骤S1、根据采动煤岩体突水微震监测系统的工程要求及实施条件，确定采动煤岩体突水微震监测的系统架构；其中，所述系统架构包括检波器阵列在煤矿回采工作面的回采巷道中的安装方式及安装位置，所述检波器阵列用于拾取被监测采动煤岩体产生的微震信号；
步骤S2、对所述检波器阵列拾取的微震信号，采用阈值算法和STA/LTA算法判断是否存在采动煤岩体突水微震触发事件，若是，保存触发事件波形；
步骤S3、滤除所述触发事件波形中的噪声信号，对滤噪后的触发事件波形，采用STA/LTA算法拾取p波到时t，并根据采动煤岩体整体最优波速vopt通过单纯形法迭代求解出所述触发事件的发生位置；
步骤S4、根据公式(1)计算所述触发事件发生时的微震辐射能量；其中，ρ为煤岩体密度；v为煤岩体弹性波速；R为震源到检波器阵列的距离；Jc为辐射能通量，通过质点速度谱在频域中积分获得；Fc为地震波辐射类型经验系数；
步骤S5、分析所述触发事件发生时微震辐射能量的时空分布特征和聚集程度，设定微震辐射能量释放阈值，以识别出采动煤岩体潜在突水灾害区域。
需要说明的是，所述步骤S5中，微震辐射能量释放阈值根据工程经验值进行设置。设定微震辐射能量释放阈值后，若某个被监测区域的采动煤岩体的微震辐射能量大于微震辐射能量释放阈值，则判定该被监测区域为潜在突水灾害区域。
由上述技术方案可知，本发明提供的这种采动煤岩体突水微震监测方法，通过对检波器阵列输出的微震信号，采用阈值算法和STA/LTA算法判断是否存在采动煤岩体突水微震触发事件，并采用STA/LTA算法拾取p波到时t，并通过单纯形法迭代求解出触发事件的发生位置，根据触发事件的辐射能量，分析触发事件发生时微震辐射能量的时空分布特征和聚集程度，通过设定微震辐射能量释放阈值，从而识别出采动煤岩体潜在突水灾害区域，相比现有技术，能实现对采动煤岩体突水的全局、实时和动态监测，监测方法简单，监测效率高，便于推广实施。
参见图2，优选地，所述采动煤岩体突水微震监测方法，还包括：
步骤S6、根据所述触发事件发生时微震辐射能量的时空分布特征及采动煤岩体潜在突水灾害区域，确定采动煤岩体潜在突水通道赋存的主控地质因素，建立煤岩体突水灾害分析的三维地质力学模型，为采动煤岩体突水灾害的监测预警提供模型支持；
步骤S7、根据所述触发事件的时空分布规律和所述触发事件发生时微震辐射能量的时空分布特征，利用能量耗散原理，将煤矿围岩微震辐射能量转化为围岩损伤弱化的力学参数，为采动煤岩体的数值分析提供力学修正参数。
优选地，所述步骤S1具体包括附图中未示出的：
步骤S10、根据采动煤岩体突水微震监测系统的工程要求及实施条件，确定采动煤岩体突水微震监测系统的系统架构；
步骤S11、通过现场波速测试试验获取可供选择的岩体弹性波波速范围(vmin,vmax)，其中，vmin和vmax分别表示岩体弹性波波速的最小值和最大值；
步骤S12、对采动煤岩体现场进行定点爆破或敲击试验，人工激发弹性波，滤除所述人工激发弹性波中的噪声信号，并对滤除噪声信号后的人工激发弹性波，采用STA/LTA算法拾取p波到时t；
步骤S13、遍历岩体弹性波波速vi∈(vmin,vmax)内的所有弹性波波速，采用单纯形法迭代求解微震震源位置Xi(xi,yi,zi,vi)，并利用公式minE(vi)＝||Xi-X0||，i＝1,2L n，n≥2(2)，计算出采用单纯形法迭代求解的震源位置Xi(xi,yi,zi,vi)与人工震源位置X0(x0,y0,z0)之间的误差最小值minE(vi)，其对应的弹性波波速为采动煤岩体整体最优波速vopt；
步骤S14、若minE(vi)的值不在预设误差范围内，返回步骤S10，调整采动煤岩体突水微震监测系统的系统架构，直至minE(vi)的值在预设误差范围；其中，调整采动煤岩体突水微震监测系统的系统架构包括调整采动煤岩体突水微震监测系统中检波器阵列的安装方式及安装位置。
需要说明的是，所述预设误差范围根据用户实际需要和微震监测工程要求进行设置。
可以理解的是，检波器的安装方式及安装位置不仅影响微震信号的监测，还影响不同微震定位算法的定位速度、精度及定位结果的唯一性。合理的检波器布置方案不仅能够更大范围地监测到更多有效微震信号，而且能使定位算法快速准确地确定震源位置和发震时间。目前，检波器布置主要是根据经验，不同经验的工作人员布置的检波器监测到的微震信号往往差异较大，不能使检波器最大程度地监测到有效微震信号，也很难保证检波器形成一个良性的阵列，致使微震定位速度及精度受到不同程度的影响。
由于采动煤岩体在变形破坏过程中会诱发微破裂的萌生和发展，常伴随着岩石弹性应变能的释放，并以弹性波的形式在岩石内传播。本发明提供的这种采动煤岩体突水微震监测方法，利用煤岩体的整体波速模型，对检波器阵列的布置位置进行优化搜索，从而使检波器阵列更精准更快速地获取微震信号。
可以理解的是，所述煤岩体整体波速模型可通过测定场区岩体弹性波的波速来获取，并根据现场试验进行调整，减小定位误差。
优选地，所述步骤S2具体包括附图中未示出的：
步骤S21、采集所述检波器阵列拾取的微震信号，判断所述微震信号初至时的幅值As是否大于第一触发阈值At；其中，At根据采动煤岩体突水微震监测系统工程现场的背景噪声值确定；
步骤S22、若As大于At，判断所述微震信号是否满足值，若是，则判定所述微震信号存在采动煤岩体突水微震触发事件；其中，STA为微震信号短时能量平均值，LTA为微震信号长时能量平均值，T为第二触发阈值。
其中，T根据采动煤岩体现场微震触发试验来确定，以保证触发事件判定的准确度。
可以理解的是，对检波器阵列输出的微震信号，采用阈值算法和STA/LTA算法共同判断是否存在采动煤岩体突水微震触发事件，可以提高触发事件判断的准确性。
优选地，所述步骤S21中At根据如下附图中未示出的步骤进行确定：
步骤S210、计算预设时段内采动煤岩体突水微震监测系统工程现场的背景噪声的平均振幅值；
步骤S211、若微震信号初至时的幅值As大于所述平均振幅值，则判定所述微震信号中可能存在采动煤岩体突水微震触发事件；
步骤S212、若被监测采动煤岩体实际状态稳定，利用步骤S211判定的可能存在的触发事件每分钟低于30个，且微震信号的信噪比在预设比值范围内，则将所述平均振幅值设为第一触发阈值At；否则，返回步骤S210，调整预设时段的长度重新计算所述平均振幅。
需要说明的是，所述预设比值范围根据用户实际需要和微震监测工程要求进行设置。
优选地，所述步骤S3中滤除所述触发事件波形中的噪声信号具体为：
对工程现场大于预设幅值的噪声信号进行快速傅里叶变换，以获取噪声信号的时域和频域特征，建立噪声信号特征库；根据噪声信号特征库中噪声信号的时域和频域特征滤除所述触发事件波形中的噪声信号。
需要说明的是，所述预设幅值根据工程经验和现场监测试验确定。
针对煤矿开采而言，检波器所得到的信号具有多样性、突发性和不确定性，微震信号在传输过程中，除了带有有用的信息外，同时还会夹杂如机械噪音、环境噪音等干扰无用信号，因此对触发事件波形进行滤波，有利于提高后续微震监测结果的准确度。
优选地，所述步骤S7具体包括：
根据公式(3)计算微震触发事件发生后岩体的损伤变量D；根据损伤变量D，对采动煤岩体的力学参数进行修正。
其中，η地震效率，E表示震源尺寸范围内的岩体单元损伤后被微震监测系统拾取的辐射能，Ue为岩体单元可释放总应变能，Ue由公式
(4)计算，其中E0为岩体初始的弹性模量，σ1、σ2和σ3分别为岩体第一、第二和第三主应力，ν为岩体泊松比。
参见图3，优选地，所述步骤S1中采动煤岩体突水微震监测的系统架构包括硬件架构和软件架构，其硬件架构包括：
检波器阵列，包括预设数量的检波器2，所述检波器阵列设置在回采工作面的回采巷道1中；
信号采集单元3，包括：集线器31、模数转换器32和存储单元33，其中，所述模数转换器32的输入端通过所述集线器31与所述检波器阵列连接，输出端与所述存储单元33连接；
信号分析单元4，包括：处理子单元41和分析子单元42，其中，所述处理子单元41的输入端与所述存储单元33连接，输出端与所述分析子单元42连接；所述处理子单元用于根据存储单元中的微震信号，计算震源参数，所述震源参数包括微震发生时刻、微震辐射能量和震源位置；所述分析子单元用于根据处理子单元输出的震源参数，分析采动煤岩体微震事件和微震能量时空分布规律；
预警决策单元，与所述分析子单元连接，用于根据所述采动煤岩体微震事件和微震能量时空分布规律进行突水灾害预警。
需要说明的是，所述预设数据根据用户实际需要和微震监测工程要求进行设置。
优选地，所述模数转换器32为24位模数转换器。
优选地，所述检波器2为单分量和/或三分量检波器，所述单分量检波器包括单分量加速度检波器和/或单分量速度检波器，所述三分量检波器包括三分量加速度检波器和/或三分量速度检波器。
优选地，当所述检波器2为单分量/三分量加速度检波器时，所述检波器2的频响范围为50～10000Hz，灵敏度为1V/g或30V/g；当所述检波器2为单分量/三分量速度检波器时，所述检波器2的频响范围为4.5～5000Hz灵敏度为27.6V/m/s或43.3V/m/s。
作为优选的技术方案，选用频响范围为15～1500Hz，灵敏度为43.3V/m/s的速度检波器。检波器能够对回采工作面的微破裂事件进行24h连续监测，实时获取微震事件的动态时空数据等多项震源参数信息。
优选地，所述信号采集单元3为多个，当所述检波器2为单分量检波器时，每六个检波器2为一组，连接一个信号采集单元3；当所述检波器2为三分量检波器时，每两个检波器2为一组，连接一个信号采集单元3。
需要说明的是，每组检波器分布在所述回采工作面的回采巷道中，并根据被监测工作面的具体情况进行检波器及线路的布设，同时对检波器进行优化布置，以获得良好的空间阵列，减小定位误差。
参见图2，本实施例共采用五组单分量速度检波器2，每组单分量速度检波器2为六个，均设置在上述回采工作面的回采巷道1中。
通过上述微震监测系统，能够实时监测煤矿回采工作面区域微破裂的萌生和演化，获取煤矿采掘过程中突水通道的分布特征，有利于提高对回采工作面突水灾害分析预警的可靠性。
由上述技术方案可知，本发明提供的这种采动煤岩体突水微震监测系统，在回采工作面的回采巷道中布置预设数量的检波器，当被监测采动煤岩体内出现微震时，检波器可将微震信号拾取。由于检波器在回采巷中是多点设置的，所以通过采集多点同步数据测定各检波器接收到微震信号的时刻，连同检波器坐标及所测微震信号的波速带入方程组求解，即可确定震源的时空参数，达到定位的目的。相比现有技术，本发明提供的这种采动煤岩体突水的微震监测系统，部署简单，监测结果可靠性高。
本发明不局限于上述最佳实施方式，任何人在本发明的启示下都可得出其他各种形式的产品，但不论在其形状或结构上作任何变化，凡是具有与本申请相同或相近似的技术方案，均落在本发明的保护范围之内。术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性。术语“多个”指两个或两个以上，除非另有明确的限定。
Claims (9)
1.一种采动煤岩体突水微震监测方法，其特征在于，包括：
步骤S1、根据采动煤岩体突水微震监测系统的工程要求及实施条件，确定采动煤岩体突水微震监测的系统架构；其中，所述系统架构包括检波器阵列在煤矿回采工作面的回采巷道中的安装方式及安装位置，所述检波器阵列用于拾取被监测采动煤岩体产生的微震信号；
步骤S2、对所述检波器阵列拾取的微震信号，采用阈值算法和STA/LTA算法判断是否存在采动煤岩体突水微震触发事件，若是，保存触发事件波形；
步骤S3、滤除所述触发事件波形中的噪声信号，对滤噪后的触发事件波形，采用STA/LTA算法拾取p波到时t，并根据采动煤岩体整体最优波速vopt通过单纯形法迭代求解出所述触发事件的发生位置；
步骤S4、根据公式计算所述触发事件发生时的微震辐射能量；其中，ρ为煤岩体密度；v为煤岩体弹性波速；R为震源到检波器阵列的距离；Jc为辐射能通量，通过质点速度谱在频域中积分获得；Fc为地震波辐射类型经验系数；
步骤S5、分析所述触发事件发生时微震辐射能量的时空分布特征和聚集程度，设定微震辐射能量释放阈值，以识别出采动煤岩体潜在突水灾害区域；
其中，所述步骤S1具体包括：
步骤S10、根据采动煤岩体突水微震监测系统的工程要求及实施条件，确定采动煤岩体突水微震监测系统的系统架构；
步骤S11、通过现场波速测试试验获取可供选择的岩体弹性波波速范围(vmin,vmax)，其中，vmin和vmax分别表示岩体弹性波波速的最小值和最大值；
步骤S12、对采动煤岩体现场进行定点爆破或敲击试验，人工激发弹性波，滤除所述人工激发弹性波中的噪声信号，并对滤除噪声信号后的人工激发弹性波，采用STA/LTA算法拾取p波到时t；
步骤S13、遍历岩体弹性波波速vi∈(vmin,vmax)内的所有弹性波波速，采用单纯形法迭代求解微震震源位置Xi(xi,yi,zi,vi)，并利用公式minE(vi)＝||Xi-X0||，i＝1，2，......n，n≥2(2)，计算出采用单纯形法迭代求解的震源位置Xi(xi,yi,zi,vi)与人工震源位置X0(x0,y0,z0)之间的误差最小值min E(vi)，其对应的弹性波波速为采动煤岩体整体最优波速vopt；
步骤S14、若min E(vi)的值不在预设误差范围内，返回步骤S10，调整采动煤岩体突水微震监测系统的系统架构，直至min E(vi)的值在预设误差范围；其中，调整采动煤岩体突水微震监测系统的系统架构包括调整采动煤岩体突水微震监测系统中检波器阵列的安装方式及安装位置。
2.根据权利要求1所述的采动煤岩体突水微震监测方法，其特征在于，还包括：
步骤S6、根据所述触发事件发生时微震辐射能量的时空分布特征及采动煤岩体潜在突水灾害区域，确定采动煤岩体潜在突水通道赋存的主控地质因素，建立煤岩体突水灾害分析的三维地质力学模型，为采动煤岩体突水灾害的监测预警提供模型支持；
步骤S7、根据所述触发事件的时空分布规律和所述触发事件发生时微震辐射能量的时空分布特征，利用能量耗散原理，将煤矿围岩微震辐射能量转化为围岩损伤弱化的力学参数，为采动煤岩体的数值分析提供力学修正参数。
3.根据权利要求1所述的采动煤岩体突水微震监测方法，其特征在于，所述步骤S2具体包括：
步骤S21、采集所述检波器阵列输出的微震信号，判断所述微震信号初至时的幅值As是否大于第一触发阈值At；其中，At根据采动煤岩体突水微震监测系统工程现场的背景噪声值确定；
步骤S22、若As大于At，判断所述微震信号是否满足值，若是，则判定所述微震信号存在采动煤岩体突水微震触发事件；其中，STA为微震信号短时能量平均值，LTA为微震信号长时能量平均值，T为第二触发阈值。
4.根据权利要求3所述的采动煤岩体突水微震监测方法，其特征在于，所述步骤S21中At根据如下步骤进行确定：
步骤S210、计算预设时段内采动煤岩体突水微震监测系统工程现场的背景噪声的平均振幅值；
步骤S211、若微震信号初至时的幅值As大于所述平均振幅值，则判定所述微震信号中可能存在采动煤岩体突水微震触发事件；
步骤S212、若被监测采动煤岩体实际状态稳定，利用步骤S211判定的可能存在的触发事件每分钟低于30个，且微震信号的信噪比在预设比值范围内，则将所述平均振幅值设为第一触发阈值At；否则，返回步骤S210，调整预设时段的长度重新计算所述平均振幅。
5.根据权利要求1所述的采动煤岩体突水微震监测方法，其特征在于，所述步骤S3中滤除所述触发事件波形中的噪声信号具体为：
对工程现场大于预设幅值的噪声信号进行快速傅里叶变换，以获取噪声信号的时域和频域特征，建立噪声信号特征库；根据噪声信号特征库中噪声信号的时域和频域特征滤除所述触发事件波形中的噪声信号。
6.根据权利要求2所述的采动煤岩体突水微震监测方法，其特征在于，所述步骤S7具体包括：
根据公式计算微震触发事件发生后岩体的损伤变量D；根据损伤变量D，对采动煤岩体的力学参数进行修正；
其中，η地震效率，E表示震源尺寸范围内的岩体单元损伤后被微震监测系统拾取的辐射能，Ue为岩体单元可释放总应变能，Ue由公式
计算，其中E0为岩体初始的弹性模量，σ1、σ2和σ3分别为岩体第一、第二和第三主应力，ν为岩体泊松比。
7.根据权利要求1所述的采动煤岩体突水微震监测方法，其特征在于，所述步骤S1中采动煤岩体突水微震监测的系统架构包括硬件架构和软件架构，其硬件架构包括：
检波器阵列，包括预设数量的检波器，所述检波器阵列设置在回采工作面的回采巷道中；
信号采集单元，包括：集线器、模数转换器和存储单元，其中，所述模数转换器的输入端通过所述集线器与所述检波器阵列连接，输出端与所述存储单元连接；
信号分析单元，包括：处理子单元和分析子单元，其中，所述处理子单元的输入端与所述存储单元连接，输出端与所述分析子单元连接；所述处理子单元用于根据存储单元中的微震信号，计算震源参数，所述震源参数包括微震发生时刻、微震辐射能量和震源位置；所述分析子单元用于根据处理子单元输出的震源参数，分析采动煤岩体微震事件和微震能量时空分布规律；
预警决策单元，与所述分析子单元连接，用于根据所述采动煤岩体微震事件和微震能量时空分布规律进行突水灾害预警。
8.根据权利要求7所述的采动煤岩体突水微震监测方法，其特征在于，
所述检波器为单分量和/或三分量检波器，所述单分量检波器包括单分量加速度检波器和/或单分量速度检波器，所述三分量检波器包括三分量加速度检波器和/或三分量速度检波器。
9.根据权利要求8所述的采动煤岩体突水微震监测方法，其特征在于，
所述信号采集单元为多个，当所述检波器为单分量检波器时，每六个检波器为一组，连接一个信号采集单元；当所述检波器为三分量检波器时，每两个检波器为一组，连接一个信号采集单元。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610513797.4A CN105891874B (zh) | 2016-06-30 | 2016-06-30 | 一种采动煤岩体突水微震监测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610513797.4A CN105891874B (zh) | 2016-06-30 | 2016-06-30 | 一种采动煤岩体突水微震监测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105891874A CN105891874A (zh) | 2016-08-24 |
| CN105891874B true CN105891874B (zh) | 2017-12-12 |
Family
ID=56718661
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610513797.4A Active CN105891874B (zh) | 2016-06-30 | 2016-06-30 | 一种采动煤岩体突水微震监测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105891874B (zh) |
Families Citing this family (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107356510B (zh) * | 2017-08-11 | 2022-04-22 | 山东科技大学 | 模拟采场围岩突水通道扩展演化试验系统 |
| CN107478725B (zh) * | 2017-08-31 | 2020-05-19 | 北京市政建设集团有限责任公司 | 一种公路隧道中夹岩隔墙稳定性评价方法 |
| CN108387948B (zh) * | 2018-01-23 | 2020-08-04 | 华北科技学院 | 煤层底板水危险源区域定位方法及装置 |
| CN108319915B (zh) * | 2018-01-31 | 2020-06-19 | 中国科学院武汉岩土力学研究所 | 一种岩爆信号阈值动态调整的多时窗简化形式识别方法 |
| CN109100787A (zh) * | 2018-06-29 | 2018-12-28 | 安徽万泰地球物理技术有限公司 | 一种等离子体脉冲谐振增透的微震监测评价方法 |
| CN109521467A (zh) * | 2018-11-26 | 2019-03-26 | 阳泉煤业（集团）股份有限公司 | 一种基于突出煤层巷道的超前探测方法 |
| CN110118991B (zh) * | 2019-05-16 | 2020-06-23 | 中国矿业大学 | 一种基于微震损伤重构的采动应力评估方法 |
| CN110687590A (zh) * | 2019-09-20 | 2020-01-14 | 贵州正厦工程检测技术咨询有限公司 | 一种用于微震监测的无线智能数据传输系统及方法 |
| CN110910613B (zh) * | 2019-12-10 | 2022-04-05 | 大连理工大学 | 一种岩体微震无线监测接收预警系统 |
| CN111816204B (zh) * | 2020-06-18 | 2022-09-13 | 山西宏安翔科技股份有限公司 | 一种三分量拾音系统 |
| CN112799129B (zh) * | 2021-01-14 | 2022-08-02 | 华北科技学院 | 一种导水通道的识别方法及系统 |
| CN112785182A (zh) * | 2021-01-29 | 2021-05-11 | 河北煤炭科学研究院有限公司 | 承压水采煤隐伏导水通道空间分析方法、装置及终端设备 |
| CN113777272B (zh) * | 2021-08-30 | 2024-01-19 | 安徽理工大学 | 用于智能加载多维相似模型试验的多场源监测及分析系统 |
| CN114137605A (zh) * | 2021-10-21 | 2022-03-04 | 陕西延长石油矿业有限责任公司 | 一种煤矿岩体中岩桥的智能定位及特征识别方法 |
| CN114154215A (zh) * | 2021-11-23 | 2022-03-08 | 中铁十四局集团第一工程发展有限公司 | 基于机理预测与误差修正的隧道突涌水预测方法及系统 |
| CN116797020A (zh) * | 2023-05-24 | 2023-09-22 | 中国矿业大学 | 一种考虑岩层结构演化的煤矿顶板离层突水微震预警方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103529488A (zh) * | 2013-09-30 | 2014-01-22 | 武汉长盛煤安科技有限公司 | 矿井顶底板突水监测预报系统及方法 |
| CN104989453A (zh) * | 2015-06-19 | 2015-10-21 | 河北煤炭科学研究院 | 煤矿突水井下全空间实时连续监测预警方法 |
| CN105652311A (zh) * | 2016-01-18 | 2016-06-08 | 陕西煤业化工技术研究院有限责任公司 | 一种监测底板突水的微震监测方法 |
-
2016
- 2016-06-30 CN CN201610513797.4A patent/CN105891874B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103529488A (zh) * | 2013-09-30 | 2014-01-22 | 武汉长盛煤安科技有限公司 | 矿井顶底板突水监测预报系统及方法 |
| CN104989453A (zh) * | 2015-06-19 | 2015-10-21 | 河北煤炭科学研究院 | 煤矿突水井下全空间实时连续监测预警方法 |
| CN105652311A (zh) * | 2016-01-18 | 2016-06-08 | 陕西煤业化工技术研究院有限责任公司 | 一种监测底板突水的微震监测方法 |
Non-Patent Citations (4)
| Title |
|---|
| THE PARTITION OF RADIATED ENERGY BETWEEN P AND S WAVES;JOHN BOATWRIGHT et al.;《Bulletin of the Seismological Society of America》;19840430;第74卷(第2期);第361-376页 * |
| 基于微震监测的岩质边坡稳定性三维反馈分析;徐奴文等;《岩石力学与工程学报》;20140131;第33卷;第3096-3100页、图5-8 * |
| 微震监测技术在大岗山水电站右岸边坡中的应用;徐奴文等;《防灾减灾工程学报》;20100930;第30卷;第217-218页，附图4 * |
| 高精度微震监测技术在煤矿突水监测中的应用;姜福兴等;《岩石力学与工程学报》;20080930;第27卷(第9期);第1933-1977页，图1、4-5 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105891874A (zh) | 2016-08-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105891874B (zh) | 一种采动煤岩体突水微震监测方法 | |
| Ghosh et al. | Application of underground microseismic monitoring for ground failure and secure longwall coal mining operation: a case study in an Indian mine | |
| CA2562005C (en) | Microseismic fracture mapping using seismic source timing measurements for velocity calibration | |
| CN104390537B (zh) | 一种基于爆破振动测试的边坡预裂爆破开挖损伤控制方法 | |
| Zeng et al. | Surface microseismic monitoring of hydraulic fracturing of a shale‐gas reservoir using short‐period and broadband seismic sensors | |
| CN106199700B (zh) | 一种地下水封储油洞库微震监测方法及系统 | |
| CN106501848B (zh) | 一种隧道掘进过程中隐性断层超前物探方法 | |
| CN104880729B (zh) | 一种基于连续跟踪槽波信号的煤巷超前探测异常构造方法 | |
| CN113189644B (zh) | 一种微震震源定位方法及系统 | |
| CN110516650A (zh) | 一种基于震动传感器阵列的浅层盲空间震源定位系统 | |
| CN113153430B (zh) | 巷道围岩破坏声发射定位与波速成像监测及灾变预警方法 | |
| CN105607040A (zh) | 一种矿区防盗采监测定位方法及系统 | |
| US20080112263A1 (en) | System and method for determining seismic event location | |
| CN104749637B (zh) | 随钻地震钻头震源侧帮地质构造探测方法 | |
| CN109001810A (zh) | 基于微震监测的重力坝坝体裂缝超前预警方法 | |
| EP2951393B1 (en) | Method for determining fracture network volume using passive seismic signals | |
| US9416641B2 (en) | Borehole microseismic systems and methods | |
| NO20141369A1 (no) | System og fremgangsmåte for overvåking av mikroseismiske hendelser i en undergrunnstruktur | |
| CN106154321B (zh) | 射孔信号的检测方法和装置 | |
| US8960280B2 (en) | Method for determining fracture plane orientation using passive seismic signals | |
| Hu et al. | Quantitative comparative analysis of monitoring effect of microseismic fracturing borehole and surface | |
| Das et al. | Overview of seismic monitoring and assessment of seismic hazard based on a decade of seismic events | |
| EA022596B1 (ru) | Способ контроля и определения опасного состояния массива горных пород при подземных горных работах | |
| Dias et al. | A microseismic monitoring pilot project of natural caves in Carajás–PA | |
| CN118519209A (zh) | 一种基于震电联合的疲劳压裂区各向异性探测方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20171109 Address after: 116024 Liaoning, Dalian, Ganjingzi Ling Road, No. 2 Applicant after: Dalian University of Technology Address before: 100000 Institute of mechanics, Chinese Academy of Sciences, 15 West Fourth Ring Road, Beijing, Haidian District Applicant before: Ma Ke Applicant before: Zhuang Duanyang Applicant before: Tang Chunan |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |