可视化储量计算普通克里格法研究.pdfVIP

中国地球物理2006 可视化储量计算的普通克里格法研究 揣媛媛范继璋肖克炎 (吉林大学综音信息矿产预测研究所吉林长春130026) 1．引言。 可视化储量计算的优势在于采用人机交互的方式智能编辑剖面，动态地圈定矿体，将生成的矿体 细分为晶胞，与普通克里格计算所划分的三维网格体一一对应，直观、准确地计算出矿体每个晶胞的 开发三维可视化储量计算模块，该模块应用在吉林省西岔金矿中取得了较理想的预测结果。 2．地质建模及可视化 (1)数据准备 数据库文件包含三个数据表文件，即①钻孔位置表，包括钻孔在三维空间的起点坐标(x，Yz)以 及钻孔的长度；②钻孔形态表。包括钻孔编号、孔段范围、钻孔在空间的倾斜方向和倾角；③采样信 息表t包括孔段起始、孔段终止、样品编号、岩性及含矿信息等。 (2)剖面定义和曲面生成 根据已有的钻孔资料定义剖面文件，描述岩性、剖面边框、钻孔段投影线、地层界线、断层线等 信息，经过投影变换和坐标变换生成三维钻孔立体图。根据地质轮廓线的表面建模法，采用遗传算法 来生成两轮廓线间的三角面片，再由面片组合成空间曲面，形成矿体的线框模型。 (3)实体生成 由编辑好的剖面文件自动生成三维地质体，再经光照模型、裁剪处理等计算机图形学算法将三维 地质体逼真地显示出来。． 3．普通克里格法储量计算 (13计算实验变差函数 主要采用半变差函数进行计算，规则网数据实验变差函数计算算法： ①要计算变差函数的方向．可能是多个，假设为m个。②要计算的步数n。⑧变化性测度类型， 并初始化为0。⑤确定数据中的一点，对参数给定的每个方向找步数为1~n的点，与给定点组成点对， 所要计算的变差函数。 (2)储量估算 地质体表示为一系列固定形状晶胞(三角形面片组成的封闭体)的集合，普通克里格法实质也是 把地质体划分为规则的晶胞单元，这样，三维建模生成的晶胞与克里格划分的晶胞一一对应起来，估 计地层内部的地质属性或矿体内部的品位值，就能准确地估算地质体的储量。 整个矿体的储量计算方法是先计算矿物在每个晶胞内的平均品位，采用晶胞储量计算公式(品位 ×比重×体积÷100)算出每个晶胞的储量，最后把每个晶胞的储量相加得到总储量。其中，计算平均 品位是最关键的步骤，其算法描述如下： ①初始化组合样数组CA，一个组合样表示为一个空间点和一个品位值，初始化块段数组BA，初 始化搜索距离D。②如果块段数组BA中的每个块段都被处理过，转④，否则在块段数组BA中取下 一个块段。③计算块段的中心点坐标。计算到该中心点的距离小于D的组合样，将这些组合样到块段 中心点距离的平方存入临时距离数组DA，将组合样相应的品位值存入一个临时品位值数组QA。如果 DA中有某个值为0，则把相应的品位值存入块段中．否则按距离平方反比加权方法计算块段的平均品 位，将结累存入块段中，转②。④返回块段数组BA，结束。 拳研究由国家863项目《地学空间三维智能分析与可视化技术》(2002AAl35160)资助。 参考文献 王家华等．克里金地质绘图技术北京：石油_[业出版社，1999 ·188·