GPS地壳形变观测及其在中亚大三角地震构造域的应用

肖根如

（中国地震局地质研究所，北京100029）

GPS地壳形变观测及其在中亚大三角地震构造域的应用

肖根如

（中国地震局地质研究所，北京100029）

纵观全球中强地震的空间分布，除一系列的洋脊、裂谷、海沟、转换断层和大陆内部的古板块边缘等大型构造活动带被清晰地勾画和显示外，在亚洲大陆东部地区相对密集的地震还清楚地勾画了一个面状的三角形区域。自1898年以来，该区域记录到的6级及以上地震占到全球地震的4.4%。在如此小的范围内发生如此多的强地震，前人将这一特殊区域称之为“中亚大三角地震构造域”。中亚大三角地震构造域位于亚洲的中东部地区，地处60°E～120°E，10°N～60°N的矩形范围内，总体形态大致呈正三角形，其西南边界为平行于喜马拉雅山脉的区域，是大三角形的底边；东边界大致沿东经105°线，即从中国境内的南北地震带沿华南地块、鄂尔多斯地块西边缘向北延伸至贝加尔湖；西北边界则从帕米尔高原开始，穿越天山、阿尔泰山至贝加尔湖地区。

自1900年以来，中亚大三角地震构造域内共发生过121次7级以上大地震。其中8级以上大地震有22次，这些大地震主要发生在青藏高原、天山造山带以及贝加尔断裂带上。最近十多年来发生的破坏性地震有：1997年11月8日玛尼MS7.9地震、2001年11月14日昆仑山MS8.1地震、2008年3月21日于田MS7.3地震、2008年5月12日汶川MS8.0地震、2010年4月14日玉树MS7.1地震等。

中亚大三角地震构造域是一个由密集的中强地震所自然呈现出的特殊区域，其内部跌宕起伏的高原、山脉、山间盆地和地堑，以及错综复杂的活动断裂直接反映了该区域强烈的地壳运动和构造形变。中亚大三角地震构造域主要包含青藏高原和西域两大构造区域。其中，青藏高原构造区域内部发育有一系列规模宏大、活动性强的弧形断裂带，主要的大型断裂带有：喜马拉雅主中央断裂带、西昆仑断裂带、喀喇昆仑—嘉黎断裂带、鲜水河断裂带、东昆仑断裂带、阿尔金断裂带和红河断裂带；西域构造区内部发育有天山断裂带、戈壁阿尔泰断裂带、蒙古—阿尔泰断裂带、Bolnay断裂带等。

本论文以高精度、高密度、大范围的GPS观测资料为主，并结合活动断裂资料和最近百余年的中强地震资料，对中亚大三角地震构造域的现今地壳形变特征以及地壳形变与地震活动性的关系进行探讨。具体研究方法与成果包括：

（1）在吸纳国内外关于高精度GPS数据处理方法和策略的基础上，完善了一套基于GIPSY和QOCA的高精度、高效率的自动化GPS数据处理系统。在USGS（United States Geological Survey）的自动化GPS处理软件包GP的基础上，吸收和采纳目前国际上GPS数据处理的最新模型与方法，利用TEQC、GIPSY、Ambizap以及QOCA等软件，完善出一套实用化的GPS前期数据处理与后期数据分析系统，实现了从GPS数据准备、GPS数据处理、联合平差及结果可视化展示等一系列步骤的自动化，实现了高效的GPS数据处理。

（2）融合国内外多渠道GPS资料，获得中亚大三角地震构造域的综合GPS速度场。主要通过两种方式获得中亚大三角地震构造域的GPS速度场：①野外流动加密观测；②将不同机构、不同区域、不同时期、不同参考框架下的GPS速度场资料，以不同数据的公共点为“桥梁”，通过球面欧拉旋转变换，融合、归化至相同的参考框架。通过一致性检核和粗差剔除等环节，最后得到1 647个GPS站点的速率。首次给出中亚大三角地震构造域及其周缘的GPS速度场。

（3）基于高精度、高密度的GPS速度场结果，并结合中亚大三角地震构造域的活动构造概要，总结出中亚大三角地震构造域GPS速度场的主要特征和区域内部的差异运动特点：

①相对于稳定的欧亚板块，中亚大三角地震构造域的北东向水平地壳运动速度的大小，从南部喜马拉雅地区的约40 mm／a向北依次梯度变小至北部贝加尔地区约3 mm／a。反映出中亚大三角在印度板块北东向俯冲、推挤作用下的缩短、隆升。

②在整个中亚大三角地震构造域，青藏高原表现出了最强烈的水平向地壳差异运动。由于塔里木块体、阿拉善块体和鄂尔多斯块体的“准刚性”阻挡，以及高原北部区域对南部区域的自然阻挡作用，青藏高原在发生北东向强烈挤压缩短的同时，产生了明显的东向挤出式“逃逸”。

③在青藏高原的东南部，由于“阿萨姆犄角”的北东向强烈插入，再加之华南“准刚性”块体在东部的阻挡，使该区域的高塑性地壳物质围绕喜马拉雅东构造结产生强烈的顺时针流滑式挤出逃逸。

④塔里木地块在整个中亚大三角地震构造域的现今地壳变形中扮演着重要的角色：由于塔里木地块的调节和转化作用，使中亚大三角地震构造域中部的水平地壳运动，呈现出从NNE到NE方向的扇状发散特征。在塔里木盆地内部，速度大小比较均匀，反映出较小的内部形变。

⑤天山地区是中亚大三角地震构造域仅次于青藏高原的强烈挤压变形区。在跨越天山前后，GPS速度场从15～20 mm／a迅速变小到1～5 mm／a。

⑥在蒙古中西部及贝加尔湖地区，站点的运动速度普遍在2～5 mm／a之间，且运动方向无明显规律性，反映出该区域相对比较微弱的地壳形变背景。

⑦在帕米尔高原地区，由于缺少足够的GPS观测站点，我们尚无法直观判定其内部的地壳运动差异性，但根据其周缘的GPS速度场状况，可推测其可能存在塔里木刚性块体阻挡下强烈的逆时针旋转。

（4）基于优选的空间内插算法，给出中亚大三角地震构造域连续分布的GPS应变率场。采用二维“高张力样条”（τ＝0.95）内插算法，对非均匀的GPS速度场进行了0.5°×0.5°的均匀预测加密，然后根据每个1°×1°网格边界及内部的9个内插速度矢量，计算出相应1°×1°网格内的均匀应变率。由此得出中亚大三角地震构造域各个区域定量、直观的地壳形变状况。

①整个喜马拉雅弧形地块，在沿印度板块与欧亚板块会聚的方向上承受着强烈的挤压缩短，典型的压缩应变率为（30～60）nanostrain／a，局部量值高达（80～90）nanostrain／a。其中，喜马拉雅地块西部挤压应变比东部更强烈。同时，在大致垂直于板块会聚的方向上，喜马拉雅地块还承受着轻微的横向拉张。

②拉萨地块沿印度板块与欧亚板块会聚的方向上，挤压缩短的典型应变率明显地减小为（20～30）nanostrain／a。但是，在该地块的中部，横向拉伸应变率变得非常明显，典型量值为（20～30）nanostrain／a，局部高达（40～60）nanostrain／a。这种明显拉张应变与该地区广泛发育的南北向正断层相一致。

③羌塘地块的中西部典型的应变是近南北向的挤压和东西向的拉伸，量值在（15～25）nanostrain／a。而在青藏高原的东部，应变场表现出明显的近南北向拉伸，典型量值为（20～30）nanostrain／a。

④祁连地块东边的四川盆地及其以北区域，并无显著的应变，典型量值在10 nanostrain／a以下。整个祁连地块和柴达木地块的应变状况比较均匀，主要表现为NE－SW向挤压和NWSE向拉伸，挤压应变的典型量值为（15～25）nanostrain／a，拉张应变的典型量值为（10～15）nanostrain／a。

⑤青藏高原东南角的川滇地区，应变强烈且方向变化较复杂，不过，大致东西向的扩散拉张和南北向的挤出式压缩特征非常清晰。

⑥天山地区主要以天山南北向挤压缩短变形为主，挤压缩短呈现非均匀分布特征。挤压的最大量值为30 nanostrain／a；

⑦贝加尔地区整体上处于NW－SE向的拉张状态，拉张的最大量值为40 nanostrain／a，在东部地区存在少量NNW向和NNE向的压缩状态，压缩的最大量值为10 nanostrain／a，面膨胀应变的结果表明，贝加尔地区整体上呈膨胀状态。

（5）中亚大三角地震构造域地壳形变块体的划分和平均应变状况的研究。中亚大三角地震构造域内主要有青藏高原－喜玛拉雅挤压区、天山挤压造山带、贝加尔拉张区以及川滇剪切区。本文在前人关于活动地块划分的基础上，以地壳形变特征为依据，划分出5类地壳形变块体，包括准刚性块体、纯拉张变形块体、纯挤压变形块体、拉张为主变形块体、挤压为主变形块体等。在中亚大三角地震构造域所划分出的具体形变块体包括：东塔里木块体、西塔里木块体、阿拉善块体、萨彦块体、阿穆尔块体、准噶尔块体、天山块体、祁连块体、巴颜喀拉块体、羌塘块体、拉萨块体、滇南块体、川滇块体。其中鄂尔多斯块体、塔里木块体、阿拉善块体为刚性块体。

（6）综合采用半无限弹性空间的断裂位错模型、准刚性块体欧拉旋转模型和弹性块体应变与旋转模型，对中亚大三角地震构造域的现今地壳形变GPS速度场进行了拟合、解释。综合与简化中亚大三角地震构造域内主要活动断裂带，建立每一条断裂在半无限空间的三维几何模型，并赋予必要的运动方式先验信息。以中亚大三角地震构造域的1 233个GPS速度矢量为约束，通过在合理的范围内约束和调整各断层段运动速率的先验值，用半无限弹性空间的深断裂位错模型最佳拟合了GPS速度场，并反演获得了所有断裂段的运动速率。同时对于无法利用该模型解释的GPS速度场残差，进一步采用准刚性块体旋转模型以及弹性块体应变与旋转模型进行了拟合、解释。亦即，尝试采用非连续形变模型和连续形变模型相结合的方法，初步解释了中亚大三角地震构造域的地壳运动与形变特征。

（7）分析、研究了中亚大三角地震构造域地震活动性与地壳形变关系，包括地震空间分布与地震矩累积率场、面膨胀率场和最大剪切应变场的高值对应关系，并首次对比了地震矩累积率与地震能量释放量之间的空间对应关系。结果表明：在中亚大三角地震构造域内，虽然应变率场、面膨胀率以及剪切应变率场与研究区域内的地震存在一定的空间相关性，但关系复杂。而采用地震能量释放量与地震矩累积率的空间对比发现，尽管两者有着一定的相关性，但并不明显，最可能的原因是我们的地震能量释放计算，仅采用了100年的地震资料，与板内地震的平均复发周期相比，这样短时间段内统计、计算的地震能量释放分布，完全不能代表一个地震周期内的能量释放分布。

中亚大三角地震构造域；GPS数据处理；准刚性块体；地壳形变；地震活动性

P315；

A；

10.3969／j.issn.0235－4975.2012.08.007

（作者电子信箱，肖根如：xgr541＠21cn.com）