Dynamic coupling of pore‐scale and reservoir‐scale models for multiphase flow. Issue 9 (23rd September 2013)

• Record Type: Journal Article
• Authors: Sheng, Qiang; Thompson, Karsten
• Is Part Of: Water resources research. Volume 49:Issue 9(2013:Sep.)
• Page Start: 5973
• Page End: 5988
• Publication Date: 2013-09-23
• Subjects: Hydrology -- Periodicals; 333.91
• DOI: 10.1002/wrcr.20430
• Languages: English
• ISSNs: 0043-1397

Abstract

<abstract abstract-type="main"> <title> <x xml:space="preserve">Abstract</x> </title> <p>[1] The concept of coupling pore‐scale and continuum‐scale models for subsurface flow has long been viewed as beneficial, but implementation has been slow. In this paper, we present an algorithm for direct coupling of a dynamic pore‐network model for multiphase flow with a traditional continuum‐scale simulator. The ability to run the two models concurrently (exchanging parameters and boundary conditions in real numerical time) is made possible by a new dynamic pore‐network model that allows simultaneous injection of immiscible fluids under either transient‐state or steady‐state conditions. Allowing the pore‐scale model to evolve to steady state during each time step provides a unique method for reconciling the dramatically different time and length scales across the coupled models. The model is implemented by embedding networks in selected gridblocks in the reservoir model. The network model predicts continuum‐scale parameters such as relative permeability or average capillary pressure from first principles, which are used in the continuum model. In turn, the continuum reservoir simulator provides boundary conditions from the current time step back to the network model to complete the coupling process. The model is tested for variable‐rate immiscible displacements under conditions in which relative permeability depends on flow rate, thus demonstrating a situation that cannot be modeled using a traditional approach. The paper discusses numerical challenges with this approach, including the fact that there is not a way to explicitly force pore‐scale phase saturation to equal the continuum saturation in the host gridblock without an artificial constraint. Hurdles to implementing this type of modeling in practice are also discussed.</p> </abstract>