利用模板匹配技术识别大岗山水库库区爆破事件

采用匹配定位方法,利用大岗山水库库区及附近（29.2°—29.9°N,101.9°—102.5°E）5次已知爆破事件对2014年3月1日至2015年12月31日内连续波形进行模板检测与定位,经过人工复核后共获取23次疑爆事件,这些事件与模板事件间有较高的相关性,相关系数达0.773 4,匹配后疑爆事件数为模板数的4.6倍,提高了对疑爆事件的识别率,同时剔除了原地震目录中人工误判的8次疑爆事件,有效地提高了地震编目质量。将匹配出来的疑爆事件及距离最近的台站记录到的NS向波形分别在Google Earth上投影,从影像图上可见扫描出的事件附近存在采石场或蓄水前修建大坝的场地,较符合爆破实际情况,重定位偏离量基本在0.5 km之内。研究认为,利用匹配定位的方法,能有效地区分不同类型的地震活动。     

Tectonic setting of the seismogeny of the 1976 Tangshan, China M =7.8 earthquake

Based on the analysis of multi-temporal and multi-spectral satellite images for North China region, we have found that one year and more before the occurrence of 1976 Tangshan earthquake, the anomalies of electro-magnetic radiation on the satellite images indicated that the NE-trending Tangshan fault zone was dissected by the NNW-trending Nantai-Tangshan fault, and the Changping-Fengnan fault was dragged to form an arcuate bending at Fengnan. All these indicate the right-lateral translation along the Tangshan fault in NE direction. In order to gain an insight into the features of these faults, a shallow seismic exploration along the Tangshan and Changping-Fengnan faults has been carried out. The results have indicated that the NE-trending Tangshan fault is a high angle right-lateral strike-slip normal fault, dipping northwest, while the NWW-trending Changping-Fengnan fault is a southwest-dipping left-lateral strike-slip normal fault. The project is supported by the National Natural Science Foundation of China (No. 49672170).     

Rupture process of the M L=4.1 earthquake in Huailai Basin on July 20, 1995

On July 20, 1995, an earthquake of M _L=4.1 occurred in Huailai basin, northwest of Beijing, with epicenter coordinates 40.326°N, 115.448°E and focal depth 5.5 km. Following the main shock, seismicity sharply increased in the basin. This earthquake sequence was recorded by Sino-European Cooperative Huailai Digital Seismograph Network (HDSN) and the hypocentres were precisely located. About 2 hours after the occurrence of the main shock, a smaller event of M _L=2.0 took place at 40.323°N, 115.447°E with a focal depth of 5.0 km, which is very close to the main shock. Using the M _L=2.0 earthquake as an empirical Green’s function, a regularization method was applied to retrieve the far-field source-time function (STF) of the main shock. Considering the records of HDSN are the type of velocity, to depress high frequency noise, we removed instrument response from the records of the two events, then integrated them to get displacement seismogram before applying the regularization method. From the 5 field stations, P phases in vertical direction which mostly are about 0.5 s in length were used. The STFs obtained from each seismic phases are in good agreement, showing that the M _L=4.1 earthquake consisted of two events. STFs from each station demonstrate an obvious “seismic Doppler effect”. Assuming the nodal plane striking 37° and dipping 40°, determined by using P wave first motion data and aftershock distribution, is the fault plane, through a trial and error method, the following results were drawn: Both of the events lasted about 0.1 s, the rupture length of the first one is 0.5 km, longer than the second one which is 0.3 km, and the rupture velocity of the first event is 5.0 km/s, larger than that of the second one which is about 3.0 km/s; the second event took place 0.06 s later than the first one; on the fault plane, the first event ruptured in the direction γ=140° measured clockwise from the strike of the fault, while the second event ruptured at γ=80°, the initial point of the second one locates at γ=−100° and 0.52 km from the beginning point of the first one. Using far-field ground displacement spectrum measurement method, the following source parameters about the M _L=4.1 earthquake were also reached: the scalar earthquake moment is 3.3×10¹³ N·m, stress drop 4.6 MPa, rupture radius 0.16 km. Contribution No. 99FE2022, Institute of Geophysics, China Seismological Bureau. This study is supported by the Chinese Joint Seismological Science Foundation (95-07-411).     

Regional seismicity triggered by the Ms=7.8 Tangshan event on July 28, 1976 and the static stress field change

We studied the seismicity before and after the MS=7.8 Tangshan event on July 28, 1976 (39(28(N, 118(11(E) and the results show that in 3 regions outside of the source zone, seismicity rate increasing were observed, which was significant in 0.99 significance level using Z-statistic test and was proposed to be triggered by the MS=7.8 Tangshan earthquake. The magnitude of the greatest triggered event was 5.5. The epicenter distances of these earthquakes were several ten kilometers to 300 km. The static stress change (CFS of Coulomb failure was calculated using an elastic dislocation model in half space and the (CFS on the major rupture directions in these three regions were positive.     

The focal mechanism solution of the M s 6.5 Wuding earthquake in Yunnan Province of October 24, 1995

The focal mechanism of Wuding earthquake with magnitude M=6.5 on October 24, 1995 is estimated by the method of inversion of seismic tensor in this paper. The two principal axes are nearly horizontal, the principal compressive axis is about N30°E, the nodal plane which strikes about N105°E (N75°W) is probably consistent with the actual rupture plane.     

Synthesis of accelerations of the 1976 Tangshan earthquake(M_S 7.8)in near-and far-field by using semi-empirical method

Ｓｙｎｔｈｅｓｉｓｏｆａｃｃｅｌｅｒａｔｉｏｎｓｏｆｔｈｅ１９７６Ｔａｎｇｓｈａｎｅａｒｔｈｑｕａｋｅ（ＭＳ７．８）ｉｎｎｅａｒａｎｄｆａｒｆｉｅｌｄｂｙｕｓｉｎｇｓｅｍｉｅｍｐｉｒｉｃａｌｍｅｔｈｏｄＱＩＦＥＮＧＬＵＯ１）（罗奇峰）ＹＵ?..     

1976年唐山MS7.8地震余震序列持续时间及对地震危险性分析的意义

1976年7月28日唐山MS7.8大地震对唐山及其周边地区造成了重大的人员伤亡和财产损失. 主震之后约15小时滦县又发生了MS7.1地震; 同年11月15日宁河也发生了MS6.9地震. 唐山MS7.8主震后的余震一直持续至今,使该区域至今保持了与主震前相比具有较高的地震活动性．如何估计余震的持续时间,并进一步将余震从主震目录中去除,一直是地震学中所关注的问题．该文通过对数线性回归和理论计算,从不同角度求取并讨论了1976年唐山MS7.8大地震的余震持续时间．结果表明,由对数线性回归计算得到的余震持续时间约为80 a．而基于Dieterich的余震触发理论所得到的余震持续时间则与区域剪应力变化率有关．区域剪应力变化率可有几种不同方法求得: ① 根据剪应力变化率和静态应力降Delta;tau;e及地震回复周期tr之间的关系求取应力变化率,该方法所得到的余震持续时间约为70——100 a;② Ziv和Rubin对Dieterich的方法进行了修正,给出了通过远场加载速率和断层宽度求取应力变化率, 该方法得到的余震持续时间约为80 a;③ 由背景场地震活动性求取远场剪应力速率, 可以得到该地区二维分布式的余震持续时间,此方法得到的研究区域内余震持续时间为130——160 a．综上,唐山地区余震持续时间约为70——140 a,据此, 该地区现今所发生的地震仍为MS7.8唐山地震所触发的余震．      

1976年唐山地震前的中源地震前兆研究

研究了唐山地震发生前的12年以上h≥33．1km地震的活动特征,得知唐山地震前确实存在中源地震前兆MDE;其最主要的特点是：大地震发生前几年,发生多次h≥60km的中源地震,特别是100km左右或大于100km的中源地震,是可能发生板内大地震的前兆信号。文中还研究了唐山MDE的一些具体特征。     

1976年唐山MS7.8地震同震及现今形变特征

基于1976年唐山M_S7.8地震的同震位移和现今GPS速度结果, 本文分析了该地震的同震破裂参数、 同震应变释放分布及其现今应变积累特征． 对唐山同震位移二维位错解析公式的拟合结果显示, 唐山地震前的断层闭锁深度约为18—23 km, 断层倾角可能介于74°—90°之间, 同震位错量约为3.1—3.4 m． 在同震过程中沿发震断层表现出显著的右旋错动特征, 在发震断层两侧沿NE向表现出左旋应变释放特征． 在本文现有的GPS测站密度和精度情况下, 尚无法识别出现今阶段沿唐山断裂的明显蠕滑特征． 唐山同震应变释放和现今GPS应变率积累结果均显示, 唐山断裂SE侧的剪切形变(速率)量值大于其NW侧, 同震与现今阶段的形变量值相差约1000倍．      

2020年7月12日唐山5.1级地震分析

对2020年7月12日唐山5.1级地震的发震特点、地震的性质、发震构造以及破裂机制进行初步分析,推测唐山断裂可能为其控震断裂.地震前唐山地区和震中所处的华北构造区的地震活动性异常以缺震和显著平静为主,表明该区域地壳应力积累到了一定程度.分析认为:此次唐山5.1级地震属于1976年7月28日唐山7.8级大震震区内的地震起...     

Analysis of foreshock sequence of the 1975 Haicheng earthquake of M S 7.3

The “earthquake nucleation” is discussed in this paper. The acceleration is a property of the nucleation phase and is a necessary condition of earthquake instability too. If the acceleration property of this nucleating process is described by the equation dΘ/dt=C/(t _f−t) ⁿ , the process can be summarized briefly that the rate of cumulative seismic release is proportional to the inverse power of the remaining time to failure. Based on this principle, the foreshock sequence of the 1975 Haicheng earthquake withM _S7.3, was analysed backward. It is stated clearly that the time-to-failure and magnitude of the mainshock can be predicted successfully if the coefficientr ² attains to the maximum. In the estimation of mainshock time, the error can generally be less than, or far less than, one-half the remaining time between the time of the last used data point and the mainshock. Contribution No. 95A0024, Institute of Geophysics, SSB, China. This study is the "Eighth-Five" contract project ofSSB.     

Tempo-spatial rupture process of the 1997 Mani, Xizang (Tibet), China earthquake of M S=7.9

An earthquake of M _S=7.4 occurred in Mani, Xizang (Tibet), China on November 8, 1997. The moment tensor of this earthquake was inverted using the long period body waveform data from China Digital Seismograph Network (CDSN). The apparent source time functions (ASTFs) were retrieved from P and S waves, respectively, using the deconvolution technique in frequency domain, and the tempo-spatial rupture process on the fault plane was imaged by inverting the azimuth dependent ASTFs from different stations. The result of the moment tensor inversion indicates that the P and T axes of earthquake-generating stress field were nearly horizontal, with the P axis in the NNE direction (29°), the T axis in the SEE direction (122°) and that the NEE-SWW striking nodal plane and NNW-SSE striking nodal plane are mainly left-lateral and right-lateral strike-slip, respectively; that this earthquake had a scalar seismic moment of 3.4×10²⁰ N·m, and a moment magnitude of M _W=7.6. Taking the aftershock distribution into account, we proposed that the earthquake rupture occurred in the fault plane with the strike of 250°, the dip of 88° and the rake of 19°. On the basis of the result of the moment tensor inversion, the theoretical seismograms were synthesized, and then the ASTFs were retrieved by deconvoving the synthetic seismograms from the observed seismograms. The ASTFs retrieved from the P and S waves of different stations identically suggested that this earthquake was of a simple time history, whose ASTF can be approximated with a sine function with the half period of about 10 s. Inverting the azimuth dependent ASTFs from P and S waveforms led to the image showing the tempo-spatial distribution of the rupture on the fault plane. From the "remembering" snap-shots, the rupture initiated at the western end of the fault, and then propagated eastward and downward, indicating an overall unilateral rupture. However, the slip distribution is non-uniform, being made up of three sub-areas, one in the western end, about 10 km deep ("western area"); another about 55 km away from the western end and about 35 km deep ("eastern area"); the third about 30 km away from the western end and around 40 km deep ("central area"). The total rupture area was around 70 km long and 60 km wide. From the "forgetting" snap-shots, the rupturing appeared quite complex, with the slip occurring in different position at different time, and the earthquake being of the characteristics of "healing pulse". Another point we have to stress is that the locations in which the rupture initiated and terminated were not where the main rupture took place. Eventually, the static slip distribution was calculated, and the largest slip values of the three sub-areas were 956 cm, 743 cm and 1 060 cm, for the western, eastern and central areas, respectively. From the slip distribution, the rupture mainly distributed in the fault about 70 km eastern to the epicenter; from the aftershock distribution, however, the aftershocks were very sparse in the west to the epicenter while densely clustered in the east to the epicenter. It indicated that the Mani M _S=7.9 earthquake was resulted from the nearly eastward extension of the NEE-SWW to nearly E-W striking fault in the northwestern Tibetan plateau. Contribution No. 99FE2016, Institute of Geophysics, China Seismological Bureau. This work is supported by SSTCC Climb Project 95-S-05 and NSFDYS 49725410.     

唐山M 4.8地震地磁谐波振幅比变化

针对2012年5月28日河北唐山M 4.8地震,采用地磁谐波振幅比方法,分析昌黎、静海、红山及承德地磁台2007年1月至2015年12月地磁观测数据,结果显示：各台站10 min周期的计算结果波动性明显,其他周期形态较为一致;距震中最近的昌黎地磁台Y_ZHx（NS）测向的10-40 min周期地磁谐波在下降至极低值后,开始转折时发生地震,而Y_ZHy（EW）测向未出现明显转折现象;位于河北平原地震带中部的静海地磁台对此次地震反应显著,2个测向在震前均出现下降-转折-上升的变化特征。通过对比位于不同地区的4个地磁台站地磁谐波振幅比变化,发现受大范围整体构造应力场作用的影响,计算结果具有相似的趋势性变化特征,具体到局部地区或断裂带内部,各台站之间显示出一定差异。     

Spatial and temporal rupture process of the January 26, 2001, Gujarat, India, M S=7.8 earthquake

The source parameters, such as moment tensor, focal mechanism, source time function (STF) and temporal-spatial rupture process, were obtained for the January 26, 2001, India, M _S=7.8 earthquake by inverting waveform data of 27 GDSN stations with epicentral distances less than 90°. Firstly, combining the moment tensor inversion, the spatial distribution of intensity, disaster and aftershocks and the orientation of the fault where the earthquake lies, the strike, dip and rake of the seismogenic fault were determined to be 92°, 58° and 62°, respectively. That is, this earthquake was a mainly thrust faulting with the strike of near west-east and the dipping direction to south. The seismic moment released was 3.5×10²⁰ Nm, accordingly, the moment magnitude M _W was calculated to be 7.6. And then, 27 P-STFs, 22 S-STFs and the averaged STFs of them were determined respectively using the technique of spectra division in frequency domain and the synthetic seismogram as Green’s functions. The analysis of the STFs suggested that the earthquake was a continuous event with the duration time of 19 s, starting rapidly and ending slowly. Finally, the temporal-spatial distribution of the slip on the fault plane was imaged from the obtained P-STFs and S-STFs using an time domain inversion technique. The maximum slip amplitude on the fault plane was about 7 m. The maximum stress drop was 30 MPa, and the average one over the whole rupture area was 7 MPa. The rupture area was about 85 km long in the strike direction and about 60 km wide in the down-dip direction, which, equally, was 51 km deep in the depth direction. The rupture propagated 50 km eastwards and 35 km westwards. The main portion of the rupture area, which has the slip amplitude greater than 0.5 m, was of the shape of an ellipse, its major axis oriented in the slip direction of the fault, which indicated that the rupture propagation direction was in accordance with the fault slip direction. This phenomenon is popular for strike-slip faulting, but rather rare for thrust faulting. The eastern portion of the rupture area above the initiation point was larger than the western portion below the initiation point, which was indicative of the asymmetrical rupture. In other words, the rupturing was kind of unilateral from west to east and from down to up. From the snapshots of the slip-rate variation with time and space, the slip rate reached the largest at the 4th second, that was 0.2 m/s, and the rupture in this period occurred only around the initiation point. At the 6th second, the rupture around the initiation point nearly stopped, and started moving outwards. The velocity of the westward rupture was smaller than that of the eastward rupture. Such rupture behavior like a circle mostly stopped near the 15th second. After the 16th second, only some patches of rupture distributed in the outer region. From the snapshots of the slip variation with time and space, the rupture started at the initiation point and propagated outwards. The main rupture on the area with the slip amplitude greater than 5 m extended unilaterally from west to east and from down to up between the 6th and the 10th seconds, and the western segment extended a bit westwards and downwards between the 11th and the 13th seconds. The whole process lasted about 19 s. The rupture velocity over the whole rupture process was estimated to be 3.3 km/s. Foundation item: 973 Project (G1998040705) from Ministry of Science and Technology, P. R. China, and the National Science Foundation of China under grant No.49904004. Contribution No. 02FE2026, Institute of Geophysics, China Seismological Bureau.     

A new technique for moment tensor inversion with applications to the 2008 Wenchuan M S8.0 earthquake sequence

A MS8.0 earthquake occurred in Wenchuan County, Sichuan Province, China, on May 12, 2008, and subsequently, numerous aftershocks followed. We obtained the moment tensor solutions and source time functions (STFs) for the Wenchuan earthquake and its seven larger aftershocks (MS5.0~6.0) by a new technique of moment tensor inversion using the broadband and long-period seismic waveform data from the Global Seismic Network (GSN). Firstly, the theoretical background and technical flow of the new technique was briefly introduced, and an aftershock of the Wenchuan earthquake sequence was employed to illustrate the real procedure for inverting the moment tensor; secondly, the moment tensor solutions and STFs of the eight events, including the main shock, were presented, and finally, the interpretation of the results was made. The agreement of our results with the GCMT results indicates the new approach is efficient and feasible. By using this approach, not only the moment tensor solution can be obtained but also the STF can be retrieved; the inverted STFs indicate that the source rupture process may be com-plicated even for the moderate earthquakes. The inverted focal mechanisms of the Wenchuan earthquake sequence show that the most of the aftershocks occurred in the main faults of the Longmenshan fault zone with predomi-nantly thrustingwith minor right-lateral strike-slip component, but some of them may have occurred in the sub-faults with strike-slip faulting in the vicinity of the main faults.     

Damage of 1975 Haicheng Ms7.3 earthquake in the Korean Peninsula and its effect on Kyushu, Japan

The 1975 Haicheng MS7.3 earthquake affected whole Korean Peninsula and even Kyushu,Japan.In this paperisoseismal map of this earthquake was drawn for these areas by sorting out and analyzing data systematically col-lected,furthermore,the characteristics of the earthquake damage were discussed.In the northern Korean Peninsula,the intensity is V in great majority of this area except local area with intensity VI,and the damage is basically inaccordance with typical characteristics of corresponding intensity grade.In the southern Korean Peninsula,theintensity is IV in large part of the area,however,in some big cities such as Seoul etc.,not only high buildingsshook strongly,but also some people appeared perceptions such as feeling dizzy and so on,and electrical servicewas interrupted on one or two districts because of transformer tripping.These phenomena could be caused by theeffect of long-period wave generated from a large earthquake on far-field.     

The middle-long term prediction of the February 3, 1996 Lijiang earthquake (MS=7) by the "criterion of activity in quiescence"

IntroductionIn the book Future CataS~ologr published in 1992, we proposed a viewpoiflt on using the"criterion of activity in quiescence" to predict big eathquake (MsZ7) (GUO, et al, 1992), and predicted in the book that in futore several years or in ten years a big earthquake (Ms27) will be possible to occur in the Zhongdian and nearby in Yunnan Province. In the 1994 nation-wide earthquake tendency consultation meeting we pointed out, once more, in the Zhongdian region of Yunnan Province…     

Seismicity anomalies before the great earth—quake of Ms=8.1 in the Kunlun Pass and its significance to earthquake prediction

A great earthquake of M _S=8.1 took place in the west of Kunlun Pass on November 14, 2001. The epicenter is located at 36.2°N and 90.9°E. The analysis shows that some main precursory seismic patterns appear before the great earthquake, e.g., seismic gap, seismic band, increased activity, seismicity quiet and swarm activity. The evolution of the seismic patterns before the earthquake of M _S=8.1 exhibits a course very similar to that found for earthquake cases with M _S≥7. The difference is that anomalous seismicity before the earthquake of M _S=8.1 involves in the larger area coverage and higher seismic magnitude. This provides an evidence for recognizing precursor and forecasting of very large earthquake. Finally, we review the rough prediction of the great earthquake and discuss some problems related to the prediction of great earthquakes.