0 引言
目前，對礦區沉降監測的主要手段有傳統的水準、GPS、EDM測量等。然而這些監測方法都需要耗費大量的人力，物力和財力，不適合用于對礦區進行快速、準確、實時的監測。合成孔徑干涉雷達InSAR）是近年來發展起來的一種先進的空間對地觀測技術。其方便快捷、覆蓋范圍廣、空間分辨率高等特點，使InSAR技術迅速成為國際研究的熱點。
近年來，InSAR技術在監測礦區沉陷和地質災害方面應用已取得了一定的成果，Canec等首先用D-lnSAR監測了法國Gardanne附近煤礦的沉降。基于前人研究的基礎上，本文以甘肅平涼礦區為研究對象，采用D-InSAR技術對地表開采沉陷進行監測研究，為更好的了解礦區地表形變特征，同時也將進一步拓展和完善D-InSAR技術在礦山領域的實際應用。
1 數據與處理方法
D-InSAR技術是通過對兩幅干涉條紋圖進行差分來消除地形相位的處理過程。根據去除地形相位采用的數據和處理方法的不同，D-InsAR的方法可分為：兩軌法、三軌法和四軌法[8]。文中采用兩軌法進行處理。
本文以甘肅華亭煤礦為例，開展利用D-InSAR技術提取煤礦沉陷形變的量級和范圍的研究。華亭煤礦位于甘肅省平涼市，從地質構造角度，該區域位于六盤山斷裂帶的中段東側�？紤]到研究區域地表植被覆蓋茂密，本文選取了日本ALOS PALSAR衛星影像資料，波長23.6cm，對植被有較好的穿透能力，影像獲取時間為2007年7月至2010年7月，共11景升軌數據。
本文采用瑞士GAMMA公司開發的GAMMA雷達干涉處理軟件進行SAR影像數據處理，具體流程如下[9]：
①獲取干涉圖：確定主、輔影像，經配準、干涉、基線估計等步驟后得到干涉相位圖；
�、谖闹械腄EM采用美國NASA公開的SRTM DEM來模擬地形相位，空間分辨率為90m（3弧秒）。將SRTM DEM與主影像進行配準，轉換到雷達坐標系下，模擬地形相位。
�、郢@取形變量：從步驟①中獲取的干涉圖去除地形相位，得到差分干涉圖，經相位解纏，地理編碼，獲得形變信息。
數據處理時，為了提高干涉圖的信噪比，對干涉圖進行了兩次基于頻譜的自適應濾波，濾波窗口分別為32×32、64×64，多視因子給定4：9，即每個像元約為60×60m2。通過自由組合干涉，共挑選出7個相干性較好的干涉圖，各干涉圖的空間垂直基線和時間間隔如圖1所示，顯然，大部分干涉圖垂直基線較短，小于1km，較大程度地削弱了外部DEM誤差對最終形變結果的影響。
2 結果分析
華亭煤礦是上個世紀甘肅省發現并開發規模最大的煤田，目前開采量依舊較大，通過對7幅干涉圖（圖2）的`對比分析可以得到如下幾點結論和認識：
�、貺波段具有較強的植被穿透能力，干涉圖相干性較好，僅在煤礦沉陷形變梯度較大的地方失相干，其他區域全部相干，顯示出L波段雷達影像在煤礦沉陷監測中應用的潛力和優勢；
�、诿旱V沉陷的時間應發生在2008年5月26日到2009年7月14日之間。首先，圖2b、c、d均選取的是非公共影像干涉圖，表明干涉圖中的條紋信號不是大氣相位延遲，也不是外部DEM誤差引起的相位，即2007年至2009年煤礦發生了沉陷，引起了地表局部范圍變形；其次，圖2a、e、g中并未顯示出明顯的形變條紋，由此可以斷定煤礦沉陷發生的時間介于2008年5月26日到2009年7月14日之間；
�、蹚膱D2b、c、d中可以看出，煤礦沉陷引起的地表最大變形約8cm；
④通過圖2b、c、d，可以圈定煤礦沉陷的形狀和范圍。沉陷區域近似南北向的矩形，其中南北長約3km，東西寬約2km。
一般煤礦沉降可以簡單的分為三個階段：緩慢沉降階段、快速沉降階段和沉降減緩并且消失階段。圖2a、e應屬緩慢沉降階段，圖2b、c、d、f應為快速沉降階段，圖g為沉降消失階段。通過監測沉降速率緩慢階段的出現并采取積極措施能夠避免危害性較大的快速沉降階段的出現，因此監測地面沉降對保證礦區安全生產和及時控制并阻止地面沉降的發生非常有用。
3 結論
本文采用兩軌法D-InSAR技術獲取了甘肅平涼華亭煤礦在2007-2010年的地表形變場，結果表明，華亭煤礦沉陷在2008年5月26日到2009年7月14日之間十分明顯，最大形變達到8cm，沉陷范圍為近似南北向矩形，南北長約3km，東西寬約2km。最后通過區域小震活動分布情況證實了本文InSAR結果的可靠性。表明D-InSAR技術監測到的沉陷信息可以較為準確的反映礦區實際沉陷狀況，顯示出了巨大的潛力和優勢，可以作為一種獲取大范圍煤礦區的地表沉陷的有效方法，為研究我國煤礦區的沉陷提供了有效的方法。
參考文獻：
.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[J].Joumal of Applied Geophysies，2000，43：43-54.
. Poland Engineering Geology，2009，106（3）：58-67.
[4]吳立新，高均海，葛大慶等.工礦區地表沉陷D-InSAR監測試驗研究[J].東北大學學報（自然科學版），2005，26（8）：775-752.
[5]董玉森，Ge Liulin， Chang Hsingchun，張志.基于差分雷達干涉測量的礦區地面沉降監測研究[J].武漢大學學報（信息科學版），2007，32（l0）：888-892.
[6]劉廣，郭華東， Ramon Hanssen，等. InSAR技術在礦區沉降監測中的應用研究[J]，國土資源遙感，2008，76（2）：51-57.
[7]范洪冬.InSAR若干關鍵算法及其在地表沉降監測中的應用研究[D].中國礦業大學，2010.
[8]Massonnet D， Rossi M， Carnlona C，etal. The displacement field of the Landers earthquake mapped by radar interferometry [J]. Nature， 1993， 364（8）：138-142.
[9]陳艷玲.星載SAR及InSAR技術在地球科學中的應用研究[D].中國科學院研究生院（上海天文臺），2007.
【基于DInSAR技術的礦區發展的研究綜述論文】相關文章：
基于創業投資發展對策建議綜述論文12-10
綠色建材的研究與發展綜述分析論文05-21
基于知識和組織學習理論的技術聯盟形成研究綜述的論文04-12
影響家庭農場發展的因素研究綜述論文05-23
網絡道德研究綜述的論文04-12
姜夔研究綜述的論文04-12