Numerical simulation of precipitation impact on fractured karst system evolution
-
摘要: 本文在构建有降雨入渗及河流补给的裂隙型岩溶含水系统的概念模型基础上,采用有限差分数值方法,利用裂隙渗流立方定律及Dreybort(1996)经过实验总结出的碳酸钙溶蚀经验公式,定量地计算了裂隙型岩溶含水系统经过10 000年溶蚀的演化形态及岩溶泉的变化。模拟发现,岩溶含水系统的溶蚀主要发生在降雨入渗面处、河流补给处和岩溶泉的出露点附近。随着溶蚀的进行,岩溶含水系统的潜水位不断下降部分岩溶泉被疏干,通过计算岩溶泉汇流管道的隙宽变化,发现岩溶含水系统会出现明显的差异性溶蚀。在模拟过程中调节降雨量的大小发现降雨入渗对有河流补给的岩溶含水系统的演化影响并不显著。Abstract: The 10,000 years evolution of fractured karst system,which has precipitation infiltration and river recharge,is calculated using finite difference method, cubic law and Dreybort's (1996) empirical formula after building on the conceptual model of fractured karst system The result shows that the dissolution of karst aquifer system mainly occurs in the point of precipitation infiltration, recharge rivers and karst springs. With the dissolution of the karst aquifer system, the water tables were declining and some of the karst springs were dry. We found that the number of precipitation infiltration does not significantly affect the evolution of fractured karst system which has river recharge.
-
Key words:
- karst evolution /
- fracture seepage /
- carbonate dissolution rate /
- numerical simulation
-
[1] 王大纯,张人权,史毅虹,等.水文地质学基础[M].北京:地质出版社,1986:131144.- [2] Dreybrodt W.The Role of Dissolution Kinetics in the Development of Karst Aquifers in Limestone: A Model Simulation of Karst Evolution[J].The Journal of Geology,1990,98(5): 639-655. [3] -Palmer A N.Origin and morphology of limestone caves[J].Geological Society of America Bulletin,1991,103(1):1-21. [4] -Groves C G,Howard A D.Early development of karst systems: 1.Preferential flow path enlargement under laminar flow[J].Water Resources Research,1994,30(10):2837-2846. [5] Howard A D,Groves C G.Early Development of Karst Systems: 2.Turbulent Flow[J].Water Resources Research,1995,31(1):19-26. [6] -Siemers J,Dreybrodt W.Early development of Karst aquifers on percolation networks of fractures in limestone[J].Water Resources Research,1998,34(3):409-420. [7] Kaufmann G,Braun J.Karst Aquifer evolution in fractured,porous rocks[J].Journal of Hydrology,2016,35(11):3223-3238. [8] Romanov D.Evolution of Karst Aquifers in Natural and Man Made Environments: A Modeling Approach[D].University of Bremen,2003. [9] Kaufmann G.Karst aquifer evolution in a changing water table environment[J].Water Resources Research,2002,38(6):26-1-26-9. [10] Kaufmann G.Karst conduit evolution[M].Evolution of karst: from prekarst to cessation,2002. [11] Kaufmann G.Numerical models for mixing corrosion in natural and artificial karst environments[J].Water Resources Research,2003,39(6):813-813. [12] Kaufmann G.A model comparison of karst aquifer evolution for different matrix-flow formulations[J].Journal of Hydrology,2003,283(1-4):281-289. [13] Kaufmann G.Structure and evolution of karst aquifers: A finite-element numerical modelling approach[M].Processes in Speleogenesis: A modelling approach,2005. [14] Kaufmann G.Modelling karst geomorphology on different time scales[J].Geomorphology,2009,106(1-2):62-77. [15] Yu Qingchun.Some investigation on early organization of Karst system[J].Jounral of China University of Geosciences,1999,10:314-321. [16] 薛亮,于青春.岩溶含水系统演化过程中河间地块分水岭消失过程的数值模拟分析[J].水文地质工程地质,2009,36(2):7-12. [17] 徐钟济.蒙特卡洛法[M].计算科学丛书,1985. [18] 刘再华,Dreybrodt W.DBL理论模型及方解石溶解/沉积速率预报[J].中国岩溶,1998,17(1):1-7. [19] Plummer L N,Wigley T M L,Parkhurst D L.The kinetics of calcite dissolution in CO2-water systems at 5 ℃ to 60 ℃and 0.0 to 1.0 atm CO2[J].American Journal of Science,1978,278(2):179-216. [20] Dreybrodt W.Principles of early development of karst conduits under natural and man- made conditions revealed by mathematical analysis of numerical models[J].Water Resources Research,1996,32(9):2923-2935. [21] 王云,于青春,薛亮,等.裂隙岩溶含水系统溢流泉演化过程的数值模拟[J].中国岩溶,2010,29(4):378-384.
点击查看大图
计量
- 文章访问数: 1631
- HTML浏览量: 342
- PDF下载量: 1302
- 被引次数: 0