岩石蠕变过程中声发射活动的实验研究

在常温下、单轴加压下对两种岩石完整样品进行蠕变破坏过程中声发射活动实验.根据实验获得的数据讨论了蠕变三个阶段中声发射活动特征.在瞬态蠕变阶段,随着加载,声发射活动逐渐增强;在稳态蠕变阶段,在较低的声发射活动背景上,可能出现较强的信号,声发射幅值-频度关系的 m 值可能稳定于某一低值,也可能下降;进入加速蠕变阶段,声发射活动水平增强,m 值下降加快,也可能表现为明显回升后再下降,并与岩石样品的性质有关.将蠕变破坏三个阶段中声发射活动与地震孕育过程中某些地震活动图象变化(如 b 值、前震活动)进行了联系,提出地壳岩石蠕变破坏机制形成前震——主震——(余震)型地震序列类型的可能性.      

The influence of stiffness of the true triaxial testing machine on the rupturing and acoustic emission of rocks and its correlation with seismic activities

This paper introduces how the ratioR of the characteristic stiffness of rock samplevs. the stiffness of testing machines would influence the rupturing process and the acoustic emission (AE) on the part of the tested rock samples. Result of the experiment shows: WhenR>0.20, the rock sample would rupture abruptly; whereas whenR<0.20, the rock samples would rupture slowly. When the samples rupture abruptly, the time-dependent variation of the AE rate takes such a pattern:peak value—stable low values—rises to the maximum value (concentration)—drops back to the minimum value (quiescence)—(rises again)—ruptures. Moreover, smallerR-value tends to be associated with longer quiescence and vice versa. WhenR>1.50, no pre-failure quiescence is detected. When the rock samples rupture slowly, the variation pattern of the AE rate (after the stress has increased to more than 50% of the rupturing stress) is as the following:stable low (or high) values—rises (or drops) to its maximum (or minimum) values and then continues for some time—ruptures. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 223–233, 1991. This study is supported by the Chinese Joint Seismological Science Foundation.     

The experimental research on b value of AE for the rock specimens with pre-existing crack or notch under uniaxial compression

A lot of researches onb value have been made in seismology. Since the 1960’s Mogi, Scholz and others have studied AE of rock specimens in laboratory and discovered that it is related to natural earthquakes. All former researchers used integral specimens to studyb value in the laboratory. However a major earthquake is usually related to a existing seismic-fault in that area. For this reason, a series of fracture experiments with rock and glass specimens having pre-existing crack or notch is performed in order to examine the effect of preexisting crack tob value. The experimental results show that theb value begins to decrease as soon as the initiation of the crack and finally drop to a very low value when the specimen breaks unstably. Based on these, a brief discussion on the possible mechanism ofb value change for natural earthquakes is given. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,9, 393–400, 1987. Projects sponsored by the Chinese Joint Seismological Science Foundation.     

Lateral density anomaly distribution in the lower mantle obtained from satellite low order geopotential coefficients

In this paper, we attempt to use satellite gravity data and a new inversion method to study the lateral density anomaly distribution in the mantle. First, density difference Δρ(τ,θ,φ) is expanded in terms of a three—dimensional orthogonal function system, the coefficients of the expansion are to be determined. Then, a set of observation equations is established from the relationship between density anomaly and disturbing geopotential. In the equations the unknown vector contains the coefficients of density anomaly expansion, the observational vector is obtained by computing geopotential perturbations using the potential coefficients of GEM10B, and a filtering process is done for the observational values by properly selecting the harmonic degrees of geopotentical. Finally, the lateral density variations in the lower mantle (670 km toCM boundary) are investigated. In this case, the degrees of disturbing geopotential are selected as 2–11, the truncated degrees of density anomaly expansion are taken asL=6 andK=4, and the damping least squares method is used to solve the observation equations. The resulting model shows the high level of density perturbations at 670 km discontinuity and core — mantle boundary, a high — density zone circumscribing the Pacific and a lower—density region under the center of Pacific. These features are in agreement with the three—dimensional seismic velocity variation features by Dziewonski (1984). In the Antarctic region and some parts of Atlantic and Indian Ocean, however, the resulting density anomalies are negatively correlated with the seismic velocity anomalies, the cause resulting in these phenomena is preliminarily analysed in this paper. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 53–65, 1991. The principle and method represented in this paper can also be suitable to study the lateral density anomaly distribution in the earth’s crust and the upper mantle.     

The experimental study of parameters of rocks under ultrahigh-pressure state

The parameters of rocks in ultrahigh-presure state are very important to understanding the geodynamics and explosion engineering of rocks. In these experiments the plane wave dynamite lenses were used to drive the flying plates. The flying plate hits the target and produces the ultrahigh-pressure shock wave. The propagation velocity and the particle velocity of shock wave in rock samples were measured by probes, and the relation between density of rocks and pressure was obtained by Rankine-Hugoniot formulas. The maximum pressure in experiments was higher than 70 GPa. The ultrahigh-pressure state formulas obtained are: Limestoneρ=3.22+0.016p; Graniteρ=3.20+0.024p. In which the unit for ρ is g/cm³, the unit forp is GPa. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 96–103, 1991.     

Time interval of earthquake recurrence for strike - slip fault

On the basis of fault’s dynamic model of Knopoffet al. (1973), this paper has finally obtained a simple approximate formula to be able to estimate the recurrence time intervalT _R of earthquake on strike-slip fault. Preliminary result holds thatμ andδ _s — δ _f have not much effect onT _R . Leta is the ratio of the coseismic displacementD _s to the total displacementD _t in whole event course, i.e.,a =D _s /D _t , thena = 1/3 may represent the standard theoretical state in whichT _R is independent onμ andδ _s — δ _f . At this time,T _R is the arithmetic average ofs ₀/v andkd/β, wheres ₀ is the long-term preseismic accumulated slippage,v is fault’s average displacement rate,d is the fracture length on the fault of seismic focal region andβ is shear wave velocity. In addition,k =υ ₀/aυ, whereυ ₀ is the initial fracture velocity of actual structure at the coseismic instant. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,14, 187–194, 1992. This paper is a part of contracted item of State Seismological Bureau — Tectonic Physical Study of Earthquake Recurrence Period and Characteristic Magnitude.     

Research on the diffusion of pore fluid and in-situ hydraulic diffusivity in the epicentral region of Xinfengjiang reservoir earthquakes

According to the fact that the Xinfengjiang reservoir earthquakes are caused mainly by water seepage, this paper using the data ofM _s⩾2. 0 earthquakes, studies the hydraulic diffusivity of the mainshock zone by the expansion of the distribution area of epicenters. It is indicated thatin-situ hydraulic diffusivity during the preshock activity of the Xinfengjiang reservoir region was about 6. 2 m²/s. However, after the main shock, thein-situ hydraulic diffusivity in the main shock region increased by fifty percent, that is to say, to 9. 7 m²/s. During the long period after main shock occurrence thein-situ hydraulic diffusivity was affected by significant anisotropy of the medium and fluctuation of water level. No regularity can be found. In addition, we compare the diffusivity found by experiment with rock samples collected with thein-situ hydraulic diffusivity estimated. It is shown that the diffusivity of rock samples with fractures is about the same as the diffusivity estimatedin-situ. However, the diffusivity of whole rock samples is 3 orders of magnitudes smaller than that estimatedin-situ. Finally, we discuss the limits on the method by the expansion of distribution area of epicenters in the study of reservoir induced seismicity. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 364–371, 1991. This subject is sponsored by the Chinese Joint Seismological Science Foundation.     

The remote sensing observation in experiments of rock failure and the beginning of remote sensing rock mechanics

The remote sensing observational study for infrared radiation of rocks was proceeded during the loading on rocks until failure. The major instruments used in experiments were transient spectrum apparatus, intelligent spectrum apparatus, infrared radiation thermometer, infrared spectrum radiometer, and infrared thermal imaging system. The experiments for 26 kinds of rocks were made. The studies show that infrared radiation temperature of rocks increases along with increasing of stress. The amplitude of infrared radiation spectrum of rocks also increases along with increasing of stress. The observational results of infrared thermal imaging of rocks are consistent with infrared radiation temperature. Before formation of major faults for some rocks, the belt-shape thermal imaging of temperature anomaly displaies in position of future major faults. This study has led the new technology of remote sensing into rock mechanics and tend to establish a new field in rock mechanics — remote sensing rock mechanics (or remote sensing rock physics). The application of remote sensing rock mechanics in prediction of earthquake and rock burst, and in measurement of stress field in rock mass is expected. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,14, Supp., 645–652, 1992. Jin-Shen HAO, Ji-Han LI, Xiao-Hong LIU, Yi-Qiao ZHI, Jin-Kai ZHANG, Yong-Hong Lü, Yi LIU, Yun-Shen YU, He ZHANG, Quan-Quan JI, Xiao-Fan ZHU and Ning CHEN took part in this work. This subject is supported by the Chinese Joint Seismological Science Foundation (91006). Work of Institute of Geophysics, SSB (93A0009).     

Exploring nonlinear characteristics in seismogenic process

Power spectral characteristics of the seismic activities before and after the seven large earthquake — Haicheng, Tangshan, Lancang, Longling, Mabian, Puer and Songpan are studied by spectrum analysis method. The results show that they all have continuous power spectrum and aperiodic behaviour. When the time away from the large earthquake occurrence, the characteristics of the seismic activities are obviously accompanied by stochastic components. About one year before main shock, the moving orbits constructed by spectrum amplitudes showed the contracting tendency. And the properties during that time may display as flicker noise. One’s another time almost mixed whole phase space. But the essence of seismic activities after main shocks is a chaotic state accompanied with much noises. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,15, 276–281, 1993. This subject is supported by the Chinese Joint Seismological Foundation.     

Regional characteristics of coda attenuation after the lancang-Gengma earthquake

In the light of the single scattering model of coda originating from local earthquakes, and based on the aftershock coda registered respectively at the 4 short period stations installed near the foci shortly after theM7.6 Lancang andM7.2 Gengma earthquakes, this paper has tentatively calculated the rate of amplitude attenuation and theQ _c-value of the coda in the Lancang and Gengma areas using a newly-founded synthetic determination method. Result of the study shows the rate of coda amplitude attenuation demonstrates remarkable regional differences respectively in the southern and northern areas. The southern area presents a faster attenuation (Q _c=114), whereas the northern area shows a slower attenuation (Q _c=231). The paper also discusses the reasons causing such differences. Result of the study also suggests a fairly good linear relation between the coda source factorA _o(f) and the seismic moment and the magnitude. Using the earthquake scaling law, the following formulas can be derived: lgM ₀=lgA ₀(f)+17.6,M _D=0.67lgA ₀(f)+1.21 and logM ₀=1.5M _D+15.79. In addition, the rates of amplitude attenuationβ _s andβ _m are respectively calculated using the single scattering and multiple scattering models, and the ratioβ _s/β_m=1.20−1.50 is found for the results respectively from the two models. Finally, the mean free pathL of the S-wave scattering in the southern and northern areas are determined to be 54 km and 122 km respectively by the relations which can distinguish between the inherentQ _i and scatteringQ _s, testify to this areas having lowQ-values correspond to stronger scatterings. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,14, 71–82, 1992. This study is partly supported by the Seismological Science Foundation of the State Seismological Bureau of China, and the present English version of the paper is translated from its Chinese original by Wenyi Xia, Seismological Bureau of Yunnan Province.     

Effect of occurrence rate of acoustic emissions on their statistical behavior

The parameterm in Ishimoto-Iida's relation was investigated for acoustic emissions (AEs) occurring in rock samples under uniaxial compression. In the experiment, we found: 1) The large AEs are counted without serious error but the number of small AEs is systematically underestimated at high AE rates, 2) the frequency distribution of maximum AE amplitudes becomes nonlinear in logarithmic scale with increasing AE rate, and 3) there exists a strong negative correlation betweenm-value and AE rate. The miscount of small AEs was interpreted as due to overlap of the large and small AEs. We call the miscount masking effect. A statistical analysis based on the masking effect showed that them-value decreases more effectively as the AE rate increases, and thus the masking effect is a possible origin both for the nonlinear frequency distribution of maximum AE amplitudes and for the negative correlation ofm-value with AE rate. We emphasize that one should be careful of the masking effect to examine correctly the change, ofm-value. In order to eliminate the masking effect, AEs should be measured by a measurement system with low sensitivity. Even if the masking effect is eliminated, them-value decreases before the main fracture of a rock sample. Them-value is a key parameter to predict the main fracture.     

The effect of stress on pore pressure in rocks and the mechanism of water table anomaly before earthquakes

Triaxial compressive experiments of porous rock samples were carried out under various confining pressures and initial pore pressures without drainage; axial strain and pore pressure were observed versus differential stress. The results of such experiments show that pore pressure increases with increase of differential stress at low differential stress; pore pressure decreases with increase of differential stress at medium and high differential stress. Pore pressure also increases with large amplitude decrease of differential stress at high differential stress. Based on such experiments, it is suggested that water table anomaly before an earthquake reflects the change of differential stress in crustal rocks. The anomalous behavior of water tables in the epicentral and peripheral areas before the great Tangshan earthquake of July 28, 1976 are explained by such suggestion. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 88–95, 1991. This study is supported by the Chinese Joint Seismological Science Foundation. Professor Yongtai Che gave us much help in applying fund support and supplying earthquake case histories.     

The revised m Lg scale in six eastern provinces of China

The influence of the four different methods of measuring Lg amplitude, and the selection of different geometrical attenuation coefficient ζ-values (=5/6 or 1) on the determination ofγ value of Lg wave are discussed.γ=0.0034±0.0001km⁻¹ (when ζ=5/6) for six eastern provinces is redetermined. The revised magnitude calibration function ofm _Lg (mxh),q _E (Δ)=(5/6)logΔ+0.00147Δ+1.81 is deduced. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 171–178, 1991. Projects sponsored by the Chinese Joint Seismological Science Foundation.     

Spatial and temporal characteristics of acoustic emission activity during deformation of rock samples with inhomogeneous fault under biaxial compression

The spatial and temporal distributions of acoustic emission （AE） during the deformation of samples containing an inhomogeneous fault have been studied under biaxial compression. The results show that the fault strength and the duration from loading to failure increase and the failure mode changes from abrupt instability to gradual failure with increase of lateral stress σ2. The pre-setting fault and its heterogeneity play an important role in controlling AE spatial distribution during the deformation. The basic pattern of AE spatial distribution is controlled by the pre-setting fault, especially the parts with inhomogeneous strength and the strong segment, and the localization of fracturing starts from the positions with inhomogeneous strength, With increase of σ5, the dense AE distribution area spreads from the positions with inhomogeneous strength to the whole strong segment gradually. AE temporal sequence is significantly affected by σ2. The fault shows abrupt instability in final failure, which occurs on a background of ＂enhancement-quiescence＂ in microfracturing activity at lower σ2. At higher σ5, the fault shows gradual failure, and AE occurrence rate increases continuously and AER increase exponentially before and after the failure, The effect of σ2 on b-value is also remarkable, b-value shows precursory decrease in the weakening stage when fault behavior is abrupt instability, but it shows balanced change when fault behavior is gradual failure,     

Source parameters of bohai earthquake, July 18, 1969 and yongshan earthquake, May 11, 1974 determined by synthetic seismograms of teleseismic P waves

The source parameters of the Bohai Sea earthquake, July 18, 1969 and Yongshan, Yunnan earthquake, May 11, 1974 were determined by full — wave theory synthetic seismograms of teleseismic P waves. P+pP+sP wereform were calculated with WKBJ approximation and real integral paths. One — dimensional unilateral, finite propagation source was also considered. By trail — and — error in comparing the theoretical seismograms with the observational ones of WWSSN stations, the source parameters were obtained as follow: for Bohai earthquake, φ=195°, δ=85°, λ=65°,M _o=0.9×10¹⁹Nm,L=59.9km.V _R=3.5km/s, ∧ _R =160°; for Yongshan earthquake, φ=240°, δ=80°, ∧=150°,M _o=1.3×10¹⁸Nm,L=48.8km,V _R=3km/s, ∧ _R =−10°, where φ is strike, δ dip angle, λ slip angle,M _o seismic moment,L rupture length,V _R rupture propagation speed. As III type fractures the faulting propagated along the fault planes, and ∧ _R is the angle from the strike to the propagation direction. Yongshan earthquake showed complexity in its focal process, having four sub—ruptures during the first 60 seconds. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 1–8, 1991.     

Variations of the geomagnetic field and the geomagnetic effects of three earthquakes ( M s>5.0) in Jiangsu Province, China

Since 1979 the repeated observations and experiments of geomagnetic total intensity and vertical component have been carried out for ten years in the geomagetic network which is located in Jiangsu Province, China. Three earthquakes aboveM _s 5.0 occurred during the decade, and some seismomagnetic effects were observed. The observation results show that the anomalies of the vertical geomagnetic component can’t be observed untill some months before the earthquake (M _s>5.0) in this area. In this paper it is suggested that a densely distributed network for continuous observation of geomagnetic vertical component may catch seismomagnetic anomalies and thus improve earthquake prediction in the light of the geomagnetic measurements of the mid — or — low latitude locations. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 80–87, 1991. This study is sponsored by the Chinese Joint Seismological Science Foundation.     

Granite deformation and behavior of acoustic emission sequence under the temperature and pressure condition at different crust depths

IntroductionFocal depths of shallow strong intraplate earthquakes are mostly distributed in the highstrength range of lithosphere controlled by the rheology property of granite and diorite (Brace,Kohlstedt, 1980; Sibson, 1982; Meissner, Strehlau, 1982; CHEN, Molnap 1983). This is obviouslyrelated to the change of rock deformation characteriStics at different crust depths. So, for earthquake stUdy, both of the rock failure types (fractUre or rock flow) and its mechanical instabilityforms (…     

The application of wave field continuation for seismic refraction

Based on wave equation a fundamental formula for plane wave propagaton is derived, the theoretical method of inverting wave volocity structure of earth interior through the observational wave field is discussed in this paper. The observational wave field can be decomposed into the surface plane wave field by using τ −p transform. The maximum amplitude curve in the plane wave field can show stably the change tendency of wave velocity with depth in the Earth. This property can be used to restrict the solution space. By using the wave field continuation method more useful information from the observational wave field can be extracted and the inverion solution not only can be obtained simply and quickly, but also is stable and less influenced by the subjective factor. The wave field continuation is a fine inversion method. Theoretical analysis and numerical modelling are carried out in the study of wave field continuation. By applying homomorphic decovolution the signal — to — noese ratio is improved. Finaly a sonar refraction profile in the northern part of the South China Sea is interpreted and computed. It is found as a result that there is a velocity interface from 1.76 km/s to 2.21 km/s at the depth of 1.4 km. The velocity gradients in the upper and lower layers are 0.54 km · s⁻¹/km and 0.63 km · s⁻¹/km respectively. A discussion of the characteristics of shallow sea structure in the view of tectonic movements is geven. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13. 41–52, 1991.     

The possible changes of tectonic stress fields in North China in 16–17th centuries

The ancient tectonic stress field in North China during 16–17th centuries were recovered by studying seismogenic faults of four great earthquakes. Three of them are dip slip events, the direction of the maximum principal stress is vertical, and that of the minimum principal stress is in the NW—SE direction and nearly horizontal tensile stress. Another earthquake is a complicated one which includes thrust and dip-slip activities. The above-stated character is quite different from that of present stress field, it is inferred that the stress field has significantly changed since 16th century, the main stress axis revolved an angle of 90 degrees nearly. The mantle upheaval in the faulted basins in North China during 16–17th centuries may cause the tensile stress field and a series of great earthquakes. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,14, 144–148, 1992.     

Metamorphic characteristics and geotectonic implications of the high-pressure granulites from Namjagbarwa, eastern Tibet

A large area of high-pressure garnet-kyanite granulite is exhumed in the Namjagbarwa area, which provides a window for observing the deep crust rocks and structures of the Tibetan Plateau. Three mineral assemblages can have been distinguished in the garnet-kyanite HP granulites by petrography, i.e. M₁. Mus+Bi+P1+Q, M₂. Gt+Ky +perphite/antiperphite+Rt+Q, M₃. Gt+Sill+Cord+Sp+Ilm ± Opx. Metamorphic conditions of the peak granulite assemblages (M₂) formatted by thickening of crusts, with available isotopic ages of 45–69 Ma, are at 1.4—1.8 Gpa and 750—850°. Their retrograde assemblages overprinted by decompressure during the uplift, with available isotopic ages of 18—23 Ma, were formed at 0.60—0.70 Gpa, 621—726°. The thermobarometric evaluation, petrogenetic grid and corresponding isotopic ages indicate a clockwise isothermal decompression metamorphic path. The HP granulite metamorphic history indicates that the collision of the Indian Plate with the Eurasian Plate had begun at 70 Ma, far earlier than the widely accepted 45 Ma. Project supported by the National Natural Science Foundation of China (Grant No. 49732100), the National Key Project for Basic Research, and the Chinese Academy of Sciences Project for Tibetan Research Project (GrantNos. KZ951-A1-204, KZ95T-06).