Efficient triple-grid multiscale finite element method for 3D groundwater flow simulation in heterogeneous porous media

In this paper, an efficient triple-grid multiscale finite element method (ETMSFEM) is proposed for 3D groundwater simulation in heterogeneous porous media. The main idea of this method is to employ new 3D linear base functions and the domain decomposition technique to solve the local reduced elliptical problem, thereby simplifying the base function construction process and improving the efficiency. Furthermore, by using the ETMSFEM base functions, this method can solve Darcy’s equation with high efficiency to obtain a continuous velocity field. Therefore, this method can considerably reduce the computational cost of solving for heads and velocities, which is crucial for large-scale 3D groundwater simulations. In the application section, we present numerical examples to compare the ETMSFEM with several classical methods to demonstrate its efficiency and effectiveness.

[1]  Xinguang He,et al.  A stochastic dimension reduction multiscale finite element method for groundwater flow problems in heterogeneous random porous media , 2013 .

[2]  Jichun Wu,et al.  Efficient Triple-Grid Multiscale Finite Element Method for Solving Groundwater Flow Problems in Heterogeneous Porous Media , 2016, Transport in Porous Media.

[3]  Björn Engquist,et al.  Multiscale methods for the wave equation , 2007 .

[4]  J. Banavar,et al.  The Heterogeneous Multi-Scale Method , 1992 .

[5]  Yalchin Efendiev,et al.  Numerical Homogenization of Nonlinear Random Parabolic Operators , 2004, Multiscale Model. Simul..

[6]  Thomas Y. Hou,et al.  A Multiscale Finite Element Method for Elliptic Problems in Composite Materials and Porous Media , 1997 .

[7]  Zhiming Chen,et al.  A mixed multiscale finite element method for elliptic problems with oscillating coefficients , 2003, Math. Comput..

[8]  Yalchin Efendiev,et al.  Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media , 2014, J. Comput. Phys..

[9]  Yalchin Efendiev,et al.  Sparse Generalized Multiscale Finite Element Methods and their applications , 2015, 1506.08509.

[10]  Marcus J. Grote,et al.  Finite Element Heterogeneous Multiscale Method for the Wave Equation , 2011, Multiscale Model. Simul..

[11]  A. Ortega-Guerrero,et al.  Evolution of long‐term land subsidence near Mexico City: Review, field investigations, and predictive simulations , 2010 .

[12]  A modified inverse procedure for calibrating parameters in a land subsidence model and its field application in Shanghai, China , 2016, Hydrogeology Journal.

[13]  T. Arbogast.  Implementation of a Locally Conservative Numerical Subgrid Upscaling Scheme for Two-Phase Darcy Flow , 2002 .

[14]  Clayton V. Deutsch,et al.  Geostatistical Software Library and User's Guide , 1998 .

[15]  Shujun Ye,et al.  Application of the multiscale finite element method to flow in heterogeneous porous media , 2004 .

[16]  Jichun Wu,et al.  The effects of artificial recharge of groundwater on controlling land subsidence and its influence on groundwater quality and aquifer energy storage in Shanghai, China , 2016, Environmental Earth Sciences.

[17]  A. Massoudieh.  Inference of long‐term groundwater flow transience using environmental tracers: A theoretical approach , 2013 .

[18]  Liangsheng Shi,et al.  Multiscale-finite-element-based ensemble Kalman filter for large-scale groundwater flow , 2012 .

[19]  Yalchin Efendiev,et al.  Generalized multiscale finite element methods (GMsFEM) , 2013, J. Comput. Phys..

[20]  Mehrdad T. Manzari,et al.  Effects of using altered coarse grids on the implementation and computational cost of the multiscale finite volume method , 2013 .

[21]  Jørg E. Aarnes,et al.  On the Use of a Mixed Multiscale Finite Element Method for GreaterFlexibility and Increased Speed or Improved Accuracy in Reservoir Simulation , 2004, Multiscale Model. Simul..

[22]  Yalchin Efendiev,et al.  An adaptive GMsFEM for high-contrast flow problems , 2013, J. Comput. Phys..

[23]  R. Rigon,et al.  A mass‐conservative method for the integration of the two‐dimensional groundwater (Boussinesq) equation , 2013 .

[24]  Dongxiao Zhang,et al.  A multiscale probabilistic collocation method for subsurface flow in heterogeneous media , 2010 .

[25]  Thomas Y. Hou,et al.  Convergence of a multiscale finite element method for elliptic problems with rapidly oscillating coefficients , 1999, Math. Comput..

[26]  Yalchin Efendiev,et al.  Adaptive multiscale model reduction with Generalized Multiscale Finite Element Methods , 2016, J. Comput. Phys..

[27]  P. Kitanidis,et al.  Fast iterative implementation of large‐scale nonlinear geostatistical inverse modeling , 2014 .

[28]  Jichun Wu,et al.  Land subsidence in China , 2005 .

[29]  Yalchin Efendiev,et al.  Mixed Generalized Multiscale Finite Element Methods and Applications , 2014, Multiscale Model. Simul..

[30]  Todd Arbogast,et al.  A two-scale numerical subgrid technique for waterflood simulations , 2002 .

[31]  Olof Runborg,et al.  Multi-scale methods for wave propagation in heterogeneous media , 2009, 0911.2638.

[32]  Mehdi Ghommem,et al.  Mode decomposition methods for flows in high-contrast porous media. Global-local approach , 2013, J. Comput. Phys..

[33]  Chunhong Xie,et al.  Cubic-Spline Multiscale Finite Element Method for Solving Nodal Darcian Velocities in Porous Media , 2015 .

[34]  Srinivasulu Ale,et al.  Long term (1960–2010) trends in groundwater contamination and salinization in the Ogallala aquifer in Texas , 2014 .

[35]  Thomas Y. Hou,et al.  Flow based oversampling technique for multiscale finite element methods , 2008 .

[36]  Yuqun Xue,et al.  Mechanical modeling of aquifer sands under long-term groundwater withdrawal , 2012 .

[37]  Yalchin Efendiev,et al.  Accurate multiscale finite element methods for two-phase flow simulations , 2006, J. Comput. Phys..

[38]  A cubic‐spline technique to calculate nodal Darcian velocities in aquifers , 1994 .

[39]  Eric T. Chung,et al.  A numerical homogenization method for heterogeneous, anisotropic elastic media based on multiscale theory , 2015 .

[40]  H. Tchelepi,et al.  Multi-scale finite-volume method for elliptic problems in subsurface flow simulation , 2003 .

[41]  Gour-Tsyh Yeh,et al.  On the computation of Darcian velocity and mass balance in the finite element modeling of groundwater flow , 1981 .

[42]  Yalchin Efendiev,et al.  Homogenization of nonlinear random parabolic operators , 2005, Advances in Differential Equations.

[43]  Li Ren,et al.  Finite volume multiscale finite element method for solving the groundwater flow problems in heterogeneous porous media , 2005 .

[44]  Jichun Wu,et al.  Modified Multiscale Finite-Element Method for Solving Groundwater Flow Problem in Heterogeneous Porous Media , 2014 .

[45]  E Weinan,et al.  The heterogeneous multi-scale method for homogenization problems , 2005 .

[46]  Nicolaos Theodossiou,et al.  Application of Non-Linear Simulation and Optimisation Models in Groundwater Aquifer Management , 2004 .