地质统计学变异函数 53页

GS-Variogram变异函数与结构分析信息管理学院王玉兰E-mail:wyl@cdut.edu.cn,wang_wyl@163.comTel:84073385(o)\n变异函数与结构分析变异函数是对区域变量结构分析的工具也是研究对象非均质性描述的手段通过计算研究对象不同方向上的变异函数可以得到不同方向的结构；套合方法是对区域变量结构分析的方法均质化描述的方法；克立格估计方法就是利用结构特征估计的。\n变异函数变异函数：variogram//semi-variogram区域化变量在某方向上相距h的增量的方差，称为区域化变量在该方向上的变异函数，记为γ(x,h)；即γ(x,h)=½Var[Z(x)-Z(x+h)]=½E[Z(x)-Z(x+h)]2–½{E[Z(x)]-E[Z(x+h)]}2在二阶平稳假设或内蕴假设下有：γ(h)=½E[Z(x)-Z(x+h)]2在二阶平稳假设有：γ(h)=C(0)-C(h)\n变异函数变异函数的估计：在二阶平稳假设或内蕴假设下区域化变量在某方向上的变异函数为γ(h)；即Z(x)-Z(x+h)只依赖于在某方向上分隔它们的向量h，而与具体位置无关；则每对预测数据[z(x),z(x+h)]都可看成是[Z(x)-Z(x+h)]的一个取样(现实)，从而可用样本方差估计总体方差：γ*(h)=Σ[z(xi)-z(xi+h)]2/2N(h)\n本章的主要内容变异函数及变异曲线：讨论变异函数及曲线描述变异函数的理论模型：介绍变异函数的理论模型及特点实验变异函数曲线的计算与拟合：学会如何利用观测样本估计计算研究对象特征的变异函数值并用理论模型进行拟合；结构分析：对研究对象的异向性进行分析的方法，主要讨论不同异向性情况下的结构套合方法结构分析的实施步骤：介绍对研究对象从观测开始到特征描述的一般过程\n变异函数及变异曲线变异函数：由于其能反映区域化变量的结构特性，又称为结构函数；γ(h)=½E[Z(x)-Z(x+h)]2由于h和x是有方向的，一般描述：γ(h,α)=½E[Z(x)-Z(x+h)]2连续时：离散时：xx+hV\n变异函数及变异曲线变异曲线：γ(h,α)与h的关系曲线（h又称为滞后距）理想的变异曲线如图：C(0):先验方差a:变程(Range)C0：块金效应(nuggeteffect)C+C0：基台(sill)hCC0C(0)aγ(h)\n变异函数及变异曲线变异函数的性质：设Z(x)是二阶平稳的，则γ(h)存在且平稳，并有下列性质：(1)γ(0)=0(2)γ(h)>=0(3)γ(-h)=γ(h)(4)[-γ(h)]是条件非负定函数(5)γ(∞)=C(0)变异函数与协方差函数的关系曲线C(h)=C(0)-γ(h)hCC(0)aγ(h)\n变异函数及变异曲线变异函数的功能：1、通过变程a反映区域化变量的影响范围跃迁现象若区域化变量Z(x)的变异函数具有一个变程a和一个基台C，则Z(x)与落在以x为中心、以a为半径的邻域内的任何其它Z(x+h)有空间相关性；且相关程度随着两点距离的增大而减弱。hCC(0)aγ(h)\n变异函数及变异曲线变异函数的功能：2、变异函数在原点的性状可反映区域化变量的空间连续性抛物线型(Parabolictype连续型):|h|->0,r(h)->A|h|2线性型lineartype:|h|->0,r(h)->A|h|随机型randomtyper(0)=0;r(h)=C0;h>0间断型discontinuoustyper(0)=0;r(h)=C0;h->0可迁型(transitiontype)\n变异函数及变异曲线变异函数的功能：3、不同方向上的变异函数可反映区域化变量的和向异性hCC(0)aγ(h)\n变异函数及变异曲线变异函数的功能：4、变异函数如果是可迁的,则基台的大小反映该方向上的变化幅度hCC(0)aγ(h)\n变异函数及变异曲线变异函数的功能：5、块金常数C0的大小可反映区域化变量随机性的大小hCC(0)aγ(h)\n变异函数的理论模型设Z(x)是满足内蕴假设的区域化变量且具有各向同性的变异函数γ(h)。变异函数的理论模型分类：\n变异函数的理论模型变异函数理论模型函数形式：有基台的模型\n变异函数的理论模型变异函数理论模型函数形式：有基台的模型\n变异函数的理论模型变异函数理论模型函数形式：无基台的模型其它模型\n变异函数的理论模型变异函数理论模型函数形式：无基台的模型\n变异函数的计算与拟合变异函数的计算计算公式\n变异函数的计算与拟合变异函数的计算计算方法在指定方向上对指定h，搜索所有相距h的点对[z(xi),z(xi+h)]，并统计点对数N(h)。计算量依赖于数据的空间构型，按构型搜索方法可分为两类：规则数据构型：已知取样数据点在空间是按规律进行的；在指定方向上可得到基本滞后距不规则数据构型：横不成行竖不成列，找不到基本滞后距\n变异函数的计算与拟合变异函数的计算计算方法规则数据构型：小样数据利用基本滞后距h有规律地直接计算，大样数据抽样计算不规则数据构型：确定基本滞后距，给出角度容差和距离容差后计算,小样数据取大容差，大样数据取小容差\n变异函数的计算与拟合变异函数的计算（加快计算速度和减少计算量和存储量）点对文件法：(序号，角度，滞后距)存储量大TheVariogramGrid:极坐标网，原点为某观测点位置，方位角为计算方向，极距为相对于某观测点的滞后距\n变异函数的计算与拟合TheVariogramGrid:方向：8个方向滞后距：4个(100,200,300,400)最大滞后距400共32个网格单元\n变异函数的计算与拟合TheVariogramGrid:有三个观测点位置{(50,50),(100,200),(500,100)}.三个点对:A(50,50),(100,200)B(50,50),(500,100)C(100,200),(500,100)则三个点对在图中的位置A71.57158.11B6.34452.77C-14.04412.31\n变异函数的计算与拟合设Z(x)是一维区域化变量满足风蕴假设。有8个观测值如图y计算变异函数值γ(1)=3.0；γ(2)=1.67；γ(3)=2.80；γ(4)=2.87；γ(5)=1；γ(6)=4；变异函数图\n变异函数的计算与拟合设Z(x)是二维区域化变量满足风蕴假设。有41个观测值如图,网格边长为a。计算4个方向变异函数方向1：γ(a)=4.1；γ(2a)=8.84；γ(3a)=12.08；方向2:γ(a)=4.25；γ(2a)=8.22；γ(3a)=10.9；方向3：γ(1.414a)=5.03；γ(2.828a)=11.91；γ(4.242a)=17.25；方向4：γ(1.414a)=6.47；γ(2.828a)=11.25；γ(4.242a)=15.44；变异函数图\n变异函数的计算与拟合计算示例数据量:164最大滞后距：1618个方向10个滞后距共180个变异计算网格RangeMidrangeMinimumMaximumAverageVarianceCoef.ofVariationX1634.826.842.834.822.4Y302053520125Z79.9440.480.5180.454.3150669.6061.93347\n变异函数的计算与拟合变异函数的拟合用最优理论变异函数模型拟合计算出的变异函数；拟合的一般步骤：选择合适的理论模型；选择恰当的拟合方法；实际拟合。单一模型的拟合、多模型套合结构的拟合\n变异函数的计算与拟合变异函数的理论模型选择任意有基台的模型都可用球状模型拟合球状模型：单一模型的拟合球状模型的组合:多模型套合结构的拟合\n变异函数的计算与拟合人工拟合方法：观察-调整；最小二乘方法：确定模型-用最小二乘法求最佳参数；加权多项式回归方法：确定模型-用对应滞后距的点对数N(hi)为权系数，在最小二乘意义下拟合。\n变异函数的计算与拟合加权多项式回归拟合单个球状模型模型选择：单个球状模型模型变换:变成多项式形式加权处理：每个点对计算结果乘以权系数最小二乘计算：\n变异函数的计算与拟合加权多项式回归方法计算示例：b0=1.02;b1=2.28;b2=-0.01C0=1.02;a=8.312;C=12.634C0+C=13.654hr(h)1.54.12.36.33.28.03.99.54.710.75.411.76.212.5713.07.713.38.513.59.213.7\n变异函数的计算与拟合加权多项式回归拟合二级套合球状模型模型选择：二级球状模型模型变换:变成多项式形式加权处理：每个点对计算结果乘以权系数最小二乘计算：分段用最小二乘拟合计算\n变异函数的计算与拟合加权多项式回归拟合二级套合球状模型模型选择：二级球状模型模型变换:变成多项式形式加权处理：每个点对计算结果乘以权系数最小二乘计算：分段用最小二乘拟合计算\n变异函数的计算与拟合加权多项式回归方法计算示例：C0=0.3353;a1=2.198;C1=0.8665a2=8.1932;C2=1.43hr(h)Ni0.50.731911.16181.51.441821.65172.51.791732163.52.071542.18154.52.251452.3145.52.441362.53136.52.511272.59127.52.6411\n变异函数的计算与拟合加权多项式回归方法计算示例：C0=0.273;a=7.383;C=5.4733hr(h)0.91.31.92.23.1344.54.85.86.24.77.96.19.45\n变异函数的计算与拟合人工拟合方法：观察-调整；最小二乘方法：确定模型-用最小二乘法求最佳参数；加权多项式回归方法：确定模型-用对应滞后距的点对数N(hi)为权系数，在最小二乘意义下拟合。\n区域化变量结构分析实验变异函数-用理论函数似合-区域化变量分析结构分析-构造变异函数模型对全部有效结构信息作定量化的概括，以表征区域化变量的主要特征。主要使用方法套合结构法（neststructure);套合结构是把分别出现在不同滞后距h上和不同方向a上同时起作用的变异性组合起来。形成一个统一的描述区域化变量特征的变异函数套合结构表示为：γ(h)=γ0(h)+γ1(h)+γ2(h)+…+γ(h)+…\n区域化变量结构分析一个方向上的套合结构套合结构表示为：γ(h)=γ0(h)+γ1(h)+γ2(h)+…+γi(h)+…套合结构中每个变异函数模型代表一个特定尺度上的变异性，可以是用不同模型表示变异函数；例如：某区域化变量在某一方向上的变异性可用三个层次的变异性表示：γ(h)=γ0(h)+γ1(h)+γ2(h)\n区域化变量结构分析不同方向的套合结构概念：各向同性：区域化变量在各个不同方向上性质相同；各向异性：不同方向上性质不同；各向异性表现为不同方向的变异函数的差异；各向异性的分类及特征：几何异向性和带状异向性\n区域化变量结构分析不同方向的套合结构几何异向性：（geometricanisotropy)区域化变量在不同方向上表现出变异程度相同（基台相同）而连续性不同（变程不同）的异向性为什么称为几何异向性？这种异向性可以通过简单几何变换变为各向同性；各向异性比K:变程比a1/a2：表示不同方向的变异性之差异大小各向异性表现为不同方向的变异函数的差异；\n区域化变量结构分析不同方向的套合结构二维几何异向性变换：水平方向球状模型1：参数C，a2;垂直方向球状模型2：参数C，a1;各向异性比K=a2/a1；\n区域化变量结构分析不同方向的套合结构带状异向性：（zonalanisotropy)区域化变量在不同方向上表现出的变异性之差不能简单用变程比得到，即无法用几何变换将它们变为各向同性为什么称为带状异向性？不同的变异性只局限在该方向带中。例：多层状矿区，由于矿层及夹层组分变化显著，垂直方向的变异比平面上的变异大。垂直方向含有重力作用因素，矿化作用包括成矿元素的自然扩散和重力双重作用。平面γ1(h)γ(hw)=γ1(h)+γ2(hw)\n区域化变量结构分析几何异向性的套合套合的目的：如何将各向异结构变换为各向同性结构；或用统一的一个变异函数表示水平方向球状模型1：参数C，a2;垂直方向球状模型2：参数C，a1;各向异性比K=a2/a1；\n区域化变量结构分析几何异向性的套合几何异向性的套合设两个方向都是球状模型变程分别是a4和a1,基台为c(b)K1=a4/a1h*=(hu2+(K1hv2)1/2\n区域化变量结构分析带状异向性的套合套合的方法：1）把各个方向的结构各自当成独立的成份进入套合结构。2）把变化最简单的某个方向的结构看成是各向同性，在其它方向上做叠加；\n区域化变量结构分析一般套合结构套合结构：基本结构的叠加组合\n区域化变量结构分析结构分析具体实施方法1）选择符合研究目的的区域化变量并获取相应数据。如品位、厚度、饱和度、孔隙度、渗透率等2）对被研究的数据进行审议。包括数据支撑等3）进行数据统计分析。了解数据的基本情况4）计算实验变异函数。要求根据数据的具体情况选择合适的计算方法；5）进行变异函数拟合和结构套合，获取区域化变量的结构模型。6）结构模型的检验。交叉检验法。离差检验法。7）对变异函数进行地质解释。\n结构分析示例水厂北山铁矿可溶铁(SFe)变异函数的套合分析背景:前寒武纪沉积变质型铁矿,22米中段经下SFe的变异函数在水平及垂直方向呈现出几何异向性,且在这两个方向上均由两个球状模型结构组成,垂直方向:15.6,53.3sph(60),41.08sph(228)水平方向:15.6,53.3sph(131.5),41.08sph(500)总基台高:15.6+53.3+41.08=109.98异性比K=a水平/a垂直=2.19线性变换矩阵A的对角线元素为[1,1,2.19]统一后的模型结构:r(h)=15.6+53.3sph(131.5)+41.08sph(500)\n结构分析示例水厂北山铁矿可溶铁(SFe)变异函数的套合分析矿区196-40m中段的SFe的变异函数表现为带状异向性垂直方向:15.6,53.3sph(60),41.08sph(228)水平方向:15.6,44,36sph(130),9.95sph(500)套合1:r(h)=15.6+53.3sph(60)+41.08sph(228)+44,36sph(130)+9.95sph(500)变换矩阵:[0,0,1];[1,1,0]垂直-水平:0,25sph(40),17sph(228)套合2:r(h)=15.6+44,36sph(130)+9.95sph(500)+25sph(40)+17sph(228)变换矩阵:[1,1,1];[0,0,1]\n结构分析示例湖北大冶铁矿尖山矿体Fe\Cu\S结构分析背景:该矿为一大型铁铜矽卡岩型矿床尖山矿体Fe变异系数20.97%,虽存在贫富两种矿,但品位总体离散不大,所以将矿体看成统一体计算变异函数并进行分析,计算r(h)时样品用6m组合样。Fe变异函数曲线呈各向同性,且可用两个球状模型拟合:r(h)=15+30sph(60)+170sph(66)Cu品位分布服从三对数正态分布：ln(Cu+0.3)服从正态分布,其变异函数曲线呈几何异向性：各向同性结构：0.015,0.08sph(36)水平结构：0.114sph(120);垂直结构：0.114sph(225)套合结构:r(h)=0.015+0.08sph(36)+0.114sph(225,kh),k=1.9S呈对数正态分布,ln(S)呈几何异向性各向同性结构：0.14sph(24)水平结构：0.256sph(100);垂直结构：0.256sph(140)套合结构:r(h)=0.14sph(24)+0.256sph(140,kh),k=1.4小结构:富矿规模和富集程度含矿率计算结果:各向同性r(h)=0.034+0.035sph(80)\n结构分析具体实施方法\n本章的主要内容变异函数及变异曲线：讨论变异函数及曲线描述变异函数的理论模型：介绍变异函数的理论模型及特点实验变异函数曲线的计算与拟合：学会如何利用观测样本估计计算研究对象特征的变异函数值并用理论模型进行拟合；结构分析：对研究对象的异向性进行分析的方法，主要讨论不同异向性情况下的结构套合方法结构分析的实施步骤：介绍对研究对象从观测开始到特征描述的一般过程