基于GF-3衛星SAR數據的鄭州及豫北地區“21·7”暴雨洪澇遙感監測

郭 鵬， 耿維成， 張翠萍

(1.中國氣象局·河南省農業氣象保障與應用技術重點實驗室，鄭州 450003； 2.河南省氣象科學研究所，鄭州 450003；3.國家衛星氣象中心，北京 100081； 4.鄭州市氣象局，鄭州450003)

引 言

隨著全球氣候變化和城市化進程的不斷發展，近年來大城市及其周邊區域暴雨洪澇災害頻繁發生[1-5]，嚴重危及人民群眾生命財產安全和城市的正常運行。衛星遙感技術具有覆蓋面廣、時效性強等特點，可對洪澇發生區域及時進行大范圍掃描，獲取洪澇水體范圍、淹沒范圍、水量水深、受災影響等多種要素[6-7]，實現對洪澇的快速響應、持續觀測和影響分析，提升防汛工作的科學決策能力。

衛星遙感監測水體主要原理是根據衛星傳感器接收的從地表水體反射的電磁波的波段反射特征，提取識別水體信息[8]。在各種遙感水體監測方法中，應用較多的是歸一化水體指數(NDWI)方法。Mcfeeters[9]在歸一化植被指數(NDVI)基礎上，構建了歸一化差異水體指數(NDWI)。該指數采用綠光波段與近紅外波段的比值來最大程度地抑制植被信息并突出水體信息。徐涵秋[10]在NDWI基礎上，提出了改進的歸一化差異水體指數(MNDWI)，并用于不同類型的遙感影像實驗。實驗結果表明，提取的大部分水體信息特別是城鎮范圍內的水體信息比NDWI的效果好。莫偉華等[11]根據MODIS傳感器波段特點，構建了混合水體指數(CIWI)，實現水體信息快速有效提取，效果和精度均優于NDWI和MNDWI的效果和精度。李博等[12]對多日MODIS數據中的典型地物進行采樣分析，在歸一化植被指數法(NDVI)基礎上，將其與MODIS第1波段數據組合得到山區水體增強指數(MEWI)，實現了圖像中山區溝壑的模糊處理，突出了水體信息。以上研究主要利用可見光-近紅外波段光譜數據進行計算得到NDWI或其他水體指數，確定水體判識閾值，確定洪澇淹沒面積，通過對比災前災后水體范圍變化，對災情特征進行分析。

不同類型的遙感數據源在洪澇災情監測中具有各自的優勢，目前應用較多的主要有風云系列氣象衛星、高分辨率系列衛星、EOS/MODIS系列衛星等[13-20]。由于降雨過后洪澇淹沒地區通常被云層覆蓋，可見光-近紅外波段的衛星傳感器均無法穿透云層，難以及時獲得洪澇信息。微波傳感器不受天氣制約，具有全天候、全天時的數據獲取能力[21]，能夠穿云透霧，迅速為災區提供第一手信息。目前利用微波雷達遙感數據開展洪澇監測較少，SSM/I等被動微波遙感數據對地表水體敏感[22-23]，但空間分辨率較低，使其應用受到一定限制。國產高分三號衛星于2016年8月發射升空，是我國首顆空間分辨率達到1 m的多極化合成孔徑雷達(SAR)成像衛星，具有12種成像模式[24-25]，可提供1 m至500 m分辨率、10 km至650 km幅寬的微波遙感數據，實現全球海洋和陸地資源的全天候、全天時監測，為洪澇監測提供了新的技術途徑。本研究擬利用高分三號主動雷達衛星數據，對河南鄭州及其周邊地區在2021年7月20日和7月21-22日豫北地區特大暴雨形成的城市內澇和河流洪澇淹沒地區進行遙感監測，計算城區內澇淹沒面積，以及賈魯河下游河道決堤淹沒面積，為決策、防汛救災提供依據。

1 研究區域與數據

鄭州市是河南省省會，位于黃淮海平原中心地區，常住人口超過1000萬，是全國重要的鐵路、航空、電力樞紐城市。2021年7月17日至22日鄭州市遭遇持續強降雨，全市184個氣象站降水量均超過200 mm，其中180個站的降水量超過250 mm。鄭州國家氣象站最大小時降水量達201.9 mm(20日16-17時)，突破中國大陸2418個國家級氣象站小時雨量歷史極值。另外，7月21日至22日，受切變線影響，河南中北部出現連續強降水天氣，新鄉市等地出現特大暴雨。受此次暴雨影響，鄭州市及豫北多條河流超過防汛保證水位，并淹沒部分城區及低洼地區，防汛形勢嚴峻。

此次“21·7”特大暴雨降水強度大、影響范圍廣，利用衛星遙感技術可第一時間對暴雨產生的洪澇災害作出響應。由于天氣原因，研究區在暴雨結束后長時間被云層覆蓋，極大限制了可見光—近紅外波段遙感傳感器的應用。為及時獲取地表洪澇信息，并為防汛救災提供決策參考，河南省生態氣象和衛星遙感中心聯合國家衛星氣象中心，分別選取暴雨前后的高分三號合成孔徑雷達(SAR)數據，對河南鄭州及豫北地區洪澇淹沒情況進行監測。

2 洪澇水體遙感識別方法和結果

由于不同地物具有不同的電磁散射特性，合成孔徑雷達(SAR)后向散射強度與地物表面粗糙度密切相關。在SAR圖像中，由于水體表面光滑，對雷達信號后向散射能力較弱，因此呈現暗色調，與其他地物類型有較大差別。水體閾值的判定主要有經驗法、雙峰法等[26]。此外，由于地物陰影的反射系數也較低，一般情況下在SAR圖像上水體和陰影難以有效區分，但本文研究區都位于平原及城市周邊地區，因此較少考慮地形起伏及陰影影響。

利用中國氣象局國家衛星氣象中心開發的高分辨率衛星資料氣象應用軟件洪澇監測專題分析工具，首先對研究區域的高分三號SAR數據進行預處理(輻射定標、多視處理、圖像濾波、地形校正)，得到散射系數強度局部圖；利用地形指數計算研究區陰影，由于研究區域地勢較為平坦，判識結果為該地區未產生陰影。最終生成HH和HV的Tiff格式產品，本研究選用HV通道，通過人機交互判識確定閾值，得到研究區的水體分布及變化情況。

2.1 鄭州及其附近地區洪澇遙感監測結果

將鄭州及其附近地區暴雨開始前的7月15日和暴雨結束后的7月21日兩幅SAR影像分別進行處理，設置閾值提取水體信息，得到暴雨前后研究區內水體變化情況(圖1和圖2)。在影像中，各種地物根據各自后向散射強度大小以不同灰度顯示，其中水體由于表面光滑后向散射最小，在圖中以黑色顯示，包括黃河、賈魯河、城郊水庫及暴雨后產生的洪澇水體。

圖1中偏上部黑色東西向河流為黃河，鄭州位于黃河以南，城區由于建筑物較多，后向散射強度大，在圖像中以淺色顯示。受地勢影響，城區積水無法排入黃河，只能通過賈魯河下泄。賈魯河發源于城區西南部丘陵地帶，繞過城北流向東南方，在河南省周口市匯入沙潁河，最終匯入淮河。賈魯河河道曲折，流速緩慢，這對鄭州及其附近地區及時排澇造成不利影響。

圖1 2021年7月15日河南鄭州及周邊區域高分衛星監測圖

圖2 2021年7月21日河南鄭州及周邊區域高分衛星監測圖

從圖2可見，7月21日暴雨后河南鄭州及周邊區域出現大范圍洪澇水體，其中城市建成區積水時間短，20日為集中降雨時段，21日大部積水已順利排出(地鐵、隧道、施工工地及老舊小區除外)。在21日的高分三號圖像上，市中心已無明顯積水，與15日暴雨前相比變化不大，但黃河沿岸、賈魯河沿岸水體范圍有了明顯擴大，城區以東平原地區也出現了多處水體。

為了更精確地反映暴雨前后研究區水體變化情況，將7月15日和7月21日兩幅高分三號雷達衛星監測圖進行處理，將暴雨后增加的水體(即洪澇水體)用紅色顯示(見圖3)。賈魯河是鄭州市主要的排水通道，城區內多條河流、管道均流入賈魯河，由于河道排水不暢，致使沿岸地區積水較重；城市以東地區地勢低平，積水明顯多于其他地區的，又以潮河匯入東風渠附近最為嚴重(圖3中紅框1所示)。同時，賈魯河流經尉氏縣境內發生決口，造成大量農田淹沒(圖3中紅框2所示)。

為便于開展災情評估，將暴雨前后水體變化與研究區土地利用圖疊加，分析淹沒范圍內土地利用類型。結果顯示(圖4)，鄭州及周邊區域洪澇水體面積約95 km2，其中淹沒耕地面積約66 km2，淹沒城鄉、工礦和居民用地面積約28 km2。淹沒的城鄉、工礦和居民用地呈點狀零散分布在鄭州城市郊區及東部平原地區，耕地主要集中在尉氏縣、中牟縣、祥符區交界處。另外，由于受黃河上游降水影響，黃河水位有所上漲，黃河灘區有大量沙洲和沿岸土地被淹沒。

2.2 豫北部分地區洪澇遙感監測結果

7月21日至22日，豫北地區遭遇特大暴雨，新鄉市及周邊地區所降雨水均需要通過衛河下泄，衛河水位迅速上漲，沿岸多地被淹。本研究利用同樣的方法，選取暴雨前(圖5，7月21日)和暴雨后(圖6，7月22日)兩幅SAR影像開展水體監測，并分析暴雨造成的洪澇影響。

圖3 2021年7月15-21日河南鄭州及周邊區域洪澇變化高分衛星監測圖

圖4 2021年7月15-21日河南鄭州及周邊區域洪澇水體與土地利用疊加圖

圖5 2021年7月21日河南新鄉及周邊區域高分衛星監測圖

從圖5中可看出，衛河是研究區內主要的排水通道，流經豫北地區多個市縣，途中有多條河流匯入。衛河河道較窄，暴雨前(7月21日)河道較細，河道內及沿岸水體面積較小，在圖像中不易區分。

暴雨后(7月22日)水體遙感監測結果(圖6)顯示，在研究區內各種水體有了明顯擴大，主要分布在新鄉市衛河沿岸，其中共產主義渠合河站附近有多條河流交匯，加上地勢低洼，水位上漲后淹沒面積較大(圖7中紅框所示)。另外，新鄉市轄區內原有湖泊池塘等水體的面積均有明顯擴大。

圖6 2021年7月22日河南新鄉及周邊區域高分衛星監測圖

圖7 2021年7月21-22日河南新鄉及周邊區域洪澇變化高分衛星監測圖

經估算，從7月21日至22日暴雨過后，新鄉市及其周邊地區新增洪澇水體面積約121 km2(圖8)，主要受災區域以耕地為主，其中淹沒耕地面積約112 km2，城鄉、工礦和居民用地受災面積約6 km2，城區水體均及時排出，受影響較小。

圖8 2021年7月21-22日河南新鄉及周邊區域洪澇水體與土地利用疊加圖

3 小 結

河南省鄭州市和新鄉市分布于黃河兩岸，兩地所降暴雨產生的洪水均無法排入黃河，只能通過賈魯河和衛河排泄。由于這兩條河流河道狹窄，排水不暢，降雨區又位于兩條河流源頭附近，洪水行程長，導致兩地大面積長時間受淹，造成工農業嚴重損失。本研究在暴雨結束后的7月21日(鄭州)和7月22日(新鄉)立即采用高分三號星載合成孔徑雷達SAR數據開展水體識別及淹沒范圍估算。估算結果證明，SAR數據不受云層影響，可有效彌補可見光-近紅外波段衛星傳感器由于天氣原因造成的數據空缺，在暴雨過后及時準確提取地表洪澇水體信息，為防汛救災提供決策依據。由于耕地地物類型簡單，淹沒范圍大，利用高分三號SAR數據對耕地上水體的監測效果要好于對城鎮、工礦和居民用地上的水體監測效果，計算的淹沒面積、分布范圍更加準確。

此次強降雨造成鄭州市內多處街道積水。本研究試圖對城區內淹沒范圍進行提取，由于受到高分三號衛星空間分辨率限制，雖然可以監測到城市街道積水，但因水體范圍受混合像元影響而難以準確提取。同時，由于市中心排水及時，地表積水時間短、面積小，地鐵、隧道處積水無法監測，從選取的21日衛星影像上看，監測提取的城市洪澇水體主要集中在郊區與賈魯河沿岸，這符合城市排水系統分布及實際淹沒情況。今后在遙感監測基礎上，開展地面與交叉數據驗證，并結合驗證結果，修正模型參數，進一步提高遙感洪澇監測精度。