摘要
由于此地区构造复杂,USGS,Harvard,台湾中央气象局等机构给出的震源机制解不尽相同,断层参数不确定性大,因此将断层的几何参数当做未知参数,以GPS同震位移资料为约束进行非线性反演。首先,假设断层同震破裂面为矩形,并且破裂面上的滑动是均匀的,在均匀弹性半空间的假定下,可以用Okada位错公式计算由断层参数与滑动量确定的地表位移场,将其与GPS观测位移场的加权残差和作为非线性反演的目标函数。然后,固定断层的几何参数,将断层沿走向和倾向划分为若干个子断层,用Okada位错公式计算格林函数矩阵,并且用二次光滑的拉普拉斯算了作为约束,线性反演滑动分布。