断层破裂传播速度与破坏性地震的烈度分布——浅源地震的断层不均匀性与短周期波动的多普勒效应

地震的等震线不仅反映了地震的大小,而且也反映了地震断层过程的类型和破裂速度.单侧破裂的断层过程呈现出蛋型的等震线,而双侧破裂的断层过程呈现出椭圆的等震线,等震线的长轴与短轴之比对于破裂传播速度十分敏感.本研究利用理论等震线与观测结果相比较的方法,确定了1964年日本新地震、1983年日本海中部地震、1975年中国海城地震和1976年中国唐山地震的破裂传播速度、断层类型以及断层走向.得到的断层破裂速度是剪切波的0.7——0.9倍.这些值比用长周期地震波所确定的要稍大.产生这种差异的原因是:对烈度起主要影响的短周期地震波强烈地依赖于小尺度的断层不均匀破裂过程,以及局部的破裂传播的开始和终止;而由长周期波所得到的破裂速度却反映了在整个断层上破裂传播的平均过程.根据等震线图所得出的断层类型及断层走向与其它独立方法的结果相一致.这意味着本方法可以应用于推断某些历史地震的断层类型、破裂传播速度及破裂传播方向.      

震源机制与烈度分布的关系（二）

除了通常熟知的震源深度和介质吸收系数对烈度分布有较大影响外,震源的破裂方式（单侧或双侧）,破裂速度的大小和方向,断层长度,断层面倾角和错动方式都对烈度分布有较大影响。尤其是破裂速度的大小和方向对烈度分布的影响更大。破裂速度大,而且方向向上,可使震中烈度增高一度多,极震区面积扩大1-2倍;反之,可使震中裂度降低一度、极震区面积减小许多。 根据震源参数对烈度分布影响的分析,可由历史地震的等震线反推其震源参数。 文中试行反推了三个历史地震的震源参数。     

Fracture characteristics of the 1997 Jiashi, Xinjiang, China, earthquake swarm inferred from source spectra

Broadband P and S waves source spectra of 12 MS5.0 earthquakes of the 1997 Jiashi, Xinjiang, China, earthquake swarm recorded at 13 GDSN stations have been analyzed. Rupture size and static stress drop of these earthquakes have been estimated through measuring the corner frequency of the source spectra. Direction of rupture propagation of the earthquake faulting has also been inferred from the azimuthal variation of the corner frequency. The main results are as follows: ①The rupture size of MS6.0 strong earthquakes is in the range of 10～20 km, while that of MS=5.0～5.5 earthquakes is 6～10 km.② The static stress drop of the swarm earthquakes is rather low, being of the order of 0.1 MPa. This implies that the deformation release rate in the source region may be low. ③ Stress drop of the earthquakes appears to be proportional to their seismic moment, and also to be dependent on their focal mechanism. The stress drop of normal faulting earthquakes is usually lower than that of strike-slip type earthquakes. ④ For each MS6.0 earthquake there exists an apparent azimuthal variation of the corner frequencies. Azimuthally variation pattern of corner frequencies of different earthquakes shows that the source rupture pattern of the Jiashi earthquake swarm is complex and no uniform rupture expanding direction exists.     

Evolution of Stress Fields and Faulting in Seismic Zones

—Measurements indicate that stress magnitudes in the crust are normally limited by the frictional equilibrium on pre-existing, optimally oriented faults. Fault zones where these limitations are frequently reached are referred to as seismic zones. Fault zones in the crust concentrate stresses because their material properties are different from those of the host rock. Most fault zones are spatially relatively stable structures, however the associated seismicity in these zones is quite variable in space and time. Here we propose that this variability is attributable to stress-concentration zones that migrate and expand through the fault zone. We suggest that following a large earthquake and the associated stress relaxation, shear stresses of a magnitude sufficient to produce earthquakes occur only in those small parts of the seismic zone that, because of material properties and boundary conditions, encourage concentration of shear stress. During the earthquake cycle, the conditions for seismogenic fault slip migrate from these stress-concentration regions throughout the entire seismic zone. Thus, while the stress-concentration regions continue to produce small slips and small earthquakes throughout the seismic cycle, the conditions for slip and earthquakes are gradually reached in larger parts of, and eventually the whole, seismogenic layer of the seismic zone. Prior to the propagation of an earthquake fracture that gives rise to a large earthquake, the stress conditions in the zone along the whole potential rupture plane must be essentially similar. This follows because if they were not, then, on entering crustal parts where the state of stress was unfavourable to this type of faulting, the fault propagation would be arrested. The proposed necessary homogenisation of the stress field in a seismic zone as a precursor to large earthquakes implies that by monitoring the state of stress in a seismic zone, its large earthquakes may possibly be forecasted. We test the model on data from Iceland and demonstrate that it broadly explains the historical, as well as the current, patterns of seismogenic faulting in the South Iceland Seismic Zone.     

2003年青海德令哈6.7级地震序列的震源机制解及其构造含义

利用8个流动数字地震台和国家数字地震台站的地震波形记录,测量了2003年4月17日青海德令哈6.7级地震及其主要余震的直达P波、SV波、SH波的初动方向和振幅比,应用Snoke(2003)的测定震源机制解的格点尝试法,测定出德令哈地震序列的48个2.4级以上地震的震源机制解.搜集分析了美国哈佛大学测定的德令哈6.7级主震和2004年二期地震活动中的7个地震的震源机制解.基于余震空间分布特征和对震源机制解特征的分析,讨论了德令哈地震序列的可能断层活动方式和地震的构造含义.结果表明,主震和大部分余震都是沿NWW-SEE走向的逆断层错动,北边的上盘可能沿低角度向北倾的断层面向南仰冲;个别正断层余震可能是震源区挤压变形弧顶区附近发生的局部张性破裂;在二期地震活动中,逆断层和走滑断层都有,走滑断层地震主要发生在震源区东侧.德令哈地震活动是青藏高原东北缘NWW-SEE向延伸的挤压带继续处于隆升活动中的表现,这一继承性新构造运动是德令哈地震序列的可能发震原因.     

基于实时观测数据快速判定地震破裂方向的方法研究

针对地震烈度速报工作中快速评估大震断层破裂方向的需求,使用我国西部多个地震的观测资料,根据有限移动源理论,利用趋势面分析方法探讨了震中一定范围内台站初始P波参数的方向差异特征,提出一种利用初始P波快速判定地震破裂断层方向的方法。并通过断层距模型和最小二乘法,统计回归最优拟合参数以确定地震断层破裂方式。通过实际震例对比分析了根据所提方法获得的判定结果与实际地震断层破裂分布的吻合度,证实了方法的可用性。本研究为利用地震实时监测资料快速进行地震断层破裂方向判定提供了新思路,可为地震烈度速报中破裂方向的快速判定和动态修正提供参考。     

郯庐断裂带鲁苏皖段地壳速度结构的分段特征及其地质意义

本文采用天然地震近震走时反演地壳三维速度结构的方法获得了郯庐断裂带鲁苏皖段及附近地壳(30°N—37°N,113°E—122°E)三维速度结构.对地壳内分层速度结构的分析发现,郯庐断裂带鲁苏皖段存在速度的分段特征.郯庐断裂带鲁苏皖段浅层35.3°N以北,34.5°N—35.3°N间,33°N—34.5°N间呈现的速度分段和地表出露地层有关,与地质上安丘段、莒县—郯城段,新沂—泗洪段三个破裂单元相对应,且和各段的地震活动相呼应,表明郯庐带新沂到泗洪段可能是断裂的闭锁段.郯庐断裂带鲁苏皖段地壳速度结构自浅至深分为三段,大体位置是:南段(32.5°N—33°N以南),中段(32.5°N—33°N至35°N—35.3°N),北段(35°N—35.3°N以北).上地壳分段与苏鲁超高压变质岩带的插入有关,中、下地壳速度分段则可能和火山岩滞留有关.地壳各层速度结构不同段的速度差异反映了构造块体的速度差异,表明各构造块体在地壳下部仍有差异,郯庐带西侧速度总体高于东侧,反映了不同构造块体的形成和组成差别,也说明了该断裂带可能延伸到莫霍面.而不同深度的分段性可能反映了不同地质演化过程.     

Synthetic isoseismals of three earthquakes in California and Nevada

Recent tests on a series of earthquakes in California and Nevada suggest that in some regions the approximate shapes of the territories with equal earthquake-induced damage (expressed in terms of macroseismic intensity) could be synthetically traced out with a simple formula. This formula takes into account some gross features of the source: depth and length, unilateral or bilateral rupture, radiation patterns, rupture velocity, and directivity. Having been formulated on an empirical basis, the formula is however compatible with the so-called asymptotic approach, in which the far-field component of the Green's function is used. This paper presents the synthetic isoseismals of the earthquakes at Cedar Mountain, Nevada, 1932; Fairview Peak-Dixie Valley, Nevada, 1954; and Coalinga, California, 1983. An overall consistency, from acceptable to remarkable, between the observed intensity patterns and the synthetically back-predicted intensity has been obtained for them. Where the detailed modelling techniques available today are inapplicable, due to insufficient information on the features of the seismic sources, or to save time and money, the new formula may be utilizable for improving seismic hazard calculations.The formula was also used inversely for back-predicting geometric-kinematic parameters of the Coalinga 1983 earthquake from macroseismic maps. This gave characteristics for its source which are in good agreement with the majority of data inferred from modelling and from analyzing modern instrumental recordings. This striking result opens new perspectives in retrieving information on the source of ancient earthquakes for which only macroseismic information is available.     

2008年汶川大地震单侧破裂过程的动力学机制研究

2008年汶川大地震的破裂过程极其发杂,向东北方向的破裂距离长达300 km,而向西南方向的破裂长度很小,呈现出单侧破裂的主要特征.尽管汶川地震破裂呈单侧传播的现象引起许多地震学家的关注,但其物理机制至今还不是十分清楚.本文利用有限单元计算方法,模拟了汶川地震的破裂过程.模型中根据龙门山断裂带两侧（东南侧为四川盆地,西北侧为川西高原）实际的地震波速度来确定模型的介质物性参数,利用目前观测的应力环境来选定初始应力条件.模拟结果表明：破裂在汶川地震的震中处成核后,先向断层两侧自发传播,但向东北方向的传播距离明显大于向西南方向;断层面上的正应力在东北方向（破裂的正方向）随着传播距离的增大而不断减小,位错速率随着破裂的传播距离而越来越大,其脉冲变得越来越尖锐,即产生了Weertman脉冲.研究结果显示：由于这种脉冲的出现,破裂在正方向上（东北方向）能够自己放大、自己愈合、自行维持,摩擦热极小,所以破裂能够沿着东北方向一直传播,直到应力场方位发生变化,不利于破裂时才最后终止.但在西南方向,破裂过程中断层面上的正应力增大,阻碍破裂继续扩展.最后就出现了汶川地震中破裂朝东北方向单侧优势传播的基本格局.模拟结果还表明：若断层面两侧介质均匀,则破裂向两侧是对称传播,且破裂距离很短,因此这种情况无法产生像汶川大地震那样的特大地震.因此,文中的模拟结果表明龙门山断裂带两侧的物性差异是造成汶川大地震单侧传播的决定性因素.断层两侧物性差异（bimaterial contrast）影响断层破裂过程的研究对于深入认识地震动力学过程、地震灾害预测及评估等有重要的科学意义.     

1988年云南澜沧-耿马地震的烈度分布及地表破裂

1988年11月6日,在云南省西南部的澜沧-耿马断裂带上发生了两次大于7级的地震.地震造成的严重破坏和人员伤亡主要是由于极震区内抗震性能极差的毛石房、砖柱土坯房的大量倒塌所致.澜沧地震的震中基本烈度可达Ⅸ度,耿马地震极震区烈度达Ⅹ度.澜沧地震构造活动的地表证据主要是出现在极震区内的张性地裂缝带和小断层陡坎.地裂缝带和小断坎主要出现在四条相对连续的北北西走向的狭窄地带内,其长度从几百 m 到6km 不等.澜沧地震地表破裂带长约35km,宽约3km,最大垂直位移量和最大右旋水平位移量分别为1.5m 和1.4m.耿马地震地表断裂活动的明显证据包括一系列北北西走向的地裂缝带和一段长约5km的地震断层陡坎.耿马地震地表破裂带长约24km,其最大垂直位移3.5m,最大右旋水平位移为3m.两次地震形成的地表破裂带均具有右旋-正断层性质.本文讨论了这两次地震的 度分布及地表破裂特征.      

Dynamic faulting under rate-dependent friction

We discuss the effects of rate-dependent friction on the propagation of seismic rupture on active faults. Several physicists using Burridge and Knopoff's box and spring model of faulting have proposed that fault complexity may arise from the spontaneous development of a self-similar stress distribution on the fault plane. If this model proves to be correct, it has important consequences for the origin of the complexity of seismic sources. In order to test these ideas on a more realistic earthquake model, we developed a new boundary integral equation method for studying rupture propagation along an antiplane fault in the presence of nonlinear rate-dependent friction. We study rupture dynamics of models with single and twin asperities. In our models, asperities are places on the fault with a higher value of prestress. Othewise all fault parameters are homogeneous. We show that for models with such asperities, a slip velocity weakening friction leads to the propagation of supersonic healing phases and to the spontaneous arrest of fracture if the prestress outside the asperities is low enough. For models with asperities, we can also observe narrow slip velocity pulses, qualitatively similar to the so-called Heaton pulses observed in some earthquake accelerograms. We also observe a complex distribution of stress after the rupture that depends on details of the initial distribution of asperities and on the details of the friction law.     

2012年江苏高邮、宝应交界MS4.9地震发震断层参数测定

基于一维单侧有限移动震源模式,根据地震波传播过程中的多普勒效应,分别利用P波和S波拐角频率的方位变化,反演2012年7月20日江苏高邮、宝应交界MS4.9地震的发震断层面参数。P波和S波拐角频率的反演结果一致显示:本次地震的断层面破裂方向为232°左右,破裂面呈NE-SW向;地震马赫数v/c为0.2左右,平均破裂速度小于S波速度,破裂长度较短,为0.2~0.3km左右。破裂面方位与震源机制解、宏观烈度调查和余震精定位的研究结果具有一致性,结合震区周边的地质构造背景,分析认为滁河断裂很可能是高邮、宝应交界MS4.9地震的发震构造。     

2008年汶川8.0级地震发震断裂的滑动速率、复发周期和构造成因

2008年5.12汶川大地震发生在中国大陆南北地震带中段.由于龙门山断裂带历史上只发生过3次6～6¹/₂级强震,而且其晚第四纪构造活动速率很低,以至于对其潜在地震危险性认识不足.为什么在龙门山地区突发大地震,该地震具有哪些特征？其成因机制是什么？本文在地震地质科学考察的基础上,利用震前的GPS观测结果,试图对上述问题进行一些初步的思考和探讨.结果表明,5.12汶川大地震是龙门山断裂带的映秀—北川断裂突发错动的结果,地表上形成200多公里长的地表破裂带;灌县—江油断裂在地震中也发生了破裂,形成的地表破裂带长达60多公里.震前的GPS观测表明,横跨整个龙门山断裂带的滑动速率不超过～2 mm/yr,单条断裂的活动速率不超过～1 mm/yr,与地震地质研究结果和历史地震记录相一致.利用地震地质考察和地震波反演得到的最大同震位移可以获得相当于5.12汶川大地震的强震复发周期为2000～6000年.龙门山断裂带发育在破裂强度很大的变质杂岩体中,断裂带本身在剖面上呈“犁形”或“铲形”结构,有利于能量积累,形成破坏性巨大的地震.所以,5.12汶川大地震是一次低滑动速率、长复发周期和高破坏强度的巨大地震,是一种值得高度重视和深入研究的新的地震类型.     

Nonuniform and Unsteady Sliding of a Plate Boundary in a Great Earthquake Cycle: A Numerical Simulation Using a Laboratory-derived Friction Law

—A numerical study is conducted to simulate complicated sliding behavior and earthquake activity on a subducting plate boundary. A 2-D model of a uniform elastic half-space with a semi-infinite thrust fault is set up, and the frictional stress prescribed by a rate- and state-dependent friction law is assumed to act on the plate boundary fault. Spatial nonuniformity of friction parameters representing rate-dependence of friction and of slip-dependence of friction are introduced in the model to obtain complicated sliding behavior in the numerical simulation. Analogs of great earthquakes that break the entire seismogenic plate boundary repeatedly occur at a constant time interval. Smaller events of seismic or aseismic sliding occur during a great earthquake cycle. Regions of rate-strengthening of friction and of a large characteristic distance in slip-dependence of friction behave as barriers or asperities. Rupture propagation is often arrested in such a region and a great earthquake occurs later when the region is broken. The variety of earthquake activity observed in many regions along real plate boundaries may be explained by similar nonuniformity in friction parameters. Conversely, the friction parameters on plate boundaries might be estimated from comparison of theoretical simulations with observations of earthquake activity. Simulation results indicate that spatiotemporal variation in stress due to aseismic sliding may play an important part in generating earthquakes.     

晋冀蒙交界地区地震活动性特征分析

采用较为完整的地震资料,分析晋冀蒙交界地区地震活动性特征,结果表明,公元1929年以后,研究区的6.0级以上地震主要分布在河套地震带和张家口-渤海地震带上,在北纬40°线及附近区域形成优势分布区域,1970年以来的中等以上地震主要分布在内蒙古和林格尔-山西大同-河北张北形成的三角形区域内。研究区2005年9月以来中等地震持续1 304天的平静现象被2009年3月山西原平4.2级地震打破,岱海断陷1997年11月以来至今中等地震持续平静,统计结果显示,研究区和岱海断陷中等地震显著平静对研究区特别是河套地震带中强地震具有明显的指示意义,为显著的震兆标志,多种迹象表明,晋冀蒙交界三角形区域和呼和浩特-包头断陷为有利发震区域,发震强度或为6.0级左右。     

地震断裂和应力场

本文讨论了地震断裂和构造应力之间的关系。首先综述了地震机制研究的结果,指出断层面解可以给出较好的断层面的取向,而得不到确切的应力场;认为地震是一个复杂的破裂过程,而P波初动只反映了破裂开始的状态。应用了库仑准则对地震断裂进行了讨论,得出在各向同性介质中,最大主应力与断层面的夹角α≤45°,断层面解只是它的一个特例,μ=0,即没有内摩擦,或相当于围压很高的岩石试验结果;在软弱面存在的情况下,断裂面与最大主应力的夹角可以分布在更大的范围内,可以大于45°,也可以小于45°,这范围取决于破裂参数S_w,μ_w以及中等主应力σ_2。讨论了震前应力变化的过程,指出了五种基本可能引起地震的过程     

利用小震分布和区域应力场确定大震断层面参数方法及其在唐山地震序列中的应用

根据成丛小震发生在大震断层面及其附近的原则,将模拟退火算法和高斯-牛顿算法结合,给出了利用小震密集程度求解主震断层面走向、倾角、位置及其误差的稳健估计方法,在此基础上考虑区域构造应力参数,给出了估计在已求得的断层面上的滑动角的方法.该方法还可用于小震活跃地区活断层走向、倾角和滑动角的确定.将这种方法用于唐山地震序列,采用2002年4月1日至2006年5月31日发生在地震破裂区的精定位地震目录,求得了唐山地震、滦县地震、宁河及卢龙断裂带的断层面走向、倾角、位置及滑动角参数.与前人给出的断层面解进行比较,发现利用小震精定位资料和区域构造应力场得到的结果与前人采用其他资料和方法得到的结果近似,验证了这种方法的有效性.另外,本研究首次发现滦县地震区东部的小震呈北东-南西向条带状成丛发生,可精确刻画为一条断裂带,较为精确地确定了此断层的走向、倾角和滑动角.该断裂及宁河断裂在唐山地震序列发生时是否破裂需要运用其他资料进行验证.     

Source parameters of central California earthquakes on February 24 and September 4, 1972

The directivity function method is combined with an earthquake sequence study to obtain the reliable estimates of rupture area and rupture velocity of two right-lateral strike-slip earthquakes occurring along the San Andreas fault zone in central California. By utilizing a significant difference in velocity structure on both sides of the main fault, a modified version of the directivity function is formulated and applied to near-field SH waves. On assuming aftershock areas following their step-wise increase with time as rupture areas, an extensive systematic search for the best fit between the observed and theoretical directivity functions is made for a combination of both source parameters. The rupture at the February 24, 1972 earthquake (M_L = 5.0) propagated unilaterally with an average velocity of about 2.3 km/sec. It produced a rectangular area, the horizontal and vertical lengths being about 5 and 2 km, respectively. The rupture at the September 4, 1972 earthquake (M_L = 4.6) was of bilateral, yielding a nearly square rupture area of which side length is about 2 km. The rupture velocity of this earthquake, though some ambiguities resulting from a lower quality of the observed directivity function, is estimated at 1.9 km/sec or less. A difference in average rupture velocity between both earthquakes might imply its dependence on the ambient tectonic environment such as represented by local stress and humidity. By taking into account a post-earthquake creep increase of about 3.0–3.5 cm observed at a creepmeter station situated over the focus of the February 24, 1972 earthquake, its stress drop and seismic moment are estimated at about 10 bar and 10²³ dyn · cm, respectively. The above procedure has a broad applicability for recovering reliable estimates of source parameters, especially when it is combined with a synthetic approach.     

2013年岷漳地震和2014年景谷地震的能矩比对比研究

2013年岷漳地震和2014年景谷地震两个地震的矩震级均为M_W6.1,但相同等级烈度区域的面积差异显著,前者均明显大于后者.地面震动的强烈程度与震源破裂过程中辐射的地震波高频成分密切相关,作者系统地测定与研究了这两个地震的震源谱和反映其震源特征的能矩比等震源参量.结果表明,虽然岷漳地震和景谷地震释放的地震矩基本相同,分别为1.6×10¹⁸Nm和1.8×10¹⁸Nm,而且前者还比后者略小些,但其释放的地震波辐射能却差异显著,分别为1.3×10¹⁴J和0.6×10¹⁴J;因而前者的能矩比（8.1×10^-5）是后者能矩比（3.3×10^-5）的大约2.5倍.相应地,岷漳地震的视应力和应力降也都明显大于景谷地震的视应力和应力降.此外,这两个地震震源谱结果均显示出其高频趋势均按角频率ω的-2幂次衰减,符合"ω²模式";采用理论震源谱公式计算结果分别约占实际结果的54%和50%,两个地震的地震波辐射能解析解之间的差异依然约2.3倍.震源破裂特征的显著差异导致岷漳地震单位地震矩释放的地震波辐射能更多,即岷漳地震能矩比更大,是导致这两个矩震级均为M_W6.1地震的地震烈度差异显著的一个重要原因.     

Simulation of Ground Motion from Strong Earthquakes of Kamchatka Region (1992–1993) at Rock and Soil Sites

To estimate the parameters of ground motion in future strong earthquakes, characteristics of radiation and propagation of seismic waves in the Kamchatka region were studied. Regional parameters of radiation and propagation of seismic waves were estimated by comparing simulations of earthquake records with data recorded by stations of the Kamchatka Strong Motion Network. Acceleration time histories of strong earthquakes (M _w = 6.8–7.5, depths 45–55 km) that occurred near the eastern coast of Kamchatka in 1992–1993 were simulated at rock and soil stations located at epicentral distances of 67–195 km. In these calculations, the source spectra and the estimates of frequency-dependent attenuation and geometrical spreading obtained earlier for Kamchatka were used. The local seismic-wave amplification was estimated based on shallow geophysical site investigations and deep crustal seismic explorations, and parameters defining the shapes of the waveforms, the duration, etc. were selected, showing the best-fit to the observations. The estimated parameters of radiation and propagation of seismic waves describe all the studied earthquakes well. Based on the waveforms of the acceleration time histories, models of slip distribution over the fault planes were constructed for the studied earthquakes. Station PET can be considered as a reference rock station having the minimum site effects. The intensity of ground motion at the other studied stations was higher than at PET due to the soil response or other effects, primarily topographic ones. At soil stations INS, AER, and DCH the parameters of soil profiles (homogeneous pyroclastic deposits) were estimated, and nonlinear models of their behavior in the strong motion were constructed. The obtained parameters of radiation and propagation of seismic waves and models of soil behavior can be used for forecasting ground motion in future strong earthquakes in Kamchatka.