基于栅格数据的局部奇异性分析迭代方法及其软件实现.pdfVIP

613—622．2014 地 质 学 刊 第 38卷 第 4期 doi：10．3969／j．issn．1674—3636．2014．04．613 基于栅格数据的局部奇异性分析迭代方法及其软件实现 陈志军 ，成秋明 ， ，陈建国 ， (1．中国地质大学地质过程与矿产资源国家重点实验室，湖北 武汉 430074；2．中国地质大学 (武 汉)资源学 院，湖北 武汉 430074；3．York大学地球空 间科学系和地理系，加拿大 多伦多 M3J1P3) 摘要：近些年，弱缓化探异常识别已成为成矿预测和勘查评价中十分关键的问题。多重分形理论 的局部奇异性分 析因其崭新的思路、简便的方法、优 良的效果而在弱缓异常识别中受到广泛关注。在深入剖析局部奇异性分析基 本思想及计算方法的基础上，引入正规化尺度参数 ，提出了迭代方法来改进局部奇异性指数的估值。给出了奇异 性迭代算法并用 c++编程实现，软件功能强劲，操作灵活，不仅适用于各向同性情形，还适用于各 向异性尺度的奇 异性计算和方向性奇异性计算。软件的动态链接库版本已在 GeoDAS矿产资源定量预测专业软件中应用。 关键词：局部奇异性分析；栅格数据；多重分形 ；迭代方法；正规化尺度参数 中图分类号：P628 文献标识码：A 文章编号：1674-3636(2014)04-0613-10 型，试图通过局部异常的有效识别来达到空间预测 0 引 言 目的，该模型正逐渐在覆盖区勘查地球化学弱缓异 常识别中表现出巨大的应用潜力，成为非线性矿产 分形与多重分形在勘查地球化学数据、勘查地 预测的有效方法之一(成秋明，2012)。应用的需求 球物理数据、矿石品位数据等地质勘查数据 中的普 也促进了局部奇异性分析 自身模型的发展，出现了 遍存在性为成矿过程奇异性与矿产预测定量化的新 多种形式的扩展模型(陈志军，2007)。主要从提高 理论与新方法研究提供了崭新的思路和新型的数学 奇异性指数计算精度方面介绍局部奇异性分析迭代 工具(成秋明，2007；Agterberg，2003)。不均匀性、 分析方法，设计算法并编程实现，同时指明准确计算 自相似性 (自放射性、广义 自相似性)、奇异性等通 奇异性指数需要注意的一些关键参数设置。 常被用来描述地球化学场特征、矿床空间分布、热液 成矿系统矿床 品位一吨位分形模型的非线性性质。 1 局部奇异性分析方法的基本思想 一 种五参数的多重级联模型最近被提出，该模型对 deWijs’s锌数据等多重分形经典数据的非对称奇异 奇异性理论指 出：奇异性是指在很小 的时间一 谱进行了成功建模，为现实世界中地球化学场多重 空间范围内具有 巨大能量释放或巨量物质形成的现 分形性的存在性提供 了理论依据。分形维数、多重 象(成秋 明，2006；Chengetal，2009)。假定在 E 分形对称度、多重分形奇异谱等参数从非线性角度 维 (拓扑)空间中有一 由质点凝聚成分维数为 D的 对描述成矿物质富集与贫化的统计规律发挥了重要 分形(如矿物的富集)，尺度上质点集合的质量 作用 (成秋明，2000，2007)，然而，对于矿产资源定 (L)与集合中质点的数 目Ⅳ()成正比，则该分形体 量预测中的 “定位”(矿产资源空间分布位置)预测 的密度可表示为(陈颐等，2005)： 功能则有所不足，需要将全局性统计转变为局部化 p： L (1) 模型。基于多重分形理论建立的局部奇异性分析模 收稿 日期：2014—08—