利用Sentinel-1A数据反演2016年台湾地震断层的同震滑动分布（三）
分享
-
2017-01-01
在第二篇中已经通过非线性反演获得均匀滑动断层几何参数（经度、纬度、走向、倾向、滑动角、深度、断层的长度、宽度以及滑动量），本篇也是基于以上结果采用线性反演获得断层精细滑动分布。其中，为了准确反演，同样加入6个轨道参数来线性估计轨道误差；采用固定滑动角，结合非负最小二乘算法求解。
1 数据导入
打开线性反演工具/SARscape/Interferometry/Displacement Modeling/ Linear Inversion，添加降采样形变文件Taiwan__reflat_disp.shp，可以看到总共得到1068个采样点，设置线性函数拟合轨道参数，估计轨道误差。创建工程文件model_linear.xml。
2 线性反演参数设置
添加非线性反演得到的结果文件model_nonlinear.xml。
将非线性反演所得几何参数结果导入线性反演模型中。
设置线性反演参数，首先将反演的断层长度、宽度适当在走向和倾向延伸，使得断层滑动信息可以很好的完全显示；设置反演为固定滑动角，采用非负最小二乘进行解算；设置走向和倾向细分数，对断层进行划分；设置阻尼系数，越大反演精细断层越平滑。
3 结果展示
得到线性反演结果为：
得出的地震矩为
Nm，震级为
6.5级。
将DInSAR形变场和模拟形变场叠加一起，在ArcScene显示。
将反演的断层精细滑动分布和模拟的形变场，在ArcScene显示。
4评定标准
（1）断层几何参数与机构给出结果是否符合。
（2）断层精细滑动分布与实际是否相符。
（3）DInSAR形变结果与模拟形变图是否吻合。
注意事项：
为尽量完全的将断层信息展示，需要调节断层反演的长度、宽度，但是随着长宽的不断变大，相应的模拟形变图的值会变小，建议尽量将断层信息完全展示，但也不应断层过大。
文章来源：http://blog.sina.com.cn/s/blog_764b1e9d0102xl7h.html
1 数据导入
打开线性反演工具/SARscape/Interferometry/Displacement Modeling/ Linear Inversion，添加降采样形变文件Taiwan__reflat_disp.shp，可以看到总共得到1068个采样点，设置线性函数拟合轨道参数，估计轨道误差。创建工程文件model_linear.xml。
2 线性反演参数设置
添加非线性反演得到的结果文件model_nonlinear.xml。
将非线性反演所得几何参数结果导入线性反演模型中。
设置线性反演参数，首先将反演的断层长度、宽度适当在走向和倾向延伸，使得断层滑动信息可以很好的完全显示；设置反演为固定滑动角，采用非负最小二乘进行解算；设置走向和倾向细分数，对断层进行划分；设置阻尼系数，越大反演精细断层越平滑。
3 结果展示
得到线性反演结果为：
得出的地震矩为
将DInSAR形变场和模拟形变场叠加一起，在ArcScene显示。
将反演的断层精细滑动分布和模拟的形变场，在ArcScene显示。
4评定标准
（1）断层几何参数与机构给出结果是否符合。
（2）断层精细滑动分布与实际是否相符。
（3）DInSAR形变结果与模拟形变图是否吻合。
注意事项：
为尽量完全的将断层信息展示，需要调节断层反演的长度、宽度，但是随着长宽的不断变大，相应的模拟形变图的值会变小，建议尽量将断层信息完全展示，但也不应断层过大。
文章来源：http://blog.sina.com.cn/s/blog_764b1e9d0102xl7h.html
