均值归一总水平导数边界识别方法.docx

均值归一总水平导数边界识别方法王赛昕，刘林静（中国地质大学地球物理与空间信息学院，湖北武汉４３００７４）摘要：分析了近年利用重磁异常进行边界识别的方法，指出Ｔｉｌｔ－ａｎｇｌｅ法、ＴｈｅｔａＭａｐ法等一类方法是通过计算重磁异常的高阶导数及其比值来增强地质边界信息，由于高通滤波对数据误差放大，在一个计算点上两种导数的比值可能会有比较大的波动，导致实际资料计算结果无法利用。为了避免在一个点计算比值造成波动的影响，笔者采用总水平导数均值归一的方法，避免了实际资料处理结果的振荡；同时，对不同强度的总水平梯度异常采用不同的均值进行归一，突出了弱地质边界信息。将该方法应用于大杉树钒钛磁铁矿远景区对辉长岩的边界识别，得到较好的地质效果。关键词：重磁异常；边界识别；归一化；总水平导数；辉长岩中图分类号：Ｐ６３１文献标识码：Ａ收稿日期：２０１１－０９－１５ＥｄｇｅＤｅｔｅｃｔｉｏｎｂｙＥｑｕａｌｉｚｅｄａｎｄＮｏｒｍａｌｉｚｅｄＡｍｐｌｉｔｕｄｅｏｆＴｏｔａｌＨｏｒｉｚｏｎｔａｌＤｅｒｉｖａｔｉｖｅｓＷａｎｇＳａｉｘｉｎ，ＬｉｕＬｉｎｊｉｎｇ（ＩｎｓｔｉｔｕｔｅｏｆＧｅｏｐｈｙｓｉｃｓ＆Ｇｅｏｍａｔｉｃｓ，ＣｈｉｎａＵｎｉｖｅｒｓｉｔｙｏｆＧｅｏｓｃｉｅｎｃｅｓ，ＷｕｈａｎＨｕｂｅｉ４３００７４，Ｃｈｉｎａ）Ａｂｓｔｒａｃｔ：Ｉｎｔｈｉｓｐａｐｅｒ，ｗｅａｎａｌｙｚｅｓｏｍｅｍｅｔｈｏｄｓｆｏｒｅｄｇｅｄｅｔｅｃｔｉｏｎｂａｓｅｄｏｎｇｒａｖｉｔｙａｎｄｍａｇｎｅｔｉｃａｎｏｍａｌｙｉｎｒｅｃｅｎｔｙｅａｒｓ．ＷｅｐｏｉｎｔｅｄｏｕｔｔｈａｔｔｈｅｍｅｔｈｏｄｓｓｕｃｈａｓＴｉｌｔａｎｇｌｅａｎｄＴｈｅｔａＭａｐｅｎｒｉｃｈｔｈｅｓｏｕｒｃｅｅｄｇｅｂｙｃａｌｃｕｌａｔｉｎｇｔｈｅｈｉｇｈｏｒｄｅｒｄｅｒｉｖａｔｉｖｅｓｏｆｇｒａｖｉｔｙｏｒｍａｇｎｅｔｉｃａｎｏｍａｌｙａｎｄｔｈｅｒａｔｉｏｏｆｔｈｅｍ．Ｂｅｃａｕｓｅｔｈｅｈｉｇｈｐａｓｓｆｉｌｔｅｒｉｎｃｒｅａｓｅｓｄａｔａｅｒｒｏｒ，ｔｈｅｒａｔｉｏｏｆｔｈｅｔｗｏｄｅｒｉｖａｔｉｖｅｓｍａｙｆｌｕｃｔｕａｔｅｓｔｒｏｎｇｌｙｏｎａｄｏｔａｎｄｔｈｅａｃｔｕａｌｄａｔａｃａｎｎｏｔｂｅｕｓｅｄ．Ｔｏａｖｏｉｄｔｈｅｅｆｆｅｃｔｏｆｔｈｅｒａｔｉｏｆｌｕｃｔｕａｔｉｏｎｏｎａｄｏｔ，ｗｅａｄｏｐｔｔｈｅｍｅｔｈｏｄｎａｍｅｄｅ－ｑｕａｌｉｚｅｄａｎｄｎｏｒｍａｌｉｚｅｄａｍｐｌｉｔｕｄｅｏｆｔｏｔａｌｈｏｒｉｚｏｎｔａｌｄｅｒｉｖａｔｉｖｅｓ，ｗｈｉｃｈａｖｏｉｄｔｈｅｆｌｕｃｔｕａ－ｔｉｏｎｏｆｔｈｅｃａｌｃｕｌａｔｅｄｒｅｓｕｌｔｏｆｔｈｅａｃｔｕａｌｄａｔａ．Ｍｏｒｅｏｖｅｒ，ｉｔｕｓｅｓｄｉｆｆｅｒｅｎｔａｖｅｒａｇｅｔｏｎｏｒ－ｍａｌｉｚｅｔｏｔａｌｈｏｒｉｚｏｎｔａｌｄｅｒｉｖａｔｉｖｅｓｔｏｗａｒｄｓｄｉｆｆｅｒｅｎｔｉｎｔｅｎｓｉｔｙ．Ａｓａｒｅｓｕｌｔ，ｉｔｅｎｈａｎｃｅｓｔｈｅｉｎｆｏｒｍａｔｉｏｎｏｆｆａｉｎｔｓｏｕｒｃｅｅｄｇｅ．ＡｐｐｌｙｉｎｇｔｈｉｓｍｅｔｈｏｄｔｏｄｅｔｅｃｔｉｎｇｔｈｅｇａｂｂｒｏｅｄｇｅｉｎＤａｓ－ｈａｎｓｈｕｖａｎａｄｉｕｍａｎｄｔｉｔａｎｉｕｍｍａｇｎｅｔｉｔｅｌｏｎｇｒａｎｇｅｒｅｇｉｏｎ，ａｇｏｏｄｇｅｏｌｏｇｉｃａｌｅｆｆｅｃｔｉｓａ－ｃｈｉｅｖｅｄ．基金项目：全国危机矿山接替资源找矿项目（编号：２００７９９０８４）；国土资源部公益性行业科研专项（编号：２００９１１０１７－０３）；中国地质大学（武汉）第八届“挑战杯”大学生课外学术作品竞赛项目资助。工程地球物理学报（ＣｈｉｎｅｓｅＪｏｕｒｎａｌｏｆＥｎｇｉｎｅｅｒｉｎｇＧｅｏｐｈｙｓｉｃｓ）第８卷７００Ｋｅｙｗｏｒｄｓ：ｇｒａｖｉｔｙａｎｄｍａｇｎｅｔｉｃａｎｏｍａｌｙ；ｅｄｇｅｄｅｔｅｃｔｉｏｎ；ｎｏｒｍａｌｉｚｅｄａｍｐｌｉｔｕｄｅ；ｔｏｔａｌｈｏｒ－ｉｚｏｎｔａｌｄｅｒｉｖａｔｉｖｅｓ；ｇａｂｂｒｏ平导数边界识别方法，其基本思想是对Ｔｉｌｔ－ａｎ－引言１ｇｌｅ再求总水平导数［１２］。Ｗｉｊｎｓ等在２００５年提出θ图法，该方法利用解析信号振幅对总水平导数进行归一化，利用极大值位置确定地质体的边界位置［１３］。由于θ图是基于导数的比值而来，所以θ图能很好的平衡高幅值异常和低幅值异常，起到边界增强的效果。分析近年发展起来的这些方法，它们都是计算重磁异常的高阶导数及其比值来增强地质边界信息的。人们知道，高阶导数是一种高通滤波，它对数据误差有很强的放大作用，在一个计算点上两种导数的比值可能会有比较大的波动，导致实际资料计算结果无法利用。为了避免在一个点计算比值造成波动的影响，人们采用总水平导数均值归一的方法，避免了实际资料处理结果的振荡；同时，由于对不同强度的总水平梯度异常采用不同的均