地质统计学方法在地下水水位估值中的应用.docVIP

应用协同-泛克立格法进行潜水位估值 陈家军( 郭乔羽 王金生 （北京师范大学环境科学研究所，环境模拟与污染控制国家重点联合实验室，北京100875） 摘 要 在许多区域水资源水环境问题中，地下水流或水质模拟往往采用数值方法，需给出各个节点的初始水位值以反映流场的初始状态。本文阐述了，以容易获取的节点地表高程信息，作为协同区域化变量，用协同-泛克立格法进行地下水位估值的原理，以河南省焦作市修武段地下水数值模拟分析区为例，给出了利用该方法时的实验变差函数获取方法，分析了这一方法的估值情况，并与未加入协同区域化变量地表高程的泛克立格法进行了比较，指出协同-泛克立格法在进行松散多孔介质区域地下水位估值时能取得更为良好的效果。 关键词 地下水，水位估值，协同-泛克立格法，变差函数，协同区域化变量. 概述 在区域水资源、水环境问题中，地下水水流模拟和水质模拟往往采用数值方法，需要在每个节点上给出初始流场值作为初始条件，反映流场的初始状态。通常情况下，剖分节点数总是多于统一测量地下水位的观测孔数，因而大部分节点上的初始地下水位要根据统一测量点的地下水位，通过插值求得。地质统计学方法(克立格法是一种最佳空间估计法，其本质是最佳无偏估计，最初起源于地质矿产储量预测并在这一领域得到迅猛发展，近年来也广泛地应用于水资源、水环境问题的相关研究领域，如沼泽中化学物质循环的空间变异规律[1]（Marnette，Stein，1993），水域污染空间模式，海洋微生物[2](Beliaeff，Marie-laure，1995)，含水层导水率的分布[3] (Chu et al., 1982)，污染物的浓度[4] (Dingman et al. 1988)，区域降水量估值[5] (Geren et al., 1994)，降水网络设计[6] (Bastin et al., 1984)，流量估值[7] (Huang et al., 1998)。利用协同-克立格在地下水位估值中也有一定的研究[8-12] (Hoeksema et al.,1984, 1989; Yeh et al., 1995; Parks et al., 1996)，Hoeksema et al.(1989) [10]在地下水位估值研究中所用的协同-克立格方法是协同-普通克立格方法，其中忽略了地下水位变化的趋势特征。本文在泛克立格法的基础上，加入协同区域化变量地表高程，用协同-泛克立格法估计地下水水位，因为对于区域地下水(潜水)其从水位高处流向水位低处且往往与地形标高相一致，这种方法在可获得的国内外文献中尚未见到。提出了通过多次迭代获取剩余的稳健实验变差函数，从而获取区域化变量的稳健实验变差函数的方法。对协同-泛克立格法和单纯的泛克立格法在估值上的准确性进行了比较，表明采用协同-泛克立格法进行地下水位估值能取得更为准确的结果。 地下水水位估值的协同-泛克立格法 在所研究的地下水水流系统中遍布观测孔是不可能的。进行水流数值模拟时，节点数总是多于统一测量水位的观测孔数，而且还有一些统一测量水位点不能与节点重合，因而大多数节点上的初始水位要根据统一测量水位通过插值求得。在节点水位值信息不足的同时，节点的地表高程信息是容易从地形图等途径获得的。利用这些补充信息，使信息量增大，会使地下水位估值的准确性提高。对于松散多孔介质中有潜水分布的区域，地表高程与地下水位具有相关性，即地下水位往往随着地形的高低起伏而同步变化，尽管这不是绝对的，但这一变化趋势是存在的，因而可采用协同-泛克立格法进行地下水位估值。在用这一方法进行估值时，地下水位、地表高程作为协同区域化变量，其中包括结构部分和随机部分，而地下水位作为主估计量。通常区域流场和区域地形都具有一定的趋势，也就是说其期望值不是常数而是空间坐标的函数，协同区域化变量地下水位和地表高程是非平稳的，存在漂移。假定协同区域化变量地下水位、地表高程满足： (1) 式中：(地下水位；(地表高程；(空间坐标（向量）；(的数学期望值；(漂移函数；(方差；(变差函数；k=1代表水位，k=2代表地表高程。 设在流场内有个统一测量地下水位点，个统一测量地表高程点，则任一点地下水位的估计值由下式计算： (2) 式中：(任一点水位的估计值；(协同-泛克立格权系数；(统一测点地下水位值；(统一测点地表高程值. 漂移可写成[13]： (3) 通常可表示成多项式。由于地下水初始流场在一般情况下水力坡度以及地形变化较为平缓，因而在本文中假定平面二维流场，线性漂移，则可写成：