張光磊，張登崤，程海星

(中煤西安設(shè)計工程有限責(zé)任公司 陜西 西安 710054)

礦井工作面為煤炭開采的第一現(xiàn)場，位于地下深處，由于礦道狹窄、作業(yè)設(shè)備多且環(huán)境溫度高，極易發(fā)生安全事故，其中在煤炭開采過程中，從煤炭本身和礦井壁散發(fā)的瓦斯在礦壓環(huán)境不佳且溫度過高環(huán)境中，易發(fā)生瓦斯爆炸情況[1-2]。因此，礦井工作面的礦壓處理極其重要，是煤礦企業(yè)管理環(huán)節(jié)的重中之重。為解決礦井工作面礦壓過程中存在的問題，楊帥等[3]通過監(jiān)測和解算礦井內(nèi)礦壓參數(shù)，實現(xiàn)礦井綜采工作面礦壓顯現(xiàn)仿真，該方法雖然可實現(xiàn)礦井綜采工作面礦壓顯現(xiàn)仿真，但受礦井內(nèi)環(huán)境因素影響，無法做到實時解算監(jiān)測到的參數(shù)，存在一定程度延遲，故此應(yīng)用效果不理想；李曼等[4]通過設(shè)置監(jiān)測點對礦井工作面礦壓進行仿真，建立礦井顯現(xiàn)仿真模型，從傳感器獲取到的礦壓分布云圖內(nèi)分析礦井礦壓值，但由于礦井內(nèi)環(huán)境復(fù)雜，監(jiān)測到礦壓據(jù)摻雜不同程度噪聲，而該方法并未對噪聲進行處理，因此應(yīng)用性較差。

虛擬現(xiàn)實技術(shù)將計算機技術(shù)、仿真技術(shù)等結(jié)合在一起通過計算機模擬虛擬環(huán)境，使用戶仿佛置身于真實環(huán)境內(nèi)，提升用戶環(huán)境沉浸感[5]。因此，本文結(jié)合虛擬現(xiàn)實技術(shù)，以建立虛擬礦井工作面環(huán)境方式，研究基于虛擬現(xiàn)實技術(shù)的礦井綜采工作面礦壓顯現(xiàn)仿真方法，為煤炭企業(yè)礦井工作面管理提供技術(shù)支持。

1 虛擬現(xiàn)實技術(shù)的綜采工作面礦壓顯現(xiàn)仿真

1.1 工作面礦壓顯現(xiàn)系統(tǒng)虛擬環(huán)境構(gòu)建

礦井綜采工作面礦壓顯現(xiàn)仿真的目的是通過模擬礦井內(nèi)煤炭開采作業(yè)環(huán)境與操作過程，研究礦井內(nèi)礦壓情況，因此礦井綜采工作面礦壓顯現(xiàn)仿真分為2個步驟：①礦井工作面虛擬場景建立；②利用礦壓變化數(shù)學(xué)模型，分析礦井內(nèi)礦壓變化情況，在礦井工作面礦壓顯現(xiàn)系統(tǒng)虛擬環(huán)境中可實現(xiàn)礦井巷道礦壓計算功能，為后續(xù)礦井綜采工作面礦壓顯現(xiàn)仿真提供數(shù)據(jù)基礎(chǔ)。

(1)虛擬現(xiàn)實技術(shù)開發(fā)工具選擇。在建立礦井工作面礦壓虛擬環(huán)境之初，需選擇相應(yīng)虛擬現(xiàn)實技術(shù)開發(fā)工具，在此選取Multigen軟件負責(zé)構(gòu)建礦井工作面礦壓虛擬環(huán)境，利用Vega軟件設(shè)計交互界面，為Multigen軟件建模提供調(diào)用函數(shù)。三維建模工具選擇Multigen軟件，該軟件為3D建模軟件[6]，其數(shù)據(jù)庫結(jié)構(gòu)為樹狀層次結(jié)構(gòu)，具備光照選擇、紋理映射以及3D聲音等功能。實時仿真軟件選擇Vega，該軟件由Lynx圖形工具箱和C語言運行環(huán)境的Vega函數(shù)庫組成[7]，為Multigen軟件提供用戶界面設(shè)計與場景建模調(diào)用函數(shù)。

圖1 礦井工作面礦壓顯現(xiàn)系統(tǒng)的仿真流程Fig.1 Simulation flow of mine pressure display system in mine working face

1.2 工作面礦壓顯現(xiàn)系統(tǒng)虛擬場景構(gòu)建

礦井工作面虛擬場景構(gòu)建主要由場景建模、幾何建模、場景形象建模、行為建模、人物建模。

(1)場景建模。礦井工作面場礦壓系統(tǒng)景建?？煽醋魇堑V井內(nèi)三維圖形對象的集合，是將礦井工作面礦壓場景以數(shù)學(xué)的方式在計算機內(nèi)呈現(xiàn)出來的過程[8]。礦井工作面礦壓場景建模為礦井巷道模型和設(shè)備模型，巷道模型包括運輸槽、礦壓槽等，設(shè)備模型包括采煤機、傳輸設(shè)備等，而運輸槽、礦壓槽等設(shè)備或場景需經(jīng)過幾何模型、形象模型和運動模型共同組成。

(2)幾何建模。對于場景建模中的礦井巷道模型和設(shè)備模型構(gòu)建，可將礦井巷道和設(shè)備當(dāng)作是由各種形狀的組成的物體，通過繪圖軟件內(nèi)的OpenGL可實現(xiàn)礦井巷道和設(shè)備輪廓構(gòu)建[9]，其步驟為利用VC++軟件依據(jù)礦井巷道和設(shè)備的特征，定義相框和設(shè)備的點、面以及形狀，經(jīng)過連接點、面和形狀，完成礦井巷道和設(shè)備繪制。

(3)場景形象建模。形象建模是將幾何建模后的礦井巷道和設(shè)備進行紋理和顏色填充以及光照渲染，其中紋理映射也稱紋理貼圖，是將場景數(shù)據(jù)庫內(nèi)存儲的礦井巷道和設(shè)備圖像貼到幾何構(gòu)建后的礦井巷道和設(shè)備內(nèi)，通過調(diào)整顏色和光照情況，使礦井巷道和設(shè)備場景與實際場景相似度最大，完成場景形象建模。紋理映射是一個放射轉(zhuǎn)換和雙線性轉(zhuǎn)換的過程[10]，面對礦井工作面復(fù)雜場景，需不同紋理映射方式完成其場景形象建模。

紋理映射方式如下：① Wrap映射方式，當(dāng)需要映射的場景或物體為立方體形狀時，利用該映射方式，從立方體的面的6個頂點出發(fā)，紋理貼圖方向依據(jù)該6個頂點呈現(xiàn)柱面環(huán)繞，完成紋理映射。②Cylindrical映射方式，該映射方式適用于柱形物體貼圖。③Spherical該紋理方式適用于有圓弧狀的物體，如礦井巷道頂端和側(cè)邊緣等。④Planai Eahc Face映射方式，該映射方式適用于大面積平面，如礦井巷道地面、煤炭斷層平面等，均采用該技術(shù)進行映射。

(4)行為建模。行為建模是礦井工作面虛擬場景的重要組成部分，通過行為模型，實現(xiàn)礦井工作面虛擬場景運動控制，井工作面虛擬場景的行為建模利用Multigen軟件內(nèi)的DOF節(jié)點建立運動設(shè)備或人物的相對坐標后，利用Vc模塊建立行為對象類別，完成行為建模過程。

(5)人物建模。用戶通過操控礦井工作面虛擬場景內(nèi)的人物，實現(xiàn)虛擬對象和實體對象之間的數(shù)據(jù)交換，其中虛擬人物利用獨立運動單元和活動關(guān)節(jié)鉸接相結(jié)合，以曲面、三角形等形狀呈現(xiàn)運動單元，關(guān)節(jié)則利用曲面呈現(xiàn)，虛擬人物的手、面部表情和口型等均使用復(fù)雜自由度驅(qū)動構(gòu)成。虛擬人物構(gòu)建完成后，利用手柄按鈕操控人物動作，用戶可通過手柄控制虛擬場景人物進入到礦井巷道內(nèi)，可操控設(shè)備運行，變更設(shè)備位置等。

孕期黃體酮和雌激素的改變 、血流加速導(dǎo)致牙齦腫脹，容易出現(xiàn)牙齦炎、牙齦腫脹，牙齦出血或者牙周病，嚴重時會影響牙齒軟組織周圍的韌帶和骨骼，這些疾病會導(dǎo)致早產(chǎn)和低出生體重兒。分娩結(jié)

1.3 虛擬環(huán)境漫游與交互設(shè)計

礦井工作面礦壓顯現(xiàn)系統(tǒng)漫游方式可劃分為自動式、查詢式和交互式三種形式漫游，用戶可自定義漫游位置與方式，迅速到達礦井工作面的任意位置。而虛擬場景與現(xiàn)實場景的交互則利用數(shù)據(jù)手套實現(xiàn)，利用計算機將數(shù)據(jù)手套內(nèi)的操作指令發(fā)送至礦壓系統(tǒng)的虛擬場景內(nèi)，通過控制虛擬人物實現(xiàn)礦壓系統(tǒng)內(nèi)設(shè)備位置變更以及礦壓操作。

1.4 數(shù)據(jù)庫設(shè)計

礦井工作面礦壓顯現(xiàn)系統(tǒng)虛擬場景構(gòu)建是通過數(shù)據(jù)運算實現(xiàn)的，因此外部數(shù)據(jù)庫與場景數(shù)據(jù)庫構(gòu)建尤為重要，在場景數(shù)據(jù)庫內(nèi)，由壓縮形式的圖像和向量、參數(shù)等形式的圖形組成，是虛擬場景建模必不可少的素材，是提升虛擬場景真實感的途徑之一，而圖像數(shù)據(jù)來源于攝像機拍攝的礦井工作面實景，視頻經(jīng)過相關(guān)軟件剪輯處理后，以O(shè)penGL格式存儲為視頻數(shù)據(jù)。圖形數(shù)據(jù)則由動態(tài)和靜態(tài)兩種圖形數(shù)據(jù)組成，其中動態(tài)圖像數(shù)據(jù)為煤炭開采活動中的相關(guān)參數(shù)，靜態(tài)圖形數(shù)據(jù)則利用高清攝像機拍攝，將礦井工作面的煤層情況、設(shè)備位置、斷層構(gòu)造以及煤炭開采場景的各種細節(jié)以圖片形式記錄下來，對畫面清晰度不夠和光照不均勻數(shù)據(jù)等經(jīng)過處理后存儲為通用圖片格式備用，利用上述數(shù)據(jù)組成外部數(shù)據(jù)庫與場景數(shù)據(jù)庫，通過輸入設(shè)備和輸入接口為虛擬場景構(gòu)建提供場景數(shù)據(jù)和參數(shù)數(shù)據(jù)。

綜上所述，礦井工作面礦壓顯現(xiàn)系統(tǒng)虛擬場景構(gòu)建的基礎(chǔ)即為數(shù)據(jù)庫設(shè)計，而場景數(shù)據(jù)庫構(gòu)建是通過攝像機拍攝即可實現(xiàn)，而外部數(shù)據(jù)庫則是礦井工作面礦壓顯現(xiàn)系統(tǒng)當(dāng)前礦壓參數(shù)，該參數(shù)依據(jù)其礦壓的基本數(shù)學(xué)模型獲取，其步驟如下。

設(shè)b為礦壓數(shù)量，礦井工作面礦壓可以表示為Gy=(gy1，gy2，…，gyb)，由此得到如下方程：

(1)

式中，Rij為i時的第j個礦壓值；Rii為i時礦壓；Pi為第i時刻的礦壓能量代數(shù)和；gyi、gyj為風(fēng)量。

(2)

式中，ci1為第1個元素；Rl、gl分別為礦壓和風(fēng)量。

G=(V，E)為礦井工作面礦壓圖，令|V|、|E|數(shù)值分別為m、n，將該礦壓圖看作一棵樹，依據(jù)余弦方式對該網(wǎng)絡(luò)圖進行排序，其順序為樹在前樹枝在后，令Bk、Uk分別為礦壓圖的基本關(guān)聯(lián)矩陣和獨立回路矩陣，由Bk=(B11B12)、Uk=(U11U12)=(IU12)表示。

其中，B11、B12和U11、U12分別為(m-1)×(n-m+1)、(m-1)×(m+1)階矩陣；I為n-m+1階單位矩陣，依據(jù)礦壓圖的基本關(guān)聯(lián)矩陣，獲取到堵路回路矩如下：

(3)

式中，T為可變參量。QN為礦壓圖分支的有風(fēng)量，則樹和樹枝的風(fēng)量由Qy、Qz表示，依據(jù)風(fēng)量平衡原理可知獨立回路矩陣與礦壓圖與風(fēng)量關(guān)系如下：

QN=QyUk

(4)

通過上述公式可推導(dǎo)出，礦井工作面礦壓圖內(nèi)只有樹的風(fēng)量向量表現(xiàn)為獨立狀態(tài)，因此某一獨立礦壓回路內(nèi)存在分支礦壓阻力代數(shù)與恒等于礦壓能量和，其表達公式如下：

(5)

式中，hj、pj分別為j個分支礦壓阻力和礦壓能量代數(shù)和。

將公式(5)轉(zhuǎn)換為矩陣，其表達式如下：

Uk(HT-PT)=0

(6)

式中，HT、PT分別為礦壓阻力和礦壓能量代數(shù)和的矩陣表達形式。

式(3)、式(4)組成礦井工作面礦壓顯現(xiàn)系統(tǒng)基本回路方程，式(5)和式(6)組成礦井工作面礦壓顯現(xiàn)系統(tǒng)基本數(shù)學(xué)方程，經(jīng)過2個方程組的反復(fù)運算，可獲取到該框架工作面礦壓系統(tǒng)的基本參數(shù)，并利用該參數(shù)構(gòu)建外部數(shù)據(jù)庫，為礦井工作面礦壓顯現(xiàn)系統(tǒng)虛擬場景構(gòu)建提供礦壓相關(guān)參數(shù)。

至此，礦井工作面礦壓顯現(xiàn)系統(tǒng)虛擬場景構(gòu)建完畢，當(dāng)?shù)V井工作面礦壓顯現(xiàn)系統(tǒng)發(fā)生變化時，通過求解礦壓基本數(shù)學(xué)模型，變更礦井工作面礦壓顯現(xiàn)系統(tǒng)虛擬場景，實現(xiàn)礦井綜采工作面礦壓顯現(xiàn)仿真。

2 仿真模擬測試

以某省大型礦區(qū)為實驗對象，使用本文方法構(gòu)建該礦區(qū)礦井工作面礦壓顯現(xiàn)系統(tǒng)虛擬仿真場景，分析該礦井內(nèi)礦壓情況。

2.1 虛擬場景性能測試

虛擬場景構(gòu)建效果是影響礦井工作面仿真結(jié)果的重要因素，從虛擬場景構(gòu)建角度進行測試，分析本文方法虛擬場景構(gòu)建能力，結(jié)果如圖2所示。分析圖2可知，在采煤機控制界面內(nèi)包含任務(wù)列表、機械認知、工具背包、井下導(dǎo)航以及采煤機開關(guān)等控制按鈕，通過該控制按鈕實現(xiàn)采煤機運作，而在圖2(b)內(nèi)，虛擬任務(wù)關(guān)節(jié)彎曲度以及人物動作較為逼真，而人物的行進方向與箭頭指示方向一致，綜上而言，本文方法構(gòu)建的礦井工作面礦壓顯現(xiàn)系統(tǒng)虛擬場景具備較強功能，且畫面感十足，可較好地呈現(xiàn)礦井工作面礦壓顯現(xiàn)系統(tǒng)當(dāng)前情況。

圖2 虛擬場景構(gòu)建測試結(jié)果Fig.2 Test results of virtual scenario construction

為更好地呈現(xiàn)本文方法構(gòu)建的礦井礦壓系統(tǒng)虛擬場景應(yīng)用效果，以礦井巷道為例，設(shè)置當(dāng)前礦井巷道礦壓為2 154 Pa，使用本文方法獲取礦井巷道當(dāng)前礦壓結(jié)果如圖3所示。分析圖3可知，在本文構(gòu)建的工作面虛擬場景內(nèi)，不僅可獲取到該礦井巷道礦壓2情況，還可實時獲取風(fēng)機當(dāng)前功率以及其工作效率。利用本文構(gòu)建的礦壓顯現(xiàn)系統(tǒng)虛擬場景可獲取工作面當(dāng)前礦壓情況，且數(shù)值不存在誤差，可有效應(yīng)用于后續(xù)該礦井礦壓仿真。

圖3 礦井巷道虛擬場景Fig.3 Virtual scene of mine roadway

2.2 礦井工作面礦壓測試

從礦壓監(jiān)測角度展開測試，以該礦井工作面的12個分支礦壓巷道作為測試對象，使用本文方法統(tǒng)計該12個分支礦壓巷道24 h內(nèi)礦壓情況，結(jié)果見表1。

表1 礦壓監(jiān)測結(jié)果Tab.1 Monitoring results of ore pressure Pa

分析表1可知，12個礦井分支巷道的最大礦壓、最小礦壓、平均礦壓以及當(dāng)前礦壓數(shù)值均不同，狀態(tài)描述的是當(dāng)前礦壓情況，在12個礦井分支巷道內(nèi)的當(dāng)前礦壓中，第8個、第9個和第12個礦井分支巷道的當(dāng)前礦壓數(shù)值較其24 h內(nèi)平均礦壓數(shù)值相差較大，被描述為當(dāng)前礦壓過小，而第6個和第10個礦井分支巷道的當(dāng)前礦壓已超過其24 h內(nèi)的平均礦壓，因此其狀態(tài)描述為當(dāng)前礦壓過大，該結(jié)果表明，本文方法可有效仿真礦井工作面內(nèi)礦壓巷道的當(dāng)前礦壓情況，并可依據(jù)24 h內(nèi)該巷道礦壓平均值判斷當(dāng)前礦壓是否符合礦壓需求，具備較高的實際應(yīng)用價值。

2.3 應(yīng)用性測試

對礦井工作面礦壓進行仿真的目的是實現(xiàn)礦井瓦斯?jié)舛缺O(jiān)測，避免礦難事故發(fā)生，對該礦井工作面某一主巷道未進行礦壓處理之前的瓦斯?jié)舛日归_仿真，而后對該主巷道進行5 min礦壓處理，并對礦壓處理后的主巷道瓦斯?jié)舛冗M行仿真，分析本文方法的實際應(yīng)用效果，結(jié)果如圖4所示。為更便捷呈現(xiàn)瓦斯?jié)舛扰c礦井埋深關(guān)系，分析圖4可知，主巷道瓦斯?jié)舛瘸尸F(xiàn)山峰狀分布，在礦井埋深為20～80 m時，其瓦斯?jié)舛茸罡撸M行礦壓前該礦井主巷道的瓦斯?jié)舛雀哌_0.27%左右，隨著礦井埋深的增加，瓦斯?jié)舛戎饾u下降，而進行5 min礦壓處理后的主巷道瓦斯?jié)舛惹€雖然與礦壓處理前整體趨勢相同，但瓦斯?jié)舛认陆得黠@，瓦斯?jié)舛茸罡邽?.22%左右，礦壓5 min后瓦斯?jié)舛认陆导s0.05%，該結(jié)果表明，本文方法可有效對礦井巷道瓦斯?jié)舛冗M行監(jiān)測仿真，也從側(cè)面印證了該方法對礦井工作面礦壓情況仿真效果好。

圖4 礦壓前后瓦斯?jié)舛惹闆rFig.4 Gas concentration before and after ore pressure

3 結(jié)語

本文依托虛擬現(xiàn)實技術(shù)，對基于虛擬現(xiàn)實技術(shù)的礦井綜采工作面礦壓顯現(xiàn)仿真方法展開研究，并從不同角度對該方法進行了測試。測試結(jié)果表明，本文方法構(gòu)建的礦井工作面礦壓顯現(xiàn)系統(tǒng)虛擬場景畫面感十足且功能眾多，可有效仿真礦井工作面礦壓情況；可依據(jù)24 h內(nèi)該巷道礦壓平均值判斷當(dāng)前礦壓是否符合礦壓需求，實際應(yīng)用性強；可有效仿真風(fēng)機礦壓，并描述礦井工作面礦壓調(diào)整過程。