MAPGIS在豫西灰池子岩体地面伽玛总量测量中的应用

郭长江　王　娟

（河南省核工业地质局地调院，河南　信阳　464000）



MAPGIS在豫西灰池子岩体地面伽玛总量测量中的应用

郭长江王娟

（河南省核工业地质局地调院，河南信阳464000）

摘要：本文以《河南省卢氏县灰池子岩体外围地面伽玛总量测量》项目1∶5万地面伽玛总量测量为例，从原始数据采集、室内资料整理、图件绘制等方面，简述伽玛测量数据处理的方法和步骤，重点对MAPGIS软件在伽玛总量测量图件绘制中的应用及各程序之间的转换进行阐述。

关键词：MAPGIS；伽玛总量测量；数据处理；DTM分析

1　项目概况及原始数据采集

1.1项目概况

《河南省卢氏县灰池子岩体外围地面伽玛总量测量》项目区位于河南省卢氏县西南部，位于豫陕两省三县交界处，与陕西省丹凤、商南两县接壤。主要任务是：通过1∶5万地面伽玛总量测量工作，探明该区内放射性高场、偏高场、异常场的分布特征；系统整理项目区内异常点、带的资料，圈定找矿远景区，为豫西成矿带远景调查工作提供依据[1，2]。

1.2数据的采集

在开展工作区地面伽玛总量测量工作时，测量人员、地质人员在全面了解测区的地质、构造、岩石、蚀变和矿化等情况，野外均采用GPS定点，测量仪器为FD-3013型 γ辐射仪。根据项目设计，由南部向北部共设测量基线4条，依次编号为A、B、C、D，测线编号依基线号依次进行编号。共完成基本测量线178条，总线长488.64km，共完成设计4 409个测点的测量任务。

2　利用MAPGIS制作伽玛总量等值图

MAPGIS数字高程模型DTM是利用一个坐标系中大量采集的已知X、Y、Z的坐标点，对连续地面的一种表示地形特征空间分布的数字定量模型，是带有空间位置特征和地形属性特征的数字描述。X、Y表示该点的平面坐标，Z值可以表示高程、坡度、温度等变量信息。在伽玛总量测量中，Z表示伽玛测值。

2.1原始数据的处理

在进行实际的野外勘测过程中，对于伽玛总量测量一般需要通过GPS技术对测量点进行定位，获得相应的采样坐标X、Y及伽玛测值，并将这些原始的数据导入到数据处理表格中进行整理，形成名称为“伽玛总量测量结果”的原始数据记录，然后将Excel格式另存为“.csv”格式。将“.csv”文件扩展名直接更改为“.txt”，用操作系统自带的“记事本”打开文件。然后将第一行（即X、Y坐标和测值符号行）改写为“notgrid”，再将“.txt”扩展名更改为“.det”文件保存，留作DTM分析时使用。

表1　伽玛场的分级标准

2.2背景值和异常下限值的确定

异常下限值的确定通常需要将高值点的数据进行剔除后，根据伽玛测值的平均值与标准值之间的差距来确定。在实际的统计计算工作中，需要利用Excel表格中的计算函数，将伽玛测值按照高低顺序进行排列，然后根据伽玛测值来计算出平均含量，将计算结果比标准含量大3倍的数据进行剔除，再次利用平均值和标准值的计算方法进行计算，直到将所有含量值较大的数据剔除之后，计算出本项目区岩石的平均值和均方差。根据《EJ/T831-94地面伽玛总量测量规范》中伽玛分级标准计算出工作区各伽玛异常场圈定值（见表1）。

2.3数据投影变换

数据的投影变换过程的完成需要利用MAPGIS6.7的投影变换功能中“用户文件投影转换”来实现，具体工作步骤：投影变换-用户数据点投影变换-打开用户文件“XX地区伽玛总量测量结果.txt”（只能是纯文本格式）-设置用户文件选项（选择“按指定分隔符”）-设置分隔符-指定数据起始位置-设置用户投影参数-结果投影参数-设置用户文件选项-设置生成点图元参数-投影转换-确定。如果需要将投影结果写到文本文件，那么点“写到文件”按系统提示输入投影结果文件名，输入完毕即进行转换。需要注意的是，保存的文件没有文件类型，需要手动输入“.txt”文件类型，才可以使用记事本打开。以上步骤也可以利用MAP⁃GIS二次开发软件Section来实现，过程更为简单和直观。

2.4离散网格数据化

在地面伽玛总量测量成图过程中，如何保证生成的地质图清晰、准确，对原始数据进行离散和网络化处理，是非常关键的一步。MAPGIS中有2种方法可以生成伽玛总量异常图。打开MAPGIS6.7空间分析中DTM分析模块，具体操作步骤分别如下。

2.4.1GRD模型。空间分析-DTM分析-打开三角剖分文件（伽玛总量测量结果.det）。离散数据网格化-生成Grd文件-平面等值线图绘制（设置等值线、等值层值、注记等参数）-保存设置-确认-在“文件”菜单下另存保存WT、WL、WP文件。

2.4.2TIN模型。空间分析-DTM分析-打开点数据文件-点数据高程提取-快速生成三角剖分网-整理三角剖分网-追踪剖分等值线（设置等值线、等值层值、注记等参数）保存设置-确认-在“文件”菜单下另存保存WT、WL、WP文件。需要说明的是，生成的等值线图，个别地方还存在比较生硬、圈定不合理的问题，需做人为修改，使之合理、美观。实际操作中还发现MAPGIS6.7无法使用DTM分析模块现象，那么可以用MAPGIS6.6来制作。

在追踪剖分等值线时，关键的是等值线参数设置，对不需要表示的等值层采取删除一层的命令来选取。对特殊的等值层采取添加一层的方法来增加。等值线光滑处理采用中精度处理，这样经追踪剖分的等值线才会比较光滑，接近真实情况。

2.5建立图件工程文件

主要步骤如下：MAPGIS主菜单-“图形处理”-“输入编辑”，新建工程，打开MAPGIS编辑子系统。选择“不生成可编辑项”，在“文件名”栏目下方空白处，右击鼠标，添加项目前面保存的伽玛总量测量WT、WL、WP文件添加到系统中，复位窗口后就可以看到添加的图形内容。文件按照WP、WL、WT排序，点击刷新按钮。还可以按照需要添加其他已编辑好的地质、地理资料图层，通过编辑有关WT、WL、WP文件，添加文字说明、比例尺、图例、图框等工作，就可以制作出一张符合规范的地面伽玛测量等值图。

总之，MAPGIS系统在伽玛总量测量图件绘制中的应用效果良好。需要注意的是，MAPGIS6.7生成的伽玛等值图存在圈定不合理的现象，部分线条显得较为生硬，这就要求地质技术人员通过地质特征进行再次编辑和加工，才能获得准确且美观的伽玛等值图。

参考文献：

［1］武汉中地信息工程有限公司.MAPGIS地理信息系统使用手册-数字制图篇［M］.北京：北京中地时代软件工程有限公司，2003.

［2］王旭东.浅议MAPGIS地质测绘数字化应用技术［J］.城市建设理论研究，2011（20）:17-18.

中图分类号：P631.6

文献标识码：A

文章编号：1003-5168（2016）01-0063-02

收稿日期：2015-12-21

作者简介：郭长江（1977-），男，本科，助理工程师，研究方向：地质测绘。

The Application of MAPGIS in Yuxi Huichizi Gamma Measurement of the Total Mass

Guo ChangjiangWang Juan
（Henan Province Nuclear Industry Geological Bureau，Xinyang Henan 464000）

Abstract:This paper taked the henan province lushi county of huichizi peripheral rock ground gamma total measure⁃ment"project 1:50 000 ground gamma total measurement as an example,the collection of raw data,and indoor data processing and plotting were described based on the gamma measurement data processing methods and steps,focus⁃ing on the transitions between the application of MAPGIS software in the gamma total measurement plotting and the program were described.

Keywords:MAPGIS;gamma ray measurement;data processing;DTM analysis