Calculation of all-time apparent resistivity of large loop transient electromagnetic method with very fast simulated annealing

Abstract

In large loop transient electromagnetic method (TEM), the late time apparent resistivity formula cannot truly reflect the geoelectric model, thus it needs to define the all-time apparent resistivity with the position information of measuring point. Utilizing very fast simulated annealing (VFSA) to fit the theoretical electromagnetic force (EMF) and measured EMF could obtain the all-time apparent resistivity of the measuring points in rectangular transmitting loop. The selective cope of initial model of VFSA could be confirmed by taking the late time apparent resistivity of transient electromagnetic method as the prior information. For verifying the correctness, the all-time apparent resistivities of the geoelectric models were calculated by VFSA and dichotomy, respectively. The results indicate that the relative differences of apparent resistivities calculated by these two methods are within 3%. The change of measuring point position has little influence on the tracing pattern of all-time apparent resistivity. The first branch of the curve of all-time apparent resistivity is close to the resistivity of the first layer medium and the last branch is close to the resistivity of the last layer medium, which proves the correctness of the arithmetics proposed.

Keywords

very fast simulated annealing (VFSA) / large loop transient electromagnetic method / rectangular loop / all-time apparent resistivity

References

[1] JiangB.-yuan.Applied near zone magnetic source transient electromagnetic exploration [M], 1998, Beijing, Geological Press

[2] XiongBin.. Inverse spline interpolation for the calculation of all-time resistivity for the larger-loop transient electromagnetic method [J]. Journal of Jilin University: Earth Science Edition, 2005, 35(4): 515-519

[3] LiJ.-p., LiT.-l., ZhaoX.-f., LiangT.-mu.. Study on the TEM all-time apparent resistivity of arbitrary shape loop source over the layered medium [J]. Progress in Geophysics, 2007, 22(6): 1777-1780

[4] WangH.-jun.. Time domain transient electromagnetic all time apparent resistivity translation algorithm [J]. Chinese J Geophys, 2008, 51(6): 1936-1942

[5] ZhangC.-f., WengA.-h., SunS.-d., DongR.-chun.. Computation of whole-time apparent resistivity of large rectangular loop [J]. Journal of Jilin University: Earth Science Edition, 2009, 39(4): 755-758

[6] WengA.-h., LiuY.-h., ChenY.-l., DongR.-c., LiaoX.-d., JiaD.-yu.. Computation of transient electromagnetic field from a rectangular loop over stratified earths [J]. Chinese J Geophys, 2010, 53(3): 646-650

[7] XiZ.-z., LiuJ., LongX., HouH.-tao.. Three-component measurement in transient electromagnetic method [J]. Journal of Central South University: Science and Technology, 2010, 41(1): 272-276

[8] FuC.-m., DiQ.-y., WangM.-yue.. Calculate electromagnetic fields in stratified medium with layer-matrix [J]. Chinese J Geophys, 2010, 53(1): 177-188

[9] LiJ.-h., LiuS.-c., ZhuZ.-q., LuG.-y., ZhaoX.-b., ZhaoS.-qiu.. Relationship between electromagnetic field and magnetic field’s symmetric excited by rectangular loop [J]. Journal of Central South University: Science and Technology, 2010, 41(2): 638-642

[10] NewmanG. A., CommerM.. New advances in three dimensional transient electromagnetic inversion [J]. Geophys J Int, 2005, 160: 5-32

[11] TangJ.-t., WangF.-y., RenZ.-y., GuoR.-wen.. 3-D direct current resistivity forward modeling by adaptive multigrid finite element method [J]. Journal of Central South University of Technology, 2010, 17: 587-592

[12] KaikkonenP., SharmaS. P.. A comparison of performances of linearized and global nonlinear 2-D inversions of VLF and VLF-R electromagnetic data [J]. Geophysics, 2001, 66(3): 462-475

[13] RydenN., ParkC. B.. Fast simulated annealing inversion of surface waves on pavement using phase-velocity spectra [J]. Geophysics, 2006, 71(4): R49-R58

[14] MaX. Q.. Simultaneous inversion of prestack seismic data for rock properties using simulated annealing [J]. Geophysics, 2002, 67(6): 1877-1885

[15] MisralS., SacchiM. D.. Global optimization with model-space preconditioning application to AVO inversion [J]. Geophysics, 2008, 73(5): R71-R82

[16] RoyL., SenM. K., DonaldD. B., StoffaP. L., RichterT. G.. Inversion and uncertainty estimation of gravity data using simulated annealing: An application over Lake Vostok, East Antarctica [J]. Geophysics, 2005, 70(1): J1-J12

[17] DafflonB., IrvingJ., HolligerK.. Simulated-annealingbased conditional simulation for the local-scale characterization of heterogeneous aquifers [J]. Journal of Applied Geophysics, 2009, 68: 60-70

[18] YinC. C., HodgesG.. Simulated annealing for airborne EM inversion [J]. Geophysics, 2007, 72(4): F189-F195

[19] KnightJ. H., RaicheA. P.. Transient electromagnetic calculations using the Gaver-Stehfest inverse laplace transform method [J]. Geophysics, 1982, 47(1): 47-50

[20] EverettM. E.. Transient electromagnetic response of a loop source over a rough geological medium [J]. Geophys J Int, 2009, 177: 421-429

[21] GuptasarmaD., SinghB.. New digital linear filters for Hankel J0 and J1 transforms [J]. Geophysical Prospecting, 1997, 45: 745-762