真三轴压机刚度对岩石破裂声发射的影响及与地震活动性的联系.

本文叙述了岩石特征刚度与压机刚度比R对岩石破裂及声发射的影响.结果表明:当R>0.20时,岩样产生突发断裂;当R1.50时,断裂前没有平静期.缓裂时,声发射率的变化形式为(应力增至大于50％破裂应力后)平稳低(高)值——上升(下降)至极大(极小)值并持续之——破裂.      

The experimental study of acoustic emission during creep of rocks

Under normal temperature, the creep experiments with complete samples of Gabbro and Marble rocks are made under uniaxial compression. It is found that at the instantaneous creeping stage, AE activities increase with loading; at the steady creeping stage, large AE signals may appear at lower background of AE action, andm—value which shows the relationship between AE amplitude and frequency keeps stable on the low value or decreases; at the accelerate creeping stage, AE activities increase andm—value decreases quickly or decreases again after recovery. These experimental results are related to the quality of the samples. In this paper, AE activities during three stages of creeping process are connected with the seismic patterns (for exampleb—value, foreshocks). Finally, it is pointed out the possibility that the foreshock—mainshock—aftershock earthquake sequence has been formed by the mechanism of creep fracture of crustal rocks. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 104–112, 1991. This subject is sponsored by the Chinese Joint Seismological Science Foundation.     

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.     

The Q value variations in the preparing process of rock rupture

In this paper, the "spectral amplitude ratio method" (SAR), "energy method" (EN) and "coda wave method" (CW) are used to calculate theQ value variations of gneiss in the preparing rupture process. The obtained results show that the variation state ofQ values by SAR features the shape of relative stability—gradual increment to the maximum—then decrement and final rupture. The variation state ofQ values by EN is just contrary to that by SAR, i. e. with the shape of stability—decrement—increment—and final rupture. The varation state ofQ values by CW is similar to that by EN, its main frequency features the shape of relatively high value—decrement to the minimum—increment—and final rupture. But to the high frequency (higher than the main frequency), the variation state ofQ values features the shape of the stable value-increment to the maximum-decrement-and final rupture. At the same time, the results by coda wave amplitude spectrum show that, when stress reaches 70% of rupture stress, the high frequency component of S wave rapidly reduces (Q _c increasing); at the time of impending the main rupture, the main frequency component reduces with a large scale (Q _c increasing again), this may be the reason which causes the different variation states of two codaQ values. The result of amplitude spectra of P, S (initial wave) waves also show that with the appearance of microcracks the frequency band of S wave turn to be narrow, the high frequency component is reduced quickly, i. e. the S wave spectra have different variation states with different frequency components. That is why theQ _s obtained by different methods have different variation characteristics.     

2018年印尼帕卢MW7.5地震——一次超剪切破裂事件

利用反投影方法,使用日本密集台网Hi-net远场垂直分量568条P波资料对2018年9月28日印尼帕卢M_W7.5地震震源破裂过程进行成像,结果显示此次地震的能量释放比较集中,主要集中在10~20 s之间.破裂有两个集中区,破裂峰值分别位于12 s和19 s,最大能量释放区域位于震中南侧约0~50 km内,另一破裂集中区覆盖了帕卢市及周边区域.破裂主要向南侧延展,破裂总长度至少100 km,平均破裂速度约4.1 km·s^-1,属于一次超剪切破裂事件.     

Average slip rate,earthquake rupturing segmentation and recurrence behavior on the Litang fault zone,western Sichuan Province,China

The Litang fault zone (LFZ) is an important active fault within the northwestern Sichuan sub-block. To-gether with the Garzê-Yushu, Xianshuihe, and An-ninghe fault zones on its northern, eastern and south-eastern sides, the LFZ constitutes the lateral extrusion tectonic system in the southeastern part of the Qing-hai-Tibetan Plateau[1,2] (Fig. 1). According to instru-mental records, historical recordings and field investi- gation, an earthquake (Ms7.3) occurred on its middle to south se…     

Behaviour of deep immersed tunnel under combined normal fault rupture deformation and subsequent seismic shaking

Immersed tunnels are particularly sensitive to tensile and compressive deformations such as those imposed by a normal seismogenic fault rupturing underneath, and those generated by the dynamic response due to seismic waves. The paper investigates the response of a future 70 m deep immersed tunnel to the consecutive action of a major normal fault rupturing in an earthquake occurring in the basement rock underneath the tunnel, and a subsequent strong excitation from a different large-magnitude seismic event that may occur years later. Non-linear finite elements model the quasi-static fault rupture propagation through the thick soil deposit overlying the bedrock and the ensuing interaction of the rupture with the immersed tunnel. It is shown that despite imposed bedrock offset of 2 m, net tension or excessive compression between tunnel segments could be avoided with a suitable design of the joint gaskets. Then, the already deformed (“injured”) structure is subjected to strong asynchronous seismic shaking. The thick-walled tunnel is modelled as a 3-D massive flexural beam connected to the soil through properly-calibrated nonlinear interaction springs and dashpots, the supports of which are subjected to the free-field acceleration time histories. The latter, obtained with 1-D wave propagation analysis, are then modified to account for wave passage effects. The joints between tunnel segments are modeled with special non-linear hyper-elastic elements, properly accounting for their 7-bar longitudinal hydrostatic pre-stressing. Sliding is captured with special gap elements. The effect of segment length and joint properties is explored parametrically. A fascinating conclusion emerges in all analysed cases for the joints between segments that were differentially deformed after the quasi-static fault rupture: upon subsequent very strong seismic shaking, overstressed joints de-compress and understressed joints re-compress—a “healing” process that leads to a more uniform deformation profile along the tunnel. This is particularly beneficial for the precariously de-compressed joint gaskets. Hence, the safety of the immersed tunnel improves with “subsequent” strong seismic shaking!     

Rupture behavior and deformation localization of the Kunlunshan earthquake (<Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript>7.8) and their tectonic implications

Earthquake surface rupture is the result of transformation from crustal elastic strain accumulation to permanent tectonic deformation. The surface rupture zone produced by the 2001 Kunlunshan earthquake (M _w 7.8) on the Kusaihu segment of the Kunlun fault extends over 426 km. It consists of three relatively independent surface rupture sections: the western strike-slip section, the middle transtensional section and the eastern strike-slip section. Hence this implies that the Kunlunshan earthquake is composed of three earthquake rupturing events, i.e. the M _w =6.8, M _w =6.2 and M _w ⩽=7.8 events, respectively. The M _w =7.8 earthquake, along the eastern section, is the main shock of the Kunlunshan earthquake, further decomposed into four rupturing subevents. Field measurements indicate that the width of a single surface break on different sections ranges from several meters to 15 m, with a maximum value of less than 30 m. The width of the surface rupture zone that consists of en echelon breaks depends on its geometric structures, especially the stepover width of the secondary surface rupture zones in en echelon, displaying a basic feature of deformation localization. Consistency between the Quaternary geologic slip rate, the GPS-monitored strain rate and the localization of the surface ruptures of the 2001 Kunlunshan earthquake may indicate that the tectonic deformation between the Bayan Har block and Qilian-Qaidam block in the northern Tibetan Plateau is characterized by strike-slip faulting along the limited width of the Kunlun fault, while the blocks themselves on both sides of the Kunlun fault are characterized by block motion. The localization of earthquake surface rupture zone is of great significance to determine the width of the fault-surface-rupture hazard zone, along which direct destruction will be caused by co-seismic surface rupturing along a strike-slip fault, that should be considered before the major engineering project, residental buildings and life line construction. Supported by the National Natural Science Foundation of China (Grant No. 40474037) and the National Basic Research Program of China (Grant No. 2004CB418401)     

On the seismicity of central Kamchatka before the 1997, M = 7.9 Kronotskii earthquake

We report results from a detailed study of seismicity in central Kamchatka for the period from 1960 to 1997 using a modified traditional approach. The basic elements of this approach include (a) segmentation of the seismic region concerned (the Kronotskii and Shipunskii geoblocks, the continental slope and offshore blocks), (b) studying the variation in the rate of M = 4.5–7.0 earthquakes and in the amount of seismic energy release over time, (c) studying the seismicity variations, (d) separate estimates of earthquake recurrence for depths of 0–50 and 50–100 km. As a result, besides corroborating the fact that a quiescence occurred before the December 5, 1997, M = 7.9 Kronotskii earthquake, we also found a relationship between the start of the quiescence and the position of the seismic zone with respect to the rupture initiation. The earliest date of the quiescence (decreasing seismicity rate and seismic energy release) was due to the M = 4.5–7.0 earthquakes at depths of 0–100 km in the Kronotskii geoblock (8–9 years prior to the earthquake). The intermediate start of the quiescence was due to distant seismic zones of the Shipunskii geoblock and the circular zone using the RTL method, combining the Shipunskii and Kronotskii geoblocks (6 years). Based on the low magnitude seismicity (M≥2.6) at depths of 0–70 km in the southwestern part of the epicentral zone (50–100 km from the mainshock epicenter), the quiescence was inferred to have occurred a little over 3 years (40 months) before the mainshock time and a little over 2 years (25 months) in the immediate vicinity of the epicenter (0–50 km). These results enable a more reliable identification of other types of geophysical precursors during seismic quiescences before disastrous earthquakes.     

Frequency analysis of the 2004 Sumatra-Andaman earthquake using spectral seismograms

Frequency analysis of the Sumatra-Andaman earthquake of 2004, one of the most significant and best-recorded earthquakes, is based on spectral seismograms obtained from their broadband seismograms. The Sumatra-Andaman earthquake is found to have a wide-range frequency content of P-wave radiation during the rupturing process. On the basis of stacking spectral seismograms we distinguished four time events of the rupturing process of a total length of about 540 s. The frequency, f _max, is the highest for the first event (0.163 Hz in time interval 0–88 s), lowest for the second — which is the strongest (0.075 Hz in time interval 88–204 s). For third and fourth events frequencies are similar (0.089 and 0.082 Hz in time intervals 204–452 and 452–537 s, respectively). The frequency also shows an azimuthal dependence (±0.02 Hz). Azimuths for which the frequency, f _max, has maximum and minimum values are 203–222° and 23–42°, respectively. These observations are discussed in relation to previously published papers on this topic.     

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 (…     

A data collection and analysis system for research on rock fracturing

We have made a data collection and analysis system for recording microcracks in rock to study the three dimensional pattern of fracture and the dynamic properties of rock during pressurization. This system, with 8 channel data detection, records in real-time, from the initial microcracks to the final vibration waveform excited by the master fracture. The detected data are recorded digitally at the hard disk of an IBM computer (or 286, 386 personal computer) as well as transmitted to disk. The sampling rate for each channel is 10 MHz. Collecting every acoustic emission events in the loaded rock sample with 8 channels by this system, mathematically modelling the AE event travel time, and taking the coordinates of the AE event hypocenter and the three components of P wave velocity as unknown parameters, we set up a set of residual equations for joint inversion, so that the three dimensional localization of AE event hypocenters can be completed under variable velocities, which will lay foundations for the research on rock fracture clearly. The Chinese version of this paper appeared in the Chinese edition ofActa Seismological Sinica,13, 489–495, 1991. The design and development of the system, and the experiments are sponsored by the Chinese Joint Seismological Science Foundation.     

Tunnel effect of fractal fault and transient S-wave velocity rupture (TSVR) of in-plane shear fault

ｎｔｒｏｄｕｃｔｉｏｎＴｒａｎｓｉｅｎｔＳｗａｖｅｖｅｌｏｃｉｔｙｒｕｐｔｕｒｅ（ＴＳＶＲ）ｍｅａｎｓｔｈｅｖｅｌｏｃｉｔｙｖｏｆｆａｕｌｔｒｕｐｔｕｒｅｉｓｂｅｔｗｅｅｎＳｗａｖｅｖｅｌｏｃｉｔｙβａｎｄＰｗａｖｅｖｅｌｏｃｉｔｙα．Ｉｔｓｅ...     

Statistical features of aftershock distribution size for moderate and large earthquakes in Chinese mainland

Based on data of earthquake sequences with MS≥5.0 in Chinese mainland from 1970 to 2004,for different se-quence types and different rupture modes of the main shock,the relationship between aftershock distribution size Rand the magnitude of the main shock M0 has been studied statistically.Considering the rupture mode of the mainshock,we give the quantitative statistical relationships between R and M0 under 95%confidence level for differentsequence types.Qualitatively,lgR,the logarithm of the aftershock distribution size,is positively correlative to theM0,but the data distribution is dispersed.Viewing from different sequence types,the correlation between R and M0is very weak for isolated earthquake type(IET)sequence,R distributes in the range from 5 to 60 km;For main-shock-aftershock type(MAT),lgR is positively correlative to M0;For multiple main shock type(MMT),the core-lation between lgR and M0 is not very obvious when M0≤6.2 and R distributes in the range from 5 to 70 km,whileit shows a linear correlation when M0≥6.3.The statistical results also show that the occupational ratios of differentsequence types for strike-slip and oblique slip are almost the same.But for dip-slip(mostly are thrust mechanisms),the ratio of MAT is higher than that of IET and MMT.Comparing with previous results,it indicates that,when M0is large enough,R is mainly determined by M0 and there is almost no relationship with the rupture mode of themain shock.     

Influence of “stress reversal” on rock resistivity during loading procedure

The properties of rock resitivity were studied under pressure, particularly with “stress reversal”, a procedure in which the pressure applied was increased and decreased. It was observed that, 1) With pressure increasing, the main feature of resistivity change was increase-steady-decrease for high-saturation rock samples (saturation 70–100%). But the main feature for low-saturation samples was different. 2) In 10 out of 11 cases of “stress reversal” for high-saturation samples the resistivity droped (about 2%). Such drop could explain the anomalies in geoelectricity terms, which are commonly observed before earthquakes in China. 3) It was also observed shortly before rock failure that, a) the resistivity drops more dramatically (about 20%) during “stress reversal” period, which is much more than ordinary drops. b) these drops occurred not only during stress decrease but also during stress increase. c) Resistivity exhibits anisotropy: the resistivity along different directions may differ by 10%. These three features may indicate that the rock is nearing failure, while ordinary resistivity drops are only connected with “stress reversal” and may not mean the imminence of rock failure. 4) Resistivity increase was observed during the “stress reversal” period for low-saturation rock samples. The results mentioned above were explained with the effect of water flowing in and out of the cracks of rock. The temporary factors which yield a reduction of the maximum main stress, may enhence the possibility of earthquake occurrence.     

DIRDOP: a directivity approach to determining the seismic rupture velocity vector

Directivity effects are a characteristic of seismic source finiteness and are a consequence of the rupture spread in preferential directions. These effects are manifested through seismic spectral deviations as a function of the observation location. The directivity by Doppler effect method permits estimation of the directions and rupture velocities, beginning from the duration of common pulses, which are identified in waveforms or relative source time functions. The general model of directivity that supports the method presented here is a Doppler analysis based on a kinematic source model of rupture (Haskell, Bull Seismol Soc Am 54:1811–1841, 1964) and a structural medium with spherical symmetry. To evaluate its performance, we subjected the method to a series of tests with synthetic data obtained from ten typical seismic ruptures. The experimental conditions studied correspond with scenarios of simple and complex, unilaterally and bilaterally extended ruptures with different mechanisms and datasets with different levels of azimuthal coverage. The obtained results generally agree with the expected values. We also present four real case studies, applying the method to the following earthquakes: Arequipa, Peru (M _w = 8.4, June 23, 2001); Denali, AK, USA (M _w = 7.8; November 3, 2002); Zemmouri–Boumerdes, Algeria (M _w = 6.8, May 21, 2003); and Sumatra, Indonesia (M _w = 9.3, December 26, 2004). The results obtained from the dataset of the four earthquakes agreed, in general, with the values presented by other authors using different methods and data.     

A new fault rupture scenario for the 2003 Mw 6.6 Bam earthquake, SE Iran: Insights from the high-resolution QuickBird imagery and field observations

The December 26, 2003 M_w 6.6 Bam earthquake is one of the most disastrous earthquakes in Iran. QuickBird panchromatic and multispectral satellite imagery with 61 cm and 2.4 m ground resolution, respectively provide new insights into the surface rupturing process associated with this earthquake. The results indicate that this earthquake produced a 2–5 km-wide surface rupture zone with a complex geometric pattern. A 10-km-long surface rupture zone developed along the pre-existing Bam fault trace. Two additional surface rupture zones, each 2–5 km long, are oblique to the pre-existing Bam fault in angles of 20–35°. An analysis of geometric and geomorphic features also shows that movement on the Bam fault is mainly right-lateral motion with some compressional component. This interpretation is consistent with field investigations, analysis of aftershocks as well as teleseismic inversion. Therefore, we suggest that the 2003 Bam earthquake occurred on the Bam fault, and that the surface ruptures oblique to the Bam fault are caused by secondary faulting such as synthetic shears (Reidel shears). Our fault model for the Bam earthquake provides a new tectonic scenario for explaining complex surface deformations associated with the Bam earthquake.     

Features and physical process of the dynamic evolution pattern of ground resistivity precursor front

ＩｎｔｒｏｄｕｃｔｉｏｎＳｉｎｃｅｔｈｅ１９６０′ｓ,ｔｈｅｄｅｖｅｌｏｐｍｅｎｔｏｆｍｏｄｅｒｎｓｃｉｅｎｃｅａｎｄｔｅｃｈｎｏｌｏｇｙｈａｓｇｒａｄｕａｌｙｍａｄｅｉｔｐｏｓｉｂｌｅｔｏｐｒｅｄｉｃｔｅａｒｔｈｑｕａｋｅｓａｎｄｈａｓｐｒｏｍｏｔ...     

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.