当前位置:文档之家› 化探数据处理步骤

化探数据处理步骤

化探数据处理步骤
化探数据处理步骤

所有数据均输入计算机、以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为异常面积)

各类系数计算和所利用的公式均符合规范要求。

化探数据处理成图过程

化探数据处理成图的过程 毕武12段新力12黄显义12袁小龙12彭仲秋12李永华12 1?乌鲁木齐金维图文信息科技有限公司,新疆,乌鲁木齐, 830091 2?新疆地矿局物化探大队计算中心,新疆,昌吉, 831100 0前言 GeolPAS 软件用户群不断扩大,由于各用户对系统的熟悉程度不同,对软件 功能 了解不够,有必要分专题将GeolPAS 处理数据及成图过程做一系统总结,下 面就化探数据处理成图的过程做一总结。 数据输入 *数据预处理 I I I II 「图件绘「输出… 图1金维地学信息处理研究应用系统 (GeolPAS )图件制作流程 1处理步骤 化探处理的成果包括:(1)参数统计表;(2)R 型聚类分析-谱系图;(3)重复 样三层套合方差分析或者重复样合格率计算结果;(4)点位数据图;(5)地球化学 图;(6)直方图;(7)组合异常图;(8)综合异常图;(9)远景区划图;(10)单元素 异常参数统计(附表册);(11)异常剖析(附图册);(12)综合异常登记卡(附表册)。 在GeolPAS 系统中,化探数据处理分为以下几个主要步骤: 1.1数据检查 数值检查,坐标检查,重复样坐标检查。 1.2分析处理 重复样三层套合方差分析、重复样合格率计算、化探特征参数统计、化探背 景值分析、R 型聚类分析、因子分析。 1.3数据分析 原测数据 重磁数据处理 电法数据处理 化探数据处理 数学地质数据处理 水文地质数据处理

数据变换;衬值、累加衬值;数据累加、累乘、比值;异常归一化。

1.4网格化 离散数据网格化、XYZ数据转网格数据 1.5成图 点位数据图、彩色等量线图、直方图、组合异常图、单元素异常图、综合异常图、剖析图。 1.6单元素异常参数统计 1.7综合异常登记卡 图2化探数据处理成图流程 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 Ti 2 矿区地球化学特征研究的数据处理(以哈拉街吐为例) ⑴作出SSPS数据文件 作出下列SSPS数据文件: ①文件1:整个矿区数据文件; ②文件2:矿区地层数据文件; ③文件3:矿区岩浆岩数据文件; ④文件4 :下二叠统大石寨组(P1d)数据文件; ⑤文件5 :下白垩统大磨拐河含煤组(K1d)数据文件;

⑥文件6 :华力西晚期侵入岩数据文件; ⑦文件7 :燕山期早期侵入岩数据文件; ⑧文件8 :燕山期晚期侵入岩数据文件; ⑨文件9:已知矿附近一定范围数据文件 每一数据文件的内容项目包括: 序号野外号 X坐标 Y坐标各元素的含量 ⑵整个矿区和各地质单元(各地层、各岩浆岩)样品各元素含量特征统计 统计的参数包括: ①元素含量平均值; ②最大值; ③最小值; ④标准离差; ⑤变化系数(标准离差/含量平均值); ⑥浓度克拉克值(元素含量平均值/该元素的克拉克值) 整个矿区和各地质单元统计结果含量平均值、最小值、最大值用表表示。 ⑶整个矿区和各地质单元样品各元素的概率分布特征统计 ①标准离差 ②峰度 ③偏度 ④概率分布曲线特征 ⑷矿区各地层样品各元素的局域丰度和蚀变-矿化叠加系数特征统计 根据地球化学过程的基本定律(A.B.Vstelius,1960),一个矿区地层中元素的“丰度”应该是沉积岩沉积成岩时的初始平均含量,而不应包括后期岩浆、蚀变、矿化作用等地质作用造成的元素含量的增赢或亏损。而矿区内局部地区地层中元素的“局域丰度”,至少应排除最后蚀变-成矿作用叠加的那一部分元素的含量。若本区各地层中元素概率分布及其偏度和峰度特征表明元素呈偏对数正态分布。这说明地层中多数元素都受到了后期不同程度的蚀变-成矿作用的叠加。据此,剔除了不服从正态分布的超差样品(即含量大

化探数据处理成图过程.

化探数据处理成图的过程 毕武1、2段新力1、2黄显义1、2袁小龙1、2彭仲秋1、2李永华1、2 1.乌鲁木齐金维图文信息科技有限公司,新疆,乌鲁木齐,830091 2.新疆地矿局物化探大队计算中心,新疆,昌吉,831100 0 前言 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、>等字符,如果有这样的情况要找实验室给予纠正。数据准备好后,我们要把数据转换成TXT后缀的文本文件,这就做好了处理前的准备数据工作。 2.2 分析处理 2.2.1 重复样三层套合方差分析 一般是从分析样品的结果中挑出来重复样的分析值,每组四个样品,按如下顺序排列: 11 第一次采样第一次分析 12 第一次采样第二次分析 21 第二次采样第一次分析 22 第二次采样第二次分析 结果: 三层套合方差分析成果--元素:Cu 三层套合方差分析成果--元素:Au 查表:Fa0.05(18,19)=1.35 Fb0.05(19,38)=1.22 Cu:11.152 > 1.35 1.482 > 1.22 即:F1>Fa F2>Fb Au:10.781 > 1.35 0.483 < 1.22 即:F1>Fa F2

化探数据处理及图件编制

化探数据处理及图件编制 第二节分析方法及质量评述 一、分析方法 本次扫面和异常查证的全部样品均交由四川省地矿局华阳地矿检测中心测试,根据任务书要求共分析测试元素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个样品密码插入一组,与样品同时分析。然后计算四个监控样实际测定值与监控样定值之间的平均对数偏差X和对数标准差λ(统计结果见表3,2),其计算公式如下: L n ,lgC,i,1 ΔlgC,lg C-lgCX,L定值测定值 n n22(,lgC)-nX,L,1iλ, n,1 6、重复分析样编码 重复采样及重复分析样的编码,两个二级监控样的密码空号均由野外采样单位确定。同时野外采样单位负责重复采样和重复分析结果的方差分析。 三、质量评述 1、报出率

Mapgis处理化探数据程序与步骤

Mapgis处理化探数据程序与步骤 一、安装系统 (一)软件拷贝 将Mapgis65版以及模拟狗拷贝到计算机某个硬盘中 (二)装模拟狗 在控制面板界面中通过添加新硬件的方式装好模拟狗,注意不同的示窗操作系统用不同的版本。 (三)装Mapgis65 系统 打开硬盘中的Mapgis65文件夹,点击Setup按钮,计算机会自动安装,安装过程根据提示进行即可。 二、数据预处理 (一)数据分类 将各元素含量与其对应采样点的坐标分别用Excel表格程序整理好,并按元素分别建立不同的文件夹。 (二)数据格式转换 1、将Excel文本格式另存为纯文本格式 2、在Word界面打开某某元素的纯文本数据文件,在编辑菜单下点替换按钮,后将所有数据间的空格在英文状态下用逗号替换。然后将第一行(即元素符号及x、y坐标符号行)改写为“notgrid”,再将数据另存为det文件。 三、数据处理与成图 (一)高程点标注

1、打开Mapgis 主菜单,选择空间分析下拉菜单的DTM分析。 2、在“文件”下拉菜单中点击“打开三角剖分文件”,后打开某元素数据的det文件。 3、在“模型应用”菜单中选择“高程点标注制图”。 4、在高程点标注显示界面中,分别选择:标注符号及其尺寸大小,标注数据位置、字体大小、颜色、小数位数等。确认后,在“文件”菜单中选择“存数据于点数据文件”保存处理好的高程点数据。 (二)离散数据网格化 1、在空间分析下拉菜单的DTM窗口中,在“文件”菜单中选择“打开数据文件”打开前面保存的某元素高程点数据文件。 2、在Grd模型菜单中,选择“离散数据网格化”按钮。 3、在离散数据网格化界面中选择:X、Y坐标的起点、终点、网格间距、网格化方法(一般用距离幂函数反比加权网格化),修改文件换名,后确定。计算机将自动进行数据网格化处理。 4、保存网格化后的数据。 (三)无效化处理 1、在DTM窗口中打开已经进行的网格化Grd数据文件。 2、在“文件”菜单中选择“打开数据文件”,打开对应元素的点数据文件。 3、在在Grd模型菜单中,选择“规则网无效化处理”按钮。 4、在裁剪区域选区方式中选择“鼠标选取任意区域”,确认后,用鼠标选取需要勾绘等值线的区域,注意:在最后一个选择点封闭时,通

MAPGIS数字高程模拟系统在化探数据处理中的应用

MAPG IS 数字高程模拟系统在化探数据处理中的应用 何明华 (甘肃省地勘局第一地质矿产勘查院,甘肃天水 741020) 摘要:地球化学普查水系沉积物测量和土壤测量数据具有三维特征。利用MAPGIS 数字高程模拟系统绘制普查区金、银、铜等元素的点位图、原始数据图、地球化学图,不仅缩短了数据处理的周期,而且保证了空间数据的精确性,具有较强的实用性。关键词:数据高程系统;处理;化探;数据 中图分类号:P 208 文献标识码:B 文章编号:1001-358X(2005)02-0013-03 随着计算机数据处理能力的不断提高,自动测量仪器广泛使用以及制图技术的发展,一种全新的数字描述地理现象的方法)))数字高程模型(DT M )日渐普及。数字高程模型(DT M )是利用已有的观测数据经过专业处理产生,利用计算机自动产生各类专业地学图件并进行专业分析。MAPGIS 地理信息系统是集数字制图、数据库管理及空间分析为一体的空间信息系统。随着该系统在地质勘查部门的进一步推广和应用,地质找矿和地质勘查技术方法和手段有了进一步的提高。在地球化学普查中,利用MAPG I S 数字地面高程模拟系统绘制普查区金、银、铜等金属元素的点位图、原始数据图、地球化学图,不仅提高了地质工作人员的工作效率,缩短了工作时间,而且保证了空间数据的精确性。1 化探数据的特征 地球化学普查的主要目的是在勘查区域内查明成矿有利地段以及与找矿有关的地球化学特征,圈出 各类区域性异常及成矿有利的远景区,寻找目标靶区,为进一步开展其它地质勘查工作提供依据。通常采用的化探方法主要有:水系沉积物测量和土壤测 量,工作精度为1:5万,采样密度一般为4-8个点/km 2 。采样点在勘查区分布比较均匀(图1) [1] 。 图1 设计采样点分布图 水系沉积物测量和土壤测量数据具有三维特征,即采样点平面坐标及每一个采样点样品经过化验得出的金、银、铜、铅、锌等金属元素的含量,平面坐标由地质人员利用GPS 全球定位系统采集。数据结构见(表1) 表1 水系沉积物测量数据一览表 采样点号 样品编号 采样点坐标Au Ag Cu Pb X Y w(10-6)w(10-6)w(10-2)w(10-2)88c 13821650440470117010812511261588d 23821080440530114010342114241289c 33821300440721110010772317221289d 43822078440789117010872415231490c 53821850440920313010492211311290d 63822230440500018011352512221589a 73822800440850110010662219211888b 83822330441240113010962517221680c 938219504412502011012582513261789b 103822870441280119010902414321890a 11 3822650 441680 113 01158 2516 3513 13 第2期2005年6月 矿 山 测 量M I NE SURVEYI NG No 12 June 12005

地球化学数据处理与图件编制方法流程

地球化学数据 处理与图件编制方法流程 编写人:刘红杰 QQ:498236930 内蒙古第三地质矿产勘查开发院

一、指导思想 成矿地质背景地球化学研究就是从地球化学特征出发,借助已建立的地球化学信息提取技术,充分利用地球化学调查所获得的海量数据信息,提取有关反应成矿地质背景条件的地球化学信息,并编制相应地球化学图及相应的推断解释图件,为资源潜力评价有关成矿地质背景的研究提供地球化学支撑。 二、工作内容 (一)基础图件 成矿地质背景条件的地球化学信息提取首先是要编制有关基础地球化学图件。主要有: 1. 39种元素(化合物)地球化学图 2. 地球化学组合异常图 3. 地球化学综合异常图 (二)解释推断图件 地球化学解释推断图件,内容包括: 1. 地球化学推断解译地质图 2. 地球化学找矿预测图 三、工作方法 (一)数据校正处理 由于区域地球化学数据受地理景观、采样介质、分析手段的影响,不可避免的产生明显的系统误差,尤其是涉及到区域性的化探数据,这种误差更为突出。因此,在各省进行数据处理与专题地球化学图编制之前,有必要分别对各元素进行系统误差的处理,以便能更好地反映地质现象和矿产信息。误差处理主要针对图幅间(包括分析批次)明显的系统分析误差(必须处理)和地质景观环境差异影响解释的效果(根据解释的需要确定)。

1. 系统误差特征及处理原则 (1)分析误差源,所展示的数据误差与周边数据值具有明显的台阶状。 (2)数据误差在空间上具有区域性特点,区域、图幅或分析批次。 (3)在数据值的分布上,掩盖了地球化学特征和地质特征展布的延续性和规律性。 (4)在数据处理方法上,尽可能地选择线性校正,通过简单的计算可以复原数据。 2. 系统误差处理步骤 (1)按原始点位采用符号分级的方式生成元素的符号图,分级方法采用累计频率方式。 (2)通过校正图示窗浏览原始数据全图,确定具有明显的数据台阶区域,区域的确定原则是由区域->图幅->批次;采用图形编辑工具,在图上直接圈定要处理的区域(用面的方式表示)。 (3)确定局部图幅和分析批次范围产生的系统误差,校正单元由系统提供的工具直接在显示窗中勾绘,确定完所有需要校正单元,各校正单元的ID需设定为唯一。 (4)建立校正单元与处理数据表空间位置索引关系。 (5)确定各单元的校正值或校正系数,主要方法是与单元周边数据进行对比分析,部分规律性较复杂的单元可以通过统计规律确定,同时还需考虑地球化学分布的整体空间分布趋势和地质背景; 计算方法推荐采用:V ai=A V i+B 其中:V ai为校正点校正后数据,A为校正系数,V i为校正点原始数据,B 校正常数。A与B值的确定参照校正单元周边数据单元(正常的数据单元)确定,可以由经验判断,也可通过统计规律确定。 (6)数据校正,可采用SQL语言操作模式或应用软件系统提供的专用工具,按确定的校正值对各校正单元逐一进行计算。 元。 (5)单元校正计算,以基准单元为基础,应用归一化方法调整其他单元的数据值。计算方法推荐采用: V aij=V ij×(V j/V r)

最新最详细化探数据处理与编图流程

一、指导思想 成矿地质背景地球化学研究就是从地球化学特征出发,借助已建立的地球化学信息提取技术,充分利用地球化学调查所获得的海量数据信息,提取有关反应成矿地质背景条件的地球化学信息,并编制相应地球化学图及相应的推断解释图件,为资源潜力评价有关成矿地质背景的研究提供地球化学支撑。 二、工作内容 (一)基础图件 成矿地质背景条件的地球化学信息提取首先是要编制有关基础地球化学图件。主要有: 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)进行了比对。计算两次分析值之间的相对偏差(RE%),具体计算采用如下公式:RE%= |(C1-C2)︳/(C1+C2)×100%,当RE<33.3% 时为合格,合格率=合格样品/总样品数×100%。总合格率大于80%。 推荐软件:GeoExpl,MapGIS, Geoipas1.64 .Suffer (二)坐标投影变换 在坐标投影变换和成图时经常出现的是将“源数据投影参数”的单位、比例尺弄错的情况,比如把数据直接转换成结果投影的单位等,这些是不需要做的,我们一般工作默认用的投影参数就是我们的地图参数,比如“投影平面直角坐标,北京54,高斯-克吕格投影坐标系”或者“投影平面直角坐标,西安80,高斯-克吕格投影坐标系”,实际工作的坐标单位一般用米,比如我们要成5万图,那参数设置就是: 源数据投影参数,比例尺:1,坐标单位:米,21度带 结果投影参数,比例尺:50000,坐标单位:毫米,21度带 推荐软件:GeoExpl,MapGIS ,Geoipas1.64 (三)数据网格化 离散数据网格化处理是空间数据插值的一种,即把无规则分布的空间数据内插为规则分布的空间数据集。网格数据是编制地球化学图件的重要数据源。 网格化处理一般包括这样几个过程:①空间几何属性的确定;②插值方法(模型)的选择;③空间数据的探索分析,包括对数据的均值、方差、协方差、独立性和变异函数的估计等;④插值方法评价;⑤重新选择内插方法,直到合理。 网格化数据处理中要确定主要参数包括: 1. 网格距:根据采样密度确定,一般网格距应与采样密度一致; 2. 数据搜索半径:一般选择网格距的2.5倍; 3. 数据计算模式:最近点或距离指数加权;

相关主题
文本预览
相关文档 最新文档