孫偉，李江衛(wèi)，白潔，童欣

(武漢市測(cè)繪研究院，湖北 武漢 430022)

1 引 言

近年來(lái)，隨著城市建設(shè)開(kāi)發(fā)力度加大，多種因素的綜合影響導(dǎo)致武漢市主城區(qū)部分區(qū)域陸續(xù)出現(xiàn)地面沉降。2013年以來(lái)，漢口香港路建設(shè)大道兩側(cè)、泛后湖區(qū)域、王家墩中央商務(wù)區(qū)及漢口火車站等區(qū)域相繼發(fā)生不同程度的地面沉降，影響了市民工作和生活。傳統(tǒng)的地面沉降監(jiān)測(cè)方法有衛(wèi)星定位測(cè)量、水準(zhǔn)測(cè)量、測(cè)量機(jī)器人等，其作業(yè)周期長(zhǎng)、空間采樣率低、接觸式測(cè)量、受天氣影響較大等劣勢(shì)導(dǎo)致沉降監(jiān)測(cè)工作效率低下，無(wú)法滿足當(dāng)前城市公共安全管理和應(yīng)急保障的需要。

針對(duì)城市地面沉降問(wèn)題，國(guó)內(nèi)外專家和大量工程技術(shù)人員開(kāi)展了很多有效的監(jiān)測(cè)工作，近20年逐步發(fā)展起的合成孔徑雷達(dá)干涉測(cè)量(Interferometric Synthetic Aperture Radar，InSAR)技術(shù)相對(duì)于傳統(tǒng)的大地監(jiān)測(cè)方法，具有非接觸式、高分辨率、高準(zhǔn)確率等優(yōu)勢(shì)，已成為地面、建筑物形變監(jiān)測(cè)應(yīng)用領(lǐng)域的研究熱點(diǎn)和發(fā)展方向。Ferretti等[1]基于差分合成孔徑雷達(dá)干涉測(cè)量(Differential InSAR，D-InSAR)理論，在二十世紀(jì)初提出了永久散射體干涉測(cè)量(Persistent Scatterer InSAR，PS-InSAR)，是利用長(zhǎng)時(shí)間序列SAR影像集進(jìn)行時(shí)間和空間域形變量估算，以提取永久散射體形變信息的干涉測(cè)量方法，該方法在監(jiān)測(cè)城市地面緩慢沉降中具有顯著優(yōu)勢(shì)。文獻(xiàn)[2]～[6]采用不同類型的SAR影響數(shù)據(jù)，利用PS-InSAR技術(shù)分別在北京、上海、天津、唐山等地面沉降監(jiān)測(cè)工作中進(jìn)行了有效試驗(yàn)，證實(shí)該技術(shù)的有效性、可靠性。為準(zhǔn)確掌握近幾年武漢市主城區(qū)地面沉降分布及動(dòng)態(tài)演化特征，本文利用武漢地區(qū)的36景COSMO-SkyMde影像數(shù)據(jù)，考慮到城市樓宇密集的特點(diǎn)，采用時(shí)序PS-InSAR技術(shù)進(jìn)行地面沉降監(jiān)測(cè)試驗(yàn)，為全面開(kāi)展自然資源監(jiān)測(cè)調(diào)查和城市公共安全管理探索新的技術(shù)手段。

2 監(jiān)測(cè)區(qū)范圍和數(shù)據(jù)處理

2.1 監(jiān)測(cè)區(qū)范圍

武漢地處江漢平原東部、長(zhǎng)江中游，長(zhǎng)江及其最大支流漢江橫穿武漢市主城區(qū)，形成了漢口、漢陽(yáng)、武昌的三鎮(zhèn)格局。屬殘丘性河湖沖積平原，黏土和亞黏土層較厚，受新構(gòu)造運(yùn)動(dòng)和河水流態(tài)的影響，特別是長(zhǎng)江三峽水利樞紐建成后沿江地區(qū)地下水位下降，改變了長(zhǎng)江流域季節(jié)性水位變化規(guī)律。近年來(lái)城市建設(shè)、開(kāi)發(fā)力度加大，高層建筑、地下工程建設(shè)等多種因素的綜合影響導(dǎo)致武漢市中心城區(qū)部分區(qū)域出現(xiàn)地面沉降。據(jù)調(diào)查，2013年以來(lái)，武漢市江岸區(qū)、江漢區(qū)、硚口區(qū)和武昌區(qū)共有40余個(gè)小區(qū)、單位相繼發(fā)生了地面沉降，造成建筑物附屬設(shè)施及市政道路不同程度開(kāi)裂、下沉，管道接頭脫節(jié)等，特別是江岸區(qū)后湖區(qū)域尤為突出[7]，調(diào)查與收集的地面沉降點(diǎn)共有157處、沉降路段3處。此外還有建設(shè)大道附近及其延長(zhǎng)線、竹葉山、江漢路、青年路、新華路、漢口火車站和武昌復(fù)興路附近的小區(qū)和單位均發(fā)生了不同程度的地面沉降現(xiàn)象。

圖1武漢市主城區(qū)監(jiān)測(cè)區(qū)域和數(shù)據(jù)源覆蓋范圍

2.2 PS-InSAR數(shù)據(jù)處理

本文采用覆蓋武漢市主城區(qū)范圍的意大利COSMO-SkyMed影像數(shù)據(jù)，監(jiān)測(cè)區(qū)域和數(shù)據(jù)源覆蓋范圍如圖1所示，基于PS-InSAR技術(shù)手段計(jì)算分析主城區(qū)范圍的地面沉降情況。SAR影像分辨率為 3 m，共計(jì)36景數(shù)據(jù)，時(shí)間區(qū)間為2012年6月～2017年6月，SAR影像數(shù)據(jù)基本參數(shù)見(jiàn)表1。COSMO-SkyMed影像數(shù)據(jù)具有覆蓋范圍廣、幾何分辨率高、入射角可選、多極化等特點(diǎn)，在城市建筑物密集區(qū)的地面沉降監(jiān)測(cè)方面效果顯著[8]。

監(jiān)測(cè)區(qū)域SAR影像數(shù)據(jù)參數(shù) 表1

根據(jù)PS-InSAR技術(shù)的基本原理[9]，可將PS-InSAR數(shù)據(jù)處理整理歸納為4個(gè)主要步驟：數(shù)據(jù)預(yù)處理、差分干涉計(jì)算、時(shí)空域形變估計(jì)和形變量計(jì)算，詳細(xì)流程如圖2所示。利用本文所述的COSMO-SkyMed影像數(shù)據(jù)源，采用意大利PSP-IFSAR軟件依據(jù)上述數(shù)據(jù)處理流程完成了數(shù)據(jù)處理工作，主要步驟包括：①公用主影像優(yōu)化選取。采用最佳相關(guān)性組合算法，計(jì)算了36景影像像對(duì)的時(shí)空基線，如圖3所示，最終選取的公用主影像為2015年12月2日的SAR雷達(dá)影像。②SAR圖像精確配準(zhǔn)。以選出的主影像為基準(zhǔn)圖像，將其他35幅影像都配準(zhǔn)到基準(zhǔn)SAR圖像的網(wǎng)格內(nèi)，配準(zhǔn)精度優(yōu)于0.1個(gè)像素。③DEM與主影像配準(zhǔn)、裁剪和組合，計(jì)算生成相位干涉圖，提取所有PS點(diǎn)的干涉相位。④時(shí)空相位解纏和回歸分析，計(jì)算分離大氣和軌道誤差相位。⑤PS點(diǎn)重新識(shí)別及形變速率反演和高程誤差估計(jì)。大氣和軌道誤差相位補(bǔ)償后，對(duì)SAR圖像的每個(gè)像素點(diǎn)進(jìn)行逐點(diǎn)分析，通過(guò)設(shè)置時(shí)間相關(guān)系數(shù)閾值，將在閾值以上的像素點(diǎn)選為最終的PS點(diǎn)，并估計(jì)其形變速率和高程誤差。⑥坐標(biāo)變換將視線向映射到垂直地面方向，得到所有高相干PS點(diǎn)的線性形變速率，采用內(nèi)插方法即得到區(qū)域線性形變速率。

圖2 PS-InSAR數(shù)據(jù)處理流程

圖3 時(shí)空基線分布圖

3 結(jié)果與分析

3.1 整體沉降監(jiān)測(cè)結(jié)果

根據(jù)計(jì)算結(jié)果，在武漢市主城區(qū)及其周邊 1 km范圍共提取到 7 571 260個(gè)PS點(diǎn)及其形變特征，包括了地理坐標(biāo)、形變發(fā)生時(shí)點(diǎn)和速度等信息。圖4和圖5分別給出了2012年6月～2017年6月期間武漢市主城區(qū)年平均形變速率圖、累計(jì)形變量圖，整體監(jiān)測(cè)的形變情況統(tǒng)計(jì)如下：

(1)通過(guò)監(jiān)測(cè)成果分析，獲取監(jiān)測(cè)區(qū)域內(nèi)主要的形變漏斗區(qū)域共計(jì)21處，分布區(qū)域詳見(jiàn)圖6，面積總計(jì)約 94.1 km2。面積較大的形變漏斗主要分布于主城區(qū)西北部，分別為：泛后湖區(qū)域、江漢區(qū)核心區(qū)、王家墩CBD連接漢西區(qū)域，對(duì)應(yīng)圖2中編號(hào)為5、6、7地塊，面積由依次為 22.3 km2、15.0 km2、10.4 km2。

(2)2012年～2013年期間，主城區(qū)整體形變不明顯，發(fā)生形變的區(qū)域主要集中在洪山區(qū)南湖周邊、青山區(qū)西部、江岸區(qū)南部、江漢區(qū)常青公園周邊、東西湖區(qū)東部，其累積形變量約 15.0 mm；武昌區(qū)域發(fā)生形變較大的青山區(qū)西部科技苑社區(qū)及周邊，形變時(shí)間區(qū)間主要為2013年6月～2015年6月；泛后湖區(qū)域形變發(fā)生時(shí)間為2014年6月左右，直到2017年6月形變趨勢(shì)仍未顯著減小。

(3)武漢市主城區(qū)監(jiān)測(cè)范圍內(nèi)PS點(diǎn)數(shù)量分布按照平均形變速率大小統(tǒng)計(jì)如表2所示，超過(guò)95%的PS點(diǎn)平均形變速率在 ±10 mm/a以內(nèi)；平均形變速率超過(guò) 27.0 mm/a的PS點(diǎn)共有 12 217個(gè)，集中分布在上述形變漏斗區(qū)域面積較大的部分。武漢市主城區(qū)監(jiān)測(cè)到的最大平均形變速率值為 -60.9 mm/a，位置在泛后湖區(qū)域的姑嫂村路與唐家墩路口處附近，累計(jì)形變量達(dá) -279.6 mm。

圖5 2012年6月～2017年6月武漢市主城區(qū)累計(jì)形變量圖

圖6 武漢市主城區(qū)21處形變漏斗區(qū)域分布示意圖

武漢市主城區(qū)InSAR監(jiān)測(cè)PS點(diǎn)平均形變速率范圍統(tǒng)計(jì)表 表2

續(xù)表2

3.2 重點(diǎn)區(qū)域沉降監(jiān)測(cè)結(jié)果

根據(jù)主城區(qū)整體監(jiān)測(cè)結(jié)果，結(jié)合現(xiàn)場(chǎng)實(shí)地巡查情況，對(duì)公眾關(guān)注的、沉降速率較大的泛后湖區(qū)域選取部分PS點(diǎn)進(jìn)行重點(diǎn)分析。

(1)江達(dá)路常青公園-金盾花園附近區(qū)域。圖7給出了該區(qū)域的PS點(diǎn)形變速率圖，形變速率較大。監(jiān)測(cè)期內(nèi)，最大形變速率為 24.7 mm/a，最大累計(jì)形變量為 -112.8 mm。圖7中A點(diǎn)位于金盾花園小區(qū)內(nèi)的建筑上，其形變歷史如圖8所示，五年間的平均形變速率為 24.7 mm/a，累計(jì)形變量為 -109.4 mm。在監(jiān)測(cè)期間內(nèi)，該點(diǎn)一直處于形變狀態(tài)，在2013年12月～2014年6月加速下滑，期間形變量為 -40.4 mm。

圖7 江達(dá)路常青公園-金盾花園附近區(qū)域PS點(diǎn)形變速率圖

圖8金盾花園小區(qū)建筑A點(diǎn)形變歷史曲線

(2)石橋一路新華書(shū)店周邊區(qū)域。圖9給出了該區(qū)域的PS點(diǎn)形變速率圖，形變速率較大。監(jiān)測(cè)期內(nèi)最大平均形變速率為 51.7 mm/a，最大累計(jì)形變量為 -242.2 mm。圖9中A點(diǎn)位于新華書(shū)店靠近文博路的出口處，其形變歷史如圖10所示，五年間的形變速率高達(dá) 51.7 mm/a，累計(jì)形變量 -242.2 mm。在整個(gè)監(jiān)測(cè)期間內(nèi)，該形變一直非常明顯，截止到2017年6月，形變趨勢(shì)仍未顯著減小。

圖9 石橋一路新華書(shū)店附近PS點(diǎn)形變速率圖

圖10新華書(shū)店出口處A點(diǎn)形變歷史曲線

3.3 地鐵沿線區(qū)域形變監(jiān)測(cè)結(jié)果

城市軌道交通工程修建、運(yùn)營(yíng)過(guò)程中，不可避免會(huì)對(duì)周邊產(chǎn)生影響，表現(xiàn)為地表或地表建(構(gòu))筑物的形變，因此對(duì)市政道路、軌道交通沿線區(qū)域進(jìn)行監(jiān)測(cè)，掌握和了解監(jiān)測(cè)區(qū)域內(nèi)地表形變速率和形變量對(duì)于保障人民生活安全意義重大。

本文研究的監(jiān)測(cè)時(shí)間段內(nèi)是武漢市地鐵建設(shè)的高峰期，依據(jù)武漢地鐵線路對(duì)監(jiān)測(cè)區(qū)域內(nèi)已運(yùn)營(yíng)和在建地鐵線路進(jìn)行矢量化，提取了已運(yùn)營(yíng)和在建地鐵線路沿線 200 m范圍內(nèi)的PS點(diǎn)，如圖11所示。結(jié)果顯示，武漢市主城區(qū)范圍內(nèi)已運(yùn)營(yíng)地鐵線路和在建地鐵線路多次穿過(guò)形變較大的區(qū)域，尤其是主城區(qū)西北側(cè)漢江、長(zhǎng)江與三環(huán)線所形成的三角形區(qū)域，如6號(hào)線石橋站周邊，3號(hào)線雙墩站、趙家條站周邊，8號(hào)線汪家墩站周邊等。其中，6號(hào)線石橋站周邊形變最為明顯，監(jiān)測(cè)期內(nèi)年平均形變速率可達(dá) 39 mm/a，累計(jì)形變量超過(guò) -180 mm。

圖11 武漢市主城區(qū)地鐵沿線2012年6月～2017年6月平均形變速率

3.4 泛后湖區(qū)域沉降成因分析

歷史上整個(gè)后湖地區(qū)均為湖灘，普遍隱埋一層一般厚 10 m以上、最厚可達(dá) 30 m的淤泥質(zhì)黏土或淤泥。由于近代人們填塘圍垸活動(dòng)，使湖泊周圍沼澤地帶的淤泥質(zhì)軟土裸露地表或被人工填土所掩蓋。該土層具有天然含水量高、孔隙比大、壓縮性高、強(qiáng)度低、滲透系數(shù)小的特點(diǎn)，具有觸變性、高壓縮性和流變性等工程地質(zhì)特征。高層建筑物的不斷加載、超深基坑降水工程手段的實(shí)施以及大氣降水補(bǔ)給的減少導(dǎo)致地下水位降低后，極易引發(fā)軟土地面沉降等地質(zhì)災(zāi)害。

隨著武漢市城市發(fā)展和經(jīng)濟(jì)建設(shè)的不斷進(jìn)行，地下水開(kāi)采造成的地下水位下降是主要外因，場(chǎng)地地表以下存在較厚的填土及軟土是主要內(nèi)因。造成地下水位下降的主要因素包括：一是近幾年，武漢段長(zhǎng)江洪峰減少，洪峰水位降低，且高水位維持時(shí)間較短，大量減少了長(zhǎng)江水的水量補(bǔ)給；二是近幾年超深基坑(開(kāi)挖深度進(jìn)入砂層)相繼開(kāi)工，深井降水措施造成上層滯水流失、水位下降，引發(fā)地表填土及軟土的固結(jié)沉降；三是近年來(lái)武漢地區(qū)降雨量也有所減少，地下水位隨之降低。

3.5 精度比較

為評(píng)估本文InSAR監(jiān)測(cè)成果精度，收集了監(jiān)測(cè)區(qū)域范圍、時(shí)間內(nèi)的精密水準(zhǔn)數(shù)據(jù)，分別是武漢地鐵6號(hào)線、8號(hào)線二等水準(zhǔn)多期監(jiān)測(cè)成果。經(jīng)過(guò)時(shí)空基準(zhǔn)統(tǒng)一，利用最鄰近點(diǎn)法的克里金插值方法對(duì)InSAR地表形變監(jiān)測(cè)成果精度進(jìn)行了分析和評(píng)估，計(jì)算統(tǒng)計(jì)了精密水準(zhǔn)值和InSAR監(jiān)測(cè)成果得數(shù)學(xué)精度一致性，如表3所示，根據(jù)平均誤差、中誤差、相關(guān)系數(shù)等參數(shù)指標(biāo)，PS-InSAR監(jiān)測(cè)結(jié)果精度優(yōu)于 ±5 mm，且與精密水準(zhǔn)結(jié)果吻合度較好，表明PS-InSAR監(jiān)測(cè)結(jié)果可靠，能實(shí)現(xiàn)對(duì)地面沉降情況的有效監(jiān)測(cè)。

精密水準(zhǔn)與InSAR監(jiān)測(cè)結(jié)果數(shù)學(xué)精度一致性統(tǒng)計(jì) 表3

4 結(jié) 語(yǔ)

本文開(kāi)展了高分辨率長(zhǎng)時(shí)間序列SAR影像監(jiān)測(cè)地表形變的研究，采用PS-InSAR技術(shù)對(duì)36景COSMO-SkyMed影像數(shù)據(jù)進(jìn)行干涉處理，獲取了武漢市主城區(qū)2012年6月～2017年6月的地面沉降信息。根據(jù)PS點(diǎn)分析，主要沉降區(qū)位于漢口的泛后湖區(qū)、江漢區(qū)核心區(qū)、王家墩CBD連接漢西區(qū)，其中泛后湖區(qū)域的沉降尤為嚴(yán)重，最大沉降速率達(dá)到 -60.9 mm/a，歷史上大面積淤泥質(zhì)黏土覆蓋層及地下水開(kāi)采是該區(qū)域發(fā)生較大沉降的主要原因。經(jīng)與高精度水準(zhǔn)觀測(cè)數(shù)據(jù)對(duì)比，驗(yàn)證了PS-InSAR技術(shù)在監(jiān)測(cè)武漢市地面沉降工作中的可行性、可靠性。

基于PS-InSAR的沉降監(jiān)測(cè)技術(shù)以其非接觸式、高分辨率、高準(zhǔn)確率等優(yōu)勢(shì)，已成為地表、建筑物形變監(jiān)測(cè)應(yīng)用領(lǐng)域的研究熱點(diǎn)和發(fā)展方向。在當(dāng)前服務(wù)自然資源“兩統(tǒng)一”和開(kāi)展自然資源監(jiān)測(cè)調(diào)查工作要求的大背景下，利用InSAR技術(shù)的優(yōu)勢(shì)和特點(diǎn)，研究和建立完整的變形監(jiān)測(cè)技術(shù)體系，對(duì)城市控制或防治地面沉降和地質(zhì)災(zāi)害工作具有重要的現(xiàn)實(shí)意義。