A secant infinite approximation algorithm for the full-time apparent resistivity in transient electromagnetic method
ZENG Qingning1,2, LUO Ying2, LIU Shuai2, LONG Chao1,2
1. Key Laboratory of Cognitive Radio and Information Processing, Ministry of Education, Guilin University of Electronic Technology, Guilin, Guangxi 541004, China; 2. College of Information and Telecommunications, Guilin University of Electronic Technology, Gui-lin, Guangxi 541004, China
Abstract:In the geology exploration with the transient electromagnetic method (TEM),the apparent resistivity calculation is quite important as the appa-rent resistivity may intuitively indicate underground geological bodies.However,conventional approximation formulas usually lead a large error in the apparent resistivity calculation.According to the characteristics of implicit functions of apparent resistivity for half-space central-loop TEM,we propose in this paper a secant infinite approximation algorithm for the full-time apparent resistivity in the transient electromagnetic method,in which the secant method is used to calculate the roots of non-linear functions.The proposed algorithm avoids complicated derivative calculations of implicit functions.Simulation and real data examples show that the proposed algorithm has faster convergence than the bisection algorithm and can calculate with arbitrary accuracy the full-time apparent resistivity for half-space central-loop TEM.
薛国强,李貅,底青云.瞬变电磁法正反演问题研究进展[J].地球物理学进展,2008,23(4):1165-1172.XUE Guoqiang,LI Xiu,DI Qingyun.Research progress in TEM forward modeling and inversion calculation[J]. Progress in Geophysics,2008,23(4):1165-1172.
[3]
Key K.1D inversion of multicomponent,multifrequency marine CSEM data:methodology and synthetic stu-dies for resolving thin resistive layers[J].Geophy-sics,2009,74(2):F9-F20.
[4]
Flores C,Peralta-Ortega S A.Induced polarization with in-loop transient electromagnetic soundings:a case study of mineral discrimination at El Arcoporphyry Copper[J].Journal of Applied Geophysics,2009,68(3):423-436.
[5]
Dennis Z R,Cull J P.Transient electromagnetic surveys for the measurement of near-surface anisotropy[J].Journal of Applied Geophysics,2012,76(1):64-73.
[6]
Sudha,Tezkan B,Siemon B.Appraisal of a new 1D weighted joint inversion of ground based and helicopter-borne electromagnetic data[J].Geophysical Prospecting,2014,62(3):597-614.
[7]
郭嵩巍,王绪本.瞬变电磁全区视电阻率数值计算方法研究[J].物探化探计算技术,2010,32(5):500-507.GUO Songwei,WANG Xuben.Study on numerical calculation methods of transient electromagnetic all-time apparent resistivity[J].Computing Techniques for Geophysical and Geochemical Exploration,2010,32(5):500-507.
[8]
李建慧,朱自强,刘树才,等.模拟退火法计算大定源瞬变电磁法的视电阻率[J].石油地球物理勘探,2011,46(1):138-142.LI Jianhui,ZHU Ziqiang,LIU Shucai,et al.Calculation of apparent resistivity in large fixed loop TEM by simulated annealing algorithm[J].Oil Geophysical Prospecting,2011,46(1):138-142.
[9]
许洋铖,吴燕清,崔靭,等.基于"双烟圈"理论的瞬变电磁全期视电阻率计算[J].物探与化探,2014,38(3):516-521.XU Yangcheng,WU Yanqing,CUI Ren,et al.The calculation of full time apparent resistivity for tran-sient electromagnetic survey based on "double smoking ring" theory[J].Geophysical and Geochemical Exploration,2014,38(3):516-521.
[10]
陈明生,许洋鋮.瞬变电磁场关断效应及全期视电阻率的普适算法[J].煤田地质与勘探,2017,45(4):131-134.CHEN Mingsheng,XU Yangcheng.Suitable algorithm of transient electromagnetic turn-off effects and whole-stage apparent resistivity[J].Coal Geology & Exploration,2017,45(4):131-134.
[11]
白登海,Maxwell A M,卢健,等.时间域瞬变电磁法中心方式全程视电阻率的数值计算[J].地球物理学报,2003,46(5):697-704.BAI Denghai,Maxwell A M,LU Jian,et al.Numerical calculation of all-time apparent resistivity for the central loop transient electromagnetic method[J].Chinese Journal of Geophysics,2003,46(5):697-704.
[12]
熊彬.大回线瞬变电磁法全区视电阻率的逆样条插值计算[J].吉林大学学报(地球科学版),2003,35(4):515-519.XIONG Bin.Inverse spline interpolation for the calculation of all-time resistivity for the large-loop transient electromagnetic method[J].Journal of Jilin University (Earth Science Edition),2003,35(4):515-519.
陈清礼,严良俊,付志红,等.长偏移距瞬变电磁法全区视电阻率的二分搜索数值算法[J].石油地球物理勘探,2009,44(6):779-782.CHEN Qingli,YAN Liangjun,FU Zhihong,et al.The all-time apparent resistivity bisearch numerical algorithm based on long offset transient electromagnetic method[J].Oil Geophysical Prospecting,2009,44(6):779-782.
[15]
张成范,翁爱华,孙世栋,等.计算矩形大定源回线瞬变电磁测深全区视电阻率.吉林大学学报(地球科学版),2009,39(4):755-758.ZHANG Chengfan,WENG Aihua,SUN Shidong,et al.Computation of whole-time apparent resistivity of large rectangular loop[J].Journal of Jilin University (Earth Science Edition),2009,39(4):755-758.
[16]
杨海燕,岳建华.地下瞬变电磁法全区视电阻率核函数算法[J].中国矿业大学学报,2013,42(1):83-87.YANGHaiyan,YUE Jianhua.Kernel function algorithm of all-time apparent resistivity used in underground transient electromagnetic method[J].Journal of China University of Mining & Technology,2013,42(1):83-87.
[17]
李文尧,晏冲为.中心回线瞬变电磁法全期视电阻率的二分法求解[J].昆明理工大学学报(自然科学版),2013,38(2):26-33.LI Wenyao,YAN Chongwei.Bisection algorithm of central loop time domain transient electromagnetism all time apparent resistivity[J].Journal of Kunming University of Science and Technology(Natural Science Edition),2013,38(2):26-33.
[18]
曾庆宁,张法全,张海如.井下瞬变电磁法中全期视电阻率的牛顿求解法[J].物探化探计算技术,2014,36(4):401-404.ZENG Qingning,ZHANG Faquan,ZHANG Hairu.Newton algorithm for solving full-time apparent resistivity of underground-mine transient electromagnetic method[J].Computing Techniques for Geophysical and Geochemical Exploration,2014,36(4):401-404.
[19]
苏朱刘,胡文宝.中心回线方式瞬变电磁测深虚拟全区视电阻率和一维反演方法[J].石油物探,2002,41(2):216-221.SU Zhuliu,HU Wenbao.Pseudo-full-region apparent resistivity and its one-dimensional inversion for center-loop-line configuration TEM data[J].Geophysical Prospecting for Petroleum,2002,41(2):216-221.
[20]
John H,Mathews,Kurtis D F. Numerical Methods Using MATLAB[M].Prentice Hall,2010,51-56.
[21]
刘玲,王正盛.数值计算方法[M].北京:科学出版社,2010,21-45.
[22]
严良俊,胡文宝,陈清礼,等.长偏移距瞬变电磁测深的全区视电阻率求取及快速反演方法[J].石油地球物理勘探,1999,34(5):532-539.YAN Liangjun,HU Wenbao,CHEN Qingli,et al.The estimation and fast inversion of all-time apparent resistivities in long-offset transient electromagnetic sounding[J].Oil Geophysical Prospecting,1999,34(5):532-539.
[23]
赵越,许枫,李貅,等.中心回线海底三维瞬变电磁响应规律分析[J].石油地球物理勘探,2017,52(5):1093-1102.ZHAO Yue,XU Feng,LI Xiu,et al.Characteristic of 3D TEM response on seafloor with the central loop configuration[J].Oil Geophysical Prospecting,2017,52(5):1093-1102.
[24]
刘晓波,李艳莉,李磊,等.地震属性和瞬变电磁技术在煤田勘探中的应用[J].石油地球物理勘探,2018,53(2):218-223.LIU Xiaobo,LI Yanli,LI Lei,et al.Joint application of seismic attributes and transient electromagnetic data in the coalfield exploration[J].Oil Geophysical Prospecting,2018,53(2):218-223.