A method and system for regulating and controlling explosion velocity of a gas explosion composition

Abstract

The present invention discloses a detonation velocity regulation and control method and system for gas explosion compositions, relating to the technical field of regulation and control. The method comprises: obtaining multiple experimental cycles, obtaining reference detonation wave arrival times, and forming a reference arrival time sequence; obtaining a stability index, and if the stability index exceeds a preset threshold, obtaining a theoretical detonation wave arrival time at a detection location and obtaining a reference average correction value; obtaining a target detonation velocity, and obtaining a target detonation wave arrival time at a detection location based on the target detonation velocity, and obtaining a target combined pressure based on the target detonation wave arrival time and the reference average correction value; conducting an experiment at the target combined pressure in the current experimental cycle, and obtaining the current detonation wave arrival time at the detection location in the current experimental cycle; obtaining a current detonation velocity based on the current detonation wave arrival time, and obtaining a regulation strategy based on the target detonation velocity and the current detonation velocity. The present invention has the advantages of precise regulation and control, reliable correction value determination, and dynamic environmental adaptation.

Description

Technical Field

[0001] The present invention relates to the field of regulation and control technology, and in particular to a method and system for regulating and controlling explosion velocity applied to a gas explosion composition. Background Art

[0002] In the field of industrial blasting, precise control of the detonation velocity of gas explosion compositions directly affects the blasting range and safety. Existing technologies mainly perform open-loop regulation by presetting the theoretical relationship between the combined gas pressure and the detonation velocity, that is, directly implementing the regulation after inferring the theoretical pressure value based on the target detonation velocity. In order to obtain a more accurate detonation velocity-gas pressure mapping relationship and improve the regulation and control accuracy of the implementation process, it is necessary to conduct multiple regulation and control experiments on the same ratio to improve the regulation and control. However, the existing detonation velocity regulation and control methods applied to gas explosion compositions have many defects.

[0003] First, due to dynamic factors such as fluctuations in the composition of explosive materials, changes in ambient temperature and humidity, and mechanical errors of equipment, the theoretical curve of detonation velocity-pressure based on ideal conditions deviates seriously from the actual response. Directly applying the theoretical value will cause the arrival time of the detonation wave to continue to shift, which not only fails to guarantee accuracy but also poses a safety hazard. Second, when obtaining the correction between the actual value and the theoretical value, the experimental stability is not taken into account. If forced adjustment is made under non-steady-state conditions such as uneven mixing of explosive materials or sensor abnormalities, the final adjustment control accuracy will be inaccurate due to invalid data. Multiple accumulated errors may even cause detonation velocity to lose control, which is typically manifested as the actual detonation velocity oscillating violently around the set value and failing to converge. Finally, the combined air pressure is directly implemented after a single correction based on the theoretical level without building a closed-loop adjustment logic. This will make it impossible to implement the verification cycle and will generate additional trial and error times for positioning the air pressure value. Summary of the Invention

[0004] In view of the defects in the prior art, the present invention provides a method and system for regulating and controlling explosion velocity applied to a gas explosion composition.

[0005] A detonation velocity regulation and control method for a gas explosion composition comprises: obtaining multiple experimental cycles under a preset combined gas pressure, and obtaining a reference detonation wave arrival time at a detection position in each experimental cycle, arranging the multiple reference detonation wave arrival times in sequence according to the order of the experimental cycles to form a reference arrival time sequence; obtaining a stability index based on the reference arrival time sequence, and if the stability index exceeds a preset threshold, obtaining a theoretical detonation wave arrival time at the detection position based on the preset combined gas pressure, and obtaining a reference average correction value based on the reference arrival time sequence and the theoretical detonation wave arrival time; obtaining a target detonation velocity, and obtaining a target detonation wave arrival time at the detection position based on the target detonation velocity, and obtaining a target combined gas pressure based on the target detonation wave arrival time and the reference average correction value; conducting an experiment at the target combined gas pressure in the current experimental cycle, and obtaining the current detonation wave arrival time at the detection position in the current experimental cycle; obtaining a current detonation velocity based on the current detonation wave arrival time, and obtaining a regulation strategy based on the target detonation velocity and the current detonation velocity.

[0006] Optionally, obtaining an adjustment strategy based on the target detonation velocity and the current detonation velocity includes: comparing the magnitude relationship between the target detonation velocity and the current detonation velocity, and determining the adjustment direction of the combined air pressure and the adjustment of the reference average correction amount based on the magnitude relationship; when the current detonation velocity is equal to the target detonation velocity, the adjustment strategy indicates that the target combined air pressure remains unchanged, and outputs a mapping relationship between the reference average correction amount and the target detonation velocity; when the current detonation velocity is less than the target detonation velocity, the adjustment strategy indicates an increase in the target combined air pressure and an increase in the reference average correction amount; when the current detonation velocity is greater than the target detonation velocity, the adjustment strategy indicates a decrease in the target combined air pressure and a decrease in the reference average correction amount.

[0007] Optionally, obtaining the theoretical detonation wave arrival time at the detection position according to the preset combined air pressure includes: obtaining the theoretical detonation velocity according to the preset combined air pressure; and obtaining the theoretical detonation wave arrival time at the detection position according to the theoretical detonation velocity.

[0008] Optionally, obtaining the target combined air pressure according to the target detonation wave arrival time and the reference average correction amount includes: obtaining a predicted detonation wave arrival time according to the target detonation wave arrival time and the reference average correction amount; obtaining a predicted detonation velocity according to the predicted detonation wave arrival time, and obtaining the target combined air pressure according to the predicted detonation velocity.

[0009] Optionally, obtaining the stability index based on the reference arrival time series includes: obtaining the average value of the arrival time of each reference detonation wave in the reference arrival time series; obtaining the difference between the arrival time of each reference detonation wave in the reference arrival time series and the average value, obtaining the sum of all the differences, and obtaining the stability index based on the sum of all the differences.

[0010] Optionally, obtaining a reference average correction value based on the reference arrival time sequence and the theoretical detonation wave arrival time includes: obtaining a difference between each reference detonation wave arrival time in the reference arrival time sequence and the theoretical detonation wave arrival time, and obtaining the reference average correction value based on the multiple differences.

[0011] Also provided is a detonation velocity regulation and control system for a gas explosion composition. The system is used to implement a detonation velocity regulation and control method for a gas explosion composition. The system includes: an acquisition module for acquiring multiple experimental cycles under a preset combined gas pressure, and obtaining a reference detonation wave arrival time at a detection position in each experimental cycle, arranging the multiple reference detonation wave arrival times in sequence according to the experimental cycles to form a reference arrival time sequence; a data processing module for obtaining a stability index based on the reference arrival time sequence, and if the stability index exceeds a preset threshold, obtaining a theoretical detonation wave arrival time at the detection position based on the preset combined gas pressure, and obtaining a reference average correction value based on the reference arrival time sequence and the theoretical detonation wave arrival time; a first regulation and control module for obtaining a target detonation velocity, obtaining a target detonation wave arrival time at the detection position based on the target detonation velocity, obtaining a target combined gas pressure based on the target detonation wave arrival time and the reference average correction value, conducting an experiment at the target combined gas pressure in the current experimental cycle, and obtaining the current detonation wave arrival time at the detection position in the current experimental cycle; and a second regulation and control module for obtaining a current detonation velocity based on the current detonation wave arrival time, and obtaining a regulation strategy based on the target detonation velocity and the current detonation velocity.

[0012] Optionally, the second regulation control module is also used to: compare the size relationship between the target detonation velocity and the current detonation velocity, and determine the adjustment direction of the combined air pressure and the adjustment of the reference average correction amount based on the size relationship; when the current detonation velocity is equal to the target detonation velocity, the regulation strategy indicates that the target combined air pressure remains unchanged, and outputs a mapping relationship between the reference average correction amount and the target detonation velocity; when the current detonation velocity is less than the target detonation velocity, the regulation strategy indicates an increase in the target combined air pressure and an increase in the reference average correction amount; when the current detonation velocity is greater than the target detonation velocity, the regulation strategy indicates a decrease in the target combined air pressure and a decrease in the reference average correction amount.

[0013] Optionally, the data processing module is further used to: obtain a theoretical detonation velocity according to a preset combined gas pressure; and obtain a theoretical detonation wave arrival time at the detection position according to the theoretical detonation velocity.

[0014] Optionally, the data processing module is further used to: obtain a predicted detonation wave arrival time based on the target detonation wave arrival time and a reference average correction amount; obtain a predicted detonation velocity based on the predicted detonation wave arrival time; and obtain a target combined air pressure based on the predicted detonation velocity.

[0015] The beneficial effects of the present invention are embodied in:

[0016] In the entire explosion velocity regulation and control method applied to gas explosion compositions, infrared radiation signals are first sent at different times, and a calculation model is constructed in combination with key parameters such as surface emissivity, surface temperature, and downward radiation brightness temperature of the atmosphere, thereby fully considering the complexity of water vapor distribution and changes in the atmosphere, providing a data basis for subsequent accurate acquisition of atmospheric transmittance and water vapor content, and at the same time improving the accuracy of water vapor detection data; further, by introducing a correction model, the water vapor content obtained from the preliminary detection is dynamically adjusted and corrected based on real-time meteorological conditions and meteorological influencing factors during the detection time period, such as wind speed differences and air pressure differences. This correction mechanism can effectively eliminate the potential impact of meteorological differences on the detection results, making the final water vapor content data more accurate and reliable; further, the entire method overcomes the limitations of traditional satellite remote sensing technology that is easily interfered with by surface characteristics and atmospheric conditions, and fully improves the stability and consistency of the detection data through refined data processing procedures and optimization algorithms. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] To more clearly illustrate the specific embodiments of the present invention or the technical solutions in the prior art, the following briefly describes the drawings required for the specific embodiments or the description of the prior art. Similar elements or parts are generally identified by similar reference numerals throughout the drawings. Elements or parts in the drawings are not necessarily drawn to scale.

[0018] Figure 1 Schematic diagram of part of the flow chart of S1 and S2 in the detonation velocity regulation and control method applied to a gas explosion composition according to the present invention;

[0019] Figure 2 Schematic diagram of a portion of the flow chart of S3 in the detonation velocity regulation and control method for a gas explosion composition according to the present invention;

[0020] Figure 3 Schematic diagram of a portion of the flow chart of S4 in the detonation velocity regulation and control method for a gas explosion composition according to the present invention;

[0021] Figure 4 Schematic diagram of the steps of the detonation velocity regulation and control method applied to a gas explosion composition according to the present invention;

[0022] Figure 5 Schematic diagram of some steps of the detonation velocity adjustment control method S4 applied to a gas explosion composition according to the present invention;

[0023] Figure 6 Schematic diagram of some steps of the detonation velocity adjustment control method S2 applied to a gas explosion composition according to the present invention;

[0024] Figure 7 Schematic diagram of some steps of the detonation velocity regulation control method S3 applied to a gas explosion composition according to the present invention. DETAILED DESCRIPTION

[0025] To make the objectives, technical solutions, and advantages of the embodiments of the present invention more clear, the technical solutions of the embodiments of the present invention will be clearly and completely described below in conjunction with the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Generally, the components of the embodiments of the present invention described and shown in the drawings herein can be arranged and designed in various different configurations.

[0026] Therefore, the following detailed description of the embodiments of the present invention provided in the accompanying drawings is not intended to limit the scope of the invention as claimed, but rather merely represents selected embodiments of the present invention. All other embodiments derived by persons of ordinary skill in the art based on the embodiments of the present invention without creative effort are intended to fall within the scope of protection of the present invention.

[0027] It should be noted that similar reference numerals and letters represent similar items in the following drawings. Therefore, once an item is defined in one drawing, it does not need to be further defined or explained in subsequent drawings. In addition, the terms "first," "second," etc. are used only to distinguish the descriptions and are not to be understood as indicating or implying relative importance.

[0028] like Figure 1 、 Figure 2 、 Figure 3 and Figure 4 As shown, a method for regulating and controlling the detonation velocity of a gas explosion composition is provided, comprising:

[0029] S1. Acquire multiple experimental cycles under a preset combined gas pressure, and obtain a reference detonation wave arrival time at a detection position in each experimental cycle, and arrange the multiple reference detonation wave arrival times in sequence according to the order of the experimental cycles to form a reference arrival time sequence;

[0030] S2. Obtaining a stability index based on a reference arrival time series. If the stability index exceeds a preset threshold, obtaining a theoretical detonation wave arrival time at the detection position based on a preset combined air pressure, and obtaining a reference average correction value based on the reference arrival time series and the theoretical detonation wave arrival time.

[0031] S3. Obtain a target detonation velocity, and obtain a target detonation wave arrival time at the detection position based on the target detonation velocity, and obtain a target combined air pressure based on the target detonation wave arrival time and a reference average correction amount, and perform an experiment according to the target combined air pressure in the current experimental cycle, and obtain the current detonation wave arrival time at the detection position in the current experimental cycle;

[0032] S4. Obtain the current detonation velocity according to the current detonation wave arrival time and the reference average correction amount, and obtain an adjustment strategy according to the target detonation velocity and the current detonation velocity.

[0033] In this embodiment, it should be noted that in S1, a time series baseline layer of experimental data is established, which is a prerequisite for achieving closed-loop control of detonation velocity. This step first requires executing a multi-cycle repeated experiment (e.g., 10 consecutive detonations) at a fixed, preset combined gas pressure. Sensors are then used to capture the physical arrival time of each detonation wave at a preset detection location (e.g., 3 meters from the detonation point). This raw data is referred to as the "reference detonation wave arrival time," and its essence reflects the dynamic characteristics of the actual propagation of the detonation wave at the current pressure. Subsequently, each arrival time is arranged into a serialized structure in strict experimental order, forming a "reference arrival time sequence." This time-series storage preserves the continuity of the experimental process, allowing dynamic fluctuation patterns (such as a trend of decreasing arrival time due to gradual improvement in gas mixing uniformity) or sudden disturbances (such as a single sensor false trigger) to be effectively identified in subsequent analysis.

[0034] Taking a mining blasting project as an example, when technicians conduct five test explosions at a fixed gas pressure of 0.5 MPa, step S1 requires recording the specific time at which each detonation wave reaches a sensor within the rock mass (e.g., 53.2 ms for the first explosion, 52.8 ms for the second explosion, 53.1 ms for the third explosion, 52.9 ms for the fourth explosion, and 53.0 ms for the fifth explosion). Once these data are sequenced, the stability information they imply is crucial: small fluctuations in the sequence (such as the data between 52.8 and 53.2 ms) indicate that the experimental conditions are under control; however, large jumps (such as a sudden jump to 60.1 ms for the third explosion) suggest a gas mixing failure or sensor malfunction during a particular experiment. This time series structure provides direct input for subsequent stability quantification in S2 (e.g., calculation of the dispersion index), ensuring that corrections are generated based solely on a reliable data set—avoiding distortion of the regulatory baseline due to single outliers or unsteady conditions (e.g., insufficiently stirred explosives).

[0035] In S2, a reliable correction benchmark is established by dynamically evaluating the stability of the experimental data. This step first quantifies the volatility of the reference arrival time series generated in S1. This involves calculating the deviation of multiple arrival times in the series from their average value, generating a quantitative stability index. This index is compared with a preset threshold. Only when the stability requirement is met (e.g., the index exceeds the preset threshold) is the correction phase permitted. This effectively creates a quality filter at the data level, eliminating invalid data interference caused by uneven gas mixing, transient sensor failures, or sudden environmental changes. Once the data is confirmed to be stable, a theoretical model is used to calculate the theoretical detonation wave arrival time at that location based on a preset combined pressure. This theoretical value is then compared with the statistical characteristics of the actual series (e.g., the historical average arrival time). Ultimately, a reference average correction value is generated that reflects the long-term deviation trend. It should be noted that the preset threshold for the stability index can be determined through experimental calibration, dynamic adaptive mechanisms, or engineering recommendations. For example, an experimental calibration process involves repeated detonation of a fixed ratio in a constant temperature and humidity laboratory. The time series of the detonation wave arrival is recorded, the deviation sum is calculated, and the upper limit of the 95% confidence interval is taken to determine the corresponding preset threshold. For example, a dynamic adaptive mechanism is employed. For high-precision blasting scenarios, the preset threshold is 0.85; for strong interference environments, the preset threshold is 0.6. For continuous loss of control scenarios, the preset threshold is self-learned, automatically relaxing the conditions (automatically reducing the threshold by a factor of 0.9). For example, engineering recommendations are used, such as a preset threshold of 0.75 for open-pit mine blasting, 0.85 for high-precision underground tunnels, and 0.65 for fractured rock formations.

[0036] Taking a tunnel blasting project as an example, consider a series of five experiments showing the detonation wave arrival time sequence of [20.5ms, 20.7ms, 20.3ms, 19.8ms, 20.4ms]. A dispersion analysis confirms that the fluctuations are within an acceptable range (the stability index meets the required standard). In this case, a theoretical model based on the air pressure values ​​calculates a theoretical arrival time of 19.2ms. By comparing the sequence mean of 20.34ms with the theoretical value of 19.2ms, a reference average correction of +1.14ms is determined (the actual delay is slower than the theoretical value). This correction essentially captures inherent errors under the current experimental conditions that cannot be accounted for by the theoretical model, such as the retarding effect of the tunnel's unique rock formations on the detonation velocity. Conversely, if the series exhibits abnormal fluctuations of [20.5ms, 35.1ms, 20.3ms] (perhaps due to a fuse delay failure), the correction is suspended due to insufficient stability indicators, and the experiment must be repeated until the data stabilizes. This mechanism fundamentally prevents occasional abnormal data from contaminating the correction reference database, ensuring the accuracy of subsequent adjustments.

[0037] In S3, the historical correction amount is used to synchronously calibrate the mapping relationship between the target reference and the air pressure. This step first converts the target detonation velocity into the theoretical detonation wave arrival time at the detection position, and then introduces the reference average correction amount generated by S2 for the first correction (for example: if the theoretical time is 50ms and the correction amount is +2ms, the predicted actual time is 52ms). This correction essentially compensates for the systematic deviations that the theoretical model cannot cover (such as inherent equipment delays or material batch differences), so that the target setting value fits the actual physical scenario. Subsequently, the required combined air pressure is inferred based on the corrected predicted arrival time - through the detonation velocity-air pressure mapping relationship library, the mapping relationship library can be obtained by a limited number of experiments on commonly used gas explosion combinations (CH4:O2:N2). Generally speaking, , where D is the detonation velocity, P is the air pressure, and It is the empirical coefficient under the predetermined ratio obtained by limited experiments in the prior art. For example, when the predetermined ratio of CH4:O2:N2 is 1:2:3, =2000, =0.5, In summary, the time parameter is converted into an accurate combined pressure value. Finally, the current cycle experiment is executed at this pressure value, and the detonation wave arrival time is collected in real time. This data will serve as a key input for closed-loop verification.

[0038] Furthermore, for example, consider the need to increase detonation velocity during tunnel blasting. If the target velocity requires the detonation wave to arrive within 40 ms at a distance of 30 meters, the system first calculates a theoretical arrival time of 38 ms (ignoring corrections). However, combined with the reference average correction obtained by S2 (e.g., +3 ms, reflecting the retarding effect of the tunnel rock on detonation velocity), the system predicts that it will actually take 41 ms to reach the target velocity. Based on this prediction, a query of the air pressure mapping database reveals that to achieve a 41 ms arrival time, the air pressure must be increased from 0.8 MPa to 1.0 MPa. After an actual experiment at 1.0 MPa, the arrival time was recorded as 40.8 ms. This process validates the effectiveness of the dual correction method: the initial theoretical value (38 ms) is corrected to approach the actual value (41 ms). The resulting air pressure (1.0 MPa) results in an actual arrival time (40.8 ms) that deviates only 0.2% from the predicted time (41 ms), far exceeding the error of the uncorrected solution and significantly improving the success rate of one-shot adjustments.

[0039] In S4, the detonation wave arrival time and detection position collected during the current experimental cycle are first used to recover the actual detonation velocity (current detonation velocity). This calculated current detonation velocity is then directly compared with the target detonation velocity to generate a three-way adjustment strategy (increase / decrease / maintain air pressure), forming the core feedback channel for closed-loop control. This allows for a precise mapping between the output reference average correction and the target detonation velocity, adaptive to the current scenario, ultimately achieving precise detonation velocity control.

[0040] Taking the later stages of tunnel blasting regulation as an example, let's assume that during the S3 phase, the measured arrival time of the detonation wave at 1.0 MPa is 40.8 ms, resulting in a calculated current detonation velocity of 783 m / s. If the target detonation velocity is 800 m / s, the current value is considered too low. The strategy then issues a "pressurization" command (rather than directly applying the theoretical value): in the next cycle, the pressure is raised to 1.05 MPa. After repeating S3-S4, the newly measured arrival time is 39.6 ms, and the detonation velocity rises to 809 m / s. Because the current value (809 m / s) exceeds the target (800 m / s), a "depressurization" command is immediately issued in the opposite direction. This bidirectional oscillation convergence mechanism (increase pressure → slightly exceed target → depressurize), combined with the continuous purification of time series data by the correction factor, ensures that the detonation velocity quickly stabilizes within the target range within 3-4 cycles, completely avoiding the persistent deviation or oscillation divergence caused by uncorrected measured values ​​in traditional open-loop regulation.

[0041] In summary, the entire detonation velocity regulation and control method applied to gas explosion compositions improves the accuracy and reliability of detonation velocity control of gas explosion compositions; among them, dynamic quality filtering of experimental data is achieved by introducing stability quantification indicators (such as time series discreteness analysis), eliminating invalid data interference under non-steady-state conditions (uneven mixing, sensor failure) from the root, ensuring that the generated reference average correction is based only on repeatable physical laws (such as the blocking effect of fixed rock formations), avoiding the correction benchmark drift caused by occasional outliers in traditional methods; further, the reference average correction is incorporated into the target detonation wave time prediction to improve the physical matching of the gas pressure setting, and after solving the current detonation velocity, the strategic decision (increase / decrease) always responds to the actual detonation velocity change rather than measurement noise; further, a negative feedback loop is formed based on the adjustment strategy of the current and target detonation velocities, and combined with the continuous purification of the data, the detonation velocity is quickly stabilized in the target range (such as convergence within a small number of cycles).

[0042] like Figure 3 and Figure 5 As shown, in one embodiment, obtaining the adjustment strategy according to the target detonation velocity and the current detonation velocity in S4 includes:

[0043] S41, comparing the magnitude relationship between the target detonation velocity and the current detonation velocity, and determining an adjustment direction of the combined gas pressure and an adjustment of a reference average correction amount based on the magnitude relationship;

[0044] S42, when the current detonation velocity is equal to the target detonation velocity, the adjustment strategy indicates that the target combined air pressure remains unchanged, and outputs a mapping relationship between a reference average correction amount and the target detonation velocity;

[0045] S43, when the current detonation velocity is less than the target detonation velocity, the adjustment strategy instructs to increase the target combined air pressure and increase the reference average correction amount;

[0046] S44: When the current detonation velocity is greater than the target detonation velocity, the adjustment strategy instructs to reduce the target combined air pressure and lower the reference average correction amount.

[0047] In this embodiment, it should be noted that in S41, the adjustment decision is triggered by real-time calculation of the purified detonation velocity deviation. In other words, the current detonation velocity obtained by the secondary correction in S4 is directly compared with the target detonation velocity to generate a discretized magnitude relationship signal (greater than / equal to / less than).

[0048] In S42, when the current detonation velocity is equal to the target value, an air pressure maintenance command is output, and the mapping relationship between the reference average correction value and the target detonation velocity is output. The significance of engineering safety lies in maintaining the air pressure unchanged after entering the steady state, preventing redundant adjustments from introducing new disturbances, achieving accurate mapping relationship determination, and thus improving subsequent control accuracy. A typical scenario is underground tunnel blasting. When the detonation velocity stabilizes to the target value (2800m / s) or thereabouts, the S42 command can avoid the disruption of rock stress balance due to fine-tuning of air pressure, eliminating the risk of inducing rock bursts; it is locked only when the purified detonation velocity is strictly equal to the target value, avoiding critical oscillations caused by fuzzy thresholds.

[0049] In S43, a boost command is triggered when the current detonation velocity is less than the target value. This can be done based on the detonation velocity-pressure mapping library described in the above-mentioned embodiment, or based on an established detonation velocity-pressure gradient relationship (e.g., every 0.1 MPa corresponds to a 50 m / s detonation velocity increment). The optimal boost amplitude can be calculated. For example, in submarine blasting, if a sudden drop in water temperature causes a 5% decrease in detonation velocity, a gradient boost of 0.12 MPa can be used to accurately compensate for environmental interference. The boost amplitude is constrained by a stability indicator: if the data fluctuates significantly (e.g., due to uneven gas mixing), the adjustment step size is automatically reduced to avoid sudden pressure changes that could cause detonation instability.

[0050] In S44, a "decompression" command is triggered when the current detonation velocity exceeds the target value. When rock fractures expand, causing a sudden increase in detonation velocity (e.g., current 2700 m / s, target 2500 m / s), S44 immediately issues a decompression command, reducing the pressure (e.g., by 0.15 MPa) based on the gradient relationship to bring the detonation velocity back to the target range for the next cycle. In conjunction with time series data stability analysis, if there are continuous oscillations from "overshoot to decompression to undershoot," the decompression range is automatically increased (e.g., from 0.15 MPa to 0.20 MPa) to accelerate convergence to a steady state.

[0051] like Figure 1 and Figure 6 As shown, in one embodiment, obtaining the theoretical detonation wave arrival time at the detection position according to the preset combined gas pressure in S2 includes:

[0052] S21. Obtaining a theoretical detonation velocity based on a preset combined air pressure;

[0053] S22. Obtain a theoretical detonation wave arrival time at the detection position according to the theoretical detonation velocity.

[0054] In this embodiment, it should be noted that in S21, the detonation velocity-pressure mapping relationship library is called based on the preset combined pressure to obtain the corresponding detonation velocity value. The core mapping relationship depends on a finite number of experimental formulas ( Its engineering significance lies in that, for example, in a methane, oxygen, and nitrogen mixture (fixed ratio 1:2:3), the increase in gas pressure will increase the frequency of molecular collisions and accelerate the propagation of the detonation reaction. The theoretical model is calibrated by the parameters and Quantify this effect; when the ambient humidity changes suddenly (such as water seepage in a tunnel), the actual reactivity of the gas explosive deviates from the ideal state. However, because the theoretical model does not incorporate adaptive environmental variables, the output explosion velocity value (such as 2500m / s) may have a systematic deviation from the physical reality (actually only 2400m / s).

[0055] In S22, the theoretical arrival time of the detonation wave at the detection location is calculated based on the propagation of the detonation wave over a fixed distance (e.g., time = distance / detonation velocity). Its core flaw is the lack of modeling of interference factors. For example, if the detection point is 30 meters from the detonation point and the theoretical detonation velocity is 2500 m / s, actual rock fractures or friction in the charge tube wall can delay the propagation of the detonation wave (e.g., the actual arrival time is 13.5 ms). This deviation is not represented in the theoretical model and its relevance to subsequent corrections. The reference average correction generated by S2 is precisely to compensate for these unmodeled factors, allowing subsequent S3 and S4 to approximate the actual physical response through dynamic correction.

[0056] Further, for example, in mine blasting, if based on the same gas explosion composition ratio, when the preset pressure is 3MPa, S21 calls the pre-stored empirical coefficient =1950, =0.48, outputting a theoretical detonation velocity of 2350 m / s. Based on this, S22 calculates the theoretical arrival time for a 30-meter process to be 12.77 ms. However, due to fluctuations in methane concentration within the cave, the actual propagation time of the detonation wave reaches 14.1 ms. This error is captured by S2's reference average correction and used to calibrate the target pressure setpoint in S3, ultimately improving control accuracy.

[0057] like Figure 2 and Figure 7 As shown, in one embodiment, obtaining the target combined pressure according to the target detonation wave arrival time and the reference average correction amount in S3 includes:

[0058] S31, obtaining a predicted detonation wave arrival time according to the target detonation wave arrival time and a reference average correction value;

[0059] S32. Obtain a predicted detonation velocity based on the predicted detonation wave arrival time, and obtain a target combined air pressure based on the predicted detonation velocity.

[0060] In this embodiment, it should be noted that in S31, the target detonation wave arrival time (derived from the target detonation velocity theory) is superimposed on a reference average correction (e.g., +3ms) to generate a predicted detonation wave arrival time. This operation essentially translates the theoretical value into physical reality, compensating for irreducible factors such as device delay and environmental resistance. For example, in a mineral blasting operation, if the theoretical target time is 40ms, but the historical correction indicates that the drilling cuttings resistance effect requires +10ms compensation, this method differs from existing methods in that it avoids directly applying the theoretical time (40ms), which would result in an excessively low air pressure setting and could lead to engineering accidents caused by a severe detonation velocity deficiency.

[0061] In S32, based on the predicted detonation wave arrival time from S31, control parameters are inferred from a pre-built detonation velocity-pressure relationship library. The predicted detonation velocity is calculated using the ratio of distance to predicted time (e.g., 30 meters / 44 milliseconds = 681 m / s). The target combined pressure (e.g., 681 m / s corresponds to 0.2 MPa) is then retrieved from the detonation velocity-pressure mapping library. For example, in a submarine blasting scenario, if the predicted detonation wave arrival time is affected by currents, a reference average correction of 4 milliseconds is applied, reducing the target detonation velocity to the predicted detonation velocity, e.g., 650 m / s. This predicted detonation velocity is then used to determine the precise target combined pressure, e.g., 1.15 MPa, and subsequently control the system. This approach, rather than the existing open-loop method that directly applies theoretical values ​​(a target detonation velocity of 700 m / s corresponds to a combined pressure of 1.25 MPa), addresses the existing failure to consider the reference average correction corresponding to environmental impacts, thereby improving control accuracy.

[0062] In one embodiment, obtaining the stability index according to the reference in-place time series in S2 includes:

[0063] Obtain the average value of each reference detonation wave arrival time in the reference arrival time series;

[0064] Obtaining the difference between each reference detonation wave arrival time in the reference arrival time series and the average value, and obtaining the sum of all the differences;

[0065] The stability index is obtained based on the sum of all differences.

[0066] It should be noted that the stability index obtained from the reference in-place time series in S2 can be obtained using the following expression, specifically expressed as: ;in,

[0067] is a stability indicator, is the number of reference detonation wave arrival times in the reference arrival time series, is the i-th reference detonation wave arrival time in the reference arrival time series, is the jth reference detonation wave arrival time in the reference arrival time series.

[0068] In this embodiment, it should be noted that Absolute deviation accumulation, used to calculate the arrival time of each reference detonation wave and the series mean Example: If the sequence is [20.0, 20.2, 20.3, 19.8] (unit: ms), the mean , the absolute deviation sum is calculated as: |20.0-20.075|+|20.2-20.075|+|20.3-20.075|+|19.8-20.075|=0.75. This quantifies volatility and directly characterizes the degree of dispersion in experimental data (small fluctuations, such as ±0.2ms, have a deviation sum of ≈0.5; large jumps, such as a 15ms spike, have a deviation sum greater than 10). It also allows for the identification of abnormal operating conditions: uneven mixing of explosive materials or sensor failure can cause a sharp shift in single-point data. In these cases, the deviation sum increases significantly, triggering stability alerts and other processes.

[0069] Furthermore, exp(-sum of all differences) is a negative exponential function transformation, mapping the absolute deviation sum to the interval (0, 1). When the deviation sum → 0 (ideal stability), exp(-sum of all differences) = 1; when the deviation sum → ∞ (severe instability), exp(-sum of all differences) approaches 0. The exponential function also implements nonlinear response characteristics, making it sensitive to small deviations. When the deviation sum increases from 0.1 to 0.3, the stability index decreases from 0.90 to 0.74 (an 18% decrease), forcing high data consistency (e.g., ±0.1ms). It also strongly suppresses large deviations, with the stability index approaching 0 when the deviation sum is greater than 2, completely shielding against interference from single faults (e.g., momentary sensor failure). In summary, the combination of linear absolute deviation and exponential compression achieves linear amplification of small deviations, ensuring reproducibility across multiple experiments. Large deviations return the exponential to zero, effectively blocking invalid data with a single click.

[0070] Obtaining stability indicators based on the reference in-place time series, converting the experimental stability quantification problem into programmable discrimination logic, and providing a data source for subsequent dynamic corrections are the core algorithm basis for overcoming problems such as detonation velocity loss of control and oscillation divergence.

[0071] In one embodiment, obtaining a reference average correction value according to the reference arrival time sequence and the theoretical detonation wave arrival time in S2 includes:

[0072] Obtain the difference between the reference detonation wave arrival time and the theoretical detonation wave arrival time in the reference arrival time series;

[0073] The sum of all the differences in this embodiment is obtained, and then the average value is calculated to form a reference average correction value.

[0074] It should be noted that the reference average correction amount obtained in S2 based on the reference arrival time series and the theoretical detonation wave arrival time can be obtained using the following expression, which is specifically expressed as: ;in,

[0075] is the reference average correction amount, is the number of reference detonation wave arrival times in the reference arrival time series, is the kth reference detonation wave arrival time in the reference arrival time series, is the theoretical detonation wave arrival time.

[0076] In this embodiment, it should be noted that, in the entire expression, Represents the deviation directional capture, calculating the reference detonation wave arrival time for each experiment minus the theoretical detonation wave arrival time Example: Assume that the theoretical time under a certain air pressure is =19.2ms, 5 measured sequences are [20.5, 20.7, 20.3, 19.8, 20.4]ms, then the single deviations are: [1.3ms, 1.5ms, 1.1ms, 0.6ms, 1.2ms]. Quantify the fixed deviation (such as rock formation blockage) and non-random noise that are not covered by the theoretical model in a single experiment; eliminate the single correction distortion defect. The existing method uses the single deviation data to correct the subsequent settings, such as taking However, if the gas explosion mixture fluctuates, the value will be invalid (such as The exception is 35ms, resulting in a deviation of 15.8ms).

[0077] Further, It represents the arithmetic mean, taking the average of n deviations to generate a long-term statistical reference average correction value. This suppresses random disturbances. If a sensor failure causes a relatively large deviation, the abnormal value will be diluted by other valid data in the mean calculation to avoid contamination of the correction benchmark; strengthen anti-interference ability: ambient temperature and humidity fluctuations can cause a single deviation of 0.5ms, but the average of n=5 times will compress the random noise to 0.22ms (standard deviation theory), making Only systematic physical deviations (such as rock block) are retained.

[0078] In summary, the core error in detonation velocity control stems from fixed propagation delays (such as equipment response lag + rock mass blockage) and non-proportional distortion. The entire expression is implemented through a three-order design consisting of single-point deviation directional capture → statistical mean filtering → linear physical superposition, quantifying the systematic deviation between the theoretical model and physical reality into a programmable constant correction. This ultimately solves the problem of detonation velocity oscillation caused by inaccurate corrections, providing an algorithmic basis for high-precision detonation velocity control.

[0079] Also provided is a detonation velocity regulation and control system for a gas explosion composition, the system being used to implement the detonation velocity regulation and control method for a gas explosion composition in any one of the above embodiments, the system comprising:

[0080] An acquisition module is used to acquire multiple experimental cycles under a preset combined gas pressure, and obtain a reference detonation wave arrival time at the detection position in each experimental cycle, and arrange the multiple reference detonation wave arrival times in sequence according to the order of the experimental cycles to form a reference arrival time sequence;

[0081] a data processing module for obtaining a stability index based on a reference arrival time series; if the stability index exceeds a preset threshold, obtaining a theoretical detonation wave arrival time at the detection position based on a preset combined gas pressure; and obtaining a reference average correction value based on the reference arrival time series and the theoretical detonation wave arrival time;

[0082] a first regulating and controlling module, configured to obtain a target detonation velocity, obtain a target detonation wave arrival time at a detection position based on the target detonation velocity, obtain a target combined air pressure based on the target detonation wave arrival time and a reference average correction amount, conduct an experiment according to the target combined air pressure in a current experimental cycle, and obtain a current detonation wave arrival time at the detection position in the current experimental cycle;

[0083] The second regulation control module is used to obtain the current detonation velocity according to the current detonation wave arrival time, and obtain the regulation strategy according to the target detonation velocity and the current detonation velocity.

[0084] In one embodiment, the second regulation control module is further used to: compare the magnitude relationship between the target detonation velocity and the current detonation velocity, and determine the adjustment direction of the combined air pressure and the adjustment of the reference average correction amount based on the magnitude relationship; when the current detonation velocity is equal to the target detonation velocity, the regulation strategy indicates that the target combined air pressure remains unchanged, and outputs a mapping relationship between the reference average correction amount and the target detonation velocity; when the current detonation velocity is less than the target detonation velocity, the regulation strategy indicates an increase in the target combined air pressure and an increase in the reference average correction amount; when the current detonation velocity is greater than the target detonation velocity, the regulation strategy indicates a decrease in the target combined air pressure and a decrease in the reference average correction amount.

[0085] In one embodiment, the data processing module is further configured to: obtain a theoretical detonation velocity according to a preset combined gas pressure; and obtain a theoretical detonation wave arrival time at the detection position according to the theoretical detonation velocity.

[0086] In one embodiment, the first regulation control module is further configured to: obtain a predicted detonation wave arrival time based on the target detonation wave arrival time and the reference average correction amount; obtain a predicted detonation velocity based on the predicted detonation wave arrival time; and obtain a target combined gas pressure based on the predicted detonation velocity.

[0087] In this embodiment, it should be noted that, regarding the above-mentioned explosion velocity regulation and control system applied to the gas explosion composition, the specific manner of performing the operation has been described in detail in the embodiment of the explosion velocity regulation and control method applied to the gas explosion composition, and will not be elaborated here.

[0088] The preferred embodiments of the present invention are described in detail above in conjunction with the accompanying drawings. However, the present invention is not limited to the specific details in the above embodiments. Within the technical concept of the present invention, various simple modifications can be made to the technical solution of the present invention, and these simple modifications all fall within the scope of protection of the present invention.

[0089] It should also be noted that the various specific technical features described in the above specific embodiments can be combined in any appropriate manner without contradiction. To avoid unnecessary repetition, the present invention will not further describe various possible combinations.

[0090] In addition, the various embodiments of the present invention may be arbitrarily combined, and as long as they do not violate the concept of the present invention, they should also be regarded as the contents disclosed by the present invention.

[0091] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the above embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the above embodiments, or make equivalent replacements for some or all of the technical features therein. These modifications or replacements do not deviate the essence of the corresponding technical solutions from the scope of the technical solutions of the embodiments of the present invention, and they should all be included in the scope of the claims and description of the present invention.

Claims

1. A method for regulating and controlling the detonation velocity of a gas explosion composition, characterized in that: include: Acquire multiple experimental cycles under a preset combined gas pressure, and obtain a reference detonation wave arrival time at the detection position in each experimental cycle, and arrange the multiple reference detonation wave arrival times in sequence according to the order of the experimental cycles to form a reference arrival time sequence; A stability index is obtained based on a reference arrival time series. If the stability index exceeds a preset threshold, a theoretical arrival time of the detonation wave at the detection position is obtained based on a preset combined gas pressure, and a reference average correction value is obtained based on the reference arrival time series and the theoretical arrival time of the detonation wave; Obtaining a target detonation velocity, and obtaining a target detonation wave arrival time at a detection position based on the target detonation velocity, and obtaining a target combined pressure based on the target detonation wave arrival time and a reference average correction amount, and conducting an experiment according to the target combined pressure in a current experimental cycle, and obtaining a current detonation wave arrival time at the detection position in the current experimental cycle; The current detonation velocity is obtained according to the current detonation wave arrival time, and the adjustment strategy is obtained according to the target detonation velocity and the current detonation velocity.

2. The detonation velocity adjustment and control method for a gas explosion composition according to claim 1, characterized in that: The adjustment strategy obtained according to the target detonation velocity and the current detonation velocity includes: comparing the magnitude relationship between the target detonation velocity and the current detonation velocity, and determining an adjustment direction of the combined gas pressure and an adjustment of a reference average correction amount based on the magnitude relationship; When the current detonation velocity is equal to the target detonation velocity, the adjustment strategy indicates that the target combined air pressure remains unchanged, and outputs a mapping relationship between a reference average correction amount and the target detonation velocity; When the current detonation velocity is less than the target detonation velocity, the adjustment strategy indicates increasing the target combined gas pressure and increasing the reference average correction amount; When the current detonation velocity is greater than the target detonation velocity, the adjustment strategy indicates reducing the target combined gas pressure and lowering the reference average correction amount.

3. The detonation velocity adjustment and control method for a gas explosion composition according to claim 1, characterized in that: The step of obtaining the theoretical detonation wave arrival time at the detection position according to the preset combined air pressure includes: Obtain theoretical detonation velocity based on preset combined air pressure; The theoretical detonation wave arrival time at the detection position is obtained according to the theoretical detonation velocity.

4. The detonation velocity adjustment and control method for a gas explosion composition according to claim 1, characterized in that: The obtaining of the target combined pressure according to the target detonation wave arrival time and the reference average correction amount includes: Obtaining the predicted detonation wave arrival time according to the target detonation wave arrival time and the reference average correction value; The predicted detonation velocity is obtained according to the predicted detonation wave arrival time, and the target combined air pressure is obtained according to the predicted detonation velocity.

5. The detonation velocity adjustment and control method for a gas explosion composition according to claim 1, characterized in that: The obtaining of the stability index according to the reference in-place time series includes: Obtain the average value of each reference detonation wave arrival time in the reference arrival time series; The difference between each reference detonation wave arrival time in the reference arrival time series and the average value is obtained, and the sum of all the differences is obtained, and the stability index is obtained according to the sum of all the differences.

6. The detonation velocity adjustment and control method for a gas explosion composition according to claim 1, characterized in that: The obtaining of a reference average correction value according to the reference arrival time sequence and the theoretical detonation wave arrival time comprises: The reference detonation wave arrival time in the reference arrival time sequence is obtained by subtracting the theoretical detonation wave arrival time to obtain a difference, and a reference average correction amount is obtained according to the multiple differences.

7. A detonation velocity regulation and control system applied to a gas explosion composition, characterized in that: The system is used to implement the detonation velocity regulation and control method for a gas explosion composition according to any one of claims 1 to 6, and the system comprises: An acquisition module is used to acquire multiple experimental cycles under a preset combined gas pressure, and obtain a reference detonation wave arrival time at the detection position in each experimental cycle, and arrange the multiple reference detonation wave arrival times in sequence according to the order of the experimental cycles to form a reference arrival time sequence; a data processing module for obtaining a stability index based on a reference arrival time series; if the stability index exceeds a preset threshold, obtaining a theoretical detonation wave arrival time at the detection position based on a preset combined gas pressure; and obtaining a reference average correction value based on the reference arrival time series and the theoretical detonation wave arrival time; a first regulating and controlling module, configured to obtain a target detonation velocity, obtain a target detonation wave arrival time at a detection position based on the target detonation velocity, obtain a target combined air pressure based on the target detonation wave arrival time and a reference average correction amount, conduct an experiment according to the target combined air pressure in a current experimental cycle, and obtain a current detonation wave arrival time at the detection position in the current experimental cycle; The second regulation control module is used to obtain the current detonation velocity according to the current detonation wave arrival time, and obtain the regulation strategy according to the target detonation velocity and the current detonation velocity.

8. The detonation velocity regulation and control system for a gas explosion composition according to claim 7, characterized in that: The second regulation control module is further configured to: comparing the magnitude relationship between the target detonation velocity and the current detonation velocity, and determining an adjustment direction of the combined gas pressure and an adjustment of a reference average correction amount based on the magnitude relationship; When the current detonation velocity is equal to the target detonation velocity, the adjustment strategy indicates that the target combined air pressure remains unchanged, and outputs a mapping relationship between a reference average correction amount and the target detonation velocity; When the current detonation velocity is less than the target detonation velocity, the adjustment strategy indicates increasing the target combined gas pressure and increasing the reference average correction amount; When the current detonation velocity is greater than the target detonation velocity, the adjustment strategy indicates reducing the target combined gas pressure and lowering the reference average correction amount.

9. The detonation velocity regulation and control system for gas explosion composition according to claim 7, characterized in that: The data processing module is further configured to: Obtain theoretical detonation velocity based on preset combined air pressure; The theoretical detonation wave arrival time at the detection position is obtained according to the theoretical detonation velocity.

10. The detonation velocity regulation and control system for gas explosion composition according to claim 7, characterized in that: The data processing module is further configured to: Obtaining the predicted detonation wave arrival time according to the target detonation wave arrival time and the reference average correction value; The predicted detonation velocity is obtained according to the predicted detonation wave arrival time, and the target combined air pressure is obtained according to the predicted detonation velocity.

Images