Inversion of lattice network structure subjected to carbonate mercury intrusion capillary pressure

Hamiltonian Monte Carlo posterior sampling

Research output: Contribution to journalArticle

Abstract

The wide range of pore sizes in carbonates (differences of three to five orders of magnitude) aggravates the already difficult problem of multiple minima in pore network model (PNM) parameters. In this paper, we propose a method to estimate the PNM structural parameters. We have applied the method to four cases: one synthetic case and three carbonate rock samples. The PNM structural parameter estimates were conditioned to the mercury intrusion capillary pressure (MICP) via a stochastic inversion algorithm that uses molecular dynamics combined with stochastic steps using Hamiltonian mechanics. This algorithm, Hamiltonian Monte Carlo (HMC), allows for large moves in the phase space making use of the periodicity developed due to the Hamiltonian formulation, thus large modifications in the model parameters are possible. Additionally, we introduced a new estimator of the pore-size distribution using MICP. The statistics of the stochastic inversion found the multiple local minima and the global minimum of the misfit function in a difficult case. The results for the synthetic case showed that the volume and fraction of small pores in the pore-size distribution can lead to multiple local minima. For the three carbonate rock samples, the model parameters clearly showed that the volume exponents are greater than zero, which contradicts the usual assumption and inferences made by other methods. As the pore-size distribution becomes multi-modal, the volume exponent decreases but is greater than zero. These results suggest that when microporosity is ignored the volume exponent will systematically tend to be strongly underestimated (values close to zero). Our findings demonstrate that lattice-type PNM parameter estimation is a difficult problem because of multiple minima. The HMC method was able to escape from those multiple minima, resulting in better estimates of the PNM parameters.

Original languageEnglish
Pages (from-to)73-106
Number of pages34
JournalTransport in Porous Media
Volume106
Issue number1
Early online date15 Nov 2014
DOIs
Publication statusPublished - Jan 2015

Fingerprint

Hamiltonians
Capillarity
Carbonates
Mercury
Sampling
Pore size
Rocks
Microporosity
Parameter estimation
Molecular dynamics
Mechanics
Monte Carlo methods
Statistics

Keywords

  • Hamiltonian
  • Lattice
  • Mercury intrusion capillary pressure
  • Pore network model
  • Pore-size distribution
  • Prediction relative permeabilities
  • Stochastic

ASJC Scopus subject areas

  • Chemical Engineering(all)
  • Catalysis

Cite this

@article{d642f636c18c4a488e1d057b79a8fc80,
title = "Inversion of lattice network structure subjected to carbonate mercury intrusion capillary pressure: Hamiltonian Monte Carlo posterior sampling",
abstract = "The wide range of pore sizes in carbonates (differences of three to five orders of magnitude) aggravates the already difficult problem of multiple minima in pore network model (PNM) parameters. In this paper, we propose a method to estimate the PNM structural parameters. We have applied the method to four cases: one synthetic case and three carbonate rock samples. The PNM structural parameter estimates were conditioned to the mercury intrusion capillary pressure (MICP) via a stochastic inversion algorithm that uses molecular dynamics combined with stochastic steps using Hamiltonian mechanics. This algorithm, Hamiltonian Monte Carlo (HMC), allows for large moves in the phase space making use of the periodicity developed due to the Hamiltonian formulation, thus large modifications in the model parameters are possible. Additionally, we introduced a new estimator of the pore-size distribution using MICP. The statistics of the stochastic inversion found the multiple local minima and the global minimum of the misfit function in a difficult case. The results for the synthetic case showed that the volume and fraction of small pores in the pore-size distribution can lead to multiple local minima. For the three carbonate rock samples, the model parameters clearly showed that the volume exponents are greater than zero, which contradicts the usual assumption and inferences made by other methods. As the pore-size distribution becomes multi-modal, the volume exponent decreases but is greater than zero. These results suggest that when microporosity is ignored the volume exponent will systematically tend to be strongly underestimated (values close to zero). Our findings demonstrate that lattice-type PNM parameter estimation is a difficult problem because of multiple minima. The HMC method was able to escape from those multiple minima, resulting in better estimates of the PNM parameters.",
keywords = "Hamiltonian, Lattice, Mercury intrusion capillary pressure, Pore network model, Pore-size distribution, Prediction relative permeabilities, Stochastic",
author = "Juri, {Juan Ernesto} and {Van Dijke}, {Marinus Izaak Jan} and Sorbie, {Kenneth Stuart}",
year = "2015",
month = "1",
doi = "10.1007/s11242-014-0391-7",
language = "English",
volume = "106",
pages = "73--106",
journal = "Transport in Porous Media",
issn = "0169-3913",
publisher = "Springer",
number = "1",

}

TY - JOUR

T1 - Inversion of lattice network structure subjected to carbonate mercury intrusion capillary pressure

T2 - Hamiltonian Monte Carlo posterior sampling

AU - Juri, Juan Ernesto

AU - Van Dijke, Marinus Izaak Jan

AU - Sorbie, Kenneth Stuart

PY - 2015/1

Y1 - 2015/1

N2 - The wide range of pore sizes in carbonates (differences of three to five orders of magnitude) aggravates the already difficult problem of multiple minima in pore network model (PNM) parameters. In this paper, we propose a method to estimate the PNM structural parameters. We have applied the method to four cases: one synthetic case and three carbonate rock samples. The PNM structural parameter estimates were conditioned to the mercury intrusion capillary pressure (MICP) via a stochastic inversion algorithm that uses molecular dynamics combined with stochastic steps using Hamiltonian mechanics. This algorithm, Hamiltonian Monte Carlo (HMC), allows for large moves in the phase space making use of the periodicity developed due to the Hamiltonian formulation, thus large modifications in the model parameters are possible. Additionally, we introduced a new estimator of the pore-size distribution using MICP. The statistics of the stochastic inversion found the multiple local minima and the global minimum of the misfit function in a difficult case. The results for the synthetic case showed that the volume and fraction of small pores in the pore-size distribution can lead to multiple local minima. For the three carbonate rock samples, the model parameters clearly showed that the volume exponents are greater than zero, which contradicts the usual assumption and inferences made by other methods. As the pore-size distribution becomes multi-modal, the volume exponent decreases but is greater than zero. These results suggest that when microporosity is ignored the volume exponent will systematically tend to be strongly underestimated (values close to zero). Our findings demonstrate that lattice-type PNM parameter estimation is a difficult problem because of multiple minima. The HMC method was able to escape from those multiple minima, resulting in better estimates of the PNM parameters.

AB - The wide range of pore sizes in carbonates (differences of three to five orders of magnitude) aggravates the already difficult problem of multiple minima in pore network model (PNM) parameters. In this paper, we propose a method to estimate the PNM structural parameters. We have applied the method to four cases: one synthetic case and three carbonate rock samples. The PNM structural parameter estimates were conditioned to the mercury intrusion capillary pressure (MICP) via a stochastic inversion algorithm that uses molecular dynamics combined with stochastic steps using Hamiltonian mechanics. This algorithm, Hamiltonian Monte Carlo (HMC), allows for large moves in the phase space making use of the periodicity developed due to the Hamiltonian formulation, thus large modifications in the model parameters are possible. Additionally, we introduced a new estimator of the pore-size distribution using MICP. The statistics of the stochastic inversion found the multiple local minima and the global minimum of the misfit function in a difficult case. The results for the synthetic case showed that the volume and fraction of small pores in the pore-size distribution can lead to multiple local minima. For the three carbonate rock samples, the model parameters clearly showed that the volume exponents are greater than zero, which contradicts the usual assumption and inferences made by other methods. As the pore-size distribution becomes multi-modal, the volume exponent decreases but is greater than zero. These results suggest that when microporosity is ignored the volume exponent will systematically tend to be strongly underestimated (values close to zero). Our findings demonstrate that lattice-type PNM parameter estimation is a difficult problem because of multiple minima. The HMC method was able to escape from those multiple minima, resulting in better estimates of the PNM parameters.

KW - Hamiltonian

KW - Lattice

KW - Mercury intrusion capillary pressure

KW - Pore network model

KW - Pore-size distribution

KW - Prediction relative permeabilities

KW - Stochastic

U2 - 10.1007/s11242-014-0391-7

DO - 10.1007/s11242-014-0391-7

M3 - Article

VL - 106

SP - 73

EP - 106

JO - Transport in Porous Media

JF - Transport in Porous Media

SN - 0169-3913

IS - 1

ER -