Method and system for positioning surrounding rock stress anomaly area in mining disturbance area

Abstract

The invention provides a method and a system for positioning a surrounding rock stress anomaly area in a mining disturbance area, and relates to the technical field of stress anomaly identification and localization. Data of a primary rock stress field is acquired, a ground stress value of a to-be-measured point of the primary rock stress field is calculated, and the stress anomaly is determined based on known measuring point data and the calculated stress value. The method comprises the following steps: carrying out gridding partition on a roadway surrounding rock area, calculating the maximum horizontal stress, the minimum horizontal stress and the vertical stress of a to-be-measured area by utilizing an interpolation method, comprehensively generating primary rock stress, collecting surrounding rock three-direction mining-induced stress data, calculating axial stress, circumferential stress and radial stress, generating a stress threshold value, and comparing with the primary rock stress. And positions and distribution of abnormal areas can be rapidly and accurately judged. Through accurate stress calculation, flexible sensor layout and a comprehensive analysis method, more effective mine disaster risk early warning and forecasting and more reliable mine safety management are realized.

Description

Technical Field

[0001] The present invention relates to the technical field of stress anomaly identification and positioning, and in particular to a method and system for positioning an abnormal stress area of ​​surrounding rocks in a mining disturbance zone. Background Art

[0002] Mining will cause redistribution of the geostress field, leading to deformation and damage of the surrounding rock. Therefore, identifying and monitoring areas of abnormal surrounding rock stress is crucial to preventing catastrophic events. The traditional method is to use surface deformation monitoring instruments, such as total stations, GPS, etc., to monitor the deformation of the surface of the mining area.

[0003] In the prior art, the publication number CN108802825B discloses a method and system for locating coal-rock dynamic disasters by infrasound monitoring, which is connected to the monitoring main station through an optical fiber network; the monitoring substation is composed of three infrasound sensors arranged in a triangle, and the infrasound sensors are connected to the infrasound monitoring instrument; the monitoring main station is equipped with a computing server; the positioning method is that the main station computing server filters the collected data, performs time-frequency analysis on the data through short-time Fourier transform (STFT), compares the frequency energy density, extracts the main frequency band signal, and locates the calculation result through statistical analysis and time delay estimation theory. The system monitors the infrasound signal generated by the fracture of loaded coal and rock, and uses this method to perform non-contact and regional testing on the stress distribution state of the surrounding rock of the mine or tunnel to determine the coal and rock fracture and stress abnormality areas.

[0004] Insufficient existing technology:

[0005] In the existing technology, although the infrasound monitoring system is capable of conducting non-contact regional testing, the data needs to undergo multiple processing steps, including filtering, short-time Fourier transform and statistical analysis, which will lead to insufficient real-time performance. In addition, the complex underground environment may affect the accuracy of the signal, especially when dealing with noise and outliers, which will affect the positioning accuracy, resulting in insufficient real-time and accuracy of data processing.

[0006] Existing technologies mainly rely on short-time Fourier transform and time delay estimation theory for signal processing and positioning calculations, but these algorithms have limited adaptability and robustness when dealing with complex and dynamically changing underground environments. When faced with strong noise interference, multipath effects or reflected signals, STFT may not be able to accurately extract the main frequency band signals, affecting the reliability of the positioning results. It is necessary to introduce more advanced signal processing algorithms to improve the performance of the system under different working conditions.

[0007] Therefore, it is necessary to provide a method and system for locating abnormal stress areas of surrounding rocks in mining disturbance areas to solve the above problems.

[0008] The above information disclosed in this Background section is only for enhancement of understanding of the background of the present disclosure and therefore it may contain information that does not constitute the prior art that is already known to one of ordinary skill in the art. Summary of the invention

[0009] The purpose of the present invention is to provide a method and system for locating abnormal stress areas of surrounding rocks in a mining disturbance zone, so as to solve the problems raised in the above-mentioned background technology.

[0010] To achieve the above object, the present invention provides the following technical solutions:

[0011] A method for locating an abnormal area of ​​surrounding rock stress in a mining disturbance area, the specific steps comprising:

[0012] Step 1: Based on the three-dimensional grid division, the tunnel surrounding rock area is divided into multiple test areas, and multiple test points are set in the test areas. The maximum horizontal stress, minimum horizontal stress and vertical force of the test points are determined based on the stress relief method;

[0013] Step 2: Calculate the maximum horizontal stress, minimum horizontal stress and vertical stress of the area to be tested based on the maximum horizontal stress, minimum horizontal stress and vertical stress of each point to be tested in the area to be tested;

[0014] Step 3: According to the calculation, the maximum horizontal stress, the minimum horizontal stress and the vertical stress of the area to be measured are obtained, and the original rock stress of the area to be measured is comprehensively generated;

[0015] Step 4: Bury a high-precision, long-lasting three-way mining stress monitoring sensor at the center of the area to be tested to collect its mining stress, and decompose the mining stress to generate axial stress, annular stress and radial stress;

[0016] Step 5: Based on the axial stress, annular stress and radial stress, a stress threshold of the area to be tested is comprehensively generated, and the stress threshold of the same area to be tested is compared with the original rock stress. If the original rock stress exceeds the stress threshold, it indicates that the area to be tested is an abnormal area.

[0017] Furthermore, the maximum horizontal stress, the minimum horizontal stress and the vertical stress of the test point are obtained according to the following method:

[0018] Using the hydraulic fracturing method, water is injected into the borehole to increase the pressure in the hole until the rock breaks. The fracture pressure is related to the ground stress. The initial fracture pressure, pore pressure and closing pressure of the rock at the test point are recorded to calculate the maximum horizontal stress, minimum horizontal stress and vertical stress. The formula is:

[0019]

[0020] σv =γ*h*g

[0021] Among them, σ H , σ h , σ v Respectively represent the maximum horizontal stress, minimum horizontal stress, and vertical stress of the test point, P hf Indicates the initial fracture pressure of rock, P c Indicates closing pressure, P k represents the pore pressure, γ represents the density of the surrounding rock, h represents the drilling depth, and g represents the deadweight of the rock.

[0022] Furthermore, the maximum horizontal stress, the minimum horizontal stress and the vertical stress at the measured area are obtained by:

[0023] Determine the coordinates of the center point of the area to be measured and mark it as (x j ,y j , z j ), calibrate the coordinates of the i-th point to be measured in the area to be measured as (x i ,y i , z i ), i represents the index of the test point in the test area, and i∈[1,n], n is the number of test points in the test area, based on the stress data of each test point in the test area and the distance between the test point and the center point, the maximum horizontal stress, minimum horizontal stress and vertical stress of the center point are calculated, and used as the maximum horizontal stress, minimum horizontal stress and vertical stress of the test area, and the maximum horizontal stress, minimum horizontal stress and vertical stress of the center point are calculated, and the formula based on which is:

[0024]

[0025] in, Respectively represent the center point (x j ,y j , z j ), the maximum horizontal stress, minimum horizontal stress and vertical stress at Respectively represent the points to be tested (x i ,y i , z i ), the maximum horizontal stress, minimum horizontal stress and vertical stress at Represents the distance from the center point to the i-th point to be measured.

[0026] Furthermore, the original rock stress is comprehensively generated based on the formula:

[0027]

[0028] Among them, σ 原represents the original rock stress in the measured area, ω H ,ω h ,ω v are the weight ratios corresponding to the maximum horizontal stress, minimum horizontal stress and vertical stress, ω H >ω h >ω v , and satisfies ω H +ω h +ω v =1.

[0029] Furthermore, the axial stress, the hoop stress and the radial stress are obtained according to the following method:

[0030] Polar coordinate method is used to decompose the dynamic stress of the center point into the dynamic stress perpendicular to the yz plane in the x direction, the dynamic stress perpendicular to the xz plane in the y axis direction, and the dynamic stress perpendicular to the xy plane in the z axis direction. The axial stress is along the axis direction and is directly expressed by the measured dynamic stress perpendicular to the xy plane. The calculation of the hoop stress and radial stress requires determining the position angle and integrating the shear stress of the xy plane. The formula is:

[0031] F z =σ zz

[0032] F θ =σ xx *cos 2 (θ)+σ yy *sin 2 (θ)+2τ xy *cos(θ)*sin(θ)

[0033] F r =σ xx *sin 2 (θ)+σ yy *cos 2 (θ)-2τ xy *cos(θ)*sin(θ)

[0034] Among them, F z 、F θ 、F r Respectively represent axial stress, hoop stress, radial stress, σ xx , σ yy , σ zz They represent the dynamic stress perpendicular to the yz plane in the x direction, the dynamic stress perpendicular to the xz plane in the y direction, and the dynamic stress perpendicular to the xy plane in the z direction, respectively. xy represents the shear stress in the xy plane, and θ represents the position angle in polar coordinates.

[0035] Furthermore, the method for generating the stress threshold is based on:

[0036] T=a*F z +b*F θ +c*F r

[0037] Where T represents the stress threshold, a, b, and c represent the weight ratios of the corresponding axial stress, hoop stress, and radial stress, respectively, and a <b=c。

[0038] Furthermore, compared with the generated original rock stress, if the original rock stress exceeds the stress threshold, it means that the area to be tested is an abnormal area. The logic is as follows:

[0039]

[0040] Among them, Q represents a logical value for judging whether the area to be tested is an abnormal area. When Q=0, it means that the original rock stress does not exceed the stress threshold, and part of the area to be tested is not an abnormal area; when Q=1, it means that the original rock stress exceeds the stress threshold, and part of the area to be tested is a stress abnormal area.

[0041] The present invention further provides a system for locating an abnormal area of ​​surrounding rock stress in a mining disturbance zone, wherein the abnormal area locating system is used to execute any one of the above-mentioned methods for locating an abnormal area of ​​surrounding rock stress in a mining disturbance zone, comprising:

[0042] A stress data acquisition module, wherein the stress data acquisition module divides the tunnel surrounding rock area into a plurality of test areas based on three-dimensional grid division, sets a plurality of test points in the test areas, and determines the maximum horizontal stress, the minimum horizontal stress, and the vertical stress of the test points based on a stress relief method;

[0043] A gridding and interpolation calculation module, which calculates the maximum horizontal stress, the minimum horizontal stress and the vertical stress of the area to be measured based on the maximum horizontal stress, the minimum horizontal stress and the vertical stress of each point to be measured in the area to be measured;

[0044] An original rock stress generation module, the original rock stress generation module is used to obtain the maximum horizontal stress, the minimum horizontal stress and the vertical stress of the area to be measured according to the calculation, and comprehensively generate the original rock stress of the area to be measured;

[0045] A mining stress acquisition module, wherein the mining stress acquisition module is used to bury a high-precision, long-term three-way mining stress monitoring sensor at the center point of the area to be measured to collect the mining stress, and decompose the mining stress to generate axial stress, annular stress and radial stress;

[0046] The stress threshold and anomaly detection module is used to comprehensively generate a stress threshold of the area to be tested based on the axial stress, annular stress and radial stress, and compare the stress threshold of the same area to be tested with the original rock stress. If the original rock stress exceeds the stress threshold, it means that the area to be tested is a stress abnormality area.

[0047] Compared with the prior art, the present invention has the following beneficial effects:

[0048] The present invention improves the accuracy of stress distribution calculation in the area to be measured by introducing grid partitioning and interpolation methods, thereby enhancing the real-time and accuracy of data processing. Through flexible sensor layout, by burying three-way mining stress monitoring sensors at the center of the grid, the method is applicable to various geological conditions and complex environments, providing a more flexible stress monitoring method, being able to more accurately monitor and identify coal and rock dynamic disasters, and reducing false alarms and missed alarms.

[0049] The present invention not only relies on infrasonic signals, but also combines a variety of stress data for comprehensive analysis. This multi-level stress analysis enhances the ability to identify abnormal areas and helps to more accurately identify potential risks. By introducing a stress threshold comparison method, the present invention can effectively determine whether the original rock stress exceeds the safety range, timely identify and warn abnormal areas, and improve the safety and reliability of the overall system.

[0050] The present invention greatly improves monitoring accuracy, adaptability and safety through precise stress calculation, flexible sensor layout and comprehensive analysis methods, and realizes effective risk management and more reliable mine safety management. BRIEF DESCRIPTION OF THE DRAWINGS

[0051] Figure 1 It is a schematic diagram of the overall method flow of the present invention.

[0052] Figure 2 It is a schematic diagram of the system module flow of the present invention. DETAILED DESCRIPTION

[0053] In order to make the objectives, technical solutions and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with specific embodiments.

[0054] It should be noted that, unless otherwise defined, the technical terms or scientific terms used in the present invention should be understood by people with ordinary skills in the field to which the present invention belongs. The words "first", "second" and similar words used in the present invention do not indicate any order, quantity or importance, but are only used to distinguish different components. "Include" or "comprise" and similar words mean that the elements or objects appearing before the word include the elements or objects listed after the word and their equivalents, without excluding other elements or objects. "Connect" or "connected" and similar words are not limited to physical or mechanical connections, but may include electrical connections, whether direct or indirect. "Up", "down", "left", "right" and the like are only used to indicate relative positional relationships. When the absolute position of the described object changes, the relative positional relationship may also change accordingly.

[0055] Example:

[0056] See also Figure 1 , a method for locating abnormal stress areas of surrounding rocks in a mining disturbance area, the specific steps comprising:

[0057] Step 1: Based on the three-dimensional grid division, the tunnel surrounding rock area is divided into multiple test areas, and multiple test points are set in the test areas. The maximum horizontal stress, minimum horizontal stress and vertical force of the test points are determined based on the stress relief method;

[0058] Step 2: Calculate the maximum horizontal stress, minimum horizontal stress and vertical stress of the area to be tested based on the maximum horizontal stress, minimum horizontal stress and vertical stress of each point to be tested in the area to be tested;

[0059] Step 3: According to the calculation, the maximum horizontal stress, the minimum horizontal stress and the vertical stress of the area to be measured are obtained, and the original rock stress of the area to be measured is comprehensively generated;

[0060] Step 4: Bury a high-precision, long-lasting three-way mining stress monitoring sensor at the center of the area to be tested to collect its mining stress, and decompose the mining stress to generate axial stress, annular stress and radial stress;

[0061] Step 5: Based on the axial stress, annular stress and radial stress, a stress threshold of the area to be tested is comprehensively generated, and the stress threshold of the same area to be tested is compared with the original rock stress. If the original rock stress exceeds the stress threshold, it indicates that the area to be tested is an abnormal area.

[0062] It should be noted that hydraulic fracturing provides direct data on ground stress by directly measuring the fracture pressure and closure pressure. The maximum and minimum horizontal stresses and vertical stresses are calculated using formulas, which can convert the measured data into usable information and incorporate factors such as initial fracture pressure, pore pressure, and closure pressure into the calculation, ensuring the comprehensiveness and accuracy of stress assessment.

[0063] Therefore, it is necessary to obtain the maximum horizontal stress, minimum horizontal stress and vertical stress of the test point, and the method is as follows:

[0064] Using the hydraulic fracturing method, water is injected into the borehole to increase the pressure in the hole until the rock breaks. The fracture pressure is related to the ground stress. The initial fracture pressure, pore pressure and closing pressure of the rock at the test point are recorded to calculate the maximum horizontal stress, minimum horizontal stress and vertical stress. The formula is:

[0065]

[0066] σ v =γ*h*g

[0067] Among them, σ H , σ h , σ v Respectively represent the maximum horizontal stress, minimum horizontal stress, and vertical stress of the test point, P hf Indicates the initial fracture pressure of rock, P c Indicates closing pressure, P k represents pore pressure, γ represents surrounding rock density, h represents drilling depth, and g represents rock deadweight; in the above formula, P c The closing pressure refers to the pressure when the pressure in the hole drops until it remains stable.

[0068] It should be noted that grid partitioning can divide complex three-dimensional space into uniform small units, making the calculation and analysis of stress distribution more detailed and accurate. Grid stress calculation for the entire area can predict the changing trend of ground stress at different locations, which is helpful to identify possible dangerous areas, prevent and deal with them. The gaps between known measuring point data are filled by interpolation methods to enhance the integrity and visualization of data. The grid unit size of 5 mm in length, width and height is selected to ensure the calculation accuracy while taking into account the operability and efficiency of the calculation. Too large or too small grid units will affect the calculation accuracy and efficiency. The use of three-dimensional grid generation algorithms can automatically divide complex three-dimensional space into regular grid units, which is helpful for the systematic and consistent processing of data. The weighting factor and distance calculation method used in the formula ensure the rationality and accuracy of the interpolation results. Through effective weighting and distance calculation, the distribution characteristics of stress in space can be more realistically reflected.

[0069] Therefore, the method for obtaining the maximum horizontal stress, the minimum horizontal stress and the vertical stress in the measured area is:

[0070] Determine the coordinates of the center point of the area to be measured and mark it as (x j ,y j , z j ), calibrate the coordinates of the i-th point to be measured in the area to be measured as (x i ,y i , z i ), i represents the index of the test point in the test area, and i∈[1,n], n is the number of test points in the test area, based on the stress data of each test point in the test area and the distance between the test point and the center point, the maximum horizontal stress, minimum horizontal stress and vertical stress of the center point are calculated, and used as the maximum horizontal stress, minimum horizontal stress and vertical stress of the test area, and the maximum horizontal stress, minimum horizontal stress and vertical stress of the center point are calculated, and the formula based on which is:

[0071]

[0072] in, Respectively represent the center point (x j ,y j , z j ), the maximum horizontal stress, minimum horizontal stress and vertical stress at Respectively represent the points to be tested (x i ,y i , z i ), the maximum horizontal stress, minimum horizontal stress and vertical stress at Represents the distance from the center point to the i-th point to be measured; in the above formula, the inverse distance weighted interpolation method is a commonly used spatial interpolation method. Its basic idea is that the closer the distance is, the greater the influence on the interpolation result. In the three formulas for calculating the stress value of the center point, the numerator is the weighted sum of the stress value of each measuring point by its distance from the center point. The closer the distance is, the greater the influence on the result; the denominator is the weighted inverse sum of the distances of all measuring points, which is used for normalization processing to make the result independent of the specific number and distance of measuring points.

[0073] It should be noted that integrating stresses in different directions into an original rock stress value can more comprehensively evaluate the overall stress state of the rock, and through the synthetic stress index, the analysis process can be simplified, which is convenient for engineering design and decision making. The setting of weights allows the reflection of the different degrees of influence of stresses in different directions on the rock mass, thereby providing a more accurate stress evaluation. Through reasonable weight distribution, the stress characteristics of the rock mass under actual geological conditions can be more realistically simulated, and the accuracy and reliability of the model prediction can be improved. In underground engineering, comprehensive stress evaluation can better guide support design and construction plans and improve safety.

[0074] Therefore, it is necessary to comprehensively generate the original rock stress, based on the formula:

[0075]

[0076] Among them, σ 原 represents the original rock stress in the measured area, ω H ,ω h ,ω v are the weight ratios corresponding to the maximum horizontal stress, minimum horizontal stress and vertical stress, ω H >ω h >ω v , and satisfies ω H +ω h +ω v =1; In the above formula, the maximum horizontal stress is given the highest proportional weight because in geological conditions, horizontal stress has a greater impact on the deformation and destruction of rock mass, especially in underground engineering; the minimum horizontal stress is given the second highest proportional weight because the minimum horizontal stress significantly affects the behavior of rock mass, but has a smaller impact on the deformation and destruction of rock mass than the maximum horizontal stress; and the vertical stress is given the smallest proportional weight because the vertical stress is generally generated by the deadweight of the stratum, and its impact on the rock mass is far less than the impact of the maximum and minimum horizontal stresses on the rock mass.

[0077] It should be noted that by calculating the axial, annular and radial stresses, the stress distribution during the mining process can be comprehensively analyzed, providing key data support for engineering design, identifying and evaluating stress concentration areas in different directions, and helping to predict potential danger points and ensure structural stability and safety. Moreover, through the polar coordinate method, complex dynamic stresses can be decomposed into more understandable components, facilitating real-time evaluation of the impact of dynamic changes on the structure.

[0078] Therefore, it is necessary to obtain the axial stress, hoop stress and radial stress, and the method is based on:

[0079] Polar coordinate method is used to decompose the dynamic stress of the center point into the dynamic stress perpendicular to the yz plane in the x direction, the dynamic stress perpendicular to the xz plane in the y axis direction, and the dynamic stress perpendicular to the xy plane in the z axis direction. The axial stress is along the axis direction and is directly expressed by the measured dynamic stress perpendicular to the xy plane. The calculation of the hoop stress and radial stress requires determining the position angle and integrating the shear stress of the xy plane. The formula is:

[0080] F z =σ zz

[0081] F θ =σ xx *cos 2 (θ)+σ yy *sin 2 (θ)+2τ xy *cos(θ)*sin(θ)

[0082] F r =σ xx *sin 2 (θ)+σ yy *cos 2 (θ)-2τ xy *cos(θ)*sin(θ)

[0083] Among them, F z 、F θ 、F r Respectively represent axial stress, hoop stress, radial stress, σ xx , σ yy , σ zz They represent the dynamic stress perpendicular to the yz plane in the x direction, the dynamic stress perpendicular to the xz plane in the y direction, and the dynamic stress perpendicular to the xy plane in the z direction, respectively. xy represents the shear stress in the xy plane, and θ represents the position angle in polar coordinates.

[0084] It should be noted that by synthesizing stresses in different directions to generate a stress threshold, the bearing capacity of the rock mass structure can be more comprehensively reflected. Calculating the threshold can help identify weak links in the mining area, thereby optimizing the design and improving the overall stability. The stress threshold provides an index for real-time monitoring and rapid assessment of stress changes, improving the response speed.

[0085] Therefore, the method for generating the stress threshold is as follows:

[0086] T = a*F z +b*F θ +c*F r

[0087] Where T represents the stress threshold, and a, b, and c respectively represent the weight ratios corresponding to the axial stress, circumferential stress, and radial stress, and a < b = c; in the above formula, the ratio weights are set as a < b = c because in various rock mass structures, the circumferential and radial stresses have a more significant impact on the bearing capacity and stability. Therefore, equal and relatively high weights are set. The role of the axial stress is often to affect the deformation of the structure in the axial direction, ensuring that it maintains the designed shape and size, and it has a less significant impact on the bearing capacity and stability compared to the other two stresses. So, the smallest weight ratio is given.

[0088] It should be noted that by establishing a logical judgment formula, potential abnormal areas can be effectively identified, preventive maintenance can be carried out, structural failure or accidents can be avoided, resources and energy can be concentrated on areas that need attention, the efficiency of detection and maintenance work can be improved, and the structure can be ensured to operate within a safe range, avoiding risks caused by exceeding the designed bearing capacity.

[0089] Therefore, the stress threshold needs to be compared with the generated in-situ stress. If the in-situ stress exceeds the stress threshold, it indicates that the area to be measured is an abnormal area. The basis for the logic is as follows:

[0090]

[0091] Where Q represents the logical value for judging whether the area to be measured is an abnormal area. When Q = 0, it indicates that the in-situ stress does not exceed the stress threshold, and the area to be measured is not an abnormal area; when Q = 1, it indicates that the in-situ stress exceeds the stress threshold, and the area to be measured is a stress abnormal area.

[0092] Please refer to Figure 2 , a system for locating stress abnormal areas in the surrounding rock of a mining disturbance area. The abnormal area location system is used to execute any one of the above-mentioned methods for locating stress abnormal areas in the surrounding rock of a mining disturbance area, and includes:

[0093] A stress data acquisition module, wherein the stress data acquisition module divides the tunnel surrounding rock area into a plurality of test areas based on three-dimensional grid division, sets a plurality of test points in the test areas, and determines the maximum horizontal stress, the minimum horizontal stress, and the vertical stress of the test points based on a stress relief method;

[0094] A gridding and interpolation calculation module, which calculates the maximum horizontal stress, the minimum horizontal stress and the vertical stress of the area to be measured based on the maximum horizontal stress, the minimum horizontal stress and the vertical stress of each point to be measured in the area to be measured;

[0095] An original rock stress generation module, the original rock stress generation module is used to obtain the maximum horizontal stress, the minimum horizontal stress and the vertical stress of the area to be measured according to the calculation, and comprehensively generate the original rock stress of the area to be measured;

[0096] A mining stress acquisition module, wherein the mining stress acquisition module is used to bury a high-precision, long-term three-way mining stress monitoring sensor at the center point of the area to be measured to collect the mining stress, and decompose the mining stress to generate axial stress, annular stress and radial stress;

[0097] The stress threshold and anomaly detection module is used to comprehensively generate a stress threshold of the area to be tested based on the axial stress, annular stress and radial stress, and compare the stress threshold of the same area to be tested with the original rock stress. If the original rock stress exceeds the stress threshold, it means that the area to be tested is a stress abnormality area.

[0098] The above formulas are all dimensionless and numerical calculations. The formula is a formula for the most recent real situation obtained by collecting a large amount of data and performing software simulation. The preset parameters in the formula are set by technicians in this field according to actual conditions.

[0099] The above embodiments may be implemented in whole or in part by software, hardware, firmware or any other combination thereof. When implemented by software, the above embodiments may be implemented in whole or in part in the form of a computer program product. Those skilled in the art may appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein may be implemented by electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are performed by hardware or software methods depends on the specific application and design constraints of the technical solution.

[0100] The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, and may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0101] The above description is only a specific implementation manner of the present application, but the protection scope of the present application is not limited thereto. Any technician familiar with the technical field can easily think of changes or substitutions within the technical scope disclosed in the present application, which should be included in the protection scope of the present application.

Claims

1. A method for locating abnormal stress areas of surrounding rocks in a mining disturbance area, characterized in that: The specific steps include: Step 1: Based on the three-dimensional grid division, the tunnel surrounding rock area is divided into multiple test areas, and multiple test points are set in the test areas. The maximum horizontal stress, minimum horizontal stress and vertical force of the test points are determined based on the stress relief method; Step 2: Calculate the maximum horizontal stress, minimum horizontal stress and vertical stress of the area to be tested based on the maximum horizontal stress, minimum horizontal stress and vertical stress of each point to be tested in the area to be tested; Step 3: According to the calculation, the maximum horizontal stress, the minimum horizontal stress and the vertical stress of the area to be measured are obtained, and the original rock stress of the area to be measured is comprehensively generated; Step 4: Bury a high-precision, long-lasting three-way mining stress monitoring sensor at the center of the area to be tested to collect its mining stress, and decompose the mining stress to generate axial stress, annular stress and radial stress; Step 5: Based on the axial stress, annular stress and radial stress, a stress threshold of the area to be tested is comprehensively generated, and the stress threshold of the same area to be tested is compared with the original rock stress. If the original rock stress exceeds the stress threshold, it indicates that the area to be tested is an abnormal area.

2. A method for locating abnormal stress areas of surrounding rocks in a mining disturbance area according to claim 1, characterized in that: The maximum horizontal stress, the minimum horizontal stress and the vertical stress of the measured point are obtained according to the following method: Using the hydraulic fracturing method, water is injected into the borehole to increase the pressure in the hole until the rock breaks. The fracture pressure is related to the ground stress. The initial fracture pressure, pore pressure and closing pressure of the rock at the test point are recorded to calculate the maximum horizontal stress, minimum horizontal stress and vertical stress. The formula is: s v =γ*h*g Among them, σ H , σ h , σ v Respectively represent the maximum horizontal stress, minimum horizontal stress, and vertical stress of the test point, P hf Indicates the initial fracture pressure of rock, P c Indicates closing pressure, P k represents the pore pressure, γ represents the density of the surrounding rock, h represents the drilling depth, and g represents the deadweight of the rock.

3. A method for locating abnormal stress areas of surrounding rocks in a mining disturbance zone according to claim 2, characterized in that: The method for obtaining the maximum horizontal stress, the minimum horizontal stress and the vertical stress in the measured area is: Determine the coordinates of the center point of the area to be measured and mark it as (x j ,y j , z j ), calibrate the coordinates of the i-th point to be measured in the area to be measured as (x i ,y i , z i ), i represents the index of the test point in the test area, and i∈[1,n], n is the number of test points in the test area, based on the stress data of each test point in the test area and the distance between the test point and the center point, the maximum horizontal stress, minimum horizontal stress and vertical stress of the center point are calculated, and used as the maximum horizontal stress, minimum horizontal stress and vertical stress of the test area, and the maximum horizontal stress, minimum horizontal stress and vertical stress of the center point are calculated, and the formula based on which is: in, Respectively represent the center point (x j ,y j , z j ), the maximum horizontal stress, minimum horizontal stress and vertical stress at Respectively represent the points to be tested (x i ,y i , z i ), the maximum horizontal stress, minimum horizontal stress and vertical stress at Represents the distance from the center point to the i-th point to be measured.

4. A method for locating abnormal stress areas of surrounding rocks in a mining disturbance zone according to claim 3, characterized in that: The formula for comprehensively generating the original rock stress is: Among them, σ 原 represents the original rock stress in the measured area, ω H ,ω h ,ω v are the weight ratios corresponding to the maximum horizontal stress, minimum horizontal stress and vertical stress, ω H >ω h >ω v , and satisfies ω H +ω h +ω v =1.

5. A method for locating abnormal stress areas of surrounding rocks in a mining disturbance zone according to claim 1, characterized in that: The axial stress, the hoop stress and the radial stress are obtained according to the method: Polar coordinate method is used to decompose the dynamic stress of the center point into the dynamic stress perpendicular to the yz plane in the x direction, the dynamic stress perpendicular to the xz plane in the y axis direction, and the dynamic stress perpendicular to the xy plane in the z axis direction. The axial stress is along the axis direction and is directly expressed by the measured dynamic stress perpendicular to the xy plane. The calculation of the hoop stress and radial stress requires determining the position angle and integrating the shear stress of the xy plane. The formula is: F z =s zz F θ =s xx *cos 2 (i)+s yy *sin 2 (θ)+2τ xy *cos(θ)*sin(θ) F r =s xx *sin 2 (i)+s yy *cos 2 (i)-2t xy *cos(θ)*sin(θ) Among them, F z 、F θ 、F r Respectively represent axial stress, hoop stress, radial stress, σ xx , σ yy , σ zz They represent the dynamic stress perpendicular to the yz plane in the x direction, the dynamic stress perpendicular to the xz plane in the y direction, and the dynamic stress perpendicular to the xy plane in the z direction, respectively. xy represents the shear stress in the xy plane, and θ represents the position angle in polar coordinates.

6. A method for locating abnormal stress areas of surrounding rocks in a mining disturbance area according to claim 5, characterized in that: The stress thresholds are generated based on the following method: T=a*F z +b*F θ +c*F r Wherein, T represents the stress threshold, a, b, and c represent the weight ratios of the corresponding axial stress, hoop stress, and radial stress, respectively, and a<b=c.

7. The method for locating abnormal stress areas of surrounding rocks in a mining disturbance zone according to claim 1, characterized in that: Compared with the generated original rock stress, if the original rock stress exceeds the stress threshold, it means that the area to be tested is a stress abnormal area, and the logic is as follows: Wherein, Q represents a logical value for judging whether the area to be tested is an abnormal area. When Q=0, it indicates that the original rock stress does not exceed the stress threshold, and part of the area to be tested is not an abnormal area; when Q=1, it indicates that the original rock stress exceeds the stress threshold, and part of the area to be tested is an abnormal area.

8. A system for locating abnormal area of ​​surrounding rock stress in a mining disturbance area, the abnormal area locating system is used to execute the method for locating abnormal area of ​​surrounding rock stress in a mining disturbance area according to any one of claims 1 to 7, comprising: A stress data acquisition module, wherein the stress data acquisition module divides the tunnel surrounding rock area into a plurality of test areas based on three-dimensional grid division, sets a plurality of test points in the test areas, and determines the maximum horizontal stress, the minimum horizontal stress, and the vertical stress of the test points based on a stress relief method; A gridding and interpolation calculation module, which calculates the maximum horizontal stress, the minimum horizontal stress and the vertical stress of the area to be measured based on the maximum horizontal stress, the minimum horizontal stress and the vertical stress of each point to be measured in the area to be measured; An original rock stress generation module, the original rock stress generation module is used to obtain the maximum horizontal stress, the minimum horizontal stress and the vertical stress of the area to be measured according to the calculation, and comprehensively generate the original rock stress of the area to be measured; A mining stress acquisition module, wherein the mining stress acquisition module is used to bury a high-precision, long-term three-way mining stress monitoring sensor at the center point of the area to be measured to collect the mining stress, and decompose the mining stress to generate axial stress, annular stress and radial stress; The stress threshold and anomaly detection module is used to comprehensively generate a stress threshold of the area to be tested based on the axial stress, annular stress and radial stress, and compare the stress threshold of the same area to be tested with the original rock stress. If the original rock stress exceeds the stress threshold, it means that the area to be tested is a stress abnormality area.

Images