论文部分内容阅读
2000年后中国大陆及其周边发生的大震,如2001年昆仑山地震、2004年苏门答腊地震、2008年汶川地震、2011年东日本大地震和2015年尼泊尔地震,均对中国内陆产生了持续性的应力扰动。大地震引起的应力场重新分配会诱发或抑制区域地震活动。因此,对这些影响范围大、持续时间长的大地震的同震及震后形变机制的精细研究对评估地震风险具有重要意义。汶川地震是发生在青藏高原东边界龙门山断裂的一次罕见的板内逆冲斜滑型特大地震,其同震和持续的震后松弛造成的应力场扰动,可能提升震中周缘断层上地震发生的概率,也为探测断层及岩石圈流变属性提供了难得的机会。为确定汶川地震震后松弛造成的持续影响,本文反演了汶川地下黏滞结构,并重点分析了汶川震后形变特征,同时还研究了震后余滑的影响。本文完成的主要工作如下:（1）实现从GPS坐标时间序列提取震间、同震、震后等构造信息以及非地震阶跃和季节性信号的自适应批处理程序Pytsfit（python）。该程序还可以对震间速度场、阶跃及季节性信号予以先验紧约束,提取震后瞬态形变。同时,还集成了图形绘制功能。（2）以网络工程1999-2007年速度场为基础,应用Greenspline算法内插得到汶川地震震区周围陆态网络及四川CORS站的震前速度场,最终获得扣除震间形变的GPS震后时间序列。震后形变信息显示汶川震后7年仍在持续松弛,震后两年松弛速度迅速衰减,之后松弛速度缓慢衰减。川西高原震后形变明显大于四川盆地,川西高原近场形变量达到100mm量级,而四川盆地量级仅10mm左右。虽然两盘震后形变量级差异明显,但两盘震后形变方向总体上与同震方向相近,呈逆冲兼右旋走滑状。（3）采用震后黏弹性松弛模拟程序PSGRN/PSCMP,分别对川西高原和四川盆地不同时间段、不同空间分布的GPS数据进行模拟,约束岩石流变参数。最终,本文得到川西高原低速层黏滞系数为2.51×101 8Pa s,而中下地壳的黏滞系数为3.98×101 8Pa s。采用四川盆地20082015期间的震后形变,约束得到四川盆地黏弹性层黏滞系数最优值为1.0×1020Pa s,相应弹性层深度为40km。黏弹性松弛模型在近场拟合较差,仅考虑黏弹性松弛单一震后机制不足以解释汶川震后累积形变。此外,不同的同震破裂模型模拟的震后黏弹性松弛形变仅对断层近场GPS测站有一定影响。（4）最优地球结构的地震形变模拟值系统性低于川西高原西南部地区以及鲜水河断裂带沿线的GPS观测值,系统性高于川西高原东北端方向的GPS观测值。本文认为川西高原本身存在地球介质物性的横向差异性,如果将川西高原分成西南部、中部和东北部三部分,川西高原以西壳幔黏滞系数可能低于中部黏滞系数,而东北部黏滞系数可能高于中部黏滞系数。（5）对Wang等（2011）同震破裂断层进行几何拓展,在北川断裂带北段增加北西倾向的滑脱带。以此为基础,采用弹性位错理论反演了20102015年期间的震后余滑分布。考虑到可用于震后余滑约束的GPS观测点稀疏以及同震滑动与震后余滑的空间互补关系,对同震滑动区域进行惩罚约束。论文采用“松约束”（同震滑动大于2m区域约束为0.1m）和“紧约束”（同震滑动大于6m区域约束为0.1m）两种约束方案,分别反演了断层面余滑空间展布。其中,“松约束”相较“紧约束”拟合误差更小,显示汶川同震破裂2m以上区域（破裂凹凸体）亦可能发生震后余滑,这说明WQ模型反演得到的凹凸体区域可能并不存在;（6）震后余滑产生的震后累积形变拟合汶川震后近场形变场较好,却对于中远场的数据拟合较差——远小于实际观测值,进一步说明汶川震后形变是多种震后机制共同作用的结果。且由于震后余滑模拟值也系统性低于龙日坝断裂带以西地区以及鲜水河断裂带沿线的GPS观测值,因此可以进一步确定川西高原西南侧存在明显的地球介质物性差异的可能性。