Document Type : Original Research Paper

Authors

1 Faculty of Mining, Petroleum and Geophysics, Shahrood University of Technology, Shahrood, Iran.

2 Shahrood University of Technology, Faculty of Mining, Petroleum and Geophysics, Shahrood, Iran.

3 Federal Institute for Geosciences and Natural Resources (BGR), Hannover, Germany.

4 Isfahan University of Technology (IUT), Faculty of Mining Engineering, Isfahan, Iran.

Abstract

It is about 30 years that Helicopter electromagnetic (HEM) surveys are being used for rapid mineral and ground water exploration, environmental investigations and also geological mapping in extensive areas. Despite this, one of the most important problems in using obtained data from the surveys is accurate interpretation of the data. Otherwise, there will be no beneficial results while spending high costs. Thus the interpretation of the data is as old as the surveys. Several experts have tried to improve the interpretation of HEM data and they have achieved great successes. Almost the results of all these surveys are presented as resistivity (or conductivity)-depth sections. To reach this target, the first step is to solve the electromagnetic induction integral equation. As solving this integral is not possible using analytical methods, several numerical methods such as Laplace transformation, Hankel transformation and Jacobi-Matrix methods have been suggested for the solution of the integral, and different approaches have been presented with each method by various authorities. One of the most important solution methods is fast Hankel transformation. In this paper, it is attempted to use this method for finally obtaining resistivity-depth sections. For solving the induction equation by this method, we need the kernel function of the integral and weighting coefficients that replace the Bessel function in the integral. For this, first we use the Guptasarma-Singh method. Then results of this method are corrected and evaluated. Then, these results will be analyzed and tested with two synthetic models in addition to presenting the results of inverse modeling. Finally, by adding new parameter named α0 to induction equation, we will clearly see an improvement in the results of inverse modeling. Meanwhile, the problem of singularity that occurs at high frequencies is almost removed.
 

Keywords

References
Fraser, D. C., 1978- Resistivity mapping with an airborne multi-coil electromagnetic system, Geophysics, 43: 144–172.
Guptasarma, D. & Singh, B., 1997- New digital filters for Hankel J0 and J1 transforms, Geophysical Prospecting, 45: 745-762.
Mundry, E., 1984- On the interpretation of airborne electromagnetic data for the two-layer case: Geophysical Prospecting, 32: 336–346.
Sengpiel, K. P. & Siemon, B., 1998- Examples of 1D inversion of multifrequency AEM data from 3D resistivity distributions, Exploration Geophysics, 29: 133–141.
Sengpiel, K. P. & Siemon, B., 2000- Advanced inversion methods for airborne electromagnetic, Exploration. Geophysics, 65: 1983–1992.
Sengpiel, K. P., 1990- Theoretical and practical aspects of ground-water exploration using airborne electromagnetic techniques. In: Fitterman, D.V. (ed.), Proceedings of the USGS Workshop on Developments and Applications of Modern Airborne Electromagnetic Surveys. 1987: Golden, Co, 149–154. USGS Bulletin 1925, Denver, Co.
Siemon, B., 2001- Improved and new resistivity –depth profiles for helicopter electromagnetic data. Journal of Applied Geophysics, 46: 65– 76.
Siemon, B., 2007- Levelling of helicopter-borne frequency-domain electromagnetic data, Journal of Applied Geophysics, xx: xxx-xxx.
Wait, J. R., 1982- Geo-Electromagnetism, Academic Press Inc, New York.