摘 要：由于地下水開采和大量工程建設(shè)的影響，不少城市地面沉降問題突出。利用2011年至2013年間獲取的41景TerraSAR-X高分辨率SAR數(shù)據(jù)，采用相干點(diǎn)目標(biāo)分析技術(shù)對常州地表形變進(jìn)行監(jiān)測，獲得了地面沉降速率圖，并與水準(zhǔn)測量結(jié)果進(jìn)行了對比驗(yàn)證，達(dá)到了毫米級精度。研究結(jié)果顯示，常州市武進(jìn)地區(qū)存在幾處明顯的沉降，在2011-2013年間最大沉降速率達(dá)到28.7 mm/a，體現(xiàn)了IPTA技術(shù)在城市地面沉降監(jiān)測中有很好的應(yīng)用前景。

關(guān)鍵詞：InSAR IPTA 地面沉降 TerraSAR-X

中圖分類號：P64 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號：1674-098X（2015）04（b）-0009-03

常州位于長三角地區(qū)，自20世紀(jì)70年代，由于大量集中開采地下水，導(dǎo)致城市地區(qū)出現(xiàn)嚴(yán)重的地面沉降，是地面沉降災(zāi)害比較有代表性的城市之一[1]。傳統(tǒng)利用水準(zhǔn)測量監(jiān)測地面沉降的方法，雖然精度高，但存在勞動(dòng)強(qiáng)度大、覆蓋密度低且更新速度慢等不足。雷達(dá)差分干涉測量技術(shù)（differential interferometric SAR， D-InSAR）具有能定期獲取大范圍地表形變的優(yōu)勢，是監(jiān)測地面沉降的熱門技術(shù)。然而，城市地面沉降是一個(gè)緩慢的漸變過程，持續(xù)時(shí)間長且形變量較小，常規(guī)D-InSAR技術(shù)往往易受到時(shí)間、空間去相干和大氣延遲等影響。為了克服這些問題，許多學(xué)者開始研究時(shí)間序列分析方法，包括永久散射體雷達(dá)干涉測量（PS-InSAR）技術(shù)、短基線（SBAS）技術(shù)，等等，已在城市地表形變監(jiān)測方面取得了豐碩的成果[2-4]。

PS-InSAR技術(shù)是2000年由意大利學(xué)者 Ferretti[5]等提出的，該方法利用相位穩(wěn)定的點(diǎn)（PS）提取形變信息，從而避免了傳統(tǒng)DInSAR技術(shù)受時(shí)間和空間失相干的限制。SAR干涉點(diǎn)目標(biāo)分析技術(shù)（interferometric point target analysis， IPTA）是PSInSAR技術(shù)的一種，其在低相干區(qū)域仍能獲取一定數(shù)量的PS點(diǎn)，并僅對提取的點(diǎn)的時(shí)間維和空間維特征進(jìn)行分析[6]。該文以常州市武進(jìn)區(qū)為實(shí)驗(yàn)區(qū)，采用IPTA技術(shù)對其地表形變進(jìn)行監(jiān)測，探討IPTA時(shí)序方法的關(guān)鍵技術(shù)，并分析高分辨率TerraSAR-X數(shù)據(jù)在城市地表形變監(jiān)測中的優(yōu)勢和存在的問題。

1 干涉點(diǎn)目標(biāo)分析（IPTA）技術(shù)

IPTA技術(shù)提取時(shí)間序列中后向散射特性較為穩(wěn)定的相干目標(biāo)，對這些目標(biāo)的干涉相位進(jìn)行處理，在相位解纏的同時(shí)，根據(jù)差分相位的各組成部分在時(shí)間維和空間維的特征不同，分離不精確DEM所引入的誤差以及大氣效應(yīng)誤差，提取對應(yīng)點(diǎn)的時(shí)間序列上地表形變信息。該技術(shù)的主要原理如下：假設(shè)有N幅不同時(shí)相的SAR影像，選取其中一幅SAR影像作為公共主影像，其余N-1幅為輔影像，將輔影像依次與主影像配準(zhǔn)、干涉，并去除地平效應(yīng)，則可以形成時(shí)間序列干涉圖；利用外部DEM數(shù)據(jù)，進(jìn)行差分干涉處理，可得到N-1幅差分干涉相位圖，從而得到每個(gè)PS點(diǎn)的N-1個(gè)時(shí)序差分干涉相位。對于第k個(gè)干涉對的某一相干點(diǎn)目標(biāo)i而言，其相位模型可以表示如下[6]：

（1）

其中，wrap{}表示相位纏繞算子，表示相干目標(biāo)i的差分干涉相位觀測值，表示DEM高程誤差引起的地形相位殘差（可以用式（2）表示），表示沿LOS（line of sight）方向的地表形變相位，表示大氣延遲引起相位，表示噪聲相位。

（2）

式（2）中，λ為雷達(dá)波長，ri為雷達(dá)傳感器到地面目標(biāo)的斜距，θ為雷達(dá)入射角；B┴為干涉對的垂直基線；?h為DEM高程改正；令，可以根據(jù)SAR數(shù)據(jù)頭文件計(jì)算獲得。

地表形變相位包括線性形變相位和非線性形變相位兩部分，可以用式（3）表示。

（3）

其中，K2=4π/λ，t為干涉對時(shí)間基線，vi為沿雷達(dá)視線方向的線性形變速率；為非線性形變相位。

IPTA技術(shù)的具體過程如下：（1）根據(jù)時(shí)空基線最優(yōu)原則選擇主影像；（2）根據(jù)振幅特性與后向散射強(qiáng)度提取干涉點(diǎn)目標(biāo)，生成PS候選點(diǎn)目標(biāo)列表。一般為盡量不在初選中漏掉點(diǎn)目標(biāo)，初選條件一般設(shè)置較為寬松；（3）根據(jù)DEM數(shù)據(jù)和初始基線模擬生成平地相位和地形相位，并從每個(gè)點(diǎn)目標(biāo)的干涉相位中除去，生成SLC點(diǎn)數(shù)據(jù)，計(jì)算差分干涉相位；（4）對差分干涉相位進(jìn)行回歸分析，提取線性形變，回歸曲線斜率則代表了高程糾正值和相對線性變形速率；（5）利用空間濾波和時(shí)間濾波對殘余相位的成分進(jìn)行分離，對相位模型進(jìn)行迭代精化，解算高程糾正值和線性變形速率的最優(yōu)解；（6）對PS候選點(diǎn)進(jìn)行質(zhì)量控制，提取高質(zhì)量的PS點(diǎn)，獲取各個(gè)點(diǎn)在時(shí)間序列上形變速率和形變總量。

2 實(shí)驗(yàn)與分析

2.1 實(shí)驗(yàn)數(shù)據(jù)簡介

本次實(shí)驗(yàn)選取了2011年至2013年間41景TerraSAR-X高分辨率條帶模式SAR數(shù)據(jù)，這些數(shù)據(jù)均為單視復(fù)數(shù)影像（Single Look Complex， SLC），DEM數(shù)據(jù)使用SRTM-3數(shù)據(jù)。根據(jù)時(shí)間基線及空間基線最優(yōu)原則，選取20120509為公共主影像，其余為輔影像，將輔影像逐個(gè)與主影像進(jìn)行干涉可形成40對干涉對。實(shí)驗(yàn)區(qū)影像覆蓋范圍大約為900 km2，主要是常州市武進(jìn)區(qū)，該區(qū)域地勢平坦，人工建筑物較多。

2.2 形變結(jié)果與分析

利用IPTA方法對41景TerraSAR-X數(shù)據(jù)進(jìn)行處理，根據(jù)振幅特性與后向散射強(qiáng)度來提取PS候選點(diǎn)，首先提取PS候選點(diǎn)，并去除不穩(wěn)定的點(diǎn)，最終得到相干點(diǎn)目標(biāo)主要分布在密集建筑、橋梁等區(qū)域，水體、植被覆蓋區(qū)較少，說明相干點(diǎn)目標(biāo)提取正確；然后生成SLC點(diǎn)數(shù)據(jù)