| Citation: |
LU Wen, LU Chuiyu, HE Xin, SUN Qingyan, ZHANG Bo, JIA Yangwen. Comparative study on two drying-rewetting algorithms of groundwater model cells[J]. Hydrogeology & Engineering Geology, 2024, 51(5): 22-34. doi: 10.16030/j.cnki.issn.1000-3665.202311026
|
Comparative study on two drying-rewetting algorithms of groundwater model cells
-
Abstract
When simulating drying-rewetting process of grid cells in numerical groundwater modeling using the block-centered finite-difference approach, the models is highly probable to run into non-convergence, which could greatly affect the applicability of groundwater models. This study used ideal case and practical simulation in Denmark to comprehensively compare the simulation capabilities and characteristics of two algorithms, namely the empirical trial (ET) method proposed by MODFLOW and the always active cell (AAC) method proposed by COMUS, in the drying-rewetting simulation. The results show that: (1) For the ET method, the selection of parameter combination has a significant influence on the model convergence and the simulation results. It is compulsory to continuously optimize the parameter to avoid model failures such as non-convergence and large simulation errors when using the ET method, which greatly increases the difficulty of groundwater model application and time cost. (2) The simulation results from the AAC method are more reliable than those from the ET method. Theoretically, the effect of the AAC method is equivalent to the optimal parameter set in the ET method. Therefore, parameter optimization is no longer needed in the AAC method, which can effectively reduce the difficulty of using groundwater models and meanwhile reduce the uncertainty of simulation results. (3) The numerical accuracy of the intercell horizontal hydraulic conductance in the AAC method is consistent with that of the classical harmonic average method, demonstrating that the AAC method can also be used in the simulation without the drying-rewetting process. In summary, the AAC method is more suitable for simulating the drying-rewetting process of groundwater model cells and is expected to be more widely used in groundwater numerical simulation.
-
-
References
[1] MERRITT M L. A rewetting approximation for a simulator of flow in a surficial aquifer overlain by seasonally inundated wetlands[J]. Groundwater,1994,32(2):286 − 292. doi: 10.1111/j.1745-6584.1994.tb00643.x [2] 郭子豪. 黄土丘陵沟壑区典型沟道土地整治工程对水系平衡影响研究[D]. 杨凌:中国科学院大学(中国科学院教育部水土保持与生态环境研究中心),2021. [GUO Zihao. The impacts of typical gully land consolidation project on water system balance in Loess Hilly and Gully Region[D]. Yangling:Research Center for Eco-environments and Soil and Water Conservation,Chinese Academy of Sciences & Ministry of Education,2021. (in Chinese with English abstract)] GUO Zihao. The impacts of typical gully land consolidation project on water system balance in Loess Hilly and Gully Region[D]. Yangling: Research Center for Eco-environments and Soil and Water Conservation, Chinese Academy of Sciences & Ministry of Education, 2021. （in Chinese with English abstract） [3] ELZEHAIRY A A,LUBCZYNSKI M W,GURWIN J. Interactions of artificial lakes with groundwater applying an integrated MODFLOW solution[J]. Hydrogeology Journal,2018,26(1):109 − 132. doi: 10.1007/s10040-017-1641-x [4] HERRERA P A,LANGEVIN C,HAMMOND G. Estimation of the water table position in unconfined aquifers with MODFLOW 6[J]. Groundwater,2023,61(5):648 − 662. doi: 10.1111/gwat.13270 [5] KOURAKOS G,HARTER T. Simulation of unconfined aquifer flow based on parallel adaptive mesh refinement[J]. Water Resources Research,2021,57(12):1 − 25. [6] PREZIOSI E,GUYENNON N,PETRANGELI A B,et al. A stepwise modelling approach to identifying structural features that control groundwater flow in a folded carbonate aquifer system[J]. Water,2022,14(16):1 − 16. [7] PAINTER S,BAŞAĞAOĞLU H,LIU Angang. Robust representation of dry cells in single-layer MODFLOW models[J]. Groundwater,2008,46(6):873 − 881. [8] KEATING E,ZYVOLOSKI G. A stable and efficient numerical algorithm for unconfined aquifer analysis[J]. Groundwater,2009,47(4):569 − 579. [9] JEPPESEN J,CHRISTENSEN S,LADEKARL U L. Modelling the historical water cycle of the Copenhagen Area 1850–2003[J]. Journal of Hydrology,2011,404(3/4):117 − 129. [10] LIN Lin,YANG Jinzhong,ZHANG Bin,et al. A simplified numerical model of 3-D groundwater and solute transport at large scale area[J]. Journal of Hydrodynamics,2010,22(3):319 − 328. doi: 10.1016/S1001-6058(09)60061-5 [11] FITTS C R. Modeling dewatered domains in multilayer analytic element models[J]. Groundwater,2018,56(4):557 − 561. doi: 10.1111/gwat.12645 [12] HUNT R J,FEINSTEIN D T. MODFLOW-NWT:Robust handling of dry cells using a Newton formulation of MODFLOW-2005[J]. Groundwater,2012,50(5):659 − 663. doi: 10.1111/j.1745-6584.2012.00976.x [13] BEDEKAR V,NISWONGER R G,KIPP K,et al. Approaches to the simulation of unconfined flow and perched groundwater flow in MODFLOW[J]. Groundwater,2012,50(2):187 − 198. doi: 10.1111/j.1745-6584.2011.00829.x [14] NISWONGER R G,PANDAY S,IBARAKI M. MODFLOW-NWT,a newton formulation for MODFLOW-2005[R]. Reston: U. S. Geological Survey Techniques and Methods Report,2011. [15] 张博. 跨流域调水下苏干湖盆地水循环响应数值模拟[D]. 北京:中国水利水电科学研究院,2018. [ZHANG Bo. Numerical simulation of water cycle response of the Sugan Lake Basin under the influence of inter-basin water transfer project[D]. Beijing:China Institute of Water Resources and Hydropower Research,2018. (in Chinese with English abstract)] ZHANG Bo. Numerical simulation of water cycle response of the Sugan Lake Basin under the influence of inter-basin water transfer project[D]. Beijing: China Institute of Water Resources and Hydropower Research, 2018. （in Chinese with English abstract） [16] LU Chuiyu. C++ object-oriented model for underground water simulation[R]. Beijing:China Institute of Water Resources and Hydropower Research,2023. [17] HARBAUGH A W. MODFLOW-2005:The U S geological survey modular ground-water model-the ground-water flow process[R]. Reston: U. S. Geological Survey Techniques and Methods,2005. [18] LU Chuiyu,HE Xin,ZHANG Bo,et al. Comparison of numerical methods in simulating lake-groundwater interactions:Lake Hampen,Western Denmark[J]. Water,2022,14(19):1 − 13. [19] KIDMOSE J,ENGESGAARD P,NILSSON B,et al. Spatial distribution of seepage at a flow-through lake:Lake Hampen,Western Denmark[J]. Vadose Zone Journal,2011,10(1):110 − 124. doi: 10.2136/vzj2010.0017 [20] OLIVEIRA OMMEN D A,KIDMOSE J,KARAN S,et al. Importance of groundwater and macrophytes for the nutrient balance at oligotrophic Lake Hampen,Denmark[J]. Ecohydrology,2012,5(3):286 − 296. doi: 10.1002/eco.213 [21] KIDMOSE J,ENGESGAARD P,OMMEN D A O,et al. The role of groundwater for lake-water quality and quantification of N seepage[J]. Groundwater,2015,53(5):709 − 721. doi: 10.1111/gwat.12281 [22] LU Chuiyu,ZHANG Bo,HE Xin,et al. Simulation of lake-groundwater interaction under steady-state flow[J]. Groundwater,2021,59(1):90 − 99. doi: 10.1111/gwat.13033 [23] REGNERY J,LI Dong,LEE J,et al. Hydrogeochemical and microbiological effects of simulated recharge and drying within a 2D meso-scale aquifer[J]. Chemosphere,2020,241:1 − 9. [24] WU Ming,WU Jianfeng,LIN Jin,et al. Evaluating the interactions between surface water and groundwater in the arid mid-eastern Yanqi Basin,northwestern China[J]. Hydrological Sciences Journal,2018,63(9):1313 − 1331. doi: 10.1080/02626667.2018.1500744 [25] SONG Jian,YANG Yun,SUN Xiaomin,et al. Basin-scale multi-objective simulation-optimization modeling for conjunctive use of surface water and groundwater in northwest China[J]. Hydrology and Earth System Sciences,2020,24(5):2323 − 2341. [26] PANDAY S,LANGEVIN C D,NISWONGER R G,et al. MODFLOW-USG version 1:An unstructured grid version of MODFLOW for simulating groundwater flow and tightly coupled processes using a control volume finite-difference formulation[R]. Reston: U. S. Geological Survey Techniques and Methods Report Book 6,2013. -
Access History
Figures(10)
Tables(3)
Export File
Citation
LU Wen, LU Chuiyu, HE Xin, SUN Qingyan, ZHANG Bo, JIA Yangwen. Comparative study on two drying-rewetting algorithms of groundwater model cells[J]. Hydrogeology & Engineering Geology, 2024, 51(5): 22-34. doi: 10.16030/j.cnki.issn.1000-3665.202311026
Format
Content
DownLoad:
-
Figure 1.
Schematic diagram of the ET method
-
Figure 2.
Schematic diagram of calculating the average effective saturated thickness between adjacent grid cells in the same layer using the AAC method
-
Figure 3.
Groundwater flow between adjacent grid cells in the same layer using the AAC method
-
Figure 4.
Model cell for which storage factors are used during different states under the AAC method
-
Figure 5.
Comparisons of simulated groundwater heads in the ideal case
-
Figure 6.
Information of the Lake Hampen basin in the practical simulation in Denmark [18]
-
Figure 7.
Model configuration of the Lake Hampen basin
-
Figure 8.
Comparisons of the numbers of the wet cells identified by the ET schemes and the AAC scheme
-
Figure 9.
Comparisons of head errors simulated by the ET schemes and the AAC scheme
-
Figure 10.
Comparison of simulated daily average water table and daily lake water level in the ET schemes and the AAC scheme