吉林大学学报(地球科学版) ›› 2015, Vol. 45 ›› Issue (3): 952-961.doi: 10.13278/j.cnki.jjuese.201503304

• 地球探测与信息技术 • 上一篇    下一篇

大地电磁测深与地震初至波走时交叉梯度反演

李桐林1, 张镕哲1, 朴英哲1,2   

  1. 1. 吉林大学地球探测科学与技术学院, 长春 130026;
    2. 金策工业综合大学资源探测工学系, 朝鲜 平壤 999093
  • 收稿日期:2014-09-13 发布日期:2015-05-26
  • 作者简介:李桐林(1962),男,教授,博士生导师,主要从事电磁法理论和应用研究,E-mail:litl@jlu.edu.cn。
  • 基金资助:

    中国地质调查局项目([2013]01-060-004);国家重大科学仪器设备开发专项项目(2011YQ05006009)

Joint Inversion of Magnetotelluric and First-Arrival Wave Seismic Traveltime with Cross-Gradient Constraints

Li Tonglin1, Zhang Rongzhe1, Pak Yongchol1,2   

  1. 1. College of Geo-Exploration Science and Technology, Jilin University, Changchun 130026, China;
    2. Resource Exploration and Engineering Department, Kimchaek University of Technology, Pyongyang 999093, Korea
  • Received:2014-09-13 Published:2015-05-26

PDF (PC)

801

摘要：

在地球物理勘探过程中,观测数据有限且存在误差,导致反演结果存在非唯一解。为了解决非唯一性问题,笔者研究了交叉梯度数学理论,实现了用交叉梯度约束2D 大地电磁和地震初至波走时的联合反演,并建立了2个模型。通过各自模型实验和不同模型间的对比,证明了交叉梯度联合反演比单独反演修复几何形态和物性参数的效果好,且交叉梯度与岩石物性无关;利用得到的交叉梯度图进一步证明了联合反演模型相似度高,也从交叉绘图中得到联合反演的物性关系更好的结论,反向证明了联合反演的正确性。

关键词: 大地电磁测深, 地震初至波走时, 交叉梯度, 联合反演

Abstract:

In the process of geophysical exploration, observation data are limited and errors exist, leading to the presence of non-unique inversion results. In order to solve the problem of non-uniqueness, the authors study the mathematical theory of cross-gradient, and achieve the joint inversion of 2D magnetotelluric and first-arrival wave seismic traveltime with cross-gradient. At the same time, we establish two models. Through test of each model and comparison between the two models, it is proven not only that the result of the cross-gradient joint inversion is better than the single inversion in remediation on geometry and physical property parameters, but also that the cross-gradient is independent from the petro-physical properties. Further, the cross-gradient figure evidences the high similarity of the joint inversion models, and the cross-plot of joint inversion presents a better correlation with physical properties. These, in turn, prove the validity of the joint inversion.

Key words: magnetotelluric, first-arrival wave seismic traveltime, cross-gradient, joint inversion

中图分类号: 

  • P631.3

[1] 杨辉,戴世坤,宋海斌,等.综合地球物理联合反演综述[J].地球物理学进展,2002,17(2):262-271. Yang Hui, Dai Shikun, Song Haibin, et al. Integrated Geophysical Joint Inversion Review[J]. Progress in Geophysics, 2002, 17(2): 262-271.

[2] 于鹏,王家林,吴健生,等.地球物理联合反演的研究现状和分析[J].勘探地球物理进展,2006,29(2):87-93. Yu Peng, Wang Jialin, Wu Jiansheng, et al. Research Status and Analysis of Joint Inversion[J]. Progress in Exploration Geophysics,2006,29(2):87-93.

[3] Heincke B, Hobbs R. Joint Inversion of MT, Gravity and Seismic Data Applied to Sub-Basalt Imaging[C]//Annual Meeting. New Orleans: SEG, 2006: 784-787.

[4] Colombo D, Stefano M D. Geophysical Modeling via Simultaneous Joint Inversion of Seismic Gravity, and Electromagnetic Data: Application to Prestack Depth Imaging[J]. The Leading Edge, 2007, 26(3): 326-331.

[5] Haber E, Oldenburg D. Joint Inversion: A Structural Approach Inverse Problems[J]. IOP Science, 1997, 13: 63-77.

[6] Zhang J, Morgan F D. Joint Seismic and Electrical Tomography[C]// Symposium on the Application of Geophysics to Engineering and Environmental Problems. [S.l.]: SEG,1997: 391-396.

[7] Gallardo L A, Meju M A. Characterization of Heterogeneous Near-Surface Materials by Joint 2D Inversion of DC Resistivity and Seismic Data[J]. Geophysical Research Letters, 2003, 30(13): 1658-1661.

[8] 彭淼,谭捍东,姜枚,等.基于交叉梯度耦合的大地电磁与地震走时资料三维联合反演[J].地球物理学报,2013,56(8):2728-2738. Peng Miao, Tan Handong, Jiang Mei, et al. Three-Dimension Joint Inversion Magnetotelluric and Seismic Travel Time Data with Cross-Gradient Constraints[J]. Chinese Journal of Geophysics, 2013, 56(8):2728-2738.

[9] 王俊,孟小红,陈朝曦,等.交叉梯度理论及其在地球物理联合反演中的应用[J].地球物理学进展,2013, 28(4):2094-2103. Wang Jun, Meng Xiaohong, Chen Zhaoxi, et al. The Theory of Cross-Gradient and Its Application in Geophysical Joint Inversion[J]. Progress in Geophysics, 2013, 28(4):2094-2103.

[10] Gallardo L A, Meju M A. Joint Two-Dimensional DC Resistivity and Seismic Travel Time Inversion with Cross-Gradients Constraints[J]. Jouranl of Geophsical Research, 2004, 109: 3311.

[11] Vidale J. Finite-Difference Calculation of Travel-Ti-mes in Three Dimenisons[J]. Geophysics, 1990, 55: 521-526.

[12] Zelt C A, Barton P J. Three-Dimensional Seismic Refraction Tomography: A Comparison of Two Methods Applied to Data from the Faeroe Basin[J]. Geophys Res, 1998, 103: 7187-7210.

[13] 刘小军,王家林,于鹏,等.基于二次场的二维大地电磁有限元法数值模拟[J].同济大学学报:自然科学版,2007,35(8):1113-1117. Liu Xiaojun, Wang Jialin, Yu Peng, et al. Numerical Simulation of Two-Dimensional Magnetotelluric Secondary Field Based on the Finite Element Method[J]. Journal of Tongji University:Natural Science, 2007, 35 (8):1113-1117.

[14] Wannamaker P E, Stodt J A, Rijo L. A Stable Finite Element Solution for Two-Dimensional Magnetotelluric Modeling[J]. Geophys J R,1987, 88: 277-296.

[15] de Lugao P P, Wannamaker P E. Calculating the Two-Dimensional Magnetotlluric Jacobian in Finite Elements Using ReciProcity[J]. Geophys J,1996, 127: 806-810.

[16] 朴英哲,李桐林,刘永亮,等.在大地电磁二维Occam反演中求取拉格朗日算子方法的改进[J].吉林大学学报:地球科学版,2014,44(2):660-667. Pak Yongchol, Li Tonglin, Liu Yongliang, et al. Improvement of Choosing Lagrange Multiplier on MT 2D Occam Inversion[J]. Journal of Jilin University : Earth Science Edition, 2014, 44(2): 660-667.

[17] de Groot-Hedlin C D, Constable S C. Occam's Inversion to Generate Smooth, Two-Dimensional Models from MagnetotelluricData[J]. Geophysics, 1990, 93: 1613-1624.

[18] Constable S C, Parker K L, Constable C G. Occam's Inversion a Practical Algorithm for Generating Smooth Models from EM Sounding Data[J]. Geophysics,1987, 52: 289-300.
[1] 高科宁, 邓居智, 刘文玉, 余辉, 陈辉, 周聪, 王彦国, 刘宝山.

黑龙江省三合屯金矿床深部电性结构与成矿模式 [J]. 吉林大学学报(地球科学版), 2022, 52(3): 684-700.
[2] 翁望飞, 王德恩, 胡召齐, 朱强, 汪启年, 吴冀明, 唐国强, 张鹏飞. 安徽省休宁—歙县整装勘查区晚侏罗世逆冲推覆构造特征及其与成矿作用关系[J]. 吉林大学学报(地球科学版), 2020, 50(5): 1518-1538.
[3] 张鹏辉, 张小博, 方慧, 王小江, 刘建勋. 地球物理资料揭示的嫩江—八里罕断裂中段深浅构造特征[J]. 吉林大学学报(地球科学版), 2020, 50(1): 261-272.
[4] 薛林福, 祝铭, 李文庆, 刘文玉, 刘正宏, 刘泽宇. 岩浆泡破裂引发地震的模式——以吉林松原2013年地震群为例[J]. 吉林大学学报(地球科学版), 2018, 48(6): 1865-1875.
[5] 韩松, 刘国兴, 韩江涛. 华南地区进贤-柘荣剖面深部电性结构[J]. 吉林大学学报(地球科学版), 2016, 46(6): 1837-1846.
[6] 张广智，杜炳毅，陈怀震，高建虎，李超，李远. 纵横波弹性阻抗联合反演方法[J]. 吉林大学学报(地球科学版), 2014, 44(5): 1695-1704.
[7] 朴英哲，李桐林，刘永亮. 在大地电磁二维Occam反演中求取拉格朗日乘子方法改进[J]. 吉林大学学报(地球科学版), 2014, 44(2): 660-667.
Viewed
Full text


Abstract

Cited
  Shared   
  Discussed   
[1] 尤敏鑫，刘建民. 同位素地球化学在峨眉山大火成岩省研究中的应用现状与进展[J]. 吉林大学学报(地球科学版), 2014, 44(4): 1231 -1243 .
[2] 杨春梅, 李洪奇,陆大卫,张方礼,高 原,邵英超. 不同驱替方式下岩石电阻率与饱和度的关系[J]. J4, 2005, 35(05): 667 -671 .
[3] 祝洪臣，张炯飞，权 恒. 大兴安岭中生代两期成岩成矿作用的元素、同位素特征及其形成环境[J]. J4, 2005, 35(04): 436 -0442 .
[4] 朱建伟, 赵刚, 刘博, 郭巍, 成俊. 油页岩测井识别技术及应用[J]. J4, 2012, 42(2): 289 -295 .
[5] 陈力，梁海安，张文娟，荣帆. 模糊数学方法在城市工程地质环境区划中的应用--以抚顺市城区为例[J]. J4, 2008, 38(5): 837 -0840 .
[6] 高桂梅,苏 克,王文颖,甘树才,刘招君. 吉林省桦甸油页岩中稀土元素和微量元素的研究[J]. J4, 2006, 36(6): 974 -0979 .
[7] 吴孔运，蒋忠诚，叶 晔. 不同植物群落对灰岩试块溶蚀速率的影响[J]. J4, 2007, 37(5): 967 -0971 .
[8] 周彦章，迟宝明，刘中培. 山东夏甸金矿床充水机理构造控制模式[J]. J4, 2008, 38(2): 255 -0260 .
[9] 张渊，刘连登，孙景贵，陈国华，张洪喜，闫复传,杨开春. 胶东西北部黄埠岭金矿床两期次叠加成矿[J]. J4, 2008, 38(1): 21 -0026 .
[10] 鲁程鹏, 束龙仓, 苑利波, 张蓉蓉, 黄币娟, 王彬彬. 基于示踪试验求解岩溶含水层水文地质参数[J]. J4, 2009, 39(4): 717 -721 .
版权所有 © 《吉林大学学报(地球科学版)》编辑部
地址：长春市西民主大街938号(130026)　电话：+86-431-88502374；13756165364　E-mail: xuebao1956@jlu.edu.cn
本系统由北京玛格泰克科技发展有限公司设计开发