地质统计学资源量估算

2 样品等长度组合
划分方法 划分范围 组合样长 最小样长 品位阈值 品位特高值 同时划分加密样品 导入矿体边界约束样品
3 组合样品数据分析
数据导入：加载离散数据点 将.cbs文件导入系统；选择进行储量估算的矿
体文件。
3 数据分析
直方图分析 了解组合样数据分布形态、是否服从正态分布、
位置图分析 了解组合样数据位置分布是否有误、数据是否
存在趋势。
• 数据分布正常 • 数据不存在趋势
4 实验变差函数计算及拟合
数据处理
正态变换 对数变换 指示变换
变差函数模型 变差函数类型
4 实验变差函数计算及拟合
实验变差函数计算
变差函数参数 变差函数图形 各向异性角度 变差函数表面图
1.矿块模型显示与查询
可以根据某个矿体内矿块内某种元素的品位大小区间，将矿块按不同的颜色 进行显示
矿块(单元素)显示
8 块体品位模型应用
1.矿块模型显示与查询
可以根据某个矿体内矿块内某种元素的品位大小区间，将矿块按不同的颜色 进行显示
显示部分品位区间的矿块
多元素集成显示
8 块体品位模型应用
2.矿块模型切割
5 克里格品位估值（示例）
克里格品位估值结果
三维图形显示，用户可以查询任意块体模型属性信息。
6 资源量计算
用户输入矿石比重，更新数据库，重新计算储量。
6 资源量分级
两种方式：一种根据估值过程中用到的样品个数、工程个 数以及搜索次数进行区分，用户可修改分级名称以及区分 范围；一种根据克里格估值方差进行分级。
最佳步长大小。
• 最佳步长为20m
4 实验变差函数计算及拟合
变差函数分析
计算变差函数（短距离变差函数） 点击上一步选择“短距离变差函数”，双击列表
控件设置步长大小，一般情况下步长大小为组合样样 长。
• 块金值为0.3
4 实验变差函数计算及拟合
变差函数分析
计算变差函数（计算三个相互垂直的变差函数）
结果保存到数据库中
7 报表生成
对数据库中矿体储量进行统计，生成统计报表。
8 块体品位模型应用
1.矿块模型显示与查询 2.矿块模型切割 3.基于矿块模型计算采空区动用储量
8 块体品位模型应用
1.矿块模型显示与查询
矿块模型可以与原始勘探工程数据、 矿体约束信息组合显示
选择数据源
空块显示
8 块体品位模型应用
块段离散化个 数
5 克里格品位估值
交叉验证
散点图 直方图 误差平均值 标准中误差 统计量比较
5 克里格品位估值（示例）
数据变换及模型选择 选择合适的克里格估值类型，并完成变差函数
的导入。
• 克里格类型：普通 克里格估值算法
• 导入变差函数文件
5 克里格品位估值（示例）
搜索邻域设置 通过一个搜索椭球定义搜索邻域，一般情况下
选中“计算第二、第三方向”输入方位角和倾角，点击“输入
第二、第三方向参数”按钮，计算三个相互垂直方向的变差函
数。
• 三个相互垂直方
向的变差函数
4 实验变差函数计算及拟合
变差函数分析
拟合变差函数 通过滑动条设置变差函数的块金值、基台值和变程。
主轴方向 • 变差函数模型类
型：球状模型 • 块金值：0.3 • 基台值：2.185 • 变程：97.8m
• 误差平均值：0.0068
• 标准中误差：1.062
• 小于中误差样品比例： 89.64%
5 克里格品位估值（示例）
克里格品位估值结果
克里格估值结果包括：块体模型中心点坐标、块体模型样 品品位值、品位估计方差，结果保存到数据库中。
• 数据库中包含： • 块体模型参数表 • 变差函数参数表 • 空块模型表 • 块体模型品位表
4 实验变差函数计算及拟合
变差函数分析
拟合变差函数 通过滑动条设置变差函数的块金值、基台值和
Hale Waihona Puke Baidu变程。
半轴方向 • 变差函数模型类
型：球状模型 • 块金值：0.3 • 基台值：2.185 • 变程：80.6m
4 实验变差函数计算及拟合
变差函数分析
拟合变差函数
通过滑动条设置变差函数的块金值、基台值和
4 实验变差函数计算及拟合 变差函数分析
计算变差函数（全方向变差函数）
步长总间距
• 步长数目*步长大小=区域长度的一半 • 方位角，倾角均为0度。 • 方位角容差，倾角容差均为90度。
4 实验变差函数计算及拟合 变差函数分析
计算变差函数（全方向变差函数） 点击“计算变差函数”按钮，观察图形，寻找
主轴：基台值为 2.185， 变程为97.8.1m
半轴：基台值为 2.185， 变程为80.6m
次轴：基台值为 2.185， 变程为31.8m
5 克里格品位估值
克里格估值类型
简单克里格 普通克里格 泛克里格 指示克里格
距离密次方反比
5 克里格品位估值
搜索椭球设置
搜索椭球半径 旋转角度 多次搜索 搜索类型
8 块体品位模型应用
2.矿块模型切割
不同方向的切割后 的模型整合显示
8 块体品位模型应用
3.基于矿块模型计算采空区动用储量
根据实际测量数据，在三维视图中显示采空区边界范围 根据采空区边界范围和采空区所属采掘面间关系自动生成采 空区三维模型 基于采空区三维边界，手工交互生成采空区三维模型功能 显示采空区与矿体表面模型相交结果并获取采空区动用储量 显示采空区与矿体块体模型相交结果并获取采空区动用储量 可对求交后的矿体表面模型进行联合显示与信息查询 可对求交后的矿块模型进行联合显示与信息查询
4 实验变差函数计算及拟合
变差函数拟合
自动拟合变差函数 变差函数图形 理论变差函数参数 各向异性结构套合
4 实验变差函数计算及拟合
变差函数分析
数据变换及模型选择
根据数据分析结果 进行数据变换，特异值 处理，选择；通过计算 全方向变差函数求得最 佳步长，通过计算短距 离变差函数求得块金值， 通过定向变差函数计算 求得各向异性角度
在矿块模型显示的基础上，可以对该模型进行不同方向的切割 ，获得切割面上的相关信息 切割方向可以为：勘探线方向、纵剖面方向、中段方向 相应块体模型剖切结果可以保存为矢量数据，与二维剖面图、 中段图成果进行整合显示 模型可以按多个元素联合显示和切割，并保存平面点线区文件 ；切割过程可以根据勘探线和不同标高方向、纵剖面方向切割 可以自定义切割起止高程和切割步长，自动获得不同中段的多 个切割面。
培训提纲
地质统计学法资源量估算 流程与操作
地质统计学法估算流程与操作
原始样品数据统计分析 样品等长度组合 数据分析 实验变差函数计算及拟合 克里格品位估值 储量分级 报表生成
1 原始样品分析
原始样品数据分析
了解原始采样数据分布形态以及数据录入是否正确。
1 原始样品分析
原始样品样长分析 了解原始采样长度的平均值、中值等信息，确
8 块体品位模型应用
3.基于矿块模型计算采空区动用储量
定最佳采样长度。
• 平均值：1.86 m • 中值： 2m
• 组合样样长选用 2m
2 样品等长度组合
划分组合样 为避免在克里格插值中引入误差，对原始样品
按一定长度进行组合。
• 组合样样长：2m • 最小样长：0.1m • 划分方式：长度加权平均 • 划分范围：忽略无样品区
点击“确定”生成一个.cbs文件
变程。
次轴方向
• 变差函数模型类 型：球状模型
• 块金值：0.3
• 基台值：2.185
• 变程：31.8m
4 实验变差函数计算及拟合
变差函数分析结果
利用球状模型对变差函数进行理论模型的拟合获得的三个方向的 理论模型。根据三个方向的理论模型拟合，可以看出矿体在各方向上呈 现变异程度相同而变异连续性不同的几何异向性特征，即三方向基台值 相同、变程不同。
进行数据变换。
• 最小值 ：0 • 最大值：11.77 • 中值：0.438 • 平均值：1.39 • 标准差：1.88 • 峰度：18.68
3 数据分析
概率图分析 了解组合样数据分布形态、数据中是否存在特
异值。
• 数据为偏态分布
• 存在特异值
• 概率值97.5%对应的 品位值6.08
3 数据分析
椭球半径与变差函数变程相一致。
• 搜索椭球半径与变 差函数变程相一致
• 搜索椭球旋转角度 与变差函数各向异 性角度相一致
• 搜索方式多次搜索
5 克里格品位估值（示例）
交叉验证结果
利用组合样对所获得的理论模型进行交叉验证验证获得的散点图 验证基本统计量，验证结构显示误差平均值趋近于0，标准中误差趋 近于1，说明获得的理论模型是可靠的。