Rock Constituents Research Articles

A comparative study of the elevated radioactivity in beach placer deposits of Bangladesh.

Beach sediments are mineral deposits formed through weathering and erosion of either igneous or metamorphic rocks. Among the rock constituent minerals are some natural radionuclides that contribute to ionizing radiation exposure on Earth. Kolatoli and Kuakata are the two major beaches with heavy mineral deposits and important tourist sites in Bangladesh. Natural radioactivity in Kolatoli and Kuakata beach sand deposits along the southern coast of Bangladesh was assessed and compared to identify the sources, causes, and possible environmental impact. Result shows most of the radionuclides have higher activity concentrations than the background level, and the activity varies with the sample locations. The dominant radionuclides were found to be the radionuclides of thorium series i.e. Th-232 and Ra-228 followed by uranium series and K-40. The radioactivity in Kolatoli beach sands was observed to be much higher than Kuakata beach due to the presence of a higher content of heavy minerals i.e. illmenite, rutile, zircon, garnet and monazite. Furthermore, monazite and zircon are the two radioactive minerals that are considered to be the main contributors to the radioactivity in Kolatoli beach sand. These minerals are dominated by the activity of thorium series radionuclides i.e. Th-232 and Ra-228 surpass the activity of all other radionuclides such as U-238, U-234, Th-230, Ra-226, Po-210, and K-40. However, major contribution of radioactivity in Kuakata beach sand comes from uranium series radionuclides such as U-238, U-234, Ra-226, and Po-210. Beach morphology, sedimentological, and geochemical evolution of those minerals might be important areas of further study for the radioactivity monitoring activity in those areas.

Coupling nanoscopic tomography and micromagnetic modelling to assess the stability of geomagnetic recorders

The recording of planetary magnetic fields is often attributed to uniformly-magnetised nanoscopic iron oxides, called single-domain. Yet, the main magnetic constituents of rocks are more complex, non-uniformly magnetised grains in single or multi-vortex states. We know little about their behaviour due to limitations in defining their precise shape and internal magnetic structure. Here we combine non-destructive Ptychographic X-ray Computed Nano-tomography with micromagnetic modelling to explore the magnetic stability of remanence-bearing minerals. Applied to a microscopic rock sample, we identified hundreds of nanoscopic grains of magnetite/maghemite with diverse morphologies. Energy barrier calculations were performed for these irregularly shaped grains. For some grains, these morphological irregularities near the transition from single-domain to the single-vortex state allow for multiple domain states, some unstable and unable to record the field for significant periods. Additionally, some other grains exhibit temperature-dependent occupancy probabilities, potentially hampering experiments to recover the intensity of past magnetic fields.

Molecular insights into dual competitive modes of CH4/CO2 in shale nanocomposites: Implications for CO2 sequestration and enhanced gas recovery in deep shale reservoirs

CO2 enhanced shale gas recovery has significant potential for improving the recovery rate of adsorbed gas and facilitating the geological sequestration of CO2. Nevertheless, the feasibility and microscopic mechanisms of CO2 enhanced gas recovery in deep shale gas reservoirs remain to be elucidated. In this study, a quartz-kerogen slit model was constructed to represent the shale composite nanopores of deep shale reservoirs. Utilizing the grand canonical Monte Carlo and molecular dynamics methods, the adsorption and diffusion behavior of CH4, CO2 and their mixtures was investigated under high-temperature and high-pressure conditions representative of deep shale reservoirs. The influences of temperature, pressure and aperture size were analyzed. The microscopic mechanisms governing the CH4/CO2 competitive adsorption in deep shale nanopores were uncovered by considering the competitive interactions between different rock constituents and quantifying the various gas storage states. The results show that the adsorption capacity of CO2 is greater than that of CH4 under the deep shale reservoir conditions. The replacement efficiency of CH4 by CO2 is higher in small shale pores. High temperature has a more significant negative impact on the adsorption of CO2 compared to that of CH4. Compared to mid-deep shale reservoirs, the nanopores of deep shale reservoirs have a higher proportion of free CH4/CO2 states. Quartz exhibits greater affinity for CH4/CO2 than kerogen, but kerogen possesses a larger specific surface area, resulting in a higher adsorption capacity per unit volume. The hierarchy of CH4/CO2 diffusion capacity in various storage states is as follows: dissolved gas in the kerogen matrix < adsorbed gas on the quartz surface < adsorbed gas on the kerogen surface < free gas in the pore center. This study improves the understanding of CO2 enhanced gas recovery in deep shale gas reservoirs.

Evaluating the role of dissolved silica for dolomite formation in evaporitic environments: Insights from prolonged laboratory experiments

The mineral Dolomite CaMg(CO3)2 is a common constituent of sedimentary rocks. Despite centuries of research, the mechanism of its formation remains elusive and debated. Recent studies have shown the presence of silica in solution promote the incorporation of Mg into the carbonate mineral, forming crystal phases that may be precursors to dolomite. The goal of this study was to evaluate with laboratory experiments whether dissolved silica may play a role for dolomite formation in sabkha (i.e., salt flats) environments. Several models for dolomite formation are based on the studies of sabkhas, which are often cited as modern analogue for ancient dolomite-rich sedimentary sequences. We performed long-incubation time (i.e., up to 600 days) laboratory precipitation experiments at 30 °C with solution mimicking the sabkha pore waters (characterized by a salinity of 23 % and Mg: Ca ratio of 15) to which we added different concentrations of Si. Our results revealed a positive correlation (p-value <0.001) between Si concentration in solution and the mol% Mg of the carbonate minerals forming in the experiment. With 2 mM of Si, the bulk precipitate was comprised of 90 % stoichiometric dolomite with possible signs or ordering. Moreover, the rhombohedral morphology of the crystals is analogue to that of natural dolomite previously described from sabkha sediments. Our results suggest that dissolved Si may play an important role for dolomite formation in evaporitic environments.

The where, when, and how of ooid formation: What ooids tell us about ancient seawater chemistry

Ooids are a common constituent of carbonate rock throughout the last 3.2 billion years. The chemical and physical characteristics of ooids record essential details of the carbonate chemistry and hydrodynamic conditions of their depositional environment. As a result, ooids may be a powerful tool for reconstructing Earth's ancient climates and environments. However, lack of constraints on fundamental aspects of ooid formation, such as the role of abrasion, underscores the potential for additional discoveries utilizing ooids as paleoenvironmental proxies. Here, we use three-dimensional models built by serially grinding and imaging ooids to show how ooid shape changes during growth. We find that ooid abrasion plays only a small role in the formation of ooids, even in giant ooids, where the mass should cause abrasion to be most intense. A lack of abrasion during ooid formation suggests a greater sensitivity of ooid distribution and size to the chemistry of depositional environments. This insight enables us to harmonize our understanding of ooid formation with longstanding questions regarding their spatial and temporal distribution, as well as size variations. With a new database of ooid occurrences, we analyze the influence of seawater chemistry, carbonate shelf area, and evolutionary developments on ooid abundance in the Modern, Last Glacial Period, and throughout the Phanerozoic. Finally, we leverage our new understanding of ooid growth to show how ooid size, including giant ooids, can shed light on relative changes in alkalinity throughout the Phanerozoic.

A critique of using epitaxial criterion to discriminate between protogenetic and syngenetic mineral inclusions in diamond

Distinguishing syngenetic from protogenetic inclusions in natural diamonds is one of the most debated issues in diamond research. Were the minerals that now reside in inclusions in diamonds born before the diamond that hosts them (protogenesis)? Or did they grow simultaneously and by the same reaction (syngenesis)? Once previously published data on periclase [(Mg,Fe)O] and magnesiochromite (MgCr2O4) inclusions in diamond have been re-analysed, we show that the main arguments reported so far to support syngenesis between diamond and its mineral inclusions, definitely failed. Hence: (a) the epitaxial relationships between diamond and its mineral inclusion should no longer be used to support syngenesis, because only detecting an epitaxy does not tell us which was the nucleation substrate (there are evidences that in case of epitaxy, the inclusion acts as a nucleation substrate); (b) the morphology of the inclusion should no longer be used as well, as inclusions could be protogenetic regardless their shapes. Finally, we advance the hypothesis that the majority of inclusions in diamonds are protogenetic, e.g., they are constituent of rocks in which diamonds were formed and not products of reactions during diamond growth.

Creep Versus Consolidation in Tunnelling Through Squeezing Ground—Part A: Basic Time Effects

Although squeezing ground may undergo rapid convergences following tunnel excavation, its behaviour is often markedly time dependent due to creep or consolidation. The effects of creep and consolidation on shield tunnelling are comparatively evaluated in two companion papers, with the aim of demonstrating their qualitative similarities and distinctive features. The present, first paper investigates the basic time effects, placing focus on the time development of ground deformations and the complex interaction between ground, tunnel boring machine (TBM) and tunnel support during excavation and during construction standstills. The presented numerical simulations indicate several qualitative similarities between the two mechanisms of time dependency, in respect of the time development of ground deformations, the counter-intuitive behaviour of increasing shield loading with increasing rate of advance under certain conditions, as well as the thoroughly adverse effect of the additional time-dependent deformations taking place during construction standstills on the shield loading. However, they also underscore two prominent differences resulting from the fundamentally different nature of creep (a purely mechanical rheological process) and consolidation (a coupled hydromechanical process): first, the consistently more extensive plastic yielding in consolidating ground, which is partially associated with the seepage forces exerted by the pore water on the solid rock constituents. Second, the role of seepage forces as a potential destabilising agent, particularly for the tunnel face, which does not happen in the case of creep and may be critical for shield and cutterhead jamming. Building upon these investigations, the companion paper compares creep and consolidation with respect to the transferability of experiences about the required thrust force to tunnels of different diameter or to adjacent tunnels.

Plutonium mobility and reactivity in a heterogeneous clay rock barrier accented by synchrotron-based microscopic chemical imaging

The long-term safe disposal of radioactive waste corresponds to a challenging responsibility of present societies. Within deep geological waste disposal concepts, host rocks correspond to the ultimate safety barrier towards the environment. To assess the performance of such barriers over extended time scales, mechanistic information on the interaction between the radiotoxic, long-lived radionuclides like plutonium and the host rock is essential. Chemical imaging based on synchrotron microspectroscopic techniques was used to visualize undisturbed reactive transport patterns of Pu within pristine Opalinus Clay rock material. Pu+V is shown to be progressively reduced along its diffusion path to Pu+IV and Pu+III due to interaction with redox-active clay rock constituents. Experimental results and modeling emphasize the dominant role of electron-transfer reactions determining the mobility of Pu in reactive barrier systems. The effective migration velocity of Pu is controlled by the kinetic rates of the reduction to Pu+IV and Pu+III and the redox capacity of the involved electron donor pools. To advance our predictive capabilities further, an improved understanding of the nature and capacity of redox-active components of the reactive barrier material is fundamental. The findings represent an essential contribution to the evaluation of the long-term safety of potential nuclear waste repositories and have implications regarding the development of effective geological disposal strategies.

The Role of the Mineralogical Composition on Wettability via Flotation Test and Surface Complexation Modeling (SCM)

Minerals are the chief constituents of rocks and have varied properties, such as the surface area, surface charge, site density, etc. Hence, numerous interactions are bound to occur in a reservoir during rock–fluid (i.e., rock, crude oil and brine) interactions. This study seeks to assess the role of the mineralogical composition in the wettability of sandstone rocks (SRs) and mineral mixture (MM) using both surface complexation modeling (SCM) and a flotation test. From the considered sandstone rocks, both the experimental results and the simulated counterparts revealed that the SRs were preferentially hydrophilic. For the MM, when the mass fraction of the hydrophobic mineral was increased, the affinity of the MM became slightly hydrophobic, and vice versa. For the dominant sandstone reservoir rock minerals with predominantly negatively charged surfaces, negligible oil adsorption took place due to the interfacial repulsive forces at the oil–brine and mineral–brine interfaces. For the MM with low calcite content, the wetting preference was influenced by the mineral with a prominent surface area. Our developed model portrayed that the main mechanism of oil adhesion onto sandstone minerals was divalent cation bridging. Nonetheless, adhesion of carboxylate (&gt;COO−) onto the illite, montmorillonite and calcite sites also took place, with the latter being more pronounced.

Hyperspectral mapping of magmatic-hydrothermal sericite, Battle Mountain mining district, Nevada

We demonstrate the utility of multi-scale imaging spectroscopy data (hyperspectral data) to differentiate magmatic-hydrothermal sericite from other white micas based on grain size, in a case study of the Battle Mountain mining district in Humboldt and Lander Counties, Nevada. The Battle Mountain mining district contains porphyry, skarn, and distal disseminated deposits that formed in magmatic-hydrothermal systems, as well as sedimentary-rock hosted gold deposits of less certain origin. We generated spectral-based mineral maps of the Battle Mountain mining district to examine the spatial distribution of spectrally dominant minerals such as kaolinite, montmorillonite, and white micas (muscovite, sericite, and illite). We then used spectral-based mineral mapping to characterize the grain size and composition of white micas, enabling us to differentiate igneous muscovite from magmatic-hydrothermal sericite. By integrating our spectral-based maps with geologic mapping, airborne aeromagnetic data, and prior studies on mines and deposits in Battle Mountain, we determined that finer grained, higher Al sericite spatially correlated with the footprints of known magmatic-hydrothermal systems and is interpreted as forming from magmatic-hydrothermal alteration. Larger grained, lower Al muscovites are distal from known magmatic-hydrothermal systems and are interpreted as being igneous derived siliciclastic sediments that became constituents of the host sedimentary rocks. Techniques described in this contribution can be applied to any appropriate hyperspectral data, regardless of sensor type used to collect the data.

Empirical determination of the effective solid modulus in organic-rich shales

Calculating the change in the saturated bulk modulus of a saturated rock with new fluid properties requires a priori selection of an effective bulk modulus of the solid constituents. When the rock constituents have similar mineral moduli, the theoretical bounds on the solid modulus are close to each other. However, when solid properties vary greatly, as in organic-rich shales, the actual effective solid modulus of a physical rock may vary significantly between the bounds which results in uncertainty in the predicted change in the saturated bulk modulus of the rock. We use a semi-empirical rock physics model utilizing the Brown–Korringa equation for mineralogically heterogenous rocks and introduce three parameters to estimate the pore space compressibility, the dry frame compressibility, and the fractional position of the effective solid modulus relative to the Reuss and Voigt bounds. We optimize for these three parameters in seven organic shale formations and find that the Reuss bound for the effective solid material modulus best fits the data when organic content is high. Furthermore, we use this model to fluid substitute to 100% brine saturation and find Gassmann’s equation using the Hill average predicts similar saturated moduli to the semi-empirical Brown–Korringa rock physics model when volume fraction of solid organic matter is less than 5%. However, at higher organic contents, we find that the error using the Gassmann–Hill approach increases, and the semi-empirical Brown–Korringa model better fits the data.

Impact of Fe(II) on 99Tc diffusion behavior in illite

A comprehensive understanding of the geochemical behavior of 99Tc is of great importance for safe disposal of radioactive waste and remediation of contaminated environmental sites. Illite is one of the most common constituents of clay rocks, and thus used in this work as a model system for studying the retention and transport of 99Tc in clay-rich systems. In this study, a through-diffusion technique was applied to investigate the diffusion behavior of Tc in compacted illite clay under oxic and anoxic conditions. Particular focus of this investigation was on the role of Fe(II) on the redox state and mobility of Tc in clay. As the diffusion cells contained stainless steel filters for confining the clay plug, 99Tc-diffusion in the filters was assessed first, followed by 99Tc diffusion in the illite column with or without Fe-loading. Two types of Fe(II) loadings in illite were considered, surface complexed Fe(II) on illite edge sites and Fe(II) added as pyrite grains. COMSOL Multiphysics was used to analyze the diffusion data. The measured filter porosity was about 0.2 and the effective diffusion coefficient, De, for Tc (as TcO4−) in the filter was 0.59 × 10−10 m2/s. Tc diffusion in illite under ambient conditions showed a typical anion diffusion behavior with De in the range 0.38-0.44 × 10−10 m2/s and the anion accessible porosity ε of approximately 0.2. In the presence of Fe(II), De for Tc migration in illite was one order of magnitude lower, showing that Fe(II) has a significant effect on the migration of 99Tc. Analysis of the Tc distribution in the system suggests that most Tc was retarded at the filters, especially the ones connected to the high concentration reservoir. The remaining Tc quantities were immobilized at sample boundaries next to the filters with higher concentrations observed in the domains close to the reservoirs. Almost no Tc was immobilized in the middle part of the sample, where the Fe(II) was preloaded. This observation contradicts the anticipation, because Tc was expected to be immobilized in the middle of the sample preloaded with Fe(II). A possible explanation is that a delocalized redox reaction may occur, which means that electrons from the oxidation of the preloaded Fe(II) in the central part of the cell are transferred to the filter via the walls of the steel diffusion cell. Tc reacts with the electrons in the filter at a relative slow rate, which results in most of the Tc retarded in the filter, while some small amounts of it may diffuse through the clay.

Uranium Speciation in a Chalky Soil.

In this article, the speciation and behavior of anthropogenic metallic uranium deposited on natural soil are approached by combining EXAFS (extended X-ray absorption fine structure) and TRLFS (time-resolved laser-induced fluorescence spectroscopy). First, uranium (uranyl) speciation was determined along the vertical profile of the soil and bedrock by linear combination fitting of the EXAFS spectra. It shows that uranium migration is strongly limited by the sorption reaction onto soil and rock constituents, mainly mineral carbonates and organic matter. Second, uranium sorption isotherms were established for calcite, chalk, and chalky soil materials along with EXAFS and TRLFS analysis. The presence of at least two adsorption complexes of uranyl onto carbonate materials (calcite) could be inferred from TRLFS. The first uranyl tricarbonate complex has a liebigite-type structure and is dominant for low loads on the carbonate surface (<10 mgU/kg(rock)). The second uranyl complex is incorporated into the calcite for intermediate (∼10 to 100 mgU/kg(rock)) to high (high: >100 mgU/kg(rock)) loads. Finally, the presence of a uranium-humic substance complex in subsurface soil materials was underlined in the EXAFS analysis by the occurrence of both monodentate and bidentate carboxylate (or/and carbonate) functions and confirmed by sorption isotherms in the presence of humic acid. This observation is of particular interest since humic substances may be mobilized from soil, potentially enhancing uranium migration under colloidal form.

Bayesian inversion for rapid multiwell interpretation of well logs and core data: Applications to unconventional organic-shale formations

Multimineral analysis is widely used to calculate in situ porosity, fluid saturation, and mineralogy of rocks penetrated by a well. It delivers weight/volumetric concentrations of rock solid/fluid constituents by combining multiple borehole geophysical measurements and often is referred to as petrophysical joint inversion. Recently, a probabilistic method was developed for improved petrophysical estimation of rock constituents and their uncertainty from well logs. This method mitigates borehole and instrument-related environmental effects present in the measurements and efficiently propagates the uncertainty from measurement noise and rock-physics models (RPMs) to compositional/petrophysical estimations. The probabilistic estimation method is extended to the challenging conditions of multiple neighboring wells penetrating similar rock formations where borehole/drilling environmental conditions, borehole instruments, and measurement noise may vary from well to well. A calibration step is performed in a few key wells with core data and/or advanced borehole measurements; it enables the same RPMs and prior models to be implemented in nearby wells but with limited measurements. In addition, a precomputed surrogate model constructed with radial basis function interpolation is implemented for accurate and efficient nuclear-property calculations. The multiwell interpretation method is verified using synthetic and field examples of organic-rich shale formations. Results find that the probabilistic method (1) improves rock petrophysical/compositional estimations by mitigating borehole environmental effects and incorporating a priori knowledge, (2) yields regionally consistent compositions among wells, and (3) quantifies the uncertainty of the estimations. As a result, the probabilistic approach is especially suitable for assessing petrophysical/compositional properties in multiwell settings with complex rock constituents and/or limited borehole measurements.

Robust and interpretable mineral identification using laser-induced breakdown spectroscopy mapping

Fast and precise identification of minerals in geological samples is of paramount importance for the study of rock constituents and for technological applications in the context of mining. However, analyzing samples based only on the extrinsic properties of the minerals such as color can often be insufficient, making additional analysis crucial to improve the accuracy of the methods. In this context, Laser-induced breakdown spectroscopy mapping is an interesting technique to perform the study of the distribution of the chemical elements in sample surfaces, thus allowing deeper insights to help the process of mineral identification. In this work, we present the development and deployment of a processing pipeline and algorithm to identify spatial regions of the same mineralogical composition through chemical information in a fast and automatic way. Furthermore, by providing the necessary labels to the results on a training sample, we can turn this unsupervised methodology into a classifier that can be used to generalize and classify minerals in similar but unseen samples. The results obtained show good accuracy in reproducing the expected mineral regions and extend the interpretability of previous unsupervised methods with a visualization tool for cluster assignment, thus paving for future applications in contexts requiring high-throughput mineral identification systems, such as mining.

Molecular Simulation of Methane Adsorption in Deep Shale Nanopores: Effect of Rock Constituents and Water

The molecular models of nanopores for major rock constituents in deep shale were constructed. The microscopic adsorption behavior of methane was simulated by coupling the grand canonical Monte Carlo and Molecular Dynamics methods and the effect of rock constituents was discussed. Based on the illite and kerogen nanopore models, the discrepancies in microscopic water distribution characteristics were elucidated, the effects of water on methane adsorption and its underlying mechanisms were revealed, and the competitive adsorption characteristics between water and methane were elaborated. The results show a similar trend in the microscopic distribution of methane between different shale rock constituents. Illite and kerogen slit pores have no significant difference in methane adsorption capacity. The adsorption capacity per unit mass of kerogen is greater than that of illite due to the smaller molar mass of the kerogen skeleton and its large intermolecular porosity. Illite has a greater affinity for water than methane. With increasing water content, water molecules preferentially occupy the high-energy adsorption sites and then overspread the entire pore walls to form water adsorption layers. Methane molecules are adsorbed on the water layers, and methane adsorption has little effect on water adsorption. Kerogen is characterized as mix-wetting. Water molecules are preferentially adsorbed on polar functional groups and gather around to form water clusters. In kerogen with high water content, methane adsorption can facilitate water cluster fusion and suppress water spreading along pore walls. In addition to adsorption, some water molecules dissolve in the kerogen matrix.

Feldspar Pb isotope evidence of cryptic impact-driven hydrothermal alteration in the Paleoproterozoic

Hypervelocity impacts throughout Earth's history have profoundly affected the evolution of the continental crust. Accessory minerals like zircon are typically used to date impact events and rock-forming minerals like quartz are routinely used as shock barometers. However, feldspar group minerals – a major constituent of most crustal rocks – are generally underutilized in the documentation of impact-induced deformation and alteration. Alkali feldspar contains appreciable amounts of Pb and analysis of Pb isotopes in feldspar may offer the opportunity to identify impact-related isotopic modifications of shocked crustal target rocks and estimate their timing. Here, we apply a combination of laser ablation inductively coupled plasma (LA-ICP) and thermal ionization mass spectrometry (TIMS) Pb isotope analysis with imaging techniques, including electron backscatter diffraction (EBSD), cathodoluminescence (CL), and time of flight secondary ion mass spectrometry (ToF-SIMS), to shocked alkali feldspar from monzogranite in the oldest confirmed terrestrial impact structure (2229 ± 5 Ma) at Yarrabubba, Western Australia. Alkali feldspar preserves microstructures such as sub-planar and irregular fractures, sets of planar deformation bands that accommodate misorientations of up to ∼20°, sets of damage lamellae, and broad domains of lattice damage that can be linked to impact-related deformation. The Pb isotope compositions in alkali feldspar correlate with variations in electron diffraction band contrast – a proxy for crystallinity – and also the degree of misorientation and CL response. Less damaged alkali feldspar yields Pb model ages similar to the igneous zircon U–Pb crystallization age of the host monzogranite (∼2650 Ma), whereas younger Pb model ages correspond to zones of damage (high relative misorientation, low crystallinity, weak CL response). The observed Pb isotope behaviour implies radiogenic ingrowth of Pb, from decay of U and Th within damaged alkali feldspar, and therefore mixing with a grain-scale Pb reservoir that formed at the time of impact. The U and Th zonation in some shock-deformed alkali feldspar is broadly similar and follows the orientation of sub-planar fractures and damage lamellae. Detailed imaging reveals the zones of U and Th enrichment are associated with trains of monazite micro-inclusions, in conjunction with magnetite and/or hematite in places, which are inferred to have precipitated during impact-induced hydrothermal circulation. Hence, the Pb isotopic data record grain-scale hydrothermal alteration in superficially weakly altered monzogranite target rocks.

Unraveling Textural and Chemical Features in Volcanic Rocks Through Advanced Image Processing: A Case Study From the 2019 Paroxysmal Eruptions of Stromboli

AbstractThe Quantitative X‐Ray Map Analyzer software, a new tool for image processing, has been tested on intertwined pumices and scoriae emitted during the two paroxysmal eruptions of Stromboli of 2019, whose textural and compositional heterogeneities reflect the coexistence of low porphyritic shoshonitic‐basalts and high porphyritic shoshonites. The procedure applied was designed to quantitatively document the complex variations in texture and composition of these products, allowing substantial time reduction of analytical and data processing. The procedure utilizes Principal Components Analysis and supervised Maximum Likelihood Classification for multivariate statistical data handling of an array of X‐ray elemental maps acquired at the millimeter scale in thin‐sections. This technique permits the production of high‐contrast colored images, which allow the classification of rock constituents, extrapolating the associated modal abundances and imaging chemical variations within the glass. Results highlight the close interconnection at the microscale of two types of magma in the erupted products, manifested in each processed image by the presence of contiguous areas preserving textural bulk properties typical of the pumice or scoria. The proportion of the two magmas feeding the eruptions is not simply represented by the proportions of scoria and pumice in individual clasts, as both scoria and pumice contain glass with the composition of shoshonitic‐basalts and shoshonites. This method also allows the recognition of important discordances between the textural and chemical features of the two fractions involved, as well as discernment of the compositions of the two interacting magmas at the microscopic scale, even in those micro‐domains showing evidence of intense interaction processes.

Interpretation of well logs and core data via Bayesian inversion

Estimating in situ petrophysical and compositional properties of rocks (e.g., porosity, mineralogy, and fluid saturation) from well logs and core measurements is critical for the evaluation of subsurface fluid resources. Traditional multimineral analysis of well logs is susceptible to abnormal borehole and geometric conditions, such as shoulder-bed effects and tool/borehole-related biases. Moreover, traditional multimineral analysis rarely includes uncertainty quantification, much less assessments of the impact of measurement noise, borehole conditions, and/or inaccurate rock-physics models (RPMs) on estimated properties. Probabilistic formation evaluation with credible-interval estimations can improve the outcome of mineral and fluid evaluations from well logs by explicitly incorporating the effects of measurement biases and a priori compositional and petrophysical constraints in the estimation. A probabilistic workflow is developed to estimate petrophysical/compositional properties from well logs and available core data. The workflow consists of two sequential steps: borehole measurements (e.g., density, resistivity, and gamma ray) are first converted into a layer-by-layer earth model, which consists of piecewise constant property values, with associated uncertainty estimations. The estimated earth model is tool/borehole independent, allowing for a common baseline to effectively compare properties among neighboring wells. Next, the estimated earth-model physical properties are used to derive petrophysical properties (e.g., volumetric concentration of fluid and solid rock constituents) via petrophysical joint inversion. In each step, Bayesian inversion is implemented with Markov chain Monte Carlo sampling in contrast to deterministic algorithms used by traditional multimineral methods. The developed method has the following advantages over traditional approaches: (1) it propagates uncertainty from well logs to petrophysical properties by assimilating measurement noise, tool/borehole-related biases, and RPM errors and (2) it enables the implementation of a priori knowledge per rock class to constrain the petrophysical estimation, effectively integrating field-specific core/regional geologic data. Examples of thinly laminated and organic shale formations successfully verify the reliability, efficiency, and accuracy of our estimation method.

ОСОБЕННОСТИ ИЗОТОПНОГО СОСТАВА УГЛЕРОДА ОРГАНИЧЕСКОГО ВЕЩЕСТВА И УГЛЕРОДА И КИСЛОРОДА КАРБОНАТНОЙ СОСТАВЛЯЮЩЕЙ НЕФТЕГАЗОПРОИЗВОДЯЩИХ ОТЛОЖЕНИЙ КЕМБРИЯ СИБИРСКОЙ ПЛАТФОРМЫ

The paper presents data on the carbon (δ13С) and oxygen (δ18О) isotopic composition of the carbonate rock constituent (CRC) and carbon (δ13C) of dispersed organic matter (SOM) of rocks of the Shumnaya, Inikan and Kuonamka Formations of the Lower-Middle Cambrian of the Siberian Platform. Dependence of the carbon and oxygen isotopic composition of CRC on the facies subtype of domanicoid deposits of the Cambrian was revealed. The interrelation of the carbon isotope composition of DOM with depositional conditions of sediments of the Lower and Middle Cambrian was considered. In the south-east of the Siberian Platform, Lower Cambrian carbonate deposits with carbon anomalous isotopic composition of CRC (δ13Сср = 6,6 ‰) and carbon of DOM (δ13С = –36,5 ‰) were found.