朱鵬程

（1. 陜西鐵路工程職業(yè)技術(shù)學(xué)院，陜西 渭南 714000）

煤炭資源作為我國主要工業(yè)能源之一，在為經(jīng)濟建設(shè)做出重要貢獻的同時也給環(huán)境帶來了嚴(yán)重破壞，主要包括工業(yè)污染[1]、開采沉陷[2]、水土流失[3]等。由于井下采煤導(dǎo)致煤層上覆巖層發(fā)生斷裂與位移，且隨著開采范圍的增大，巖層移動波及地表，引起地表位移，不可避免地對采動影響區(qū)域內(nèi)的建（構(gòu)）筑物、耕地、人類生產(chǎn)活動產(chǎn)生影響。為定量分析采動對地表的影響以及二者之間的對應(yīng)關(guān)系，相關(guān)學(xué)者開展了大量研究，如李昱昊[4]等以無人機攝影測量技術(shù)為基礎(chǔ)，利用BP 神經(jīng)網(wǎng)絡(luò)去除了沉降盆地內(nèi)的噪點，進而提高了下沉盆地監(jiān)測精度；廉旭剛[5]等研究了免像控?zé)o人機攝影測量技術(shù)監(jiān)測開采沉陷的精度；楊緒霆[6]等通過地形跟隨模式提高了無人機攝影測量方法監(jiān)測礦區(qū)沉陷的精度。除無人機監(jiān)測外，InSAR 技術(shù)也被廣泛應(yīng)用于地表移動監(jiān)測，如譚志祥[7]、沙永蓮[8]等采用SBAS-InSAR 技術(shù)實現(xiàn)了沉陷盆地的毫米級監(jiān)測；張童康[9]等將InSAR 與支持向量機相結(jié)合，改進了開采沉陷預(yù)測模型。三維激光掃描技術(shù)也是常見的開采沉陷監(jiān)測方法，如王磊[10]、賈秉松[11]等分別將其應(yīng)用于礦區(qū)建筑物的變形監(jiān)測和采動損害評估以及受采動影響的山區(qū)邊坡沉陷監(jiān)測。已有研究幾乎均是通過某種技術(shù)手段監(jiān)測地表某點的縱向位移，即單點的沉降量，但地表某點受采動影響的位移是一個三維移動的過程[12]，既包含縱向位移，又包含橫向位移（水平位移），二者疊加才是其真實的位移量。為了彌補已有監(jiān)測方案的缺陷，更加全面地描述地表受采動影響的變化過程，本文通過無人機攝影測量技術(shù)與計算機視覺技術(shù)相結(jié)合的方法，實現(xiàn)了采動影響區(qū)域內(nèi)地表的水平移動監(jiān)測，可為井下開采、沉陷預(yù)計、人員安置等提供更全面的數(shù)據(jù)。

1 無人機監(jiān)測礦區(qū)地表水平移動方法

本文提出的無人機礦區(qū)地表水平移動監(jiān)測方法的技術(shù)流程見圖1，首先利用開采沉陷預(yù)計方法預(yù)計采動影響范圍，并根據(jù)監(jiān)測精度確定無人機作業(yè)的航高、重疊度、航線等參數(shù)；然后采集預(yù)計影響范圍內(nèi)不同時期的時序影像，通過Pix4Dmapper處理影像獲得時序DOM；再利用Global Mapper 對DOM 進行網(wǎng)格化處理，并進行特征點匹配，計算特征點在不同時序影像上的位置；最后計算特征點的水平移動距離與方向，繪制水平移動矢量圖。

圖1 無人機礦區(qū)地表水平移動監(jiān)測方法流程圖

2 水平移動提取關(guān)鍵技術(shù)

2.1 網(wǎng)格化處理

原始影像經(jīng)Pix4Dmapper 處理后生成DOM，再利用Global Mapper 進行網(wǎng)格化處理，具體流程見圖2，網(wǎng)格大小由沉陷預(yù)計的最大水平移動量ΔDmax確定，假設(shè)ΔDmax=500 mm，則網(wǎng)格邊長為1.2ΔDmax即600 mm。劃分網(wǎng)格大小的目的：①減小特征點匹配時的運算量，只需在對應(yīng)網(wǎng)格內(nèi)匹配同名特征點即可；②確定時序移動增量監(jiān)測范圍，隨著采動的影響，地表發(fā)生水平移動與沉降后，其特征會與初期DOM 特征差異巨大，導(dǎo)致無法匹配，劃定網(wǎng)格后只需匹配當(dāng)前網(wǎng)格與上期對應(yīng)網(wǎng)格，通過計算移動增量并累加增量的方法獲得當(dāng)期與初期之間的水平移動大小與方向。

圖2 網(wǎng)格化處理

2.2 特征點匹配

常見的特征點匹配算法包括Harris 角點檢測算法、SIFT 算法、SURF 算法、FAST 角點檢測算法、BRIEF 描述子以及基于FAST 角點檢測算法和BRIEF描述子相結(jié)合的ORB算法等。大量學(xué)者已將上述算法應(yīng)用于各場景中并做出了相關(guān)改進，如王丞[13]等利用Harris 角點檢測算法實現(xiàn)了三維點云的自適應(yīng)特征描述、提取與粗配準(zhǔn)；張占平[14]等將Harris 算子與SIFT描述子相結(jié)合，并輔以無人機影像的POS數(shù)據(jù)與公開的SRTM 數(shù)據(jù)，提高了無人機傾斜攝影影像的匹配效率；徐啟文[15]等利用一種具有動態(tài)閾值的改進SURF算法提升了圖像拼接效率；陳偉[16]等結(jié)合FAST提取的特征點和SURF 算法的描述子實現(xiàn)了圖像的快速拼接；陶卓[17]等利用稀疏光流法改進了傳統(tǒng)ORB 算法，使算法的匹配精度更高、魯棒性更強。本文對同一場景進行了常見特征檢測與提取算法的測試，并統(tǒng)計了5種特征點提取算法的提取效果，結(jié)果見表1。通過對比提取效率、特征點提取數(shù)量以及特征點匹配準(zhǔn)確性，最終確定本文的特征點提取與匹配采用經(jīng)典SIFT算法。

表1 常見特征檢測與提取算法比較

2.3 水平移動大小與方向的計算

利用SIFT算子匹配兩幅影像中同名特征點后，需將特征點像素坐標(biāo)轉(zhuǎn)換為地理坐標(biāo)，本文采用GDAL空間數(shù)據(jù)轉(zhuǎn)換坐標(biāo)庫進行坐標(biāo)轉(zhuǎn)換，轉(zhuǎn)換公式為：

式中，trans(0)、trans(3)為影像左上角像素點左上角位置的地理坐標(biāo)；trans(1)、trans(5)為影像上單個像素的寬度與高度；如若該影像是指北的，則trans(2)、trans(4)均為0；m、n為特征點在影像中的像素坐標(biāo)；x、y為該特征點轉(zhuǎn)換后的地理坐標(biāo)。

將同名特征點的像素坐標(biāo)轉(zhuǎn)換為地理坐標(biāo)后，便可計算該特征點的水平移動量和方向，假設(shè)同一網(wǎng)格匹配n對特征點，各對特征點水平移動量的計算公式為：

為使水平移動監(jiān)測更加可視化，可計算水平移動方向，假定某點水平移動是由(x1'y1)指向(x2y2)，且水平移動矢量與東方向形成的夾角為α（圖3），則有：

圖3 水平移動矢量示意圖

2.4 基于水平移動大小與方向的誤匹配剔除

由于礦區(qū)地物特征的特殊性和復(fù)雜性，采用SIFT算子提取的匹配特征點存在誤差（圖4a），為提高水平移動監(jiān)測精度，本文提出了一種基于位移大小和方向的誤匹配剔除方法。通常在同一網(wǎng)格內(nèi)，各特征點間水平移動量差別不大，且水平移動方向具有一致性，根據(jù)這兩個特性可剔除誤匹配和不符合地表水平移動規(guī)律的匹配。首先計算同一網(wǎng)格內(nèi)所有水平移動量的平均值ΔA，設(shè)置0.9ΔA～1.1ΔA為滿足誤匹配剔除的大小閾值范圍；再計算同一網(wǎng)格內(nèi)所有水平移動方向的平均值Δα，設(shè)置0.9ΔA～1.1ΔA為滿足誤匹配剔除的角度閾值范圍；最后結(jié)合兩個閾值范圍，可基本剔除誤匹配，剔除后效果見圖4b。

圖4 誤匹配剔除前后對比

3 實測案例分析

3.1 案例概況

監(jiān)測案例中采動區(qū)域為山東省濟寧市南屯煤礦33上08工作面（圖5），測區(qū)中心位于35°22′43″N、116°52′01″E，平均采深為300 m、走向長680 m、傾向長170 m、平均煤厚5.2 m。該工作面采動前已對地表村莊易地搬遷，但仍保留部分工業(yè)建筑，且該區(qū)域為高潛水位地區(qū)，采動后地下水會出露地表，因此對該區(qū)域的地表移動進行監(jiān)測非常必要。

圖5 監(jiān)測區(qū)域

基于概率積分法的開采沉陷預(yù)計方法是目前較成熟的礦山開采沉陷預(yù)計手段，主要基于已有變形監(jiān)測數(shù)據(jù)與對應(yīng)工作面參數(shù)，根據(jù)隨機介質(zhì)理論，以統(tǒng)計學(xué)的視角將采動區(qū)域分解為無限個微小開采單元，進而反演計算工作面參數(shù)與變形監(jiān)測數(shù)據(jù)之間的相互關(guān)系，主要通過下沉系數(shù)、水平移動系數(shù)、拐點偏移距、主要影響角正切以及開采影響傳播角等預(yù)計參數(shù)進行描述和計算，并以此預(yù)計同地質(zhì)采礦條件下的其他待采動工作面的變形情況[18]。由3308工作面參數(shù)與歷史預(yù)計參數(shù)計算可知，該工作面預(yù)計最大下沉值為4 549 mm，東西最大水平移動為1 396 mm，南北最大水平移動為1 455 mm（圖6）。

圖6 33上08工作面開采地表下沉等值線

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

為監(jiān)測3308工作面采動對地表的影響，檢驗本文方法的有效性，2021 年3—4 月對3308 工作面采動影響區(qū)域采集3 期無人機影像數(shù)據(jù)與像控點坐標(biāo)數(shù)據(jù)。無人機采用大疆精靈RTK，作業(yè)高度為100 m，影像分辨率為2.73 cm/pixel，具體飛行數(shù)據(jù)見表2。像控點采用60 cm×60 cm 倒三角形像控點（圖7），并采用GPS-RTK+DS3 級水準(zhǔn)測量獲得坐標(biāo)與高程，其中平面誤差為2 cm，高程誤差為1 mm，共布設(shè)35 個像控點（圖8），包括10 個像控點和25 個檢查點。

表2 無人機飛行數(shù)據(jù)

圖7 33上08工作面像控點布設(shè)與量測

圖8 像控點布設(shè)位置分布圖

3.3 數(shù)據(jù)處理與精度分析

利用本文方法對3 期數(shù)據(jù)進行處理，獲得水平移動大小與方向圖（圖9），并將實測值與提取值進行比較，統(tǒng)計本文方法的水平移動提取誤差（表3），首期水平移動由二期DOM 與初期DOM 比較獲得，二期水平移動由三期DOM 與二期DOM 比較獲得，經(jīng)計算首期水平移動提取中誤差為1.9 cm，二期水平移動提取中誤差為3.3 cm，綜合誤差為3.8 cm。

表3 水平移動提取誤差統(tǒng)計

圖9 兩期地表移動監(jiān)測結(jié)果

4 結(jié) 語

本文基于無人機攝影測量技術(shù)提出了一種礦區(qū)地表水平移動監(jiān)測方法，并通過實測案例檢驗了該方法的可行性。

1）地表水平移動作為受采動影響的地表變形的重要組成部分，前人的相關(guān)研究重點關(guān)注地表沉降監(jiān)測，而對水平移動監(jiān)測的研究存在缺失，本文是對相關(guān)研究的補充。

2）本文方法可基本完成地表水平移動的自動化監(jiān)測，較少需要人為干預(yù)，且精度可滿足礦區(qū)日常監(jiān)測要求。在影像分辨率為2.73 cm/pixel 的前提下，水平移動提取精度為3.8 cm，且隨著硬件設(shè)備性能的提升，該方法可獲得更高精度的結(jié)果。

本文為受采動影響的地表水平移動監(jiān)測提供了一種更全面、便利、可視化的監(jiān)測手段，可廣泛應(yīng)用于采動區(qū)域的變形監(jiān)測，為井下開采、耕地保護和搬遷移民等提供前瞻性和實時的數(shù)據(jù)支持，對于井下安全開采、采動區(qū)治理以及采動影響范圍內(nèi)的生命財產(chǎn)安全等具有重大意義。