基于磁法考古数据数据处理系统MATLAB实现.pdfVIP

基于磁法考古数据的数据处理系统的MATLAB 实现 王春庆，陈辉，杨进 中国地质大学（北京）地球物理与信息技术学院，北京（100083 ） 中国地质大学（北京）地下信息探测技术与仪器教育部重点实验室，北京（100083 ） E-mail ：chunqingicy732@ 摘 要：磁法用于考古探测是地球物理考古中的一种常用方法。本文基于磁法考古的数据， 采用MATLAB 语言，自行编制了一套数据处理系统，并在实际考古数据中验证了该系统的 可靠性和实用性。 关键词：磁法，考古，数据处理系统，MATLAB 1. 引 言 在考古勘查中，数十年来，应用地球物理勘探方法(简称物探) 由于受到各种原因的影响， 始终没有得到十分广泛应用，虽然也取得一些成绩，但总不尽人意，究其原因，除考古与物 探人员的交流不够、考古经费有限外，物探勘查效果及精度、勘察方法技术、数据处理软件 是主要影响考古物探顺利开展的因素。我国是一个文明古国，地下埋藏文物众多，尚未探查 的地下遗迹等无计其数，对我们来说是一笔巨大的财富，更具有很高的科学研究价值。因此， 我国今后很长的一段时间内，考古的任务将是很重的。研究、改进探测古遗址的地球物理探 测技术，对我国具有特殊的重要意义；研究开发简便，快捷，有效的数据处理软件，在物探 界具有举足轻重的地位。 磁测技术是一种有效的地球物理勘探方法，用于考古研究开始于 1958 年[1] 。当时主要 是用来确定古代炼铁炉和砖瓦窑的具体位置，以及通过探测古墓葬的石头地基来寻找古墓。 由于该方法具有速度快、成本低、效果好和无损伤的优点，几十年来一直是地球物理考古中 的一种常用方法，并越来越受到世界各国的重视。特别是由于近年来各类高灵敏度磁力仪和 磁梯度仪的出现，则使高精度磁测技术更广泛地用于考古调查研究，为进一步探测磁性弱， 埋藏深的古墓葬、古遗迹提供了得力的工具。高精度磁力测量广泛应用于固体矿产勘查、地 质填图、石油资源勘查及工程勘查等领域。磁法用于考古探测主要表现在两方面[2 3 4] ：一是 直接探测或填绘出地下埋藏的古物（如被掩埋的石质古建筑地基、古墓、陶器等等），这些 被探测对象的磁性一般不同于周围的背景；二是通过上部土壤的磁化率测量来寻找古人类活 动的遗址。 埋藏在地下的古代砖、瓦、陶器、窑、炉、墓葬、城墙、道路和居民点等古文物和古遗 迹，都具有一定的磁性，并与其周围的物质形成了一定的磁性差异。如果这些磁性差异明显， 并在磁力仪灵敏度所能分辨的范围内，则高精度磁力仪在地面所观测到的磁异常就成为人们 进行考古研究的重要依据。另外，磁法较地质雷达法和高密度电法轻便、快速、廉价和具有 干扰易识别的优势。 为了提高由原始数据到成图的工作效率，避免反复运算带来的麻烦，特为磁法考古数据 运算和处理过程设计了数据处理系统。该系统使用计算机MATLAB语言[5]，包括数据的处理、 网度分析、图件显示等几方面的内容。 2. 系统概述 本文借助在陕西某古遗址采用高精度磁法进行考古探测的数据，自行编制了一套数据处 - 1 - 理系统，并在实际应用中验证了其可靠性和实用性。 本处理系统图形用户界面采用MATLAB 语言开发也是一个尝试，大多搞计算机开发的 同行都是采用VC ，VB 等语言，但 MATLAB 语言也不次之。被称作为第四代计算机语言 的MATLAB ，利用其丰富的函数资源，使编程人员从繁琐的程序代码中解放出来。MATLAB 给用户带来的是最直观，最简洁的程序开发环境，符合人们思维习惯。除具备卓越的数值计 算能力外，MATLAB 的图形功能强大，具有较强的编辑图形界面的能力。 数据为采用高精度磁法野外采集原始数据进行日变改正后的数据，输入文件格式主要针 对点号、线号、坐标、磁场改正值形式的原*.txt 格式文件，经处理后输出格式文件是*.grd 和*.txt 格式文件。 本处理系统中的图形用户主界面根据用途主要分为数据处理、磁测图件和网度分析三部 分，如图1。“数据处理”包括五个子界面，分别为网格化、圆滑、噪声