多孔介质非均质性空间变异结构分析.pdfVIP

第18卷 第4期 地质灾害与环境保护 Vo1．18， No．4 2 0 0 7年1 2月 Journal of Geological Hazards and Environment Preservation December 2007 文章编号： 1006—4362(2007)O4～0074一O5 多孑L介质非均质性空间变异结构分析 谷石岩，张征，赵娟，徐永利 (北京林业大学资源与环境学院，北京 100083) 摘要： 非均质多孔介质的空间分布具有随机性与结构性的双重性质。地统计学方法在研究具有结构 性和随机性双重性变量的空间变异性方面具有极强的科学性。本文在地理信息系统平台ArcGIS上， 采用Geostatistical Analyst地统计学分析模块中的kriging插值方法研究非均质性多孔介质的空间变 异性。以河北省廊坊市地下水污染物的现场数据为研究对象，利用普通克立格法对氟化物、硫酸盐、氯 化物和汞浓度的空间变异性进行分析。 关键词： 普通克立格；空间变异性；地统计学；多孔介质 中图分类号： 0213 文献标识码： A 多孑L介质是由固体骨架和流体组成的一类复合 最大年降水量1 087 mm，最小年降水量273 mm，年 介质。天然土体、含水层、岩层以及油藏都属于多孑L 平均无霜期197 d，最大冻土深度70 cm。区域地貌 介质[1 ]。自然界的多孔介质大多具有非均质 属于河流冲积平原，具有多层水文地质结构，地下水 性[4]，被视为同一质地的区域内，在同一时刻各点的 多以大气降水、灌溉入渗补给、侧向径流补给为主。 均值存在着明显差异，天然土壤和含水层及岩层沉 1．2 区域化变量及变异函数 积过程具有随机性，它们都是不均匀、具有高度空间 地统计学(Geostatistic)是 Matheron于 1963 异质性的复合体，与其内在的母质、水文、生物等都 年创立并发展的一门学科。该方法是以区域化变量 有一定的关系，这使得含水层的水文地质特征以及 为核心和理论基础，以空间相关和变异函数为基本 土壤的结构、构造和土壤各种矿物的组成等土壤特 工具的一种数学地质方法，一般以半方差图作为描 性具有非常明显的空间变异性。地统计学正是研究 述空间变异性的工具 ]。 分布于空间中并显示出一定随机性和结构性的自然 区域变量是指以空间点的三个直角坐标(X ， 现象。与此同时，污染物可以通过土壤、岩体等多孔 x ，X )为自变量的随机场Z(x ，x ，x )一Z(z)。 介质的孔隙渗透到地下，多孑L介质的裂隙是地下水 当对它进行了一次观测后，就得到了它的一个现实 赋存的空间和运动通道，也是污染物迁移的通道 ]。 Z(z)，它是一个普通的三元实值函数或空间点函 因此利用地统计学方法，研究地下水污染物浓度的 数。当空间点z沿方向变化为z+h时，就把区域 空间变异性，可以很好的表达多孔介质非均质性的 变化量在z与z+h处的值Z(z)与Z(x+ )的差 空间结构。 的方差的一半定义为区域变化量Z(z)在z方向上 1 材料与方法 的变异函数，记为7(x， )： 1 7(x， )一÷El Z(z)～Z(z+ )1 0 1．1 研究区概况 ’ 廊坊市位于河北省中部偏东，属