ECMOR 2026#
Publication#
This page describes how to reproduce the results in:
Landa-Marbán, D., Sandve, T. H., and Gasda, S. E. (2026). Improving pressure communication in coarsened aquifer models for CO2 storage via explicit non-net cell treatment. ECMOR 2026.
The workflows use Bash scripts for most tasks and Python scripts for the more comprehensive processing steps.
Prerequisites#
The main preprocessing and postprocessing tools are pycopm and plopm. Install their current development versions with:
pip install git+https://github.com/cssr-tools/pycopm.git
pip install git+https://github.com/cssr-tools/plopm.git
The Troll workflow also requires OPM Flow. Some figures require PowerPoint, ParaView, or ResInsight for manual assembly or screenshots, as described in the corresponding sections.
Method: Figure 1#
Run figure1.sh from the
publications/ecmor2026/method directory:
. ./figure1.sh
Note
PowerPoint is used to assemble the subfigures and add the colored connections.
2D synthetic model#
Run figure2.sh and figure3_table1.py from the
publications/ecmor2026/results_2d_synthetic_model directory:
. ./figure2.sh
python3 figure3_table1.py
Note
PowerPoint is used to add the sensor, region, and well annotations to Figure 2.
3D synthetic model#
Run figure4_table2.py from the
publications/ecmor2026/results_3d_synthetic_model directory:
python3 figure4_table2.py
Note
The graphics in Figure 4 are generated from screenshots using ParaView.
Troll aquifer model#
The MODEL.DATA and FILLED_MODEL.DATA decks have the same dimensions and
total number of cells, 12,450,809, as the Troll aquifer model. Their remaining
properties use common homogeneous literature values, while the well locations
match those described in the ECMOR 2026 paper.
The repository does not distribute the actual Troll aquifer model. Contact the Norwegian Offshore Directorate to obtain the model, then adapt the supplied files in the results_troll/results directory.
Warning
Do not run MODEL.DATA or FILLED_MODEL.DATA before adapting them to
the actual Troll aquifer model. The supplied layouts contain many active
cells and are not substitutes for the authority-provided model.
Prepare the filled model#
Using the Petrel model provided by the Norwegian Offshore Directorate, export
the non-net cell identifiers to FLUXNUM and assign those entries a value of
zero. Then use EQUALREG to assign PORO and permeability values in the
non-net cells. See FILLED_MODEL.DATA for the implemented keyword pattern.
Run the Troll simulations#
Create the dual models and run all simulations with run_results.sh from the
publications/ecmor2026/results_troll directory:
. ./run_results.sh
Generate Table 3#
Generate Table 3 with table3.py:
python3 table3.py
Generate Figures 5 to 9#
Generate Figures 5 to 9 with figure5-9.sh:
. ./figure5-9.sh
Note
Figures 5c and 5d are generated from screenshots using
ResInsight. PowerPoint is used to add the wells
in Figure 7 and to assemble Figure 9 from figure9a.png,
figure9b.png, and figure9c.png. The color bar is taken from
to_extract_colorbar_for_figure9.png.