不确定性是客观世界的固有属性,普遍存在于各个学科领域。对地下水模型进行不确定性分析,一直是水文地质学的研究热点和难点。最近十年,多模型分析(Multiple Model Analysis, MMA)的方法开始受到学者们的重视,并得到了广泛的应用。本文从岩性分布、渗透系数随机分布和岩体复杂构造三方面开展了多模型分析。首先,选取常见的河间地块问题,以渗流区岩性分布的不确定性为研究对象,构建了多模型方案,模拟得到了预测潜水位及其置信区间。岩性分布不确定性由置信区间的变化得到了很好的体现,验证了多模型分析方法的实用性。基于非稳定井流溶质运移问题,评价了建模过程中网格单元尺度选择和渗透系数场随机生成过程中抽样、插值等不确定因素对模拟的影响。分析了不同位置观测孔的浓度变化特征以及污染羽形态差异,结果表明网格单元尺度对模拟结果影响最大,抽样误差次之,克里格插值方法最小。研究还指出了AICc准则评价网格尺度、渗透系数随机分布等不确定性的局限性。针对雅砻江锦屏二级水电站工程区的隧洞涌水问题,在MODFLOW环境下建立了5种工况下的三维稳定流模型,讨论了复杂构造所引起的参数方案不确定性,利用UCODE程序优化了参数,使用多模型方法预测了防渗前后隧洞群的总涌水量分别为40.35 m~3/s和31.19 m~3/s。包括渗透结构在内的各种水文地质条件的不确定性,都应该作为多模型分析的对象。在系统分析概念模型不确定性的基础之上,选用分布式地下水数值模拟程序MODFLOW,建立了太原地区岩溶地下水系统的多模型分析框架,全面、系统地研究了1956年至1994年气候变化和人类活动影响下的晋祠泉和兰村泉泉流量衰减过程。研究取得了几点认识:1)将东山岩溶水系统纳入研究区的模拟效果更好;2)晋祠、兰村泉域地下水系统不存在所谓的可变边界;3)渗透系数垂向上衰减程度的差异所引起的不确定性大于补给方式;4)晋祠泉对降雨量变化较敏感,降水量减少和开采增加都是致使晋祠断流的重要原因。网格尺度不确定性对模拟有着重要影响,但往往被研究者忽视,为了探讨评价此类不确定性的方法,进行了数值试验。试验分析了模型网格大小不同所导致的渗透系数分布的差异,以及布置井、河流和观测孔等要素的变化,基于AICc的多模型分析结果表明,传统的信息量准则并不适于评价网格尺度因素引起的模拟不确定性。为此,提出了基于网格尺度的信息量准则GIC (Grid Information Criterion),实现了对AICc准则的改进。