地震前震源机制解一致性参数空间动态特征分析.docxVIP

地震前震源机制解一致性参数空间动态特征分析 0 震区周边地区地震波场变化 结果表明，强震活动受到结构应力场的限制，源机制通常与应力场的方向一致（李钦祖等人，1982）。运用大量中小地震震源机制解既可以描述震源区应力场状态,也能够探讨大范围构造应力场的作用,从动力成因和运动过程两方面分析大震区数十年来应力场演化、不同阶段的特征。前苏联学者Sadovsky等(1972)报道过加尔姆地区中强地震前一年半到两年时间内,小震P轴取向由紊乱趋于一致。陈颙(1980)根据海城前震序列的研究,提出用震源机制的一致性作为判别前震的新参数;刁桂苓和于新昌(1980,1982)发现在唐山7.8级地震和北京延庆5.4 级地震前,由震区附近小震P波初动资料得到的综合断层面解矛盾符号占全部符号的比例降低的现象;刁桂苓等(1994)通过分析夏威夷Kaoiki地区的震源机制变化,提出考虑震源区应力场的改变进行地震预报的思路。刁桂苓等(2004)对大同地震序列,赵英萍等(2004)对张北地震序列震源机制一致性参数进行时间扫描发现,在强余震前一致性参数都出现低值异常。美国南加州地区两个M7.0大震(Landers和Hector Mine地震)前,该区小震震源机制存在趋于主震震源机制的现象(万永革,2008)。张致伟等(2010)反演汶川地震序列829个震源机制,发现余震区南、北段的震源机制解一致性参数随时间变化总体平稳。 2006年7月4日11时56分,在河北文安(38°53.85′N,116°16.88′E)发生5.1级地震,震源深度19.3 km(刁桂苓等,2006)。文安地震发生在河北平原强震构造带上,从公元777年至1977年该构造带内发生6级以上强震19次(中国地震局震害防御司,1995,1999)。在强震构造带的河间至宁河之间,是6级以上地震缺失时间最长的一段,能量释放最不充分,文安位于缺震段的中南部。文安地震前,2001年10月～2006年6月形成ML≥1.5小震空区,计算该空区边缘区域ML≥2.0地震视应力,发现在2004年出现明显的视应力升高过程(陈学忠,李艳娥,2007);地震波速比在文安和唐山附近地区出现正常—降低—恢复—发震的异常演变过程(王林瑛等,2008)。本文利用平面插值的方法对文安地震前后震源机制解一致性参数进行平面插值扫描,讨论地震前后一致性参数空间动态的变化特征,提取有预测意义的信息。 1 震源机制解的稳定性 本研究利用河北省遥测台网记录的河北及邻区ML≥2.0地震的数字化波形数据,并补充周边省份波形数据,震源位置取自中国地震台网中心汇编C01目录,震源机制计算采用梁尚鸿等(1984)提出的利用区域地震台网的垂直向PˉˉˉΡˉ、SˉSˉ最大振幅比资料测定小震震源参数的方法。胡新亮等(2004)利用首都圈记录到的数字化波形反演震源机制,并与格点尝试法得到的震源机制进行比较,二者基本一致,证实了该方法的可靠性。付虹等(2009)通过对云南地区19个地震的实际对比分析,发现P波和S波最大振幅比值与地震的记录形式无关,与振幅的绝对大小无关,与震级大小也基本无关。 为了保证单个震源机制解的稳定性,参与计算的台站数目一般要求6个以上,而且在震中周围均匀分布,震中距控制在150 km之内,避免出现Pn和P11混淆的情况;记录振幅要求大于0.5 mm,而且不能限幅。从图1中的台网分布来看,上述条件是容易达到的。反演计算采用的速度结构模型(表1)综合了人工地震测深和天然地震层析成像的结果(孙若昧,刘福田,1995;张成科,张先康,1997;王椿镛等,1993)。计算得到从2002年1月～2008年4月河北及邻区的震源机制解共673个,覆盖了113.0°～119.0° E,36.0°～42.0°N的范围(图1)。 一个区域的应力场可分为区域构造应力场和震源应力场。区域构造应力场主要反映整个区域内的应力状态,对区域构造应力场主要研究是否存在一个较为稳定的空间取向,认为这个稳定的空间取向对中强地震的发生起控制作用;震源应力场反映震源区附近短时间内的应力场变化、扰动,通过研究地震震源应力场取向的变化,就能发现构造运动的情况及强震的前兆信息。如果背景构造应力场的空间方向由3个相互正交的应力主轴σ1,σ2,σ3表示,震源机制的取向由3个相互正交的应力轴P、B、T表示(图2),那么就可以定义震源机制和构造应力场的一致性参数(王俊国,刁桂苓,2005):a=α+β+γ,其中:α=∠σ1P,β=∠σ2B,γ=∠σ3T。一致性参数的最大值为270°,最小值为0°,它反映了单个震源机制解的力轴与平均应力张量的差异程度,强震前孕震区内大量中小地震的震源机制一致性参数反映了孕震区内地震震源应力场与构造区域应力场的关系。 2 震源机制的一致性参数 对震源机制解的各个参数