宋繼德，邸志眾

(1．陜西省地質(zhì)礦產(chǎn)勘查開發(fā)局 第二綜合物探大隊，陜西 西安 710016;2．北京中勘邁普科技有限公司，北京 100086)

0 引言

InSAR(SAR Interferometry 或Interferometric SAR)，即合成孔徑雷達干涉測量，是指利用SAR 數(shù)據(jù)中的相位信息進行干涉測量處理，結(jié)合雷達參數(shù)和衛(wèi)星位置信息反演地表三維及其變化信息的遙感手段［1，2］。D-InSAR 即差分合成孔徑雷達干涉測量，是指對干涉相位進行差分處理提取變化信息的干涉測量手段，是InSAR 技術(shù)的延伸。InSAR 技術(shù)已被廣泛用于大區(qū)域地形測圖和地表信息提取。常見的InSAR 方式有順軌(Along track)，交軌(Across track)和重復(fù)軌(Repeat pass)。由于地表形變監(jiān)測是測量一定時間間隔內(nèi)的地表形變量，因而重復(fù)軌是目前應(yīng)用的主要方式，具體是指利用同一地區(qū)不同時刻獲取的SAR 數(shù)據(jù)進行差分干涉處理以提取地表形變信息。由于InSAR具有全天候、全天時、覆蓋面廣和高精度獲取地表信息的能力，因而在地球科學(xué)的諸多領(lǐng)域如地震、火山、礦區(qū)塌陷、地面沉降、滑坡、冰川等災(zāi)害中逐漸得到了廣泛的應(yīng)用［1，5］。

本文以山東省重點地區(qū)礦山開采的遙感動態(tài)監(jiān)測為研究內(nèi)容，展開濟寧、兗州地區(qū)地面沉降InSAR 監(jiān)測試驗研究，圍繞濟寧和兗州地區(qū)煤礦開采、地下水開采等多種因素作用下地面沉降(陷)的分布與發(fā)展?fàn)顩r，介紹以D-InSAR 技術(shù)為主要手段的地面沉降調(diào)查與監(jiān)測工作的研究情況與成果。

1 差分合成孔徑雷達干涉測量原理與數(shù)據(jù)處理方法

1．1 基本原理

重復(fù)軌干涉測量的幾何原理［2］，如圖1 所示。圖中O1、O2分別為兩次獲取SAR 信號的天線位置，B 為兩天線間空間距離，稱為空間基線。

圖1 重復(fù)軌D-InSAR 幾何示意圖Fig．1 Geometric schematic diagram of repeat pass D-InSAR

設(shè)基線與水平方向的夾角為α，H 為平臺高度，地面點P到兩天線間的距離分別為R1和R2，θ 為第一副天線的側(cè)視角，P 點高程為h。將基線沿雷達視線向分解，垂直于視線向分量為B⊥，稱為垂直基線。雷達對地物觀測時，接收到地物的后向散射回波信號可以表示為:

兩次觀測所獲取的SAR 影像分別稱為主影像(Master)和輔影像(Slave)，接收到P 點的雷達信號分別為:

雷達后向散射回波信號中的相位信息包含兩部分信息，一是與往返路徑相關(guān)的相位，另一部分是目標(biāo)后向散射相位，可以表示為:

兩次獲取的相位分別為:

式中:ΔR 為路程差，與雷達與目標(biāo)的距離相關(guān)。若兩次雷達獲取的相位中的散射相位保持穩(wěn)定，即:

則經(jīng)過共軛相乘處理后得到的干涉相位為:

干涉相位φint的貢獻值(Contributions)由以下5 個因素構(gòu)成:平地相位δφflat;地形相位δφtopo;兩次觀測期間觀測目標(biāo)沿雷達視線向(LOS)移動引起的相位變化δφmov;電磁波傳輸過程中大氣波動引起的相位延遲δφatm;其他噪聲相位δφnoise。則干涉相位可以表示為:

1．2 數(shù)據(jù)處理方法

根據(jù)式(8)，要得到反映地形信息的干涉相位δφtopo，應(yīng)當(dāng)首先從干涉相位和中消除平地相位的影響，其次從中減去其他各項得到地形相位，再進行相位解纏和相位到高程的轉(zhuǎn)換，將生成的高程信息轉(zhuǎn)換到特定的坐標(biāo)系統(tǒng)下，生成DEM 數(shù)據(jù)［1，6］。由于地形形變相位、熱噪聲相位、大氣波動相位以及軌道偏差引起的相位波動可以通過合理的處理方法來消除或者降低。對于ERS-1/2，其特殊的“串行模式”(Tandem Mode)SAR 數(shù)據(jù)是生成DEM 的最佳選擇［6］。因而，利用InSAR 獲取DEM 數(shù)據(jù)的處理流程，見圖2。

圖2 InSAR 數(shù)據(jù)處理流程Fig．2 Data processing flow of InSAR

1．3 D-InSAR 地表形變監(jiān)測技術(shù)及應(yīng)用

隨著D-InSAR 技術(shù)研究的不斷深入，利用其進行地表形變信息獲取的數(shù)據(jù)處理方法也在不斷進步。針對不同形變體的變形特征，產(chǎn)生了不同的監(jiān)測方法。按照監(jiān)測目標(biāo)空間分布的不同和數(shù)據(jù)處理方式的差別，可將InSAR 信息提取方式分為兩類:基于點目標(biāo)的D-InSAR 方法和基于連續(xù)目標(biāo)的D-InSAR 監(jiān)測方法。

基于單個點目標(biāo)的方法以角反射器方法為代表，主要針對的是形變區(qū)域無相干目標(biāo)或存在極少數(shù)相干目標(biāo)的情況，監(jiān)測區(qū)域較小，適合于局部變形監(jiān)測，如滑坡體、危巖活動監(jiān)測等［1，7，8］。

連續(xù)目標(biāo)方法指的是在空間大面積范圍內(nèi)，長時間條件下依然存在較多相干目標(biāo)的情況［1］。常規(guī)差分干涉測量方法難以解決多時相動態(tài)監(jiān)測問題。緩慢連續(xù)變形體的移動監(jiān)測依賴于長時間連續(xù)監(jiān)測，對于重復(fù)軌監(jiān)測，長時間周期會引起失相干，這種情況下，對于大區(qū)域連續(xù)形變體的移動監(jiān)測則比較困難。由于監(jiān)測對象多是與人們生產(chǎn)和生活緊密相關(guān)的城市或城鎮(zhèn)，這些區(qū)域在長時間間隔下依然保持了很好的相干性，因而，可以通過多時相分析，獲取這些區(qū)域的地表形變信息。由此，則產(chǎn)生了累積紋圖法(Stacking Interferograms)方法和永久性散射體(Permanent Scatterers)技術(shù)。這兩種方法通過長時間序列分析，彌補了常規(guī)重復(fù)軌干涉測量的不足，提高了時間采樣率，能更準(zhǔn)確的獲取形變隨時間變化而產(chǎn)生的形變?？梢詫⑿巫兲卣鞣譃榫€性和非線性分別進行處理［3，4，8］。但由于永久性散射體技術(shù)對SAR 數(shù)據(jù)數(shù)目要求較高，一般情況下難以完成大面積形變。

2 研究區(qū)域地表形變InSAR 監(jiān)測的技術(shù)方法

山東省兗州地區(qū)為典型的煤礦開采區(qū)，開采沉陷誘發(fā)的礦區(qū)環(huán)境破壞不僅出現(xiàn)在井田開采上覆地區(qū)，同時也影響到了城市和城鎮(zhèn)的地面穩(wěn)定性，因而礦山沉陷InSAR 監(jiān)測的首要任務(wù)是明確2004 ～2010 年開采沉陷的變化特征，進而選擇合理的數(shù)據(jù)處理方法和相關(guān)技術(shù)參數(shù)。

針對緩慢變形和快速沉降兩種不同的形變類型，采用常規(guī)D-InSAR 技術(shù)與永久散射體干涉測量相結(jié)合的工作方法、利用永久性散射體干涉測量(PSI)技術(shù)相較其他技術(shù)，在剔除地形誤差和大氣波動對相位信號的影響方面的優(yōu)勢，且目前可獲得的SAR 數(shù)據(jù)在精度和分辨率方面也足以滿足PSI 對數(shù)據(jù)的要求。

2．1 永久散射體干涉測量(PSInSAR)技術(shù)

永久散射體干涉測量(Permanent Scatterers Interferometry SAR，PSInSAR)的核心思想是對永久散射體干涉相位進行時間序列分析，根據(jù)各相位分量的時空特征，估算大氣波動，DEM 誤差以及噪聲等，將其從差分干涉相位中逐個分離，最終獲取每個PS 的線性和非線形形變速率、大氣延遲量(Atmosphere Phase Screen)以及DEM 誤差。經(jīng)PS 方法處理，獲取的年度形變率的精度可以達到毫米級。其基本思想是基于大量的SAR 數(shù)據(jù)(20-30 景甚至更多)，從中篩選出具有穩(wěn)定散射特性的相干點目標(biāo)，構(gòu)成離散點觀測網(wǎng)絡(luò)(較之常規(guī)的變形監(jiān)測網(wǎng)密度更高)，通過分析PS 點目標(biāo)相位變化獲取地表形變狀況。由于將永久散射體作為觀測對象，降低了空間基線對相干性的影響，即使在臨界基線的條件下，仍然可以通過分析PS 差分干涉相位的變化反演形變信息、提取地表形變速率和累積形變量［8］(見圖3)。

PS 技術(shù)一般采用線性形變模型提取點目標(biāo)對應(yīng)的形變量，非線性模擬的處理過程復(fù)雜，而且非常耗時，限制了其用于大面積的形變測量，但隨著處理技術(shù)的進步，處理時間將逐步縮小，處理的范圍也可再進一步擴大。

PSInSAR 對數(shù)據(jù)量要求較高，只有SAR 影像個數(shù)達到一定的程度才能篩選出在整個時間跨度內(nèi)具有穩(wěn)定信號的PS 點，對于相干目標(biāo)較多的地區(qū)，如城區(qū)等，由于地物在長時間間隔下保持了較高的相干性，數(shù)據(jù)量較少的情況下也可進行處理。一般情況下，至少必須滿足5 個/km2才能完成相位解纏。

圖3 PSInSAR 數(shù)據(jù)處理的基本方法Fig．3 The basic method of data processing of PSInSAR

PS 方法在處理過程中通過分析每個點目標(biāo)對應(yīng)的相位，將其中包含的因地表形變，高程以及大氣波動引起的相位分量加以區(qū)分，最終得到各個點對應(yīng)的形變相位。通過分析PS 數(shù)據(jù)集，估算每個干涉紋圖中由大氣波動引起的相位分量，此外可以確定出每一點精確的高程信息。完成上述過程后，在干涉處理時將二者對相位的貢獻量加以去除。經(jīng)PS 方法處理，獲取的年度形變率的精度可以達到±0．1 mm/a。

2．2 可獲取的主要雷達數(shù)據(jù)及參數(shù)分析

當(dāng)前用于InSAR 技術(shù)研究和應(yīng)用較多的是ERS -1/2、ENVISAT、RADARSAT -1/2 以及ALOS PALSAR 獲取的SAR 數(shù)據(jù)，其空間分辨多為20 ～30 m，可歸為中等分辨率雷達數(shù)據(jù)。與之相比較而言，近幾年發(fā)射的雷達衛(wèi)星以高分辨率為主，包括Cosmo- Skymed 及TerraSAR - X，其最高分辨率可以達到1 ～3 m。表1 所列為當(dāng)前幾種雷達傳感器所提供的SAR 數(shù)據(jù)。

表1 SAR 數(shù)據(jù)及相關(guān)參數(shù)Tab．1 SAR data and relevant parameters

從雷達波段而言，目前主要的中等分辨率數(shù)據(jù)為C 波段數(shù)據(jù)(ENVISAT)。由于這些衛(wèi)星具有較高的定位精度(可以達到5 cm)，空間基線解算結(jié)果可靠，數(shù)據(jù)相干性較好，因此，在多數(shù)地表形變監(jiān)測中均可達到較好的監(jiān)測結(jié)果。高分辨率雷達衛(wèi)星TerraSAR-X 重復(fù)周期為11 d，采用X 波段，對地表形變等變化信息更為敏感，可以應(yīng)用于形變速率極快的監(jiān)測，礦山開采沉陷等非常適用。Cosmo-Skymed 與TerraSAR -X 相似，也屬于高分辨率雷達數(shù)據(jù)，最高分辨率達1 m，掃描模式獲取數(shù)據(jù)為3 m地面分辨率。在空間基線不大的條件下，L 波段的SAR 數(shù)據(jù)能夠在很長的時間內(nèi)保持較高的相干性，能滿足變化幅度較大的形變監(jiān)測。特別對于快速的地表移動監(jiān)測，應(yīng)用L 波段數(shù)據(jù)也較為合適。同時，由于這些數(shù)據(jù)分辨率的提高，空間采樣單元的個數(shù)將會增大，能更為準(zhǔn)確的描述地表形狀。因此，在很大程度上滿足了從短時間，變化速率高以及長時間，變化幅度大的形變監(jiān)測。

2．3 工作方法與技術(shù)

針對兗州地區(qū)地表形變調(diào)查和監(jiān)測的工作重點，主要目標(biāo)為查明工作區(qū)地面沉降(陷)的歷史狀況和現(xiàn)狀，調(diào)查兗州、濟寧和鄒城下屬的多個礦區(qū)開采沉陷的基本特征和分布范圍，開展兗州市地面沉降InSAR 監(jiān)測。因而，采用常規(guī)D-InSAR 技術(shù)與永久散射體干涉測量相結(jié)合的工作方法進行。

2．3．1 工作區(qū)SAR 數(shù)據(jù)獲取

采用的雷達數(shù)據(jù)為歐空局自2004 年11 月至2010 年10 月間接收的ENVISAT SAR 數(shù)據(jù)，其有效工作波段是C 波段，數(shù)據(jù)獲取周期為35 d，空間分辨率為20 m，是目前常用且應(yīng)用效果顯著的SAR 數(shù)據(jù)。獲取ASAR 數(shù)據(jù)共計24 景，如表2 所示。覆蓋范圍，見圖4。

表2 ENVISAT 衛(wèi)星數(shù)據(jù)列表Tab．2 List of ENVISAT satellite data

圖4 ENVISAT 衛(wèi)星數(shù)據(jù)覆蓋范圍Fig．4 Coverage of ENVISAT satellite data

區(qū)域性地面沉降監(jiān)測實際使用的數(shù)據(jù)從2004 年12 月到2010 年10 月計20 景，根據(jù)區(qū)域地面沉降長時序監(jiān)測需求形成干涉像對。

2．3．2 數(shù)據(jù)處理和分析

數(shù)據(jù)處理和分析系統(tǒng)以專業(yè)InSAR 數(shù)據(jù)處理軟件，如GAMMA、Doris 和SARSCAPE 為主，選擇利用ArcGIS 平臺為數(shù)據(jù)后處理系統(tǒng)，針對工作區(qū)環(huán)境特點及礦區(qū)開采沉陷關(guān)鍵In-SAR 處理分析技術(shù)，分別利用常規(guī)D-InSAR 技術(shù)和永久散射體InSAR 技術(shù)分別針對快速沉陷和緩慢沉降開展監(jiān)測。InSAR 數(shù)據(jù)處理過程主要包括最優(yōu)干涉數(shù)據(jù)選擇、干涉像對的基線估算、平地效應(yīng)去除、相干性分析、大氣效應(yīng)分析、數(shù)據(jù)配準(zhǔn)、離散點相位解纏，時間序列分析方法等。

2．3．3 地表形變信息提取

采用永久散射體InSAR 時間序列分析技術(shù)提取地表形變信息，并結(jié)合若干實地觀測結(jié)果獲取不同時間間隔的地面沉降速率圖。針對礦山開采沉陷的發(fā)生特征，利用常規(guī)D-InSAR 技術(shù)開展數(shù)據(jù)處理，城市地表緩慢變化和活動斷裂移動監(jiān)測利用PSInSAR 技術(shù)開展數(shù)據(jù)處理與分析。

2．3．4 監(jiān)測結(jié)果后處理

在多期地面沉降信息分析與綜合的基礎(chǔ)上編制工作區(qū)地面沉降(陷)速率圖和累積沉降量圖。

3 區(qū)域地面沉降與開采沉陷InSAR 監(jiān)測與分析

3．1 煤礦區(qū)開采沉陷歷史狀況InSAR 動態(tài)監(jiān)測

考慮到礦山塌陷的大變形過程對相干性的影響，以及季節(jié)性植被覆蓋變化，選擇冬季干涉像對，且間隔時間較短的數(shù)據(jù)進行D-InSAR 數(shù)據(jù)處理，共獲取8 個能清晰反應(yīng)礦山沉陷的干涉紋圖，用于反映工作區(qū)礦區(qū)開采沉陷2011 年之前歷史狀況的動態(tài)過程。

全區(qū)總的開采沉陷區(qū)主要分布于兗州、曲阜和鄒城沿市泗河沿岸以及濟寧市區(qū)周邊。這些地區(qū)不同時段內(nèi)差分干涉條紋密集，表明了各個煤礦在監(jiān)測時段內(nèi)開采沉陷的發(fā)生位置和范圍。在2004 年11 月到2005 年2 月3 個月的時間內(nèi)，共有4 個差分干涉紋圖(圖5a、b、c、d)顯示了該時間段的塌陷變化過程，各個采區(qū)的干涉條紋密度隨著時間間隔的增加而整體上增大。由041121 -050306 差分干涉紋圖(圖5b)可知兗州市城區(qū)東部與曲阜交界(序號2)的開采沉陷漏斗干涉條紋達5 個色周，沉陷量達到15 cm，表明在105 d內(nèi)由于煤礦開采引起了最大15 cm的沉降，沉降呈現(xiàn)出盆地形狀的分布特征，與傳統(tǒng)開采沉陷的沉陷盆地分布相同。位于濟寧南張鎮(zhèn)西北(序號1)的沉陷漏斗在041121 -050130(圖5a)的干涉紋圖超過2 個色周，沉陷量超過6 cm。050103 和050306 兩景SAR 數(shù)據(jù)相干性好，差分干涉紋圖(圖5d)各沉陷中心清晰，可以明顯看到該沉陷漏斗中心向北偏移。

2005 年12 月到2006 年2 月有1 個相干性好的干涉紋圖(圖5e)，該干涉紋圖表現(xiàn)了該時間段內(nèi)各沉陷中心的沉陷狀況，曲阜與兗州在泗河交界處的沉陷漏斗干涉條紋有5 -6 個色周，中心沉陷量達15 ～18 cm;濟寧南張鎮(zhèn)-劉堤頭間(序號3)的沉陷漏斗條紋達4 個色周，中心沉陷量超過12 cm。2009 年12 月到2010 年10 月相干性良好的3 個干涉紋圖(圖5f、g、h)清晰地表現(xiàn)了該時間段內(nèi)的開采沉陷狀況。由091220 -100124(圖5f)和100124 -100228(圖5g)干涉紋圖可知，位于曲阜與兗州在泗河交界處張家村煤礦一帶(序號4)的沉陷范圍在35 d的時間間隔內(nèi)達到3 個色周的沉降變化，表明在091220 -100228共70 d的時間內(nèi)沉陷量達到近18 cm，屬于快速沉降階段。濟寧主城區(qū)西北前王村南沉陷中心(序號5)在091220 -100124 和100124 -100228 干涉紋圖中分別達到了6 個和3 個以上色周，共計約10 個條紋，表明最大沉陷量超過30 cm。圖5h 為直接利用091220 和100228 兩景數(shù)據(jù)干涉后的條紋圖，反映了105 d內(nèi)整個沉降區(qū)的變化狀況，是圖5f 和圖5g 反映的沉降量的總和。

總體上，干涉圖的空間位置和分布形態(tài)一致，但局部干涉圖的條紋清晰程度降低。這種結(jié)果產(chǎn)生的原因是由于時間間隔進一步增大，相干性下降引起的。同時由于沉陷量增大，局部出現(xiàn)非相干的移動變化，這種快速變化限制了利用InSAR 干涉圖之間確定沉降量的精準(zhǔn)程度。同時，位于樓莊和梁家營一帶出現(xiàn)了局部回彈，干涉條紋與其它沉降引起的形態(tài)相反，變化量為1個色周，影響范圍約2 ～3 km。

3．2 全區(qū)2004 ～2010 年區(qū)域性地面沉降InSAR 連續(xù)監(jiān)測

利用長時間序列進行PSInSAR 處理，提取相干目標(biāo)的沉降速率。研究區(qū)內(nèi)的濟寧、兗州等城市，地面沉降主要分為緩慢沉降以及煤礦開采沉陷。前者為緩慢微小變化，主要發(fā)生于濟寧市城區(qū)以及城區(qū)周邊，而后者則多發(fā)生于煤礦區(qū)，與采礦條件、煤層厚度與賦存條件，開采方式等因素相關(guān)。因而，在煤礦區(qū)開采沉陷動態(tài)監(jiān)測的基礎(chǔ)上，實現(xiàn)全區(qū)地面沉降的連續(xù)監(jiān)測。

圖5 工作區(qū)煤礦開采沉陷InSAR 動態(tài)監(jiān)測結(jié)果圖Fig．5 Dynamic monitoring results for applying InSAR in subsidence area caused by coal mining

利用2004 年12 月～2006 年02 月和2008 年9 月～2010 年10 月兩個時間段的ENVISAT ASAR 數(shù)據(jù)進行PSInSAR 時序分析，得到兩個時間段的平均沉降速率，如圖6、圖7 所示。監(jiān)測結(jié)果表明:工作區(qū)總的沉降范圍與礦區(qū)分布一致，分別位于濟寧市任城區(qū)北部、兗州市東南部、曲阜市西部以及鄒城市西北部;工作區(qū)內(nèi)有兩大明顯的沉降區(qū)域，兗州市東南、曲阜市西部以及鄒城市西北部沉降連接成片，沉降速率大于20 mm/a的區(qū)域超過180 km2，另一個明顯的沉降區(qū)域為濟寧市任城區(qū)北部。

2008 年09 月～2010 年10 月InSAR 監(jiān)測地面沉降結(jié)果(見圖7)表明兩大沉降區(qū)沉降范圍有所擴大，兗州—曲阜—鄒城沉降區(qū)內(nèi)沉降中心沉降速率有所減緩，但西部整體沉降加劇，東部有所減緩;濟寧任城區(qū)沉降區(qū)與汶上縣、嘉祥縣交界處沉降明顯減緩，主城區(qū)西北郊南張鎮(zhèn)一帶有所加劇，北部岱莊—栗鄉(xiāng)一帶有所減緩，但任城區(qū)東部和南部出現(xiàn)明顯沉降中心。

3．3 全區(qū)累積地面沉降量分析

通過2004 年12 月至2010 年10 月共24 景ASAR 數(shù)據(jù)的PSInSAR 時序分析處理，得到6 年間工作區(qū)累積地面沉降量，如圖8 所示。由InSAR 監(jiān)測累積地面沉降量所示可知，在近6 年的時間內(nèi)，全區(qū)最大累積沉降量超過250 mm，兩個沉降區(qū)分別為濟寧市北部煤礦區(qū)和曲阜、兗州交界沿泗河一線。

圖6 2004．12 ～2006．02 InSAR 監(jiān)測地面沉降速率圖Fig．6 Land subsidence rate by applying InSAR from December 2004 to February 2006

濟寧、兗州轄區(qū)約4 060 km2的面積，去除大面積水域共約3 922 km2，通過對累積地面沉降量空間插值分析統(tǒng)計可得各沉降范圍影響面積，如表3 所示。

圖7 2008．09 ～2010．10 InSAR 監(jiān)測地面沉降速率圖Fig．7 Land subsidence rate by applying InSAR from September 2008 to October 2010

圖8 2004．12 ～2010．10 InSAR 監(jiān)測累積地面沉降量Fig．8 The cumulative land subsidence amount by applying In-SAR from December 2004 to October 2010

表3 累積地面沉降影響范圍統(tǒng)計表Tab．3 The statistical table of influence range about cumulative land subsidence

對濟寧、兗州各沉降范圍影響面積進行單獨統(tǒng)計，各地區(qū)主要沉降主要分布，如表4、表5 所示。

表4 濟寧市累積地面沉降影響范圍統(tǒng)計Tab．4 The statistical table of influence range about cumulative land subsidence in Jining city

表5 兗州市累積地面沉降影響范圍統(tǒng)計Tab．5 The statistical table of influence range about cumulative land subsidence in Yanzhou city

3．4 礦區(qū)范圍與沉降范圍比較

通過對累積地面沉降量空間插值分析得到該區(qū)域累積沉降等值線圖，將濟寧和兗州地區(qū)煤礦區(qū)礦權(quán)范圍數(shù)據(jù)與實際開采沉陷范圍進行疊加，得到如圖9 所示的累積沉降范圍與開采范圍。圖中左下角深灰色范圍為主要煤礦區(qū)礦權(quán)劃定范圍，比較發(fā)現(xiàn)，兩個主要沉降區(qū)分別對應(yīng)著兗州和濟寧地區(qū)兩個主要礦區(qū)范圍，其它無礦權(quán)范圍地區(qū)均為較小沉降區(qū)。對比沉降量的大小可知，大于60 mm沉降量的地區(qū)為煤礦開采區(qū)、整個沉降量為60 mm等值線所包圍的區(qū)域均為煤礦礦權(quán)覆蓋范圍;東山王樓煤礦沉降突出，中心區(qū)大于80 mm，影響范圍較廣。

圖9 工作區(qū)累積地面沉降范圍與礦區(qū)開采范圍Fig．9 The mining range and the cumulative land subsidence range in study area

4 結(jié)論

本次研究工作針對兗州和鄒城地區(qū)區(qū)域性地面沉降和煤礦開采區(qū)快速沉陷兩種類型地表形變，利用2004 ～2010 年多期ENVISAT ASAR 雷達數(shù)據(jù)，應(yīng)用PSInSAR 處理技術(shù)與常規(guī)D-In-SAR 差分干涉處理方法，獲取了全區(qū)2 個不同時間段的區(qū)域性地面沉降速率和6 年時間內(nèi)全區(qū)累積沉降量結(jié)果，并得到了多個時段礦區(qū)動態(tài)沉陷監(jiān)測成果。監(jiān)測結(jié)果查明了濟寧、兗州地區(qū)地面沉降分布狀況，統(tǒng)計得到了全區(qū)累積沉降量和影響地區(qū)范圍，全面掌握了該區(qū)域目前地面沉降的總體分布態(tài)勢和局部沉降特征。區(qū)域性地面沉降結(jié)果表明，濟寧、兗州的地面沉降主要與礦山開采相關(guān)，城市活動等其他因素帶來的地面沉降甚微;工作區(qū)內(nèi)地面沉降顯著區(qū)有兩個，一個為兗州市東南、曲阜市西部以及鄒城市西北部沉降區(qū)，另一個明顯的沉降區(qū)域為濟寧市任城區(qū)北部。比較兩期監(jiān)測成果，雖然最大沉降速率有所減小，但沉降嚴(yán)重區(qū)域(年沉降速率＞20 mm/a)有擴大的趨勢。

煤礦開采沉陷是上述區(qū)域沉降的主要因素，由于煤礦開采和開采過程的輸排水等，改變了工作區(qū)局部地質(zhì)環(huán)境，出現(xiàn)地面破壞。這種破壞與大規(guī)模開采地下水引發(fā)的區(qū)域性大范圍沉降漏斗不同，從空間分布上沉降中心比較集中，沉降速率比較大，局部甚至出現(xiàn)地裂縫和塌陷。因而，從監(jiān)測手段上InSAR 技術(shù)實現(xiàn)面上控制和城市全覆蓋，對于單個點上的連續(xù)快速沉降，還需進一步利用其它地面觀測的方式開展工作。

為了保證本研究工作的連續(xù)性和發(fā)展性，下一步工作中擬采用RADARSAT-2 衛(wèi)星獲取的精細模式SAR 數(shù)據(jù)開展相應(yīng)的監(jiān)測工作。該數(shù)據(jù)分辨率為8 m，優(yōu)于ENVISAT 數(shù)據(jù)，覆蓋范圍為50 ×50 km2，獲取周期為24 d，應(yīng)用PSInSAR時序分析方法，連續(xù)獲取數(shù)據(jù)既可滿足礦區(qū)動態(tài)監(jiān)測，也可滿足項目全區(qū)緩慢沉降監(jiān)測的需要。

［1］ 劉國祥．利用雷達干涉技術(shù)監(jiān)測區(qū)域地表形變［M］．北京:測繪出版社，2006:148．

［2］ 舒寧．雷達影像干涉測量原理［M］．武漢:武漢大學(xué)出版社，2003:99．

［3］ 程滔，單新建． CR、PS 干涉測量聯(lián)合解算算法研究［J］． 地震，2007，27(2):64 -71．

［4］ 范景輝，李梅，郭小方，等．基于PSInSAR 方法和ASAR 數(shù)據(jù)監(jiān)測天津地面沉降的試驗研究［J］．國土資源遙感，2007(4):23 -27．

［5］ 張潔．干涉雷達數(shù)據(jù)處理技術(shù)與應(yīng)用［D］． 武漢:中國地質(zhì)大學(xué)(武漢)，2005．

［6］ 韓春林．SAR 成像技術(shù)和SAR 圖像處理算法研究［D］．成都:電子科技大學(xué)，2004．

［7］ 楊成生． 基于D-InSAR 技術(shù)的煤礦沉陷監(jiān)測［D］． 西安:長安大學(xué)，2008．

［8］ 董玉森．干涉雷達中的永久散射體和交叉干涉技術(shù)的研究與應(yīng)用［D］．武漢:中國地質(zhì)大學(xué)(武漢)，2007．