Table 2.

Details of the iTOUGH2-EOS7C simulations

Simulation no.HeterogeneityEnd timeInitial conditionsBoundary conditionsMethane gas source
LocationParameter adjustment (s)Water pressureSaturationLocationRate (kg s−1)Stop time
1n/an/a10 myrHydrostaticSl = 1SteadyTight2.3 × 10−11100 kyr
2High Ctz2(1/P0)10 myrHydrostaticSl = 1SteadyTight2.3 × 10−11100 kyr
3High Ctz0.5(1/P0)10 myrHydrostaticSl = 1SteadyTight2.3 × 10−11100 kyr
4*OutliersFitted10 myrHydrostaticSl = 1SteadyTight2.3 × 10−11100 kyr
5n/an/a300 kyrOverpressureSl = 1TransientTight2.3 × 10−11100 kyr
6n/an/a100 kyrPg,tight = Pw, CambrianSg,tight = 0.01SteadyNone (methane gas present initially)
  • *For simulation 4, the van Genuchten formula was fitted to the outlier mercury porosimetry datasets by adjusting both λ and 1/P0. These fitted parameters were then applied at the locations (depths) from which the outlier samples were collected in the field.