Coal mine blasting longitudinal wave velocity measurement method and device and storage medium

Abstract

This application discloses a method and apparatus for measuring the longitudinal wave velocity in coal mine blasting. The method includes: acquiring the first arrival time of the blast wave to the source detector and the second arrival time to the target detector after the blasting of a pre-set explosive in a coal mine roadway, wherein there are multiple target detectors; calculating the distance between each target detector and the source detector to obtain the source-detector distance corresponding to each target detector; calculating the blast wave travel time corresponding to each target detector based on the first and second arrival times; calculating the initial longitudinal wave velocity corresponding to each target detector based on the blast wave travel time and source-detector distance corresponding to each target detector; constructing an initial longitudinal wave velocity sample set containing the initial longitudinal wave velocity corresponding to each target detector; removing error samples from the initial longitudinal wave velocity sample set according to the Laida criterion to obtain an effective longitudinal wave velocity sample set; and determining the blasting longitudinal wave velocity based on the effective longitudinal wave velocity sample set, thereby realizing the measurement of the longitudinal wave velocity in coal mine blasting.

Description

Technical Field

[0001] This application relates to the field of mine vibration monitoring technology, and in particular to a method and device for measuring longitudinal wave velocity in coal mine blasting, a storage medium, and a computer device. Background Technology

[0002] During coal mining, the extraction of coal leads to the redistribution of stress in the surrounding rock and the collapse and movement of overburden, generating microseismic events and releasing elastic vibration waves. The accuracy of microseismic source location is a crucial indicator of the quality of coal mine microseismic network monitoring, and largely determines the effectiveness of microseismic monitoring.

[0003] The accuracy of microseismic source location is significantly affected by the P-wave velocity. A reasonable P-wave velocity can ensure accurate source location, while an inaccurate velocity can lead to large location errors. However, in the field of microseismic monitoring in coal mines, the solution for the P-wave velocity has not been well resolved. Therefore, a method for measuring the P-wave velocity in coal mine blasting is needed. Summary of the Invention

[0004] In view of this, this application provides a method and apparatus for measuring longitudinal wave velocity in coal mine blasting, a storage medium, and a computer device, which realizes the measurement of longitudinal wave velocity in coal mine blasting.

[0005] According to one aspect of this application, a method for measuring longitudinal wave velocity in coal mine blasting is provided, the method comprising:

[0006] After the pre-set explosives are detonated in a coal mine roadway, the first arrival time of the blast wave to the source detector and the second arrival time to the target detector are obtained, wherein there are multiple target detectors.

[0007] The distance between each target detector and the source detector is calculated to obtain the source-detector distance corresponding to each target detector. Based on the first arrival time and the second arrival time, the explosion wave travel time corresponding to each target detector is calculated.

[0008] Based on the blast wave travel time corresponding to each target detector and the source-detector distance, the initial P-wave velocity corresponding to each target detector is calculated, and an initial P-wave velocity sample set containing the initial P-wave velocity corresponding to each target detector is constructed.

[0009] Error samples are removed from the initial P-wave velocity sample set according to the Laida criterion to obtain an effective P-wave velocity sample set, and the blasting P-wave velocity is determined based on the effective P-wave velocity sample set.

[0010] Optionally, the source detector is horizontally arranged on the coal wall in the coal mine roadway and at a predetermined distance from the blasting position of the predetermined explosive; the target detector is vertically arranged on the roof and / or floor of the coal mine roadway and within a predetermined range corresponding to the blasting position; the step of calculating the distance between each target detector and the source detector to obtain the source-detector distance corresponding to each target detector includes:

[0011] Obtain the three-dimensional coordinates of the source detector and the target detector;

[0012] Based on the three-dimensional coordinates of the source detector and the three-dimensional coordinates of the target detector, the source-detector distance between each target detector and the source detector is calculated.

[0013] Optionally, calculating the blast wave travel time corresponding to each target detector based on the first arrival time and the second arrival time includes:

[0014] The difference between the second duration and the first duration corresponding to each of the target detectors is calculated to obtain the blast wave travel time corresponding to each of the target detectors.

[0015] Optionally, calculating the initial P-wave velocity for each target detector based on the blast travel time and the source-detector distance for each target detector includes:

[0016] The ratio of the source-detector distance to the travel time of the blast wave is calculated for each target detector to obtain the initial longitudinal wave velocity for each target detector.

[0017] Optionally, the step of removing error samples from the initial P-wave velocity sample set according to the Laida criterion to obtain an effective P-wave velocity sample set includes:

[0018] Calculate the arithmetic mean, standard deviation, and absolute difference between each initial P-wave velocity and the arithmetic mean for all initial P-wave velocities in the initial P-wave velocity sample set.

[0019] If the ratio of the absolute difference to the standard deviation is greater than or equal to a preset multiple, the corresponding initial P-wave velocity is marked as an error sample, and the error sample is removed from the initial P-wave velocity sample set to obtain an effective P-wave velocity sample set.

[0020] Optionally, determining the blasting P-wave velocity based on the effective P-wave velocity sample set includes:

[0021] The minimum and maximum source-detection distances are determined based on the source-detection distances corresponding to each effective P-wave velocity in the effective P-wave velocity sample set.

[0022] Based on the minimum source-receiver distance, the maximum source-receiver distance, and the source-receiver distance corresponding to the effective P-wave velocity, determine the source-receiver distance normalization coefficient for each effective P-wave velocity;

[0023] The ratio of each effective P-wave velocity to the source-detector distance normalization coefficient is calculated, and the average value of the ratios is used to obtain the blasting P-wave velocity.

[0024] Optionally, before removing error samples from the initial P-wave velocity sample set according to the Laida criterion to obtain an effective P-wave velocity sample set, and determining the blasting P-wave velocity based on the effective P-wave velocity sample set, the method further includes:

[0025] The blasting was repeated multiple times, and after each blast, the calculated initial P-wave velocity was added to the initial P-wave velocity sample set.

[0026] According to another aspect of this application, a coal mine blasting longitudinal wave velocity measuring device is provided, the device comprising:

[0027] The information acquisition module is used to acquire the first arrival time of the blast wave to the source detector and the second arrival time to the target detector after the blasting of a preset explosive in a coal mine roadway. The target detectors are multiple.

[0028] The data calculation module is used to calculate the distance between each target detector and the source detector respectively, to obtain the source-detector distance corresponding to each target detector, and to calculate the blast wave travel time corresponding to each target detector based on the first arrival time and the second arrival time;

[0029] The set construction module is used to calculate the initial P-wave velocity corresponding to each target detector based on the blast wave travel time and the source-detector distance, and construct an initial P-wave velocity sample set containing the initial P-wave velocity corresponding to each target detector.

[0030] The wave velocity determination module is used to remove error samples from the initial P-wave velocity sample set according to the Laida criterion, obtain an effective P-wave velocity sample set, and determine the blasting P-wave velocity based on the effective P-wave velocity sample set.

[0031] Optionally, the data calculation module is further configured to:

[0032] The source detector is horizontally arranged on the coal wall in the coal mine roadway and is at a preset distance from the blasting position of the preset explosive.

[0033] The target detector is vertically arranged on the roof and / or floor of the coal mine roadway and within a preset range corresponding to the blasting position;

[0034] Obtain the three-dimensional coordinates of the source detector and the target detector;

[0035] Based on the three-dimensional coordinates of the source detector and the three-dimensional coordinates of the target detector, the source-detector distance between each target detector and the source detector is calculated.

[0036] Optionally, the data calculation module is further configured to:

[0037] The difference between the second duration and the first duration corresponding to each of the target detectors is calculated to obtain the blast wave travel time corresponding to each of the target detectors.

[0038] Optionally, the collection building module is further configured to:

[0039] The ratio of the source-detector distance to the travel time of the blast wave is calculated for each target detector to obtain the initial longitudinal wave velocity for each target detector.

[0040] Optionally, the collection building module is further configured to:

[0041] Calculate the arithmetic mean, standard deviation, and absolute difference between each initial P-wave velocity and the arithmetic mean for all initial P-wave velocities in the initial P-wave velocity sample set.

[0042] If the ratio of the absolute difference to the standard deviation is greater than or equal to a preset multiple, the corresponding initial P-wave velocity is marked as an error sample, and the error sample is removed from the initial P-wave velocity sample set to obtain an effective P-wave velocity sample set.

[0043] Optionally, the wave velocity determination module is further configured to:

[0044] The minimum and maximum source-detection distances are determined based on the source-detection distances corresponding to each effective P-wave velocity in the effective P-wave velocity sample set.

[0045] Based on the minimum source-receiver distance, the maximum source-receiver distance, and the source-receiver distance corresponding to the effective P-wave velocity, determine the source-receiver distance normalization coefficient for each effective P-wave velocity;

[0046] The ratio of each effective P-wave velocity to the source-detector distance normalization coefficient is calculated, and the average value of the ratios is used to obtain the blasting P-wave velocity.

[0047] Optionally, the collection building module is further configured to:

[0048] The blasting was repeated multiple times, and after each blast, the calculated initial P-wave velocity was added to the initial P-wave velocity sample set.

[0049] According to another aspect of this application, a storage medium is provided that stores a computer program thereon, which, when executed by a processor, implements the above-described method for measuring longitudinal wave velocity in coal mine blasting.

[0050] According to another aspect of this application, a computer device is provided, including a storage medium, a processor, and a computer program stored on the storage medium and executable on the processor, wherein the processor executes the program to implement the above-described method for measuring longitudinal wave velocity in coal mine blasting.

[0051] By employing the above technical solution, this application provides a method and apparatus for measuring the longitudinal wave velocity of coal mine blasting, a storage medium, and a computer device. After the blasting of a pre-set explosive in a coal mine roadway, the method acquires the first arrival time of the blasting longitudinal wave at the source detector and the second arrival time at the target detector. Multiple target detectors are used. The distance between each target detector and the source detector is calculated to obtain the source-detector distance for each target detector. The blasting longitudinal wave travel time for each target detector is calculated based on the first and second arrival times. The initial longitudinal wave velocity for each target detector is calculated based on the blasting longitudinal wave travel time and the source-detector distance. An initial longitudinal wave velocity sample set containing the initial longitudinal wave velocity for each target detector is constructed. Error samples in the initial longitudinal wave velocity sample set are removed according to the Laida criterion to obtain an effective longitudinal wave velocity sample set. The blasting longitudinal wave velocity is determined based on the effective longitudinal wave velocity sample set, thus realizing the measurement of the longitudinal wave velocity of coal mine blasting.

[0052] The above description is only an overview of the technical solution of this application. In order to better understand the technical means of this application and to implement it in accordance with the contents of the specification, and to make the above and other objects, features and advantages of this application more obvious and understandable, the following are specific embodiments of this application. Attached Figure Description

[0053] The accompanying drawings, which are included to provide a further understanding of this application and form part of this application, illustrate exemplary embodiments and are used to explain this application, but do not constitute an undue limitation of this application. In the drawings:

[0054] Figure 1 A schematic flowchart of a method for measuring longitudinal wave velocity in coal mine blasting, provided in an embodiment of this application, is shown.

[0055] Figure 2 This illustration shows a schematic diagram of the installation positions of the source detector and the target detector provided in an embodiment of this application;

[0056] Figure 3 A flowchart illustrating another method for measuring longitudinal wave velocity in coal mine blasting provided in an embodiment of this application is shown.

[0057] Figure 4 A schematic diagram of the source detector installation provided in an embodiment of this application is shown;

[0058] Figure 5 A schematic diagram of the structure of a coal mine blasting longitudinal wave velocity measuring device provided in an embodiment of this application is shown.

[0059] Among them, 21-target detector, 22-coal mine roadway, 23-source detector, 24-blasting location, 25-main coal mine roadway, 46-blast hole charging section, 47-blast hole sealing section, 48-anchor bolt, 49-source detector, 410-coal mine roadway, 411-coal wall. Detailed Implementation

[0060] The present application will be described in detail below with reference to the accompanying drawings and embodiments. It should be noted that, unless otherwise specified, the embodiments and features described in the embodiments of the present application can be combined with each other.

[0061] This embodiment provides a method for measuring longitudinal wave velocity in coal mine blasting, such as... Figure 1 As shown, the method includes:

[0062] Step 101: After the pre-set explosive material in the coal mine roadway is detonated, the first arrival time of the blast wave to the source detector and the second arrival time to the target detector are obtained, wherein there are multiple target detectors.

[0063] In the above embodiments of this application, firstly, the user installs a preset explosive in the aforementioned coal mine roadway, and then ignites the preset explosive. After the preset explosive is ignited and detonated, it will generate a blast wave. By obtaining the relevant parameters of the blast wave, the wave velocity of the coal mine blast wave can be measured.

[0064] In one specific embodiment, before measuring the wave velocity of a coal mine blasting wave, the coal mine where the aforementioned coal mine roadway is located needs to cease production to avoid noise interference. Simultaneously, the target detector and source detector are tuned to the same seismic monitoring system, and the clocks of each monitoring channel are precisely synchronized. Specifically, the method for installing the pre-set explosive is as follows: a blast hole is excavated in the aforementioned coal mine roadway; then, a pre-set explosive is installed in the blast hole, and the blast hole is sealed. Subsequently, the source detector is installed 0.5m from the blast hole opening, i.e., the blasting location, and multiple target detectors are simultaneously installed at different distances around the blast hole. Figure 2As shown, the aforementioned target detector 21 can be located in the coal mine roadway 22 and the main coal mine roadway 25. The aforementioned target detector 21 is arranged around the blasting location 24 at a certain distance, and is arranged in a combination of near and far distances. The aforementioned source detector 23 and target detector 21 can be passive vibration sensors. Then, the aforementioned preset explosive is ignited, and the arrival time of the blast wave caused by the preset explosive is obtained by professional software, which shows the first arrival time of the blast wave to the aforementioned source detector 23 and the second arrival time of the blast wave to the aforementioned target detector 21.

[0065] Step 102: Calculate the distance between each target detector and the source detector to obtain the source-detector distance corresponding to each target detector, and calculate the explosion wave travel time corresponding to each target detector based on the first arrival time and the second arrival time.

[0066] Next, the distance between each of the aforementioned target detectors and the aforementioned source detectors is calculated to obtain the source-detector distance corresponding to each of the aforementioned target detectors, i.e., the propagation distance of the aforementioned blast wave. Based on the aforementioned first arrival time and the aforementioned second arrival time, the blast travel time corresponding to each of the aforementioned target detectors is calculated. The aforementioned first arrival time is the time when the aforementioned blast wave arrives at the aforementioned source detector, i.e., the blast time. Based on the gap between the second arrival time of the aforementioned target detector and the aforementioned blast time, the waveform travel time of the aforementioned blast wave can be obtained. Therefore, after the aforementioned preset explosive material detonates, based on the distance between the aforementioned target detectors and the aforementioned source detectors, i.e., the propagation distance of the aforementioned blast wave, and the waveform travel time of the aforementioned blast wave, subsequent calculations and measurements of the blast longitudinal wave velocity can be performed.

[0067] Step 103: Calculate the initial P-wave velocity corresponding to each target detector based on the explosion wave travel time and the source-detector distance, and construct an initial P-wave velocity sample set containing the initial P-wave velocity corresponding to each target detector.

[0068] Next, based on the blast wave travel time corresponding to each of the aforementioned target detectors and the aforementioned source-detector distance, the initial P-wave velocity corresponding to each of the aforementioned target detectors is calculated, and then an initial P-wave velocity sample set containing the initial P-wave velocity corresponding to each of the aforementioned target detectors is constructed.

[0069] Step 104: Eliminate error samples from the initial P-wave velocity sample set according to the Laida criterion to obtain an effective P-wave velocity sample set, and determine the blasting P-wave velocity based on the effective P-wave velocity sample set.

[0070] Finally, error samples in the initial P-wave velocity sample set were removed according to the Laida criterion to obtain an effective P-wave velocity sample set. Based on this effective P-wave velocity sample set, the blasting P-wave velocity was determined. Therefore, by obtaining accurate and reasonable coal mine blasting P-wave velocities as fundamental parameters for coal mine microseismic source location, the accuracy of source location and the quality of coal mine microseismic monitoring can be improved, providing data support for safe coal mining.

[0071] In the above embodiments of this application, as the coal mining location changes, the geology and mining conditions will also change accordingly, so the blasting wave velocity will also be different. Therefore, when measuring the longitudinal wave velocity of coal mine blasting at different coal mining locations, it is necessary to remeasure the blasting longitudinal wave velocity using the above method.

[0072] By applying the technical solution of this embodiment, after the blasting of a pre-set explosive in a coal mine roadway, the first arrival time of the blast wave to the source detector and the second arrival time to the target detector are obtained. Multiple target detectors are used, and the distance between each target detector and the source detector is calculated to obtain the source-detector distance for each target detector. The blast wave travel time for each target detector is calculated based on the first and second arrival times. The initial P-wave velocity for each target detector is calculated based on the blast wave travel time and the source-detector distance, constructing an initial P-wave velocity sample set containing the initial P-wave velocity for each target detector. Error samples in the initial P-wave velocity sample set are removed according to the Laida criterion to obtain an effective P-wave velocity sample set. The blasting P-wave velocity is determined based on the effective P-wave velocity sample set, thus realizing the measurement of the blasting P-wave velocity in coal mines. By using the aforementioned blasting P-wave velocity as the basic parameter for locating the microseismic source in coal mines, the accuracy of source location and the quality of microseismic monitoring in coal mines are improved, providing data support for safe coal mining.

[0073] Furthermore, as a refinement and extension of the specific implementation of the above embodiments, and to fully illustrate the specific implementation process of this embodiment, another method for measuring longitudinal wave velocity in coal mine blasting is provided, such as... Figure 3 As shown, the method includes:

[0074] Step 301: After the pre-set explosive material in the coal mine roadway is detonated, the first arrival time of the blast wave to the source detector and the second arrival time to the target detector are obtained, wherein there are multiple target detectors.

[0075] Step 302: Obtain the three-dimensional coordinates of the source detector and the target detector, and calculate the source-detector distance between each target detector and the source detector based on the three-dimensional coordinates of the source detector and the target detector.

[0076] The source detector is horizontally arranged on the coal wall in the coal mine roadway and at a predetermined distance from the blasting position of the predetermined explosive, while the target detector is vertically arranged on the roof and / or floor of the coal mine roadway and within a predetermined range corresponding to the blasting position.

[0077] In the above embodiments of this application, after the preset explosive in the coal mine roadway is detonated, the first arrival time of the blast wave to the source detector and the second arrival time to the target detector are obtained, wherein there are multiple target detectors.

[0078] The aforementioned source detector is horizontally arranged on the coal wall in the aforementioned coal mine roadway and is at a predetermined distance from the blasting position of the aforementioned preset explosive, for example, 0.5m away. When installing the aforementioned source detector, the anchor bolt of the aforementioned source detector needs to be fully anchored, for example, the anchor bolt length can be greater than 2m. The aforementioned target detector can be installed on the top or bottom plate of the coal seam and within the predetermined range corresponding to the aforementioned blasting position. When installing the aforementioned target detector, the aforementioned target detector needs to be kept in a vertical state.

[0079] Specifically, the aforementioned source detector installation diagram is as follows: Figure 4 As shown, the aforementioned pre-set explosive is... Figure 4 The borehole charging section 46 is constructed with the aforementioned source detector 49 mounted on a horizontal anchor bolt 48 and positioned on the coal wall 411 of the coal mine roadway 410. The borehole direction is parallel to the anchor bolt direction. When the borehole is being charged, the length of the borehole sealing section 47 is greater than the length of the anchor bolt 48 to prevent the anchor bolt 48 from collapsing when the aforementioned pre-set explosives detonate.

[0080] Next, the three-dimensional coordinates of the source detector and the target detector are obtained. Based on the aforementioned three-dimensional coordinates of the source detector and the target detector, the source-detector distance between each of the aforementioned target detectors and the source detector is calculated.

[0081] Step 303: Calculate the difference between the second duration and the first duration for each target detector to obtain the blast wave travel time for each target detector.

[0082] Next, the difference between the second duration and the first duration corresponding to each of the aforementioned target detectors is calculated to obtain the blast wave travel time corresponding to each of the aforementioned target detectors.

[0083] Step 304: Calculate the ratio of the source-detector distance to the travel time of the blast wave for each target detector to obtain the initial P-wave velocity for each target detector, and construct an initial P-wave velocity sample set containing the initial P-wave velocity for each target detector.

[0084] Next, the ratio of the source-receiver distance to the travel time of the blast wave is calculated for each of the aforementioned target detectors to obtain the initial P-wave velocity for each of the aforementioned target detectors. Then, an initial P-wave velocity sample set containing the initial P-wave velocities for each of the aforementioned target detectors is constructed.

[0085] Step 305: Repeat the blasting multiple times, and after each blast, add the calculated initial P-wave velocity to the initial P-wave velocity sample set.

[0086] In the above embodiments of this application, the blasting is repeated multiple times, and the calculated initial P-wave velocity is added to the initial P-wave velocity sample set after each blast.

[0087] Step 306: Calculate the arithmetic mean, standard deviation, and absolute difference between each initial P-wave velocity and the arithmetic mean of all initial P-wave velocities in the initial P-wave velocity sample set. If the ratio of the absolute difference to the standard deviation is greater than or equal to a preset multiple, the corresponding initial P-wave velocity is marked as an error sample, and the error sample is removed from the initial P-wave velocity sample set to obtain a valid P-wave velocity sample set.

[0088] Next, the arithmetic mean, standard deviation, and absolute difference between each initial P-wave velocity and the arithmetic mean are calculated for all the initial P-wave velocities in the aforementioned initial P-wave velocity sample set. If the ratio of the absolute difference to the standard deviation is greater than or equal to a preset multiple, the corresponding initial P-wave velocity is marked as an error sample and removed from the aforementioned initial P-wave velocity sample set to obtain a valid P-wave velocity sample set.

[0089] Step 307: Based on the source-receiver distance corresponding to each effective P-wave velocity in the effective P-wave velocity sample set, determine the minimum and maximum source-receiver distances. Based on the minimum, maximum, and source-receiver distances corresponding to the effective P-wave velocities, determine the source-receiver distance normalization coefficient for each effective P-wave velocity. Calculate the ratio of each effective P-wave velocity to the source-receiver distance normalization coefficient, and calculate the average value of the ratios to obtain the blasting P-wave velocity.

[0090] Finally, based on the source-receiver distance corresponding to each effective P-wave velocity in the aforementioned effective P-wave velocity sample set, the minimum and maximum source-receiver distances are determined. Based on the aforementioned minimum, maximum, and source-receiver distances corresponding to the aforementioned effective P-wave velocities, the source-receiver distance normalization coefficient for each of the aforementioned effective P-wave velocities is determined. The ratio of each of the aforementioned effective P-wave velocities to the aforementioned source-receiver distance normalization coefficient is calculated, and the average value of the aforementioned ratios is calculated to obtain the blasting P-wave velocity.

[0091] In one specific embodiment, for example: assuming there are m target detectors, the aforementioned blasting test is repeated n times at different locations within the installation range of the m target detectors, where mn > 50. This ensures a sufficiently large sample size to guarantee the accuracy of the blasting longitudinal wave velocity measurement results. Specifically:

[0092] After the pre-set explosives are detonated in the coal mine roadway, step 301 obtains the first arrival time of the blast wave to the source detector and the second arrival time to the target detector, thus obtaining the aforementioned first arrival times of n source detectors, where the first arrival time of the source detector for the i-th blast is denoted as t. 0i (1≤i≤n), the aforementioned second arrival time of nm target detectors is obtained, where the second arrival time of the i-th blast on the j-th target detector is denoted as t. ij (1≤i≤n,1≤j≤m).

[0093] Assume the three-dimensional coordinates of the source detector for the i-th blast are (x... 0i y 0i , z 0i The three-dimensional coordinates of the detector of the j-th target in the i-th blast are (x, y, y). ij y ij , z ij ), where 1≤i≤n, 1≤j≤m.

[0094] The source-detector distance S between the detector and the source detector for the j-th target during the i-th blast in step 302 is obtained. ij The formula for calculating the three-dimensional distance between the target detector and the source detector is as follows:

[0095]

[0096] The blast wave travel time T corresponding to the detector of the j-th target during the i-th blast in step 303 is used. ij The calculation formula is as follows:

[0097] T ij =t ij -t 0i

[0098] The travel time T of the blast wave is calculated in step 304. ij Source-detection distance S ij The ratio is used to obtain mn initial P-wave velocity samples. An initial P-wave velocity sample set containing the aforementioned mn initial P-wave velocity samples is constructed, where the initial P-wave velocity of the j-th target detector in the i-th blast is represented as v. ij (1≤i≤n,1≤j≤m).

[0099] The initial longitudinal wave velocity V of the detector of the j-th target during the i-th blast. ijThe calculation formula is as follows:

[0100]

[0101] Step 306 calculates all the aforementioned initial P-wave velocities V in the aforementioned initial P-wave velocity sample set. ij arithmetic mean for:

[0102]

[0103] Calculate all the aforementioned initial P-wave velocities V in the aforementioned initial P-wave velocity sample set. ij The standard deviation M is:

[0104]

[0105] like Then it is assumed that the aforementioned initial longitudinal wave velocity V ij Samples with large errors are considered coarse error samples and are deleted; otherwise, they are retained. After deleting the aforementioned coarse error samples, the effective P-wave velocity sample set is obtained.

[0106] The source-detection distance normalization coefficient g is calculated in step 307. h Specifically:

[0107] Assuming that the effective P-wave velocity sample set obtained through step 306 contains k samples, and the h-th sample is denoted as v h The source-detection distance corresponding to the h-th sample is S. h The source-detector distances corresponding to all samples form a set Q, where h≤k≤mn.

[0108] Then, the source-receiver distance normalization coefficient is used as a weighting factor to calculate the weighted average of the effective P-wave velocity samples, which is the final required blasting P-wave velocity. The calculation method is as follows:

[0109] The source-detection distance normalization coefficient corresponding to the h-th sample Where MinQ is the minimum value of the elements in set Q, and MaxQ is the maximum value of the elements in set Q;

[0110] Final blast longitudinal wave velocity

[0111] Furthermore, as Figure 1 In specific implementation of the method, this application provides a coal mine blasting longitudinal wave velocity measuring device, such as... Figure 5 As shown, the device includes:

[0112] The information acquisition module 51 is used to acquire the first arrival time of the blast wave to the source detector and the second arrival time to the target detector after the blasting of the preset explosive in the coal mine roadway. The target detector is multiple.

[0113] The data calculation module 52 is used to calculate the distance between each target detector and the source detector respectively, to obtain the source-detector distance corresponding to each target detector, and to calculate the blast wave travel time corresponding to each target detector based on the first arrival time and the second arrival time;

[0114] The set construction module 53 is used to calculate the initial P-wave velocity corresponding to each target detector based on the explosion wave travel time corresponding to each target detector and the source-detector distance, and construct an initial P-wave velocity sample set containing the initial P-wave velocity corresponding to each target detector.

[0115] The wave velocity determination module 54 is used to remove error samples from the initial P-wave velocity sample set according to the Laida criterion, obtain an effective P-wave velocity sample set, and determine the blasting P-wave velocity based on the effective P-wave velocity sample set.

[0116] Optionally, the data calculation module 52 is further configured to:

[0117] The source detector is horizontally arranged on the coal wall in the coal mine roadway and is at a preset distance from the blasting position of the preset explosive.

[0118] The target detector is vertically arranged on the roof and / or floor of the coal mine roadway and within a preset range corresponding to the blasting position;

[0119] Obtain the three-dimensional coordinates of the source detector and the target detector;

[0120] Based on the three-dimensional coordinates of the source detector and the three-dimensional coordinates of the target detector, the source-detector distance between each target detector and the source detector is calculated.

[0121] Optionally, the data calculation module 52 is further configured to:

[0122] The difference between the second duration and the first duration corresponding to each of the target detectors is calculated to obtain the blast wave travel time corresponding to each of the target detectors.

[0123] Optionally, the set construction module 53 is further configured to:

[0124] The ratio of the source-detector distance to the travel time of the blast wave is calculated for each target detector to obtain the initial longitudinal wave velocity for each target detector.

[0125] Optionally, the set construction module 53 is further configured to:

[0126] Calculate the arithmetic mean, standard deviation, and absolute difference between each initial P-wave velocity and the arithmetic mean for all initial P-wave velocities in the initial P-wave velocity sample set.

[0127] If the ratio of the absolute difference to the standard deviation is greater than or equal to a preset multiple, the corresponding initial P-wave velocity is marked as an error sample, and the error sample is removed from the initial P-wave velocity sample set to obtain an effective P-wave velocity sample set.

[0128] Optionally, the wave velocity determination module 54 is further configured to:

[0129] The minimum and maximum source-detection distances are determined based on the source-detection distances corresponding to each effective P-wave velocity in the effective P-wave velocity sample set.

[0130] Based on the minimum source-receiver distance, the maximum source-receiver distance, and the source-receiver distance corresponding to the effective P-wave velocity, determine the source-receiver distance normalization coefficient for each effective P-wave velocity;

[0131] The ratio of each effective P-wave velocity to the source-detector distance normalization coefficient is calculated, and the average value of the ratios is used to obtain the blasting P-wave velocity.

[0132] Optionally, the set construction module 53 is further configured to:

[0133] The blasting was repeated multiple times, and after each blast, the calculated initial P-wave velocity was added to the initial P-wave velocity sample set.

[0134] It should be noted that other corresponding descriptions of the functional units involved in the coal mine blasting longitudinal wave velocity measuring device provided in this application embodiment can be found in the following references. Figure 1 and Figure 3 The corresponding descriptions in the method will not be repeated here.

[0135] Based on the above, Figure 1 and Figure 3 Accordingly, this application also provides a storage medium storing a computer program, which, when executed by a processor, implements the above-described method. Figure 1 and Figure 3 The method for measuring longitudinal wave velocity in coal mine blasting is shown.

[0136] Based on this understanding, the technical solution of this application can be embodied in the form of a software product. This software product can be stored in a non-volatile storage medium (such as a CD-ROM, USB flash drive, or portable hard drive), and includes several instructions to cause a computer device (such as a personal computer, server, or network device) to execute the methods described in the various implementation scenarios of this application.

[0137] Based on the above, Figure 1 and Figure 3 The method shown, and Figure 5 To achieve the above objectives, the present application also provides a computer device, specifically a personal computer, server, network device, etc., as shown in the virtual device embodiment. This computer device includes a storage medium and a processor; the storage medium stores a computer program; the processor executes the computer program to achieve the above-described objectives. Figure 1 and Figure 3 The method for measuring longitudinal wave velocity in coal mine blasting is shown.

[0138] Optionally, the computer device may also include a user interface, a network interface, a camera, radio frequency (RF) circuitry, sensors, audio circuitry, a Wi-Fi module, etc. The user interface may include a display screen, input units such as a keyboard, etc., and optional user interfaces may also include USB interfaces, card reader interfaces, etc. The network interface may optionally include standard wired interfaces, wireless interfaces (such as Bluetooth interfaces, Wi-Fi interfaces), etc.

[0139] Those skilled in the art will understand that the computer device structure provided in this embodiment does not constitute a limitation on the computer device, and may include more or fewer components, or combine certain components, or have different component arrangements.

[0140] The storage medium may also include an operating system and a network communication module. The operating system is a program that manages and stores the hardware and software resources of a computer device, supporting the operation of information processing programs and other software and / or programs. The network communication module is used to enable communication between the various components within the storage medium, as well as communication with other hardware and software within the physical device.

[0141] Through the above description of the embodiments, those skilled in the art can clearly understand that this application can be implemented by means of software plus necessary general-purpose hardware platform, or by hardware. After the blasting of a preset explosive in a coal mine roadway, the first arrival time of the blast wave to the source detector and the second arrival time to the target detector are obtained. There are multiple target detectors. The distance between each target detector and the source detector is calculated to obtain the source-detector distance corresponding to each target detector. The blast wave travel time corresponding to each target detector is calculated based on the first and second arrival times. The initial P-wave velocity corresponding to each target detector is calculated based on the blast wave travel time and the source-detector distance. An initial P-wave velocity sample set containing the initial P-wave velocity corresponding to each target detector is constructed. Error samples in the initial P-wave velocity sample set are removed according to the Laida criterion to obtain an effective P-wave velocity sample set. The blasting P-wave velocity is determined based on the effective P-wave velocity sample set, thus realizing the measurement of the P-wave velocity of coal mine blasting.

[0142] Those skilled in the art will understand that the accompanying drawings are merely schematic diagrams of a preferred embodiment, and the modules or processes shown in the drawings are not necessarily essential for implementing this application. Those skilled in the art will understand that the modules in the apparatus of the embodiment can be distributed within the apparatus of the embodiment as described, or can be modified to be located in one or more apparatuses different from this embodiment. The modules of the above-described embodiment can be combined into one module, or further divided into multiple sub-modules.

[0143] The serial numbers in this application are for descriptive purposes only and do not represent the superiority or inferiority of any particular implementation scenario. The above disclosures are merely a few specific implementation scenarios of this application; however, this application is not limited thereto, and any variations conceived by those skilled in the art should fall within the protection scope of this application.

Claims

1. A method for measuring longitudinal wave velocity in coal mine blasting, characterized in that, The method includes: After the pre-set explosives in the coal mine roadway are detonated, the first arrival time of the blast wave to the source detector and the second arrival time to the target detector are obtained, wherein there are multiple target detectors. The distance between each target detector and the source detector is calculated to obtain the source-detector distance corresponding to each target detector. Based on the first arrival time and the second arrival time, the explosion wave travel time corresponding to each target detector is calculated. Based on the blast wave travel time corresponding to each target detector and the source-detector distance, the initial P-wave velocity corresponding to each target detector is calculated, and an initial P-wave velocity sample set containing the initial P-wave velocity corresponding to each target detector is constructed. Repeat the blasting process multiple times, and after each blast, add the calculated initial P-wave velocity to the initial P-wave velocity sample set; The error samples in the initial P-wave velocity sample set are removed according to the Laida criterion to obtain the effective P-wave velocity sample set. The minimum and maximum source-detection distances are determined based on the source-detection distances corresponding to each effective P-wave velocity in the effective P-wave velocity sample set. Based on the minimum source-receiver distance, the maximum source-receiver distance, and the source-receiver distance corresponding to the effective P-wave velocity, determine the source-receiver distance normalization coefficient for each effective P-wave velocity; The ratio of each effective P-wave velocity to the source-detector distance normalization coefficient is calculated, and the average value of the ratios is used to obtain the blasting P-wave velocity.

2. The method according to claim 1, characterized in that, The source detector is horizontally arranged on the coal wall in the coal mine roadway and at a predetermined distance from the blasting position of the predetermined explosive; the target detector is vertically arranged on the roof and / or floor of the coal mine roadway and within a predetermined range corresponding to the blasting position; the step of calculating the distance between each target detector and the source detector to obtain the source-detector distance corresponding to each target detector includes: Obtain the three-dimensional coordinates of the source detector and the target detector; Based on the three-dimensional coordinates of the source detector and the three-dimensional coordinates of the target detector, the source-detector distance between each target detector and the source detector is calculated.

3. The method according to claim 1, characterized in that, The step of calculating the explosion wave travel time corresponding to each target detector based on the first arrival time and the second arrival time includes: The difference between the second arrival time and the first arrival time for each target detector is calculated to obtain the blast wave travel time for each target detector.

4. The method according to claim 1, characterized in that, The step of calculating the initial P-wave velocity for each target detector based on the explosion wave travel time and the source-detector distance includes: The ratio of the source-detector distance to the travel time of the blast wave is calculated for each target detector to obtain the initial longitudinal wave velocity for each target detector.

5. The method according to claim 1, characterized in that, The step of removing error samples from the initial P-wave velocity sample set according to the Laida criterion to obtain an effective P-wave velocity sample set includes: Calculate the arithmetic mean, standard deviation, and absolute difference between each initial P-wave velocity and the arithmetic mean for all initial P-wave velocities in the initial P-wave velocity sample set. If the ratio of the absolute difference to the standard deviation is greater than or equal to a preset multiple, the corresponding initial P-wave velocity is marked as an error sample, and the error sample is removed from the initial P-wave velocity sample set to obtain an effective P-wave velocity sample set.

6. A device for measuring longitudinal wave velocity in coal mine blasting, characterized in that, The device includes: The information acquisition module is used to acquire the first arrival time of the blast wave to the source detector and the second arrival time to the target detector after the blasting of a preset explosive in a coal mine roadway. The target detectors are multiple. The data calculation module is used to calculate the distance between each target detector and the source detector to obtain the source-detector distance corresponding to each target detector, and to calculate the blast wave travel time corresponding to each target detector based on the first arrival time and the second arrival time. The set construction module is used to calculate the initial P-wave velocity corresponding to each target detector based on the blast wave travel time corresponding to each target detector and the source-detector distance, and construct an initial P-wave velocity sample set containing the initial P-wave velocity corresponding to each target detector. The wave velocity determination module is used for repeated blasting, and after each blast, it adds the calculated initial P-wave velocity to the initial P-wave velocity sample set; it removes error samples from the initial P-wave velocity sample set according to the Laida criterion to obtain the effective P-wave velocity sample set; it determines the minimum and maximum source-receiver distances based on the source-receiver distances corresponding to each effective P-wave velocity in the effective P-wave velocity sample set; it determines the source-receiver distance normalization coefficient for each effective P-wave velocity based on the minimum, maximum, and source-receiver distances corresponding to the effective P-wave velocities; it calculates the ratio of each effective P-wave velocity to the source-receiver distance normalization coefficient, and calculates the average of the ratios to obtain the blasting P-wave velocity.

7. A storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it implements the method for measuring longitudinal wave velocity in coal mine blasting as described in any one of claims 1 to 5.

8. A computer device, comprising a storage medium, a processor, and a computer program stored on the storage medium and executable on the processor, characterized in that, When the processor executes the computer program, it implements the method for measuring longitudinal wave velocity in coal mine blasting as described in any one of claims 1 to 5.

Images