基于D—InSAR技術(shù)的礦區(qū)發(fā)展的研究綜述

1、基于DInSAR技術(shù)的礦區(qū)發(fā)展的研究綜述0 引言目前，對(duì)礦區(qū)沉降監(jiān)測(cè)的主要手段有傳統(tǒng)的水準(zhǔn)、GPS、EDM測(cè)量等。然而這些監(jiān)測(cè)方法都需要耗費(fèi)大量的人力，物力和財(cái)力，不適合用于對(duì)礦區(qū)進(jìn)行快速、準(zhǔn)確、實(shí)時(shí)的監(jiān)測(cè)。合成孔徑干涉雷達(dá)InSAR）是近年來發(fā)展起來的一種先進(jìn)的空間對(duì)地觀測(cè)技術(shù)。其方便快捷、覆蓋范圍廣、空間分辨率高等特點(diǎn)，使InSAR技術(shù)迅速成為國(guó)際研究的熱點(diǎn)。近年來，InSAR技術(shù)在監(jiān)測(cè)礦區(qū)沉陷和地質(zhì)災(zāi)害方面應(yīng)用已取得了一定的成果，Canec等首先用D-lnSAR監(jiān)測(cè)了法國(guó)Gardanne附近煤礦的沉降【1】。2000年，Canec和Delaeourt又用相同的方法得到了該煤礦1992到1

2、995年的地面沉陷變化，和地面的水準(zhǔn)測(cè)量結(jié)果相似【2】；2009年，Perski等對(duì)波蘭的wieliezka鹽礦，采用InsAR技術(shù)進(jìn)行監(jiān)測(cè)與分析【3】。目前，我國(guó)學(xué)者也己開展了用D-InSAR技術(shù)監(jiān)測(cè)礦區(qū)沉降的研究?；谇叭搜芯康幕A(chǔ)上，本文以甘肅平?jīng)龅V區(qū)為研究對(duì)象，采用D-InSAR技術(shù)對(duì)地表開采沉陷進(jìn)行監(jiān)測(cè)研究，為更好的了解礦區(qū)地表形變特征，同時(shí)也將進(jìn)一步拓展和完善D-InSAR技術(shù)在礦山領(lǐng)域的實(shí)際應(yīng)用。1 數(shù)據(jù)與處理方法D-InSAR技術(shù)是通過對(duì)兩幅干涉條紋圖進(jìn)行差分來消除地形相位的處理過程。根據(jù)去除地形相位采用的數(shù)據(jù)和處理方法的不同，D-InsAR的方法可分為：兩軌法、三軌法和四軌法

3、。文中采用兩軌法進(jìn)行處理。本文以甘肅華亭煤礦為例，開展利用D-InSAR技術(shù)提取煤礦沉陷形變的量級(jí)和范圍的研究。華亭煤礦位于甘肅省平?jīng)鍪校瑥牡刭|(zhì)構(gòu)造角度，該區(qū)域位于六盤山斷裂帶的中段東側(cè)。考慮到研究區(qū)域地表植被覆蓋茂密，本文選取了日本ALOS PALSAR衛(wèi)星影像資料，波長(zhǎng)23.6cm，對(duì)植被有較好的穿透能力，影像獲取時(shí)間為2007年7月至2010年7月，共11景升軌數(shù)據(jù)。本文采用瑞士GAMMA公司開發(fā)的GAMMA雷達(dá)干涉處理軟件進(jìn)行SAR影像數(shù)據(jù)處理，具體流程如下：獲取干涉圖：確定主、輔影像，經(jīng)配準(zhǔn)、干涉、基線估計(jì)等步驟后得到干涉相位圖；文中的DEM采用美國(guó)NASA公開的SRTM DEM來模

4、擬地形相位，空間分辨率為90m（3弧秒）。將SRTM DEM與主影像進(jìn)行配準(zhǔn)，轉(zhuǎn)換到雷達(dá)坐標(biāo)系下，模擬地形相位。獲取形變量：從步驟中獲取的干涉圖去除地形相位，得到差分干涉圖，經(jīng)相位解纏，地理編碼，獲得形變信息。數(shù)據(jù)處理時(shí)，為了提高干涉圖的信噪比，對(duì)干涉圖進(jìn)行了兩次基于頻譜的自適應(yīng)濾波，濾波窗口分別為32x32、64x64，多視因子給定4：9，即每個(gè)像元約為60x60m2。通過自由組合干涉，共挑選出7個(gè)相干性較好的干涉圖，各干涉圖的空間垂直基線和時(shí)間間隔如圖1所示，顯然，大部分干涉圖垂直基線較短，小于1km，較大程度地削弱了外部DEM誤差對(duì)最終形變結(jié)果的影響。2 結(jié)果分析華亭煤礦是上個(gè)世紀(jì)甘肅省

5、發(fā)現(xiàn)并開發(fā)規(guī)模最大的煤田，目前開采量依舊較大，通過對(duì)7幅干涉圖（圖2）的對(duì)比分析可以得到如下幾點(diǎn)結(jié)論和認(rèn)識(shí)：L波段具有較強(qiáng)的植被穿透能力，干涉圖相干性較好，僅在煤礦沉陷形變梯度較大的地方失相干，其他區(qū)域全部相干，顯示出L波段雷達(dá)影像在煤礦沉陷監(jiān)測(cè)中應(yīng)用的潛力和優(yōu)勢(shì)；煤礦沉陷的時(shí)間應(yīng)發(fā)生在2008年5月26日到2009年7月14日之間。首先，圖2b、c、d均選取的是非公共影像干涉圖，表明干涉圖中的條紋信號(hào)不是大氣相位延遲，也不是外部DEM誤差引起的相位，即2007年至2009年煤礦發(fā)生了沉陷，引起了地表局部范圍變形；其次，圖2a、e、g中并未顯示出明顯的形變條紋，由此可以斷定煤礦沉陷發(fā)生的時(shí)間介

6、于2008年5月26日到2009年7月14日之間；從圖2b、c、d中可以看出，煤礦沉陷引起的地表最大變形約8cm；通過圖2b、c、d，可以圈定煤礦沉陷的形狀和范圍。沉陷區(qū)域近似南北向的矩形，其中南北長(zhǎng)約3km，東西寬約2km。一般煤礦沉降可以簡(jiǎn)單的分為三個(gè)階段：緩慢沉降階段、快速沉降階段和沉降減緩并且消失階段。圖2a、e應(yīng)屬緩慢沉降階段，圖2b、c、d、f應(yīng)為快速沉降階段，圖g為沉降消失階段。通過監(jiān)測(cè)沉降速率緩慢階段的出現(xiàn)并采取積極措施能夠避免危害性較大的快速沉降階段的出現(xiàn)，因此監(jiān)測(cè)地面沉降對(duì)保證礦區(qū)安全生產(chǎn)和及時(shí)控制并阻止地面沉降的發(fā)生非常有用。3 結(jié)論本文采用兩軌法D-InSAR技術(shù)獲取了

7、甘肅平?jīng)鋈A亭煤礦在2007-2010年的地表形變場(chǎng)，結(jié)果表明，華亭煤礦沉陷在2008年5月26日到2009年7月14日之間十分明顯，最大形變達(dá)到8cm，沉陷范圍為近似南北向矩形，南北長(zhǎng)約3km，東西寬約2km。最后通過區(qū)域小震活動(dòng)分布情況證實(shí)了本文InSAR結(jié)果的可靠性。表明D-InSAR技術(shù)監(jiān)測(cè)到的沉陷信息可以較為準(zhǔn)確的反映礦區(qū)實(shí)際沉陷狀況，顯示出了巨大的潛力和優(yōu)勢(shì)，可以作為一種獲取大范圍煤礦區(qū)的地表沉陷的有效方法，為研究我國(guó)煤礦區(qū)的沉陷提供了有效的方法。參考文獻(xiàn)：【1】Carnec C， Massonnet D， and King C. Two examples of the use of

8、 SAR interferometry on displacement fields of small spatial extent .Geophysical Research Letters， 1996， 23（24）：3579-3582.【2】Carnec C and Delaeourt C. Three years of mining subsidence monitored by SAR inierferometry near Gardanne Franee.Joumal of Applied Geophysies，2000，43：43-54.【3】Perski Z，Hanssen R

9、，Wojeik A， etal. InSAR analyses of terrain deformation near the Wieliczka Salt Mine. Poland Engineering Geology，2009，106（3）：58-67.【4】吳立新，高均海，葛大慶等.工礦區(qū)地表沉陷D-InSAR監(jiān)測(cè)試驗(yàn)研究.東北大學(xué)學(xué)報(bào)（自然科學(xué)版），2005，26（8）：775-752.【5】董玉森，Ge Liulin， Chang Hsingchun，張志.基于差分雷達(dá)干涉測(cè)量的礦區(qū)地面沉降監(jiān)測(cè)研究.武漢大學(xué)學(xué)報(bào)（信息科學(xué)版），2007，32（l0）：888-892.【6】劉廣，郭華東， Ramon Hanssen，等. InSAR技術(shù)在礦區(qū)沉降監(jiān)測(cè)中的應(yīng)用研究，國(guó)土資源遙感，2008，76（2）：51-57.【7】范洪冬.InSAR若干關(guān)鍵算法及其在地表沉降監(jiān)測(cè)中的應(yīng)用研究.中國(guó)礦業(yè)大學(xué)，2010.Massonnet D， Rossi M， Carnlona C，etal. The displace