我国岩溶分布十分广泛,类型复杂繁多,约占全国陆地面积的三分之一以上。岩溶地质体中发育的裂隙、管道、溶洞等岩溶形态不仅为地下水、矿产等资源提供了赋存的空间,也是发生岩溶塌陷、岩溶水突水等地质环境问题的场所。虽然国内有许多学者对岩溶含水层采用数值模拟的方法进行了研究,但对层流-紊流共存流场中岩溶裂隙网络演化规律的研究较少。随着溶蚀过程的不断进行,裂隙网络中的隙宽不断被扩宽,地下水流的流态逐渐从层流向紊流过渡,当裂隙中的地下水流为紊流时,裂隙网络中的水头分布、流量、以及裂隙扩宽速率与层流相比有很大的差别。本研究利用有限差分数值模拟技术,提出了一个层流-紊流共存流场中岩溶裂隙网络演化数值模拟的方法。改进了课题组开发的岩溶裂隙数值模拟软件FractureToKarst,利用C++语言在水头计算部分新增了紊流模块,实现了该软件在层流-紊流共存流场中进行岩溶裂隙网络演化模拟的功能。采用非连续介质分析方法,基于二维光滑平行板模型假定,构建岩溶含水系统二维裂隙网络渗流模型。对裂隙网络进行离散化处理后,将渗流模型与溶解动力学模型进行耦合。首先,利用水均衡原理建立裂隙网络各节点水头的方程组。当裂隙段中的水流为层流时,流量采用立方定律表示;当裂隙段中的水流为紊流时,流量采用Lomize G(1951)提出的紊流公式表示,并采用迭代法求得各节点水头值和各裂隙段的流量。之后,依据Dreybrodt W(1990,1996)提出的碳酸盐岩溶解速率经验公式求解各裂隙段的溶解速率,根据质量守恒原理求得该时间步长内裂隙的扩宽量,最终求得溶蚀后的隙宽。求出新的隙宽后,便可重复以上步骤,重新计算水头分布,求得经过下一时间步长后的隙宽,从而得到岩溶裂隙网络的演化过程。本研究分别设计了承压岩溶含水层与潜水岩溶含水层紊流模型对软件模拟的正确性进行了验证,并且模拟了岩溶含水层溢流泉和古岩溶盆地两个模型中裂隙网络的演化过程。研究发现,只接受降雨补给的溢流泉模型和古岩溶盆地模型直到分水岭消失时也没有紊流的发生。接受河流与降雨补给的溢流泉模型在1200年时,在泉的排泄点处首先发生了紊流,并且紊流随着演化时间逐渐向系统内部发展。在整个岩溶系统演化的过程中,在潜水位附近岩溶发育程度最为强烈。