化探数据处理
化探数据处理的一般性方法
化探数据处理的一般性方法一、分析质量：1．看技术报告中质量评述部分，看各项指标是否合格。
主要有：检出限、报出率、异常检查率、异检合格率、内检率、内检合格率、国家一级或二级标准物质的准确度和精密度等。
看是否符合标准（设计）。
这些数据由化验室提供。
这些一般不会超差的，否则化验室不能给出化验报告。
2．RE计算RE=abs（C1-C2）/（C1+C2）×200要求：小于3倍检出限时，RE≤85%，大于3倍检出限时，RE≤66%为合格。
总合格率一般要求大于70%。
需要说明的是：这种方法适用于简单对比，就是一个点取2个样时使用。
如果用三重套合分析，就不能用这种方法了。
所以写设计时，一定要用简单对比方法。
还有，这些年的化探中的RE还没有超标的，如果真的超标了，我也不知道怎样处理。
二、数据处理：1．剔除一级及二级标准样；2．剔除重复样；3．剔除0值及化验室输入错误的值，或小于检出限的值。
主要是指≤0.3这样数据的≤号；4．做原始数据图；5．计算异常下限，做单元素异常图，圈定单元素异常。
单点异常，只有外带的不圈，有中、内带的圈出。
外带用黄色，中带用浅红色，内带用深红色。
单元素异常编号为Au-1，Ag-1等。
需要指出的是，圈定异常时应该形成数据异常图，但交报告时，必须分开，就是形成一张数据图，再形成一张异常图；需要说明的是，如果面积较大（这个没有标准），总体说是水系面积超过一个5万图幅，就要分子区计算下限。
分子区的原则是不同年代、不同地质体都要划分成子区。
如果各个子区的异常下限接近，就采用总的，否则分别确定子区的异常下限，然后分别圈定子区的异常。
6．填单元素异常评序表。
异常点数、面积、平均值、极大值、标准离差、衬度、规模、浓度分带等。
评序有5参数和多参数两个评序，5参数不参与的参数有标准离差、浓度分带及异常点数。
需要说明的是，排序时，单项值高的给1，以下类推2、3等，一样的值给一样的排序。
不产生空的值。
浓度分带有内带的给1，中带的给2，外带的给3；7．做组合异常图，先把所有的元素做成1张组合异常图，只要外带，圈定组合异常。
地球化学化探数据处理与成图
（三）位置的相对性 无论是土壤测量还是水系沉积物测量所获得的异常， 往往与异常源都会发生不同程度的位移。这种位移与表生 介质本身的位移程度和采样的布局有关。特别是水系沉积 物异常的位移更为明显，可达几公里甚至更大的距离。因 此，查明异常与异常源的空间关系，就成为异常查证中的 首要任务。 （四）表生作用带来的复杂性 不同的景观条件下，表生地球化学作用会有很大的 差异，制约了元素在表生环境中的分散和富集。因此，只 有在同一景观内，异常才有较好的可对比性。地理景观不 同，表生地球化学作用也就不同，元素在地表迁移、分散、 富集的规律也就不同。在异常对比和解释上，除考虑引起 异常的原生因素（地质背景、矿床类型）以外，在一定程 度上必须注重异常所处的地理景观条件及表生地球化学环 境。一般来说，只有同一地理、地质景观区的区域化探异 常才有对比研究的基础。
2、化探数据处理解决的主要问题：
①研究采样和分析中的误差，优化采样布局
②抑制数据噪音，突出主体趋势
③揭示多种数据的内在联系，提取隐蔽的有用信息 ④显示数据空间分布模式，编制地球化学图件
⑤异常对比、分类、评序，等等。
3、化探数据处理中应该注意的问题： 1、地球化学数据通常蕴含多种有用信息并伴随某些 不规律的变化，同时在数据获取过程中还存在分析测定误 差，这些使化探数据的复杂性增加了，在化探数据处理中 要将这些不规律成分和分析误差除去。 2、找矿信息总是同地球化学异常相联系的。最普通 的化探数据处理是对一组化探数据计算出背景值和变化范 围（如用平均值和标准离差来衡量），据此确定出地球化 学异常的下限值。当地球化学背景随着地理位置出现趋势 变化时，要相应地采取适当的处理方法以便获得随地理位 置而变的背景值和异常下限。
化探讲座数据处理20120212
（4）在数据校正方法上，尽可能地选择线性校正，通过简
单的计算(jì suàn)可以复原数据。。
共八十页
二、常规(chángguī)数据预处理
▲ 元素分析系统误差校正校正步骤
(1) 按原始点位采用符号分级的方式生成元素的符号图或（累计频率）地球化
学初图。
(2) 确定具有明显的数据台阶区域，区域的确定原则是由区域->图幅->批次
低于检出限的数据取检出限的1/2, 如：<3.0=1.5
高于检出最高极限的数据取125%,如：>1000=1250
特异值（高出周边很多的但点异常值）
（1）实验室处理；（2）野外处理；（3）实际处理（背景(bèijǐng)剔除、统计及异常
保留）；4、空白值（单空点——多点平均，区空保留）
共八十页
三、数据(shùjù)统计分析
▲根据各省及示范区确定的地理投影模型，对准备各类空间数据（数据表和图形
数据）转换为统一的坐标系统(xìtǒng)，椭球参数建议采用西安80/IUUG 1975坐
标系(即长半轴：6378140米，短半轴：6356755.3041米)。
▲推荐软件系统(xìtǒng)：GeoExpl,MapGIS . MapGIS .Geoipas1.64 等
--2012年院培训(péixùn)讲座
化探资料
综合整理常规方法技术
(zīliào)
湖北省地质(dìzhì)调查院
共八十页
讲课(jiǎng kè)内容
一、化探数据来源及特点
二、常规(chángguī)数据预处理
三、常规数据统计分析
四、地球化学背景及异常确定与评价表述
五、地球化学数据处理新方法
化探数据处理成图过程
化探数据处理成图的过程毕武1、2段新力1、2黄显义1、2袁小龙1、2彭仲秋1、2李永华1、21.乌鲁木齐金维图文信息科技有限公司，新疆，乌鲁木齐，8300912.新疆地矿局物化探大队计算中心，新疆，昌吉，8311000 前言GeoIPAS软件用户群不断扩大，由于各用户对系统的熟悉程度不同，对软件功能了解不够，有必要分专题将GeoIPAS处理数据及成图过程做一系统总结，下面就化探数据处理成图的过程做一总结。
1 处理步骤化探处理的成果包括：(1)参数统计表；(2)R型聚类分析-谱系图；(3)重复样三层套合方差分析或者重复样合格率计算结果；(4)点位数据图；(5)地球化学图；(6)直方图；(7)组合异常图；(8)综合异常图；(9)远景区划图；(10)单元素异常参数统计(附表册)；(11)异常剖析(附图册)；(12)综合异常登记卡(附表册)。
在GeoIPAS系统中，化探数据处理分为以下几个主要步骤：1.1 数据检查数值检查，坐标检查，重复样坐标检查。
1.2 分析处理重复样三层套合方差分析、重复样合格率计算、化探特征参数统计、化探背景值分析、R型聚类分析、因子分析。
1.3 数据分析数据变换；衬值、累加衬值；数据累加、累乘、比值；异常归一化。
1.4 网格化离散数据网格化、XYZ数据转网格数据1.5 成图点位数据图、彩色等量线图、直方图、组合异常图、单元素异常图、综合异常图、剖析图。
1.6 单元素异常参数统计1.7 综合异常登记卡图2 化探数据处理成图流程2 具体处理过程2.1 数据检查我们的数据处理工作从化验室提供的样品分析报告开始，项目要提供坐标和样品对应的分析数据，坐标我们一般取实际工作中的米单位，系统中默认东西向横坐标为X坐标，不加带号，南北向纵坐标为Y坐标，需要时还要提供样品对应的地质编码，我们拿到这个数据后首先进行数据检查，以确保数据中不出现写错、漏填、负数、0、>等字符，如果有这样的情况要找实验室给予纠正。
化探数据处理原理及方法
(二)单元划分基本类型
1.规则单元
网格单元 2.自然单元 (1)地质体单元 定性划分;
矿体往往位于地质体内,接触带和地质体外.
(2)地质异常单元 以网格单元为样品单元,样品单元大小根据预测尺度而定; 综合致矿信息定量标度的样品异常单元集合.
（三）地球化学数据分析
1. 单变量 （1）针对要解决的地学问题，分析多源信息来源、原理、 精度、准确度 、用途及意义。
（四）主要成矿元素的统计意义
（四） 图示分析
（直方图、点阵图、三角图、玫瑰图、曲线图、投影图等） 1、直方图（单元素含量－频数直方图） 2、散点图（只管展示元素间、介质间相互关系）
3、饼图（含量分布百分比）
4、研究统计分布特征 （1）正态分布（分布比较均匀或样本密度大） （2）对数正态分布（元素含量变化大或多因素叠加） （3）二项分布（不确定性大，地质体复杂极不均匀）
i 1
n
n 1
式中:i=1，2，…,n；为样本数；j=1,2,…,m为变量 数； xij 为原始观测值； Sj 为标准偏差； xj 为平均值； 处理后的xij值yij为无量纲数据。
二、多源地学信息分析与变换
4、极差化
y ij
( xij x j min ) ( x j max x j min )
（一）地球化学找矿分类
地球化学找矿——地球化学探矿，化探（Geochemical Prospecting） 1、方法分类
测量方式： 航空化探——放射性、气体
海洋化探——海水、海底沉积物、生植物、放射性、气体 地表化探——岩石、土壤、水系沉积物、水、生植物、气
地下（井中）化探——岩矿石、水、放射性、气体
（二）背景与异常的概念
化探数据处理与解释评价资料课件
CATALOGUE
目 录
• 化探数据处理概述 • 化探数据处理方法 • 化探数据解释 • 化探评价资料 • 化探数据处理与解释的实践案例 • 总结与展望
01
CATALOGUE
化探数据处理概述
化探数据的特点
多元性
化探数据通常包含多种元素或化 合物的浓度信息，呈现出多元性
数据质量挑战
原始化探数据可能存在 采集、传输、存储等方 面的误差，影响处理与 解释的准确性。需要采 取合适的质量控制措施
以提高数据可靠性。
多源性数据融合
在实际工作中，化探数 据通常需要与其他地质 、地球物理、地球化学 等多源性数据进行融合 解释。如何实现多源性 数据的有效融合与协同 解释是一个重要挑战。
便于后续解释评价
将处理后的数据用于后续 的地质解释和资源评价， 提高工作效率。
化探数据处理的基本流程
1. 数据收集与整理
收集原始化探数据，并进行必要的格式转换和 整理。
01
3. 特征提取与选择
利用统计方法、图像处理等手段提取 与地质目标相关的特征，并选择重要
特征。
03
5. 数据可视化与表达
将处理后的数据通过图表、图像等方式进行 可视化表达，便于后续解释评价。
指导找矿方向
数据解释可以揭示地质构造和成矿规律，从而指 导找矿工作的方向。
数据解释的方法
统计分析法
通过对化探数据进行统计分析，可以了解元素的分布特征、异常 形态等，进而推断地质背景和成矿可能性。
地质解释法
结合地质资料，对化探数据进行解释，从地质角度揭示成矿规律和 矿产分布。
地球化学模型法
利用地球化学模型对化探数据进行处理和解释，可以更深入地了解 元素迁移、富集规律。
化探数据处理步骤
所有数据均输入计算机、以MAPGIS 制图系统为平台，以原始数据筛选替换特高值后，转计算成对数值后，用 X +2S 求出异常下限，分别以X －2S 、X －0.5S 、X +0.5S 、X +2S 分出色区，绘制各元素地球化学图，以上做图过程均在计算机上用MAPGIS 软件完成。
对化验室的样品分析结果取对数分组作直方图，证明所有元素均符合对数正态分布。
元素异常参数的确定：首先对原始数据进行假设正态检验，再作X －
+3S 特高值逐步剔除，然后进行各参数统计。
Au 元素含量为W×10－9，其它元素含量为W×10－6。
（1）背景平均值：f
fxc X ∑∑= （2）对数标准离差：1)(22-∑-∑=
n n fxc fxL S （3）对数异常下限：T0=X －L+2S
（4）变异系数：%1001%2230285.2⨯-=⋅S e Cv
e －自然对数，2.30285为常用对数与自然对数模数的倒数
S －对数标准离差，1为常数
（5）衬度：To
Xa K =（Xa 为异常平均值） （6）异常规模：P=k×km 2（km 2为异常面积）
各类系数计算和所利用的公式均符合规范要求。
化探数据处理说明
化探数据处理说明平顺项⽬化探数据处理⽅法⼀、基本概念1、异常⾯积：⽤GeoCIPS 软件直接统计异常⾯积，计量单位为km 2。
2、异常强度：异常区内⼤于下限数据的算术平均值。
3、最⼤值：异常区内数据最⼤值。
4、异常下限：根据作图结果调整异常下限。
5、异常衬度：异常均值/异常下限。
6、⾯⾦属量：元素剩余含量（异常均值减去异常下限）与异常⾯积的乘积。
7、NAP 值：异常衬度×异常⾯积。
8、异常排序：各异常按NAP 值⼤⼩排序，⼤的在前。
9、⾦计量单位⽤×10－9（ppb ），其余元素⽤×10－6（ppm ）。
⼆、单元素异常的圈定1、异常下限的确定表5-1 各元素异常下限⼀览表单位：Au 为ppb ，其它元素为ppm 。
逐步剔除法剔除⼤于+3S 的值、⼩于-3S 的值后求元素的平均值和标准差S ，选择+2S 定为计算下限，计算下限作为参考，根据表4-1分析结果对各分析指标的下限适当调整，对于有找矿可能的分析指标根据成图效果稍降低了下限，保留了较多的异常，对于找矿指⽰意义⼩的分析指标则提⾼了异常下限，仅保留了异常强度⾼的异常，将部分异常⾯积⼩、异常极⼤值/下限⼩、异常点数为1（少数为2）的异常删除，⼒求异常图可以直观的反x x x x映预查区的元素特征。
具体采⽤的异常下限及浓度分级见表5-1。
2、浓度分级预查区各分析指标尽量采⽤1、2、4分带。
由于预查区各分析指标整体含量低，仅Au、W采⽤了1、2、4分带，部分变异系数⼩和采⽤下限较⼤的分析指标，灵活调整了浓度分带。
三、综合异常的圈定与类别划分1、综合异常的圈定根据预查区内单元素异常分布及组合特征，以主要成矿元素的异常为主，把在空间上分布基本⼀致，相互重合的多个单元素异常圈定为⼀个综合异常。
共圈定以Ag、As、Au、Bi、Cu、Hg、Mo、Pb、Sb、Sn、W、Zn为主要异常元素的综合异常7个。
按所处地质环境、找矿意义和已有资料现阶段的认识⽔平，将各综合异常按下列标准进⾏分类，其中⼄2类异常1个，⼄3类异常4个，丁类异常2个。