(Luận văn thạc sĩ) đánh giá tác động dòng chảy hạ lưu sông ba đến sự bồi lấp cửa sông đà 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 - PHẠM QUỐC SỸ ĐÁNH GIÁ TÁC ĐỘNG DỊNG CHẢY HẠ LƢU SƠNG BA ĐẾN SỰ BỒI LẤP CỬA SÔNG ĐÀ DIỄN, TỈNH PHÚ YÊN LUẬN VĂN THẠC SĨ KHOA HỌC Hà Nội – Năm 2018 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - PHẠM QUỐC SỸ ĐÁNH GIÁ TÁC ĐỘNG DỊNG CHẢY HẠ LƢU SƠNG BA ĐẾN SỰ BỒI LẤP CỬA SÔNG ĐÀ DIỄN, TỈNH PHÚ YÊN Chuyên ngành: Thủy văn học Mã số: 60440224 LUẬN VĂN THẠC SĨ KHOA HỌC NGƯỜI HƯỚNG DẪN KHOA HỌC: PGS TS NGYỄN THỌ SÁO Hà Nội – Năm 2018 LỜI CẢM ƠN Sau thời gian nghiên cứu, em hoàn thành luận văn thạc sĩ khoa học “Đánh giá tác động dòng chảy hạ lưu sông Ba đến bồi lấp cửa sông Đà Diễn, tỉnh Phú Yên” Em xin chân thành gửi lời cảm ơn sâu sắc đến PGS.TS Nguyễn Thọ Sáo, thầy giáo, cán Khoa Khí tượng Thủy văn Hải dương học - Trường Đại học Tự nhiên, Đại học Quốc Gia Hà Nội tận tình hướng dẫn, ân cần truyền đạt kiến thức, tạo điều kiện thuận lợi để em có hội học tập, nghiên cứu, phát triển kỹ năng, hoàn thành chương trình học Em xin cảm ơn thầy cô Ban chủ nhiệm đề tài cấp nhà nước “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” mã số ĐTĐL.CN.15/15 Trường Đại học Khoa học Tự nhiên, Đại học Quốc gia Hà Nội chủ trì, cung cấp số liệu, tài liệu hỗ trợ em nhiều trình học tập, nghiên cứu hoàn hoàn thành luận văn Em xin cảm ơn anh chị em nhóm “G‟Group” giúp đỡ q trình hồn thành luận văn Xin cảm ơn gia đình, bạn bè người thân ln bên cạnh, giúp đỡ động viên em vượt qua khó khăn thời gian qua Luận văn có nhiều cố gắng khó tránh khỏi thiếu sót Em mong nhận nhận xét, góp ý quý báu quý thầy cô, bạn bè, đồng nghiệp, nhà khoa học độc giả để luận văn em hoàn thiện Xin trân trọng cảm ơn! Tác giả Phạm Quốc Sỹ MỤC LỤC Chƣơng ĐIỀU KIỆN ĐỊA LÝ TỰ NHIÊN VÀ KINH TẾ XÃ HỘI KHU VỰC NGHIÊN CỨU 1.1 Điều kiện tự nhiên khu vực nghiên cứu 1.1.1 Điều kiện tự nhiên .4 1.1.2 Đặc điểm khí tượng khí hậu 12 1.1.3 Đặc điểm thủy văn 21 1.1.4 Đặc điểm hải văn…………… ……………………………………… 22 1.2 Điều kiện kinh tế xã hội 24 1.2.1 Tổ chức hành lưu vực sông Ba 24 1.2.2 Dân số .25 1.2.3 Đặc điểm kinh tế 25 1.3 Đánh giá tổng hợp yếu tố động lực ảnh hƣởng đến diễn biến cửa Đà Rằng 26 Chƣơng CƠ SỞ LÝ THUYẾT MƠ HÌNH MIKE 11 VÀ MƠ HÌNH MIKE 21……………… 28 2.1 Giới thiệu chung 28 2.2 Mơ hình Mike 11 28 2.2.1 Giới thiệu phần mềm Mike 11 28 2.2.2 Hệ phương trình sử dụng phần mềm Mike 11 30 2.3 Mơ hình Mike 21 31 2.3.1 Giới thiệu phần mềm Mike 21 31 2.3.2 Hệ phương trình sử dụng phần mềm Mike 21 32 Chƣơng ĐÁNH GIÁ TÁC ĐỘNG DỊNG CHẢY HẠ LƢU SƠNG BA ĐẾN SỰ BỒI LẤP CỬA SÔNG ĐÀ DIỄN, TỈNH PHÚ YÊN 40 3.1 Phạm vi tính tốn 40 3.2 Kết hiệu chỉnh kiểm định mơ hình Mike 11 41 3.2.1 Thiết lập thơng số mơ hình Mike 11 41 i 3.2.2 Kết hiệu chỉnh mơ hình chiều 45 3.2.3 Kết kiểm định mơ hình chiều: 46 3.3 Kết hiệu chỉnh kiểm định mơ hình Mike 21 49 3.3.1 Số liệu phục vụ nghiên cứu .49 3.3.2 Thiết lập thơng số mơ hình Mike 21 51 3.3.3 Kết hiệu chỉnh mơ hình hai chiều năm 2015 52 3.3.4 Kết kiểm định mơ hình hai chiều năm 2016 .52 3.4 Kết tính tốn theo kịch 53 3.4.1 So sánh trường dòng chảy pha triều xuống thời kì tháng khác nhau.…… .55 3.4.2 So sánh trường sóng theo hướng sóng khác Bắc, Đông từ tháng đến tháng 58 3.4.3 So sánh thay đổi đáy 15 ngày với hướng sóng khác nhau…………… 59 3.4.4 So sánh thay đổi đáy với hướng sóng 65 TÀI LIỆU THAM KHẢO 73 Tài liệu Tiếng Việt 73 Tài liệu Tiếng Anh 73 ii DANH MỤC BẢNG Bảng 1.1: Các đặc trưng sơng Ba số sông lưu vực[2] 11 Bảng 1.2: Tốc độ gió trung bình tháng năm (Đơn vị: m/s)[2] .13 Bảng 1.3: Một số đặc trưng mưa năm (Đơn vị: mm) [2] 16 Bảng 1.4: Lượng mưa trung bình nhiều năm tháng (Đơn vị: mm)[2] 18 Bảng 1.5: Độ ẩm tuyệt đối trung bình tháng năm(mb) [2] 20 Bảng 1.6: Độ ẩm tương đối trung bình tháng năm (Đơn vị: %)[2] 20 Bảng 1.7: Lưu lượng lũ lớn số trạm lưu vực sông Ba [2] 21 Bảng 1.8 Bảng tính tốn cao độ biên độ thủy triều dựa số liệu toàn cầu[2]25 Bảng 1.9: Độ cao sóng bình qn (m) theo tháng mùa năm[2]…… 26 Bảng 3.1: Các kịch theo hướng sóng khác 53 iii DANH MỤC HÌNH Hình 1.1: Sơ đồ hành khu vực tỉnh Phú Yên .5 Hình 1.2: Địa hình vùng cửa sông Đà Rằng .6 Hình 1.3: Bản đồ mạng lưới sơng lưu vực sông Ba 10 Hình 1.4: Sơ đồ vùng hạ lưu sơng Ba từ Củng Sơn đến cửa Đà Rằng .12 Hình 1.5: Phân bố lượng mưa theo mùa mưa mùa khơ 19 Hình 2.1 Sơ đồ lưới tính so le theo hai chiều x, y………………………………….34 Hình 2.2 Sơ đồ quét thời gian trung tâm………………………………………… 35 Hình 3.1: Phạm vi tính tốn mơ hình Mike 11 40 Hình 3.2: Phạm vi tinh tốn mơ hình Mike 21 40 Hình 3.3: Lưới tính toán luận văn 41 Hình 3.4: Mạng lưới sơng 42 Hình 3.5: Mặt cắt ngang sơng 43 Hình 3.6: Biên lưu lượng Củng Sơn .45 Hình 3.7: Biên mực nước dựa theo tương quan thủy triều 45 Hình 3.8: So sánh mực nước tính tốn thực đo trạm Phú Lâm .46 Hình 3.9 Biên lưu lượng Củng Sơn……………………………………… 47 Hình 3.10: Biên mực nước dựa theo tương quan thủy triều 48 Hình 3.11: Kết kiểm định mơ hình trạm Phú Lâm 49 Hình 3.12: Địa hình khu vực tháng 3/2016 khu vực tháng 9/2016 50 Hình 3.13: Tọa độ yếu tố đo đạc trạm đo cụ thể .50 Hình 3.14: Điều kiện biên mực nước 51 Hình 3.15: So sánh mực nước trạm D 52 Hình 3.16: So sánh mực nước trạm F 52 Hình 3.17: Hoa sóng tính từ số liệu gió đo trạm Tuy Hịa 54 Hình 3.18 : KB1 Trường dịng chảy từ tháng đến tháng pha triều xuống theo hướng sóng Bắc 55 Hình 3.19: KB4 Trường dòng chảy tháng đến tháng pha triều xuống theo hướng sóng Bắc 56 iv Hình 3.20: KB7 Trường dịng chảy tháng 10 đến tháng 12 pha triều xuống theo hướng sóng Bắc 56 Hình 3.21: KB 3Trường sóng theo hướng sóng Đơng từ tháng đến tháng 58 Hình 3.22: KB2 Trường sóng theo hướng sóng Đơng Bắc từ tháng đến tháng 58 Hình 3.23: KB1Trường sóng theo hướng sóng Bắc từ tháng đến tháng .59 Hình 3.24: KB Biến đổi đáy theo hướng sóng Đơng từ ngày 6/4/2016 đến ngày 21/4/2016 60 Hình 3.25: KB Biến đổi đáy theo hướng sóng Đơng Bắc từ ngày 6/4/2016 đến ngày 21/4/2016 60 Hình 3.26: KB Biến đổi đáy theo hướng sóng Bắc từ ngày 6/4/2016 đến ngày 21/4/2016 61 Hình 3.27: KB Biến đổi đáy theo hướng sóng Đơng từ ngày 9/5/2016 đến ngày 23/5/2016 62 Hình 3.28: KB Biến đổi đáy theo hướng sóng Đơng Bắc từ ngày 9/5/2016 đến ngày 23/5/2016 62 Hình 3.29: KB Biến đổi đáy theo hướng sóng Bắc từ ngày 9/5/2016 đến ngày 23/5/2016 63 Hình 3.30: KB9 Biến đổi đáy theo hướng sóng Đơng từ ngày 12/10/2016 đến ngày 26/10/2016 .64 Hình 3.31: KB8 Biến đổi đáy theo hướng sóng Đơng Bắc từ ngày 12/10/2016 đến ngày 26/10/2016 64 Hình 3.32: KB Biến đổi đáy theo hướng sóng Bắc từ ngày 12/10/2016 đến ngày 26/10/2016 .65 Hình 3.33: Biến đổi đáy theo hướng sóng Đơng sau 10 ngày KB3 (ngày 16/4/2016) 66 Hình 3.34: Biến đổi đáy theo hướng sóng Đơng sau 10 ngày KB6 ( ngày 23/5/2016 66 Hình 3.35: Biến đổi đáy theo hướng sóng Đơng sau 10 ngày KB 9(ngày 22/10/2016) .67 v Hình 3.36: Biến đổi đáy theo hướng sóng Đơng Bắc sau 10 ngày KB (ngày 16/4/2016) 67 Hình 3.37: Biến đổi đáy theo hướng sóng Đông Bắc sau 10 ngày KB5 ( ngày 23/5/2016) 68 Hình 3.38: Biến đổi đáy theo hướng sóng Đơng Bắc sau 10 ngày KB8 (ngày 22/10/2016) .68 Hình 3.39: Biến đổi đáy theo hướng sóng Bắc sau 10 ngày KB1 (ngày 16/4/2016) 69 Hình 3.40: Biến đổi đáy theo hướng sóng Bắc sau 10 ngày KB4 (ngày 23/5/2016) 69 Hình 3.41: Biến đổi đáy theo hướng sóng Bắc sau 10 ngày KB7 (ngày 22/10/2016) .70 vi MỞ ĐẦU Việt Nam nằm khu vực nhiệt đới gió mùa, có đường bờ biển khoảng 3260 km (chưa kể bờ đảo) 29 tỉnh thành phố, có mật độ lưới sông vào loại lớn giới, với mức trung bình km/km2, 20 km bờ biển có cửa sơng Diện tích vùng ven biển Việt Nam 136.700 km2 (chiếm 41% diện tích nước), với dân số tính đến năm 2010 khoảng 42 triệu người (chiếm 48% nước) Các ngành kinh tế biển ven biển khai thác dầu khí, khống sản biển, đánh bắt ni trồng thủy sản, vận tải đường thủy, du lịch,… có vai trò quan trọng chiếm tỷ lệ đáng kể kinh tế Việt Nam Hàng năm kinh tế vùng ven biển đóng góp khoảng 30% GDP 50% tổng thu nhập xuất nước [2] Dải ven biển nước ta phong phú tài nguyên, nơi tập trung dân cư đông đúc so với vùng khác nơi có nhiều cơng trình dân sinh kinh tế, quốc phịng quan trọng Tuy nhiên, vùng hứng chịu nhiều thiên tai bão, nước dâng, sạt lở bồi tụ Bồi tụ bờ biển, cửa sông nước ta tạo nên bãi bồi quý giá cho nhiều vùng, nhiều nơi lại gây hậu nghiêm trọng bồi lấp luồng tàu, bến cảng Các cửa sông ven biển miền Trung phần lớn không ổn định, thường xuyên bị phá xói lở (mở ra) vào mùa mưa, mùa khơ lại bị bồi lấp (đóng lại) có xu hướng dịch chuyển theo chu kỳ cửa Thuận An, cửa Tư Hiền, cửa Lở, cửa Đà Nông… Hậu nghiêm trọng bồi lấp cửa sơng làm giảm khả lũ, gây ngập lụt diện rộng, làm hóa đầm phá, gây ách tắc tầu thuyền ngư dân vào cửa sơng Cịn xói lở gây nhiều thiệt hại đến cơng trình bến cảng, đê kè nơi người dân Có thể nói dạng thiên tai có diễn biến phức tạp, gây thiệt hại lớn để lại hậu lâu dài kinh tế, xã hội môi trường Cửa sơng đóng vai trị quan trọng việc thơng thương với giới bên ngồi Từ xa xưa, cửa sơng đóng vai trị quan trọng đời sống người Nhiều cửa sông lớn trở thành nôi văn minh cổ đại Chính có vai trị quan trọng mà từ lâu cửa sông đối tượng nghiên Hình 3.23: KB1 Trƣờng sóng theo hƣớng sóng Bắc từ tháng đến tháng Với tháng đến tháng 10 đến 12 tương tự Ta nhận thấy theo hướng sóng Bắc dịch xuống phía cửa, cịn theo hướng sóng Đơng Bắc hướng sóng Đơng thẳng vào cửa gây nguy hiểm Trong hướng sóng Đơng Bắc chiếm ưu hướng sóng Bắc Đông độ cao lẫn tần suất xuất Do hướng sóng Đơng Bắc gần vng góc với đường bờ nên vận chuyển bùn cát ngang bờ tác động hướng sóng lớn Đây nguyên nhân gây biến động cửa Đà Rằng Qua kết mơ trường sóng cho thấy, vùng ven biển cửa Đà Rằng chịu ảnh hưởng mạnh mẽ thuỷ triều lên xuống vào mùa kiệt Khu vực nghiên cứu ngày, triều dâng sóng có khả tiến sâu vào cửa triều xuống sóng cách bờ Do cửa Đà Rằng có hướng vng góc với hướng Đơng Bắc nên mùa lũ, sóng có hướng tác động trực tiếp vào cửa sơng, chiều cao sóng trung bình cửa sông khoảng 0,5m 3.4.3 So sánh thay đổi đáy 15 ngày với hƣớng sóng khác 3.4.3.1 Sự biến đổi đáy sau 15 ngày từ ngày 6/4/2016 đến ngày 21/4/2016: Biến đổi đáy đặc trưng tháng đến tháng 59 Ngày 21/4/2016 Hình 3.24: KB Biến đổi đáy theo hƣớng sóng Đơng từ ngày 6/4/2016 đến ngày 21/4/2016 Hình 3.25: KB Biến đổi đáy theo hƣớng sóng Đơng Bắc từ ngày 6/4/2016 đến ngày 21/4/2016 60 Hình 3.26: KB Biến đổi đáy theo hƣớng sóng Bắc từ ngày 6/4/2016 đến ngày 21/4/2016 Qua hình ảnh cho ta thấy tháng đến tháng 4: - Sóng theo hướng Bắc làm cho doi cát bị chia cắt mạnh bị phá vỡ mảng theo hướng sóng Đơng Bắc Đơng doi cát chạy song song chắn ngang cửa khiến cửa bị đóng thu hẹp lại - Độ biến động đáy phía bờ trái 0.24-0.28m, Độ biến động đáy phía bờ phải 0.20-0.28m Độ biến động đáy luồng cửa 0.04-0.08m Ta nhận thấy với ngày sóng theo hướng Bắc làm xói lở đáy lớn - Vào mùa kiệt, trường dòng chảy nhỏ, cấp lưu lượng nhỏ nên biến đổi đáy thời kì nhỏ 3.4.3.2 Sự biến đổi đáy sau 15 ngày từ ngày 9/5/2016 đến ngày 23/5/2016: Biến đổi đáy đặc trưng tháng đến tháng Ngày 23/5/2016 61 Hình 3.27: KB Biến đổi đáy theo hƣớng sóng Đơng từ ngày 9/5/2016 đến ngày 23/5/2016 Hình 3.28: KB Biến đổi đáy theo hƣớng sóng Đơng Bắc từ ngày 9/5/2016 đến ngày 23/5/2016 62 Hình 3.29: KB Biến đổi đáy theo hƣớng sóng Bắc từ ngày 9/5/2016 đến ngày 23/5/2016 - Như ngày 23/5/2016 bồi lớn 0.24m xói thấp -0,32m theo hướng sóng - Cũng tương tự tháng đến tháng 4, biến đổi đáy theo hướng sóng Bắc xiên thẳng vào xuống phía cửa làm doi cát bị chia cắt mạnh phân mảng cịn biến đổi đáy theo hướng sóng Đơng Đông Bắc xiên thẳng vào cửa doi cát chạy song song chắn ngang cửa - Thời kì mùa kiệt, cấp lưu lượng nhỏ, trường dòng chảy tương đối nhỏ nên biến đổi đáy nhỏ không lớn 3.4.3.3 Sự biến đổi đáy sau 15 ngày từ ngày 12/10/2016 đến ngày 26/10/2016 Biến đổi đáy đặc trưng tháng 10 đến tháng 12 Ngày 26/10/2016 63 Hình 3.30: KB9 Biến đổi đáy theo hƣớng sóng Đơng từ ngày 12/10/2016 đến ngày 26/10/2016 Hình 3.31: KB8 Biến đổi đáy theo hƣớng sóng Đơng Bắc từ ngày 12/10/2016 đến ngày 26/10/2016 64 Hình 3.32: KB Biến đổi đáy theo hƣớng sóng Bắc từ ngày 12/10/2016 đến ngày 26/10/2016 Qua hình ảnh trên, tương tự tháng đến tháng tháng đến tháng 9, ta thấy biến đổi đáy 15 ngày lớn với hướng sóng bắc doi cát bị chia cắt mạnh, tác động sóng tới trình vận chuyển bùn cát đáng kể cịn hướng sóng đơng đơng bắc doi cát song song chạy chắn ngang cửa Với tháng 10 đến tháng 12 mùa lũ, cấp lưu lượng lớn, trường dòng chảy lớn khiến cho tương bồi xói lớn so với tháng đến tháng tháng đến tháng mùa kiệt Như ngày 26/10/2016 đáy bị bồi lớn 0.8m xói nhỏ -0.4m - Độ biến động đáy bờ trái 0.72-0.8m, Độ biến động đáy bờ phải 0.00-0.08m - Độ biến động đáy luồng cửa 0.08-0.16m 3.4.4 So sánh thay đổi đáy với hƣớng sóng 3.4.4.1 Cùng theo hướng sóng Đơng 65 Hình 3.33: Biến đổi đáy theo hƣớng sóng Đơng sau 10 ngày KB3 (ngày 16/4/2016) Hình 3.34: Biến đổi đáy theo hƣớng sóng Đơng sau 10 ngày KB ( ngày 23/5/2016) 66 Hình 3.35: Biến đổi đáy theo hƣớng sóng Đơng sau 10 ngày KB (ngày 22/10/2016) Qua hình ảnh ta thấy rằng: Ngày 22/10/2016 theo hướng sóng Đơng bị xói mạnh, xói sâu vào với ngày 16/4/2016 ngày 23/5/2016 3.4.4.2 Cùng theo hướng sóng Đơng Bắc Theo hướng Đông Bắc thay đổi đáy tương tự Hình 3.36: Biến đổi đáy theo hƣớng sóng Đơng Bắc sau 10 ngày KB (ngày 16/4/2016) 67 Hình 3.37: Biến đổi đáy theo hƣớng sóng Đơng Bắc sau 10 ngày KB (ngày 23/5/2016) Hình 3.38: Biến đổi đáy theo hƣớng sóng Đơng Bắc sau 10 ngày KB (ngày 22/10/2016) 68 Qua hình ảnh (hình 3.36, 3.37, 3.38) chụp vùng cửa sông trên, ta thấy vùng cửa sông ngày 22/10/2016 bị xói so với ngày 16/4/2016 ngày 23/5/2016 3.4.4.2 Cùng theo hướng sóng Bắc Hình 3.39: Biến đổi đáy theo hƣớng sóng Bắc sau 10 ngày KB (ngày 16/4/2016) Hình 3.40: Biến đổi đáy theo hƣớng sóng Bắc sau 10 ngày KB (ngày 23/5/2016) 69 Hình 3.41: Biến đổi đáy theo hƣớng sóng Bắc sau 10 ngày KB (ngày 22/10/2016) Qua hình ảnh trên, ta rút nhận xét: với hướng sóng Bắc ngày 22/10/2016 xói lở mạnh hơn, xói sâu so với ngày 23/5/2016 ngày 16/4/2016 Ta kết luận theo hướng sóng Bắc theo tháng 10 đến tháng 12 mùa lũ xói mạnh, xói sâu so với tháng đến tháng tháng đến tháng mùa kiệt, tượng bồi xói vào tháng 10 đến tháng 12 phức tạp Như theo hướng sóng Bắc gây tượng bồi xói nghiêm trọng so với hướng sóng Đơng Đơng Bắc 70 KẾT LUẬN Nghiên cứu diễn biến cửa sông vấn đề khó phức tạp, đặc biệt cửa sông ven biển miền Trung, nơi thường xuyên có thiên tai bão, lũ Cửa Đà Rằng cửa sông lớn miền Trung với đầy đủ yếu tố tự nhiên kinh tế, xã hội tác động đến diễn biến Quá trình nghiên cứu tổng quan điều kiện tự nhiên, điều kiện khí tượng thủy hải văn kinh tế xã hội, trạng khu vực cửa sông Đà Diễn, tỉnh Phú Yên nghiên cứu, tính tốn xây dựng mơ hình tốn diễn biến hình thái cửa sơng Đà Diễn việc xác định nguyên nhân chế bồi lấp cửa sông cho phép rút số đánh giá tác động dịng chảy hạ lưu sơng Ba đến bồi lấp cửa sông Đà Diễn, tỉnh Phú Yên sau: Đánh giá trường dịng chảy sơng: Dịng chảy từ sơng đóng vai trị để phân bố lại bùn cát từ sông cửa tải bùn cát sóng bứt khỏi bờ đáy khu gần bờ, dịng dư (dịng từ sơng đổ dịng sóng) tạo nên dịng chảy ven bờ tổng hợp vận chuyển bùn cát dọc bờ Đánh giá trường dòng chảy kết hợp với điều kiện thủy động lực ven biển - Vào mùa kiệt dòng chảy sơng ngịi khơng đáng kể, ngược lại nhân tố động lực biển giữ vai trò chủ đạo trình biến động phát triển bồi tụ xói lở cửa sơng Vùng ven biển cửa sơng Đà Rằng chịu ảnh hưởng mạnh mẽ thủy triều lên xuống Ở 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 nên dòng triều lên lượng bùn cát mang từ biển vào bồi lắng lại vùng cửa sơng Vì vậy, thời kỳ này, cửa sơng bị bồi lấp thu hẹp lại, bãi bồi cửa sơng xuất hiện, cửa Đà Rằng có thời gian bị bồi gần lấp kín, bồi tụ nhiều khu vực cửa sông vào khoảng 0,25 – 0.8 m, lớn lên tới 0.8 m - Vào mùa lũ, dòng triều xuống kết hợp với dòng lũ từ sông đổ mang lượng lớn bùn cát phía ngồi cửa sơng, lắng đọng chủ yếu khu vực ngồi cửa sơng, bãi cát sơng dịng lũ đẩy phía biển Tuy nhiên, phía sơng có nhiều nơi bị xói mạnh, dịng chảy lũ sóng vào sâu gây ngập lụt khu vực cửa sông Đà Rằng, cửa sông mở rộng Khu vực cửa sơng xói lở, độ sâu xói lớn lên tới 71 -0.056m Phía ngồi cửa sơng xuất bãi bồi Lượng bùn cát lắng đọng chủ yếu khu vực ngồi cửa sơng vào mùa lũ (tháng 11) Ta nhận thấy lịng sơng xói sâu mạnh vào mùa lũ Đánh giá trường sóng: Với cấp lưu lượng, hướng sóng Bắc gây nguy hiểm nhất, hướng sóng thẳng xuống phía cửa, doi cát vỡ mảng; cịn hướng sóng Đơng, Đơng Bắc doi cát song song chắn ngang cửa Nghiên cứu sở cho nghiên cứu sâu việc xác định nguyên nhân chế bồi lấp cửa sông Đà Diễn, tỉnh Phú yên đề xuất giải pháp chỉnh trị Luận văn đưa hướng nghiên cứu lĩnh vực xây dựng mơ hình tốn sau: Mơ hình tốn luận văn sử dụng làm sở để xây dựng mơ hình mơ diễn biến hình thái cửa sơng Đà Diễn, tỉnh Phú Yên Kiến nghị: Trong luận văn bước đầu thiết lập mơ hình mơ đánh giá diễn biến vùng cửa sơng Đà Diễn, có xét đến tác động dịng chảy hạ lưu sơng Ba nhiên kết đánh giá mơ hình mơ cịn nhiều hạn chế cần cải tiến thêm nghiên cứu tiếp sau 72 TÀI LIỆU THAM KHẢO Tài liệu Tiếng Việt Phạm Duy Huy Bình (2017), Ứng dụng phương pháp hồi quy đa biến logistic xây dựng mơ hình nhận thức diễn biến hình thái khu vực cửa sông Đà Diễn, tỉnh Phú Yên, Luận văn Thạc sĩ khoa học chuyên ngành Thủy văn học, trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội Bùi Minh Hòa (2012), Nghiên cứu xây dựng đồ ngập lụt hạ lưu sông Ba, Luận văn Thạc sĩ ngành Thủy văn học, trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội 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 chuyên ngành Chỉnh trị sông bờ biển, trường Đại học Thủy lợi - Bộ Nông nghiệp &PTNT Nguyễn Tiền Giang (2015), 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, Thuyết minh đề tài cấp nhà nước mã số ĐTĐL.CN.15/15 Lê Đình Thành, Nguyễn Bá Quỳ 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, Hà nội 2010 Tài liệu Tiếng Anh MIKE - 11 reference/ scientific document MIKE - 21 reference/ scientific document 73 ... NHIÊN - PHẠM QUỐC SỸ ĐÁNH GIÁ TÁC ĐỘNG DỊNG CHẢY HẠ LƢU SƠNG BA ĐẾN SỰ BỒI LẤP CỬA SÔNG ĐÀ DIỄN, TỈNH PHÚ YÊN Chuyên ngành: Thủy văn học Mã số: 60440224 LUẬN VĂN THẠC SĨ KHOA HỌC NGƯỜI... trình chuyển động hạt bùn cát giải s (từ nhận đến điểm lắng đọng) đánh giá xác 39 Chƣơng ĐÁNH GIÁ TÁC ĐỘNG DỊNG CHẢY HẠ LƢU SƠNG BA ĐẾN SỰ BỒI LẤP CỬA SÔNG ĐÀ DIỄN, TỈNH PHÚ YÊN 3.1 Phạm vi tính... luận văn thạc sĩ khoa học ? ?Đánh giá tác động dịng chảy hạ lưu sơng Ba đến bồi lấp cửa sông Đà Diễn, tỉnh Phú Yên? ?? Em xin chân thành gửi lời cảm ơn sâu sắc đến PGS.TS Nguyễn Thọ Sáo, thầy giáo,