地质统计学5变差函数及结构分析cjg2011.pptVIP

变差图的拟合 拟合方法：如何加权？ 模型优选：球状？ 结构套合：各向异性？ 加权最小二乘拟合 加权系数的选择 加权系数的选择 模型优选 在所有的建立变差函数模型的过程中，我们都面临着一项任务，从众多基本的理论模型中选择适当的类型。 依靠实验变差在原点附近的性质 线性型——球状模型或者指数模型； 抛物线型——高斯模型。 不同方向上的结构套合 5. 比例效应 （1）比例效应的概念 设两个不同中心点位于x0，x0? 的准平稳邻域V(x0), V(x0 ?) ， γ(h, x0), γ(h, x0 ?) 分别为两邻域的变差函数。二者具有一定的比例关系： 这种现象称为结构的比例效应。即：存在一个不依赖于空间点x0及其邻域V(x0)的平稳结构γ0 (h) ，使得下式成立： 　又由于Z(x)在邻域内平稳，故在V(x0)邻域内有： 从而 　　　　也近于一个常数。 因此，当比例效应存在时，就存在一个不依赖于x0的平稳结构γ0（h）差一个依赖于x0的常数比例因子 　　　 。 正比例效应 反比例效应 （２）比例效应的类型 正比例效应：实验变差函数曲线 　　　　 为升函数，基台值随　　　　的上 　　　　　　升而上升； 反比例效应：实验变差函数曲线 　　　　 为降函数，基台值随　　　　的上 　　　　　　升而下降； 对数正态型分布的区域化变量 反对数正态型分布的区域化变量 （３）比例效应的消除 　　分析比例效应的目的就是为了如何消除这种比例效应，使之平稳。当存在以　　　　表示的正比例效应（即Z(x)具有对数正态型的直方图）时，消除比例效应的的方法有二： 　① 将原始数据取对数，然后再求其变差函数：即作对数变差函数； 　② 作相对变差函数：即以某邻域内的变差函数除以该邻域内观测数据的平均值的平方，从而得到与位置x0无关的相对变差函数γ０(h)，即： 　对于反比例效应，则可将上公式中的分母换成　　　　　　，从而得到消除反比例的相对变差函数。 （４）套合结构中的比例效应 　在准平稳性假设下，套合结构 　　　　　 的每一组成结构的比例效应可以是不同的，即对每一个模型 　　 可能有不同的 　　　　 。作为一个特例，当正比例效应仅对块金常数C0起作用时，套合结构为： 　在实际工作中，要想知道有无比例效应，是正还是反，以及这种比例效应属于哪一种结构，是否既作用于微观结构又作用于宏观结构，实用的方法就是在透明纸上绘制出一系列（对不同x0点及其邻域V(x0)）相对变差函数曲线，把它们套在一起进行直观比较，得出结论 三、变差函数的拟合 １.直观拟合球状模型 　（１） 求变程a （见图） 　　１）计算所用数据（用于计算实验变差函数γ*(h)的实验方差σ*2； 　　２）在实验变差函数图的纵坐标轴上过σ*点作一条平行于横坐标轴的直线； 　　３）以直线连接实验变差函数γ*(h)的头两至三个点，此直线与过σ*2点的直线相交，交点横坐标为 ，即假定其值为h0，则a = （２）求块金值C0：连接γ*(h) 的头两三个点的直线与纵坐标轴的交点的纵坐标为C0（见图）。若C00，则取C0＝0。 （３）求拱高C：C= σ*2- C0， 2. 用加权多项式回归法拟合球状模型 已知数据：　　 hi Ni γ*(hi) ? h1 N1 γ*(h1) ? h2 N2 γ*(h2) ?… … … ? hn Nn γ*(hn) ?变差函数： ? ? ? ?因h=0　和　ha　的情况都很简单，所以仅讨论0h≤a 的情况： ? 令：y= γ (h)，x1=h，x2=h3，b0=C0， 则上式变为：y= b0+b2x+b2x2 这样对球模型变差函数的拟合问题就变成了多元线性回归问题。 3. 用加权多项式回归法拟合二级套合球状模型 　　公式： ? ? 　　(1)??? 将实验变差函数分两个区分另进行拟合 　　　先取前区0h≤a１ 来拟合： ? ? ? 　　　按前述方法，令：y= γ (h)，x1=h，x2=h3，b0=C0， ? 　　(2)用后半区数据点a１h≤a２ 拟合第二段 ? ? ? 　　　设： y= γ (h)，x1=h，x2=h3，b0=C0＋C１， ?? 自动拟合 + 人机交互修正