基于数据变化分析的重力局部异常识别方法.pdfVIP

·492· 中国地球物理2011 (15)信息技术与地球物理 基于数据变化分析的重力局部异常识别方法 Methodsof local identificationbasedon gravityanomaly ofdata analysis changes 谢永茂姚长利张晨郑元满 Xie Yao Chen Yuanman YongmaoChangliZhang Zheng 地下信息探测技术与仪器教育部重点实验室(中国地质大学·北京) 北京100083 1．前言 通常，地球物理场是地下不同深度、不同尺度的各种地质体所引起的异常总和，其中局部地质体形成 的异常叠加在区域背景场上。在地球物理勘探作业，为了能够通过反演计算得到局部地质体的物性及几何 形态，达到找矿的目的，通常需要从实测数据中对局部异常进行识别与提取，然后再进行进一步的处理计 算及分析。但是，局部异常与区域背景场的叠加使局部异常难以辨认，通常依赖于解释人员的专业水平和 主观认识。本文基于滑动窗口的数据变化程度，和图形图像处理与识别相关技术思想，对重力局部异常的 识别和提取进行了相应的方法技术研究，提出了一种简单且相对有效的重力局部异常识别方法，形成自动 化的识别算法，为大范围数据异常的自动识别打下基础。 2．研究方法及技术要点 相对于区域背景场而言，局部异常区域，特别是局部异常的边界位置，其数据变化的程度较大。因此， 首先要通过建立对数据变化程度的评判方法，对局部异常区域进行初步的识别；然后利用图形图像处理与 识别相关技术，对初步识别的结果进行再识别与提取。其中主要的方法步骤与技术要点如下： 1)初步识别局部异常区域 基于选定的的滑动窗口内数据的均方差衡量数据的变化程度，当变化程度大于设置的阈值时，认 为该点数据位于局部异常区域内。为了能够反映整个数据区域的相对变化，该阈值设置为：舻∥，其中 ∥为所有点的窗口均方差的均值，A为系数。窗口较小，能较好地识别出局部异常的边界，但是位于局 部异常中间的平缓区域需进～步识别。 2)识别局部异常的平缓区域 第一步是利用数字图像处理中的贴标签法区分初步识别的各个局部异常，即对各个异常贴上连续的标 签号，不同标签号代表不同异常，为后续对各个局部异常的平缓区域进行识别，以及可能需要的对其特征 差异进行描述与反演计算并获得所需要的结果等工作打下基础。 第二步是基于计算机图形学中的射线法思想，以及对于一些特例情况加以针对性分析与处理，针对每 个局部异常，识别其边界区域内的平缓区域。 3)提取已识别的局部异常 先从原始数据中将已识别的局部异常区域剥离，再对余下的数据重新进行网格化处理得到区域背景 场，最后将原始数据减去区域背景场，即可得到需要的局部异常。 3．结语与建议 通过多个模型数据计算检验，证明了该方法技术的有效性，以及图形图像处理与识别技术在识别局部 异常中具有一定的价值和可行性。但方法技术仍需进一步深化、完善。仅通过均方差这单一的评判尺度， 对实际资料的复杂组合异常，还达不到理想的识别结果。需要进一步挖掘更多可用的数据信息，如区域内 的均值、斜率等，以及局部的趋势等，将这些参数进行合理组合，从而形成一套多参数、多尺度的分析异 常、识别异常的方法技术，为重力异常的处理解释服务。 (1212010660502)的支持。 参考文献 【l】程方道，刘东甲，姚汝信．划分重力区域场与局部场的研究【J】．物探化探计算技术，1987，9(1)：I~9． 【2】杨文采，施志群，候遵泽，等．离散小波变换与重力异常多重分解叨．地球物理学报，2001，44(4)：534--541 【3】阮秋绮．数字图像处理学【M】．北京：电子工业出版社，2000．