基于MODFLOW-CFP的岩溶水模型降雨非线性入渗补给研究——以湖南省香花岭地区为例

杨郑秋; 杨杨; 邵景力; 苏春田; 崔亚莉; 罗飞

doi:10.11932/karst20190504

基于MODFLOW-CFP的岩溶水模型降雨非线性入渗补给研究——以湖南省香花岭地区为例

doi: 10.11932/karst20190504

1.
北京师范大学，北京 100875
2.
中国地质科学院岩溶地质研究所/自然资源部、广西岩溶动力学重点实验室，广西桂林 541004
3.
中国地质大学（北京），北京 100083

基金项目: 国家自然科学基金项目（41502252）；中国地质调查局地质调查项目（DD20160303、DD20160301）

计量
- 文章访问数: 1977
- HTML浏览量: 611
- PDF下载量: 388
- 被引次数: 0
出版历程
- 发布日期: 2019-10-25

Study on non-linear rainfall infiltration recharge of numerical karst water model based on MODFLOW-CFP：A case study of Xianghualing area, Hunan Province

1.
College of Water Sciences，Beijing Normal University，Beijing 100875, China
2.
Institute of Karst Geology,CAGS/Key Laboratory of Karst Dynamics, MNR&GZAR, Guilin, Guangxi 541004, China
3.
China University of Geosciences(Beijing), Beijing 100083, China

摘要

摘要: 利用MODFLOW-CFP建立湖南省香花岭渗流-管道流耦合模型，并使用降雨量系数法实现非线性入渗过程，以探讨此方法在模型中的适用性。通过人工试错反演参数，得出6个不同降雨强度区间的降雨量系数；非线性入渗处理后，地下河流量模拟结果的纳什系数提高至0.91。结果表明：降雨量系数法可使该模型更好地模拟地下河流量变化，此方法也适用于一些岩溶小流域的数值模拟工作。
- 岩溶水 /
- 数值模拟 /
- MODFLOW-CFP /
- 非线性降雨入渗
Abstract: Carbonate rocks are widely distributed in southwest China, and these strata constitute sophisticated karst aquifers, which are local important water and tourism resources. However, groundwater system has strong heterogeneity there, which makes it difficult to characterize subterranean runoff. At the same time, the characteristic that the karst aquifer can be recharged by rainfall infiltration through different ways has not been considered in most of the numerical simulation work. The study of karst groundwater and solute transport plays a key role in the development and utilization of karst groundwater and in the prevention and control of pollution. The purpose of this study is to discuss the relationship between rainfall intensity and rainfall recharge infiltration coefficient in the numerical model, and to use practical method to deal with the variation of coefficient to make the simulation more accurate.The Xianghualing area is located in Chenzhou City of Hunan Province, with low latitude. The aquifer in the study area is composed of Carboniferous limestone. It is a relatively independent karst basin surrounded by mountains with an area of about 14km2, and the surface is covered by luxuriant vegetation. The average annual temperature is 17.9 °C. The annual rainfall ranges from 1，000 mm to 2，400 mm, but most of the rainfall events usually occurs in summer, and there are often rainstorms between July and October. The groundwater level fluctuates with the change of rainfall conditions, and the groundwater table is too low in most areas that the evaporation barely exists. The main artificial extraction of groundwater are mine drainage, the drainage tunnel intersects with the underground river in the discharge area.According to the binary theory of karst system, the karst aquifer can be generalized into two subsystems. The fissure can be generalized into the equivalent porous media system recharged by diffuse infiltration. It is assumed that the groundwater moves in the form of seepage. Underground rivers and karst large-size passages are generalized as conduit system recharged by point infiltration, where is the main discharge pathway of karst groundwater system. MODFLOW-CFPM1 was used to develop the equivalent pore-conduit coupling model in Xianghualing, and its CRCH module was used to simulate the centralized point infiltration. Then the nonlinear infiltration process was represented by rainfall coefficient method and its applicability in the model was well discussed. Depend on data size, six rainfall coefficients of different rainfall intensity intervals were set up and obtained by artificial trial-and-error inversion parameters. After the nonlinear infiltration treatment, the Nash coefficient of the simulation results of underground river flow is increased to 0.91. The results turn out that the rainfall coefficient method can make the model better simulate the change of underground river discharge, and the method is also suitable for the numerical simulation of some karst small watershed.
- karst water /
- numerical simulation /
- MODFLOW-CFP /
- non-linear rainfall infiltration recharge

HTML

参考文献(14)

[1]	赵良杰,夏日元,杨杨,等.基于MODFLOW的岩溶管道水流模拟方法探讨与应用［J］.中国岩溶,2017,36(3):346-351.
[2]	Shoemaker W B, Kuniansky E L, Birk S, et al .Documentation of a Conduit Flow Process (CFP) for MODFLOW-2005：U.S. Geological Survey Techniques and Methods, Book6, Chapter A 24, 50 p.
[3]	Reimann T, Hill M E. MODFLOW‐CFP: A New C-onduit Flow Process for MODFLOW–2005［J］. Groundwater, 2009, 47(3):321-325.
[4]	Hill M E. An evaluation of conduit conceptualizations and model performance［D］.Tampa Bay：University of South Florida, 2008.
[5]	Gallegos J J, Hu B X, Davis H. Erratum: Simulating flow in karst aquifers at laboratory and sub-regional scalesusing MODFLOW-CFP［J］. Hydrogeology Journal, 2013, 21(8):1911-1912.
[6]	陆海建.广西平果铝赤泥堆场对岩溶含水层污染风险评价［D］. 桂林：桂林理工大学,2015.
[7]	洪淑娜, 蓝俊康, 陈丽娜,等. 利用GMS模拟预测崇左市拟建稀土分离项目影响区岩溶水污染［J］. 桂林理工大学学报, 2014,34(1):99-106.
[8]	张鑫, 徐世光, 黄建国,等. GMS在地下水数值模拟中的应用：以云南文山州某水库为例［J］. 河南科学, 2015,33(11):1994-1997.
[9]	易连兴,夏日元,王喆,等.岩溶峰丛洼地区降水入渗系数：以寨底岩溶地下河流域为例［J］.中国岩溶, 2017,36(4)，512-517.
[10]	黄丹红, 成建梅, 刘军,等. 岩溶含水层降雨非线性入渗补给的处理方法［J］.地下水,2006,28(2):23-25.
[11]	陈崇希.滞后补给权函数：降雨补给潜水滞后性处理方法［J］. 水文地质工程地质,1998(6):22-24.
[12]	刘廷玺, 朝伦巴根, 马龙,等. 规定时段降雨入渗补给系数的确定［J］. 内蒙古农业大学学报(自然科学版), 2002, 23(2):44-48.
[13]	常勇. 裂隙—管道二元结构的岩溶泉水文过程分析与模拟［D］. 南京：南京大学, 2015.
[14]	Ghasemizadeh R, Hellweger F, Butscher C, et al. Re-view: Groundwater flow and transport modeling of karst aquifers, with particular reference to the North Coast Lime-stone aquifer system of Puerto Rico［J］. Hydrogeology Journal, 2012, 20(8):1441-1461.