基于SBAS技術(shù)的輸電線路地表形變監(jiān)測(cè)

馬 儀，耿 浩，黃 然，張 輝，馬御棠，劉 靖

(1.電力遙感技術(shù)聯(lián)合實(shí)驗(yàn)室(云南電網(wǎng)有限責(zé)任公司電力科學(xué)研究院)，昆明 650217；2.云南電網(wǎng)有限責(zé)任公司，昆明 650000；3.北京空間飛行器總體設(shè)計(jì)部，北京 100094)

云南地處我國(guó)西南邊陲，位于亞歐板塊與印度洋板塊的交界處，地形復(fù)雜，西北部為深谷的橫斷山區(qū)，東部和南部為云貴高原。云南緯度低，海拔高，受季風(fēng)影響，形成了年溫差小、日溫差大、河流縱橫交錯(cuò)的地形地貌[1]，省內(nèi)地質(zhì)災(zāi)害頻發(fā)，僅2021年4月發(fā)生地震2009次，按震級(jí)統(tǒng)計(jì)，其中2.0～2.9級(jí)20次、3.0～3.9級(jí)5次。云南省主要遭受的地質(zhì)災(zāi)害為滑坡、泥石流、崩塌、地面塌陷等[2]。地質(zhì)災(zāi)害會(huì)對(duì)該省的輸電系統(tǒng)造成安全隱患。

輸電線路覆蓋區(qū)域廣闊，沿途地理環(huán)境復(fù)雜，容易受到地質(zhì)災(zāi)害的影響，輸電桿塔及輸電網(wǎng)長(zhǎng)期在這種環(huán)境中會(huì)受到擾動(dòng)，主要影響包括地質(zhì)災(zāi)害造成的地表劇烈形變或地表緩慢形變引起的桿塔本體形變。形變會(huì)對(duì)電網(wǎng)安全產(chǎn)生影響，例如引發(fā)桿塔傾斜、倒塌等事故。高成本的人工巡查可以發(fā)現(xiàn)形變明顯的桿塔，但不能很好地監(jiān)測(cè)到桿塔地基的緩慢形變，且人工巡查效率較低，巡檢周期長(zhǎng)，不能很好地達(dá)到巡檢目的。久而久之，桿塔區(qū)域地表的形變導(dǎo)致桿塔形變或傾斜，因未及時(shí)發(fā)現(xiàn)而形成安全隱患，對(duì)電網(wǎng)系統(tǒng)造成破壞，會(huì)導(dǎo)致大面積斷電，給社會(huì)帶來(lái)經(jīng)濟(jì)損失[3]。因此，對(duì)輸電網(wǎng)絡(luò)沿線地表沉降的周期性監(jiān)測(cè)是十分必要的，有效的監(jiān)測(cè)能及時(shí)發(fā)現(xiàn)輸電網(wǎng)絡(luò)中存在的安全隱患，及時(shí)對(duì)隱患區(qū)域做出處理，減少地質(zhì)災(zāi)害給電網(wǎng)系統(tǒng)帶來(lái)的威脅。

隨著遙感技術(shù)的發(fā)展，SAR衛(wèi)星數(shù)據(jù)在地質(zhì)方面的研究越來(lái)越廣泛，關(guān)于SAR衛(wèi)星測(cè)量地表形變的技術(shù)越來(lái)越成熟，將SAR衛(wèi)星數(shù)據(jù)測(cè)量技術(shù)應(yīng)用于電網(wǎng)線路的監(jiān)測(cè)，能夠?qū)崿F(xiàn)輸電線路沿線地表形變的監(jiān)測(cè)，降低各類地質(zhì)災(zāi)害引起的輸電網(wǎng)絡(luò)故障發(fā)生率。衛(wèi)星運(yùn)用合成孔徑雷達(dá)(SAR)對(duì)觀測(cè)區(qū)域發(fā)射微波，并接收后向散射回波完成對(duì)地成像，無(wú)需依靠太陽(yáng)輻射，全天候觀測(cè)成像，對(duì)雨雪穿透力強(qiáng)[4]。將兩幅SAR影像進(jìn)行干涉處理，即干涉雷達(dá)測(cè)量技術(shù)(Interferometric SAR, InSAR)，可用于對(duì)地表形變的測(cè)量。InSAR技術(shù)有著高精度、高分辨率、短周期、覆蓋范圍廣、全天候等特點(diǎn)，在形變監(jiān)測(cè)中應(yīng)用廣泛[5]。在InSAR技術(shù)的運(yùn)用過程中，研究人員發(fā)現(xiàn)受時(shí)間去相干和空間去相干的影響，無(wú)法獲得形變隨時(shí)間的演化過程，于是Berandino等提出了短基線集(SBAS-InSAR)技術(shù)(SBAS)[6]。SBAS能夠有效地獲取長(zhǎng)時(shí)序的地表形變信息，且形變測(cè)量結(jié)果精度高，達(dá)到mm級(jí)，還能消除大氣效應(yīng)和失相干的影響。綜上，文中利用SBAS技術(shù)，基于Sentinel-1A SAR衛(wèi)星影像，對(duì)影像覆蓋區(qū)域的3條輸電線路沿線區(qū)域進(jìn)行地表形變監(jiān)測(cè)，并對(duì)該區(qū)域存在隱患的輸電桿塔加以分析，給出維護(hù)建議。

1 SBAS的技術(shù)原理及其數(shù)據(jù)處理流程

1.1 InSAR技術(shù)的發(fā)展

1969年Rogers等首次應(yīng)用InSAR技術(shù)獲取金星的地面高程[7]；1980年Zebker等利用InSAR技術(shù)提取了地面的高程信息并應(yīng)用于監(jiān)測(cè)地表形變[8],1993年Massonnent等利用InSAR技術(shù)處理 ERS-1數(shù)據(jù)獲取了1992年美國(guó)蘭德斯地震的形變場(chǎng)[9-10]，此后，InSAR獲取形變技術(shù)得到了全世界的關(guān)注，被廣泛用于地質(zhì)災(zāi)害監(jiān)測(cè)，例如地震、火山活動(dòng)、滑坡、凍土等。在InSAR的基礎(chǔ)上，F(xiàn)erretti提出永久性散射差分干涉測(cè)量(PS-InSAR)技術(shù)(PS)，用于監(jiān)測(cè)意大利安科納地區(qū)的滑坡[11]；2002年，Berardino提出了短基線差分干涉測(cè)量SBAS技術(shù)，并用于意大利的Campo Flegrei火山口的時(shí)序地表形變監(jiān)測(cè)，取得了理想的地表形變監(jiān)測(cè)結(jié)果[6]。2004年，Trasatti E基于ENVISAT數(shù)據(jù)，利用SBAS-DInSAR技術(shù)成功監(jiān)測(cè)了Campi Flegrei火山口抬升現(xiàn)象[12]。PS和SBAS都是以消除InSAR技術(shù)處理過程中遇到的大氣干擾、時(shí)間去相干和空間去相干等問題提出來(lái)的處理方法，并且都可以有效地抑制去相干和大氣延遲效應(yīng)以及噪聲。

PS能夠獲取大范圍的高精度形變速率，該技術(shù)要求25景以上的SAR影像[13]才可獲得可靠的測(cè)量結(jié)果，且該方法通常要求具有穩(wěn)定散射特性的點(diǎn)，采用差值內(nèi)插出整個(gè)監(jiān)測(cè)區(qū)域的地表形變，犧牲空間點(diǎn)的密度來(lái)獲取地表長(zhǎng)時(shí)間的形變序列。SBAS技術(shù)可通過少量的SAR影像完成對(duì)地表形變的監(jiān)測(cè)，結(jié)果可靠，該方法以增加時(shí)間采樣率來(lái)提高時(shí)間相干性和空間相干性，相比于PS技術(shù)減少了對(duì)高相干點(diǎn)數(shù)量和空間分布的約束條件，因此SBAS技術(shù)在滿足精度要求的條件下，更適用于山區(qū)、礦區(qū)、流域、滑坡等復(fù)雜地形或地質(zhì)災(zāi)害的形變監(jiān)測(cè)。孫偉等人利用SBAS技術(shù)，結(jié)合遙感影像、降雨數(shù)據(jù)和坡度坡向數(shù)據(jù)對(duì)萬(wàn)安縣輸油管道重點(diǎn)區(qū)段進(jìn)行形變監(jiān)測(cè)，結(jié)果表明降雨天數(shù)越多累積形變?cè)酱蟮慕Y(jié)論[14]；戴可人等人利用SBAS技術(shù)對(duì)雅礱江流域雅江縣-木里縣段的高山峽谷區(qū)域進(jìn)行滑坡災(zāi)害的識(shí)別和預(yù)警，并成功找到了8處隱患區(qū)域，有效減少了滑坡災(zāi)害的發(fā)生[15]；那靜則對(duì)迪慶州北部區(qū)域進(jìn)行時(shí)間跨度兩年的地表形變周期性監(jiān)測(cè)，通過監(jiān)測(cè)結(jié)果對(duì)災(zāi)害易發(fā)區(qū)域進(jìn)行排查，為后續(xù)形變監(jiān)測(cè)研究和地質(zhì)災(zāi)害防治提供技術(shù)參考[16]；馮文凱等人同樣利用SBAS技術(shù)將金沙江流域沃達(dá)村巨型老滑坡進(jìn)行劃定，分成了強(qiáng)烈形變區(qū)和均勻形變區(qū)，結(jié)果表明SBAS技術(shù)在復(fù)雜山區(qū)地質(zhì)災(zāi)害監(jiān)測(cè)預(yù)警領(lǐng)域與其他監(jiān)測(cè)方式具有較好的重合度，也印證了SBAS技術(shù)監(jiān)測(cè)的準(zhǔn)確性[17]。綜上所述，文中采用SBAS技術(shù)對(duì)輸電桿塔沿線進(jìn)行形變監(jiān)測(cè)。

1.2 SBAS的基本原理

SAR衛(wèi)星影像是復(fù)數(shù)影像，影像的像元記錄地面后向散射信息，后向散射信息包括每一地面單元的后向散射強(qiáng)度以及雷達(dá)天線距離相關(guān)的相位信息，像元的散射強(qiáng)度與相位共用一個(gè)復(fù)數(shù)表示。通過對(duì)SAR影像進(jìn)行干涉處理，可以得到反映地表起伏和形變的干涉相位條紋，再借助復(fù)數(shù)影像對(duì)的成像幾何關(guān)系和已知的控制點(diǎn)信息，即可獲得每個(gè)像元的高程或形變量[18-19]。

SBAS是對(duì)InSAR差分處理的形變值進(jìn)行時(shí)序分析的技術(shù)，該方法將SAR影像依據(jù)時(shí)間基線和空間基線閾值分成若干短基線子集，生成差分干涉圖，提高相干性[20]。該技術(shù)的基本算法是將SAR影像數(shù)據(jù)分成若干小集合，然后根據(jù)干涉條件等設(shè)置時(shí)間基線和空間基線的閾值，對(duì)每個(gè)小集合利用最小二乘法求解小集合內(nèi)的時(shí)間形變序列，再聯(lián)合各個(gè)小集合進(jìn)行奇異值分解(SVD)，反演出整個(gè)地表在影像時(shí)間段內(nèi)完整沉降序列。使用SVD的方法可以將分開的小基線之間聯(lián)合起來(lái)求解，也可以通過獲取時(shí)間和空間信息分離大氣延遲相位，小基線分析了SAR影像間所有可用的干涉圖，大大增加了形變差分干涉圖的空間密度和時(shí)間采樣率[19,21-27]。

假設(shè)有同一區(qū)域的N+1幅SAR影像按時(shí)間t排序，任選一幅SAR影像作為主影像，然后將其他N幅SAR影像在同一坐標(biāo)系下與主影像配準(zhǔn)，并進(jìn)行多視差分處理得到M幅多視后的差分干涉圖，干涉圖的時(shí)間基線與空間基線均小于設(shè)定的閾值，并且M滿足

(1)

對(duì)干涉圖進(jìn)行濾波以及解纏處理，并在解纏后的干涉圖中篩選相干性較好的點(diǎn)，對(duì)這些點(diǎn)用式(2)計(jì)算。假設(shè)第i幅解纏后的差分干涉圖分別由主圖像和輔影像tA和tB時(shí)刻獲得的SAR影像干涉差分生成的，那么距離向坐標(biāo)為r和方位向坐標(biāo)為x的高相干性相干點(diǎn)的像元干涉相位可以寫成：

δφi(x,r)=φA(x,r)-φB(x,r)≈

(2)

忽略地形殘余相位、大氣延遲相位以及噪聲相位，式(2)簡(jiǎn)化為：

δφi(x,r)=φA(x,r)-φB(x,r)≈

(3)

待求的像元點(diǎn)的形變量所對(duì)應(yīng)的N個(gè)未知相位值可表示為：

ΦT=[φ(t1),φ(t2),…,φ(tN)].

(4)

通過M個(gè)差分干涉圖得出的干涉相位表示為：

ΔΦT=[δφ1,δφ2,…,δφM].

(5)

主影像和從影像的時(shí)間序列分別為：

IE=[IE1,…IEJ,…,IEM],

(6)

IS=[IS1,…IEJ,…,ISM].

(7)

第j幅差分干涉圖對(duì)應(yīng)的干涉相位為：

δφ1=φ(tIEj)-φ(tISj),j=1,2,…,M.

(8)

將所有干涉圖的干涉相位用方程表示，并擴(kuò)展為矩陣形式：

AΦ=ΔΦ.

(9)

A是M×N的一個(gè)矩陣，每一行對(duì)應(yīng)一幅干涉圖，每一列對(duì)應(yīng)一個(gè)時(shí)間序列的SAR影像。

(10)

矩陣中-1表示從影像所在列，+1表示主影像所在列。對(duì)于式(9)，當(dāng)矩陣Φ為N×1的矩陣時(shí)，M≥N，A的秩為N，根據(jù)最小二乘法求得形變相位：

Φ=(ATA)-1ATΔΦ.

(11)

在實(shí)際中，SAR影像數(shù)據(jù)往往分為不同的子集矩陣，因此ATA的秩小于列數(shù)N，即ATA為奇異矩陣，這使得式(11)有無(wú)窮的解。為求得合適的解，需通過SVD方法求得矩陣A的廣義逆矩陣，大致步驟為利用SVD分解矩陣A，在最小二乘約束下，求解相位，對(duì)相位在時(shí)間上量化獲得形變速率，用于形變量的計(jì)算和形變分析。

1.3 處理流程

SBAS技術(shù)用于形變量估計(jì)的步驟大致如下：選定主影像和從影像，設(shè)定空間基線閾值，對(duì)配準(zhǔn)完成的SAR數(shù)據(jù)干涉處理，生成M幅干涉圖，計(jì)算相干系數(shù)；根據(jù)相干系數(shù)，選擇高相干點(diǎn)作為相位解纏參考點(diǎn)；利用干涉圖基線和外部DEM進(jìn)行差分干涉處理，去除平地效應(yīng)和地形相位的影響，再對(duì)差分干涉圖進(jìn)行濾波；利用SVD和最小二乘法解算小基線幾何，獲得大氣延遲差和地形相位；最后去噪，獲取時(shí)間序列地表形變信息。

具體處理流程如圖1所示。

圖1 SBAS技術(shù)處理流程

2 輸電桿塔線路形變監(jiān)測(cè)

2.1 數(shù)據(jù)介紹

文中的SAR衛(wèi)星數(shù)據(jù)來(lái)源于哨兵1號(hào) A(Sentinel-1A)衛(wèi)星。哨兵一號(hào)為歐洲航天局研制的SAR衛(wèi)星，通過C波段對(duì)地全天時(shí)全天候觀測(cè)，該衛(wèi)星重訪周期為12 d，具有雙極化和重訪周期短的優(yōu)點(diǎn)，可實(shí)現(xiàn)全球高分辨率監(jiān)測(cè)和對(duì)同一地區(qū)的長(zhǎng)時(shí)間序列監(jiān)測(cè)。哨兵1A具有4種成像模式，每種成像模式分辨率與幅寬不一致。文中SAR衛(wèi)星數(shù)據(jù)覆蓋的3條輸電線，即位于云南省楚雄彝族自治州元謀縣北部的RT甲線、RT乙線、YF直流線，針對(duì)3條輸電線路的分布和沿線的形變規(guī)律，選取40景升軌Sentinel-1A影像作為處理數(shù)據(jù)。區(qū)域中心坐標(biāo)為101°46′E，25°50′N，影像覆蓋面積為1 000 km2，影像采集時(shí)間范圍為2018年11月3日至2020年5月8日。該組影像成像模式為IW模式，分辨率為5 m×12 m，采用VV+VH極化方式，數(shù)據(jù)的成像日期間隔12 d/景，數(shù)據(jù)覆蓋范圍如圖2所示。

圖2 Sentinel-1A影像覆蓋范圍

2.2 輸電線路周邊地質(zhì)沉降分析

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

對(duì)40景SAR影像干涉處理，得到的干涉對(duì)時(shí)空基線如圖3所示。由圖知，數(shù)據(jù)中未生成連續(xù)的干涉對(duì)，因此采用SBAS技術(shù)進(jìn)行年平均形變速率反演。經(jīng)InSAR技術(shù)處理得到的形變速率是沿雷達(dá)視線方向上的速率，為更加直觀地看到地面點(diǎn)在垂直方向上的變化情況，將數(shù)據(jù)平均形變速率投影到垂直地面方向，結(jié)果如圖4所示。數(shù)據(jù)結(jié)果表明，影像中大部分區(qū)域呈現(xiàn)出±40 mm/a的形變速率波動(dòng)，形變速率為正的代表地表抬升，形變速率為負(fù)表示地表下沉形變。

圖3 40景Sentinel-1A數(shù)據(jù)時(shí)空基線

圖4 數(shù)據(jù)覆蓋范圍內(nèi)年平均形變速率

2.2.2 基于桿塔線路的形變分析

數(shù)據(jù)范圍內(nèi)存在3條不同的桿塔線路，分別為RT甲線、RT乙線和YF直流線，線路分布情況如圖4所示，RT甲線和RT乙線分布一致，為東西走向，YF直流線為南北走向。分別對(duì)3條線路進(jìn)行沿線形變速率分析，結(jié)果如圖5所示。

圖5顯示，RT甲線、RT乙線和YF直流線3條線路上的點(diǎn)在形變速率上均存在不同程度的波動(dòng)，其中，RT甲線的波動(dòng)幅度較小，其波動(dòng)范圍在±40 mm/a之間，最大抬升速率與最大下沉速率的差值最大為80 mm/a；RT乙線沿線上的形變速率極差最大，該線大部分點(diǎn)的形變速率波動(dòng)范圍在±30 mm/a之間，沿線最大下沉速率點(diǎn)超過了-60 mm/a，最大抬升速率點(diǎn)大概為50 mm/a；YF直流線的波動(dòng)幅度受到其線路上出現(xiàn)的沉降漏斗影響，形變速率極差達(dá)到了90 mm/a。數(shù)據(jù)結(jié)果顯示，RT甲線和RT乙線在沿線上沒有明顯的下沉量逐漸變大或變小的趨勢(shì)，YF直流線則呈現(xiàn)下沉速率逐漸變小的趨勢(shì)，綜上所述，RT甲線、RT乙線和YF直流線3條線路沿線地表狀態(tài)整體較為穩(wěn)定。

3 重點(diǎn)桿塔區(qū)域形變分析

3條輸電線路桿塔沿線均呈現(xiàn)出較為穩(wěn)定的波動(dòng)，沉降速率大多不超過±40 mm/a。而圖5(b)、圖5(c)顯示，在RT乙線和YF直流線上存在明顯的沉降漏斗，為更好地了解每條線路上出現(xiàn)的形變異常點(diǎn)和區(qū)域內(nèi)的時(shí)間序列形變情況，文中對(duì)RT乙線和YF直流線的部分區(qū)域重點(diǎn)分析。

3.1 RT乙線重點(diǎn)觀察區(qū)域

RT乙線桿塔沿線大部分形變點(diǎn)速率處于-40～40 mm/a之間，局部監(jiān)測(cè)到“漏斗”式的形變速率變化段。根據(jù)圖5(b)知，漏斗靠近#84桿塔，因此對(duì)#83～#85桿塔線路進(jìn)一步分析。將#83～#85桿塔線路視為剖面線，分析剖面線年平均形變速率變化情況，如圖6所示，剖面線上的形變速率變化情況如圖7所示，在#84桿塔東側(cè)出現(xiàn)明顯的沉降漏斗，最大下沉速率達(dá)到了-70 mm/a。#84處于沉降漏斗的邊緣，容易受到沉降漏斗面積逐漸擴(kuò)大的的影響。

圖5 輸電線路沿線形變速率折線圖

圖6 RT乙線重點(diǎn)觀察區(qū)域干涉影像

圖7 RT乙線#83～#85形變速率

對(duì)#84桿塔和#85桿塔所在位置和沉降漏斗中心點(diǎn)(25°50′36″N，101°49′19.76″E)作為采樣點(diǎn)分別進(jìn)行時(shí)序分析，結(jié)果如圖8所示，結(jié)果表明，#84桿塔形變速率為負(fù)，處于沉降趨勢(shì)，#85形變速率為正，為抬升趨勢(shì)。與沉降中心對(duì)比，漏斗中心下沉速率明顯，兩桿塔的形變速率則顯得平穩(wěn)。2018年11月3日至2020年5月8日期間，漏斗中心累計(jì)下沉量達(dá)到了110 mm，且呈現(xiàn)繼續(xù)下沉的趨勢(shì)，下沉速率趨向于增大的狀態(tài)。

圖8 #84、#85桿塔及下沉中心時(shí)序形變分析圖

3.2 YF直流線重點(diǎn)觀察區(qū)域

圖5(c)顯示，沿桿塔線路由北(#51)至南(#70)距#51桿塔2.1～2.3 km范圍內(nèi)形變速率驟降，該區(qū)域形變速率與周邊地區(qū)速率明顯不同。沉降漏斗區(qū)域位置如圖9所示，圖中顯示，下沉區(qū)域位于桿塔#55～#56之間，為西南-東北的走向，此外，在#55～#58桿塔間存在離散的下沉點(diǎn)。對(duì)形變區(qū)域進(jìn)一部研究，分析形變情況，將#55～#58桿塔線路作為剖面線，如圖10所示。

圖9 YF直流線重點(diǎn)觀察區(qū)域干涉圖

圖10 YF直流線#55～#58形變速率

圖10顯示，4座桿塔之間存在3個(gè)明顯的沉降漏斗，其中，#55～#57之間的形變速率處于-15～-45 mm/a之間，兩線桿間形變采樣點(diǎn)最大與最小形變速率差值小于30 mm/a，最大下沉速率為-45 mm/a，兩桿塔之間的形變速率變化較小，下沉速率不大，故將#55～#57桿塔線路視為穩(wěn)定變化的線路。#58桿塔位于沉降漏斗的邊緣部分，漏斗區(qū)域最大下沉速率大于-60 mm/a。

對(duì)#58桿塔和沉降漏斗中心(25°56′23.47″N，101°47′35.07°E)進(jìn)行時(shí)間序列形變分析，兩點(diǎn)形變速率變化如圖11所示。2019年6月之前，兩點(diǎn)的形變趨勢(shì)及累計(jì)形變量未出現(xiàn)較大差異；自2019年6月起，#58桿塔開始趨于平穩(wěn)，而下沉中心仍然在緩慢下沉，且截至2020年5月8日，該中心點(diǎn)仍存在下沉趨勢(shì)，考慮到逐漸下沉帶來(lái)的累計(jì)形變?cè)龃?、形變差異的擴(kuò)大，需將#58桿塔列為重點(diǎn)觀察對(duì)象，繼續(xù)觀察該桿塔的形變規(guī)律和周邊區(qū)域累計(jì)沉降量的變化。

圖11 YF直流線#58和下沉中心時(shí)間序列形變分析圖

3.3 區(qū)域形變分析

不同方向上累計(jì)形變差異會(huì)致使桿塔塔基在不同方向上出現(xiàn)形變，而隨著形變差異的不斷擴(kuò)大，桿塔會(huì)出現(xiàn)傾斜、桿塔區(qū)域出現(xiàn)地面裂縫等影響桿塔穩(wěn)定性及壽命的情況。出現(xiàn)形變嚴(yán)重的“沉降漏斗”時(shí)，漏斗范圍內(nèi)及范圍外1 km區(qū)域內(nèi)的桿塔均作為重點(diǎn)觀察的對(duì)象，以防范和應(yīng)對(duì)地面下沉和抬升引起的裂縫、地面塌陷、地面沉降等不同程度的地質(zhì)災(zāi)害給桿塔帶來(lái)的影響。

4 結(jié) 論

為監(jiān)測(cè)輸電桿塔沿線區(qū)域的地表形變，消除地質(zhì)災(zāi)害帶來(lái)的輸電安全隱患，確保輸電系統(tǒng)的正常運(yùn)作，文中提出運(yùn)用SBAS技術(shù)監(jiān)測(cè)3條輸電桿塔沿線形變情況。將40景Sentinel-1A數(shù)據(jù)運(yùn)用SBAS技術(shù)處理，得出數(shù)據(jù)覆蓋范圍內(nèi)的3條桿塔線路在2018年11月至2020年5月期間的年平均形變速率和時(shí)間序列形變情況。RT甲線、RT乙線和YF直流線均呈現(xiàn)為穩(wěn)定的形變規(guī)律，多在±40 mm/a之間，此外，YF直流線沿線的年平均形變速率呈現(xiàn)下沉速率逐漸減小的情況。

鑒于RT甲線波動(dòng)幅度小的特點(diǎn)，僅在RT乙線以及YF直流線中挑選出重點(diǎn)分析區(qū)域進(jìn)一步分析。結(jié)果表明，RT乙線上的#84和#85桿塔之間存在沉降漏斗，且#84 桿塔位于漏斗邊緣地區(qū)，易受到該漏斗的影響。在YF直流線沿線的#55～#58 線路上，存在大小程度不一的沉降漏斗，目前影響范圍較小，但形變速率為-70 mm/a，且靠近#58桿塔，因此#58號(hào)桿塔容易受到該沉降漏斗的影響。綜上，需對(duì)RT乙線#84桿塔和YF直流線#58桿塔加強(qiáng)監(jiān)測(cè)，避免因地表形變引發(fā)桿塔形變隱患。

通過運(yùn)用SBAS技術(shù)，有效地完成對(duì)3條線路輸電桿塔沿線的地表形變監(jiān)測(cè)，并對(duì)監(jiān)測(cè)結(jié)果進(jìn)行分析，給出輸電桿塔在地表形變上的維護(hù)建議。下一步考慮擴(kuò)大SBAS用于輸電桿塔沿線監(jiān)測(cè)范圍，普及到其他各個(gè)輸電線路，降低維護(hù)成本，提高監(jiān)測(cè)效率。