Method and device for evaluating mine earthquake intensity of coal mine tunnel

Abstract

The invention provides a method and a device for evaluating mine earthquake intensity of a coal mine tunnel, and relates to the technical field of coal mine safety. The method comprises the following steps: acquiring mine seismic waveform data of a coal mine tunnel area; and determining the mine earthquake intensity of the coal mine tunnel area according to the peak particle velocity PPV value of the mine earthquake waveform data and a predetermined mapping relation between the PPV value and the mine earthquake intensity. The mapping relation between the PPV value and the mine earthquake intensity is determined by performing numerical simulation on various mine earthquake waveform sample data and a numerical model. The mine earthquake intensity can be accurately predicted according to the mine earthquake waveform data, and the obtained mine earthquake intensity can accurately reflect the influence and damage of the mine earthquake waveform data on the roadway, so that a worker can evaluate potential threats of the roadway structure, the effectiveness of mine safety management is remarkably improved, and the safety of the mine is improved. Therefore, economic losses caused by mine earthquakes are effectively controlled and reduced.

Description

Technical Field

[0001] The present invention relates to the technical field of coal mine safety, and in particular to a method and device for evaluating earthquake intensity in a coal mine tunnel. Background Art

[0002] Earthquakes caused by mining are called mine tremors, which are a type of induced earthquakes. Mine tremors often occur along with mining, and their focal depths are often at or above the mining surface. Therefore, the focal depth of mine tremors is often only tens of meters to thousands of meters. Compared with natural earthquakes, they are extremely shallow-source earthquakes, with the characteristics of shallow focal points, slow attenuation, and strong destructive power. After a mine tremor occurs, it will have an impact both underground and on the ground in the mine, and the intensity of a mine tremor is an important indicator for dividing the degree of damage and danger of a mine tremor. These impacts play an important role in safe and green production.

[0003] At present, the assessment of mine earthquake intensity mainly refers to the earthquake intensity scale, that is, the impact of earthquakes on the ground to evaluate the strength of its impact. This method combines macroscopic surveys and instrument measurements to obtain intensity levels, taking into account factors such as house damage, human perception, object reactions, lifeline engineering damage, and other earthquake damage phenomena. According to the breaking of thick and hard key layers and the law of energy propagation, the concept of "vibration damage boundary" was proposed, and a preliminary assessment method for ground vibration damage induced by mine earthquakes was established. However, the mine earthquake intensity determined by the above method lacks systematicity and pertinence in assessing the impact and damage of coal mine tunnels. No reasonable basis for division has been proposed, and the impact of mine earthquakes on coal mine tunnels cannot be reflected. Summary of the invention

[0004] The present disclosure aims to solve one of the technical problems in the related art at least to some extent.

[0005] To this end, the first embodiment of the present disclosure proposes a method for evaluating the intensity of mine earthquakes in coal mine tunnels, comprising:

[0006] Obtain mine earthquake waveform data in coal mine tunnel areas;

[0007] Determine the mine earthquake intensity of the coal mine roadway area according to the peak particle velocity PPV value of the mine earthquake waveform data and the mapping relationship between the predetermined PPV value and the mine earthquake intensity;

[0008] The mapping relationship between the PPV value and the mine earthquake intensity is determined by the following method:

[0009] Acquire mechanical parameters of coal seams and surrounding rocks in the coal mine tunnel area, as well as geological data of the coal mine tunnel area;

[0010] Constructing a numerical model of the coal mine tunnel area according to the mechanical parameters of the coal seam, the mechanical parameters of the surrounding rock and the geological data;

[0011] embedding a cohesive unit into the numerical model, wherein the cohesive unit is used to describe the damage initiation and development process of the material interface in the coal mine tunnel area;

[0012] Inputting a variety of mining earthquake waveform sample data into the numerical model embedded with the cohesive unit respectively, simulating the damage evolution process of the coal mine roadway area under the various mining earthquake waveform sample data, and obtaining the roadway stress information, roadway deformation information, roadway damage morphology information and the damage information of the cohesive unit corresponding to the PPV value of each of the mining earthquake waveform sample data after simulation;

[0013] According to the tunnel stress information, tunnel deformation information, tunnel damage morphology information and the damage information of the cohesive unit corresponding to the PPV value of each of the mine earthquake waveform sample data, the mapping relationship between the PPV value and the mine earthquake intensity is determined.

[0014] In some embodiments of the present disclosure, the numerical model of the coal mine tunnel area is constructed according to the mechanical parameters of the coal seam, the mechanical parameters of the surrounding rock and the geological data, including: constructing a three-dimensional geometric model of the coal mine tunnel area based on the geological data; after meshing the three-dimensional geometric model, constructing the numerical model of the coal mine tunnel area in combination with the mechanical parameters of the coal seam and the mechanical parameters of the surrounding rock.

[0015] In some embodiments of the present disclosure, the tunnel damage morphology information includes at least one of the following: tunnel crack distribution information and broken area information.

[0016] In some embodiments of the present disclosure, the mapping relationship between the PPV value and the mine earthquake intensity is determined based on the tunnel stress information, tunnel deformation information, tunnel damage morphology information and damage information of the cohesive unit corresponding to the PPV value of each of the mine earthquake waveform sample data, including: evaluating the degree of tunnel damage corresponding to the PPV value of each of the mine earthquake waveform sample data, and obtaining the mapping relationship between the PPV value and the degree of tunnel damage; establishing a mapping relationship between the degree of tunnel damage and the mine earthquake intensity according to a mine earthquake classification standard; determining the mapping relationship between the PPV value and the mine earthquake intensity according to the mapping relationship between the PPV value and the degree of tunnel damage, and the mapping relationship between the degree of tunnel damage and the mine earthquake intensity.

[0017] The second aspect of the present disclosure provides a device for evaluating the intensity of earthquake in a coal mine tunnel, comprising:

[0018] An acquisition module, used to acquire mine earthquake waveform data in the coal mine tunnel area;

[0019] The first determination module is used to determine the mine earthquake intensity of the coal mine tunnel area according to the peak particle velocity PPV value of the mine earthquake waveform data and the mapping relationship between the predetermined PPV value and the mine earthquake intensity.

[0020] In some embodiments of the present disclosure, a second determination module is also included; the second determination module includes: an acquisition unit, used to acquire the mechanical parameters of the coal seam and the mechanical parameters of the surrounding rock in the coal mine tunnel area, as well as the geological data of the coal mine tunnel area; a first construction unit, used to construct a numerical model of the coal mine tunnel area according to the mechanical parameters of the coal seam, the mechanical parameters of the surrounding rock and the geological data; a second construction unit, used to embed a cohesive unit into the numerical model, the cohesive unit is used to describe the damage initiation and development process of the material interface in the coal mine tunnel area; a simulation unit, used to integrate multiple mine shock waves The shape sample data are respectively input into the numerical model after the cohesive unit is embedded, so as to simulate the damage evolution process of the coal mine roadway area under the multiple mine earthquake waveform sample data, and obtain the roadway stress information, roadway deformation information, roadway damage morphology information and the damage information of the cohesive unit corresponding to the PPV value of each of the mine earthquake waveform sample data after simulation; and the determination unit is used to determine the mapping relationship between the PPV value and the mine earthquake intensity according to the roadway stress information, roadway deformation information, roadway damage morphology information and the damage information of the cohesive unit corresponding to the PPV value of each of the mine earthquake waveform sample data.

[0021] In some embodiments of the present disclosure, the first construction unit is specifically used to: construct a three-dimensional geometric model of the coal mine tunnel area based on the geological data; after meshing the three-dimensional geometric model, construct a numerical model of the coal mine tunnel area in combination with the mechanical parameters of the coal seam and the mechanical parameters of the surrounding rock.

[0022] In some embodiments of the present disclosure, the tunnel damage morphology information includes at least one of the following: tunnel crack distribution information and broken area information.

[0023] In some embodiments of the present disclosure, the determination unit is specifically used to: evaluate the degree of tunnel damage corresponding to the PPV value of each type of mine earthquake waveform sample data according to the tunnel stress information, tunnel deformation information, tunnel damage morphology information and damage information of the cohesive unit corresponding to the PPV value of each type of mine earthquake waveform sample data, and obtain a mapping relationship between the PPV value and the degree of tunnel damage; establish a mapping relationship between the degree of tunnel damage and the intensity of mine earthquake according to a mine earthquake classification standard; determine the mapping relationship between the PPV value and the degree of tunnel damage, and the mapping relationship between the degree of tunnel damage and the intensity of mine earthquake according to the mapping relationship between the PPV value and the degree of tunnel damage, and the mapping relationship between the degree of tunnel damage and the intensity of mine earthquake.

[0024] A third aspect of the present application provides an electronic device, comprising: a processor, and a memory communicatively connected to the processor;

[0025] The memory stores computer-executable instructions;

[0026] The processor executes the computer-executable instructions stored in the memory to implement the method described in the first aspect.

[0027] The method for assessing the mine earthquake intensity in coal mine tunnels provided by the present invention can accurately predict the mine earthquake intensity based on the mapping relationship between the predetermined PPV value and the mine earthquake intensity and the mine earthquake waveform data of the coal mine tunnel area. The obtained mine earthquake intensity can accurately reflect the impact and damage of the mine earthquake waveform data on the tunnel, so that coal mine workers can assess the potential threat of the mine earthquake waveform data to the tunnel structure, which significantly improves the effectiveness of mine safety management and can effectively control and reduce the economic losses caused by mine earthquakes.

[0028] Additional aspects and advantages of the present disclosure will be given in part in the following description and in part will be obvious from the following description or learned through practice of the present disclosure. BRIEF DESCRIPTION OF THE DRAWINGS

[0029] The above and / or additional aspects and advantages of the present disclosure will become apparent and easily understood from the following description of the embodiments in conjunction with the accompanying drawings, in which:

[0030] Figure 1 A flow chart of a method for determining a mapping relationship between a PPV value and a mine earthquake intensity provided in an embodiment of the present disclosure;

[0031] Figure 2 A schematic diagram of a numerical model of a coal mine tunnel area provided by an embodiment of the present disclosure;

[0032] Figure 3 A schematic diagram of embedding a cohesive unit in a numerical model of a coal mine tunnel area provided by an embodiment of the present disclosure;

[0033] Figure 4 A schematic diagram of a mixed mode cohesive traction response provided by an embodiment of the present disclosure;

[0034] Figure 5 A schematic diagram of a flow chart of a method for assessing earthquake intensity in a coal mine tunnel provided by an embodiment of the present disclosure;

[0035] Figure 6 A comparison diagram of original waveform data and pre-processed mining earthquake waveform data provided by an embodiment of the present disclosure;

[0036] Figure 7 A schematic diagram of a device for assessing earthquake intensity in a coal mine tunnel provided by an embodiment of the present disclosure;

[0037] Figure 8 A schematic diagram of another device for assessing earthquake intensity in coal mine tunnels provided in an embodiment of the present disclosure. DETAILED DESCRIPTION

[0038] Embodiments of the present disclosure are described in detail below, examples of which are shown in the accompanying drawings, wherein the same or similar reference numerals throughout represent the same or similar elements or elements having the same or similar functions. The embodiments described below with reference to the accompanying drawings are exemplary and are intended to be used to explain the present disclosure, and should not be construed as limiting the present disclosure.

[0039] Specifically, the following describes the method and device for evaluating the seismic intensity of a coal mine tunnel according to an embodiment of the present disclosure with reference to the accompanying drawings. In order to better understand the method for evaluating the seismic intensity of a coal mine tunnel according to an embodiment of the present disclosure, the following first describes a method for determining the mapping relationship between the PPV value and the seismic intensity. Figure 1 The following is a flow chart of a method for determining a mapping relationship between a PPV value and a mine earthquake intensity provided by an embodiment of the present disclosure. Figure 1 As shown, the method for determining the mapping relationship between the PPV value and the mine earthquake intensity may include the following steps:

[0040] Step 101, obtaining mechanical parameters of the coal seam and the surrounding rock in the coal mine tunnel area, as well as geological data of the coal mine tunnel area.

[0041] In some embodiments of the present disclosure, the geological data of the coal mine tunnel area may include data such as stratum lithology, stratigraphic sequence, thickness, fault distribution, groundwater conditions, and the actual layout and dimensions of the mine tunnels, and the mechanical parameters of the coal seam and surrounding rock may include parameters such as density, elastic modulus, Poisson's ratio, compressive strength, etc. obtained through geological exploration and field measurements.

[0042] Step 102, constructing a numerical model of the coal mine tunnel area according to the mechanical parameters of the coal seam, the mechanical parameters of the surrounding rock and the geological data.

[0043] It should be noted that the numerical model of the coal mine tunnel area can reflect the topographic and geological characteristics of the coal mine tunnel area through geological data, such as the length, width, height, curves, intersections, etc. of the tunnel. The numerical model sets boundary conditions through the actual action relationship in the coal mine tunnel area to reflect the mechanical properties of the material in the coal mine tunnel area.

[0044] In some embodiments of the present disclosure, a three-dimensional geometric model of a coal mine tunnel area can be constructed based on geological data using dedicated geological modeling software, and the geometric model can be meshed using software such as ANSYS, Abaqus, and HyperMesh to ensure a balance between mesh size and computational efficiency. After meshing the three-dimensional geometric model, boundary conditions are set in combination with actual interaction relationships in the coal mine tunnel area to construct a numerical model of the coal mine tunnel area. Figure 2 A schematic diagram of a numerical model of a coal mine tunnel area provided in an embodiment of the present disclosure.

[0045] Step 103, embedding the cohesive unit into the numerical model, the cohesive unit is used to describe the damage initiation and development process of the material interface in the coal mine tunnel area.

[0046] It should be noted that the cohesive unit is specially designed to simulate the bonding behavior within a material or between different materials. Figure 3 A schematic diagram of embedding a cohesive unit in a numerical model of a coal mine tunnel area provided in an embodiment of the present disclosure. In some embodiments of the present disclosure, the cohesive unit can be automatically generated by writing code. Each cohesive unit has a specific force-displacement relationship, namely, a traction-separation constitutive model, which is used to describe the damage initiation and development process of the material interface.

[0047] Among them, the initial linear elastic behavior of the traction separation law can be expressed as follows:

[0048]

[0049] where t is the nominal traction stress vector, which in three-dimensional problems consists of three components (t n ,t s and t t ) composition; n , ε s and ε t are the three components of nominal strain; δ n , δ s and δt is the corresponding spacing; T0 is the original thickness of the viscous unit; the matrix containing parameters such as Enn is the elastic matrix, which provides the complete coupling behavior between all components of the traction vector and the separation vector. A schematic diagram of the commonly used mixed mode cohesive traction response can be found in Figure 4 .like Figure 4 As shown, t no ,t so ,t to and t mo are the relative normal direction, relative shear force, relative shear force and effective traction at the beginning of the crack; δ no , δ so , δ to and δ mo are the relative normal, relative shear, relative shear force and effective displacement at the beginning of damage; δ nf , δ sf , δ tf and δ mf They are pure normal displacement, pure shear displacement, pure shear displacement and effective displacement respectively. no , k so and k to They are initial tensile stiffness, initial shear stiffness and initial shear stiffness respectively.

[0050] In some embodiments of the present disclosure, material properties can be set for the cohesive unit in the numerical model according to the cohesive unit material parameter library of the rock material designed and trained in advance, including parameters such as the modulus of the normal (tensile) and tangential (shear), the bond strength and the fracture energy. Among them, the cohesive unit material parameter library of the rock material may include the property information of a variety of rock materials, and according to the actual rock materials in the coal mine roadway area, the matching material properties are determined in the cohesive unit material parameter library of the rock material.

[0051] After the cohesive unit used to describe the damage initiation and development process of the material interface in the coal mine tunnel area is embedded in the numerical model, the mine shock waveform can be input into the numerical model after the cohesive unit is embedded, and the damage evolution process of the tunnel under different mine shock waveforms can be simulated by explicit iterative solution. When the stress in the cohesive unit reaches the strength limit of the material, damage occurs until the failure limit is reached and the unit fails. The damage evolution criterion can refer to formula (2)-formula (6):

[0052] Due to damage accumulation, the changes in the normal stress and shear stress components of the traction-separation model can be referred to as follows:

[0053]

[0054]

[0055] The tensile stiffness and shear stiffness can be expressed as:

[0056] k n =(1-D)k no

[0057] k s =(1-D)k so

[0058] k t =(1-D)k to (4)

[0059] In order to describe the damage evolution of a crack under the combination of normal deformation and shear deformation on the crack surface, the effective displacement is introduced as follows:

[0060]

[0061] Among them, δ m is an equivalent pure displacement, <> is a Macaulay bracket, when <δ n > is greater than zero, <δ n > equal to δ n , otherwise <δ n > is equal to zero.

[0062] The damage variable D can be evolved according to the following expression:

[0063]

[0064] Among them, δ mm It represents the maximum pure displacement during loading. The damage variable D, or stiffness reduction rate, is a parameter used to describe the degree of damage.

[0065] Step 104, input various mining earthquake waveform sample data into the numerical model embedded with the cohesive unit respectively, simulate the damage evolution process of the coal mine tunnel area under various mining earthquake waveform sample data, and obtain the tunnel stress information, tunnel deformation information, tunnel damage morphology information and cohesive unit damage information corresponding to the PPV value of each mining earthquake waveform sample data after simulation.

[0066] It should be noted that the multiple mining earthquake waveform sample data are multiple different mining earthquake waveform sample data, which can be obtained by preprocessing the original collected mining earthquake waveform. The multiple mining earthquake waveform sample data are respectively input into the numerical model embedded with the cohesive unit to perform numerical simulation to simulate the damage evolution process of the coal mine tunnel area under the multiple mining earthquake waveform sample data. The stress distribution map, deformation cloud map and damage cloud map of the numerical model under the action of each mining earthquake waveform sample data (i.e., different PPV values) can be output through the numerical simulation software. The damage information of the cohesive unit is obtained through the damage cloud map. The damage information of the cohesive unit may include the stiffness reduction rate D and the damage state value (both values ​​are: 1 represents complete failure, 0 represents no damage at all, and between 0-1 means damage but not complete failure, corresponding to Figure 4 The damage curve is 0 is the part before the vertex of the curve, 1 represents the end point where the curve finally intersects with the x-axis, and 0-1 is the descending curve). The stress distribution diagram can be used to analyze the stress information at various locations inside the tunnel, such as the maximum principal stress and shear stress of the surrounding rock around the tunnel. The deformation cloud map can be used to determine the tunnel deformation information of the numerical model, especially the deformation features such as the sinking of the tunnel roof, the displacement of the two side walls, and the uplift of the tunnel bottom. The tunnel damage morphology information can be obtained through any one of the stress distribution diagram, deformation cloud map, and damage cloud map. In some embodiments of the present disclosure, the tunnel damage morphology information may include at least one of the following: tunnel crack distribution information and broken area information.

[0067] Step 105, determining the mapping relationship between the PPV value and the mine earthquake intensity according to the tunnel stress information, tunnel deformation information, tunnel damage morphology information and cohesive unit damage information corresponding to the PPV value of each mine earthquake waveform sample data.

[0068] After numerical simulation of the numerical model based on the mine earthquake waveform sample data, the PPV value of each mine earthquake waveform sample data corresponds to a set of roadway stress information, roadway deformation information, roadway damage morphology information and cohesive unit damage information. The damage condition of the coal mine roadway is evaluated based on each set of roadway stress information, roadway deformation information, roadway damage morphology information and cohesive unit damage information. The damage condition of the coal mine roadway is matched with the damage condition in the existing mine earthquake intensity table to determine the mine earthquake intensity corresponding to the PPV value of each mine earthquake waveform sample data.

[0069] In some embodiments of the present disclosure, the degree of tunnel damage corresponding to the PPV value of each mine-seismic waveform sample data can be evaluated based on the tunnel stress information, tunnel deformation information, tunnel damage morphology information and cohesive unit damage information corresponding to the PPV value of each mine-seismic waveform sample data, and the mapping relationship between the PPV value and the tunnel damage degree can be obtained.

[0070] Optionally, the overall damage degree of the tunnel can be evaluated based on the tunnel stress information, tunnel deformation information, tunnel damage morphology information and damage information of the cohesive unit, such as safety factor calculation, probabilistic damage model, etc. The tunnel damage degree corresponding to different PPV values ​​is counted, and a mapping relationship between the PPV value and the tunnel damage degree is established. Then, based on the existing earthquake intensity scale, mine earthquake intensity scale and mine earthquake classification and grading standards, the description of earthquake damage to houses, human feelings, object reactions and instrument-measured earthquake intensity is integrated to establish a mapping relationship between the tunnel damage degree and the mine earthquake intensity. Based on the mapping relationship between the PPV value and the tunnel damage degree, and the mapping relationship between the tunnel damage degree and the mine earthquake intensity, the mapping relationship between the PPV value and the mine earthquake intensity is determined.

[0071] By implementing the embodiments of the present disclosure, a numerical model of the coal mine tunnel area constructed based on geological data and mechanical parameters of materials in the tunnel is more in line with actual on-site conditions. By numerically simulating the numerical model, the damage to the coal mine tunnel under the influence of different mine earthquake waveforms is obtained, and then a mapping relationship between the PPV value of the mine earthquake waveform and the mine earthquake intensity is established, the scope and degree of the impact of the mine earthquake on the coal mine tunnel is accurately restored and predicted, and a theoretical basis and technical support are provided for the evaluation of the mine earthquake intensity. In subsequent applications, the mine earthquake intensity can be accurately evaluated based on the PPV value of the mine earthquake waveform, and the potential threat of mine earthquakes to the tunnel structure can be predicted, which significantly improves the effectiveness of mine safety management and effectively controls and reduces the economic losses caused by mine earthquakes.

[0072] Figure 5 The following is a flow chart of a method for evaluating the intensity of earthquake in a coal mine tunnel provided by an embodiment of the present disclosure. Figure 5 As shown, the method for evaluating the seismic intensity of a coal mine tunnel may include the following steps:

[0073] Step 501, obtaining mine earthquake waveform data in the coal mine tunnel area.

[0074] In some embodiments of the present disclosure, high-sensitivity seismographs and microseismic monitoring systems can be used to capture real-time vibration signals in coal mines. These detection devices can cover major mining surfaces and potential high-risk areas to ensure comprehensive capture of mine earthquake activities. Seismic waveform signals are collected by the microseismic monitoring system and transmitted to the central processing unit. During the data acquisition process, time synchronization and data integrity must be guaranteed to facilitate subsequent accurate analysis. The collected raw waveform data can be pre-processed, such as denoising, signal amplification, and data compression, and organized and optimized into waveform data files that can be input into numerical simulations for analysis to obtain mine earthquake waveform data (i.e., waveform files corresponding to different PPV values).

[0075] Figure 6A comparison diagram of the original waveform data and the pre-processed mining earthquake waveform data provided by the embodiment of the present disclosure, wherein Figure 6 (Left) is the original waveform data, Figure 6 (Right) The mine earthquake waveform data is obtained after preprocessing the original waveform data.

[0076] Step 502, determining the mine earthquake intensity in the coal mine tunnel area according to the peak particle velocity PPV value of the mine earthquake waveform data and the mapping relationship between the predetermined PPV value and the mine earthquake intensity.

[0077] Among them, the determination process of the mapping relationship between the PPV value and the predetermined PPV value and the mine earthquake intensity has been Figure 1 The above is described in the embodiment shown and will not be described in detail here.

[0078] By implementing the embodiments of the present disclosure, the mine earthquake intensity is accurately predicted based on the mapping relationship between the predetermined PPV value and the mine earthquake intensity and the mine earthquake waveform data of the coal mine tunnel area. The obtained mine earthquake intensity can accurately reflect the impact and damage of the mine earthquake waveform data on the tunnel, so that coal mine workers can evaluate the potential threat of the mine earthquake waveform data to the tunnel structure, which significantly improves the effectiveness of mine safety management and effectively controls and reduces the economic losses caused by mine earthquakes.

[0079] Figure 7 A schematic diagram of a device for evaluating the intensity of earthquakes in a coal mine tunnel provided by an embodiment of the present disclosure. Figure 7 As shown, the device for evaluating the seismic intensity of a coal mine tunnel includes: an acquisition module 701 and a first determination module 702 .

[0080] Among them, the acquisition module 701 is used to obtain the mine earthquake waveform data in the coal mine tunnel area.

[0081] The first determination module 702 is used to determine the mine earthquake intensity in the coal mine tunnel area according to the peak particle velocity PPV value of the mine earthquake waveform data and the mapping relationship between the predetermined PPV value and the mine earthquake intensity.

[0082] Optionally, in some embodiments of the present disclosure, Figure 8 As shown, the device for assessing earthquake intensity in coal mine tunnels further includes a second determination module 803. The second determination module 803 includes: an acquisition unit 804, a first construction unit 805, a second construction unit 806, a simulation unit 807 and a determination unit 808.

[0083] An acquisition unit 804 is used to acquire mechanical parameters of the coal seam and the surrounding rock in the coal mine tunnel area, as well as geological data of the coal mine tunnel area;

[0084] The first construction unit 805 is used to construct a numerical model of the coal mine tunnel area according to the mechanical parameters of the coal seam, the mechanical parameters of the surrounding rock and the geological data;

[0085] A second construction unit 806 is used to embed a cohesive unit into the numerical model, where the cohesive unit is used to describe the damage initiation and development process of the material interface in the coal mine tunnel area;

[0086] The simulation unit 807 is used to input various mining earthquake waveform sample data into the numerical model embedded with the cohesive unit, simulate the damage evolution process of the coal mine roadway area under various mining earthquake waveform sample data, and obtain the roadway stress information, roadway deformation information, roadway damage morphology information and cohesive unit damage information corresponding to the PPV value of each mining earthquake waveform sample data after simulation;

[0087] The determination unit 808 is used to determine the mapping relationship between the PPV value and the mine earthquake intensity according to the tunnel stress information, tunnel deformation information, tunnel damage morphology information and cohesive unit damage information corresponding to the PPV value of each mine earthquake waveform sample data.

[0088] In some embodiments of the present disclosure, the first construction unit 805 is specifically used to: construct a three-dimensional geometric model of the coal mine tunnel area based on geological data; after meshing the three-dimensional geometric model, construct a numerical model of the coal mine tunnel area by combining the mechanical parameters of the coal seam and the mechanical parameters of the surrounding rock.

[0089] In some embodiments of the present disclosure, the tunnel damage morphology information includes at least one of the following: tunnel crack distribution information and broken area information.

[0090] In some embodiments of the present disclosure, the determination unit 808 is specifically used to: evaluate the degree of tunnel damage corresponding to the PPV value of each type of mine earthquake waveform sample data according to the tunnel stress information, tunnel deformation information, tunnel damage morphology information and damage information of the cohesive unit corresponding to the PPV value of each type of mine earthquake waveform sample data, and obtain the mapping relationship between the PPV value and the degree of tunnel damage; establish a mapping relationship between the degree of tunnel damage and the intensity of mine earthquake according to the mine earthquake classification standard; determine the mapping relationship between the PPV value and the degree of tunnel damage, and the mapping relationship between the degree of tunnel damage and the intensity of mine earthquake according to the mapping relationship between the PPV value and the degree of tunnel damage, and the mapping relationship between the degree of tunnel damage and the intensity of mine earthquake.

[0091] in, Figure 8 801-802 and Figure 7 701-702 have the same function and structure.

[0092] Regarding the device in the above embodiment, the specific manner in which each module performs operations has been described in detail in the embodiment of the method, and will not be elaborated here.

[0093] In order to implement the above embodiments, the present disclosure also proposes an electronic device, including: a processor, and a memory communicatively connected to the processor; the memory stores computer-executable instructions; the processor executes the computer-executable instructions stored in the memory to implement the method provided by the above embodiments.

[0094] In order to implement the above embodiments, the present disclosure further proposes a computer-readable storage medium, in which computer-executable instructions are stored. When the computer-executable instructions are executed by a processor, they are used to implement the methods provided by the above embodiments.

[0095] In order to implement the above embodiments, the present disclosure also proposes a computer program product, including a computer program, which implements the methods provided by the above embodiments when executed by a processor.

[0096] In the description of the aforementioned embodiments, the description with reference to the terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" etc. means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present disclosure. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any one or more embodiments or examples in a suitable manner. In addition, those skilled in the art may combine and combine the different embodiments or examples described in this specification and the features of the different embodiments or examples, unless they are contradictory.

[0097] In addition, the terms "first" and "second" are used for descriptive purposes only and should not be understood as indicating or implying relative importance or implicitly indicating the number of the indicated technical features. Thus, a feature defined as "first" or "second" may explicitly or implicitly include at least one of the features. In the description of the present disclosure, "plurality" means at least two, such as two, three, etc., unless otherwise clearly and specifically defined.

[0098] Any process or method description in a flowchart or otherwise described herein may be understood to represent a module, segment or portion of code that includes one or more executable instructions for implementing the steps of a custom logical function or process, and the scope of the preferred embodiments of the present disclosure includes additional implementations in which functions may not be performed in the order shown or discussed, including performing functions in a substantially simultaneous manner or in reverse order depending on the functions involved, which should be understood by technicians in the technical field to which the embodiments of the present disclosure belong.

[0099] The logic and / or steps represented in the flowchart or otherwise described herein, for example, can be considered as an ordered list of executable instructions for implementing logical functions, and can be embodied in any computer-readable medium for use by an instruction execution system, device or apparatus (such as a computer-based system, a system including a processor, or other system that can fetch instructions from an instruction execution system, device or apparatus and execute the instructions), or in combination with these instruction execution systems, devices or apparatuses. For the purpose of this specification, "computer-readable medium" can be any device that can contain, store, communicate, propagate or transmit a program for use by an instruction execution system, device or apparatus, or in combination with these instruction execution systems, devices or apparatuses. More specific examples of computer-readable media (a non-exhaustive list) include the following: an electrical connection with one or more wires (electronic device), a portable computer disk box (magnetic device), a random access memory (RAM), a read-only memory (ROM), an erasable and programmable read-only memory (EPROM or flash memory), a fiber optic device, and a portable compact disk read-only memory (CDROM). In addition, the computer-readable medium may even be paper or other suitable medium on which the program is printed, since the program may be obtained electronically, for example, by optically scanning the paper or other medium and then editing, interpreting or processing in other suitable ways if necessary, and then stored in a computer memory.

[0100] It should be understood that the various parts of the present disclosure can be implemented in hardware, software, firmware or a combination thereof. In the above-mentioned embodiments, multiple steps or methods can be implemented in software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented in hardware, as in another embodiment, any one of the following technologies known in the art or a combination thereof can be used to implement: a discrete logic circuit having a logic gate circuit for implementing a logic function for a data signal, a dedicated integrated circuit having a suitable combination of logic gate circuits, a programmable gate array (PGA), a field programmable gate array (FPGA), etc.

[0101] A person skilled in the art may understand that all or part of the steps in the method for implementing the above-mentioned embodiment may be completed by instructing related hardware through a program, and the program may be stored in a computer-readable storage medium, which, when executed, includes one or a combination of the steps of the method embodiment.

Claims

1. A method for evaluating the intensity of earthquake in a coal mine tunnel, characterized in that: The following steps are involved: Obtain mine earthquake waveform data in coal mine tunnel areas; Determine the mine earthquake intensity of the coal mine roadway area according to the peak particle velocity PPV value of the mine earthquake waveform data and the mapping relationship between the predetermined PPV value and the mine earthquake intensity; The mapping relationship between the PPV value and the mine earthquake intensity is determined by the following method: Acquire mechanical parameters of coal seams and surrounding rocks in the coal mine tunnel area, as well as geological data of the coal mine tunnel area; Constructing a numerical model of the coal mine tunnel area according to the mechanical parameters of the coal seam, the mechanical parameters of the surrounding rock and the geological data; embedding a cohesive unit into the numerical model, wherein the cohesive unit is used to describe the damage initiation and development process of the material interface in the coal mine tunnel area; Inputting a variety of mining earthquake waveform sample data into the numerical model embedded with the cohesive unit respectively, simulating the damage evolution process of the coal mine roadway area under the various mining earthquake waveform sample data, and obtaining the roadway stress information, roadway deformation information, roadway damage morphology information and the damage information of the cohesive unit corresponding to the PPV value of each of the mining earthquake waveform sample data after simulation; According to the tunnel stress information, tunnel deformation information, tunnel damage morphology information and the damage information of the cohesive unit corresponding to the PPV value of each of the mine earthquake waveform sample data, the mapping relationship between the PPV value and the mine earthquake intensity is determined.

2. The method according to claim 1, characterized in that The step of constructing a numerical model of the coal mine tunnel area according to the mechanical parameters of the coal seam, the mechanical parameters of the surrounding rock and the geological data includes: Constructing a three-dimensional geometric model of the coal mine tunnel area based on the geological data; After meshing the three-dimensional geometric model, a numerical model of the coal mine tunnel area is constructed by combining the mechanical parameters of the coal seam and the mechanical parameters of the surrounding rock.

3. The method according to claim 1, characterized in that The tunnel damage morphology information includes at least one of the following: tunnel crack distribution information and broken area information.

4. The method according to claim 1, characterized in that Determining the mapping relationship between the PPV value and the mine earthquake intensity according to the tunnel stress information, tunnel deformation information, tunnel damage morphology information and the damage information of the cohesive unit corresponding to the PPV value of each mine earthquake waveform sample data includes: According to the tunnel stress information, tunnel deformation information, tunnel damage morphology information and the damage information of the cohesive unit corresponding to the PPV value of each of the mine earthquake waveform sample data, the tunnel damage degree corresponding to the PPV value of each of the mine earthquake waveform sample data is evaluated to obtain the mapping relationship between the PPV value and the tunnel damage degree; Establish a mapping relationship between the degree of tunnel damage and the intensity of mine earthquake according to the mine earthquake classification standard; The mapping relationship between the PPV value and the mine earthquake intensity is determined based on the mapping relationship between the PPV value and the degree of tunnel damage, and the mapping relationship between the degree of tunnel damage and the mine earthquake intensity.

5. A device for evaluating the intensity of earthquake in a coal mine tunnel, characterized in that: include: An acquisition module, used to acquire mine earthquake waveform data in the coal mine tunnel area; The first determination module is used to determine the mine earthquake intensity of the coal mine tunnel area according to the peak particle velocity PPV value of the mine earthquake waveform data and the mapping relationship between the predetermined PPV value and the mine earthquake intensity.

6. The device according to claim 5, characterized in that Also includes a second determination module; the second determination module includes: An acquisition unit, used to acquire mechanical parameters of the coal seam and the surrounding rock in the coal mine tunnel area, as well as geological data of the coal mine tunnel area; A first construction unit is used to construct a numerical model of the coal mine tunnel area according to the mechanical parameters of the coal seam, the mechanical parameters of the surrounding rock and the geological data; A second construction unit is used to embed a cohesive unit into the numerical model, wherein the cohesive unit is used to describe the damage initiation and development process of the material interface in the coal mine tunnel area; A simulation unit, used to input a plurality of mining earthquake waveform sample data into a numerical model embedded with a cohesive unit, simulate the damage evolution process of the coal mine roadway area under the plurality of mining earthquake waveform sample data, and obtain the roadway stress information, roadway deformation information, roadway damage morphology information and the damage information of the cohesive unit corresponding to the PPV value of each mining earthquake waveform sample data after simulation; A determination unit is used to determine the mapping relationship between the PPV value and the mine earthquake intensity according to the tunnel stress information, tunnel deformation information, tunnel damage morphology information and damage information of the cohesive unit corresponding to the PPV value of each of the mine earthquake waveform sample data.

7. The device according to claim 6, characterized in that The first building block is specifically used for: Constructing a three-dimensional geometric model of the coal mine tunnel area based on the geological data; After meshing the three-dimensional geometric model, a numerical model of the coal mine tunnel area is constructed by combining the mechanical parameters of the coal seam and the mechanical parameters of the surrounding rock.

8. The device according to claim 6, characterized in that The tunnel damage morphology information includes at least one of the following: tunnel crack distribution information and broken area information.

9. The device according to claim 6, characterized in that The determining unit is specifically used for: According to the tunnel stress information, tunnel deformation information, tunnel damage morphology information and the damage information of the cohesive unit corresponding to the PPV value of each of the mine earthquake waveform sample data, the tunnel damage degree corresponding to the PPV value of each of the mine earthquake waveform sample data is evaluated to obtain the mapping relationship between the PPV value and the tunnel damage degree; Establish a mapping relationship between the degree of tunnel damage and the intensity of mine earthquake according to the mine earthquake classification standard; The mapping relationship between the PPV value and the mine earthquake intensity is determined based on the mapping relationship between the PPV value and the degree of tunnel damage, and the mapping relationship between the degree of tunnel damage and the mine earthquake intensity.

10. An electronic device, characterized in that: include: A processor, and a memory communicatively connected to the processor; The memory stores computer-executable instructions; The processor executes the computer-executable instructions stored in the memory to implement the method according to any one of claims 1 to 4.

Images