地质统计学课件.ppt

地质统计学一、地质统计学的发展历史和现状地质统计学的发展历史和现状 什么是地质统计学？地质统计学(Geostatistics)包含经典统计学与空间统计学，按其基本原理可定义为：地质统计学是以区域化变量理论为基础，以变异函数为主要工具，研究那些在空间分布上既有随机性，又有结构性的自然现象的科学。地质统计学诞生过程上世纪40年代后期，当南非统计学家H.S西奇尔（Sichel）判明南非各金矿的样品品位呈对数正态分布以后，才真正确立了地质统计学的开端。1951年，南非的矿山工程D.G.克立格Daniel Krige）在H.S西奇尔研究的基础上提出一个论点：“可以预计，一个矿山总体中的金品位的相对变化要大于该矿山某一部分中的金品位的相对变化”。换句话说，以较近距离采集的样品很可能比以较远距离采集的样品具有更近似的品位。这一论点是描述在多维空间内定义的数值特征的空间统计学据以建立的基础。到上世纪60年代，才认识到需要把样品值之间的相似性作为样品间距离的函数来加以模拟，并且得出了半变异函数。法国概率统计学家马特隆（Matheron）创立了一个理论框架，为克立格作出的经验论点提供了精确而简明的数学阐释。马特隆创造了一个新名词“克立格法”（Kriging），藉以表彰克立格在矿床的地质统计学评价工作中所起到的先驱作用。即1962年，马特隆在克立格和西奇尔研究的基础上，将他们的成果理论化、系统化，并首先提出了区域化变量（Regionalized variable）的概念，为了更好地研究具有随机性及结构性的自然现象，提出了地质统计学（Geostatistics）一词，发表了应用地质统计学，该著作的出版标志着地质统计学作为一门新兴边缘学科而诞生。地质统计学开始进入了学术界。在法国枫丹白露成立了地质统计学中心（Centre de Geostatistiques），培养了一大批学员，不仅为地质统计学的研究而且为它的传播起到了巨大的作用。地质统计学是在1977年由美国福禄尔采矿金属有限公司（Flour Mining&Meta Incorporation）H.M.Parker博士随美中贸易全国委员会矿业代表团来华访问，传入我国，继而得到进一步的发展。1989年11月召开的全国第一届地质统计学学术讨论会，地质统计学在我国的发展进入了一个新的阶段，理论研究更加深入，涉及的方法原理更加广泛。地质统计学已经被广泛地承认是矿床评价的必要部分，在我国已经认可用地质统计学对矿床进行评价的地质报告。1、区域化变量理论区域化变量理论 区域化变量G.马特隆定义区域化变量是：一种在空间上具有数值的实函数，它在空间的每一个点取一个确定的数值，即当由一个点移到下一个点时，函数值是变化的。从地质及矿业角度来看，区域化变量具有如下性质：（1）空间局限性：即它被限制在一个特定的空间（如一个矿体内）；该空间称为区域化的几何域；区域化变量是按几何支撑定义的。（2）连续性：不同的区域化变量具有不同的连续性，这种连续性是通过相邻样品之间的变异函数来描述的。（3）异向性：当区域化变量在各个方向上具有相同的性质时称各向同性，否则称各向异性。（4）相关性：一定范围内、一定程度上的空间相关性，当超出这一范围后相关性减弱以至消失。（5）对于任一区域化变量而言，特殊的变异性可以叠加在一般规律之上。相关关系相关系数(计算公式)样本相关系数的计算公式相关系数(取值及其意义)-1.0+1.00-0.5+0.5完全负相关完全负相关完全负相关完全负相关无线性相关无线性相关无线性相关无线性相关完全正相关完全正相关完全正相关完全正相关负相关程度增加负相关程度增加负相关程度增加负相关程度增加r正相关程度增加正相关程度增加正相关程度增加正相关程度增加 协方差函数(Covariance)内蕴假设(intrinsic assumption)估计方差估计方差(Estimation variance)1 估计方差估计方差(Estimation variance)3 离差方差离差方差(Dispersion variance)2、变异函数及结构分析、变异函数及结构分析 为表征一个矿床金属品位等特征量的变化，经典统计学通常采用均值、方差等一类参数，这些统计量只能概括该矿床中金属品位等特征量的全貌，却无法反映局部范围和特定方向上地质特征的变化。地质统计学引入变异函数这一工具，它能够反映区域化变量的空间变化特征相关性和随机性，特别是透过随机性反映区域化变量的结构性，故变异函数又称结构函数。变异函数及变异曲线 变异函数及变异曲线变异函数与协方差之间的关系 存在趋势的变异函数具有空穴效应的变异函数 变异函数在原点处的性状 变异函数的理论模型Co-块金常数a变程Co+C基台三种有基台值标准模型比较具有空穴效应的变异函数 无基台值标准模型 变异函数结构分析 几何异向性结构的套合带状异向性结构的套合比例效应 相对变异函数改正 变异函数的套合 普通克立格法(Ordinary Kriging)普通克立格方程组正态分布的误差图示 x95%95%的概率的概率的概率的概率 -1.96-1.96 x x +1.96+1.96 x x99%99%的概率的概率的概率的概率 -2.58-2.58 x x +2.58+2.58x x90%90%的概率的概率的概率的概率 -1.65-1.65 x x +1.65+1.65 x x泛克立格法(Universal Kriging)指示克立格法(Indicator Kriging)指示克立格法(Indicator Kriging)协同克里格法（协同克里格法（Co-Kriging）协同克里格法（协同克里格法（Co-Kriging）协同克里格法（协同克里格法（Co-Kriging）协同克里格法（协同克里格法（Co-Kriging）