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

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

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

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

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

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

（一）研究區(qū)域

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

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

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

二、研究方法

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

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

（二）小基線集時(shí)序InSAR處理

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

三、 結(jié)果與分析

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

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

（二）大藤峽主壩時(shí)序變形分析

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

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

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

（三）沿岸斜坡形變分析

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

四、結(jié)論與展望

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

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

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

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

作者單位：貴州省第一測(cè)繪院