Hệ thống sông Hồng là một mạng lưới các con sông, tập hợp quanh con sông chính là sông Hồng, góp nước cho sông Hồng hoặc nhận nước của con sông này đổ ra Biển Đông. Hệ thống sông Hồng tạo nên phần lớn diện tích đồng bằng Bắc Bộ, một vùng bình nguyên tam giác châu thổ lớn thứ hai của Việt Nam. Cửa Ba Lạt là kết thúc của sông Hồng đổ ra biển nằm ở 20019’ độ vĩ Bắc, 106031’độ kinh Đông. Cửa Ba Lạt là nơi tiếp giáp địa giới hành chính giữa hai huyện Giao Thủy (Nam Định) và Tiền Hải (Thái Bình). Đây là khu vực ngập nước cửa sông mang ý nghĩa rất quan trọng về mặt kinh tế xã hội, sinh học cũng như nghiên cứu khoa học. 2.2. Khí hậu Khí hậu mang tính chất chung của khí hậu đồng bằng Bắc Bộ, thuộc khí hậu nhiệt đới gió mùa ẩm, có 4 mùa rõ rệt (xuân, hạ, thu, đông), mùa đông khí hậu khô do chịu tác động của gió mùa đông bắc. Khí hậu nhiệt đới gió mùa ẩm có mùa đông lạnh từ tháng 11 đến tháng 2. Mùa hè từ tháng 5 đến tháng 9. Nhiệt độ trung bình năm khoảng 240C. Lượng mưa trung bình 1175mm với số ngày mưa khoảng 133 ngày. Hai hướng gió chính trong năm ở đây là hướng Đông Bắc từ tháng 10 đến tháng 3 năm sau. Hướng Đông Nam từ tháng 4 đến tháng 9. 2.3. Thủy văn cửa sông Vùng ven biển tỉnh tỉnh Nam Định, Thái Bình có 4 cửa sông lớn, đó là cửa sông Ba Lạt (sông Hồng), cửa sông Ninh Cơ , cửa sông Đáy và cửa sông Trà Lý, sông Luộc. Mật độ sông trong vùng không cao (0,33kmkm2) nên khi lũ xảy ra vẫn có hiện tượng ngập úng tạm thời tại một số vùng, đặc biệt là đối với vùng ven biển nhu cầu rửa mặn rất lớn, do đó hệ thống sông này cần phải được tăng cường bằng các kênh mương tưới tiêu.
Báo cáo nghiên cứu khoa học NGHIÊN CỨU TÍNH TỐN LAN TRUYỀN MẶN KHU VỰC CỬA SÔNG BA LẠT SVTH: Đồng Thị Dung, lớp 54B2 GVHD: PGS.TS Nghiêm Tiến Lam SV: Đồng Thị Dung _ Lớp:54B2 Contents Đặt vấn đề 1.1 Mục tiêu phạm vi nghiên cứu 1.2 Phương pháp nghiên cứu Giới thiệu khu vực nghiên cứu 2.1 Vị trí địa lý 2.2 Khí hậu 2.3 Thủy văn cửa sông 2.4 Hải văn biển 2.5 Độ mặn 2.6 Hiện trạng xâm nhập mặn Ứng dụng mơ hình EFDC tính toán thủy động lực xâm nhập mặn vùng cửa sông 3.1 Tổng quan mơ hình tốn EFDC 3.2 Thiết lập mơ hình EFDC cho khu vực Cửa Ba Lạt 3.2.1 (1) Tài liệu địa hình (2) Số liệu thủy, hải văn 3.3 Thiết lập mơ hình thủy động lực Cửa Ba Lạt 3.3.1 Thiết lập miền tính tốn 3.3.2 Thiết lập điều kiện biên 3.3.3 Thiết lập điều kiện ban đầu 10 3.3.4 Hiệu chỉnh mơ hình thủy động lực 11 a Thời kỳ tính tốn 12 b Thiết lập thơng số mơ hình thủy lực 12 3.3.5 3.4 Các số liệu Kiểm định mơ hình thủy động lực 15 Thiết lập mơ hình tính tốn lan truyền mặn cho khu vực Cửa Ba Lạt 17 3.4.1 Phạm vi miền tính tốn 17 3.4.2 Thời gian tính tốn 17 3.4.3 Điều kiện ban đầu 17 3.4.4 Điều kiện biên 17 3.4.5 Hiệu chỉnh mơ hình lan truyền mặn 18 3.4.6 Kiểm định mơ hình lan truyền mặn 22 Kết luận kiến nghị 27 SV: Đồng Thị Dung _ Lớp:54B2 Đặt vấn đề Hiện nay, mực nước biển dâng cao khí hậu thay đổi nguy nghiêm trọng có tính tồn cầu, nguy trở nên nghiêm trọng nước có mật độ dân cư dày đặc vùng đất thấp ven biển Việt Nam Mực nước biển dâng cao làm trình xâm nhập mặn vùng cửa sơng thuộc dải ven biển đồng Bắc Bộ diễn biến phức tạp lấn sâu vào đất liền, ảnh hưởng đến trình lấy nước phục vụ ngành kinh tế Vì cần phải có tính tốn xác định diễn biến trình thủy lực xâm nhập mặn hệ thống vùng cửa sông ven biển Trong nghiên cứu ta tính tốn vận chuyển chất, lan truyền mặn cửa Ba Lạt, nơi Sông Hồng đổ biển 1.1 Mục tiêu phạm vi nghiên cứu Nghiên cứu thực nhằm đánh giá lan truyền, xâm nhập mặn khu vực cửa sơng Ba Lạt nằm hai tình Nam Định Thái Bình Phạm vi nghiên cứu tốn thủy lực lan truyền mặn khu vực cửa Ba Lạt với ảnh hưởng chế độ thủy triều dòng chảy từ sông Vùng nghiên cứu khu vực cửa sơng chịu tác động sóng, gió nên bỏ qua tác động q trình sóng, gió 1.2 Phương pháp nghiên cứu Nghiên cứu thực với phương pháp nghiên cứu phương pháp tiếp cận kế thừa, phương pháp phân tích tổng hợp phương pháp mơ hình tốn Nghiên cứu kế thừa kết nghiên cứu khu vực thực trước đây, kế thừa, áp dụng có chọn lọc kiến thức cơng cụ mơ hình thủy động lực học, lan tryền mặn có giới nước Đề tài nghiên cứu khai thác, sử dụng phần mềm mơ hình EFDC để mô thủy chế độ thủy động lực học, chất lượng nước khu vực nghiên cứu với công cụ xử lý số liệu khác Theo phương pháp mơ hình tốn, vào số liệu thu thập mục đích nghiên cứu, đồ án thực bước: v Thiết lập miền tính tốn, lưới tính, trạm kiểm tra SV: Đồng Thị Dung _ Lớp:54B2 v Thiết lập điều kiện ban đầu, điều kiện biên (biên trên, biên dưới) v Hiệu chỉnh, kiểm định mơ hình Sau hiệu chỉnh kiểm định, mơ hình dùng để tính tốn mơ phục vụ cho việc phân tích trạng chế độ thủy động lực, lan truyền mặn khu vực nghiên cứu Giới thiệu khu vực nghiên cứu 2.1 Vị trí địa lý Hệ thống sơng Hồng mạng lưới sông, tập hợp quanh sơng sơng Hồng, góp nước cho sông Hồng nhận nước sông đổ Biển Đông Hệ thống sông Hồng tạo nên phần lớn diện tích đồng Bắc Bộ, vùng bình nguyên tam giác châu thổ lớn thứ hai Việt Nam Cửa Ba Lạt kết thúc sông Hồng đổ biển nằm 20019’ độ vĩ Bắc, 106031’độ kinh Đông Cửa Ba Lạt nơi tiếp giáp địa giới hành hai huyện Giao Thủy (Nam Định) Tiền Hải (Thái Bình) Đây khu vực ngập nước cửa sông mang ý nghĩa quan trọng mặt kinh tế xã hội, sinh học nghiên cứu khoa học Hình 1: Khu vực nghiên cứu cửa sông Ba Lạt SV: Đồng Thị Dung _ Lớp:54B2 2.2 Khí hậu Khí hậu mang tính chất chung khí hậu đồng Bắc Bộ, thuộc khí hậu nhiệt đới gió mùa ẩm, có mùa rõ rệt (xuân, hạ, thu, đơng), mùa đơng khí hậu khơ chịu tác động gió mùa đơng bắc Khí hậu nhiệt đới gió mùa ẩm có mùa đơng lạnh từ tháng 11 đến tháng Mùa hè từ tháng đến tháng Nhiệt độ trung bình năm khoảng 240C Lượng mưa trung bình 1175mm với số ngày mưa khoảng 133 ngày Hai hướng gió năm hướng Đông Bắc từ tháng 10 đến tháng năm sau Hướng Đông Nam từ tháng đến tháng 2.3 Thủy văn cửa sông Vùng ven biển tỉnh tỉnh Nam Định, Thái Bình có cửa sơng lớn, cửa sông Ba Lạt (sông Hồng), cửa sông Ninh Cơ , cửa sông Đáy cửa sông Trà Lý, sông Luộc Mật độ sông vùng không cao (0,33km/km2) nên lũ xảy có tượng ngập úng tạm thời số vùng, đặc biệt vùng ven biển nhu cầu rửa mặn lớn, hệ thống sơng cần phải tăng cường kênh mương tưới tiêu Bảng 1: Một số trạm đo thủy văn STT Trạm Ba Lạt Dương Liễu Vũ Thuận Hành Thiện Ngô Xá Phú Hào Nhật Tảo Tiến Đức Triều Dương Kinh độ ’ 20 19 20020’ 20023’ 20020’ 20023’ 20025’ 20035’ 20035’ 20038’ SV: Đồng Thị Dung _ Lớp:54B2 Vĩ độ Ghi ’ 106 31 106024’ 106018’ 106018’ 106016’ 106013’ 106007’ 106024’ 106005’ Thủy Văn Thủy Văn Thủy Văn Thủy Văn Thủy Văn Thủy Văn Thủy Văn Thủy Văn Thủy Văn 2.4 Hải văn biển Chế độ sóng khu vực thay đổi theo mùa Vào mùa lạnh, hướng sóng ngồi khơi Đơng Bắc, Đơng, ven bờ hướng Đông, Đông Bắc Đông Nam Thủy triều vùng biển Nam Định-Thái Bình thuộc chế độ nhật triều, ngày có đỉnh triều chân triều, với biên độ triều thuộc loại lớn nước ta, biên độ triều trung bình từ 1,6 – 1,7m, lớn đạt 3,3m, nhỏ 0,1m Ảnh hưởng thủy triều đến sông vùng lớn 2.5 Độ mặn Ở ngồi khơi Biển Đơng không đổi, mùa mưa độ mặn khoảng 32‰ mùa khơ 33‰ Ở vùng ven biển, độ nhiễm mặn thay đổi theo mùa ảnh hưởng nước từ sông đổ vào Chiều dài xâm nhập mặn trung bình 1‰ 4‰ sơng Hồng tương ứng 12 km 10 km Chiều dài xâm nhập mặn 1‰ sâu xảy sông Hồng 14 km Độ mặn sông ven biển tăng từ đầu mùa đến mùa khơ sau giảm dần đến cuối mùa vụ Sự thay đổi có liên quan tới dòng nước từ thượng nguồn đổ 2.6 Hiện trạng xâm nhập mặn Cùng với tốc độ phát triển theo hướng cơng nghiệp hóa, đại hóa nước, tỉnh lân cận tỉnh Thái Bình, Nam Định bước chuyển dần từ tỉnh canh tác nơng nghiệp sang xây dựng hình thành khu, cụm cơng nghiệp Nguồn nước cửa Ba Lạt có biến động tác động khí tượng thủy văn hoạt động người Bên cạnh hai bên bờ nhu cầu sử dụng nước ngày tăng tăng dân số phát triển kinh tế - xã hội Vì cơng tác điều tra, dự báo, xác định vùng xâm nhập mặn, mức độ lan truyền mặn cần thiết công phát triển kinh tể xã hội Ứng dụng mơ hình EFDC tính tốn thủy động lực xâm nhập mặn vùng cửa sông 3.1 Tổng quan mô hình tốn EFDC Mơ hình EFDC (Environmental Fluid Dynamics Code) phần mềm mơ hình tốn có khả dự báo, tính tốn mơ q trình dòng chảy, lan truyền vật chất có tính đến q trình sinh - địa - hóa sơng, suối, hồ, cửa sông, ven biển, biển đại dương Mơ hình quan bảo vệ mơi trường Mỹ US EPA SV: Đồng Thị Dung _ Lớp:54B2 phát triển từ năm 1980, đến 1994 nhà khoa học Viện Khoa Học Biển Virgina tiếp tục xây dựng Mơ hình xây dựng dựa phương trình động lực học (bảo tồn vật chất, bảo tồn động lượng bảo toàn lượng), nguyên tắc bảo tồn khối lượng bảo tồn thể tích Mơ hình mơ hình đa chiều (1 chiều, chiều, chiều) nên có khả đạt độ xác cao việc mơ hình hóa hệ thống đầm lầy, đất ngập nước, kiểm sốt dòng chảy, dòng sinh sóng gần bờ q trình vận chuyển trầm tích Mơ hình EFDC gồm modul chính: mơ hình thủy động lực học, mơ hình chất lượng nước, mơ hình vận chuyển trầm tích, mơ hình lan truyền – phân hủy chất độc môi trường nước mặt Mơ hình thủy động lực học EFDC lại gồm modul lan truyền vận chuyển, bao gồm: động lực học, màu sắc, nhiệt độ, độ mặn… Mơ hình EFDC gồm modul chính: Mơ hình thủy động lực học Mơ hình chất lượng nước Mơ hình vận chuyển trầm tích Mơ hình lan truyền, phân hủy chất độc môi trường nước mặt Mô-đun thủy động lực mơ hình EFDC đựa phương trình xấp xỉ thủy tĩnh chiều cho hệ tọa độ theo phương thẳng đứng tọa độ cong trực giao nằm ngang 3.2 Thiết lập mơ hình EFDC cho khu vực Cửa Ba Lạt 3.2.1 Các số liệu Với mục tiêu nghiên cứu tính tốn vận chuyển chất sông Ba Lạt Dựa vào trạng hệ thống trạm quan trắc, tình hình số liệu thu thập khu vực tính tốn, lựa chọn sử dụng số liệu địa hình, thủy hải văn (1) Tài liệu địa hình Tài liệu địa hình sơng Hồng thu thập từ hai nguồn: SV: Đồng Thị Dung _ Lớp:54B2 - Địa hình sông lấy từ số liệu điều tra mặt cắt sông Hồng hàng năm Bộ nông nghiệp năm 2009 - Địa hình vùng cửa sơng lấy dựa vào Đề án 47 “ Về việc phê duyệt đề án tổng thể điều tra quản lý tài nguyên – môi trường biển đến năm 2010, tầm nhìn đến năm 2020” (2) Số liệu thủy, hải văn Số liệu thủy văn sử dụng bao gồm số liệu mực nước đo đạc theo năm 2002, 2003 trạm thủy văn Ba Lạt Dương Liễu 3.3 Thiết lập mơ hình thủy động lực Cửa Ba Lạt 3.3.1 Thiết lập miền tính tốn Miền tính tốn mơ hình bao gồm phần hạ du lưu vực sơng Hồng – Thái Bình, bao gồm nhánh sơng Đào, Ninh Cơ, Trà Lý, sông Luộc Biên khu vực gần thành phố Hưng Yên, biên lấy biên thủy triều ngồi biển khơi xa khơng chịu ảnh hưởng dòng chảy sơng Do khơng thu thập số liệu địa hình nhánh sơng nhỏ mạng lưới kênh rạch lưu vực ảnh hưởng tới kết tính tốn khơng nhiều nên khuôn khổ đồ án sông nhỏ khơng mơ Miền mơ hình xây dựng thuộc dạng lưới cong trực giao Đây dạng lưới mơ hình phù hợp với vùng nghiên cứu đáp ứng đặc điểm địa hình dòng chảy sơng có độ xác cao so với dòng chảy thực tế Trong đồ án sử dụng phần mềm Delft3D tạo lưới cho lưu vực nghiên cứu để xây dựng miền lưới tính tốn mơ cho vùng nghiên cứu Trên miền tính có tất 3859 lưới, đưa vào mơ hình EFDC theo lựa chọn Inport Grid Vùng lưới mơ hình EFDC thể Hình SV: Đồng Thị Dung _ Lớp:54B2 Hình 2: Lưới tính tốn mơ hình EFDC 3.3.2 Thiết lập điều kiện biên Điều kiện biên sử dụng giá trị mực nước thực đo trạm quan trắc lưu vực sông Hồng năm 2002 SV: Đồng Thị Dung _ Lớp:54B2 Hình 3: Vị trí biên tính tốn 3.3.3 Thiết lập điều kiện ban đầu Khi thiết lập điều kiện ban đầu cần khai báo liệu: - File mực nước ban đầu (Surface Elevations): dựa vào mực nước thực đo trạm phía thượng lưu hạ lưu lấy độ dốc mặt nước số ta có cao trình mặt nước theo đường mặt cắt dọc sơng Các điểm lưới lại tồn miền mơ hình EFDC có khả tự nội suy số liệu mực nước toàn vùng nghiên cứu dạng (X, Y, Z) - File cao trình đáy (Bottom Elevations) Kết số hóa đồ nội suy địa hình ta đươc số liệu địa hình mơ hình EFDC Hình 19 Cốt cao độ địa hình lấy theo cao độ chuẩn quốc gia Địa hình toàn lưu vực nghiên cứu dao động từ cao độ -28,035 m +2,576 m SV: Đồng Thị Dung _ Lớp:54B2 10 Hình 7: Phân bố độ sâu mực nước sông Hồng pha triều rút ngày 16/3/2002 Hình 8: Phân bố độ sâu mực nước sơng Hồng pha triều dâng ngày 25/3/2002 Hình 9: Phân bố vận tốc khu vực cửa Ba Lạt lúc 03:00 ngày 25/3/2002 SV: Đồng Thị Dung _ Lớp:54B2 14 3.3.5 Kiểm định mơ hình thủy động lực Sau hiệu chỉnh, thông số thủy lực chọn sử dụng để kiểm định mơ hình Mơ hình kiểm định trạm Ba Lạt, thời đoạn kiểm định mơ hình tháng năm 2003 Bước thời gian thực mô kiểm định lựa chọn 5s, độ rối 10m2/s Hình 10: Đường q trình mực nước thực đo tính tốn tháng năm 2003 cửa Ba Lạt thời kỳ kiểm định Hình 11: Phân bố mực nước khu vực cửa Ba Lạt triều dâng ngày 28/03/2003 SV: Đồng Thị Dung _ Lớp:54B2 15 Hình 12: Phân bố vận tốc khu vực cửa Ba Lạt lúc 03:00 ngày 28/3/2003 Nhận xét kết kiểm định mơ hình thủy lực: Kết kiểm định đường trình mực nước thực đo tính tốn trạm Ba Lạt tốt Dao động mực nước trạm tính từ mơ hình với số liệu thực đo có trùng khớp pha triều độ lớn (Hình 29) Kết tính tốn tiêu Nash RMSE Ba Lạt đạt yêu cầu cho phép với giá trị trình bày Bảng Bảng 3: Kết đánh giá sai số kiểm định mô hình Trạm Nash (%) RMSE Ba Lạt 91 0.198 Kết tương đối tốt hiệu chỉnh kiểm định mơ hình thơng qua tiêu đánh giá Nash RMSE Vì vậy, kết luận thơng số thủy lực ta chọn có đủ độ tin cậy để xây dựng mơ hình tính tốn lan truyền mặn SV: Đồng Thị Dung _ Lớp:54B2 16 3.4 Thiết lập mơ hình tính tốn lan truyền mặn cho khu vực Cửa Ba Lạt Mơ hình xây dựng nhằm mục đích đánh giá so sánh biến đổi trình lan truyền mặn thay đổi số tính chất mơ hình So sánh thay đổi chạy số liệu thực tế mơ hình từ ta đánh giá mức độ tin cậy, mơ hình từ áp dụng mơ hình vào để dự tính trường hợp xảy tương lai Các số liệu miền tính tốn, điều kiện biên mực nước, điều kiện ban đầu mực nước ta lấy mơ hình thủy lực 3.4.1 Phạm vi miền tính tốn Mơ hình mơ lan truyền mặn dùng kết mơ hình thủy động lực làm điều kiện cho tính tốn Mơ hình tính tốn cho thủy lực có kết tương đối tốt, nên tất thông số thủy lực chuyển sang sử dụng tính tốn mơ hình lan truyền mặn 3.4.2 Thời gian tính tốn Mơ hình tính toán lan truyền mặn thiết lập chạy với thời gian tháng, bước thời gian chạy mô hình 5s 3.4.3 Điều kiện ban đầu Dữ liệu địa hình lấy từ mơ hình thủy lực Hình 13: Vị trí trạm đo độ mặn 3.4.4 Điều kiện biên SV: Đồng Thị Dung _ Lớp:54B2 17 Biên mực nước: giữ ngun thơng số mơ hình thủy lực Biên mặn: Đối với liệu biên mặn ta thiết lập vào cửa sổ Hình 31 cách nhập trực tiếp số liệu: - Ta vào mục Salinity => time series Data chọn Edit sau nhập số liệu hình: Hình 14: Cửa sổ thiết lập biên mặn mơ hình 3.4.5 Hiệu chỉnh mơ hình lan truyền mặn Kết mơ hình chạy độ mặn cho tháng năm 2002 trạm Ba Lạt Hình 15: Đường q trình độ mặn thực đo tính toán tháng năm 2002 tầng mặt cửa Ba Lạt thời kỳ hiệu chỉnh SV: Đồng Thị Dung _ Lớp:54B2 18 Hình 16: Đường trình độ mặn thực đo tính tốn tháng năm 2002 tầng cửa Ba Lạt thời kỳ hiệu chỉnh Hình 17: Đường q trình độ mặn thực đo tính tốn tháng năm 2002 tầng đáy cửa Ba Lạt thời kỳ hiệu chỉnh Từ kết Hình 15, Hình 16 Hình 17 thấy vào tháng (mùa nước cạn) diễn biến độ mặn định dòng triều Khi ta xét đến phân bố nồng độ mặn kỳ triều lên rút tháng, dịch chuyển độ mặn chiều với lên xuống dòng triều Trong ba lớp ta thấy độ mặn lớp mặt thấp lớp đáy, lớp đáy có độ mặn cao Nguyên nhân dẫn tới chênh lệch tỷ trọng nước biển nằm khoảng 1.020 tới 1.030 kg/m³ bề mặt sâu lòng đại dương, áp suất cao, nước biển đạt tỷ trọng riêng tới 1.050 kg/m³ hay cao Như nước biển nặng nước (nước tinh khiết đạt tỷ trọng riêng tối đa 1.000 g/ml nhiệt độ °C) trọng lượng bổ sung muối tượng điện giải Trạm Ba Lạt gần cửa Ba Lạt (nơi sông Hồng đổ biển) nên thay đổi thủy lực dòng chảy sơng chịu ảnh hưởng lớn dòng chiều SV: Đồng Thị Dung _ Lớp:54B2 19 Kết mơ hình chạy độ mặn cho tháng năm 2002 trạm Dương Liễu Hình 18: Đường q trình độ mặn thực đo tính tốn tháng năm 2002 tầng mặt trạm Dương Liễu thời kỳ hiệu chỉnh Hình 19: Đường trình độ mặn thực đo tính tốn tháng năm 2002 tầng trạm Dương Liễu thời kỳ hiệu chỉnh Hình 20: Đường trình độ mặn thực đo tính tốn tháng năm 2002 tầng đáy trạm Dương Liễu thời kỳ hiệu chỉnh Từ kết Hình 18, Hình 19 Hình 20 thấy Trạm Dương Liễu không bị xâm nhập mặn nhiều (độ mặn vào khoảng 0.050/00), nguyên nhân Trạm Dương Liễu nằm sâu bên sông, trạm cách cửa sông Ba Lạt 27km phía thượng nguồn SV: Đồng Thị Dung _ Lớp:54B2 20 nên Dương Liễu khơng chịu nhiều tác động thủy triều mà chịu tác động dòng chảy thủy lực sơng Hình 21: Phân bố độ mặn triều lên (25/3/2002) Hình 22: Trường dòng chảy triều lên lúc 02:00 (9/3/2002) Hình 23: Phân bố độ mặn triều xuống (20/3/2002) SV: Đồng Thị Dung _ Lớp:54B2 21 Hình 24: Diễn biến xâm nhập mặn dọc sông pha triều dâng ngày 25/3/2002 Hình 25: Diễn biến xâm nhập mặn dọc sông pha triều rút ngày 16/3/2002 Từ hai hình Hình 24 Hình 25 ta nhận thấy yếu tố thủy lực, vận tốc, lưu lượng vùng cửa sông chịu tác động ảnh hưởng nhiều thủy triều Cụ thể triều dâng, vận tốc sông nhỏ, nước từ biển vào sông làm mực nước tăng, xâm nhập nước mặn vào tầng đáy sông Khi triều cao, lưu lượng vào sông nhiều dẫn tới biên mặn sâu lên thượng nguồn gây tổn hại đến kinh tế, ôi nhiễm mặn nguồn nước ngọt, ảnh hưởng tới đời sống, sinh hoạt người dân vùng 3.4.6 Kiểm định mơ hình lan truyền mặn Sau hiệu chỉnh, thơng số tính tốn lan truyền mặn chọn sử dụng để kiểm định mô hình Mơ hình kiểm định trạm Ba Lạt Dương Liễu thời đoạn kiểm định mơ hình tháng năm 2003 Bước thời gian thực mô kiểm định lựa chọn 5s, độ rối 10m2/s Kết mô kiểm định diễn biến mặn trạm thủy văn Ba Lạt, Dương Liễu thể hình SV: Đồng Thị Dung _ Lớp:54B2 22 •Trạm Ba Lạt Hình 26: Đường q trình độ mặn thực đo tính tốn tháng năm 2003 tầng mặt trạm Ba Lạt thời kỳ kiểm định Hình 27: Đường trình độ mặn thực đo tính tốn tháng năm 2003 tầng trạm Ba Lạt thời kỳ kiểm định Hình 28: Đường trình độ mặn thực đo tính tốn tháng năm 2003 tầng đáy trạm Ba Lạt thời kỳ kiểm định • Trạm Dương Liễu SV: Đồng Thị Dung _ Lớp:54B2 23 Hình 29: Đường q trình độ mặn thực đo tính tốn tháng năm 2003 tầng mặt trạm Dương Liễu thời kỳ kiểm định Hình 30: Đường trình độ mặn thực đo tính tốn tháng năm 2003 tầng trạm Dương Liễu thời kỳ kiểm định Hình 31: Đường q trình độ mặn thực đo tính toán tháng năm 2003 tầng đáy trạm Dương Liễu thời kỳ kiểm định SV: Đồng Thị Dung _ Lớp:54B2 24 Hình 32: Diễn biến độ mặn triều dâng ngày 28/3/2003 Hình 33: Diễn biến độ mặn triều rút ngày 10/3/2003 Nhận xét: Kết thực đo tính tốn tương đối tốt số liệu thực đo đo vài thời điểm nên biết thị hình có vài điểm, số liệu chạy mơ hình số liệu đo theo nên biểu đồ so sánh mang tính chất tương đối Tuy nhiên phần thể kết chạy mơ hình so với thực tế tương đồng nhau, thể qua việc biểu thị pha triều dâng pha triều rút tang giảm khoảng thời gian Cụ thể pha triều dâng vào ngày 14,15, 16,17,18,27,28,29,30 Trạm Dương Liễu cách xa cửa sơng phía thượng nguồn nên bị tác động, ảnh hưởng thủy triều nên có độ mặn thấp SV: Đồng Thị Dung _ Lớp:54B2 25 Hình 34: Diễn biến xâm nhập mặn dọc sông pha triều rút ngày 10/3/2003 Hình 35: Diễn biến xâm nhập mặn dọc sông pha triều dâng ngày 28/3/2003 SV: Đồng Thị Dung _ Lớp:54B2 26 Kết luận kiến nghị Báo cáo trình bày khái quát đặc điểm tự nhiên trạng xâm nhập mặn khu vực đồng sơng Hồng nói chung khu vực Cửa Ba Lạt nói riêng thiết lập mơ hình tốn thủy động lực lan truyền mặn cho Cửa Ba Lạt Mơ hình tốn thủy động lực lan truyền mặn cho Cửa Ba Lạt thiết lập dựa phần mềm mơ hình EFDC hiệu chỉnh kiểm định với số liệu thực đo khu vực Một số kết thu nghiên cứu bao gồm: • Từ điều kiện tự nhiên số liệu địa hình, thủy triều, số liệu độ mặn thu thập nghiên cứu xác định thơng số hợp lý cho mơ hình thủy lực chiều EFDC cho khu vực cửa Ba Lạt • Đã đánh giá ảnh hưởng thủy triều tới xâm nhập mặn vào sông quan trọng • Đã thiết lập mơ hình mơ q trình xâm nhập mặn cho khu vực cửa Ba Lạt Tuy nhiên, kết nghiên cứu đề tài dừng lại nghiên cứu toán chiều biến đổi xâm nhập mặn vùng cửa sông Ba Lạt Kiến nghị quan quản lý cần nghiên cứu đầy đủ toán lan truyền mặn chiều thu thập bổ sung thêm số liệu làm tăng độ tin cậy độ xác mơ hình để làm sở cho toán đánh giá xâm nhập mặn phục vụ cho sản xuất đề xuất giải pháp nhằm khắc phục, hạn chế rủ ro tượng xâm nhập mặn gây Tài liệu tham khảo Số liệu địa hình: Số liệu điều tra mặt cắt sơng Hồng hàng năm Bộ nông nghiệp năm 2009 Đề án 47 Số liệu thủy văn: Mạng lưới trạm đo đạc, số liệu mực nước trạm (thu thập tài liệu đo theo giờ) (trang waterdata.vn www.hymetdata.gov.vn) Số liệu hải văn: Mực nước thực đo theo ven biển), số liệu độ mặn Số liệu dân sinh, kinh tế, xã hội (thu thập) Luận văn thạc sĩ Nguyễn Thanh Tâm “Nghiên cứu vân chuyển bùn cát khu vực lạch huyện” Phần mềm EFDC, hướng dẫn sử dụng, sở lý thuyết phương pháp số http://www.ds-intl.biz Trang wed “Viện nước, tưới tiêu môi trường”: iwe.vn SV: Đồng Thị Dung _ Lớp:54B2 27 SV: Đồng Thị Dung _ Lớp:54B2 28 ... vi nghiên cứu Nghiên cứu thực nhằm đánh giá lan truyền, xâm nhập mặn khu vực cửa sơng Ba Lạt nằm hai tình Nam Định Thái Bình Phạm vi nghiên cứu toán thủy lực lan truyền mặn khu vực cửa Ba Lạt. .. đồng sơng Hồng nói chung khu vực Cửa Ba Lạt nói riêng thiết lập mơ hình tốn thủy động lực lan truyền mặn cho Cửa Ba Lạt Mơ hình tốn thủy động lực lan truyền mặn cho Cửa Ba Lạt thiết lập dựa phần... EFDC cho khu vực cửa Ba Lạt • Đã đánh giá ảnh hưởng thủy triều tới xâm nhập mặn vào sông quan trọng • Đã thiết lập mơ hình mơ trình xâm nhập mặn cho khu vực cửa Ba Lạt Tuy nhiên, kết nghiên cứu đề