基于MAPGIS6.7化探数据处理方法及异常图绘制.docVIP

基于MAPGIS6.7化探数据处理方法及异常图绘制[摘 要] 本文以苏丹红海州捷拜拉吐哈度57区金矿区水系沉积物测量中Au元素化探数据为例，介绍了MAPGIS6.7环境中化探数据的处理方法和流程，以及化探异常图的绘制。 [关键词] MAPGIS6.7 化探数据 处理方法 化探异常图 1. 前言 我队在苏丹红海州捷拜拉吐哈度57区金矿区开展了1∶5万～1∶20万水系沉积物测量，所有采样点位均使用GPS精确定位，共采集样品数量2200个，并由化验单位完成对化探样品Au元素含量的分析。 MAPGIS6.7是一个大型基础地理信息系统软件平台，它可以方便、快捷、有效的对经过化验分析的数据进行处理并自动圈定元素异常[1]。 2. 数据处理 2.1 原始数据整理 水系沉积物测量通过GPS对采样点进行定位，得到采样点坐标（X、Y）。在Excel中导入原始数据字段，包括采样点编号、X坐标、Y坐标、Au元素分析结果，然后另存为制表符分隔的文本文件，命名为“化探结果.txt”。 2.2 背景值和异常下限值的确定 异常下限主要根据剔除高值点后样品Au元素含量的算术平均值和标准离差来确定。其计算公式为： 其中， ——Au元素含量的异常下限值； ——剔除高值点剩余样品Au元素含量的算术平均值； ——剔除高值点剩余样品Au元素含量的标准离差值； 为经验系数，一般取1.65～3。本文 的取值为： =3。 统计计算时，利用Excel相关函数将Au元素含量所在的列进行从大到小排列，然后计算全部样品中Au元素含量的算术平均值和标准离差值，将大于 的含量数据剔除后，再次统计计算全部样品中Au元素含量的算术平均值和标准离差，直至将大于 的含量数据全部剔除后，用剩余样品中Au元素含量的算术平均值来度量Au元素的背景值，并以剩余样品的算术平均值加3倍的标准离差值来度量Au元素的异常下限值。表1为剔除高值点后Au元素的相关统计数据。 表1 剔除高值点后Au元素的相关统计数据 最大值 最小值 平均值 标准离差 异常下限 3.55 0.10 1.36 0.73 3.56 本文在绘制化探异常图时，将背景值 设定为1.36，异常下限值 设定为4.0。 3. 利用MAPGIS绘制化探异常图 3.1 数据的投影变换处理 投影变换处理就是将化探采样点坐标通过MAPGIS处理生成化探数据点位图。投影变换在MAPGIS6.7中的操作步聚： （1）在MAPGIS6.7下单击实用服务→投影变换→用户文件投影变换，如图1； （2）打开文本文件“化探结果.txt”，点击按指定分隔符→设置分隔符，将分隔符号中的“TAB键”打勾，“空格”中的勾取消,结果显示如图2； （3）设置X坐标、Y坐标所对应的的列号； （4）将用户投影参数、结果投影参数分别根据需求进行设置； （5）设置投影点图元参数，然后投影变换→确定，投影结果如图3； （6）保存点文件，命名为“化探结果.wt”。 3.2 数据属性连接 数据属性连接是将投影点位图与化探数据相关联，实现了将外部数据库的属性数据输入到MAPGIS图形文件的属性数据中[1]。属性连接在MAPGIS6.7中的操作步聚： （1）在Access中打开“化探结果.txt”→按步骤（默认选项）生成“化探结果.mdb” →在Access表文件中右击“化探结果”→导出（保存类型：dBASEⅢ（*.dbf）），生成“化探结果.DBF”； （2）在MAPGIS6.7下单击库管理→属性库管理； （3）在MapGis属性管理子系统中单击属性→连接属性； （4）打开连接文件“化探结果.wt”，打开被连文件“化探结果.DBF”（注：均按序号连接），结果如图4； 3.3 离散网格数据化 化探数据都是离散数据，为了绘制平面等值线图，需要对原始数据进行网格化处理，即生成GRD模型。MAPGIS6.7环境下生成GRD模型的步骤： （1）在MAPGIS6.7下单击空间分析→DTM管理，即打开MapGis数字地面模型子系统； （2）单击“文件”→打开数据文件→点数据文件，打开“化探结果.txt”； （3）单击“处理点线”→点数据高程点提取→选择Au元素含量值为高程值（即字段4）； （4）单击“Grd模型”→离散网格数据化，如图5； （5）根据需求设置相应参数，生成“化探数据.Grd”。 3.4 绘制平面等值线图 化探采样就是为了形成该地区各元素的化学异常图，圈定工作区的矿产异常范围，为后期地质找矿提供依据[2]。绘制平面等值线图就可以达到这个效果。MAPGIS6.7环境下生成平面等值线图的步骤： （1）单击“Grd模型