基于粘弹性球体地球模型的震后位移与重力变化计算软件

在考虑地球曲率、成层结构、可压缩性与自重的前提下,Tanaka等[1-2]提出一套较为完备的粘弹性球体位错理论,可计算全球任意位置由地震产生的同震与震后形变(含位移、重力变化、大地水准面变化)。Gao等[3]给出了与上述理论相匹配、界面友好的计算软件,能计算30个震后时间点对应的震后形变。本研究针对Gao等的软件进行改进,可计算震后任意时间点对应的震后形变。新软件由3个部分组成:1)与32个震后时间点相对应的32套离散格林函数数值框架;2)格林函数插值计算程序,可针对上述32套格林函数数值框架进行插值运算,输出任意震后时间点对应的格林函数数值结果;3)积分计算程序,调用上述格林函数数值结果,计算任意类型地震在地表任意位置产生的同震与震后形变。一般情况下,使用者只需按要求准备辅助文件,提供发震断层模型和观测站位置信息,以及震中周围地区地幔粘滞性因子,先后运行格林函数插值计算程序和积分计算程序,即可计算出目标地震在地表任意位置产生的同震与震后形变。本文基于粘弹球体位错理论与弹性球体位错理论,分别计算2011年日本MW9.0地震引起的远场同震位移,2套结果的高度一致性证明了新程序的正确性。最后,介绍需要注意的若干事项,便于使用者掌握该软件。