地质统计学(2)_区域化变量的理论cjg2011.pptVIP

第二章 地质统计学理论基础 一、随机场与区域化变量 1．定义：以空间点x的三个直角坐标xu,xv,xw为自变量的随机场 Z(xu,xv,xw）=Z(x）称为一个区域化变量。 [区域化变量具有两重性]： 观测前，将Z(x）看作随机场；观测后，将Z(x）看作一个普通的三元实值函数。即空间点函数，一次观测后，就得到它的一个实观Z(x)。 第一节 区域化变量的理论 2.功能 能同时反映地质变量的结构性与随机性。 ①当空间点x固定后， Z(x）即为一个随机变量； ②x与x+h两点处的Z(x）具有某种程度的相关性（因随机场有相关函数R（x，x+h））即为一个随机变量； 3.物理学或地质学特征 ①空间局限性；②不同程度的连续性；③不同类型的各向异性。 1. 协方差函数 若Z(x）是随机场，在空间两点x和x+h 处两个随机变量Z(x）和 Z（x+h）的二阶中心混合矩 称为随机场的Z(x）自协方差函数，简称协方差函数。一般地讲，它是依赖于点x和向量h 的函数。 特殊地：当h =0时， 就等于方差函数: 当其不依赖于x时简称方差，故有: 二、协方差函数与变差函数 2. 变差函数与变差图 假设空间点x只在一维的x轴上变化，我们把区域化变量Z（x）在x，x+h两点处的值之差的方差之半定义为Z（x）在x 轴方向上的变差函数，记为r (x,h),即: 在二阶平稳和本征假设条件下： 于是变差函数的计算公式变为： x x+h 基本公式 在二维、三维情况下定义时，以一维变差函数为基础，需考虑各向异性，结构套合等问题。 当r（x，h）与x的取值无关时，r（x，h）只依赖与h（滞后、间隔、步长），则可将r（x，h）写成r（h），此时以h为横坐标，r（h）为纵坐标作出图形谓之变差图。 变程a 拱高C 基台值C+C0 块金常数C0 h r (h) [问题]：由数理统计知：要估计变差函数值 就要估计数学期望值 这必须有若干对Z（ x ）和Z（ x+h ）的值才可通过求 平均数的办法来估计上述数学期望。而这在实际地质，采矿工作中是不可实现的，因为不可能恰在空间同一点上重复直接取得二个样品。这就使统计陷入困境。需借助假设来解决。 三．平稳假设与本征假设 两个重要的假设条件: 1. 平稳假设 2. 本征假设 1．? 平稳假设 ① 严格的平稳假设 区域化变量Z（x）的任意n维分布函数不因空间点 x发生位移h而改变。 即： 这种要求是Z（x）的各阶矩存在，且平稳，这在实际中不能满足，且不好验证。所以实用上采用的只需一、二阶矩且平稳就够了。→ 二阶平稳（弱平稳）。 ? ② 二阶平稳假设 满足下列两个条件 1）整个研究区内，Z（x）的数学期望存在，且等于常数， 2）整个研究区内，Z（x）的协方差函数存在且平稳（即只依赖于滞后h，而与x无关） 特殊地：当h=0时 =C（0） 即方差存在且为常数。当上述条件仍不能满足时，条件进一步放宽，导致本征假设。 2. 本征假设 区域化变量Z（x）的增量[Z（x）- Z（x+h）]满足下列两个条件： 1）? 在整个研究区内有:  2）增量[Z（x）– Z（x+h）]的方差函数存在且平稳（不依赖于x）即： = =2r(h), 即Z（x）的变差函数存在且平稳。 3 .二阶平稳假设与本征假设的比较 总的结论：二阶平稳假设较强，本征假设较弱 1）? 由二阶平稳假设的第一个条件可推出本征假设条件一。 如：设 y为一服从柯西分布的随机变量，其概率密度为 则：