2013年前郭MS5.8震群矩张量研究

2013年吉林前郭M_S5.8震群为爆发性震群,目前余震活动仍然在持续.基于吉林、辽宁、黑龙江、内蒙古四省地震台网记录的前郭震群波形资料,利用波形信噪比、震源类型、台站及速度模型组合的指标选择最佳的反演方案,应用矩张量的三种反演模式,对序列中5个M_S≥5.0地震进行矩张量反演研究,获得了全矩张量、偏矩张量和纯双力偶的矩张量.使用F-test对地震的三种模式的矩张量反演结果进行显著性检验来确定最佳反演模式.结果显示,5个地震的最优矩张量解均为全矩张量模式反演获得的结果,其双力偶分量仅有20%~65%,矩心深度位于地下3~4 km处,地震在Hudson震源类型图上的投影远离双力偶震源类型区域.这些结果表明,震源类型并非典型的构造地震,推断前郭地震可能是与人类活动有关的诱发地震.     

Decomposition of seismic moment tensors for underground nuclear explosions

Generally the decomposition of a seismic moment tensor is not unique. However, to favorably view the characteristics of a certain seismic source, one must decompose a seismic moment tensor into parts according to assumptions about the properties of the seismic source. Different from natural earthquakes in which the shear dislocation component plays a predominant role in the source process, and the seismic moment tensor can be separated into an isotropic component, a double couple, and a compensated linear vector dipole (CLVD), underground nuclear explosions have three major components in their source process, i.e., the explosion, the tensional spalling, and the tectonic strain release associated with the explosion. In such a situation the conventional moment tensor decomposition for earthquakes is not convenient to estimate the yield of the explosion and to characterize the tectonic strain release. In this paper, an alternative decomposition scheme is proposed to deal with the moment tensor of underground nuclear explosions, which might benefit the approach to study the tectonic strain release induced by underground nuclear detonations.     

非对称地震矩张量时间域反演:理论与方法

本文在对称地震矩张量反演的基础上,进一步研究了非对称地震矩张量时间域反演的理论与方法,结果表明:非对称地震矩张量反演与对称地震矩张量反演类似,只需将对称地震矩张量反演方法略加改动,即增加3个待解参数,便可实现非对称地震矩张量反演.为了判断非对称地震矩张量反演相对于对称地震矩张量反演是否存在过度拟合,运用了AIC准则 (赤池信息准则).为了定量地描述地震矩张量之间的差异,引入了地震矩张量的矢量表示法.通过分析格林函数与地震矩张量各分量之间的关系,得出:在非对称地震矩张量反演时,若仅用垂直向地动位移数据,将无法区分M_xy与M_yx这两个分量, 需要同时运用垂直向与水平向地动位移数据进行联合反演才能区分M_xy与M_yx; 若采用不同的速度结构模型或不同的格林函数计算方法,则需重新评估地震矩张量各分量的分辨度问题.为检验非对称地震矩张量反演方法的可行性, 利用合成地震图进行了一系列数值试验.数值试验结果表明,在非对称地震矩张量反演中,有必要引入S波进行P波与S波联合反演以提高反演的准确性和判定断层面的能力.      

On the observability of isotropic seismic sources: The July 31, 1970 Colombian earthquake

Theoretical calculations are made to study the observability of isotropic components of seismic sources. In particular we consider the 1970 deep Colombian earthquake, for which a precursory isotropic component was previously reported by Gilbert and Dziewonski.We compare an ultra-long period vertical record at Pasadena of the 1970 event to synthetic seismograms calculated both for Gilbert and Dziewonski's source model and for the pure double-couple source of Furumoto and Fukao, and obtain better overall agreement for the latter. The amplitude of the long-period synthetic for the isotropic source is about 5–15 times smaller than the synthetic for the deviatoric source, suggesting that the data may be relatively insensitive to the presence of a small isotropic source. When this possibility was tested, the overall agreement was found to be almost completely insensitive to the presence of even a reasonably large isotropic component.However, the isotropic source was derived from multi-station moment tensor inversion, rather than from single-station studies. A numerical experiment on the effect of lateral heterogeneity of eigenfrequencies and of Q on the inversion for the moment tensor shows that even relatively small amounts of heterogeneity can produce spurious isotropic sources from moment tensor inversion.     

Ambiguous Moment Tensors and Radiation Patterns in Anisotropic Media with Applications to the Modeling of Earthquake Mechanisms in W-Bohemia

Anisotropic material properties are usually neglected during inversions for source parameters of earthquakes. In general anisotropic media, however, moment tensors for pure-shear sources can exhibit significant non-double-couple components. Such effects may be erroneously interpreted as an indication for volumetric changes at the source. Here we investigate effects of anisotropy on seismic moment tensors and radiation patterns for pure-shear and tensile-type sources. Anisotropy can significantly influence the interpretation of the source mechanisms. For example, the orientation of the slip within the fault plane may affect the total seismic moment. Also, moment tensors due to pure-shear and tensile faulting can have similar characteristics depending on the orientation of the elastic tensor. Furthermore, the tensile nature of an earthquake can be obscured by near-source anisotropic properties. As an application, we consider effects of inhomogeneous anisotropic properties on the seismic moment tensor and the radiation patterns of a selected type of micro-earthquakes observed in W-Bohemia. The combined effects of near-source and along-path anisotropy cause characteristic amplitude distortions of the P, S1 and S2 waves. However, the modeling suggests that neither homogeneous nor inhomogeneous anisotropic properties alone can explain the observed large non-double-couple components.The results also indicate that a correct analysis of the source mechanism, in principle, is achievable by application of anisotropic moment tensor inversion.     

用应变地震观测求解震源矩张量的基本原理

震源机制解,即对地震矩张量的推断,对于地震研究具有至关重要的意义.应变地震观测是张量观测,与摆式地震仪的位移矢量观测不同,可以为地震研究提供新的数据源.本文讨论用应变地震观测求解震源矩张量的基本原理.在距离震源足够远的地方,地震波可以看成平面波,其性质决定于震源矩张量.假设平面地震波的应变张量可以由震源矩张量通过坐标变换计算得到,就可以通过观测应变地震波求解震源机制.这个假设至少对于双力偶震源机制是成立的.由此可以证明,在理想的无限介质中,只要有两个以上不同地点的应变地震波观测,就可以解出震源矩张量.这为解决震源机制问题提供了新的方法.目前的地震矩张量求解方法需要两方面的条件:或者需要很多观测点(例如体波反演),或者需要长周期地震波资料(例如面波反演).这些方法只适用于分析比较大的地震.对于小震,因为通常其震中周围不会有足够多的摆式地震仪观测点观测到其地震波,而地震波周期又短,难以利用传统方法给出可靠的震源机制解,所以只需少数观测点就能求解震源矩张量的新方法特别有意义.用应变地震观测求解震源机制,可以给出更为精确的结果.     

断层厚度的地震效应和非对称地震矩张量

本文导出了具有厚度和滑动弱化区域的断层的非对称地震矩张量表示式,指出要求地震矩张量具有对称性不是一个绝对必要的限制.在非对称地震矩张量中,位错项对应于对称地震矩张量,拉力项对应于非对称地震矩张量.由于拉力项等效于单力偶,所以在非对称地震矩张量解的两个节面上,沿滑动矢量方向的力偶强度不再相同,与较大力偶相联系的节面为断层面,与较小力偶相联系的节面为辅助面.这一性质可用以从两个正交的节面中判断哪一个节面是断层面.如果忽略拉力项,会高估与位错对应的标量地震矩.只有满足相应的约束条件的非对称地震矩张量才能表示具有厚度和滑动弱化区域的断层模型,并从中分离出与位错和拉力对应的地震矩张量.     

Isotropic and Nonisotropic Components of Earthquakes and Nuclear Explosions on the Lop Nor Test Site, China

—?We test the hypothesis that the existence of an observable non-zero isotropic component of seismic moment can be used as a discriminant to distinguish nuclear explosions from shallow earthquakes. We do this by applying the method described herein to a small set of data recorded between 1990 and 1996 following events (seven nuclear explosions, three earthquakes) that occurred on the Lop Nor test site in Western China. We represent each source as a sum of an isotropic component at the surface and a nonisotropic, double-couple component at an estimated depth. The explosions all possess a significant non-zero isotropic component and the estimated depth of the double-couple component of the moment tensor, presumably the result of tectonic release, lies between about 0 and 3?km. For the earthquakes studied, the isotropic component is indistinguishable from zero and the depths of the sources are estimated at 3, 17 and 31?km. The data set we have studied, although still very small, suggests that certain source characteristics (namely, double-couple depth and the ratio of the isotropic to nonisotropic components of seismic moment) may prove useful in discriminating explosions from shallow earthquakes. Further work is needed to determine whether these observations hold for explosions at other test sites, to investigate a much larger set of shallow earthquakes located in regions of interest, and to study the robustness of the estimated source parameters as source magnitude and the number of observing stations decrease.     

Inversion of waveform from the data of six moderate strong earthquakes in Chinese mainland to seismic moment tensors and source mechanism

Based on Generalized Seismic Ray Theory (Helmberger, 1968), a new quickly linear inversion method from the data of seismic waveform to seismic moment tensor and source mechanism for domestic earthquake is studied in this paper. Six moderately strong earthquakes which occurred in Chinese mainland in the past few years are studied. The seismic source parameters of these earthquakes, seismic moment tensors, scalar seismic moments, fault plane solutions and source time functionsetc, are obtained. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,15, 261–268, 1993.     

A preliminary empirical model for frequency-dependent attenuation of Fourier amplitude spectra In Serbia from the Vrancea earthquakes

We present the frequency-dependent attenuation model for Fourier amplitude spectra of strong earthquake ground motion in Serbia from intermediate depth earthquakes in the Vrancea source zone in Romania. The development of this type of scaling is the essential first step toward developing the corresponding attenuation and scaling equations for pseudo relative velocity spectra (PSV), which are necessary for seismic macro- and microzoning in the territory of Serbia. Such scaling is necessary because the Vrancea source zone produces large earthquakes with shaking that attenuates differently from the local earthquakes in Serbia. Development of such a scaling model is associated with several difficulties, the principal one being the lack of recorded strong motion accelerograms at epicentral distances exceeding 300 km. To reduce uncertainties with such scaling, we require our preliminary scaling equations to be consistent with independent estimates of seismic moment, stress drop, and radiated wave energy. In the future, when the recorded strong motion data from Vrancea earthquakes increases several-fold of what it is today, it will become possible to perform this analysis again, thus leading to more reliable and permanent scaling estimates.     

Moment Tensor Inversion of Regional Phases: Application to a Mine Collapse

—?The aim of our study consists of analyzing potentially non-double-couple seismic events recorded at regional distances. In order to define the nature of the seismic source, a moment tensor inversion is carried out as this method is general enough not to initially constrain the source mechanism. In this paper we present an application to a seismic event induced by a mine collapse which occurred near the town of Halle in Germany. Because of its induced nature, many parameters such as the location and geometry of this seismic source are known. This information allows us to test the influence of inadequate propagation modeling on the moment tensor obtained from the inversion. Green's functions have been computed with the reflectivity method in a flat layered medium, using the European model EurID (Du et? al., 1998; Dufumier et al., 1997). From the inversion of P-wave seismograms recorded by the German Regional Seismic Network will, we obtained a source time function which can be decomposed into two subevents. The first one has a large isotropic part and a deviatoric mechanism with near vertical nodal planes. No volume change is observed for the second subevent, but a deviatoric component opposite of the first one. The addition of S-waves does not change the results of the inversion which are stable. Surface waves were not used because of their poor dispersion curves. Based on the moment tensor obtained from these inversions, the physical process at the source is compatible with a large cavity collapse.     

用强震动观测资料估算2007年宁洱M_S6.4地震震源参数

通过对2007年6月3日宁洱M_S 6.4、M_S 5.1两次地震的强震动台站观测记录的谱分析,估算了两次地震的辐射能量、地震矩、矩震级等震源参数.结果表明,虽然对拐角频率、应力降等参数各台站的估算结果离散性较大,效果不甚理想,但利用近场强震动观测记录,可以对地震辐射能量、地震矩、矩震级3个参数作出较好的估计.     

区域台网地震矩张量快速反演系统研究

基于美国内华达大学地震实验室的矩张量反演程序,采用JAVA和HTML语言,研发了一套近实时全波形矩张量快速反演系统,用于区域台网地震矩张量的快速产出.该系统由区域地震速报系统触发,能自动实现中等强度地震矩张量的快速反演,同时也包含人机交互模块,可进行人工校正,最终以网页形式对外服务.本文将系统运行以来成功反演的结果与USGS和中国地震局地球物理研究所的计算结果进行对比,验证了该系统的可靠性.      

地震震源运动学参数获取方法研究进展

本文基于从简单到复杂的各向同性点源模型、质心矩张量模型、有限矩张量模型和有限震源滑动分布模型4种震源运动学模型,概述了相应的研究进展.     

Moment tensor inversion of small to moderate earthquakes in the Capital Region in 2004

The moment tensor solutions of 51 small to moderate earthquakes occurred in the Capital Region in the year of 2004 are obtained by inverting the broadband waveform data. Accordingly, other source parameters, such as scalar seismic moments, moment magnitudes, double-couple （DC） components and compensated-linear-vector-dipole （CLVD） components, are determined as well as fault parameters and stress-axis parameters. The inverted results are evaluated by groups of numerical tests.     

The determination of the best fitting and the extreme source models by a global inversion method

In this paper, a global inversion method is developed for seismic moment tensor inversion by using the body wave forms. The algorithm depends on neither the selection of starting model nor the forms of objective function and constraints. When the error function, measure of the difference between the observed and synthetic waveforms, is chosen as the objective function, the best fitting source model is found; when a certain combination in seismic moment tensor elements is selected as the objective function and the values of error function are constrained in a suitable bound, the extreme source models can be produced by minimizing or maximizing this combination. By changing the form of the combination of moment tensor elements, a variety of different source characteristics can be considered. Therefore the extreme solution provides an estimation of the uncertainty in the best fitting source model. The seismic waveform data was used to evaluate the effectiveness of this algorithm. This research was supported by the National Natural Science Foundation of China.     

Studies on seismic source

During the period of 1999~2002, the Chinese seismologists made a serious of developments in the study on seismic sources including observations, experiments and theory. In the field of observation, the methods of the accuracy location of earthquake sources, the inversion of seismic moment tensor and the mechanism of earthquake source are improved and developed. A lot of important earthquake events are studied by using these methods. The rupture processes of these events are inverted and investigated combined with the local stress fields and the tectonic moment by using the measurements of surface deformation. In the fields of experiments and theory, many developments are obtained in cause of seismic formation, condition of stress and tectonics, dynamics of earthquake rupture, rock fracture and nucleation of strong earthquakes.     

Studies on seismic source

During the period of 1999-2002, the Chinese seismologists made a serious of developments in the study on seis-mic sources including observations, experiments and theory. In the field of observation, the methods of the accu-racy location of earthquake sources, the inversion of seismic moment tensor and the mechanism of earthquake source are improved and developed. A lot of important earthquake events are studied by using these methods. The rupture processes of these events are inverted and investigated combined with the local stress fields and the tec-tonic moment by using the measurements of surface deformation. In the fields of experiments and theory, many developments are obtained in cause of seismic formation, condition of stress and tectonics, dynamics of earthquake rupture, rock fracture and nucleation of strong earthquakes.     

2004年首都圈地区中小地震的矩张量反演

通过首都圈地区宽频带数字地震波形反演, 获得了2004年发生在该地区的51次中小地震的矩张量解, 由此确定了这些地震的标量地震矩、 矩震级、 双力偶分量和补偿线性矢量偶极分量的大小以及断层面参数和应力轴参数, 并通过数值试验对反演结果进行了评价.