A method and system for positioning stress abnormal areas of surrounding rock in a mining disturbance area

Abstract

The present application provides a kind of mine exploitation disturbed zone surrounding rock stress abnormal region positioning method and system, it is related to stress anomaly identification and positioning technical field, the present application obtains the data of original rock stress field, calculates the ground stress value of original rock stress field test point, based on the known test point data and the stress value calculated, roadway surrounding rock area is zoned by grid, and the maximum horizontal stress, minimum horizontal stress and vertical stress of the region to be measured are calculated using interpolation method, the original rock stress is generated in a comprehensive manner, and three-way mining stress data of surrounding rock are collected, including axial stress, hoop stress and radial stress, stress threshold is generated, and compared with original rock stress, the position and distribution of abnormal region are quickly and accurately judged.The present application realizes more effective mine disaster risk early warning and more reliable mine safety management through accurate stress calculation, flexible sensor layout and comprehensive analysis method.

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 region of surrounding rock in a mining disturbance zone. Background Art

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

[0003] Prior art publication CN108802825B discloses a method and system for locating coal and rock dynamic disasters using infrasound monitoring. The system is connected to a monitoring master station via a fiber optic network. The monitoring substation consists of three infrasound sensors arranged in a triangle, connected to an infrasound monitor. The master station is equipped with a computing server. The computing server at the master station filters the collected data, performs time-frequency analysis on the data using a short-time Fourier transform (STFT), compares the frequency energy density, extracts the main frequency band signal, and calculates the location using statistical analysis and time-delay estimation theory. The system monitors the infrasound signals generated by the fracture of loaded coal and rock. This method allows for non-contact, regional testing of the stress distribution of the surrounding rock in mines or tunnels, identifying areas of coal and rock fracture and stress anomalies.

[0004] Insufficient existing technology:

[0005] In existing technologies, although infrasound monitoring systems are capable of conducting non-contact regional tests, the data must undergo multiple processing steps, including filtering, short-time Fourier transform, and statistical analysis, which results in 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 system's performance under different working conditions.

[0007] Therefore, it is necessary to provide a method and system for locating abnormal stress areas in 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 form the prior art that is already known to a person 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 area, 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 stress area of ​​surrounding rock in a mining disturbance area, comprising the following steps:

[0012] Step 1: Based on 3D 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 stress 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 test area based on the maximum horizontal stress, minimum horizontal stress, and vertical stress of each test point in the test area;

[0014] Step 3: Based on the calculated maximum horizontal stress, minimum horizontal stress, and vertical stress of the test area, the in-situ rock stress of the test area is comprehensively generated;

[0015] Step 4: Bury a high-precision, long-lasting three-dimensional mining stress monitoring sensor at the center of the area to be measured to collect mining stress. Decompose the mining stress to generate axial stress, hoop stress, and radial stress.

[0016] Step 5: Based on the axial stress, hoop stress, and radial stress, a stress threshold of the area to be tested is comprehensively generated. 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, minimum horizontal stress and vertical stress of the test point are obtained according to the following method:

[0018] Using hydraulic fracturing, 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 closure 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 rock weight.

[0022] Furthermore, the maximum horizontal stress, minimum horizontal stress and 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 the maximum horizontal stress, minimum horizontal stress and vertical stress of the center point are used as the maximum horizontal stress, minimum horizontal stress and vertical stress of the test area. The maximum horizontal stress, minimum horizontal stress and vertical stress of the center point are calculated based on the following formula:

[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 measured (x i ,y i , z i ), the maximum horizontal stress, minimum horizontal stress and vertical stress at Indicates the distance from the center point to the i-th point to be measured.

[0026] Furthermore, the in-situ rock stress is comprehensively generated based on the following formula:

[0027]

[0028] Among them, σ 原represents the in-situ 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 satisfy ω H +ω h +ω v =1.

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

[0030] The polar coordinate method is used to decompose the dynamic stress at 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 axial direction and is directly expressed using 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 in 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 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] Among them, 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, when compared with the original rock stress, if the original rock stress exceeds the stress threshold, it means that the tested area is an abnormal area. The logic is as follows:

[0039]

[0040] Wherein, Q represents the 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 the 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 the 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 area. The abnormal area locating system is used to execute any of the above-mentioned methods for locating an abnormal area of ​​surrounding rock stress in a mining disturbance area, comprising:

[0042] A stress data acquisition module, which divides the tunnel surrounding rock area into multiple test areas based on three-dimensional grid division, sets multiple test points in the test areas, and determines the maximum horizontal stress, minimum horizontal stress, and 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, minimum horizontal stress, and vertical stress of the area to be measured based on the maximum horizontal stress, minimum horizontal stress, and vertical stress of each point to be measured in the area to be measured;

[0044] An in-situ rock stress generation module is used to comprehensively generate the in-situ rock stress of the area to be measured based on the calculated maximum horizontal stress, minimum horizontal stress, and vertical stress of the area to be measured;

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

[0046] The stress threshold and anomaly detection module is used to comprehensively generate a stress threshold of the test area based on the axial stress, hoop stress, and radial stress, and compare the stress threshold of the same test area with the original rock stress. If the original rock stress exceeds the stress threshold, it indicates that the test area is a stress anomaly 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 measured area by introducing grid partitioning and interpolation methods, thereby enhancing the real-time and accuracy of data processing. Through flexible sensor layout and the embedding of three-dimensional 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, capable of more accurately monitoring and identifying 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 multiple 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 in-situ rock stress exceeds the safe range, promptly identify and warn of abnormal areas, and improve the safety and reliability of the overall system.

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

[0051] Figure 1 Schematic diagram of the overall method 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 with reference to specific embodiments.

[0054] It should be noted that, unless otherwise defined, the technical or scientific terms used in the present invention should have the usual meanings understood by people with ordinary skills in the field to which the present invention belongs. The "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 position relationships. When the absolute position of the object being described changes, the relative position 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, comprising the following steps:

[0057] Step 1: Based on 3D 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 stress 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 test area based on the maximum horizontal stress, minimum horizontal stress, and vertical stress of each test point in the test area;

[0059] Step 3: Based on the calculated maximum horizontal stress, minimum horizontal stress, and vertical stress of the test area, the in-situ rock stress of the test area is comprehensively generated;

[0060] Step 4: Bury a high-precision, long-lasting three-dimensional mining stress monitoring sensor at the center of the area to be measured to collect mining stress. Decompose the mining stress to generate axial stress, hoop stress, and radial stress.

[0061] Step 5: Based on the axial stress, hoop stress, and radial stress, a stress threshold of the area to be tested is comprehensively generated. 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 by directly measuring the fracture pressure and closure pressure, hydraulic fracturing provides direct data on ground stress. Using formulas to calculate the maximum and minimum horizontal stresses and vertical stresses can convert the measured data into usable information, and factors such as initial fracture pressure, pore pressure, and closure pressure are included in 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 hydraulic fracturing, 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 closure 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 at which 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 cells, making the calculation and analysis of stress distribution more detailed and accurate. Gridding stress calculations for the entire region can predict the changing trends of ground stress at different locations, helping to identify potential hazardous areas for prevention and treatment. Interpolation fills gaps between known measurement points, enhancing data integrity and visualization. A grid cell size of 5 mm in length, width, and height is chosen to ensure both computational accuracy and operability and efficiency. Grid cells that are too large or too small will affect both accuracy and efficiency. Using a three-dimensional grid generation algorithm can automatically divide complex three-dimensional space into regular grid cells, facilitating systematic and consistent data processing. The weighting factors and distance calculation methods used in the formula ensure the rationality and accuracy of the interpolation results. Effective weighting and distance calculations can more realistically reflect the distribution characteristics of stress in space.

[0069] Therefore, the method for obtaining the maximum horizontal stress, minimum horizontal stress and 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 the maximum horizontal stress, minimum horizontal stress and vertical stress of the center point are used as the maximum horizontal stress, minimum horizontal stress and vertical stress of the test area. The maximum horizontal stress, minimum horizontal stress and vertical stress of the center point are calculated based on the following formula:

[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 measured (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 measured point; 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 stress value of each measuring point weighted by its distance to the center point. The closer the measuring point 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. The synthetic stress index can simplify the analysis process and facilitate 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 in-situ rock stress, based on the following formula:

[0075]

[0076] Among them, σ 原 represents the in-situ 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 satisfy ω 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 damage 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 damage 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 that of the maximum and minimum horizontal stresses on the rock mass.

[0077] It should be noted that by calculating the axial, circumferential 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, 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 according to the following method:

[0079] The polar coordinate method is used to decompose the dynamic stress at 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 axial direction and is directly expressed using 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 in 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 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 that can be used 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, 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 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, which is conducive to concentrating resources and energy on areas that need attention, improving the efficiency of detection and maintenance work, and ensuring that the structure operates within a safe range, avoiding risks caused by exceeding the design bearing capacity.

[0089] Therefore, it is necessary to compare the stress threshold with the in-situ stress. If the in-situ stress exceeds the stress threshold, it indicates that the待测区域 (to-be-detected area) is an abnormal area. The basis logic is as follows:

[0090]

[0091] Where, Q represents the logical value for judging whether the to-be-detected area is an abnormal area. When Q = 0, it indicates that the in-situ stress does not exceed the stress threshold, and the to-be-detected area is not an abnormal area; when Q = 1, it indicates that the in-situ stress exceeds the stress threshold, and the to-be-detected area 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, which divides the tunnel surrounding rock area into multiple test areas based on three-dimensional grid division, sets multiple test points in the test areas, and determines the maximum horizontal stress, minimum horizontal stress, and 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, minimum horizontal stress, and vertical stress of the area to be measured based on the maximum horizontal stress, minimum horizontal stress, and vertical stress of each point to be measured in the area to be measured;

[0095] An in-situ rock stress generation module is used to comprehensively generate the in-situ rock stress of the area to be measured based on the calculated maximum horizontal stress, minimum horizontal stress, and vertical stress of the area to be measured;

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

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

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

[0099] The above embodiments can be implemented in whole or in part by software, hardware, firmware, or any other combination thereof. When implemented using software, the above embodiments can be implemented in whole or in part in the form of a computer program product. Those skilled in the art will appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented by electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are performed by hardware or software 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 separate, and the components shown as units may or may not be physical units, and may be located in one place or distributed across multiple network units. Some or all of these units may be selected to achieve the purpose of this embodiment as needed.

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

Claims

1. A method for locating abnormal stress areas in surrounding rocks in a mining disturbance area, characterized in that: The specific steps include: Step 1: Based on 3D 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 stress 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 test area based on the maximum horizontal stress, minimum horizontal stress, and vertical stress of each test point in the test area; Step 3: Based on the calculated maximum horizontal stress, minimum horizontal stress, and vertical stress of the test area, the in-situ rock stress of the test area is comprehensively generated; Step 4: Bury a high-precision, long-lasting three-dimensional mining stress monitoring sensor at the center of the area to be measured to collect mining stress. Decompose the mining stress to generate axial stress, hoop stress, and radial stress. Step 5: Based on the axial stress, hoop stress, and radial stress, a stress threshold of the area to be tested is comprehensively generated. 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. The method for locating abnormal stress areas of surrounding rocks in a mining disturbance zone according to claim 1, characterized in that: The maximum horizontal stress, minimum horizontal stress and vertical stress of the measured point are obtained according to the following method: Using hydraulic fracturing, 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 closure 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 rock weight.

3. The method for locating abnormal stress areas of surrounding rocks in a mining disturbance zone according to claim 2, characterized in that: The maximum horizontal stress, minimum horizontal stress and vertical stress at the measured area are obtained as follows: 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 the maximum horizontal stress, minimum horizontal stress and vertical stress of the center point are used as the maximum horizontal stress, minimum horizontal stress and vertical stress of the test area. The maximum horizontal stress, minimum horizontal stress and vertical stress of the center point are calculated based on the following formula: 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 measured (x i ,y i , z i ), the maximum horizontal stress, minimum horizontal stress and vertical stress at Indicates the distance from the center point to the i-th point to be measured.

4. A method for locating abnormal stress areas in 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 in-situ 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 satisfy ω H +ω h +ω v =1.

5. The method for locating abnormal stress areas of surrounding rocks in a mining disturbance zone according to claim 1, characterized in that: The axial stress, hoop stress and radial stress are obtained according to the following method: The polar coordinate method is used to decompose the dynamic stress at 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 axial direction and is directly expressed using 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 in 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 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 in surrounding rocks in a mining disturbance zone 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 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。 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 original rock stress, if the original rock stress exceeds the stress threshold, it means that the tested area is a stress abnormal area. The logic is as follows: Wherein, Q represents the 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 the 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 the part of the area to be tested is an abnormal area.

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

Images