化探数据处理

化探数据处理的一般性方法一、分析质量：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张组合异常图，只要外带，圈定组合异常。

化探数据处理成图的过程毕武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、>等字符，如果有这样的情况要找实验室给予纠正。

所有数据均输入计算机、以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个。

内蒙古扎赉特旗东芒合矿和哈拉街吐矿化探数据处理及图件编制方法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坐标各元素的含量⑵整个矿区和各地质单元（各地层、各岩浆岩）样品各元素含量特征统计统计的参数包括：①元素含量平均值；②最大值；③最小值；④标准离差；⑤变化系数（标准离差/含量平均值）；⑥浓度克拉克值（元素含量平均值/该元素的克拉克值）整个矿区和各地质单元统计结果含量平均值、最小值、最大值用表表示。

化探数据处理及图件编制第二节分析方法及质量评述一、分析方法本次扫面和异常查证的全部样品均交由四川省地矿局华阳地矿检测中心测试，根据任务书要求共分析测试元素14种。

样品从加工到测试到质量监控均按中华人民共和国地质矿产行业标准DZ0130-1994《地质矿产实验室测试质量管理规范》、ISSN-1870《1?5万区域地质调查及地球化学样品分析方法及质量管理指导性规程》和2002年新疆地勘局试验管理科《1?5万化探样品分析质量过程管理规则报告》进行。

14种元素的分析方法见表3,3。

二、技术要求1、报出率十四种元素的总报出率应大于95%。

2、外检样对已测试样品，测试单位按照3%的比率进行外检。

3、分析质量检查及质量监控方案为了有重点地监控元素的分析质量，实验室在送样单位确定的分析元素中，要再选择若干种主要监控元素并根据这些元素在本省制备的全部GRS二级标样中选择四个在元素含量范围及基体组成均为合适的GRD二级标样作为本图幅质量检查监控之用。

主要监控元素和二级标样的选择均应和送样单位协商进行。

每一大批样品测定完毕后，应将数据交给质量管理人员，对每一小批中插入的四个二级标样及四个重复分析(内部检查)样进行统计计算，并及时绘制日常质量监控图，在日常金的分析工作中，必须进行不小于10%的内检抽查。

为满足在一个较大范围的成矿远景区带内的1?5万图幅的拼接，应对分析的准确度进行检查和考核，为此实验室应在每一个1?5万普查化探项目完成后，分析8个GSD一级标样一次，痕金分析也应用金标样作准确度检查。

准确度和精密度计算结果应符合表3,2的要求。

4、微量金由于金元素在自然界中的均匀度和赋存状态对分析检测影响比较大，为确保金元素的分析质量，化验室特采用两种监控措施:第一，在每一分析批次的50个样品中插入两个国家?级标准物质GBW系列，用以计算实测值与推荐值之间的对数偏差:ΔlgC,lgC,lgC; 定值实测值第二，该地区样品分析结果结束或阶段性结束后，再对高、低异常点进行随机抽样检查约20%.5、?级标样为严格监控各元素的分析质量，实验室选取了四个不同含量的GRD系列监控样，每批次50个样品密码插入一组，与样品同时分析。