Geophysical inversion contributions to mineral exploration: lessons from the Footprints project

2019 ◽  
Vol 56 (5) ◽  
pp. 525-543 ◽  
Author(s):  
Marc A. Vallée ◽  
William A. Morris ◽  
Stéphane Perrouty ◽  
Robert G. Lee ◽  
Ken Wasyliuk ◽  
...  

Magnetic and gravity inversions are used to create 2D or 3D models of the magnetic susceptibility and density, respectively, using potential field data. Unconstrained inversions generate an output based on mathematical constraints imposed by the inversion algorithm. Constrained inversions integrate lithological, structural, and petrophysical information in the inversion process to produce more geologically meaningful results. This study analyses the validity of this assertion in the context of the NSERC-CMIC Mineral Exploration Footprints project. Unconstrained and constrained geophysical inversions were computed for three mining sites: a gold site (Canadian Malartic, Québec), a copper site (Highland Valley, British Columbia), and a uranium site (Millennium – McArthur River, Saskatchewan). After initially computing unconstrained inversions, constrained inversions were developed using physical property measurements, which directly link geophysics to geology, and lithological boundaries extracted from an interpreted geological model. While each derived geological model is consistent with the geophysical data, each site exhibited some magnetic complexity that confounded the inversion. The gold site includes regions with a strong magnetic signature that masks the more weakly magnetic zone, thereby hiding the magnetic signature associated with the ore body. Initial unconstrained inversions for the copper site yielded solutions with invalid depth extent. A consistency between the constrained model and the geological model is reached with iterative changes to the depth extent of the model. At the uranium site, the observed magnetic signal is weak, but the inversion provided some insights that could be interpreted in terms of an already known complexly folded geological model.

Minerals ◽  
2021 ◽  
Vol 11 (11) ◽  
pp. 1195
Author(s):  
Daniel Bombardieri ◽  
Mark Duffett ◽  
Andrew McNeill ◽  
Matthew Cracknell ◽  
Anya Reading

Over the last two decades, Mineral Resources Tasmania has been developing regional 3D geological and geophysical models for prospective terranes at a range of scales and extents as part of its suite of precompetitive geoscience products. These have evolved in conjunction with developments in 3D modeling technology over that time. Commencing with a jurisdiction-wide 3D model in 2002, subsequent modeling projects have explored a range of approaches to the development of 3D models as a vehicle for the better synthesis and understanding of controls on ore-forming processes and prospectivity. These models are built on high-quality potential field data sets. Assignment of bulk properties derived from previous well-constrained geophysical modeling and an extensive rock property database has enabled the identification of anomalous features that have been targeted for follow-up mineral exploration. An aspect of this effort has been the generation of uncertainty estimates for model features. Our experience is that this process can be hindered by models that are too large or too detailed to be interrogated easily, especially when modeling techniques do not readily permit significant geometric changes. The most effective 3D modeling workflow for insights into mineral exploration is that which facilitates the rapid hypothesis testing of a wide range of scenarios whilst satisfying the constraints of observed data.


Geophysics ◽  
2020 ◽  
Vol 85 (5) ◽  
pp. K13-K26 ◽  
Author(s):  
Esmaeil Eshaghi ◽  
Anya M. Reading ◽  
Michael Roach ◽  
Mark Duffett ◽  
Daniel Bombardieri ◽  
...  

The Heazlewood-Luina-Waratah area is a prospective region for minerals in northwest Tasmania, Australia, associated with historically important ore deposits related to the emplacement of granite intrusions and/or ultramafic complexes. The geology of the area is poorly understood due to the difficult terrain and dense vegetation. We have constructed an initial high-resolution 3D geologic model of this area using constraints from geologic maps and geologic and geophysical cross sections. This initial model is improved upon by integrating results from 3D geometry and physical property inversion of potential field (gravity and magnetic) data, petrophysical measurements, and updated field mapping. Geometry inversion reveals that the Devonian granites in the south are thicker than previously thought, possibly connecting to deep sources of mineralization. In addition, we identified gravity anomalies to the northeast that could be caused by near-surface granite cupolas. A newly discovered ultramafic complex linking the Heazlewood and Mount Stewart Ultramafic Complexes in the southwest also has been modeled. This implies a greater volume of ultramafic material in the Cambrian successions and points to a larger obducted component than previously thought. The newly inferred granite cupolas and ultramafic complexes are targets for future mineral exploration. Petrophysical property inversion reveals a high degree of variation in these properties within the ultramafic complexes indicating a variable degree of serpentinization. Sensitivity tests suggest maximum depths of 2–3 km for the contact aureole that surrounds major granitic intrusions in the southeast, whereas the Heazlewood River complex is likely to have a deeper source up to 4 km. We have demonstrated the value of adding geologic and petrophysical constraints to 3D modeling for the purpose of guiding mineral exploration. This is particularly important for the refinement of geologic structures in tectonically complex areas that have lithology units with contrasting magnetic and density characteristics.


2019 ◽  
Vol 219 (3) ◽  
pp. 1698-1716 ◽  
Author(s):  
M Malovichko ◽  
A V Tarasov ◽  
N Yavich ◽  
M S Zhdanov

SUMMARY This paper presents a feasibility study of using the controlled-source frequency-domain electromagnetic (CSEM) method in mineral exploration. The method has been widely applied for offshore hydrocarbon exploration; however, nowadays this method is rarely used on land. In order to conduct this study, we have developed a fully parallelized forward modelling finite-difference (FD) code based on the iterative solver with contraction-operator preconditioner. The regularized inversion algorithm uses the Gauss–Newton method to minimize the Tikhonov parametric functional with the Laplacian-type stabilizer. A 3-D parallel inversion code, based on the iterative finite-difference solver with the contraction-operator preconditioner, has been evaluated for the solution of the large-scale inverse problems. Using the computer simulation for a synthetic model of Sukhoi Log gold deposit, we have compared the CSEM method with the conventional direct current sounding and the CSEM survey with a single remote transmitter. Our results suggest that, a properly designed electromagnetic survey together with modern 3-D inversion could provide detailed information about the geoelectrical structure of the mineral deposit.


2014 ◽  
Vol 644-650 ◽  
pp. 2670-2673
Author(s):  
Jun Wang ◽  
Xiao Hong Meng ◽  
Fang Li ◽  
Jun Jie Zhou

With the continuing growth in influence of near surface geophysics, the research of the subsurface structure is of great significance. Geophysical imaging is one of the efficient computer tools that can be applied. This paper utilize the inversion of potential field data to do the subsurface imaging. Here, gravity data and magnetic data are inverted together with structural coupled inversion algorithm. The subspace (model space) is divided into a set of rectangular cells by an orthogonal 2D mesh and assume a constant property (density and magnetic susceptibility) value within each cell. The inversion matrix equation is solved as an unconstrained optimization problem with conjugate gradient method (CG). This imaging method is applied to synthetic data for typical models of gravity and magnetic anomalies and is tested on field data.


2021 ◽  
Author(s):  
Mark Jessell

<p>In geological settings characterised by folded and faulted strata, and where good field data exist, we have been able to automate a large part of the 3D modelling process directly from the raw geological database (maps, bedding orientations and drillhole data). The automation is based upon the deconstruction of the geological maps and databases into positional, gradient and spatial and temporal topology information, and the combination of deconstructed data into augmented inputs for 3D geological modelling systems, notably LoopStructural and GemPy.</p><p>When we try to apply this approach to more complex terranes, such as greenstone belts, we come across two types of problem:</p><ul><li>1) Insufficient structural data, since the more complexly deformed the geology, the more we need to rely on secondary structural information, such as fold axial traces and vergence to ‘solve’ the structures. Unfortunately these types of data are not always stored in national geological databases. One approach to overcoming this is to analyse the simpler (i.e. bedding) data to try and estimate the secondary information automatically.</li> </ul><p> </p><ul><li>2) The available information is unsuited to the logic of the modelling system. Most modern modelling platforms assume the knowledge of a chronostratigraphic hierarchy, however, especially in more complexly deformed regions, a lithostratigraphy may be all that is available. Again a pre-processing of the map and stratigraphic information may be possible to overcome this problem.</li> </ul><p>This presentation will highlight the progress that has been made, as well as the road-blocks to universal automated 3D geological model construction.</p><p> </p><p>We acknowledge the support of the MinEx CRC and the Loop: Enabling Stochastic 3D Geological Modelling (LP170100985) consortia. The work has been supported by the Mineral Exploration Cooperative Research Centre whose activities are funded by the Australian Government's Cooperative Research Centre Programme. This is MinEx CRC Document 2020/xxx.</p><p> </p>


Minerals ◽  
2020 ◽  
Vol 10 (8) ◽  
pp. 663 ◽  
Author(s):  
Sebeom Park ◽  
Yosoon Choi

Over the past decade, unmanned aerial vehicles (UAVs) have been used in the mining industry for various applications from mineral exploration to mine reclamation. This study aims to review academic papers on the applications of UAVs in mining by classifying the mining process into three phases: exploration, exploitation, and reclamation. Systematic reviews were performed to summarize the results of 65 articles (June 2010 to May 2020) and outline the research trend for applying UAVs in mining. This study found that UAVs are used at mining sites for geological and structural analysis via remote sensing, aerial geophysical survey, topographic surveying, rock slope analysis, working environment analysis, underground surveying, and monitoring of soil, water, ecological restoration, and ground subsidence. This study contributes to the classification of current UAV applications during the mining process as well as the identification of prevalent UAV types, data acquired by sensors, scales of targeted areas, and styles of flying control for the applications of UAVs in mining.


Geophysics ◽  
2012 ◽  
Vol 77 (3) ◽  
pp. E213-E225 ◽  
Author(s):  
Gianluca Fiandaca ◽  
Esben Auken ◽  
Anders Vest Christiansen ◽  
Aurélie Gazoty

Time-domain-induced polarization has significantly broadened its field of reference during the last decade, from mineral exploration to environmental geophysics, e.g., for clay and peat identification and landfill characterization. Though, insufficient modeling tools have hitherto limited the use of time-domain-induced polarization for wider purposes. For these reasons, a new forward code and inversion algorithm have been developed using the full-time decay of the induced polarization response, together with an accurate description of the transmitter waveform and of the receiver transfer function, to reconstruct the distribution of the Cole-Cole parameters of the earth. The accurate modeling of the transmitter waveform had a strong influence on the forward response, and we showed that the difference between a solution using a step response and a solution using the accurate modeling often is above 100%. Furthermore, the presence of low-pass filters in time-domain-induced polarization instruments affects the early times of the acquired decays (typically up to 100 ms) and has to be modeled in the forward response to avoid significant loss of resolution. The developed forward code has been implemented in a 1D laterally constrained inversion algorithm that extracts the spectral content of the induced polarization phenomenon in terms of the Cole-Cole parameters. Synthetic examples and field examples from Denmark showed a significant improvement in the resolution of the parameters that control the induced polarization response when compared to traditional integral chargeability inversion. The quality of the inversion results has been assessed by a complete uncertainty analysis of the model parameters; furthermore, borehole information confirm the outcomes of the field interpretations. With this new accurate code in situ time-domain-induced polarization measurements give access to new applications in environmental and hydrogeophysical investigations, e.g., accurate landfill delineation or on the relation between Cole-Cole and hydraulic parameters.


Geophysics ◽  
2010 ◽  
Vol 75 (1) ◽  
pp. L1-L11 ◽  
Author(s):  
Yaoguo Li ◽  
Sarah E. Shearer ◽  
Matthew M. Haney ◽  
Neal Dannemiller

Three-dimensional (3D) inversion of magnetic data to recover a distribution of magnetic susceptibility has been successfully used for mineral exploration during the last decade. However, the unknown direction of magnetization has limited the use of this technique when significant remanence is present. We have developed a comprehensive methodology for solving this problem by examining two classes of approaches and have formulated a suite of methods of practical utility. The first class focuses on estimating total magnetization direction and then incorporating the resultant direction into an inversion algorithm that assumes a known direction. The second class focuses on direct inversion of the amplitude of the magnetic anomaly vector. Amplitude data depend weakly upon magnetization direction and are amenable to direct inversion for the magnitude of magnetization vector in 3D subsurface. Two sets of high-resolution aeromagnetic data acquired for diamond exploration in the Canadian Arctic are used to illustrate the methods’ usefulness.


Sign in / Sign up

Export Citation Format

Share Document