陈仕炎
(长沙有色冶金设计研究院有限公司,长沙410011)
MapGIS地理信息系统是武汉中地信息工程有限公司研制的,集数字制图、数据库管理以及空间分析为一体的空间地理信息系统。它以其功能强大而广泛应用于城市规划、测绘、地质勘查、资源管理等领域,并已成为我国各领域数字化建设的首选软件之一。在地质找矿方面,该系统在提升地质勘查技术手段上发挥了重要作用,大大提高了工作效率,而且保证了空间数据的精确性。
本文以新疆恰尔隆一带1∶5万水系沉积物地球化学测量为例,从背景值和异常下限的确定、不同地球化学图件绘制等方面简述了化探数据处理的方法和步骤,通过MapGIS软件与Excel电子表格结合使用在化探数据处理过程中取得了良好的效果。
工作区位于新疆西昆仑山西段北坡,塔里木盆地南缘。地势南高北低,海拔一般为3 000~5 500m,最大相对高差在1 500~2 000m左右。区内水系较发育,均属内陆河的源头水系,主干水系为塔什库尔干河和恰尔隆河。支流水系比较发育,水系组合格局多呈现为树枝状(调查区北部)及网格状(调查区中南部),河流水量比较充沛,为开展水系沉积物测量的有利地段。
本次化探取样共计2 706件,化探样品测试由新疆地矿局中心实验室测试完成,分析项目包括:Au、Ag、As、Sb、Cu、Pb、Zn、W、Sn、Mo、Bi、Ni、Co、U、Th等15个元素。
成图前逐点采集样点的平面坐标。数据排列顺序为X坐标、Y坐标、元素含量(见表1),图面坐标与实际坐标要进行转换。分别建立各元素含量数据的电子表格文件,在Excel电子表格中另存文件,文件类型为*.txt。
考虑到方法的实用性、有效性和易操作性,通过几种方法在工作区的试验对比,迭代法确定的背景值及异常下限较低,更有利于突出弱异常。因此,本工作区背景值和异常下限的确定选用迭代法。
应用迭代法在Excel中计算异常下限的过程:(1)计算全区各元素原始数据的均值(X1)和标准偏差(Sd1);(2)按X1+nSd1的条件剔除一批高值后获得一个新数据集,再计算此数据集的均值(X2)和标准偏差(Sd2);(3)重复第二步,直至无特高值点存在,求出最终数据集的均值(X)和标准偏差(Sd),则X做为背景值C0,X+nSd(n根据情况选1.5或2,3)做为异常下限Ca。
表1 测区内部分水系沉积物样品原始数据
表2 工作区元素地球化学特征值及参数表
(1)进入 MapGIS 6.7主菜单,点击空间分析模块中“DTM分析”子模块。
(2)在“Grd模型”下拉菜单中单击“离散数据网格化”命令,根据需要修改网格参数调整网格疏密程度,选择网格化方法,导入*.txt文件,生成网格数据(*.Grd文件)并保存(图1)。
(3)在平面等值线图绘制命令中导入*.Grd文件,进行等值线图的绘制。
图1 离散数据网格化显示视窗
(4)打开生成的网格数据文件,在设置等值线参数对话框中(图2),选择等值线套区、保留边界线、绘制色阶和等值线光滑处理等,并分别对等值层的起始值、终止值、步长值和线参数等项进行设置,点击“确定”按钮即可显示元素地球化学异常分布图,结果如图3所示。
(5)保存所生成的点文件、线文件和区文件。
图2 设置等值线参数视窗
图3 测区部分地段金元素地球化学异常分布图
依照上述生成元素地球化学异常分布图的步骤,并且根据工作区地球化学特征,确定元素组合,把元素组合中各元素的单元素异常添加在一起,然后在MapGIS编辑子系统创建工程,并在工程管理器窗口依次添加各元素的地球化学异常分布图等相关文件;最后点击窗口中“1∶1”按钮,在窗口位置便可看到各元素地球化学等值线图叠加在一起,形成地球化学综合异常分布图。测区部分地段综合异常图见图4。
图4 测区部分地段化探综合异常图
Excel电子表格具有强大的数理统计功能,可以快速统计计算出元素地球化学特征参数及异常下限,而MapGIS软件在化探数据处理中具有独到的功用,如利用它自动生成各元素的地球化学等值线图,直接编制综合异常图等,具有操作简单、图件精确、快速高效的特点。
[1]袁义生,刘应忠,罗明学,等.应用 MAPGIS制作地球化学图单元素异常图及综合异常图[J].贵州地质,2007,24(4):156.
[2]张钊,韦龙明,陈三明,等.MAPGIS在地质填图及化探数据处理中的应用[J].物探化探计算技术,2010,32(2):221-224.
[3]赵荣军.不同方法在栾川北部化探数据处理中的应用[J].地质与勘探,2006,42(3):67-71.
[4]吴锡生.化探数据处理方法[M].北京:地质出版社,1992:25-26.
[5]李明,崔永文,李俊,等.MapGIS与Excel相结合在辽北化探数据处理中的应用[J].科研探索与知识创新,2009,12(下):83-84.
[6]纪宏金,林瑞庆,周永昶.关于若干化探数据处理方法的讨论[J].地质与勘探,2001,37(4):56-59.
[7]宋运红,贾大成,李健,等.河北承德三沟地区化探数据处理与异常特征[J].吉林大学学报:地球科学版,2008,38(增):167-169.
[8]戴慧敏,宫传东,鲍庆中,等.区域化探数据处理中几种异常下限确定方法的对比——以内蒙古查巴奇地区水系沉积物为例[J].物探与化探,2010,34(6):782-786.
[9]春乃芽.如何利用Excel处理化探数据[J].物探化探计算技术,2006,28(3):272-276.
[10]杨世莹.Excel数据统计与分析范例[M].北京:中国青年电子出版社,2005:18-27.