肖波，趙振峰，喜文飛，周定義，姜乃齊

(1.云南交通職業(yè)技術(shù)學(xué)院，云南 昆明 650500； 2.昆明理工大學(xué)國(guó)土資源工程學(xué)院，云南 昆明 650093；3.云南師范大學(xué)地理學(xué)部，云南 昆明 650500； 4.云南省核工業(yè)地質(zhì)調(diào)查院，云南 昆明 650106)

1 引 言

近年來(lái)，為主動(dòng)融入國(guó)家提出的“一帶一路”和“長(zhǎng)江經(jīng)濟(jì)帶”發(fā)展戰(zhàn)略，落實(shí)習(xí)近平總書(shū)記視察云南的重要指示，云南省調(diào)整《云南省道網(wǎng)規(guī)劃(2014-2030年)》中的高速公路網(wǎng)布局方案，提出了“五縱五橫一邊兩環(huán)二十聯(lián)”的高速公路布網(wǎng)方案，為“做強(qiáng)滇中、搞活沿邊、聯(lián)動(dòng)廊帶、多點(diǎn)支撐、雙向開(kāi)放”的戰(zhàn)略做好交通運(yùn)輸服務(wù)。截至2020年底年，云南高速公路通車?yán)锍桃呀?jīng)超過(guò) 9 000 km。高速公路工程建設(shè)完成投入運(yùn)營(yíng)后，不可避免地會(huì)引發(fā)區(qū)域性的路面沉降，產(chǎn)生裂縫、裂隙、橋梁橋頭跳車等病害問(wèn)題，破壞路基路面的穩(wěn)定性，影響行車安全，危及人們的生命和財(cái)產(chǎn)安全[1]。滇中城市區(qū)，包括昆明、曲靖、楚雄和玉溪4個(gè)州市，位于云南的發(fā)展核心區(qū)，密集的路網(wǎng)規(guī)劃和建設(shè)，更是高速公路路面沉降的高發(fā)區(qū)。因此，迫切需要監(jiān)測(cè)滇中城市區(qū)高速公路的形變發(fā)展趨勢(shì)，掌握其路基路面結(jié)構(gòu)的穩(wěn)定性，保障高速公路的正常運(yùn)營(yíng)。

目前，高速公路的變形監(jiān)測(cè)已經(jīng)明顯地滯后于建設(shè)，如果繼續(xù)再利用傳統(tǒng)監(jiān)測(cè)方法，如水準(zhǔn)儀、經(jīng)緯儀、全站儀、GNSS、填埋降標(biāo)等，對(duì)大范圍的高速公路網(wǎng)進(jìn)行安全排查，不僅耗時(shí)、耗力，養(yǎng)護(hù)人員的人身安全也得不到保障，且容易出現(xiàn)排查不到位、遺漏等問(wèn)題。隨著空間衛(wèi)星的不斷升空，越來(lái)越多遙感數(shù)據(jù)的免費(fèi)使用，新興的InSAR技術(shù)以無(wú)須大量的外業(yè)工作，就能獲得大范圍內(nèi)毫米級(jí)的形變監(jiān)測(cè)精度而備受人們的青睞。從2010年開(kāi)始，陸續(xù)有國(guó)內(nèi)外學(xué)者將InSAR技術(shù)用于城市道路[2～4]、高速公路[5～8]、地鐵[9～12]等線性工程的形變監(jiān)測(cè)，驗(yàn)證了InSAR技術(shù)對(duì)線性工程形變探測(cè)的可靠性及應(yīng)用潛力。筆者檢索文獻(xiàn)發(fā)現(xiàn)，現(xiàn)有的研究中，眾多的學(xué)者選取的研究對(duì)象相對(duì)單一，僅對(duì)一條高速公路或是一個(gè)標(biāo)段進(jìn)行形變監(jiān)測(cè)與成因分析，針對(duì)大規(guī)模城市群高速公路網(wǎng)沿線地面沉降探測(cè)的研究不多，利用獲取的監(jiān)測(cè)成果進(jìn)行預(yù)測(cè)的更是少見(jiàn)。為此，本文采用SBAS-InSAR技術(shù)和BP神經(jīng)網(wǎng)絡(luò)模型相結(jié)合的方法，對(duì)滇中地區(qū)(包括昆明、楚雄、玉溪、曲靖)的6條高速公路、昆明主城區(qū)的一、二、三環(huán)3條環(huán)線和繞城高速沿線災(zāi)害進(jìn)行識(shí)別與預(yù)測(cè)。首先利用69景2018年7月20日～2020年11月18日的Sentinel-1A SAR影像，借助SBAS-InSAR技術(shù)獲取高速公路沿線地表沉降范圍和沉降值，探測(cè)沿線的災(zāi)害隱患點(diǎn)，再將其沉降值作為BP算法的訓(xùn)練樣本，建立已訓(xùn)練樣本和未來(lái)沉降值之間的函數(shù)模型；然后探討高速公路沿線識(shí)別出的21個(gè)潛在災(zāi)害點(diǎn)的沉降特征；最后從沉降速率場(chǎng)中，隨機(jī)提取 1 000個(gè)沉降點(diǎn)值，作為BP神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)樣本，建立InSAR監(jiān)測(cè)值和沉降預(yù)測(cè)值之間的函數(shù)模型，預(yù)測(cè)未來(lái)賓川斷陷盆地的地面沉降趨勢(shì)。在輸入學(xué)習(xí)樣本的訓(xùn)練刺激下，BP神經(jīng)網(wǎng)絡(luò)不斷改變網(wǎng)絡(luò)權(quán)值，逐漸使沉降預(yù)測(cè)值接近期望值[13]。該算法的學(xué)習(xí)能力較強(qiáng)，能夠處理復(fù)雜條件下的非線性形變問(wèn)題，被廣泛地應(yīng)用于沉降預(yù)測(cè)中，可為今后大范圍的道路管養(yǎng)和維護(hù)提供技術(shù)支持。

2 研究區(qū)與研究數(shù)據(jù)

2.1 研究區(qū)概況

滇中地區(qū)屬于滇東高原盆地，以山地和山間盆地地形為主，地勢(shì)起伏和緩。滇中地區(qū)包括：昆明、曲靖、楚雄和玉溪4個(gè)城市。滇中地區(qū)多盆地，集中了云南全省近一半的山間平地。區(qū)域位于長(zhǎng)江、珠江和紅河上游，屬亞熱帶氣候，日照充足，四季如春，氣候宜人，干濕季分明。土壤類型以紅壤為主。植被類型多樣，多為次生植被和人工植被。范圍總面積 94 558 km2，占全省土地面積的24%。滇中地區(qū)是云南省經(jīng)濟(jì)社會(huì)最發(fā)達(dá)、資源稟賦和基礎(chǔ)條件最優(yōu)越、區(qū)位優(yōu)勢(shì)最顯著的區(qū)域，是聯(lián)系東南亞和南亞各個(gè)國(guó)家的重要物資集散通道，是銜接印度洋、太平洋的陸上樞紐。本文選取昆明市主城區(qū)一環(huán)、二環(huán)、三環(huán)3條典型的城市道路，及機(jī)場(chǎng)高速、內(nèi)外繞城高速、汕昆高速(昆明-陸良段)、晉紅高速(昆明-晉寧-紅塔區(qū))、昆玉高速(昆明-玉溪)、昆曲高速(昆明-曲靖)、杭瑞高速(昆明-祿豐恐龍山段)8條滇中地區(qū)的高速公路作為研究對(duì)象。研究區(qū)地理位置及高速公路分布如圖1所示。

圖1 研究區(qū)地理位置及高速公路分布圖

2.2 研究數(shù)據(jù)

本文從歐空局官網(wǎng)下載了69景覆蓋研究區(qū)范圍的sentinel-1A SAR，對(duì)應(yīng)的精密軌道數(shù)據(jù)由ESA網(wǎng)站平臺(tái)獲取(https：//qc.sentinel1.eo.esa.int/)。SAR影像的具體參數(shù)如表1所示。為了去除地形效應(yīng)的影響，本文獲取了日本宇宙航空研究開(kāi)發(fā)機(jī)構(gòu)(JAXA)2015年5月免費(fèi)提供的全球數(shù)字地表模型數(shù)據(jù)，水平分辨率為 30 m(1弧秒)，高程精度 5 m。

表1 SAR影像的具體參數(shù)

3 研究方法

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

SBAS-InSAR技術(shù)是繼D-InSAR和PS-InSAR技術(shù)后，由Berardino和Lanari等人[14，15]提出的一種時(shí)間序列InSAR技術(shù)。該方法首先在獲取的多景SAR影像中挑選出一景作為超級(jí)主影像，在此基礎(chǔ)上，通過(guò)多主影像組合成一系列短基線干涉圖，有效地克服時(shí)間和空間失相干影響；然后根據(jù)干涉圖結(jié)果，選取沒(méi)有殘余地形條紋、形變條紋和相位躍變的GCP(Ground Control Points)點(diǎn)，對(duì)所有數(shù)據(jù)對(duì)進(jìn)行重去平和軌道精煉；最后，聯(lián)合最小二乘法和奇異分解法(Singular Value Decomposition，SVD)求解時(shí)序地表形變。

時(shí)序SBAS-InSAR技術(shù)獲取的地表形變?yōu)榈乇砣S形變投影到雷達(dá)視線方向(Line of Sight，LOS)的形變速率，而高速公路沿線地表沉降主要考慮的是垂直形變，因此，需要將LOS向的形變速率投影到垂向上。

LOS向的形變量[16]為：

DLOS=(UNsinφ-UEcosφ)sinθ+UVcosθ

(1)

式(1)中，UN、UE和UV表示南北、東西和垂直方向的地表形變分量，φ表示方位角，θ表示入射角?？紤]到InSAR技術(shù)對(duì)南北、東西方向的形變不敏感，不管是降軌還是升軌，垂直形變對(duì)LOS向形變的貢獻(xiàn)都在90%以上[17]，本文垂直方向上的形變計(jì)算公式為：

UV≈DLOS/cosθ

(2)

3.2 BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型建立

BP(Back Propagation)神經(jīng)網(wǎng)絡(luò)是由Rumelhart和McCelland等研究人員1986年提出的一種誤差后向傳播算法[18]。如圖2所示，BP模型由輸入層、隱含層、輸出層3層組成，其原理[19]為：

圖2 BP神經(jīng)網(wǎng)絡(luò)示意圖

設(shè)時(shí)間序列X(i)，i=1，2，…，N，N為數(shù)據(jù)個(gè)數(shù)，非線性的時(shí)間序列預(yù)測(cè)模型為：

X(t)=φ[X(t-1)，…，X(t-p)]

(3)

式中，φ[·]為非線性函數(shù)，p為模型的階數(shù)。由時(shí)間序列可以構(gòu)造(N-p)個(gè)樣本，樣本號(hào)輸入層和輸出層為：

神經(jīng)網(wǎng)絡(luò)建模的本質(zhì)就是用神經(jīng)網(wǎng)絡(luò)的隱式來(lái)表達(dá)輸入與輸出的函數(shù)關(guān)系，將所構(gòu)造樣本集代入網(wǎng)絡(luò)進(jìn)行訓(xùn)練，建立從輸入到輸出之間的非線性的映射關(guān)系，并將“知識(shí)信息”儲(chǔ)存在連接權(quán)上。因而，可用于非線性分類、預(yù)測(cè)等問(wèn)題，作為函數(shù)計(jì)算器，能以任意精度逼近任何非線性函數(shù)[20]。

誤差后向傳播過(guò)程是通過(guò)使能量函數(shù)最小化來(lái)完成輸入輸出的映射。將能量函數(shù)E定義為輸出層單元的均方誤差，設(shè)有(N-p)個(gè)學(xué)習(xí)樣本，則：

(4)

式中，x(t)，yi分別為樣本i的實(shí)際輸出和希望輸出。

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

4.1 滇中地區(qū)高速公路沿線沉降分布特征

本文使用SARscape軟件，對(duì)研究區(qū)內(nèi)的69景Sentinel-1A SAR影像進(jìn)行處理，為了更清晰、準(zhǔn)確地探討道路沿線地表形變對(duì)行車安全的影響，本文在城區(qū)道路(一環(huán)、二環(huán)) 200 m范圍內(nèi)，三環(huán)及其余各高速公路 1 000 m范圍內(nèi)建立緩沖區(qū)，提取道路沿線的形變特征及時(shí)空分布，如圖3所示。圖3中，形變速率值小于0的區(qū)域(紅色)，表示遠(yuǎn)離衛(wèi)星運(yùn)動(dòng)，實(shí)際地表表現(xiàn)為沉降；形變速率大于0的區(qū)域(綠色)，表示朝著衛(wèi)星運(yùn)動(dòng)，實(shí)際地表表現(xiàn)為抬升。研究區(qū)內(nèi)3條城市主干道及8條高速公路分布范圍內(nèi)，2018年7月～2020年11月期間，累積沉降速率為 -64.54 mm/a～+44.68 mm/a，地表形變廣泛分布，共探測(cè)出22處明顯沉降區(qū)域，分別編號(hào)為1～22。

圖3 研究區(qū)形變速率

為了更清晰、準(zhǔn)確地描述SBAS-InSAR技術(shù)探測(cè)識(shí)別的明顯沉降點(diǎn)，將22個(gè)沉降漏斗的沉降特征及危險(xiǎn)性統(tǒng)計(jì)如表2所示。滇中地區(qū)的高速網(wǎng)絡(luò)明顯沉降點(diǎn)最多的為昆明繞城高速公路，一共識(shí)別出6個(gè)沉降點(diǎn)，編號(hào)分別為6、7、11、12、15、16；沉降最嚴(yán)重的為三環(huán)線的4號(hào)沉降點(diǎn)，最大沉降速率為 -76.50 mm/a，位于靠近滇池附近的草海隧道一帶，涉及草海左岸的董家村、小圍村，右岸從西向東沿線的新河村、河尾村、楊家塘、小河尾、太河村等地，沉降面積約為 6.7 km2；研究區(qū)內(nèi)較為穩(wěn)定的高速公路為杭瑞高速昆明至楚雄段，除跨越昆明主城區(qū)的3號(hào)沉降點(diǎn)外，基本處于穩(wěn)定狀態(tài)。本文按照探測(cè)的22個(gè)沉降點(diǎn)的沉降速率，采用自然斷點(diǎn)法，將其對(duì)高速公路運(yùn)營(yíng)的可能影響程度分為輕微、一般、嚴(yán)重三個(gè)等級(jí)，由表2可知，可能致災(zāi)的沉降點(diǎn)為12個(gè)，占識(shí)別總個(gè)數(shù)的54.5%；影響一般的為8個(gè)，占比為36.4%；輕微的僅為2個(gè)，占比為9.1%。

表2 滇中地區(qū)高速公路網(wǎng)沿線沉降情況

續(xù)表2

4.2 滇中地區(qū)高速公路網(wǎng)絡(luò)地表沉降預(yù)測(cè)

從上述的監(jiān)測(cè)結(jié)果可以知道，滇中地區(qū)高速網(wǎng)絡(luò)不同路段均出現(xiàn)了不均勻的沉降現(xiàn)象，涉及的沉降區(qū)域，范圍大、零星分布。為了縮小監(jiān)測(cè)范圍，準(zhǔn)確、及時(shí)地進(jìn)行管養(yǎng)，本文在沉降監(jiān)測(cè)的基礎(chǔ)上，隨機(jī)選取研究區(qū)的高速公路網(wǎng)沿線緩沖區(qū)范圍內(nèi)的沉降監(jiān)測(cè)點(diǎn)，利用BP神經(jīng)網(wǎng)絡(luò)模型對(duì)未來(lái)的沉降趨勢(shì)進(jìn)行預(yù)測(cè)分析。試驗(yàn)中，從研究區(qū)內(nèi)隨機(jī)抽取 1 000組InSAR監(jiān)測(cè)數(shù)據(jù)，選出980組數(shù)據(jù)作為計(jì)算樣本，在980組數(shù)據(jù)中，選取80%作為訓(xùn)練樣本，20%作為測(cè)試樣本，得到滇中地區(qū)高速公路沉降趨勢(shì)預(yù)測(cè)的BP神經(jīng)網(wǎng)絡(luò)模型，利用剩余的20組數(shù)據(jù)對(duì)該模型進(jìn)行驗(yàn)證。根據(jù)選取的數(shù)據(jù)，利用Matlab軟件構(gòu)建BP神經(jīng)網(wǎng)絡(luò)模型訓(xùn)練時(shí)，設(shè)置的參數(shù)如下：學(xué)習(xí)速率0.03，訓(xùn)練目標(biāo)0.001、訓(xùn)練次數(shù) 20 000。驗(yàn)證比對(duì)結(jié)果和最佳迭代次數(shù)如圖4所示。圖4(a)中預(yù)測(cè)輸出值用藍(lán)色星號(hào)表示，期望輸出值(SBAS監(jiān)測(cè)值)用紅色圓圈表示；圖4(b)中橫坐標(biāo)表示BP神經(jīng)網(wǎng)絡(luò)訓(xùn)練時(shí)的迭代收斂次數(shù)，縱坐標(biāo)表示均方誤差。實(shí)驗(yàn)結(jié)果表明，BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)的均方誤差為 0.167 mm，平均絕對(duì)誤差為 0.379 mm。進(jìn)一步統(tǒng)計(jì)BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)的精度，如表3所示，SBAS監(jiān)測(cè)值與BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)值之間的最大誤差為 0.666 mm/a。說(shuō)明BP神經(jīng)網(wǎng)絡(luò)模型對(duì)滇中地區(qū)高速公路地表沉降值預(yù)測(cè)效果較好，可為今后的道路災(zāi)害識(shí)別普查提供必要的技術(shù)手段支持。

圖4 BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)試驗(yàn)結(jié)果

表3 BP神經(jīng)網(wǎng)絡(luò)模型預(yù)測(cè)精度統(tǒng)計(jì)表

5 結(jié) 語(yǔ)

本文運(yùn)用SBAS-InSAR技術(shù)和SARScape軟件處理68景滇中地區(qū)的Sentinel-1A影像，獲取昆明主城區(qū)道路(一環(huán)、二環(huán))200 m范圍內(nèi)，三環(huán)及滇中地區(qū)高速公路網(wǎng) 1 000 m范圍內(nèi)，道路沿線的地表形變特征及時(shí)空分布，對(duì)研究范圍內(nèi)的地表沉降特征進(jìn)行分析；同時(shí)，利用BP神經(jīng)網(wǎng)絡(luò)模型，預(yù)測(cè)研究區(qū)內(nèi)的沉降趨勢(shì)，得出以下結(jié)論：

(1)研究區(qū)內(nèi)3條城市主干道及8條高速公路分布范圍內(nèi)，2018年7月～2020年11月期間，累積沉降速率為 -76.50 mm/a～+44.68 mm/a，地表形變廣泛分布，共探測(cè)出22個(gè)災(zāi)害隱患點(diǎn)；

(2)BP神經(jīng)網(wǎng)絡(luò)模型的學(xué)習(xí)能力和容錯(cuò)能力較強(qiáng)，對(duì)滇中地區(qū)高速公路地表沉降值預(yù)測(cè)效果接近InSAR監(jiān)測(cè)值，可為今后的滇中地區(qū)高速公路致災(zāi)范圍圈定提供依據(jù)。