Gridding analysis method for micro-seismic activity energy plane evolution

Abstract

The invention provides a gridding analysis method for micro-seismic activity energy plane evolution, and belongs to the technical field of coal mine safety. The method comprises the following steps: dynamically determining a grid division range based on a plane coordinate extreme value of a micro-seismic event in a target area, and dividing the area into rectangular grids with equal sizes through a preset grid side length parameter; counting an energy accumulated value of a micro-seismic event in each grid, and performing energy level conversion by adopting a logarithmic function to eliminate magnitude difference; and finally, a two-dimensional energy level distribution thermodynamic diagram is generated through an interpolation algorithm, and energy space aggregation characteristics and dynamic evolution laws are visually displayed. According to the method, by establishing the nonlinear mapping relation between the grid coordinates and the energy values and combining the thermodynamic diagram visualization technology, the problems that in a traditional method, energy distribution characterization is not sensitive, and the spatio-temporal evolution trend is difficult to recognize are effectively solved. According to the method, the quantitative precision and the dynamic monitoring efficiency of micro-seismic activity analysis are remarkably improved, and reliable technical support is provided for mine safety early warning.

Description

Technical Field

[0001] This invention relates to the field of coal mine safety mining technology, specifically a gridded analysis method for the planar evolution of microseismic activity energy. Background Technology

[0002] Rockbursts are often accompanied by microseismic events. When the seismic waves generated by these microseismic events propagate to the vicinity of the working face or roadway, they disturb the surrounding rock in the near-field of the mining space. Once the mechanical conditions for rockburst to occur are met, a rockburst will be induced. Therefore, understanding the occurrence patterns and spatial distribution of microseismic events is beneficial for taking preventative measures to avoid rockbursts.

[0003] Because numerous microseismic events are generated during the mining process, directly projecting them onto a base map can, to some extent, reflect their distribution characteristics, but it cannot intuitively reflect the concentrated areas and evolution trends of microseismic energy. Therefore, a characterization method that can reflect the evolution of planar microseismic activity is needed. Summary of the Invention

[0004] To address the aforementioned technical shortcomings, the purpose of this invention is to provide a gridded analysis method for the planar evolution of microseismic activity energy, so as to reflect the concentrated areas and evolution trends of microseismic energy.

[0005] To achieve the above objectives, the present invention adopts the following technical solution:

[0006] A gridded analysis method for the planar evolution of microseismic activity energy includes the following steps:

[0007] S1. Determine the grid division range based on the extreme values ​​of the microseismic event coordinates within the target area and time period, and divide the area into m rows and n columns of equal-sized rectangular grids, with the length and width of each grid being l.

[0008] S2. For each grid cell, count all microseismic events falling within that grid cell and calculate the total energy value E of that grid cell. ij =∑E k E k This represents the energy released by the k-th microseismic event within the grid.

[0009] S3. Perform energy level transformation on the total energy value of each grid, and calculate the energy level value e using a nonlinear mapping function. ij Energy level conversion involves logarithmic processing of energy values ​​to reduce the difference in numerical magnitude.

[0010] S4. Generate a two-dimensional energy distribution heat map based on the energy level values ​​of each grid cell, and visualize the spatial accumulation characteristics and temporal evolution of microseismic energy.

[0011] Preferably, in step S1, the method for determining the mesh division range specifically includes:

[0012] Calculate the minimum value x along the x-axis of all microseismic events within the target area. min and maximum value x max The minimum value of y in the y-axis direction min and maximum value y max This determines the grid division range of the target area; according to the formula m=(x max -x min ) / l+1, n=(y max -y min ) / l+1 determines the number of grid rows and columns, where l is the preset grid side length parameter.

[0013] Preferably, in step S2, the condition for determining whether a microseismic event falls within a specified grid is as follows: For a microseismic event with coordinates (x, y), it is determined to fall within the grid with center point coordinates (x0, y0) when the following formula is satisfied simultaneously:

[0014] x∈[x0-l / 2,x0+l / 2]

[0015] y∈[y0-l / 2,y0+l / 2]

[0016] Where l represents the length and width of a single grid cell; x0 = x min +(i-0.5)*l, y0=y min +(j-0.5)*l; i,j are the row and column indices of the grid, respectively.

[0017] Preferably, in step S3, the energy level conversion specifically adopts the following formula:

[0018] e ij =log 10 (E ij )

[0019] Among them, E ij This represents the total energy value of the target grid.

[0020] Preferably, step S4 includes:

[0021] A heat map of energy level distribution is generated using the difference method, and the energy release status of different regions is reflected by color mapping.

[0022] Compared with the prior art, the beneficial effects of the present invention are as follows:

[0023] This invention uses grid division to statistically analyze total energy and reduces numerical magnitude differences through logarithmic energy level conversion, making microseismic events of different energy scales easier to distinguish in visualization and solving the problem of insensitivity in energy distribution characterization.

[0024] The generated two-dimensional energy distribution heat map can clearly show the spatial accumulation characteristics of microseismic energy. At the same time, combined with data from different time periods, it can intuitively present the evolution law over time, overcoming the drawback of the difficulty in identifying the spatiotemporal evolution trend in traditional methods.

[0025] By establishing a nonlinear mapping relationship between grid coordinates and energy values, quantitative analysis of microseismic activity energy was achieved, significantly improving the accuracy of the analysis. The visualization method of the heat map makes the results more intuitive and easier to interpret quickly, improving the efficiency of dynamic monitoring.

[0026] The method has clear steps, from grid division, energy statistics, energy level conversion to heat map generation. The process is well-defined and highly operable, and can be directly applied to microseismic monitoring and analysis in actual mines, providing effective data support for preventing safety accidents such as rockbursts. Attached Figure Description

[0027] Figure 1 This is a flowchart of the identification process of the present invention;

[0028] Figure 2 This is a schematic diagram of microseismic meshing in this invention;

[0029] Figure 3 This is the micro-vibration projection in the embodiment of the present invention;

[0030] Figure 4 This is a micro-vibration energy level thermogram in an embodiment of the present invention. Detailed Implementation

[0031] The invention will now be further described with reference to the accompanying drawings.

[0032] like Figure 1 , Figure 2 As shown, a gridded analysis method for the planar evolution of microseismic activity energy includes the following steps:

[0033] S1. Determine the grid division range based on the extreme values ​​of the microseismic event coordinates within the target area and time period, and divide the area into m rows and n columns of equal-sized rectangular grids, with the length and width of each grid being l.

[0034] The specific methods for determining the grid division range include:

[0035] Calculate the minimum value x along the x-axis of all microseismic events within the target area. min and maximum value x max The minimum value of y in the y-axis direction min and maximum value y max This determines the grid division range of the target area; according to the formula m=(x max -x min) / l+1, n=(y max -y min ) / l+1 determines the number of grid rows and columns, where l is the preset grid side length parameter.

[0036] S2. For each grid cell, count all microseismic events falling within that grid cell and calculate the total energy value E of that grid cell. ij =∑E k E k This represents the energy released by the k-th microseismic event within the grid.

[0037] The condition for determining whether a microseismic event falls within a specified grid is as follows: For a microseismic event with coordinates (x, y), it is determined to fall within the grid with center point coordinates (x0, y0) if the following formula is satisfied simultaneously:

[0038] x∈[x0-l / 2,x0+l / 2]

[0039] y∈[y0-l / 2,y0+l / 2]

[0040] Where l represents the length and width of a single grid cell; x0 = x min +(i-0.5)*l, y0=y min +(j-0.5)*l; i,j are the row and column indices of the grid, respectively.

[0041] S3. Perform energy level transformation on the total energy value of each grid, and calculate the energy level value e using a nonlinear mapping function. ij Energy level conversion involves logarithmic processing of energy values ​​to reduce the difference in numerical magnitude.

[0042] The energy level conversion uses the following formula:

[0043] e ij =log 10 (E ij )

[0044] Among them, E ij This represents the total energy value of the target grid.

[0045] S4. A two-dimensional energy distribution heatmap is generated based on the energy level values ​​of each grid cell, visually displaying the spatial accumulation characteristics and temporal evolution of microseismic energy. The energy level distribution heatmap is generated using the difference method, and color mapping reflects the energy release state of different regions.

[0046] Example:

[0047] This method was used to implement a gridded analysis method for the planar evolution of microseismic activity energy in the 250108-1 working face of a certain mine. The specific steps are as follows:

[0048] 1. Collect microseismic data for the 250108-1 working face in February 2025, such as... Figure 3 As shown.

[0049] 2. Calculate the minimum value x along the x-axis for all microseismic events within the target area. min and maximum value x max The minimum value of y in the y-axis direction min and maximum value y max This determines the grid division range of the target area. Figure 3 It can be seen that in this embodiment, the minimum value x in the x-axis direction is... min The value is 374516.53, and the maximum value is x. max The value is 374939.24; the minimum value of y in the y-axis direction is... min The value is 3902592.85, and the maximum value is y. max The value is 3903491.42.

[0050] 3. Collect all microseismic events falling within each grid and calculate the total energy value E of that grid. ij And according to formula e ij =log 10 (E ij Calculate the energy levels of the grid.

[0051] 4. Generate a heat map of energy level distribution using the difference method. Figure 4 The color chart reflects the energy release state in different regions; as shown in the figure... Figure 4 Compare Figure 3 It is more intuitive and better reflects the characteristics and concentration areas of energy release.

[0052] Through innovative gridded analysis and visualization technology, the shortcomings of traditional microseismic energy analysis methods are effectively addressed. It has significant advantages in terms of sensitivity, identification ability, accuracy, efficiency and practicality, providing strong technical support for safe mining in fields such as coal mining.

Claims

1. A gridded analysis method for the planar evolution of microseismic activity energy, characterized in that, Includes the following steps: S1. Determine the grid division range based on the extreme values ​​of the microseismic event coordinates within the target area and time period, and divide the area into m rows and n columns of equal-sized rectangular grids, with the length and width of each grid being l. S2. For each grid cell, count all microseismic events falling within that grid cell and calculate the total energy value E of that grid cell. ij =∑E k E k This represents the energy released by the k-th microseismic event within the grid. S3. Perform energy level transformation on the total energy value of each grid, and calculate the energy level value e using a nonlinear mapping function. ij Energy level conversion involves logarithmic processing of energy values ​​to reduce the difference in numerical magnitude. S4. Generate a two-dimensional energy distribution heat map based on the energy level values ​​of each grid cell, and visualize the spatial accumulation characteristics and temporal evolution of microseismic energy.

2. The gridded analysis method for the planar evolution of microseismic activity energy as described in claim 1, characterized in that, In step S1, the method for determining the mesh division range specifically includes: Calculate the minimum value x along the x-axis of all microseismic events within the target area. min and maximum value x max The minimum value of y in the y-axis direction min and maximum value y max This determines the grid division range of the target area; according to the formula m=(x max -x min ) / l+1, n=(y max -y min ) / l+1 determines the number of grid rows and columns, where l is the preset grid side length parameter.

3. The gridded analysis method for the planar evolution of microseismic activity energy as described in claim 1, characterized in that, In step S2, the condition for determining whether a microseismic event falls within a specified grid is as follows: For a microseismic event with coordinates (x, y), it is determined to fall within the grid with center point coordinates (x0, y0) when the following formula is satisfied simultaneously: x∈[x0-l / 2,x0+l / 2] y∈[y0-l / 2,y0+l / 2] Where l represents the length and width of a single grid cell; x0 = x min +(i-0.5)*l, y0=y min +(j-0.5)*l; i,j are the row and column indices of the grid, respectively.

4. The gridded analysis method for the planar evolution of microseismic activity energy as described in claim 1, characterized in that, In step S3, the energy level conversion specifically adopts the following formula: and ij log 10 (AND ij ) Among them, E ij This represents the total energy value of the target grid.

5. The gridded analysis method for the planar evolution of microseismic activity energy as described in claim 1, characterized in that, Step S4 includes: A heat map of energy level distribution is generated using the difference method, and the energy release status of different regions is reflected by color mapping.

Images