总第22期成矿预测中非参数蒙特卡洛模拟方法研究.pdfVIP

1g91年第 3期 矿 产 与 地 质 第 5 卷 1991年 6 月 M INERAL RESOURCES AND GEOLOGY 总第 22期 成矿预测中非参数蒙特卡洛模拟方法研究 赵玉琛 安徽省地质矿产局 322队地质研究所) 摘要 非参数方法是在成矿预谢中．当统计数据不服从或不完全服从规则分布时的一种蒙 特卡洛模拟形式 非毒数方法采用拉格朗 日插值函数和三次条样函数．并均满足密度函 数的条悻．可 以对统计数据进行拙台．借 以建立各种抽样模型 即便是离散 型数据。也 可用牛顿插值法为其建立非参数抽样模型。因而。它具有算法统一 。编程简单和有利于 计算机 自动化模拟等优点，具有广 阔的应用前景 关键词 ：成矿预测，非参数方法，抽样模型 蒙特卡洛模拟是在概率意义下 的一种合理的抽样统计方法。在地学领域 内． 目前主要用 于矿产资源预测 ，其抽样模型通常是靠选配经典分布的方法来建立 ，即通过统计样本直方 图 的峰度 、偏度检验 ．用形状近似的已知分布去拟合 ，再用特征参数 (如平均值 、均方差等) 经过理论计算构塑出分布函数：这种依赖统计参数的拟合形式称参数方法 。 实际经验证明，一些统计数据并不根 从l或不完全眼从经典概率的规则分布，而有些数据 ． 如成矿有利度等控制性参数，究竟服从l何种 分布，还不能从l理论上给予说明；这种情况下． 参数方法就很棘手。为此，建议用拉格朗El插值函数或样条函数等非参数方法，通过数值逼 近的办法建立其抽样模型。它具有无需论证数据理论分布型式，适于各种情况，算法统一和 编程简单等优点。 1 非参数方法的理论依据 矿产资源预测工作中的抽样模型的建立 ，往往依赖少数统计数据，甚至无法画出直方图． 而 只是一批数 组 (如文献 ：1的储量模型)=非参数方法在于寻找一个合适 的函数 f(x)-它 能拟台这批数组或直方图，且满足密度函数的条件，构塑出抽样摸型。这在数学上归结为插 值 问题。 用插值函数拟台数组或直方图有多种 算法，比较流行和便于编程的是本文建议的拉格朗 日函数和样条函数。下面首先讨论它们是否满足密度函数的两个先决条件，即 f(x) 0和 lf： )dx=l： 为讨论方便．以下均考虑 为等距节点的情况．这并不破坏论述的一般性 1．1 拉格朗日插值函数 (Pn(x)1 Pn(x)是n次多项式．一般形式为PnIx)=Co+C +…… +CnXn或写成求和形式：Pn(x)= ykL (x) 式 中 ii为模拟的直方图柱的个数， 为各拄的高度 ．Lk(x)为与柱宽有关的函数。由于y 0．Lk(x)≥O，故Pn(x) 0，满足f(x) 0的要求。 直方图具有I： f(x)dx=L的性质，若在其一)。到)。区间内插几个节点，则要求 下式 成立 ： l9l J’：f(dx)dxJ：Pn(x)dx=高AkY：l 式中A=旦n垦_Jf0： {_二dt，其中a、b分别相当直方图的两个端点位置，按题意令 I；0 k — I 。 i± k Q }1 a一。(．b一+。c·由于已经证咀壹0Joi。_}dt=l[23，因而，Ak相当于各