EXPERIMENTAL STUDY ON THE CHANGES OF ULTRASONIC CODA WAVE AND ACOUSTIC EMISSION DURING ROCK LOADING AND DEFORMATION

The coda wave propagation path has received extensive attention as it is more sensitive to small changes in the medium than the direct wave. During the process of loading, the wave velocity, medium or source changes may cause the coda wave to change. The physical mechanism of change in the ultrasonic coda wave varies during different deformation stages. Meanwhile, there exist local damages in the rock sample during the deformation, and it will be accompanied by acoustic emission. Combining the ultrasonic coda wave and acoustic emission is beneficial to characterize the coda wave characteristics and damage degree of the sample at different deformation stages. In this paper, three kinds of rocks, including granodiorite, marble and sandstone with the sizes of 50mm×50mm×150mm, are used to carry out observations of ultrasonic coda wave and acoustic emission during the whole process of loading so as to study characteristics of the coda wave at different deformation stages. The major results are given below: 1)There is a good correspondence between the coda wave variation and the acoustic emission evolution process. When the acoustic emission frequency increases, the coda wave changes accordingly. In particular, the coda wave changes in the early stages of increased acoustic emission frequency, which indicates that the early damage information of rock can be obtained by analysis of the coda wave. 2)The physical mechanism of the coda wave change is different in different deformation stages. At the initial stage of loading, there are obvious scatterer changes in the coda wave change; then, in the linear elastic deformation stage, the wave velocity change is dominating; in the late-stage of loading, the scatterer change increases and coexists with the wave velocity change, the scatterer change effect is related with the rock micro-fracture degree, the rock will locally be damaged before rupturing, and the role of the scatterer will be enhanced. 3)With the increase of loading, the amplitude of increase of the wave velocity generally decreases gradually, which is basically consistent with the understanding obtained through the direct wave. The interference of acoustic emission can be eliminated because of the Kaiser effect when analyzing the coda wave. The consistency of the wave velocity change and stress loading and unloading is further verified. 4)The micro-fracture generated during rock deformation will change the physical mechanism of the coda wave change, and the scatterer effect will be significantly enhanced. At the same time, the acoustic emission waveform will cause interference to the ultrasonic coda wave. This means that attention needs to be paid when analyzing rock damage using only coda wave data. In short, the ultrasonic coda wave and acoustic emission can reflect the damage inside the rock, and the change mechanism of the coda wave in different deformation stages is different. The joint observation of the two can play a mutual verification role, which is conducive to improving the reliability of the observation results.     

基于理论地震图方法的波导结构对 地下核爆炸Lg波传播特性研究

基于频率-波数域算法的理论地震波形图方法, 可以数值模拟频率达到10 Hz、 震中距达1000 km的区域理论地震波形图. 该算法适用于计算大量分层地壳结构中激发的导波Lg波. 本文在前苏联东哈萨克斯坦地下核试验场至我国乌鲁木齐台站间的地球介质速度模型中, 引入速度梯度结构、 速度扰动分布的薄叠加层结构、 降低Q值结构以及速度扰动与Q值变化的综合结构来模拟实际地壳波导结构, 较好地模拟出东哈萨克斯坦地下核爆炸地震在乌鲁木齐台站记录的宽频带地震波形图, 模拟出完整的Lg波序列, 该序列符合Lg波能量分布特征, 且能够解释Lg波波尾的特征. 结果表明, Lg波的形状和峰值结构均依赖于地壳的不同波导结构.      

利用两种模型对宁夏及邻区尾波 Q 值进行对比研究

分别利用Aki模型和Sato模型,选取2008年1月-2009年12月宁夏地震台网记录的ML≥2.0地震的数字地震波资料,计算了宁夏及邻区平均尾波Q值,并拟合了Q值对频率的依赖关系.结果为:Aki模型为Q(f)=(212±87.62)f0.7584±0.19;Sato模型为Q(f)=(44.81±15.87)f0.9491±0.0805.与国内其他区域相比,两种模型计算结果均显示本区域Q值较低,对频率依赖性较高.研究结果认为两种模型均适合对宁夏及邻区尾波Q值的研究.     

Scattered Surface Wave Energy in the Seismic Coda

One of the many important contributions that Aki has made to seismology pertains to the origin of coda waves (Aki, 1969; Aki and Chouet, 1975). In this paper, I revisit Aki's original idea of the role of scattered surface waves in the seismic coda. Based on the radiative transfer theory, I developed a new set of scattered wave energy equations by including scattered surface waves and body wave to surface wave scattering conversions. The work is an extended study of Zeng et al. (1991), Zeng (1993) and Sato (1994a) on multiple isotropic-scattering, and may shed new insight into the seismic coda wave interpretation. The scattering equations are solved numerically by first discretizing the model at regular grids and then solving the linear integral equations iteratively. The results show that scattered wave energy can be well approximated by body-wave to body wave scattering at earlier arrival times and short distances. At long distances from the source, scattered surface waves dominate scattered body waves at surface stations. Since surface waves are 2-D propagating waves, their scattered energies should in theory follow a common decay curve. The observed common decay trends on seismic coda of local earthquake recordings particular at long lapse times suggest that perhaps later seismic codas are dominated by scattered surface waves. When efficient body wave to surface wave conversion mechanisms are present in the shallow crustal layers, such as soft sediment layers, the scattered surface waves dominate the seismic coda at even early arrival times for shallow sources and at later arrival times for deeper events.     

The complexity of the crust from refraction/wide-angle reflection data

ThePmP wide-angle reflected signal from the Moho shows a wide variation in its characteristics when cornpared from region to region. If the earth's crust is simple and homogeneous and the Moho is a sharp discontinuity, supercriticalPmP wavelets are large and isolated in the 100 km to 300 km distance range when compared to the preceding signals on the seismic traces. If the crust itself has numerous short reflectors and the Moho transition zone is thick or badly disrupted from past tectonic disturbances, thenPmP is often poorly defined and difficult to identify in the coda of earlier arriving signals. These signals are poorly correlated because the reflectors themselves are in general not only discontinuous but tend to be distributed at various depths within the crust. The effect of the vertical velocity gradient in the crust is to make the reflective field for downgoing waves much greater than for upgoing waves. The large variations in reflection coefficients wiht angle of incidence has the effect of making the coda generated by a reflective machanism distance dependent. The reflected signals are also contaminated with scattered signals from smaller scale heterogeneities which may be distributed more uniformly within the crust.In conventional seismic inversion methods applied to crustal refraction experiments, emphasis is usually placed on obtaining a velocity model of the subsurface structure. These models are relatively simple and are limited in the amount of complexity which can be uniquely inferred from the data. The main problem arises because of the difficulties in the identification of the origin of the signals which are themselves often incoherent from trace to trace.In this paper it is shown how the conventional record section may be complemented with a normal moveout corrected intensity section which emphasizes areas of large signal complexities. Data from this intensity section is then used as input to obtain a quantitative measurement of a complexity parameter. These measurements may be used to infer or compare differences in crustal heterogeneity from one region to another. The discussion is illustrated with both numerical modeling data as well as data from recent crustal experiments which were conducted over the Canadian Shield.