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

REGULARIZED SYNCHRONOUS JOINT INVERSION OF MT AND SEISMIC DATA

作　　者：陈晓[1] 于鹏[1] 张罗磊[1] 李洋[1] 王家林[1] 

机构地区：[1]同济大学海洋地质国家重点实验室,上海200092

出　　处：《地震地质》2010年第3期402-408,共7页Seismology and Geology

基　　金：国家高技术研究发展计划项目(2008AA093001);国家自然科学基金(40674063);国家重点基础研究发展计划项目(2007CB41170602;2007CB411702);同济大学海洋地质国家重点实验室重点项目(MG20080103);国家科技重大专项(2008ZX05005005010HZ)共同资助

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

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

分 类 号：P318.62[天文地球—固体地球物理学]