利用ers12sar影像獲取西藏瑪尼地震的部分差分干涉圖

利用ers12sar影像獲取西藏瑪尼地震的部分差分干涉圖

0地表變形檢測(cè)合成孔雷達(dá)干擾測(cè)量技術(shù)（imar）是利用從合成孔中提取的相位信息為信息源獲得三維信息和變更信息的技術(shù)。它是利用空間上分開的兩副天線或同一天線在重復(fù)飛行的軌道上對(duì)同一區(qū)域進(jìn)行兩次成像,得到的兩幅復(fù)圖像(包括強(qiáng)度信息和相位信息)經(jīng)配準(zhǔn)后生成相位差圖像,利用相位差圖像來提取地面目標(biāo)的三維信息。雷達(dá)差分干涉測(cè)量(D-InSAR)由InSAR引申而來,它通過同一地區(qū)不同時(shí)相的兩幅干涉紋圖進(jìn)行差分組合或采取外部DEM模擬,消除地形影響,探測(cè)地表形變信息?？捎糜跈z測(cè)雷達(dá)視線方向厘米級(jí)或毫米級(jí)的地球表面形變,能夠全天候長(zhǎng)期獲取大面積、高精度、毫米級(jí)的緩慢地表變化信息,可以得到大范圍(數(shù)百千米×數(shù)百千米)的沉降趨勢(shì),在地震形變監(jiān)測(cè)、地面沉降監(jiān)測(cè)等方面已顯示了強(qiáng)大的優(yōu)越性。D-InSAR技術(shù)由于觀測(cè)結(jié)果的準(zhǔn)確性和觀測(cè)方法的特殊性,如其得到的觀測(cè)結(jié)果是整個(gè)區(qū)域面的變形信息(而其他相同精度的觀測(cè)手段得到的只是一些點(diǎn)和線的形變觀測(cè)信息),它能對(duì)人員不能到達(dá)的危險(xiǎn)地區(qū)的變形進(jìn)行高精度監(jiān)測(cè)等,引起了科技界和工程界的廣泛關(guān)注。國(guó)際上研究利用D-InSAR技術(shù)測(cè)量地震形變,開始于20世紀(jì)90年代初。1993年,Massonnet等利用ERS1SAR數(shù)據(jù)成功地應(yīng)用InSAR技術(shù)監(jiān)測(cè)到美國(guó)Landers地震的地表形變;1994年,Zebker又對(duì)Massonnet的方法提出了改進(jìn);Fujiwara等利用JERS-1星載InSAR技術(shù)監(jiān)測(cè)到日本伊豆半島的地殼形變。國(guó)內(nèi),王超等人于2000年利用InSAR技術(shù)獲取了張北—尚義地震的同震形變;單新建等利用InSAR技術(shù)對(duì)2001年昆侖大地震同震形變場(chǎng)進(jìn)行了分析;單新建等人還對(duì)西藏瑪尼地震震源的斷層參數(shù)進(jìn)行了研究。目前,D-InSAR技術(shù)已得到廣泛應(yīng)用,但是其需要解決的關(guān)鍵技術(shù)問題還有很多。同時(shí),由于受到數(shù)據(jù)源的限制,其用于地震形變的提取在國(guó)內(nèi)例子還較少,所以有必要展開更多的研究和應(yīng)用。本文在前人研究的基礎(chǔ)上,利用ERS-1/2SAR影像進(jìn)行了兩軌法D-InSAR提取地震形變的嘗試研究,并對(duì)獲得的結(jié)果進(jìn)行了分析。1d-smar數(shù)據(jù)處理地震前后獲取的雷達(dá)影像生成的干涉圖包含地形相位和形變相位。為了獲取形變相位,必須去除地形相位,具體方法有兩軌法、三軌法、四軌法。兩軌法是利用震前和震后兩景影像,進(jìn)行干涉處理,獲取干涉相位(包含地形相位和形變相位),利用外部DEM模擬在上述干涉像對(duì)的空間幾何條件下的地形相位,從干涉相位中減去之,就得到了形變相位。然后經(jīng)過相位解纏,解纏相位到形變量的轉(zhuǎn)換計(jì)算,地理編碼,得到地理坐標(biāo)下的形變量。圖1是兩軌法D-InSAR數(shù)據(jù)處理流程。兩軌法的關(guān)鍵是外部DEM與SAR主影像的配準(zhǔn)。本文采用瑞士GAMMA公司開發(fā)的GAMMA雷達(dá)干涉處理軟件,具體流程如下:1)獲取干涉圖:以震前1996年4月16日獲取的ERS-2影像為主影像,震后1998年4月21日獲取的ERS-2影像為輔影像。經(jīng)配準(zhǔn)、干涉、基線估計(jì)等步驟后得到干涉相位圖;2)SRTMDEM模擬地形相位:由于缺乏該地區(qū)高分辨率DEM,本文的外部DEM采用美國(guó)NASA公開的SRTMDEM來模擬地形相位,空間分辨率為3弧度秒,約為90m。首先將SRTMDEM與主影像進(jìn)行配準(zhǔn),轉(zhuǎn)換到雷達(dá)坐標(biāo)系下,模擬地形相位。3)獲取形變量:從步驟1)中獲取的干涉圖中減去地形相位,得到差分干涉圖,經(jīng)過相位解纏,解纏相位到形變的計(jì)算,地理編碼,最終獲得形變量。2地表同震變形場(chǎng)及數(shù)據(jù)處理1997年11月8日,在西藏自治區(qū)那曲地區(qū)瑪尼鄉(xiāng)北約150km處發(fā)生了MS7.9地震,震中位于φN35.2°,λE87.3°?，斈岬卣鸢l(fā)生在藏北高原北部長(zhǎng)達(dá)270km的瑪爾蓋茶卡斷層附近。野外考察資料表明,瑪尼地震的地表同震破裂帶長(zhǎng)120km,以左旋剪切為主,最大左旋位移量達(dá)4.5m,其東西兩側(cè)水平位移約2～3m,最大地裂縫帶寬度達(dá)300～400m。許力生等1999年的研究也表明,斷層性質(zhì)為左旋—逆走滑斷層,破裂方式為由西向東。西藏瑪尼地區(qū)氣候干燥,地表裸露,植被稀少,去相干因素少,非常適合利用D-InSAR技術(shù)提取地表同震位移場(chǎng)。本次實(shí)驗(yàn)選取了震前1996年4月16日和震后1998年4月21日獲取的ERS-2SAR影像進(jìn)行差分干涉,獲取同震形變信息。外部DEM采用美國(guó)NASA免費(fèi)提供的SRTM3弧度秒分辨率的DEM(見圖2),數(shù)據(jù)處理采用瑞士GAMMA遙感公司開發(fā)的GAMMA軟件。表1是所用SAR影像的基本參數(shù)。3區(qū)分變形條紋的意義由于西藏瑪尼地區(qū)植被稀少,地表裸露,雖然時(shí)間間隔約兩年,但大部分地區(qū)的相干系數(shù)均在0.5以上(見圖3),總體干涉情況很好。圖4給出了差分干涉條紋圖,從中可以看出,形變條紋清晰可見,地震引起的斷裂位置十分明顯,距離斷層越遠(yuǎn)條紋越稀疏,距離斷層越近則條紋越密集,并且差分干涉條紋呈現(xiàn)以震中為中心的平行向外拓展的模式,與地震形變規(guī)律相吻合。通過形變條紋的周期變化趨勢(shì)可以看出,斷裂帶以北的條紋向北東方向發(fā)散,斷裂帶以南向東南方向發(fā)散,形成左旋扭動(dòng)的態(tài)勢(shì),這一結(jié)果與野外考察所得結(jié)論一致。由于地震形變較大,所以差分干涉圖條紋清晰可辨,故根據(jù)條紋的數(shù)量定量的獲取視線向的形變量,也可根據(jù)衛(wèi)星的視角將視線向的形變量換算為水平形變量,這樣避免了相位解纏,不僅節(jié)省時(shí)間,而且避免了新的誤差引入。我們直接通過圖4中差分干涉條紋的色周數(shù)來定量的判斷由地震引起的雷達(dá)視線向形變量北盤最大約為0.7m,南盤最大約為0.8m。表明了兩軌法D-InSAR技術(shù)探測(cè)地震大范圍形變的巨大潛力和優(yōu)勢(shì)。需要指出的是,由于大氣效應(yīng)、基線估計(jì)誤差、外部DEM的誤差等的存在,我們認(rèn)為圖4中遠(yuǎn)離斷裂的上半部分為誤差,不是形變條紋。由于受數(shù)據(jù)的限制,本文只包含了瑪尼地震引起的部分?jǐn)嗔?研究具有區(qū)域性。而單新建等的研究結(jié)果覆蓋了整個(gè)地震斷裂,是全局的。4斷層上盤和下盤計(jì)算本文嘗試?yán)脙绍壏―-InSAR技術(shù)獲取了1997年瑪尼地震的部分形變,認(rèn)為由地震引起的雷達(dá)視線向形變量北盤最大約為0.7m,南盤最大約為0.8m。表明了利用D-InSAR技術(shù)獲取地震形變的巨大潛力和優(yōu)勢(shì)。研究得到的是視線向相對(duì)形變量,即假定離斷裂帶遠(yuǎn)處的某點(diǎn)形變量為0,得到的斷層上盤和下盤的最大形變量。如果有控制點(diǎn),則可以得到整個(gè)地區(qū)的絕對(duì)形變量。另外,研究得到的形變量是雷達(dá)視線向的,而實(shí)際的形變可能是任意方向的,所以需要根據(jù)經(jīng)驗(yàn)和地表實(shí)測(cè)數(shù)據(jù)建立形變模型模擬反演,才能確定形變機(jī)制。我國(guó)幅員遼闊,自然條件復(fù)雜,很多地區(qū)如青藏地區(qū)自然條件極其惡劣,很難也不太可能通過常規(guī)手段進(jìn)行監(jiān)測(cè),致使這些地區(qū)的監(jiān)測(cè)資料缺少。而InSAR技術(shù)不受這些因素的限制,具有全天候、全天時(shí)作業(yè)的優(yōu)勢(shì),并且這些地區(qū)植被稀少,更利于發(fā)揮InSAR技術(shù)的優(yōu)勢(shì)。InSAR技術(shù)的應(yīng)用使我