Global Buckley--Leverett for Multicomponent Flow in Fractured Media: Isothermal Equation-of-State Coupling and Dynamic Capillarity

This paper presents an isothermal Global Buckley--Leverett framework for multicomponent multiphase flow in fractured media that integrates equation-of-state phase behavior, Maxwell--Stefan diffusion, and dynamic capillarity to resolve hyperbolicity issues and ensure a well-posed problem while maintaining interpretability for applications like carbon storage and geothermal exchange.

The Gist

Imagine you are trying to predict how a complex mixture of fluids (like oil, water, gas, and dissolved chemicals) moves through a giant, sponge-like rock formation deep underground. This rock isn't just a simple sponge; it's a fractured sponge, meaning it has tiny pores and large, jagged cracks running through it.

This paper presents a new, upgraded "rulebook" for simulating that movement. The authors call it the Global Buckley–Leverett (GBL-N) framework.

Here is the breakdown of what they did, using simple analogies:

1. The Old Rulebook vs. The New One

The Old Way (Classical Buckley–Leverett):
Think of the old method like a traffic light system for a simple two-lane road (just water and oil). It's fast, easy to understand, and works great when the traffic is light and the road is smooth.

• The Problem: As soon as you add a third lane (a third fluid phase), or if the road is bumpy, cracked, or the cars are changing their shape (chemical reactions), the old traffic lights break down. The math gets "stuck" (it loses hyperbolicity), meaning the computer can't decide which way the fluids should go, leading to chaotic or impossible results.

The New Way (GBL-N):
The authors built a smart, adaptive navigation system that keeps the simplicity of the old traffic lights but adds the brains to handle complex, real-world chaos. It's like upgrading from a paper map to a GPS that knows about traffic jams, road closures, and weather.

2. The Five "Superpowers" Added

To make this work for fractured, complex reservoirs, they added five essential physics "superpowers" that the old model ignored:

• The Shape-Shifter (Equation of State):

  • Analogy: Imagine the fluids aren't just water and oil; they are like shapeshifters. Depending on the pressure and temperature, a gas might turn into a liquid, or a liquid might dissolve into a gas.
  • What the paper does: It tracks these transformations in real-time. If the pressure drops, the gas expands; if it rises, it compresses. The model knows exactly how the "shapeshifters" behave.

• The Social Mixer (Maxwell–Stefan Diffusion):

  • Analogy: In the old model, fluids moved in separate lanes. In reality, molecules are like people at a crowded party. They bump into each other, swap places, and push one another around. A heavy molecule might slow down a light one, or a fast one might drag a slow one.
  • What the paper does: It accounts for this "crowded party" effect. It calculates how different chemical components push and pull each other as they move, which is crucial for things like carbon storage where gases mix with brine.

• The Sticky Sponge (Dynamic Capillarity):

  • Analogy: Think of the rock pores as tiny straws. When you pour water into a dry straw, it doesn't just flow instantly; it has to "wet" the walls first. This takes time. The old model assumed this happened instantly.
  • What the paper does: It adds a "time delay" to the wetting process. This small delay is actually a mathematical magic trick that stops the simulation from crashing. It smooths out the jagged edges of the math, ensuring the computer always finds a unique, stable answer.

• The Stress-Sensitive Rock (Stress-Sensitive Permeability):

  • Analogy: Imagine the rock is made of squishy foam. If you squeeze the foam (increase pressure from above), the holes get smaller, and fluid can't flow through as easily.
  • What the paper does: It acknowledges that the rock changes shape under pressure. As you pump fluids in or out, the rock's ability to let fluids pass (permeability) changes dynamically.

• The Fracture Highway (Non-Darcy Flow):

  • Analogy: In the tiny pores, fluid moves like a slow, polite line of people (Darcy flow). But in the big cracks (fractures), the fluid is a rushing river. If the water moves too fast, it creates turbulence and bumps into the rough walls, losing energy.
  • What the paper does: It treats the cracks differently. It adds a "friction penalty" for high-speed flow in fractures, ensuring the model doesn't overestimate how fast fluids can race through the cracks.

3. How It Solves the "Traffic Jam" Problem

The biggest achievement of this paper is solving the mathematical crash that happens with three or more fluids.

• The Issue: In the old math, when three fluids meet, the equations sometimes become "elliptic" (like a circle instead of a line). This means the future state of the system depends on the future in a confusing way, making it impossible to predict what happens next.
• The Fix: By combining the "Social Mixer" (diffusion) and the "Sticky Sponge" (dynamic capillarity), the authors turned the math into a pseudo-parabolic system.
  • Simple Translation: They added just enough "smoothing" and "friction" to the equations to ensure they always have a clear, unique path forward. It's like adding a shock absorber to a car; the ride is smoother, and the car never flips over.

4. The "Global Pressure" Trick

To keep the simulation fast, they use a Global Pressure concept.

• Analogy: Imagine a team of runners (the different fluid phases). Instead of calculating the speed of every single runner individually, you calculate the speed of the team leader (Global Pressure).
• Once you know where the leader is going, you can easily figure out where the rest of the team is, based on their specific strengths (mobility) and weaknesses (gravity, capillarity). This allows the computer to solve the problem much faster without losing accuracy.

Why Does This Matter?

This model is a practical backbone for some of the most important energy challenges today:

1. Carbon Capture: Storing CO2 underground involves injecting gas into rock that already has water and oil. The gas changes phase and mixes with everything. This model predicts exactly where that CO2 will go.
2. Geothermal Energy: Pumping hot water through fractured rock requires understanding how the rock cracks open or close under stress.
3. Contaminant Cleanup: If toxic chemicals leak into the ground, this model helps predict how they spread through complex, cracked soil.

The Bottom Line

The authors took a classic, simple tool (Buckley–Leverett) and upgraded it with modern physics to handle the messy, complex reality of fractured rocks and chemical mixtures. They managed to keep the tool easy to interpret (you can still see the "traffic flow") while making it mathematically robust (it won't crash) and physically accurate (it accounts for stress, diffusion, and phase changes).

It's the difference between using a ruler to measure a winding river versus using a drone with a laser scanner: the result is the same shape, but the drone gives you a much more accurate, reliable, and useful picture of the whole system.

Technical Summary

1. Problem Statement

The classical Buckley–Leverett (BL) framework is the standard for modeling two-phase, immiscible, incompressible flow in porous media. However, it fails to address the complexities of modern unconventional reservoirs (shales, tight formations) and carbon storage applications, which involve:

• Multicomponent and Multiphase Flow: Systems with $N_c$ components and $N_p$ phases (often $N_p \geq 3$).
• Loss of Hyperbolicity: In three-phase systems, the governing equations often lose strict hyperbolicity, leading to non-unique solutions (umbilic points, elliptic regions) and numerical instability.
• Complex Physics: Real-world scenarios require Equation-of-State (EOS) phase behavior, multicomponent diffusion (Maxwell–Stefan), dynamic capillarity, stress-sensitive permeability, and non-Darcy (inertial) flow in fractures.
• Fractured Media: The need to model matrix-fracture exchange and high-velocity inertial losses in fractures without losing the physical interpretability of the BL fractional-flow concept.

The authors aim to develop a Global Buckley–Leverett (GBL-N) framework that retains the intuitive fractional-flow split of the classical BL model while rigorously incorporating these advanced physical mechanisms to ensure mathematical well-posedness.

2. Methodology

The authors formulate a thermodynamically consistent, isothermal, conservative transport system for $N_c$ components and $N_p$ phases in both matrix and fracture continua.

A. Governing Equations and State Variables

• Conservation: The system uses strictly conservative molar balances for each component $i$, accounting for mobile pore storage, adsorbed mass on the solid, and matrix-fracture exchange.
• EOS Coupling: An isothermal Equation-of-State (EOS) flash map relates the overall composition ($z$), pressure ($p$), and temperature ($T$) to phase properties (saturation $S_\alpha$, mole fractions $x_{i\alpha}$, density $\rho_\alpha$, viscosity $\mu_\alpha$). This ensures thermodynamic consistency across phase transitions.
• Diffusion: Multicomponent diffusion is modeled using Maxwell–Stefan (MS) theory, which accounts for cross-coupling between species and reduces to Fickian diffusion only in limiting cases. A constraint ensures zero net molar diffusive flux relative to the phase average.

B. Momentum and Flow Laws

• Dynamic Capillarity: The capillary pressure is decomposed into an equilibrium term (based on pore-morphology models like Alonso-Marroquín) and a rate-dependent dynamic term ($\tau_\alpha \partial_t S_\alpha$). This term is crucial for regularization.
• Non-Darcy Flow: In fractures (and potentially high-rate matrix zones), the model includes Forchheimer inertial corrections. This introduces a per-phase damping factor ($\chi_\alpha$) that reduces the effective mobility at high velocities.
• Stress Sensitivity: Permeability and porosity are modeled as exponential functions of effective stress, capturing the closure of fractures and pores under pressure.

C. Global Pressure and Fractional-Flow Split

• Global Pressure ($p_g$): The authors define a scalar global pressure driven by a mobility-weighted sum of phase potentials. This decouples the pressure equation (elliptic/weakly parabolic) from the transport equations.
• Fractional-Flow Decomposition: The phase velocity is split into:
  1. A fractional-flow term (mobility-weighted share of the total flux).
  2. Drift terms representing buoyancy (density differences) and capillary forces.
  3. Damping: Inertial effects are folded into "apparent mobilities" ($\tilde{\lambda}_\alpha$), preserving the BL structure while correcting flux magnitudes.

D. Matrix-Fracture Coupling

The model supports both Dual-Porosity (Warren-Root) and Embedded Discrete Fracture Models (EDFM/pEDFM). Exchange fluxes are driven by phase potential differences (including gravity and dynamic capillarity) and are scaled by the fracture-side damping factor to maintain strict conservation.

3. Key Contributions

1. Restoration of Well-Posedness: The paper mathematically demonstrates that the combination of Maxwell–Stefan diffusion and dynamic capillarity renders the transport system pseudo-parabolic. This resolves the loss of strict hyperbolicity inherent in classical three-phase BL models, ensuring unique and stable solutions even in regions where eigenvalues coalesce.
2. Unified Framework: It integrates EOS phase behavior, multicomponent diffusion, dynamic capillarity, stress sensitivity, and non-Darcy flow into a single "Global BL" structure.
3. Exact Fractional-Flow Split with Inertia: Unlike previous attempts that struggle with non-Darcy flow, this formulation derives an exact operational split for phase velocities that includes Forchheimer damping, maintaining the interpretability of the BL fractional-flow curve.
4. Generalized Total Differential (gTD) Condition: The authors clarify the conditions under which a global pressure decoupling is fully equivalent to the original system. When the strict TD condition fails, they propose a conservative $H^1$-projection surrogate to approximate the global pressure potential.
5. Coordinate-Free and Cylindrical Forms: The equations are presented in a coordinate-free vector form for generality, with specific derivations for axisymmetric (cylindrical) radial injection, making it directly applicable to CO2 storage and geothermal scenarios.

4. Results and Analysis

• Classical Limit: The model rigorously reduces to the classical Buckley–Leverett equations when diffusion, dynamic capillarity, and inertia are disabled, validating its consistency.
• Energy Estimates: An energy inequality is derived showing that MS diffusion controls composition gradients ($\nabla z$) and dynamic capillarity controls saturation time-derivatives ($\nabla \partial_t S$). This proves the system is strictly well-posed.
• Fracture Modeling: The inclusion of Forchheimer damping in fractures is shown to be critical for high-rate flows, preventing unphysical flux predictions while maintaining conservation at matrix-fracture interfaces.
• Numerical Implementation Strategy: The paper outlines a sequential iterative algorithm:
  1. EOS Flash (update properties).
  2. Solve Global Pressure equation.
  3. Reconstruct Phase Fluxes (using the split).
  4. Advance Conservative Transport (including diffusion and adsorption).

5. Significance

This work provides a practical, thermodynamically consistent backbone for simulating complex flow in fractured, compositional reservoirs.

• For Carbon Capture and Storage (CCS): It offers a robust tool for modeling CO2 injection, where multicomponent mixing, phase changes, and dynamic capillarity are dominant, and where the loss of hyperbolicity in standard simulators can lead to numerical failure.
• For Unconventional Reservoirs: It bridges the gap between pore-scale physics (slip, surface diffusion) and Darcy-scale simulation, allowing for more accurate predictions in shale and tight gas formations.
• Theoretical Impact: By proving that dynamic capillarity and MS diffusion are the necessary regularizers for multicomponent multiphase flow, it resolves a long-standing mathematical issue in reservoir simulation theory.

In summary, the paper successfully modernizes the Buckley–Leverett framework, transforming it from a simplified two-phase tool into a rigorous, general-purpose engine for next-generation energy and environmental applications.