基于Sentinel-1數(shù)據(jù)和SBAS-InSAR技術(shù)的大藤峽水利樞紐變形監(jiān)測與分析

馮瀚 汪娟 嚴(yán)浩元 楊雋 向淇文

變形監(jiān)測是大壩安全監(jiān)測的主要內(nèi)容之一，也是大壩健康穩(wěn)定運行的主要依據(jù)。針對大藤峽水利樞紐設(shè)施的變形監(jiān)測，基于2020年2月至2023年1月的Sentinel-1數(shù)據(jù)，利用星載時序SBAS-InSAR方法，開展了大藤峽大壩及其周邊區(qū)域的形變監(jiān)測。由變形監(jiān)測結(jié)果可知，大藤峽左壩主體部分保持穩(wěn)定，右岸有1處斜坡累積變形較大，水庫上游兩岸存在4處變形斜坡。應(yīng)用研究表明，Sentinel-1數(shù)據(jù)和SBAS-InSAR技術(shù)可以實現(xiàn)大藤峽大壩低成本、大面積的時序變形監(jiān)測，為大壩安全監(jiān)測預(yù)警提供一種有效的技術(shù)手段，在水庫大壩形變監(jiān)測領(lǐng)域有著較為廣闊的應(yīng)用前景。

水利樞紐工程在防洪、蓄水、灌溉、航運、發(fā)電等方面發(fā)揮著巨大的社會效益和經(jīng)濟效益，水利樞紐的安全直接關(guān)系到人民群眾的生命財產(chǎn)安全，進行變形監(jiān)測是大壩安全的重要保障。針對大壩的變形監(jiān)測，傳統(tǒng)監(jiān)測技術(shù)主要采用全站儀、水準(zhǔn)儀和GPS等獲取大壩形變信息，需要耗費大量的人力、物力和財力，且只能獲取局部形變信息，無法實現(xiàn)大壩全方位高密度的變形監(jiān)測。星載合成孔徑雷達干涉測量（Interferometric Synthetic Aperture Radar， InSAR）技術(shù)具有全天時、全天候、高精度、大范圍獲取地面變形信息的優(yōu)勢，為水庫大壩的安全監(jiān)測提供了一種新型工具和監(jiān)測手段。

國內(nèi)外已有學(xué)者利用InSAR技術(shù)進行大壩的變形監(jiān)測研究，如王騰等人（2011）采用40幅SAR影像，結(jié)合PS-InSAR與QPS-InSAR方法獲取了三峽大壩的時序變形信息，詳細分析了混凝土重力壩在建設(shè)及運行過程中的穩(wěn)定性及形變規(guī)律；Pietro Milillo，Daniele Perissin等人（2015）利用198幅COSMO-SkyMed和TerraSAR-X 影像數(shù)據(jù)，對意大利的Pertusillo大壩進行五年的持續(xù)監(jiān)測，并使用地面測量數(shù)據(jù)進行了驗證，結(jié)果表明時序InSAR結(jié)果與實際相符，可以用于獲取大壩的時序變形監(jiān)測?？傮w而言，時序InSAR可以獲取水庫大壩的高精度變形信息，是一種高精度、大范圍大壩安全監(jiān)測的有力技術(shù)手段。

大藤峽水利樞紐作為珠江流域的關(guān)鍵控制性工程，對保障流域防洪和供水安全、促進地區(qū)經(jīng)濟社會可持續(xù)發(fā)展具有重要作用。本文基于Sentinel-1衛(wèi)星SAR數(shù)據(jù)，使用小基線集InSAR（Small Baseline Subset InSAR，SBAS-InSAR）技術(shù)，開展大藤峽水利樞紐工程的時序InSAR變形監(jiān)測，分析了大藤峽左壩主體及正在建設(shè)的右壩主體的時序變形情況，并對大藤峽黔江上游兩岸潛在的滑坡風(fēng)險隱患進行了識別。

一、研究區(qū)域和數(shù)據(jù)

（一）研究區(qū)域

大藤峽水利樞紐工程地處廣西壯族自治區(qū)桂平市，位于珠江流域西江干流黔江河段的大藤峽出口弩灘附近，是紅水河（黔江）梯級規(guī)劃中最末一個梯級，大壩整體呈現(xiàn)Y字形，黔江主壩壩型為混凝土重力壩，最大壩高80.01m，壩長1343.098m，左岸黔江副壩為粘土心墻石渣壩，壩長為1239m。大藤峽水庫所在區(qū)域的平均高程約為23～43m，在正常投入生產(chǎn)使用時，水庫蓄水水位為61.00m，設(shè)計總庫容約為34.3億m3，防洪庫容約15億m3，控制流域面積約19.86萬km2，約占西江流域面積的56.4%。

（二）SAR影像數(shù)據(jù)

雷達數(shù)據(jù)來源于歐空局新一代中高分辨率的 Sentinel-1 A衛(wèi)星。Sentinel-1屬于歐空局哥白尼計劃（GMES）中的地球觀測衛(wèi)星，由A/B兩顆相距180°的衛(wèi)星星座組成（B星于2021年12月發(fā)生故障，已無法獲取數(shù)據(jù)），單星重訪周期為12天，雙星串飛可達到6天的重訪周期。本文以2021年2月17日的SAR數(shù)據(jù)為主影像完成配準(zhǔn)預(yù)處理，選定130米空間基線與62天時間基線標(biāo)準(zhǔn)組成基線集，總共生成222個干涉對，干涉對組成的時空基線分布情況見圖1，其中空間基線最大的長度為129m。

二、研究方法

（一） 數(shù)據(jù)預(yù)處理

Sentinel-1數(shù)據(jù)預(yù)處理流程主要包括以下4個關(guān)鍵步驟：（1）添加精密軌道星歷信息；（2）提取研究區(qū)域的條帶與Burst；（3）影像配準(zhǔn)；（4）研究區(qū)域裁剪與多視處理。

（二）小基線集時序InSAR處理

小基線集時序InSAR數(shù)據(jù)處理流程見圖2。SBASInSAR技術(shù)優(yōu)選較短時空基線的影像集進行干涉，提高了影像干涉相干性，在形變求解過程中，通過DEM誤差改正、三維解纏、殘余大氣分離等流程，利用奇異值分解（Singular Value Decomposition，SVD）得到監(jiān)測點的形變速率。

三、 結(jié)果與分析

（一）研究區(qū)域平均形變速率

按照第2節(jié)處理流程，利用SBAS-InSAR技術(shù)獲取了研究區(qū)域雷達視線向（Line of Sight， LOS）的平均形變速率結(jié)果。結(jié)果顯示，沿江兩岸監(jiān)測點稀少（Sentinel衛(wèi)星波長較短，對植被穿透性較差，在沿岸植被覆蓋山區(qū)的干涉相干性低），而城區(qū)監(jiān)測點密集，監(jiān)測效果好。研究區(qū)域整體較為穩(wěn)定，城市群及線狀高鐵線路的形變量在±10mm之間，無明顯的地面沉降災(zāi)害。

（二）大藤峽主壩時序變形分析

根據(jù)2021年1月、2022年3月及2022年10月大藤峽大壩的SAR影像結(jié)果顯示，其中黑色低值區(qū)域為水體，高亮線狀目標(biāo)為已建成的大壩。結(jié)合2021年9月的光學(xué)影像可知，研究監(jiān)測時段內(nèi)大藤峽水利樞紐工程左壩已經(jīng)建設(shè)完工，而右壩于2022年年底才基本完成修建，地表變化較大。

分析可知，大藤峽左壩上分布有較多的高相干監(jiān)測點，而右壩上高相干監(jiān)測點則較為稀少，造成右壩區(qū)域監(jiān)測點稀少的原因是由于在監(jiān)測時段內(nèi)，右壩圍堰區(qū)建設(shè)施工造成地表快速變化，從而導(dǎo)致該區(qū)域出現(xiàn)干涉失相干，進而造成了高相干監(jiān)測形變點的缺失。但整體而言，利用Sentinel-1影像開展大藤峽大壩時序InSAR變形監(jiān)測是有效可行的。

選取大藤峽大壩上8個重點特征區(qū)域，進一步分析大壩的時序變形狀況，結(jié)果見圖3。p1、p2、p3特征區(qū)域位于左壩主體，3個區(qū)域的累積時序變形見圖3（a）（b）（c），結(jié)果顯示左壩主體上3個區(qū)域的累積變形量較小，左壩主體基本處于穩(wěn)定狀態(tài)。p4、p5區(qū)域位于大藤峽船閘區(qū)域，其時序變形見圖3（d）（e），p4、p5區(qū)域的變形處于-20～20mm之間，變形波動較大，該現(xiàn)象主要是由大量超強散射信號干擾引起。p6、p7區(qū)域位于右壩的圍堰攔水設(shè)施上，其時序變形見圖3（f）（g），結(jié)果顯示P6區(qū)域累積形變量較大，約為65mm，p7區(qū)域變形量基本在0mm附近。p8區(qū)域為右岸的一處斜坡，時序變形見圖3（h），累積變形約為60mm，線性變形趨勢明顯，具有一定危險性。

（三）沿岸斜坡形變分析

大藤峽庫區(qū)黔江上游沿江兩岸的平均形變速率如圖4。該區(qū)域發(fā)現(xiàn)了4處明顯的變形區(qū)：A1、A2、A3、A4區(qū)域，其時序變形圖顯示四處區(qū)域的變形基本呈線性趨勢，疑似為滑坡體活動變形。

四、結(jié)論與展望

本文基于2020年2月至2023年1月年期間大藤峽庫區(qū)的Sentinel-1時序數(shù)據(jù)，利用星載時序SABS-InSAR監(jiān)測方法，開展了大藤峽水利樞紐及其沿江兩岸區(qū)域的變形監(jiān)測研究，研究結(jié)論如下：

（1）基于Sentinel-1影像的SBAS-InSAR技術(shù)可以實現(xiàn)大藤峽區(qū)域的時序變形監(jiān)測，為大藤峽大壩的變形監(jiān)測提供了一種新思路。

（2）大藤峽左壩主體基本處于穩(wěn)定狀態(tài)，大藤峽沿江兩岸發(fā)現(xiàn)4處具有較明顯線性形變趨勢的斜坡。

（3）監(jiān)測結(jié)果對于大藤峽水利樞紐的健康運營具有一定的參考價值，后期聯(lián)合GNSS、精密水準(zhǔn)測量及實地勘察結(jié)果等對大壩安全監(jiān)測分析，將更有助于確保大壩安全。

作者單位：貴州省第一測繪院