第四、五讲 地质统计学理论基础
地质学基础PPT课件全文

2024/8/31
8
第89页/共45页
第一节 地球概况
二、地球的物理性质
(一)地球的密度和重力
布格校正
2024/8/31
地壳均衡说
9
第第190页页/共/共4455页页
第一节 地球概况
二、地球的物理性质 (二)地磁
2024/8/31
10
第110页/共45页
第一节 地球概况
二、地球的物理性质 (二)地磁
2024/8/31
26
第267页/共45页
第二节 地球的结构
二、地球的内部圈层 (一)地壳
2.地壳的厚度和结构
陆壳:约占地壳面积的 1/3多一点,陆壳具有 明显的双层 结构,即存 在上、下地壳。
洋壳:位于海洋之下, 约占地壳面积2/3少一 点,其上为 约4km厚的 海水,洋壳缺失上地壳
2024/8/31
第一节 地球概况 一、地球的形状和大小
大地水准体” (Geoid)
2024/8/31
1
第12页/共45页
第一节 地球概况 一、地球的形状和大小
2024/8/31
2
第23页/共45页
第一节 地球概况
一、地球的形状和大小
地球赤道半经(α):6378137m 地球极半经(с):6356752m
2024/8/31
用以计算岩石的年龄的公式为:
t=ln(1+D/N)/λ
λ为衰变常数;D为子体同位素含量;N为母体同位素含量。 目前用放射性同位素方法测得地球上最古老岩石的年龄为40- 43亿年; 对来自外星球的陨石及月岩的测定,获得的最大年龄为45- 47亿年。 据此,确定地球的年龄为至少有45亿年。
2024/8/31
地质统计学

地统计(Geostatistics)又称地质统计,是在法国著名统计学家G. Matheron大量理论研究的基础上逐渐形成的一门新的统计学分支。
它是以区域化变量为基础,借助变异函数,研究既具有随机性又具有结构性,或空间相关性和依赖性的自然现象的一门科学。
凡是与空间数据的结构性和随机性,或空间相关性和依赖性,或空间格局与变异有关的研究,并对这些数据进行最优无偏内插估计,或模拟这些数据的离散性、波动性时,皆可应用地统计学的理论与方法。
地统计学与经典统计学的共同之处在于:它们都是在大量采样的基础上,通过对样本属性值的频率分布或均值、方差关系及其相应规则的分析,确定其空间分布格局与相关关系。
但地统计学区别于经典统计学的最大特点即是:地统计学既考虑到样本值的大小,又重视样本空间位置及样本间的距离,弥补了经典统计学忽略空间方位的缺陷。
地统计分析理论基础包括前提假设、区域化变量、变异分析和空间估值。
第一章品位与储量计算第一节概述投资一个矿床开采项目,首先必须估算其品位和储量。
一个矿床的矿量、品位及其空间分布是对矿床进行技术经济评价、可行性研究、矿山规划设计以及开采计划优化的基础,是矿山投资决策的重要依据。
因此,品位估算、矿体圈定和储量计算是一项影响深远的工作,其质量直接影响到投资决策的正确性和矿山规划及开采计划的优劣。
从一个市场经济条件下的矿业投资者的角度看,这一工作做不好可能导致两种对投资者不利的决策:(1)矿体圈定与品位、矿量估算结果比实际情况乐观,估计的矿床开采价值在较大程度上高于实际可能实现的最高价值,致使投资者投资于利润远低于期望值,甚至带来严重亏损的项目。
(2)与第一种情况相反,矿床的矿量与品位的估算值在较大程度上低于实际值,使投资者错误地认为在现有技术经济条件下,矿床的开采不能带来可以接受的最低利润,从而放弃了一个好的投资机会。
然而,准确地估算出一个矿床的矿量、品位绝非易事。
大部分矿体被深深地埋于地下,即使有露头,也只能提供靠近地表的局部信息。
地质统计学

第一章绪论一、历史背景与产生地质统计学是二十世纪六七十年代发展起来的一门新兴的数学地质学科的分支。
它开始主要是为解决矿床从普查勘探、矿山设计到矿山开采整个过程中各种储量计算和误差估计问题而发展起来的。
它是由法国著名学者G. 马特隆教授于1962年创立的。
其核心即所谓的“克立格”。
它是一种无偏的最小误差的储量计算方法。
该方法按照样品与待估块段的相对空间位置和相关程度来计算块段品位及储量,并使估计误差为最小。
这是南非采矿工程师D. G. Krige 根据南非金矿的具体情况与1952年提出的,故命名为克立格法。
后来法国学者G. 马特隆(Matheron)对克立格提出的方法进行研究,认为克立格提出的方法是在考虑了空间分布特征的基础上,合理地改进了统计学,是一种传统方法与统计学方法结合起来的新方法。
同时为了解决具二重型(结构型与随机性)的地质变量的条件下使用统计方法的问题。
马特隆教授提出了区域化变量的概念(Regionalized Variable),从而创立了地质统计学。
根据地质统计学理论,地质特征可以用区域化变量的空间分布特征来表征。
而研究区域化变量的空间分布特征分布的主要数学工具是变差函数(Variogram)。
到七十年代中后期,马特隆的学生JOURENL等在研究其它地质变量的基础上,认为某些地质变量并不是一成不变的,而是有一定波动的,这样使用克立格法就不能很好再现地质变量的分布特征。
因此他们采样模拟的方法,将克立格估计的离散方差的波动性模拟出来,从而产生了随机模拟法。
因此,从二十世纪八十年代以来,地质统计学分为两派:一派以法国的马特隆教授等人为主,仍致力于克立格估计的研究;一派以美国JOURENL等人为主,主要致力于随机模拟方法的研究。
地质统计学的产生是在经典统计学的基础上发展起来的。
在此前,为了反映地质变量的空间变化性,一些地质学家曾经使用一些经典的概率统计方法来研究地质变量。
但由于地质变量并不是纯粹的随机变量,因此,直接用简单的统计方法解决复杂的地质问题,有一定的局限性。
地质统计学及其应用介绍PPT

第二节 地质统计学的研究现状及优点
一、研究现状
(1)线性平稳地质统计学
(2)线性非平稳地质统计学 (3)条件模拟 (4)平稳非线性地质统计学 (5)储量参数确定
1.初步形成了一直完整的理论体系
基本概念:区域化变量 基本工具:变差函数 基本假设:二阶平稳假设和本征假设 基本公式:估计方差,离散方差 基本方法:克里格法
问题2:品位空间变化问题:矿化的空间结构。如:走向上变化小 ,倾向变化大,权值不一样。
问题3:矿化强度的空间变化问题:离散度。这与问题2相关联,离 散度是衡量经济开采可行度的重要因素。
问题4:缺乏估计精度的方1:不考虑样品的空间分布
40
40
30
30
20
20
地 质 统 计 学
地质统计学简介及其应用ppt课件

(Simple Kriging)
的 2、普通克里金
(Ordinary Kriging)
几 3、非稳态克里金 (Nostationary Kriging)
种
克 4、内在趋势克里金 (Universal Kriging)
里
(泛克里金)
金
5、外在趋势克里金 (Kriging with an External Drift)
头
尾
滞后距(Lag)
经营者提供商品或者服务有欺诈行为 的,应 当按照 消费者 的要求 增加赔 偿其受 到的损 失,增 加赔偿 的金额 为消费 者购买 商品的 价款或 接受服 务的费 用
散 点 图
经营者提供商品或者服务有欺诈行为 的,应 当按照 消费者 的要求 增加赔 偿其受 到的损 失,增 加赔偿 的金额 为消费 者购买 商品的 价款或 接受服 务的费 用
经营者提供商品或者服务有欺诈行为 的,应 当按照 消费者 的要求 增加赔 偿其受 到的损 失,增 加赔偿 的金额 为消费 者购买 商品的 价款或 接受服 务的费 用
4、内在趋势克里金
(Universal Kriging)
经营者提供商品或者服务有欺诈行为 的,应 当按照 消费者 的要求 增加赔 偿其受 到的损 失,增 加赔偿 的金额 为消费 者购买 商品的 价款或 接受服 务的费 用
经营者提供商品或者服务有欺诈行为 的,应 当按照 消费者 的要求 增加赔 偿其受 到的损 失,增 加赔偿 的金额 为消费 者购买 商品的 价款或 接受服 务的费 用
几种克里金算法之间的关系
简单克里金 (Simple)
非稳态克里金 (Nostationary) 同位协克里金(Collocated Cokriging)
地质统计学讲义.doc

地质统计学讲义第1章地质统计学的发展历史和现状1.1地质统计学的发展历史地质统计学是根据英文单词Geostatistics的字面意思翻译过来的,从词源学上讲,按照韦氏(N .Webster)大词典对于“geo”(地球、土地)和“Statistics”(统计学)两词的释义,地质统计学(Geostatistics)的定义便是:“关于取自地球的大量数据的收集、分析、解释和表达的一个数学分支”。
就矿山地质统计学的内容范围来说,这一定义是十分恰当的。
地质统计学包含经典统计学与空间统计学,其重点是地球状况,也就是说着重于地质特征的分析。
按其基本原理可定义为:地质统计学是以区域化变量理论为基础,以变异函数为主要工具,研究那些在空间分布上既有随机性,又有结构性的自然现象的科学。
早在上世纪10年代里,传统的统计学方法就已用于分析地质数据。
在地质矿产方面最初也是利用传统的统计学作为分析数据的工具,直到上世纪40年代后期,当南非统计学家H.S西奇尔(Sichel)判明南非各金矿的样品品位呈对数正态分布以后,才真正确立了地质统计学的开端。
1951年,南非的矿山工程D.G.克立格(Daniel Krige)在H.S西奇尔研究的基础上提出一个论点:“可以预计,一个矿山总体中的金品位的相对变化要大于该矿山某一部分中的金品位的相对变化”。
换句话说,以较近距离采集的样品很可能比以较远距离采集的样品具有更近似的品位。
这一论点是描述在多维空间内定义的数值特征的空间统计学据以建立的基础。
到上世纪60年代,才认识到需要把样品值之间的相似性作为样品间距离的函数来加以模拟,并且得出了半变异函数。
法国概率统计学家马特隆(Matheron)创立了一个理论框架,为克立格作出的经验论点提供了精确而简明的数学阐释。
马特隆创造了一个新名词“克立格法”(Kriging),藉以表彰克立格在矿床的地质统计学评价工作中所起到的先驱作用。
即1962年,马特隆在克立格和西奇尔研究的基础上,将他们的成果理论化、系统化,并首先提出了区域化变量(Regionalized variable)的概念,为了更好地研究具有随机性及结构性的自然现象,提出了地质统计学(Geostatistics)一词,发表了《应用地质统计学》,该著作的出版标志着地质统计学作为一门新兴边缘学科而诞生。
地质统计学基本原理

Z(x 差h)的方差之半定义为区域化变量 的Z(变x)差函数,记为
(x, h)
(x, h) 1 Var[Z (x) Z (x h)]
2
变差函数定义
• 定义:在任一方向 a ,相距 | h |的两个区域 化变量 Z(x) 和 Z(x h) 的增量的方差的一半。
• 公式: (h) 1 E[Z (x) Z (x h)]2
几点注意内容
• 变差函数参数
• 块金值:块金值越小,距离越近的点越重要,这样会导 致权值的变化范围变大(从负值到大于1的值变化),使 数据出现异常。块金值越大,估值结果越平滑。
当时h 0,上式变成:
Var[Z(x)] C(0) x
即它有有限先验方差。
本征假设
当区域化变量Z(x) 的增量 Z(x) Z(x h) 满足下列两个条 件时,称该区域化变量满足本征假设: (1)在整个研究区内,区域化变量Z(x的) 增量 Z(x) Z(x 的h)
期望为0: E[Z(x) Z(x h)] 0 x,h
滞后距
实验变差函数计算实例
• 相距为200米的样本点对。
实验变差函数计算实例
• 滞后距为200米的变差函数值。
变差函数计算实例
• 变差函数图:滞后距200米的变差函数点
变差函数
20 18 16 14 12 10
8 6 4 2 0
0
100
200
300
400
500
滞后距
变差函数计算实例
• 变差函数图:滞后距300米、400米的变差函数点
几何各向异性
• 基台值相同 • 变程不同
在不同的方向具有相同的变异程 度(基台值相同)但具有不同的 连续程度(变程不同)为几何各 向异性。
地质统计学基础知识

信息管理学院 王玉兰 E-mail:wyl@,wang_wyl@ Tel:84073385 (o)
概率论基础
• 随机事件:
– 在一定的条件下可能发生也可不发生的事件; – 结果具有不确定性的事件;
• 统计概率(随机事件的概率):随机事件出现的可能性;N次 重复试验中,事件A出现了n次,则事件A发生的概率n/N • 随机事件之间的关系:包含、相等、互斥、互逆; • 概率的运算:
– – – – – 加法定理:P(A+B)=P(A)+P(B) 乘法定理:P(AB)=P(A|B)P(B)=P(B|A)P(A)--P(A|B)=P(AB)/P(B) 全概率公式:P(A)=∑P(Hi)P(A|Hi) Hi:独立事件完备群 Bayes公式:P(Hi|A)= P(Hi)P(A|Hi)/P(A)
随机函数、随机过程与随机场
• 随机函数:
– 具有n个参数的随机变量族Z(x1,x2,…,xn;w); – 或随机函数所有实现的集合
• 随机过程:
– 只有一个(时间)变量的随机函数Z(t;w) – 对每个固定t, Z(t;w)是随机变量; – 或Z(于多个变量的随机函数Z(x1,x2,…;w)
统计推断(总体估计)
• 总体与样本、个体及样本观测值:
– – – – – 总体:研究对象的全体; 样本:总体的一部分; 个体:组成总体的基本单元; 样品:从总体中按某种方式随机抽取的个体; 样本观测值:对样本中所有样品的观测、测量或分析值 的集合。
• 理论分布与经验分布:
– 理论分布:总体的客观存在的分布特征; – 经验分布或实验分布:由从总体中随机抽取的样本所对 应的分布;
统计推断(总体估计)
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
当空间一点x固定之后,Z(x)(表示x点处的矿石 品位)就是一个随机变量,体现了其随机性。
在空间两个不同点x及x+h(此处h也是个三维向量
(hu,hv,hw)。它的模
h
h h h 2 2 2
u
v
w
表示x点与(x+h)点
的距离)处的品位Z(x)与Z(x+h)具有某种程度的
相关性,这就体现了其结构性的一面。
6/88
区域化变量的属性
• 1、空间局限性 • 2、连续性 • 3、异向性 • 4、相关性 • 5、叠加性
7/88
1、空间局限性
区域化变量被限制 于一定的空间,该 区间称为区域化变 量的几何域。例如, 矿体的范围,油藏 的范围,断块的范 围等都可以看成是 区域化变量的几何 域。
Z(xi)
Z(xk) Z(xj)
第三章 区域化变量与变差函数
区域化变量及其基本特征 变差函数的定义 变差函数曲线 变差函数的理论模型 变差函数的结构分析
1/88
第一节 区域化变量
区域化变量(Regionalized Variable) 是地质统计学研究的对象,它是一种在空
间上具有数值的实函数(G Matheron),也就 是说,它在空间的每一个点取一个确定的数 值,即当由一个点移到下一个点时,函数值 是变化的
2 (x, h) Var[Z(x) Z(x h)]
E[Z (x) Z (x h)]2 {E[Z (x)] E[Z(x h)]}2
Z(x1)
观测前是一个随机场,
Z(x2)Z(x7) NhomakorabeaZ(x3) Z(x6) Z(x8)
依赖于坐标(xu,xv,xw)
Z(x4) Z(x5)
观测前随机变量的集合
观测后是一个空间点 函数,在具体的坐标 上有一个具体的值
z(x1)
z(x2)
z(x7)
z(x3) z(x6) z(x8)
z(x4) z(x5)
确定性
当h=0时,上式变为 Var[Z(x)]=C(0) x 此式说明:方差函数也存在,且为常数C(0)
16/88
本征假设(内蕴假设)
在实际应用中,有时连二阶平稳假设的要求也不能满足, (如协方差函数不存在或方差函数不存在等)。这时,可 以再放宽条件,得到本征假设 当区域化变量Z(x)的增量[ Z(x)- Z(x+h)]满足下列两条 件时,称其满足本征假设: 1、在整个研究区内有
4/88
观测后实数(实现)的集合
区域化变量举例
在地质、采矿领域中许多变量都可看成是 区域化变量:资源储量、储层厚度、地形 标高、矿石内有害组分含量、岩石破碎程 度、孔隙度、渗透率、泥质含量等。有的 是二维的,有的是三维的。区域化变量正 是地质统计学研究的对象。
5/88
区域化变量的功能
由于区域化变量是一种随机函数,因而能同时反 映地质变量的结构性与随机性。
P{Z(x1 h) z1, Z(x2 h) z2,, Z(xn h) zn}
Fx1h,x2 h,,xn h (z1, z2 ,, zn )
n,h,x1, x2,, xn
这种假设要求的条件太强了,实际上很难满足。在地质统计 学中,只需要假设Z(x)的一阶、二阶矩存在且平稳就够了, 也就是二阶平稳假设
17/88
三、变差函数
变差函数是地质统计学所特有的基本工具,它 既能描述区域化变量的结构性变化,又能描述 其随机性变化。是地质统计学计算的基础。
18/88
1、变差函数的定义
设区域化变量Z(x)定 义在一维数轴x上, 把Z(x)在x,x+h两点 处的值之差的方差之 半定义为Z(x)在x轴方 向上的一维变差函数, 记为:
8/88
2、连续性
不同的区域化变量具有不 同程度的连续性,这种连 续性是通过区域化变量的 变差函数来实现的
9/88
3、异向性
区域化变量 在各个方向 具有不同的 性质时称为 各向异性, 否则称为各 向同性。在 地质上,各 向异性是绝 对的,而各 向同性是相 对的
10/88
区域化变量 在一定的范 围之内呈现 一定程度的 空间相关性, 当超出这一 范围之后, 相关性变弱 以至消失
15/88
二阶平稳假设
当区域化变量Z(x)满足下列两条件时,称其为二阶平稳的: 1、在整个研究区内Z(x)的数学期望均存在,且等于常数,即
E(Z(x)=m(常数) x 2、在整个研究区内Z(x)的协方差函数存在且平稳(即只依赖 于基本步长h,而与x无关,即:
Cov{Z(x),Z(x+h)}=E[Z(x)Z(x+h)-E[Z(x)]E[Z(x+h)] =E[Z(x)Z(x+h)]-m2 =C(h) x h
2/88
区域化变量图示
区域化变量是以空间点x 的三个直角坐标(xu,xv,xw) 为自变量的随机场
Z(xu,xv,xw)=Z(x) 当对它进行一次观测后观 测后,就得到了他的一个 实现z(x),它是一个普通的 三元实值函数或空间点函 数
w
v Z(xu,xv,xw) u
3/88
区域化变量的两重性
随机性
但是在许多情况下,在某一个位置(u)只有一个样品,那 么z(u)是已知的。也就是说,区域化变量的取值是唯一 的,不能重复,为了克服这个困难,提出了如下的平稳 假设。
14/88
平稳假设
假设区域化变量Z(x)的任意n维分布函数均不因空间点x发生 位移h而改变,即:
Fx1,x2,,xn (z1, z2,, zn ) P{Z (x1) z1, Z(x2 ) z2,, Z(xn zn}
4、相关性
11/88
5、叠加性
对于任一区域化变量而言,特殊的 变异性可以叠加在一般的规律之上
12/88
二、地质统计学的若干基本假设
平稳假设 内蕴假设
13/88
平稳假设(stationary assumption)
任何统计的推断,不论是单变量的累积概率分布函数(cdf) 或是它的任何阶矩(均质、方差),还是多变量的cdf 及其任何阶矩(协方差),都需要重复取样。
E[Z(x)- Z(x+h)]=0 , x ,h 若Z(x) ,x存在,则此条件等于
E[Z(x)]=E[Z(x+h)]=m(常数) x ,h 2、增量[ Z(x)- Z(x+h)]的方差函数存在且平稳(即方差 函数不依赖于x)
Var[Z(x)- Z(x+h)]= E[Z(x)- Z(x+h)]2-{E[Z(x)- Z(x+h)]}2 = E[Z(x)- Z(x+h)]2=2(x,h)=2 (h) x ,h