張永滎，張生鵬,2，貢寶扎西，李水平，李亞剛，省天琛

（ 1. 青海省基礎(chǔ)測繪院, 西寧 810001；2. 青海師范大學(xué) 地理科學(xué)學(xué)院, 西寧 810008；3. 合肥工業(yè)大學(xué) 土木與水利工程學(xué)院, 合肥 230009 ）

0 引 言

水電站是通過調(diào)控水資源時空分布成為蓄水防旱的重要樞紐工程，其運(yùn)行直接關(guān)系到庫區(qū)及下游居民的重大生命財(cái)產(chǎn)安全[1]. 中國已建成各類水庫近10萬座，但是多數(shù)水庫由于建造年代久遠(yuǎn)，受限于當(dāng)時的條件，病險(xiǎn)情況普遍存在，最嚴(yán)重時其數(shù)量占比達(dá)到30%. 形變監(jiān)測是水庫安全運(yùn)營評估的基礎(chǔ)，目前水庫形變監(jiān)測主要采用傳統(tǒng)全站儀、水準(zhǔn)測量等手段在重點(diǎn)部位進(jìn)行監(jiān)測，傳統(tǒng)大地測量的工作特點(diǎn)決定其監(jiān)測結(jié)果空間采樣率不足及觀測成果具有周期性，導(dǎo)致在大范圍監(jiān)測結(jié)果中無法及時有效地發(fā)現(xiàn)新的形變區(qū)域，造成重大安全隱患的遺漏[2-3].

全球衛(wèi)星導(dǎo)航系統(tǒng)(GNSS)和合成孔徑雷達(dá)干涉測量(InSAR)技術(shù)作為大地測量新技術(shù)，為水庫形變監(jiān)測帶來了新的方案. 王騰等[4]采用InSAR永久散射體(PS)技術(shù)主要對三峽水庫大壩進(jìn)行了穩(wěn)定性分析，得出利用此種技術(shù)可以獲得遠(yuǎn)高于傳統(tǒng)手段監(jiān)測點(diǎn)的結(jié)論. 裴媛媛等[5]采用PS技術(shù)對上海部分堤壩進(jìn)行了穩(wěn)定性分析，得出此方法可以得到毫米級精度，更好地用在人工線狀地物的結(jié)論. 張明晶[6]對于橫泉水庫GNSS自動化監(jiān)測進(jìn)行研究，得出此類系統(tǒng)可以提高監(jiān)測的實(shí)時性和可靠性的結(jié)論. MILILLO等[7]利用InSAR時序方法研究了伊拉克摩蘇爾大壩2004—2016失穩(wěn)過程，找出了各個時段詳細(xì)的形變趨勢. 上述研究都是單獨(dú)利用GNSS或InSAR進(jìn)行監(jiān)測，其中InSAR基本采用PS技術(shù)，聚焦點(diǎn)在水庫人工建筑，GNSS監(jiān)測未單獨(dú)驗(yàn)證北斗衛(wèi)星導(dǎo)航系統(tǒng)(BDS)的可用性，也未融合現(xiàn)有連續(xù)運(yùn)行參考站系統(tǒng)(CORS)技術(shù)的優(yōu)勢，其實(shí)時定位精度停留在厘米級. 龍羊峽水電站作為黃河上游第一座大型梯級電站，對下游水電站的水量調(diào)節(jié)起主導(dǎo)作用，由于其特殊的地形、地貌、巖土體類型及氣象因素的誘發(fā)，成為大型滑坡的高發(fā)育區(qū)，并且迄今為止還完全采用傳統(tǒng)測量方式進(jìn)行形變監(jiān)測[8]. 本文綜合考慮GNSS的高時間分辨率和InSAR的高空間分辨率，結(jié)合青海省連續(xù)運(yùn)行參考站系統(tǒng)(QH-CORS)數(shù)據(jù)成果采用InSAR小基線集技術(shù)對龍羊峽庫區(qū)進(jìn)行了地表形變監(jiān)測分析，驗(yàn)證并采用InSAR技術(shù)進(jìn)行大面積形變區(qū)排查，基于CORS成果提升BDS實(shí)時動態(tài)(RTK)差分精度，進(jìn)行重點(diǎn)區(qū)域形變RTK差分監(jiān)測方法的可行性及精度.

1 研究方法

1.1 BDS RTK差分監(jiān)測

基于CORS成果進(jìn)行BDS實(shí)時動態(tài)差分計(jì)算的步驟如下：

步驟一：將基準(zhǔn)站的觀測值和實(shí)測氣象元素按特定頻率匯集到數(shù)據(jù)中心，建立新的對流層延遲改正模型進(jìn)行RTK定位；

步驟二：選擇經(jīng)驗(yàn)?zāi)Ｐ?，?jì)算所有站前一段時間每天的全球經(jīng)驗(yàn)?zāi)Ｐ秃蛯?shí)測氣象元素的天頂對流層延遲(ZTD)；

步驟三：建立高程、基于全球經(jīng)驗(yàn)?zāi)Ｐ秃突趯?shí)測氣象元素對流層延遲的模型之間的關(guān)系，進(jìn)行新的修正模型；

步驟四：利用改正后的基準(zhǔn)站和監(jiān)測站ZTD模型得到基準(zhǔn)站觀測值的改正量；

步驟五：處理后的基準(zhǔn)站觀測值播發(fā)給監(jiān)測站，監(jiān)測站進(jìn)行實(shí)時差分解算出三維(3D)坐標(biāo).

1.2 短基線子集技術(shù)

短基線子集(SBAS)技術(shù)于2002被提出，是當(dāng)前較主流的InSAR時間序列分析方法[9]. 該方法依據(jù)時空基線的條件，自由組合基線較短的影像對，組合成若干個短基線集合，再利用奇異值分解方法，將組合成的所有短基線集聯(lián)合起來求解，從而得到研究區(qū)域內(nèi)整個觀測時間的形變時間序列和平均沉降速率.SBAS的基本步驟如下.

采集同一區(qū)域按照時間順序t0， ···，tN排列的N+1幅雷達(dá)影像，選取其中一幅影像作為主影像進(jìn)行配準(zhǔn)，則形成M幅差分干涉圖.

對于從影像tA和主影像tB(tB＞tA) 時刻獲取的合成孔徑雷達(dá)(SAR)影像生成的第j幅差分干涉圖，在不考慮大氣延遲相位、殘余地形相位和噪聲位的情況下，其中任意像元 (x,y) 的干涉相位為

將式(1)中相位采用兩個獲取時間之間平均相位速度和時間的乘積的方式表示，則有

式中： Φ 為殘余地形相位，得第j幅干涉圖的相位值為

利用矩陣形式表示為

B是M×N矩陣，對B進(jìn)行奇異值求解，獲取最小范圍下的視線向的相位平均速率，對各時間段的沉降速率進(jìn)行時間域上的積分，得到視線向上時間序列的形變量，通過三角函數(shù)，將它們轉(zhuǎn)為垂直向上監(jiān)測時間段內(nèi)的累計(jì)沉降量和形變速率[10-12].

2 研究區(qū)域概況

龍羊峽水電站位于青海省共和縣與貴德縣之間的黃河干流上，黃河穿越其中，右岸是險(xiǎn)峻的茶納山，左岸是莽原. 屬高原干旱-半干旱大陸性氣候，庫區(qū)植被稀少，水位基本每年7月開始上升，11月份達(dá)到最高值，然后開始逐漸回落. 水庫右岸發(fā)育北東東向活動斷裂和北西向斷裂，巖土體類型為雙層結(jié)構(gòu)巖土體，下部軟弱的層狀半固結(jié)巖土體，上部為松散沖洪積一般粘性土[8，13].

研究區(qū)及GNSS監(jiān)測點(diǎn)布設(shè)位置如圖1所示，其中紅色框內(nèi)為龍羊峽水庫庫區(qū)，紅色三角點(diǎn)為QHCORS龍羊峽站及GNSS監(jiān)測點(diǎn).

圖1 研究區(qū)及BDS站點(diǎn)布設(shè)圖

3 成果獲取及分析

3.1 GNSS數(shù)據(jù)處理

在龍羊峽大壩下方布設(shè)了5個GNSS監(jiān)測點(diǎn)，將數(shù)據(jù)通過無線網(wǎng)絡(luò)傳實(shí)時輸?shù)絈H-CORS核心服務(wù)器上. 選取QH-CORS龍羊峽站為基準(zhǔn)站，按照1.1節(jié)中的研究方法對采集到的BDS數(shù)據(jù)進(jìn)行RTK差分定位解算. 其步驟一中特定頻率為1 Hz，步驟二中全球?qū)α鲗幽Ｐ蜑镚PT2w，利用薩斯塔莫寧(Saastamoinen)模型計(jì)算龍羊峽站每天各時刻的ZTD. 圖2為流程及具體方法解算出2019年BDS數(shù)據(jù)RTK差分定位形變結(jié)果.

圖2 BDS監(jiān)測點(diǎn)實(shí)時動態(tài)解算流程圖

利用GNSS定位GAMIT軟件解算收集的與SAR和BDS數(shù)據(jù)時間相對應(yīng)的13天的靜態(tài)數(shù)據(jù)[14]，為了最大程度的削弱水域?qū)馑愠晒挠绊?，利用商業(yè)軟件預(yù)算數(shù)據(jù)得到高精度的先驗(yàn)坐標(biāo)，進(jìn)行潮汐改正.采用目前模擬對流層狀況精度最高的模型進(jìn)行大氣改正，其主要設(shè)置為：

衛(wèi)星軌道：國際GNSS服務(wù)(IGS)快速精密星歷；

衛(wèi)星截止高度角：15°；

數(shù)據(jù)采樣間隔：10 s；

坐標(biāo)約束：CORS水平方向給予5 cm，垂直方向給予10 cm的約束，監(jiān)測點(diǎn)給予10 m的約束；

潮汐改正：otl_FES2004.grid；

氣象模型：VMF1[15].

基線重復(fù)性可以反映出基線解內(nèi)部精度，是評定數(shù)據(jù)處理質(zhì)量的重要指標(biāo). 表1顯示基線重復(fù)性在各方向的固定誤差均小于2 mm，相對精度都達(dá)到10-7水平，滿足形變監(jiān)測毫米級精度的要求[14].

表1 基線向量重復(fù)性統(tǒng)計(jì)表

3.2 SAR數(shù)據(jù)處理

選擇覆蓋研究區(qū)域的2017-10-15—2019-10-17獲取的25景Sentinel-1SAR影像，C波段，分辨率5 m×20 m，精密軌道文件，30 m分辨率的數(shù)字高程模型數(shù)據(jù).

采用SBAS技術(shù)解譯獲取的全部影像，構(gòu)建了一個所有時相配對、連接均勻的SBAS數(shù)據(jù)集合，如圖3所示，主要參數(shù)：

圖3 解纏像對連接圖

極化方式：VV；

臨界基線最大百分比：3%；

最大時間距離：60；

多視視數(shù)：1∶5；

解纏方式：3D，Delaunay MCF；

相關(guān)系數(shù)閾值：0.35；

大氣低通：1 200；

大氣高通：365；

形變投影：垂直.

采用的解纏方法可以很好地處理兩個相干性高的孤立區(qū)域，更好地抑制水汽對解譯的影響.

圖4通過解譯獲得了龍羊峽庫區(qū)2017—2019年在垂直方向上的形變速率，區(qū)域內(nèi)大面積沉降區(qū)主要在入庫區(qū)、查西、旅游公司區(qū)域. 最嚴(yán)重沉降區(qū)域位于入庫區(qū)，最大沉降速率52.48 mm/a；最大抬升點(diǎn)在查西，最大抬升速率43.60 mm/a. 其全部分布于庫區(qū)的南岸，與其他學(xué)者研究龍羊峽滑坡集中發(fā)育在南岸的論述相符合[16-17].

3.3 成果分析

3.3.1 地表形變原因推斷

入庫區(qū)和旅游公司所在區(qū)域巖土主要以粘土和湖積粉砂為主. 黃河在龍羊峽內(nèi)直流而下進(jìn)入入庫區(qū)，在其沖擊作用下將部分泥沙沖入庫區(qū)，并在黃河下切影響下，沖刷南岸造成沉降現(xiàn)象. 旅游公司地勢平緩，在庫區(qū)水的浸潤、浪蝕作用下形成一個沉降量較小的沉降形變. 查西在黃河下切的影響下形成臨空面，馬蹄形地勢造成沉降和抬升的同時存在水位下降，山體內(nèi)地下水位下降幅度相對滯后，地下水向庫區(qū)滲流，產(chǎn)生動水壓力，增大山體的下滑力[8，13].

3.3.2 成果精度驗(yàn)證

所得成果的精度采用外符合精度評定. 外符合精度是以外部提供的參考值為比對基準(zhǔn)，主要反映觀測值與參考值之間的偏差程度. 外符合精度計(jì)算公式如下：

式中：m外是 外符合精度； Δ 是參考值減去觀測值后的殘差；n為監(jiān)測點(diǎn)參與計(jì)算數(shù)據(jù)期數(shù)[15].

將靜態(tài)解算成果與BDS RTK監(jiān)測成果相減得到的殘差值代入式(5)，得到各個點(diǎn)的外符合精度值.然后對各點(diǎn)外符合精度值求平均值得到監(jiān)測精度. 同理，得到InSAR監(jiān)測成果的精度.

由表2～3可知，BDS RTK差分監(jiān)測成果的精度為7.1 mm，InSAR監(jiān)測成果的精度為4.4 mm. 以GNSS靜態(tài)解算成果為參考值，兩種技術(shù)的監(jiān)測精度都達(dá)到毫米級.

表3 InSAR監(jiān)測成果精度計(jì)算表 mm

4 結(jié)論和討論

通過InSAR提取高空間分辨率的地表形變信息，GNSS獲取高時間分辨率的空間3D信息. 近年來國家對BDS系統(tǒng)的建設(shè)投入越來越大，而CORS的迅猛發(fā)展推動了其廣闊的應(yīng)用前景. 本文以龍羊峽水庫為研究區(qū)域，通過兩種技術(shù)獲取了整個庫區(qū)垂直方向的地表形變信息，分析結(jié)果表明：在黃河影響下，2017—2019年龍羊峽庫區(qū)存在大面積的地表形變區(qū)，最嚴(yán)重沉降區(qū)域位于入庫區(qū)，最大沉降速率為52.48 mm/a；最大抬升點(diǎn)在查西，最大抬升速率為43.60 mm/a. 以GNSS靜態(tài)解算成果為參考值進(jìn)行驗(yàn)證，BDS RTK差分監(jiān)測成果的精度為7.1 mm，In-SAR監(jiān)測成果的精度為4.4 mm，都滿足形變監(jiān)測的精度要求. 采用InSAR進(jìn)行大面積形變區(qū)篩選，再結(jié)合CORS利用BDS進(jìn)行重點(diǎn)區(qū)域?qū)崟r監(jiān)測的模式對于水庫全方位形變監(jiān)測是可行的. 由于InSAR對平面信息不敏感，本文只對垂直方向進(jìn)行了驗(yàn)證分析，但是庫區(qū)全方位形變監(jiān)測在進(jìn)行垂直形變監(jiān)測的同時進(jìn)行平面形變監(jiān)測也十分重要[18]，今后結(jié)合GNSS空間3D信息及水汽反演成果，對深入研究InSAR水汽干擾消除及3D形變場的建立很有意義.