DGSS-REInfo地质统计学资源量估算流程

1、培培 训训 提提 纲纲 地质统计学法估算流程与操作地质统计学法估算流程与操作 原始样品数据分析原始样品数据分析 了解原始采样了解原始采样数据分布形态数据分布形态以及以及数据录入是否正确数据录入是否正确。1 1 原始样品分析原始样品分析原始样品样长分析原始样品样长分析 了解原始采样长度的了解原始采样长度的平均值平均值、中值中值等信息，等信息，确确定最佳采样长度定最佳采样长度。 平均值：平均值：1.86 m 中值：中值： 2m 组合样样长选用组合样样长选用2m1 1 原始样品分析原始样品分析划分组合样划分组合样 为为避免在克里格插值中引入误差避免在克里格插值中引入误差，对原始样品，对原始样品按一定

2、长度进行组合。按一定长度进行组合。 组合样样长：组合样样长：2m 最小样长：最小样长：0.1m 划分方式：划分方式：长度加权平均长度加权平均 划分范围：划分范围：忽略无样品区忽略无样品区 点击点击“确定确定”生成一个生成一个.cbs文件文件2 2 样品等长度组合样品等长度组合划分方法划分范围组合样长最小样长品位阈值品位特高值同时划分加密样品导入矿体边界约束样品2 2 样品等长度组合样品等长度组合数据导入：加载离散数据点数据导入：加载离散数据点 将将.cbs文件文件导入系统；选择进行储量估算的矿导入系统；选择进行储量估算的矿体文件。体文件。3 3 组合样品数据分析组合样品数据分析直方图分析直方图

3、分析 了解组合样了解组合样数据分布形态、是否服从正态分布、数据分布形态、是否服从正态分布、进行进行数据变换数据变换。 最小值最小值 ：0 最大值：最大值：11.77 中值：中值：0.438 平均值：平均值：1.39 标准差：标准差：1.88 峰度：峰度：18.683 3 数据分析数据分析概率图分析概率图分析 了解组合样了解组合样数据分布形态、数据分布形态、数据中是否存在数据中是否存在特特异值异值。 数据为数据为偏态分布偏态分布 存在特异值存在特异值 概率值概率值97.5%对应对应的品位值的品位值6.083 3 数据分析数据分析位置图分析位置图分析 了解组合样数据了解组合样数据位置分布是否有误、

4、位置分布是否有误、数据数据是否是否存在趋势存在趋势。 数据分布正常数据分布正常 数据不存在趋势数据不存在趋势3 3 数据分析数据分析数据处理 正态变换正态变换 对数变换对数变换 指示变换指示变换变差函数模型变差函数类型4 4 实验变差函数计算及拟合实验变差函数计算及拟合实验变差函数计算变差函数参数变差函数参数变差函数图形变差函数图形各向异性角度各向异性角度变差函数表面图变差函数表面图4 4 实验变差函数计算及拟合实验变差函数计算及拟合变差函数拟合自动拟合变差函数自动拟合变差函数变差函数图形变差函数图形理论变差函数参数理论变差函数参数各向异性结构套合各向异性结构套合4 4 实验变差函数计算及拟合

5、实验变差函数计算及拟合数据变换及模型选择数据变换及模型选择 根据数据分析结果根据数据分析结果进行进行数据变换数据变换，特异值特异值处理处理，选择；通过计算，选择；通过计算全方向变差函数全方向变差函数求得最求得最佳步长，通过计算佳步长，通过计算短距短距离变差函数离变差函数求得块金值，求得块金值，通过通过定向变差函数定向变差函数计算计算求得各向异性角度求得各向异性角度4 4 实验变差函数计算及拟合实验变差函数计算及拟合计算变差函数（计算变差函数（全方向变差函数全方向变差函数） 步长数目步长数目*步长大小步长大小=区域长度的一半区域长度的一半 方位角，倾角均为方位角，倾角均为0度度。 方位角容差，倾

6、角方位角容差，倾角容差均为容差均为90度度。4 4 实验变差函数计算及拟合实验变差函数计算及拟合步长总间距步长总间距计算变差函数（计算变差函数（全方向变差函数全方向变差函数） 点击点击“计算变差函数计算变差函数”按钮，观察图形，寻找按钮，观察图形，寻找最佳步长大小。最佳步长大小。 最佳步长为最佳步长为20m4 4 实验变差函数计算及拟合实验变差函数计算及拟合计算变差函数（计算变差函数（短距离变差函数短距离变差函数） 点击上一步选择点击上一步选择“短距离变差函数短距离变差函数”，双击列表，双击列表控件设置步长大小，一般情况下步长大小为控件设置步长大小，一般情况下步长大小为组合样样组合样样长。长。

7、 块金值为块金值为0.34 4 实验变差函数计算及拟合实验变差函数计算及拟合计算变差函数（计算变差函数（计算三个相互垂直的变差函数计算三个相互垂直的变差函数） 选中选中“计算第二、第三方向计算第二、第三方向”输入方位角和倾角，点击输入方位角和倾角，点击“输入输入第二、第三方向参数第二、第三方向参数”按钮，计算三个相互垂直方向的变差函按钮，计算三个相互垂直方向的变差函数。数。 三个相互垂直方三个相互垂直方向的变差函数向的变差函数4 4 实验变差函数计算及拟合实验变差函数计算及拟合拟合变差函数拟合变差函数 通过滑动条设置变差函数的通过滑动条设置变差函数的块金值块金值、基台值基台值和和变程变程。主轴

8、方向主轴方向 变差函数模型类变差函数模型类型：型：球状模型球状模型 块金值：块金值：0.3 基台值：基台值：2.185 变程：变程：97.8m4 4 实验变差函数计算及拟合实验变差函数计算及拟合拟合变差函数拟合变差函数 通过滑动条设置变差函数的通过滑动条设置变差函数的块金值块金值、基台值基台值和和变程变程。半轴方向半轴方向 变差函数模型类变差函数模型类型：型：球状模型球状模型 块金值：块金值：0.3 基台值：基台值：2.185 变程：变程：80.6m4 4 实验变差函数计算及拟合实验变差函数计算及拟合拟合变差函数拟合变差函数 通过滑动条设置变差函数的通过滑动条设置变差函数的块金值块金值、基台值

9、基台值和和变程变程。次轴方向次轴方向 变差函数模型类变差函数模型类型：型：球状模型球状模型 块金值：块金值：0.3 基台值：基台值：2.185 变程：变程：31.8m4 4 实验变差函数计算及拟合实验变差函数计算及拟合4 4 实验变差函数计算及拟合实验变差函数计算及拟合克里格估值类型简单克里格简单克里格普通克里格普通克里格泛克里格泛克里格指示克里格指示克里格距离密次方反比距离密次方反比5 5 克里格品位估值克里格品位估值搜索椭球设置搜索椭球半径搜索椭球半径旋转角度旋转角度多次搜索多次搜索搜索类型搜索类型块段离散化个数5 5 克里格品位估值克里格品位估值交叉验证散点图散点图直方图直方图误差平均值

10、误差平均值标准中误差标准中误差统计量比较统计量比较5 5 克里格品位估值克里格品位估值数据变换及模型选择数据变换及模型选择 选择合适的克里格估值类型，并完成变差函数选择合适的克里格估值类型，并完成变差函数的导入。的导入。 克里格类型：克里格类型：普通普通克里格估值算法克里格估值算法 导入变差函数文件导入变差函数文件5 5 克里格品位估值（示例）克里格品位估值（示例）搜索邻域设置搜索邻域设置 通过一个通过一个搜索椭球搜索椭球定义搜索邻域，一般情况下定义搜索邻域，一般情况下椭球半径椭球半径与变差函数变程相一致与变差函数变程相一致。 搜索椭球半径搜索椭球半径与变与变差函数变程相一致差函数变程相一致

11、搜索椭球旋转角度搜索椭球旋转角度与变差函数各向异与变差函数各向异性角度相一致性角度相一致 搜索方式搜索方式多次搜索多次搜索5 5 克里格品位估值（示例）克里格品位估值（示例）交叉验证结果交叉验证结果 误差平均值：误差平均值：0.0068 标准中误差：标准中误差：1.062 小于中误差样品比例：小于中误差样品比例：89.64%5 5 克里格品位估值（示例）克里格品位估值（示例）克里格估值结果包括：块体模型克里格估值结果包括：块体模型中心点坐标中心点坐标、块体模型、块体模型样样品品位值品品位值、品位估计方差品位估计方差，结果保存到数据库中。，结果保存到数据库中。 数据库中包含：数据库中包含： 块体

12、模型参数表块体模型参数表 变差函数参数表变差函数参数表 空块模型表空块模型表 块体模型品位表块体模型品位表5 5 克里格品位估值（示例）克里格品位估值（示例）三维图形显示，用户可以三维图形显示，用户可以查询任意块体模型属性信息查询任意块体模型属性信息。5 5 克里格品位估值（示例）克里格品位估值（示例）用户用户输入矿石比重输入矿石比重，更新数据库更新数据库，重新计算储量重新计算储量。6 6 资源量计算资源量计算两种方式：一种根据估值过程中用到的两种方式：一种根据估值过程中用到的样品个数样品个数、工程个工程个数数以及以及搜索次数搜索次数进行区分，用户可修改分级名称以及区分进行区分，用户可修改分级

13、名称以及区分范围；一种根据克里格估值方差进行分级。范围；一种根据克里格估值方差进行分级。结果保存到数据库中结果保存到数据库中6 6 资源量分级资源量分级对数据库中矿体储量进行统计，对数据库中矿体储量进行统计，生成统计报表生成统计报表。7 7 报表生成报表生成8 8 块体品位模型应用块体品位模型应用 8 8 块体品位模型应用块体品位模型应用矿块模型可以与原始勘探工程数据、矿块模型可以与原始勘探工程数据、矿体约束信息组合显示矿体约束信息组合显示 8 8 块体品位模型应用块体品位模型应用可以根据某个矿体内矿块内某种元素的品位大小区间，将矿块按不同的颜可以根据某个矿体内矿块内某种元素的品位大小区间，将

14、矿块按不同的颜色进行显示色进行显示 矿块矿块( (单元素单元素) )显示显示8 8 块体品位模型应用块体品位模型应用可以根据某个矿体内矿块内某种元素的品位大小区间，将矿块按不同的颜色可以根据某个矿体内矿块内某种元素的品位大小区间，将矿块按不同的颜色进行显示进行显示 多元素集成显示多元素集成显示显示部分品位区间的矿块显示部分品位区间的矿块8 8 块体品位模型应用块体品位模型应用l在矿块模型显示的基础上，可以对该模型进行不同方向的切割在矿块模型显示的基础上，可以对该模型进行不同方向的切割，获得切割面上的相关信息，获得切割面上的相关信息l切割方向可以为：勘探线方向、纵剖面方向、中段方向切割方向可以为

15、：勘探线方向、纵剖面方向、中段方向l相应块体模型剖切结果可以保存为矢量数据，与二维剖面图、相应块体模型剖切结果可以保存为矢量数据，与二维剖面图、中段图成果进行整合显示中段图成果进行整合显示l模型可以按多个元素联合显示和切割，并保存平面点线区文件模型可以按多个元素联合显示和切割，并保存平面点线区文件；切割过程可以根据勘探线和不同标高方向、纵剖面方向切割；切割过程可以根据勘探线和不同标高方向、纵剖面方向切割l可以自定义切割起止高程和切割步长，自动获得不同中段的多可以自定义切割起止高程和切割步长，自动获得不同中段的多个切割面。个切割面。 8 8 块体品位模型应用块体品位模型应用 l根据实际测量数据，在三维视图中显示采空区边界范围根据实际测量数据，在三维视图中显示采空区边界范围l根据采空区边界范围和采空区所属采掘面间关系自动生成采根据采空区边界范围和采空区所属采掘面间关系自动生成采空区三维模型空区三维模型l基于采空区三维边界，手工交互生成采空区三维模型功能基于采空区三维边界，手工交互生成采空