Speaker: Paul Marchal

Date: Thursday 29th of September 2022, 1:20 pm.

Abstract:
Large volume of scientific and technical data have been acquired within research & development projects carried out by Orano and the Université de Lorraine. Throughout this long-term collaboration, top-ranked articles have been published based on genuine datasets resulting from specialized analyses, in particular mineralogy, isotope geochemistry and fluid composition. Samples are stored in a well referenced lithothèque and regularly used for new batches of analysis. These results have largely contributed to the knowledge of the unconformity related deposits and are considered as key parameters for defining new projects, constraining any modelling and defining machine learning analysis strategies. In order to provide access to these numerical datasets through Otelo and to connect to the Orano standards under acQuire, a joint-project for data valorization has been set up in order to develop semi-automated methods for historical data integration and design of database templates. This project aims first to provide standards methods and templates to help new data integration into Otelo’s data repository and secondly to build tools to format and integrate historical datasets. The application of these methods and standards contribute to making datasets easier to find, accessible, interoperable and reusable (FAIR Principles). For dataset, we focused especially on drillhole data analysis, i.e. data associated with sampling on cores. Concerning new data integration, initial work was the redaction of good practice guidelines1 for data management based with focus on: Explicit & coherent data naming; access & centralization of all the necessary data; use of international standards; minimization of useless duplicates & blanks. To implement, we elaborated a standardized data template containing all required metadata data fields, in order to generalize for each new project, the realization of an integrable data file collecting all samples analysis data. Then, based on established models and templates, we designed a JavaScript web-application, on the one hand to help new file edition & also to assess historical data files quality and on the other hand extract information from sparse files. We had also implemented a geographical visualization tool based on drillholes. It allows to compare data from different projects & to extract samples lists according to available attributes, either it’s about project metadata or analysis data values.

Speaker: Luyu Wang

Date: Thursday 12th of May 2022, 1:15 pm.

Abstract:

In a geological repository, percolation of gas through clay rock may generate an expansion of the conductive pathways, which will lead to a localized displacement of water away from these paths. To this end, this work investigates the hydromechanical coupling during gas migration in saturated porous media. A mesoscopic model based on the multiple relaxation time Lattice Boltzmann Method (MRT-LBM) and the spring model is proposed, in which the interpolated bounce-back (IBB) boundary is applied to capture the exact displacement of the deformable pore walls with irregular shapes. We propose a numerical assessment of this non-conforming method, coupling IBB and LBM, by comparing with the reference solution obtained on the boundary-conforming grid. The impact of grid resolution on accuracy is discussed for single- and two-phase flow conditions.

Speaker: Nathan Amrofel

Date: Thursday 15th of September 2022, 1:15 pm.

Abstract:
A full understanding of the migration behavior of corrosion gases in clay rock is of fundamental importance for the reliability of scenarios predicting the long-term evolution of geological repositories. Due to the low permeability of host clay rock, the produced gas will accumulate as a distinct phase until the pressure becomes large enough. The high pressure generated will desaturate the surrounding clay rock by displacement of pore water far along gas paths, but also by the diffusion of water vapor through the gas. In order to better understand the impact of key transport processes occurring in gas flow in clay material, a pore-scale numerical study taking into account the capillary-dominated two-phase flow, the diffusion of water vapor in the gas phase and specific features of nanoporous materials such as kelvin effect is proposed. The work has been carried out using the Smoothed Particle Hydrodynamics (SPH) method, a Lagrangian and meshless method which has emerged as an efficient and reliable tool for simulating complex fluid flows. A drying algorithm with Kelvin effect, which drives the thermodynamic equilibrium between the fluid phase and the gas phase at nanoscale, has been implemented in a two-phase flow SPH code, initially developed at IRSN.

Speaker: Guillaume Caumon

Date: Thursday 28th of April 2022, 1:15 pm.

Abstract:

This seminar will mainly consist of a review and explanation of some aspects of a recent paper by Miltenberger, Mukerji, Hariharan, Passalacqua & Nesvold (Mathematical Geoscience 2022), entitled: "A Graph-Theoretic Monte Carlo Framework for Comparing Delta Surface Dynamics and Subsurface Structure in Numerical Models and Physical Experiments".

Speaker: Narges Dashtbesh

Date: Thursday 23rd of June 2022, 1:15 pm.

Abstract:

Understanding and modelling contaminant transport is necessary to assess the lifetime of pollution sources and their severity and optimize the remediation strategies. The transfer of contaminants from the NAPL (Non-Aqueous Phase Liquid) phase to the aquifer is a multi-scale problem driven by mass transfer between both phases and is generally described by local non-equilibrium models. While the mass exchange coefficient, an important property in such models, play a key role in the fate of the pollution source, it is usually approximated by a constant value estimated from empirical correlations. However, it generally shows a transient behaviour and can evolve with NAPL phase composition and relative solubilities, which remains poorly studied. In this work, we upscale numerically this effective property from 2D pore-scale numerical simulations. We study the impact of different factors on the form and behaviour of the mass exchange coefficient in the presence of a multi-component NAPL source depending on whether the internal mass transfer is limiting. The potential implications of replacing this time-and-space-dependent mass transfer coefficient with a constant and unique value are discussed.

 

Speaker: Xiaodong Zhang

Date: Thursday 21st of April 2022, 1:15 pm.

Abstract:

The prediction of fracture mechanisms in engineering materials and structure components is of great interest. The regularized damage modeling named as phase field approach has shown great efficiencies by approximating sharp crack surface with a diffusive variable to automatically study crack initiation and propagation without additional considerations of sharp crack discontinuities. This presentation proposes a phase field model for fracture in poroelastic media based on total energy minimization. We model the fluid flow by modifying Darcy's law considerig an additive part because of the crack and the fluid flow is controlled by the elastic energy and dissipation energy not the mass conservation balance as most papers. We use a staggered scheme and implement our approach in FEniCS.  Several benchmark examples are first realized for the assessment and validation of the proposed phase field model.

Speaker: Julien Herrero

Date: Thursday 16th of June 2022, 1:15 pm.

Abstract:

Evaluation of georesources involves to appropriately manage the level of detail needed in geomodels for subsurface porous flow and transport problems. The non-uniqueness of the problems considered in aquifer or reservoir modeling calls for using methods such as stochastic Bayesian inversion to appropriately estimate rock property parameters and quantify uncertainty. These approaches often consider a fixed number of model parameters. In geomodeling, however, some model components are discrete at the scale of concern (e.g., minerals, facies, fractures, layers), hence lead to unknown number of parameters in the inverse problem. To address this issue, we propose to use transdimensional Monte Carlo methods (also known as reversible jump Markov chain Monte Carlo) which are a way to solve the inverse problem with a suitable geological parameterization when the number of model parameters is an unknown. We consider as a first application example the 2D case of implicit layer interfaces defined by a level-set in a porous reservoir model. Using a Voronoi diagram parameterization, the number of discrete layers becomes an unknown parameter. Rather than randomly perturbing the petrophysical field of interest, we build the prior model by an averaging process of log data located inside layers. Hence, each layer is defined by thickness and interface depth information, and a constant petrophysical value. Through a history matching problem, a set of flow simulations is performed to generate production data used in the Markov chain as an acceptance criterion. These numerical simulations can only be solved on a model discretization conformal to discontinuities. To address this challenge, we capitalize on the local mesh updating strategy presented by Legentil et al. (2022). First results on a set of horizontal layers demonstrate that the algorithm is capable to capture the main geological discontinuities from a prior permeability model and flow data, suggesting that this transdimensional tool could be applied for more complex geometries such as anticlines or faulted reservoir models.

Speaker: Guillaume Modeste

Date: Thursday 14th of April 2022, 1:15 pm.

Abstract:

Malgré la forte émission de gaz à effet de serre liée à sa combustion, le charbon demeure une ressource stratégique. Produit lors du processus d'houillification, du méthane (ou gaz de charbon) est présent dans les couches de charbon, adsorbé à la surface des pores de la matrice charbonneuse. En diminuant la pression de fluide dans la couche de charbon par le pompage continu des eaux naturellement présentes, le méthane se désorbe puis est extrait. Comme toute activité souterraine, la question sur l'impact de l'extraction de gaz de charbon sur les déplacements de surface et la sismicité de la région se pose. La question est d'autant plus pertinente que la désorption induit des effets physico-mécaniques supplémentaires. La transposition des connaissances concernant les réservoirs conventionnels aux couches de charbon n'est donc pas triviale. En premier lieu, une cartographie des déplacements de surface sur la période 2014-2019 est réalisée par la technique d'interférométrie radar. La carte des résultats se révèle être le reflet des exploitations souterraines passées et présentes du bassin houiller lorrain. Le soulèvement dû à l'ennoyage des travaux miniers est par ailleurs nettement détecté. Ensuite, après une synthèse des propriétés géomécaniques des éléments géologiques et des contraintes régionales, un modèle géomécanique est élaboré pour étudier les déplacements de surface. En reproduisant la déformation volumique de la couche de charbon, la propagation de la perturbation à travers le recouvrement et ses répercussions en surface sont modélisées. La déformation induite par l'extraction est estimée à partir d'une loi empirique. Des déplacements cumulés de l'ordre de quelques centimètres peuvent être attendus. En dernier lieu, les travaux investiguent le potentiel impact de l'extraction sur la sismicité induite et la réactivation de failles. Toutefois, cette partie des travaux n'est que fraichement abordée et aucun résultat préliminaire n'est disponible. Ainsi, à travers cette étude, des informations sont apportées concernant l'état actuel des déplacements de surface dans le bassin houiller lorrain et des éléments de réponse concernant le potentiel impact de l'extraction de gaz de charbon sur les déplacements et la sismicité induite.

Speaker: Gloria Arienti

Date: Thursday 19th of May 2022, 1:15 pm.

Abstract:

Introduction to the geomodelling workflow that is being used to build the new 3D structural model of the Italian North-Western Alps. Input data are represented by structural surveys and detailed geological mapping, defining a truly 3D dataset that compensates for the absence of subsurface data thanks to important differences in elevation of 3-4 km, from valley floors to mountain summits. Our modelling workflow is based on a first step of conceptual modelling in vertical cross-sections, based on classical and sound structural concepts, followed by interpolation with implicit and explicit surface algorithms. In our area, geological complexity is given by several finite faults, internal to the model, that intersect each other as the result of multiple brittle deformation events. The basement units are also interested by polyphase ductile tectonics: isoclinal and cuspate-lobate folds, polyphasic folds, and other structures that produce great thickness variations of the lithological bodies.