張桂芳
(中國地震局地質(zhì)研究所,北京 100029)
時(shí)序InSAR技術(shù)是近些年發(fā)展起來的一種新的InSAR測量技術(shù)。其中PS-InSAR技術(shù)是通過分析高相干的永久散射體相位變化信息提取形變,受時(shí)間和空間基線的限制小,該方法需要假設(shè)形變特征為勻速,同時(shí)需要大量的SAR數(shù)據(jù)(25景)來估算大氣效應(yīng)貢獻(xiàn)值;短基線技術(shù)(SBAS)可以降低幾何去相干的影響,同時(shí)較短基線使DEM誤差的敏感性大大降低,但該方法通過一定大小的窗體計(jì)算相干性,并據(jù)此來選擇時(shí)間序列穩(wěn)定點(diǎn),窗體操作降低了被選像元的分辨率,但抑制了噪聲相位的影響;干涉圖疊加方法(interferogram stacking)是將多幅差分干涉圖進(jìn)行疊加,提高結(jié)果中形變信息對大氣干擾的相對精度,可以獲取大區(qū)域的連續(xù)形變場;CR-InSAR技術(shù)將人工角反射器作為觀測對象,通過建立CR點(diǎn)上的形變模型,實(shí)現(xiàn)對重點(diǎn)目標(biāo)如滑坡、水庫和斷層的形變監(jiān)測。從以上的敘述可以看出,幾種時(shí)序InSAR技術(shù)各有優(yōu)缺點(diǎn),而本文最終的研究目標(biāo)是監(jiān)測京西北試驗(yàn)區(qū)的地表微小形變,所以論文選用以上幾種時(shí)序InSAR技術(shù)作為地表形變監(jiān)測手段,獲取覆蓋試驗(yàn)區(qū)不同尺度的形變速率場。在沒有GPS和水準(zhǔn)監(jiān)測資料的情況下,同一區(qū)域不同監(jiān)測技術(shù)手段得到的形變場之間可以相互檢驗(yàn),對比其可靠性和精度。圍繞這一目標(biāo),論文的研究工作主要聚焦到以下內(nèi)容。
在理想狀態(tài)下,雷達(dá)散射截面與有效散射面積和雷達(dá)波長有關(guān)。為更清楚地認(rèn)識(shí)在不同環(huán)境背景下,不同尺寸的角反射器對雷達(dá)信號實(shí)際的反應(yīng)特性,論文對角反射器的輻射強(qiáng)度進(jìn)行分析。以三角形角反射器為例,在完成角反射器設(shè)計(jì)制作、選點(diǎn)安裝和定標(biāo)以后,利用R-D模型和鄰域搜索算法對角反射器像素進(jìn)行了精確定位,然后統(tǒng)計(jì)各個(gè)角反射器對應(yīng)的輻射強(qiáng)度,并對角反射器安裝前的背景輻射強(qiáng)度與不同尺寸角反射器輻射強(qiáng)度進(jìn)行對比分析。通過分析可知,反射器安裝后輻射強(qiáng)度較背景輻射強(qiáng)度有明顯增加,而且輻射強(qiáng)度隨角反射器尺寸的增加而增大。角反射器安裝前,背景的平均輻射強(qiáng)度為214.117 1,即46.6 dB;邊長為1 m×1 m角反射器的平均輻射強(qiáng)度為646.314 6,即56.2 dB;邊長為1.2 m×1.2 m角反射器的平均輻射強(qiáng)度為902.059 0,即59.1 dB;減去輻射定標(biāo)常數(shù)45.1 dB,則背景和兩個(gè)不同邊長角反射器的輻射強(qiáng)度分別為1.5 dB、11.1 dB和14.0 dB。角反射器安裝前后輻射強(qiáng)度離散指數(shù)則隨角反射器尺寸的增大而減小,安裝前背景、邊長為1 m×1 m角反射器和邊長為1.2 m×1.2 m角反射器的輻射強(qiáng)度平均離散指數(shù)分別為0.305 7、0.133 7和0.116 8,這說明角反射器尺寸越大,角反射器輻射強(qiáng)度越穩(wěn)定。但角反射器相位在安裝后沒有表現(xiàn)出穩(wěn)定的特性,這可能與在SAR數(shù)據(jù)獲取和raw數(shù)據(jù)處理過程引起的相位誤差有關(guān)。論文還分析了角反射器所在像素及其周圍8鄰域像素的相位,發(fā)現(xiàn)在3×3窗口內(nèi)9個(gè)像素間的相位不具相關(guān)性,彼此之間的相位差很大,而且沒有任何規(guī)律,這說明角反射器相位受噪聲影響仍然嚴(yán)重。
角反射器的優(yōu)勢之一是能長時(shí)間保持相干性,而且相干性也不受空間基線的影響。所以可以充分利用數(shù)據(jù),在長時(shí)間大空間基線條件下分析CR點(diǎn)的形變;同時(shí)論文也在小基線集的情況下對CR點(diǎn)的形變做了分析。由第一部分內(nèi)容的分析可知,角反射器像素的相位受噪聲影響嚴(yán)重。為抑制噪聲相位的影響,在進(jìn)行干涉處理前,對角反射器像素及其周圍8鄰域像素的實(shí)部和虛報(bào)分別作了平均處理,然后對具有平均后實(shí)部和虛部的CR點(diǎn)集復(fù)共軛相乘。針對離散CR點(diǎn)差分干涉相位解纏的問題,在分析總結(jié)已有的角反射器相位解纏方法的基礎(chǔ)上,參考基于相對高程的相位解纏方法,提出了基于高程差的相位解纏方法。前者是基于CR點(diǎn)之間相對高程差進(jìn)行相位解纏,而后者是基于CR點(diǎn)實(shí)測高程與其對應(yīng)的DEM高程之間的差進(jìn)行相位解纏。試驗(yàn)區(qū)CR點(diǎn)差分干涉相位的解纏應(yīng)用證明了本論文提出的基于高程差的相位解纏方法的實(shí)用性。大氣水汽空間上的垂直分層使得大氣效應(yīng)的影響與地形相關(guān),大氣效應(yīng)在1~2 km的范圍內(nèi)是空間相關(guān)的,而試驗(yàn)區(qū)內(nèi)部分角反射器的距離超過了這個(gè)空間相關(guān)的范圍,針對這種情況,研究過程中用線性模型模擬大氣相位與地形的關(guān)系,從而減弱CR點(diǎn)上大氣相位誤差。通過CR點(diǎn)集相位分析模型得到了長時(shí)間大空間基線和小基線集兩種條件下的CR點(diǎn)形變速率場,并對其進(jìn)行對比分析。兩種條件下CR點(diǎn)形變趨勢基本相同,即CR點(diǎn)形變以視線向下降為主,形變速率為-2.0~0.0 mm/a。在沒有GPS和水準(zhǔn)形變資料的情況下,用兩種條件下分析得到的CR點(diǎn)形變速率相互比對,即可驗(yàn)證CR點(diǎn)形變的正確性,也有力地證明了角反射器的優(yōu)勢:在時(shí)序數(shù)據(jù)集中不受時(shí)間和空間基線的限制,能夠保持良好的相關(guān)性,為長時(shí)間大空間基線數(shù)據(jù)的利用提供了更好的條件。
論文用PS-InSAR技術(shù)和干涉圖疊加技術(shù)在京西北試驗(yàn)區(qū)做了應(yīng)用研究,分別獲取PS-InSAR覆蓋了70 km×45 km范圍的離散點(diǎn)形變速率場,和覆蓋盆地內(nèi)部(不包含水庫區(qū)域)及盆地內(nèi)邊緣的干涉圖疊加形變速率場。兩種技術(shù)獲取的形變速率場在重疊區(qū)域結(jié)果基本一致,即盆地內(nèi)部視線向下降,形變速率為-2.0~0.0 mm/a;盆地內(nèi)部的鎮(zhèn)級居住區(qū)比縣級居住區(qū)下降速率大,對應(yīng)的形變速率分別為-3.0~-2.0 mm/a和-2.0~-1.0 mm/a;盆地邊緣區(qū)的形變速率是-1.0~1.0 mm/a,說明斷裂帶的形變速率梯度范圍是0~2.0 mm/a;盆地外邊緣到山區(qū)的運(yùn)動(dòng)趨勢是視線向上升,PS-InSAR時(shí)序形變場顯示上升速率為1.0~4.0 mm/a,PS-InSAR技術(shù)得到的形變速率略大于干涉圖疊加形變速率,這可能是由PS-InSAR技術(shù)誤差傳播引起的。據(jù)方仲景等人(1993)對延礬盆地北緣斷裂的研究,其垂直形變速率為0.17~0.91 mm/a,本文給出延礬盆地北緣斷裂的形變速率與這一研究結(jié)果基本一致。