張珂， 楊應迪， 劉學通， 蔣成龍

(安徽理工大學 能源與安全學院， 安徽 淮南 232001)

0 引言

隨著煤礦的深入開采，對煤礦的三維建模及可視化的需求越來越迫切，巷道的三維顯示已經(jīng)成為國內(nèi)外礦山數(shù)字化建設的一個方向。計算機仿真技術的發(fā)展，為建立礦井通風系統(tǒng)的三維模型提供了有效的手段[1]。

礦井通風系統(tǒng)圖是對巷道情況最直觀的反映，很多學者研究了三維仿真技術，開發(fā)出了一系列礦井通風系統(tǒng)來實現(xiàn)巷道的三維顯示。如賈慶仁等[2]利用二維巷道雙邊線計算巷道三線的方法進行建模，開發(fā)了虛擬礦井通風系統(tǒng)。蘇楠[3]利用三維仿真技術結合Solidworks建模，實現(xiàn)了礦井通風網(wǎng)絡的三維仿真功能。陳海林等[4]通過三維坐標計算建立巷道模型，生成了三維通風仿真系統(tǒng)。王建民[5]采用線框模型進行巷道建模，生成了巷道模型。李剛等[6]利用軸投影等原理進行巷道建模。以上模型雖然都實現(xiàn)了巷道的三維可視化，但建模算法復雜，且以巷道兩端端點坐標直接確定巷道的走勢，不能較好地體現(xiàn)巷道的起伏。

針對以上問題，本文以煤礦工程平面圖為基礎數(shù)據(jù)，利用巷道間的拓撲關系進行定位點采集，在中心線加載斷面算法實現(xiàn)的基礎上，構建了三維巷道仿真模型?；谌S巷道模型，將通風設施模型同步到三維系統(tǒng)中，生成了礦井通風系統(tǒng)三維模型，生成的立體圖能清晰地展示礦井巷道的整體走勢，結合數(shù)據(jù)庫信息管理可實現(xiàn)對礦井模型的添加、修改等智能操作。

1 三維巷道模型構建

1.1 構建流程

三維巷道模型構建流程如圖1所示。首先在煤礦工程平面圖上進行導線點數(shù)據(jù)采集和巷道斷面及參數(shù)的選擇；然后由導線點逼近中心線處理，在巷道底板中心線上加載巷道斷面圖元，利用巷道間的拓撲關系確定定位點坐標，生成三維巷道模型；最后進行圖形的連通性判斷。

1.2 基礎數(shù)據(jù)來源

煤礦工程平面圖是煤礦生產(chǎn)實體在圖紙上的真實反映，涵蓋了煤礦的基本信息，包括井下巷道、工作面參數(shù)、工作面狀況等[7-9]。三維巷道模型的構建主要通過巷道斷面形狀和巷道物理位置來體現(xiàn)，通過采集巷道底板中心線點的三維坐標，然后沿中心線軌跡所形成的斷面進行加厚處理即形成礦井巷道模型，而這些數(shù)據(jù)都可以從煤礦工程平面圖上獲取。所以，以煤礦工程平面圖作為基礎數(shù)據(jù)的來源繪制三維通風系統(tǒng)圖是可靠的。

圖1 三維巷道模型構建流程Fig.1 Construction flow of 3D roadway model

1.3 巷道網(wǎng)絡結構模型

1.3.1 定位點約定

為了更加清晰地展示巷道三維模型，巷道數(shù)據(jù)模型中的點需要進行分類處理[10]。在空間中，三維建模的基礎就是以(x，y，z)為坐標的三維定位點，定位點的精度高，才能反映巷道整體的形態(tài)，采集對應的定位點的坐標保存到定位點數(shù)據(jù)表。定位點的分類如圖2所示。

圖2 定位點分類Fig.2 Location point classification

定位點：為了確定巷道模型的參考位置和巷道方位等設置的點，可分為基準點和巷道點。

基準點：為了更換圖紙時圖形不改變方位，以方便采集設置的點，一般以不輕易發(fā)生移動和改變的井口作為參考點，以指北針作為方位點。

巷道點：巷道點分為端點和偽點，端點為巷道的起點和終點，偽點為標高差距較大、斷面改變和通風設施存在的地方，包括變坡點、變形點和設施點。

1.3.2 定位點采集原理

在AutoCAD軟件中，利用煤礦工程平面圖進行交互式采集，識別巷道底板中線點三維坐標，根據(jù)定位點的設置及巷道的拓撲關系采集巷道線上的定位點。巷道元素的拓撲關系如圖3所示。

在圖3中，1、2、3表示巷道線，D1,D2,D3,D5,D6,D7分別為巷道1、2、3的端點，其中D2,D5也是巷道1的偽點，表示巷道的交點。D4,D8,D9為巷道中的偽點，可表示巷道的變坡點和變形點等。圖元根據(jù)巷道形狀選擇拱形或者是梯形，采集巷道實際尺寸與之相匹配，與此同時生成三維立體巷道模型，并能在三維立體圖形上展示井巷參數(shù)信息。

圖3 巷道各元素的拓撲關系Fig.3 Topological relation of each element of roadway

1.4 中心線加載斷面算法的實現(xiàn)

巷道三維建模算法為中心線加載斷面算法，巷道斷面沿著中心線移動，然后進行加厚處理生成三維巷道模型。礦井巷道形態(tài)多樣，主要體現(xiàn)在巷道的物理位置關系上。把巷道抽象成一條線，即為巷道中線[11-12]。由于導線點相連接時會形成一條折線，這樣生成的三維立體圖形會有彎曲現(xiàn)象。為了避免這種現(xiàn)象，利用現(xiàn)有的導線點的坐標、巷道寬度、左邊距、高程校正量等數(shù)據(jù)，采用導線點逼近中心線的方法用巷道底板中線代替導線。依次將巷道中所有的導線點逼近到巷道中心線上，減小由導線點生成巷道帶來的誤差，從而把巷道抽象成一條線進行建模。

在礦井中，各條巷道交錯，形態(tài)各異，以半圓拱形巷道為例研究斷面的加載。巷道斷面的加載是巷道斷面沿中心線進行移動，此過程中巷道可以看作是連貫整體，半圓拱形巷道移動軌跡如圖4所示。加載后的圖形為線條圖，為了形成三維立體巷道模型，還需要對巷道模型進行加厚處理。加載形成的巷道實體模型如圖5所示。基于真實的三維坐標進行巷道模型的構建，能清晰地識別立體巷道間的交錯關系，可通過放大、縮小和旋轉進行模型展示。

圖4 巷道斷面移動軌跡模型Fig.4 Movement trajectory model of roadway section

圖5 巷道實體模型Fig.5 Entity model of roadway

1.5 拓撲關系識別

三維巷道模型的生成過程、各分支之間的拓撲關系是通風系統(tǒng)研究的核心問題。模型生成之后需要對建好的模型進行連通性判斷，在圖形拓撲關系的研究中，主要研究分支與巷道端點之間的關系。將通風系統(tǒng)抽象為空間有向圖，有向圖的各端點就是井巷工程三維空間結構模型中的各空間關鍵點，有向圖各邊的方向就是風流方向。采用深度優(yōu)先搜索算法[13-14]計算通風網(wǎng)絡圖中2個端點之間所有的通路。

深度優(yōu)先搜索算法識別步驟如下：

Step1：從數(shù)據(jù)庫中獲得巷道端點、支路數(shù)據(jù)。

Step2：提取出數(shù)據(jù)中端點和邊的關系，填充鄰接矩陣。

Step3：標記和儲存所有度為0的端點(度表示以該點作為一個端點的邊的個數(shù)，度為0表示該點為通路的最后一個點)。對連通支路進行顏色標記。

Step4：若存在度不為0的端點，則繼續(xù)執(zhí)行深度優(yōu)先搜索算法。執(zhí)行深度優(yōu)先搜索算法后，找出所有的連通分支，并存儲在內(nèi)存中。

巷道中風流的方向應從巷道的始端點A1指向末端點A2，當通風網(wǎng)絡改變或者繪制發(fā)生錯誤時，會導致巷道的方向和實際風流方向不一致。在經(jīng)過通風網(wǎng)絡解算后，系統(tǒng)會自動識別風流方向改變的巷道，然后通知巷道對象改變方向，將巷道始端點變?yōu)锳2，末端點變?yōu)锳1。向風流標注符號對象發(fā)送旋轉180°的信息，然后進行顯示。當進行巷道對象刪除操作時，向風流方向?qū)ο蟀l(fā)送消息要求隨之進行刪除操作。

當通風網(wǎng)絡的通路判定完成后，系統(tǒng)將根據(jù)標注及存儲數(shù)據(jù)的情況自動輸出結果圖，并且彈出對話框，給出判定結果。如果輸出結果是連通圖，則全圖會使用一種顏色對點和相連的邊進行顏色標記。如果不是連通圖，會分為幾個連通子集，不同的連通子集系統(tǒng)會自動標記出不同的顏色，并以對話框的形式提醒。

2 礦井通風系統(tǒng)三維模型構建

礦井通風系統(tǒng)三維模型的構建主要包括三維巷道模型和通風設施模型的構建。在三維巷道模型的基礎上，采集通風機和風門等參數(shù)數(shù)據(jù)，構建通風設施模型，生成完整的礦井通風系統(tǒng)三維模型及礦井資源數(shù)據(jù)庫。構建流程如圖6所示。

(1) 基于三維巷道模型，分析通風網(wǎng)絡的拓撲結構，在煤礦工程平面圖上采集通風機、風門等參數(shù)數(shù)據(jù)，建立與通風系統(tǒng)密切相關的風速、風流、風壓、風阻等參數(shù)數(shù)據(jù)表。

(2) 二維平面圖中布置的通風構筑物可基于巷道建模原理，采用巷道繪制類似的方法，利用四邊形來構造實體模型，通過二三維一體化功能，同步到三維系統(tǒng)中，生成對應的三維模型。然后將礦井通風系統(tǒng)圖抽象為有向網(wǎng)絡圖，為確定定位點之間的通風網(wǎng)絡連通性，在進行通路判斷時遍尋每一個定位點[15]。

圖6 礦井通風系統(tǒng)三維模型構建流程Fig.6 Construction process of 3D model of mine ventilation system

風門模型：建模時將風門的三維空間位置、所處巷道、類型及開關狀態(tài)等數(shù)據(jù)存儲到風門屬性表中。

通風機模型：建模時將通風機的三維空間位置、所在巷道及風量等數(shù)據(jù)存儲到通風機屬性表中。

(3) 礦井資源數(shù)據(jù)庫生成。數(shù)據(jù)庫中包括巷道定位點表、風門屬性表、通風機屬性表和通風機巷道關聯(lián)表等。結合礦井資源數(shù)據(jù)庫，可以對礦井通風系統(tǒng)三維模型進行刪除、重繪等處理，如圖7所示。在對話框中輸入始點和末點的節(jié)點號即可查詢建立的巷道分支，并可對該分支進行相應的修改。

圖7 智能操作Fig.7 Intelligent operation

3 工程應用

平頂山煤業(yè)集團有限責任公司六礦通風方式為分區(qū)與中央并列混合式通風，于1970年投產(chǎn)，通風方法為抽出式通風。共有一水平主井、一水平副井、北一副井、北二副井4個進風井進風，北一風井、北二風井2個回風井回風。其中北一風井擔負戊二上山采區(qū)回風，北二風井擔負丁一采區(qū)、丁二采區(qū)、北二下山采區(qū)回風。

將上述礦井通風系統(tǒng)三維模型的構建方法應用到該礦。首先在AutoCAD中打開該礦工程平面圖，在平面圖上選定參考點和井田邊界，然后進行巷道的端點和偽點的采集，端(偽)點號自動生成并且保存到信息采集數(shù)據(jù)庫中，進入點集合，生成的巷道會進入巷道集合。根據(jù)定位點點號在對應三維圖上生成三維巷道模型，軟件中定位點采集界面如圖8所示，在該礦主要采集了3個水平的運輸大巷、北一風井、北一副井等6個井巷工程及風門、通風機等參數(shù)數(shù)據(jù)，建立了相應的巷道模型、風門模型和通風機模型，六礦立體圖基本信息見表1。

圖8 定位點采集界面Fig.8 Acquisition interface of location points

表1 六礦立體圖基本信息
Table 1 Basic information of stereogram of No.6 Coal Mine

項目信息巷道分支/條428采集定位點/個629巷道類型半圓拱形、矩形、梯形等通風設施風門、擋風墻、通風機等

以六礦明斜井、-100 m大巷局部通風系統(tǒng)網(wǎng)絡圖來展示礦井基本的拓撲結構，如圖9所示。圖中為北一副井到北二風井的回風路線，編號為測點，可進行風速、斷面等參數(shù)的測試，提供立體圖繪制所需的參數(shù)。

圖9 局部通風網(wǎng)絡Fig.9 Local ventilation network

在煤礦工程平面圖的基礎上，利用相關建模原理及數(shù)據(jù)的采集生成三維巷道圖，然后根據(jù)巷道的拓撲關系判斷由進風巷到回風巷是否存在串并聯(lián)巷道等。建立的巷道通風立體圖如圖10所示。局部巷道立體圖如圖10(a)所示，對節(jié)點及拐彎處進行平滑處理，使立體圖呈現(xiàn)得更加清晰?？蛇M行風流方向的標注及巷道顏色等的修改，如圖10(b)所示。最終繪制出六礦通風系統(tǒng)整體圖，如圖10(c)所示，該圖能表達出巷道的走勢及工作面的分布狀況，在風井處添加工業(yè)廣場的建模，使立體圖的呈現(xiàn)更加清晰。同時，繪圖過程所建立的基礎資源數(shù)據(jù)庫如圖11所示，收錄了巷道節(jié)點坐標參數(shù)、巷道斷面及巷道名稱等數(shù)據(jù)信息，數(shù)據(jù)庫可以隨著巷道改變而實時更新。

(a) 局部巷道立體圖 (b) 風流方向標注

(c) 通風系統(tǒng)整體

圖11 資源數(shù)據(jù)庫Fig.11 Resource database

將此建模方法應用到平頂山煤業(yè)集團有限責任公司六礦，改變了以往手工繪制圖紙的方式，實現(xiàn)了二維圖到三維圖的快速轉換，更加真實地展示出了六礦的巷道信息及工作面狀況；巷道資源數(shù)據(jù)庫提供了井巷模型的屬性數(shù)據(jù)，使得數(shù)據(jù)既相互獨立又相互兼容，簡化了信息的查詢過程，數(shù)據(jù)庫的信息可實時更新，大大提高了建模效率。

4 結論

(1) 以煤礦工程平面圖為基礎數(shù)據(jù)，設置定位點進行數(shù)據(jù)采集，利用中心線加載斷面的方法建立礦井通風系統(tǒng)三維模型，生成的巷道立體圖對巷道形狀和走勢的表達更加立體和清晰，具有放大、縮小和旋轉等展示功能。

(2) 實際應用表明：該通風系統(tǒng)三維模型的構建方法大大提高了建模效率，可更加真實地展示煤礦巷道信息及工作面狀況，且巷道資源數(shù)據(jù)庫提供了井巷模型的屬性數(shù)據(jù)，可結合建立的資源數(shù)據(jù)庫對礦井巷道信息進行添加、修改和刪除等操作，為下一步研究通風網(wǎng)絡解算和災變模擬提供了優(yōu)良的平臺。