视频配套操作步骤地质统计学

1、SURPAC 软件地质统计学（操作步骤）GEMCOM 国际软件公司SURPAC 中国办事处目录第 1章地质域11.1 查看地质域1样品组合32.1 在 SURPAC 软件中提取地质域中的数据3基础统计73.1 显示柱状图7特异值104.1 特异值104.2 确定并处理特异值104.2.1 特异值下限的确定方法10第 2章第 3章第 4章4.2.2特异值的处理方法104.2.3用 SURPAC 软件处理特异特异值114.2.4验证处理结果11第 5章方差图155.1 全方向向变异函数155.2 查看主方差图的倾斜平面175.3 计算主方向方差图195.3.1 查看已有的主方向方差图结果195.3

2、.2计算主方向方差图225.4 计算次主方向变差图265.4.1 查看已有的次主方向方差图结果265.4.2计算次主方向方差图285.5 异向性椭球参数325.5.1 查看已有的椭球体325.5.2异向性椭球参数33第 6章交叉验证376.16.1 交叉验证37创建 3 维块模型397.1 创建一个三维块模型39通克立格法428.1 用异向性椭球在 3 维块模型上进行普通克立格估值428.2 各向同性椭球在 3 维块模型上进行普通克立格估值45第 7章第 8章I第9 章部分百分比计算489.1 计算块模型中的块体在 3DM 实体中的百分比48第10章 模型检验4910.1 勘探线数据与模型比较

3、4910.1.1 实体内块模型百分比与矿体的实体剖面比较4910.1.2 块体模型的克立格估值与组合样品的对比5110.2 品位吨位曲线5410.2.1 创建品位吨位曲线5510.3 块模型值的基础统计5810.3.1 显示块模型值的基本统计5910.4 趋势分析6010.4.1 导出块质心6010.4.2 趋势分析61II第 1 章地质域1.1 查看地质域1.打开 all_composites2.str。2.选择显示>全部隐藏。3.选择显示>点>标记。4. 按下面的窗体输入信息，点击执行。5. 选择显示> 3D 网格.6. 按下面的窗体输入信息，点击执行。7. 在图形

4、工作区中进行点击和拖动，将看到如下的旋转的数据。1线串文件中的点所有钻孔以 2 米为组合样长的样品品位。D1 描述字段包含金的组合品位，根据品位的高低为分不同的线串，见下表。8.再打开 ore_solid1.dtm 文件。实体模型表征了一个单独的地质域，只有落在域内部的组合样品才能用于估计位于地质域内的块体模型。2第 2 章样品组合2.1 在 SURPAC 软件中提取地质域中的数据域ore_solid1.dtm表现的是QV1地质带，下面将完成中仅在QV1域中提取组合样品的操作：1.数据库>打开/新建2.数据库>组合>根据勘探工程点击执行定义地质带：在下一个窗体，输入“geol

5、ogy”表中的字段“rock”被定义为地质带，“geology”表是区间型表，“rock”字段中的“QV1”标明了每个钻孔穿过矿体“ore_solid1.dtm”的区间。33.文件>打开>线/DTM4. 显示> 全部隐藏5. 显示> 点> 标记46. 显示> 点> 属性7. 数据库>显示>钻孔8. 查看> 观察平面> 创建图形剖面59.查看>观察平面>后一个平面执行 7 次上述可由运行宏文件_01_创建钻孔组合.tcl来实现。6第 3 章基础统计3.1显示柱状图块体模型>地质统计>基础统计1.文件>

6、 从线文件导入数据2.点击执行。3.74. 统计> 报告5.点击执行。显示输出文件“raw_gold.not”内容。89第 4 章特异值4.1 特异值“特异值”是一个统计名词，是指在数据母体中那些比全部数值的平均值或中位数高得多的数值。例如，下面的数据母体，236 将被认为是一个特异值。1 3 5 5 8 8 8 2 3 6特异值在变异函数计算时，能产生问题。另外，在估值中有特异值的存在，将产生不切合实际的结果。用于减少特异值的影响的办法是对特异值用一个边界品位替代，在上面的例子中，236 的数值将被“截取”， 用数值“9”替代。1 3 5 5 8 8 8 9另外一种办法就是去掉特异值。

7、4.2 确定并处理特异值查看一下 gold_comp2.str 的柱状图和输出的报告，将会发现，多数的品位数据都集中在 0 至 10g/t之间，同时也发现几个比 10g/t 大的特异值。4.2.1特异值下限的确定方法1.传统方法:用平均品位的 68 倍做为特异值的下限2.数理统计法:对组合的样品进行基本统计分析，按 95%的置信区间确定特异值的下限.即:特异值下限=均值+1.96× 标准差3.本实例的特异值下限=3.828+1.96×6.831=17.2174.2.2特异值的处理方法1.传统方法:用平均品位替代特异值102.截取法:将高于特异值的下限的部分截取掉,相当于用特

8、异值的下限替代特异值3.本实例的特异值下限为了避免特异值对估值的影响，采取了“截取”的办法，即将超过 17g/t 的品位值替换为 17g/t。4.2.3用 SURPAC 软件处理特异特异值1. 线文件工具>输入下面信息：线串运算D1 描述字段将表达式的结果：iif(d1>17，17，d1)可表述为：如果D1的值大于17， 则D1的值等于17， 否则等于D1的原来值。2.点击执行。为检验线串运算的结果，你将在基础统计窗口中对数据进行分析，4.2.4验证处理结果1.块体模型> 地质统计> 基础统计112. 文件>从线文件导入数据3. 点击执行。12一个柱状图和一条描述

9、累计频数曲线。注意这时数据中的最大值是 17。4. 统计>报告5.点击执行。显示输出文件“gold_cut17.not”的内容。13上述操作可由宏文件03_截取特异值.tcl.实现。14第 5 章方差图5.1 全方向向变异函数1.块体模型>地质统计> 方差模型2.变异函数>新建线文件方差图3.点击高级选项卡154.点击执行。5.变异函数>模型。把滞后距修改为7.6165.2 查看主方差图的倾斜平面1.文件>打开>线/DTM2.显示>为 DTM 着色173.文件>打开>线/DTM4. 显示> 线> 带颜色充填18上述操作可由

10、运行宏文件 05_data_geometry.tcl 来实现。5.3 计算主方向方差图5.3.1查看已有的主方向方差图结果1. 将display_major_axis1.str文件拖到图形工作区。2. 将primary_variogram_map1.dtm文件拖到图形工作区。193.点击按钮。4.显示>为DTM着色。5. 将strike_dip_plane1.str文件拖到图形工作区。6. 显示> 线> 带颜色充填207.显示> 3D坐标网主方向方差图显示在图形工作区中215.3.2计算主方向方差图1.块体模型> 地质统计>方差模型2.方差图> 新建方

11、差图3.点击高级选项卡224.点击执行。5.移动滞后距滑块，并点击选项卡6. 变异函数> 模型237. 点击并拖动基台值和变程的标记点，进行拟合8. 方差图> 选择具有最大的连续性的方向249. 方差图>保存为DTM10. 文件> 试验方差图11. 文件> 方差图模型12. 差图窗口也可以平铺，点击右键，弹出窗口，选TileWindows25上述操作可由运行宏文件_06_primary_variogram_map.tcl.tcl来实现5.4 计算次主方向变差图5.4.1查看已有的次主方向方差图结果1.选择主程序窗口。2. 将display_semi_major_a

12、xis1.str文件拖到图形工作区。3. 将secondary_variogram_map1.dtm文件拖到图形工作区。264.显示>为DTM着色。显示> 3D 坐标网5.27主方向和次主方向方差图显示在图形工作区中5.4.2计算次主方向方差图1. 选择方差图窗口2. 方差图> 次变差图283.移动滞后距滑块，并选第10个选项卡变异函数> 模型4.295.点击并拖动基台值和变程的标记点6.方差图>保存为最有最大的连续性方向7.方差图>保存为DTM308.文件> 保存> 试验方差图9.文件> 保存> 方差图模型10. 方差图窗口也可以平

13、铺，点击右键，弹出窗口，选 TileWindows。31上述操作可由运行宏文件07_secondary_variogram_map.tcl来实现。5.5 异向性椭球参数5.5.1查看已有的椭球体1.将 ellipsoid2.str 文件拖到图形工作区中325.5.2异向性椭球参数1.方差图>沿轴提取试验方差图2.文件>打开>方差模型33调节滞后距3.4.选择 semi_major 选项卡，调整滞后距为 32，点击并拖动基台值和变程的标记点。345.选择 minor 选项卡，调整滞后距为 24，点击并拖动基台值和变程的标记点。上述操作可由运行宏文件08_anisotropy_e

14、llipsoid.tcl.6.方差图> 创建异向性椭球报告。357.方差图>椭球体观察仪。8.选择文件> 关闭36第 6 章交叉验证6.1 交叉验证1.重置图形工作区。2.块体模型> 地质统计方差验证3.输入下面信息，点击执行。3738第 7 章创建 3 维块模型7.1 创建一个三维块模型1. 删除 geostats_example.mdl 和 QV1.CON 两个文件。2. 块体模型>新建/打开3. 块体模型> 显示394. 块体模型>约束>新建图形约束405. 块体模型> 属性> 新建上述操作可由宏文件_09_create_blo

15、ck_m.tcl来实现41第 8 章通克立格法8.1 用异向性椭球在 3 维块模型上进行普通克立格估值1 块体模型>估值>普通克立格法4243块体模型显示如下：448.2 各向同性椭球在 3 维块模型上进行普通克立格估值1. 块体模型>估值>普通克立格法4546上述操作可用宏文件_10_ordinary_kriging.tcl来实现47第 9 章部分百分比计算9.1 计算块模型中的块体在 3DM 实体中的百分比重置图形工作区。1.2.打开 geostats_example.mdl，并进行显示.将约束文件 qv1.con 拖到图形工作区中估值> 部分百分比3.输入下

16、面信息，点击执行。4.块体模型> 保存5.6.块体模型>关闭上述操作可由运行宏文件_11_partial_percentage.tcl 来实现。48第 10 章 模型检验10.1勘探线数据与模型比较10.1.1 实体内块模型百分比与矿体的实体剖面比较1.将 ore_solid.dtm 文件拖入到图形工作区2.将 geostats_example_final.mdl 文件拖入到图形工作区3.将 qv1.con 文件拖入到图形工作区块体模型块体模型显示4.块体模型显示根据属性为模型着色5.496.查看观察平面创建图形剖面507.显示点属性10.1.2 块体模型的克立格估值与组合样品的对

17、比1.重置图形工作区，关闭块模型2. 将 ore_solid1.dtm 拖入图形工作区3. 将 surpac.ddb 拖入图形工作区4.数据库显示515.将 section_example.mdl 文件拖入图形工作区中6.将 文件拖入图形工作区中块体模型块体模型显示7.块体模型显示根据属性为模型着色8.529.查看观察平面创建图形剖面53显示点属性10.10.2品位吨位曲线另一个方法是报告品位、吨位，并构造一条品位吨位曲线。5410.2.1 创建品位吨位曲线1.将 geostats_example_final.mdl 文件拖入图形工作区中2.将 文件拖入图形工作区中块体模型块体模型显示3.块体

18、模型显示根据属性为模型着色4.555.块体模型块体模型报告565710.3块模型值的基础统计块模型值的基础统计也是检验块模型的法。58Anisotr opic Kriged Gold体积吨位Anisotropic Kriged GoldAnisotropic Kriged Gold0.0 -> 1.000001.0 -> 2.02371744743481.82861450.92.0 -> 3.0129266425853282.466366763.693.0 -> 4.059767811953553.454119692.534.0 -> 5.02523075046134.432233498.265.0 -> 6.0900041800085.4971777.366.0 -> 7.026998539966.44347821.57.0 -> 8.0.44196105.38.0 -&