勘查地球化学实习报告021111班.docx

勘查地球化学实习报告 秭归5万水系沉积物地球化学异常图绘制 021111班 -- 一．离群样品（剔除）及数据统计 首先Grafer软件打开.text格式的秭归五万原始数据，对Cu所在列数据进行降序排列 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异常点异常参数统计 　异常1号异常2号异常3号异常4号异常5