-
摘要: 由于传统的阻尼最小二乘法只适合于模型较少的简单模型,因此当介质的层数较多时,反演就会受到多解性的影响,有时甚至出现不收敛的情况,并且反演十分耗时。为此,使用正则化思想引入模型约束进行反演,且正则化因子通过计算每次迭代的数据目标函数和模型目标函数自适应得到,使反演能够稳定地进行;引入拟牛顿法来更新雅可比矩阵,大大缩短反演所需要的时间,通过典型的3层与多层理论模型的反演试算,证明了拟牛顿法自适应正则化反演算法对初始模型的要求不高,拟合效果好,收敛速度快,适应性强,体现了良好的稳定性和可靠性。Abstract: Because the damping least square method is only suitable for a simple model, when the media has multiple layers, inversion is effected by multiple solutions, sometimes even no convergence occurs, and the inversion is very time-consuming. For this purpose, regularization idea is used to introduce the model-constrained inversion, in the cause of iteration, the adaptive regular factor of the regularized inversion algorithm is calculated adaptively based on the relation of the data objective function and the model objective function; it can result in a stable convergence in the iteration course of the inversion. In this paper, the inversion used quasi-Newton method to update Jacobian matrix, it reduced greatly the time required for the inversion, and the typical spreadsheet three and multi-theoretical model as an example to inversion is proved less demanding on the initial model of the adaptive quasi-Newton regularized inversion algorithm. The inversion result also show that the adaptive regularized inversion is steady and reliable and has good fitting effect, fast convergence and strong adaptability.
-
-
[1] 牛之琏. 时间域电磁法原理[M]. 长沙:中南大学出版社,2007. [2] 闫述,薛国强,陈明生. 大回线源瞬变电磁响应理论研究回顾及展望[J]. 地球物理学进展,2011,26(3):941-947. YAN Shu,XUE Guoqiang,CHEN Mingsheng. Review and perspective of theoretical study on large-loop TEM response[J]. Progress in Geophysics,2011,26(3):941-947.
[3] 彭仲秋,段新力,毕武,等. 瞬变电磁“烟圈”反演实现[J]. 物探化探计算技术,2012,34(5):573-575. PENG Zhongqiu,DUAN Xinli,BI Wu,et al. TEM "smoking rings" inversion soft-wave ware development[J]. Computing Techniques for Geophysical and Geochemical Exploration,2012,34(5):573-575.
[4] 陆桂福,张强,黄力军. 回线源瞬变电磁测深快速反演方法[J]. 物探化探计算技术,2005,27(4):289-291. LU Guifu,ZHANG Qiang,HUANG Lijun. A rapid inversion technique for loop source transient electromagnetic method[J]. Computing Techniques for Geophysical and Geochemical Exploration,2005,27(4):289-291.
[5] 李锋平,杨海燕,邓居智,等. 地面瞬变电磁法一维烟圈反演技术研究[J]. 地球物理学进展,2016,31(2):688-694. LI Fengping,YANG Haiyan,DENG Juzhi,et al. One-dimensional smoke ring inversion technology of ground transient electromagnetic method[J]. Progress in Geophysics,2016,31(2):688-694.
[6] 蒋邦远. 磁源TEM测深一维反演经验式及应用效果[J]. 物探与化探,2000,24(2):99-104. JIANG Bangyuan. A rapid 1-D empirical inversion formula and its application results in TEM sounding[J]. Geophysical and Geochemical Exploration,2000,24(2):99-104.
[7] 郭文达,朱希安. 烟圈反演和视纵向电导解释方法对比[J]. 煤田地质与勘探,2014,42(6):93-96. GUO Wenda,ZHU Xi'an. Comparison between the smoke loop inversion and apparent longitudinal conductivity interpretation method[J]. Coal Geology & Exploration,2014,42(6):93-96.
[8] 黄皓平,王维中. 时间域航空电磁数据的反演[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.
[9] 杨云见,王绪本,何展翔. 考虑关断时间效应的瞬变电磁一维反演[J]. 物探与化探,2005,29(3):234-236. YANG Yunjian,WANG Xuben,HE Zhanxiang. 1D inversion of transient electromagnetic data in consideration of ramp time effect[J]. Geophysical and Geochemical Exploration,2005,29(3):234-236.
[10] 王园园,刘斌,王晨. 基于阻尼最小二乘法的瞬变电磁反演算法研究[J]. 电子测试,2013,1:1-4. WANG Yuanyuan,LIU Bin,WANG Chen. Study on transient electromagnetic inversion based on damped least square method[J]. Electronic Test,2013,1:1-4.
[11] 张维. 大定源瞬变电磁法一维正反演研究[D]. 长沙:中南大学,2013. [12] 唐宝山. 瞬变电磁法中心回线装置一维正反演研究[D]. 北京:中国地质大学(北京),2008. [13] 郭嵩巍. 中心回线瞬变电磁一维正反演算法研究[D]. 成都:成都理工大学,2010. [14] 徐玉聪,赵宁,秦策,等. 大定源瞬变电磁一维自适应正则化反演[J]. 地质与勘探,2015,51(2):360-365. XU Yucong,ZHAO Ning,QIN Ce,et al. One-dimensional adaptive regularization inversion of transient electromagnetic sounding with a large fixed source[J]. Geology and Exploration,2015,51(2):360-365.
[15] 毛立峰,王绪本,李文杰. 飞行高度同时反演的固定翼航空瞬变电磁一维反演[J]. 地球物理学报,2011,54(8):2136-2147. MAO Lifeng,WANG Xuben,LI Wenjie. Research on 1D inversion method of fix-wing airborne transient electromagnetic record with flight altitude inversion simultaneously[J]. Chinese Journal of Geophysics,2011,54(8):2136-2147.
[16] 毛立峰,王绪本,陈斌. 直升机航空瞬变电磁自适应正则化一维反演方法研究[J]. 地球物理学进展,2011,26(1):300-305. MAO Lifeng,WANG Xuben,CHEN Bin. Study on an adaptive regularized 1D inversion method of helicopter TEM data[J]. Progress in Geophysics,2011,26(1):300-305.
[17] 毛立峰,陈斌,吕东伟. 测高数据不准时的直升机航空瞬变电磁一维反演方法理论研究[J]. 物探与化探,2011,35(3):402-405. MAO Lifeng,CHEN Bin,LYU Dongwei. The 1D inversion method of airborne TEM for inaccurate altimeter data[J]. Geophysical and Geochemical Exploration,2011,35(3):402-405.
[18] 覃庆炎,罗威,周洪生. 中心回线瞬变电磁自适应正则化反演[J]. 科学技术与工程,2014,14(13):7-10. QIN Qingyan,LUO Wei,ZHOU Hongsheng. Inversion display of aria used in in-loop transient electromagnetic sounding[J]. Science Technology and Engineering,2014,14(13):7-10.
[19] 翁爱华,高玲玲,刘云鹤. 矩形大定源瞬变测深数据模型约束反演[J]. 物探化探计算技术,2008,30(4):314-316. WENG Aihua,GAO Lingling,LIU Yunhe. Inversion of TEM data from rectangular loop on a layeredearth model[J]. Computing Techniques for Geophysical and Geochemical Exploration,2008,30(4):314-316.
[20] 陈小斌,赵国泽,汤吉,等. 大地电磁自适应正则化反演算法[J]. 地球物理学报,2005,48(4):937-946. CHEN Xiaobin,ZHAO Guoze,TANG Ji,et al. An adaptive regularized inversion algorithm for magnetotelluric data[J]. Chinese Journal of Geophysics,2005,48(4):937-946.
[21] CHRISTIANSEN A V,AUKEN E. Optimizing a layered and laterally constrained 2D inversion of resistivity data using Broyden's update and 1D derivatives[J]. Journal of Applied Geophysics,2004,56(4):247-261.
[22] 胡祖志,胡祥云,何展翔. 大地电磁非线性共轭梯度拟三维反演[J]. 地球物理学报,2006,49(4):1226-1234. HU Zuzhi,HU Xiangyun,HE Zhanxiang. Pseudo-Three-Dimensional magnetotelluric inversion using nonlinear conjugate gradients[J]. Chinese Journal of Geophysics,2006,49(4):1226-1234.
[23] 胡祖志. 大地电磁拟三维反演研究[D]. 武汉:中国地质大学(武汉),2005. [24] 米萨克N.纳比吉安. 勘查地球物理电磁法[M]. 赵经祥,等,译. 北京:地质出版社,1992. [25] GUPTASARMA D,SINGH B. New digital linear filters for Hankel J0 and J1 transforms[J]. Geophysical Prospecting,1997,45(5):745-762.
[26] 朴化荣,殷长春. 利用G-S逆拉氏变换法计算瞬变测深正演问题[J]. 物化探计算技术,1987,9(4):295-302. PIAO Huarong,YIN Changchun. Calculation of transient EM sounding using the Gaver-Stehfest inverse Laplace transform method[J]. Computing Techniques for Geophysical and Geochemical Exloration,1987,9(4):295-302.
[27] 罗延钟,昌彦君. G-S变换的快速算法[J]. 地球物理学报,2000,43(5):684-690. LUO Yanzhong,CHANG Yanjun. A rapid algorithm for G-S transform[J]. Chinese Journal of Geophysics,2000,43(5):684-690.
[28] FARQUHARSON C G,OLDENBURG D W. Inversion of time-domain electromagnetic data for a horizontally layered earth[J]. Geophysical Journal International,1993,114(3):433-442.
计量
- 文章访问数: 131
- HTML全文浏览量: 58
- PDF下载量: 23