勘查地球化学实习报告021111班 - 图文

勘查地球化学实习报告

秭归5万水系沉积物地球化学异常图绘制

021111班 20111000---

一．离群样品（剔除）及数据统计

首先Grafer软件打开.text格式的秭归五万原始数据，对Cu所在列数据进行降序排列

1） box plot盒形图方法

箱线图作为描述统计的工具之一，其功能有独特之处，能直观明了地识别数据批中的异常值。异常值常被定义为大于外部上线Q3+1.5IQR或小于外部下限Q1-1.5IQR的值（Q1为内部下限、Q3为外部上线、IQR（Q3-Q1）为四分位距），本作业只考虑上部离群点作为异常，具体做法如下：

对Cu所在列数据进行box-plot绘制，由上至下不断删减列表中的数据，逐步剔除落在box图外部上线的离散异常点直至完全。此法所得标准方差7.76、中位数22.4、背景值23.4、异常下限为43.4、Q3=28.2、Q1=18、IQR=10.2；也可对剔除离群值后的样品统计取K0+2S（23.4+2*7.76=38.92）为异常下限，可知比前一种方法计算所得的异常下限偏小。（本次实习按39为异常下限进行作图，另外也可采用excel宏VBA编程进行箱图分析，直接完成数据统计）

2）调整坐标值：恢复原始数据的排列方式，调整坐标值适应相应比例尺，对X,Y列数据进行变换X=X/50,Y=Y/50.另存为秭归5万2.txt,导入excel调整后到处存为秭归5万3.txt

二. Surfer数据网格化

a.) SURFER→网格化→数据→打开数据文件秭归5万3.txt→调整X，Y，Z列分别为X,Y,Cu列→克里克法→.grid文件 b). 过滤：网格化→过滤→改输出文件名为Cu3，输出时存为BINARY(二进制格式)

三. MAPGIS异常图制作

a).MAPGIS→空间分析→DTM分析→平面等值线绘制→加入Cu3.grid文件(调整参数：等值线套区，光滑程度高，原始数据范围，外39、中78和内带值156及最大值和相应充填颜色(内178，中175，外190) → 另存数据于（点Cu3.wt、线Cu3.wl、面文件Cu3.wp）

b).建立工程文件：MAPGIS→输入编辑→新建工程文件→添加点、线、面文件(由上到下面线点排列 →另存工程文件Cu.MPJ

c) 建立采样点文件:MAPGIS→实用服务→投影变换→P投影转换→U用户文件投影转换→打开秭归五万3.txt→选定行(第2行) →按指定分隔符→(设置分格符→ TAB键,属性名称所在行设为第一行,序号1坐标XY不选→确定)→选定X第3，Y坐标第2列→选定点图元参数(子图号184,高度3) →不投影变换→数据生成→1：1选定NONAME.WT确定→ 另存采样点文件为Sampling.WT

d）异常编号：加入 Sampling.WT给异常编号从Cu1-Cu5（字大小为7）

e)制作图框：实用服务→投影变换→系列标准图框 →键盘生成矩形图框 →（任意公里矩形分幅→ 绘制实线坐标线 →比例尺1：50000→图廓参数X：488.05-501.95 Y：3402.05-3412.99 →XY网起始值489,3403,网间隔1km →点参数高宽4→图框文件名：图框5万→移动放大图例（整块移动+变换线XY，基点坐标10067，68036，XY约束X设为12+区变换弧段输入同样参数）→加图名，比例尺及责任表→工程输出页面设置为系统自动检测→保存，输出，生成JEP文件

f)修改点文件sampling。wt-编辑点属性结构(字符串改浮点，字段20，小数点1保存WIN7系统要先关闭TASKHOST进程) →CU.WP区文件区合并更改ID为1-5 →工作区换名另存区文件CU33.WP

四．统计异常参数

a).异常面积=图上面积×50×50/1000000(平方千米) b).异常参数

i.空间分析→空间分析→装入SAMPLING点文件和区文件CU33.WP(合并后的面文件)

ii.点对区相交分析→新点文件

iii.文件转换→O输出输出点属性 点属性文件C33.DBF) IV.EXCEL统计相应ID0的变量统计值。

现在以2号异常为例计算水系沉积物Cu异常参数： 异常下限：39

异常样品数：由数据筛选所得，为49 异常峰值（极大值,=MAX（C11：C59））：381.62 平均异常强度（平均值，=AVERAGE（C11：C59））：65.81 异常衬度（本次取：异常强度/异常下限）：65.81/39=1.69 异常面积（图上面积×50×50/1000000）：5487.44*2500/1000000=13.72km2 异常规模（异常衬度×异常面积）：1.69×13.72=23.15 km2

表一：秭归5万水系沉积物Cu异常点异常参数统计 异常峰值CMAX 异常强度C 异常样品数N 异常下限Ca 异常衬度 异常面积S 异常规模 异常1号 52.70 43.35 2.00 39.00 1.11 0.39 0.43 异常2号 381.60 65.81 49.00 39.00 1.69 13.72 23.15 异常3号 74.00 52.80 4.00 39.00 1.35 0.73 0.99 异常4号 异常5号 51.10 54.00 46.66 45.27 7.00 3.00 39.00 39.00 1.20 1.16 1.20 0.45 1.44 0.53 注：导出数据见附件C33.DBF/CU33.xlsx，水系沉积物Cu异常参数统计和计算过程见附件Cu异常参数统计.xlsx。

五．心得体会

本次实习是化探的第一次实习，相对于资源信息工程和矿床统计等的实习来说工作量显得比较大，作图过程中涉及Excel，Grapher,Sufer,Mapgis等软件的操作；过程比较复杂，软件制图也显的比较专业化。对于还没经过Mapgis专业训练的我们来说难度的，先是经过老师的细微讲解和请教同学，基本弄清了每一步的大致流程。虽然一开始并不顺利，许多操作按键不熟悉，也有漏掉一些步骤，但经过反复调试不断琢磨，最后顺利的完成了水系沉积物的制作，算出了异常参数，独立地完成了作业和报告编写。、、、、、、最后，谢谢老师的认真讲解。

编辑：021111班 20111000---