研究地下水流相关问题时,地下水模型构建中的一个突出问题是勘探资料分布不均或缺乏,以及含水系统本身的非均质性,造成了水文地质参数的空间变异性。在通常的地下水数值模拟中只能用简单的参数分区来描述参数的非均质性,但由于分区过大而忽略分区内部参数的差异性,往往导致模拟结果的不确定性。寻求一种尽可能利用有限的勘探资料,对未知区域内的含水层参数进行合理估值的方法,是目前大区域地下水流模拟中的关键问题。自然情况下,含水层参数不仅具有随机性,也具有一定的结构性。地质统计学就是研究这种具有"二重性"的区域化变量的数学工具,其以变异函数为基本工具来研究分布于空间并呈现一定结构性与随机性的自然现象。作为地质统计学的基本工具,变异函数不仅影响区域化变量的结构分析,还将决定插值结果的精度。。为探究在空间变异分析过程中变异函数模型的选取问题,以及变异函数模型的不同对克里格插值或模拟的影响,本文分别基于球状模型、指数模型和高斯模型,利用非条件模拟生成随机场,通过Monte Carlo法分别获得样本数为50、100、150和200的样本,再利用球状模型、指数模型和高斯模型估计不同样本的变异函数模型参数,最后基于这3种模型及对应的参数估计值进行条件模拟,将模拟结果与原始值对比以评价其模拟精度。同时,本文以新疆焉耆盆地和静县内一农用水源地第三含水层渗透系数采样数据为例,对比分析了高斯模型、球状模型、指数模型在拟合实验变异函数时的差异,以及这种差异对克里格插值的影响。研究结果显示:(1)采用非条件模拟生成随机场的方法存在非遍历性问题,即单次模拟实现的数据统计特征与期望值会有所偏差,但多次实现的平均值与期望值接近。(2)对于球状模型和指数模型,非条件模拟实现的变异函数模型曲线绝大部分平均分布在标准模型曲线上下两侧,且基台值的平均值等于原始模型,但块金常数和变程值要略大,高斯模型非条件模拟实现的变异函数模型曲线全部位于标准模型曲线左侧,变程值相较于原始模型总体明显偏小,块金常数略微偏大。(3)采用最小二乘法或者GS+软件自动拟合实验变异函数时,得到的变程值指数模型>球状模型>高斯模型,块金常数高斯模型>球状模型>指数模型,基台值指数模型>球状模型=高斯模型。(4)采用与原始模型相同函数形式的模型进行变异函数拟合,拟合得到的变异函数模型参数与初始场变异函数模型参数最为接近,球状模型次之。(5)当采样点数较少时,基于指数模型的条件模拟结果总是具有最高的精度、球状模型其次、高斯模型最差;当样本点数较多时,总体上,基于与原始模型相同函数形式的变异函数模型的条件模拟结果精度最高,基于球状模型的条件模拟结果精度仅次于原始模型。(6)在选用多种模型拟合实验变异函数时,变异函数模型的差异主要体现在变程值的不同:变程值通过影响实测点对待估点的作用大小从而影响克里格插值结果:变程值过小,用来进行估值的实测点之间以及实测点与估值点之间相关性降低甚至消失,克里格法退化为简单的算术平均;变程值过大,参与克里格计算的实测点多,插值结果趋近于平稳;当变程值相对适中,克里格插值结果与实测值吻合较好,插值精度也达到最大。