A method and device for calculating mine earthquake energy considering stress wave propagation attenuation

Abstract

The present invention provides a method and device for calculating the energy of a mine earthquake that takes into account the attenuation of stress wave propagation, the method comprising: obtaining the distance between a plurality of sensors arranged in a monitoring area and a target source in the monitoring area; obtaining waveform data obtained by each sensor monitoring the target source, and calculating the initial total energy of the mine earthquake and the peak vibration velocity corresponding to each sensor according to the waveform data; obtaining the material attenuation coefficient of the stress wave of the mine earthquake according to each peak vibration velocity and the distance between each sensor and the target source; and obtaining the corrected total energy of the mine earthquake according to the material attenuation coefficient and the initial total energy of the mine earthquake. The embodiment of the present application calculates the material attenuation coefficient of the stress wave of the mine earthquake according to the waveform data monitored by a plurality of sensors at different distances from the source by considering the attenuation of stress wave propagation, and takes the material attenuation coefficient into consideration in the calculation of the energy of the mine earthquake, so as to more accurately calculate the true energy value of the target source.

Description

Technical Field

[0001] The present invention relates to the technical field of mine earthquake monitoring, and in particular to a mine earthquake energy calculation method and device taking into account stress wave propagation attenuation. Background Art

[0002] Mine tremors, also known as mining-induced earthquakes, are vibration phenomena caused by the rapid release of elastic energy of rock masses disturbed by mining that accompanies the development of underground solid resources. When a mine tremor occurs, it releases energy and produces vibration effects on underground or surface structures. Accurately assessing the energy generated by mine tremors is of great significance for further understanding of mine tremors and mine safety prevention and control. However, when calculating the energy of mine tremors in the field of mine tremors, energy parameter correction is required. The energy parameters are set based on experience and there is no scientific reference basis, which results in an error between the calculated energy and the actual energy.

[0003] Therefore, the prior art has defects and needs to be improved and developed. Summary of the invention

[0004] The present application provides a method for calculating mine earthquake energy taking into account the attenuation of stress wave propagation, so as to solve the technical problem that when calculating mine earthquake energy in related technologies, energy parameters need to be corrected, while the energy parameters are set based on experience without scientific reference basis, resulting in an error between the calculated energy and the actual energy.

[0005] To achieve the above objectives, this application adopts the following technical solutions:

[0006] The first aspect of the present application provides a method for calculating mining earthquake energy taking into account stress wave propagation attenuation, including:

[0007] Obtaining the distances between a number of sensors deployed in the monitoring area and the target seismic source in the monitoring area;

[0008] Acquire waveform data obtained by each sensor monitoring the target earthquake source, and calculate the total energy and peak vibration velocity of the initial mining earthquake corresponding to each sensor according to the waveform data;

[0009] The material attenuation coefficient of the mine earthquake stress wave is obtained according to each peak vibration velocity and the distance between each sensor and the target seismic source;

[0010] The corrected total energy of the mine tremor is obtained according to the material attenuation coefficient and the initial total energy of the mine tremor.

[0011] In one implementation, obtaining the distances between a plurality of sensors arranged in the monitoring area and a target seismic source in the monitoring area includes:

[0012] Obtaining location information of several sensors deployed in the monitoring area and source location information of target sources in the monitoring area;

[0013] Establishing a three-dimensional spatial coordinate system, obtaining the spatial position coordinates of each sensor in the three-dimensional spatial coordinate system according to the position information, and obtaining the source spatial coordinates of the target source in the three-dimensional spatial coordinate system according to the source position information;

[0014] Obtaining the distance between each sensor and the target seismic source according to each of the spatial position coordinates and the seismic source spatial coordinates;

[0015] The sensors are arranged on the ground within the monitoring area, and the distances between each sensor and the target seismic source are different.

[0016] In one implementation, the step of acquiring waveform data obtained by each sensor monitoring the target earthquake source and calculating the total energy and peak vibration velocity of the initial mining earthquake corresponding to each sensor according to the waveform data includes:

[0017] Acquire waveform data obtained by each sensor monitoring the target seismic source, and obtain a function of the waveform data changing with time according to the waveform data;

[0018] The energy density corresponding to each sensor is calculated according to the function of the waveform data corresponding to each sensor changing with time, and the energy density is calculated by spherical integration to obtain the total energy of the initial mine earthquake corresponding to each sensor;

[0019] The peak vibration velocity of each sensor is obtained according to the function of the waveform data corresponding to each sensor changing with time.

[0020] In one implementation, the sensor is a three-component sensor; the step of acquiring waveform data obtained by each sensor monitoring the target seismic source and obtaining a function of the waveform data changing with time according to the waveform data includes:

[0021] Acquire waveform data obtained by each sensor monitoring the target earthquake source according to a pre-constructed three-dimensional coordinate system, wherein the xy plane of the three-dimensional coordinate system is on a horizontal plane, the x-axis direction points to the earthquake source direction, the y-axis direction is perpendicular to the x-axis direction, and the z-axis direction is a vertically upward direction, and the waveform data includes waveform data in the x-axis direction, waveform data in the y-axis direction, and waveform data in the z-axis direction;

[0022] According to the waveform data, an x-axis direction variation function of the x-axis direction waveform data varying with time, a y-axis direction variation function of the y-axis direction waveform data varying with time, and a z-axis direction variation function of the z-axis direction waveform data varying with time are obtained.

[0023] In one implementation, the energy density corresponding to each sensor is calculated according to the function of the waveform data corresponding to each sensor changing with time, and the energy density is spherically integrated to obtain the total energy of the initial mine tremor corresponding to each sensor, including:

[0024] Obtain the rock density and the longitudinal wave velocity of the stress wave in the rock mass in the monitoring area, and obtain the energy density corresponding to each sensor according to the function of the rock density, the propagation velocity of the stress wave in the rock mass and the waveform data corresponding to each sensor changing with time;

[0025] According to the distance between each sensor and the target earthquake source, the energy density is calculated by spherical integration to obtain the total energy of the initial mine earthquake corresponding to each sensor.

[0026] In one implementation, obtaining the material attenuation coefficient of the mine earthquake stress wave according to each peak vibration velocity and the distance between each sensor and the target seismic source includes:

[0027] The spherical wave geometric attenuation coefficient is obtained according to the distance between each sensor and the target source;

[0028] According to each of the peak vibration velocities, the distance between each sensor and the target seismic source, and the spherical wave geometric attenuation coefficient, a pre-constructed attenuation formula is used for fitting to obtain the material attenuation coefficient of the mine earthquake stress wave.

[0029] In one implementation, obtaining the corrected total energy of the mine tremor according to the material attenuation coefficient and the initial total energy of the mine tremor includes:

[0030] Calculate the average value of the distance between each sensor and the target earthquake source to obtain the average distance, and calculate the average value of the initial mining earthquake total energy of each sensor to obtain the average energy;

[0031] The preset earthquake source equivalent radius is obtained, and the corrected total energy of the mine earthquake is obtained according to the material attenuation coefficient, the average distance, the earthquake source equivalent radius and the average energy.

[0032] A second aspect of the present application provides a mining earthquake energy calculation device taking into account stress wave propagation attenuation, including:

[0033] A distance acquisition module is used to acquire the distance between a number of sensors deployed in the monitoring area and the target source in the monitoring area;

[0034] A calculation module is used to obtain waveform data obtained by each sensor monitoring the target earthquake source, and calculate the total energy and peak vibration velocity of the initial mining earthquake corresponding to each sensor according to the waveform data;

[0035] A coefficient acquisition module is used to obtain the material attenuation coefficient of the mine earthquake stress wave according to each peak vibration velocity and the distance between each sensor and the target seismic source;

[0036] The correction module is used to obtain the corrected total energy of the mine tremor according to the material attenuation coefficient and the initial total energy of the mine tremor.

[0037] A third aspect of the present application provides a terminal, which includes a memory, a processor, and a mine tremor energy calculation program that takes into account the attenuation of stress wave propagation, which is stored in the memory and can be run on the processor. When the processor executes the mine tremor energy calculation program that takes into account the attenuation of stress wave propagation, the steps of the mine tremor energy calculation method that takes into account the attenuation of stress wave propagation as described above are implemented.

[0038] A fourth aspect of the present application provides a computer-readable storage medium, on which is stored a mine tremor energy calculation program that takes into account the attenuation of stress wave propagation. When the mine tremor energy calculation program that takes into account the attenuation of stress wave propagation is executed by a processor, the steps of the mine tremor energy calculation method that takes into account the attenuation of stress wave propagation as described above are implemented.

[0039] Beneficial effects of the present invention: The embodiment of the present application obtains the distance between several sensors arranged in the monitoring area and the target earthquake source in the monitoring area; obtains the waveform data obtained by each sensor monitoring the target earthquake source, and calculates the initial mine earthquake total energy and peak vibration velocity corresponding to each sensor according to the waveform data; obtains the material attenuation coefficient of the mine earthquake stress wave according to each peak vibration velocity and the distance between each sensor and the target earthquake source; obtains the corrected mine earthquake total energy according to the material attenuation coefficient and the initial mine earthquake total energy. The embodiment of the present application calculates the material attenuation coefficient of the mine earthquake stress wave according to the waveform data monitored by several sensors at different distances from the earthquake source by considering the attenuation of stress wave propagation, and takes the material attenuation coefficient into account in the mine earthquake energy calculation, so as to more accurately calculate the true energy value of the target earthquake source. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0041] Figure 1 It is a flow chart of a preferred embodiment of the method for calculating mine earthquake energy taking into account stress wave propagation attenuation in the present invention.

[0042] Figure 2 It is a schematic diagram of laying out three-component sensors on the ground.

[0043] Figure 3 This is an example diagram of waveforms detected by a three-component sensor in three directions.

[0044] Figure 4 It is a schematic diagram of calculating kinetic energy by spherical integration in three-dimensional space.

[0045] Figure 5 It is a schematic diagram of the vibration velocity of particles in three directions on the monitoring surface when stress waves pass through the rock mass.

[0046] Figure 6 This is an example diagram of the peak vibration velocity attenuation fitting curve.

[0047] Figure 7 It is a functional principle block diagram of a preferred embodiment of a mine earthquake energy calculation device that takes stress wave propagation attenuation into consideration in the present invention.

[0048] Figure 8 It is a functional principle block diagram of a preferred embodiment of the terminal in the present invention. DETAILED DESCRIPTION

[0049] In order to make the purpose, technical solution and advantages of the present invention clearer and more specific, the present invention is further described in detail below with reference to the accompanying drawings and examples. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not used to limit the present invention.

[0050] At present, there are still the following problems in calculating the energy of mine earthquakes:

[0051] First, the microseismic monitoring technology commonly used in mine earthquakes is deployed in the surrounding rock or on the ground for a long time, and the sensor response lacks real-time correction conditions, so it is impossible to accurately measure the true value of the vibration parameters of the mine earthquake stress wave;

[0052] Second, energy parameters need to be set when calculating source energy, and the selection of parameter values ​​is mostly determined by experience without scientific reference, resulting in errors between the calculated energy and the actual energy;

[0053] Third, at present, when calculating the energy of mine tremors, the influence of the material attenuation coefficient suitable for the geological conditions of coal mines on mine tremors is not taken into account. The calculation results of mine tremor energy have large errors, which seriously affect our accurate understanding and assessment of mine tremors.

[0054] See also Figure 1 The method for calculating mining earthquake energy considering stress wave propagation attenuation according to an embodiment of the present invention comprises the following steps:

[0055] Step S100: Obtain the distances between a number of sensors deployed in the monitoring area and the target seismic source in the monitoring area.

[0056] In one embodiment of the present application, the step S100 specifically includes:

[0057] Step S110, obtaining location information of a number of sensors deployed in the monitoring area and source location information of a target source in the monitoring area;

[0058] Step S120, establishing a three-dimensional spatial coordinate system, obtaining the spatial position coordinates of each sensor in the three-dimensional spatial coordinate system according to the position information, and obtaining the source spatial coordinates of the target source in the three-dimensional spatial coordinate system according to the source position information;

[0059] Step S130: Obtain the distance between each sensor and the target seismic source according to the spatial position coordinates and the seismic source spatial coordinates.

[0060] The sensors are arranged on the ground within the monitoring area, and the distances between each sensor and the target seismic source are different.

[0061] Specifically, n (n is not less than 5) three-component sensors are deployed in the monitoring area; the longitude and latitude of each sensor can be obtained by using GPS positioning and other methods, and the longitude and latitude are used as the location information, and the target source is used as the known source, that is, the source location information of the known target source, which can also be the longitude and latitude of the source. A spatial three-dimensional coordinate system is established, and the location information is converted into the spatial three-dimensional coordinate system. The spatial position coordinates of each sensor are A and B respectively. 1 , A 2 …A n , the spatial coordinate of the earthquake source is B, and the distance between each sensor and the earthquake source is calculated as r 1 、r 2 …r n .

[0062] The embodiment of the present application can realize the calculation of mine earthquake energy by installing sensors on the ground. Specifically, the calculation of mine earthquake energy can be realized by arranging sensors at different distances from the earthquake source. The sensor arrangement is easy, the disassembly is convenient, and the monitoring cost is low. In addition, the real-time calibration of the sensor can be realized. Since the sensor that is easy to install on the ground is adopted, such as the vibration pickup, the problem of the inconsistency of the coupling of the underground surrounding rock arrangement sensor for a long time cannot be corrected. The vibration response of the sensor designed by the present invention can be corrected at any time to ensure that the vibration parameters monitored by the sensor are true, and provide a reliable basis for the calculation of mine earthquake energy.

[0063] like Figure 1 As shown, the method for calculating mining earthquake energy taking into account stress wave propagation attenuation also includes the following steps:

[0064] Step S200, acquiring waveform data obtained by each sensor monitoring the target earthquake source, and calculating the total energy and peak vibration velocity of the initial mining earthquake corresponding to each sensor according to the waveform data.

[0065] In one embodiment of the present application, step S200 specifically includes:

[0066] Step S210, acquiring waveform data obtained by each sensor monitoring the target seismic source, and obtaining a function of waveform data changing with time according to the waveform data;

[0067] Step S220, calculating the energy density corresponding to each sensor according to the function of the waveform data corresponding to each sensor changing with time, performing spherical integral calculation on the energy density, and obtaining the total energy of the initial mine earthquake corresponding to each sensor;

[0068] Step S230: Obtain the peak vibration velocity of each sensor according to a function of the waveform data corresponding to each sensor varying with time.

[0069] Specifically, the mine earthquake source generated by the coal mine is monitored, and the waveform data monitored by each group of sensors are obtained as f 1 、f 2 …f n ; Calculate different distances r from the earthquake source i The total energy of the mine earthquake corresponding to the sensor E i .

[0070] The embodiment of the present application assumes that the earthquake source propagates as a spherical wave. The total energy of the earthquake source propagated to a distance r from the earthquake source can be calculated by spherical integration of the unit density energy to obtain the total energy of the earthquake source.

[0071] In one embodiment, the sensor is a three-component sensor; the step S210 specifically includes:

[0072] Step S211, acquiring waveform data obtained by each sensor monitoring the target seismic source according to a pre-constructed three-dimensional coordinate system, wherein the xy plane of the three-dimensional coordinate system is on a horizontal plane, the x-axis direction points to the seismic source direction, the y-axis direction is perpendicular to the x-axis direction, and the z-axis direction is a vertically upward direction, and the waveform data includes waveform data in the x-axis direction, waveform data in the y-axis direction, and waveform data in the z-axis direction;

[0073] Step S212, obtaining, according to the waveform data, an x-axis direction variation function of the x-axis direction waveform data varying with time, a y-axis direction variation function of the y-axis direction waveform data varying with time, and a z-axis direction variation function of the z-axis direction waveform data varying with time.

[0074] The embodiment of the present application obtains the peak vibration velocity of each sensor through the x-axis direction change function, the y-axis direction change function and the z-axis direction change function.

[0075] In one embodiment, the step S220 specifically includes:

[0076] Step S221, obtaining the rock density and the longitudinal wave velocity of the stress wave in the rock mass in the monitoring area, and obtaining the energy density corresponding to each sensor according to the function of the rock density, the propagation velocity of the stress wave in the rock mass and the waveform data corresponding to each sensor changing with time;

[0077] Step S222: Perform spherical integral calculation on the energy density according to the distance between each sensor and the target earthquake source to obtain the total energy of the initial mine earthquake corresponding to each sensor.

[0078] Specifically, the total energy density e of the stress wave passing through a unit rock mass is equal to the kinetic energy e carried by the stress wave. k and stress wave propagation process, particle vibration process, elastic potential energy stored in the rock e u , that is e=e k +e u .

[0079] The kinetic energy expression of the stress wave at a certain time t during the propagation of the rock mass is:

[0080]

[0081] Where ρ represents the rock mass density, C represents the propagation speed of stress waves in the rock mass, and f x (t) represents the function of the waveform data of the x direction monitored by the sensor changing with time, f y (t) represents the function of the waveform data of the y direction monitored by the sensor changing with time, f z (t) represents the function of the waveform data monitored by the sensor in the z direction changing with time; according to the basic theory of stress waves, the kinetic energy caused by stress waves in the rock mass is equal to the strain energy stored in the rock mass, that is, e u =e k , the total energy density is:

[0082] e=e u +e k =ρdtC(f x 2 (t)+f y 2 (t)+f z 2 (t)) (2)

[0083] Assuming that the earthquake source propagates as a spherical wave, the total energy of the earthquake source propagated to a distance r from the earthquake source can be calculated by spherical integration of the unit density energy to obtain the total energy of the earthquake source E:

[0084]

[0085] Where r is the distance between the sensor and the target source, and then the total energy E of the mine earthquake corresponding to each sensor is calculated. 1、E 2 …E n .

[0086] When calculating the peak vibration velocity v of the earthquake source monitored by each sensor i When , the calculation method is as follows:

[0087] The waveform data monitored by the three-component velocity sensor over time is f x (t),f y (t),f z (t), the calculation formula of the combined velocity waveform function is:

[0088]

[0089] The peak vibration velocity of each sensor is v i =max(f i (t)).

[0090] like Figure 1 As shown, the method for calculating mining earthquake energy taking into account stress wave propagation attenuation also includes the following steps:

[0091] Step S300: obtaining the material attenuation coefficient of the mine earthquake stress wave according to each peak vibration velocity and the distance between each sensor and the target seismic source.

[0092] In one embodiment of the present application, step S300 specifically includes:

[0093] Step S310, obtaining a spherical wave geometric attenuation coefficient according to the distance between each sensor and the target source;

[0094] Step S320: fitting a pre-constructed attenuation formula according to each of the peak vibration velocities, the distance between each sensor and the target seismic source, and the spherical wave geometric attenuation coefficient to obtain the material attenuation coefficient of the mine earthquake stress wave.

[0095] Specifically, according to the stress wave propagation attenuation theory, the geometric attenuation caused by the diffusion of the wavefront in the vibration velocity propagation attenuation process and the material attenuation caused by the damping of the rock material, for spherical wave propagation, the theoretical formula for the attenuation of the vibration velocity with the propagation distance is:

[0096]

[0097] In the formula, v is the peak vibration velocity, A is a constant, α is the material attenuation coefficient, and r is the distance from the earthquake source. is the geometric attenuation coefficient of spherical wave.

[0098] The calculated v i With r i Based on the formula The formal fitting is performed to obtain A and α.

[0099] The embodiment of the present application takes into account the attenuation characteristics of stress waves, and then inverts the energy at the source of the mine earthquake. That is, the material attenuation characteristics of stress waves in rock materials are taken into account, which improves the scientific nature of the calculation of mine earthquake energy, and then can accurately calculate the true energy value of the source of the mine earthquake.

[0100] like Figure 1 As shown, the method for calculating mining earthquake energy taking into account stress wave propagation attenuation also includes the following steps:

[0101] Step S400: obtaining a corrected total energy of the mine tremor according to the material attenuation coefficient and the initial total energy of the mine tremor.

[0102] In one embodiment of the present application, step S400 specifically includes:

[0103] Step S410, calculating the average value of the distance between each sensor and the target seismic source to obtain the average distance between the sensor and the target seismic source, and calculating the average value of the initial mining earthquake total energy of each sensor to obtain the average energy;

[0104] Step S420, obtaining a preset earthquake source equivalent radius, and obtaining a corrected total mine earthquake energy according to the material attenuation coefficient, the average distance, the earthquake source equivalent radius and the average energy.

[0105] Specifically, the equivalent radius r of the earthquake source is pre-set according to the size of the earthquake source. 0 , and r 0 < <r i The total energy calculation formula of mine earthquake considering attenuation is:

[0106]

[0107] In the formula, E represents the total energy of the modified mine earthquake, α represents the material attenuation coefficient, represents the average distance between the sensor and the target source, represents the average energy.

[0108] Specific embodiments are listed below for illustration.

[0109] Figure 2 This is a schematic diagram of sensor deployment. 10 sensors are deployed on the ground. Figure 2 Where S1 and S2 represent different propagation wave vibration planes, and R1 and R2 represent the distances between different propagation wave vibration planes and the earthquake source O. The three-dimensional coordinates of each sensor are measured, and the coordinates of each sensor are shown in Table 1;

[0110]

[0111] Table 1

[0112] The spatial coordinates of the earthquake source are obtained as (4323033, 382420, -140), and the distance between each sensor and the earthquake source is calculated. The results are shown in Table 1.

[0113] The stress waves generated by the mine earthquake source are monitored. Each group of sensors obtains the three-component waveform data of the source. The functions of the waveform data changing with time are f 1 -f 10 , Figure 3 It is the waveform data monitored by sensor 1.

[0114] like Figure 4 As shown, the total energy E of the mine earthquake is calculated according to formula (3): 1 =1.72×10 7 J, calculate the combined velocity function of each sensor according to formula (4), and calculate the peak vibration velocity monitored by each sensor, such as Figure 5 shown.

[0115] According to formula (5), the calculated peak vibration velocity is fitted with the distance from the earthquake source. The fitting result is as follows: Figure 6 As shown. The calculated values ​​of A and α are 146 and -2.455×10 -4 .

[0116] The effective radius of the earthquake source is 5m, the average distance between the sensor and the earthquake source is 665.274m, and the average energy is 1.68×10 7 J. According to formula (6), the total energy of the mine earthquake with attenuation considered in the earthquake source monitored this time is 1.98×10 7 J.

[0117] In one embodiment, if Figure 7 As shown, based on the above-mentioned method for calculating the energy of a mine earthquake taking into account the attenuation of stress wave propagation, the present invention also provides a device for calculating the energy of a mine earthquake taking into account the attenuation of stress wave propagation, including:

[0118] The distance acquisition module 100 is used to acquire the distance between a plurality of sensors arranged in the monitoring area and the target source in the monitoring area;

[0119] A calculation module 200 is used to obtain waveform data obtained by each sensor monitoring the target earthquake source, and calculate the total energy and peak vibration velocity of the initial mining earthquake corresponding to each sensor according to the waveform data;

[0120] A coefficient acquisition module 300, for obtaining a material attenuation coefficient of a mine earthquake stress wave according to each peak vibration velocity and the distance between each sensor and a target seismic source;

[0121] The correction module 400 is used to obtain the corrected total energy of the mine tremor according to the material attenuation coefficient and the initial total energy of the mine tremor.

[0122] It should be noted that the aforementioned explanation of the embodiment of the method for calculating the energy of a mine tremor taking into account the attenuation of stress wave propagation is also applicable to the device for calculating the energy of a mine tremor taking into account the attenuation of stress wave propagation in this embodiment, and will not be repeated here.

[0123] Figure 8 A schematic diagram of the structure of a terminal provided in an embodiment of the present application. The terminal may include:

[0124] A memory 501 , a processor 502 , and a computer program stored in the memory 501 and executable on the processor 502 .

[0125] When the processor 502 executes the program, the method for calculating the energy of a mine earthquake taking into account the attenuation of stress wave propagation provided in the above embodiment is implemented.

[0126] Furthermore, the terminal further includes:

[0127] The communication interface 503 is used for communication between the memory 501 and the processor 502 .

[0128] The memory 501 is used to store computer programs that can be executed on the processor 502 .

[0129] The memory 501 may include a high-speed RAM memory, and may also include a non-volatile memory (non-volatile memory), such as at least one disk memory.

[0130] If the memory 501, the processor 502 and the communication interface 503 are implemented independently, the communication interface 503, the memory 501 and the processor 502 can be connected to each other through a bus and communicate with each other. The bus can be an Industry Standard Architecture (ISA) bus, a Peripheral Component (PCI) bus or an Extended Industry Standard Architecture (EISA) bus. The bus can be divided into an address bus, a data bus, a control bus, etc. For ease of representation, only one line is used in the figure, but it does not mean that there is only one bus or one type of bus.

[0131] Optionally, in a specific implementation, if the memory 501, the processor 502 and the communication interface 503 are integrated on a chip, the memory 501, the processor 502 and the communication interface 503 can communicate with each other through an internal interface.

[0132] The processor 502 may be a central processing unit (CPU), or an application specific integrated circuit (ASIC), or one or more integrated circuits configured to implement the embodiments of the present application.

[0133] This embodiment also provides a computer-readable storage medium on which a computer program is stored. When the program is executed by a processor, the above-mentioned method for calculating mine earthquake energy taking into account the attenuation of stress wave propagation is implemented.

[0134] In the description of this specification, 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 application. 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 N 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, without contradiction.

[0135] 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 technical features indicated. Therefore, a feature defined as "first" or "second" may explicitly or implicitly include at least one of the features. In the description of this application, "N" means at least two, such as two, three, etc., unless otherwise clearly and specifically defined.

[0136] Any process or method description in a flowchart or otherwise described herein may be understood to represent a module, fragment or portion of code comprising one or N executable instructions for implementing the steps of a custom logical function or process, and the scope of the preferred embodiments of the present application includes additional implementations, in which the order shown or discussed may not be followed, 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 application belong.

[0137] 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 read instructions from an instruction execution system, device or apparatus and execute 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 N wirings (electronic devices), 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 by optically scanning the paper or other medium and then editing, interpreting or otherwise processing in a suitable manner if necessary and then storing it in a computer memory.

[0138] It should be understood that the various parts of the present application can be implemented by hardware, software, firmware or a combination thereof. In the above embodiment, the N steps or methods can be implemented by software or firmware stored in a memory and executed by a suitable instruction execution system. If implemented by hardware, as in another embodiment, it can be implemented by any one of the following technologies known in the art or their combination: 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.

[0139] 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. When the program is executed, it includes one or a combination of the steps of the method embodiment.

[0140] In addition, each functional unit in each embodiment of the present application may be integrated into a processing module, or each unit may exist physically separately, or two or more units may be integrated into one module. The above-mentioned integrated module may be implemented in the form of hardware or in the form of a software functional module. If the integrated module is implemented in the form of a software functional module and sold or used as an independent product, it may also be stored in a computer-readable storage medium.

[0141] The storage medium mentioned above may be a read-only memory, a magnetic disk or an optical disk, etc. Although the embodiments of the present application have been shown and described above, it can be understood that the above embodiments are exemplary and cannot be understood as limiting the present application. A person of ordinary skill in the art may change, modify, replace and modify the above embodiments within the scope of the present application.

[0142] In summary, the present invention discloses a method and device for calculating the energy of a mine earthquake that takes into account the attenuation of stress wave propagation, the method comprising: obtaining the distance between a number of sensors arranged in a monitoring area and a target source in the monitoring area; obtaining waveform data obtained by each sensor monitoring the target source, and calculating the initial total energy of the mine earthquake and the peak vibration velocity corresponding to each sensor according to the waveform data; obtaining the material attenuation coefficient of the stress wave of the mine earthquake according to each peak vibration velocity and the distance between each sensor and the target source; and obtaining the corrected total energy of the mine earthquake according to the material attenuation coefficient and the initial total energy of the mine earthquake. The embodiment of the present application calculates the material attenuation coefficient of the stress wave of the mine earthquake according to the waveform data monitored by a number of sensors at different distances from the source by considering the attenuation of stress wave propagation, and takes the material attenuation coefficient into account in the calculation of the energy of the mine earthquake, so as to more accurately calculate the true energy value of the target source.

[0143] It should be understood that the application of the present invention is not limited to the above examples. For ordinary technicians in this field, improvements or changes can be made based on the above description. All these improvements and changes should fall within the scope of protection of the claims attached to the present invention.

Claims

1. A method for calculating mining earthquake energy taking into account stress wave propagation attenuation, characterized in that: include: Obtaining the distances between a number of sensors deployed in the monitoring area and the target seismic source in the monitoring area; Acquire waveform data obtained by each sensor monitoring the target earthquake source, and calculate the total energy and peak vibration velocity of the initial mining earthquake corresponding to each sensor according to the waveform data; The material attenuation coefficient of the mine earthquake stress wave is obtained according to each peak vibration velocity and the distance between each sensor and the target seismic source; Obtaining a corrected total energy of the mine tremor according to the material attenuation coefficient and the initial total energy of the mine tremor; The obtaining of the distances between a plurality of sensors arranged in the monitoring area and the target seismic source in the monitoring area comprises: Obtaining location information of several sensors deployed in the monitoring area and source location information of target sources in the monitoring area; Establishing a three-dimensional spatial coordinate system, obtaining the spatial position coordinates of each sensor in the three-dimensional spatial coordinate system according to the position information of the plurality of sensors, and obtaining the spatial coordinates of the target seismic source in the three-dimensional spatial coordinate system according to the seismic source position information; Obtaining the distance between each sensor and the target seismic source according to each of the spatial position coordinates and the seismic source spatial coordinates; The sensors are arranged on the ground within the monitoring area, and the distances between each sensor and the target seismic source are different; The step of acquiring waveform data obtained by each sensor monitoring the target earthquake source and calculating the total energy and peak vibration velocity of the initial mining earthquake corresponding to each sensor according to the waveform data includes: Acquire waveform data obtained by each sensor monitoring the target seismic source, and obtain a function of the waveform data changing with time according to the waveform data; The energy density corresponding to each sensor is calculated according to the function of the waveform data corresponding to each sensor changing with time, and the energy density is calculated by spherical integration to obtain the total energy of the initial mine earthquake corresponding to each sensor; The peak vibration velocity of each sensor is obtained according to the function of the waveform data corresponding to each sensor changing with time; The sensor is a three-component sensor; the waveform data obtained by each sensor monitoring the target seismic source is obtained, and a function of the waveform data changing with time is obtained according to the waveform data, including: Acquire waveform data obtained by each sensor monitoring the target earthquake source according to a pre-constructed three-dimensional coordinate system, wherein the xy plane of the three-dimensional coordinate system is on a horizontal plane, the x-axis direction points to the earthquake source direction, the y-axis direction is perpendicular to the x-axis direction, and the z-axis direction is a vertically upward direction, and the waveform data includes waveform data in the x-axis direction, waveform data in the y-axis direction, and waveform data in the z-axis direction; According to the waveform data, an x-axis direction variation function of the x-axis direction waveform data varying with time, a y-axis direction variation function of the y-axis direction waveform data varying with time, and a z-axis direction variation function of the z-axis direction waveform data varying with time are obtained; The energy density corresponding to each sensor is calculated according to the function of the waveform data corresponding to each sensor changing with time. The energy density is calculated by spherical integration to obtain the total energy of the initial mine earthquake corresponding to each sensor, including: Obtain the rock density and the longitudinal wave velocity of the stress wave in the rock mass in the monitoring area, and obtain the energy density corresponding to each sensor according to the function of the rock density, the propagation velocity of the stress wave in the rock mass and the waveform data corresponding to each sensor changing with time; According to the distance between each sensor and the target earthquake source, the energy density is calculated by spherical integration to obtain the total energy of the initial mine earthquake corresponding to each sensor; The corrected total energy of the mine tremor is obtained according to the material attenuation coefficient and the initial total energy of the mine tremor, including: Calculate the average value of the distance between each sensor and the target earthquake source to obtain the average distance, and calculate the average value of the initial mining earthquake total energy of each sensor to obtain the average energy; Obtaining a preset earthquake source equivalent radius, and obtaining a corrected total energy of the mine earthquake according to the material attenuation coefficient, the average distance, the earthquake source equivalent radius and the average energy; The calculation formula of the energy density is: ,in, It represents the elastic potential energy stored in the rock during the particle vibration process during stress wave propagation. represents the kinetic energy carried by the stress wave, represents the rock mass density, It represents the propagation speed of stress wave in rock mass. It represents the function of the waveform data in the x direction monitored by the sensor changing with time. It represents the function of the waveform data in the y direction detected by the sensor changing with time. It represents the function of the waveform data in the z direction monitored by the sensor changing with time; Assuming that the earthquake source propagates as a spherical wave, the calculation formula for the total energy of the initial mining earthquake is: ， Where r is the distance between the sensor and the target source; The calculation formula of the function of the waveform data corresponding to each sensor changing with time is: ; The equivalent radius r0 of the earthquake source is pre-set according to the size of the earthquake source, and , the calculation formula of the modified total energy of mining earthquake is: , where represents the total energy of the modified mine earthquake, represents the material attenuation coefficient, represents the average distance between the sensor and the target source, Represents the average energy.

2. The method for calculating mining earthquake energy considering stress wave propagation attenuation according to claim 1 is characterized in that: The material attenuation coefficient of the mine earthquake stress wave is obtained according to each peak vibration velocity and the distance between each sensor and the target seismic source, including: The spherical wave geometric attenuation coefficient is obtained according to the distance between each sensor and the target source; According to each of the peak vibration velocities, the distance between each sensor and the target seismic source, and the spherical wave geometric attenuation coefficient, a pre-constructed attenuation formula is used for fitting to obtain the material attenuation coefficient of the mine earthquake stress wave.

3. A mining earthquake energy calculation device taking into account stress wave propagation attenuation, characterized in that: include: A distance acquisition module is used to acquire the distance between a number of sensors deployed in the monitoring area and the target source in the monitoring area; A calculation module is used to obtain waveform data obtained by each sensor monitoring the target earthquake source, and calculate the total energy and peak vibration velocity of the initial mining earthquake corresponding to each sensor according to the waveform data; A coefficient acquisition module is used to obtain the material attenuation coefficient of the mine earthquake stress wave according to each peak vibration velocity and the distance between each sensor and the target seismic source; A correction module, used for obtaining a corrected total energy of the mine tremor according to the material attenuation coefficient and the initial total energy of the mine tremor; The obtaining of the distances between a plurality of sensors arranged in the monitoring area and the target seismic source in the monitoring area comprises: Obtaining location information of several sensors deployed in the monitoring area and source location information of target sources in the monitoring area; Establishing a three-dimensional spatial coordinate system, obtaining the spatial position coordinates of each sensor in the three-dimensional spatial coordinate system according to the position information of the plurality of sensors, and obtaining the spatial coordinates of the target seismic source in the three-dimensional spatial coordinate system according to the seismic source position information; Obtaining the distance between each sensor and the target seismic source according to each of the spatial position coordinates and the seismic source spatial coordinates; The sensors are arranged on the ground within the monitoring area, and the distances between each sensor and the target seismic source are different; The step of acquiring waveform data obtained by each sensor monitoring the target earthquake source and calculating the total energy and peak vibration velocity of the initial mining earthquake corresponding to each sensor according to the waveform data includes: Acquire waveform data obtained by each sensor monitoring the target seismic source, and obtain a function of the waveform data changing with time according to the waveform data; The energy density corresponding to each sensor is calculated according to the function of the waveform data corresponding to each sensor changing with time, and the energy density is calculated by spherical integration to obtain the total energy of the initial mine earthquake corresponding to each sensor; The peak vibration velocity of each sensor is obtained according to the function of the waveform data corresponding to each sensor changing with time; The sensor is a three-component sensor; the waveform data obtained by each sensor monitoring the target seismic source is obtained, and a function of the waveform data changing with time is obtained according to the waveform data, including: Acquire waveform data obtained by each sensor monitoring the target earthquake source according to a pre-constructed three-dimensional coordinate system, wherein the xy plane of the three-dimensional coordinate system is on a horizontal plane, the x-axis direction points to the earthquake source direction, the y-axis direction is perpendicular to the x-axis direction, and the z-axis direction is a vertically upward direction, and the waveform data includes waveform data in the x-axis direction, waveform data in the y-axis direction, and waveform data in the z-axis direction; According to the waveform data, an x-axis direction variation function of the x-axis direction waveform data varying with time, a y-axis direction variation function of the y-axis direction waveform data varying with time, and a z-axis direction variation function of the z-axis direction waveform data varying with time are obtained; The energy density corresponding to each sensor is calculated according to the function of the waveform data corresponding to each sensor changing with time. The energy density is calculated by spherical integration to obtain the total energy of the initial mine earthquake corresponding to each sensor, including: Obtain the rock density and the longitudinal wave velocity of the stress wave in the rock mass in the monitoring area, and obtain the energy density corresponding to each sensor according to the function of the rock density, the propagation velocity of the stress wave in the rock mass and the waveform data corresponding to each sensor changing with time; According to the distance between each sensor and the target earthquake source, the energy density is calculated by spherical integration to obtain the total energy of the initial mine earthquake corresponding to each sensor; The corrected total energy of the mine tremor is obtained according to the material attenuation coefficient and the initial total energy of the mine tremor, including: Calculate the average value of the distance between each sensor and the target earthquake source to obtain the average distance, and calculate the average value of the initial mining earthquake total energy of each sensor to obtain the average energy; Obtaining a preset earthquake source equivalent radius, and obtaining a corrected total energy of the mine earthquake according to the material attenuation coefficient, the average distance, the earthquake source equivalent radius and the average energy; The calculation formula of the energy density is: ,in, It represents the elastic potential energy stored in the rock during the particle vibration process during stress wave propagation. represents the kinetic energy carried by the stress wave, represents the rock mass density, It represents the propagation speed of stress wave in rock mass. It represents the function of the waveform data in the x direction monitored by the sensor changing with time. It represents the function of the waveform data in the y direction detected by the sensor changing with time. It represents the function of the waveform data in the z direction monitored by the sensor changing with time; Assuming that the earthquake source propagates as a spherical wave, the calculation formula for the total energy of the initial mining earthquake is: ， Where r is the distance between the sensor and the target source; The calculation formula of the function of the waveform data corresponding to each sensor changing with time is: ; The equivalent radius r0 of the earthquake source is pre-set according to the size of the earthquake source, and , the calculation formula of the modified total energy of mining earthquake is: , where represents the total energy of the modified mine earthquake, represents the material attenuation coefficient, represents the average distance between the sensor and the target source, Represents the average energy.

4. A terminal, characterized in that: The terminal includes a memory, a processor, and a mine tremor energy calculation program that takes into account stress wave propagation attenuation, which is stored in the memory and can be run on the processor. When the processor executes the mine tremor energy calculation program that takes into account stress wave propagation attenuation, the steps of the mine tremor energy calculation method that takes into account stress wave propagation attenuation as described in any one of claims 1-2 are implemented.

5. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a mine tremor energy calculation program that takes into account the attenuation of stress wave propagation. When the mine tremor energy calculation program that takes into account the attenuation of stress wave propagation is executed by a processor, the steps of the mine tremor energy calculation method that takes into account the attenuation of stress wave propagation as described in any one of claims 1-2 are implemented.

Images