Mô phỏng trường thủy động lực khu vực cửa đà diễn, tỉnh phú yên

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - Bùi Thị Hạnh MÔ PHỎNG TRƢỜNG THỦY ĐỘNG LỰC KHU VỰC CỬA ĐÀ DIỄN, TỈNH PHÚ YÊN LUẬN VĂN THẠC SĨ KHOA HỌC Hà Nội - 2016 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - Bùi Thị Hạnh MÔ PHỎNG TRƢỜNG THỦY ĐỘNG LỰC KHU VỰC CỬA ĐÀ DIỄN, TỈNH PHÚ YÊN Chuyên ngành: Thủy văn học Mã số: 60 44 02 24 LUẬN VĂN THẠC SĨ KHOA HỌC NGƢỜI HƢỚNG DẪN KHOA HỌC: PGS TS Trần Ngọc Anh TS Đặng Thanh Mai Hà Nội - 2016 MỤC LỤC DANH MỤC HÌNH DANH MỤC BẢNG LỜI CẢM ƠN MỞ ĐẦU 10 Chƣơng - TỔNG QUAN 11 1.1 Tổng quan điều kiện địa lý tự nhiên & KTXH lưu vực sông Ba cửa sông Đà Diễn 11 1.1.1 Đặc điểm tự nhiên lưu vực sông Ba cửa Đà Diễn 11 1.1.2 Tổng quan điều kiện kinh tế xã hội lưu vực sông Ba 26 1.2 Tổng quan mô hình lựa chọn mơ hình 27 1.2.1 Tổng quan 27 1.2.2 Lựa chọn mơ hình 28 1.2.3 Giới thiệu module tốn mơ hình MIKE21/3 FM COUPLE 29 CHƢƠNG - THIẾT LẬP MƠ HÌNH 30 2.1 Cơ sở lý thuyết mơ hình MIKE 21 30 2.2 Cơ sở liệu 34 2.2.1 Số liệu địa hình 34 2.2.2 Sóng, gió 34 2.2.3 Số liệu thủy, hải văn 35 2.3 Thiết lập mơ hình 36 2.3.1 Miền tính lưới tính 36 2.3.2 Điều kiện biên 38 2.4 Hiệu chỉnh mơ hình 38 2.5 Kiểm định mơ hình 40 CHƢƠNG - ỨNG DỤNG MƠ HÌNH MƠ PHỎNG TRƢỜNG THỦY ĐỘNG LỰC CỬA SÔNG ĐÀ DIỄN 43 3.1 Xây dựng kịch 43 3.2 Kết tính tốn trạng 46 3.2.1 Kết tháng 11/2015 46 3.2.2 Kết tháng 5/2016 49 3.3 Kết tính tốn với kịch 52 3.3.1 Kịch sóng Đơng Bắc (KB_NE) 52 3.3.2 Kịch sóng Đông (KB_E) 56 3.3.3 Sóng Nam (S) 61 KẾT LUẬN VÀ KIẾN NGHỊ 65 Kết luận 65 Kiến nghị 67 TÀI LIỆU THAM KHẢO 68 DANH MỤC HÌNH Hình 1.1 Bản đồ lƣu vực sơng Ba 11 Hình 1.2 Bản đồ khu vực tỉnh Phú Yên (tỉ lệ 1:600.000) [2] 12 Hình 1.3 Ảnh chụp cửa Đà Diễn từ vệ tinh (nguồn Google Earth) 12 Hình 1.4 Bản đồ phân bố lƣợng mƣa mùa mƣa (tỉ lệ 1:100.000) 16 Hình 1.5 Bản đồ phân bố lƣợng mƣa mùa khơ (tỉ lệ 1:100.000) 16 Hình 1.6 Sơ đồ vùng hạ lƣu sông Ba từ Củng Sơn đến cửa Đà Diễn .22 Hình 1.7 Hoa sóng tính từ số liệu sóng tồn cầu 25 Hình 1.8 Các module mơ hình MIKE21/3 FM COUPLED 29 Hình 2.1 Sơ đồ lƣới tính so le theo hai chiều x y .31 Hình 2.2 Sơ đồ quét thời gian trung tâm 32 Hình 2.3 Địa hình khu vực nghiên cứu tháng 11/2015 34 Hình 2.4 Địa hình khu vực nghiên cứu tháng 5/2016 .34 Hình 2.5 Sơ đồ trạm đo đợt khảo sát tháng 11/2015 .35 Hình 2.6 Sơ đồ trạm đo đợt khảo sát tháng 5/2016 35 Hình 2.7 Miền tính khu vực nghiên cứu 36 Hình 2.8 Lƣới Mesh II cho miền tính lớn 37 Hình 2.9 Lƣới tính Mesh I cho khu vực nghiên cứu 37 Hình 2.10 Mực nƣớc trạm C mực nƣớc tính tốn đƣợc từ mơ hình 39 Hình 2.11 So sánh dịng chảy tính toán với thực đo trạm C 39 Hình 2.12 So sánh độ cao sóng trung bình tính tốn với thực đo trạm C 39 Hình 2.13 So sánh mực nƣớc tính tốn với thực đo trạm D 40 Hình 2.14 Dịng chảy tính tốn với thực đo trạm E .40 Hình 2.15 Độ cao sóng tính tốn với thực đo trạm E 41 Hình 2.16 Mực nƣớc triều tính tốn với thực đo trạm F 41 Hình 2.17 Độ cao sóng tính tốn với thực đo trạm F 41 Hình 2.18 Dịng chảy tính toán với thực đo trạm G 42 Hình 2.19 Độ cao sóng tính tốn với thực đo trạm G 42 Hình 3.1 Chế độ sóng tháng lũ kiệt khu vực cửa sông Đà Diễn 43 Hình 3.2 Đƣờng tần suất chiều cao sóng Đơng Bắc lớn 44 Hình 3.3 Đƣờng tần suất chiều cao sóng Đơng lớn 45 Hình 3.4 Đƣờng tần suất chiều cao sóng Nam lớn 45 Hình 3.5 Kết trƣờng sóng tháng 11/2015 46 Hình 3.6 Kết trƣờng sóng tháng 11/2015 thời kì triều cao 47 Hình 3.7 Kết trƣờng sóng tháng 11/2015 thời kì triều thấp .47 Hình 3.8 Trƣờng dịng chảy cửa Đà Diễn thời điểm đỉnh triều 11/2015 47 Hình 3.9 Trƣờng dịng chảy cửa Đà Diễn thời điểm chân triều 11/2015 .48 Hình 3.10 Trƣờng dịng chảy phía cửa Đà Diễn tháng 11/2015 triều rút 48 Hình 3.11 Trƣờng dịng chảy phía cửa Đà Diễn tháng 11/2015 triều lên .48 Hình 3.12 Kết trƣờng sóng tháng 5/2016 49 Hình 3.13 Trƣờng sóng cửa Đà Diễn thời điểm đỉnh triều 5/2016 50 Hình 3.14 Trƣờng sóng cửa Đà Diễn thời điểm chân triều 5/2016 50 Hình 3.15 Trƣờng dòng chảy cửa Đà Diễn thời điểm đỉnh triều 5/2016 51 Hình 3.16 Trƣờng dịng chảy cửa Đà Diễn thời điểm chân triều 5/2016 .51 Hình 3.17 Trƣờng sóng tồn miền tính tốn KB_NE1% 52 Hình 3.18 Trƣờng sóng tồn miền tính tốn KB_NETB 52 Hình 3.19 Trƣờng sóng khu vực cửa sơng KB_NE1% (đỉnh triều) 53 Hình 3.20 Trƣờng sóng khu vực cửa sông KB_NE1% (chân triều) 53 Hình 3.21 Trƣờng dịng chảy vùng cửa sơng KB_NE1% (đỉnh triều) .53 Hình 3.22 Trƣờng dịng chảy vùng cửa sơng KB_NETB (đỉnh triều) .54 Hình 3.23 Trƣờng sóng vùng cửa sơng KB_NETB thời kì đỉnh triều 54 Hình 3.24 Trƣờng sóng vùng cửa sơng KB_NETB thời kì chân triều .54 Hình 3.25 Vùng xốy trƣớc cửa sơng KB_NE1% (triều lên) 55 Hình 3.26 Vùng xốy trƣớc cửa sông KB_NETB (triều rút) 55 Hình 3.27 Trƣờng sóng miền tính tốn KB_E1% 57 Hình 3.28 Trƣờng sóng miền tính tốn KB_ETB .57 Hình 3.29 Trƣờng sóng khu vực cửa sơng KB_E1% (đỉnh triều) 57 Hình 3.30 Trƣờng sóng khu vực cửa sông KB_E1% (chân triều) 58 Hình 3.31 Trƣờng sóng khu vực cửa sơng KB_ETB (chân triều) 58 Hình 3.32 Trƣờng sóng khu vực cửa sông KB_ETB (đỉnh triều) 58 Hình 3.33 Trƣờng dịng chảy khu vực cửa sơng KB_E1% (đỉnh triều) 59 Hình 3.34 Trƣờng dịng chảy khu vực cửa sông KB_E1% (chân triều) 59 Hình 3.35 Trƣờng dịng chảy khu vực cửa sơng KB_ETB (đỉnh triều) 60 Hình 3.36 Trƣờng dịng chảy khu vực cửa sông KB_ETB (chân triều) 60 Hình 3.37 Trƣờng sóng KB_S1% .61 Hình 3.38 Trƣờng sóng KB_S1% (chân triều) 62 Hình 3.39 Trƣờng sóng KB_S1% (đỉnh triều) 62 Hình 3.40 Trƣờng dịng chảy KB_S1% (đỉnh triều) 62 Hình 3.41 Trƣờng dịng chảy KB_S1% (chân triều) 63 Hình 3.42 Trƣờng sóng KB_STB (đỉnh triều) 63 Hình 3.43 Trƣờng sóng KB_STB (chân triều) 63 Hình 3.44 Trƣờng dịng chảy KB_STB (đỉnh triều) 64 Hình 3.45 Trƣờng dịng chảy KB_STB (chân triều) 64 DANH MỤC BẢNG Bảng 1.1 Tần suất (%) hƣớng gió thịnh hành khu vực Phú Yên .15 Bảng 1.2 Lƣợng mƣa trung bình nhiều năm tháng (mm) 17 Bảng 1.3 Nhiệt độ trung bình tháng năm (oC) 18 Bảng 1.4 Cao độ biên độ thủy triều 24 Bảng 1.5 Tần suất hƣớng sóng thịnh hành theo số liệu sóng tồn cầu 25 Bảng 3.1 Kịch sóng 44 Bảng 3.2 Các tham số thống kê đƣờng tần suất chiều cao sóng 45 Bảng 3.3 Các đặc trƣng sóng 46 LỜI CẢM ƠN Nghiên cứu đƣợc thực khuôn khổ đề tài luận văn cao học chuyên ngành Thủy văn học, khóa 2013 - 2015 Khoa Khí tƣợng, Thủy văn Hải dƣơng học, Trƣờng Đại học Khoa học Tự nhiên Luận văn đƣợc hoàn thành dƣới hƣớng dẫn tận tình, PGS TS Trần Ngọc Anh TS Đặng Thanh Mai Tác giả xin bày tỏ lời cảm ơn chân thành tới tập thể ngƣời hƣớng dẫn khoa học bảo, hỗ trợ, động viên suốt trình học tập, nghiên cứu hoàn thành luận văn Tác giả xin gửi lời cảm ơn tới ThS Đặng Đình Đức cán Trung tâm Động lực học thủy khí mơi trƣờng ln dành thời gian giải đáp, thảo luận vấn đề khúc mắc liên quan đến ứng dụng mơ hình q trình thực luận văn Tác giả xin bày tỏ biết ơn sâu sắc đến PGS TS Nguyễn Tiền Giang, Ban Chủ nhiệm thành viên cho phép hỗ trợ tác giả sử dụng nguồn số liệu quý giá khu vực cửa Đà Diễn, tỉnh Phú Yên, mà thiếu chắn luận văn khơng thể đƣợc hồn thành Cuối cùng, tác giả xin chân thành cảm ơn thầy cô Bộ môn Thủy văn Tài ngun nƣớc, Văn phịng Khoa Khí tƣợng Thủy văn Hải dƣơng học quan tâm động viên tạo điều kiện thuận lợi để tác giả hồn thành nhiệm vụ Do thời gian kiến thức hạn chế, luận văn không tránh khỏi thiếu sót Tác giả mong nhận đƣợc ý kiến đóng góp q thầy cơ, bạn bè, đồng nghiệp, nhà khoa học độc giả để luận văn đƣợc hoàn thiện Xin trân trọng cảm ơn! Hà Nội, ngày tháng Học viên năm 2016 MỞ ĐẦU Cửa sông Đà Diễn cửa sông hệ thống sơng Ba, cửa sơng quan trọng thành phố Tuy Hịa, tỉnh Phú Yên Cửa sông vừa cảng cá vừa nơi neo đậu 900 tàu thuyền, nhƣng địa hình lòng dẫn khu vực thƣờng xuyên bị biến đổi sau mùa bão lũ gây ảnh hƣởng không nhỏ đến giao thơng thủy, lũ phát triển kinh tế - xã hội Hiện nay, tỉnh Phú Yên xây kè chắn sóng, tổ chức nạo vét cửa sơng liên tục nhƣng thời gian sau tình trạng bồi lấp lại tái diễn Hiện tƣợng xói lở xảy nghiêm trọng bờ phía Bắc nƣớc thƣợng nguồn đổ nhƣng bị doi cát phía Đơng cản trở dịng chảy Bờ sơng bị sạt lở đoạn lớn uy hiếp đến hồ nuôi tôm, nhân dân tạo mỏ hàn tạm thời để nắn dịng chảy sang phía Nam Cho đến có số nghiên cứu cửa Đà Diễn có số giải pháp nhƣ: xây dựng kè, nạo vét luồng ngành giao thông, số hoạt động bảo vệ bờ chống sạt lở nhân dân quyền địa phƣơng Các giải pháp sở thực tế sạt lở, bồi lấp chƣa có nghiên cứu sâu định hƣớng hay quy hoạch tổng thể nên hiệu thấp Do đó, đề tài “Mơ trƣờng thủy động lực khu vực cửa Đà Diễn, tỉnh Phú Yên” đƣợc chọn làm nội dung luận văn thạc sĩ nhằm mô trƣờng thủy động lực mong nguồn tài liệu tham khảo góp phần vào việc xác định nguyên nhân, để từ định hƣớng, đƣa giải pháp nhƣ thiết kế công trình chỉnh trị cửa sơng Đà Diễn 10 Hình 3.22 Trường dịng chảy vùng cửa sơng KB_NETB (đỉnh triều) Hình 3.23 Trường sóng vùng cửa sơng KB_NETB thời kì đỉnh triều Hình 3.24 Trường sóng vùng cửa sơng KB_NETB thời kì chân triều 54 Khi triều lên, dịng triều sâu vào sông khoảng km với kịch KB_NE1% (Hình 3.22) gần km với kịch KB_NETB tính từ cửa (Hình 3.23) Dịng chảy lớn triều xuống đạt 1,9 m/s với kịch KB_NE1% 1,6 m/s kịch KB_NETB Với kịch sóng KB_NE1%, khu xốy khu vực trƣớc cửa sơng dịch vào bờ Bắc kích thƣớc khu nhỏ so với trạng tháng 11/2015 Vận tốc dịng chảy khu xốy trƣờng hợp (0,04 m/s) nhỏ lần so với vận tốc khu xoáy trƣờng hợp trạng tháng 11/2015 Xoáy sinh đƣợc sinh dòng triều kết hợp với địa hình (Hình 3.25) Hình 3.25 Vùng xốy trước cửa sơng KB_NE1% (triều lên) Hình 3.26 Vùng xốy trước cửa sông KB_NETB (triều rút) 55 Đối với kịch sóng KB_NETB, khu xốy xuất có vị trí gần nhƣ trạng tháng 11/2015 nhƣng vận tốc lớn khu xoáy 0,06 m/s nhỏ so với trạng 0,07 m/s lớn kịch sóng KB_NE1% 0,02 m/s Xốy sinh trƣờng hợp dòng lũ từ lạch Chùa đổ biển (Hình 3.26) Tại cửa sơng, so sánh hai kịch KB_NETB KB_NE1% nhận thấy, dòng chảy thấp so 0,02 m/s, mực nƣớc 0,05 m độ cao sóng có nghĩa 0,3 m Hƣớng sóng cửa khơng có thay đổi lớn hƣớng NE chiếm chủ đạo Ở phía bờ Bắc, so sánh kịch KB_NETB với kịch KB_NE1%, nhận thấy dòng chảy thấp 0,03 m/s, mực nƣớc thấp 0,08 m độ cao sóng có nghĩa thấp 0,6 m Hƣớng sóng bờ Bắc có biến động rõ rệt kịch KB_NE1% KB_NETB hƣớng sóng gần nhƣ khơng thay đổi Ở bờ Nam, so sánh kịch KB_NETB với kịch KB_NE1% nhận thấy dòng chảy thấp 0,01 m/s, mực nƣớc thấp 0,04 m, độ cao sóng có nghĩa thấp 0,17 m Cũng giống nhƣ bờ Bắc, hƣớng sóng bờ Nam có biến động rõ rệt kịch KB_NE1% KB_NETB hƣớng sóng gần nhƣ khơng thay đổi Ảnh hƣởng biến động sóng độ lớn sóng tới dịng chảy, mực nƣớc độ cao sóng bờ Bắc lớn bờ Nam 3.3.2 Kịch sóng Đơng (KB_E) Đối với nhóm kịch sóng Đơng, sóng từ biển vào gần nhƣ không đổi hƣớng vào bờ với kịch KB_E1% (Hình 3.27) Trong đó, với kịch KB_ETB sóng dần chuyển hƣớng từ Đơng sang Tây Bắc, vào khu vực cửa sơng sóng có hƣớng gần nhƣ vng góc với cửa (Hình 3.28) Chiều cao sóng lớn khu vực cửa sông 2,7 m với kịch KB_E1% 1,5 m với kịch KB_ETB Độ cao sóng từ biển vào với kịch KB_ETB biến đổi khơng sóng trƣờng hợp chịu tác động địa hình khu vực cửa Khi biên độ triều thay đổi, sóng khu vực có thay đổi độ lớn khu vực cửa Cụ thể với kịch sóng KB_E1%, độ cao sóng biến đổi ± 0,18 m, KB_ETB ± 0,06 m (Hình 3.29 - Hình 3.33) 56 Hình 3.27 Trường sóng miền tính tốn KB_E1% Hình 3.28 Trường sóng miền tính tốn KB_ETB Hình 3.29 Trường sóng khu vực cửa sơng KB_E1% (đỉnh triều) 57 Hình 3.30 Trường sóng khu vực cửa sơng KB_E1% (chân triều) Hình 3.31 Trường sóng khu vực cửa sơng KB_ETB (chân triều) Hình 3.32 Trường sóng khu vực cửa sơng KB_ETB (đỉnh triều) 58 Hình 3.33 Trường dịng chảy khu vực cửa sơng KB_E1% (đỉnh triều) Trong hai kịch sóng này, dịng chảy có giá trị nhỏ, xuất dịng ven bờ có hƣớng Tây Nam Cửa sơng chịu tác động chủ yếu dịng triều Khi triều lên, dòng ven kết hợp dòng triều tiến sâu vào khu vực cửa sơng phía bờ trái tới cầu Đà Rằng (cũ) Khi triều rút, dòng chảy hƣớng ngồi cửa sơng kết hợp với dịng ven bờ hƣớng phía Bắc, vận tốc lớn cửa sơng 0,8 m/s kịch sóng KB_E1% KB_ETB 0,6 m/s (Hình 3.34, Hình 3.36) Hình 3.34 Trường dịng chảy khu vực cửa sơng KB_E1% (chân triều) 59 Hình 3.35 Trường dịng chảy khu vực cửa sơng KB_ETB (đỉnh triều) Hình 3.36 Trường dịng chảy khu vực cửa sông KB_ETB (chân triều) Tại cửa sông, so sánh giá trị dòng chảy, mực nƣớc độ cao sóng có nghĩa hai kịch sóng KB_E1% KB_ETB Nhận thấy giá trị kịch sóng KB_E1% cao KB_ETB: dịng chảy 0,019 m/s, mực nƣớc 0,049 m, độ cao sóng có nghĩa 0,2 m Ở đây, khơng có thay đổi lớn hƣớng sóng Ở phía bờ Bắc, so sánh kịch KB_ETB với kịch KB_E1%, nhận thấy dòng chảy thấp 0,019 m/s, mực nƣớc thấp 0,05 m độ cao sóng có nghĩa thấp 0,3 m Hƣớng sóng bờ Bắc gần nhƣ khơng thay đổi hai kịch KB_E1% KB_ETB 60 Ở bờ Nam, so sánh kịch KB_NETB với kịch KB_NE1% nhận thấy dòng chảy thấp 0,02 m/s, mực nƣớc thấp 0,05 m, độ cao sóng có nghĩa thấp 0,29 m Hƣớng sóng hai kịch KB_NETB KB_NE1% khơng có thay đổi lớn Sự biến động sóng tới dịng chảy, mực nƣớc độ cao sóng bờ Nam lớn bờ Bắc 3.3.3 Sóng Nam (S) Trƣờng hợp kịch sóng KB_S1%, dịng chảy triều rút tiến sâu vào sơng qua lạch bờ Bắc (Hình 3.37) Sóng S từ từ đổi hƣớng tiến vào khu vực cửa Đà Diễn có hƣớng gần vng góc với bờ Ở khu vực cửa sơng, chiều cao sóng lớn 1,7 m Độ cao sóng chênh lệch khoảng ± 0,15 m triều lên triều rút (Hình 3.38, Hình 3.39) Đối với kịch sóng KB_STB tƣơng tự nhƣ kịch sóng KB_ S1% trạng Ở kịch này, dòng triều chiếm ƣu so với dòng chảy sơng (Hình 3.42, Hình 3.43), triều lên dịng triều vào sơng đến chân cầu Đà Rằng (cũ) Chênh lệch vận tốc khu vực cửa triều lên so với KB_ S1% 0,06 m/s Độ cao sóng lớn khu vực cửa gần m Hình 3.37 Trường sóng KB_S1% 61 Hình 3.38 Trường sóng KB_S1% (chân triều) Hình 3.39 Trường sóng KB_S1% (đỉnh triều) Hình 3.40 Trường dịng chảy KB_S1% (đỉnh triều) 62 Hình 3.41 Trường dịng chảy KB_S1% (chân triều) Hình 3.42 Trường sóng KB_STB (đỉnh triều) Hình 3.43 Trường sóng KB_STB (chân triều) 63 Hình 3.44 Trường dịng chảy KB_STB (đỉnh triều) Hình 3.45 Trường dịng chảy KB_STB (chân triều) Tại cửa sông, so sánh giá trị dịng chảy, mực nƣớc độ cao sóng có nghĩa hai kịch sóng KB_S1% KB_STB Nhận thấy giá trị kịch sóng KB_S1% cao KB_STB: dòng chảy 0,017 m/s, mực nƣớc 0,045 m, độ cao sóng có nghĩa 0,19 m Ở đây, khơng có thay đổi lớn hƣớng sóng Ở phía bờ Bắc, dịng chảy, mực nƣớc, độ cao sóng bờ bắc trƣờng hợp kịch sóng KB_STB thấp so với KB_S1% lần lƣợt 0,016 m/s, 0,04 m 0,4 m Hƣớng sóng hai kịch KB_S1% KB_STB có biến đổi lớn Ở phía bờ Nam, dịng chảy, mực nƣớc độ cao sóng có nghĩa trƣờng kịch KB_NETB thấp so với kịch KB_NE1% lần lƣợt 0,02 m/s, 0,05 m 0,29 m Có thay đổi lớn hƣớng sóng với KB_S1% KB_STB Sự biến đổi sóng có tác động lớn đến phía bờ Nam phía bờ Bắc kịch KB_S1% KB_STB 64 KẾT LUẬN VÀ KIẾN NGHỊ Kết luận Kết nghiên cứu “Mô trƣờng thủy động lực khu vực cửa Đà Diễn, tỉnh Phú Yên” thu đƣợc nhƣ sau: - Tổng hợp phân tích đƣợc điều kiện địa lý tự nhiên, khí tƣợng thủy văn kinh tế xã hội lƣu vực sông Ba khu vực cửa Đà Diễn Nhận thấy rằng, khu vực nơi quan trọng góp phần phát triển kinh tế xã hội cho tỉnh Phú Yên Nhƣng vùng cửa sông Đà Diễn bị biến động theo mùa (mùa lũ cửa sông đƣợc mở rộng ra, mùa kiệt cửa sông bị bồi lấp) dẫn đến hoạt động đánh bắt thủy sản giao thƣơng bị đình trệ Do đó, ảnh hƣởng trực tiếp đến đời sống ngƣời dân nhƣ mục tiêu phát triển kinh tế tỉnh - Mơ hình MIKE21 mơ hình có cở sở lý thuyết rõ ràng, bao gồm số module tính dịng chảy thơng dụng mơ tính tốn trƣờng dịng chảy, trƣờng sóng, vận chuyển bùn cát, chất lƣợng nƣớc cho vùng sông lục địa, cửa sông, vịnh vùng ven biển Mơ hình lƣới phi cấu trúc MIKE21FM tính vận chuyển bùn cát bồi xói vùng nghiên cứu diễn địa hình đáy bờ phức tạp, bị chia cắt - Số liệu sử dụng luận văn hoàn toàn đáng tin cậy, đƣợc cung cấp Đề tài ĐTĐL.CN 15/15 “Nghiên cứu sở khoa học để xác định chế bồi lấp, sạt lở đề xuất giải pháp ổn định cửa sông Đà Diễn Đà Nông tỉnh Phú Yên phục vụ phát triển bền vững sở hạ tầng kinh tế xã hội” - Kết hiệu chỉnh mơ hình tốt Các mơ dịng chảy cho thấy có phù hợp giá trị pha dao động số sai khác định việc phân tỷ lệ lƣu lƣợng biên chƣa sát với điều kiện thực tiễn - Kết kiểm định mơ hình đạt kết tốt, đƣờng trình mực nƣớc thực đo tính tốn tƣơng đối bám sát Tuy nhiên, yếu tố dịng chảy sóng, cịn sai khác tính tốn thực đo Địa hình khu vực cửa biến động mạnh mẽ sau mùa mƣa lũ, nguyên nhân sai khác 65 nghiên cứu sử dụng số liệu đo địa hình 11/2015 làm điều kiện tính tốn cho đợt 5/2016 Nhƣ vậy, chấp nhận đƣợc để sử dụng cho mô trƣờng thủy động lực vùng cửa Đà Diễn theo kịch tính tốn - Khu vực cửa Đà Diễn chịu ảnh hƣởng mạnh mẽ thủy triều Vào mùa khơ, triều lên sóng có khả tiến sâu vào cửa sơng, dịng chảy sóng có sức ảnh hƣớng so với dịng triều Dịng chảy sơng ngịi hầu nhƣ khơng có vai trị đáng kể, ngƣợc lại nhân tố động lực biển giữ vai trị chủ đạo Khu vực cửa sơng phía biển chủ yếu chịu ảnh hƣởng tác động dòng ven bờ từ phía Nam lên, kết hợp đồng pha với dịng triều Ở khu vực cửa phía sơng, dòng chảy khu vực chủ yếu dòng triều chi phối Trƣờng dịng chảy tồn vùng nghiên cứu đƣợc mơ có xét đến tƣơng tác sóng với dịng chảy cửa sơng dao động mực nƣớc triều Vào mùa lũ, dòng triều, dòng chảy lũ kết hợp với dịng chảy sóng ven bờ có hƣớng Tây Bắc - Đơng Nam Vùng cửa sơng chủ yếu chịu ảnh hƣởng mạnh mẽ dòng chảy sơng - Kịch sóng đƣợc xây dựng dựa vào chế độ sóng khu vực cửa Đà Diễn, với hƣớng sóng chủ đạo mùa lũ mùa kiệt, bao gồm sóng Nam (S), sóng Đơng (E) sóng Đơng Bắc (NE) - Kịch sóng thiết kế gồm sóng tần suất 1% sóng chế độ, tổ hợp với hƣớng sóng tạo thành kịch - Hƣớng sóng vào cửa sơng gần nhƣ vng góc với cửa kịch KB_NE, KB_E KB_S - Chiều cao sóng cửa sơng với kịch 1% lớn với sóng Đơng Bắc (NE) nhỏ với sóng Nam (S) - Dịng triều lên tới cầu Đà Rằng cũ kịch sóng mùa kiệt (KB_E, KB_S) Đối với kịch sóng mùa lũ (KB_NE), triều lên dịng triều sâu vào sông gần km kịch KB_NE1% gần km kịch KB_NETB 66 - Trƣớc cửa sông tồn khu xốy dịng triều vào kịch KB_NE1% triều lên Đối với kịch KB_NETB, khu xốy tồn dịng lũ từ Lạch Chùa đổ kết hợp với địa hình, phạm vị khu xốy vận tốc dịng chảy lớn KB_NE1% Các kịch sóng mùa kiệt, hồn lƣu trƣớc cửa sơng nhỏ không tồn - Khi triều lên triều xuống, không nhận thấy thay đổi hƣớng sóng kịch - Chênh lệch chiều cao sóng triều lên triều xuống với phƣơng án KB_NE1% 0,25 m; KB_NETB 0,15; KB_E1% 0,18 m; KB_ETB 0,06 m; KB_S1% 0,15 m KB_STB 0,04 m Kiến nghị Tuy nghiên cứu đạt đƣợc số kết định, xong kết mô chƣa thực tốt yếu tố biên đƣợc lấy từ kết tính tốn nghiên cứu khác, đồng thời khu vực nghiên cứu có địa hình biến đổi lớn Do đó, giai đoạn nghiên cứu cần phải xem xét lại việc sử dụng biên nghiên cứu rộng (tính tốn cụ thể tác động trình nạo vét đến chế độ thủy động lực cửa sông, ảnh hƣởng việc xây dựng công trình ) 67 TÀI LIỆU THAM KHẢO Tiếng Việt 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, Nhà xuất Tài Nguyên - Môi trƣờng đồ Việt Nam, Hà Nội Bùi Minh Hòa (2012), Nghiên cứu xây dựng đồ ngập lụt hạ lưu lưu vực sông Ba, Luận văn thạc sĩ khoa học, Trƣờng Đại học Khoa học Tự nhiên Phạm Thu Hƣơng (2012), Nghiên cứu sở khoa học cho việc đề xuất giải pháp ổn định cửa Đà Rằng, tỉnh Phú Yên, Luận án Tiến sĩ kĩ thuật, Trƣờng Đại học Thủy Lợi Hoàng Thị Hằng Nga (2013), Sử dụng mơ hình Mike 21 đánh giá q trình lan truyền nhiệt vùng cửa sông Trà Lý ảnh hưởng nhà máy nhiệt điện Thái Bình, Luận văn Thạc sĩ, Trƣờng Đại học Khoa học Tự nhiên Nguyễn Thọ Sáo (2003), Dự báo tượng xói lở, bồi tụ bờ biển cửa sông Đà Rằng, Báo cáo đề mục thuộc đề tài KC09.05/01 - 05 Lê Đình Thành, nnk (2010), Nghiên cứu đề xuất giải pháp ổn định cửa sông ven biển miền Trung, Đề tài KC08.07/06 - 10 Cấn Thu Văn (2010), Ứng dụng mơ hình MIKE - FLOOD tính tốn ngập lụt hạ lưu sông Ba, Luận văn thạc sĩ, Trƣờng Đại học Khoa học Tự nhiên Tiếng Anh DHI, (2011), MIKE21/3 Coupled model FM, User guide 68 ... trƣờng thủy động lực khu vực cửa Đà Diễn, tỉnh Phú Yên? ?? đƣợc chọn làm nội dung luận văn thạc sĩ nhằm mô trƣờng thủy động lực mong nguồn tài liệu tham khảo góp phần vào việc xác định nguyên nhân,... TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - Bùi Thị Hạnh MÔ PHỎNG TRƢỜNG THỦY ĐỘNG LỰC KHU VỰC CỬA ĐÀ DIỄN, TỈNH PHÚ YÊN Chuyên ngành: Thủy văn học Mã số: 60 44 02 24 LUẬN VĂN THẠC SĨ KHOA HỌC... HÌNH MƠ PHỎNG TRƢỜNG THỦY ĐỘNG LỰC CỬA SƠNG ĐÀ DIỄN 3.1 Xây dựng kịch Qua khảo sát thực tế q trình nghiên cứu, dịng chảy sơng Ba có tác động lớn tới chế độ thủy động lực khu vực cửa sông Đà Diễn