地质建模型培训步骤.doc
地质建模型培训步骤 GEMCOM国际软件公司 SURPAC中国办事处 目 录 第1章 地质数据库结构1 1.1 分析实例地质数据1 第2章 创建一个Surpac地质数据库5 2.1 创建数据库结构5 2.2 编辑转换表8 2.3 导入数据9 第3章 三维空间显示钻孔13 3.1 设置钻孔显示风格13 3.2 显示钻孔16 第4章 地质解译和剖面品位计算20 4.1 创建剖面20 4.2 地质解译(圈定矿体)24 4.3 计算线段品位31 第5章 形成实体模型33 5.1 建立实体模型33 第6章 数据提取38 6.1 提取钻孔平面图/剖面图38 6.2 提取样品表中数据,获得取样长度41 6.3 统计样长,确定组合样长43 第7章 组合46 7.1 样品组合原理47 7.2 根据勘探工程组合47 第8章 基本统计与分析54 8.1 基本统计54 8.2 处理特高品位59 第9章 创建块体模型62 9.1 创建块模型63 9.2 增加属性及背景值66 9.3 约束68 9.4 块体模型的显示70 第10章 块体模型赋值74 10.1 直接赋值法74 10.2 距离幂次反比法78 10.3 根据属性为模型着色83 10.4 第二次估值85 第11章 资源储量报告93 11.1 块模型的报告93 II Error No text of specified style in document. 第1章 地质数据库结构 1.1 分析实例地质数据 该实例为一套培训数据,作为初学者学习所用,此例为仅有钻孔数据的金矿床。 1. 打开Surpac,将培训文档中的地质建模设为当前工作目录。 2. 右键编辑开孔坐标.txt、测斜数据.txt、岩性数据.txt和化验数据.txt,这四个表已经按照SURPAC数据库的要求所整理,本章所做的所有工作都源于并且仅源于这四张表。现分析四个文件的内容 开孔坐标.txt内容 该地质数据与Surpac地质数据库的对应关系如下表示。 开孔坐标.txt 数据库中强制表Collar 注释 孔号 hole_id 钻孔的编号,不能重复 Y y 孔口的北坐标 X x 孔口的东坐标 Z z 孔口的标高 最大孔深 max_depth 终孔深度 孔迹线类型 hole_path 见说明 钻机类型 钻机类型(需手动添加) 钻机类型 说明 1. 前六项均为强制字段,软件自动创建。只需手工添加钻机类型这一选项字段。 2. 孔迹线类型curved或linear或vertical 钻孔为curved; 若是探槽或坑道应为 linear; 若工程是铅直向下的为vertical,倾角为-90,方位角为0 测斜数据.txt内容 该地质数据与Surpac地质数据库的对应关系如下表示。 测斜数据.txt 数据库中强制表Survey 注释 孔号 hole_id 钻孔的编号 测斜深度 depth 测斜深度 倾角 dip 该深度至下点的倾角 方位角 azimuth 该深度至下点的方位角 注意 1. 倾角向上为正,向下为负 岩性数据.txt 该地质数据与Surpac地质数据库的对应关系如下表示。 岩性数据.txt 数据库中选项表(需手工创建) 注释 孔号 hole_id 钻孔的编号 样号 sample_id 样号,可以为空 深度自 depth_from 深度至 depth_to 岩性代号 岩性(需手工添加) 岩性代号,可用中文 注意 1. 深度自和深度至不能超过最大孔深,如有重叠,会报错。 2. 共有四种岩性VOL,OVB,QV1,QV2。 化验数据.txt内容 该地质数据与Surpac地质数据库的对应关系如下表示。 化验数据.txt 数据库中选项表(需手工创建) 注释 孔号 hole_id 钻孔的编号 样号 sample_id 样号,可以为空 深度自 depth_from 深度至 depth_to Au au(需手工添加) 化验的au品位 说明 1. 化验的Au品位单位为g/t,因化验精度的限制,当品位打开/新建”。 2.输入建立的数据库名称一般情况下,应用项目名称来命名。 3.选择数据库类型信息一般是根据用户计算机安装的数据库类型来决定。 后台数据库为access2000,注意如果使用access数据库,在本机中必须安装Access数据库,可以通过MicroSoft产品Office中安装。如果不想安装其它数据库产品,Surpac建议使用Paradox数据库,使用这个数据库不需要其它产品的支持。 存储 – 是将各样品中心点的坐标,存储在数据库,这样会增加数据库的大小。 计算y x z测量和采样坐标 – 是当需要的时候,软件自动计算样品中心点的坐标,推荐使用该选项。 4. 定义要创建的选项表 如前所说,Surpac会自动创建3个强制表 collar, survey 和translation。 在该实例中,因数据库还缺少岩性表和化验表,因此创建这两个选项表。 5. 为collar(孔口) 表添加选项字段钻机类型 强制字段表中即为对应表中各强制字段的一些格式。 在该实例中,数据库中还缺少钻机类型这一字段,因此在这添加一项。其它设置如图。 6. 为岩性表添加选项字段岩性 7.为化验表添加选项字段au 注意类型选“实数”型。 8. 这样,我们定义的4张表和各自包含的字段,执行后,出现以下信息 说明Surpac已建立了一个以Access20000为后台的数据库,数据库的表和字段全部定义好了,并且数据库已经打开,但数据库里没有内容。 在左侧文件浏览器中,发现两个新文件 其中“培训数据库.ddb”为数据库定义文件,它是Surpac和数据库之间的桥梁,“培训数据库.mdb”为access2000数据库。至此,数据库已经建立起来了,但没有任何数据。打开 “培训数据库.mdb” 文件。用户在Access数据库中查看数据库格式和内容,确认建立的数据库格式正确。 以上步骤已录制为宏文件m1_创建数据库.tcl。 2.2 编辑转换表 Translation(转换)表是一个强制表,每个地质数据库中都有软件自动创建的该表,该表所包含的字段有 Field 字段名 table_name 表名 field_name 字段名 code 字段内容 num_equiv 转换后的字段内容 description 描述 在本例中,必须编辑转换表,否则化验表中一些非数字的记录将无法倒入到数据库中。下面是对转换表添加记录的步骤。若实际情况中所有化验结果均是数字,则无需编辑该表。 1. 编辑插入记录,为translation表插入记录。 2.选择translation表。 3.如下输入两条记录。 在该例中,translation转换表的作用 将化验表中au的化验值为“数据库中的表名称; 字段名称数据库表中的字段; 包含如果某字段不需导入,则不选。例如Survey中y,x,z自动计算,所以不需导入; 列该字段对应的txt文本文件中的列数 ; 这里要求反复分析text文本中的记录,一定要对应好,而且需要将所有的表中的字段都对应完成后方可“执行”。 5. 选择需要加载的txt文本文件 “进行样品重叠检验”在载入数据过程中,自动检查样品记录表中取样间隔是否有重叠 “最多错误数目”加载过程中,有可能text中数据格式与数据库中不一致,或其它错误引起导入错误,Surpac会自动报告错误,如果错误数目大于指定值(用户自定义),则停止加载,本例选择50个。 载入类型插入,更新,插入并更新,我们选择插入。有的情况,需要覆盖以前的数据,或追加新的数据,则需要用到更新,插入/更新功能。 执行后,出现一个文本文件,显示导入结果。 以上步骤已录制为宏文件m3_导入数据.tcl。 97 第3章 三维空间显示钻孔 一旦将地质数据库建立起来,我们就可以利用Surpac强大的图形显示系统,在三维空间显示地质数据,包括钻孔的轨迹线,品位值、岩性及代码、岩层走向等,总之,几乎所有的地质信息都可以以字符、图表、图案方式显示出来。打开菜单数据库显示,得到下面的子菜单 3.1 设置钻孔显示风格 显示钻孔时,我们必须做一些准备工作。如果希望沿着钻孔的方向,对不同的岩性着不同的颜色、不同的品位区间显示不同的风格,则需要设置钻孔显示风格, 任务 设置钻孔显示风格 1. 运行菜单数据库显示钻孔显示风格,得到如下界面 2. 在岩性表处,点击右键,获取字段代码,将不同的岩性设置为不同的颜色, 各岩性颜色设置如下, OVB periwinkle QV1 red QV2 red VOL blue 3. 在化验表处点击右键,获取最小-最大范围,以在不同的范围区间设置不同的颜色, 4. 区间的确定要参考对应项目的工业指标,比如边界品位、块段品位、矿区品位等指标。 各品位区间颜色设置如下 -990 periwinkle 01 cyan 13 pink 35 red 5999 plum 执行,则对于不同的岩性和品位,会显示不同的颜色。 3.2 显示钻孔 Surpac数据库中,对于钻孔的显示有多种效果,比如钻孔根据品位的不同以圆柱体的方式显示孔迹线,将孔号和终孔深度等孔口表中的信息显示到钻孔上,在孔迹线左侧显示岩性图案,标注任意数据库中有的数据等等。下面只介绍几种最常用的显示方式,其余可自行尝试应用。 1. 运行菜单数据库显示钻孔,在孔迹线风格里,填写如下 2. 在孔口风格里,选hole_id,以在孔口位置标注孔号, 3. 在地质图案中,如下填写, 4. 在标注栏中如下填写, 5. 定义查询约束为空,即对所有钻孔进行显示,点击执行, 6. 在图形工作区中即可以以上定义的方式显示所有钻孔,如图, 其余选项卡中的内容是, 标签,是指在轨迹线的一侧,显示化验品位值或岩性代码字符串等。本例中,在轨迹线的右侧显示au(金)的品位,左侧显示tau的品位,分别是左右对齐,偏移量为1。组合是指如果相邻的样本如果值相同,则可以合并为较长样长的单一值来显示。如果使用风格,则本例会将au品位字符串显示为在01之间为兰色,13之间为黄色,35之间为红色,5以上为绿色,如果不显示风格,则需给定颜色,所有字符串都同一种颜色。 图表,是指在钻孔的一侧,绘制出品位曲线图或直方图,来直观的查看品位分布情况,本例选择“化验表”,字段为au,图形类型为曲线(line),不使用风格,采用红颜色,长度/单位比值为11,超过50则不显示。 深度标记会在钻孔轨迹线的一侧,描述钻孔的长度标记。本例中,会沿钻孔轨迹线,30米间隔,用红色字符来记录钻孔长度。 视倾斜指示器,是用来描述岩层的走向和倾角的,一般在煤矿中,钻孔不仅要记录煤层代码,而且还要记录煤层在此处的走向和倾角,以便判断矿体的发展趋势。 以上步骤已录制为宏文件m4_显示钻孔.tcl。 第4章 地质解译和剖面品位计算 前面已介绍了地质数据库的建立和三维空间显示,下面介绍任何利用地质数据库和三维显示功能,来圈定矿体。本节内容是地质工程师所必须掌握的,是建立矿体实体模型和块体的基础,是地质建模的基础部分。 一般而言,钻孔会按照勘探线来布置,勘探线的间距和钻孔在同一勘探线的距离,称勘探网度,这些参数一般由地质勘察工程师确定。我们将沿着勘探线,建立一系列的剖面,在各剖面内确定矿体边界,各剖面线的矿体连成三角网,就生成矿体的实体。本节的内容就是如何定义剖面、圈矿体、报告各剖面矿体的品位,为实体模型和块体模型做准备。 4.1 创建剖面 在Surpac中,利用切剖面的工具,可在任意位置切剖面。但是,由于要进行剖面解译,因此应该在勘探线上创建剖面。 任务 创建勘探线剖面 1. 连接培训数据库.ddb,显示钻孔,或者直接点击m4_显示钻孔.tcl。双击地表.dtm,将地表也显示出来。 2. 显示二维坐标网。点击显示二维坐标网。如下填写。 3. 在图形工作区中,调整大小,可以看到本例中有7条勘探线,东西向,分别在y6750,6800,6850,6900,6950,7000,7050分布。 4. 创建6750剖面。点击剖面定义。如下填写,注意确认端点坐标。 5. 在实体剖面选项卡中,选上实体剖面,从而可将地表这个DTM进行剖切,得到的地形线存到剖面图层地形线中。 6. 执行后,在图形工作区中按住左键拖拽出一条线。松开左键。 7. 在随后出现的对话框中,将Y坐标都改为6750,从而沿Y6750 这条线创建剖面。 8. 在图形工作区中显示Y6750 剖面。 以上步骤已录制为宏文件m5_创建6750剖面.tcl。 4.2 地质解译(圈定矿体) 地质解译,也就是在勘探线剖面上圈定矿体,是地质工程师的一项重要工作。也是在Surpac 中建立实体模型的基础。地质解译工作,不能随意,应该参照不同的矿床的工业指标和解译方式。 在本例中,主要的解译原则是见矿点(大于边界品位1g/t)间直线连接;有限外推为工程间距的1/2;边部1/4相邻工程间距尖灭。 任务 解译6800剖面 1. 连接培训数据库,打开地表.dtm,显示钻孔(或直接双击m4_显示钻孔.tcl)。 2. 创建6750剖面,或直接双击m5_创建6750剖面.tcl。 3. 点击后一剖面,得到6800剖面, 图形工作区中显示 4. 因见矿点间直线连接,因此需要捕捉到见矿点。在工具栏中,点击捕捉到点。 5. 在图层窗口中点击,将各剖面解译线写在不同的图层,从而分别保存。 6. 通过显示的au品位,判断见矿点,将见矿点间采用画线工具圈定出矿体。 得到如下结果 注意按顺时针依次画线。 7. QV1矿体的右上应做有限外推,其中地表作为一种限制。 1 改变线号,点击状态栏中,编号设为2. 2 应用两点间取中点的工具,根据操作提示,从而得到两个中点。 3 连接两个中点延伸至地表。 4 再应用取中点工具,找到外推点。如图示。 5 点击编辑点插入,根据提示,在见矿点附近拖拽至外推点。 6 如上步骤,将QV2的外推点也找到,并调整矿体边界线。 8. QV1矿体的左下应做无限外推, 1 先找到中点,如上述方法, 2 报告两个中点的方位和距离,注意由于是向左下外推,因此应先选右上点,再沿左下方向,选另一点,报告。 报告的结果显示在消息窗口中, 3 编辑移动特定线段,点击左下角的点,在弹出窗口中,如下填写。 注意方位和倾角都等于第2)步中报告出的结果,而倾斜距离应该等于报告的距离的1/4。结果如下, 4 编辑点插入,将无限外推点插入到矿体边界线中。 5 如上步骤,将QV2的外推点找到。 9. 将2号线删除,得到结果如下, 10. 将矿体解译线保存, 11. 重复310步,得到所有剖面的解译线,如图。 12. 对于矿体的两端,楔形外推, 2. 新建图层“矿体解译线7050侧”。 3. 画2号线。 4. 找到中心点的方位和距离,以相邻剖面间距离的一半移动,得到如下结果 5. 重复以上步骤,得到矿体两端的界线。 在文件导航器中,得到矿体解译线6750侧.str,矿体解译线6750.str,矿体解译线6800.str,矿体解译线6850.str,矿体解译线6900.str,矿体解译线6950.str,矿体解译线7000.str,矿体解译线7050.str,矿体解译线7050侧.str等9个线文件。 4.3 计算线段品位 应用所画线段的品位这一功能可确定剖面上所画的线段的品位。 任务 确定所画线段的品位 1. 连接培训数据库。 2. 运行m4_显示钻孔.tcl。 3. 运行m5_创建6750剖面.tcl。 4. 打开矿体解译线6750.str。 5. 选择 剖面所画线段品位,如下输入,点击执行。 6. 执行后,根据提示,选择一个线段。可以在消息窗口中看到结果, 7. 按ESC退出该命令。 8. 点击查询点属性,选择刚才操作的那条线段上的点,可以在消息窗口中,看到其d1字段上即是该线段的au品位。 注意有两种计算线段内平均品位的方法,长度加权和面积加权。 长度加权利用样长为权,来计算闭合线的平均品位。针对上图,公式如下 平均品位A*d1B*d2C*d3/ABC 其中A、B、C为样品长度,d1,d2,d3为样的平均品位,grad为最终的平均品位。 面积加权在两个钻孔的中间,定义一条边界线(如图兰线部分),样品在水平方向、矿体界限内及边界线内的公共投影面积(蓝色部分),即为该段样品的权值。公式如下 S1 平均品位S1*d1S2*d2S3*d3/S1S2S3 S2 S3 第5章 形成实体模型 5.1 建立实体模型 标准化分割线段 实体建模先处理两端 加入辅助线 利用相连段间-结束端连接(注意选择顺序,不要封闭,先矿体的闭合段,后是尖灭到的线) 其他部分 检验实体 第6章 数据提取 在Surpac中,可提取地质数据库中的数据,形成线文件,用于绘图、统计等方面。 6.1 提取钻孔平面图/剖面图 钻孔数据存储于地质数据库中,但可从数据库中提取出需要的相关数据用于绘图和估值,下表中是几个命令对应的功能及目的, 菜单 功能 目的 用于打印的平面图 提取某一(或一系列)水平的钻孔,存为线文件,其中钻孔的各组成部分用固定的编号表示; 用于绘制或者查看平面图 用于打印的剖面图 提取某一(或一系列)剖面的钻孔,存为线文件,其中钻孔的各组成部分用固定的编号表示; 用于绘制或者查看剖面图 用于计算的平面图 提取某一(或一系列)水平的钻孔,存为线文件,其中化验品位作为描述字段; 用于估值 用于计算的剖面图 提取某一(或一系列)剖面的钻孔,存为线文件,其中化验品位作为描述字段; 用于估值 1. 设置好当前工作目录,打开培训数据库,点击数据库数据提取用于打印的剖面图。 2. 定义剖面方向,如图示,填写。 3. 选择化验表中要提取的品位,如图示填写。 4. 选择岩性字段,如图示填写。 5. 定义查询约束,即定义要提取那些钻孔,因提取所有,所以默认空,执行。 6. 生成7个钻孔剖面,打开钻孔剖面6800.str,如下图示。 7. 右键编辑,查看该线文件,各线号含义。 以上步骤已录制为宏文件m6_提取所有钻孔剖面.tcl。 6.2 提取样品表中数据,获得取样长度 任务 提取样品表中数据,获得取样长度 目的提取出取样长度以进行统计,确定组合样长度 1.数据提取 样品表数据,如下填写, 文件名处可自行定义。 样品坐标点可选底部、孔口、中间、顶部。意为,记录提取出的数据的点所要在的位置。 2.选择化验表。 3.选择要提取的样品表中的字段,可如下填写,执行。 4.选择所有样品的方式,执行。 5.对所有数据进行提取,即无约束,执行。 6.得到结果,点击文件导航器中生成的提取化验表中全部样品的数据.str ,点击右键,查看其结果。可以看到D5 即为所要的样长。另外,D1为样号,D2为depth_from, D3为depth_to, D4为au,D5为样长,D6为提取出的样长与原样长之比。 以上步骤已录制为宏文件m7_提取全部样品样长.tcl。 6.3 统计样长,确定组合样长 任务统计全部样品样长,以确定组合样长 1.分析 基础统计窗口, 2.选择上节中的结果线文件,由于要对样长进行统计,选择D5。 3.执行后,可选择多种方式,比如柱状图,频率分布曲线等方式查看样长的分布。 3. 点击 统计 报告,对样长进行报告。 4.报告结果如下所示。 可以看出,样长的均值是2.006米,绝大多数是2米,因此可将组合样长定为2米,这在之后的组合中会使用到。 第7章 组合 到目前为止,我们已经建立了数据库和进行了剖面解译,下一步,要将数据库中au的品位提取出来,保存在线文件内,以便于分析样品中au品位分布规律和块体插值。 组合的含义是,将空间不等长的样长和品位,量化到一些离散点上,如au品位,在该点附近,最有可能是这个值。组合样最终产生一些离散的点,除了三维坐标外,在它的描述字段中,存放该点最有可能的品位值。 组合的方式有多种,但是只有根据勘探工程和根据台阶高程这两种方式组合得到的线文件,能够在工程的方向上,产生均匀(等距离)的离散点,从而才可用于地质统计和作为块体模型估值的数据库进行估计插值。 Surpac提供了以下6种组合样的方法 Composite downhole根据勘探工程 Composite by grade根据品位约束 Composite by geology根据地质约束 Composite by elevation根据台阶高程 Composite From End Of Hole自钻孔末端 Multiple elements多种元素 在Surpac6.0中,组合命令有如下, 7.1 样品组合原理 如何组合样品一般采用长度加权法如下图 右侧是每个样的长度和品位值,在左侧,按照2米间距,可选择在中间位置产生一个点,来描述该点品位组合后的值 在此图中01米的样,品位为1.00;12.5米的样,品位为2.00;2.55米的样,品位为3.00; 按照组合样定义,在02米的中心位置,该点的品位为 grade1*1.001*2.00/111.50 同理,在24米的中心位置,该点的品位为 grade0.5*2.001.5*3.00/112.75 在1213米处,出现负值,这在实际情况是不存在的,我们有不同的处理方法,本例只是简单忽略未当零处理 grade1*7/17.00 在1415米处,此处只有1米长,占组合样长2米的一半,次此处依然包含 grade1*7/17.00 注意出现负值或样长不够,Surpac用另外一种线号,来区分正常的样品。对待负值或样长不够的情况,Surpac提供了可选的处理方案,在后面会介绍到。 7.2 根据勘探工程组合 任务 根据勘探工程组合 1. 选择组合根据勘探工程, 2. 如下填写,点击执行。 3. 因为无约束,所以不需填写,点击执行。 4.重置图形工作区。 5.打开根据勘探工程组合_无区带.str到图形工作区中。 6.点击显示隐藏线在图层中。 7.显示点标记。 8. 结果如图示, 9.在文件导航器中,在该线文件处点击右键编辑,可查看各描述字段的含义,其中注意到D1字段即为组合后的品位。 任务 在一地质带内根据勘探工程组合 应用此步骤,可实现根据一定的取样长度对矿体内部的样品进行组合。 1.确定组合样长,如7.4节中所述,统计得到的平均取样长度为2.006m,因此组合样长可设为2m。 2.由于要对矿体内部组合,所以组合应该约束在一个地质带内。从第2步到第9步,是为了得到矿体内部这样一个地质带。 新建一个表用于存放钻孔与矿体相交的结果。数据库数据库管理创建数据表,如下图填写, 3.定义一个选项字段,叫地质带名称, 4. 打开矿体.dtm, 5. 数据库分析钻孔与3dm相交,无约束,默认为空,执行, 6. 如下填写, 7. 弹出一个报告, 8. 同样,对于2号体,应为qv2矿体,如下填写, 9.结果如下, 10.对qv1矿体内部的样品进行组合,组合根据勘探工程,如下填写, 11.定义地质带,由于我们进行了钻孔与矿体相交工作,得到了一个名为钻孔与矿体相交的间隔型表。如下填写, 上述操作能够得到qv1矿体内的组合样品。 12.同样,对qv2矿体内部的样品进行组合,组合根据勘探工程,如下填写, 13.定义地质带,如下填写, 上述操作能够得到qv2矿体内的组合样品。 第8章 基本统计与分析 产生组合样品点后,我们要统计样品的值与空间分布情况,以了解矿体分布规律,为确定特高值和块体模型估值作准备。 8.1 基本统计 在上章中得到了组合后的样品点,我们要对样品品位的离散值进行一些基本统计,以更好的掌握这些原始数据,并确定特高值,进行处理。 1.对qv1矿体内组合样的Au品位进行基本统计,运行菜单数据库分析基本统计窗口,弹出对话框 2.可以看到如下结果, 这是基本统计窗口,其中 X轴表示au的品位 Y轴表示频次 直方图表示该品位出现的频次 曲线表示积分频率曲线 以上是自然分布状态,如果需要对数分布的曲线,按动,输入自然对数, 注意,如果遇到下图双峰分布曲线,则不用按照以下常规的方式处理,需要先划分域,以将两部分分开,具体处理方式在本书中不做叙述。 3.统计报告,如下填写, 得到如下结果, 样本数参与统计的样本数目 91 最小、最大值样本的最小、最大值 均值样品的平均值,本例为15.107 方差值描述样品围绕均值变化的大小。公式为 标准离差方差值的开方, 变量系数 偏度,描述样品对称分布的情况,如果完全对称分布,则值为零 峰度,描述样品分布峰值的情况 其它属于高级统计部分,可参照相关的帮助文档 4.对qv2矿体内组合样的Au品位进行基本统计,运行菜单数据库分析基本统计窗口,弹出对话框 2.可以看到如下结果, 这是基本统计窗口,其中 X轴表示au的品位 Y轴表示频次 直方图表示该品位出现的频次 曲线表示积分频率曲线 3.统计报告,如下填写, 得到如下结果, 8.2 处理特高品位 特高品位指的是极个别的与平均值偏离很多的值。通过前面的基本统计,发现以下问题 为了消除特高品位的影响,必须将其处理,处理方式有多种,比如使特高值等于均值,或赋一个其它的数值,或者直接删掉这些值。我们在此选择第一种方式。 首先我们确定特高值的下限,这也有多种方式,比如通过统计、类比法等等。一般我们选取分布大于97.5的样品值,即截取2.5的高品位,或均值的68倍。 95 CI 均值 1.96 * 标准差 处理特高品位时通常是分矿体来进行。 在本例中,qv1矿体即为15.107 1.96 * 65.142 142.785,可取为143; qv2矿体即为 5.375 1.96 * 3.731 12.688,可取为13 1.处理特高值qv1矿体内的特高品位 打开根据勘探工程组合_qv1矿体内.str. 2.编辑图层运算,如下填写, 其中 iifd1143,143,d1,其含义为 如果d1中的值大于143,则变成143,如果d1小于143,d1保持原值不变。iif是个判断语句,判断第一部分成立,选择第2部分,否则选第3部分。 3.保存。 4.处理qv2矿体内的特高品位 打开根据勘探工程组合_qv2矿体内.str. 5.编辑图层运算,如下填写, 其中 iifd113,13,d1,其含义为 如果d1中的值大于13,则变成13,如果d1小于13,d1保持原值不变。iif是个判断语句,判断第一部分成立,选择第2部分,否则选第3部分。 6.保存。 得到的此结果文件经过了特高值处理,我们可以对它再进行基本统计。这个文件即是之后要估值的数据源。 创建块体模型 第9章 创建块体模型 建立块模型和品位模型流程 9.1 创建块模型 1. 在图形窗口中,调入文件矿体.dtm 得到如下显示图像 2. 选择查询报告层范围得出矿体范围及走向。 在信息窗中可以得出X、Y、Z的范围 3. 选择块体模型块体模型新建/打开。出现选择模型的表,输入模型的名字,长度不超过40个字符。 4. 输入如下所示的信息,点击执行。 注意如果输入模型的名字不存在,就创造一个新的模型。如果是已有的模型,则直接调用。当创建了约束文件,载入此文件,则会在约束文件的基础上,创建一个新的块体模型。这是非常有用的一点,一方面可以快速显示约束范围的块模型,另一方面,可以精简模型属性,提高运算速度。 5. 出现下面的对话框,点击执行。 6. 输入如下所示的信息,点击执行。 参数说明 l 描述 模型的描述是可选的,但是它的纪录对于描述模型的目的是很有用的。 l 定义块模型使用 通常情况下使用,如果使用,则原点坐标即为最小坐标。 l 定义块模型使用范围 在这里可以选择三种方法中的一种即用“最小/最大坐标”或“原点坐标/范围”,二者是相互运算的,也可以选择“自线文件得到范围”。这里选择坐标最小值Y6720、 X2360、Z90;最大值Y7080、X2530、Z160。通常最小值和最大值一般取整。 l 用户块大小 以前讨论过,根据矿体的形态和工程控制网度,选择用户块尺寸10X10X5。通常在XY平面块尺寸一般为勘探线间距的1/31/5,在Z方向上块尺寸一般为组合样长度的23倍或是台阶高度的整数分之一。但是块在块模型范围内必须是整数个块。 l 旋转 当选择“旋转”时,要求对模型的旋转方位、倾角和倾伏角进行定义。这里暂都设定为零“0”,表明不旋转。 l 次级分块 有三个选项标准的(Standard)、可变的(Variable)和无(None)。这里选择标准的(Standard)。 l 最小块尺寸 软件中已经设置了不同数级(2n)进行缩分,用户只需根据要求选择即可。这里选择5X5X2.5。 l 保持审计跟踪 将来在“块体模型块体模型显示模型审核记录”中记录模型的操作、编辑和修改时间。 7. 得到下面的窗口,点击执行。 显示模型确认表,可以检查模型的大小,旋转等参数,如果你需要更改任意参数,点击取消按钮;确认参数,则点击执行按钮。 8. 选择块体模型块体模型保存,保存刚才新建的块体模型。 注意在块模型中操作,所有的数据存在内存中,意味着你总是工作在模型的复印件中,任何改变都没有存储, 直到你退出模型的时候,在退出之前一定要存盘,因而,随时存盘,就能避免意外发生(如突然停电),在填充操作后(手册后面会讲到),所有的值发生了变化,存盘将覆盖了以前的模型。 至此块体模型已经创建完成,状态栏里会出现按钮表示块体模型已经创造成功了,并且已经完成连接。 模型建好之后,就不能改变范围,旋转,块大小,每条边的最大限度,要改变模型的参数就只好重新创造一个新的模型了。 9. 选择块体模型块体模型关闭,关闭块体模型。 当完成相应的工作后,需要退出块体模型,选择“关闭”功能,弹出关闭对话框,用户可以根据操作的情况选择退出方式“保存并退出”、“保存为新模型并退出”或“不保存退出”。 如果你想要看这个任务的所有操作过程,请运行 _m01_创建块体模型.tcl 。 注意无论任何时候当宏暂停时,在消息窗口出现“单击图形工作区继续”提示信息时,你需要在单击图形工作区。当然你也需要单击弹出窗口的执行按钮。 9.2 增加属性及背景值 前面已经描述过,块体模型实际上是一个空间数据库,它的目的就是用来存储相关地质信息(包括岩石类型、品位分布、比重等等)。而这些属性是通过块体进行反映出来,创建块模型后,开始增加属性,属性是模拟模型空间的道具,这些属性可以是数字,字符代表顺序,间隔,比率等,属性也可以是通过别的属性计算得出。 1. 打开文件块体模型.mdl。 2. 选择块体模型属性新建。 3. 输入如下所示的信息,点击执行。 参数说明 对于新创建的块体模型,需要“新建”属性并且必须命名。 l 属性名字 属性名字的长度不超过30个字符,属性名字中允许空格,不过我们并不推荐这样,这样会使块的数学计算的功能复杂化。在左边数字或下面灰白处点击右键可以“插入、增加和删除”属性栏。 l 属性类型 属性类型可以是字符(Character)、实数(Real)、整数