Using Causal Inference in Field Development Optimization: Application to Unconventional Plays

Mathematical Geosciences - In the current era of big data and machine learning, a strong focus exists on prediction and classification. In industrial applications, however, many important questions...     

Spectroscopic investigation and theoretical modeling of kaolinite-group minerals and other low-temperature phases

This article summarizes some recent results obtained on the physical properties of environmental minerals, mostly kaolinite-group minerals and Fe- and Al-(hydr)oxides occurring in lateritic soils. The defective structure of these minerals, including impurities, stacking faults and radiation-induced defects, is probed using infrared spectroscopy and electron paramagnetic resonance. Resulting information bears on models of soil formation and transformation mechanisms of minerals in low-temperature environments. We underline the increasing impact of quantum chemical modeling in this field, providing straightforward interpretations of spectroscopic signals and overcoming the limits of fingerprint approaches. Importantly, the first-principles modeling of isotopic fractionation factors provides new links between mineralogical and geochemical investigations of secondary minerals.     

Barium and strontium uptake into larval protoconchs and statoliths of the marine neogastropod Kelletia kelletii

The trace elemental composition of calcified larval hard parts may serve as useful tags of natal origin for invertebrate population studies. We examine whether the trace metal barium (Ba) deposits into the calcium carbonate matrix of molluscan larval statolith and protoconch in proportion to seawater Ba concentration at two temperatures (11.5 and 17°C). We also examine strontium (Sr) uptake as a function of temperature. Using encapsulated larvae of the marine gastropod, Kelletia kelletii, reared in the laboratory under controlled conditions, we demonstrate a significant inverse effect of temperature and a positive effect of seawater Ba/Ca ratio on Ba incorporation into larval carbonates. Ba/Ca partition coefficients (D_Ba) in protoconch were 1.13 at 11.4°C and 0.88 at 17.1°C, while D_Ba in larval statolith measured 1.58 at 11.4°C and 1.29 at 17.1°C. Strontium incorporation into statoliths is also inversely affected by temperature, but there was a significant positive effect of temperature on Sr incorporation into protoconch. These data suggest larval statoliths and protoconchs can meaningfully record variation in seawater physical and chemical properties, and, hence, have potential as natural tags of natal origin.     

New igneous zircon Pb/Pb and metamorphic Rb/Sr ages in the Yaounde Group (Cameroon, Central Africa): implications for the Central African fold belt evolution close to the Congo Craton

Three meta-igneous bodies from the Yaounde Group have been analyzed for their petrography, geochemistry, and ²⁰⁷Pb/²⁰⁶Pb zircon ages. According to their geochemical patterns, they represent meta-diorites. The meta-plutonites yielded identical zircon ages with a mean of 624?±?2?Ma interpreted as their intrusion age. This age is in agreement with previously published zircon ages of meta-diorites from the Yaounde Group. The meta-diorites derived mainly from crustal rocks with minor contribution from mantle material. The ⁸⁷Rb/⁸⁶Sr isochron ages of one meta-diorite sample and three meta-sedimentary host rocks are significantly younger than the obtained intrusion age. Therefore, they are not related to igneous processes. ⁸⁷Rb/⁸⁶Sr isochron ages differ from sample to sample (599?±?3, 572?±?4, 554?±?5, 540?±?5?Ma) yielding the oldest Neoproterozoic age (~600?Ma) for a paragneiss sample at a more northern location. The youngest Rb/Sr isochron age (~540?Ma) was obtained for a mica schist sample at a more southern location closer to the border of the Congo Craton. The ⁸⁷Rb/⁸⁶Sr whole rock-biotite ages are interpreted as cooling ages related to transpressional processes during exhumation. Therefore, several discrete metamorphic events related to the exhumation of the Yaounde Group were dated. It could be shown by Rb/Sr dating for the first time that these late tectonic processes occurred earlier at more distant northern locations of the Yaounde Group and lasted at least until early Cambrian (~540?Ma) more closely to the border of the Congo Craton.     

Gradual Deformation and Iterative Calibration of Sequential Stochastic Simulations

Based on the algorithm for gradual deformation of Gaussian stochastic models, we propose, in this paper, an extension of this method to gradually deforming realizations generated by sequential, not necessarily Gaussian, simulation. As in the Gaussian case, gradual deformation of a sequential simulation preserves spatial variability of the stochastic model and yields in general a regular objective function that can be minimized by an efficient optimization algorithm (e.g., a gradient-based algorithm). Furthermore, we discuss the local gradual deformation and the gradual deformation with respect to the structural parameters (mean, variance, and variogram range, etc.) of realizations generated by sequential simulation. Local gradual deformation may significantly improve calibration speed in the case where observations are scattered in different zones of a field. Gradual deformation with respect to structural parameters is necessary when these parameters cannot be inferred a priori and need to be determined using an inverse procedure. A synthetic example inspired from a real oil field is presented to illustrate different aspects of this approach. Results from this case study demonstrate the efficiency of the gradual deformation approach for constraining facies models generated by sequential indicator simulation. They also show the potential applicability of the proposed approach to complex real cases.     

Pan-African metamorphic evolution in the southern Yaounde Group (Oubanguide Complex,Cameroon) as revealed by EMP-monazite dating and thermobarometry of garnet metapelites

Garnet-bearing micaschists and paragneisses of the Yaounde Group in the Pan-African Central African Orogenic Belt in Cameroon underwent a polyphase structural evolution with the deformation stages D₁–D₂, D₃ and D₄. The garnet-bearing assemblages crystallized in course of the deformation stage D₁–D₂ which led to the formation of the regional main foliation S₂. In X_Ca–X_Mg coordinates one can distinguish several zonation trends in the garnet porphyroblasts. Zonation trends with increasing X_Mg and variably decreasing X_Ca signalize a garnet growth during prograde metamorphism. Intermineral microstructures provided criteria for local equilibria and a structurally controlled application of geothermobarometers based on cation exchange and net transfer reactions. The syndeformational P–T path sections calculated from cores and rims of garnets in individual samples partly overlap and align along clockwise P–T trends. The P–T evolution started at ～450 °C/7 kbar, passed high-pressure conditions at 11–12 kbar at variable temperatures (600–700 °C) and involved a marked decompression toward 6–7 kbar at high temperatures (700–750 °C). Th–U–Pb dating of metamorphic monazite by electron microprobe (EMP-CHIME method) in eight samples revealed a single period of crystallization between 613 ± 33 Ma and 586 ± 15 Ma. The EMP-monazite age populations between 613 ± 33 Ma enclosed in garnet and 605 ± 12 Ma in the matrix apparently bracket the high temperature–intermediate pressure stage at the end of the prograde P–T path. The younger monazites crystallized still at amphibolite-facies conditions during subsequent retrogression. The Pan-African overall clockwise P–T evolution in the Yaounde Group with its syndeformational high pressure stages and marked pressure variations is typical of the parts of orogens which underwent contractional crustal thickening by stacking of nappe units during continental collision and/or during subduction-related accretionary processes.     

Stable isotope geochemistry of the Archean Val-d‘Or (Canada) orogenic gold vein field

A systematic study of the auriferous quartz veins of the Val-dOr vein field, Abitibi, Quebec, Canada, demonstrates that the C, O, S isotope composition of silicate, carbonate, borate, oxide, tungstate and sulphide minerals have a range in composition comparable to that previously determined for the whole Superior Province. The oxygen isotope composition of quartz from early quartz–carbonate auriferous veins ranges from 9.4 to 14.4 whereas later quartz-tourmaline-carbonate veins have ¹⁸O_quartz values ranging from 9.2 to 13.8 . Quartz-carbonate veins have carbonate (¹⁸O: 6.9–12.5 ; ¹³C: –6.2– –1.9 ) and pyrite (³⁴S: 1.2 and 1.9 ) isotope compositions comparable to those of quartz-tourmaline-carbonate veins (¹⁸O: 7.9–11.7 ; ¹³C: –8.0 – –2.4 ; ³⁴S: 0.6–6.0 ). ¹⁸O_quartz values in quartz-tourmaline-carbonate veins have a variance comparable to analytical uncertainty at the scale of one locality, irrespective of the type of structure, the texture of the quartz or its position along strike, across strike or down-dip a vein. In contrast, the oxygen isotope composition of quartz in quartz-tourmaline-carbonate veins displays a regional distribution with higher ¹⁸O values in the south-central part of the vein field near the Cadillac Tectonic Zone, and which ¹⁸O values decrease regularly towards the north. Another zone of high ¹⁸O values in the northeast corner of the region and along the trace of the Senneville Fault is separated by a valley of lower ¹⁸O values from the higher values near the Cadillac Tectonic Zone. Oxygen isotope isopleths cut across lithological contacts and tectonic structures. This regional pattern in quartz-tourmaline-carbonate veins is interpreted to be a product of reaction with country rocks and mixing between (1) a deep-seated hydrothermal fluid of metamorphic origin with minimum ¹⁸O=8.5 , ¹³C=0.6 and ³⁴S=–0.4 , and (2) a supracrustal fluid, most likely Archean seawater with a long history of water-rock exchange and with maximum ¹⁸O=3.9 , ¹³ C=–5.6 and ³⁴S=5.0 .     

Role of microorganisms in oceanic anoxic events (OAEs)

The accumulation of organic matter (OM) during oceanic anoxic events (OAEs) that took place in the Mesozoic may reflect different biogeochemical pathways and feedbacks related to palaeoenvironmental changes. It originates from a combination of nutrient cycling, metal concentrations and position of the chemocline and/or pycnocline controlling biological production. Because ocean stagnation cannot induce euxinia, an expanded oxygen minimum zone (OMZ) associated with upwelling and improved hydrological cycles is more commonly accepted. In some depositional settings, it permitted the evolution of anoxia/dysoxia to euxinia (sulphidic waters). The latter is recorded by the presence of anoxygenic phototrophs and involved denitrification, anaerobic oxidation of methane and sulphate reduction. In order to better constrain different OAE palaeoenvironments, models are suggested here that show the evolution from main microbial communities thriving through the water column to bacterial and more complex microbial mats. Therefore, microbial processes have a key role in OM production and metabolization. They directly influence the position of the chemocline. Detailed study of the OM preserved in black shales will be key to understanding OAEs.     

Trace elements in magnetite as petrogenetic indicators

We have characterized the distribution of 25 trace elements in magnetite (Mg, Al, Si, P, Ca, Sc, Ti, V, Cr, Mn, Co, Ni, Cu, Zn, Ga, Ge, Y, Zr, Nb, Mo, Sn, Hf, Ta, W, and Pb), using laser ablation ICP-MS and electron microprobe, from a variety of magmatic and hydrothermal ore-forming environments and compared them with data from the literature. We propose a new multielement diagram, normalized to bulk continental crust, designed to emphasize the partitioning behavior of trace elements between magnetite, the melt/fluid, and co-crystallizing phases. The normalized pattern of magnetite reflects the composition of the melt/fluid, which in both magmatic and hydrothermal systems varies with temperature. Thus, it is possible to distinguish magnetite formed at different degrees of crystal fractionation in both silicate and sulfide melts. The crystallization of ilmenite or sulfide before magnetite is recorded as a marked depletion in Ti or Cu, respectively. The chemical signature of hydrothermal magnetite is distinct being depleted in elements that are relatively immobile during alteration and commonly enriched in elements that are highly incompatible into magnetite (e.g., Si and Ca). Magnetite formed from low-temperature fluids has the lowest overall abundance of trace elements due to their lower solubility. Chemical zonation of magnetite is rare but occurs in some hydrothermal deposits where laser mapping reveals oscillatory zoning, which records the changing conditions and composition of the fluid during magnetite growth. This new way of plotting all 25 trace elements on 1 diagram, normalized to bulk continental crust and elements in order of compatibility into magnetite, provides a tool to help understand the processes that control partitioning of a full suit of trace elements in magnetite and aid discrimination of magnetite formed in different environments. It has applications in both petrogenetic and provenance studies, such as in the exploration of ore deposits and in sedimentology.