CN114879256B - Rock burst monitoring method and device, medium and equipment - Google Patents
Rock burst monitoring method and device, medium and equipment Download PDFInfo
- Publication number
- CN114879256B CN114879256B CN202210315730.5A CN202210315730A CN114879256B CN 114879256 B CN114879256 B CN 114879256B CN 202210315730 A CN202210315730 A CN 202210315730A CN 114879256 B CN114879256 B CN 114879256B
- Authority
- CN
- China
- Prior art keywords
- vibration
- target area
- computer tomography
- source
- shock 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 239000011435 rock Substances 0.000 title claims abstract description 50
- 238000012544 monitoring process Methods 0.000 title claims abstract description 28
- 230000035939 shock Effects 0.000 claims abstract description 142
- 238000002591 computed tomography Methods 0.000 claims abstract description 133
- 230000004927 fusion Effects 0.000 claims abstract description 56
- 238000009826 distribution Methods 0.000 claims abstract description 35
- 239000003245 coal Substances 0.000 claims description 26
- 238000005422 blasting Methods 0.000 claims description 25
- 238000004891 communication Methods 0.000 claims description 22
- 238000012545 processing Methods 0.000 claims description 8
- 238000003860 storage Methods 0.000 claims description 7
- 230000001133 acceleration Effects 0.000 claims description 6
- 238000012806 monitoring device Methods 0.000 claims description 4
- 238000001514 detection method Methods 0.000 abstract description 9
- 238000005065 mining Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000005284 excitation Effects 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 230000010365 information processing Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000000149 penetrating effect Effects 0.000 description 2
- 238000004590 computer program Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 238000011022 operating instruction Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000003068 static 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
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
-
- 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/34—Displaying seismic recordings or visualisation of seismic data or attributes
- G01V1/345—Visualisation of seismic data or attributes, e.g. in 3D cubes
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The application discloses a rock burst monitoring method, a device, a medium and equipment, relates to the technical field of rock burst detection, and mainly aims to solve the technical problems that the accuracy of an inversion result of a shock wave computer tomography obtained based on a single seismic source is low, the real-time performance is poor, and therefore rock burst cannot be effectively monitored and early-warned. Comprising the following steps: respectively acquiring vibration wave propagation ray distribution data generated by various vibration sources in a target area, and respectively determining vibration wave computer tomography inversion results of the various vibration sources in the target area; based on the number of vibration wave propagation rays passing through the target area, acquiring weight ratios of vibration wave computer tomography inversion results of the vibration sources in various vibration source fusion vibration wave computer tomography inversion results; and determining a plurality of seismic source fusion shock wave computer tomography inversion results of the target area based on the weight ratio. The method is mainly used for detecting rock burst.
Description
Technical Field
The application relates to the technical field of rock burst monitoring, in particular to a rock burst monitoring method and device, a medium and equipment.
Background
With the increasing contradiction between the demand of coal resources and the gradual exhaustion of shallow coal resources, a plurality of coal mines will fully enter the deep mining stage. In the environment of deep mining 'three high and one disturbance', rock burst disasters become one of the important difficulties faced by mine safety production. The rock burst is monitored and early-warned, so that disaster loss caused by the rock burst can be effectively prevented. Traditional rock burst monitoring and early warning mainly focuses on static evaluation on impact dangers, dynamic load monitoring based on real-time monitoring data, embedding stress meters in rock mass and other modes, and cannot effectively monitor rock burst.
At present, a CT detection method (shock wave computer tomography) is widely applied to rock burst monitoring, namely a detection imaging technology for measuring the internal density and structure of an opaque object, mainly utilizing the positive correlation relation between the propagation wave speed of shock waves and stress, and calculating the stress distribution by detecting the propagation speed of the shock waves in a coal rock mass. However, on one hand, the inversion of the passive shock wave computer tomography based on the microseismic event has lower inversion precision due to the influence of the positioning precision of the microseismic event; on the other hand, the active shock wave computer tomography detection based on blasting is high in accuracy, but is influenced by underground blasting operation, and can be performed only in individual sections and time periods, so that the real-time performance is poor.
Disclosure of Invention
In view of this, the present application provides a rock burst monitoring method and device, and mainly aims to solve the technical problem that the accuracy of the obtained shock wave computer tomography result based on a single seismic source is low and the real-time performance is poor, so that the rock burst cannot be effectively monitored and pre-warned.
According to one aspect of the present application, there is provided a rock burst monitoring method, comprising:
respectively acquiring vibration wave propagation ray distribution data generated by various vibration sources in a target area, and respectively determining vibration wave computer tomography inversion results of the various vibration sources in the target area;
based on the number of vibration wave propagation rays passing through the target area, acquiring weight ratios of vibration wave computer tomography inversion results of the vibration sources in various vibration source fusion vibration wave computer tomography inversion results;
and determining a plurality of seismic source fusion shock wave computer tomography inversion results of the target area based on the weight ratio.
Preferably, the determining the inversion result of the multiple seismic source fusion shock wave computer tomography of the target area based on the weight ratio specifically includes:
and based on the weight ratio, performing superposition processing on the shock wave computer tomography inversion results of various shock sources in the target area to obtain various shock source fusion shock wave computer tomography inversion results of the target area.
Preferably, the weight ratio of the shock wave computer tomography inversion result of the seismic source in the shock wave computer tomography inversion result of the seismic source fusion is obtained, which specifically includes:
in the target area, if the number of the vibration wave propagation rays generated by the artificial vibration source is greater than or equal to a preset vibration wave propagation ray number threshold, weight ratio is distributed according to the proportion that the number of the vibration wave propagation rays generated by the vibration source is the sum of the vibration wave propagation rays generated by various vibration sources except the artificial vibration source;
and in the target area, if the number of the vibration wave propagation rays generated by the artificial vibration source is smaller than a preset vibration wave propagation ray number threshold, weight ratio is distributed according to the proportion of the number of the vibration wave propagation rays generated by various vibration sources to the sum of the number of the vibration wave propagation rays generated by various vibration sources.
Preferably, the method for determining the shock wave computed tomography inversion results of the various seismic sources in the target area specifically includes:
generating vibration wave propagation ray distribution data generated in the target area by the vibration source based on the position coordinate data of the vibration source, the vibration starting time and the vibration wave receiving time;
and acquiring a shock wave computer tomography inversion result of the seismic source in the target area based on the shock wave propagation ray distribution data.
Preferably, before generating the vibration wave propagation ray distribution data generated in the target area by the vibration source based on the position coordinate data, the starting time and the vibration wave receiving time of the vibration source, the method further includes:
establishing a data connection with the vibration sensor;
and based on the data connection, indicating the vibration sensor to acquire position coordinate data, the starting moment and the vibration wave receiving moment of the vibration source so as to generate vibration wave propagation ray distribution data generated by the vibration source in the target area.
Preferably, the method further comprises:
the vibration sensor comprises a speed sensor and/or an acceleration sensor, so that position coordinate data, a starting moment and a vibration wave receiving moment of the vibration source are acquired based on the vibration sensor.
Preferably, the seismic source comprises: at least two of a coal cutter focus, a blasting focus, a manual focus and a microseismic event focus.
According to another aspect of the present application, there is provided a rock burst monitoring device comprising:
the first acquisition module is used for respectively acquiring vibration wave propagation ray distribution data generated by various vibration sources in a target area and respectively determining vibration wave computer tomography inversion results of the various vibration sources in the target area;
the second acquisition module is used for acquiring the weight ratio of the shock wave computer tomography inversion result of the seismic source in a plurality of seismic source fusion shock wave computer tomography inversion results based on the number of the shock wave propagation rays passing through the target area;
and the determining module is used for determining a plurality of seismic source fusion shock wave computer tomography inversion results of the target area based on the weight ratio.
Preferably, the determining module is specifically configured to:
and based on the weight ratio, performing superposition processing on the shock wave computer tomography inversion results of various shock sources in the target area to obtain various shock source fusion shock wave computer tomography inversion results of the target area.
Preferably, the second obtaining module is specifically configured to:
in the target area, if the number of the vibration wave propagation rays generated by the artificial vibration source is greater than or equal to a preset vibration wave propagation ray number threshold, weight ratio is distributed according to the proportion that the number of the vibration wave propagation rays generated by the vibration source is the sum of the vibration wave propagation rays generated by various vibration sources except the artificial vibration source;
and in the target area, if the number of the vibration wave propagation rays generated by the artificial vibration source is smaller than a preset vibration wave propagation ray number threshold, weight ratio is distributed according to the proportion of the number of the vibration wave propagation rays generated by various vibration sources to the sum of the number of the vibration wave propagation rays generated by various vibration sources.
Preferably, the first obtaining module specifically includes:
the generation unit is used for generating vibration wave propagation ray distribution data generated by the vibration source in the target area based on the position coordinate data of the vibration source, the vibration starting time and the vibration wave receiving time;
the first acquisition unit is used for acquiring the shock wave computer tomography inversion result of the seismic source in the target area based on the shock wave propagation ray distribution data.
Preferably, before the generating unit, the module further includes:
the establishing unit is used for establishing data connection with the vibration sensor;
and the second acquisition unit is used for indicating the vibration sensor to acquire position coordinate data, starting time and vibration wave receiving time of the vibration source based on the data connection so as to generate vibration wave propagation ray distribution data generated by the vibration source in the target area.
Preferably, the apparatus further comprises:
the vibration sensor comprises a speed sensor and/or an acceleration sensor, so that position coordinate data, a starting moment and a vibration wave receiving moment of the vibration source are acquired based on the vibration sensor.
Preferably, the seismic source comprises: at least two of a coal cutter focus, a blasting focus, a manual focus and a microseismic event focus.
According to yet another aspect of the present application, there is provided a storage medium having stored therein at least one executable instruction for causing a processor to perform operations corresponding to the rock burst monitoring method described above.
According to still another aspect of the present application, there is provided a terminal including: the device comprises a processor, a memory, a communication interface and a communication bus, wherein the processor, the memory and the communication interface complete communication with each other through the communication bus;
the memory is used for storing at least one executable instruction, and the executable instruction enables the processor to execute the operation corresponding to the rock burst monitoring method.
By means of the technical scheme, the technical scheme provided by the embodiment of the application has at least the following advantages:
the application provides a rock burst monitoring method and device, firstly, vibration wave propagation ray distribution data generated by various vibration sources in a target area are respectively obtained, and vibration wave computer tomography inversion results of the various vibration sources in the target area are respectively determined; secondly, based on the number of vibration wave propagation rays passing through the target area, obtaining the weight ratio of the vibration wave computer tomography inversion result of the vibration source in the multiple vibration source fusion vibration wave computer tomography inversion results; and finally, determining a plurality of seismic source fusion shock wave computer tomography inversion results of the target area based on the weight ratio. Compared with the prior art, the method and the device for determining the weight ratio of each seismic source in the multi-seismic-source fusion shock wave computer tomography inversion result are determined by obtaining the quantity of the shock wave propagation rays generated by the plurality of seismic sources and passing through the target area, the multi-seismic-source fusion shock wave computer tomography inversion result in the current target area is further determined according to the weight ratio, the accuracy of the shock wave computer tomography inversion result is effectively improved, meanwhile, the problem of poor real-time performance caused by a single seismic source is solved, and therefore effective detection and early warning of rock burst are achieved.
The foregoing description is only an overview of the technical solutions of the present application, and may be implemented according to the content of the specification in order to make the technical means of the present application more clearly understood, and in order to make the above-mentioned and other objects, features and advantages of the present application more clearly understood, the following detailed description of the present application will be given.
Drawings
Various other advantages and benefits will become apparent to those of ordinary skill in the art upon reading the following detailed description of the preferred embodiments. The drawings are only for purposes of illustrating the preferred embodiments and are not to be construed as limiting the application. Also, like reference numerals are used to designate like parts throughout the figures. In the drawings:
FIG. 1 shows a flow chart of a rock burst monitoring method provided in an embodiment of the present application;
FIG. 2 shows a shock wave computed tomography schematic of a shearer source provided in an embodiment of the present application;
FIG. 3 shows a shock wave computed tomography schematic of an artificial seismic source provided in an embodiment of the present application;
FIG. 4 shows a shock wave computed tomography schematic of a blasted seismic source provided in an embodiment of the present application;
FIG. 5 shows a shock wave computed tomography schematic of a microseismic event source provided by an embodiment of the present application;
FIG. 6 shows a block diagram of an apparatus for rock burst monitoring provided in an embodiment of the present application;
fig. 7 shows a schematic structural diagram of a terminal according to an embodiment of the present application.
Detailed Description
Exemplary embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art.
Meanwhile, it should be understood that the sizes of the respective parts shown in the drawings are not drawn in actual scale for convenience of description.
The following description of at least one exemplary embodiment is merely exemplary in nature and is in no way intended to limit the application, its application, or uses.
Techniques, methods, and apparatus known to one of ordinary skill in the relevant art may not be discussed in detail, but are intended to be part of the specification where appropriate.
It should be noted that: like reference numerals and letters denote like items in the following figures, and thus once an item is defined in one figure, no further discussion thereof is necessary in subsequent figures.
Embodiments of the present application may be applied to computer systems/servers that are operational with numerous other general purpose or special purpose computing system environments or configurations. Examples of well known computing systems, environments, and/or configurations that may be suitable for use with the computer system/server include, but are not limited to: personal computer systems, server computer systems, thin clients, thick clients, hand-held or laptop devices, microprocessor-based systems, set-top boxes, programmable consumer electronics, network personal computers, small computer systems, mainframe computer systems, and distributed cloud computing technology environments that include any of the foregoing, and the like.
A computer system/server may be described in the general context of computer-system-executable instructions, such as program modules, being executed by a computer system. Generally, program modules may include routines, programs, objects, components, logic, data structures, etc., that perform particular tasks or implement particular abstract data types. The computer system/server may be implemented in a distributed cloud computing environment in which tasks are performed by remote processing devices that are linked through a communications network. In a distributed cloud computing environment, program modules may be located in both local and remote computing system storage media including memory storage devices.
The embodiment of the application provides a rock burst monitoring method, as shown in fig. 1, which comprises the following steps:
101. and respectively acquiring vibration wave propagation ray distribution data generated by various vibration sources in the target area, and respectively determining vibration wave computer tomography inversion results of the various vibration sources in the target area.
In this embodiment of the present application, the current execution end may be a rock burst disaster detection system, configured to monitor rock burst during deep mining, so as to predict occurrence of disasters, thereby ensuring safe performance of mining tasks. The vibration wave computer tomography is an imaging technology for measuring the internal density and structure of an opaque object, and mainly utilizes the positive correlation relation between the propagation wave speed of the vibration wave and the stress to calculate the stress distribution by detecting the propagation speed of the vibration wave in a coal rock mass; the inversion results are monitoring results obtained based on such imaging techniques. The prior art generally analyzes the propagation rays of shock waves generated by a certain type of seismic source by using a shock wave computer tomography method, for example, analyzes the propagation rays of shock waves generated by a microseismic event or analyzes the propagation rays of shock waves generated by artificial blasting. However, when analyzing the vibration wave propagation rays generated by the microseismic event, the problem of lower accuracy of the inversion result can be caused due to the influence of the positioning accuracy of the microseismic event; the method is characterized in that the method analyzes vibration wave propagation rays generated by artificial blasting, is influenced by underground blasting operation, can be performed only in individual sections and time periods, and causes the problem of poor real-time inversion results. Based on the above, the embodiments of the present application overcome the above-mentioned drawbacks in the prior art by fusing the shock wave propagation ray distribution data generated by multiple seismic sources in the target area. Firstly, vibration wave propagation ray distribution data generated by various vibration sources in a target area are obtained, and vibration wave computer tomography inversion results of the various vibration sources in the target area are respectively determined.
It should be noted that the number of the multiple seismic sources may be two or more, for example, any two or three of the blasting sources, artificial sources, and microseismic events.
102. Based on the number of shock wave propagation rays passing through the target area, the weight ratio of the shock wave computer tomography inversion result of the seismic source in the multiple seismic source fusion shock wave computer tomography inversion results is obtained.
In the embodiment of the application, it can be understood that, due to the adoption of the rock burst monitoring method of fusion of multiple seismic sources, the shock wave computer tomography inversion results of all the seismic sources need to be fused, so as to obtain the shock wave computer tomography inversion results after fusion of all the seismic sources. Specifically, the shock wave computer tomography inversion results of all the seismic sources are fused based on the weight ratio of the shock wave computer tomography inversion results of all the seismic sources in the shock wave computer tomography inversion results of various seismic source fusion.
In the target area, the weight ratio of the shock wave computer tomography inversion result of each seismic source in the shock wave computer tomography inversion result of multiple seismic source fusion has an association relation with the number of rays of the shock wave propagation rays generated by each seismic source passing through the target area, so that the weight ratio of the shock wave computer tomography inversion result of each seismic source in the shock wave computer tomography inversion result of multiple seismic source fusion can be obtained based on the association relation and the number of rays of the shock wave propagation rays generated by each seismic source passing through the target area. The target area is used for representing a current area monitored by a current executing end.
103. And determining the inversion result of the multiple seismic source fusion shock wave computer tomography of the target area based on the weight ratio.
In the embodiment of the application, based on the weight ratio of the shock wave computer tomography inversion result of each shock source obtained in step 102 in the shock wave computer tomography inversion result of multiple shock source fusion, the shock wave computer tomography inversion results of each shock source are fused to obtain multiple shock source fusion shock wave computer tomography inversion results in a target area, so that the monitoring of rock burst of the target area is realized.
Compared with the prior art, the method and the device for determining the weight ratio of each seismic source in the multi-seismic-source fusion shock wave computer tomography inversion result are determined by obtaining the quantity of the shock wave propagation rays generated by the plurality of seismic sources and passing through the target area, the multi-seismic-source fusion shock wave computer tomography inversion result in the current target area is further determined according to the weight ratio, the accuracy of the shock wave computer tomography inversion result is effectively improved, meanwhile, the problem of poor real-time performance caused by a single seismic source is solved, and therefore effective detection and early warning of rock burst are achieved.
For further explanation and limitation, in the embodiments of the present application, determining a plurality of seismic source fusion shock wave computed tomography inversion results based on weight ratios specifically includes: and (3) based on the weight ratio, carrying out superposition processing on the shock wave computer tomography inversion results of various seismic sources in the target area to obtain the multi-seismic source fusion shock wave computer tomography inversion results of the target area.
Specifically, because the weight ratio of the shock wave computer tomography inversion results of each seismic source in the shock wave computer tomography inversion results of multiple seismic source fusion shock waves is used for representing the ray density of the shock wave propagation rays generated by each seismic source and passing through the target area, the shock wave computer tomography inversion results of the various seismic sources in the target area are subjected to superposition processing based on the weight ratio, and the accurate multiple shock wave computer tomography inversion results of the multiple seismic source fusion shock waves in the target area can be obtained, so that effective monitoring and early warning of rock burst in the target area can be realized. The superposition formula is as follows:
V=a*V 1 +b*V 2 +c*V 3 +d*V 4
wherein V represents inversion results of multiple seismic source fusion shock wave computer tomography; v (V) 1 A represents the weight ratio of the shock wave computer tomography inversion result of the coal cutter seismic source in the shock wave computer tomography inversion result of the various seismic source fusion shock wave computer tomography inversion results; v (V) 2 B represents the weight ratio of the shock wave computer tomography inversion result of the artificial seismic source in the shock wave computer tomography inversion result of the various seismic source fusion shock waves; v (V) 3 C represents the weight ratio of the shock wave computer tomography inversion result of the blasting seismic source in the shock wave computer tomography inversion result of the various seismic source fusion shock waves; v (V) 4 And d represents the weight ratio of the shock wave computer tomography inversion result of the micro-seismic event source in the multi-source fusion shock wave computer tomography inversion result.
For further explanation and limitation, in the embodiment of the present application, the weight ratio of the shock wave computed tomography inversion result of the seismic source to the shock wave computed tomography inversion result of the plurality of seismic source fusion shock waves is specifically: in the target area, if the number of the vibration wave propagation rays generated by the artificial vibration source is larger than or equal to a preset vibration wave propagation ray number threshold value, weight ratios are distributed according to the proportion that the number of the vibration wave propagation rays generated by the vibration source accounts for the sum of the numbers of the vibration wave propagation rays generated by various vibration sources except the artificial vibration source; and in the target area, if the number of the vibration wave propagation rays generated by the artificial vibration source is smaller than a preset vibration wave propagation ray number threshold value, weight ratios are distributed according to the proportion that the number of the vibration wave propagation rays generated by various vibration sources accounts for the sum of the number of the vibration wave propagation rays generated by various vibration sources.
Specifically, firstly, a target area is divided into grids with preset sizes (such as grids with 20m x 20 m), the number of vibration wave propagation rays generated by a current seismic source passing through the target area is obtained, and the distribution of weight ratios can be divided into two cases based on the difference of the number of vibration wave propagation rays generated by an artificial seismic source.
(one) if N2 i When the number of the propagation rays of the shock waves is more than or equal to a preset threshold value (e.g. 10);
b=M
(II) if N2 i <Presetting a vibration wave propagation ray quantity threshold (e.g. 10);
wherein N1 i -the number of rays in the ith grid during computer tomography of the seismic source shock wave of the shearer;
N2 i -the number of rays in the ith grid during computer tomography of the artificial source shock wave;
N3 i -the number of rays in the ith grid during computer tomography of the blasts source shock wave;
N4 i -the number of rays in the ith grid during computer tomography of the source shock wave of the microseismic event;
m-weight ratio of artificial seismic source vibration wave computer tomography inversion result.
For further explanation and limitation, in the embodiments of the present application, the shock wave computed tomography inversion results of various seismic sources in the target area are determined respectively, and specifically include: generating vibration wave propagation ray distribution data generated in a target area by a vibration source based on position coordinate data of the vibration source, a vibration starting time and a vibration wave receiving time; and acquiring a shock wave computer tomography inversion result of the seismic source in the target area based on the shock wave propagation ray distribution data.
For example, taking four seismic sources as examples, the four seismic sources are respectively a coal cutter seismic source, a manual seismic source, a blasting seismic source and a micro-seismic event seismic source. As shown in fig. 2, the vibration source coordinates and the vibration time of the mechanical vibration of the coal cutter are calculated in real time based on the starting position coordinates of the coal cutter and the travel record of the coal cutter, so as to obtain the vibration wave propagation rays from the coal cutter to the roadway vibration sensors on two sides of the working surface, and the vibration wave computer tomography inversion result 1 of the vibration source of the coal cutter is obtained based on the vibration wave propagation rays. As shown in fig. 3, the artificial seismic sources installed in the lanes on two sides of the working surface are used for active excitation (on-line or portable), so as to obtain the vibration wave propagation rays from the artificial seismic sources to the vibration sensors on the opposite sides of the lane, and the vibration wave computed tomography inversion result 2 of the artificial seismic sources is obtained based on the vibration wave propagation rays. In the area inconvenient to use artificial seismic source, the active excitation is carried out by using a blasting mode, the blasting can be carried out by using broken roof blasting, bottom plate blasting, coal seam blasting and the like which are carried out frequently by using a coal mine, and the blasting excitation can also be carried out by special vibration wave computer tomography, as shown in fig. 4, the blasting point and the vibration sensors for receiving vibration signals are respectively positioned in the roadways on two sides of the working face, the blasting charge of the special blasting for the vibration wave computer tomography is not less than 200g of coal mine allowed explosive, the blasting hole is constructed by dividing the waist line of the upper part of the roadway by the vertical coal wall, and the vibration wave propagation rays from the blasting seismic source to the vibration sensors on the roadway on two sides of the working face are obtained, and the vibration wave computer tomography inversion result 3 of the blasting seismic source is obtained based on the vibration wave propagation rays. As shown in fig. 5, using a microseismic event frequently occurring in a coal mine as a seismic source, obtaining a vibration wave propagation ray from the microseismic event to roadway vibration sensors on two sides of a working surface, and obtaining a vibration wave computer tomography inversion result 4 of the microseismic event based on the vibration wave propagation ray.
Further, in the embodiment of the present application, before generating the vibration wave propagation ray distribution data generated in the target area by the vibration source based on the position coordinate data, the starting time and the vibration wave receiving time of the vibration source, the method of the embodiment further includes: establishing a data connection with the vibration sensor; based on the data connection, the vibration sensor is instructed to acquire position coordinate data, a starting moment and a vibration wave receiving moment of the vibration source so as to generate vibration wave propagation ray distribution data generated by the vibration source in the target area.
Specifically, data connection between the current execution end and each seismic source is established, and a vibration sensor is indicated to capture the starting moment and the receiving moment of vibration waves of the seismic source in real time through the data connection. And meanwhile, the position coordinate data of the vibration sensor and each seismic source are acquired through a computer program or a manual recording mode. Further, shock wave propagation ray distribution data generated in the target area by each seismic source is generated based on the acquired data.
Optionally, in an embodiment of the present application, the embodiment method further includes: the vibration sensor comprises a speed sensor and/or an acceleration sensor to acquire position coordinate data of a seismic source, a starting time and a vibration wave receiving time based on the vibration sensor.
Specifically, the sensor may be, but is not limited to, a speed sensor, an acceleration sensor, or the like. It may be previously arranged at a work site according to a preset rule. Illustratively, vibration sensors are arranged at the same time in the roadways on two sides of the stope face; 4-8 sensors are arranged on each roadway, the distance between the sensors is 20-50 m, and the sensors are arranged on the side parts of the roadway or the anchor rod end parts of the top plate; the length of the anchor rod exposed out of the coal rock mass is generally 20mm; the length of the anchor rod penetrating into the coal rock body is not less than 1.5m, and the anchor rod is anchored in the whole course.
Optionally, in an embodiment of the present application, the seismic source includes: at least two of a coal cutter focus, a blasting focus, a manual focus and a microseismic event focus.
When the artificial seismic source is adopted, real-time on-line or portable artificial seismic sources can be installed in the roadways on two sides of the stoping working face. When the real-time on-line state is adopted, 4-8 artificial seismic sources are installed on each side of the roadway, and the distance between the artificial seismic sources is 20-50 m; when used for portability, a single artificial seismic source is used to fire every 10m in the roadway. The artificial seismic source is installed at the end part of an anchor rod at the side part or the bottom plate of a roadway when in use, the length of the anchor rod exposed out of a coal rock body is generally 20mm, the length of the anchor rod penetrating into the coal rock body is not less than 1.5m, and the anchor rod is anchored in the whole course.
The application provides a rock burst monitoring method, which comprises the steps of firstly, respectively acquiring vibration wave propagation ray distribution data generated by various vibration sources in a target area, and respectively determining vibration wave computer tomography inversion results of the various vibration sources in the target area; secondly, based on the number of vibration wave propagation rays passing through the target area, obtaining the weight ratio of the vibration wave computer tomography inversion result of the vibration source in the multiple vibration source fusion vibration wave computer tomography inversion results; and finally, determining a plurality of seismic source fusion shock wave computer tomography inversion results of the target area based on the weight ratio. Compared with the prior art, the method and the device for determining the weight ratio of each seismic source in the multi-seismic-source fusion shock wave computer tomography inversion result are determined by obtaining the quantity of the shock wave propagation rays generated by the plurality of seismic sources and passing through the target area, the multi-seismic-source fusion shock wave computer tomography inversion result in the current target area is further determined according to the weight ratio, the accuracy of the shock wave computer tomography inversion result is effectively improved, meanwhile, the problem of poor real-time performance caused by a single seismic source is solved, and therefore effective detection and early warning of rock burst are achieved.
Further, as an implementation of the method shown in fig. 1, an embodiment of the present application provides a rock burst monitoring device, as shown in fig. 6, where the device includes:
the first acquisition module 21, the second acquisition module 22, the determination module 23.
A first obtaining module 21, configured to obtain vibration wave propagation ray distribution data generated by various vibration sources in a target area, and determine vibration wave computed tomography inversion results of the various vibration sources in the target area;
a second obtaining module 22, configured to obtain a weight ratio of the shock wave computed tomography inversion result of the seismic source in a plurality of seismic source fusion shock wave computed tomography inversion results based on a number of shock wave propagation rays passing through the target area;
a determining module 23, configured to determine a plurality of seismic source fusion shock wave computed tomography inversion results of the target area based on the weight ratio.
In a specific application scenario, the determining module is specifically configured to:
and based on the weight ratio, performing superposition processing on the shock wave computer tomography inversion results of various shock sources in the target area to obtain various shock source fusion shock wave computer tomography inversion results of the target area.
In a specific application scenario, the second obtaining module is specifically configured to:
in the target area, if the number of the vibration wave propagation rays generated by the artificial vibration source is greater than or equal to a preset vibration wave propagation ray number threshold, weight ratio is distributed according to the proportion that the number of the vibration wave propagation rays generated by the vibration source is the sum of the vibration wave propagation rays generated by various vibration sources except the artificial vibration source;
and in the target area, if the number of the vibration wave propagation rays generated by the artificial vibration source is smaller than a preset vibration wave propagation ray number threshold, weight ratio is distributed according to the proportion of the number of the vibration wave propagation rays generated by various vibration sources to the sum of the number of the vibration wave propagation rays generated by various vibration sources.
In a specific application scenario, the first obtaining module specifically includes:
the generation unit is used for generating vibration wave propagation ray distribution data generated by the vibration source in the target area based on the position coordinate data of the vibration source, the vibration starting time and the vibration wave receiving time;
the first acquisition unit is used for acquiring the shock wave computer tomography inversion result of the seismic source in the target area based on the shock wave propagation ray distribution data.
In a specific application scenario, before the generating unit, the module further includes:
the establishing unit is used for establishing data connection with the vibration sensor;
and the second acquisition unit is used for indicating the vibration sensor to acquire position coordinate data, starting time and vibration wave receiving time of the vibration source based on the data connection so as to generate vibration wave propagation ray distribution data generated by the vibration source in the target area.
In a specific application scenario, the apparatus further includes:
the vibration sensor comprises a speed sensor and/or an acceleration sensor, so that position coordinate data, a starting moment and a vibration wave receiving moment of the vibration source are acquired based on the vibration sensor.
In a specific application scenario, the seismic source includes: at least two of a coal cutter focus, a blasting focus, a manual focus and a microseismic event focus.
The application provides a rock burst monitoring device, which comprises the steps of firstly, respectively acquiring vibration wave propagation ray distribution data generated by various vibration sources in a target area, and respectively determining vibration wave computer tomography inversion results of the various vibration sources in the target area; secondly, based on the number of vibration wave propagation rays passing through the target area, obtaining the weight ratio of the vibration wave computer tomography inversion result of the vibration source in the multiple vibration source fusion vibration wave computer tomography inversion results; and finally, determining a plurality of seismic source fusion shock wave computer tomography inversion results of the target area based on the weight ratio. Compared with the prior art, the method and the device for determining the weight ratio of each seismic source in the multi-seismic-source fusion shock wave computer tomography inversion result are determined by obtaining the quantity of the shock wave propagation rays generated by the plurality of seismic sources and passing through the target area, the multi-seismic-source fusion shock wave computer tomography inversion result in the current target area is further determined according to the weight ratio, the accuracy of the shock wave computer tomography inversion result is effectively improved, meanwhile, the problem of poor real-time performance caused by a single seismic source is solved, and therefore effective detection and early warning of rock burst are achieved.
According to one embodiment of the present application, there is provided a storage medium storing at least one executable instruction for performing the rock burst monitoring method of any of the method embodiments described above.
Based on such understanding, the technical solution of the present application may be embodied in the form of a software product, which may be stored in a non-volatile storage medium (may be a CD-ROM, a U-disk, a mobile hard disk, etc.), and includes several instructions for causing a computer device (may be a personal computer, a server, or a network device, etc.) to perform the methods described in various implementation scenarios of the present application.
Fig. 7 is a schematic structural diagram of a terminal according to an embodiment of the present application, and the specific embodiment of the present application is not limited to a specific implementation of the terminal.
As shown in fig. 7, the terminal may include: a processor (processor) 302, a communication interface (Communications Interface) 304, a memory (memory) 306, and a communication bus 308.
Wherein: processor 302, communication interface 304, and memory 306 perform communication with each other via communication bus 308.
A communication interface 304 for communicating with network elements of other devices, such as clients or other servers.
The processor 302 is configured to execute the program 310, and may specifically perform the relevant steps in the rock burst monitoring method embodiment described above.
In particular, program 310 may include program code including computer-operating instructions.
The processor 302 may be a central processing unit CPU, or a specific integrated circuit ASIC (Application Specific Integrated Circuit), or one or more integrated circuits configured to implement embodiments of the present application. The one or more processors included in the terminal may be the same type of processor, such as one or more CPUs; but may also be different types of processors such as one or more CPUs and one or more ASICs.
Memory 306 for storing programs 310. Memory 306 may comprise high-speed RAM memory or may also include non-volatile memory (non-volatile memory), such as at least one disk memory.
Program 310 may be specifically operable to cause processor 302 to:
respectively acquiring vibration wave propagation ray distribution data generated by various vibration sources in a target area, and respectively determining vibration wave computer tomography inversion results of the various vibration sources in the target area;
based on the number of vibration wave propagation rays passing through the target area, acquiring weight ratios of vibration wave computer tomography inversion results of the vibration sources in various vibration source fusion vibration wave computer tomography inversion results;
and determining a plurality of seismic source fusion shock wave computer tomography inversion results of the target area based on the weight ratio.
The storage medium may also include an operating system, a network communication module. The operating system is a program for managing the physical equipment hardware and software resources based on rock burst monitoring of multiple seismic source fusion, and supports the operation of information processing programs and other software and/or programs. The network communication module is used for realizing communication among all components in the storage medium and communication with other hardware and software in the information processing entity equipment.
In this specification, each embodiment is described in a progressive manner, and each embodiment is mainly described in a different manner from other embodiments, so that the same or similar parts between the embodiments are mutually referred to. For system embodiments, the description is relatively simple as it essentially corresponds to method embodiments, and reference should be made to the description of method embodiments for relevant points.
The methods and systems of the present application may be implemented in a number of ways. For example, the methods and systems of the present application may be implemented by software, hardware, firmware, or any combination of software, hardware, firmware. The above-described sequence of steps for the method is for illustration only, and the steps of the method of the present application are not limited to the sequence specifically described above unless specifically stated otherwise. Furthermore, in some embodiments, the present application may also be implemented as programs recorded in a recording medium, the programs including machine-readable instructions for implementing the methods according to the present application. Thus, the present application also covers a recording medium storing a program for executing the method according to the present application.
It will be appreciated by those skilled in the art that the modules or steps of the application described above may be implemented in a general purpose computing device, they may be centralized on a single computing device, or distributed across a network of computing devices, or they may alternatively be implemented in program code executable by computing devices, such that they may be stored in a memory device for execution by the computing devices and, in some cases, the steps shown or described may be performed in a different order than what is shown or described, or they may be implemented as individual integrated circuit modules, or as individual integrated circuit modules. Thus, the present application is not limited to any specific combination of hardware and software.
The foregoing description is only of the preferred embodiments of the present application and is not intended to limit the same, but rather, various modifications and variations may be made by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principles of the present application should be included in the protection scope of the present application.
Claims (10)
1. A method of rock burst monitoring comprising:
respectively acquiring vibration wave propagation ray distribution data generated by various vibration sources in a target area, and respectively determining vibration wave computer tomography inversion results of the various vibration sources in the target area;
based on the number of vibration wave propagation rays passing through the target area, acquiring weight ratios of vibration wave computer tomography inversion results of the vibration sources in various vibration source fusion vibration wave computer tomography inversion results;
and determining a plurality of seismic source fusion shock wave computer tomography inversion results of the target area based on the weight ratio.
2. The method of claim 1, wherein determining a plurality of source fusion shock wave computed tomography inversion results for the target region based on the weight ratios, comprises:
and based on the weight ratio, performing superposition processing on the shock wave computer tomography inversion results of various shock sources in the target area to obtain various shock source fusion shock wave computer tomography inversion results of the target area.
3. The method of claim 1, wherein the obtaining the weight ratio of the shock wave computed tomography inversion result of the seismic source in the multiple seismic source fusion shock wave computed tomography inversion result specifically comprises:
in the target area, if the number of the vibration wave propagation rays generated by the artificial vibration source is greater than or equal to a preset vibration wave propagation ray number threshold, weight ratio is distributed according to the proportion that the number of the vibration wave propagation rays generated by the vibration source is the sum of the vibration wave propagation rays generated by various vibration sources except the artificial vibration source;
and in the target area, if the number of the vibration wave propagation rays generated by the artificial vibration source is smaller than a preset vibration wave propagation ray number threshold, weight ratio is distributed according to the proportion of the number of the vibration wave propagation rays generated by various vibration sources to the sum of the number of the vibration wave propagation rays generated by various vibration sources.
4. The method of claim 1, wherein the determining shock wave computed tomography inversion results for each seismic source within the target area comprises:
generating vibration wave propagation ray distribution data generated in the target area by the vibration source based on the position coordinate data of the vibration source, the vibration starting time and the vibration wave receiving time;
and acquiring a shock wave computer tomography inversion result of the seismic source in the target area based on the shock wave propagation ray distribution data.
5. The method of claim 4, wherein the generating shock wave propagation ray distribution data generated by the seismic source within the target area based on the location coordinate data, the time of onset, and the time of shock wave reception of the seismic source further comprises:
establishing a data connection with the vibration sensor;
and based on the data connection, indicating the vibration sensor to acquire position coordinate data, the starting moment and the vibration wave receiving moment of the vibration source so as to generate vibration wave propagation ray distribution data generated by the vibration source in the target area.
6. The method according to any of claims 1-5, wherein the method further comprises:
the vibration sensor comprises a speed sensor and/or an acceleration sensor, so that position coordinate data, a starting moment and a vibration wave receiving moment of the vibration source are acquired based on the vibration sensor.
7. The method of claim 1, wherein the seismic source comprises: at least two of a coal cutter focus, a blasting focus, a manual focus and a microseismic event focus.
8. A rock burst monitoring device, comprising:
the first acquisition module is used for respectively acquiring vibration wave propagation ray distribution data generated by various vibration sources in a target area and respectively determining vibration wave computer tomography inversion results of the various vibration sources in the target area;
the second acquisition module is used for acquiring the weight ratio of the shock wave computer tomography inversion result of the seismic source in a plurality of seismic source fusion shock wave computer tomography inversion results based on the number of the shock wave propagation rays passing through the target area;
and the determining module is used for determining a plurality of seismic source fusion shock wave computer tomography inversion results of the target area based on the weight ratio.
9. A storage medium having stored therein at least one executable instruction that causes a processor to perform operations corresponding to the rock burst monitoring method of any one of claims 1-7.
10. An electronic device, comprising: the device comprises a processor, a memory, a communication interface and a communication bus, wherein the processor, the memory and the communication interface complete communication with each other through the communication bus;
the memory is configured to store at least one executable instruction that causes the processor to perform operations corresponding to the rock burst monitoring method of any one of claims 1-7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210315730.5A CN114879256B (en) | 2022-03-29 | 2022-03-29 | Rock burst monitoring method and device, medium and equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210315730.5A CN114879256B (en) | 2022-03-29 | 2022-03-29 | Rock burst monitoring method and device, medium and equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114879256A CN114879256A (en) | 2022-08-09 |
| CN114879256B true CN114879256B (en) | 2024-02-02 |
Family
ID=82670669
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210315730.5A Active CN114879256B (en) | 2022-03-29 | 2022-03-29 | Rock burst monitoring method and device, medium and equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114879256B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118409361B (en) * | 2024-06-27 | 2024-09-17 | 徐州华东机械有限公司 | Method and system for monitoring rock burst of ultrathin coal seam |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102667529A (en) * | 2009-10-27 | 2012-09-12 | 雪佛龙美国公司 | Method and system for seismic imaging and earth modeling using beam tomography |
| CN107817524A (en) * | 2017-09-28 | 2018-03-20 | 安徽理工大学 | The method and apparatus of three-dimensional seismic tomography |
| CN111413735A (en) * | 2020-05-11 | 2020-07-14 | 安徽理工大学 | Coal face rapid earthquake transmission chromatography method capable of simultaneously exciting multiple seismic sources |
| CN112114359A (en) * | 2020-08-13 | 2020-12-22 | 中南大学 | Dangerous area detection method, system and terminal based on active and passive seismic source signals and readable storage medium |
| CN114219211A (en) * | 2021-11-12 | 2022-03-22 | 中国矿业大学 | Rock burst risk assessment method based on theory and data-driven fusion |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130258810A1 (en) * | 2012-03-30 | 2013-10-03 | Wenyi Hu | Method and System for Tomographic Inversion |
-
2022
- 2022-03-29 CN CN202210315730.5A patent/CN114879256B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102667529A (en) * | 2009-10-27 | 2012-09-12 | 雪佛龙美国公司 | Method and system for seismic imaging and earth modeling using beam tomography |
| CN107817524A (en) * | 2017-09-28 | 2018-03-20 | 安徽理工大学 | The method and apparatus of three-dimensional seismic tomography |
| CN111413735A (en) * | 2020-05-11 | 2020-07-14 | 安徽理工大学 | Coal face rapid earthquake transmission chromatography method capable of simultaneously exciting multiple seismic sources |
| CN112114359A (en) * | 2020-08-13 | 2020-12-22 | 中南大学 | Dangerous area detection method, system and terminal based on active and passive seismic source signals and readable storage medium |
| CN114219211A (en) * | 2021-11-12 | 2022-03-22 | 中国矿业大学 | Rock burst risk assessment method based on theory and data-driven fusion |
Non-Patent Citations (2)
| Title |
|---|
| 蔡武 ; 窦林名 ; 李振雷 ; 巩思园 ; 韩荣军 ; 刘军 ; .矿震震动波速度层析成像评估冲击危险的验证.地球物理学报.2016,(第01期),全文. * |
| 贺日政 ; 高锐 ; 郑洪伟 ; 管烨 ; 李秋生 ; 李文辉 ; 熊小松 ; 邓攻 ; .地震波速层析成像方法研究进展.地质学报.2010,(第06期),全文. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114879256A (en) | 2022-08-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Dong et al. | Three-dimensional analytical solution of acoustic emission source location for cuboid monitoring network without pre-measured wave velocity | |
| PL230226B1 (en) | Method and the system for the assessment of the occurrence of a hazard of high-power shocks, generated by the underground exploitation | |
| Feng et al. | Performance and feasibility analysis of two microseismic location methods used in tunnel engineering | |
| CN112114359B (en) | Dangerous area detection method, system and terminal based on active and passive seismic source signals and readable storage medium | |
| CN114879256B (en) | Rock burst monitoring method and device, medium and equipment | |
| Huang et al. | Source identification of microseismic events in underground mines with interferometric imaging and cross wavelet transform | |
| CN113189644B (en) | Microseismic source positioning method and system | |
| CN105607040A (en) | Mining area illegal mining prevention monitoring and positioning method and system | |
| Dong et al. | Dynamic stability analysis of rockmass: a review | |
| Sherki et al. | Design of real time sensor system for detection and processing of seismic waves for earthquake early warning system | |
| Yan et al. | On the proper estimation of the confidence interval for the design formula of blast-induced vibrations with site records | |
| Song et al. | Determination of RVE size based on the 3D fracture persistence | |
| Chen et al. | Research on in situ stress inversion of deep‐buried tunnel based on pressure/tension axis mechanism and geological structure | |
| Huang et al. | Stability analysis of deep-buried hard rock underground laboratories based on stereophotogrammetry and discontinuity identification | |
| Chen et al. | Replacing RQD and discontinuity spacing with the modified blockiness index in the rock mass rating system | |
| Asming et al. | LORS2 software package for monitoring the seismicity of the Khibiny massif: Construction principles and basic algorithms | |
| CN109521221A (en) | A kind of drill bursting construction hard rock tunnel microseism wave velocity of wave real time acquiring method | |
| Tzou et al. | Enhancing the safety management of NATM using the tunnel seismic prediction method: a case study | |
| Spathis | Innovations in blast measurement: Reinventing the past | |
| Vinoth et al. | Status and developments of slope monitoring techniques in opencast mines | |
| Hadjigeorgiou et al. | Drift reinforcement design based on Discrete Fracture Network (DFN) modelling | |
| Gangrade et al. | Investigating seismicity surrounding an excavation boundary in a highly stressed dipping underground limestone mine | |
| Ghaedi Ghalini et al. | Estimation of in-situ block size distribution in jointed rock masses using combined photogrammetry and discrete fracture network | |
| CN103389522A (en) | Method and system for forecasting tunnel geology | |
| CN109100800A (en) | A kind of surrounding rock stability method for early warning and system based on micro-seismic imaging |
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 |