维普资讯 http://www.cqvip.com 2007年24卷 贵 州 地 质 Vo1.24 No.2(To1.91)2007 ?156- GUIZH0U GE0L0GY 第2期(总第91期) 应用MAPG l S制作地球化学图 单元素异常图及综合异常图 袁义生 ,刘应忠 ,罗明学 ,何彦南2 (1.贵州省地质调查院,贵州贵阳550005;2.贵州省地质矿产局105地质大队,贵州贵阳550018) [摘要]利用MAPGIS离散数据网格化模型绘制地球化学图,通过四种不同的网格化方法 效果对比,在元素含量分析数据极差相对较小时,采用Kring泛克立格法网格化模型绘制的 地球化学图效果佳。分析数据极差相对较大时,采用距离幂函数反比加权网格化模型绘制的 地球化学图效果较好。两种数据网格模型绘制的地球化学图均达到化探规范中对地球化学图 的要求,因此提出了应用MAPGIS制作地球化学等含量线图,利用地球化学等含量线图的属 性,绘制单元素异常图及综合异常图。 [关键词]GRD模型;等量线;色块图;属性;MAPGIS [中图分类号]TP302.4 [文献标识码]A [文章编号]1000—5943(2007)01—0156—05 地球化学图是化探工作中最重要的基础图 坐标、第三列为对应坐标采样点的元素含量) 件,特点是数据多,制图工作量大。如果能用计 如图1。 算机软件自动绘制等含量线色块图,且生成的图 形文件具有属性,将会大大的提高工作效率, 2离散数据网格化 MAPGIS就是这样的软件。 MAPGIS是武汉中地信息工程有限公司研制 在化探工作中水系沉积物测量得到的数据为 的大型基础地理信息系统软件平台,是一个集当 离散数据。微机绘制等值线图最关键的是对原始 代最先进的图形、图像、地质、地理、遥感、测 数据进行网格化(生成GRD数学模型),即采 绘、人工智能、计算机科学于一体的大型智能软 用一定的网格方法对不规则分布的原始数据进行 件系统,是集数字制图、数据库管理及空间分析 插值,生成原始数据分布范围内规则间距的数据 为一体的空间信息系统。 点分布,因此数学模型是绘制等值线图的核心。 MAPGIS离散数据网格化生成GRD模型操作 1原始数据整理 如下: (1)在MAPGIS程序中,进入空间分析下 水系沉积物测量,利用手持GPS地面卫星 的DTM分析即可启动MAPGIS数学模型系统, 定位系统对野外采样点进行定位,得出采样点的 选择离散数据网格化如图2。 坐标(x、Y),通过实验室对化探样品进行分 (2)在离散数据网格化下,本菜单允许用 析,得出各采样点元素含量,由坐标和元素含 户对离散数据进行网格化,系统将弹出如下对话 量,通过Excel生成绘制地球化学图的数据文件 框如图3,该对话框显示了原始数据在x—Y平 .TXT,文件格式为x、Y、z(第一列为绘制 面的范围,用户如果需要扩大或缩小范围,可以 地球化学图的x坐标、第二列为绘制地化图的Y 修改网格化参数中的有关项;通过修改网格间 ①[收稿日期】2007—03—13 [作者简介】袁义生(1965一),女,贵州贵阳人,所学专业:地球化学勘查,工程师,长期从事区域化探工作。 维普资讯 http://www.cqvip.com 第2期 袁义生,等:应用MAPGIS制作地球化学图、单元素异常图及综合异常图 ?157? 图1在Excel中生成制图数据 Fig.1 Data generation in excel 图2对离散数据进行网格 Fi.g2 Meshingfor scatter data 距,用户可以调整网格的疏密程度;在网格类型 栏中,用户有四种不同的网格方法选择,分别为 距离幂函数反比加权网格化、Kring泛克立格法 网格化、稠密数据中值选取网格化、稠密数据高 斯距离权网格化。以上参数设置好后,选择文件 换名按钮系统将弹出标准的文件对话框,输入文 件名。最后按确定,计算机开始对原始数据进行 网格化,并以用户输入的文件名保存网格化后的 结果( .GRD)。 3平面等值线图绘制 选中MAPGIS平面等值线图绘制菜单,打开 网格好的 .GRD文件,系统将弹出平面等值 线图绘制对话框,各参数选择如图4。 在设置中选择等值线套区、绘制色阶、保留 边界线;在光滑中选择等值线光滑处理,光滑 度:高程度;选择轴向标尺;制图注记。 在等值线定层中: (1)等值层值 根据中华人民共和国地质矿产行业标准1: 5万地球化学普查规范(DZ/T0011—91)及1: 20万区域地球化学勘查规范(DZ/T0167 95) 规定,地球化学等含量线的间隔一般可采用 0.1Lg g(或ng/g)间隔取等值线值(见表 1)。等值层值参数的确定,也可采用累积频率 法确定含量线间隔,绘制的地球化学图也能得到 很好的效果,如绘制1:20万贵州省水系沉积物 测量地球化学图,图面既美观,又能表现每个色 标所代表的统计参数含义。 表1 地球化学图等含量线间隔表 Table 1 Interval of Geochemical Isoline 等量线值经过圆滑图上标等量线值经过圆滑图上标 LS g/g 注的的 /g值Lg g/g注的的 g/g值 (2)线参数 在线参数中定义好绘制地球化学图的线型、 线颜色、线宽、线类型(一般用光滑)、等值线 所在图层。 (3)区参数 根据规范,色区值的确定应剔除元素数据特 高值后,计算出数据平均值(x)及标准离差 (s),色区参数确定如表2: 表2地球化学图色区参数表 Table 2 Parameter of Geochemical Mapping Color Zone 维普资讯 http://www.cqvip.com ?158? 贵州地质 2007年第24卷 图3数据网格化选择的参数 Fig.3 Parameter on data meshing choice 元素含X一2S X一0.5 X+o.5 X+2S 量范围<X一2S ~ ~ ~ ~ (p,g/g)X一0.5SX+0.5S X+2S 2(X+2S) 注:X为背景值;S为标准离差 (4)注记参数 在需要标注的等含量线上用Yes,不需要标 注的等含量线上用No。 以上参数设定完成后,按存设置按钮系统将 弹出标准的文件对话框,输入文件名,按确定, 系统将按输入的文件名保存为 .VEL文件。如 果需要重新绘制等值线图时,可以按装入设置按 钮,输入文件名 .VEL文件即可。最后按确定 地球化学图就自动绘制完毕。在文件菜单下 ———存数据于一点数据文件,线数据文 件,区数据文件。 4应用实例 不同的网格方法生成的GRD模型绘制的地 球化学图有一定差异。数据极差相对较小时,采 用Kring泛克立格法网格化模型绘制的地球化学 图效果好,如绘制某区1:5万水系沉积物测量 地球化学图。用同一数据,四种不同网格方法, 生成的GRD模型绘制的地球化学图如图5所示: 相比之下,Kring泛克立格法网格化绘制的地球 图4平面等值线圈参数设置选择 Fi.g4 Selection of parameter setup of plan Isoline graph 化学图线条圆滑,效果较好。元素含量分析数据 极差相对较大时,采用距离幂函数反比加权网格 化模型绘制的地球化学图效果好,如绘制贵州省 地球化学图。用1:20万水系沉积物测量中分析 数据极差相对较大的同一数据,用四种不同网格 方法,生成的GRD模型绘制的地球化学图如图 6所示:相比之下,采用距离幂函数反比加权网 格化模型绘制的地球化学图效果较好。 图5数据极差相对较小时同一数据不同网格方法对比图 Fi.g5 Comparison of diferent meshing way with small data range and same data 1、距离幂函数反比加权网格化 2、Kring泛克立格法网格化 3、稠密数据中值选取网格化 4、稠密数据高斯距离权网格化 剔除元素数据特高值后,计算出数据平均值 (X)、标准离差(S)。 异常下限(N)=平均值+1.65标准离差 维普资讯 http://www.cqvip.com 第2期 袁义生,等:应用MAPGIS制作地球化学图、单元素异常图及综合异常图 ?159? 图6数据极差相对较大时同—数据不同网格方法对比 Fig.6 Comparison ofdiferent meshing way with large data range and same data 确定元素的异常下限N为元素的一级浓度 带,2N为元素的二级浓度带,4N为元素的三级 浓度带。在图形处理的输入编辑下打开绘制好的 地球化学图,在L线编辑菜单——一x参数编辑 ——一根据属性赋参数,系统弹出标准的文件对 话框如图6,选择高度=N,按确定,系统弹出 如图7的文件对话框,选择线型、线颜色、线宽 度、线类型、x系数、Y系数、图层号(异常 层)。以上参数设置完成后,按确定线含量为N 的等含量线即在异常层图层里。用同样的方法把 2N、4N的等含量线放在异常层图层里。在A图 层下关所有层c改层开关N改线L打开异常层; 在G工作区下——一P部份存文件——一存部份 线L得出元素的一、二、三级异常线,将异常 线添加套合在工作区地质图上,绘制好图示图 例,即得单元素异常图。 5 综合常图 根据工作区地球化学特征,确定元素组合, 把元素组合中各元素的单元素异常添加在一起, 套合在工作区地质图上,绘制好图示图例,即得 综合异常图。 圆 豫像磁亨罐辫 }矗 { 口辑件缸忭8 躬.祀瀚嘲旺哪蚺-曩雌|I:疆芑 髓 哟I D叠瘩 鞋甓 杂感 ?“ {l 国 一f{:e 霞鱼¨鼍.f.} 曩 旨瘩傍 簟 i卜 l‘ 图7根据属性提取异常线 Fig.7 Abnormal line abstraction based on attribution 6 结束语 MAPGIS是一个功能强大的制图软件,在数 据极差相对较小时,用Kring泛克立格法网格化 模型绘制的地球化学图效果较好;数据极差相对 较大时,采用距离幂函数反比加权网格化模型绘 制的地球化学图效果佳。两种网格化模型都可以 绘制地球化学图。网格间距的确定一般可以根据 工作区的地质背景,只要绘制的地球化学图能客 观反映元素的地球化学特征即可。 由于目前地质调查项目要求提交的图形文件 为MAPGIS格式,因此在MAPGIS软件下绘制地 球化学图是十分必要的。 本文得到贵州地调院冯济舟、陈光荣两位高 级工程师的指导,在此深表感谢! [参考文献] [1] 中地软件丛书编委会,MAPGIS地理信息系统使用手册 (空间分析篇)[M].武汉中地数码科技有限公司、武 汉中地信息工程有限公司. [2] 中华人民共和国地质矿产行业标准1:57/地球化学普查 规范(DZ/T0011__91)[s].中国地调局. 维普资讯 http://www.cqvip.com ?160? 贵州地质 2007年第24卷 Application of MAPGIS in Creating Geochemical, Single Element Anomaly and Comprehensive Anomaly Maps YUAN Yi—sheng ,LIU Ying—zhong ,LUO Ming—xue ,HE Yan~nan (1 Gu ̄hou Academy ofGeology Surveying,550004 Guiyang, Gu ̄hou,China;2 Geologic 啊105,Guizhou Bureau ofGeology and Minerals Exploration&Development,Guiyang 550018,Guizhou,China) [Abstract] Using MAPGIS meshing model of disperse data for plotting geochemicla map,compared with effects from four types of diferent meshing way,when analytic data range of elementary content is relative small,the effect is good by Kring logarihmitc meshing model for geochemical map.When extreme diference of analytical data is large,the geochemical map is favorable by meshing model of inverse distance power function weighting.Geochemical maps made by two kinds of data both come up to the required standard,therefore it is presented that creates geochemical isoline map wih MAPGIS and pltots single element anomaly and integrated anomaly maps by geochemical isoline. [Key Word] cola quantity;coal core assay;check and naalysis on abnormal data (上接第100页) (4)物探异常呈带状分布,与构造和地 层、岩性套合较好,展现出较大的资源潜力。 因此,龙吟地区具有较好的找矿前景,通过进 —队,20O5. [3] 黄林。赵征,刘金海,等.黔西北五指山背斜矿集区铅 锌矿成矿远景浅析[J].贵州地质,2006.23(3):203 —2o5. 步的地质勘查工作,有望找到中一大型铅锌矿床。 [参考文献] [4] 贵州一O/k地质大队一分队.1:20万盘县幅区域地质调 查报告(下册)[R].贵州省地矿局,1973. [5] 贵州省地质矿产局.贵州省区域矿产志[M].北京:地质出版社,1987. [1] 朱华,杨光龙.贵州省普安县龙吟铅锌银矿普查地质报 告[R].贵州省地矿局106地质大队,1997. [2] 杨光龙,朱华,王智谦,等.贵州省普安县龙吟铅锌矿 预杏地质报告『R1.贵州省地矿局106地质大 Analyses on Geologic Characteristic and Ore——controllng iFactors in the Longyin Lead—Zinc Deposit,Pulan County YANG Guang—long , (1 Guizhou Academy of eolGogy Surveying,550004 Guiyang, Guizhou,China;2 Geologic Party 106,Guizhou Bureau of Geology and Miernals Exploratoin&Developmerit,Zunyi 563000,Guizhou,China) [Abstract]The ore bodies near surface in the Longyin lead—zinc deposit,Pu"an County in Southwest Guizhou Prefecture are characterized by vein and stratoid in occurrence,small ore body but high qualiy,prte— dominately OCCUrring in the space of tensile and tension—torsional faults,as well its derived secondary small faults and interbedding slides,and signiicantly appearfing to be dominated by structure.By analyses on geo— logic characteristics and ore—controlling factors in the Longyin lead—zinc deposit,this area is considered a favorable prospect potential。and is expected to find out middle—large—size lead—zinc deposit by in depth geologic works. [Key Word]lead—zinc deposit;geologic characteirstics;ore—controlling factos;prrospects potentil;a I ̄ngyin mining area;Puhn;Guizhou
应用MAPGIS制作地球化学图单元素异常图及综合异常图
应用MAPGIS制作地球化学图单元素异常图及综合异常图 .doc
将本文的Word文档下载到电脑