化探数据处理成图过程讲解.docx

PAGE PAGE 10 化探数据处理成图的过程 毕武 1、2 段新力 1、2 黄显义 1、2 袁小龙 1、2 彭仲秋 1、2 李永华 1、2 乌鲁木齐金维图文信息科技有限公司，新疆，乌鲁木齐，830091 新疆地矿局物化探大队计算中心，新疆，昌吉，831100 前言 GeoIPAS 软件用户群不断扩大，由于各用户对系统的熟悉程度不同，对软件功能了解不够，有必要分专题将 GeoIPAS 处理数据及成图过程做一系统总结，下面就化探数据处理成图的过程做一总结。 数据输入 数据预处理数据检查数据整理 原测数据 投影变换 数据网格化格式转换 数据处理 重磁数据处理电法数据处理 化探数据处理 数学地质数据处理水文地质数据处理 图件绘制 点位数据图点位符号图 剖面图 …… 等值线图场强图 三元图 三维立体图 输 出 打印输出 图像输出 文本文件图形文件 图 1 金维地学信息处理研究应用系统(GeoIPAS)图件制作流程 处理步骤 化探处理的成果包括：(1)参数统计表；(2)R 型聚类分析-谱系图；(3)重复样三层套合方差分析或者重复样合格率计算结果；(4)点位数据图；(5)地球化学图；(6)直方图；(7)组合异常图；(8)综合异常图；(9)远景区划图；(10)单元素异常参数统计(附表册)；(11)异常剖析(附图册)；(12)综合异常登记卡(附表册)。 在 GeoIPAS 系统中，化探数据处理分为以下几个主要步骤： 数据检查 数值检查，坐标检查，重复样坐标检查。 分析处理 重复样三层套合方差分析、重复样合格率计算、化探特征参数统计、化探背景值分析、R 型聚类分析、因子分析。 数据分析 数据变换；衬值、累加衬值；数据累加、累乘、比值；异常归一化。 网格化 离散数据网格化、XYZ 数据转网格数据 成图 点位数据图、彩色等量线图、直方图、组合异常图、单元素异常图、综合异常图、剖析图。 单元素异常参数统计 综合异常登记卡 开始 开始 数据输入 分析样品与坐标对应 TXT 数据预处理 TXT 数值检查 坐标检查 重复样提取 采样编码统计 TXT 数据分析处理 TXT 矢量 三层套合方差分析或 重复样合格率计算化探特征参数统计 R 型聚类分析 TXT 网格化 GRD 离散数据网格化 XYZ 数据转网格数据 TXT 点位数据图 直方图 成图 矢量 地球化学图 GRD 组合异常图 单元素异常图 综合异常图剖析图 TXT 矢量 异常统计 单元素异常参数统计 CSV 结束 图 2 化探数据处理成图流程 具体处理过程 数据检查 我们的数据处理工作从化验室提供的样品分析报告开始，项目要提供坐标和样品对应的分析数据，坐标我们一般取实际工作中的米单位，系统中默认东西向横坐标为X 坐标，不加带号，南北向纵坐标为Y 坐标，需要时还要提供样品对应的地质编码，我们拿到这个数据后首先进行数据检查，以确保数据中不出现写错、漏填、负数、 0、等字符，如果有这样的情况要找实验室给予纠正。数据准备好后，我们要把数据转换成 TXT 后缀的文本文件，这就做好了处理前的准备数据工作。 分析处理 重复样三层套合方差分析 一般是从分析样品的结果中挑出来重复样的分析值，每组四个样品，按如下顺序排列： 第一次采样第一次分析 第一次采样第二次分析 第二次采样第一次分析 第二次采样第二次分析结果： 三层套合方差分析成果--元素:Cu 方差来源 离差平方和 自由度 平均离差平方和 F-检验值 地点间 13.224 18 0.735 11.152 样品间 1.252 19 0.066 1.482 分析间 1.689 38 0.044 三层套合方差分析成果--元素:Au 方差来源 离差平方和 自由度 平均离差平方和 F-检验值 地点间 0.573 18 0.032 10.781 样品间 0.056 19 0.003 0.483 分析间 0.232 38 0.006 查 表 ：Fa0.05(18,19)=1.35 Fb0.05(19,38)=1.22 Cu：11.152 1.35 1.482 1.22 即：F1Fa F2Fb Au：10.781 1.35 0.483 1.22 即：F1Fa F2Fb 说明：一般是 F1＞临界值，F2＜临界值。 Faf0.05()理论值,说明采样与分析误差不显著，此次分析结果可靠， Fbf0.05()理论值,说明采样与分析误差显著，且以采样误差为主。 例如报告中的描述： 先列一个表格说明所有元素的 F1 和 F2； 重复样三层套合方差分析，39 种元素 F1 值远超过临界值，说明采样与分析误差变化均低于区内地球化学自然变化，数据是可用的。F2 值表现为两种情况， Ag、Pb、P、U、MgO、La 等 6 个元素 F2 值小于临界值，说明采样误差不大于分析误差，