矿体几何制图数据的稳健处理-中国科技论文在线.PDFVIP

第2o卷第2期 测 绘 学 报 VoI．20，No．2 1991年5月 AGTA GEODAETIGA etGARTOGRAPHIGA SINICA a1 — L 矿体几何制图数据的稳健处理 汪 云 甲 (中国矿业大学) 提 要 本文讨论了矿体几何制图数据的特点及常用处理方法存在的问题，比较分析了游动 中值与游动平均的修匀效果及优缺点，研究了基于 Huber-M 估计的稳健游动平均修匀 与稳健趋势面分析方法。理论分析、实例计算、蒙特卡洛模拟结果表明：这两种方法对 异值影响有较好的抑制能力，它们既保留了传统数据处理方法的优点，又克服了某些不 足，在剩余不服从正态分布时亦能取得较好的效果。 一 、 引 言 矿体几何学是矿山测量专业学科的重要组成部分，是测量、地质、采矿之间的边缘学 科，它主要研究矿体形状和矿产特性的空间分布规律，利用图解模型和数学模型解决与矿体 形状及矿产特性有关的地质采矿问题。矿体几何制图是矿体几何学的中心内容，包括绘制各 种矿体几何图件、建立矿床模型等等[1]。 矿体几何制图数据系由对矿体取样而得。由于试样的代表性与不可重复性，同时由于矿 体指标本身的特性，使得矿体指标观测值既含有规律性变化，又含有局部的随机误差。因此 在进行矿体几何制图前应对原始数据进行处理，以抑制随机部分，将规律部分有效地显示出 来，这对某些矿体性质的几何制图尤为重要。本文首先讨论矿体几何制图数据的特点及常用 数据处理方法存在的问题，在此基础上研究矿体几何制图数据的稳健处理方法。 二、矿体几何制图数据的特点及常用处理方法中的问题 常用的矿体几何制图数据处理方法有两种：图解法与解析法，通常前者称为修匀，后者 称为趋势面分析，而最一般的修匀方法即为游动平均。显然，这两种方法均要求剩余服从正 态分布，只有在这样的条件下它们才可能取得较好的效果。但是，矿体几何制图数据 (即矿 体指标观测值)有以下特点： (1)总体分布的多样性与不规贝U性。矿体某一特定指标大小及其空间展布规律是漫长地 质年代多次地质作用而形成的结果，由于诸种地质条件及因素的影响，在不同矿体甚至同一 矿体的不同区域，矿体指标的总体常呈不同的分布，而且这些分布往往难以用理论的函数来 拟合。同时，样本中往往含有其它总体的样本，甚至就是两个或更多总体的混合物，因而其属 性的分布不是某种单一分布，而是多重分布的复合。如研究物采化采数据时的基本假设就是 数据中至少包含两个总体——背景总体与矿化总体，物采与化采的目的正是从背景总体中鉴 ‘本文1990年8月l8日收到。 转载 中国科技论文在线 92 测 绘 学 报 2O卷 别出属矿化总体的样品(1111异点)。 (2)矿体指标观测值既包含规律成分，又包含随机误差，异值或 风“暴样品 往往较多。 如化采数据，异值个数常常大于5％，有时甚至大于20 。 (3)矿体指标观测值既存在技术误差，又存在代表性误差，一般不能或不作重复测量。 上述这些特点，决定了矿体几何制图数据不完全服从正态分布，对某些金属品位及地球 化学数据，一般不服从正态分布，对煤矿床，以前人们认为煤厚服从正态分布，但作者最近 对平顶山矿区七个矿的煤厚进行了统计检验 ，发