張曉巖，路 鑫

(潞安化工集團(tuán)有限公司 古城煤礦，山西 長治 046000)

古城煤礦自投產(chǎn)至今，全面的地表沉陷監(jiān)測數(shù)據(jù)相對缺乏，工作面開采所引起的地表變形規(guī)律不夠明確，對厚煤層綜放開采條件下地面沉陷機理的研究不夠深入，從而導(dǎo)致適用于本礦區(qū)的地表移動角值參數(shù)及開采沉陷預(yù)計參數(shù)尚不完備。高精度的地面水準(zhǔn)和位移監(jiān)測是研究地表移動變形規(guī)律的重要手段，但傳統(tǒng)的“十字線”地表移動觀測站只能獲取主斷面上的“點狀”信息，不能完全反映整個盆地的移動變形特征，導(dǎo)致所求取的角量參數(shù)具有不確定性。

本文以古城煤礦S1306工作面為觀測工作面，對地表沉陷“天-地”(“天”指無人機低空雷達(dá)測量；“地”指傳統(tǒng)地面測量)協(xié)同監(jiān)測的技術(shù)方案進(jìn)行設(shè)計，以獲取高精度、面狀的地表變形信息。研究成果可為后續(xù)沉陷規(guī)律研究和預(yù)計參數(shù)反演提供數(shù)據(jù)基礎(chǔ)，亦對指導(dǎo)保護(hù)煤柱留設(shè)和“三下”壓煤開采、對礦區(qū)生態(tài)修復(fù)及環(huán)境后評價[1-3]。

1 工程概況

S1306工作面位于古城煤礦南部一采區(qū)，地勢平坦，地表多為農(nóng)田；周邊村莊和公路較多，東以常廟線農(nóng)村公路為界，南臨西酪余村，西側(cè)有王史村，北臨東李高村。工作面走向長約2 100 m，傾向?qū)捈s360 m，煤層平均厚度6.15 m，傾角約為5°，平均采深536 m；松散層厚度約90 m.綜合放頂煤開采，全部垮落法管理頂板。工作面井上下對照情況如圖1所示。

圖1 S1306工作面井上下對照圖及測點布置圖

2 地表巖移觀測方案

2.1 地表巖移觀測站設(shè)計

2.1.1 設(shè)計參數(shù)選取

根據(jù)礦區(qū)生產(chǎn)資料，結(jié)合本礦S1306工作面具體情況，并參考周邊煤礦(五陽煤礦)的地表實測沉陷參數(shù)，選取S1306工作面巖移觀測站的設(shè)計參數(shù)如表1所示。

表1 S1306工作面地表移動觀測站設(shè)計參數(shù)

為確保觀測線沿煤柱一側(cè)的邊界不受開采影響，觀測線應(yīng)有足夠的長度。因此，按煤礦測量規(guī)程取各種移動角的調(diào)整值：Δγ=Δβ=Δδ=20°.

2.1.2 觀測站位置設(shè)計

1) 傾向觀測線位置設(shè)計。為保證傾向觀測線通過充分采動區(qū)，觀測線到開切眼的距離D必須滿足公式(1)[4]：

D≥(H0-h)cot(δ-Δδ)+hcotφ

(1)

式中：H0為工作面平均采深；h為松散層厚度，取90 m；δ為走向移動角；Δδ為走向移動角修正值；φ為松散層移動角。

經(jīng)計算得：D≥399 m，即傾向觀測線與開切眼的距離必須大于399 m.

2) 走向觀測線位置設(shè)計。S1306工作面煤層傾角α=5°，為盡可能保證走向觀測線通過最大下沉點，走向觀測線應(yīng)從工作面中間向下山方向偏移。根據(jù)五陽煤礦實測最大下沉角，S1306工作面下山方向偏移量P應(yīng)滿足公式(2)：

P=(H0-h)cotθ

(2)

式中：H0為回采工作面平均開采深度；h為松散層厚度；θ為最大下沉角。

將參數(shù)值帶入公式(2)計算得：P=39 m.此外，由于S1306工作面傾向?qū)?60 m<1.4H0，傾向方向未達(dá)到充分采動，所以走向觀測線設(shè)置為一條。

2.1.3 觀測線長度設(shè)計

1) 傾向觀測線長度設(shè)計。傾向觀測線長度必須保證覆蓋上山盆地半長、下山盆地半長和工作面的寬度，并有一定的安全距離，因回采工作面傾向方向未達(dá)到充分采動，所以布置整條觀測線，傾向觀測線最小長度可按公式(3)計算[5]:

L傾向=2hcotφ+(H1-h)cot(β-Δβ)+(H2-h)cot(γ-Δγ)+Lcosα

(3)

式中：φ為松散層移動角；H1、H2分別為上、下山方向采深；Lcosα為工作面水平投影寬度；其它符號含義與前述相同。

經(jīng)計算，得到傾向觀測線長度約為1 224 m.從與走向觀測線交叉點開始，上山方向長度約為640 m，下山方向長度為584 m.

2) 走向觀測線長度設(shè)計走向觀測線要保證停采線方向延伸出沉陷盆地邊界范圍之外，工作面之內(nèi)要測到最大下沉點，全長走向觀測線長度可用公式(4)計算[5]：

L走向=2hcotφ+2(H0-h)cot(δ-Δδ)+l1

(4)

式中：h為松散層厚度；φ為松散層移動角；H0為平均采深；δ和Δδ分別為走向移動角及其修改值；l1為工作面走向長度。

經(jīng)計算得到，走向觀測線長度約為2 898 m，其中，工作面外側(cè)約782 m，工作面內(nèi)側(cè)約2 111 m.

2.1.4 觀測點數(shù)目及密度的確定

觀測線上的測點數(shù)目及其密度，主要取決于開采深度和設(shè)站的目的。根據(jù)《煤礦測量規(guī)程》，采深大于400 m時，測點間距為35 m；控制點間距為50 m.經(jīng)計算，需布置工作測點120個，傾向觀測線上布置40個測點，走向觀測線上布置80個測點；其中，控制點9個(走向3個，傾向6個)。測點具體位置如圖1所示。

2.2 觀測內(nèi)容

地表移動觀測站的觀測工作主要包括：連接測量、全面觀測、日常觀測等。具體如下：

2.2.1 連接測量

與測區(qū)內(nèi)已知控制點進(jìn)行連測，以獲得觀測站控制點的平面位置和高程。在觀測點埋好10～15 d、點位固結(jié)之后，獨立進(jìn)行兩次，點位中誤差需小于7 cm，高程測量按三等水準(zhǔn)的精度要求進(jìn)行。

2.2.2 全面測量

全面觀測是確定各測點的初始位置和高程。需連接測量后、地表開始移動之前，獨立進(jìn)行兩次全面觀測，兩次觀測時間間隔小于5 d.兩次觀測的同一點高程差小于10 mm，同一邊的長度差不大于4 mm時，取平均值作為觀測站的原始觀測數(shù)據(jù)(又稱初次觀測)。

2.2.3 日常觀測

日常觀測工作，指的是首次和末次全面觀測之間適當(dāng)增加的水準(zhǔn)測量工作。在回采工作面推進(jìn)0.2～0.5H0后，每隔4～5 d進(jìn)行一次水準(zhǔn)測量。在移動過程中，視地表下沉速度， 1～3個月觀測一次。實測按四等水準(zhǔn)測量的精度要求進(jìn)行。

觀測站的各項觀測，如表2所示。

表2 觀測站觀測程序

3 無人機激光雷達(dá)監(jiān)測方案設(shè)計

無人機(unmanned aerial vehicle,UAV)激光雷達(dá)監(jiān)測原理與地面三維激光掃描原理相同，是以無人機作為觀測平臺，對同一位置地表在不同時段進(jìn)行兩次掃描，獲得兩個時刻地表的數(shù)字高程模型(Digital Elevation Model, DEM)，用兩期DEM相減，可以得到監(jiān)測區(qū)域的地表的實測下沉值，同時根據(jù)動態(tài)求參原理，獲得地表移動變形參數(shù)[6]。其原理如圖2所示。

圖2 UAV-Lidar沉陷盆地獲取原理圖[6]

無人機飛行平臺擬選擇大疆M600Pro；激光雷達(dá)選擇智喙ARS-450i激光測量系統(tǒng)，它集成高精度激光掃描儀(測距精度5 mm)、GPS、IMU等傳感器，可獲取三維激光點云和定位定姿數(shù)據(jù)；通過配備的數(shù)據(jù)處理和應(yīng)用軟件，可以快速生成DEM，制作3D模型。

3.1 飛行范圍確定

工作面開采引起地表沉降的范圍即為無人機激光雷達(dá)測區(qū)飛行范圍。利用測區(qū)地質(zhì)資料及相關(guān)參數(shù)，結(jié)合收集的采掘工程平面圖、井上下對照圖等，規(guī)劃無人機地表監(jiān)測范圍為3 761 m×1 623 m，約6.1 km2.

3.2 飛行參數(shù)設(shè)置

飛行參數(shù)主要有重疊度、飛行高度。根據(jù)攝影測量規(guī)范，并結(jié)合工作面的地形條件和精度要求，設(shè)計航向重疊率和旁向重疊率為均30%.地面分辨率設(shè)置為6 cm，利用相機的焦距(15 mm)和像元尺寸，得到本次項目設(shè)計的相對飛行高度為65 m.S1306工作面無人機激光雷達(dá)設(shè)計航線如圖3所示[7-9]。

圖3 古城煤礦S1306工作面無人機激光雷達(dá)飛行范圍及航線

3.3 數(shù)據(jù)后處理

無人機外業(yè)結(jié)束后的數(shù)據(jù)后處理主要包括：POS解算、點云融合、點云濾波和DEM構(gòu)建等內(nèi)容，分述如下：

3.3.1 POS解算

主要是用基站收集的GNSS靜態(tài)數(shù)據(jù)和POS數(shù)據(jù)(移動站GNSS數(shù)據(jù)+IMU數(shù)據(jù)+里程計數(shù)據(jù))組合解算，以輸出后期數(shù)據(jù)出來所需的高精度定位定姿數(shù)據(jù)[8]。

3.3.2 點云融合及精度驗證

點云融合是將機載激光雷達(dá)系統(tǒng)采集到的原始數(shù)據(jù)進(jìn)行融合處理，生成平面坐標(biāo)的點云文件，以進(jìn)行后續(xù)的點云濾波及DEM構(gòu)建。進(jìn)行精度驗證，可保證點云數(shù)據(jù)的真實可用性。

3.3.3 點云濾波

其目的是提出三維點云數(shù)據(jù)中的非地面點，如植被點、建筑物點、車輛點、橋梁點和電力線點等，以構(gòu)建生成DEM。目前主流的點云濾波算法有不規(guī)則三角網(wǎng)漸進(jìn)加密、分層穩(wěn)健線性內(nèi)插、漸進(jìn)窗口數(shù)學(xué)形態(tài)學(xué)濾波、基于坡度的濾波方法等。

3.3.4 構(gòu)建DEM

利用前期處理得到的優(yōu)質(zhì)點云數(shù)據(jù)，采用規(guī)則格網(wǎng)或不規(guī)則三角網(wǎng)進(jìn)行測區(qū)DEM的構(gòu)建；兩期DEM相減即得到測區(qū)下沉盆地。

4 結(jié) 語

1) 本文通過設(shè)計觀測線位置、長度、數(shù)量，以及觀測方法等內(nèi)容，建立了完整且可靠的傳統(tǒng)地表巖移觀測站，有利于后續(xù)獲取各種巖移參數(shù)，總結(jié)研究區(qū)地表的移動變形規(guī)律。

2) UAV激光雷達(dá)監(jiān)測地表沉陷有重量輕、攜帶方便、成果處理效率高的優(yōu)點，因此為了更全面地研究地表沉陷規(guī)律，本文詳細(xì)設(shè)計了UAV激光雷達(dá)監(jiān)測方案，包括無人機飛行范圍劃定、飛行參數(shù)設(shè)置、實測步驟介紹以及最后數(shù)據(jù)處理流程等內(nèi)容，其中后期數(shù)據(jù)處理的精確與否是獲得巖移規(guī)律和巖移參數(shù)的關(guān)鍵。

3) 監(jiān)測礦區(qū)地表沉陷有多種方法，除了上述兩種之外，仍可以考慮利用InSAR大規(guī)模監(jiān)測不同時間跨度的地表沉陷情況。關(guān)于綜合運用“空-天-地”三種監(jiān)測技術(shù)來獲取巖移規(guī)律和地表移動參數(shù)，也是未來開采沉陷領(lǐng)域的重要趨勢。