徐志強(qiáng),王龍輝

（核工業(yè)二〇八大隊(duì)，內(nèi)蒙古 包頭 014010）

我國(guó)砂巖型鈾礦儲(chǔ)量計(jì)算主要使用礦體水平投影圖，現(xiàn)有礦體水平投影圖制作，首先是利用地質(zhì)技術(shù)人員提供的礦段揭穿點(diǎn)深度，結(jié)合測(cè)斜成果，計(jì)算得到揭穿點(diǎn)的坐標(biāo)，把揭穿點(diǎn)坐標(biāo)通過(guò)測(cè)繪軟件CASS 7.0展繪到圖形里，按照礦化信息，把工業(yè)孔、礦化孔、異?？椎姆?hào)手工添加進(jìn)圖形；然后，將平米鈾含量、厚度、品位等參數(shù)按一定格式，通過(guò)計(jì)算機(jī)輸入對(duì)應(yīng)鉆孔。依據(jù)基點(diǎn)基線(xiàn)位置，繪出勘探網(wǎng)，結(jié)合鉆孔的坐標(biāo)值范圍生成圖框；最后圈定儲(chǔ)量范圍，傳統(tǒng)方法是根據(jù)鉆孔坐標(biāo)及礦化信息、勘探線(xiàn)位置，通過(guò)AutoCAD中畫(huà)圓及offset命令畫(huà)輔助線(xiàn)，按照礦體圈定的范圍，通過(guò)輔助線(xiàn)的交點(diǎn)來(lái)確定礦體邊界的拐點(diǎn)，并通過(guò)Pline命令手動(dòng)連接礦體拐點(diǎn)，同時(shí)給每一個(gè)礦體加上參數(shù)標(biāo)示框，完成礦體水平投影圖的繪制。

在實(shí)際工作中，需及時(shí)繪制礦體水平投影圖，實(shí)時(shí)更新估算資源儲(chǔ)量，在制作期間應(yīng)注意礦段劃分原則，塊段外推、塊段圈定及圖件繪制等各種事項(xiàng)［1］。礦體投影圖不僅是為了更好的服務(wù)于工作［2］，也是進(jìn)行下一步勘查、礦床評(píng)價(jià)的重要依據(jù)［3］，為儲(chǔ)量的計(jì)算起到重要的作用［4］。在實(shí)際工作中，生產(chǎn)一線(xiàn)工作者使用CAD進(jìn)行礦體水平投影圖的繪制，電子的圖件有助于提高效率、減少人為誤差、方便后期保存及處理，但通過(guò)CAD制圖依然需要手動(dòng)鼠標(biāo)繪制工作效率仍然不高。礦產(chǎn)勘查、儲(chǔ)量核實(shí)中塊段多、數(shù)據(jù)量大，容易使部分塊段選漏，操作過(guò)程單一，出錯(cuò)機(jī)率大，容易造成眼誤或手誤，質(zhì)量自檢性差，檢查過(guò)程費(fèi)時(shí)費(fèi)力，難以從不同角度分析礦體并繪制礦體水平投影圖。因此，開(kāi)發(fā)一種快速、自動(dòng)生成砂巖型鈾礦礦體水平投影圖的程序，能夠進(jìn)一步提高鈾礦地質(zhì)人員的工作效率，對(duì)鈾礦勘查工作具有重要意義。

1 應(yīng)用開(kāi)發(fā)環(huán)境介紹

1.1 應(yīng)用開(kāi)發(fā)環(huán)境的確定

程序以Microsoft office的Excel為數(shù)據(jù)平臺(tái)，用Excel VBA對(duì)數(shù)據(jù)進(jìn)行處理，以AutoCAD為圖形平臺(tái)，通過(guò)AutoCAD VBA進(jìn)行圖形繪制?，F(xiàn)有鈾礦生產(chǎn)數(shù)據(jù)的電子數(shù)據(jù)以E xcel數(shù)據(jù)格式為主，如坐標(biāo)數(shù)據(jù)、測(cè)斜數(shù)據(jù)、礦體數(shù)據(jù)、鈾礦礦化數(shù)據(jù)等，故在Excel平臺(tái)上用Excel VBA進(jìn)行數(shù)據(jù)處理即是最實(shí)用也是最直接快捷，Auto CAD圖形平臺(tái)圖形繪制方便、使用廣泛、繪制圖形精度高，特別是圖形面積計(jì)算精確，所以礦體水平投影圖的圖形平臺(tái)選用Auto CAD。

1.2 Excel VBA簡(jiǎn)介

Excel軟件是0ffice辦公系列軟件之一，具有強(qiáng)大的數(shù)據(jù)錄入和數(shù)據(jù)處理能力［5］。雖然Excel可以滿(mǎn)足通常大部分?jǐn)?shù)據(jù)的處理需求，但要對(duì)一些數(shù)據(jù)做專(zhuān)業(yè)性較強(qiáng)的處理，就要在Excel中使用VBA。VBA誕生于VB之后，繼承了VB語(yǔ)言的很多對(duì)象、屬性和方法，語(yǔ)法結(jié)構(gòu)與VB基本相同。通過(guò)VBA能夠編寫(xiě)程序?qū)崿F(xiàn)可以數(shù)據(jù)處理的自動(dòng)化，極大提升數(shù)據(jù)處理效率［6］。

Excel應(yīng)用程序的開(kāi)發(fā)需要在Visual Basic編輯器（即VBE）中進(jìn)行。VBE能夠與Excel無(wú)縫結(jié)合，但VBE環(huán)境不能獨(dú)立打開(kāi)，必須在啟動(dòng)Excel或其它宿主程序后才能打開(kāi)運(yùn)行，使用Visual Basic編輯器可以完成VBA過(guò)程創(chuàng)建、VBA用戶(hù)窗體創(chuàng)建、對(duì)象屬性查看或修改以及VBA程序調(diào)試等任務(wù)。

1.3 Auto CAD介紹

AutoCAD是目前工程設(shè)計(jì)中應(yīng)用最廣泛的軟件之一，被廣泛應(yīng)用于機(jī)械、建筑、管線(xiàn)施工、家裝、地質(zhì)等行業(yè)［7］。AutoCAD具有強(qiáng)大的二維繪圖功能，如繪制平面圖、剖面圖、進(jìn)行尺寸標(biāo)注等，并可以進(jìn)行二次開(kāi)發(fā)。

1.3.1 Auto CAD Active X Automation技術(shù)

AutoCAD Active X Automation將AutoCAD對(duì)象顯示到“外部世界”，使得許多編程語(yǔ)言或應(yīng)用程序（例如VB、C++、C#、Java、Word VBA、Excel VBA等）可以訪(fǎng)問(wèn)和編輯CAD對(duì)象。VBA通過(guò)Auto CAD Active X Automation接口將消息發(fā)送到AutoCAD，并通過(guò)接口對(duì)AutoCAD進(jìn)行控制。Auto CAD ActiveX Automation技術(shù)簡(jiǎn)化了程序設(shè)計(jì)過(guò)程，便于對(duì)AutoCAD進(jìn)行二次開(kāi)發(fā)，同時(shí)為其它Windows應(yīng)用程序（如Excel和Word）之間的通信提供了公共平臺(tái)［8］。

1.3.2 Auto CAD對(duì)象模型

Application對(duì)象位于A(yíng)utoCAD對(duì)象模型最頂層，其它層狀結(jié)構(gòu)的對(duì)象是它的子對(duì)象。Preferences是系統(tǒng)參數(shù)對(duì)象，可用于A(yíng)utoCAD的參數(shù)設(shè)置。Documents、Document分別是圖形文件對(duì)象集合、圖形文件對(duì)象，通過(guò)Document對(duì)象可實(shí)現(xiàn)對(duì)AutoCAD文件的新建、打開(kāi)、輸入、輸出等操作。Layer、Layers是圖層對(duì)象、圖層集合對(duì)象，用于控制圖層文件對(duì)象，如添加、刪除圖層，設(shè)置圖層的線(xiàn)型、尺寸和顏色等。

2 基本思路與方法

2.1 基本思路

本方法流程如圖1所示，利用鉆孔孔口坐標(biāo)表、勘探線(xiàn)基線(xiàn)起止點(diǎn)坐標(biāo)、方向、基礎(chǔ)孔距、基礎(chǔ)線(xiàn)距通過(guò)VBA圖框及勘探線(xiàn)計(jì)算程序得到圖框、坐標(biāo)及勘探網(wǎng)數(shù)據(jù)表，并通過(guò)VBA自動(dòng)繪制到圖形空間［9］；利用鉆孔孔口坐標(biāo)表、鉆孔測(cè)斜數(shù)據(jù)表、礦體揭穿深度表計(jì)算揭穿點(diǎn)坐標(biāo)，再結(jié)合揭穿點(diǎn)坐標(biāo)、勘探類(lèi)型、鉆孔礦化信息表，將鉆孔位置及礦化信息展上CAD圖，并由上述數(shù)據(jù)計(jì)算礦體拐點(diǎn)坐標(biāo)；將礦體拐點(diǎn)坐標(biāo)上圖，由程序進(jìn)行初步連接，得到礦體水平投影圖。

圖1 基于VBA技術(shù)自動(dòng)生成鈾礦礦體水平投影圖基本流程圖Fig.1 The basic flow chart of automatically generating horizontal projection map of uranium ore body based on VBA technology

2.2 程序?qū)崿F(xiàn)方法

2.2.1 EXCEL VBA數(shù)據(jù)預(yù)處理模塊實(shí)現(xiàn)

2.2.1.1 礦體揭穿點(diǎn)計(jì)算

由測(cè)斜數(shù)據(jù)、揭穿礦體中心深度及孔口坐標(biāo)

計(jì)算礦體揭穿點(diǎn)平面坐標(biāo)。先從測(cè)斜成果中取得揭穿中點(diǎn)深度的上一測(cè)斜點(diǎn)和下一測(cè)斜點(diǎn)成果數(shù)據(jù)，用距離加權(quán)法計(jì)算出X和Y值，再根據(jù)孔口坐標(biāo)計(jì)算出礦體揭穿中點(diǎn)的坐標(biāo)。距離加權(quán)的插值算法（圖2）為：

圖2 距離加權(quán)示意圖Fig.2 Diagram of distance weighting

式中：X、Y—揭穿點(diǎn)的坐標(biāo)，X1、Y1、X2、Y2分別為揭穿點(diǎn)上下兩個(gè)測(cè)點(diǎn)的坐標(biāo)；D—揭穿點(diǎn)的深度，D1、D2分別為上下兩測(cè)點(diǎn)的深度。

2.2.1.2 礦體水平投影圖的圖框和坐標(biāo)格網(wǎng)及勘探線(xiàn)的計(jì)算

1）圖框計(jì)算

利用Excel公式計(jì)算出所有點(diǎn)X、Y坐標(biāo)的極大和極小值，分別在X和Y坐標(biāo)的極大值、極小值坐標(biāo)值再加、減一個(gè)固定值，得到圖框四個(gè)角點(diǎn)的坐標(biāo)。

2）坐標(biāo)格網(wǎng)的計(jì)算

根據(jù)作圖比例尺，計(jì)算出坐標(biāo)網(wǎng)格的坐標(biāo)間距，對(duì)極值坐標(biāo)按坐標(biāo)間距取整，計(jì)算出每條網(wǎng)格線(xiàn)的坐標(biāo)。

3）勘探網(wǎng)的計(jì)算

以基點(diǎn)坐標(biāo)為中心，根據(jù)勘探線(xiàn)方位角及基線(xiàn)方位角、孔距和線(xiàn)距，計(jì)算出每條勘探線(xiàn)的起點(diǎn)和終點(diǎn)坐標(biāo)。

2.2.1.3 礦體拐點(diǎn)的計(jì)算

1）判斷每一個(gè)工業(yè)孔在勘探線(xiàn)方向、基線(xiàn)方向共四個(gè)方向距離最近鉆孔的礦化類(lèi)型。

以工業(yè)孔為中心計(jì)算每個(gè)鉆孔與該工業(yè)孔的距離及方位角，判斷其與勘探線(xiàn)兩個(gè)方向和基線(xiàn)兩個(gè)方向距離最近的鉆孔孔位，每個(gè)方向的方向角范圍設(shè)為30°，距離按地質(zhì)塊段法要求，根據(jù)礦化類(lèi)型判斷，每個(gè)方向有五種可能：無(wú)鉆孔、無(wú)礦孔、異常孔、礦化孔及工業(yè)孔。

推算算法如下：

兩點(diǎn)的方位角及距離計(jì)算方法如下：

K＝（Y2－Y1）/（X2－X1）

如果X2－X1＜0，方位角α＝arctan（K）＋180如 果X2－X1＞0，且Y2－Y1＞0，α＝arctan（K）＋180

如 果X2－X1＞0，且Y2－Y1＜0，α＝arctan（K）＋360

算法的VBA代碼如下：

dis=（（x1-x2）*（x1-x2）+（y1-y2）*（y1-y2））^（0.5）’兩點(diǎn)間距離

k=（y2-y1）/（x2-x1）’求斜率

Ifx2-x1＜0 Then angle=Application.WorksheetFunction.Degrees（Atn（k））+180

Ifx2-x1＞0 Andy2-y1＞0 Then

angle = Application.Worksheet Function.Degrees（Atn（k））’求方位角

Ifx2-x1＞0 Andy2-y1＜0 Then

angle=Application.Worksheet Function.Degrees（Atn（k））+360

If angle＞K D-anglefw And angle＜K D+anglefw Then

If dis＜disKN And dis＜disfw Then……

2）根據(jù)四個(gè)方向關(guān)系計(jì)算礦體拐點(diǎn)坐標(biāo)

根據(jù)塊段法推算拐點(diǎn)坐標(biāo)。每一個(gè)方向?qū)?yīng)五種不同的可能，它們所對(duì)應(yīng)的距離d值不同。若對(duì)應(yīng)方向無(wú)鉆孔，d為基本工程間距的四分之一；若對(duì)應(yīng)方向最近鉆孔為無(wú)礦孔，d取值為兩孔間距的四分之一或基本工程間距的四分之一；若對(duì)應(yīng)方向最近鉆孔為礦化孔，d取值為兩孔間距的二分之一或基本工程間距的二分之一；對(duì)應(yīng)方向最近鉆孔為工業(yè)孔，取d值為0；若距離最近鉆孔為異常孔，按無(wú)礦孔處理（圖3-圖7，紅、藍(lán)、黃點(diǎn)分別表示工業(yè)孔、礦化孔、異?？祝?/p>

圖3 對(duì)應(yīng)鉆孔d值Fig.3 The d value assignment diagram for counterpart borehole

圖4 對(duì)應(yīng)無(wú)礦孔d值Fig.4 The d value assignment diagram for barren bore hole

圖5 對(duì)應(yīng)礦化孔d值Fig.5 The d value assignment diagram for mineralized hole

圖6 對(duì)應(yīng)異?？譫值Fig.6 The d value assignment diagram for abnormal borehole

圖7 對(duì)應(yīng)工業(yè)孔d值為0Fig.7 The d value assignment diagram for the industrial bore

然后編程由已知d1與d2計(jì)算拐點(diǎn)坐標(biāo)（圖8，紅、藍(lán)點(diǎn)分別表示工業(yè)孔、礦化孔），算法如下：

圖8 由兩個(gè)d值計(jì)算拐點(diǎn)坐標(biāo)Fig.8 Calculating the inflection point coordinates from two d values

式中：X、Y—計(jì)算拐點(diǎn)坐標(biāo)；Y0、X0—工業(yè)孔的坐標(biāo)；ΔY、ΔX為向Y方向和X方向推算距離。

推算出相互關(guān)聯(lián)的工業(yè)孔組，根據(jù)每一個(gè)工業(yè)孔四個(gè)方向的鉆孔類(lèi)型，把相互連接的工業(yè)孔挑選出來(lái)，并且生成一張礦體分類(lèi)表，并且運(yùn)用EXCEL的公式把計(jì)算出來(lái)的相應(yīng)工業(yè)孔的礦體拐點(diǎn)坐標(biāo)復(fù)制到表格相應(yīng)的工業(yè)孔的行里。

2.2.2 Auto CAD VBA繪圖模塊實(shí)現(xiàn)

2.2.2.1 VBA實(shí)現(xiàn)鉆孔和平米鈾含量展圖

1）圖符號(hào)的制作及鉆孔揭穿點(diǎn)的展圖

按照現(xiàn)行鈾礦規(guī)范制作鈾礦化信息的圖塊，將其放在圖形文件目錄下，用礦化信息類(lèi)別命名，如“工業(yè)鈾礦孔.dwg”，在自動(dòng)繪制時(shí)直接作為塊展繪到圖形空間。循環(huán)讀取揭穿點(diǎn)表的坐標(biāo)和孔號(hào)，按照表內(nèi)的礦化信息讀取對(duì)應(yīng)的圖塊，將圖塊和孔號(hào)插入到圖形空間中。

關(guān)鍵代碼如下：

2）平米鈾含量展圖

在鉆孔號(hào)后，需要按照特定格式標(biāo)明鉆孔的所屬塊段號(hào)碼、水平厚度、平均品位、米百分含量，按照已知比例尺算得每個(gè)字的長(zhǎng)、寬進(jìn)，按照規(guī)范要求自動(dòng)繪制。

關(guān)鍵代碼如下：

2.2.2.2 VBA圖框和坐標(biāo)格網(wǎng)及勘探線(xiàn)的自動(dòng)生成

1）圖框和坐標(biāo)格網(wǎng)的自動(dòng)生成

利用EXCEL計(jì)算出的區(qū)域內(nèi)X的最大及最小值和Y的最大最小值，再加上圖框的增加值得到圖框的四個(gè)角點(diǎn)，并對(duì)四個(gè)角坐標(biāo)取整公里網(wǎng)格，繪制出圖框的內(nèi)外框，并根據(jù)Y的極值和X最大值，繪制圖名并放到中間；根據(jù)圖形比例尺繪制公里網(wǎng)格并標(biāo)注。

關(guān)鍵代碼如下：

2）勘探網(wǎng)的自動(dòng)生成

根據(jù)勘探網(wǎng)的基點(diǎn)、基線(xiàn)方向、勘探線(xiàn)方向、線(xiàn)號(hào)距、孔號(hào)距及勘探網(wǎng)號(hào)繪制勘探網(wǎng)。首先，根據(jù)基點(diǎn)和基線(xiàn)方向繪制0號(hào)基線(xiàn)，再根據(jù)基點(diǎn)、勘探線(xiàn)方向和線(xiàn)號(hào)距計(jì)算每條勘探線(xiàn)的起點(diǎn)和終點(diǎn)坐標(biāo)（算法見(jiàn)下）以及勘探線(xiàn)號(hào)，然后繪制勘探線(xiàn)，并在勘探線(xiàn)起點(diǎn)和終點(diǎn)標(biāo)明勘探線(xiàn)號(hào)。

勘探線(xiàn)起點(diǎn)和終點(diǎn)計(jì)算方法（圖9）為：

圖9 勘探線(xiàn)端點(diǎn)計(jì)算示意Fig.9 Calculation of end point of exploration line

式中：X1、Y1—計(jì)算勘探線(xiàn)起點(diǎn)坐標(biāo)；Y0、X0—基點(diǎn)的坐標(biāo)，ΔY、ΔX—向Y方向和X方向推算距離；β—勘探線(xiàn)方位角。

勘探線(xiàn)自動(dòng)生成的主要代碼如下：

2.2.2.3 VBA實(shí)現(xiàn)礦體邊界的初步圈定

礦體邊界由拐點(diǎn)坐標(biāo)連接構(gòu)成，并且需要確定拐點(diǎn)連接方向及順序。使用GIS點(diǎn)在線(xiàn)段左邊判斷原理算法結(jié)合兩點(diǎn)間距離的限差進(jìn)行初步圈定。GIS左手算法，已知兩點(diǎn)A（X1，Y1）和B（X2，Y2），線(xiàn)外一點(diǎn)P（X0，Y0），（X2-X1）（Y0-Y1）-（Y2-Y1）（X0-X1）＞0則點(diǎn)P在線(xiàn)段AB的左側(cè)（圖10）。

圖10 點(diǎn)在線(xiàn)段左邊判斷Fig.10 Judgment diagram of point on the left side of line segment

3 礦體水平投影圖自動(dòng)生成實(shí)例應(yīng)用及與手工圖對(duì)比分析

3.1 實(shí)例應(yīng)用

本研究方法已經(jīng)用于實(shí)際生產(chǎn)中，現(xiàn)以某鈾礦詳查項(xiàng)目為例，對(duì)本程序進(jìn)行實(shí)例驗(yàn)證。

3.1.1 數(shù)據(jù)預(yù)處理模塊的驗(yàn)證

3.1.1.1 揭穿中點(diǎn)計(jì)算

在實(shí)際的野外生產(chǎn)工作中，鉆探施工期間，鉆孔往往不是直孔，均會(huì)在一定程度上朝一定的方位進(jìn)行偏移，所以孔口坐標(biāo)并不能代表鉆孔礦體的所在位置，一般根據(jù)孔口坐標(biāo)及鉆孔的偏斜方位及偏斜距計(jì)算鉆孔礦體的實(shí)際坐標(biāo)。本文由EXCEL表格的礦體頂板和底板數(shù)據(jù)及揭穿中點(diǎn)的上一點(diǎn)和下一點(diǎn)的測(cè)斜坐標(biāo)數(shù)據(jù)計(jì)算出揭穿中點(diǎn)的坐標(biāo)（圖11，數(shù)據(jù)單位：m）。

圖11 揭穿點(diǎn)計(jì)算示意圖Fig.11 Schematic diagram of intersection point calculation

3.1.1.2 圖框和坐標(biāo)格網(wǎng)及勘探線(xiàn)的計(jì)算

在制作礦體水平投影圖時(shí)，必須先計(jì)算好圖框和坐標(biāo)網(wǎng)格及勘探線(xiàn)，才能開(kāi)展下一步工作，可以由EXCEL公式統(tǒng)計(jì)出X坐標(biāo)、Y坐標(biāo)的極大和極小值，并輸入礦床的勘探網(wǎng)的基點(diǎn)坐標(biāo)、勘探線(xiàn)方位角、孔距、線(xiàn)距、勘探網(wǎng)號(hào)等參數(shù)，點(diǎn)擊按鈕直接生成圖框和勘探線(xiàn)（圖12，數(shù)據(jù)單位：m）。

圖12 圖框及坐標(biāo)網(wǎng)格計(jì)算示意圖Fig.12 Schematic diagram of Frame and coordinate grid calculation

3.1.1.3 礦體拐點(diǎn)的計(jì)算

圖框和坐標(biāo)格網(wǎng)及勘探線(xiàn)生成后，便可開(kāi)展下一步工作，對(duì)礦體的拐點(diǎn)進(jìn)行計(jì)算，才能對(duì)礦體的范圍進(jìn)行圈定。由揭穿點(diǎn)坐標(biāo)、勘探線(xiàn)方位角、基線(xiàn)方位角、勘探線(xiàn)方向基本工程間距、基線(xiàn)方向基本工程間距，利用VBA編寫(xiě)的函數(shù)推算出工業(yè)孔四個(gè)方向所對(duì)應(yīng)情況（圖13，數(shù)據(jù)單位：m）。再由推算的情況計(jì)算每一個(gè)工業(yè)孔的礦體拐點(diǎn)（圖14，數(shù)據(jù)單位：m）。

圖13 四個(gè)方向點(diǎn)判斷示意圖Fig.13 Schematic diagram of four directional points

圖14 拐點(diǎn)計(jì)算示意Fig.14 Indication of inflection point calculation

3.1.2 Auto CAD繪圖模塊的驗(yàn)證

EXCEL模塊對(duì)需要繪制圖形數(shù)據(jù)處理好了后，通過(guò)程序的執(zhí)行按扭自動(dòng)生成礦體水平投影圖。

3.1.2.1 鉆孔和平米鈾含量展圖

通過(guò)程序進(jìn)行計(jì)算后，對(duì)數(shù)據(jù)進(jìn)行預(yù)處理，生成每個(gè)鉆孔的礦體揭穿點(diǎn)（圖15，數(shù)據(jù)單位：m），在利用揭穿點(diǎn)數(shù)據(jù)在A(yíng)uto CAD中進(jìn)行揭穿點(diǎn)的展布，形成揭穿點(diǎn)展圖效果圖，可以確定揭穿點(diǎn)的位置及鉆孔的礦化類(lèi)型，與以往模式一致，但速度較快，較為便捷（圖16，數(shù)據(jù)單位：m）。

圖15 揭穿點(diǎn)表意圖Fig.15 List of intersection the point of intent

圖16 揭穿點(diǎn)展圖效果圖Fig.16 Intersection point rendering

3.1.2.2 圖框和坐標(biāo)格網(wǎng)及勘探線(xiàn)的自動(dòng)生成

在前期數(shù)據(jù)比較齊全的情況下，點(diǎn)擊按鈕直接生成圖框和勘探線(xiàn)，圖框和坐標(biāo)網(wǎng)格會(huì)根據(jù)圖面的大小進(jìn)行適當(dāng)?shù)钠ヅ?，不需要再進(jìn)行人工計(jì)算，勘探線(xiàn)也會(huì)根據(jù)鉆孔的位置進(jìn)行生成，如圖17所示，圖框及坐標(biāo)格網(wǎng)基本覆蓋鉆孔范圍，勘探線(xiàn)及基線(xiàn)也進(jìn)行了自動(dòng)繪制，勘探線(xiàn)基本為北西向南東，勘探線(xiàn)數(shù)量為13條，走向?yàn)?35°，基線(xiàn)基本為北東向南西，基線(xiàn)為1條，走向?yàn)?5°，為下一步礦體的圈定做好了準(zhǔn)備工作。

圖17 圖框和坐標(biāo)網(wǎng)格及勘探線(xiàn)效果圖Fig.17 Overview of frame and coordinate grid of exploration line

3.1.2.3 礦體邊界的初步圈定

程序可根據(jù)區(qū)內(nèi)每個(gè)鉆孔的實(shí)際情況，在計(jì)算出礦體拐點(diǎn)坐標(biāo)的基礎(chǔ)上，圈定出礦體的邊界范圍，對(duì)無(wú)法連接的礦體進(jìn)行分割，圈定出相應(yīng)的礦體塊段，如圖18所示，礦體基本走向?yàn)楸睎|南向西向，且被劃分為多個(gè)塊段，初步圈定了礦體的邊界，與人工制作圖件對(duì)比基本無(wú)較大的差異，圈定的邊界準(zhǔn)確可靠。

圖18 礦體初步圈定效果圖Fig.18 Overview of preliminary delineated orebody

3.2 程序自動(dòng)成圖與手工成圖對(duì)比分析

3.2.1 效率對(duì)比

礦體水平投影圖制作，程序繪制效率要遠(yuǎn)比傳統(tǒng)方法繪制的效率高，程序成圖所需時(shí)間基本為15 min左右，而人工用時(shí)基本為665 min，為程序成圖時(shí)間的44倍，這也體現(xiàn)了程序成圖的優(yōu)勢(shì)，所需時(shí)間較短，提高了工作效率，能夠?qū)D件進(jìn)行及時(shí)的更新（表1）。

表1 程序自動(dòng)成圖與手工成圖效率對(duì)比Table 1 Comparison of efficiency between automatic mapping and manual mapping

3.2.2 質(zhì)量對(duì)比

因程序是以Excel輸入的數(shù)據(jù)為基礎(chǔ)的，所以只需要檢查輸入數(shù)據(jù)有沒(méi)有錯(cuò)誤，程序自動(dòng)生成錯(cuò)誤率很低。而手工繪制礦體水平投影圖，繪制過(guò)程避免了人為誤差，錯(cuò)誤率較高。但從礦體水平投影圖的內(nèi)容來(lái)看，礦體的圈定上，手工圈定的相比程序圈定的精準(zhǔn)，程序是初步圈定，可以先給地質(zhì)研究區(qū)圈定大致的礦體范圍，進(jìn)而精確圈定，隨著程序在日后與生產(chǎn)實(shí)踐結(jié)合逐步完善，最終一定能實(shí)現(xiàn)礦體精確圈連。

4 結(jié)論

1）利用Excel方便的數(shù)據(jù)輸入、輸出與編輯功能和強(qiáng)大的公式運(yùn)算能力，結(jié)合AutoCAD完善的圖形平臺(tái)，應(yīng)用VBA技術(shù)編程已實(shí)現(xiàn)礦體水平投影圖自動(dòng)繪制。

2）本應(yīng)用程序所用數(shù)據(jù)平臺(tái)和圖形平臺(tái)都是生產(chǎn)中的常用軟件，具有適用性和易操作性。

3）在使用程序進(jìn)行礦體水平投影圖的生成，用時(shí)較短，效率較高，且繪制的準(zhǔn)確性較人工的高，基本做到避免人為誤差。