ISSN 1000-3665 CN 11-2202/P
    焦友军,黄奇波,王旭升,等. 等水位河间地块岩溶裂隙和流场演变模拟初探[J]. 水文地质工程地质,2024,51(1): 1-11. DOI: 10.16030/j.cnki.issn.1000-3665.202307009
    引用本文: 焦友军,黄奇波,王旭升,等. 等水位河间地块岩溶裂隙和流场演变模拟初探[J]. 水文地质工程地质,2024,51(1): 1-11. DOI: 10.16030/j.cnki.issn.1000-3665.202307009
    JIAO Youjun, HUANG Qibo, WANG Xusheng, et al. Preliminary modeling investigation on changes in Karst fractures and seepage field in rockmass between streams with equal water levels[J]. Hydrogeology & Engineering Geology, 2024, 51(1): 1-11. DOI: 10.16030/j.cnki.issn.1000-3665.202307009
    Citation: JIAO Youjun, HUANG Qibo, WANG Xusheng, et al. Preliminary modeling investigation on changes in Karst fractures and seepage field in rockmass between streams with equal water levels[J]. Hydrogeology & Engineering Geology, 2024, 51(1): 1-11. DOI: 10.16030/j.cnki.issn.1000-3665.202307009

    等水位河间地块岩溶裂隙和流场演变模拟初探

    Preliminary modeling investigation on changes in Karst fractures and seepage field in rockmass between streams with equal water levels

    • 摘要: 经典地下水动力学根据Dupuit假设得到了河渠间潜水含水层水位分布的解析解,常用于实际工程。该解析模型是否适用于具有裂隙网络的河间地块岩溶含水层,这个问题远未得到充分论证。考虑降雨入渗强度为100~800 mm/a的情形,建立剖面二维离散裂隙网络渗流模型,对等水位河间地块开展数值模拟研究。裂隙网络包含陡倾角和缓倾角2组裂隙,具有随机分布特征,隙宽均值为0.01 cm。根据稳定流场模拟结果划分包气带和饱水带裂隙界面,分析水位分布特征,计算饱水带裂隙网络的等效渗透系数,并与经典解析模型反算的等效渗透系数进行了对比,发现采用解析模型的误差小于25%。进一步模拟岩溶作用下裂隙演变的情景,对10 ka内河间地块裂隙状态和准稳态渗流进行了预测,发现隙宽最大值达到0.07 cm,经典解析模型仍然能够通过水位反算出数量级相符的等效渗透系数。模拟结果表明潜水面具有不规则形状,而且在河岸会出现明显的渗出面泉点,裂隙泉点数量在岩溶演变过程中逐渐减少。经典解析模型虽然能够估算河间地块裂隙网络的等效渗透系数,但不能刻画潜水面的不规则形态,也不能预测岩溶裂隙渗流场的演变趋势。

       

      Abstract: The analytical solution of distributed water table in an unconfined aquifer between streams has been obtained in classic groundwater hydraulics based on Dupuit’s assumption, and it is frequently used in engineering practice. However, the applicability of this analytical model to groundwater in karst fractures needs a further verification. A 2D profile model of discrete fractures in rockmass between streams with equal water levels was simulated, considering a rainfall infiltration rate of 100−800 mm/a. The random fracture network includes two fracture groups, steeply or slightly inclining with the same statistic average aperture of 0.01 cm. Phreatic surface were identified by checking the interface between saturated and unsaturated fractures. The equivalent hydraulic conductivity of saturated fractures was obtained from seepage simulations and compared to the inversely estimated result from the classical analytical solution. It indicates that the relative error of inverse estimation is smaller than 25%. Further, the change in fractures aperture due to dissolution of Karst rocks was simulated with predictions on the state of fractures and quasi-steady state seepage in a 10 ka period. It finds that the maximum of the fracture aperture reaches 0.07 cm, and the classical model can yield an equivalent hydraulic conductivity from water table distribution with the same order of magnitude as the real hydraulic conductivity. The simulation shows an irregular shape of water table and seepage springs on river sides, where the number of springs reduced gradually following the karstification process. Though the classical analytical model can be used to estimate the equivalent hydraulic conductivity of the fractured aquifer between streams, it is unable to reveal the irregular water table shape and predict the change in seepage in Karst fractures.

       

    /

    返回文章
    返回