地震地质 ›› 2010, Vol. 32 ›› Issue (3): 402-408.DOI: 10.3969/j.issn.0253-4967.2010.03.006

• 大地电磁法研究与应用 • 上一篇    下一篇

大地电磁与地震正则化同步联合反演

陈晓, 于鹏, 张罗磊, 李洋, 王家林   

  1. 同济大学海洋地质国家重点实验室, 上海, 200092
  • 收稿日期:2010-01-17 修回日期:2010-07-12 出版日期:2010-09-30 发布日期:2011-01-14
  • 作者简介:陈晓,男,1986年生,2004年毕业于江汉石油学院(现长江大学),现于同济大学攻读固体地球物理学博士学位.主要研究方向为大地电磁测深法以及多种地球物理方法的联合反演,电话:021-65982207,E-mail:dwjhtj@hotmail.com.
  • 基金资助:
    国家高技术研究发展计划项目(2008AA093001);国家自然科学基金(40674063);国家重点基础研究发展计划项目(2007CB41170602,2007CB411702);同济大学海洋地质国家重点实验室重点项目(MG20080103);国家科技重大专项(2008ZX05005005010HZ)共同资助

REGULARIZED SYNCHRONOUS JOINT INVERSION OF MT AND SEISMIC DATA

CHEN Xiao, YU Peng, ZHANG Luo-lei, LI Yang, WANG Jia-lin   

  1. State Key Laboratory of Marine Geology, Tongji University, Shanghai 200092, China
  • Received:2010-01-17 Revised:2010-07-12 Online:2010-09-30 Published:2011-01-14

摘要: 文中在于鹏等提出的电阻率和速度随机分布的大地电磁与地震联合反演方法的基础上,将Tikhonov正则化思想引入到联合反演中,加入先验信息进行模型约束,以最小模型为稳定器,采用L曲线方法来确定近似最佳的正则化因子。考虑到线性寻优算法容易陷入局部极小,文中采用非线性的模拟退火方法来实现大地电磁与地震的同步联合反演。通过模型试验的对比分析,我们认为加入有效模型约束的正则化联合反演可以比单纯考虑数据拟合的联合反演和单独反演方法更有效地提高解的稳定性和计算效率,获得更接近实际而且稳定的解。

关键词: 正则化, L曲线, 联合反演, 模拟退火, 大地电磁, 地震

Abstract: All the geophysical inversion problems are ill-posed.The application of regularization theory and joint inversion provide us ideas to solve the ill-posed inverse geophysical problems.Regularization theory has been widely used in single geophysical inversion problem which has achieved certain effects;however,what geophysicists concerned in the previous studies on joint inversion just was data misfit with little consideration about regularization theory.In this paper,based on the joint inversion method with random resistivity and velocity distributions proposed by YU Peng,the conception of regularization is introduced into synchronous joint inversion of magnetotelluric(MT) and seismic data.We use priori information to improve the stability of solution,and set stabilizer as the smallest model.How to meet the best trade-off between data misfit and model constraint,and what's intellectual choice for the regularization factor have always been the hotspots in the geophysical inversion problem.We use L-curve method to determine the quasi-optimal value of the regularization parameter,and apply it into joint inversion.To avoid immerging in local minimum of linear inversion method,nonlinear simulated annealing method is used to realize joint inversion.The forward modeling of MT data is obtained by using finite element method and the seismic travel-time computing is obtained by using improved ray-tracing method.Testing the regularized and conventional synchronous joint inversion on synthetic data sets shows that regularized joint inversion with effective model constraint has powerful availability to improve stability of solution and efficiency of computation,and to help us to obtain a stable solution which is close to the fact.

Key words: regularization, L-curve, joint inversion, simulated annealing, MT, earthquake

中图分类号: