地震同震滑动分布反演平滑因子的确定

doi:10.11947/j.AGCS.2018.20170724

测绘学报 ›› 2018, Vol. 47 ›› Issue (12): 1571-1580.doi: 10.11947/j.AGCS.2018.20170724

地震同震滑动分布反演平滑因子的确定

王乐洋^1,2,3, 赵雄^1,2

1. 东华理工大学测绘工程学院, 江西南昌 330013;
2. 流域生态与地理环境监测国家测绘地理信息局重点实验室, 江西南昌 330013;
3. 江西省数字国土重点实验室, 江西南昌 330013

收稿日期:2017-12-20 修回日期:2018-08-26 出版日期:2018-12-20 发布日期:2018-12-24
作者简介:王乐洋(1983-),男,博士,副教授,研究方向为大地测量反演及大地测量数据处理。E-mail:wleyang@163.com
基金资助:
国家自然科学基金（41664001；41204003；41874001）；江西省杰出青年人才资助计划（20162BCB23050）；国家重点研发计划（2016YFB0501405）

Determination of Smoothing Factor for the Co-seismic Slip Distribution Inversion

WANG Leyang^1,2,3, ZHAO Xiong^1,2

1. Faculty of Geomatics, East China University of Technology, Nanchang 330013, China;
2. Key Laboratory of Watershed Ecology and Geographical Environment Monitoring, NASG, Nanchang 330013, China;
3. Key Laboratory for Digital Land and Resources of Jiangxi Province, Nanchang 330013, China

Received:2017-12-20 Revised:2018-08-26 Online:2018-12-20 Published:2018-12-24
Supported by:
The National Natural Science Foundation of China (Nos. 41664001;41204003;41874001);The Support Program for Outstanding Youth Talents in Jiangxi Province (No. 20162BCB23050);The National Key Research and Development Program(No. 2016YFB0501405)

摘要/Abstract

摘要： 针对大地测量地震同震滑动分布反演中平滑因子（又称正则化参数）的确定问题，一般采用模型粗糙度与数据拟合残差之间的折中曲线来确定（为便于区分，本文将该方法称为"L曲线"）。本文在L曲线的基础上提出一种确定平滑因子的新方法——折中相交曲线法。模拟试验反演结果表明，利用折中相交曲线法确定的平滑因子反演地震滑动分布各参数精度要优于L曲线法。利用折中相交曲线法确定平滑因子反演拉奎拉与台湾美浓实际地震，并与L曲线法反演结果进行对比分析。分析结果表明，利用折中相交曲线法确定的平滑因子反演拉奎拉与台湾美浓实际同震滑动分布各参数结果均在国内外其他学者研究的范围内，并且利用折中相交曲线法确定平滑因子较L曲线法具有计算效率高、无需依赖数据拟合度、确定平滑因子大小更为合适等优点。

关键词: 地震同震滑动分布, 平滑因子, 折中相交曲线法, L曲线法, 拉奎拉地震, 美浓地震

Abstract: For the determination of the smoothing factor (also known as the regularization parameter) in the co-seismic slip distribution inversion, the compromise curve between the model roughness and the data fitting residual is generally used to determine (in order to distinguish the method proposed in this paper, the method is called "L curve" according to its shape). Based on the L-curve, the eclectic intersection curve as a new method are proposed to determine the smoothing factor in this paper. The results of the simulated experiment show that the inversion accuracy of the parameters of the seismic slip distribution with the smoothing factor determined by the eclectic intersection curve method is better than that of the L curve method. Moreover, the eclectic intersection curve method and the L curve method are used to determine the smoothing factor of L'Aquila and Taiwan Meinong earthquake slip distribution inversion respectively, and the inversion results are compared and analyzed. The analysis results show that the L'Aquila and Meinong of Taiwan actual earthquake slip distribution results are in the range of other scholars at home and abroad, and compared with the L curve method, the eclectic intersection curve method has advantages of high computation efficiency, no need to depend on data fitting degree and more appropriate of smoothing factor and so on.

Key words: co-seismic slip distribution, smoothing factor, eclectic intersection curve method, L curve method, the L'Aquila earthquake, the Meinong earthquake

中图分类号:

P228

王乐洋, 赵雄. 地震同震滑动分布反演平滑因子的确定[J]. 测绘学报, 2018, 47(12): 1571-1580.

WANG Leyang, ZHAO Xiong. Determination of Smoothing Factor for the Co-seismic Slip Distribution Inversion[J]. Acta Geodaetica et Cartographica Sinica, 2018, 47(12): 1571-1580.

参考文献

[1] 许才军, 邓长勇, 周力璇. 利用方差分量估计的地震同震滑动分布反演[J]. 武汉大学学报(信息科学版), 2016, 41(1):37-44. XU Caijun, DENG Changyong, ZHOU Lixuan. Coseismic Slip Distribution Inversion Method Based on the Variance Component Estimation[J]. Geomatics and Information Science of Wuhan University, 2016, 41(1):37-44.
[2] FUNNING G J, PARSONS B, WRIGHT T J, et al. Surface Displacements and Source Parameters of the 2003 Bam (Iran) Earthquake from Envisat Advanced Synthetic Aperture Radar Imagery[J]. Journal of Geophysical Research:Solid Earth, 2005, 110(B9):B09406.
[3] JIANG Zaisen, WANG Min, WANG Yanzhao, et al. GPS Constrained Coseismic Source and Slip Distribution of the 2013M_w6.6 Lushan, China, Earthquake and Its Tectonic Implications[J]. Geophysical Research Letters, 2014, 41(2):407-413.
[4] 许才军, 刘洋, 温扬茂. 利用GPS资料反演汶川M_w7.9级地震滑动分布[J]. 测绘学报, 2009, 38(3):195-201, 215. XU Caijun, LIU Yang, WEN Yangmao. M_w7.9 Wenchuan Earthquake Slip Distribution Inversion from GPS Measurements[J]. Acta Geodaetica et Cartographica Sinica, 2009, 38(3):195-201, 215.
[5] 张国宏, 屈春燕, 宋小刚, 等. 基于InSAR同震形变场反演汶川M_w7.9地震断层滑动分布[J]. 地球物理学报, 2010, 53(2):269-279. ZHANG Guohong, QU Chunyan, SONG Xiaogang, et al. Slip Distribution and Source Parameters Inverted from Co-seismic Deformation Derived by InSAR Technology of Wenchuan M_w7.9 Earthquake[J]. Chinese Journal of Geophysics, 2010, 53(2):269-279.
[6] 温扬茂, 何平, 许才军, 等. 联合Envisat和ALOS卫星影像确定L'Aquila地震震源机制[J]. 地球物理学报, 2012, 55(1):53-65. WEN Yangmao, HE Ping, XU Caijun, et al. Source Parameters of the 2009 L'Aquila Earthquake, Italy from Envisat and ALOS Satellite SAR Images[J]. Chinese Journal of Geophysics, 2012, 55(1):53-65.
[7] 李志才, 张鹏, 金双根, 等. 基于GPS观测数据的汶川地震断层形变反演分析[J]. 测绘学报, 2009, 38(2):108-113, 119. LI Zhicai, ZHANG Peng, JIN Shuanggen, et al. Wenchuan Earthquake Deformation Fault Inversion and Analysis Based on GPS Observations[J]. Acta Geodaetica et Cartographica Sinica, 2009, 38(2):108-113, 119.
[8] WANG Leyang, ZHAO Xiong, GAO Hua. A Method for Determining the Regularization Parameter and the Relative Weight Ratio of the Seismic Slip Distribution with Multi-source Data[J]. Journal of Geodynamics, 2018(118):1-10.
[9] OKADA Y. Surface Deformation due to Shear and Tensile Faults in a Half-space[J]. Bulletin of the Seismological Society of America, 1985, 75(4):1135-1154.
[10] OKADA Y. Internal Deformation due to Shear and Tensile Faults in a Half-space[J]. Bulletin of the Seismological Society of America, 1992, 92(2):1018-1040.
[11] SUN Wenke, OKUBO S, VANÍČEK P. Global Displacements Caused by Point Dislocations in a Realistic Earth Model[J]. Journal of Geophysical Research:Solid Earth, 1996, 101(B4):8561-8577.
[12] 孙建宝, 徐锡伟, 沈正康, 等. 基于线弹性位错模型及干涉雷达同震形变场反演1997年玛尼M_w7.5级地震参数-Ⅰ. 均匀滑动反演[J]. 地球物理学报, 2007, 50(4):1097-1110. SUN Jianbao, XU Xiwei, SHEN Zhengkang, et al. Parameter Inversion of the 1997 Mani Earthquake from InSAR Co-seismic Deformation Field Based on Linear Elastic Dislocation Model-I. Uniform Slip Inversion[J]. Chinese Journal of Geophysics, 2007, 50(4):1097-1110.
[13] 刘洋, 许才军, 温扬茂, 等. 2008年大柴旦M_w6.3级地震的InSAR同震形变观测及断层参数反演[J]. 测绘学报, 2015, 44(11):1202-1209. DOI:10.11947/j.AGCS.2015.20140628. LIU Yang, XU Caijun, WEN Yangmao, et al. The InSAR Coseismic Deformation Observation and Fault Parameter Inversion of the 2008 Dachaidan M_w6.3 Earthquake[J]. Acta Geodaetica et Cartographica Sinica, 2015, 44(11):1202-1209. DOI:10.11947/j.AGCS.2015.20140628.
[14] TIKHONOV A N. Regularization of Incorrectly Posed Problems[J]. Soviet Mathematics Doklady, 1963, 4(1):1624-1627.
[15] GOLUB G H, HEATH M, WAHBA G. Generalized Cross-Validation as a Method for Choosing a Good Ridge Parameter[J]. Technometrics, 1979, 21(2):215-223.
[16] 王振杰. 大地测量中不适定问题的正则化解法研究[D]. 武汉:中国科学院测量与地球物理研究所, 2003. WANG Zhenjie. Research on the Regularization Solutions of Ill-posed Problems in Geodesy[D]. Wuhan:Institute of measurement and Geophysics of the Chinese Academy of Sciences, 2003.
[17] 王振杰, 欧吉坤. 用L-曲线法确定岭估计中的岭参数[J]. 武汉大学学报(信息科学版), 2004, 29(3):235-238. WANG Zhenjie, OU Jikun. Determining the Ridge Parameter in a Ridge Estimation Using L-curve Method[J]. Geomatics and Information Science of Wuhan University, 2004, 29(3):235-238.
[18] HANSEN P C. Analysis of Discrete Ill-posed Problems by Means of the L-curve[J]. SIAM Review, 1992, 34(4):561-580.
[19] HANSEN P C, O'LEARY D P. The Use of the L-curve in the Regularization of Discrete Ill-posed Problems[J]. SIAM Journal on Scientific Computing, 1993, 14(6):1487-1503.
[20] FAN Qingbiao, XU Caijun, YI Lei, et al. Implication of Adaptive Smoothness Constraint and Helmert Variance Component Estimation in Seismic Slip Inversion[J]. Journal of Geodesy, 2017, 91(10):1163-1177.
[21] 鲁洋为, 王振杰. 用U曲线法确定岭估计中的岭参数[J]. 导航定位学报, 2015(3):132-134, 138. LU Yangwei, WANG Zhenjie. Determining the Ridge Parameter in Ridge Estimation Using U-curve Method[J]. Journal of Navigation and Positioning, 2015(3):132-134, 138.
[22] 袁强强, 沈焕锋, 李平湘, 等. 自适应正则化多幅影像超分辨率重建[J]. 中国图象图形学报, 2010, 15(12):1720-1727. YUAN Qiangqiang, SHEN Huanfeng, LI Pingxiang, et al. Adaptively Regularized Muti-frame Image Super-resolution Reconstruction[J]. Journal of Image and Graphics, 2010, 15(12):1720-1727.
[23] XU Peiliang. Truncated SVD Methods for Discrete Linear Ill-posed Problems[J]. Geophysical Journal International, 1998, 135(2):505-514.
[24] JONSSON, ZEBKER H, SEGALL P, et al. Fault Slip Distribution of the 1999M_w 7.1 Hector Mine, California, Earthquake, Estimated from Satellite Radar and GPS Measurements[J]. Bulletin of the Seismological Society of America, 2002, 92(4):1377-1389.
[25] 王乐洋, 许才军, 张朝玉. 一种确定联合反演中相对权比的两步法[J]. 测绘学报, 2012, 41(1):19-24. WANG Leyang, XU Caijun, ZHANG Chaoyu. A Two-step Method to Determine Relative Weight Ratio Factors in Joint Inversion[J]. Acta Geodaetica et Cartographica Sinica, 2012, 41(1):19-24.
[26] WANG Leyang, XU Guangyu. Variance Component Estimation for Partial Errors-in-variables Models[J]. Studia Geophysica et Geodaetica, 2016, 60(1):35-55.
[27] WALTERS R J, ELLIOTT J R, D'AGOSTINO N, et al. The 2009 L'Aquila Earthquake (Central Italy):A Source Mechanism and Implications for Seismic Hazard[J]. Geophysical Research Letters, 2009, 36(17):L17312.
[28] 王乐洋, 李海燕, 温扬茂, 等. 地震同震滑动分布反演的总体最小二乘方法[J]. 测绘学报, 2017, 46(3):307-315. DOI:10.11947/j.AGCS.2017.20160212. WANG Leyang, LI Haiyan, WEN Yangmao, et al. Total Least Squares Method Inversion for Coseismic Slip Distribution[J]. Acta Geodaetica et Cartographica Sinica, 2017, 46(3):307-315. DOI:10.11947/j.AGCS.2017.20160212.
[29] 李海燕. 震源位错模型参数反演方法研究[D]. 南昌:东华理工大学, 2016. LI Haiyan. Research on Inversion Method for Hypocenter Parameter Base on Dislocation Model[D]. Nanchang:East China University of Technology, 2016.
[30] ANZIDEI M, BOSCHI E, CANNELLI V, et al. Coseismic Deformation of the Destructive April 6, 2009 L'Aquila Earthquake (Central Italy) from GPS Data[J]. Geophysical Research Letters, 2009, 36(17):L17307.
[31] CHELONI D, D'AGOSTINO N, D'ANASTASIO E, et al. Coseismic and Initial Post-seismic Slip of the 2009M_w6.3 L'Aquila Earthquake, Italy, from GPS Measurements[J]. Geophysical Journal International, 2010, 181(3):1539-1546.
[32] 王乐洋, 高华, 冯光财. 2016年台湾美浓M_w 6.4地震震源参数的InSAR和GPS反演[J]. 地球物理学报, 2017, 60(7):2578-2588. WANG Leyang, GAO Hua, FENG Guangcai. InSAR and GPS Inversion for Source Parameters of the 2016M_w6.4 Meinong, Taiwan Earthquake[J]. Chinese Journal of Geophysics, 2017, 60(7):2578-2588.
[33] FENG Guangcai, LI Zhiwei, SHAN Xinjian, et al. Geodetic Model of the 2015 April 25M_w 7.8 Gorkha Nepal Earthquake and M_w7.3 Aftershock Estimated from InSAR and GPS Data[J]. Geophysical Journal International, 2015, 203(2):896-900.
[34] HUANG M H, TUNG H, FIELDING E J, et al. Multiple Fault Slip Triggered Above the 2016 M_w 6.4 Meinong Earthquake in Taiwan[J]. Geophysical Research Letters, 2016, 43(14):7459-7467.
[35] LEE S J, YEH T Y, LIN Y Y. Anomalously Large Ground Motion in the 2016 ML 6.6 Meinong, Taiwan, Earthquake:A Synergy Effect of Source Rupture and Site Amplification[J]. Seismological Research Letters, 2016, 87(6):1319-1326.

地震同震滑动分布反演平滑因子的确定

Determination of Smoothing Factor for the Co-seismic Slip Distribution Inversion

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 4

编辑推荐

Metrics

本文评价

[1]	王乐洋, 陈涛, 邹传义. 病态乘性误差模型的加权最小二乘正则化迭代解法及精度评定[J]. 测绘学报, 2021, 50(5): 589-599.
[2]	王乐洋, 熊露雲. 一种Partial EIV半参数模型的系统误差处理方法[J]. 测绘学报, 2018, 47(1): 25-34.
[3]	王乐洋, 李海燕, 温扬茂, 许才军. 地震同震滑动分布反演的总体最小二乘方法[J]. 测绘学报, 2017, 46(3): 307-315.
[4]	刘晓刚李迎春肖云翟振和. 重力与磁力测量数据向下延拓中最优正则化参数确定方法研究[J]. 测绘学报, 2014, 43(9): 881-887.