吉林大学学报(地球科学版) ›› 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)

摘要:

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

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

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]. 吉林大学学报(地球科学版), 2026, 56(2): 673-683.
[2] 赵天硕, 宋超, , 刘财, , , , 徐雨歆.

物理信息神经网络地震走时层析成像程函方程因式分解方法 [J]. 吉林大学学报(地球科学版), 2026, 56(2): 694-702.

[3] 黄兴国, 翁央央, 韩丽. 基于广义递归卷积的孔隙黏弹地震波正演模拟[J]. 吉林大学学报(地球科学版), 2026, 56(1): 377-385.
[4] 蒋川东, 孟璐瑶, 杨麒宇, 王云志, 徐洋, 王言章, 林君. 深部金属矿产地下电磁探测技术进展与展望[J]. 吉林大学学报(地球科学版), 2025, 55(6): 2100-2119.
[5] 邵逸飞, 龚启舟, 许权 张越, 王言章, 王世隆. 基于八叉树网格的时间域航空电磁深切割地形三维正演模拟[J]. 吉林大学学报(地球科学版), 2025, 55(6): 2120-2131.
[6] 丁乾龙, 沈金松, 陈双全, 冉尚, 龙刚. 频率域传播矩阵法地震AVA多参数反演[J]. 吉林大学学报(地球科学版), 2025, 55(3): 957-969.
[7] 晋浩颖, 程大伟, 詹红兵, 杨胜科, 张晓斐, 张琳, .

考虑毛细饱水带的非承压含水层中仿Theis井流模型及解析解 [J]. 吉林大学学报(地球科学版), 2025, 55(2): 563-574.

[8] 邢慧婷, 冯晅, 刘财, 郭智奇, 逄硕, 乔汉青, .

基于贝叶斯框架的裂缝型储层频变AVO反演及多参数预测 [J]. 吉林大学学报(地球科学版), 2025, 55(2): 657-669.

[9] 龙刚, 沈金松, 苏朝阳, 冉尚. 电各向异性海底多金属硫化矿地层时域电磁响应模拟与特征分析[J]. 吉林大学学报(地球科学版), 2025, 55(1): 274-288.
[10] 吕华星, 陈兆明, 张振波, 姜大朋, 李克成, 郭伟. 机器学习高分辨融合反演在地层对比中的应用——以珠江口盆地开平凹陷开平A构造带为例[J]. 吉林大学学报(地球科学版), 2025, 55(1): 289-297.
[11] 郑轩, 宗兆云, 付亚群, 骆坤. 海域油气盖层叠前地震反演预测方法及应用[J]. 吉林大学学报(地球科学版), 2024, 54(6): 2128-2141.
[12] 张科, 窦松江, 高欢欢, 石石, 李涛, 汤戈, 杨冰.  频谱恢复和波形反演在大宁—吉县区块薄储层预测中的应用[J]. 吉林大学学报(地球科学版), 2024, 54(6): 2142-2153.
[13] 雷松达, 王显祥, 刘遂明. 复杂地质条件下浅海区电性源瞬变电磁法三维响应特征[J]. 吉林大学学报(地球科学版), 2024, 54(3): 1016-1030.
[14] 李瑞友, 白细民, 张勇, 汪靖, 朱亮, 丁小辉, 李广. 基于小波包分解与GA优化BP神经网络的瞬变电磁反演[J]. 吉林大学学报(地球科学版), 2024, 54(3): 1003-1015.
[15] 顾观文, 王顺吉, 李桐林, 武晔, 许志河. 不同供电方向激励对多极距组合中梯装置三维激电探测效果的影响#br#[J]. 吉林大学学报(地球科学版), 2024, 54(1): 292-309.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 吴远坤, 刘成林, 于春勇. 松辽盆地双城断陷深层原油成藏模式[J]. 吉林大学学报(地球科学版), 2024, 54(5): 1443 -1456 .
[2] 马荣,石建省,刘继朝. 人工内分泌网络模型在水文地质参数研究中的应用[J]. 吉林大学学报(地球科学版), 2013, 43(3): 914 -921 .
[3] 李宁, 王成文. 东北及邻区晚古生代地层接触关系与佳-蒙地块的形成和演化[J]. 吉林大学学报(地球科学版), 2017, 47(5): 1331 -1340 .
[4] 胡大千,初凤友,姚 杰. 中太平洋YJA海山富钴结壳矿物组成与元素地球化学[J]. J4, 2006, 36(01): 32 -0037 .
[5] 姜 雪, 程日辉,于民凤. 裂谷地层的气候和构造控制:Zscape模型分析与在松辽盆地北安断陷的应用[J]. J4, 2006, 36(01): 54 -0059 .
[6] 李春柏,张新涛,刘 立,任延广,孟 鹏. 布达特群热流体活动及其对火山碎屑岩的改造作用--以海拉尔盆地贝尔凹陷为例[J]. J4, 2006, 36(02): 221 -0226 .
[7] 孟宪纲,薄万举,刘志广,刘勇,畅柳,李朝柱,王子平. 芦山7.0级地震与巴颜喀拉块体中东段的活动性[J]. 吉林大学学报(地球科学版), 2014, 44(5): 1705 -1711 .
[8] 陈欢庆, 梁淑贤, 舒治睿, 邓晓娟, 彭寿昌. 冲积扇砾岩储层构型特征及其对储层开发的控制作用——以准噶尔盆地西北缘某区克下组冲积扇储层为例[J]. 吉林大学学报(地球科学版), 2015, 45(1): 13 -24 .
[9] 贾大成,邢立新, 潘 军, M. J. van Bergen, H. van Roermund. 伊通上地幔剪切带捕虏体中富铝尖晶石的地球化学特征[J]. J4, 2006, 36(04): 497 -502 .
[10] 谢忠雷,杨佰玲,包国章,董德明. 茶园土壤不同形态镍的含量及其影响因素[J]. J4, 2006, 36(04): 599 -604 .