Abstract:Understanding fluid flows in heterogeneous porous media is fundamental to applied geosciences. The wide connectivity variations in the natural aquifer or oil reservoirs make the equivalent permeability have strong spatial variations. When performing the simulations for subsurface flows, the permeabilities may have strong discontinuities across the interfaces between different grid cells. Utilizing the traditional numerical schemes to simulate flows in strong heterogeneous media, the refinement ratio for the grid cell needs to increase dramatically to get an accurate result. Recently, we proposed a finite analytic numerical scheme to solve the two-dimensional fluid flows in heterogeneous porous media. With only 2×2 or 3×3 subdivisions, this scheme can provide rather accurate solutions. In this paper, we develop the finite analytic numerical method for solving the three-dimensional fluid flows in heterogeneous porous media. For the rectangular grid system, it is generally proposed that the pressure gradient in a plane normal to the edge joining different permeability regions will tend to infinite as approaching the edge according to a typical power-law solution and the tangential derivate of the pressure along the edge must be of limited value due to the pressure continuity. Consequently, the three-dimensional flow will reduce to the two-dimensional one in the neighborhood around each edge. Such quasi-two-dimensional behavior is then applied to construct a finite analytic numerical scheme. Numerical examples show that the proposed scheme can provide rather accurate solutions with only 2×2×2 or 3×3×3 subdivisions and the convergent speed is independent of the permeability heterogeneity. Due to its high calculation efficiency, the proposed scheme is utilized to test the well known LLM (Landau, Lifshitz and Matheron) conjecture, which provides keq/kG=exp⁡(16σln⁡k2) for the isotropic log-normal porous medium. The numerical results do not support this conjecture for large σln⁡k, but strongly suggest the linear relation keq/kG=1+16σln⁡k2.

Numerical simulation for porous medium equation by local discontinuous Galerkin finite element method

Local discontinuous galerkin method for radial porous flow with dispersion and adsorption

Finite Analytic Numerical Method for Two-Dimensional Fluid Flow in Heterogeneous Porous Media

A locally mass-conservative enriched Petrov-Galerkin method without penalty for the Darcy flow in porous media

Discontinuous Galerkin Method for Level Set Equation on Unstructured Grids

Numerical analysis of a class of penalty discontinuous Galerkin methods for nonlocal diffusion problems

Fully Discrete Local Discontinuous Galerkin Approximation for Time-Space Fractional Subdiffusion/Superdiffusion Equations

A Local Discontinuous Galerkin Methods with Local Lax-Friedrichs Flux and Modified Central Flux for One Dimensional Nonlinear Convection-Diffusion Equation

Two Finite Element Approaches For The Porous Medium Equation That Are Positivity Preserving And Energy Stable

Numerical methods for non－linear coupled problems in porous media

Analysis Of A Local Discontinuous Galerkin Method For Time-Fractional Advection-Diffusion Equations

Numerical Solutions of Convection-Diffusion Equations by Hybrid Discontinuous Galerkin Methods

Local Discontinuous Galerkin Methods for Diffusive–viscous Wave Equations

Numerical Simulation of Nonlinear Flow Based on Element-Free Galerkin Method

Positivity-preserving and energy-dissipating discontinuous Galerkin methods for nonlinear nonlocal Fokker-Planck equations

Convergence analysis of a numerical scheme for the porous medium equation by an energetic variational approach

A Discontinuous Galerkin Finite Element Method for Dynamic and Wave Propagation Problems in Non‐linear Solids and Saturated Porous Media

A Numerical Method Based on Fully Discrete Direct Discontinuous Galerkin Method for the Time Fractional Diffusion Equation

Finite Element Approximation of the Fractional Porous Medium Equation

A Positivity‐preserving Lagrangian Discontinuous Galerkin Scheme with Exact Riemann Solver for Gas‐water Compressible Flows