摘要
本文尝试利用GRACE观测资料,结合白重力、粘弹性、球体分层位错理论反演该地区粘滞性构造。但是,在应用时发现该理论模型在计算时存在一些问题,即对格林函数积分时未将距离计算点过近的了断层重新分块,进而导致最终的计算结果错误。为解决这个问题,我们利用弹性球体位错理论积分程序对该粘弹性位错理论的格林函数进行积分。根据该x弹位错理论与球体位错理论格林函数定义时的差异及在计算格林函数程序时的数据调整发现弹性球体位错的重力格林函数是x弹球体位错重力格林函数的2.6倍,垂直位移格林函数4.059×102倍,大地水准面格林函数4.059×10-1倍。此外,勃弹球体位错理论的水平走滑断层和垂直倾滑断层的格林函数与弹性球体位错理论相应断层的格林函数符号相反。调整后模拟了2004年苏门答腊Mw9.3地震引起的同震重力变化,与弹性球体位错理论模型模拟的结果具有很好的一致性;并将结果与GRACE重力变化点位时间序列进行比较,结果显示该地区断层两侧粘滞性系数存在差异,下盘所在区域粘滞性系数在1018 Pas量级,上盘所在区域粘滞性系数在lO19Pas量级。