Thomas Ramstad

Relative Permeability Calculations from Two-Phase Flow Simulations Directly on Digital Images of Porous Rocks

We present results from a systematic study of relative permeability functions derived from two-phase lattice Boltzmann (LB) simulations on X-ray microtomography pore space images of Bentheimer and Berea sandstone. The simulations mimic both unsteady- and steady-state experiments for measuring relative permeability. For steady-state flow, we reproduce drainage and imbibition relative permeability curves that are in good agreement with available experimental steady-state data. Relative permeabilities from unsteady-state displacements are derived by explicit calculations using the Johnson, Bossler and Naumann method with input from simulated production and pressure profiles. We find that the nonwetting phase relative permeability for drainage is over-predicted compared to the steady-state data. This is due to transient dynamic effects causing viscous instabilities. Thus, the calculated unsteady-state relative permeabilities for the drainage is fundamentally different from the steady-state situation where transient effects have vanished. These effects have a larger impact on the invading nonwetting fluid than the defending wetting fluid. Unsteady-state imbibition relative permeabilities are comparable to the steady-state ones. However, the appearance of a piston-like front disguises most of the displacement and data can only be determined for a restricted range of saturations. Relative permeabilities derived from unsteady-state displacements exhibit clear rate effects, and residual saturations depend strongly on the capillary number. We conclude that the LB method can provide a versatile tool to compute multiphase flow properties from pore space images and to explore the effects of imposed flow and fluid conditions on these properties. Also, dynamic effects are properly captured by the method, giving the opportunity to examine differences between steady and unsteady-state setups.

Cluster evolution in steady-state two-phase flow in porous media.

We report numerical studies of the cluster development of two-phase flow in a steady-state environment of porous media. This is done by including biperiodic boundary conditions in a two-dimensional flow simulator. Initial transients of wetting and nonwetting phases that evolve before steady state has occurred, undergo a crossover where every initial pattern is broken up. For flow dominated by capillary effects with capillary numbers in order of 10(-5), we find that around a critical saturation of nonwetting fluid the nonwetting clusters of size have a power-law distribution ns similar to s(-tau)with the exponent tau=1.92+/-0.04 for large clusters. This is a lower value than the result for ordinary percolation. We also present scaling relation and time evolution of the structure and global pressure.

Correlation length exponent in the three-dimensional fuse network

We present numerical measurements of the critical correlation length exponent nu in the three-dimensional fuse model. Using sufficiently broad threshold distributions to ensure that the system is the strong-disorder regime, we determine nu to be nu=0.83+/-0.04 based on analyzing the fluctuations of the survival probability. This value is different from that of ordinary percolation, which is 0.88.

Self-affinity in the gradient percolation problem.

We study the scaling properties of the solid-on-solid front of the infinite cluster in two-dimensional gradient percolation. We show that such an object is self-affine with a Hurst exponent equal to 23 up to a cutoff length approximately g{-4/7}, where g is the gradient. Beyond this length scale, the front position has the character of uncorrelated noise. Importantly, the self-affine behavior is robust even after removing local jumps of the front. The previously observed multiaffinity is due to the dominance of overhangs at small distances in the structure function. This is a crossover effect.

Capillary-driven instability of immiscible fluid interfaces flowing in parallel in porous media.

When immiscible wetting and nonwetting fluids move in parallel in a porous medium, an instability may occur at sufficiently high capillary numbers so that interfaces between the fluids initially held in place by the porous medium are mobilized. A boundary zone containing bubbles of both fluids evolves, which has a well-defined thickness. This zone moves at constant average speed toward the nonwetting fluid. A diffusive current of bubbles of nonwetting fluid into the wetting fluid is set up.