Full moment tensor inversion constrained by double-couple focal mechanism for induced seismicity*

In this study, we propose a new method to determine full moment tensor solution for induced seismicity. This method generalizes the full waveform matching algorithm we have developed to determine the double-couple (DC) focal mechanism based on the neighbourhood algorithm. One major difference between the new method and the former one is that we adopt a new misfit function to constrain the candidate moment tensor solutions with respect to a reference DC solution in addition to other misfit terms characterizing the waveform matching. Through synthetic tests using a real passive seismic survey geometry, the results show the new constraint can help better recover the DC components of inverted moment tensors. We further investigate how errors in the velocity model and source location affect the moment tensor solution. The synthetic test results indicate that the constrained inversion is robust in recovering both the DC and non-DC components. We also test the proposed method on several real induced events in an oil/gas field in Oman using the same observation system as synthetic tests. While it is found that the full moment tensor solutions without using the DC constraints have much larger non-DC components than solutions with the DC constraints, both solutions are able to fit the observed waveforms at similar levels. The synthetic and real test results suggest the proposed DC constrained inversion method can reliably retrieve full moment tensor solutions for the induced seismicity.     

长江三峡工程周边地区的采矿诱发地震

地震观测和实地调查表明,长江三峡工程周缘广泛存在采矿诱发地震。采矿诱发地震主要是由于矿山采空区的出现使浅部地壳差应力增大岩体失稳而产生的,它们不同于天然地震。水库蓄水前应对这些地震的背景进行详细研究。在进行区域地震震情分析和预报、区域应力场及活动断裂研究时应严格加以区别     

甘肃平凉地区地震活动特征及其发震机理研究

针对2001年以来甘肃省平凉地区中小地震活动增强现象比较突出问题,通过遗传算法重新定位地震深度,利用矩张量反演方法反演震源机制解,并根据平凉地区地震时空分布特征、震相特征以及区域地质构造环境,认为平凉地区地震主要为煤矿诱发地震,可能发震机理为煤体压缩、顶底板拉伸、断层矿震等.     

Source Parameters and Scaling Relations for Reservoir Induced Seismicity in the Longtan Reservoir Area

We investigate the relationship between the impoundment and seismicity in the Longtan reservoir, southwestern China and find evidence that the seismicity was reservoir induced. After the reservoir impoundment, a pronounced increase in seismicity was observed in five clusters mainly concentrated in the areas where few earthquakes had occurred before the first filling. The observed induced seismicity shows a strong correlation with the filling cycles. The activity levels in the five clusters are different due to differences in the structures and permeabilities of the faults. Source parameters for 1,616 earthquakes with M _L 0.1–4.2 recorded by 24 fixed and temporary stations deployed around the reservoir were calculated after applying corrections for geometrical spreading, frequency-dependent Q, and site effects. The static stress drop and apparent stress in this area both appear to increase with increasing seismic moment over the entire magnitude range. Our results show that reservoir induced earthquakes have ten times lower average stress drop than natural tectonic earthquakes. These results may indicate that the reservoir induced seismicity can occur with a lower tectonic stress due to the high pore pressures of the underground medium, and that the effect of the water decreases the coefficient of friction.     

Evidence that earthquakes have been triggered by reservoir in the Song Tranh 2 region,Vietnam

In the Song Tranh 2 (ST2) hydropower reservoir located in the Quang Nam province, Central Vietnam, earthquakes started occurring soon after impoundment of the reservoir in late 2010. Earthquakes continue to occur in the region, and two earthquakes of M 4.6 and 4.7 on October 22, 2012 and November 15, 2012, respectively, have been reported (Trieu et al. 2014; Giang et al. 2015) in the vicinity of the reservoir. In the present study, b-value has been estimated, and focal mechanism solutions have been computed for the first time using moment tensor inversion approach. Also, the influence of impoundment of reservoir on the occurrence of earthquakes has been computed for the ST2 region based on Coulomb stress. A quality data set of 595 earthquakes recorded for the period of October 2012 to April 2014 at ten stations of the seismic network operated by the Institute of Geophysics (IGP) has been used to calculate b-values for the northern and southern seismicity clusters of the region. In general, the b-values associated with reservoir-triggered seismicity (RTS) are found to be higher than the regional b-values in the frequency-magnitude relation of earthquakes. For the ST2 region, it is found that the b-values for the northern and southern clusters are 0.94 ± 0.04 and 0.90 ± 0.04, respectively. Focal mechanism solutions obtained for the two earthquakes close to the reservoir have a right-lateral strike-slip mechanism, with the preferred planes trending NW-SE. These results are concurrent with the orientation of the nearby local surface faults, which we confirm as the active faults in this region. Influence of the stresses due to reservoir water load on the local seismicity is computed based on the obtained focal mechanism by using the concept of fault stability. It is found that most of the earthquakes occur in the positive Coulomb stress region, which shows the influence of reservoir impoundment on earthquake occurrence in the vicinity. Our results suggest that the local earthquakes are triggered by the impoundment of the ST2 reservoir.     

Modelling short-period crustal phases at regional distances for the seismic source parameter inversion

Waveforms of short-period crustal phases observed at regional distances are modelled in an attempt to retrieve the source parameters of small earthquakes. For a simple shield type crustal structure, the synthetic seismograms calculated for elementary source mechanisms show that the waveforms of crustal phases are quite distinct for different source mechanisms, and suggest that these features can be used to constrain the source mechanisms. Waveforms of the crustal phases are also strongly affected even by small variations in the crustal thickness as well as in the focal depth.A generalized, linearized waveform inversion technique is used to invert the regional crustal phases for a deviatoric moment tensor source. The inversion technique is applied to regional short-period records from a small earthquake in central Sweden for determining its source parameters. The waveform inversion, which mainly used Lg waves, produced a satisfactory fit between the observed and synthetic records, suggesting that the source mechanism and focal depth estimates were reasonable. The focal depth of 40 km obtained for the event studied suggests the occurrence of a deep crustal earthquake and demonstrates the usefulness of crustal phases in studying small intraplate earthquakes.     

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

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

Moment tensor inversion of nine events in Iran using INSN data

In this research, moment tensor (MT) inversion procedure in frequency domain is applied to nine events in order to evaluate the data efficiencies in computing source parameters of earthquakes in Iran. The events are recorded by Iranian National Seismological Network (INSN) stations, currently comprising 22 broadband stations. For the purpose of this study, first, the inversion procedure is applied to obtain source parameters. Then, the results obtained here are compared with those of Harvard global centroid moment tensors (GCMT) in order to estimate their reliabilities. The results show reasonable consistencies with those of Harvard GCMT within 0.075 uncertainties in logarithmic scalar moment or 0.05 in moment magnitude and 20° in Kagan angle for focal angles. Furthermore, according to the present research, the calculated source parameters are not significantly sensitive to the azimuthal gap.     

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

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

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

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

Earthquake detection probability within a seismically quiet area: application to the Bruchsal geothermal field

In applications such as oil and gas production, deep geothermal energy production, underground storage, and mining, it is common practice to implement local seismic networks to monitor and to mitigate induced seismicity. For this purpose, it is crucial to determine the capability of the network to detect a seismic event of predefined magnitude in the target area. The determination of the magnitude of completeness of a network is particularly required to properly interpret seismic monitoring results. We propose a method to compute the detection probability for existing local seismic networks, which (i) strictly follows the applied detection sequence; (ii) estimates the detection capability where seismicity has not yet occurred; and (iii) delivers the results in terms of probabilities. The procedure includes a calibration of a local magnitude scale using regional earthquakes recorded by the network and located outside the monitored area. It involves pre‐processing of the seismograms recorded at each station as performed during the triggering sequence, which is assumed based on amplitude thresholds. Then, the calibrated magnitude–distance–amplitude relations are extrapolated at short distances and combined to reproduce the network detection sequence. This generates a probability to detect a seismic event of a given magnitude at a specified location. This observation‐based approach is an alternative to a fully theoretical detection capability modelling and includes field conditions. Seismic wave attenuation by geometrical spreading and intrinsic attenuation, site effect, and instrumental responses are partly accounted for by the calibration. We apply this procedure on the seismic network deployed in the Bruchsal geothermal field (Germany). Although the system was in good working order, no induced seismicity was identified in the area between June 2010, when monitoring started, and November 2012. The recording of distant seismicity during this time period, however, allowed the application of the proposed procedure. According to the applied network detection parameters, the results indicate that the absence of seismicity can be interpreted as a 95% probability that no seismic event with M_L ≥ 0.7 occurred below the network at 2.4‐km depth, i.e., in the geothermal reservoir.     

汶川MS8.0级地震发生背景与过程的研究

本文首先阐明汶川M_S8.0级地震发生在由区域布格重力异常和地震震中分布所确定的武都—松潘—茂汶—汶川—泸定地震带上.汶川地震所在地段是地震前兆和中小地震（M≤7.0）的空白区,震前出现明显的孕震空区,M_S8.0级地震发生在空区周围区域中小地震活动峰值之后的减少段里.地震的破裂超出孕震空区范围,空区内、外余震活动呈现出不同的衰减特征,依此将余震活动分为WS和NE两个区段.地震破裂过程、4级以上余震矩张量及震区应力场反演和余震应力降的测定结果表明,两个区域的位错、余震机制解和应力降及最大主应力的方〖JP2〗向等明显有别.根据这些特征和地震应力触发的研究,推测NE段地震的发生可能是〖JP〗由WS段主破裂的发生所触发.     

应用近场速度记录快速估算川滇地区中强地震强度

采用Wu等(2004)提出的地震强度测定方法,应用近场速度波形记录,测定地震参数M_(ew),实现了川滇地区中强地震(M_W≥5.7)强度近实时快速测定。M_(ew)是一个与矩震级M_W对接的地震参数,与传统的应用振幅测定震级的方法相比,不会存在震级饱和问题,与全球矩心矩张量(GCMT)的矩震级M_W具有较好的一致性。该方法解决了应用近场波形资料和矩张量反演方法求取大地震矩震级比较困难的问题,可为区域中强地震速报提供近实时矩震级提供参考,从而提高地震速报的时效性和准确性,同时为强震加速度数据的应用奠定基础。     

Source parameters determination for earthquakes in Kushiro,Japan considering source time function

This paper applies a new formulation to do moment tensor inversion for earthquakes in the Kushiro area of Japan. Comparing with conventional moment tensor inversion method, the new one takes the effect of source time function into consideration. For the inversion, best solution is obtained by minimizing the difference between the observed seismograms and the synthetic ones. And the best-fitting focal depth is determined from the variance reduction. The results indicate that half duration of source time function is proportional to the magnitude of earthquakes. Large earthquakes have long half duration, whereas that of moderate-small earthquakes is comparatively shorter. The focal mechanisms of all three earthquakes are of thrust fault type, which is mainly ascribed to the collision of the North American plate with the Eurasia plate in the late Cretaceous or Paleogene.     

二滩水库蓄水前天然地震活动本底背景

了解和掌握水库蓄水前库坝区及其周围天然地震背景水平,不仅可以对水库诱震的危险性进行前期评价,而且还可以为蓄水后的诱发地震活动监测提供可靠的鉴别依据。二滩水库遥测地震台网在水库蓄水前已连续运行了６年,积累了世界水库地震监测史上极为难得的资料。为此,本文利用这批资料提供了表征库坝区蓄水前天然地震活动水平确切实用的各项指标,同时也使用四川台网的地震资料,给出了包括库区的较大区域地震活动背景,讨论了对库区可能的影响。     

中国台湾以东地区地震矩张量研究及其构造意义

本文利用长周期Ｐ波波形反演方法，推断了１９８６年１０月至１９９１年期间发生在中国台湾岛东部及其附近海域中８次中强地震的矩张量．矩张量解揭示了发生在这一区域边界处的地震震源性质的规律性变化：沿西部边界的主旋错动，沿北部边界的右旋错动，西北角花莲地区的逆冲断错及显著的体积收缩分量．这一结果表明，台湾地区是菲律宾海板块与欧洲大陆板块边界中的一特殊地段，海洋板块的俯冲在这里受阻．随着菲律宾海板块向西北方向的相对运动，在花莲地区发生了强烈碰撞．     

Moment Tensor Solutions and Triggering Environment for Earthquakes in Koyna-Warna Water Reservoirs Region,India

We consider nine earthquakes in the Koyna-Warna reservoir region on the western side of the Peninsular India. The deviatoric moment tensors of these earthquakes have been evaluated by minimizing the least-squares misfit between observed and synthetic seismograms. We use broadband seismograms of observatories at KARD and PUNE which are at distances of nearly 50 and 150 km, respectively, from the epicenters. Both surface wave inversion and the difference between the arrival times of SH and SV show the presence of an anisotropic crust. However, we have obtained an equivalent isotropic structure by improving the published crustal structures of this area through inversion of surface wave group velocity data. The deviatoric moment tensors of the earthquakes are decomposed into two components: double-couple and compensated linear vector dipoles (CLVD). The double-couple components of all the nine earthquakes show normal faulting with minor strike slip; the T axis is consistently subhorizontal with an average azimuth of 260.6° and the P axis is nearly vertical. The fault planes of six events give average strike direction and dip, respectively as 194.0° and 51.8° and are associated with the main fault of the area. The other three events lie in the southern part of this area and have strike direction between SSE and SE which is parallel to the tectonic features in this part. The CLVD component is generally within 20 percent of the total moment tensor. Recent studies show that anisotropy can produce source mechanism with CLVD up to 30 percent and can also cause high pore fluid pressure leading to fault instability more rapidly compared to conventional mechanism in an isotropic medium. It appears that the anisotropic crust, noted in the present work, is generating the CLVD component and also gives the proper environment to trigger earthquakes by reservoirs through pore fluid pressure.     

三峡库区秭归与巴东交界Ms4.0地震矩张量及应力环境研究

三峡库区湖北秭归与巴东交界地区2013年以来发生5次Ms4.0以上地震.本文基于湖北、重庆等省属13个测震台站宽频带波形数据,采用kiwi软件包和rapidinv12程序反演了这些地震的全矩张量解,结果显示所有地震双力偶解均为走滑兼逆断性质,两个节面走向、倾角、滑动角与P、T轴走向和倾俯角的一致性均较好,主压应力P轴N...     

水口水库地震对南平地区的影响

水口水库建成蓄水后诱发了频繁的地震活动，延续时间已达8年之久，,学在继续活动。近年来，南平地区的邵武，政和等地也相继有中小地震活动，为了摸清南平地区地震活动与水口水库诱发地震是否关联，本文对该地区及水口水库诱发地震前后的地震活动性及时空特征，区域构造背景及水库诱发地震的活动趋势进行综合分析研究。     

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.