摘要

A finite layer approach for the general problem of three-dimensional (3D) flow to horizontal wells in multilayered aquifer systems is presented, in which the unconfined flow can be taken into account. The flow is approximated by an integration of the standard finite element method in vertical direction and the analytical techniques in the other spatial directions. Because only the vertical discretization is involved, the horizontal wells can be completely contained in one specific nodal plane without discretization. Moreover, due to the analytical eigenfunctions introduced in the formulation, the weighted residual equations can be decoupled, and the formulas for the global matrices and flow vector corresponding to horizontal wells can be obtained explicitly. Consequently, the bandwidth of the global matrices and computational cost rising from 3D analysis can be significantly reduced. Two comparisons to the existing solutions are made to verify the validity of the formulation, including transient flow to horizontal wells in confined and unconfined aquifers. Furthermore, an additional numerical application to horizontal wells in three-layered systems is presented to demonstrate the applicability of the present method in modeling flow in more complex aquifer systems.