化探数据处理方法

内蒙古扎赉特旗东芒合矿和哈拉街吐矿

化探数据处理及图件编制方法

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 Ti

2 矿区地球化学特征研究的数据处理(以哈拉街吐为例)

⑴作出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由区域构造、地层、岩浆岩等因素引起，反映了区域因素引起的元素的含量变化，这种变化往往代表来自多重母体不同背景的起伏变化，实际上反映了元素含