HỘI THẢO ỨNG DỤNG GIS TOÀN QUỐC 2011
302
THEO DÕIDIỄNTIẾNLŨLƯUVỰCSÔNGMEKONG
LÀM CƠSỞDỰBÁOLŨỞĐỒNGBẰNGSÔNGCỬULONG
SỬ DỤNGẢNHVIỄNTHÁMMODIS
(MONITORING THE PROGRESS OF FLOOD IN THE MEKONG RIVER BASIN AS A
BASIS FOR FLOOD FORECASTING IN MEKONG DELTA)
Phan Thanh Nhàn, Võ Quang Minh
Bộ môn Tài nguyên đất đai, Đại học Cần Thơ
Email: ptnhan18students@gmail.com
; vqminh@ctu.edu.vn
Abstract: This study used MODIS satellite image (MOD09) for monitoring the progress of flood
in Mekong River Basin and testing the flood forecasting method for the Mekong Delta in flooding
stage of 2011. The results showed that (80% reliability): MODIS image can be used to monitor
the progress of flood in large areas of the Mekong River Basin. There was a close relationship
between enhanced vegetation index EVI, land surface water index LSWI with growth status of
plants and surface water of the flood. Risk flood maps during the flood season of the study area
were established as the basis for developing the flood forecasting method applied to the Mekong
Delta. With accuracy about 91%, this flood forecast method weekly opened a new direction for
researching about environmental disasters using the resource of satellite image at low cost.
Therefore, it should use these images for monitoring the process, forecasting flood capability and
other related fields in combination with other types of vegetation indices.
Keywords: MODIS, EVI, LSWI, flood, Mekong River Basin, Mekong delta.
1. ĐẶT VẤN ĐỀ
Những năm gần đây do ảnh hưởng của biến đổi khí hậu, đã có nhiều trận lũ lớn xảy ra
tại vùng ĐồngbằngsôngCửuLong (ĐBSCL) mà phần lớn nguyên nhân là do lũ thượng
nguồn sôngMeKong đổ về. Do đó, các cơ quan dựbáo khí tượng thủy văn cần có những giải
pháp thiết thực, đưa ra những dựbáo về khả năng ngập lũ, giúp người dân có thể tránh hoặc
hạn chế những thiệt hại về người và tài sản. Công nghệ ảnhviễnthám là một trong những
công cụ có thể giúp theo dõi, giám sát hỗ trợ cho dựbáo ngập lũ. Hiện nay, có nhiều nghiên
cứu ứng dụng các mô hình dựbáo ngập lũở vùng ĐBSCL. Tuy nhiên, các phương pháp giám
sát và mô hình này đòi hỏi các số liệu và phương pháp tính toán phức tạp và bị giới hạn bởi
nhiều y
ếu tố như chu kỳ dòng chảy, hướng dòng chảy và lũ thượng nguồn. Với phương pháp
sử dụngảnhviễnthámMODIS đa phổ và đa thời gian của NASA kết hợp mô hình số độ cao
DEM (Digital Elevation Model) có thể giúp theodõi được diễntiến của lũ vùng hạ lưusông
MeKong làmcơsở cho dựbáolũ vùng ĐBSCL.
Nghiên cứu thực hiện với mục tiêu: (1) Đánh giá khả năng sửdụng ả
nh MODIS đa thời
gian theodõidiễntiếnlũ vùng lưuvựcsông Mekong; (2) Nghiên cứu đề xuất giải pháp dự
báo lũ cho vùng ĐBSCL giai đoạn mùa lũ năm 2011.
2. DỮ LIỆU VÀ PHƯƠNG PHÁP NGHIÊN CỨU
2.1. Dữ liệu
2.1.1. Dữ liệu ảnh MODIS/TERRA đa thời gian
Sử dụngảnh vệ tinh MODIS/TERRA SURFACE REFLECTANCE 8-DAY L3
HỘI THẢO ỨNG DỤNG GIS TOÀN QUỐC 2011
303
GLOBAL 500 M SIN GRID V005 (MOD09) với độ phân giải thời gian là 8 ngày và độ phân
giải không gian là 500m. Ảnh được thu thập từ hệ thống dữ liệu quan sát Trái đất của NASA
(EOS, 2006) từ năm 2009 đến nửa đầu năm 2011. Sửdụng các Band 1, 2, 3 và 6 (đỏ, hồng
ngoại gần và sóng ngắn hồng ngoại…) để tính toán các chỉ số EVI, LSWI và DVEL làmcơsở
cho việc phân tích, đánh giá kết quả.
2.1.2. Dữ liệu ảnh SPOT
Ảnh SPOT 5 có độ phân giải 80 x 80m, 3 band Red – Green – Blue thu thập từ trung
tâm viễnthám Singapore (CRIPS), được sửdụng để đối chiếu kiểm tra kết quả giải đoán từ
ảnh MODIS.
2.1.3. Bản đồ số độ cao DEM của khu vực nghiên cứu
Mô hình số độ cao DEM khu vực hạ lưusôngMeKong thu thập từ Trung Tâm Phân
tích dữ liệu trái đất ERSDAC của Nhật (http://www.gdem.aster.ersdac.or.jp/) dùng để đánh
giá hướng dòng chảy của các hệ thống sông trên những dạng địa hình của khu vực nghiên
cứu. Mô hình được tính toán trên cơsở các chỉ số về cao độ, độ dốc và hướng dòng chảy của
nước để định hướng dòng chảy của các sông.
2.2. Phương pháp nghiên cứu
2.2.1. Phương pháp thành lập bản đồ ngập lũ
Sử dụng phương pháp của Sakamoto
et al.,
(2007) đề nghị cho đánh giá lũ. Nếu giá trị
EVI ≥ 0,2 được coi như điểm ảnh đục che phủ bởi mây và sẽ được gỡ bỏ khỏi ảnh.
(Thenkabail
et al.,
2005; Xiao
et al,
,2006). Theo Xiao
et al.
(2005, 2006), nếu EVI > 0,3 đối
tượng được phân loại là pixel không lũ. Nếu EVI ≤ 0,05 và LSWI nhỏ hơn hoặc bằng 0, điểm
ảnh sẽ được xác định là pixel nước liên quan.
Hình 1: Phương pháp theodõidiễntiến và dựbáolũ cho ĐBSCL bằngảnh
MODIS (Nguồn: Sakamoto et al, 2007)
Ảnh chỉ số EVI, LSWI, DVEL của MOD09A1
Xây dựng bản đồ ngập lũ, theodõi về diễntiến và dựbáolũ cho khu
vực ĐBSCL năm 2011
ρ
BLUE
> 0,2
EVI ≤ 0,05
LSWI ≤ 0,0
EVI > 0,3
0,1< EVI ≤ 0,3 EVI ≤ 0,1
Không lũ
Pixel lẫn lộn
DVEL≤ 0,05
EVI ≤ 0,3
Các pixel liên quan đến nước
Chỉ số không Mây: EVI, LSWI, DVEL
N
gập suốt >120 ngày
Pixel lũ Sông, hồ, biển
SAI
Ảnh MOD09 phản xạ 8 ngày
HỘI THẢO ỨNG DỤNG GIS TOÀN QUỐC 2011
304
Sau đó cần phân loại nơi nào là pixel lũ hoặc hỗn hợp hoặc các đối tượng ngập nước dài
hạn. Nếu pixel liên quan đến nước có EVI < 0,1, được coi là pixel lũ. Nếu EVI > 0,1 nhưng <
0,3, pixel liên quan đến nước được xác định là điểm ảnh hỗn hợp. Những khu vực ngập nước
liên tục được tách ra từ các pixel lũ và hỗn hợp. Các pixel liên quan đến nước có thời gian
ngập lũ > 120 ngày sẽ được phân loại là đối tượng ngập nước dài hạn.
2.2.2. Phương pháp dựbáodiễn biến ngập lũ
1
Hệ số tăng (giảm) trung bình k = DT nguy cơ ngập tuần 1/ DT tăng (giảm) tuần 2 so với tuần 1 = 0.02 được
chọn là giá trị có tần suất xuất hiện cao nhất trong các kết quả tính toán k
1
, k
2
,k
3
…k
n
từ kết quả giải đoán các bộ
ảnh trong mùa lũ các năm 2009 đến 2010. Dấu (+) tương ứng giai đoạn lũ bắt đầu và dấu (-) tương ứng giai
đoạn kết thúc mùa lũ.
Để dựbáodiễn biến ngập lũ cho một tuần sau đó sửdụng một bộ ảnh chụp ở thời điểm
hiện tại, kết quả giải đoán được phân tích và đánh giá kết hợp với các dữ liệu độ đốc, hướng
dòng chảy khai thác từ DEM, dữ liệu khí tượng thủy văn để dựbáo xác định về diện tích và
sự phân bố
không gian ngập lũ của tuần sau đó. Kết quả dựbáo đưa ra ở tuần trước sẽ được
kiểm chứng bởi kết quả giải đoán bộ ảnh của tuần sau đó. Với kết quả giải đoán bộ ảnh hiện
tại về diện tích ngập lũ, diện tích vùng có nguy cơ ngập kết hợp nguồn dữ liệu khí tượng thủy
văn thu thập đưa ra những dựbáo cho diễn biến ngập lũ của tuần tiếp theo.
3. KẾT QUẢ THẢO LUẬN
3.1. Diễntiến ngập lũ vùng nghiên cứu
3.1.1. Sự phân bố không gian ngập vùng nghiên cứu
Trên cơsở phân loại các đối tượng bằngsự kết hợp các chỉ số EVI, LSWI và DVEL cho
được loạt ảnh thể hiện sự phân bố không gian ngập lũ xuyên suốt mùa lũ các năm 2009, 2010
và 2011. Kết quả so sánh diện tích ngập nước ở các năm cho thấy được xu thế diễn biến ngập
lũ làmcơsở cho việc xác định các thời điểm bắt đầu cósự gia tăng một cách đột ngột (lũ xuất
hiện) về diện tích ngập, tăng đến mức tối đa (lũ đạt đỉnh) và thời điểm bắt đầu giảm (lũ rút).
k
1
Kết quả giải đoán lũ vùng hạ lưusôngMeKong tuần 1
DT ngập lũ x DT có nguy cơ ngập
Diện tích n
g
ậ
p
lũdựbáo cho tu
ầ
n 2
+
Ki
ể
m chứn
g
Hình 2: Phương pháp dựbáo khả năng ngập lũ khu vực ĐBSCL
Dữ liệu KTTV
Kết quả giải đoán lũ vùng hạ lưusôngMeKong tuần 2
DT ngập lũ
Diện tích ngập lũdựbáo cho tuần 3
+
Kiểm chứng
Dữ li
ệ
u KTTV
Kết quả giải đoán lũ vùng hạ lưusôngMeKong tuần 3
k x DT có nguy cơ ngập
HỘI THẢO ỨNG DỤNG GIS TOÀN QUỐC 2011
305
Hình 3: Diễn biến ngập vùng lưuvựcsôngMeKong qua các năm
Kết quả giải đoán cho thấy vào thời gian từ 20/7/2009 lũ đã tiến về các tỉnh đầu nguồn
ĐBSCL. Diễn biến cùng thời điểm đối với năm 2010 là ngày 21/8 và năm 2011 là ngày 05/8.
Sau ngày 05/8/2011 lũ bắt đầu tăng nhanh do ảnh hưởng của mưa, bãoở Trung và Hạ lưu kết
hợp với triều cường sông MeKong. Lũ 2009 đạt đỉnh lũ trong giai đoạn nửa đầu tháng 10, đối
với năm 2010 thời điểm này vào khoảng cuối tháng 10 đến đầu tháng 11. Lũ năm 2011 tăng
nhanh và đạt đỉnh vào khoảng cuối tháng 10. Lũ 2009 và 2010 bắt đầu rút vào khoảng tuần
cuối tháng 11 hoặc đầu tháng 12 đến nữa cuối tháng 12 lũ rút hết hoàn toàn.
3.1.2. Sự thay đổidiện tích ngập nước vùng hạ lưusôngMeKong
Nhìn chung lũ các năm có diễ
n biến khá giống nhau, diện tích mặt nước khoảng 10.000
km
2
, tương đối ổn định trong mùa khô và bắt đầu có nhiều biến đổi từ đầu tháng 7, tuy nhiên
về mức độ ngập cósự khác biệt nhất định giữa các năm. Nếu như mùa lũ năm 2010 có dấu
hiệu giảm so với năm 2009 thì lũ năm 2011 có phần tăng lên bằng hoặc hơn năm 2009. Tại
ĐBSCL, lũ năm 2009 bắt đầu vào khoảng 20/7/2009 và đến đầu tháng 12/2009 thì rút dần,
tổ
ng diện tích ngập vào mùa lũ từ đầu tháng 7 đến hết tháng 12 ước đạt 460.254km
2
, đối với
năm 2010 vào khoảng trên 317.053 km
2
trong khi năm 2011 đạt 248.643 km
2
tính từ đầu mùa
lũ đến cuối tháng 9 của năm và tăng dần đến cuối tháng 10.
2009
2010
201
21/
14/9
11/
9
16/10
14/20/7
11
9/11
16
14/9 21/8
16/10
8/10
1/1
13/
04/
2
1
2
1 2
3
3
19/12
2
7/
12
1
12/
HỘI THẢO ỨNG DỤNG GIS TOÀN QUỐC 2011
306
Hình 4: Biến đổidiện tích ngập nước vùng hạ lưu và một số tỉnh ĐBSCL qua các năm
3.1.3. So sánh kết quả giải đoán với số liệu quan trắc
Có sự tương quan khá tốt giữa mực nước thực đo tại trạm Tân Châu và Châu Đốc so với
diện tích ngập cùng thời điểm tại khu vực tỉnh An Giang. Kết quả kiểm chứng cho hệ số R
2
các năm lần lượt là 0,84 và 0,84 (năm 2009), R
2
=0,87 và R
2
=0,84 (năm 2010) và R
2
=0,95 và
R
2
=0,94 (năm 2011). Với hệ số R
2
khá cao chứng tỏ rằng hai giá trị này tỉ lệ thuận với nhau,
khi mực nước sông càng tăng cao thì diện tích ngập nước cũng sẽ tăng và ngược lại. Điều này
càng có ý nghĩa trong việc đánh giá độ sâu ngập thông qua diện tích ngập nước của khu vực
nghiên cứu tại một thời điểm nhất định.
3.1.4. Ngày bắt đầu, ngày kết thúc và thời gian ngập xuyên suốt của lũ
S
ự phân bố không gian của những ngày bắt đầu thay đổi từng năm. Ngày bắt đầu của
trận lũ năm 2009 sớm hơn so với trận lũ năm 2010 (Hình 6). Đối với năm 2011, ngày bắt đầu
dọc theosông Hậu và sông Mê Kông tương tự như năm 2009 tuy nhiên ở khu vực cận biên
của vùng đồngbằng ngập lũ, đặc biệt là khu vực ven biển gần vịnh Thái Lan, ngày bắt đầu
của năm 2011 sớm hơn so với năm 2010.
Hình 5. Tương quan giữa diện tích ngập tỉnh An Giang và mực nước thực đo
tại trạm Tân Châu và Châu Đốc năm 2011
HỘI THẢO ỨNG DỤNG GIS TOÀN QUỐC 2011
307
NĂM 2009 NĂM 2010 NĂM 2011
Hình 6: Ngày bắt đầu, ngày kết thúc, thời gian ngập suốt của lũ
Chu kỳ ngập các năm thể hiện thời gian kéo dài của lũ được xác định dựa vào ngày bắt
đầu và ngày kết thúc của lũ. Đối với trận lũ tương đối lớn như 2009 thì thời gian bắt đầu diễn
ra sớm hơn lũ năm 2010 và thời gian kết thúc cũng muộn hơn và do đó tình trạng ngập úng
kéo dài gây khó khăn trong tiêu nước và ảnh hưởng lịch thời vụ trong gieo sạ. Lũ năm 2011
có thời gian bắt đầu cũng khá sớm và theo ước đoán thì thời gian kết thúc cũng sẽ muộn hơn
năm 2010 đồng nghĩa với một diện tích khá lớn đất sản xuất nông nghiệp sẽ bị ngập úng.
Ngày
bắt đầu
ngập lũ
Ngày
kết thúc
ngập lũ
Ngập lũ
xuyên
suốt
HỘI THẢO ỨNG DỤNG GIS TOÀN QUỐC 2011
308
Hình 7: So sánh kết quả giải đoán từ hai ảnhMODIS và SPOT khu vực An Giang
(ngày 04/12/2010)
3.1.5. So sánh kết quả giải đoán ảnhMODIS với ảnh SPOT
Kết quả giải đoán ảnhMODIS cần phải được kiểm chứng với kết quả giải đoán từ ảnh
SPOT độ phân giải cao. Nghiên cứu chọn hai ảnh đại diện thuộc khu vực hai tỉnh An Giang
(ngày 04/12/2010) và Đồng Tháp (ngày 04/12/2009) để làmcơsởđối chiếu với ảnh MODIS.
Kết quả kiểm chứng ảnh chụp khu vực tỉnh An Giang được thể hiện ở các hình 7.
Kết quả so sánh hai loại ảnh bản khu vực tỉnh An Giang cho thấy mức độ tương thích
khá cao (khoảng 80%). Khu vực màu xanh đậm được xác định là ngập lũcó mức độ tương
thích lớn giữa hai loại ảnh. Tỉ lệ % diện tích sai lệch (chỉ cóở
ảnh MODIS hoặc SPOT) là
không đáng kể. Ngoài phần diện tích chung hiện diệnở hai loại ảnh phần diện tích chỉ có
riêng ởảnhMODIS là tương đốicao hơn so với phần diện tích lũ giải đoán chỉ cóởảnh
SPOT. Điều này phần nào được giải thích bởi sự khác biệt về độ phân giải cũng như phương
pháp phân loại ảnh là khác nhau dẫn đến những kết quả không giống nhau.
3.2. Kết quả dựbáo ngập lũ khu vựcĐồngBằngSôngCửuLong
Quá trình thực hiện được chia làm 11 đợt dựbáo tính từ ngày 04/7/2011 đến ngày
30/9/2011, tương ứng với 11 kết quả kiểm chứng từ thực tế. Kết quả tương quan giữa diện
tích ngập lũdựbáo và diện tích ngập thực tế giải đoán từ ảnhMODIS cho thấy hệ số tương
quan là khá cao (R
2
= 0,91). Hình 8 cho thấy mức độ sai lệch trong các dựbáo đưa ra là không
quá lớn và hoàn toàn có thể chấp nhận được. Điều này chứng tỏ rằng bằng việc áp dụng
phương pháp dựbáo từ ảnhMODIS kết hợp với các dữ liệu khác tương đối chính xác về diễn
biến lũ.
CHÚ THÍCH
Vùng không ngập ở hai ảnh SPOT và MODIS Vùng ngập lũ chỉ cóở SPOT
Vùng ngập lũ chung ở hai ảnh
Vùng ngập lũ chỉ cóởMODIS
HỘI THẢO ỨNG DỤNG GIS TOÀN QUỐC 2011
309
Hình 8: Tương quan và sai lệch giữa diện tích dựbáo và diện tích thực tế
4. KẾT LUẬN
Từ kết quả trên cho thấy khả năng ứng dụngảnhviễnthámMODIS (MOD09) trong
theo dõidiếntiếnlũ vùng lưuvựcsôngMeKong là khá tốt. Có mối liên hệ giữa đặc điểm ảnh
đến sự phân bố không gian và thời gian ngập lũ vùng hạ lưusông MeKong. Với độ chính xác
của kết quả giải đoán khá cao cho thấy khả năng có thể ứng dụngảnh MOD09 rộng rãi cho
việc giám sát, theodõidiễn tiến, dựbáolũ trong một khoảng thời gian dài và liên tiếp nhau
cho một khu vựccó phạm vi rộng lớn như quốc gia, vùng lãnh thổ. Với độ chính xác khoảng
91%, phương pháp dựbáocó thể được áp dụng cho dựbáo khả năng diễn biến ngập của lũ
khu vực ĐBSCL với giới hạn khả năng dựbáo về thời gian là một tuần.
Với nguồn ảnhMODIS không tốn nhiều chi phí như hiện nay cùng với việc hoàn thiện
dần phương pháp dựbáo thì đây có thể xem là một giải pháp hiệu quả trong theodõi và dự
báo, lũ lụt, thiên tai.
Tài liệu tham khảo
EOS. (2006), NASA earth observing system data gateway, from http://edcimswww.cr.usgs.gov
Sakamoto, T., Nguyen, N. V., Kotera, A., Ohno, H., Ishitsuka, N., and Yokozawa, M. (2007),
Detecting temporal changes in the extent of annual flooding within the Cambodia and the
Vietnamese Mekong Delta from MODIS time-series imagery, Remote Sensing of Environment,
109(3), 295-313.
Thenkabail, P. S., Schull, M., and Turral, H. (2005), Ganges and Indus river basin land use/land cover
(LULC) and irrigated area mapping using continuous streams of MODIS data, Remote Sensing of
Environment, 95, 317-341.
Xiao, X., Bole S., Liu J., Zhuang D. (2005), “Mapping paddy rice agriculture in southern China using
multi-temporal MODIS images”, Remote Sensing of Environment , 95, pp. 480–492.
Xiao, X., Boles, S., Frolking, S., Li, C., Bau, J. Y., and Salas, W. (2006). Mapping paddy rice
agriculture in South and Southeast Asia using multitemporal MODIS images. Remote Sensing of
Environment, 100, 95-113.
. HỘI THẢO ỨNG DỤNG GIS TOÀN QUỐC 2011 302 THEO DÕI DIỄN TIẾN LŨ LƯU VỰC SÔNG MEKONG LÀM CƠ SỞ DỰ BÁO LŨ Ở ĐỒNG BẰNG SÔNG CỬU LONG SỬ DỤNG ẢNH VIỄN THÁM MODIS (MONITORING THE. giúp theo dõi được diễn tiến của lũ vùng hạ lưu sông MeKong làm cơ sở cho dự báo lũ vùng ĐBSCL. Nghiên cứu thực hiện với mục tiêu: (1) Đánh giá khả năng sử dụng ả nh MODIS đa thời gian theo dõi. pháp theo dõi diễn tiến và dự báo lũ cho ĐBSCL bằng ảnh MODIS (Nguồn: Sakamoto et al, 2007) Ảnh chỉ số EVI, LSWI, DVEL của MOD09A1 Xây dựng bản đồ ngập lũ, theo dõi về diễn tiến và dự báo lũ