摘要：许多研究人员利用GPS测量的速度资料计算了地应变率场,但其结果差异较大.本文将地质统计学中的Kriging方法引入到GPS观测的速度场研究中,通过Kriging插值得到青藏高原及邻区均匀网格节点上的速度值,然后运用有限单元中形函数(Lagrange插值函数)的求导方法,计算每个网格单元积分点处的地应变率分量,从而获得青藏高原及邻区的地应变率场的分布.计算结果显示,青藏高原主体处在南北向受挤压、东西向被拉张的应变状态之中,但高原东部地区则正好相反,即南北向拉张、东西向出现挤压.青藏高原及邻区主应变率的方位与震源机制解中P轴、T轴的方向基本一致;最大主压应变率的高值区分布在喜马拉雅主边界冲断带及附近地区,高原内部出现主张应变率大于压应变率的现象,且高原内部处在拉张应变状态.面膨胀率结果也表明,喜马拉雅山及附近地区为面收缩区,而高原内部其他地区主要为膨胀区;最大剪应变率分布清晰地显示出青藏高原周边的主要断裂带轮廓.文中的应变率计算结果预示青藏高原及周边地区现今的地应变与较长期的地质活动之间有一定的继承关系.

关键词：kriging方法 gps观测 地应变率场 青藏高原 

单位：中国地震局地壳应力研究所; 北京; 100085; 北京大学地空学院地球物理学系; 北京; 100871; 北京大学地空学院地球物理学系; 北京; 100871; 中国科学院研究生院计算地球动力学实验室; 北京; 100039