Mine seismic intensity and supporting structure performance evaluation method, device, equipment and medium

Abstract

The application discloses a mine shock intensity and supporting structure performance evaluation method and device, equipment and medium, relates to the technical field of supporting structure evaluation, and comprises the following steps: determining the peak particle vibration velocity based on the coal mine geological survey data and the particle peak vibration velocity determination model; determining the coal mine shock intensity grade, determining the load amplification coefficient based on the coal mine shock intensity grade, correcting the mine shock dynamic load vector through the load amplification coefficient, and obtaining the corrected mine shock dynamic load vector; correcting the stiffness matrix, and determining the target supporting-rock coupling dynamics model according to the corrected stiffness matrix; determining the structure response information of the target part based on the corrected mine shock dynamic load vector and the target supporting-rock coupling dynamics model; determining the comprehensive damage degree through the structure response information of the target part and the coal mine shock intensity grade, and determining the performance safety grade of the supporting structure according to the comprehensive damage degree. The mine shock intensity is accurately calculated, and the dynamic response quantitative analysis of the supporting structure is realized.

Description

Technical Field

[0001] This invention relates to the field of support structure assessment technology, and in particular to methods, devices, equipment and media for assessing the relationship between seismic intensity and support structure performance. Background Technology

[0002] Mining tremors are non-natural earthquake phenomena induced by mining activities. They are essentially a dynamic process in which the elastic energy of rock mass is rapidly released and propagates in the form of seismic waves. Strong mining tremors can easily induce disasters such as underground rockbursts, tunnel collapses, and support failures, seriously threatening the lives of underground workers and even causing public safety problems such as damage to surface buildings and surface subsidence.

[0003] Existing methods for assessing mine seismic intensity largely follow the natural earthquake intensity classification system, combining macroscopic surveys and instrumental monitoring results to determine the intensity level. The assessment dimensions focus on building damage, surface deformation, and human perception, without specifically adapting to the stress characteristics and failure mechanisms of underground mine support structures. Some mine seismic assessment techniques use peak particle velocity (PPV) as the core indicator, determining the intensity level through pre-defined mapping relationships, but fail to consider the dynamic coupling effect of seismic vibrations with the support structure and surrounding rock, thus failing to accurately quantify the real-time dynamic response and remaining bearing capacity of the support structure.

[0004] Furthermore, existing technologies rely solely on a single vibration velocity parameter for assessment. Limited by the monitoring equipment's range, this leads to data distortion under strong mining seismic conditions, resulting in significant errors in the assessment results. Moreover, the lack of linkage analysis between seismic intensity and support structure damage hinders support design optimization and safe mining decisions, failing to meet the safety management needs of deep mining and seismically active areas. Therefore, accurately calculating seismic intensity and quantitatively analyzing the dynamic response of support structures to provide a scientific basis for optimizing mine support structures and preventing seismic disasters is a pressing issue that needs to be addressed. Summary of the Invention

[0005] In view of this, the purpose of this invention is to provide a method, device, equipment, and medium for evaluating the intensity of mine seismic events and the performance of support structures, enabling accurate calculation of seismic intensity and quantitative analysis of the dynamic response of support structures, thus providing a scientific basis for optimizing mine support structures and preventing mine seismic disasters. The specific solution is as follows: In the first aspect, this application discloses a method for evaluating the intensity of mine earthquakes and the performance of support structures, including: Collect geological survey data and support structure design parameters of the coal mine area, and determine the peak vibration velocity of the particles based on the geological survey data of the coal mine area and the pre-constructed peak vibration velocity determination model. Based on the pre-constructed mapping relationship, the coal mine seismic intensity level corresponding to the peak vibration velocity of the particle is determined. Based on the coal mine seismic intensity level, the load amplification factor is determined. The mine seismic load vector in the initial support-surrounding rock coupled dynamic model constructed based on the support structure design parameters is corrected by the load amplification factor to obtain the corrected mine seismic load vector. The stiffness matrix in the initial support-surrounding rock coupled dynamic model is modified so as to determine the target support-surrounding rock coupled dynamic model based on the modified stiffness matrix. The structural response information of the target location is determined based on the modified mine vibration load vector and the target support-surrounding rock coupled dynamic model; the target location includes anchor bolts / cables, shotcrete, steel supports and surrounding rock; the structural response information includes any one or a combination of stress, strain, plastic deformation, crack width, damage factor, buckling degree, residual deformation, plastic zone range and displacement. The comprehensive damage degree of the support structure is determined by the structural response information of the target location and the seismic intensity level of the coal mine, and the performance safety level of the support structure is determined based on the comprehensive damage degree.

[0006] Optionally, the collection of geological survey data and support structure design parameters for the coal mine area includes: The geological survey data and support structure design parameters of the coal mine area are collected in real time by an array of vibration sensors deployed in underground roadways, mining areas and on the surface. The geological survey data of the coal mine area includes the horizontal distance from the monitoring point to the seismic source and the equivalent source depth of the mine earthquake. The support structure design parameters include the specifications of the anchor bolts / cables, the strength of the sprayed layer, the model of the steel support, the support density and the anchoring force.

[0007] Optionally, before determining the peak particle vibration velocity based on the geological exploration data of the coal mine area and the pre-constructed peak particle vibration velocity determination model, the method further includes: Screening for mine tremor events that meet preset conditions that occur in coal mining areas within a target time period; Based on the parameter information corresponding to the mining earthquake event, a mathematical model is fitted using the least squares method. The accuracy of the mathematical model is verified using the leave-one-out cross-validation method. Based on the corresponding verification results, the peak vibration velocity of the mass point is determined to confirm the model. Accordingly, determining the peak particle vibration velocity based on the geological exploration data of the coal mine area and a pre-constructed peak particle vibration velocity determination model includes: Based on the geological exploration data of the coal mine area and the pre-constructed peak particle vibration velocity determination model, the peak particle vibration velocity is determined; the peak particle vibration velocity determination model is as follows: ; Where Y is the peak velocity of the particle; M is the magnitude of the seismic event; r is the horizontal distance from the monitoring point to the seismic source; h0 is the equivalent source depth of the seismic event; S is the site category coefficient; G is the geological structure influence coefficient; a, b, c, d and e are undetermined coefficients for model fitting; and ε is the random error term.

[0008] Optionally, before determining the coal mine seismic intensity level corresponding to the peak particle vibration velocity based on the pre-constructed mapping relationship, the method further includes: Based on enterprise standards for monitoring and evaluating seismic intensity in coal mines, types of underground coal roadway support, damage characteristics, and measured parameters from instruments, a mapping relationship between seismic intensity levels in coal mines and peak particle vibration velocities is established.

[0009] Optionally, before correcting the mine vibration load vector in the initial support-surrounding rock coupled dynamic model constructed based on the support structure design parameters using the load amplification factor, the method further includes: An initial support-surrounding rock coupled dynamic model is constructed based on the support structure design parameters; the initial support-surrounding rock coupled dynamic model is as follows: ; Among them, [M] s For the support-surrounding rock coupling consistent quality matrix; [C] s Let K be the Rayleigh damping matrix; s (u) is the stiffness matrix; F(t) is the mine vibration load vector; , , These are the nodal acceleration, velocity, and displacement vectors, respectively. Accordingly, the load amplification factor is determined based on the seismic intensity level of the coal mine, and the seismic load vector in the initial support-surrounding rock coupled dynamic model constructed based on the support structure design parameters is corrected using the load amplification factor to obtain the corrected seismic load vector, including: Determine the first product between the earthquake intensity level of the coal mine and the first value; The sum of the first product and the second value is determined as the load amplification factor; The product of the load amplification factor and the third value is determined as the corrected mine vibration load vector.

[0010] Optionally, the modification of the stiffness matrix in the initial support-surrounding rock coupled dynamic model includes: The comprehensive geological condition coefficient is determined based on the weighted average of the rock strength coefficient and the structural complexity coefficient. Determine the second product between the comprehensive geological condition coefficient and the fourth value; Determine the target sum between 1 and the second product; The product of the target and the stiffness matrix in the initial support-surrounding rock coupled dynamic model is determined as the corrected stiffness matrix.

[0011] Optionally, determining the comprehensive damage degree of the support structure based on the structural response information of the target location and the seismic intensity level of the coal mine includes: Determine the comprehensive geological condition coefficient and the support structure characteristic coefficient based on the structural response information of the target location; The damage degree of the target location is determined using a damage degree determination formula based on the comprehensive geological condition coefficient, the support structure characteristic coefficient, and the seismic intensity level of the coal mine; the damage degree determination formula is as follows: ; Wherein, D represents the degree of damage; I represents the intensity level of the coal mine seismic event; G represents the comprehensive geological condition coefficient; S represents the characteristic coefficient of the support structure; and P represents the effective duration of the seismic event. The weighted average of the damage levels at each target location is determined as the overall damage level of the support structure.

[0012] Optionally, determining the performance safety level of the support structure based on the comprehensive damage degree includes: The performance safety level of the support structure is determined based on the comparison between the comprehensive damage degree and the preset threshold range; wherein, different preset threshold ranges correspond to different performance safety levels.

[0013] Secondly, this application discloses a device for evaluating the intensity of mine tremors and the performance of support structures, comprising: The peak vibration velocity determination module is used to collect geological exploration data of the coal mine area and design parameters of the support structure, and determine the peak vibration velocity of the particles based on the geological exploration data of the coal mine area and the pre-constructed peak vibration velocity determination model. The correction module is used to determine the coal mine seismic intensity level corresponding to the peak vibration velocity of the mass point based on the pre-constructed mapping relationship, determine the load amplification factor based on the coal mine seismic intensity level, and correct the mine vibration load vector in the initial support-surrounding rock coupled dynamic model based on the support structure design parameters through the load amplification factor to obtain the corrected mine vibration load vector. The model determination module is used to correct the stiffness matrix in the initial support-surrounding rock coupled dynamic model, so as to determine the target support-surrounding rock coupled dynamic model based on the corrected stiffness matrix. The structural response information determination module is used to determine the structural response information of the target location based on the modified mine vibration load vector and the target support-surrounding rock coupled dynamic model; the target location includes anchor bolts / cables, shotcrete layer, steel support and surrounding rock; the structural response information includes any one or a combination of stress, strain, plastic deformation, crack width, damage factor, buckling degree, residual deformation, plastic zone range and displacement. The performance and safety level determination module is used to determine the comprehensive damage degree of the support structure by using the structural response information of the target location and the seismic intensity level of the coal mine, and to determine the performance and safety level of the support structure based on the comprehensive damage degree.

[0014] Thirdly, this application discloses an electronic device, including: Memory, used to store computer programs; A processor for executing computer programs to implement the steps of the aforementioned method for evaluating the intensity of seismic events and the performance of support structures.

[0015] Fourthly, this application discloses a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the steps of the aforementioned method for evaluating the intensity of seismic events and the performance of support structures.

[0016] This application collects geological survey data and support structure design parameters of coal mine areas. Based on the geological survey data and a pre-constructed peak particle vibration velocity determination model, it determines the peak particle vibration velocity. Based on a pre-constructed mapping relationship, it determines the coal mine seismic intensity level corresponding to the peak particle vibration velocity. Based on the coal mine seismic intensity level, it determines the load amplification factor. The load amplification factor is used to correct the seismic load vector in the initial support-surrounding rock coupled dynamic model constructed based on the support structure design parameters, resulting in a corrected seismic load vector. The stiffness matrix in the initial support-surrounding rock coupled dynamic model is also corrected according to... The modified stiffness matrix determines the coupled dynamic model of the target support-surrounding rock. Based on the modified mine seismic load vector and the coupled dynamic model, the structural response information of the target location is determined. The target location includes anchor bolts / cables, shotcrete, steel supports, and surrounding rock. The structural response information includes any one or a combination of stress, strain, plastic deformation, crack width, damage factor, buckling degree, residual deformation, plastic zone range, and displacement. The comprehensive damage degree of the support structure is determined by the structural response information of the target location and the mine seismic intensity level. The performance safety level of the support structure is determined based on the comprehensive damage degree. Therefore, this application achieves deep coupling between mine seismic intensity and the dynamic response of the support structure, breaking through the limitations of traditional single-parameter assessment, accurately quantifying the damage degree of the support structure, and significantly improving assessment accuracy. It realizes accurate calculation of mine seismic intensity, quantitative analysis of the dynamic response of the support structure, and comprehensive damage degree classification, providing a scientific basis for mine support structure optimization and mine seismic disaster prevention. Attached Figure Description

[0017] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on the provided drawings without creative effort.

[0018] Figure 1 This application discloses a flowchart for evaluating the intensity of mine earthquakes and the performance of support structures. Figure 2 This is a schematic diagram of a device for evaluating the intensity of mine earthquakes and the performance of support structures disclosed in this application. Figure 3 This is a structural diagram of an electronic device disclosed in this application. Detailed Implementation

[0019] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0020] Existing technologies rely solely on a single vibration velocity parameter for assessment. Limited by the monitoring equipment's range, this leads to data distortion under strong mining seismic conditions, resulting in significant errors in the assessment results. Furthermore, the lack of linkage analysis between seismic intensity and support structure damage hinders support design optimization and safe mining decisions, failing to meet the safety management needs of deep mining and seismically active areas. To address these technical problems, this application discloses a method, apparatus, equipment, and medium for assessing seismic intensity and support structure performance. This method enables accurate calculation of seismic intensity and quantitative analysis of the dynamic response of support structures, providing a scientific basis for optimizing mine support structures and preventing seismic disasters.

[0021] See Figure 1 As shown in the figure, an embodiment of the present invention discloses a method for evaluating the intensity of mine earthquakes and the performance of support structures, including: Step S11: Collect geological survey data of the coal mine area and design parameters of the support structure. Based on the geological survey data of the coal mine area and the pre-constructed peak vibration velocity determination model, determine the peak vibration velocity of the particles.

[0022] In this embodiment, a vibration sensor array deployed in underground roadways, mining areas, and on the surface collects real-time geological survey data and support structure design parameters of the coal mine area. The geological survey data includes the horizontal distance from the monitoring point to the seismic source and the equivalent seismic source depth. The support structure design parameters include anchor bolt / cable specifications, shotcrete intensity, steel support type, support density, and anchoring force. Specifically, the vibration sensor array deployed in underground roadways, mining areas, and on the surface collects real-time seismic acceleration, velocity, and displacement time-history data; simultaneously, it collects geological survey data of the mining area, including surrounding rock mechanical parameters, geological structure distribution, and stratum attitude; and collects support structure design parameters and measured data, covering anchor bolt / cable specifications, shotcrete intensity, steel support type, support density, and anchoring force. The collected seismic data undergoes filtering, noise reduction, and baseline correction preprocessing to remove environmental noise and equipment interference signals, retaining effective vibration waveform data to provide an accurate data source for subsequent model calculations.

[0023] Before determining the peak particle vibration velocity, seismic events occurring in the coal mine area within the target time period that meet preset conditions are screened. Based on the parameter information corresponding to the seismic events, a mathematical model is fitted using the least squares method. The accuracy of the mathematical model is verified using leave-one-out cross-validation. The peak particle vibration velocity determination model is determined based on the verification results. Correspondingly, the peak particle vibration velocity is determined based on the geological exploration data of the coal mine area and the pre-constructed peak particle vibration velocity determination model, including: determining the peak particle vibration velocity based on the geological exploration data of the coal mine area and the pre-constructed peak particle vibration velocity determination model; the peak particle vibration velocity determination model is: ; Wherein, Y is the peak vibration velocity of the particle, in cm / s, and is the core assessment parameter for the destructive effect of the mine earthquake; M is the magnitude of the mine earthquake, characterizing the amount of energy released; r is the horizontal distance from the monitoring point to the source; h0 is the equivalent source depth of the mine earthquake; S is the site category coefficient, calibrated based on mechanical test data such as the elastic modulus and wave velocity of the rock strata in the mining area; G is the geological structure influence coefficient, assigned based on the exploration results such as the scale of faults, the degree of development of fracture zones, and the density of fractures in the mining area; a, b, c, d, and e are undetermined coefficients for model fitting, determined through regression analysis of measured data; ε is the random error term, following a normal distribution N(0, σ). 2 ).

[0024] In one specific embodiment, valid seismic samples from the mining area over a period of 5-10 years are selected. The selection criteria are: magnitude M ≥ 1.0, complete and anomaly-free monitoring data, and matching with corresponding geological and site parameters. Parameters such as magnitude and distance are standardized and converted to eliminate dimensional differences. The least squares method is used to fit the model coefficients, and irrelevant variables with a significance level α > 0.05 are eliminated through stepwise regression analysis. The model accuracy is verified using leave-one-out cross-validation, requiring the relative error between predicted and measured values ​​to be ≤ 15%. If this is not met, additional samples are added and the model is refitted. The applicable range of the model is defined as: magnitude 0.5 ≤ M ≤ 4.0, monitoring distance 0.1 km ≤ r ≤ 5.0 km, and focal depth 0.05 km ≤ h0 ≤ 1.5 km.

[0025] Step S12: Determine the coal mine seismic intensity level corresponding to the peak vibration velocity of the mass point based on the pre-constructed mapping relationship, determine the load amplification factor based on the coal mine seismic intensity level, and correct the mine vibration load vector in the initial support-surrounding rock coupled dynamic model constructed based on the support structure design parameters through the load amplification factor to obtain the corrected mine vibration load vector.

[0026] In this embodiment, a mapping relationship between the seismic intensity level and the peak velocity of mass points in a coal mine is established based on enterprise standards for coal mine seismic intensity monitoring and evaluation, underground coal roadway support types, damage characteristics, and instrument-measured parameters. Simultaneously, the criteria for personnel perception and support damage assessment are clarified, with specific classifications shown in Table 1. Table 1. Mapping Relationship Diagram

[0027] In this embodiment, before correcting the mine vibration load vector in the initial support-surrounding rock coupled dynamic model constructed based on the support structure design parameters using the load amplification factor, an initial support-surrounding rock coupled dynamic model is constructed based on the support structure design parameters; the initial support-surrounding rock coupled dynamic model is as follows: ; Among them, [M] s For the support-surrounding rock coupling consistent quality matrix; [C] s Let K be the Rayleigh damping matrix; s (u) is the stiffness matrix; F(t) is the mine vibration load vector; , , These represent the nodal acceleration, velocity, and displacement vector, respectively.

[0028] Then, determine the first product between the coal mine seismic intensity level and the first value; determine the sum between the first product and the second value as the load amplification factor; and determine the product between the load amplification factor and the third value as the corrected mine seismic load vector.

[0029] In one specific embodiment, a coupled dynamic model of the support-surrounding rock is established using the finite element method. This equation characterizes the instantaneous dynamic equilibrium relationship of the support-surrounding rock coupled system, i.e., the resultant force of inertial force, damping force, and elastic restoring force equals the mine vibration load. The modeling and solution requirements are as follows: Element selection: C3D8R (three-dimensional eight-node reduced integral solid element) solid elements are used for the surrounding rock; T3D2 truss elements are used for the anchor bolts / cables; S4R shell elements are used for the shotcrete; and B31 beam elements are used for the steel support. Material constitutive model: Drucker-Prager elastoplastic model is used for the surrounding rock; bilinear elastoplastic model is used for the anchor bolts / cables; and concrete damage plastic model is used for the shotcrete. Contact boundary: The support and surrounding rock are in surface-to-surface contact; Cohesive elements are used for anchor bolt bonding; the model is fixed at the bottom, horizontally constrained on the sides, and the self-weight stress of the overburden is applied at the top. Solution process: After completing geometric modeling and mesh generation, the mine vibration load is applied, and explicit dynamic solution is used to extract stress, strain, and displacement data to analyze the structural damage evolution law. Intensity-load coupling: Calculating the load amplification factor k I=0.2I+0.1, corrected dynamic load Where I represents the seismic intensity level of the coal mine.

[0030] Step S13: Correct the stiffness matrix in the initial support-surrounding rock coupled dynamic model, so as to determine the target support-surrounding rock coupled dynamic model based on the corrected stiffness matrix.

[0031] In this embodiment, a comprehensive geological condition coefficient is determined based on a weighted average of the rock strength coefficient and the structural complexity coefficient; a second product is determined between the comprehensive geological condition coefficient and a fourth value; a target sum is determined between 1 and the second product; and the product of the target sum and the stiffness matrix in the initial support-surrounding rock coupled dynamic model is determined as the corrected stiffness matrix. The worse the geological conditions (e.g., fault development, fractured surrounding rock), the greater the stiffness reduction, and the more easily the structure deforms. In a specific embodiment, geological-structural coupling: corrected stiffness matrix. The model is adapted to the geological fracture characteristics. Here, G: comprehensive geological condition coefficient, G = 0.6G1 + 0.4G2; G1: rock strength coefficient, assigned a value of 0~1 based on test data such as the uniaxial compressive strength and elastic modulus of the surrounding rock, according to industry grading standards; G2: structural complexity coefficient, assigned a value of 0~1 based on the degree of development of faults, fracture zones, and fissures. Based on this, the target support-surrounding rock coupled dynamic model is determined according to the modified stiffness matrix.

[0032] Step S14: Determine the structural response information of the target location based on the corrected mine vibration load vector and the target support-surrounding rock coupled dynamic model; the target location includes anchor bolts / cables, shotcrete layer, steel support and surrounding rock; the structural response information includes any one or a combination of stress, strain, plastic deformation, crack width, damage factor, buckling degree, residual deformation, plastic zone range and displacement.

[0033] In this embodiment, the corrected dynamic load F'(t) is a time-varying vector; it is not a fixed force, but rather a series of force values ​​at different times. This load is not applied to a single point on the model, but rather to all vibration-related nodes. The common practice is to start from the pre-processed acceleration time history; multiply by the mass corresponding to the node; obtain the load value for each node at each time step; and apply these loads as "external forces" to the corresponding nodes of the model. This is equivalent to "pushing" and "pulling" the model at each node on the bottom or side according to the time sequence, simulating the force of the vibration wave propagating to the tunnel. The solver, following the explicit dynamics method, cuts the vibration process into very small time steps (e.g., 0.1 milliseconds per step), pushing forward step by step. After the calculation, the results of each time step are recorded. The post-processing stage extracts useful information from this massive amount of data: the response of the anchor bolts / cables: extracting the maximum principal stress σ1_max of each anchor bolt; extracting the plastic strain of each anchor bolt. (If entering the plastic phase); check the stress-time curve to determine if yielding has occurred. Shotcrete response: extract the concrete damage factor j (automatically output by the damage plasticity model), the maximum value j_max represents the microcrack density; extract the crack width w (calculated through the displacement difference between adjacent nodes); check if any elements have been deleted (representing crushing or cracking). Steel support response: extract the maximum compressive stress σ3_max; determine if buckling has occurred through displacement (sudden increase in lateral displacement); calculate the residual deformation δ_residual after the vibration ends (final displacement minus the elastic recovery portion). Surrounding rock response: extract the equivalent plastic strain ε_eq_plastic, the plastic zone is the area where ε_eq_plastic>0; extract the maximum displacement u_max of key nodes; check if the plastic zone is continuous (extending from the tunnel surface to depth). In summary, the modified dynamic load F'(t) is applied as an "external force" to the support-surrounding rock model with the modified stiffness matrix [K]'s. The dynamic equations are solved frame by frame in time step by explicit dynamic method, thereby calculating the structural response data such as stress, strain, displacement, and damage factor at each moment.

[0034] Step S15: Determine the comprehensive damage degree of the support structure through the structural response information of the target location and the seismic intensity level of the coal mine, and determine the performance safety level of the support structure based on the comprehensive damage degree.

[0035] In this embodiment, a comprehensive geological condition coefficient and a support structure characteristic coefficient are determined based on the structural response information of the target location; the damage degree of the target location is determined using a damage degree determination formula based on the comprehensive geological condition coefficient, the support structure characteristic coefficient, and the coal mine seismic intensity level; the damage degree determination formula is: ; Wherein, D represents the damage degree, ranging from 0 to 1, with a larger value indicating more severe damage; I represents the seismic intensity level of the coal mine, with levels I to X corresponding to 1 to 10; G represents the comprehensive geological condition coefficient; S represents the support structure characteristic coefficient, S = 0.7S1 + 0.3S2; S1: support strength coefficient, assigned a value of 0 to 1 based on the measured values ​​of anchor bolt / anchor cable yield strength, spray layer grade, and support stiffness; S2: structural integrity coefficient, assigned a value of 0 to 1 based on the support failure rate and the number of failed components. P represents the effective vibration duration of the seismic event, extracting the duration of the vibration signal exceeding the threshold; the weighted average of the damage degrees of each target location is determined as the comprehensive damage degree of the support structure.

[0036] Damage is calculated separately for four components: anchor bolts / cables, shotcrete layer, steel support, and surrounding rock. The damage degree of each component is determined by comparing the calculated stress / strain with the material's ultimate limit value—0 for intact, 1 for complete failure, and intermediate values ​​are taken proportionally. A weighted sum is then used to obtain the overall damage degree D. z The damage of the four sub-items is summed by weight. The default weights are 0.3 for anchor cables, 0.2 for shotcrete, 0.3 for steel supports, and 0.2 for surrounding rock. These weights can be adjusted according to the actual support type.

[0037] Finally, the performance safety level of the support structure is determined based on the comparison between the comprehensive damage level and the preset threshold range; different preset threshold ranges correspond to different performance safety levels. Safety level determination criteria: Level I (Green Safety): 0 ≤ D z <0.3: The support system is intact with no obvious damage, and normal mining operations can proceed; Level II (Yellow Alert): 0.3≤D z <0.6: Localized damage to the support structure, decreased load-bearing capacity, requiring increased monitoring and localized reinforcement; Level III (Red Control): D z ≥0.6: The support has seriously failed, posing a major safety hazard. Production must be stopped immediately, personnel evacuated, and comprehensive repairs carried out.

[0038] In addition, the following recommendations are made for support optimization and mining: Safety level: Carry out preventive support optimization, increase the spacing of anchor cables, add linkage components, and improve the overall seismic resistance; Warning level: For damaged and weak areas, install additional anchor cables, repair the shotcrete layer, grout to reinforce the surrounding rock, and slow down the mining speed; Control level: Remove failed support, adopt a high-strength combined support system, adjust the mining layout to avoid high-risk areas, and resume operations only after acceptance.

[0039] It should be noted that in this application, distributed fiber optic sensing technology can be used to replace traditional sensors in the data acquisition stage to achieve distributed monitoring of mine seismic activity; mine seismic intensity assessment can employ machine learning algorithms such as deep learning and random forests, relying on massive amounts of data to build predictive models; and the support structure response analysis can use the discrete element method, adapted to the simulation of rock masses with well-developed joints and fractures. The above solutions should be selected comprehensively based on the mine site conditions, cost budget, and computational efficiency.

[0040] A coal mine transport roadway prone to frequent seismic activity was selected as an engineering case study. An assessment was conducted according to the technical process of this invention: 1) Six vibration sensors were deployed to collect and preprocess seismic data; 2) The seismic intensity assessment model coefficients for the mining area were fitted, the PPV value was calculated, and the intensity was determined; 3) A finite element model of roadway support and surrounding rock was established, and seismic loads were imported to conduct dynamic analysis; 4) The comprehensive damage degree D was calculated. z =0.42, judged as a yellow warning level; 5) Take control measures such as localized anchor cable installation, repair of sprayed layer, and slowing down tunneling speed to effectively eliminate safety hazards. In this way, it adapts to the geological and support characteristics of underground mines, constructs a dedicated assessment model and grading standard, which is highly targeted and applicable to various mines with frequent mine earthquakes; it forms a closed-loop control system of "monitoring-assessment-disposal-optimization", which can directly guide the optimization of support design and safe mining decisions, effectively reduce the risk of mine earthquake disasters, and ensure safe production in mines.

[0041] In summary, this application collects geological exploration data and support structure design parameters of the coal mine area. Based on the geological exploration data and a pre-constructed peak particle vibration velocity determination model, the peak particle vibration velocity is determined. Based on the pre-constructed mapping relationship, the corresponding coal mine seismic intensity level is determined. Based on the coal mine seismic intensity level, a load amplification factor is determined. The load amplification factor is used to correct the seismic load vector in the initial support-surrounding rock coupled dynamic model constructed based on the support structure design parameters, resulting in a corrected seismic load vector. The stiffness matrix in the initial support-surrounding rock coupled dynamic model is also corrected. The target support-surrounding rock coupled dynamic model is determined based on the modified stiffness matrix. The structural response information of the target location is determined based on the modified mine seismic load vector and the target support-surrounding rock coupled dynamic model. The target location includes anchor bolts / cables, shotcrete, steel supports, and surrounding rock. The structural response information includes any one or a combination of stress, strain, plastic deformation, crack width, damage factor, buckling degree, residual deformation, plastic zone range, and displacement. The comprehensive damage degree of the support structure is determined through the structural response information of the target location and the mine seismic intensity level. The performance safety level of the support structure is then determined based on the comprehensive damage degree. Therefore, this application achieves deep coupling between mine seismic intensity and the dynamic response of the support structure, breaking through the limitations of traditional single-parameter assessment, accurately quantifying the damage degree of the support structure, and significantly improving assessment accuracy. It realizes accurate calculation of mine seismic intensity, quantitative analysis of the dynamic response of the support structure, and comprehensive damage degree classification, providing a scientific basis for the optimization of mine support structures and the prevention and control of mine seismic disasters.

[0042] See Figure 2 As shown in the figure, an embodiment of the present invention discloses a device for evaluating the intensity of mine earthquakes and the performance of support structures, comprising: The peak vibration velocity determination module 11 is used to collect geological exploration data of the coal mine area and design parameters of the support structure, and determine the peak vibration velocity of the particles based on the geological exploration data of the coal mine area and the pre-constructed peak vibration velocity determination model. The correction module 12 is used to determine the coal mine seismic intensity level corresponding to the peak vibration velocity of the mass point based on the pre-constructed mapping relationship, determine the load amplification factor based on the coal mine seismic intensity level, and correct the mine vibration load vector in the initial support-surrounding rock coupled dynamic model based on the support structure design parameters through the load amplification factor to obtain the corrected mine vibration load vector. Model determination module 13 is used to correct the stiffness matrix in the initial support-surrounding rock coupled dynamic model, so as to determine the target support-surrounding rock coupled dynamic model based on the corrected stiffness matrix. The structural response information determination module 14 is used to determine the structural response information of the target part based on the modified mine vibration load vector and the target support-surrounding rock coupled dynamic model; the target part includes anchor bolts / cables, shotcrete layer, steel support and surrounding rock; the structural response information includes any one or a combination of stress, strain, plastic deformation, crack width, damage factor, buckling degree, residual deformation, plastic zone range and displacement. The performance safety level determination module 15 is used to determine the comprehensive damage degree of the support structure by using the structural response information of the target location and the seismic intensity level of the coal mine, and to determine the performance safety level of the support structure based on the comprehensive damage degree.

[0043] Since the embodiments of the device part correspond to the embodiments described above, please refer to the embodiments described in the method part for the embodiments of the device part, and will not be repeated here.

[0044] In summary, this application collects geological exploration data and support structure design parameters of the coal mine area. Based on the geological exploration data and a pre-constructed peak particle vibration velocity determination model, the peak particle vibration velocity is determined. Based on the pre-constructed mapping relationship, the corresponding coal mine seismic intensity level is determined. Based on the coal mine seismic intensity level, a load amplification factor is determined. The load amplification factor is used to correct the seismic load vector in the initial support-surrounding rock coupled dynamic model constructed based on the support structure design parameters, resulting in a corrected seismic load vector. The stiffness matrix in the initial support-surrounding rock coupled dynamic model is also corrected. The target support-surrounding rock coupled dynamic model is determined based on the modified stiffness matrix. The structural response information of the target location is determined based on the modified mine seismic load vector and the target support-surrounding rock coupled dynamic model. The target location includes anchor bolts / cables, shotcrete, steel supports, and surrounding rock. The structural response information includes any one or a combination of stress, strain, plastic deformation, crack width, damage factor, buckling degree, residual deformation, plastic zone range, and displacement. The comprehensive damage degree of the support structure is determined through the structural response information of the target location and the mine seismic intensity level. The performance safety level of the support structure is then determined based on the comprehensive damage degree. Therefore, this application achieves deep coupling between mine seismic intensity and the dynamic response of the support structure, breaking through the limitations of traditional single-parameter assessment, accurately quantifying the damage degree of the support structure, and significantly improving assessment accuracy. It realizes accurate calculation of mine seismic intensity, quantitative analysis of the dynamic response of the support structure, and comprehensive damage degree classification, providing a scientific basis for the optimization of mine support structures and the prevention and control of mine seismic disasters.

[0045] Furthermore, embodiments of this application also disclose an electronic device, Figure 3This is a structural diagram of an electronic device 20 according to an exemplary embodiment. The content of the diagram should not be construed as limiting the scope of this application.

[0046] Figure 3 This is a schematic diagram of the structure of an electronic device 20 provided in an embodiment of this application. Specifically, the electronic device 20 may include: at least one processor 21, at least one memory 22, a power supply 23, a communication interface 24, an input / output interface 25, and a communication bus 26. The memory 22 stores a computer program, which is loaded and executed by the processor 21 to implement the relevant steps in the method for evaluating the intensity of seismic events and the performance of support structures disclosed in any of the foregoing embodiments. Alternatively, the electronic device 20 in this embodiment may specifically be a computer.

[0047] In this embodiment, the power supply 23 is used to provide operating voltage for each hardware device on the electronic device 20; the communication interface 24 can create a data transmission channel between the electronic device 20 and external devices, and the communication protocol it follows can be any communication protocol applicable to the technical solution of this application, and is not specifically limited here; the input / output interface 25 is used to acquire external input data or output data to the outside world, and its specific interface type can be selected according to specific application needs, and is not specifically limited here.

[0048] In addition, the memory 22, as a carrier for resource storage, can be a read-only memory, random access memory, disk or optical disk, etc. The resources stored thereon can include operating system 221, computer program 222, etc., and the storage method can be temporary storage or permanent storage.

[0049] The operating system 221 is used to manage and control the various hardware devices on the electronic device 20 and the computer program 222, which may be Windows Server, Netware, Unix, Linux, etc. In addition to including a computer program capable of performing the method for evaluating seismic intensity and support structure performance disclosed in any of the foregoing embodiments, the computer program 222 may further include computer programs capable of performing other specific tasks.

[0050] Furthermore, this application also discloses a computer-readable storage medium for storing a computer program; wherein, when the computer program is executed by a processor, it implements the aforementioned method for evaluating the intensity of seismic events and the performance of support structures. Specific steps of this method can be found in the corresponding content disclosed in the foregoing embodiments, and will not be repeated here.

[0051] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on its differences from other embodiments. Similar or identical parts between embodiments can be referred to interchangeably. For the apparatus disclosed in the embodiments, since it corresponds to the method disclosed in the embodiments, the description is relatively simple; relevant parts can be referred to in the method section.

[0052] Those skilled in the art will further recognize that the units and algorithm steps of the various examples described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of both. To clearly illustrate the interchangeability of hardware and software, the components and steps of the various examples have been generally described in terms of functionality in the foregoing description. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of this application.

[0053] The steps of the methods or algorithms described in conjunction with the embodiments disclosed herein can be implemented directly by hardware, a software module executed by a processor, or a combination of both. The software module can be located in random access memory (RAM), main memory, read-only memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, hard disk, removable disk, CD-ROM, or any other form of storage medium known in the art.

[0054] Finally, it should be noted that in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Without further limitations, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes said element.

[0055] The technical solutions provided in this application have been described in detail above. Specific examples have been used to illustrate the principles and implementation methods of this application. The descriptions of the above embodiments are only for the purpose of helping to understand the methods and core ideas of this application. At the same time, for those skilled in the art, there will be changes in the specific implementation methods and application scope based on the ideas of this application. Therefore, the content of this specification should not be construed as a limitation of this application.

Claims

1. A method for evaluating the intensity of mine earthquakes and the performance of support structures, characterized in that, include: Collect geological survey data and support structure design parameters of the coal mine area, and determine the peak vibration velocity of the particles based on the geological survey data of the coal mine area and the pre-constructed peak vibration velocity determination model. Based on the pre-constructed mapping relationship, the coal mine seismic intensity level corresponding to the peak vibration velocity of the particle is determined. Based on the coal mine seismic intensity level, the load amplification factor is determined. The mine seismic load vector in the initial support-surrounding rock coupled dynamic model constructed based on the support structure design parameters is corrected by the load amplification factor to obtain the corrected mine seismic load vector. The stiffness matrix in the initial support-surrounding rock coupled dynamic model is modified so as to determine the target support-surrounding rock coupled dynamic model based on the modified stiffness matrix; The structural response information of the target location is determined based on the corrected mine vibration load vector and the target support-surrounding rock coupled dynamic model. The target area includes anchor bolts / cables, shotcrete layer, steel support, and surrounding rock; the structural response information includes any one or a combination of stress, strain, plastic deformation, crack width, damage factor, buckling degree, residual deformation, plastic zone range, and displacement. The comprehensive damage degree of the support structure is determined by the structural response information of the target location and the seismic intensity level of the coal mine, and the performance safety level of the support structure is determined based on the comprehensive damage degree. Before determining the peak particle vibration velocity based on the geological exploration data of the coal mine area and the pre-constructed peak particle vibration velocity determination model, the method further includes: Screening for mine tremor events that meet preset conditions that occur in coal mining areas within a target time period; Based on the parameter information corresponding to the mining earthquake event, a mathematical model is fitted using the least squares method. The accuracy of the mathematical model is verified using the leave-one-out cross-validation method. Based on the corresponding verification results, the peak vibration velocity of the mass point is determined to confirm the model. Accordingly, determining the peak particle vibration velocity based on the geological exploration data of the coal mine area and the pre-constructed peak particle vibration velocity determination model includes: Based on the geological exploration data of the coal mine area and the pre-constructed peak particle vibration velocity determination model, the peak particle vibration velocity is determined; the peak particle vibration velocity determination model is as follows: ； Where Y is the peak velocity of the particle; M is the magnitude of the mine earthquake; r is the horizontal distance from the monitoring point to the source; h0 is the equivalent source depth of the mine earthquake; S is the site category coefficient; G is the geological structure influence coefficient; a, b, c, d and e are undetermined coefficients for model fitting; ε is the random error term; The stiffness matrix in the modified initial support-surrounding rock coupled dynamic model includes: The comprehensive geological condition coefficient is determined based on the weighted average of the rock strength coefficient and the structural complexity coefficient. Determine the second product between the comprehensive geological condition coefficient and the fourth value; Determine the target sum between 1 and the second product; The product of the target and the stiffness matrix in the initial support-surrounding rock coupled dynamic model is determined as the corrected stiffness matrix.

2. The method for evaluating the intensity of mine seismic events and the performance of support structures according to claim 1, characterized in that, The collected geological survey data of the coal mine area and the design parameters of the support structure include: The geological survey data and support structure design parameters of the coal mine area are collected in real time by an array of vibration sensors deployed in underground roadways, mining areas and on the surface. The geological survey data of the coal mine area includes the horizontal distance from the monitoring point to the seismic source and the equivalent source depth of the mine earthquake. The support structure design parameters include the specifications of the anchor bolts / cables, the strength of the sprayed layer, the model of the steel support, the support density and the anchoring force.

3. The method for evaluating the intensity of mine seismic events and the performance of support structures according to claim 1, characterized in that, Before determining the coal mine seismic intensity level corresponding to the peak particle vibration velocity based on the pre-constructed mapping relationship, the process also includes: Based on enterprise standards for monitoring and evaluating seismic intensity in coal mines, types of underground coal roadway support, damage characteristics, and measured parameters from instruments, a mapping relationship between seismic intensity levels in coal mines and peak particle vibration velocities is established.

4. The method for evaluating the intensity of mine earthquakes and the performance of support structures according to claim 1, characterized in that, Before correcting the mine vibration load vector in the initial support-surrounding rock coupled dynamic model constructed based on the support structure design parameters using the load amplification factor, the process also includes: An initial support-surrounding rock coupled dynamic model is constructed based on the support structure design parameters; the initial support-surrounding rock coupled dynamic model is as follows: ； Among them, [M] s For the support-surrounding rock coupling consistent quality matrix; [C] s Let K be the Rayleigh damping matrix; s (u) is the stiffness matrix; F(t) is the mine vibration load vector; , , These are the nodal acceleration, velocity, and displacement vectors, respectively. Accordingly, the load amplification factor is determined based on the seismic intensity level of the coal mine, and the seismic load vector in the initial support-surrounding rock coupled dynamic model constructed based on the support structure design parameters is corrected using the load amplification factor to obtain the corrected seismic load vector, including: Determine the first product between the earthquake intensity level of the coal mine and the first value; The sum of the first product and the second value is determined as the load amplification factor; The product of the load amplification factor and the third value is determined as the corrected mine vibration load vector.

5. The method for evaluating the intensity of mine seismic events and the performance of support structures according to claim 1, characterized in that, The determination of the comprehensive damage degree of the support structure based on the structural response information of the target location and the seismic intensity level of the coal mine includes: Determine the comprehensive geological condition coefficient and the support structure characteristic coefficient based on the structural response information of the target location; The damage degree of the target location is determined using a damage degree determination formula based on the comprehensive geological condition coefficient, the support structure characteristic coefficient, and the seismic intensity level of the coal mine; the damage degree determination formula is as follows: ； Wherein, D represents the degree of damage; I represents the intensity level of the coal mine seismic event; G represents the comprehensive geological condition coefficient; S represents the characteristic coefficient of the support structure; and P represents the effective duration of the seismic event. The weighted average of the damage levels at each target location is determined as the overall damage level of the support structure.

6. The method for evaluating the intensity of mine seismic events and the performance of support structures according to any one of claims 1 to 5, characterized in that, Determining the performance safety level of the support structure based on the comprehensive damage degree includes: The performance safety level of the support structure is determined based on the comparison between the comprehensive damage degree and the preset threshold range; wherein, different preset threshold ranges correspond to different performance safety levels.

7. A device for evaluating the intensity of mine earthquakes and the performance of support structures, characterized in that, include: The peak vibration velocity determination module is used to collect geological exploration data of the coal mine area and design parameters of the support structure, and determine the peak vibration velocity of the particles based on the geological exploration data of the coal mine area and the pre-constructed peak vibration velocity determination model. The correction module is used to determine the coal mine seismic intensity level corresponding to the peak vibration velocity of the mass point based on the pre-constructed mapping relationship, determine the load amplification factor based on the coal mine seismic intensity level, and correct the mine vibration load vector in the initial support-surrounding rock coupled dynamic model based on the support structure design parameters through the load amplification factor to obtain the corrected mine vibration load vector. The model determination module is used to correct the stiffness matrix in the initial support-surrounding rock coupled dynamic model, so as to determine the target support-surrounding rock coupled dynamic model based on the corrected stiffness matrix. The structural response information determination module is used to determine the structural response information of the target location based on the corrected mine vibration load vector and the target support-surrounding rock coupled dynamic model. The target area includes anchor bolts / cables, shotcrete layer, steel support, and surrounding rock; the structural response information includes any one or a combination of stress, strain, plastic deformation, crack width, damage factor, buckling degree, residual deformation, plastic zone range, and displacement. The performance safety level determination module is used to determine the comprehensive damage degree of the support structure by using the structural response information of the target location and the intensity level of the coal mine seismic event, and to determine the performance safety level of the support structure based on the comprehensive damage degree. The device is also used to screen coal mine seismic events that meet preset conditions and occur in the coal mine area within a target time period; to fit a mathematical model by least squares method based on the parameter information corresponding to the seismic events; to verify the accuracy of the mathematical model by leave-one-out cross-validation method; and to determine the peak vibration velocity of the mass point and determine the model based on the corresponding verification results. The peak particle vibration velocity determination module is used to determine the peak particle vibration velocity based on the geological exploration data of the coal mine area and a pre-constructed peak particle vibration velocity determination model; the peak particle vibration velocity determination model is as follows: ； Where Y is the peak velocity of the particle; M is the magnitude of the mine earthquake; r is the horizontal distance from the monitoring point to the source; h0 is the equivalent source depth of the mine earthquake; S is the site category coefficient; G is the geological structure influence coefficient; a, b, c, d and e are undetermined coefficients for model fitting; ε is the random error term; The model determination module is used to determine the comprehensive geological condition coefficient based on the weighted average result of the rock strength coefficient and the structural complexity coefficient. Determine the second product between the comprehensive geological condition coefficient and the fourth value; determine the target sum between 1 and the second product; and determine the corrected stiffness matrix by multiplying the target sum with the stiffness matrix in the initial support-surrounding rock coupled dynamic model.

8. An electronic device, characterized in that, include: Memory, used to store computer programs; A processor for executing a computer program to implement the steps of the method for evaluating the intensity of seismic events and the performance of support structures as described in any one of claims 1 to 6.

9. A computer-readable storage medium, characterized in that, A computer-readable storage medium stores a computer program that, when executed by a processor, implements the steps of the method for evaluating the intensity of seismic events and the performance of support structures as described in any one of claims 1 to 6.

Images