葉萍萍

(甘肅省基礎(chǔ)地理信息中心，甘肅 蘭州 730000)

地面沉降作為一種常見的地質(zhì)災(zāi)害，具有持續(xù)時(shí)間長(zhǎng)、范圍大和影響因素復(fù)雜等特點(diǎn)。目前我國(guó)已有多個(gè)城市發(fā)生不同程度的地面沉降，對(duì)城市基礎(chǔ)設(shè)施造成重大安全隱患。

傳統(tǒng)的地表形變監(jiān)測(cè)主要是水準(zhǔn)測(cè)量、全站儀監(jiān)測(cè)等，雖然測(cè)量精度高，但是人工作業(yè)強(qiáng)度大、成本高，并且難以實(shí)現(xiàn)大范圍監(jiān)測(cè)[1-2]。合成孔徑雷達(dá)干涉測(cè)量技術(shù)作為新興的空間對(duì)地觀測(cè)技術(shù)，憑其全天時(shí)作業(yè)、監(jiān)測(cè)范圍廣、成本低、精度高等優(yōu)點(diǎn)，在城市沉降監(jiān)測(cè)等方面得到了廣泛的應(yīng)用。常規(guī)InSAR技術(shù)往往受限于時(shí)空失相關(guān)和大氣延遲的影響，形變探測(cè)精度較低，為克服以上問題，在此基礎(chǔ)上發(fā)展形成了多種時(shí)序InSAR技術(shù)，主要包括PS-InSAR、SBAS-InSAR。

蘭州新區(qū)作為西北第一個(gè)國(guó)家級(jí)新區(qū)，戰(zhàn)略地位十分重要，是西北地區(qū)重要的經(jīng)濟(jì)增長(zhǎng)極、國(guó)家重要的產(chǎn)業(yè)基地、承接產(chǎn)業(yè)轉(zhuǎn)移示范區(qū)[3]。新區(qū)近年為了促進(jìn)城市發(fā)展，一直推行平山造地，為城市儲(chǔ)備土地。因此，對(duì)新區(qū)，尤其是平山造地區(qū)域進(jìn)行地表形變監(jiān)測(cè)具有十分重要的意義。本次監(jiān)測(cè)選擇蘭州新區(qū)城區(qū)作為研究區(qū)域，選取2016年1月～2018年6月期間共28景Sentinel-1A雷達(dá)數(shù)據(jù)，采用SBAS-InSAR技術(shù)進(jìn)行差分干涉處理，生成區(qū)域形變圖,并根據(jù)結(jié)果分析形變?cè)?，為?zāi)害防治提供參考和依據(jù)。

1 研究區(qū)概述

蘭州新區(qū)坐落于甘肅省蘭州市北部的秦王川盆地，位于東經(jīng)103°26′20″～103°58′18″，北緯36°12′30″～36°52′27″之間[4]。本次研究區(qū)域范圍為蘭州新區(qū)城區(qū)，包含新區(qū)建成區(qū)以及新區(qū)的平山造地區(qū)域。具體位置如圖1所示。

圖1 研究區(qū)位置圖

2 數(shù)據(jù)來源

2.1 雷達(dá)影像數(shù)據(jù)

雷達(dá)影像為C波段的Sentinel-1A影像，來源于歐空局網(wǎng)站。共有28景，時(shí)間跨度為2016年1月～2018年6月，采用升軌模式，VV單極化，干涉寬幅(IW)模式的斜距單視復(fù)數(shù)(SLC)[5]，空間分辨率為5 m×20 m。雷達(dá)影像數(shù)據(jù)采集時(shí)間如表1所示。

表1 雷達(dá)影像數(shù)據(jù)采集時(shí)間

2.2 DEM數(shù)據(jù)

DEM數(shù)據(jù)采用國(guó)際熱帶農(nóng)業(yè)中心CIAT發(fā)布的90 m分辨率的SRTM DEM數(shù)據(jù)，用以消除由地形產(chǎn)生的相位。DEM的范圍必須大于雷達(dá)影像的范圍。

2.3 軌道數(shù)據(jù)

軌道信息是SBAS-InSAR數(shù)據(jù)處理中非常重要的信息,從最初的圖像配準(zhǔn)到最后的形變圖像生成都有著重要的作用，含有誤差的軌道信息造成基線誤差以殘差條紋的形式存在于干涉圖中[6]，使用衛(wèi)星精密軌道數(shù)據(jù)對(duì)軌道信息進(jìn)行修正，可有效去除因軌道誤差引起的系統(tǒng)性誤差。

雷達(dá)影像數(shù)據(jù)的衛(wèi)星軌道信息采用歐空局發(fā)布的POD精密定軌星歷數(shù)據(jù)。這是最精確的軌道數(shù)據(jù)，但該數(shù)據(jù)距離GNSS下行21天之后才可以使用，定位精度優(yōu)于5 cm。

3 技術(shù)方法和數(shù)據(jù)處理

3.1 技術(shù)方法

本次監(jiān)測(cè)采用短基線集差分干涉測(cè)量技術(shù)，即Small Baseline Subset InSAR(SBAS-InSAR)技術(shù)，它是在常規(guī)InSAR技術(shù)的基礎(chǔ)上發(fā)展的一種差分干涉測(cè)量技術(shù)，由Berardino等人在2002年提出[7]。為了克服常規(guī)InSAR技術(shù)由于時(shí)空基線過長(zhǎng)引起的失相干，SBAS-InSAR技術(shù)通過將雷達(dá)影像進(jìn)行適當(dāng)?shù)慕M合，選擇空間基線、時(shí)間基線都較小的InSAR干涉對(duì)，保證每幅干涉圖的高相干性。在保證相干性的同時(shí)，差分短基線集中包含的數(shù)據(jù)增加了時(shí)間采樣率，實(shí)現(xiàn)了連續(xù)觀測(cè)[8-9]。

SBAS-InSAR技術(shù)主要是利用奇異值分解(SVD)的方法將多個(gè)子集聯(lián)合進(jìn)行最小二乘法求解，獲得微小形變的信息[8-9]，在保留了傳統(tǒng) InSAR方法中大面積覆蓋和高時(shí)間分辨率的優(yōu)點(diǎn)的同時(shí)，提高了形變監(jiān)測(cè)的精度。

SBAS-InSAR具體技術(shù)流程如圖2所示。

圖2 SBAS-InSAR技術(shù)流程圖

3.2 數(shù)據(jù)處理

3.2.1 數(shù)據(jù)準(zhǔn)備

Sentinel-1A數(shù)據(jù)解壓后，導(dǎo)入軟件，生成強(qiáng)度圖數(shù)據(jù)(_pwr)、slc索引文件(.slc_list)等數(shù)據(jù)。由于下載的雷達(dá)影像范圍遠(yuǎn)大于蘭州新區(qū)的范圍，數(shù)據(jù)量較大，處理耗時(shí)，需將雷達(dá)影像進(jìn)行裁切。

3.2.2 生成連接圖

對(duì)輸入的數(shù)據(jù)進(jìn)行干涉像對(duì)的配對(duì)，輸入N景數(shù)據(jù)，能得到的最大配對(duì)數(shù)是(N×(N-1))/2，軟件會(huì)選擇最優(yōu)的組合方式進(jìn)行配對(duì)，也會(huì)自動(dòng)選擇超級(jí)主影像作為參考影像，所有的像對(duì)都會(huì)配準(zhǔn)到超級(jí)主影像上。

通過反復(fù)試驗(yàn)，調(diào)整參數(shù)，最后確定空間基線臨界基線最大百分比為2%，時(shí)間基線閾值為120 d，最后得到連接圖如圖3所示。

圖3 生成連接圖

3.2.3 干涉工作流

對(duì)所有配對(duì)的干涉像對(duì)進(jìn)行干涉處理，從相干性生成，去平、濾波和相位解纏，所有的數(shù)據(jù)對(duì)都配準(zhǔn)到超級(jí)主影像上。

依次查看各個(gè)像對(duì)的相干性圖和相位解纏圖。相干性圖顏色越黑，相干性越低，顏色越白，相干性越高。相位解纏圖也是一樣，白色區(qū)域越多越好。相干性圖、濾波后的干涉圖、相位解纏圖，只要其中一個(gè)結(jié)果較好，都可以保留。若有相干性低的，就用連接圖編輯工具移除該像對(duì)。

查看所有像對(duì)的干涉結(jié)果，相干性圖、濾波后的干涉圖和相位解纏圖的效果都比較好。

3.2.4 軌道精煉和重去平

估算和去除殘余的恒定相位和解纏后還存在的相位坡道。

選擇20171005-20170911這一像對(duì)作為輸入數(shù)據(jù)?？刂泣c(diǎn)要選擇沒有形變、穩(wěn)定的點(diǎn)，不能選擇有殘余地形條紋、形變條紋、相位躍變等地方。在進(jìn)行控制點(diǎn)選擇的時(shí)候，可以多選擇一些，軌道精煉和重去平做完后，查看控制點(diǎn)的誤差，刪除不好的控制點(diǎn)。

3.2.5 SBAS反演STEP1

這是SBAS反演的核心，第一次估算形變速率和殘余地形。二次解纏，用來對(duì)輸入的干涉圖進(jìn)行優(yōu)化?；谀Ｐ陀?jì)算出所有像對(duì)的形變和高程，提供線性、二次方、無位移、三次方四個(gè)模型，線性模型最穩(wěn)定，其他需要密集的連接圖和高相干性才能得到可靠的結(jié)果[10]。

3.2.6 SBAS反演STEP2

核心是計(jì)算時(shí)間序列上的位移，以反演第一步計(jì)算的形變速率為基礎(chǔ)，通過大氣濾波去除大氣相位，得到更加純凈的時(shí)間序列上的最終位移結(jié)果。大氣高通、大氣低通兩個(gè)選項(xiàng)，對(duì)大氣影響進(jìn)行估計(jì)，最后每個(gè)時(shí)間都從測(cè)量的位移中減去這些大氣部分。控制點(diǎn)使用軌道精煉時(shí)使用的控制點(diǎn)文件，用來移除恒定相位或者斜坡相位。

3.2.7 地理編碼

對(duì)SBAS-InSAR的結(jié)果進(jìn)行地理編碼，可以將地表形變結(jié)果投影到自定義方向上。

4 結(jié)果分析

圖4為2016～2018年蘭州新區(qū)城區(qū)地表形變速率圖。監(jiān)測(cè)結(jié)果表明：2016～2018年，研究區(qū)域形變速率為-54.18～18.23 mm/a，大部分區(qū)域形變速率在-4.39～4.66 mm/a之間，處于穩(wěn)定狀態(tài)，形變不明顯，主要形變區(qū)域集中在新區(qū)東部職教園區(qū)和東南部平山造地區(qū)域，且整體呈現(xiàn)下降趨勢(shì)。

圖4 蘭州新區(qū)城區(qū)地表形變速率圖

通過分析平山造地相關(guān)專題資料、查看遙感衛(wèi)星影像以及實(shí)地調(diào)查，可以得出新區(qū)地表沉降的主要原因是在新區(qū)東部和東南部進(jìn)行平山造地等工程開發(fā)，使得地面荷載加重，導(dǎo)致地面下沉。

5 結(jié) 論

本文通過對(duì)蘭州新區(qū)城區(qū)進(jìn)行時(shí)間序列的形變監(jiān)測(cè)研究，形成了SBAS-InSAR地表形變監(jiān)測(cè)的技術(shù)流程和監(jiān)測(cè)成果。

SBAS-InSAR技術(shù)監(jiān)測(cè)范圍廣、成本低、精度高、全天候作業(yè)等特點(diǎn)對(duì)城鎮(zhèn)地區(qū)地表形變監(jiān)測(cè)有很強(qiáng)的示范作用。使用SBAS-InSAR技術(shù)進(jìn)行時(shí)間序列監(jiān)測(cè)，能夠?yàn)榈刭|(zhì)災(zāi)害監(jiān)測(cè)和防治提供技術(shù)手段，實(shí)現(xiàn)城市地區(qū)地表形變動(dòng)態(tài)監(jiān)測(cè)，幫助有關(guān)部門動(dòng)態(tài)掌握城市地表形變的發(fā)生與變化規(guī)律，避免城市地質(zhì)災(zāi)害發(fā)生，對(duì)城市建設(shè)有著十分重要的意義。

在非城市地區(qū)，尤其是山區(qū)，由于地面實(shí)測(cè)數(shù)據(jù)缺少，還需進(jìn)一步加強(qiáng)SBAS-InSAR技術(shù)研究，爭(zhēng)取達(dá)到更好的監(jiān)測(cè)精度。