地球物理学报  2010, Vol. 53 Issue (9): 2091-2101   PDF    
四川紫坪铺水库库区地震剪切波分裂研究
张永久1 , 高原2 , 石玉涛2 , 太龄雪2     
1. 四川省地震局, 成都 610041;
2. 中国地震局地震预测研究所, 北京 100036
摘要: 本研究利用四川紫坪铺水库数字地震台网2004年8月17日~2008年5月11日的地震观测波形资料, 使用剪切波分裂SAM系统分析方法, 获得了四川紫坪铺水库库区8个数字地震台站的快剪切波偏振结果.结果表明, 紫坪铺库区台站的快剪切波偏振优势方向主要为NE或NW方向; 台站的快剪切波偏振优势方向与区域主压应力方向或活动断裂走向一致; 快剪切波偏振方向变化可能与汶川大地震前区域应力场的增加和龙门山断裂带微破裂增加有关, 慢剪切波时间延迟的变化与四川紫坪铺水库水位升降变化相关.
关键词: 四川紫坪铺水库      剪切波分裂      快剪切波偏振      主压应力      活动断裂     
The shear-wave splitting study of Sichuan Zipingpu reservoir region
ZHANG Yong-Jiu1, GAO Yuan2, SHI Yu-Tao2, TAI Ling-Xue2     
1. Earthquake Administration of Sichuan Province, Chengdu 610041, China;
2. Institute of Earthquake Science, China Earthquake Administration, Beijing 100036, China
Abstract: Using seismic data recorded by Zipingpu reservoir Digital Seismic Network from Aug.17, 2004 to May 11, 2008, the dominant polarization directions of fast shear-waves are obtained at 8 digital seismic stations in this region by SAM technique.The results show that dominant directions of polarizations of fast shear-waves are mainly in nearly NE or NW direction in Zipingpu reservoir region. The dominant polarization directions of fast shear-waves are consistent with the strikes of active faults, or regional principal compressive stress.The direction variety of polarization of fast shear-wave are related with the enhancement of regional stress and micro-cracking of Longmenshan fault before Wenchuan earthquake, and there are correlation between delay time of slow shear-wave and water level of Zipingpu reservoir..
Key words: Sichuan Zipingpu reservoir      Shear-wave splitting      Polarization of fast shear-wave      Principal compressive stress      Active faults     
1 引言

地壳介质的地震各向异性是一种比较常见的现象[1],利用地壳中的剪切波分裂可以研究地壳介质的地震各向异性特性,分析地壳应力场的状态,对有关各向异性参数进行静态和动态的描述[2].剪切波分裂对各向异性比较敏感,国内已有诸多学者利用剪切波分裂来研究地壳介质的各向异性[3~5],结果显示,复杂地质构造会造成剪切波偏振方向的不同[2, 6].由于剪切波偏振方向反映了地震台站下方地壳的主压应力方向,因而可以用来研究地壳应力场特征[7].Gao和Crampin的研究表明[8],剪切波分裂参数可以反映震前的应力积累和临震前的应力释放过程,并可用于研究地震预测问题.最新的研究还发现,慢剪切波的时间延迟对裂隙的水压变化具有很好的反映[9].

紫坪铺水库位于都江堰市区北西9km处的岷江上游.水库是坝高156 m总库容11.12亿m3的大型水利工程.紫坪铺水利枢纽建有包括7个子台的数字遥测地震台网进行实时监测,同时四川区域台网的油榨坪台也位于库区范围内.自2004年8月紫坪铺水库台网投入试运行以来,除了个别被盗子台一段时间无数据外,记录资料基本完整.

与大桥水库、二滩水库蓄水后的地震活动相比,紫坪铺水库的弱震活动相对活跃[10, 11],尤其是2005年9月下闸蓄水后在2006年10月和2008年2月发生的两次震群活动表现更为明显.因此,四川紫坪铺水库台网从2004年8月17日开始运行到2008年汶川大地震前记录了丰富的中小地震资料.本文旨在利用这些数字地震波形资料,对四川紫坪铺水库库区在汶川大地震前剪切波分裂特征进行初步研究,并探讨水库蓄放水对剪切波分裂特征的影响.

2 地质构造背景、台站及库区地震活动概况

位于扬子准地台与松潘-甘孜地槽之间的构造过渡带是中国大陆以地台为主的稳定区和中国西部以地槽为主的活动区的分界带,是新生带以来中国西部强烈的褶皱隆起区之一,构造运动导致推覆逆掩的巨型断裂系,它就是龙门山推覆构造带[12].

四川紫坪铺水库库区位于龙门山断裂带中段.龙门山断裂由一系列压性、压扭性断裂及褶皱组成,断裂带深部有北北东向地壳厚度陡变带或重力梯度带,浅部由多条断层及强烈的构造变形带组成,属控制青藏高原东缘的I级活动地块边界构造带,部分断裂全新世活动性强烈.断裂以青川-茂汶断裂、北川-保兴断裂、江油-灌县断裂规模较大,断裂总体走向N40°E左右,倾向NW,倾角50°~80°[13].在2008年5月12日汶川8.0级大地震之前,龙门山断裂带地震活动强度不高,历史上最大地震震级为6.5级,地震分布呈NE向,与龙门山断裂带走向一致.龙门山断裂带位于重力梯度带上,也是地壳厚度的分界线,其西侧地壳厚度为60~70km,东侧在50km以下.根据地震震源机制解资料得到龙门山断裂带所在地区的现代应力场的主压应力轴方位为北西西向[14],水压致裂地应力的测量结果显示龙门山断裂最大水平主应力在各段略有差异,总体上呈NW-SE的方向[15].

本研究使用了两个台网的波形数据资料,即四川紫坪铺水库台网和四川区域台网.四川紫坪铺水库数字化地震台网从2004年8月开始运行,紫坪铺水库台网共包括八角(BAJ)、白岩(BAY)、桂花树(GHS)、灵隐寺(LYS)、庙子坪(MZP)、桃子坪(TZP)和钻洞子(ZDZ)7个台站,这7个台站均匀地展布在库区周边,平均台距为10km.台站使用FBS-3B型数字拾震器,配备EDAS-24L型24位数据采集器,数据采用超短波段的单工数字无线传输与DDN公用数据网接力传输、中心统一记录的方式.观测系统的动态范围优于100dB,观测频带为0.05~40Hz,采样率秒100点.采用了GPS统一授时和定位系统,数据采集的时间服务精度优于1 ms.台网的7个子台全部坐落于较完整的岩层上,其平均速度干扰水平分布在5×10-8~2× 10-7m·s-1间.除了这7个台站外,四川区域数字台网的油榨坪(YZP)台也位于紫坪铺库区范围内(图 1).

图 1 紫坪铺库区地震台站、断裂及地震(ML>1.0)分布 左边阴影处为水磨震群位置,右边阴影处为董家山震群位置. Fig. 1 Distribution of stations, faults and earthquakes (ML>l. 0) in Zipingpu regio The left shadow is Shuimo shock cluster, and right shadow is Dongjiashan shock cluster.

2008年5月12日汶川8.0级地震前,四川紫坪铺水库库坝区地震活动强度不高,自2004年水库台网运行以来,没有记录到4级以上地震.但与省内其他的水库库区地震活动比较,比如与大桥水库、二滩水库蓄水后的地震活动相比,紫坪铺水库的弱震活动相对活跃[10, 11],尤其是2005年9月下闸蓄水后分别在2006年10月和2008年2月发生的两次震群活动表现特别明显.从四川紫坪铺水库台网2004年8月到2008年5月11日记录的库区地震目录分析,库区弱震活动频度受到水库的水位大幅度升降变化的明显影响(图 2).2005年9月开始下闸蓄水后,水位从9月30日的760 m迅速上升到10月10日的821 m,在缓慢上升到12月3日的835m后,水位开始下降,地震活动出现逐渐增加的趋势;2006年8月13日开始,库水位再次出现大幅度上升,到10月7日累计上升54m,此后水位开始下降,结果10月12日出现了水磨震群(图 1中左边阴影位置),该震群距离2008年5月12日汶川8.0级地震不到5km;2007年12月14日库水位再次达到高值后开始下降,2008年2月14日又在大坝下游10km左右的地方发生了董家山震群(图 1中右边阴影位置).

图 2 四川紫坪铺库区0.1级以上地震M-T(a)、N-T(b)及水库水位变化(c)N表示地震频次数,T表示时间. Fig. 2 The M-T (a), N-T(b) and water level (c) of Zipingpu region
3 方法和数据处理

地震波在各向异性介质中传播时,由于传播速度的不同,剪切波会分裂成快、慢两个剪切波列,这两个剪切波具有近似垂直的偏振特性.地壳岩石中普遍分布有近似直立、定向排列的平行裂隙,即Extensive-dilatancy Anisotropy(EDA)裂隙[16]. EDA裂隙结构是各向异性的,剪切波在这种结构中传播时,会分裂成快剪切波和慢剪切波.快剪切波的偏振方向与应力诱发的定向排列裂隙面的走向一致:即快剪切波偏振方向与原地最大主压应力方向一致;慢剪切波的时间延迟反映了地壳介质的各向异性程度,它受到介质中裂隙的物理特性和所含的流体特性的影响.剪切波入射到自由表面时,若入射角大于临界角,则会发生全反射现象,这是剪切波窗口的限制条件.对于泊松介质,临界角约为35°,但由于地壳表面沉积层的影响,简单地按照单层介质计算,选择入射角小于45°范围内的地震射线,即可满足剪切波窗口的要求[17].

本研究使用剪切波分裂系统分析法---SAM方法[18].SAM方法是在相关函数的基础上提出的一种分析方法,包括相关函数计算、时间延迟校正和偏振分析检验3个部分,具有自我检验的特点[18, 19].新的SAM方法软件系统升级为SAM(2007)[20, 21],具体的计算原理和分析方法可参阅相关文献[2, 6, 19, 21~24].

图 3显示的分别是桃子坪台、八角台和灵隐寺台记录到的三次不同地震事件剪切波的质点运动轨迹,即剪切波偏振图.三次地震事件分别是TZP台记录的2006年1月18日1.7级地震,震中距13km,震源深度12km;BAJ台记录的2007年2月12日1.3级地震,震中距3km,震源深度6km;LYS台记录的2008年2月14日3.7级地震,震中距9km,震源深度12km.由于剪切波已经发生了分裂,在水平方向上两个分量的剪切波的质点运动轨迹并不是线性的.如果快剪切波偏振方向为α角,把剪切波按α角进行顺时针旋转,得到快、慢剪切波,再进行时间延迟校正,把慢剪切波提前,扣除时间延迟的影响,得到偏振图,此时的两列剪切波的偏振图应该表现为更接近线性的偏振.若快剪切波偏振方向的计算结果不正确,由于快、慢剪切波相互混在一起,则得不到分离的快、慢剪切波.若时间延迟的计算结果不正确,则扣除时间延迟影响的剪切波偏振图不会是线性的.如果发现结果有偏差,则需要根据波形和偏振特征进行修正,这也正是SAM方法的优势所在[6, 18, 22].

图 3 剪切波分裂分析结果 从左至右分别为TZP、BAJ、LYS.上为南北(NS)、东西(EW)分量剪切波形及质点运动轨迹,下为快(F)、慢(S)剪切波及经时间延迟校正的质点运动轨迹. Fig. 3 The results of shear-wave splitting analysi The left is TZP; the middle is BAJ; and the right is LYS. The first is shear waveforms of north-south (NS), east-west (EW) and the trail of particle of NS and EWrespectively, the second is fast shear-wave (F), slow shear-wave (S) and the trail of particle of fast and slow shear-wave, which have eliminated the effect of time delay.
4 计算结果

根据四川紫坪铺水库数字地震台网提供的地震目录,本文收集整理了四川紫坪铺水库数字地震台网监测范围内自2004年8月到2008年5月11日的1.0级以上并有深度记录的地震460次.通常震中定位误差有几千米,剪切波分裂的时间延迟在相当程度上会受到震源定位精度以及震源机制的影响.如果震源深度误差是5~12km,延时误差可达到30%~40%[25, 26].因此,本研究在进行剪切波分裂分析前运用双差定位方法[27]对所有地震进行了重新精确定位.

如前所述,剪切波窗理论入射角的临界值为35°,但射线在接近地表时,由于低速层或沉积层的影响,射线传播的轨迹会逐渐趋向垂直方向.根据龙门山地区的速度结构[28]并结合石油部门的浅层勘探速度结果,在该地区给出厚度分别为1km、5km、20km,速度为3.7km/s、4.3km/s、5.9km/s的条件.对于一个深度为12km,震中距为19.5km的地震,简单计算入射角为58.5°,但实际入射角为35. 3°.本研究中使用单层模型概念下的58°入射角作为剪切波窗口,根据这里给出的速度模型,本研究使用的数据都在剪切波窗口内.根据剪切波窗口的约束条件,得到剪切波窗口内的有效地震记录520条.在此基础上挑选观测质量较好的波形记录,最后得到用于计算剪切波分裂的地震记录226条,其中八角台70条,其余各台的记录在20~25条之间.图 4给出了8个台的剪切波分裂结果.

图 4 剪切波偏振方向等面积极射投影与等面积玫瑰图 图中的短线方向是相应台站记录的每个地震事件的快剪切波偏振方向,长短代表慢剪切波延迟时间大小. Fig. 4 Lower hemispherical project and equal-area rose diagram of fast shear-wave polarizations The directions of the short lines in the figure is each polarization of the fast shear-wave which is the station corresponded, and the length represent the delay tme level of slow shear-wave.

快剪切波偏振方向与台站的原地主压应力方向一致,受区域应力场的控制,与断裂分布和地质构造等有关[6, 19, 23, 24, 29].图 5是紫坪铺库区8个台站的快剪切波偏振方向优势分布.从图 5中可以看出,除GHS、MZP和TZP三个台快剪切波偏振优势方向仅一个方向突出外,其余的5个台站的快剪切波偏振都显示出了两个相对明显的优势方向,一个为NE向,另一个为NW向,而YZP台的快剪切波偏振优势方向为近N-S.

图 5 紫坪铺库区快剪切波偏振方向 具有两个优势方向的台站中,深色的表示优势方向Ⅰ,浅色的表示优势方向Ⅱ. Fig. 5 Polarization of fast shear-wave for Zipingpu region The dark show dominant direction 1, undertone showdominant direction Ⅱ.

由于同一台站接收到不同地震事件的快剪切波偏振方向不完全一致.因此,采用求角度平均的方法[30~33]分别计算各个台站的平均快剪切波偏振优势方向θR值及慢波时间延迟(表 1图 5),R值在0~1.0之间,其大小反映了快剪切波偏振方向的集中程度,越接近1.0说明快剪切波偏振方向越集中.

表 1 四川紫坪铺水库台网7个台站及油榨坪剪切波分裂参数 Table 1 The shear-wave split parameters of Sichuan Zipingpu reservoir network and YZP

将整个研究区域8个台共计226个剪切波分裂事件的数据综合起来,得到全部记录的快剪切波偏振方向等面积投影玫瑰图(图 6).所有剪切波分裂事件的平均快波偏振方向为126.2°±52.2°,R值为0.64,慢波的平均延迟时间为5.38±3.62 ms/km.同时,图 6也显示出整个区域的快剪切波偏振方向具有两个主要优势方向,两个优势方向分别为51.2°±24.7°和135.9°±21.1°,R值分别为0.91和0.93.即研究区域的快剪切波偏振方向一个优势是北西方向,与研究区的主压应力方向一致;另一个优势为近北东方向,与研究区里的龙门山断裂带的走向一致.

图 6 紫坪铺库区快剪切波偏振方向等面积投影玫瑰图 Fig. 6 Project and equal-area rose diagram of fast shear-wave polarizations in Zipingpu region
5 分析与讨论

赵珠等[28]的研究结果表明,龙门山及附近地区剪切波速度结构具横向不均匀性.这与龙门山断裂带位于重力梯度带上,也是地壳厚度的分界线的构造上的复杂性相关.本研究范围尽管只涉及龙门山中段的一部分,但其结果也显示出这种构造上及应力场双重影响的复杂性.

紫坪铺库区8个台站的台间距最大的不到20km,但剪切波分裂分析结果已显示出明显差异. 8个台的平均快剪切波偏振方向除GHS的与龙门山后山断裂走向一致外,其余7个台的结果没有显示出与龙门山断裂走向或区域主压应力方向的简单一致(表 1).除GHS、MZP和TZP三个台快剪切波偏振方向有很清楚的一个优势方向外,其余5个台站的快剪切波偏振都显示出了两个相对明显的优势方向.在三个有单个优势方向的台站中,TZP和MZP台的快剪切波偏振的优势方向为NW,GHS台快剪切波偏振的优势方向为NE.在5个有复杂偏振优势方向分布的台站中,BAJ、LYS、ZDZ和BAY等4个台快剪切波偏振的两个优势方向分别为NW和NE方向,而YZP台的快剪切波偏振的主要优势方向为NNW近N-S方向,而另一个较离散的方向相对集中指向约为NE(图 5).

YZP台的快剪切波偏振优势方向为近N-S向,与原地主压应力方向、GPS观测结果及龙门山断裂走向都存在较大的差异,这与原来使用更长时间、仅包含区域台网波形资料的研究结果是一致的[33].这从一个方面也证明我们的研究结果可能真实地反映了研究台站附近的应力场或断层情况.

2008年5月12日汶川8.0级地震后,石玉涛等[34]研究了汶川余震序列地震的各向异性变化特征.其研究所用台站分布于整个龙门山,其中有三个台距离本研究所用台站较近.三个台中,L105台的平均快剪切波为近N-S,和本研究YZP台结果一致;L106、L5501台的快剪切波都有两个优势方向,与本研究BAJ、LYS、ZDZ和BAY等台结果一致.进一步分析台站分布与龙门山断裂的分布,L105、YZP两个台都分布在龙门山断裂的东南侧,且距离相对远.而L106、L5501、BAJ、LYS、ZDZ和BAY等台都在龙门山断裂上或距离断裂很近.可见,剪切波分裂参数真实地反映了构造和应力场特征.

从快剪切波偏振方向随时间变化过程分析(图 7),四川紫坪铺库区的快剪切波的偏振方向在2006年10月之前比较分散,在近NS、EW、NE和NW向都有明显分布.自2006年11月开始,快剪切波偏振方向出现了新的变化趋势,近NS向的逐渐消失,近EW向的自2007年下半年开始也逐渐减少,NE和NW两个方向逐渐突出,而这两个方向一个是区域应力场方向,另一个是汶川大地震发震断裂---龙门山断裂方向.因此,四川紫坪铺库区快剪切波偏振方向随时间的逐渐集中这一变化可能与汶川大地震前区域应力场的增加和龙门山断裂带微破裂增加有关.

图 7 紫坪铺库区快剪切波偏振方向随时间的变化 Fig. 7 The variety of fast shear-wave polarizations in Zipingpu region

慢剪切波时间延迟参数是度量各向异性程度大小的物理量,它对地壳应力场的微小变化所引起的地壳裂隙的变化非常敏感.紫坪铺库区8个台站共计226条记录的平均慢剪切波时间延迟为3.92±2.72 ms/km,最大值是八角台(BAJ),两个不同优势方向的时间延迟分别为7.52±4.25 ms/km和7.51±4.52 ms/km,全部记录的平均值为7.41±4.32 ms/km.其次是灵隐寺台(LYS),两个不同优势方向的时间延迟分别为4.97±2.73 ms/km和8.72±4.11 ms/km,全部记录的平均值为6.22±3.64ms/km.最小的是油榨坪台,两个不同优势方向的时间延迟分别为4.12±1.85ms/km和2.89±1.64 ms/km,全部记录的平均值为3.42±1.80ms/km.其他5个台的时间延迟相差较小,都在4ms/km左右.引人注目的是,慢剪切波延迟时间较大的灵隐寺台和八角台正好位于四川紫坪铺水库的库坝区和库尾区,而慢剪切波延迟时间最小的油榨坪台距离四川紫坪铺水库水域最远.根据雷兴林等[35]的研究结果,灵隐寺台和八角台正好处于水库蓄放水引起的应力变化较大的区域两端,而油榨坪台所在的位置应力基本不受水位变化影响.

从慢剪切波时间延迟随时间变化过程分析,各个台站表现出了较好的一致性.由于地震记录样本数量的关系,仅对所有台站的记录和有效记录最多的八角台(表 1)慢剪切波时间延迟进行具体分析.在四川紫坪铺水库2005年9月底开始蓄水后,其水位变化出现了三次比较大的升降过程,升降幅度在50~75m之间,水位的大幅度升降变化必然导致库区应力场的同步变化.从整个库区和八角台慢剪切波时间延迟的变化过程分析,慢剪切波时间延迟与四川紫坪铺水库水位的升降变化显示出相关性(图 8),每当水位大幅度上升到最大值后,慢剪切波延迟时间随之增加,水位下降后慢剪切波延迟时间随之缩短.这个关系也基本符合Gao和Crampin[36]实验研究的剪切波分裂时间延迟变化与压力的关系.另外,整个库区和八角台的慢剪切波延迟时间数据比较离散(图 8a8b),根据APE理论,时间延迟可能会由于应力的微小变化而产生显著的变化[37],高孔隙流体压力的剧烈变化会导致时间延迟数据的发散[38].因而认为,造成整个库区和八角台时间延迟数据发散可能是高孔隙流体压力的剧烈变化所致.

图 8 库区(a)、八角台(b)慢剪切波时间延迟变化与水库水位变化(c) (a),(b)中单点线是原始计算值,实线是5点滑动滤波值. Fig. 8 The variety of slow shear-wave delay time of reservoir area (a), BAJ (b) and water level (c) The dot line is originality value, and real line is sliding filter value with five point of (a), (b).

裂隙密度与剪切波速度各向异性直接相关,是描述地壳介质物理性质的基本参数,并随着应力场的变化而变化.根据关于具有稀疏、充水平行直立排列裂隙的岩石波速异常公式[39],可以大致估计岩石裂隙密度.对于泊松介质,可以推出如下公式:

(3)

式中ρ为裂隙密度,β为不含裂隙介质的剪切波速度,τ为时间延迟.公式(3)中假设了裂隙充满传播路径,入射波与裂隙方向平行.根据(3)式计算得到紫坪铺库区的裂隙密度平均值为0.013±0.009,八角台附近和灵隐寺台附近的裂隙密度值分别为0.025±0.015和0.021±0.012.可见紫坪铺库区的平均裂隙密度与吴晶等[24]计算的京西北地区的结果比较接近,而处于四川紫坪铺水库库坝和库尾附近的灵隐寺台和八角台附近的裂隙密度比整个库区的平均值高60%以上.

6 结论

龙门山断裂带位于重力梯度带上,也是地壳厚度的分界线.四川紫坪铺水库台网位于龙门山三大主干断裂上,其地形地貌及地下结构都较为复杂.本文使用四川紫坪铺水库库区8个数字地震台站2004年8月~2008年5月11日的波形资料,研究了四川紫坪铺水库库区的剪切波分裂现象,给出了该地区各台站的快波偏振方向,各个台站数据处理的结果呈现出较好的一致性.

快剪切波偏振优势方向代表了观测区域的原地(台站及其附近)主压应力方向或观测区域活动断裂走向.本研究计算得到的四川紫坪铺水库台网以及油榨坪台共计226个剪切波分裂事件的平均快波偏振方向为126.2°±52.2°,慢波的平均延迟时间为5.38±3.62 ms/km.同时,结果显示紫坪铺库区大部分台站有两个快剪切波偏振优势方向,一个是NW方向,另一个是NE方向.这些结果充分体现了四川紫坪铺水库库区北西向的区域主压应力场和北东向的龙门山活动断裂的共同作用.

YZP台的快剪切波偏振优势方向尽管比较离散,但主优势方向为近N-S向,与原来使用更早时间记录波形资料的研究结果是一致的[33],与库区主压应力方向及龙门山断裂走向都存在较大的差异,与紫坪铺库区其他7个台站都有两个优势方向的结果也不一致.其他学者[4, 5, 23, 24]的研究结果也显示,根据快剪切波偏振优势方向的分布不仅可以判断地下断层是否为活动断层,而且还可能揭示新的构造关系,比如在一定深度下存在的隐伏断层等.因此,YZP的剪切波分裂结果可能真实地反映了该台附近的应力场或断层情况.

其他一些剪切波分裂研究结果[2, 22, 29, 36, 40~45]显示,强震前后剪切波分裂参数会出现变化.本文的研究结果表明,库区快剪切波偏振方向随时间的变化可能与汶川大地震前区域应力场的增加和龙门山断裂带微破裂增加有关.慢剪切波时间延迟分析结果显示,紫坪铺库区8个台站中,位于库坝附近的LYS台和位于库尾附近的BAJ台处于水库蓄放水引起的应力变化较大的区域两端,其平均慢波时间延迟明显大于其他台站.而整个库区和BAJ台的慢波时间延迟与水库水位的变化相关.根据时间延迟得到的处于四川紫坪铺水库库坝和库尾附近的LYS台和BAJ台附近的裂隙密度比整个库区的平均值高60%以上.

本文的研究结果仅使用了汶川地震前的资料,得到的结果也是初步的.2008年5月12日汶川8.0级地震后,石玉涛等[34]利用汶川地震序列资料进行了剪切波分裂研究,但其研究范围远超出了四川紫坪铺水库库区的范围,而且汶川地震使得四川紫坪铺水库台网的多个台站损毁,他们的研究仅使用了区域台网及震后流动台站资料,而没有使用四川紫坪铺水库台网资料.尽管本文与这些研究的具体范围和使用台站资料都不同,但对于相同区域的研究结果是有一致性的.

致谢

感谢四川省地震局水库研究所和监测研究所提供各种帮助和支持.

参考文献
[1] Crampin S. An introduction to wave propagation in anisotropic media. Geophys. J. R. astr. Soc. , 1984, 76: 17-28. DOI:10.1111/j.1365-246X.1984.tb05018.x
[2] 高原, 郑斯华, 周蕙兰. 唐山地区快剪切波偏振图像及其变化. 地球物理学报 , 1999, 42(2): 228–232. Gao Y, Zheng S H, Zhou H L. Polarization patterns of fast shear wave in Tangshan region and their variations. Chinese J. Geophys. (in Chinese) , 1999, 42(2): 228-232.
[3] 姚陈, 王培德, 陈运泰. 卢龙地区S波偏振与上地壳裂隙各向异性. 地球物理学报 , 1992, 35(3): 305–315. Yao C, Wang P D, Chen Y T. Shear-wave polarization and crack induced anisotropy of upper crust in Lulong, north China. Chinese J. Geophys. (in Chinese) , 1992, 35(3): 305-315.
[4] 赖院根, 刘启元, 陈九辉, 等. 新疆伽师强震群区的横波分裂与应力场特征. 地球物理学报 , 2002, 45(1): 83–92. Lai Y G, Liu Q Y, Chen J H, et al. Features of the shear-wave splitting and stress field in the Xinjiang-Jiashi strong earthquake region. Chinese J. Geophys. (in Chinese) , 2002, 45(1): 83-92.
[5] 赖院根, 刘启元, 陈九辉, 等. 首都圈地区横波分裂与地壳应力场特征. 地球物理学报 , 2006, 49(1): 189–196. Lai Y G, Liu Q Y, Chen J H, et al. Shear wave splitting and the features of the crustal stress field in the Capital circle. Chinese J. Geophys. (in Chinese) , 2006, 49(1): 189-196.
[6] 高原, 郑斯华, 孙勇. 唐山地区地壳裂隙各向异性. 地震学报 , 1995, 17(3): 283–293. Gao Y, Zheng S H, Sun Y. Crustal crack anisotropy in Tangshan. Acta Seismologica Sinica (in Chinese) , 1995, 17(3): 283-293.
[7] Crampin S. Seismic-wave propagation through a cracked solid:polarization as a possible dilatancy diagnostic. Geophys. J. R. astr. Soc. , 1978, 53: 467-496. DOI:10.1111/gji.1978.53.issue-3
[8] Gao Y, Crampin S. Observations of stress relaxation before earthquakes. Geophys. J. Int. , 2004, 157(2): 578-582. DOI:10.1111/gji.2004.157.issue-2
[9] Tang C, Rial J A, Lees J M. Shear-wave splitting:a diagnostic tool to monitor fluid pressure in geothermal fields. Geophys. Res. Lett. , 2005, 32: L21317. DOI:10.1029/2005GL023551
[10] 胡先明, 谢蓉华. 紫坪铺水库地震台网记录地震的复杂性. 四川地震 , 2006(2): 17–20. Hu X M, Xie R H. Seismic recording complexity of Zipingpu Dam Telemetered Seismic Network. Earthquake Research in Sichuan (in Chinese) , 2006(2): 17-20.
[11] 胡先明. 紫坪铺水库蓄水前天然地震活动. 四川地震 , 2007(2): 16–21. Hu X M. Seismicity background around Zipingpu Dam before impounding. Earthquake Research in Sichuan (in Chinese) , 2007(2): 16-21.
[12] 邓起东, 张培震. 中国活动构造基本特征. 中国科学(D辑) , 2002, 32(12): 1020–1030. Deng Q D, Zhang P Z. Basic characteristics of active tectonics of China. Science in China (Ser. D) (in Chinese) , 2002, 32(12): 1020-1030.
[13] 唐荣昌, 文德华, 黄祖智, 等. 松潘-龙门山地区主要活动断裂带第四纪活动特征. 中国地震 , 1991, 7(3): 64–71. Tang R C, Wen D H, Huang Z Z, et al. The quaternary activity characteristics of several major active faults in the Songpan-Longmenshan region. Earthquake Research in China (in Chinese) , 1991, 7(3): 64-71.
[14] 邓起东, 陈社发, 赵小麟. 龙门山及其邻区的构造和地震活动及动力学. 地震地质 , 1994, 16(4): 389–403. Deng Q D, Chen S F, Zhao X L. Tectonics seismicity and dynamics of Longmenshan Mountains and its adjacent regions. Seismology and Geology (in Chinese) , 1994, 16(4): 389-403.
[15] 安其美, 丁立丰, 王海忠, 等. 龙门山断裂带的性质与活动性研究. 大地测量与地球动力学 , 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.
[16] Crampin S, Atkinson B K. Microcracks in the earth's crust. First Break , 1985, 3(3): 16-20.
[17] Crampin S. A review of shear-wave splitting in the compliant crack-critical anisotropic earth. Wave Motion , 2005, 41: 59-77. DOI:10.1016/j.wavemoti.2004.05.006
[18] 高原, 刘希强, 梁维, 等. 剪切波系统分裂方法(SAM)软件系统. 中国地震 , 2004, 20(1): 101–107. Gao Y, Liu X Q, Liang W, et al. Systematic analysis method of shear wave splitting:SAM software system. Earthquake Research in China (in Chinese) , 2004, 20(1): 101-107.
[19] 高原, 郑斯华. 唐山地区剪切波分裂研究(Ⅱ)--相关函数分析法. 中国地震 , 1994, 10(Suppl.): 11–21. Gao Y, Zheng S H. On shear wave splitting in Tangshan region (Ⅱ)-Correlation Function Analysis Method. Earthquake Research in China (in Chinese) , 1994, 10(Suppl.): 11-21.
[20] 高原, 石玉涛, 梁维, 等. 剪切波分裂分析系统SAM (2007)--软件系统. 中国地震 , 2008, 24(4): 345–353. Gao Y, Shi Y T, Liang W, et al. Systematic analysis method of shear-wave splitting SAM (2007):software system. Earthquake Research in China (in Chinese) , 2008, 24(4): 345-353.
[21] 石玉涛, 高原, 吴晶, 等. 剪切波分裂分析系统SAM (2007)--区域地震台网资料应用实例. 中国地震 , 2008, 24(4): 354–361. Shi Y T, Gao Y, Wu J, et al. Systematic analysis method of shear-wave splitting SAM (2007):examples of application in the regional seismograph networks. Earthquake Research in China (in Chinese) , 2008, 24(4): 354-361.
[22] 高原, 郑斯华, 王培德. 海南省东方地区1992年小震群剪切波分裂研究. 地球物理学报 , 1996, 39(2): 221–232. Gao Y, Zheng S H, Wang P D. Shear wave splitting of small earthquakes in 1992 in Hainan Dongfang area. Chinese J. Geophys. (in Chinese) , 1996, 39(2): 221-232.
[23] 石玉涛, 高原, 吴晶, 等. 云南地区地壳介质各向异性--快剪切波偏振特性. 地震学报 , 2006, 28(6): 574–585. Shi Y T, Gao Y, Wu J, et al. Seismic anisotropy in crust in Yunnan-preliminary study on polarization of fast shear wave. Acta Seismologica Sinica (in Chinese) , 2006, 28(6): 574-585.
[24] 吴晶, 高原, 陈运泰, 等. 首都圈西北部地区地壳介质地震各向异性特征初步研究. 地球物理学报 , 2007, 50(1): 209–220. Wu J, Gao Y, Chen Y T, et al. Seismic anisotropy in the crust in northwestern capital area of China. Chinese J. Geophys. (in Chinese) , 2007, 50(1): 209-220.
[25] Crampin S, Volti T, Stefansson R. A successfully stress forecast earthquake. Geophys. J. Int. , 1999, 138: FI-F5.
[26] 张尉, 戴恒昌, 李向阳, 等. 水库气枪震源产生的S波及其分裂. 地震地质 , 2009, 31(1): 22–33. Zhang W, Dai H C, Li X Y, et al. Shear wave generated by airguns in a reservoir and its splitting. Seismology and Geology (in Chinese) , 2009, 31(1): 22-33.
[27] Waldhauser F.HypeDD-A program to compute double-difference hypocenter locations.US Geological Survey Open-File Report, 2001.113
[28] 赵珠, 范军, 郑斯华, 等. 龙门山断裂带地壳速度结构和震源位置的精确修定. 地震学报 , 1997, 19(6): 615–622. Zhao Z, Fan J, Zheng S H, et al. Precision determination of the crustal structure and hypocentral locations in the Longmenshan thrust nappe belt. Acta Seismologica Sinica (in Chinese) , 1997, 19(6): 615-622.
[29] Gao Y, Crampin S. A further stress-forecast earthquake (with hindsight), where migration of source earthquakes causes anomalies in shear-wave polarizations. Tectonophysics , 2006, 426: 253-262. DOI:10.1016/j.tecto.2006.07.013
[30] Davis J C. Statistics and Data Analysis in Geology. Hoboken N J: John Wiley, 1986 : 649 .
[31] Mardia K V, Jupp P E. Directional Statistics. Hoboken N J: John Wiley, 2000 : 429 .
[32] 马宏生, 张国民, 刘杰, 等. 云南及邻区应力应变场分区耦合特性初步研究. 地震学报 , 2007, 29(2): 130–141. Ma H S, Zhang G M, Liu J, et al. Coupling characteristics of stress and strain at different layers of different sub-regions in Yunnan and its adjacent areas. Acta Seismologica Sinica (in Chinese) , 2007, 29(2): 130-141.
[33] 张永久, 高原, 石玉涛, 等. 四川区域地震台网的剪切波分裂研究. 地震学报 , 2008, 30(2): 123–134. Zhang Y J, Gao Y, Shi Y T, et al. A study on S wave splitting using waveform data from Sichuan Regional Seismic Network. Acta Seismologica Sinica (in Chinese) , 2008, 30(2): 123-134.
[34] 石玉涛, 高原, 赵翠萍, 等. 汶川地震余震序列的地震各向异性. 地球物理学报 , 2009, 52(2): 398–407. Shi Y T, Gao Y, Zhao C P, et al. A study of seismic anisotropy of Wenchuan earthquake sequence. Chinese J.Geophys. (in Chinese) , 2009, 52(2): 398-407.
[35] 雷兴林, 马胜利, 闻学泽, 等. 地表水体对断层应力与地震时空分布影响的综合分析--以紫坪铺水库为例. 地震地质 , 2008, 30(4): 1046–1064. Lei X L, Ma S L, Wen X Z, et al. Integrated analysis of stress and regional seismicity by surface loading-a case study of Zipingpu reservoir. Seismology and Geology (in Chinese) , 2008, 30(4): 1046-1064.
[36] Gao Y, Crampin S. Temporal variations of shear-wave splitting in field and laboratory studies in China. J. Appl. Geophys. , 2003, 54: 279-287. DOI:10.1016/j.jappgeo.2003.01.002
[37] Crampin S, Zatsepin S V. Modelling the compliance of crustal rock-Ⅱ. Response to temporal changes before earthquakes. Geophys. J. Int. , 1997, 129: 495-506.
[38] Crampin S, Peacock S, Gao Y, et al. The scatter of time-delays in shear-wave splitting above small earthquakes. Geophys. J. Int. , 2004, 156(1): 39-44. DOI:10.1111/gji.2004.156.issue-1
[39] Hudson J A. Wave speeds and attenuation of elastic waves in material containing cracks. Geophys.J.R.astr.Soc. , 1981, 64: 133-150. DOI:10.1111/j.1365-246X.1981.tb02662.x
[40] 郑治真. S波分裂在地震趋势估计中的应用. 地震研究 , 1991, 14(1): 1–8. Zheng Z Z. The application of shear wave splitting in earthquake tendency evaluation. Journal of Seismologica Research (in Chinese) , 1991, 14(1): 1-8.
[41] 郑秀芬, 陈朝辉, 张春贺. 区域数字地震台网开展近场源剪切波分裂的实用化研究. 地震地磁观测与研究 , 2006, 27(2): 1–3. Zheng X F, Chen C H, Zhang C H. The introduction to the practical near-source shear wave splitting study in regional digital seismograph network. Seismological and Geomagnetic Observation and Research (in Chinese) , 2006, 27(2): 1-3.
[42] 郑秀芬, 陈朝辉, 张春贺. 1999年台湾集集地震余震区--嘉义地区地震的剪切波分裂参数随时间变化的研究. 地球物理学报 , 2008, 51(1): 149–157. Zheng X F, Chen C H, Zhang C H. Study on temporal variations of shear wave splitting in the Chiayi area, aftershock zone of 1999 Chichi earthquake, Taiwan. Chinese J. Geophys. (in Chinese) , 2008, 51(1): 149-157.
[43] Crampin S, Both D C, Evans R, et al. Comment on "Quantitative measurements of shear wave potarizations at the Anza seismic network, southern California: implications for shear wave splitting and earthquake prediction, by R. C. Aster, P.M. Shearer and J. Berger".J. Geophys. Res. , 1991, 96: 6403-6414.
[44] Liu Y, Crampin S, Main L. Shear-wave anisotropy:spatial and temporal variations in time delays at Parkfield, Central California. Geo phys. J. Int. , 1997, 130: 771-785. DOI:10.1111/gji.1997.130.issue-3
[45] 陈天长, 李志君, 张润生, 等. 1989年6月9日石棉5.2级地震后横波分裂的时间变化. 地震学报 , 1995, 17(1): 52–61. Chen T C, Li Z J, Zhang R S, et al. The time variety of transverse wave splitting after Shimian 5.2 earthquake on 6/9, 1989. Acta SeismologicaS inica (in Chinese) , 1995, 17(1): 52-61.