A method for evaluating damage degree of rock mass

Abstract

The application relates to the technical field of surrounding rock damage detection, and relates to a rock mass damage degree evaluation method.The method comprises the following steps: collecting microseismic data corresponding to each microseismic event through a microseismic monitoring system arranged on a rock mass to be monitored and recording microseismic moments.According to the microseismic data corresponding to the microseismic event occurring at each microseismic moment, the rock mass fracture form and the rock mass damage amount of the microseismic event occurring at the microseismic moment are determined.According to the rock mass fracture form of the microseismic event occurring at the microseismic moment and a preset fractal dimension algorithm, the microseismic time dimension value corresponding to the microseismic event occurring at the microseismic moment is determined.According to the microseismic time dimension value corresponding to the microseismic event occurring at each microseismic moment and the rock mass damage amount, a rock mass damage degree evaluation index is determined.The application can analyze and evaluate the whole process of gradual accumulation of rock mass damage and can determine the signs of underground disasters.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of surrounding rock damage detection, and particularly relates to a rock mass damage degree evaluation method. BACKGROUND

[0002] At present, rock burst and rock pressure disaster are the main reasons for restricting roadway excavation. Rock mass shows a deterioration trend under load, often accompanied by energy accumulation and crack propagation, thereby inducing underground disasters.

[0003] The existing rock mass damage degree evaluation method mainly evaluates from the changes of rock mass elastic modulus, dissipation energy and damage energy. These methods mostly need to measure the laboratory reference parameters, and cannot be accurately and quickly applied to the actual engineering site. At the same time, the stability analysis of surrounding rock damage on site mostly revolves around the monitored data, and cannot form a relatively unified damage discrimination mode, and cannot effectively analyze and evaluate the whole process of rock mass damage accumulation, thereby delaying the best opportunity for identifying the signs of underground disasters.

[0004] Therefore, there is an urgent need for a rock mass damage degree evaluation method. SUMMARY

[0005] Therefore, it is necessary to provide a rock mass damage degree evaluation method in view of the above technical problems.

[0006] In a first aspect, a rock mass damage degree evaluation method is provided, and the method comprises:

[0007] Through a microseismic monitoring system arranged on the rock mass to be monitored, microseismic data corresponding to each microseismic event is collected, and a microseismic moment is recorded;

[0008] For each microseismic moment, according to the microseismic data corresponding to the microseismic event occurring at the microseismic moment, the rock mass fracture form and the rock mass damage amount of the microseismic event occurring at the microseismic moment are determined;

[0009] According to the rock mass fracture form of the microseismic event occurring at the microseismic moment and a preset fractal dimension algorithm, a microseismic time dimension value corresponding to the microseismic event occurring at the microseismic moment is determined;

[0010] According to the microseismic time dimension value corresponding to the microseismic event occurring at each microseismic moment and the rock mass damage amount, a rock mass damage degree evaluation index is determined, the rock mass damage degree evaluation index is used for analyzing and evaluating the whole process of rock mass damage accumulation, and the signs of underground disasters are identified.

[0011] As an optional implementation, the microseismic data comprises P waves, and the rock mass fracture form of the microseismic event occurring at each microseismic moment is determined according to the microseismic data corresponding to the microseismic event occurring at the microseismic moment, and the rock mass fracture form of the microseismic event occurring at each microseismic moment comprises:

[0012] The P wave corresponding to the microseismic event occurring at each microseismic moment is interpreted to obtain the moment tensor of the microseismic event occurring at the microseismic moment;

[0013] The moment tensor of the microseismic event occurring at each microseismic moment is decomposed into an inflation center, a pure double force couple and a compensated linear vector doublet based on a moment tensor decomposition method;

[0014] The tensile failure value and the shear failure value of the microseismic event occurring at each microseismic moment are determined according to the inflation center, the pure double force couple and the compensated linear vector doublet corresponding to the moment tensor of the microseismic event occurring at the microseismic moment;

[0015] The rock mass fracture form of the microseismic event occurring at each microseismic moment is determined according to the tensile failure value, the shear failure value, a preset first failure threshold and a preset second failure threshold of the microseismic event occurring at the microseismic moment.

[0016] As an optional implementation, the microseismic data comprises elastic strain energy released by a microseismic event, and maximum principal stress, intermediate principal stress, minimum principal stress, Poisson's ratio, volume and damage elastic modulus of each rock mass unit within a damage radius, and the rock mass damage amount of the microseismic event occurring at each microseismic moment is determined according to the microseismic data corresponding to the microseismic event occurring at the microseismic moment, and the rock mass damage amount of the microseismic event occurring at each microseismic moment comprises:

[0017] For each rock mass unit in the rock mass to be monitored, the elastic strain energy released when the rock mass unit is damaged is determined according to the maximum principal stress, the intermediate principal stress, the minimum principal stress, Poisson's ratio, the volume and the damage elastic modulus of the rock mass unit;

[0018] For each microseismic moment, the rock mass damage amount of the microseismic event occurring at the microseismic moment is determined according to the elastic strain energy released by the microseismic event occurring at the microseismic moment and the total elastic strain energy released when each rock mass unit within the damage radius of the microseismic event is damaged.

[0019] As an optional implementation, the formula for determining the elastic strain energy released when each rock mass unit in the rock mass to be monitored is damaged is:

[0020] ;

[0021] wherein, represents the elastic strain energy released when the rock mass unit within the damage radius of the microseismic event is damaged.i elastic strain energy released by the i-th rock unit within the damage radius of the microseismic event, σ 1 is the maximum principal stress of the rock unit, σ 2 is the intermediate principal stress of the rock unit, σ 3 is the minimum principal stress of the rock unit, v is the Poisson's ratio; V i the volume of the rock unit, E Di is the damage elastic modulus of the rock unit.

[0022] As an optional implementation, the formula for determining the rock damage amount of the microseismic event occurring at each microseismic moment according to the elastic strain energy released by the microseismic event occurring at the microseismic moment and the total elastic strain energy within the damage radius of the microseismic event is:

[0023] ;

[0024] wherein, ∆U represents the elastic strain energy released by the microseismic event, represents the elastic strain energy released by the i-th rock unit within the damage radius of the microseismic event when the rock unit is damaged, and m is the number of rock units within the damage radius of the microseismic event, represents the total elastic strain energy released by the rock unit within the damage radius of the microseismic event when the rock unit is damaged.

[0025] As an optional implementation, the formula for determining the tensile failure value and the shear failure value of the microseismic event occurring at each microseismic moment according to the explosion center, the pure double force couple and the compensation linear vector doublet corresponding to the moment tensor of the microseismic event is:

[0026] ;

[0027] ;

[0028] wherein, S1 is the tensile failure value, S2 is the shear failure value, M ISO is the explosion center, M DC is the pure double force couple, and M CLVD is the compensation linear vector doublet.

[0029] As an optional implementation, the formula for determining the rock fracture form of the microseismic event occurring at each microseismic moment according to the tensile failure value, the shear failure value, the preset first failure threshold and the preset second failure threshold is:

[0030] The tensile failure value and the shear failure value of the microseismic event occurring at each microseismic moment, if the tensile failure value is less than the preset first failure threshold, the rock mass failure form of the microseismic event occurring at the microseismic moment is tensile failure form;

[0031] If the shear failure value is greater than the preset second failure threshold, the rock mass failure form of the microseismic event occurring at the microseismic moment is shear failure form.

[0032] As an optional implementation, the formula for determining the microseismic time dimension value corresponding to the microseismic event occurring at the microseismic moment according to the rock mass failure form of the microseismic event occurring at the microseismic moment and the preset fractal dimension algorithm is:

[0033] ;

[0034] ;

[0035] Wherein, d i The microseismic event dimension value is represented, A represents an arbitrary non-empty bounded subset of Rn space, N(Ai) represents the i-th microseismic moment, for any r>0, the minimum number of microseismic events in space distribution is measured by a positive cube of length r, and n is the dimension of the microseismic event distribution graph in space.

[0036] As an optional implementation, the rock mass damage degree evaluation index is determined according to the microseismic time dimension value corresponding to each microseismic event occurring at the microseismic moment and the rock mass damage amount, comprising:

[0037] For each microseismic moment, the sum of the rock mass damage amounts of all microseismic events occurring from the first microseismic moment to the microseismic moment is determined as the rock mass cumulative damage amount corresponding to the microseismic moment;

[0038] The microseismic time dimension value and the rock mass cumulative damage amount of the microseismic moment form a rock mass damage matrix;

[0039] For each microseismic moment, the rock mass damage degree evaluation index corresponding to the microseismic event occurring at the microseismic moment is determined according to the rock mass damage matrix corresponding to the microseismic moment.

[0040] As an optional implementation, the formula for determining the rock mass damage degree evaluation index corresponding to the microseismic event occurring at the microseismic moment according to the rock mass damage matrix corresponding to the microseismic moment for each microseismic moment is:

[0041] ;

[0042] ;

[0043] wherein, M x is a rock mass damage matrix, microseismic time x corresponds to a time vector, is a microseismic time x corresponds to a microseismic time dimension value of a shear failure form, is a microseismic time x corresponds to a microseismic time dimension value of a tensile failure form, D x is a microseismic time x corresponds to a rock mass cumulative damage amount, C is a microseismic time x corresponds to a rock mass damage degree evaluation index of a microseismic event.

[0044] In a second aspect, a computer device is provided, comprising a memory and a processor, wherein the memory has stored thereon a computer program capable of running on the processor, and the processor implements the method steps of the first aspect when executing the computer program.

[0045] In a third aspect, a computer readable storage medium is provided, having stored thereon a computer program, and the computer program is executed by a processor to implement the method steps of the first aspect.

[0046] The application provides a rock mass damage degree evaluation method, and the technical scheme provided by the embodiment of the application at least brings the following beneficial effects: through a microseismic monitoring system arranged on a rock mass to be monitored, microseismic data corresponding to each microseismic event is collected, and a microseismic moment is recorded. For each microseismic moment, according to the microseismic data corresponding to the microseismic event occurring at the microseismic moment, the rock mass fracture form of the microseismic event occurring at the microseismic moment and the rock mass damage amount are determined. According to the rock mass fracture form of the microseismic event occurring at the microseismic moment and a preset fractal dimension algorithm, the microseismic time dimension value corresponding to the microseismic event occurring at the microseismic moment is determined. According to the microseismic time dimension value corresponding to the microseismic event occurring at each microseismic moment and the rock mass damage amount, a rock mass damage degree evaluation index is determined, the rock mass damage degree evaluation index is used for analyzing and evaluating the whole process of rock mass damage progressive accumulation, and the occurrence signs of underground disasters are distinguished. The application is based on microseismic monitoring, the microseismic data of different types of microseismic events corresponding to the destruction of the rock mass to be monitored is analyzed, and the rock mass fracture form is determined. Then the microseismic time dimension value corresponding to the microseismic event occurring at the microseismic moment is determined, and according to the microseismic time dimension value corresponding to the microseismic event occurring at each microseismic moment and the rock mass damage amount, the rock mass damage degree evaluation index is determined. Thus, the damage degree of the rock mass can be comprehensively obtained. The error judgment caused by a single microseismic event index is avoided, the accuracy of the stability judgment of the roadway surrounding rock is ensured, the damage condition of the roadway surrounding rock is comprehensively and accurately evaluated, and the design and implementation work of the roadway supporting scheme are guided.

[0047] It should be understood that the foregoing general description and the following detailed description are only exemplary and explanatory and are not restrictive of the application. BRIEF DESCRIPTION OF DRAWINGS

[0048] In order to more clearly illustrate the technical solutions in the embodiments of the application or the prior art, hereinafter, a brief introduction will be given to the drawings needed to be used in the embodiments or the prior art description. Obviously, the drawings in the following description are only some embodiments of the application, and for those skilled in the art, other drawings can also be obtained without creative labor on the basis of these drawings.

[0049] Figure 1 A flowchart of a rock mass damage degree evaluation method provided by the embodiment of the application;

[0050] Figure 2 A flowchart of a method for determining a rock mass fracture form provided by the embodiment of the application;

[0051] Figure 3 A flowchart of another method for determining a rock mass fracture form provided by the embodiment of the application;

[0052] Figure 4A flowchart of a method for determining rock mass damage amount provided by an embodiment of the present application is shown in FIG. 1.

[0053] Figure 5 A flowchart of a method for determining rock mass damage degree evaluation index provided by an embodiment of the present application is shown in FIG. 2.

[0054] Figure 6 A flowchart of an example of a rock mass damage degree evaluation method provided by an embodiment of the present application is shown in FIG. 3.

[0055] Figure 7 A structural schematic diagram of a computer device provided by an embodiment of the present application is shown in FIG. 4. DETAILED DESCRIPTION

[0056] In order to make the purposes, technical solutions and advantages of the present application clearer, the present application is further described in detail below with reference to the drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and do not limit the present application.

[0057] A rock mass damage degree evaluation method provided by an embodiment of the present application will be described in detail below with reference to the specific embodiments, Figure 1 A flowchart of a rock mass damage degree evaluation method provided by an embodiment of the present application is shown in FIG. 2. Figure 1 As shown in FIG. 2, the specific steps are as follows.

[0058] Step 101, through a microseismic monitoring system arranged on a rock mass to be monitored, microseismic data corresponding to each microseismic event is collected, and a microseismic time is recorded.

[0059] In implementation, a high-precision microseismic monitoring system is prearranged in a roadway, and through a plurality of sensors of the microseismic monitoring system preembedded in a rock mass to be measured of surrounding rock, microseismic data corresponding to each microseismic event is collected in real time, and a microseismic time is recorded through a computer.

[0060] Step 102, for each microseismic time, according to microseismic data corresponding to a microseismic event occurring at the microseismic time, a rock mass fracture form and a rock mass damage amount of the microseismic event occurring at the microseismic time are determined.

[0061] In implementation, for each microseismic time, according to microseismic data corresponding to a microseismic event occurring at the microseismic time, a rock mass fracture form and a rock mass damage amount of the microseismic event occurring at the microseismic time are determined. Specifically, a method of inversion calculation of a focal mechanism of the collected microseismic event can be used to determine the rock mass fracture form of the microseismic event. The rock mass damage amount corresponding to the microseismic event can be determined by using a fractal dimension calculation method in combination with a rock mass damage mechanism and a damage model.

[0062] Optionally, Figure 2A flowchart of a method for determining a rock mass rupture form according to an embodiment of the present application is shown in FIG. 1. As shown in FIG. 1, the microseismic data includes P waves. In step 102, for each microseismic moment, the rock mass rupture form of the microseismic event occurring at the microseismic moment is determined according to the microseismic data corresponding to the microseismic event occurring at the microseismic moment. The specific steps are as follows: Figure 2

[0063] In step 201, for each microseismic moment, the P waves corresponding to the microseismic event occurring at the microseismic moment are interpreted to obtain the moment tensor of the microseismic event occurring at the microseismic moment.

[0064] In implementation, each microseismic event can be regarded as a small earthquake event. The seismic waves generated during the microseismic event include S waves and P waves. If the earthquake is regarded as a point source, the moment tensor can be used to describe the deformation occurring at the source when the earthquake occurs. The seismic moment tensor is usually obtained by inverting the observed waveforms, so as to determine the source mechanism solution, and then analyze the tectonic stress field when the earthquake occurs. The inversion calculation of the microseismic event can be realized by interpreting the seismic waves. Therefore, in the embodiment of the present application, for each microseismic moment, the P waves corresponding to the microseismic event occurring at the microseismic moment are interpreted to obtain the moment tensor M of the microseismic event occurring at the microseismic moment.

[0065] In step 202, the moment tensor of the microseismic event occurring at the microseismic moment is decomposed into an inflation center, a pure double-couple, and a compensated linear vector doublet based on the moment tensor decomposition method.

[0066] In implementation, based on the moment tensor decomposition method, the moment tensor M of the microseismic event occurring at the microseismic moment can be decomposed into an inflation center M ISO , a pure double-couple M DC , and a compensated linear vector doublet M CLVD , that is, M=M ISO +M DC +M CLVD .

[0067] Alternatively, the specific steps of decomposing the moment tensor M into the inflation center M ISO , the pure double-couple M DC , and the compensated linear vector doublet M CLVD are as follows:

[0068] The moment tensor M is defined as M=M , and the moment tensor M is decomposed as follows: Therefore, the eigenvalues M1, M2, and M3 are respectively M1=(σ1+σ2+σ3) / 3, M2=(σ1-σ3) / 2, and M3=(2σ2-σ1-σ3) / 6.

[0069] wherein the three eigenvalues M1>M 2, ​> M3, M1 corresponding to the principal axis is T axis (tension axis), M2 corresponding to the principal axis is B axis (zero axis), M3 corresponding to the principal axis is P axis (compression axis).

[0070] Further, the expansion center M decomposed by the moment tensor M ISO , pure double force couple M DC And the compensation linear vector dipole M CLVD Can be expressed as follows:

[0071] ;

[0072] ;

[0073] .

[0074] Step 203, according to the expansion center, pure double force couple and compensation linear vector dipole corresponding to the moment tensor of the microseismic event occurred at the microseismic moment, determine the tensile failure value and shear failure value of the microseismic event occurred at the microseismic moment.

[0075] In the implementation, the tensile failure value and shear failure value of the microseismic event occurred at the microseismic moment can be determined according to the expansion center, pure double force couple and compensation linear vector dipole corresponding to the moment tensor of the microseismic event occurred at the microseismic moment.

[0076] Optionally, in step 203, for the moment tensor of each microseismic event occurred at the microseismic moment, the formula for determining the tensile failure value and shear failure value of the microseismic event occurred at the microseismic moment according to the expansion center, pure double force couple and compensation linear vector dipole corresponding to the moment tensor is:

[0077] ;

[0078] ;

[0079] Wherein, S1 is the tensile failure value, S2 is the shear failure value, M ISO Is the expansion center, M DC Is the pure double force couple, M CLVD Is the compensation linear vector dipole.

[0080] Step 204, according to the tensile failure value, shear failure value, first failure threshold and second failure threshold preset of the microseismic event occurred at the microseismic moment, determine the rock mass fracture form of the microseismic event occurred at the microseismic moment.

[0081] In the implementation, the rock mass fracture form of the microseismic event occurred at the microseismic moment can be determined according to the tensile failure value, shear failure value, first failure threshold and second failure threshold preset of the microseismic event occurred at the microseismic moment.

[0082] Optionally, Figure 3 A flowchart of another method for determining the rock mass failure mode provided by the embodiments of the present application is shown in Figure 3 As shown in FIG. 4, in step 204, for the tensile failure value and the shear failure value of the microseismic event occurring at each microseismic time, the specific steps of determining the rock mass failure mode of the microseismic event occurring at the microseismic time according to the tensile failure value, the shear failure value, the preset first failure threshold and the preset second failure threshold are as follows:

[0083] In step 301, for the tensile failure value and the shear failure value of the microseismic event occurring at each microseismic time, if the tensile failure value is less than the preset first failure threshold, the rock mass failure mode of the microseismic event occurring at the microseismic time is the tensile failure mode.

[0084] In implementation, different rock mass failure modes can correspond to different time instants during the microseismic event, and the rock mass failure modes include the tensile failure mode and the shear failure mode. For example, the first failure threshold can be 0.4, and when the tensile failure value is less than 0.4, the rock mass failure mode of the microseismic event is the tensile failure mode.

[0085] In step 302, if the shear failure value is greater than the preset second failure threshold, the rock mass failure mode of the microseismic event occurring at the microseismic time is the shear failure mode.

[0086] In implementation, if the shear failure value is greater than the preset second failure threshold, the rock mass failure mode of the microseismic event occurring at the microseismic time is the shear failure mode. For example, the second failure threshold can be 0.6, and when the tensile failure value is greater than 0.6, the rock mass failure mode of the microseismic event is the shear failure mode.

[0087] Optionally, Figure 4 A flowchart of a method for determining the rock mass damage provided by the embodiments of the present application is shown in Figure 4 As shown in FIG. 4, in step 102, the microseismic data includes the elastic strain energy released by the microseismic event, and the maximum principal stress, the intermediate principal stress, the minimum principal stress, the Poisson's ratio, the volume and the damage elastic modulus of each rock mass unit within the damage radius, for each microseismic time, according to the microseismic data corresponding to the microseismic event occurring at the microseismic time, the specific steps of determining the rock mass damage of the microseismic event occurring at the microseismic time are as follows:

[0088] In step 401, for each rock mass unit in the rock mass to be monitored, the elastic strain energy released by the rock mass unit when damaged is determined according to the maximum principal stress, the intermediate principal stress, the minimum principal stress, the Poisson's ratio, the volume and the damage elastic modulus of the rock mass unit.

[0089] In implementation, for each rock mass unit in the rock mass to be monitored, the elastic strain energy released when the rock mass unit is damaged can be determined according to the maximum principal stress, the intermediate principal stress, the minimum principal stress, the Poisson's ratio, the volume and the damage elastic modulus of the rock mass unit.

[0090] Optionally, in step 401, for each rock mass unit in the rock mass to be monitored, the formula for determining the elastic strain energy released when the rock mass unit is damaged is:

[0091] ;

[0092] wherein, represents the elastic strain energy released when the i th rock mass unit within the damage radius of the microseismic event is damaged, i σ 1 is the maximum principal stress of the rock mass unit, σ 2 is the intermediate principal stress of the rock mass unit, σ 3 is the minimum principal stress of the rock mass unit, v is the Poisson's ratio; V i is the volume of the rock mass unit, E Di is the damage elastic modulus of the rock mass unit.

[0093] In step 402, for each microseismic moment, the rock mass damage amount of the microseismic event occurring at the microseismic moment is determined according to the elastic strain energy released by the microseismic event occurring at the microseismic moment and the total elastic strain energy released when each rock mass unit within the damage radius of the microseismic event is damaged.

[0094] In implementation, the microseismic events of different failure forms can be analyzed to calculate the dimension number, to obtain the dimension number changes of shear, tension and mixed failure, and to determine the rock mass damage degree and the damage degree thereof in combination with the rock mass damage mechanism and the damage model. Therefore, for each microseismic moment, the rock mass damage amount of the microseismic event occurring at the microseismic moment can be determined according to the elastic strain energy released by the microseismic event occurring at the microseismic moment and the total elastic strain energy released when each rock mass unit within the damage radius of the microseismic event is damaged.

[0095] Optionally, in step 402, for each microseismic moment, the formula for determining the rock mass damage amount of the microseismic event occurring at the microseismic moment is:

[0096] ;

[0097] ​wherein, ΔU represents the elastic strain energy released by the microseismic event, represents the elastic strain energy released by the i-th rock mass unit within the damage radius of the microseismic event, and m is the number of rock mass units within the damage radius of the microseismic event, represents the total elastic strain energy released by the rock mass units within the damage radius of the microseismic event.

[0098] In step 103, according to the rock mass fracture form of the microseismic event occurring at the microseismic moment and the preset fractal dimension algorithm, the microseismic time dimension value corresponding to the microseismic event occurring at the microseismic moment is determined.

[0099] In implementation, after determining the rock mass fracture form of the microseismic event, fractal dimension calculation can be performed on the microseismic events of different fracture forms in time dimension. For the microseismic events in three-dimensional space, box dimension is used to perform fractal calculation on the range of the microseismic event at the moment, so as to determine the microseismic time dimension value corresponding to the microseismic event occurring at the microseismic moment.

[0100] Optionally, in step 103, according to the rock mass fracture form of the microseismic event occurring at the microseismic moment and the preset fractal dimension algorithm, the formula for determining the microseismic time dimension value corresponding to the microseismic event occurring at the microseismic moment is:

[0101] ;

[0102] ;

[0103] wherein, d i represents the microseismic event dimension value, A represents an arbitrary non-empty bounded subset of Rn space, N(Ai) represents the i-th microseismic moment, for any one r>0, the minimum number of microseismic events in space distribution is measured by using a positive cube with length r, and n is the dimension of the microseismic event distribution graph in space.

[0104] In implementation, the logr can be taken as the abscissa, the logN(Ai) can be taken as the ordinate, and the set A can be plotted in the logarithmic coordinate system, that is, the dimension variation image of the microseismic event. The dimension calculation of the microseismic event distribution graph in space takes time as the node, which can be divided into intervals of 8h. Along the time axis of the development of the microseismic event, eight hours of data are collected for summary calculation, so that the calculated dimension value represents the development degree of the microseismic time and space within 8h.

[0105] In step 104, according to the microseismic time dimension value corresponding to the microseismic event occurring at each microseismic moment and the rock mass damage amount, the rock mass damage degree evaluation index is determined, which is used to analyze and evaluate the whole process of the gradual accumulation of rock mass damage, and to judge the signs of underground disasters.

[0106] In practice, the sensitivity of the measurement matrix will be an intuitive reflection of the damage variables. Therefore, the sensitivity of the measurement matrix can be determined based on the microseismic time dimension and rock mass damage amount corresponding to the microseismic events occurring at each microseismic moment, and this sensitivity can be determined as an indicator for assessing the degree of rock mass damage.

[0107] Optionally, Figure 5 A flowchart illustrating a method for determining rock mass damage assessment indices, as provided in this application embodiment, is shown below. Figure 5 As shown, in step 104, the specific steps for determining the rock mass damage assessment index based on the microseismic time dimension and rock mass damage amount corresponding to the microseismic events occurring at each microseismic moment are as follows:

[0108] Step 501: For each microseismic moment, the sum of the rock mass damage amounts of all microseismic events that occurred from the first microseismic moment to that microseismic moment is determined as the cumulative rock mass damage amount corresponding to that microseismic moment.

[0109] In practice, for each microseismic moment, the sum of the rock mass damage amounts of all microseismic events that occurred from the first microseismic moment to that microseismic moment can be determined as the cumulative rock mass damage amount corresponding to that microseismic moment.

[0110] Step 502: Combine the microseismic time dimension at the microseismic moment with the cumulative rock mass damage to form a rock mass damage matrix.

[0111] In practice, the time dimension of the microseismic event at that moment and the cumulative damage to the rock mass can be combined to form a rock mass damage matrix. For example: Rock Mass Damage Matrix , Micro-earthquake moment x The corresponding microseismic time dimension of the shear failure mode, Micro-earthquake moment x The corresponding microseismic time dimension of the tension failure mode, D x Micro-earthquake moment x The corresponding cumulative damage to the rock mass.

[0112] Step 503: For each microseismic moment, determine the rock mass damage assessment index corresponding to the microseismic event that occurred at that microseismic moment based on the rock mass damage matrix corresponding to that microseismic moment.

[0113] In practice, for each microseismic moment, the rock mass damage assessment index corresponding to the microseismic event occurring at that moment can be determined based on the rock mass damage matrix corresponding to that microseismic moment.

[0114] Optionally, in step 503, for each microseismic moment, the formula for determining the rock mass damage degree evaluation index corresponding to the microseismic event occurring at the microseismic moment is determined according to the rock mass damage matrix corresponding to the microseismic moment.

[0115] ;

[0116] ;

[0117] wherein, M x is the rock mass damage matrix, is the microseismic moment x corresponding to the time vector, is the microseismic time dimension value of the microseismic event occurring at the microseismic moment x corresponding to the shear failure form, is the microseismic time dimension value of the microseismic event occurring at the microseismic moment x corresponding to the tensile failure form, D x is the rock mass cumulative damage amount of the microseismic moment x corresponding to the rock mass damage degree evaluation index of the microseismic event occurring at the microseismic moment C x .

[0118] Optionally, Figure 6 is the flow chart of an example of the rock mass damage degree evaluation method provided by the embodiments of the present application, as shown in Figure 6 , the specific steps are as follows:

[0119] Step 601, obtaining the microseismic data (such as coordinate, moment tensor, energy, frequency, source radius and other parameters) of the microseismic event;

[0120] Step 602, decomposing the moment tensor and further analyzing and determining the rock mass failure form;

[0121] Step 603, using the box dimension to calculate the dimension of the microseismic event of different rock mass failure forms, and obtaining the dimension change of the shear and tensile failure forms;

[0122] Step 604, calculating the sensitive value (rock mass damage degree evaluation index) of the rock mass damage matrix composed of the microseismic time dimension value and the rock mass cumulative damage amount, measuring the damage degree change according to the sensitive value of different moments, and then evaluating the damage.

[0123] ​The embodiment of the application provides a rock mass damage degree evaluation method, and the technical scheme provided by the embodiment of the application at least brings the following beneficial effects: through the microseismic monitoring system arranged on the rock mass to be monitored, microseismic data corresponding to each microseismic event is collected, and a microseismic moment is recorded. For each microseismic moment, according to the microseismic data corresponding to the microseismic event occurring at the microseismic moment, the rock mass fracture form of the microseismic event occurring at the microseismic moment and the rock mass damage amount are determined. According to the rock mass fracture form of the microseismic event occurring at the microseismic moment and the preset fractal dimension algorithm, the microseismic time dimension value corresponding to the microseismic event occurring at the microseismic moment is determined. According to the microseismic time dimension value corresponding to the microseismic event occurring at each microseismic moment and the rock mass damage amount, a rock mass damage degree evaluation index is determined, the rock mass damage degree evaluation index is used for analyzing and evaluating the whole process of rock mass damage progressive accumulation, and an underground disaster occurrence sign is judged. The application is based on microseismic monitoring, the microseismic data of different types of microseismic events corresponding to the destruction of the rock mass to be monitored is analyzed, the fracture form of the rock mass is determined. Then the microseismic time dimension value corresponding to the microseismic event occurring at the microseismic moment is determined, and the rock mass damage degree evaluation index is determined according to the microseismic time dimension value corresponding to the microseismic event occurring at each microseismic moment and the rock mass damage amount. Thus, the damage degree of the rock mass can be comprehensively obtained. The error judgment caused by a single microseismic event index is avoided, the accuracy of the stability judgment of the roadway surrounding rock is ensured, the damage condition of the roadway surrounding rock is comprehensively and accurately evaluated, and the design and implementation work of the roadway supporting scheme are guided.

[0124] It should be understood that although Figures 1 to 6 The steps in the flowchart of the method are shown in sequence according to the arrows, but these steps are not necessarily executed in the order indicated by the arrows. Unless otherwise specified herein, the execution of these steps is not strictly limited in sequence, and these steps can be executed in other orders. Moreover, Figures 1 to 6 At least part of the steps in the method can include multiple steps or multiple stages, which are not necessarily executed at the same time, but can be executed at different times, and the execution order of these steps or stages is not necessarily sequential, but can be executed alternately or alternately with at least part of other steps or steps or stages in other steps.

[0125] It can be understood that the same / similar parts of each embodiment of the method in the specification can be mutually referred to, and each embodiment mainly explains the difference from other embodiments, and the related part can be referred to the description of other method embodiments.

[0126] In one embodiment, a computer device is provided, such as Figure 7As shown, it comprises a memory and a processor, the memory has a computer program stored thereon, which can be run on the processor, and the processor implements the method steps of the rock mass damage degree evaluation method when the computer program is run.

[0127] In one embodiment, a computer readable storage medium has a computer program stored thereon, which, when executed by a processor, implements the steps of the rock mass damage degree evaluation method.

[0128] A person of ordinary skill in the art can understand that all or part of the processes in the above-mentioned embodiments can be completed by a computer program instructing related hardware, and the computer program can be stored in a non-volatile computer readable storage medium. When the computer program is executed, it can include the processes of the above-mentioned embodiments. Any reference to memory, storage, database or other medium used in the embodiments provided by the present application can include non-volatile and / or volatile memory. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM) or flash memory. Volatile memory can include random access memory (RAM) or external cache memory. As an illustration but not limitation, RAM is available in various forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDR SDRAM), enhanced SDRAM (ESDRAM), synchronous link (Synchlink) DRAM (SLDRAM), Rambus direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM) and the like.

[0129] It should be noted that, in this document, relational terms such as first and second and the like can only be used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply that there is any such actual relationship or order between these entities or operations. Moreover, the terms "comprise", "comprise" or any other variant thereof are intended to cover non-exclusive inclusion, so that a process, method, article or apparatus including a series of elements includes not only those elements, but also other elements not explicitly listed or inherent to such a process, method, article or apparatus. Without more limitations, the element defined by the statement "comprises a" does not exclude the presence of additional identical elements in the process, method, article or apparatus including the element.

[0130] It should be further noted that the user information (including but not limited to user equipment information, user personal information, etc.) and data (including but not limited to data for display, analyzed data, etc.) involved in the present application are all information and data authorized by the user or authorized by all parties.

[0131] Each of the embodiments in the specification is described in a related manner, and the same and similar parts between the embodiments can be referred to each other. Each embodiment focuses on the difference from other embodiments. In particular, for the system embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and the related parts can be referred to the part of the method embodiment.

[0132] The technical features of the above embodiments can be combined in any manner. In order to make the description simple, not all possible combinations of the technical features in the above embodiments are described, but as long as the combination of the technical features does not exist contradictory, it should be considered as the scope of the specification.

[0133] The above embodiments only express several implementation manners of the present application, and the description is relatively specific and detailed, but it should not be understood as a limitation on the scope of the patent. It should be pointed out that for those skilled in the art, without departing from the concept of the present application, some modifications and improvements can be made, which are all within the protection scope of the present application. Therefore, the protection scope of the patent of the present application should be subject to the appended claims.

Claims

1. A method for evaluating the damage degree of a rock mass, characterized by, The method comprises: Collecting microseismic data corresponding to each microseismic event through a microseismic monitoring system arranged on the rock mass to be monitored, and recording the microseismic time; For each microseismic time, determining the rock mass fracture form and rock mass damage amount of the microseismic event occurring at the microseismic time according to the microseismic data corresponding to the microseismic event occurring at the microseismic time; According to the rock mass fracture form of the microseismic event occurring at the microseismic time and the preset fractal dimension algorithm, the microseismic time dimension value corresponding to the microseismic event occurring at the microseismic time is determined; According to the microseismic time dimension value corresponding to the microseismic event occurring at each microseismic time and the rock mass damage amount, a rock mass damage degree evaluation index is determined, which is used to analyze and evaluate the whole process of rock mass damage progressive accumulation and judge the signs of underground disasters.

2. The method of claim 1, wherein, The microseismic data includes P wave, and the rock mass fracture form of the microseismic event occurring at each microseismic time is determined according to the microseismic data corresponding to the microseismic event occurring at the microseismic time, which comprises: For each microseismic time, the P wave corresponding to the microseismic event occurring at the microseismic time is interpreted to obtain the moment tensor of the microseismic event occurring at the microseismic time; Based on the moment tensor decomposition method, the moment tensor of the microseismic event occurring at the microseismic time is decomposed into an inflation center, a pure double force couple and a compensating linear vector doublet; According to the inflation center, the pure double force couple and the compensating linear vector doublet corresponding to the moment tensor of the microseismic event occurring at the microseismic time, the tensile failure value and the shear failure value of the microseismic event occurring at the microseismic time are determined; According to the tensile failure value, the shear failure value, the first damage threshold and the second damage threshold of the microseismic event occurring at the microseismic time, the rock mass fracture form of the microseismic event occurring at the microseismic time is determined.

3. The method of claim 1, wherein, The microseismic data includes the elastic strain energy released by the microseismic event, and the maximum principal stress, the intermediate principal stress, the minimum principal stress, the Poisson's ratio, the volume and the damage elastic modulus of each rock mass unit within the damage radius, and the rock mass damage amount of the microseismic event occurring at each microseismic time is determined according to the microseismic data corresponding to the microseismic event occurring at the microseismic time, which comprises: For each rock mass unit in the rock mass to be monitored, the elastic strain energy released when the rock mass unit is damaged is determined according to the maximum principal stress, the intermediate principal stress, the minimum principal stress, the Poisson's ratio, the volume and the damage elastic modulus of the rock mass unit; For each microseismic time, the rock mass damage amount of the microseismic event occurring at the microseismic time is determined according to the elastic strain energy released by the microseismic event and the total elastic strain energy released when each rock mass unit within the damage radius of the microseismic event is damaged.

4. The method of claim 3, wherein, The formula for determining the elastic strain energy released when each rock mass unit in the rock mass to be monitored is damaged according to the maximum principal stress, the intermediate principal stress, the minimum principal stress, the Poisson's ratio, the volume and the damage elastic modulus of the rock mass unit is: ； wherein, the elastic strain energy released by the damage of the i-th rock mass element within the damage radius of a microseismic event, i σ 1 is the maximum principal stress of the i-th rock mass element, σ 2 is the intermediate principal stress of the i-th rock mass element, σ 3 is the minimum principal stress of the i-th rock mass element, v is the Poisson's ratio; V i the volume of the i-th rock mass element, E Di is the damaged elastic modulus of the i-th rock mass element.​ 5. The method of claim 3, wherein, The formula for determining the rock mass damage quantity of the microseismic event occurring at each microseismic moment according to the elastic strain energy released by the microseismic event occurring at the microseismic moment and the total elastic strain energy within the damage radius of the microseismic event is: ； where ΔU represents the elastic strain energy released by the microseismic event, where ΔU represents the elastic strain energy released by the microseismic event, where ΔU represents the elastic strain energy released by the microseismic event, 6. The method of claim 2, wherein, The formula for determining the tensile failure value and the shear failure value of the microseismic event occurring at each microseismic moment according to the expansion center, the pure double force couple and the compensation linear vector doublet corresponding to the moment tensor of the microseismic event is: ； ； where S1 is the tensile failure value, S2 is the shear failure value, M ISO is the expansion center, M DC is the pure double force couple, M CLVD is the compensating linear vector dipole.

7. The method of claim 2, wherein, The formula for determining the rock mass failure form of the microseismic event occurring at each microseismic moment according to the tensile failure value, the shear failure value, the preset first failure threshold and the preset second failure threshold is: If the tensile failure value of the microseismic event occurring at each microseismic moment is less than the preset first failure threshold, the rock mass failure form of the microseismic event occurring at the microseismic moment is a tensile failure form; If the shear failure value of the microseismic event occurring at each microseismic moment is greater than the preset second failure threshold, the rock mass failure form of the microseismic event occurring at the microseismic moment is a shear failure form.

8. The method of claim 1, wherein, The formula for determining the microseismic time dimension value corresponding to the microseismic event occurring at each microseismic moment according to the rock mass failure form of the microseismic event occurring at the microseismic moment and the preset fractal dimension algorithm is: ； ； wherein d i wherein, represents the microseismic event dimension value, A represents an arbitrary non-empty bounded subset of Rn space, N(Ai) represents the i-th microseismic moment, for any one r>0, the minimum number of microseismic events in the spatial distribution is measured by a positive cube of length r, and n is the dimension of the spatial distribution pattern of microseismic events.

9. The method of claim 1, wherein, The formula for determining the rock mass damage degree evaluation index according to the microseismic time dimension value corresponding to the microseismic event occurring at each microseismic moment and the rock mass damage quantity is: For each microseismic moment, the sum of the rock mass damage quantities of all microseismic events occurring from the first microseismic moment to the microseismic moment is determined as the rock mass cumulative damage quantity corresponding to the microseismic moment; The microseismic time dimension value and the rock mass cumulative damage quantity of the microseismic moment form a rock mass damage matrix; For each microseismic moment, the rock mass damage degree evaluation index corresponding to the microseismic event occurring at the microseismic moment is determined according to the rock mass damage matrix corresponding to the microseismic moment.

10. The method of claim 9, wherein, The formula for determining the rock mass damage degree evaluation index corresponding to the microseismic event occurring at each microseismic moment according to the rock mass damage matrix corresponding to the microseismic moment is: ； ； wherein, M x is the rock mass damage matrix, microseismic time x corresponds to the time vector, is the microseismic time x corresponds to the microseismic time dimension value of the shear failure form, is the microseismic time x corresponds to the microseismic time dimension value of the tensile failure form, D x is the microseismic time x corresponds to the rock mass cumulative damage amount, C is the microseismic time x corresponds to the rock mass damage degree evaluation index of the microseismic event.

Images