Rockburst determination method, rockburst determination device and electronic device

Abstract

The present application provides a rockburst determination method, rockburst determination device, and electronic device. The rockburst determination method includes: obtaining the elastic energy density of a current region of a rock mass, and calculating rockburst parameters of the current region based at least on the elastic energy density of the current region; determining whether a rockburst has occurred in the current region based on the rockburst parameters of the current region; if a rockburst has occurred in the current region, obtaining the elastic energy density of a next region of the rock mass, and calculating rockburst parameters of the next region based at least on the elastic energy density of the next region; determining whether a rockburst has occurred in the next region based on the rockburst parameters of the next region; if a rockburst has occurred in the next region, repeating the calculation step and the second determination step at least once in sequence until no rockburst has occurred in the next region. This method solves the problem in the prior art of being unable to accurately and quantitatively predict rockbursts.

Description

Technical Field

[0001] The present application relates to the field of coal mines, and in particular, to a rockburst determination method, a rockburst determination device, a computer-readable storage medium, and an electronic device. Background Art

[0002] As underground mining deepens, the high ground stress in deep tunnels leads to an increased frequency of rock burst disasters. The rock burst dynamic disasters induced by deep engineering excavation are extremely destructive to underground engineering, causing serious consequences such as damage to construction equipment and casualties of construction personnel, and resulting in huge economic losses.

[0003] In order to provide a scientific basis and safe and reasonable theoretical and technical support that is more in line with the actual conditions of the engineering design and construction of deep underground projects, systematic research on the prediction and assessment methods of rockburst dynamic disasters has become an inevitable development trend and has important engineering significance. Summary of the Invention

[0004] The main purpose of the present application is to provide a rockburst determination method, a rockburst determination device, a computer-readable storage medium and an electronic device, so as to at least solve the problem in the prior art that rockburst cannot be accurately and quantitatively predicted.

[0005] According to one aspect of the present application, a rockburst determination method is provided, comprising: an acquisition step of acquiring an elastic energy density of a current region of a rock mass, and calculating rockburst parameters of the current region based at least on the elastic energy density of the current region, wherein the rock mass is a rock mass for underground engineering and includes multiple regions, the current region is any one of the multiple regions, and the elastic energy density is the elastic energy of the rock mass per unit volume; a first determination step of determining whether a rockburst has occurred in the current region based on the rockburst parameters of the current region; a calculation step of acquiring an elastic energy density of a next region of the rock mass if a rockburst has occurred in the current region, and calculating rockburst parameters of the next region based at least on the elastic energy density of the next region if a rockburst has occurred in the current region, wherein the next region is an adjacent region to the current region; a second determination step of determining whether a rockburst has occurred in the next region based on the rockburst parameters of the next region; and if a rockburst has occurred in the next region, repeating the calculation step and the second determination step in sequence at least once until a rockburst has not occurred in the next region, wherein during the repetitions, the next region of the current time is an adjacent region to the next region of the previous time.

[0006] Optionally, calculating the rockburst parameters of the corresponding area at least according to the elastic energy density of the corresponding area, wherein the corresponding area is the current area or the next area, includes: obtaining the constraint pressure of the corresponding area; according to the constraint pressure p and the formula Calculate the rockburst parameter URLERI of the corresponding area, where U i is the elastic energy density of the rock mass at the i-th calculation step, U i+1 is the elastic energy density of the rock mass in the (i+1)th calculation step, and dt is the time step for calculating the rockburst parameters of the corresponding area.

[0007] Optionally, determining whether a rockburst occurs in the corresponding area based on the rockburst parameters of the corresponding area, wherein the corresponding area is the current area or the next area, includes: when the rockburst parameters of the corresponding area are greater than 1, determining that a rockburst occurs in the corresponding area.

[0008] Optionally, in the case of rock burst in the current area, the elastic energy density of the next area of ​​the rock mass is obtained, including: according to the seismic moment M0 of the rock burst in the current area, the magnitude M of the rock burst in the current area, the amplitude-energy relationship coefficients α1 and α2, the stiffness μ of the fault medium where the rock burst occurs in the current area, the shear strength τ of the rock mass y , seismic radiation efficiency η and Calculate the elastic energy ε of the next area of ​​the rock mass T ; According to the elastic energy of the next area of ​​the rock mass, the elastic energy of the next area within a unit volume is calculated to obtain the elastic energy density of the next area of ​​the rock mass.

[0009] Optionally, the method also includes: obtaining engineering parameters of the underground project, wherein the engineering parameters include at least excavation method parameters and hard rock mechanical behavior parameters, the excavation method parameters are used to characterize the excavation method of the underground project, and the hard rock mechanical behavior parameters are used to characterize the hard rock mechanical behavior of the underground project; establishing an underground engineering model based on the engineering parameters; when it is determined that a rock burst occurs in the corresponding area, the time when the rock burst occurs in the corresponding area, the location where the rock burst occurs in the corresponding area, and the degree of damage caused by the rock burst in the corresponding area are obtained according to the underground engineering model and displayed on a display interface, wherein the corresponding area is the current area or the next area.

[0010] Optionally, obtaining the excavation method parameters includes: obtaining the excavation method of the underground project, the density of the explosive, the explosion speed of the explosive, the isentropic index of the explosive, the diameter of the explosion hole, the distance between the centers of two adjacent explosion holes, the depth of the explosion hole, the expansion speed of the cracks generated by the explosion, the propagation speed of the rarefaction wave in the explosion, and the average speed of the rarefaction wave reflection, wherein the rarefaction wave is generated by the gas from the top to the bottom of the explosion hole after the explosion; determining the excavation method of the underground project, when the excavation method of the underground project is the drilling and blasting method, according to the formula Calculate the equivalent explosion stress peak P, where ρ e is the density of the explosive, V d is the explosion velocity of the explosive, γ is the isentropic index of the explosive; when the excavation method of the underground project is mechanical method, according to the formula Calculate the equivalent explosion stress peak value P, where r b is the diameter of the explosion hole, S is the distance between the centers of two adjacent explosion holes; according to the formula Calculate the load rise time t r , where L is the depth of the explosion hole; according to the formula Calculate the load duration t d , where c f is the expansion speed of the cracks generated by the explosion, c u1 is the propagation speed of the rarefaction wave in the explosion, c u2 is the average velocity of the rarefaction wave reflection; according to the formula Calculate the dynamic unloading duration t du .

[0011] Optionally, obtaining the hard rock mechanical behavior parameters includes: according to the formula: E(κ, I1) = B1 + B2κ + B3κ 2 +B4ln(I1), calculate the elastic-plastic coupling eigenvalue E, where κ is the plastic internal variable, I1 is the first stress invariant, B1, B2, B3 and B4 are the first fitting parameters; when the plastic internal variable is less than or equal to the first threshold, according to the formula Calculating strain hardening eigenvalues Among them, χ1, ζ1 and ξ1 are the first calculation parameters; when the plastic internal variable is greater than the first threshold, according to the formula Calculate the strain hardening eigenvalue in, is the residual internal friction angle; when the plastic internal variable is less than or equal to the second threshold, according to the formula c(κ)=χ2κ 2+ζ2κ+ξ2, calculate the strain softening characteristic value c, where χ2, ζ2 and ξ2 are the second calculation parameters; when the plastic internal variable is greater than the second threshold and less than or equal to the third threshold, according to the formula c(κ)=χ3κ 2 +ζ3κ+ξ3, calculate the strain softening characteristic value c, where the third threshold is greater than the second threshold, χ3, ζ3 and ξ3 are third calculation parameters; when the plastic internal variable is greater than the third threshold, according to the formula c(κ)=c r , calculate the strain softening characteristic value c, where c r is the residual cohesion; according to the formula The dilution characteristic value ψ is calculated, where C1, C2, C3, C4, C5, and C6 are second fitting parameters.

[0012] According to another aspect of the present application, a rockburst determination device is provided, comprising: an acquisition unit, configured to acquire an elastic energy density of a current region of a rock mass, and calculate rockburst parameters of the current region based at least on the elastic energy density of the current region, wherein the rock mass is a rock mass for underground engineering and includes a plurality of regions, the current region is any one of the plurality of regions, and the elastic energy density is the elastic energy of the rock mass per unit volume; a first determination unit, configured to determine, in a first determination step, whether a rockburst occurs in the current region based on the rockburst parameters of the current region; and a calculation unit, configured to calculate, in the current region, whether a rockburst occurs in the current region. In the case that a rock burst occurs in a region, the elastic energy density of a next region of the rock mass is obtained, and the rock burst parameters of the next region are calculated based on at least the elastic energy density of the next region, wherein the next region is an adjacent region of the current region; a second determination unit is used for determining, in a second determination step, whether a rock burst occurs in the next region based on the rock burst parameters of the next region; and a repetition unit is used for, in the case that a rock burst occurs in the next region, sequentially repeating the calculation step and the second determination step at least once until no rock burst occurs in the next region, wherein, during the repetition process, the current next region is an adjacent region of the previous next region.

[0013] According to another aspect of the present application, a computer-readable storage medium is provided, wherein the computer-readable storage medium includes a stored program, wherein when the program is executed, the device where the computer-readable storage medium is located is controlled to execute any one of the methods described.

[0014] According to another aspect of the present application, an electronic device is provided, comprising a memory and a processor, wherein a computer program is stored in the memory, and the processor is configured to execute any one of the methods described above through the computer program.

[0015] Applying the technical solution of the present application, first, the elastic energy density of the current region of the rock mass is obtained, and the rockburst parameters of the current region are calculated based on at least the elastic energy density of the current region; then, based on the rockburst parameters of the current region, it is determined whether a rockburst has occurred in the current region; if a rockburst has occurred in the current region, the elastic energy density of the next region of the rock mass is obtained, and the rockburst parameters of the next region are calculated based on at least the elastic energy density of the next region; then, based on the rockburst parameters of the next region, it is determined whether a rockburst has occurred in the next region; finally, if a rockburst has occurred in the next region, the calculation step and the second determination step are repeated at least once in sequence until no rockburst has occurred in the next region. By calculating the rockburst parameters of the current region and the next region, the precise rockburst phenomenon can be quantitatively determined, especially in unconstructed areas, so that the risk of rockburst can be predicted in advance, solving the problem of the inability to accurately and quantitatively predict rockburst in the prior art. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] The drawings that constitute part of this application are used to provide a further understanding of this application. The illustrative embodiments of this application and their descriptions are used to explain this application and do not constitute an improper limitation on this application. In the drawings:

[0017] Figure 1 A hardware structure block diagram of a mobile terminal for executing a rockburst determination method provided in an embodiment of the present application is shown;

[0018] Figure 2 A schematic flow chart of a rockburst determination method provided according to an embodiment of the present application is shown;

[0019] Figure 3 A schematic diagram showing a rock burst result provided according to an embodiment of the present application is shown;

[0020] Figure 4 A schematic diagram of equivalent explosion stress provided in an embodiment of the present application is shown;

[0021] Figure 5 A schematic diagram showing a load rise time and load duration provided according to an embodiment of the present application is shown;

[0022] Figure 6 A schematic diagram showing a dynamic unloading duration provided according to an embodiment of the present application is shown;

[0023] Figure 7 A schematic diagram of the overall process of a rockburst determination method provided in accordance with an embodiment of the present application is shown;

[0024] Figure 8 A structural block diagram of a rockburst determination device provided according to an embodiment of the present application is shown.

[0025] The above drawings include the following reference numerals:

[0026] 102. Processor; 104. Memory; 106. Transmission device; 108. Input / output device. DETAILED DESCRIPTION

[0027] It should be noted that, in the absence of conflict, the embodiments and features of the embodiments in this application can be combined with each other. The present application will be described in detail below with reference to the accompanying drawings and in combination with the embodiments.

[0028] In order to enable those skilled in the art to better understand the present invention, the following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments in the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts should fall within the scope of protection of this application.

[0029] It should be noted that the terms "first", "second", etc. in the specification and claims of the present application and the above-mentioned drawings are used to distinguish similar objects and are not necessarily used to describe a specific order or sequential order. It should be understood that the data used in this way can be interchanged where appropriate, so that the embodiments of the present application described here. In addition, the terms "including" and "having" and any of their variations are intended to cover non-exclusive inclusions. For example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0030] As introduced in the background technology, it is impossible to achieve accurate quantitative prediction of rockburst in the existing technology. To solve the above problem, the embodiments of the present application provide a rockburst determination method, a rockburst determination device, a computer-readable storage medium and an electronic device.

[0031] The technical solutions in the embodiments of the present invention will be described clearly and completely below with reference to the accompanying drawings in the embodiments of the present invention.

[0032] The method embodiments provided in the embodiments of the present application can be executed in a mobile terminal, a computer terminal or a similar computing device. Taking running on a mobile terminal as an example, Figure 1 FIG. 1 is a hardware structure block diagram of a mobile terminal of a rockburst determination method according to an embodiment of the present invention. Figure 1 As shown, the mobile terminal may include one or more ( Figure 1Only one is shown) a processor 102 (the processor 102 may include but is not limited to a microprocessor MCU or a programmable logic device FPGA and other processing devices) and a memory 104 for storing data, wherein the mobile terminal may also include a transmission device 106 and an input and output device 108 for communication functions. It will be understood by those skilled in the art that Figure 1 The structure shown is only for illustration and does not limit the structure of the mobile terminal. Figure 1 More or fewer components than shown, or with Figure 1 Different configurations shown.

[0033] The memory 104 can be used to store computer programs, such as software programs and modules of application software, such as the computer program corresponding to the device information display method in the embodiment of the present invention. The processor 102 executes various functional applications and data processing by running the computer program stored in the memory 104, that is, implementing the above-mentioned method. The memory 104 may include a high-speed random access memory and may also include a non-volatile memory, such as one or more magnetic storage devices, flash memory, or other non-volatile solid-state memory. In some examples, the memory 104 may further include a memory remotely located relative to the processor 102, and these remote memories can be connected to the mobile terminal via a network. Examples of the above-mentioned networks include but are not limited to the Internet, corporate intranets, local area networks, mobile communication networks, and combinations thereof. The transmission device 106 is used to receive or send data via a network. Specific examples of the above-mentioned network may include a wireless network provided by the mobile terminal's communication provider. In one example, the transmission device 106 includes a network adapter (Network Interface Controller, abbreviated as NIC), which can be connected to other network devices via a base station to communicate with the Internet. In one example, the transmission device 106 may be a radio frequency (RF) module, which is used to communicate with the Internet wirelessly.

[0034] In this embodiment, a rockburst determination method running on a mobile terminal, a computer terminal, or a similar computing device is provided. It should be noted that the steps shown in the flowcharts of the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions, and although a logical sequence is shown in the flowcharts, in some cases, the steps shown or described can be executed in an order different from that shown.

[0035] Figure 2 FIG. 1 is a flow chart of a rockburst determination method according to an embodiment of the present application. Figure 2 As shown, the method includes the following steps:

[0036] Step S201, an acquisition step, obtaining an elastic energy density of a current region of a rock mass, and calculating rockburst parameters of the current region based at least on the elastic energy density of the current region, wherein the rock mass is a rock mass for an underground engineering project and includes multiple regions, the current region is any one of the multiple regions, and the elastic energy density is the elastic energy of the rock mass per unit volume;

[0037] Specifically, a rockburst is a dynamic crack expansion process in a rock mass accompanied by energy transfer and release. The critical crack expansion rate is used as the threshold for rockburst occurrence. When the energy release rate of the rock mass exceeds the threshold, a rockburst event is considered to have occurred.

[0038] Step S202, a first determination step, determining whether a rock burst occurs in the current area based on the rock burst parameters of the current area;

[0039] Specifically, the relative energy release rate of the rock mass per unit time and the elastic energy density of the above-mentioned current area are used as rockburst criterion to achieve a reasonable description of the rockburst process and further track the change law of the elastic energy density of the rock mass during the rockburst process.

[0040] Step S203, a calculation step, in which, when a rockburst occurs in the current area, elastic energy density of a next area of ​​the rock mass is obtained, and rockburst parameters of the next area are calculated based at least on the elastic energy density of the next area, wherein the next area is an adjacent area of ​​the current area;

[0041] Specifically, before a rockburst occurs, a considerable amount of elastic strain energy has already accumulated in the rock mass. During a rockburst, most of this strain energy is converted into dissipated energy, represented by the plastic properties and surface energy consumed by the rock mass. The remaining energy is converted into electromagnetic radiation energy, acoustic energy, and radiation energy represented by kinetic energy, and released to the outside world, namely the next area.

[0042] Step S204, a second determination step, determining whether a rock burst occurs in the next area according to the rock burst parameters of the next area;

[0043] Specifically, the calculation method of the rockburst parameters may be the same as the calculation method of the rockburst parameters of the current area.

[0044] Step S205: When rock burst occurs in the next area, repeat the calculation step and the second determination step at least once in sequence until no rock burst occurs in the next area, wherein, during the repetition process, the next area of ​​the current time is an adjacent area of ​​the next area of ​​the previous time.

[0045] Specifically, the basis for judging whether the rock burst does not occur in the next area may be that the surrounding rock of the rock mass maintains mechanical balance.

[0046] Through this embodiment, first, the elastic energy density of the current region of the rock mass is obtained, and the rockburst parameters of the current region are calculated based on at least the elastic energy density of the current region; then, based on the rockburst parameters of the current region, it is determined whether a rockburst has occurred in the current region; if a rockburst has occurred in the current region, the elastic energy density of the next region of the rock mass is obtained, and the rockburst parameters of the next region are calculated based on at least the elastic energy density of the next region; then, based on the rockburst parameters of the next region, it is determined whether a rockburst has occurred in the next region; finally, if a rockburst has occurred in the next region, the calculation step and the second determination step are repeated at least once in sequence until no rockburst has occurred in the next region. By calculating the rockburst parameters of the current region and the next region, the precise rockburst phenomenon can be quantitatively determined, especially in unconstructed areas, allowing for early prediction of the risk of rockburst, thereby resolving the problem of the inability to accurately and quantitatively predict rockburst in the prior art.

[0047] In a specific implementation process, the above steps S201 and S203 can be implemented by the following steps: Step S2011, obtaining the constraint pressure of the corresponding area; Step S2012, according to the constraint pressure p and the formula: Calculate the rockburst parameter URLERI of the corresponding area above, where U i is the elastic energy density of the rock mass at the i-th calculation step, U i+1 is the elastic energy density of the rock mass in the i+1th calculation step. This method can further calculate the precise elastic energy density value of the current region or the next region, and dt is the time step for calculating the rockburst parameters of the corresponding region.

[0048] Specifically, a rockburst is a dynamic crack expansion process in a rock mass accompanied by energy transfer and release. The critical crack expansion velocity is used as the threshold for rockburst occurrence. When the energy release rate of the rock mass exceeds the threshold, a rockburst event is considered to have occurred. Using the relative energy release rate of the rock mass per unit time as the rockburst criterion allows for a reasonable description of the rockburst process and further tracks the changes in the elastic energy density of the rock mass during the rockburst process.

[0049] In order to further accurately determine whether a rock burst has occurred, the above-mentioned steps S202 and S204 of the present application can be implemented by the following steps: Step S2021, when the rock burst parameter of the above-mentioned corresponding area is greater than 1, determine that a rock burst has occurred in the above-mentioned corresponding area.

[0050] Specifically, the rockburst parameter is realized by detecting the dynamic change of the elastic energy density of each yield element. When its value is greater than 1, it means that a rockburst has occurred in the element.

[0051] The above step S203 can also be implemented in other ways, for example: Step S2031, according to the earthquake moment M0 of the rock burst in the above current area, the magnitude M of the rock burst in the above current area, the amplitude-energy relationship coefficients α1 and α2, the stiffness μ of the fault medium of the rock burst in the above current area, the shear strength τ of the above rock mass, y , seismic radiation efficiency η and Calculate the elastic energy ε of the next area of ​​the above rock mass T Step S2032: Calculate the elastic energy of the next region per unit volume based on the elastic energy of the next region of the rock mass to obtain the elastic energy density of the next region of the rock mass. This method can quickly obtain the elastic energy density of the next region.

[0052] Specifically, before a rockburst occurs, a considerable amount of elastic strain energy has already accumulated in the rock mass. During a rockburst, most of the strain energy is converted into dissipated energy, represented by the plastic properties and surface energy consumed by the rock mass. The remaining energy is converted into electromagnetic radiation, acoustic energy, and radiation energy represented by kinetic energy, and released to the outside world. The impact of a rockburst on the surrounding rock mass primarily includes the release of radiation energy and the dynamic unloading of geostress at the new free surface. Due to the similarities between rockbursts and earthquakes, the ratio of rockburst radiation energy to total energy is further calculated to be equal to that of an earthquake. Finally, the seismic radiation rate is approximately 1% to 5%. Therefore, the radiation energy release process of a mesh element in the numerical calculation is equated with the explosive detonation process. After a rockburst occurs in a mesh element, it acts on the adjacent element nodes. Based on the stress continuity condition, the new boundary conditions during the rockburst process are further determined.

[0053] In some embodiments, the method further includes: step S206, obtaining engineering parameters of the underground project, wherein the engineering parameters include at least excavation method parameters and hard rock mechanical behavior parameters, the excavation method parameters are used to characterize the excavation method of the underground project, and the hard rock mechanical behavior parameters are used to characterize the hard rock mechanical behavior of the underground project; step S207, establishing an underground project model based on the engineering parameters; step S208, when it is determined that a rockburst has occurred in the corresponding area, obtaining the time when the rockburst occurred in the corresponding area, the location where the rockburst occurred in the corresponding area, and the degree of damage caused by the rockburst in the corresponding area based on the underground project model and displaying them on a display interface, wherein the corresponding area is the current area or the next area. This method can further obtain a three-dimensional numerical model of the underground project.

[0054] Specifically, the above engineering parameters may include parameters such as engineering geology, ground stress, and geological engineering structure. According to the above underground engineering model, the time when the rock burst occurred in the above corresponding area, the location where the rock burst occurred in the above corresponding area, and the degree of damage caused by the rock burst in the above corresponding area are obtained and displayed on the display interface, so that the following can be obtained: Figure 3 Schematic diagram of the simulation results shown.

[0055] The above step S207 can be specifically implemented by the following steps: Step S2071, obtaining the excavation method of the underground project, the density of the explosive, the explosion speed of the above explosive, the isentropic index of the above explosive, the diameter of the explosion hole, the distance between the centers of two adjacent explosion holes, the depth of the explosion hole, the expansion speed of the cracks generated by the explosion, the propagation speed of the rarefaction wave in the explosion, and the average speed of the reflection of the above rarefaction wave, wherein the above rarefaction wave is generated by the gas from the top to the bottom of the above explosion hole after the explosion; Step S2072, determining the excavation method of the above underground project. In the case where the excavation method of the above underground project is the drilling and blasting method, according to the formula Calculate the equivalent explosion stress peak P, where ρ e is the density of the above explosive, V d is the explosion velocity of the explosive, γ is the isentropic index of the explosive; Step S2073, when the excavation method of the underground project is mechanical method, according to the formula Calculate the equivalent explosion stress peak P, where r b is the diameter of the explosion hole, S is the distance between the centers of two adjacent explosion holes; Step S2074, according to the formula Calculate the load rise time t r , where L is the depth of the explosion hole; Step S2075, according to the formula Calculate the load duration t d , where c f is the expansion speed of the cracks generated by the explosion, c u1 is the propagation speed of the rarefaction wave in the explosion, c u2 is the average velocity of the rarefaction wave reflection; Step S2076, according to the formula Calculate the dynamic unloading duration t du This method can be further used to calculate the precise excavation method parameters.

[0056] Specifically, when excavating underground projects, the blast load pressure on the blast hole wall is physically measured during the explosion, taking into account the combined effects of blast load and dynamic unloading. Using multi-hole millisecond blasting technology, the blast stress near the blast hole and the excavation load near the boundary of the underground project excavation are treated equivalently. The stress wave generated by the explosion near the blast hole exhibits a negative power function attenuation law as the propagation distance increases, and the equivalent blast stress peak value at the excavation boundary is obtained, such as Figure 4 As shown, the horizontal axis represents time, including load rise time and load duration, and the vertical axis represents stress. Figure 5 As shown. In addition, rock blasting is an instantaneous and complex process. The entire process lasts only a few milliseconds. In the numerical analysis, it is assumed that the explosion load waveform is a triangular shock wave, and the load rise time and load duration of the blasting process are obtained; after excavation, a free surface is formed near the contour of the underground engineering. At the same time, the radial stress and shear stress of the surface surrounding rock are instantly zero. The boundaries of the radial and shear stress unloading values ​​of the foundation pit are calculated by numerical simulation during the dynamic process. The boundaries of the unbalanced force at each node of the excavation are directly determined during the excavation process, and the dynamic unloading of the ground stress is linear unloading. The dynamic unloading duration of the surrounding rock is approximately equal to the explosion fracture duration between adjacent blasting holes in the same blasting circle, and then the dynamic unloading duration of the blasting is determined, such as Figure 6 shown.

[0057] The above step S207 can also be implemented by the following steps: Step S2077, according to the formula: E(κ, I1) = B1 + B2κ + B3κ 2 +B4ln(I1), calculate the elastic-plastic coupling eigenvalue E, where κ is the plastic internal variable, I1 is the first stress invariant, B1, B2, B3 and B4 are the first fitting parameters; step S2078, when the above plastic internal variable is less than or equal to the first threshold, according to the formula Calculating strain hardening eigenvalues Wherein, χ1, ζ1 and ξ1 are the first calculation parameters; Step S2079, when the above plastic internal variable is greater than the above first threshold, according to the formula Calculate the above strain hardening eigenvalues in, is the residual internal friction angle; Step S2080, when the above plastic internal variable is less than or equal to the second threshold, according to the formula c(κ)=χ2κ 2 +ζ2κ+ξ2, calculate the strain softening characteristic value c, where χ2, ζ2 and ξ2 are the second calculation parameters; Step S2081, when the above-mentioned plastic internal variable is greater than the above-mentioned second threshold and less than or equal to the third threshold, according to the formula c(κ)=χ3κ 2+ζ3κ+ξ3, calculate the strain softening characteristic value c, where the third threshold is greater than the second threshold, χ3, ζ3 and ξ3 are third calculation parameters; step S2082, when the plastic internal variable is greater than the third threshold, according to the formula c(κ)=c r , calculate the strain softening characteristic value c, where c r is the remaining cohesion; step S2083, according to the formula: Calculate the dilution characteristic value ψ, where C1, C2, C3, C4, C5, and C6 are the second fitting parameters. This method can further calculate the accurate mechanical behavior parameters of hard rock.

[0058] Specifically, the accurate description of the mechanical behavior of high-stress hard rock under rockburst environment includes: significant elastic-plastic coupling characteristics, strain hardening, softening characteristics, and dilution characteristics during deformation. In the rockburst mechanics model, the relationship between the elastic modulus, the first stress invariant, and the plastic internal variable is used to characterize the elastic-plastic coupling characteristics of high-in-situ stress hard rock and the influence of strain rate on rock strength. The non-associative flow rule is further used to express the expansion characteristics of hard rock, and the relationship between the expansion angle, constraint pressure, and internal factors can be plotted. In the rockburst mechanics model, the strain hardening and softening characteristics are further reflected through the relationship between the yield function and the plastic internal variable. In the rockburst mechanics model, the strain hardening and softening characteristics are further reflected through the relationship between the yield function and the plastic internal variable.

[0059] In order to enable those skilled in the art to more clearly understand the technical solution of the present application, the implementation process of the rockburst determination method of the present application will be described in detail below with reference to specific embodiments.

[0060] This embodiment relates to a specific rockburst determination method, such as Figure 7 As shown, the following steps are included:

[0061] Step S1: In order to establish a rockburst process model and calculation parameters, facilitate dynamic calculation of the rockburst process simulation, and collect and organize underground engineering condition information;

[0062] Step S2: Adopting appropriate excavation methods in response to different underground engineering environments, including mechanical excavation or drilling;

[0063] Step S3: Selecting a rockburst mechanical model and calculation parameters;

[0064] Step S4: dynamically calculating the parameters;

[0065] Step S5: Using the rockburst criterion to determine whether a rockburst has occurred, if a rockburst has occurred, determining the new boundary conditions after the rockburst, and returning to step S4 to perform dynamic calculation again;

[0066] Step S6: If no rock burst occurs, determine whether the surrounding rock has reached a mechanical equilibrium state. If not, return to step S4 to continue dynamic calculation.

[0067] Step S7: When the equilibrium state is reached, the simulation ends.

[0068] The embodiments of the present application also provide a rockburst determination device. It should be noted that the rockburst determination device of the embodiments of the present application can be used to execute the rockburst determination method provided in the embodiments of the present application. The device is used to implement the above-mentioned embodiments and preferred implementation methods, and the details that have been described will not be repeated. As used below, the term "module" can be a combination of software and / or hardware that implements a predetermined function. Although the devices described in the following embodiments are preferably implemented in software, implementation in hardware, or a combination of software and hardware, is also possible and contemplated.

[0069] The rockburst determination device provided in the embodiments of the present application is introduced below.

[0070] Figure 8 Schematic diagram of a rockburst determination device according to an embodiment of the present application. Figure 8 As shown, the device includes:

[0071] an acquisition unit 10, configured to acquire, in an acquisition step, an elastic energy density of a current region of the rock mass, and calculate rockburst parameters of the current region based at least on the elastic energy density of the current region, wherein the rock mass is a rock mass for underground engineering and includes multiple regions, the current region is any one of the multiple regions, and the elastic energy density is the elastic energy of the rock mass per unit volume;

[0072] Specifically, a rockburst is a dynamic crack expansion process in a rock mass accompanied by energy transfer and release. The critical crack expansion rate is used as the threshold for rockburst occurrence. When the energy release rate of the rock mass exceeds the threshold, a rockburst event is considered to have occurred.

[0073] A first determining unit 20 is configured to determine, in a first determining step, whether a rockburst occurs in the current area based on the rockburst parameters of the current area;

[0074] Specifically, the relative energy release rate of the rock mass per unit time and the elastic energy density of the above-mentioned current area are used as rockburst criterion to achieve a reasonable description of the rockburst process and further track the change law of the elastic energy density of the rock mass during the rockburst process.

[0075] a calculation unit 30 configured to, in a calculation step, obtain, when a rockburst occurs in the current region, an elastic energy density of a next region of the rock mass, and calculate rockburst parameters of the next region based at least on the elastic energy density of the next region, wherein the next region is an adjacent region to the current region;

[0076] Specifically, before a rockburst occurs, a considerable amount of elastic strain energy has already accumulated in the rock mass. During a rockburst, most of this strain energy is converted into dissipated energy, represented by the plastic properties and surface energy consumed by the rock mass. The remaining energy is converted into electromagnetic radiation energy, acoustic energy, and radiation energy represented by kinetic energy, and released to the outside world, namely the next area.

[0077] A second determining unit 40 is used for the second determining step of determining whether a rock burst occurs in the next area according to the rock burst parameters of the next area;

[0078] Specifically, the calculation method of the rockburst parameters may be the same as the calculation method of the rockburst parameters of the current area.

[0079] The repeating unit 50 is used to repeat the above calculation step and the above second determination step at least once in sequence when a rock burst occurs in the above next area, until no rock burst occurs in the above next area, wherein, during the repetition process, the current next area is an adjacent area of ​​the previous next area.

[0080] Specifically, the basis for judging whether the rock burst does not occur in the next area may be that the surrounding rock of the rock mass maintains mechanical balance.

[0081] According to this embodiment, the acquisition unit acquires the elastic energy density of the current region of the rock mass and calculates the rockburst parameters of the current region based at least on the elastic energy density of the current region; the first determination unit determines whether a rockburst has occurred in the current region based on the rockburst parameters of the current region; the calculation unit acquires the elastic energy density of the next region of the rock mass if a rockburst has occurred in the current region and calculates the rockburst parameters of the next region based at least on the elastic energy density of the next region; the second determination unit determines whether a rockburst has occurred in the next region based on the rockburst parameters of the next region; and the repetition unit repeats the calculation step and the second determination step in sequence at least once if a rockburst has occurred in the next region until no rockburst has occurred in the next region. By calculating the rockburst parameters of the current region and the next region, the precise rockburst phenomenon can be quantitatively determined, especially in unconstructed areas, allowing for early prediction of the risk of rockburst, thereby resolving the problem of the inability to accurately and quantitatively predict rockburst in the prior art.

[0082] In a specific implementation process, the acquisition unit or calculation unit includes a first acquisition module and a first calculation module, wherein the first acquisition module is used to obtain the constraint pressure of the corresponding area; the first calculation module is used to calculate the constraint pressure p according to the constraint pressure p and the formula Calculate the rockburst parameter URLERI of the corresponding area above, where U i is the elastic energy density of the rock mass at the i-th calculation step, U i+1 is the elastic energy density of the rock mass in the i+1th calculation step, and dt is the time step for calculating the rockburst parameters of the corresponding area. The device can further calculate the accurate elastic energy density value of the current area or the next area.

[0083] Specifically, a rockburst is a dynamic crack expansion process in a rock mass accompanied by energy transfer and release. The critical crack expansion velocity is used as the threshold for rockburst occurrence. When the energy release rate of the rock mass exceeds the threshold, a rockburst event is considered to have occurred. Using the relative energy release rate of the rock mass per unit time as the rockburst criterion allows for a reasonable description of the rockburst process and further tracks the changes in the elastic energy density of the rock mass during the rockburst process.

[0084] In order to further achieve accurate judgment of whether rock burst occurs, the first determination module and the second determination module of the present application are respectively used to determine that rock burst occurs in the corresponding area when the rock burst parameter of the corresponding area is greater than 1.

[0085] Specifically, the rockburst parameter is realized by detecting the dynamic change of the elastic energy density of each yield element. When its value is greater than 1, it means that a rockburst has occurred in the element.

[0086] The calculation unit includes a second calculation module and a third calculation module, wherein the second calculation module is used to calculate the rock burst moment M0 of the current region, the magnitude M of the rock burst, the amplitude-energy relationship coefficients α1 and α2, the stiffness μ of the fault medium of the rock burst, the shear strength τ of the rock mass, and the rock burst magnitude M of the current region. y , seismic radiation efficiency η and Calculate the elastic energy ε of the next area of ​​the above rock mass T The third calculation module is used to calculate the elastic energy of the next region per unit volume based on the elastic energy of the next region of the rock mass, thereby obtaining the elastic energy density of the next region of the rock mass. This device can quickly obtain the elastic energy density of the next region.

[0087] Specifically, before a rockburst occurs, a considerable amount of elastic strain energy has already accumulated in the rock mass. During a rockburst, most of the strain energy is converted into dissipated energy, represented by the plastic properties and surface energy consumed by the rock mass. The remaining energy is converted into electromagnetic radiation, acoustic energy, and radiation energy represented by kinetic energy, and released to the outside world. The impact of a rockburst on the surrounding rock mass primarily includes the release of radiation energy and the dynamic unloading of geostress at the new free surface. Due to the similarities between rockbursts and earthquakes, the ratio of rockburst radiation energy to total energy is further calculated to be equal to that of an earthquake. Finally, the seismic radiation rate is approximately 1% to 5%. Therefore, the radiation energy release process of a mesh element in the numerical calculation is equated with the explosive detonation process. After a rockburst occurs in a mesh element, it acts on the adjacent element nodes. Based on the stress continuity condition, the new boundary conditions during the rockburst process are further determined.

[0088] In some embodiments, the device further includes a processing unit, an establishment unit, and a display unit, wherein the processing unit is used to obtain engineering parameters of the underground project, wherein the engineering parameters include at least excavation method parameters and hard rock mechanical behavior parameters, wherein the excavation method parameters are used to characterize the excavation method of the underground project, and the hard rock mechanical behavior parameters are used to characterize the hard rock mechanical behavior of the underground project; the establishment unit is used to establish an underground project model based on the engineering parameters; and the display unit is used to, when it is determined that a rock burst has occurred in the corresponding area, obtain the time when the rock burst has occurred in the corresponding area, the location of the rock burst in the corresponding area, and the degree of damage caused by the rock burst in the corresponding area based on the underground project model and display them on a display interface, wherein the corresponding area is the current area or the next area. The device can further obtain a three-dimensional numerical model of the underground project.

[0089] Specifically, the above engineering parameters may include parameters such as engineering geology, ground stress, and geological engineering structure. According to the above underground engineering model, the time when the rock burst occurred in the above corresponding area, the location where the rock burst occurred in the above corresponding area, and the degree of damage caused by the rock burst in the above corresponding area are obtained and displayed on the display interface, so that the following can be obtained: Figure 3 Schematic diagram of the simulation results shown.

[0090] The processing unit includes a second acquisition module, a fourth calculation module, a fifth calculation module, a sixth calculation module, a seventh calculation module and an eighth calculation module, wherein the second acquisition module is used to obtain the excavation method of the underground project, the density of the explosive, the explosion speed of the explosive, the isentropic index of the explosive, the diameter of the explosion hole, the distance between the centers of two adjacent explosion holes, the depth of the explosion hole, the expansion speed of the cracks generated by the explosion, the propagation speed of the rarefaction wave in the explosion, and the average speed of the reflection of the rarefaction wave, wherein the rarefaction wave is generated by the gas from the top to the bottom of the explosion hole after the explosion; the fourth calculation module is used to determine the excavation method of the underground project. When the excavation method of the underground project is the drilling and blasting method, according to the formula Calculate the equivalent explosion stress peak P, where ρ e is the density of the above explosive, V d is the explosion velocity of the explosive, γ is the isentropic index of the explosive; the fifth calculation module is used for calculating the excavation method of the underground engineering according to the formula when the excavation method of the underground engineering is the mechanical method. Calculate the equivalent explosion stress peak P, where r b is the diameter of the explosion hole, S is the distance between the centers of two adjacent explosion holes; the sixth calculation module is used to calculate according to the formula Calculate the load rise time t r , where L is the depth of the explosion hole; the seventh calculation module is used to calculate the explosion hole according to the formula Calculate the load duration t d , where c f is the expansion speed of the cracks generated by the explosion, c u1 is the propagation speed of the rarefaction wave in the explosion, c u2 is the average velocity of the above-mentioned rarefaction wave reflection; the eighth calculation module is used to calculate the average velocity of the rarefaction wave reflection according to the formula Calculate the dynamic unloading duration t du The device can further calculate the precise excavation method parameters.

[0091] Specifically, when excavating underground projects, the blast load pressure on the blast hole wall is physically measured during the explosion, taking into account the combined effects of blast load and dynamic unloading. Using multi-hole millisecond blasting technology, the blast stress near the blast hole and the excavation load near the boundary of the underground project excavation are treated equivalently. The stress wave generated by the explosion near the blast hole exhibits a negative power function attenuation law as the propagation distance increases, and the equivalent blast stress peak value at the excavation boundary is obtained, such as Figure 4 As shown, the horizontal axis represents time, including load rise time and load duration, and the vertical axis represents stress. Figure 5As shown. In addition, rock blasting is an instantaneous and complex process. The entire process lasts only a few milliseconds. In the numerical analysis, it is assumed that the explosion load waveform is a triangular shock wave, and the load rise time and load duration of the blasting process are obtained; after excavation, a free surface is formed near the contour of the underground engineering. At the same time, the radial stress and shear stress of the surface surrounding rock are instantly zero. The boundaries of the radial and shear stress unloading values ​​of the foundation pit are calculated by numerical simulation during the dynamic process. The boundaries of the unbalanced force at each node of the excavation are directly determined during the excavation process, and the dynamic unloading of the ground stress is linear unloading. The dynamic unloading duration of the surrounding rock is approximately equal to the explosion fracture duration between adjacent blasting holes in the same blasting circle, and then the dynamic unloading duration of the blasting is determined, such as Figure 6 shown.

[0092] The processing unit further includes a ninth calculation module, a tenth calculation module, an eleventh calculation module, a twelfth calculation module, a thirteenth calculation module, a fourteenth calculation module, and a fifteenth calculation module, wherein the ninth calculation module is configured to calculate the value of the first and second kappa ... 2 +B4ln(I1), calculate the elastic-plastic coupling eigenvalue E, where κ is the plastic internal variable, I1 is the first stress invariant, B1, B2, B3 and B4 are the first fitting parameters; the tenth calculation module is used to calculate the elastic-plastic coupling eigenvalue E according to the formula when the above plastic internal variable is less than or equal to the first threshold. Calculating strain hardening eigenvalues Among them, χ1, ζ1 and ξ1 are the first calculation parameters; the eleventh calculation module is used to calculate the plastic internal variable according to the formula when it is greater than the first threshold. Calculate the above strain hardening eigenvalues in, is the residual internal friction angle; the twelfth calculation module is used for, when the above plastic internal variable is less than or equal to the second threshold, according to the formula c(κ)=χ2κ 2 +ζ2κ+ξ2, calculate the strain softening characteristic value c, where χ2, ζ2 and ξ2 are the second calculation parameters; the thirteenth calculation module is used to calculate the strain softening characteristic value c according to the formula c(κ)=χ3κ when the above-mentioned plastic internal variable is greater than the above-mentioned second threshold and less than or equal to the third threshold. 2 +ζ3κ+ξ3, calculate the strain softening characteristic value c, wherein the above-mentioned third threshold value is greater than the above-mentioned second threshold value, χ3, ζ3 and ξ3 are third calculation parameters; the fourteenth calculation module is used to calculate the strain softening characteristic value c according to the formula c(κ)=c when the above-mentioned plastic internal variable is greater than the above-mentioned third threshold value. r , calculate the strain softening characteristic value c, where c r is the residual cohesion; the fifteenth calculation module is used according to the formula: The dilution characteristic value ψ is calculated, where C1, C2, C3, C4, C5, and C6 are the second fitting parameters. The device can further calculate the accurate mechanical behavior parameters of hard rock.

[0093] Specifically, the accurate description of the mechanical behavior of high-stress hard rock under rockburst environment includes: significant elastic-plastic coupling characteristics, strain hardening, softening characteristics, and dilution characteristics during deformation. In the rockburst mechanics model, the relationship between the elastic modulus, the first stress invariant, and the plastic internal variable is used to characterize the elastic-plastic coupling characteristics of high-in-situ stress hard rock and the influence of strain rate on rock strength. The non-associative flow rule is further used to express the expansion characteristics of hard rock, and the relationship between the expansion angle, constraint pressure, and internal factors can be plotted. In the rockburst mechanics model, the strain hardening and softening characteristics are further reflected through the relationship between the yield function and the plastic internal variable. In the rockburst mechanics model, the strain hardening and softening characteristics are further reflected through the relationship between the yield function and the plastic internal variable.

[0094] The rockburst determination device includes a processor and memory. The acquisition unit, first determination unit, calculation unit, second determination unit, and repetition unit are all stored as program units in the memory. The processor executes the program units stored in the memory to implement the corresponding functions. All of the modules are located in the same processor; alternatively, the modules can be located in different processors in any combination.

[0095] The processor contains a kernel, which retrieves the corresponding program unit from the memory. One or more kernels can be set, and the presence of a rock burst can be determined by adjusting kernel parameters.

[0096] The memory may include non-permanent memory in a computer-readable medium, random access memory (RAM) and / or non-volatile memory, such as read-only memory (ROM) or flash RAM, and the memory includes at least one memory chip.

[0097] An embodiment of the present invention provides a computer-readable storage medium, which includes a stored program. When the program is executed, the device where the computer-readable storage medium is located is controlled to execute the rockburst determination method.

[0098] Specifically, the rockburst determination method includes:

[0099] Step S201, an acquisition step, obtaining an elastic energy density of a current region of a rock mass, and calculating rockburst parameters of the current region based at least on the elastic energy density of the current region, wherein the rock mass is a rock mass for an underground engineering project and includes multiple regions, the current region is any one of the multiple regions, and the elastic energy density is the elastic energy of the rock mass per unit volume;

[0100] Specifically, a rockburst is a dynamic crack expansion process in a rock mass accompanied by energy transfer and release. The critical crack expansion rate is used as the threshold for rockburst occurrence. When the energy release rate of the rock mass exceeds the threshold, a rockburst event is considered to have occurred.

[0101] Step S202, a first determination step, determining whether a rock burst occurs in the current area based on the rock burst parameters of the current area;

[0102] Specifically, the relative energy release rate of the rock mass per unit time and the elastic energy density of the above-mentioned current area are used as rockburst criterion to achieve a reasonable description of the rockburst process and further track the change law of the elastic energy density of the rock mass during the rockburst process.

[0103] Step S203, a calculation step, in which, when a rockburst occurs in the current area, elastic energy density of a next area of ​​the rock mass is obtained, and rockburst parameters of the next area are calculated based at least on the elastic energy density of the next area, wherein the next area is an adjacent area of ​​the current area;

[0104] Specifically, before a rockburst occurs, a considerable amount of elastic strain energy has already accumulated in the rock mass. During a rockburst, most of this strain energy is converted into dissipated energy, represented by the plastic properties and surface energy consumed by the rock mass. The remaining energy is converted into electromagnetic radiation energy, acoustic energy, and radiation energy represented by kinetic energy, and released to the outside world, namely the next area.

[0105] Step S204, a second determination step, determining whether a rock burst occurs in the next area according to the rock burst parameters of the next area;

[0106] Specifically, the calculation method of the rockburst parameters may be the same as the calculation method of the rockburst parameters of the current area.

[0107] Step S205: When rock burst occurs in the next area, repeat the calculation step and the second determination step at least once in sequence until no rock burst occurs in the next area, wherein, during the repetition process, the next area of ​​the current time is an adjacent area of ​​the next area of ​​the previous time.

[0108] Specifically, the basis for judging whether the rock burst does not occur in the next area may be that the surrounding rock of the rock mass maintains mechanical balance.

[0109] An embodiment of the present invention provides a processor, which is used to run a program, wherein the rockburst determination method is executed when the program is run.

[0110] Specifically, the rockburst determination method includes:

[0111] Step S201, an acquisition step, obtaining an elastic energy density of a current region of a rock mass, and calculating rockburst parameters of the current region based at least on the elastic energy density of the current region, wherein the rock mass is a rock mass for an underground engineering project and includes multiple regions, the current region is any one of the multiple regions, and the elastic energy density is the elastic energy of the rock mass per unit volume;

[0112] Specifically, a rockburst is a dynamic crack expansion process in a rock mass accompanied by energy transfer and release. The critical crack expansion rate is used as the threshold for rockburst occurrence. When the energy release rate of the rock mass exceeds the threshold, a rockburst event is considered to have occurred.

[0113] Step S202, a first determination step, determining whether a rock burst occurs in the current area based on the rock burst parameters of the current area;

[0114] Specifically, the relative energy release rate of the rock mass per unit time and the elastic energy density of the above-mentioned current area are used as rockburst criterion to achieve a reasonable description of the rockburst process and further track the change law of the elastic energy density of the rock mass during the rockburst process.

[0115] Step S203, a calculation step, in which, when a rockburst occurs in the current area, elastic energy density of a next area of ​​the rock mass is obtained, and rockburst parameters of the next area are calculated based at least on the elastic energy density of the next area, wherein the next area is an adjacent area of ​​the current area;

[0116] Specifically, before a rockburst occurs, a considerable amount of elastic strain energy has already accumulated in the rock mass. During a rockburst, most of this strain energy is converted into dissipated energy, represented by the plastic properties and surface energy consumed by the rock mass. The remaining energy is converted into electromagnetic radiation energy, acoustic energy, and radiation energy represented by kinetic energy, and released to the outside world, namely the next area.

[0117] Step S204, a second determination step, determining whether a rock burst occurs in the next area according to the rock burst parameters of the next area;

[0118] Specifically, the calculation method of the rockburst parameters may be the same as the calculation method of the rockburst parameters of the current area.

[0119] Step S205: When rock burst occurs in the next area, repeat the calculation step and the second determination step at least once in sequence until no rock burst occurs in the next area, wherein, during the repetition process, the next area of ​​the current time is an adjacent area of ​​the next area of ​​the previous time.

[0120] Specifically, the basis for judging whether the rock burst does not occur in the next area may be that the surrounding rock of the rock mass maintains mechanical balance.

[0121] An embodiment of the present invention provides a device, comprising a processor, a memory, and a program stored in the memory and executable on the processor. When the processor executes the program, at least the following steps are performed:

[0122] Step S201, an acquisition step, obtaining an elastic energy density of a current region of a rock mass, and calculating rockburst parameters of the current region based at least on the elastic energy density of the current region, wherein the rock mass is a rock mass for an underground engineering project and includes multiple regions, the current region is any one of the multiple regions, and the elastic energy density is the elastic energy of the rock mass per unit volume;

[0123] Step S202, a first determination step, determining whether a rock burst occurs in the current area based on the rock burst parameters of the current area;

[0124] Step S203, a calculation step, in which, when a rockburst occurs in the current area, elastic energy density of a next area of ​​the rock mass is obtained, and rockburst parameters of the next area are calculated based at least on the elastic energy density of the next area, wherein the next area is an adjacent area of ​​the current area;

[0125] Step S204, a second determination step, determining whether a rock burst occurs in the next area according to the rock burst parameters of the next area;

[0126] Step S205: When rock burst occurs in the next area, repeat the calculation step and the second determination step at least once in sequence until no rock burst occurs in the next area, wherein, during the repetition process, the next area of ​​the current time is an adjacent area of ​​the next area of ​​the previous time.

[0127] The devices in this article can be servers, PCs, PADs, mobile phones, etc.

[0128] The present application also provides a computer program product, which, when executed on a data processing device, is adapted to execute a program for initializing at least the following method steps:

[0129] Step S201, an acquisition step, obtaining an elastic energy density of a current region of a rock mass, and calculating rockburst parameters of the current region based at least on the elastic energy density of the current region, wherein the rock mass is a rock mass for an underground engineering project and includes multiple regions, the current region is any one of the multiple regions, and the elastic energy density is the elastic energy of the rock mass per unit volume;

[0130] Step S202, a first determination step, determining whether a rock burst occurs in the current area based on the rock burst parameters of the current area;

[0131] Step S203, a calculation step, in which, when a rockburst occurs in the current area, elastic energy density of a next area of ​​the rock mass is obtained, and rockburst parameters of the next area are calculated based at least on the elastic energy density of the next area, wherein the next area is an adjacent area of ​​the current area;

[0132] Step S204, a second determination step, determining whether a rock burst occurs in the next area according to the rock burst parameters of the next area;

[0133] Step S205: When rock burst occurs in the next area, repeat the calculation step and the second determination step at least once in sequence until no rock burst occurs in the next area, wherein, during the repetition process, the next area of ​​the current time is an adjacent area of ​​the next area of ​​the previous time.

[0134] Obviously, those skilled in the art will appreciate that the various modules or steps of the present invention described above can be implemented using a general-purpose computing device, can be centralized on a single computing device, or can be distributed across a network of multiple computing devices. They can be implemented using program code executable by the computing device, and thus, can be stored in a storage device and executed by the computing device. In some cases, the steps shown or described herein can be performed in a different order than that shown, or can be fabricated as separate integrated circuit modules, or multiple modules or steps can be fabricated as a single integrated circuit module. Thus, the present invention is not limited to any particular combination of hardware and software.

[0135] Those skilled in the art will appreciate that the embodiments of the present application can be provided as methods, systems, or computer program products. Therefore, the present application can adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application can adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to magnetic disk storage, CD-ROM, optical storage, etc.) that contain computer-usable program code.

[0136] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 a process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0137] These computer program instructions may also be stored in a computer readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 a process or multiple processes and / or boxes Figure 1 The function specified in one or more boxes.

[0138] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operational steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing the instructions executed on the computer or other programmable device for implementing the process. Figure 1 a process or multiple processes and / or boxes Figure 1 A step that specifies a function in one or more boxes.

[0139] In a typical configuration, a computing device includes one or more processors (CPUs), input / output interfaces, network interfaces, and memory.

[0140] The memory may include non-permanent memory in a computer-readable medium, random access memory (RAM) and / or non-volatile memory in the form of read-only memory (ROM) or flash RAM. The memory is an example of a computer-readable medium.

[0141] Computer-readable media includes permanent and non-permanent, removable and non-removable media that can be implemented by any method or technology to store information. The information can be computer-readable instructions, data structures, program modules or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disc read-only memory (CD-ROM), digital versatile disc (DVD) or other optical storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices or any other non-transmission media that can be used to store information that can be accessed by a computing device. As defined herein, computer-readable media does not include transitory computer-readable media (transitory media), such as modulated data signals and carrier waves.

[0142] It should also be noted that the terms "comprises," "includes," or any other variations thereof are intended to encompass non-exclusive inclusion, such that a process, method, commodity, or apparatus that includes a series of elements includes not only those elements but also other elements not explicitly listed, or includes elements inherent to such process, method, commodity, or apparatus. In the absence of further limitations, an element defined by the phrase "comprises a ..." does not exclude the presence of other identical elements in the process, method, commodity, or apparatus that includes the element.

[0143] From the above description, it can be seen that the above embodiments of the present application achieve the following technical effects:

[0144] 1) The rockburst determination method of the present application first obtains the elastic energy density of the current area of ​​the rock mass, and calculates the rockburst parameters of the current area based on at least the elastic energy density of the current area; then, based on the rockburst parameters of the current area, determines whether a rockburst has occurred in the current area; in the event that a rockburst has occurred in the current area, obtains the elastic energy density of the next area of ​​the rock mass, and calculates the rockburst parameters of the next area based on at least the elastic energy density of the next area; then, based on the rockburst parameters of the next area, determines whether a rockburst has occurred in the next area; finally, in the event that a rockburst has occurred in the next area, repeats the calculation step and the second determination step at least once in sequence until no rockburst occurs in the next area. By calculating the rockburst parameters of the current area and the next area, the precise rockburst phenomenon can be quantitatively determined, especially in unconstructed areas, and the risk of rockburst can be predicted in advance, solving the problem that the existing technology cannot accurately and quantitatively predict rockburst.

[0145] 2) The rockburst determination device of the present application comprises an acquisition unit that acquires the elastic energy density of the current region of the rock mass and calculates the rockburst parameters of the current region at least based on the elastic energy density of the current region; a first determination unit determines whether a rockburst has occurred in the current region based on the rockburst parameters of the current region; a calculation unit acquires the elastic energy density of the next region of the rock mass in the case of a rockburst in the current region and calculates the rockburst parameters of the next region at least based on the elastic energy density of the next region; a second determination unit determines whether a rockburst has occurred in the next region based on the rockburst parameters of the next region; and a repetition unit repeats the calculation step and the second determination step in sequence at least once in the case of a rockburst in the next region until no rockburst occurs in the next region. By calculating the rockburst parameters of the current region and the next region, the precise rockburst phenomenon can be quantitatively determined, especially in unconstructed areas, so that the risk of rockburst can be predicted in advance, solving the problem that the prior art cannot accurately and quantitatively predict rockburst.

[0146] The above description is merely a preferred embodiment of the present application and is not intended to limit the present application. Various modifications and variations are possible for those skilled in the art. Any modifications, equivalent substitutions, or improvements made within the spirit and principles of the present application shall be included within the scope of protection of the present application.

Claims

1. A rockburst determination method, characterized in that: include: an obtaining step of obtaining an elastic energy density of a current region of the rock mass, and calculating a rockburst parameter of the current region based at least on the elastic energy density of the current region, wherein the rock mass is a rock mass for underground engineering and includes multiple regions, the current region is any one of the multiple regions, and the elastic energy density is the elastic energy of the rock mass per unit volume; A first determining step is to determine whether a rock burst occurs in the current area according to rock burst parameters of the current area; a calculating step of obtaining, when a rockburst occurs in the current area, an elastic energy density of a next area of ​​the rock mass, and calculating a rockburst parameter of the next area based at least on the elastic energy density of the next area, wherein the next area is an adjacent area of ​​the current area; A second determining step is to determine whether a rock burst occurs in the next area according to the rock burst parameters of the next area; In the event that rock burst occurs in the next region, the calculating step and the second determining step are sequentially repeated at least once until rock burst does not occur in the next region, wherein during the repetition, the next region of the current time is an adjacent region of the previous next region, In the case of rock burst in the current area, the elastic energy density of the next area of ​​the rock mass is obtained, including: according to the seismic moment M0 of the rock burst in the current area, the magnitude M of the rock burst in the current area, the amplitude-energy relationship coefficients α1 and α2, the stiffness μ of the fault medium where the rock burst occurs in the current area, and the shear strength τ of the rock mass. y , seismic radiation efficiency η and Calculate the elastic energy ε of the next area of ​​the rock mass T ; According to the elastic energy of the next area of ​​the rock mass, the elastic energy of the next area per unit volume is calculated to obtain the elastic energy density of the next area of ​​the rock mass, The method further includes: obtaining engineering parameters of the underground project, wherein the engineering parameters include at least excavation method parameters and hard rock mechanical behavior parameters, the excavation method parameters are used to characterize the excavation method of the underground project, and the hard rock mechanical behavior parameters are used to characterize the hard rock mechanical behavior of the underground project; establishing an underground engineering model based on the engineering parameters; and when it is determined that a rockburst occurs in a corresponding area, obtaining, based on the underground engineering model, a time when the rockburst occurs in the corresponding area, a location where the rockburst occurs in the corresponding area, and a degree of damage caused by the rockburst in the corresponding area, and displaying the information on a display interface, wherein the corresponding area is the current area or the next area; Obtaining the excavation method parameters includes: obtaining an excavation method for the underground project, a density of explosives, an explosion velocity of the explosives, an isentropic index of the explosives, a diameter of an explosion hole, a distance between centers of two adjacent explosion holes, a depth of the explosion hole, a propagation velocity of a crack generated by the explosion, a propagation velocity of a rarefaction wave in the explosion, and an average velocity of reflection of the rarefaction wave, wherein the rarefaction wave is generated by gas from the top to the bottom of the explosion hole after the explosion; Determine the excavation method of the underground project. If the excavation method of the underground project is the drilling and blasting method, according to the formula Calculate the equivalent explosion stress peak P, where ρ e is the density of the explosive, V d is the explosion velocity of the explosive, γ is the isentropic index of the explosive; when the excavation method of the underground project is mechanical method, according to the formula Calculate the equivalent explosion stress peak value P, where r b is the diameter of the explosion hole, S is the distance between the centers of two adjacent explosion holes; according to the formula Calculate the load rise time t r , where L is the depth of the explosion hole; according to the formula Calculate the load duration t d , where c f is the expansion speed of the cracks generated by the explosion, c u1 is the propagation speed of the rarefaction wave in the explosion, c u2 is the average velocity of the rarefaction wave reflection; according to the formula Calculate the dynamic unloading duration t du , Obtaining the hard rock mechanical behavior parameters includes: according to the formula E(κ, I1) = B1 + B2κ + B3κ 2 +B4ln(I1), calculate the elastic-plastic coupling eigenvalue E, where κ is the plastic internal variable, I1 is the first stress invariant, and B1, B2, B3, and B4 are the first fitting parameters. When the plastic internal variable is less than or equal to the first threshold, according to the formula: Calculating strain hardening eigenvalues Among them, χ1, ζ1 and ξ1 are the first calculation parameters; when the plastic internal variable is greater than the first threshold, according to the formula Calculate the strain hardening eigenvalue in, is the residual internal friction angle; when the plastic internal variable is less than or equal to the second threshold, according to the formula: c(κ)=χ2κ 2 +ξ2κ+ξ2, calculate the strain softening characteristic value c, where χ2, ζ2 and ξ2 are the second calculation parameters; when the plastic internal variable is greater than the second threshold and less than or equal to the third threshold, according to the formula: c(κ)=χ3κ 2 +ξ3κ+ξ3, calculate the strain softening characteristic value c, where the third threshold is greater than the second threshold, χ3, ξ3 and ξ3 are third calculation parameters; when the plastic internal variable is greater than the third threshold, according to the formula c(κ)=c r , calculate the strain softening characteristic value c, where c r is the residual cohesion; according to the formula The dilution characteristic value ψ is calculated, where C1, C2, C3, C4, C5, and C6 are second fitting parameters.

2. The method according to claim 1, characterized in that Calculating rockburst parameters of the corresponding area based at least on the elastic energy density of the corresponding area, wherein the corresponding area is the current area or the next area, comprising: Obtaining the constraint pressure of the corresponding area; According to the constraint pressure p and the formula Calculate the rockburst parameter URLERI of the corresponding area, where U i is the elastic energy density of the rock mass at the i-th calculation step, U i+1 is the elastic energy density of the rock mass in the (i+1)th calculation step, and dt is the time step for calculating the rockburst parameters of the corresponding area.

3. The method according to claim 1, characterized in that Determining whether a rock burst occurs in a corresponding area according to rock burst parameters of the corresponding area, wherein the corresponding area is the current area or the next area, includes: When the rockburst parameter of the corresponding area is greater than 1, it is determined that a rockburst occurs in the corresponding area.

4. A rockburst determination device according to the rockburst determination method according to any one of claims 1 to 3, characterized in that: include: an acquisition unit, configured to acquire, in an acquisition step, an elastic energy density of a current region of the rock mass, and calculate a rockburst parameter of the current region based at least on the elastic energy density of the current region, wherein the rock mass is a rock mass for underground engineering and includes multiple regions, the current region is any one of the multiple regions, and the elastic energy density is the elastic energy of the rock mass per unit volume; A first determining unit is configured to determine, in a first determining step, whether a rockburst occurs in the current area according to rockburst parameters of the current area; a calculation unit configured to, in a calculation step, obtain the elastic energy density of a next region of the rock mass when a rock burst occurs in the current region, and calculate the rock burst parameters of the next region based at least on the elastic energy density of the next region, wherein the next region is an adjacent region to the current region; a second determining unit, configured to determine, in a second determining step, whether a rockburst occurs in the next area according to the rockburst parameters of the next area; A repeating unit is used to repeat the calculation step and the second determination step in sequence at least once when a rock burst occurs in the next area, until no rock burst occurs in the next area, wherein, during the repetition process, the current next area is an adjacent area of ​​the previous next area.

5. A computer-readable storage medium, characterized in that The computer-readable storage medium includes a stored program, wherein when the program is executed, the device where the computer-readable storage medium is located is controlled to execute the method according to any one of claims 1 to 3.

6. An electronic device comprising a memory and a processor, characterized in that: A computer program is stored in the memory, and the processor is configured to execute the method according to any one of claims 1 to 3 through the computer program.

Images