DOI: 10.5281/zenodo.13710228
title: Hydro-mechanical simulation of CO2 Injection into a faulted aquifer

description: Summary

This dataset contains the results of geomechanical simulations conducted on a faulted aquifer under conditions of CO2 injection. The primary focus of the simulations is the pressure evolution within the rock matrix and along the fault, as well as the associated changes in the mechanical state, including rock deformation and fault slip. Additionally, the simulations explore the sensitivity of fault stability under varying orientations of far-field stress.

The dataset includes raw data in VTK format, as well as an illustrative Jupyter notebook that provides a comprehensive explanation of the problem's geometry, boundary and initial conditions, and an interpretation of the observed physical phenomena. The Jupyter notebook is designed to be run both online and locally.

These simulations were performed using an open-source FEM-based geomechanical simulator. Detailed instructions for running the notebook, along with a link to the geomechanical simulator, are provided in the description below.

Contributions



Emil Gallyamov did contribute to the production of the dataset and its visualisation.

Ismaël Gomes did contribute to the development of the visualisation interface through Jupyter Notebooks and Streamlit.

Guillaume Anciaux did contribute to the development of the visualisation interface and data curation.


Data collection: period and details



From 20 April, 2024 to 30 August, 2024, the .pvt and .npy files were generated, curated, and visualisation routines were developed.


Funding sources



ENAC Interdisciplinary Cluster Grant project OSGEOCGS.


Notebook demonstration

Online

An interactive notebook showcasing visualisations of the dataset is available on RenkuLab. Clicking on this link will launch an interactive session with Jupyter notebook in a separate tab. Please permit a couple of minutes for this session to load. Later, go through this notebook and update cells by pressing Shift + Enter (normal Jupyter notebook usage).  

Running locally

Alternatively, you can launch the notebook on your computer. Download the dataset, install dependencies, and launch Jupyter notebook:

pip install -r requirements_freeze.txt

jupyter notebook

Then, open notebooks/DataVisualisation.ipynb.

Reproducing the dataset

To recreate the results found in this dataset, install the solver and go through the example at examples/injection_fault.

Data structure and information

The repository has the following structure:



data: simulation results



reservoir_vs_time: bulk, solid, cohesive and fault fields for the duration of the simulation



paraview: vtk files

*.npy: numpy arrays to store numerical data at specified locations



fault_vs_angle: only cohesive fields for 180 degrees of rotation



cohesive_0*.vtu: vtu files containing absolute values of the fields

cohesive_init_0*.vtu: all the fields in these files are frozen to the initial (at rest) system state

cohesive_parsed.pvd

cohesive_init_parsed.pvd: postprocessing of Ismaël is substracting one set of data from the other one and plots the difference - increase of slip





notebooks: data visualisation through Jupyter Notebooks



images: images used for illustration in notebooks (e.g. schema of stress rotation, model geometry, etc.)



stress_rotation: contains images with scheme of rotation for angles from 0 to 180



DataVisualisation.ipynb: main notebook with visualized DataVisualisation



library: all the scripts needed to visualize DataVisualisation