指示克立格估值计算的理论和方法.docxVIP

指示克立格估值计算的理论和方法 近十年来，随着人口的增加和经济的发展，许多地区的地下水开发利用程度日益严重，导致了一系列区域地质环境问题。地下水数值模拟作为一种常规研究工具,规模和范围也日趋扩大。此类区域性模型应用中的一个突出问题是勘探资料特别是含水层的渗透性参数资料分布不均或缺乏,在通常的地下水数值模拟中只能用简单的参数分区来描述渗透性的非均质性。由于分区过大而忽略分区内部参数的差异性,导致模拟误差的出现。寻求一种尽可能利用有限的勘探资料,对未知区域内的含水层参数进行合理估值的方法,是目前大区域地下水流模拟中的关键性问题。 地质统计学中的克立格方法是线性无偏最优的估值方法。20世纪70年代末之后,克立格方法在水文地质学领域得到了广泛应用。Delhomme (1979 )用克立格法估计法国下诺曼底含水层的导水系数并绘制了导水系数的最优等值线图; Dagan等(1982)用相同方法绘制了意大利威尼斯市地下水头分布的等水头线并以此作为研究地面沉降的依据;Neuman和 Cliftun (1982)用克立格法给出了美国亚利桑拉Avra河谷含水层渗透系数的估计;Ahmed等(1987)也结合实例对使用克立格法求渗透系数与其他几种方法进行了比较。由于自然界大量的数据样本不服从正态分布,用高斯克立格法进行估值比较困难,为此Journel于1983年创立了指示克立格法,它不要求数据总体服从正态分布,且不受特异值的影响。指示克立格法自创立以来,在矿产(Fytas和Chaouai,1990)、土壤(Goovaerts,Webster和Dubois,1997)、土地管理(傅佩红,2004)、石油(Data-Gupta等,1999和Hohn等,1994)等领域有广泛应用。然而,使用指示克立格法对第四纪含水层渗透性空间分布进行估值国内外还鲜有报道。 本文以华北某第四系地层为例,以指示变异函数为结构分析工具,运用指示克立格法对该区域含水层渗透性的空间分布进行三维估计,结果表明指示克立格法能很好地描述第四系含水层渗透性的空间分布规律。 1 拉格朗日函数中的克立格估计 在研究含水层渗透性时,把含水层的渗透系数作为区域化变量,记作Z(u),它同时具有随机性和结构性的特点。建立变量Z(u)的模型时,设渗透系数在区域范围内的变化为二阶平稳过程,即其均值函数和协方差函数存在且平稳,则所有取样点所形成的分布可以代替总体分布,而区域内任一点处变量取值的概率分布服从这一总体。作所有取样点的累积频率曲线,将累积频率曲线由小到大截成K+1段,则有K个截断值(阈值)Zk,k=1,2,…,K,通常以其十分位数作为阈值(如果是类型变量,则以其类型作为阈值)。为求出位置u处的渗透系数Zu小于等于阈值Zk的概率,引进二值指示函数,其定义为 Ι(u?Ζk)={1当位置u处的渗透系数Ζu小于等于阈值Ζk时0其他(1)I(u?Zk)=?????1当位置u处的渗透系数Zu小于等于阈值Zk时0其他(1) 指示变异函数可描述为 γΙ(h)=1Ν(h)Ν(h)∑α=1[i(uα+h,Ζk)-i(uα,Ζk)]2(2)γI(h)=1N(h)∑α=1N(h)[i(uα+h,Zk)?i(uα,Zk)]2(2) 它与协方差函数有如下关系: γI(h)=CI(0)-CI(h),CI(0) =Var[I(u,z)]为方差。 在u处出现渗透系数Zu小于等于阈值Zk的数学期望是: E{Ι(u,Ζk)}=1?Ρrob{Ζ(u)≤Ζk}+0?Ρrob{Ζ(u)Ζk}=Ρrob{Ζ(u)≤Ζk}=F(Ζk)(3)E{I(u,Zk)}=1?Prob{Z(u)≤Zk}+0?Prob{Z(u)Zk}=Prob{Z(u)≤Zk}=F(Zk)(3) 上式表明在u处出现渗透系数Zu小于等于阈值Zk的数学期望等于它出现的概率。所以,指示克立格法估计量能够用于估计特征出现的概率,这个估计值可通过实测渗透系数的线性函数表示: F(u,Ζk)*=E{Ι(u,Ζk)ㄧ(n)}*=Ρrob*{Ζ(u)≤Ζkㄧ(n)}=n∑α=1λα(u,Ζk)Ι(uα,Ζk)(4)F(u,Zk)?=E{I(u,Zk)ㄧ(n)}?=Prob?{Z(u)≤Zkㄧ(n)}=∑α=1nλα(u,Zk)I(uα,Zk)(4) 式中F(u,Zk)*为渗透系数在位置u处真实概率F(u,Zk)的克立格估计;uα是第α个实测渗透系数的位置,α=1,2,…,n;λα(u,Zk)是相应的截断值Zk的简单克立格法权重。要确定每个权,必须使估计方差达到最小,即: min{E[F(u,Ζk)*-F(u,Ζk)]2}min{E[F(u,Zk)??F(u,Zk)]2} 且受无偏性条件n∑α=1λα(u,Ζk)=1的限制。构造拉格朗日函数L (λ1,…,λn,μ),μ为拉格朗日乘数。用L对