The Parabolic Variational Inequalities for Variably Saturated Water Flow in Heterogeneous Fracture Networks

Abstract

Fractures are ubiquitous in geological formations and have a substantial influence on water seepage flow in unsaturated fractured rocks. While the matrix permeability is small enough to be ignored during the partially saturated flow process, water seepage in heterogeneous fracture systems may occur in a non-volume-average manner as distinguished from a macroscale continuum model. This paper presents a systematic numerical method which aims to provide a better understanding of the effect of fracture distribution on the water seepage behavior in such media. Based on the partial differential equation (PDE) formulations with a Signorini-type complementary condition on the variably saturated water flow in heterogeneous fracture networks, the equivalent parabolic variational inequality (PVI) formulations are proposed and the related numerical algorithm in the context of the finite element scheme is established. With the application to the continuum porous media, the results of the numerical simulation for one-dimensional infiltration fracture are compared to the analytical solutions and good agreements are obtained. From the application to intricate fracture systems, it is found that water seepage flow can move rapidly along preferential pathways in a nonuniform fashion and the variably saturated seepage behavior is intimately related to the geometrical characteristics orientation of fractures.