GPS观测研究汶川8.0级强震前地壳内部应力场变化特征.pdfVIP

《测绘通报》测绘科学前沿技术论坛论文集 1 GPS观测研究汶川8.0级强震前 地壳内部应力场变化特征∗ 张永志 吴小利 王卫东 刘宁 长安大学地质工程与测绘学院，中国陕西，西安，710054 摘要 汶川大地震造成了许多的人员伤亡和财产损失，造成强地震的主要原因是龙门山地区地壳中局部应力场的释放和地壳内本身 存在断层等地质缺陷的滑动。本文通过格林第二内部积分公式导出了地面位移与地壳内应力变化的关系，分析和讨论了利用地面位移变化 速率来计算地球内应力场变化的可能性和计算方法; 最后, 采用2001-2005年间在中国四川地区观测的GPS位移场速率数据, 对汶川地震区 震前地壳内深度为15公里处的最大主应力、最大剪应力变化进行了计算分析。结果表明，四川地区地壳内的局部最大主压应力变化在震中 附近变化较大，而最大剪应力变化不明显。 关键词 GPS速率 反演 应力场变化 汶川地震 一、 引言 GPS(global positioning system ) 观测可精确地获得地球表面的三维位移场变化速率. 而GPS 技术观测到的地球表面位移场变化可能是地球内部应力场变化在地面的直接响应。因此，利用地面观 测到的位移场速率来研究地壳内部的应力场变化与地震活动等地球物理现象的关系无疑具有重要的 科学意义。GPS数据在研究强地震活动引起的地面变形在国外已取得许多成功的例子。如Snay等(1991) 采用GPS资料,研究了与Loma Prieta 地震有关的水平形变[1] .Ching （2005）利用GPS观测的位移场研 究了中国台湾地区的构造逃逸变形[2] 。Peyrat等（2006）采用测震、强地面运动、GPS 和 InSAR 数 据对发生在开罗北面Tarapaca 地区的 Mw7.7级地震的震源参数进行了反演[3] 。Johnson等（2006）采 用GPS数据和震后应力变化数据对断层的滑动摩擦参数进行了反演计算研究[4] 。我国学者在利用GPS 资料研究地壳形变和地质块体间的相对运动中也取得了一些结果.如黄立人(黄立人等, 1998)用GPS 资料研究了我国华北地区和青藏高原的地壳形变和地质块体间的相对运动[5] 、王琪等利用GPS资料研 究我国云南丽江7. 0 级地震前的地壳形变和地质块体间的相对运动(王琪等, 1998)[6] ，黄立人等 (1999) 利用华北地区GPS 资料,采用有限元方法计算了华北地区地壳中的应力场分布[7] .本文在格林 作者简介：张永志，男,1961年1月生，博士，教授（博士生导师），长安大学地质工程与测绘学院工作，主要从事地壳形变、 重力测量、INSAR、GPS测量等的教学与研究。 《测绘通报》测绘科学前沿技术论坛论文集 2 [8] 第内部积分公式 的基础上，导出了地面位移场与地球内应力场变化的近似表示关系，并利用四川地 区（2001-2005）观测的GPS地面位移场速率,对我国汶川强震活动区震前的地壳深部最大主应力场和 剪应力变化分布进行了计算研究，并分析了该地应力场空间变化分布与强震活动位置、断层分布的关 系。得出了一些有意义的结论和认识。 二、 地面位移场与地球内部应力场的关系 为了利用地面观测的位移场计算地球内部的应力场，本文首先建立地面位移场速率与地球内部应 力的关系。如图1所示,n为曲面S （S S1 +S 2 ） 的单位外法线向量,设体积内部任意两点P、Q的位移 [8] 向量α、β在体积V 内与S 面上二阶连续、可微。由格林第二内部积分公式 （管泽霖等，1981）可得 β o s1