变异函数及结构分析课件

——地统计学的工具第四章变异函数及结构分析——地统计学的工具第四章变异函数及结构分析1一、协方差函数的计算公式第一节协方差函数和变异函数的性质设区域化变量Z(x)满足(准)二阶平稳假设，h为两样本点空间分隔距离，Z(xi)与Z(xi+h)分别是Z(x)在空间位置xi和xi+h上的观测值(i=1,2,…N(h)),则计算协方差的公式为：一、协方差函数的计算公式第一节协方差函数和变异函数的性质2变异函数及结构分析课件3协方差函数曲线图：以h为横坐标，C#(h)为纵坐标作图协方差函数曲线图：以h为横坐标，C#(h)为纵坐标作图4二、协方差函数的性质区域化变量Z(x)在二阶平稳假设下，其协方差函数存在且平稳，定义为1.C(0)=Var[Z(x)]≥0，即先验方差不能小于零2.C(h)=C(-h)，即C(h)对h=0的直线对称，是一个偶函数二、协方差函数的性质区域化变量Z(x)在二阶平稳假设下，其协5证：令x-h=y，则x=y+h，带入上式得图形特征及含义证：令x-h=y，则x=y+h，带入上式得图形特征及含义63.|C(h)|≤C(0)，即协方差函数绝对值小于等于先验方差证：3.|C(h)|≤C(0)，即协方差函数绝对值小于等于先74.|h|→∞时,C(h)→0，或写作C(∞)=0，即当空间距离很大时，协方差函数值很小意义(空间局限性)：当距离很大时，Z(x)和Z(x+h)之间的线性相关基本不存在5.C(h)必须是一个非负定函数，由C(xi-xj）构成的协方差函数矩阵必须是非负定矩阵正定条件(positivedefinitecondition)区域化变量Z(x)二阶平稳，其数学期望为m，协方差为C(h),变异函数为γ(h)，令Y是该类型区域化变量的任意有限线性组合，即：4.|h|→∞时,C(h)→0，或写作C(∞)=0，即当8较难理解则由C(xi-xj）(i,j=1,2…n)构成的协方差函数矩阵是非负定矩阵，即C(h)为非负定函数二阶平稳区域化变量的协方差函数是有条件的较难理解则由C(xi-xj）(i,j=1,2…n)构成的协方9三、实验(经验)变异函数(experimentalvariogram)的计算公式设区域化变量Z(x)满足(准)二阶平稳条件或(准)本征假设，h为两样本点空间分隔距离，Z(xi)与Z(xi+h)分别是Z(x)在空间位置xi和xi+h上的观测值(i=1,2,…N(h)),则计算实验变异函数的公式为：变异函数曲线图：以h为横坐标，γ#(h)为纵坐标作图三、实验(经验)变异函数(experimentalvari10变异函数计算实例（1）一维变异函数的计算x1x2x3x4x5x6x7x8x9x104

3

4

5797877以下为一研究对象在水平方向上的采样数据，满足二阶平稳或本征假设，采样值如图所示，点间分隔距离h=1米，计算γ#(h)变异函数计算实例（1）一维变异函数的计算x1x211变异函数及结构分析课件12变异函数及结构分析课件13两方面理解：变异性的理解与相关性的理解两方面理解：变异性的理解与相关性的理解14作业：x1x2x3x4x5x6x7x8

2

4

3

15364以下为一研究对象在水平方向上的采样数据，满足二阶平稳或本征假设，采样值如图所示，点间分隔距离h=1米，计算γ#(1),γ#(2),γ#(3)作业：x1x2x3x4x515（2）二维变异函数的计算下图为正方形网格状的采样数据，*号处为无数据点，点间距离h为100米，请分别计算南北、东西、西北和东南方向上的变异函数值。（2）二维变异函数的计算下图为正方形网格状的采样数据，*号处16变异函数及结构分析课件17变异函数及结构分析课件18变异函数及结构分析课件19西北和东南方向上的变异函数值的计算，注意分隔距离h的确定和样本数据对的查找西北和东南方向上的变异函数值的计算，注意分隔距离h的确定和样20变异函数及结构分析课件21作业：下图为正方形网格状的采样数据，网格交叉空白处为无数据点，点间距离h为a米，请分别计算南北方向γ#(a),

西北—东南方向上γ#(a)。作业：下图为正方形网格状的采样数据，网格交叉空白处为无数据点22四、变异函数的性质区域化变量Z(x)满足二阶平稳或本征假设条件，则变异函数存在且平稳，计算公式为1.γ(0)=0，即在h=0时，变异函数为零2.γ(h)=γ(-h)，即γ(h)对h=0的直线对称，是一个偶函数四、变异函数的性质区域化变量Z(x)满足二阶平稳或本征假设条233.γ(h)≥0,即研究现象的变异性只能大于或等于零4.|h|→∞时,γ(h)→C(0)，或写作γ(∞)=C(0)，即当空间距离很大时，变异函数值接近先验方差5.[-γ(h)]必须是一个条件非负定函数，即由[-γ(xi-xj)]构成的变异函数矩阵必须是条件非负定矩阵。3.γ(h)≥0,即研究现象的变异性只能大于或等于零4.24区域化变量Z(x)二阶平稳，其数学期望为m，协方差为C(h),变异函数为γ(h)，令Y是该类型区域化变量的任意有限线性组合，即：区域化变量Z(x)二阶平稳，其数学期望为m，协方差为C(h)25区域化变量Z(x)的变异函数γ(h)是有条件的，即需满足条件非负定条件区域化变量Z(x)的变异函数γ(h)是有条件的，即需满足条26五、协方差函数与变异函数的关系五、协方差函数与变异函数的关系27协方差函数和变异函数的曲线图问题：为什么只画出了h>0的关系图？当h足够大(即存在a>0，当h≥a)时，可以使C(h)=0,γ(h)=C(0),a称为变程(range)1、变程a表示区域化变量从存在空间相关状态(当|h|<a时）转向不存在空间相关状态(当|h|>a时)的转折点2、变程a的大小反映区域化变量影响范围的大小，或说反映该变量自相关范围的大小。也可说变程a是区域化变量空间变异尺度或空间自相关尺度变程a的意义：协方差函数和变异函数的曲线图问题：为什么只画出了h>0的关系28第二节变异函数的功能一、变异函数通过“变程”反映变量的影响范围——变异函数的跃迁现象变异函数γ(h)是一个单调递增函数，当h超过某一数值(变程a)后，γ(h)不再继续单调地增大，而往往稳定在一个极限值γ(∞)附近，这种现象称为“跃迁现象”(transitionphenomena)γ(∞)极限值称为基台值(sill),即C(0)【二阶平稳条件】，基台值的大小反映变量变化幅度的大小凡具有一个变程a和一个基台值的变异函数，称为“跃迁型”的变异函数“变程”反映变量的影响范围(图示）第二节变异函数的功能一、变异函数通过“变程”反映变量的影29二、不同方向上的变异函数图可反映区域化变量的各向异性——变异函数表示的各向异性如果在各个方向上区域化变量的变异性相同或相近，则称区域化变量是各向同性的，反之称为各向异性通过作出各个方向上的变异函数图，并放到一起来比较、分析、研究，就可以确定区域化变量的各向异性（包括有无各向异性，及各向异性的类型等）二、不同方向上的变异函数图可反映区域化变量的各向异性30变异函数及结构分析课件31三、块金常数C0的大小可反映区域化变量的随机性大小——变异函数的块金效应当h=0时，变异函数γ(h)≠0，而等于一个常数C0，这种现象称为“块金效应”(nuggeteffect),C0称为块金常数或块金方差(nuggetvariance)块金效应的图形表示三、块金常数C0的大小可反映区域化变量的随机性大小32“块金效应”

主要有两种来源：1、区域化变量在小于抽样尺度h时所具有的变异性2、采样分析误差当样点间的距离大于微域结构的范围，或样点的大小大于微域结构的范围就会出现块金效应（Webester，1985）“块金效应”主要有两种来源：1、区域化变量在小于抽样尺度h33四、变异函数在原点处的性状可反映区域化变量的空间连续性变异函数在原点处的性状主要有五种类型，每种类型反映了变量的不同程度的空间连续性

1、抛物线型（parabolictype）当|h|→0时,γ(h)→A|h|2(A为常数)，即变异函数曲线在原点处趋向一条抛物线，反映区域化变量是具有高度连续性的，如矿层厚度四、变异函数在原点处的性状可反映区域化变量的空间连续性变34

2、线性型（lineartype）当|h|→0时,γ(h)→A|h|(A为常数)，即变异函数曲线在原点处趋向一条直线，或说在原点处有斜向的切线存在，反映区域化变量是具有平均的连续性，如金属品位

3、间断型（discontinuoustype）当|h|→0时,γ(h)→C0，即变异函数曲线在原点处间断，说明块金效应存在，又称“块金效应型”，反映区域化变量的连续性很差，但当h增大时，γ(h)又变的较为连续了，如金品位2、线性型（lineartype）当|h|→0时,γ35

4、随机型（randomtype）这种变异函数可看成具有基台值C0和无穷小变程a的跃迁型变异函数，则无论h多小，h总大于a，故Z(x)与Z(x+h)总是互不相关又称纯块金效应型，反映了区域化变量完全不存在空间相关的情况，则本质上此区域化变量为普通随机变量此时，C0=C(0)4、随机型（randomtype）这种变异函数可看成具36

5、过渡型：介于抛物线型和随机型间当|h|→0时,γ(h)→C0，即有块金效应；当|h|=a时,γ(a)=C(0)，即有基台值(C0+C)和变程a,C称为“拱高”过渡型是实际研究工作中最常遇到的一种类型5、过渡型：介于抛物线型和随机型间当|h|→0时,γ(37第三节变异函数的理论模型

思考：是否有了采样数据及变异函数计算公式就可以获知任意距离h的区域化变量变异性？设Z(x)具有各向同性的变异函数

γ(h)，则常见的变异函数模型如下：变异函数的理论模型有基台值模型无基台值模型可以有或无基台值模型：孔穴效应模型球状模型、指数模型高斯模型线性有基台模型纯块金效应模型幂函数模型对数模型线性无基台模型第三节变异函数的理论模型思考：是否有了采样数据及变异函38一、有基台值模型

1、球状模型（sphericalmodel）若模型满足二阶平稳假设，且有有限先验方差，γ(h)值随h的变大而增大，当h达一定值(h>a)时，γ(h)达到一定值——基台值，则称此类模型为有基台值模型式中：C0为块金常数，(C0+C)为基台值，C为拱高，a为变程一、有基台值模型1、球状模型（sphericalmode39当C0=0，C=1，称为标准球状模型，其图形为：原点处切线的斜率为3/2a,与基台值线交点的横坐标为2a/3球状模型是地统计学应用最广的理论模型，许多区域化变量的理论模型都可以用球状模型来拟合当C0=0，C=1，称为标准球状模型，其图形为：原点处切线的40

2、指数模型（exponentialmodel）式中：C0，C意义同前，但a不是变程当C0=0，C=1，称为标准指数模型，其图形为：由于1-e-3=1-0.05=0.95≈1，则变程为3a2、指数模型（exponentialmodel）式中：41

3、高斯模型（gaussianmodel）式中：C0，C意义同前，但a不是变程由于1-e-3=1-0.05=0.95≈1，则变程为√3

a当C0=0，C=1，称为标准高斯函数模型，其图形为：3、高斯模型（gaussianmodel）式中：C0，42模型通过原点切线与基台值线交点的横坐标变程原点处的性状球状2a/3a直线指数a3a直线高斯无交点√3

a抛物线三种模型的比较模型通过原点切线与基台值线交点的横坐标变程原点处的性状球状243

4、线性有基台值模型（linearwithsillmodel）式中：C0，C意义同前，A为常数，表示直线的斜率，变程为a4、线性有基台值模型（linearwithsillm44

5、纯块金效应模型（purenuggeteffectmodel）此时，C0=C(0)此种模型意味着区域化变量为随机分布，样点间的协方差函数对于所有距离h均等于0，即变量不存在空间相关性5、纯块金效应模型（purenuggeteffect45变异函数及结构分析课件46二、无基台值模型

1、幂函数模型（powermodel）若与模型相应的区域化变量不满足二阶平稳假设，仅满足本征假设，γ(h)值随h的变大而增大，但不能达到一定值，即无基台值，则称此类模型为无基台值模型当改变参数θ时，可以表示原点处的各种性状二、无基台值模型1、幂函数模型（powermodel）47

2、线性无基台值模型（linearwithoutsillmodel）2、线性无基台值模型（linearwithoutsil48

3、对数模型（powermodel）3、对数模型（powermodel）49三、孔穴效应模型(holeeffectmodel)当变异函数γ(h)在大于一定距离后，并非单调递增，而具有一定周期波动，此种模型称为孔穴效应模型有基台值无基台值三、孔穴效应模型(holeeffectmodel)当变50第四节变异函数的结构分析一、结构分析、套合结构概念采样数据计算γ#(h)试验变异函数曲线对区域化变量进行分析合适的理论模型实际中区域化变量的变化性很复杂：(1)可能在不同方向上有不同的变异性；(2)在同一方向上包含不同尺度上的多层次的变异性这么复杂！？第四节变异函数的结构分析一、结构分析、套合结构概念采样数51矿床或矿体的变异性往往由多种原因引起采样、样品制备及分析等过程所产生的误差原因矿物成分的变化，如金矿等品位变化剧烈的矿床上尤为明显矿层与夹层的交替变化矿床分布引起的变异≈01～ncm米至百米公里尺度不同原因引起的变异特性，其变异尺度的大小不同显然，大尺度的变异总是包含着小尺度的变异，小尺度的变异在大尺度变异曲线上只能作为“块金效应”出现矿床或矿体的变异性往往由多种原因引起采样、样品制备及分析等过52土壤的空间变异性与土壤母质、气候、水文、地形和生物等因素相关Z(x)Z(x+h）hh≈0h≤1mh≤100m取样和测定误差+其它因素，如水分+地形影响合适的理论模型！？结构分析土壤的空间变异性与土壤母质、气候、水文、地形和生物等因素相关53结构分析:就是构造一个变异函数模型对于全部有效结构信息作定量化的概括，以表征区域化变量的主要特征。结构分析的主要方法:套合结构套合结构(nestedstructure):就是把分别出现在不同距离h上和(或）不同方向α上同时起作用的变异性组合起来套合结构表达式:套和结构可以表示为多个变异函数之和，每一个变异函数代表一种特定尺度上的变异性，表达式为：γi(h)可以是相同或不同的理论模型结构分析:就是构造一个变异函数模型对于全部有效结构信息作定量54二、单一方向上的套合

设区域化变量Z(x)在某一方向上的变异性是由γ0(h)、γ1(h)和γ2(h)组成微观尺度，纯块金效应模型变程为a1的球状模型变程为a2的球状模型(a2>a1)二、单一方向上的套合设区域化变量Z(x)在某一方向上的55变异函数及结构分析课件56变异函数及结构分析课件57作业：

设区域化变量Z(x)在某一方向上的变异性是由γ0(h)、γ1(h)和γ2(h)组成，请求出其套合结构表达式并写出计算过程微观尺度，纯块金效应模型变程为a的球状模型变程为3a的指数模型作业：设区域化变量Z(x)在某一方向上的变异性是由γ058变异函数及结构分析课件59——地统计学的工具第四章变异函数及结构分析——地统计学的工具第四章变异函数及结构分析60一、协方差函数的计算公式第一节协方差函数和变异函数的性质设区域化变量Z(x)满足(准)二阶平稳假设，h为两样本点空间分隔距离，Z(xi)与Z(xi+h)分别是Z(x)在空间位置xi和xi+h上的观测值(i=1,2,…N(h)),则计算协方差的公式为：一、协方差函数的计算公式第一节协方差函数和变异函数的性质61变异函数及结构分析课件62协方差函数曲线图：以h为横坐标，C#(h)为纵坐标作图协方差函数曲线图：以h为横坐标，C#(h)为纵坐标作图63二、协方差函数的性质区域化变量Z(x)在二阶平稳假设下，其协方差函数存在且平稳，定义为1.C(0)=Var[Z(x)]≥0，即先验方差不能小于零2.C(h)=C(-h)，即C(h)对h=0的直线对称，是一个偶函数二、协方差函数的性质区域化变量Z(x)在二阶平稳假设下，其协64证：令x-h=y，则x=y+h，带入上式得图形特征及含义证：令x-h=y，则x=y+h，带入上式得图形特征及含义653.|C(h)|≤C(0)，即协方差函数绝对值小于等于先验方差证：3.|C(h)|≤C(0)，即协方差函数绝对值小于等于先664.|h|→∞时,C(h)→0，或写作C(∞)=0，即当空间距离很大时，协方差函数值很小意义(空间局限性)：当距离很大时，Z(x)和Z(x+h)之间的线性相关基本不存在5.C(h)必须是一个非负定函数，由C(xi-xj）构成的协方差函数矩阵必须是非负定矩阵正定条件(positivedefinitecondition)区域化变量Z(x)二阶平稳，其数学期望为m，协方差为C(h),变异函数为γ(h)，令Y是该类型区域化变量的任意有限线性组合，即：4.|h|→∞时,C(h)→0，或写作C(∞)=0，即当67较难理解则由C(xi-xj）(i,j=1,2…n)构成的协方差函数矩阵是非负定矩阵，即C(h)为非负定函数二阶平稳区域化变量的协方差函数是有条件的较难理解则由C(xi-xj）(i,j=1,2…n)构成的协方68三、实验(经验)变异函数(experimentalvariogram)的计算公式设区域化变量Z(x)满足(准)二阶平稳条件或(准)本征假设，h为两样本点空间分隔距离，Z(xi)与Z(xi+h)分别是Z(x)在空间位置xi和xi+h上的观测值(i=1,2,…N(h)),则计算实验变异函数的公式为：变异函数曲线图：以h为横坐标，γ#(h)为纵坐标作图三、实验(经验)变异函数(experimentalvari69变异函数计算实例（1）一维变异函数的计算x1x2x3x4x5x6x7x8x9x104

3

4

5797877以下为一研究对象在水平方向上的采样数据，满足二阶平稳或本征假设，采样值如图所示，点间分隔距离h=1米，计算γ#(h)变异函数计算实例（1）一维变异函数的计算x1x270变异函数及结构分析课件71变异函数及结构分析课件72两方面理解：变异性的理解与相关性的理解两方面理解：变异性的理解与相关性的理解73作业：x1x2x3x4x5x6x7x8

2

4

3

15364以下为一研究对象在水平方向上的采样数据，满足二阶平稳或本征假设，采样值如图所示，点间分隔距离h=1米，计算γ#(1),γ#(2),γ#(3)作业：x1x2x3x4x574（2）二维变异函数的计算下图为正方形网格状的采样数据，*号处为无数据点，点间距离h为100米，请分别计算南北、东西、西北和东南方向上的变异函数值。（2）二维变异函数的计算下图为正方形网格状的采样数据，*号处75变异函数及结构分析课件76变异函数及结构分析课件77变异函数及结构分析课件78西北和东南方向上的变异函数值的计算，注意分隔距离h的确定和样本数据对的查找西北和东南方向上的变异函数值的计算，注意分隔距离h的确定和样79变异函数及结构分析课件80作业：下图为正方形网格状的采样数据，网格交叉空白处为无数据点，点间距离h为a米，请分别计算南北方向γ#(a),

西北—东南方向上γ#(a)。作业：下图为正方形网格状的采样数据，网格交叉空白处为无数据点81四、变异函数的性质区域化变量Z(x)满足二阶平稳或本征假设条件，则变异函数存在且平稳，计算公式为1.γ(0)=0，即在h=0时，变异函数为零2.γ(h)=γ(-h)，即γ(h)对h=0的直线对称，是一个偶函数四、变异函数的性质区域化变量Z(x)满足二阶平稳或本征假设条823.γ(h)≥0,即研究现象的变异性只能大于或等于零4.|h|→∞时,γ(h)→C(0)，或写作γ(∞)=C(0)，即当空间距离很大时，变异函数值接近先验方差5.[-γ(h)]必须是一个条件非负定函数，即由[-γ(xi-xj)]构成的变异函数矩阵必须是条件非负定矩阵。3.γ(h)≥0,即研究现象的变异性只能大于或等于零4.83区域化变量Z(x)二阶平稳，其数学期望为m，协方差为C(h),变异函数为γ(h)，令Y是该类型区域化变量的任意有限线性组合，即：区域化变量Z(x)二阶平稳，其数学期望为m，协方差为C(h)84区域化变量Z(x)的变异函数γ(h)是有条件的，即需满足条件非负定条件区域化变量Z(x)的变异函数γ(h)是有条件的，即需满足条85五、协方差函数与变异函数的关系五、协方差函数与变异函数的关系86协方差函数和变异函数的曲线图问题：为什么只画出了h>0的关系图？当h足够大(即存在a>0，当h≥a)时，可以使C(h)=0,γ(h)=C(0),a称为变程(range)1、变程a表示区域化变量从存在空间相关状态(当|h|<a时）转向不存在空间相关状态(当|h|>a时)的转折点2、变程a的大小反映区域化变量影响范围的大小，或说反映该变量自相关范围的大小。也可说变程a是区域化变量空间变异尺度或空间自相关尺度变程a的意义：协方差函数和变异函数的曲线图问题：为什么只画出了h>0的关系87第二节变异函数的功能一、变异函数通过“变程”反映变量的影响范围——变异函数的跃迁现象变异函数γ(h)是一个单调递增函数，当h超过某一数值(变程a)后，γ(h)不再继续单调地增大，而往往稳定在一个极限值γ(∞)附近，这种现象称为“跃迁现象”(transitionphenomena)γ(∞)极限值称为基台值(sill),即C(0)【二阶平稳条件】，基台值的大小反映变量变化幅度的大小凡具有一个变程a和一个基台值的变异函数，称为“跃迁型”的变异函数“变程”反映变量的影响范围(图示）第二节变异函数的功能一、变异函数通过“变程”反映变量的影88二、不同方向上的变异函数图可反映区域化变量的各向异性——变异函数表示的各向异性如果在各个方向上区域化变量的变异性相同或相近，则称区域化变量是各向同性的，反之称为各向异性通过作出各个方向上的变异函数图，并放到一起来比较、分析、研究，就可以确定区域化变量的各向异性（包括有无各向异性，及各向异性的类型等）二、不同方向上的变异函数图可反映区域化变量的各向异性89变异函数及结构分析课件90三、块金常数C0的大小可反映区域化变量的随机性大小——变异函数的块金效应当h=0时，变异函数γ(h)≠0，而等于一个常数C0，这种现象称为“块金效应”(nuggeteffect),C0称为块金常数或块金方差(nuggetvariance)块金效应的图形表示三、块金常数C0的大小可反映区域化变量的随机性大小91“块金效应”

主要有两种来源：1、区域化变量在小于抽样尺度h时所具有的变异性2、采样分析误差当样点间的距离大于微域结构的范围，或样点的大小大于微域结构的范围就会出现块金效应（Webester，1985）“块金效应”主要有两种来源：1、区域化变量在小于抽样尺度h92四、变异函数在原点处的性状可反映区域化变量的空间连续性变异函数在原点处的性状主要有五种类型，每种类型反映了变量的不同程度的空间连续性

1、抛物线型（parabolictype）当|h|→0时,γ(h)→A|h|2(A为常数)，即变异函数曲线在原点处趋向一条抛物线，反映区域化变量是具有高度连续性的，如矿层厚度四、变异函数在原点处的性状可反映区域化变量的空间连续性变93

2、线性型（lineartype）当|h|→0时,γ(h)→A|h|(A为常数)，即变异函数曲线在原点处趋向一条直线，或说在原点处有斜向的切线存在，反映区域化变量是具有平均的连续性，如金属品位

3、间断型（discontinuoustype）当|h|→0时,γ(h)→C0，即变异函数曲线在原点处间断，说明块金效应存在，又称“块金效应型”，反映区域化变量的连续性很差，但当h增大时，γ(h)又变的较为连续了，如金品位2、线性型（lineartype）当|h|→0时,γ94

4、随机型（randomtype）这种变异函数可看成具有基台值C0和无穷小变程a的跃迁型变异函数，则无论h多小，h总大于a，故Z(x)与Z(x+h)总是互不相关又称纯块金效应型，反映了区域化变量完全不存在空间相关的情况，则本质上此区域化变量为普通随机变量此时，C0=C(0)4、随机型（randomtype）这种变异函数可看成具95

5、过渡型：介于抛物线型和随机型间当|h|→0时,γ(h)→C0，即有块金效应；当|h|=a时,γ(a)=C(0)，即有基台值(C0+C)和变程a,C称为“拱高”过渡型是实际研究工作中最常遇到的一种类型5、过渡型：介于抛物线型和随机型间当|h|→0时,γ(96第三节变异函数的理论模型

思考：是否有了采样数据及变异函数计算公式就可以获知任意距离h的区域化变量变异性？设Z(x)具有各向同性的变异函数

γ(h)，则常见的变异函数模型如下：变异函数的理论模型有基台值模型无基台值模型可以有或无基台值模型：孔穴效应模型球状模型、指数模型高斯模型线性有基台模型纯块金效应模型幂函数模型对数模型线性无基台模型第三节变异函数的理论模型思考：是否有了采样数据及变异函97一、有基台值模型

1、球状模型（sphericalmodel）若模型满足二阶平稳假设，且有有限先验方差，γ(h)值随h的变大而增大，当h达一定值(h>a)时，γ(h)达到一定值——基台值，则称此类模型为有基台值模型式中：C0为块金常数，(C0+C)为基台值，C为拱高，a为变程一、有基台值模型1、球状模型（sphericalmode98当C0=0，C=1，称为标准球状模型，其图形为：原点处切线的斜率为3/2a,与基台值线交点的横坐标为2a/3球状模型是地统计学应用最广的理论模型，许多区域化变量的理论模型都可以用球状模型来拟合当C0=0，C=1，称为标准球状模型，其图形为：原点处切线的99

2、指数模型（exponentialmodel）式中：C0，C意义同前，但a不是变程当C0=0，C=1，称为标准指数模型，其图形为：由于1-e-3=1-0.05=0.95≈1，则变程为3a2、指数模型（exponentialmodel）式中：100

3、高斯模型（gaussianmodel）式中：C0，C意义同前，但a不是变程由于1-e-3=1-0.05=0.95≈1，则变程为√3

a当C0=0，C=1，称为标准高斯函数模型，其图形为：3、高斯模型（gaussianmodel）式中：C0，101模型通过原点切线与基台值线交点的横坐标变程原点处的性状球状2a/3a直线指数a3a直线高斯无交点√3

a抛物线三种模型的比较模型通过原点切线与基台值线交点的横坐标变程原点处的性状球状2102

4、线性有基台值模型（linearwithsillmodel）式中：C0，C意义同前，A为常数，表示直线的斜率，变程为a4、线性有基台值模型（linearwithsillm103

5、纯块金效应模型（purenuggeteffectmodel）此时，C0=C(0)此种模型意味着区域化变量为随机分布，样点间的协方差函数对于所有距离h均等于0，即变量不存在空间相关性5、纯块金效应模型（purenuggeteffect104变异函数及结构分析课件105二、无基台值模型

1、幂函数模型（powermodel）若与模型相应的区域化变量不满