基于三維激光掃描技術(shù)的某金礦采場體積測量與計(jì)算

李鵬宇 姜 岳 宗 琪 牟占武 袁 鑫 姜 巖

(1.山東科技大學(xué)測繪科學(xué)與工程學(xué)院，山東 青島 266590；2.山東科技大學(xué)地球科學(xué)與工程學(xué)院，山東 青島 266590；3.山東黃金礦業(yè)股份有限公司，山東 煙臺 264000)

在黃金礦山開采中，采場驗(yàn)收測量與計(jì)算是一項(xiàng)重要的日常工作。傳統(tǒng)采場驗(yàn)收測量一般采用支距丈量法，即首先利用全站儀或激光測距儀進(jìn)行觀測，然后采用一定的間距截取平行斷面，計(jì)算各斷面的面積，進(jìn)而利用梯形公式計(jì)算總體積[1]。該方法通過在礦房頂板上布設(shè)控制點(diǎn)采用點(diǎn)上對中方式進(jìn)行觀測，測量難度較大[2]，導(dǎo)致每次驗(yàn)收測量耗時較長。此外，礦房頂板不穩(wěn)固也在很大程度上增加了測量工作的安全隱患。三維激光掃描技術(shù)突破了傳統(tǒng)單點(diǎn)測量方法，具有高效率、高精度、非接觸三維測量等優(yōu)勢，在諸多領(lǐng)域得到了廣泛應(yīng)用[3]。山東某黃金礦山為提高采場驗(yàn)收測量精度與工作效率，將三維激光掃描技術(shù)應(yīng)用于采場驗(yàn)收測量，對采場空間體積進(jìn)行測量和計(jì)算。

1 三維激光掃描測量原理

本研究應(yīng)用的三維激光掃描儀型號為Trimble FX，首先由激光極管對目標(biāo)發(fā)射出激光脈沖信號，然后由接收器接收目標(biāo)物表面反射回的脈沖信號，根據(jù)每個激光脈沖從發(fā)出到目標(biāo)物表面再返回到接收器所經(jīng)過的時間來計(jì)算距離，同時測量每個脈沖的角度，經(jīng)過數(shù)據(jù)處理可以得到目標(biāo)物的三維真實(shí)坐標(biāo)[4-5]。

三維激光掃描儀發(fā)射器發(fā)射1個激光脈沖信號，經(jīng)物體表面漫反射后，沿幾乎相同的路徑反向傳回至接收器，據(jù)此可以計(jì)算出目標(biāo)點(diǎn)P與掃描儀O的斜距S，控制編碼器同步測量每個激光脈沖橫向、縱向掃描角度觀測值θ、α。三維激光掃描坐標(biāo)系一般為儀器自定義坐標(biāo)系，X軸在橫向掃描面內(nèi)，Y軸在橫向掃描面內(nèi)與X軸垂直，Z軸與橫向掃描面垂直構(gòu)成了右手坐標(biāo)系(圖1)[6-7]。

圖1 三維激光掃描坐標(biāo)系Fig.1 Coordinate system of 3D laser scanning

三維激光掃描坐標(biāo)系統(tǒng)屬于極坐標(biāo)測量系統(tǒng)，圖1中P點(diǎn)的三維坐標(biāo)計(jì)算公式為

(1)

2 采場體積測量與計(jì)算方法

2.1 點(diǎn)云空間配準(zhǔn)

為解決井下復(fù)雜采場空間單次測量點(diǎn)云數(shù)據(jù)與整體坐標(biāo)的統(tǒng)一問題，需要對點(diǎn)云數(shù)據(jù)進(jìn)行空間配準(zhǔn)。現(xiàn)實(shí)中采場空間體積測量工作是由一系列單次測量工作組成，其點(diǎn)云數(shù)據(jù)需要通過配準(zhǔn)統(tǒng)一至整體坐標(biāo)系中[8]。本研究采用基于標(biāo)靶的配準(zhǔn)方法，以轉(zhuǎn)換成符合要求的初始測站數(shù)據(jù)建立的坐標(biāo)系為基準(zhǔn)，通過同名標(biāo)靶的不同坐標(biāo)采用基于羅德里格矩陣坐標(biāo)轉(zhuǎn)換方法進(jìn)行點(diǎn)云數(shù)據(jù)配準(zhǔn)、平差，得到配準(zhǔn)后的拼接站點(diǎn)云數(shù)據(jù)，其中每隔8站加入黑白標(biāo)靶的配準(zhǔn)、平差數(shù)據(jù)，以此類推，即可得到整個采場的點(diǎn)云拼合數(shù)據(jù)[5]。點(diǎn)云數(shù)據(jù)處理流程如圖2所示。

圖2 點(diǎn)云數(shù)據(jù)處理流程Fig.2 Processing flow of point clouds data

羅德里格矩陣坐標(biāo)轉(zhuǎn)換方法的基本原理是將點(diǎn)云數(shù)據(jù)在空間坐標(biāo)系O0-X0Y0Z0中的坐標(biāo)轉(zhuǎn)換至坐標(biāo)系O0-XiYiZi中，經(jīng)過坐標(biāo)軸旋轉(zhuǎn)和平移便可實(shí)現(xiàn)點(diǎn)云數(shù)據(jù)坐標(biāo)系統(tǒng)一。坐標(biāo)系之間的轉(zhuǎn)換公式為

(2)

式中，xi，yi，zi為經(jīng)過轉(zhuǎn)化后的點(diǎn)云數(shù)據(jù)坐標(biāo)；x0，y0，z0為點(diǎn)云數(shù)據(jù)在三維激光掃描儀坐標(biāo)系中的坐標(biāo)； [Δx,Δy,Δz]T為平移參數(shù)矩陣；λ為平移尺度；R為坐標(biāo)旋轉(zhuǎn)矩陣。

2.2 點(diǎn)云清洗

由于采場空間表面粗糙度大、干濕不均，導(dǎo)致在三維激光掃描過程中不可避免地會獲得大量的噪聲點(diǎn)云。該類噪聲點(diǎn)(包括漂移點(diǎn)、孤立點(diǎn)、冗余點(diǎn)、混雜點(diǎn)等)的存在不僅增加了數(shù)據(jù)量，而且嚴(yán)重影響了點(diǎn)云質(zhì)量和后續(xù)三維模型構(gòu)建，會降低數(shù)據(jù)處理效率，因而有必要進(jìn)行點(diǎn)云濾波，剔除該類冗余數(shù)據(jù)。將經(jīng)過濾波后的點(diǎn)云數(shù)據(jù)單獨(dú)生成點(diǎn)云集合，即可完成每一測站的點(diǎn)云數(shù)據(jù)濾波工作。

2.3 三維建模

通過數(shù)據(jù)拼接和數(shù)據(jù)清洗得到了采場空間表面信息，在此基礎(chǔ)上利用曲面擬合、特征面粘貼方法將采場空間表面細(xì)部數(shù)據(jù)進(jìn)行特征化處理，并選用適宜的建模方法對點(diǎn)云數(shù)據(jù)進(jìn)行三維建模。在儀器自身因素、入射角和測距的雙重幾何因素、采場空間環(huán)境和巖體表面特性等一系列因素影響下，本研究采集的點(diǎn)云數(shù)據(jù)在邊墻和拱頂交接部位質(zhì)量較差，存在點(diǎn)云缺失，在拱頂部位，部分區(qū)域出現(xiàn)空洞，需要通過建模方法進(jìn)行處理。選用基于TIN的建模方法(將抽稀后的點(diǎn)云數(shù)據(jù)生成密集且不透明的三角網(wǎng)來生成模型)對點(diǎn)云數(shù)據(jù)進(jìn)行三維建模。

2.4 體積計(jì)算

三維建模后的數(shù)據(jù)為規(guī)則的格網(wǎng)狀數(shù)據(jù)，利用格網(wǎng)將采場頂部的點(diǎn)云切割成n個小方形，通過投影面投影形成n個長方體。取其中任意1個小長方體，假設(shè)采樣間隔為d，其在數(shù)值上等于小投影面的邊長dx，即

d=dx=dy,

高為hi，那么整個采場的體積即為n個小長方體體積的總和，即：

(3)

如圖3所示，采樣間隔越小，對整個采場的劃分便越接近實(shí)際，體積計(jì)算精度也越高，更能凸顯出三維激光掃描技術(shù)的優(yōu)越性。

圖3 采場體積計(jì)算原理示意Fig.3 Calculation theory schematic of stope volume

3 精度分析

影響三維激光掃描空間點(diǎn)位精度的因素有儀器誤差、目標(biāo)物體反射面引起的誤差、外界環(huán)境條件引起的誤差等[9-11]。在進(jìn)行采場體積測量時，井下溫度、濕度、能見度、風(fēng)力以及氣壓等條件變化較大，特別是采場空間有限，激光束在采場表面的入射角較大，采場巖石表面粗糙、顏色灰暗，導(dǎo)致反射回的激光較少，(一般情況下反射率小于50%)，導(dǎo)致設(shè)備的最大掃描距離顯著減小。為提高掃描精度，有必要控制掃描距離，一般在水平采場，掃描距離不宜超過30%標(biāo)稱掃描距離；在傾斜采場不超過20%標(biāo)稱掃描距離。本研究使用的Trimble FX型三維激光掃描儀的中誤差m0=±2 mm。

本研究采用三維激光掃描儀結(jié)合全站儀對采場進(jìn)行掃描測量，首先利用全站儀進(jìn)行控制測量，獲得黑白靶標(biāo)的三維坐標(biāo)，然后根據(jù)靶標(biāo)掃描結(jié)果進(jìn)行坐標(biāo)轉(zhuǎn)換和匹配。靶標(biāo)的平面與高程控制測量精度以及空間點(diǎn)位中誤差的計(jì)算公式為

(4)

式中，mS為測距誤差；mθ、mα分別為測角誤差；S為水平測量距離；mi、mv分別為儀器高和靶標(biāo)高量取中誤差；ρ為S在XY平面的投影值。

本研究控制測量采用拓普康GTS-332N型全站儀，其標(biāo)稱精度分別為

(5)

mα=±2″,mθ=±5″,S=15 m,θ=45°,mi=mv= ±2 mm,m0=±2 mm。經(jīng)式(4)計(jì)算，Mxyz=±4 mm，則三維激光掃描點(diǎn)P的空間點(diǎn)位中誤差

綜上分析，三維激光掃描測量精度與5″級全站儀的測量精度基本相等，由于Trimble FX型掃描儀一測回耗時約5 min，最大掃描速度為800 000 pts/s，可見采用三維激光掃描技術(shù)能夠代替?zhèn)鹘y(tǒng)采場驗(yàn)收測量方法，可以對采場體積進(jìn)行高精度計(jì)算。

4 工程實(shí)例

本研究以山東某金礦為例，選取該礦某采空區(qū)的部分區(qū)段進(jìn)行試驗(yàn)(圖4)。采場控制點(diǎn)由脈外運(yùn)輸巷開始布設(shè)，由外至內(nèi)連續(xù)設(shè)置測站直至采空區(qū)。

圖4 采場平面Fig.4 Planar graph of the stope

4.1 數(shù)據(jù)采集

本研究試驗(yàn)在采場和脈外運(yùn)輸巷分別設(shè)站Scan01、Scan02進(jìn)行采空區(qū)三維激光掃描(圖5)。在掃描前，需在各測站附近設(shè)置一定數(shù)量并在各測站都可以掃描出的小球或目標(biāo)靶，以確保2次掃描后可以進(jìn)行點(diǎn)云數(shù)據(jù)配準(zhǔn)。為確保掃描數(shù)據(jù)一次性配準(zhǔn)成功，往往會多放置小球和目標(biāo)靶，便于在點(diǎn)云數(shù)據(jù)配準(zhǔn)后可以將采場和脈外運(yùn)輸巷進(jìn)行有效連接[12]。

圖5 三維激光掃描站點(diǎn)分布Fig.5 Distribution of the measurementstations of 3D laser scanning

4.2 數(shù)據(jù)處理

點(diǎn)云數(shù)據(jù)采集完畢后首先進(jìn)行預(yù)處理[13]，而后利用Trimble Realworks 6.5軟件進(jìn)行進(jìn)一步處理。利用目標(biāo)球和目標(biāo)靶進(jìn)行配準(zhǔn)實(shí)現(xiàn)點(diǎn)云數(shù)據(jù)拼接，將原本不在一個坐標(biāo)系內(nèi)的點(diǎn)云數(shù)據(jù)統(tǒng)一于同一坐標(biāo)系內(nèi)。以轉(zhuǎn)換成大地坐標(biāo)的第一站數(shù)據(jù)為基準(zhǔn)，首先采用羅德里格矩陣坐標(biāo)轉(zhuǎn)換方法將下一站的點(diǎn)云數(shù)據(jù)進(jìn)行配準(zhǔn)、平差，得到配準(zhǔn)后的兩站的點(diǎn)云數(shù)據(jù)；然后將每組黑白目標(biāo)靶所控制的多站點(diǎn)云數(shù)據(jù)進(jìn)行拼合，得到采場整體點(diǎn)云數(shù)據(jù)(圖6)。

圖6 采場整體點(diǎn)云數(shù)據(jù)Fig.6 Point clouds data of the whole stope

由于三維激光掃描點(diǎn)云數(shù)據(jù)量非常大，若直接采用原始數(shù)據(jù)進(jìn)行建模，則極大降低了建模效率。因此在構(gòu)建巷道模型之前需對模型進(jìn)行抽稀處理。本研究試驗(yàn)選取空間采樣的點(diǎn)云抽稀方法以及選擇基于TIN的建模方法進(jìn)行采場三維建模[14]，結(jié)果如圖7所示。

圖7 采場三維模型Fig.7 3D model of stope

利用Trimble Realworks 6.5軟件對采場模型體積進(jìn)行了計(jì)算，不同采樣間隔對應(yīng)的采場體積如表1所示。分析表1可知：不同采樣間隔對于采場體積計(jì)算精度的影響較小，因此可以根據(jù)采場空間的實(shí)際需求設(shè)置不同的采樣間隔來計(jì)算采場體積。

表1 不同采樣間隔對應(yīng)的采場體積Table 1 Stope volume calculation results under theconditions of different sampling interval

4.3 SURPAC軟件中三維激光掃描數(shù)據(jù)共享

SURPAC軟件擁有完整的三維立體建模體系，是一款集地質(zhì)、測量、采礦的信息共享平臺，已被廣泛應(yīng)用于金屬礦山生產(chǎn)管理工作中。因此，將三維激光掃描數(shù)據(jù)應(yīng)用于SURPAC軟件中有助于實(shí)現(xiàn)采場空間數(shù)據(jù)共享，便于礦山企業(yè)生產(chǎn)管理。由于SURPAC軟件無法批量處理點(diǎn)云數(shù)據(jù)，故本研究將處理后的點(diǎn)云數(shù)據(jù)進(jìn)一步進(jìn)行了抽稀處理，并以“*.dxf”格式導(dǎo)出(圖8)。

圖8 抽稀處理后的點(diǎn)云數(shù)據(jù)Fig.8 Point clouds after tin processing

將抽稀處理后的點(diǎn)云數(shù)據(jù)導(dǎo)入SURPAC軟件中進(jìn)行采場三維建模(圖9)，而后進(jìn)行了采場體積計(jì)算，結(jié)果見表2。

圖9 SURPAC軟件構(gòu)建的采場三維模型Fig.9 3D model of stope established by SURPAC software表2 不同計(jì)算參數(shù)對應(yīng)的采場空間體積Table 2 Stope volume calculation results under theconditions of different calculation parameters

抽稀間距/m不同等高線間距對應(yīng)的采場體積/m30.51.01.52.00.52073.201900.911608.511252.411.02071.161884.461584.121239.221.52024.131824.431538.451210.072.01943.461791.121468.401179.55

分析表2可知：在其他條件固定的情況下，SURPAC軟件計(jì)算出的采場體積隨著等高線間距和抽稀間距的增大而減??；等高線間距對于采場體積計(jì)算精度的影響較大，等高線間距越大，體積計(jì)算結(jié)果便越小，其與采場真實(shí)體積(2 039.14 m3)的誤差也越大[15]。在不影響SURPAC軟件采場體積計(jì)算效率和精度的前提下，為確保采場體積計(jì)算精度，等高線間距應(yīng)小于0.5 m，抽稀間距不宜大于1.5 m。

5 結(jié) 論

(1)應(yīng)用Trimble FX型三維激光掃描儀可以快速、準(zhǔn)確、高效地獲取礦山采場三維空間數(shù)據(jù)，平均每次掃描測量耗時不超過5 min，點(diǎn)位中誤差為±4.5 mm，與5″級全站儀的測量精度相當(dāng)，可見三維激光掃描儀能夠滿足采場驗(yàn)收測量的精度要求。

(2)基于激光掃描的點(diǎn)云數(shù)據(jù)可以構(gòu)建采場三維空間精細(xì)模型，點(diǎn)云數(shù)據(jù)采樣間隔對于采場體積計(jì)算精度的影響較小，合理采樣間隔為50～200 mm。

(3)SURPAC軟件的不同參數(shù)選取對于采場體積的計(jì)算精度影響較大，為確保采場體積計(jì)算精度，等高線間距應(yīng)小于0.5 m，抽稀間距應(yīng)不大于1.5 m。

[1] 賴東林.采空區(qū)測量新方法探討[J].有色金屬：礦山部分，2001(5)：42-43.

Lai Donglin.Discussion on a new method of goaf measurement[J].Nonferrous Metals:Mine Section,2001(5):42-43.

[2] 徐安泰，李傳晗，于 超，等.距離交會法在礦房驗(yàn)收的應(yīng)用及誤差分析[J].計(jì)算機(jī)光盤軟件與應(yīng)用，2012(11)：102-103.

Xu Antai,Li Chuanhan,Yu Chao,et al.Application of intersection law in mine acceptance and error analysis[J].Computer CD Software and Applications,2012(11):102-103.

[3] 吳 侃，汪云甲，王歲權(quán)，等.礦山開采沉陷監(jiān)測及預(yù)測新技術(shù)[M].北京：中國環(huán)境科學(xué)出版社，2012．

Wu Kan,Wang Yunjia,Wang Suiquan,et al.New Technology of Monitoring and Prediction in Mining Subsidence[M].Beijing:China Environmental Science Press,2012.

[4] 閆騰飛，李元輝，徐 帥，等.采空區(qū)三維激光掃描系統(tǒng)地理坐標(biāo)定位方法[J].金屬礦山，2015(5)：135-139.

Yan Tengfei,Li Yuanhui,Xu Shuai,et al.Research on geographical coordinate positioning method of three-dimensional laser scanning system for mine goaf[J].Metal Mine,2015(5):135-139.

[5] 張卓峰.井下激光掃描測量方法及精度分析[D].青島：山東科技大學(xué)，2013.

Zhang Zhuofeng.Laser Scanning Measurement Method and Precision Analysis in Mining[D].Qingdao:Shandong University of Science and Technology,2013.

[6] 柏雯娟.用三維激光掃描技術(shù)監(jiān)測礦山開采沉陷[J].金屬礦山，2017(1)：132-135.

Bai Wenjuan.Mining subsidence monitoring method based on laser scanning technique[J].Metal Mine,2017(1):132-135.

[7] 馮曉剛，劉明星，撒利偉，等.三維激光掃描在土石方量精計(jì)算中的應(yīng)用[J].地理空間信息，2015(4)：49-51.

Feng Xiaogang,Liu Mingxing,Sa Liwei,et al.Application of 3D laser scanning in earth volume calculation[J].Journal of Geomatics,2015(4):49-51.

[8] 韓 亞，王衛(wèi)星，李 雙，等.基于三維激光掃描技術(shù)的礦山滑坡變形趨勢評價方法[J].金屬礦山，2014(8)：103-107.

Han Ya,Wang Weixing,Li Shuang,et al.An evaluation method of mine landslide and deformation trend based on three-dimensional laser scanning technology[J].Metal Mine,2014(8):103-107.

[9] 冒愛泉，朱益虎，郝思寶，等.地面三維激光掃描儀精度測評方法和誤差改正模型研究[J].測繪通報(bào)，2014(2)：72-75.

Mao Aiquan,Zhu Yihu,Hao Sibao,et al.The research of accuracy evaluation method and error correction model based on 3D laser scanner[J].Bulletin of Surveying and Mapping,2014(2):72-75.

[10] 鄭德華，沈云中，劉 春.三維激光掃描儀及其測量誤差影響因素分析[J].測繪工程，2005(2)：32-34.

Zheng Dehua,Shen Yunzhong,Liu Chun.3D Laser scanner and its effect factor analysis of surveying error[J].Engineering of Surveying and Mapping,2005(2):32-34.

[11] 王俊杰.三維激光掃描數(shù)據(jù)采集誤差來源及精度分析[J].科協(xié)論壇，2011(5)：97-98.

Wang Junjie.Sources and accuracy analysis of 3D laser scanning data acquisition error[J].Technology Association Forum,2011(5):97-98.

[12] 王仲鋒，王 渝，陳 雷.三維激光掃描技術(shù)在城市地下空間測量中的應(yīng)用研究[J].測繪與空間地理信息，2017，40(5)：1-4.

Wang Zhongfeng,Wang Yu,Chen Lei.Application of 3D laser scanning technology measure urban underground space[J].Geomatics and Spatial Information Technology,2017,40(5):1-4.

[13] 劉文龍，趙小平.基于三維激光掃描技術(shù)在滑坡監(jiān)測中的應(yīng)用研究[J].金屬礦山，2009(2)：131-133.

Liu Wenlong,Zhao Xiaoping.Study on the application of 3D Laser Scanning Technology in Landslide Monitoring[J].Metal Mine,2009(2):131-133.

[14] 臧 偉，錢 林，孫寶軍，等.地面三維激光掃描技術(shù)在工程測量中的應(yīng)用[J].北京測繪，2015(3)：130-135.

Zang Wei,Qian Lin,Sun Baojun,et al.Application of ground 3D laser scanning technology in engineering survey[J].Beijing Surveying and Mapping,2015(3):130-135.

[15] 沈 寧.三維激光掃描技術(shù)在隧道變形監(jiān)測中的應(yīng)用研究[D].青島：山東科技大學(xué)，2017.

Shen Ning.Research on Laser Scanner Applied for Deformation Monitoring of Tunnel[D].Qingdao:Shandong University of Science and Technology,2017.