Rock blasting hole filling quality assessment method based on infrared thermal imaging technology

Abstract

The present invention discloses a method for evaluating the filling quality of rock blasting boreholes based on infrared thermal imaging technology, comprising the steps of determining monitoring points, setting up instruments, capturing thermal images, processing video data, determining and calculating a filling quality evaluation index, and evaluating the filling quality. The present invention has the beneficial effect of dynamically displaying the explosive products formed by the high-temperature, high-pressure gas released by explosives at the moment of blasting on a display through an infrared thermal imager, thereby compensating for the movement process of the explosive products rushing out of the free surface that cannot be displayed by visible light camera equipment based on naked eye observation. Furthermore, by mapping the thermal image into a three-dimensional volume calculation, the total volume of the explosive products of the mixed gas rushing out is obtained, thereby quantifying the total amount of the explosive products that escaped. Furthermore, by creatively proposing a quantitative index ξ for judging the quality of filling, the present invention provides a clear target orientation for filling optimization and a basic guarantee for the scientific and rational optimization method.

Description

Technical Field

[0001] The invention belongs to engineering blasting technology, in particular to a rock blasting hole filling quality assessment method based on infrared thermal imaging technology. Background Art

[0002] In blasting projects in fields such as transportation, water conservancy, mining, and tunnel blasting, the quality of blasthole filling significantly impacts the blasting results. Good filling quality not only ensures the full reaction of the explosive within the blasthole, maximizing the conversion of the explosive's chemical energy into mechanical energy, but also reduces the temperature and pressure of the explosive gas, alters the stress pattern in the rock, and improves the thermal efficiency of the explosive, converting more thermal energy into mechanical work. This increases the utilization rate of the explosive energy and ultimately improves the blasting results. However, in actual construction, poor filling quality is unavoidable due to various reasons. This causes the energy of the explosive within the blasthole to rush out along the hole wall after detonation. This premature energy escape creates a lumen in the filling section at the top of the blasthole, depriving the surrounding rock of energy, preventing effective rock fragmentation, and resulting in the formation of oversized fragments on the surface after blasting, increasing secondary processing costs. Therefore, it is necessary to quantitatively analyze the energy diffusion state of explosives during blasting to assess filling quality and provide a reference for blasting optimization. However, the existing technology lacks quantitative research on the impact of the filling section on the blasting effect, making it difficult to take targeted improvement measures. Summary of the Invention

[0003] The purpose of the present invention is to address the problem that in the existing blasting engineering field, there is a lack of quantitative research on the impact of filling quality on blasting effect, which results in the improvement measures being not targeted enough. The present invention provides a rock blasting borehole filling quality assessment method based on infrared thermal imaging technology. By introducing infrared thermal imaging technology into engineering blasting to detect borehole filling quality, the method can effectively capture and present the free surface diffusion trajectory of the explosive products escaping when the explosive explodes in the borehole, and display the size of the gas cloud in different colors through post-processing. By comparing the size of the temperature color block, the filling quality can be identified. This provides technicians with a new method and means for scientifically analyzing blasting quality, and also provides a basis for targeted improvements to improve filling quality.

[0004] To achieve the aforementioned objectives, the present invention adopts the following technical solutions.

[0005] A rock blasting hole filling quality assessment method based on infrared thermal imaging technology comprises the following steps:

[0006] S1, Determination of monitoring points: Based on the blasting design parameters and surrounding environmental data of the blasting area, determine the photography boundary and boundary control points of the blasting area, and measure the position coordinates of the control points;

[0007] S2, instrument installation: calculating the location of the instrument installation point based on the parameter characteristics of the thermal imager and the position coordinates of the boundary control point, installing the thermal imager according to the calculation result, and completing the shooting parameter setting;

[0008] S3, thermal imaging: during the detonation process in the target blasting area, infrared thermal imaging video is obtained by using the thermal imager;

[0009] S4, video data processing: first use analysis software to decompose the infrared thermal image video into i The corresponding single-frame picture; then, the color cloud image that is different from the ambient temperature in the single-frame picture is mapped into a three-dimensional body to obtain the time t i Volume V of the gas product of the explosion i Finally, after screening according to the principle of non-repetitive calculation, the total volume of explosive gas escaping from the blasthole during the blasting process is obtained by summing up V 红 ;

[0010] S5. Determination and calculation of the filling quality evaluation index: The filling quality evaluation index is defined as ξ, which is calculated as follows:

[0011]

[0012] Where V 红 is the total volume of infrared thermal imaging of explosive gas products escaping from the free surface of the blasthole top; V is the theoretical volume of gas released by the detonation reaction of the total amount of explosives in the blasting area;

[0013] S6, filling quality evaluation: according to the pre-determined evaluation index ξ cr The evaluation index ξ is evaluated. If it does not exceed the index value, it is qualified; if it exceeds the index value, the filling quality is improved.

[0014] The method of the present invention adopts the above-mentioned scheme, and clearly displays the movement process of the explosive products formed by the high-temperature and high-pressure gas released by the explosives at the moment of blasting from the free surface at the top of the blasthole on the display through the infrared thermal imager, making up for the movement process of the explosive products rushing out of the free surface that cannot be displayed by the visible light camera equipment based on naked eye observation, and by calculating the volume mapped by the thermal image at different temperatures, the total volume of the mixed gas explosive products rushing out of the free surface is obtained, and the total amount of the explosive products escaping from the free surface is quantified; in addition, the method also creatively proposes a quantitative index ξ for judging the quality of filling, which provides a clear target direction for filling optimization and also provides a basic guarantee for the scientific and rational optimization method. Among them, the evaluation index ξ crThe value range is usually 5% to 15%, determined by the type of blasting project. The method of the present invention can also make use of the thermal image distribution uniformity resolution feature of the analysis software to qualitatively judge the quality of the filling. If the filling quality is good, the cloud image distribution is uniform. Otherwise, the quality is poor and needs to be improved. The improvement is mainly in the filling material and filling method.

[0015] Among them, the evaluation index is proposed based on the following basic logic.

[0016] After detonation, explosives produce high-temperature, high-pressure explosive gas products. During this process, the explosive energy is converted into the internal energy and kinetic energy of the explosive gas products. The explosive gas then acts on the rock, causing it to break and throw. Some of these explosive products do not contribute to the rock-breaking process, but instead escape from the blocked area of ​​the blasthole, resulting in energy loss. Therefore, the quality of the blasthole packing is directly related to the amount of explosive gas that escapes within a given time after the blast. This can be characterized by the amount of explosive gas that escapes from the upper blocked section of the blasthole, that is, the volume of gas that escapes within a given time:

[0017] According to the basic theory in this field, the total energy released by the explosive gas products produced by the explosion of explosives can be calculated by the following formula:

[0018] Q 总 =C·m·T=C·ρ t ·V·T;

[0019] In the formula, T is the relative temperature of the explosive gas product at time t, m is the mass of the explosive gas product, C is the specific heat of the explosive gas product, ρ t is the density of the explosive gas products at time t, and V is the total volume of the explosive gas products produced by the explosive explosion. The calculation formula for the volume of the explosive gas products produced by the explosive explosion is:

[0020] V=M e v e ;

[0021] In the formula, M e is the total mass of explosives in the blasting area, V e The volume of explosive products produced by the complete reaction of unit mass of explosives is related to the composition of the explosives and is a constant under certain conditions of explosives.

[0022] At time t, the energy of the explosive gas lost due to escaping from the upper plugging section of the blasthole is:

[0023] Q ξ =C·ρ t ·V 红 T;

[0024] From this, we can calculate the percentage of energy lost from the blocked part of the blasthole at time t:

[0025]

[0026] Preferably, in step S2, the instrument installation point is determined according to the direction of the rock mass thrown away from the blasting area and is higher than the elevation H of the step top. top The horizontal distance L between the installation point and the center of the blasting area is greater than the horizontal distance R of the earthquake source when the maximum blasting vibration velocity is 1 cm / s. v , and is less than the object distance R at the longest focal length of the target burst area d ; Elevation of erection point H 红 Determined based on the installation point location and the instrument's field of view angle β. This ensures instrument safety and obtains clear infrared video data covering the explosion area.

[0027] Further preferably, the infrared thermal imager is installed on the top of the step, and the installation point height H 红 Calculate as follows:

[0028] H 红 =H top +L*tanα, where α is the angle between the corresponding side of the field of view of the infrared thermal imager and the horizontal plane of the target explosion area, where 0≤α<90°.

[0029] This fully considers the field of view of the instrument itself to meet the requirements of infrared video data shooting range and clarity.

[0030] Preferably, in step S4, the infrared thermal imaging video is decomposed into i The corresponding single-frame picture includes importing the infrared thermal imaging video stream data into the analysis software that matches the format, and using the analysis software to decompose the infrared thermal imaging video stream into single-frame pictures, and making each frame picture correspond to a time t i ;

[0031] The color cloud image that is different from the ambient temperature in a single frame image is mapped into a three-dimensional volume to obtain the temperature at the time t i Volume V of the gas product of the explosion i , including importing the color cloud images that are different from the ambient temperature and presented in different temperature gradients in all single-frame images into digital image processing software in jpg format; then combining the position coordinates of the boundary control points measured in step S1, using the polyline function of the drawing tool to draw a closed contour line along the edge of the color block, and based on the shooting angle, the explosion area range and the time segmentation based on the delay parameter, mapping the surface formed by the closed contour line of each segment into a three-dimensional body, and obtaining the time t at that time i Volume V corresponding to different temperature gradients i ; and calculate the volume sum of various temperature gradients that are different from the ambient temperature at a single moment in the corresponding partition;

[0032] After screening according to the principle of non-repetitive calculation, the total volume of explosive gas escaping from the blasthole during the blasting process in the blasting area is obtained by summing up. 红 , including removing the part of the explosive gas generated by the previous moment that has not yet dissipated naturally; the maximum principle screening means only selecting the maximum value V in the sum of the explosive gas volumes in each blasting period jmaz As the basis for total volume calculation;

[0033] Total volume of explosive gas V 红 Calculate as follows: In the formula, j is the specific time segment, and n is the number of segments which is a natural number.

[0034] This ensures accurate total volume results and evaluation results. The process of mapping the surface formed by the closed contour lines into a three-dimensional volume involves mapping the surface as the full cross-section of a solid of revolution in the case of a single-hole explosion. For example, the full cross-section of a sphere is a circle with the sphere's diameter as its diameter, and the full cross-section of a cylinder is a rectangle with the cylinder's length and diameter as its two adjacent sides. In the case of simultaneous multi-hole explosions, the amount of explosive product released at that moment can be determined by multiplying the single-hole mapping by the number of simultaneously exploded holes based on the shooting angle. When the resulting image is a composite multi-hole image, empirical data from single-hole distribution patterns can be combined to map it into a wavy three-dimensional volume. The most recent explosive product and the previous explosive product that has not yet naturally dissipated can be distinguished based on temperature differences reflected by color differences in the infrared image; alternatively, differences in diffusion rates can be used to distinguish them, thereby eliminating any remaining, undissipated, portions from the previous explosion.

[0035] Preferably, in step S5, V is calculated as follows:

[0036] V=M e v e ;

[0037] Where M e is the total mass of explosives in the blasting area, v e The volume of explosive products produced by the complete reaction of a unit mass of explosives is a constant based on the explosive composition. This is used to accurately determine the total volume of explosive products. The total explosive mass in the blasting area is determined based on the size of the blasting area, the arrangement of blastholes, the rock properties in the blasting area, and the blasting design criteria. Given these factors, the total explosive mass in the blasting area is the corresponding fixed quantity.

[0038] Preferably, in step S6, the principle for evaluating the quality of packing is that the smaller ξ is, the better the packing quality is; conversely, the worse the packing quality is, the more necessary it is to optimize the blasting design; wherein, ξ≤ξ cr , then qualified; ξ>ξ cr , then it is unqualified and the blasting design needs to be improved, where ξ crThe value range is 5% to 15%. By rationally determining evaluation indicators and process control, blasting efficiency can be improved and resources can be saved. The evaluation indicators can be adjusted appropriately based on the type of blasting project, combined with historical experience, the level of refinement of on-site management, and cost targets.

[0039] The beneficial effect of the present invention is that, through the infrared thermal imager, the movement process of the explosive products formed by the high-temperature and high-pressure gas released from the explosive at the moment of blasting, escaping from the free surface at the top of the blasthole, is clearly and dynamically displayed on the display, thereby compensating for the movement process of the explosive products rushing out of the free surface that cannot be displayed by visible light camera equipment based on naked eye observation. By calculating the volume mapped by the thermal image at different temperatures, the total volume of the mixed gas explosive products rushing out of the free surface is obtained, and the total amount of explosive products escaping from the free surface is quantified. In addition, the method also creatively proposes a quantitative index ξ for judging the quality of the packing, providing a clear target for packing optimization and a basic guarantee for the scientific and rational optimization method. The method can also use the thermal image distribution uniformity resolution characteristics of the analysis software to qualitatively judge the quality of the packing. If the packing quality is good, the cloud map distribution is uniform. Otherwise, the quality is poor and needs improvement. The improvement is mainly in the packing material and packing method. BRIEF DESCRIPTION OF THE DRAWINGS

[0040] Figure 1 It is a flow chart of the present invention.

[0041] Figure 2 It is a schematic diagram of the step surface parameters, instrument installation points and punching in the present invention.

[0042] Figure 3 It is an infrared thermal imaging diagram of step blasting in the present invention.

[0043] In the attached figure: 1-infrared thermal imager; 2-tripod; 3-control point; 4-contour of blast gas products; 5-blasting flying rocks; 6-open surface; 7-filling material; 8-explosive; 9-blast hole. Symbol H 红 - elevation of the instrument installation point; L - horizontal distance between the instrument installation point and the center of the blasting area; β - instrument field of view angle; α - angle between the corresponding side of the infrared thermal imager's field of view and the horizontal plane of the target blasting area; L 药 -Blasthole charge length; L 堵 -Blasthole filling length; H top -step surface elevation; W-line of least resistance; H-step height. DETAILED DESCRIPTION

[0044] The present invention will be further described below with reference to the accompanying drawings, but the present invention is not limited to the scope of the embodiments described.

[0045] See also Figure 1 、 Figure 2 and Figure 3 A rock blast hole filling quality assessment method based on infrared thermal imaging technology includes the following steps:

[0046] S1, Determination of monitoring points: Based on the blasting design parameters and surrounding environmental data of the blasting area, determine the photography boundary and boundary control points of the blasting area, and measure the position coordinates of the control points;

[0047] S2, instrument installation: Calculate the location of the instrument installation point based on the parameter characteristics of the thermal imager and the coordinates of the boundary control points, install the thermal imager according to the calculation results, and complete the shooting parameter setting;

[0048] S3, thermal imaging: during the detonation process in the target blasting area, infrared thermal imaging video is obtained by using the thermal imager;

[0049] S4, video data processing: first use analysis software to decompose the infrared thermal image video into i The corresponding single-frame picture; then, the color cloud image that is different from the ambient temperature in the single-frame picture is mapped into a three-dimensional body to obtain the time t i Volume V of the gas product of the explosion i Finally, after screening according to the principle of non-repetitive calculation, the total volume of explosive gas escaping from the blasthole during the blasting process is obtained by summing up V 红 ;

[0050] S5. Determination and calculation of the filling quality evaluation index: The filling quality evaluation index is defined as ξ, which is calculated as follows:

[0051]

[0052] Where V 红 is the total volume of infrared thermal imaging of explosive gas products escaping from the free surface of the top of the blasthole; V is the theoretical volume of gas released by the detonation reaction of the total amount of explosives in the blasting area.

[0053] S6, filling quality evaluation: according to the pre-determined evaluation index ξ cr The evaluation index ξ is evaluated. If it does not exceed the index value, it is qualified; if it exceeds the index value, the filling quality is improved.

[0054] In step S2, the instrument installation point is determined according to the direction of the rock mass thrown away from the blasting area and is higher than the elevation H of the top plate of the step. top The horizontal distance L between the installation point and the center of the blasting area is greater than the horizontal distance R of the earthquake source when the maximum blasting vibration velocity is 1 cm / s. v , and is less than the object distance R at the longest focal length of the target burst area d ; Elevation of erection point H 红Determined based on the installation point location and the instrument field of view angle β. The infrared thermal imager is installed on the top of the step, and the installation point elevation H 红 Calculate as follows:

[0055] H 红 =H top +L*tanα, where α is the angle between the corresponding side of the field of view of the infrared thermal imager and the horizontal plane of the target explosion area, where 0≤α<90°.

[0056] In step S4, the infrared thermal imaging video is decomposed into the time t i The corresponding single-frame picture includes importing the infrared thermal imaging video stream data into the analysis software that matches the format, and using the analysis software to decompose the infrared thermal imaging video stream into single-frame pictures, and making each frame picture correspond to a time t i ;

[0057] The color cloud image that is different from the ambient temperature in a single frame image is mapped into a three-dimensional volume to obtain the temperature at the time t i Volume V of the gas product of the explosion i , including importing the color cloud images that are different from the ambient temperature and presented in different temperature gradients in all single-frame images into digital image processing software in jpg format; then combining the position coordinates of the boundary control points measured in step S1, using the polyline function of the drawing tool to draw a closed contour line along the edge of the color block, and based on the shooting angle, the explosion area range and the time segmentation based on the delay parameter, mapping the surface formed by the closed contour line of each segment into a three-dimensional body, and obtaining the time t at that time i Volume V corresponding to different temperature gradients i ; and calculate the volume sum of various temperature gradients that are different from the ambient temperature at a single moment in the corresponding partition;

[0058] After screening according to the principle of non-repetitive calculation, the total volume of explosive gas escaping from the blasthole during the blasting process in the blasting area is obtained by summing up. 红 , including removing the part of the explosive gas generated by the previous moment that has not yet dissipated naturally; the maximum principle screening means only selecting the maximum value V in the sum of the explosive gas volumes in each blasting period jmax As the basis for total volume calculation;

[0059] Total volume of explosive gas V 红 Calculate as follows: In the formula, j is the specific time segment, and n is the number of segments which is a natural number.

[0060] In step S5, V is calculated as follows:

[0061] V=M e v e ;

[0062] Where M e is the total mass of explosives in the blasting area, v e It is the volume of explosive products produced by complete reaction of unit mass of explosive, which is a constant based on the composition of the explosive.

[0063] In step S6, the principle of evaluating the quality of packing is that the smaller ξ is, the better the packing quality is; conversely, the worse the packing quality is, the more necessary it is to optimize the blasting design; where ξ≤ξ cr , then qualified; ξ>ξ cr , then it is unqualified and the blasting design needs to be improved, where ξ cr The value is generally between 5% and 15%, which is determined according to the type of blasting project.

[0064] The following is an application example of deep hole bench blasting in an open pit mine, and combined with Figure 1 、 Figure 2 and Figure 3 , further illustrating the aforementioned rock blasting hole filling quality assessment method based on infrared thermal imaging technology.

[0065] The step height H=15m, the step surface elevation EL.1200m, the instrument used is a handheld infrared thermal imager 1, which is mounted on a tripod 2, has an infrared resolution of 384*288, a field of view of 25°*19°, a temperature measurement range of -20°C to 650°C, and a focal length of f15.

[0066] The evaluation process consists of the following steps.

[0067] The first step is to determine the photography boundary and boundary control points of the blasting area based on the blasting design parameters and surrounding environmental data, and measure the position coordinates of the control points;

[0068] The blasting design parameters and surrounding environment data of the blasting area include: the blast hole diameter in the blasting design is 152 mm, the hole spacing is 7.5 m × 4 m, the blast hole depth is 15.5 m to 16.5 m, the filling length is 4.5 m to 5.5 m, the average single hole charge is 285 kg, the explosive type is mixed emulsion explosive, the total number of blast holes is 34, and they are arranged in 3 rows, with 12 in the front row and the middle row respectively, and 10 in the back row; the environment around the blasting area is relatively simple, with the background rock being limestone. During the filming period, the atmospheric temperature was 26°C, the rock temperature was 31°C, and the relative humidity was 43%.

[0069] The theoretical volume of gas released by the detonation reaction of explosives in the blasting area is calculated based on the total mass of explosives in the blasting area determined by the blasting design. The total mass of explosives in the blasting area is 34×285Kg, and the volume of explosive products produced by the complete reaction of unit mass of explosives is 1.67m 3 Kg -1 , then V=M e v e=34×285Kg×1.67m 3 Kg -1 =14625m 3 .

[0070] The second step is instrument setup: The instrument setup point location is calculated based on the thermal imager's parameter characteristics and the coordinates of the boundary control points. The thermal imager is then set up according to the calculated results, and the shooting parameters are set. This includes calculating a safe distance of 209 meters from the center of the blast zone, the setup elevation being located at the upper step (the EL.1215m platform), the infrared thermal imager's field of view angle β, and the angle α = 4° between the corresponding edge of the field of view of the infrared thermal imager 1 and the horizontal plane of the target blast zone. The horizontal distance L of the infrared thermal imaging lens setup point from the blast zone center is calculated to be 210 meters. The tripod is set at the same height as the blast zone horizontal plane, with its center 210 meters from the blast zone center. The infrared thermal imager 1 is adjusted in height. A 50cm×50cm red flag 3 is used to mark the control point at the blast zone boundary. The instrument is set up at the designated location, and the ambient temperature parameters and the distance from the target blast zone to the shooting point are input. The acquisition time is set to 30 seconds, and the shooting frame rate is set to 300 frames / s.

[0071] The third step is thermal imaging: during the detonation process in the target blasting area, the thermal imager is used to obtain infrared thermal imaging video; including detonation after confirmation with the blaster and taking real-time thermal imaging pictures of the explosion process in the blasting area.

[0072] Step 4: Video data processing: First, use analysis software to decompose the infrared thermal image video into i The corresponding single-frame picture; then, the color cloud image that is different from the ambient temperature in the single-frame picture is mapped into a three-dimensional body to obtain the time t i Volume V of the gas product of the explosion i Finally, after screening according to the principle of non-repetitive calculation, the total volume of explosive gas escaping from the blasthole during the blasting process is obtained by summing up V 红 The process includes importing the acquired image data into professional thermal imaging analysis software, adjusting the temperature difference color band to display the corresponding relationship between the color and temperature of the current thermal image, and marking the blast holes with inconsistent temperature with the ambient temperature on the upper free surface according to the color difference range. There are 32 blast holes with a temperature of 42.6℃ on the top of the blast holes and the color blocks are uniform. The pictures are exported to digital image processing software, the contour lines are drawn and mapped into a three-dimensional body, and then the volume V is calculated. i , after screening according to the principle of non-repetitive calculation, calculate the total volume V of the explosive gas escaping the blast hole 红 =144m 3 .

[0073] Step 5: Calculation of the filling quality evaluation index: The total volume of the explosive gas escaping the blasthole V 红Substitute the volume V of the explosion zone into the formula We obtain ξ=0.98%.

[0074] Step 6: Filling quality evaluation: According to the requirements of the blasting project, take cr =6%. The principle of evaluating the quality of filling is that the smaller ξ is, the better the filling quality is; conversely, the worse the filling quality is, the more necessary it is to optimize the blasting design. cr , then qualified; ξ>ζ cr , then it is unqualified and the blasting design needs to be improved. The calculated result based on ξ=0.98% is lower than ζ cr , indicating that the filling quality is good.

[0075] At the same time, it can be seen from the uniformity of the distribution of the thermal imaging images that it is evenly distributed, proving that the filling materials and filling operation methods meet the blasting design requirements.

[0076] The above describes in detail the preferred embodiments of the present invention. It should be understood that those skilled in the art can make numerous modifications and variations based on the concepts of the present invention without inventive effort. Therefore, any technical solutions that can be derived by those skilled in the art through logical analysis, reasoning, or limited experimentation based on the concepts of the present invention and the prior art should be within the scope of protection defined by the claims.

Claims

1. A rock blast hole filling quality assessment method based on infrared thermal imaging technology, characterized in that: The following steps are involved: S1, Determination of monitoring points: Based on the blasting design parameters and surrounding environmental data of the blasting area, determine the photography boundary and boundary control points of the blasting area, and measure the position coordinates of the control points; S2, instrument installation: calculating the location of the instrument installation point based on the parameter characteristics of the thermal imager and the position coordinates of the boundary control point, installing the thermal imager according to the calculation result, and completing the shooting parameter setting; S3, thermal imaging: during the detonation process in the target blasting area, infrared thermal imaging video is obtained by using the thermal imager; S4, video data processing: first use analysis software to decompose the infrared thermal image video into i The corresponding single-frame picture; then, the color cloud image that is different from the ambient temperature in the single-frame picture is mapped into a three-dimensional body to obtain the time t i Volume V of the gas product of the explosion i Finally, after screening according to the principle of non-repetitive calculation, the total volume of explosive gas escaping from the blasthole during the blasting process is obtained by summing up V 红 ; The infrared thermal imaging video is decomposed into the time t using analysis software i The corresponding single-frame picture includes importing the infrared thermal imaging video stream data into the analysis software that matches the format, and using the analysis software to decompose the infrared thermal imaging video stream into single-frame pictures, and making each frame picture correspond to a time t i ; The color cloud image that is different from the ambient temperature in a single frame image is mapped into a three-dimensional volume to obtain the temperature at the time t i Volume V of the gas product of the explosion i , including importing the color cloud images that are different from the ambient temperature and presented in different temperature gradients in all single-frame images into digital image processing software in jpg format; then combining the position coordinates of the boundary control points measured in step S1, using the polyline function of the drawing tool to draw a closed contour line along the edge of the color block, and based on the shooting angle, the explosion area range and the time segmentation based on the delay parameter, mapping the surface formed by the closed contour line of each segment into a three-dimensional body, and obtaining the time t at that time i Volume V corresponding to different temperature gradients i ; and calculate the volume sum of various temperature gradients that are different from the ambient temperature at a single moment in the corresponding partition; After screening according to the principle of non-repetitive calculation, the total volume of explosive gas escaping from the blasthole during the blasting process in the blasting area is obtained by summing up. 红 ; This includes removing the part of the explosive gas generated by the previous blast that has not yet dissipated naturally; the maximum principle screening means only selecting the maximum value V in the sum of the explosive gas volumes in each blasting period. jmax As the basis for total volume calculation; Total volume of explosive gas V 红 Calculate as follows: In the formula, j is the specific time segment, and n is the number of segments as a natural number; S5. Determination and calculation of the filling quality evaluation index: The filling quality evaluation index is defined as ξ, which is calculated as follows: Where V 红 is the total volume of infrared thermal imaging of explosive gas products escaping from the free surface of the blasthole top; V is the theoretical volume of gas released by the detonation reaction of the total amount of explosives in the blasting area; S6, filling quality evaluation: according to the pre-determined evaluation index ξ cr The evaluation index ξ is evaluated. If it does not exceed the index value, it is qualified; if it exceeds the index value, the filling quality is improved.

2. The method according to claim 1, characterized in that In step S2, the instrument installation point is determined according to the direction of the rock mass thrown away from the blasting area and is higher than the elevation H of the top plate of the step. top The horizontal distance L between the installation point and the center of the blasting area is greater than the horizontal distance R of the earthquake source when the maximum blasting vibration velocity is 1 cm / s. v , and is less than the object distance R at the longest focal length of the target burst area d ; Elevation of erection point H 红 Determined based on the installation point location and the instrument field of view angle β.

3. The method according to claim 2, characterized in that The infrared thermal imager is installed on the top of the step, and the installation point height H 红 Calculate as follows: H 红 =H top +L*tanα, where α is the angle between the corresponding side of the field of view of the infrared thermal imager and the horizontal plane of the target explosion area, where 0≤α<90°.

4. The method according to any one of claims 1 to 3, characterized in that In step S5, V is calculated as follows: V=M e in e ; Where M e is the total mass of explosives in the blasting area, v e It is the volume of explosive products produced by complete reaction of unit mass of explosive, which is a constant based on the composition of the explosive.

5. The method according to any one of claims 1 to 3, characterized in that In step S6, the principle of evaluating the quality of packing is that the smaller ξ is, the better the packing quality is; conversely, the worse the packing quality is, the more necessary it is to optimize the blasting design; where ξ≤ξ cr , then qualified; ξ>ξ cr , then it is unqualified and the blasting design needs to be improved, ξ cr The value range is 5% to 15%.

Images