化探数据处理成图过程精要.doc

化探数据处理成图的过程 毕武1、2　段新力1、2　黄显义1、2　袁小龙1、2　彭仲秋1、2　李永华1、2 1.乌鲁木齐金维图文信息科技有限公司，新疆，乌鲁木齐，830091 2.新疆地矿局物化探大队计算中心，新疆，昌吉，831100 0 前言 GeoIPAS软件用户群不断扩大，由于各用户对系统的熟悉程度不同，对软件功能了解不够，有必要分专题将GeoIPAS处理数据及成图过程做一系统总结，下面就化探数据处理成图的过程做一总结。 1 处理步骤 化探处理的成果包括：(1)参数统计表；(2)R型聚类分析-谱系图；(3)重复样三层套合方差分析或者重复样合格率计算结果；(4)点位数据图；(5)地球化学图；(6)直方图；(7)组合异常图；(8)综合异常图；(9)远景区划图；(10)单元素异常参数统计(附表册)；(11)异常剖析(附图册)；(12)综合异常登记卡(附表册)。 在GeoIPAS系统中，化探数据处理分为以下几个主要步骤： 1.1 数据检查 数值检查，坐标检查，重复样坐标检查。 1.2 分析处理 重复样三层套合方差分析、重复样合格率计算、化探特征参数统计、化探背景值分析、R型聚类分析、因子分析。 1.3 数据分析 数据变换；衬值、累加衬值；数据累加、累乘、比值；异常归一化。 1.4 网格化 离散数据网格化、XYZ数据转网格数据 1.5 成图 点位数据图、彩色等量线图、直方图、组合异常图、单元素异常图、综合异常图、剖析图。 1.6 单元素异常参数统计 1.7 综合异常登记卡 2 具体处理过程 2.1 数据检查 我们的数据处理工作从化验室提供的样品分析报告开始，项目要提供坐标和样品对应的分析数据，坐标我们一般取实际工作中的米单位，系统中默认东西向横坐标为X坐标，不加带号，南北向纵坐标为Y坐标，需要时还要提供样品对应的地质编码，我们拿到这个数据后首先进行数据检查，以确保数据中不出现写错、漏填、负数、0、等字符，如果有这样的情况要找实验室给予纠正。数据准备好后，我们要把数据转换成TXT后缀的文本文件，这就做好了处理前的准备数据工作。 2.2 分析处理 2.2.1 重复样三层套合方差分析 一般是从分析样品的结果中挑出来重复样的分析值，每组四个样品，按如下顺序排列： 11 第一次采样第一次分析 12 第一次采样第二次分析 21 第二次采样第一次分析 22 第二次采样第二次分析 结果： 三层套合方差分析成果--元素: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值小于临界值，说明采样误差不大于分析误差，以分析误差为主；其余33种元素F2值均大于临界值，说明误差应以采样误差为主。 2.2.2 地层编码统计 从样品相对应的地质图中自动读取地层的颜色号作为编码，实现每个样品的地层自动编码。 在GeoIPAS系统中的化探-化探特征参数统计-地层编码统计 其处理结果是在输入的数据文件中增加一列，过程是将坐标投影到地质图相应的投影参数上，以相同地层颜色来套取数据，取地层的颜色号为落入此地层的采样点的地层编码，所以要想一次做好，要先对地质图做个预处理，最好是按编码顺序，将颜色号改成相应的顺序号，比如Q的颜色改为1，N的颜色改为2，E的颜色改为3，等等。相对应的出来的结果就是按顺序的地层编码了，可以直接用了，绘制的直方图也是这个顺序了。 2.2.3 化探特征参数统计 化探参数统计成果--元素:Cu 地质单元 样品数 (N) 平均值(X) 中位数(M) 标准离差(S) 变化系