Mine air leakage rate measuring method, system, device and medium

Abstract

The application provides a mine air leakage rate measuring method and system, equipment and medium, and belongs to the technical field of mine ventilation and safety. The method comprises the following steps: determining a to-be-monitored area based on a mine grid map, dividing the to-be-monitored area into at least one target monitoring area; for each target monitoring area, performing an air leakage rate detection operation to obtain the air leakage rate corresponding to each target monitoring area; and determining the air leakage rate of the to-be-monitored area based on the air leakage rates corresponding to the target monitoring areas. The mine air leakage rate measuring method and system, equipment and medium provided by the application can accurately locate the position of an air leakage branch, accurately obtain the air leakage rate of the mine, and improve the measurement accuracy.

Description

Technical Field

[0001] This application belongs to the field of mine ventilation and safety technology, and more specifically, it relates to a method, system, equipment and medium for measuring mine air leakage rate. Background Technology

[0002] As the core site for coal resource extraction, coal mines are a key supporting unit of the energy supply system. Through a rational layout of roadways and the construction of ventilation systems, mines ensure a safe underground working environment and achieve efficient coal extraction and transportation.

[0003] Mine air leakage rate refers to the proportion of air volume that leaks out of the mine's ventilation system through gaps in roadways and air doors, bypassing designated ventilation areas (areas where airflow is not effectively utilized according to ventilation design requirements). An excessively high leakage rate can disrupt the stability of the ventilation system, leading to insufficient airflow in the ventilation areas, ineffective dilution of harmful gases, and an increased risk of poisoning for underground workers. Current methods for measuring underground air leakage rate use the traditional tracer gas method, injecting a single tracer gas into the air inlet and collecting its concentration at the return air outlet and branches. The leakage rate is then determined based on the gas concentration. While existing methods can calculate the leakage rate, they cannot accurately pinpoint the location of leaks. Since different leak locations significantly affect the accuracy of gas concentration collection, the current technology, lacking accurate location determination, also results in low accuracy in determining the leakage rate based on gas concentration. Summary of the Invention

[0004] The purpose of this application is to provide a method, system, equipment, and medium for measuring the air leakage rate in mines, so as to accurately locate the position of the air leakage branch, accurately obtain the air leakage rate of the mine, and improve the measurement accuracy.

[0005] A first aspect of this application provides a method for measuring mine air leakage rate, including:

[0006] Based on the mine grid map, the area to be monitored is determined, and the area to be monitored is divided into at least one target monitoring area;

[0007] For each target monitoring area, an air leakage rate detection operation is performed to obtain the air leakage rate corresponding to each target monitoring area. Based on the air leakage rate corresponding to each target monitoring area, the air leakage rate of the area to be monitored is determined.

[0008] The air leakage rate detection operation includes:

[0009] By uniformly injecting a composite tracer of a preset concentration into the air inlet of the target monitoring area, sensor data corresponding to each branch of the target monitoring area during the target time period are obtained. The sensor data includes: magnetic signal intensity data, fluorescence concentration data, sound wave intensity data, and airflow velocity data. The composite tracer is a combination of gaseous tracer and fluorescent nano-tracer. The sensor data of each branch is obtained based on the sensor group set at that branch.

[0010] The initial air leakage area was determined based on fluorescence concentration data and sound wave intensity data.

[0011] The target air leakage area is determined based on the initial air leakage area, magnetic signal intensity data, and airflow velocity data. The target air leakage branch is then determined based on the target air leakage area.

[0012] The leakage rate of the target air leakage branch is determined based on the concentration data of the composite tracer at the target air leakage branch and the leakage air volume data at the target air leakage branch. The leakage air volume data at the target air leakage branch is determined based on the airflow velocity data and the cross-sectional area of ​​the target air leakage branch.

[0013] Based on the leakage rate of each target's leakage branch, determine the leakage rate of the target's monitoring area.

[0014] A second aspect of this application provides a mine air leakage rate measurement system, comprising:

[0015] The target monitoring area determination module is used to determine the area to be monitored based on the mine grid map and divide the area to be monitored into at least one target monitoring area.

[0016] The air leakage rate detection module is used to perform air leakage rate detection for each target monitoring area, obtain the air leakage rate corresponding to each target monitoring area, and determine the air leakage rate of the area to be monitored based on the air leakage rate corresponding to each target monitoring area.

[0017] Specifically, the air leakage rate detection module, when performing the air leakage rate detection operation, is used for:

[0018] By uniformly injecting a composite tracer of a preset concentration into the air inlet of the target monitoring area, sensor data corresponding to each branch of the target monitoring area during the target time period are obtained. The sensor data includes: magnetic signal intensity data, fluorescence concentration data, sound wave intensity data, and airflow velocity data. The composite tracer is a combination of gaseous tracer and fluorescent nano-tracer. The sensor data of each branch is obtained based on the sensor group set at that branch.

[0019] The initial air leakage area was determined based on fluorescence concentration data and sound wave intensity data.

[0020] The target air leakage area is determined based on the initial air leakage area, magnetic signal intensity data, and airflow velocity data. The target air leakage branch is then determined based on the target air leakage area.

[0021] The leakage rate of the target air leakage branch is determined based on the concentration data of the composite tracer at the target air leakage branch and the leakage air volume data at the target air leakage branch. The leakage air volume data at the target air leakage branch is determined based on the airflow velocity data and the cross-sectional area of ​​the target air leakage branch.

[0022] Based on the leakage rate of each target's leakage branch and the concentration decay law of the composite tracer, the leakage rate of the target's monitoring area is determined.

[0023] A third aspect of this application provides an electronic device, including a memory, a processor, and a computer program stored in the memory and running on the processor, wherein the processor executes the computer program to implement the steps of the above-described method for measuring the leakage rate of a mine.

[0024] A fourth aspect of this application provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the steps of the above-described method for measuring mine air leakage rate.

[0025] The beneficial effects of the mine air leakage rate measurement method, system, equipment, and medium provided in this application embodiment are as follows: This application embodiment determines the initial air leakage area based on the fluorescence concentration data and sound wave intensity data of the target monitoring area, and then determines the target air leakage area by combining the magnetic signal intensity data and airflow velocity data, thereby determining the target air leakage branch, achieving precise positioning of the air leakage location; After the target air leakage branch (air leakage location) is accurately located, the measurement accuracy of the concentration of the composite tracer corresponding to the determined target air leakage branch is improved, thereby improving the accuracy of the air leakage rate of the target air leakage branch determined based on the concentration of the composite tracer and the cross-sectional area of ​​the branch, and thus improving the accuracy of the air leakage rate measurement of the area to be monitored. Attached Figure Description

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

[0027] Figure 1 A schematic flowchart of a method for measuring mine air leakage rate provided in an embodiment of this application;

[0028] Figure 2A schematic diagram of a rectangular area on a horizontal plane of a mine provided in an embodiment of this application;

[0029] Figure 3 This is a flowchart illustrating a method for performing air leakage rate detection according to an embodiment of this application.

[0030] Figure 4 A schematic diagram illustrating the determination of the length of the overlapping line segment between the projection of the first air leakage area and the projection of the second air leakage area, provided as an embodiment of this application;

[0031] Figure 5 A structural block diagram of a mine air leakage rate measurement system provided in an embodiment of this application;

[0032] Figure 6 This is a schematic block diagram of an electronic device provided in an embodiment of this application. Detailed Implementation

[0033] In the following description, specific details such as particular system architectures and techniques are set forth for illustrative purposes and not for limitation, in order to provide a thorough understanding of the embodiments of this application. However, those skilled in the art will understand that this application may also be implemented in other embodiments without these specific details. In other instances, detailed descriptions of well-known systems, apparatuses, circuits, and methods have been omitted so as not to obscure the description of this application with unnecessary detail.

[0034] To make the objectives, technical solutions, and advantages of this application clearer, the following description will be provided in conjunction with the accompanying drawings and specific embodiments.

[0035] Please refer to Figure 1 , Figure 1 This is a flowchart illustrating a method for measuring mine air leakage rate according to an embodiment of this application. The method may include:

[0036] S101: Determine the area to be monitored based on the mine grid map, and divide the area to be monitored into at least one target monitoring area.

[0037] In this embodiment, based on the mine grid map, the monitoring area is divided by selecting continuous grid cells. In this embodiment, workers can select all grid cells containing each target monitoring area on the grid map by drawing a box, and then combine all the selected grid cells into the monitoring area. Alternatively, grid cells can be selected based on the mine's internal roadways and their numbers, and then the selected grid cells can be combined into the monitoring area. Specifically, as... Figure 2 The diagram shown is a schematic representation of a rectangular area on a certain horizontal plane of the mine. Figure 2 The rectangular area is divided using mutually perpendicular dashed lines, forming several square grids with sides of 5 meters. The corners of each grid cell can be represented by coordinates (…). xi , y j The area is indicated by a solid line frame to clearly define its location in space. Several consecutive square grid cells are selected and enclosed with a solid line frame to form the target monitoring area, for example... Figure 2 In the diagram, A represents target monitoring area 1, and B represents target monitoring area 2.

[0038] Based on the size, shape, and monitoring range of the area to be monitored, the area to be monitored is divided into at least one target monitoring area. If the shape of the area to be monitored is regular and its range is smaller than a preset range, it can be directly used as a target monitoring area. If the shape of the area to be monitored is irregular and its range is larger than the preset range, it can be divided into multiple target monitoring areas according to the branch boundaries of the tunnel. For example, a long tunnel can be divided into target monitoring areas at preset lengths according to its distance from the mine shaft. The preset range can be adjusted according to the specific application scenario. In this embodiment, the maximum length range in three-dimensional space (length, width, height) is set to no more than 20 meters, the maximum width range to no more than 10 meters, and the maximum height range to no more than 10 meters. The area to be monitored within this three-dimensional space range can be considered a small range, and in this case, it can be used as a target monitoring area. If the range of the area to be monitored exceeds the three-dimensional space range, it can be divided into multiple target monitoring areas.

[0039] S102: For each target monitoring area, perform an air leakage rate detection operation to obtain the air leakage rate corresponding to each target monitoring area.

[0040] Before performing the air leakage rate detection operation, this embodiment first acquires parameters such as the area type, number and location of branch entrances, and mine ventilation volume of the target monitoring area. Environmental interference factors are then investigated. After confirming that there are no electromagnetic interferences, dust, or other interference factors within the target monitoring area, corresponding sensor groups are deployed according to the area type and calibrated. In this embodiment, the area type can include the mine production ventilation room and shaft area, the underground ventilation roadway area, the mining area ventilation area, and the underground sealed wall area. For example, when deploying sensor groups in the underground sealed wall area, four sensor groups can be evenly spaced within a 1-3 meter radius around the sealed wall, with two more sensor groups placed at the junction of the sealed wall and the roadway wall. Each point corresponds to one sensor group to obtain sensor data from various parts of the sealed wall. Similarly, a sensor group can be placed at each branch entrance of the underground ventilation roadway area to monitor the sensor data at each branch entrance.

[0041] like Figure 3 As shown, for each target monitoring area (e.g., for target monitoring area A), the leakage rate detection operation in S102 may include steps Sa-Sd, wherein,

[0042] Sa: A composite tracer of a preset concentration is injected at a constant speed into the air inlet of the target monitoring area to obtain sensor data corresponding to each branch of the target monitoring area during the target time period.

[0043] In this embodiment, each target monitoring area has multiple branch openings. For each target monitoring area, corresponding sensors are set according to the spatial characteristics of the target monitoring area to obtain the corresponding sensor data of each branch opening of the target monitoring area. The sensor data includes: magnetic signal intensity data, fluorescence concentration data, sound wave intensity data, and airflow velocity data. The composite tracer is a combination of gas tracer and fluorescent nano tracer. The sensor data of each branch opening is obtained based on the sensor group set at that branch opening.

[0044] Sb: Determine the initial air leakage area based on fluorescence concentration data and sound wave intensity data.

[0045] In this embodiment, for each target monitoring area, fluorescence concentration data at each branch outlet is acquired. This data is then filtered to remove outliers (e.g., abrupt changes in concentration). The trend of fluorescence concentration variation is determined based on this data. Acoustic intensity data at each branch outlet is also acquired and processed to remove noise interference. Finally, the initial air leakage area within the target monitoring area is determined based on both the fluorescence concentration data and the acoustic intensity data.

[0046] Sc: Determine the target air leakage area based on the initial air leakage area, magnetic signal intensity data, and airflow velocity data, and determine the target air leakage branch point based on the target air leakage area.

[0047] This embodiment obtains magnetic signal intensity data from each branch in the initial leakage area to determine branches with abnormal magnetic signals. A branch with abnormal magnetic signals is defined as one where the absolute difference between the magnetic signal intensity data and the mean of the standardized magnetic signal intensity values ​​is greater than a preset standard deviation. Airflow velocity data is obtained from each branch. If the absolute difference between this airflow velocity data and the mean of the standard airflow velocity is greater than a preset standard deviation, the branch is determined to have abnormal airflow velocity. Based on the branches with abnormal magnetic signals and abnormal airflow velocity, branches with multi-parameter coordinated anomalies are determined. These branches simultaneously satisfy both magnetic signal anomalies and airflow velocity anomalies. The target leakage area is determined based on the branches with multi-parameter coordinated anomalies and the initial leakage area. The branches with multi-parameter coordinated anomalies in the target leakage area are then designated as target leakage branches.

[0048] Sd: Determine the leakage rate of the target air leakage branch based on the concentration data of the composite tracer at the target air leakage branch and the leakage air volume data at the target air leakage branch. Based on the leakage rate of each target air leakage branch, determine the leakage rate of the target monitoring area.

[0049] In this embodiment, for each target monitoring area, a composite tracer of a preset concentration is injected at a uniform speed into the air inlet of the target monitoring area to obtain the concentration data of each target leakage branch. The leakage air volume data of the target leakage branch is determined based on the airflow velocity data and the cross-sectional area of ​​the target leakage branch. The cross-sectional area and airflow velocity data of each target leakage branch are obtained to determine the leakage air volume of the target leakage branch; combined with the total mine air volume, the leakage rate of the target leakage branch is determined based on the leakage air volume of the target leakage branch and the concentration of the composite tracer at the target leakage branch; then, the leakage rate of the target monitoring area is determined based on the leakage rate of each target leakage branch in the target monitoring area.

[0050] S103: Determine the air leakage rate of the area to be monitored based on the air leakage rate corresponding to each target monitoring area.

[0051] In this implementation, when the area to be monitored includes multiple target monitoring areas (two or more), the leakage rate of each target monitoring area is assigned a corresponding weight according to the ventilation volume to determine the leakage rate of the area to be monitored. When the area to be monitored includes only one target monitoring area, the leakage rate of that target monitoring area is the leakage rate of the area to be monitored.

[0052] As can be seen from the above, the embodiments of this application acquire multiple types of sensor data, such as magnetic signal intensity data, fluorescence concentration data, sound wave intensity data, and airflow velocity data, at each branch opening, and monitor the air leakage area in the mine from different dimensions. Compared with single parameter detection, it can more accurately determine the target air leakage area and the target air leakage branch opening, reducing misjudgments. This embodiment uses a composite tracer combining a gas tracer and a fluorescent nano-tracer to monitor the location of the air leakage branch opening, improving the accuracy of determining the air leakage location and achieving precise positioning of the air leakage location. After accurately locating the target air leakage branch opening (air leakage location), the measurement accuracy of the concentration of the composite tracer corresponding to the determined target air leakage branch opening is improved, thereby improving the accuracy of the air leakage rate of the target air leakage branch opening determined based on the concentration of the composite tracer and the cross-sectional area of ​​the branch opening, and further improving the accuracy of the air leakage rate measurement of the monitored area.

[0053] In one embodiment of this application, determining the initial air leakage area based on fluorescence concentration data and sound wave intensity data includes:

[0054] The sub-region corresponding to the branch that meets the first preset condition is determined as the first air leakage region. The branch that meets the first preset condition is the branch where the fluorescence concentration data shows an upward trend. One branch corresponds to one sub-region.

[0055] The sub-region corresponding to the branch that meets the second preset condition is determined as the second air leakage region. The branch that meets the second preset condition is the branch where the attenuation of the sound wave intensity data is higher than the preset attenuation level.

[0056] A spatial coordinate system is determined based on the first and second air leakage areas, and the two-dimensional and three-dimensional overlap between the first and second air leakage areas is determined based on the spatial coordinate system.

[0057] If the two-dimensional overlap is greater than or equal to the preset two-dimensional overlap threshold, and the three-dimensional overlap is greater than or equal to the preset three-dimensional overlap threshold, then the first air leakage area and the second air leakage area are determined to overlap.

[0058] The first and second air leakage areas are merged into an initial air leakage area.

[0059] In this embodiment, a first preset condition and a second preset condition are set, wherein the first preset condition is that the trend of fluorescence concentration data is an upward trend, and the second preset condition is that the attenuation of sound wave intensity data is higher than a preset attenuation level.

[0060] In this embodiment, fluorescence concentration data of each branch point within the target monitoring area at each time point is acquired, the concentration difference between adjacent fluorescence concentration data of each branch point is determined, and the branch point whose concentration difference is higher than a preset concentration threshold for a consecutive preset number of times is determined as the branch point whose fluorescence concentration data shows an upward trend. The sub-regions corresponding to all branches whose fluorescence concentration data shows an upward trend are determined as the first air leakage area.

[0061] For example, a preset acquisition time interval for fluorescence concentration data is set to 1 minute. According to this preset acquisition time interval, fluorescence concentration data for each branch in the target monitoring area is acquired at each acquisition time. For each branch, if the branch meets a first preset condition, the sub-region corresponding to that branch is determined as the first air leakage sub-region. All first air leakage sub-regions in the target monitoring area are determined as first air leakage regions. The criterion for determining whether the fluorescence data shows an upward trend is: based on the fluorescence concentration data acquired at each time, the concentration difference between adjacent times at that branch is calculated. If, within a preset number of consecutive times, this concentration difference is higher than a preset concentration threshold, then the fluorescence concentration data is determined to show an upward trend.

[0062] A sound wave transmitter is installed in the target monitoring area to emit initial sound waves. Sound wave intensity data for each branch in the target monitoring area is acquired. For each branch, if the branch meets a second preset condition, the sub-area corresponding to that branch is designated as a second air leakage sub-area. All second air leakage sub-areas in the target monitoring area are designated as second air leakage sub-areas. The criterion for determining whether a branch meets the second preset condition is: determining the sound wave attenuation degree of each branch based on the sound wave intensity data and the initial sound wave intensity data; the formula for calculating the sound wave attenuation degree is: ,in, For the first i The degree of sound wave attenuation at each branch port The initial sound wave intensity, For the first i The sound wave intensity at each branch port.

[0063] The preset attenuation level is set to 20%. For each branch, if the sound wave attenuation level of the branch is greater than the preset attenuation level, the sub-region corresponding to the branch is determined to be a sub-region of the second air leakage region. The second air leakage region is determined based on the sub-regions corresponding to all branches in the target monitoring region that meet the condition that the sound wave attenuation level is greater than the preset attenuation level.

[0064] This embodiment uses fluorescence concentration data and sound wave intensity data to detect air leakage areas, and analyzes from both chemical and physical dimensions to more comprehensively and accurately locate air leakage areas, reduce misjudgments caused by a single detection method, and thus improve the measurement efficiency of air leakage rate in the target area.

[0065] This embodiment takes the starting point of the overlapping part of the first and second air leakage areas in the mine roadway (for example, Figure 2 A spatial coordinate system is constructed with point O as the origin; the horizontal extension direction of the tunnel is set as the X-axis, representing the lateral distance in meters; the direction perpendicular to the X-axis in the horizontal direction is set as the Y-axis, representing the longitudinal distance in meters; the direction perpendicular to the horizontal plane formed by the X-axis and Y-axis is set as the Z-axis, with the positive direction of the Z-axis pointing upwards, representing the height (or depth) in meters. The Z-axis is used to characterize the vertical position of the air leakage area.

[0066] Based on the spatial coordinate system, the two-dimensional overlap and three-dimensional overlap of the first and second air leakage areas are determined respectively. If the two-dimensional overlap is greater than or equal to a preset two-dimensional overlap threshold and the three-dimensional overlap is greater than or equal to a preset three-dimensional overlap threshold, then the first and second air leakage areas are determined to overlap. The first and second air leakage areas are then merged into a whole area, and this whole area is taken as the initial air leakage area.

[0067] If the two-dimensional overlap is less than a preset two-dimensional overlap threshold or the three-dimensional overlap is less than a preset three-dimensional overlap threshold, then check whether the sensor collecting fluorescence concentration data is faulty and whether the sensor collecting sound intensity data is misaligned. If the fluorescence sensor is faulty or the sound sensor is misaligned, repair the fluorescence sensor and the sound sensor, and re-collect fluorescence concentration data and sound intensity data. Determine the first leakage area based on the fluorescence concentration data collected by the repaired fluorescence sensor, and determine the second leakage area based on the sound intensity data collected by the repaired sound sensor. Obtain the two-dimensional and three-dimensional overlap of the first and second leakage areas to determine the initial leakage area. If the fluorescence sensor is not faulty or the sound sensor is not misaligned, mark the first leakage area as a suspected fluorescence leakage area and the second leakage area as a suspected sound leakage area. Add these two suspected leakage areas to the list awaiting review for manual verification. The historical fluorescence concentration data and historical acoustic intensity data of the target monitoring area can be manually verified. If the fluorescence concentration data of the current first leakage area shows a continuous and slow upward trend compared with the historical fluorescence concentration data, then the suspected fluorescence leakage area is identified as the first leakage area. If the acoustic intensity data of the current second leakage area shows a slight increase compared with the historical acoustic data, then the suspected acoustic leakage area is identified as the second leakage area.

[0068] In one embodiment of this application, determining the two-dimensional overlap and three-dimensional overlap between the first and second air leakage areas according to a spatial coordinate system includes:

[0069] Determine the first spatial outline model based on the first air leakage area;

[0070] The second spatial outline model is determined based on the second air leakage area;

[0071] Based on the spatial coordinate system, the first spatial contour model and the second spatial contour model are projected and transformed to determine the projection of the first air leakage area and the projection of the second air leakage area.

[0072] The intersection area of ​​the first and second air leakage area projections is determined based on the first air leakage area projection and the second air leakage area projection, and the two-dimensional overlap between the first and second air leakage areas is determined based on the intersection area.

[0073] The target monitoring area is divided into a voxel grid according to a preset step size;

[0074] The voxels covered by the first air leakage area and the voxel covered by the second air leakage area are determined based on the first air leakage area, the second air leakage area, and the voxel grid.

[0075] Based on the voxels covered by the first air leakage area and the voxels covered by the second air leakage area, determine the number of overlapping voxels in the first air leakage area and the second air leakage area.

[0076] The three-dimensional overlap between the first and second air leakage regions is determined based on the number of overlapping voxels.

[0077] In this embodiment, fluorescence concentration data of different branch entrances in the mine are acquired using a fluorescence sensor. The fluorescence concentration data of each branch entrance is processed using an edge detection algorithm, and regions where the fluorescence concentration transitions from high to normal are selected. The boundary points of these transition regions are used as feature points. The three-dimensional coordinates of the feature points of the first air leakage area boundary are collected, and a first spatial contour model is formed by splicing spatial polygons. The first spatial contour model contains the spatial geometry of the feature points of all boundaries of the first air leakage area.

[0078] Acquire acoustic intensity data from different branch points using acoustic sensors, filter the data, and extract acoustic intensity within a preset frequency range (e.g., 100-1000Hz). Select the point with the largest rate of change in acoustic intensity (e.g., the point that rapidly drops from peak value to normal level) as a feature point. Acquire the 3D coordinates of the feature points at the boundary of the second air leakage area, and form a second spatial contour model by splicing spatial polygons. The second spatial contour model contains the spatial geometry of all feature points at the boundary of the second air leakage area.

[0079] Based on the spatial coordinate system, the horizontal plane (the plane formed by the X-axis and Y-axis) is selected as the projection plane. The coordinates of all boundary points in the first spatial contour model and the second spatial contour model are projected and transformed respectively. The X-axis and Y-axis coordinates are retained, and the Z-axis coordinate is 0 after projection, thus obtaining the projection of the first air leakage area and the projection of the second air leakage area.

[0080] A scan-line algorithm is used to determine the intersection region between the projections of the first and second air leakage regions. The intersection area of ​​this region is then determined; this intersection region is the two-dimensional coordinate system intersection region after projecting the first and second air leakage regions. For example... Figure 4This diagram illustrates the determination of the length of the overlapping line segment between the projections of the first and second leaky areas, where m represents the scan line. Specifically, the coordinates of all vertices in the first and second leaky area projections are extracted. These coordinates are then sorted according to the x-axis. The x-coordinates of all vertices are used as the sorting criterion, arranging them in ascending order of their x-axis values ​​to ensure that subsequent scan lines can systematically cover the x-axis range of both projections, avoiding scanning chaos. A preset scan step size of 1 is set, and the scan line is perpendicular to the x-axis. This scan line scans along the x-axis. If the scan line m passes through the projections of the two areas, there will be an overlap in the y-axis coverage of each projection along this scan line. Figure 4 In the diagram, y1 represents the y value corresponding to target monitoring area A, and y2 represents the y value corresponding to target monitoring area B. The overlapping part between y1 and y2 is taken as the overlapping line segment. The overlapping line segments are summed to determine the length of the intersection area. The length of the intersection area is then multiplied by the scanning step size to determine the intersection area.

[0081] In this embodiment, all vertex coordinates are sorted according to the x-axis in order to efficiently determine the order of intersections between the scan line (which is a straight line parallel to the y-axis) and each side of the first and second leaky area projections. This allows for accurate determination of the overlap between the first and second leaky area projections during each scan, and then calculation of the length of the overlapping line segment.

[0082] The areas of the first and second leaky regions are determined using the Gaussian area formula. The two-dimensional overlap between the first and second leaky regions is then determined based on their intersection area, the area of ​​the first leaky region's projection, and the area of ​​the second leaky region's projection. The formula for calculating the two-dimensional overlap is: ,in, The degree of overlap between the first and second air leakage areas is defined as follows: The area of ​​intersection. Let be the area projected from the first air-leaking area. This represents the projected area of ​​the second air leakage area.

[0083] This embodiment uses a spatial coordinate system as a reference, dividing the smallest cubic range of the target monitoring area, which includes the first and second air leakage areas, into several voxel grids of equal size. The side length of the voxel grid can be set according to actual accuracy requirements. For each voxel, it is determined whether its center point belongs to the first spatial contour. If the center point of a voxel belongs to the first spatial contour, it is determined that the voxel belongs to the voxels covered by the first air leakage area. All voxels in the voxel grid whose center points are inside the first spatial contour are traversed, and these voxels are determined as the set of voxels covered by the first air leakage area. For each voxel, it is determined whether its center point is inside the second spatial contour. The voxel grid is traversed, and all voxels whose center points are inside the second spatial contour are determined as the set of voxels covered by the second air leakage area.

[0084] The number of overlapping voxels in the voxel sets covered by the first and second air leakage areas is counted, along with the total number of voxel sets covered by the first and second air leakage areas. Based on the number of overlapping voxels, the total number of voxel sets covered by the first and second air leakage areas, the 3D overlap degree between the first and second air leakage areas is determined. The formula for calculating the 3D overlap degree is: ,in, The three-dimensional overlap between the first and second air leakage areas. The number of overlapping voxels. The total number of voxel sets covered by the first air leakage area. The total number of voxel sets covered by the second air leakage area.

[0085] For example, the boundary coordinates of the first spatial contour model are determined to be X∈[20,30]m, Y∈[5,8]m, and Z∈[0,3]m; the boundary coordinates of the second spatial contour model are determined to be X∈[25,35]m, Y∈[6,9]m, and Z∈[0,3]m.

[0086] Projecting the first spatial contour model and the second spatial contour model respectively, in this example, the coordinates of the first air leakage area projection are X∈[20,30]m, Y∈[5,8]m, and thus the area S1 of the first air leakage area projection is determined to be 30m². 2 The coordinates of the projection of the second air leakage area are X∈[25,35]m, Y∈[6,9]m, thus determining the area S2 of the projection of the second air leakage area to be 30m². 2 Based on the coordinates of the projections of the first and second leaky areas, the coordinates of the intersection region are determined to be X∈[25,30]m and Y∈[6,8]m, thus determining the intersection area S to be 10m². 2Based on the two-dimensional overlap calculation formula, the two-dimensional overlap between the first and second air leakage areas is determined to be 20%.

[0087] Based on the first and second spatial contour models, the minimum cube range containing the first and second air leakage areas within the target monitoring area is determined. This minimum cube range includes the first and second air leakage areas, and may also include non-air leakage areas. The range of this minimum cube is: X∈[20,35]m, Y∈[5,9]m, Z∈[0,3]m. Setting the voxel step size to 1m×1m×1m, this range is divided into 15×4×3=180 voxels. The first air leakage area covers 10×3×3=90 voxels, and the second air leakage area covers 10×3×3=90 voxels. The number of overlapping voxels between the first and second air leakage areas is 5×2×3=30 voxels. Based on the formula for calculating the three-dimensional overlap, the three-dimensional overlap between the first and second air leakage areas is determined to be 20%.

[0088] The preset two-dimensional overlap threshold is set to 50%, and the preset three-dimensional overlap threshold is also set to 50%. In this example, the two-dimensional overlap is 20%, and the three-dimensional overlap is 20%. Based on the two-dimensional and three-dimensional overlap, it can be confirmed that the two-dimensional overlap is less than the preset two-dimensional overlap threshold, and the three-dimensional overlap is less than the preset three-dimensional overlap threshold. Therefore, it is determined that the first air leakage area and the second air leakage area do not overlap. The first air leakage area is then designated as a suspected fluorescent air leakage area, and the second air leakage area is designated as a suspected acoustic air leakage area. These two suspected air leakage areas are added to the list of areas awaiting review, pending manual verification.

[0089] In one embodiment of this application, determining the target air leakage area based on the initial air leakage area, magnetic signal intensity data, and airflow velocity data includes:

[0090] The branch that meets the third preset condition is identified as the branch with magnetic signal anomaly. The branch that meets the third preset condition is the branch where the absolute difference between the magnetic signal intensity data and the mean of the standardized magnetic signal intensity value is greater than the preset intensity standard deviation.

[0091] The branch of multi-parameter coordinated anomaly is determined based on the branch of magnetic signal anomaly and the corresponding airflow velocity data. The branch of multi-parameter coordinated anomaly is the branch of magnetic signal anomaly and the branch of airflow velocity data anomaly. The branch of airflow velocity anomaly is the branch of airflow velocity data where the absolute difference between the mean of airflow velocity data and the mean of standard airflow velocity data is greater than the preset airflow standard deviation.

[0092] The spatial gradient difference of the initial air leakage region is determined based on the branch point of the multi-parameter collaborative anomaly and the initial air leakage region.

[0093] The target air leakage area is determined based on the spatial gradient difference.

[0094] In this embodiment, magnetic signal intensity data of all branch points in the target monitoring area are acquired, and the mean and standard deviation of the magnetic signal intensity data are determined based on the magnetic signal intensity data of all branch points. The magnetic signal intensity data of each branch point is standardized using the following formula: ,in, This is the normalized value of the magnetic signal intensity. Magnetic signal intensity data for each branch port, The mean of the magnetic signal intensity data. This represents the standard deviation of the magnetic signal intensity data. The mean of the standardized magnetic signal intensity values ​​is determined. If the mean of the standardized magnetic signal intensity values ​​is 0, the mean of the original magnetic signal intensity data can be used for calculation. If the mean of the standardized magnetic signal intensity values ​​is not 0, the arithmetic mean of the standardized magnetic signal intensity values ​​at all branch points can be determined, and this arithmetic mean is used as the standardized mean of the magnetic signal intensity. This standardized mean reflects the average level of the standardized magnetic signal.

[0095] Set a preset intensity standard deviation, which is used to determine whether the magnetic signal data of each branch port is abnormal.

[0096] The absolute difference between the magnetic signal intensity data of each branch and the mean of the standardized magnetic signal intensity values ​​is compared with the preset intensity standard deviation. If the absolute difference of the branch is greater than the preset intensity standard deviation, the branch is determined to be a branch with abnormal magnetic signal; if the absolute difference of the branch is less than or equal to the preset intensity standard deviation, the branch is determined to be a branch with normal magnetic signal.

[0097] Obtain airflow velocity data for all branch outlets in the target monitoring area that meet the third preset condition, determine the number of all branch outlets that meet the third preset condition, add up the airflow velocity data of all branch outlets, and then divide by the number of all branch outlets to determine the average value of the airflow velocity data.

[0098] Calculate the square of the difference between the airflow velocity data at each branch and the mean airflow velocity data. Sum these squares, divide by the number of all branches, and take the square root of the result to obtain the standard deviation. For the airflow velocity data at each branch, calculate the absolute difference between the airflow velocity data at that branch and the mean of the standard airflow velocity data. Compare this absolute difference with the preset airflow standard deviation. If the absolute difference of the branch is greater than the preset airflow standard deviation, the branch is determined to have abnormal airflow velocity; if the absolute difference of the branch is less than or equal to the preset airflow standard deviation, the branch is determined to have normal airflow velocity.

[0099] Traverse all branches with abnormal magnetic signals, and filter out branches with both abnormal magnetic signals and abnormal airflow speeds. If a branch is both a branch with abnormal magnetic signals and a branch with abnormal airflow speeds, then the branch is determined to be a branch with multiple parameters co-existing abnormalities. If a branch is a branch with abnormal magnetic signals but not a branch with abnormal airflow speeds, then the branch is determined to be a branch with a single abnormal magnetic signal.

[0100] Based on the spatial coordinate system constructed above, the spatial coordinates of the region corresponding to the multi-parameter collaborative anomaly branch are determined, as well as the set of boundary spatial coordinates corresponding to the initial air leakage region. The difference between the spatial coordinates of the multi-parameter collaborative anomaly branch and the initial air leakage region is determined based on these spatial coordinates and the spatial coordinates of the initial air leakage region. This difference represents the offset of the coordinates corresponding to the multi-parameter collaborative anomaly branch relative to the initial air leakage boundary. The spatial gradient difference is obtained by weighted summation of the differences between the spatial coordinates of all multi-parameter collaborative anomaly branches and the initial air leakage region.

[0101] For example, taking the center coordinates X0 (10,20,5) of the initial air leakage area as the spatial coordinates of the initial air leakage area, the three multi-parameter collaborative anomaly branches obtained are C, D and E, and the spatial coordinates of the three multi-parameter collaborative anomaly branches are (12,23,7), (8,18,6) and (15,25,4).

[0102] The difference between the spatial coordinates of the multi-parameter collaborative anomaly branch C and the spatial coordinates of the initial air leakage area is Δ. x C The value is 2, △ y C The value is 3, △ z C The difference between the spatial coordinates of the multi-parameter collaborative anomaly branch D and the spatial coordinates of the initial air leakage area is △. x D =-2, △ y D =-2, △ z D The difference between the spatial coordinates of the multi-parameter collaborative anomaly branch E and the spatial coordinates of the initial air leakage area is Δ. x E It is 5, △ y E It is 5, △ z E It is -1.

[0103] Weights are assigned to each coordinate direction based on the degree of influence of different coordinate directions on air leakage. For example, the weight of the X coordinate direction is 0.4, the weight of the Y coordinate direction is 0.4, and the weight of the Z coordinate direction is 0.2.

[0104] Based on the absolute values ​​of the differences between the spatial coordinates of each multi-parameter coordinated anomaly branch and the spatial coordinates of the initial air leakage area in each coordinate direction, and the corresponding weights in each coordinate direction, the weighted difference values ​​for each multi-parameter coordinated anomaly branch are determined. Based on the above data, the weighted difference values ​​for the three multi-parameter coordinated anomaly branches are determined to be 2.4, 1.8, and 4.2, respectively. Summing the weighted difference values ​​of the three multi-parameter coordinated anomaly branches, the spatial gradient difference is determined to be 8.4.

[0105] If a preset spatial gradient difference threshold is set, and the spatial gradient difference between the multi-parameter collaborative anomaly branch and the initial leakage area is less than or equal to the preset spatial gradient difference threshold, then the area where the multi-parameter collaborative anomaly branch is located is merged with the initial leakage area to obtain the target leakage area. If the spatial gradient difference between the multi-parameter collaborative anomaly branch and the initial leakage area is greater than the preset spatial gradient difference threshold, then it is further determined whether the branch of the multi-parameter collaborative anomaly belongs to a new leakage area, or the magnetic signal intensity data and airflow velocity data of each branch are reacquired, and the judgment is made again based on the reacquired data, and finally the target leakage area containing all areas that meet the conditions is determined.

[0106] In one embodiment of this application, if the target monitoring area includes multiple target air leakage branches, before determining the air leakage rate of the target air leakage branch based on the concentration data of the composite tracer at the target air leakage branch and the air leakage volume data of the target air leakage branch, the method further includes:

[0107] Acquire the airflow velocity data corresponding to each target air leakage branch opening;

[0108] The air leakage volume data of each target air leakage branch is determined based on the airflow velocity data corresponding to each target air leakage branch and the cross-sectional area of ​​the target air leakage branch.

[0109] Specifically, for each target leakage branch, the leakage rate of the target leakage branch is determined based on the concentration data of the composite tracer and the leakage air volume data of the target leakage branch, including:

[0110] Obtain the total ventilation volume of the mine;

[0111] The leakage rate of the target leakage branch is determined based on the total mine air volume, the concentration data of the composite tracer, and the leakage air volume data.

[0112] In this embodiment, airflow velocity data corresponding to each target air leakage branch is obtained using an airflow velocity measuring device (e.g., a thermal anemometer or an ultrasonic anemometer). The cross-sectional area of ​​each target air leakage branch is measured. For regularly shaped branches, the cross-sectional area can be determined using set formulas by measuring parameters such as side length and diameter; for irregularly shaped branches, the cross-sectional area can be determined using a segmentation method or a scanner.

[0113] The leakage air volume of each target leakage branch is determined based on the airflow velocity data of that branch and its corresponding cross-sectional area. The formula for calculating the leakage air volume is: ,in, For the first i The air leakage volume at each target air leakage branch outlet. For the first i Airflow velocity data at each target leaking branch opening; For the first i The cross-sectional area of ​​each target air leakage branch.

[0114] Furthermore, in this embodiment, an air volume measuring instrument is installed at the location of the mine's main intake airway or the mine's main return airway, and the total air volume of the mine is obtained based on the air volume measuring instrument.

[0115] The concentration data of the composite tracer at each target air leakage branch is obtained. In this embodiment, the concentration of the composite tracer is assumed to be uniformly distributed during the air leakage process and to follow the law of conservation of mass. The leakage rate is determined based on the total mine air volume, the concentration data of the composite tracer at each target air leakage branch, and the leakage air volume at the target air leakage branch. The formula for calculating the leakage rate is: ,in, Let be the air leakage rate of the i-th target air leakage branch. For the first i The air leakage volume at each target air leakage branch outlet. For the first i The concentration of the composite tracer at each target air leakage branch. This represents the total ventilation volume of the mine. The concentration of the composite tracer released at the source of the air leak.

[0116] This embodiment measures the airflow velocity data for each target leakage branch individually, determines the leakage volume of each target leakage branch by combining the cross-sectional area of ​​each target leakage branch, and then determines the leakage rate of each target leakage branch by combining the concentration data of the composite tracer. This helps mine workers understand the leakage distribution of the ventilation system, thereby optimizing and adjusting the ventilation system, improving its efficiency, and ensuring the ventilation safety of the mine.

[0117] For example, the measured airflow velocity at the target leakage branch 1 isv =1 = 5 m / s, the measured cross-sectional area of ​​the target air leakage branch 1 is =6m 2 The air leakage volume of the target air leakage branch outlet 1 is determined to be... =30m 3 / s; The measured airflow velocity at the target leakage branch 2 is v 2 = 4 m / s, the target air leakage branch 2 is measured to be circular with a diameter of 2 m, then the cross-sectional area of ​​the target air leakage branch 2 is... ≈3.14m 2 The air leakage volume of the target leakage branch outlet 2 is determined to be... ≈12.56m 3 / s.

[0118] The measured total ventilation volume in the mine is =500m 3 / s, assuming the concentration of the composite tracer released from the air leakage source is =10mg / m 3 The concentration of the composite tracer collected at the target air leakage branch 1 was [missing information]. =8mg / m 3 The concentration of the composite tracer collected at the target air leakage branch 2 was [value missing]. =6mg / m 3 Therefore, the leakage rate of target leakage branch 1 was determined to be 4.8%, and the leakage rate of target leakage branch 2 was approximately 1.51%.

[0119] In one embodiment of this application, the target air leakage branch includes either an outward leakage branch or an inward supply branch. Before determining the air leakage rate of the target monitoring area based on the air leakage rate of each target air leakage branch, the method further includes:

[0120] The first target air leakage branch and the second target air leakage branch are determined based on the air leakage branch of each target. The first target air leakage branch is the branch for outward leakage, and the second target air leakage branch is the branch for inward replenishment.

[0121] Based on the leakage rate of each target's leakage branch, determine the leakage rate of the target's monitoring area, including:

[0122] The total outward leakage rate within the monitoring area of ​​the target is determined based on the leakage rate of the first target's leakage branch.

[0123] The total inward supply leakage rate within the monitoring area of ​​the target is determined based on the leakage rate of the second target's leakage branch.

[0124] The total leakage rate of the target monitoring area is determined based on the total outward leakage rate and the total inward supply leakage rate.

[0125] In this embodiment, the target air leakage branch port includes either an outward leakage branch port or an inward supply branch port. In this embodiment, the first target air leakage branch port is an outward leakage branch port, and the second target air leakage branch port is an inward supply branch port. The outward leakage branch port is the target air leakage branch port from the inside of the target monitoring area to the outside of the target monitoring area, and the inward supply branch port is the target air leakage branch port from the outside of the target monitoring area to the inside of the target monitoring area.

[0126] For all target leakage branches within the target monitoring area, the first and second target leakage branches are determined based on their airflow direction. The concentration difference of the composite tracer is used to determine whether the target leakage branch is an outward leakage branch or an inward supply branch. In this embodiment, the composite tracer is located inside the target monitoring area. If the concentration of the composite tracer inside the target leakage branch is much higher than that outside, it indicates that the airflow is from the inside of the target leakage branch to the outside, and the target leakage branch is determined to be an outward leakage branch. If the concentration of the composite tracer inside the target leakage branch is much lower than that outside, it indicates that the airflow is from the outside of the target leakage branch to the inside, and the target leakage branch is determined to be an inward supply branch.

[0127] Obtain the leakage rate of each first target leakage branch, and sum the leakage rates of all first target leakage branches to obtain the total outward leakage rate of the target monitoring area. Obtain the leakage rate of each second target leakage branch, and sum the leakage rates of all second target leakage branches to obtain the total inward supply leakage rate of the target monitoring area.

[0128] The total leakage rate of the target monitoring area is determined by comparing its total outward leakage rate with its total inward supply leakage rate. The total leakage rate can be characterized by the difference between these two rates. The formula for calculating the total leakage rate of the target monitoring area is as follows: ,in, The total air leakage rate of the target monitoring area. This represents the total outward leakage rate of the target monitoring area. This represents the total inward air supply leakage rate for the target monitoring area. If... > This indicates that the entire monitoring area of ​​the target is leaking air outwards; if < This indicates that the entire monitoring area of ​​the target is supplied with external air.

[0129] This embodiment can clearly identify whether each target monitoring area is mainly leaking air outwards or mainly receiving air inwards. It can plug the branches that leak outwards and provide reasonable fresh air supply to the branches that receive air inwards, thereby optimizing the mine's ventilation system, improving ventilation efficiency, and ensuring production safety.

[0130] Corresponding to the mine air leakage rate measurement method in the above embodiment, Figure 5 This is a structural block diagram of a mine air leakage rate measurement system provided in one embodiment of this application. For ease of explanation, only the parts relevant to the embodiment of this application are shown. References Figure 5 The mine air leakage rate measurement system 50 includes: a target monitoring area determination module 51 and an air leakage rate detection module 52.

[0131] Among them, the target monitoring area determination module 51 is used to determine the area to be monitored based on the mine grid map and divide the area to be monitored into at least one target monitoring area.

[0132] The air leakage rate detection module 52 is used to perform air leakage rate detection operation for each target monitoring area, obtain the air leakage rate corresponding to each target monitoring area, and determine the air leakage rate of the area to be monitored based on the air leakage rate corresponding to each target monitoring area.

[0133] Specifically, when performing the air leakage rate detection operation, the air leakage rate detection module 52 is used for:

[0134] By uniformly injecting a composite tracer of a preset concentration into the air inlet of the target monitoring area, sensor data corresponding to each branch of the target monitoring area during the target time period are obtained. The sensor data includes: magnetic signal intensity data, fluorescence concentration data, sound wave intensity data, and airflow velocity data. The composite tracer is a combination of gaseous tracer and fluorescent nano-tracer. The sensor data of each branch is obtained based on the sensor group set at that branch.

[0135] The initial air leakage area was determined based on fluorescence concentration data and sound wave intensity data.

[0136] The target air leakage area is determined based on the initial air leakage area, magnetic signal intensity data, and airflow velocity data. The target air leakage branch is then determined based on the target air leakage area.

[0137] The leakage rate of the target air leakage branch is determined based on the concentration data of the composite tracer at the target air leakage branch and the leakage air volume data at the target air leakage branch. The leakage air volume data at the target air leakage branch is determined based on the airflow velocity data and the cross-sectional area of ​​the target air leakage branch.

[0138] Based on the leakage rate of each target's leakage branch and the concentration decay law of the composite tracer, the leakage rate of the target's monitoring area is determined.

[0139] In one embodiment of this application, when determining the initial leakage area based on fluorescence concentration data and sound wave intensity data, the leakage rate detection module 52 is specifically used for:

[0140] The sub-region corresponding to the branch that meets the first preset condition is determined as the first air leakage region. The branch that meets the first preset condition is the branch where the fluorescence concentration data shows an upward trend. One branch corresponds to one sub-region.

[0141] The sub-region corresponding to the branch that meets the second preset condition is determined as the second air leakage region. The branch that meets the second preset condition is the branch where the attenuation of the sound wave intensity data is higher than the preset attenuation level.

[0142] A spatial coordinate system is determined based on the first and second air leakage areas, and the two-dimensional and three-dimensional overlap between the first and second air leakage areas is determined based on the spatial coordinate system.

[0143] If the two-dimensional overlap is greater than or equal to the preset two-dimensional overlap threshold, and the three-dimensional overlap is greater than or equal to the preset three-dimensional overlap threshold, then the first air leakage area and the second air leakage area are determined to overlap.

[0144] The first and second air leakage areas are merged into an initial air leakage area.

[0145] In one embodiment of this application, when the air leakage rate detection module 52 determines the two-dimensional overlap and three-dimensional overlap between the first air leakage area and the second air leakage area according to the spatial coordinate system, it is specifically used for:

[0146] Determine the first spatial outline model based on the first air leakage area;

[0147] The second spatial outline model is determined based on the second air leakage area;

[0148] Based on the spatial coordinate system, the first spatial contour model and the second spatial contour model are projected and transformed to determine the projection of the first air leakage area and the projection of the second air leakage area.

[0149] The intersection area of ​​the first and second air leakage area projections is determined based on the first air leakage area projection and the second air leakage area projection, and the two-dimensional overlap between the first and second air leakage areas is determined based on the intersection area.

[0150] The target monitoring area is divided into a voxel grid according to a preset step size;

[0151] The voxels covered by the first air leakage area and the voxel covered by the second air leakage area are determined based on the first air leakage area, the second air leakage area, and the voxel grid.

[0152] Based on the voxels covered by the first air leakage area and the voxels covered by the second air leakage area, determine the number of overlapping voxels in the first air leakage area and the second air leakage area.

[0153] The three-dimensional overlap between the first and second air leakage regions is determined based on the number of overlapping voxels.

[0154] In one embodiment of this application, when the air leakage rate detection module 52 determines the target air leakage area based on the initial air leakage area, magnetic signal intensity data, and airflow velocity data, it is specifically used for:

[0155] The branch that meets the third preset condition is identified as the branch with magnetic signal anomaly. The branch that meets the third preset condition is the branch where the absolute difference between the magnetic signal intensity data and the mean of the standardized magnetic signal intensity value is greater than the preset intensity standard deviation.

[0156] The branch of multi-parameter coordinated anomaly is determined based on the branch of magnetic signal anomaly and the corresponding airflow velocity data. The branch of multi-parameter coordinated anomaly is the branch of magnetic signal anomaly and the branch of airflow velocity data anomaly. The branch of airflow velocity anomaly is the branch of airflow velocity data where the absolute difference between the mean of airflow velocity data and the mean of standard airflow velocity data is greater than the preset airflow standard deviation.

[0157] The spatial gradient difference of the initial air leakage region is determined based on the branch point of the multi-parameter collaborative anomaly and the initial air leakage region.

[0158] The target air leakage area is determined based on the spatial gradient difference.

[0159] In one embodiment of this application, if the target monitoring area includes multiple target air leakage branch outlets, the air leakage rate detection module 52, before determining the air leakage rate of the target air leakage branch outlet based on the concentration data of the composite tracer at the target air leakage branch outlet and the air leakage volume data of the target air leakage branch outlet, is further configured to:

[0160] Acquire the airflow velocity data corresponding to each target air leakage branch opening;

[0161] The air leakage volume data of each target air leakage branch is determined based on the airflow velocity data corresponding to each target air leakage branch and the cross-sectional area of ​​the target air leakage branch.

[0162] Specifically, for each target leakage branch, the leakage rate detection module 52, when determining the leakage rate of the target leakage branch based on the concentration data of the composite tracer and the leakage air volume data of the target leakage branch, is used for:

[0163] Obtain the total ventilation volume of the mine;

[0164] The leakage rate of the target leakage branch is determined based on the total mine air volume, the concentration data of the composite tracer, and the leakage air volume data.

[0165] In one embodiment of this application, the target air leakage branch port includes an outward leakage branch port or an inward supply branch port. Before determining the air leakage rate of the target monitoring area based on the air leakage rate of each target air leakage branch port, the air leakage rate detection module 52 is further configured to:

[0166] The first target air leakage branch and the second target air leakage branch are determined based on the air leakage branch of each target. The first target air leakage branch is the branch for outward leakage, and the second target air leakage branch is the branch for inward replenishment.

[0167] When determining the air leakage rate of the target monitoring area based on the air leakage rate of each target air leakage branch, the air leakage rate detection module 52 is specifically used for:

[0168] The total outward leakage rate within the monitoring area of ​​the target is determined based on the leakage rate of the first target's leakage branch.

[0169] The total inward supply leakage rate within the monitoring area of ​​the target is determined based on the leakage rate of the second target's leakage branch.

[0170] The total leakage rate of the target monitoring area is determined based on the total outward leakage rate and the total inward supply leakage rate.

[0171] In one embodiment of this application, when the air leakage rate detection module 52 determines the sub-region corresponding to the branch opening that meets the first preset condition as the first air leakage region, it is specifically used for:

[0172] Obtain the fluorescence concentration data of each branch in the target monitoring area at each time, determine the concentration difference between the fluorescence concentration data of each branch at adjacent times, and determine the branch whose concentration difference is higher than the preset concentration threshold for a consecutive preset number of times as the branch whose fluorescence concentration data shows an upward trend. Determine the sub-area corresponding to all branches whose fluorescence concentration data shows an upward trend as the first air leakage area.

[0173] When the air leakage rate detection module 52 determines the sub-region corresponding to the branch port that meets the second preset condition as the second air leakage region, it is specifically used for:

[0174] Acquire the acoustic intensity data of each branch in the target monitoring area, determine the acoustic attenuation degree of each branch based on the acoustic intensity data of each branch and the initial acoustic intensity data, and determine the sub-area corresponding to all branches with attenuation degrees higher than the preset attenuation degree as the second air leakage area.

[0175] See Figure 6 , Figure 6 This is a schematic block diagram of an electronic device provided according to an embodiment of this application. Figure 6 The electronic device 600 in this embodiment may include one or more processors 601, one or more input devices 602, one or more output devices 603, and one or more memories 604. The processors 601, input devices 602, output devices 603, and memories 604 communicate with each other via a communication bus 605. The memory 604 stores computer programs, including program instructions. The processor 601 executes the program instructions stored in the memory 604. The processor 601 is configured to invoke the program instructions to perform the functions of each module / unit in the above system embodiments, for example... Figure 5 The functions of the target monitoring area determination module 51 and the air leakage rate detection module 52 are shown.

[0176] It should be understood that, in the embodiments of this application, the processor 601 may be a central processing unit (CPU), but it may also be other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The general-purpose processor may be a microprocessor or any conventional processor.

[0177] Input device 602 may include a touchpad, a fingerprint sensor (for collecting the user's fingerprint information and fingerprint orientation information), a microphone, etc., and output device 603 may include a display (LCD, etc.), a speaker, etc.

[0178] The memory 604 may include read-only memory and random access memory, and provides instructions and data to the processor 601. A portion of the memory 604 may also include non-volatile random access memory. For example, the memory 604 may also store sensor data, composite tracer concentration data, target leakage areas, and leakage rate data of target leakage branches in each target monitoring area during the target time period.

[0179] In specific implementations, the processor 601, input device 602, and output device 603 described in the embodiments of this application can execute the implementation method described in the mine air leakage rate measurement method provided in the embodiments of this application, or they can execute the implementation method of the electronic device described in the embodiments of this application, which will not be repeated here.

[0180] In another embodiment of this application, a computer-readable storage medium is provided. This computer-readable storage medium stores a computer program, which includes program instructions. When executed by a processor, the program instructions implement all or part of the processes in the methods described above. Alternatively, the computer program can instruct related hardware to implement these processes. The computer program can be stored in a computer-readable storage medium, and when executed by a processor, it can implement the steps of the various method embodiments described above. The computer program includes computer program code, which can be in the form of source code, object code, executable files, or certain intermediate forms. The computer-readable medium can include any entity or system capable of carrying computer program code, a recording medium, a USB flash drive, a portable hard drive, a magnetic disk, an optical disk, a computer memory, a read-only memory (ROM), a random access memory (RAM), an electrical carrier signal, a telecommunication signal, and a software distribution medium, etc.

[0181] The computer-readable storage medium can be an internal storage unit of the electronic device in any of the foregoing embodiments, such as a hard disk or memory of the electronic device. The computer-readable storage medium can also be an external storage device of the electronic device, such as a plug-in hard disk, smart media card (SMC), secure digital card (SD), flash card, etc., provided on the electronic device. Furthermore, the computer-readable storage medium can include both internal and external storage units of the electronic device. The computer-readable storage medium is used to store computer programs and other programs and data required by the electronic device. The computer-readable storage medium can also be used to temporarily store data that has been output or will be output.

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

[0183] Those skilled in the art will clearly understand that, for the sake of convenience and brevity, the specific working process of the electronic devices and units described above can be referred to the corresponding process in the foregoing method embodiments, and will not be repeated here.

[0184] In the embodiments provided in this application, it should be understood that the disclosed electronic devices and methods can be implemented in other ways. For example, the system embodiments described above are merely illustrative. For instance, the division of modules / units is only a logical functional division, and in actual implementation, there may be other division methods. For example, multiple modules, units, or components may be combined or integrated into another system, or some features may be ignored or not executed. In addition, the mutual coupling or direct coupling or communication connection shown or discussed may be indirect coupling or communication connection through some interfaces or modules / units, or it may be an electrical, mechanical, or other form of connection.

[0185] The modules / units described as separate components may or may not be physically separate. Similarly, the components shown as modules / units may or may not be physical modules / units; they may be located in one place or distributed across multiple network modules / units. Some or all of the modules / units can be selected to achieve the purpose of the embodiments of this application, depending on actual needs.

[0186] Furthermore, the functional modules in the various embodiments of this application can be integrated into one processing unit, or each module can exist physically separately, or two or more modules can be integrated into one unit. The integrated modules / units described above can be implemented in hardware or in the form of software functional modules / units.

[0187] The above are merely specific embodiments of this application, but the scope of protection of this application is not limited thereto. Any person skilled in the art can easily conceive of various equivalent modifications or substitutions within the technical scope disclosed in this application, and these modifications or substitutions should all be covered within the scope of protection of this application. Therefore, the scope of protection of this application should be determined by the scope of the claims.

Claims

1. A method for measuring air leakage rate in mines, characterized in that, include: Based on the mine grid map, the area to be monitored is determined, and the area to be monitored is divided into at least one target monitoring area; For each target monitoring area, an air leakage rate detection operation is performed to obtain the air leakage rate corresponding to each target monitoring area. Based on the air leakage rate corresponding to each target monitoring area, the air leakage rate of the area to be monitored is determined. The air leakage rate detection operation includes: By uniformly injecting a composite tracer of a preset concentration into the air inlet of the target monitoring area, sensor data corresponding to each branch of the target monitoring area during the target time period are obtained. The sensor data includes: magnetic signal intensity data, fluorescence concentration data, sound wave intensity data, and airflow velocity data. The composite tracer is a combination of gaseous tracer and fluorescent nano-tracer. The sensor data of each branch is obtained based on the sensor group set at that branch. The sub-region corresponding to the branch that meets the first preset condition is determined as the first air leakage region. The branch that meets the first preset condition is the branch whose fluorescence concentration data shows an upward trend. One branch corresponds to one sub-region. The sub-region corresponding to the branch that meets the second preset condition is determined as the second air leakage region. The branch that meets the second preset condition is the branch where the attenuation of the sound wave intensity data is higher than the preset attenuation level. A spatial coordinate system is determined based on the first air leakage area and the second air leakage area, and the two-dimensional overlap and three-dimensional overlap between the first air leakage area and the second air leakage area are determined based on the spatial coordinate system. If the two-dimensional overlap is greater than or equal to a preset two-dimensional overlap threshold, and the three-dimensional overlap is greater than or equal to a preset three-dimensional overlap threshold, then the first air leakage area and the second air leakage area are determined to overlap. The first and second leak areas are merged into an initial leak area; The branch that meets the third preset condition is identified as the branch with abnormal magnetic signal. The branch that meets the third preset condition is the branch where the absolute difference between the magnetic signal intensity data and the mean of the standardized magnetic signal intensity value is greater than the preset intensity standard deviation. The branch of multi-parameter coordinated anomaly is determined based on the branch of magnetic signal anomaly and the airflow velocity data corresponding to the branch of magnetic signal anomaly. The branch of multi-parameter coordinated anomaly is the branch of magnetic signal anomaly and the branch of airflow velocity data anomaly. The branch of airflow velocity anomaly is the branch where the absolute difference between the mean of the airflow velocity data and the mean of the standard airflow velocity data is greater than the preset airflow standard deviation. The spatial gradient difference of the initial air leakage region is determined based on the branch point of the multi-parameter collaborative anomaly and the initial air leakage region. The target air leakage area is determined based on the spatial gradient difference; Determine the target air leakage branch point based on the target air leakage area; The leakage rate of the target air leakage branch is determined based on the concentration data of the composite tracer at the target air leakage branch and the leakage air volume data at the target air leakage branch. The leakage air volume data at the target air leakage branch is determined based on the airflow velocity data and the cross-sectional area of ​​the target air leakage branch. Based on the leakage rate of each target's leakage branch, determine the leakage rate of the target's monitoring area.

2. The method for measuring mine air leakage rate as described in claim 1, characterized in that, Determining the two-dimensional and three-dimensional overlap between the first and second air leakage areas according to the spatial coordinate system includes: Determine the first spatial outline model based on the first air leakage area; The second spatial outline model is determined based on the second air leakage area; Based on the spatial coordinate system, the first spatial contour model and the second spatial contour model are projected to determine the projection of the first air leakage area and the projection of the second air leakage area. The intersection area of ​​the first air leakage area projection and the second air leakage area projection is determined based on the first air leakage area projection and the second air leakage area projection, and the two-dimensional overlap between the first air leakage area and the second air leakage area is determined based on the intersection area. The target monitoring area is divided into a voxel grid according to a preset step size; The voxels covered by the first air leakage area and the voxels covered by the second air leakage area are determined based on the first air leakage area, the second air leakage area, and the voxel grid. Based on the voxels covered by the first air leakage area and the voxels covered by the second air leakage area, determine the number of overlapping voxels among the voxels covered by the first air leakage area and the voxels covered by the second air leakage area. The three-dimensional overlap between the first and second air leakage regions is determined based on the number of overlapping voxels.

3. The method for measuring mine air leakage rate as described in claim 1, characterized in that, If the target monitoring area contains multiple target air leakage outlets, before determining the air leakage rate of the target air leakage outlet based on the concentration data of the composite tracer at the target air leakage outlet and the air leakage volume data of the target air leakage outlet, the following steps are also included: Acquire the airflow velocity data corresponding to each target air leakage branch opening; The air leakage volume data of each target air leakage branch is determined based on the airflow velocity data corresponding to each target air leakage branch and the cross-sectional area of ​​the target air leakage branch. Specifically, for each target leakage branch, the leakage rate of the target leakage branch is determined based on the concentration data of the composite tracer and the leakage air volume data of the target leakage branch, including: Obtain the total ventilation volume of the mine; The leakage rate of the target leakage branch is determined based on the total mine air volume, the concentration data of the composite tracer, and the leakage air volume data.

4. The method for measuring mine air leakage rate as described in claim 3, characterized in that, The target air leakage branch includes an outward leakage branch or an inward supply branch. Before determining the air leakage rate of the target monitoring area based on the air leakage rate of each target air leakage branch, the method further includes: The first target air leakage branch and the second target air leakage branch are determined based on each target air leakage branch. The first target air leakage branch is the branch for outward leakage, and the second target air leakage branch is the branch for inward replenishment. The step of determining the air leakage rate of the target monitoring area based on the air leakage rate of each target air leakage branch includes: The total outward leakage rate within the monitoring area of ​​the target is determined based on the leakage rate of the first target leakage branch. The total inward supply leakage rate within the monitoring area of ​​the target is determined based on the leakage rate of the second target leakage branch. The total leakage rate of the target monitoring area is determined based on the total outward leakage rate and the total inward supply leakage rate.

5. The method for measuring mine air leakage rate as described in claim 1, characterized in that, The step of defining the sub-region corresponding to the branch that meets the first preset condition as the first air leakage region includes: Obtain the fluorescence concentration data of each branch in the target monitoring area at each time, determine the concentration difference between the fluorescence concentration data of each branch at adjacent times, and determine the branch whose concentration difference is higher than the preset concentration threshold for a consecutive preset number of times as the branch whose fluorescence concentration data shows an upward trend. Determine the sub-area corresponding to all branches whose fluorescence concentration data shows an upward trend as the first air leakage area. The step of defining the sub-region corresponding to the branch that meets the second preset condition as the second air leakage region includes: Acquire the acoustic intensity data of each branch in the target monitoring area, determine the acoustic attenuation degree of each branch based on the acoustic intensity data of each branch and the initial acoustic intensity data, and determine the sub-area corresponding to all branches with attenuation degrees higher than the preset attenuation degree as the second air leakage area.

6. A mine air leakage rate measurement system, characterized in that, include: The target monitoring area determination module is used to determine the area to be monitored based on the mine grid map and divide the area to be monitored into at least one target monitoring area. The air leakage rate detection module is used to perform air leakage rate detection operation for each target monitoring area, obtain the air leakage rate corresponding to each target monitoring area, and determine the air leakage rate of the area to be monitored based on the air leakage rate corresponding to each target monitoring area. Specifically, when performing the air leakage rate detection operation, the air leakage rate detection module is used for: By uniformly injecting a composite tracer of a preset concentration into the air inlet of the target monitoring area, sensor data corresponding to each branch of the target monitoring area during the target time period are obtained. The sensor data includes: magnetic signal intensity data, fluorescence concentration data, sound wave intensity data, and airflow velocity data. The composite tracer is a combination of gaseous tracer and fluorescent nano-tracer. The sensor data of each branch is obtained based on the sensor group set at that branch. The sub-region corresponding to the branch that meets the first preset condition is determined as the first air leakage region. The branch that meets the first preset condition is the branch whose fluorescence concentration data shows an upward trend. One branch corresponds to one sub-region. The sub-region corresponding to the branch that meets the second preset condition is determined as the second air leakage region. The branch that meets the second preset condition is the branch where the attenuation of the sound wave intensity data is higher than the preset attenuation level. A spatial coordinate system is determined based on the first air leakage area and the second air leakage area, and the two-dimensional overlap and three-dimensional overlap between the first air leakage area and the second air leakage area are determined based on the spatial coordinate system. If the two-dimensional overlap is greater than or equal to a preset two-dimensional overlap threshold, and the three-dimensional overlap is greater than or equal to a preset three-dimensional overlap threshold, then the first air leakage area and the second air leakage area are determined to overlap. The first and second leak areas are merged into an initial leak area; The branch that meets the third preset condition is identified as the branch with abnormal magnetic signal. The branch that meets the third preset condition is the branch where the absolute difference between the magnetic signal intensity data and the mean of the standardized magnetic signal intensity value is greater than the preset intensity standard deviation. The branch of multi-parameter coordinated anomaly is determined based on the branch of magnetic signal anomaly and the airflow velocity data corresponding to the branch of magnetic signal anomaly. The branch of multi-parameter coordinated anomaly is the branch of magnetic signal anomaly and the branch of airflow velocity data anomaly. The branch of airflow velocity anomaly is the branch where the absolute difference between the mean of the airflow velocity data and the mean of the standard airflow velocity data is greater than the preset airflow standard deviation. The spatial gradient difference of the initial air leakage region is determined based on the branch point of the multi-parameter collaborative anomaly and the initial air leakage region. The target air leakage area is determined based on the spatial gradient difference; Determine the target air leakage branch point based on the target air leakage area; The leakage rate of the target air leakage branch is determined based on the concentration data of the composite tracer at the target air leakage branch and the leakage air volume data at the target air leakage branch. The leakage air volume data at the target air leakage branch is determined based on the airflow velocity data and the cross-sectional area of ​​the target air leakage branch. Based on the leakage rate of each target's leakage branch and the concentration decay law of the composite tracer, the leakage rate of the target's monitoring area is determined.

7. An electronic device comprising a memory, a processor, and a computer program stored in the memory and running on the processor, characterized in that, When the processor executes the computer program, it implements the steps of the method as described in any one of claims 1 to 5.

8. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed by a processor, it implements the steps of the method as described in any one of claims 1 to 5.

Images