Search

Mechanistic induced polarization (IP) models describe the relationships between the intrinsic properties of geomaterials and their frequency-dependent complex conductivity spectra. However, the uncertainties associated with estimating petrophysical properties from IP data are still poorly understood. Therefore, practitioners rarely use mechanistic models to interpret actual IP data. We propose a framework for critically assessing any IP model's sensitivity and parameter estimation limitations. The framework consists of a conditional variational autoencoder (CVAE), an unsupervised Bayesian neural network specializing in data dimension reduction and generative modeling. We train the CVAE on the IP signatures of synthetic mixtures of metallic mineral inclusions in electrolyte-filled host geomaterials and describe the effect of data transformations on the model. First, the CVAE's Jacobian reveals the relative importance of each petrophysical property for generating spectral IP data. The most critical parameters are the conductivity of the host, the volumetric content of the inclusions, the characteristic length of the inclusions, and the permittivity of the host. The inclusions' diffusion coefficient, permittivity, and conductivity, as well as the host's diffusion coefficient, only have marginal importance for generative IP modeling. A parameter estimation experiment yields the standardized accuracy of petrophysical properties using various model constraints scenarios and corroborates the sensitivity analysis results. Finally, we visualize the effects of data transformations and model constraints on the petrophysical parameter space. We conclude that a common logarithm data transformation yields optimal parameter estimation results and that constraining the electrochemical properties of the geomaterial improves estimates of the characteristic length of its metallic inclusions and vice versa., 31 pages, 12 figures, 1 table. Original manuscript submitted to SEG Geophysics for review. For associated code, see https://doi.org/10.5281/zenodo.6948651

Electrical impedance spectroscopy (EIS), also known as Spectral Induced Polarization (SIP) has proven to be a successful non-destructive technique able to characterize the chemical and physical properties of a complex structure (concrete, rocks, and soil) under various environmental conditions. Because this technique aims to measure small amplitude signals, it can easily be influenced by parasitic effects which are not representative of the intrinsic properties of the investigated material. The purpose of this paper is to validate the resolution of our experimental methodology as a preliminary step to ensure an accurate measurement of the bulk complex resistivity response of concrete samples within a wide frequency range (1.43 mHz–20 kHz). A cement mortar sample is first used as a homogeneous (isotropic) material to improve the sample holder design for reducing the errors and controlling the external parameters (such as coupling effect, geometric factor, contact impedance, electrode polarization) during the laboratory measurements. Also, the performance and installation of the measuring components such as the electrodes (current and potential) and the sample holder design are assessed. The specimen dimensions and designs are selected based on a sample standard size that will be used to determine the physico-mechanical properties of concrete (compressive strength, modulus of elasticity, ultrasonic pulse velocity, length variation). The examinations are performed in parallel with the SIP measurement on the same sample. The accuracy of our measurement setup is then validated using aluminum-bronze plates bonded to the sample using a conductive gel to transmit electric current with a density less than 10 - 3 A m−2. The results showed that by controlling the sponge moisture in the Ag/AgCl non-polarizable voltage electrodes, our SIP measurement system is able to measure the phase lag with an error smaller than 1 mrad over a frequency range from 1.43 mHz to 20 kHz.

Structurally complex zones within orogenic terranes typically correspond to areas where there is interference between multiple fold generations and are known to be favorable pathways for fluid flow because of their higher permeability. In the Canadian Malartic district, gold anomalies have been linked with zones of structural complexity that have been quantified by outcrop bedding orientation measurements and calculation of bedding variance maps. In this work, historical apparent resistivity and induced polarization data in the Canadian Malartic district were reprocessed and combined with new surveys to create a compilation of inverted chargeability and resistivity, which were then interpreted together with airborne electromagnetics and outcrop structural data. The results indicate chargeability anomalies, up to five times the background value, associated with the sulfide mineral content in monzodioritic dikes that are thickened in folds and hydrothermally altered. Although the airborne apparent half-space resistivity is mostly sensitive to conductive surficial cover, the inverted ground resistivity method is sensitive to deeper structure and likely represents bedrock signal at depths greater than 25 m. Inverted ground resistivity exhibits strong anisotropy in areas of subvertical bedding, where measured resistivities can vary by up to a factor of two, over the same location, depending on whether the survey lines are perpendicular or parallel to the strike of bedding. This result is observed at scales of 50 cm up to 100 m. Analysis of inverted ground resistivity together with bedding variance indicates a strong correlation between structurally complex zones with high bedding variance and a decrease in resistivity at depths greater than 25 m. This suggests that in places where the presence of disseminated gold cannot be directly detected, or where the outcrop exposure is limited due to overburden cover, geophysical data may still succeed in identifying structural complexity zones that could potentially host mineralization.

Applications of the spectral induced polarization (SIP) method to mineral exploration are limited by our knowledge of the relationships among rock texture, mineral composition, and electrical properties. Laboratory SIP responses were measured on rock samples from the Canadian Malartic gold deposit. Field SIP responses were also measured at the outcrop scale, along a profile that intersects a well-studied mineralized zone. The mineralogy and the texture of sedimentary rocks from this deposit were quantitatively determined with mineral liberation analysis. A systematic decrease (Pearson [Formula: see text]) in total chargeability with increasing fraction of the sulfide mineral interfaces associated with feldspar minerals (namely, K-feldspar and albite) was observed. On the other hand, total chargeability increased with the fraction of sulfide mineral interfaces associated with carbonates and micas (Pearson [Formula: see text]). At Canadian Malartic, proximal alteration in the mineralized zones is marked by rocks that lack a foliation plane and that were subjected to pervasive K-feldspar, albite, and pyrite alteration. In contrast, distal alteration in sedimentary rocks is marked by biotite, albite, carbonate, and pyrite that are oriented along the regional [Formula: see text] foliation. In the least-altered (LA) sedimentary rocks, quartz and biotite are associated with pyrrhotite and ilmenite as the main sulfide and oxide mineral phases, respectively. SIP measurements conducted at district and outcrop scales and along a drill core indicated that proximally altered sedimentary rocks were characterized by low total chargeability values ([Formula: see text] to [Formula: see text] in the laboratory and [Formula: see text] in the field). In contrast, the LA sedimentary rocks were characterized by total chargeability values up to [Formula: see text] in the laboratory and [Formula: see text] in the field. We conclude that mineralized zones associated with this type of ore deposit are characterized by low chargeability anomalies.

The Canadian Malartic deposit is a world class intrusion-related Archean gold deposit hosted in the Pontiac Subprovince, Superior Province, in Quebec, Canada. Laboratory petrophysical properties measurements were performed on 824 rock samples collected from the various rock types observed within the ore body and peripheral host rocks. The various rock types present in the Malartic District, mainly meta-sedimentary rocks, felsic-intermediate intrusive rocks and mafic dykes have contrasting grain densities and magnetic susceptibilities. Using support vector machines, it is shown that these two physical properties can be used to predict the rock type of a sample with an average precision and recall rate of 89%. Within the meta-sedimentary rocks class, variations in magnetic susceptibility are due to the changes in mineralogy associated with hydrothermal alteration. These are caused by the destruction of iron-bearing silicate minerals and magnetite in unaltered rocks (10−4 to 10−3 SI) to form pyrite, carbonates, K-feldspar and Fe-depleted hydrothermal biotite in altered rocks (10−5 to 10−4 SI). Within the felsic-intermediate intrusive rocks, grain densities below 2.7 g/cm3 and magnetic susceptibilities in the 10−6 to 10−4 SI range yield the highest probabilities that a rock has been submitted to carbonate and pyrite alteration. However, magnetic susceptibility and grain density of these intrusive rocks are also dependent on their Fe2O3/Al2O3 and TiO2/Al2O3 ratios, which are not related to the hydrothermal alteration footprint, but are rather due to distinct protolith compositions. Within the mafic dykes, grain density is the best indicator of hydrothermal alteration. Unaltered mafic dykes (2.95 to 3.10 g/cm3) are mostly composed of amphibole, whereas altered mafic dykes (2.70 to 2.95 g/cm3) have reduced amphibole contents and higher abundances of carbonates, pyrite, quartz, and biotite alteration. The support vector machine classifier is extended to predict if meta-sedimentary rocks, felsic-intermediate intrusive rocks and mafic dykes have undergone hydrothermal alteration with average F1 scores of 73%, 69% and 93%, respectively. In altered meta-sedimentary rocks, the integration of grain density and magnetic susceptibility allows the identification of altered but unmineralized samples. The classifier is further extended to predict if the gold content of meta-sedimentary rocks is above or below threshold values of 0.01, 0.1 and 1 ppm with average F1 scores of 83%, 80%, and 76%, respectively. Using conceptual models of the rock physical properties at the Malartic District scale, it is shown that ground magnetic surveys are the most promising geophysical tool for early-stage greenfield exploration of this type of deposit. However, depending on the scale at which the surveys are conducted, magnetic susceptibility contrasts between the various investigated rock types can overshadow the specific signatures of hydrothermally altered rocks. This in part explains why past airborne geophysical exploration campaigns for this type of deposit in the Malartic District were inconclusive. Finally, the machine learning process used in this case study can be applied in advanced exploration stages, during which drilling and either subsequent laboratory petrophysical analyses of core samples or downhole geophysical surveys produce large amounts of data that can be used to train prediction models.

Spectral induced polarization (SIP) measurements are now widely used to infer mineralogical or hydrogeological properties from the low-frequency electrical properties of the subsurface in both mineral exploration and environmental sciences. We present an open-source program that performs fast multi-model inversion of laboratory complex resistivity measurements using Markov-chain Monte Carlo simulation. Using this stochastic method, SIP parameters and their uncertainties may be obtained from the Cole-Cole and Dias models, or from the Debye and Warburg decomposition approaches. The program is tested on synthetic and laboratory data to show that the posterior distribution of a multiple Cole-Cole model is multimodal in particular cases. The Warburg and Debye decomposition approaches yield unique solutions in all cases. It is shown that an adaptive Metropolis algorithm performs faster and is less dependent on the initial parameter values than the Metropolis-Hastings step method when inverting SIP data through the decomposition schemes. There are no advantages in using an adaptive step method for well-defined Cole-Cole inversion. Finally, the influence of measurement noise on the recovered relaxation time distribution is explored. We provide the geophysics community with a open-source platform that can serve as a base for further developments in stochastic SIP data inversion and that may be used to perform parameter analysis with various SIP models.

The structural setting of the Pontiac Subprovince in the vicinity of the world-class Canadian Malartic gold deposit has been revisited by combining and reinterpreting airborne geophysical surveys together with a century of structural observations. Felsic-intermediate intrusive bodies are a key component of this deposit. Defining the regional and local favorable structural setting for intrusive rock emplacement within the clastic meta-sedimentary rocks of the Pontiac Group, south of the Cadillac Larder Lake Deformation Zone, may further assist gold exploration in similar tectonic environment. Three structural domains are interpreted in the area based on the geometry of the bedding, folds and structural fabrics related to the three major phases of deformation. During these events four phases of magmatism and one metamorphic episode occurred. The North domain, which hosts the Canadian Malartic deposit, is characterized by highly variable bedding orientations produced by the interference patterns of isoclinal F1 folds overprinted by open to tight, steeply dipping, F2 folds. The bedding in this domain is cut by a penetrative S2 biotite foliation, which possibly built on rheological changes enhanced by metasomatism in the footprint of the Canadian Malartic deposit. By contrast, the Central and South domains display homogeneous bedding orientations cross-cut by a discrete S2 biotite foliation and syn- to late- D2 metamorphic porphyroblasts. In all domains the L2 stretching lineation consistently plunges at about 60 degrees toward the east. D3 is a minor deformation event in the Pontiac Subprovince, which possibly correlates with late dextral transcurrent movement along the Cadillac Larder Lake Deformation Zone. Intrusive bodies were emplaced throughout the first and second deformation events. Phase 1 monzonite, quartz-monzodiorite and granodiorite plutons (ca. 2683–2680 Ma) intruded into consolidated Pontiac sedimentary rocks during D1. Phase 2 quartz-monzodiorite bodies (ca. 2679–2676 Ma) predominantly formed in the North Domain along F1 fold axial surfaces prior to or at the onset of D2. Phase 3 basic dykes (ca. 2675–2673 Ma) cross-cut earlier felsic-intermediate intrusions across all domains and subsequently developed an S2 metamorphic hornblende foliation. Phase 4 magmatism (ca. 2672–2662 Ma) is related to the Decelles Batholith S-type granite and pegmatite, which is interpreted to be contemporaneous with the peak of regional metamorphism. The Decelles Batholith may project at depth underneath the Canadian Malartic deposit and could be associated with magmatic-hydrothermal mineralizing fluids in the Canadian Malartic deposit. Gold mineralization at Canadian Malartic is spatially located on the contact of Phase 2 quartz-monzodiorite bodies. The proximity to the Cadillac Larder Lake Deformation Zone combined with the rheological contrast between steeply dipping clastic meta-sedimentary rocks and quartz-monzodiorite intrusions favored the protracted failure of the contacts between these two rock masses, thus forming favorable conduits for hydrothermal fluids (e.g., the Sladen Fault Zone). Such specific rheological behavior is demonstrated by domains of structural complexity, emphasized by the variance of the bedding dip. These structurally complex zones systematically host gold mineralization in the Canadian Malartic district of the Pontiac Subprovince proximal to the Cadillac Larder Lake Deformation Zone.

Alkali-silica reaction (ASR) is one of the most damaging problems with chemical alteration causing mechanical and physical changes in the concrete properties. Therefore in order to control concrete durability, a fast, non-destructive method aiming at detecting ASR is needed before an advanced stage (external visible signs) especially in field investigation. This paper focuses on the validation of the complex electrical measurements using spectral induced polarization (SIP) to monitor ASR development in concrete. Significant relationships were studied between the electro-chemical and physico-mechanical properties of concrete affected by ASR over time (22 weeks in accelerated conditions). Non-destructive methods (SIP, ultrasonic pulse velocity, expansion, and mass variation) and destructive testing (compressive strength, modulus of elasticity, petrographic examination) were used to monitor concrete mixtures (reactive and reference) submitted to accelerated laboratory test conditions. A significant evolution of SIP responses (phase lag, total chargeability, mean relaxation time) was found in the concrete specimens affected by ASR over the low frequency range. This observation was confirmed by other measurements specially expansion, UPV, modulus of elasticity, and petrographic analysis. Finally, this paper suggests SIP method as an effective non-destructive technique to monitor ASR development.

The distress caused by alkali-silica reaction (ASR) to concrete structures can occur long time after the reaction has been triggered. We propose to use spectral induced polarization (SIP) as a non-destructive method for early detection of ASR reaction before the damage is apparent. Our research deals with monitoring the difference between the electrical behaviour (phase lag, bulk resistivity, relaxation time, total chargeability) of non-reactive (NR) and reactive (RC) concrete samples affected by ASR. Laboratory measurements of complex resistivity were done in the frequency range 1.43 mHz-20 kHz. A Debye decomposition (DD) model was used to determine the DC bulk resistivity ( ρ 0 ), and the distribution of the chargeability as a function of the relaxation time. The total chargeability ( Σ m ) and the mean relaxation time ( τ mean ) were computed as DD parameters to characterize the ASR development. While ρ 0 was strongly dependent on the electrical conductivity of the solution (water, NaOH) with a constant change with time, both concrete mixtures saturated with the same solution showed an increase ρ 0 due to the development of the matrix skeleton over time. In this work, ρ 0 did not play a significant role to characterize the ASR development. For the samples affected by ASR, a significant evolution of Σ m and τ mean was found compared to the non-affected mixtures (RC, NRC-H2O), specially over the low frequency range [1.48 mHz τ > 10 2 s that could be associated with the polarization of aggregates affected by ASR.

Canadian Malartic is a large-tonnage, low-grade Archean gold deposit (16.3 Moz, 1.08 g/t Au) located in the Abitibi region of Quebec, Canada. A large part of the mineralization is hosted in the Pontiac Group metasedimentary rocks, which consist of mudstones to greywackes at upper greenschist to amphibolite facies. In exploration and production environments, these lithologies are challenging to characterize by conventional core logging methods, while in a research setting thin section-sized samples (2 cm × 4 cm) do not capture the full extent of mineralogical variability, which can extend from centimeters to meters away from mineralized zones. Here, high-resolution hyperspectral imagery (0.2–1.0 mm/pixel) in both shortwave infrared (SWIR, 1000–2500 nm) and longwave infrared (LWIR, 8000–12000 nm) is acquired for over two thousand meters of drill core, and is used to visualize changes in mineralogy and mineral chemistry related to metamorphism and hydrothermal alteration. Unaltered metasedimentary rocks contain metamorphic white mica with AlVI contents varying between 1.90 and 1.75 apfu (2195 to 2203 nm), depending on metamorphic grade. Hydrothermal alteration is characterized by white mica which becomes progressively more phengitic with increasing alteration intensity, with AlVI contents ranging from 1.70 to 1.50 apfu (2204 to 2212 nm). Phengitic white mica extends from meters to tens of meters away from major mineralized zones, and can be used as a vector towards mineralization in an exploration setting. White mica composition is correlated to Au content, and can be used to discriminate between unmineralized ( 0.3 g/t Au, 2205–2208 nm), and highly mineralized (>1.0 g/t Au, >2208 nm) samples, which is a simple metric that can directly be applied to sort ore in a production environment. The Mg# (molar Mg/[Mg + Fe]) of biotite, on the other hand, is unaffected by metamorphic grade, and consistently is Mg# 55–60 (2251–2250 nm) in unaltered samples. In mineralized samples, biotite is Mg-rich (Mg# >65