Fundamental Research of Wave Phenomenon in Frozen Rock Massif during Shaft Sinking

Shaft sinking in frozen soils with adoption of blasting for rocks ripping in some cases leads destruction of freezing columns. Such damages are coincided to rocks interface. At the same time, in case rock stratification is absent in geological section, rapture of freezing columns never has a place. Physical explanation of this phenomenon is possibly on the base of wave processes, which having place in specific conditions of shaft sinking with adoption of blasting works.


Introduction
Review of theoretical studies, concerning to a problem of negative development's investigation of concomitant to the ground freezing process in shafts, and shows, that deformation of freezing columns often coincided with a regions of interbedded layers of frozen rocks with different acoustic stiffness, such as frozen sands and clays.But reason of this phenomenon is not explained [1].This article will show the results of researches which explain appear of thickness resonance inside the vertical nonhomogeneous layer, abutting to later face of built shaft.And it will gives appraisal of feature distribution of normal waves in the waveguide generated between the free surface of pit bottom and a horizontal interface throughout the height for the heterogeneous waveguides.Critical thickness of the waveguide that leads to deformations of the frozen columns has been discovered.

Problem Formulation
Presence of horizontal interface of soils with different acoustic stiffness, lower pit bottom, creates prerequisites for beginnings of waveguide zone with drastic horizontal interfaces.Lower bound of this waveguide is an interface with different type soils, such as our example between frozen water-logged sand and clay.And the upper bound is a rock face, it is absolutely soft wall with reflection factor = 0 on pressure V [3].
Length of waveguide is correspond to the rock face (L = D = 7 m).At the both ends of waveguide, it mates with vertical layer with thickness h v , with dissimilar surface.Altitude of horizontal waveguide h g depending on shaft sinking technology and can range between 2 and 4 m [2].Thus is can be given a ration of waveguide and layer thicknesses: 1 v hg h ≥ .Physical analog of an object in the case being considered in horizontal profile presented in Figure 1.
To make a characteristic of low bound, and consequently a classification of its type, we will estimate reflection index from this bound.
Under variation of longitudinal waves speed in frozen clay over the range 1 2800 3200 m / s and in frozen sand 2 5600 m / s , refraction coefficient of P-waves, is: Consequently, low bound can be taken as a "strong" bound, and this waveguide, by ratio of its compounded mediums acoustic characteristics, is similar to so-called near surface sound channel, with a free upper bound (where upper bound-air) and lower drastic bound , on which n ˂ 1.
During the blasting of explosive charge, on each waveguide's end (solution of the problem is in plane) inside of it, appears composite interference wave field, which theoretically, composed of infinite sum of normal waves with different numbers of l.Each normal wave number l is running lengthwise waveguide axis (lengthwise axis X) with its oven phase velocity C l , which depends on the number of wave l, frequency of the source f, height of waveguide h and medium properties inside the waveguide.
At the same time each normal wave number l is stagnant along the thickness of waveguide h (lengthwise axis Z), i.e. amplitude of each normal wave number l congruence 0. For example, in the simplest case, when reflectance from low absolutely rigid bound V 2 = 1, and reflectance from upper interface V 1 = −1, as in the current task, then wave potential Ψ inside of the waveguide [1] [4] Figure 1.An object of physical analog assumes the presence of horizontal interface.
( ) ( ) where 2 2 1 1 π 2 ; ; Z 0 -position of a point source; Z-position of observation point; ( ) ( ) Each normal wave number l propagates lengthwise the layer with a speed: where C-speed of waves in medium of the waveguide; λ -wavelength, C f λ = .
In the case under consideration, lower bound of waveguide is characterized by the absolute number of reflectance V 1 , which in general can be calculated so: θ -is an angel of incidence on the bound.Similar formula can be given for reflectance from the second bound V 2 .In case if waveguide has two bounds with arbitrary index of reflectance V 1 and V 2 , then more generally normal wave field inside the waveguide can be written so: where b = ikcosϴ; ϴ l -roots of dispersion equation When the upper bound is absolutely soft (V 2 = −1), in a rather big distances from the radiator ( 1 r λ > ), when is useful asymptotic path of Hankel transform, then Formula (2) for normal waves will be written so [2]: where Z 0 , Z-radiator's and detector's location.
cos ; sin 1 It being understood that wave field inside the waveguide springs up by point (spherical) source: ( ) where A-vibration altitude of a source; radiating area, V-particle velocity on a sphere.Hereinafter amplitude of primary field will equals to unit (A 0 = 1) X l is the root of dispersion equation: ; ( ) Each normal wave number l is characterized by her distribution of waves amplitude through the thickness of the waveguide (lengthwise axis Z), her phase velocity C l and absorption factor ( ) ( ) To assess traits of normal waves distribution, which are created by explosion source inside the waveguide between the free surface of pit bottom and horizontal interface, we made calculations of first free normal waves amplitude distribution (l = 1, 2, 3) with waveguide height ( ( ) ) for waveguides with different height h ( : 0, 5; 0, 75;1; 2;3; 4. h Also we took into account that, that coordinate of radiator Z 0 is the center of each layer.Value of characteristics for each waveguide and appropriate normal waves are shown in the Table 1 To assess the possibility of excitation by each acoustic mode of thickness resonance in vertical layer with heterogeneous interface, where calculated ankles:

Model and Simulation Result
As a physical analog of item involved, we choose so-cold "acoustic medium", in which are considered just longitudinal acoustic waves.Task solution is given in the plane same to profile of an object.To assess the possibility of thickness resonance excitation by each acoustic mode in the vertical layer with heterogeneous interface were calculated amplitudes of first three waves and angles of incidence on the heterogeneous layer.Calculations were made in software Mathcad by ranking method.
Results of calculations are shown in Figures 2-7.
From derived distributions A l = f(Z) it is clear, that the biggest amplitude in the waveguide's sectional view will be in all instances of main normal wave with l = 1.Besides, in all instances, first mode amplitude's splitting form in sectional view of a waveguide is approximate to piston valve gear.Value of ankle ϴ l is shown in Table 2.
That is why hereinafter we will analyze behavior just of first two modes (l = 1 and l = 2) under changing of waveguide's height.Values of maximum amplitude for the first two modes of a waveguide, with different height h, are shown in Table 3.It is apparent from Table 3 that then less is the height of waveguide h, then bigger is the maximum amplitude of the main mode.On increasing of a waveguide's height up to h = 3, amplitude A 1 (h) will slight decrease.And just for a waveguide with h = 4 there is a sudden amplitude decrease.Thus, ratio of amplitude A l of a waveguide with h = 2 m to the same amplitude for a waveguide with h = 4 m is more than 5 times.Besides, waveguide with h = 4 m has a different main mode amplitude pattern over a cross section-it changes into a double-humped distribution.

Discussion and Conclusion.
From this article, it may be inferred that the pronounced canalization of the wave energy from the radiator which inside the waveguide can be seen when the varying of distance from the pit bottom to a vertical interface in the range of 0.5 to 3 m and the distance h > 4 m this wave effect is practical insignificant.
From the results of this study I suppose to ensure practical guidance on the betterment of the shaft sinking technology with adoption of ground freezing method.On application of blasting workings during shaft sinking, in case of existence of interface between frozen soils in front of a pit bottom, to avoid possible destruction of freezing columns, we recommend to use blasting just in a distance not less than 4 meters between the pit bottom and horizontal interface [6].
To increase the probability of trouble-free operation of blasting works I suppose that the stowage of explosive should at the depth of 0.5 m from pit bottom face.When shaft sinking in a distance is less than 4 m to the rock interface it is important to use another shaft sinking method in frozen rocks, such as create of a relieve slot by the help of compensation holes.
The value of dispersion equation root for normal waves number l = 1, 2, 3 was calculated on the base of graphic solutions by congruent curve.Calculations were made for waveguides with height : 0,5;0, 75;1; 2;3; 4.

Figure 2 .
Figure 2. First normal wave amplitude distribution according to a height of a waveguide ( 0 0.5 , 0.25 h m Z m = =).

Figure 5 .Figure 6 .
Figure 5. First three normal wave's amplitude distribution according to a height of a waveguide ( ) 0

Table 1 .
. Characteristics of the waveguide and normal waves with different number l.

Table 2 .
Angle ϴ l according to h.

Table 3 .
Dependence of maximum amplitude A l from a waveguide's thickness.