《化探数据处理成图过程》由会员分享，可在线阅读，更多相关《化探数据处理成图过程（16页珍藏版）》请在金锄头文库上搜索。
1、化探数据处理成图的过程毕武1、2段新力1、2黄显义1、2袁小龙1、2彭仲秋1、2李永华1、21.乌鲁木齐金维图文信息科技有限公司，新疆，乌鲁木齐，8300912.新疆地矿局物化探大队计算中心，新疆，昌吉，8311000 前言原测数据水文地质数据处理数据预处理重磁数据处理化探数据处理数学地质数据处理点位数据图打印输出数据检查数据整理数据网格化投影变换格式转换电法数据处理数据处理图件绘制输出数据输入图1 金维地学信息处理研究应用系统(GeoIPAS)图件制作流程图图像输出文本文件图形文件剖面图等值线图场强图三元图三维立体图点位符号图GeoIPAS软件用户群不断扩大，由于各用户对系统的熟悉程度不同，
2、对软件功能了解不够，有必要分专题将GeoIPAS处理数据及成图过程做一系统总结，下面就化探数据处理成图的过程做一总结。1 处理步骤化探处理的成果包括：(1)参数统计表；(2)R型聚类分析-谱系图；(3)重复样三层套合方差分析或者重复样合格率计算结果；(4)点位数据图；(5)地球化学图；(6)直方图；(7)组合异常图；(8)综合异常图；(9)远景区划图；(10)单元素异常参数统计(附表册)；(11)异常剖析(附图册)；(12)综合异常登记卡(附表册)。在GeoIPAS系统中，化探数据处理分为以下几个主要步骤：1.1 数据检查数值检查，坐标检查，重复样坐标检查。1.2 分析处理重复样三层套合方差分
3、析、重复样合格率计算、化探特征参数统计、化探背景值分析、R型聚类分析、因子分析。1.3 数据分析数据变换；衬值、累加衬值；数据累加、累乘、比值；异常归一化。1.4 网格化离散数据网格化、XYZ数据转网格数据1.5 成图点位数据图、彩色等量线图、直方图、组合异常图、单元素异常图、综合异常图、剖析图。1.6 单元素异常参数统计1.7 综合异常登记卡开始数据输入数据预处理网格化成图异常统计数据分析处理结束分析样品与坐标对应数值检查坐标检查重复样提取采样编码统计三层套合方差分析或重复样合格率计算化探特征参数统计R型聚类分析离散数据网格化XYZ数据转网格数据点位数据图直方图地球化学图组合异常图单元素异常
4、图综合异常图剖析图TXTGRD单元素异常参数统计图2 化探数据处理成图流程TXTTXTTXTTXTTXT 矢量GRD矢量TXT 矢量CSV2 具体处理过程2.1 数据检查我们的数据处理工作从化验室提供的样品分析报告开始，项目要提供坐标和样品对应的分析数据，坐标我们一般取实际工作中的米单位，系统中默认东西向横坐标为X坐标，不加带号，南北向纵坐标为Y坐标，需要时还要提供样品对应的地质编码，我们拿到这个数据后首先进行数据检查，以确保数据中不出现写错、漏填、负数、0、等字符，如果有这样的情况要找实验室给予纠正。数据准备好后，我们要把数据转换成TXT后缀的文本文件，这就做好了处理前的准备数据工作。2.2
5、 分析处理2.2.1 重复样三层套合方差分析一般是从分析样品的结果中挑出来重复样的分析值，每组四个样品，按如下顺序排列：11 第一次采样第一次分析12 第一次采样第二次分析21 第二次采样第一次分析22 第二次采样第二次分析结果：三层套合方差分析成果-元素:Cu方差来源离差平方和自由度平均离差平方和F-检验值地点间13.224180.73511.152样品间1.252190.0661.482分析间1.689380.044三层套合方差分析成果-元素:Au方差来源离差平方和自由度平均离差平方和F-检验值地点间0.573180.03210.781样品间0.056190.0030.483分析间0.23
6、2380.006查表：Fa0.05(18,19)=1.35 Fb0.05(19,38)=1.22Cu：11.152 1.35 1.482 1.22 即：F1Fa F2FbAu：10.781 1.35 0.483 Fa F2f0.05()理论值,说明采样与分析误差不显著，此次分析结果可靠，Fbf0.05()理论值,说明采样与分析误差显著，且以采样误差为主。例如报告中的描述：先列一个表格说明所有元素的F1和F2；重复样三层套合方差分析，39种元素F1值远超过临界值，说明采样与分析误差变化均低于区内地球化学自然变化，数据是可用的。F2值表现为两种情况，Ag、Pb、P、U、MgO、La等6个元素F2值
7、小于临界值，说明采样误差不大于分析误差，以分析误差为主；其余33种元素F2值均大于临界值，说明误差应以采样误差为主。2.2.2 地层编码统计从样品相对应的地质图中自动读取地层的颜色号作为编码，实现每个样品的地层自动编码。在GeoIPAS系统中的化探-化探特征参数统计-地层编码统计其处理结果是在输入的数据文件中增加一列，过程是将坐标投影到地质图相应的投影参数上，以相同地层颜色来套取数据，取地层的颜色号为落入此地层的采样点的地层编码，所以要想一次做好，要先对地质图做个预处理，最好是按编码顺序，将颜色号改成相应的顺序号，比如Q的颜色改为1，N的颜色改为2，E的颜色改为3，等等。相对应的出来的结果就是
8、按顺序的地层编码了，可以直接用了，绘制的直方图也是这个顺序了。2.2.3 化探特征参数统计化探参数统计成果-元素:Cu地质单元样品数(N)平均值(X)中位数(M)标准离差(S)变化系数(Cv)最大值最小值浓集克拉克值(C)致矿系数(Z)1.算术全区210541.91934.7328.6280.683250.492.491.1038.6162.对数全区21051.5241.5410.3070.2012.3990.3960.042.2553.算术剔除高值全区203939.01234.0222.9730.589108.722.491.02710.7344.对数剔除高值全区21051.5241.541
9、0.3070.2012.3990.3960.042.255注：浓集克拉克值可以选用系统提供默认的，也可以直接输入用户所在不同地区的浓集克拉克值。这个值可以保存为文件，以后导入使用。2.2.4 R型聚类分析从谱系图中可以看到，Cr、Ni、Co相关性最好，我们在做组合异常图时可以参考谱系图来做元素的组合，但这个也不能完全作为组合参考，我们通常还是以常用组合方式来做组合异常图，如：Cu-Pb-Zn；W-Sn-Mo；Cr-Ni-Co；Au-Ag-As-Sb。2.3 网格化2.3.1 工作比例尺与测网密度DZT 0145-94土壤地球化学测量规范 表1 工作比例尺与测网密度工作阶段比例尺矩形网格正方形网
10、格点km2线距(m) 点距(m)点线距(m)区域调查1：200 000 2 000500 1 000 11100 000 1 000250500 5001 000 14普查1：50 000 500100250 250500 4201：25 000 25050100 125250 1680详查1：10 0001002050 50100 1005001：5 000501025 25501：2 000 20510从上表可以看到，化探在实际工作中一般和物探工作一起工作的常采用矩形网格，但在实际成图中却是用正方形网格的好，不用经过网格化处理，没有数据损失，也符合绘图的正方形网。1:20万化探在全国范围已
11、基本完成，只在新疆等高山区有少量的没有覆盖全，最常采用的是组合样，即：1点4km2，少量的有单点样的，这主要用转网格数据的方式来将数据做网格化处理。单点样的网格化一般采用500500到10001000的网格距，搜索半径30003000到50005000。组合样做转网格数据，注意起始坐标为奇数公里，间距为2000，然后要进行一次补空扩边处理，这主要是处理数据成图到图框边。1:5万化探是现在化探普查的主要工作，一般在山区采用自由网的水系样，土壤测量一般用正方形网格采样，1:1万化探作为1:5万化探异常检查的主要工作手段，常采用10040的矩形网格，少量有10020、20040、15040等加密或者
12、抽稀采样。也有采用物探布设的测网非南北方向的测线。2.3.2 下面用几个问题来具体说明。问题1：1:5万土壤测量，网距330300的怎么设置搜索半径和间距？ 5万化探，网距500*500的，用GeoIPAS做地球化学图的时候，怎么设置搜索半径?答：a、关于网格距，最好采用正方形网格，网格距为野外采样的点线距的一半，比如野外采样为500500，网格化的网格距选250250就可以了，10040的选为5040或者5050。b、网格化的搜索半径选择需要一定的经验，建议你首先看看点位数据图，选择一个既能将所有数据联系起来，又不至于网格后把未采样区填充的参数，一般是网格距的2-3倍。问题2：在进行1:5万
13、化探采样时，乙类异常区采样密度设计为12个点km2，丙、丁类异常区及其它基岩区采样密度设计为个8个点km2，网格化时X Y 坐标的间距、行列数如何选？答：一般1:5万化探点线距250250米或者500500米；搜索半径以满足丙、丁类异常区网格后只有个别空白区或者无空区，搜索半径一般10001000到20002000米；如果对网格化参数设置不太清楚，可以选择“计算网格化参数”；网格化后如果异常比较细碎零乱，可用数据预处理二维数据滤波处理周围点平均圆滑处理，X和Y方向点数选择3，原点倍数选择2-4倍。问题3：请教如下问题 1:2.5万土壤测量，线距250m ，点距40怎么设置搜索半径和间距？答：搜
14、索半径500*500，间距40*40或者50*40，注意最好与原始点位重叠.问题4：1:1万化探，南北向测线，10040网格距。答：这样的数据本身符合网格数据的格式，只是以TXT文本格式保存，所以我们可以借助GeoIPAS软件的“数据预处理”“数据格式转换”“XYZ数据转为网格数据”功能，这里注意数据坐标的起始值、终止值和网格间距，网格间距就选择10040，起始值和终止值要保证数据的原始点位都与网格点重合。这里有一个问题，就是有些项目用的是GPS坐标，不一定是与网格点位完全重合，一般我们建议用规则网的坐标。圆滑方法：周围点平均圆滑；圆滑参数：X方向点数：3，Y方向点数：3，原点倍数：3问题5：我的测网不是正南北向是45度方向，用咱们的软件进行网格化后变成正南北的，而且图幅形状不是原来测网的形状。请问这样的问题怎样解决?1:1万
