Skip to content


Switch branches/tags

Name already in use

A tag already exists with the provided branch name. Many Git commands accept both tag and branch names, so creating this branch may cause unexpected behavior. Are you sure you want to create this branch?

Latest commit


Git stats


Failed to load latest commit information.
Latest commit message
Commit time

Carbon Capture and Storage

License: MIT


This is my BSc thesis which I defended in 9 August 2019. The repo includes my BSc works with few additions post-BSc that I made to enhance my methodology. You are recommended to read first my presented papers to understand what is this about, in the link given at the end of this material.



This is a complete modelling workflow that integrates reservoir simulation data, rock physics, time-lapse seismic, and geomechanics of CO2 injection in carbonate rock. The result of reservoir simulation is CO2 saturation and pressure changes. These information are used in rock physics modelling to (1) characterize pore structure of the carbonate rock based on Xu and Payne (2009) and (2) compute seismic response (Vp, Vs, and density) of carbonate rock as CO2 saturation and pressure increases, based on Kuster and Toksöz (1974) using the previous result of pore structure. After seismic response is computed, synthetic seismic record is simulated (seismic modelling). However, this part is worked separately using NORSAR software, not in Python, so that is not covered in this repo. Last part is geomechanical modelling, that simulates Mohr-Coulomb stress path due to change in pressure.

Case study is used from an Indonesian CCS pilot project in Gundih. Several key points of this CCS are:

  • Injection is done in carbonate reservoir,
  • CO2 is injected into 100% water-saturated zone,
  • Depth 3 km,
  • High temperature (158oC), high pressure (273 bar),
  • 10 years of injection,
  • Rate of injection 800 ton per day.

Content of this repository

Notebooks are in the main folder:

  • 01_reservoir_simulation.ipynb: displays the result of reservoir simulation done in Petrel Eclipse, compute the reservoir pressure increase due to CO2 injection
  • 02_rock_physics_modelling.ipynb: computes Vp, Vs, and density change due to CO2 saturation and increasing pressure
  • 03_xu_payne.ipynb: Xu-Payne (2009) carbonate pore characterization
  • geomechanics.m: geomechanical modelling, built in Matlab language. Will be re-written in Python language when having leisure time :)

Folder lib includes all functions used in each of the notebooks.

Folder data includes the input data.

Next task and others in mind

  • Re-writing geomechanical modelling from MATLAB code to Python.
  • Computational modelling of increasing porosity effect to change of porosity and pore compressibility.
  • Exploration of PyFEHM, PyTOUGH, Open Porous Media applied in CO2 injection.


CoolProp: to calculate CO2 property

Read my papers

My BSc research had been presented, see the extended abstracts in my ResearchGate here:


Integration of reservoir simulation, rock physics, seismic modeling, & geomechanics for CCS monitoring. Originally my BSc thesis and will be presented in EAGE 2020 Amsterdam.








No releases published


No packages published