The airborne electromagnetic discovery of the Detour zinc‐copper‐silver deposit, northwestern Québec

Geophysics ◽  
1981 ◽  
Vol 46 (9) ◽  
pp. 1278-1290 ◽  
Author(s):  
L. E. Reed

In June 1974, a diamond drill operated for Selco Mining Corp. intersected zinc‐copper sulfides in Brouillan Township in northwestern Québec. To date, two bodies have been outlined. These bodies were discovered during a ground follow‐up of a Mark VI Input® electromagnetic (EM) survey. The Input survey covered an area selected on the basis of regional geology and local outcrops of acid volcanic rocks. Conductors were identified that appeared to be associated with potentially favorable geology. They were selected for ground follow‐up. One was the discovery zone. The airborne responses over the zone were less encouraging than those often observed over highly conductive massive sulfides. The low apparent conductivity‐thickness (5 mhos) was suggestive of conductive overburden. However, the character of the profiles suggested a bedrock source. Ground geophysical confirmation identified a drill target. Subsequent to the discovery, more intensive geophysical surveys, both ground and airborne, were carried out. The best EM response suggested a confined source within a much larger mineralized halo. Weaker ground EM response from the halo correlated with the early channel response of the Input system. An airborne EM survey conducted in 1958 over the same area identified both conductive zones. However, they were not followed up. Only with later advances in exploration philosophy, geologic appreciation, and instrumentation were the conductive zones recognized as viable exploration targets.

2013 ◽  
Vol 1 (1) ◽  
pp. T101-T112 ◽  
Author(s):  
Nasreddine Bournas ◽  
Ellen Clements ◽  
Rob Hearst

Airborne electromagnetic, ground direct current and induced polarization (DCIP), and magnetotelluric (MT) surveys have extensively been used in mining exploration and more particularly for the exploration of base metal mineralization. The continuous development of geophysical techniques with advances in the instrumentation and signal processing and the recent development of robust 3D inversion algorithms make possible the detection and accurate delineation of deep-seated mineralization of economic interest. Recently, a deep-penetrating TITAN-24 DCIP and MT survey was conducted over the Silver Queen project area, located in British Columbia, Canada in two phases (2011 and 2012, respectively) by Quantec Geoscience Ltd. on behalf of New Nadina Explorations Ltd. for the exploration of porphyry-style polymetallic mineralization. The ground survey was carried out as a follow-up to the helicopter-borne z-axis tipper electromagnetic survey flown during the spring of 2011 by Geotech Ltd. with the aim to delineate favorable areas for the exploration of porphyry-style deposits. A deep-seated significant zone of anomalously high chargeability occurring in coincidence with a conductive zone was detected by the ground DCIP and MT survey. Drill-testing based on the 3D inversion results of the data led to the discovery of a new significant deep-seated porphyry-style mineralization. The discovery drillhole contained visible disseminated to semimassive sulphide mineralization, gold and molybdenite over 350 m for a total drillhole depth of approximately 800 m and occurs in association with a significant mineralized stockwork zone open at depth.


Geophysics ◽  
1969 ◽  
Vol 34 (5) ◽  
pp. 729-738 ◽  
Author(s):  
P. H. Nelson ◽  
D. B. Morris

The secondary magnetic field induced by a time‐domain, airborne EM system is calculated by transforming the tabulated mutual impedances of two magnetic dipoles above an earth of homogeneous or layered resistivity structure. The computational procedure is extended to produce response curves useful in interpreting data from a particular system, the Barringer Input system. It is demonstrated that the apparent resistivity can be estimated through use of the receiver channel ratios, a method which is independent of absolute system calibration. Layered earth calculations indicate to what extent conductive overburden cases can be readily distinguished, in terms of the conductivity‐thickness parameter, but separate interpretation of layer resistivity and thickness will require an amplitude‐calibrated flight system.


Geophysics ◽  
2007 ◽  
Vol 72 (4) ◽  
pp. F189-F195 ◽  
Author(s):  
Changchun Yin ◽  
Greg Hodges

The traditional algorithms for airborne electromagnetic (EM) inversion, e.g., the Marquardt-Levenberg method, generally run only a downhill search. Consequently, the model solutions are strongly dependent on the starting model and are easily trapped in local minima. Simulated annealing (SA) starts from the Boltzmann distribution and runs both downhill and uphill searches, rendering the searching process to easily jump out of local minima and converge to a global minimum. In the SA process, the calculation of Jacobian derivatives can be avoided because no preferred searching direction is required as in the case of the traditional algorithms. We apply SA technology for airborne EM inversion by comparing the inversion with a thermodynamic process, and we discuss specifically the SA procedure with respect to model configuration, random walk for model updates, objective function, and annealing schedule. We demonstrate the SA flexibility for starting models by allowing the model parameters to vary in a large range (far away from the true model). Further, we choose a temperature-dependent random walk for model updates and an exponential cooling schedule for the SA searching process. The initial temperature for the SA cooling scheme is chosen differently for different model parameters according to their resolvabilities. We examine the effectiveness of the algorithm for airborne EM by inverting both theoretical and survey data and by comparing the results with those from the traditional algorithms.


Geophysics ◽  
1996 ◽  
Vol 61 (1) ◽  
pp. 66-73 ◽  
Author(s):  
Richard S. Smith ◽  
Jan Klein

Airborne induced‐polarization (IP) measurements can be obtained with standard time‐domain airborne electromagnetic (EM) equipment, but only in the limited circumstances when the ground is sufficiently resistive that the normal EM response is small and when the polarizability of the ground is sufficiently large that the IP response can dominate the EM response. Further, the dispersion in conductivity must be within the bandwidth of the EM system. One example of what is hypothesized to be IP effects are the negative transients observed on a GEOTEM® survey in the high arctic of Canada. The dispersion in conductivity required to explain the data is very large, but is not inconsistent with some laboratory measurements. Whether the dispersion is caused by an electrolytic or dielectric polarization is not clear from the limited ground follow‐up, but in either case the polarization can be considered to be induced by eddy currents associated with the EM response of the ground. If IP effects are the cause of the negative transients in the GEOTEM data, then the data can be used to estimate the polarizabilities in the area.


Geophysics ◽  
2002 ◽  
Vol 67 (2) ◽  
pp. 492-500 ◽  
Author(s):  
James E. Reid ◽  
James C. Macnae

When a confined conductive target embedded in a conductive host is energized by an electromagnetic (EM) source, current flow in the target comes from both direct induction of vortex currents and current channeling. At the resistive limit, a modified magnetometric resistivity integral equation method can be used to rapidly model the current channeling component of the response of a thin-plate target energized by an airborne EM transmitter. For towed-bird transmitter–receiver geometries, the airborne EM anomalies of near-surface, weakly conductive features of large strike extent may be almost entirely attributable to current channeling. However, many targets in contact with a conductive host respond both inductively and galvanically to an airborne EM system. In such cases, the total resistive-limit response of the target is complicated and is not the superposition of the purely inductive and purely galvanic resistive-limit profiles. Numerical model experiments demonstrate that while current channeling increases the width of the resistive-limit airborne EM anomaly of a wide horizontal plate target, it does not necessarily increase the peak anomaly amplitude.


Geophysics ◽  
1984 ◽  
Vol 49 (8) ◽  
pp. 1301-1312 ◽  
Author(s):  
G. T. DeMoully ◽  
A. Becker

Recent improvements in equipment quality make it possible to increase the usefulness of airborne electromagnetic (EM) systems in areas of moderate electrical conductivity for the purpose of constructing simple electrical property maps which can be related to surficial geology. This application of airborne electromagnetics may be demonstrated and evaluated using Barringer/Questor Mark VI Input® survey results in places where independent verifications of the airborne data interpretation are available. For this purpose we have developed a set of computer algorithms which read digitally recorded Input data and interpret them automatically in terms of a simple electrical section that is defined by a single conductive layer whose thickness, conductivity, and subsurface depth are determined from the data. Because this technique is formally based on a one‐dimensional, three‐layer, three‐parameter, horizontally stratified earth model, it is only applicable in regions where the surficial formations are mildly dipping and the conductive layer is covered by, and rests on, highly resistive materials. The interpretation method is illustrated by three field examples. At the first field survey site, in Alberta, Canada, airborne EM survey data are used to map the depth of the interface between coarse and clayey sands. Data from a second survey site, this time in the Western USA, are interpreted to yield the section of a subsurface valley filled with conductive clay. The final example, taken from British Columbia, Canada, involves the mapping of all the three parameters for a weathered volcanic unit.


Geophysics ◽  
2021 ◽  
pp. 1-56
Author(s):  
Aaron Davis

Airborne geophysical surveys routinely collect data along traverse lines at sample spacing distances that are two or more orders of magnitude less than between line separations. Grids and maps interpolated from such surveys can produce aliasing; features that cross flight lines can exhibit boudinage or string-of-beads artefacts. Boudinage effects can be addressed by novel gridding methods. Following developments in geostatistics, a non-stationary nested anisotropic gridding scheme is proposed that accommodates local anisotropy in survey data. Computation is reduced by including anchor points throughout the interpolation region that contain localised anisotropy information which is propagated throughout the survey area with a smoothing kernel. Additional anisotropy can be required at certain locations in the region to be gridded. A model selection scheme is proposed that employs Laplace approximations for determining whether increased model complexity is supported by the surrounding data. The efficacy of the method is shown using a synthetic data set obtained from satellite imagery. A pseudo geophysical survey is created from the image and reconstructed with the method above. Two case histories are selected for further elucidation from airborne geophysical surveys conducted in Western Australia. The first example illustrates improvement in gridding the depth of palaeochannels interpreted from along-line conductivity-depth models of a regional airborne electromagnetic survey in the Mid-West. The second example shows how improvements can be made in producing grids of aeromagnetic data and inverted electrical conductivity from an airborne electromagnetic survey conducted in the Pilbara. In both case histories, nested anisotropic kriging reduces the expression of boudinage patterns and sharpens cross-line features in the final gridded products permitting increased confidence in interpretations based on such products.


Geophysics ◽  
2010 ◽  
Vol 75 (4) ◽  
pp. WA163-WA177 ◽  
Author(s):  
Julian Vrbancich ◽  
Peter Fullagar ◽  
Richard Smith

A floating transient electromagnetic (TEM) system (“sea ring”) simulating a low-altitude helicopter airborne electromagnetic (AEM) system was constructed to test the accuracy of the AEM method for measuring water depth and estimating sediment thickness in shallow coastal waters. A square transmitter loop [Formula: see text], plus concentric inner and outer receiver loops, was strung from masts supported by the circular sea-ring base. Data were stacked over periods from 1 to approximately [Formula: see text] and with loop heights ranging between approximately 5.5 and [Formula: see text] above sea level. The towed sea ring provides a stable platform at a known fixed altitude in calm waters. We have undertaken modeling to investigate the effect of vertical and horizontal displacements of the loops, and to compare circular and square loopgeometries, in proximity to the sea surface. With relatively long stacking times, as long as approximately [Formula: see text], the uncertainty in altitude can be reduced to very low levels. The sea ring has been deployed near Port Lincoln, Australia, in an area with known bathymetry, seawater conductivity [Formula: see text], and sediment thickness (from marine seismic). Initial 1D inversion of raw sea-ring data highlighted significant instrument calibration errors. Empirical correction factors were defined to achieve agreement between measured and modeled TEM responses at selected control points. These corrections were then applied to the entire data set. The 1D inversion of corrected sea-ring data predicted seawater depths within approximately [Formula: see text] of known depths down to [Formula: see text]. These results provide an upper limit to the accuracy that can be expected from airborne EM systems for shallow bathymetry using current technology.


Geophysics ◽  
2012 ◽  
Vol 77 (5) ◽  
pp. E317-E327 ◽  
Author(s):  
Terence Kratzer ◽  
James C. Macnae

A major impediment in the path toward airborne induced polarization (IP) is an effective method to quantify data from inductive sources, such as those used in airborne electromagnetic systems. We modeled inductive IP using a combination of Warburg and exponential decay models as a basis for fitting electromagnetic data from ground time-domain electromagnetic (TEM) and airborne versatile TEM (VTEM) surveys. Observed decays were deconvolved into electromagnetic and IP constituents by constrained least-squares fitting of basis functions modified to account for transmitter waveforms. The method was confirmed through synthetic modeling of 2D and 3D structures, and when applied to ground TEM or airborne TEM data, obtained an estimate of apparent chargeability at each station or fiducial. In the case of a VTEM survey in Africa, the apparent chargeabilities mapped graphitic sediments and provided spatially consistent indications of clay concentrations. A limitation on this airborne IP for airborne applications is motion noise, which places a lower limit on usable base frequency and begins to significantly affect the signal at the later delay times, when IP effects are most visible.


2020 ◽  
Vol 224 (1) ◽  
pp. 543-557
Author(s):  
Thomas M Hansen

SUMMARY Probabilistic inversion methods, typically based on Markov chain Monte Carlo, exist that allow exploring the full uncertainty of geophysical inverse problems. The use of such methods is though limited by significant computational demands, and non-trivial analysis of the obtained set of dependent models. Here, a novel approach, for sampling the posterior distribution is suggested based on using pre-calculated lookup tables with the extended rejection sampler. The method is (1) fast, (2) generates independent realizations of the posterior, and (3) does not get stuck in local minima. It can be applied to any inverse problem (and sample an approximate posterior distribution) but is most promising applied to problems with informed prior information and/or localized inverse problems. The method is tested on the inversion of airborne electromagnetic data and shows an increase in the computational efficiency of many orders of magnitude as compared to using the extended Metropolis algorithm.


Sign in / Sign up

Export Citation Format

Share Document