摘要：地球介质的黏度是地球动力学研究中的重要参数.在考虑震后震源区介质物性变化的情况下,本文提出利用GPS观测的时间序列资料反演地壳、地幔黏度的方法.计算中,将遗传算法与Marquardt法相结合,利用遗传算法搜索结果作为Marquardt法反演的初始模型,运用微扰法形成Jacobian矩阵.根据余震分布及地震波速度变化确定震源区的范围及介质的物性变化情况,利用黏弹性有限单元模型（Maxwell体）及台湾集集地震后GPS观测的时间序列数据,反演了台湾地区地壳、地幔的黏性系数.反演计算及数值实验表明,该方法反演黏度的效率高,稳定性好.反演结果初步显示台湾地区的地壳、地幔的黏度分别为1.2×10^18Pa·s和3.6×10^19Pa·s.

关键词：黏度系数反演 黏弹性有限元 gps观测 集集地震 

单位：中国地震局地壳应力研究所; 北京100085; 北京大学地球物理学系; 北京100871