A parallelizable method for two-phase flows in naturally-fractured reservoirs

Jim Douglas*, Felipe Pereira, Li-Ming Yeh

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

17 Scopus citations


A parallelizable, semi-implicit numerical method is proposed for the study of naturally-fractured reservoir systems. It has proved to be computationally efficient in producing accurate numerical solutions for the dual-porosity model for immiscible, two-phase flow in such reservoirs. The method combines hybridized mixed finite elements, a new version of the modified method of characteristics, a sophisticated operator-splitting procedure for separating the pressure calculation in the fractures from that of the saturation, another operator splitting to handle the interaction of the matrix blocks and the fractures, and domain decomposition iterative procedures for both the pressure and the saturation. It permits moderately long time steps for the pressure and the saturation in the fractures and matrix blocks by using short, inexpensive microsteps to treat the transport portion of the saturation equation in the fractures. This paper is devoted to the formulation of the method and a discussion of numerical results for five-spot and vertical cross-section examples.

Original languageEnglish
Pages (from-to)333-368
Number of pages36
JournalComputational Geosciences
Issue number3-4
StatePublished - 1 Sep 1997


  • Mathematical Modeling
  • Industrial Mathematic
  • Iterative Procedure
  • Domain Decomposition
  • Method Combine

Fingerprint Dive into the research topics of 'A parallelizable method for two-phase flows in naturally-fractured reservoirs'. Together they form a unique fingerprint.

Cite this