There is a clear demand to increase detection depths in the context of raw material exploration programs. Semi-airborne electromagnetic (semi-AEM) methods can address these demands by combining the advantages of powerful transmitters deployed on the ground with efficient helicopter-borne mapping of the magnetic field response in the air. The penetration depth can exceed those of classic airborne EM systems because low frequencies and large transmitter-receiver offsets can be realized in practice. A novel system has been developed that combines high-moment horizontal electric bipole transmitters on the ground with low-noise three-axis induction coil magnetometers, a three-axis fluxgate magnetometer, and a laser gyro inertial measurement unit integrated within a helicopter-towed airborne platform. The attitude data are used to correct the time series for motional noise and subsequently to rotate into an earth-fixed reference frame. In a second processing step, and as opposed to existing semi-AEM systems, we transform the data into the frequency domain and estimate the complex-valued transfer functions between the received magnetic field components and the synchronously recorded injection current by regression analysis. This approach is similar to the procedure used in controlled-source EM. For typical source bipole moments of 20–40 kAm and for rectangular current waveforms with a fundamental frequency of approximately 10 Hz, we can estimate reliable three-component (3C) transfer functions in the frequency range from 10 to 5000 Hz over a measurement area of 4×5 km2 for a single source installation. The system has the potential to be used for focused exploration of deep targets.

We implemented a novel multi-resolution grid approach to direct current resistivity (DCR) modeling in 3-D. The multi-resolution grid was initially developed to solve the electromagnetic forward problem and helped to improve the modeling efficiency. In the DCR forward problem, the distribution of the electric potentials in the subsurface is estimated. We consider finite-difference staggered grid discretization, which requires fine grid resolution to accurately model electric potentials around the current electrodes and complex model geometries near the surface. Since the potential variations attenuate with depth, the grid resolution can be decreased correspondingly. The conventional staggered grid fixes the horizontal grid resolution that extends to all layers. This leads to over-discretization and therefore unnecessary high computational costs (time and memory). The non-conformal multi-resolution grid allows the refinement or roughening for the grid’s horizontal resolution with depth, resulting in a substantial reduction of the degrees of freedom, and subsequently, computational requirements. In our implementation, the coefficient matrix maintains its symmetry, which is beneficial for using the iterative solvers and solving the adjoint problem in inversion. Through comparison with the staggered grid, we have found that the multi-resolution grid can significantly improve the modeling efficiency without compromising the accuracy. Therefore, the multi-resolution grid allows modeling with finer horizontal resolutions at lower computational costs, which is essential for accurate representation of the complex structures. Consequently, the inversion based on our modeling approach will be more efficient and accurate.

The classical radio-magnetotelluric (RMT) method is nowadays routinely applied to various environmental, engineering, and exploration problems. The technique uses radio transmitters broadcasting in the frequency range of 10 kHz to 1 MHz, and the measurements are carried out in the far field. The well-known disadvantages of RMT are a lack of robust radio transmitters in remote areas; the absence of transmitters broadcasting below 10 kHz, which limits the penetration depth; and a possible low signal-to-noise ratio. To overcome these difficulties, controlled sources can be used (controlled-source RMT [CSRMT]). We extend the CSRMT method to perform measurements not only in the far field but also in the transition zone. In CSRMT practice, it often is challenging to maintain far-field conditions for logistical reasons. Therefore, part of the measured data contains signatures of the source field, which cannot be interpreted with magnetotelluric software. In addition, the source placed directly in the survey area allows us to increase the signal-to-noise ratio and resolution. Such CSRMT in the transition zone is, in fact, a controlled-source electromagnetic method but with full impedance tensor and tipper vector transfer functions. We develop new procedures for the 3D modeling and inversion of the tensor radio-frequency data measured in the transition zone of two perpendicular horizontal electric dipole sources. In this case, the geometry of the source must be considered in the forward modeling. The developed modeling and inversion software is tested on a synthetic 3D model. The 3D resistivity models derived from the real data confirm the geologic settings and are consistent with the available borehole information. Therefore, we conclude that the CSRMT approach extended to include the source field is feasible and that the developed procedures are reliable.

This study compares the efficiency of 3-D transient electromagnetic forward modeling schemes on the multi-resolution grid for various modeling scenarios. We developed time-domain finite-difference modeling based on the explicit scheme earlier. In this work, we additionally implement 3-D transient electromagnetic forward modeling using the backward Euler implicit scheme. The iterative solver is used for solving the system of equations and requires a proper initial guess that has significant effect on the convergence. The standard approach usually employs the solution of a previous time step as an initial guess, which might be too conservative. Instead, we test various initial guesses based on the linear extrapolation or linear combination of the solutions from several previous steps. We build up the implicit scheme forward modeling on the multi-resolution grid, which allows for the adjustment of the horizontal resolution with depth, hence improving the performance of the forward operator. Synthetic examples show the implicit scheme forward modeling using the linearly combined initial guess estimate on the multi-resolution grid additionally reduces the run time compared to the standard initial guess approach. The result of comparison between the implicit scheme developed here with the previously developed explicit scheme shows that the explicit scheme modeling is more efficient for more conductive background models often found in environmental studies. However, the implicit scheme modeling is more suitable for the simulation with highly resistive background models, usually occurring in mineral exploration scenarios. Thus, the inverse problem can be solved using more efficient forward solution depending on the modeling setup and background resistivity.

We discuss the implementation of multi-resolution framework to 3-D Direct Current (DC) problem. Commonly used staggered (SG) grid fixes the horizontal grid resolution for all depths. Thus, employing the fine horizontal resolution may lead to an over-discretised forward problem, subsequently affecting the performance of the inversion. We implemented a novel multi-resolution (MR) grid approach to the 3-D DC modeling and inversion problem, which allows adjustment of the horizontal resolution with depth. By using finer resolution for the near-surface regions, MR grid can ensure the modeling accuracy and describe the shallow features in the inversion model as well. The ability to use relatively coarser horizontal resolution for the deeper regions reduces the computation costs compare to the SG grid modeling. As a result, modeling and inversion can be accelerated several times by solving a smaller problem. Our grid resembles non-conformal rectangular grid, which commonly used in finite-elements modelling.