基于栅格数据模型的局部奇异性分析迭代算法.pdfVIP

基于栅格数据模型的局部奇异性分析迭代算法 陈志军。2成秋明“¨陈建国“2 1．中国地质大学(武汉)地质过程与矿产资源国家重点实验室，湖北武汉430074 2．中国地质大学(武汉)资源学院，湖北武汉430074 3．York大学地球空问科学系和地理系，加拿大多伦多M3JIP3 摘要：地球化学样品往往受到古代和现代众多复杂的相-百作用过程的影响，需对采集样品 的统计分布律、空间相关及变异性、异常几何特征(丛聚形状、方向性)、模式的尺度独立 性等进行多方面的考察。基于多重分形理论的局部奇异性分析为弱缓异常的识别提供了一种 崭新的研究思路和方法技术，受到国内外同行的广泛关注。本文对奇异性分析的基本思想进 行了介绍，局部奇异性指数作为指示异常的关键指标，如何减少其估算偏差在实际应用中尤 为重要。方形窗口滑动平均方法是目前流行的最简单化的算法，本文在剖析常规算法的基础 上提出了一种迭代算法来获取最终的局部奇异性指数及其相伴的局部系数，该系数具有非奇 异的特点。 关键词：局部奇异性分析；栅格数据模型；迭代算法；多重分形； 近年来，国内外许多地质学家从不同角度认识到地球化学数据场具有多霞分形的特点 (Cheng，etal，1994；李长江等，1999；成秋明，2000)。地球化学样品往往受到古代和现代众多 复杂的相互作用过程的影响，元素地球化学含餐的空间分布具有极其复杂的不规则性，我们 需对采集样品的统计分布律、空间相关及变异性、异常几何特征(丛聚形状、方向性)、模 式的尺度独立性等进行多方面的考察。才能更好地发现异常，特别在高背景巾提取矿致异常、 在低背景中发现弱缓异常。研究表明，成矿异常范围内元素的浓度分布与背景区的分布服从 不同的统计分布：前者服从分形分布(Preato分布)，直方图一卜望明显的正向偏斜甚至很长的 et 拖尾，而后者服从对数正态分布(Cheng 这十年间的主要进展在GEEA杂志进行了评述，他指出：模式识别和异常检测(Pattem and detecti∞)的数据分析方法技术取得了重要进展，多种形式的分形和 recognitionanomaly et 相关分析方法(Chenga1．1994；ChengAgterberg1996；Gon9alves200l：ShenCohen et 2005；Li a1．2004)对单变量或多变量的多总体分离得到了极大关注。我国学者成秋明提出 的C．A(含量．面积模型)。S．A(能谱密度．面积模型)等最新的分析和多重分形方法己被荷 andMineral in 兰学者Carranza(2009)编著的(GeochemicalAnomalY ProspectivityMapping 在对该书的书评中建议还应该引用另一个强有力的新的水系沉积物和其他地球化学数据处 理的方法一局部奇异性分析方法。 l局部奇异性分析方法的基本思想 奇异性理论指出：奇异性是指在很小的时间．空间范围内具有巨大能量释放或巨量物质 分形(如矿物的富集)，其分维数为D，工尺度上质点集合的质量肘弛)与集合中质点的数目 Ⅳ弛)成正比，该分形体的密度为： ‘基金项目：国家自然科学基金青年科学基金项目(No；国家杰出青年科学研究基金资助项Iil(No． 中国地质调查局工作项目(资[20wW评ol·18-01)：中央高校基本科研业务费专项资金资助项II资助；地 质过程与矿产资源国家重点实验室科技部专项经费资助。 作者简介：陈志军(1978-)．男，江苏常熟人。副教授，主要从事数学地质的教学和研究。E-mail： ‘ p：罂芘掣芘矿F E r (1) 当DE时，物体体积增长要比它的质量增长快，随着三的增大，密度会逐渐变小，物 体也会变得越来越稀疏：反之从三的刁i断减少方向考察，物体的密度将越来越高直至无穷， 呈现强烈的奇异性。模型(1)是一个简单的分形模型，只需一个分维数就能完全确定密度， 而对于复杂情况，需从多重分形分布来考察，用多个参数才能完整地描述这个富集过程的密 度分布。