The Lanczos Reduction Method For Long-Term Simulation Of Groundwater Flow In Discretely-Fractured Porous Media

K Zhang,Ad Woodbury

2000-01-01

Abstract:One of the more advanced approaches for simulating groundwater flow in fractured porous media is the discrete fracture approach. However, this approach is limited by the large computational overheads associated with traditional modeling methods. One primary reason is because of the vastly different flow velocities between the fractures and the matrix blocks, which results in sharp groundwater head gradients near the fracture-matrix interface. In addition, for most field problems, the formations contain complex randomly distributed networks of discrete fractures. In order to simulate the sharp hydraulic head gradients at the interface between fractures and the matrix, a fine time and spatial discretization is needed. The discretization can easily involve many tens to hundreds of thousands of unknowns and time steps. In this work, we extend the Lanczos reduction method to the modeling of groundwater now using the discrete fracture approach. The Lanczos reduction method reduces a finite element equation system to a much smaller tridiagonal system of first-order differential equations. The reduced system can be solved by a standard tridiagonal algorithm with little computational effort. Because solving the reduced system is more efficient compared to solving the original system; this method is especially suitable for the problem of large-scale and long-term simulation. In this paper, we develop an iterative version algorithm, in which the preconditioned conjugate gradient solver is employed within the Lanczos reduction process. In addition, a "shift" method is introduced, which can greatly accelerate the iterative solution convergent speed. The model is verified by comparison with dual-porosity approach. The efficiency and accuracy of the method are demonstrated on a field scale problem and compared to the performance of classic time marching using an iterative solver on the original system.