A stability analysis method based on low-vibration blasting in cultural relic protection zone

Abstract

The present application relates to the technical field of stability analysis, and specifically discloses a kind of stability analysis method based on low vibration blasting of cultural relics protection zone, comprising the following steps: S1: monitoring each to be blasted area, S2: first blasting influence analysis, S3: second blasting influence analysis, S4: comprehensive blasting influence analysis, S5: blasting vibration influence analysis, S6: blasting pre-crack width design and S7: low vibration blasting stability analysis;The present application analyzes the first and second blasting influence coefficient, and then analyzes the comprehensive blasting influence coefficient;Use historical blasting data to analyze the blasting vibration influence coefficient;Based on the comprehensive blasting influence coefficient and the blasting vibration influence coefficient, design target pre-crack width;And the stability of low vibration blasting is evaluated.The present application realizes the accurate analysis of low vibration blasting of cultural relics protection zone, improves the scientific nature of blasting vibration stability performance evaluation, and guarantees the safety of cultural relics in the process of low vibration blasting of cultural relics protection zone.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of stability analysis technology, and particularly relates to a stability analysis method based on low-vibration blasting in a cultural relic protection zone. BACKGROUND

[0002] With the rapid development of urbanization and the continuous advancement of infrastructure construction, blasting engineering is increasingly carried out around the cultural relic protection zone. Cultural relics carry valuable historical and cultural value, and proper protection of them is an important mission to inherit national culture. Therefore, in the blasting operation, low-vibration blasting analysis is a core technical challenge to protect cultural relics, coordinate the relationship between engineering construction and cultural protection.

[0003] The traditional low-vibration blasting stability analysis method in the cultural relic protection zone mainly includes: calculating and limiting the maximum charge of a single section by an empirical formula to reduce the blasting vibration energy; using pre-splitting blasting technology to build a pre-splitting crack between the blasting area and the cultural relic protection zone to try to block the propagation of vibration waves; relying on engineering experience to optimize the initiation sequence to reduce the vibration superposition effect.

[0004] However, these methods have significant defects. First, the traditional method does not comprehensively analyze geological data and cultural relic stress parameters, making it difficult to fully assess the potential impact of blasting on geological conditions and cultural relic structures. Second, there is a lack of in-depth analysis of historical blasting data, resulting in inaccurate understanding of the characteristics of blasting vibration. In addition, the pre-splitting crack width design is not optimized in combination with the comprehensive blasting influence coefficient and the blasting vibration influence coefficient, and the low-vibration blasting stability performance evaluation lacks a systematic analysis of the blasting control parameters and the pre-splitting crack quality coefficient, making it difficult to accurately output the blasting stability performance evaluation index, unable to achieve accurate analysis of the whole process from blasting impact analysis, parameter design to effect evaluation, and difficult to meet the strict requirements of low-vibration blasting in the cultural relic protection zone. SUMMARY

[0005] In order to overcome the above-mentioned defects of the prior art, the embodiments of the present application provide a stability analysis method based on low-vibration blasting in a cultural relic protection zone to solve the problems raised in the background art.

[0006] To achieve the above-mentioned purpose, the present application provides the following technical scheme: a stability analysis method based on low-vibration blasting in a cultural relic protection zone, comprising the following steps: S1: monitoring each blasting area, S2: first blasting impact analysis, S3: second blasting impact analysis, S4: comprehensive blasting impact analysis, S5: blasting vibration impact analysis, S6: blasting pre-splitting crack width design, and S7: low-vibration blasting stability analysis.

[0007] S1: monitoring each blasting area: a monitoring network is arranged around the cultural relic protection zone, each blasting area of the cultural relic protection zone is simply denoted as each blasting area, and each blasting area is monitored in real time.

[0008] S2: First blasting influence analysis: Obtain the geological data around each blasting area, and analyze the first blasting influence coefficient of each blasting area based on the geological data;

[0009] S3: Second blasting influence analysis: Obtain the cultural relic body stress parameters around each blasting area, and analyze the second blasting influence coefficient of each blasting area based on the cultural relic body stress parameters;

[0010] S4: Comprehensive blasting influence analysis: Comprehensive analysis of the first blasting influence coefficient and the second blasting influence coefficient of each blasting area, and obtain the comprehensive blasting influence coefficient of each blasting area;

[0011] S5: Blasting vibration influence analysis: Obtain the historical particle vibration velocity, vibration duration and vibration acceleration of each blasting area, and analyze the blasting vibration influence coefficient of each blasting area;

[0012] S6: Blasting pre-crack width design: Analysis of the comprehensive blasting influence coefficient and the blasting vibration influence coefficient of each blasting area, and obtain the target pre-crack width that meets the requirements of each blasting area;

[0013] S7: Low-vibration blasting stability analysis: Based on the blasting control parameters and the pre-crack quality coefficient of each blasting area after blasting, the stability of low-vibration blasting of each blasting area is analyzed, the blasting stability performance evaluation index of each blasting area is obtained, and the performance evaluation is carried out, and the low-vibration blasting stability evaluation result is output.

[0014] Preferably, the monitoring of each blasting area is performed as follows:

[0015] First, a monitoring network composed of vibration sensors, displacement sensors and acoustic sensors is arranged around the cultural relic protection area, and the sensors are distributed around the cultural relic protection area in a grid shape, forming a density of 50m 2 A monitoring network of a measuring point;

[0016] Secondly, each blasting area of the cultural relic protection area is simply recorded as each blasting area, and each area is assigned a unique area code, and the each blasting area is sequentially labeled as 1, 2,..., i,..., n, i is the number of each blasting area.

[0017] Preferably, the first blasting influence analysis is performed as follows:

[0018] Obtain the geological data around each blasting area, including rock compressive strength, joint density and groundwater depth;

[0019] The actual compressive strength of the rock of the i-th to-be-blasted area Ct i The compressive strength deviation Sd of the rock of the i-th to-be-blasted area is calculated by subtracting the preset reference rock compressive strength Ct0 and taking the absolute value i ;

[0020] Based on the groundwater level depth, the groundwater level depth influence coefficient is analyzed, and the calculation formula is specifically as follows: Wherein, kw i represents the groundwater level depth influence coefficient of the i-th to-be-blasted area, Hw i represents the groundwater level depth of the i-th to-be-blasted area, Hw0 represents the critical groundwater level depth, and β represents the water pressure influence coefficient.

[0021] The correlation model of the geological data and the blasting influence is established, and the first blasting influence coefficient of each to-be-blasted area is calculated, and the calculation formula is specifically as follows:

[0022] Wherein, FIC i represents the first blasting influence coefficient of the i-th to-be-blasted area, Sd max represents the preset maximum allowable rock compressive strength deviation, Jd i represents the joint density of the i-th to-be-blasted area, and Jd0 represents the preset reference joint density.

[0023] Preferably, the execution mode of the second blasting influence analysis is specifically as follows:

[0024] The cultural relic body is comprehensively detected, and the cultural relic body stress parameters around each to-be-blasted area are obtained, and the cultural relic body stress parameters include cultural relic body structure allowable tensile stress, material elastic modulus and structure stress concentration coefficient.

[0025] Based on the cultural relic body stress parameter analysis, the second blasting influence coefficient of each to-be-blasted area is obtained, and the calculation formula is specifically as follows:

[0026]

[0027] Wherein, SIC i represents the second blasting influence coefficient of the i-th to-be-blasted area, St iz represents the i-th to-be-blasted area corresponding to the z-th cultural relic body structure allowable tensile stress, St iz ′ represents the i-th to-be-blasted area corresponding to the z-th cultural relic body preset reference structure allowable tensile stress, Em iz represents the i-th to-be-blasted area corresponding to the z-th cultural relic material elastic modulus, Em iz ′ represents the i-th to-be-blasted area corresponding to the z-th cultural relic preset reference material elastic modulus, and Cfiz represents the structural stress concentration coefficient of the zth cultural relic body corresponding to the ith blasting area, λ1, λ2, and λ3 represent the weight coefficients of the allowable tensile stress, the material elastic modulus, and the structural stress concentration coefficient of the cultural relic body structure respectively, z represents the number of each cultural relic body, and g represents the total number of cultural relic bodies.

[0028] Preferably, the execution mode of the comprehensive blasting influence analysis is specifically as follows:

[0029] The first blasting influence coefficient and the second blasting influence coefficient of each blasting area are obtained, and a comprehensive blasting influence coefficient of each blasting area is obtained through comprehensive analysis, and the calculation formula is specifically as follows:

[0030] IC i = w1xFIC i + w2xSIC i , wherein IC i represents the comprehensive blasting influence coefficient of the ith blasting area, FIC i represents the first blasting influence coefficient of the ith blasting area, SIC i represents the second blasting influence coefficient of the ith blasting area, w1 and w2 represent the weight coefficients of the first blasting influence coefficient and the second blasting influence coefficient respectively, and w1+w2=1.

[0031] Preferably, the execution mode of the blasting vibration influence analysis is specifically as follows:

[0032] The historical blasting data of each blasting area is obtained, including the particle vibration velocity, the vibration duration, and the vibration acceleration, the blasting vibration influence coefficient of each blasting area is calculated, and the calculation formula is as follows:

[0033] , wherein IV i represents the blasting vibration influence coefficient of the ith blasting area, v ij represents the particle vibration velocity measured in the jth historical blasting corresponding to the ith blasting area, v th represents a preset vibration velocity safety threshold, T ij represents the vibration duration measured in the jth historical blasting corresponding to the ith blasting area, T th represents a preset vibration duration safety threshold, a ij represents the vibration acceleration measured in the jth historical blasting corresponding to the ith blasting area, a th represents a preset vibration acceleration safety threshold, and kf ijThe frequency correction coefficient of the jth historical blasting corresponding to the ith blasting area is represented, m represents the total number of historical blasting data, j represents the number of each historical blasting data, j=1, 2, 3, …, m, i is the number of each blasting area, and a1, a2, and a3 represent the weight indexes of the particle vibration speed, vibration duration, and vibration acceleration, respectively.

[0034] Preferably, the execution mode of the blasting pre-crack width design is specifically as follows:

[0035] The comprehensive blasting influence coefficient and the blasting vibration influence coefficient of each blasting area are obtained, analyzed, and the target pre-crack width of each blasting area is obtained. The calculation formula of the target pre-crack width of each blasting area is specifically as follows:

[0036] Wherein, Wp i represents the target pre-crack width of the ith blasting area, W0 represents the preset reference pre-crack width, IC i represents the comprehensive blasting influence coefficient of the ith blasting area, IV i represents the blasting vibration influence coefficient of the ith blasting area, k represents a constant, q i represents the linear charge density of the ith blasting area, q0 represents the preset reference linear charge density, Sp i represents the pre-crack hole spacing of the ith blasting area, S0 represents the preset reference pre-crack hole spacing, represents a correction coefficient.

[0037] Preferably, the execution mode of the low-vibration blasting stability analysis is specifically as follows:

[0038] The blasting control parameters after blasting of each blasting area are obtained, and the blasting control parameters include the actual vibration speed and the main vibration frequency;

[0039] The blasting stability performance evaluation index of each blasting area is calculated, and the calculation formula is specifically as follows:

[0040] Wherein, EEI i represents the blasting stability performance evaluation index of the ith blasting area, V max represents the preset safe maximum vibration speed, Vv i represents the actual vibration speed of the ith blasting area, f i represents the main vibration frequency of the ith blasting area, fn represents the natural frequency of the protection target, fr represents the preset frequency threshold, wherein fr=2fn, Qf i represents the pre-crack quality coefficient of the ith blasting area, and e represents a natural constant.

[0041] Based on the blasting stability performance evaluation index of each to be blasted area, the stability performance evaluation is carried out, and the specific content of the stability performance evaluation is as follows:

[0042] The blasting stability performance evaluation index of each to be blasted area is read, which is compared with the preset stability performance evaluation index threshold value, if the blasting stability performance evaluation index of a to be blasted area is greater than the preset stability performance evaluation index threshold value, it is determined that the blasting stability state of the to be blasted area is normal, otherwise, it is determined that the blasting stability state of the to be blasted area is abnormal, the blasting stability state abnormal data of the to be blasted area is marked as the low vibration blasting stability performance evaluation result of the to be blasted area, and the low vibration blasting stability performance evaluation result is output.

[0043] Preferably, the calculation formula of the pre-crack quality coefficient of the i th to be blasted area is specifically as follows:

[0044] Wherein, Qf i represents the pre-crack quality coefficient of the i th to be blasted area, Le i represents the effective pre-crack length of the i th to be blasted area, Le i ' represents the preset theoretical total length of pre-crack of the i th to be blasted area, Dd i represents the preset allowable pre-crack surface deviation angle of the i th to be blasted area, Da i represents the actual pre-crack surface deviation angle of the i th to be blasted area, Wa i represents the actual pre-crack width of the i th to be blasted area, Wp i represents the target pre-crack width of the i th to be blasted area.

[0045] As described above, the low vibration blasting stability analysis method provided by the application has at least the following beneficial effects:

[0046] The application provides a stability analysis method based on low-vibration blasting in a cultural relic protection area. BRIEF DESCRIPTION OF DRAWINGS

[0047] The application is further described by using the drawings, but the embodiments in the drawings do not constitute any limitation on the application, and other drawings can be obtained by those skilled in the art without creative labor on the basis of the following drawings.

[0048] Figure 1 FIG. 1 is a flowchart of the stability analysis method based on low-vibration blasting in a cultural relic protection area. DETAILED DESCRIPTION

[0049] The technical solutions in the embodiments of the application will be clearly and completely described below with reference to the drawings in the embodiments of the application. Obviously, the described embodiments are only part of the embodiments of the application, rather than all the embodiments of the application. Based on the embodiments in the application, all other embodiments obtained by those skilled in the art without creative labor fall within the protection scope of the application.

[0050] Please refer to Figure 1 The application provides a stability analysis method based on low-vibration blasting in a cultural relic protection area, which comprises the following steps: S1: monitoring each blasting area, S2: first blasting impact analysis, S3: second blasting impact analysis, S4: comprehensive blasting impact analysis, S5: blasting vibration impact analysis, S6: blasting pre-crack width design, and S7: low-vibration blasting stability analysis.

[0051] S1: monitoring each blasting area: a monitoring network is arranged around the cultural relics protection area, each blasting area of the cultural relics protection area is simply recorded as each blasting area, and each blasting area is monitored in real time;

[0052] In this embodiment, it needs to be specifically pointed out that the execution mode of monitoring each blasting area is specifically as follows:

[0053] Firstly, a monitoring network composed of vibration sensors, displacement sensors and acoustic sensors is arranged around the cultural relics protection area; specifically, the vibration sensor adopts a three-way acceleration sensor with a sensitivity of 100 mV / g and a frequency response range of 0.5 Hz-5 kHz; the displacement sensor adopts a laser displacement sensor with an accuracy of 0.01 mm; the acoustic sensor adopts a sound emission sensor with a frequency range of 50 kHz-300 kHz; the sensors are distributed around the cultural relics protection area in a grid shape, forming a density of 50 m 2 The monitoring network of one measuring point;

[0054] Secondly, each blasting area of the cultural relics protection area is simply recorded as each blasting area, and each area is assigned a unique area code, and the each blasting area is sequentially labeled as 1, 2,..., i,..., n, i is the number of each blasting area;

[0055] Finally, the warning threshold is set: when the monitored vibration speed exceeds 0.5 cm / s or the displacement exceeds 0.2 mm, the system automatically sends a warning signal, the monitoring system operates all day long, records data once per second, and stores the data in the database for subsequent analysis.

[0056] S2: first blasting influence analysis: obtaining geological data around each blasting area, and analyzing the first blasting influence coefficient of each blasting area based on the geological data;

[0057] In this embodiment, it needs to be specifically pointed out that the execution mode of the first blasting influence analysis is specifically as follows:

[0058] Obtaining geological data around each blasting area, the geological data including rock compressive strength, joint density and groundwater depth;

[0059] The actual rock compressive strength Ct i of the i-th blasting area is obtained by subtracting the preset reference rock compressive strength Ct0 and taking the absolute value, and the rock compressive strength deviation Sd i of the i-th blasting area is calculated.

[0060] Based on the groundwater depth, the groundwater depth influence coefficient is analyzed, and the calculation formula is specifically as follows: Where, kw iHw represents the influence coefficient of groundwater level depth in the i-th area to be blasted. i Hw0 represents the groundwater level depth of the i-th area to be blasted, Hw0 represents the critical groundwater level depth, and β represents the water pressure influence coefficient, which takes a value between 0.005 and 0.015.

[0061] It should be noted that the depth of the groundwater level affects the performance of explosives and the blasting effect through water pressure. Increased water depth reduces the detonation velocity and intensity of the explosives.

[0062] A correlation model between the geological data and the blasting impact is established, and the first blasting impact coefficient for each area to be blasted is calculated. The specific calculation formula is as follows:

[0063] Among them, FIC i Sd represents the first blasting influence coefficient of the i-th area to be blasted. max Jd represents the preset maximum allowable deviation in rock compressive strength. i Jd0 represents the joint density of the i-th region to be blasted, and Jd0 represents the preset reference joint density.

[0064] In this embodiment, it should be specifically noted that the rock compressive strength deviation Sd of the i-th region to be blasted in the formula... i Larger joint density deviation | Jd i -Jd0| The larger the value, the greater the influence coefficient of groundwater level depth (kw) i The larger the value, the greater the first blasting influence coefficient FIC of the i-th blasting area. i The larger the value, the worse the overall adverse impact of the geological conditions in the area to be blasted on the blasting. i It reflects the degree of influence of groundwater level depth on blasting; the larger the value, the more significant the influence.

[0065] S3: Second blasting impact analysis: Obtain the stress parameters of the cultural relics around each area to be blasted, and analyze the second blasting impact coefficient of each area to be blasted based on the stress parameters of the cultural relics.

[0066] In this embodiment, it should be specifically explained that the execution method of the second blasting impact analysis is as follows:

[0067] A comprehensive inspection of the cultural relic was conducted to obtain the stress parameters of the cultural relic around each area to be blasted. The stress parameters of the cultural relic included the allowable tensile stress of the cultural relic structure, the elastic modulus of the material, and the stress concentration factor of the structure.

[0068] The detection methods used for comprehensive testing include: non-contact strain measurement, ultrasonic testing, and infrared thermography.

[0069] A second blasting influence coefficient of each to-be-blasted region is obtained based on the cultural relic body stress parameter analysis, and the calculation formula is specifically as follows:

[0070]

[0071] wherein, SIC i represents the second blasting influence coefficient of the i-th to-be-blasted region, St iz represents the z-th cultural relic body structure allowable tensile stress corresponding to the i-th to-be-blasted region, St iz ' represents the z-th cultural relic body preset reference structure allowable tensile stress corresponding to the i-th to-be-blasted region, Em iz represents the z-th cultural relic material elastic modulus corresponding to the i-th to-be-blasted region, Em iz ' represents the z-th cultural relic preset reference material elastic modulus corresponding to the i-th to-be-blasted region, Cf iz represents the z-th cultural relic body structure stress concentration coefficient corresponding to the i-th to-be-blasted region, λ1, λ2 and λ3 represent weight coefficients of the cultural relic body structure allowable tensile stress, the material elastic modulus and the structure stress concentration coefficient respectively, z represents the number of each cultural relic body, and g represents the total number of the cultural relic bodies;

[0072] In this embodiment, it needs to be specifically explained that the structure stress concentration coefficient of the cultural relic body is obtained in the following manner:

[0073] The three-dimensional laser scanning and photographic measurement technology is used to obtain the cultural relic geometric data, a three-dimensional digital model is constructed, the material elastic modulus, Poisson's ratio and compressive / tensile strength parameters of the cultural relic are obtained through micro-damage detection (such as ultrasonic testing and point load testing), the finite element software (such as ANSYS and ABAQUS) is input, the finite element analysis is run, the maximum principal stress of the stress concentration part is determined, the nominal stress of the region is calculated through the integral path, the maximum principal stress of the stress concentration part is divided by the nominal stress of the region, and the structure stress concentration coefficient of the cultural relic body is calculated, wherein the greater the structure stress concentration coefficient value of the cultural relic body is, the more significant the local stress concentration is.

[0074] It needs to be specifically explained that in one specific embodiment, λ1 can be set to 0.5, λ2 can be set to 0.3, and λ3 can be set to 0.2.

[0075] S4: Comprehensive blasting influence analysis: the first blasting influence coefficient and the second blasting influence coefficient of each to-be-blasted region are comprehensively analyzed to obtain the comprehensive blasting influence coefficient of each to-be-blasted region;

[0076] In this embodiment, it needs to be specifically explained that the execution manner of the comprehensive blasting influence analysis is specifically as follows:

[0077] The first blasting influence coefficient and the second blasting influence coefficient of each to-be-blown region are acquired, and comprehensive analysis is performed thereon to obtain the comprehensive blasting influence coefficient of each to-be-blown region, and the calculation formula is as follows:

[0078] IC i =w1×FIC i +w2×SIC i , wherein, IC i represents the comprehensive blasting influence coefficient of the i-th to-be-blown region, FIC i represents the first blasting influence coefficient of the i-th to-be-blown region, SIC i represents the second blasting influence coefficient of the i-th to-be-blown region, w1 and w2 represent the weight coefficients of the first blasting influence coefficient and the second blasting influence coefficient respectively, and w1+w2=1.

[0079] It should be specifically noted that, in a specific embodiment, w1 is 0.4 and w2 is 0.6.

[0080] S5: Blasting vibration influence analysis: the historical particle vibration velocity, vibration duration and vibration acceleration of each to-be-blown region are acquired, and the blasting vibration influence coefficient of each to-be-blown region is obtained through analysis;

[0081] In this embodiment, it should be specifically noted that the execution mode of the blasting vibration influence analysis is as follows:

[0082] The historical blasting data of each to-be-blown region, including the particle vibration velocity, vibration duration and vibration acceleration, are acquired, and the blasting vibration influence coefficient of each to-be-blown region is calculated, and the calculation formula is as follows:

[0083] , wherein, IV i represents the blasting vibration influence coefficient of the i-th to-be-blown region, v ij represents the particle vibration velocity measured in the j-th historical blasting of the i-th to-be-blown region, v th represents a preset vibration velocity safety threshold, T ij represents the vibration duration measured in the j-th historical blasting of the i-th to-be-blown region, T th represents a preset vibration duration safety threshold, a ij represents the vibration acceleration measured in the j-th historical blasting of the i-th to-be-blown region, a th represents a preset vibration acceleration safety threshold, kf ijindicates the frequency correction coefficient of the jth historical blasting corresponding to the ith blasting area to be blasted, m indicates the total number of historical blasting data, j indicates the number of each historical blasting data, j = 1, 2, 3, …, m, i is the number of each blasting area, and α1, α2, and α3 respectively indicate the weight index of the particle vibration velocity, the vibration duration, and the vibration acceleration;

[0084] It needs to be specifically pointed out that in a specific embodiment, α1 indicates the weight index of the particle vibration velocity in the formula, generally taking a value range of 1.5-2.0, here taking 1.8, for highlighting the contribution degree of the particle vibration velocity to the blasting vibration influence coefficient; α2 indicates the weight index of the vibration duration in the formula, generally taking a value range of 1.2-1.8, here taking 1.5; and α3 indicates the weight index of the vibration acceleration in the formula, generally taking a value range of 1.0-1.5, here taking 1.2.

[0085] For the frequency correction coefficient kf ij of the jth historical blasting corresponding to the ith blasting area to be blasted, ij when kf 文物 is within the range of ±20% of the cultural relic natural vibration frequency f ij , kf ij =1.5, otherwise, kf i =1.0.

[0086] S6: blasting pre-crack width design: analyzing the comprehensive blasting influence coefficient and the blasting vibration influence coefficient of each blasting area to be blasted to obtain the target pre-crack width meeting the each blasting area to be blasted;

[0087] In the embodiment, it needs to be specifically pointed out that the execution mode of the blasting pre-crack width design is specifically as follows:

[0088] The comprehensive blasting influence coefficient and the blasting vibration influence coefficient of each blasting area to be blasted are obtained, and are analyzed to obtain the target pre-crack width meeting the each blasting area to be blasted, and the calculation formula of the target pre-crack width of the each blasting area to be blasted is specifically as follows:

[0089] Wherein, Wp i indicates the target pre-crack width of the ith blasting area to be blasted, W0 indicates a preset reference pre-crack width, IC i indicates the comprehensive blasting influence coefficient of the ith blasting area to be blasted, IV i indicates the blasting vibration influence coefficient of the ith blasting area to be blasted, k indicates a constant, q i indicates the linear charge density of the ith blasting area to be blasted, q0 indicates a preset reference linear charge density, and Sp iS0 represents a preset reference presplit hole spacing, represents a correction coefficient.

[0090] S7: Low-vibration blasting stability analysis: based on the blasting control parameters and the presplit crack quality coefficient of each blasting area after blasting, the low-vibration blasting stability of each blasting area is analyzed, the blasting stability performance evaluation index of each blasting area is obtained, and the performance evaluation is performed, and the low-vibration blasting stability evaluation result is output.

[0091] In this embodiment, it needs to be specifically pointed out that the execution mode of the low-vibration blasting stability analysis is as follows:

[0092] The blasting control parameters after blasting of each blasting area are obtained, and the blasting control parameters include the actual vibration speed and the main vibration frequency;

[0093] The blasting stability performance evaluation index of each blasting area is calculated, and the calculation formula is as follows:

[0094] EEI i represents the blasting stability performance evaluation index of the i-th blasting area, V max represents a preset safety allowable maximum vibration speed, Vv i represents the actual vibration speed of the i-th blasting area, f i represents the main vibration frequency of the i-th blasting area, fn represents the natural frequency of the protection target, and fr represents a preset frequency threshold, wherein fr = 2fn, Qf i represents the presplit crack quality coefficient of the i-th blasting area, and e represents a natural constant;

[0095] It needs to be specifically pointed out that in the formula, the smaller the actual vibration speed Vv i of the i-th blasting area, the greater the difference |f i -fn| of the main vibration frequency of the i-th blasting area and the natural frequency of the protection target, and the greater the presplit crack quality coefficient Qf i of the i-th blasting area, the greater the blasting stability performance evaluation index of the i-th blasting area, indicating that the blasting of the blasting area is more stable.

[0096] It needs to be specifically pointed out that the natural frequency fn of the protection target refers to the inherent vibration frequency (unit: Hz) of the protection target (such as ancient buildings, sculptures), which is determined by the structural material and geometric shape (such as the natural frequency of wood structure is usually 1-5 Hz), and when the main vibration frequency f i of blasting is farther away from fn, resonance damage is less likely to occur.

[0097] Based on the blasting stability performance evaluation index of each to-be-blasted region, the stability performance of each to-be-blasted region is evaluated, and the specific content of the evaluation is as follows:

[0098] The blasting stability performance evaluation index of each to-be-blasted region is read, and compared with the preset stability performance evaluation index threshold. If the blasting stability performance evaluation index of a to-be-blasted region is greater than the preset stability performance evaluation index threshold, it is determined that the blasting stability state of the to-be-blasted region is normal. Otherwise, it is determined that the blasting stability state of the to-be-blasted region is abnormal, the blasting stability state abnormal data of the to-be-blasted region is marked as the low-vibration blasting stability performance evaluation result of the to-be-blasted region, and the low-vibration blasting stability performance evaluation result is output.

[0099] In this embodiment, it needs to be specifically pointed out that the calculation formula of the pre-crack quality coefficient of the i-th to-be-blasted region is as follows:

[0100] Wherein, Qf i represents the pre-crack quality coefficient of the i-th to-be-blasted region, Le i represents the effective pre-crack length of the i-th to-be-blasted region, which refers to the pre-crack length that is actually formed and can effectively block the propagation of blasting vibration, Le i ′ represents the preset theoretical total pre-crack length of the i-th to-be-blasted region, Dd i represents the preset allowable pre-crack surface deviation angle of the i-th to-be-blasted region, Da i represents the actual pre-crack surface deviation angle of the i-th to-be-blasted region, Wa i represents the actual pre-crack width of the i-th to-be-blasted region, Wp i represents the target pre-crack width of the i-th to-be-blasted region.

[0101] Finally: the above only for the preferred embodiments of the present application, and not for limiting the present application, any modification, equivalent replacement, improvement, etc. made within the spirit and principles of the present application, should be included in the protection scope of the present application.

[0102] The above is only a specific embodiment of the present application, but the protection scope of the present application is not limited to this. Any person skilled in the art can easily think of changes or replacements within the technical range disclosed in the present application, which should be included in the protection scope of the present application. Therefore, the protection scope of the present application should be subject to the protection scope of the claims.

Claims

1. A stability analysis method based on low-vibration blasting in a cultural relic protection zone, characterized by, Comprising the following steps: S1: monitoring each blasting area: a monitoring network is laid around the cultural relics protection area, each blasting area of the cultural relics protection area is simply recorded as each blasting area, and each blasting area is monitored in real time; S2: first blasting influence analysis: obtain the geological data around each blasting area, and analyze the first blasting influence coefficient of each blasting area based on the geological data; S3: second blasting influence analysis: obtain the cultural relics body stress parameters around each blasting area, and analyze the second blasting influence coefficient of each blasting area based on the cultural relics body stress parameters; S4: comprehensive blasting influence analysis: comprehensively analyze the first blasting influence coefficient and the second blasting influence coefficient of each blasting area to obtain the comprehensive blasting influence coefficient of each blasting area; S5: blasting vibration influence analysis: obtain the historical particle vibration velocity, vibration duration and vibration acceleration of each blasting area, and analyze the blasting vibration influence coefficient of each blasting area; S6: blasting pre-crack width design: analyze the comprehensive blasting influence coefficient and the blasting vibration influence coefficient of each blasting area to obtain the target pre-crack width meeting the each blasting area; S7: low-vibration blasting stability analysis: based on the blasting control parameters and the pre-crack quality coefficient of each blasting area after blasting, the stability of low-vibration blasting of each blasting area is analyzed to obtain the blasting stability performance evaluation index of each blasting area, and the performance is evaluated, and the low-vibration blasting stability evaluation result is outputted; The calculation formula of the pre-crack quality coefficient of the blasting area is as follows: wherein, Qf i a preset theoretical pre-splitting total length of the i-th blasting area to be blasted, Le i an effective pre-splitting length of the i-th blasting area to be blasted, a preset theoretical pre-splitting total length of the i-th blasting area to be blasted, Dd i a preset allowable pre-splitting surface deviation angle of the i-th blasting area to be blasted, Da i an actual pre-splitting surface deviation angle of the i-th blasting area to be blasted, Wa i an actual pre-splitting width of the i-th blasting area to be blasted, Wp i a target pre-splitting width of the i-th blasting area to be blasted.

2. The method according to claim 1, characterized in that: The execution mode of monitoring each blasting area is as follows: First, a monitoring network composed of vibration sensors, displacement sensors and acoustic sensors is laid around the cultural relics protection area, and the sensors are distributed in a grid shape around the cultural relics protection area to form a monitoring network with a density of one measuring point per 50 square meters; Secondly, each blasting area of the cultural relics protection area is simply recorded as each blasting area, and each area is assigned a unique area code, and the each blasting area is sequentially labeled as 1, 2,..., i,..., n, i is the number of each blasting area.

3. The method according to claim 1, characterized in that: The execution mode of the first blasting influence analysis is as follows: Obtain the geological data around each blasting area, the geological data including rock compressive strength, joint density and groundwater depth; The actual compressive strength of the rock in the i-th blasting area is obtained The preset reference rock compressive strength The rock compressive strength deviation of the i-th blasting area is calculated by subtracting and taking the absolute value ; Based on the groundwater level depth, the groundwater level depth influence coefficient is analyzed, and the calculation formula is specifically as follows: wherein, kw i denotes the groundwater level depth influence coefficient of the i-th to-be-blasted area, Hw i denotes the groundwater level depth of the i-th to-be-blasted area, Hw 0 denotes the critical groundwater level depth, denotes the water pressure influence coefficient; An association model of the geological data and blasting influence is established, and the first blasting influence coefficient of each blasting area is calculated, and the calculation formula is as follows: wherein, FIC i represents a first blasting influence coefficient of the i-th blasting area to be blasted, represents a preset allowed maximum deviation of rock compressive strength, Jd i represents a joint density of the i-th blasting area to be blasted, Jd 0 represents a preset reference joint density.

4. The method according to claim 1, characterized in that: The execution mode of the second blasting influence analysis is as follows: Comprehensively detect the cultural relics body, obtain the cultural relics body stress parameters around each blasting area, and the cultural relics body stress parameters include cultural relics body structure allowable tensile stress, material elastic modulus and structure stress concentration coefficient; Based on the cultural relics body stress parameters, the second blasting influence coefficient of each blasting area is analyzed, and the calculation formula is as follows: wherein, SIC i a second blasting influence coefficient of the i-th region to be blasted, a structure allowable tensile stress of the z-th cultural relic body corresponding to the i-th region to be blasted, a preset reference structure allowable tensile stress of the z-th cultural relic body corresponding to the i-th region to be blasted, an elastic modulus of the z-th cultural relic material corresponding to the i-th region to be blasted, a preset reference elastic modulus of the z-th cultural relic material corresponding to the i-th region to be blasted, a structure stress concentration coefficient of the z-th cultural relic body corresponding to the i-th region to be blasted, respectively represent weight coefficients of the structure allowable tensile stress, the elastic modulus, and the structure stress concentration coefficient of the cultural relic body, z represents the number of each cultural relic body, and g represents the total number of cultural relic bodies.

5. The method according to claim 1, wherein: The execution mode of the comprehensive blasting influence analysis is as follows: The first blasting influence coefficient and the second blasting influence coefficient of each blasting area are obtained, and comprehensive analysis is performed thereon to obtain the comprehensive blasting influence coefficient of each blasting area, and the calculation formula is as follows: wherein, IC i represents a first blast impact coefficient of the i-th region to be blasted, FIC i represents a first blast impact coefficient of the i-th region to be blasted, SIC i represents a second blast impact coefficient of the i-th region to be blasted, represent weight coefficients of the first blast impact coefficient and the second blast impact coefficient, respectively, and .

6. The method according to claim 1, wherein: The execution mode of the blasting vibration influence analysis is as follows: The historical blasting data of each blasting area are obtained, including the particle vibration velocity, vibration duration and vibration acceleration, the blasting vibration influence coefficient of each blasting area is calculated, and the calculation formula is as follows: wherein, IV i represents the blasting vibration influence coefficient of the i-th blasting area, represents the particle vibration velocity of the i-th blasting area corresponding to the j-th historical blasting measurement, represents the preset vibration velocity safety threshold, represents the vibration duration of the i-th blasting area corresponding to the j-th historical blasting measurement, represents the preset vibration duration safety threshold, represents the vibration acceleration of the i-th blasting area corresponding to the j-th historical blasting measurement, represents the preset vibration acceleration safety threshold, represents the frequency correction coefficient of the i-th blasting area corresponding to the j-th historical blasting, m represents the total number of historical blasting data, j represents the number of each historical blasting data, j = 1, 2, 3, …, m, i is the number of each blasting area, respectively represent the weight index of particle vibration velocity, vibration duration, and vibration acceleration.

7. The method according to claim 1, wherein the method is characterized by: The execution mode of the blasting pre-crack width design is as follows: The comprehensive blasting influence coefficient and the blasting vibration influence coefficient of each blasting area are obtained, and analysis is performed thereon to obtain the target pre-crack width of each blasting area, and the calculation formula of the target pre-crack width of each blasting area is as follows: wherein, Wp i a target pre-splitting width of the i-th blasting area, W 0 a preset reference pre-splitting width, IC i a comprehensive blasting influence coefficient of the i-th blasting area, IV i a blasting vibration influence coefficient of the i-th blasting area, k represents a constant, a linear charge density of the i-th blasting area, a preset reference linear charge density, a pre-split hole spacing of the i-th blasting area, a preset reference pre-split hole spacing, a correction coefficient.

8. The method according to claim 1, characterized in that: The execution mode of the low-vibration blasting stability analysis is as follows: The blasting control parameters after blasting of each blasting area are obtained, and the blasting control parameters include the actual vibration velocity and the main vibration frequency; The blasting stability performance evaluation index of each blasting area is calculated, and the calculation formula is as follows: wherein, EEI i represents the blasting stability performance evaluation index of the i-th blasting area, represents the preset safety-allowed maximum vibration velocity, Vv i represents the actual vibration velocity of the i-th blasting area, f i represents the main vibration frequency of the i-th blasting area, fn represents the natural frequency of the protection target, fr represents the preset frequency threshold, wherein fr =2 fn , Qf i represents the pre-crack quality coefficient of the i-th blasting area, and e represents the natural constant; Based on the blasting stability performance evaluation index of each blasting area, the stability performance evaluation is performed, and the specific content of the stability performance evaluation is as follows: The blasting stability performance evaluation index of each blasting area is read, and is compared with the preset stability performance evaluation index threshold value, if the blasting stability performance evaluation index of a certain blasting area is greater than the preset stability performance evaluation index threshold value, it is determined that the blasting stability state of the blasting area is normal, otherwise, it is determined that the blasting stability state of the blasting area is abnormal, the blasting stability state abnormal data of the blasting area is marked as the low-vibration blasting stability performance evaluation result of the blasting area, and the low-vibration blasting stability performance evaluation result is output.

Images