地球化学化探数据处理与成图剖析
- 格式:ppt
- 大小:4.60 MB
- 文档页数:38
最新最详细化探数据处理与编图流程⼀、指导思想成矿地质背景地球化学研究就是从地球化学特征出发,借助已建⽴的地球化学信息提取技术,充分利⽤地球化学调查所获得的海量数据信息,提取有关反应成矿地质背景条件的地球化学信息,并编制相应地球化学图及相应的推断解释图件,为资源潜⼒评价有关成矿地质背景的研究提供地球化学⽀撑。
⼆、⼯作内容(⼀)基础图件成矿地质背景条件的地球化学信息提取⾸先是要编制有关基础地球化学图件。
主要有: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. 传统方法,均值加标准差 (1)2. 直方图解法 (2)3. 概率格纸图解法. (3)4. 多重分形法。
(4)5. 85%累计频率法。
(5)小结 (6)1.传统方法,均值加标准差在excel中用过函数,求均值,求标准差,先对数据中的极大/极小值进行剔除,大于/小于三倍标准差的剔除掉,直到无剔除点。
然后用均值加2倍标准差求异常下限。
图,D列中的函数,E列中的结果。
异常下限位114.86.。
2.直方图解法图1图2首先,做频率直方图,(图1的数据是某化探区数据)含量频率分布图上呈现双峰曲线,左边是背景部分,右边是异常部分,双峰间谷底处(0.7)为异常下限。
求真值得 5.所以,异常下限位5。
图2另一个化探区的数据,是单峰曲线,在频率极大值的0.6倍处画一条平行直线,与曲线一侧相交,其横坐标长度即为σ。
用Ca=Co+2*σ=0.16+2*0.665=1.49,求得为真值为31。
3.概率格纸图解法.图4是概率格纸。
发现纵坐标(累计频率)是不均匀的。
把样本值小于或等于某个样本n i的数据频率累加,即得到小于或等于n i的累积频率。
概率格纸用excel能轻松的做出来。
制造方法如下。
图4.图4显示了概率格纸的制造过程。
原理就是把标准正态分布曲线投影到纵坐标上。
首先确定纵坐标数值,如B列,0.1、1、5、10、20、30、40、50、60、70、80、90、95、99、99.9.。
如果想要纵坐标线密一点,也可以插入更多的数。
然后在C列中用NORMSINV函数,求对应频率的分位数(如果把标准正态分布,正着放,分位数就是横坐标)。
这时的原点(0)在50%处,我们想要原点在0处,那么把C列的数统一加-03.090232(C5),---(处理化探数据的时候,加的也是相同的数)。
即输入公式”D5”=C5-$C$5…。
E列为x 值,根据实际化探数据,设定最大和最小值。