矿床体视化的一种空间规则体数据生成方法

D0I:10.13374/i.issm1001053x.2002.06.003 第24卷第6期 北京科技大学学报 Vol.24 No.6 2002年12月 Journal of University of Science and Technology Beijing Dec.2002 矿床体视化的一种空间规则体数据生成方法 李翠平李仲学胡乃联 北京科技大学矿业研究所，北京100083 摘要系统地论述了矿床体视化过程中由原始采样数据生成空间规则体数据的实现方法， 分析了能够形成地矿体空间体数据的主要插值方法的优缺点，在充分考虑矿石品位的空间相 关性和变异性的前提下，利用Shepard法实现了空间体数据的生成. 关键词体视化；体数据；Shepard法 分类号TP319.4:TD853 矿床及其开采设计的可视化技术研究正在 钻孔的相应位置.在处理地勘钻孔数据时，若用 成为地矿工程数字化的一个重要发展方向，. 曲线形态来表达钻孔，即用微分的形式建立钻 在矿床体视化仿真过程中，为了实现矿体的真 孔的参数方程，自然是比较精确的.但是，这种 三维显示及矿床属性（品位、储量）的计算与处 做法计算量大，而且一般情况下钻杆具有很强 理，必须对有关地质数据库中的原始采样数据 的刚度，钻孔的弯曲程度不会很大，所以用微分 进行预处理，使之转化为满足需要的三维体数 方法所提高的精确程度相对要确定矿岩属性的 据67刀 空间分布和空间变化来说其意义不大，因此，本 根据对某铁矿部分数据的分析，本文给出 研究中将钻孔处理为折线形式. 了一种由原始地质采样数据生成空间规则体数 对于生产勘探而言，因为勘探目的是探明 据的方法，为建立矿床体视化仿真系统奠定了 相应生产水平的矿岩分布，故一般钻探深度都 基础. 不大，最深也很少超过两个水平.所以，钻孔质 量一般较好，测斜次数由于钻孔深度的不同一 1样品数据到空间点数据的转化 般不多于2次. 将钻孔样品数据转化为空间点数据是矿床 在数据处理的过程中，当钻孔的测斜次数 n>1时，近似认为钻孔在空间呈现n段的折线 体视化中数据预处理的首要工作，由于采样中 形式.各折点分别处在相应的钻孔测斜点处，这 的样品数据是线数据（样品具有长度），而不是 点数据，为实现数据的规则化处理，使之得到规 样任意折点以及折线上的任意一点的位置坐标 则网格的三维体数据，需将样品数据转化为三 可通过相应测斜点的孔深、倾角和方位角确定 如图1所示，设钻孔在任意深度上的孔深 维空间点数据 在地质勘探工作中，钻探深度都比较大，一 般为数百米，有时达上千米.同时，由于钻孔所 穿过岩层的工程条件及钻杆刚性的影响，所形 (x.y,z) 成的钻孔在三维空间实际上是以曲线的形态存 在的，即任意位置上的钻孔深度实际上是孔口 Xo.Yo.Ze) 到该点的曲线长度，故地质工作人员在实际勘 测中经常进行不同深度的测斜工作，用以标识 收稿日期20010604李翠平女，28岁，博士研究生 *国家自然科学基金资助项目QN0.79870059)和教育部高校骨 图1钻孔在任意深度的空间坐标示意图 干教师资助计划项目 Fig.1 Coordinates of any point along a drillhole

第 2 4 卷 第 6 期 2 0 0 2 年 1 2 月 北 京 科 技 大 学 学 报 J o u r n a l o f U n i v e r s iyt o f s e ic n e e a n d 介 e h n o l o yg B e ij i n g V b l 一 2 4 N O 一 6 D e C . 2 0 0 2 矿床体视化 的一种 空 间规则体数据生成 方法 李翠平 李仲 学 胡 乃 联 北京科技大学矿业研究所 , 北京 10 0 0 83 摘 要 系统地 论述 了矿床体视化过 程 中由原 始采样 数据 生成空 间规则体数据 的实现方 法 , 分析 了能够形 成地矿 体空 间体数据 的主要 插值方 法的优缺点 , 在充 分考虑 矿石 品位 的空 间相 关性 和变异性 的前提 下 , 利用 Sh eP ar d 法 实现 了空间体数 据 的生 成 . 关键 词 体视化 ; 体数据 ; s h叩adr 法 分类 号 T P 3 1 9 . 4 ; T D 8 5 3 矿床及其开采设计的 可视化技术研究正在 成为地矿工程数字化 的一个重要发展方 向 `1一 .s] 在矿床体视化仿真过程 中 , 为了 实现矿体 的真 三 维显示及矿 床属性 (品位 、 储量 ) 的计算与处 理 , 必须对有 关地质数据库 中的原始采样 数据 进行预处理 , 使之转化为满足需要的三 维体数 据 “ , 7 , . 根据 对某铁矿部分数据 的分析 , 本文给 出 了一种由原始地质采样数据生成空 间规则体数 据 的方法 , 为建立矿床体视化仿真系统奠定 了 基础 . 1 样品数据到空 间点数据的转化 将钻孔样品数据转化为空 间点数据是矿床 体视化 中数据预处理 的首要工作 . 由于采样 中 的样 品数据是线数据 (样 品具有长度) , 而不 是 点数据 , 为实现数据的规则化处理 , 使之得到规 则 网格的三维体 数据 , 需将样品 数据转化为三 维空 间点数据 . 在地质勘探工作 中 , 钻探深度都 比较大 , 一 般为数百 米 , 有时达上千米 . 同时 , 由于 钻孔所 穿过岩层 的工程条件及钻杆刚性的影响 , 所 形 成 的钻孔在三 维空 间实际上是 以 曲线的形态存 在的 , 即 任意位置上的钻 孔深度 实际上是孔 口 到该 点的曲线长度 , 故地 质工作人员在实 际勘 测 中经常进行不 同深度 的测 斜工作 , 用 以标识 钻孔的相应位置 . 在处理地勘钻孔数据时 , 若用 曲线形态来表达钻孔 , 即 用 微分的 形式建立钻 孔 的参数方程 , 自然是 比较精确 的 . 但是 , 这种 做法计算 量大 , 而 且一般情况下钻杆具有很强 的 刚度 , 钻孔 的弯曲程度不会很大 , 所以用微分 方法所提高的精确程度相对要确定矿岩属 性的 空 间分布和 空 间变化来说其意义 不大 . 因此 , 本 研究 中将 钻孔 处理为折线形 式 . 对 于 生 产勘探 而言 , 因为勘探 目的是探 明 相应生产水平 的矿岩分布 , 故一般钻探深度都 不大 , 最深也很少 超过两个水平 . 所 以 , 钻孔质 量一般较好 , 测斜次数 由于钻孔深度 的不 同一 般不 多于 2 次 . 在数据处理的 过程 中 , 当钻孔 的测斜 次数 n > 1 时 , 近似认 为钻孔 在空 间呈 现 n 段 的折线 形式 . 各折点分别处在相应 的钻孔测 斜点处 , 这 样任意折点以及折线上的任意一点的位置坐标 可通 过相应测斜点的孔深 、 倾 角和方位角确定 . 如 图 1 所示 , 设 钻孔在任意深度上 的 孔深 收稿 日期 2 0 01 刃 6 es o 4 李翠 平 女 , 28 岁 , 博士研究生 * 国家 自然科 学基金资 助项 目(N .0 7 9 8 7 00 59) 和教育部高校骨 干教 师资助计划 项 目 图 1 钻 孔在任 意深 度的空 间坐标示 意图 F i g · 1 C o o r d in a te s o f a n y P o i n t a l o n g a d r i l l h o l e DOI: 10. 13374 /j . issn1001 -053x. 2002. 06. 003

·590 北京科技大学学报 2002年第6期 为1，方位角为a,倾角为B,此深度对应折线的起 层的序号；△x,△y,△z分别为网格点在三个方向 点坐标为(，，z),起点处的孔深为l,钻孔在此 上的间距，三者之间未必相等，当相等时为均匀 深度的相对孔深为，'=l-lo,1"=cos3.则此点 网格结构化数据，否则为规则网格结构化数据； 的坐标(xy,z)为： f八x,y,z)为体数据在网格点(，，)处的函数值. [x=x+/'cosa=xo+/cosBcosa 数据的空间插值是地矿工程计算机体视化 y=yo+l'sina=yo+lcosB sina. (1) 过程中必不可少的一项工作，插值算法的精度 z=z0+l sinB 和效率将直接影响地矿工程的三维显示效果 这样，对所有钻孔均从孔口水平（就实例而 通过对现有插值方法的分析，可以用来形成地 言，为-288m)开始，按垂直深度1m由样品数 矿体空间体数据的插值方法主要有Shepard:法、 据求得相应的空间点坐标.而对应的各空间点 四面体剖分法和克立格法 的属性（品位和岩性）值的确定方法见图2所示。 Shepard法的优点是计算方便，易于理解； 缺点是没有考虑研究对象属性的空间相关性和 位于相邻样品间 变异性，同时插值半径的确定常凭经验，没有确 位于同一样品中 切的方法.四面体剖分法的插值精度较高，可逐 位于下边界 层细化；但随细化的深人，计算及处理时间会成 位于间断的样品间 倍的增加，即精度的提高是以绘制时间为代价； 同时该方法的数据结构描述复杂，也没有考虑 研究对象属性的空间相关性和变异性.克立格 法的突出优点是充分考虑了属性的空间相关性 图2空间点属性确定示意图 和变异性，易于构造规则体数据；但该方法整 Fig.2 Determination of attributes at any point along a 个插值过程比较繁琐，计算量大，效率低 drillhole 通过对插值方法的分析，同时考虑系统的 若所求点刚好位于一个样品中间，则属性 硬件实现环境，拟采用计算方便的Shepard插值 同该样品；若所求点刚好位于两相邻样的交界 方法.但对于插值半径的确定充分考虑了矿石 处，则该点的属性与后一个样品的属性相同；若 品位的空间相关性和变异性的影响，从而克服 所求点位于两间断的样品之间，则该点的品位 了该方法自身的缺点.计算步骤为： 为0，即为岩石；若所求点位于下边界（样品间断 (1)对研究范围内所得空间点的品位值进行 处或孔底)，则该点的属性同该样品 统计，建立品位频度直方图（见图3），分析品位 通过上述原始数据的转化方法，共得到空 的随机性分布.对所得的空间点数据按不同的 间，点数据9159个，但这些空间点数据的分布仍 品位区段（就本文的实例而言，按品位介于24% 呈现散乱、不规则的体数据形式，有待于进一步 26%,26%一28%…，68%-70%等不同区段)分别 处理才能实现矿床的体视化 统计各区段的出现频率 (2)对研究范围内所得空间点数据的矿石品 2 三维规则体数据的生成 位的变异性进行分析，建立TFe品位在垂直和 由于各种勘测数据一般只能是对三维连续 地魂惠信 数据场的有限采样，不论是原始的采样数据还 是通过其转化所得的空间点数据，都具有明显 的离散特征：同时由于地质体自身的大尺度的 特性，又具有明显的稀疏特征.为了实现矿体的 真三维显示，需要通过数据的空间插值处理，形 成均匀网格结构化数据或规则网格结构化数 据.这样处理得到的体数据可表示为： fx,y,z)=fi△x,jAy,k△z) (2) 图3品位频率直方图 其中，，，k分别为任一个采样点所在的行、列、 Fig.3 Frequency diagram of ore grades

北 京 科 技 大 学 学 报 2 0 0 2 年 第 6 期 为 l , 方位 角为a , 倾角为刀 , 此深度对应折线 的起 点坐 标为x(0 , y0 , 局) , 起点处的孔深为10 , 钻孔在此 深度 的相对孔深为 lr, ’l = l一 10 , l 护` = ’.l co 叨 . 则此点 的坐标 x( ,y 力为 : 卜 一 x0丫 c ?as 一 局月 “ osf aso } , 一 y0 材 『 {!ln 一 +y0 ` “ osP ’ 叫 夕= 几+l sl 丫 ( l ) 这样 , 对所有钻孔均从孔 口 水平 (就 实例 而 言 , 为一 2 8 m )开始 , 按垂直深度 l m 由样 品数 据求得 相应 的空 间点坐标 . 而对应 的各 空 间点 的属性 (品位和 岩性 )值的确定方法见 图 2 所示 . ) 位于相邻样 品间 位 于同一样品 中 位于下边界 位于间断的样 品间 二曰鬓翼馨1 图 2 空间点属 性确 定示 意图 iF .g Z D e et 彻i n a it o n Of a t r i b u et s a t a n y pio n t a l o n g a d r il lk o l e 若所求点 刚好位 于一个样 品中间 , 则属 性 同该样 品 ; 若 所求点 刚好位 于两 相邻样 的交界 处 , 则该点 的属性 与后一个样 品 的属性 相同 ; 若 所求点位 于两间断的样品之 间 , 则 该点的 品位 为 0 , 即为岩石 ; 若所求点位于下边界(样 品 间断 处或孔底 ) , 则该点 的属性 同该 样品 . 通过上述原始数据 的转化方法 , 共得到空 间点数据 9 159 个 , 但这些空间点数 据的分布仍 呈现散乱 、 不规则 的体数据形式 , 有待于进一步 处理才 能实现矿床 的体视化 . 层 的 序号 ; 山 , 妙 , 山分别 为 网格点 在三个方 向 上的 间距 , 三者之间未必相 等 , 当相等时为均匀 网格结构化数据 , 否 则为规则网格结构化数据 ; fx( , y , z) 为体数据 在网格点i(, j , k) 处 的函 数值 . 数据 的空 间插值是地矿工程计算机体视化 过程 中必不 可少 的一项工作 , 插值算法 的精度 和 效率将直 接影响地矿 工程的三维显示效果 . 通过对现有插值方法 的分析 , 可 以用来形成 地 矿体空 间体 数据的插值 方法主要有 hs eP ar d 法 、 四面体剖分法和克立格 法 . hs eP ar d 法 的优点是计算方便 , 易于 理解 ; 缺点是 没有考虑研究对象属 性的空 间相关性和 变异性 , 同时插值半径 的确定常凭经验 , 没有确 切的 方法 . 四面体剖分法的插值精度较 高 , 可 逐 层细 化 ; 但 随细化的深人 , 计算及处理时间会成 倍 的增加 , 即精度的提高是 以绘制时间为代价 ; 同时该方法 的数据 结构描述复杂 , 也没有考虑 研究对象属 性的空间相关性和变异性 . 克立格 法 的突 出优 点是充分考虑 了属性 的空 间相关性 和 变异性 , 易于构造规则体数 据 ; 但该方法整 个插值过程 比较 繁琐 , 计算量大 , 效率低 . 通过对插 值方法的 分析 , 同时考虑系统 的 硬件实现环境 , 拟采用计算方便 的 S h叩ar d 插值 方法 . 但对于插值半径 的确 定充分考虑 了矿石 品位 的空 间相关 性和变异性 的影 响 , 从而克服 了该 方法 自身的缺点 . 计算步骤 为 : ( l) 对研究范 围内所得空 间点 的 品位值进行 统计 , 建立 品位频度直方 图 (见 图 3) , 分析 品位 的随机性分布s[] . 对所得的空 间点数据按不 同的 品位区 段(就本文 的实例而言 , 按品位介于 24 % 一 26 % , 26 % 一28 % , … , 68 % 一 70 % 等不 同区段 )分别 统计各 区 段的 出现频率 . 但)对研究范 围内所得空 间点数据的 矿石 品 位的变异 性进行分析 , 建立 TF e 品位在垂直 和 硒哪脚娜娜娜斌娜脚雄娜娜舞娜黝卿卿卿脚寒稗)爵摹攀馨缪弊辫鬓拳 ~ 娜粼娜攀舞麟 篇黑豁蕊麟咖翔剐即70D 糊删剐姗.20 均健.匀”,吸阶叼俪朋晰姗却殉知10栩曲比,lJ 。 ùùWoMóLdMet 2 三维规则体数据的生成 由于各种勘测数据一般 只 能是对三维连续 数据场 的 有限采样 , 不论 是原始 的采样数据还 是通过其转 化所得 的空 间点数据 , 都具有 明显 的离散特 征 ; 同时 由于地 质体 自身 的大尺度 的 特性 , 又具有 明显的 稀疏 特征 . 为 了实现矿 体的 真三维显示 , 需要通过数据的空 间插值处理 , 形 成 均匀 网 格结 构 化数 据或 规则 网格 结构 化数 据 . 这样处理 得到 的体数据 可表示 为 : 刀元y, 习 = f( i公 ,j 匆沂山 ) (2 ) 其 中 , i,j , k分别 为任一个采样点所在 的行 、 列 、 念 图 3 品位频 率直 方图 F ig ·飞 F er q u e n cy d i a g r a m o f o er g r a d e s

VoL24 李翠平等：刊矿床体视化的一种空间规则体数据生成方法 591· 水平方向上的实验变异函数曲线（见图4）.垂直 结合勘探网度的大小，确定空间插值半径（就本 方向上的变异性由5个方向（沿勘探线剖面上 文的实例而言，空间插值半径为15m). 的30°，60°，90°，120°和150)上TFe品位的变异 (4)利用Shepard法进行空间插值，形成体素 函数及其平均变异函数曲线来确定；水平方向 大小为10m×10m×1m的三维规则体数据（就本 上的变异性由x,y两水平方向上TFe品位的变 文的实例而言，共包含体素105350个，经数据 异函数及其平均变异函数曲线来确定， 压缩保存在数据库中的体素为23397个) (3)根据TFe品位的变异性分析结果，同时 (a) (b) 图4变异函数曲线，(a)垂直方向，()水平方向 Fig.4 Variogram of ore grades 3定义空间体数据的数据结构 x=nodes[n].i*xStep +xStartCord; y=nodes[n]j*yStep+yStartCord; 对于规则体数据fx,y,z)=fi△x,j△y,k△z)的 z=nodes[n].k*zStep zStartCord; 表示，无需给出各采样点的空间位置，只要给出 该网格点品位值为：nodes[n].Grade.. 起始点的坐标、空间任意网格点的序号及x,y,z 三个方向上的增量，即可确定该点的空间位置. 4 结论 这种表示方式既方便数据的存取，又节省存储 数据的预处理是体视化技术的一个重要方 空间.系统中体数据的数据结构形式为： 面，不同的体视化对象和相应的数据来源，要求 const int NUMBER=10000; 对原始数据采用不同的预处理方法.本文面向 const double xStartCord =3355700.0; 地矿工程体视化技术开发与系统建立，提出了 const double yStartCord=87850.0; const double zStartCord =-288.0; 一种对矿床地质勘测数据进行预处理的方法， 包括首先由原始样品数据转化为空间点数据， const double xStep =10.0; 再通过空间的插值处理生成系统所需的规则体 const double yStep=10.0; 数据等两个基本步骤.本方法充分考虑了地矿 const double zStep=1.0; 工程体视化对象及数据取样特征、数据空间相 struct tagNode 关性及变异性等特点，且结合某地下铁矿的实 int i; ∥该网格点在x轴上的序号 际数据将方法予以实现，验证了方法的可行性 intj; ∥该网格点在y轴上的序号 和有效性. int k; ∥该网格点在z轴上的序号 double Grade; ∥该网格点品位值 参考文献 double xNextGrade;I∥沿x轴的下一点的品位值 1毛善君，许友志，张海荣，等.空间地质模型及其可视 double yNextGrade;∥沿y轴的下一点的品位值 化系统[A].中国数据地质M小.北京：地质出版社， double zNextGrade;∥沿z轴的下一点的品位值 1996.186 nodes[NUMBER]; 2胡金星，吴立新，高卫贞，等.三维地学模拟体视化 则任一网格点n的实际坐标为： 技术的应用研究.煤炭学报，1999,244：345

4 4 ) 水 平方向 上的实验 变异函数 李翠平等曲线见 矿床体视图垂 化的一种空直结 间规则体合勘探 数据生成网度的大 方法小确定 空间插 值半径(就 二本 9 ( L o 1 V 2 5 方的 向上的306 变异性由9012 5个方0和1 向(沿勘)上TFe 探线剖面品位的变 上文异 的实例(4)利 而言空用Shepard 间插值半法进行 径为15空间插值 m).形成体 素 函上异 数及其的变异函 平均变异性由x,y其平均变 函数曲两水平方异函数 线来确定向上Te曲线来确 ;水平方F品位的定 向大变文压 小为ro的实例缩保存 m10而言共在数据库 lm的三包含体素中 维规则1053素为2 体数据(就个经97) 本据 (3)根 据TFe品 位的变 异性分析 结果同 时 瑟万图4变 异函数曲F啥4Va …瑟线(a)垂直rigm 方向伪fregd )水平方向 引 , 一 一~ 二`~ `出“二 ~一 x =nodes fnlixste 。+xstar Cord: 表 -对于示无需 一规则体数给出各 拒封卜y采样点的 z)=f(i酞空间位置 一,j’勿kAz)只要给 y的z出该 一nodes=网格点 [n].j*ysteP[n]k*zst品位值 +ys以ePztar为nodes[ cord;Cn].Grade 起三 始点的个方向 坐标空上的增量 间任意网即可确 格点的定该 序号及x,空间位 yz置4 结论 汀do ubleGro。xN ade;xtGrade; 刀该网格/沿x轴 点品位的下一 值点的品位 …值 辫参考文1毛善君 献许友志 张海荣筹 空间地质 模型及其 可视 d} oubleyN水nods协 xtGrade[UMBER] ;刀沿夕轴/ 的下一 点的品位 值 化系统19%2胡金星 [A]中国86昊立新 数据地质高卫贞等 [M]北京三维地 地质出版学模拟体 社视化 则 任一网 格点,的 实际坐标 为 技术的 应用研究 [J]煤炭学 报192 4():35

·592 北京科技大学学报 2002年第6期 3李杰，毛善君，张海荣，等.层状矿床的三维仿真与6李仲学，李翠平.矿床仿真及体视化技术).计算机 动态显示).中国矿业大学学报，1999,28(4)：335 仿真，2000,17(5)6 4朱大培，牛文杰，杨钦，等.地质构造的三维可视化7李翠平，李仲学，赵文广矿床的体视化及仿真系统 [.北京航空航天大学学报，2001,27(4)：448 框架[).金属矿山，2001(10)：44 5黄文静，唐龙，唐泽圣.体绘制及三维交互技术在地 8侯景儒，郭光裕，矿床统计预测及地质统计学的理 质数据可视化中的应用)，工程图学学报，1998(3)：60 论与应用M.北京：冶金工业出版社，1993 A Method for Spatialization and Regularization of Geological Data in the Volume Visualization of Mineral Deposits LI Cuiping,LI Zhongxue,HU Nailian Mining Research Institute,UST Beijing,Beijing 100083,China ABSTRACT A methodology is presented for the spatialization and regularization of geological data in de- veloping a volume visualization system of mineral deposits.In the process of spatialization,drill holes are treated as a series of broken !ines with the spatial direction,and the grade at a point is obtained through layering of the deposit with one meter depth.The regularization of geological data is realized using the Shepard method with the regionalization of mineralization taken into account through geostatistical variograms.A case of underground iron ore mine is studied to verify and validate the methodology. KEY WORDS volume visualization;data spatialization and regularization:Shepard method (上接第584页) Outburst-prone Tectonophisical Environment and Its Applications GUO Deyong",HAN Dexin',Wang Xinyi 1)China University of Mining Technology (Beijing Campus),Beijing 100083,China 2)Pingdingshan Coal (Group)Co Ltd,Pingdingshan 467000,China ABSTRACT Based on the research of the control mechanism of geology structure for the outburst and the distribution rule,the theory and method of tectonophisics were applied to the forecast of outburst.It was pro- posed that the outburst-prone tectonophisical environment is composed of 4 factors,that is tectonic combina- tion characteristic,tectonic stress field,and tectonic coal and gas.The control mechanism of all factors of the outburst-prone tectonophisical environment was discussed.It was shown that the synthetic action of multi fac- tors of the outburst-prone tectonophisical environment controls the outburst in the geology tectonic area,and this provides the theoretic basement for forecasting the outburst proneness due to the geology structure.En- gineering practice proves that the outburst of coal and gas can be forecasted accurately with the analysis of the outburst-prone tectonophisical environment. KEY WORDS outburst of coal and gas;tectonophysical environment;outburst forecast

9北 京2 科 技 大 学 学 报 5 年 第 期 2 2 00 6 李 杰 , 毛 善君 , 张海荣 , 等 . 层状 矿床的三维 仿真 与 动态显示 [J ] . 中 国矿业 大学 学报 , 19 9 9 , 2 8 ( 4 ) : 3 3 5 朱大培 , 牛文杰 , 杨钦 , 等 . 地 质构造 的三维 可视 化 【J ] . 北京航 空航 天大学 学报 , 2 0 0 1 , 2 7 ( 4 ) : 4 4 8 黄文静 , 唐 龙 , 唐泽圣 . 体绘制 及三维 交互 技术在 地 质 数据可视化 中的应用 明 . 工程 图学 学报 , 1 9 9 8( 3) : 60 6 李 仲学 , 李翠 平 . 矿 床仿 真及体视化技 术 [J] . 计算机 仿真 , 2 00 0 , 1 7 ( 5 ) : 6 7 李翠平 , 李 仲学 , 赵 文广 . 矿床的体视 化及 仿真系统 框架 I J ] . 金属 矿 山 , 2 0 0 1 ( 1 0 ) : 4 4 8 侯景 儒 , 郭光裕 . 矿 床统计预 测及 地质统计学 的理 论与应 用 【M ] . 北京 : 冶金 工业 出版社 , 19 93 A M e ht o d fo r S P at i a li z at i o n an d R e g u l ar i z at i o n o f G e o l o g i c a l D at a in ht e Vo l u 工n e Vi s ua li z at i o n o f M ine r a l D e P o s it s LI uC iP i 刀g, IL hZ o n g x ue, H U N改il ia n M i n i n g eR s e ar hc ls t i t u t e , U S T B e ij l n g , B e ij i n g 10 0 0 8 3 , C l l i n a A B S T R A C T A m e ht o d o l o gy 1 5 P r e s e in e d for ht e s Piat a li z at i o n an d r e gu l iar z at i o n o f g e o l o g i c a l d at a in d e - v e l o Pi n g a v o 】u m e v i s u a li z at i o n sy s t e m o f m i n e r a l d e P o s it s . I n ht e P r o e e s s o f s P at i a li z at ion , dr ill h o l e s aer tr e a t e d a s a s e r i e s o f br o ke n lin e s w iht ht e s Pat i a l d i r e e it o n , an d ht e gr a d e at a P o int i s o b t a i n e d htr o u g h lay e r i n g o f ht e d eP o s it w iht o n e m e t e r d e Pht . The r e g u l ar i z at i o n o f g e o l o g i c a l d at a 1 5 r e a li z e d u s i n g ht e s h e P ar d m e ht o d w iht t h e r e g i o n a li z at i o n o f m i n e r a li z at i o n tak e n int o a e c o u n t htr o u g h g e o s t at i s t i e a l v iar o gr am s . A e a s e o f un d e gr o un d i r o n o er m i n e 1 5 s ut d i e d t o v e r i fy an d v a lid at e ht e m e ht o d o l o gy . K E Y W O R D S v o l um e v i s u a li z at i o n : d at a s Pat i a li z at i o n an d r e gu l ar i z at i o n ; Sh e Par d m e ht o d 拭材派洲C餐来今芳孕冷快来来弓{ ( 来 来来来拭(令冷矢谓e 矢舞弓装来白卜 (上接 第 5 84 页 ) O u ht ur s t 一 P r o n e eT e t o n o P h i s i c a l E n v i r o n m e n t an d It s AP P li e at i o n s G U O D ey o 心 ’ , 儿魂N D ex in] ’ , 肠gn iX yn 扩’ l ) C h in a U in v e rs iyt o f M in l ll g & eT e hn o l o gy (B e ij in g C田 n P u s ) , B e ij i n g 10 0 0 83 , C h in a 2 ) P in g d i n g sh an C o a l ( G or uP ) C o L dt , P I ll g di gn sh an 4 6 7 0 0 0 , C h i n a A B S T R A C T B a s e d o n ht e r e s e ar e h o f ht e e o ntr o l m e c h an i s m o f g e o l o gy s trU e trU e of r ht e o u t b ur s t an d 小e d i s tr ib ut i o n ru l e , ht e ht e o yr an d m e ht o d o f t e c ot n o Ph i s i c s w e r e a P P1i e d t o ht e fo er e a s t o f o ut b ur s t . l w a s P r o - Po s e d ht a t ht e o u t b ur s t . P r o n e et e t o n o hP i s i e al e n v i r o mn e nt 1 5 e o m P o s e d o f 4 fa c t o r s , ht at 1 5 t e c t o n i e c o n i b in a ` t ion e h ar e t e r i s t i c , t e e t o 正e s etr s s if e ld , an d t e e t o n i e e o a l a n d g a s . T h e e o ntr o l m e c h an i s m o f a ll fa c t o r s o f ht e o u t b ur s t 一 P r o n e t e e t o n op hi s i e a l e n v i r o mn e nt w a s d i s e u s s e d . It w a s s h o w n ht a t ht e s y n th e t i e a e t i o n o f mu lt i fa e - t o r s o f ht e o ut b ur s t 一 P r o n e t e e t o on Ph i s i e a l e vn i r o mn e in e o ntr o l s ht e o ubt ur s t i n ht e g e o l o gy t e c t on i e ar e a , an d ht i s P r o v ide s ht e ht e o r e t i c b a s e m e nt for fo r e c a s t i n g ht e o u t b u r s t P r o n e n e s s du e t o ht e g e o l o gy s utr e trU e . 劲 - g i n e e r i n g P r a c t i e e Por v e s ht at ht e o u ht ur s t o f c o a l an d g a s e an b e fo r e e a s t e d a e e ur at e l y w it h ht e an a ly s i s o f ht e o u t b ur s t 一 rP o n e t e e t o n op h i s i e a l e vn ir o nr e nt . K E Y W O R D S o Ll t b ur s t o f e o a l a xl d g as : et c t o n o P场 s i c a l e vn i r o mn e nt : o t l t b ur s t for e c a s t