黃友菊,覃怡婷,吳 慧,韋 強(qiáng)

(廣西壯族自治區(qū)自然資源遙感院,廣西 南寧 530000)

隨著中國(guó)社會(huì)經(jīng)濟(jì)的飛速發(fā)展,國(guó)家對(duì)水利基礎(chǔ)設(shè)施的建設(shè)和維護(hù)投入力度加大,且建造了各型水庫(kù),以進(jìn)行水土保持、水利安全維護(hù)和水務(wù)管理。形變監(jiān)測(cè)分析和預(yù)測(cè)是建立水庫(kù)安全評(píng)價(jià)機(jī)制的基礎(chǔ)。目前常基于以點(diǎn)概面、以局部情況代替整體形變的原則,對(duì)全站儀水準(zhǔn)測(cè)量、GPS離散點(diǎn)所獲取的形變信息進(jìn)行處理與分析[1]。這種方式無(wú)法獲取實(shí)際的空間連續(xù)形變信息,除了耗費(fèi)大量人力物力資源外,可能還會(huì)遺漏潛在的安全隱患。因此,采用覆蓋范圍廣、可實(shí)現(xiàn)連續(xù)觀測(cè)的遙感技術(shù)監(jiān)測(cè)是發(fā)展的必然趨勢(shì)。

差分雷達(dá)干涉測(cè)量技術(shù)(differential interferometric synthetic aperture radar,D-InSAR)是將同一研究區(qū)不同時(shí)間的兩幅SAR影像通過(guò)差分干涉處理得到地表形變信息。作為近年發(fā)展迅速的對(duì)地觀測(cè)衛(wèi)星遙感技術(shù),D-InSAR具有全天候、大面積同步獲取地形信息及形變信息的特點(diǎn),且具有一定的地表穿透能力[2-4],為水利設(shè)施的建設(shè)和維護(hù)提供了新穎思路和技術(shù)手段。但普通D-InSAR技術(shù)的精度會(huì)受地表散射特性的影響,導(dǎo)致時(shí)間、空間的失相干,同時(shí)大氣延遲也會(huì)引入噪聲[5]。時(shí)序InSAR技術(shù)的提出極大克服了失相干現(xiàn)象,針對(duì)時(shí)間基線較長(zhǎng)的數(shù)據(jù)集和多普勒質(zhì)心頻移變化較大的數(shù)據(jù)具有一定適用性,同時(shí)減緩了大氣相位影像及地形效應(yīng)對(duì)差分的干擾[6]。利用時(shí)序InSAR技術(shù)對(duì)水利設(shè)施監(jiān)測(cè)進(jìn)行研究和分析在國(guó)內(nèi)外得到了廣泛的應(yīng)用。文獻(xiàn)[7]探討了利用InSAR技術(shù)識(shí)別三峽工程蓄水過(guò)程中地表形變的可行性和有效性。文獻(xiàn)[8]對(duì)上海市江口南岸和杭州灣北側(cè)堤壩的沉降速率進(jìn)行監(jiān)測(cè),證明了時(shí)間序列SAR影像在堤壩形變監(jiān)測(cè)上的有效性和應(yīng)用潛力。文獻(xiàn)[9]利用SBAS-InSAR技術(shù)對(duì)挪威斯瓦特湖壩進(jìn)行形變監(jiān)測(cè),探測(cè)到壩體的局部形變情況,精度可達(dá)毫米級(jí)。

現(xiàn)有研究多側(cè)重于對(duì)大型城市群區(qū)域內(nèi)的水利設(shè)施、大型樞紐工程的監(jiān)測(cè),對(duì)位于植被茂盛、地形起伏較大的山區(qū)的大中型水庫(kù)大壩缺乏關(guān)注。運(yùn)用時(shí)序InSAR技術(shù)提取形變時(shí),多數(shù)研究基于先驗(yàn)經(jīng)驗(yàn),只使用單基線集提取形變情況,未考慮不同基線集對(duì)提取形變信息的影響,從而找到最適合研究區(qū)實(shí)際情況的基線條件。本文綜合考慮山區(qū)水庫(kù)大壩特殊的環(huán)境特征和自身特點(diǎn),以貴州省貴陽(yáng)市平寨水庫(kù)及周邊區(qū)域?yàn)檠芯繀^(qū)開(kāi)展試驗(yàn),運(yùn)用時(shí)序InSAR技術(shù),在不同時(shí)間基線集合條件下提取庫(kù)區(qū)形變,從而分析獲取最符合研究區(qū)實(shí)際情況的形變信息。

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

1.1 地理位置

研究區(qū)位于貴州省貴陽(yáng)市南部,范圍約75 km2,平均海拔約1100 m。研究區(qū)地屬亞熱帶季風(fēng)氣候,冬無(wú)嚴(yán)寒,夏無(wú)酷暑,全年域內(nèi)長(zhǎng)存大面積植被。水庫(kù)位于黔南州壩王河上游惠水縣雅水鎮(zhèn)與擺金鎮(zhèn)交接處,壩址在庫(kù)區(qū)北偏東的位置。研究區(qū)域如圖1所示。

圖1 研究區(qū)地理位置

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

選取歐空局長(zhǎng)時(shí)間序列Sentinel-1A升軌雷達(dá)影像數(shù)據(jù)59景,時(shí)間跨度為2019年1月—2020年12月,分辨率為25 m,C波段(波長(zhǎng)5.4 cm),中心入射角為39.17°,數(shù)據(jù)處理采用SBAS-InSAR技術(shù)。成果為時(shí)序InSAR地面沉降結(jié)果,所獲取的地表形變信息包括監(jiān)測(cè)目標(biāo)的地理位置、平均沉降速率、沉降歷史等。

2 基于山區(qū)水庫(kù)研究的時(shí)序InSAR方法

時(shí)序InSAR技術(shù)多應(yīng)用于城市區(qū)域的沉降監(jiān)測(cè)。地處山區(qū)的水庫(kù)因特殊的地域特征給數(shù)據(jù)處理帶來(lái)一定困難,主要影響因素為:①山區(qū)地表多以植被覆蓋,且地形起伏較大,導(dǎo)致影像之間存在一定的時(shí)空失相干,從而使影像干涉效果不佳;②研究區(qū)以大面積水體覆蓋,人工建筑區(qū)域較少且分散,雷達(dá)的側(cè)視成像在含水體量大的范圍內(nèi)后向散射強(qiáng)度低,未合理選取基線條件會(huì)導(dǎo)致試驗(yàn)結(jié)果與研究區(qū)實(shí)際情況不匹配,影響形變結(jié)果提取。

短基線集干涉測(cè)量(small baseline subset interferometric synthetic aperture radar,SBAS-InSAR)技術(shù)[10]作為新型時(shí)序InSAR干涉測(cè)量技術(shù)之一,需要以一定數(shù)量的SAR數(shù)據(jù)作為基礎(chǔ)數(shù)據(jù)源,完成長(zhǎng)時(shí)間內(nèi)對(duì)地持續(xù)、有效的形變監(jiān)測(cè)。其基本原理可表述為:假定在時(shí)間t1,t2,…,tn內(nèi)獲取同一區(qū)域的N幅SAR影像,根據(jù)干涉組合條件﹐形成M幅干涉圖像﹐M與N滿足N/2≤M≤N(N-1)/2。對(duì)于任意干涉圖i,在去除平地及地形相位影響后,第i個(gè)像素的干涉相位可表示為

δφi(x)=φ(tB,x)-φ(tA,x)

(1)

式中,λ為發(fā)射信號(hào)中心波長(zhǎng);tA、tB(tA

Aφ=Δφ

(2)

式中,φ為待求時(shí)段的相位矩陣;Δφ表示M幅差分干涉圖的相位矩陣;A為M×N的系數(shù)矩陣,其中矩陣的行對(duì)應(yīng)一景干涉圖,列對(duì)應(yīng)該時(shí)刻的SAR影像,主影像所在位置為列數(shù)+1,輔影像所在位置為列數(shù)-1,其余為0。A為列滿秩矩陣,且N為A的秩,通過(guò)最小二乘法計(jì)算可得

φ=(ATA)-1ATΔφ

(3)

式(3)為所有SAR影像在一個(gè)集合的理想情況下的最優(yōu)解,而在實(shí)際情況中存在多個(gè)集合,且需要考慮大氣和軌道相位誤差的影響。因此,對(duì)矩陣A進(jìn)行奇異值(singular value decomposition,SVD)分解,解算出該方程式最小范數(shù)意義上的最小二乘解[11];估算地形和大氣延遲相位的影響,將模型相位值與真實(shí)形變相位值作差去除殘差影響;通過(guò)大氣高通、低通濾波減弱大氣相位誤差的影響[12];最后求得目標(biāo)點(diǎn)的形變速率值,反演出地表的形變相位值。由于SBAS-InSAR方法屬于多視成像觀測(cè),降低了形變圖的空間分辨率,因此該方法適用于大范圍面狀或線狀的形變區(qū)域監(jiān)測(cè)[13]。技術(shù)流程如圖2所示。

圖2 SBAS-InSAR技術(shù)流程

在D-InSAR技術(shù)中,時(shí)間基線、空間基線和多普勒質(zhì)心頻移是影響干涉像對(duì)質(zhì)量的重要因素[14]。時(shí)間基線越長(zhǎng),時(shí)間失相干越嚴(yán)重;空間垂直基線越大,空間失相干越嚴(yán)重;多普勒質(zhì)心頻率基線越大,干涉相干性越小。時(shí)間基線為兩幅影像成像的時(shí)間間隔。時(shí)間基線的長(zhǎng)短影響時(shí)間去相干程度,主要是指地表和大氣介電常數(shù)在一段時(shí)間內(nèi)發(fā)生的變化[15]。季節(jié)的交替更迭、氣候的變化等自然與人類(lèi)活動(dòng)都會(huì)影響地表介電常數(shù)[16]。時(shí)間基線對(duì)相干程度的影響遠(yuǎn)遠(yuǎn)復(fù)雜于空間基線對(duì)相干程度的影響。多數(shù)研究表明,時(shí)間基線越長(zhǎng),不可控的隨機(jī)因素造成的影響越大[17]。

3 結(jié)果分析

文獻(xiàn)[18]提出了以最佳時(shí)間基線作為公共主影像選取的依據(jù),即時(shí)間基線和最短[18]。但本文是對(duì)山區(qū)水庫(kù)大壩的形變監(jiān)測(cè)進(jìn)行研究,不能一味選擇最短時(shí)間基線。在氣候濕潤(rùn)或植被茂密的區(qū)域,時(shí)間基線過(guò)短會(huì)導(dǎo)致最后提取形變的區(qū)域過(guò)小,而無(wú)法進(jìn)行結(jié)果研究和分析。本文綜合考慮山區(qū)特殊的環(huán)境特征和自身特點(diǎn),在Sentinel-1A雷達(dá)衛(wèi)星的重訪周期為12 d的基礎(chǔ)上,分別建立時(shí)間基線為36、60、90 d的數(shù)據(jù)集合作對(duì)比,探究時(shí)間基線對(duì)干涉相干性的影響,以及時(shí)間基線對(duì)最后獲取的形變情況的影響。選擇的時(shí)間基線集分別對(duì)應(yīng)重訪周期12的公倍數(shù),同時(shí)考慮了1～3個(gè)月的季節(jié)因素,分析2019—2020年水庫(kù)周?chē)鷧^(qū)域沉降,以及2019、2020年整年水庫(kù)大壩形變情況。

本文針對(duì)Sentinel-1A衛(wèi)星空間基線較短的特點(diǎn),將空間基線閾值統(tǒng)一設(shè)置為臨界基線的45%,使空間基線的長(zhǎng)度在300 m以下,從而保證較好的相干性[19]。本次處理中空間基線閾值為臨界基線的30%。選取2019年時(shí)間基線60 d集合,以獲取時(shí)間為20190317的SAR影像,作為主影像的其中一組像對(duì)進(jìn)行展示。小基線集主從影像對(duì)及參數(shù)見(jiàn)表1。

表1 小基線主從影像對(duì)及時(shí)空基線參數(shù)

如圖3所示,空間基線為±12 d時(shí)相位信息最清晰,其成像季節(jié)氣候大致相同,地表植被反射特性相似,庫(kù)區(qū)水體部分、高反射率的城市建筑區(qū)域位置清晰。-24 d的相位圖像相較-12 d圖像噪聲點(diǎn)增多,除了城市建筑區(qū)域及水體有明顯的信息外,其他區(qū)域像素離散程度增加,多片狀區(qū)域信息殘缺。對(duì)比±12 d相干性(如圖3(d)、(e)所示)及±36 d相干性(如圖3(b)、(f)所示),雖然時(shí)間基線絕對(duì)值相同,但由于季節(jié)變化,植被覆蓋率增加,研究區(qū)整體含水量上升,干涉的相干性受影響更大。圖3(f)中研究區(qū)主要的水體及周邊信息區(qū)域已經(jīng)因?yàn)橥耆喔蔁o(wú)法目視判別,原區(qū)域變?yōu)殡S機(jī)的噪點(diǎn),相干性相較-48 d相干性(如圖3(a)所示)更難識(shí)別。通過(guò)以上分析可以看出,時(shí)間基線過(guò)長(zhǎng),受地表植被覆蓋影響及季節(jié)氣候變化因素造成的時(shí)間失相干越嚴(yán)重,對(duì)干涉的相干性影響越大。

因此,運(yùn)用SBAS-InSAR技術(shù)獲取山區(qū)水庫(kù)研究區(qū)域形變情況,需要對(duì)不同時(shí)間基線數(shù)據(jù)集合進(jìn)行數(shù)據(jù)處理,通過(guò)對(duì)比分析確定獲取形變結(jié)果的最佳時(shí)間基線集。不同時(shí)間基線集合生成的干涉像對(duì)數(shù)量及剔除相干性較差的像對(duì)數(shù)量見(jiàn)表2。

表2 時(shí)間基線集合生成的干涉像對(duì)情況

通過(guò)SBAS-InSAR技術(shù)處理,獲取平寨水庫(kù)研究區(qū)在2019年1月28日—12月30日、2020年1月11日—12月24日兩個(gè)時(shí)間段時(shí)間基線集分別為36、60、90、120 d的年平均形變速率,如圖4—圖5所示。

圖4 2019年平均形變速率

圖5 2020年平均形變速率

分析圖4(a)、圖5(a)可知,時(shí)間基線集為36 d時(shí),SBAS-InSAR技術(shù)只能在圖上提取較少的形變情況,且?guī)靺^(qū)附近的重點(diǎn)研究區(qū)域形變提取量也較少,無(wú)法提取到壩體的具體形變信息。分析圖4(d)、圖5(d)可知,時(shí)間基線集為120 d時(shí),研究區(qū)內(nèi)存在大片形變區(qū)域,但由于時(shí)間基線過(guò)長(zhǎng),干涉流受時(shí)間失相關(guān)的影響,導(dǎo)致存在大量形變信息誤判及部分區(qū)域形變量與現(xiàn)實(shí)情況不符的問(wèn)題,主要體現(xiàn)在研究區(qū)西北山區(qū)及庫(kù)區(qū)南部、東南部山區(qū),監(jiān)測(cè)的形變區(qū)域起伏變化較大且雜亂無(wú)章。采用60 d與90 d的時(shí)間基線對(duì)研究區(qū)進(jìn)行SBAS-InSAR監(jiān)測(cè)提取的形變區(qū)域基本一致,但時(shí)間基線為90 d時(shí),仍存在部分區(qū)域形變突然增大或減少的問(wèn)題。結(jié)合圖3時(shí)序相干性系數(shù)圖分析可知,對(duì)于貴州山區(qū)內(nèi)植被覆蓋率大、含水量高的研究區(qū),時(shí)間基線集選定越大,時(shí)間失相干越嚴(yán)重,對(duì)于形變提取存在的形變量驟降及驟升問(wèn)題影響越大;而60 d時(shí)間基線集提取的形變量更為詳細(xì),因此選定60 d時(shí)間基線集對(duì)平寨水庫(kù)庫(kù)區(qū)形變進(jìn)行SBAS-InSAR監(jiān)測(cè)得到的結(jié)果較為合理。

根據(jù)時(shí)間基線集為60 d的平寨水庫(kù)研究區(qū)2019、2020年形變速率圖所提取的形變信息,選取壩體西側(cè)特征點(diǎn)A與東側(cè)特征點(diǎn)B,如圖6所示。

圖6 壩體形變特征點(diǎn)

將A、B兩點(diǎn)衛(wèi)星視線方向的具體形變量,疊加庫(kù)區(qū)大壩上兩點(diǎn)于2019—2020年內(nèi)所獲取的幾何水準(zhǔn)測(cè)量數(shù)據(jù),并作對(duì)比研究,然后將SBAS-InSAR技術(shù)提取的形變結(jié)果與水準(zhǔn)數(shù)據(jù)作差取絕對(duì)值制圖進(jìn)行分析,結(jié)果如圖7—圖8所示。

圖7 A點(diǎn)水準(zhǔn)數(shù)據(jù)和InSAR數(shù)據(jù)形變差值分析

圖8 B點(diǎn)水準(zhǔn)數(shù)據(jù)和InSAR數(shù)據(jù)形變差值分析

疊加距研究區(qū)最近的貴陽(yáng)站氣象站2019—2020年月平均降雨量數(shù)據(jù)(見(jiàn)表3)進(jìn)行分析[20]。

表3 2019—2020年貴陽(yáng)市月均降水量 mm

由圖7和圖8可知,堤壩上典型點(diǎn)位的最大形變值約為30 mm,形變的發(fā)生主要是受降雨量及水庫(kù)水位變化的外界作用。在2019年4—7月、2020年5—7月、9—10月的3個(gè)時(shí)間段內(nèi),月平均降雨量增大,點(diǎn)位形變速率有明顯抬升,且受降雨導(dǎo)致的庫(kù)區(qū)水位增長(zhǎng)的影響,抬升持續(xù)一段時(shí)間后緩慢回落。60 d和90 d時(shí)間基線集提取的形變信息與水準(zhǔn)測(cè)量數(shù)據(jù)的變化趨勢(shì)具有一致性,但當(dāng)時(shí)間基線集為90 d時(shí),提取的形變量整體偏大,且存在驟升或驟降的問(wèn)題,與水準(zhǔn)測(cè)量數(shù)據(jù)擬合性不佳。對(duì)兩種數(shù)據(jù)的差值進(jìn)行分析,在降雨量較大的月份(5、6、7、9月),地表含水量的增加對(duì)干涉的相干性產(chǎn)生一定的影響,差值也存在一定程度的增大。

4 結(jié) 論

時(shí)序InSAR提供了一種全新的衛(wèi)星遙感監(jiān)測(cè)方法,可以在低成本條件下提取較為精確的形變信息。本文針對(duì)水系豐富、降水充足、山勢(shì)險(xiǎn)峻的研究區(qū)地表形變情況提取困難的問(wèn)題,利用長(zhǎng)時(shí)間序列Sentinel-1A升軌雷達(dá)影像數(shù)據(jù),運(yùn)用時(shí)序InSAR技術(shù),采用SBAS-InSAR方法,提取平寨水庫(kù)周邊形變信息,分析了不同時(shí)間基線集對(duì)庫(kù)區(qū)形變提取的適應(yīng)性,主要結(jié)論如下。

(1)對(duì)于植被覆蓋率及含水量較高的西南山地區(qū)域,時(shí)間基線集為60 d時(shí)提取研究區(qū)形變信息最為準(zhǔn)確,既具有一定形變提取量供后續(xù)分析與研究,同時(shí)也符合實(shí)際的形變情況,存在誤判和缺漏較其他時(shí)間基線集少。

(2)平寨水庫(kù)研究區(qū)內(nèi)水工建筑物及周邊城鎮(zhèn)等人工建筑的形變受降雨與庫(kù)區(qū)蓄水量的綜合影響,存在±30 mm的形變情況。根據(jù)壩體上提取的特征點(diǎn)形變情況,壩體形變受降雨量影響,降雨量大的月份存在明顯的形變抬升。

(3)與幾何水準(zhǔn)測(cè)量數(shù)據(jù)進(jìn)行對(duì)比分析發(fā)現(xiàn),兩種數(shù)據(jù)監(jiān)測(cè)結(jié)果在趨勢(shì)上具有一致性,60 d時(shí)間基線集SBAS-InSAR提取的形變信息與水準(zhǔn)監(jiān)測(cè)數(shù)據(jù)差值的絕對(duì)值,除了別時(shí)間點(diǎn)偏差較大外,其他時(shí)段均保持在5 mm以內(nèi),偏差較小。本文試驗(yàn)結(jié)果驗(yàn)證了SBAS-InSAR技術(shù)監(jiān)測(cè)結(jié)果的準(zhǔn)確性。

由于缺少精細(xì)的大壩結(jié)構(gòu)數(shù)據(jù),本文數(shù)據(jù)分析僅基于幾何水準(zhǔn)監(jiān)測(cè)及SBAS-InSAR技術(shù)成果,同時(shí)考慮了降雨量、氣候等自然環(huán)境因素對(duì)地表形變提取的影響,未加入人為施工開(kāi)發(fā)導(dǎo)致的水土流失及自然災(zāi)害等因素引起的形變,后續(xù)研究將加入人文環(huán)境因素進(jìn)行綜合分析,實(shí)現(xiàn)對(duì)山區(qū)水庫(kù)的維護(hù)監(jiān)測(cè)及災(zāi)害預(yù)警。