An efficient hybrid scheme for fast and accurate inversion of airborne transient electromagnetic data
Anders Vest Christiansen 1 4 Esben Auken 1 Casper Kirkegaard 1 Cyril Schamper 2 Giulio Vignoli 3
+ Author Affiliations
- Author Affiliations
1 Department of Geoscience, Aarhus University, Aarhus 8000, Denmark.
2 Unité Mixte de Recherche 7619 Sisyphe, Université Pierre et Marie Curie, 75252 Paris, France.
3 Geological Surveys of Denmark and Greenland, Aarhus 8000, Denmark.
4 Corresponding author. Email: anders.vest@geo.au.dk
Exploration Geophysics 47(4) 323-330 https://doi.org/10.1071/EG14121
Submitted: 25 November 2014 Accepted: 8 July 2015 Published: 13 August 2015
Abstract
Airborne transient electromagnetic (TEM) methods target a range of applications that all rely on analysis of extremely large datasets, but with widely varying requirements with regard to accuracy and computing time. Certain applications have larger intrinsic tolerances with regard to modelling inaccuracy, and there can be varying degrees of tolerance throughout different phases of interpretation. It is thus desirable to be able to tune a custom balance between accuracy and compute time when modelling of airborne datasets. This balance, however, is not necessarily easy to obtain in practice. Typically, a significant reduction in computational time can only be obtained by moving to a much simpler physical description of the system, e.g. by employing a simpler forward model. This will often lead to a significant loss of accuracy, without an indication of computational precision.
We demonstrate a tuneable method for significantly speeding up inversion of airborne TEM data with little to no loss of modelling accuracy. Our approach introduces an approximation only in the calculation of the partial derivatives used for minimising the objective function, rather than in the evaluation of the objective function itself. This methodological difference is important, as it introduces no further approximation in the physical description of the system, but only in the process of iteratively guiding the inversion algorithm towards the solution. By means of a synthetic study, we demonstrate how our new hybrid approach provides inversion speed-up factors ranging from ~3 to 7, depending on the degree of approximation. We conclude that the results are near identical in both model and data space. A field case confirms the conclusions from the synthetic examples: that there is very little difference between the full nonlinear solution and the hybrid versions, whereas an inversion with approximate derivatives and an approximate forward mapping differs significantly from the other results.
Key words: AEM, approximate Jacobian, hybrid minimisation, large dataset inversion.
References
Allard, M., 2007, On the origin of the HTEM species, in B. Milkereit, ed., Proceedings of Exploration 07: Fifth Decennial International Conference on Mineral Exploration: Decennial Mineral Exploration Conferences, 355–374.Auken, E., and Christiansen, A. V., 2004, Layered and laterally constrained 2D inversion of resistivity data: Geophysics, 69, 752–761
| Layered and laterally constrained 2D inversion of resistivity data:Crossref | GoogleScholarGoogle Scholar |
Auken, E., Christiansen, A. V., Jacobsen, L., and Sørensen, K. I., 2008, A resolution study of buried valleys using laterally constrained inversion of TEM data: Journal of Applied Geophysics, 65, 10–20
| A resolution study of buried valleys using laterally constrained inversion of TEM data:Crossref | GoogleScholarGoogle Scholar |
Auken, E., Christiansen, A. V., Kirkegaard, C., Fiandaca, G., Schamper, C., Behroozmand, A. A., Binley, A., Nielsen, E., Effersø, F., Christensen, N. B., Sørensen, K. I., Foged, N., and Vignoli, G., 2014, An overview of a highly versatile forward and stable inverse algorithm for airborne, ground-based and borehole electromagnetic and electric data: Exploration Geophysics, ,
| An overview of a highly versatile forward and stable inverse algorithm for airborne, ground-based and borehole electromagnetic and electric data:Crossref | GoogleScholarGoogle Scholar |
Broyden, C. G., 1965, A class of methods for solving nonlinear simultaneous equations: Mathematics of Computation, 19, 577–593
| A class of methods for solving nonlinear simultaneous equations:Crossref | GoogleScholarGoogle Scholar |
Christensen, N. B., 1997, Electromagnetic subsurface imaging. A case for adaptive Born approximation: Surveys in Geophysics, 18, 477–510
| Electromagnetic subsurface imaging. A case for adaptive Born approximation:Crossref | GoogleScholarGoogle Scholar |
Christensen, N. B., Reid, J. E., and Halkjaer, M., 2009, Fast, laterally smooth inversion of airborne time-domain electromagnetic data: Near Surface Geophysics, 7, 599–612
| Fast, laterally smooth inversion of airborne time-domain electromagnetic data:Crossref | GoogleScholarGoogle Scholar |
Christiansen, A. V., and Auken, E., 2004, Optimizing a layered and laterally constrained 2D inversion of resistivity data using Broyden’s update and 1D derivatives: Journal of Applied Geophysics, 56, 247–261
| Optimizing a layered and laterally constrained 2D inversion of resistivity data using Broyden’s update and 1D derivatives:Crossref | GoogleScholarGoogle Scholar |
Christiansen, A. V., and Auken, E., 2012, A global measure for depth of investigation: Geophysics, 77, WB171–WB177
| A global measure for depth of investigation:Crossref | GoogleScholarGoogle Scholar |
Christiansen, A. V., Auken, E., and Viezzoli, A., 2011, Quantification of modeling errors in airborne TEM caused by inaccurate system description: Geophysics, 76, F43–F52
| Quantification of modeling errors in airborne TEM caused by inaccurate system description:Crossref | GoogleScholarGoogle Scholar |
Cox, L. H., Wilson, G. A., and Zhdanov, M. S., 2010, 3D inversion of airborne electromagnetic data using a moving footprint: Exploration Geophysics, 41, 250–259
| 3D inversion of airborne electromagnetic data using a moving footprint:Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BC3cXhsFyjt7jN&md5=74793b40de54b147501e8ce64633db02CAS |
Haber, E., Ascher, U. M., and Oldenburg, D., 2000, On optimization techniques for solving nonlinear inverse problems: Inverse Problems, 16, 1263–1280
| On optimization techniques for solving nonlinear inverse problems:Crossref | GoogleScholarGoogle Scholar |
Harbo, M. S., Pedersen, J., Johnsen, R., and Petersen, K., eds., 2011, Groundwater in a future climate: Danish Ministry of the Environment Nature Agency.
Hestenes, M. R., and Stiefel, E., 1952, Methods of conjugate gradients for solving linear systems: Journal of Research of the National Bureau of Standards, 49, 409–436
| Methods of conjugate gradients for solving linear systems:Crossref | GoogleScholarGoogle Scholar |
Huang, H., and Fraser, D. C., 1996, The differential parameter method for multifrequency airborne resistivity mapping: Geophysics, 61, 100–109
| The differential parameter method for multifrequency airborne resistivity mapping:Crossref | GoogleScholarGoogle Scholar |
Huang, H., and Palacky, G. J., 1991, Damped least-squares inversion of time-domain airborne EM data based on singular value decomposition: Geophysical Prospecting, 39, 827–844
| Damped least-squares inversion of time-domain airborne EM data based on singular value decomposition:Crossref | GoogleScholarGoogle Scholar |
Jørgensen, F., Scheer, W., Thomsen, S., Sonnenborg, T. O., Hinsby, K., Wiederhold, H., Schamper, C., Roth, B., Kirsch, R., and Auken, E., 2012, Transboundary geophysical mapping of geological elements and salinity distribution critical for the assessment of future sea water intrusion in response to sea level rise: Hydrology and Earth System Sciences, 16, 1845–1862
| Transboundary geophysical mapping of geological elements and salinity distribution critical for the assessment of future sea water intrusion in response to sea level rise:Crossref | GoogleScholarGoogle Scholar |
Kirkegaard, C., and Auken, E., 2015, A parallel, scalable and memory efficient inversion code for very large scale airborne EM surveys: Geophysical Prospecting, 63, 495–507
| A parallel, scalable and memory efficient inversion code for very large scale airborne EM surveys:Crossref | GoogleScholarGoogle Scholar |
Kirkegaard, C., Sonnenborg, T. O., Auken, E., and Flemming, J., 2011, Salinity distribution in heterogeneous coastal aquifers mapped by airborne electromagnetics: Vadose Zone Journal, 10, 125–135
| Salinity distribution in heterogeneous coastal aquifers mapped by airborne electromagnetics:Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BC3MXltVaisrc%3D&md5=ad5dd0702f57313a8497c728148dc539CAS |
Kirkegaard, C., Foged, N., Auken, E., Christiansen, A. V., and Sørensen, K. I., 2012, On the value of including x-component data in 1D modeling of electromagnetic data from helicopterborne time domain systems in horizontally layered environments: Journal of Applied Geophysics, 84, 61–69
| On the value of including x-component data in 1D modeling of electromagnetic data from helicopterborne time domain systems in horizontally layered environments:Crossref | GoogleScholarGoogle Scholar |
Legault, J. M., Witter, J. B., Berardelli, P., Lombardo, S., and Orta, M., 2011, Recent ZTEM airborne AFMAG EM survey results over Reese River and other geothermal test areas: Geothermal Resources Council 2011 Annual Meeting, 23–26 October 2011, San Diego, California, 35, 879–884.
Ley-Cooper, A. Y., Viezzoli, A., Guillemoteau, J., Vignoli, G., Macnae, J., Cox, L., and Munday, T., 2015, Airborne electromagnetic modelling options and their consequences in target definition: Exploration Geophysics, 46, 74–84
| Airborne electromagnetic modelling options and their consequences in target definition:Crossref | GoogleScholarGoogle Scholar |
Loke, M. H., and Barker, R. D., 1996, Rapid least squares inversion of apparent resistivity pseudosections by a quasi-Newton method: Geophysical Prospecting, 44, 131–152
| Rapid least squares inversion of apparent resistivity pseudosections by a quasi-Newton method:Crossref | GoogleScholarGoogle Scholar |
Macnae, J., King, A., Stolz, N., Osmakoff, A., and Blaha, A., 1998, Fast AEM data processing and inversion: Exploration Geophysics, 29, 163–169
| Fast AEM data processing and inversion:Crossref | GoogleScholarGoogle Scholar |
Marquardt, D., 1963, An algorithm for least squares estimation of nonlinear parameters: SIAM Journal on Applied Mathematics, 11, 431–441
| An algorithm for least squares estimation of nonlinear parameters:Crossref | GoogleScholarGoogle Scholar |
Minsley, B., Cox, L. H., Brodie, R., Wilson, G. A., Abraham, J. D., and Zhdanov, M. S., 2012, A comparison of AEM inversion methods for discontinuous permafrost at Fort Yukon, Alaska: Symposium on the Application of Geophysics to Environmental and Engineering Problems, SAGEEP, Tuscon, Arizona.
Oldenburg, D. W., and Ellis, R. G., 1991, Inversion of geophysical data using an approximate inverse mapping: Geophysical Journal International, 105, 325–353
| Inversion of geophysical data using an approximate inverse mapping:Crossref | GoogleScholarGoogle Scholar |
Pfaffhuber, A. A., Monstad, S., and Rudd, J., 2009, Airborne electromagnetic hydrocarbon mapping in Mozambique: Exploration Geophysics, 40, 237–245
| Airborne electromagnetic hydrocarbon mapping in Mozambique:Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BD1MXhtFCktb%2FN&md5=2cb93cbb1c950dc41209bc05e0615c77CAS |
Plessix, R. E., 2006, A review of the adjoint-state method for computing the gradient of a functional with geophysical applications: Geophysical Journal International, 167, 495–503
| A review of the adjoint-state method for computing the gradient of a functional with geophysical applications:Crossref | GoogleScholarGoogle Scholar |
Romdhane, A., Grandjean, G., Brossier, R. M., Rejiba, F., Operto, S., and Vineux, J., 2011, Shallow structure characterization by 2D elastic full-waveform inversion: Geophysics, 76, R81–R93
| Shallow structure characterization by 2D elastic full-waveform inversion:Crossref | GoogleScholarGoogle Scholar |
Sengpiel, K. P., and Siemon, B., 2000, Advanced inversion methods for airborne electromagnetic exploration: Geophysics, 65, 1983–1992
| Advanced inversion methods for airborne electromagnetic exploration:Crossref | GoogleScholarGoogle Scholar |
Siemon, B, Christiansen, A. V., and Auken, E, 2009, A review of helicopter-borne electromagnetic methods for groundwater exploration: Near Surface Geophysics, 7, 629–646
Sørensen, K. I., and Auken, E., 2004, SkyTEM – a new high-resolution helicopter transient electromagnetic system: Exploration Geophysics, 35, 191–199
Torres-Verdín, C., Druskin, V. L., Fang, S., Knizhnerman, L. A., and Malinverno, A., 2000, A dual-grid nonlinear inversion technique with applications to the interpretation of dc resistivity data: Geophysics, 65, 1733–1745
| A dual-grid nonlinear inversion technique with applications to the interpretation of dc resistivity data:Crossref | GoogleScholarGoogle Scholar |
Viezzoli, A., Christiansen, A. V., Auken, E., and Sørensen, K. I., 2008, Quasi-3D modeling of airborne TEM data by spatially constrained inversion: Geophysics, 73, F105–F113
| Quasi-3D modeling of airborne TEM data by spatially constrained inversion:Crossref | GoogleScholarGoogle Scholar |
Viezzoli, A., Munday, T., Auken, E., and Christiansen, A. V., 2010, Accurate quasi 3D versus practical full 3D inversion of AEM data – the Bookpurnong case study: Preview, 149, 23–31
Ward, S. H., and Hohmann, G. W., 1988, Electromagnetic theory for geophysical applications, in M. N. Nabighian, ed., Electromagnetic methods in applied geophysics: Society of Exploration Geophysicists (SEG), 131–311.
Wilson, G. A., Raiche, A., and Sugeng, F., 2006, 2.5D inversion of airborne electromagnetic data: Exploration Geophysics, 37, 363–371
Zhdanov, M. S., 2002, Geophysical inverse theory and regularization problems: Elsevier.