0 引言

龙门山断裂带地处南北地震带中段，介于青藏高原东缘巴颜喀拉与华南地块、四川盆地之间，是1条NE向断裂带。以映秀和北川附近为界，可将龙门山断裂带分为南、中、北3段(邓起东等，1994；陈国光等，2007)。2008年5月12日在龙门山断裂带中段发生了MW7.9汶川大地震，主震震源破裂过程研究(王卫民等，2008；张勇等，2008)表明汶川地震是1次以逆冲为主兼具少量右旋走滑的事件。主要的动力来自巴颜喀拉块体朝SEE方向的运动受到华南地块阻挡而在龙门山断裂带的长期应力积累。龙门山断裂带能孕育出如此大的地震可归结于3大特征(Burchfiel et al.，2008；张培震等，2008；李延兴等，2009；Zhang et al.，2010)，组成龙门山断裂带的是强度大的杂岩体，断裂带剖面为铲型的构造结构，处于低活动速率的闭锁状态下。

汶川地震发生后，在龙门山断裂带上产生了大量的余震。丰富的余震资料以及大尺度破裂带的存在为研究龙门山断裂带动力学提供了重要的数据基础和模型参数。从震源机制解推断应力场是现今研究应力场主要思路之一。许多学者及机构对汶川中强余震进行了震源机制求解，获得了其时空变化特征及蕴含的构造意义(胡幸平等，2008；张勇等，2009；郑勇等，2009；王勤彩等，2009；Chen et al.，2011；崔效锋等，2011；易桂喜等，2012；Yuan et al.，2013)。通过震源机制反演应力张量，王勤彩等(2009)发现最大主应力方向沿主破裂方向发生改变，总体呈现由EW到NW-SE再到EW的特征；盛书中等(2012)从应力场角度对汶川地震断层进行南、北段划分，发现在北川附近是构造应力场及发震断层的主要分界点；Long等(2012)获得了汶川余震带内各分段应力场时空特征；Zhang等(2015)分析了汶川地震前后应力场特征及主震对周边断层的影响。这些结果提供了对龙门山断裂带构造应力场主要特征的认识，但由于现有研究主要针对4级以上的强余震，对局部应力场研究则稍显不足。

中强地震所需的驱动力数量级远大于小震，这决定了其震源机制主要受区域构造应力的控制，对局部应力场变化不敏感。相对而言，小震不仅数量上占优，而且除受构造应力场影响外，还对大地震触发的库仑应力变化以及周边断层的影响敏感。小震的缺点在于信噪比低导致的有效观测数据少，随着近年来密集台网的架设及区域精细速度结构的研究，小震定位结果及震源参数求解的质量有了显著提高。为了获得龙门山断裂带中段应力场更精细的时空变化，需要采用小震震源机制解资料。冯梅等(2013)通过对汶川地震后震源区部分小震的分析，发现震后水平最大主应力方向转变为SN向。但他们所用的台站和地震数量较少，要对汶川地震后区域应力场变化进行细致研究，仍然需要开展进一步的工作。

本文利用川西密集流动地震台阵的布局优势，从小震(震级范围为ML2.5～3.5)震源机制解入手，对龙门山断裂带中段局部应力场的时空变化特征进行研究，以便更好地了解龙门山断裂带中段应力场的震后调整变化过程，希望从中获得新的认识。

1 数据与方法

1.1 震源机制求解

对于中强地震通常采用宽频带波形数据进行矩张量反演；对微小震而言，由于整体信噪比过小，利用完整的某段波形资料进行反演就不那么适合，只能从中提取出信噪比相对高的离散数据，如P波初动极性、以及P和S振幅等。对于震源机制解近似为双力偶模型的情形下，P波初动极性物理性质清晰，求解快捷方法成熟。在观测数据较少时，如果能补充些信噪比高的振幅数据可增强对节面的约束。本文拟通过P波初动极性结合SH/P振幅比反演计算震源机制解(Hardebeck et al.，2002，2003)。

震源机制解的结果除了受所选用方法以及数据影响外，还与地震定位结果和区域速度结构相关，影响体现在回溯到震源球上的方位角和离源角。为了提高小震震源机制解的质量，我们采用邓文泽等(2014)用双差层析成像联合反演的精定位震源位置和区域三维速度模型结果计算方位角和离源角，P和SH理论振幅是根据双力偶位错点源在半无限空间远场的位移响应计算得来的。

互联网时代，机遇与危险并存。其为管理会计的发展提供了很多有利的条件，但同时也面临很多问题。诸如：对现有的认知存在误区、技术水平不达标、缺乏系统化的理论以及缺乏专业型人才等等。

1.2 应力场反演

从震源机制解反演应力场，虽不能测定地震孕育处完整的应力张量，但可以得到应力主轴的方向以及它们之间的相对大小(Gephart et al.，1984; Michael，1987)。

在应力场均匀的前提下，假设断层面上的剪应力方向与滑动方向一致

临床病理特征 两组患者性别差异无统计学意义，肿瘤疾病相关资料有可比性，但老年组的合并症明显高于年轻组，包括糖尿病(P=0.045)、心脏病(P=0.039)、呼吸系统疾病(P=0.041)、高血压病(P=0.035)等(表1)。

(1)

式(1)中， 是单位滑动矢量， 是断层面上的剪应力矢量，σ是应力张量，是断层面单位法向矢量。

式(2)中，S1、S2和S3分别是最大、中间和最小主压应力。

R=(S2-S3)/(S1-S3)

(2)

从式(1)可以确定3个应力主轴的方向和1个描述不同主应力相对大小的应力比R，一般定义其为

我抢在别人之前拾起传票，心里明白是怎么回事了。我得镇住混乱的现场。我说，李书记闹肚子，这两天一直带病坚守岗位呢。接着，我即兴发挥把欢迎词致了，又请嘉宾讲话，我应变能力还强吧，临时决定将行长一军，请他作重要指示。行长满脸尴尬，对着麦说，谢谢同学们，我没有重要指示，我只有美好的祝愿，祝愿同学们好好学习天天向上！我用喊声压住掌声说，这个指示还不重要啊！这曾经是毛主席的伟大指示啊！

为了评估本文小震震源机制解方法的稳定性，本文随机选取崔效锋等(2011)文章中提及的在本文研究区域内的4级以上地震20个，采取与计算小震相同的步骤方法反演计算这些震源机制解。为了使强余震与小余震的反演结果具有可比性，只选取了强余震周边10—13个台站的数据进行震源机制反演。图2 给出了两者震源机制解对比图。为了量化2个震源机制解的相似程度，一般用以P、T、B为轴的坐标轴之间的最小旋转角来描述。最小旋转角在30°以内，两者的震源机制解相似性高。最小旋转角对比统计结果显示，55%的测量结果最小旋转角都在30°以内，95%都在50°以内，这说明两者给出的震源机制解总体上具有很好的一致性，表明了该算法的可靠性。

(3)

使应力场反演变得线性化。并通过bootstrap技术消除从震源机制2参考节面挑选正确断层面的技术难题。

图 1 川西台站分布以及地震余震分布 Fig. 1 Map showing seismic stations of western Sichuan and aftershocks of the 2008 Wenchuan earthquake. 红色三角代表地震台站；黄色的五星代表汶川主震；绿色和黑色小点代表重定位后的余震分布(邓文泽等，2014)，其中绿色小点为本文所选用的数据；白色矩形框为本文研究区域；黑色实线代表断裂，主要有： ①汶川-茂县断裂， ②映秀-北川断裂，③灌县-江油断裂

历朝统治者为了自身及国家安全，往往对海路采取过于谨慎的态度，大规模的商贸大多难以正常进行。但即使是在这种较为封闭的环境下，海上贸易因其便捷性仍然缓慢发展着。例如古代胶东半岛，作为南方沿海与京城联系的海上枢纽，长久以来广东、福建等地商人乘船在此贩运，商贾云集，之后更成为海上丝绸之路的起点。胶东传统民居具有浓厚的南北交融特色，山墙与潮汕地区类似，墙面抹灰与闽粤民居雷同，而硬山顶屋面屋脊处理则又具有北方民居风格。明清以来，随着沿海口岸的开放，不止闽粤沿海，通商口岸从宁波开始不断扩散，海外贸易逐渐发展。沿海聚落中最为有统一特色的就是妈祖文化，通过每村设置的公共建筑进行祭祀可以看出其同源性。

1.3 研究区域和数据

本文主要研究区域为龙门山中段映秀至北川这一片余震带(图1)。选取的数据为汶川地震后2008年5—10月发生在龙门山中段的余震事件，经邓文泽等(2014)重定位后震级在2.5～3.5之间的小震事件675个，观测台站为中国地震局在川西地区架设的流动台阵。为了提高震源机制解质量，我们限制震中距为120km，最少有效记录台站数为8个，最后符合要求的余震事件共有638个。

1.4 震源机制解质量评估

此外，由式(1)进行的反演是非线性的。Michael(1984)在上述框架内，额外增添了所有断层面上剪应力大小一致的假设

本文采用Hardebeck等(2006)的阻尼线性反演方法分析龙门山断裂带中段应力场的时空特征。在反演过程中，需挑选1个节面作为实际断层面；小震的特性决定了其震源机制的多样性，为了不人为干预反演结果，将每个节面被选到的概率都设为0.5。将阻尼因子设为1，并进行2i000次bootstrap重采样，以获得相对稳定的95%置信区间。

图 2 不同方法获得的震源机制解 Fig. 2 Focal mechanism solutions from different methods. a 相邻的2个震源机制右边为崔效锋等(2011)的结果，左边为本文结果，震源机制序号与崔效锋等(2011)文章中的地震序列序号相同；b 2个震源机制解最小旋转角的直方统计图

2 结果与讨论

2.1 汶川小余震的震源机制解

步骤1:初始化,t=0,t为当前迭代次数,N为种群大小。若k≥N,则根据适应度值从小到大选取前N个作为初始值将(xP1,yP1),(xP2,yP2),…,(xPi,yPi),…,(xPN,yPN)赋值给(xW1(t),yW1(t)),…,(xWi(t),yWi(t)),…,(xWN(t),yWN(t)),i的取值为1～N;若k

图 3 不同类型震源机制解分布 Fig. 3 Spatial distribution of different kinds of focal mechanism solutions. a 逆冲、走滑和正断分别用红色、蓝色和绿色表示；b 震后早、晚期各类型震源机制数目所占比例直方图； c 根据Frohlich(1992)画出的颜色比例三角图

图 4 震源机制解沿垂直龙门山断裂带方向的投影 Fig. 4 Projections of focal mechanism solutions normal to the Longmenshan fault zone. a 震后早期，b 震后晚期；AA′、BB′、CC′、DD′分别对应图1 中的子区域；WMF、YBF、GJF分别代表汶川-茂汶断裂、映秀-北川断裂和灌县-江油断裂

参考陈九辉等(2009)给出的余震空间形态特征将龙门山中段分为4个区域(图1)。图4 给出震源机制在各区域垂直龙门山断裂带剖面上的投影。

对比震后不同时期震源机制空间形态的分布特征发现，与主震性质差异较大的正断型事件主要集中在汶川地震破裂带南端的2个区域，一是主震起始破裂区西北侧，二是小鱼洞-理县次级破裂西北侧。震源区主要受局部地质构造影响。小鱼洞-理县断裂左旋走滑，破裂东北侧存在明显的NW向运动，受中央断裂切割出现解耦，并在局部表现出较强的拉张性质。强余震震源机制研究(郑勇等，2009)表明震后早期在小鱼洞-理县NW-SE向左旋走滑断裂区域含有一定的拉张运动。相对而言，破裂西南侧沿SE运动则会在断层交界处积累一定的应变能，震后晚期正断型聚集延伸下方出现逆冲型集中可能就是应力释放的表现。4个分区上、下盘震源机制类型随时间分布变化明显，表明震后应力场调整、变化显著，与时间密切相关。

首先对汶川地震的余震进行震源机制解测量。在638个参与反演计算的余震事件中，最终求得质量相对较好的震源机制解488个。图3 给出了这些震源机制解的空间分布图。图3 显示，汶川地震的余震震源机制以逆冲型和走滑型为主，但相比前人基于强余震的同类研究，我们获得的小震震源机制解结果中同时具有相当比例的正断型地震。强余震具有较大的破裂长度，受较大范围应力场的控制，而小震则主要受局部应力场控制，这可能是我们获得的结果与强余震结果(胡幸平等，2008；王勤彩等，2009；易桂喜等，2012)存在差异的主要原因。为了进一步了解震后小震震源机制分布特征，对小震震源机制解随时间的变化趋势进行了研究。以2008年6月12日为界，将震后时间分成早期和晚期2个时段。时间界限的划分以6月9日和6月11日在龙门山中段发生的2个强余震为主要依据；这2个强余震发生之后，汶川地震余震数量开始逐渐减少，认为震后应力调整进入了新的阶段。统计结果显示，早期余震中逆冲、走滑、正断与其他型所占比例分别为37.6%、19.3%、18.7%和24.4%，晚期则分别为47.7%、22.8%、12.9%和16.6%。随时间最明显的趋势是逆冲型比例大幅度上升，正断及其他型比重下降。

到了震后晚期(图4b)，经过一段时间的应力调整和释放过程之后，AA′区域中央断裂和灌县-江油断裂上仍然显示清晰的逆冲分量，而在中央断裂NW侧余震数目减少，原先的逆冲型震源机制多被走滑和正断所代替，后山断裂上盘浅部逆冲活动明显。BB′区域断层分界线变得模糊，中央断裂下盘显著的逆冲性质消失，而在原先正断型聚集延伸下方出现逆冲型集中，表明局部有较强的逆冲活动。CC′区域最明显的变化是整体上逆冲型所占比重大幅度上升，其他性质小震比重降低。DD′剖面断裂下盘及断层附近逆冲型数目增多，上盘NW侧震源机制变化明显，由早先的逆冲转变为走滑。震后晚期CC′和DD′剖面上均出现1条主要为逆冲性质的向四川盆地倾斜的余震分布带，这可能是汶川地震后应力调整造成的1个次级逆冲破裂，破裂进入四川盆地长达20km。

图4a给出了早期余震震源机制分布的剖面图。AA′剖面穿过主震起始破裂区，震后初期AA′剖面小震主要以逆冲为主。但在震中NW侧10km处震源机制含有较多的拉张分量，推测是由于主震的发生造成的应力释放导致的结果。BB′剖面所属区域是小鱼洞—理县NW向余震分支起始段，震源机制类型显示其沿NW变化明显，中央断裂与前山断裂之间小震活动以逆冲为主，中央断裂上盘靠近断层处聚集大量正断型事件，再往NW则走滑型数目开始增多。CC′区域是主震由纯逆冲型破裂向逆冲加走滑破裂转换的起始区域，余震震源机制各类型混杂，几乎看不出分布规律。DD′区域震源机制显示断裂上盘靠近断层处震源性质多为走滑和正断，沿NW方向移动逆冲型比例大幅度增加，显示上盘仍然具有较强的残存应力场。震后初期地表垂直位错和断层面上破裂滑动量较大(徐锡伟等，2008；Shen et al.，2009)的AA′和DD′地区，上盘逆冲活动明显。

缺钙发生原因：通常在土壤酸度较大（pH较低）的情况下易导致缺钙。柑橘对钙的需求量较大，但由于在柑橘生产中常施石灰、过磷酸钙、钙镁磷肥及喷含钙的肥料、农药等，使得柑橘园缺钙的现象并不多见。仅山坡地，或土质差、土壤有机质含量低时，钙素易流失，会导致土壤缺钙。其次，大量施用生理酸性的化肥，易使土壤酸化，并加速钙的流失。另外，在干旱年份土壤水分不足时，氮和其他盐离子（如钾）浓度相应提高，影响根系对钙的吸收，也会发生暂时缺钙。近年发现柑橘裂果趋势越来越严重，与土壤营养失衡密切相关，土壤过量的磷和钾是造成柑橘缺钙裂果的主要原因。

2.2 震后应力场

图5 给出了震后早、晚期各区域应力场的反演结果，并根据Lund等(2007)的方法计算了水平最大主压应力SH方位。为了解不同深度间的应力场差异，在空间分段的基础上，以15km为界，按深度划分为10～15km和15～20km 2个区间。

图 5 早、晚期各分区应力张量结果 Fig. 5 Stress tensors from inversion results in each area for two different periods. a 早期10～15km深度；b 早期15～20km深度；c 晚期10～15km深度；d 晚期15～20km深度。左边给出各个主压应力轴在下半球等面积投影的95%置信区间分布，最佳拟合主压应力轴用空心符号表示，A、B、C、D分别对应图1 中的不同区域；中间给出基于95%置信区间应力比的频数直方图；右边给出最大水平主应力SH的方位分布，应力结构根据最佳拟合应力张量结果

从图 5a可以看出，在10～15km深度范围内，AA′、CC′和DD′区域S1近水平，方位十分相似，都是NWW-SEE，计算的SH方位结果类似。BB′区域SH方位与其他地区差异较大，应该是受小鱼洞-理县走滑破裂影响，SH方位转变为NE-SW，几乎与断层走向一致，与相邻区域的应力场方向有近60°的差异。AA′与CC′区域S2和S3各有一定的倾角，说明应力结构是走滑与逆冲的混合体，BB′区S2较S3更接近于垂直，表明BB′区域为走滑型应力结构，与NW向走滑断裂相吻合。DD′区域S2和S3相互覆盖，最佳拟合应力张量更支持逆冲型断层。

在15～20km深度(图5b)，AA′和BB′区域SH方位相似近EW，CC′和DD′则是NWW-SEE。沿破裂面有呈现出由EW顺时针旋转至NWW-SEE的微弱变化趋势。AA′和DD′区域S3虽范围广但接近于垂直，反映出逆冲型应力结构。BB′区域S2和S3应力结构则是走滑与逆冲的混合体。CC′区域震源机制类型较零散(图5a)，但应力场反演结果显示走滑特征十分显著。

震后初期，不同深度间应力场差别最大的就是小鱼洞—理县地区(图5 BB′剖面)。速度结构研究显示(邓文泽等，2014)小鱼洞地区在15km以上两侧存在速度差异，与本文SH方位陡变结果类似，结合震源机制投影图(图5a)可推测小鱼洞-理县NW向走滑型地震的震源深度在15km以上。15km以下仍受主震NE向破裂传播特征的控制。从SH方位差异可看出，与应力结构相比，SH方位对破裂面特征更敏感。

从图5c可见，震后晚期10～15km深度间SH方位基本在EW和NWW-SEE之间。分布规律与图5b类似，猜测浅部SH方位的变化可能追随深部的变化，并有一定的延迟。小鱼洞区SH方位不再沿断层走向，而是逐渐与其他区域趋向一致，表明经过应力调整过程，NW向走滑断裂引起的动态应力对SH方位的影响已经弱化，不再占据主导地位。从应力张量看，这种变化表现在各应力主轴置信区间均很大，但从其格局分布上可看出该区断层形态仍以走滑为主。

随着国家文化事业的改革和发展，高校图书馆的服务功能正在发生深刻变化，如今的高校图书馆已经突破了传统的以图书借阅为主的单一服务功能，正在向多元化服务的方向发展，除了简单的校内文献信息资源服务外，还具有社会教育、社会信息服务等职能.与其它类型的图书馆相比，高校图书馆在硬件设施和软件条件方面都具有明显的优势.

在15～20km深度(图5d)，SH方位变化趋势不再一致，AA′区十分离散，优势方位不明显，BB′顺时针旋转至NW-SE，CC′仍保持NWW-SEE不变，DD′方位离散度增加并有逆时针旋转趋势。整个区域S3近垂直，说明应力结构为逆冲。BB′区SH方位比早期顺时针旋转明显，显著的变化可能归结于后期的应力释放。Long等(2013)发现北川附近S1方位为NEE-SWW并且沿NW向方位有逆时针旋转变化。本文的结果与其一致，沿NW向S1方位变化原因体现在震源机制由逆冲向走滑转变(图4b DD′剖面)。从SH方位及应力结构变化可以看出深部应力有了一定的调整和恢复，可能向着区域背景应力场转变。

王凯英等(2012)得到的汶川地震前应力场结果提出汶川地震震中区是1个应力梯度带。从震后应力场研究中发现这种应力突变区域主要在15km深度以下。应力的陡变体现在SH方位的模糊及应力椭球的形态上，应力比R数值较大，在0.7左右，说明S1和S2十分接近。在此情况下，应力结构仍为逆冲型，推测震源区深部彭灌杂岩存在高强度的局部隆升凸起，导致彭灌杂岩NW、SE和SW边界侧SH方位指向中心，造成SH方位的模糊。震后晚期应力场特征也观察到了这一点。速度结构揭示的震源区下方中央断裂和后山断裂之间包夹着显著的高速体(邓文泽等，2014)可能就是局部隆升的产物。由于彭灌杂岩强度大，在其NE侧会出现由局部隆升导致的物质横向变化，并在边界积累巨大的剪应变，小鱼洞-理县NW向走滑断裂成因可能与之相关。

3 结论

本文利用川西流动地震台阵记录的汶川地震余震波形数据，求解了小震震源机制解，并通过震源机制解结果反演了龙门山中段应力场分布，展示了其时空变化特征。获得的主要认识有：

“雪龙”号11月2日从上海出发，计划2019年4月12日返回上海，将进行为期162天、总航程约3.7万海里（约合6.85万公里）的南极综合科学考察。本次航行，“雪龙”号将6次穿越以狂风巨浪著称的“魔鬼西风带”，科考队员需同步建设考察站并完成科学研究计划，困难与挑战空前。

(1)对汶川小余震进行震源机制求解，发现震源机制以逆冲型为主，同时具有相当比例的正断型地震。各类型比重随时间变化明显，逆冲型随时间比例大幅度上升，正断及其他型比例下降。震源机制投影发现余震性质的空间分布与断层深部展布及区域速度结构具有相关性，时空形态揭示了局部应力场的复杂性。

(2)应力张量的反演结果显示，震后初期小鱼洞浅部地区受左旋走滑破裂影响，SH方位为NE-SW向，其他区域SH方位基本在EW和NWW-SEE范围内，表明震后初期应力场主要受主震破裂传播引起的动态应力影响。破裂传播引起的动态应力使小震破裂面走向与传播方向趋于一致，表现在SH方位趋向一致。随着时间的推移，应力场趋于平衡的过程中，深部应力结构由逆冲、走滑混合型向纯逆冲型转变，表明深部区域通过走滑来平衡破裂传播造成的应力失稳，并向背景应力场转变。浅部区域的应力调整追随深部的变化特征，并有一定的滞后。

(3)应力场反演揭示主震震中区深部SH方位模糊且为逆冲型应力结构，表明局部存在高强度的隆升凸起。先前发现的应力异常可能来自于20km深度左右彭灌杂岩的隆升，杂岩体内NE侧会出现由局部隆升导致的物质横向变化，小鱼洞-理县走滑断裂的成因可能与之有关。

分别准确称取12.5 mg苦杏仁苷与5.0 mg野黑樱苷标准品，甲醇溶解，分别定溶于25 mL、10 mL容量瓶中，即分别得到质量浓度均为500 μg/mL的苦杏仁苷与野黑樱苷甲醇标准储备液，4 ℃储存备用。分别吸取2 mL各标准储备液至不同10 mL容量瓶中，甲醇稀释至刻度，即得100 μg/mL苦杏仁苷标准溶液和100 μg/mL野黑樱苷标准溶液，4 ℃储存备用。

参考文献

陈国光，计凤桔，周荣军，等. 2007. 龙门山断裂带晚第四纪活动性分段的初步研究 [J]. 地震地质，29(3): 657— 673. doi: 10.3969/j.issn.0253-4967.2007.03.019.

CHEN Guo-guang，JI Feng-ju，ZHOU Rong-jun，et al. 2007. Primary research of activity segmentation of Longmenshan fault zone since late-quaternary [J]. Seismology and Geology，29(3): 657— 673(in Chinese).

陈九辉，刘启元，李顺成，等. 2009. 汶川MS8.0地震余震序列重新定位及其地震构造研究 [J]. 地球物理学报，52(2): 390—397.

CHEN Jiu-hui，LIU Qi-yuan，LI Shun-cheng，et al. 2009. Seismotectonic study by relocation of the Wenchuan MS8.0 earthquake sequence [J]. Chinese Journal of Geophysics，52(2): 390—397(in Chinese).

崔效锋，胡幸平，俞春泉，等. 2011. 汶川地震序列震源机制解研究 [J]. 北京大学学报(自然科学版)，47(6): 1063—1072.

CUI Xiao-feng，HU Xing-ping，YU Chun-quan，et al. 2011. Research on focal mechanism solutions of Wenchuan earthquake sequence [J]. Acta Scientiarum Naturalium Universitatis Pekinensis，47(6): 1063—1072(in Chinese).

邓起东，陈社发，赵小麟. 1994. 龙门山及其邻区的构造和地震活动及动力学 [J]. 地震地质，16(4): 389— 403.

DENG Qi-dong，CHEN She-fa，ZHAO Xiao-lin. 1994. Tectonics，seismicity and dynamics of Longmenshan mountains and its adjacent regions [J]. Seismology and Geology，16(4): 389— 403(in Chinese).

邓文泽，陈九辉，郭飚，等. 2014. 龙门山断裂带精细速度结构的双差层析成像研究 [J]. 地球物理学报，57(4): 1101—1110.

DENG Wen-ze，CHEN Jiu-hui，GUO Biao，et al. 2014. Fine velocity structure of the Longmenshan fault zone by double-difference tomography [J]. Chinese Journal of Geophysics，57(4): 1101—1110(in Chinese).

冯梅，张纪中，安美建. 2013. 2008年汶川地震前后震源区构造应力场变化 [J]. 地震地质，35(4): 701—720. doi: 10.3969/j.issn.0253-4967.2013.04.002.

FENG Mei，ZHANG Ji-zhong，AN Mei-jian. 2013. Tectonic stress changes around the epicenter before and after the 2008 Wenchuan earthquake [J]. Seismology and Geology，35(4): 701—720(in Chinese).

胡幸平，俞春泉，陶开，等. 2008. 利用P波初动资料求解汶川地震及其强余震震源机制解 [J]. 地球物理学报，51(6): 1711—1718.

HU Xing-ping，YU Chun-quan，TAO Kai，et al. 2008. Focal mechanism solutions of Wenchuan earthquake and its strong aftershocks obtained from initial P wave polarity analysis [J]. Chinese Journal of Geophysics，51(6): 1711—1718(in Chinese).

李延兴，张静华，周伟，等. 2009. 汶川MS8.0地震孕育发生的机制与动力学问题 [J]. 地球物理学报，52(2): 519—530.

LI Yan-xing，ZHANG Jing-hua，ZHOU Wei，et al. 2009. The mechanism and dynamics of the generation and occurrence for Wenchuan MS8.0 earthquake [J]. Chinese Journal of Geophysics，52(2): 519—530(in Chinese).

盛书中，万永革. 2012. 由构造应力场研究汶川地震断层的分段性 [J]. 地震学报，34(6): 741—753.

SHENG Shu-zhong，WAN Yong-ge. 2012. Segmentation of the Wenchuan earthquake fault derived from tectonic stress analysis [J]. Acta Seismologica Sinica，34(6): 741—753(in Chinese).

王凯英，Rebetsky Y，刁桂苓，等. 2012. 2008年汶川地震前震源区的异常应力状态 [J]. 地震地质，34(1): 11—16. doi: 10.3969/j.issn.0253-4967.2012.01.002.

WANG Kai-ying，Rebetsky Y，DIAO Gui-ling，et al. 2012. Stress state of Wenchuan area before the 2008 strong earthquake [J]. Seismology and Geology，34(1): 11—16(in Chinese).

王勤彩，陈章立，郑斯华. 2009. 汶川大地震余震序列震源机制的空间分段特征 [J]. 科学通报，54(16): 2348—2354.

WANG Qin-cai，CHEN Zhang-li，ZHENG Si-hua. 2009. Spatial segmentation characteristic of focal mechanism of aftershock sequence of Wenchuan earthquake [J]. Chinese Science Bulletin，54(13): 2263—2270.

王卫民，赵连锋，李娟，等. 2008. 四川汶川8.0级地震震源过程 [J]. 地球物理学报，51(5): 1403—1410.

WANG Wei-min，ZHAO Lian-feng，LI Juan，et al. 2008. Rupture process of the MS8.0 Wenchuan earthquake of Sichuan，China [J]. Chinese Journal of Geophysics，51(5): 1403—1410(in Chinese).

徐锡伟，闻学泽，叶建青，等. 2008. 汶川MS8.0地震地表破裂带及其发震构造 [J]. 地震地质，30(3): 597— 629. doi: 10.3969/j.issn.0253-4967.2008.03.003.

XU Xi-wei，WEN Xue-ze，YE Jian-qing，et al. 2008. The MS8.0 Wenchuan earthquake surface ruptures and its seismogenic structure [J]. Seismology and Geology，30(3): 597— 629(in Chinese).

易桂喜，龙锋，张致伟. 2012. 汶川MS8.0地震余震震源机制时空分布特征 [J]. 地球物理学报，55(4): 1213—1227.

YI Gui-xi，LONG Feng，ZHANG Zhi-wei. 2012. Spatial and temporal variation of focal mechanisms for aftershocks of the 2008 MS8.0 Wenchuan earthquake [J]. Chinese Journal of Geophysics，55(4): 1213—1227(in Chinese).

张培震，徐锡伟，闻学泽，等. 2008. 2008年汶川8.0级地震发震断裂的滑动速率、复发周期和构造成因 [J]. 地球物理学报，51(4): 1066—1073.

ZHANG Pei-zhen，XU Xi-wei，WEN Xue-ze，et al. 2008. Slip rates and recurrence intervals of the Longmen Shan active fault zone，and tectonic implications for the mechanism of the May 12 Wenchuan earthquake，2008，Sichuan，China [J]. Chinese Journal of Geophysics，51(4): 1066—1073(in Chinese).

张勇，冯万鹏，许力生，等. 2008. 2008年汶川大地震的时空破裂过程 [J]. 中国科学 (D辑)，38(10): 1186—1194.

ZHANG Yong，FENG Wan-peng，XU Li-sheng，et al. 2009. Spatio-temporal rupture process of the 2008 great Wenchuan earthquake [J]. Science in China(Ser D)，52(2): 145—154.

张勇，许力生，陈运泰. 2009. 2008年汶川大地震震源机制的时空变化 [J]. 地球物理学报，52(2): 379—389.

ZHANG Yong，XU Li-sheng，CHEN Yun-tai. 2009. Spatio-temporal variation of the source mechanism of the 2008 great Wenchuan earthquake [J]. Chinese Journal of Geophysics，52(2): 379—389(in Chinese).

郑勇，马宏生，吕坚，等. 2009. 汶川地震强余震(MS≥5.6)的震源机制解及其与发震构造的关系 [J]. 中国科学 (D辑)，39(4): 413— 426.

ZHENG Yong，MA Hong-sheng，LÜ Jian，et al. 2009. Source mechanism of strong aftershocks(MS≥5.6)of the 2008/05/12 Wenchuan earthquake and the implication for seismotectonics [J]. Science in China(Ser D)，52(6): 739—753.

Burchfiel B C，Royden L H，Van Der Hilst R D，et al. 2008. A geological and geophysical context for the Wenchuan earthquake of 12may 2008，Sichuan，People’s Republic of China [J]. GSA Today，18(7): 4—11.

Chen C，Yu C Q，Tao K，et al. 2011. Spatial distribution and focal mechanism solutions of the Wenchuan earthquake series: results and implications [J]. Earthquake Science，24(1): 115—125.

Frohlich C. 1992. Triangle diagrams: ternary graphs to display similarity and diversity of earthquake focal mechanisms [J]. Physics of the Earth and Planetary Interiors，75(1-3): 193—198. doi: 10.1016/0031-9201(92)90130-N.

Gephart J W，Forsyth D W. 1984. An improved method for determining the regional stress tensor using earthquake focal mechanism data: application to the San Fernando earthquake sequence [J]. Journal of Geophysical Research，89(B11): 9305—9320.

Hardebeck J L，Michael A J. 2006. Damped regional-scale stress inversions: methodology and examples for southern California and the Coalinga aftershock sequence [J]. Journal of Geophysical Research，111(B11): B11310. doi: 10.1029/2005JB004144.

Hardebeck J L，Shearer P M. 2002. A new method for determining first-motion focal mechanisms [J]. Bulletin of the Seismological Society of America，92(6): 2264—2276.

Hardebeck J L，Shearer P M. 2003. Using S/P amplitude ratios to constrain the focal mechanisms of small earthquakes [J]. Bulletin of the Seismological Society of America，93(6): 2434—2444.

Long F，Yi G X，Wen X Z，et al. 2012. Spatio-temporal variation of the stress field in the Wenchuan aftershock region [J]. Earthquake Science，25(5-6): 517—526.

Lund B，Townend J. 2007. Calculating horizontal stress orientations with full or partial knowledge of the tectonic stress tensor [J]. Geophysical Journal International，170(3): 1328—1335. doi: 10.1111/j.1365-246X.2007.03468.x.

Michael A J. 1984. Determination of stress from slip data: faults and folds [J]. Journal of Geophysical Research，89(B13): 11517—11526.

Michael A J. 1987. Use of focal mechanisms to determine stress: A control study [J]. Journal of Geophysical Research，92(B1): 357—368.

Shen Z K，Sun J B，Zhang P Z，et al. 2009. Slip maxima at fault junctions and rupturing of barriers during the 2008 Wenchuan earthquake [J]. Nature Geoscience，2(10): 718—724.

Yuan T，Ning J Y，Yu C Q，et al. 2013. Focal Mechanism Solutions of the 2008 Wenchuan earthquake sequence from P-wave polarities and SH/P amplitude ratios: new results and implications [J]. Earthquake Science，26(6): 357—372.

Zhang P Z，Wen X Z，Shen Z K，et al. 2010. Oblique，high-angle listric-reverse faulting and associated development of strain: The Wenchuan earthquake of May 12，2008，Sichuan，China [J]. Annual Review of Earth and Planetary Sciences，38:353—382.

Zhang X，Zha X J，Dai Z Y. 2015. Stress changes induced by the 2008 MW7.9 Wenchuan earthquake，China [J]. Journal of Asian Earth Sciences，98:98—104.