Ultrasonic precision wind measurement and resistance online monitoring method and related equipment

Abstract

The present application belongs to the technical field of mine ventilation and resistance monitoring, and specifically provides an ultrasonic precision wind measurement and resistance online monitoring method and related equipment. The method mainly includes: obtaining the original sound wave characteristics collected by the ultrasonic transducer array, integrating and generating an initial feature data set; based on the initial feature data set and the underground propagation calibration value, generating an anti-dust propagation time difference set by fusion; parsing the tunnel space grid according to the anti-dust propagation time difference set, generating a resistance hot zone distribution map; solving the graded resistance state index by comparing the dynamic flow field of the resistance hot zone distribution map; generating a spatial coordinate safety instruction based on the positioning result of the graded resistance state index. The present application can accurately correct the interference of multiple factors such as dust on sound wave propagation in complex environments, significantly improve the accuracy of propagation time difference calculation, optimize the tunnel resistance distribution analysis and the spatial coordinate generation of safety instructions, and effectively ensure the accuracy and reliability of mine ventilation monitoring.

Description

Technical Field

[0001] The present application belongs to the technical field of mine ventilation and resistance monitoring, and particularly relates to an ultrasonic precise wind measurement and resistance online monitoring method and related equipment. Background Art

[0002] In the field of mine ventilation and resistance monitoring, various sensing methods are often used to collect data, and wind measurement and resistance monitoring are achieved by analyzing physical quantities such as sound waves and airflow. Generally, the propagation characteristics of sound waves are used to collect original sound wave information and perform preliminary processing. The propagation time difference is calculated in combination with simple calibration data to analyze the relevant parameters of the tunnel space, and then a resistance distribution map is drawn. The flow field is then compared based on the map, the resistance status is solved, and finally a safety instruction is generated.

[0003] However, the mine environment is complex, with many interference factors such as dust. When sound wave propagation is affected by dust, the lack of precise correction of propagation characteristics and anti-interference fusion processing can easily lead to deviations in the calculation of propagation time differences, resulting in inaccurate tunnel spatial analysis and resistance map drawing, which in turn affects the accuracy of resistance state solution, resulting in limited accuracy in the spatial coordinate generation of safety instructions, making it difficult to reliably meet the needs of safe and efficient mine ventilation monitoring. Summary of the Invention

[0004] This application effectively solves the problem of limited spatial coordinate accuracy of safety instruction generation in the existing technology during mine ventilation and resistance monitoring, which is caused by the lack of deep adaptive processing of multi-factor interference in complex environments in the existing technology, by providing an ultrasonic precision wind measurement and resistance online monitoring method and related equipment. It can accurately correct complex environmental interference, improve the accuracy of propagation time difference calculation, optimize resistance distribution analysis and instruction generation, and ensure the accuracy and reliability of mine ventilation monitoring.

[0005] In order to achieve the above objectives, this application adopts the following technical solutions:

[0006] In a first aspect, the present application provides a method for accurate ultrasonic wind measurement and online monitoring of resistance, comprising:

[0007] The original sound wave features collected by the ultrasonic transducer array are obtained and integrated to generate an initial feature data set.

[0008] Based on the initial feature data set and the underground propagation calibration value, the anti-dust propagation time difference set is generated by fusion.

[0009] The spatial grid of the roadway is analyzed based on the set of resistance dust propagation time differences to generate a distribution map of resistance hotspots.

[0010] The graded resistance state index is calculated by comparing the dynamic flow field of the resistance hot zone distribution map.

[0011] Based on the positioning results of the graded resistance status indicators, spatial coordinate safety instructions are generated.

[0012] Furthermore, the original sound wave features collected by the ultrasonic transducer array are obtained and integrated to generate an initial feature data set, including:

[0013] Obtain the absolute completion time of the measured effective echo signal in the opposite propagation direction; calculate the absolute value of the frequency difference between the carrier frequency of the effective echo signal and the carrier frequency of the transmitted signal, and record it as the frequency offset characteristic value.

[0014] The absolute completion time and frequency offset feature value of the same path are bound and integrated into the initial feature data set.

[0015] The time delay of the effective echo signal is within the range of the preset maximum propagation time of the lane.

[0016] Furthermore, based on the initial feature data set and the downhole propagation calibration value, a set of anti-dust propagation time difference values ​​is generated by fusion, including:

[0017] The frequency offset equivalent time correction amount is calculated based on the frequency offset characteristic value of each path in the initial characteristic data set and the calibration value of the downhole ultrasonic propagation velocity; the basic propagation time difference is calculated based on the absolute completion time in the opposite direction of each path in the initial characteristic data set.

[0018] The basic propagation time difference is superimposed with the frequency offset equivalent time correction to generate a fused anti-interference propagation time difference.

[0019] Generate a fused propagation time difference value set including the fused interference-rejected propagation time difference values ​​of all paths.

[0020] Furthermore, the roadway spatial grid is analyzed based on the set of anti-dust propagation time differences to generate a resistance hotspot distribution map, including:

[0021] Based on the pre-built three-dimensional digital model of the tunnel, a tunnel space grid model is constructed.

[0022] The wind speed vector value at the grid node position is calculated by using the fused propagation time difference set in combination with the position of the ultrasonic path in the tunnel space grid model.

[0023] Calculate the spatial variation gradient of the wind speed vector values ​​of adjacent grid nodes.

[0024] The resistance impact hotspot area and the benchmark flow field reference area are marked according to the spatial variation gradient value.

[0025] Furthermore, the drag impact hotspot area and the benchmark flow field reference area are marked according to the spatial variation gradient value, including:

[0026] A first threshold and a second threshold are set.

[0027] Grid nodes whose spatial variation gradient values ​​exceed the first threshold are identified and marked as resistance impact hotspots; grid nodes whose spatial variation gradient values ​​are lower than the second threshold are identified and marked as baseline flow field reference areas.

[0028] Furthermore, by comparing the dynamic flow field of the resistance hot zone distribution map, the graded resistance state indicators are calculated, including:

[0029] Extract the wind speed vector values ​​of the grid nodes in the reference area of ​​the benchmark flow field and calculate the average benchmark wind speed vector value.

[0030] For each grid node in the resistance-affecting hotspot area: calculate the wind speed vector difference between the wind speed vector value and the average reference wind speed vector value; and solve the dynamic resistance factor based on the wind speed vector difference and the node ventilation cross-sectional area.

[0031] Aggregate dynamic resistance factors to generate a comprehensive resistance state assessment value.

[0032] When the comprehensive resistance status evaluation value exceeds the preset warning threshold and the resistance impact hotspot coincides with the coordinates of the gas sensor, the resistance-gas coupling warning mark is triggered.

[0033] Furthermore, based on the positioning results of the graded resistance state indicators, spatial coordinate safety instructions are generated, including:

[0034] The physical blockage point is located according to the spatial distribution coordinates of the resistance hotspot area and associated with the roadway equipment database.

[0035] Set the third threshold. If the rate of change of the dynamic resistance factor reaches the third threshold and the damper coordinates are located, a damper opening adjustment instruction is generated. If the resistance-gas coupling warning mark is activated, a strong ventilation instruction + personnel evacuation instruction is generated. If the resistance-affected hotspot area persists for a specified time, a patrol robot cleaning instruction is generated.

[0036] The dynamic resistance factor is mapped to the three-dimensional digital model of the tunnel to generate a real-time resistance heat map.

[0037] In a second aspect, the present application provides an ultrasonic precision wind measurement and drag online monitoring system, which includes:

[0038] Original sound wave feature acquisition module: obtains the original sound wave features collected by the ultrasonic transducer array and integrates them to generate an initial feature data set.

[0039] Anti-dust propagation feature fusion module: Based on the initial feature data set and the downhole propagation calibration value, it fuses and generates an anti-dust propagation time difference set.

[0040] Tunnel space grid analysis module: Analyze the tunnel space grid based on the anti-dust propagation time difference set to generate the resistance hot zone distribution map.

[0041] Resistance state dynamic solution module: solves the graded resistance state index through dynamic flow field comparison of resistance hot zone distribution map.

[0042] Safety instruction coordinate generation module: Generates spatial coordinate safety instructions based on the positioning results of the graded resistance state indicators.

[0043] In a third aspect, the present application provides an ultrasonic precision wind measurement and online resistance monitoring device, which includes a memory and a processor; the memory is used to store computer programs; the processor is used to implement the steps of the ultrasonic precision wind measurement and online resistance monitoring method described in the first aspect when executing the computer program.

[0044] In a fourth aspect, a readable storage medium stores computer program instructions. When the computer program instructions are read and executed by a processor, the steps of the ultrasonic precision wind measurement and resistance online monitoring method described in the first aspect are executed.

[0045] Beneficial effects of this application:

[0046] This application forms an initial data set based on the acquisition of acoustic wave characteristics by an ultrasonic array, fuses anti-dust time difference data with underground calibration values, analyzes and generates a roadway resistance hot zone map, and solves the graded resistance index through flow field comparison. It outputs spatial coordinate safety instructions based on the index positioning results, effectively solving the problem of the lack of deep adaptive processing of multi-factor interference in complex environments in mine ventilation and resistance monitoring in the existing technology, resulting in limited spatial coordinate accuracy of safety instruction generation. It can accurately correct the interference of multiple factors such as dust in complex environments on sound wave propagation, significantly improve the calculation accuracy of propagation time difference, optimize the roadway resistance distribution analysis and the spatial coordinate generation of safety instructions, and effectively ensure the accuracy and reliability of mine ventilation monitoring.

[0047] Other features and advantages of the present application will be described in the following description, and in part will become apparent from the description, or will be understood by practicing the present application. The purpose and other advantages of the present application can be realized and obtained by the structures indicated in the description and the drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0048] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, a brief introduction will be given below to the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without any creative work.

[0049] Figure 1 A schematic diagram of the process of an ultrasonic precise wind measurement and resistance online monitoring method of the present application is shown;

[0050] Figure 2 A module schematic diagram of an ultrasonic precision wind measurement and resistance online monitoring system of the present application is shown. DETAILED DESCRIPTION

[0051] In order to solve the problems raised by the background technology, this application forms an initial data set based on the acquisition of acoustic wave characteristics by an ultrasonic array, combines the anti-dust time difference data with the underground calibration value, analyzes and generates the tunnel resistance hot zone map and solves the graded resistance index through flow field comparison, and outputs spatial coordinate safety instructions based on the indicator positioning results. It can accurately correct complex environmental interference, improve the calculation accuracy of propagation time difference, optimize resistance distribution analysis and instruction generation, and ensure the accuracy and reliability of mine ventilation monitoring.

[0052] To make the purpose, technical solutions, and advantages of the embodiments of this application more clear, the technical solutions in the embodiments of this application will be clearly and completely described below in conjunction with the drawings in the embodiments of this application. Obviously, the described embodiments are part of the embodiments of this application, not all of the embodiments. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of this application.

[0053] In some embodiments, as Figure 1 As shown, the present application provides an ultrasonic precise wind measurement and drag online monitoring method, comprising:

[0054] S1. Obtain the original sound wave features collected by the ultrasonic transducer array and integrate them to generate an initial feature data set.

[0055] S2. Based on the initial feature data set and the downhole propagation calibration value, generate a set of anti-dust propagation time difference values ​​by fusion.

[0056] S3. Analyze the tunnel space grid based on the anti-dust propagation time difference set to generate a resistance hot zone distribution map.

[0057] S4. Calculate the graded resistance state index by comparing the dynamic flow field of the resistance hot zone distribution map.

[0058] S5. Generate spatial coordinate safety instructions based on the positioning results of the graded resistance state indicators.

[0059] In some embodiments, in S1, original acoustic wave features collected by the ultrasonic transducer array are obtained and integrated to generate an initial feature data set, including:

[0060] S11. Obtain the absolute completion time of the measured effective echo signal in the opposite propagation direction; calculate the absolute value of the frequency difference between the carrier frequency of the effective echo signal and the carrier frequency of the transmitted signal, and record it as the frequency offset characteristic value; wherein, the time delay of the effective echo signal is within the preset maximum propagation time range of the channel.

[0061] Ultrasonic transducer arrays are deployed at key nodes in mine tunnels, covering main tunnels, branch tunnels and closed wall areas.

[0062] The ultrasonic transducer is controlled to transmit an ultrasonic transmission signal in the direction of the target lane. The receiving end captures the ultrasonic echo signal on the propagation path and verifies the validity of the received ultrasonic echo signal. Only the ultrasonic echo signals with a delay within the preset maximum propagation time range of the lane are screened as valid echo signals.

[0063] The preset maximum propagation time range of the tunnel is calculated based on the nominal propagation speed of ultrasound in mine air (340m / s) and the maximum straight-line distance of the tunnel (e.g., 1500m), with an additional 20% safety margin. The specific formula is: ;in, Represents the maximum propagation time range threshold, Represents the maximum straight-line distance of the roadway, Represents the nominal propagation velocity.

[0064] Obtain the absolute completion time of the measured effective echo signal in the opposite propagation direction. For example, the absolute completion time of the effective echo signal from the transmitting end A to the receiving end B is The absolute completion time of the effective echo signal from the transmitter B to the receiver A is .

[0065] The main frequency of the effective echo signal can be analyzed by fast Fourier transform The carrier frequency of the transmitted signal , according to the formula Calculate and obtain the frequency offset characteristic value .

[0066] S12. Bind the absolute completion time and frequency offset feature value of the same path and integrate them into an initial feature data set.

[0067] The initial feature data set includes the path number, transmitter, receiver, and frequency offset feature value. .

[0068] In some embodiments, in S2, based on the initial feature data set and the downhole propagation calibration value, a set of anti-dust propagation time difference values ​​is generated by fusion, including:

[0069] S21. Calculate the frequency offset equivalent time correction value based on the frequency offset characteristic value of each path in the initial characteristic data set, combined with the calibration value of the downhole ultrasonic propagation velocity; calculate the basic propagation time difference based on the absolute completion time in the opposite direction of each path in the initial characteristic data set.

[0070] Dust causes the sound wave energy to attenuate, making the actual propagation speed Lower than ideal (340m / s), quantify the speed deviation by the frequency offset characteristic value, and calculate the frequency offset equivalent time correction according to the linear mapping relationship : ; where k represents the mine tunnel attenuation coefficient, which can be determined by underground calibration tests, and the typical value is 0.018 .

[0071] According to the formula Calculate the difference in travel time in two opposite directions ,like , ,but .

[0072] S22. Superimpose the basic propagation time difference and the frequency offset equivalent time correction amount to generate a fused anti-interference propagation time difference.

[0073] like , representing the wind flow from A to B, then the fusion anti-interference propagation time difference ;like , representing the wind flow from B to A, then the fusion anti-interference propagation time difference , thereby eliminating the asymmetric interference of dust on the measurement with and against the wind.

[0074] S23. Generate a fused propagation time difference value set including the fused anti-interference propagation time differences of all paths.

[0075] The fused propagation time difference set includes the path number and the fused anti-interference propagation time difference .

[0076] In some embodiments, S3 analyzes the roadway space grid based on the anti-dust propagation time difference set to generate a resistance hot zone distribution map, including:

[0077] S31. Construct a tunnel spatial grid model based on the pre-built three-dimensional digital model of the tunnel.

[0078] Based on the pre-built three-dimensional digital model of the tunnel (including tunnel centerline coordinates, cross-sectional dimensions, and slope data), the physical space is discretized into structured grid units.

[0079] During discretization, reference points are generated along the tunnel centerline. The spacing between reference points is dynamically adjusted based on the tunnel shape. For example, the spacing between reference points on straight sections is 2 meters, while reference points are increased to 0.5 meters in curved areas (radius of curvature ≤ 30 meters). At the reference points, a cross-section perpendicular to the tunnel centerline is created, and mesh nodes are generated based on the cross-sectional dimensions.

[0080] Each grid node is associated with the position coordinates (x, y, z) and the ventilation cross-sectional area A. The grid topological connection relationship is stored as an adjacency matrix to obtain the tunnel space grid model.

[0081] S32. Utilize the fused propagation time difference set and the position of the ultrasonic path in the tunnel spatial grid model to solve the wind speed vector value at the grid node position.

[0082] Call the fused anti-interference propagation time difference in the fused propagation time difference set , combining the spatial topological relationship between the ultrasonic path and the tunnel space grid model, an overdetermined set of equations is established: ; The wind speed vector value at the grid node position is solved by the least squares method, which includes the magnitude and direction components.

[0083] S33. Calculate the spatial variation gradient of the wind speed vector values ​​of adjacent grid nodes.

[0084] Perform differential calculation on the wind speed vector values ​​of adjacent nodes in the tunnel space grid model, such as With downstream nodes , spatially varying gradient value ;in, Represents the modulus difference of the wind speed vector values ​​of adjacent nodes A and B, Represents the spatial distance measured in wind flow direction between adjacent nodes A and B.

[0085] S34. Mark the drag impact hotspot area and the baseline flow field reference area according to the spatial variation gradient value.

[0086] In some embodiments, marking the resistance impact hotspot area and the baseline flow field reference area according to the spatially varying gradient value in S34 includes:

[0087] S341. Set a first threshold and a second threshold.

[0088] The first threshold can be 1.8 times the average value of the spatial variation gradient values ​​of all nodes in the roadway space grid model, and the second threshold can be 0.3 times the average value of the spatial variation gradient values ​​of all nodes in the roadway space grid model.

[0089] S342. Identify the grid nodes whose spatial variation gradient values ​​exceed the first threshold and mark them as resistance impact hotspots; identify the grid nodes whose spatial variation gradient values ​​are lower than the second threshold and mark them as baseline flow field reference areas.

[0090] Scan all nodes, if the node space changes the gradient value If the node spatial gradient value exceeds the first threshold, it is marked as a resistance impact hotspot. If the value is lower than the second threshold, it is marked as the reference area of ​​the baseline flow field, and a hot zone distribution map is finally generated.

[0091] In some embodiments, in S4, the graded resistance state index is calculated by comparing the dynamic flow field of the resistance hot zone distribution map, including:

[0092] S41. Extract the wind speed vector values ​​of the grid nodes in the reference area of ​​the reference flow field and calculate the average reference wind speed vector value.

[0093] Extract the wind speed vector values ​​of all nodes in the reference area of ​​the benchmark flow field and mark them as ,in Represents the wind speed vector value of the mth node.

[0094] Calculate the magnitude of the average reference wind speed vector value : ; Among them, represents the node The modulus of the wind speed vector, m represents the number of grid nodes, Represents the wind speed vector value of the kth node, Represents the modulus of the wind speed vector value at the kth node.

[0095] Calculate the direction of the average reference wind speed vector value : ;in, Represents the direction angle of the wind speed vector value of the kth node, Represents the sum of the sine values ​​of the wind speed direction angles at all nodes, Represents the sum of the cosine values ​​of the wind speed direction angles at all nodes.

[0096] S42. For each grid node in the resistance-affecting hotspot, calculate the wind speed vector difference between the wind speed vector value and the average reference wind speed vector value; and calculate the dynamic resistance factor based on the wind speed vector difference and the node's ventilation cross-sectional area.

[0097] The wind speed vector difference includes the velocity component and directional components , , ;in, Represents the modulus of the wind speed vector value at the node, Represents the direction angle of the wind speed vector value at the node.

[0098] The dynamic resistance factor F is: ; Among them, K represents the mine air damping constant, A represents the ventilation cross-sectional area, and 90 is a dimensionless value representing the base ratio for angle conversion.

[0099] S43. Aggregate dynamic resistance factors to generate a comprehensive resistance state assessment value.

[0100] The comprehensive resistance state evaluation value R is: ; Where n represents the total number of nodes in the hotspot area, Represents the dynamic resistance factor of the i-th node.

[0101] S44. When the comprehensive resistance status evaluation value exceeds the preset warning threshold and the resistance impact hotspot area coincides with the gas sensor coordinates, the resistance-gas coupling warning mark is triggered.

[0102] Based on the historical resistance data of the mine under normal ventilation conditions, the maximum value of the comprehensive resistance state evaluation value R in the last 90 natural days is extracted, and the warning threshold is set as Set to 1.2 times the maximum value.

[0103] Based on the three-dimensional coordinate system of the mine, the Euclidean distance between the coordinates of the center point of the resistance impact hotspot area and the coordinates of the mine gas sensor is calculated. If the distance is less than or equal to the preset distance threshold, such as 5 meters, it is determined to be spatial overlap.

[0104] At the same time, the comprehensive resistance state assessment value is greater than the warning threshold If the distance between the hotspot area and the gas sensor is less than or equal to the distance threshold, the resistance-gas coupling warning mark will be activated.

[0105] In some embodiments, generating a spatial coordinate safety instruction based on the positioning result of the graded resistance state indicator in S5 includes:

[0106] S51. Locate the physical blockage point according to the spatial distribution coordinates of the resistance impact hotspot area and associate it with the roadway equipment database.

[0107] Based on the spatial distribution coordinate set of the resistance impact hotspot area, the geometric center point algorithm can be used to locate the physical blockage point.

[0108] Specifically, the three-dimensional coordinate data of all nodes in the resistance-affecting hotspot area are extracted, the arithmetic mean of these coordinates is calculated as the coordinate of the blockage center point, the coordinate of the center point is associated and matched with the tunnel equipment database, and the associated tunnel equipment, such as air doors, sensors, robot docking stations, etc., are scanned within a first specified range (such as a 15-meter radius).

[0109] S52. Set the third threshold. If the rate of change of the dynamic resistance factor reaches the third threshold and the damper coordinates are located, a damper opening adjustment instruction is generated; if the resistance-gas coupling warning mark is activated, a strong ventilation instruction + personnel evacuation instruction is generated; if the resistance-affected hotspot area persists for a specified time, a patrol robot cleaning instruction is generated.

[0110] The third threshold is set as the trigger benchmark for the change rate of the dynamic resistance factor, and the value of the third threshold is 80% of the historical maximum change rate.

[0111] Trigger three-level security instructions based on real-time status:

[0112] When the rate of change of the dynamic resistance factor exceeds a third threshold and the blockage point is within a second specified range (eg, a range of 5 meters) of the damper, a damper opening adjustment instruction is generated.

[0113] When the resistance-gas coupling warning mark is activated, strong ventilation instructions and personnel evacuation instructions are generated synchronously.

[0114] When the resistance hotspot persists for a preset time threshold (e.g., 2 hours), a cleaning instruction is generated for the inspection robot.

[0115] The command execution priority is: personnel evacuation command > air door adjustment command > robot cleaning command.

[0116] S53. Map the dynamic resistance factor to the three-dimensional digital model of the roadway to generate a real-time resistance heat map.

[0117] The dynamic resistance factor values ​​of all nodes in the tunnel spatial grid model are mapped to the three-dimensional digital model of the tunnel, and a real-time resistance heat map is generated through color-scale rendering.

[0118] For example, the area with a dynamic resistance factor less than 50N is displayed in dark blue and marked as a low resistance area; the area of ​​50-100N is displayed in yellow and marked as a medium resistance area; the area greater than 100N is displayed in red and marked as a high resistance area; at the same time, a flashing red star mark is marked at the center of the blockage.

[0119] In some embodiments, as Figure 2 As shown, the present application provides an ultrasonic precision wind measurement and resistance online monitoring system, which includes:

[0120] Original sound wave feature acquisition module: obtains the original sound wave features collected by the ultrasonic transducer array and integrates them to generate an initial feature data set.

[0121] Anti-dust propagation feature fusion module: Based on the initial feature data set and the downhole propagation calibration value, it fuses and generates an anti-dust propagation time difference set.

[0122] Tunnel space grid analysis module: Analyze the tunnel space grid based on the anti-dust propagation time difference set to generate the resistance hot zone distribution map.

[0123] Resistance state dynamic solution module: solves the graded resistance state index through dynamic flow field comparison of resistance hot zone distribution map.

[0124] Safety instruction coordinate generation module: Generates spatial coordinate safety instructions based on the positioning results of the graded resistance state indicators.

[0125] In some embodiments, the present application provides an ultrasonic precision wind measurement and online resistance monitoring device, which includes a memory and a processor; the memory is used to store computer programs; the processor is used to implement the steps of the ultrasonic precision wind measurement and online resistance monitoring method when executing the computer program.

[0126] In some embodiments, a readable storage medium stores computer program instructions. When the computer program instructions are read and executed by a processor, the steps of the ultrasonic precision wind measurement and resistance online monitoring method are executed.

[0127] Among them, any reference to memory, storage, database or other media used in the embodiments provided in this application may include non-volatile and / or volatile memory. Non-volatile memory may include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM) or flash memory. Volatile memory may include random access memory (RAM) or external cache memory.

[0128] It should be noted that, in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any actual relationship or order between these entities or operations. Moreover, the terms "comprises," "comprising," or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, article, or device that includes a series of elements includes not only those elements, but also other elements that are not explicitly listed, or elements that are inherent to such process, method, article, or device. In the absence of further limitations, an element defined by the phrase "comprising a ..." does not exclude the presence of other identical elements in the process, method, article, or device that includes the element.

[0129] Although the present application has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein; and these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present application.

Claims

1. An ultrasonic precise wind measurement and resistance online monitoring method, characterized in that: include: Obtaining the original sound wave features collected by the ultrasonic transducer array and integrating them to generate an initial feature data set; Based on the initial feature data set and the downhole propagation calibration value, the anti-dust propagation time difference set is generated by fusion; The roadway space grid is analyzed based on the set of resistance dust propagation time differences to generate a resistance hot zone distribution map; By comparing the dynamic flow field of the resistance hot zone distribution map, the graded resistance state index is calculated; Generate spatial coordinate safety instructions based on the positioning results of the graded resistance status indicators; The original sound wave features collected by the ultrasonic transducer array are obtained and integrated to generate an initial feature data set, including: Obtain the absolute completion time of the measured effective echo signal in the opposite propagation direction; calculate the absolute value of the frequency difference between the effective echo signal carrier frequency and the transmitted signal carrier frequency, and record it as the frequency offset characteristic value; the time delay of the effective echo signal is within the preset maximum propagation time range of the lane; Bind the absolute completion time and frequency offset feature value of the same path and integrate them into the initial feature data set; Based on the initial feature data set and the underground propagation calibration value, a set of anti-dust propagation time difference values ​​is generated by fusion, including: Calculate the frequency offset equivalent time correction value based on the frequency offset characteristic value of each path in the initial characteristic data set in combination with the downhole ultrasonic propagation velocity calibration value; calculate the basic propagation time difference based on the absolute completion time in the opposite direction of each path in the initial characteristic data set; The basic propagation time difference is superimposed with the frequency offset equivalent time correction to generate a fused anti-interference propagation time difference; Generate a fused propagation time difference value set including the fused interference-rejected propagation time difference values ​​of all paths.

2. The ultrasonic precise wind measurement and resistance online monitoring method according to claim 1 is characterized in that: The roadway spatial grid is analyzed based on the set of resistance dust propagation time differences to generate a resistance hot zone distribution map, including: Construct a tunnel spatial grid model based on the pre-built 3D digital model of the tunnel; Utilizing the fused propagation time difference set and combining the position of the ultrasonic path in the tunnel spatial grid model, the wind speed vector value at the grid node position is solved; Calculate the spatial variation gradient of wind speed vector values ​​at adjacent grid nodes; The resistance impact hotspot area and the benchmark flow field reference area are marked according to the spatial variation gradient value.

3. The ultrasonic precise wind measurement and resistance online monitoring method according to claim 2 is characterized in that: Mark the resistance impact hotspot area and the reference area of ​​the benchmark flow field according to the spatial variation gradient value, including: Setting a first threshold and a second threshold; Grid nodes whose spatial variation gradient values ​​exceed the first threshold are identified and marked as resistance impact hotspots; grid nodes whose spatial variation gradient values ​​are lower than the second threshold are identified and marked as baseline flow field reference areas.

4. The ultrasonic precise wind measurement and resistance online monitoring method according to claim 2 is characterized in that: By comparing the dynamic flow field of the resistance hot zone distribution map, the graded resistance state indicators are calculated, including: Extract the wind speed vector values ​​of the grid nodes in the reference area of ​​the benchmark flow field and calculate the average benchmark wind speed vector value; For each grid node in the resistance-affecting hotspot area: calculate the wind speed vector difference between the wind speed vector value and the average reference wind speed vector value; and calculate the dynamic resistance factor based on the wind speed vector difference and the node's ventilation cross-sectional area. Aggregate dynamic resistance factors to generate comprehensive resistance status assessment values; When the comprehensive resistance status evaluation value exceeds the preset warning threshold and the resistance impact hotspot coincides with the coordinates of the gas sensor, the resistance-gas coupling warning mark is triggered.

5. The ultrasonic precise wind measurement and resistance online monitoring method according to claim 4 is characterized in that: Based on the positioning results of the graded resistance status indicators, spatial coordinate safety instructions are generated, including: Locate physical blockage points based on the spatial distribution coordinates of resistance hotspots and associate them with the roadway equipment database; A third threshold is set. If the rate of change of the dynamic resistance factor reaches the third threshold and the damper coordinates are located, a damper opening adjustment instruction is generated. If the resistance-gas coupling warning flag is activated, a strong ventilation instruction and personnel evacuation instruction are generated. If the resistance-affected hotspot persists for a specified period of time, a patrol robot cleaning instruction is generated. The dynamic resistance factor is mapped to the three-dimensional digital model of the tunnel to generate a real-time resistance heat map.

6. An ultrasonic precision wind measurement and resistance online monitoring system, characterized in that: It includes: Original sound wave feature acquisition module: obtains the original sound wave features collected by the ultrasonic transducer array and integrates them to generate an initial feature data set; Anti-dust propagation feature fusion module: Based on the initial feature data set and the downhole propagation calibration value, it generates an anti-dust propagation time difference set; Tunnel space grid analysis module: Analyzes the tunnel space grid based on the anti-dust propagation time difference set and generates a resistance hot zone distribution map; Resistance state dynamic calculation module: calculates the graded resistance state index through dynamic flow field comparison of resistance hot zone distribution map; Safety instruction coordinate generation module: Generates spatial coordinate safety instructions based on the positioning results of the graded resistance state indicators; The original sound wave features collected by the ultrasonic transducer array are obtained and integrated to generate an initial feature data set, including: Obtain the absolute completion time of the measured effective echo signal in the opposite propagation direction; calculate the absolute value of the frequency difference between the effective echo signal carrier frequency and the transmitted signal carrier frequency, and record it as the frequency offset characteristic value; the time delay of the effective echo signal is within the preset maximum propagation time range of the lane; Bind the absolute completion time and frequency offset feature value of the same path and integrate them into the initial feature data set; Based on the initial feature data set and the underground propagation calibration value, a set of anti-dust propagation time difference values ​​is generated by fusion, including: Calculate the frequency offset equivalent time correction value based on the frequency offset characteristic value of each path in the initial characteristic data set in combination with the downhole ultrasonic propagation velocity calibration value; calculate the basic propagation time difference based on the absolute completion time in the opposite direction of each path in the initial characteristic data set; The basic propagation time difference is superimposed with the frequency offset equivalent time correction to generate a fused anti-interference propagation time difference; Generate a fused propagation time difference value set including the fused interference-rejected propagation time difference values ​​of all paths.

7. An ultrasonic precision wind measurement and resistance online monitoring device, characterized in that: It includes a memory and a processor; the memory is used to store computer programs; the processor is used to implement the steps of the ultrasonic precise wind measurement and resistance online monitoring method described in any one of claims 1 to 6 when executing the computer program.

8. A readable storage medium, characterized in that: The readable storage medium stores computer program instructions. When the computer program instructions are read and executed by a processor, the steps of the ultrasonic precise wind measurement and resistance online monitoring method according to any one of claims 1 to 6 are executed.

Images