内蒙古扎赉特旗东芒合矿和哈拉街吐矿化探数据处理及图件编制方法1 化探数据质量评价的数据处理(分矿区)⑴统计重采样和重分析抽查样所占样品总数的比例比例 = （重采样和重分析抽查样数/工作样总数）100%⑵作出SSPS数据文件将重采样和重分析样分别作成SSPS数据文件。
文件中列出项目为：①重采抽查样重采样号元素含量相应的工作样号元素含量②重分析抽查样重分析样号元素含量相应的工作样号元素含量⑶计算各元素相对误差重采样和重分析抽查样相对误差均按RE(%) = |C1-C2|/0.5×(C1+C2)×100%计算。
C1为重采样或重分析抽查样的分析含量C2为重采样或重分析抽查样的相应的工作样的分析含量| |为绝对值RE(%)≤30%为合格，＞30为超差（不合格）；（Au：RE(%)≤50%为合格，＞50为超差）⑷计算各元素的合格率η= （抽查样品中合格的样品数/抽查样品的总数）100%合格率（η）应＞80%，即这批样品的分析结果是可信的。
⑸列表表示检查或分析质量结果表××化探重采样抽查各元素的合格率(%)Cu Pb Zn Cr Ni Co Sn V Ag Ti2 矿区地球化学特征研究的数据处理(以哈拉街吐为例)⑴作出SSPS数据文件作出下列SSPS数据文件：①文件1：整个矿区数据文件；②文件2：矿区地层数据文件；③文件3：矿区岩浆岩数据文件；④文件4 ：下二叠统大石寨组（P1d）数据文件；⑤文件5 ：下白垩统大磨拐河含煤组（K1d）数据文件；⑥文件6 ：华力西晚期侵入岩数据文件；⑦文件7 ：燕山期早期侵入岩数据文件；⑧文件8 ：燕山期晚期侵入岩数据文件；⑨文件9：已知矿附近一定范围数据文件每一数据文件的内容项目包括：序号野外号 X坐标 Y坐标各元素的含量⑵整个矿区和各地质单元（各地层、各岩浆岩）样品各元素含量特征统计统计的参数包括：①元素含量平均值；②最大值；③最小值；④标准离差；⑤变化系数（标准离差/含量平均值）；⑥浓度克拉克值（元素含量平均值/该元素的克拉克值）整个矿区和各地质单元统计结果含量平均值、最小值、最大值用表表示。
⑶整个矿区和各地质单元样品各元素的概率分布特征统计①标准离差②峰度③偏度④概率分布曲线特征⑷矿区各地层样品各元素的局域丰度和蚀变-矿化叠加系数特征统计根据地球化学过程的基本定律(A.B.Vstelius,1960),一个矿区地层中元素的“丰度”应该是沉积岩沉积成岩时的初始平均含量，而不应包括后期岩浆、蚀变、矿化作用等地质作用造成的元素含量的增赢或亏损。
而矿区内局部地区地层中元素的“局域丰度”，至少应排除最后蚀变-成矿作用叠加的那一部分元素的含量。
若本区各地层中元素概率分布及其偏度和峰度特征表明元素呈偏对数正态分布。
这说明地层中多数元素都受到了后期不同程度的蚀变-成矿作用的叠加。
据此，剔除了不服从正态分布的超差样品(即含量大于或等于元素的平均值加上2倍标准离差(c+2δ)的那些样品)后，再求出的元素含量平均数即为本区各地层中元素的“局域丰度”(x)。
该局域丰度可被认为是本区各地层在区域沉积作用（包括海底喷气喷流、热水沉积作用等）的元素平均含量。
而该区各地层元素的平均值（c）则是在局域丰度的基础上叠加了岩浆、蚀变-矿化作用后的平均含量。
用克拉克值除以局域丰度(x/克值)所得的浓集系数称为初始浓集系数（k 0），反映了沉积作用对元素的初步聚集程度；用局域丰度除以平均值(c/x) 所得的浓集系数称为局域蚀变-矿化叠加系数(k)，反映了本区地层形成后遭受的矿化强度。
因此k 0和k可作为评价地层含矿性的参数。
统计结果可用下表表示：表××哈拉街吐矿区下二叠统中微量元素的平均值（c）、局域丰度(x)、初始浓集系数（k 0）和蚀变-矿化叠加系数(k)-9-6表××哈拉街吐矿区下白垩统中微量元素的平均值（c）、局域丰度(x)、初始浓集系数（k 0）和蚀变-矿化叠加系数(k)-9-6⑸整个矿区和各地质单元样品各元素的相关关系特征统计作出整个矿区和各地质单元样品各元素的相关矩阵⑹整个矿区和各地质单元样品各元素的聚类特征统计作出整个矿区和各地质单元样品各元素的R型聚类图3 矿区单元素地球化学异常的数据处理和成图⑴背景值和异常下限值的确定地质地球化学观察法和统计法相结合确定。
①在全区数据表上观察，一般含量的中间值确定为背景值（C0）；②高于背景值3～5倍的数据确定为异常下限值(C a)；③统计某一元素≥异常下限值的样品的个数占总样品数的百分数，一般主要成矿元素异常下限值的样品所占的百分数为30%左右。
④根据元素的概率分布曲线特征，若分布曲线为对称对数正态分布，取正态分布曲线的中值作为背景值，取背景值＋2×标准离差(C0＋2δ)作为异常下限；若分布曲线为不对称对数正态分布，则将全区数据去掉大于标准离差的数据后，再作概率分布曲线，直至出现对称对数正态分布，这时再正态分布曲线的中值作为背景值，取背景值＋2×标准离差作为异常下限。
⑤根据地质地球化学观察法和统计法相结合正确确定背景值和异常下限值。
⑵编制全区单元素含量等值线图等值线间隔按实际情况适宜确定⑶编制全区单元素异常图将等值线的间隔设为异常下限值（1 C a）、2倍异常下限值（2 C a）、4倍异常下限值（4 C a），每间隔间充填适当的颜色，即为异常图。
⑷单元素异常图的整饰每张图可放置2～3个单元素异常图。
4 矿区综合地球化学异常的数据处理和成图⑴共生元素累加和异常图的编制①共生累加元素的选择综合矿区元素相关分析、聚类分析、因子分析等获得的信息及元素地球化学习性，选择有较密切共生关系元素的含量相加，作为综合指标圈定异常图，此时的背景值、异常下限值要重新求取。
②编图方法以求得的几个元素含量的和成图，方法同单元素成图，可作出综合等值线图或综合异常图。
⑵对抗元素（累加）商异常图的编制①对抗（累加）商元素的选择综合矿区元素相关分析、聚类分析、因子分析等获得的信息及元素地球化学习性，选择有明显对抗关系的元素，求其商，作为综合指标圈定异常图，此时的背景值、异常下限值要重新求取。
②编图方法以求得的元素含量的商成图，方法同单元素成图，可作出综合等值线图或综合异常图。
⑶单元素趋势剩余值异常图的编制①趋势剩余值数字化处理过程趋势剩余值数字化及成图方法简述如下：经趋势分析处理的采样点上元素含量Z将被分解成3个基本分量：趋势分量Z i、异常分量A I和随机分量L i，即Z = Z i＋T i或Z = Z i＋A I＋L i趋势分量Z i由区域构造、地层、岩浆岩等因素引起，反映了区域因素引起的元素的含量变化，这种变化往往代表来自多重母体不同背景的起伏变化，实际上反映了元素含量的背景部分。
异常分量A i由局域范围的矿化作用、特殊的岩性或地球化学异常引起。
在趋势分析中，用其剩余值（残差值）T i表示局部因素引起的元素含量变化。
剩余值T i 和趋势值不同，它的变化是不连续的，代表了元素的异常部分。
有时为了进一步突出异常，又将剩余值T i分解为异常分量A I和随机分量L i。
随机分量L i有可能由采样、样品加工、分析等偶然误差因素造成。
本次处理的数据主要在矿田范围内及其周边，目的是要突出地球化学异常，故采用趋势分析的剩余值成图。
趋势分析通常是在二维空间进行的。
每一采样点上元素的含量Z和地理坐标之间存在着函数关系，即Z = f(x,y)趋势分析的任务就是导出函数f(x,y)的一个多项式去拟合或逼近元素含量的观察值。
函数f(x,y)的一次、二次、三次等多项式（趋势面方程）为：Z i = a0+a1x i+a2y i+a3x i y iZ i = a0+a1x i+a2y i+a3x i y i+a4x i2+a5y i2Z i = a0+a1x i+a2y i+a3x i y i+a4x i2y+a5xy i2+a6x i3+a7y i3式中，Z i是各采样点或元素含量Z的预测值，或者说用Z i去拟合实际的元素的含量Z，Z i 就是所谓趋势分量。
a0、a1、a2、a3、a4、a5……等待定系数须用最小二乘法原理求取，最小二乘法可使多项式中的误差波动达到极小，也就是使趋势值Z i更好地逼近元素的含量Z。
求取各采样点元素含量的趋势值Z i后，即可用每点的实际元素含量Z减去该点的趋势值Z i求得每点的剩余值T i。
即T i = Z - Z i剩余值为负值的，计为零。
根据每点所求得的剩余值T i即可绘出剩余值等值线图。
趋势面方程阶数的选择在数学理论上是应以求出的趋势值Z i最为逼近观察值Z，即拟合度最大为宜。
拟合度以趋势面上的变差占总变差的百分比来表示。
但在实际成图过程中，拟合度的选择要充分考虑化探工作的目的和工作区地质构造及矿化分布等情况，拟合度过高，剩余异常就少了。
一般区域化探拟合度选择为20%。
由于本次工作的目的是突出矿区的剩余异常，拟合度的选择为15%左右。
②编图方法以求得的每个元素含量的趋势剩余值T i成图，方法同单元素成图，作出综合等值线图即可。
⑷单元素值因子得分异常图的编制①因子得分数字化处理过程因子得分数字化及成图方法简述如下：化探样品的多元素分析数据（变量）之间存在一定的相关关系，因此有可能用较少的综合指标反映各变量中的各类信息，而各综合指标之间彼此不相关，就是说其所代表的信息不重叠，这种数理统计方法称为因子分析。
代表各类信息的综合指标称为因子或主成分。
具体的程序是：①通过对原始数据变量（元素含量）的相关分析结果，进行初始因子的提取并获取碎石图信息和因子负荷矩阵；②在此基础上选择一定的方法进行因子旋转并获取旋转后的因子负荷矩阵和因子转换矩阵；③建立起因子得分系数矩阵和回归因子分数的协方差矩阵（以考察旋转后的因子是否仍然正交）；④根据因子得分系数导出每个因子得分的计算方程式；⑤将每一采样点的各个原始变量（元素含量）代入因子得分的计算式，求出每一采样点各因子得分值；⑥根据每一采样点各因子得分值编制因子得分等值线图。
如某矿区次生晕数据因子分析根据原始数据相关分析结果，初始因子采用主成分分析法提取公因子，初始分析结果和因子碎石图（图1）表明：第１初始因子和第２初始因子之间、第６初始因子和第７初始因子之间特征值差值较大，第3、第4、第5、第６初始因子相互之间特征值差值较小，因此可以初步得出提取３个综合因子将能得出概括出各元素数据所反映的大部分的信息的认识。
因子旋转使用Varimax法（方差最大旋转法），经３次迭代收敛，得出的旋转后的因子负荷系数矩阵（表1）表明：经旋转后３个综合因子负荷系数对各原始变量（元素）的绝对值有较大的差别。
第1因子（F1）对Pb、Zn有绝对值较大的负荷因子，第2因子（F２）对As、F有绝对值较大的负荷因子,第3因子（F3）对Sn有绝对值较大的负荷因子。