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 CO₂ injection in carbonate rock. The result of reservoir simulation is CO₂ 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 CO₂ 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 (158^oC), high pressure (273 bar),
• 10 years of injection,
• Rate of injection 800 ton per day.

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.

• 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 CO₂ injection.

CoolProp: to calculate CO₂ property

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

• Paper 1: Time-lapse Seismic Modelling
• Paper 2: Geomechanical Modelling