CN114185084A - Method and device for determining influence range of bottoming operation, electronic equipment and medium - Google Patents

Method and device for determining influence range of bottoming operation, electronic equipment and medium Download PDF

Info

Publication number
CN114185084A
CN114185084A CN202210143734.XA CN202210143734A CN114185084A CN 114185084 A CN114185084 A CN 114185084A CN 202210143734 A CN202210143734 A CN 202210143734A CN 114185084 A CN114185084 A CN 114185084A
Authority
CN
China
Prior art keywords
line
determining
microseismic events
partition
influence range
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.)
Granted
Application number
CN202210143734.XA
Other languages
Chinese (zh)
Other versions
CN114185084B (en
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.)
BGRIMM Technology Group Co Ltd
Original Assignee
BGRIMM Technology Group Co Ltd
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 BGRIMM Technology Group Co Ltd filed Critical BGRIMM Technology Group Co Ltd
Priority to CN202210143734.XA priority Critical patent/CN114185084B/en
Publication of CN114185084A publication Critical patent/CN114185084A/en
Application granted granted Critical
Publication of CN114185084B publication Critical patent/CN114185084B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/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/30—Analysis
    • G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Acoustics & Sound (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 provides a method and a device for determining the influence range of bottoming operation, electronic equipment and a medium, wherein the method comprises the following steps: for each bottoming period, acquiring a vibration signal in the bottoming period through a micro-seismic monitoring system, and determining the seismic source position and the seismic magnitude of a micro-seismic event in each bottoming period; dividing the bottoming operation surface into a plurality of partitions based on the seismic source position of the microseismic event, and calculating the distance between each partition and a bottoming line; wherein, the subregion includes: at least one pull-bottom line front partition and at least one pull-bottom line rear partition; determining the number of microseismic events in each partition, and determining the bottom pulling operation influence range of the bottom pulling period based on the number of the microseismic events in each partition and the magnitude of the microseismic events; and determining the average value of the influence ranges of the bottom pulling operation in a plurality of bottom pulling periods as the target influence range of the bottom pulling operation. The invention can improve the accuracy of dividing the influence range of the bottom pulling operation.

Description

Method and device for determining influence range of bottoming operation, electronic equipment and medium
Technical Field
The invention relates to the technical field of mine safety, in particular to a method and a device for determining an influence range of bottom pulling operation, electronic equipment and a medium.
Background
The metal mine faces the ground pressure problem in the mining process, the stability of surrounding rocks is influenced by blasting and ore removal operation, and the safety of underground personnel and equipment can be threatened by ground pressure disasters. The mine mined by the natural caving method is subjected to bottom-drawing operation to break the stress balance state formed by surrounding rocks of the bottom structure, and the stress redistribution stage is accompanied with the fracture deformation of a rock body. Meanwhile, the single explosive loading of the bottom-pulling blasting is large, the high-energy blasting shock wave can aggravate the cracking degree of the surrounding rock, and once the large-scale through damage is formed, the ground pressure disaster can be caused. Therefore, the influence range of the bottoming operation on the stability of the surrounding rock is analyzed, the cracking and development degree of the surrounding rock is evaluated, and the method has great significance for safe production of mines. At present, a method for analyzing the stability of surrounding rocks by using natural caving bottom-drawing operation mainly comprises numerical simulation, building a numerical simulation model by using mine engineering parameters, analyzing the stress state of the surrounding rocks of a bottom structure after each step of bottom drawing, and summarizing the time-space evolution characteristics and rules of stress at each stage. However, the analysis result of the numerical simulation method is difficult to be consistent with the actual ground pressure display rule of the mine site, so that the accuracy of dividing the influence range of the bottom pulling operation is influenced.
Disclosure of Invention
In view of the above, an object of the present invention is to provide a method, an apparatus, an electronic device and a medium for determining an influence range of a pull-down operation, so as to divide the influence range of the pull-down operation more accurately.
In order to achieve the above purpose, the embodiment of the present invention adopts the following technical solutions:
in a first aspect, an embodiment of the present invention provides a method for determining an influence range of a bottom pulling operation, including: for each bottoming period, acquiring a vibration signal in the bottoming period through a micro-seismic monitoring system, and determining the seismic source position and the seismic magnitude of a micro-seismic event in each bottoming period; dividing the bottoming operation surface into a plurality of partitions based on the seismic source position of the microseismic event, and calculating the distance between each partition and a bottoming line; wherein, the subregion includes: at least one pull-bottom line front partition and at least one pull-bottom line rear partition; determining the number of microseismic events in each partition, and determining the bottom pulling operation influence range of the bottom pulling period based on the number of the microseismic events in each partition and the magnitude of the microseismic events; and determining the average value of the influence ranges of the bottom pulling operation in a plurality of bottom pulling periods as the target influence range of the bottom pulling operation.
In one embodiment, dividing the bottoming job surface into a plurality of partitions based on the source locations of the microseismic events comprises: determining areas at two sides of a preset reference point and a first preset distance from the reference point along the bottom pulling tendency as an inclination range; wherein, the datum point is the midpoint of the bottom pulling line when the bottom pulling period is finished; based on the seismic source position of the microseismic event, in the inclination range, the reference point is used as a starting point, the front area of the bottom pulling line is divided into at least one front subarea of the bottom pulling line according to a preset distance, and the rear area of the bottom pulling line is divided into at least one rear subarea of the bottom pulling line.
In one embodiment, the pull-down operation impact range includes: the influence range in front of the bottom drawing line and the influence range behind the bottom drawing line; determining the number of microseismic events in each partition, and determining the influence range of the bottom pulling operation of the bottom pulling period based on the number of the microseismic events in each partition and the magnitude of the microseismic events, wherein the influence range comprises the following steps: determining the number of microseismic events in each partition and the total number of microseismic events in the tendency range; calculating the ratio of the number of the microseismic events in each subarea to the number of the first microseismic events in the total number of the microseismic events, and determining the front influence range of the first bottoming line and the rear influence range of the first bottoming line based on the ratio of the number of the first microseismic events in each subarea; determining the number of microseismic events with the magnitude greater than a magnitude threshold in each zone based on the magnitude of the microseismic events; wherein, the magnitude threshold is the median of the magnitude of the microseismic event in the tendency range; calculating the second microseismic event number ratio of the number of microseismic events with the seismic level larger than the seismic level threshold value in each subarea to the total number of the microseismic events, and determining the front influence range of the second bottom pulling line and the rear influence range of the second bottom pulling line based on the second microseismic event number ratio in each subarea; and determining the maximum value in the first bottom pulling wire front influence range and the second bottom pulling wire front influence range as the bottom pulling wire front influence range, and determining the maximum value in the first bottom pulling wire rear influence range and the second bottom pulling wire rear influence range as the bottom pulling wire rear influence range.
In one embodiment, calculating a first ratio of the number of microseismic events in each zone to the total number of microseismic events and determining a first pull-bottom front impact range and a first pull-bottom rear impact range based on the first ratio of the number of microseismic events in each zone comprises: calculating a first percentage of the number of microseismic events of the partition in front of each bedding line in the total number of microseismic events and a second percentage of the number of microseismic events of the partition in back of each bedding line in the total number of microseismic events; determining the front partition of the bottom pulling line with the first percentage larger than a first threshold value as a first front target area, and determining the maximum distance between the first front target area and the bottom pulling line as the front influence range of the first bottom pulling line; and determining the rear partition of the bottom pulling line with the second percentage larger than the first threshold as a first rear target area, and determining the maximum distance between the first rear target area and the bottom pulling line as the rear influence range of the first bottom pulling line.
In one embodiment, calculating a second ratio of the number of microseismic events with a magnitude greater than a magnitude threshold in each zone to the total number of microseismic events, and determining a second pull-bottom line front impact range and a second pull-bottom line rear impact range based on the second ratio of the number of microseismic events in each zone comprises: calculating a third percentage of the number of microseismic events with the seismic level larger than the seismic level threshold value in the partition in front of each bottom pulling line in the total number of the microseismic events and a fourth percentage of the number of the microseismic events with the seismic level larger than the seismic level threshold value in the partition behind each bottom pulling line in the total number of the microseismic events; determining the front partition of the bottom pulling line with the third percentage larger than the second threshold as a second front target area, and determining the maximum distance between the second front target area and the bottom pulling line as the front influence range of the second bottom pulling line; and determining the rear partition of the bottom pulling line with the fourth percentage larger than the second threshold as a second rear target area, and determining the maximum distance between the second rear target area and the bottom pulling line as the front influence range of the second bottom pulling line.
In one embodiment, determining the average value of the pull-down operation influence ranges of a plurality of pull-down cycles as the pull-down operation target influence range comprises: determining the average value of the front influence ranges of the bottom pulling lines in a plurality of bottom pulling periods as the front target influence range of the bottom pulling operation; and determining the average value of the influence ranges behind the bottom pulling lines of the plurality of bottom pulling periods as the target influence range behind the bottom pulling operation.
In one embodiment, prior to determining the source location and magnitude of the microseismic event within each bottoming cycle, the method further comprises: and eliminating blasting signals and noise signals in the vibration signals.
In a second aspect, an embodiment of the present invention provides an apparatus for determining an influence range of a bottoming operation, including: the microseismic event determining module is used for acquiring a vibration signal in each bottoming period through the microseismic monitoring system and determining the seismic source position and the seismic magnitude of the microseismic event in each bottoming period; the partition dividing module is used for dividing the bottoming operation surface into a plurality of partitions based on the seismic source position of the microseismic event and calculating the distance between each partition and a bottoming line; wherein, the subregion includes: at least one pull-bottom line front partition and at least one pull-bottom line rear partition; the first determining module is used for determining the number of microseismic events in each partition, and determining the influence range of the bottom pulling operation in the bottom pulling period based on the number of the microseismic events in each partition and the magnitude of the microseismic events; and the second determining module is used for determining the average value of the influence ranges of the bottom pulling operation in the plurality of bottom pulling periods as the target influence range of the bottom pulling operation.
In a third aspect, an embodiment of the present invention provides an electronic device, which includes a processor and a memory, where the memory stores computer-executable instructions capable of being executed by the processor, and the processor executes the computer-executable instructions to implement the steps of any one of the methods provided in the first aspect.
In a fourth aspect, the present invention provides a computer-readable storage medium, on which a computer program is stored, where the computer program is executed by a processor to perform the steps of any one of the methods provided in the first aspect.
The embodiment of the invention has the following beneficial effects:
according to the method, the device, the electronic equipment and the medium for determining the influence range of the bottom pulling operation, provided by the embodiment of the invention, firstly, for each bottom pulling period, a vibration signal in the bottom pulling period is obtained through a micro-seismic monitoring system, and the seismic source position and the seismic magnitude of a micro-seismic event in each bottom pulling period are determined; then, dividing the bottoming operation surface into a plurality of partitions (at least one partition in front of the bottoming line and at least one partition behind the bottoming line) based on the seismic source position of the microseismic event, and calculating the distance between each partition and the bottoming line; then, determining the number of microseismic events in each partition, and determining the bottom pulling operation influence range of the bottom pulling period based on the number of the microseismic events in each partition and the magnitude of the microseismic events; and finally, determining the average value of the influence ranges of the bottom pulling operation in a plurality of bottom pulling periods as the target influence range of the bottom pulling operation. The method can divide areas according to the seismic source position of the microseismic event, and can simultaneously determine the influence ranges of the bottom pulling operation of the natural caving method on the bottom pulling operation in front of and behind the bottom pulling line, thereby accurately describing the actual ground pressure display rule of the mine site; secondly, the method determines the influence range of the bottom pulling operation by the natural caving method based on the number of the microseismic events and the distribution of the seismic magnitude, and eliminates the interference of blasting vibration on the total energy of the regional vibration signals, so that the division of the influence range of the bottom pulling operation is more accurate.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
In order to make the aforementioned and other objects, features and advantages of the present invention comprehensible, preferred embodiments accompanied with figures are described in detail below.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.
Fig. 1 is a flowchart of a method for determining an influence range of a bottoming operation according to an embodiment of the present invention;
FIG. 2 is a schematic diagram illustrating partition division of a bottoming operation plane according to an embodiment of the present invention;
FIG. 3 is a schematic structural diagram of an apparatus for determining an influence range of a bottoming operation according to an embodiment of the present invention;
fig. 4 is a schematic structural diagram of an electronic device according to an embodiment of the present invention.
Detailed Description
To make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings, and it is apparent that the described embodiments are some, but not all embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
At present, numerical simulation analysis is a common method for analyzing the stability of surrounding rocks by bottoming operation, and the method utilizes numerical simulation software to carry out excavation simulation in steps according to the actual bottoming operation sequence, and analyzes the stress state and the time-space evolution characteristics of a bottom structure after each step of bottoming. The model established by the numerical simulation method simplifies the actual engineering conditions, the rock parameters obtained by the laboratory test of the rock test piece have certain errors, and the difference between the initial stress and the boundary condition and the actual complex stress state of the mine is large, so the analysis result is difficult to be consistent with the actual ground pressure display rule of the mine site, the practical guiding significance is lacked for the safety of mining, and the accuracy of dividing the influence range of the bottoming operation is influenced; secondly, rock mass parameters of the numerical simulation method are acquired for one time in an experiment, after initial stress and boundary conditions are set, the rock mass parameters cannot be changed along with the stable state of the ore rock in the mining process, and the production operation arrangement of a mine site is changed at any time, so that the model cannot reflect the real-time change of deformation and fracture of the ore rock, and an analysis result does not have the function of monitoring and early warning.
Based on this, the method, the device, the electronic device and the medium for determining the influence range of the pull-down operation provided by the embodiment of the invention enable the division of the influence range of the pull-down operation to be more accurate.
To facilitate understanding of the embodiment, a detailed description will be given to a method for determining an influence range of a pull-down operation disclosed in the embodiment of the present invention, which can be executed by an electronic device, such as a computer, a smart phone, and the like. Referring to the flowchart of the method for determining the influence range of the bottom pulling operation shown in fig. 1, it is shown that the method mainly includes the following steps S101 to S104:
step S101: and for each pull-down period, acquiring a vibration signal in the pull-down period through a microseismic monitoring system, and determining the seismic source position and the seismic magnitude of a microseismic event in each pull-down period.
In one embodiment, one bottoming cycle is from the beginning of a broached face blast penetration to before the next blast penetration of the face. Specifically, the microseismic monitoring system can be used for collecting vibration signals in a bottoming period, eliminating bottoming blasting signals and noise signals, processing microseismic events and determining the seismic source positions and the seismic levels of all the microseismic events in the bottoming period. In particular, it can be obtained by known geophysical seismic source parameter calculation methods, such as: after the microseismic monitoring system records microseismic signals, the arrival time of P \ S waves of seismic waves is calibrated manually or automatically by using professional processing software, the position of a seismic source is calculated through sensor coordinates and sensor trigger time, and the seismic level is calculated through seismic waveforms recorded by the sensors.
Step S102: dividing the bottoming operation surface into a plurality of partitions based on the seismic source position of the microseismic event, and calculating the distance between each partition and a bottoming line; wherein, the subregion includes: at least one pull-bottom line front partition and at least one pull-bottom line rear partition.
In one embodiment, the seismic source positions of the microseismic events can be determined, and the microseismic events before and after the microseismic events can be partitioned by taking the bottom line as a reference, the partition needs to cover the microseismic events at the farthest ends before and after the bottom line, and then the distance from each partition to the bottom line is calculated.
Step S103: and determining the number of microseismic events in each partition, and determining the pull-down operation influence range of the pull-down period based on the number of microseismic events in each partition and the magnitude of the microseismic events.
In one embodiment, the bottoming operation influence range may be determined according to the number of microseismic events in each partition, specifically, the number of microseismic events in each partition and the total number of all microseismic events in the trend range are counted, the percentage of the number of microseismic events in each partition in the total number of the microseismic events is calculated, and the area where the value of the percentage exceeds the threshold is determined as the bottoming operation influence range; then determining the bottom pulling operation influence range of the bottom pulling period according to the magnitude of the microseismic events, specifically, counting the magnitude of the microseismic events in each partition, calculating the magnitude median of all the microseismic events in the inclination range, counting the number of the microseismic events of which the magnitude is greater than the magnitude median in each partition, calculating the percentage of the number of the microseismic events of which the magnitude is greater than the magnitude median in each partition to the total number of the microseismic events in the inclination range, and determining the area of which the percentage value exceeds the threshold as the bottom pulling operation influence range; and finally, taking the maximum value in the two bottoming operation influence ranges as the bottoming operation influence range of the bottoming period.
Step S104: and determining the average value of the influence ranges of the bottom pulling operation in a plurality of bottom pulling periods as the target influence range of the bottom pulling operation.
In an embodiment, an average value of the pull-down operation influence ranges of a plurality of continuous pull-down cycles may be taken as a final pull-down operation influence range, that is, a pull-down operation target influence range, in this embodiment, the number of continuous pull-down cycles is not less than 2.
The method for determining the influence range of the bottom pulling operation provided by the embodiment of the invention can be used for carrying out region division according to the seismic source position of the microseismic event and simultaneously determining the influence ranges of the bottom pulling operation of the natural caving method on the front and the back of the bottom pulling line, thereby accurately describing the actual ground pressure display rule of the mine site; secondly, the method determines the influence range of the bottom pulling operation by the natural caving method based on the number of the microseismic events and the distribution of the seismic magnitude, and eliminates the interference of blasting vibration on the total energy of the regional vibration signals, so that the division of the influence range of the bottom pulling operation is more accurate.
In one embodiment, for the aforementioned step S102, i.e., when dividing the bottoming plane into a plurality of partitions based on the source location of the microseismic event, the following methods may be adopted, including but not limited to:
firstly, determining areas at two sides of a preset reference point and a first preset distance from the reference point along the bottom pulling tendency as an inclined range; wherein, the reference point is the midpoint of the bottom pulling line at the end of the bottom pulling period.
In specific application, the midpoint position of the pull-bottom line at the end of the pull-bottom period is determined as a reference point, and the distance between two sides of the reference point along the pull-bottom trendαThe area within the range is defined as the tendency range.
Then, based on the seismic source position of the microseismic event, in the inclination range, taking the reference point as a starting point, dividing the front area of the bottom line into at least one front partition of the bottom line according to a preset distance, and dividing the rear area of the bottom line into at least one rear partition of the bottom line.
In a specific application, see the bottoming job face partition shown in FIG. 2Dividing the schematic diagram, based on the location of the seismic source of the microseismic event, taking the reference point as a starting point in the inclination range according to a preset distanceβRespectively carrying out region division on the front and the rear of the bottom pulling line until the last microseismic event at the front and the rear ends of the bottom pulling line is covered by the regions, and obtaining the front and the rear of the bottom pulling linenThe front partition of the bottom drawing line and the rear partition of the bottom drawing line are sharedmA zone behind the pull bottom line; distance between two sides of reference point along bottom-pulling inclinationαThe value can be 15 m-20 m, and the preset intervalβThe value can be 15 m-20 m.
Further, calculating the distance from the partition in front of each bottom pulling line to the bottom pulling lineX fi (ii) a Calculating the distance from each zone behind the bottom drawing line to the bottom drawing lineX bj (ii) a The specific calculation method is as follows:
Figure P_220215170831107_107238001
Figure F_220215170830047_047150001
wherein,ithe serial number of the partition in front of the bottom drawing line is shown,nindicates the total number of partitions in front of the pull-down line,jthe serial number of the partition behind the pull bottom line is shown,mrepresenting the total number of partitions behind the pull bottom line,X fi is shown asiThe distance from the front subarea of each bottom pulling line to the bottom pulling line,X bj is shown asiThe distance from the rear subarea of each bottom drawing line to the bottom drawing line.
Considering that surrounding rocks of the bottom structure in front of and behind the natural caving faraday working surface have long service life requirements, stability monitoring is needed to realize monitoring and early warning of mine ground pressure disasters. Therefore, the embodiment of the invention can simultaneously determine the influence range of the natural caving method bottom pulling operation on the stability of surrounding rocks in front of and behind the pushing line, monitor, analyze and count microseismic events in front of and behind the bottom pulling working face, and divide the influence range of the stability of the surrounding rocks in front of and behind the pushing line. In one embodiment, the pull-down operation impact range includes: the influence range in front of the bottom drawing line and the influence range behind the bottom drawing line; for the aforementioned step S103, that is, when determining the number of microseismic events in each partition and determining the bottom pulling operation influence of the bottom pulling period based on the number of microseismic events in each partition and the magnitude of the microseismic events, the following methods, including but not limited to, may be adopted, and mainly include the following steps 1 to 5:
step 1: the number of microseismic events within each zone is determined, as well as the total number of microseismic events within the range of dip.
In one embodiment, the number of microseismic events in each partition is counted according to the partitions and the trend range determined in step S102q k And the total number of microseismic events for all microseismic events in the range of the trendQ. Specifically, the total number of microseismic events may be calculated according to the following formula:
Figure F_220215170830172_172673002
wherein,q k is shown askThe number of microseismic events for each zone,kthe sum of the number of the subareas in front of the bottom drawing line and the subareas behind the bottom drawing line.
Step 2: and calculating the ratio of the number of the microseismic events in each subarea to the number of the first microseismic events in the total number of the microseismic events, and determining the front influence range of the first bottoming line and the rear influence range of the first bottoming line based on the ratio of the number of the first microseismic events in each subarea.
In one embodiment, the number of microseismic events per zone may be calculatedq k Account for the total number of microseismic eventsQIn percentage (b)δWill beδIs greater than a first threshold valueλIs determined as the bottoming operation influence range (i.e., the first bottoming line front influence range and the first bottoming line rear influence range). In particular applications, the first pull-bottom line front and first pull-bottom line rear impact ranges may be determined in a manner including, but not limited to:
first, a first percentage of the number of microseismic events per bedding line forward partition to the total number of microseismic events and a second percentage of the number of microseismic events per bedding line rearward partition to the total number of microseismic events are calculated.
Specifically, the number of microseismic events per zone may be calculated according to the following formulaq k Account for the total number of microseismic eventsQIn percentage (b)δ：
Figure M_220215170831138_138465001
The percentage of the number of the microseismic events of the partition in front of the pull-bottom line to the total number of the microseismic events is a first percentage, and the percentage of the number of the microseismic events of the partition in back of the pull-bottom line to the total number of the microseismic events is a second percentage.
Then, the front partition of the pull-bottom line with the first percentage larger than the first threshold is determined as a first front target area, and the maximum distance between the first front target area and the pull-bottom line is determined as a first pull-bottom line front influence range.
In a particular application, the first threshold valueλThe value of (b) can be 10% -15%. For each of the pull-bottom line front sections, if the first percentage thereof is greater than the first threshold, the pull-bottom line front section is determined to be the first front target area. If the first percentage of only one front partition of the pull-bottom line is greater than the first threshold, the front partition of the pull-bottom line is the first front target area, and the distance between the partition and the pull-bottom line can be directly determined as the front influence range of the first pull-bottom lineL f (ii) a If the first percentage of the plurality of front partition areas of the pulling bottom line is larger than the first threshold value, the plurality of front partition areas of the pulling bottom line are determined as first front target areas, and then the maximum distance between the plurality of front partition areas of the pulling bottom line and the pulling bottom line is determined as the front influence range of the first pulling bottom lineL f . For example, referring to FIG. 2, for the front part of the pull-bottom line, the front sections of the pull-bottom line are F1, F2, F3, F4, F5 and F6, and the distance from each section to the reference point (i.e. the pull-bottom line) is 15m, 30m, 45m, 60m, 75m and 90m, respectively, if the pull-bottom line is pulledThe first percentages of the in-front bay sections F1, F2 are both greater than the first threshold, then the first front target zones are the in-front pull-bottom bay sections F1 and F2, and the first in-front pull-bottom influence zoneL f I.e., the distance of the lower wire front section F2 from the reference point is 30 m.
And finally, determining the rear partition of the bottom pulling line with the second percentage larger than the first threshold as a first rear target area, and determining the maximum distance between the first rear target area and the bottom pulling line as the rear influence range of the first bottom pulling line.
For each of the pullup line rear sections, if the second percentage thereof is greater than the first threshold, the pullup line rear section is determined to be the first rear target area. If the first percentage of only one partition behind the pull-bottom line is greater than the first threshold, the partition behind the pull-bottom line is the first rear target area, and the distance between the partition and the pull-bottom line can be directly determined as the influence range behind the first pull-bottom lineL b (ii) a If the first percentage of the plurality of back drawn line subareas is larger than the first threshold value, the plurality of back drawn line subareas are determined as first back target areas, and then the maximum distance between the plurality of back drawn line subareas and the drawn line is determined as the first back influence rangeL b . For example, referring to fig. 2, for the area behind the pull-bottom line, the rear partition of each pull-bottom line is B1, B2, B3, B4 and B5, and the distance from each partition to the reference point is 15m, 30m, 45m, 60m and 75m respectively; if the first percentages of the backfire rear partitions B1, B2 are both greater than the first threshold, then the first rear target zones are the backfire rear partitions B1 and B2, the first backfire rear reachL b I.e., the distance of the rear draw line section B2 from the reference point is 30 m.
And step 3: determining the number of microseismic events with the magnitude greater than a magnitude threshold in each zone based on the magnitude of the microseismic events; wherein the magnitude threshold is the median of magnitudes of microseismic events within the trend range.
In specific application, the magnitude of the microseismic events in each partition is counted, and the magnitudes of all the microseismic events in the inclination range are takenMedian numberM mad As magnitude threshold, then count magnitude greater than magnitude in each partitionM mad Number of microseismic eventse k 。
And 4, step 4: calculating the second microseismic event number ratio of the number of microseismic events with the seismic level larger than the seismic level threshold value in each subarea to the total number of the microseismic events, and determining the front influence range of the second bottom pulling line and the rear influence range of the second bottom pulling line based on the second microseismic event number ratio in each subarea;
in a specific application, the magnitude of the earthquake in each partition is calculated to be larger thanM mad Number of microseismic eventse k Total number of all microseismic events in trend rangeQIn percentage (b)ζPercent ofζIs greater than a second threshold value
Figure SYM_220215170830001
Is determined as the bottoming work influence range (i.e., the second bottoming line front influence range and the second bottoming line rear influence range). In particular applications, the second pull-bottom line front and second pull-bottom line rear impact ranges may be determined in a manner including, but not limited to:
first, the third percentage of the total number of microseismic events with the magnitude greater than the magnitude threshold in each bottom-line front partition and the fourth percentage of the total number of microseismic events with the magnitude greater than the magnitude threshold in each bottom-line rear partition are calculated.
Specifically, the magnitude of each partition greater than the magnitude of each partition can be calculated according to the following formulaM mad Number of microseismic eventse k Account for the total number of microseismic eventsQIn percentage (b)ζ：
Figure M_220215170831249_249898001
The percentage of the number of the microseismic events with the seismic level of the partition in front of the bottom pulling line larger than the seismic level threshold value in the total number of the microseismic events is a third percentage, and the percentage of the number of the microseismic events with the seismic level of the partition in back of the bottom pulling line larger than the seismic level threshold value in the total number of the microseismic events is a fourth percentage.
And then, determining the front partition of the bottom pulling line with the third percentage larger than the second threshold value as a second front target area, and determining the maximum distance between the second front target area and the bottom pulling line as the front influence range of the second bottom pulling line.
In a particular application, the second threshold value
Figure SYM_220215170831001
The value of (b) can be 5% -10%. For each of the pull-bottom line front sections, if the third percentage thereof is greater than the second threshold, the pull-bottom line front section is determined to be the second front target area. If the third percentage of only one front partition of the bottom drawing line is greater than the second threshold, the front partition of the bottom drawing line is the second front target area, and the distance between the partition and the bottom drawing line can be directly determined as the front influence range of the second bottom drawing lineD f (ii) a If the third percentage of the front subareas of the plurality of bottom drawing lines is larger than the second threshold value, the front subareas of the plurality of bottom drawing lines are determined as a second front target area, and then the maximum distance between the front subareas of the plurality of bottom drawing lines and the bottom drawing lines is determined as the front influence range of the second bottom drawing lineD f . For example, referring to fig. 2, for the front pull-bottom line area, the front pull-bottom line partitions are F1, F2, F3, F4, F5 and F6, the distances from the reference point (i.e., the pull-bottom line) of each partition are 15m, 30m, 45m, 60m, 75m and 90m, respectively, if the third percentages of the front pull-bottom line partitions F1 and F2 are greater than the second threshold, the second front target area is the front pull-bottom line partitions F1 and F2, and the second front pull-bottom line influence range is the second front pull-bottom line influence rangeD f I.e., the distance of the lower wire front section F2 from the reference point is 30 m.
And finally, determining the rear partition of the bottom pulling line with the fourth percentage larger than the second threshold as a second rear target area, and determining the maximum distance between the second rear target area and the bottom pulling line as the front influence range of the second bottom pulling line.
For eachAnd determining the back-drawn subarea as a second back target area if the fourth percentage of the back-drawn subarea is greater than the second threshold. If the fourth percentage of only one partition behind the pull bottom line is greater than the second threshold, the partition behind the pull bottom line is the second rear target area, and the distance between the partition and the pull bottom line can be directly determined as the influence range behind the second pull bottom lineD b (ii) a If the fourth percentage of the plurality of back partitions of the pull bottom line is larger than the second threshold value, the plurality of back partitions of the pull bottom line are determined as a second back target area, and then the maximum distance between the plurality of back partitions of the pull bottom line and the pull bottom line is determined as the influence range behind the second pull bottom lineD b . For example, referring to fig. 2, for the area behind the pull-bottom line, the rear partition of each pull-bottom line is B1, B2, B3, B4 and B5, and the distance from each partition to the reference point is 15m, 30m, 45m, 60m and 75m respectively; if the fourth percentage of the backfire rear subareas B1 and B2 is greater than the second threshold, the second rear target areas are the backfire rear subareas B1 and B2, and the second backfire rear influence rangeD b I.e., the distance of the rear draw line section B2 from the reference point is 30 m.
And 5: and determining the maximum value in the first bottom pulling wire front influence range and the second bottom pulling wire front influence range as the bottom pulling wire front influence range, and determining the maximum value in the first bottom pulling wire rear influence range and the second bottom pulling wire rear influence range as the bottom pulling wire rear influence range.
In one embodiment, the first pull-bottom line front influence range may beL f And the front influence range of the second pull bottom lineD f The maximum value of the two is determined as the influence range of the bottom pulling period bottom pulling operation on the front part of the bottom pulling wire in front of the bottom pulling wireR f1The influence range behind the first pull bottom lineL b Influence range behind the second pull bottom lineD b The maximum value of the two is determined as the influence range of the bottoming cycle bottoming operation on the rear part of the bottoming lineR b1。
Further, can be takenlThe average value of the influence ranges of the continuous bottoming operations is used as the influence range of the final bottoming operation (namely the influence range of the target in front of the bottoming operation and the influence range of the target behind the bottoming operation), and the number of the continuous bottoming periods is not less than 2. Specifically, the average value of the influence ranges in front of the bottom pulling lines in a plurality of bottom pulling periods is determined as the influence range of a target in front of the bottom pulling operation; and determining the average value of the influence ranges behind the bottom pulling lines of the plurality of bottom pulling periods as the target influence range behind the bottom pulling operation.
In a specific application, the following formula can be used for calculation:
Figure M_220215170831297_297174001
Figure M_220215170831375_375291001
wherein,R f showing the influence range of the target in front of the bottom pulling operation,R b and s represents the sequence number of the bottom pulling period.
For convenience of understanding, the embodiment of the present invention further provides a specific example of a method for determining an influence range of a bottoming operation, which mainly includes the following steps (1) to (7):
step (1): and (3) acquiring a vibration signal in a certain natural caving bottom pulling period by using a micro-seismic monitoring system, eliminating bottom pulling blasting signals and noise signals, processing micro-seismic events, and determining the seismic source positions and the seismic levels of all the micro-seismic events.
Step (2): determining the midpoint position of the bottom pulling line at the end of the bottom pulling period as a reference point, and taking the area within 15m of the distance between the two sides of the reference point along the bottom pulling tendency as an tendency analysis range; according to the position distribution of the periodic microseismic event, the front and the back of the bottom line are divided into regions according to the distance 15m by taking the reference point as the starting point in the inclination analysis range, the front of the bottom line is divided into 6 subareas, and the back of the bottom line is divided into 5 subareas, which can be seen in fig. 2.
And (3): the distances from the reference points to the sections F1, F2, F3, F4, F5 and F6 in front of the pull-down line are respectively 15m, 30m, 45m, 60m, 75m and 90m, and the distances from the sections B1, B2, B3, B4 and B5 in back of the pull-down line are respectively 15m, 30m, 45m, 60m and 75 m.
And (4): calculating to obtain the total number of all microseismic events in the tendency rangeQ278 microseismic events per zone, as shown in table 1, are calculated to obtain the percentage of the number of microseismic events in each zone to the total number of microseismic events.λTaking 10%, calculating to obtain the influence range of the front part of the periodic bottom line to be 30m and the influence range of the rear part of the periodic bottom line to be 30 m.
TABLE 1 number and percentage of microseismic events for each zone
Partition name F1 F2 F3 F4 F5 F6 B1 B2 B3 B4 B5
qk 64 48 26 19 14 3 47 29 8 13 5
δ 23% 17% 9% 7% 5% 1% 17% 10% 3% 5% 2%
And (5): calculating to obtain the seismic median of all the microseismic events in the trend analysis range as-2.9, and as shown in Table 2, obtaining the number of the microseismic events with the seismic magnitude larger than the seismic median in each subarea according to the seismic magnitude distribution of the microseismic events in each subarea, andthe number of microseismic events having a magnitude greater than the median magnitude in each zone is a percentage of the total number of microseismic events.
Figure SYM_220215170831001
And 5% is taken, and the influence range in front of the periodically drawn bottom line is 45m, and the influence range in back of the periodically drawn bottom line is 30 m.
TABLE 2 number and percentage of microseismic events with magnitude greater than the median of magnitude for each zone
Partition name F1 F2 F3 F4 F5 F6 B1 B2 B3 B4 B5
ek 33 25 13 10 9 1 24 16 5 9 3
ζ 12% 9% 5% 4% 3% 0% 9% 6% 2% 3% 1%
And (6): and (5) determining that the influence range in front of the periodic bottom drawing line is 45m and the influence range in back of the periodic bottom drawing line is 30m according to the statistical results of the step (4) and the step (5).
And (7): the average value of the influence ranges of the bottom pulling operation in 3 continuous cycles is taken as the influence range of the final bottom pulling operation, which is shown in table 3, namely the influence range of the front part of the bottom pulling line is 50m, and the influence range of the rear part of the bottom pulling line is 40 m.
TABLE 3 average of the impact Range of the bottoming operation for 3 consecutive periods
Cycle of bottom pulling operation T1 T2 T3 Mean value of
Rfs 45 45 60 50
Rfs 30 45 45 40
In summary, the method for determining the influence range of the natural caving method bottoming operation based on the real-time monitoring of the microseismic signal by the microseismic monitoring system provided by the embodiment of the invention has the following technical effects:
(1) the acquired microseismic signals can be accurately positioned, the time-space evolution and energy change of microseismic events can reflect the deformation and fracture state of the ore rock in real time, the actual mine ground pressure appearance rule on the site of the mine can be accurately described, and the analysis result can be used for monitoring and early warning of mine ground pressure disasters.
(2) The mining blasting signals collected by the micro-seismic monitoring system are screened and removed, the interference of blasting vibration on the total energy of the regional vibration signals is removed, the calculation of regional energy synthesis is replaced by counting the magnitude of each event, the major-magnitude micro-seismic events in each region are analyzed in an important mode, the interference of the major-magnitude blasting events in an energy calculation method on the total energy of the regions is avoided, and the stability evaluation of surrounding rocks of each region is more accurate.
(3) The invention aims at the influence range of the bottom pulling operation in the process of the bottom pulling operation by the natural caving method, gives consideration to the surrounding rock structures in front of and behind the working face, carries out all-around monitoring and analysis on the stability of the rock mass, takes the average value of the influence ranges of the bottom pulling operation for a plurality of continuous periods as the influence range of the final bottom pulling operation, and avoids the accidental division of the influence ranges.
(4) The invention mainly considers the earthquake disaster precursor of the major-earthquake-level microseismic event, and uses the ratio of the number of the microseismic events in the subarea to the major-earthquake-level microseismic event as a judgment basis, so that the influence range division of the bottoming operation is more accurate.
It should be noted that any particular value in all examples shown and described herein should be construed as merely exemplary and not limiting, and thus other examples of exemplary embodiments may have different values.
As to the method for determining the influence range of the backing-up operation provided in the foregoing embodiment, an embodiment of the present invention further provides a device for determining the influence range of the backing-up operation, and referring to a schematic structural diagram of the device for determining the influence range of the backing-up operation shown in fig. 3, the device may include the following components:
and the microseismic event determining module 301 is configured to, for each bottoming cycle, acquire a vibration signal in the bottoming cycle through the microseismic monitoring system, and determine a seismic source position and a seismic magnitude of the microseismic event in each bottoming cycle.
A partition dividing module 302, configured to divide the bottoming operation plane into multiple partitions based on the source position of the microseismic event, and calculate a distance between each partition and a bottoming line; wherein, the subregion includes: at least one pull-bottom line front partition and at least one pull-bottom line rear partition.
The first determining module 303 is configured to determine the number of microseismic events in each partition, and determine the pull-down operation influence range of the pull-down period based on the number of microseismic events in each partition and the magnitude of the microseismic events.
The second determining module 304 is configured to determine an average value of the pull-down operation influence ranges of the plurality of pull-down cycles as a pull-down operation target influence range.
The device for determining the influence range of the bottom pulling operation provided by the embodiment of the invention can perform region division according to the seismic source position of the microseismic event, and can simultaneously determine the influence ranges of the bottom pulling operation of the natural caving method on the front and the rear of the bottom pulling line, thereby accurately describing the actual ground pressure display rule of a mine site; secondly, the device determines the influence range of the bottom-pulling operation by the natural caving method based on the number of the microseismic events and the distribution of the seismic magnitude, and eliminates the interference of blasting vibration on the total energy of the regional vibration signals, so that the division of the influence range of the bottom-pulling operation is more accurate.
In an embodiment, the partitioning module 302 is further configured to: determining areas at two sides of a preset reference point and a first preset distance from the reference point along the bottom pulling tendency as an inclination range; wherein, the datum point is the midpoint of the bottom pulling line when the bottom pulling period is finished; based on the seismic source position of the microseismic event, in the inclination range, the reference point is used as a starting point, the front area of the bottom pulling line is divided into at least one front subarea of the bottom pulling line according to a preset distance, and the rear area of the bottom pulling line is divided into at least one rear subarea of the bottom pulling line.
In one embodiment, the pull-down operation impact range includes: the influence range in front of the bottom drawing line and the influence range behind the bottom drawing line; the first determining module 303 is further configured to: determining the number of microseismic events in each partition and the total number of microseismic events in the tendency range; calculating the ratio of the number of the microseismic events in each subarea to the number of the first microseismic events in the total number of the microseismic events, and determining the front influence range of the first bottoming line and the rear influence range of the first bottoming line based on the ratio of the number of the first microseismic events in each subarea; determining the number of microseismic events with the magnitude greater than a magnitude threshold in each zone based on the magnitude of the microseismic events; wherein, the magnitude threshold is the median of the magnitude of the microseismic event in the tendency range; calculating the second microseismic event number ratio of the number of microseismic events with the seismic level larger than the seismic level threshold value in each subarea to the total number of the microseismic events, and determining the front influence range of the second bottom pulling line and the rear influence range of the second bottom pulling line based on the second microseismic event number ratio in each subarea; and determining the maximum value in the first bottom pulling wire front influence range and the second bottom pulling wire front influence range as the bottom pulling wire front influence range, and determining the maximum value in the first bottom pulling wire rear influence range and the second bottom pulling wire rear influence range as the bottom pulling wire rear influence range.
In an embodiment, the first determining module 303 is further configured to: calculating a first percentage of the number of microseismic events of the partition in front of each bedding line in the total number of microseismic events and a second percentage of the number of microseismic events of the partition in back of each bedding line in the total number of microseismic events; determining the front partition of the bottom pulling line with the first percentage larger than a first threshold value as a first front target area, and determining the maximum distance between the first front target area and the bottom pulling line as the front influence range of the first bottom pulling line; and determining the rear partition of the bottom pulling line with the second percentage larger than the first threshold as a first rear target area, and determining the maximum distance between the first rear target area and the bottom pulling line as the rear influence range of the first bottom pulling line.
In an embodiment, the first determining module 303 is further configured to: calculating a third percentage of the number of microseismic events with the seismic level larger than the seismic level threshold value in the partition in front of each bottom pulling line in the total number of the microseismic events and a fourth percentage of the number of the microseismic events with the seismic level larger than the seismic level threshold value in the partition behind each bottom pulling line in the total number of the microseismic events; determining the front partition of the bottom pulling line with the third percentage larger than the second threshold as a second front target area, and determining the maximum distance between the second front target area and the bottom pulling line as the front influence range of the second bottom pulling line; and determining the rear partition of the bottom pulling line with the fourth percentage larger than the second threshold as a second rear target area, and determining the maximum distance between the second rear target area and the bottom pulling line as the front influence range of the second bottom pulling line.
In an embodiment, the second determining module 304 is further configured to: determining the average value of the front influence ranges of the bottom pulling lines in a plurality of bottom pulling periods as the front target influence range of the bottom pulling operation; and determining the average value of the influence ranges behind the bottom pulling lines of the plurality of bottom pulling periods as the target influence range behind the bottom pulling operation.
In one embodiment, the microseismic event determination module 301 is further configured to: and eliminating blasting signals and noise signals in the vibration signals.
The device provided by the embodiment of the present invention has the same implementation principle and technical effect as the method embodiments, and for the sake of brief description, reference may be made to the corresponding contents in the method embodiments without reference to the device embodiments.
The embodiment of the invention also provides electronic equipment, which specifically comprises a processor and a storage device; the storage means has stored thereon a computer program which, when executed by the processor, performs the method of any of the above embodiments.
Fig. 4 is a schematic structural diagram of an electronic device 100 according to an embodiment of the present invention, where the electronic device 100 includes: a processor 40, a memory 41, a bus 42 and a communication interface 43, wherein the processor 40, the communication interface 43 and the memory 41 are connected through the bus 42; the processor 40 is arranged to execute executable modules, such as computer programs, stored in the memory 41.
The Memory 41 may include a high-speed Random Access Memory (RAM) and may also include a non-volatile Memory (non-volatile Memory), such as at least one disk Memory. The communication connection between the network element of the system and at least one other network element is realized through at least one communication interface 43 (which may be wired or wireless), and the internet, a wide area network, a local network, a metropolitan area network, etc. may be used.
The bus 42 may be an ISA bus, PCI bus, EISA bus, or the like. The bus may be divided into an address bus, a data bus, a control bus, etc. For ease of illustration, only one double-headed arrow is shown in FIG. 4, but that does not indicate only one bus or one type of bus.
The memory 41 is used for storing a program, the processor 40 executes the program after receiving an execution instruction, and the method executed by the apparatus defined by the flow process disclosed in any of the foregoing embodiments of the present invention may be applied to the processor 40, or implemented by the processor 40.
The processor 40 may be an integrated circuit chip having signal processing capabilities. In implementation, the steps of the above method may be performed by integrated logic circuits of hardware or instructions in the form of software in the processor 40. The Processor 40 may be a general-purpose Processor, and includes a Central Processing Unit (CPU), a Network Processor (NP), and the like; the device can also be a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field-Programmable Gate Array (FPGA), or other Programmable logic devices, discrete Gate or transistor logic devices, discrete hardware components. The various methods, steps and logic blocks disclosed in the embodiments of the present invention may be implemented or performed. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like. The steps of the method disclosed in connection with the embodiments of the present invention may be directly implemented by a hardware decoding processor, or implemented by a combination of hardware and software modules in the decoding processor. The software module may be located in ram, flash memory, rom, prom, or eprom, registers, etc. storage media as is well known in the art. The storage medium is located in a memory 41, and the processor 40 reads the information in the memory 41 and completes the steps of the method in combination with the hardware thereof.
The computer program product of the readable storage medium provided in the embodiment of the present invention includes a computer readable storage medium storing a program code, where instructions included in the program code may be used to execute the method described in the foregoing method embodiment, and specific implementation may refer to the foregoing method embodiment, which is not described herein again.
The functions, if implemented in the form of software functional units and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
Finally, it should be noted that: the above-mentioned embodiments are only specific embodiments of the present invention, which are used for illustrating the technical solutions of the present invention and not for limiting the same, and the protection scope of the present invention is not limited thereto, although the present invention is described in detail with reference to the foregoing embodiments, those skilled in the art should understand that: any person skilled in the art can modify or easily conceive the technical solutions described in the foregoing embodiments or equivalent substitutes for some technical features within the technical scope of the present disclosure; such modifications, changes or substitutions do not depart from the spirit and scope of the embodiments of the present invention, and they should be construed as being included therein. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims (10)

1. A method for determining the influence range of a bottom pulling operation is characterized by comprising the following steps:
for each bottoming period, acquiring a vibration signal in the bottoming period through a microseismic monitoring system, and determining the seismic source position and the seismic magnitude of a microseismic event in each bottoming period;
dividing a bottoming operation surface into a plurality of partitions based on the seismic source position of the microseismic event, and calculating the distance between each partition and a bottoming line; wherein the partitioning comprises: at least one pull-bottom line front partition and at least one pull-bottom line rear partition;
determining the number of microseismic events in each partition, and determining the pull-down operation influence range of the pull-down period based on the number of microseismic events in each partition and the magnitude of the microseismic events;
and determining the average value of the influence ranges of the bottom pulling operation in the plurality of bottom pulling periods as a target influence range of the bottom pulling operation.
2. The method of claim 1, wherein partitioning the bottoming plane into a plurality of partitions based on the source location of the microseismic event comprises:
determining areas at two sides of a preset reference point and a first preset distance from the reference point along the bottom pulling tendency as an inclination range; wherein, the reference point is the midpoint of the bottom pulling line at the end of the bottom pulling period;
based on the source position of the microseismic event, in the inclination range, with the reference point as a starting point, dividing the front area of the bottom line into at least one front partition of the bottom line according to a preset distance, and dividing the rear area of the bottom line into at least one rear partition of the bottom line.
3. The method of claim 2, wherein the bottoming operation impact range comprises: the influence range in front of the bottom drawing line and the influence range behind the bottom drawing line;
determining the number of microseismic events in each partition, and determining the pull-down operation influence range of the pull-down period based on the number of microseismic events in each partition and the magnitude of the microseismic events, wherein the method comprises the following steps:
determining the number of microseismic events in each of the zones and the total number of microseismic events in the dip range;
calculating the ratio of the number of the microseismic events in each partition to the number of the first microseismic events in the total number of the microseismic events, and determining the front influence range of a first bottoming line and the rear influence range of the first bottoming line based on the ratio of the number of the first microseismic events in each partition;
determining, based on the magnitude of the microseismic events, a number of microseismic events for which the magnitude in each of the zones is greater than a magnitude threshold; wherein the magnitude threshold is the median of magnitudes of microseismic events within the trend range;
calculating the second microseismic event number ratio of the number of microseismic events with the seismic magnitude larger than the seismic magnitude threshold value in each subarea to the total number of the microseismic events, and determining the front influence range of a second bottoming line and the rear influence range of the second bottoming line based on the second microseismic event number ratio in each subarea;
determining the maximum value in the first bottom pulling line front influence range and the second bottom pulling line front influence range as a bottom pulling line front influence range, and determining the maximum value in the first bottom pulling line rear influence range and the second bottom pulling line rear influence range as a bottom pulling line rear influence range.
4. A method according to claim 3 wherein calculating a first ratio of the number of microseismic events in each of the zones to the total number of microseismic events and determining a first pull-bottom line front impact range and a first pull-bottom line back impact range based on the first ratio of the number of microseismic events in each of the zones comprises:
calculating a first percentage of the number of microseismic events of each bottom wire front section to the total number of microseismic events and a second percentage of the number of microseismic events of each bottom wire back section to the total number of microseismic events;
determining the front partition of the bottom pulling line with the first percentage larger than a first threshold value as a first front target area, and determining the maximum distance between the first front target area and the bottom pulling line as a first front influence range of the bottom pulling line;
and determining the rear partition of the bottom pulling line with the second percentage larger than the first threshold value as a first rear target area, and determining the maximum distance between the first rear target area and the bottom pulling line as a first bottom pulling line rear influence range.
5. A method according to claim 3 wherein calculating a second ratio of the number of microseismic events having a magnitude greater than a magnitude threshold in each of the zones to the total number of microseismic events and determining a second bedding line front impact range and a second bedding line rear impact range based on the second ratio of the number of microseismic events in each of the zones comprises:
calculating a third percentage of the number of microseismic events with the seismic level larger than the seismic level threshold value in each bottom line front partition to the total number of the microseismic events and a fourth percentage of the number of the microseismic events with the seismic level larger than the seismic level threshold value in each bottom line rear partition to the total number of the microseismic events;
determining the front partition of the bottom pulling line with the third percentage larger than a second threshold value as a second front target area, and determining the maximum distance between the second front target area and the bottom pulling line as a second front influence range of the bottom pulling line;
and determining the rear partition of the bottom pulling line with the fourth percentage larger than the second threshold as a second rear target area, and determining the maximum distance between the second rear target area and the bottom pulling line as a front influence range of the second bottom pulling line.
6. The method of claim 3, wherein determining an average of pull-down operation target influence ranges of a plurality of pull-down cycles comprises:
determining the average value of the front influence ranges of the bottom pulling lines of the plurality of bottom pulling periods as the front target influence range of the bottom pulling operation;
and determining the average value of the influence ranges behind the bottoming lines of the plurality of bottoming cycles as the target influence range behind the bottoming operation.
7. The method of claim 1, wherein prior to determining the source location and magnitude of the microseismic event within each of the pull-down periods, the method further comprises:
and eliminating blasting signals and noise signals in the vibration signals.
8. An apparatus for determining the range of influence of a bottoming operation, comprising:
the microseismic event determining module is used for acquiring a vibration signal in each bottoming period through a microseismic monitoring system and determining the seismic source position and the seismic magnitude of a microseismic event in each bottoming period;
the partition dividing module is used for dividing the bottoming operation surface into a plurality of partitions based on the seismic source position of the microseismic event and calculating the distance between each partition and a bottoming line; wherein the partitioning comprises: at least one pull-bottom line front partition and at least one pull-bottom line rear partition;
the first determination module is used for determining the number of microseismic events in each partition, and determining the influence range of the bottom pulling operation of the bottom pulling period based on the number of the microseismic events in each partition and the magnitude of the microseismic events;
and the second determining module is used for determining the average value of the influence ranges of the bottom pulling operation in the plurality of bottom pulling periods as the target influence range of the bottom pulling operation.
9. An electronic device comprising a processor and a memory, the memory storing computer-executable instructions executable by the processor, the processor executing the computer-executable instructions to perform the steps of the method of any one of claims 1 to 7.
10. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the method according to any one of the claims 1 to 7.
CN202210143734.XA 2022-02-17 2022-02-17 Method and device for determining influence range of bottoming operation, electronic equipment and medium Active CN114185084B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210143734.XA CN114185084B (en) 2022-02-17 2022-02-17 Method and device for determining influence range of bottoming operation, electronic equipment and medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210143734.XA CN114185084B (en) 2022-02-17 2022-02-17 Method and device for determining influence range of bottoming operation, electronic equipment and medium

Publications (2)

Publication Number Publication Date
CN114185084A true CN114185084A (en) 2022-03-15
CN114185084B CN114185084B (en) 2022-05-17

Family

ID=80546109

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210143734.XA Active CN114185084B (en) 2022-02-17 2022-02-17 Method and device for determining influence range of bottoming operation, electronic equipment and medium

Country Status (1)

Country Link
CN (1) CN114185084B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008124759A1 (en) * 2007-04-09 2008-10-16 Baker Hughes Incorporated Location of a microseismic event using the combination of a vsp array and a surface array
CN104297449A (en) * 2014-10-20 2015-01-21 北京矿冶研究总院 Method for testing similar model of cemented jointed rock mass by natural caving method
CN106777772A (en) * 2017-01-09 2017-05-31 辽宁工程技术大学 A kind of Rock Burst hazard prediction method based on coal or rock dynamic system
CN108798785A (en) * 2018-06-06 2018-11-13 中煤能源研究院有限责任公司 A kind of coal mine rock burst prediction and warning method
CN113338924A (en) * 2021-05-11 2021-09-03 紫金矿业集团股份有限公司 Control method for surface subsidence range by natural caving method
CN113657044A (en) * 2021-07-30 2021-11-16 深圳市中金岭南有色金属股份有限公司凡口铅锌矿 Ore body large-seismic-level event early warning method and device, terminal equipment and medium

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008124759A1 (en) * 2007-04-09 2008-10-16 Baker Hughes Incorporated Location of a microseismic event using the combination of a vsp array and a surface array
CN104297449A (en) * 2014-10-20 2015-01-21 北京矿冶研究总院 Method for testing similar model of cemented jointed rock mass by natural caving method
CN106777772A (en) * 2017-01-09 2017-05-31 辽宁工程技术大学 A kind of Rock Burst hazard prediction method based on coal or rock dynamic system
CN108798785A (en) * 2018-06-06 2018-11-13 中煤能源研究院有限责任公司 A kind of coal mine rock burst prediction and warning method
CN113338924A (en) * 2021-05-11 2021-09-03 紫金矿业集团股份有限公司 Control method for surface subsidence range by natural caving method
CN113657044A (en) * 2021-07-30 2021-11-16 深圳市中金岭南有色金属股份有限公司凡口铅锌矿 Ore body large-seismic-level event early warning method and device, terminal equipment and medium

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
彭张 等: "自然崩落法放矿量与底部结构应力分布关系研究", 《有色金属（矿山部分）》 *
杨云琛 等: "普朗铜矿地压综合监测数据分析与研究*", 《采矿技术》 *

Also Published As

Publication number Publication date
CN114185084B (en) 2022-05-17

Similar Documents

Publication Publication Date Title
RU2605192C2 (en) Updating microseismic histogram data
CN109271738B (en) Numerical inversion method for acquiring Weibull distribution parameters of roadway surrounding rock
EP2924471B1 (en) Systems and methods for modeling fracture networks in reservoir volumes from microseismic events
Yarahmadi et al. Discontinuity modelling and rock block geometry identification to optimize production in dimension stone quarries
CN115343757A (en) Ground pressure risk early warning method and device, electronic equipment and medium
CN115390130B (en) Method and device for predicting large-energy microseismic event in coal mining
CN117250670B (en) Detection capability assessment method, system and equipment of mine earthquake integrated monitoring station network
CN114814940B (en) Microseism monitoring station network evaluation method and device, electronic equipment and medium
CN113189644B (en) Microseismic source positioning method and system
CN114185084B (en) Method and device for determining influence range of bottoming operation, electronic equipment and medium
US20150355354A1 (en) Method of analyzing seismic data
Guzev et al. Algorithm of potentially burst-hazard zones dynamics representation in massif of rocks by results of seismic-acoustic monitoring
CN116166988A (en) Coal-bearing stratum abnormal structure identification and classification method
CN115688428A (en) Method, device and equipment for calculating critical value of explosive quantity of blasting near zone supporting structure
CN114879256B (en) Rock burst monitoring method and device, medium and equipment
CN114723863B (en) Vibration speed cloud picture prediction method and system for open bench blasting and electronic equipment
CN115310686A (en) Method, device, equipment and medium for predicting blasting qualified rate of strip mine
Palleske A study of biases, assumptions and practical consideration for the use of discrete fracture networks in geomechanical practice
CN113919500A (en) Microseism event space-time fractal dimension calculation method and device and electronic equipment
CN110824559B (en) Method and device for generating co-coordinate point gather and computer-readable storage medium
Kluckner et al. Estimation of the in situ block size in jointed rock masses using three-dimensional block simulations and discontinuity measurements
US20180306939A1 (en) Correcting Biases In Microseismic-Event Data
CN114518595B (en) Method and device for positioning micro-seismic source, electronic equipment and medium
CN112305603B (en) Quality monitoring method and system for node instrument
RU2815465C1 (en) System for semi-deterministic modeling of fracturing based on random matrix

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