基于地统计学的年降水量与蒸发量空间变异规律研究.docxVIP

基于地统计学的年降水量与蒸发量空间变异规律研究 0 无偏内插估计方法在水文、法学及环境领域的应用 基于区域变量理论，区域统计是最重要的函数。它以区块变量理论为基础，以矩阵函数为主要工具。它用于研究空间分布、结构、空间相关和依赖的自然现象。它所提供的克立格法是一种最优的无偏内插估计方法。它在水文地质、气象学及环境方面有着广泛的应用。在环境科学领域,刘瑞民等应用地统计学方法对太湖水质参数空间分布特性进行了研究;在气象学中,孟健等应用Kriging的空间分析方法分析了城市大气中SO2浓度的空间分布;在水文地质学中,梁天刚等对甘肃省河西走廊以东地区多年平均降水量进行了空间分布模拟研究;李丽娟等将Kriging法应用于土壤水、径流模数、降水等的估值工作中。本文尝试利用地统计学理论和克立格插值方法对研究区域年降水量和蒸发量的空间变异规律进行分析研究。 1 地统计学研究方法的基本理论 1.1 区域化现象的概念 当一个变量呈现为空间分布时,就称之为区域化变量(Regionalized Variable)。这种变量常常反映某种空间现象的特征,用区域化变量来描述的现象称之为区域化现象。区域化变量,亦称区域化随机变量,G.Matheron(1963)将它定义为以空间点x的三个直角坐标为自变量的随机场。区域化变量具有两个最显著,而且也是最重要的特征,即随机性和结构性。区域化随机变量之间的差异,可以用空间协方差来表示。 1.2 地统计学函数 协方差函数和变异函数是以区域化变量理论为基础建立起来的地统计学的两个最基本的函数。地统计学的主要方法之一,克立格法就是建立在变异函数理论和结构分析基础之上的。 1.2.1 变换函数的估计 设Z(x)是系统某属性Z在空间位置x处的值,Z(x)为一区域化随机变量,并满足二阶平稳假设,h为两样本点空间分隔距离,Z(xi)和Z(xi+h)分别是区域化变量Z(x)在空间位置xi和xi+h处的实测值[i=1,2,…,N(h)],那么,变异函数γ(h)的离散计算公式为: γ(h)=12Ν(h)Ν(h)∑i=1[Ζ(xi)-Ζ(xi+h)]2γ(h)=12N(h)∑i=1N(h)[Z(xi)?Z(xi+h)]2 对不同的空间分隔距离h,计算出相应的协方差函数值c(h)和变异函数值γ(h),如果分别以h为横坐标,c(h)或γ(h)为纵坐标,画出协方差函数和变异函数曲线图,就可以直接展示区域化变量Z(x)的空间变异特点。可见,变异函数能同时描述区域化变量的随机性和结构性,从而在数学上对区域化变量进行严格分析,是空间变异规律分析和空间结构分析的有效工具。 1.2.2 变异函数理论模型 地统计学将变异函数理论模型分为三大类:第一类是有基台值模型,包括球状模型、指数模型、高斯模型、线性有基台值模型和纯块金效应模型;第二类是无基台值模型,包括幂函数模型、线性无基台值模型、抛物线模型;第三类是孔穴效应模型。考虑到本研究的实际需要,下面只简单介绍两种常见的变异函数理论模型:球状模型和指数模型。 球状模型,其一般公式为: γ(h)={0h=0c0+c(3h2a-h32a3)0＜h≤ac0+ch＞aγ(h)=???????0c0+c(3h2a?h32a3)c0+ch=00＜h≤ah＞a 式中:c0为块金(效应)常数,c为拱高,c0+c为基台值,a为变程。当c0=0,c=1时,称为标准球状模型。球状模型是地统计分析中应用最广泛的理论模型,许多区域化变量的理论模型都可以用该模型去拟合。 指数模型,其一般公式为: γ(h)={0h=0c0+c(1-e-ha)h＞0γ(h)={0c0+c(1?e?ha)h=0h＞0 式中:c0和c意义与前相同,但a不是变程。当h=3a时,1-e-ha=1-e-3≈0.95≈11?e?ha=1?e?3≈0.95≈1,即γ(3a)≈c0+c,从而指数模型的变程a′约为3a。当c0=0,c=1时,称为标准指数模型。 1.3 地统计学的空间插值方法,主要有 克立格(Kriging)插值法,又称空间局部估计或空间局部插值法,是根据变异函数模型而发展起来的一系列地统计的空间插值方法,是地统计学的主要内容之一。克立格法适用的条件是:变异函数和相关分析的结果表明区域化变量存在空间相关性。其实质是利用区域化变量的原始数据和变异函数的结构特点,对未采样点的区域化变量的取值进行线性无偏和最优估计。 2 区域化变量及其特征 年降水量和蒸发量,既服从地带性规律,同时又受随机性因素的影响,因此它们是典型的区域化变量。文章以甘肃省53个气象台站多年平均降水量和蒸发量数据为实测值,拟合了年降水量和蒸发量的半变异函数理论模型,并采用普通克立格法做了空间插值计算,结论如下。 2.1 年降水量的空间变异拟合 半变异函数模型,是克立格空间插值的前提条件,同时它也