地球化学数据处理与图件编制方法流程编写人:刘红杰QQ:498236930内蒙古第三地质矿产勘查开发院一、指导思想成矿地质背景地球化学研究就是从地球化学特征出发，借助已建立的地球化学信息提取技术，充分利用地球化学调查所获得的海量数据信息，提取有关反应成矿地质背景条件的地球化学信息，并编制相应地球化学图及相应的推断解释图件，为资源潜力评价有关成矿地质背景的研究提供地球化学支撑。
二、工作内容（一）基础图件成矿地质背景条件的地球化学信息提取首先是要编制有关基础地球化学图件。
主要有：1. 39种元素（化合物）地球化学图2. 地球化学组合异常图3. 地球化学综合异常图（二）解释推断图件地球化学解释推断图件，内容包括：1. 地球化学推断解译地质图2. 地球化学找矿预测图三、工作方法（一）数据校正处理由于区域地球化学数据受地理景观、采样介质、分析手段的影响，不可避免的产生明显的系统误差，尤其是涉及到区域性的化探数据，这种误差更为突出。
因此，在各省进行数据处理与专题地球化学图编制之前，有必要分别对各元素进行系统误差的处理，以便能更好地反映地质现象和矿产信息。
误差处理主要针对图幅间（包括分析批次）明显的系统分析误差（必须处理）和地质景观环境差异影响解释的效果（根据解释的需要确定）。
1. 系统误差特征及处理原则（1）分析误差源，所展示的数据误差与周边数据值具有明显的台阶状。
（2）数据误差在空间上具有区域性特点，区域、图幅或分析批次。
（3）在数据值的分布上，掩盖了地球化学特征和地质特征展布的延续性和规律性。
（4）在数据处理方法上，尽可能地选择线性校正，通过简单的计算可以复原数据。
2. 系统误差处理步骤（1）按原始点位采用符号分级的方式生成元素的符号图，分级方法采用累计频率方式。
（2）通过校正图示窗浏览原始数据全图，确定具有明显的数据台阶区域，区域的确定原则是由区域－＞图幅－＞批次；采用图形编辑工具，在图上直接圈定要处理的区域（用面的方式表示）。
（3）确定局部图幅和分析批次范围产生的系统误差，校正单元由系统提供的工具直接在显示窗中勾绘，确定完所有需要校正单元，各校正单元的ID需设定为唯一。
（4）建立校正单元与处理数据表空间位置索引关系。
（5）确定各单元的校正值或校正系数，主要方法是与单元周边数据进行对比分析，部分规律性较复杂的单元可以通过统计规律确定，同时还需考虑地球化学分布的整体空间分布趋势和地质背景；计算方法推荐采用：V ai=A V i+B其中：V ai为校正点校正后数据，A为校正系数，V i为校正点原始数据，B 校正常数。
A与B值的确定参照校正单元周边数据单元（正常的数据单元）确定，可以由经验判断，也可通过统计规律确定。
（6）数据校正，可采用SQL语言操作模式或应用软件系统提供的专用工具，按确定的校正值对各校正单元逐一进行计算。
元。
（5）单元校正计算，以基准单元为基础，应用归一化方法调整其他单元的数据值。
计算方法推荐采用：V aij=V ij×(V j/V r)其中V aij为第j个单元的第i个数据校正结果值，V ij为第j个单元的第i个原始数据，V j为第j个单元的平均值，V r基准单元的平均值。
V j和V r也可分别换为计算单元和基准单元的中位数。
（6）（7）利用校正计算结果重新生成符号分级图。
（8）观察全图，对部分校正结果不理想的单元，可通过上述步骤，对单元和校正值进行调整，并重新计算，直到校正数据和成图效果符合全局规律为止。
3. 地质与景观区域校正步骤（1）确定校正单元，选择地质单元（根据解释的需要，原则上可划分到系或组）或景观单元图层，编制单元唯一编码(ID)。
（2）建立校正单元图层与处理数据表空间位置索引关系。
（3）按单元统计计算数据特征值，主要应包括如下项：数据点数，最大值、最小值、平均值、中位数、几何均值，方差，累频分段值（25%、75%、85%、95%可选项）。
当各单元数据离散程度较高，尤其是突变点影响较大时，统计中应进行异点剔除处理，采用剔除方法：Vs+2.5×Dev，其中Vs为单元平均值，Dev单元标准差。
（4）确定基准单元，根据解释或综合异常编图的需要，选择一个单元作为校正的基准单利用校正后的数据制作符号分级图，（7）观察全图，对部分校正结果不理想的单元，可通过上述步骤，调整校正参数，并重新计算。
（二）坐标投影变换空间坐标转换是区域地球化学数据库建立和数据处理分析及编图重要环节。
各省应参考成图范围、成图比例尺及数据模型规范中《空间坐标系统及其参数规定分册》的相关规定选择投影坐标系统。
各投影模型椭球参数建议采用西安80/IUUG 1975坐标系(即长半轴：6378140米，短半轴：6356755.3041米)。
所收集的各类数据和基础图件（底图）均应按上述转换后，方能进行数据处理分析与编图。
转换方式分为表数据转换和图形数据转换，表数据转换可直接通过操作纵横坐标项（如经度、纬度）应用数字计算方法进行，图形数据的坐标转换需要借助于图形处理或GIS软件。
基本方法包括：（1）确定数据源空间坐标类型，收集投影参数（如中央经线、标准纬线，参考零纬线，椭球参数等）。
（2）确定转换的空间坐标类型及投影参数。
（3）对于表数据转换需要正确地选定横坐标与纵坐标数据项（或字段），建议转换后的坐标保存到新的字段中。
（4）注意原坐标和转换后的坐标单位推荐软件：GeoExpl，MapGIS（三）数据网格化离散数据网格化处理是空间数据插值的一种，即把无规则分布的空间数据内插为规则分布的空间数据集。
网格数据是编制地球化学图件的重要数据源。
网格化处理一般包括这样几个过程：①空间几何属性的确定；②插值方法（模型）的选择；③空间数据的探索分析，包括对数据的均值、方差、协方差、独立性和变异函数的估计等；④插值方法评价；⑤重新选择内插方法，直到合理。
网格化数据处理中要确定主要参数包括：1. 网格距：根据采样密度确定,一般网格距应与采样密度一致；2. 数据搜索半径：一般选择网格距的2.5倍；3. 数据计算模式：最近点或距离指数加权；参数数据的选择可根据不同研究目的改变。
推荐软件系统：GeoExpl, Surfer（四）基础图件编制1. 单元素地球化学图编图要素（1）简化地理图（以素图的方式表示）；（2）地质矿产简化图（以素图的方式表示）；（3）单元素按含量高低分级，建议采用累频分级方式，分19级，分级频率：0.5-1.2-2-3-4.5-8-15-25-40-60-75-85-92-95.5-97-98-98.8-99.5-100。
分级色阶选取：以冷色调（兰色）作为低值区，随着数据增大，颜色变暖，由兰－绿－黄－红－深红。
各省可根据本省元素地球化学特征和推断解释的需要作适当调整。
推荐软件系统：GeoExpl，Surfer图件样例：图一单元素地球化学图样式1422. 组合异常图编图要素（1）简化地理图（以素图的方式表示）；（2）地质矿产简化图（以素图的方式表示）；（3）元素按统计规律确定异常下限，分为外、中、内级，分别按累频的85-95.5-98划分。
一张图上可选择3-5 个元素，以异常下限为基准，可划分2-3条浓度带，每种元素通过不同颜色区别。
也可采用采用39种元素参与的主因子分析方法得到的因子得分图作为组合异常图。
3．因子分析及因子得分图（1）因子分析利用元素原始数据集对元素进行因子分析，取特征根较大且累计贡献率达80%～85%以上的前若干因子选做主要因子。
对选定前若干因子模型做最大方差旋转，确定主因子模型。
推荐软件系统：GeoExpl，SPSS，SASGeoExpl软件应用示范：在GeoExpl软件的“数据处理分析”模块中，选择“多元分析”—＞“因子分析”功能（如图二）。
在多元素数据表中，选择分析元素（不少于3个，可根据研究目的选择元素组合。
不要选择元素以外的变量，如坐标、ID等），设定好横、纵坐标，点击“特征值输入”设置分析结果的保存位置，最后点击“因子分析”，系统将自动计算因子元素特征值，即元素特征根和特征根百分比。
运算完成后将弹出“确定因子个数”的窗口，推荐按特征根累计百分比大于80%的个数输入（如图三）。
图二因子分析模块图三确定因子个数确定后软件将自动进行因子分析处理计算，计算结果将输出相关矩阵、特征根、特征向量、初始因子矩阵和旋转因子（载荷）矩阵，以及初始因子得分和旋转因子得分表。
（旋转因子与初始因子的差别主要在于旋转因子使复杂的矩阵变得更简单，分组更明显，易于地质解释，因此一般采用旋转因子的结果进行解释分析）。
分析结果将保存为ftr文件，该文件可以用记事本等文本工具直接打开（如图四）。
图四因子分析参数结果将“旋转因子矩阵”导入Excel表中得到更为直观的结果（如图五）。
图五旋转因子矩阵利用因子载荷矩阵，确定各因子的组合，一般来说，对应元素因子载荷越大，对这个因子的贡献也越大，并说明包含该元素的信息量越多。
按照载荷值大小对元素排序，一般取载荷绝对值>0.5或>0.6的元素作为该因子元素组合。
（2）编制主因子得分图应用对应每一个数据点的因子得分制作各因子二维空间得分图，以便于进行地质解释。
推荐软件系统：GeoExpl，MapGIS，SPSS，SASGeoExpl软件应用示范：GeoExpl中，进行因子分析后将会自动生成两个数据表，其中，后缀为“_FI”的，是对应每一个数据点的初始因子得分表，后缀为“_FR”的，是对应每一个数据点的旋转因子得分表（如图六）。
图六因子得分表可利用旋转因子得分表进行离散数据网格化后编制等值区图（如图七），以得到各主因子的因子得分图，作为组合异常图（如图八），该异常图的元素组合即是各因子中载荷较高的元素组合。
图七因子得分表数据网格化4. 综合异常图编图要素（1）简化地理图（以素图的方式表示）；（2）地质矿产简化图；（3）按多元素的综合方法确定综合异常。
综合异常图是通过应用数据处理结果（如多元素叠加、空间分析）和对多元素异常综合解释按照相关规则圈定的综合性图件。
综合地球化学异常图是圈定预测区的基础。
可以根据不同的组合方式，构成不同表达形式的综合异常图（图九）。
·多元素异常空间逻辑与的叠加，确定异常元素的异常下限，按元素最大空间范围圈定异常。
·多元素异常空间逻辑和的叠加，确定异常元素的异常下限，按元素交积圈定异常。
·累加、累乘元素组合叠加，根据元素的组合关系，按累加累乘的结果，圈定异常，异常的表达方式可以参考单元素异常的表示方法。
图八因子得分图148图九综合异常图149·多元素归一（标准化）加权叠加，选择相关元素分别进行归一化处理，根据元素的主次关系，分别赋以权重计算。
·多元统计分析处理综合，应用多元统计数据处理方法，如因子分析的因子得分，二维聚类分析的分类值，表达方式可以根据不同的处理方法和要突出的地质矿产特征，参考组合异常图的表达方式。