NGHIÊN CỨU & TRAO ĐỔI THỬ NGHIỆM ĐÁNH GIÁ TÁC ĐỘNG CỦA BIẾN ĐỔI KHÍ HẬU ĐẾN NGẬP LỤT KHU VỰC HẠ LƯU SÔNG LAM PGS.TS Nguyễn Thanh Sơn, PGS.TS Trần Ngọc Anh, ThS Đặng Đình Khá Trường Đại học Khoa học Tự nhiên, Đại học quốc gia Hà Nội ThS Nguyễn Xuân Tiến - Đài Khí tượng Thủy văn khu vực Bắc Trung Bộ CN Lê Viết Thìn - Viện Khoa học Khí tượng Thủy văn Biến đổi khí hậu B ài báo giới thiệu số kết đánh giá tác động biến đổi khí hậu (BĐKH) đến ngập lụt khu vực hạ lưu lưu vực sơng Lam Nhằm mơ tình trạng ngập lụt khu vực nghiên cứu, mơ hình thủy lực kết nối 1-2 chiều Mike Flood xây dựng với biên dịng chảy kết mơ từ mưa sử dụng mơ hình Mike NAM Hiện trạng ngập lụt khu vực hạ lưu tái thông qua tổ hợp kết mô ngập lụt từ trận lũ lịch sử (1978, 1988, 2010) Trên sở kịch phát thải trung bình Bộ Tài nguyên Môi trường khuyến nghị, biến đổi lượng mưa ba tháng lớn khu vực Bắc Trung Bộ mực nước biển dâng tương ứng Cửa Hội sử dụng làm kịch biên đầu vào cho hệ thống mơ hình xây dựng tương ứng với trận mưa điển hình lựa chọn mơ tình hình ngập lụt khu vực hạ lưu lưu vực sông Lam Kết tổ hợp đồ ngập lụt tương ứng cung cấp mô tả ngập lụt khu vực nghiên cứu điều kiện BĐKH sở so sánh với đồ trạng nhằm xác định tác động BĐKH đến diện tích ngập lụt, độ sâu ngập lụt thời gian ngập lụt Mở đầu Đối với nghiên cứu tác động BĐKH đến vấn đề ngập lụt, thơng thường có yếu tố nhấn mạnh gia tăng lượng mưa nước biển dâng (NBD) Các kịch NBD xây dựng cho vùng biển khác với mốc thời gian tương lai theo kịch BĐKH, nhiều trường hợp sử dụng trực tiếp làm biên đầu vào cho mơ hình thủy lực Riêng việc đánh giá gia tăng lượng mưa ảnh hưởng đến ngập lụt mơ hình thủy văn - thủy lực cần phải làm rõ tác động BĐKH việc so sánh trường hợp trạng (chưa xét đến BĐKH) trường hợp có BĐKH Trong thực tế, có số cách tiếp cận sau: a) sử dụng toán, làm đầu vào cho hệ thống mơ hình thủy văn – thủy lực; b) sử dụng trận mưa điển hình khứ (thường trận gây ngập lụt lịch sử khu vực nghiên cứu) để mô trạng ngập lụt (chưa xét đến BĐKH), sau định lượng hóa gia tăng lượng mưa từ kịch có sẵn theo ngun tắc thu phóng trận mưa điển hình ứng dụng mơ hình thủy văn - thủy lực nhằm thu tranh ngập lụt có tác động BĐKH; c) sử dụng trận mưa theo tần suất thiết kế (1%, 2%, 5%, 10%,…) tương ứng theo nhu cầu quy hoạch thiết kế cơng trình thu phóng tác động BĐKH trận mưa thiết kế lựa chọn ứng dụng mơ hình tương tự mơ hình khí hậu khu vực, ứng dụng kịch Cách tiếp cận c) có nhiều thuận lợi ứng phát thải khác nhau, với điều kiện ban đầu dụng trực tiếp vào công tác quy hoạch, thiết kế, mô trạng thái thời tiết để nhận phù hợp với nhà thực hành kết có chuỗi số liệu mưa ngày 100 năm tính phổ biến cao, việc sử dụng trận lũ phân tích, sử dụng để tính tốn giá trị mang thực xuất theo (cách tiếp cận b) cho phép tính đại diện lượng mưa ứng với thời đoạn hình dung tốt tác động trạng 10 năm 20 tương lai lưới tính tương lai có BĐKH Mặt khác, theo cách tiếp Người đọc phản biện: TS Nguyễn Kiên Dũng TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 09 - 2014 13 NGHIÊN CỨU & TRAO ĐỔI cận b) việc thu thập số liệu, liệu có nhiều thuận nhánh cấp I lợi dễ dàng kiểm chứng thực tế Theo cách tiếp Những sông nhánh lớn sơng Lam Nậm cận a) mang tính đồng cao, số liệu mưa nhận Mô, Huổi Nguyên, sông Hiếu, sông Giăng sông chi tiết cho ô lưới, phản ánh biến La Các sông đóng góp lượng dịng chảy đáng đổi tác động yếu tố cục đến phân bố kể vào sơng Lam mưa, có độ phân giải thời gian tốt (mưa ngày) Đối với lưu vực sông nhỏ, mưa phân bố tương đối đồng theo không gian việc sử dụng + Sông La hợp lưu hai nhánh Ngàn Phố Ngàn Sâu bắt nguồn từ vùng núi cao phía tây Hà Tĩnh, có tổng diện tích 3.210 km2 đổ vào hạ lưu sông Cả Chợ Tràng trận lũ (lịch sử theo tần suất thiết kế) làm sở + Sông Hiếu bắt nguồn từ dãy núi cao Phu Hoạt để mơ điều kiện trạng có xét đến có độ cao đỉnh núi 2.452 m huyện Quế Phong, BĐKH hợp lý dễ dàng ứng dụng thực Quỳ Châu, Nghĩa Đàn Tân Kỳ đổ vào sông Cả tiễn [1, 2] Đối với lưu vực sông lớn, việc lựa ngã ba Cây Chanh Diện tích tồn lưu vực chọn trận lũ thiết kế theo tần suất khó 5.340 km2, chiều dài sông 228 km khăn (do phải tính tốn tổ hợp lựa chọn điểm + Các sông nhánh lớn Nậm Mô, Huổi khống chế tính tần suất cho tồn lưu vực,…), Ngun, sơng Giăng có tổng lượng dịng chảy năm sử dụng trận lũ lịch sử chiếm tới 62,3% lượng dòng chảy năm tới Yên hết mức độ ngập lụt cho toàn lưu vực Thượng sơng Cả Việc sử dụng mơ hình khí hậu khu vực có nhiều ưu điểm lại hàm chứa tính bất định cao, đặc biệt đợt mưa lớn gây lũ thời đoạn mô dài (dự tính cho hàng thập kỷ sau) khó khăn việc so sánh với ngập lụt thời kỳ Do vậy, cần có hướng đánh giá khắc phục nhược điểm nêu Giới thiệu vùng nghiên cứu Lưu vực sông Lam vị trí từ 18015'50" đến 20010'30" vĩ độ Bắc, từ 103045'10" đến 105015'20" kinh độ Đơng Phía bắc giáp lưu vực sơng Chu, phía tây giáp lưu vực sơng Mêkơng, phía nam giáp lưu Hình Sơ đồ lưu vực sơng Lam vực sơng Gianh, phía đơng giáp biển Đơng Tổng diện tích lưu vực 27.200 km2, phần diện tích Việt Nam 17.730 km2, chiếm 65,2% diện tích lưu vực Diện tích thuộc Lào 9.470 km2 chiếm 34,8% diện tích lưu vực Dịng sơng Cả có chiều dài 531km, chảy qua Lào 170 km qua Nghệ An - Hà Tĩnh 361km (hình 1) [4] Hệ thống sơng Lam có mật độ lưới sơng 0,6 km/km2 Các sơng suối đổ vào dịng ngắn dốc bắt nguồn từ vùng núi cao tỉnh Xiêm Khoảng, Nghệ An, Hà Tĩnh Tổng số có 44 sơng 14 TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 09 - 2014 Mùa lũ lưu vực sông Cả hạ du từ tháng - 11, lũ lớn thường xuất vào tháng 9, 10 Thiết lập mơ hình Mike FLOOD Nhằm mô ngập lụt lưu vực sông Lam, mơ hình MIKE FLOOD lựa chọn [1, 2, 5, 6] Trên lưu vực nghiên cứu, mơ hình gồm có thành phần: mạng lưới sơng mơ mơ hình MIKE 11 với biên dịng chảy đầu vào gia nhập khu mô từ mưa mơ hình MIKE NAM, khu vực trũng ven sông NGHIÊN CỨU & TRAO ĐỔI khu vực hạ lưu, xuất dòng chảy tràn từ sơng ngập úng nội đồng mơ hình MIKE 21 sử dụng nhằm mơ dịng chảy chiều diện tích ngập lụt Cụ thể thành phần mô tả chi tiết a Xây dựng mạng lưới thủy lực chiều Mô hình thủy lực chiều MIKE 11 sử dụng để mô lại chế độ thủy động lực sơng Mạng lưới thủy lực sử dụng để tính tốn bao gồm sơng lưu vực với tổng chiều dài 461km bao gồm 292 mặt cắt thể bảng hình 2, Bảng Thông tin đặc trưng mạng thủy lực chiều Sông Cả Sông Giăng 225 22 Số mặt cắt 164 11 Sông Hiếu 105 48 Tân Kỳ Sông La Sông Lam Sông Ngàn Phố Sông Ngàn Sâu Sông Hào 13 35 30 25 6.25 12 23 16 14 Linh Cảm Ngã ba Chợ Tràng Sơn Diệm Hòa Duyệt Sơng La TT Tên sơng Chiều dài (km) Hình Bản đồ phân chia tiểu lưu vực vùng nghiên cứu Điểm đầu Cửa Rào Thác Muối Điểm cuối Nhập lưu vào sông Lam Nhập lưu vào sông Cả Nhập lưu vào sông Cả(Ngã ba Cây Chanh) Ngã ba Chợ Tràng Cửa Hội Linh Cảm Linh Cảm Sông Lam Hình Sơ đồ mạng thủy lực 1D vùng nghiên cứu b Xây dựng mạng lưới thủy lực hai chiều có thay đổi nhiều địa hình lưới tính nhỏ Mơ hình MIKE 21 sử dụng để tính tốn hơn, từ 30 – 100 m (hình 4) dịng chảy bãi ngập lũ, vùng tính tốn chiều Sau xây dựng mạng lưới thủy lực Mike vùng nghiên cứu xác định sở 11 Mike 21 nghiên cứu tiến hành Coupling cà đồ địa hình kết hợp số liệu điều tra khảo sát mạng lưới thủy lực chiều chiều, liên kết trận lũ lịch sử nhằm đảm bảo vùng tính tốn bao bên lựa chọn để kết nối mơ hình trùm vùng ngập lưu vực Từ đồ địa hình 1:10.000 Bộ Tài nguyên Môi trường cấp, c Hiệu chỉnh kiểm định mơ hình nghiên cứu tiến hành xây dựng lưới tính cho Kết hiệu chỉnh kiểm định mơ hình mưa miền tính chiều Khu vực nghiên cứu rời rạc dòng chảy MIKE NAM cho trận lũ 10/2010 hóa theo lưới phần tử hữu hạn (FEM) với kích thước 10/2013 trạm Nghĩa Khánh Sơn Diệm cạnh ô lưới từ 100 – 200 m cho khu vực có địa sơng Hiếu sơng Ngàn Phố cho kết tốt, hình tương đối phẳng, cịn với khu vực đánh giá theo tiêu Nash đạt 70% (hình 5, 6) TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 09 - 2014 15 NGHIÊN CỨU & TRAO ĐỔI Hình Miền tính chiều vùng nghiên cứu Hình So sánh đường trình lưu lượng tính tốn thực đo trạm Nghĩa Khánh (trận lũ 10/2010) Hình So sánh đường trình lưu lượng tính tốn thực đo trạm Sơn Diệm (trận lũ 10/2013) Mơ hình thủy lực hiệu chỉnh kiểm định với trận lũ 10/2010 10/2013 cho lưu vực sông Lam Số liệu dùng để hiệu chỉnh kiểm định mực nước thực đo hệ thống sông số liệu diện ngập chụp từ vệ tinh ngày 20/10/2010 UNOSAT công bố Trận lũ từ 14 - 25/10/2010 sử dụng để hiệu chỉnh thơng số mơ hình Kết so sánh mực nước tính tốn thực đo trạm Linh Cảm sông La Yên Thượng, Nam Đàn sông Cả tốt với tiêu Nash 75% (hình 7-10) Nghiên cứu tiến hành hiệu chỉnh thông số thủy lực bãi ngập lũ thơng qua việc so sánh diện ngập tính tốn diện ngập thu từ vệ tinh vào ngày 20/10/2010 Kết cho thấy giá trị tính tốn mơ hình phù hợp với giá trị quan trắc diện ngập độ sâu ngập lụt (hình 11) Hình So sánh đường trình mực nước Hình So sánh đường trình mực nước tính tốn thực đo trạm Dừa (trận lũ tính tốn thực đo trạm n Thượng 10/2010) (trận lũ 10/2010) 16 TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 09 - 2014 NGHIÊN CỨU & TRAO ĐỔI Hình So sánh đường q trình mực nước tính tốn Hình 10 So sánh đường q trình mực nước tính thực đo trạm Linh Cảm (trận lũ 10/2010) toán thực đo trạm Nam Đàn (trận lũ 10/2010) Hình 11 So sánh đường diện ngập tính tốn ảnh vệ tinh ngày 20/ X/2010 Mơ hình thủy lực kiểm định với trận lũ từ ngày 13 - 25/10/2013 cho kết tốt với tiêu Nash 70% (hình 12 -15) Kết hiệu chỉnh kiểm định mơ hình cho kết tốt với tiêu Nash đạt >70% Do vậy, sơ đồ mạng lưới thủy lực xây dựng cho lưu vực sông thông số hiệu chỉnh kiểm định đáng tin cậy việc sử dụng mơ hình để mơ kịch ngập lụt khứ tương lai tác động BĐKH NBD 12.sánh S sánh So đườn ng quátrình trình mực mực nước nh HìnhHình 12 So đường q nướctín tính th hực đo trạm trạ ạm Dừa (trận nlũlũ10/2013) 10/2013)) tốntốn thực đo Dừa (trận Hình h 13 13 So So sánh sánh h đường ttrìnhmực ước tính tốn Hình đườngq qtrình mựcnư nước tính tốn thực c đoYên trạm Yên (trận Thượng (trận thực đovàtại trạm Thượng lũg10/2013) 0/2013) lũ 10 Hình 14 So sánh đường q trình mực nước tính Hình 15 So sánh đường q trình mực nước tính tốn thực đo trạm Linh Cảm (trận lũ 10/2013) toán thực đo trạm Nam Đàn (trận lũ 10/2013) TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 09 - 2014 17 NGHIÊN CỨU & TRAO ĐỔI Xây dựng đồ ngập lụt khu vực nghiên cứu điều kiện trạng a Lựa chọn trận mưa tính tốn Nhằm đánh giá điều kiện ngập lụt trạng chưa có tác động BĐKH lưu vực, phân tích tài liệu trận mưa gây lũ điển hình lưu vực (bảng 2, hình 16), nhận thấy có trận mưa có phân bố rộng đặc trưng lưu vực (1978, 1998 2010) sử dụng làm trận mưa tính tốn Mặt khác, phân bố mưa khác nên diễn biến diện tích ngập lụt lưu vực tương đối khác tổ hợp để xây dựng đồ ngập lụt bao quát tối đa khu vực có nguy bị ngập lũ Bảng Tổng lượng mưa trận lũ lưu vực sông Lam Trạm Mường Xén n Nghĩa Khán nh Quỳ Châu Dừa Đơ Lương Hịa Duyệt Linh Cảm Cửa Hội 9/1978 135,3 421,5 178,1 809,2 957,9 554,2 660,8 516,0 9/1988 457,2 355,0 424,9 667,8 584,2 377,9 484,4 612,5 10/2002 10/2007 10/2010 10 0/2013 58,4 310,0 104,0 68,0 241,2 310,0 343,0 280,0 146,8 426,0 196,0 164,6 282,5 339,0 275,0 234,0 312,5 285,3 355,7 327,0 519,6 356,8 940,0 482,0 467,6 285,6 652,0 506,0 319,2 186,7 1089 212,0 Hình 16 Tổng lượng mưa trạm trận lũ điển hình 18 b Xây dựng đồ ngập lụt Từ đầu vào lượng mưa trạm mưa ứng với trận lũ năm 1978, 1988 2010, sử dụng mơ hình hiệu chỉnh kiểm định thu kết độ sâu ngập lụt dạng ASCII, sau xử lý phần mềm Mapinfo để xây dựng vùng ngập lụt với độ sâu khác cho trận lũ lớp thông tin GIS (hình 17-19) Trên sở đồ ngập lụt có, tiến hành tổ hợp phương pháp chồng xếp đồ để tạo đồ ngập lụt tối đa cho khu vực lưu vực sơng Lam (hình 20) trạng ngập lụt diễn lưu vực điều kiện chưa có tác động BĐKH Hình 17 Bản đồ ngập lụt hạ lưu lưu vực sông Lam trận lũ 9/1978 Hình 18 Bản đồ ngập lụt hạ lưu lưu vực sơng Lam trận lũ 10/1988 TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 09 - 2014 NGHIÊN CỨU & TRAO ĐỔI Hình 19 Bản đồ ngập lụt hạ lưu lưu vực sơng Lam trận lũ 10/2010 Hình 20 Bản đồ ngập lụt hạ lưu lưu vực sông Lam điều kiện trạng chưa có tác động Xây dựng đồ ngập lụt khu vực nghiên cứu điều kiện BĐKH & NBD a Kịch BĐKH NBD Theo Kịch BĐKH&NBD năm 2012 [6] khuyến nghị Bộ Tài ngun Mơi trường kịch phát thải trung bình B2 sử dụng để triển khai, xây dựng thực giải pháp ứng phó với BĐKH&NBD Theo kịch phát thải B2, lượng mưa tỉnh Nghệ An – Hà Tĩnh qua thập kỷ vào mùa khô (tháng 3, 4, 5) có xu hướng giảm dần qua thập kỷ, ngược lại vào mùa mưa (tháng 9, 10, 11) lại có xu hướng tăng Tuy nhiên, tốc độ tăng mùa mưa nhanh so với tốc độ giảm mùa mưa Vào cuối kỷ 21, mức giảm mùa mưa giảm khoảng 8,5%, cịn mức tăng mùa mưa 10,8% [6] Giá trị lượng mưa tăng lớn mùa mưa thời đoạn sử dụng để tính tốn khả ngập lụt cho lưu vực sơng Lam Ngồi giá trị mực nước biển dâng thời đoạn vùng biển Cửa Hội đưa vào tính tốn nhằm đánh giá khả bất lợi xảy tương lai b Xây dựng đồ ngập lụt 22.Bản Bản ngập lụtlũtrận lũ 10/1988 21.đồ Bản đồ ngập trận 9/1978 Hình 22 đồđồ ngập lụt trận 10/1988 Hình 21.Hình Bản ngập lụt lụt trận lũlũ9/1978 Hình tác động BĐKH NBD đến năm 2100 tác động BĐKH NBD đến năm 2100 tác động BĐKH NBD đến năm 2100 tác động BĐKH NBD đến năm 2100 đồ ngập trận 10/2010 Hình 24 đồđồ ngập lụt hạ vực sông 24.Bản Bản ngập lụtlưu hạlưu lưu lưu vực sơng Lam Hình 23.Hình Bản23 đồBản ngập lụt lụt trận lũlũ10/2010 Hình tác động BĐKH NBD đến năm 2100 Lam tác động BĐKH NBD đến năm tác động BĐKH NBD đến năm 2100 tác động BĐKH NBD đến năm 2100 2100 TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 09 - 2014 19 NGHIÊN CỨU & TRAO ĐỔI Các trận mưa khứ tương ứng với trận lũ lịch sử mô sử dụng để thu phóng với điều kiện gia tăng tác động BĐKH đến 2100 theo kịch B2 đưa vào mơ hình thủy văn thủy lực có Kết mơ ngập lụt xây dựng thành đồ tương ứng (hình 21-23), sau tiến hành tổ hợp thu tranh ngập lụt khu vực hạ lưu tác động BĐKH NBD (hình 24) So sánh với điều kiện ngập lụt trạng (hình 20) nhận thấy nhìn chung ngập lụt gia tăng đáng kể diện tích ngập lụt độ sâu ngập lụt tương ứng Các tính tốn chi tiết từ đồ cho thấy diện tích ngập trạng 109.368 tăng lên BĐKH&NBD thành 115.331 (5,45%) so với nguy ngập trạng chủ yếu tập trung khu vực huyện Hưng Nguyên Nghệ An huyện Đức Thọ tỉnh Hà Tĩnh Một số vùng thuộc Hưng Nguyên độ sâu ngập lụt gia tăng đáng kể từ 0,5 m điều kiện trạng đến 1,5 m có tác động BĐKH Kết luận Các kết nghiên cứu đề nghị phương pháp sử dụng tổ hợp ngập lụt trận lũ lịch sử điển hình lưu vực làm điều kiện trạng (nền) sở mơ tác động BĐKH thông qua gia tăng lượng mưa tương ứng (thu phóng từ trận mưa thực gây lũ) NBD để xây dựng đồ ngập lụt điều kiện có BĐKH&NBD Cách làm có sở khoa học ứng dụng dễ dàng thực tiễn có tiềm ứng dụng cho lưu vực sơng lớn khác có điều kiện tương tự Lời cảm ơn: Nghiên cứu nhận hỗ trợ số liệu tài đề tài BĐKH-19 “Đánh giá mức độ tổn thương kinh tế - xã hội lũ lụt số lưu vực sơng miền Trung bối cảnh biến đổi khí hậu khai thác cơng trình thủy điện, thủy lợi” thuộc Chương trình “Khoa học cơng nghệ phục vụ Chương trình mục tiêu quốc gia ứng phó với biến đổi khí hậu” Các tác giả trân trọng cám ơn hỗ trợ quý giá Tài liệu tham khảo Trần Ngọc Anh (2011) Xây dựng đồ ngập lụt hạ lưu sông Bến Hải Thạch Hãn, tỉnh Quảng Trị Tạp chí Khoa học ĐHQGHN, Khoa học Tự nhiên Công nghệ Tập 27, số 1S, tr 1-8 Hồng Thái Bình, Trần Ngọc Anh Đặng Đình Khá (2010) Ứng dụng mơ hình MIKE FLOOD tính tốn ngập lụt hệ thống sơng Nhật Lệ tỉnh Quảng Bình, Tạp chí Khoa học ĐHQGHN, Khoa học Tự nhiên Công nghệ 26, số 3S,285-294 Bộ Tài nguyên Môi trường (2012) Kịch biến đổi khí hậu, nước biển dâng cho Việt Nam Viện Quy hoạch Thủy lợi (2006) Quy hoạch sử dụng tổng hợp nguồn nước lưu vực sông Cả”, Báo cáo tổng hợp Denmark Hydraulic Institute (DHI), 2007, “MIKE FLOOD Reference Manual”, DHI Denmark Hydraulic Institute (DHI), 2007, “MIKE FLOOD User Guide”,DHI, 514 pp 20 TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 09 - 2014 ... khảo Trần Ngọc Anh (2011) Xây dựng đồ ngập lụt hạ lưu sông Bến Hải Thạch Hãn, tỉnh Quảng Trị Tạp chí Khoa học ĐHQGHN, Khoa học Tự nhiên Công nghệ Tập 27, số 1S, tr 1-8 Hồng Thái Bình, Trần Ngọc. .. Khoa học ĐHQGHN, Khoa học Tự nhiên Công nghệ Tập 27, số 1S, tr 1-8 Hồng Thái Bình, Trần Ngọc Anh Đặng Đình Khá (2010) Ứng dụng mơ hình MIKE FLOOD tính tốn ngập lụt hệ thống sơng Nhật Lệ tỉnh Quảng... Chu, phía tây giáp lưu vực sơng Mêkơng, phía nam giáp lưu Hình Sơ đồ lưu vực sơng Lam vực sơng Gianh, phía đơng giáp biển Đơng Tổng diện tích lưu vực 27.200 km2, phần diện tích Việt Nam 17.730