The Effect of Harmonic Oscillations on the Deformation and Acoustic Regime of the Fault Zone Model

G. A. Sobolev and A. V. Ponomarev

Institute of Seismology, UIPE RAS, Moscow, Russia


Electronic version of the paper published in:
Seismology in Europe, Papers presented at the ESC XXV General Assembly, Sep. 9-14, 1996, pp. 94-99.
© Copyright 1996 by G. Sobolev and A. Ponomarev
© Copyright 1997 by Geophysical Center RAS (electronic version only)

Figures 5


Abstract

Experiments were carried out on concrete models using biaxial compression. The models were composed of the outer layers with higher durability and the middle layer with lower durability and including rock debris of different sizes. It was assumed the model simulated the fault zone. The construction of the models and the selected type of load allowed the realization of the classical type of failure by accumulation of subparallel echelon-like cracks in the inner layer with subsequent displacement macrofault. The location of the acoustic emission (AE) events permitted the tracing of the geometry and the stages of fissuring. The experiments were carried out at the average rate of longitudinal deformation of 10-6 s-1. The major compressing load was modulated by sinusoidal vibration with 2, 10, 30, or 100 s periods. This treatment resulted in smoother transition of the model into post-failure stress state and prolonged the stage of deformation from the peak load to macrodestruction. It was observed that in the spectra of the AE time series the number of acoustic events and amount of energy release are concentrated in 20 - 30% interval near the cyclic load maximums. It is shown that the distortion of the form of cyclic changes in the load and deformation is a prognostic symptom of pending instability. It could be provide a possibility to reveal the unstable stage without knowledge of the absolute load and deformation values. It is suggested such approach can be used for estimation of transition of rock massif into post-peak stage under Earth crust conditions.

Go to Contents


Introduction

It is generally recognized, that the Earth's crust is periodically influenced by forces of cosmic and meteorological origin.The rhythmic occurrences of many terrestrial processes have been also established in a wide range of periods [Gamburtsev et al., 1994]. The seismicity of the Earth is within the scope of these processes by its diurnal, seasonal and even longer periodic variations.

The earth tides is a powerful source of periodic stresses and deformations in the Earth's crust. Since the publication by Nishimura [1950], there were many attempts at using earth tides for the assessment of seismic hazards [Seymour, 1972; Mikumo et al., 1978; Burton, 1986; Yaramanci et al., 1988; Hartzell and Heaton, 1989; Nikolaev and Nikolaev, 1993; Xiang-chu Yin et al., 1994; et al.] The effect of mechanical vibrations in addition to the main load, however, may alter the destruction process itself. In their paper, Sadovski et al. [1981] have shown that this influence on the samples of different rocks and man-made materials provokes transition from brittle destruction to the release of elastic energy, accumulated in the sample, by plastic deformation.

A paper by Sobolev et al. [1996] establishes that additional vibration reduces the time period between the successive movements of the stick-slip type on the contact between the blocks of rock.

The present research is aimed at studying the effect of vibration on deformation, destruction, and acoustic regime of complex samples with a lower strength layer as a model in the first approximation of a fault zone.

Go to Contents


Methods

Figure 1 The experiments were carried out on models composed of concrete. Figure 1 shows the general view of the model. The two lateral layers were stronger than the central layer owing to the 400-brand cement that was used to make them, whereas the 50-brand cement was the binding component in the central layer. Besides quartz sand, the central layer also contained rock chips of different size ranging from 2 to 20 mm. The length of the model on its side facet was 300 mm, the width along the perpendicular was 220 mm. The size of the model in its third dimension was 65 mm. The thickness of the lateral and central layers was the same. Owing to the jutting out angles of the chips, the contacts between the lateral and central layers were rough, and the average angle of their dip towards the vertical axis coincided with the main load F and was 25o.

The lateral compressing stress G was 3 tons and was maintained constant throughout the experiment with the accuracy of 1% by roller padding between the side facets of the model and the pistons of the press. The vertical load F, also applied by roller padding, was gradually increased to its peak value and then reduced in the post-peak area of deformation to correspond to the condition of the constant rate of relative deformation \varepsilon = 10-6. The experiments were conducted on a press "Inova" with the servocontrol at the Borok Observatory of the Institute of Physics of the Earth. The load F and the movement of the pistons in the vertical direction D were recorded with sampling rate 0.2 s or 1 s on PC IBM. The conclusions are based on the results of ten experiments, seven of which had F load modulated by additional vibration with periods 2, 10, 30 and 100 s. The amplitude of vibration was about 12% of the maximum load.

The choice of construction of the models and the type of loading have caused formation of fissures and their development in the weaker central layer. A system of echelonlike fissures appeared, the typical case of which is shown in Figure 1. The macrodestruction occurred as a rupture of the shear type cutting across the system of echelonlike fissures and through the middle and along the central layer.

The appearing acoustic signals were recorded on PC IBM by two systems. The first system had six piezoreceivers and was tuned on registration of all signals above noise level. The output parameters of the system, i.e., amplitude and duration of impulses, allowed to estimate their energy. The second system with eight piezoreceivers recorded the waveforms of signals, and was used to locate the strongest of them. The main energybearing frequencies of signals were confined to the range above 30 kHz. In the course of the experiment, usually lasting several hours, the model was periodically sounded in different directions to measure the elastic waves velocities. The initial average velocities of compressional waves in the lateral layers were 4 km s-1 and in the central layer 2 km s-1. As soon as a system of echelonike fissures was formed, the elastic waves velocity in the central layer was reduced to 1.3 - 1.5 km s-1, whereas in the lateral layers it remained practically unchanged.

Go to Contents


Results

Figure 2 Figure 2 represents a typical example of a comparison between the experiments without vibration (a) and those with additional vibration (b). The difference was that in the first case the destruction was more brittle with a macrorupture near the peak load. The additional vibrations allowed the experiment to reach far into the post-peak area with gradual reduction of the model's strength. Cessation of vibration partially restored the strength at the conclusive stage of the experiment (Figure 2b).

The location of AE events indicates that fissuring always started in the lower part of the central layer, near the mobile piston of the press, and migrated upwards. The arrows in Figure 2b show the moments of three acoustic signals corresponding to the location of their sources 1, 2, 3 in Figure 1.

The acoustic signals with large amplitudes (Figure 2) were recorded only near the peak load and in the post-peak area of deformation. The distribution of the weaker signals during vibration period was also studied. For this purpose their numbers in the subsequent cycles of vibration-relative additional load release - were summed up to obtain "seasonal variations". The spectral analysis of the distribution in time of the number of acoustic signals distinctly reveals the periodicity corresponding to the vibration period. Figure 3 The amount of signals, however, does not strictly follow the sinusoidal law. Most of them are grouped in the 20-30% interval before and near F maximum. Figure 3 represents the mode of seasonal variations as compared with the sinusoidal period of F variation.


One of the principal problems of prognostic importance is to determine the point of maximum rheological curve F - D. The experience of destruction mechanics implies that, on reaching this point, the loaded body passes to the stage of unstable deformation that initiates macrodestruction without outer energy supply. In his paper, Sobolev [1993] emphasized that the discovery of this point in the focus of a pending earthquake should improve prediction methods, because at that stage the short-term precursors become apparent. The difficulty is that we do not know the absolute values of stresses and deformations under conditions in the Earth's interior, though we easily obtain them in the laboratory. Therefore, an idea occurs to use cyclic variations of stresses and deformations to find the point of F - D maximum.

Figure 4 In the present paper, we used the following procedures to study this problem. We obtained a nonlinear low-frequency trend of F and D plots by averaging data from ten vibration periods. The residual plots were then standardized to standard deviation. The next step was to construct the F - D hysteresis from individual cycles with vibration period T. Figure 4 shows the loops of hysteresis averaged from ten cycles on the linear part of load (a) and near the maximum F- D points (b). It is apparent that hysteresis in Figure 4b has a peculiar feature in the part covering about 1/4 of the duration of the cycle at large F and D values. This feature shows that the reduction of deformation lags against the reduction of load on the descending arm of the F - D curve. This is an indication that the deformation process is unstable. Proceeding from this result, an attempt was made to determine the instability stages without the absolute F and D values (Figure 5). Figure 5 The values of \DeltaF and \DeltaD cyclic changes, relative to the lines of the nonlinear trend, were calculated. After standardization the \DeltaF and \DeltaD values, we calculated the functions

FD(j) = \Sigmani =1[(\DeltaFi+n - \DeltaDi+n) - (\DeltaFi-n - \DeltaDi-n)],

where \DeltaFi, \DeltaDi) are load and displacement values at the points of maximums of vibration cycles, while n changes from 1 to n and corresponds to \DeltaF and \DeltaD values within 1/4 of the duration of the cycle, j corresponds to the amount of vibration cycles.

The FD(j) function demonstrates the deviation from the synchronous changes of the F and D curves on the 1/4 portions of the cycle, where the largest differences from the linear law were observed in the hysteresis loops. A correlation of the F loading plot with the F-D function plot in Figure 5 shows that in transition of the model to the instability stage (the peak load on F plot) the FD values increase. A prognostic aspect is consequently revealed. A tendency to deviation from the linear law on approach to instability is also apparent in the study of asymmetry in the F and D cycles separately. The evaluation of stability of this prognostic feature, however, requires supplementary experiments.

There is an apparent opportunity for studying the instability of crustal blocks in the course of an analysis of distortions in the earth tides or in the signals from vibrosources close to peak loads.

Go to Contents


Conclusions

A study has been conducted of the process of deformation and destruction of a model with a weaker layer simulating the fault zone under conditions of biaxial load with additional vibration.

The application of vibration reduces the possible brittle destruction and promotes a smooth transition of the model into the post-peak state.

A distinct periodicity in the appearance of acoustic signals corresponding to the vibration period reveals the appearance of impulses primarily in the 20-30% interval near the cyclic load maximums.

The approach of the model to instability is manifested by the expansion of the hysteresis area in the load-displacement loop.

The distortion of the form of cyclic changes in the load and in the deformation is a prognostic symptom of pending instability that provides a possibility to reveal this stage without knowledge of the absolute force and deformation values.

There is an obvious opportunity to observe an instability in the generating earthquake focus by studying similar changes in the oscillations from earth tides or vibrosources.

Go to Contents


Acknowledgments

We extend our thanks for the help rendered during the research to V. Babichev, A. Koltzov, A. Mostrioukov, A. Patonin - our colleagues at the United Institute of Physics of the Earth, and to D. Frolov, V. Pikulin, S. Stanchits at the Physical-Technical Institute. The work was supported by the Russian Foundation for Basic Research, grant 95-05-14295.
Go to Contents

References

1. Gamburtsev A. G. et al. Atlas of natural processes. Order and chaos in lithosphere and other spheres. Vol. 1. UIPE RAS, Moscow, 1994, p. 176.

2. Nishimura E. On Earth tides. Trans. Am. Geophys. Union, 31, 1950, pp. 357-376.

3. Seymour S. Earthquake-tide correlation. Geophys. J. Roy. Astr. Soc., v.28, 1972, pp.27-84.

4. Mikumo T., M. Rato, H. Doi et al., Possibility of temporal variations in Earth tidal strain amplitudes associated with major earthquakes. Earthquake Precursors. Japan Scientific Societies Press, 1978, pp. 123-136.

5. Burton P. W. Is there coherence between Earth tides and earthquakes? Nature, v.321, 1986 ( 6066, p.115.

6. Yaramanci U., J. Zschau, M. Westerhaus, Tidal tilt admittance at Mudurnu Valley, Turkey. Anhang zum DFG Gesamtbericht uber das Vorhaben Geodatisch-Geophysikal. Untersuchungen im seismoaktiven Bereich der nord-anatolischen Storzone. 1988, 1, 26.

7. Hartzell S., T. Heaton, The fortnightly tide and tidal triggering of earthquakes. Bull. Seismol. Soc. Amer., v.79, 1989, pp.1282-1286.

8. Nikolaev A. V. and V. A. Nikolaev. Earth tides triggering of earthquakes. Continental Earthquakes. Seismological Press, Beijing, 1993, pp. 319-327.

9. Xiang-chu Yin, Xue-zhong Chen, Zhi-ping Song, Can Yin. The load-unload response ratio (LURR) theory and its application to earthquake prediction. Journal of Earthquake Prediction Research, v.3, 1994, p.325-333.

10. Sadovski M. A., K. M. Mirzoev, S. H. Negmatullaev, N. G. Salomov, Effect of mechanical vibration on the characteristics of plastic deformation of materials. Fizika Zemli, 1981, No.6, pp.32-42 ( in Russian).

11. Sobolev G. A., A. V. Ponomarev, A. V. Koltzov and V. B. Smirnov. Simulation of triggering earthquakes in the laboratory. PAGEOPH, v.147, 1996, No. 2, pp1230-1241.

12. Sobolev G. A. Fundamentals of earthquake prediction. Nauka, Moscow, 1993, p.313 (in Russian).


Go to Contents

Load: [HTML document with graphics for local use]  [MS Word file prepared by authors]

HTML version prepared and loaded by V. Nechitailenko on January 23, 1997


Go to Content of GPO, No. 1