周 宇 聂兆生 贾治革 陈正松 杨少敏

(中国地震局地震研究所(地震大地测量实验室)，武汉 430071)

基于非均匀地壳模型的格林函数法反演2001年昆仑山口西Ms8.1地震的同震滑移*

周 宇 聂兆生 贾治革 陈正松 杨少敏

(中国地震局地震研究所(地震大地测量实验室)，武汉 430071)

建立基于地壳不均匀的三维有限元模型，将昆仑山口西Ms8.1地震产生的长约400 km的破裂带在走向上分为16段，生成各个破裂段滑移相对于地表位移的格林函数，并以GPS观测的水平位移为约束，用阻尼最小二乘法反演破裂滑移分布，结果表明反演结果与地质考察的地表破裂数据接近。用反演结果重构的地表位移场与观测值一致，并体现了断裂带南北两侧不同的形变模式。

有限元;最小二乘法;昆仑山口西Ms8.1地震;同震滑移;GPS

1 引言

昆仑山口西Ms8.1地震发生在东昆仑断裂带上，该断裂带是青藏高原一条古老的缝合带，其左旋走滑非常活跃。震后，美国哈佛大学、美国地质调查局(USGS)、日本东京大学地震研究所和中国地震局地球物理研究所先后给出了此次地震的震中位置及矩张量震源机制解，结果略有差异，但多认为此次主震的破裂方式以走滑破裂为主，兼有少量的逆冲滑动分量［1，2］。野外地质考察发现，地震沿东昆仑断裂带西段形成了长约350～420 km的地表破裂带，位错以走滑为主，最大左旋水平错距 6～7 m［3-5］。

文献［1，2］运用全球地震波资料反演了该地震的时空破裂过程。文献［1］的反演结果为:断层的破裂长度约380 km、宽约30 km、最大滑动量5.8 m、平均错动(左旋走滑)1.8 m;文献［2］假定断层破裂深度为40 km，得到的断层最大滑动量为2.2 m，平均滑动量仅1.2 m。……