CN114779331A - Mine earthquake risk area prediction method based on accumulated microseismic response - Google Patents

Mine earthquake risk area prediction method based on accumulated microseismic response Download PDF

Info

Publication number
CN114779331A
CN114779331A CN202210450383.7A CN202210450383A CN114779331A CN 114779331 A CN114779331 A CN 114779331A CN 202210450383 A CN202210450383 A CN 202210450383A CN 114779331 A CN114779331 A CN 114779331A
Authority
CN
China
Prior art keywords
roadway
ppv
microseismic
risk area
predicting
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
Application number
CN202210450383.7A
Other languages
Chinese (zh)
Inventor
刘耀琪
曹安业
白贤栖
白金正
周坤友
阚吉亮
田鑫元
杨旭
王常彬
薛成春
郭文豪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Mining and Technology CUMT
Original Assignee
China University of Mining and Technology CUMT
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Mining and Technology CUMT filed Critical China University of Mining and Technology CUMT
Priority to CN202210450383.7A priority Critical patent/CN114779331A/en
Publication of CN114779331A publication Critical patent/CN114779331A/en
Pending legal-status Critical Current

Links

Images

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/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/01—Measuring or predicting earthquakes

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Remote Sensing (AREA)
  • Acoustics & Sound (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Business, Economics & Management (AREA)
  • Emergency Management (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention discloses a mine earthquake risk area prediction method based on accumulated microseismic response, which comprises the following steps: the mining earthquake risk area prediction method comprises the steps of firstly determining an area to be predicted, collecting microseismic monitoring data of the area, further counting roadway accumulated PPV distribution and determining a danger threshold value by constructing a roadway simplified model and calculating the PPV distribution of a roadway, and predicting the mining earthquake risk area according to the roadway accumulated PPV response frequency distribution. The method provided by the invention has the advantages of clear calculation model, universality and strong operability, and can effectively predict the impact display area.

Description

Mine earthquake risk area prediction method based on accumulated microseismic response
Technical Field
The invention belongs to the technical field of coal mining and coal mine safety, and particularly relates to a mine earthquake risk area prediction method based on accumulated microseismic response.
Background
Mine earthquake is generated along with coal mine production activities, and can only avoid the harm caused by the mine earthquake to the greatest extent under the limited possibility. Because the underground coal mine workers and the mechanical operation places are all located in the roadway, how to evaluate the influence of the continuously generated mine earthquake on the roadway and then predict the risk area of the mine earthquake can be taken the first place.
In recent years, microseismic systems are widely deployed in rock burst mines for capturing mine seismic signals, and clear waveforms of the seismic signals can be obtained through the microseismic systems. The microseismic waveform contains abundant mineral earthquake information, such as the amplitude, duration and vibration frequency of the waveform, which can reflect the influence of mineral earthquake on the roadway from some sides. However, how to effectively utilize the information in the waveform and fully consider the influence of the continuously generated mine earthquake on the roadway to predict the mine earthquake risk area is not solved well. In addition, a micro-seismic system is often provided with a certain number of stations around a coal mine excavation space for collecting waveforms, so that waveform information obtained from the micro-seismic system is only coal and rock mass vibration information near the stations; and when the waveform information is used for risk assessment of the roadway, vibration information of the whole roadway is required, so that how to obtain the vibration information of the whole roadway by a scientific method according to the limited microseismic station information becomes more important.
Disclosure of Invention
The invention aims to provide a mine earthquake risk area prediction method based on accumulated microseismic response, so as to solve the problems in the prior art.
In order to achieve the purpose, the invention provides a mine earthquake risk area prediction method based on accumulated microseismic response, which comprises the following steps:
acquiring microseismic monitoring data of a region to be predicted, and constructing a roadway simplified model;
and predicting a mine earthquake risk area based on the microseismic monitoring data and the roadway simplified model.
Optionally, the microseismic monitoring data includes the coordinate of a microseismic monitoring system station, the coordinate of a microseismic source, and the peak value of a PPV picked up by different microseismic stations, and the PPV is the particle peak velocity of the roadway.
Optionally, the method for constructing the simplified roadway model includes: simplifying the roadway into straight lines, and dispersing the straight lines into scattered points.
Optionally, the process of predicting the mine earthquake risk area based on the microseismic monitoring data and the roadway simplified model includes: and calculating the radius of a seismic source, calculating the PPV of the scatter points based on the size relationship between the scatter points of the simplified roadway model and the positioning of the seismic source and the radius of the seismic source, and predicting the mine earthquake risk area based on the PPV of the simplified roadway model and the scatter points.
Optionally, the process of predicting the mineral earthquake risk area based on the roadway simplified model and the scattered point PPV includes:
calculating the PPV distribution of the roadway based on the roadway simplified model;
counting the cumulative PPV distribution of the roadway;
calculating the distribution of the response times of the accumulated PPV of the roadway and determining a danger threshold;
and predicting a mine earthquake risk area based on the distribution of the roadway accumulated PPV response times.
Optionally, the source radius is calculated using the following equation:
Figure BDA0003617016750000021
wherein r is0Is the focal radius, VAIs the apparent volume of the seismic source; g is the shear modulus of the coal rock mass; m0As seismic sourcesThe seismic moment of (a); esIs the source energy.
Optionally, the calculating the PPV value of the scatter point based on the magnitude relationship between the scatter point and the source location and the source radius includes: when the linear distance between the scattered point and the seismic source is less than or equal to the radius of the roadway, calculating the PPV value of the scattered point by adopting the following formula:
ppvnear field＝1.28(Cs/G)ρRppa
Wherein CSIn terms of the shear wave velocity, ρ is the density of the propagation medium, R is the distance from the seismic source to the station, and ppa is the peak acceleration of the particle at the station, which is obtained by deriving the velocity waveform.
Optionally, the process of calculating the PPV value of the scatter based on the relationship between the scatter and the source location and the source radius further includes: and when the straight-line distance between the scattered point and the seismic source is greater than the radius of the roadway, counting the distances from the seismic source to all stations and the corresponding PPV to obtain the attenuation relation between the PPV and the distances, and obtaining the PPV value of the scattered point based on the attenuation relation.
Optionally, the method for determining the risk threshold is: selecting 80% of the PPV distribution as a threshold.
Optionally, the hazard threshold is determined taking into account a spatial relationship of a source size to the roadway.
The invention has the technical effects that:
the mine earthquake risk area is predicted by acquiring the micro-earthquake monitoring data, constructing the roadway simplified model and based on the micro-earthquake monitoring data and the roadway simplified model. The method provided by the invention has the advantages of clear calculation model, universality and strong operability, and can effectively predict the impact display area.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this application, illustrate embodiments of the application and, together with the description, serve to explain the application and are not intended to limit the application. In the drawings:
fig. 1 is a flowchart of a method for predicting a mine earthquake risk area according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of the impact development position during recovery of a working surface LW250106-1 in an embodiment of the present invention;
FIG. 3 is a schematic diagram of a research area and microseismic station arrangement in an embodiment of the present invention;
FIG. 4 is a diagram of a typical shock waveform for a microseismic event within a region of interest in accordance with an embodiment of the present invention;
FIG. 5 is a schematic diagram of the peak velocity pickup of particles for a channel (probe) in a region of interest in an embodiment of the present invention;
FIG. 6 is a schematic diagram illustrating the peak velocity and energy attenuation of the particle propagated by the mine seismic waves within and outside the range of the source radius according to an embodiment of the present invention;
FIG. 7 is a schematic view of a fit of particle peak velocity attenuation coefficients for microseismic events at a time outside the radius of the seismic source within the area of interest in accordance with an embodiment of the present invention;
FIG. 8 is a schematic diagram of counting roadway cumulative PPV distributions and determining a hazard threshold in an embodiment of the present invention;
fig. 9 is a schematic diagram of comparison between an actual impact-exhibiting area and a predicted mineral earthquake risk area according to an embodiment of the present invention.
Detailed Description
It should be noted that the embodiments and features of the embodiments in the present application may be combined with each other without conflict. The present application will be described in detail below with reference to the accompanying drawings in conjunction with embodiments.
It should be noted that the steps illustrated in the flowcharts of the figures may be performed in a computer system such as a set of computer-executable instructions and that, although a logical order is illustrated in the flowcharts, in some cases, the steps illustrated or described may be performed in an order different than here.
Example one
As shown in fig. 1, the present embodiment provides a method for predicting a mine earthquake risk area based on cumulative microseismic responses, including:
firstly, determining a region to be predicted and collecting microseismic monitoring data of the region, further counting the cumulative PPV (particle peak velocity) distribution of the roadway and determining a danger threshold by constructing a roadway simplified model and calculating the PPV (particle peak velocity) distribution of the roadway, and predicting a mine earthquake risk region according to the distribution of the cumulative PPV response times of the roadway; the method specifically comprises the following steps: the acquired microseismic data mainly comprises microseismic monitoring system station coordinates, microseismic source coordinates and PPV peak values picked up by different microseismic station; simplifying the roadway into a straight line when constructing a roadway simplified model; and when the PPV threshold value of the roadway microseismic damage is counted, the spatial relation between the size of the seismic source and the roadway is considered.
As a preferred technical scheme of the invention: when the roadway simplification model is constructed and the PPV (particle peak velocity) distribution of the roadway is calculated, the roadway is simplified into a straight line, then the straight line is further dispersed into scattered points, the radius of the seismic source is calculated, and the PPV of the scattered points is calculated according to the size relationship between the dispersed scattered points and the positioning and radius of the seismic source.
As a preferred technical scheme of the invention: the seismic source radius r0Calculated using the formula:
Figure BDA0003617016750000051
in the formula, VAIs the apparent volume of the seismic source; g is the shear modulus of the coal rock mass; m0Seismic moment as a seismic source; esIs the source energy.
As a preferred technical scheme of the invention: and when the linear distance between the scattered point and the seismic source is not more than the radius of the roadway, calculating the PPV value of the scattered point by adopting the following formula:
ppvnear field＝1.28(Cs/G)ρRppa
Wherein CSIn terms of shear wave velocity, ρ is the density of the propagation medium, R is the distance from the seismic source to the station, and ppa is the peak acceleration of the particles at the station, which can be derived from the velocity waveform.
As a preferred technical scheme of the invention: when the straight-line distance between the scattered point and the seismic source is larger than the radius of the roadway, the distances from the seismic source to all the stations and the corresponding PPV need to be counted to obtain the attenuation relation between the PPV and the distances, and the PPV value of the scattered point is obtained according to the attenuation relation.
As a preferred technical scheme of the invention: and when the statistical roadway accumulates the PPV distribution and determines the danger threshold value, 80 percent of all the PPV distributions are used as the threshold value.
Example two
As shown in fig. 2 to 9, the present embodiment provides a prediction example of a method for predicting a mine earthquake risk area based on cumulative microseismic responses, which includes:
as shown in figure 2, in 2017, 14 months and 5 days, when a working face LW250106-1 of a certain mine of Huating pavilion in Gansu is mined, an impact display event occurs together, so that the bottom drum damage is caused within a range of 197-257 m in front of a working face transportation crossheading, the height of the bottom drum reaches 0.2-0.3 m, equipment in an inner roadway is damaged, and no casualties exist. As shown in figure 3, in order to monitor the mine earthquake events during the working face mining, a microseismic station is installed in the whole mine range, particularly near the LW250106-1 working face, the microseismic station can effectively monitor and record the vibration waveform generated by the mine earthquake events (figure 4), and the origin time and space coordinates of a seismic source can be calculated through the microseismic post-processing degree.
In order to test the effectiveness of the mine earthquake risk area prediction method provided by the invention, the accident case is adopted for verification. Firstly, the LW250106-1 working face is determined to be a target area predicted by the method, and microseismic data are acquired when the LW250106-1 working face is recovered to establish a microseismic database. The microseismic data comprises microseismic data in a month before the LW250106-1 working surface is impacted, and mainly comprises seismic source coordinates, microseismic waveforms recorded by each microseismic station, coordinates of the microseismic stations and the like.
Because the area of the cross section of the underground roadway of the coal mine is negligible relative to the whole mine, the roadway can be simplified to be provided with a spatial straight line when the dimension space problem of the mine is analyzed; for the convenience of subsequent prediction calculation processing, the straight line may be further discretized into scattered points, for example, one lane has a length of 1000m, and after being simplified into the straight line, the straight line may be further simplified into scattered points consisting of 1001 scattered points and having an interval of 1 m. In practicing the present invention, the waveforms recorded at different stations in each microseismic event are first picked up and the PPV value is picked up, for example, fig. 5 is a schematic diagram of the PPV peak value of a waveform recorded at a station for a certain microseismic event. And after the PPV values of all the data in the microseismic database are picked up, calculating the cumulative effect of the microseismic event on the roadway. Specifically, the correlation relationship between the scattered points of the roadway after dispersion and the positions of the seismic sources and the radii of the seismic sources is firstly determined, for example, as shown in fig. 6, if the scattered points after dispersion (left dark regions in the figure) are located within the range of the radii of the seismic sources, the PPV values of the scattered points are calculated by using the formula (1).
ppvNear field＝1.28(Cs/G)ρRppa
Wherein C isSIn terms of shear wave velocity, ρ is the density of the propagation medium, R is the distance from the seismic source to the station, and ppa is the peak acceleration of the particle at the station, which can be derived from the velocity waveform.
In addition, as a preferable embodiment of the present invention: the seismic source radius r0Calculated using the formula:
Figure BDA0003617016750000071
in the formula, VAIs the apparent volume of the seismic source; g is the shear modulus of the coal rock mass; m is a group of0Seismic moments of a seismic source; esIs the source energy.
If the scattered point is located outside the radius of the seismic source, firstly, the attenuation relations between the PPV of all the received microseismic waveforms of the microseismic event and the propagation distance are counted, and after an attenuation relation expression (shown in figure 7) is obtained, the PPV numerical values of all the scattered points are calculated by adopting the attenuation relation expression according to the distance between the scattered points and the center of the seismic source.
Further, the risk threshold is determined by the PPV caused by all microseismic events in the database to the nearest point of the roadway, wherein after a lot of experimental tests by the inventor, as shown in fig. 8, the prediction effect obtained when 80% of all statistics is taken as the critical threshold is the best, so that the PPV value corresponding to 80% of the ratio in the cumulative frequency is taken as the critical threshold, and the critical PPV threshold is 0.01m/s in the embodiment.
Further, the statistical distribution of the times of PPV (pulse per square volt) larger than 0.01m/s caused by the microseismic events to roadway scattered points in the microseismic database is counted. Fig. 9 shows the statistical dispersion distribution of the roadway with PPV greater than 0.01m/s in this embodiment, which is used as the predicted mineral earthquake risk area. From fig. 9, it can be found that the 250106-1 working face impact display area (dark area in fig. 9) falls into the predicted mineral shock risk area and is obviously located at a strong shock frequency (PPV greater than 0.01m/s), which indicates that the prediction method provided by the invention completely predicts the impact display area and has good prediction effect.
The above description is only for the preferred embodiment of the present application, but the scope of the present application is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present application should be covered within the scope of the present application. Therefore, the protection scope of the present application shall be subject to the protection scope of the claims.

Claims (9)

1. A mine earthquake risk area prediction method based on accumulated microseismic response is characterized by comprising the following steps:
acquiring microseismic monitoring data of a region to be predicted, and constructing a roadway simplified model;
and predicting a mine earthquake risk area based on the microseismic monitoring data and the roadway simplified model.
2. The method for predicting the mine earthquake risk area based on the cumulative microseismic response of claim 1, wherein the method comprises the following steps: the microseismic monitoring data comprise microseismic monitoring system station coordinates, microseismic source coordinates and the peak value of PPV (particle per volt) picked up by different microseismic stations, and the PPV is the particle peak value speed of the roadway.
3. The method for predicting the mine earthquake risk area based on the cumulative microseismic response as recited in claim 1, wherein: the construction method of the roadway simplified model comprises the following steps: simplifying the roadway into straight lines, and dispersing the straight lines into scattered points.
4. The method for predicting the mine earthquake risk area based on the cumulative microseismic response of claim 1, wherein the method comprises the following steps: the process for predicting the mine earthquake risk area based on the microseismic monitoring data and the roadway simplified model comprises the following steps: and calculating the radius of a seismic source, calculating the PPV of the scattered points based on the size relationship between the scattered points of the roadway simplified model and the positioning of the seismic source and the radius of the seismic source, and predicting the mine earthquake risk area based on the PPV of the roadway simplified model and the scattered points.
5. The method for predicting the mine earthquake risk area based on the cumulative microseismic response of claim 1, wherein the method comprises the following steps: the process for predicting the mine earthquake risk area based on the roadway simplified model and the PPV of the scattered points comprises the following steps:
calculating the PPV distribution of the roadway based on the roadway simplified model;
counting the cumulative PPV distribution of the roadway;
calculating the distribution of the response times of the accumulated PPV of the roadway and determining a danger threshold;
and predicting a mine earthquake risk area based on the distribution of the roadway accumulated PPV response times.
6. The method for predicting the mine earthquake risk area based on the cumulative microseismic response as recited in claim 4, wherein: the source radius is calculated using the following equation:
Figure FDA0003617016740000021
wherein r is0Is the focal radius, VAIs the apparent volume of the seismic source; g is the shear modulus of the coal rock mass; m0Seismic moment as a seismic source; esIs the source energy.
7. The method for predicting the mine earthquake risk area based on the cumulative microseismic response as recited in claim 4, wherein: the process of calculating the PPV value of the scatter based on the size relation between the scatter and the seismic source location and the seismic source radius comprises the following steps: when the linear distance between the scattered point and the seismic source is less than or equal to the radius of the roadway, calculating the PPV value of the scattered point by adopting the following formula:
ppvnear field＝1.28(Cs/G)ρRppa
Wherein CSIn terms of the shear wave velocity, ρ is the density of the propagation medium, R is the distance from the seismic source to the station, and ppa is the peak acceleration of the particle at the station, which is obtained by deriving the velocity waveform.
8. The method for predicting the mine earthquake risk area based on the cumulative microseismic response as recited in claim 4, wherein: the process of calculating the PPV value of the scatter point based on the size relation between the scatter point and the seismic source location and the seismic source radius further comprises the following steps: and when the straight-line distance between the scattered point and the seismic source is greater than the radius of the roadway, counting the distances from the seismic source to all the stations and the corresponding PPV to obtain the attenuation relation between the PPV and the distances, and obtaining the PPV value of the scattered point based on the attenuation relation.
9. The method for predicting the mine earthquake risk area based on the cumulative microseismic response as recited in claim 5, wherein: the method for determining the danger threshold comprises the following steps: and determining a danger threshold value based on the spatial relationship between the size of the seismic source and the roadway, and selecting 80% of the PPV distribution as the danger threshold value.
CN202210450383.7A 2022-04-26 2022-04-26 Mine earthquake risk area prediction method based on accumulated microseismic response Pending CN114779331A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210450383.7A CN114779331A (en) 2022-04-26 2022-04-26 Mine earthquake risk area prediction method based on accumulated microseismic response

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210450383.7A CN114779331A (en) 2022-04-26 2022-04-26 Mine earthquake risk area prediction method based on accumulated microseismic response

Publications (1)

Publication Number Publication Date
CN114779331A true CN114779331A (en) 2022-07-22

Family

ID=82433886

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210450383.7A Pending CN114779331A (en) 2022-04-26 2022-04-26 Mine earthquake risk area prediction method based on accumulated microseismic response

Country Status (1)

Country Link
CN (1) CN114779331A (en)

Similar Documents

Publication Publication Date Title
US20230168406A1 (en) Rock burst hazard prediction method based on seismic wave energy attenuation characteristics of mine earthquake cluster
Wang et al. Acoustic emission/microseismic source location analysis for a limestone mine exhibiting high horizontal stresses
Yang et al. Frequency comparison of blast-induced vibration per delay for the full-face millisecond delay blasting in underground opening excavation
Shen et al. Displacement, stress and seismicity in roadway roofs during mining-induced failure
CN111648826B (en) Coal mine rock burst prediction early warning system and method
CN115755185B (en) Method for judging disaster causing performance of high-energy ore earthquake based on microseism monitoring
Li et al. Rockburst occurrences and microseismicity in a longwall panel experiencing frequent rockbursts
CN108594295B (en) Method for evaluating coal bed blasting pressure relief effect
Wang et al. Location error based seismic cluster analysis and its application to burst damage assessment in underground coal mines
Potvin Strategies and tactics to control seismic risks in mines
Shirzadegan et al. Large scale dynamic testing of rock support system at Kiirunavaara underground mine
CN112324506B (en) Dynamic early warning method for preventing rock burst of coal mine based on microseism
CN114412573B (en) Mine wall breakage monitoring and early warning system and method
Jiang et al. Combined early warning method for rockburst in a Deep Island, fully mechanized caving face
CN114994791B (en) Method for evaluating monitoring capability of well-ground integrated micro-seismic monitoring system
CN114109508A (en) Rock burst monitoring method based on energy system analysis and application
CN116842411A (en) Coal mine microseismic platform network monitoring capability evaluation method based on vibration wave detection probability
CN109239775B (en) Method for tracking and positioning stolen mining of mineral resources
Wang et al. Ground motion characteristics and their cumulative impacts to burst risks in underground coal mines
Chen et al. Analysis of ground vibrations caused by open pit production blasts–a case study
CN114779331A (en) Mine earthquake risk area prediction method based on accumulated microseismic response
Vinoth et al. Slope stability monitoring by quantification and behavior of microseismic events in an opencast coal mine
Westley-Hauta et al. Characterisation of seismic activity at a kimberlite block caving operation in a complex geological setting in Quebec, Canada
Giacomini et al. In situ experiments of rockfall in an open pit coal mine
CN114704329B (en) Method for monitoring earthquake waves of deep tunnel tunneling mine and rapidly detecting advanced stress

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