In digital rock physics, surface relaxation of NMR decay on a 3D porous media is often simulated with Random-Walker methods, where Brownian walkers representing the magnetized nuclei are absorbed by the solid surface under an interfacial absorption probability. Previous simulations usually neglected the blocky pixel appearance of 3D pore space images, which results in over estimation of the sur...

Backscatter SEM imaging and small-angle neutron scattering (SANS) data are combined within a statistical framework to quantify the microstructure of a porous solid in terms of a continuous pore-size distribution spanning over five orders of magnitude of length scale, from 10 A to 500 microm. The method is demonstrated on a sample of natural sandstone and the results are tested against mercury p...

Colloidal Gas Aphron as a mobility control in enhanced oil recovery is becoming attractive; it is also designed to block porous media with micro-bubbles. In this paper, the effects of surfactant concentration, polymer concentration, temperature and salinity on the bubble size of the Colloidal Gas Aphron were studied. Effects of injection rates, Colloidal Gas Aphron fluid composition, heterogene...

We developed two new EOS additions to the TOUGH+ family of codes, the RealGasH2O and RealGas. The RealGasH2O EOS option describes the non-isothermal two-phase flow of water and a real gas mixture in a gas reservoir (including a tight/shale gas one). The gas mixture is treated as either a single-pseudocomponent having a fixed composition, or as a multicomponent system composed of up to 9 individ...

Due to the intricate structure of porous rocks, relationships between porosity or saturation and petrophysical transport properties classically used for reservoir evaluation and recovery strategies are either very complex or nonexistent. Thus, the pore network model extracted from the natural porous media is emphasized as a breakthrough to predict the fluid transport properties in the complex m...

Modeling important engineering problems related to flow-induced damage (in the context of hydraulic fracturing among others) depends critically on characterizing the interaction of porous media and interstitial fluid flow. This work presents a new formulation for incorporating the effects of pore pressure in a nonlocal representation of solid mechanics. The result is a framework for modeling fl...

A general volume-averaging technique is used to derive equations satis ed by the average scalar stresses and strains in multicomponent porous rock. The resulting equations are combined with general thought experiments to produce the e ective-stress rules that determine the volumetric changes of the rock induced by changes in the con ning and uid pressures. The composite porous material speci ca...

Microbial plugging has been proposed as an effective low cost method of improving sweep efficiency by reducing permeability of water-swept regions of a reservoir. To optimize the field implementation of this method, better mechanistic and volumetric understanding of biofilm growth within a porous medium is needed. In particular, the engineering design hinges upon a quantitative relationship bet...

Truly multidimensional methods for hyperbolic equations use flow-based information to determine the computational stencil, as opposed to applying one-dimensional methods dimension by dimension. By doing this, the numerical errors are less correlated with the underlying computational grid. This can be important for reducing bias in flow problems that are inherently unstable at simulation scale, ...

We are concerned with central differencing schemes for solving scalar hyperbolic conservation laws. We compare the Kurganov-Tadmor (KT) two-dimensional (Kurganov and Tadmor, 2000) second order semi-discrete central scheme in dimension by dimension formulation with a new two-dimensional approach introduced here and applied in numerical simulations for two-phase, two-dimensional flows in heteroge...