陸燕燕，何 敏，何秀鳳

（河海大學(xué) 衛(wèi)星及空間信息應(yīng)用研究所，江蘇 南京 210098）

我國(guó)煤礦資源豐富，煤礦的開(kāi)采和利用，在產(chǎn)生巨大經(jīng)濟(jì)和社會(huì)效益的同時(shí)，也造成了大規(guī)模的地表沉降［1］。徐州作為我國(guó)重要的煤礦產(chǎn)地，近幾年由于煤炭資源的開(kāi)采以及地下水的過(guò)度抽取，已造成不同程度的地表沉降。徐州張雙樓煤礦是徐州重要的煤礦產(chǎn)地之一，位于徐州市西北的沛縣境內(nèi)，地理坐標(biāo)為116°46′～116°53′E，34°46′～34°51′N，該礦井田面積約34km2，礦區(qū)地形平坦，地勢(shì)西高東低。1979年1月開(kāi)工建設(shè)，1986年12月投入生產(chǎn)，礦井原設(shè)計(jì)年生產(chǎn)能力120萬(wàn)t，到目前為止核定年生產(chǎn)能力為225萬(wàn)t，是徐州礦務(wù)集團(tuán)的主力礦井之一。經(jīng)過(guò)30多年的開(kāi)采，張雙樓煤礦區(qū)及其周圍發(fā)生了不同程度的地表形變、民房開(kāi)裂、路表塌陷、礦區(qū)局部地震等地質(zhì)災(zāi)害問(wèn)題，嚴(yán)重影響了當(dāng)?shù)鼐用竦纳?，阻礙當(dāng)?shù)亟?jīng)濟(jì)的發(fā)展。因此，急需采取適當(dāng)?shù)姆椒▽?duì)張雙樓煤礦地面形變?yōu)暮M(jìn)行監(jiān)測(cè)、治理和修復(fù)。

傳統(tǒng)的礦區(qū)地表沉降監(jiān)測(cè)方法主要有水準(zhǔn)測(cè)量和GPS測(cè)量，這兩種方法都屬于基于點(diǎn)的監(jiān)測(cè)，其離散點(diǎn)監(jiān)測(cè)信息只能獲取空間離散點(diǎn)形變。對(duì)于礦區(qū)大面積微小形變監(jiān)測(cè)而言，需要布設(shè)數(shù)量龐大的監(jiān)測(cè)點(diǎn)，因而作業(yè)成本高，變形監(jiān)測(cè)工作量大，觀測(cè)周期長(zhǎng)，觀測(cè)結(jié)果也只反映了測(cè)站點(diǎn)的形變量，難以得到大范圍內(nèi)的整體形變結(jié)果。相比較而言，合成孔徑雷達(dá)差分干涉測(cè)量（DInSAR）作為一種新型的空間對(duì)地觀測(cè)技術(shù)，具有高精度（毫米至厘米級(jí)）、高分辨率（數(shù)十米）［2］、低成本、全天候、全天時(shí)等獨(dú)特優(yōu)點(diǎn)，且能夠?qū)ρ芯繀^(qū)域進(jìn)行大尺度連續(xù)覆蓋的面狀觀測(cè)，這使其在煤礦區(qū)地表形變的監(jiān)測(cè)中具有巨大的優(yōu)勢(shì)，無(wú)疑將成為一種未來(lái)煤礦區(qū)地面沉陷監(jiān)測(cè)極具潛力的空間對(duì)地觀測(cè)技術(shù)。

本文采用DInSAR技術(shù)對(duì)徐州張雙樓煤礦地區(qū)的地表沉降進(jìn)行監(jiān)測(cè)，實(shí)驗(yàn)數(shù)據(jù)為2011－01－16和2011－03－03獲得的兩景ALOS PALSAR影像。在數(shù)據(jù)處理過(guò)程中，采用DInSAR中的兩軌法，其中忽略大氣效應(yīng)對(duì)試驗(yàn)結(jié)果的影響，并采用SRTM DEM去除地形相位。

1 利用D－InSAR兩軌法監(jiān)測(cè)礦區(qū)的基本原理

DInSAR的監(jiān)測(cè)礦區(qū)基本原理是利用衛(wèi)星在形變前后對(duì)同一礦區(qū)分別成像，對(duì)兩景圖像進(jìn)行干涉處理，得到包含地形相位和形變相位的干涉圖，再利用已知的該地區(qū)的DEM，基于已有的成像參數(shù)模型模擬干涉紋圖，從而達(dá)到消除地形因素的效果［3］。DInSAR得到的相位可以用以下公式表示：式中：φf(shuō)lat為地球曲面引起的平地效應(yīng)相位，采用精密軌道數(shù)據(jù)可以去除；φtopo為地形起伏引起的地形相位，利用已知的DEM反演模擬地形相位消除；φatm為大氣延遲相位，φnosie為噪聲相位，這兩項(xiàng)可以通過(guò)濾波的方法給予削弱或去除；φdef為地表形變相位。圖1為兩軌法DInSAR的數(shù)據(jù)流程圖。

圖1 兩軌法DInSAR流程圖

在兩軌法DInSAR數(shù)據(jù)處理中，地形相位的去除對(duì)形變的結(jié)果影響較大，因此，獲得高精度的DEM至關(guān)重要。SRTM（Shuttle Radar Topography Mission）DEM數(shù)據(jù)由美國(guó)太空總署（NASA）和國(guó)防部國(guó)家測(cè)繪局（NIMA）聯(lián)合測(cè)量，獲取北緯60°至南緯60°之間總面積超過(guò)1.19億km2的雷達(dá)影像數(shù)據(jù)，覆蓋地球80%以上的陸地表面［4］。本文采用的SRTM3是SRTM DEM的產(chǎn)品之一，標(biāo)稱的絕對(duì)高程精度是±16m，相對(duì)高程精度是10m，標(biāo)稱的絕對(duì)平面精度是±20m，相對(duì)平面精度是15m［5－6］，目前中國(guó)地區(qū)獲得的SRTM3的分辨率是90m。

2 實(shí)驗(yàn)結(jié)果及分析

本文選取了2011－03－03和2011－01－16獲取的兩景ALOS PALSAR的影像進(jìn)行差分干涉。為了更好地考察46d間張雙樓礦區(qū)采礦引起的地表變化，截取了張雙樓煤礦所在區(qū)域?yàn)檠芯繉?duì)象，此區(qū)域包含七堡村、八堡村、陳莊、卓莊、南張莊及后周甲在內(nèi)的十多個(gè)村莊。圖2所示的矩形區(qū)域即為此次實(shí)驗(yàn)所截取的區(qū)域，試驗(yàn)區(qū)域面積約為6.089×8.238km2。表1為所選用的SAR影像的基本參數(shù)，46d的時(shí)間基線以及480m的垂直基線保證了圖像的良好相干性，從圖3中可以得出試驗(yàn)區(qū)域80%以上的區(qū)域相干值大于0.5，符合差分干涉測(cè)量的要求。

圖2 實(shí)驗(yàn)區(qū)域強(qiáng)度圖

圖3 相干系數(shù)圖

表1 SAR影像數(shù)據(jù)資料

本次實(shí)驗(yàn)中忽略大氣延遲的影響，主要原因有：①由表1的SAR數(shù)據(jù)獲取時(shí)間，查閱這一時(shí)期徐州地區(qū)的氣象資料顯示選取的這兩景數(shù)據(jù)時(shí)間段內(nèi)，徐州地區(qū)氣候干燥、少降雨、空氣清新、大氣能見(jiàn)度高、天氣以晴好為主；②本文的研究區(qū)域較小，為6.089×8.238km2，在這么小的范圍內(nèi)大氣狀態(tài)一般不會(huì)發(fā)生較大的變化。

兩幅圖像經(jīng)過(guò)配準(zhǔn)、重采樣、干涉處理后得到的干涉相位圖如圖4所示。此時(shí)的干涉相位包括噪聲相位、參考面相位、地形相位、形變相位。假定本次試驗(yàn)數(shù)據(jù)處理過(guò)程足夠精確，可忽略噪聲相位的影響；采用軌道數(shù)據(jù)和干涉條紋頻率相結(jié)合的方法估計(jì)基線，并根據(jù)成像幾何，反演并去除參考面相位。

圖4 未去除地形相位干涉圖

為了消除干涉相位中地形相位的影響，利用SRTM DEM來(lái)反演地形相位，通過(guò)兩軌差分法來(lái)消除［7］。本次實(shí)驗(yàn)中忽略了DEM誤差引入的相位誤差，這是因?yàn)楸疚氖褂玫母缮鎸?duì)的垂直基線為480m，SRTM DEM的相對(duì)高程誤差為10m，由高程模糊度可知，由DEM誤差引入的誤差相位達(dá)1/10個(gè)干涉條紋周期，即1.18cm。一般煤礦短時(shí)間內(nèi)沉降會(huì)達(dá)到幾十到上百厘米，再加上實(shí)驗(yàn)區(qū)地勢(shì)平坦（見(jiàn)圖5），所以地形誤差不會(huì)影響礦區(qū)沉降監(jiān)測(cè)［8］。

圖5 研究區(qū)域SRTM DEM（單位：m）

經(jīng)過(guò)一系列精細(xì)的差分處理后得到形變圖（見(jiàn)圖6），從圖6中可以看出，從2011－01－16到2011－03－03共46d的時(shí)間跨度內(nèi)，實(shí)驗(yàn)區(qū)有3處發(fā)生明顯的沉降，這3處區(qū)域分別位于七堡村、卓莊和后周甲村，此3處區(qū)域均位于張雙樓礦區(qū)內(nèi)，其分布與張雙樓礦井的分布一致。其中，七堡村地區(qū)干涉條紋明顯，沉降量較大，發(fā)生約420mm的形變量，形變面積約為0.49km2。其余兩處位于礦區(qū)附近的卓莊以及后周甲村附近，卓莊附近發(fā)生了約為95mm的沉降，沉降面積約為0.243km2，后周甲村發(fā)生約為92mm的沉降量，形變面積約為0.245km2。

圖6 形變圖（單位：m）

通過(guò)分析實(shí)驗(yàn)沉降結(jié)果，查閱張雙樓礦區(qū)近幾年煤礦開(kāi)采情況以及礦區(qū)地質(zhì)構(gòu)造等相關(guān)資料，得出張雙樓礦區(qū)是采空區(qū)塌陷的易發(fā)區(qū)［9］。由此推斷張雙樓礦區(qū)發(fā)生沉降的主要原因是由于礦區(qū)長(zhǎng)期的開(kāi)采，形成了大量的采空區(qū)，由于采空區(qū)的頂板坍塌以及煤柱破壞造成地表沉降導(dǎo)致這3處的沉降。

3 結(jié)束語(yǔ)

本文以徐州張雙樓礦區(qū)為研究區(qū)域，利用DIn－SAR兩軌法監(jiān)測(cè)張雙樓煤礦及其附近地區(qū)的沉降情況，獲得了張雙樓礦區(qū)的地表形變結(jié)果。結(jié)果表明，2011－01－16到2011－03－03期 間，張 雙 樓 礦 區(qū) 共出現(xiàn)3處沉降漏斗區(qū)域，漏斗中心最大沉降量達(dá)到420mm，其附近的卓莊沉降量達(dá)到95mm，后周甲沉降量92mm，并且這3處沉降漏斗區(qū)域與礦區(qū)及其采空區(qū)分布一致。

由于本次實(shí)驗(yàn)數(shù)據(jù)比較少，只能采用DInSAR方法來(lái)監(jiān)測(cè)礦區(qū)沉降變化，并未考慮時(shí)空去相關(guān)以及大氣非均勻性影響，使得地表形變的精度受到限制。隨著以后實(shí)驗(yàn)數(shù)據(jù)的增多，將采用永久散射體的雷達(dá)差分測(cè)量方法（PS－InSAR）進(jìn)一步研究，從而得到高精度的地表形變。

［1］何敏，陸曉燕，何秀鳳.利用D－InSAR二軌法監(jiān)測(cè)徐州大屯中心區(qū)地表形變［J］.地理空間信息，2011，9（5）：3－5.

［2］劉國(guó)祥，丁曉利.沿海地區(qū)DInSAR形變探測(cè)：精度與可應(yīng)用性分析［J］.測(cè)繪通報(bào)，2006（9）：9－13.

［3］王超，張紅，劉智.星載合成孔徑雷達(dá)干涉測(cè)量［M］.北京科學(xué)出版社，2002.

［4］何敏，何秀鳳.利用星載InSAR技術(shù)提取鎮(zhèn)江地區(qū)DEM及其精度分析［J］.計(jì)算機(jī)應(yīng)用，2010，30（2）：537－539.

［5］SRTM Documentation［EB/OL］.［2009－08－01］.ftp：//e0srp01u.ecs.nasa.gov/srtm/.

［6］陳俊勇.對(duì)SRTM3和GTOPO30地形數(shù)據(jù)質(zhì)量的評(píng)估［J］.武漢大學(xué)學(xué)報(bào)：信息科學(xué)版，2005，30（11）：941－944.

［7］何敏，何秀鳳.利用D－InSAR技術(shù)監(jiān)測(cè)鹽城地區(qū)地表形變［J］.測(cè)繪通報(bào)，2010（11）：1－3.

［8］黃寶偉，宋小剛，王振杰，等.基于D－InSAR技術(shù)的葛亭煤礦地面沉降監(jiān)測(cè)研究［J］.工程勘察，2012（4）：55－60.

［9］童立元，劉松玉，方磊.徐州礦區(qū)采空塌陷區(qū)公路工程地質(zhì)災(zāi)害與防治對(duì)策分析［J］.工程地質(zhì)學(xué)報(bào)，2005，13（2）：165－168.