Nghiên cứu chế độ thủy động lực cửa tiên châu tỉnh phú yên

LỜI CAM ĐOAN Tôi - Mai Duy Khánh xin cam đoan cơng trình nghiên cứu thân học viên Các kết nghiên cứu kết luận luận văn trung thực, không chép từ nguồn hình thức Việc tham khảo nguồn tài thực trích dẫn ghi nguồn tài liệu tham khảo quy định Tác giả luận văn Mai Duy Khánh i LỜI CÁM ƠN Với tất tình cảm mình, tơi xin bày tỏ lịng kính trọng biết ơn sâu sắc tới PGS.TS Trần Thanh Tùng, TS Lê Tuấn Hải, hai người thầy tận tình bảo, dìu dắt hướng dẫn cho tơi suốt trình nghiên cứu thực luận văn Tôi xin chân thành cảm ơn thầy giáo Kỹ thuật cơng trình biển đường thủy tạo điều kiện cho suốt q trình học tập nghiên cứu Tơi xin trân trọng cảm ơn Ban Giám hiệu, Phòng Đào tạo (Trường Đại học Thủy lợi), tạo điều kiện thuận lợi cho tơi hồn thành luận văn Đồng thời, tơi xin gửi lời cảm ơn tới Ban chủ nhiệm đề tài Độc lập cấp Nhà nước “Nghiên cứu giải pháp chỉnh trị chống sa bồi luồng tàu cho cảng cá khu neo đậu tàu thuyền tỉnh Phú Yên vùng lân cận, áp dụng cho cửa Tiên Châu” mã số ĐTĐL.CN.33/18 cho phép sử dụng số liệu, kết nghiên cứu đề tài Đây nguồn tài liệu quý giá giúp thực luận văn Xin gửi lời cảm ơn chân thành tới cha mẹ-những người sinh tôi, điểm tựa vững cho Tôi xin gửi lời cảm ơn tới vợ – người bạn đồng hành lúc khó khăn Cuối tơi xin chân thành cảm ơn Ban lãnh đạo anh, chị bạn đồng nghiệp Trung tâm Tư vấn kỹ thuật biển Phát triển cảng - Viện Kỹ thuật công trình ln ủng hộ, động viên giúp đỡ tơi suốt q trình thực luận văn Xin trân trọng cảm ơn! Học viên Mai Duy Khánh ii MỤC LỤC DANH MỤC HÌNH ẢNH v DANH MỤC BẢNG viii DANH MỤC TỪ VIẾT TẮT ix MỞ ĐẦU CHƯƠNG TỔNG QUAN KHU VỰC NGHIÊN CỨU 1.1 Giới thiệu khu vực nghiên cứu 1.1.1 Vị trí địa lý, đặc điểm địa hình, địa mạo khu vực nghiên cứu 1.1.2 Đặc điểm khí tượng, khí hậu [2] 1.1.3 Đặc điểm thủy văn, hải văn, thủy động lực khu vực nghiên cứu [2] 10 1.1.4 Đặc điểm dân sinh, kinh tế – xã hội [4] 16 1.2 Hiện trạng cơng trình khu vực nghiên cứu 19 1.2.1 Cảng Tiên Châu khu nước quanh cảng 19 1.2.2 Các cơng trình đập dâng lưu vực sông Kỳ Lộ 22 1.2.3 Các cơng trình bảo vệ bờ sơng Kỳ Lộ 23 1.3 Khái quát chung nghiên cứu cửa sông giới Việt Nam 32 1.3.1 Phân loại cửa sông giới Việt Nam 32 1.3.2 Khái quát chung nghiên cứu cửa sông Việt Nam miền Trung 35 1.4 Kết luận chương 39 CHƯƠNG NGHIÊN CỨU QUY LUẬT DIỄN BIẾN DOI CÁT VEN BỜ PHÍA BẮC CỬA TIÊN CHÂU BẰNG TƯ LIỆU ẢNH VIỄN THÁM [28] 40 2.1 Cơ sở liệu 40 2.2 Giải đoán đường bờ biển từ ảnh viễn thám 41 2.3 Phân tích quy luật phát triển doi cát 47 2.3.1 Phân tích lịch sử diễn biến doi cát cửa Tiên Châu 48 2.3.2 Phân tích tương quan diễn biến doi cát động lực sông 50 2.4 Kết luận chương 55 CHƯƠNG THIẾT LẬP MƠ HÌNH THỦY ĐỘNG LỰC KHU VỰC CỬA TIÊN CHÂU… 56 3.1 Giới thiệu mơ hình thủy động lực Delft-3D 56 3.2 Cơ sở lý thuyết phương pháp tính 58 3.2.1 Mô đun tính tốn thủy động lực học [32] 58 3.2.2 Mơ đun tính tốn sóng [33] 61 3.3 Tài liệu sử dụng nghiên cứu 64 3.3.1 Tài liệu địa hình đường bờ 64 3.3.2 Tài liệu mực nước, sóng lưu lượng 65 3.4 Thiết lập mơ hình thủy động lực học 67 3.4.1 Xác định miền tính lưới tính 67 3.4.2 Thiết lập địa hình 69 iii 3.4.3 Thiết lập điều kiện biên mơ hình thủy động lực học 70 3.4.4 Thiết lập điều kiện ban đầu mơ hình thủy động lực học 72 3.4.5 Thiết lập thông số thủy lực 72 3.4.6 Hiệu chỉnh kiểm định mô hình thủy động lực 74 3.5 Thiết lập mơ hình sóng 79 3.5.1 Xác định miền tính tốn 79 3.5.2 Xây dựng lưới tính 80 3.5.3 Thiết lập địa hình 81 3.5.4 Thiết lập điều kiện biên mơ hình sóng 81 3.5.5 Thiết lập điều kiện ban đầu mơ hình sóng 82 3.5.6 Hiệu chỉnh mô hình sóng 82 3.6 Kết luận chương 85 CHƯƠNG NGHIÊN CỨU DIỄN BIẾN THỦY ĐỘNG LỰC CỬA TIÊN CHÂU………………………………………………… 86 4.1 Khái quát quy luật diễn biến cửa Tiên Châu 86 4.2 Xây dựng tập kịch tính tốn mơ 89 4.2.1 Thiết lập tham số mơ nhóm kịch lũ 90 4.2.2 Thiết lập dựng tham số mô nhóm kịch sóng 95 4.2.3 Thiết lập tham số mô cho kịch bão đổ khu vực nghiên cứu 97 4.2.4 Tổng hợp kịch mô 100 4.3 Phân tích kết mơ 102 4.3.1 Kết mơ nhóm kịch lũ 102 4.3.2 Kết mơ nhóm kịch sóng 108 4.3.3 Kết mô kịch bão 111 4.4 Kết luận chương 115 KẾT LUẬN VÀ KIẾN NGHỊ 118 Kết luận 118 Kiến nghị 120 DANH MỤC CƠNG TRÌNH ĐÃ CƠNG BỐ CỦA TÁC GIẢ 122 TÀI LIỆU THAM KHẢO 123 PHỤ LỤC 127 iv DANH MỤC HÌNH ẢNH Hình 1.1 Bản đồ khu vực nghiên cứu Hình 1.2 Lưu vực sơng Kỳ Lộ Hình 1.3 Vị trí khu vực cửa Tiên Châu Hình 1.4 Vị trí điểm trích xuất sóng T1 12 Hình 1.5 Biến trình độ cao sóng điểm T1 giai đoạn 2011 – 2014 13 Hình 1.6 Hoa sóng điểm T1 giai đoạn 2011 – 2014 .13 Hình 1.7 Mặt bến cảng Tiên Châu khu nước trước bến 20 Hình 1.8 Vị trí số cơng trình bảo vệ bờ sơng Kỳ Lộ tiêu biểu khu vực .24 Hình 1.9 Vị trí kè An Ninh Đơng 26 Hình 1.10 Mặt cắt điển hình kè An Ninh Đơng 26 Hình 1.11 Đê, kè Bình Thạnh, xã Xuân Thọ 2, thị xã Sông Cầu 28 Hình 1.12 Mặt cắt ngang đại diện đê 28 Hình 1.13 Mặt cắt ngang đại diện kè 29 Hình 1.14 Kè sơng Vét xã An Dân, huyện Tuy An, tỉnh Phú Yên .30 Hình 1.15 Cụm cơng trình mỏ hàn, xã An Ninh Tây, huyện Tuy An, tỉnh Phú Yên 31 Hình 1.16 Một số cửa sơng miền Trung điển hình có dạng lưỡi cát cua 35 Hình 2.1 Quy trình thu thập, chiết tách đường bờ từ ảnh viễn thám .41 Hình 2.2 Ảnh quang học trước sau tăng cường chấp lượng ảnh 42 Hình 2.3 Ranh giới khu vực cắt ảnh 45 Hình 2.4 Sơ họa cách xác định chiều rộng cửa Tiên Châu 46 Hình 2.5 Sơ họa cách xác định góc doi cát so với hướng Bắc .47 Hình 2.6 Diễn biến chiều rộng doi cát giai đoạn từ 1988 đến 2020 48 Hình 2.7 Diễn biến góc doi cát thời kì nhiều năm 49 Hình 2.8 Tỉ lệ xuất góc doi cát từ 1988–2020 .49 Hình 2.9 Biểu đồ trận lũ bão có ảnh hưởng tới khu vực nghiên cứu 50 Hình 2.10 Cửa Tiên Châu sau trận lũ lịch sử năm 2009 51 Hình 2.11 Diễn biến chiều rộng cửa Tiên châu điều kiện lũ .52 Hình 2.12 Tương quan chiều rộng cửa lưu lượng đỉnh lũ Qmax 52 Hình 2.13 Diễn biến góc doi cát điều kiện lũ giai đoạn 1995–2020 53 Hình 2.14 Tương quan góc doi cát chiều rộng cửa Tiên Châu .54 Hình 3.1 Dữ liệu địa hình khu vực nghiên cứu .64 Hình 3.2 Vị trí trạm đo mực nước trạm đo sóng, dịng chảy 66 Hình 3.3 Sơ đồ vị trí trạm thủy văn Hà Bằng trạm quan trắc tài ngun nước 67 Hình 3.4 Miền tính mơ hình 68 Hình 3.5 Lưới tính mơ hình thủy lực .69 Hình 3.6 Các biên mơ hình .71 Hình 3.7 Các vị trí trích xuất kết mơ 74 v Hình 3.8 Vị trí trạm đo mực nước sử dụng để hiệu chỉnh kiểm định mơ hình thủy lực 75 Hình 3.9 Hiệu chỉnh mực nước từ 02/10/2019 tới 14/10/209 77 Hình 3.10 Kiểm định mực nước từ 5/1/2020 đến 16/1/2020 79 Hình 3.11 Miền tính mơ hình sóng 79 Hình 3.12 Lưới tính khu vực 80 Hình 3.13 Lưới tính chi tiết lưới tính khu vực 81 Hình 3.14 Các biên mơ hình sóng 81 Hình 3.15 Vị trí trạm đo sóng sử dụng để hiệu chỉnh kiểm mơ hình sóng 82 Hình 3.16 Kết hiệu chỉnh chiều cao sóng 84 Hình 3.17 Kết hiệu chỉnh chu kỳ sóng 84 Hình 3.18 Kết hiệu chỉnh hướng sóng 84 Hình 4.1 Diễn biển chiều rộng cửa Tiên Châu điều kiện lũ 86 Hình 4.2 Cửa Tiên Châu trước bão (ngày 21/8/2009) sau bão Mirinae (ngày 28/1/2010) 88 Hình 4.3 Cửa Tiên Châu trước bão (24/09/2020) sau bão Molave (18/11/2020) 89 Hình 4.4 Đường tần suất lũ trạm An Thạnh 91 Hình 4.5 Đường biến trình lưu lượng lũ trạm An Thạnh sử dụng để mô kịch 91 Hình 4.6 Đường trình mực nước triều tháng 11/2019 95 Hình 4.7 Hoa sóng trung bình nhiều năm 2005–2018 96 Hình 4.8 Tương quan chiều cao – chu kì sóng giai đoạn 2005–2018 96 Hình 4.9 Đường bão Damrey 2017 98 Hình 4.10 Diễn biến chiều cao sóng bão DAMREY 11/2017 [34] 99 Hình 4.11 Diễn biến chu kỳ sóng bão DAMREY 11/2017 [34] 99 Hình 4.12 Diễn biến hướng sóng bão DAMREY 11/2017 [34] 100 Hình 4.13 Kết mơ trường dịng chảy KB1 KB2 khu vực nghiên cứu thời điểm lưu lượng lũ đạt đỉnh 102 Hình 4.14 Diễn biến mực nước, lưu tốc vị trí cửa Tiên Châu KB1 KB2 103 Hình 4.15 Diễn biến lưu tốc dịng chảy vị trí trích xuất kết mô KB1, KB2 103 Hình 4.16 Vị trí trạm trích xuất kết tính tốn 104 Hình 4.17 Đường trình mực nước số vị trí KB1, KB2 105 Hình 4.18 Kết mơ trường vận tốc dịng chảy KB1 KB2 vị trí đoạn sơng co hẹp khúc sông cong 105 Hình 4.19 Kết mơ trường vận tốc dịng chảy KB1 KB2 khu vực cửa Tiên Châu thời điểm lũ lưu lượng lũ đạt đỉnh 106 Hình 4.20 Bồi lắng phân bố bùn khu vực nghiên cứu 107 Hình 4.21 Bề rộng hướng mũi tên cát cửa Tiên Châu trước sau lũ 2016 108 Hình 4.22 Kết mơ trường sóng KB3, KB4 108 vi Hình 4.23 Trường sóng khu vực cửa Tiên Châu 109 Hình 4.24 Trường dịng chảy KB3 KB4 .110 Hình 4.25 Đường bờ khu vực nghiên cứu thời kì khơng có lũ 111 Hình 4.26 Trường sóng KB5 .111 Hình 4.27 So sánh trường sóng KB3 KB5 .112 Hình 4.28 Trường dịng chảy bão KB5 113 Hình 4.29 So sánh trường dòng chảy KB3 KB5 114 Hình 4.30 Đường bờ khu vực nghiên cứu trước sau trận bão Damrey 2017 115 vii DANH MỤC BẢNG Bảng 1.1 Đặc trưng hình thái sơng Kỳ Lộ Bảng 1.2 Đặc trưng dịng chảy trung bình sơng Kỳ Lộ 14 Bảng 1.3 Mực nước lưu lượng dịng chảy năm theo tần suất sơng Kỳ Lộ trạm Hà Bằng 14 Bảng 1.4 Tóm tắt thơng tin số trận lũ lớn lưu vực sông Kỳ Lộ 16 Bảng 1.5 Diện tích, dân số mật độ dân số khu vực nghiên cứu năm 2017 16 Bảng 1.6 Dân số trung bình phân theo nam nữ giai đoạn 2012 – 2016 17 Bảng 1.7 Dân số trung bình phân theo thành thị, nông thôn giai đoạn 2012 – 2016 17 Bảng 1.8 Thông tin cơng trình đập Tam Giang 22 Bảng 1.9 Các thơng tin cơng trình đập Hà Yến 23 Bảng 1.10 Các thơng tin cơng trình đập Đồng Kho 23 Bảng 1.11 Tóm tắt thơng tin số cơng trình bảo vệ bờ sơng Kỳ Lộ khu vực nghiên cứu 24 Bảng 1.12 Bảng thông số kỹ thuật kè An Ninh Đông 27 Bảng 1.13 Các thơng số kỹ thuật chủ yếu tuyến đê, kè Bình Thạnh 29 Bảng 1.14 Các cách phân loại cửa sông giới [8] [9],[10], [11], [12], [13] ,[14] 33 Bảng 2.1 Thống kê loại ảnh, số ảnh độ phân giải, thời gian thu thập 40 Bảng 3.1 Tổng hợp vị trí trạm đo việc xây dựng hiệu chỉnh mơ hình 67 Bảng 3.2 Bộ số điều hòa 70 Bảng 3.3 Các thơng số mơ hình 72 Bảng 3.4 Tọa độ điểm trích xuất kết tính tốn 73 Bảng 3.5 Vị trí mặt cắt trích xuất kết tính toán 73 Bảng 3.6 Đánh giá mức độ phù hợp mơ hình số Nash–Sutcliffe 76 Bảng 3.7 Bảng thông số hiệu chỉnh mơ hình thủy lực 77 Bảng 3.8 Kết kiểm định mơ hình 78 Bảng 3.9 Kết kiểm định mơ hình thủy lực 78 Bảng 3.10 Bảng thông số hiệu chỉnh mơ hình sóng 83 Bảng 3.11 đánh giá sai số độ cao sóng 85 Bảng 4.1 Lưu lượng lũ tháng 11/2016 trạm An Thạnh 92 Bảng 4.2 Thống kê bão đổ vào Phú Yên lân cận 98 Bảng 4.3 Tập kịch tính tốn mơ 101 viii DANH MỤC TỪ VIẾT TẮT BĐ BIAS BTCT Delft-3D ĐTĐL E, N, W, S GIS HĐND KB KHCN Mike 21 FM NOAA RMSE SW SWAN UBND WAVEWATCH III Báo động Sai số trung bình Bê tơng cốt thép Mơ hình thủy lực vận chuyển bùn cát mã nguồn mở Đề tài độc lập Hướng Đông, Bắc, Tây, Nam Hệ thống thông tin đại lý Hội đồng Nhân dân Kịch Khoa học cơng nghệ Mơ hình thủy động lực bùn cát chiều Cơ quan Quản lý Khí Đại dương Quốc gia Mỹ Sai số trung bình quân phương Hướng Tây Nam Mơ hình sóng nước sơng Ủy ban Nhân dân Mơ hình sóng tồn cầu ix DANH MỤC CƠNG TRÌNH ĐÃ CƠNG BỐ CỦA TÁC GIẢ Tùng T T and Khánh M D., “Nghiên cứu quy luật diễn biến doi cát ven bờ khu vực cửa Tiên Châu ảnh vệ tinh Landsat,” Tạp Chí Khoa Học Kỹ Thuật Thủy Lợi Môi Trường, no 71, Art no 71, 2020 122 TÀI LIỆU THAM KHẢO [1] V A Kiệt and L T Hương, “Đặc trưng hình thái sơng ngịi Nam Trung Bộ,” Tập San Quy Hoạch Tài Nguyên Nước Miền Trung, 2015 [2] Sở Khoa học Cơng nghệ tỉnh Phú n, “Báo cáo khí hậu - thủy văn Phú Yên,” Phú Yên, 2015 [3] V A Kiệt, “Tình hình bão lũ tỉnh Phú Yên năm gần đây, xác định nguyên nhân trận lũ lụt lịch sử năm 2009 địa bàn tỉnh Phú Yên,” Hội thảo khoa học BĐKH Phú Yên, Báo cáo hội thảo, Apr 2010 [4] Cục thống kê tỉnh Phú Yên, “Niên giám thống kê tỉnh Phú Yên 2017,” Phú Yên, 2017 [5] D W Pritchard, “What is an estuary: Physical Viewpoint,” 1967, Accessed: May 15, 2021 [Online] Available: https://tamug-ir.tdl.org/handle/1969.3/24383 [6] R W Dalrymple, B A Zaitlin, and R Boyd, “Estuarine facies models; conceptual basis and stratigraphic implications,” J Sediment Res., vol 62, no 6, pp 1130–1146, Nov 1992, doi: 10.1306/D4267A69-2B26-11D7-8648000102C1865D [7] P V Giáp and L P Hậu, Chỉnh trị cửa sông ven biển, vol 215 Hà Nội: Nhà xuất Xây dựng, 1996 [8] N V Giáp, “Nghiên cứu chế bồi, xói đề xuất giải pháp chỉnh trị khu vực Cửa Lấp, tỉnh Bà Rịa Vũng Tàu,” Luận án Tiến sĩ, Đại học Thủy lợi, Hà Nội, 2015 [9] L P Hậu, Động lực học cơng trình cửa sơng Hà Nội: NXB Xây dựng, 2005 [10] Thuỷ V T T., “Giới thiệu phương pháp việc phân loại vùng cửa sông ven biển,” 2014, Accessed: May 16, 2021 [Online] Available: http://tailieuso.tlu.edu.vn/handle/DHTL/8408 [11] J M Coleman and L D Wright, “Modern river deltas: variability of processes and sand bodies,” 1975 123 [12] A T Ippen, Estuary and Coastline Hydronamics McGraw-Hill, New York, 1966 [13] H B Simmons, “View of SALINITY PROBLEMS | Coastal Engineering Proceedings,” 1951 https://journals.tdl.org/icce/index.php/icce/article/view/1473/759 (accessed May 17, 2021) [14] G M Friedman and J E Sanders, Principles of sedimentology Wiley, New York, 1978 [15] L P Hậu, T V An, and L P Hợp, Diễn biến cửa sông vùng đồng Bắc Bộ Nhà xuất Xây dựng, 2002 [16] N T Lam, “Thủy động lực động lực hình thái hệ thống cửa sơng chịu tác động theo mùa, cửa Thuận An Tư Hiền, Thừa Thiên Huế,” Luận án Tiến sĩ, Đại học Công nghệ Delft, Hà Lan [17] P T Hương, N B Qùy, and N L Long, “Ứng dụng mơ hình Mike 21 FM nghiên cứu ảnh hưởng sóng dịng chảy đến cửa sơng Đà Rằng tỉnh Phú n,” Tạp Chí Khoa Học Công Nghệ Hàng Hải, no 27, pp 42–46, Aug 2011 [18] T N Anh, Đ Đ Đức, and Đ Đ Khá, “Xây dựng phương pháp đánh giá tác động biến đổi khí hậu nước biển dâng tới cơng trình hạ tầng kỹ thuật-thử nghiệm cho khu vực ven biển tỉnh Khánh Hịa,” Tạp Chí Khoa Học ĐHQGHN Các Khoa Học Trái Đất Và Môi Trường, no Tập 29, số 4, pp 1–12, 2013 [19] T N Anh and N B Tùng, “Ứng dụng mô Mike mô xâm nhập mặn sông Ninh Cơ điều kiện biến đổi khí hậu nước biển đâng,” Tạp Chí KHCN Thủy Lợi, no 38, pp 21–32, 2020 [20] “Nghiên cứu chế bồi lấp cửa An Dũ giải pháp khắc phục,” Bình Định, Việt Nam, Dự án [21] “Nghiên cứu đánh giá giải pháp công nghệ chống xói lở bờ biển áp dụng đề xuất giải pháp công nghệ phù hợp cho số vùng trọng điểm miền 124 Trung,” Trường Đại học Thủy lợi, Hà Nội, Đề tài khoa học công nghệ cấp Bộ, 2019 2017 [22] “Nghiên cứu ứng dụng giải pháp nuôi bãi nhân tạo cho đoạn bờ biển bị xói lở khu vực miền Trung Việt Nam,” Trường Đại học Thủy lợi, Hà Nội, Đề tài nghiên cứu tiềm Nhà nước KC.08.TN03/11–15, 2013 [23] Luân N T and Sơn N H., “Nghiên cứu biến động vùng cửa sông Cái, Nha Trang qua tư liệu viễn thám (giai đoạn 1999-2013),” Tạp Chí Khoa Học Kỹ Thuật Thủy Lợi Và Môi Trường, no 45, Art no 45, 2014 [24] V A Kiệt and B V Chanh, “Nghiên cứu ứng dụng mơ hình Marine Muskingum dự báo thủy văn lưu vực sông Kỳ Lộ tỉnh Phú Yên,” Tạp Chí Khoa Học Tài Ngun Và Mơi Trường, no 20, pp 66–72, 2018 [25] N V Bình, “Nghiên cứu đề xuất mơ hình quản lý chất lượng xây dựng cơng trình tuyến lũ cứu nạn vượt sơng Kỳ Lộ, huyện Đồng Xuân, tỉnh Phú Yên,” Luận văn thạc sĩ, Đại học Thủy lợi, 2016 [26] T C Danh, “Nghiên cứu, bổ sung đặc điểm khí hậu - thủy văn Phú Yên xây dựng đồ nguy ngập lụt lưu vực sông Kỳ Lộ đoạn từ Phú Mỡ đến hạ lưu,” Phú Yên, Đề tài NCKH cấp tỉnh, 2012 [27] “Nghiên cứu giải pháp chỉnh trị chống sa bồi luồng tàu cho cảng cá khu neo đậu tàu thuyển tỉnh Phú Yên vùng lân cận, áp dụng cho cửa Tiên Châu,” Hà Nội, Đề tài khoa học công nghệ độc lập cấp Quốc gia mã số ĐTĐL.CN.33/18 [28] Tùng T T and Khánh M D., “Nghiên cứu quy luật diễn biến doi cát ven bờ khu vực cửa Tiên Châu ảnh vệ tinh Landsat,” Tạp Chí Khoa Học Kỹ Thuật Thủy Lợi Và Mơi Trường, no 71, Art no 71, 2020 [29] “EarthExplorer.” https://earthexplorer.usgs.gov/ (accessed May 19, 2021) [30] T T Tùng, Đ Đ Đức, M D Khánh, N B Tùng, and L T Hải, “Báo cáo khảo sát thủy hải văn khu vực cửa Tiên Châu lưu vực sông Kỳ Lộ,” Đại học Thủy lợi, Hà Nội, Đề tài KHCN cấp Nhà nước, mã số ĐTĐLCN.33/18, 2019 125 [31] E Boak and I Turner, “Shoreline Definition and Detection: A Review,” J Coast Res J Coast Res., vol 21, pp 688–703, Jul 2005, doi: 10.2112/03-0071.1 [32] Deltares, User Manual Delft2D-FLOW The Netherlands, 2011 [33] Deltares, User Manual Delft2D-Wave The Netherlands, 2011 [34] D C Điển, Đ V Mạnh, N T Việt, N N Thế, and N V Định, “Tính tốn đặc trưng sóng biển Đơng mơ hình SWAN sử dụng trường gió tái phân tích NCEP,” Tuyển Tập Cơng Trình Hội Nghị Khoa Học Cơ Học Thủ Khí Tồn Quốc, 2017 126 PHỤ LỤC Phụ lục 1: Danh sách ảnh viễn thám thu thập TT Loại ảnh Độ phân giải Thời gian chụp STT Loại ảnh Độ phân giải Thời gian chụp Landsat 30m 1/24/2020 3:01 35 Landsat 30m 4/7/2003 2:49 Landsat 30m 12/23/2019 3:01 36 Landsat 30m 11/14/2002 2:49 Landsat 30m 5/13/2019 3:00 37 Landsat 30m 4/4/2002 2:49 Landsat 30m 11/2/2018 3:01 38 Landsat 30m 7/30/2001 2:41 Landsat 30m 5/10/2018 3:00 39 Landsat 30m 3/24/2001 2:41 Landsat 30m 9/12/2017 3:01 40 Landsat 30m 11/8/2000 2:51 Landsat 30m 5/23/2017 3:00 41 Landsat 30m 4/14/2000 2:53 Landsat 30m 9/9/2016 3:01 42 Landsat 30m 9/11/1999 2:37 Landsat 30m 4/18/2016 3:00 43 Landsat 30m 3/19/1999 2:40 10 Landsat 30m 9/7/2015 3:00 44 Landsat 30m 10/26/1998 2:39 11 Landsat 30m 1/26/2015 3:01 45 Landsat 30m 3/16/1998 2:37 12 Landsat 30m 9/4/2014 3:01 46 Landsat 30m 11/24/1997 2:35 13 Landsat 30m 1/7/2014 3:02 47 Landsat 30m 5/16/1997 2:29 14 Landsat 30m 11/4/2013 3:02 48 Landsat 30m 9/2/1996 2:18 15 Landsat 30m 5/28/2013 3:03 49 Landsat 30m 8/17/1996 2:17 16 Landsat 30m 11/9/2012 2:56 50 Landsat 30m 10/18/1995 2:01 17 Landsat 30m 7/4/2012 2:55 51 Landsat 30m 5/27/1995 2:08 18 Landsat 30m 10/22/2011 2:54 52 Landsat 30m 10/15/1994 2:16 19 Landsat 30m 6/8/2011 2:50 53 Landsat 30m 4/6/1994 2:21 20 Landsat 30m 9/25/2010 2:51 54 Landsat 30m 9/26/1993 2:23 21 Landsat 30m 1/28/2010 2:51 55 Landsat 30m 4/3/1993 2:23 127 TT Loại ảnh Độ phân giải Thời gian chụp STT Loại ảnh Độ phân giải Thời gian chụp 22 Landsat 30m 8/21/2009 2:50 56 Landsat 30m 12/28/1992 2:22 23 Landsat 30m 10/5/2008 2:45 57 Landsat 30m 10/7/1991 2:25 24 Landsat 30m 7/17/2008 2:47 58 Landsat 30m 4/30/1991 2:23 25 Landsat 30m 2/8/2008 2:51 59 Landsat 30m 11/5/1990 2:20 26 Landsat 30m 9/1/2007 2:54 60 Landsat 30m 3/26/1990 2:21 27 Landsat 30m 4/10/2007 2:55 61 Landsat 30m 12/4/1989 2:24 28 Landsat 30m 6/10/2006 2:53 62 Landsat 30m 10/1/1989 2:26 29 Landsat 30m 9/6/2006 2:50 63 Landsat 30m 4/8/1989 2:30 30 Landsat 30m 11/14/2005 2:49 64 Landsat 30m 9/12/1988 2:31 31 Landsat 30m 6/23/2005 2:48 65 Landsat 30m 6/8/1988 2:31 32 Landsat 30m 11/11/2004 2:46 66 Landsat 30m 8/21/1975 2:19 33 Landsat 30m 8/23/2004 2:44 67 Landsat 30m 9/9/1973 2:31 34 Landsat 30m 7/28/2003 2:49 68 Landsat 30m 2/23/1973 2:33 Phụ lục 2: Thống kê tham số đặc trưng hình học doi cát cửa Tiên Châu Năm 1988 1989 1990 1991 Thời gian Chiều rộng cửa B (m) Góc mũi tên cát (độ) (so hướng Bắc) 08/06/1988 112 97 12/09/1988 112 97 08/04/1989 395 80 01/10/1989 330 99 04/12/1989 134 55 26/03/1990 123 86 05/11/1990 278 82 03/04/1991 164 95 128 Năm 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 Thời gian Chiều rộng cửa B (m) Góc mũi tên cát (độ) (so hướng Bắc) 07/10/1991 185 95 28/12/1992 250 82 03/04/1993 233 87 26/09/1993 233 87 06/04/1994 220 68 15/10/1994 237 79 27/05/1995 161 82 18/10/1995 161 91 17/08/1996 82 131 02/09/1996 82 131 16/05/1997 82 140 24/11/1997 72 110 16/03/1998 82 115 26/10/1998 124 121 19/03/1999 93 103 11/09/1999 93 103 14/04/2000 113 106 08/11/2000 71 106 21/03/2001 96 97 04/04/2002 98 89 04/11/2002 91 89 07/04/2003 78 124 23/08/2004 93 81 11/11/2004 100 81 129 Thời gian Chiều rộng cửa B (m) Góc mũi tên cát (độ) (so hướng Bắc) 23/06/2005 74 70 14/11/2005 80 105 10/06/2006 110 103 10/04/2007 100 89 01/09/2007 85 127 08/02/2008 181 127 17/07/2008 218 96 05/10/2008 206 160 21/08/2009 191 63 28/01/2010 328 28 25/09/2010 328 28 2011 08/06/2011 303 70 2013 28/05/2013 96,3 125 07/01/2014 103,5 150 15/05/2014 44,2 128 31/05/2014 72 128 04/09/2014 56,2 128 26/01/2015 92,3 148 11/02/2015 92,3 148 02/05/2015 57,5 141 07/09/2015 57,5 138 18/04/2016 59,8 140 04/05/2016 60,7 140 23/07/2016 61,7 140 Năm 2005 2006 2007 2008 2009 2010 2014 2015 2016 130 Năm 2017 2018 2019 2020 Thời gian Chiều rộng cửa B (m) Góc mũi tên cát (độ) (so hướng Bắc) 09/09/2016 45,7 140 23/05/2017 166,8 55 11/08/2017 166,8 55 12/09/2017 153,7 55 07/03/2018 148 41 10/05/2018 165 48 02/11/2018 160 78 13/05/2019 97,5 97 23/12/2019 85,5 128 24/01/2020 53,7 115 Phụ lục 3: Kết mơ trường sóng KB3, KB4,KB5 Hình PL Kết mơ trường vận tốc dịng chảy vị trí đoạn sơng co hẹp khúc sơng cong KB1 131 Hình PL Kết mơ trường vận tốc dòng chảy khu vực cửa Tiên Châu thời điểm lũ lưu lượng lũ đạt đỉnh KB1 Hình PL Kết mơ trường vận tốc dịng chảy vị trí đoạn sơng co hẹp khúc sơng cong KB2 132 Hình PL Kết mơ trường vận tốc dịng chảy khu vực cửa Tiên Châu thời điểm lũ lưu lượng lũ đạt đỉnh KB2 Hình PL Trường sóng miền tính KB3 133 Hình PL Trường sóng khu vực cửa Tiên Châu KB3 Hình PL Trường sóng miền tính KB4 134 Hình PL Trường sóng khu vực cửa Tiên Châu KB4 Hình PL Trường sóng miền tính KB5 135 Hình PL 10 Trường sóng khu vực cửa Tiên Châu KB5 136