化探数据处理方法

化探数据处理的一般性方法一、分析质量：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. 地球化学组合异常图3. 地球化学综合异常图（⼆）解释推断图件地球化学解释推断图件，内容包括：1. 地球化学推断解译地质图2. 地球化学找矿预测图三、⼯作⽅法（⼀）数据校正处理1|数据检查的必要性，因为实验室的分析报告还是⼿⼯输⼊的，还是存在录⼊错误的，我们重点检查的是“>”，数据中间的空格等录⼊错误问题；另外还有畸变检查，数据的特⼤值,⽐如超过10倍变差，⼀般对这样的分析值实验室会很重视的,你也可以提出让他们再确认⼀下,做到⼼中有数。

另⼀类错误可能会是我们录⼊样号或者坐标时出现的错误，如：“56b” 写成“56 b”,程序是以空格分开数据的，数据如果写成这样就会产⽣错误结果，有时在完成处理后才可能发现，这样⼀来我们前⾯的⼯作就作废了。

所以数据检查是⾮常必要的。

2|异常下限值的确定采⽤逐渐剔除法：①计算全区各元素原始数据的均值(X)和标准偏差(S)；②按X1+3S1的条件剔除⼀批⾼值后获得⼀个新数据集,再计算此数据集的均值(X2)和标准偏差(S2)；③重复第⼆步，直⾄⽆特⾼值点存在，求出最终数据集的均值(X)和标准偏差(S),则X做为背景平均值，S为标准离差，T(异常下限值)= X (背景平均值)+2S(标准离差)求出理论异常下限值，再结合地球化学等量线、地质背景及圈定效果确定出实⽤异常下限值。

3|重复样样品合格率统计野外重采样品以密码样形式插⼊样品中进⾏了分析，结果(C2)与第⼀次分析结果(C1)进⾏了⽐对。

所有数据均输入计算机、以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:1万化探的数据处理成图方法探讨毕武1、2段新力1、2袁小龙1、2黄显义1、2彭仲秋1、21.乌鲁木齐金维图文信息科技有限公司，新疆，乌鲁木齐，8300912.新疆地矿局物化探大队信息中心，新疆，昌吉，831100摘要：对1:1万化探的数据处理成图方法，在实际工作中的一些应用经验。

关键词：1:1万化探数据处理成图方法0前言随着1:5万区域化探的实施，在1:5万区域化探异常范围内开展的1:1万化探勘查也越来越普遍，对1:1万化探目前没有具体的规范，一般是参照1:5万区域化探规范执行，可是成图效果却有时不尽理想，下面就两个例子举例说明1:1万化探的数据处理成图方法。

1编制的步骤及要求1.1图件编制的要求：a土壤地球化学测量图件分两部分：基础图与推断解释图。

b区域调查和普查工作的图件编制按化探区域调查有关规范执行。

c图件编制必须符合地球化学勘查图式图例及用色标准规定。

d成果报告需做交通测区位置图、实际材料图、等值线图、综合异常图及其它推断解释图。

1.2例子一：南北向测线，100×40网格距。

这样的数据本身符合网格数据的格式，只是以TXT文本格式保存，所以我们可以借助GeoIPAS软件的“数据预处理”→“数据格式转换”→“XYZ数据转为网格数据”功能，这里注意数据坐标的起始值、终止值和网格间距，网格间距就选择100×40，起始值和终止值要保证数据的原始点位都与网格点重合。

这里有一个问题，就是有些项目用的是GPS坐标，不一定是与网格点位完全重合，一般我们建议用规则网的坐标。

下面先看一下图1-1点位数据图和图1-2点位符号图；图1-3原始数据转为网格数据绘制的地球化学图以及其转换参数：X最小值320000，X最大值32500网格化0，间距100，列数51；Y最小值4121000，Y最大值4127000，间距40，行数151。

图1-4距离平方导数加权网格化绘制的地球化学图，GeoIPAS系统→数据预处理→离散数据网格化→距离平方导数；搜索范围：圆域，R=300米；网格化参数：X最小值319880，X最大值325120，间距40，列数132；Y最小值41209200，Y最大值4127080，间距40，行数155。

化探数据处理说明平顺项⽬化探数据处理⽅法⼀、基本概念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坐标各元素的含量⑵整个矿区和各地质单元（各地层、各岩浆岩）样品各元素含量特征统计统计的参数包括：①元素含量平均值；②最大值；③最小值；④标准离差；⑤变化系数（标准离差/含量平均值）；⑥浓度克拉克值（元素含量平均值/该元素的克拉克值）整个矿区和各地质单元统计结果含量平均值、最小值、最大值用表表示。