Different geophysical methods provide information about various physical properties of rock formations and mineralization. In many cases, this information is mutually complementary. At the same time, inversion of the data for a particular survey is subject to considerable uncertainty and ambiguity as to causative body geometry and intrinsic physical property contrast. One productive approach to reducing uncertainty is to jointly invert several types of data. Non-uniqueness can also be reduced by incorporating additional information derived from available geological and/or geophysical data in the survey area to reduce the searching space for the solution. This additional information can be incorporated in the form of a joint inversion of multiphysics data. This paper presents an overview of the main ideas and principles of novel methods of joint inversion, developed over the last decade, which do not require a priori knowledge about specific empirical or statistical relationships between the different model parameters and/or their attributes. These approaches are designated as follows: (1) Gramian constraints; (2) Gramian-based structural constraints; (3) localized Gramian constraints; and (4) joint focusing constraints. We provide a short description of the mathematical foundations of each of these approaches and discuss the practical aspects of their applications in mineral exploration.

This paper discusses the physical and mathematical principles of the airborne induced polarization (IP) method. The possibility of extracting information about the IP properties of rocks from airborne survey data has become a subject of active research recently. We introduce a method for the joint inversion of the airborne EM data into the electrical conductivity and IP parameters based on the generalized effective-medium theory of induced polarization (GEMTIP). We also present the results of the inversion of the airborne EM data collected over the Echum Project Area, in Northwestern Ontario, Canada, into 3D conductivity and chargeability models. Obtaining IP physical property models from an airborne geophysical survey may result in a paradigm change in mineral exploration by pulling more information and value from airborne EM surveys.

The Saudi Arabian Glass Earth Pilot Project is a geophysical exploration program to explore the upper crust of the Kingdom for minerals, groundwater, and geothermal resources as well as strictly academic investigations. The project began with over 8000 km² of green-field area. Airborne geophysics including electromagnetic (EM), magnetics, and gravity were used to develop several high priority targets for ground follow-up. Based on the results of airborne survey, a spectral induced polarization (SIP) survey was completed over one of the prospective targets. The field data were collected with a distributed array system, which has the potential for strong inductive coupling. This was examined in a synthetic study, and it was determined that with the geometries and conductivities in the field survey, the inductive coupling effect may be visible in the data. In this study, we also confirmed that time domain is vastly superior to frequency domain for avoiding inductive coupling, that measuring decays from 50 ms to 2 s allow discrimination of time constants from 1 ms to 5 s, and the relaxation parameter C is strongly coupled to intrinsic chargeability. We developed a method to fully include all 3D EM effects in the inversion of induced polarization (IP) data. The field SIP data were inverted using the generalized effective-medium theory of induced polarization (GEMTIP) in conjunction with an integral equation-based modeling and inversion methods. These methods can replicate all inductive coupling and EM effects, which removes one significant barrier to inversion of large bandwidth spectral IP data. The results of this inversion were interpreted and compared with results of drill hole set up in the survey area. The drill hole intersected significant mineralization which is currently being further investigated. The project can be considered a technical success, validating the methods and effective-medium inversion technique used for the project.

Recent developments in large-scale geophysical inversions made it possible to invert the results of entire airborne geophysical surveys over large areas into 3D models of the subsurface. This paper presents the methods for and results of the interpretation of the data acquired by a multiphysics airborne geophysical survey in Saudi Arabia. The project involved the acquisition, processing, and interpretation of airborne electromagnetic, gravity, and magnetic geophysical data over an 8000 square kilometer area. All the collected data were carefully analyzed and inverted in 3D models of the corresponding physical properties of the subsurface, including 3D density, magnetization vector, and conductivity models. This paper summarizes the interpretation of all geophysical data sets collected during the field airborne survey. The goal of the paper is to demonstrate how the advanced 3D modeling and inversion methods can be effectively used for interpretation of multiphysics airborne survey data and to study and analyze the potential of the survey area for natural resource exploration in Saudi Arabia.

This paper develops a novel method of 3D inversion of induced polarization (IP) survey data, based on a generalized effective-medium model of the IP effect (GEMTIP). The electrical parameters of the effective-conductivity model are determined by the intrinsic petrophysical and geometrical characteristics of composite media, such as the mineralization and/or fluid content of rocks and the matrix composition, porosity, anisotropy, and polarizability of formations. The GEMTIP model of multiphase conductive media provides a quantitative tool for evaluation of the type of mineralization, and the volume content of different minerals using electromagnetic (EM) data. The developed method takes into account the nonlinear nature of both electromagnetic induction and IP phenomena and inverts the EM data in the parameters of the GEMTIP model. The goal of the inversion is to determine the electrical conductivity and the intrinsic chargeability distributions, as well as the other parameters of the relaxation model simultaneously. The recovered parameters of the relaxation model can be used for the discrimination of different rocks, and in this way may provide an ability to distinguish between uneconomic mineral deposits and zones of economic mineralization using geophysical remote sensing technology.

This paper introduces an original approach to the joint inversion of airborne electromagnetic (EM) data for three-dimensional (3D) conductivity and chargeability models using hybrid finite difference (FD) and integral equation (IE) methods. The inversion produces a 3D model of physical parameters, which includes conductivity, chargeability, time constant, and relaxation coefficients. We present the underlying principles of this approach and an example of a high-resolution inversion of the data acquired by a new active time domain airborne EM system, TargetEM, in Ontario, Canada. The new TargetEM system collects high-quality multicomponent data with low noise, high power, and a small transmitter–receiver offset. This airborne system and the developed advanced inversion methodology represent a new effective method for mineral resource exploration.