# Asymptotic Structure of the Seismic Radiation from an Explosive Column.

1. IntroductionThe radiation of elastic waves produced by a cylindrical charge in a borehole has important applications in seismic exploration, mining geophysics, and blast vibration prediction. Today it is a well-known fact that, in the blasting engineering, the curves of blasting vibration in the near field and the far field show that vibration amplitudes induced by P wave and S wave are generally in the same magnitude, so the effect of S wave could not be ignored.

A pioneering theoretical work is that developed by Heelan (1953) [1, 2], where he study the radiation from a cylindrical source of finite length, the walls of which are subjected to symmetric lateral and tangential stresses. His papers were credited with demonstrating the generation of shear waves by artificial sources in a borehole; at the time this was a fundamentally surprising result. Subsequently, several authors have been making improvements to the results of Heelan, principally by Jordan [3], Hazebroek [4], and Abo-Zena [5]. One of the main advantages of the Heelan solution is its algebraic simplicity and clear physical interpretation, which allows easy implementation in more complex models.

An exact analytical solution of the problem solved by Heelan [2] was developed by Meredith [6]; Meredith et al. [7] and Meredith et al. [8] analyze the exact full-field solution for the axisymmetric wave radiation caused by any time-varying pressure load source acting over a finite length of a borehole. Tubman [9] and Tubman et al. [10] derive a general expression for the dispersion relations and the impulse response of a radially layered borehole.

Blair and Minchinton [11] analyze the near field of vibration from an explosive column; Blair [12] makes a comparison between Heelan solution and exact solutions for seismic radiation from a short cylindrical charge. Analytical solutions for describing Mach waves produced by infinitely long explosive columns were studied by Rossmanith and Kouzniak [13]. Blair [14] studied the blast vibration dependence on charge length, velocity of detonation, and layered media. Trivino et al. [15] studied seismic radiation patterns from cylindrical explosive charges by analytical and combined finite-discrete element methods. Kumar et al. [16] analyze the blast-induced ground vibration equations for rocks using mechanical and geological properties. Lee and Balch [17] generalize the results of Heelan to a fluid-filled borehole. The interaction between two adjacent blast holes has been studied by Yi et al. [18]. Trivino et al. [19] have studied the seismic waveforms from explosive sources located in boreholes and initiated in different directions. Using the Heelan results, Chen et al. [20] have studied the influence of millisecond time, charge length, and detonation velocity on blasting vibration.

In this work we study the structure of the seismic radiation in the far field produced by an explosive column. Using an asymptotic solution for the far field of vibration (Heelan's solution), we find analytical expressions to the peak particle velocity (PPV) diagrams. The analysis is extended to the case of a charge with finite velocity of detonation, obtaining a first correction to the Heelan's solution in the case of finite VOD. The results are compared with those existing in the literature.

2. Mathematical Model

2.1. Short Cylindrical Charge (Heelan's Solution). Consider a cylindrical explosive charge length 21 and radius a, centered at the origin of coordinates and embedded in an infinite elastic medium (Figure 1). The vertical walls of this cylinder are subjected to stresses which are symmetric about the vertical axis, uniform in the sense that they have the same instantaneous value at all points, and finite in time duration.

It has been shown (Heelan 1953) that for a pressure transient acting on the walls of a short section of an infinitely long cylindrical hole of radius a, the far field radial and tangential displacements are

[mathematical expression not reproducible]. (1)

To simplify the notation we will define the coefficients [F.sub.1]([phi]) and [F.sub.2]([phi]):

[mathematical expression not reproducible], (2)

where 2l is the total length of the loaded section of the cylinder, [phi] is the polar angle measured from the cylinder's axis, R is the distance to the observer, [c.sub.s] = [square root of [mu]/[rho]] and [c.sub.p] = [square root of ([lambda] + 2[mu])/[rho]] are the shear and compressional wave speeds in the medium, and [mu] and [rho] are the rock's rigidity and density, respectively. In this analysis, the source pressure function is based on the experimental data of Grady et al. [21] and Larson [22] and is given by

P(t) = [P.sub.VN][(e[gamma]/n).sup.n][t.sup.n][e.sup.-[gamma]t]H(t), (3)

where "e" is the Euler's constant, [P.sub.VN] is the von Neumann borehole pressure at the detonation front produced by the explosive, H(t) is the Heaviside unit step function, n is an integer, and [gamma] is a pressure decay parameter. The values of the parameter n are in the range of 3 to 9.

Equations (1) show that, in the far field, radial and tangential displacements travel as waves in the solid (S wave and P wave), with speeds [c.sub.s] and [c.sub.p], respectively, where [c.sub.p] > [c.sub.s]. Another important feature of the solution of Heelan is the angular variation of the amplitude of the radiated P and S waves, where the displacement induced by shear wave could be larger than that caused by P wave. Figure 2 shows the angular distribution of P and S waves.

You can see that there are regions where the P waves dominate over the S waves and vice versa. Clearly the S waves predominate around [[phi].sub.n] = (2n + 1)([pi]/4), with n = 0, 1, 2, 3. The P waves predominate around [[phi].sub.n] = n([pi]/2). The region where the P waves predominate over S waves can be estimated from the relationship [F.sub.1]([phi]) [greater than or equal to] [F.sub.2]([phi]); then

[mathematical expression not reproducible]. (4)

This equation leads to the following expression:

[mathematical expression not reproducible], (5)

where [bar.[lambda]] = [c.sub.s]/[c.sub.p] < 1. The angular widths for windows radiation P waves are given by

[DELTA][[phi].sub.1] = 2 [cos.sup.-1] ([A.sub.+])

[DELTA][[phi].sub.2] = [pi] - 2 [cos.sup.-1] ([A.sub.-]), (6)

where the parameters [A.sub.+,-] are given by

[mathematical expression not reproducible]. (7)

Figure 3 shows the angular width for the P wave radiation as a function of the parameter [bar.[lambda]]; it can be seen that the angular aperture for the radiation field of P waves in general is less than that of S waves, which shows the importance of S waves in the fracture process material by detonating cylindrical explosive charges.

2.2. Particle Velocity. The velocity time histories [v.sub.R](r, z, t) and [v.sub.[phi]](r, z, t) predicted by the Heelan solution for the pressure function equation are

[mathematical expression not reproducible]. (8)

The particle velocity is defined as [absolute value of [??]] = [square root of [v.sup.2.sub.R] + [v.sup.2.sub.[phi]]],

[mathematical expression not reproducible]. (9)

In general, (9) is a complex function of the coordinates and time, but we will see that it maybe susceptible to considerable simplifications given the angular structure discussed above.

Besides, P wave and S wave produced by blasting will separate gradually because of their velocity differences when the distance from measured point to blasting source increases. Another important aspect is that the width of the wave fronts are given approximately by [[delta].sub.p] ~ ([c.sub.p]/[gamma]) and [[delta].sub.s] ~ ([c.sub.s]/[gamma]) for P wave and S wave, respectively. Thus, the Heelan model shows that the P wave and the S wave do not overlap in time.

Figure 4 shows the advancing wave fronts as a function of time and the spatial coordinates. The values of the parameters used are [c.sub.p] = 5000 m/s, [c.sub.s] = 3000 m/s, a = 2l = 0.1 m, and three different values of the parameter [gamma].

It can be seen that as the value of the parameter [gamma] increases, the width of the wave fronts decrease, making more evident uncoupling between the P and S waves. However, the angular distribution of the radiation remains intact following the distribution given by (4). Another important aspect is that the regions near the axis of the explosive charge, vibration induced by the explosive charge, are very low; this feature is also explained by (4), where the angular opening for waves in this region is very small.

2.3. The Peak Particle Velocity. The peak particle velocity is given by

[mathematical expression not reproducible]. (10)

As demonstrated above, the P wave and the S wave do not overlap in time. Thus, the peak particle velocity can be expressed approximately as follows:

[mathematical expression not reproducible]. (11)

Note that the maximum value of peak particle velocity in time depends only on the maximum value of the function ([d.sup.2]/d[t.sup.2]){p(t - R/[c.sub.p,s])}; that is,

[mathematical expression not reproducible], (12)

where [G.sub.n] is a numerical value that depends on n, for n = 4, [G.sub.n] = 1.85.

Finally, taking into account the above considerations, the following approximate expression for the PPV is obtained:

[mathematical expression not reproducible]. (13)

Figure 5 shows the PPV diagrams numerically obtained from (10) (continuous line) compared with the analytic solution (13) (dashed line). You can see that the analytical solution predicts quite well the numerical results, showing that this analytical solution allows a useful tool for simple and rapid assessments under field conditions.

3. Large Cylindrical Charge a/L [much less than] 1

In the case when the length (L) of the cylindrical charge is not small, the angular distribution of the velocity field is very different from the case of a charge of small length. Mainly because when a/L [much less than] 1, the detonation of the explosive cannot be considered instant and it becomes necessary to consider the finite velocity of detonation VOD. A large cylindrical charge can be modeled as a superposition of an array of n Heelan's elements considering a delay time for the detonation of each element (Figure 6). The length of each element is defined as [DELTA]L = L/2n.

3.1. Particle Velocity. The field generated by the superposition of the charges is given by

[mathematical expression not reproducible], (14)

where the field generated by the ith charge located in position [z.sub.i] is given by

[mathematical expression not reproducible] (15)

with

[mathematical expression not reproducible]. (16)

Figure 7 shows the advancing wave fronts as a function of time and the spatial coordinates. The values of the parameters used are [c.sub.p] = 5000 m/s, [c.sub.s] = 3000 m/s, L = 6m, and three different values of the parameter [gamma]. From Figure 7 it can be seen that (14) reproduces well the fronts of P and S waves, where cones formed by these fronts are appreciated. For smaller times (t < L/[c.sub.d]), represented fronts coincide quite well with those found by Rossmanith and Kouzniak [13].

3.2. The Peak Particle Velocity. In this case, the peak particle velocity is given by

[mathematical expression not reproducible]. (17)

Figure 8 shows the PPV diagrams numerically obtained from (17). The values of the parameters used are [c.sub.p] = 5000 m/s and [c.sub.s] = 3000 m/s. To analyze the effect of the charge length L and attenuation coefficient on the PPV diagrams, three different values for these parameters were considered. The results show that, for moderate charge lengths (L ~ 1m) and moderates values of [gamma], the PPV diagrams (Figure 8(a)) are very similar to those delivered by Heelan solution for a charge of small length. As gamma increases (at constant L), the symmetry of the PPV diagrams with respect to the plane z = 0 is broken, showing a higher intensity of vibration in opposite to the initiation of detonation regions (Figures 8(b)-8(c)). The same effect is found if [gamma] is kept fixed and the value of the length of the load increases (Figures 8(e)-8(h)), which shows that there is a relationship between [gamma] and L governing structure diagrams PPV.

3.3. Approximate Analytical Solution. To find an analytical solution of the field of vibration, we can take the following approach in the far field:

[R.sub.i] = R - [z.sub.i] cos [[phi].sub.i]

[[phi].sub.i] = [phi]. (18)

In the case n = 0, the field generated by the ith charge located in position [z.sub.i] is given by

[mathematical expression not reproducible]. (19)

Introducing ((18)-(19)) into (14), it is not difficult to show that the fields of vibration can be approximated by

[mathematical expression not reproducible]. (20)

Finally, the following expressions for the velocity fields are found:

[mathematical expression not reproducible]. (21)

In the limit [DELTA]L [right arrow] 0 and [gamma]L/[c.sub.d] [much less than] 1, the following expression is obtained:

[mathematical expression not reproducible]. (22)

Figure 9 shows the angular distribution of P and S waves in the case of finite velocity of detonation VOD obtained from (22), where dashed line represented Heelan's solution and continuous line represented the case of VOD finite.

The velocity field obtained from (22) is identical to that of Heelan, except for the geometric factor: 1 - ([gamma]L/2[c.sub.p,s])(cos [phi] - [c.sub.p,s]/[c.sub.d]). This geometric factor represents a first correction to Heelan's solution in the case of finite length of the explosive and finite VOD. In the case [gamma]L/[c.sub.p,s] [right arrow] 0, the velocity field is reduced to the solution of Heelan.

4. Conclusions

This paper demonstrates that the seismic radiation in the far field produced by an explosive column of finite length can be modeled by overlapping of asymptotic solutions in the far field of vibration (Heelan's solution); the results agree quite well with exact solutions [12, 13]. The asymptotic solutions found in this work allow reducing the calculation time and facilitating the understanding of the phenomena of vibration in the elastic medium. In addition, it is shown that the geometry and structure of the peak particle velocity (PPV) diagrams depend strongly on the parameters [gamma]L/[c.sub.p,s] and [c.sub.p,s]/[c.sub.d]. The asymptotic model developed in this paper can be extended to more complex problems, such as the interaction between two or more columns of explosives.

https://doi.org/10.1155/2018/9681465

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

[1] P. A. Heelan, "On the theory of head waves," Geophysics, vol. 18, pp. 871-893, 1953.

[2] P. A. Heelan, "Radiation from a cylindrical source of finite length," Geophysics, vol. 18, pp. 685-696, 1953.

[3] D. W. Jordan, "The stress wave from a finite cylindrical explosive source," Journal of Mathematics and Mechanics, vol. 11, pp. 503-551, 1962.

[4] P. Hazebroek, "Elastic waves from a finite line source," Proceedings of the Royal Society A Mathematical, Physical and Engineering Sciences, vol. 294, no. 1436, pp. 38-65, 1966.

[5] A. M. Abo-Zena, "Radiation from a finite cylindrical explosive source," Geophysics, vol. 42, no. 7, pp. 1384-1393, 1977.

[6] J. A. Meredith, Numerical and analytical modeling of downhole seismic sources: the near and far field, Massachusetts Institute of Technology, 1990.

[7] J. A. Meredith, C. H. Cheng, and M. N. Toksoz, Modelling of downhole seismic sources ii: an analysis of the heelan/brekhovskikh results and comparison of point source radiation to radiation from boreholes, Massachusetts Institute of Technology. Earth Resources Laboratory, 1990.

[8] J. A. Meredith, M. N. Toksoz, and C. H. Cheng, "Secondary shear waves from source boreholes," Geophysical Prospecting, vol. 41, no. 3, pp. 287-312, 1993.

[9] K. M. Tubman, Full waveform acoustic logs in radially layered boreholes, Massachusetts Institute of Technology, 1984.

[10] K. M. Tubman, C. H. Cheng, and M. N. Toksoz, "Synthetic full waveform acoustic logs in cased boreholes," Geophysics, vol. 49, no. 7, pp. 1051-1059, 1984.

[11] D. P. Blair and A. Minchinton, "Near Field Blast Vibration Models," in Proceedings of The 8th International Symposium on Rock Fragmentation by Blasting--Fragblast, vol. 8, pp. 152-159, Santiago, Chile, 2006.

[12] D. P. Blair, "A comparison of Heelan and exact solutions for seismic radiation from a short cylindrical charge," Geophysics, vol. 72, no. 2, pp. E33-E41, 2007.

[13] H. P. Rossmanith and N. Kouzniak, "Supersonic detonation in rock: part 2 particle displacements and velocity fields for single and multiple non-delayed and delayed detonating blast holes," International Journal of Blasting and Fragmentation, vol. 8, pp. 95-117, 2004.

[14] D. P. Blair, "Blast vibration dependence on charge length, velocity of detonation and layered media," International Journal of Rock Mechanics and Mining Sciences, vol. 65, pp. 29-39, 2014.

[15] L. Trivino, B. Mohanty, and A. Munjiza, "Investigation of seismic radiation patterns from cylindrical explosive charges by analytical and combined finite-discrete element methods," in Proceedings of The 9th International Symposium On Rock Fragmentation by Blasting, pp. 415-426, Granada, Spain, 2009.

[16] R. Kumar, D. Choudhury, and K. Bhargava, "Determination of blast-induced ground vibration equations for rocks using mechanical and geological properties," Journal of Rock Mechanics and Geotechnical Engineering, vol. 8, no. 3, pp. 341-349, 2016.

[17] M. W. Lee and A. H. Balch, "Theoretical seismic wave radiation from a fluid-filled borehole," Geophysics, vol. 47, no. 9, pp. 1308-1314, 1982.

[18] C. Yi, D. Johansson, U. Nyberg, and A. Beyglou, "Stress wave interaction between two adjacent blast holes," Rock Mechanics and Rock Engineering, vol. 49, no. 5, pp. 1803-1812, 2016.

[19] L. F. Trivino, B. Mohanty, and B. Milkereit, "Seismic waveforms from explosive sources located in boreholes and initiated in different directions," Journal of Applied Geophysics, vol. 87, pp. 81-93, 2012.

[20] S.-H. Chen, J. Wu, and Z.-H. Zhang, "Influence of millisecond time, charge length and detonation velocity on blasting vibration," Journal of Central South University, vol. 22, no. 12, pp. 4787-4796, 2015.

[21] D. E. Grady, M. E. Kipp, and C. S. Smith, "Explosive fracture studies on oil shale," SPE Journal, vol. 20, pp. 349-356, 1980.

[22] D. B. Larson, "Explosive energy coupling in geologic materials," International Journal of Rock Mechanics and Mining Sciences, Sciences and Geomechanics Abstracts, vol. 19, pp. 157-166, 1982.

Marco Rosales-Vera (ID), (1) Ivan Pineyro, (1) and Gaston Pinilla (2)

(1) CODELCOTECH, Los Jardines 279, Santiago, Chile

(2) CODELCO CHILE, Gerencia Corporativa de Innovacion y Tecnologia, Huerfanos 1270, Santiago, Chile

Correspondence should be addressed to Marco Rosales-Vera; mrosa001@codelcotech.cl

Received 27 September 2017; Revised 13 February 2018; Accepted 18 February 2018; Published 29 March 2018

Academic Editor: Emil Manoach

Caption: Figure 1: Schematic of a cylindrical charge.

Caption: Figure 2: Angular distribution of P and S waves.

Caption: Figure 3: Angular width for the P wave radiation as a function of the parameter [bar.[lambda]].

Caption: Figure 4: P and S wave radiation patterns for source in a borehole wave fronts as a function of time and the spatial coordinates. The values of the parameters used are [c.sub.p] = 5000 m/s, [c.sub.s] = 3000 m/s, a = 2l = 0.1 m, and three different values of the parameter [gamma].

Caption: Figure 5: PPV diagrams for different values of [bar.[lambda]] = [c.sub.s]/[c.sub.p].

Caption: Figure 6: Large cylindrical charge modeled as a superposition of an array of elements of Heelan.

Caption: Figure 7: Wave fronts as a function of time and the spatial coordinates. The values of the parameters used are [c.sub.p] = 5000 m/s, [c.sub.s] = 3000 m/s, L = 6 m, and three different values of the parameter [gamma].

Caption: Figure 8: PPV diagrams numerically obtained from (17).

Caption: Figure 9: Angular distribution of P and S waves in the case of finite velocity of detonation VOD obtained from (22) (filled line).

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Rosales-Vera, Marco; Pineyro, Ivan; Pinilla, Gaston |

Publication: | Advances in Acoustics and Vibration |

Article Type: | Report |

Date: | Jan 1, 2018 |

Words: | 3341 |

Previous Article: | A Laser-Based Noncontact Vibration Technique for Health Monitoring of Structural Cables: Background, Success, and New Developments. |

Next Article: | Spectral Representation of Uncertainty in Experimental Vibration Modal Data. |

Topics: |