王 旋，鄭文青，苗小利

（中煤（西安）地下空間科技發(fā)展有限公司，陜西 西安710199）

1 概述

近年來(lái)，隨著我國(guó)城鎮(zhèn)化進(jìn)程的不斷加快，城市常住人口激增，城市占地面積急速擴(kuò)張，而在城市飛速發(fā)展中，不可避免帶來(lái)一系列的“城市病”。人地關(guān)系緊張，大中型城市發(fā)展迅速，而城市的擴(kuò)張路線(xiàn)大多是以老城區(qū)為中心，向四周地區(qū)開(kāi)發(fā)。而城市高層建筑群在自然地表上建設(shè)，會(huì)受到各種地質(zhì)災(zāi)害的影響，如山體滑坡和塌陷等。

為應(yīng)對(duì)可能發(fā)生的地質(zhì)災(zāi)害，應(yīng)當(dāng)建立常態(tài)化的監(jiān)測(cè)機(jī)制，以便對(duì)災(zāi)害高發(fā)地區(qū)提出預(yù)警。傳統(tǒng)的測(cè)量方法，如水準(zhǔn)測(cè)量和GPS測(cè)量，覆蓋密度低，監(jiān)測(cè)頻率低，空間分辨率小，不能對(duì)地表變化進(jìn)行高效率的監(jiān)測(cè)。InSAR技術(shù)在近十幾年發(fā)展迅猛，合成孔徑雷達(dá)差分干涉測(cè)量(DInSAR)技術(shù)的應(yīng)用更為廣泛，其探測(cè)地表形變精度高，范圍廣。具有全天時(shí)、全天候、分辨率高、以及主動(dòng)探測(cè)等優(yōu)點(diǎn)，在地震、地面沉降、火山和冰川活動(dòng)監(jiān)測(cè)中應(yīng)用非常廣泛[1]。本文選取Sentinel-1衛(wèi)星數(shù)據(jù)對(duì)蘭州市地表形變進(jìn)行監(jiān)測(cè)，Sentinel-1是歐空局發(fā)射的地球觀測(cè)衛(wèi)星，是哥白尼環(huán)境監(jiān)測(cè)計(jì)劃的重要組成部分，Sentinel-1A在2014年成功入軌運(yùn)行，此后對(duì)全球用戶(hù)提供數(shù)據(jù)服務(wù)，是InSAR技術(shù)重要的數(shù)據(jù)源之一。

2 研究區(qū)概況及數(shù)據(jù)選擇

2.1 研究區(qū)概況

蘭州市地處中國(guó)西北部的甘肅省，是西北地區(qū)的經(jīng)濟(jì)、政治、教育、文化中心。其地理位置十分特殊，黃河穿城而過(guò)，在較為平坦的河谷兩岸，蘭州市城關(guān)區(qū)開(kāi)始成形，此后再向其他地區(qū)擴(kuò)張。就形成了主城區(qū)沿狹長(zhǎng)的黃河河谷地帶分布發(fā)展的態(tài)勢(shì)，如圖1所示。近幾十年來(lái)，由于人口激增，城區(qū)面積由原來(lái)解放之初的16km2擴(kuò)張到超過(guò)220km2，但受特殊地形條件的限制，城市發(fā)展只能向河谷兩側(cè)山地和臺(tái)塬地區(qū)尋求空間。蘭州地區(qū)地表大部分被第四紀(jì)松散沉積物所覆蓋[2]。主城區(qū)所在的黃河河谷地帶大部分為濕陷性黃土，受水浸后極易發(fā)生下陷，給城市建筑帶來(lái)極大損害。

近年來(lái)，蘭州市城市建設(shè)工程對(duì)地質(zhì)災(zāi)害的影響急劇增加，主要就在于人地關(guān)系緊張，由此而引發(fā)的地質(zhì)災(zāi)害頻率愈高，范圍更廣，受特殊地質(zhì)條件影響，這些災(zāi)害問(wèn)題將一直存在[3-4]。據(jù)相關(guān)統(tǒng)計(jì)，在蘭州市，地質(zhì)災(zāi)害造成的人員傷亡和財(cái)產(chǎn)損失90%以上由人類(lèi)的建設(shè)工程引起，在施工區(qū)域200m的范圍內(nèi)，集中分布著該地區(qū)96%的地質(zhì)災(zāi)害，由此可見(jiàn)，蘭州市地質(zhì)災(zāi)害高發(fā)的關(guān)鍵因素就在于建筑工程對(duì)自然地表的破壞[5-7]。所以，建立長(zhǎng)效的地表變化監(jiān)測(cè)機(jī)制，刻不容緩。傳統(tǒng)的變形監(jiān)測(cè)方式雖然精度高，但是工程項(xiàng)目耗時(shí)耗力，所以采用時(shí)序InSAR技術(shù)對(duì)蘭州市主城區(qū)進(jìn)行長(zhǎng)時(shí)間序列的監(jiān)測(cè)，為地質(zhì)災(zāi)害提供動(dòng)態(tài)、實(shí)時(shí)、高精度的監(jiān)測(cè)成果。

2.2 數(shù)據(jù)選擇

本文選取2017年1月至2018年8月間的46景升軌Sentinel-1影像，軌道編號(hào)55，數(shù)據(jù)模式選擇干涉寬幅模式（IW），極化方式為同極化（VV），影像范圍如圖1所示。以SRTM-30數(shù)據(jù)作為參考DEM，精密軌道數(shù)據(jù)由歐空局軌道數(shù)據(jù)發(fā)布網(wǎng)站下載，所有影像相對(duì)于主影像的垂直基線(xiàn)及時(shí)間基線(xiàn)見(jiàn)表1。

表1 Sentinel-1A影像參數(shù)

3 SBAS-InSAR實(shí)驗(yàn)

3.1 SBAS-InSAR技術(shù)原理

傳統(tǒng)DInSAR技術(shù)受去相干因素及大氣延遲影響，且未能計(jì)算各影響因子對(duì)相位誤差的貢獻(xiàn)值，測(cè)量精度較低。2002年，Berardino和Lanari等人在分析了去相干因素對(duì)DInSAR精度影響的基礎(chǔ)上，總結(jié)出小基線(xiàn)集方法(Small Baseline Subset，SBAS)，其原理為設(shè)置時(shí)間與空間基線(xiàn)的閾值，選取小基線(xiàn)干涉影像對(duì)，將符合閾值的小基線(xiàn)集聯(lián)立，再利用奇異值分解(SVD)的算法求解，可得研究區(qū)長(zhǎng)時(shí)間序列地表形變的平均形變速率和時(shí)間序列形變量[8]。

實(shí)現(xiàn)的基本方法為，選擇N景SAR影像，按時(shí)間順序排列為 t0，…，ti，…，tN，以其中一景為主影像，做影像配準(zhǔn)，并設(shè)置時(shí)間和空間基線(xiàn)閾值，根據(jù)閾值范圍可將影像劃分為若干集合，分別在集合做差分干涉處理，設(shè)有N景影像分為L(zhǎng)個(gè)集合，共有M幅差分干涉圖，則M滿(mǎn)足：

所有的差分干涉干涉相位可組成觀測(cè)方程：

式中，IEk與ISk分別為主影像序列與從影像序列，是含有N個(gè)未知數(shù)的M個(gè)方程組，即：

式中，A為M×N階矩陣，；Φ為影像中高相干點(diǎn)對(duì)應(yīng)相位組成的向量；δΦ為各干涉解纏后相位組成向量。

當(dāng)SAR影像只有一組時(shí)，矩陣A的秩為N，對(duì)式（3）采用最小二乘法求解：

將式（5）代入式（2）中，可得：

即：

式中，B為M×N階矩陣，表達(dá)式為：

求解方程（8），首先對(duì)矩陣B做奇異值分解，求方程的最小范數(shù)最小二乘解，可得對(duì)應(yīng)地表在LOS方向的平均位移速率，即可得地表在LOS方向，時(shí)間序列上的形變量。

必須說(shuō)明的是，在式（7）中，δΦ還包含有大氣和其他噪聲引起的誤差，可利用時(shí)空濾波的方法來(lái)計(jì)算大氣對(duì)電磁波傳播的影響[9-10]。其基本原理為：大氣對(duì)電磁波的傳播的影響在空間域連續(xù)，表現(xiàn)為低頻信號(hào)；在時(shí)間域隨機(jī)，表現(xiàn)為高頻信號(hào)。通過(guò)時(shí)空濾波的方法進(jìn)行計(jì)算。

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

Sentinel-1A衛(wèi)星數(shù)據(jù)的IW模式，每景影像幅寬為250km，針對(duì)蘭州市主城區(qū)的地表形變監(jiān)測(cè)，選取中心城區(qū)約900km2區(qū)域。根據(jù)SAR影像的時(shí)間和空間基線(xiàn)報(bào)告，選擇獲取日期2018年1月9日的為主影像。空間基線(xiàn)閾值設(shè)為200m，時(shí)間基線(xiàn)閾值設(shè)為60d，符合基線(xiàn)閾值的影像組成干涉像對(duì)。影像連接圖與時(shí)間基線(xiàn)如圖2所示。

圖2 影像連接圖與時(shí)間基線(xiàn)

根據(jù)小基線(xiàn)集數(shù)據(jù)處理流程，（1）影像配準(zhǔn)，將所有影像配準(zhǔn)到20180109超級(jí)主影像上；（2）差分干涉，根據(jù)基線(xiàn)閾值劃分的集合進(jìn)行差分干涉處理，并以最小費(fèi)用流（MCF）法對(duì)結(jié)果進(jìn)行相位解纏，計(jì)算各個(gè)像對(duì)的相干性圖與解纏圖；（3）根據(jù)參考DEM和軌道參數(shù)，計(jì)算并去除平地相位及地形相位；（4）形變反演，利用SVD方法估算地表形變速率及殘余地形，利用時(shí)空濾波去除大氣延遲相位，獲得時(shí)間序列上地表平均形變速率；（5）地理編碼，將影像與成果由SAR坐標(biāo)系轉(zhuǎn)換到地理坐標(biāo)系下，分析實(shí)驗(yàn)結(jié)果。

經(jīng)過(guò)SBAS-InSAR數(shù)據(jù)處理，獲得了蘭州市主城區(qū)在2017年到2018年地表的年平均沉降速率結(jié)果，如圖3所示，其中紅色代表地面沉降，綠色代表地面抬升。

圖3 蘭州市區(qū)2017-2018年地表平均形變速率

4 實(shí)驗(yàn)結(jié)果分析

由圖3可得蘭州市主城區(qū)地表存在大范圍的沉降變形，在主城區(qū)所在的河谷地帶向南北兩側(cè)山坡擴(kuò)張的區(qū)域，地表沉降尤為顯著。在北部山區(qū)地帶形成了四大沉降區(qū)，由東到西分為A、B、C、D四個(gè)區(qū)域。其分布情況如地表形變速率圖所示。選擇各區(qū)域內(nèi)居民聚居地和典型地表形變區(qū)域，繪制其在時(shí)間序列上累積沉降曲線(xiàn)圖，如圖4所示。

圖4 ABCD四個(gè)區(qū)域內(nèi)沉降曲線(xiàn)圖

A區(qū)域內(nèi)，大浪溝村累積沉降量較小，年平均沉降速率為27.7mm/a，但是蘭州碧桂園沉降速率為104.2mm/a，直岔溝沉降速率為96.2mm/a，白道坪村沉降速率為72.1mm/a，石溝新村沉降速率為98.3mm/a，累積沉降量很大，其中尤以蘭州碧桂園的沉降速率和累積沉降量為甚。

B區(qū)域內(nèi)，恒大綠洲的年平均沉降速率為44.7mm/a，甘肅廣播電視大學(xué)職業(yè)技術(shù)學(xué)院的沉降速率為48.4mm/a，金城御景園的沉降速率為64.1mm/a，馨和園的沉降速率為 63.4mm/a，甘肅警察職業(yè)學(xué)院的沉降速率為82.9mm/a。

C區(qū)域內(nèi)，皋蘭縣第五中學(xué)的年平均沉降速率為31.8mm/a，大沙溝的沉降速率為44.6mm/a，安寧滑雪場(chǎng)的沉降速率為57.9mm/a，保利領(lǐng)秀山的沉降速率為64.7mm/a，蘭州萬(wàn)科城沉降速率為79.3mm/a。

D區(qū)域內(nèi)，李麻沙溝的年平均沉降速率為28.4mm/a，蘭州北車(chē)站的沉降速率為31mm/a，蘭州西機(jī)務(wù)段的沉降速率為37mm/a，白家鋪?zhàn)拥某两邓俾蕿?7.4mm/a，鳳凰山的沉降速率為52.1mm/a。

根據(jù)各區(qū)域內(nèi)的累計(jì)沉降曲線(xiàn)圖可以更清楚的觀察到，城區(qū)北部山區(qū)地帶沉降十分嚴(yán)重，且沉降區(qū)域在時(shí)間和空間上都有連續(xù)發(fā)展的趨勢(shì)。綜合來(lái)看，沉降速率及累積沉降量大的區(qū)域集中在新式居民區(qū)周?chē)驗(yàn)檫@些小區(qū)往往都是高層、密集的建筑物，而且都是近幾年新建成的居住區(qū)，所以對(duì)地表形變影響較大，如圖5所示。

圖5 城區(qū)內(nèi)地表沉降曲線(xiàn)圖

城區(qū)內(nèi)七里河區(qū)的地表沉降情況最為嚴(yán)重，城區(qū)向南部山區(qū)擴(kuò)張較多，其方式多為局部地形改造，或沿山谷地帶向兩側(cè)山脊方向延伸，所以沉降速率，累積沉降量及地表沉降范圍都是最大的。其中蘭州西站年平均沉降速率為19.1mm/a，黃河樓的沉降速率為24.9mm/a，明珠家園的沉降速率為24.8mm/a，彭家坪沉降速率為 30.7mm/a，五星坪街道沉降速率為35.1mm/a，而伏龍坪街道的沉降速率為54.2mm/a，沉降范圍最廣，沉降速率和累積沉降量最大，應(yīng)對(duì)可能發(fā)生的地質(zhì)災(zāi)害準(zhǔn)備應(yīng)急預(yù)案。

城關(guān)區(qū)是蘭州市的中心區(qū)，是各級(jí)黨、政、軍機(jī)關(guān)駐地，也是蘭州市和甘肅省的政治、科技、文化、教育、經(jīng)濟(jì)中心，人口數(shù)量最多。由于城區(qū)發(fā)展歷史長(zhǎng)遠(yuǎn)，建筑變遷及各種城市建設(shè)的影響，城區(qū)的地表形變也最為復(fù)雜，這一點(diǎn)由城關(guān)區(qū)地表沉降曲線(xiàn)圖可以看出。其中，華鼎中央都會(huì)的年平均沉降速率為24.2mm/a，蘭州高新開(kāi)發(fā)區(qū)總部經(jīng)濟(jì)二號(hào)園區(qū)沉降速率為19.2mm/a，城關(guān)區(qū)人民檢察院沉降速率為19mm/a，貢元巷周?chē)两邓俾蕿?8.2mm/a，蘭州東站調(diào)度站沉降速率為23mm/a。

5 結(jié)論

本文利用46景升軌Sentinel-1A影像數(shù)據(jù)，通過(guò)時(shí)序InSAR技術(shù)—短基線(xiàn)集干涉測(cè)量(SBAS-In-SAR)的方法來(lái)監(jiān)測(cè)蘭州市的地表形變情況。獲得了2017年1月到2018年8月的蘭州市主城區(qū)地表形變結(jié)果，并分析了形變區(qū)域的時(shí)空分布，總結(jié)了地表形變的規(guī)律，經(jīng)分析等到結(jié)論：

1）城區(qū)的南北部山地和臺(tái)塬地帶是新城區(qū)擴(kuò)張的主要方向，城市建設(shè)的大型建筑往往聚集在一個(gè)區(qū)域，工程活動(dòng)對(duì)自然地表的改造對(duì)地表形變有一定的影響，大型建筑物建成以后會(huì)對(duì)地面施加壓力，這兩者是造成新城區(qū)地表快速沉降的主要因素，區(qū)域內(nèi)地表恢復(fù)穩(wěn)定需要一個(gè)較長(zhǎng)的時(shí)期。

2）蘭州市主城區(qū)整體存在沉降變形，這與蘭州市的地質(zhì)環(huán)境有關(guān)，沿黃河谷地分布著濕陷性黃土，對(duì)城市建筑有很大的影響，且城區(qū)內(nèi)人口密集，城市建設(shè)和改造活動(dòng)比較活躍，多個(gè)方面共同作用使得地表形變異常復(fù)雜，但城市發(fā)展歷史較長(zhǎng)，大型建筑對(duì)地表的壓力趨于穩(wěn)定，所以相對(duì)于新開(kāi)發(fā)的區(qū)域沉降速率較小。

本文的實(shí)驗(yàn)成果分析與文獻(xiàn)內(nèi)蘭州市地質(zhì)災(zāi)害區(qū)域的劃分結(jié)果相符，可以對(duì)蘭州市內(nèi)的地質(zhì)災(zāi)害和城市地表監(jiān)測(cè)提供支持，可作為常態(tài)化監(jiān)測(cè)城市地表形變的技術(shù)手段，為防災(zāi)減災(zāi)、城市規(guī)劃與建設(shè)以及智慧城市發(fā)展提供科學(xué)依據(jù)。