A method and system for evaluating time-delay rockburst induced by blasting disturbance

Abstract

The present invention provides an evaluation method and system for time-delay rockburst induced by blasting disturbance, which relates to the technical field of tunnel engineering, including selecting basic dimensions from multiple physical factors based on dimensional analysis, and constructing a general model of blasting impact stress by using the basic dimensions and basic models; establishing a numerical analysis model of the tunnel according to parameter information of the tunnel to be blasted, simulating vibration acceleration data and blasting impact stress data of surrounding rock particles in adjacent tunnels facing blasting under different charge amounts and different blasting center distances in the numerical analysis model; substituting the acceleration data and blasting impact stress data into the general model to fit a blasting impact stress calculation model of the tunnel to be blasted; constructing an evaluation model for the probability of inducing time-delay rockburst in adjacent tunnels when the tunnel to be blasted is blasted by using the blasting impact stress calculation model. The present invention provides a set of systematic evaluation systems to evaluate the risk of time-delay rockburst, and provide scientific guidance and early warning for construction.

Description

Technical Field

[0001] The present invention relates to the technical field of tunnel engineering, and in particular to an evaluation method and system for time-delay rockburst induced by blasting disturbance. Background Art

[0002] In the field of tunnel construction, especially in the construction of deep tunnels, drilling and blasting is a common technical means of excavation. However, the strong shock wave generated by this method when blasting the subsequent tunnel often causes serious damage to the surrounding rock of the preceding tunnel, which may trigger a disaster called time-delay rockburst. This type of rockburst has a time-delay characteristic, that is, it will not occur until a period of time after the excavation area, which makes prediction and prevention extremely difficult. Time-delay rockburst not only poses a direct threat to the safety of construction workers, but may also lead to construction delays and increased costs. Despite its increasing importance, the current prediction method for this complex phenomenon is still immature, and there is a lack of a systematic evaluation system to assess the risk of rockburst, which makes it impossible to provide scientific guidance and early warning for construction. Summary of the invention

[0003] The purpose of the present invention is to provide a method and system for evaluating time-delay rockburst induced by blasting disturbance, so as to improve the above-mentioned problem. In order to achieve the above-mentioned purpose, the technical solution adopted by the present invention is as follows:

[0004] In a first aspect, the present application provides a method for evaluating a time-delay rockburst induced by a blasting disturbance, comprising:

[0005] Acquire multiple physical factors related to blasting impact stress, and construct a basic model of blasting impact stress using the multiple physical factors;

[0006] Selecting basic dimensions from multiple physical factors based on dimensional analysis, and constructing a general model of blasting impact stress using the basic dimensions and basic models;

[0007] Acquiring parameter information of a tunnel to be blasted, and establishing a numerical analysis model of the tunnel according to the parameter information of the tunnel to be blasted, wherein the tunnel to be blasted is a double-hole single-line tunnel;

[0008] The vibration acceleration data and blasting impact stress data of surrounding rock mass points in adjacent tunnels were simulated in the numerical analysis model under different charge amounts and different blasting center distances;

[0009] Substituting the acceleration data and the blasting impact stress data into the general model to obtain a blasting impact stress calculation model of the tunnel to be blasted;

[0010] An evaluation model for the probability of inducing a time-delayed rockburst in an adjacent tunnel when the tunnel to be blasted is blasted is constructed using the blasting impact stress calculation model.

[0011] In a second aspect, the present application also provides an evaluation system for blasting disturbance-induced time-delay rockburst, comprising:

[0012] Acquisition module: acquiring multiple physical factors related to blasting impact stress, and constructing a basic model of blasting impact stress using the multiple physical factors;

[0013] Selection module: selecting basic dimensions from multiple physical factors based on dimensional analysis, and constructing a general model of blasting impact stress using the basic dimensions and basic models;

[0014] Model building module: obtaining parameter information of the tunnel to be blasted, and building a numerical analysis model of the tunnel according to the parameter information of the tunnel to be blasted, wherein the tunnel to be blasted is a double-hole single-line tunnel;

[0015] Simulation module: simulates the vibration acceleration data and blasting impact stress data of the surrounding rock mass points in adjacent tunnels under different charge amounts and different blast center distances in the numerical analysis model;

[0016] Fitting module: using acceleration data and blasting impact stress data to substitute into the general model to fit the blasting impact stress calculation model of the tunnel to be blasted;

[0017] Calculation module: An evaluation model for the probability of inducing a time-delay rockburst in an adjacent tunnel when the tunnel to be blasted is blasted is constructed based on the blasting impact stress calculation model.

[0018] In a third aspect, the present application also provides an evaluation device for a time-delay rockburst induced by a blasting disturbance, comprising:

[0019] Memory for storing computer programs;

[0020] A processor is used to implement the steps of the method for evaluating time-delay rockburst induced by blasting disturbance when executing the computer program.

[0021] In a fourth aspect, the present application further provides a readable storage medium having a computer program stored thereon, and when the computer program is executed by a processor, the steps of the above-mentioned method for evaluating time-lag rockburst induced by blasting disturbance are implemented.

[0022] The beneficial effects of the present invention are:

[0023] The present invention constructs a dimensionless expression of blasting impact stress through dimensional analysis, which improves the universality and adaptability of the model. Combined with the parameters of the tunnel to be blasted, a numerical simulation is performed to obtain the vibration acceleration and blasting impact stress data under different charges and blasting center distances. Combined with these data, a calculation model for blasting impact stress is fitted, providing a quantitative tool for rock burst risk assessment. The present invention significantly improves the safety and reliability of tunnel engineering blasting operations, optimizes blasting design, reduces rock burst risks, ensures the safety of construction personnel, improves construction efficiency, and has important engineering application value.

[0024] Other features and advantages of the present invention will be described in the following description, and partly become apparent from the description, or be understood by implementing the embodiments of the present invention. The purpose and other advantages of the present invention can be realized and obtained by the structures particularly pointed out in the written description, claims, and drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0025] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required for use in the embodiments are briefly introduced below. It should be understood that the following drawings only show certain embodiments of the present invention and therefore should not be regarded as limiting the scope. For ordinary technicians in this field, other related drawings can be obtained based on these drawings without creative work.

[0026] Figure 1 It is a schematic flow chart of a method for evaluating a time-delay rockburst induced by blasting disturbance described in an embodiment of the present invention;

[0027] Figure 2 It is a schematic diagram of blast disturbance of adjacent tunnels facing blast in an embodiment of the present invention;

[0028] Figure 3 It is a schematic diagram of the structure of the evaluation system for the time-delay rockburst induced by blasting disturbance described in an embodiment of the present invention;

[0029] Figure 4 It is a schematic diagram of the structure of the evaluation equipment for the time-delay rockburst induced by blasting disturbance described in an embodiment of the present invention.

[0030] Markings in the figure:

[0031] 800. Evaluation equipment for time-delay rockburst induced by blasting disturbance; 801. Processor; 802. Memory; 803. Multimedia component; 804. I / O interface; 805. Communication component. DETAILED DESCRIPTION

[0032] In order to make the purpose, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, rather than all of the embodiments. The components of the embodiments of the present invention generally described and shown in the drawings here can be arranged and designed in various different configurations. Therefore, the following detailed description of the embodiments of the present invention provided in the drawings is not intended to limit the scope of the claimed invention, but merely represents selected embodiments of the present invention. Based on the embodiments in the present invention, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present invention.

[0033] It should be noted that similar reference numerals and letters represent similar items in the following drawings, so once an item is defined in one drawing, it does not need to be further defined and explained in the subsequent drawings. At the same time, in the description of the present invention, the terms "first", "second", etc. are only used to distinguish the description and cannot be understood as indicating or implying relative importance.

[0034] Embodiment 1:

[0035] This embodiment provides a method for evaluating time-delay rockburst induced by blasting disturbance.

[0036] See also Figure 1 , the figure shows that the method comprises:

[0037] S1. Acquire multiple physical factors related to blasting impact stress, and construct a basic model of blasting impact stress using multiple physical factors;

[0038] Specifically, detailed information of the tunnel project is collected, including tunnel design drawings, geological survey reports, etc., to determine that the propagation of blasting shock waves in the surrounding rock is related to physical factors such as particle vibration acceleration, charge amount, propagation path, and geotechnical medium conditions (such as geotechnical properties, joints, and geological structures). The dimensional parameters of the relevant physical factors are shown in Table 1.

[0039] Table 1 Dimensional parameters

[0040]

[0041] In Table 1, M represents the mass dimension, L represents the length dimension, and T represents the time dimension.

[0042] In this embodiment, the basic model of blasting impact stress is:

[0043] ;

[0044] Based on the above embodiments, the method further includes:

[0045] S2. Selecting a basic dimension from a plurality of physical factors based on a dimensional analysis method, and generating a general model of blasting impact stress using the basic dimension and the basic model;

[0046] Specifically, step S2 includes:

[0047] S21. According to Buckingham's theorem, charge, distance from the explosion center and vibration wave velocity are selected from multiple physical factors as basic dimensions. The three basic dimensions are independent and can form a full-rank matrix, that is, they can be combined into the dimensions of all other related physical quantities;

[0048] use To represent dimensionless, we can get:

[0049] ;

[0050] in, , , , , Respectively , , , , , represents a dimensionless function;

[0051] S22. Using the three basic dimensions, the remaining physical factors are converted into corresponding dimensionless parameters;

[0052] Specifically, The transformation process is as follows, assuming:

[0053] ;

[0054] In the formula, , , is the unknown coefficient, , , , Substituting the dimensionless into the above formula, we get:

[0055] ;

[0056] In order to make the above expression dimensionless, the dimensions of the numerator and denominator must be equal, that is:

[0057] ;

[0058] therefore, , , , and substitute it into In, get The dimensionless parameter is ;

[0059] Similarly, , , , , The dimensionless parameters are: , , , , ;

[0060] S23. Generate a dimensionless expression of blasting impact stress using each dimensionless parameter and the basic model;

[0061] Specifically, step S23 includes:

[0062] S231. Substitute the dimensionless parameters into the basic model to generate an initial dimensionless expression of the blasting impact stress:

[0063] ;

[0064] S232. Selecting key dimensionless parameters from each dimensionless parameter, taking the remaining dimensionless parameters as constants, simplifying the initial dimensionless expression of the blasting impact stress, and obtaining the dimensionless expression of the blasting impact stress;

[0065] In this embodiment, through sensitivity analysis, and As the key dimensionless parameter, the dimensionless expression of blasting impact stress is obtained as follows:

[0066] ;

[0067] S24. Use key dimensionless parameters in dimensionless expressions to construct new dimensionless parameters :

[0068] make ;

[0069] In the formula, and is the unknown coefficient. For a given particle in the tunnel surrounding rock, and is a constant.

[0070] S25. A general model of blasting impact stress is constructed by the new dimensionless parameters;

[0071] Specifically, the S25 includes:

[0072] S251. Establishing a first relationship between blasting impact stress and a new dimensionless parameter, performing polynomial rearrangement on the natural logarithm of the first relationship to obtain a logarithmic expression;

[0073] Specifically, due to and is a constant, and There is a functional relationship between them, namely:

[0074] ;

[0075] Taking ln of the formula, we get:

[0076] ;

[0077] In the formula, , are the coefficients obtained by function transformation;

[0078] After polynomial rearrangement of the formula, we get a logarithmic expression:

[0079] ;

[0080] In the formula, represents the effect of accompanying distance on the attenuation of explosion stress wave, is the decay index; and They represent the effects of particle vibration velocity and charge amount on the attenuation of explosion stress wave, and are the corresponding attenuation indexes respectively;

[0081] S252. Simplify the logarithmic expression to obtain a second relationship between blasting impact stress and charge amount, explosion center distance and vibration acceleration, that is, a general model of blasting impact stress;

[0082] Specifically, , , removing the logarithm from the logarithmic expression gives:

[0083] ;

[0084] In the formula, , , , , are all constants;

[0085] Command Parameters , , , , then the second relationship between blasting impact stress and charge amount, explosion center distance and vibration acceleration can be established:

[0086] ;

[0087] Based on the above embodiments, the method further includes:

[0088] S3. Obtain parameter information of the tunnel to be blasted, and establish a numerical analysis model of the tunnel according to the parameter information of the tunnel to be blasted, wherein the tunnel to be blasted is a double-hole single-line tunnel;

[0089] Specifically, step S3 includes:

[0090] S31. Obtain the size of the subsequent tunnel, the size of the preceding tunnel, the burial depth, the mechanical parameters of the surrounding rock and the geological conditions from the data of the tunnel to be blasted and the geological survey report;

[0091] S32. According to the size of the subsequent tunnel, the size of the preceding tunnel, the burial depth, the mechanical parameters of the surrounding rock and the geological conditions, the double-hole single-line tunnel model is established using ANSYS / LS-DYNA finite element software, and the fluid (such as air) and solid (such as surrounding rock, explosive) involved in the tunnel excavation process are coupled and analyzed;

[0092] The double-hole single-line tunnel model is divided into solid unit hexahedral grids; non-reflection boundary conditions are set on the front, back, left, right and bottom of the double-hole single-line tunnel model, and the top surface of the double-hole single-line tunnel model is set as a free boundary to prevent the blasting seismic wave from reflecting on the boundary.

[0093] S33. In the double-hole single-line tunnel model, surrounding rock, air and explosive materials are set respectively to obtain a numerical analysis model of the tunnel;

[0094] Specifically, the Johnson-Holmquist model is used for rock materials. This model is suitable for use under large strain, high strain rate and high pressure conditions. The equivalent strength of rock is related to pressure, strain rate and damage. The compressive strength and failure principal stress of rock are defined as rock failure criteria, so that when the compressive strength of rock reaches the set value or the unit principal stress reaches the set value, it fails, thereby simulating rock blasting failure. Its yield stress and strain rate The relationship is as follows:

[0095] ;

[0096] In the formula, is the initial yield stress, is the strain rate, C, is the Cowper Symonds strain rate parameter, is the effective plastic strain, is the plastic hardening modulus, is the pressure sensitivity coefficient.

[0097] Specifically, the explosive material is *MAT_HIGH_EXPLOSIVE_BURN, which is a material model keyword used in finite element analysis software such as LS-DYNA. The keyword is used to define the behavior of the explosive material, especially to simulate the combustion and explosion process of the explosive.

[0098] The detonation pressure is calculated using the JWL state equation:

[0099] ;

[0100] In the formula, is the detonation pressure, is the internal energy of the explosive detonation products, is the relative volume of the detonation products, , , , , are all property constants of the selected explosives.

[0101] Specifically, the air pressure state equation of the model is expressed as:

[0102] ;

[0103] In the formula, is the internal energy of the material, , , , , , , is a coefficient related to the properties of the gas. Specifically, = = = = =0, = =γ-1, where γ is the adiabatic index of an ideal gas, γ = 1.4;

[0104] is the dimensionless form of relative volume, , where V is the current relative volume of the gas.

[0105] Based on the above embodiments, Figure 2 As shown, the method also includes:

[0106] S4. Simulate the vibration acceleration data and blasting impact stress data of the surrounding rock mass points in adjacent tunnels under different charge amounts and different blast center distances in the numerical analysis model;

[0107] In this embodiment, simulation parameters need to be set during the simulation process, including:

[0108] Determine the charge range to be simulated, for example, from 100 kg to 500 kg, with a step size of 50 kg;

[0109] Determine the range of explosion center distance to be simulated, for example, from 1m to 10m, and simulate with a step length of 1m;

[0110] For each combination of charge amount and distance from the detonation center, corresponding simulation conditions are set, including the location, type and distribution of the explosives. The location of the explosives includes the cutout holes, auxiliary holes or peripheral holes of the subsequent tunnel.

[0111] Based on the above embodiments, the method further includes:

[0112] S5. Substituting the acceleration data and blasting impact stress data into the general model to obtain the blasting impact stress calculation model of the tunnel to be blasted;

[0113] Specifically, the acceleration data, the blasting impact stress data, the charge amount and the distance from the explosion center are substituted into the second relational expression for fitting to obtain the blasting impact stress calculation model.

[0114] Substitute the acceleration data, blasting impact stress data, charge amount and explosion center distance into the second relationship to solve , , The specific value of can be used to obtain the blasting impact stress calculation model.

[0115] Specifically, step S5 includes:

[0116] Based on the above embodiments, the method further includes:

[0117] S6. constructing an evaluation model of the probability of a time-delay rockburst in an adjacent tunnel when the tunnel to be blasted is induced by the blasting impact stress calculation model;

[0118] Specifically, step S6 includes:

[0119] S61. Obtaining the uniaxial compressive strength of the surrounding rock, obtaining the uniaxial compressive strength of the surrounding rock, and generating an evaluation model based on the ratio of the uniaxial compressive strength of the surrounding rock to the blasting impact stress calculation model;

[0120] ;

[0121] In the formula, represents the rockburst probability, It represents the uniaxial compressive strength of the surrounding rock. The uniaxial compressive strength of the surrounding rock needs to be obtained by sampling the surrounding rock at the engineering site, obtaining rock samples for blasting in adjacent tunnels before blasting, and conducting indoor uniaxial compression tests on the rock samples.

[0122] S62. Collect the charge amount, distance from the blasting center and acceleration of the tunnel to be blasted, substitute the charge amount, distance from the blasting center and acceleration into the evaluation model to calculate the probability of time-delay rockburst in adjacent tunnels:

[0123] when When it is greater than the rockburst assessment threshold N, there is no risk of rockburst in the adjacent tunnels. When it is less than the rockburst assessment threshold N, there is a risk of time-delay rockburst in the adjacent tunnel explosion detection.

[0124] Embodiment 2:

[0125] like Figure 2 As shown, this embodiment provides an evaluation system for time-delay rockburst induced by blasting disturbance, the system comprising:

[0126] Acquisition module: acquiring multiple physical factors related to blasting impact stress, and constructing a basic model of blasting impact stress using the multiple physical factors;

[0127] Selection module: selecting basic dimensions from multiple physical factors based on dimensional analysis, and constructing a general model of blasting impact stress using the basic dimensions and basic models;

[0128] Model building module: obtaining parameter information of the tunnel to be blasted, and building a numerical analysis model of the tunnel according to the parameter information of the tunnel to be blasted, wherein the tunnel to be blasted is a double-hole single-line tunnel;

[0129] Simulation module: simulates the vibration acceleration data and blasting impact stress data of the surrounding rock mass points in adjacent tunnels under different charge amounts and different blast center distances in the numerical analysis model;

[0130] Fitting module: using acceleration data and blasting impact stress data to substitute into the general model to fit the blasting impact stress calculation model of the tunnel to be blasted;

[0131] Calculation module: An evaluation model for the probability of inducing a time-delay rockburst in an adjacent tunnel when the tunnel to be blasted is blasted is constructed based on the blasting impact stress calculation model.

[0132] Based on the above embodiment, the selection module includes:

[0133] The first selection unit: According to Buckingham's theorem, the charge, the distance from the explosion center and the vibration wave speed are selected from multiple physical factors as basic dimensions, and a full-rank matrix is ​​formed by the three basic dimensions;

[0134] Conversion unit: using the three basic dimensions to convert the remaining physical factors into corresponding dimensionless parameters;

[0135] Generation unit: Generates dimensionless expression of blasting impact stress using dimensionless parameters and basic model;

[0136] Dimensionless parameter construction unit: construct new dimensionless parameters using key dimensionless parameters in dimensionless expressions;

[0137] General model building unit: building a general model of blasting impact stress by using the new dimensionless parameters.

[0138] Based on the above embodiment, the generating unit includes:

[0139] Substitution unit: Substituting the dimensionless parameters into the basic model to generate an initial dimensionless expression of the blasting impact stress;

[0140] The first simplification unit: select key dimensionless parameters from various dimensionless parameters, take the remaining dimensionless parameters as constants, simplify the initial dimensionless expression of blasting impact stress, and obtain the dimensionless expression of blasting impact stress.

[0141] Based on the above embodiments, the general model building unit includes:

[0142] A relationship establishing unit: establishing a first relationship between the blasting impact stress and the new dimensionless parameter, performing polynomial rearrangement on the natural logarithm of the first relationship to obtain a logarithmic expression;

[0143] The second simplification unit: simplifies the logarithmic expression to obtain a general model of blasting impact stress.

[0144] Based on the above embodiments, the model building module includes:

[0145] Acquisition unit: obtains the size of the subsequent tunnel, the size of the preceding tunnel, the burial depth, the mechanical parameters of the surrounding rock and the geological conditions from the data of the tunnel to be blasted and the geological survey report;

[0146] Mesh division unit: According to the size of the subsequent tunnel, the size of the preceding tunnel, the burial depth, the mechanical parameters of the surrounding rock and the geological conditions, a double-hole single-line tunnel model is established, and the double-hole single-line tunnel model is divided into solid unit hexahedral meshes;

[0147] Material setting unit: In the double-hole single-line tunnel model, the materials of surrounding rock, air and explosives are set respectively to obtain the numerical analysis model of the tunnel.

[0148] Based on the above embodiment, the calculation module includes:

[0149] The first calculation unit: obtains the uniaxial compressive strength of the surrounding rock, and generates an evaluation model based on the ratio of the uniaxial compressive strength of the surrounding rock to the blasting impact stress calculation model;

[0150] The second calculation unit: collects the charge amount, explosion center distance and acceleration when the tunnel to be blasted is blasted, and substitutes the charge amount, explosion center distance and acceleration into the evaluation model to calculate the probability of time-delay rockburst in adjacent tunnels.

[0151] It should be noted that, regarding the system in the above embodiment, the specific manner in which each module performs operations has been described in detail in the embodiment of the method, and will not be elaborated here.

[0152] Embodiment 3:

[0153] Corresponding to the above method embodiment, this embodiment also provides an evaluation device for blasting disturbance induced time-lag type rockburst. The evaluation device for blasting disturbance induced time-lag type rockburst described below and the evaluation method for blasting disturbance induced time-lag type rockburst described above can be referenced to each other.

[0154] Figure 3 FIG. 8 is a block diagram of an evaluation device 800 for blasting disturbance-induced time-delay rockburst according to an exemplary embodiment. Figure 3 As shown, the blasting disturbance induced time-delay rockburst assessment device 800 may include: a processor 801 and a memory 802. The blasting disturbance induced time-delay rockburst assessment device 800 may also include one or more of a multimedia component 803, an I / O interface 804, and a communication component 805.

[0155] The processor 801 is used to control the overall operation of the blasting disturbance induced time-delay type rockburst assessment device 800 to complete all or part of the steps in the above-mentioned blasting disturbance induced time-delay type rockburst assessment method. The memory 802 is used to store various types of data to support the operation of the blasting disturbance induced time-delay type rockburst assessment device 800, and these data may include, for example, instructions for any application or method operating on the blasting disturbance induced time-delay type rockburst assessment device 800, and application-related data, such as contact data, sent and received messages, pictures, audio, video, etc. The memory 802 can be implemented by any type of volatile or non-volatile storage device or a combination thereof, such as static random access memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic memory, flash memory, disk or optical disk. The multimedia component 803 may include a screen and an audio component. The screen may be, for example, a touch screen, and the audio component is used to output and / or input audio signals. For example, the audio component may include a microphone for receiving external audio signals. The received audio signal may be further stored in the memory 802 or sent via the communication component 805. The audio component also includes at least one speaker for outputting audio signals. The I / O interface 804 provides an interface between the processor 801 and other interface modules, and the above-mentioned other interface modules can be keyboards, mice, buttons, etc. These buttons can be virtual buttons or physical buttons. The communication component 805 is used for wired or wireless communication between the blasting disturbance induced time-delay rockburst evaluation device 800 and other devices. Wireless communication, such as Wi-Fi, Bluetooth, Near Field Communication (NFC), 2G, 3G or 4G, or a combination of one or more of them, so the corresponding communication component 805 may include: Wi-Fi module, Bluetooth module, NFC module.

[0156] In an exemplary embodiment, the evaluation device 800 for time-lag rockburst induced by blasting disturbance can be implemented by one or more application-specific integrated circuits (ASIC), digital signal processors (DSP), digital signal processing devices (DSPD), programmable logic devices (PLD), field programmable gate arrays (FPGA), controllers, microcontrollers, microprocessors or other electronic components to execute the above-mentioned evaluation method for time-lag rockburst induced by blasting disturbance.

[0157] In another exemplary embodiment, a computer-readable storage medium including program instructions is also provided, and when the program instructions are executed by a processor, the steps of the above-mentioned blasting disturbance induced time-delay type rockburst assessment method are implemented. For example, the computer-readable storage medium may be the above-mentioned memory 802 including program instructions, and the above-mentioned program instructions may be executed by the processor 801 of the blasting disturbance induced time-delay type rockburst assessment device 800 to complete the above-mentioned blasting disturbance induced time-delay type rockburst assessment method.

[0158] Embodiment 4:

[0159] Corresponding to the above method embodiment, a readable storage medium is also provided in this embodiment. The readable storage medium described below and the evaluation method for blasting disturbance-induced time-lag rockburst described above can refer to each other.

[0160] A readable storage medium stores a computer program, which, when executed by a processor, implements the steps of the method for evaluating a time-delay rockburst induced by blasting disturbance in the above method embodiment.

[0161] The readable storage medium may specifically be a U disk, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk, or other readable storage medium that can store program codes.

[0162] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. For those skilled in the art, the present invention may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

[0163] The above is only a specific embodiment of the present invention, but the protection scope of the present invention is not limited thereto. Any person skilled in the art can easily think of changes or substitutions within the technical scope disclosed by the present invention, which should be included in the protection scope of the present invention. Therefore, the protection scope of the present invention should be based on the protection scope of the claims.

Claims

1. A method for evaluating time-delay rockburst induced by blasting disturbance, characterized in that: include: Acquire multiple physical factors related to blasting impact stress, and construct a basic model of blasting impact stress using the multiple physical factors; Selecting basic dimensions from multiple physical factors based on dimensional analysis, and constructing a general model of blasting impact stress using the basic dimensions and basic models; Acquiring parameter information of a tunnel to be blasted, and establishing a numerical analysis model of the tunnel according to the parameter information of the tunnel to be blasted, wherein the tunnel to be blasted is a double-hole single-line tunnel; The vibration acceleration data and blasting impact stress data of surrounding rock mass points in adjacent tunnels were simulated in the numerical analysis model under different charge amounts and different blasting center distances; Substituting the acceleration data and the blasting impact stress data into the general model to obtain a blasting impact stress calculation model of the tunnel to be blasted; An evaluation model for the probability of inducing a time-delay rockburst in an adjacent tunnel when the tunnel to be blasted is blasted is constructed using the blasting impact stress calculation model; Among them, based on the dimensional analysis method, basic dimensions are selected from multiple physical factors, and the general model of blasting impact stress is constructed using the basic dimensions and basic models, including: According to Buckingham's theorem, the charge, the distance from the explosion center and the vibration wave speed are selected from multiple physical factors as basic dimensions, and a full-rank matrix is ​​formed by the three basic dimensions. The remaining physical factors are converted into corresponding dimensionless parameters respectively using the three basic dimensions; The dimensionless expression of blasting impact stress is generated by using dimensionless parameters and basic models; Construct new dimensionless parameters using key dimensionless parameters in dimensionless expressions; A general model of blasting impact stress is constructed by using the new dimensionless parameters; The general model of blasting impact stress is constructed by the new dimensionless parameters, including: Establishing a first relationship between the blasting impact stress and the new dimensionless parameter, performing polynomial rearrangement on the natural logarithm of the first relationship to obtain a logarithmic expression; The logarithmic expression is simplified to obtain a general model of blasting impact stress; The process of obtaining parameter information of the tunnel to be blasted and establishing a numerical analysis model of the tunnel according to the parameter information of the tunnel to be blasted includes: Obtain the size of the subsequent tunnel, the size of the preceding tunnel, the burial depth, the mechanical parameters of the surrounding rock and the geological conditions from the data of the tunnel to be blasted and the geological survey report; According to the size of the trailing tunnel, the size of the leading tunnel, the burial depth, the mechanical parameters of the surrounding rock and the geological conditions, a double-hole single-line tunnel model is established, and the double-hole single-line tunnel model is divided into solid unit hexahedral grids; In the double-hole single-line tunnel model, the materials of surrounding rock, air and explosives are set respectively to obtain the numerical analysis model of the tunnel.

2. The method for evaluating time-delay rockburst induced by blasting disturbance according to claim 1, characterized in that: Substituting each dimensionless parameter into the basic model to generate the dimensionless expression of blasting impact stress includes: Substituting the dimensionless parameters into the basic model to generate an initial dimensionless expression of the blasting impact stress; The key dimensionless parameters are selected from various dimensionless parameters, and the remaining dimensionless parameters are taken as constants. The initial dimensionless expression of blasting impact stress is simplified to obtain the dimensionless expression of blasting impact stress.

3. An evaluation system for time-delay rockburst induced by blasting disturbance, characterized in that: include: Acquisition module: acquiring multiple physical factors related to blasting impact stress, and constructing a basic model of blasting impact stress using the multiple physical factors; Selection module: selecting basic dimensions from multiple physical factors based on dimensional analysis, and generating a general model of blasting impact stress using the basic dimensions and basic models; Model building module: obtaining parameter information of the tunnel to be blasted, and building a numerical analysis model of the tunnel according to the parameter information of the tunnel to be blasted, wherein the tunnel to be blasted is a double-hole single-line tunnel; Simulation module: simulates the vibration acceleration data and blasting impact stress data of the surrounding rock mass points in adjacent tunnels under different charge amounts and different blast center distances in the numerical analysis model; Fitting module: using acceleration data and blasting impact stress data to substitute into the general model to fit the blasting impact stress calculation model of the tunnel to be blasted; Calculation module: constructing an evaluation model for the probability of inducing a time-delay rockburst in an adjacent tunnel when the tunnel to be blasted is blasted based on the blasting impact stress calculation model; Wherein, the selection module includes: The first selection unit: According to Buckingham's theorem, the charge, the distance from the explosion center and the vibration wave speed are selected from multiple physical factors as basic dimensions, and a full-rank matrix is ​​formed by the three basic dimensions; Conversion unit: using the three basic dimensions to convert the remaining physical factors into corresponding dimensionless parameters; Generation unit: Generates dimensionless expression of blasting impact stress using dimensionless parameters and basic model; Dimensionless parameter construction unit: construct new dimensionless parameters using key dimensionless parameters in dimensionless expressions; A general model building unit: building a general model of blasting impact stress by using the new dimensionless parameters; Wherein, the general model building unit includes: A relationship establishing unit: establishing a first relationship between the blasting impact stress and the new dimensionless parameter, performing polynomial rearrangement on the natural logarithm of the first relationship to obtain a logarithmic expression; The second simplification unit is used to simplify the logarithmic expression to obtain a general model of blasting impact stress; Wherein, the model building module includes: Acquisition unit: obtains the size of the subsequent tunnel, the size of the preceding tunnel, the burial depth, the mechanical parameters of the surrounding rock and the geological conditions from the data of the tunnel to be blasted and the geological survey report; Grid division unit: According to the size of the subsequent tunnel, the size of the preceding tunnel, the burial depth, the mechanical parameters of the surrounding rock and the geological conditions, Establish a double-hole single-line tunnel model and divide the double-hole single-line tunnel model into solid unit hexahedral meshes; Material setting unit: In the double-hole single-line tunnel model, the materials of surrounding rock, air and explosives are set respectively to obtain the numerical analysis model of the tunnel.

4. The blasting disturbance induced time-delay rockburst assessment system according to claim 3, characterized in that: The generating unit comprises: Substitution unit: Substituting the dimensionless parameters into the basic model to generate an initial dimensionless expression of the blasting impact stress; The first simplification unit: select key dimensionless parameters from various dimensionless parameters, take the remaining dimensionless parameters as constants, simplify the initial dimensionless expression of blasting impact stress, and obtain the dimensionless expression of blasting impact stress.

Images