基于TOUGH2的柴达木盆地诺木洪剖面地下水流模拟

干旱区地理 ›› 2016, Vol. 39 ›› Issue (3): 548-554.

基于TOUGH2的柴达木盆地诺木洪剖面地下水流模拟

顾晓敏¹, 张戈², 郝奇琛¹, 邵景力¹, 崔亚莉¹, 张秋兰¹, 肖勇¹

1 中国地质大学(北京)水资源与环境学院, 北京 100083;
2 中国地质调查局西安地质调查中心, 陕西西安 710001

收稿日期:2015-12-11 修回日期:2016-02-27 出版日期:2016-05-25
通讯作者: 邵景力(1959-),男,教授,水文地质专业,主要从事水文地质教学与科研工作.Email:jshao@cugb.edu.cn
作者简介:顾晓敏(1991-),女,博士研究生,水文地质专业,主要从事地下水研究.Email:cugbgxm@163.com
基金资助:
中国地质调查局干旱-半干旱区地下水与生态重点实验室开放基金项目(1212011220974)

Application of TOUGH2 in groundwater numerical simulation of Qaidam Basin

GU Xiao-min¹, ZHANG Ge², HAO Qi-chen¹, SHAO Jing-li¹, CUI Ya-li¹, Zhang Qiu-lan¹, XIAO Yong¹

1 China University of Geosciences(Beijing), School of Water Resource and Environment, Beijing 100083, China;
2 Xi'an Center of Geological Survey, CGS, Xi'an 710001, Shaanxi, China

Received:2015-12-11 Revised:2016-02-27 Online:2016-05-25

摘要/Abstract

摘要： 基于TOUGH2数值模拟软件的EOS9模块建立柴达木盆地典型剖面饱和-非饱和地下水流数值模型。模型中对上边界条件的处理采用考虑潜在蒸发能力和土壤含水率的表土含水量计算方法。利用试估校正法对模型参数进行识别验证,模拟结果得出实测水位与模拟水位拟合相关系数为0.9764;均衡计算结果表明地下水补给量主要为南侧山前补给,占补给总量的90.1%,排泄量主要为溢出排泄和蒸发排泄,分别占64.3%和35.7%。地下水流系统发育三级流动层次,其中局部水流系统排泄是地下水排泄的主要方式,循环量主要为泉溢出排泄,占58.8%;中间流循环量主要为溢出排泄和蒸发排泄,占27.6%;区域流循环量主要为蒸发排泄,占13.6%。

关键词: 柴达木盆地, 饱和-非饱和, 地下水流系统, Tough2

Abstract: In order to accurately and efficiently study the transforming relationship between surface water and groundwater, and the renewal capacity of groundwater system in Qaidam Basin, northwest China, a two-dimensional variably saturated numerical model of a typical profile from Nuomuhong River to Amunike Mountain was developed using TOUGH2 simulator. On the basis of previous studies, the evaporation was a function of soil saturation given as an upper boundary to characterize water transport near surface through iterative calculation. The numerical model was calibrated with 41 groundwater observation data based on trial-and-error method. The simulated results fitted well with observation data. The correlation coefficient of simulated groundwater heads was 0.9764 m, and the absolute error was 1.6 m. The simulation results showed that the typical profile of Qaidam Basin could be divided into three orders of groundwater flow system:local groundwater flow system accounted for 58.8% of total amount of groundwater flow with spring as the main discharge; the middle flow system accounted for 27.6% of the total amount of water, with evaporation and river as the main discharge; regional flow system accounted for 13.6% of total amount of water, with evaporation as the most important discharge. The results showed that TOUGH2 with modified boundary condition can accurately describe water transport in variable saturated zone under the conditions of phreatic fluctuation. The study of the renewal capacity of groundwater and the characteristic of groundwater flow system has great significance for the areas lack of intensive studies. In the future research, different supply sources and forms, as well as dispersion and density changes can be taken into consideration to simulate groundwater ages and further study the groundwater systems.

Key words: Qaidam Basin, saturated-unsaturated, groundwater flow system, TOUGH2

中图分类号:

P641.8

顾晓敏, 张戈, 郝奇琛, 邵景力, 崔亚莉, 张秋兰, 肖勇. 基于TOUGH2的柴达木盆地诺木洪剖面地下水流模拟[J]. 干旱区地理, 2016, 39(3): 548-554.

GU Xiao-min, ZHANG Ge, HAO Qi-chen, SHAO Jing-li, CUI Ya-li, Zhang Qiu-lan, XIAO Yong. Application of TOUGH2 in groundwater numerical simulation of Qaidam Basin[J]. , 2016, 39(3): 548-554.

参考文献

[1] 陈亚宁,杨青,罗毅,等. 西北干旱区水资源问题研究思考[J]. 干旱区地理,2012,35(1):1-9.[CHEN Yaning,YANG Qing, LUO Yi,et al. Ponder on the issues of water resources in the arid region of northwest China[J]. Arid Land Geography,2012,35(1):1-9.]

[2] ANDERSON M P,WOESSNER W W. Applied groundwater modeling:simulation of flow and advectivetransport[M]. Oxford:Gulf Professional Publishing,1992:203-205.

[3] 夏雨波,谢海澜,王冰,等. 地下水循环演化模式及研究方法综述[J]. 地质调查与研究,2012,35(4):299-303.[XIA Yubi, XIE Hailan,WANG Bing,et al. Review on regional groundwater cycle and evolvement[J]. Geological Survey and Research, 2012,35(4):209-303.]

[4] SYKES J F,NORMANI S D,JENSEN M R,et al. Regional-scale groundwater flow in a Canadian shield setting[J]. Canadian Geotechnical Journal,2009,46(7):813-827.

[5] CANDELA L,ELORZA F J,TAMOH K,et al. Groundwater modelling with limited data sets:the Chari-Logone area(Lake Chad Basin,Chad)[J]. Hydrological Processes,2014,28(11):3714-3727.

[6] 魏国孝,王德军,王刚,等. GIS技术和FEFLOW在酒泉东盆地地下水系统数值模拟中的应用[J]. 兰州大学学报:自然科学版,2008,43(6):1-6.[WEI Guoxiao,WANG Dejun,WANG Gang,et al. Application of GIS technique and FEFLOW to numerical simulation of groundwater system in Jiuquan East Basin[J]. Journal of Lanzou University(Natural Sciences),2007,43(6):1-6.]

[7] 张亚宁,张勃,王亚敏,等. 1960-2008年黑河流域参考作物蒸发蒸腾量的时空变化[J]. 干旱区地理,2011,34(5):779-787.[ZHANG Yaning,ZHANG Bo,WANG Yamin,et al. Temporal and spatial changes of reference crop evapotranspiration in the Heihe River Basin from 1960 to 2008[J]. Arid Land Geography, 2011,34(5):779-787.]

[8] 黄金廷,侯光才,尹立河,等. 干旱半干旱区天然植被的地下水水文生态响应研究[J]. 干旱区地理,2011,34(5):788-793.[HUANG Jinting,HOU Guangcai,YIN Lihe,et al. Eco-hydrological response of natural vegetation in arid and semi-arid:a review[J]. Arid Land Geography,2011,34(5):788-793.]

[9] 邵景力,崔亚莉,张德强. 基于包气带水分运移数值模型的黄河三角洲蒸发量研究[J]. 地学前缘,2005,12(4):95-100.[SHAO Jingli,CUI Yali,ZHANG Deqiang. Evapotranspiration from modeling of soil water flow in the Yellow River Delta[J]. Earth Science Frontiers,2005,12(4):95-100.]

[10] 张殿君,张学霞,武鹏飞. 黄土高原典型流域土地利用变化对蒸散发影响研究[J]. 干旱区地理,2011,34(3):400-408.[ZHANG Dianjun,ZHANG Xuexia,WU Pengfei. Relationship between ET and LUCC in a typical watershed of Loess Plateau over the past 20 years[J]. Arid Land Geography,2011,34(3):400-408.]

[11] 施小清,张可霓,吴吉春. TOUGH2软件的发展及应用[J]. 工程勘察,2009,37(10):29-34.[SHI Xiaoqing,ZHANG Keni, WU Jichun. The history and application of TOUGH2 code[J]. Geotechnical Investigation & Surveying,2009,37(10):29-34.]

[12] 尚松浩,毛晓敏,潜水蒸发研究进展[J]. 水利水电科技进展, 2010,30(4):85-89.[SHANG Songhao,MAO Xiaomin. Research progress on evaporation from phreatic water[J]. Advances in Science and Technology of Water Resources,2010,30(4):85-89.]

[13] 童新,刘廷玺,杨大文,等. 半干旱沙地-草甸区水面蒸发模拟及其影响因子辨识[J]. 干旱区地理,2015,38(1):10-17.[TONG Xin,LIU Tingxi,YANG Dawen,et al. Simulating evaporation from a water surface for the sand-meadow ecotone of the semiarid region in North China[J]. Arid Land Geography, 2015,38(1):10-17.]

[14] RICHARDS L A. Capillary conduction of liquids through porous mediums[J]. Journal of Applied Physics,2004,1(5):318-333.

[15] VAN GENUCHTEN M T. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils[J]. Soil Science Society of America Journal,1980,44(5):892-898.

[16] NARASIMHAN T N,WITHERSPOON P A. An integrated finite difference method for analyzing fluid flow in porous media[J].Water Resources Research,1976,12(1):57-64.

[17] 芮孝芳. 水文学原理[M]. 北京:中国水利水电出版社,2012:104-126.[RUI Xiaofang. Hydrology theory[M]. Beijing:China Water & Power Press,2012:104-126.]

[18] GONÇALVES T D,FISCHER T,GRÄBE A,et al. Groundwater flow model of the Pipiripau watershed,Federal District of Brazil[J]. Environmental Earth Sciences,2013,69(2):617-631.

[19] TOTH J. A theoretical analysis of groundwater flow in small drainage basins[J]. Journal of Geophysical Research,1963,68(16):4795-4812.

[20] 周仰效,李文鹏. 区域地下水流模拟[J]. 水文地质工程地质, 2009,36(1):1-10.[ZHOU Yangxiao,LI Wenpeng. Modeling of regional groundwater flow[J]. Hydrogeology & Engineering Geology,2009,36(1):1-10.]

[21] FREEZE R A,WITHERSPOON P A.Theoretical analysis of regional groundwater flow:1. Analytical and numerical solutions to the mathematical model[J]. Water Resources Research,1966, 2(4):641-656.

[22] FREEZE R A,WITHERSPOON P A. Theoretical analysis of regional groundwater flow:2. Effect of water-table configuration and subsurface permeability variation[J]. Water Resources Research, 1967,3(2):623-634.

[23] JIANG X W,WAN L,WANG X S,et al. Effect of exponential decay in hydraulic conductivity with depth on regional groundwater flow[J]. Geophysical Research Letters,2009,36(24):88-113.

[24] 梁杏,牛宏,张人权,等. 盆地地下水流模式及其转化与控制因素[J]. 地球科学:中国地质大学学报,2012,37(2):269-275.[LIANG Xing,NIU Hong,ZHANG Renquan,et al. Basinal groundwater flow patterns and their transformation and do-minant factors[J]. Earth Science-Journal of China University of Geosciences,2012,37(2):269-275.]