Discussion and application of simulation methods for karst conduit flow based on MODFLOW
-
摘要: 针对岩溶管道水流特征,选取广西寨底岩溶流域为研究对象,分别采用三维有限差分地下水流动模型(MODFLOW)中的Drain和River模块概化模拟岩溶管道特征,结合管道内上游、中游及下游三个观测孔对比分析两种方法的适用性,结果表明:末期等水位线图中位于管道两侧处流线均有形状突变,符合实际管道水流动态特征,且能模拟出管道内水位变化趋势。最后分析两种方法模拟管道水流的原理,Drain模块中岩溶管道仅起排水作用,不允许管道内水流流向周围含水层,而River模块允许岩溶管道水与周围含水层水流交换,因而相对Drain模块而言,应用River模块概化模拟岩溶管道更加精确。 Abstract: The Zhaidi karst underground river system is located in Guangxi of southwest China and is a typical karst spring catchment with a mean annual discharge of around 1.62 m3/s. The catchment area is about 31.05 km2, with altitudes ranging from 190 m to 820 m a.s.l, and is developed mainly in the Devonian limestones with the landform type of karst peak-cluster depression. Precipitation mainly takes place in spring and summer and mean annual precipitation recorded from 1971 to 2000 is about 1,601.1 mm. The Zhaidi karst underground river system is characterized by extremely anisotropic fractured carbonate rocks and well interconnected karst conduits. There are eight karst underground river sub-systems and many sky-windows, pools and sinkholes. According to tracer tests, these eight subsystems are also well connected each other, where the groundwater runs from the sinkhole G37 through the karst conduits toward the outlet G47. Although numerical technique is a powerful tool for evaluating karst water resources, but it is difficult to accurately depict the pipe flow in the karst conduits in the model due to its anisotropic property. In this paper, the conduit flow of the Zhaidi karst basin was studied, by using Drain module and River module in MODFLOW to simulate the karst conduit, for which the applicability of the modules was examined by the data derived from three observation boreholes respectively situated in the upstream, midstream and downstream of the basin. To observe the groundwater level fluctuation and to obtain the time series data, two boreholes (ZK7 and ZK8) were drilled on the conduit with intention. From upstream to downstream, the distances between sinkhole G37 and borehole ZK7, borehole ZK7 and ZK8, borehole ZK8 to outlet G47 are 1,400 m, 900 m, 140 m, respectively. In addition, the outlet G47 was equipped with a rectangle sharp-crested weir with a level logger installed. The results showed that there was a sudden change of the water table contour near conduit and the water level variation trend was simulated, which meant both of method could conform to the actual conduit features. At last, the principle of two simulation methods were analyzed which showed karst conduit in Drain module only played a drainage role without exchanging water, while in River module it could exchange water with aquifer. It can be concluded that River module is better that Drain module to simulate karst conduit flow.
-
Key words:
- conduit flow /
- karst aquifer /
- numerical simulation /
- MODFLOW
-
[1] 杨杨,唐建生,苏春田,等. 岩溶区多重介质水流模型研究进展[J]. 中国岩溶, 2014,33(4):419-424. [2] 陈皓锐,高占义,王少丽,等. 基于Modflow的潜水位对气候变化和人类活动改变的响应[J]. 水利学报, 2012, 43(3): 96-105, 114. [3] Rahnama M B, Zamzam A. Quantitative and Qualitative Simulation of Groundwater by Mathematical Models in Rafsanjan Aquifer using MODFLOW and MT3DMS[J]. Arabian Journal of Geosciences, 2013, 6(3): 901-912. [4] 崔光中,于浩然,朱远峰. 我国岩溶地下水系统中的快速流[J]. 中国岩溶, 1986, 5(4): 297-305. [5] Loper D E, Chicken E. A leaky-Conduit Model of Transient Flow in Karstic Aquifers[J]. Mathematical Geosciences, 2011, 43(8): 995-1009. [6] 成建梅,陈崇希. 广西北山岩溶管道—裂隙—孔隙地下水流数值模拟初探[J]. 水文地质工程地质, 1998(4): 52-56. [7] 陈崇希,胡立堂. 渗流—管流耦合模型及其应用综述[J]. 水文地质工程地质, 2008(3):76-81. [8] 孙晨,束龙仓,鲁程鹏,等. 裂隙-管道介质泉流量衰减过程试验研究及数值模拟[J]. 水利学报, 2014, 45(1): 54-61, 68. [9] Cheng J M, Chen C X. An integrated linear/non-linear flow model for the conduit-fissure-pore media in the karst triple void aquifer system[J]. Environmental Geology, 2005, 47(2):163-174. [10] Saller S P, Ronayne M J, Long A J. Comparison of a karst groundwater model with and without discrete conduit flow[J]. Hydrogeology Journal, 2013, 21(7): 1555-1566. [11] 陶小虎,赵坚,陈孝兵,等. 岩溶含水层水流模型研究进展[J]. 水利水电科技进展, 2014, 34(2): 80-88. [12] 罗明明,尹德超,张亮,等. 南方岩溶含水系统结构识别方法初探[J]. 中国岩溶, 2015, 34(6): 543-550. [13] Gallegos J J, Hu B X, Davis H. Simulation flow in karst aquifers at laboratory and sub-regional scales using Modflow-CFP[J]. Hydrogeology Journal, 2013, 21(8): 1749-1760. [14] 肖斌,许模,曾科,等. 基于Modflow的岩溶管道概化与模拟探讨[J]. 地下水, 2014, 36(1): 59-61, 84. [15] McDonald M G, Harbaugh A W. A Modular Three Dimensional Finite-difference Groundwater Flow Model [M]. U.S. Geological Survey, 1984. [16] Harbaugh A W, Banta E R, Hill M C, et al. MODFLOW-2000, The U.S. Geological Survey Modular Groundwater Model [M]. U. S. Geological Survey, 2000. [17] 赵良杰,夏日元,易连兴,等. 岩溶地下河浊度来源及对示踪试验影响的定量分析[J]. 地球学报, 2016, 37(2): 241-246.
点击查看大图
计量
- 文章访问数: 2336
- HTML浏览量: 450
- PDF下载量: 1266
- 被引次数: 0