Tunnel rock burst pretreatment method and system based on while-drilling multi-physics field parameters

Abstract

The invention mainly relates to the technical field of underground engineering safety, and provides a tunnel rockburst pretreatment method and system based on multi-physics field parameters while drilling in order to improve the rockburst risk recognition and pressure relief effect through the multi-physics field parameters in the deep hard rock tunnel drilling process. The core of the method is that a surrounding rock three-dimensional geologic model is established based on multi-source physical field parameters in the surrounding rock while-drilling process; establishing a rockburst risk index prediction model of the surrounding rock three-dimensional geologic model based on an improved D-S evidence theory to calculate a rockburst risk index; and determining simulated pressure relief parameters according to the rockburst risk indexes, performing simulated pressure relief, judging whether the change rate of the multi-source physical field parameters meets the set requirements or not, if so, performing pressure relief according to the simulated pressure relief parameters, otherwise, re-determining the pressure relief parameters, and performing pressure relief until the change rate of the multi-source physical field parameters meets the set requirements.

Description

Technical Field

[0001] The present invention mainly relates to the technical field of underground engineering safety, and in particular to a tunnel rockburst pretreatment method and system based on multi-physical field parameters while drilling. Background Art

[0002] During the excavation of deep hard rock tunnels, rock burst disasters have become a key problem restricting project safety due to their suddenness and destructive power.

[0003] Current rockburst pretreatment methods, such as local blasting unloading and shotcrete support, have the following defects: ① The advanced prediction method is single, and only single stress or acoustic wave monitoring is performed through borehole stress gauges, TSP seismic wave methods, etc., which cannot fully capture the multi-physical field coupling information such as surrounding rock stress, acoustic wave velocity, and resistivity. As a result, it is impossible to accurately identify the critical conditions of rockburst under complex stress states, the prediction accuracy is insufficient, and it is difficult to accurately judge the rockburst risk in advance; ② Unloading measures are often implemented passively after the signs of rockburst appear. The layout parameters of the unloading holes are mostly determined by engineering experience, with fixed spacing and uniform depth. They are not linked to the advanced prediction system and cannot be dynamically adjusted according to real-time geological changes to form a closed-loop control of "forecast-adjustment-implementation". There is a problem of lagging and passive unloading process. Summary of the Invention

[0004] The technical problem to be solved by the present invention is to provide a tunnel rockburst pretreatment method and system based on multi-physical field parameters while drilling, with the aim of jointly predicting the tunnel rockburst risk through multi-physical field parameters during the drilling process of deep-buried hard rock tunnels, thereby improving the prediction effect of rockburst risk and the pressure relief effect.

[0005] The technical solution adopted by the present invention to solve the above technical problems is: In one aspect, the present invention provides a tunnel rockburst pretreatment method based on multi-physics field parameters while drilling, the method comprising: Establish a three-dimensional geological model of the surrounding rock based on multi-source physical field parameters during the surrounding rock drilling process; A rockburst risk index prediction model for the surrounding rock 3D geological model is established based on the improved DS evidence theory; During the drilling process, the multi-source physical field parameters are collected according to the set drilling distance, and the rock burst risk index is calculated based on the rock burst risk index prediction model; The simulated pressure relief parameters are determined according to the rock burst risk index, and simulated pressure relief is performed to determine whether the change rate of the multi-source physical field parameters meets the set requirements. If so, pressure relief is performed according to the simulated pressure relief parameters. Otherwise, the pressure relief parameters are re-determined until the change rate of the multi-source physical field parameters meets the set requirements, and then pressure relief is performed according to the corresponding pressure relief parameters.

[0006] Furthermore, the three-dimensional geological model of the surrounding rock is: ,in 、 and are the acoustic wave velocity, electromagnetic response intensity and stress value during drilling, 、 、 are the weight coefficients of acoustic wave velocity, electromagnetic response intensity and stress value, , weight coefficient 、 、 The method to determine is: Calculate rock burst tendency based on ground stress or rock damage criterion ; According to the initial weight coefficient set 、 、 and the acoustic wave velocity during drilling , electromagnetic response intensity and stress values Calculate and predict rock burst tendency ; Based on the calculation of rock burst tendency and predicting rock burst tendency Initial weight coefficient 、 、 Make corrections and determine the corrected weight coefficient 、 、 .

[0007] Furthermore, the initial weight coefficient is 、 、 Make corrections to obtain the corrected weight coefficient 、 、 value.

[0008] Furthermore, the rockburst risk index prediction model is ,in is the rockburst risk index, The maximum principal stress detected during tunnel boring, The delayed uniaxial compressive strength is tested during tunnel boring. The apparent resistivity detected during tunnel boring process, is the baseline longitudinal wave velocity of the intact rock mass, is the base apparent resistivity.

[0009] Furthermore, the simulation unloading parameters are determined based on the rock burst risk index, including: Determine the rock burst risk level based on the rock burst risk index; The scope of the pressure relief area, spacing, diameter, depth and pressure relief pressure of the pressure relief holes are determined based on the rock burst risk level.

[0010] Furthermore, judging whether the rate of change of the multi-source physical field parameters meets the set requirements includes: obtaining the attenuation rate of the ultrasonic wave velocity and the decrease in apparent resistivity of the tunnel after simulated pressure unloading, judging whether the decrease in ultrasonic wave velocity and apparent resistivity reaches the set value, if it reaches the set value, then performing pressure unloading according to the simulated pressure unloading parameters, otherwise, calculating whether the rock burst risk index after simulated pressure unloading is less than the set rock burst risk index, if so, directly increasing the pressure unloading pressure, otherwise it is necessary to recalculate the rock burst risk index according to the multi-source physical field parameters collected during the surrounding rock drilling process, and reset the pressure unloading hole spacing, pressure unloading hole spacing, hole diameter, depth and pressure unloading pressure.

[0011] Furthermore, the spacing between the pressure relief holes is calculated based on the critical crack extension length calculated by Griffith fracture theory.

[0012] Furthermore, the method also includes: when the rock burst level is high risk, setting a rupture pressure safety factor, setting an actual relief pressure based on the rupture pressure safety factor and the rupture pressure, and when the rock burst level is medium risk, using the rupture pressure as the relief pressure.

[0013] Furthermore, the method also includes: decompression based on hydraulic fracturing, gas pressure fracturing or laser induced fracture.

[0014] On the other hand, the present invention also provides a tunnel rockburst preprocessing system based on multi-physical field parameters while drilling, the system comprising a data acquisition module, a data processing module, a guidance control system and an execution module; The data acquisition module is used to collect multi-source physical field parameters during the surrounding rock drilling process; The data processing module is used to calculate the rockburst risk index based on the multi-source physical field parameters collected during the surrounding rock drilling process; The guidance control system is used to determine the rock burst risk level according to the rock burst risk index, generate simulated pressure relief parameters, verify whether the simulated pressure relief parameters meet the set requirements, and generate a pressure relief instruction based on the pressure relief parameters that meet the set requirements and transmit it to the execution module; The execution module is used to perform pressure relief according to the pressure relief instruction to eliminate the risk of rock burst.

[0015] Beneficial effects of the present invention: (1) The tunnel rockburst pretreatment method based on multi-physical field parameters while drilling described in the present invention establishes a three-dimensional geological model of the surrounding rock through multiple physical field parameters during the tunnel drilling process, and establishes a rockburst risk index prediction model based on the three-dimensional geological model of the surrounding rock to predict the rockburst risk index during the drilling process. On the one hand, it can achieve quantitative evaluation and accurate prediction of the rockburst risk of the surrounding rock, and avoid the potential fracture development risk under normal circumstances being missed by monitoring with a single prediction parameter; on the other hand, based on the critical extension length of the fracture and the annular pressure relief zone theory, the pressure relief process is simulated and the pressure relief parameters are selected to accurately match the pressure relief parameters with the stress state of the surrounding rock and the law of fracture extension, which can improve the accuracy and safety of the engineering execution process; (2) During the simulated pressure relief process, the hydraulic fracturing pressure formula is used to accurately control the expansion of the cracks along the direction of the maximum principal stress, which helps the cracks form effective stress release channels during the actual pressure relief process. Through real-time monitoring of ultrasonic and transient electromagnetic, it is ensured that the maximum principal stress reduction in the target area reaches the target and the probability of rock burst is reduced. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] Figure 1 This is a side view of the arrangement of annular pressure relief holes in front of the tunnel face; Figure 2 This is the front view of the arrangement of the annular pressure relief holes in front of the tunnel face; Figure 3 This is a schematic diagram of the hydraulic fracturing sealing device and crack expansion principle; Figure 4 This is a schematic diagram of multi-source physical field parameter acquisition; Figure 5 This is a schematic diagram for calculating the rockburst risk index; Figure 6 It is a schematic diagram for setting pressure relief parameters; Figure 7 Schematic diagram of rockburst pretreatment.

[0017] 1-Tunnel face contour line, 2-Annular pressure relief area contour line, 3-Tunnel face pressure relief hole, 4-Upper pressure relief hole in the annular pressure relief area, 5-Hole sealer, 6-Hydraulic splitting fracture, 7-Surrounding rock. DETAILED DESCRIPTION

[0018] The core of the tunnel rockburst pretreatment method based on multi-physics field parameters while drilling described in the present invention to solve the above technical problems is: Based on the establishment of a three-dimensional geological model of the surrounding rock, the types of physical field parameters for measuring the rockburst risk during the surrounding rock drilling process and the weights of each parameter are determined. Based on the types of physical field parameters and the weights of each parameter, a rockburst risk index prediction model is established to predict the rockburst risk index. Based on the rockburst risk index, the rockburst risk level is determined. Different pressure relief parameters are matched for different rockburst risk levels. Pressure relief is simulated based on the pressure relief parameters to determine whether the simulated pressure relief results meet the set requirements. If so, the pressure relief parameters are executed according to the simulated pressure relief parameters.

[0019] like Figure 4-Figure 7 As shown, the tunnel rockburst pretreatment method based on multi-physics field parameters while drilling described in the present invention mainly includes four stages: multi-source physical field parameter acquisition, rockburst risk index calculation, pressure relief parameter setting and rockburst pretreatment.

[0020] Multi-source physical field parameter acquisition A circular downhole sensor array is set at the front end of the drill rod of the tunnel boring machine. The downhole sensor array is installed within 30 cm of the front end of the drill rod and advances synchronously with the drilling direction. The downhole sensor array collects physical field data every time it advances a set distance during the drilling process.

[0021] like Figure 4 As shown in the figure, the LWD sensor array specifically includes: borehole stress sensors for real-time monitoring of the three-dimensional stresses in the surrounding rock ahead of the tunnel face; ultrasonic transducers for transmitting and receiving electromagnetic signals to measure the longitudinal wave velocity of the surrounding rock; and transient electromagnetic probes for collecting apparent resistivity. The LWD sensor array can simultaneously collect physical field parameters such as stress, longitudinal wave velocity, and apparent resistivity of the surrounding rock 0-30 meters ahead of the tunnel face during tunnel drilling.

[0022] Calculation of rockburst risk index like Figure 5 As shown, first, a three-dimensional geological model of the surrounding rock is constructed based on the collected multi-source physical field parameters. The three-dimensional geological model is: ,in 、 and are the acoustic wave velocity, electromagnetic response intensity and stress value collected during the drilling process, 、 、 are the weight coefficients of acoustic wave velocity, electromagnetic response intensity and stress value in the three-dimensional geological model.

[0023] Weight coefficient 、 、 Determined through BP neural network training, it is used to reflect the contribution of different physical field parameters in different surrounding rock types to rock burst risk. 、 、 The acquisition specifically includes: Define the BP neural network structure and parameters: The BP neural network input layer is used to receive the acoustic wave velocity, electromagnetic response intensity and stress value. The physical field parameters of the input layer are vector Form input, number of neurons in the input layer ; Number of neurons in the hidden layer , the activation function is Sigmoid function ; Output layer outputs three-dimensional geological model indicators , the number of neurons , the activation function is a linear function , the output expression is the predicted value , To predict rock burst tendency.

[0024] Prepare training data: The collection group contains Geological sample data, building a sample set , Used to identify the group of sample data, is the total number of samples, To calculate the rock burst tendency of rock mass.

[0025] It can be obtained by two methods: ground stress and rock damage; Ground stress determination and acquisition Specifically include: when the maximum principal stress Exceeds the uniaxial compressive strength of rock mass 70% of the corresponding ; Maximum principal stress Uniaxial compressive strength of rock mass For every 10% decrease, Decrease linearly by 0.2, refer to the following formula: ; Rock damage determination and acquisition Specifically include: When the ultrasonic wave velocity attenuation rate , apparent resistivity drop hour, , each decrease The wave velocity attenuation rate, Decrease by 0.3.

[0026] If the in-situ stress criterion and the rock damage criterion If the difference is greater than 0.3, microseismic monitoring is initiated for additional verification, and the microfracture event density N (times / m³) is used for correction. ,in is the ground stress criterion value, is the rock damage criterion value.

[0027] Weight initialization and hyperparameter setting: 、 、 Assign random initial value 、 、 ~ ; Initial learning rate The value is 0.01 (adaptive adjustment); the maximum number of iterations ; Error threshold .

[0028] Forward propagation and loss calculation: For the sample set Group data and calculate predicted values , using mean square error as the loss function , is the true value.

[0029] Back propagation and weight update: Weighting of sound wave velocity , partial derivatives , weight update formula .

[0030] Weight of electromagnetic intensity , partial derivatives , weight update formula .

[0031] Corresponding stress value weight , partial derivatives , update formula .

[0032] Iteration termination and result output: Repeat the forward propagation, loss calculation, and back propagation process until the number of iterations reaches Or loss function , and get the optimized weight coefficient 、 、 .

[0033] In actual engineering, considering issues such as time progress, the weight coefficient 、 、 The initial value can be directly selected based on engineering experience and dynamically adjusted according to the rock type, such as the typical value of granite 、 、 .

[0034] Based on the determined acoustic wave velocity, electromagnetic response intensity, and three-dimensional stress value weight coefficient, an improved DS evidence theory is used to perform spatiotemporal coupling analysis and construct a three-dimensional rockburst risk index prediction model: ,in is the rockburst risk index, is the maximum principal stress, is the delayed uniaxial compressive strength, is the baseline longitudinal wave velocity of the intact rock mass, is the base apparent resistivity.

[0035] During the drilling process, multi-source physical field parameters are collected according to the set drilling distance, and the rock burst risk index is calculated based on the rock burst risk index prediction model; the simulated pressure relief parameters are determined according to the rock burst risk index, and simulated pressure relief is performed to determine whether the change rate of the multi-source physical field parameters meets the set requirements. If so, pressure relief is performed according to the simulated pressure relief parameters; otherwise, the pressure relief parameters are re-determined until the change rate of the multi-source physical field parameters meets the set requirements.

[0036] Pressure relief parameter setting When the rockburst risk index When the rockburst risk index is , it indicates that the surrounding rock in front of the tunnel face is in a state of high stress concentration, and the rock mass integrity is high, and the rockburst risk level is high; when the rockburst risk index is , indicating that the stress state of the rock mass at the tunnel face is relatively mild, the unloading parameters are based on the principle of "moderate intervention + cost control", and the rock burst risk level is medium risk. , it means that the rock burst risk level is low and pressure relief pretreatment is not considered for the time being.

[0037] For high-risk areas, such as Figure 6 As shown in the figure, the selection of simulation decompression parameters includes: Annular pressure relief area contour line 2 range: Figure 1 and Figure 2 As shown, with the tunnel face contour line 1 as the center, the longitudinal range is 10~15m (covering the risk area + 5m safety margin), and the radial radius is , = is the tunnel diameter, ensuring that all potential stress concentration areas are covered. Spacing of pressure relief holes: Calculate the critical crack extension length based on Griffit fracture theory ( is the fracture toughness of rock (MPa・m¹ / ²), is the maximum principal stress (MPa) and the spacing between pressure relief holes is Depth and angle of pressure relief hole: The depth of the pressure relief hole needs to penetrate the stress concentration area, which is 15~20m (for example, when the burial depth is 1200m, the stress influence depth is about 12m, and the hole depth is designed to be 15m). The drilling angle is at a certain depression angle with the tunnel axis to ensure that the cracks expand along the high stress surface. Specifically, the drilling angle is designed to be at a depression angle of 15°-30° with the tunnel axis. The pressure relief pressure is determined according to the fracture pressure. And the set burst pressure safety factor a (a>1) is set to MPa.

[0038] For medium-risk areas, the parameter settings for simulated pressure relief include: the pressure relief range is reduced to 8~12m in the longitudinal direction and 8~12m in the radial direction relative to the high-risk area. , covering only the core risk area. Spacing and depth of pressure relief holes: Calculation of critical crack extension length , take the hole spacing (Compared The drilling volume is reduced by 50%). The hole depth is 10m (thickness of the penetrating stress concentration area + 2m safety margin). Calculation of pressure relief pressure: If the calculated ,Since the rock pressure state in the medium risk area is relatively mild, there is no need to set a ,fracture pressure safety factor, and the fracture pressure of the surrounding rock can be ,directly used as the pressure relief pressure.

[0039] The pressure relief method can be selected according to the actual engineering situation, such as hydraulic fracturing, gas pressure fracturing or laser induced fracturing to create cracks.

[0040] In this embodiment, hydraulic fracturing technology is used. The pressure relief range, pressure relief hole spacing and other pressure relief parameters are set according to the rock burst risk. After the water injection equipment is connected and before the water injection begins, Figure 3 As shown, a sealing device 5 is used to seal one end of the pressure relief hole 3 at the tunnel face to ensure that no hydraulic loss occurs during water injection, and then water injection is started. Figure 3 As shown in Figure 1, when the water pressure reaches the fracture pressure, the hydraulic fracture 6 starts to expand from the upper pressure relief hole 4 in the annular pressure relief zone along the direction of the maximum principal stress, thereby effectively reducing the concentrated stress of the surrounding rock 7 in the target area. The water injection time t is calculated using the fracture expansion model: ,in is the dynamic viscosity of the rock mass, The target fracture radius is 0.001mm. After the water injection fracturing is completed, the U-TEM ultrasonic-transient electromagnetic generator is started to induce the expansion of rock microcracks through acoustic-electric coupling, which can further disperse the stress concentration area. 100kHz ultrasonic waves and transient electromagnetic signals with an off time of less than or equal to 1μs are emitted to monitor the ultrasonic wave velocity and apparent resistivity before pressure relief in real time, and calculate the ultrasonic wave velocity attenuation rate before and after simulated pressure relief. and apparent resistivity drop ,when and When the cracks are fully developed and the stress release meets the standard, the water pressure will be automatically increased to Repeat the cracking until the index requirements are met.

[0041] Rockburst pretreatment Verify whether the simulated pressure relief parameters meet the set requirements, and generate pressure relief instructions based on the pressure relief parameters that meet the set requirements and transmit them to the execution equipment for rock burst pretreatment.

[0042] The present invention also provides a tunnel rockburst preprocessing system based on multi-physical field parameters while drilling, the system comprising a data acquisition module, a data processing module, a guidance control system and an execution module; The data acquisition module is used to collect multi-source physical field parameters during the surrounding rock drilling process; The data processing module is used to calculate the rockburst risk index based on the multi-source physical field parameters collected during the surrounding rock drilling process; The guidance control system is used to determine the rock burst risk level according to the rock burst risk index, generate simulated pressure relief parameters, verify whether the simulated pressure relief parameters meet the set requirements, and generate a pressure relief instruction based on the pressure relief parameters that meet the set requirements and transmit it to the execution module; The execution module is used to perform pressure relief according to the pressure relief instruction to eliminate the risk of rock burst; Figure 7 As shown, the execution module includes the drilling machine arm system and drilling equipment. The drilling machine arm system is driven by a six-axis servo motor and supports 360-degree rotation and 45-degree pitch adjustment. The drilling equipment is a 76mm diameter hydraulic rock drill. The depth control accuracy is controlled within plus or minus 2cm through real-time calibration of the laser rangefinder. After the execution module completes the single-cycle decompression, the tunnel boring machine continues to drill 2m, and the drilling sensor is used to verify whether the rock burst risk index has dropped to If the standard is not met, the spacing between the pressure relief holes should be increased to , ensuring that risks are controllable.

[0043] Flexible sprayed anchor support is promptly implemented behind the pressure relief area, spraying C25 fiber concrete (thickness 15cm), Φ22 mortar anchor rods, length 3m, spacing 1.2m×1.2m, forming a coordinated protection system with the advance pressure relief.

Claims

1. A tunnel rockburst pretreatment method based on multi-physics field parameters while drilling, characterized by: The method comprises: Establish a three-dimensional geological model of the surrounding rock based on multi-source physical field parameters during the surrounding rock drilling process; A rockburst risk index prediction model for the surrounding rock 3D geological model is established based on the improved DS evidence theory; During the drilling process, the multi-source physical field parameters are collected according to the set drilling distance, and the rock burst risk index is calculated based on the rock burst risk index prediction model; The simulated pressure relief parameters are determined according to the rock burst risk index, and simulated pressure relief is performed to determine whether the change rate of the multi-source physical field parameters meets the set requirements. If so, pressure relief is performed according to the simulated pressure relief parameters. Otherwise, the pressure relief parameters are re-determined until the change rate of the multi-source physical field parameters meets the set requirements, and then pressure relief is performed according to the corresponding pressure relief parameters.

2. The tunnel rockburst pretreatment method based on multi-physics field parameters while drilling according to claim 1, characterized in that: The three-dimensional geological model of the surrounding rock is: ,in 、 and are the acoustic wave velocity, electromagnetic response intensity and stress value during drilling, 、 、 are the weight coefficients of acoustic wave velocity, electromagnetic response intensity and stress value, , weight coefficient 、 、 The method to determine is: Calculate rock burst tendency based on ground stress or rock damage criterion ; According to the initial weight coefficient set 、 、 and the acoustic wave velocity during drilling , electromagnetic response intensity and stress values Calculate and predict rock burst tendency ; Based on the calculation of rock burst tendency and predicting rock burst tendency Initial weight coefficient 、 、 Make corrections and determine the corrected weight coefficient 、 、 .

3. The tunnel rockburst pretreatment method based on multi-physics field parameters while drilling according to claim 2, characterized in that: Based on the BP neural network model, the initial weight coefficient 、 、 Make corrections to obtain the corrected weight coefficient 、 、 value.

4. The tunnel rockburst pretreatment method based on multi-physics field parameters while drilling according to claim 2, characterized in that: The rockburst risk index prediction model is: ,in is the rockburst risk index, The maximum principal stress detected during tunnel boring, The delayed uniaxial compressive strength is tested during tunnel boring. The apparent resistivity detected during tunnel boring process, is the baseline longitudinal wave velocity of the intact rock mass, is the base apparent resistivity.

5. The tunnel rockburst pretreatment method based on multi-physics field parameters while drilling according to claim 3, characterized in that: The simulation unloading parameters determined according to the rockburst risk index include: Determine the rock burst risk level based on the rock burst risk index; The scope of the pressure relief area, spacing, diameter, depth and pressure relief pressure of the pressure relief holes are determined based on the rock burst risk level.

6. The tunnel rockburst pretreatment method based on multi-physics field parameters while drilling according to claim 5, characterized in that: Determining whether the rate of change of the multi-source physical field parameters meets the set requirements includes: obtaining the attenuation rate of the ultrasonic wave velocity and the decrease in apparent resistivity of the tunnel after simulated pressure unloading, and determining whether the decrease in ultrasonic wave velocity and apparent resistivity reaches the set value. If so, pressure unloading is performed according to the simulated pressure unloading parameters. Otherwise, calculating whether the rock burst risk index after simulated pressure unloading is less than the set rock burst risk index. If so, directly increase the pressure unloading pressure. Otherwise, it is necessary to recalculate the rock burst risk index based on the multi-source physical field parameters collected during the surrounding rock drilling process, and reset the pressure unloading hole spacing, pressure unloading hole spacing, hole diameter, depth and pressure unloading pressure.

7. The tunnel rockburst pretreatment method based on multi-physics field parameters while drilling according to claim 5, characterized in that: The spacing between the pressure relief holes is calculated based on the critical crack extension length calculated by Griffith fracture theory.

8. The tunnel rockburst pretreatment method based on multi-physics field parameters while drilling according to claim 5, characterized in that: The method further includes: when the rock burst level is high risk, setting a rupture pressure safety factor, setting an actual relief pressure based on the rupture pressure safety factor and the rupture pressure, and when the rock burst level is medium risk, using the rupture pressure as the relief pressure.

9. The tunnel rockburst pretreatment method based on multi-physics field parameters while drilling according to claim 5, characterized in that: The method further includes: decompression based on hydraulic fracturing, gas fracturing or laser induced fracture.

10. A tunnel rockburst pretreatment system based on multi-physics field parameters while drilling, for implementing the tunnel rockburst pretreatment method based on multi-physics field parameters while drilling according to any one of claims 1 to 9, characterized in that: The system includes a data acquisition module, a data processing module, a guidance control system and an execution module; The data acquisition module is used to collect multi-source physical field parameters during the surrounding rock drilling process; The data processing module is used to calculate the rockburst risk index based on the multi-source physical field parameters collected during the surrounding rock drilling process; The guidance control system is used to determine the rock burst risk level according to the rock burst risk index, generate simulated pressure relief parameters, verify whether the simulated pressure relief parameters meet the set requirements, and generate a pressure relief instruction based on the pressure relief parameters that meet the set requirements and transmit it to the execution module; The execution module is used to perform pressure relief according to the pressure relief instruction to eliminate the risk of rock burst.

Images