周 呂 李佳豪 王 成 李 淑 朱子林 魯金金

1 桂林理工大學(xué)測繪地理信息學(xué)院,桂林市雁山街319號,541006 2 城市空間信息工程北京市重點實驗室,北京市羊坊店路15號,100038 3 廣西壯族自治區(qū)自然資源調(diào)查監(jiān)測院,南寧市建政路5號,530023 4 北京交通大學(xué)城市地下工程教育部重點實驗室,北京市上園村3號,100044 5 南寧理工學(xué)院,桂林市雁山街317號,541006

上海地區(qū)人口密度較大，大范圍地面沉降監(jiān)測難度大、任務(wù)重。目前針對上海地區(qū)沉降監(jiān)測的研究主要集中于局部監(jiān)測，且對基礎(chǔ)設(shè)施的系統(tǒng)性監(jiān)測較少；此外，各種實驗數(shù)據(jù)均較為久遠(yuǎn)，缺少時效性[1-4]。本文以PSInSAR技術(shù)為基礎(chǔ)，對上海地區(qū)2018～2020年35景Sentinel-1A影像數(shù)據(jù)進行處理與分析，獲取該地區(qū)及基礎(chǔ)設(shè)施的累積沉降量及沉降速率場，并結(jié)合城市化進程及基礎(chǔ)設(shè)施建成時間等信息，分析研究區(qū)及其內(nèi)部多條地鐵沿線的沉降時空特征，為基礎(chǔ)設(shè)施的維護提供科學(xué)的數(shù)據(jù)支撐。

1 研究區(qū)概況及數(shù)據(jù)源

1.1 研究區(qū)概況

上海(120°52′～122°12′E，30°40′～31°53′N)位于長江三角洲沖積平原，平均海拔為2.19 m左右，最高點海拔為103.70 m。截至2019年底，上海擁有軌道交通車站415座，線路共有17條(含磁懸浮線1條)；道路長度達(dá)4 227 km，道路面積為7 985萬m2；共有大小橋梁1 717座，公路橋梁8 482座；公路里程達(dá)10 392 km，其中高速公路通車?yán)锍虨?81 km。本文具體研究區(qū)域如圖1所示。

圖1 研究區(qū)概況

1.2 數(shù)據(jù)源

選取2018-01～2020-03覆蓋上海主城區(qū)的35景C波段升軌Sentinel-1A影像獲取沉降速率場及累積沉降量，Sentinel-1A具體參數(shù)見表1。利用 NASA 提供的 SRTM去除地形相位，采用ESA發(fā)布的精密軌道數(shù)據(jù)對影像數(shù)據(jù)進行軌道精化與相位重去平。

表1 影像源基本參數(shù)

2 研究方法

本實驗基于SARscape軟件，采用PSInSAR技術(shù)對數(shù)據(jù)進行處理，技術(shù)流程如圖2所示，具體處理流程如下[3,5-6]：

圖2 技術(shù)流程

1)將SAR影像集按時間序列排列，并綜合每景影像的時間基線、空間基線、多普勒質(zhì)心頻率差等3個不同因素，選取滿足最大相關(guān)系數(shù)條件的主影像作為本次處理的主影像。經(jīng)綜合計算，選取2019-04-06的SAR影像作為主影像。

2)將其他輔影像與主影像進行配準(zhǔn)。經(jīng)差分干涉處理，獲得差分干涉相位φi,j：

φi,j=φdef+φDEM+φorbit+φatm+φnoise

(1)

式中，φi,j為差分干涉相位；φdef為形變相位；φDEM為地形起伏造成的相位；φorbit為軌道誤差引起的相位；φatm為大氣延遲產(chǎn)生的相位差；φnoise為噪聲引起的相位。DEM引起的相位可由外部參考DEM進行消除，PSInSAR技術(shù)中較為穩(wěn)定的PS點可以降低φnoise的影響。

3)本文采用振幅離差指數(shù)法選取相干性較高的候選PS點，基于這些候選PS點的相位，利用空間域、時間域相位解纏方法分析相位變化過程，進而對相位信息進行解纏。對解纏后的差分相位進行分類，以消除其中的相位誤差。

4)基于候選PS點進行構(gòu)網(wǎng)，解算出整幅影像的高程殘差值、線性形變及殘余相位。對殘余相位進行時間域高通濾波和空間域低通濾波，得到每幅差分干涉圖的大氣相位與非線性形變估計值。

和普通輸尿管硬鏡相比較，輸尿管軟鏡可以到達(dá)許多硬鏡不能到達(dá)的腎盂及腎盞，與PCNL手術(shù)相比較又避免了通道建立對腎臟的損傷，因此近年來輸尿管軟鏡越來越多地應(yīng)用于腎臟及輸尿管上段結(jié)石的治療，其適應(yīng)證也越來越擴大，從結(jié)石大小、結(jié)石部位都有不同程度的放寬[7，8]。除此之外，術(shù)者的經(jīng)驗也越來越影響手術(shù)的效果，相同情況下，經(jīng)驗越豐富，手術(shù)時間越短，結(jié)石清除率越高。但輸尿管軟鏡的使用過程中，仍有一定比例的病例無法上鏡，不能順利進行輸尿管軟鏡下碎石取石，限制了輸尿管軟鏡的更大范圍應(yīng)用[9]。

5)去除大氣相位的影響后，將線性形變與非線性形變進行求和，得到上海地區(qū)形變速率場。

3 結(jié)果分析

3.1 研究區(qū)沉降速率分析

經(jīng)過PSInSAR技術(shù)處理后，得到上海地區(qū)2018-01～2020-03 LOS方向(遠(yuǎn)離雷達(dá)視線方向為負(fù)，靠近雷達(dá)視線方向為正)沉降速率圖(圖3)。從圖中可以看出，上海市主城區(qū)大部分區(qū)域沉降速率為-1～-10 mm/a，但在上海站附近區(qū)域出現(xiàn)較明顯的地面下沉現(xiàn)象，沉降速率在-1.9～-6.9 mm/a之間。在上海市南部，沉降現(xiàn)象呈現(xiàn)大范圍分布，且有多處較為明顯的沉降區(qū)域，并在閔行區(qū)滬金高速與申嘉湖高速附近出現(xiàn)較明顯的沉降漏斗，沉降速率均在-3.5 mm/a以上，最大沉降速率為-22.9 mm/a。

圖3 上海市2018～2020年地表形變速率

3.2 研究區(qū)地表沉降時空變化分析

為更好地分析上海地區(qū)沉降場的時空分布特征，選取研究區(qū)內(nèi)A、B、C、D四個特征區(qū)域進行沉降時間序列分析。特征區(qū)域空間分布、沉降情況分別如圖4、5所示，各區(qū)域時序點累積沉降量如圖6所示(時序點為圖5中綠色星形所在區(qū)域)。

圖4 特征區(qū)分布情況

圖5 特征區(qū)沉降情況

圖6 各區(qū)域累積沉降量

區(qū)域A為中山北路與南北高架橋交叉處，圖5(a)為該地區(qū)形變圖，從圖中可以看出，該區(qū)域沉降現(xiàn)象較為明顯，最大累積沉降量達(dá)-27.81 mm。該區(qū)域為上海市中心城區(qū)，具有眾多寫字樓與居民樓，還有鐵路及多條地鐵線路。水位監(jiān)測井?dāng)?shù)據(jù)顯示，2018-03地下水平均水位為-0.56 m，而2020-03該監(jiān)測井平均水位下降至-1.03 m，地下水水位呈現(xiàn)明顯下降趨勢。過度開采地下水及鐵路和地鐵線路的運行均會加大地面負(fù)荷，成為地面沉降的主要因素。

區(qū)域B位于上海市閔行區(qū)梅隴鎮(zhèn)許涇村附近，該區(qū)域沉降較為明顯，并形成沉降漏斗，最大累積沉降量達(dá)-47.3 mm，具體沉降速率如圖5(b)所示。2018～2020年，閔行區(qū)經(jīng)濟建設(shè)一直處于高速發(fā)展?fàn)顟B(tài)，而隨著城市化進程的加快，地面負(fù)荷也大大增加，隨之出現(xiàn)較為明顯的沉降現(xiàn)象。

區(qū)域C位于上海市浦東新區(qū)S32申嘉湖高速與滬南公路交叉處，該區(qū)域最大沉降量為-29.81 mm，沉降較為明顯, 具體沉降速率如圖5(c)所示。據(jù)資料顯示，該地區(qū)具有多個大型小區(qū)，居民生活用水主要來源于地下水。區(qū)域內(nèi)水位數(shù)據(jù)顯示，該地區(qū)2018-01地下水平均水位為-2.82 m，2020年地下水平均水位下降至-3.46 m，地下水的過度開采直接導(dǎo)致地面下沉。

區(qū)域D位于浦東新區(qū)川沙附近，該地區(qū)發(fā)展較為穩(wěn)定，累積沉降量也較為穩(wěn)定，主要集中在-5～3 mm之間, 具體沉降速率如圖5(d)所示?！镀謻|新區(qū)總體規(guī)劃暨土地利用總體規(guī)劃(2017～2035)草案》對該區(qū)人口規(guī)模及用地規(guī)模作出規(guī)定，要求嚴(yán)格控制人口總量、擴大生態(tài)空間和鎖定城鎮(zhèn)空間。該政策緩解了地下水的過度開發(fā)及城鎮(zhèn)化飛速發(fā)展造成的地面沉降現(xiàn)象，因此該地區(qū)未出現(xiàn)過大沉降現(xiàn)象。

3.3 地下水變化與地表沉降變化相關(guān)性分析

為分析上海地區(qū)地表沉降與地下水變化的相關(guān)性，獲取上海地區(qū)3個地下水水位監(jiān)測井的數(shù)據(jù)(水位監(jiān)測井分布如圖1所示)，對地下水水位監(jiān)測井周邊100 m內(nèi)PS點的累積沉降量求取平均值作為對比依據(jù)。其中，G1點地下水水位呈現(xiàn)持續(xù)下降趨勢，G2點地下水水位呈現(xiàn)平穩(wěn)變化趨勢，G3點地下水水位呈現(xiàn)穩(wěn)步上升趨勢，3處地下水水位呈現(xiàn)不同變化趨勢，可使對比結(jié)果更具代表性，地下水水位與累積沉降量對比結(jié)果如圖7所示。為進一步分析地表沉降變化與地下水水位變化的關(guān)系，去除地表沉降變化與水位變化中的線性趨勢后進行分析。計算結(jié)果表明，平均累積沉降量變化與地下水水位變化峰值存在一定相似性，并具有滯后特性，即地下水水位對地表沉降量存在直接影響。

圖7 地下水水位與累積沉降量對比

3.4 地鐵沿線地面沉降分析

截至2019年底，上海地區(qū)擁有地鐵線路16條，運營里程705 km，日均客運量1 065.03萬人次。面對如此大的客運流量，對上海地區(qū)地鐵線路及其沿線建筑物進行沉降監(jiān)測必不可少。本文研究區(qū)內(nèi)包含上海地區(qū)多條地鐵線路，具體如圖8所示。選取其中3條地鐵線路(5、8、16號)作為研究對象，將其周圍500 m作為緩沖區(qū)，分析地鐵線路及其周圍的形變場。

圖8 研究區(qū)地鐵線路

3.4.1 地鐵5號線(莘莊站-東川路站)

上海地鐵5號線是上海軌道交通第4條建成通車的線路,全長37.376 km。本文研究區(qū)內(nèi)包含地鐵5號線莘莊站至東川路站附近區(qū)域，提取得到兩站之間500 m范圍內(nèi)的沉降速率場(圖9)。由圖可知，該段沉降速率均大于-2.8 mm/a；其中，顓橋站至北橋站之間的沉降現(xiàn)象最為明顯，最大沉降速率達(dá)-12.0 mm/a，最大累積沉降量達(dá)-29.07 mm。該條線路于2003年開通運營，屬于營運時間較長的地鐵線路，但仍在莘莊站至東川路站之間出現(xiàn)較為明顯的沉降現(xiàn)象，這可能與“4改6”擴建改造工程有關(guān)。2018年后，上海地鐵5號線“4改6”擴建改造工程相繼完工，而站臺改建完成后的自然沉降也是該段發(fā)生明顯沉降的因素之一。此外，閔行區(qū)2018～2020年工業(yè)飛速發(fā)展及城市化進程的加快和隨之帶來的地下水過度開采也是該地區(qū)發(fā)生沉降的主要原因。

圖9 上海地鐵5號線沉降速率

3.4.2 地鐵8號線(市光路站-沈杜公路站)

上海地鐵8號線是上海市第6條建成運營的地鐵線路，呈東北-西南走向。截至2017-03，上海地鐵8號線全長37.5 km，共設(shè)30座車站，2019年日乘客量最高達(dá)122.1萬人次。

本文研究區(qū)內(nèi)包含地鐵8號線全線，具體沉降情況如圖10所示。從圖中可以看出，地鐵8號線北部相對較為穩(wěn)定，北部的青藏北路站與中興路站沉降較為明顯，最大沉降量為-15.65 mm。這兩站之間為上海市商業(yè)區(qū)，具有眾多大廈及工廠，而該地區(qū)的用水主要來源于地下水，這成為該區(qū)域地面下沉的因素之一。地鐵8號線南部沉降現(xiàn)象較為明顯，主要體現(xiàn)在凌兆新村站至沈杜公路站之間，其沉降速率均在-4.1 mm/a以上，最大累積沉降量達(dá)-24.44 mm。該地區(qū)具有眾多工業(yè)園區(qū)及居民區(qū)，其中以紡織業(yè)為主的浦江工業(yè)園區(qū)近年來飛速發(fā)展，基礎(chǔ)設(shè)施的建設(shè)成為該地區(qū)沉降的主要原因之一。

圖10 上海地鐵8號線沉降速率

3.4.3 地鐵16號線(龍陽路站-航頭東站)

上海地鐵16號線全長58.96 km，其中地下線路長13.69 km，共設(shè)站13座。本文研究區(qū)內(nèi)包含地鐵16號線部分線路(龍陽路站至航頭東站)，從圖11可以看出，該線路沉降主要集中在鶴沙航城站與航頭東站之間，最大累積沉降量達(dá)-35.96 mm。該地區(qū)位于特征區(qū)域C處，地下水的過度開采是該區(qū)沉降的主要原因之一。

圖11 上海地鐵16號線沉降速率

4 結(jié) 語

本文基于PSInSAR技術(shù)，利用35景Sentinel-1A影像獲取上海地區(qū)2018～2020年的沉降速率場與累積沉降量，并提取上海市地鐵5、8、16號線周圍500 m的沉降場，結(jié)合城市化進程數(shù)據(jù)與地下水開采等影響因素，對研究區(qū)地表沉降原因進行討論分析，主要研究結(jié)果如下：

1)上海市沉降現(xiàn)象不均勻性較為明顯，北部沉降相對平穩(wěn)，南部區(qū)域則出現(xiàn)較為明顯的大范圍沉降現(xiàn)象。閔行區(qū)許涇村沉降較為明顯，最大累積沉降量達(dá)-47.3 mm；浦東新區(qū)S32申嘉湖高速與122省道交叉處最大沉降量為-29.81 mm，兩地均存在較大形變現(xiàn)象。

2)以閔行區(qū)許涇村附近為中心的區(qū)域具有明顯的沉降漏斗，沉降速率均在-3.5 mm/a以上，最大沉降速率為-22.86 mm/a，該地區(qū)的沉降現(xiàn)象與區(qū)域基礎(chǔ)設(shè)施建設(shè)存在一定關(guān)系。

3)采用灰色關(guān)聯(lián)度分析方法對比分析地下水變化與地表沉降變化之間的相關(guān)性發(fā)現(xiàn)，上海地區(qū)地表沉降與地下水水位變化呈現(xiàn)一定的正相關(guān)性。

4)上海地鐵5號線顓橋站至北橋站之間的沉降現(xiàn)象最為明顯，最大沉降速率達(dá)-12.0 mm/a，該線出現(xiàn)明顯沉降與“4改6”擴建改造工程和閔行區(qū)的發(fā)展存在一定關(guān)系；地鐵8號線北部的青藏北路站與中興路站沉降較為明顯，最大沉降量為-15.65 mm，該線南部的凌兆新村站至沈杜公路站之間的沉降速率均在-4.1 mm/a以上，最大累積沉降量達(dá)-24.44 mm；地鐵16號線沉降主要集中在鶴沙航城站與航頭東站之間，最大沉降量達(dá)-35.96 mm。