Method and system for improving prediction of the effect of a blast based on the length of the blast hole obstruction

Abstract

The application relates to the technical field of impact effect prediction improvement, and discloses an impact effect prediction improvement method and system based on a blast hole blockage length, which comprises the following steps: according to pre-acquired blast hole blockage detection data, analyzing the stress distribution in the blast hole to obtain a blast hole blockage analysis result; combining the blast hole blockage analysis result and preset impact parameters to construct an impact effect analysis model, predicting the dynamic distribution of energy in the blast hole to obtain an energy distribution result; based on the energy distribution result, taking each blast hole as a node and the position relationship between the blast holes as an edge to construct a blast area action graph; simulating the propagation process of impact stress waves and cracks, analyzing the node state after impact to obtain a first impact effect; analyzing the improvement direction through a preset impact improvement model to generate an improvement scheme; and the application can improve the measurement precision and accuracy of blast hole blockage parameters, enhance the authenticity and comprehensiveness of the overall blasting area fragmentation effect simulation process, and improve the impact effect.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of impact effect prediction improvement, and more particularly to an impact effect prediction improvement method and system based on borehole plugging length. BACKGROUND

[0002] At present, in the fields of mine exploitation, tunnel excavation, geotechnical engineering and the like, borehole plugging is one of the key links affecting the impact effect, economic benefits and safety in blasting impact operation. The plugging length and the state of the plugging material, such as the compactness and the particle size distribution, directly affect the propagation of the explosion stress wave, the expansion and working process of the explosion gas, and the expansion of the explosion-generated cracks, and further determine the rock crushing size, the blast pile shape, the fly rock distance and the impact vibration intensity.

[0003] The prior art has the following problems: relying on manual measurement or simple instrument detection, the precision is low and it is difficult to reflect the compactness state and internal structure of the plugging material, the judgment accuracy of the plugging quality is low; based on a single empirical formula to analyze the plugging section, ignoring the complex reflection and transmission process of the explosion stress wave in the plugging section, the prediction of the dynamic distribution of the explosion energy in the axial and radial directions of the borehole is not accurate; analyzing a simple single hole or single row of holes, ignoring the interaction between multiple boreholes in the impact area, it is difficult to truly reflect the impact process under complex geological conditions and charging structure; in order to solve at least one of the above problems, the present application proposes an impact effect prediction improvement method and system based on borehole plugging length. SUMMARY

[0004] In view of the deficiencies in the prior art, the purpose of the present application is to provide an impact effect prediction improvement method and system based on borehole plugging length, which can effectively solve the problems in the background art. The specific technical solutions of the present application are as follows:

[0005] The impact effect prediction improvement method based on borehole plugging length comprises:

[0006] According to the pre-acquired borehole plugging material detection data, the stress distribution inside the borehole is analyzed, the borehole plugging length and the plugging state are determined, and the borehole plugging analysis result is obtained;

[0007] Combined with the borehole plugging analysis result and the preset impact parameters, an impact effect analysis model is constructed, the stress wave propagation and gas expansion coupling in the plugging section are analyzed, the dynamic distribution of the energy in the borehole is predicted, and the energy distribution result is obtained;

[0008] Based on the energy distribution result, each borehole is taken as a node, and the positional relationship between the boreholes is taken as an edge, to construct a blast zone action graph;

[0009] According to the blast zone effect diagram, the propagation process of the shock stress wave and the crack is simulated through a preset shock effect prediction model, the node state after the shock is analyzed, and a first shock effect is obtained;

[0010] Based on the first shock effect, an improvement direction is analyzed through a preset shock improvement model, an improvement scheme is generated, and the shock effect is predicted and improved.

[0011] Specifically, the stress distribution inside the blast hole is analyzed according to the pre-acquired detection data of the plug in the blast hole, the plug length and the plug state of the blast hole are determined, and a plug analysis result of the blast hole is obtained, including:

[0012] According to the pre-acquired detection data of the plug in the blast hole, the stress variation of the stress signal in the direction perpendicular to the hole wall and the direction parallel to the hole wall is analyzed, a stress distribution map is constructed, the position distance of the plug and the blast hole mouth is identified, and the plug length of the blast hole is obtained.

[0013] The acoustic data characteristics corresponding to the plug length of the blast hole are analyzed, the plug profile along the plug length of the blast hole is constructed based on the acoustic data characteristics, and the plug state is analyzed and obtained.

[0014] In combination with the plug length and the plug state of the blast hole, the plug analysis result of the blast hole is obtained.

[0015] Specifically, the stress distribution map is constructed according to the pre-acquired detection data of the plug in the blast hole, the position distance of the plug and the blast hole mouth is identified, and the plug length of the blast hole is obtained, including:

[0016] According to the pre-acquired detection data of the plug in the blast hole, the stress variation of the stress signal in the direction perpendicular to the hole wall and the direction parallel to the hole wall is analyzed, the normal stress sequence and the tangential stress sequence are calculated, and the stress distribution map is constructed.

[0017] According to a preset distance interval, the ratio of the normal stress to the tangential stress at each position point in the stress distribution map is calculated, the variation of the ratio along the axial direction of the blast hole is analyzed, and a stress coupling ratio curve is obtained.

[0018] The inflection points in the stress coupling ratio curve are identified, and the blast hole cross sections corresponding to the inflection points are taken as a first interface set.

[0019] According to the pre-acquired radar signal in the blast hole, the mutation of the echo signal intensity is analyzed, the radar mutation position is selected from the first interface set, the corresponding first interface is taken as the position distance of the plug and the blast hole mouth, and the plug length of the blast hole is calculated.

[0020] Specifically, the analysis borehole plug length corresponds to the acoustic data characteristics, based on the acoustic data characteristics inversion along the borehole plug length distribution of the plug profile, analysis of the plug state, including:

[0021] According to the borehole plug length corresponds to the acoustic data characteristics, the energy attenuation ratio of each position point acoustic signal in the preset high frequency sub-band relative to the energy of the acoustic emission signal is calculated, and the energy attenuation feature sequence is obtained;

[0022] The shift amount of the center frequency of each position point acoustic signal along the axial direction is calculated, and the shift feature sequence is obtained;

[0023] The energy attenuation feature sequence and the shift feature sequence are respectively interpolated to map the density sequence and the granularity sequence;

[0024] The mutation point in the density sequence is identified, the jump of the granularity sequence at the mutation point position is analyzed, and the first plug region is divided, the average density in the region is calculated, the plug profile along the borehole plug length distribution is constructed, and the plug state is analyzed.

[0025] Specifically, the analysis results of the borehole plug and the preset impact parameters are combined to construct an impact effect analysis model, the stress wave propagation and gas expansion coupling in the plug section are analyzed, the dynamic distribution of energy in the borehole is predicted, and the energy distribution result is obtained, including:

[0026] According to the density sequence in the analysis results of the borehole plug, the impedance change and the gas permeation are analyzed, and the impact effect analysis model is constructed combined with the preset impact parameters;

[0027] The stress wave propagation and gas expansion coupling in the plug section are analyzed through the impact effect analysis model, the dynamic distribution of energy in the borehole is predicted, and the energy distribution result is obtained.

[0028] Specifically, the stress wave propagation and gas expansion coupling in the plug section are analyzed through the impact effect analysis model, the dynamic distribution of energy in the borehole is predicted, and the energy distribution result is obtained, including:

[0029] The stress wave propagation process is simulated through the impact effect analysis model, the energy accumulation condition is obtained, the expansion of the gas and the movement process of the plug are simulated based on the energy accumulation condition, and the gas expansion efficiency is calculated;

[0030] The energy accumulation condition and the gas expansion efficiency are coupled to calculate the energy accumulation distribution value;

[0031] Based on the energy accumulation distribution value, the dynamic distribution of energy in the borehole is predicted, and the energy distribution result is obtained.

[0032] Specifically, the first impact effect is obtained by simulating the propagation process of the impact stress wave and the crack through a preset impact effect prediction model according to the blast zone effect diagram, and analyzing the state of the node after the impact, and the first impact effect comprises:

[0033] According to the blast zone effect diagram, the breaking performance and impact time of each node are analyzed, and the dynamic stress intensity is calculated;

[0034] Based on the dynamic stress intensity, the propagation process of the impact stress wave and the crack is simulated through a preset impact effect prediction model, the state of the node after the impact is analyzed, and the first impact effect is obtained.

[0035] Specifically, the first impact effect is obtained by simulating the propagation process of the impact stress wave and the crack through a preset impact effect prediction model according to the blast zone effect diagram, and analyzing the state of the node after the impact, and the first impact effect comprises:

[0036] Based on the dynamic stress intensity, the propagation process of the impact stress wave and the crack is simulated through a preset impact effect prediction model, the state of the node after the impact is analyzed, and the first impact effect is obtained.

[0037] For each node, the impact efficiency is calculated by combining the dynamic stress intensity of the node and the corresponding instantaneous stress, the state of the node after the impact is analyzed, and the impact damage state is obtained.

[0038] According to the impact damage state, the impact effect index of the corresponding node position is mapped, and the first impact effect is analyzed.

[0039] Specifically, based on the first impact effect, the improvement direction is analyzed through a preset impact improvement model, and an improvement scheme is generated to predict and improve the impact effect, comprising:

[0040] Based on the first impact effect, the impact effect and the improvement direction are analyzed through a preset impact improvement model, and the effect defect node and the effect excess node are identified.

[0041] According to the effect defect node and the effect excess node, the associated adjustment node set is screened out;

[0042] For each node in the adjustment node set, the improvement parameter is analyzed, and the corresponding improvement scheme is generated to predict and improve the impact effect.

[0043] The impact effect prediction and improvement system based on the borehole plugging length is used to realize the impact effect prediction and improvement method based on the borehole plugging length, comprising:

[0044] The borehole plugging analysis module analyzes the stress distribution inside the borehole according to the pre-acquired borehole plugging detection data, determines the borehole plugging length and the plugging state, and obtains the borehole plugging analysis result.

[0045] An energy distribution analysis module combines the blast hole plugging analysis result and the preset impact parameter to construct an impact effect analysis model, predicts the dynamic distribution of energy in the blast hole by analyzing the stress wave propagation and gas expansion coupling in the plugging section, and obtains an energy distribution result;

[0046] A blast zone action graph construction module constructs a blast zone action graph based on the energy distribution result, taking each blast hole as a node and the positional relationship between blast holes as an edge;

[0047] An impact effect prediction module simulates the propagation process of impact stress waves and cracks according to the blast zone action graph through a preset impact effect prediction model, analyzes the node state after impact, and obtains a first impact effect;

[0048] An impact improvement module analyzes the improvement direction through a preset impact improvement model based on the first impact effect, generates an improvement scheme, and predicts the improvement of the impact effect.

[0049] The beneficial effects of the present application are as follows: through multi-source data fusion, the blast hole plugging length and internal density and particle size profile can be accurately identified; the dynamic distribution of explosion energy in the axial direction of the blast hole is quantified by combining the coupling effect of stress wave propagation and gas expansion; the synergistic effect of stress wave and crack propagation between multiple holes is simulated by constructing a blast zone action graph, and the impact effect is dynamically predicted; the defect area is identified based on the impact effect prediction result, and a targeted improvement method is developed; the measurement precision and accuracy of the blast hole plugging parameters can be improved, the accuracy of the impact energy dynamic distribution prediction result can be improved, the authenticity and comprehensiveness of the overall blast zone fragmentation effect simulation process can be enhanced, the impact scheme can be automatically optimized, and the impact effect can be improved. BRIEF DESCRIPTION OF DRAWINGS

[0050] Figure 1 A work flow chart of the impact effect prediction and improvement method based on the blast hole plugging length in the embodiments of the present application;

[0051] Figure 2 A schematic diagram of the blast zone action graph in the embodiments of the present application;

[0052] Figure 3 A work flow chart of the plugging state analysis process in the embodiments of the present application;

[0053] Figure 4 A structural schematic diagram of the impact effect prediction and improvement system based on the blast hole plugging length in the embodiments of the present application. DETAILED DESCRIPTION

[0054] The present application will be further described in detail below with reference to the drawings and embodiments.

[0055] In the embodiments of the present application, the word "exemplary" or "for example" is used to mean serving as an example, instance, or illustration. Any implementation or design solution described as "exemplary" or "for example" in the embodiments of the present application should not be construed as being preferred or advantageous over other implementations or design solutions. Rather, the use of "exemplary" or "for example" is intended to present concepts in a concrete manner.

[0056] Hereinafter, the terms "first", "second", and the like are used generically and are only intended for the purpose of description, and should not be construed as indicating or implying relative importance or implicitly indicating the number of the technical features indicated. Therefore, the features defined with "first", "second" can explicitly or implicitly include one or more of the features. In the description of the embodiments of the present application, the meaning of "a plurality of" is two or more, unless otherwise specified.

[0057] Reference Figure 1 As shown in the specific embodiments of the present application, the impact effect prediction improvement method based on the borehole plugging length includes:

[0058] S101, according to the pre-acquired borehole plugging detection data, analyze the stress distribution inside the borehole, determine the borehole plugging length and the plugging state, and obtain the borehole plugging analysis result;

[0059] S102, combine the borehole plugging analysis result and the preset impact parameter, construct an impact effect analysis model, analyze the stress wave propagation and gas expansion coupling in the plugging section, predict the dynamic distribution of energy in the borehole, and obtain the energy distribution result;

[0060] S103, based on the energy distribution result, taking each borehole as a node and the position relationship between the boreholes as an edge, construct a blast area action graph;

[0061] S104, according to the blast area action graph, simulate the propagation process of the impact stress wave and the crack through the preset impact effect prediction model, analyze the node state after the impact, and obtain the first impact effect;

[0062] S105, based on the first impact effect, analyze the improvement direction through the preset impact improvement model, generate an improvement scheme, and predict the improvement of the impact effect.

[0063] In the embodiment, the original detection data distributed along the axis of the blast hole is obtained by a sensor pre-installed in the blast hole, the sensor includes but is not limited to a miniature stress meter, an acoustic wave emission / receiving probe, and a radar antenna array. According to the collected plug detection data in the blast hole, the stress distribution inside the blast hole is analyzed by fusing multi-source data, the length and state of the blast hole plug are determined, and the blast hole plug analysis result is obtained. Through multi-source data fusion analysis, accurate initial boundary conditions and physical parameters are provided for subsequent impact energy analysis, the input error of the impact effect analysis model caused by the uncertainty of the plug parameters is avoided, and the calculation error of the explosion gas action time caused by the length estimation deviation can be avoided by accurately determining the length of the plug. By identifying the accurate plug state, accurate state support is provided for analyzing the stress wave propagation impedance matching, gas permeation leakage channel and other processes, and the authenticity and reliability of the impact effect analysis process are improved.

[0064] Specifically, an impact effect analysis model is constructed in combination with the blast hole plug analysis result and the preset impact parameters, the propagation process of the initial stress wave generated after the explosion of the explosive in the blast hole is simulated, the energy distribution result is obtained by analyzing the stress wave propagation and gas expansion coupling in the plug section, and predicting the dynamic distribution of energy in the blast hole. By constructing the impact effect analysis model, the energy dynamic distribution in the blast hole plug section can be accurately analyzed, the part of the explosion energy that is effectively used to break the rock and the part of the explosion energy that is consumed to push or compress the plug without effect are identified, the problems of early gas leakage caused by the too loose top of the plug section or the excessive attenuation of the stress wave caused by the too dense middle section are identified, the energy distribution is accurately analyzed, accurate data support is provided for predicting the impact effect, and the accuracy of the impact effect analysis result is improved.

[0065] As shown in Figure 2 Based on the energy distribution result, each blast hole is taken as a node, and the positional relationship between the blast holes is taken as an edge to construct a blast zone action graph. The energy distribution result of each blast hole calculated is taken as the core attribute of the corresponding blast hole node, including but not limited to the distribution information of energy along its own axial direction. According to the three-dimensional spatial positional relationship between the blast hole nodes, the attenuation degree of the stress wave when propagating from one node to another node is calculated as the weight of the edge based on the positional relationship, and the connection edge is constructed between the corresponding nodes to obtain the blast zone action graph taking the blast holes as nodes and the interaction relationship as edges. By constructing the blast zone action graph, the inter-hole synergistic effect and interference effect can be quantitatively analyzed. The inter-hole synergistic effect includes but is not limited to stress wave superposition and enhanced fragmentation, and the interference effect includes but is not limited to the effect of the post-blast hole affected by the cracks or stress field formed by the pre-blast hole. By simulating the overall impact process through analyzing the inter-hole synergistic effect, the accuracy of the prediction results of the overall fragmentation size distribution, the blast pile shape, and the contour effect of the blast zone can be improved.

[0066] Specifically, according to the blast zone effect diagram, the propagation process of the shock stress wave and the crack is simulated through a preset shock effect prediction model, the propagation of the stress wave field generated by the node to the adjacent node area is analyzed according to the initiation time sequence of the node, the state of the node after the shock is analyzed, and the first shock effect is obtained; through the dynamic numerical simulation based on the blast zone effect diagram structure, the spatial distribution difference of the shock effect in the blast zone can be predicted, the area where the stress wave superposition is insufficient due to improper inter-hole delay can be accurately identified, and the accuracy of the shock effect prediction result is improved.

[0067] Specifically, based on the first shock effect, the improvement direction is analyzed through a preset shock improvement model, the area where the effect does not meet the standard and the area where energy is wasted or over damaged are automatically identified, and an improvement scheme is generated to predict and improve the shock effect. By automatically identifying the improvement area and generating the corresponding improvement scheme, the improvement method can be quickly obtained, the optimization efficiency is improved, the blindness and partiality of experience parameter adjustment are avoided, the obtained improvement scheme has clear pertinence and operability, and the overall shock effect, safety and economy of the shock process are improved.

[0068] The present application can accurately identify the borehole plugging length and internal density, and the particle size profile through multi-source data fusion; quantize the dynamic distribution of explosion energy in the borehole axial direction by combining the coupled stress wave propagation and gas expansion effect; simulate the synergistic effect of stress wave and crack propagation among multiple holes by constructing a blast zone effect diagram to dynamically predict the shock effect; identify the defect area based on the shock effect prediction result and develop a targeted improvement method; the measurement accuracy and accuracy of the borehole plugging parameter can be improved, the accuracy of the shock energy dynamic distribution prediction result can be improved, the authenticity and comprehensiveness of the overall blasting effect simulation process can be enhanced, the shock scheme can be automatically optimized, and the shock effect can be improved.

[0069] Further, according to the pre-acquired borehole plugging object detection data, the stress distribution in the borehole is analyzed, the borehole plugging length and the plugging state are determined, and the borehole plugging analysis result is obtained, including:

[0070] S201, according to the pre-acquired borehole plugging object detection data, the stress signal is analyzed in the stress change in the direction perpendicular to the hole wall and the direction parallel to the hole wall, a stress distribution diagram is constructed, the position distance of the plugging object and the borehole mouth is identified, and the borehole plugging length is obtained;

[0071] S202, analyze the acoustic data characteristics corresponding to the borehole plugging length, and construct the plugging profile along the borehole plugging length based on the acoustic data characteristics, and analyze the plugging state;

[0072] S203, combined with the borehole plugging length and the plugging state, the borehole plugging analysis result is obtained.

[0073] In the embodiment, according to the pre-acquired borehole plug detection data, including but not limited to the signals recorded by the stress sensor array arranged along the borehole axial direction, the stress signal changes in the direction perpendicular to the hole wall and the direction parallel to the hole wall are analyzed respectively, the synthetic stress vector detected by the stress sensor at each measuring point is decomposed into the normal stress component perpendicular to the hole wall and the tangential stress component parallel to the hole wall according to the axial and radial geometric coordinate system of the borehole, the stress distribution map is constructed, the ratio of the normal stress to the tangential stress at each measuring point is calculated, the stress coupling ratio curve is constructed, the position distance of the plug to the borehole mouth is identified, and the borehole plug length is obtained.

[0074] It should be noted that the borehole plug length can be accurately calculated through data analysis and data fusion, avoiding the misjudgment problem caused by single stress mutation or radar reflection in complex hole environment, and the position distance of the plug to the borehole mouth is identified through analysis of the stress coupling ratio curve combined with the radar signal, and the accuracy and reliability of the calculation result of the borehole plug length can be improved through double analysis and judgment of mechanics and electromagnetism.

[0075] Specifically, for the determined plug length section, the corresponding acoustic wave detection data is obtained by the acoustic wave emission source and the receiver array arranged along the hole, the acoustic wave data characteristics corresponding to the borehole plug length are analyzed, the plug profile along the borehole plug length distribution is constructed based on the acoustic wave data characteristics inversion, and the plug state is obtained by analysis; the borehole plug state can be converted into a specific quantitative relationship through acoustic wave characteristic inversion analysis of the borehole plug state, and whether the internal structure of the plug is uniform, whether there is a weak zone or gap can be accurately analyzed through the plug profile, which provides accurate data support for predicting whether the explosion gas can be effectively sealed and whether the stress wave propagation impedance is matched, and improves the accuracy of the analysis result of the borehole plug.

[0076] Specifically, the borehole plug length and the plug state are associated, the plug length is taken as the core metadata, and the plug profile is taken as the detailed attribute under the length. Through data association, the spatial coordinate origin of the plug profile can be aligned with the starting point defined by the length, ensuring the consistency and integrity of the information, and obtaining the borehole plug analysis result. Through data association and integration, accurate data support is provided for the impact effect prediction and improvement process, the consistency and accuracy of the data in the transmission process are improved, the calculation deviation caused by non-uniform or incomplete input parameters is eliminated, and the logic consistency of the impact effect analysis and prediction process and the accuracy and reliability of the prediction result are improved.

[0077] Further, according to the pre-acquired borehole plug detection data, the stress signal changes in the direction perpendicular to the hole wall and the direction parallel to the hole wall are analyzed respectively, the stress distribution map is constructed, the position distance of the plug to the borehole mouth is identified, and the borehole plug length is obtained, including:

[0078] S301, according to the pre-acquired borehole plug detection data, analyze the stress signal respectively in the direction perpendicular to the hole wall and the direction parallel to the hole wall, calculate the normal stress sequence and the tangential stress sequence, and construct the stress distribution map;

[0079] S302, according to the pre-set distance interval, calculate the ratio of the normal stress and the tangential stress at each position point in the stress distribution map, analyze the change of the ratio along the borehole axial direction, and obtain the stress coupling ratio curve;

[0080] S303, identify the inflection point in the stress coupling ratio curve, and take the borehole cross section corresponding to the inflection point as the first interface set;

[0081] S304, according to the pre-acquired radar signal in the borehole, analyze the mutation of the echo signal intensity, select the radar mutation position from the first interface set, take the corresponding first interface as the position distance between the plug and the borehole mouth, and calculate the borehole plug length.

[0082] In this embodiment, according to the pre-acquired borehole plug detection data, including but not limited to the synthetic stress vector collected by the three-axis stress sensor arranged at a certain interval along the borehole axial direction, the stress signal is analyzed respectively in the direction perpendicular to the hole wall and the direction parallel to the hole wall, the normal stress perpendicular to the hole wall reflects the direct compression or support of the medium to the hole wall, and the tangential stress parallel to the hole wall reflects the friction, shear or viscous effect between the medium and the hole wall; according to the geometric axial direction of the borehole, a space coordinate system is established, in which the coordinate axis Z is parallel to the borehole axis, and the X and Y axes are located in the tangent plane of the hole wall; the original stress vector measured by the sensor is decomposed into the space coordinate system through coordinate transformation; the normal stress component perpendicular to the hole wall is obtained by calculating the projection of the stress vector in the unit normal vector direction of the hole wall, and the tangential component in the X and Y axial directions is obtained by calculating the synthetic stress vector in the tangent plane, and the normal stress value and the tangential stress value are obtained; sorting and connecting according to the axial position, forming the normal stress sequence and the tangential stress sequence; taking the axial position as the horizontal coordinate and the normal stress and the tangential stress as the vertical coordinate, a stress distribution map reflecting the change trend of the stress component is constructed.

[0083] It should be noted that by decomposing the detected combined stress into normal component and tangential component and constructing the distribution sequence map, the independent physical quantity for judging the medium constraint state can be separated out, the normal stress sequence reflects the pressure distribution of the medium to the hole wall, and the tangential stress sequence reflects the shear constraint condition between the medium and the hole wall, the specific relationship between the two components can be analyzed by independent analysis, the mechanical mutation condition between the air section and the plug section is identified, and the medium interface is identified, thereby improving the accuracy of the identification result of the borehole plug length.

[0084] Specifically, in the air section in the blast hole, the medium has little effect on the hole wall, the normal stress and the tangential stress are close to zero, and the corresponding ratio presents irregular oscillation or cannot be defined. When entering the plug section, the solid particles contact the hole wall, and under the action of the vertical pressure corresponding to the normal stress, the lateral friction resistance corresponding to the tangential stress is generated, and the ratio between them presents a relatively stable proportion or coupling relationship. According to the preset distance interval, the distance interval can be set according to the sensor arrangement interval, the normal stress value and the tangential stress value of each axial position point are extracted, the ratio of the absolute value of the normal stress to the absolute value of the tangential stress is calculated for each position point, the stress coupling ratio is obtained, and the stress coupling ratios calculated by all position points are arranged according to the corresponding axial positions to obtain a stress coupling ratio curve.

[0085] It needs to be emphasized that by constructing and analyzing the stress coupling ratio curve, the independent stress component information is fused, the interference condition of the absolute value of a single stress component can be effectively inhibited, the characteristics reflecting the essential difference of the medium constraint state are amplified, the complex two-dimensional stress distribution information is compressed into a trend line, the curve shape mutation point can be quickly screened out by analyzing the trend line, the complexity of the analysis process is reduced, and the analysis efficiency and accuracy of the blast hole plugging length are improved.

[0086] Specifically, when the detection point enters the plug from the air, the stress coupling relationship will change discontinuously, and at the corresponding inflection point position in the stress coupling ratio curve, the first difference sequence of the stress coupling ratio curve is calculated, the points where the difference value changes in sign or the absolute value exceeds the preset absolute value threshold are screened out, the preset absolute value threshold can be set according to the accuracy requirement of the blast hole plugging length calculation, each identified inflection point corresponds to a cross section of the blast hole, and the cross sections are integrated to form a first interface set. Through the inflection point screening, the mutation point can be quickly analyzed, which provides a reference for the plug interface judgment process, improves the accuracy and reliability of the interface detection, and improves the efficiency of the analysis process.

[0087] Specifically, according to the radar signal pre-acquired in the blast hole, the radar data along the axial direction of the blast hole collected at the same position as the stress detection is called, the change of the echo signal strength along the axial direction in the radar data is analyzed, the position where the echo strength suddenly and significantly increases is screened out as the radar mutation position, the position corresponds to the physical interface where the dielectric constant changes, the first interface set is compared with the radar mutation position, the coordinate points appearing in both the first interface set and the radar mutation position list are identified, the position points determined through the double verification of the stress mechanical response and the radar electromagnetic response can reflect the real interface position between the air and the plug, the interface closest to the blast hole mouth and passing the double verification is taken as the starting interface of the plug, the axial coordinate value of the interface reflects the distance from the upper surface of the plug to the blast hole mouth, and the blast hole plugging length is obtained.

[0088] It should be noted that the determination of the borehole blockage length by stress and radar information double detection can reduce the misjudgment risk caused by single detection method. The multiple candidate interfaces provided by stress analysis are verified and screened by radar detection, which can effectively filter out the misjudgment results caused by the limitation of single method, improve the accuracy and reliability of the final blockage length measurement result, and provide accurate initial length parameters for the impact effect analysis process.

[0089] As shown in Figure 3 , the acoustic wave data characteristics corresponding to the borehole blockage length are analyzed, and the blockage profile along the borehole blockage length is constructed based on the acoustic wave data characteristics, and the blockage state is analyzed, including:

[0090] S401, according to the acoustic wave data characteristics corresponding to the borehole blockage length, the energy attenuation ratio of the acoustic wave signal at each position point in the preset high frequency sub-band relative to the energy of the acoustic wave emission signal is calculated, and the energy attenuation feature sequence is obtained;

[0091] S402, the shift amount of the center frequency of the acoustic wave signal at each position point along the axial direction is calculated, and the shift feature sequence is obtained;

[0092] S403, the energy attenuation feature sequence and the shift feature sequence are respectively interpolated to map the density sequence and the granularity sequence;

[0093] S404, the mutation point in the density sequence is identified, the jump condition of the granularity sequence at the mutation point position is analyzed, and the first blockage region is divided, the average density in the region is calculated, the blockage profile along the borehole blockage length is constructed, and the blockage state is analyzed.

[0094] In this embodiment, for the borehole blockage length section, the acoustic wave signal excited by the transmitting transducer and collected by the receiving transducer array arranged along the axial direction is obtained, for each receiving point position, the time-frequency analysis is performed on the received time-domain acoustic wave signal, and the different frequency components are separated; according to the detection system characteristics and the typical blockage response, the high frequency sub-band is set, the total energy of the received signal in the high frequency sub-band is calculated, the original waveform of the transmission source signal is obtained, and the reference energy of the original waveform in the same high frequency sub-band is calculated; for each position point, the ratio of the high frequency sub-band energy of the received signal to the high frequency sub-band reference energy of the transmission signal is calculated to obtain the energy attenuation ratio, and the smaller the value, the more severe the energy attenuation; along the borehole blockage length direction, the energy attenuation ratios of all sampling position points are arranged to obtain the energy attenuation feature sequence.

[0095] It should be noted that by calculating the energy attenuation ratio in a specific high-frequency sub-band, the medium density can be reflected, the full-band acoustic signal can be converted into a single characteristic quantity reflecting the medium structure state, the response of the signal to the mechanical properties of the medium can be enhanced, the influence of low-frequency interference can be suppressed, accurate data support can be provided for predicting the state of the plug, and the accuracy of the analysis result of the plug state can be improved.

[0096] Specifically, for the acoustic signal of each receiving point position, the power spectral density is obtained by Fourier transform analysis, and the center frequency is determined by identifying the frequency corresponding to the peak point. For each position point, the difference between the center frequency of the received signal and the center frequency of the transmitted signal is calculated to obtain the center frequency offset of the point. The difference corresponding to each sampling point is calculated along the borehole axis and arranged in order of position to obtain the offset feature sequence. By calculating the offset of the center frequency of the acoustic wave, the distribution of the medium can be reflected. The offset feature sequence can distinguish between attenuation caused by changes in density and scattering effects caused by changes in particle composition, reflect the internal structural heterogeneity of the plug, and enhance the accuracy of the state prediction result.

[0097] Specifically, the energy attenuation feature sequence and the offset feature sequence are respectively interpolated to generate two continuous curves, which represent the continuous functions of the energy attenuation ratio and the center frequency offset with depth, respectively. By using a neural network model pre-trained using a large amount of historical energy attenuation ratio data, the interpolated energy attenuation ratio of each position is mapped to the corresponding density value. By using a random forest model pre-trained using a large amount of historical center frequency offset data, the interpolated center frequency offset of each position is mapped to the corresponding particle size value. After mapping all axial positions, the density sequence and the particle size sequence are combined to reflect the compaction degree and the change in particle composition of the plug along the axial direction.

[0098] It should be noted that by data interpolation and mapping, parameter sequences for mechanical analysis and judgment can be constructed. The density sequence reflects where the plug is well compacted and where it is loose. The particle size sequence reflects where the plug has coarse particles and where it has fine particles. Interpolation processing can ensure the continuity of the profile in space, eliminate information gaps between discrete sampling points, and improve the accuracy of the identification result of the plug state.

[0099] Specifically, by analyzing the extreme points of the first derivative, the mutation points in the compactness sequence are identified. At each mutation point position, it is analyzed whether the granularity sequence also has a significant change at the position. When the compactness and granularity both show mutation characteristics near the same position, the position is taken as the layering interface. Based on the confirmed interface, the entire plug length is divided into several continuous first plug regions. For each divided region, the average value of the compactness of all points inside the region is calculated as the average compactness of the region. Taking the borehole axial depth as the horizontal axis and the compactness and granularity as the vertical axis attributes, the divided region boundaries and the average compactness values of each region are labeled, and a plug profile is constructed. According to the profile, the plug state is analyzed. For example, the plug is divided into three layers, the top 0-0.5 meters is a loose coarse particle zone, the middle 0.6-1.2 meters is a dense fine particle zone, and the bottom 1.3 meters to the hole bottom is a medium dense mixed zone.

[0100] It should be noted that by combining multi-parameter joint mutation point recognition and region division, the internal structure of the plug can be layered, and a corresponding comprehensive plug state report can be generated. The state report can determine whether there is layering, how many layers, and the compactness and particle characteristics of each layer. The impact energy action can be accurately predicted, and the accuracy of the impact effect prediction result can be improved.

[0101] Further, combined with the blast hole plug analysis result and the preset impact parameter, an impact effect analysis model is constructed. By analyzing the stress wave propagation and gas expansion coupling in the plug section, the dynamic distribution of energy in the blast hole is predicted, and the energy distribution result is obtained, including:

[0102] S501, according to the compactness sequence in the blast hole plug analysis result, the impedance change and gas permeation are analyzed, and combined with the preset impact parameter, an impact effect analysis model is constructed;

[0103] S502, by analyzing the stress wave propagation and gas expansion coupling in the plug section through the impact effect analysis model, the dynamic distribution of energy in the blast hole is predicted, and the energy distribution result is obtained.

[0104] In this embodiment, according to the compactness sequence in the blast hole plugging analysis result, by using the nonlinear correlation model pre-trained by a large number of historical compactness sequences, the impedance change and gas permeation situation are mapped, the compactness value of each axial position point is mapped to the corresponding dynamic wave impedance value and equivalent gas permeation coefficient, and the preset impact parameters are combined, including but not limited to explosive type and performance parameters, charge amount, charge structure, initiation point position, etc., to construct an impact effect analysis model. The model includes a wave equation describing the propagation of stress waves in a variable impedance medium, and a fluid dynamics equation describing the expansion, seepage and movement of the plugging object of the explosion gas in a variable permeability porous medium. The existing technology includes a large number of model construction processes, which will not be repeated here.

[0105] It should be noted that by mapping the compactness to the impedance and permeation parameters and constructing the impact effect analysis model, the energy folding back situation caused by the strong reflection of the initial stress wave due to the loose top of the plugging section, the easy and rapid leakage of the explosion gas can be quantitatively analyzed; and the basis for accurately predicting the dynamic distribution of energy at different axial positions in the blast hole is provided.

[0106] Specifically, the coupling of stress wave propagation and gas expansion in the plugging section is analyzed by the impact effect analysis model, the dynamic distribution of energy in the blast hole is predicted, and the energy distribution result is obtained. By simulating and analyzing the coupling of stress wave propagation and gas expansion in the plugging section, the part of the total explosion energy that is converted into effective mechanical energy for breaking rocks and moving rock blocks, and the part that is consumed in invalid work can be identified, thereby improving the accuracy of the energy distribution result.

[0107] Further, the coupling of stress wave propagation and gas expansion in the plugging section is analyzed by the impact effect analysis model, the dynamic distribution of energy in the blast hole is predicted, and the energy distribution result is obtained, including:

[0108] S601, simulating the stress wave propagation process by the impact effect analysis model to obtain the energy accumulation situation, based on the energy accumulation situation, simulating the expansion of the gas and the movement process of the plugging object, and calculating the gas expansion efficiency;

[0109] S602, coupling the energy accumulation situation and the gas expansion efficiency, and calculating the energy accumulation distribution value;

[0110] S603, based on the energy accumulation distribution value, predicting the dynamic distribution of energy in the blast hole, and obtaining the energy distribution result.

[0111] In the embodiment, the stress wave propagation process is simulated by the shock effect analysis model, the initial conditions are determined according to the charge parameters and the impedance distribution of the blockage along the axial direction, the stress wave propagation process is simulated based on the initial conditions, for the blockage section, the change of the density along the axial direction is mapped as the spatial distribution of the wave impedance, the model calculates the transmission and reflection coefficients of the stress wave at each impedance change interface, the input energy and the output energy at each position are integrated to obtain the corresponding energy accumulation, and the gas expansion process is further simulated, the calculated energy accumulation is taken as the initial state input, the work done by the gas pressure to drive the blockage movement and the work done by the gas pressure directly acting on the hole wall rock are calculated, and the corresponding gas expansion efficiency is obtained by summation.

[0112] It should be noted that by simulating the energy accumulation and the gas expansion process, the enhancement effect of the stress wave on the subsequent gas expansion efficiency can be reflected, the accuracy of the calculated gas expansion efficiency can be improved, and through correlation analysis, the accuracy and authenticity of the calculation result of the gas expansion efficiency can be improved.

[0113] Specifically, the energy accumulation distribution value is calculated by coupling the energy accumulation and the gas expansion efficiency by using the deep neural network model pre-trained by a large amount of historical energy accumulation and gas expansion efficiency data; by coupling the energy accumulation and the gas expansion efficiency, the synergistic or competitive effect of the stress wave and the gas expansion energy in breaking the rock can be reflected, the coupling calculation can avoid overestimating the contribution of the gas energy in the low-efficiency area or underestimating its contribution in the high-efficiency area, and the calculated energy accumulation distribution value can accurately reflect the actual impact physical process, thereby providing accurate data support for the prediction and optimization process of the impact effect.

[0114] Specifically, based on the energy accumulation distribution value, the relative energy distribution proportion corresponding to each axial position is obtained by calculating the ratio between the energy accumulation distribution value of each position and the total sum of the energy accumulation distribution values; in time sequence, the energy accumulation change curve at different depths over time is generated to obtain the energy distribution result, reflecting the dynamic process of energy action. By analyzing the energy distribution result, the dynamic characteristics and structural defects of energy distribution can be analyzed, thereby providing accurate data support for the impact effect prediction and optimization adjustment, and improving the accuracy of the impact effect prediction and analysis result.

[0115] Further, according to the blast area effect diagram, the propagation process of the shock stress wave and the crack is simulated by a preset impact effect prediction model, the state of each node after the impact is analyzed, and a first impact effect is obtained, including:

[0116] S701, according to the blast area effect diagram, the breaking performance and the impact time of each node are analyzed, and the dynamic stress intensity is calculated;

[0117] S702, simulate the propagation process of the impact stress wave and the crack based on the dynamic stress intensity through a preset impact effect prediction model, analyze the state of the node after the impact, and obtain a first impact effect.

[0118] In this embodiment, according to the blast area action diagram, the attribute information corresponding to each node is analyzed, based on the energy distribution characteristics of the node, the energy and the rock mass are analyzed to map the corresponding static fragmentation potential index to obtain the fragmentation performance; the impact time is obtained from the preset initiation time sequence table; in the blast area action diagram, all initiated nodes are traversed, the spatial distance, relative orientation, initiation time difference and energy characteristics of the nodes are calculated, the instantaneous stress state of the stress wave generated by them when propagating to the node position is calculated, and vector superposition is performed, the superposed stress field is coupled and superimposed with the reference stress field generated by the node itself explosion to obtain the dynamic stress intensity, reflecting the real mechanical load parameter of the node under the specific time sequence environment. The dynamic stress intensity of each node is calculated to fuse the time sequence coupling effect, which provides an accurate data basis for subsequent simulation of stress wave superposition, crack penetration and other dynamic processes.

[0119] Specifically, based on the dynamic stress intensity, the propagation process of the impact stress wave and the crack is simulated through a preset impact effect prediction model, the state of the node after the impact is analyzed, and a first impact effect is obtained. Through simulation and analysis of the impact effect, the impact blind area and other positions caused by improper inter-hole delay can be identified in the case of multi-hole delay impact, and the first impact effect obtained reflects the overall fragmentation effect and the specific damage distribution, which provides an accurate data basis for formulating improvement schemes.

[0120] Further, based on the dynamic stress intensity, the propagation process of the impact stress wave and the crack is simulated through a preset impact effect prediction model, the state of the node after the impact is analyzed, and a first impact effect is obtained, including:

[0121] S801, based on the dynamic stress intensity, simulate the propagation process of the impact stress wave and the crack through a preset impact effect prediction model according to the impact time sequence, calculate the instantaneous stress caused by the stress wave generated by the node to the adjacent node;

[0122] S802, for each node, calculate the impact efficiency combining the dynamic stress intensity of the node and the corresponding instantaneous stress, analyze the state of the node after the impact, and obtain the impact damage state;

[0123] S803, according to the impact damage state, map the impact effect index of the corresponding node position, and analyze to obtain a first impact effect.

[0124] In this embodiment, based on the dynamic stress intensity, the propagation process of the impact stress wave and the crack is simulated in the order of the impact time by a preset impact effect prediction model, which includes but is not limited to a deep neural network model pre-trained by a large amount of historical dynamic stress intensity data. The model sets the time step according to the initiation timing and the impact effect analysis accuracy requirement. In each time step, for the node currently being initiated, the model generates an outward propagating stress wave according to its dynamic stress intensity, calculates the propagation process of the stress wave in three-dimensional space, including but not limited to geometric attenuation with distance, internal friction attenuation due to material inelasticity, and reflection and refraction when encountering geological interfaces. For each target node that has not been initiated or has not completed energy release, the model monitors and calculates the stress tensor component instantaneous value generated at the corresponding target node coordinate position by the stress wave propagated from all initiated nodes, and performs vector superposition. The superimposed stress tensor reflects the instantaneous stress that the target node position bears due to initiation of other blastholes at the simulation time.

[0125] It should be noted that by dynamically calculating the mutual interference process of stress waves between nodes in time sequence, the inter-hole dynamic interference can be quantified in impact effect prediction, the instantaneous stress can be accurately analyzed, accurate data can be provided for accurately calculating the real impact efficiency of each node in the dynamic environment, and the prediction deviation caused by independent impact analysis of multi-hole impact can be avoided.

[0126] Specifically, for each node, the dynamic stress intensity and the corresponding instantaneous stress are superimposed and mapped to calculate the impact efficiency by using a neural network model pre-trained by a large amount of historical dynamic stress intensity and instantaneous stress data, the impact efficiency is superimposed on the initial state of the node to analyze the state of the node after impact, and the impact damage state is obtained. By analyzing the impact damage state of the node, the real damage efficiency of each blasthole can be quantified, the accuracy and authenticity of the impact damage state can be improved, accurate data support can be provided for analyzing the spatial distribution of the crushing effect in the blast area, and accurate direction can be provided for the effect improvement process.

[0127] Specifically, according to the impact damage state, it is mapped to one or more impact effect indicators by a preset mapping rule library, the mapping rule library is established based on a large amount of experimental data, field calibration or theoretical analysis, and the mapped indicators include but are not limited to qualitative grades and quantitative descriptors. The mapping results of all blasthole nodes are integrated and analyzed to obtain the first impact effect. By analyzing the first impact effect, the impact effect situation can be predicted, potential quality defects and over-crushing areas can be warned, accurate optimization direction can be provided for the impact effect improvement process, and the controllability and safety of the impact engineering are improved.

[0128] Further, based on the first impact effect, an improvement direction is analyzed through a preset impact improvement model, an improvement scheme is generated to predict and improve the impact effect, including:

[0129] S901, based on the first impact effect, the impact effect and the improvement direction are analyzed through a preset impact improvement model, and effect defect nodes and effect excess nodes are identified;

[0130] S902, according to the effect defect nodes and the effect excess nodes, a set of associated adjustment nodes is screened out;

[0131] S903, for each node in the set of adjustment nodes, the improvement parameters are analyzed, and the corresponding improvement scheme is generated to predict and improve the impact effect.

[0132] In this embodiment, based on the first impact effect, the impact effect and the improvement direction are analyzed through a preset impact improvement model, and the impact improvement model includes but is not limited to a support vector machine model pre-trained through a large amount of impact effect data. The model compares the predicted effect index of each node with the preset target threshold one by one, identifies the nodes whose predicted effect index does not meet the target requirement as effect defect nodes, reflects that the energy input at this place is insufficient or the action efficiency is low, and cannot achieve effective fragmentation; the nodes whose predicted effect index shows that the fragmentation degree far exceeds the target requirement are identified as effect excess nodes, reflecting that there is energy waste at this place; through the rapid screening of effect defect nodes and effect excess nodes, accurate improvement direction is provided for the subsequent effect improvement process, and the improvement effect and efficiency are improved.

[0133] Specifically, according to the effect defect nodes and the effect excess nodes, the blast zone action diagram is analyzed, for each effect defect node, the upstream nodes which are detonated before it and whose stress waves can propagate to the position and produce preloading effect are screened out; for each effect excess node, the strengthening source nodes which provide over-stress wave superposition are screened out; combined with the screened upstream nodes and strengthening source nodes, a set of adjustment nodes is constructed. By combining the blast zone action diagram to screen out the adjustment nodes, the effective adjustment entry point can be quickly screened out, the set of adjustment nodes screened out provides accurate improvement direction for the impact effect improvement process, improves the efficiency of the subsequent parameter optimization stage and the overall coordination of the improvement measures, avoids new negative reactions caused by local adjustment, and improves the improvement effect of the improvement scheme.

[0134] Specifically, for each node in the set of adjustment nodes, the improvement parameters are analyzed through a neural network model pre-trained through a large amount of historical improvement parameters, and the corresponding improvement scheme is generated to predict and improve the impact effect. By generating the improvement parameters corresponding to each node, the corresponding improvement scheme can be quickly obtained, the scheme development time of the impact optimization process is reduced, and the improvement effect and efficiency are improved.

[0135] As Figure 4 shown, the blast hole blockage length-based impact effect prediction improvement system for implementing the blast hole blockage length-based impact effect prediction improvement method comprises:

[0136] A blast hole blockage analysis module analyzes the stress distribution inside the blast hole according to the pre-acquired blast hole blockage detection data, determines the blast hole blockage length and the blockage state, and obtains the blast hole blockage analysis result.

[0137] An energy distribution analysis module, in combination with the blast hole blockage analysis result and the preset impact parameters, constructs an impact effect analysis model, predicts the dynamic distribution of energy in the blast hole by analyzing the stress wave propagation and gas expansion coupling in the blockage section, and obtains the energy distribution result.

[0138] A blast area action graph construction module, based on the energy distribution result, takes each blast hole as a node and the positional relationship between the blast holes as an edge to construct a blast area action graph.

[0139] An impact effect prediction module, according to the blast area action graph, simulates the propagation process of the impact stress wave and the crack by using the preset impact effect prediction model, analyzes the node state after the impact, and obtains the first impact effect.

[0140] An impact improvement module, based on the first impact effect, analyzes the improvement direction by using the preset impact improvement model, generates an improvement scheme, and predicts the impact effect.

[0141] In this embodiment, the blast hole blockage analysis module realizes accurate quantitative diagnosis of the blast hole blockage length and the internal density, particle size profile by fusing stress, radar, and acoustic wave and other multi-source detection data, provides high-precision and high-reliability initial conditions for subsequent analysis, and solves the problem of impact energy analysis error source caused by unknown blockage state; the energy distribution analysis module constructs a physical model coupling stress wave propagation and gas expansion effect according to the real physical properties of the blockage, dynamically predicts the effective distribution proportion of explosion energy in the axial direction of the blast hole, quantitatively predicts the energy regulation effect of the blockage section, can identify the enrichment area with high energy utilization rate and the leakage area with serious energy loss, and provides accurate data support for predicting the single-hole impact efficiency.

[0142] Specifically, the blast zone interaction graph construction module abstracts the entire blast zone as a network graph model with blast holes as nodes and interaction relationships as edges, and integrates the energy distribution characteristics of each hole, realizing the dimensionality of the perspective from single-hole analysis to multi-hole collaborative analysis, and providing data support for the quantitative study of stress wave superposition and crack penetration between holes; the impact effect prediction module dynamically simulates the whole process of stress wave propagation, superposition and rock mass progressive damage and crack propagation under multi-hole delay blasting based on the numerical model of the blast zone interaction graph, predicts the final fragmentation state of each region, can accurately analyze the spatial distribution of the overall fragmentation effect of the blast zone after impact, the fragmentation characteristics and potential problem areas, and improve the predictability and controllability of the impact process; the impact improvement module automatically diagnoses problem nodes based on the predicted effect, generates an impact parameter optimization adjustment scheme, and can quickly and accurately guide design iteration to improve impact effect and safety.

[0143] The above only describes the preferred embodiments of the present application, and the protection scope of the present application is not limited to the above-described embodiments. Any technical solutions falling within the scope of the present application should be considered as falling within the protection scope of the present application. It should be noted that, for ordinary skilled persons in the art, some improvements and refinements without departing from the principles of the present application should also be considered as falling within the protection scope of the present application.

Claims

1. A method for predicting and improving the impact effect based on borehole plugging length, characterized in that, include: Based on the pre-acquired detection data of blockages inside the borehole, the stress distribution inside the borehole is analyzed to determine the length and state of the blockage, and the borehole blockage analysis results are obtained. Based on the analysis results of borehole plugging and the preset impact parameters, an impact effect analysis model is constructed. By analyzing the stress wave propagation and gas expansion coupling in the plugged section, the dynamic distribution of energy in the borehole is predicted, and the energy distribution results are obtained. Based on the energy distribution results, each blast hole is treated as a node, and the positional relationship between blast holes is treated as an edge to construct a blast zone effect graph. Based on the explosion zone action diagram, the propagation process of impact stress wave and crack is simulated by a preset impact effect prediction model, and the node state after impact is analyzed to obtain the first impact effect. Based on the first impact effect, the impact effect and improvement direction are analyzed through a preset impact improvement model, and the defective and excessive effect nodes are identified. Based on the nodes with defective effects and the nodes with excessive effects, filter out the set of related adjustment nodes; For each node in the set of adjustment nodes, analyze the improvement parameters and generate corresponding improvement schemes to predict and improve the impact effect.

2. The method for predicting and improving the impact effect based on the borehole plugging length according to claim 1, characterized in that, The process involves analyzing the stress distribution inside the borehole based on pre-acquired data on borehole blockage detection, determining the borehole blockage length and state, and obtaining borehole blockage analysis results, including: Based on the pre-acquired detection data of blockages inside the borehole, the stress signal is analyzed to determine the stress changes in the direction perpendicular to the borehole wall and the direction parallel to the borehole wall, a stress distribution map is constructed, the position and distance between the blockage and the borehole opening are identified, and the length of the borehole blockage is obtained. The acoustic data characteristics corresponding to the borehole blockage length are analyzed, and a blockage profile along the borehole blockage length is constructed based on the acoustic data characteristics to obtain the blockage status. Based on the borehole blockage length and blockage status, the borehole blockage analysis results are obtained.

3. The method for predicting and improving the impact effect based on the borehole plugging length according to claim 2, characterized in that, The process involves analyzing the stress changes in the direction perpendicular to and parallel to the borehole wall based on pre-acquired data on borehole blockage detection, constructing a stress distribution map, identifying the location and distance between the blockage and the borehole opening, and obtaining the borehole blockage length. This includes: Based on the pre-acquired detection data of blockages inside the borehole, the stress signal is analyzed in the direction perpendicular to the borehole wall and in the direction parallel to the borehole wall, the normal stress sequence and the tangential stress sequence are calculated, and a stress distribution map is constructed. According to the preset distance interval, calculate the ratio of normal stress to tangential stress at each location point in the stress distribution diagram, analyze the change of the ratio along the borehole axis, and obtain the stress coupling ratio curve; Identify the inflection point in the stress coupling ratio curve and take the cross-section of the borehole corresponding to the inflection point as the first interface set; Based on the radar signals pre-acquired inside the borehole, the abrupt changes in the echo signal intensity are analyzed. The locations of radar abrupt changes are selected from the first interface set. The corresponding first interface is used as the distance between the blockage and the borehole opening, and the length of the borehole blockage is calculated.

4. The method for predicting and improving the impact effect based on the borehole plugging length according to claim 3, characterized in that, The analysis involves analyzing the acoustic data characteristics corresponding to the borehole blockage length, constructing a blockage profile along the borehole blockage length based on these characteristics, and analyzing the blockage state, including: Based on the acoustic data characteristics corresponding to the borehole plugging length, the energy attenuation ratio of the acoustic signal at each location point in the preset high-frequency sub-band relative to the energy of the acoustic emission signal is calculated to obtain the energy attenuation characteristic sequence. Calculate the offset of the center frequency of the acoustic signal along the axial direction at each location point to obtain the offset feature sequence; Interpolation processing is performed on the energy decay feature sequence and the offset feature sequence respectively to map them into a density sequence and a granularity sequence; Abrupt changes in the density sequence are identified, and the abrupt changes in the particle size sequence are analyzed at the abrupt change points. These abrupt changes are then divided into the first blockage region. The average density within the region is calculated, and a blockage profile is constructed along the blockage length of the borehole to analyze the blockage status.

5. The method for predicting and improving the impact effect based on the borehole plugging length according to claim 1, characterized in that, The impact effect analysis model is constructed by combining the borehole plugging analysis results and preset impact parameters. By analyzing the stress wave propagation and gas expansion coupling within the plugged section, the dynamic distribution of energy in the borehole is predicted, and the energy distribution results are obtained, including: Based on the density sequence in the borehole plugging analysis results, the impedance change and gas permeation are mapped and analyzed. Combined with the preset impact parameters, an impact effect analysis model is constructed. By analyzing the stress wave propagation and gas expansion coupling within the blocked section using an impact effect analysis model, the dynamic distribution of energy in the borehole is predicted, and the energy distribution results are obtained.

6. The method for predicting and improving the impact effect based on the borehole plugging length according to claim 5, characterized in that, The analysis of stress wave propagation and gas expansion coupling within the blocked section using an impact effect analysis model predicts the dynamic energy distribution within the borehole, yielding energy distribution results, including: The stress wave propagation process is simulated by an impact effect analysis model to obtain the energy accumulation. Based on the energy accumulation, the gas expansion and the process of pushing the blockage are simulated, and the gas expansion efficiency is calculated. Couple the energy accumulation situation with the gas expansion efficiency to calculate the energy accumulation distribution value; Based on the energy accumulation distribution value, the dynamic distribution of energy in the borehole is predicted, and the energy distribution result is obtained.

7. The method for predicting and improving the impact effect based on the borehole plugging length according to claim 1, characterized in that, The first impact effect is obtained by simulating the propagation process of impact stress wave and cracks using a preset impact effect prediction model based on the explosion zone action diagram, analyzing the post-impact node state, and including: Based on the explosion zone action diagram, analyze the crushing performance and impact time of each node, and calculate the dynamic stress intensity; Based on dynamic stress intensity, the propagation process of impact stress wave and crack is simulated by a preset impact effect prediction model, and the state of nodes after impact is analyzed to obtain the first impact effect.

8. The method for predicting and improving the impact effect based on the borehole plugging length according to claim 7, characterized in that, The first impact effect is obtained by simulating the propagation process of impact stress waves and cracks through a preset impact effect prediction model based on dynamic stress intensity, analyzing the nodal state after impact, and including: Based on dynamic stress intensity, the propagation process of impact stress waves and cracks is simulated by a preset impact effect prediction model according to the order of impact time, and the instantaneous stress caused by the stress wave generated by the node to the adjacent node is calculated. For each node, the impact effectiveness is calculated by combining the dynamic stress intensity of the node and the corresponding instantaneous stress, and the node state after impact is analyzed to obtain the impact damage state. Based on the impact damage state, the impact effect index is mapped to the corresponding node location, and the first impact effect is obtained through analysis.

9. A system for predicting and improving the impact effect based on borehole plugging length, characterized in that, The method for improving the prediction of impact effects based on borehole plugging length as described in any one of claims 1 to 8 includes: The borehole blockage analysis module analyzes the stress distribution inside the borehole based on the pre-acquired detection data of blockages inside the borehole, determines the length and state of the blockage, and obtains the borehole blockage analysis results. The energy distribution analysis module combines the results of borehole blockage analysis with preset impact parameters to construct an impact effect analysis model. By analyzing the stress wave propagation and gas expansion coupling within the blockage section, it predicts the dynamic distribution of energy in the borehole and obtains the energy distribution results. The blast zone effect graph construction module constructs a blast zone effect graph based on the energy distribution results, treating each blast hole as a node and the positional relationship between blast holes as edges. The impact effect prediction module simulates the propagation process of impact stress wave and cracks through a preset impact effect prediction model based on the explosion zone action diagram, analyzes the node state after impact, and obtains the first impact effect. The impact improvement module, based on the first impact effect, analyzes the improvement direction through a preset impact improvement model and generates an improvement plan to predict and improve the impact effect.

Images