Speaker: Sylvain Favier

Date: Thursday 20th of May 2021, 1:20 pm.

Abstract:

La formation des gisements de nickel latéritique en Nouvelle-Calédonie est classiquement supposée être régie par des processus supergènes, ainsi que par l’écoulement gravitaire d’eau de pluie favorisant l’enrichissement en nickel à la base du profil latéritique. Cependant, l'hétérogénéité de la distribution du nickel semble avoir été favorisée par des processus secondaires contrôlés par les effets combinés de la tectonique, de l'évolution géomorphologique et des systèmes hydrologiques depuis la formation de la latérite primaire. Les processus de circulation des fluides et de transfert de masse ne sont pas purement descendants dans des conditions de basse température, mais peuvent également être liés à des circulations latérales de fluides, et à un drainage local le long des zones endommagées à proximité des failles (Cathelineau et al., 2016a ; 2016b ; Myagkiy et al., 2019). Cette étude vise à étudier par la modélisation du transport réactif l'impact de fractures discrètes sur la distribution du nickel.

La dissolution d’une colonne d’olivine fracturée est ainsi simulée. Cette dissolution est supposée être cinétiquement contrôlée tandis que la précipitation des produits d'altération secondaires est considérée comme se produisant selon un équilibre local. Les résultats de deux approches numériques différentes sont présentés et discutés. La première est basée sur un modèle 1D à double porosité d'une colonne d'olivine, utilisant PhreeqC associé à la base de données thermodynamiques llnl.dat. La seconde est une modélisation 2D du transport réactif dans un milieu poreux fracturé basée sur le couplage de PhreeqC et Comsol Multiphysics par ICP.

Alors que le modèle 1D vise à décrire la tendance générale de la progression du front d'altération et la redistribution globale des minéraux, le modèle 2D se concentre sur la géométrie des fractures pour mettre l’impact d’un réseau de fracturation sur la redistribution du nickel.

Les résultats des deux modélisations sont ensuite comparés aux données de terrain observées en Nouvelle-Calédonie et à une modélisation précédente d'un profil d'olivine sans fractures (Myagkiy et al., 2017) pour valider les modèles et mettre en évidence les différences induites par le réseau de fractures.

 

Speaker: Augustin Gouy

Date: Thursday 29th of April 2021, 1:30 pm.

Abstract:

De nos jours, la demande croissante en énergie, ainsi que les préoccupations grandissantes concernant le changement climatique, ont conduit à la recherche de sources d’énergie alternatives. Le gaz de houille, ou méthane de houille (“Coal Bed Methane” en anglais), peut contribuer à cette diversification et à un futur mix énergétique plus respectueux de l'environnement. Actuellement, le CBM est une source d’énergie précieuse aux États-Unis ou en Chine. Les réserves significatives de gaz estimées en Lorraine (3,71x10 11 m 3 contre 2x10 13 m 3 pour la totalité des ressources aux USA) représentent une cible prometteuse pour l'exploitation du CBM en France. En effet, sur la base de la réutilisation des mines de charbon inondées abandonnées, le méthane peut être récupéré par simple pompage de l’eau saturant le massif. La récupération du CBM peut également être améliorée par l’injection de CO2 . Dans ce qu'on appelle la récupération améliorée du méthane de houille (ECBM), le méthane est remplacé par le CO2 qui sera ainsi stocké dans la roche de manière permanente par le biais d'un processus préférentiel d'adsorption / désorption in situ. Outre l’augmentation attendue de la production de méthane, cette technique présente l’avantage supplémentaire de permettre le stockage de grandes quantités de CO2 et de contribuer ainsi à la réduction des émissions de gaz à effet de serre. L'objectif du projet est de mettre en œuvre des simulations numériques reproduisant les mécanismes mis en jeu lors de la récupération du CBM dans le charbon. Le charbon est un milieu à double porosité caractérisé par un réseau de fractures orthogonales, connues sous le nom de cleats (dont l'ouverture est inférieure à 0,1 mm) à travers lesquelles l'eau et le gaz percolent, entourant une matrice de charbon où le méthane est principalement piégé sous la forme d'une couche adsorbée dans les micropores où la diffusion prévaut. Les processus d’adsorption / désorption entraînent des déformations de la matrice rocheuse qui peut gonfler ou se contracter. L’efficacité de la récupération du méthane est ainsi déterminée par le couplage entre la dynamique d’écoulement diphasique, les mécanismes de sorption / désorption et les processus de déformation résultant des variations de la contrainte effective. Le projet se concentre sur les phénomènes en jeu dans les fractures de charbon.

Speaker: Franck Amoih

Date: Thursday 16th of April 2021, 1:00 pm.

Abstract:

Les gisements de charbon lorrain renferment une ressource en gaz non conventionnel (CH4 majoritairement) estimée à 370 milliards de m3. Ils constituent de ce fait un atout important dans la politique d’indépendance énergétique de la France et de l’Europe à l’horizon 2020. En outre, le stockage de CO2 dans les veines de charbon en vue de réduire l’émission des gaz à effet de serre, est présenté comme l’une des solutions les plus prometteuses du fait de la grande affinité du CO2 avec le charbon.

Le présent sujet de thèse inscrit dans le cadre du projet REGALOR (REssources GAzières de LORraine), vise à améliorer la compréhension des échanges gazeux sur les sites d’adsorption du charbon lorrain dans un contexte d’injection de CO2 et de récupération assistée de CH4. Il vise aussi à mieux caractériser l’influence de l’eau sur le phénomène d’adsorption. Il s’agit d’un sujet qui vient compléter et améliorer les anciens travaux de recherche réalisés dans le cadre des projets ANR CHARCO et VELCO qui ont été menés en condition statique et sur du gaz pur. Ces travaux ont montré une influence négative de l’humidité sur la capacité d’adsorption du CO2 sur le charbon à partir d’un nombre d’essais réduit. Les expériences de cette thèse seront réalisées en tenant compte des paramètres pression, température, contraintes in situ des veines de charbon. Des essais mécaniques en cellules triaxiales seront réalisés sur du charbon recompacté en vue de s’affranchir des soucis d’hétérogénéité liés au charbon ainsi que sur des échantillons naturellement fracturés. Le charbon lorrain sera caractérisé au moyen des méthodes de pétrophysique classique ; et d’autres paramètres tels que la fracturation, la nature et la maturité seront aussi déterminés. L’idée est de trouver les éventuels liens entre la capacité d’adsorption, l’humidité et les propriétés physiques du charbon lorrain.

Speaker: Christophe Antoine

Date: Thursday 8th of April 2021, 1:20 pm.

Speaker: Paul Baville

Date: Thursday 1st of April 2021, 1:20 pm.

Abstract:

To evaluate the difference between a three-dimensional Bézier triangle cubic interpolation and a parametric three-dimensional surface, I compute the absolute volume in between. The Bézier triangle is defined within the Bézier space (u,v,w are the barycentric coordinates) and the parametric surface is defined in the three-dimensional space (x,y are the Cartesian coordinates).
The Bézier triangle is defined by three points in between a three-dimensional cubic interpolation is computed. This triangle is meshed in 100 sub-triangles (each edge is divided into 10 parts) defined by 66 nodes. These nodes have barycentric coordinates which have a Cartesian equivalent and the absolute volume is computed by multiplying the vertical difference between nodes and the parametric surface and the planar surface of the Bézier triangle.
The goal of this seminar is to compare two ways to compute the vertical difference (at nodes or at sub-triangle barycenters) and to discuss another way to obtain similar results by limiting the computational time.

Speaker: Marius Huber

Date: Thursday 18th of March 2021, 1:20 pm.

Abstract:

Deep seated gravitational slope deformations (DSGSDs) are a commonly observed type of slope instability process found in mountain ranges all over the world. Their rates of movement differ fundamentally from more superficial rapid mass movements known in high relief terrain, with displacements characterized by long periods of dubious activity or time spans of inactivity. If some DSGSDs do not induce catastrophic failures but instead cease activity, others might be in a preparatory stage before a large volume catastrophic landslide. This last evolutionary stage could be favored by geological and fracturing conditions, as well as topographical, climatic, or seismo-tectonic conditions. However, due to the scarcity of very large rockslides, this possible catastrophic evolution and the reasons for it are still poorly documented and understood.

In order to explore the conditions that lead or not to an unfavorable catastrophic failure, we propose to use discrete element models to investigate DSGSDs from a mechanical viewpoint. After validating our modeling approach based on comparisons with Limit Equilibrium Analysis of simple slope geometries, we explore preconditions and mechanisms acting towards failure considering the influence of mechanical rock properties, slope topography, and pre-existing structures. Furthermore, we can implement our method on real existing case studies by shaping our simulations with digital terrain models (DTMs).

Comparison of our modeling results with those DSGSDs for which long-term and recent activity has been documented (e.g. La Clapière case in the South Western Alps of France) provides insight into their prevailing patterns of deformation (continuous or stepwise over time) and the conditions (geological, fracturing or topographic) for a catastrophic landslide or progressive stabilization. This may lead to better understanding and prediction of behavior for large scale deep seated deformation in mountainous areas, especially relevant for estimating natural hazards.


Speaker: Mustapha Zakari

Date: Thursday 25th of February 2021, 1:20 pm.

Abstract:

Unstructured grids can offer flexible and accurate geometrical approximations of reservoir models. Upscaling methods like Darcy flow-based numerical approaches can produce realistic and accurate reservoir grids. Initial flow-based upscaling approaches were developed for cartesian grids. Flow simulations were computed on groups of fine-scale grid elements covering one or many coarser elements (local, extended local, global methods). The main difficulty in flow-based methods is to define accurate fine-scale boundary conditions allowing to recover realistic coarse-scale fluxes. While cartesian grids use orthogonal directions to define vertical or horizontal pressure gradients across the simulation domain in unstructured meshes no more vertical or horizontal coarse grid alignment can be used. We investigate here fluid flow immersed boundary methods to allow incorporating cartesian-like vertical and horizontal boundary conditions for Cartesian like upscaling on unstructured meshes. I will briefly present cartesian-like flow-based upscaling approaches and then give a short review of immersed boundary methods from cartesian grids to unstructured mesh ones.

Speaker: Kou Du

Date: Thursday 18th of February 2021, 1:20 pm.

Abstract:

The aim of this presentation is to extend recent works devoted to the study of the effect of 3D pores of concave shape embedded in isotropic matrix to the case of transversely isotropic (TI) matrix. In the first part, we'll present a brief review of numerical calculation of compliance contribution tensors (which have been introduced in this seminar last year). Then, the approximate relations for the compliance contribution tensor of pores of two reference shapes, supersphere, and axisymmetrical superspheroid, are developed on the basis of 3D Finite Element Modelling, and known exact solutions for the limiting cases of spherical pores and circular crack. In the last part application to effective elastic coefficients of transversely isotropic materials in the frame of homogenization theory is presented to illustrate the impact of concavity parameter on overall properties.

Speaker: Guillaume Caumon

Date: Thursday 11th of February 2021, 1:20 pm.

Abstract:

In this seminar, I will present some work recently done in collaboration with Equinor on the use of magnetostratigraphic data to better constrain the ages of formations in the Central Graben. For this, we used individual correlation to 5 wells to the geomagnetic polarity time scale (GPTS) using the Cupydon code (Lallier et al, EPSL 2013). As the recorded series are relatively short, we also started to use magnetic chrons to constrain log-based multi-well correlation. The ultimate goal is to use these correlation results to create a composite magnetostratigraphic series that could be matched to the GPTS. Significant uncertainty yet exists about the results, as no external dating constraints are available. I will suggest some avenues for improving the methods and discuss some sources of uncertainty in the process.

 

Speaker: Ahmad Mostafa

Date: Thursday 4th of February 2021, 1:20 pm.

Abstract:

Coal beds are dual permeability systems characterized by a porous matrix enclosed within sets of orthogonal fractures known as cleats. Production of coalbed methane (CBM) consists of desorbing methane from the low permeable coal matrix to the high permeable cleat system. Unlike in conventional reservoir exploitation, sorption mechanisms cause shrinkage and swelling of the matrix which increase the complexity of the phenomena at stake, leading to complex reservoir behaviors in terms of production. A 3D discrete element method (DEM) coupled to a pore scale finite volume method (PFVM) is used here to better understand the different mechanisms at stake. The model, implemented in the open-source software Yade Open DEM (Smilauer et al., 2015), is an offspring of the hydro-mechanical model proposed by Catalano et al. (2014). The coal matrix is treated as an assembly of bonded particles interacting one with another through elastic-brittle contact laws. The pore space is discretized into tetrahedra, generated from a regular triangulation of the particle assembly. Both Knudsen and surface diffusion, as well as sorption processes, are modeled considering the coal matrix as a microporous material. The method is hydro-mechanically coupled in the sense that changes in pore pressure produce hydrostatic forces that deform the solid skeleton, while deformation of the pore space induces pore pressure changes that promote interporal flow. In addition, sorption induced deformations are taken into account by considering an additional pressure term related to the concentration of gas within the medium (the so-called solvation pressure). In this work, we first present the model and its constitutive equations. We assess its capabilities by comparing its predictions to well established solutions describing diffusive flow in porous media as well as to classic poroelasticity concepts. In particular, we focus on the influence of sorption induced deformations on the Biot coefficient estimation. Finally, we compare the model predictions to swelling test data from the literature to illustrate its consistency.

Speaker: Enrico Scarpa

Date: Thursday 28th of January 2021, 1:20 pm.

Abstract:

Nowadays, karst aquifers are strategic groundwater resources for the provision of drinking water all around the world. One-fourth of the population worldwide drinks water stored in karst aquifers, which are extremely vulnerable to various pollution sources. However, due to the complexity of karstification patterns and geometries, it is not easy to simulate groundwater flow in karstic aquifers. Over the last years, different techniques have been proposed to model karstic conduit networks (e.g., MODFLOW-CFP, GROUNDWATER). Those techniques include methods based on purely statistical approaches, object-based methods, genetic algorithms that mimic the physico-chemical processes that bring to conduit formation. Consequently, new possibilities have been introduced to evaluate the impact of karstic network architectures on groundwater flow accurately.
This work focuses on modeling the karst network stochastically using a geostatistical technique known as Multiple-Point Statistics (MPS) to characterize the karst architecture in the Dogger aquifer at the hydrogeological experimental site (HES) of Poitiers, France. In particular, here we apply the Direct Sampling (DS), an MPS simulation method, with the recent implementations proposed from the literature. Field data include core description, hydrogeological tests, and seismic surveys. A critical review of the available data set suggests the occurrence of cave patterns and several, well-identified, karst layers. The first step consists of reviewing all the (available data) collected in the last ten years and then studying the correlation between the geological, geophysical, and hydrogeological information. In the second step, DS simulations are performed to obtain several equiprobable realizations. These different DS outcomes can be compared and can estimate a global 3D karstic framework architecture of the Dogger aquifer. The resulting representation of the geological heterogeneity can support flow modeling at the HES and explore the impact of the karstic framework on the dynamic flow. 

Speaker: Fabrice Gbewade

Date: Thursday 21st of January 2021, 1:20 pm.

Abstract:

Prévoir l’endommagement hydraulique dans les systèmes géologiques constitue un challenge majeur dans l’ingénierie de subsurface. Ces processus très couplés hydro-mécaniquement (HM) qui naissent à différentes échelles, de l’échelle de la fissure (centimètre) à celle de la fracture ou de la faille (mètre ou décimètre), constituent des problèmes majeurs en géo-ingénierie qui ne sont pas encore complètement compris. L’objectif principal de ce projet est donc d’améliorer notre compréhension des processus HM qui prennent place en subsurface pour une meilleure évaluation des risques associés et de proposer des méthodes non intrusives permettant d’améliorer l’évaluation de ces processus. La présentation portera sur une synthèse bibliographique des essais de fracturations hydrauliques avec écoute acoustiques présents  dans la littérature en fonction de paramètres hydromécaniques.

Speaker: Capucine Legentil, Paul Baville & Zoé Renat

Date: Thursday 14th of January 2021, 1:20 pm.

Abstract:

Ce séminaire sera l'occasion de faire découvrir aux 3As l'IAMG (International Association for Mathematical Geosciences) et son Student Chapter. L'IAMG est l'association scientifique la plus proche des thématiques de recherche de l'équipe, et possède un Student Chapter à Nancy, qui est géré par les doctorants de l'équipe.

Le bureau de cette année est composé des doctorants de et 3ème année (Capucine Legentil, Zoé Renat et Paul Baville).

Speaker: Xiaodong Zhang

Date: Thursday 7th of January 2021, 1:20 pm.

Abstract:

It is about a diffusive crack modeling based on the introduction of a phase field. A thermodynamically consistent framework has been primarily studied aiming at developing phase field models for elastics solids by considering the hydromechanical coupling.

Speaker: Fabrice Taty

Date: Thursday 17th of December 2020, 1:20 pm.

Abstract:

Le but de cet article était donc d’utiliser la théorie des graphes pour mieux comprendre les incertitudes liées à l’interprétation structurale. En effet, compte tenu des observations limitées et de la qualité médiocre des images sismiques, définir une architecture structurale est cependant fastidieux et est sujet à des incertitudes. Le formalisme de graphe nous permettra de trouver les associations possibles entre les différentes preuves de failles spatialement tout en utilisant des règles géologiques numériques dans le but de contraindre ces associations.