田 輝,孫岐發(fā),王福剛,馬詩敏,都基眾
1.中國地質調查局沈陽地質調查中心(沈陽地質礦產研究所),遼寧沈陽110026;
2.吉林大學環(huán)境與資源學院,吉林長春130026
D-InSAR在盤錦濕地地面沉降監(jiān)測中的應用
田 輝1,2,孫岐發(fā)1,王福剛2,馬詩敏1,都基眾1
1.中國地質調查局沈陽地質調查中心(沈陽地質礦產研究所),遼寧沈陽110026;
2.吉林大學環(huán)境與資源學院,吉林長春130026
采用二軌差分的方法,對濕地區(qū)的SAR數據進行配準、濾波、去平地效應、相位解纏、差分處理等處理,最后得到垂直形變圖,實現了地表形變的監(jiān)測.通過計算,得到東郭葦廠、歡喜嶺和西八千鄉(xiāng)3個沉降中心.3處沉降量分別為-169、-78和-105 mm.沉降面積方面,東郭葦廠(A)沉降面積為5.14 km2,橢圓形沉降區(qū)的長軸方向為北東-南西向.歡喜嶺(B)沉降區(qū)面積為0.42 km2,平面形態(tài)近似圓形.西八千鄉(xiāng)(C)沉降區(qū)面積為5.28 km2,橢圓形沉降區(qū)長軸方向為北東-南西向.
GAMMA軟件;二軌差分;合成孔徑雷達干涉測量;盤錦濕地;遼寧省
盤錦濕地區(qū)發(fā)現地面沉降現象,沉降造成葦場泵站報廢,給蘆葦種植、河蟹養(yǎng)殖造成了嚴重的影響,損失過億.利用D-InSAR技術,快速、準確的獲取沉降信息,從而為防災、減災工作提供依據,對濕地的利用與保護服務.
盤錦濕地位于下遼河平原最南端遼河入海口處,占地面積1000 hm2.盤錦濕地是我國最美的六大濕地之一,同時也是世界上植被類型保護最好的蘆葦沼澤地.保護區(qū)內有丹頂鶴、黑嘴鷗等國家一、二類保護動物43種,棲息的珍稀鳥類多達263種.因此,盤錦又有“鶴鄉(xiāng)”和“黑嘴鷗之鄉(xiāng)”之稱.每年9月,紅海灘美麗壯觀的景色吸引著大量國內外游客來此觀光[5-7].
研究區(qū)在大地構造上屬于華北臺地的東北部,在區(qū)域構造上位于遼河斷陷地帶.人類活動主要影響第四系、新近系地層,詳細論述如下.
第四系(Q):第四系地層平行不整合于新近系,并于裂谷兩側超覆角度不整合于前古近系巖層之上.其內部沉積連續(xù),地層一般厚65~400 m,水平結構特點由層次簡單向復雜過渡至層次不明顯,巖性由砂礫石層向細砂、粉細砂過渡.
新近系(N):研究區(qū)新近系平行不整合或微角度沉積不整合于古近系及前古近系裂谷基底地層之上,可分為明化鎮(zhèn)組、館陶組2個巖石地層單位.明化鎮(zhèn)組(Nm):下段為灰綠、黃褐、雜色泥巖、砂質泥巖夾灰白色砂巖、砂礫巖.半成巖,層理較發(fā)育、含植物殘屑.上段為灰白色,中粗?!屑毩I皫r、砂礫巖、含礫砂巖夾黃綠、灰綠色泥巖、泥質粉砂巖薄層或透鏡體.顆粒分選好,礫徑一般0.5~20 mm,大者可為70~100 mm,由北向南顆粒變細.頂部分布泥巖,厚度4~8 m,橫向變化不大,局部地段如高升、紅旗等地缺失.該組地層埋深一般在100~400 m之間,總體趨勢由東西兩側向中間,由北東向南西埋藏由淺至深,埋藏最深處位于趙圈河葦場以南.明化鎮(zhèn)組地層厚度一般為100~800 m,地層厚度變化趨勢與埋深一致.館陶組(Ng):巖性為灰、灰白色厚層狀含漂礫砂礫巖、細砂巖、薄層含礫砂巖及灰綠色、淺灰綠色、黃綠色泥巖、泥質粉砂巖,偶夾紫紅色泥巖.黏土質半膠結,礫石成分主要為石英、石英砂巖、火山巖、變質巖及少量花崗巖.礫石分選較差,北部、西部曙光、歡喜嶺地區(qū)礫徑較大,一般3~50 mm,大者可達150 mm以上,東部礫徑變細,一般3~10 mm.田莊臺、太平房地層中夾有玄武巖和玄武質凝灰?guī)r.該組地層在研究區(qū)西部部分缺失.缺失線南起艾家,經石新鎮(zhèn)、胡家延伸到北部的后王回窩堡.在缺失線附近,館陶組地層覆于花崗片麻巖、片麻巖、砂礫巖之上.地層由東西兩側向中間,由北東向南西厚度加大(50~370 m)、埋藏變深(350~1200 m).厚度最厚位于雙臺河口一帶可達370 m,埋藏最深處位于大洼縣西南海外河一帶,可達1200 m.
本文利用GAMMA軟件的MSP、ISP、DIFF&GEO模塊,采用的數據為ALOS PALSAR衛(wèi)星雷達數據,DEM分辨率為90 m.
二軌差分數據處理如圖1所示.利用主、從兩幅SAR影像,經過配準、干涉圖生成、濾波、去平地效應、相位解纏等步驟,得到干涉圖,再與DEM模擬的干涉圖進行差分處理,得到差分干涉干涉圖,最終得到地表形變圖[8-14].
(2) 巖棉復合型保溫模板體系圍護結構的平均傳熱系數K0為0.355 W/(m2·K),其高效的保溫性能,可滿足我國現行建筑墻體最高節(jié)能率要求。
圖1 兩軌差分干涉數據處理流程圖Fig.1 Two-pass differential interferometric data processing flowchart
3.1 影像的配準
在配準過程中查看配準誤差,按照經驗,只有小于0.2個象元才能繼續(xù)生成干涉圖.如果配準誤差過大,就需要繼續(xù)改進偏移多項式.通過改變影像塊的大小、影像塊的數目、閾值大小、配準多項式系數等來減小配準誤差.
3.2 干涉圖的生成
利用配準后的SLC影像,外加SLC的參數文件和偏移量的參數文件,生成研究區(qū)干涉圖.圖2a為2009年12月24日影像與2010年12月27日影像生成的研究區(qū)干涉圖.
3.3 干涉圖去平地效應
干涉成像時,平坦的地面產生干涉條紋,這些條紋的存在增加了解纏的難度.因此,應將平地相位從干涉圖中去除.
如圖2b所示,在生成的原始干涉圖(圖2a)中,最直觀的看只有細的彩色條紋.通過去平地效應處理,即能看出地形變化的趨勢.圖中存在顆粒狀斑點,即雷達噪聲斑點,需要進一步處理.
3.4 干涉圖濾波
干涉圖濾波可以提高解纏精度,有效地減少雷達噪聲帶來的誤差.具體計算過程如下:首先,利用復數干涉圖計算局部坡度,沿著斜坡通過計算得到干涉相位的平均值.然后,局部計算干涉圖像的功率譜,依據功率譜設計濾波器,對干涉圖進行濾波.圖2b、c是干涉圖濾波前、后對比圖.濾波前(圖2b)干涉圖面昏暗,有許多麻點;濾波后(圖2c)干涉圖面鮮艷,相位清晰,此時表明噪聲相位已經被有效去除.
圖2 2009-12-24—2010-12-27時段計算結果圖Fig.2 Calculation during 2009-12-24—2010-12-27a—干涉圖(interferogram);b—去平地相位(flat-earth phase removal);c—濾波(filtering);d—差分干涉圖(differential interferogram);e—相干圖(coherence map);f—垂直形變圖(vertical deformation map)
3.5 差分干涉圖的生成
具體步驟如下:先利用外部DEM數據模擬地形相位,然后采用最小二乘擬合方法優(yōu)化模擬地形相位和真實解纏相位影像之間的相位縮放因子,最后生成差分相位圖,即為去地形相位圖.圖2d為去地形相位的干涉圖.
3.6 相干圖的生成
相干圖的質量,直接決定了InSAR測量結果的質量.整個區(qū)域的相干系數越高,相干性越好,相干圖質量也就越高.圖2e為研究區(qū)的相干圖.其中黃色、紅色的區(qū)域為高亮區(qū),相干系數較高,能達到0.8左右,表明相干性較好.
3.7 形變圖生成
圖2f為2009-12-24~2010-12-27時段研究區(qū)垂直形變圖,該圖顯示出3個明顯的沉降區(qū),分別為東郭葦廠(A)、歡喜嶺(B)和西八千鄉(xiāng)(C).圖中一個顏色周期表示200 mm垂直形變量,A、B、C三處沉降量分別為-169、-78和-105 mm,平均沉降速率分別為167.6、77.3和104.1 mm/a.經過計算,東郭葦廠(A)沉降面積為5.14 km2,橢圓形沉降區(qū)的長軸方向為北東-南西向;歡喜嶺(B)沉降區(qū)面積為0.42 km2,平面形態(tài)近似圓形;西八千鄉(xiāng)(C)沉降區(qū)面積為5.28 km2,橢圓形沉降區(qū)長軸方向為北東-南西向.
3.8 形變產生的危害
1)對于泵站的危害:現有泵站持續(xù)下沉,完全淹沒在水中,已無法正常使用.
2)對于灌渠的危害:灌渠橫穿沉降中心區(qū),沉降導致灌渠局部下沉,靠天然的水力梯度,河水已無法輸入到蘆葦生長區(qū),必須要加高護堤,并利用大功率水泵,將河水輸送到地勢高的地方.
3)對于蘆葦生產的影響:蘆葦生長區(qū)大部分區(qū)域靠人工灌溉,河水到不了地勢高(相對于沉降區(qū))的地方,對于蘆葦生長極為不利,直接導致蘆葦減產,葦廠收入降低.
綜上所述,在沉降區(qū)必須投入更多的人力、物力、財力,才能維持葦廠生產及局部的生態(tài)平衡.
基于二軌法模型,對盤錦濕地地區(qū)ALOS PALSAR影像數據進行差分干涉測量數據處理,并從差分干涉圖中提取了地面沉降信息.得到東郭葦廠(A)、歡喜嶺(B)、西八千鄉(xiāng)(C)3個沉降中心的沉降速率及沉降面積信息.在濕地中,人類生產、建筑、經濟活動,應遠離沉降區(qū),減少沉降對人類的損失.
[1]穆冬.干涉合成孔徑雷達成像技術研究[M].南京:南京航空航天大學出版社,2001:32—35.
[2]張紅.基于相干目標的D-InSAR方法研究[M].北京:測繪出版社, 2009:1—20.
[3]周春霞.星載SAR干涉測量技術及其在南極冰川地貌地形研究中的應用[D].武漢:武漢大學,2004:23—30.
[4]單新建,馬瑾,宋曉宇,等.利用星載D-InSAR技術獲取的地表形變場研究張北-尚義地震震源破裂特征[J].中國地震,2002,18(1):119—126.
[5]陶梅,金寶丹,鈕勁濤.盤錦濕地功能開發(fā)與合理利用[J].遼寧工程技術大學學報,2009,28(3):84—86.
[6]李琦.盤錦濕地的價值及其開發(fā)利用[J].現代農業(yè)科技,2010(6):316—318.
[7]田輝.D-InSAR技術在盤錦地區(qū)地面沉降監(jiān)測中的應用研究[D].長春:吉林大學,2014:30—54.
[8]WANG Gui-jie,XIE Mo-wen,CHAI Xiao-qing.D-InSAR-based landslide location and monitoring at Wudongde hydropower reservoir in China[J]. Environmental Earth Sciences,2013,69(8):2763—2777.
[9]YANG Cheng-sheng,ZHANG Qin,ZHAO Chao-ying,et al.Monitoring mine collapse byD-InSAR[J].MiningScience and Technology,2010,20:696—700.
[10]FAN Hong-dong,DENG Ka-zhong,JU Cheng-yu.Land subsidence monitoring by D-InSAR technique[J].Mining Science and Technology(China),2011,21:869—872.
[11]何秀鳳.InSAR對地觀測數據處理方法與綜合測量[M].北京:測繪出版社,2012:110—150.
[12]王超,張紅,劉智.星載合成孔徑雷達干涉測量[M].北京:科學出版社,2002:55—78.
[13]廖明生.雷達干涉測量-原理與信號處理基礎[M].北京:測繪出版社,2003:43—96.
[14]郭華東.雷達對地觀測理論與應用[M].北京:科學出版社,2000:51—88.
APPLICATION OF D-InSAR IN THE LAND SUBSIDENCE MONITORING FOR PANJIN WETLAND
TIAN Hui1,2,SUN Qi-fa1,WANG Fu-gang2,MA Shi-min1,DU Ji-zhong1
1.ShenyangInstitute ofGeology and Mineral Resources,CGS,Shenyang 110034,China;
2.College of Environment and Resources,Jilin University,Changchun 130026,China
The two-pass D-InSAR method is adopted for the registration,filtering,flat-earth phase removal,phase unwrapping,differential processing of the SAR data in Panjin wetland,to obtain the vertical deformation map and realize the land subsidence monitoring.By calculation,the subsidence amounts of three subsidence centers,namely Dongguoweichang, Huanxiling and Xibaqianxiang,are-169,-78 and-105 mm,respectively.In respect of subsidence areas,Dongguoweichang(A)is 5.14 km2,in oval settlement area with NE-SW trending long axis;Huanxiling(B)is 0.42 km2,nearly circular plane shape;Xibaqianxiang(C)is 5.28 km2,in oval settlement area with NE-SW trending long axis.
GAMMA software;two-pass differential;interferometric synthetic aperture radar(InSAR);Panjin wetland; LiaoningProvince
1671-1947(2015)05-0507-04
P642.26
A
2014-06-25;
2015-01-04.編輯:李蘭英.
中國地質調查局“東北平原地面沉降調查”項目(編號:12120113011500)資助.
田輝(1984—),男,碩士,主要從事水文地質、環(huán)境地質調查研究工作,通信地址遼寧省沈陽市皇姑區(qū)黃河北大街280號,E-mail// 359585977@qq.com