SPE Journal
Volume 14, Number 3, September 2009, pp. 532-542

SPE-110448-PA

Physics-Driven Interface Modeling for Drainage and Imbibition in Fractures

View full textPDF ( 1,640 KB )

DOI  More information 10.2118/110448-PA http://dx.doi.org/10.2118/110448-PA

Citation

  • Prodanovic, M. and Bryant, S.L. 2009. Physics-Driven Interface Modeling for Drainage and Imbibition in Fractures. SPE J.  14 (3): 532-542. SPE-110448-PA. doi: 10.2118/110448-PA.

Discipline Categories

  • 6.3.2 Multi-phase Flow
  • 6.1.2 Faults and Fracture Characterization

Keywords

  • drainage, imbibition, rough-wall fractures, capillarity, level set method

Summary

The geometric distribution of immiscible fluid phases in fractures is not readily accessible experimentally, so aperture-scale simulations of drainage and imbibition in realistic fractures can provide valuable insight. We implement a level set method for computing the location within a fracture of the interface between two fluids controlled by capillary forces. The movement of the interface in response to changes in capillary pressure is approximated as quasistatic displacement. Fluid interfaces are thus constant mean curvature surfaces, satisfying the Young-Laplace equation. We apply a progressive quasistatic (PQS) algorithm to determine when spontaneous pore-level events occur during fluid displacement. The algorithm captures reversible and irreversible behavior. We illustrate the approach in two types of rough-walled fractures. One type is a 3D crack between irregular, impermeable surfaces; the other type is a gap between irregular 2D and 3D grain packs. We focus on the disconnected (defending) fluid volumes and (advancing) fluid main pathway as the geometric properties of the fracture are varied, notably aperture and the number of contact points between the upper and lower fracture surfaces. The curvature of the fluid/fluid interface in the plane of the fracture is often ignored by invasion percolation simulation techniques, yet it is known to influence strongly the fluid cluster properties. Our simulations establish the exact position and shape of the interface in realistic fracture geometries, from which fluid volumes, contact areas, and interface curvatures can be obtained. This establishes a new, mechanistic basis for predicting relative permeabilities in fractures and for evaluating transfer functions in dual-porosity flow models.

© 2009. Society of Petroleum Engineers

View full textPDF ( 1,640 KB )

History

  • Original manuscript received: 27 July 2007
  • Meeting paper published: 11 November 2007
  • Revised manuscript received: 18 November 2008
  • Manuscript approved: 18 November 2008
  • Published online: 27 July 2009
  • Version of record: 28 September 2009