BÀI BÁO KHOA HỌC NGHIÊN CỨU ĐÁNH GIÁ LƯỢNG DÒNG CHẢY SÔNG ĐÀ TỪ TRUNG QUỐC VÀO VIỆT NAM PHỤC VỤ CHO BÀI TOÁN QUY HOẠCH VÀ QUẢN LÝ TÀI NGUYÊN NƯỚC LƯU VỰC SÔNG ĐÀ Trương Vân Anh1, Nguyễn Thu Hiền1, Đặng Quốc Khánh2 Tóm tắt: Do mạng lưới trạm khí tượng thủy văn nước ta khơng đủ dày để kiểm soát lượng dịng chảy biên giới đổ vào lãnh thổ Việt Nam nên việc tính tốn dịng chảy ngồi biên giới lưu vực sơng suối lớn có ý nghĩa đặc biệt quan trọng cho toán quy hoạch quản lý tài nguyên nước lưu vực sông liên quốc gia Trong nghiên cứu này, phương pháp thống kê kết hợp với mơ hình tốn sử dụng để tính tốn lượng dịng chảy sơng Đà từ lãnh thổ Trung Quốc vào Việt Nam khơng có số liệu thực đo Kết cho thấy lượng dịng chảy sơng Đà phía Trung Quốc chiếm khoảng 40 - 80% Lai Châu, 25 - 65% Tạ Bú 22 - 55% Hịa Bình Con số có ý nghĩa cho việc quản lý tổng hợp nguồn nước sơng Đà nói riêng tồn hệ thống sơng Hồng nói chung sơng Đà chiếm đến 37% tổng lượng dịng chảy sơng Hồng tại Sơn Tây Từ khóa: Lưu vực sông liên quốc gia, quy hoạch quản lý tài nguyên nước, MIKE NAM, MIKE 11 -Muskingum, lưu vực sông Đà Ban Biên tập nhận bài: 24/4/2017 54 Mở đầu Quản lý tài nguyên nước lưu vực sông liên quốc gia vấn đề mang nhiều thách thức Nguyên nhân khơng có hợp tác thực chất quốc gia liên quan trình khai thác sử dụng nguồn nước lưu vực sơng Các quốc gia thượng nguồn ln có ưu việc khai thác sử dụng nguồn nước; quốc gia hạ du lại phải đối mặt với cạn kiệt, suy thoát nguồn nước không chủ động biết nguồn nước sẵn có để lên kế hoạch sử dụng cách hợp lý Trong số 13 lưu vực sông lớn Việt Nam có đến 10 lưu vực sơng liên quốc gia với Lào, Campuchia, Thái Lan Trung Quốc với phần diện tích nằm ngồi lãnh thổ Việt Nam lớn gấp 3,3 lần diện tích lưu vực nước [1] Theo ước tính, tổng lượng nước lãnh thổ Việt Nam khoảng 830 - 840 tỷ m3/năm có khoảng 310 - 330 tỷ m3 (khoảng 37%) nước Khoa Khí tượng - Thủy văn, Đại học Tài nguyên Môi trường Hà Nội Tạp chí Khí tượng Thủy văn Email: truongvananh.vn@gmail.com TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 06 - 2017 Ngày phản biện xong: 15/5/2017 nội sinh, phần lại chảy từ địa phận nước vào Việt Nam [2] Do mạng lưới trạm khí tượng thủy văn khơng đủ để quan trắc nguồn nước đổ vào lãnh thổ Việt Nam với việc thiếu/hoặc khơng có hợp tác quốc gia chia sẻ nguồn nước dẫn đến khó khăn việc quy hoạch quản lý tổng hợp tài nguyên nước lưu vực sông Cùng với phát triển khoa học cơng nghệ, có nhiều phương pháp giúp tính tốn nguồn nước phần lưu vực khơng có số liệu đo đạc ứng dụng cơng nghệ viễn thám để tính tốn biến đầu vào (mưa bốc hơi) cho mơ hình tập trung mưa - dòng chảy [3], hay ứng dụng GIS để tính tốn biến đầu vào (sử dụng đất) cho mơ hình phân bố mưa - dịng chảy [4] Tuy nhiên nghiên cứu thích hợp với lưu vực tự nhiên, chưa có cơng trình khai thác lớn sơng Ở lưu vực sông liên quốc gia lớn nước ta sông Hồng, sơng Mekong, phía thượng nguồn Trung Quốc xây dựng nhiều hồ chứa để phân bổ lại nguồn nước Quy mơ quy trình vận hành hồ chứa phía Việt Nam ta hầu BÀI BÁO KHOA HỌC khơng nắm được, việc sử dụng loại mơ hình mưa rào - dịng chảy khơng thích hợp Mục tiêu nghiên cứu đánh giá lượng dòng chảy từ Trung Quốc vào Việt Nam lưu vực sông Đà Để đạt mục tiêu đó, nghiên cứu sử dụng phương pháp thống kê đơn giản kết hợp với phương pháp mơ hình tốn Chuỗi số liệu thực đo trạm Lai Châu, Tạ Bú Hịa Bình giai đoạn trước xây dựng hệ thống hồ chứa lưu vực sơng Đà phía Việt Nam sử dụng để kiểm định lại độ xác phương pháp Kết sau phục vụ giải toán quy hoạch quản lý tài ngun nước lưu vực sơng Đà nói riêng tồn hệ thống sơng Hồng nói chung Tài liệu thu thập phương pháp nghiên cứu 2.1 Tổng quan khu vực nghiên cứu Sông Đà sông lớn hệ thống sông Hồng Sông chảy qua địa phận ba quốc gia Trung Quốc, Lào Việt Nam Với tổng diện tích lưu vực 52,900 km2, phần lãnh thổ phía Trung Quốc chiếm khoảng 47%, Lào chiếm khoảng 2% lại Việt Nam khoảng 51%, sơng Đà đóng góp khoảng 37% tổng lượng nước sơng Hồng Sơn Tây Tồn lưu vực trải dài từ 20040’N - 25000’N 100022’ 105024’W (Hình 1) Trên dịng sơng Đà (Phần lãnh thởViệt Nam), Lai Châu trạm đo lưu lượng thượng nguồn sơng, phía sau biên giới Việt Trung có diện tích khống chế 33,800 km2 26,800 km2 nằm địa phận Trung Quốc; tiếp đến trạm thủy văn Tạ Bú, với diện tích khống chế 45900 km2 trạm thủy văn Hịa Bình 51800 km2 Ngồi sơng nhánh Nậm Mức, Nậm Giàng, Nậm Mu có trạm thủy văn Nậm Mức, Nậm Giàng Bản Củng tương ứng Riêng trạm Bản Củng ngừng đo từ năm 1987 Như nhận thấy số liệu quan trắc đo đạc lưu vực khơng nhiều Hình Lưu vực sơng Đà 2.2 Phương pháp nghiên cứu Để tính tốn lượng dịng chảy sơng Đà phía Trung Quốc chảy Việt Nam, nghiên cứu sử dụng phương pháp thống kê đơn giản kết hợp với mơ hình mưa - dịng chảy mơ hình diễn tốn dịng chảy MIKE 11 - NAM MIKE 11 MUSKINGUM Theo đó, chênh lệch dịng chảy thực TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 06- 2017 55 BÀI BÁO KHOA HỌC đo trạm khống chế dọc sông Đà (Lai Châu, Tạ Bú Hịa Bình) lượng dịng chảy nội sinh lãnh thổ Việt Nam đóng góp trạm coi lượng dòng chảy đến từ Trung Quốc Để biểu diễn hàm toán học, ta gọi lưu lượng thực đo tương ứng Lai Châu, Tạ Bú Hịa Bình QtLC, QtTB, QtHB; lưu lượng đóng góp cho dịng chảy tương ứng Lai Châu, Tạ Bú Hịa Bình từ phần lãnh thổ Việt Nam Qt(LC_VN), Qt(TB_VN), Qt(HB_VN); lưu lượng đóng góp cho dịng chảy tương ứng Lai Châu, Tạ Bú Hịa Bình từ phần lãnh thổ Trung Quốc Qt(LC_TQ), Qt(TB_TQ), Qt(HB_TQ); lượng dịng chảy sơng Đà đóng góp từ lãnh thổ Trung Quốc QtTQ Từ ta có: ̴்ொ ܳ௧ ൌ ܳ௧̴ே ܳ௧ (1) ̴்்ொ (2) ܳ௧் ൌ ܳ௧̴்ே ܳ௧ (3) ு̴்ொ ு̴ே ܳ௧ு ൌ ܳ௧ (4) Do đó, chuỗi số liệu Q , Q , Q thu thập từ trạm đo thủy văn tương ứng; chuỗi số liệu Qt(LC_VN), Qt(TB_VN), Qt(HB_VN) tính tốn từ mơ hình mưa rào dịng chảy MIKE NAM mơ hình thủy văn diễn tốn dịng chảy MIKE 11 MUSKINGUM Theo đó, lưu vực sông Đà địa phận Việt Nam chia sáu tiểu lưu vực Dòng chảy cửa tiểu lưu vực tính tốn mơ hình MIKE NAM vàdiễn tốn dịng chảy cửa tiểu lưu vực đến sông Đà dọc dịng sơng Đà tính tốn mơ hình MIKE 11 MUSKINGUM 2.2.1 Giới thiệu mơ hình MIKE NAM MIKE NAM mơ hình nhận thức thơng số tập trung mưa - dòng chảy xem xét lưu vực đơn vị không gian với bốn bể chứa theo chiều thẳng đứng bao gồm bể chứa tuyết tan, bể chứa mặt, bể chứa tầng rễ bể chứa nước ngầm [7] Các thông số mô hình gói gọn thơng số thể đặc tính bể chứa trình bày bảng Mục tiêu hiệu chỉnh kiểm định mơ hình tìm tập hợp thơng số giúp tạo ்ொ ܳ௧ 56 ܳ௧ ̴்ொ ൎ ܳ௧ ̴்்ொ ൎ ܳ௧ ு̴்ொ ൎ ܳ௧ LC t TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 06 - 2017 TB t HB t phù hợp mô đo đạc Mơ hình cho phép xem xét bốn mục tiêu: i) phù hợp dịng chảy mơ thực đo; ii) phù hợp hình dạng đường trình thủy văn; iii) phù hợp đỉnh lũ theo thời gian, độ lớn tổng lượng iv) phù hợp cho dòng chảy kiệt Ở cần ý có đánh đổi mục tiêu khác Ví dụ, tìm thơng số tốt để mô đỉnh lũ lại mơ dịng chảy kiệt ngược lại Đối với tiểu lưu vực sơng khơng có số liệu đo đạc thủy văn mượn số liệu lưu vực tương tự Lưu vực tương tự xem xét dựa phân tích tương tự đặc tính lưu vực liên quan đến thơng số nêu Đầu tiên lưu vực nên có kiểu khí hậu Tất đặc tính khác lưu vực hiểu phản hồi thủy văn lưu vực [5] kích thước, hình dạng, mật độ sông suối, sử dụng đất, thổ nhưỡng 2.2.2 Giới thiệu mơ hình MIKE 11- MUSKINGUM Để diễn tốn dịng chảy đến trạm đo thủy văn dọc dịng sơng Đà, phương pháp diễn tốn động học (Kinematic Routing Method) Muskingum MIKE 11HD sử dụng.Phương pháp phương pháp diễn toán thủy văn sử dụng để kiểm soát mối quan hệ lưu lượng - dung tích trữ, mơ tổng lượng trữ lũ sông thông qua việc kết hợp trữ dạng nêm trữ dạng Trong suốt thời kỳ lên sóng lũ, dịng chảy đến lớn dòng chảy sản sinh trữ dạng nêm Trong thời kỳ lũ đi, dòng chảy lớn dịng chảy đến, tạo nêm âm Phương trình đặc trưng phương pháp là: ାଵ ାଵ (5) ܳ ൌ ܥଵ ܳ ܥଶ ܳ ܥଷ ܳ ܥସ ାଵ ାଵ Trong số i, j đại diện cho điểm lưới xem xét mức thời gian Các biến số, C1C4 cho bởi: ܥଵ ൌ ܥଷ ൌ ο௧ିଶ௫ ଶሺଵି௫ሻାο௧ ଶሺଵି௫ሻିο௧ ଶሺଵି௫ሻାο௧ ο௧ାଶ௫ (6) ܥଶ ൌ ଶሺଵି௫ሻାο௧ (8) ܥସ ൌ ଶሺଵି௫ሻାο௧ ொೌ ο௧ (7 (9 BÀI BÁO KHOA HỌC Theo phương pháp Muskingum, thông số cần xác định bao gồm K (thời gian chảy truyền) x (trọng số), cho không đổi theo không gian thời gian Kết thảo luận Như trình bày trên, lưu vực sông Đà địa phận Việt Nam chia thành tiểu lưu vực: Tả thượng Lai Châu, Hữu Thượng - Lai Châu, Nậm Mức, Nậm Mu, khu Lai Châu Tạ Bú, Tạ Bú - Hịa Bình với trạm thủy văn khống chế tương ứng sau: Lai Châu cho tiểu lưu vực (Tả - Hữu Thượng Lai Châu), Nậm Mức, Nậm Mu, Tạ Bú Hịa Bình Số liệu dịng chảy thực đo khống chế lưu vực có Nậm Mức, Bản Củng với trạm khống chế tương ứng Nậm Mức Nậm Mu Do thơng số mơ hình NAM xây dựng cho hai tiểu lưu vực có đủ số liệu Nậm Mức Nậm Mu thơng qua tính tốn hiệu chỉnh (giai đoạn 1974 - 1984) kiểm định (giai đoạn 1985 - 1987) Việc lựa chọn giai đoạn hiệu chỉnh kiểm định xa so với số liệu đo đạc cho trạm Bản Củng kết thúc năm 1987 Các lưu vực lại phân tích xem xét khả mượn thơng số lưu vực tương tự 3.1 Hiệu chỉnh kiểm định mơ hình MIKE NAM cho hai tiểu lưu vực đủ số liệu Nghiên cứu xây dựng hai mơ hình mơ tốt cho hai tiểu lưu vực có đủ số liệu quan trắc Nậm Mức Nậm Mu Giá trị số đánh giá khả mơ dịng chảy lưu vực thể bảng cho kết cao với phù hợp hai đường q trình lưu lượng thực đo tính tốn hai giai đoạn hiệu chỉnh kiểm định (Hình 3) Như kết luận mơ hình NAM sử dụng mơ dịng chảy lưu vực sơng Nậm Mức Nậm Mu với thông số kiểm định tương ứng cho lưu vực trình bày bảng Hình Các tiểu lưu vực lưu vực sơng Đà Bảng Giá trị số NASH mô hình NAM cho lưu vực Nậm Mức Nậm Mu ChӍ sӕ NASH Quá trình HiӋu chӍnh KiӇm ÿӏnh PBIAS RSR NұmMӭc Nұm Mu Nұm Mӭc Nұm Mu Nұm Mӭc Nұm Mu 0.84 0.89 4.12 -44.51 0.45 0.66 0.75 0.71 -3.00 -30.20 0.62 0.53 Bảng Bộ thông số mô hình MIKE NAM hai lưu vực sơng Nậm Mức Nậm Mu Thông sӕ Ĉѫn vӏ Mӕi liên hӋ vӟi ÿһc tính lѭu vӵc Nұm Mu Nұm Mӭc Umax Lmax CQOF CKIF CK1,2 TOF TIF TG CKBF CK2 Mm Mm [-] Giӡ Giӡ [-] [-] [-] Giӡ Ĉѫn Loҥi thҧm phӫ mұt ÿӝ thҧm phӫ Loҥi ÿҩt, loҥi thҧm phӫ mұt ÿӝ thҧm phӫ Loҥi ÿҩt, ÿӝ dӕc lѭu vӵc, thҧm phӫ, mұt ÿӝ lѭӟi sông Ĉӝ che phӫ bӅ mһt Ĉӝ c lѭudӕ vӵc, loҥi ÿҩt kích thѭӟc lѭu vӵc Loҥi ÿҩt thҧm phӫ Loҥi ÿҩt thҧm phӫ Loҥi ÿҩt thҧm phӫ Ĉӝ dӕc lѭu vӵc, kích thѭӟc lѭu vӵc, hình dҥng lѭu vӵc Mӕi liên hӋ vӟi ÿһc tính lѭu vӵc 1.06 100 0.638 860.4 24.1 0.171 0.00161 0.01 1200 24.1 138 0.85 201 34 0.88 0.28 0.04 1100 34.2 TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 06 - 2017 57 BÀI BÁO KHOA HỌC (a) (b) (c) (d) Hình a) đường q trình dịng chảy thực đo tính tốn giai đoạn hiệu chỉnh b) giai đoạn kiểm định lưu vực Nậm Mức; c) d) lưu vực sông Nậm Mu 58 3.2 Tính tốn dịng chảy tiểu lưu vực khơng có số liệu Dịng chảy bốn tiểu lưu vực cịn lại tính tốn cách mượn số liệu lưu vực tương tự So sánh điều kiện tự nhiên chọn trạm tương tự cho mô hình MIKE NAM (Bảng 3) cho thấy lưu vực Nậm Mu lưu vực tương tự với tiểu lưu vực tả Lai Châu; lưu vực Nậm Mức lưu vực tương tự Hữu Lai Châu, khu Lai Châu - Tạ Bú Tạ Bú - Hòa Bình Kiểm định độ xác tính tốn thực phần tính tốn dịng chảy đến sơng Đà từ Trung Quốc 3.3 Diễn tốn dịng chảy tiểu lưu vực trạm khống chế dọc sơng Đà Phần trình bày thiết lập mơ hình diễn tốn dịng chảy MIKE 11 Muskingum cho sông Đà Hệ thống chia thành hai đoạn: Đoạn đầu Lai Châu Tạ Bú với chiều dài 170 km Đoạn có hai nhập lưu khu TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 06 - 2017 Nậm Mức Nậm Mu với chiều dài tương ứng 23,6 km 38 km Độ dốc đáy sơng trình qn 0,038 - 0,04% hệ sống nhám Manning 0,035 Đoạn sau Tạ Bú Hồ Bình với chiều dài 167 km độ dốc đáy bình quân 0,06% Hệ số nhám Manning toàn đoạn 0,035 Như vậy, có bốn đoạn sơng tính diễn tốn dịng chảy bao gồm Lai Châu - Tạ Bú, Tạ Bú Hịa Bình, Nậm Mức - sơng Đà Bản Củng sông Đà chuỗi số liệu thực đo hai trạm khống chế Tạ Bú Hịa Bình sử dụng để hiệu chỉnh mơ hình Q trình hiệu chỉnh giúp xác định thơng số mơ hình diễn tốn dịng chảy Muskingum bao gồm: x trọng số K thời gian chảy truyền trung bình đoạn thơng qua phân tích biểu đồ dịng chảy thực đo tính tốn Mơ hình diễn toán đánh giá sử dụng hệ số tương quan R hệ số NASH -Sut- BÀI BÁO KHOA HỌC cliffe (EI) Giá trị hai hệ số gần chứng tỏ mơ hình xác ngược lại Hình Sơ đồ diễn tốn dịng chảy vùng nghiên cứu 3.4 Hiệu chỉnh kiểm định mơ hình diễn tốn dịng chảy MIKE 11 - MUSKINGUM cho hệ thống sông Đà Trong phần này, kết mơ so sánh với dịng chảy thực đo trạm khống chế Trong phương pháp Muskingum, có hai thơng số phân tích K x Hai thông số đặc trưng cho tốc độ hình thành đỉnh lũ hình dạng đường q trình dịng chảy Lưu lượng thực đo trạm khống chế Tạ Bú Hịa Bình sử dụng để hiệu chỉnh mơ hình giai đoạn 1974 -1984 sau: - Biên của mạng tính toán bao gồm Lai Châu, Bản Củng Nậm Mức biên khu mạng tính tốn bao gồm Lai Châu - Tạ Bú Tạ Bú - Hòa Bình: Lấy chuỗi số liệu tính tốn từ mơ hình MIKE NAM tương ứng cho tiểu lưu vực - Hiệu chỉnh thông số K x đoạn sơng mạng tính tốn đến hai trạm khống chế Tạ Bú Hịa Bình Kết thể hình bảng a) Hiệu chỉnh mơ hình: Mơ hình hiệu chỉnh cách so sánh kết tính tốn với ch̃i sớliệu thực đo thời kỳ 1974 -1984 hai trạm Tạ Bú Hòa Bình (Hình 5) Hệ sớtương quan (R2) Chỉ sớNash (EI) rất cao (Bảng 3) chứng tỏ mơ hình mơ tốt mạng tính tốn b) Kiểm định mơ hình:Với thơng sớ lựa chọn từ q trình hiệu chỉnh, tiến hành kiểm định lại mơ hình với sớ liệu dòng chảy thực đo năm 1986 với xuất trận lũ lớn so với trận lũ xảy giai đoạn kiểm định Kết kiểm định trình bày hình Mơ hình tính xác thời gian xuất đỉnh giá trị đỉnh lũ năm kiểm định Kết kiểm định mơ hình tốt thể giá trị hệ số tương quan số Nash hai trạm khống chế lớn 0,95 (Bảng 3), đường trình thực đo tính tốn khá sát chân, đỉnh lũ thời gian xuất đỉnh (Hình 6) Như vậy, mơ hình sử dụng diễn tốn dịng chảy dọc sơng Đà Hình Q trình dịng chảy thực đo (đỏ) tính tốn (xanh) Tạ Bú (a) Hịa Bình (b) giai đoạn hiệu chỉnh từ 1974 -1984 TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 06 - 2017 59 BÀI BÁO KHOA HỌC Hình Đường q trình dịng chảy thực đo (đỏ) tính tốn (xanh) trạm Tạ Bú (a) Hịa Bình (b) năm 1986 (a) (b) Bảng Chỉ số đánh giá giai đoạn hiệu chỉnh Giai ÿoҥn HiӋu chӍnh KiӇm ÿӏnh Trҥm khӕng chӃ Tҥ Bú Hịa Bình Tҥ Bú Hịa Bình 3.5 Tính tốn lượng dịng chảy đến từ Trung Quốc Như đề cập trên, khơng có số liệu từ Trung Quốc sử dụng Do dịng chảy phía Trung Quốc tính tốn từ kết diễn tốn dọc sơng Đà với giai đoạn trùng với giai đoạn hiệu chỉnh (1974 - 1984) Quy trình thực sau: Tại trạm kiểm soát Lai Châu, biết tiểu lưu vực thượng Lai Châu, khác thực đo mô gọi Qt(LC_TQ) Tại trạm đo lưu lượng Tạ Bú, khác thực đo mô gọi Qt(TB_TQ) Tại trạm đo lưu lượng Hịa Bình, khác thực R2 0.971 0.988 0.955 0.981 EI 0.981 0.973 0.988 0.961 đo mô gọi Qt(HB_TQ) So sánh nhận thấy tương đồng giá trị dạng đường trình Qt(LC_TQ), Qt(TB_TQ), Qt(HB_TQ) (Hình 7), sai khác nhỏ lỗi mơ hình sai số từ kết đo đạc, bỏ qua Lượng dịng chảy sơng Đà phía bên Trung Quốc đóng góp trạm khống chế Lai Châu, Tạ Bú Hòa Bình dọc dịng sơng Đà tính tốn tỷ lệ giá trị dịng chảy phía Trung Quốc QtTQ với lưu lượng dòng chảy trạm thể hình bảng Hình Sự khác lưu lượng thực đo (Qthực đo) lưu lượng đóng góp từ nội biên (Qnội sinh) đến trạm khống chế dọc sông Đà 60 TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 06 - 2017 BÀI BÁO KHOA HỌC lai chau ta bu % dịng chҧy ÿóng góp tҥi trҥm 0.9 hoa binh 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 10 11 12 Hình Tỷ lệ đóng góp (%) dịng chảy sơng Đà từ Trung Quốc đến ba trạm khống chế Lai Châu (xanh đen), Tạ Bú (xanh da trời) Hịa Bình (vàng) dọc sơng Đà Tháng Bảng Tỷ lệ đóng góp dịng chảy sơng Đà phía Trung Quốc đến trạm Lai Châu, Tạ Bú Hịa Bình Tháng Tӹ lӋ ÿóng góp dịng chҧy sơng Ĉà phía Trung Quӕc Lai Châu Tҥ Bú Hịa Bình 0.82 0.62 0.56 0.81 0.63 0.57 0.84 0.65 0.56 0.61 0.42 0.37 0.40 0.26 0.24 0.42 0.29 0.26 Kết phần trăm dịng chảy Trung quốc đóng góp cho trạm Hịa Bình thấp mùa lũ (26 - 48% từ tháng đến tháng 10) cao mùa khô (52 - 57% từ tháng 11 đến tháng 3) Tháng tháng giao mùa, dòng chảy đến từ Trung Quốc chiếm khoảng 38% Điều gây nên bất lợi quy hoạch quản lý nguồn nước cho hệ thống hồ chứa lưu vực sông Đà nói riêng tồn lưu vực sơng Hồng chung lũ xảy đến mùa mưa hạn hán mùa khô, đặc biệt Trung Quốc xây dựng nhiều bậc thang thủy điện thượng nguồn sông Đà Kết luận kiến nghị Phương pháp phân tích thống kê phương pháp mơ hình tốn sử dụng đểtính tốn lượng dịng chảy ngoại sinh đóng góp cho trạm khống chế dọc sơng Đà Theo lượng dịng chảy phía Trung Quốc đóng góp lớn trạm Lai Châu, khoàng từ 40 - 80% tống Tháng 10 11 12 Tӹ lӋ ÿóng góp dịng chҧy sơng Ĉà phía Trung Quӕc Lai Châu Tҥ Bú Hịa Bình 0.52 0.38 0.35 0.58 0.44 0.40 0.69 0.52 0.43 0.73 0.58 0.48 0.73 0.59 0.52 0.69 0.55 0.49 lượng dòng chảy đến Lai Châu nên việc quản lý vận hành hồ Lai Châu (hồ chứa bậc thang hồ chứa sơng Đà) gặp nhiều khó khăn Tiếp đến trạm Tạ Bú (nằm chân hồ Sơn La) chiếm khoảng từ 25 - 65% trạm Hịa Bình từ 25 - 55% Vấn đề cần đối mặt vào mùa khơ tỷ lệ dịng chảy phía bên Trung Quốc lại cao thấp mùa lũ Do việc cấp nước mùa khơ gặp nhiều khó khăn phía Trung Quốc vận hành loạt cơng trình điều tiết lại dịng chảy phía thượng nguồn Nghiên cứu làm tiền đề cho nghiên cứu quy hoạch quản lý lưu vực sơng Đà nói riêng lưu vực sơng Hồng nói chung Kết hiệu chỉnh kiểm định cho thấy mơ hình sử dụng nghiên cứu đáng tin cậy có thể áp dụng cho lưu vực sông liên quốc gia khác có điều kiện tương tự sơng Đà TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 06 - 2017 61 BÀI BÁO KHOA HỌC Tài liệu tham khảo Đào Trọng Tứ công sự, (2011), Tổ chức Quản lý lưu vực sông Việt Nam “Quyền lực Thách thức” Thanh Tâm, (2015), Nguồn nước Việt Nam mối liên quan đến nước có chung nguồn nước Irena Ymeti, (2007), Rainfall estimation by remote sensing for conceptual rainfall-runoff modeling in the Upper Blue Nile basin, International Institute for Geoinformation science and earth observation enschede, the Netherlands Manoj K.Jain, (2004), AGIS based distributed rainfall - runoff model, Elservier Publisher Howard S Wheater, Thorsten Wagener, Hoshin, (2004), Rainfall-runoff modeling for gauged and ungauged catchment, Imperial College Press DHI, (2011), MIKE 11 user guide DHI, (2011), MIKE NAM user guide D N Moriasiet all, (2007), Model evaluation guidelines for Systematic quantification of accuracy in watershed simulations, Trans ASABE, 50, 885-900 ESTIMATINGTHE FOREIGN FLOW FROM CHINA TO VIETNAM SUPPORTING WATER RESOURCES PLANNING AND MANAGEMENT IN DARIVER BASIN Truong Van Anh1, Nguyen Thu Hien1, Dang Quoc Khanh2 Hanoi University of Natural Resources and Environment Scientific and Technical Hydrometeorological Journal Abstract: Due to the meteo-hydrological monitoring network is not dense enough to observe the foreign flows in the transboundary river basins in Vietnam, the estimation of these flows is especially important for water resources planning and management.This research briefly describes a method for estimating Da River’s flow which is coming from China part The results show that this flow accounts for about 40 up to 80% in Lai Chau, similarly 25 to 65% in Ta Bu and approximately 22 to 55% in Hoa Binh These numbers have significant meaning for water resources planning and management in Da River basin in particular and Red river system in general Keywords: Transbounday river basins, water resources planning and management, MIKE NAM, MIKE 11 - Muskingum, Da River basin 62 TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 06 - 2017