吉林大学学报(地球科学版) ›› 2016, Vol. 46 ›› Issue (1): 254-261.doi: 10.13278/j.cnki.jjuese.201601302

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

时间域航空电磁数据加权横向约束反演

殷长春, 邱长凯, 刘云鹤, 蔡晶   

  1. 吉林大学地球探测科学与技术学院, 长春 130026
  • 收稿日期:2015-05-19 出版日期:2016-01-26 发布日期:2016-01-26
  • 通讯作者: 邱长凯(1991),男,研究生,主要从事电磁法正反演和方法技术研究,E-mail:qiuck13@mails.jlu.edu.cn E-mail:qiuck13@mails.jlu.edu.cn
  • 作者简介:殷长春(1965),男,教授,国家"千人计划"特聘专家,主要从事电磁勘探理论特别是航空和海洋电磁方面的研究,E-mail:yinchangchun@jlu.edu.cn
  • 基金资助:

    国家重大科研装备研制项目(ZDYZ2012-1-03);国家自然科学基金项目(41274121);国土资源部老矿山复核项目(中地调研合同[2013]第214号)

Weighted Laterally-Constrained Inversion of Time-Domain Airborne Electromagnetic Data

Yin Changchun, Qiu Changkai, Liu Yunhe, Cai Jing   

  1. College of GeoExploration Science and Technology, Jilin University, Changchun 130026, China
  • Received:2015-05-19 Online:2016-01-26 Published:2016-01-26
  • Supported by:

    Supported by Projects on the Development of the Key Equipment of Chinese Academy of Sciences (ZDYZ2012-1-03),China Natural Science Foundation (41274121) and Ministry of Land and Resources of the Peoples Republic of China Old Mine Review Proiect(China Geological Survey Project Agreement[2013] No.214)

PDF (PC)

493

摘要：

传统的单点一维反演技术处理时间域航空电磁数据得到的电阻率和层厚度容易呈现横向不连续性,导致层界面不光滑。即使相邻测点的反演结果也会出现突变现象,与实际地质情况不符,给解释工作带来困难。笔者针对时间域航空电磁数据空间采样密集、在沉积岩地区工作时地下电性结构横向连续性较好的特点,研究了时间域航空电磁数据参数加权横向约束反演方法。反演时同时考虑数据拟合、相邻测点电阻率和层厚度横向约束以及深度横向约束,利用参数加权来调整对各层模型参数和深度横向光滑度的约束强度。通过对理论和实测数据反演并与传统的单点一维反演结果进行对比,验证了拟二维横向加权约束反演是处理时间域航空电磁数据的有效方法。横向约束反演结果电阻率和层厚度横向光滑连续,层界面清晰。引入参数加权实现对不同参数横向光滑度和连续性的制约,可进一步提高反演解的稳定性。加权因子对反演结果有较大影响,在实际应用中,应该综合考虑地质资料选择合适的加权因子。

关键词: 航空电磁, 时间域, 加权横向约束, 正演, 反演

Abstract:

The resistivity and thickness inverted by traditional single-site 1D methods for time-domain airborne electromagnetic (AEM) data frequently demonstrate sharp variations along a profile. This phenomenon occurs even between neighbor stations, resulting in non-smooth layer interfaces, which can seriously complicate the interpretation. Considering that AEM systems generally have a high spatial sampling rate, and the underground structures in areas like sedimentary formations are continuous, we use the weighted laterally-constrained inversion (WLCI) algorithm to invert the time-domain AEM data. Data fitting,lateral constraints on resistivity,thickness and depth are considered at the same time, and further the constraint strength is adjusted by using the weighting factors (chosen on geological information). Through testing the above algorithm on both synthetic and survey data, and comparing the results with those from conventional inversions, the WLCI algorithm is proved to be a very reliable and effective tool to process the time-domain AEM data. The inverted results by WLCI algorithm are smooth and continuous with clear layer interfaces. Weighting for different constrained parameters can also improve the stability of inversion. The weighting factors have an important impact on the inversion results. To choose the appropriate weighting factors, the field geological information should be taken into account.

Key words: airborne electromagnetic, time-domain, weighted laterally-constrained inversion, forward modeling, inversion

中图分类号: 

  • P631.3

[1] 雷栋,胡祥云,张素芳.航空电磁法的发展现状[J].地质找矿论丛,2006,21(1):40-44. Lei Dong, Hu Xiangyun, Zhang Sufang. Development Status of Airborne Electromagnetic[J]. Contributions to Geology and Mineral Resources Research, 2006, 21(1):40-44.

[2] 殷长春,张博,刘云鹤,等.航空电磁勘查技术发展现状及展望[J].地球物理学报,2015,58(8):2637-2653. Yin Changchun, Zhang Bo, Liu Yunhe, et al. Review on Airborne EM Technology and Developments[J]. Chinese Journal of Geophysics, 2015, 58(8):2637-2653.

[3] Macnae J, King A, Stolz N, et al. Fast AEM Data Processing and Inversion[J]. Exploration Geophysics, 1998, 29(1/2):163-169.

[4] Wolfgram P, Karlik G. Conductivity-Depth Transform of GEOTEM Data[J]. Exploration Geophysics, 1995, 26(2/3):179-185.

[5] 陈小红,段奶军.时间域航空电磁快速成像研究[J].地球物理学进展,2012,27(5):2123-2127. Chen Xiaohong, Duan Naijun. Study on Fast Imaging of Airborne Time-Domain Electromagnetic Data[J]. Progress in Geophysics, 2012, 27(5):2123-2127.

[6] Huang H, Rudd J. Conductivity-Depth Imaging of He-licopter-Borne TEM Data Based on a Pseudolayer Half-Space Model[J]. Geophysics, 2008, 73(3):F115-F120.

[7] 朱凯光,林君,韩悦慧,等.基于神经网络的时间域直升机电磁数据电导率深度成像[J].地球物理学报,2010,53(3):743-750. Zhu Kaiguang, Lin Jun, Han Yuehui, et al. Research on Conductivity Depth Imaging of Time Domain Helicopter-Borne Electromagnetic Data Based on Neural Network[J]. Chinese Journal of Geophysics, 2010, 53(3):743-750.

[8] 毛立峰.中心回线直升机TEM资料的电导率-深度成像方法[J].CT理论与应用研究,2013,22(3):429-437. Mao Lifeng. Conductivity-Depth Imaging Algorithm for Central-Loop Helicopter TEM[J]. CT Theory and Applications, 2013, 22(3), 429-437.

[9] 嵇艳鞠, 冯雪, 于明媚, 等. 基于多元线性回归的HTEM三维异常体电导率-深度识别[J]. 吉林大学学报(地球科学版), 2014, 44(5):1687-1694. Ji Yanju, Feng Xue, Yu Mingmei, et al. Conductivity-Depth Identification of HTEM 3D Anomalies Based on Multiple Linear Regression[J]. Journal of Jilin University (Earth Science Edition), 2014, 44(5):1687-1694.

[10] 黄皓平,王维中.时间域航空电磁数据的反演[J].地球物理学报,1990,33(1):87-97. Huang Haoping, Wang Weizhong. Inversion of Time-Domain Airborne Electromagnetic Data[J]. Chinese Journal of Geophysics, 1990, 33(1):87-97.

[11] 罗勇,陆从德,王宇航.时间域航空电磁一维阻尼特征参数反演方法[J].地球物理学进展,2014,29(6):2723-2729. Luo Yong, Lu Congde, Wang Yuhang. 1D Inverse Method Based on Damped Eigenparameter for Airborne Time Domain Electromagnetic Data[J]. Progress in Geophysics, 2014, 29(6):2723-2729.

[12] 强建科,李永兴,龙剑波.航空瞬变电磁数据一维Occam反演[J].物探化探计算技术,2013,35(5):501-505. Qiang Jianke, Li Yongxing, Long Jianbo. 1D Occam Inversion Method for Airborne Transient Electromagnetic Data[J]. Computing Techniques for Geophysical and Geochemical Exploration, 2013, 35(5):501-505.

[13] Constable S C, Parker R L, Constable C G. Occam's Inversion:A Practical Algorithm for Generating Smooth Models from Electromagnetic Sounding Data[J]. Geophysics, 1987, 52(3):289-300.

[14] Sattel D. Inverting Airborne Electromagnetic (AEM) Data with Zohdy's Method[J]. Geophysics, 2005, 70(4):G77-G85.

[15] 刘云鹤,殷长春.三维频率域航空电磁反演研究[J].地球物理学报,2013,56(12):4278-4287. Liu Yunhe, Yin Changchun. 3D Inversion for Frequency-Domain HEM Data[J]. Chinese Journal of Geophysics, 2013, 56(12):4278-4287.

[16] Auken E, Foged N, Sørensen K I. Model Recognition by 1-D Laterally Constrained Inversion of Resistivity Data[C]//8th EEGS-ES Meeting. Aveiro:Universidade de Aveiro, 2002.

[17] Auken E, Christiansen A V. Layered and Laterally Constrained 2D Inversion of Resistivity Data[J]. Geophysics, 2004, 69(3):752-761.

[18] Santos F A M. 1-D Laterally Constrained Inversion of EM34 Profiling Data[J]. Journal of Applied Geophysics, 2004, 56(2):123-134.

[19] Auken E, Christiansen A V, Jacobsen L, et al. La-terally Constrained 1D-Inversion of 3D TEM Data[C]//10th European Meeting of Environmental and Engineering Geophysics. Utrecht:[s. n.], 2004.

[20] Siemon B, Auken E, Christiansen A V. Laterally Constrained Inversion of Helicopter-Borne Frequency-Domain Electromagnetic Data[J]. Journal of Applied Geophysics, 2009, 67(3):259-268.

[21] 蔡晶,齐彦福,殷长春,等.频率域航空电磁数据的加权横向约束反演[J].地球物理学报,2014,57(1):953-960. Cai Jing, Qi Yanfu, Yin Changchun, et al. Weighted Laterally-Constrained Inversion of Frequency-Domain Airborne EM Data[J]. Chinese Journal of Geophy-sics, 2014, 57(1):953-960.

[22] Vallée M A, Smith R S. Inversion of Airborne Time-Domain Electromagnetic Data to a 1D Structure Using Lateral Constraints[J]. Near Surface Geophysics, 2009, 7(1):63-71.

[23] 殷长春,黄威,贲放.时间域航空电磁系统瞬变全时响应正演模拟[J].地球物理学报,2013,56(9):3153-3162. Yin Changchun, Huang Wei, Ben Fang. The Full-Time Electromagnetic Modeling for Time-Domain Airborne Electromagnetic Systems[J]. Chinese Journal of Geophysics, 2013, 56(9):3153-3162.

[24] Yin C, Smith R S, Hodges G, et al. Modeling Results of On-and Off-Time B and dB/dt for Time-Domain Airborne EM Systems[C]//70th Annual EAGE Conference and Exhibition. Rome:[s. n.], 2008:1-4.

[25] Smith R S, Lee T. Asymptotic Expansion for Calculation of Transient EM Fields Induced by a Vertical Magnetic Dipole Above a Conductive Half-Space[J]. Pure and Applied Geophysics, 2004, 161:385-397.

[26] Marquardt D. An Algorithm for Least-Squares Estimation of Nolinear Parameters[J]. Journal of the Society for Industrial and Applied Mathematics, 1963, 11(2):431-441.

[27] Chen T, Hodges G, Miles P. MULTIPULSE-High Resolution and High Power in One TDEM System[J]. Exploration Geophysics, 2015, 46(1):49-57.
[1] 郑国磊, 徐新学, 李世斌, 袁航, 马为, 叶青. 天津市重力数据反演解释[J]. 吉林大学学报(地球科学版), 2018, 48(4): 1221-1230.
[2] 王泰涵, 黄大年, 马国庆, 李野, 林松. 基于并行预处理算法的三维重力快速反演[J]. 吉林大学学报(地球科学版), 2018, 48(2): 384-393.
[3] 李建平, 翁爱华, 李世文, 李大俊, 李斯睿, 杨悦, 唐裕, 张艳辉. 基于球坐标系下有限差分的地磁测深三维正演[J]. 吉林大学学报(地球科学版), 2018, 48(2): 411-419.
[4] 曾昭发, 霍祉君, 李文奔, 李静, 赵雪宇, 何荣钦. 任意各向异性介质三维有限元航空电磁响应模拟[J]. 吉林大学学报(地球科学版), 2018, 48(2): 433-444.
[5] 刘新彤, 刘四新, 孟旭, 傅磊. 低频缺失下跨孔雷达包络波形反演[J]. 吉林大学学报(地球科学版), 2018, 48(2): 474-482.
[6] 殷长春, 卢永超, 刘云鹤, 张博, 齐彦福, 蔡晶. 多重网格准线性近似技术在三维航空电磁正演模拟中的应用[J]. 吉林大学学报(地球科学版), 2018, 48(1): 252-260.
[7] 陈辉, 尹敏, 殷长春, 邓居智. 大地电磁三维正演聚集多重网格算法[J]. 吉林大学学报(地球科学版), 2018, 48(1): 261-270.
[8] 安振芳, 张进, 张建中. 海洋三维VC观测系统优化设计[J]. 吉林大学学报(地球科学版), 2018, 48(1): 271-284.
[9] 翁爱华, 李斯睿, 杨悦, 李大俊, 李建平, 李世文. 磁电法基本原理、发展现状及前景展望[J]. 吉林大学学报(地球科学版), 2017, 47(6): 1838-1854.
[10] 王通, 王德利, 冯飞, 程浩, 魏敬轩, 田密. 三维稀疏反演多次波预测及曲波域匹配相减技术[J]. 吉林大学学报(地球科学版), 2017, 47(6): 1865-1874.
[11] 李大俊, 翁爱华, 杨悦, 李斯睿, 李建平, 李世文. 地-井瞬变电磁三维交错网格有限差分正演及响应特性[J]. 吉林大学学报(地球科学版), 2017, 47(5): 1552-1561.
[12] 陈帅, 李桐林, 张镕哲. 考虑激发极化效应的瞬变电磁一维Occam反演[J]. 吉林大学学报(地球科学版), 2017, 47(4): 1278-1285.
[13] 张代磊, 黄大年, 张冲. 基于遗传算法优化的BP神经网络在密度界面反演中的应用[J]. 吉林大学学报(地球科学版), 2017, 47(2): 580-588.
[14] 高秀鹤, 黄大年, 孙思源, 于平. 重力梯度数据协克里金三维反演确定岩脉倾向[J]. 吉林大学学报(地球科学版), 2017, 47(2): 589-596.
[15] 翁爱华, 刘佳音, 贾定宇, 杨悦, 李建平, 李亚彬, 赵祥阳. 有限长导线源频率测深有限内存拟牛顿一维反演[J]. 吉林大学学报(地球科学版), 2017, 47(2): 597-605.
Viewed
Full text


Abstract

Cited
  Shared   
  Discussed   
No Suggested Reading articles found!
版权所有 © 《吉林大学学报(地球科学版)》编辑部
地址：长春市西民主大街938号(130026)　电话：+86-431-88502374；13756165364　E-mail: xuebao1956@jlu.edu.cn
本系统由北京玛格泰克科技发展有限公司设计开发