郭利民

（中國(guó)地震局地質(zhì)研究所，北京100029）

基于InSAR與多源數(shù)據(jù)的三維形變場(chǎng)獲取研究與應(yīng)用

郭利民

（中國(guó)地震局地質(zhì)研究所，北京100029）

中圖分類號(hào)：P315.72＋5；

文獻(xiàn)標(biāo)識(shí)碼：A；

doi：10.3969/j.issn.0235－4975.2015.08.008

合成孔徑雷達(dá)干涉測(cè)量技術(shù)（Interferometry Synthetic Aperture Radar，InSAR）近幾十年來得到迅猛發(fā)展，在獲取地面形變信息方面具有獨(dú)特的優(yōu)勢(shì)，使得人類可以在太空對(duì)地表進(jìn)行長(zhǎng)時(shí)間大范圍的監(jiān)測(cè)。但是，InSAR獲取的形變不是地表的真實(shí)變形，而是變形的垂直向、東西向、南北向分量在衛(wèi)星視線方向上的投影，即存在視線向（Line of Sight，LOS）模糊問題。對(duì)于地震、地殼運(yùn)動(dòng)、工程建設(shè)等變形特征的認(rèn)識(shí)，需要提取地表運(yùn)動(dòng)的三維特征，這不僅對(duì)深刻理解斷層運(yùn)動(dòng)方式、形變特征，分析地震發(fā)生機(jī)制有十分重要的作用，而且對(duì)人們有效監(jiān)測(cè)及預(yù)報(bào)地面沉降、塌陷等災(zāi)害同樣具有重要的實(shí)用價(jià)值。

本文詳細(xì)分析了InSAR形變監(jiān)測(cè)的視線向模糊問題及可能的三維形變場(chǎng)解決方案。根據(jù)InSAR和GPS的技術(shù)特點(diǎn)，分別提出和改進(jìn)了InSAR與GPS結(jié)合插值處理模型、多波段InSAR與多孔徑雷達(dá)干涉測(cè)量技術(shù)（Multiple－aperture SAR Interferometry，MAI）結(jié)合直接解算模型和聯(lián)合PSInSAR與GPS震間微小三維場(chǎng)模型。在此基礎(chǔ)上，以汶川地震、玉樹地震及海原斷裂帶為研究對(duì)象，獲取了同震和震間三維形變場(chǎng)，并詳細(xì)分析和論述了形變場(chǎng)特征，取得了一些新的認(rèn)識(shí)。

本文研究工作主要取得了如下成果和認(rèn)識(shí)：

（1）分析GPS三維形變分量和InSAR視線向觀測(cè)值之間的投影關(guān)系，探討了重軌In－SAR形變測(cè)量的局限性，即視線向模糊問題。針對(duì)各自的觀測(cè)特點(diǎn)，采用多波段InSAR、GPS、MAI等觀測(cè)手段及之間的結(jié)合，能夠在空間域和觀測(cè)優(yōu)勢(shì)方向上互相補(bǔ)充，不但可提高觀測(cè)精度，而且可獲取地表豐富的三維形變場(chǎng)信息。通過分析和研究，本文提出和改進(jìn)了3個(gè)三維形變場(chǎng)獲取模型。

①提出基于InSAR與GPS結(jié)合獲取三維形變場(chǎng)的插值處理模型。選取部分質(zhì)量較好的GPS觀測(cè)值，通過投影至衛(wèi)星觀測(cè)視線向，對(duì)同名點(diǎn)InSAR視線向觀測(cè)量進(jìn)行校正；利用GPS水準(zhǔn)觀測(cè)精度高這一優(yōu)勢(shì)，利用盡可能多的GPS水平觀測(cè)數(shù)據(jù)，通過插值方法獲取地表水平形變場(chǎng)；最后結(jié)合InSAR一維視線向觀測(cè)值，獲取同震三維形變場(chǎng)垂直方向分量?？紤]到斷層上盤或者下盤的形變矢量的大小和方向變化呈有序穩(wěn)定狀態(tài)，插值過程中先對(duì)地表

水平形變矢量的大小和方向進(jìn)行插值，然后再進(jìn)行EW和NS向分解，這比直接對(duì)EW和NS向形變進(jìn)行插值更接近真實(shí)的形變特征。

②改進(jìn)多波段InSAR與MAI結(jié)合直接解算模型。直接利用不同SAR衛(wèi)星的3種不同視線向的形變場(chǎng)，構(gòu)建三維形變解算模型進(jìn)行解算。其優(yōu)勢(shì)是具有較高的可信度，避免了模擬、假設(shè)等方法帶來的不確定因素。缺點(diǎn)是需要不同SAR衛(wèi)星，在一些研究區(qū)數(shù)據(jù)無法保障。

③提出PSInSAR與GPS震間微小三維場(chǎng)插值模型。針對(duì)地表微小形變運(yùn)動(dòng)提出PSIn－SAR與GPS聯(lián)合獲取三維震間形變場(chǎng)模型。采用PSInSAR方法獲取研究區(qū)的地表微小形變量，結(jié)合GPS水平形變觀測(cè)，依據(jù)插值算法，提出了獲取震間三維形變場(chǎng)模型。

（2）通過InSAR技術(shù)獲取了汶川、玉樹地震的同震形變場(chǎng)，分析了地震同震地表形變特征。

首先，通過InSAR獲取的汶川地震發(fā)震斷層周圍約450km×500km區(qū)域的同震形變場(chǎng)，形變場(chǎng)反映出汶川地震地表同震形變場(chǎng)的分布特征，即干涉條紋以北東向發(fā)震斷層（映秀—北川斷裂）為中心，且呈包絡(luò)狀分布。形變場(chǎng)南側(cè)甚至到達(dá)樂山、重慶一帶，北側(cè)到南坪、武都、康縣。從干涉條紋的分布格局可以看出，發(fā)震斷層兩側(cè)干涉條紋有多次起伏，且有一定的兩盤不對(duì)稱性，這說明位錯(cuò)量沿?cái)鄬臃植际遣痪鶆虻?。在發(fā)震斷層南西端映秀鎮(zhèn)附近，出現(xiàn)了排列細(xì)密的同心干涉條紋，說明這一帶的形變強(qiáng)度非常大。北川縣城南東側(cè)，有平行于斷層密集的干涉條紋包絡(luò)，表明該段發(fā)震斷層有較大的地表位錯(cuò)量。在汶川地震震中位置，條紋開始包絡(luò)震中，表明震中附近是發(fā)震斷層的西端點(diǎn)。在青川附近，條紋已包絡(luò)了發(fā)震斷層，顯示出了斷層?xùn)|端形成的形變場(chǎng)格局。發(fā)震斷層附近，斷層兩側(cè)均表現(xiàn)為局部抬升，且沿?cái)鄬有巫兞糠植己懿痪鶆?，表現(xiàn)出較強(qiáng)的分段性。

玉樹地震獲取的InSAR形變場(chǎng)條紋圍繞甘孜—玉樹斷裂北西向展布的橢圓狀干涉條紋，而且條紋在斷層兩側(cè)基本呈對(duì)稱分布。反映斷層兩盤形變梯度和形變量大體相當(dāng)，沒有明顯形變突出的主動(dòng)盤，這也正是剪切走滑斷層運(yùn)動(dòng)所具有的形變場(chǎng)特征。在形變場(chǎng)的西部和東部圍繞斷層跡線有兩個(gè)局部形變條紋環(huán)，而在這兩個(gè)局部形變之間有一小段未發(fā)生地表破裂的弱形變區(qū)。由形變條紋分布，我們可初步判定發(fā)震斷層長(zhǎng)約76km；東段約32km的斷層滑動(dòng)到達(dá)地表，形成較為明顯的地表破裂分布；西段則形成較為連續(xù)的形變條紋。InSAR形變場(chǎng)整體分布特征與野外調(diào)查及偏移量追蹤的結(jié)果均較為吻合。

（3）應(yīng)用基于InSAR與GPS三維插值處理模型，獲得汶川地震三維同震形變場(chǎng)，分析了三維形變特征，特別是斷層垂直形變特征。

通過三維形變場(chǎng)可以看出，發(fā)震斷層南段向北運(yùn)動(dòng)明顯，有一定的右旋走滑分量，龍門山發(fā)震斷層的逆沖分量由南向北逐步遞減，其主要逆沖形變量集中在映秀—連山坪、茶坪—北川縣城—南壩區(qū)段。同時(shí)注意到，在映秀鎮(zhèn)—連山坪一帶出現(xiàn)了東向位移極大區(qū)域，斷層上盤東向位移量達(dá)到了5.85m，斷層下盤向西位移量大約為1m。而在茶坪—北川縣城—南壩區(qū)段，斷層下盤向西位移量達(dá)到了－2.80m，斷層上盤向東位移量為2.9m。發(fā)震斷層南段上盤向東運(yùn)動(dòng)占絕對(duì)優(yōu)勢(shì)，而發(fā)震斷層中北段，則轉(zhuǎn)換為發(fā)震斷層上下盤對(duì)沖位移為主。斷層兩側(cè)垂直形變衰減較快，橫跨斷裂帶形變量大于30cm的寬度不超過50km；沿發(fā)震斷層附近垂直形變高值區(qū)分布不均勻，主要集中分布在發(fā)震斷裂的南段、中段和北端。

（4）應(yīng)用多波段InSAR與MAI結(jié)合直接解算模型，獲得玉樹地震斷層分布特征及震區(qū)三維形變特征。

利用ALOS和ENVISAT數(shù)據(jù)差分干涉得到的形變場(chǎng)及ALOS數(shù)據(jù)采用MAI方法得到