曼家寨锌矿体储量基于克里格法的计算
曼家寨露天采场西帮边坡原位直剪试验研究
关 系数 ｒ 。试验表 明 ， 西 帮边坡 原位 直 剪 试验 研 究 对 露 天采 场 边 坡稳 定性 评 价 及 工程 治理 优化 具
有重要 价值 。 关键 词 边坡 原位 直剪 试验 抗 剪 强度 变形
Ｉ ｎ － ｓ ｉ ｔ ｕ Ｄ ｉ ｒ ｅ ｃ ｔ Ｓ ｈ ｅ ａ ｒ Ｔ ｅ ｓ ｔ Ｓ ｔ ｕ ｄ ｙ ｏ ｆ Ｗｅ ｓ ｔ Ｓ ｌ ｏ ｐ ｅ ｏ ｆ Ｍａ ｎ ｊ ｉ ａ ｚ ｈ ａ ｉ Ｏ ｐ ｅ ｎ － ｐ ｉ ｔ
现
代
矿
业
总 第５ ３ １ 期
Ｍ ０ＲＤＥＮ ＭＩ ＮＩ ＮＧ
２ ０ １ ３ 年 ７月 第 ７期
・
采选工程 ・
曼家 寨 露天 采 场 西帮 边 坡原 位 直 剪试 验研 究
李 诚 项 宏海 段蔚 平
（ １ ． 云南华联锌铟股份有 限公 司； ２ ． 中钢 集团马鞍 山矿山研 究院有限公 司； ３ ． 金属矿 山安全与健康 国家重点 实验室）
ａ ｎ ｄ ｃ ｏ ｒ ｒ ｅ ｌ ａ ｔ ｉ ｏ ｎ ｃ ｏ ｅ ｉｃ ｆ ｉ ｅ ｎ ｔ ｒ ｗｅ ｒ ｅ ｏ ｂ ｔ ａ ｉ ｎ ｅ ｄ ｂ ｙ ｕ ｓ ｉ ｎｇ ｌ ｅ ａ ｓ ｔ ｓ ｑ ｕ ａ ｒ ｅ ｓ ｒ ｅ ｇ ｒ ｅ ｓ ｓ ｉ ｏ ｎ ．Ｔｈ ｅ ｔ ｅ ｓ ｔ ｉ ｎ ｄ ｉ ｃ ａ ｔ ｅ ｓ ｔ ｈ ａ ｔ ｉ ｎ— ｓ ｉ ｔ ｕ
不同储量估算方法在离子型稀土资源储量估算中的应用探索
- 143 -技 术 经 济 与 管 理0 引言离子吸附型稀土矿是1970年江西地质局第七地质大队（前身江西908大队）在龙南足洞矿区首次发现，在我国华南地区不断取得找矿突破，目前已成为我国优势战略矿产资源，是世界中重稀土元素主要来源，具有短期内其他国家无可替代的地位[1-3]。
离子吸附型稀土矿体多呈似层状沿全风化层分布，平面形态受风化壳形态的控制，呈阔叶状随地形而变化，边界一般受沟谷展布的控制，矿体倾角由山顶至山坡不断变陡。
矿体的厚度、品位受地貌位置影响较大，从山顶、山腰到山脚厚度、品位一般逐渐变小。
任何一个矿块的资源量都有唯一的实际值，在计算范围和参数相同的条件下，使用不同的计算方法，得出的结果均可接近实际值。
目前，国内地勘单位一般采用平面投影地质块段法，该方法适用范围广泛，但是受取样工程密度、地形起伏和人为选择矿块边界工程点影响，多数学者认为该方法存在较大的误差。
该文对多个不同情况的矿块使用多种方法估算稀土资源量进行对比，分析各种估算方法之间的差异，以便选择更加适合的资源量估算方法，为今后进行资源储量估算工作提供新的选择。
1 估算方法原理对比平面投影地质块段法是将矿块的单工程平均厚度、平均品位计算出来后，按矿块块段分别计算单块段平均厚度、平均品位，再根据投影的平面面积与平均厚度计算块段体积，体积与矿石密度计算矿石量，最后以矿石量和平均品位得到估算资源量。
这种方法是目前普遍采用的估算方式，然而，离子型稀土矿的矿体赋存形态与其他金属矿床存在较大的差异。
离子型稀土矿是稀土母岩在长期风化和化学腐蚀等作用下不断解离，在随水流向下迁移的过程中被黏土矿物吸附而形成的一种淋积型矿床，其矿体形态常受风化壳形态的控制[4-5]，当使用平面投影地质块段法时存在如下问题：1）受地形影响较大，基于其计算原理，地形变化越大，平面投影后的误差值越大。
2）矿体分块段后，相邻块段会不可避免地反复使用相同钻孔样品进行计算，受这种人工选择块段分界线的影响，误差值可达10%以上。
矿量计算方法
矿量计算方法LG GROUP system office room 【LGA16H-LGYY-LGUA8Q8-LGA162】资源量与储量计算方法储量（包括资源量，下同）计算方法的种类很多，有几何法（包括算术平均法、地质块段法、开采块段法、断面法、等高线法、线储量法、三角形法、最近地区法/多角形法），统计分析法（包括距离加权法、克里格法），以及SD法等等。
（一）地质块段法计算步骤：首先，在矿体投影图上，把矿体划分为需要计算储量的各种地质块段，如根据勘探控制程度划分的储量类别块段，根据地质特点和开采条件划分的矿石自然（工业）类型或工业品级块段或被构造线、河流、交通线等分割成的块段等；然后，主要用算术平均法求得各块段储量计算基本参数，进而计算各块段的体积和储量；所有的块段储量累加求和即整个矿体（或矿床）的总储量。
地质块段法储量计算参数表格式如表下所列。
表地质块段法储量计算表块段编号资源储量级别块段面积（m2）平均厚度(m）块段体积（m3）矿石体重（t/m3）矿石储量（资源量）平均品位（%）金属储量(t）备注需要指出，块段面积是在投影图上测定。
一般来讲，当用块段矿体平均真厚度计算体积时，块段矿体的真实面积S需用其投影面积S′及矿体平均倾斜面与投影面间的夹角α进行校正。
在下述情况下，可采用投影面积参加块段矿体的体积计算：①急倾斜矿体，储量计算在矿体垂直纵投影图上进行，可用投影面积与块段矿体平均水平（假）厚度的乘积求得块段矿体体积。
图在矿体垂直投影图上划分开采块段(a)、(b)—垂直平面纵投影图； (c)、(d)—立体图1—矿体块段投影； 2—矿体断面及取样位置②水平或缓倾斜矿体，在水平投影图上测定块段矿体的投影面积后，可用其与块段矿体的平均铅垂（假）厚度的乘积求得块段矿体体积。
优点：适用性强。
地质块段法适用于任何产状、形态的矿体，它具有不需另作复杂图件、计算方法简单的优点，并能根据需要划分块段，所以广泛使用。
矿产资源储量估算方法
几种常见的矿产资源储量估算方法固体储量估算方法主要是几何法和统计分析法。
一、几何法（一）断面法（剖面法）原理就是当矿体被一系列勘查断面横切为若干块段，就可以以这些断面图为基础，估算相邻两断面间的矿块储量乃至整个矿床储量。
分为垂直断面法和水平断面法。
第一步：计算体积1、当相邻两断面的矿体形状相似，且其相对面积差（S1－S2）÷S1小于40%时，用梯形体积公式V=(S1+S2)×L÷2。
其中V为两断面间的矿体体积；L为相邻两剖面间的距离；S1、S2为相邻两端面上的矿体面积。
2、当相邻两断面的矿体形状相似，且其相对面积差（S1－S2）/S1大于40%时，选用截锥体积公式，即V=（S1+S2+√S1×S2）×L÷3。
其中V为两断面间的矿体体积；L为相邻两剖面间的距离；S1、S2为相邻两端面上的矿体面积。
3、当相邻两断面的矿体形状不同，不论面积相差多少，除油一对应边相等时，可用梯形体积公式外，其余均应选用似角柱体（辛浦生）公式，即V=[（S1+S2）÷2+2S m]×L÷3 =（S1+S2+4S m）×L÷6。
其中V为两断面间的矿体体积；L为相邻两剖面间的距离；S1、S2为相邻两端面上的矿体面积。
S m为似角柱体的平均断面面积。
4、当在相邻的两剖面中只有一个剖面有面积，而另一剖面上矿体已尖灭，或矿体两段边缘部分的块段只有一个断面控制时，其体积计算可根据剖面上的矿体面积形状或矿体尖灭特点不同选择不同公式。
（1）当矿体作楔尖灭时，块段体积用楔形公式计算。
V=L×S÷2（2）当矿体作锥形尖灭时，块段体积可用锥形公式计算。
V=L×S÷3第二步，计算两剖面间块段的矿石储量Q=V×d。
其中Q为块段矿石储量，V为块段的矿体体积，d为块段矿石平均体重。
第三步，计算出两剖面间块段的金属储量P=Q×C。
克里金法 案例
克里金法案例【最新版】目录1.克里金法的定义和原理2.克里金法的应用案例3.克里金法的优缺点分析正文【克里金法案例】克里金法是一种插值方法，主要用于空间数据的预测和模拟。
该方法基于距离衰减原理，通过计算周围已知数据点的加权平均值，预测或模拟未知数据点的值。
克里金法具有较强的理论基础和实用性，被广泛应用于地理信息系统、环境科学、地质勘探等领域。
本文将通过具体案例，介绍克里金法的应用及其优缺点。
一、克里金法的定义和原理克里金法（Kriging）是一种插值方法，其名称来源于南非的克里金（Kruger）金矿。
该方法是由南非矿业工程师丹尼斯·格里高利（Dennis G.Krige）于 1951 年提出的，用于预测金矿中的黄金含量。
后来，克里金法逐渐被应用于其他领域，成为一种重要的空间数据插值方法。
克里金法的原理是基于距离衰减的，即一个数据点对预测值的贡献与其距离成反比。
具体来说，克里金法通过计算已知数据点到预测点的距离，然后根据距离的大小赋予不同的权重，最后计算权重的加权平均值，得到预测值。
二、克里金法的应用案例1.地理信息系统：在地理信息系统中，克里金法可以用于空间数据的预测和模拟，如地形高程、土地利用、土壤类型等地理信息的预测。
2.环境科学：在环境科学领域，克里金法可以用于预测污染物的分布，如水质中的重金属含量、大气中的污染物浓度等。
3.地质勘探：在地质勘探领域，克里金法可以用于预测矿产资源的分布，如金矿、铜矿等。
三、克里金法的优缺点分析1.优点：克里金法具有较强的理论基础，可以较好地模拟空间数据的变化规律；同时，该方法具有较高的计算效率，适用于大规模空间数据的预测和模拟。
2.缺点：克里金法的预测结果受输入数据质量的影响较大，如果输入数据存在噪声或误差，预测结果也可能出现偏差；此外，克里金法对于数据点的分布和密度有一定的要求，当数据点分布稀疏或不规律时，预测结果可能不准确。
综上所述，克里金法作为一种插值方法，在空间数据的预测和模拟方面具有广泛的应用。
曼家寨采场东部边坡地下水同位素测试研究
具有 明 显承压 性 ， 地下 水 主要 为上层 滞水 或潜 水 ， 可
以利用 疏 水工 程疏 干边 坡 内地 下 水 ， 以此 增 强 边 坡
稳 定 。见 图 ２ 。
１ 垂 向流 测 试
１ １ 测试 原理 ．
７ ９ Ｏ
２
Ｏ
式（ ） ８ 即为 单 孔稀 释 指 数 定 律 ， 即渗 透 速 度 的
Ｉ ｌ
● ●
．
ｆ ４１
计 算公 式 。
９ ５
ｌ
■
５ ２
又 因为
Ｖ ：ａ ｆ ｗ Ｖ， （） ５
２ ２ 地 下水 流速 测试结 果 ．
ｌ
（ ．云 南华联锌铟 股份有限公 司；．中钢集 团马鞍 山矿 山研 究院有限公 司； １ ２
３ ．金属矿 山安 全与健 康 国家重点 实验 室 ；． ４ 河海 大学土木 ＿程 学院） Ｔ －
摘
要
通过 现 场试验 ， 研 究 了边坡 垂 向流 、 平流速 、 下水 流 向的基 础 上 ， 出 了地 下水 在 水 地 得
因此
２ａ Ｖ ｆ Ｂ ：
— —
深度上 的地 下水 流速 ， 试验 点流 速均 大于 ０ ００１ 各 ．０
Ｏ
●
Ｏ
●
ｍ ｄ 其 中在 １ ｍ 以上深 度流速 大 于 １０ ｍ ｄ 见 表 ／， ５ ． ／ （
３ ９
９
ｌ
１， ）通过 结果 分析 表 明边坡 地 层 在 １ ５ｍ深 度 为 ～１
孔 １ ７个 ， 其 中 的 ５个 钻孔 进行 了地 下水 流 速和 流 对 向的测试 。流 速 、 向测试 在天 然 条件下 进 行 ， 孔 流 钻 编号 分 别 为 Ｂ Ｋ ３４、 Ｓ ４ — 、 Ｓ ４ — Ｂ Ｋ ５２ Ｓ ４ － Ｂ Ｋ ３３ Ｂ Ｋ ３２、 Ｓ ３ — 和 以前 的一 个 老 钻 孔 ， 老 钻 孔 靠 近 １ 该 ９线 ， 在
kriging基础知识ppt课件
（应用随机函数理论）
最新版整理ppt
井眼 地震
4
第一节 基本原理
一、随机变量与随机函数 1. 随机变量
为一个实值变量，可根据概率分布取不同的值。 每次取值（观测）结果z为一个确定的数值，称为 随机变量Z的一个实现。
P
最新版整理ppt
5
连续变量：
累积分布函数（cdf）
Z (u)
cumulative distribution function
F (u ;z ) P o{ rZ b (u ) z }
P
条件累积分布函数（ccdf）后验 conditional cumulative distribution function
F ( u ;z |( n ) ) P o { Z r( b u ) z |( n )}
第二讲
克里金插值
克里金方法（Kriging）, 是以南非矿业 工程师D.G.Krige (克里格)名字命名的一项 实用空间估计技术，是地质统计学 的重要
组成部分，也是地质统计学的核心。
最新版整理ppt
1
地质统计学
由法国巴黎国立高等矿业学院G．马特隆教授于 1962年所创立。 主要是为解决矿床储量计算和误差估计问题而 发展起来的
（将空间位置作为随机函数的自变量）
•空间一点处的观测值可解释为一个随机变量在该点
处的一个随机实现。
• 空间各点处随机变量的集合构成一个随机函数。
（可以应用随机函数理论解决插值和模拟问题）
最新版整理ppt
14
考虑邻近点，推断待估点 ----空间统计推断要求平稳假设
严格平稳
F ( u 1 , , u K ; z 1 , , z K ) F ( u 1 h , , u K h ; z 1 , , z K )
克里金法 案例
克里金法案例【原创版】目录1.克里金法的定义与原理2.克里金法的应用案例3.克里金法的优缺点正文【克里金法】克里金法，全称克里金插值法（Kriging Interpolation），是一种基于随机场理论的插值方法，主要应用于空间数据的预测和模拟。
它是由南非矿业工程师丹尼尔·克里金（Daniel Krige）于 1951 年提出的，用于解决矿产资源勘探中的空间数据预测问题。
克里金法的基本原理是：假设空间数据由一个或多个随机场构成，通过构建随机场模型并求解其协方差矩阵，从而实现对未知数据的预测和模拟。
【应用案例】克里金法在许多领域都有广泛应用，如地质勘探、环境监测、气象预报等。
这里举一个地质勘探的案例：假设我们在某地区进行矿产资源勘探，已经获得了一系列钻孔的矿产品位数据。
我们需要预测该区域内其他位置的矿产品位。
这时，我们可以使用克里金法来解决这个问题。
首先，根据已有的钻孔数据，构建矿产品位的随机场模型。
然后，通过求解协方差矩阵，可以得到任意位置的矿产品位预测值。
这样，我们就可以预测该区域内其他位置的矿产品位，为矿产资源的开发和利用提供科学依据。
【优缺点】克里金法具有以下优点：1.可以处理空间相关数据，考虑数据的空间变异特性；2.具有较强的理论依据，预测结果具有较高的可信度；3.可以处理不完全数据，适用于资料匮乏的情况。
然而，克里金法也存在一些缺点：1.计算过程较为复杂，需要求解协方差矩阵；2.对输入数据的质量要求较高，数据质量会影响预测结果；3.克里金法假设数据遵循特定的统计模型，当实际数据不符合假设时，预测结果可能会出现偏差。
综上所述，克里金法是一种强大的空间数据预测方法，在许多领域具有广泛的应用前景。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
５ 变 异 函数 的交 叉 验 证
克 里 格估 值 是 通 过 矿 体 内 各 组 合 样 数 据 建 立 理 论 变 异
¨ 。 ． －． 一 … ～ 。一 。一 －． －． 。 一 。一 一 ． － ． ＋ ’
Ｏｒ ｅ ｒ ｈｉ ｇ ａ ｒ ｂ ｄｙ Ｓ ｅ ｏ ｅ ｏ ｈ ｅ ｓ ａ ｃ ｎ ｔｏ ｅ ｏ ’ ｎｄ ｚ ｎ ｎ ｔ ｅ ｂａ ｉ ｆ ａ ｌ ｚ ｎ ｅ ｌ ｇ ｃ ｌｐｒ ｐｅ ｔｅ ｓｓ ｏ ｎａ ｙ ｉ ｇ ｇ ｏ ｏ ｉ ａ ｏ ｒ ｉ ｓ
频率
瓣
％）
一
ｌ 一 — — — Ｏ Ｏ ８ ６ ４ ２
２
４
６
图 １ 曼家寨矿段 ｚ ｎ品 位 组 合 样 对 数 统 计 值
Ｓ ７ １ ●
图
格 残差 直 方 图
攮 率瓣
囊 ｝ ：嚣 量 墨 要耋
：
能§ 嚣 荟 Ｍ 耋 。
● ２ Ｓ ７ １■
图６ Ｚ ｎ克 里 格 残 差／ 准差 直方 图 标
Ｕ ● ■ ２ ● Ｔ ｓ ， ● １ Ｓ ● ２ Ｈ ４ ７ ｊ ａ １ ● １ ０ Ｈ ２ ２
图 ２ Ｚ 走 向方 向 实验 变 异 函 数 图 ｌ｛｜． ｎ 蓥垂萎蓥ｌ ； 室耋兰嚣蚕 Ｉｌｌ１
； ：
．
： ．６ 豫ｌｌ｜一 的建 立 ｚ ｎ ．Ｉ｜ｌ Ｉ；ｉｉ ． 品位 模 型 ｌｌｊ
．
：
兰
ｌ ｉ
． 克
； 里
品 位 建 模 就 是 对 建 模 范 围 内 各 个 单 元 块 的 品 位 根 据 其 ｌｌ ｌ一
兰篓～
周 围 的 地 质 信 息 进 行 估 值 ， 前 述 的 一 切 工 作 ， 其 目 的 都 是
为 了准 确 地 建 立 品位 模 型 。本 次 采 用 Ｓ ｒ ａ 提 供 的 克 里 格 ｕｐ ｃ
矿 区 为单 斜 构 造 ，地 层 走 向南 北 ， 向 西 倾 斜 ，倾 角 一 般 １。 ５。 区 内 出露 地 层 为 中 寒 武 统 田蓬 组 区 域 变 质 岩 ， Ｏ～３ 。 变 质 程 度 中等 ，锡 、锌 、铜 、铅 矿 化 背 景 较 高 ， 具 矿 源 层
特征 。
４ 变 异 函数 的 计 算 和 拟 合
Ｗ ＡＮＧ ａ － ｎ ，ＺＨ０Ｕ Ｂ ｎ ， Ｌ ｉ ，ＱＩＰ ｎ ｗｅ ｃｈ ｏ ｗｅ ｉ ＩＬ ａ－ ｎ
（ ． ａ ｓ ｉ ｉ ｉ ｇ Ｃｏ ａ ｙ ｏ ｎ ｈ ａ Ｉ ｏ ＆ Ｓ ｅ ｌＧｒ ｕ ｒ ｏ ａ ｉ ｎ Ｂ ｉｈ ｎ １ ４ ０ １ Ｂ ｎ ｈ ｎ ｎ ｍｐ ｎ ｆＴｏ ｇ ｕ ｒ ｎ Ｍ ｔ ｅ ｏ ｐ Ｃｏ ｐ ｒ ｔ ， ａ ｓ ａ ３ ３ ４： ｏ
法 进 行 估 值 。 以建 立 的 马 关 都 龙 曼 家 寨 锡 锌 多 金 属 矿 床 三 维 实 体 模 型 为 约 束 条 件 ，对 矿 体 进 行 品 位 克 里 格 估 值 。 综 合 考 虑 矿 山 现 有 采 矿 方 法 、矿 区 勘 探 网 度 及 元 素 变 异 函 数
倾 伏 角 ７，倾 向方 向方 位 角 为 ２ ０，倾 伏 角 ２ 。 。 ７。 ８ ，垂 向方 向 方位 角为 ９。 Ｏ ，倾 伏 角 ６ 。 ２，得 到 的 实 验 变 异 函 数 曲 线 较 好 （ ２ ４ 。 在此 基 础 上 ，计 算 了各 方 向 的 理 论 变 异 函数 参 图 ～ ） 数 （ １。 表 ）
０” ０ ，北 纬 ２ 。３ ０ ２ 。 ５０ ｒ 。 ２５ ０ ～ ２５ ０， ｒ
立 变 量 的理 论 模 型 ；按 开 采 要 求 划 分 体 积 及 形 状 相 同 的 待 估 块 段 ，利 用 待 估 块 段 的 内 外 信 息 ，计 算 待 估 块 段 平 均 参
数和矿石及金属储量 ｌ。 ２ ］
２ Ｃ ｉａＧｅ ｌｇ ｉｅ ｓｔ ． ｈ ｎ ｏｏ ｙＵｎｖ ｒ ｉ ｙ，Ｂｅｊ ｇ １ ０ ８ ，Ｃｈｎ ） ｉｎ ０ ０ ３ ｉ ｉａ
Ａｂ ｔ ａ ｔ Ｂ ｓ ｄ ｏ ｎ ｌｓ ｓ ｏ ． ｒ ｂ ｄ ｅ ｌ ｇｃ ｌｐ ｏ ｅ ｔｅ ｎ ｅ ｔ ｎ ｃ ｆ ｒ ， ｔ ｅ ｔ ｅ ｒ ｆ“ ａ ｌ ｄ ｓｏ ａ ｉｎ ｓｒ ｃ ： ａ ｅ ｎ ａ ａ ｙ ｉ ｆ Ｎｏ ８ ｏ ｅ ｏ ｙＳ ｇ ｏ ｏ ｉａ ｒ ｐ ｒｉｓ ａ ｄ ｔ ｃｏ ｉ ｏ ｍ ｈ ｈ ｏ ｙ ｏ ｆ ｕ ｔ ｉｌ ｃ ｔ ｏ
２ 克 里 格 法 概 述
克里 格 法 是地 质 统 计 学 的 一 个 主 要 方 法 ， 是 一 种 最 优 无 偏 估 量 的储 量计 算 方 法 。该 方 法 以变 异 函 数 为 工 具 ，建
收 稿 日期 ：２ １— ０ —２ ０１ ３ １ 作 者 简 介 ：伍 伟 （ ９ １ ） １ ８ 一 ，男 （ 族 ） 汉 ，云 南 楚 雄 人 ，昆 明冶 金 研 究 院博 士 ，工 程师 ，主 要研 究 方 向 ：矿产 经 济 ，工 艺 矿 物 学研 究 。
ｌ ｃ ｔ ｎ ｏ ｌ ｄ ｏ ｅ ｂ ｄ ． Ｔｈ ｓ ｏ ｅ ｂ ｄ ｓ ｉ ｇ ｉ ａ ｏｄ ｄ， ｗｈｃ ｅｐ ｖ ｒ ｏ ｈ ｏ ｔａ ｉｔｏ ｆ“ ｎ ｒ ａ ｉ ｇ ｏ ｅ ｏ ａ ｉ ｆｂ ｉ ｒ ｏ ｙ ｏ ｎ ｕ ， ｒ ｏ ｙＳ ｍｉｓｎ ｓ ｖ ｉ ｅ ｉ ｈ ｈ ｌ ｓ ｏ ｅ ｃ ｍｅ ｔ ｅ ｃ ｎ ｒ ｄｃ ｉ ｎ ｏ ｉ ｃ ｅ ｓｎ ｒ ｐ ｏ ｕ ｔ ｎ ｗｉｈ ｄ ｃ ｅ ｓ ｄ ｒ ｓ ｕ ｃ ”ａ ｄ ｂ ｉｇ ｃ ｎ ｍｉ ｅ ｅ ｉ． ｒ ｄ ｃ ｉ ｔ ｅ ｒ ａ ｅ ｅ ｏ ｒ ｅ ｎ ｒｎ ｓ ｅ ｏ ｏ ｃ ｂ ｎ ｆｔ ｏ
０ ２ ４ ６ ８ １ ０ １ ０ １ ０ １ ０ １ ０ ０ ｍ ０ ０ ０ ０ Ｏ ２ ４ ６ ８ ２ ０／
的特 征 等 因 素 ，确 定 品 位 模 型 单 元 块 尺 寸 为 行 宽 ×歹 宽 × ０ 层 厚 ５ｍ×５ｒ×２ ５ｒ， 采 用 球 状 变 异 函 数 模 型 ，分 别 沿 ｎ ． ｎ 矿 体 走 向 、倾 向 、厚 度 ３个 方 向 进 行 变 异 函 数 分 析 ，最 终 得 到 了 马关 都 龙 曼 家 寨 锡 锌 多 金 属 矿 床 的 品 位模 型 。
摘 要 ：阐 述 了 应 用 克 里 格 法 建 立 了曼 家 寨 锌 矿 体 的 品 位 模 型 ，并 进 行 了 储 量 计 算 。用 该 计 算 结 果 与 传 统 方 法 的 计算 结 果 进 行 对 比 ，证 明 了该 方 法 的准 确 性 。 关 键 词 ：曼 家 寨 锌 矿 体 ；克 里 格 法 ；品 位 模 型 ；储 量 计 算 中 图 分 类 号 ：Ｐ６ ４ ７ ２ ． 文 献 标 识 码 ：Ａ 文 章 编 号 ：１ ７ — ８ ５ （０ １ ５ ０ ３ ３ ６ １ ５ ０ ２ １ ）０ —０ １ —０
１ 矿 区地 质 概 况
都 龙 曼 家 寨 锡 锌 多 金 属 矿 区 ，位 于 云 南 省 东 南 部 ，隶 属 于 文 山壮 族 、 苗 族 自治 州 马 关 县 都 龙 镇 ，矿 区 面 积 约 ５ｋ ，呈 南 北 向展 布 ，位 于 老 君 山 花 岗 岩 体西 南 侧 外 接 触 ｍ 带 ，北 面 和 东 面 以 Ｆ 断 裂 为 界 ，宽 约 １ ５ｋ 。 ． ｍ， 向南 延 至 中 越 国境 线 ，长 约 ８ｋ ｍ，构 成 南 北 向 展 布 的 锡 、 锌 、 铜 、银 等 多 金 属 矿 带 。其 地 理 坐 标 为 ：东 经 １４３ ０ ～ １ ４３ ０ 。２０ ０ 。 ３
１ ４
矿 业 工 程
第９ 第 ５ 卷 期
８ 函数 模 型 的基 础 上 ，应 用 ７ 建 ５ 理 ３变 ２ 函 模 型 来 对 已 ６立 的 ４ 论 异 ● 数
间 的误 差 小 于 ｌ ；两 标 准差 的 残 差 参 数 接 近 ９ ； ｎ的 ５ ５ ｚ
异 函数 进 行 交 叉 验 证 ，得 到 验 证 结 果 见 表 ２ 图 ５ 。 及 ～６
从 正 态 分 布 ］ 。说 明 所 选 择 的 变 异 函 数 理 论 参 数 很 接 近 实 际 ，能 用 于 Ｚ ｎ品 位 的 克 里 格 估 值 。 表 ２ 方 差 检 验 统 计 结 果
ｐａ ｉｙ ｉ ｓ ｄ ｔ ｉ ｅ ｓ ａ ｃ ｎ ａ ｅ ｄｙＳ ｎ ｏ ｄ Ｏ ｎｖ ｓ ｉ ｔ ｈｅ ｓｚ ｒｔ ” ｓ ｕ ｅ ｏ ｇｕｄｅ ｏｒ ｅ ｒ ｈｉｇ ｔｏｒ ｂｏ ｅ ｄ ｚ ｎｅ ａｎ ｔ ｉ ｅ ｔｇａ ｅ ｔ ｉｅ， ｏｃ ｕ ｒ ｎｃ ｃ ｒ ｅ ｅ， ３一 Ｄ ｏｒ ｎ ｓ ａ ｉ ｆ ｍ ａ ｄ ｐ ｔａｌ
３ 样 品数 据 的 统计 分 析
地 质 统 计 学 的计 算 要 求 有 效 数 据 必 须 确 定 在 同 等 大 小 的承 载 上 ，以保 持 其 均 匀 性 ，也 就 是 要 保 证 每 个 样 品 的 权 大小一致 ，以达到 等效 的支撑¨ ］ ２ ，为 此 ，必 须 进 行 样 品 的组 合 。根 据 该 矿 山 的 采 矿 技 术 ，可 将 组 合 样 的 样 长 确 定 为 １ｍ，对 Ｚ ｎ进 行 统 计 分 析 ，从 图 １ 看 出 ，Ｚ 可 ｎ组 合 样 基 本 服 从 对 数 正 态分 布 ，故 可 采 用 对 数 克 里 格 法 对 矿 体 进 行 品位 估 值 。
第 ９ 卷 第 ５ 期 ２１ ０ １年 １ Ｏ月
矿 业 工 程
Ｍ ｉ ｉ ｇ ｎｅ ｒ ｎ ｎ ｎｇ Ｅｎ ｉ ｅ ｉ ｇ １ ３