李如仁，楊 震，余 博

(1. 沈陽建筑大學(xué)交通工程學(xué)院，遼寧 沈陽 110168； 2. 中國礦業(yè)大學(xué)(北京)地球科學(xué)與測繪工程學(xué)院，北京 100083)

GB-InSAR集成GIS的露天煤礦邊坡變形監(jiān)測

李如仁1，楊 震2，余 博1

(1. 沈陽建筑大學(xué)交通工程學(xué)院，遼寧 沈陽 110168； 2. 中國礦業(yè)大學(xué)(北京)地球科學(xué)與測繪工程學(xué)院，北京 100083)

對露天礦邊坡進(jìn)行實時、精確測量可為邊坡提供可靠的預(yù)警和治理信息，直接關(guān)系到露天礦的生產(chǎn)連續(xù)性和效益。地基SAR干涉測量技術(shù)(GB-InSAR)是最近10年發(fā)展起來的一種相對較新的變形監(jiān)測技術(shù)，具有全天候、全天時的觀測特點。而國內(nèi)針對露天礦邊坡變形監(jiān)測大多處于試驗階段，缺乏長時間連續(xù)觀測的能力，本文詳細(xì)地闡述了GB-InSAR的基本原理，并將GB-InSAR得到的變形數(shù)據(jù)與GIS有機(jī)地連接起來，利用GIS強(qiáng)大的空間顯示、分析能力，高效、實時、直觀地展示變形信息，為露天礦生產(chǎn)提供了有力的決策支持。

GB-InSAR；露天煤礦邊坡；GIS；變形監(jiān)測

隨著露天礦生產(chǎn)的不斷進(jìn)行，其邊坡角也會逐漸增加，增加了露天礦邊坡災(zāi)害發(fā)生的概率。因此，對露天礦邊坡進(jìn)行實時的變形監(jiān)測，掌握邊坡災(zāi)害的過程，利用一定的數(shù)據(jù)處理方法分析和建立預(yù)報、預(yù)警系統(tǒng)是至關(guān)重要的[1]。

全站儀、測量機(jī)器人和水準(zhǔn)儀等傳統(tǒng)測量技術(shù)在變形監(jiān)測中已經(jīng)得到廣泛的應(yīng)用，但容易受氣候、時間、通視等條件的限制；GPS克服了傳統(tǒng)觀測技術(shù)的眾多缺點，在變形監(jiān)測中取得了廣泛的應(yīng)用，但露天礦邊坡的坡度往往比較大，可視衛(wèi)星的數(shù)目及幾何位置嚴(yán)重影響了GPS的精度甚至可行性。GPS和全站儀等傳統(tǒng)測量技術(shù)能獲得相對較高精度的單點位移[2]，但難以整體反映邊坡變形，更重要的是這些獲取變形信息的方式都是接觸式的，對于危險的邊坡不宜進(jìn)行觀測。

合成孔徑雷達(dá)干涉測量技術(shù)(synthetic aperture radar interferometry，InSAR)是近些年發(fā)展起來的極具潛力的微波遙感新技術(shù)，它已在冰川漂移、地震變形、火山運動等領(lǐng)域表現(xiàn)出很好的前景[3-6]。GB-InSAR是在星載SAR的基礎(chǔ)上發(fā)展起來的一項新技術(shù)，其具有靈活多變、分辨率高、平臺穩(wěn)定、造價相對低廉等優(yōu)點。GB-InSAR不僅能夠監(jiān)測山體滑坡、冰川移動、雪崩等無法接觸或較為危險的領(lǐng)域，也能在大型橋梁、高層建筑物的變形監(jiān)測中發(fā)揮重大的作用[7]，國內(nèi)利用GB-InSAR進(jìn)行露天礦變形監(jiān)測僅僅停留在試驗階段，缺乏長時間連續(xù)監(jiān)測的能力，因此本文以錫林浩特勝利東二號露天礦為研究對象，利用GB-InSAR技術(shù)獲取礦區(qū)邊坡的變形數(shù)據(jù)，隨后利用GIS強(qiáng)大的空間數(shù)據(jù)處理、分析能力進(jìn)行二次開發(fā)，實時、直觀地進(jìn)行露天礦邊坡變形監(jiān)測。

1 GB-InSAR與GIS集成

1.1 GB-InSAR原理

GB-InSAR采用的是正側(cè)視條帶的工作模式，系統(tǒng)采用步進(jìn)頻率連續(xù)波，同時通過控制天線的運動實現(xiàn)空間高分辨率，其距離向分辨率和方位向分辨率表達(dá)式分別為

Δr=c/2B

(1)

Δφ=λ/2L

(2)

式(1)—式(2)中，B為雷達(dá)脈沖寬度；c為光速；L為合成后的天線孔徑；λ為雷達(dá)波長。

1.2D-InSAR

D-InSAR即差分干涉測量技術(shù)，利用D-InSAR可以實現(xiàn)目標(biāo)物位移變化的精準(zhǔn)測量[8]。設(shè)GB-InSAR在不同時間獲取的同一目標(biāo)區(qū)域的兩幅SAR復(fù)圖像為I1和I2，則將復(fù)圖像的對應(yīng)像素進(jìn)行共軛相乘即可得到相位差(前提是兩幅復(fù)圖像匹配)，并形成干涉相位圖。其干涉相位計算方式為

(3)

式中，*表示共軛，且∠(I(q))=φ(q)，干涉相位圖中任意像素點q的相位差為

(4)

在實際處理過程中求得的相位差φ(q)只是(-π,π]內(nèi)的主值，若要得到真實相位差φ(q)就必須進(jìn)行相位解纏，由相位差可求得目標(biāo)的視向變形值，其計算公式為

(5)

式中，fc代表雷達(dá)的工作頻率；c為光速。

1.3 GIS

GIS以其強(qiáng)大的空間數(shù)據(jù)存儲、管理、分析等功能，已被廣泛應(yīng)用于國土、遙感、國防、資源環(huán)境等領(lǐng)域[9]。目前GIS已經(jīng)成為一項服務(wù)于信息化的技術(shù)工具，如何開發(fā)出快速、低成本的GIS應(yīng)用系統(tǒng)成為企事業(yè)單位、高校、科研院所關(guān)注的熱點。GIS已滲透到很多行業(yè)，市場上活躍著許多開發(fā)商推出的GIS二次開發(fā)產(chǎn)品，用戶可以根據(jù)實際需求進(jìn)行選擇。由Esri公司開發(fā)的GIS二次開發(fā)產(chǎn)品ArcGIS Engine成為目前的主流，它具有靈活、易用、簡潔、可移植性強(qiáng)、成本低等優(yōu)點，深受GIS從業(yè)人員的喜愛[10]。

1.4 GB-InSAR與GIS集成

通過GB-InSAR技術(shù)可以全天候、全天時地獲取露天礦邊坡的變形數(shù)據(jù)，而利用GIS可以將變形數(shù)據(jù)快速、直觀地展現(xiàn)出來。因此，GB-InSAR與GIS集成可以充分利用兩者的優(yōu)勢，實時、高效地對露天礦邊坡進(jìn)行變形監(jiān)測，集成原理如圖1所示。

圖1 GB-InSAR與GIS集成原理

如圖1所示，GB-InSAR與GIS集成是利用GB-InSAR技術(shù)獲取露天礦邊坡變形數(shù)據(jù)，以ArcSDE為存儲變形數(shù)據(jù)的數(shù)據(jù)庫與GIS之間橋梁，并利用GIS二次開發(fā)組件開發(fā)露天礦邊坡監(jiān)測系統(tǒng)，進(jìn)而實時、直觀地對露天礦邊坡進(jìn)行變形監(jiān)測。這里需要說明的是，ArcSDE(SDE即spatial database engine，空間數(shù)據(jù)庫引擎)是ArcGIS與數(shù)據(jù)庫之間通道。

2 GB-InSAR與GIS集成在露天礦邊坡變形監(jiān)測中的應(yīng)用

2.1 錫林浩特勝利露天礦概況

內(nèi)蒙古大唐國際錫林浩特礦業(yè)有限公司東二號露天礦位于勝利煤田中東部，礦區(qū)面積為49.63 km2，東西長約8.0 km，南北寬約6.6 km，主體構(gòu)造形態(tài)為NE-SW走向，兩翼不對稱的寬緩向斜。2013年初東2號露天礦已形成南北寬約2.3 km、東西長約3.5 km、面積約7.5 km2的采坑，邊坡高度202 m。隨著東二露天煤礦進(jìn)一步的開拓和延伸，南幫邊坡穩(wěn)定性成為制約其安全生產(chǎn)的重要因素，因此需要對其邊坡進(jìn)行實時的變形監(jiān)測。

2.2 GB-InSAR數(shù)據(jù)處理

GB-InSAR數(shù)據(jù)處理的原理與星載和機(jī)載SAR的相似，由于GB-InSAR采用零空間基線工作的模式，并且軌道參數(shù)已知，因此無需考慮去基線估計和地平效應(yīng)等步驟，且GB-InSAR影像獲取時間間隔短，可以忽略大氣延遲的影響[11]。GB-InSAR差分干涉測量的處理步驟主要包括GB-InSAR復(fù)圖像的配準(zhǔn)、干涉圖生成和相位噪聲濾波、相位解纏及變形值提取，GB-InSAR高精度的變形監(jiān)測需要精準(zhǔn)的干涉相位來保證，復(fù)圖像的精確配準(zhǔn)是成功進(jìn)行干涉測量的前提[12]，因此在礦區(qū)穩(wěn)固的地點安置若干個三面角反射器，并以此對GB-InSAR影像進(jìn)行配準(zhǔn)[13]。

本文結(jié)合某一天特定時段的7幅雷達(dá)監(jiān)測數(shù)據(jù)對邊坡變形提取的方法及結(jié)果進(jìn)行說明。在數(shù)據(jù)處理中利用累積干涉相位法對所有相鄰兩幅圖像相位差進(jìn)行累加，這樣可以消除相位纏繞和降低噪聲影響，設(shè)I1,I2,…IN分別為在時間t1,t2,…tN時獲得的目標(biāo)區(qū)域SAR圖像，則任意一個像元q點處目標(biāo)的累積干涉相位由式(6)可得，結(jié)合式(5)即可得q的累計位移，結(jié)果如圖2所示。

(6)

圖2 變形圖像序列

圖2為一個時間序列的6副差分后變形圖像，時長為308 min。變形方向?qū)?yīng)著坡度在視線方向的移動，負(fù)值表示距離減小，即朝著觀測者運動，從結(jié)果可以看出，礦體對應(yīng)圖中的中心區(qū)域變形較為嚴(yán)重，有些區(qū)域變形甚至呈振蕩趨勢。

2.3 GIS在露天礦邊坡變形監(jiān)測中的應(yīng)用

上文對GB-InSAR技術(shù)提取露天礦邊坡變形數(shù)據(jù)的方法進(jìn)行了詳細(xì)的闡述，并對變形數(shù)據(jù)的精度進(jìn)行了驗證。本文以ArcSDE為橋梁，將存儲、管理變形數(shù)據(jù)的Oracle數(shù)據(jù)庫與GIS有機(jī)地連接起來，并在Windows Form開發(fā)環(huán)境下，以C#為開發(fā)語言，實現(xiàn)ArcGIS Engine控件的集成，進(jìn)而完成露天礦邊坡變形監(jiān)測系統(tǒng)的開發(fā)。根據(jù)露天礦邊坡變形監(jiān)測的需求與組件式結(jié)構(gòu)和分層次開發(fā)的思想，將系統(tǒng)分為數(shù)據(jù)庫、數(shù)據(jù)層、邏輯層和表現(xiàn)層，如圖3所示。

圖3 露天礦變形監(jiān)測系統(tǒng)框架

為了盡可能滿足露天礦邊坡變形監(jiān)測的需求，本文結(jié)合GIS的功能特點，在作變形分析的同時，附加了很多其他輔助功能，如視圖管理、數(shù)據(jù)導(dǎo)出、警報設(shè)置等，系統(tǒng)功能框架如圖4所示。

從圖4可以看出，本文利用ArcSDE將數(shù)據(jù)庫與露天礦邊坡變形監(jiān)測系統(tǒng)連接起來，而GB-InSAR得到的變形數(shù)據(jù)存儲在數(shù)據(jù)庫中，從而實現(xiàn)了兩者的集成。從功能框架可以看出，該系統(tǒng)不僅包括漫游、縮放等基本功能，還可以加載、顯示礦區(qū)的實景及三維視圖，最重要的是，作為一個露天礦邊坡監(jiān)測系統(tǒng)，其可以任意設(shè)置重點監(jiān)測點和區(qū)域，并對點和區(qū)域的變形數(shù)據(jù)進(jìn)行顯示、輸出，可以根據(jù)變形信息及礦區(qū)的實際情況進(jìn)行相應(yīng)預(yù)警。

2.4 GB-InSAR與GIS集成的露天礦變形邊坡變形監(jiān)測系統(tǒng)的實現(xiàn)

GB-InSAR與GIS集成的露天礦邊坡變形監(jiān)測系統(tǒng)開發(fā)的過程繁復(fù)，且內(nèi)容較多，本文結(jié)合圖文僅對其中的視圖管理、變形分析兩大功能模塊進(jìn)行簡單的介紹，具體示例如圖5、圖6所示。

圖5(a)為雷達(dá)視圖，即為渲染后的變形數(shù)據(jù)；圖5(b)是實景視圖，其中上半?yún)^(qū)域是礦區(qū)實景圖，下半?yún)^(qū)域?qū)?yīng)的是經(jīng)過渲染的實時變形圖，兩者具有相同的坐標(biāo)系且有聯(lián)動作用；圖5(c)顯示的是監(jiān)測區(qū)域，可以清晰地看出監(jiān)測區(qū)域在三維視圖與變形圖中的位置；圖5(d)是實景圖與實時變形圖的疊加，通過該視圖可以直觀地了解露天礦邊坡整體變形情況；圖5(e)是實時變形圖與三維圖的疊加，在視圖的左側(cè)有三維視圖工具欄，可以進(jìn)行各個角度瀏覽；圖5(f)是圖5(e)的局部，圖中的深色區(qū)域為選中的監(jiān)測區(qū)域。在每個視圖的左下側(cè)是監(jiān)測區(qū)域的實時信息，主要包括數(shù)據(jù)的時間、區(qū)域信息、區(qū)域面積、位移、位移加速度、土方量等，右側(cè)為變形信息的走勢圖，如位移走勢圖、位移速度走勢圖、位移加速度走勢圖等，且可多個走勢圖同時顯示，如圖6所示，通過走勢圖可直觀地得到監(jiān)測區(qū)域的變形信息，并可以通過下方的兩個滑條自主選擇走勢圖的時間區(qū)間。

圖4 系統(tǒng)功能框架

圖5 系統(tǒng)視圖

圖6 走勢圖

3 結(jié) 論

(1) 利用GB-InSAR掃描數(shù)據(jù)，通過差分干涉測量技術(shù)獲取了露天礦邊坡變形數(shù)據(jù)。試驗結(jié)果表明，GB-InSAR技術(shù)可以快速、大范圍地獲取露天礦邊坡變形信息。

(2) 將GB-InSAR得到的變形數(shù)據(jù)通過Oracle數(shù)據(jù)庫和ArcSDE與GIS結(jié)合起來，并充分利用GIS強(qiáng)大的地圖顯示、分析功能，將變形數(shù)據(jù)通過計算、分析、渲染高效、直觀地顯示出來，為露天礦生產(chǎn)提供了有力的支持。

(3) 盡管GB-InSAR技術(shù)在地表變形監(jiān)測方面具有傳統(tǒng)技術(shù)手段難以比擬的優(yōu)勢，但是時空失相關(guān)和大氣效應(yīng)這兩項因素嚴(yán)重影響到邊坡變形信息的精度，將來有必要結(jié)合永久散射體技術(shù)對GB-InSAR技術(shù)進(jìn)行深入的研究。

[1] 姜晨光,逄曉周,董勤景,等. 露天礦邊坡自動監(jiān)測系統(tǒng)的開發(fā)與應(yīng)用[J]. 中國地質(zhì)災(zāi)害與防治學(xué)報,2006,17(4):123-128.

[2] 劉超,高井祥,王堅,等. GPS/偽衛(wèi)星技術(shù)在露天礦邊坡監(jiān)測中的應(yīng)用[J]. 煤炭學(xué)報,2010,35(5):755-759.

[3] MONSERRAT O, CROSETTO M, LUZI G. A Review of Ground-based SAR Interferometry for Deformation Measurement[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2014, 93(7): 40-47.

[4] BROUSSOLLE J, KYOVTOROV V, BASSO M, et al. MELISSA, A New Class of Ground Based InSAR System: An Example of Application in Support to the Costa Concordia Emergency[J].ISPRS Journal of Photogrammetry and Remote Sensing, 2014, 91(5):50-57.

[5] CADUFF R,KOS A,SCHLUNEGGER F,et al.Terrestrial Radar Interferometric Measurement of Hillslope Deformation and Atmospheric Disturbances in the Illgraben Debris-flow Catchment, Switzerland[J]. IEEE Geoscience and Remote Sensing Letters, 2014, 11(2): 434-438.

[6] CADUFF R, SCHLUNEGGER F, KOS A, et al. A Review of Terrestrial Radar Interferometry for Measuring Surface Change in the Geosciences[J]. Earth Surface Processes and Landforms, 2015, 40(2): 208-228.

[7] 蘇晨. 地基合成孔徑雷達(dá)變形監(jiān)測技術(shù)研究[D].重慶：重慶大學(xué),2014.

[8] 龍四春,陳鵬琦,袁英,等.GB-InSAR技術(shù)及其形變監(jiān)測[J]. 測繪通報,2015(9):1-5.

[9] 范文瑜,張榮群,高玲玲. 基于ArcGIS Engine的銀川平原濕地管理信息系統(tǒng)開發(fā)研究[J]. 測繪科學(xué),2010(S1):212-214.

[10] 吳建華. 基于ArcGIS Engine的GIS軟件開發(fā)方法[J]. 測繪通報,2010(11):54-57.

[11] 楊紅磊,彭軍還,崔洪曜. GB-InSAR監(jiān)測大型露天礦邊坡變形[J]. 地球物理學(xué)進(jìn)展,2012,27(4):1804-1811.

[12] 張祥. 地基SFCW SAR差分干涉測量技術(shù)研究[D].長沙：國防科技大學(xué),2011.

[13] 任月龍,李如仁,張信.基于多傳感器網(wǎng)的露天礦邊坡變形監(jiān)測[J]. 煤炭學(xué)報,2014,39(5):767-773.

[14] ZEBKER H A,VILLASENOR J.Decorrelation in Interferometric Radar Echoes[J]. IEEE Transactions on Geoscience and Remote Sensing, 1992, 30(5): 950-959.

[15] ZEBKER H A,ROSEN P A,HENSLEY S. Atmospheric Effects in Interferometric Aperture Radar Surface Deformation and Topographic Maps[J]. Journal of Geophysical Research(Solid Earth),1997,102 ( B4) :7547 -7563.

[16] IGLESIAS R, FABREGAS X, AGUASCA A, et al. Atmospheric Phase Screen Compensation in Ground-based SAR with a Multiple-regression Model over Mountainous Regions[J]. IEEE Transactions on Geoscience & Remote Sensing, 2013, 52(5): 2436-2449.

[17] AGRAM P. Persistent Scatterer Interferometry in Natural Terrain[D]. Stanford: Stanford University, 2010.

Slope Deformation Monitoring of Open Pit Coal Mine by GIS Integrated GB-InSAR

LI Ruren1，YANG Zhen2，YU Bo1

(1. School of Transportation Engineering, Shenyang Jianzhu University, Shenyang 110168, China;2. School of Geoscience and Surveying Engineering, China University of Mining & Technology, Beijing 100083, China)

The real-time and accurate monitoring for open pit coal mine slope can provide reliable information for slope stability and early warning, which is directly related to the economic benefit and sustained production. The ground-based InSAR (GB-InSAR), a relatively new technology for deformation monitoring, has been developed in the last ten years, which has the characteristics of all-weather and all day. However, most of the open pit slope deformation monitoring in China are in the experimental stage, which lacks the ability of continuous observation for a long time. In this paper, the basic principle of GB-InSAR is described in detail and the deformation data obtained by GB-InSAR is organically connected with GIS. Deformation information is efficiently, real-time and intuitively displayed by using of GIS powerful spatial analysis and display ability, providing powerful decision support for open-pit mine production.

GB-InSAR; open pit coal mine; GIS; deformation monitoring

李如仁，楊震，余博.GB-InSAR集成GIS的露天煤礦邊坡變形監(jiān)測[J].測繪通報，2017(5)：26-30.

10.13474/j.cnki.11-2246.2017.0147.

2016-10-08

國家重點研發(fā)計劃(2016YFC0501101)

李如仁(1966—)，男，教授，研究方向為GB-InSAR與GIS的集成。E-mail: rurenli@163.com

楊 震。E-mail: 593797596@qq.com

P237

A

0494-0911(2017)05-0026-05