设计测点在地形图上完成，在1：50000地形图上按0.25km2划出的采样格子（称小格）中的绝大多数（95%以上）小格内都有采样点分布，应做到基本不出现或很少出现连续3个以上空白小格。当采样密度平均为4-5个点/km2时，小格内样品数一般不要多于2个。采样点分布均匀并不意味着将采样点布置在采样格子中央，而应将采样点布置在每一个格子中能最大限度控制汇水面积处。样品采集按照设计坐标采用手持GPS严格按照相关规范和确定的方案进行采样，最终完成水系沉积物样品采集247件。
2.2 预处理
对采集的样品送往国土资源部昆明矿产资源监督检验中心分析，分析项目包括Au、Ag、Cu、Pb、Zn、W、Sn、Bi、Mo、As、Sb、Hg共12种元素（表1）。将野外手持GPS获得的采样点的坐标与样品12种元素测试分析得到的结果一一对应整理到一个文本文件（*.txt）中。
通过MAPGIS 6.7主菜单→实用服务→投影变换→P投影转换→用户投影转换打开第一步骤的文本文件，设置X（坐标）位于与Y（坐标）位于文件中所在的列。设置分隔符。勾选分隔符号中的Tab键，设置属性名称所在列，同时设置数据类型（一般为双精度），这样就可以在生成的图元文件中带有属性信息。设置输入投影参数与输出投影参数，需要设置到当地的投影中心点经度，同时设置好点图元参数，点击投影变换即完成了含有属性信息点文件（化探.WT）的制作。
3 统计分析和异常图件的制作
3.1 统计分析
3.1.1 数字特征统计
将前述经过投影变换的点文件（化探.WT）添加到数字填图的工程系统中，点击综合数据处理→数字特征统计→数字特征统计完成对老挝波乔水系沉积物12种元素分析测试结果的统计，按照对数平均值+2倍数标准离差来确定异常下限，同时可以得到12中元素的数字特征统计直方图及变差系数。
测试的12种元素分布显示，仅As、Au、Bi、Cu、Pb含量分布不均匀，属强分异型，Sb、Zn为中等分异常型，其余5种元素为弱分异或者无明显分异类型。即从地球化学观点看，矿区多金属矿成矿具有一定潜力，可以进一步开展工作。
3.1.2 聚类分析[14]
利用数字填图系统中综合数据处理→多元统计→聚类分析可以实现化探数据处理中的聚类分析谱系图制作。
由聚类分析谱系图（图1）看出，相关系数>0.5时，分为Ag、Sn、Sb和As、Cu、Zn、Bi、W、Mo及Au、Pb三个组合。反映该异常存在As、Cu、Zn、Bi、W、Mo中高温热液成矿作用和Ag、Sn、Sb与Au、Pb的低温热液成矿过程，从该多金属矿区的岩性条件来看，具备寻找Cu、Pb、Zn、Au的成矿条件。
3.1.3 因子分析
利用数字填图系统中综合数据处理→多元统计→因子分析可以实现化探数据处理中正交因子解的获得（表2）。
从表2可以看出，Cu、Pb、Zn三种元素主要受第1因子的影响，Ag主要受第2因子的影响，Au主要受第3因子的影响，Hg主要受第4因子的影响。
从因子分析的结果来看，矿区寻找Cu、Pb、Zn、Au是具备潜力的，同时再从正交因子计量表格中对应第1因子，第3因子可以找出最具备找矿潜力的坐标区域。
3.2 地球化学图件制作
3.2.1 点位数据图制作
将该点文件添加到MAPGIS编辑子系统中将点文件处于输入状态，点击N点编辑→根据属性标注释既可以完成12种元素的原始数据图制作。
3.2.2 单元素异常图制作
利用MAPGIS中的DTM分析或者数字填图软件中的等值线图均可以实现单元素异常图的绘制，级别以统计的异常下限为依据，按照1、2、4、8定级，用颜色的渐变反映异常的外、中、内带（图2）。
3.2.3 组合元素异常图制作
依据其元素地球化学场分布形态、位置及亲合性，及聚类分析、因子分析成果将12种元素分为3种组，①Ag-Sn-Sb-Hg组合；②As-Cu-Zn-Bi-W-Mo-Hg组合；③Au-Pb-Hg组合，同时为了更好的分析，制作了12种元素组合异常图。
从综合异常图来分析，可以圈定矿区东北部Cu、Pb、Zn成矿远景区，矿区中部及西部Au成矿远景区。
4 结论
MAPGIS与数字填图软件都具有一定的空间数据处理及成图的功能，特别是数字填图软件具备较强的空间数据分析能力，这也是越来越多的地学工作者开始喜欢使用这两个软件完成地球化学数据处理的最大优点。
笔者采用数字填图软件对老挝波乔的水系沉积物采样的测试结果进行了特征值统计，得到了相关元素的异常下限及变差系数，为开展进一步的化探数据处理打下了基础。同时开展了聚类分析，得到了相关元素的谱系图，为组合元素异常图绘制及成矿远景区预测奠定了基础。开展了因子分析，也为开展成矿远景区预测奠定了基础。
在地球化学图件的绘制方面采用MAPGIS完成，形成了点位数据图、单元素异常图、组合元素异常图，为地球化学异常区的评价及解释提供了基础资料，同时也为下一步开展进一步的地质工作打下了基础。
参考文献：
[1]蔡朝阳，孙德有，吴国学.mapgis与surfer相结合在黑河市化探工作中的应用[J].物探与化探，2009，33（2）：220-223.
[2]陈志军，陈建国.基于MapGIS的图切地质剖面系统面向对象设计与实现[J].地质学刊，2012，36（3）：256-264.
[3]汪新庆，冯磊.基于数据字典的MapGIS属性数据逻辑检查——以全国矿产资源潜力评价基础成果数据检查为例[J].地质学刊，2013，37（3）：478-481.
[4]胡雄伟，吴良士.老挝人民民主共和国地质特征与区域成矿[J].矿床地质，2009，28（1）：104-106.
[5]杨加庆，李强，左琼华，等.老挝布劳铁矿地质特征及成因探讨[J].云南地质，2013，32（2）：160-163.
[6]陈载林，左琼华，王建飞.老挝波乔省某铁矿区磁法勘探[J].地球物理学进展，2013，28（3）：1447-1452.
[7]孙建中，程之牧，姜志祥，等.老挝会晒蓝宝石资源评价[J].上海地质，1998，2：59-62.
[8]莫运东，孙明祥，王国洪，等.景洪南林一大勐龙地区铁多金属矿资源潜力[J].云南地质，2009，28（3）：59-62.
[9]赵含军，范建辉，廖锦川.南勐腊新山铁多金属矿矿化富集规律[J].云南地质，2011，30（2）：157-159.
[10]王爱萍，高李阳，张吉寿.云南勐海关双铁矿矿床成因[J].云南地质，2011，30（2）：160-163.
[11]罗先熔.勘查地球化学[M].北京：冶金工业出版社，2007.
[12]杨小峰.地球化学找矿方法[M].北京：地质出版社，2007.
[13]王进进，陈玉梁，李映海，等.云南维西山加铜矿床地质特征及成因浅析[J].价值工程，2016，35（8）：160-163.
[14]纪宏金，孙丰月，曾道明，等.对应聚类分析与盲矿体的地球化学预测[J].物探化探计算技术，2007，29（2）：138-142.
