化探讲座数据处理2012湖北

1、注：多目标样注：多目标样16或或4平方千米组合样分析，其它一般分析原始单点样平方千米组合样分析，其它一般分析原始单点样 选择研究对象中一切能提供地质、找矿选择研究对象中一切能提供地质、找矿信息的，并能直接或间接测定的地球化学变量均可作为化探信息的，并能直接或间接测定的地球化学变量均可作为化探指标。一般几种至数十种（含指标。一般几种至数十种（含经数据处理后的各种组合变量）经数据处理后的各种组合变量）其中，其中，Cij为第为第j个单元第个单元第i个样品校正个样品校正后的值，后的值，Cij为第为第j个单元第个单元第i个样品原个样品原始值，始值，nj为第为第j个单元的样品数，个单元的样品数，n为为全区

2、样品数全区样品数总平均值或标准总平均值或标准单元平均值单元平均值J单元平均值单元平均值浓集克拉克值某地质体平均值/克拉克值相对浓集系数某地质体平均值/区域背景 网格距：根据采样密度确定网格距：根据采样密度确定, ,一般网格距应与采样密度一一般网格距应与采样密度一致致 数据搜索半径：一般选择网格距的数据搜索半径：一般选择网格距的2.52.5倍倍 数据计算模式：最近点或距离指数加权数据计算模式：最近点或距离指数加权 参数数据的选择可根据不同研究目的改变参数数据的选择可根据不同研究目的改变 网格结点网格结点网格距网格距原始数据点原始数据点网格化的几何关系网格化的几何关系网格化处理网格化处理 （1）连

3、续性、均匀性、渐变性（2）峰值（极大值）、平均强度、衬度KCa/C0(C/Ca、Cmax/Ca) （3）线金属量Ml=x (Cx-C0 )面金属量Ms=s (Cx-C0 ) (4) 形态描述（线状、带状、等轴状、不规则状）与地质体关系 (5) 异常的元素组合特征指示意义，例如：AsSbBiHgAg，WMoSnCu 组合等（6）分带性特征：元素性质分带、浓度分带（水平与垂直）（7）异常剥蚀度系数E=（Kele21Kele22Kele2n）/（Kele11Kele12Kele1n）其中Kele1n、Kele2n分别为典型的前缘元素和尾晕元素的异常衬度。（8）不同介质中某元素含量（平均值）比值系数，

4、如：Ki=某元素在土壤中的平均含量/该元素在岩石中的平均含量反映在表生地球化学作用下元素相对于不同地层的贫化与富集程度，用于异常评价和研究元素在不同表生条件下（不同地质背景、不同气候、不同氧化还原环境等）的迁移、富集规律。（8）金属元素特征比值，如：Au/Au+Ag、Cu/Pb+Zn、Ni/Co、Fe/Mn、Ba/Sr、K+Na/Ca+Mg、K+Na/Al等等，表征元素在不同地质背景条下的相对富集与分配规律(9) 地球化学元素特征组合参数及其意义规格化面金属量（NPA）法（谢学锦，1980）NPA值=异常区某元素的平均含量/该元素的异常下限*大于异常下限的异常面异常面积NPA值越大找矿希望越大

5、。还可以计算元素组或主要成矿元素的总的NPA值：分带序列参数 其中：PK1为主成矿元素NPA值之和；PK2为前缘元素NPA值之和；PK3为尾晕元素NPA值之和；在无已知分带模式时，前缘、后尾元素的确定可依据邵跃（1984）的方案或经验分带序列。即：前缘元素为As2、Sb、Ba、Hg、Sr，尾部元素为W 、Mo、Sn、Bi、Ni、Be。该方法仅用于热液矿床或热液再造型矿床。均值（原始数据直接计算）均值（原始数据直接计算）平均含量大小平均含量大小高差（原始数据直接计算）高差（原始数据直接计算）相对平均值的离散程度（反映成矿的可能）相对平均值的离散程度（反映成矿的可能）背景值（剔除所有野值后计算，一

6、般此时为正态分布，符合概率统计概念）背景值（剔除所有野值后计算，一般此时为正态分布，符合概率统计概念）- -真实背景大小真实背景大小背景离差（剔除所有野值后计算）背景离差（剔除所有野值后计算）计算异常下限需要计算异常下限需要变异系数变异系数离差离差/ /平均值，越大更易成矿，一般用大于平均值，越大更易成矿，一般用大于1 1判别，如大于判别，如大于5 5肯定可肯定可成矿。成矿。衬值衬值原数据原数据/ /（背景值、异常下限、同类岩石（背景值、异常下限、同类岩石）-比较值比较值异常强度异常强度最大值最大值/ /异常下限异常下限面金属量面金属量平均值平均值* *面积，成矿规模大小面积，成矿规模大小NA

7、PNAP值值衬值衬值* *面积，不同异常间相加或比较面积，不同异常间相加或比较外、中、内带外、中、内带异常下限的异常下限的2n(n=02n(n=0、1 1、2 2或其它或其它) )，平面分带性，平面分带性相关系数相关系数相关程度，用临界相关系数判别相关程度，用临界相关系数判别1、元素之间的线性累加（线性相关分析、极差化处理）累加、乘积、比值等2、R、Q型聚类3、因子分析4、判别分析。灰关关联滤灰关关联滤相似类比相似类比人工经网网人工经网网自适应学习识别自适应学习识别小滤波技小滤波技高频若信息提取高频若信息提取模模算法模模算法隶属分类隶属分类遗算法遗算法特征归属分类特征归属分类1 1、采样点位图

8、、异常查证实际材料图（原始数据图）、采样点位图、异常查证实际材料图（原始数据图） 2 2、地球化学图、地球化学图 等值线勾绘：等值线勾绘：0.1lg0.1lg间隔（方差）法间隔（方差）法 累积频率分级法累积频率分级法3 3、单元素异常或组合（地球化学）异常图、单元素异常或组合（地球化学）异常图 单元素异常单元素异常3 3级，不同异常元素的套叠（每元素一条异常下限，重点成矿级，不同异常元素的套叠（每元素一条异常下限，重点成矿元素可做面色分级）形成组合异常元素可做面色分级）形成组合异常4 4、综合异常图、综合异常图 不同异常元素并（也可用交）集的综合（一般一组有成因联系的一组元不同异常元素并（也可

9、用交）集的综合（一般一组有成因联系的一组元素圈一条线，同组异常间无相交），可手工勾绘素圈一条线，同组异常间无相交），可手工勾绘 ，也可计算机绘制。，也可计算机绘制。 必须对异常进行编号并做地质与价值分类必须对异常进行编号并做地质与价值分类5 5、推断解释图、推断解释图 推断地质构造图、成矿预测图、潜力评价图推断地质构造图、成矿预测图、潜力评价图 （规范未要求的，有条件或（规范未要求的，有条件或必要时做）必要时做） 采样点位图采样点位图 主要内容包括水系、主要居民点、主要地物标志、交通道路、高斯方主要内容包括水系、主要居民点、主要地物标志、交通道路、高斯方里网格、经纬度坐标、采样点位置及采样信息

10、和分析数据信息属性。重里网格、经纬度坐标、采样点位置及采样信息和分析数据信息属性。重复采样、质量检查点以不同颜色或符号标注。比例尺为复采样、质量检查点以不同颜色或符号标注。比例尺为1异常查证实际材料图异常查证实际材料图 是在异常查证工作完成后编制。一般宜在地质矿产底图或综合异常图是在异常查证工作完成后编制。一般宜在地质矿产底图或综合异常图投放完整的异常查证实际材料，包括可能的面积性加密样点，异常查证投放完整的异常查证实际材料，包括可能的面积性加密样点，异常查证地质地质- -化探剖面，查证时随机采集的各类样品点位，槽探、坑探工程位置化探剖面，查证时随机采集的各类样品点位

11、，槽探、坑探工程位置等。比例尺为等。比例尺为11、地球化学图等量线间隔选取 a、采用0.1lg（g/g和ng/g）含量间隔(表7）。当数据为异常含量时，可适当将等量线抽稀为0.2 lgg/g (或ng/g)或更大，使等量线图面上间距不小于0.7mm。 b、采用累积频率的分级方法，推荐： 0.5、1.2、2、3、4.5、8、15、25、40、60、75、85、92、95.5、97、98、98.8、99.5分级间隔对应的含量作等量线勾绘，可根据人作区数据的分布特征进行调整，使地球化学图能反映区内地质、矿产分布特征。采用累积频率分级法确定等值线间隔应对元素的测试精度分析，以

12、保证等值线的可靠性2、编图要素： 简化地理图 、地质矿产简化图 、单元素按含量高低分级、全区或分地质构造、景观单元数据统计特征值3、表达方式：图名： xxx预测区x元素（或化合物）地球化学图 、地理图以素图的方式表示、 地质矿产图以素图的方式表示4、 等值区图例, 元素特征值表或图 XXXXXX灰区设置灰区设置: : 分级灰阶的选取方式为：以冷灰调（兰灰）作为低值区，随着数据分级灰阶的选取方式为：以冷灰调（兰灰）作为低值区，随着数据的增大，颜灰变暖，即由兰绿黄红深红变化，的增大，颜灰变暖，即由兰绿黄红深红变化，各灰区内不同含量各灰区内不同含量线间隔还可以用过渡灰阶表示。线间隔还可以用过渡灰阶表

13、示。含量坐标取对数，其组距可采用0.1 lgg/g或ng/g。，组端值规定小数点后第二位数字为7(对数负值为3) 责任表中应增加样品分析单位、分析栏 第一次第一次修改后修改后等值线等值线加密，加密，背景区背景区范围合范围合理加大理加大甲地地球化学图实例乙地地球化学图的绘制实例成图美观勾绘合理，但数据位未预处理成图美观勾绘合理，但数据位未预处理（ng/gng/g较合理较合理) )，并以一定的含量等级划分出异常浓度外带、，并以一定的含量等级划分出异常浓度外带、中带和内带。当数据为单一地质单元或数据基本符合正态分布时，按照全区中带和内带。当数据为单一地质单元或数据基本符合正态分布时，按照全区异常下限

14、值的异常下限值的1 1、2424、3838倍划分三个浓度带，勾绘异常外、中、内带。也倍划分三个浓度带，勾绘异常外、中、内带。也可采用累频方式，取可采用累频方式，取85%85%频数的值作为异常下限值，采用频数的值作为异常下限值，采用90%90%、95%95%频数值将频数值将异常划分为外、中、内异常划分为外、中、内3 3级浓度分带。异常内、中、外带含量等级含量可根级浓度分带。异常内、中、外带含量等级含量可根据人作区实际情况，进行调整。据人作区实际情况，进行调整。，当数据为两个以上的多个母体分布时，可依多，当数据为两个以上的多个母体分布时，可依多个母体分别统计出的衬值，统一确定衬度异常下限，勾绘异常

15、。并以衬度值个母体分别统计出的衬值，统一确定衬度异常下限，勾绘异常。并以衬度值1 1、2424、3838划分三个异常浓度分带。划分三个异常浓度分带。，以求得普查区区域背景趋势面，然后用原始数，以求得普查区区域背景趋势面，然后用原始数据减去该点的背景趋势值所获得的剩余异常值作图。当求得剩余异常值后，据减去该点的背景趋势值所获得的剩余异常值作图。当求得剩余异常值后，以正剩余异常值确定剩余异常下限，并以一定的倍值量勾绘正剩余异常及浓以正剩余异常值确定剩余异常下限，并以一定的倍值量勾绘正剩余异常及浓度分区。在应用剩余异常值制作异常图时，应对求取地球化学背景趋势面的度分区。在应用剩余异常值制作异常图时，

16、应对求取地球化学背景趋势面的方法进行试验、研究，以确保圈出的异常具有明确的找矿指示意义。方法进行试验、研究，以确保圈出的异常具有明确的找矿指示意义。无论采用上述哪种方案制作异常图，都应在异常图上表明异常图作图方无论采用上述哪种方案制作异常图，都应在异常图上表明异常图作图方法、作图参数和异常图例。法、作图参数和异常图例。 内蒙甲内蒙甲鄂西乙鄂西乙不足：未确定一个主要元素，选择不足：未确定一个主要元素，选择高于异常下限的高于异常下限的3 3个浓度带，用面个浓度带，用面灰表示，分别为橙灰、浅红、红灰。灰表示，分别为橙灰、浅红、红灰。 1、根据地球化学子区分别确定异常下限所圈定根据地球化学子区分别确定

17、异常下限所圈定的单元素异常，将空间上密切相伴的单元素异常，将空间上密切相伴( (叠合在一起叠合在一起) )、同种成因的所有元素的正异常，归并为一个综合异同种成因的所有元素的正异常，归并为一个综合异常。综合异常区域内形成正异常的所有元素的集合常。综合异常区域内形成正异常的所有元素的集合称为该综合异常的元素组合称为该综合异常的元素组合 研究异常元素组合类型及标型元素组合中各元研究异常元素组合类型及标型元素组合中各元素所占的地位，可以判断异常的成因、区分矿与非素所占的地位，可以判断异常的成因、区分矿与非矿异常、推断矿床类型、判断矿体剥蚀程度等重要矿异常、推断矿床类型、判断矿体剥蚀程度等重要的地质问题

18、。因此，异常元素组合类型的划分具有的地质问题。因此，异常元素组合类型的划分具有重要的实践和理论意义。重要的实践和理论意义。根据单元素异常在图面上的位置关系，构成元素组合的条件：根据单元素异常在图面上的位置关系，构成元素组合的条件：(1)(1)与主要元素异常位置重合或部分重合；与主要元素异常位置重合或部分重合；(2)(2)与主要元素异常紧密相邻，有其它同类元素异常同时与相邻与主要元素异常紧密相邻，有其它同类元素异常同时与相邻异常重合或部分重合异常重合或部分重合( (即可能处于同一矿田中，是元素水平分带即可能处于同一矿田中，是元素水平分带或次生作用产物或次生作用产物) )；(3)(3)这些元素能够

19、在相同地球化学作用过程中共同富集，即成因这些元素能够在相同地球化学作用过程中共同富集，即成因基本一致；基本一致；(4)(4)与已知矿床与已知矿床( (矿田、矿集区、成矿带矿田、矿集区、成矿带) )的地球化学异常元素组的地球化学异常元素组合一致或相近；合一致或相近；(5)(5)与不同侵入岩类、或特殊岩性、特殊地质作用的典型元素组与不同侵入岩类、或特殊岩性、特殊地质作用的典型元素组合一致或相近。合一致或相近。A A 元素按相对规模由大到小依次排列。元素按相对规模由大到小依次排列。B B 在顺序排列的元素符号间以英文连接符在顺序排列的元素符号间以英文连接符( (短横短横 杠杠) )相连。相连。C C

20、 主要成矿元素加下划线，并以推断的矿床主要成矿元素加下划线，并以推断的矿床( (组组) ) 类型的颜灰表示。类型的颜灰表示。D D 将元素浓度分带数以上标方式标注于元素符号右上角。将元素浓度分带数以上标方式标注于元素符号右上角。 如：如：Cr3-Co2-Ni3-Cu2Cr3-Co2-Ni3-Cu2将同一元素组合的所有元素异常共同覆盖的范围作为该综合异将同一元素组合的所有元素异常共同覆盖的范围作为该综合异常的范围常的范围, ,该范围的外边界线稍做圆滑作为综合异常的界线该范围的外边界线稍做圆滑作为综合异常的界线 综合异常边界线和范围用主成矿元素综合异常边界线和范围用主成矿元素( (或推断的矿床组类

21、型或推断的矿床组类型) )的花纹、颜灰来表示。的花纹、颜灰来表示。综合异常的面积按主异常元素中面积最大者计算。当综合异常综合异常的面积按主异常元素中面积最大者计算。当综合异常中某种元素的异常不止一处时，统计该元素异常地球化学参数中某种元素的异常不止一处时，统计该元素异常地球化学参数时这应对几处异常合并处理。时这应对几处异常合并处理。 综合异常按其中心的地理位置上从西到东、由北到南统一顺综合异常按其中心的地理位置上从西到东、由北到南统一顺序编号序编号( (即北部异常为小号即北部异常为小号) )。按人作区或人作图幅编号（连幅。按人作区或人作图幅编号（连幅也可统一顺序编号）。也可统一顺序编号）。A

22、A 元素相对规模大；元素相对规模大；异常平均强度异常平均强度 异常平均强度异常平均强度相对规模相对规模该元素异常面积。该元素异常面积。 该区该元素异常下限该区该元素异常下限B B 是推断矿化类型的主要成矿（指示）元素；是推断矿化类型的主要成矿（指示）元素；C C 是推断侵入岩类或特殊岩性、特殊地质作用的典型是推断侵入岩类或特殊岩性、特殊地质作用的典型指标元素。指标元素。原则上把所有综合异常编制在一张图上，若异常特别密集，图原则上把所有综合异常编制在一张图上，若异常特别密集，图面负担过重，可分编面负担过重，可分编CuCu多金属、贵金属多金属、贵金属2 2张（非金属组）张（非金属组） 1 1 图名

23、：图名： xxxxxx地区（幅）地球化学综合异常图地区（幅）地球化学综合异常图 2 2 地质矿产底图以简化素图的方式表示。地质矿产底图以简化素图的方式表示。 3 3 把与研究确定的矿种把与研究确定的矿种( (组组) )有关的综合异常标示于底图上；有关的综合异常标示于底图上； 4 4 将综合异常编号将综合异常编号- -异常类别标注于综合异常范围内，二者间异常类别标注于综合异常范围内，二者间用英文连接符连接；当异常区内无法标注时，可以在其外适当用英文连接符连接；当异常区内无法标注时，可以在其外适当位置标示，用引线引向异常区；位置标示，用引线引向异常区； 5 5 将元素组合标注于综合异常范围内将元素

24、组合标注于综合异常范围内, ,在编号下方。元素组合在编号下方。元素组合可换行表示可换行表示 6 6 在地质矿产图例的基础上增加地球化学异常的图例。包括在地质矿产图例的基础上增加地球化学异常的图例。包括主要异常元素及划纹、颜灰、异常编号主要异常元素及划纹、颜灰、异常编号- -类别、元素组合、主要类别、元素组合、主要元素、推断的主要矿种组颜灰、引线等。元素、推断的主要矿种组颜灰、引线等。 7 7 责任表（编图责任单位和责任等）。责任表（编图责任单位和责任等）。主要矿类的颜色主要矿类的线条异常编号-类别元素组合范围小的压盖范围大的异常甲地两次异常甲地两次异常效果对比效果对比异常方向性、与地质异常方向

25、性、与地质背景吻合程度大为改背景吻合程度大为改观（合理分割、取舍）观（合理分割、取舍）鄂西乙地综合异常效果图鄂西乙地综合异常效果图特点：反映组合异常范围、元素组合、组合异常元特点：反映组合异常范围、元素组合、组合异常元素规模相对大小，增加了综合异常范围。可考虑改素规模相对大小，增加了综合异常范围。可考虑改进：反映综合异常规模相对大小、反映主找矿矿种进：反映综合异常规模相对大小、反映主找矿矿种元素、反映综合异常找矿价值类别。元素、反映综合异常找矿价值类别。编制原则编制原则地球化学找矿预测图在综合异常分类、评序、评价的基地球化学找矿预测图在综合异常分类、评序、评价的基础上，结合地质矿产资料、异常查

26、证结果和找矿模型等础上，结合地质矿产资料、异常查证结果和找矿模型等评价要素，分矿种分别编制。要素：评价要素，分矿种分别编制。要素：图名图名地质矿产简图（素图）地质矿产简图（素图）综合异常图（素图）综合异常图（素图）找矿预测区边界、矿类、级别及编号找矿预测区边界、矿类、级别及编号 地质矿产图例地质矿产图例地球化学找矿预测区图例地球化学找矿预测区图例找矿预测区划分原则找矿预测区划分原则 地球化学找矿预测区主要依据地球化学异常分类评价地球化学找矿预测区主要依据地球化学异常分类评价结果、地球化学异常组合与空间分布（分带）规律等综结果、地球化学异常组合与空间分布（分带）规律等综合因素进行划分。合因素进行

27、划分。 以地球化学异常的平面分布（分区）、元素组合、成因分类为以地球化学异常的平面分布（分区）、元素组合、成因分类为依据，参考成矿区（带）、地球化学区（带）、地质构造区（带）依据，参考成矿区（带）、地球化学区（带）、地质构造区（带）的划分成果和地球化学推断地质构造成果，圈定找矿预测区。的划分成果和地球化学推断地质构造成果，圈定找矿预测区。 以同类异常的数量和找矿意义分类结果（或每个异常资源量定以同类异常的数量和找矿意义分类结果（或每个异常资源量定量预测结果）为依据对找矿预测区进行分级。量预测结果）为依据对找矿预测区进行分级。 以地球化学异常的平面分布（分区）、元素组合、成因分类为以地球化学异常

28、的平面分布（分区）、元素组合、成因分类为依据，参考成矿区（带）、地球化学区（带）、地质构造区（带）依据，参考成矿区（带）、地球化学区（带）、地质构造区（带）的划分成果和地球化学推断地质构造成果，圈定找矿预测区。的划分成果和地球化学推断地质构造成果，圈定找矿预测区。 以同类异常的数量和找矿意义分类结果（或每个异常资源量定以同类异常的数量和找矿意义分类结果（或每个异常资源量定量预测结果）为依据对找矿预测区进行分级。量预测结果）为依据对找矿预测区进行分级。找矿预测区全国统一划分找矿预测区全国统一划分A A、B B、C C三级三级矿预测区划分原则矿预测区划分原则 A A级：级：存在一个以上甲、乙类异常

29、，经与区内或附近已知矿床（田）找存在一个以上甲、乙类异常，经与区内或附近已知矿床（田）找矿模型对比分析，认为有希望找到（或新增储量）达大型以上规模的矿床矿模型对比分析，认为有希望找到（或新增储量）达大型以上规模的矿床或矿田的区域；或矿田的区域； 已知成矿类型有利，根据地球化学定量模型计算，预测资源量远已知成矿类型有利，根据地球化学定量模型计算，预测资源量远远超出已探明储量，有望找到（或新增储量）达大型以上规模的矿床（田）远超出已探明储量，有望找到（或新增储量）达大型以上规模的矿床（田）的区域；的区域； 异常特征和查证结果表明具有找到新矿种（接替资源）的巨大潜异常特征和查证结果表明具有找到新矿种

30、（接替资源）的巨大潜力，有望找到中型以上规模矿床的区域力，有望找到中型以上规模矿床的区域 B B级：级：根据成矿区（带）以外所建的找矿模型或理论模型判断，有一个以根据成矿区（带）以外所建的找矿模型或理论模型判断，有一个以上的乙类异常存在，有希望找到大型以上规模的矿床；根据地球化学定量上的乙类异常存在，有希望找到大型以上规模的矿床；根据地球化学定量模型计算，预测区预测总资源量巨大，有希望找到大型以上规模的矿床。模型计算，预测区预测总资源量巨大，有希望找到大型以上规模的矿床。 C C级：级：多个丙类异常存在，已知地质条件有利或一般，未进行异常查证或多个丙类异常存在，已知地质条件有利或一般，未进行异

31、常查证或查证后未获得重要突破，但推测有希望找到人业矿体或小型以上矿床；有查证后未获得重要突破，但推测有希望找到人业矿体或小型以上矿床；有甲、乙类异常存在，但人作或人程控制程度已经很高（包括深部控制），甲、乙类异常存在，但人作或人程控制程度已经很高（包括深部控制），深、边部找矿有一定潜力，但重大找矿突破的可能性较小。深、边部找矿有一定潜力，但重大找矿突破的可能性较小。找矿靶区划分原则找矿靶区划分原则相当于矿产成矿预测的相当于矿产成矿预测的、级预测区。从地球化学找矿预级预测区。从地球化学找矿预测区中划分出来的，与同成矿区（带）内的典型矿床（模型）测区中划分出来的，与同成矿区（带）内的典型矿床（模型）十分相似的，或通过三级查证发现有利找矿线索的，或与其十分相似的，或通过三级查证发现有利找矿线索的，或与其它预测方法高度吻合的，具有明确找矿方向和目标的甲、乙它预测方法高度吻合的，具有明确找矿方向和目标的甲、乙类异常分布区。类异常分布区。 一般根据一般根据1515万比例尺化探人作资料圈定。万比例尺化探人作资料圈定。 地球化学找矿靶区一般为可直接部署二级查证或普查以上地球化学找矿靶区一般为可直接部署二级查证或普查以上找矿评价人作的区域。找矿评价人作的区域。找矿预测图编图方法找矿预测图