kriging基础知识PPT.ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,克里金插值,第二讲,克里金方法（,Kriging,）,是以南非矿业工程师,D.G.Krige(,克里格,),名字命名的一项实用空间估计技术，是,地质统计学,的重要组成部分，也是地质统计学的核心。,1,地质统计学,主要是为解决矿床储量计算和误差估计问题而发展起来的,由法国巴黎国立高等矿业学院,G,马特隆教授于,1962,年所创立。,2,H.S.Sichel(1947),D.G.Krige(1951),Kriging,法,（克里金法，克立格法）,:“,根据样品空间位置不同、样品间相关程度的不同，对每个样品品位赋予不同的权，进行滑动加权平均，以估计中心块段平均品位”,G.,Materon(1962),提出了,“地质统计学”概念,(,法文,Geostatistique),发表了专著,应用地质统计学论,。阐明了一整套区域化变量的理论，为地质统计学奠定了理论基础。,区域化变量理论,克里金估计,随机模拟,应用统计学方法研究金矿品位,1977,年我国开始引入,3,克里金插值方法,井眼,地震,（普通克里金）,（应用,随机函数,理论）,不仅考虑待估点位置与,已知数据位置的相互关,系，而且还考虑变量的,空间相关性。,4,为一个实值变量，可根据概率分布取不同的值。每次取值（观测）结果,z,为一个确定的数值，称为随机变量,Z,的,一个实现。,P,一、随机变量与随机函数,第一节 基本原理,1.,随机变量,5,连续变量：,累积分布函数（,cdf,）,cumulative distribution function,条件累积分布函数（,ccdf,）,后验,conditional cumulative distribution function,离散变量（类型变量）：,Z(u),P,P,不同的取值方式：,估计（,estimation,）,模拟（,simulation,）,6,连续型地质变量,构造深度,砂体厚度,有效厚度,孔隙度,渗透率,含油饱和度,离散型地质变量,（范畴变量）,砂体,相,流动单元,隔夹层,断层,类型变量,7,设,离散型随机变量,的所有可能取值为,x,1,，,x,2,，,，其相应的概率为,P(=x,k,)=p,k,k=1,2,.,随机变量的特征值：,(1),数学期望,是随机变量,的整体代表性特征数。,则当级数 绝对收敛时，称此级数的和为,的数学期望，记为,E(),，或,E,。,E()=,8,设,连续型随机变量,的可能取值区间为,(-,，,+),，,p(x),为其概率密度函数，若无穷积分,绝对收敛，则称它为,的数学期望，记为,E(),。,E()=,数学期望是随机变量的最基本的数字特征，,相当于随机变量以其取值概率为权的加权平均数。,从矩的角度说，数学期望是,的一阶原点矩。,对于一组样本：,9,为随机变量,的离散性特征数。若数学期望,E-E(),2,存在，则称它为,的方差，记为,D(),，或,Var(),，或,2,。,=,从矩的角度说，方差是,的二阶中心矩。,(2),方差,其简算公式为,D()=E(,2,)E(),2,D()=E-E(),2,方差的平方根为标准差，,记为,10,研究范围内的一组随机变量。,简记为,随机场：,当随机函数依赖于多个自变量时，称为随机场。,如具有三个自变量,(,空间点的三个直角坐标,),的随机场,2.,随机函数,条件累积分布函数（,ccdf,）,P,11,二个随机变量,，,的协方差为二维随机变量,(,，,),的二阶混合中心矩,11,，记为,Cov(,，,),，或,。,协方差,(Variance),：,Cov(,)=,=E-E()-E(),其简算公式为,Cov(,)=E()-E()E(),随机函数的特征值,12,二、统计推断与平稳要求,P,任何统计推断（,cdf,，数学期望等）均要求重复取样。,但在储层预测中，一个位置只能有一个样品。,同一位置重复取样，得到,cdf,，不现实,13,考虑邻近点，推断待估点,空间一点处的观测值可解释为一个随机变量在该点,处的一个随机实现。,空间各点处随机变量的集合构成一个随机函数。,区域化变量：,能用其空间分布来表征一个自然现象的变量。,（将空间位置作为随机函数的自变量）,（可以应用随机函数理论解决插值和模拟问题）,14,考虑邻近点，推断待估点,-,空间统计推断要求平稳假设,严格平稳,对于单变量而言：,可从研究区内所有数据的累积直方图推断而得,（将邻近点当成重复取样点）,太强的假设，不符合实际,P,15,当区域化变量,Z(u),满足下列二个条件时，则称其为二阶平稳或弱平稳：,EZ(u)=EZ(u+h),=m(,常数,),x,h,随机函数在空间上的变化没有明显趋势，,围绕,m,值上下波动。,在整个研究区内有,Z(u),的数学期望存在，,且等于常数，即：,二阶平稳,16,在整个研究区内，,Z(u),的协方差函数存在且平稳,(,即只依赖于滞后,h,，而与,u,无关,),，即,CovZ(u),Z(u+h),=EZ(u)Z(u+h)-EZ(u)EZ(u+h),=EZ(u)Z(u+h)-,=C(h),特殊地，当,h=0,时，上式变为,VarZ(u)=C(0),即方差存在且为常数。,协方差不依赖于空间绝对位置，而依赖于相对位置，,即具有空间的平稳不变性。,u,u+h,17,在整个研究区内有,EZ(u)-Z(u+h)=0,本征假设,当区域化变量,Z(u),的增量,Z(u)-Z(u+h),满足下列二条件时，称其为满足本征假设或内蕴假设。,可出现,EZ(u),不存在，,但,EZ(u)-Z(u+h),存在并为零的情况,intrinsic hypothese,EZ(u),可以变化,但,EZ(u)-Z(u+h)=0,（比二阶平稳更弱的平稳假设）,18,增量,Z(u)-Z(u+h),的方差函数,(,变差函数,Variogram),存在且平稳,(,即不依赖于,u),，即：,VarZ(u)-Z(u+h),=EZ(u)-Z(u+h),2,-,EZ(u)-Z(u+h),2,=EZ(u)-Z(u+h),2,=2(u,h),=2(h),相当于要求：,Z(u),的变差函数存在且平稳。,19,例：物理学上的著名的布朗运动是一种呈现出无限离散性的物理现象，其随机函数的理论模型就是维纳,-,勒维,(Wiener-Levy),过程,(,或随机游走过程,),。,布朗运动,:,可出现协方差函数不存在，但变差函数存在的情况。,既不能确定验前方差，也不能确定协方差函数。,但是其增量却具有有限的方差：,VarZ(x)-Z(x+h)=2 =A|h|(,其中，,A,是个常数,),，,变差函数,=|h|,，且随着,|h|,线性地增大。,20,若区域化变量,Z(x),在整个区域内不满足二阶平稳,(,或本征假设,),，但在有限大小的邻域内是二阶平稳,(,或本征,),的，则称,Z(x),是准二阶平稳的,(,或准本征的,),。,准二阶平稳假设及准本征假设,21,设 为区域上的一系列观测点，为相应的观测值。区域化变量在 处的值 可采用一个线性组合来估计：,三、克里金估计,(,基本思路,Z*(x,0,),无偏,最优,无偏性,和,估计方差最小,被作为 选取的标准,-,以普通克里金为例,22,从本征假设出发,可知 为常数,有,可得到关系式,:,（,1,）无偏条件,Z*(x,0,),（在搜寻邻域内为常数，不同邻域可以有差别）,23,（,2,）估计方差最小,应用拉格朗日乘数法求条件极值,Z*(x,0,),24,进一步推导，可得到,n+1,阶的线性方程组,即克里金方程组,当随机函数不满足二阶平稳,而满足内蕴（本征）假设时,可用变差函数来表示克里金方程组如下,:,Z*(x,0,),25,最小的估计方差,即克里金方差可用以下公式求解,:,Z*(x,0,),26,变差函数,(,或叫,变程方差函数,，或,变异函数,),是地质统计学所特有的基本工具。它既能描述区域化变量的空间结构性变化，又能描述其随机性变化。,跃迁现象,1.,变差函数的概念与参数,四、变差函数及其结构分析,27,假设空间点,x,只在一维的,x,轴上变化，则将区域化变量,Z(x),在,x,，,x+h,两点处的,值之差,的方差之半定义为,Z(x),在,x,轴方向上的变差函数，记为,一维情况下的定义：,VarZ(x)-Z(x+h),EZ(x)-Z(x+h),2,-EZ(x)-Z(x+h),2,=,=,半变差函数,(,或半变异函数,),28,在,二阶平稳假设，或作本征假设,，此时：,地质统计学中最常用的基本公式之一。,EZ(x)-Z(x+h)=0,h,VarZ(x)-Z(x+h),EZ(x)-Z(x+h),2,-EZ(x)-Z(x+h),2,=,=,EZ(x)-Z(x+h),2,=,则：,29,（二阶平稳假设条件下边查函数与写防查的关系）,30,变程,(Range),：,指区域化变量在空间上具有相关性的范围。在变程范围之内，数据具有相关性；而在变程之外，数据之间互不相关，即在变程以外的观测值不对估计结果产生影响。,31,具不同变程的克里金插值图象,32,块金值,(Nugget),：变差函数如果在原点间断，在地质统计学中称为“块金效应”，表现为在很短的距离内有较大的空间变异性，无论,h,多小，两个随机变量都不相关。它可以由测量误差引起，也可以来自矿化现象的微观变异性。在数学上，块金值,c0,相当于变量纯随机性的部分。,33,如果品位完全是典型的随机变量，则不论观测尺度大小，所得到的实验变差函数曲线总是接近于纯块金效应模型。,当采样网格过大时，将掩盖小尺度的结构，而将采样尺度内的变化均视为块金常数。这种现象即为块金效应的尺度效应。,块金效应的尺度效应,1,2,1,1,1,3,3,3,3,34,基台值,(Sill),：,代表变量在空间上的总变异性大小。即为变差函数在,h,大于变程时的值，为,块金值,c,0,和,拱高,cc,之和。,拱高,为在取得有效数据的尺度上，可观测得到的变异性幅度大小。当块金值等于,0,时，基台值即为拱高。,=C(0)C(h),35,几何各向异性：,变差函数在空间各个方向上的,变程不同,，但,基台值不变,（即变化程度相等）。这种情况能用一个简单的几何坐标变换将各向异性结构变换为各向同性结构。,带状各向异性：,不同方向的变差函数具有,不同的基台值,，其中,变程可以不同，也可以相同,。这种情况不能通过坐标的线性变换转化为各向同性，因而结构套合是比较复杂的。,地质变量相关性的各向异性,1,2,1,1,1,3,3,3,3,(2),36,2.,变差函数的理论模型,设,Z(x),为满足本征假设的区域化变量，则常见的理论变差函数有以下几类：,球状模型,指数模型,高斯模型,幂函数模型,空洞效应模型,37,接近原点处，变差函,数呈线性形状，在变,程处达到基台值。,原点处变差函数的切,线在变程的,2/3,处与,基台值相交。,球状模型：,c,为基台值，,a,为变程，,h,为滞后距。,38,指数模型：,变差函数渐近地逼近,基台值。,在实际变程处，变差,函数为,0.95c,。,模型在原点处为直线。,39,高斯模型：,变差函数渐近地逼近,基台值。,在实际变程处，变差函,数为,0.95c,。,模型在原点处为抛物线。,40,幂函数模型：,幂函数模型为一种无基台值的变差函数模型。这是一种特殊的模型。,当,=1,时，变差函数为一直线，即为线性模型，这一模型即为著名的,布朗运动（随机行走过程）,的变差函数模型；,当,1,时，变差函数为抛物线形状，为,分数布朗运动,(fBm),的变差函数模型。,布朗运动,分数布朗运动,分数布朗运动,h,41,空洞效应模型,(Hole Effect),：,变差函数并非单调增加，而显示出一定周期性的波动。,模型可以有基台值，也,可以无基台值；可以有,块金值，也可以无块金,值。,空洞效应在地质上多沿,垂向上出现，如富矿层,与贫矿层互层、砂岩与,泥岩频繁薄互层等等。,（,b,为富矿化带重复距离）,h,42,通过区域化变量的空间观测值来,构建相应的变差函数模型,以表征该变量的主要结构特征。,(,求变差,),(1),数据准备,区域化变量的选取,、,数据质量检查及校正,、,数据的变换,（如对渗透率进行对数变换）、,数据的统计,（如分相对储层参数计算平均值、,方差，作直方图、相关散点图等）、,丛聚数据的解串,等。,3.,区域化变量的,结构分析,43,(2),实验变差函数的计算,实验变差函数是指应用观测值计算的变差函数。对于不同的滞后距,h,，可算出相应的实验变差函数,。,=,一维实验变差函数的计算公式,(i=1,N(h),Z(x,i,)-Z(x,i,+h),2,的算术平均值一半即为一个,h,的变差函数值,44,对不同的滞后,h,，进行计算，得出各个,h,的变差函数值,=,h,3h,5h,h,45,设,Z(x),为一维区域化变量，满足,本征假设,，又已知,Z(1)=2,，,Z(2)=4,，,Z(3)=3,，,Z(4)=1,，,Z(5)=5,，,Z(6)=3,，,Z(7)=6,，,Z(8)=4,，,，,例：,试求：,=,46,=,=,=,2,2,+1,2,+2,2,+4,2,+2,2,+3,2,+2,2,=,=3.00,1,2,+3,2,+2,2,+2,2,+1,2,+1,2,=,=1.67,1,2,+1,2,+0,2,+5,2,+1,2,=,=2.80,47,2D,情况,（,1,）分不同方向,，进行,1D,变差函数计算,3D,情况,：,增加垂向方向,（,2,）确定主变程方向,次变程方向,角度容限,步长容限,h,3h,5h,h,四方向试算,（考虑主变程方向的 走向、倾向和倾角）,48,(3),理论变差函数的最优拟合与结构套合,选择合适的理论变差函数模型，同时还需进行结构套合，从而得到一条反映不同层次（或不同空间规模）结构的、统一的、最优的变差函数曲线。,球状模型,指数模型,高斯模型,幂函数模型,空洞效应模型,49,复杂的区域化变量往往包含各种尺度上的多层次、多方向的变化性，反映在变差函数上即为多层次结构。将不同结构组合为统一结构的过程称为“结构套合”,结构套合,各层次套合,例如，对于,200,米宽的河道，在,h,50m,的观测尺度上可以将,其与河道间的变化性,区分出来，但却无法区分,层理和矿物成分的变化性,(,即无法找出更细微的结构来,),，它们在,50m,尺度得到的结构上只能作为“块金效应”出现。若观测尺度为,500,米，河道的变化也只能作为“块金效应”。,1,2,1,1,1,3,3,3,3,大尺度的变化性总是包含着小尺度的变化性，但却不能从大尺度的变化性中区分出小尺度的变化性。,50,=,=,=,代表微观变化性的变程极小的球状模型，可近似地看作纯块金效应型,球状模型，没有块金常数，基台值为,C,1,，变程为,a,1,，反映了小规模范围的变化,球状模型，没有块金常数，基台值为,C,2,，变程较大，为,a,2,，反映了大规模范围的变化,可以用反映各种不同尺度变化性的多个变差函数之和来表示一个套合结构。,（各层次理论模型可以不一样）,可以是不同模型的变差函数,51,其中,则套合结构的表达式为,=,=,=,=,52,对于,几何各向异性,，先根据异向比压缩,距离轴，使之成为各向同性的模型；,对于,带状各向异性,，运用模型叠加的方法加以处理。先用压缩距离轴的办法，使其变程变为相同，然后再把具有相同变程的两个球状模型叠加起来，构成一个新的球状模型,各方向套合,（,将各向异性套合为各向同性,，以便于,在克里金估计时，不同方向均可用统一,的结构模型计算实际的变差函数值）,53,(4),变差函数参数的最优性检验：,变差函数是否符合实际，应该进行检验。一种实用的检验方法为,“,交叉验证法,”,（,Cross-validation,），检验标准是在各实测点，根据周围点计算的,克里金估计值与该实测值的误差平方,平均最小。,估计误差的平方,与,克里金估计方差,之比越接近,1,，则说明变差函数与实际的符合程度越高。实际上，这种方法在检验变差函数的同时，也在检验所使用的克里金估计方法的适用性。,Z*(x,0,),54,五、克里金插值中权系数的确定,（以普通克里金为例）,i,求取变差函数（或协方差）；,解克里金方程组,在结构分析的基础上,55,设有一个油藏，在平面上,S,1,，,S,2,，,S,3,，,S,4,处有四个井点，其孔隙度值分别为,Z,1,，,Z,2,，,Z,3,，,Z,4,。据此估计,S,0,点处的孔隙度值,Z,0,设孔隙度,Z(x),是二阶平稳的。其在平面上的二维变差函数是一个各向同性的球状模型，其参数为：块金值,C,0,2,，变程,a,200,，拱高,C,20,，即：,实例,56,Z,0,的估计量为,普通克里金方程组的矩阵形式为,K =M,2,(,求解）,57,求解,:,C,ij,C11,C12,C01,试求,?,58,C,11,=C,22,=C,33,=C,44,=C(0)=,2,=C,0,+C=22,？,由于,C(h)=C(0)-(h)=22-(h),当,ij,时，,C,ij,=C(|S,i,-S,j,|)=22-(|S,i,-S,j,|).,于是，,C,12,=C,21,=C,04,=22-,59,将以上数值代入普通克里金方程组解的矩阵形式中，得,经计算得,:=0.5182,=0.0220,=0.0886,=0.3712,。,Z,0,0.5182Z,1,+0.0220Z,2,+0.0886Z,3,+0.3712Z,4,60,六、克里金插值中搜索策略,搜索邻域,注意,1,：,搜索邻域中的数据点才参加估计,节省,CPU,和内存,局域平稳,搜索椭圆或椭球的选择方法与选择变差函数椭圆或椭球相同。,61,注意,2,：,参与计算的数据点不能太多，否则计算太慢,一般软件中都,内置或可选,最大的 数据点数目（与待估点最近的数据点），如,10,。,注意,3,：,防止数据丛聚带来的数据代表性不强,井眼,井眼垂向数据太密，若待估点与该井近，则可能忽视邻井数据,八分搜寻,，保证各象限均有代表数据,62,若搜寻范围无数据，则应用边际概率。,63,x,0,第二节,克里金插值方法,简单克里金,(SK),普通克里金,(OK),泛克里金,(UK),协同克里金,(CK),贝叶斯克里金（,BK,）,指示克里金,(IK),64,一、简单克里金,(SK),(Simple Kriging),所有克里金估计都应用线性回归算法，形式为：,m,为期望,求取权系数的克里金方程组的非平稳形式,求,(n+1),个,m(u),求,(n+1),(n+1),个,C(u,u,),65,二阶平稳假设,EZ(u)=EZ(u+h)=m(,常数,),C(u,u+h)=C(h),简单克里金估计的平稳形式：,EZ(u)=EZ(u+h)=m(,常数,),66,应用条件：,随机函数二阶平稳,随机函数的期望值,m,为常数并,已知,不能用于具有局部趋势的情况,简单克里金方程组的平稳形式：,C(u,u+h)=C(h),（,C,与位置有关）,（,C,与位置无关）,67,二、普通克里金,(OK),(Ordinary Kriging),68,应用要求：,随机函数二阶平稳或符合内蕴假设,随机函数的期望值,m,在搜寻邻域内稳定但未知,协方差平稳,与简单克里金相比，普通克里金相当于在每一个位置,u,，重新估计,m,。,由于普通克里金估计常使用滑动数据邻域，相当于均值,m,随位置可变，即,Z*(u),，此时，实际上是一种非平稳算法，对应于变化的均值和平稳的协方差。,69,三、泛克里金,(UK),(Universal Kriging),非平稳随机函数的漂移函数,(drift),简称为漂移或趋势,随机函数,=,趋势,+,残差,区域化变量,Z(X),是非平稳的，即,EZ(x)=m(x),Kriging with a trend model (KT),具有趋势的克里金,70,用光滑的确定性函数来模拟，或用拟合方法,趋势函数,一维的线性趋势,二维的二次趋势,:,71,用均值为,0,、,协方差函数为 的平稳随机函数来模拟。,泛克里金估计值,:,残差,72,为权值,是与,(,K,+1),个权值的限制条件相对应的,(,K,+1),个拉格朗日参数,泛克里金方程组,为残差协方差函数,73,具有外部漂移的克里金,(Kriging with an external Drift),估计值,当,K,=1,时，线性趋势函数为,趋势函数,可理解为二级变量,74,（,1,）外部变量必须在空间光滑,地变化，否则可能导致,KT,线性系统不稳定；,（,2,）在主变量的所有数据点,u,处和待估计的,位置,u,处，外部变量都必须是已知的。,克里金,方程组：,可理解为地震,数据（如深度）,（,K=0,时，？）,75,利用几个变量之间的空间相关性，对其中的一个或几个变量进行空间估计，尤其适用于被估计变量的观察数据较少的情况。,协同克里金估计值（初始变量和二级变量）,四、协同克里金（,CK,）,(Cokriging),-,随机变量在位置,0,处的估计值；,-,初始变量的,n,个样本数据；,-,二级变量的,m,个样本数据；,-,需要确定的协同克里金加权系数。,及,76,协同克里金方程组,传统普通协克里金,77,标准化普通协克里金,m,X,=Ex(u),m,Y,=Ey(u),78,为协同克里金的简化形式，即如果二级变量密集取样时，只保留与估计点同位的二级变量。,对应的协同克里金方程组只要求知道,Z-,协方差函数,以及,Z-Y,互协方差函数,C,Z,（,h,）,（同位两种数据的,相关系数）,（方差函数）,同位协同克里金,Collocated Cokriging,79,五、贝叶斯克里金（,BK,）,H.Omre,在（,1987,）把线性贝叶斯理论用于克里金估计技术，提出了贝叶斯克里金估计技术。他构想了一个模型，把用于空间估计的数据分为两类：,观察数据：,是指那些精度比较高，但数量比较少的数据,猜测数据：,是指那些精度比较低，但分布广泛的数据,在观测数据比较多的地方，估计结果主要受观测数据的影响；在观测数据比较少的地方，则主要受猜测数据的影响。,显然，井数据和地震数据的关系符合贝叶斯估计中观测数据和猜测数据的关系。,80,设,Z(x),，,xA,，是观察数据的区域化变量。,设,M(x),，,xA,，是猜测数据的区域化变量。,Z*(x,0,),=,EZ(x,0,),=a,0,+,M,(x,0,),EM(x),M,(x),，,xA,M(x),是对,Z(x),的一种猜测，误差为,a,0,x,0,a,0,？,81,设已得到,Z(x),，,xA,的一组,(N,个,),观察值,Z(x,i,);i=1,2,N,。,定义一个新的随机函数：,Z,T,(x),Z(x)-,M,(x),，,xA,Z,T,(x,i,),=Z(x,i,)-,M,(x),，,i=1,2,N,Z(x,0,),的贝叶斯克里金估计量为,x,0,对这个,N,个观察值有,（相当于误差,a,0,）,（误差的随机函数）,82,基于无偏性和估计方差最小两个条件：,利用拉格朗日乘数法可得到贝叶斯克里金方程组：,Z,(x,x,)=,Z,M,(x,-x,)+,M,(x,x,),83,将数据按照不同的门槛值编码为,1,或,0,的过程。,对于模拟目标区内的每一类相，当它出现于某一位置时，指示变量为,1,，否则为,0,。,A,(100),B,(010),A,(100),C,(001),类型变量的指示变换：,变量,u,属于范畴,A,其它,六、指示克里金,(IK),Indicator Kriging,指示变换,1982,年由,AGJournel(,儒尔奈耳,),教授提出,84,(00111),(00001),(01111),(00011),首先将连续变量截断为类型变量，然后进行指示变换。,如：,z=10,，,15,，,20,，,25,30,连续变量的指示变换,85,设沿空间某一方向，在间距为,h,的,5,对样品点处观测了,Z(x,),及,Z(x,+h),的值,(=1,2,，,，,5),。,设指示,X,I(x,;z)=,设指示,Y=I(x,+h;z)=,86,假定只选定了,5,个门限值：,0,，,2,，,3,，,6,，,9,87,指示,(,函数,),的数学期望,当,x,固定时，若,z,再给定，则,I(x,；,z),就是个随机变量，就有数学期望：,EI(x,；,z),1PI(x,；,z)=1+0PI(x,；,z)=0,=PI(x,；,z)=1=P =F(x,；,z),，,z+,在,x,点处区域化变量,Z(x),的先验分布函数,F(x,；,z),，,z+,就是,x,点处指示函数的数学期望,E,I(x,；,z),i,*,(x,；,z)=,88,(1)z=3,时，设指示随机变量,X,I (x;3),E(X)=EI(x;3)=,V,ar,(X)=V,ar,I(x;3)=m,x,(1-m,x,)=0.60.4=0.24=,试求指示随机变量的数学期望和方差：,89,(2)z=6,时，设指示随机变量,Y,I(x+h;6),E(Y)=EI(x+h;6)=,V,ar,(Y)=V,ar,I(x+h;6)=m,Y,(1-m,Y,)=0.40.6=0.24=,90,设,Z(X),是个一维区域化变量，在等间距的,10,个点 处有,10,个观测值：,Z(1)=3,，,Z(2)=5,，,Z(3)=6,，,Z(4)=2,，,Z(5)=7,，,Z(6)=1,，,Z(7)=4,，,Z(8)=8,，,Z(9)=9,，,Z(10)=7,，设门限值,Z,分别等于,2,3,4,5,6,7,8,时，求指示变差函数的估计值,(h;z),，,h=1,，,2,，,3,，,4,，,5,。,计算,91,I,(h;5)=,(h;5)=,首先，计算,i(x,;5),：,i(x,1,;5),1,i(x,2,;5),1,i(x,3,;5),0,i(x,4,;5)=1,i(x,5,;5),0,i(x,6,;5)=1,i(x,7,;5),1,i(x,8,;5),0,i(x,9,;5),0,i(x,10,;5),0,(1;5)=,92,列表计算各指示值,i(x,;z),根据,i(x,;z),值算出不同,z,值的,(h;z),值,Z(X),3,5,6,2,7,1,4,8,9,7,93,根据,i(x,;z),值算出的不同,z,值的,(h;z),值,计算的各指示值,i(x,;z),94,指示克里金,线性估计量,i,*,(x;z)=,其中，,x,(=1,2,，,n),点处的,Z(x,),值已知，,(x;z)(=1,2,，,n),为,IK,权系数,(00111),(00001),(01111),(00011),（普通指示克里金）,95,C,I,(h;z),为指示协方差，,为指示变差函数,是拉格朗日乘数,有多少个门限值,z,，就有多少个,IK,方程组,IK,方程组：,96,从,IK,方程组中解出,i,*,(x;z)=,求出,i(x;z),的估计值,i,*,(x;z),实际上是,I(x;z),的条件数学期望,EI(x;z)|(n),的估计值。,EI(x;z)|(n)=0PI(x;z)=0|(n)+1PI(x;z)=1|(n),=PI(x;z)=1|(n),=PZ(x,）,z)|(n),=Fx;z|(n),为后验概率分布函数，或条件概率分布函数，等于在已知,n,个信息样品的条件下，概率,P Z(x)z,的大小。,i,*,(x;z)=,EI(x;z)|(n),即：,i,*,(x;z)=,EI(x;z)|(n),又：,=PZ(x,）,z)|(n),97,(00111),(00001),(01111),(00011),Z,*,(x)=,各类的,局部均值,各类的,后验概率,估计的指示值,98,Sunbeson,金矿，有,3500,个样品，边界品位值分别取,z=0.009,，,0.015,，,0.020,，,0.033,克,/,吨，用这,4,个边界品位可将金矿化分为,5,类：,应用,IK,法估计出以下的条件概率：,P,*,x(1)|(n)=18.2%,P,*,x(2)|(n)=6.5%,P,*,x(3)|(n)=0%,P,*,x(4)|(n)=57.8%,P,*,x(5)|(n)=17.5%,各概率的总和为,1,例：,99,各类的,(,局部,),均值为：,Z(x)|x(1),*,=0.004g/t,Z(x)|x(2),*,=0.0125g/t.,Z(x)|x(3),*,=0.020g/t,Z(x)|x(4),*,=0.0276g/t.,Z(x)|x(5),*,=0.125g/t,按,IK,法可得总平均品位的估计值为,Z,*,(x)=,x(k),(n)Z(x),x(k),*,=0.0393683g/t.,100,A,(100),B,(010),A,(100),C,(001),离散变量的指示克里金估计,分别对各类型进行指示克里金估计，得出各类型的概率估计。,（最大概率的类型,即为待估点的类型）,101,作为一种非参数统计方法，指示克里金估计方法在处理特高值和特低值的分布方面，具有明显的优势。,特色一：,(00111),(00001),(01111),(00011),离散变量（类型变量）,连续变量,102,应用,贝叶斯理论,，可综合各种软信息（与硬信息一起）进行指示克里金估计。,（如地震信息、试井解释、地质推理和解释等）,硬信息：能准确标定所研究地质,变量取值的信息。,软信息：不能准确标定所研究地质,变量取值，而只能提供其,概率分布的信息。,特色二：,103,使用贝叶斯理论，将先验的局部累计分布函数更新为后验的条件局部累计分布函数。,整体,先验概率,硬指示数据,软指示数据,A,(100),B,(010),A,(100),C,(001),104,l,估计的无偏性,克里金方法的优点,l,反映了变量的空间结构性,l,能得到估计精度,105,（,1,）克里金插值为局部估计方法，,对估计值的整体空间相关性考虑不够,，,它保证了数据的估计局部最优，却不能保证数据的总体最优,，因为克里金估值的方差比原始数据的方差要小。因此，当井点较少且分布不均时可能会出现较大的估计误差，特别是在井点之外的无井区误差可能更大。,克里金方法的局限性,106,（,2,）,克里金插值法为光滑内插方法,，为减小估计方差而对真实观测数据的离散性进行了平滑处理，虽然可以得到由于光滑而更美观的等值线图或三维图，但,一些有意义的异常带也可能被光滑作用而“光滑”掉了,。所以，有时，克里金方法被称为一种“移动光滑窗口”。,（用于,CCDF,的求取，应用于随机建模）,(3),107,