储量计算自动化,提高地质工作水平

jokonem

储量计算自动化,提高地质工作水平

 

张文宽  徐洪恩  谯小平  杨本锦

（四川省地质矿产局403地质队，四川峨眉 614200）

摘要：文章介绍的《储量计算自动化程序》是作者90年代DOS版程序的全面升级。该程序按照目前地质工作普遍采用的公式和方法，符合现行规范规定和要求。程序对储量计算中矿体圈定、公式判别、数据输入、对矿体有不同主张的表示方法等难点解决较好。该程序已经多个矿区应用检验，最近又得南非英铂公司采用推广。

关键词：固体矿产 储量计算 矿体圈定 应用程序

 

十多年前，计算机操作系统还是DOS的年代，我队进行超大型呷村银多金属矿床详勘时，由于矿体多，样品多，传统的手工储量计算已不适应要求，我们编写了一个DOS版的《储量计算自动化程序》,为呷村矿区详勘报告保质按时提交，发挥了重要作用。该报告于1995年6月获中华人民共和国地质矿产部地质报告成果奖。DOS版的《储量计算自动化程序》获中华人民共和国地质矿产部科技进步奖。

十多年了，DOS版《储量计算自动化程序》虽经几次升级，但总赶不上电脑的飞速发展、时代的进步。2005年我队因为夏塞铅锌银矿、淌塘铜矿、拉拉铜矿、峨水石灰岩矿的工作需要，对DOS版《储量计算自动化程序》进行全面升级。升级后操作系统为Windows XP，编程语言用Visual Basic6.0（可视化编程工具），数据库用Microsoft Access2000（微软办公系统数据库）。储量计算方法完全按照有关地质勘探规范规定。升级后的储量计算自动化程序，只需将矿床储量计算有关原始数据，如样品编号、取样长度、轴夹角（或矿体产状）、分析结果、矿体块段面积（或剖面面积）、体重等按程序表格样式输入电脑。再输入工业指标后，按电脑提示点击按钮，即可完成工程中矿体圈定与厚度、平均品位计算、块段储量计算和各类型、各级别储量的合计、总计。现将升级的有关情况介绍如下：

一、数据输入可按电脑提示手动输入，也可从电脑其它文件中复制粘贴进入样品原始数据表和储量计算原始数据表。样品原始数据表包含：序号(rd)、工程编号（ZK）、样品编号（NO）、样品长度（L）、样品品位（C）、轴夹角（A）、矿体倾向（AL）、矿体倾角（A2）、地面坡向（B1）、地面坡角（B2）等字段。只要按此输入数据，电脑会根据岩心轴夹角数据和工程代号自动识别是钻孔样品还是地表工程样品，岩心轴夹角项为“空”是地表样品，反之为钻孔样品。电脑自动分别采用不同的真厚度计算公式；

钻孔样品真厚度计算公式：M=L×Sin（A）

地表样品真厚度计算公式：

M=L{Sin（A2）×Cos（B2）×Sin（γ）±Cos（B2）×Cos（A2）}

式中±号按地面坡向与矿体倾向相反时为正。地面坡向与矿体倾向一致，并且矿体倾角小于地面坡角时为负。

γ——样品方向与岩层走向的夹角

工程中矿体加权平均品位计算公式：

PC=Σ{（M1×C1）+（M2×C2）+（M3×C3）……}÷Σ（M1+M2+M3……）

二、数据输入按地质工作者确定的顺序录入。由于补采样品等原因，常常需要向数据表添加样品或调整样品位置。新数据添放在表的最后，在序号字段“rd”栏里写正确位置的序号。程序计算时电脑会自动采用升序排列，将样品归入正确位置。序号可写成小数形式，例：要在2、3之间加入样品，则新加入样品“rd”序号可写成2.1、2.2……；若还需在2.1、2.2之间加样品，则新加入样品可写成2.1.1、2.1.2……，以此类推。

三、矿体圈定涉及边界品位、块段平均品位、最小可采厚度、最大夹石剔除厚度等多种因素，人脑圈定亦是斟酌再三。目前发表的储量计算程序多绕道而行，人工圈定后由电脑执行计算。这样降低了自动化程度，失去使用电脑的意义。我们的矿体圈定，发挥电脑功能，逐步筛选，逐步合并优化，得出较满意的结果。经多个矿区的手动圈定与电脑圈定对比，两者符合率在97%以上，分歧之点亦是各自的侧重点不同。地质人员间由于对矿床成因、类型认识有差异，对矿体圈定可能有不同的主张。电脑自动圈定矿体后，为了反映不同看法，程序中增加了手动修改一节。修改后的数据，电脑会自动读入，并据此计算矿体平均品位、平均厚度。多种圈矿方案择善而行，可以提高对矿床的认识。

四、程序设计了不同工业指标圈定矿体及计算储量的功能。输入新的工业指标后，点击电脑即可获得新指标计算的储量。这对矿床经济评价，科学合理的开发矿床很有意义。

五、程序完全按照目前地质工作普遍采用的公式和做法，符合现行规范规定和要求。使用电脑进行储量计算，除有提高精度，加快速度的好处外。还有如下优点：

1、动态计算，方便修改。储量计算完成后。审查时若有数据需要改动。牵一发而动全身。手动重算，耗时费力。电脑储量计算自动化，数据相互关联，改动数据后，电脑会对涉及项目自动重算改正。

2、使用不同工业指标，自动圈定矿体，加上手工修改，形成多种方案对比，加深了矿床认识，有利于科学决策。

3、电脑储量计算客观公正。矿体厚、品位变化大的矿区，矿体圈定涉及因素多。那个样品该圈入矿体，那个样品可做夹石剔除，手动圈定时常有操作者主观意愿影响空间。不同的人计算，会有不同的矿床平均品位和储量，这样会给有心人以可乘之机。电脑圈定，不因人而异。在矿床储量有偿使用的今天，客观公正有其特别意义。

本程序除在我队地质大调查项目和社会勘探项目中应用取得较好效果外。2006年4月17日，四川省地矿局和南非英铂公司在峨眉山召开的地质工作数字化培训班上讲授推广获得好评。

在软件编制和本文写作过程中，承西南科研所张建龙老师、我局文辉老师以及杨昌银、曾凡勇、刘云霞、万霞、刘丽琴等同事热心支持和帮助，特此致谢！

我队除编制储量计算软件外，还利用Visual Basic6.0和AutoCAD 2004编制了《钻孔、探槽数字化地质编录、钻孔柱状图和探槽素描图自动绘制与输出、表格打印、剖面上钻孔投影作图软件》。《滑坡计算软件》等，欢迎有兴趣的同行来人（函）索取。

 

[第一作者简介] 张文宽，（1954-），男，1978年毕业于成都地质学院，高级工程师，从事地质勘探与技术管理工作。

 

xinghesnak

楼主，发这个东东什么意思？

jokonem

现在很多人在找矿，他们都希望有这类型的软件来提高工作效率，所以我在这里抛砖引玉。

jokonem

我们自已在进行这方面的编程工作，盼有志之士一起来完成这项工作，当然我们是业余的，不是专业的，所以技术上不是很好，不敢和你们专业的编程人员相比。我们这帮子人都只懂一点CAD,WORD、ACCESS及VB6的皮毛，还望不吝赐教。

 

zhang007

早知道有这么一个东西我也省一点时间了，前一段时间我也编写了一个类似的东西，是用来自动读取MAPGIS图中指定矿体面积，按控制矿体的工程加权来计算单个矿体平均品位、同一剖面上的同一类矿体按面积加权计算平均品位及总面积，只是这类软件的使用太有限了，只适合比较复杂的大型矿床，在实际工作中往往都是小买卖，手工做也是飞快的。

另外，我认为储量计算的难点不在这里，而在于将图纸信息关联到数据库中相对应的工程品位中去，有时间我们可以讨论一下.

zoubo604

我想问一下各位:文章中的一个公式

地表样品真厚度计算公式：M=L{Sin（A2）×Cos（B2）×Sin（γ）±Cos（B2）×Cos（A2）}

           我在地质技术标准中看到是的这样的,前面3项与上相同,后面是±Cos（A2）×Sin（B2）.是公式变形吗.是怎么变成这个形式的,请大家把变形的过程给我写一下吧,谢谢了.

附:地质技术标准DD2006-01中的地表样品真厚度计算公式：M=L{Sin（A2）×Cos（B2）×Sin（γ）±COS（A2）×Sin（B2）}