Xâm nhập mặn sông La Hà Tĩnh 2013 sử dụng mô hình Mike 11. Trạm bến thủy, trạm Linh Cảm. Độ mặn, mùa kiệt. Các sông đổ ra biển biển luôn có sự tương tác với biển về động lượng cũng như vật chất, nhất là tại khu vực cửa sông. Sự xáo trộn giữa hai dong nước ngọt và mặn diễn ra đều đặn trong mỗi chu kỳ triều. Thời kỳ triều lên, nước mặn được đưa từ biển vào sâu trong sông và rút đi trong thời kỳ triều xuống tạo thành sự mặn hoá đều đặn trong không gian, theo thời gian dưới tác động của hai yếu tố cơ bản: lưu lượng nước ngọt từ nguồn xuống và thủy triều thể hiện qua biên độ và cường suất.
VIỆN HÀN LÂM KHOA HỌC VÀ CÔNG NGHỆ VIỆT NAM VIỆN KHOA HỌC VẬT LIỆU BÁO CÁO ĐỀ TÀI NGHIÊN CỨU KHOA HỌC CÔNG NGHỆ CẤP CƠ SỞ CHO CÁN BỘ TRẺ NĂM 2013 ĐỀ TÀI “Nghiên cứu thử nghiệm mơ hình MIKE 11 để dự báo xâm nhập mặn Sông La tỉnh Hà Tĩnh” Chủ nhiệm đề tài: CN Nguyễn Mạnh Dưỡng Đơn vị : Phịng Cơng nghệ Vật liệu Mơi trường Hà Nội, 12/2013 DANH MỤC CÁC TỪ VIẾT TẮT XNM: Xâm nhập mặn BĐKH: Biến đổi khí hậu NASH: Độ phù hợp q trình thực đo tính tốn Obs: Đo theo ca Q: Lưu lượng trạm đo H, S: Mực nước độ mặn trạm KT – XH: Kinh tế xã hội PSU: Độ mặn ‰ DANH MỤC HÌNH Hình Vị trí địa lý sơng La…………………………………………………… ……2 Hình Sơ đồ sai phân hữu hạn điểm ẩn Abbot…………………………… …….12 Hình Sơ đồ sai phân điểm ẩn Abbott mặt phẳng x~t……………… …….12 Hình Sơ đồ mạng lưới thủy lực 1D………………………………………… ……17 Hình Mạng lưới 1D Mike 11………………………………………… ……17 Hình Mặt cắt thể Mike 11……………………………………… …….18 Hình Mực nước tính tốn thực đo thời kỳ hiệu chỉnh…………… …….19 Hình Mực nước tính tốn thực đo theo số liệu quan trắc tháng 4/2008 trạm Bến Thủy……………………………………………………………….……….……19 Hình Mực nước tính tốn thực đo thời kỳ kiểm định năm 2009……… 20 Hình 10 Mực nước tính tốn thực đo theo số liệu quan trắc tháng 4/2008 trạm Bến Thủy …………………………………………… ………………………… …20 Hình 11 Mực nước tính tốn thực đo thời kỳ kiểm định năm 2011… … 21 Hình 12 Mực nước thực đo tính tốn 2011 trạm Bến Thủy………………… …21 Hình 13 Độ mặn thực đo tính tốn 2008 trạm Bến Thủy………………… ……22 Hình 14 Độ mặn thực đo tính tốn 2009 trạm Bến Thủy…………………… …23 Hình 15 Độ mặn thực đo tính tốn năm 2011 trạm Bến Thủy…………… …24 Hình 16.Trắc mặn dọc sơng Cả năm 2008……………………………………… …25 Hình 17 Trắc mặn dọc sơng Cả năm 2011………………………………………… 26 Hình 18: Độ mặn theo tính tốn trạm Chợ Tràng tháng năm 2011………… 27 Hình 19 Độ mặn tính tốn trạm Linh Cảm tháng năm 2011……………….….28 Hình 20 Diễn biến XNM Sông La theo kịch B2 (2020-2039 – 15cm)……… 31 Hình 21 Diễn biến XNM Sơng La theo kịch B2 (2040-2059 – 25cm)…… ….32 Hình 22 Diễn biến XNM Sông La theo kịch B2 (2060-2079 – 50cm)……… 33 Hình 23 Diễn biến XNM Sơng La theo kịch B2 (2080-2099 – 75cm)……… 34 Hình 24 Diễn biến XNM Sông La theo kịch B2 (2080-2099 – 100cm)…….…35 DANH MỤC BẢNG Bảng 1.Tóm tắt số mơ hình tốn thường sử dụng Việt Nam…… … …9 Bảng Đặc điểm mạng lưới thủy lực 1D ….………………………… ….16 Bảng Số liệu dùng mơ hình……………………………………… …18 Bảng Mực nước biển dâng theo kịch phát thải B2 [9]……………….….29 Bảng Diễn biến xâm nhập mặn lớn dọc sông La (Chợ Tràng đến Linh Cảm) theo kịch B2……………………………………………………………… 30 MỞ ĐẦU Các sông đổ biển biển ln có tương tác với biển động lượng vật chất, khu vực cửa sông Sự xáo trộn hai dong nước mặn diễn đặn chu kỳ triều Thời kỳ triều lên, nước mặn đưa từ biển vào sâu sông rút thời kỳ triều xuống tạo thành mặn hoá đặn không gian, theo thời gian tác động hai yếu tố bản: lưu lượng nước từ nguồn xuống thủy triều thể qua biên độ cường suất Diễn biến phức tạp xâm nhập mặn hệ thống sơng Cả nói chung, sơng La nói riêng vấn đề cần phải quan tâm diễn biến gắn liền với nhu cầu sử dụng nước cho mục tiêu kinh tế - xã hội khu vực Cần thiết phải nghiên cứu trình xâm nhập mặn với trợ giúp cơng cụ tính tốn cơng nghệ cao để bước khám phá quy luật chung đặc thù xâm nhập mặn vùng cửa sông, trước mắt áp dụng cho vùng sông La Kết nghiên cứu làm sở cho việc khai thác sử dụng tài nguyên nước mặt lợi hạn chế tác động có hại q trình xâm nhập mặn Xâm nhập mặn trình tự nhiên nên nắm quy luật dự báo q trình phục vụ cho việc lấy nước tuới theo mùa vụ trồng thời đoạn dài bố trí thời vụ gieo trồng hợp lý để hạn chế tối đa tác động xâm nhập mặn Đồng thời, mặn điều kiện thuận lợi cho khu nuôi trồng thuỷ sản, hệ sinh thái ngập nước ven sông Đây mục tiêu việc nghiên cứu xâm nhập triều mặn phục vụ hoạt động kinh tế-xã hội khu vực sông La CHƯƠNG I TỔNG QUAN KHU VỰC NGHIÊN CỨU VÀ TÌNH HÌNH XÂM NHẬP MẶN Đặc điểm địa lý tự nhiên khu vực 1.1 Vị trí địa lý Sông La phụ lưu sông Lam, dài khoảng 12,5 km chảy qua huyện Đức Thọ, tỉnh Hà Tĩnh Sông La hợp lưu hai sông Ngàn Phố (từ Hương Sơn) sông Ngàn Sâu (từ Hương Khê Vũ Quang) bến Tam Soa (Linh Cảm, Đức Thọ) Sau đó, sơng La hợp lưu với Sơng Cả (từ Nghệ An) đổ biển Cửa Hội Hình Vị trí địa lý sơng La 1.2 Địa hình Khu vực sơng La có dạng địa hình đồi núi, gò đồi ven trà sơn, thung lũng đồng với khơng gian hẹp Địa hình thấp dần từ Tây sang Đơng bị chia cắt mạnh, phía Tây Nam huyện Đức Thọ chủ yếu núi thoải chạy dọc ven trà sơn, vùng núi dốc vùng giáp địa hình hành huyện lân cận, xen địa hình đồi núi thung lũng nhỏ hẹp 1.3 Khí hậu - Chế độ nhiệt: nhiệt độ trung bình hàng năm từ 22 - 25 0C Trong năm, khí hậu chia thành mùa rõ rệt: + Mùa nắng: kéo dài từ tháng đến tháng 10, khí hậu khơ nóng, từ tháng đến tháng Nhiệt độ trung bình thàng từ 24,7 0C (tháng 4) đến 32,90 (tháng 6) Nhiệt độ cao lên đến 38,5 - 400C + Mùa mưa: kéo dài từ tháng 10 đến tháng năm sau, nhiệt độ trung bình tháng từ 18,30C (tháng 1) đến 21,80C (tháng 11) với nhiều ngày số khu vực có nhiệt độ 70C (tháng 11, 12) - Độ ẩm: nhìn chung độ ẩm khơng khí tương đối cao (trung bình khoảng từ 84 87%), độ ẩm trung bình cao khoảng 92 - 96% vào tháng 1, 2, độ ẩm trung bình thấp khoảng 55 - 70% vào tháng 6, 7, - Bốc hơi: độ bốc trung bình năm đạt 800mm Lượng bốc lớn thường xảy vào tháng với mức trung bình tháng đạt 180 - 200mm Tháng có lượng bốc nhỏ 27 - 34mm - Số nắng: 1.400 - 1.600 giờ/năm - Chế độ gió: tỉnh Hà Tĩnh nằm khu vực chịu ảnh hưởng Khí hậu nhiệt đới gió mùa gió Lào hướng gió mang tính chất phân mùa khơng rõ rệt số địa phương khác thuộc khu vực Bắc Trung Bộ Bắc Bộ Tốc độ gió trung bình đạt từ 1,7 - 2,4m/s Mùa Đơng hướng gió chủ đạo gió Tây Bắc đến gió Bắc Đơng Bắc, tần suất tổng cộng tới 50 - 60% Mùa hè hướng gió chủ đạo hướng nam, tần suất 40 - 50% Bão thường xuất tháng kết thúc vào tháng 11, 12 Bình quân năm có từ - bão có ảnh hưởng trực tiếp đến Hà Tĩnh - Chế độ mưa: Hà Tĩnh có lượng mưa lớn, trung bình 2.000 mm/năm, cá biệt có vùng lên đến 3.500mm/năm Những mưa lớn thường xảy thượng nguồn sông Ngàn Phố, Ngàn Sâu, Rào Trổ, khu vực Hoành Sơn, Kỳ Lạc - Kỳ Anh Lượng mưa thường phân bố không đồng năm: Mùa Đông Mùa Xuân lượng mưa thấp, đạt 25% lượng mưa hàng năm, chủ yếu mưa phùn kết hợp gió mùa Đơng Bắc Mưa lớn tập trung vào mùa Hạ mùa Thu; chiếm 75% lượng mưa năm, đặc biệt cuối mùa Thu thường xuất nhiều trận mưa lớn gây ngập lụt, sạt lở đất, lũ quét [8] 1.4 Chế độ thủy văn Dịng chảy sơng La hàng năm dao động theo mùa, phù hợp với mùa mưa mùa khơ có lưu lượng phân phối khơng đều.Mùa lũ bắt đầu vào tháng V đơi tháng IV kết thúc vào tháng IX tháng X, khoảng tháng (tháng VI-IX) lượng nước chiếm khoảng 75-85% lượng nước năm Mùa kiệt vào khoảng từ tháng XI đến tháng IV năm sau Giai đoạn ổn định mùa kiệt có lượng dòng chảy nhỏ, giai đoạn kéo dài tháng (tháng I-IV), lượng dòng chảy tháng 7-8% lượng dòng chảy năm, thời kỳ mặn có điều kiện xâm nhập sâu vào lục địa Sông Ngàn Sâu chi lưu sơng La, có chiều dài khoảng 131 km, bắt nguồn từ vùng núi Ông Giao Thừa núi Cũ Lân thuộc dãy núi Trường Sơn nằm địa bàn giáp ranh hai tỉnh Hà Tĩnh Quảng Bình Diện tích lưu vực khoảng 3.214 km2, độ cao trung bình 360 m, độ dốc trung bình 28,2%, mật độ sơng suối 0,87 km/km2 Tổng lượng nước trung bình nhiều năm đạt 6,15 km tương ứng với lưu lượng trung bình năm 195 m3/s modul dịng chảy năm 47 l/s km2 Sông Ngàn Phố phụ lưu sông La chảy chủ yếu địa phận huyện Hương Sơn, tỉnh Hà Tĩnh Sông Ngàn Phố bắt nguồn dòng suối nhỏ từ vùng núi Giăng Màn địa phận xã Sơn Hồng, Sơn Kim Sơn Kim huyện Hương Sơn, ven biên giớiViệt-Lào, độ cao khoảng 700 m Sông Ngàn Phố chảy gần theo hướng Tây-Đông tới ngã ba Tam Soa (xưa có bến Tam Soa), nơi giáp ranh xã Sơn Tân, Sơn Long (huyện Hương Sơn) với xã Trường Sơn, Tùng Ảnh (huyện Đức Thọ) Chiều dài khoảng 87 km Diện tích lưu vực 1.091 km², độ cao trung bình 331 m, độ dốc trung bình 25,2% Mật độ sông suối 0,91 km/km² Tổng lượng nước khoảng 1,40 km³ tương ứng với lưu lượng trung bình 45,6 m³/s Một số cơng trình nghiên cứu XNM Việt Nam Biến đổi khí hậu (BĐKH) phạm vi toàn cầu làm cho thiên tai Việt Nam ngày gia tăng số lượng, cường độ mức độ ảnh hưởng, ảnh hưởng lớn đến hoạt động sản xuất, phát triển kinh tế xã hội Lĩnh vực chịu tác động mạnh mẽ BĐKH nông nghiệp, thuỷ lợi, thuỷ sản, diêm nghiệp, lâm nghiệp, an ninh lương thực; vùng đồng dải ven biển mực nước biển dâng [5] Cụ thể sau: Nhiệt độ: Trong 50 năm qua (1958 – 2007), nhiệt độ Việt Nam tăng lên từ 0,5 -0,70C Nhiệt độ mùa đông tăng nhanh nhiệt độ mùa hè Các vùng khí hậu phía Bắc tăng nhanh phía Nam Nhiệt độ trung bình thập kỷ gần (19612000) cao trung bình năm thập kỷ trước Lượng mưa: Trên địa điểm, xu biến đổi lượng mưa hàng năm thập kỷ qua (1911-2000) không rõ rệt theo thời kỳ vùng khác Lượng mưa vùng phí Bắc giảm tăng vùng phía Nam Tính trung bình năm, lượng mưua 50 năm qua giảm khoảng 2% Khơng khí lạnh: Số đợt khơng khí lạnh ảnh hưởng tới Việt Nam giảm rõ rệt thập kỷ qua Tuy nhiên, biểu dị thường lại xuất mà gần đợt khơng khí lạnh gây rét đậm kéo dài 38 ngày tháng tháng năm 2008 Bão: Năm gần đây, bão có cường độ mạnh xuất nhiều Quỹ đạo bão có dấu hiệu dịch chuyển dần phía Nam mùa bão kết thúc muộn nhiều bão có đường dị thường Mực nước biển: Theo số liệu quan trắc trạm hải văn học ven biển Việt Nam cho thấy tốc độ dâng lên mực nước biển trung bình Việt Nam khoảng 3mm/năm tương đương với tốc độ tăng trung bình giới [9] Một vấn đề cấp thiết tượng xâm nhập mặn gia tăng nhiều nơi, làm thu hẹp diện tích đất canh tác ảnh hưởng lớn đến tài nguyên thiên nhiên, đặc biệt làm giảm trữ lượng nước nhạt đất gây nhiễm mặn nước sông Nhận thức rõ tầm quan trọng trình XNM, nhiều ngành sớm triển khai chương trình, dự án nghiên cứu tác động BĐKH, lồng ghép yếu tố BĐKH vào chiến lược, quy hoạch phát triển KT - XH, hình thành chương trình hành động thích ứng với biến đổi khí hậu Theo nghiên cứu Ngân hàng Thế giới, nước ta với bờ biển dài hai vùng đồng lớn, mực nước biển dâng cao từ 0,2 - 0,6 m, có từ 100.000 - 200.000 đất bị ngập làm thu hẹp diện tích sản xuất nông nghiệp Nếu nước biển dâng lên 1m làm ngập khoảng từ 0,3 đến 0,5 triệu Đồng sống Hồng năm lũ lớn khoảng 90% diện tích Đồng Sơng Cửu Long bị ngập từ - tháng, vào mùa khơ khoảng 70% diện tích bị xâm nhập mặn với nồng độ lớn 4g/l Ước tính Việt Nam khoảng triệu đất trồng lúa tổng số triệu nay, đe dọa nghiêm trọng đến an ninh lương thực quốc gia ảnh hưởng đến hàng chục triệu người dân [5] Nước biển dâng làm nước mặn xâm nhập sâu vào nội địa, cống hạ lưu ven sông khả lấy nước vào đồng ruộng Các thành phố ven biển bị ngập úng triều như: Thành phố Hồ Chí Minh, Cần Thơ, Cà Mau, Hải Phòng, Trà Vinh Khu vực thấp ven biển bị ngập triều gây mặn nặng như: Bến tre, Cà Mau Chế độ dịng chảy sơng suối thay đổi theo hướng bất lợi, cơng trình thuỷ lợi hoạt động điều kiện khác với thiết kế, làm cho lực phục vụ cơng trình giảm Riêng khu vực ven biển Hà Tĩnh, kịch biến đổi khí hậu (Bộ Tài nguyên Môi trường, 2009) bao gồm kịch biến đổi nhiệt độ, lượng mưa kịch nước biển dâng, Hà Tĩnh địa phương chịu nhiều biến động mạnh thời tiết Cụ thể, vào kỷ 21 nhiệt độ Hà Tĩnh tăng lên 1,5°C lượng mưa tăng lên 4,0% vào cuối kỷ 21 nhiệt độ tăng lên 2,8°C lượng mưa tăng lên 7,7% so với thời kỳ 1980-1999 Tương tự, mực nước biển dâng 10 Độ mặn ‰ Ngày Độ mặn ‰ Hình 14 Độ mặn thực đo tính tốn 2009 trạm Bến Thủy Độ mặn ‰ Ngày Ngày Độ mặn ‰ Hình 15 Độ mặn thực đo tính tốn 2011tại trạm Bến Thủy 28 Ngày Qua kết hiệu chỉnh kiểm định mơ hình cho thấy hợp lý thơng số nhám, từ ta kết luận mơ hình với mạng lưới thủy lực xây dựng thông số hiệu chỉnh kiểm định đủ tin cậy để tính tốn tồn lan truyền mặn 29 [meter] [PSU] 1-4-2008 07:00:00 32.0 20.0 30.0 18.0 28.0 16.0 26.0 12.0 24.0 38516 42800 18.0 41300 46500 51040 57466 48600 44901 52480 55200 58500 62070 64000 65595 66600 67800 69000 22.0 20.0 16.0 71300 2.0 72100 4.0 73327 8.0 60500 10.0 6.0 35800 14.0 14.0 0.0 12.0 -2.0 10.0 -4.0 8.0 -6.0 6.0 -8.0 4.0 -10.0 2.0 -12.0 SONG CA 73327 - 44901 -14.0 0.0 2000.0 4000.0 6000.0 8000.0 10000.0 12000.0 14000.0 16000.0 SONG CA 44901 - 38516 18000.0 20000.0 22000.0 24000.0 26000.0 28000.0 30000.0 32000.0 Hình 16.Trắc mặn dọc sông Cả năm 2008 Đường màu đỏ thể khuy ếch tán mặn sông Ranh giới mặn 1‰ năm 2008 tiến sâu vào sông Cả khoảng 29km tính từ Cửa Hội 34000.0 SONG CA 38516 - 36000.0 38000.0 [m] 0.0 Ranh giới mặn 4‰ năm 2008 tiến sâu vào sơng Cả khoảng 25km tính từ Cửa Hội [meter] [PSU] 1-4-2011 07:00:00 32.0 20.0 30.0 18.0 28.0 16.0 26.0 12.0 24.0 38516 42800 41300 44901 46500 51040 57466 48600 52480 55200 58500 62070 64000 65595 66600 67800 69000 22.0 20.0 18.0 16.0 71300 2.0 72100 4.0 73327 8.0 60500 10.0 6.0 35800 14.0 14.0 0.0 12.0 -2.0 10.0 -4.0 8.0 -6.0 6.0 -8.0 4.0 -10.0 2.0 -12.0 SONG CA 73327 - 44901 -14.0 0.0 2000.0 4000.0 6000.0 8000.0 10000.0 12000.0 14000.0 16000.0 SONG CA 44901 - 38516 18000.0 20000.0 22000.0 24000.0 26000.0 28000.0 30000.0 Hình 17 Trắc mặn dọc sông Cả năm 2011 Ranh giới mặn 1‰ năm 2011 tiến sâu vào sông Cả khoảng 30.5km tính từ Cửa Hội Ranh giới mặn 4‰ năm 2011 tiến sâu vào sơng Cả khoảng 25,8km tính từ Cửa Hội 32000.0 34000.0 SONG CA 38516 - 36000.0 38000.0 [m] 0.0 Hình 18: Độ mặn theo tính tốn trạm Chợ Tràng tháng năm 2011 Hình 19 Độ mặn tính tốn trạm Linh Cảm tháng năm 2011 Như tính năm 2011 độ mặn sơng La khơng đáng kể nước hồn tồn khai thác cấp cho nơng nghiệp ni trồng thủy sản 3.4.Tác động nước biển dâng đến xâm nhập sông Các kịch phát thải khí nhà kính lựa chọn để tính tốn, xây dựng kịch nước biển dâng cho Việt Nam kịch phát thải thấp (kịch B1), kịch phát thải trung bình nhóm kịch phát thải trung bình (kịch B2) kịch phát thải cao nhóm kịch phát thải cao (kịch A1FI) Kịch phát thải thấp (B1) mơ tả giới phát triển tương đối hồn hảo theo hướng phát thải khí nhà kính nhất, tốc độ tăng dân số thấp, cấu kinh tế thay đổi nhanh theo hướng dịch vụ thông tin, thỏa thuận quốc tế nhằm giảm thiểu phát thải khí nhà kính thực đầy đủ nghiêm túc phạm vi toàn cầu Tuy nhiên, với cấu kinh tế không đồng khu vực giới nay, cộng với nhận thức khác biến đổi khí hậu quan điểm khác nước phát triển nước phát triển, đàm phán quốc tế biến đổi khí hậu nhằm ổn định khí nhà kính, hạn chế mức độ gia tăng nhiệt độ mức 0C gặp nhiều trở ngại, kịch phát thải thấp có khả thành thực kỷ 21 Các kịch phát thải cao (A2, A1FI) mô tả giới không đồng quy mơ quy mơ tồn cầu, có tốc độ tăng dân số cao, chậm đổi công nghệ (A2) sử dụng tối đa lượng hóa thạch (A1FI) Đây kịch xấu mà nhân loại cẩn nghĩ đến Với nỗ lực phát triển cơng nghệ thân thiện với khí hậu, đàm phán giảm phát thải khí nhà kính, chung tay, chung sức toàn nhân loại “liên kết chống lại biến đổi khí hậu”, hi vọng kịch phát thải cao có khả xảy Hơn nữa, cịn nhiều điểm chưa chắn việc xác định kịch phát triển kinh tế - xã hội kèm theo lượng phát thải khí nhà kính tương lai Với tồn điểm chưa chắn kịch biến đổi khí hậu, nước biển dâng ứng với kịch phát thải khí nhà kính cận cận có mức độ tin cậy thấp so với kịch mức trung bình Vì lý nêu trên, kịch biến đổi khí hậu, nước biển dâng Việt Nam khuyến nghị sử dụng thời điểm kịch ứng với mức phát thải trung bình (B2) Bảng Mực nước biển dâng theo kịch phát thải B2 [9] Kịch 2020 2030 Các mốc thời gian kỷ 21 2040 2050 2060 2070 2080 2090 2100 Phát thải trung bình (B2) 7-8 11-12 15-17 2-24 25-31 31-38 36-47 42-55 49-64 Bảng Diễn biến xâm nhập mặn lớn dọc sông La (Chợ Tràng đến Linh Cảm) theo kịch B2 Khoảng cách từ Chợ Tràng đến Linh Cảm (km) 7,3 8,3 10 11 12 12,5 Độ mặn theo kịch (0/00) B2 0.15m 0.25m 0.50m 0.75m 1m 1,9 1,7 1,55 1,37 1,3 1,15 1,03 2,2 1,95 1,77 1,62 1,5 1,38 1,25 1,05 2,6 2,37 2,13 1,9 1,8 1,6 1,42 1,23 3,6 3,4 3,1 2,9 2,75 2,53 2,31 1,9 1,73 1,7 1,47 1,3 1,22 1,2 4,8 4,5 4,19 3,7 3,4 3,3 2,5 2,3 2,12 2,1 1,9 Sơng La Hình 20 Diễn biến XNM Sông La theo kịch B2 (2020-2039 – 15cm) Ranh giới mặn 1‰ tiến sâu vào sông La khoảng 6,2km Sơng Sơng La Hình 21 Diễn biến XNM Sông La theo kịch B2 (2040-2059 – 25cm) Ranh giới mặn 1‰ tiến sâu vào sông La khoảng 7,3km Sơng La Sơng La Hình 22 Diễn biến XNM Sơng La theo kịch B2 (2060-2079 – 50cm) Ranh giới mặn 1‰ tiến sâu vào sông La khoảng 8,7km Sông La Sơng La Ngàn Sâu Hình 23 Diễn biến XNM Sông La theo kịch B2 (2080-2099 – 75cm) Ranh giới mặn 1‰ tiến sâu vào tồn Sơng La tiến vào sông Ngàn Sâu khoảng 2km Sông La Ngàn Sâu Sơng Hình 24 Diễn biến XNM Sơng La theo kịch B2 (2080-2099 – 100cm) Ranh giới mặn 1‰ năm tiến sâu vào tồn sơng La vào sơng Ngàn Sâu Ranh giới mặn 4‰ tiến sâu vào sông La khoảng 3,7km KẾT LUẬN Qua kết cho thây việc áp dụng Mike 11 vào việc dự báo XNM sông khu vực ven biển cần thiết, dự báo tương lai qua có biện pháp khai thác phịng tránh XNM cách hợp lý hiệu Hiện nhiễm mặn sông La chưa đáng kể Xâm nhập mặn Sông La dự đốn kịch khác nhau, hình thành từ giá trị mực nước biển dâng theo kịch biến đổi khí hậu B2 lưu lượng thượng nguồn mùa kiệt suy giảm theo kịch sử dụng nước từ thượng nguồn Kết từ kịch cho thấy khoảng cách xâm nhập mặn khu vực sông Cả cụ thể Sông La vào tháng kiệt tương lai Đặc biệt kịch nước biển dâng 75cm 100cm, khoảng cách xâm nhập mặn 1‰ tiến vào toàn khu vực Sơng La tiến sâu vào sông Ngàn Sâu Kết nghiên cứu hữu ích cho việc qui hoạch tài nguyên nước khu vực Hà Tĩnh nói chung khu vực Đức Thọ nói riêng KIẾN NGHỊ Tiếp tục mở rộng nghiên cứu q trình XNM sơng vùng ven biển năm TÀI LIỆU THAM KHẢO http://daitudien.net/dia-li-hoc/dia-li-hoc-ve-song-ca-nghe-an.html Nguyễn Tất Đắc, Nguyễn Văn Điệp, Nguyễn Minh Sơn, 1988, “Mơ hình tính tốn dòng chảy chất lượng nước hệ thống kênh, sơng (WFQ87) kỹ thuật chương trình”, Uỷ ban Quốc gia Chương trình Thuỷ văn Quốc tế Việt Nam Viện Khoa học Khí tượng Thủy văn Mơi trường, “Xây dựng chương trình dự báo xâm nhập mặn cho khu vực đồng sông Hồng – Thái Bình”, Báo cáo tổng kết đề tài Atlas Việt Nam http://baohatinh.vn/home/xa-hoi/bien-doi-khi-hau-o-ha-tinh-thuc-trang-va-giaiphap/1k54003.aspx Trần Ngọc Ánh, Nguyễn Thanh Sơn sử dụng mơ hình MIKE11 để đánh giá dự báo nhiễm mặn số sơng tỉnh Quảng Bình, Quảng Trị 2007 Liên giám thống kê tỉnh Hà Tĩnh 2012 Kịch Biến đổi Khí hậu, nước biển dâng cho Việt Nam” BTNMT 2012 Phạm Văn Vỵ, 2009, “Ảnh hưởng số yếu tố khí tượng tới q trình thủy lực vùng cửa sơng Hồng – Thái Bình”, Tạp chí Khoa học Đại học Quốc gia Hà Nội, Khoa học Tự nhiên Công nghệ 25, Số 1S (2009) 140‐144 10 Lê Sâm, 2004, “Nghiên cứu xâm nhập mặn phục vụ phát triển kinh tế - xã hội vùng ven biển đồng sông Cửu Long”, Đề tài NCKH cấp Nhà nước, KC08.18 thuộc Chương trình Bảo vệ Mơi trường Phịng tránh thiên tai, mã số KC08 11 Lã Thanh Hà, “Nghiên cứu khả dự báo xâm nhập vùng đồng sông Hồng – sơng Thái Bình mơ hình tốn”, Tạp chí KTTV tháng số 523 năm 2004 12 Trần Quốc Đạt “Mô xâm nhập mặn đồng sông Cửu Long tác động nước biển dâng suy giảm lưu lượng thượng nguồn”, Tạp chí Khoa học 2012:21b 141-150, Trường Đại học Cần Thơ 13 User guide, “A Modelling System for Rivers and Channels” 14 Reference, “Manual A Modelling System for Rivers and Channels” ... giới mặn 1‰ tiến sâu vào sông La khoảng 7,3km Sông La Sơng La Hình 22 Diễn biến XNM Sơng La theo kịch B2 (2060-2079 – 50cm) Ranh giới mặn 1‰ tiến sâu vào sông La khoảng 8,7km Sông La Sơng La Ngàn... thấy khoảng cách xâm nhập mặn khu vực sông Cả cụ thể Sông La vào tháng kiệt tương lai Đặc biệt kịch nước biển dâng 75cm 100cm, khoảng cách xâm nhập mặn 1‰ tiến vào tồn khu vực Sơng La tiến sâu vào... thiệu trên, nên việc lựa chọn MIKE1 1 để dự báo trình xâm nhập mặn nước sông vùng ven biển hợp lý 2.2 Cơ sở lý thuyết mơ hình Mike 11 2.2.1 Giới thiệu mơ hình MIKE 11 phần mềm kỹ thuật chuyên dụng