Rock mass damage degree classification analysis method, device, equipment and storage medium

Abstract

The present invention discloses a rock mass damage degree classification analysis method, device, equipment and storage medium, and relates to the field of data processing. The method comprises: establishing a three-dimensional geometric model of the disturbance area according to the construction disturbance center area, and determining the source information corresponding to multiple microseismic positioning event points according to the three-dimensional geometric model of the disturbance area; dividing the three-dimensional geometric model of the disturbance area into grids, and assigning grid disturbance range level, event residual level, grid residual level, and grid space density level according to the source information, and obtaining a rock mass damage degree determination index according to the assignment result; based on the principle of the four-section method, determining the rock mass damage degree classification threshold according to the rock mass damage degree determination index, and classifying the rock mass damage degree according to the rock mass damage degree classification threshold. In this way, using the geometric model and the index threshold classification based on multi-source source data, accurate body damage degree classification analysis results can be effectively obtained.

Description

Technical Field

[0001] The present invention relates to the field of data processing, and in particular to a rock mass damage degree classification analysis method, device, equipment and storage medium. Background Art

[0002] At present, the analysis of rock mass damage results based on microseismic monitoring mostly stays at the statistical parameter level, focusing on the temporal evolution characteristics of microseismic parameters to predict the precursors of surrounding rock rupture. Although this method can provide early warning information to a certain extent, it is insufficient in interpreting microseismic location results, and mostly relies on the qualitative analysis of single information, lacks comprehensive consideration of multi-source microseismic source information, and thus it is difficult to establish a quantitative relationship between microseismic source results and rock mass damage degree. Summary of the invention

[0003] In view of this, the purpose of the present invention is to overcome the shortcomings of the prior art and provide a rock damage degree classification analysis method, device, equipment and storage medium for integrating multi-source microseismic data, combining multi-dimensional parameters such as spatial density, energy, positioning accuracy of microseismic events, which can more comprehensively reflect the stress changes and damage characteristics of the rock mass during the construction process, thereby accurately evaluating the degree of rock damage.

[0004] The present invention provides the following technical solutions:

[0005] In the first aspect, the present invention proposes a method for grading and analyzing the degree of rock damage, comprising: establishing a three-dimensional geometric model of the disturbance area according to the construction disturbance center area, and determining the source information corresponding to multiple microseismic positioning event points according to the three-dimensional geometric model of the disturbance area; the source information includes the three-dimensional coordinates of the source event, the residual of the source event and the energy of the source event; dividing the three-dimensional geometric model of the disturbance area into grids to obtain multiple geometric model grids; determining the total number of event-containing grids according to the three-dimensional coordinates of the source event and each geometric model grid; calculating the distance between the grid center of each geometric model grid and the disturbance center of the three-dimensional geometric model of the disturbance area; assigning a grid disturbance range level to each geometric model grid according to a preset disturbance range classification threshold and each distance, and obtaining a grid disturbance range index according to the disturbance range assignment result; and The source event residual is assigned an event residual grade and a grid residual grade to obtain a grid source residual index; the event energy grade and the grid energy grade are assigned according to the source event energy to obtain a grid source energy index; the grid space density grade is assigned according to the number of grid events in each event-containing grid to obtain a grid source density index; the grid source energy-density comprehensive index is obtained according to the grid source energy index and the grid source density index; the rock mass damage degree determination index is obtained according to the preset index weight, the grid disturbance range index, the grid source residual index and the grid source energy-density comprehensive index; based on the four-point method principle, the rock mass damage degree classification threshold is determined according to the rock mass damage degree determination index, and the rock mass damage degree grade is assigned according to the rock mass damage degree classification threshold.

[0006] In the second aspect, the present invention proposes a rock damage degree classification analysis device, including: a construction module, used to establish a three-dimensional geometric model of the disturbance area according to the construction disturbance center area, and determine the source information corresponding to multiple microseismic positioning event points according to the three-dimensional geometric model of the disturbance area; the source information includes the three-dimensional coordinates of the source event, the residual of the source event and the energy of the source event; a division module, used to divide the three-dimensional geometric model of the disturbance area into grids to obtain multiple geometric model grids; a calculation module, used to determine the total number of event-containing grids according to the three-dimensional coordinates of the source event and each geometric model grid; a first assignment module, used to calculate the distance between the grid center of each geometric model grid and the disturbance center of the three-dimensional geometric model of the disturbance area; assign a grid disturbance range level to each geometric model grid according to a preset disturbance range classification threshold and each distance, and obtain a grid disturbance range index according to the disturbance range assignment result; a second assignment module, used to The first module is used to assign event residual levels and grid residual levels according to the residual of the earthquake source event to obtain the grid earthquake source residual index; the third assignment module is used to assign event energy levels and grid energy levels according to the energy of the earthquake source event to obtain the grid earthquake source energy index; the fourth assignment module is used to assign grid space density levels according to the number of grid events in each event-containing grid to obtain the grid earthquake source density index; the determination module is used to obtain the grid earthquake source energy-density comprehensive index according to the grid earthquake source energy index and the grid earthquake source density index; the rock mass damage degree judgment index is obtained according to the preset index weight, the grid disturbance range index, the grid earthquake source residual index and the grid earthquake source energy-density comprehensive index; the fifth assignment module is used to determine the rock mass damage degree grading threshold according to the rock mass damage degree judgment index based on the four-point method principle, and assign the rock mass damage degree level according to the rock mass damage degree grading threshold.

[0007] In a third aspect, the present invention proposes a computer device, including a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the rock damage degree classification analysis method as in the first aspect is implemented.

[0008] In a fourth aspect, the present invention provides a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the rock damage degree classification analysis method as in the first aspect.

[0009] The rock damage degree classification analysis method, device, equipment and storage medium disclosed in the present invention establish a geometric model and divide the grid based on the disturbance center, establish a mapping relationship between microseismic event points and three-dimensional rock masses, realize the transformation from point to body, and make the analysis of rock damage more refined; at the same time, fully consider the multi-source microseismic source information, and establish a quantitative evaluation index between the microseismic source and rock damage based on the physical meaning of the microseismic source information, effectively avoid the influence of subjective factors, make the evaluation result more objective and accurate, thereby helping to improve the reliability of rock damage prediction and provide powerful guidance for engineering practice. BRIEF DESCRIPTION OF THE DRAWINGS

[0010] In order to more clearly illustrate the technical solution of the present invention, the following is a brief introduction to the drawings required for use in the embodiments. It should be understood that the following drawings only illustrate certain embodiments of the present invention and should not be regarded as limiting the scope of protection of the present invention. In each of the drawings, similar components are numbered similarly.

[0011] Figure 1 A schematic diagram of a process of the rock mass damage degree classification analysis method proposed in this embodiment is shown;

[0012] Figure 2 A schematic diagram of three-dimensional geometric modeling proposed in this embodiment is shown;

[0013] Figure 2 (a) shows a schematic diagram of three-dimensional geometric modeling of a tunneling tunnel proposed in this embodiment;

[0014] Figure 2 (b) shows a schematic diagram of the spatial distribution results of microseismic event points proposed in this embodiment;

[0015] Figure 2 (c) shows a schematic diagram of the interception position of the three-dimensional geometric modeling proposed in this embodiment;

[0016] Figure 2 (d) shows a projection diagram of the three-dimensional geometric modeling proposed in this embodiment after the interception range is divided into grids;

[0017] Figure 2 (e) shows a schematic diagram of the assignment result proposed in this embodiment;

[0018] Figure 3 A schematic diagram of the grid disturbance range proposed in this embodiment is shown;

[0019] Figure 3 (a) shows a schematic diagram of grid disturbance range level assignment proposed in this embodiment;

[0020] Figure 3(b) shows a schematic diagram of the microseismic event point and grid disturbance range threshold surface proposed in this embodiment;

[0021] Figure 3 (c) shows a schematic diagram of the grid disturbance range level division proposed in this embodiment;

[0022] Figure 3 (d) shows a schematic diagram of the grid disturbance range index proposed in this embodiment;

[0023] Figure 4 A residual schematic diagram proposed in this embodiment is shown;

[0024] Figure 4 (a) shows a schematic diagram of the event residual determination result proposed in this embodiment;

[0025] Figure 4 (b) shows a schematic diagram of the event residual projection result proposed in this embodiment;

[0026] Figure 4 (c) shows a schematic diagram of the grid source residual index proposed in this embodiment;

[0027] Figure 5 A schematic diagram showing the energy proposed in this embodiment is shown;

[0028] Figure 5 (a) shows a schematic diagram of event energy determination results proposed in this embodiment;

[0029] Figure 5 (b) shows a schematic diagram of grid source energy index proposed in this embodiment;

[0030] Figure 6 A schematic diagram showing the number and density of grid events proposed in this embodiment is shown;

[0031] Figure 6 (a) shows a schematic diagram of the grid event quantity statistics result proposed in this embodiment;

[0032] Figure 6 (b) shows a schematic diagram of the statistical results of the number of hierarchical grid events proposed in this embodiment;

[0033] Figure 6 (c) shows a schematic diagram of the grid source density index proposed in this embodiment;

[0034] Figure 7 shows a schematic diagram of indicators proposed in this embodiment;

[0035] Figure 7 (a) shows a schematic diagram of the grid calculation result of the earthquake source energy-earthquake source density superposition index proposed in this embodiment;

[0036] Figure 7 (b) shows a schematic diagram of the calculation results of the earthquake source energy dominant-non-earthquake source density dominant index proposed in this embodiment;

[0037] Figure 7 (c) shows a schematic diagram of the calculation results of the non-seismic source energy-dominated-seismic source density-dominated index proposed in this embodiment;

[0038] Figure 7 (d) shows a schematic diagram of the calculation results of the non-seismic energy-dominated-non-seismic density-dominated index proposed in this embodiment;

[0039] Figure 7 (e) shows a schematic diagram of the calculation results of the grid source energy-density comprehensive index proposed in this embodiment;

[0040] Figure 8 A schematic diagram of the rock mass damage degree result proposed in this embodiment is shown;

[0041] Figure 8 (a) shows a first result schematic diagram of the rock mass damage degree proposed in this embodiment;

[0042] Figure 8 (b) shows a second result schematic diagram of the rock mass damage degree proposed in this embodiment;

[0043] Figure 8 (c) shows a third result schematic diagram of the rock mass damage degree proposed in this embodiment;

[0044] Fig. 9 A schematic diagram of the structure of the rock mass damage degree classification analysis device proposed in this embodiment is shown.

[0045] Explanation of the accompanying drawings: 1-microseismic source point; 2-excavation tunnel; 3-critical interface of disturbance range; 4-disturbance range slice; 5-geometric model slice; 900-rock damage degree classification analysis device; 901-construction module; 902-division module; 903-calculation module; 904-first assignment module; 905-second assignment module; 906-third assignment module; 907-fourth assignment module; 908-determination module; 909-fifth assignment module. DETAILED DESCRIPTION

[0046] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions 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. The following detailed description of the embodiments of the present invention provided in the drawings is not intended to limit the scope of the invention claimed for protection, but only represents the selected embodiments of the present invention. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without making creative work belong to the scope of protection of the present invention. Hereinafter, the terms "including", "having" and their cognates that can be used in various embodiments of the present invention are only intended to indicate specific features, numbers, steps, operations, elements, components or combinations of the foregoing items, and should not be understood as first excluding the existence of one or more other features, numbers, steps, operations, elements, components or combinations of the foregoing items or increasing the possibility of one or more features, numbers, steps, operations, elements, components or combinations of the foregoing items. In addition, the terms "first", "second", "third" and the like are only used to distinguish the description and cannot be understood as indicating or implying relative importance. Unless otherwise defined, all terms (including technical terms and scientific terms) used herein have the same meaning as those generally understood by ordinary technicians in the field to which the various embodiments of the present invention belong. The terms, such as those defined in generally used dictionaries, will be interpreted as having the same meaning as the contextual meaning in the relevant technical field and will not be interpreted as having an idealized meaning or an overly formal meaning unless clearly defined in various embodiments of the present invention.

[0047] Example 1

[0048] The present disclosure provides a method for analyzing the degree of rock mass damage. Figure 1 A rock mass damage degree classification analysis method includes steps S101 to S109, and each step is described in detail below.

[0049] Step S101, establish a three-dimensional geometric model of the disturbance area according to the construction disturbance center area, and determine the source information corresponding to multiple microseismic location event points according to the three-dimensional geometric model of the disturbance area; the source information includes the three-dimensional coordinates of the source event, the source event residual and the source event energy.

[0050] In this embodiment, a three-dimensional geometric model of the disturbance area is established according to the construction disturbance center area. For example, the three-dimensional geometric model of the disturbance area is established with the tunneling head as the disturbance center, and the three-dimensional geometric model of the disturbance area is established with the tunneling direction as the Y axis. The three-dimensional dimensions of the model are {1500 m, 1500 m, 500}, and the coordinates of the tunneling head in the geometric model are {750 m, 750 m, 250}. Figure 2 As shown in (a).

[0051] Furthermore, the source information corresponding to the multiple microseismic location event points is determined according to the three-dimensional geometric model of the disturbance area. The source information includes the three-dimensional coordinates of the source event, the source event residual and the source event energy.

[0052] Step S102, gridding the three-dimensional geometric model of the disturbance area to obtain a plurality of geometric model grids; determining the total number of grids containing events according to the three-dimensional coordinates of the source event and each geometric model grid.

[0053] In this embodiment, the three-dimensional geometric model of the disturbance area is meshed to obtain multiple geometric model meshes. At the same time, the total number of event-containing meshes is determined based on the three-dimensional coordinates of the source event of the microseismic location event point and each geometric model mesh. .

[0054] Based on the disturbance center, a geometric model is established and a grid is divided, and a mapping relationship from microseismic event points to three-dimensional rock mass is established, realizing the transformation from point to body, making the analysis of rock mass failure more refined.

[0055] Demonstratively, the three-dimensional grid accuracy is divided into 25 m grids, the total number of units in the three-dimensional geometric model of the disturbance area is 72,000, and the total number of event points in the three-dimensional model unit is 516.

[0056] like Figure 2 As shown in (b), the spatial distribution results of microseismic event points within the Z-axis range [225 m 275 m] are intercepted and projected on the XOY plane. Figure 2 As shown in the gray area in (c), a 2D plane is projected on the XOY plane in the intercepted area, and the projection plane is divided into 25 m grids along the X and Y axis directions. The projection results are shown in Figure 2 As shown in (d), there are 105 event points in the interception range, and the number of projection result grids in the interception area is 3600. All geometric model grids in the projection plane are determined one by one to determine whether there are event points. The grids containing event points are assigned a value of 1. The result of the assignment is the grid containing event index IG, as shown in Figure 2 As shown in (e), the total number of event grids in the projection plane is counted There are 77 of them.

[0057] It should be noted that when reducing the three-dimensional space to a two-dimensional projection plane, the same horizontal plane as the excavation disturbance is selected through the Z-axis direction. Generally speaking, this plane has the largest disturbance and is also the plane of most concern in the construction process. It is also convenient for visualizing the drawing results, and the intercepted area can be changed artificially.

[0058] Step S103, calculating the distance between the grid center of each geometric model grid and the disturbance center of the three-dimensional geometric model of the disturbance area; assigning a grid disturbance range level to each geometric model grid according to a preset disturbance range classification threshold and each distance, and obtaining a grid disturbance range index according to the disturbance range assignment result.

[0059] In this embodiment, the mesh center of each geometric model mesh is calculated. The disturbance center of the 3D geometric model of the disturbance area The distance between , k =1,2,3… and set the disturbance range classification threshold As the radius, the grid disturbance range level is assigned to the 3D modeling area, such as Figure 3 As shown in (a), the distribution results of microseismic event points and grid disturbance range threshold surface in three-dimensional space are as follows: Figure 3 As shown in (b), the disturbance range assignment result is summarized as the grid disturbance range index Among them, the threshold surface is the critical threshold boundary spherical interface of each distance classification.

[0060] Specifically, the disturbance range classification threshold is preset ,in , and then calculate the distance to the disturbance center for all geometric model grids , the calculation formula is: ; In the formula, For the The mesh center coordinates of the geometric model mesh, where , , are the coordinates of the disturbance center.

[0061] Further, the distance is determined in sequence and classification threshold Size, when When the disturbance range is determined to be level 1, the value of the mesh of this part of the geometric model is assigned to 4. When the disturbance range is determined to be level 2, the value of this part of the geometric model grid is assigned to be 3. When the disturbance range is determined to be level 3, the value of this part of the geometric model grid is assigned to 2. When the disturbance range is determined to be level 4, the value of this part of the geometric model grid is assigned to 1. When the disturbance range level is determined to be level 5, the value of this part of the geometric model grid is assigned to 0, and the results of all grid disturbance range assignments in the three-dimensional grid are summarized as the grid disturbance range index. .

[0062] Exemplarily, the preset disturbance range classification thresholds are 100 m, 200 m, 300 m, and 400 m. The disturbance center coordinates {750 m, 750 m, 250 m} are taken as the sphere center. The distance from each geometric model grid to the disturbance center coordinates is calculated, and the preset disturbance range classification threshold is used as the radius to divide the grid disturbance range level of the three-dimensional modeling area, such as Figure 3 As shown in (c).

[0063] When the distance ≤100 m, the value of the geometric model grid is 4. ≤ 200 m, the geometric model is assigned a value of 3. ≤300 m, the grid is assigned a value of 2. When the distance is 300 m< When ≤400 m, the geometric model grid is assigned a value of 1, and the statistical grid assignment result is the grid disturbance range index ID ,like Figure 3 As shown in (d).

[0064] Step S104, assigning event residual levels and grid residual levels according to the source event residuals to obtain grid source residual indicators.

[0065] In this embodiment, each source event residual is assigned an event residual grade according to the source event residual, and a grid residual grade is further assigned based on the event residual grade assignment, and the grid source residual index is obtained by combining the assignment results.

[0066] In a specific embodiment, step S104 includes: determining the first residual maximum value and the first residual minimum value from the residuals of each earthquake source event; based on the quartering principle, calculating the event residual grading threshold according to the first residual maximum value and the first residual minimum value; assigning event residual grades and grid residual grades according to the event residual grading threshold to obtain a grid earthquake source residual index.

[0067] In this embodiment, from all event residual values Determine the first residual maximum value in and the first residual minimum is the interval endpoint; according to the principle of quartering, according to the maximum value of the first residual and the first residual minimum Calculate event residual classification threshold ,in , the calculation formula of event residual classification threshold is: , , .

[0068] Furthermore, event residual grades and grid residual grades are assigned according to the event residual grading threshold, and the grid source residual index is comprehensively obtained.

[0069] In a specific embodiment, an event residual grade is assigned to each earthquake source event residual according to an event residual grade threshold, and the number of events corresponding to each event residual grade is obtained according to the event residual grade assignment result; the grid residual of each event-containing grid is obtained according to the preset grid event residual weight and the number of events of each event residual grade; the second residual maximum value and the second residual minimum value are determined from each grid residual, and the grid residual grade threshold is calculated based on the second residual maximum value and the second residual minimum value based on the quartering principle; grid residual grade is assigned to each geometric model grid according to the grid residual grade threshold and the grid residual, and the grid earthquake source residual index is obtained according to the grid residual grade assignment result.

[0070] In this embodiment, each event residual is determined based on the event residual classification. Since the larger the source event residual is, the lower the positioning accuracy is, the When the residual level of the event is determined to be 4, the residual value of the source event is assigned to 1. When the residual level of the source event is determined to be level 3, the residual value of the source event is assigned to 2. When , the residual level of the source event is determined to be 2, and the residual value of the source event is assigned to 3. When , the residual level of the source event is determined to be 1, and the residual value of the source event is assigned to 4.

[0071] Count the number of events corresponding to the residual level of each source event in all geometric model grids , according to the residual level of the source event from level 1 to level 4, multiple preset event residual weights are determined , and then calculate the grid residuals of each event grid , the calculation formula is: ; In the formula, For the The event residual weights corresponding to the event grids, For the The number of events corresponding to the jth level of source event residuals in the event-containing grid, i = 1, 2, 3… v , v is the total number of event grids ,and =1, 2, 3, 4.

[0072] Further, the second residual maximum value is determined from each grid residual and the second residual minimum , according to the principle of quartering, according to the second residual maximum value and the second residual minimum Confirm the grid residual classification threshold ,in , the classification threshold calculation formula is: , , .

[0073] Furthermore, the grid residual classification is determined for each geometric modeling grid. When the grid residual level is determined to be 1, the grid value of the geometric model is assigned to 4. When the grid residual level is determined to be 2, the value of the geometric model grid is assigned to be 3. When the grid residual level is determined to be 3, the grid value of the geometric model is assigned to 2. When the grid residual level is determined to be 4, the geometric model grid is assigned a value of 1, and the grid residual level assignment results of all the geometric model grids are summarized as the grid source residual index. .

[0074] Exemplarily, the first residual maximum and the first residual minimum 6e respectively -5 and 1e -6 The threshold value of event residual classification calculated by the four-point method is 4.525 e -5 , 3.05e -5 , 1.575 e -5 , in order to determine the residual level of the source event for all events in turn.

[0075] when When , the residual level of the source event is determined to be 1, and the residual value of the source event is assigned to 4. When , the residual level of the source event is determined to be 2, and the residual value of the source event is assigned to 3. When the residual level of the source event is determined to be level 3, the residual value of the source event is assigned to 2. When the residual level of the earthquake source event is determined to be 4, the residual of the earthquake source event is assigned a value of 1. The residual determination result of the event is as follows: Figure 4 As shown in (a).

[0076] The intercept plane projection analysis is performed on the three-dimensional space. The intercept area is the Z axis range [225 m 275 m]. The intercept plane is projected on the XOY plane and a grid is drawn on the projection surface. The event residual projection results are as follows: Figure 4As shown in (b), the residual levels of all source events in all geometric model grids are accumulated. After accumulation, the maximum residual value of all grids (the second maximum residual value) is 9, and the minimum residual value of the grids (the second minimum residual value) is 1; the grid residual classification thresholds are determined to be 7, 5, and 3 according to the four-point method.

[0077] According to the grid residual classification threshold, all geometric model grids in the intercepted plane projection surface are judged for grid residual classification. When the grid residual level is determined to be 1, the grid value of the geometric model is assigned to 4. When the grid residual level is determined to be 2, the value of the geometric model grid is assigned to be 3. When the grid residual level is determined to be 3, the grid value of the geometric model is assigned to 2. When the grid residual level is determined to be 4, the geometric model grid is assigned a value of 1, and the result of the residual level assignment of all grids in the intercepted area is the grid source residual index. ,like Figure 4 As shown in (c).

[0078] Step S105, assigning event energy levels and grid energy levels according to the source event energy to obtain a grid source energy index.

[0079] In this embodiment, event energy levels are assigned according to the energy of the earthquake source event, and grid energy levels are assigned based on the event energy level assignments, thereby obtaining a grid earthquake source energy index.

[0080] In a specific embodiment, step S105 includes: determining a first energy maximum value and a first energy minimum value from the energy of each seismic source event; based on the quartering principle, calculating an event energy grading threshold value according to the first energy maximum value and the first energy minimum value; assigning event energy levels and grid energy levels according to the event energy grading threshold value to obtain a grid seismic source energy index.

[0081] In this embodiment, the energy of each source event Determine the first energy maximum and the first energy minimum is the interval endpoint; according to the principle of quartering, according to the first energy maximum and the first energy minimum Determine event energy classification threshold , and use the event energy classification threshold Event energy levels and grid energy levels are assigned to obtain grid source energy indicators.

[0082] Among them, , the classification threshold calculation formula is: , , .

[0083] In a specific embodiment, an event energy level is assigned to each earthquake source event energy according to an event energy grading threshold, and the number of events corresponding to each event energy level is obtained according to the event energy level assignment result; the grid event energy of each event-containing grid is obtained according to the preset grid event energy weight and the number of events corresponding to each event energy level; the second energy maximum value is determined from each grid event energy, and the grid energy grading threshold is calculated based on the quartering principle according to the second energy maximum value; the grid energy level is assigned to each geometric model grid according to the grid energy grading threshold and each grid event energy, and the grid earthquake source energy index is obtained according to the grid energy level assignment result.

[0084] In this embodiment, according to the event energy classification threshold The energy of each source event is determined by the event energy classification. When the event energy level is determined to be 1, the energy of the source event is assigned a value of 4. When the event energy level is determined to be 2, the energy of the source event is assigned a value of 3. When the event energy level is determined to be 3, the energy of the source event is assigned a value of 2. When , the event energy level is determined to be 4, and the energy of the source event is assigned a value of 1.

[0085] Count the number of events of each event energy level in all geometric model grids , according to the event energy level from level 1 to level 4, the preset event energy weight is determined , calculate the energy of each grid event , the calculation formula is: ; In the formula, For the The event energy weight of a grid containing events, For the The number of events corresponding to the jth event energy level in the event grid.

[0086] Further, the second energy maximum is determined from the energy of each grid event and the second energy minimum ,by As the grid energy classification threshold Determine the grid event level. When the grid event energy , assign a value of 1 to the geometric model grid, otherwise assign a value of 0, and summarize the energy level assignment results of all grids as the grid source energy index .

[0087] Demonstratively, the energy numerical sequences of all earthquake source events in the intercepted area are statistically analyzed, and the maximum value of the event energy sequence (the first energy maximum value) is 26.61 kJ, the minimum value of the event residual sequence (the first energy minimum value) is 3.59 kJ, and the event energy classification thresholds calculated by the four-point method are 20.855 kJ, 15.1 kJ, and 9.345 kJ.

[0088] According to the event energy classification threshold, the event energy level is judged in turn. When the grid energy level is determined to be 1, the grid is assigned a value of 4. When , the grid energy level is determined to be 2, and the grid is assigned a value of 3. When the grid energy level is determined to be 3, the grid is assigned a value of 2. When the grid energy level is determined to be 4, the grid is assigned a value of 1. The event energy determination result is as follows: Figure 5 As shown in (a).

[0089] For all event grids, the energy classification results of all events in each event grid are accumulated, the maximum grid event energy (second energy maximum) is 7, and the minimum grid event energy (second energy minimum) is 1; and the grids with event energy greater than 7 in all geometric model grids are counted and the grid energy level is assigned to 1, otherwise it is assigned to 0; the energy level results of all grids in the intercepted area are counted as the grid source energy index ,like Figure 5 As shown in (b).

[0090] Step S106, assigning grid space density levels according to the number of grid events in each event-containing grid to obtain a grid source density index.

[0091] In this embodiment, the number of grid events in each event-containing grid is obtained to assign a grid space density level, thereby obtaining a grid source density index.

[0092] In a specific embodiment, step S106 includes: obtaining the number of grid events in each event-containing grid; grading the number of grid events based on a preset arithmetic interval to obtain multiple grid event number levels; and determining the sum of the number of grids corresponding to each grid event number level; determining a grid event number grading threshold based on the sum of the number of event-containing grids and the sum of the number of grids corresponding to each grid event number level; assigning a grid space density level to each geometric model grid based on the grid event number grading threshold, and obtaining a grid source density index based on the grid space density level assignment result.

[0093] In this embodiment, get The number of grid events in the event-containing grid ; Then classify the number of all grid events with the preset arithmetic interval as 1, and count the total number of grids under the number of grid events of the jth classification , j =1, 2, 3…m, where m is the maximum number of events in a single grid.

[0094] The grid event quantity threshold calculation formula is: ; In the formula, is the number of grid events. u = 1, judge whether the above formula is satisfied, if not, , if the output is satisfied , and let the grid event number threshold . Determine the grid space density based on the number of grid events; if the threshold is exceeded The grid with negative grid density is assigned a value of 1, otherwise it is assigned a value of 0; the results of assigning the spatial density levels of all grids are summarized as the grid source density index. .

[0095] For example, the number of grid events in all geometric model grids in the intercepted area is counted, and the statistical results are as follows: Figure 6 As shown in (a), the grid numbers are summarized as a numerical sequence of the number of grid events. The number of all grid events is classified with an arithmetic interval of 1. The number of grids under each classification of grid event numbers is counted. The number of grids with a grid event number of 4 is 1, the number of grids with a grid event number of 3 is 6, the number of grids with a grid event number of 2 is 13, and the number of grids with a grid event number of 1 is 57. The statistical results are shown in Figure 6 As shown in (b).

[0096] It can be seen that the number of grids with a grid event number of 1 accounts for more than 70% of all grids containing events. Therefore, the grid event number classification threshold is confirmed to be 1. The grid with a grid event number greater than 1 is assigned a value of 1, and the grid with a grid event number equal to 1 is assigned a value of 0. The results of the grid event number classification assignments in the intercepted area are summarized as the grid source density index. ,like Figure 6 As shown in (c).

[0097] Step S107, obtaining a grid seismic source energy-density comprehensive index according to the grid seismic source energy index and the grid seismic source density index.

[0098] In this embodiment, according to the source energy index and grid source density index , calculate the Source energy-source density stacking index on event grid , source energy dominant-non-source density dominant index , non-seismic energy-dominated-seismic density-dominated indicators and non-seismic energy-dominated-non-seismic density-dominated indicators , respectively assign weights to the above indicators and further calculate the grid source energy-density comprehensive index .

[0099] In a specific embodiment, step S107 includes: obtaining a source energy-source density superposition index according to a grid source energy index and a grid source density index; obtaining a source energy dominant-non-source density dominant index according to a grid source energy index and a source energy-source density superposition index; obtaining a non-source energy dominant-source density dominant index according to a grid source density index and a source energy-source density superposition index; obtaining a non-source energy dominant-non-source density dominant index according to the source energy dominant-non-source density dominant index and the non-source energy dominant-source density dominant index; obtaining a grid source energy-density comprehensive index according to the source energy-source density superposition index, the source energy dominant-non-source density dominant index, the non-source energy dominant-source density dominant index and the non-source energy dominant-non-source density dominant index.

[0100] In this embodiment, A comprehensive index of grid source energy-density on a grid containing events The calculation formula is: , , , , .

[0101] In the formula, For the grid source energy index on the event grid, For the The grid source density index on the event grid is For the Mesh event metrics on an event-containing mesh.

[0102] Among them, grid calculation source energy-source density superposition index The calculation results are as follows Figure 7 As shown in (a), the source energy dominated - non-source density dominated index The calculation results are as follows Figure 7 As shown in (b), the non-seismic energy-dominated-seismic density-dominated index The calculation results are as follows Figure 7 As shown in (c), the non-seismic energy-dominated-non-seismic density-dominated index The calculation results are as follows Figure 7As shown in (d), the grid source energy-density comprehensive index The calculation results are as follows Figure 7 As shown in (e), the grid with a value of 4 represents that both energy and quantity are the main influencing factors in the grid, the grid with a value of 3 represents that event energy is the dominant influencing factor in the grid, the grid with a value of 2 represents that the number of events is the dominant influencing factor in the grid, and the grid with a value of 1 indicates that both the number of events and energy in the grid are relatively small.

[0103] Step S108, obtaining a rock mass damage degree determination index according to a preset index weight, a grid disturbance range index, a grid source residual index and a grid source energy-density comprehensive index.

[0104] In this embodiment, according to the grid disturbance range index , grid source residual index and grid source energy-density comprehensive index And the preset indicator weight A 1 , A 2 , A 3 , calculate the rock mass damage degree determination index , where the rock mass damage degree determination index of each event grid is The calculation formula is: .

[0105] Step S109, based on the principle of the four-part method, the rock mass damage degree classification threshold is determined according to the rock mass damage degree judgment index, and the rock mass damage degree grade is assigned according to the rock mass damage degree classification threshold.

[0106] In this embodiment, the rock mass damage degree index in each event grid is obtained. The maximum value of and minimum value As the endpoint of the four-point value interval, the rock mass damage degree classification threshold is confirmed according to the four-point principle ,in , the classification threshold calculation formula is: , , .

[0107] The rock mass damage degree is assigned according to the rock mass damage degree classification threshold. When , the rock mass damage level of the grid is determined to be level 1, and the grid containing the event is assigned a value of 4. When the grid rock mass damage level is determined to be 2, the value of the event grid is assigned to 3. When the grid rock mass damage level is determined to be level 3, the grid containing the event is assigned a value of 2. When , the grid rock mass damage level is determined to be level 4, the grid containing the event is assigned a value of 1, and the results of the assignment of the damage level of all grid rock masses are summarized.

[0108] Finally, based on the results of the rock damage degree classification assignment for all grids, rock damage degree level 1 is judged as large-scale damage, rock damage degree level 2 is judged as medium-scale damage, rock damage degree level 3 is judged as small-scale damage, rock damage degree level 4 is judged as micro-scale damage, and rock damage degree level 5 is judged as non-scale damage.

[0109] It should be noted that the weight is set as [0.6, 0.2, 0.2] with the disturbance range as the main influencing factor. Therefore, the rock mass damage degree evaluation index is The calculation formula is: .

[0110] The rock mass damage degree is calculated by calculating all the geometric model grids on the projection plane of the intercepted area. Figure 8 As shown in (a), when the disturbance range is the main influencing factor, the rock mass damage degree in the area closer to the disturbance center is higher, and the rock mass damage degree in the area farther from the disturbance center is lower. Setting the weight with the disturbance range as the main influencing factor can amplify the fracture behavior near the disturbance, and the evaluation result of the area directly disturbed by tunneling is more accurate.

[0111] Taking the magnitude residual as the main influencing factor, the index weight is set to [0.2, 0.6, 0.2]. Therefore, the rock mass damage degree evaluation index The calculation formula is: .

[0112] The rock mass damage degree is calculated for all grids on the projection plane of the intercepted area. Figure 8 As shown in (b), when the magnitude residual is the main influencing factor, it is found that the large-scale damage of the rock mass is still concentrated in the disturbance center, but there is medium-scale damage in the excavated area, indicating that the stability of the surrounding rock in this area has weakened and the potential risk of surrounding rock damage has increased. From the perspective of construction safety, attention should be paid to this area. High-precision positioning can more finely characterize rock mass fractures, with more accurate locations and higher credibility.

[0113] Taking the magnitude event energy-density as the main influencing factor, the index weight is set to [0.2, 0.2, 0.6]. Therefore, the rock mass damage degree evaluation index The calculation formula is: .

[0114] The rock mass damage degree is calculated for all grids on the projection plane of the intercepted area. Figure 8As shown in (c), when the energy-density of the magnitude event is the main influencing factor, compared with the evaluation result with the disturbance range as the main influencing factor, it can be clearly observed that there is medium-scale damage on the left side of the tunnel at the rear end of the excavation. Compared with the evaluation result with the magnitude residual as the main influencing factor, the large-scale damage of the rock mass near the disturbance center on the tunnel excavation face is significantly increased.

[0115] In the process of engineering application, the rock damage degree evaluation results of the three dynamic schemes can be integrated by changing the weights of the influencing factors to carry out a multi-dimensional evaluation of the rock damage caused by construction disturbance.

[0116] The rock damage degree classification analysis method proposed in this embodiment establishes a geometric model and divides the grid based on the disturbance center, establishes a mapping relationship between microseismic event points and three-dimensional rock masses, realizes the transformation from point to body, and makes the analysis of rock damage more refined; at the same time, it fully considers the multi-source microseismic source information, and establishes a quantitative evaluation index between the microseismic source and rock damage based on the physical meaning of the microseismic source information, effectively avoiding the influence of subjective factors, making the evaluation result more objective and accurate, thereby helping to improve the reliability of rock damage prediction and providing strong guidance for engineering practice.

[0117] Example 2

[0118] In addition, the present disclosure provides a rock mass damage degree classification analysis device 900, see Fig. 9, including: a construction module 901, used to establish a three-dimensional geometric model of the disturbance area according to the construction disturbance center area, and determine the source information corresponding to multiple microseismic positioning event points according to the three-dimensional geometric model of the disturbance area; the source information includes the three-dimensional coordinates of the source event, the residual of the source event and the energy of the source event; a division module 902, used to divide the three-dimensional geometric model of the disturbance area into grids to obtain multiple geometric model grids; a calculation module 903, used to determine the total number of event-containing grids according to the three-dimensional coordinates of the source event and each geometric model grid; a first assignment module 904, used to calculate the distance between the grid center of each geometric model grid and the disturbance center of the three-dimensional geometric model of the disturbance area; assign a grid disturbance range level to each geometric model grid according to a preset disturbance range classification threshold and each distance, and obtain a grid disturbance range index according to the disturbance range assignment result; a second assignment module 905, used to assign a grid disturbance range level according to the residual of the source event The event residual level is assigned and the grid residual level is assigned to obtain the grid source residual index; the third assignment module 906 is used to assign event energy level and grid energy level according to the energy of the source event to obtain the grid source energy index; the fourth assignment module 907 is used to assign grid space density level according to the number of grid events in each event-containing grid to obtain the grid source density index; the determination module 908 is used to obtain the grid source energy-density comprehensive index according to the grid source energy index and the grid source density index; the rock mass damage degree determination index is obtained according to the preset index weight, the grid disturbance range index, the grid source residual index and the grid source energy-density comprehensive index; the fifth assignment module 909 is used to determine the rock mass damage degree classification threshold according to the rock mass damage degree determination index based on the four-point method principle, and assign the rock mass damage degree level according to the rock mass damage degree classification threshold.

[0119] The device provided in the embodiment of the present disclosure can execute the steps of the rock damage degree classification analysis method provided in Example 1, which will not be described again to avoid repetition.

[0120] Example 3

[0121] In addition, the embodiment of the present disclosure provides a computer device, including a memory and a processor, the memory stores a computer program, and when the computer program is executed by the processor, the rock damage degree classification analysis method of embodiment 1 is implemented. The device provided by the embodiment of the present disclosure can execute the steps of the rock damage degree classification analysis method provided in embodiment 1, and to avoid repetition, they are not repeated.

[0122] Example 4

[0123] The embodiment of the present disclosure proposes a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, the rock damage degree classification analysis method of the present embodiment 1 is implemented. In the present embodiment, the computer-readable storage medium can be a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk, etc. The computer-readable storage medium provided in the present embodiment can implement the rock damage degree classification analysis method provided in the embodiment 1, and in order to avoid repetition, it will not be repeated here.

[0124] In all examples shown and described herein, any specific value should be interpreted as being merely exemplary and not limiting, and therefore, other examples of the exemplary embodiments may have different values. It should be noted that similar reference numerals and letters represent similar items in the following figures, and therefore, once an item is defined in one figure, it does not need to be further defined and explained in subsequent figures. The above embodiments only express several implementation methods of the present invention, and their descriptions are more specific and detailed, but they cannot be understood as limiting the scope of the present invention. It should be pointed out that, for those of ordinary skill in the art, without departing from the concept of the present invention, several variations and improvements can also be made, and these all belong to the protection scope of the present invention.

Claims

1. A method for analyzing the degree of rock mass damage by grading, characterized in that: include: Establishing a three-dimensional geometric model of the disturbance area according to the construction disturbance center area, and determining the source information corresponding to multiple microseismic location event points according to the three-dimensional geometric model of the disturbance area; The earthquake source information includes three-dimensional coordinates of earthquake source events, earthquake source event residuals and earthquake source event energy; Dividing the three-dimensional geometric model of the disturbance area into grids to obtain a plurality of geometric model grids; Determine the total number of grids containing events according to the three-dimensional coordinates of the earthquake source event and each of the geometric model grids; Calculating the distance between the grid center of each of the geometric model grids and the disturbance center of the three-dimensional geometric model of the disturbance area; assigning a grid disturbance range grade to each of the geometric model grids according to a preset disturbance range classification threshold and each of the distances, and obtaining a grid disturbance range index according to the disturbance range assignment result; According to the earthquake source event residual, an event residual grade and a grid residual grade are assigned to obtain a grid earthquake source residual index; According to the energy of the earthquake source event, an event energy level is assigned and a grid energy level is assigned to obtain a grid earthquake source energy index; Assigning grid space density levels according to the number of grid events in each of the event-containing grids to obtain a grid source density index; Obtaining a grid seismic source energy-density comprehensive index according to the grid seismic source energy index and the grid seismic source density index; Obtaining a rock mass damage degree determination index according to a preset index weight, the grid disturbance range index, the grid source residual index, and the grid source energy-density comprehensive index; Based on the principle of the four-part method, the rock mass damage degree classification threshold is determined according to the rock mass damage degree judgment index, and the rock mass damage degree is classified according to the rock mass damage degree classification threshold.

2. The rock mass damage degree classification analysis method according to claim 1 is characterized in that: The step of assigning event residual levels and grid residual levels according to the seismic source event residuals to obtain a grid seismic source residual index includes: determining a first residual maximum value and a first residual minimum value from each of the seismic source event residuals; Based on the four-part principle, calculating the event residual classification threshold according to the first residual maximum value and the first residual minimum value; The event residual grade and the grid residual grade are assigned according to the event residual grade threshold to obtain the grid source residual index.

3. The rock mass damage degree classification analysis method according to claim 2 is characterized in that: The step of assigning event residual levels and grid residual levels according to the event residual grading threshold to obtain the grid source residual index includes: Assigning an event residual grade to each of the earthquake source event residuals according to the event residual grade threshold, and obtaining the number of events corresponding to each event residual grade according to the event residual grade assignment result; Obtaining the grid residual of each event-containing grid according to the preset grid event residual weight and the number of events of each event residual level; Determine a second residual maximum value and a second residual minimum value from each of the grid residuals, and calculate a grid residual classification threshold according to the second residual maximum value and the second residual minimum value based on the four-section principle; Grid residual grading is assigned to each of the geometric model grids according to the grid residual grading threshold and the grid residual, and the grid source residual index is obtained according to the grid residual grade assignment result.

4. The rock mass damage degree classification analysis method according to claim 1 is characterized in that: The assigning of event energy levels and grid energy levels according to the seismic source event energy to obtain a grid seismic source energy index includes: Determining a first energy maximum and a first energy minimum from each of the source event energies; Based on the four-division principle, calculating an event energy classification threshold according to the first energy maximum value and the first energy minimum value; Event energy levels and grid energy levels are assigned according to the event energy classification threshold to obtain the grid source energy index.

5. The rock mass damage degree classification analysis method according to claim 4 is characterized in that: The assigning of event energy levels and grid energy levels according to the event energy classification threshold to obtain the grid source energy index includes: Assigning event energy levels to the energy of each of the seismic source events according to the event energy grading threshold, and obtaining the number of events corresponding to each event energy level according to the event energy level assignment result; Obtaining the grid event energy of each event-containing grid according to the preset grid event energy weight and the number of events corresponding to each event energy level; Determine a second energy maximum value from each of the grid event energies, and calculate a grid energy classification threshold according to the second energy maximum value based on the four-division principle; A grid energy grade is assigned to each of the geometric model grids according to the grid energy grading threshold and the energy of each grid event, and a grid source energy index is obtained according to the grid energy grade assignment result.

6. The rock mass damage degree classification analysis method according to claim 1 is characterized in that: The grid space density level is assigned according to the number of grid events in each of the event-containing grids to obtain a grid source density index, including: Obtaining the number of grid events in each of the event-containing grids; Based on a preset arithmetic progression, the number of grid events is graded to obtain a plurality of grid event number levels; and the sum of the number of grids corresponding to each grid event number level is determined; Determine a grid event quantity classification threshold according to the sum of the number of grids containing events and the sum of the number of grids corresponding to the level of each grid event quantity; A grid space density grade is assigned to each of the geometric model grids according to the grid event quantity classification threshold, and the grid source density index is obtained according to the grid space density grade assignment result.

7. The rock mass damage degree classification analysis method according to claim 1 is characterized in that: The grid seismic source energy-density comprehensive index is obtained according to the grid seismic source energy index and the grid seismic source density index, including: Obtaining a superposition index of seismic source energy and seismic source density according to the grid seismic source energy index and the grid seismic source density index; Obtaining a source energy dominant-non-source density dominant index according to the grid source energy index and the source energy-source density superposition index; Obtaining a non-seismic source energy-dominant-seismic source density-dominant index according to the grid seismic source density index and the seismic source energy-seismic source density superposition index; According to the earthquake source energy dominant-non-seismic source density dominant index and the non-seismic source energy dominant-seismic source density dominant index, a non-seismic source energy dominant-non-seismic source density dominant index is obtained; The grid seismic source energy-density comprehensive index is obtained according to the seismic source energy-seismic source density superposition index, the seismic source energy dominant-non-seismic source density dominant index, the non-seismic source energy dominant-seismic source density dominant index and the non-seismic source energy dominant-non-seismic source density dominant index.

8. A rock mass damage degree classification analysis device, characterized in that: include: A construction module is used to establish a three-dimensional geometric model of the disturbance area according to the construction disturbance center area, and determine the source information corresponding to multiple microseismic location event points according to the three-dimensional geometric model of the disturbance area; the source information includes the three-dimensional coordinates of the source event, the source event residual and the source event energy; A partitioning module, used for partitioning the three-dimensional geometric model of the disturbance area into grids to obtain a plurality of geometric model grids; A calculation module, used for determining the total number of grids containing events according to the three-dimensional coordinates of the earthquake source event and each of the geometric model grids; The first assignment module is used to calculate the distance between the grid center of each of the geometric model grids and the disturbance center of the three-dimensional geometric model of the disturbance area; assign a grid disturbance range level to each of the geometric model grids according to a preset disturbance range classification threshold and each of the distances, and obtain a grid disturbance range index according to the disturbance range assignment result; A second assignment module is used to assign event residual levels and grid residual levels according to the source event residuals to obtain a grid source residual index; A third assignment module is used to assign event energy levels and grid energy levels according to the energy of the earthquake source event to obtain a grid earthquake source energy index; A fourth assignment module is used to assign a grid space density level according to the number of grid events in each of the event-containing grids to obtain a grid source density index; A determination module, used to obtain a grid seismic source energy-density comprehensive index according to the grid seismic source energy index and the grid seismic source density index; Obtaining a rock mass damage degree determination index according to a preset index weight, the grid disturbance range index, the grid source residual index, and the grid source energy-density comprehensive index; The fifth assignment module is used to determine the rock mass damage degree classification threshold based on the four-part principle and the rock mass damage degree judgment index, and to assign the rock mass damage degree grade based on the rock mass damage degree classification threshold.

9. A computer device, characterized in that: The method comprises a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the method for grading and analyzing the degree of rock damage as claimed in any one of claims 1 to 7 is implemented.

10. A computer-readable storage medium, characterized in that: It stores a computer program, which, when executed by a processor, implements the rock mass damage degree classification analysis method as described in any one of claims 1 to 7.

Images