excelvba在矿山测量中的应用.pdf

I S S N 1 6 712 90 0 CN 431 3 4 7／ TD 采矿技术第 7卷第 4期 Mi ni ng Te c hn o l o g y，Vo 1 ． 7， No． 4 2 0 0 7年 1 2月 De C． 2 00 7 E x c e l VB A 在 矿 山 测 量 中 的 应 用 褚建民 辰洲矿业股份有限公 司， 湖南 怀化市4 1 9 6 0 7 摘要 以矿山测量工作 中经常进行的井下经纬仪 导线测量计算为实例 ， 详细介绍 了‘ E x - c e l V B A的使 用方法及其优越性 ， 并列 出该计算全部有效的程序编码 。 关键词 E x c e l ； V B A编程； 矿山测量 ； 井下导线 用 E x c e l 进行数据处理 ， 工作直观方 便 ， ， 许 多 人都乐于使用。然而， 其鲜为人知 或者说是有待 深度开发 的另一强大功能还在于它提供 的一种新 兴的程序设计语言一一宏语言 V B A。 笔者性情所致 ， 一度潜心 于 E x c e l V B A在矿 山 测量中的应用实践 ， 经不断进行调试与排错后 ， 最终 如获所愿， 成功开发出了 地面控制导线测量计算 与成果整理程序 、 井下经纬仪导线测量计算与成 果整理程序 和 井下罗盘仪导线测量计算 与成果 整理程序 等一系列矿 山测量计算程序。现不揣浅 陋， 试以矿山测量工作 中经常进行 的井下经纬仪导 线测量计算为实例， 将编程精要予以和盘托出， 以期 有抛砖引玉之效 。 1 用窗体录入测量观测数据 众所周知 ， E x c e l 是一个 电子制作表格 的软件 。 直接利用它来输人大量 的测量观测数据 ， 就会发现 很容易 出错 ， 因为在输入过程中可能会被纵横交错 的大量测 量 观测 数 据弄 得 眼花 嘹 乱 而导 至 张 冠 李戴。 V B A给了 E x c e l 很大 的支持 。它 能让 E x c e l 变 得如同魔法一般任意调用任何想要的控件来增强它 本身的功能。在 E x c e l中可用 V B A建立一个 测量 观测数据输入窗体 ， 通过它来输入每一条导线的测 量观测数据 ， 如测线名称 测站点 一测量 点 、 水平 角、 垂直角 、 斜仪高 、 觇高 的录入均从此 窗 口进行 输入和确认 ； 然后再通过模块 安全而快捷地将窗体 中的测量 观测数据交换 到 E x c e l 工作簿 中的 井下 经纬仪导线测量计算 工作表的单元格中 见图 1 。 显然 ， 用 V B A来实现输入数据要 比手工输入方 便、 可靠， 而且简捷许多。 使用窗体录入测量观测数据的具体操作步骤。 图 1 观测数据录入窗体 1 在 E x c e l 下打开一个新的工作表 ， 单击 E x c e l 主菜单上的“ 工具 宏 V i s u a l B a s i c编辑器” 选 项打开 Vi s u a l B a s i c编辑器 ， 再单击该编辑器主菜单 中的“ 视 图 工程资源管理器 ” 选项 以打开“ 工程资 源管理器” 。 2 在显示的树形图上右击 ， 单击右键菜单上 的“ 插入 窗体” 选项 以插入一个 “ 窗体” 。然后按 同样 的方式插入一个“ 模块 ” 以交换数据 ， 并且更改 这两个新建的资源的属性 。 ． 3 在添加控件前 ， 应对使用控件 的摆放位置 有一个大概的构思。窗体 中各控件的名字为各 自默 认的名字。 4 双击新建的模块 ， 编辑新建模块 中的代码 Opt i o n Ex pl i c i t Pu b l i c C X As S t r i n g P u b l i c s p j A s S t r i n g P u b l i c c z j A s S t r i n g P ub l i c x p As S t r i ng P ub l i c x y g As St r i n g P u b l i c s g As S t r i n g 采矿技术 2 0 0 7 ， 7 4 S u b I n i t ex『 1『 1 s p j ” ” c z j ” ” x p”” x y g ”” s g ”” En d S u b S u b I n p u t V al u e Di m Te mp Ms g Bo x As Vb Ms g Bo x Re s u l t Di m Co un t As I n t e g e r I f T ri m U s e r F o r m1 ． T e x t B o x 1 ． T e x t ” ” O r T ri m U s e r F o r m1 ． T e x t B o x 2 ． T e x t ” ” O r T ri m U s e r F o r m1 ． T e x t B o x 3 ． T e x t ” ” O r T ri m U s e r F o r m1 ． T e x t B o x 4 ． T e x t ” ” O r T ri m U s e r F o r m1 ． T e x t B o x 5 ． T e x t ” ” O r T ri m U s e r F o r m1 ． T e x t B o x 6 ． T e x t ” ” T e m p Ms g B o x M s -g S O X” 请 输 入 完 整 的观 测 数 据 ” ， v b O K O n l y ， ” 提示” En d En d I f Co u n t 0 R a n g e ” A 1 ” ． S e l e c t D o U n t i l S e l e c t i o n ． O f f s e t C o u n t ， 0 ． V a l u e ” ” I f T ri m U s e r F o r m1 ． T e x t B o x 1 ． T e x t S e l e c t i o n ． Of f - s e t C o u n t ， 0 ． V a l u e T h e n T e m p Ms g B o x Ms g B o x” 本 测 线 已 录 人 ， v b O K O n l y ，” 提示” End En d I f Co u n t Co u nt 1 Lo o p C X T ri m U s e r F o r m1 ． T e x t B o x 1 ． T e x t s p jT ri m U s e r F o r m1 ． T e x t B o x 2 ． T e x t c z jT ri m U s e r F o r m1 ． T e x t B o x 3 ． T e x t x p T ri m U s e r F o r m1 ． T e x t B o x 4 ． T e x t x y g T ri m U s e r F o r m1 ． T e x t B o x 5 ． T e x t s g T ri m U s e r F o r m1 ． T e x t B o x 6 ． T e x t Co u n t 0 R a n g e ” A 1 ” ． S e l e c t D o U n t i l S e l e c t i o n ． O f f s e t C o u n t ， 0 ． V al u e ” ” Co u n t Co u nt 1 Lo o p S e l e c t i o n ． O f f s e t C o u n t ， 0 ． V al u e C X S e l e c t i o n ． O f f s e t C o u n t ，1 ． V al u e s p j S e l e c t i o n ． O f f s e t C o u n t ， 2 ． V a l u e c z j S e l e c t i o n ． O f f s e t C o u n t ， 3 ． V alu e x p S e l e c t i o n ． O f f s e t C o u n t ， 4 ． V al u e x y g S e l e c t i o n ． O f f s e t C o u n t ， 5 ． V a l u e s g C all I n i t If Ms g B o x ” 还要继续录入吗 ” ， v b O K C anc e l ，” 继续录 入 ” v b C a n c e l T h e n Us e r F o r m1． Hi de En d I f En d S u b 5 在表单上新建一个 “ 按钮 ” 控件并将名字 改为“ 确认” ， 在“ 设计模式” 下双击该按钮重新进入 “ V i s u a l B a s i c编辑该按钮单击响应代码如下 P ri v a t e S u b C o mm and B u t t o n l C l i c k Us e r Fo r m1 ． S h o w Ca l l I n i t En d S u b 6 保存 后退 出“ V i s u a l B a s i c编辑器 ” ， 单击 “ 输入观测数据” 按钮 ， 则可弹 出刚刚制作 的窗体。 输入观测数据进入窗体 ， 单击 “ 确认 ” 按钮 即可以将 每条导线边 测线 测量观测数据一次性导入“ 经纬 仪导线测量计算” 工作表中的单元格中。 2 使用 Ma s g B o x函数查询测量成果 M a s g B o x函数是 E x c e l V B A 的一个库 函数 ， 在 实际中使用非常广泛 ， 几乎在任何 V B A程序 中都有 出现。使用 M a s g B o x函数查询测量成果的具体步骤 如下。 1 单击 E x c e l 主菜单上的“ 工具 自定义” 选 项 ， 在弹出的“ 自定义” 对话框中勾选“ 控件工具箱” 前面的复选框即出现浮动的“ 控件工具箱” 。 2 在“ 控件工具箱” 中单击 “ 按钮” 控件 ， 在井 下经纬仪导线测量计算表的空 白区域新拖出一个按 钮控件并在“ 设计模式” 下在该控件下右击， 在弹出 的菜单上单击 “ 属性” 弹出“ 属性 ” 对话 框。将 “ 属 性” 对话框 中的 C a p t i o n属性 值改 为 “ 测量 成果查 询” 。 3 在“ 设计模式” 下双击“ 测量成果查询” 按 钮进入编辑该按钮单击响应代码的 S u b过程。在该 过程 中编辑如下代码 P r i v a t e S ub Co mma ndBu t t o nlC l i c k Di m Te mpRa g As Ra ng e Di m Te mpOf f s e t As I nt e g e r Di m Te mpMs g Bo x As VbMs g Bo x Re s ul t S e t Te mpRa g App l i c a t i on ． Se l e c t i o n I f f T e mp Ra g ． Ro w 1 O r T e mp Ra g ． R o w 2 Or T ri m T e m p R a g ． V a l u e ” ” T h e n 褚 建民 E x c e l V B A 在矿 山测量 中的应用 9 3 T e m p Ms g B o xM s g B o x ” 请用 鼠标选择测线所在 行” ， v b O K O n l y ，” 错误选择” En d End I f Ⅱ C e l l s 2 ， T e m p R a g ． C o l u m n ． V a l u e” 测线 测站 一 测点 ” T h e n T e mp O f f s e t 0 E l s e l f C e l l s 2 ，T e mp R a g ． C o l u mn ． V a l u e ” 水平 角” T h e n T e mp O f f s e t 一1 E l s e l f C e l l s 2 ，T e mp R a g ． C o l u mn ． V a l u e ” 垂直 角” The n T e mp O f f s e t 一2 E l s e I f C e l l s 2 ， T e m p R a g ． C o l u mn ． V a l u e” 斜边” The n Te mpOf f s e t 一3 E l s e I f C e l l s 2 ，T e mp R a g ． C o l u m n ． V a l u e ” 斜仪 高” T h e n Te mpOf f s e t 一4 E l s e I f C e l l s 2 ，T e mp R a g ． C o l u mn ． V a l u e” 视高” Th e n Te mpOf f s e t 一5 E l s e I f C e l l s 2 ， T e mp R a g ． C o l u mn ． V a l u e” 平距” Th e n T e mpOf f s e t 一6 End I f T e m p Ms g B o x Ms g B o x ” 测站 一 测点 ”一 ＆T e mp R a g ． O ff s e t ， T e m p O f f s e t ． V a l u e＆C h r 1 3 ＆C h r 1 3 ＆” 平 距 D ”一 ＆ T e m p R a g ． O f f s e t ， T e m p Off s e t7 ． V a l u e＆ C h r 1 3 ＆ ” 方位角 a ”一 ＆ T e m p R a g ． Off s e t ， T e m p O f f s e t8 ． V a l u e＆ C h r 1 3 ＆ C h r 1 3 ＆ ” 纵坐标 X ”一 ＆ T e mp R a g ． O ff s e t ， T e m p O ff s e t1 2 ． V a l u e＆C h r 1 3 ＆ ” 横坐标 Y ”一 ＆ T e mp R a g ． O ff s e t ， T e m p O f f s e t1 3 ． V a l u e＆C h r 1 3 ＆ ” 高 程 H ”一 ＆ T e m p R a g ． O ff s e t ， T e m p O ff s e t 1 4 ． V a l u e ，v b O K O n l y ， ” 测量成果查询” E n d S u b 在编辑完代码后保存 ， 退出“ 设计模式” ， 使用 E x c e l 固有的查询功能在 “ 经纬仪导线测量计算 ” 工 作表中 自动查询 出“ 测线” ， 在 大表 中看着肯定不方 便 ， 这时单击“ 测线成果查询” 按钮 ， 即出现简练的 测量成果 见 图 2 ， 其 中包括测线 测站 一测点 、 边长、 方位角、 纵座标值、 横座标值、 高程值。 图 2 测量成果查询 3 测量公 式输入 E x c e l 有很丰富的函数功能 ， 经纬仪导线计算只 需要用到其 中有关的数学 函数 、 三角函数 、 以及少量 的信息 函数和文本函数 ， 按照导线计算过程的公式 ， 即输 即算。由于 E x c e l 具 有拖动功能 ， 所 以公式不 必重输 。在 E x c e l 中进行经纬仪导线计算 的关键在 于首先要解决六十进位制与十位制互换的问题 ， 否 则在 E x c e l 工作表中进行测量计算连输入与输出都 成问题 。针对这一需求可 以利用 E x c e l V B A建立两 个模块 ， 它们各 自的编程代码如下。 1 模块 1 六十进制化十进制 F u n c t i o n d e g D E s i g n S g n D E D E A b s D E 0 ． 0 0 0 0 0 0 0 0 0 1 d l I n t D E d 2 I n t D E d 1 1 0 0 d 3 DE dl d 2／ 1 0 0 d e g s i g n d l d 2／6 0 d 3／0． 3 6 En d F un c t i o n 2 模块 2 十进制化六十进制 F u n c t i o n d in s D E s i gn S g n DE D E A b s D E 0 ． 0 0 0 0 0 0 0 0 0 1 d l I n t D E d 2 D E d 1 3 6 0 0 d 3 I n t d 2／6 0 d 4 d 2 一 d 3 6 0 d in s s i gn f d l d 3／1 0 0 d 4／l O O O O En d F un c t i o n 经纬仪导线测量计算完成后 ， 即可以打印报表 形式输 出， 作硬拷贝保存。 采矿技术 4 结 论 矿山测量计算非常适合在 E x c e l 中进行 。应用 E x c e l 进行测量计算与平差 的过程实际上就是输人 数据和输人公式的过程。在 E x c e l 中测量计算的过 程与传统的手算表格 近似 ， 过程清楚 ， 界 面直观 ， 检 查修改方便。特别是即输即算， 有较大的灵活性。 上接第 8 8页 矿 、 方铅矿、 黄铁矿 ， 并充填在岩石孔隙中 ， 形成浸染 状矿。由于矿液沿构造裂隙运移 比较方便 、 成矿作 用较充分， 使充填成矿作用形成的浸染状矿体规模 较小。由于 S n 、 C u 、 P b 、 z n “、 F e 离 子地 球 化学性质和迁移能力的差异 ， 发生 了沉积分带。 4深部成矿预测 根据矿床地质特征 、 成矿模式 ， 结合区域成矿背 景分析， 预测在曼家寨矿床倾斜延伸方向仍赋存有 矿体 ， 曼家寨矿床深部具有“ 第二赋矿空 间” 。依据 如下 。 1 曼家寨西部为一单斜褶皱带 ， 地表有南北 向花 岗斑岩展布 ， 南北向构造断裂发育。曼家寨 矿 带横 向上含矿层大致呈 2 0 。～ 4 0 。 的倾斜 角继续 向 西部延伸 ， 其岩性特征为大理岩、 石英云母 片岩 、 矽 卡岩复合岩性组合 ， 与曼家寨矿段含矿岩性组合基 本一致， 对应矿体呈叠瓦状 向深部延 伸。波状起伏 发育的纵向褶皱构造伴随含矿层继续沿倾斜方 向向 深部延伸， 褶皱部位形成一系列的层间剥离空间， 它 是含矿热液运移的通道 ， 也是矿质沉淀富集的场所。 结合横向上金属分带及叠瓦状分布规律， 预测曼家 寨西延深部赋存锌铜工业矿体。 2 曼家寨矿区纵向上由北 向南 1 1 0 0～ 8 0 0 m 标高以上 ， 是主要锡锌矿体赋存部位 ， 矿体厚大、 品 位较富， 是历次勘查工作 的重 点， 储量基本 控制清 楚 ， 称之为第一含矿空 间。在其底部 由北向南 1 0 0 0 7 0 0 I l l_ 标高与 8 0 0 4 0 0 m标高范围之 间通过 以 往稀疏普查孔控制 ， 具有与第一含矿空间相同的含 矿层 ， 岩性组合仍是片岩、 矽卡岩、 大理岩组合 ， 构造 剥离空间、 构造破碎带发育 ， 普查孔均控制到不同厚 度的岩脉 、 岩株 ， 钻孔底部个别孔已揭露凸起的花 岗 岩体。其成矿地质条件与第一含矿空 间相似 ， 称之 为第二含矿空 间。物探 重力异常显示深部 7 0 0 m 标高左右 有 隐伏花 岗岩体隆起 ， 而且带状磁 化规 由于 E x c e l 具有矩阵计算功能， 对于各级各类工程 控制网的严密 平差 ， 使用 E x c e l 来进行更显得举重 若轻 ， 直接应 用严密平差公式 和 E x c e l 的功能就可 以完成平 差 ， 在速 度上不 逊 于应用 软件程 序。而 V B A的辅助应用则更使得 E x c e l 如虎添翼 。 作者简介 褚建民 1 9 5 8一 ， 男， 工程师， 主要从事矿山测量 技 术工作 ， Em a i l s u g u i l 1 1 2 g ma i l ． c o mo 模大 ， 形态完整 ， 呈南北 向带状分布； 化探锡锌铜砷 异常与磁异常吻合 ， 局部还 出现钹 、 铋 、 银 、 钨异常。 结合垂 向上金属分带规律 ， 预测该位置赋存锡锌工 业矿体 。 5 成矿预测效果 通过曼家寨矿区钻探工程验证 ， 在曼家寨深部 第二含矿空间发现了多个隐伏的工业矿体 ， 其 中以 c 1 、 c 2矿体规模最大。c 1号锌铜矿体形态呈条带 状， 赋存于 1 0 5 0 ～ 8 0 0 m标高， 9 5～ 7 1 号剖面线问， 近南北走 向， 走向延伸长大于 3 2 0 m， 倾 向延伸 8 0～ 1 2 0 11 3_ 由于矿权 因素西边未控制尖灭 ， 水平宽度 8 0～1 2 0 m， 矿体收缩膨胀现象明显 ， 厚度变化较大 ， 厚 2～2 3 ． 9 7 m， 锌 品 位 2 ． 3 1 ％ ～6 ． 8 2 ％ ， 平 均 5 ． 1 2 ％ ， 部分地段出现单铜矿或锌铜共生 ， 平均铜 品 位 0 ． 3 2 ％； C 2锡 锌矿体 ， 矿体形 态呈条带状 ， 赋存 于 6 5 0～ 5 2 0 m标高， 2 3～ 7号线间。 近南北走 向， 走 向延伸长大于 3 0 0 m， 倾 向延伸 8 0～1 2 0 m， 水平宽 度 8 0～1 6 0 m， 厚 2 ． 0 5～3 0 ． 4 8 i n ， 锌品位 2 ． 4 5 ％ 一 4 ． 0 3 ％ ， 平均 3 ． 3 8 ％ ， 锡往往以单锡形式存在 ， 锡锌 共生较少。另外， 曼家寨西延深部揭露 了数十个规 模较小的工业矿体。上述矿体产出位置上下盘不远 处多见侵入岩脉、 岩株或深部花岗岩凸起。 6 今 后找矿方 向 上述找矿实践成果证明了曼家寨矿区长期勘探 实践总结的成矿规律对本区的深部找矿预测有重要 的指导作用 。根据成矿规律 ， 预测曼家寨深部 Y坐 标 2 6 0 0 m以西及矿 区南西部位的辣子寨深部仍有 很好的找矿前 景， 且 辣子寨 深部 应存在 多个 赋矿 空 间。 参考文献 [ 1 ] 乐 昌硕 ， 等． 岩石学[ M] ． 北 京 地质 出版社 ， 1 9 8 4 ． 收稿 日期 2 0 0 7 0 81 0