回顾了岩溶区地下水流运动数值模拟的发展历程，概述了岩溶水运动的基本理论，系统全面地总结了国内外有关岩溶区地下水运动数值模拟方法，分析比较了各种计算方法的优缺点，在此基础上给出了各种计算方法的适用条件。 在仔细分析中梁水利水电工程库区水文地质条件的基础上，结合地表水测流资料、水化学分析资料及示踪连通试验资料，得出水库区岩溶渗漏条件，并在此基础上提出水库渗漏概念模型，采用水均衡法计算出库区的岩溶水渗漏量。 结合库区具体水文地质条件，讨论了三种岩溶区地下水流数值模型：线性渗流模型，非线性渗流模型和线性—非线性渗流模型。 线性渗流模型结合了中梁水库库区水文地质条件，重点进行正反问题的计算分析，利用研究区若干观测孔资料，采用有限单元法反演库区的水文地质参数，并预测了库区岩溶水渗漏量。 推导出整体坐标系下岩溶管道区非线性渗流模型渗透张量计算公式，基于Galerkin有限元法，导出了地下水非线性渗流的有限元方程，在此基础上编制了计算程序KARCON.FOR，将该程序用于典型模型的计算分析，结果表明该方法是合理的。 考虑到岩溶区渗流往往都是线性流与非线性流并存，首次提出了基于区域分解法理论的线性—非线性渗流模型，利用该方法模拟出理想模型中的水流运动，并编制了计算程序DAR-NONDAR.FOR，然后将该成果应用于中梁水利水电工程中，预测出库区岩溶水渗漏量。 在上述三种计算库区岩溶水渗漏量的模型中，由于线性—非线性渗流模型同时考虑了与中梁水库区实际情况相符合的达西流和非达西流区域并存的地下水运动，因此，利用该模型计算出的库区岩溶水渗漏量更合理。