基于MODFLOW的地下水流模型前处理优化

吉林大学学报(地球科学版)

基于MODFLOW的地下水流模型前处理优化

韩忠¹，邵景力¹，崔亚莉¹，程汤培^2,3，李玲¹，杨程¹

1.中国地质大学（北京）水资源与环境学院，北京100083;
2.北京应用物理与计算数学研究所计算物理重点实验室，北京100088;
3.北京应用物理与计算数学研究所高性能计算中心，北京100094

收稿日期:2013-10-05 出版日期:2014-07-26 发布日期:2014-07-26
通讯作者: 邵景力（1959-），男，教授，主要从事水文学及水资源专业的教学和科研工作 E-mail:jshao@cugb.edu.cn
作者简介:韩忠(1988-)，男，硕士，主要从事地下水数值模拟方面的研究，E-mail:hanzhonggl@163.com
基金资助:
国家“973”计划项目（2010CB428804）

Preprocessing Optimization of Groundwater Flow Model Based on MODFLOW

Han Zhong¹, Shao Jingli¹, Cui Yali¹, Cheng Tangpei^2,3, Li Ling¹,Yang Cheng¹

1.School of Water Resources and Environment, China University of Geosciences, Beijing100083, China;
2.Key Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Beijing100088, China;
3.High Performance Computing Center, Institute of Applied Physics and Computational Mathematics, Beijing100094, China

Received:2013-10-05 Online:2014-07-26 Published:2014-07-26

摘要/Abstract

摘要：

MODFLOW中现有的降水补给数据和井流数据是基于每个剖分网格输入的，导致基于MODFLOW建立大区域水流数值模型前处理的降水补给文件和井流文件存储开销过大、读取效率低。为此，通过改进现有的降水补给子程序包RCH及开发新的子程序包RAW，给出了一种基于面状补排项的数据输入新方法。改进后的MODFLOW程序通过读取RCH文件或RAW文件中每个分区的补排量数据，以及每个网格对应的分区编号，在程序内部实现了补排量向每个模型网格的分配。RAW子程序包实现了多层面状补排量的表达，可用于面状的地下水开采、农业回灌等源汇项的处理。相对于原始的源汇项数据存储方式，基于新方法建立的华北平原地下水流模型，RCH及RAW文件大小分别减小为原来的1/145和1/255，整个模型数据的读取时间的加速比为7.46。

关键词: MODFLOW, RCH, RAW, 大区域地下水模拟, 华北平原

Abstract:

Reading and writing precipitation and well data are all based on each grid in the existing MODFLOW. So building regional groundwater flow numerical models with refined spatial and temporal discretization often involves storing and reading large quantities of well and precipitation data. However, the storage requirement can be prohibitive and the program efficiency can be quite low. To address this problem, we have modified the existing package RCH and also developed a new package named RAW (recharge and well), which is based on the existing packages in MODFLOW. The improved MODFLOW program can only read values of every subarea and the subarea number of each grid which is stored in the RCH or RAW file, and then allocates the values to corresponding grids. With the RAW package, the area recharge and discharge could occur in all layers. Therefore, it could be used to simulate the source and sink term, like groundwater exploitation and agricultural irrigation return. The correctness has been verified and the efficiency has been tested through simulating a theoretical groundwater flow model and a real North China Plain model respectively. Compared with the traditional method of storing data, the size of RCH and RAW file decreases 1/145 and 1/255 respectively, and the time ratio for reading data is 7.46.

Key words: MODFLOW, RCH, RAW, large-scale groundwater simulation, North China Plain

中图分类号:

P641.1

韩忠，邵景力，崔亚莉，程汤培，李玲，杨程. 基于MODFLOW的地下水流模型前处理优化[J]. 吉林大学学报(地球科学版), 2014, 44(4): 1290-1296.

Han Zhong, Shao Jingli, Cui Yali, Cheng Tangpei, Li Ling,Yang Cheng. Preprocessing Optimization of Groundwater Flow Model Based on MODFLOW[J]. Journal of Jilin University(Earth Science Edition), 2014, 44(4): 1290-1296.

参考文献

Metrics

Viewed

Full text

Abstract

Cited

Shared

Discussed

[1]	束龙仓, 许杨, 吴佩鹏. 基于MODFLOW参数不确定性的地下水水流数值模拟方法[J]. 吉林大学学报(地球科学版), 2017, 47(6): 1803-1809.
[2]	吴鸣, 吴剑锋, 施小清, 刘杰, 陈干, 吴吉春. 基于谐振子遗传算法的高效地下水优化管理模型[J]. 吉林大学学报(地球科学版), 2015, 45(5): 1485-1492.
[3]	孙建国. Kirchhoff型偏移理论的研究历史、研究现状与发展趋势展望——与光学绕射理论的类比、若干新结果、新认识以及若干有待于解决的问题[J]. J4, 2012, 42(5): 1521-1552.
[4]	张兆吉, 费宇红, 郭春艳, 钱永, 李亚松. 华北平原区域地下水污染评价[J]. J4, 2012, 42(5): 1456-1461.
[5]	李文运, 崔亚莉, 苏晨, 张伟, 邵景力. 天津市地下水流-地面沉降耦合模型[J]. J4, 2012, 42(3): 805-813.
[6]	卢文喜，李平，王福林，郭龙珠. 挠力河流域三维地下水流数值模拟[J]. J4, 2007, 37(3): 541-0545.