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. These results are extended to the case of a charge with finite velocity of detonation.
1. Introduction
The 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. Triviño 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]. Triviño 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 Model2.1. Short Cylindrical Charge (Heelan’s Solution)
Consider a cylindrical explosive charge length 2l 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.
Schematic of a cylindrical charge.
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(1)uRR,ϕ,t=la22μcpR1-2cscp2cos2ϕddtpt-RcpuϕR,ϕ,t=la22μcsRsin2ϕddtpt-Rcs.To simplify the notation we will define the coefficients F1(ϕ) and F2(ϕ):(2)F1ϕ=la22μcp1-2cscp2cos2ϕF2ϕ=la22μcssin2ϕ,where 2l is the total length of the loaded section of the cylinder, ϕ is the polar angle measured from the cylinder’s axis, R is the distance to the observer, cs=μ/ρ and cp=(λ+2μ)/ρ are the shear and compressional wave speeds in the medium, and μ and ρ 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(3)Pt=PVNeγnntne-γtHt,where “e” is the Euler’s constant, PVN 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 γ 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 cs and cp, respectively, where cp>cs. 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.
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 ϕn=(2n+1)(π/4), with n=0,1,2,3. The P waves predominate around ϕn=n(π/2). The region where the P waves predominate over S waves can be estimated from the relationship F1(ϕ)≥F2(ϕ); then(4)1-2cscp2cos2ϕ∗≥cpcssin2ϕ∗.This equation leads to the following expression:(5)cos2ϕ∗=1211+λ-61+λ-4±1-λ-2+2λ-4,where λ-=cs/cp<1. The angular widths for windows radiation P waves are given by(6)Δϕ1=2cos-1A+Δϕ2=π-2cos-1A-,where the parameters A+,- are given by (7)A+,-=1211+λ-61+λ-4±1-λ-2+2λ-4.Figure 3 shows the angular width for the P wave radiation as a function of the parameter λ-; 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.
Angular width for the P wave radiation as a function of the parameter λ-.
2.2. Particle Velocity
The velocity time histories vR(r,z,t) and vϕ(r,z,t) predicted by the Heelan solution for the pressure function equation are(8)vRR,ϕ,t=la22μcpR1-2cscp2cos2ϕd2dt2pt-RcpvϕR,ϕ,t=la22μcsRsin2ϕd2dt2pt-Rcs.The particle velocity is defined as u→=vR2+vϕ2, (9)u→=F1ϕRd2dt2pt-Rcp2+F2ϕRd2dt2pt-Rcs2.In general, (9) is a complex function of the coordinates and time, but we will see that it may be 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 δp~(cp/γ) and δs~(cs/γ) 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 cp=5000 m/s, cs=3000 m/s, a=2l=0.1 m, and three different values of the parameter γ.
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 cp=5000 m/s, cs=3000 m/s, a=2l=0.1 m, and three different values of the parameter γ.
It can be seen that as the value of the parameter γ 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(10)u→PP=maxF1ϕRd2dt2pt-Rcp2+F2ϕRd2dt2pt-Rcs2.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:(11)u→PP≈maxF1ϕRd2dt2pt-Rcp,F2ϕRd2dt2pt-Rcs.Note that the maximum value of peak particle velocity in time depends only on the maximum value of the function (d2/dt2){p(t-R/cp,s)}; that is,(12)maxd2dt2pt-Rcp,s=γ2aPVNl2μcp,sennGn=αn,where Gn is a numerical value that depends on n, for n=4,Gn=1.85.
Finally, taking into account the above considerations, the following approximate expression for the PPV is obtained:(13)u→PP≈αnF1ϕR0≤ϕ≤12Δϕ112π-12Δϕ2≤ϕ≤12π-12Δϕ2π-12Δϕ1≤ϕ≤π+12Δϕ1αnF2ϕR12Δϕ1≤ϕ≤12π-12Δϕ212π+12Δϕ2≤ϕ≤π-12Δϕ1.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.
PPV diagrams for different values of λ-=cs/cp.
3. Large Cylindrical Charge a/L≪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≪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 ΔL=L/2n.
Large cylindrical charge modeled as a superposition of an array of elements of Heelan.
3.1. Particle Velocity
The field generated by the superposition of the charges is given by(14)uRr,z,t=∑i=1nuRir,z,tuϕr,z,t=∑i=1nuϕir,z,t,where the field generated by the ith charge located in position zi is given by(15)uRir,z,t=F1ϕiRid2dt2pt-Ricp-zicduϕir,z,t=F2ϕiRid2dt2pt-Rics-zicdwith (16)cosϕi=-z-ziRisinϕi=rRiRi=r2+z-zi2F1ϕi=a2ΔL2μcp1-2cscp2cos2ϕiF2ϕi=a2ΔL2μcssin2ϕi.Figure 7 shows the advancing wave fronts as a function of time and the spatial coordinates. The values of the parameters used are cp=5000 m/s, cs=3000 m/s, L=6 m, and three different values of the parameter γ. 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/cd), represented fronts coincide quite well with those found by Rossmanith and Kouzniak [13].
Wave fronts as a function of time and the spatial coordinates. The values of the parameters used are cp=5000 m/s, cs=3000 m/s, L=6 m, and three different values of the parameter γ.
3.2. The Peak Particle Velocity
In this case, the peak particle velocity is given by(17)u→PP=max∑i=1nuRir,z,t2+∑i=1nuϕir,z,t2.Figure 8 shows the PPV diagrams numerically obtained from (17). The values of the parameters used are cp=5000 m/s and cs=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~1 m) and moderates values of γ, 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 γ 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 γ and L governing structure diagrams PPV.
PPV diagrams numerically obtained from (17).
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:(18)Ri=R-zicosϕiϕi=ϕ.In the case n=0, the field generated by the ith charge located in position zi is given by(19)pt-Ricp-zicd=p0e-γt-Ri/cp-zi/cdHt-Ricp-zicd.Introducing ((18)-(19)) into (14), it is not difficult to show that the fields of vibration can be approximated by(20)∑i=1nuRir,z,t=p0γ2F1ϕRe-γt-R/cp+zncosϕ/cp-1/cd∑i=0neγiΔLcosϕ/cp-1/cd∑i=1nuϕir,z,t=p0γ2F2ϕRe-γt-R/cs+zncosϕ/cs-1/cd∑i=0neγiΔLcosϕ/cs-1/cd.Finally, the following expressions for the velocity fields are found:(21)uRr,z,t~p0γ2F1ϕRe-γt-R/cp+Lcosϕ/cp-1/cd1-eγLcosϕ/cp-1/cd1-eγΔLcosϕ/cp-1/cduϕr,z,t~p0γ2F2ϕRe-γt-R/cs+Lcosϕ/cs-1/cd1-eγLcosϕ/cs-1/cd1-eγΔLcosϕ/cs-1/cd.In the limit ΔL→0 and γL/cd≪1, the following expression is obtained:(22)uRr,z,t=La22μcpR1-2cscp2cos2ϕ1-γL2cpcosϕ-cpcdd2dt2pt-Rcpuϕr,z,t=La22μcsRsin2ϕ1-γL2cscosϕ-cscdd2dt2pt-Rcs.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.
Angular distribution of P and S waves in the case of finite velocity of detonation VOD obtained from (22) (filled line).
The velocity field obtained from (22) is identical to that of Heelan, except for the geometric factor: 1-(γL/2cp,s)(cosϕ-cp,s/cd). 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 γL/cp,s→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 γL/cp,s and cp,s/cd. 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.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
HeelanP. A.On the theory of head waves195318871893MR0059768HeelanP. A.Radiation from a cylindrical source of finite length19531868569610.1190/1.1437923JordanD. W.The stress wave from a finite cylindrical explosive source196211503551MR0148297HazebroekP.Elastic waves from a finite line source19662941436386510.1098/rspa.1966.0193Zbl0151.39801Abo-ZenaA. M.Radiation from a finite cylindrical explosive source19774271384139310.1190/1.14407992-s2.0-0017633523MeredithJ. A.1990Massachusetts Institute of TechnologyMeredithJ. A.ChengC. H.ToksozM. N.1990Massachusetts Institute of Technology. Earth Resources LaboratoryMeredithJ. A.ToksözM. N.ChengC. H.Secondary shear waves from source boreholes199341328731210.1111/j.1365-2478.1993.tb00571.x2-s2.0-0027333649TubmanK. M.1984Massachusetts Institute of TechnologyTubmanK. M.ChengC. H.ToksözM. N.Synthetic full waveform acoustic logs in cased boreholes19844971051105910.1190/1.1441720BlairD. P.MinchintonA.Near Field Blast Vibration Models8Proceedings of The 8th International Symposium on Rock Fragmentation by Blasting—Fragblast2006Santiago, Chile152159BlairD. P.A comparison of Heelan and exact solutions for seismic radiation from a short cylindrical charge2007722E33E41RossmanithH. P.KouzniakN.Supersonic detonation in rock: part 2 particle displacements and velocity fields for single and multiple non-delayed and delayed detonating blast holes2004895117BlairD. P.Blast vibration dependence on charge length, velocity of detonation and layered media20146529392-s2.0-8489038827410.1016/j.ijrmms.2013.11.007TriviñoL.MohantyB.MunjizaA.Investigation of seismic radiation patterns from cylindrical explosive charges by analytical and combined finite-discrete element methodsProceedings of The 9th International Symposium On Rock Fragmentation by Blasting2009Granada, Spain415426KumarR.ChoudhuryD.BhargavaK.Determination of blast-induced ground vibration equations for rocks using mechanical and geological properties2016833413492-s2.0-8500621578410.1016/j.jrmge.2015.10.009LeeM. W.BalchA. H.Theoretical seismic wave radiation from a fluid-filled borehole19824791308131410.1190/1.14413912-s2.0-0020330410YiC.JohanssonD.NybergU.BeyglouA.Stress wave interaction between two adjacent blast holes20164951803181210.1007/s00603-015-0876-x2-s2.0-84945156617TriviñoL. F.MohantyB.MilkereitB.Seismic waveforms from explosive sources located in boreholes and initiated in different directions20128781932-s2.0-8486723516610.1016/j.jappgeo.2012.09.004ChenS.-H.WuJ.ZhangZ.-H.Influence of millisecond time, charge length and detonation velocity on blasting vibration20152212478747962-s2.0-8494995503310.1007/s11771-015-3030-8GradyD. E.KippM. E.SmithC. S.Explosive fracture studies on oil shale198020349356LarsonD. B.Explosive energy coupling in geologic materials198219157166