第九讲地统计分析方法.ppt

《第九讲地统计分析方法.ppt》由会员分享，可在线阅读，更多相关《第九讲地统计分析方法.ppt（54页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第九讲地统计分析方法现在学习的是第1页，共54页2基础知识基础知识地统计（Geostatistics）又称地质统计，是在法国著名统计学家G.Matheron大量理论研究的基础上逐渐形成的一门新的统计学分支。它是以区域化变量为基础，借助变变异异函函数数，研究既具有随机性又具有结构性，或空间相关性和依赖性的自然现象的一门科学。凡是与空间数据的结构性和随机性，或空间相关性和依赖性，或空间格局与变异有关的研究，并对这些数据进行最优无偏内插估计，或模拟这些数据的离散性、波动性时，皆可应用地统计学的理论与方法。地统计分析的核心就是通过对采样数据的分析、对采样区地理特征的认识选择合适的空间内插方法创建表面。

2、现在学习的是第2页，共54页3n前提假设 n随机过程随机过程 地统计学认为研究区域中的所有样本值都是随机过程地统计学认为研究区域中的所有样本值都是随机过程的结果，即所有样本值都不是相互独立的，它们的结果，即所有样本值都不是相互独立的，它们是遵循一定的内在规律的。因此地统计学就是要是遵循一定的内在规律的。因此地统计学就是要揭示这种内在规律，并进行预测。揭示这种内在规律，并进行预测。n正态分布正态分布 若不符合正态分布的假设，应对数据进行变换，转若不符合正态分布的假设，应对数据进行变换，转为符合正态分布的形式，并尽量选取可逆的变换为符合正态分布的形式，并尽量选取可逆的变换形式。形式。现在学习的是第

3、3页，共54页4n平稳性平稳性 包括两种平稳性：一类是包括两种平稳性：一类是均值平稳均值平稳；另一类是与协方差；另一类是与协方差函数有关的函数有关的二阶平稳二阶平稳和与半变异函数有关的和与半变异函数有关的内蕴平稳内蕴平稳。均值平稳均值平稳，即假设均值是不变的并且与位置无关；，即假设均值是不变的并且与位置无关；二阶平稳二阶平稳是假设具有相同的距离和方向的任意两点的是假设具有相同的距离和方向的任意两点的协方差是相同的，协方差只与这两点的值相关而协方差是相同的，协方差只与这两点的值相关而与它们的位置无关；与它们的位置无关；内蕴平稳内蕴平稳假设是指具有相同距离和方向的任意两点的假设是指具有相同距离和方

4、向的任意两点的方差（即变异函数）是相同的。二阶平稳和内蕴方差（即变异函数）是相同的。二阶平稳和内蕴平稳都是为了获得基本重复规律而作的基本假设，平稳都是为了获得基本重复规律而作的基本假设，通过协方差函数和变异函数可以进行预测和估计通过协方差函数和变异函数可以进行预测和估计预测结果的不确定性。预测结果的不确定性。现在学习的是第4页，共54页5n区域化变量n当一个变量呈现一定的当一个变量呈现一定的空间分布空间分布时，称之为区时，称之为区域化变量，它反映了区域内的某种特征或现象。域化变量，它反映了区域内的某种特征或现象。n区域化变量与一般的随机变量不同之处在于，区域化变量与一般的随机变量不同之处在于，

5、一般的随机变量取值符合一定的概率分布，而一般的随机变量取值符合一定的概率分布，而区域化变量根据区域内位置的不同而取不同的区域化变量根据区域内位置的不同而取不同的值。而当区域化变量在区域内确定位置取值时，值。而当区域化变量在区域内确定位置取值时，表现为一般的随机变量，也就是说，它是与位表现为一般的随机变量，也就是说，它是与位置有关的随机变量。置有关的随机变量。n区域化变量具有两个显著特征：即区域化变量具有两个显著特征：即随机性和结随机性和结构性构性。现在学习的是第5页，共54页 地统计学是以区域化变量理论为基础，以变异函数为主要工具，研究那些在空间分布上既有随机性又有结构性，或空间相关和依赖性的

6、自然现象的科学。协方差函数和变异函数是以区域化变量理论为基础建立起来的地统计学的两个最基本的函数。地统计学的主要方法之一，克立格法就是建立在变异函数理论和结构分析基础之上的。现在学习的是第6页，共54页（一）协方差函数 n协方差函数的概念 区域化随机变量之间的差异，可以用空间协方差来表示。在概率论中,随机向量X与Y的协方差被定义为 E(.)为期望值.区域化变量在空间点x和x+h处的两个随机变量和的二阶混合中心矩定义为Z(x)的自协方差函数，即（4.2.24.2.2）(4.2.1)(4.2.1)现在学习的是第7页，共54页n协方差函数的计算公式 式中：h为两样本点空间分隔距离或距离滞后；为 在空

7、间位置 处的实测值；是 在 处距离偏离h的实测值i=1，2，是分隔距离为h时的样本点对（paris）总数，和 分别为 和 的样本平均数,即(4.2.3)(4.2.3)(4.2.4)(4.2.4)(4.2.5)(4.2.5)现在学习的是第8页，共54页 若 =m（常数），则上式可以改写为 式中：m为样本平均数，可由一般算术平均数公式求得，即(4.2.6(4.2.6)现在学习的是第9页，共54页（二）变异函数 n变异函数的概念 变异函数variograms），又称变差函数、变异矩，是地统计分析所特有的基本工具。在一维条件下变异函数定义为，当空间点x在一维x轴上变化时，区域化变量Z(x)在点x和x+

8、h处的值Z(x)与Z(x+h)差的方差的一半为区域化变量Z(x)在x轴方向上的变异函数，记为(h)，即 (4.2.7(4.2.7)方差等于平方均值减去均值的平方方差等于平方均值减去均值的平方现在学习的是第10页，共54页 在二阶平稳假设条件下，对任意的h有 因此，公式可以改写为 从上式可知，变异函数依赖于两个自变量 x和h，当变异函数 仅仅依赖于距离h而与位置x无关时，可改写成 ，即(4.2.9)(4.2.9)(4.2.8)(4.2.8)现在学习的是第11页，共54页n变异函数的性质 设Z(x)是区域化变量，在满足二阶平稳假设条件下，变异函数式具有如下性质：(1)=0，即在h=0处，变异函数为

9、0；(2)=，即 关于直线h=0是对称的，它是一个偶函数；(3)0，即 只能大于或等于0；现在学习的是第12页，共54页n变异函数的计算公式 设 是系统某属性Z在空间位置x处的值，为一区域化随机变量，并满足二阶平稳假设，h为两样本点空间分隔距离，和 分别是区域化变量 在空间位置 和 处的实测值i=1,2,N(h)，那么，变异函数 的离散计算公式为(4.2.10)(4.2.10)现在学习的是第13页，共54页 这样对不同的空间分隔距离h，计算出相应的 和 值。如果分别以h为横坐标，或 为纵坐标，画出协方差函数和变异函数曲线图，就可以直接展示区域化变量Z(x)的空间变异特点。可见，变异函数能同时描

10、述区域化变量的随机性和结构性，从而在数学上对区域化变量进行严格分析，是空间变异规律分析和空间结构分析的有效工具。现在学习的是第14页，共54页例如：假设某地区降水量Z(x)（单位：mm）是二维区域化随机变量，满足二阶平稳假设，其观测值的空间正方形网格数据如图4.2.1所示（点与点之间的距离为h=1 km）。试计算其南北方向及西北和东南方向的变异函数。现在学习的是第15页，共54页图4.2.1 空间正方形网格数据（点间距h=1 km）从图4.2.1可以看出，空间上有些点，由于某种原因没有采集到。如果没有缺失值，可直接对正方形网格数据结构计算变异函数；在有缺失值的情况下，也可以计算变异函数。只要“

11、跳过”缺失点位置即可（图4.2.2）。现在学习的是第16页，共54页 首先计算南北方向上的变异函数值，由变异函数的计算公式可得 =385/72=5.35 图4.2.2 缺失值情况下样本数对的组成和计算过程为缺失值为缺失值 现在学习的是第17页，共54页 同样计算出 最后，得到南北方向和西北东南方向上的变异函数计算结果见下表。同样可以计算东西方向上的变异函数。方向 南北 方向 西北东南 h12345h1.412.824.245.657.07N(h)36 27 21 13 5 N(h)322113825.359.2617.5525.6922.907.0612.9530.8558.1350.00现在

12、学习的是第18页，共54页n变异函数的参数 变异函数有个非常重要的参数，即基台值（sill）、变程（range）或称空间依赖范围（range of spatial dependence）、块金值（nugget）或称区域不连续性值（localized discontinuity）和分维数（fractal dimension）。前3个参数可以直接从变异函数图中得到。它们决定变异函数的形状与结构。变异函数的形状反映自然现象空间分布结构或空间相关的类型，同时还能给出这种空间相关的范围。现在学习的是第19页，共54页 当变异函数随着间隔距离h的增大，从非零值达到一个相对稳定的常数时，该常数称为基台值C0

13、+C。当间隔距离h=0时，(0)=C0，该值称为块金值或块金方差（n nugget variance）。基台值是系统或系统属性中最大的变异，变异函数达到基台值时的间隔距离a称为变程。变程表示在ha以后，区域化变量Z(x)空间相关性消失。块金值表示区域化变量在小于抽样尺度时非连续变异，由区域化变量的属性或测量误差决定。现在学习的是第20页，共54页 上述个参数可从变异函数曲线图直接得到，或通过估计曲线回归参数得到。第4个参数，即分维数用于表示变异函数的特性，由变异函数 和间隔距离h之间的关系确定 分维数D为双对数直线回归方程中的斜率，它是一个无量纲数。分维数D的大小，表示变异函数曲线的曲率，可以

14、作为随机变异的量度。但该随机分维数D与形状分维数有本质的不同。现在学习的是第21页，共54页n变异函数的理论模型 地统计学将变异函数理论模型分为3大类：第1类是有基台值模型，包括球状模型、指数模型、高斯模型、线性有基台值模型和纯块金效应模型；第2类是无基台值模型，包括幂函数模型、线性无基台值模型、抛物线模型；第3类是孔穴效应模型。下面有代表性地介绍几种常见的变异函数理论模型。现在学习的是第22页，共54页 纯块金效应模型:其一般公式为 式中：c00，为先验方差。该模型相当于区域化变量为随机分布，样本点间的协方差函数对于所有距离h均等于0，变量的空间相关不存在。(4.2.11)(4.2.11)现

15、在学习的是第23页，共54页 球状模型:其一般公式为 式中：c0为块金（效应）常数;c为拱高;c0+c为基台值;a为变程。当c0=0，c=1时，称为标准球状模型。球状模型是地统计分析中应用最广泛的理论模型，许多区域化变量的理论模型都可以用该模型去拟合。(4.2.12)(4.2.12)现在学习的是第24页，共54页 指数模型:其一般公式为 式中：c0和c意义与前相同，但a不是变程。当h=3时，即 ，从而指数模型的变程 约为 。当c0=0，c=1时，称为标准指数模型。(4.2.13(4.2.13)现在学习的是第25页，共54页 高斯模型:其一般公式为 式中：c0和c意义与前相同，a也不是变程。当

16、时，即 ，因此高斯模型的变程 约为 。当 时，称为标准高斯函数模型。(4.2.14)(4.2.14)现在学习的是第26页，共54页 幂函数模型:其一般公式为 式中：为幂指数。当变化时，这种模型可以反映在原点附近的各种性状。但是必须小于2，若 ，则函数 就不再是一个条件非负定函数了，也就是说它已经不能成为变异函数了。(4.2.15)(4.2.15)现在学习的是第27页，共54页 对数模型:其一般公式为 显然，当 ，这与变异函数的性质 不符 。因此，对数模型不能描述点支撑上的区域化变量的结构。(4.2.16)(4.2.16)现在学习的是第28页，共54页线性有基台值模型:其一般公式为 式中:该模型

17、的变程为a，基台值为 。线性无基台值模型:其一般公式为 从式中可以看出，该模型没有基台值，也没有变程。(4.2.18)(4.2.18)(4.2.17)(4.2.17)现在学习的是第29页，共54页例如:某地区降水量是一个区域化变量，其变异函数 的实测值及距离h的关系见下表，下面我们试用回归分析方法建立其球状变异函数模型。实测值(h)距离h实测值(h)距离h2.10.69.24.94.31.110.35.15.72.210.56.26.52.510.97.57.83.111.29.58.83.812.49.8现在学习的是第30页，共54页 从上面的介绍和讨论，我们知道，球状变异函数的一般形式为

18、当 时，有现在学习的是第31页，共54页 如果记 ，则可以得到线性模型 根据表中的数据，对上式进行最小二乘拟合，得到 (4.2.20)计算可知，上式的显著性检验参数F=114.054，R2=0.962，可见模型的拟合效果是很好的。(4.2.19)现在学习的是第32页，共54页 比较(4.2.20)式与(4.2.19)式，并做简单计算可知：c0=2.048，c=1.154，a=8.353，所以，球状变异函数模型为(4.2.21)现在学习的是第33页，共54页 (三)克立格插值方法 克立格（Kriging）插值法，又称空间局部估计或空间局部插值法，是地统计学的主要内容之一。克立格法是建立在变异函数

19、理论及结构分析基础之上的，它是在有限区域内对区域化变量的取值进行无偏最优估计的一种方法。克立格法适用的条件是，如果变异函数和相关分析的结果表明区域化变量存在空间相关性。其实质是利用区域化变量的原始数据和变异函数的结构特点，对未采样点的区域化变量的取值进行线性无偏、最优估计。现在学习的是第34页，共54页 克立格插值（riging interpolation)是根据变异函数模型而发展起来的一系列地统计的空间插值方法，包括：普通克立格法（ordinary riging）;泛克立格法（universal riging）;指示克立格法（indicator riging）;析取克立格法（disjunct

20、ive riging）;协同克立格法（cokriging）等。下面仅对普通克立格法作一些简单介绍。现在学习的是第35页，共54页 首先假设区域化变量 满足二阶平稳假设和本征假设，其数学期望为m，协方差函数 及变异函数 存在。即 假设在待估计点（x）的临域内共有n个实测点，即x1，x2，xn，其样本值为 。那么，普通克里格法的插值公式为(4.2.22)(4.2.22)现在学习的是第36页，共54页 其中 为权重系数，表示各空间样本点 处的观测值 对估计值 的贡献程度。可见，克立格插值的关键就是计算权重系数 。显然，权重系数的求取必须满足两个条件：一是使 的估计是无偏的，即偏差的数学期望为零；二是

21、最优的，即使估计值 和实际值 之差的平方和最小。为此，需要满足以下两个条件：现在学习的是第37页，共54页 (1)无偏性。要使 成为 的无偏估计量，即 当 时，也就是当 时，则有 这时，为 的无偏估计量。（2）最优性。在满足无偏性条件下，估计方差为(4.2.23(4.2.23)现在学习的是第38页，共54页 使用协方差函数表达，它可以进一步写为 (4.2.24)为使估计方差最小，根据拉格朗日乘数原理，令 (4.2.25)求F对 和 的偏导数，并令其为0，得克立格方程组 (4.2.26)现在学习的是第39页，共54页 (4.2.27)(4.2.28)整理后得 解线性方程组(4.2.27)式,求出

22、权重系数i和拉格朗日系数,代入公式(4.2.24),可得克立格估计方差现在学习的是第40页，共54页 在变异函数存在的条件下，根据协方差与变异函数的关系：或 ，也可以用变异函数表示普通克立格方程组和克立格估计方差，即 (4.2.29)解线性方程组(4.2.27)式，求出权重系数 和拉格朗日乘数，代入公式(4.2.24)，可得克立格估计方差 ，即(4.2.30)现在学习的是第41页，共54页上述过程也可用矩阵形式表示，令 则普通克立格方程组为 (4.2.31)解方程组(4.2.31)式，可得 (4.2.32)其估计方差为 (4.2.33)现在学习的是第42页，共54页 也可以将克立格方程组和估计

23、方差用变异函数写成上述矩阵形式。令 在以上的介绍中，区域化变量 的数学期望 可以是已知或未知的。如果m是已知常数，称为简单克立格法；如果m是未知常数，称为普通克立格法。不管是哪一种方法，均可根据方法计算权重系数和克立格估计量。(4.2.34)(4.2.35)(4.2.36)现在学习的是第43页，共54页 以图4.2.1为例，4个观测点x1，x2，x3，x4的观测值分别为Z(x1)=37、Z(x2)=42、Z(x3)=36、Z(x4)=35，如果假设降水量的变异函数是向同性（即变异函数在各个方向的变化都相同）的二维球状模型，其具体形式为（4.2.21）式。现在，我们用普通克立格法估计观测点x0的

24、降水量值Z(x0)。根据普通克立格法的基本原理，我们知道，Z(x0)估计的基本公式应该是 现在学习的是第44页，共54页 根据公式(4.2.32)，可知 (4.2.37)根据协方差与变异函数的关系以及（4.2.21）式，可得协方差函数 现在学习的是第45页，共54页 当 时，根据克立格矩阵的对称性，当 时，由此计算可得 现在学习的是第46页，共54页现在学习的是第47页，共54页 将以上计算结果代入克立格方程组（4.2.31），得 现在学习的是第48页，共54页 即 克 立 格 权 重 系 数 分 别 为：1=0.287，2=0.210,3=0.202,4=0.301，=-0.473，所以观测

25、点的降水量的克立格估计值为：根据普通克立格法的基本原理，我们知道，Z(x0)估计的基本公式应该是 37.2537.25（mmmm）。克立格估计方差为 现在学习的是第49页，共54页二、应用实例 年降水量和蒸发量，既服从地带性规律，同时又受随机性因素的影响，因此它们是典型的区域化变量。我们以甘肃省53个气象台站多年平均降水量和蒸发量数据（见教材表3.1.2）为实测值，拟合了年降水量和蒸发量的半变异函数理论模型，并采用普通克立格法和双变量协同克里格法，做了空间插值计算，结论如下。现在学习的是第50页，共54页 （一）半变异函数 半变异函数模型，是克立格空间插值的前提条件，同时它也决定着空间插值的精

26、度。一般情况下，半变异函数模型是根据半变异函数云图的分布，选择合适的理论模型，按照估计方差最小的原则，运用最小二乘法求得。图4.2.4和图4.2.5分别给出了年降水量和年蒸发量的半变异函数云图。图4.2.4 年降水量的半变异函数云图 现在学习的是第51页，共54页 图4.2.5 年蒸发量的半变异函数云图 从图4.2.4和图4.2.5可以看出，年降水量和年蒸发量的块金效应都不明显，这是因为样本点是各个气象站点的实测值，空间分辨率可以忽略不计，另外实验误差和人为误差基本上都很小。我们选择各种不同的半变异函数理论模型，经过多次拟合计算和对比分析，发现指数模型比较好地描述了年降水量的空间变异规律。其变

27、异函数的具体形式如下：现在学习的是第52页，共54页 (4.2.38)(4.2.38)式拟合的适度系数为 。我们选择各种不同的半变异函数理论模型，经过多次拟合计算和对比分析，发现球状模型比较好地描述了年蒸发量的空间变异规律。其变异函数的具体形式如下 (4.2.39)(4.2.39)式拟合的适度系数为 。现在学习的是第53页，共54页 （二）空间插值结果 基 于 半 变 异 函 数 的 理 论 模 型(4.2.38)和(4.2.39)，对甘肃省范围内的年降水量和蒸发量，用普通克立格法进行空间插值计算，得到的结果分别如图4.2.4和图4.2.5。（三）结果讨论 从图4.2.6可以看出，在甘肃省范围内，年降水量的空间分布格局总体上是东南多西北少，并且呈现从东南方向到西北方向逐渐过渡，梯度变化明显；山地多，平地少，南北方向从南部祁连山脉向北部的沙漠戈壁逐渐减少。现在学习的是第54页，共54页