地球化学化探数据处理与成图
用Excel与Mapgis处理化探数据与直方图编程
用Excel及Mapgis处理化探数据1、软件要求：Microsoft Office Excel 2003及MAPGIS672、参考资料：如何利用Excel处理化探数据，应用MAPGIS制作地球化学图单元素异常图及综合异常图3、Excel准备工作：①在化探找矿工作中 ,大量数据的处理是一项重要的前期基础工作 ,数据处理的及时、准确与否 ,将直接涉及到下一步勘查工作的部署。
作者在本文中以 Microsoft Office XP Excel工具～数据分析功能为工具 ,简述了这些功能在化探数据处理当中的一些应用。
缺省的 Windows Professional XP一般不安装用来数据处理所需要的数理统计功能 ,故需重新加载 ,步骤如下:工具 -加载宏 -分析工具库 -确定。
之后 ,还需检查在工具菜单下有无数据分析选项。
在大量数据的录入过程中 ,可以设定工具 -语音选项。
在某个数据输入完毕之后 ,按 Enter键即可语音朗读 (需要音响或耳麦 ) ,这样可以实现数据录入的同步检查 ,确保录入数据的准确性。
②数据的处理正态性检验数据是否服从正态分布(或对数正态分布)是化探数据处理的一个重要前提条件，因此必须先进行数据正态性检验。
4、Excel操作步骤：1）、Excel2003准备：把测试数据及样号坐标整理为下图所示2）、打开：①工具——数据分析——描述统计（如下图）②再点确定如下图（Au的数据统计）③再点确定如下图（Au的数据统计）这样就可以得出一系列所需数据，但经验算该组数据为非正态分布，不能做进一步计算，但此为计算流程。
5、MAPGIS操作步骤：1）、首先将Excel数据复制粘贴为记事本格式（.TXT格式）2）打开Mapgis如下①、选择②、Grd模型——离散数据网格化③、打开（做Au图）选择数据列D（X、Y、Z）—网格化方法M（按不同数据选择不同网格化）—输出网格文件名F（存储数据的位置）—确定。
④、选择打开上面保存的⑤、打开⑥、确定（如下图）⑦、选择文件—数据存储于—点数据文件、线数据文件、区数据文件（以上点线区文件保存完毕，以下就是用Mapgis正常编辑）⑧修改要编辑的项目化探直方图制图软件的编程方法2008-12-21 23:36软件采用VB+MAPGIS SDK 开发数据载入格式： CSV文本文件，内容为包含分区编码的元素分析数据分区编码与元素代号均采用MAPGIS的书写方式如C1表示为C#-1，以便成图时正确显示。
化探数据处理的一般性方法
化探数据处理的一般性方法一、分析质量：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）高程点标注显示框中的“符号”是指采样点位置符号，可以点击缺省符号区域自定义为“ｏ”符号；标注是指元素点、线号。
（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段新力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、>等字符，如果有这样的情况要找实验室给予纠正。
地质数据处理常用方法介绍
一、原始数据准备与转换
•数据分布特征：离散点 或 网格化数据 •文件格式：图幅库中加入的点图层文件（*） •实验数据：图幅
一、原始数据准备与转换
•文件转换成所需的文本文件 •将文本投影成 点文件
二、数字特征计算
.数字特征：通过最大值、最小值、方差等基本统计量以及 分级统计频数，以专题图件（直方图、饼图）形式表现 .方式：统计整个图层的数字特征、按照某个指定区域的数 字特征、多个不同子区的数字特征 .专题图件的编辑
等值线方法：方法、方法、方法 可对等值线极值点参数设置 可根据“标注格式参数”，选择输出标注：使用
“固定小数位数”还是“科学计数法”；小数部分 是几位 若计算数字特征时选择了对数计算方式（数据已对 数），可还原显示 具有“等值线套区”、“绘制色阶”、“保留边界 线”；等多参数控制绘图
三、地球化学图制作
数字特征 等值线方法选择，及其相关参数的设置（搜索方法、网格数 等） 极值点参数设置：子图、标注字体、标注格式 等值线追踪
注意点：数字特征 － 对数计算方式时，是否已经选对数
自动计算数值特征：三种方式
：绘制的等值线经过点 ：绘制的等值线可能不经过点
注意：若计算数字特征时选 择了对数计算方式（数据已 对数），则需要进行还原显 示
判别分析：
判别分析的前提是已经存在几种已知的分组类型， 需要对未知样品判别其归属于哪一类。
神经网络:
可用于离散数据变量间的非线性关系拟合（对多， 多对，多对多），达到预测的作用；用于离散数 据的模式识别，即先建立样本属性与模式之间的 关系，通过回判达到识别样本模式的功能，属于 有监督的模式识别。
选择刻划对象间两两接近程度的要素和具体标定方法，是聚类分析的关键性基础 工作。它主要依靠地理专业知识，具体标定方法很多，常用的有相关系数法、欧 式距离法、海明距离法、夹角余弦法及数量积法等。把彼此接近的归为一类的聚 类方法也有多种，如一次性完成的系统聚类和逐步完成的动态聚类；基于分明 （硬）划分的聚类和基于模糊（软）划分的聚类；具体实施上有通过矩阵运算实 现的，也有应用图论方法实施的。
最新最详细化探数据处理与编图流程
最新最详细化探数据处理与编图流程⼀、指导思想成矿地质背景地球化学研究就是从地球化学特征出发，借助已建⽴的地球化学信息提取技术，充分利⽤地球化学调查所获得的海量数据信息，提取有关反应成矿地质背景条件的地球化学信息，并编制相应地球化学图及相应的推断解释图件，为资源潜⼒评价有关成矿地质背景的研究提供地球化学⽀撑。
⼆、⼯作内容（⼀）基础图件成矿地质背景条件的地球化学信息提取⾸先是要编制有关基础地球化学图件。
主要有：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万化探的数据处理成图方法探讨毕武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。
Tin法制作地球化学图
首先把EXCEL表Y坐标放第一列，另存为CSV格式。
1：50000地化图绘制：1、用记事本打开Ag.csv格式数据，将x（坐标符号），y（坐标符号），Ag （元素符号）改为notgrid，另存为Ag.det格式。
对于1：5万化探，x、y单位为m，成图比例尺为1：50000万，那么将x、y同除以50，就变成mm。
2、依次打开MAPGIS程序主菜单\空间分析\DTM分析（MAPGIS数字地面模型自系统），打开Ag.det文件。
3、打开Tin模型\快速生成三角剖分网。
4、打开Tin模型\删除三角剖分网边。
5、打开Tin模型\三角网内插网格化：左下角坐标设为最小整公里坐标，右上角坐标设为最大整公里坐标；1：50000地化图网格化间距实验证明设为X：125，Y：125最为适当，另存栅格化结果文件Ag.Grd.6、打开Grd模型\平面等值线图绘制：选择：等值线套区、绘制色阶、保留边界线、等值线光滑处理（高程度）、示坡线、线型。
线参数设定：默认。
注记参数设定：标注字体尺寸：2mm，小数部分位数：与该元素线值位数相同，其它默认。
(.VEL)。
以后重新绘图时，可直接装入该文件。
7、确定成图，保存点、线、面文件。
8、极值点标注：等值线图极值点标注，只能在“三角剖分-追踪剖分等值线后”，在“模型应用-高程点标注制图”子菜单下标注。
在标注输出设置中，选择“等值区内部极值点标注”，然后按自己的要求设网上收集资料一、建立数据表1、用电子表格建立成图元素表（不要表头，建议将平面坐标按Y、X、数据…列次序移到表前1、2…列，Y、X坐标分别为7位、6位整数）；2、另存为文本文件（制表符分隔）类型的新文件并取名(如文件名为茨竹箐原始数据表)。
二、Mapgiss投影转换1、打开Mapgiss程序主菜单→点击″实用服务″→点击″投影变换″→点击″P投影转换″→点击″U用户文件投影转换″→打开文本文件(茨竹箐原始数据表)；2、选″指定分隔符″→指定坐标X、Y列（注意X、Y列可能倒置）；3、点击″设置分隔符″→选″Tab″及″空格″→点击″数据类型″→将所有的字符串点为″双精度″→在″属性名称所在行″选框中选择数据（注意数据不能选″无″）→在″加入″列选定所需的坐标及元素值→确定；4、点击″点图元参数″→选择子图参数（建议选编号35的圆点的子图号，高度宽度一般均为1，颜色为1）→确定；5、选″不需要投影″→点击″数据生成″→确定→关闭；6、另取点文件名（如茨竹箐数据投影)→保存→退出。
6_化探数据处理与异常查证方法_201411
二、多源地学信息分析与变换
5、取对数 Xi’ = Ln (Xi+C)
式中:i=1，2，…,n；为样本数； C 为常数，防止接近于0的元素含量太负； 功能：对原始观测值xij 含量极差比较大的元素（如： Au 、 Ag 、 Hg 等成矿元素）进行变换后，使其服从对数 正态分布。
二、多源地学信息分析与变换
•
R或Q -型因子（Factor）分析
利用维变换或转换方法减少有效变量的数目或找到数据 的不变式。因子分析法或主成分分析（PCA）是依据变量空间 （R模式）或样本空间(Q模式)的主因子解，进行正交或斜交 变换，通过线性组合优化，在信息量小于15%的前提下，相关 变量x1，x2，…，xm（或相关样本y1，y2，…，yn）可以用数 目较少且本身互不相关的因子计量f1，f2，…，fj（j＜m或j ＜n）来取代，达到特征分类提取、优化组合的目的。 在R模式因子分析中，fj = βj1x1+βj2x2+ …+βjpxm 或写成 矩阵形式： F’=X’B 式中，B是回归系数矩阵，阶数为p×m。按最小二乘法求解， B满足如下正规方程组 RB = A 或 B =R-1A 因子计量（因子得分）表达式为：
2、综合异常图
①将相关元素叠臵形成多元素组合异常图，元素组合 根据点群分析、因子分析确定， ②按矿化类型元素组合编制
③多元统计综合异常图
因子计量编图：反映某一因子元素组合异常信息。
累加晕、累乘晕：按矿化类型元素组合衬值或标
准比值累加或累乘，计算异常下限，圈定综合异常。
累加晕、累乘晕的作用是强化低缓异常。
（三）特异数据处理
1、低于捡出限的数据取捡出限1/2 如：<3.0 = 1.50
2、高于捡出最高极限的数据取2—3倍值
地球化学学习主要流程步骤
地球化学实习主要流程步骤实习四概率分布型式的检验及背景值和异常下限的确定1.将实习用的数据新建一个文件夹，以自己的学号命名2.复制表4-1的数据，选择性粘贴→文本，然后保存。
文件→另存为→一共3个要点：保存的位置：自己建的文件夹里，必须和运行程序在同一级目录；文件名：推荐4-1；保存的文件类型：csv逗号分隔 (后缀必须是csv格式)；然后以同样的操作方式复制表4-2的数据进行相同操作。
3.运行BASIC程序，得出结果4.以写字板方式（EXCEL）打开运行结果.5.简单解释一下相关参数的意思6.以相同的方式运行4-2的数据，得出运行结果.7.演示做作业。
复制实习4的作业到另外一个空白的WORD文档里，单独以自己的学号命名。
具体方式为：姓名+学号+实习内容的顺序。
比如：200701010326朱国勋No.18.开始填表：首先求出均值和方差，用各组组上限值（Xi），计算，以4-1为例：均值为1.46，方差为0.26，用EXCEL编辑公式=ROUND((A2-1.46)/0.26,2)9. 据值查标准正态概率分布函数表，得到各Xi点的累积概率F（）。
10.计算Dn值，在F（）的基础上加减13.6分组间隔频数f频率f累计频率∑fμ值概率F(μ)%临界曲线（上）临界曲线（下）0.77-0.97 3 3 3 -1.88 3.01 16.61 -10.590.97-1.17 9 9 12 -1.12 13.14 26.74 -0.461.17-1.37 22 22 34 -0.35 36.32 49.92 22.721.37-1.57 32 32 66 0.42 66.28 79.88 52.681.57-1.77 20 20 86 1.19 88.3 101.9 74.71.77-1.97 12 12 98 1.96 97.5 111.1 83.91.97-2.17 2 2 100 2.73 99.68 113.28 86.0811.建立画图的文件。
地球化学图和原始图件的编制规范(正式)
a.采样点要分布均匀、合理。要使1：50000地形图上按1km2划出的采样格子（称小格）中的绝大多数（95%以上）小格内都有采样点分布，应做到基本不出现或很少出现连续3个以上空白小格。采样密度平均为1~2个点/km2时，小格内样品数一般不要超过2个。采样点分布均匀并不意味着将采样点布置在采样格子中央，而应将采样点布置在每一个格子中能最大限度控制汇水面积处（见图2）；
GB/T14839地球化学勘查技术符号
3总则
3.1区域化探是基础性、公益性、战略性基础地质调查工作之一。工作范围应覆盖我国全部可工作的国土面积。其调查成果，主要为区域矿产勘查中成矿远景区和中大型找矿靶区的圈定、资源潜力预测、地质填图与基础地质研究、环境与土地质量评价以及社会其它行业和领域的应用提供数十种元素含量空间分布的地球化学资料。此项工作按国际分幅部署测区。基本工作比例尺为1：250000，交通很不发达、工作条件困难的边远地区也可采用1：500000工作比例尺。它们的基本采样密度（以水系沉积物测量为例）分别为1点/1~4km2和1点/16km2。成果图件视不同行业、不同领域和不同群体对地球化学资料的需求，可制作任意比例尺和不同图示表达形式的图件（册）提供应用。
d.重要异常的查证：包括拟采用的查证方法、使用的现场分析手段和设备、查证程度、提交异常查证的成果图件、简报等。
e.技术、经济指标：包括设备、材料计划、经济预算、组织编制和工作进度安排等；
f.提交成果的内容与时间：包括计划上交的报告、图件等资料，提交时间等。
4.5区域化探项目都应单独编写设计书。设计书由项目负责人负责编写。
3.2按地理特点和区域化探工作条件，可把我国大致分为以下两类地区（见图1）：
a.内地及沿海中低山、丘陵地区；
利用MAPGIS处理化探数据文档
地球化学探矿工作方法及注意事项异常下限计算公式：都采用Excle表格计算(原始数据参照表三)平均值：=A VERAGE(C2:C36)方差：STDEV(C2:C36)异常下限：平均值+K值*方差计算过程中自动剔除大于平均值加3倍离差的高值和小于平均值减3倍离差的低值。
K值自选（1-3）剔除次数的多少取决于没有大于和小于这些数据。
表一表二表三等值线图生成方法步骤：利用Excel先做一个三角部分文件，格式如图：A1：要打上notgrid.A2为生成点文件的X坐标 B2为Y坐标 C2为高程代表化学分析结果。
文件需要另存格式如图：点击保存，关闭表格。
打开MAPGIS6.7空间分析-DTM分析-点击文件打开三角部分文件。
如图：文件打开后在菜单栏中找到Grd模型打开里面有一栏离散数据网格化打开。
如图：网格化方法用：Kring泛克立格法网格化，文件换名保存文件。
“选择打开”，选三阶漂移如图：其它内容不用改动，改完后按确定。
“搜索打开”搜索类型改为八方向，其它内容不用改动，改完后按确定。
如图：全部改完后按确定。
再一次在Grd模型里找到平面等值线图绘制点击后打开离散数据网格化后保存的文件，如图：把文件打开如图：设置：等值线套区、绘制色阶、保留边界线全部选上。
光滑：等值线光滑处理选上，光滑度选择高程度。
等值线定层一般定四层，第一层为背景值，第二层为异常下限值，第三层为中带值（异常下限值的2倍），第四层为内带值（异常下限值的4倍）。
注记参数都改为YES。
区、线参数参照图例规范设定。
制图幅面：改为数据投影变换如图：原数据投影参数与结果数据投影参数都要改变，设置结果与投影点位图一致。
全部设置完后点击确定。
生成图像如图：生成图像后保存点、线、区。
保存好以后再一次打开MAPGIS6.7图形处理-输入编辑打开刚才已保存好的点、线、区文件进行图式修饰。
点位图制作方法步骤：利用Excel表格把所有化探数据先处理一下如图：注：ID代表总序号、YH代表样号、Y代表X轴坐标、X代表Y轴坐标、Au、Ag 代表元素化学分析结果。
地球化学+地质填图教程总结
地球化学+地质填图常用教程总结一、1:5万土壤采样点位设计(水产沟幅为例)1、标准图框生成首先利用Mapgis—实用服务—投影变换—系列标准图框—生成1:5万图框.图框模式：高斯坐标实线公里网.投影参数：起点纬度491000. 起点经度1210000. 网间间距1KM.图框文件名：水产沟图框.图框参数输入：图框内容：图幅名称：水产沟幅.图框参数选择：输入并绘制接图表+绘制图幅比例尺+图框外图廓线.2、图形校正图片JPG转换成msi文件： Mapgis—图像处理—图像分析—文件—数据输入—添加文件(图片JPG)—转换—转化成msi文件.校正到标准图框内：Mapgis—图像处理—图像分析—文件—打开影像—msi文件—镶嵌融合--打开参照文件—参照线文件(标准框). 镶嵌融合—添加控制点—(先点击带校正图片某点，然后点击对应的标准框点，最后按两次空格键，选择是，以此类推，尽量多的找几处控制点，通常不少于13个控制点，要求残差<1为准)—镶嵌融合—校正预览—透明度调节--回车--镶嵌融合—影像几何校正—几何校正--重命名保存文件.3、图片数字化对地形线、河流、公路、村庄等地形要素进行数字化过程. 同时，准备大格号、图例、网格线、采样位置、重复样和监控样分析表格等文件.4、采样点位设计四幅图框：据前述方法，生成5万四幅合体的图框（注意修饰）。
将四幅图之间的界线用红线标出，根据阵列复制线，画出网格线，在网格线的交汇处定出大格号标注，然后建立采样点位文件，在上面点击采样点位置，辅助工具1—导入导出功能—点位置转属性. 辅助工具2，水系沉积物自动编号，0001大格号的左上角坐标为5488000,355000，一行上的大格号从0001—0004，点击确定，重新显示文件，出来小格号标注文件.5、重复样、监控样在每个单幅的采样点位图中，设置重复样和监控样，按从上到下，从左到右的原则，即a—c—b —d，查出每50个样品，其中1个是重复样，4个是监控样，45个是普通样品，导出excel表格之后排序，对其重复样和监控样设置，然后在单幅采样点位图中标出来.二、1:1万土壤采样点位设计1、基点基线准备在地形图基础上，可以在图框文件里输入一个点，作为基点（待设计完毕可以删除，例如基点为200/100，点号为200，线号为100），图框外部边缘作为基线。
详查阶段的化探数据成图处理探讨
１ 中的测 网密度 。 这主要 是考 虑与物 探 网布 设相 同 ， 但 从化探元 素分 布规 律和成 图方 面来说 ，正方 形 网格 应 该更合适 。
表 １ 工作 比例 尺 与测 网密度
图 １ 原始 数据转网格数据 绘制的地球化学图
我们使用的网格数据绘制等值线图一般要求正方
形 网格 ， 就需 要用 到网格化 处理 工具 ， 这 如距离平 方导
随着 ｌ： ５万 区域化探 普查 的实施 ，在 ｌ： 万 区 ５
域 化 探异 常 范 围内开 展 的 １：１ 万化 探详 查 也越 来 越 普 遍 ，目前 １： 万 化探处 理一般 是参 照 １： １ ５万 区域 化 探规 范执行 ，可是 由于 多种原 因成 图效果有 时却不 尽 理想 ，下 面就利用 乌鲁木 齐金 维图文信 息科 技有 限
格化 。
图 ４ Ｋ ｉｎ ｒ ｉｇ网 格 化 绘 制 的地 球 化 学 图 ｇ 图 ２ 距 离 平方 导 数 加权 网格 化 经 圆滑 处 理地 球 化 学 图
３ 结束语
（） 查 阶段化探 的测 网宜 采用 正方形 网格 。 １详
（ 网格化遵循 的一个原则是尽量使 网格化后的 ２ ） 坐标点 位与 原始数 据坐 标点 位接 近 ，网格 化搜 索 的范
图 ３ 三 角剖 分 方 法 绘 制 的地 球 化 学 图
验 ， 中选 出满 意 的 ， 从 然后再 统一 批处 理 。
参 考文献
［］ Ｐ １ Ｇ Ｓ和 Ｇ Ｓ技术在 １： Ｉ １万土壤地球化学测量 中的应
对这 种非规则 的数 据也 可 以采 用 三角剖 分方 法绘 制地 球化学 图 ， 图 ３ 如 。
围越小 数据损 失越 小 。
（ 网格间距要小于等于最小 的原始点距 ， ３ ） 如果要 保 留更 多 的细节 ， 也可 以将 网格距减 小一倍 。 如果细 节
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
3、化探中常用的数据处理方法：
化探数据处理过程中面对大量原始数据主要依靠现代 数学一统计方法和信息技术，分析方法多种多样，有单变量 分析：（如趋势面分析，稳健统计等）和多变量分析（如判 别分析，簇群分析，相关分析，因子分析等）之别；就处理 所依据的统计分布律而言，有依赖于正态律的参数统计法和 不依赖于它的非参数统计法，属于后者的模式识别技术已有 很大发展；就处理所针对的数据变化特征而言，有研究线性 关系和非线性关系的两大类方法。近些年来，根据地球化学 数据的变异特点而采用的克里格法等研究非线性变化的方法 迅速发展，混沌、分形和智能神经网络等新型数据处理方法 也得到广泛应用。
6、多元地学空间数据管理与分析系统 （GeoExpl）应用
主要步骤： 建立数据工程、应用工程→建立基础 数据库（数据的导入、备份）→数据分析 检索→离散数据网格化→网格数据二维空 间分析（图件制作、生成）
谢谢大家！
3、化探数据是以多元素或多变量为特征的。如果参加分 析的数据含有已知因素（如矿或非矿的作用）时，数据处 理的结果可给出明确的地质解释，否则所做的地质解释就 含有较大程度的推测性。 4、应该充分注意地球化学数据中少数高含量值对统计计 算结果(如背景和异常的数值)的扰乱，尽量使用以中位数 为基础的稳健统计学代替以平均值为基础的传统统计学处 理数据。 5、保证地球化学数据处理质量的关键，是所用数字模型 的原理和特性，搞清所用参数和处理结果的地质一地球化 学意义，避免朦朦胧胧地做数字游戏
3、异常下限（背景上限）：指划分异常与背景的临界值, 大于此值为异常范围.特征元素背景含量不是一个特定值， 一定范围内元素含量变化的一系列数值，其中最大值即为 背景上限值，超过临界限值，为异常范围，异常下限异常 范围内最低值。 4、异常衬度：又称为异常清晰度.异常衬度定义为异常强 度与平均背景或异常下限的比值.有多种表示方法: C＝xa/xb C:异常衬度; xa和xb相应为异常平均值和背景 平均值。 在矿区化探工作中,在一个小范围内出露的岩性差异较大, 使得元素的背景值差异增大,此时的岩石测量就不宜使用 统一的异常下限圈定异常.而采用异常衬度圈定异常可以 更为确切地反映不同背景岩石中的地球化学异常.
5、主要化探图件类型
原始数据图：包括采样点位图、分析元素原始数据图。 等值线图：包括主要成矿和伴生元素的等值线图。 地球化学图：包括主要成矿和伴生元素的地球化学图。 异常图：包括单元素、综合元素异常图 解释推断图：包括各类解释、推断、异常剖析图 各种辅助专题图件：包括数据直方图、聚类分析图等。 其它：移动标准化地球化学图、综合异常分区图、综 合异常结构简化图等
第二部分 化探数据处理及成图
1、地球化学数据处理 2、化探数据处理解决的主要问题 3、化探中常用的数据处理方法 4、化探数据处理中应注意的几个问题 5、求背景值和异常下限的常用方法 6、主要化探图件类型 7、多元地学空间数据管理与分析系统（GeoExpl） 应用
1、地球化学数据处理
化探工作中对所获得的多种元素分析数据以及量化的相 关学科数据进行加工、分析、成图和解释的工作。其核心内 容是有关地球化学背景与地球化学异常的合理区分问题,从大 量数据中发现与矿有关的地球化学信息是化探数据处理的目 的。
（三）位置的相对性 无论是土壤测量还是水系沉积物测量所获得的异常， 往往与异常源都会发生不同程度的位移。这种位移与表生 介质本身的位移程度和采样的布局有关。特别是水系沉积 物异常的位移更为明显，可达几公里甚至更大的距离。因 此，查明异常与异常源的空间关系，就成为异常查证中的 首要任务。 （四）表生作用带来的复杂性 不同的景观条件下，表生地球化学作用会有很大的 差异，制约了元素在表生环境中的分散和富集。因此，只 有在同一景观内，异常才有较好的可对比性。地理景观不 同，表生地球化学作用也就不同，元素在地表迁移、分散、 富集的规律也就不同。在异常对比和解释上，除考虑引起 异常的原生因素（地质背景、矿床类型）以外，在一定程 度上必须注重异常所处的地理景观条件及表生地球化学环 境。一般来说，只有同一地理、地质景观区的区域化探异 常才有对比研究的基础。
2、化探数据处理解决的主要问题：
①研究采样和分析中的误差，优化采样布局
②抑制数据噪音，突出主体趋势
③揭示多种数据的内在联系，提取隐蔽的有用信息 ④显示数据空间分布模式，编制地球化学图件
⑤异常对比、分类、评序，等等。
3、化探数据处理中应该注意的问题： 1、地球化学数据通常蕴含多种有用信息并伴随某些 不规律的变化，同时在数据获取过程中还存在分析测定误 差，这些使化探数据的复杂性增加了，在化探数据处理中 要将这些不规律成分和分析误差除去。 2、找矿信息总是同地球化学异常相联系的。最普通 的化探数据处理是对一组化探数据计算出背景值和变化范 围（如用平均值和标准离差来衡量），据此确定出地球化 学异常的下限值。当地球化学背景随着地理位置出现趋势 变化时，要相应地采取适当的处理方法以便获得随地理位 置而变的背景值和异常下限。
区域化探异常特点 （一）信息的直接性或直观性 区域化探所圈定的异常一般对矿（化）体均具直观性。 因为有什么矿（化）体，就会反映与之相关的成矿或伴生 元素，例如：金矿体上会出现一组与Au有关的元素；不同 元素组合特征反映在地球化学图上，可以直观地判别各元 素地异常特征，如：浓集中心、浓度分带、元素组合特征、 异常规模、形态和变化趋势等。 （二）组合与分带 矿（化）体的地球化学异常是由一些或一系列指示元 素或组分组成的，这是由矿（化）体的化学组成决定的； 这些元素或组分在空间上呈现一定的浓度分带和组份分带， 这种分带是成矿时的地质和地球化学条件决定的。 在地球化学找矿工作中，根据指示元素或组份的组合和含 量特征，可以判定引起异常的可能矿种和矿床类型；根据 指示元素或组份的相对强度和分带，可以判定矿（化）体 的剥蚀程度，进一步还可以研究空间上矿产的成矿系列等 问题。
最为常用的方法有：回归分析方法、判别分析方法、 簇群分析方法、因子分析方法、对应分析方法、典型相关分 析方法、方差分析方法等。
4、求背景值和异常下限的常用方法：
(1)剖面图解法 (2)直方图解法 (3)概率格纸法确定背景值和异常下限 (4)计算法 计算正态分布母体的背景值C0和标准离差方差σ,按 公式: Ca=C0±Kσ K=1-3 K为信度系数,选择信度α=0.05,K=1.96,α=0.01,K＝ 3,α=1,K=1.65.一般选择2。 即平均值加2倍标准差确定异常下限。
如何理解： 1. 背景和异常的相对性-正常与异常相对 2. 背景多介质性-对任一特定元素,背景值随它所赋 存的介质性质有变化。 3. 多指标性 4. 空间上的不均匀性-不是固定值,正常状态也有变 化.元素在任何一种介质中的分布不均匀，背景范围 内元素含量围绕一常见值呈不规则起伏。 地球化学异常消失，并非不存在引起异常的元素，指 它的变化被背景起伏隐蔽。
化探数据处理与成图
第一部分地球化学勘查中的几个概念 第二部分化探数据处理及成图
第一部分 地球化学勘查中的几个概念
第一部分地球化学勘查中的几个概念
1、地球化学背景：指无矿化地球介质中元素的 正常含量分布.在一定范围内的同类介质中广泛存 在的地球化学指标(环境)特征.可理解为调查对象 影响以外的地区,地球化学指标在更大范围内呈正 常状态的情况. 2、地球化学异常：在给定地区对正常地球化学 模式的偏离，即量值超过背景值-地球化学异常, 简称异常.成矿过程或成矿后风化过程中,成矿元 素及其相关元素组合在矿体周围/附近形成相对 富集的情况.地球化学异常最早用于勘查地球化学 中,指某一区段的地球化学特征(元素含量高低和 存在形式)明显不同于周围无矿背景区的现象.