CN118688716A - A method for improving positioning accuracy of coal mine microseismic monitoring system based on wave velocity zone division - Google Patents
A method for improving positioning accuracy of coal mine microseismic monitoring system based on wave velocity zone division Download PDFInfo
- Publication number
- CN118688716A CN118688716A CN202410850990.1A CN202410850990A CN118688716A CN 118688716 A CN118688716 A CN 118688716A CN 202410850990 A CN202410850990 A CN 202410850990A CN 118688716 A CN118688716 A CN 118688716A
- Authority
- CN
- China
- Prior art keywords
- wave velocity
- sensor
- zone
- monitoring system
- wave
- 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.)
- Pending
Links
- 238000012544 monitoring process Methods 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 title claims abstract description 39
- 239000003245 coal Substances 0.000 title claims abstract description 33
- 238000005065 mining Methods 0.000 claims abstract description 10
- 238000005422 blasting Methods 0.000 claims abstract description 9
- 239000007787 solid Substances 0.000 claims description 8
- 238000005259 measurement Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 5
- 239000002360 explosive Substances 0.000 claims description 3
- 230000035939 shock Effects 0.000 claims description 3
- 239000002689 soil Substances 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 10
- 239000011435 rock Substances 0.000 description 20
- 101001121408 Homo sapiens L-amino-acid oxidase Proteins 0.000 description 2
- 102100026388 L-amino-acid oxidase Human genes 0.000 description 2
- 230000000052 comparative effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 101100012902 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) FIG2 gene Proteins 0.000 description 1
- 239000002817 coal dust Substances 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- -1 fire Substances 0.000 description 1
- 238000011017 operating method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/18—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/18—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves
- G01S5/22—Position of source determined by co-ordinating a plurality of position lines defined by path-difference measurements
-
- 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
- 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/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Geophysics (AREA)
- Radar, Positioning & Navigation (AREA)
- Business, Economics & Management (AREA)
- Emergency Management (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明属于煤矿安全技术领域,具体涉及一种基于波速区划分的煤矿微震监测系统定位精度提高方法,包括以下步骤:步骤1:在目标矿井布置微震监测系统,在井下和地面布置拾振传感器,并测量出各传感器坐标;步骤2:根据矿井地质和开采技术条件以及震动波在不同区域波速值的不同进行区域划分,根据区域划分结果构建三维波速区划分模型;步骤3:采用放炮震源的方式计算出各波速区波速值;步骤4:首先采用实体区波速V1作为系统各传感器初始波速值对震源进行定位,计算出初始位置0(X0、Y0、Z0);根据震源初始位置和各拾振传感器位置,按照其两者中间所穿过的不同波速区的距离,计算出各传感器波速;步骤5:采用新波速值重新计算震源位置。本发明震源定位更加精确。
The present invention belongs to the field of coal mine safety technology, and specifically relates to a method for improving the positioning accuracy of a coal mine microseismic monitoring system based on wave velocity zone division, comprising the following steps: step 1: arranging a microseismic monitoring system in a target mine, arranging vibration pickup sensors underground and on the ground, and measuring the coordinates of each sensor; step 2: performing regional division according to the geological and mining technical conditions of the mine and the different wave velocity values of the vibration wave in different regions, and constructing a three-dimensional wave velocity zone division model according to the regional division results; step 3: calculating the wave velocity values of each wave velocity zone by using a blasting source method; step 4: first using the entity zone wave velocity V1 as the initial wave velocity value of each sensor of the system to locate the source, and calculating the initial position 0 ( X0 , Y0 , Z0 ); according to the initial position of the source and the position of each vibration pickup sensor, according to the distance of the different wave velocity zones passed between the two, the wave velocity of each sensor is calculated; step 5: recalculating the source position by using the new wave velocity value. The source positioning of the present invention is more accurate.
Description
技术领域Technical Field
本发明属于煤矿安全技术领域,具体涉及一种基于波速区划分的煤矿微震监测系统定位精度提高方法。The invention belongs to the technical field of coal mine safety, and in particular relates to a method for improving the positioning accuracy of a coal mine microseismic monitoring system based on wave velocity zone division.
背景技术Background Art
冲击地压是一种严重的煤矿灾害,已继水、火、瓦斯、煤尘、顶板之后成为煤矿第六大灾害,且随着煤矿开采深度的不断增加,有成为煤矿最主要灾害的趋势,近些年来,冲击地压灾害频发,冲击地压监测技术是防治冲击地压灾害的首要和关键环节,《防治煤矿冲击地压细则》第四十六条规定:冲击地压矿井必须建立区域与局部相结合的冲击危险性监测制度,区域监测应当覆盖矿井采掘区域,区域监测可采用微震监测法等;煤岩体在受力达到一定程度时会产生裂隙,裂隙在产生和扩展以及岩块运动的过程中会释放震动波,微震监测技术就是一种通过监测震动波来定位煤岩体破裂或运动位置及计算释放能量大小的技术;微震监测技术是冲击地压矿井采用的最重要的监测技术,其中震源定位精度对冲击危险监测起着关键影响;微震监测技术通过一定的定位算法来定位震源位置,其中影响定位精度的关键是波速值的确定。当前微震监测技术均采用固定波速,但地层本身就是非均质体,且采矿活动会使地层赋存状态发生重大变化,当震动波穿过实体区、采空区、地质异常区、表土层等区域时波速差异较大,若采用固定波速定位,震源误差必然较大,将会严重影响微震监测技术对冲击危险的监测。Rock burst is a serious coal mine disaster. It has become the sixth largest disaster in coal mines after water, fire, gas, coal dust and roof. With the continuous increase in the depth of coal mining, it has the tendency to become the most important disaster in coal mines. In recent years, rock burst disasters have occurred frequently. Rock burst monitoring technology is the primary and key link in preventing and controlling rock burst disasters. Article 46 of the "Detailed Rules for Preventing and Controlling Rock Burst in Coal Mines" stipulates that rock burst mines must establish a regional and local impact hazard monitoring system. Regional monitoring should cover the mining area of the mine. Regional monitoring can adopt microseismic monitoring methods, etc.; when the force reaches a certain level, cracks will be generated in the coal rock mass. During the generation and expansion of the cracks and the movement of the rock blocks, vibration waves will be released. Microseismic monitoring technology is a technology that locates the fracture or movement position of the coal rock mass and calculates the amount of released energy by monitoring vibration waves; microseismic monitoring technology is the most important monitoring technology used in rock burst mines, among which the accuracy of earthquake source positioning plays a key role in monitoring impact hazard; microseismic monitoring technology locates the earthquake source through a certain positioning algorithm, among which the key to affecting the positioning accuracy is the determination of the wave velocity value. Current microseismic monitoring technologies all use a fixed wave velocity, but the strata themselves are heterogeneous, and mining activities will cause significant changes in the stratum's state of occurrence. When the vibration wave passes through solid areas, goaf areas, geological anomaly areas, topsoil layers and other areas, the wave velocity difference is large. If fixed wave velocity positioning is used, the source error will inevitably be large, which will seriously affect the microseismic monitoring technology's monitoring of impact hazards.
例如专利号为CN105022031A的中国专利公开了一种区域岩体微震震源的分层速度定位方法,步骤如下:①将待测震源的岩体区域划分为不同的波速层,建立三维直角坐标系,计算出各波速层的交界面的解析方程;②在岩体区域安装传感器并测量各传感器的空间坐标;③在所述岩体区域进行爆破实验,计算各波速层的岩体波速;④计算微震震源的空间坐标(x0,y0,z0)。本发明所述方法对微震震源的定位精度高,能够促进微震监测技术在工程实践中更好地发挥预测预警作用;但是现有技术中,将岩体划分为相互平行的不同波速层,且后续方法仅适用于层状波速区划分;该发明中要求必须在第一和最后一个波速层布置传感器;使得适用性能相对局部;为此,研究一种科学、合理的波速确定和定位方法对提高震源定位精度和冲击危险的监测能力具有重要意义。For example, the Chinese patent with patent number CN105022031A discloses a layered velocity positioning method for regional rock microseismic sources, and the steps are as follows: ① Divide the rock mass area of the source to be measured into different wave velocity layers, establish a three-dimensional rectangular coordinate system, and calculate the analytical equation of the interface of each wave velocity layer; ② Install sensors in the rock mass area and measure the spatial coordinates of each sensor; ③ Carry out blasting experiments in the rock mass area to calculate the rock mass wave velocity of each wave velocity layer; ④ Calculate the spatial coordinates (x0, y0, z0) of the microseismic source. The method of the present invention has high positioning accuracy for the microseismic source, which can promote the microseismic monitoring technology to better play a predictive and early warning role in engineering practice; However, in the prior art, the rock mass is divided into different wave velocity layers parallel to each other, and the subsequent method is only applicable to the layered wave velocity zone division; The invention requires that sensors must be arranged in the first and last wave velocity layers; The applicability is relatively local; Therefore, it is of great significance to study a scientific and reasonable wave velocity determination and positioning method to improve the accuracy of source positioning and the monitoring capability of impact danger.
发明内容Summary of the invention
本发明的目的是提供一种基于波速区划分的煤矿微震监测系统定位精度提高方法,本方法首先测量不同区域波速值,然后根据传感器及其与震源之间所穿过波速区的距离来自动计算各传感器波速值,而后根据新波速值进行震源定位并不断迭代,可以消除当前微震监测系统采用固定波速定位所带来的误差,进而大幅度提高煤矿微震系统定位精度。The purpose of the present invention is to provide a method for improving the positioning accuracy of a coal mine microseismic monitoring system based on wave velocity zone division. The method first measures the wave velocity values of different areas, and then automatically calculates the wave velocity value of each sensor based on the distance between the sensor and the wave velocity zone passed by the sensor and the source, and then locates the source according to the new wave velocity value and continuously iterates. This can eliminate the errors caused by the current microseismic monitoring system using fixed wave velocity positioning, thereby greatly improving the positioning accuracy of the coal mine microseismic system.
本发明采取的技术方案具体如下:The technical solution adopted by the present invention is as follows:
一种基于波速区划分的煤矿微震监测系统定位精度提高方法,包括以下步骤:A method for improving the positioning accuracy of a coal mine microseismic monitoring system based on wave velocity zone division comprises the following steps:
步骤1:在目标矿井布置微震监测系统,在井下和地面布置拾振传感器,其中地面传感器编号T1、T2、T3…Tn,井下传感器编号S1、S2、S3…Sn,所述地面传感器以及所述井下传感器总和大于等于四;并测量出各传感器坐标;Step 1: Arrange a microseismic monitoring system in the target mine, and arrange vibration pickup sensors underground and on the ground, where the ground sensors are numbered T 1 , T 2 , T 3 ..T n , and the underground sensors are numbered S 1 , S 2 , S 3 ..S n , and the total number of the ground sensors and the underground sensors is greater than or equal to four; and measure the coordinates of each sensor;
步骤2:根据矿井地质和开采技术条件,进行区域划分,将区域划分为呈实体区A1、垮落区A2、裂隙区A3、表土区A4以及地质构造区A5…An的不同波速区,在区域划分过程中,按照矿井实测进行划分,若无实测按矿压理论取值进行划分,根据区域划分结果构建三维波速区划分模型;Step 2: According to the geological and mining technical conditions of the mine, the region is divided into different wave velocity zones, namely, solid zone A1 , collapse zone A2 , fissure zone A3 , surface soil zone A4 and geological structure zone A5 ... An. In the process of regional division, the division is carried out according to the actual measurement of the mine. If there is no actual measurement, the division is carried out according to the theoretical value of the mine pressure. According to the regional division results, a three-dimensional wave velocity zone division model is constructed;
步骤3:采用放炮震源的方式计算出各波速区波速值,在波速区An一侧安放炸药,在另一侧安放传感器来测试波速,若受井下巷道布置,另一侧无紧挨波速区An的巷道来布置传感器,则需要选择距此区域最近的巷道布置传感器;例如,传感器与爆破点之间穿过了实体区距离D1和地质构造区距离D5,测出其波速为V15,其中实体区波速已知为V1,则地质构造区波速V5=[V15×(D1+D2)-D1×V1]/D5;Step 3: Calculate the wave velocity value of each wave velocity zone by blasting source, place explosives on one side of the wave velocity zone A n , and place sensors on the other side to test the wave velocity. If there is no roadway close to the wave velocity zone A n on the other side to arrange the sensor due to the underground tunnel layout, it is necessary to select the roadway closest to this area to arrange the sensor; for example, the distance between the sensor and the blasting point passes through the physical area distance D 1 and the geological structure area distance D 5 , and its wave velocity is measured to be V 15 , where the wave velocity in the physical area is known to be V 1 , then the wave velocity in the geological structure area V 5 = [V 15 × (D 1 +D 2 )-D 1 ×V 1 ]/D 5 ;
步骤4:首先采用实体区波速V1作为系统各传感器初始波速值对震源进行定位,计算出初始位置0(X0、Y0、Z0);根据震源初始位置和各拾振传感器位置,按照其两者中间所穿过的不同波速区的距离,计算出各传感器波速,震源与地面传感器Tn之间依次穿过A1、A2、A3、…、An;则传感器T1波速应设置为式中,Li为震动波穿过对应区域的距离,Vi是对应区域的速率;Step 4: First, use the entity area wave velocity V1 as the initial wave velocity value of each sensor in the system to locate the earthquake source and calculate the initial position 0 ( X0 , Y0 , Z0 ); according to the initial position of the earthquake source and the position of each vibration pickup sensor, according to the distance between the two passing through different wave velocity zones, calculate the wave velocity of each sensor. The earthquake source and the ground sensor Tn pass through A1 , A2 , A3 , ..., An in turn; then the wave velocity of sensor T1 should be set to Where Li is the distance that the shock wave passes through the corresponding area, and Vi is the velocity of the corresponding area;
步骤5:采用新波速值重新计算震源位置,不断迭代直至震源位置与上一次震源位置误差小于1m计算终止。Step 5: Use the new wave velocity value to recalculate the source position, and continue iterating until the error between the source position and the previous source position is less than 1m. The calculation is terminated.
本发明取得的技术效果为:The technical effects achieved by the present invention are:
由于微震震动波速穿过垮落带、裂隙带、表土带、地质构造带、煤岩实体等不同区域时波速差异较大,但现行微震监测系统的采用固定波速进行定位,必然存在很大的误差;本发明首先测量不同区域波速值,然后根据传感器及其与震源之间所穿过波速区的距离来自动计算各传感器波速值,而后根据新波速值进行震源定位并不断迭代,可以消除当前微震监测系统采用固定波速定位所带来的误差,进而大幅度提高煤矿微震系统定位精度。Since the wave velocity of microseismic vibration varies greatly when passing through different areas such as collapse zones, fissure zones, topsoil zones, geological structural zones, coal and rock entities, the current microseismic monitoring system uses a fixed wave velocity for positioning, which inevitably results in large errors. The present invention first measures the wave velocity values of different areas, and then automatically calculates the wave velocity value of each sensor based on the distance between the sensor and the wave velocity zone passed through by the sensor and the source, and then locates the source based on the new wave velocity value and continuously iterates, which can eliminate the errors caused by the current microseismic monitoring system using a fixed wave velocity for positioning, thereby greatly improving the positioning accuracy of the coal mine microseismic system.
本发明与现有技术相比,波速区划分不同;现有技术中将岩体划分为相互平行的不同波速层,且后续方法仅适用于层状波速区划分;本发明可根据矿井地质及开采技术条件划分为任意形状的波速区,且后续方法适用该划分,与现有技术相比传感器布置方式不同;现有技术要求必须在第一和最后一个波速层布置传感器;而本发明对传感器布置位置无特殊要求,只需要在目标矿井布置微震监测系统,在井下和地面布置拾振传感器,其中地面传感器编号T1、T2、T3…Tn,井下传感器编号S1、S2、S3…Sn,所述地面传感器以及所述井下传感器总和大于等于四;并测量出各传感器坐标;与现有技术相比,本发明的定位方式不同;现有技术中对比文件通过假设震源位于某一波速层中,通过求解目标函数最小值确定震源位置;而本发明首先采用实体煤波速计算震源位置,而后根据震源和传感器之间所穿过波速区位置不断重新计算传感器波速值,而后不断迭代直至与上次震源位置误差小于1m为止,即为真实震源位置,使得最终得到的真实震源位置更加精确。Compared with the prior art, the present invention has a different wave velocity zone division; in the prior art, the rock mass is divided into different wave velocity layers parallel to each other, and the subsequent method is only applicable to layered wave velocity zone division; the present invention can be divided into wave velocity zones of any shape according to the mine geology and mining technical conditions, and the subsequent method is applicable to the division, and the sensor arrangement method is different from the prior art; the prior art requires that sensors must be arranged in the first and last wave velocity layers; while the present invention has no special requirements for the sensor arrangement position, and only needs to arrange a microseismic monitoring system in the target mine, and arrange vibration pickup sensors underground and on the ground, wherein the ground sensors are numbered T 1 , T 2 , T 3 ... T n , and the underground sensors are numbered S 1 , S 2 , S 3 ... S n , the sum of the ground sensors and the downhole sensors is greater than or equal to four; and the coordinates of each sensor are measured; compared with the prior art, the positioning method of the present invention is different; the comparative document in the prior art assumes that the earthquake source is located in a certain wave velocity layer, and determines the earthquake source position by solving the minimum value of the objective function; and the present invention first uses the wave velocity of the physical coal to calculate the earthquake source position, and then continuously recalculates the sensor wave velocity value according to the position of the wave velocity zone passed between the earthquake source and the sensor, and then continuously iterates until the error with the last earthquake source position is less than 1m, that is, the real earthquake source position, so that the final real earthquake source position is more accurate.
附图说明BRIEF DESCRIPTION OF THE DRAWINGS
图1是本发明一种基于波速区划分的煤矿微震监测系统定位精度提高方法的流程图;FIG1 is a flow chart of a method for improving positioning accuracy of a coal mine microseismic monitoring system based on wave velocity zone division according to the present invention;
图2是本发明一种基于波速区划分的煤矿微震监测系统定位精度提高方法中波速区划分及震源定位示意图。2 is a schematic diagram of wave velocity zone division and earthquake source positioning in a method for improving positioning accuracy of a coal mine microseismic monitoring system based on wave velocity zone division according to the present invention.
附图中,各标号所代表的部件列表如下:In the accompanying drawings, the components represented by the reference numerals are listed as follows:
1、表土区;2、地面微震传感器T1;3、实体区;4、震源初始位置0;5、裂隙区;6、震源最终位置i+1;7、垮落区;8、井下传感器S1;9、井下传感器S2;10、井下传感器S3;11、地质构造区;12、井下巷道。1. Topsoil area; 2. Ground microseismic sensor T1; 3. Solid area; 4. Initial earthquake source position 0; 5. Fracture area; 6. Final earthquake source position i+1; 7. Collapse area; 8. Downhole sensor S1 ; 9. Downhole sensor S2 ; 10. Downhole sensor S3 ; 11. Geological structure area; 12. Underground tunnel.
具体实施方式DETAILED DESCRIPTION
为了使本发明的目的及优点更加清楚明白,以下结合实施例对本发明进行具体说明。应当理解,以下文字仅仅用以描述本发明的一种或几种具体的实施方式,并不对本发明具体请求的保护范围进行严格限定。In order to make the purpose and advantages of the present invention more clearly understood, the present invention is specifically described below in conjunction with embodiments. It should be understood that the following text is only used to describe one or several specific embodiments of the present invention, and does not strictly limit the scope of protection of the specific claims of the present invention.
如图1所示,一种基于波速区划分的煤矿微震监测系统定位精度提高方法,包括以下步骤:As shown in FIG1 , a method for improving the positioning accuracy of a coal mine microseismic monitoring system based on wave velocity zone division includes the following steps:
步骤1:如图2所示,图为本发明波速区划分及震源定位示意图,在目标矿井布置微震监测系统,在井下和地面布置拾振传感器,其中地面传感器编号T1、T2、T3…Tn,井下传感器编号S1、S2、S3…Sn,所述地面传感器以及所述井下传感器总和大于等于四;并测量出各传感器坐标;Step 1: As shown in FIG2 , which is a schematic diagram of wave velocity zone division and earthquake source location of the present invention, a microseismic monitoring system is arranged in the target mine, and vibration pickup sensors are arranged underground and on the ground, wherein the ground sensors are numbered T 1 , T 2 , T 3 ..T n , and the underground sensors are numbered S 1 , S 2 , S 3 ..S n , and the total number of the ground sensors and the underground sensors is greater than or equal to four; and the coordinates of each sensor are measured;
步骤2:根据矿井地质和开采技术条件,进行区域划分,将区域划分为呈实体区A1、垮落区A2、裂隙区A3、表土区A4以及地质构造区A5的不同波速区,在区域划分过程中,按照矿井实测进行划分,若无实测按矿压理论取值进行划分,根据区域划分结果构建三维波速区划分模型;Step 2: According to the geological and mining technical conditions of the mine, the region is divided into different wave velocity zones, namely, solid zone A1 , collapse zone A2 , fissure zone A3 , surface soil zone A4 and geological structure zone A5. In the process of regional division, the division is carried out according to the actual measurement of the mine. If there is no actual measurement, the division is carried out according to the theoretical value of the mine pressure. According to the regional division results, a three-dimensional wave velocity zone division model is constructed;
步骤3:采用放炮震源的方式计算出各波速区波速值,在地质构造区一侧安放炸药,在另一侧安放传感器来测试波速,若受井下巷道布置,另一侧无紧挨地质构造区的巷道来布置传感器,则需要选择距此区域最近的巷道布置传感器,传感器与爆破点之间穿过了实体区距离D1和地质构造区距离D5,测出其波速为V15,其中实体区波速已知为V1,则地质构造区波速V5=[V15×(D1+D2)-D1×V1]/D5;Step 3: Calculate the wave velocity value of each wave velocity zone by blasting source, place explosives on one side of the geological structure zone, and place sensors on the other side to test the wave velocity. If there is no tunnel close to the geological structure zone on the other side to arrange the sensor due to the underground tunnel layout, it is necessary to select the tunnel closest to this area to arrange the sensor. The distance between the sensor and the blasting point passes through the physical zone distance D1 and the geological structure zone distance D5 . The wave velocity is measured to be V15 , where the wave velocity of the physical zone is known to be V1 , then the wave velocity of the geological structure zone V5 = [ V15 × ( D1 + D2 ) - D1 × V1 ] / D5 ;
步骤4:首先采用实体区波速V1作为系统各传感器初始波速值对震源进行定位,计算出初始位置0(X0、Y0、Z0);根据震源初始位置和各拾振传感器位置,按照其两者中间所穿过的不同波速区的距离,计算出各传感器波速,震源与地面传感器Tn之间依次穿过A1、A2、A3、…、An;则传感器T1波速应设置为式中,Li为震动波穿过对应区域的距离,Vi是对应区域的速率;Step 4: First, use the entity area wave velocity V1 as the initial wave velocity value of each sensor in the system to locate the earthquake source and calculate the initial position 0 ( X0 , Y0 , Z0 ); according to the initial position of the earthquake source and the position of each vibration pickup sensor, according to the distance between the two passing through different wave velocity zones, calculate the wave velocity of each sensor. The earthquake source and the ground sensor Tn pass through A1 , A2 , A3 , ..., An in turn; then the wave velocity of sensor T1 should be set to Where Li is the distance that the shock wave passes through the corresponding area, and Vi is the velocity of the corresponding area;
步骤5:采用新波速值重新计算震源位置,不断迭代直至震源位置与上一次震源位置误差小于1m计算终止;具体的当时,计算终止,输出震源位置,(Xi+1,Yi+1,Zi+1)。Step 5: Use the new wave velocity value to recalculate the source position, and continue to iterate until the error between the source position and the previous source position is less than 1m. When , the calculation terminates and the source position is output, (X i+1 , Yi +1 , Zi +1 ).
本方法首先测量不同区域波速值,然后根据传感器及其与震源之间所穿过波速区的距离来自动计算各传感器波速值,而后根据新波速值进行震源定位并不断迭代,可以消除当前微震监测系统采用固定波速定位所带来的误差,进而大幅度提高煤矿微震系统定位精度。This method first measures the wave velocity values in different areas, and then automatically calculates the wave velocity value of each sensor based on the distance between the sensor and the wave velocity zone passed by the sensor and the source. Then, the source is located according to the new wave velocity value and it is continuously iterated. This can eliminate the error caused by the current microseismic monitoring system that uses fixed wave velocity positioning, thereby greatly improving the positioning accuracy of the coal mine microseismic system.
本发明的工作原理为:The working principle of the present invention is:
步骤1:例如某矿为冲击地压矿井,矿井布置了井地一体微震监测系统。在某工作面周边共有4个传感器,其中井下传感器三个,编号分别为S1、S2、S3,地面传感器1个,编号为T1,测量出了传感器三维坐标,具体位置如图2所示;Step 1: For example, a mine is a rock burst mine, and a mine-ground integrated microseismic monitoring system is deployed in the mine. There are 4 sensors around a certain working face, including three underground sensors, numbered S 1 , S 2 , and S 3 , and one ground sensor, numbered T 1 . The three-dimensional coordinates of the sensors are measured, and the specific positions are shown in Figure 2;
步骤2:根据该区域地质和开采技术条件,划分了不同波速区,分别为实体区A1、垮落区A2、裂隙区A3、表土区A4,地质构造区A5;根据矿井实测和地质勘探结果对波速区范围进行了确定,并进一步构建了三维波速区模型,如图2所示;Step 2: According to the geological and mining technical conditions of the region, different wave velocity zones were divided, namely solid zone A1 , collapse zone A2 , fissure zone A3 , topsoil zone A4 , and geological structure zone A5 ; the range of the wave velocity zone was determined according to the actual measurement of the mine and the geological exploration results, and a three-dimensional wave velocity zone model was further constructed, as shown in Figure 2;
步骤3:采用放炮震源的方式测出了各波速区波速值,实体区V1为4130m/s、垮落区V2为2470m/s、裂隙区V3为3420m/s、表土区V4为2650m/s、地质构造区V5为3750m/s;Step 3: The velocity values of each velocity zone were measured by blasting source, with V1 in the solid zone being 4130 m/s, V2 in the collapse zone being 2470 m/s, V3 in the fissure zone being 3420 m/s, V4 in the topsoil zone being 2650 m/s, and V5 in the geological structure zone being 3750 m/s.
步骤4:首先采用实体区波速V1作为系统各传感器初始波速值对震源进行定位,计算出初始位置0(X0、Y0、Z0),此处根据多个波速计算震源初始位置的计算方法为本领域人员公知常识,在此不做赘述,位置如图2所示;根据震源初始位置和各传感器位置,按照其两者中间所穿过的不同波速区的距离,计算出各传感器波速;如初始位置0与地面传感器T1之间依次穿过裂隙区、实体区、表土区,穿过的距离分别为L1=60m、L2=250m、L3=50m,则传感器T1波速应设置为VT1-1=(60×3420+250×4130+50×2650)/(60+250+50)=3806.11m/s,其他传感器同理。Step 4: First, the solid area wave velocity V1 is used as the initial wave velocity value of each sensor in the system to locate the earthquake source, and the initial position 0 ( X0 , Y0 , Z0 ) is calculated. The calculation method for calculating the initial position of the earthquake source based on multiple wave velocities is common knowledge in this field and will not be repeated here. The position is shown in Figure 2; according to the initial position of the earthquake source and the positions of each sensor, the wave velocity of each sensor is calculated according to the distance between the initial position of the earthquake source and the positions of each sensor and the different wave velocity zones passed between them; if the distance between the initial position 0 and the ground sensor T1 is L1 = 60m, L2 = 250m, and L3 = 50m respectively, then the wave velocity of sensor T1 should be set to VT1-1 = (60×3420+250×4130+50×2650)/(60+250+50) = 3806.11m/s, and the same applies to other sensors.
步骤5:采用新波速值重新计算震源位置1(X1、Y1、Z1),根据新震源位置与传感器位置重新计算传感器波速,而后再次定位,不断迭代直至震源位置与上一次震源位置误差小于1m计算终止,流程如2所示。Step 5: Use the new wave velocity value to recalculate the source position 1 (X 1 , Y 1 , Z 1 ), recalculate the sensor wave velocity based on the new source position and the sensor position, and then reposition, and continue to iterate until the error between the source position and the previous source position is less than 1m. The calculation is terminated. The process is shown in 2.
本发明中,微震震动波速穿过垮落带、裂隙带、表土带、地质构造带、煤岩实体等不同区域时波速差异较大,但现行微震监测系统的采用固定波速进行定位,必然存在很大的误差。本方法首先测量不同区域波速值,然后根据传感器及其与震源之间所穿过波速区的距离来自动计算各传感器波速值,而后根据新波速值进行震源定位并不断迭代,可以消除当前微震监测系统采用固定波速定位所带来的误差,进而大幅度提高煤矿微震系统定位精度。In the present invention, the velocity of microseismic vibration waves varies greatly when passing through different areas such as collapse zones, fracture zones, topsoil zones, geological structural zones, coal and rock entities, but the current microseismic monitoring system uses fixed wave velocities for positioning, which inevitably results in large errors. This method first measures the wave velocity values of different areas, and then automatically calculates the wave velocity values of each sensor based on the distance between the sensor and the wave velocity zone passed through by the source, and then locates the source based on the new wave velocity value and continuously iterates, which can eliminate the errors caused by the current microseismic monitoring system using fixed wave velocity positioning, thereby greatly improving the positioning accuracy of the coal mine microseismic system.
本发明中,与现有技术相比波速区划分不同;现有技术中将岩体划分为相互平行的不同波速层,且后续方法仅适用于层状波速区划分;本发明可根据矿井地质及开采技术条件划分为任意形状的波速区,且后续方法适用该划分,与现有技术相比传感器布置方式不同;现有技术要求必须在第一和最后一个波速层布置传感器;而本发明对传感器布置位置无特殊要求,只需要在目标矿井布置微震监测系统,在井下和地面布置拾振传感器,其中地面传感器编号T1、T2、T3…Tn,井下传感器编号S1、S2、S3…Sn,所述地面传感器以及所述井下传感器总和大于等于四;并测量出各传感器坐标;与现有技术相比,本发明的定位方式不同;现有技术中对比文件通过假设震源位于某一波速层中,通过求解目标函数最小值确定震源位置;而本发明首先采用实体煤波速计算震源位置,而后根据震源和传感器之间所穿过波速区位置不断重新计算传感器波速值,而后不断迭代直至与上次震源位置误差小于1m为止,即为真实震源位置,使得最终得到的真实震源位置更加精确。In the present invention, the velocity zone division is different from that in the prior art; in the prior art, the rock mass is divided into different velocity layers parallel to each other, and the subsequent method is only applicable to the layered velocity zone division; the present invention can be divided into velocity zones of any shape according to the mine geology and mining technical conditions, and the subsequent method is applicable to the division, and the sensor arrangement method is different from that in the prior art; the prior art requires that sensors must be arranged in the first and last velocity layers; while the present invention has no special requirements for the sensor arrangement position, and only needs to arrange a microseismic monitoring system in the target mine, and arrange vibration pickup sensors underground and on the ground, wherein the ground sensors are numbered T 1 , T 2 , T 3 ... T n , and the underground sensors are numbered S 1 , S 2 , S 3 ... S n , the sum of the ground sensors and the downhole sensors is greater than or equal to four; and the coordinates of each sensor are measured; compared with the prior art, the positioning method of the present invention is different; the comparative document in the prior art assumes that the earthquake source is located in a certain wave velocity layer, and determines the earthquake source position by solving the minimum value of the objective function; and the present invention first uses the wave velocity of the physical coal to calculate the earthquake source position, and then continuously recalculates the sensor wave velocity value according to the position of the wave velocity zone passed between the earthquake source and the sensor, and then continuously iterates until the error with the last earthquake source position is less than 1m, that is, the real earthquake source position, so that the final real earthquake source position is more accurate.
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以作出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。本发明中未具体描述和解释说明的结构、装置以及操作方法,如无特别说明和限定,均按照本领域的常规手段进行实施。The above is only a preferred embodiment of the present invention. It should be noted that, for those skilled in the art, several improvements and modifications can be made without departing from the principles of the present invention, and these improvements and modifications should also be considered as the protection scope of the present invention. The structures, devices and operating methods not specifically described and explained in the present invention shall be implemented according to the conventional means in the art unless otherwise specified and limited.
Claims (8)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410850990.1A CN118688716A (en) | 2024-06-28 | 2024-06-28 | A method for improving positioning accuracy of coal mine microseismic monitoring system based on wave velocity zone division |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410850990.1A CN118688716A (en) | 2024-06-28 | 2024-06-28 | A method for improving positioning accuracy of coal mine microseismic monitoring system based on wave velocity zone division |
Publications (1)
Publication Number | Publication Date |
---|---|
CN118688716A true CN118688716A (en) | 2024-09-24 |
Family
ID=92767755
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410850990.1A Pending CN118688716A (en) | 2024-06-28 | 2024-06-28 | A method for improving positioning accuracy of coal mine microseismic monitoring system based on wave velocity zone division |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN118688716A (en) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105022031A (en) * | 2015-07-03 | 2015-11-04 | 四川大学 | Layered speed positioning method for regional rock microseismic source |
CN107884822A (en) * | 2017-11-13 | 2018-04-06 | 北京矿冶研究总院 | Method for improving positioning precision of mining micro-seismic source |
CN115356767A (en) * | 2022-07-21 | 2022-11-18 | 北京科技大学 | Method and device for improving microseismic positioning and wave velocity inversion accuracy based on region segmentation |
-
2024
- 2024-06-28 CN CN202410850990.1A patent/CN118688716A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105022031A (en) * | 2015-07-03 | 2015-11-04 | 四川大学 | Layered speed positioning method for regional rock microseismic source |
CN107884822A (en) * | 2017-11-13 | 2018-04-06 | 北京矿冶研究总院 | Method for improving positioning precision of mining micro-seismic source |
CN115356767A (en) * | 2022-07-21 | 2022-11-18 | 北京科技大学 | Method and device for improving microseismic positioning and wave velocity inversion accuracy based on region segmentation |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Konicek et al. | Stress changes and seismicity monitoring of hard coal longwall mining in high rockburst risk areas | |
Cheng et al. | Numerical modelling of strata movement at footwall induced by underground mining | |
Li et al. | Deformation and failure analyses of large underground caverns during construction of the Houziyan Hydropower Station, Southwest China | |
Pan et al. | Effect of premining on hard roof distress behavior: a case study | |
CN112302722B (en) | Coal mine roadway multidirectional stress and deformation wireless monitoring and early warning method and system | |
Fan et al. | Advanced stability analysis of the tunnels in jointed rock mass based on TSP and DEM | |
Deng et al. | Influence of blasting vibrations generated by tunnel construction on an existing road | |
CN112285802B (en) | Tunnel earthquake and transient electromagnetic combined detection method | |
Zhu et al. | Overburden movement characteristics of top-coal caving mining in multi-seam areas | |
Meng et al. | In situ investigation and numerical simulation of the failure depth of an inclined coal seam floor: a case study | |
Xiao et al. | Microseismic monitoring and deformation early warning of the underground caverns of Lianghekou hydropower station, Southwest China | |
Ma et al. | An early warning method for water inrush in Dongjiahe coal mine based on microseismic moment tensor | |
Li et al. | Cause Analysis of Soft and Hard Rock Tunnel Collapse and Information Management. | |
Gai et al. | Microseismic response difference and failure analysis of roof and floor strata under dynamic load impact | |
Liu et al. | Nondestructive testing on cumulative damage of watery fractured rock mass under multiple cycle blasting | |
CN110761841B (en) | Mine group working face mining mutual interference range calculation method based on microseismic monitoring | |
CN118688716A (en) | A method for improving positioning accuracy of coal mine microseismic monitoring system based on wave velocity zone division | |
Liu et al. | The spatio-temporal evolution law of overlying rock fractures in an experimental working face (N00 mining method) based on microseismic monitoring technology | |
Li et al. | Numerical Simulation of Overburden Deformation Mechanism and Surface Settlement Characteristics Induced by Underground Coal Mining: A Case Study | |
Wang et al. | Blasting Vibration Law and Prediction in the Near‐Field of Tunnel | |
Guan et al. | Vibration characteristics, attenuation law and prediction method in the near field of tunnel blasting | |
Fan et al. | Application of seismic channel wave technology on small structure exploration in coal mine | |
Wang et al. | Failure height and fracture evolution pattern of overburden rock in fully mechanized cave mining | |
Luo et al. | Numerical analysis about blasting effect of digital detonator and millisecond detonator in upper and lower cross tunnels | |
Fu et al. | Cross-term suppression method for time-frequency spectrum of engineering blasting signals |
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 |