Abstract

The temporal discretization scheme is one important ingredient of efficient simulator for two‐phase flow in the fractured porous media. The application of single‐scale temporal scheme is restricted by the rapid changes of the pressure and saturation in the fractured system with capillarity. In this paper, we propose a multi‐scale time splitting strategy to simulate multi‐scale multi‐physics processes of two‐phase flow in fractured porous media. We use the multi‐scale time schemes for both the pressure and saturation equations; that is, a large time‐step size is employed for the matrix domain, along with a small time‐step size being applied in the fractures. The total time interval is partitioned into four temporal levels: the first level is used for the pressure in the entire domain, the second level matching rapid changes of the pressure in the fractures, the third level treating the response gap between the pressure and the saturation, and the fourth level applied for the saturation in the fractures. This method can reduce the computational cost arisen from the implicit solution of the pressure equation. Numerical examples are provided to demonstrate the efficiency of the proposed method.

Highlights

  • The fractured porous media is composed of two distinct pore spaces: the fracture and the matrix

  • We focus on two-phase fluid flow in Journal of Applied Mathematics fractured porous media with capillarity

  • A multiscale time splitting strategy is introduced for simulating two-phase flow in fractured porous media

Read more

Summary

Introduction

The fractured porous media is composed of two distinct pore spaces: the fracture and the matrix. The explicit scheme used for updating the saturation suffers from Courant number stability constraints, and as a result, we should use the smallest time-step size to match the rapid changes in the fracture network if using the single-scale time-stepping To resolve this problem, we should use the multiscale time-stepping for the pressure and the saturation in the fracture and the matrix. In the subtiming methodology for implicittype time-stepping schemes proposed in 39 , the fully implicit schemes with multiple time-step sizes are utilized to the different portions of a domain, which have different accuracy requirements This technique has been applied to the time-dependent flow and transport simulations in fractured porous media 40 , which uses the implicit subtimestepping in the locations within and near the fractures. Taking account into the different physical properties of the matrix blocks and the fracture network, and the multiphysics processes of two-phase flow, we apply multiple time scales for the pressure and saturation equations.

Governing Equations of Two-Phase Flow
Discrete-Fracture Model
Multi-Scale Time Splitting Strategy
Multi-Scale Explicit Time-Stepping for the Saturation Equation
Numerical Tests
Example 1
Example 2
Example 3
Example 4
Conclusions
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call