[发明专利]基于多源SAR数据的同震震后时空滑动分布联合反演方法

权利要求书

1.一种基于多源数据附加对数约束的时空断层滑动联合反演方法，其特征在于，包括以下步骤：

步骤一：利用InSAR技术处理多源SAR数据得到地震区域的同震和震后雷达视线向观测值，并通过地理编码将InSAR测量值从雷达坐标系下转换到通用的横轴墨卡托投影坐标系；

步骤二：通过四叉树降采样InSAR观测值，利用多峰值粒子群算法求解包括发震断层倾角、发震断层走向、均一滑动量的断层几何参数，然后固定断层倾角和走向并将其细分为若干有限元断层；

步骤三：基于弹性位错理论构建地表形变和断层时空滑动分布之间的格林函数模型，其表示如下

其中X为断层时空滑动向量，D为InSAR技术得到的地表形变，G为通过弹性位错模型得到的格林函数矩阵，B为设计系数矩阵，ε为模型残差；T表示获取的总计震后SAR数据数量，是断层有限元i在震后时刻j时的同震滑动向量，是断层有限元i在震后时刻j时的时空余滑累积向量，其中，i＝1,2,...,P；j＝1,2,...,T，P为断层有限元个数；表示不同时间节点的同震形变场，其中，j＝0,1,2,...,T，表示不同时间节点的震后形变场，其中，i＝1,2,...,T；j＝0,1,2,...,T-1；

而震后余滑遵循对数式衰减的特征，由下式表征：

x_post＝A·log(1+(t-t₀)/τ) (2)

其中，A为描述震后余滑大小的幅度系数，t为震后SAR数据获取的时间节点，t₀为地震发生时的时间节点，τ为描述震后余滑衰减速度的时间常数；

联合公式(1)和(2)将该公式(1)模型中待求的震后时空余滑转化为对数函数模型，可得如下矩阵形式的附加对数约束的同震震后时空滑动分布联合反演模型：

其中，B₁＝[1…1]_1×P和B₀＝[0…0]_1×P；

步骤四：求解联合反演模型，得到所有断层有限元上的同震滑动分布以及每个断层有限元上的余滑幅度系数和衰减时间常数，进而得到研究时间段内整个发震断层上的时空滑动分布。

2.根据权利要求1所述的基于多源数据附加对数约束的时空断层滑动联合反演方法，其特征在于，使用附加约束的非线性最小二乘求解所述联合反演模型。

3.根据权利要求1所述的基于多源数据附加对数约束的时空断层滑动联合反演方法，其特征在于，公式(3)所示联合反演模型参数仅有2P+3，其中包括P个断层有限元上的同震滑动和P个余滑幅度系数，一个均一断层滑动衰减时间常数τ以及嵌套在格林矩阵中的两个断层滑动角。

4.根据权利要求1-3中任一项所述的基于多源数据附加对数约束的时空断层滑动联合反演方法，其特征在于，所述地震形变为任意区域任意地震类型造成的形变。