基于协同泛克立格法的地下水位估值方法研究.docxVIP

基于协同泛克立格法的地下水位估值方法研究 在区域水资源和水环境方面，地下水流量模拟和水质模拟通常采用数值法。每个节点都需要给出初始流场值作为初始条件，以反映流量初始状态。通常，横截面上的节点数通常大于统一测量的地下水位观测孔的数量。因此，在大多数节点中，初始地下水位应该通过统一测量点的插值获得。地质统计学方法是最好的空间估算方法，其本质是最好的无偏差估计。它最初是在矿产资源储量预测和水环境研究领域进行的。近年来，它也被广泛应用于水资源和水环境的研究领域，如沼泽中物质循环的空间变异规律、水面污染的空间模型、海洋微生物、包括水层渗透系数的分布、污染物浓度、区域降水量估算、降水网络设计、流量估算、地下水位评价、地下水嵌入深度分布（mart、stie=beiefef、cohard；zhanghai；zhengyen；zhengyeetal.；zhangyeetal.；zhongyeetal.）在这项工作中，我们提出了一种基于泛克立方法的稳定实验变量模型，并使用泛克立方法估计各层地下水的水位。这一方法在所获得的国家文献中是罕见的。提出了一种稳定实验变量差分法，该方法可以通过多次迭代获得区域变量的函数。采用泛克立方法和单泛克立方法的估计精度，结果表明，通过泛克立方法进行地下水位评价可以取得更准确的结果。 1 地下水流场的地表高程 在所研究的地下水水流系统中遍布观测孔是不可能的.进行水流数值模拟时, 节点数总是多于统一测量水位的观测孔数, 而且还有一些统一测量水位点不能与节点重合, 因而大多数节点上的初始水位要根据统一测量水位通过插值求得.在节点水位值信息不足的同时, 节点的地表高程信息是容易从地形图等途径获得的.利用这些补充信息, 使信息量增大, 会使地下水位估值的准确性提高.对于松散多孔介质中有潜水分布的区域, 地表高程与地下水位具有相关性, 即地下水位往往随着地形的高低起伏而同步变化, 尽管这不是绝对的, 但这一变化趋势是存在的, 因而可采用协同-泛克立格法进行地下水位估值.在用这一方法进行估值时, 地下水位H1、地表高程H2作为协同区域化变量, 其中包括结构部分和随机部分, 而地下水位H1作为主估计量.通常区域流场和区域地形都具有一定的趋势, 也就是说其期望值不是常数而是空间坐标的函数, 协同区域化变量地下水位和地表高程是非平稳的, 存在漂移.假定协同区域化变量地下水位、地表高程满足: {E[Ηk(x)-Ηk(y)]=mk(x)-mk(y)12D2[Ηk(x)-Ηk(y)]=γk(x,y)k=1,2(1) 式中H1为地下水位;H2为地表高程;x为空间坐标 (向量) ;E[Hk]为Hk的数学期望值;mk(x) 为漂移函数;D2为方差;γk(x,y) 为变差函数;k=1代表水位,k=2代表地表高程. 设在流场内有n1个统一测量地下水位点,n2个统一测量地表高程点, 则任一点地下水位的估计值H*由下式计算: Η*=n1∑i=1λiΗ1i+n2∑j=1ωjΗ2j(2) 式中H*为任一点水位的估计值;λi,ωj为协同-泛克立格权系数;H1i为统一测点地下水位值;H2j为统一测点地表高程值. 漂移mk(x) 可写成: mk(x)=t∑l=0ak,lfk,l(x)(3) 通常mk(x) 可表示成多项式.由于地下水初始流场在一般情况下水力坡度以及地形变化较为平缓, 因而在本文中假定平面二维流场, 线性漂移 (l=2) , 则可写成: mk(x)=ak,0+ak,1x+ak,2y(4) 即,fk, 0(x) =1,fk, 1(x) =x,fk, 2(x) =y. 根据无偏性条件:E[H*(x) -H1(x) ]=0, 经推导可得下面的无偏条件组: {n1∑i=1λif1,l(xi)=f1,l(x)n2∑i=1ωjf2,l(xj)=0l=0,1,2(5) 而估计方差为: σ2E=D2[Η*-Η1(x)]=D2[n1∑i=1λiΗ1i+n2∑j=1ωjΗ2j-Η1(x)]=-n1∑i=1n1∑i′=1λiλi′γ(xi,xi′)-n1∑i=1n2∑j=1λjωjγ12(xi,xj)-n2∑j=1n1∑i=1ωjλiγ21(xj,xi)-n2∑j=1n2∑j′=1ωjωj′γ2(xj,xj′)+2n1∑i=1λiγ1(xi,x)+2n2∑j=1ωjγ2(xj,x)(6) 要在无偏条件组下使估计方差达到极小, 可用求条件极值的拉格朗日乘数法, 即要求使: F=σ2E-22∑l=0μl[n1∑i=1λif1,l(xi)-f1,l(x)]-22∑l=0ηl[n2∑j=1ωjf2,l(xj)](7) 达到极小时的诸λi、ωj、μl和ηl(i=1, 2, …,n1;j=1, 2, …,n2;l=0, 1, 2) .上式中μl、ηl(l=0, 1, 2) 为拉格朗日