化探数据处理原理及方法
- 格式:pdf
- 大小:2.39 MB
- 文档页数:59
化探数据处理的一般性方法一、分析质量: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 化探数据质量评价的数据处理(分矿区)⑴统计重采样和重分析抽查样所占样品总数的比例比例 = (重采样和重分析抽查样数/工作样总数)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坐标各元素的含量⑵整个矿区和各地质单元(各地层、各岩浆岩)样品各元素含量特征统计统计的参数包括:①元素含量平均值;②最大值;③最小值;④标准离差;⑤变化系数(标准离差/含量平均值);⑥浓度克拉克值(元素含量平均值/该元素的克拉克值)整个矿区和各地质单元统计结果含量平均值、最小值、最大值用表表示。
1数据准备(1)成图前逐点采集样点的平面坐标。
数据排列顺序为x坐标(米)、y坐标(米)、元素含量或采样点线号,图面坐标与实际坐标要进行转换。
分别建立各元素含量数据及点位点线号的电子表格文件,在Excel电子表格中另存文件,文件类型为.csv,逗号分隔(*.csv)。
(2)用记事本软件将*.csv文件打开,在首行加入notgrid,将文件类型"*.csv"改为"*.det"保存。
2取样点位图成图(1)进入MapGIS6.5主菜单,点击空间分析—DTM分析。
(2)进入文件菜单打开三角剖分文件…,即上述的*.det文件。
(3)进入模型应用菜单下选择高程点标注制图,弹出下列菜单。
(4)高程点标注显示框中的“符号”是指采样点位置符号,可以点击缺省符号区域自定义为“o”符号;标注是指元素点、线号。
(5)在标注位置框中设置数据的“对齐方式”选择“右边”,“角度”为“0”。
也可以自定义设置x、y位移量。
(6)在符号尺寸框中一般选择固定尺寸,为0.5~1。
(7)在“标注字体…”框中设置字体大小、颜色、类型。
(8)在“标注格式…”框中设置小数位数。
(9)设置完成点击“确认”。
(10)保存点文件、线文件。
3单元素异常图绘制(1)打开单元素高程数据三角剖分文件…,即单元素*.det文件。
(2)在TIM模型菜单下,进行快速生成三角剖分网,对于弹出的“*.tin”数据已修改是否先保存,点击否。
利用删除三角剖分网删除不必要的三角网。
(3)在TIN模型菜单下追踪剖分等值线。
(4)在设置等值线参数对话框中以前面计算的各元素异常下限为等值线外层,选择等值线套区、保留边界线、等值线光滑处理—高程度、轴向标尺。
(5)在等值线定层中设置所确定的等值线间隔、线参数类型、注记参数类型及修改区参数。
(6)在所有参数设置完毕,则点击确定,异常等值线自动生成。
(7)保存所生成的三角剖分文件、点文件、线文件、区文件。
化探数据处理及图件编制第二节分析方法及质量评述一、分析方法本次扫面和异常查证的全部样品均交由四川省地矿局华阳地矿检测中心测试,根据任务书要求共分析测试元素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个样品密码插入一组,与样品同时分析。
最新最详细化探数据处理与编图流程⼀、指导思想成矿地质背景地球化学研究就是从地球化学特征出发,借助已建⽴的地球化学信息提取技术,充分利⽤地球化学调查所获得的海量数据信息,提取有关反应成矿地质背景条件的地球化学信息,并编制相应地球化学图及相应的推断解释图件,为资源潜⼒评价有关成矿地质背景的研究提供地球化学⽀撑。
⼆、⼯作内容(⼀)基础图件成矿地质背景条件的地球化学信息提取⾸先是要编制有关基础地球化学图件。
主要有: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)进⾏了⽐对。
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、采样点坐标生成
如果为GPS测定,直接把坐标转成文本文档。
在投影变换中用户文件投影转换把坐标数据投成点文件。
(具体操作可问孔德成,注意:该步骤只有在带虚拟狗的MAPGIS61或67下运行)
如果为手工图件,则先把采样点做成点文件(注意:点的顺序要与化验结果一致),参照以下步骤操作:
(1)、在投影变换中打开点文件,在工具中选择编辑属性结构,按下图编辑后,点OK。
(2)在工具中选择点位置转换为属性,弹出下图,依次点击转换、保存、确定。
确定。
通过以上过程,在桌面上生成坐标文本。
2、异常下限的确定:从化验结果中统计下图中的各个值。
(用的是EXCEL中的公式,不明之处可问吴邦友),下图中异常下限之下的各值仅为判别推断用,可以不要。
3、异常图制作:先把点坐标及各元素化验结果复制成文本文档。
在DTM分析中选择GRD模型——离散数据网格化,打开文本文档,
参照下图设置:
确定后自动生成Tmpgrid.Grd文件。
选择GRD模型——平面等值线图绘制,打开Tmpgrid.Grd文件,根据异常下限设置等值层值及其它项,参照下图:
确定后生成等值线图。