地球物理学报  2012, Vol. 55 Issue (4): 1213-1227   PDF    
汶川MS8.0地震余震震源机制时空分布特征
易桂喜1,2 , 龙锋2 , 张致伟2     
1. 成都理工大学地球探测与信息技术教育部重点实验室, 成都 610059;
2. 四川省地震局, 成都 610041
摘要: 本文利用CAP波形反演方法,获取了汶川MS8.0地震序列中312个具有较高信噪比波形资料的4级以上余震的震源机制解和震源深度. 基于震源深度空间分布与震源机制时空分布,分析了主震后余震区断层行为特征与应力场时空变化,并对龙门山断裂带中北段的发震断层面几何形态进行了初步探讨. 获得的主要认识如下:(1)余震震源深度分布存在显著的空间分段差异. 绵竹以西的余震区南段与平武以东的北段余震深度范围大于中段(绵竹—平武段),但深度小于5 km的5级以上超浅源地震主要分布在明显偏离龙门山断裂带走向的理县NW向分支与余震区北端NNE向分支,而中段余震主要分布在7~19 km深度. (2)余震机制类型存在明显的时空差异. 余震区中段逆冲型地震占绝对优势,理县NW向分支余震则以走滑型为主,机制类型随时间变化不显著. 沿龙门山断裂带走向的余震区南段,早期(2008年8月底前)逆冲型地震比例高于走滑型、晚期走滑型地震比例显著升高并超过逆冲型;而余震区北段早期走滑型地震占绝对优势、晚期逆冲型地震比例大幅上升且超过走滑型. 南、北两段余震机制类型比例的显著变化,可能是余震区两端断层调整性运动的表现. (3)节面走向及P轴方位优势方向均存在显著的空间差异. 南段NWW向P轴方位与区域应力场一致,中段及理县NW向分支P轴优势方向NEE,而北段具NWW和NEE两个优势方向,这种差异反映了余震活动除了受区域应力场控制外,还受到主震引发的局部应力场的控制. 节面走向的多方位分布则反映不同走向的构造参与了主震后的余震活动. (4)沿龙门山断裂带走向,余震区南段具深部缓倾角、浅部高倾角的铲形断面特征;中段深部倾角均值较稳定、浅部倾角均值随深度减小而增大;北段倾角均值相对稳定,显示其断面几何形态相对简单. 上述不同区段倾角均值随深度的变化揭示龙门山断裂带中北段断层面几何形态复杂.
关键词: 汶川地震      余震      震源机制      时空变化     
Spatial and temporal variation of focal mechanisms for aftershocks of the 2008 MS8.0 Wenchuan earthquake
YI Gui-Xi1,2, LONG Feng2, ZHANG Zhi-Wei2     
1. Key Lab. of Earth Exploration and Information Technique of Education Ministry of China, Chengdu University of Technology, Chengdu 610059, China;
2. Earthquake Administration of Sichuan Province, Chengdu 610041, China
Abstract: Focal mechanism solutions and depths of 312 M≥4 aftershocks with relatively high signal-to-noise ratio waveforms of the 2008 Wenchuan MS8.0 earthquake have been derived by CAP method. Following this, we analyzed the aftershock faulting behavior and spatial-temporal variation of stress field, and further discussed preliminarily the geometry of the central-northern segment of the Longmen Shan fault zone. The major results are as follows. (1) The range of focal depths is significantly different between various segments of aftershock area. The depth ranges of aftershocks on the southern segment west of Mianzhu and northern segment east of Pingwu are greater than that on the central segment from Mianzhu to Pingwu. Relatively shallow (<5 km) aftershocks mainly occurred along NW-striking Lixian branch on the southern segment and NNE-striking branch at the northern end of the aftershock area. On the central segment, aftershocks are concentrated at the depths of 7~19 km. (2) The proportion of different-type aftershocks varies with time on different segments, especially on the southern and northern segments. Thrust slip and strike-slip motions are dominant faulting behavior for the central segment and NW-striking Lixian branch, respectively. However, on the southern segment along the Longmen Shan fault zone, thrust-slip aftershocks are more than those of strike-slip in the early stage before Aug. 31, 2008, and slightly less than those of strike-slip in the later stage. Meanwhile, for the northern segment, the proportion of aftershocks of right-lateral strike-slip type is higher than that of reverse type with minor strike-slip component in the early stage, but it is less than that of reverse type in the later stage. The significant variation of the proportion of different-type aftershocks on the southern and northern segments may be resulted from readjustment of faulting behavior on both ends of the seismogenic fault zone. (3) Variation of dominant strike of fault planes and orientation of P-axes is significant on different segments as rose diagrams indicate. The dominant orientations of P-axes are NWW with sub-horizontal plunges on the southern segment, NEE on the central segment and NW-striking Lixian branch, and NWW and NEE on the northern segment. The above finding demonstrates that the aftershock activity was influenced by both regional stress field and local stress field from the Wenchuan mainshock. (4) Along the central-northern Longmen Shan fault zone, the mean dip of the fault plane on the southern segment of the aftershock area gradually increases as the depth decreases from 18 to 10 km, and becomes relatively stable for depths shallower than 10 km, which means that the fault plane changes from gentle to steep upwards as a listric fault plane. While on the central segment, the mean dips at depths >10 km are relatively steady, and obviously increase as depths decrease upwards (<10 km). However, the relatively stable mean dips within the range of 50°~65° on the northern segment imply that the fault plane there is relatively steep and therefore more favorable for strike-slip faulting. The above results indicate the very complicated geometry of the fault plane of the central-northern Longmen Shan fault zone..
Key words: Wenchuan earthquake      Aftershocks      Focal mechanism      Spatial-temporal variation     
1 引 言

龙门山断裂带位于青藏高原东缘、南北地震带中段,主要由近于平行的茂汶-汶川断裂(后山断裂)、映秀-北川断裂(中央断裂)、灌县-彭县断裂(前山断裂)等多条主干断裂组成[1-11],宽约30~40km,长约500km,总体走向NE(N40°-50°E),倾向北西,倾角不定[4-11].2008 年5 月12 日汶川MS8.0地震就发生在龙门山断裂带的中-北段,地震破裂沿龙门山断裂带中央断裂呈NE 向单侧扩展,同震破裂还局部扩展到龙门山断裂带前山断裂中段,形成总长约235km 的地表破裂带[12],而余震密集分布带则长达近330km,强余震空间分布表现出显着的空间差异[13].汶川MS8.0 地震被认为是巴颜喀拉块体朝东南东方向的挤出运动受到华南地块阻挡、沿巴颜喀拉块体与四川盆地之间的龙门山断裂带长期积累应力应变并突然释放的结果[14-19].

震源机制解可以直观反映地震破裂的几何特征与运动学特征[20].地震序列震源机制与震源深度分布对了解震源区应力状态、勾画断层形态、厘清主余震与发震构造的关系,理解主余震孕震机理等具有重要意义[20-21].汶川地震后,不同学者通过地震资料[22-28]、GPS 与InSAR 资料[29-30]反演主震破裂过程与主震震源机制,获得主震发震断层为龙门山断裂带中央断裂-映秀-北川断裂的共识.与此同时,关于汶川余震震源机制及震源区应力场的相关研究也陆续展开[31-35],但研究主要集中于序列早期,所得结果仅能反映余震区早期的构造活动与应力场特征.

截止到2011年6月,由四川区域台网记录的汶川余震近8.8 万次,其中M4 级以上余震超过400次.丰富的余震为深入研究汶川8.0级地震后余震区断层行为与震源区应力场分布特征提供了资料保障.本文采用Zhao 等[36]提出、并经Zhu等[37]发展的CAP(CutandPaste)波形反演方法,利用2008年5月至2011年6 月四川台网记录的宽频带地震波形资料,同时求取汶川地震序列中M≥4.0 余震的震源机制解与震源深度,试图根据余震震源机制的时空分布,探讨余震区断层行为与应力场的时空变化特征,并基于余震震源机制与震源深度结果,对龙门山断裂带中北段断层面展布特征进行初步分析.

2 资料与方法 2.1 资 料

地震波形资料来源于四川区域地震台网.为了保证研究结果的可靠性,本文依据下列原则选取波形资料:(1)余震震级M≥4.0;(2)记录台站震中距小于250km、波形连续且信噪比较高.经过筛选,汶川地震序列中,满足上述条件的余震共312 个.图 1给出了所用地震的震中(圆点)与台站(三角形)分布,其中,黑色三角形代表波形反演所用台站,蓝边灰色三角形为求取P波初动解所用补充台站.为便于后续对比分析,在该图中还同时标示出了由中国地震台网测定的汶川主震震中位置、震源机制[22]及地表破裂分布[12].

图 1 区域主要活动断裂、所用余震震中与台站分布、汶川主震震中与震源机制及地表破裂展布 Fig. 1 Tectonic units, distribution of epicenters of aftershocks and stations used, epicenter and focal mechanism as well as surface rupture of Wenchuan mainshock
2.2 计算方法

本文利用CAP 波形反演方法求取余震震源机制及震源深度,具体方法参见相关文献[20-2136-40].该方法将宽频带数字地震波形记录分为体波Pnl与面波两部分,分别对Pnl波、面波进行带通滤波,并计算其理论地震波形与观测波形之间的误差函数,通过网格搜索,获取给定参数空间中误差函数达到最小的最佳解.

与其他求解震源机制解的方法相比,CAP 方法具有所需台站少、反演结果对地壳横向变化不敏感、对速度模型依赖性相对较小等优点[20-2138-40],从而保证了结果的稳定性,该方法尤其适用于横向速度变化较大的龙门山构造带.

本文计算时,体波Pnl与面波截取波形窗长分别设置为30s与60s,其带通滤波频带宽度分别为0.05~0.2Hz、0.05~0.1 Hz.计算理论地震图时所使用的速度模型为郑勇等[20]的龙门山地区分层地壳速度结构(见图 2中模型1).

图 2 反演所用速度模型(模型1)及测试模型(模型2) Fig. 2 Velocity model 1 for inversion and model 2 for test

图 3给出了2011 年6 月5 日茂县M4.3 地震波形拟合实例.不同台站的理论地震波形与实际地震波形拟合(图 3a)结果显示:观测波形与理论地震图的相位和振幅均拟合较好.图 3b 则给出了不同震源深度搜索的机制解所对应的拟合误差,当震源深度为15km 时,拟合误差达到最小,即为最佳解,此处的震源机制解显示发震断层为逆冲型,断层走向θ 、倾角δ 与滑动角λ 分别为348°、59°、102°.

图 3 2011年6月5日茂县4.3级地震波形拟合与震源机制 (a)理论地震图与观测地震图拟合.红线表示理论地震图,黑线表示观测地震图,波形下方的两行数字分别表示理论地震图相对观测地震图的移动时间及二者的相关系数(用百分比表示);b)反演误差随深度的变化.给出不同震源深度下搜索的震源机制解所对应的误差. 当深度为15 km时,误差最小,即为最佳解. Fig. 3 Waveform fitting and focal mechanism for Maoxian 4. 3 earthquake on Jun. 5,2011 (a) Comparison of synthetic waveforms with observations. Red and black lines are synthetics and observations, respectively. Numbers below the traces are the time shifts of the synthetics relative to the observations and the corresponding cross-correlation coefficients (n percentage) ; (b) Fitting error as a function of focal depth. Itshows the fitting error of focal mechanism ateach searched depth. The error i s the m i n i mum and the focal mechan i sm result i s the best at 15 km depth.
2.3 结果可靠性分析

罗艳等[21]的研究表明,速度结构差异10%所引起的深度误差小于15%,而10km 左右的水平定位误差对震源深度结果的影响则可以忽略.汶川余震区位于台网监测能力较强的区域,4 级以上余震水平定位误差均小于5km,因此,水平定位误差对震源深度的影响可以忽略.本文利用不同速度模型(见图 2)进行反演结果的稳定性测试,结果显示不同模型间深度差异小于2km,与郑勇等[20]的研究结果一致.

为了确定震源机制解结果的可靠性,本文以2008年8月1日青川6.1级、2009年10月20日北川4.9级、2010年5月25日汶川5.0级、2010年11月15日青川4.5级与2011年6月5日茂县4.3级地震为例,分别给出了利用CAP 波形反演和P 波初动解求取的震源机制解,结果列于表 1中,可以看出,2 种方法获得的结果尽管在具体数字上存在少许差异,但获得的震源力学性质是一致的.

表 1 不同方法获得的5个余震震源机制解 Table 1 Focal mechanisms derived from different methods for five aftershocks

同时,我们将利用CAP方法获得的震源机制解与能够收集到的35次由远震波形得到的汶川M5.0以上地震的HarvardCMT 结果进行对比,显示两组结果给出的所有地震震源力学性质一致.表 1列出了其中2 个地震的Harvard CMT 结果(节面解).由表 1图 4不难发现,利用不同资料与不同方法获得的结果极为接近,表明本文结果是可信的.

图 4 表 1中所列5次余震不同方法获得的震源机制解对比 Fig. 4 Comparison of focal mechanisms derived from different methods for five aftershocks listed in Table 1
3 结果分析与讨论 3.1 余震深度分布特征

利用CAP 波形反演获得的312个余震沿龙门山断裂带走向的震源深度分布见图 5a.图像显示,余震分布在1~19km 深度区间,最大深度与利用深度震相获得的定位结果一致[41],深度范围小于双差定位结果[42-45],但深度优势分布区间4~17km,与已有研究结果[41-45]相近.

图 5 余震震源深度分布 (a)全部地震;(b)剔除南段理县NW向分支与北端NNE向分支后的地震. Fig. 5 Distribution of focal depths of aftershocks (a) All aftershocks; (b) Aftershocks excepting those on NW-striking Lixian branch of southern segment and NNE-striking branch at the northern end of the aftershock area.

图 5a表明,余震震源深度分布存在显着的空间分段差异:位于余震区中段的龙门山断裂带绵竹至平武段,4级以上余震主要分布在7~19km 深度,鲜有深度小于7km 的4 级余震活动;而绵竹以西与平武以东的段落,余震深度分布范围较大,尤其是青川以东段落,具有深度小于5km 的超浅源6 级余震活动[21].图 5b为剔除了图 1中余震区南段与北端明显偏离龙门山断裂带走向的余震后的震源深度分布.对比图 5a图 5b,不难发现,深度小于5km的5级以上浅源地震主要分布在余震区南段理县NW 向分支以及余震区北端与平武-青川断裂斜交的NNE 向分支(图 1),这些偏离龙门山断裂带且明显呈带状分布的超浅源地震揭示了不同走向分支断层的存在.而中段大部分区域6 级余震空缺[13],可能代表了最大位移附近的应力释放较彻底,后期应力调整幅度相对较小,因而余震强度(震级)低于两侧.

3.2 震源机制时空分布特征

本文利用滑动角判别余震震源机制类型(即发震断层的错动性质),判别标准如下:45°~135° 为逆冲型;-135°~-45°为正断型;其余为走滑型.统计结果显示,汶川余震以逆冲型与走滑型为主,312次地震中,逆冲型、走滑型、正断型地震分别为143、149和20次,各占总数的45.8%、47.8%、6.4%.

为了更好地分析余震震源机制类型时空分布特征,根据余震空间分布特征,划分出南段明显偏离龙门山断裂带走向的理县NW 向分支(图 6中黑色椭圆圈定区域),并分别以绵竹、平武为界,将近330km长的余震密集区沿龙门山断裂带划分为南、中、北3段(见图 6,黑色虚线标示分段界限).同时,根据6级以上强余震集中发生在2008年5-8月的时间活动特征,以2008年8月31日为界,将统计时间分成早期与晚期2个统计时段.表 2列出了不同时段各段落不同机制类型余震的数量与比例,统计结果表明,机制类型比例存在显着的时空差异.

图 6 312个余震的震源机制分布(黑色虚线与椭圆标示分段边界) Fig. 6 Focal mechanisms of 312 aftershocks
表 2 不同区域、不同时段余震类型统计 Table 2 Types of aftershocks on different for two different periods

(1)理县NW 向分支该区域系余震区南段理县附近明显偏离龙门山断裂带的NW 向余震分支(见图 1图 6).参与计算的余震32次(表 2),其中,早期余震28次,逆冲、走滑、正断型地震分别为3、23 和2 次;晚期余震4次,走滑型3 次,正断型1 次,无逆冲型地震活动.可见,无论在余震活动早期还是晚期,该区域震源机制类型均以走滑型占绝对优势,与已有研究一致[2031-32].

余震节面走向统计显示,节面主要集中在NNE与NWW 两个优势方向(图 7a).根据余震的空间分布,可判断NWW 走向的节面应为其发震断层面;而滑动角与倾角则显示发震断层以高角度的左旋走滑运动为主,兼具一定的正断分量.至于该发震断层是先存的NW 向米亚罗断裂[20]因汶川主震所引起的应力调整分配的表现还是同震过程所产生的新生断层,尚需进一步研究.

图 7 理县NW向分支节面走向、滑动角、倾角(a)与P轴方位角、仰角统计(b) Fig. 7 Strikes, rakes and dips of fault planes (a) , azimuths and plunges of P axes (b) for the Lixian NW-striking branch

与沿龙门山断裂带展布的3 个段落相比,该NW 向分支的P轴仰角略大,均值近30°(图 7b).P轴方位集中在NE-NEE 向,明显有别于NW-NWW 向的区域应力场方向[46-51]和汶川主震最佳双力偶解P轴方向[22],但与NE 向主震破裂传播方向[23-25]接近.大地震可能引起震源区主应力方向发生变化[52-54].由此推测,该区域的余震应是汶川主震牵动了高角度理县NW 向分支断裂活动的结果.同时,考虑断裂节面(约290°~300°)与P轴方位(40°~80°)斜交,该断裂的运动性质应为左旋走滑,这与根据滑动角和倾角特征表现出发震断层为高角度的左旋走滑断裂的推论一致.

(2)余震区南段绵竹以西沿龙门山断裂带走向的余震区南段长约100km,参与统计的地震共116 次,其中,逆冲型、走滑型、正断型分别为59、48、9 次(表 2),显示南段余震以逆冲与走滑型为主.但在余震活动的不同时段,不同机制类型地震比例具有显着变化:早期91次余震中,逆冲型地震所占比例为51.6%,走滑型地震38.5%,逆冲型地震比例明显高于走滑型,与已有研究结果相吻合[2031-32].晚期25次地震中,逆冲型地震比例为48%,走滑型52%,且无正断型地震活动.与早期相比,晚期走滑型地震比例升幅达13.5%.这种走滑型地震比例大幅上升、并超过逆冲型地震比例的特征,显示出该段落断层运动特性的调整变化.汶川主震后的较短时期内,余震特征沿袭了主震的特征,表现为逆冲作用为主[55].而在后期震后调整过程中,走滑型地震比例增大,以此来补充和协调断裂带沿走向方向上的应力变化.此外,南段早期接近整个余震区一半数量的正断型地震,可能是在大规模岩体逆冲运动过程中,由重力作用引起的岩体正断滑移调整运动所致[56].

南段节面优势方向较单一(图 8a),余震早期节面优势方向呈NEE-NE 向,与龙门山断裂带走向大体一致;在余震活动晚期,节面优势方向除NE 向外,可看出近NS 向的节面明显增多.绝大部分近90°的滑动角表明断层性质以逆冲为主,断层面倾角在30°~90°范围,反映了断层面空间形态的复杂性.

图 8 南段节面走向、滑动角、倾角(a)与P轴方位角、仰角统计(b) Fig. 8 Strikes, rakes and dips of fault planes (a) , azimuths and plunges of P axes (b) on the southern segment

图 8b显示,该段P轴仰角主要集中在30°以内,小于理县NW 向分支,尤其是在余震晚期,P轴仰角更接近水平,与震前区域主压应力近水平的特征一致[46-48].2个时段的结果均显示该段P轴方位角优势方向单一,为NWW 向,与区域应力场方向[46-51]一致,也与汶川主震P轴方位[22]以及张勇等[24]通过波形反演获得的主震最大子事件P轴方位基本吻合,表明南段的余震断层活动主要受区域应力场控制.

(3)余震区中段

位于绵竹与平武之间的余震区中段长约120km,参与计算的地震共68次(表 2图 9).其中,早期余震44 次,逆冲型、走滑型、正断型所占比例分别为61.4%、29.5%、9.1%;晚期余震24 次,逆冲型、走滑型、正断型比例各为66.6%、29.2%、4.2%,晚期逆冲型地震比例小幅上升,而正断型地震比例有所下降;可见,中段余震震源机制类型以逆冲型占绝对优势(图 9),机制类型比例随时间变化特征不如南段显着.

图 9 余震区中段不同时期的余震震源机制分布 Fig. 9 Focal mechanisms of aftershocks on the central segment of the aftershock area during different stages

该段节面走向优势方向在不同时段存在明显差异(图 10a).早期节面走向以NE-NNE 和NNW-NW 为主,其中NNW-NW 略占优势(图 10a);而晚期节面走向明显以垂直于龙门山断裂带走向的NNW 向占优势,表明晚期以NNW 向逆冲断层活动为主(图 9b图 10a).与南段(图 8a)相比,倾角变化范围相对较小,主要分布在40°~70°之间(图 10a).滑动角分布则显示断层运动方式以逆冲为主,但相对于南段而言,其走滑分量有所增大(图 10a).

图 10 中段节面走向、滑动角、倾角(a)与P轴方位角、仰角统计(b) Fig. 10 Strikes, rakes and dips of fault planes (a) , azimuths and plunges of P axes (b) on the central

中段P轴仰角均值小于15°,较南段更接近水平.同时,与南段相比,中段早期P轴方位分布相对分散(图 10b),NE、NEE 向相对占优势,晚期则以NEE 向为主,表明中段P轴优势方向也明显有别于主震P轴方向[22]与区域应力场方向[46-51],而与主震破裂传播方向[23-25]接近.根据节面走向与P轴分布,可以推测,中段的余震一部分来源于主震后龙门山断裂带的调整性运动,另一部分则是因汶川主震引发余震区应力场的改变,并在应力场调整过程中,触发不同走向的先存或新生次级断裂、尤其是与龙门山断裂带走向近乎垂直的NNW 向断裂活动的结果.除此之外,近SN 向的岷江断裂带和虎牙断裂带在余震区中段与龙门山断裂带相交,也可能是上述余震的发震断裂之一.

(4)余震区北段

平武以东的余震区北段长约110km,共有96次地震参与统计(表 2),其中,早期余震77次,逆冲型、走滑型与正断型所占比例分别为36.4%、59.7%、3.9%,走滑型地震比例远高于逆冲型,这一特征也与已有研究结果一致[2031-32].然而,在19次晚期余震中,逆冲型地震比例达52.6%,比早期上升了16.2%;走滑型地震比例为47.4%,下降了12.3%;无正断型地震活动,由此可见,北段余震震源机制类型比例随时间变化极为显着.刁桂苓等[50]认为,震源机制类型的转换可以理解为构造应力场的改变.我们认为,在余震活动晚期所出现的南段走滑型地震比例显着上升、而北段逆冲型地震比例显着上升的现象,可能是余震区两端断层调整性运动的表现.图 11为北段早期和晚期余震机制分布,为缩小绘图范围,图 11b中未包含图 6 中远离余震密集区的宁强地震.由该图可见,与南段刚好相反,汶川主震后的较短时期内,余震表现为走滑作用为主(图 11a);在后期震后调整过程中,逆冲型地震比例大增(图 11b),由此补充和协调断裂带沿垂向方向上的应力变化.

图 11 余震区北段不同时期的余震震源机制分布 Fig. 11 Focal mechanisms of aftershocks on the northern segment of the aftershock are aduring different stages

图 1图 6可以看出,北段北端余震分布呈现为沿龙门山断裂带走向的NE 向分支与偏离龙门山断裂带走向的NNE 向分支,可能反映了汶川地震破裂带尾端破裂扩展特征[57].因此,与前3个区域相比,北段节面走向分布更为复杂,早期NE-NNE 向节面略占优势,晚期优势方向不明显(图 12a),反映了该段复杂的发震构造,其节面分布方位的离散性反映出不同走向的先存或新生构造参与了主震后的余震活动.滑动角和断层倾角显示该段早期以较大倾角的断层走滑活动为主,晚期则以中等倾角断层逆冲活动为主(图 12a).

图 12 北段节面走向、滑动角、倾角(a)与p轴方位角、仰角统计(b) Fig. 12 Strikes,rakes and dips of fault planes (a) ,azimuths and plunges of P axes(b)on the northern segment

与南段和中段一样,北段P轴仰角也呈近水平,晚期P轴仰角更小,均值约10°(图 12b).P轴方位总体分布呈现NWW 和NEE 两个优势方位,分别对应于主震前由原地应力测量获得的该段水平最大主压应力方向(NW-NWW)[58]与主震破裂传播方向[23-25],但P轴方位分布在不同时段仍然存在明显差异:余震早期P轴方位分布相对简单,集中在NWW 和NEE 两个优势方向,余震晚期P轴方位分布则相对混乱,优势方向不突出.大震后震区主压应力方向与区域应力场方向发生较大偏离、且随时间变化的特征,与李方全等[59-61]通过研究邢台、海城、龙陵、唐山等大震后震区地应力测量结果所获得的认识是一致的.这种变化特征反映了震后震源区应力调整过程.因此,北段P轴方位分布特征反映该区域余震除受区域应力场控制外,还受到汶川主震所引起的局部应力场的控制.

尤其值得关注的是,位于北段余震密集区末端的余震,早期基本为走滑型,自2008年7月24日开始,余震走滑分量明显减小而逆冲分量显着增大,且在空间位置上更偏离平武-青川断裂,北段末端余震错动方式及空间位置的变化,或许代表了汶川地震破裂和应变过程由同震及震后短期内走滑,后期调整为逆冲,这种特征符合走滑断裂位移应变向端部递减,逐渐转变为端部逆冲的规律.

2011年11 月1 日,汶川余震区北段青川与文县间再次发生5.4级强余震,CAP 波形反演结果与P波初动解均显示该地震为纯逆冲型,支持我们所获得的北段晚期余震特征.

4 龙门山断裂带中北段断面结构初步

探讨根据余震深度及震源机制解所给出的节面参数,可推测发震断层面几何形态.假定龙门山断裂带的几条主干断裂在不同深度上平行展布,选取节面走向与龙门山断裂带走向差异小于20°的地震,统计龙门山断裂带中北段不同段落倾角随深度的分布特征.

统计结果显示,沿龙门山断裂带走向,不同段落的断层面倾角范围均较大(图 13),且与深度的关系较离散.然而,倾角均值在不同段落、不同深度的变化却存在一定的规律.图 13a显示,余震区南段在10~18km 深度范围内,断层面倾角均值随深度减小而增大,由18km 时的大约40°上升至10km 时的70°;当深度小于10km 时,断层面倾角均值在60°左右,断层倾角均值的深度变化显示出深部缓倾角、浅部高倾角的铲形断面特征[101826],从而验证了Zhu等[27]的拟合结果.中段余震主要集中在7~19km深度(图 5),断面倾角均值随深度的变化特征与南段和北段存在明显差异,深度大于10km 时的倾角均值相对稳定(图 13b),在50°左右;深度小于10km 时,倾角均值存在随深度减小而增大的趋势.余震区北段不同深度上的断面倾角均值相对稳定(见图 13c),大约50°~65°,表明北段断层几何结构相对简单,断面较陡且均一,有利于走滑型破裂过程,这一结果与Shen 等[29]的研究相吻合.余震区不同段落断面几何参数的变化,揭示了龙门山断裂带中北段断面几何形态的复杂性.

图 13 余震区沿龙门山断裂带不同段落断层面倾角随深度分布 图中水平线段标示特定深度上的倾角标准差.深色影区标示倾角均值随深度的变化趋势. Fig. 13 Distribution of dip of fault plane with focal depth on different segments along Longmen Shan fault zone The horizontal line marks standard error of dips at special depth. The horizontal line marks standard error of dips atspecial depth.Grey shadow area marks the tendency of mean dip with depth.
5 结 论

本文利用CAP波形反演方法,求取了汶川地震序列中312个4级以上余震的震源机制解和震源深度,分析了余震震源深度分布与震源机制时空分布特征,探讨了余震区断层行为与应力场的时空变化,并对龙门山断裂带中北段的断层面几何形态进行了初步探讨.获得的主要结果如下:

(1)余震震源深度分布存在显着的空间差异.绵竹以西的余震区南段与平武以东的余震区北段余震深度范围大于余震区中段,但深度小于5km 的5级以上超浅源地震主要分布在明显偏离龙门山断裂带走向的余震区南段理县NW 向分支与北端NNE向分支.而余震区中段绵竹-平武段4级以上余震主要分布在7~19km 深度,鲜有深度小于7km 的4级余震活动.

(2)余震机制类型比例存在明显的时空差异.余震区南段理县NW 向分支余震以走滑型为主,中段绵竹-平武段逆冲型地震占绝对优势,时间变化不显着.沿龙门山断裂带走向的余震区南段早期逆冲型地震比例高于走滑型、晚期走滑型地震比例显着升高并超过逆冲型,而平武以东的余震区北段早期走滑型地震占绝对优势、晚期逆冲型地震比例大幅上升且超过走滑型,可能是余震区两端断层调整性运动的表现.

(3)节面走向与P轴方位优势方向存在显着的空间分段差异.这种差异反映了余震活动除了受区域应力场控制外,还受到主震引发的局部应力场的控制,而节面走向的多方位分布也反映了不同走向的构造参与了主震后的余震活动.

(4)龙门山断裂带中北段断层面几何形态复杂.余震区南段大于10km 深度时的断面倾角均值随深度减小而增大,深度小于10km 时的断面倾角均值相对稳定,显示出深部缓倾角、浅部高倾角的铲形断面特征;中段深度大于10km 时的倾角均值相对稳定,深度小于10km 时倾角均值存在随深度减小而明显增大的趋势;北段断面倾角均值在不同深度上的变化相对较小,约50°~65°,表明北段断面几何形态相对简单.上述不同段落断面几何参数的变化,揭示了龙门山断裂带中北段断面几何形态的复杂性.

余震区P轴优势方向的时空分段差异及其与区域应力场方向较大的偏离现象,以及大震破裂区两端机制类型比例的显着变化特征,可能反映了大震所引起的震源区应力场调整过程,有必要对该过程开展更进一步的深入研究,为强余震发震机理研究及强余震预测提供依据.

致谢

本文CAP 波形反演程序来源于美国圣路易斯大学的朱露培博士;图形绘制采用了GMT[62]程序;在本文研究过程中,张培震研究员就结果的合理分析进行了悉心指导,受益匪浅;审稿人给出了很好的修改建议.在此一并致以衷心感谢!

参考文献
[1] 王二七, 孟庆任, 陈智樑, 等. 龙门山断裂带印支期左旋走滑运动及其大地构造成因. 地学前缘 , 2001, 8(2): 375–384. Wang E Q, Meng Q R, Chen Z L, et al. Early Mesozoic left-lateral movement along the Longmen Shan fault belt and its tectonic implication. Earth Science Frontiers (in Chinese) , 2001, 8(2): 375-384.
[2] Densmore A L, Ellis M A, Li Y, et al. Active tectonics of the Beichuan and Pengguan faults at the eastern margin of the Tibetan Plateau. Tectonics , 2007, 26: TC4005. DOI:10.1029/2006TC001987
[3] Zhang H P, Zhang P Z, Kirby E, et al. Along-strike topographic variation of the Longmen Shan and its significance for landscape evolution along the eastern Tibetan Plateau. Journal of Asian Earth Sciences , 2011, 40(4): 855-864. DOI:10.1016/j.jseaes.2010.05.015
[4] 周荣军, 李勇, DensmoreA L, 等. 青藏高原东缘活动构造. 矿物岩石 , 2006, 26(2): 40–51. Zhou R J, Li Y, Densmore A L, et al. Active tectonics of the eastern margin of the Tibet plateau. Journal of Mineralogy and Petrology (in Chinese) , 2006, 26(2): 40-51.
[5] 邓起东, 陈社发, 赵小麟. 龙门山及其邻区的构造和地震活动及动力学. 地震地质 , 1994, 16(4): 389–403. Deng Q D, Chen S F, Zhao X L. Tectonics, scismisity and dynamics of Longmenshan Mountains and its adjacent regions. Seismology and Geology (in Chinese) , 1994, 16(4): 389-403.
[6] 许志琴, 李海兵, 吴忠良. 汶川地震和科学钻探. 地质学报 , 2008, 82(12): 1613–1622. Xu Z Q, Li H B, Wu Z L. Wenchuan earthquake and scientific drilling. Acta Geologica Sinica (in Chinese) , 2008, 82(12): 1613-1622.
[7] 李勇, 黄润秋, 周荣军, 等. 龙门山地震带的地质背景与汶川地震的地表破裂. 工程地质学报 , 2009, 17(1): 3–16. Li Y, Zhou R Q, Zhou R J, et al. Geological background of Longmenshan seismic belt and surface ruptures in Wenchuan earthquake. Journal of Engineering Geology (in Chinese) , 2009, 17(1): 3-16.
[8] 刘树根, 李智武, 曹俊兴, 等. 龙门山陆内复合造山带的四维结构构造特征. 地质科学 , 2009, 44(4): 1151–1180. Liu S G, Li Z W, Cao J X, et al. 4-D textural and structural characteristics of Longmen intracontinental composite orogenic belt, southwest China. Chinese Journal of Geology (in Chinese) , 2009, 44(4): 1151-1180.
[9] Zhou R J, Li Y, Densmore A L. Active tectonics of the Longmen Shan region on the eastern margin of the Tibetan Plateau. Acta Geologica Sinica: English Edition , 2007, 81(4): 593-604. DOI:10.1111/acgs.2007.81.issue-4
[10] 付碧宏, 时丕龙, 王萍, 等. 2008年汶川地震断层北川段的几何学与运动学特征及地震地质灾害效应. 地球物理学报 , 2009, 52(2): 485–495. Fu B H, Shi P L, Wang P, et al. Geometry and kinematics of the 2008 Wenchuan earthquake surface ruptures around the Qushan Town of Beichuan County, Sichuan: Implications for mitigation of seismic and geologic disasters. Chinese J. Geophys (in Chinese) , 2009, 52(2): 485-495.
[11] Burchfiel B C, Royden L H, van der Hilst R D, et al. A geological and geophysical context for the Wenchuan earthquake of 12 May 2008, Sichuan, People's Republic of China. GSA Today , 2008, 18(7): 4-11. DOI:10.1130/GSATG18A.1
[12] Xu X W, Wen X Z, Yu G H, et al. Coseismic reverse- and oblique-slip surface faulting generated by the 2008 Mw7|9 Wenchuan earthquake, China. Geology , 2009, 37(6): 515-518. DOI:10.1130/G25462A.1
[13] 易桂喜, 闻学泽, 辛华, 等. 2008年汶川Ms8.0地震前龙门山—岷山构造带的地震活动性参数与地震视应力分布. 地球物理学报 , 2011, 54(6): 1490–1500. Yi G X, Wen X Z, Xin H, et al. Distributions of seismicity parameters and seismic apparent stresses on the Longmenshan-Minshan tectonic zone before the 2008 Ms8.0 Wenchuan earthquake. Chinese J. Geophys (in Chinese) , 2011, 54(6): 1490-1500.
[14] 张培震, 闻学泽, 徐锡伟, 等. 2008年汶川8.0 级特大地震孕育和发生的多单元组合模式. 科学通 , 2009, 54(7): 944–953. Zhang P Z, Wen X Z, Xu X W, et al. Tectonic model of the great Wenchuan earthquake of May 12, 2008, Sichuan, China (in Chinese). Chinese Sci. Bull (in Chinese) , 2009, 54(7): 944-953.
[15] 朱介寿. 汶川地震的岩石圈深部结构与动力学背景. 成都理工大学学报 , 2008, 35(4): 348–356. Zhu J S. The Wenchuan earthquake occurrence background in deep structure and dynamics of lithosphere. Journal of Chengdu University of Technology (Science & Technology Edition) (in Chinese) , 2008, 35(4): 348-356.
[16] 滕吉文, 白登海, 杨辉, 等. 2008汶川Ms8.0地震发生的深层过程和动力学响应. 地球物理学报 , 2008, 51(5): 1385–1402. Teng J W, Bai D H, Yang H, et al. Deep processes and dynamic responses associated with the Wenchuan Ms8.0 earthquake of 2008. Chinese J. Geophys (in Chinese) , 2008, 51(5): 1385-1402.
[17] 江在森, 方颖, 武艳强, 等. 汶川8.0级地震前区域地壳运动与变形动态过程. 地球物理学报 , 2009, 52(2): 505–518. Jiang Z S, Fang Y, Wu Y Q, et al. The dynamic process of regional crustal movement and deformation before Wenchuan Ms8.0 earthquake. Chinese J. Geophys (in Chinese) , 2009, 52(2): 505-518.
[18] Zhang P Z, Wen X Z, Shen Z K, et al. Oblique, high-angle, listric-reverse faulting and associated development of strain: the Wenchuan earthquake of May 12, 2008, Sichuan, China. Annu. Rev. Earth Planet. Sci , 2010, 38: 353-382. DOI:10.1146/annurev-earth-040809-152602
[19] 闻学泽, 张培震, 杜方, 等. 2008年汶川8.0级地震发生的历史与现今地震活动背景. 地球物理学报 , 2009, 52(2): 444–454. Wen X Z, Zhang P Z, Du F, et al. The background of historical and modern seismic activities of the occurrence of the 2008 Ms8.0 Wenchuan, Sichuan Earthquake. Chinese J. Geophys (in Chinese) , 2009, 52(2): 444-454.
[20] 郑勇, 马宏生, 吕坚, 等. 汶川地震强余震(Ms≥5.6)的震源机制解及其与发震构造的关系. 中国科学D辑 , 2009, 52(6): 739–753. Zheng Y, Ma H S, Lü J, et al. Source machanism of strong aftershocks (Ms≥5. Science in China (Series D) (in Chinese) , 2009, 52(6): 739-753. DOI:10.1007/s11430-009-0074-3
[21] 罗艳, 倪四道, 曾祥方, 等. 汶川地震余震区东北端一个余震序列的地震学研究. 中国科学D辑 , 2010, 53(11): 1655–1664. Lu Y, Ni S D, Zheng X F, et al. A shallow aftershock sequence in the north-eastern end of the Wenchuan earthquake aftershock zone. Science in China (Series D) (in Chinese) , 2010, 53(11): 1655-1664. DOI:10.1007/s11430-010-4026-8
[22] 陈运泰, 许力生, 张勇等. 2008年5月12日汶川特大地震震源特性分析报告. 中国地震信息网, http://www.csi.ac.cn/sichuan/chenyuntai.pdf.2008 Chen Y T, Xu L S, Zhang Y, et al. The report of source characteristics of the May 12, 2008 Wenchuan earthquake. Earthquake Information Center, China Earthquake Administration. http://www.csi.ac.cn/sichuan/chenyuntai.pdf.2008.
[23] 王卫民, 赵连锋, 李娟, 等. 四川汶川8.0级地震震源过程. 地球物理学报 , 2008, 51(5): 1403–1410. Wang W M, Zhao L F, Li J, et al. Rupture process of the Ms 8.0 Wenchuan earthquake of Sichuan, China. Chinese J. Geophys (in Chinese) , 2008, 51(5): 1403-1410.
[24] 张勇, 许力生, 陈运泰. 2008年汶川大地震震源机制的时空变化. 地球物理学报 , 2009, 52(2): 379–389. Zhang Y, Xu L S, Chen Y T. Spatio-temporal variation of the source mechanism of the 2008 great Wenchuan earthquake. Chinese J. Geophys (in Chinese) , 2009, 52(2): 379-389. DOI:10.1002/cjg2.v52.2
[25] 赵翠萍, 陈章立, 周连庆, 等. 汶川Mw8.0级地震震源破裂过程研究: 分段特征. 科学通报 , 2009, 51(22): 3475–3482. Zhao C P, Chen Z L, Zhou L Q, et al. Rupture process of the 8. Chinese Sci. Bull (in Chinese) , 2009, 51(22): 3475-3482. DOI:10.1007/s11434-009-0425-7
[26] 张竹琪, 张培震, 王庆良. 龙门山高倾角逆断层结构与孕震机制. 地球物理学报 , 2010, 53(9): 2068–2082. Zhang Z Q, Zhang P Z, Wang Q L. The structure and seismogenic mechanism of Longmenshan high dip-angle reverse fault. Chinese J. Geophys (in Chinese) , 2010, 53(9): 2068-2082.
[27] Zhu S B, Zhang P Z. Numeric modeling of the strain accumulation and release of the 2008 Wenchuan, Sichuan, China, Earthquake. Bull. Seismol. Soc. Amer , 2010, 100(5B): 2825-2839. DOI:10.1785/0120090351
[28] Chen J, Hayes G. Preliminary result of the May 12, 2008 Mw7.9 eastern Sichuan, China Earthquake. http://earthquake.usgs.gov/earthquakes/eqinthenews/2008/us2008ryan/finite_fault.php.
[29] Shen Z K, Sun J B, Zhang P Z, et al. Slip maxima at fault junctions and rupturing of barriers during the 2008 Wenchuan earthquake. Nature Geoscience , 2009, 2(10): 718-724. DOI:10.1038/ngeo636
[30] 张国宏, 屈春燕, 宋小刚, 等. 基于InSAR同震形变场反演汶川Mw7|9地震断层滑动分布. 地球物理学报 , 2010, 53(2): 269–279. Zhang G H, Qu C Y, Song X G, et al. Slip distribution and source parameters inverted from co-seismic deformation derived by InSAR technology of Wenchuan Mw7|9 earthquake. Chinese J. Geophys (in Chinese) , 2010, 53(2): 269-279.
[31] 胡幸平, 俞春泉, 陶开, 等. 利用P波初动资料求解汶川地震及其强余震震源机制解. 地球物理学报 , 2008, 51(6): 1711–1718. Hu X P, Yu C Q, Tao K, et al. Focal mechanism solutions of Wenchuan earthquake and its strong aftershocks obtained from initial P wave polarity analysis. Chinese J. Geophys (in Chinese) , 2008, 51(6): 1711-1718.
[32] 王勤彩, 陈章立, 郑斯华. 汶川大地震余震序列震源机制的空间分段特征. 科学通报 , 2009, 54(16): 2348–2354. Wang Q C, Chen Z L, Zheng S H. Spatial segmentation characteristic of focal mechanism of aftershock sequence of Wenchuan earhtuqake. Chinese Sci.Bull (in Chinese) , 2009, 54(16): 2348-2354.
[33] Cheng W Z, Zhang Z W, Ruan X. Spatio-temporal variation and focal mechanism of the Wenchuan Ms8.0 earthquake sequence. Earthquake Science , 2009, 22(2): 109-117. DOI:10.1007/s11589-009-0109-z
[34] 郭祥云, 陈学忠, 李艳娥. 2008年5月1 2日四川汶川8.0级地震与部分余震的震源机制解. 地震 , 2010, 30(1): 50–60. Guo X Y, Chen X Z, Li Y E. Focal mechanism solutions for the Ms8.0 Wenchuan earthquake and part of its aftershocks. Earthquake (in Chinese) , 2010, 30(1): 50-60.
[35] 张致伟, 张永久, 程万正. 汶川8.0级地震序列的小震震源机制及应力场特征. 地震研究 , 2010, 33(1): 43–49. Zhang Z W, Zhang Y J, Cheng W Z, et al. Focal mechanisms and stress field of small earthquakes of the Ms8.0 Wenchuan earthquake sequence. Journal of Seismological Research (in Chinese) , 2010, 33(1): 43-49.
[36] Zhao L S, Helmberger D V. Source estimation from broadband regional seismograms. Bull. Seismol. Soc. Amer , 1994, 84(1): 91-104.
[37] Zhu L P, Helmberger D V. Advancement in source estimation techniques using broadband regional seismograms. Bull. Seismol. Soc. Amer , 1996, 86(5): 1634-1641.
[38] 吕坚, 郑勇, 倪四道, 等. 2005年11月26日九江一瑞昌Ms5.7、Ms4.地震的震源机制解与发震构造研究. 地球物理学报 , 2008, 51(1): 158–164. Lü J, Zheng Y, Ni S D, et al. Focal mechanisms and seismogenic structures of the Ms5|7 and Ms4|8 Jiujiang -Ruichang earthquakes of Nov.26, 2005. Chinese J. Geophys (in Chinese) , 2008, 51(1): 158-164.
[39] Tan Y, Zhu L P, Helmberger D V, et al. Locating and modeling regional earthquakes with two stations. J Geophys. Res , 206, 111: B01306.
[40] 龙锋, 张永久, 闻学泽, 等. 2008年8月30日攀枝花—会理6.1级地震序列ML≥4.0事件的震源机制解. 地球物理学报 , 2010, 53(12): 2852–2860. Long F, Zhang Y J, Wen X Z, et al. Focal mechanism solutions of ML≥4.0 events in the Ms6|1 Panzhihua-Huili earthquake sequence of Aug 30, 2008. Chinese J. Geophys (in Chinese) , 2010, 53(12): 2852-2860.
[41] 张瑞青, 吴庆举, 李永华, 等. 汶川中强余震震源深度的确定及其意义. 中国科学D辑 , 2008, 38(10): 1234–1241. Zhang R Q, Wu Q J, Li Y H, et al. Focal depths for moderate-sized aftershocks of the Wenchuan Ms8.0 earthquake and their implication. Science in China (Series D) (in Chinese) , 2008, 38(10): 1234-1241.
[42] 吕坚, 苏金蓉, 靳玉科, 等. 汶川8.0级地震序列重新定位及其发震构造初探. 地震地质 , 2008, 30(4): 917–925. Lü J, Su J R, Jin Y K, et al. Discussion on relocation and seismo-tectonics of the Ms8.0 Wenchuan earthquake sequences. Seismology and Geology (in Chinese) , 2008, 30(4): 917-925.
[43] 陈九辉, 刘启元, 李顺成, 等. 汶川Ms8.0地震余震序列重新定位及其地震构造研究. 地球物理学报 , 2009, 52(2): 390–397. Chen J H, Liu Q Y, Li S C, et al. Seismotectonic study by relocation of the Wenchuan Ms8.0 earthquake sequence. Chinese J. Geophys (in Chinese) , 2009, 52(2): 390-397. DOI:10.1002/cjg2.v52.2
[44] 黄媛, 吴建平, 张天中, 等. 汶川8.0级大地震及其余震序列重定位研究. 中国科学D辑 , 2008, 51(12): 1703–1712. Huang Y, Wu J P, Zhang T Z, et al. Relocation of the Ms8.0 Wenchuan earthquake and its aftershock sequence. Science in China (Series D) (in Chinese) , 2008, 51(12): 1703-1712. DOI:10.1007/s11430-008-0135-z
[45] 朱艾斓, 徐锡伟, 刁桂苓, 等. 汶川Ms 8.0地震部分余震重新定位及地震构造初步分析. 地震地质 , 2008, 30(3): 759–767. Zhu A L, Xu X W, Diao G L, et al. Relocation of the Ms8.0 Wenchuan earthquake sequence in part: preliminary seismotectonic analysis. Seismology and Geology (in Chinese) , 2008, 30(3): 759-767.
[46] 成尔林. 四川及其邻区现代构造应力场和现代构造运动特征. 地震学报 , 1981, 3(3): 231–241. Cheng E L. Recent tectonic stress field and tectonic movement of the Sichuan Province and its vicinity. Acta Seismological Sinica (in Chinese) , 1981, 3(3): 231-241.
[47] 许忠淮, 汪素云, 黄雨蕊, 等. 由多个小震推断的青甘和川滇地区地壳应力场的方向特征. 地球物理学报 , 1987, 30(5): 476–486. Xu Z H, Wang S Y, Huang Y R, et al. Directions of mean stress axes in southwestern China deduced from microearthquake data. Chinese J. Geophys (in Chinese) , 1987, 30(5): 476-486.
[48] 陈天长, 崛内茂木, 郑斯华. 利用P波初动和短周期P、S波振幅测定川滇地区地震震源机制解和应力场. 地震学报 , 2001, 23(4): 436–440. Chen T C, Horiuchi S, Zheng S H. Earthquake focal mechanisms and stress field in Sichuan-Yunnan area determined using P wave polarity and short period P and S waveform data. Acta Seismological Sinica (in Chinese) , 2001, 23(4): 436-440.
[49] 安其美, 丁立丰, 王海忠, 等. 龙门山断裂带的性质与活动性研究. 大地测量与地球动力学 , 2004, 24(2): 115–119. An Q M, Ding L F, Wang H Z, et al. Research of property and activity of Longmen Mountain fault zone. Journal of Geodesy and Geodynamics (in Chinese) , 2004, 24(2): 115-119.
[50] 刁桂苓, 徐锡伟, 陈于高, 等. 汶川Mw7.9和集集Mw7.6地震前应力场转换现象及其可能的前兆意义. 地球物理学报 , 2011, 54(1): 128–138. Diao G L, Xu X W, Chen Y G, et al. The precursory significance of tectonic stress field transformation before the Wenchuan Mw7.9 earthquake and the Chi-Chi Mw7.6 earthquake. Chinese J. Geophys (in Chinese) , 2011, 54(1): 128-138.
[51] 李延兴, 张静华, 周伟, 等. 汶川Ms8.0地震孕育发生的机制与动力学问题. 地球物理学报 , 2009, 52(2): 519–530. Li Y X, Zhang J H, Zhou W, et al. The mechanism and dynamics of the generation and occurrence for Wenchuan Ms8.0 earthquake. Chinese J. Geophys (in Chinese), (in Chinese) , 2009, 52(2): 519-530.
[52] 许忠淮, 汪素云, 黄雨蕊, 等. 由大量的地震资料推断的我国大陆构造应力场. 地球物理学报 , 1989, 32(6): 636–647. Xu Z H, Wang S Y, Huang Y R, et al. The tectonic stress field of Chinese Continent deduced from a great number of earthquakes. Chinese J. Geophys (in Chinese) , 1989, 32(6): 636-647.
[53] Yamakawa N. Stress fields in focal regions. J. Phys Eart , 1971, 19(4): 347-353. DOI:10.4294/jpe1952.19.347
[54] Xu Z H, Wang S Y. A possible change in stress field orientation due to the 1976 Tangshan earthquake. Pure and Applied Geophysics , 1986, 124(4-5): 941-955. DOI:10.1007/BF00879620
[55] Hao K X, Si H J, Fujiwara H, et al. Coseismic surface-rupture and crustal deformations of the 2008 Wenchuan earthquake Mw7.9, China. Geophys. Res. Lett , 2009, 36: L11303. DOI:10.1029/2009GL037971
[56] 赵珠, 汪碧澜, 龙思胜. 1999年四川绵竹清平5.0地震序列构造活动特征. 地震地磁观测与研究 , 2001, 22(6): 30–37. Zhao Z, Wang B L, Long S S, et al. The tectonic active characteristic of the earthquake sequence with magnitude 5.0 in Mianzhu Qingping, Sichuan in 1999. Seismological and Geomagnetic Observation and Research (in Chinese) , 2001, 22(6): 30-37.
[57] 邓起东, 陈桂华, 朱艾斓. 关于2008年汶川Ms8.0地震震源断裂破裂机制几个问题的讨论. 中国科学(D辑) , 2011, 54(9): 1360–1377. Deng Q D, Chen G H, Zhu A L. Discussion of rupture mechanisms on the seismogenic fault of the 2008 Ms8.0 Wenchuan earthquake. Science in China (Series D) (in Chinese) , 2011, 54(9): 1360-1377. DOI:10.1007/s11430-011-4230-1
[58] 郭啟良, 王成虎, 马洪生, 等. 汶川Ms8.0级大震前后的水压致裂原地应力测量. 地球物理学报 , 2009, 52(5): 1395–1401. Guo Q L, Wang C H, Ma H S, et al. In-situ hydro-fracture stress measurement before and after the Wenchuan Ms8.0 earthquake of China. Chinese J. Geophys (in Chinese) , 2009, 52(5): 1395-1401.
[59] 李方全, 王连捷. 华北地区地应力测量. 地球物理学报 , 1979, 22(1): 1–8. Li F Q, Wang L J. Stress measurements in north China. Chinese J. Geophys (in Chinese) , 1979, 22(1): 1-8.
[60] 李方全, 孙世宗, 李立球. 华北及郯庐断裂带地应力测量. 岩石力学与工程学报 , 1982, 1(1): 73–86. Li F Q, Sun S Z, Li L Q. In-situ stress measurements in north China and Tancheng-Lujiang faulz zone. Chinese Journal of Rock Mechanics and Engineering (in Chineae) (in Chinese) , 1982, 1(1): 73-86.
[61] 李方全, 刘光勋. 地应力测量、地壳上部应力状态与地震. 中国地震 , 1986, 2(1): 50–55. Li F Q, Liu G X. Stress measurement, stress state of upper crust and earthquake research. Earthquake Research in China (in Chinese) , 1986, 2(1): 50-55.
[62] Wessel P, Smith W H F. Free software helps map and display data. Eos Trans| AGU , 1991, 72(41): 441. DOI:10.1029/90EO00319