Two-phase flow and transport in a single fracture-porous medium system

Abstract

The prognosis for the remediation of contaminated fractured media is much worse than that for more homogeneous units. Fractures act as conduits for the flow of dense non-aqueous phase liquids (DNAPLs), while diffusion is responsible for the storage of dissolved mass in the surrounding matrix. A numerical model incorporating aqueous phase transport in a variable-aperture fracture and its surrounding matrix is developed and coupled with an existing two-phase flow model. The processes of transient two-phase flow, non-equilibrium dissolution, advective–dispersive transport in the fracture, and three-dimensional matrix diffusion are included in the model. Results from various investigations show that the DNAPL distribution is very sensitive to variations in aperture within a single fracture. Diffusion-controlled mass removal from both the matrix and from the hydraulically inaccessible zones within the fracture itself result in extremely large time frames for significant mass removal from these systems. Success in aqueous phase mass removal from the matrix is very sensitive to the effective fracture spacing. The hydraulic gradient in the fracture only affects the amount of water removed from the system, and does not greatly affect the amount of time required to remove the contaminant mass from the matrix. The ability to remove mass is somewhat sensitive to the porosity and effective matrix diffusion coefficient over the range of expected values.