Nghiên cứu ảnh hưởng của nước dâng do sóng đến biến động bãi biển và giải pháp bảo vệ bãi biển cửa đại, hội an

BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT TRƯỜNG ĐẠI HỌC THỦY LỢI NGUYỄN NGỌC THẾ NGHIÊN CỨU ẢNH HƯỞNG CỦA NƯỚC DÂNG DO SÓNG ĐẾN BIẾN ĐỘNG BÃI BIỂN VÀ GIẢI PHÁP BẢO VỆ BÃI BIỂN CỬA ĐẠI, HỘI AN LUẬN ÁN TIẾN SĨ KỸ THUẬT HÀ NỘI, NĂM 2021 BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT TRƯỜNG ĐẠI HỌC THỦY LỢI NGUYỄN NGỌC THẾ NGHIÊN CỨU ẢNH HƯỞNG CỦA NƯỚC DÂNG DO SÓNG ĐẾN BIẾN ĐỘNG BÃI BIỂN VÀ GIẢI PHÁP BẢO VỆ BÃI BIỂN CỬA ĐẠI, HỘI AN Ngành: Kỹ thuật xây dựng cơng trình Biển Mã số : 958 02 03 NGƯỜI HƯỚNG DẪN KHOA HỌC: PGS.TS Trần Thanh Tùng GS TS Nguyễn Trung Việt HÀ NỘI, NĂM 2021 LỜI CAM ĐOAN Tác giả xin cam đoan cơng trình nghiên cứu riêng tác giả Các kết nghiên cứu kết luận luận án trung thực, không chép từ nguồn hình thức Việc tham khảo nguồn tài liệu thực trích dẫn ghi nguồn tài liệu tham khảo theo quy định Tác giả luận án Nguyễn Ngọc Thế i LỜI CẢM ƠN Lời tác giả xin bày tỏ lòng cảm ơn sâu sắc tới PGS.TS Trần Thanh Tùng GS.TS Nguyễn Trung Việt tận tình hướng dẫn tác giả suốt trình nghiên cứu để hoàn thành luận án Tác giả xin bày tỏ lòng biết ơn Ban Giám hiệu, phòng Đào tạo, khoa Cơng trình, mơn Kỹ thuật Cơng trình Biển, trường Đại Học Thủy Lợi Trường Cao đẳng Công nghệ - Kinh tế Thủy lợi miền Trung, giúp đỡ tạo điều kiện để tác giả hoàn thành luận án Tác giả xin cảm ơn Ban chủ nhiệm đề tài: “Nghiên cứu q trình xói lở bờ biển Hội An đề xuất giải pháp bảo vệ bờ biển cách bền vững” cho phép tác giả thành viên chính, trực tiếp tham gia nghiên cứu sử dụng số liệu dự án Nhân dịp này, tác giả xin bày tỏ lời cảm ơn sâu sắc tới nhà khoa học, nhóm cộng tác nghiên cứu, thầy đồng nghiệp với tình cảm iiong chân thành động viên, giành nhiều thời gian công sức giúp đỡ đóng góp ý kiến quý báu suốt trình nghiên cứu thực luận án Cuối cùng, tác giả xin bày tỏ lòng biết ơn đến gia đình quan tâm, động viên, khích lệ, ủng hộ tạo điều kiện thuận lợi để tác giả yên tâm thực hoàn thành luận án ii MỤC LỤC MỞ ĐẦU 1 Tính cấp thiết đề tài .1 Mục tiêu nghiên cứu Đối tượng phạm vi nghiên cứu .4 Cách tiếp cận phương pháp nghiên cứu Ý nghĩa khoa học thực tiễn .5 Cấu trúc luận án CHƯƠNG TỔNG QUAN NGHIÊN CỨU VỀ NƯỚC DÂNG DO SĨNG VÀ CƠNG TRÌNH BẢO VỆ BỜ BIỂN 1.1 Khái quát chung nước dâng sóng ảnh hưởng nước dâng sóng 1.1.1 Các thuật ngữ liên quan đến nước dâng sóng .7 1.1.2 Phạm vi ảnh hưởng tác động nước dâng sóng tới bãi biển 1.2 Tổng quan nghiên cứu nước dâng sóng giới 10 1.3 Tổng quan nghiên cứu Việt Nam nước dâng bão, nước dâng sóng cơng trình bảo vệ bờ biển 15 1.3.1 Các nghiên cứu nước dâng, nước dâng sóng 15 1.3.2 Các giải pháp cơng trình bảo vệ bờ biển áp dụng Việt Nam 17 1.4 Các nghiên cứu liên quan đến biển Cửa Đại, Hội An 19 1.5 Đặt vấn đề nghiên cứu luận án 22 1.6 Kết luận Chương 22 CHƯƠNG CƠ SỞ KHOA HỌC XÁC ĐỊNH NƯỚC DÂNG DO SÓNG VÀ BIẾN ĐỘNG BÃI BIỂN 24 2.1 Lựa chọn phương pháp xác định nước dâng sóng luận án 24 2.2 Các số liệu phục vụ nghiên cứu 26 2.2.1 Tài liệu địa hình 26 2.2.2 Tài liệu thủy, hải văn 29 2.2.3 Số liệu bùn cát đáy 35 2.3 Phương pháp phân tích nước dâng sóng từ số liệu đo đạc 36 2.3.1 Cơ sở khoa học phân tích nước dâng sóng từ hình ảnh camera 36 iii 2.3.2 Phân tích ảnh, trích xuất dao động mực nước từ camera cọc tiêu .43 2.3.3 Tính tốn nước dâng sóng từ chuỗi số liệu đo đạc 44 2.4 Phương pháp mơ hình tốn mơ nước dâng sóng .48 2.4.1 Giới thiệu mơ hình sử dụng nghiên cứu luận án .48 2.4.2 Tính tốn lan truyền sóng từ nước sâu vào vùng ven bờ 49 2.4.3 Tính tốn nước dâng sóng biến động bãi biển 54 2.5 Kết luận Chương 58 CHƯƠNG NGHIÊN CỨU ẢNH HƯỞNG CỦA NƯỚC DÂNG DO SÓNG ĐẾN BIẾN ĐỘNG BÃI BIỂN CỬA ĐẠI, HỘI AN 59 3.1 Đặt vấn đề 59 3.2 Phân tích khái quát số quy luật biến động mặt cắt ngang bãi 61 3.2.1 Phân tích biến động mặt cắt ngang bãi biển qua giai đoạn 61 3.2.2 Khái quát số quy luật biến động mặt cắt ngang bãi biển 67 3.3 Nghiên cứu ảnh hưởng nước dâng sóng đến biến động bãi biển 69 3.3.1 Trường hợp tính tốn vị trí mặt cắt tính tốn 69 3.3.2 Kết mơ nước dâng sóng bão ven biển Cửa Đại 71 3.3.3 Nghiên cứu ảnh hưởng nước dâng sóng đến biến động bãi biển 78 3.3.4 Phân tích ảnh hưởng nước dâng sóng đến biến động bãi biển 84 3.4 Kết luận Chương 90 CHƯƠNG NGHIÊN CỨU ĐỀ XUẤT GIẢI PHÁP CƠNG TRÌNH BẢO VỆ BÃI BIỂN CỬA ĐẠI, HỘI AN 93 4.1 Giới thiệu chung 93 4.1.1 Vị trí địa lý, địa hình địa mạo khu vực nghiên cứu 93 4.1.2 Hiện trạng bãi biển Cửa Đại, Hội An 94 4.2 Mục tiêu giải pháp cơng trình bảo vệ bãi biển Cửa Đại, Hội An 98 4.3 Các nghiên cứu đề xuất phương án bố trí cơng trình bảo vệ bãi biển 98 4.4 Phương án bố trí cơng trình ngăn cát, giảm sóng 99 4.4.1 Phương án chung bố trí mặt cơng trình 99 4.4.2 Các phương án bố trí cơng trình 100 4.5 Đánh giá hiệu kỹ thuật phương án bố trí cơng trình 101 4.5.1 Phương pháp đánh giá 101 iv 4.5.2 Độ cao sóng sau bố trí phương án cơng trình bảo vệ 101 4.5.3 Biến động bãi biển theo phương án bố trí cơng trình 103 4.6 Lựa chọn phương án cơng trình 106 4.7 Kết luận Chương .107 KẾT LUẬN VÀ KIẾN NGHỊ 108 I Những kết đạt cửa luận án 108 II Những đóng góp luận án 111 III Những tồn kiến nghị hướng tiếp tục nghiên cứu 111 DANH MỤC CƠNG TRÌNH ĐÃ CƠNG BỐ 112 TÀI LIỆU THAM KHẢO .114 PHẦN PHỤ LỤC 119 PHỤ LỤC 119 PHỤ LỤC 154 v DANG MỤC CÁC HÌNH ẢNH Hình Tình trạng xói lở nghiêm trọng khu vực biển Cửa Đại, Hội An Hình Phạm vi khu vực nghiên cứu Hình 1.1 Sơ đồ mơ tả thuật ngữ nước dâng bão, nước dâng sóng Hình 1.2 Hình thái bãi biển Hình 1.3 Cơ chế phá hủy đụn cát ven bờ Hình 1.4 Hình ảnh minh họa tác động NDDS đến xói lở chân đụn cát bãi Tân Mỹ, Cửa Đại đợt bão tháng 10/2016 10 Hình 1.5 Ứng suất xạ sóng 10 Hình 1.6 Trắc diện mực nước sóng (a) độ cao sóng (b) theo liệu thực nghiệm 13 Hình 1.7 Một số hình ảnh tiêu biểu cơng trình bảo vệ bờ biển Việt Nam 19 Hình 2.1 Sơ đồ mô tả liên kết, bổ trợ phương pháp nghiên cứu 25 Hình 2.2 Tổng hợp địa hình khu vực biển Hội An 26 Hình 2.3 Vị trí đo đạc mặt cắt ngang địa hình bãi biển MCN-01, MCN-02 .27 Hình 2.4 Mặt cắt ngang địa hình đợt khảo sát 28 Hình 2.5 Mặt cắt ngang địa hình thời điểm trước, sau đợt ảnh hưởng bão SARIKA năm 2016 MCN-02, bãi Tân Mỹ, Cửa Đại 28 Hình 2.6 Đo đạc khảo sát địa hình khu vực nghiên cứu sau đợt bão 3/2017 .29 Hình 2.7 Biểu đồ thống kê số lượng bão xuất tháng từ 1979 – 2020 30 Hình 2.8 Biểu đồ thống kê số lượng bão phạm vi kính ảnh hưởng 30 Hình 2.9 Đường đi, đồ sóng tái phân tích đồ trường gió bão DOKSURI di chuyển ngày 14/9/2017 15/9/2017 31 Hình 2.10 Vị trí bố trí trạm đo đạc thủy, hải văn khu vực nghiên cứu 32 Hình 2.11 Diễn biến độ cao sóng 02 lần đo, đợt hai trạm SMS01 (phía Bắc Cửa Đại) trạm SMS02 (phía trước cửa sông Cửa Đại) 33 Hình 2.12 Diễn biến độ cao sóng đo đạc đợt 2, tháng 3/2017 trạm SMS01 33 Hình 2.13 Đo đạc hải văn biển Cửa Đại, Hội An từ ngày 10/3/2017 đến 26/3/2017 33 Hình 2.14 Biến trình tốc độ hướng dịng chảy đợt khảo sát đợt trạm SCR1 34 Hình 2.15 Biến trình tốc độ hướng dòng chảy đợt khảo sát đợt trạm SCR1 34 Hình 2.16 Diễn biễn mực nước từ ngày 12/1/2016 đến 4/3/2017 trạm đo 34 vi Hình 2.17 Vị trí lấy mẫu bùn cát đáy 35 Hình 2.18 Thành phần hạt đưa vào mơ hình 35 Hình 2.19 Sơ đồ tổng hợp màu 37 Hình 2.20 Sơ đồ mặt cắt dọc theo cọc tiêu 37 Hình 2.21 Sơ đồ phân tích sóng theo đường mực nước trung bình 37 Hình 2.22 Thiết lập, bố trí camera hệ thống cọc tiêu bãi tắm KS Agribank .38 Hình 2.23 Bố trí hệ thống cọc tiêu MCN địa hình bãi tắm KS Agribank .39 Hình 2.24 Một số hình ảnh lắp đặt thiết bị quan trắc thu nhận dao động mặt nước biển bão DOKSURI (từ ngày 12 đến 15//9/2017) bãi tắm KS Agribank .39 Hình 2.25 Các vị trí, độ dài quy ước kiểm tra cọc tiêu 40 Hình 2.26 Biến dạng ảnh vị trí cọc tiêu 40 Hình 2.27 Sơ đồ khối chương trình trích xuất tín hiệu từ ảnh camera cọc tiêu .42 Hình 2.28 Định dạng file liệu ảnh thu nhận từ camera 42 Hình 2.29 Trích xuất khung hình đại diện 42 Hình 2.30 Tín hiệu màu R-G-B mặt cắt chương trình đọc xử lý 42 Hình 2.31 Dao động mực nước cọc tiêu (từ C1 đến C6) thời đoạn đo đạc 43 Hình 32 Vị trí trạm đo mương dẫn nước vào giếng đo trạm hải văn Sơn Trà 44 Hình 2.33 So sánh mực nước dâng trạm Sơn Trà trạm Cửa Đại, Hội An 45 Hình 2.34 Tương quan mực nước trạm hải văn Sơn Trà trạm đo mực nước Cửa Đại, Hội An giai đoạn từ 12/2015 đến 3/2017 45 Hình 2.35 Đồ thị MN dâng trạm hải văn Sơn Trà (12/9/2016 đến 16/9/2016) 46 Hình 2.36 Phân bố độ cao sóng quan trắc nước dâng sóng theo mặt cắt ngang 47 Hình 2.37 Sơ đồ khối tính tốn nước dâng sóng biến động bãi biển 48 Hình 2.38 Miền tính tốn vùng 1, vùng 49 Hình 2.39 Miền tính tốn sóng chi tiết 49 Hình 2.40 Lưới tính tốn vùng 1, vùng địa hình đáy biển 50 Hình 2.41 Độ cao sóng biên nước sâu, biển Cửa Đại (10/2016) 50 Hình 2.42 Hiệu chỉnh độ cao sóng trạm SMS01 51 Hình 2.43 Hiệu chỉnh độ cao sóng trạm SMS02 51 Hình 2.44 Hiệu chỉnh chu kỳ sóng trạm SMS01 51 Hình 2.45 Hiệu chỉnh chu kỳ sóng trạm SMS02 51 vii Hình 2.46 Hiệu chỉnh hướng sóng trạm SMS01 51 Hình 2.47 Hiệu chỉnh hướng sóng trạm SMS02 51 Hình 2.48 Kiểm định độ cao sóng Trạm SMS01 (3/2017) 53 Hình 2.49 Kiểm định độ cao sóng Trạm SMS02 (3/2017) 53 Hình 2.50 Kiểm định chu kỳ sóng Trạm SMS01 (3/2017) 53 Hình 2.51 Kiểm định chu kỳ sóng Trạm SMS02 (3/2017) 53 Hình 2.52 Kiểm định hướng sóng trạm SMS01 (3/2017) 53 Hình 2.53 Kiểm định hướng sóng trạm SMS02 (3/2017) 53 Hình 2.54 Trường sóng 16h ngày 14/9/2017 54 Hình 2.55 So sánh độ cao sóng đo đạc tính tốn mặt cắt thuộc bãi biển KS AGRIBANK, Cửa Đại, Hội An 57 Hình 2.56 So sánh nước dâng sóng đo đạc tính tốn mặt cắt thuộc bãi biển KS AGRIBANK, Cửa Đại, Hội An 57 Hình 2.57 So sánh biến động địa hình mặt cắt ngang đo đạc tính tốn mặt cắt thuộc bãi biển KS AGRIBANK, Cửa Đại, Hội An 57 Hình 3.1 Bãi biển KS Agribank 10/2016 59 Hình 3.2 Bãi biển KS Agribank 59 Hình 3.3 Sơ đồ phân tích đánh giá ảnh hưởng NDDS tới biến động bãi biển .60 Hình 3.4 Vị trí trích xuất sóng khu vực nước sâu khu vực biển Cửa Đại, Hội An 61 Hình 3.5 Tương quan nước dâng độ cao sóng nước sâu thời kỳ trước mùa gió Tây Nam sau thời kỳ mùa gió Đơng Bắc 2016 62 Hình 3.6 Diễn biễn địa hình mặt cắt ngang MCN-01 từ 23/3/016 đến 17/8/2016 63 Hình 3.7 Hướng trường gió bão ngày 12/9/2016, 17/10/2016 [49] 63 Hình 3.8 Biến động địa hình mặt cắt ngang MCN-01 từ 17/8/2016 đến 26/10/2016 64 Hình 3.9 Biến động bãi biển MCN-02 đợt ảnh hưởng bão SARIKA 65 Hình 3.10 Tương quan mực nước dâng độ cao sóng bão số 7) năm 2016 .65 Hình 3.11 Diễn biến mặt cắt ngang địa hình bãi biển giai đoạn 66 Hình 3.12 Hình ảnh xói lở khu vực biển từ Palm Garden đến Agribank 67 Hình 3.13 Vị trí điểm tính tốn vùng ven bờ biển phía bắc Cửa Đại, Hội An 70 Hình 3.14 Phân bố trường sóng theo tần suất bão khu vực biển Cửa Đại, Hội An .71 Hình 3.15 Biểu đồ so sánh chiều cao sóng điểm biên trường hợp 72 viii 2.2 Cơ sở lý thuyết mơ hình thủy lực chiều XBEACH [56] 2.2.1 Giới thiệu tóm tắt mơ hình mơ hình XBEACH XBeach mơ hình số mã nguồn mở phát triển ban đầu để mô trình thủy động lực học hình thái động lực học tác động lên bờ biển cát với kích thước miền hàng km theo quy mơ thời gian bão Mơ hình bao gồm q trình thủy động lực học biến đổi sóng ngắn (khúc xạ, vỡ), biến đổi sóng dài (tạo ra, lan truyền tiêu tán), thiết lập sóng gây dịng chảy khơng ổn định, ngập nước ngập tràn Các trình hình thái học bao gồm tải trọng đáy vận chuyển trầm tích lơ lửng, sạt lở mặt đụn cát, tác động phá vỡ tầng đáy XBeach có hai chế độ: chế độ thủy tĩnh chế độ không thủy tĩnh Trong chế độ thủy tĩnh, biến đổi biên độ sóng ngắn giải tách biệt với sóng dài, dịng chảy thay đổi hình thái Chế độ khơng thủy tĩnh giải tất trình bao gồm chuyển động sóng ngắn, với nhu cầu tính tốn nhiều 2.2.2 Chế độ khơng thủy tĩnh Lưu lượng trung bình theo độ sâu sóng dịng chảy tính phương trình nước nơng phi tuyến tính, bao gồm áp suất khơng thủy tĩnh Phương trình cân hoạt động sóng ngắn Sóng cưỡng phương trình động lượng vùng nước nơng nhận từ phiên phụ thuộc vào thời gian phương trình cân tác dụng sóng Tương tự mơ hình HISWA [63] Đại học Delft, phân bố theo hướng mật độ hoạt động tính đến, phổ tần số biểu thị tần số, biểu thị tốt tham số phổ fm-1,0 Phương trình cân hoạt động sóng đưa ra: (2 31) Tác động sóng A tính sau: (2 32) đó:  biểu thị góc tới so với trục x, Sw biểu ô định hướng số thị mật độ lượng sóng tần số sóng nội Tần nội từ quan hệ phân tán tuyến tính σ = √gktanhkh vận tốc nhóm cg thu 161 Tốc độ truyền sóng khơng gian x, y theo phương cho bởi: Cx (x,y,t,) =cg cos() (2 33) Cy (x,y,t,) =cg sin() (2 34) (2 35) c0(x, y, t, θ) = đó: h đại diện cho độ sâu nước địa phương k số sóng Tần số sóng nội  xác định mà khơng có tương tác dịng sóng với tần số xuyên tâm tuyệt đối  Tương tác Sóng – dịng chảy Trong XBeach, điều tính đến cách hiệu chỉnh số sóng k với việc sử dụng phương trình Eikonal, điều có tác động đến nhóm tốc độ truyền sóng (x, y khơng gian hướng) Số sóng ngang sóng dọc bờ, k x ky, xác định theo phương trình (2.36; 2.37) Trong cơng thức này, số phụ đề cập đến hướng thành phần vectơ sóng kx = kxn−1 + kx n−1 ky = ky đó: số phụ n-1 đề cập đến số sóng bước thời gian trước đó, k x ky hiệu chỉnh số sóng kx ky số sóng hiệu chỉnh có tính đến diện dịng chảy Các số hạng hiệu chỉnh xác định phương trình thứ hai, phương trình Eikonal: ∂k ∂k x + ∂t y + ∂t Số sóng sau cho: ∂ ∂ w =0 ∂x w (2 38) =0 ∂y k = √k2x + k2y (2 39) Tần số bán kính tuyệt đối  tính bằng: ω = σ + kxuL + kyvL + (2 40) (2 41) L L đó: u v vận tốc Lagrangian trung bình qua bờ độ sâu dọc bờ Tốc độ truyền sóng (theo phương x y) cho bởi: 162 cx(x, y, t, θ) = cg cos(θ) + u L cy(x, y, t, θ) = cg sin(θ) + vL Tốc độ lan truyền không gian có hướng , nơi tính đến khúc xạ đáy (số hạng đầu) khúc xạ dòng điện (hai số hạng cuối), nhận từ: (2 43) cθ(x, y, t, θ) = (2 44) + sinθ (sinθ Sự tiêu tán lượng sóng (2 42) Năm cơng thức phá sóng khác thực XBeach Các cơng thức chọn cách sử dụng ngắt từ khóa (Bảng PL2.1) Bảng PL 2.1 Các cơng thức phá sóng khác thực Công Roe Roe Daly Baldo Janssen Tổng tiêu tán lượng sóng, tức tích hợp theo hướng, sóng phá vỡ mơ hình hóa theo Roelvink (1993) [64] Trong cơng thức tiêu tán sóng vỡ, ý tưởng tính tốn tiêu tán với phần nhỏ sóng đứt Q b nhân với tiêu tán kiện vỡ Trong công thức này,α áp dụng làm hệ số tản sóng O, T rep chu kỳ sóng đại diện E w lượng sóng Phần phá sóng xác định với chiều cao sóng bình phương gốc H mrs chiều cao sóng lớn Hmax Chiều cao sóng tối đa tính tỷ số độ sâu mực nước (h) cộng với phần nhỏ chiều cao sóng  Hrms cách sử dụng số ngắt  Trong cơng thức cho Hrmss  đại diện cho mật độ nước g số hấp dẫn Tổng lượng sóng E w tính cách tích phân thùng hướng sóng ̅η̅ Dw = Qb = − exp (− ( Trong biến thể (2.45), người ta sử dụng cơng thức phá sóng khác, trình bày (2.46) Công thức khác so với công thức Roelvink (1993) chọn cách sử dụng từ khóa: roelvink2 Sự khác biệt với cơng thức ban đầu tiêu tán sóng với roelvink2 tỷ lệ với H /h thay H ̅η̅ Dw = Ngồi ra, cơng thức [65] cho biết sóng hồn tồn vỡ độ cao sóng vượt ngưỡng () ngừng đứt độ cao sóng giảm xuống ngưỡng khác (2) Công thức chọn break = roelvink_daly ngưỡng thứ hai, (2), đặt cách sử dụng từ khóa (2) { Q b Trong trường hợp sóng tĩnh [66] áp dụng, trình bày (2.48, 2.48a) Trong cơng thức phá vỡ này, sóng vỡ phân đoạn Q b chiều cao sóng vỡ H b tính tốn khác so với công thức phá vỡ sử dụng cho tình khơng đứng n, α áp dụng làm hệ số tiêu tán sóng, f rep biểu thị tần số nội đại diện y hệ số hiệu chuẩn ̅η̅ Dw = Qb = exp [− ( Cuối cùng, sử dụng cơng thức [JB07] để phá sóng sóng tĩnh Cơng thức sửa đổi công thức Baldock R= H b H D w rms 164 (2 49) (2.49a) (2.49b) Trong trường hợp cố định đứng yên, tổng tiêu tán sóng phân bố tỷ lệ thuận hướng sóng với cơng thức (2.50) Dw(x, y, t, θ) = Ứng suất xạ Với phân bố khơng gian tác động sóng, ứng suất xạ đánh giá cách sử dụng lý thuyết sóng tuyến tính mơ tả bởi: Sxx,r(x, y, t) = fcos2θSrdθ Sxy,r(x, y, t) = Sỹ,r(x, y, t) = fsinθcosθSrdθ Syy,r(x, y, t) = fsin2θSrdθ Sự nhiễu loạn sóng Sự rối loạn sóng gây bề mặt nước phải chuyển phía đáy để ảnh hưởng đến khuấy lên trầm tích [RS89] sử dụng mơ hình phân rã theo cấp số nhân với độ dài trộn tỷ lệ với Hrms để ước tính thời gian lượng nhiễu loạn trung bình lớp đáy từ nhiễu động mặt nước: k kb = (2 54) −1 h exp( ) H rms kb phương sai nhiễu loạn đáy k phương sai nhiễu loạn trung bình theo thời gian mặt nước kb = exp ( Phương trình nước nơng Đối với sóng tần số thấp dịng chảy trung bình, chúng tơi sử dụng phương trình nước nơng Để tính đến thơng lượng khối lượng sóng gây dịng chảy (trở lại) tiếp theo, chúng đúc thành công thức trung bình theo chiều sâu Lagrangian trung bình (GLM) ([AM78], [WRG00]) Trong khn khổ vậy, phương trình động lượng liên tục xây dựng theo vận tốc Lagrangian định nghĩa khoảng cách mà hạt nước di chuyển chu kỳ sóng, chia cho chu kỳ Vận tốc liên quan đến vận tốc Eulerian (vận tốc trung bình sóng ngắn quan sát điểm cố định) cách: 165 L E u =u +u S S L E vàv = v + v S S u v đại diện cho Stokes trôi theo hướng x- y tương ứng ([Phi77]) Sự trơi dạt Strokes tính tốn lượng E w hướng sóng ngắn thay đổi nhóm sóng thu từ cân hoạt động sóng (2 57) uS = Các phương trình động lượng GLM kết đưa bởi: (2 58) ∂uL = uL (2 59) ∂t ∂v L =u (2 60) L ∂t Độ nhớt ngang Theo mặc định, độ nhớt ngang (vh) tính cách sử dụng mơ hình [Sma63] để tính đến trao đổi động lượng ngang quy mô không gian nhỏ kích thước lưới tính tốn, cho là: v =c h 2 1/2 √( S cS số Smagorinsky đặt 0,1 tất mơ mơ hình Cũng sử dụng giá trị người dùng xác định cho độ nhớt ngang cách tắt mơ hình Smagorinsky định giá trị trực tiếp Phương trình khuếch tán Nồng độ trầm tích cột nước lập mơ hình cách sử dụng sơ đồ khuếch tán-đối lưu trung bình theo độ sâu với thuật ngữ chìm nguồn dựa nồng độ trầm tích cân ([GV83]): ∂hC ∂t ∂hCuE + ∂x In đại diện cho nồng độ trầm tích trung bình theo độ sâu thay đổi theo thang thời gian nhóm sóng hệ số khuếch tán trầm tích Sự theo trầm tích biểu thị thời gian thích ứng, cho phép gần đơn giản dựa độ sâu nước địa phương vận tốc rơi trầm tích Một giá trị nhỏ tương ứng với phản 166 ứng trầm tích gần tức thời Hệ số hệ số hiệu chỉnh hiệu chuẩn có tính đến thực tế xác định liệu trung bình theo độ sâu Ts = max(fTs ws , Ts, min) Vận chuyển trầm tích Trong phiên XBeach, có sẵn hai cơng thức vận chuyển trầm tích Cơng thức hai cơng thức trình bày phần sau Đối với hai phương pháp, tổng nồng độ cặn cân tính theo phương trình (2.64) Trong phương trình này, giá trị nhỏ nồng độ bùn cát cân (cho tải trọng đáy tải trọng lơ lửng) so với nồng độ bùn cát tối đa cho phép Ceq = max(min(Ceq,b, Cơng thức vận chuyển trầm tích có phương trình Soulsby-Van Rijn Nồng độ cặn cân tính tốn theo: (2 65) C (2 66) ) 1.2 C (2 67) Hệ số tải trọng đáy tải trọng treo tính với: Asb = 0.005h ( Trong đường kính cặn khơng thứ ngun (D*) tính theo cơng thức sau Là độ nhớt động học dựa biểu thức Van Rijn hàm nhiệt độ nước XBeach giả định nhiệt độ không đổi 20 độ C, điều dẫn đến độ nhớt động học không đổi D∗=( Nồng độ cặn cân tính tốn theo: C eq,b = 167 A C = eq,s Công thức vận chuyển trầm tích thứ ba có phương trình Van Rijn (1993 phân biệt vận chuyển trầm tích độ cao tham chiếu trầm tích coi vận chuyển tải trọng đáy độ cao tham chiếu coi tải trọng lơ lửng Vận chuyển đáy nằm tính toán với: Sb=0.006swsD50M0.5Me0.7 Đối với tải trọng lơ lửng, trước tiên, nồng độ chuẩn tính theo: c a D50Ta1.5 = 0.015ρ s αD0∗.3 Hiệu ứng độ dốc bãi Ảnh hưởng mái dốc lên thủy động lực cục không xem xét XBeach Hai biểu thức thực để thay đổi cường độ vận chuyển bùn cát Phương thức phương thức mặc định XBeach: q x,slope =q y,slope =q q x y − αhC√ − αhC√ Phương pháp áp dụng cho vận chuyển tổng lượng bùn cát cho vận chuyển tải trọng đáy Phương pháp thứ hai dựa công thức kỹ thuật [Sou97]: qslope = q(1 − α Ngoài ra, phương pháp áp dụng cho vận chuyển toàn vận chuyển theo đáy Để thay đổi hướng vận chuyển tải trọng đáy, sử dụng biểu thức: tan(αψ,new) = f(θ) = 9( qb.x = |qb|cos(αψ,new) qb.y = |qb|sin(αψ,new) 168 Cuối cùng, điều chỉnh việc bắt đầu tiêu chí chuyển động cho vận chuyển toàn vận chuyển tải trọng đáy thơng qua ([Sou97]): adjused θ =θ cr cr Xói lở đụn cát Để giải thích cho sụt giảm vật liệu cát từ mặt đụn cát đến bờ biển q trình xói mịn cồn cát bão gây (từ khóa: tuyết lở) giới thiệu để cập nhật tiến hóa đáy Avalanching giới thiệu thông qua việc sử dụng độ dốc quan trọng cho khu vực khô ướt Người ta cho khu vực ngập nước dễ bị sụt trượt nhiều đó, hai mái dốc tới hạn riêng biệt cho điểm khô ướt sử dụng Giá trị mặc định 0,3 Khi vượt độ dốc tới hạn này, vật liệu trao đổi ô liền kề với lượng cần thiết để đưa mái dốc trở lại độ dốc tới hạn ∂ zb | > mcr | ∂x Để ngăn chặn việc tạo sóng xung kích lớn thay đổi đột ngột mức đáy, việc cập nhật đáy sạt lở giới hạn tốc độ tối đa Công thức (2.81, 2.82) cho thấy thay đổi mức đáy kết bước thời gian ∆ = ((| zb ∆zb= max (− (| 169 ... định nước dâng sóng biến động bãi biển Chương 3: Nghiên cứu ảnh hưởng nước dâng sóng đến biến động bãi biển Cửa Đại, Hội An Chương 4: Đề xuất giải pháp cơng trình bảo vệ bãi biển Cửa Đại, Hội An. .. 3.3.2 Kết mô nước dâng sóng bão ven biển Cửa Đại 71 3.3.3 Nghiên cứu ảnh hưởng nước dâng sóng đến biến động bãi biển 78 3.3.4 Phân tích ảnh hưởng nước dâng sóng đến biến động bãi biển 84 3.4... DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT TRƯỜNG ĐẠI HỌC THỦY LỢI NGUYỄN NGỌC THẾ NGHIÊN CỨU ẢNH HƯỞNG CỦA NƯỚC DÂNG DO SÓNG ĐẾN BIẾN ĐỘNG BÃI BIỂN VÀ GIẢI PHÁP BẢO VỆ BÃI BIỂN CỬA ĐẠI, HỘI AN Ngành: