The Laplace Transform Galerkin Technique for large‐scale simulation of mass transport in discretely fractured porous formations

The ability to simulate contaminant migration in large‐scale porous formations containing a complex network of discrete fractures is limited by traditional modeling approaches. One primary reason is because of vastly different transport time scales in different regions due to rapid advection along the discrete fractures and slow but persistent diffusion in the porous matrix. In addition to time‐related complexities, standard numerical methods require a fine spatial discretization in the porous matrix to represent sharp concentration gradients at the interface between the fractures and the matrix. In order to circumvent these difficulties, the Laplace transform Galerkin method is extended for application to discretely fractured media with emphasis on large‐scale modeling capabilities. The technique avoids time stepping and permits the use of a relatively coarse grid without compromising accuracy because the Laplace domain solution is relatively smooth and devoid of discontinuities even in advection‐dominated problems. Further computational efficiency for large‐grid problems is achieved by employing a preconditioned, ORTHOMIN‐accelerated iterative solver. A unique feature of the method is that each of the several needed p space solutions are independent, thus making the scheme highly parallel. Other features include the accommodation of advective‐dispersive transport in the porous matrix and the straightforward inclusion of dual‐porosity theory to represent matrix diffusion in regions where microfractures exist below the modeling scale. An example problem involving contaminant transport through an aquitard into an underlying aquifer leads to the conclusion that deep, essentially undetectable fractures in clayey aquitards can greatly compromise the quality of groundwater in the impacted aquifer.