Skip to content

DRP-10 Estaillades Verification Overview

This report documents a controlled DRP-10 benchmark where the current voids image-to-network workflow is compared against the published Estaillades reference values reported by Muljadi et al. (2016).

The reproducible artifact for this report is notebook notebooks/31_mwe_drp10_estaillades_raw_porosity_perm.ipynb.

Goal

The benchmark answers the following question:

Given the same DRP-10 Estaillades binary image, how closely does the current voids extraction-plus-PNM workflow reproduce the paper-reference porosity and absolute permeability?

In this documentation split, DRP-10 is classified as verification because the paper reference values were obtained from a numerical OpenFOAM workflow, not from a laboratory core-flood experiment.

Studied Notebook

  • 31_mwe_drp10_estaillades_raw_porosity_perm

Sources

Current Notebook Setup

The DRP-10 notebook currently uses:

  • full-sample analysis (500 x 500 x 500 voxels), no ROI cropping
  • RAW decoding with C-style voxel ordering (order='C')
  • void convention raw == 0
  • optional pre-trimming to axis-percolating paths (trim_nonpercolating_paths = True)
  • PoreSpy extraction with geometry_repairs = "imperial_export"
  • conductance model valvatne_blunt
  • pressure-dependent water viscosity from thermo, 298.15 K
  • reference outlet pressure 5.0 MPa
  • imposed pressure gradient 10 kPa/m

Figures

DRP-10 Estaillades slices

Representative orthogonal slices of the segmented Estaillades volume used in the benchmark.

DRP-10 directional permeability comparison

Directional permeability from the current voids run and the paper flow-axis reference (Kx, Table 2).

DRP-10 extracted-network statistics

Extracted-network diagnostics for pore/throat counts, size distribution, and topological quality indicators.

Results

The full CSV outputs are available here:

Porosity and Permeability Against Paper Reference

Quantity voids estimate Paper reference Relative error
Full-image porosity [%] 10.8180 10.8 +0.17%
Network absolute porosity [%] 10.9801 10.8 +1.67%
Permeability Kx [mD] 188.5889 172.0 +9.64%

Directional Permeability (Current Run)

Axis Permeability [mD]
Kx 188.59
Ky 20.44
Kz 100.16

Network Statistics (Current Run)

Metric Value
Np 4704
Nt 8520
Mean coordination 3.62
Max coordination 33
Mean pore diameter-equivalent [um] 45.84
Mean throat diameter-equivalent [um] 31.33
Connected components 1
Giant component fraction 1.00
Dead-end pore fraction 0.204

Interpretation

The DRP-10 case shows strong agreement in porosity (sub-2% relative difference) and moderate positive bias in the paper flow-direction permeability (+9.64% in Kx).

The directional spread (Kx >> Kz > Ky) indicates anisotropy in the extracted network. That anisotropy is expected to affect which scalar permeability is comparable to a paper-reported directional value.

For this reason, the current DRP-10 benchmark is best interpreted as numerical-reference consistency evidence under the current extraction and closure assumptions, not as universal physical calibration.

Limits Of This Verification

Important limits and assumptions:

  • single DRP-10 sample (Estaillades v2) under one current workflow setup
  • no ROI/subvolume sensitivity analysis in this notebook
  • permeability comparison uses paper directionality assumptions from Table 2
  • extraction and constitutive choices (snow2 + valvatne_blunt) can shift Kabs independently of linear solver correctness
  • agreement with one OpenFOAM-referenced paper does not imply agreement with all carbonate datasets or all simulator setups

Reproducible Artifacts

  • Notebook: notebooks/31_mwe_drp10_estaillades_raw_porosity_perm.ipynb
  • Outputs:
  • examples/data/drp-10/Estaillades_v2_estimated_properties.csv
  • examples/data/drp-10/Estaillades_v2_kabs_directional.csv
  • examples/data/drp-10/Estaillades_v2_network_stats.csv
  • examples/data/drp-10/Estaillades_v2_slices.png
  • examples/data/drp-10/Estaillades_v2_kabs_comparison.png
  • examples/data/drp-10/Estaillades_v2_network_stats.png
  • examples/data/drp-10/Estaillades_v2_network_plotly.html