河套灌区地下水位的空问分布规律研究.pdfVIP

第二届黄河国际论坛论文集 3 河套灌区地下水位的空问分布规律研究* 王卫光1张仁铎1，2 (1．武汉大学水资源与水电工程科学国家重点实验室； 2．Wyoming大学，IJaramie) 摘要：运用地质统计学中的普通克里金方法，对1997。2001年河套灌区203个观测井的年均 井水位高程进行了插值及空问场变异分析。分析发现地下水位高程分布符合J卜态分布，采用 克里金插值方法为计算工具进行分析，结果显示试验半变异蛹数符合一般的球状模型。从获 得的参数町知，1997—2001年期间，水位高程的相关尺度各年均存在着表现相似的儿何各向 异性，在搜索方向为310。左右时为变程长轴，变程大约134．98km，在40。左右为相应的短轴， 变程大约为60 km。同时各个年份计算结果均硅示有块金值，说明在小于变程的空间尺度上 仍然存在随机效应。最后，运用克里金方法对该河套灌医的地下水位高程进行了空间插值成 图。 关键词：地下水位地质统计学空间变异克里金方法 1引言 为了定量描述随机变量的空间变异，最初的研究根据传统统计学原理，不考 虑测定参数测定位置之间的空间关系，计算测定区域变量某一性质的统计学参 数，由此考察在这一区域变量性质的变异程度，这一方法的缺点是明湿的，它只 能估计测定区域所考察变量变异程度的强弱，不能估计测定区域所考察变量的 空间分布。为了克服这一缺点，以后的研究着重于变量空问分布的估计。克里 金方法是由南非的地质学家、采矿工程师克里金于1951年提出的，后来法国学 者马特隆于1962年提出“区域化变量”的概念后，得到不断发展和完善。它是一 种特定的滑动加权平均法，这种方法假定采样点间的距离和方向可反映一定的 空间关联，并利用它们来解释空间的变异性。它主要是利用一定的数学函数对 特点或是给定搜索半径内的所有点进行拟合来估计每个点的值，对处理含有距 离和方向上的偏差的数据尤为适用。地质统计学是以区域化变量等理论为基 础，以变异函数为基本工具来研究展布于空间并呈现一定结构性与随机性的自 * 国家自然科学基金和黄河基金联合资助项目(502a9090)。 4 然现象的科学，发展到现今已不单单应用于地质找矿和矿产估量中，它的使用范 围由逐步涉及环境科学、水科学、海洋学、农田水利、土壤学、森林及渔业等领域， 发展成为用来研究自然界具有随机性和规律性双重特征变量的具有普遍意义的 科学方法。在水科学中，克里金方法已应用于土壤水、径流模数的估值，并且在 布设水文站点与气象站点时，克里金方法也提供了强有力的理论依据。地下水 位是具有随机性与规律性双重特征的区域化变量，本文以河套灌区1997．2001 年的年均地下水位高程为对象，运用克里金方法分析了该灌区降地下水位的空 间变异特征。 2资料和方法 2．1地质统计学理论 2．1．1半方差函数 半方差函数的定义是地质统计学方法的理论基础，半方差函数是描述随机 变量空间变异结构的一个函数。假设随机函数均值稳定，方差存在且有限，该值 仅和间距有关，则半方差函数用以下的数学形式来表达 y(h)={E{[z(石+h(一z(戈)]2}(1) 计算式为 1 ％ r(^)=疣蓦[z(”¨一z(戡)]2 (2) 式中^——分离距离； ％——在间距为h的样本“对”数，它的下标h表示％是分离距离的函 数。 常用的半方差模型包括球形模型、指数模型、高斯模型和线性模型。根据样 本数据计算样本半方差值后，可使用“交叉证实法”来选择半方差函数模型。具 体的步骤为，轮流将每一点数据从样本中暂时去除，用其他剩余的样本值以及克 里金法(在克里金法中使用初步选择的变异函数模型)来估计该点的值，直到对 所有样本点进行了估值。然后，用统计方法对一一对应的观测值和估计值进行 比较，选择模型。容易看出，实测样本数据和估计值之间的误差式是所选半方差 模型是否合理的一种量度。用来评估半方差模型合理性的一些指标有：克里金 的平均误差尽可能接近于0，均方差尽可能的小，约