王 帅 张永志 姜永涛 姚志军 刘国仕

1)长安大学地质工程与测绘学院，西安 710054

2)新疆工程学院，乌鲁木齐830001

汶川地震发生后，国内外学者基于不同观测数据对汶川地震断层的滑动分布特征进行了较为深入的研究［1－4］。利用远场体波波形记录的地震数据，纪晨［5］对滑动分布进行反演，其结果显示，断层的最大滑动量为9 m，地震矩Mo=0.76 ×1021N·m。王卫民［6］利用远场地震波和近场GPS 实测同震位移数据，以野外地质调查结果作为约束联合反演，得到断层滑动量最大值为12.49 m，地震矩为1.04 ×1021N·m 。刁发启［7］利用近场GPS 同震位移数据反演得到的滑动量最大值与纪晨的研究相吻合，但地震矩略大，为0.92 ×1021N·m。沈正康利用GPS和InSAR 数据［8］得到的最大滑动量和地震矩均小于前3 种模型，分别为5 ～6 m、0.73 ×1021N·m。

为研究滑动分布和位错模型对同震位移的影响，解释震区大地测量、地球物理等观测数据所蕴涵的信息，本文以汶川地震为例，利用上述4 种断层滑动模型，分别采用Okada 经典的弹性半空间位错模型［9－10］和孙文科球形位错模型［11－13］，对同震位移场进行模拟计算，并将模拟结果与GPS 实测值进行对比，得出适用于该区域最优的断层滑动分布和位错模型。最后基于统计学中二元方差分析方法，研究了断层滑动模型和位错模型对正演结果的显著性影响，并对其显著性影响的空间特征进行分析。

1 断层滑动分布模型和位错模型的对比分析

选择上述4 种断层滑动模型和两种位错模型进行组合，模拟龙门山断裂带区域同震位移场(图1)，并计算相应的位移场残差(观测值与模拟值之差)、残差均值和方差，见表1。……