计算地球化学实习3

计算地球化学实习3

1、统计异常参数（异常面积、异常强度、衬度、极大值、异常规

模）

a)用实习1和实习2已作好的地球化学异常文件（Mapgis点、

线、区文件及采样点文件）或新产生之。

b)异常面积＝区面积×100/1000/1000平方公里

c)异常强度＝异常范围内Mo的平均值

d)异常衬度＝异常强度/异常下限值

e)异常规模＝异常面积×异常衬度

2、用Surfer的grid→convert将Mo.grd转换为二进制文件(Binary)

3、MAPGIS→空间分析→DTM分析→GRD模型→平面等值线绘制→

设置参数：等值线绘制＋绘制色阶＋等值线光滑（高）＋等值线定层（层值0、1.5、3、6、）

“确定”

保存文件：文件→另存数据文件于→（分别存点、线和区文件为Mon.wt、Mon.wl 和Mon.wp）

新建工程文件：MAPGIS→图形处理→输入编辑→添加点、线和区文件

删除背景区的内容（图中非红色部分）

保存工程文件（Mon.mpj）、保存点、线和区文件

打开工程文件 Mon.mpj

添加采样点文件，并根据点分布特征新建裁剪框区文件（Monbox.wp）

点平差使多边形封闭（第1点与最后1点重合）→输入区

保存MonBox.wp文件后将其删除之。

工程裁剪：其它→工程裁剪→

选择要剪裁的文件→添加→点击其中1文件→生成被裁剪工程→装入裁剪框→开始裁剪→右下框内按鼠标右键→复位窗口→关闭。

在“我的文档”下打开工程文件

添加图框文件和采样点文件

删除单点异常的区并保存区文件

合并区的外、中和内带：区编辑→区编辑→合并区（不保存区文件）

区面积：

ID 面积（图面）实际面积（km2）

51 7381 0.7381

16 71.6 0.00716

15 97.3 0.00973

25 203.5 0.02035

30 849.2 0.0849

将Mon.wp另存为Mon1.wp

修改采样点文件的属性结构，将Mo的“字符串”→“浮点型”：

“OK”

存采样点文件

空间分析→空间分析→装入区文件“Mon1.wp”和点文件“sampling.wt”

“点对区”相交分析

“确定”

将相交分析产生的“Mon2.wt”转换为“dbf”文件：

“文件转换”→装入点文件Mon2.wt→输出→输出点属性(dbf)→

“确定”

在Excel下打开“Mon2.dbf”

“数据”→“筛选”→“自动筛选”→“ID0”→“51”

从中可以看出有171个点位于该异常范围内：用average算出异常强度为7.32×10－6；

用Max找出极大值为85×10－6；

异常衬度＝7.32/1.5=4.88

异常规模＝0.738×4.88＝3.60

同样的方法可以算出其它异常参数。