利用衛星測高監測高郵湖水位變化

梁子亮，岳建平

（河海大學地球科學與工程學院，江蘇南京 210098）

一、引 言

衛星測高以衛星為載體，以海面作為遙測靶，由衛星上裝載的雷達測高儀向海面發射微波信號，該雷達脈沖信號傳播到達海面后，經過海面反射，再返回到雷達測高儀［1］。通過計算信號往返的時間可以得到衛星高度的測量值［2］。

相較于常規方法，衛星測高具有實時、全天候的特點，尤其在偏遠的內陸湖地區，利用衛星測高可以近實時監測水域的水位變化，建立水域水位的多年時間序列［3］。近年來利用衛星測高技術監測內陸湖泊水域變化，國內已經做了很多研究。褚永海等利用Jason-1數據分析了呼倫湖水位變化［3］；李建成等利用ENVISAT的數據分析了長江中下游湖泊水位的變化［4］；高永剛等利用波形重定后的T/P和Jason-1衛星測高數據對呼倫湖水位再次進行了分析［5］；姜衛平等利用ENVISAT數據分析了青海湖水位變化［6］；李景剛等利用Jason-2衛星測高數據分析了南洞庭湖的水位變化［7］。本文利用 Jason-1、Jason-2衛星測高數據分析了高郵湖2003—2012年的水位變化情況。

二、研究區與數據

1.高郵湖概況

高郵湖是江蘇省第三大淡水湖、中國第六大淡水湖，總面積約780 km2。它的湖盆高出東部里下河平原1.0～2.5 m。湖水較淺，最大水深2.40 m，平均水深僅1.44 m［8］。高郵湖位于淮河下游地區，為淺水湖泊，處于淮河入江水道的中段，是淮水入江的主要通道。高郵湖主要承泄淮河70%以上的來水，是典型的過水型湖泊，行洪期間河湖連為一片［9］。圖1為 Google Earth環境下顯示的Jason-1、Jason-2衛星通過洞庭湖的地面軌跡。在每個Cycle內，Jason-1、Jason-2總共有254個 Pass文件，只有Pass153弧段經過高郵湖。

圖1 Jason-1、Jason-2衛星通過高郵湖的軌跡

2.測高數據

本文采用 OSTST發布的 Jason-1、Jason-2的GDR數據，兩類數據均從AVISO網站下載。GDR數據中，每個 Cycle包含154個 Pass文件。每個Pass文件中所包含的數據記錄個數與衛星通過的地面狀況有關。每個記錄中包含了時間、位置、高度、姿態、測高、環境改正、有效波高、海況改正、后向散射系數等數據項，每個數據項的定義、單位、缺省值和取值范圍可以參考衛星測高產品手冊。Jason-1、Jason-2衛星只有Pass153弧段經過高郵湖。Jason-1測高數據為從2003年1月3日至2009年1月22日（Cycle036～Cycle259）的GDR-C數據。Jason-2測高數據為從2009年2月1日至2012年12月29日（Cycle000～Cycle128）的 GDR-T數據。Jason-2為Jason-1的后續衛星，基本特征與Jason-1相同，一個軌道周期為10 d，地面軌跡在赤道上間距約為315 km。

三、數據分析

1.湖水位計算

精確測定高度及發射的脈沖往返于衛星與湖面的雙程傳播時間，便可以確定衛星到水面的距離。該觀測值與衛星相對于參考橢球面的高度均可以從GDR數據文件中獲取。將這兩個觀測值相減，可以得到湖水面相對于參考橢球的高度，即

式中，H表示湖面橢球高；Halt表示衛星相對于參考橢球的距離；Hran表示衛星到湖面的距離；δcor表示誤差改正。

在開闊海域，GDR數據中提供了包括固體潮、電離層、對流層（干、濕）、海潮、極潮、逆氣壓和海況偏差等在內的多項改正。對于內陸湖泊，測高數據內的最優改正模型并不完全有效和適用，濕對流層和電離層改正經常缺損，這是將測高數據用于監測湖泊精度受到限制的主要因素［3］。本文改正包括對流層（干、濕）、電離層、海況、固體潮、極潮5項改正。式（1）中誤差改正為

式中，δwet為濕對流層改正；δdry為干對流層改正；δlono為電離層改正；δset為固體潮改正；δpol為極潮改正；δssb為海況改正。

2.數據編輯

在內陸水域，由于其面積較小，并且相對于開闊海面，衛星測高的數據量相對較少。介于此，海洋上的數據編輯標準并不能完全適用。在本文研究中，采用BRAT軟件對讀取的GDR數據進行編輯，標準為:① 數據位于湖泊水面范圍內，根據surface_type==1進行判斷；② 測高儀所觀測的高度有效，根據qual_alt_1hz_range_ku==0進行判斷；③各項改正的數值在有效的范圍內；④同一Cycle中的GDR數據與同一Cycle平均值的差值在3倍中誤差的限差范圍內；⑤剔除經上述判斷后所得數據中的粗差［5］。

由于Jason-2衛星發射升空后，Jason-1衛星仍在軌。從2008年7月8日至2009年1月22日兩衛星存在重疊數據，為了選擇質量較好的數據，本文對兩衛星在此期間的數據進行了分析，結果見表1。分析得出Jason-1的數據質量優于Jason-2數據，其原因可能為Jason-2衛星此時仍處于驗證階段，數據質量相對較低。選擇Jason-1數據作為2008年7月8日至2009年1月22日的湖面高計算數據。

表1 2008年7月8日至2009年1月22日間Jason-1、Jason-2湖面高數據統計 m

根據前文數據編輯標準，對測高數據進行編輯，得到時間跨度為2003年3月1日至2012年12月29日共310個湖面高數據。分析發現，每年的1月、2月及3月存在部分數據為負數，與實際不符，原因可能為冬季湖面結冰，衛星高度計跟蹤模式不同所造成的，應 剔 除 該 異 常 數 據。在 （32°58＇42″N，119°23＇55″E）附近的數據同樣存在部分不符值，原因可能為該地區有部分陸地延伸進湖泊內，衛星高度計在測量的過程中產生了誤差，以也應將其剔除。表2是將所有數據進行統計的結果，可以看出，如果剔除粗差，計算湖面平均橢球高度的標準差約0.799 4 m，精度明顯優于剔除前。

表2 高郵湖2003—2012年湖面高數據統計表 m

四、計算結果

瞬時湖面高由兩個分量之和構成:①相對于參考橢球面的大地水準面高；②相對于大地水準面的湖面高。由于內陸湖泊地區的大地水準面通常很難確定，本文將每一個瞬時湖面高減去平均湖面高，從而得到時間序列，并以此來分析水位變化，該過程同樣可以消除部分誤差的影響［3］。本文使用的平均湖面高是將9年（2003—2012年）的瞬時湖面高求取平均值而獲得的。

1.年水位異常

高郵湖位于淮河下游地區，其水位受淮河影響嚴重。將每一年內所有瞬時湖面高取平均值，減去10年平均湖面高，得到年水位異常（如圖2所示）。年水位較高處出現在2003年、2005年和2007年，這與2003年、2005年和2007年淮河發生的洪澇災害相符［10］。其中2003年水位明顯高于其余各年，這也與當年在高郵湖發生的特大洪澇災害相符，該年淮河水位達到了歷史第二高［11］。從圖2中可以看出，10年內高郵湖水位呈下降趨勢，自2009年之后下降趨勢尤為明顯，這也與近年來降水量呈減少趨勢相符［12］。

圖2 高郵湖2003—2012年年水位時間序列

2.月水位異常

高郵湖位于淮河下游地區，處于淮河入江水道的中段，其水位受淮河影響嚴重，每當洪水季節，水位升高；枯水季節，水位下降。從圖3可以看出，高郵湖水位具有明顯的季節性變化特征，每年的最高水位出現在7月。從5月開始，由于梅雨季節來臨，水位逐漸上升；7、8月達到最高峰，這個時期也是淮河流域洪峰時期。9月以后水位逐漸降低，進入平水期，11月至翌年的4月水位為年內最低值。每年的冬季和春季水位較低且穩定，變幅不大，這與淮河流域水位變化有明顯的相關性。

2003年6 月29 日至2003年7月4日和2005年7月4日至2005年7月11日淮河流域（31°～36°N、112°～36°E）出現了連續性的強降水［13］。2007年江淮流域從6月19日開始入梅，降雨增加，特別是2007年6月29日至2007年7月10日，淮河流域出現了集中強降水［11］，受淮河及降水影響，高郵湖水位連續超警戒水位。上述實際情況與圖中反映出的水位起伏相符。

圖3 高郵湖2003—2012年月水位時間序列

3.水位異常序列的Fourier分析

首先取各個半年的水位異常數據作為輸入數據，利用 Fourier變換對高郵湖湖泊水位（2003—2012年）的時間序列進行分析。在分析中由頻率和幅度構成頻率-幅度圖，用于確定水位變化，從中發現湖泊水位變化周期項。圖4給出了高郵湖水位變化序列的頻譜分析，振幅圖中有一個較大的峰值，由此推斷高郵湖水位具有較強的周期性，且周期為1.5年。

圖4 高郵湖水位變化的頻譜分析圖

五、結束語

本文利用Jason-1測高衛星6年（2003—2009年）和 Jason-2測高衛星3年（2009—2012年）的GDR數據，經數據編輯，得到了高郵湖年水位異常和月水位異常的時間序列，并利用Fourier變換分析了水位變化的周期。從兩個時間序列來看，高郵湖水位呈現下降趨勢，且自2009年之后尤為嚴重。每年的7、8月高郵湖水位達到峰值，且在9月以后水位逐漸降低。分析發現:高郵湖水位呈現出顯著的周期變化，周期為1.5年。高郵湖位于淮河下游地區，其水位受淮河影響嚴重。本文利用衛星測高資料所得出的高郵湖水位變化序列，與近幾年高郵湖及淮河水位實際情況相符。

［1］ 高永剛，岳建平，郭金運，等.衛星測高波形重定方法的改進研究［J］.測繪通報，2008（11）:7-9.

［2］ 翟國君，黃謨濤，謝錫君，等.衛星測高數據處理的理論與方法［M］.北京:測繪出版社，2000.

［3］ 褚永海，李建成，姜衛平，等.利用Jason-1數據監測呼倫湖水位變化［J］.大地測量與地球動力學，2005，25（4）:11-16.

［4］ 李建成，褚永海，姜衛平，等.利用衛星測高資料監測長江中下游湖泊水位變化［J］.武漢大學學報:信息科學版，2007，32（2）:144-147.

［5］ 高永剛，郭金運，岳建平.衛星測高在陸地湖泊水位變化監測中的應用［J］.測繪科學，2008，33（6）:73-75.

［6］ 姜衛平，褚永海，李建成，等.利用ENVISAT測高數據監測青海湖水位變化［J］.武漢大學學報:信息科學版，2008，33（1）:64-67.

［7］ 李景剛，李紀人，阮紅勛，等.Jason-2衛星測高數據在陸地水域水位變化監測中的應用——以南洞庭湖為例［J］.自然資源學報，2010，25（3）:502-510.

［8］ 孫羊林，姚志剛.江蘇省十大湖泊［J］.森林與人類，2008（9）:54-57.

［9］ 張煒，周雪晴.湖泊保護規劃中行洪通道的研究——利用水力計算確定高郵湖行洪通道［J］.水利規劃與設計，2009（5）:16-17.

［10］ 王東勇，張嬌，朱紅芳，等.2007年等淮河洪澇年7月上旬環流特征與多年平均環流特征差異分析［C］∥第四屆全國災害性天氣預報技術研討會論文集.北京:中國氣象學會，2007.

［11］ 矯梅燕，金榮花，齊丹.2007年淮河暴雨洪澇的氣象水文特征［J］.應用氣象學報，2008，19（3）:257-264.

［12］ 袁喆，楊志勇，鄭曉東，等.近50年來淮河流域降水時空變化特征分析［J］.南水北調與水利科技，2012，10（2）:24-29.

［13］ 張嬌，王東勇，田紅，等.2003年與2005年淮河流域強降水過程環流特征的對比分析［J］.氣象科學，2008，28（4）:402-408.