高銀貴，周大偉，安士凱，王 玲，張德民，詹紹奇

(1.鄂爾多斯市華興能源有限責(zé)任公司，內(nèi)蒙古 鄂爾多斯 017000；2.中國(guó)礦業(yè)大學(xué) 環(huán)境與測(cè)繪學(xué)院，江蘇 徐州 221116；3.平安煤炭開(kāi)采工程技術(shù)研究院有限責(zé)任公司，安徽 淮南 232033；4.新能礦業(yè)有限公司，內(nèi)蒙古 鄂爾多斯 017000)

0 引 言

“富煤、貧油、少氣”是我國(guó)的能源現(xiàn)狀，這就決定了煤炭資源在較長(zhǎng)時(shí)間內(nèi)始終是我國(guó)的主要消費(fèi)能源和能源戰(zhàn)略安全的保證[1]。我國(guó)煤炭資源在地理分布上呈現(xiàn)“西多東少”趨勢(shì)，隨著國(guó)家供給側(cè)結(jié)構(gòu)性改革、能源戰(zhàn)略調(diào)整及東部煤礦逐漸關(guān)閉，西部地區(qū)必將成為我國(guó)未來(lái)主要煤炭資源產(chǎn)地[2-3]。煤炭開(kāi)采在滿足國(guó)家能源需求的同時(shí)也導(dǎo)致嚴(yán)重的礦區(qū)地質(zhì)環(huán)境損害問(wèn)題，煤炭開(kāi)采引起巖層及地表變形是產(chǎn)生礦山地質(zhì)環(huán)境損害的根源。西部礦區(qū)地質(zhì)構(gòu)造簡(jiǎn)單，煤層賦存穩(wěn)定，適宜大規(guī)模機(jī)械化開(kāi)采，且由于煤層淺埋，大采高、高強(qiáng)度快速開(kāi)采及地質(zhì)環(huán)境脆弱，西部煤炭開(kāi)采引發(fā)的地表沉陷大多呈現(xiàn)出變形速度快、損害程度深、波及范圍廣的特點(diǎn)[4]。由煤炭開(kāi)采引起的礦區(qū)地表沉陷會(huì)造成地表建構(gòu)筑物以及生態(tài)環(huán)境的破壞，嚴(yán)重時(shí)甚至?xí)l(fā)滑坡、泥石流等災(zāi)害，對(duì)煤礦安全生產(chǎn)以及人們的正常生活產(chǎn)生嚴(yán)重的損害和威脅。礦區(qū)開(kāi)采沉陷監(jiān)測(cè)對(duì)合理指導(dǎo)煤炭資源開(kāi)采活動(dòng)，保護(hù)生態(tài)環(huán)境不受破壞以及人民生命財(cái)產(chǎn)安全具有重大意義。然而，西部礦區(qū)地表沉陷監(jiān)測(cè)數(shù)據(jù)較少，沉陷規(guī)律研究不深入，更沒(méi)有適合的沉陷參數(shù)。尤其國(guó)家對(duì)環(huán)境保護(hù)要求越來(lái)越嚴(yán)格，礦區(qū)環(huán)境損害監(jiān)測(cè)及評(píng)估已經(jīng)成為西部礦山企業(yè)面臨的巨大難題。

現(xiàn)有沉陷監(jiān)測(cè)技術(shù)是以全站儀、水準(zhǔn)測(cè)量、GPS等大地測(cè)量技術(shù)為主的常規(guī)地表沉陷監(jiān)測(cè)技術(shù)。該監(jiān)測(cè)方法是通過(guò)在采空區(qū)上方地表布設(shè)走向和傾向“十字”觀測(cè)線進(jìn)行沉陷監(jiān)測(cè)[5-6]，根據(jù)各觀測(cè)站獲取的平面坐標(biāo)和高程信息。該技術(shù)監(jiān)測(cè)精度高，是目前獲取礦區(qū)地表形變最常用、最有效的手段，在礦山開(kāi)采沉陷研究及實(shí)際工程中發(fā)揮了重要的作用。然而，該監(jiān)測(cè)手段有如下特點(diǎn)：①監(jiān)測(cè)范圍和尺度小，一般以單工作面為監(jiān)測(cè)對(duì)象；②僅以地表移動(dòng)變形為監(jiān)測(cè)內(nèi)容，無(wú)法獲得環(huán)境損害信息(比如植被、水體等)；③成本高、周期長(zhǎng)、工作量大、需要埋設(shè)測(cè)點(diǎn)且難以長(zhǎng)期保存等缺點(diǎn)；④只“點(diǎn)狀”觀測(cè)，數(shù)據(jù)及信息量少，不能完整反映沉陷盆地特征。這些不足導(dǎo)致其已無(wú)法適應(yīng)新形勢(shì)下西部高強(qiáng)度開(kāi)采地表?yè)p害監(jiān)測(cè)任務(wù)，因此如何快速、準(zhǔn)確、全面地監(jiān)測(cè)西部高強(qiáng)度大規(guī)模開(kāi)采引起的地表沉陷與環(huán)境損害是解決問(wèn)題的關(guān)鍵。

隨著近年低成本、輕型無(wú)人機(jī)與技術(shù)的發(fā)展，以無(wú)人機(jī)為搭載平臺(tái)的無(wú)人機(jī)低空攝影測(cè)量技術(shù)(Unmanned Aerial Vehicle,UAV)成為當(dāng)前應(yīng)用較廣泛的測(cè)繪新技術(shù)。UAV攝影測(cè)量是以無(wú)人機(jī)為平臺(tái)，集成高分辨率(多/高光譜)低空光學(xué)傳感器、GPS和IMU等技術(shù)，可同時(shí)獲取地表沉陷變形和地物環(huán)境信息，適應(yīng)于多尺度的監(jiān)測(cè)任務(wù)，具有機(jī)動(dòng)靈活、高效精準(zhǔn)、作業(yè)成本低(突出的時(shí)效性和性價(jià)比)等優(yōu)勢(shì)[7-8]，極大推動(dòng)了攝影測(cè)量技術(shù)在礦區(qū)的應(yīng)用[7,9-11]?；诖耍P者以內(nèi)蒙古鄂爾多斯某煤礦為研究對(duì)象，提出了UAV攝影測(cè)量監(jiān)測(cè)礦區(qū)地表沉陷的思路、關(guān)鍵技術(shù)流程和方法、采用UAV攝影測(cè)量技術(shù)對(duì)該礦區(qū)地表沉陷進(jìn)行監(jiān)測(cè)，采用常規(guī)全站儀/水準(zhǔn)測(cè)量數(shù)據(jù)進(jìn)行對(duì)比研究，分析了UAV攝影測(cè)量監(jiān)測(cè)精度，討論分析該技術(shù)在西部礦區(qū)高強(qiáng)度開(kāi)采地表沉陷監(jiān)測(cè)的可行性。

1 無(wú)人機(jī)監(jiān)測(cè)思路及方法

1.1 基本思路

UAV-攝影測(cè)量可以定期快速獲取監(jiān)測(cè)區(qū)域地表光學(xué)影像數(shù)據(jù)，該影像含有地表三維坐標(biāo)和豐富的地物光譜信息數(shù)據(jù)，通過(guò)對(duì)影像數(shù)據(jù)的處理可以得到監(jiān)測(cè)區(qū)域的地表數(shù)字高程模型(Digital Elevation Model,DEM)和生態(tài)環(huán)境信息(如植被、土壤及水體)，前者通過(guò)多期DEM相減可以進(jìn)一步處理得到地表沉陷變形，后者通過(guò)多期地物對(duì)比分析，可用于研究礦區(qū)生態(tài)環(huán)境精細(xì)演變。

圖1為UAV-攝影測(cè)量監(jiān)測(cè)地表沉陷的示意，當(dāng)工作面推進(jìn)到位置1時(shí)，利用UAV觀測(cè)一次地表，解算得到當(dāng)時(shí)的數(shù)字高程模型DEM1。當(dāng)工作面推進(jìn)到位置2時(shí)，再用對(duì)同一位置地表進(jìn)行第二次觀測(cè)，獲得該時(shí)刻的地表數(shù)字高程模型DEM2。通過(guò)兩期DEM相減，即DEM1減去DEM2，可以得到監(jiān)測(cè)區(qū)域的地表沉陷盆地，如若2次監(jiān)測(cè)是采動(dòng)過(guò)程中進(jìn)行的，獲得的即為動(dòng)態(tài)沉陷盆地；如若為工作面開(kāi)采首末2次監(jiān)測(cè)，獲得的即為穩(wěn)態(tài)沉陷盆地，再根據(jù)工作面開(kāi)采的相關(guān)信息，便可以反演出地表移動(dòng)變形預(yù)測(cè)參數(shù)等。

圖1 無(wú)人機(jī)攝影測(cè)量技術(shù)監(jiān)測(cè)地表沉陷思路Fig.1 Monitoring surface subsidence by UAV photogrammetry technology

1.2 DEM構(gòu)建方法及流程

UAV-攝影測(cè)量技術(shù)通過(guò)對(duì)航測(cè)數(shù)據(jù)進(jìn)行處理得到同一監(jiān)測(cè)區(qū)域不同監(jiān)測(cè)時(shí)期的DEM，多期DEM相減獲得監(jiān)測(cè)區(qū)域的地表沉陷。DEM的構(gòu)建是生成地表沉陷盆地的關(guān)鍵步驟，DEM構(gòu)建流程如圖2所示。首先，利用無(wú)人機(jī)控制平臺(tái)規(guī)劃航飛路線并獲取監(jiān)測(cè)區(qū)域的影像數(shù)據(jù)；再進(jìn)行相機(jī)標(biāo)定獲得相機(jī)畸變參數(shù)，并利用該參數(shù)對(duì)原始影像進(jìn)行畸變校正；然后，對(duì)校正后的影像進(jìn)行相對(duì)定向，使所有像片統(tǒng)一至選定的像空間輔助坐標(biāo)系中；隨后，再利用外業(yè)測(cè)量獲取的像控點(diǎn)數(shù)據(jù)進(jìn)行絕對(duì)定向，使像片由像空間輔助坐標(biāo)系轉(zhuǎn)換至地面坐標(biāo)系；最后，通過(guò)光束法區(qū)域網(wǎng)平差以及影像密集匹配加密點(diǎn)，獲取整個(gè)監(jiān)測(cè)區(qū)域的點(diǎn)云數(shù)據(jù)，再通過(guò)空間內(nèi)插構(gòu)建DEM[12-13]。其中，構(gòu)建DEM的關(guān)鍵技術(shù)如下：

圖2 DEM構(gòu)建流程Fig.2 Flow chart of DEM construction

1)影像畸變校正：無(wú)人機(jī)低空攝影測(cè)量攜帶的相機(jī)大多為非量測(cè)型相機(jī)，獲取的影像存在畸變差，因此在進(jìn)行影像數(shù)據(jù)處理前必須消除畸變差。通過(guò)調(diào)整像主點(diǎn)位置和糾正畸變參數(shù)來(lái)實(shí)現(xiàn)影像畸變校正。

2)相對(duì)定向：在2張具有一定重疊度的相鄰影像上提取相似的特征點(diǎn)，然后利用基于灰度或者特征的匹配算法進(jìn)行同名點(diǎn)匹配，相對(duì)定向后監(jiān)測(cè)區(qū)域所有影像的坐標(biāo)系統(tǒng)一至以像空間輔助坐標(biāo)系為基準(zhǔn)的自由坐標(biāo)系中。

3)絕對(duì)定向：引入像控點(diǎn)數(shù)據(jù)，解算每張像片在控制點(diǎn)坐標(biāo)系下的外方位元素和待定點(diǎn)的地面坐標(biāo)，完成影像由自由坐標(biāo)系統(tǒng)一至地面坐標(biāo)系。

4)光束法區(qū)域網(wǎng)平差：以每張像片組成的光束為平差基本單元，按照共線條件方程列出誤差方程，對(duì)測(cè)區(qū)整體進(jìn)行平差處理獲取每張像片的外方位元素以及待求點(diǎn)的地面坐標(biāo)。

5)影像匹配：尋找2幅或多幅影像之間的同名點(diǎn)，主要方法為基于點(diǎn)特征的匹配，包括SIFT、SURF、Harris-Laplace算子等算法。

2 礦區(qū)工程案例

2.1 研究區(qū)域及數(shù)據(jù)采集

試驗(yàn)區(qū)域?yàn)槲挥趦?nèi)蒙古自治區(qū)鄂爾多斯市的某煤礦2S201工作面。該工作面傾向長(zhǎng)度260 m，走向長(zhǎng)度1 253 m，煤層厚2.9～3.7 m，實(shí)際平均采厚3.17 m，煤層總體傾向西南，平均傾角2°，煤層平均采深約200 m。

采用天寶UX5無(wú)人機(jī)攝影測(cè)量系統(tǒng)進(jìn)行數(shù)據(jù)采集，無(wú)人機(jī)相關(guān)參數(shù)見(jiàn)表1；外業(yè)像控點(diǎn)數(shù)據(jù)由GPS RTK測(cè)量完成，像控點(diǎn)測(cè)量坐標(biāo)系為北京54坐標(biāo)系，數(shù)據(jù)采集現(xiàn)場(chǎng)如圖3所示。

圖3 現(xiàn)場(chǎng)采集數(shù)據(jù)Fig.3 Field data acquisition

表1 天寶UX5無(wú)人機(jī)航測(cè)系統(tǒng)參數(shù)Table 1 Parameters of aerial survey system of Tianbao UX5 UAV

兩期數(shù)據(jù)分別采集于2018年6月9日、2018年9月4日，第一期和第二期無(wú)人機(jī)影像均為560張，相對(duì)地面平均航高為250 m左右，設(shè)計(jì)為東西航向，航帶為28條，考慮到本試驗(yàn)所用無(wú)人機(jī)平臺(tái)較輕、穩(wěn)定性差和抗風(fēng)能力弱的特點(diǎn)，設(shè)計(jì)的航向和旁向重疊率均為80%，見(jiàn)表2。

表2 無(wú)人機(jī)影像數(shù)據(jù)信息Table 2 Information of UAV images

2.2 數(shù)據(jù)處理及精度分析

2.2.1 數(shù)據(jù)處理及結(jié)果

根據(jù)前述無(wú)人機(jī)攝影測(cè)量監(jiān)測(cè)地表沉陷的基本原理，首先構(gòu)建DEM，分別將外業(yè)無(wú)人機(jī)平臺(tái)采集的2期影像數(shù)據(jù)以及GPS RTK采集的像控點(diǎn)數(shù)據(jù)導(dǎo)入PhotoScan軟件中，進(jìn)行初始預(yù)處理剔除掉部分扭曲或模糊的照片，然后對(duì)齊照片；然后進(jìn)行影像匹配獲取影像的同名點(diǎn)坐標(biāo)，并剔除誤匹配點(diǎn)；利用少量的外業(yè)像控點(diǎn)數(shù)據(jù)在像片上刺同名點(diǎn)，進(jìn)行輔助GPS光束法區(qū)域網(wǎng)平差以獲取高精度的姿態(tài)數(shù)據(jù)；最后進(jìn)行密集匹配、SfM攝影測(cè)量等處理，建立實(shí)驗(yàn)區(qū)域的密集點(diǎn)云數(shù)據(jù)，再經(jīng)過(guò)規(guī)則格網(wǎng)內(nèi)插獲取數(shù)字正射影像(DOM)和數(shù)字高程模型(DEM)，如圖4和圖5所示。

圖4 第一期DOM圖像(左)和DEM圖像(右)Fig.4 DOM (left)/ DEM (right)of phase I

圖5 第二期DOM圖像(左)和DEM圖像(右)Fig.5 DOM (left)/ DEM (right)of phase Ⅱ

然后對(duì)利用ArcGIS軟件對(duì)2期DEM進(jìn)行相減，得到無(wú)人機(jī)攝影測(cè)量監(jiān)測(cè)的某煤礦2S201工作面從2018-06-09至2018-09-04間的地表沉陷，最大下沉值為2 487 mm，如圖6所示。

2.2.2 DEM及沉陷盆地精度分析

為分析無(wú)人機(jī)攝影測(cè)量技術(shù)監(jiān)測(cè)礦區(qū)地表沉陷的精度，采用全站儀和水準(zhǔn)儀于2018年6月和9月分別對(duì)布設(shè)的15個(gè)地面測(cè)站點(diǎn)進(jìn)行測(cè)量(地面測(cè)站點(diǎn)分布情況如圖6所示)，得到兩個(gè)時(shí)期的各地面觀測(cè)站的平面坐標(biāo)和水準(zhǔn)數(shù)據(jù)，并與無(wú)人機(jī)數(shù)據(jù)生成的DEM和沉陷盆地上對(duì)應(yīng)點(diǎn)的高程和沉降值數(shù)據(jù)分別進(jìn)行比較。

圖6 無(wú)人機(jī)攝影測(cè)量監(jiān)測(cè)的沉陷盆地Fig.6 Subsidence basin monitored by UAV photogrammetry

1)DEM精度評(píng)價(jià)。根據(jù)全站儀測(cè)量的地面測(cè)站點(diǎn)的平面坐標(biāo)，利用利用ArcGIS軟件提取DEM上對(duì)應(yīng)坐標(biāo)點(diǎn)的高程。兩期DEM高程與地面觀測(cè)站水準(zhǔn)高程對(duì)應(yīng)比較，見(jiàn)表3和表4。

表3 第一期DEM的誤差統(tǒng)計(jì)Table 3 Error statistics of first period DEM

表4 第二期DEM的誤差統(tǒng)計(jì)Table 4 Error statistics of second period DEM

通過(guò)對(duì)地面觀測(cè)點(diǎn)和DEM對(duì)應(yīng)點(diǎn)的高程差值的均方根誤差進(jìn)行計(jì)算，以此評(píng)估DEM的精度，均方根誤差計(jì)算公式為

(1)

式中：M為均方根誤差(即中誤差)；Δi為高程差值；n為實(shí)測(cè)點(diǎn)數(shù)量。

由上表數(shù)據(jù)計(jì)算可得，第一期和第二期的DEM的高程均方根誤差分別為0.236 m和0.220 m，兩期DEM的平均高程均方根誤差為0.228 m。

2)沉陷盆地精度評(píng)價(jià)。將2期水準(zhǔn)高程相減得到的下沉值與2期DEM相減得到的沉陷盆地上相應(yīng)的下沉值作比較，見(jiàn)表5，得到下沉中誤差為0.081 m。

由表3—表5數(shù)據(jù)可以看出，無(wú)人機(jī)監(jiān)測(cè)的下沉值在14、15號(hào)點(diǎn)的相對(duì)誤差較大，主要由于該位置處下沉量小為厘米級(jí)，無(wú)人機(jī)監(jiān)測(cè)精度無(wú)法達(dá)到要求，因此誤差較大；而7～11號(hào)點(diǎn)處則由于一些非采動(dòng)因素(農(nóng)民挖沙子)的影響，造成較大誤差。

表5 沉陷盆地誤差統(tǒng)計(jì)Table 5 Error statistics of subsidence basin

對(duì)下沉差值的均方根誤差進(jìn)行計(jì)算，以此評(píng)估沉陷盆地的精度，均方根誤差計(jì)算公式同式(1)。

根據(jù)表5中數(shù)據(jù)計(jì)算可得，沉陷盆地的下沉值中誤差為0.081 m。由于無(wú)人機(jī)DEM誤差中，系統(tǒng)誤差影響較大，造成無(wú)人機(jī)DEM的精度較低，而無(wú)人機(jī)監(jiān)測(cè)的沉陷盆地是由2期DEM相減獲得，消除了系統(tǒng)誤差的影響，使其只含隨機(jī)誤差，故無(wú)人機(jī)攝影測(cè)量解算的沉陷盆地的精度相對(duì)較高。

2.3 UAV-攝影測(cè)量監(jiān)測(cè)礦區(qū)沉陷可行性

由上述精度分析可知，本次試驗(yàn)中的UAV求取的下沉盆地精度為81 mm，就單個(gè)點(diǎn)的精度而言，不滿足礦區(qū)開(kāi)采沉陷變形監(jiān)測(cè)精度10 mm的要求而無(wú)法直接用于礦區(qū)地表沉陷監(jiān)測(cè)，尤其在對(duì)建筑物損害鑒定時(shí)，需要觀測(cè)沉陷盆地邊界的微小變形(監(jiān)測(cè)精度為毫米級(jí))。然而，礦區(qū)沉陷監(jiān)測(cè)的另一個(gè)核心作用是求取概率積分參數(shù)，與傳統(tǒng)沉陷監(jiān)測(cè)技術(shù)的線狀觀測(cè)站不同，UAV攝影測(cè)量技術(shù)獲取的是整個(gè)面域的沉陷值，10 mm的精度要求是針對(duì)常規(guī)的單個(gè)點(diǎn)而言的，盡管單點(diǎn)監(jiān)測(cè)精度不高，但利用面域內(nèi)大量的沉陷數(shù)據(jù)求取的預(yù)計(jì)參數(shù)具有較高的精度，以此對(duì)整個(gè)地表沉陷進(jìn)行預(yù)測(cè)，可以得到比較可靠的結(jié)果。

根據(jù)文獻(xiàn)[14]中的測(cè)量誤差對(duì)開(kāi)采沉陷預(yù)計(jì)參數(shù)的影響分析結(jié)果可知，以測(cè)量中誤差和最大下沉值之比作為衡量標(biāo)準(zhǔn)，比值<7%時(shí)求參結(jié)果可靠，比值>10%時(shí)求參結(jié)果不可靠，比值在兩者之間時(shí)求參結(jié)果部分可靠。根據(jù)計(jì)算可知，無(wú)人機(jī)監(jiān)測(cè)的某煤礦2S201工作面的測(cè)量中誤差81 mm與最大下沉值2 487 mm的比值為3.3%，因此認(rèn)為此次反演求取的沉陷參數(shù)是可靠的。由于2S201工作面在2018-6-9—2018-9-4的開(kāi)采引起的地表沉陷需持續(xù)一段時(shí)間才能停止，尚未形成穩(wěn)定的沉陷盆地，且此后工作面仍需繼續(xù)推進(jìn)，因此此期間形成的沉陷盆地為“動(dòng)態(tài)沉陷盆地”。故對(duì)該期間的沉陷數(shù)據(jù)利用概率積分法的穩(wěn)態(tài)預(yù)計(jì)模型和時(shí)間影響函數(shù)(Knothe時(shí)間影響函數(shù))進(jìn)行動(dòng)態(tài)求參。對(duì)UAV監(jiān)測(cè)數(shù)據(jù)與傳統(tǒng)水準(zhǔn)測(cè)量獲取的沉陷數(shù)據(jù)的求參結(jié)果進(jìn)行對(duì)比，見(jiàn)表6，可得，下沉系數(shù)相對(duì)誤差為1.4%，主要影響角正切相對(duì)誤差為20%。

表6 開(kāi)采沉陷預(yù)計(jì)參數(shù)反演結(jié)果Table 6 Inversion results of mining subsidence prediction parameters

利用求參結(jié)果進(jìn)行預(yù)計(jì)得到擬合下沉值，實(shí)測(cè)下沉值與擬合下沉值的中誤差(即擬合中誤差)應(yīng)滿足小于10%W0(W0為最大下沉值)的要求，否則求出的參數(shù)結(jié)果將會(huì)不可靠。此外，擬合中誤差與測(cè)量中誤差成正比，測(cè)量中誤差越大，則擬合中誤差也越大。擬合中誤差計(jì)算公式為

(2)

式中：v為實(shí)測(cè)下沉值與擬合下沉值的差。

根據(jù)上表中數(shù)據(jù)計(jì)算得，無(wú)人機(jī)數(shù)據(jù)求參結(jié)果的擬合中誤差為181 mm，占最大下沉值2 487 mm的百分比為7.3%，故認(rèn)為該求參結(jié)果相對(duì)可靠；而水準(zhǔn)測(cè)量的測(cè)量精度高于無(wú)人機(jī)測(cè)量，故水準(zhǔn)數(shù)據(jù)求參結(jié)果也可靠。

由表6中的對(duì)比結(jié)果可以看出，由無(wú)人機(jī)攝影測(cè)量獲得的沉陷數(shù)據(jù)反演得到的下沉系數(shù)是可靠的；而由于無(wú)人機(jī)監(jiān)測(cè)沉陷盆地邊界誤差較大，導(dǎo)致求取的主要影響角正切值偏小。因此無(wú)人機(jī)攝影測(cè)量技術(shù)可應(yīng)用于礦區(qū)開(kāi)采沉陷監(jiān)測(cè)，并可求取可靠的下沉系數(shù)。

礦區(qū)生態(tài)環(huán)境監(jiān)測(cè)以及后續(xù)的生態(tài)修復(fù)和土地復(fù)墾等也是礦區(qū)沉陷監(jiān)測(cè)的主要目的之一，無(wú)人機(jī)攝影測(cè)量技術(shù)以其低成本、高靈活度、高效率等優(yōu)點(diǎn)，通過(guò)數(shù)據(jù)處理獲取高分辨率地表DEM和全盆地沉陷數(shù)據(jù)以及大比例尺地形圖，為復(fù)墾規(guī)劃提供了基礎(chǔ)信息[15]。此外，獲取的海量的光學(xué)影像數(shù)據(jù)可以提取礦區(qū)地表植被信息，為研究煤礦開(kāi)采引起的礦區(qū)生態(tài)環(huán)境變化提供支撐，彌補(bǔ)常規(guī)監(jiān)測(cè)手段僅能獲得地表沉降數(shù)據(jù)的不足。

3 討論及建議

1)煤炭開(kāi)采引起巖層及地表變形問(wèn)題在生態(tài)環(huán)境極其脆弱的西部地區(qū)尤為突出，嚴(yán)重影響和制約了煤礦企業(yè)的可持續(xù)發(fā)展和人民的健康生活[2]。目前常規(guī)地表移動(dòng)觀測(cè)站布設(shè)在工作面主斷面上，由“點(diǎn)”組成“十字形”觀測(cè)線，具有尺度小、工作量大、數(shù)據(jù)及信息量少等不足。在西部淺埋深、大采高、快速綜合機(jī)械化及智能化的大規(guī)模、高強(qiáng)度開(kāi)采已成為主流下，開(kāi)采引起的巖層及地表移動(dòng)呈現(xiàn)出損害面積大、程度深、波及范圍廣的特點(diǎn)。常規(guī)地表移動(dòng)觀測(cè)站已無(wú)法適應(yīng)新形勢(shì)下的開(kāi)采損害監(jiān)測(cè)任務(wù)；UAV攝影測(cè)量為解決上述問(wèn)題提供了新的技術(shù)途徑。UAV攝影測(cè)量以無(wú)人機(jī)為平臺(tái)，可快速獲取含有豐富的地物光譜信息的光學(xué)影像，具有機(jī)動(dòng)靈活、高效精準(zhǔn)、作業(yè)成本低(突出的時(shí)效性和性價(jià)比)等優(yōu)勢(shì)；通過(guò)研究，UAV攝影測(cè)量可以快速獲得“面狀”全盆地沉陷數(shù)據(jù)，并可反演可靠的沉陷參數(shù)，可為礦區(qū)生態(tài)環(huán)境監(jiān)測(cè)以及后續(xù)的生態(tài)修復(fù)和土地復(fù)墾提供豐富數(shù)據(jù)。

2)對(duì)UAV攝影測(cè)量獲得的DEM以及沉陷盆地的精度進(jìn)行評(píng)估，得到UAV監(jiān)測(cè)的礦區(qū)地表DEM高程中誤差為228 mm，主要誤差來(lái)源于無(wú)人機(jī)平臺(tái)較輕易受天氣因素影響而導(dǎo)致影像質(zhì)量差，影像畸變大以及無(wú)人機(jī)飛行姿態(tài)不穩(wěn)定，這些誤差包括系統(tǒng)和偶然誤差；沉陷盆地通過(guò)DEM相減，消除了系統(tǒng)誤差的影響，沉陷盆地的精度為81 mm(即為厘米級(jí)精度)，比DEM高程精度提高了64.5%；因此，無(wú)人機(jī)監(jiān)測(cè)礦區(qū)大變形具有顯著優(yōu)勢(shì)，但無(wú)法準(zhǔn)確監(jiān)測(cè)沉陷盆地邊界區(qū)域小變形。就單個(gè)點(diǎn)的精度而言，無(wú)人機(jī)攝影測(cè)量監(jiān)測(cè)精度不滿足要求；然而，從開(kāi)采沉陷參數(shù)反演的角度來(lái)說(shuō)，與水準(zhǔn)測(cè)量求參結(jié)果相比，無(wú)人機(jī)攝影測(cè)量獲得沉陷數(shù)據(jù)求取的下沉系數(shù)相對(duì)誤差較小，僅為1.4%，由于UAV獲得下沉盆地邊界誤差較大，求取的主要影響角正切相對(duì)誤差為20%，誤差偏大。UAV攝影測(cè)量的不足之處在于：①無(wú)人機(jī)平臺(tái)較輕、攝影姿態(tài)不穩(wěn)定(影像旋偏較大)；②飛行高度低、影像幾何畸變大；③航攝影像重疊度不規(guī)則等方面的缺陷，導(dǎo)致UAV攝影測(cè)量的高程測(cè)量精度相對(duì)較低。根據(jù)實(shí)測(cè)結(jié)果，平均高程精度為228 mm(即為分米級(jí)精度)。UAV監(jiān)測(cè)礦區(qū)開(kāi)采沉陷大變形具有優(yōu)勢(shì)，卻無(wú)法準(zhǔn)確獲得沉陷邊界區(qū)域數(shù)據(jù)，這也導(dǎo)致其反演的參數(shù)主要影響角正切誤差較大的原因。

3)盡管無(wú)人機(jī)攝影測(cè)量技術(shù)在礦區(qū)地表沉陷中存在些許不足，但其監(jiān)測(cè)得到的高分辨率影像以及處理生成的各數(shù)字產(chǎn)品對(duì)于礦區(qū)生態(tài)環(huán)境監(jiān)測(cè)具有重要作用。為推廣無(wú)人機(jī)攝影測(cè)量技術(shù)在礦區(qū)中的應(yīng)用，需提高其監(jiān)測(cè)精度，可以通過(guò)傳統(tǒng)觀測(cè)站與無(wú)人機(jī)攝影測(cè)量技術(shù)聯(lián)合監(jiān)測(cè)的方法；根據(jù)文獻(xiàn)[16]中測(cè)點(diǎn)缺失對(duì)地表移動(dòng)參數(shù)的影響的研究可知，重要觀測(cè)點(diǎn)(如最大下沉處、拐點(diǎn)附近等)對(duì)概率積分法參數(shù)的確定起決定性作用，為彌補(bǔ)無(wú)人機(jī)攝影測(cè)量技術(shù)存在的缺陷，可以將觀測(cè)站布設(shè)在礦區(qū)地表最大下沉處、拐點(diǎn)附近以及盆地邊界；將觀測(cè)站數(shù)據(jù)與無(wú)人機(jī)數(shù)據(jù)聯(lián)合處理，既可得到完整可靠的沉陷盆地，又能可靠地求得開(kāi)采沉陷預(yù)計(jì)的各參數(shù)。此外，隨著無(wú)人機(jī)制造工藝、機(jī)載定位、視覺(jué)算法等技術(shù)的創(chuàng)新，無(wú)人機(jī)攝影測(cè)量技術(shù)正趨向于智能化和自動(dòng)化發(fā)展[17]；研制高效鋰電池、太陽(yáng)能電池、無(wú)人機(jī)群體化協(xié)同作業(yè)等方法提高了無(wú)人機(jī)的續(xù)航能力，使無(wú)人機(jī)攝影測(cè)量技術(shù)更能適應(yīng)于礦山的大規(guī)模監(jiān)測(cè)；無(wú)人機(jī)與相機(jī)設(shè)備的進(jìn)步，使由無(wú)人機(jī)平臺(tái)不穩(wěn)定、影像畸變等問(wèn)題造成的誤差可以減小。

4)除了提到的常規(guī)測(cè)量技術(shù)和UAV攝影測(cè)量外，目前地面三維激光掃描技術(shù)，InSAR技術(shù)也被應(yīng)用礦區(qū)地表沉陷監(jiān)測(cè)中。然而，地面三維激光掃描技術(shù)雖然可以獲得大量空間點(diǎn)云信息，從而構(gòu)建三維模型[18-21]，但其監(jiān)測(cè)成本高且需多次遷站，不具有靈活性；InSAR技術(shù)可全天時(shí)、全天候連續(xù)監(jiān)測(cè)，且監(jiān)測(cè)成本低[22-24]，但時(shí)空相位失相干阻礙了它在礦區(qū)地表沉陷監(jiān)測(cè)中的推廣應(yīng)用[25]。而無(wú)人機(jī)攝影測(cè)量技術(shù)以其成本低、監(jiān)測(cè)范圍廣、靈活度高、可連續(xù)性監(jiān)測(cè)等優(yōu)勢(shì)，逐漸在礦區(qū)地表沉陷監(jiān)測(cè)中被廣泛應(yīng)用。隨著測(cè)繪新技術(shù)的發(fā)展，多種技術(shù)融合，多源數(shù)據(jù)融合是開(kāi)采沉陷監(jiān)測(cè)的發(fā)展趨勢(shì)，通過(guò)多技術(shù)融合可以充分發(fā)揮各自優(yōu)勢(shì)，揚(yáng)長(zhǎng)避短獲得更高精度的全盆地沉陷數(shù)據(jù)，可以精細(xì)化研究地表盆地變形。通過(guò)多源數(shù)據(jù)融合即可以獲得地表變形，又可以精確獲得礦區(qū)生態(tài)環(huán)境及其變化，為礦區(qū)生態(tài)修復(fù)、環(huán)境評(píng)價(jià)及后評(píng)價(jià)提供高精度豐富的基礎(chǔ)數(shù)據(jù)。

4 結(jié) 論

1)分析了西部大多數(shù)礦區(qū)煤層淺埋，大采高、高強(qiáng)度快速開(kāi)采引發(fā)的地表沉陷呈現(xiàn)變形速度快、損害程度深、波及范圍廣的特點(diǎn)。常規(guī)觀測(cè)站的不足導(dǎo)致其已無(wú)法適應(yīng)高強(qiáng)度開(kāi)采地表?yè)p害監(jiān)測(cè)任務(wù)。在這種形式下，如何快速、準(zhǔn)確、全面地監(jiān)測(cè)西部高強(qiáng)度大規(guī)模開(kāi)采引起的地表沉陷與環(huán)境損害是解決問(wèn)題的關(guān)鍵。提出采用無(wú)人機(jī)攝影測(cè)量技術(shù)，給出了無(wú)人機(jī)攝影測(cè)量監(jiān)測(cè)礦區(qū)開(kāi)采損害的思路和方法，并以內(nèi)蒙古鄂爾多斯某煤礦為例進(jìn)行了應(yīng)用研究，結(jié)果表明無(wú)人機(jī)攝影測(cè)量技術(shù)可以獲取“面狀”全盆地?cái)?shù)據(jù)并反演出穩(wěn)定可靠的參數(shù)，具有簡(jiǎn)單、高效、快速的特點(diǎn)，彌補(bǔ)了常規(guī)監(jiān)測(cè)方法的不足。

2)以某煤礦為研究對(duì)象，利用無(wú)人機(jī)攝影測(cè)量技術(shù)獲取地表影像數(shù)據(jù)，得到監(jiān)測(cè)區(qū)域高分辨率的DEM數(shù)據(jù)和沉陷盆地。通過(guò)與全站儀/水準(zhǔn)測(cè)量的地面觀測(cè)站數(shù)據(jù)對(duì)比，對(duì)無(wú)人機(jī)DEM和沉陷盆地的精度進(jìn)行評(píng)估，精度分別為228 mm和81 mm，基于精度分析結(jié)果對(duì)無(wú)人機(jī)攝影測(cè)量技術(shù)監(jiān)測(cè)礦區(qū)地表沉陷的可行性進(jìn)行分析，得出其在開(kāi)采沉陷預(yù)計(jì)參數(shù)反演中，可以求得較為可靠的下沉系數(shù)；獲取的大量影像數(shù)據(jù)可以為礦區(qū)生態(tài)環(huán)境監(jiān)測(cè)以及后續(xù)的生態(tài)修復(fù)和土地復(fù)墾提供基礎(chǔ)數(shù)據(jù)。

3)無(wú)人機(jī)攝影測(cè)量技術(shù)作為一種新技術(shù)應(yīng)用于煤礦開(kāi)采地表沉陷中，具有監(jiān)測(cè)成本低，耗時(shí)短等優(yōu)勢(shì)，結(jié)合某煤礦研究實(shí)例來(lái)看，其監(jiān)測(cè)精度81 mm達(dá)厘米級(jí)，但對(duì)于單點(diǎn)精度來(lái)說(shuō)尚待提高，可對(duì)其誤差源進(jìn)行分析以改善其精度，使之在煤礦開(kāi)采沉陷監(jiān)測(cè)中得到推廣應(yīng)用。

4)討論分析了無(wú)人機(jī)攝影測(cè)量及現(xiàn)有監(jiān)測(cè)手段的優(yōu)缺點(diǎn)，建議采用傳統(tǒng)觀測(cè)站與無(wú)人機(jī)攝影測(cè)量技術(shù)聯(lián)合監(jiān)測(cè)的方法，可以得到更為精準(zhǔn)的沉陷盆地信息和參數(shù)。隨著測(cè)繪新技術(shù)的發(fā)展，多種技術(shù)融合發(fā)揮各自優(yōu)勢(shì)，多源數(shù)據(jù)融合進(jìn)行精準(zhǔn)礦區(qū)生態(tài)環(huán)境損害監(jiān)測(cè)是開(kāi)采沉陷監(jiān)測(cè)的發(fā)展趨勢(shì)。