利用Sentinel-1A数据反演2016年台湾地震断层的同震滑动分布（二）
分享
-
2017-01-01
基于降采样得到的InSAR视线向形变场结果，采用非线性反演均匀滑动断层几何参数（经度、纬度、走向、倾向、滑动角、深度、断层的长度、宽度以及滑动量），为了准确反演，加入6个轨道参数来线性估计轨道误差。其中，采用Levemberg-Marquardt最小二乘优化算法迭代，进行9个几何参数和6个轨道参数求解。
1 数据导入
打开非线性反演工具/SARscape/Interferometry/Displacement Modeling/Non-Linear Inversion，添加降采样形变文件Taiwan__reflat_disp.shp，可以看到总共得到1068个采样点，设置线性函数拟合轨道参数，估计轨道误差。创建工程文件model_nonlinear.xml。
2 CMT几何参数区间查询
输入CMT网址：www.globalcmt.org/CMTsearch.html，设置好查询开始时间和查询天数。
下图为得到查询结果，复制查询得到的地震编号。
3 模型区间设置
点击new进入Okada模型参数设置界面。
根据在CMT中查询得到的地震编号，查询CMT给出的建议区间参数。
设置好区间参数，点击Start，软件将根据预设的区间进行最优解搜索。
搜索结束，获得搜索结果，此时，对比得到结果与设置区间参数，如果结果在区间的上限或者下限，需要适当扩大搜索范围，不断调节参数，直至得到合适搜索区间。
经过不断调整参数区间，最终确定区间为：
得到非线性反演结果为：
得出的地震矩为
Nm，震级为
6.4级
注意事项：
区间设置是比较麻烦的过程，需要不断调节，首先建议设置大的搜索范围，待最优解趋于稳定，再将区间适当调小，以免造成多解情况。
致谢：
感谢中南大学 高华硕士在Okada模型算法中的交流指导。
文章来源：http://blog.sina.com.cn/s/blog_764b1e9d0102xl55.html
1 数据导入
打开非线性反演工具/SARscape/Interferometry/Displacement Modeling/Non-Linear Inversion，添加降采样形变文件Taiwan__reflat_disp.shp，可以看到总共得到1068个采样点，设置线性函数拟合轨道参数，估计轨道误差。创建工程文件model_nonlinear.xml。
2 CMT几何参数区间查询
输入CMT网址：www.globalcmt.org/CMTsearch.html，设置好查询开始时间和查询天数。
下图为得到查询结果，复制查询得到的地震编号。
3 模型区间设置
点击new进入Okada模型参数设置界面。
根据在CMT中查询得到的地震编号，查询CMT给出的建议区间参数。
设置好区间参数，点击Start，软件将根据预设的区间进行最优解搜索。
搜索结束，获得搜索结果，此时，对比得到结果与设置区间参数，如果结果在区间的上限或者下限，需要适当扩大搜索范围，不断调节参数，直至得到合适搜索区间。
经过不断调整参数区间，最终确定区间为：
得到非线性反演结果为：
得出的地震矩为
注意事项：
区间设置是比较麻烦的过程，需要不断调节，首先建议设置大的搜索范围，待最优解趋于稳定，再将区间适当调小，以免造成多解情况。
致谢：
感谢中南大学 高华硕士在Okada模型算法中的交流指导。
文章来源：http://blog.sina.com.cn/s/blog_764b1e9d0102xl55.html
