区化数据预处理

化探资料数据准备与预处理向运川区域地球化学研究课题组

2007年6月全国矿产资源潜力评价主要内容资料与数据准备数据预处理应用目标

应用区域地球化学进行矿产资源预测研究是以1:20万和1:50万水系沉积物资料为基础，收集局部研究区域的1:5万水系沉积物和1:1万土壤测量数据，以及岩石地球化学资料等，以成矿地质理论和地球化学理论为指导，对数据进行处理、分析、成图、推断解释、预测靶区圈定、及定量预测分析等。针对《全国矿产资源潜力评价》项目目的，结合区域地球化学的特点，重点对铜、铅、锌、镍、钨、锡、金、铬、钼、锑、稀土、银、硼等矿种进行预测分析、并初步估算单矿种资源量。数据源1:20万或1:50万水系沉积物资料与数据1:5万水系沉积物资料与数据1:5万－1:50万岩石测量资料与数据1:5000－1:2.5万土壤测量资料与数据化探异常检查与查证报告与数据其他相关地质、矿产、地理景观等资料数据分类采样信息数据元素分析数据质量监控数据图形数据文档报告

数据准备采样信息数据▲采样记录卡->数据表（Access或Excel）数据结构与内容参考《区域地球化学勘查规范》及相关标准

▲组合样信息->数据表▲采样点位布置图->数字化（MAPGIS）

水系沉积物采样点分布信息组合样信息数据准备元素分析数据▲数据表存储，采用Access或Excel

▲空间坐标信息，经度、纬度（十进制度），直角坐标（公里）

▲地质等编码信息（地层、岩石、景观区、构造单元、成矿单元等）

▲分析元素值，单位：Aung/g,氧化物％，其他ug/g

数据准备质量监控数据

▲分析元素信息表－记录元素分析方法

▲标样推荐值

▲标样分析结果

▲标样控制样质量控制参数

数据准备图形数据

▲基础地质、地理

▲采样点位分布图

▲异常图

▲异常查证分布图

采用MAPGIS或其他数字化软件分层数字化

数据准备文档资料

▲地球化学勘查报告

▲异常查证成果报告

▲省内有关地球化学与矿产研究报告及文章

▲填写各类报告主要内容信息表，包括：报告名称，提交单位，提交时间，报告编写主要人员，涉及范围，内容摘要，主要成果

条件许可的情况下将各类报告转换成Word文档保存

数据预处理空间坐标转换异点检查分析系统误差校正景观区校正

数据表数据预处理空间坐标转换▲根据各省确定的地理投影模型，对准备各类空间数据（数据表和图形数据）转换为统一的坐标系统，椭球参数建议采用西安80/IUUG1975坐标系(即长半轴：6378140米，短半轴：6356755.3041米)。

▲采用专业软件系统或GIS系统

1.表数据空间坐标转换

(1)确定源坐标系统、参数，、单位，以及坐标项

(2)确定目标坐标系统、参数、单位，给定新坐标项

源坐标系统目标坐标系统源坐标项目标坐标项数据预处理空间坐标转换

2.图形数据空间坐标转换建议在数字化完成后，进行此项转换

(1)确定源坐标系统、参数、坐标单位

(2)确定目标坐标系统、参数，坐标单位

源坐标系统目标坐标系统数据预处理异点检查▲检查空间坐标和分析元素的奇异值

▲推荐采用AccessSQL查询语言

如：SELECT索引编码,经度,纬度,AuFROMHHLLJJWHERE(Au<=0OrAu>2000）AndAu>-1E+30▲异点值替换，直接在表中替换或用相关软件

低于检出限的数据取检出限的1/2,如：<3.0=1.5

高于检出最高极限的数据取125%,如：>1000=1250

表名数据项条件数据预处理元素分析系统误差校正

由于区域地球化学数据受地理景观、采样介质、分析手段的影响，不可避免的产生明显的系统误差，尤其是涉及较大区域的化探数据，这种误差更为突出。因此，在数据处理之前，有必要分别对各元素进行系统误差的处理。▲校正由分析或采样等引起的系统误差

▲推荐采用GeoExpl

可采用系统的数据校正模块或二维空间分析子系统结合SQL数据分析功能可分别进行离散数据和网格数据的校正

▲建议直接对原始离散数据校正

数据预处理元素分析系统误差校正系统误差校正的原则：（1）分析误差源，所展示的数据误差与周边数据值具有明显的台阶状；（2）数据误差在空间上具有区域性特点，区域、图幅或分析批次；（3）在数据值的分布上，掩盖了地球化学特征和地质特征展布的延续性和规律性。（4）在数据校正方法上，尽可能地选择线性校正，通过简单的计算可以复原数据。▲误差识别方法

数据表浏览图形方式浏览数据点位图数据预处理元素分析系统误差校正▲校正步骤

(1)按原始点位采用符号分级的方式生成元素的符号图，分级方法采用累计频率方式。(2)通过校正图示窗浏览原始数据全图，确定具有明显的数据台阶区域，区域的确定原则是由区域->图幅->批次；采用图形编辑工具，在图上直接圈定要处理的区域（用面的方式表示）。(3)确定局部图幅和分析批次范围产生的系统误差，校正单元由系统提供的工具直接在显示窗中勾绘，确定完所有需要校正单元，各校正单元的ID需设定为唯一。

数据预处理元素分析系统误差校正▲校正步骤

(4)建立校正单元与处理数据表空间位置索引关系。

(5)确定各单元的校正值或校正系数，主要方法是与单元周边数据进行对比分析，部分规律性较复杂的单元可以通过统计规律确定，同时还需考虑地球化学分布的整体空间分布趋势和地质背景；计算方法推荐采用：Vai=AVi+B

其中：Vai为校正点校正后数据，A为校正系数，Vi为校正点原始数据，B校正常数。A与B值的确定参照校正单元周边数据单元（正常的数据单元）确定，可以由经验判断，也可通过统计规律确定。

(6)数据校正，可采用SQL语言操作模式或应用软件系统提供的专用工具，按确定的校正值对各校正单元逐一进行计算。

(7)利用校正计算结果重新生成符号分级图。

(8)观察全图，对部分校正结果不理想的单元，可通过上述步骤，对单元和校正值进行调整，并重新计算，直到校正数据和成图效果符合全局规律为止。

数据预处理景观区及地质构造分区校正

由于不同的地理景观和地质因素的影响，在进行大区域成图或异常研究中，使得低背景或弱小异常被掩盖，失去了许多有用成矿信息，这种校正的目的是突出弱小异常。

▲校正方法

通过圈定研究区内景观或构造单元，建立单元与数据表的空间关系，利用标准单元单元层为基准，对其他校正单元的数据进行归一化校正处理。▲方法原理

设有m个校正单元，第t个校正单元为标准单元，则：

其中，Cij’为第j个单元第i个样品校正后的值，Cij为第j个单元第i个样品原始值，nj为第j个单元的样品数，n为全区样品数J单元平均值总平均值或标准单元平均值数据处理与分析

景观区及地