Nghiên cứu ảnh hưởng của sóng phản xạ đến dòng phản hồi và xói chân đê biển mái nghiêng khu vực Bắc Bộ

Nghiên cứu ảnh hưởng của sóng phản xạ đến dòng phản hồi và xói chân đê biển mái nghiêng khu vực Bắc Bộ.Nghiên cứu ảnh hưởng của sóng phản xạ đến dòng phản hồi và xói chân đê biển mái nghiêng khu vực Bắc Bộ.Nghiên cứu ảnh hưởng của sóng phản xạ đến dòng phản hồi và xói chân đê biển mái nghiêng khu vực Bắc Bộ.Nghiên cứu ảnh hưởng của sóng phản xạ đến dòng phản hồi và xói chân đê biển mái nghiêng khu vực Bắc Bộ.Nghiên cứu ảnh hưởng của sóng phản xạ đến dòng phản hồi và xói chân đê biển mái nghiêng khu vực Bắc Bộ.Nghiên cứu ảnh hưởng của sóng phản xạ đến dòng phản hồi và xói chân đê biển mái nghiêng khu vực Bắc Bộ.Nghiên cứu ảnh hưởng của sóng phản xạ đến dòng phản hồi và xói chân đê biển mái nghiêng khu vực Bắc Bộ.Nghiên cứu ảnh hưởng của sóng phản xạ đến dòng phản hồi và xói chân đê biển mái nghiêng khu vực Bắc Bộ.Nghiên cứu ảnh hưởng của sóng phản xạ đến dòng phản hồi và xói chân đê biển mái nghiêng khu vực Bắc Bộ.

HÀ NỘI, NĂM 2022

TRƯỜNG ĐẠI HỌC THỦY LỢI

NGUYỄN THỊ PHƯƠNG THẢO

NGHIÊN CỨU ẢNH HƯỞNG CỦA SÓNG PHẢN XẠ TỚI DÒNG PHẢN HỒI VÀ XÓI CHÂN ĐÊ BIỂN MÁI

NGHIÊNG KHU VỰC BẮC BỘ

LUẬN ÁN TIẾN SĨ KỸ THUẬT

TRƯỜNG ĐẠI HỌC THỦY LỢI

LỜI CAM ĐOAN

Tác giả xin cam đoan đây là công trình nghiên cứu của bản thân tác giả Các kết quảnghiên cứu và các kết luận trong luận án là trung thực, không sao chép từ bất kỳ mộtnguồn nào và dưới bất kỳ hình thức nào Việc tham khảo các nguồn tài liệu (nếu có) đãđược thực hiện trích dẫn và ghi nguồn tài liệu tham khảo đúng quy định

Tác giả luận án

Nguyễn Thị Phương Thảo

LỜI CÁM ƠN

Tác giả xin trân trọng biết ơn sâu sắc tới thầy GS.TS Thiều Quang Tuấn đã tận tìnhhướng dẫn, chỉ bảo, động viên trong suốt thời gian thực hiện Luận án Tác giả xinđược trân trọng cảm ơn Trường Đại học Thủy lợi, Phòng Thí nghiệm thủy lực tổnghợp, Phòng Đào tạo, Khoa, Bộ môn về những hỗ trợ, tạo điều kiện làm việc tốt nhấtcho tác giả trong quá trình nghiên cứu

Tác giả cũng xin trân trọng cảm ơn các thầy cô, đồng nghiệp, bạn bè và gia đình về sựgiúp đỡ quý giá và luôn đồng hành cùng tác giả trong quá trình thực hiện và hoànthành Luận án của mình

MỤC LỤC

DANH MỤC CÁC HÌNH, ẢNH vi

DANH MỤC BẢNG BIỂU ix

DANH MỤC CÁC TỪ VIẾT TẮT VÀ GIẢI THÍCH THUẬT NGỮ ………… x

MỞ ĐẦU 1

1 Tính cấp thiết của luận án 1

2 Mục tiêu nghiên cứu 4

3 Đối tượng và phạm vi nghiên cứu 4

4 Nội dung nghiên cứu 4

5 Cách tiếp cận và phương pháp nghiên cứu 5

6 Ý nghĩa khoa học và thực tiễn của luận án 6

7 Cấu trúc luận án 7

CHƯƠNG 1 TỔNG QUAN NGHIÊN CỨU DÒNG PHẢN HỒI DO SÓNG VÀ XÓI CHÂN ĐÊ BIỂN 8

1.1 Giới thiệu chung 8

1.1.1 Dòng phản hồi và các quá trình vận chuyển bùn cát ngang bờ 8

1.1.2 Xói chân đê biển 11

1.2 Tổng quan nghiên cứu về dòng phản hồi 12

1.3 Tổng quan mô hình vận chuyển bùn cát và xói lở ngang bờ 17

1.3.1 Tổng quan mô hình sóng 18

1.3.2 Mô hình vận chuyển bùn cát và hình thái 21

1.4 Tổng quan nghiên cứu về sóng phản xạ 25

1.5 Tổng quan các nghiên cứu về dòng phản hồi và xói chân đê biển trong bão ở

Việt Nam 28

1.6 Kết luận chương 1 30

1.6.1 Những vấn đề còn tồn tại 30

1.6.2 Định hướng nghiên cứu của luận án 31

CHƯƠNG 2 CƠ SỞ KHOA HỌC MÔ HÌNH HOÁ DÒNG PHẢN HỒI DO SÓNGVÀ XÓI CHÂN ĐÊ BIỂN 33

2.1 Giới thiệu chung 33

2.2 Ảnh hưởng của phản xạ sóng do công trình đến dòng phản hồi và vận chuyển bùn cát 33

2.2.1 Cơ sở xác định hệ số phản xạ và biến đổi chiều cao sóng trước chân công

trình ……… 34

2.2.2 Ảnh hưởng của sóng phản xạ đến dòng phản hồi 37

2.2.3 Ảnh hưởng của sóng phản xạ đến nồng độ bùn cát 40

2.3 Xây dựng mô hình vật lý máng sóng nghiên cứu dòng phản hồi và xói chân đê biển ……… 42

2.3.1 Mục tiêu và điều kiện thực hiện thí nghiệm 42

2.3.2 Lựa chọn tiêu chuẩn tương tự và tỉ lệ mô hình 43

2.3.3 Thiết kế thí nghiệm và bố trí thiết bị đo đạc 46

2.3.4 Kịch bản thí nghiệm 49

2.3.5 Thực hiện mô hình 51

2.4 Phát triển cập nhật mô hình toán mô phỏng dòng phản hồi và xói bồi chân đê biển ……… 56

2.4.1 Giới thiệu chung 56

2.4.2 Những vấn đề cập nhật trong mô hình Wadibe-TC 56

2.5 Kết luận chương 2 57

CHƯƠNG 3 KẾT QUẢ NGHIÊN CỨU DÒNG PHẢN HỒI VÀ XÓI CHÂN ĐÊ BIỂN MÁI NGHIÊNG 58

3.1 Giới thiệu chung 58

3.2 Phân tích các kết quả nghiên cứu trên mô hình vật lý 58

3.2.1 Nghiên cứu cấu trúc dòng phản hồi 58

3.2.2 Phân tích đánh giá ảnh hưởng của sóng phản xạ đến dòng phản hồi 63

3.2.3 Phân tích kết quả quan trắc xói chân đê trên mô hình lòng động 66

3.2.4 Nhận xét kết quả nghiên cứu trên mô hình vật lý máng sóng 72

3.3 Kết quả phát triển cập nhật mô hình toán mô phỏng cấu trúc dòng phản hồi và xói bồi chân đê biển 74

3.3.1 Phát triển cập nhật mô hình toán Wadibe-TC 74

3.3.2 Kết quả hiệu chỉnh và kiểm định module sóng 74

3.3.3 Kết quả hiệu chỉnh và kiểm định module dòng chảy 78

3.3.4 Kết quả kiểm định module vận chuyển bùn cát và xói chân đê biển 91

3.4 Kết luận chương 3 98

CHƯƠNG 4 ỨNG DỤNG KẾT QUẢ NGHIÊN CỨU TÍNH TOÁN XÓI CHÂN ĐÊ BIỂN NAM ĐỊNH 100

4.1 Giới thiệu chung 100

4.2 Tổng quan khu vực ven biển Nam Định 100

4.2.1 Điều kiện biên tự nhiên khu vực Nam Định 100

4.2.2 Đặc điểm hệ thống đê biển Nam Định 101

4.3 Kiểm định mô hình toán cho vùng ven biển Nam Định 104

4.4 Mô phỏng một số kịch bản xói chân đê biển Nam Định 105

4.4.1 Xây dựng kịch bản 106

4.4.2 Kết quả mô phỏng xói chân đê 107

4.5 Một số đề xuất giảm xói chân đê trong bão 112

4.6 Kết luận chương 4 113

KẾT LUẬN VÀ KIẾN NGHỊ 114

1 Những kết quả đạt được 114

2 Những đóng góp mới của luận án 115

3 Tồn tại và hướng phát triển 116

4 Kiến nghị 116

DANH MỤC CÔNG TRÌNH ĐÃ CÔNG BỐ 118

TÀI LIỆU THAM KHẢO 119

PHỤ LỤC 1

DANH MỤC CÁC HÌNH, ẢNH

Hình 1 Sự cố của đoạn đê Nam Định và Đồ Sơn do bão Doksuri 9.2017 (nguồn: tác

giả,NamDinh TV) 2

Hình 2 Xói lở trước chân đê và sạt lở mái đê, mái và chân gốc kè mỏ hàn sau trận bão Damrey 7/2005 tại Nam Định [3] [4] 3

Hình 3 Cấu trúc của Luận án 7

Hình 1-1 Ảnh hưởng của đê tường đứng đến bãi [5] 8

Hình 1-2 Sơ họa phân bố dòng chảy và nồng độ bùn cát vùng sóng vỡ [18] 11

Hình 1-3 Phân bố hệ số xáo trộn bùn cát [1] 23

Hình 1-4 Bảo toàn thể tích bùn cát [77] 24

Hình 1-5 Nội dung nghiên cứu luận án 32

Hình 2-1Mật độ phương sai của sóng, dòng chảy và hiệp phương sai của sóng và dòng chảy 36

Hình 2-2 Phân bố ứng suất tiếp theo độ sâu [108] 38

Hình 2-3 Phân bố ứng suất trung bình và vận tốc [108] 38

Hình 2-4 Các thành phần trong cấu trúc dòng chảy [19] 39

Hình 2-5 Nguyên lý tính toán vận chuyển bùn cát lơ lửng [1] 40

Hình 2-6 Hệ số xáo trộn [19] 41

Hình 2-7 Mô hình thí nghiệm lòng cứng 47

Hình 2-8 Mô hình lòng động với mái đê m=4 và độ dốc bãi 1/100 49

Hình 2-9 Thực hiện thí nghiệm mô hình kết cấu đê cao không tràn 52

Hình 2-10 Thực hiện mô hình kết cấu đê thấp không thường đỉnh 52

Hình 2-11 Thực hiện thí nghiệm mô hình kết cấu đê thấp có tường đỉnh 52

Hình 2-12 Hiển thị kết quả đo sóng và vận tốc dòng chảy tại một điểm 53

Hình 2-13 Mô hình thực hiện trong máng sóng với các kịch bản kết cấu đê, độ đôc bãi 1/100 54

Hình 2-14 Mô hình thực hiện trong máng sóng với các kịch bản kết cấu đê, độ dốc bãi 1/40 55

Hình 2-15 Địa hình trước và sau khi chạy sóng – độ dốc bãi 1/100 55

Hình 2-16 Địa hình trước và sau khi chạy sóng – độ dốc bãi 1/40 55 Hình 3-1 Phân bố tốc độ dòng phản hồi của các trường hợp thí nghiệm D65H15T1959 Hình 3-2 Phân bố tốc độ dòng phản hồi của các trường hợp thí nghiệm D65H17T1659 Hình 3-3 Phân bố tốc độ dòng phản hồi của các trường hợp thí nghiệm D70H15T1959

Hình 3-4 Phân bố tốc độ dòng phản hồi của các trường hợp thí nghiệm D70H19T165

60

Hình 3-5 Ảnh hưởng của công trình đến dòng chảy trung bình- các thí nghiệm D65H15T19 61

Hình 3-6 Ảnh hưởng của công trình đến dòng chảy trung bình- các thí nghiệm D65H17T16 62

Hình 3-7 Ảnh hưởng của công trình đến dòng chảy trung bình- các thí nghiệm D70H15T19 62

Hình 3-8 Ảnh hưởng của công trình đến dòng chảy trung bình- các thí nghiệm D70H19T165 62

Hình 3-9 Phân bố hệ số phản xạ K r theo khoảng cách tương đối x/L ở lân cận phía trướccông trình 63

Hình 3-10 Tỷ số giữa chiều cao sóng tổng và chiều cao sóng tới H m0,x/H m0i,x thay đổi theo khoảng cách tương đối x/L ở lân cận phía trước công trình 63

Hình 3-11 Phân bố hệ số phản xạ K r cục bộ theo khoảng cách tương đối x/L ở lân cận phía trước công trình 64

Hình 3-12 Kết quả phân tích hồi quy hàm phân bố F x (x/L) .65

Hình 3-13 So sánh giữa kết quả tính toán và đo đạc hệ số phản xạ sóng K r,x 66

Hình 3-14 So sánh độ sâu lớn nhất tương đối giữa đê thấp có và không có tường đỉnh .69

Hình 3-15 Ảnh hưởng của hệ số phản xạ đến độ sâu hố xói lớn nhất tương đối 70

Hình 3-16 Ảnh hưởng của lưu lượng sóng tràn đến độ sâu hố xói 70

Hình 3-17 Ảnh hưởng của độ sâu nước tương đối đến độ sâu hố xói 71

Hình 3-18 Tương quan chỉ số sóng vỡ với độ sâu hố xói lớn nhất 72

Hình 3-19 Dữ liệu đầu vào module thuỷ lực của mô hình Wadibe-TC 74

Hình 3-20 Kết quả mô phỏng và đo đạc chiều cao sóng - kịch bản đê cao không tràn .76

Hình 3-21 Kết quả mô phỏng và đo đạc chiều cao sóng - kịch bản đê thấp không tườngđỉnh 76

Hình 3-22 Kết quả mô phỏng và đo đạc chiều cao sóng - kịch bản đê thấp có tường đỉnh 77

Hình 3-23 Vận tốc dòng phản hồi trung bình – đê cao không tràn 80

Hình 3-24 Vận tốc dòng phản hồi trung bình – đê thấp không tường đỉnh 81

Hình 3-25 Vận tốc dòng phản hồi trung bình – đê thấp có tường đỉnh 81

Hình 3-26 Kết quả mô phỏng kịch bản 1- D65H15T19_CW0 83

Hình 3-27 Kết quả mô phỏng kịch bản 2 - D65H17T16_CW0 84

Hình 3-28 Kết quả mô phỏng kịch bản 3- D70H15T19_CW0 84

Hình 3-29 Kết quả mô phỏng kịch bản 4 – D70H19T165_CW0 85

Hình 3-30 Kết quả mô phỏng kịch bản 5- D65H15T19_TW0 86

Hình 3-31 Kết quả mô phỏng kịch bản 6- D65H17T16_TW0 86

Hình 3-32 Kết quả mô phỏng kịch bản 7- D70H15T19_TW0 87

Hình 3-33 Kết quả mô phỏng kịch bản 8- D70H19T165_TW0 87

Hình 3-34 Kết quả mô phỏng kịch bản 9- D65H15T19_TW10 88

Hình 3-35 Kết quả mô phỏng kịch bản 10 - D65H17T16_TW10 89

Hình 3-36 Kết quả mô phỏng kịch bản 11- D70H15T19_TW10 89

Hình 3-37 Kết quả mô phỏng kịch bản 12 – D70H19T165_TW10 90

Hình 3-38 Mô hình hóa các kịch bản 92

Hình 3-39 Kết quả mô phỏng xói chân đê cao không tràn 93

Hình 3-40 Kết quả mô phỏng xói chân đê thấp không tường đỉnh 94

Hình 3-41 Kết quả mô phỏng xói chân đê thấp có tường đỉnh 95

Hình 3-42 Kết quả mô phỏng biến đổi lòng dẫn khi có và chưa có sóng phản xạ - Đê cao .96

Hình 3-43 Kết quả mô phỏng biến đổi lòng dẫn khi có và chưa có sóng phản xạ-Đê thấpkhông tường đỉnh 97

Hình 3-44 Kết quả mô phỏng biến đổi lòng dẫn khi có và chưa có sóng phản xạ-Đê thấpcó tường đỉnh 98

Hình 4-1 Thống kê bão đến khu vực nghiên cứu [118] 100

Hình 4-2 Thống kê chiều cao sóng Wavewatch III 101

Hình 4-3 Kết cấu đê biển điển hình ở khu vực Nam Định [89] 103

Hình 4-4 Thiết lập mô hình xói chân đê Thịnh Long 105

Hình 4-5 Kết quả kiểm định xói chân đê biển Thịnh Long – bão Damrey 9/2005 105

Hình 4-6 Ảnh hưởng của chiều cao sóng đến kích thước hố xói chân đê 108

Hình 4-7 Ảnh hưởng của độ sâu nước đến kích thước hố xói chân đê 108

Hình 4-8 Ảnh hưởng của chu kỳ sóng đến kích thước hố xói chân đê 108

Hình 4-9 So sánh chiều sâu hố xói lớn nhất ứng với sự thay đổi của K r0 109

Hình 4-10 Kết quả mô phỏng hố xói chân đê với giải pháp thảm đá rộng 3-6m 111

Hình 4-11 Kết quả mô phỏng hố xói chân đê với giải pháp kè mỏ hàn dài 35-70m 112

DANH MỤC BẢNG BIỂU

Bảng 2-1 Tỉ lệ của một số đại lượng vật lý cơ bản theo tiêu chuẩn Froude [102] 44

Bảng 2-2 Các kịch bản của mô hình lòng cứng 50

Bảng 2-3 Kịch bản thí nghiệm với mô hình lòng động 50

Bảng 3-1 Số liệu biên sóng và mực nước 75

Bảng 3-2 Các tham số được hiệu chỉnh 75

Bảng 3-3 Các kịch bản mô phỏng 92

Bảng 4-1 Thống kê hiện trạng đê biển Nam Định [119] 102

Bảng 4-2 Thống kê chiều cao sóng theo chu kỳ sóng từ số liệu Wavewatch III 106

Bảng 4-3 Điều kiện thuỷ lực ảnh hưởng chiều sâu hố xói lớn nhất 108

Bảng 4-4 Ảnh hưởng của kết cấu đê biển đến chiều sâu hố xói lớn nhất 110

DANH MỤC CÁC TỪ VIẾT TẮT VÀ GIẢI THÍCH THUẬT NGỮ

1 Danh mục các từ viết tắt

N L tỉ lệ dài

SWL Still water level – mực nước tĩnh

CG Current gauge – Vị trí đo dòng chảy

WG Wave gauge – Vị trí đo sóng

DxHyTz tên kịch bản ứng với độ sâu nước D là x (cm), chiều cao sóng H là y(cm)và

chu kỳ đỉnh phổ T là z(10-1s)

CW0 - Đê cao, chiều cao tường đỉnh =0cm

TW0 - Đê thấp, chiều cao tường đỉnh = 0

TW10 - Đê thấp, chiều cao tường đỉnh =10cm

Wadibe -TC (Wave, dike and beach – Time dependent Cross-shore) – Mô hình sóng,

đê và bãi – ngang bờ theo thời gian

2 Giải thích các thuật ngữ

Dòng phản hồi (undertow) là dòng chảy vuông góc với bờ hướng ra phía biển ở dưới

chân sóng

Độ sâu nước tương đối là tỷ số giữa độ sâu nước cục bộ với chiều dài sóng (d/L)

Mô hình lòng cứng – fixed bed model là mô hình vật lý thực hiện mô phỏng chuyển

động của chất lỏng mà không làm thay đổi địa hình lòng dẫn, thường dùng để nghiêncứu sóng, dòng chảy

Mô hình lòng động – mobile bed model là mô hình vật lý thường dùng để nghiên cứu

vận chuyển bùn cát và sự thay đổi địa hình lòng dẫn

MỞ ĐẦU

1 Tính cấp thiết của luận án

Vùng ven biển Việt Nam thường xuyên chịu tác động của thiên tai như bão, áp thấpnhiệt đới, gió mùa, triều cường, nước biển dâng… đặc biệt trong những năm gần đây

xu hướng các hiện tượng cực đoan này xảy ra ngày càng gia tăng cả về tần suất vàcường độ làm ảnh hưởng đến sự phát triển kinh tế, xã hội, cuộc sống của những ngườidân ven biển Do đó hệ thống đê biển luôn được chú trọng xây dựng, duy tu, kiên cố

và nâng cấp hàng năm nhằm bảo vệ vùng đất sau đê khỏi bị ngập lụt dưới tác động từphía biển Đây là giải pháp công trình biển khá phổ biến trên thế giới cũng như ở ViệtNam, đã và đang thu hút nhiều nhà khoa học nghiên cứu để hệ thống đê kè biển làmviệc ổn định lâu dài, hiệu quả, nhất là trong điều kiện cực trị nhằm hạn chế những tácđộng bất lợi dẫn đến các cơ chế hư hỏng đê Một trong các cơ chế hư hỏng đó là xói lởbãi trước và xói chân đê cần được đặc biệt nghiên cứu của khoa học biển hiện nay.Khi bão lớn đổ bộ vào vùng ven bờ, do gió mạnh và chênh lệch áp suất trong bão làmmực nước dâng cao, sóng lớn tiến vào gần đê hơn, đào bới bùn cát và đồng thời tạo ratrường dòng chảy có tốc độ lớn vận chuyển bùn cát đi, gây xói lở bờ biển Trườngdòng chảy ven bờ và các quá trình năng lượng trong đới sóng vỡ có ảnh hưởng quyếtđịnh đến các hiện tượng như suy giảm chiều cao sóng, vận chuyển bùn cát và khuếchtán các vật chất xảy ra do sóng vỡ Khi sóng vỡ, chuyển động sóng biến thành chuyểnđộng rối, đặc trưng bởi các xoáy cuộn có kích thước khác nhau Dưới ảnh hưởng củachuyển động rối do sóng vỡ tạo nên, chuyển động của chất lỏng trong vùng sóng vỡ vôcùng phức tạp, biến đổi mạnh theo cả không gian và thời gian Trong đó, có sự xuấthiện của dòng chảy vuông góc với bờ hướng ra phía biển ở dưới chân sóng, được gọi

là dòng phản hồi Dựa trên những nghiên cứu định lượng về trường dòng chảy ngang

bờ có thể giải thích định lượng được hướng vận chuyển bùn cát, nguyên nhân dẫn đến

sự xuất hiện các cồn ngầm ở bãi trước và xói chân công trình đê biển

Theo các tài liệu nghiên cứu, điều tra khảo sát thực địa cho thấy hầu hết các trườnghợp sự cố ven bờ biển như: xói chân đê, sạt lở, hỏng mái, vỡ đê kè biển thường xảy ratrong thời kỳ bão hoạt động Xói lở xảy ra do ảnh hưởng của bão được xem là xóingắn hạn

là xói chân đê gây hư hỏng, mất an toàn của công trình đê biển như vỡ đê và có thểgây rủi ro ngập lụt [2] Hình 1 và Hình 2 là ví dụ điển hình về hư hỏng mái đê và xóichân đê biển ở Nam Định, Hải Phòng do ảnh hưởng của sóng lớn trong bão Doksuri 9-

2017 và bão Damrey 7-2005 Từ Hình 1a,b có thể thấy ngoài sóng tràn gây xói lở máitrong của đê còn thấy rõ được sự tương tác giữa sóng và tường đỉnh tạo ra sóng bắnlên cao Dưới tác dụng của gió bão làm cho sóng bắn hướng về phía mái trong đê tạo

ra lượng gia tăng của sóng tràn qua đê Ảnh hưởng của sóng bão làm mất ổn địnhtường đỉnh và bị vỡ như Hình 1c

Ngoài những con sóng lớn tràn qua đê còn có những con sóng bị vỡ, trút năng lượngsóng ngay sát chân đê hay trên mái ngoài công trình đê, kè gây mất ổn định cấu kiệnbảo vệ mái ngoài (Hình 1d,e, Hình 2b,c), sự hình thành dòng chảy mạnh do sóng cóthể mang những vật liệu trên mái kè đi nơi khác làm hư hỏng sạt lở mái kè như Hình

1b, có rất nhiều cấu kiện bê tông đúc sẵn bảo vệ mái đê bị sóng vận chuyển lên mặt đê.Hoạt động của sóng bão lớn tạo ra sóng vỡ ở bãi trước và ngay sát chân đê, kè gây racác chuyển động rối làm khuấy động bùn cát lơ lửng, thậm chí các vật liệu chèn bảo vệống buy ở chân kè biển cũng bị trôi đi như Hình 1d,e,f Trên Hình 2a là kết quả đo đạcthực tế mặt cắt ngang trước và sau bão Damrey 2005 tại Nam Định của Viện cơ học vàhình ảnh bãi trước chân đê bị hạ thấp, độ sâu xói chân đê sau cơn bão lên tới 1m, gâynên mất ổn định cho toàn bộ công trình

a

b

Hình 2 Xói lở trước chân đê và sạt lở mái đê, mái và chân gốc kè mỏ hàn sau trận

bão Damrey 7/2005 tại Nam Định [3] [4]

Từ những phân tích trên cho thấy ảnh hưởng của sóng bão đến sự an toàn của đê biển

là rất nghiêm trọng bởi trong bão, mực nước dâng cao, sóng tới lớn, vận chuyển bùncát nhiều và đặc biệt xảy ra vào thời kỳ triều cường là các yếu tố động lực quan trọnggây uy hiếp sự an toàn của đê biển đặc biệt là vấn đề xói chân đê biển Tùy thuộcvào sự

tương tác giữa tải trọng của các yếu tố thủy động lực học với địa hình và kết cấu côngtrình đê biển mà năng lượng sóng bão và dòng chảy tác động gây ra mức độ xói bãitrước và chân đê cũng như mức độ hư hỏng của đê khác nhau Rất nhiều nghiên cứu vàthực tế đã chỉ ra rằng mỗi loại kết cấu công trình có ảnh hưởng nhất định đến kíchthước hố xói chẳng hạn như tường đứng cho độ sâu hố xói lớn hơn so với kè máinghiêng mà một trong những nguyên nhân gây ra sự khác nhau của kích thước hố xói

đó là bởi phản xạ sóng từ các dạng kết cấu mái khác nhau Kết cấu đê biển máinghiêng điển hình ở vùng ven biển Bắc Bộ ảnh hưởng như thế nào đến xói chân đêbiển thông qua sóng phản xạ cho đến nay vẫn chưa được làm sáng tỏ và đây là vấn đềnghiên cứu có ý nghĩa khoa học và thực tiễn giúp cho các nhà thiết kế, quản lý có thểtính toán dự báo được các tác động và đưa ra những giải pháp hỗ trợ phù hợp cho sự

an toàn của đê biển Do đó, đề tài “Nghiên cứu ảnh hưởng của sóng phản xạ đến dòng

phản hồi và xói chân đê biển mái nghiêng khu vực Bắc Bộ” được lựa chọn để nghiên

cứu trong Luận án này

2 Mục tiêu nghiên cứu

+ Đánh giá được ảnh hưởng của sóng phản xạ từ mái đê biển đến dòng phản hồi và xóichân đê biển;

+ Mô phỏng xói chân đê biển do dòng phản hồi với ảnh hưởng của sóng phản xạ từ mái

đê và ứng dụng tính xói chân đê biển tại Nam Định

3 Đối tượng và phạm vi nghiên cứu

3.1 Đối tượng nghiên cứu

Dòng phản hồi và xói chân công trình đê biển mái nghiêng dưới tác động của sự tươngtác giữa sóng bão với mái đê thông qua sóng phản xạ

3.2 Phạm vi nghiên cứu

Xói chân đê do cơ chế dòng phản hồi cho hệ thống đê biển mái nghiêng vùng ven biểnBắc Bộ

4 Nội dung nghiên cứu

Nội dung chính của Luận án bao gồm các phần như sau:

 Nghiên cứu ảnh hưởng của kết cấu mái công trình đê biển đến phân bố sóng phản

xạ và đặc trưng cấu trúc của dòng phản hồi theo mặt cắt ngang trong bão bằng việcthiết lập mô hình vật lý lòng cứng trong máng sóng;

 Nghiên cứu ảnh hưởng của kết cấu mái công trình đê biển đến xói chân đê biểnbằng mô hình vật lý lòng động;

 Phát triển cập nhật mô hình toán mô phỏng dòng phản hồi và xói chân đê có xét đếnảnh hưởng của kết cấu mái công trình đê biển và kiểm định bằng số liệu thí nghiệmứng với các điều kiện nghiên cứu trong máng sóng và ứng dụng nghiên cứu đối vớivùng biển Nam Định

Với nội dung nghiên cứu trên, Luận án trả lời 4 câu hỏi:

1) Hệ số phản xạ sóng thay đổi như thế nào trên mặt cắt ngang bãi ứng với các cấu trúc

đê biển mái nghiêng?

2) Sóng phản xạ ảnh hưởng như thế nào đến cấu trúc của dòng phản hồi và vận tốcdòng chảy trung bình?

3) Ảnh hưởng của sóng phản xạ đến vận chuyển bùn cát và xói lở chân đê diễn ra nhưthế nào?

4) Ứng dụng các nội dung nghiên cứu cho khu vực cụ thể ở Nam Định cho kết quảnhư thế nào?

5 Cách tiếp cận và phương pháp nghiên cứu

5.1 Cách tiếp cận

Cách tiếp cận được áp dụng trong Luận án: Quan sát hiện trạng thực tế vấn đề xói chân

đê dọc ven biển Bắc Bộ, tiến hành tổng hợp, phân tích và kế thừa những công trìnhnghiên cứu đã thực hiện trên thế giới và trong nước Từ đó đưa ra những khía cạnhkhoa học chưa được nghiên cứu kỹ và ứng dụng nghiên cứu trong điều kiện thực tếcủa Việt Nam Sau đó áp dụng các phương pháp có thể để tìm hiểu sâu hơn, giải thích

rõ hơn và tính toán mô phỏng định lượng đối tượng nghiên cứu

5.2 Phương pháp nghiên cứu

Để đạt được mục tiêu trên, Luận án sử dụng tổng hợp một số phương pháp nghiên cứusau đây:

- Phương pháp thống kê: Tổng hợp và phân tích các tài liệu đã có nhằm tổng kết và

kế thừa các kết quả nghiên cứu trong và ngoài nước về dòng chảy trong vùng sóng

vỡ nói chung và dòng phản hồi nói riêng, vận chuyển bùn cát ngang bờ và sự thayđổi địa hình đáy trong bão

- Phương pháp thực nghiệm: Thí nghiệm mô hình vật lý trong máng sóng ứng với cácđiều kiện khác nhau để mô phỏng cấu trúc dòng phản hồi và quá trình vận chuyểnbùn cát khu vực chân kè biển ứng với các kịch bản khác nhau

- Phương pháp mô hình toán: Phát triển có kế thừa mô hình toán học Wadibe-TC môphỏng dòng phản hồi và diễn biến xói chân kè trong bão trong điều kiện nghiên cứu,kết hợp với kết quả mô hình vật lý để kiểm nghiệm

- Phương pháp chuyên gia: Thực hiện công tác xin ý kiến đóng góp của các chuyêngia trong lĩnh vực biển về hướng nghiên cứu, phương pháp thực hiện và kết quảnghiên cứu trong suốt quá trình thực hiện

6 Ý nghĩa khoa học và thực tiễn của luận án

Mặc dù đã có rất nhiều công trình nghiên cứu về vấn đề dòng chảy và xói lở ngang bờ,nhưng diễn biến cơ chế của hiện tượng này rất phức tạp và rất khó để dự tính, dự báochính xác được với những điều kiện biên luôn biến động và hiếm số liệu đo đạc thực

tế Chính vì vậy đây là hướng nghiên cứu có ý nghĩa khoa học và có khả năng ứngdụng vào thực tế cao

6.1 Ý nghĩa khoa học

- Đóng góp về mặt khoa học đối với các nghiên cứu ảnh hưởng của công trình đê biểnmái nghiêng đến phân bố sóng phản xạ và dòng phản hồi trên mặt cắt ngang trong điềukiện bão;

- Đóng góp về mặt khoa học đối với nghiên cứu diễn biến xói lở chân kè dưới ảnhhưởng của kết cấu công trình đê biển mái nghiên điển hình ở Việt Nam thông qua việctính toán ảnh hưởng của sóng phản xạ;

- Đóng góp về mặt khoa học đối với công cụ mô hình toán nhằm dự tính, phân tíchthủy động lực học và xói chân đê, kè biển do sóng bão tốt hơn

Phần mở đầu

C1 Tổng quan nghiên cứu dòng phản hồi do sóng và xói chân đê biển

Đê biển mái

nghiêng Mô hình Điều kiện sóng bão

Sóng phản xạ C2 Cơ sở khoa học mô hình hoá dòng phản hồi do sóng và xói chân kè

Dòng phản hồi

C3 Kết quả nghiên cứu ảnh hưởng của sóng phản xạ đến dòng phản hồi và xói chân đê biển mái nghiêngTừ khóa: Hệ số xáo trộn

bùn cát C4 Ứng dụng kết quả nghiên cứu tính toán xói chân đê biển Nam Định Xói chân đê

Kết luận và kiến nghị

6.2 Ý nghĩa thực tiễn

- Ứng dụng nghiên cứu này cho bài toán thực tế có thể cho phép đánh giá, dự tính địnhlượng xác thực hơn lượng vận chuyển bùn cát ngang bờ và xói lở chân kè trong cơnbão do dòng phản hồi chi phối Nội dung nghiên cứu này có thể là cơ sở cho việc thiết

kế, xây dựng các giải pháp bảo vệ bờ bãi, chân công trình biển phù hợp nhất;

- Là cơ sở cho các nghiên cứu ứng dụng tiếp theo về vận chuyển bùn cát ven bờ

7 Cấu trúc luận án

Với những mục tiêu và phương pháp như trên, cấu trúc nội dung của Luận án đượctrình bày trong các chương như Hình 3

Ngoài phần mở đầu giới thiệu tổng quát về Luận án và phần kết luận, kiến nghị, Luận

án được trình bày trong 4 chương Chương 1 trình bày tổng quan tình hình nghiên cứu

về dòng phản hồi và xói chân đê biển từ đó đề xuất ra hướng nghiên cứu cho Luận án

Cơ sở khoa học và thiết lập mô hình vật lý cũng như mô hình toán nghiên cứu ảnhhưởng của sóng phản xạ đến dòng phản hồi và vận chuyển bùn cát, biến đổi địa hìnhngang bờ được thể hiện trong chương 2 Những kết quả nghiên cứu của Luận án vềsóng phản xạ, dòng phản hồi và xói chân đê biển mái nghiêng được nêu trong chương

3 Chương 4 là phần ứng dụng kết quả nghiên cứu vào tính toán mô phỏng cho trườnghợp cụ thể xói chân đê biển Nam Định

Hình 3 Cấu trúc của Luận án

CHƯƠNG 1 TỔNG QUAN NGHIÊN CỨU DÒNG PHẢN HỒI DO

SÓNG VÀ XÓI CHÂN ĐÊ BIỂN1.1 Giới thiệu chung

Chuyển động của các phần tử nước trong sóng làm khuấy động bùn cát, dịch chuyểntheo dòng chảy và lượng vận chuyển bùn cát tăng lên do sự có mặt của sóng Chươngnày trình bày tổng quan các nghiên cứu đã thực hiện về sóng, dòng phản hồi và vậnchuyển bùn cát, xói bồi trước chân đê biển Từ đó đưa ra những vấn đề nghiên cứu tiếptheo trong Luận án

1.1.1 Dòng phản hồi và các quá trình vận chuyển bùn cát ngang bờ

Có rất nhiều yếu tố chi phối sự hình thành dòng chảy ven bờ nói chung và dòng phảnhồi nói riêng cũng như vận chuyển bùn cát gồm thủy triều, sóng, gió, đặc điểm thủyvực, đặc điểm của bão, đặc điểm bùn cát… Trong phần này giới thiệu tổng quát về cơchế vật lý hình thành dòng phản hồi và vận chuyển bùn cát ngang bờ dẫn đến xói chân

đê biển do sóng

1.1.1.1 Dòng phản hồi

Dòng phản hồi là dòng chảy trung bình hướng ra biển, ở phía dưới mực nước chânsóng trong vùng sóng vỡ [7] Dòng phản hồi được xem như là một trong những lựcchính chi phối vận chuyển bùn cát ra xa bờ, hình thành các cồn cát, xói lở bãi biển vàxói lở đê kè biển khi có bão Để giải thích cơ chế hình thành dòng này cần xuất phát từcác lực làm dịch chuyển các phần tử nước trong chuyển động sóng như dưới đây.Theo Stive and Wind (1996) và Svendsen (1984), dòng phản hồi được hình thành do

sự mất cân bằng theo phương thẳng đứng giữa thông lượng động lượng biến đổi theo

độ sâu và gradient áp suất đồng nhất theo độ sâu do sự dâng nước bởi sóng Hay nóicách khác thì dòng phản hồi là dòng chảy hướng ra biển ở phần dưới cột nước khi đền

bù cho thông lượng khối lượng hướng vào bờ ở phía trên mực nước chân sóng [8] [9][10] Trước khi vào vùng sóng vỡ, các phần tử nước chuyển động tiến và lui theochuyển động quỹ đạo sóng, năng lượng sóng bị tiêu tán ít, nhưng đến độ sâu nướcbằng nhỏ hơn hoặc bằng khoảng 1,25 lần chiều cao sóng thì quá trình sóng vỡ diễn

ra làm tiêu hao

năng lượng sóng rất nhiều, thông lượng năng lượng sóng chuyển thành năng lượng rối

và nhiệt [11] Các loại sóng vỡ gồm vỡ cuộn (spilling), vỡ đổ (plunging) và vỡ dâng(surging) xuất hiện và dịch chuyển tùy thuộc vào đặc trưng quan trọng nhất đó là độdốc của sóng tới và của bãi biển [12] [13] Một phần năng lượng sóng vỡ sẽ chuyểnhóa thành động lượng gia tăng gọi là ứng suất bức xạ của sóng và tạo thành sóng cuộn(sóng cồn) nhiễu loạn có mặt sóng gần như dốc đứng vận chuyển nước vào bờ [14].Lượng nước do sóng cuộn mang vào khu vực ven bờ sẽ làm cho mực nước trung bìnhven bờ dâng cao và sinh ra dòng chảy ngang bờ do sự cân bằng về khối lượng

Cho đến nay đã có rất nhiều nghiên cứu được phát triển nhằm tính toán định lượngdòng phản hồi do sóng bằng mô hình mô phỏng Tất cả các nghiên cứu đều cho rằngdòng phản hồi được tạo ra do sự chênh lệch cục bộ giữa ứng suất bức xạ biến đổi theo

độ sâu gây ra bởi sóng và gradient áp suất đồng nhất theo độ sâu do độ dốc bề mặtnước (Hình 1-1 (A) ) hoặc do sự mất cân bằng lực giữa thông lượng động lượng sóng

và nước dềnh do sóng [15]… Chính sự mất cân bằng lực này đã cung cấp lực tạo rahoàn lưu theo phương thẳng đứng

Hình 1-1 (A) Sơ họa về phân bố gradient ứng suất bức xạ (Sxx/x) và áp lực (/x), (B) kết

quả phân bố theo phương đứng của dòng chảy ngang bờ [16]

Khi dòng phản hồi ổn định được thiết lập chính là sự cân bằng của ba lực quan trọnggồm ứng suất bức xạ, gradient áp lực do nước dềnh và ứng suất rối được tạo ra bởi sựbiến đổi theo phương thẳng đứng của vận tốc dòng phản hồi kết hợp với chuyển động

rối xuất hiện do sóng vỡ Các yếu tố chi phối cân bằng động lượng gồm: a) bảo toànkhối lượng cục bộ theo phương thẳng đứng, b) ứng suất tiếp tại vị trí chân sóng, c) ứngsuất tiếp tại đáy, d) biến đổi của hệ số nhớt rối (turbulence) theo phương thẳng đứng.Trong số những yếu tố này thì bảo toàn khối lượng cục bộ xác định đến độ lớn củadòng phản hồi trong sóng vỡ Thông lượng khối lượng hướng về phía bờ do sóng trongvùng từ đỉnh sóng đến chân sóng, bao gồm cả sự đóng góp của sóng vỡ cuộn (breakingwave rollers), được bù bởi dòng phản hồi tổng cộng hướng ra xa bờ Ba nhân tố còn lạiquyết định đến cấu trúc của dòng phản hồi theo phương thẳng đứng trên mặt cắt.

Thông lượng khối lượng tịnh của dòng hướng bờ thường cân bằng với dòng phản hồihướng ra biển, tuy nhiên trong trường hợp sóng tràn trong bão, cân bằng khối lượngnày bị thay đổi do lượng nước tràn mất trong đất liền Ảnh hưởng này thường đượctính gộp với sóng phản xạ trong quá trình tính toán bằng các mô hình

1.1.1.2 Các quá trình vận chuyển bùn cát ngang bờ

Vận chuyển bùn cát phụ thuộc rất nhiều yếu tố gồm chuyển động của các phần tử nướcdao động dưới điều kiện sóng, lực tạo dòng chảy, đặc điểm thủy vực và bùn cát [1].Quá trình dịch chuyển của các phần tử bùn cát bao gồm 3 giai đoạn: (1) Các hạt bùncát đáy bị khuấy động và mang lên khỏi đáy ở trạng thái lơ lửng; (2) Những thànhphần hạt này bị dịch chuyển theo phương ngang theo dòng nước; (3) Các hạt bùn cát bịlắng đọng xuống Theo Dean (1973), nếu vận tốc lắng chìm cao, hạt sẽ chìm xuốngtrước khi chân sóng xuất hiện và sẽ dịch chuyển hướng vào bờ, trong khi nếu vận tốclắng chìm thấp, nó vẫn sẽ lơ lửng và bị dịch chuyển ra xa bờ Điều này dẫn đến kết quả

là các hạt thô hơn bị lọc về phía bờ Diễn biến của mỗi giai đoạn tùy thuộc vào sự dịchchuyển của nước/ dòng chảy và đặc điểm của hạt cát

Vận chuyển bùn cát được xác định là lượng bùn cát (S) bị dịch chuyển đi theo dòngchảy trên đơn vị diện tích trong đơn vị thời gian, nghĩa là giả định rằng vận tốc chuyểnđộng của hạt bùn cát tương ứng với vận tốc của hạt chất lỏng và công thức tổng quát

để tính vận chuyển bùn cát ngang bờ bằng tích của nồng độ bùn cát (c) với vận tốc (u)[17] Như vậy để tính toán vận chuyển bùn cát ngang bờ trong bão cần phải xác địnhđược phân bố theo phương thẳng đứng của cả nồng độ bùn cát và dòng chảy Hình 1-2

trình bày sơ họa mô phỏng phân bố vận tốc dòng chảy và nồng độ bùn cát trong sóng

vỡ Có

thể thấy rằng dòng vận chuyển bùn cát đi vào bờ ít hơn dòng vận chuyển bùn cát ra xa

bờ bởi rõ ràng rằng nồng độ bùn cát phân bố theo độ sâu chịu ảnh hưởng của hiệntượng rối loạn và nhiễu động ở mức độ mạnh do hiện tượng sóng vỡ, kích thước hạtcát, tốc độ lắng chìm Do tác động của yếu tố trọng lực mà nồng độ bùn cát ở phíadưới sẽ cao hơn phía trên mặt nước, dẫn tới vận chuyển bùn cát tịnh sẽ có hướng raphía ngoài khơi với sự có mặt của dòng phản hồi

Hình 1-2 Sơ họa phân bố dòng chảy và nồng độ bùn cát vùng sóng vỡ [18]

Khi tính toán các quá trình vận chuyển ngang bờ thường được chia thành các lớp khácnhau: trên mặt nước là dòng trôi do sóng tiến vào bờ (surface drift), tiếp xuống dưới đó

là dòng phản hồi và phía đáy là dòng biên sát đáy (streamming) Vận chuyển bùn cátcũng được chia thành vận chuyển bùn cát lơ lửng (suspended load) và vận chuyển bùncát đáy (bed load) Trong điều kiện bão, ở vùng sóng vỡ gần bờ, vận chuyển bùn cát lơlửng chiếm ưu thế [19] [1] [20]

1.1.2 Xói chân đê biển

Xói chân đê biển là các hiện tượng chỉ sự hạ thấp cao độ đáy khu vực chân đê dưới tácdụng của các yếu tố thuỷ lực biển Sự thay đổi mặt cắt ngang, các hố xói hay các cồncát ngầm dọc bờ biển có thể được hình thành và các cồn cát ngầm có thể bị di chuyển

vị trí hoặc bị phá hủy tùy thuộc vào gradient vận chuyển bùn cát ven bờ

Hình dạng mặt cắt ngang trước chân đê biển có thể biến đổi trong một trận bão, mộtmùa hay một năm và nhiều năm Các yếu tố ảnh hưởng đến hố xói chân đê gồm hiệntượng dâng lên của mực nước biển trong thời kì bão, tác động của sóng lớn trong bão,dòng chảy ven bờ, tác động thường xuyên của thủy triều, đặc biệt triều cường, địa hình

phía trước chân đê thoải hay dốc, ngoài ra cần kể đến sự gia tăng xói lở do các sóngphản xạ từ kết cấu công trình đê, sự gia tăng năng lượng sóng do thềm bãi bị hạ thấp,

do ảnh hưởng của quá trình vận chuyển bùn cát dọc bờ, do sự gia tăng bão biển (khiliên tiếp xảy ra bão trong thời gian ngắn thì tình hình xói lở trở nên nghiêm trọng, docơn bão thứ hai và cơn bão kế tiếp xảy ra trên bãi biển chưa kịp phục hồi do xói lở sautrận bão đầu tiên Ví dụ trên mặt cắt (Hình 2) ở bãi biển Nam Định năm 2005 là một ví

dụ điển hình, sau ba trận bão liên tiếp, chiều sâu hố xói tại chân đê lên đến khoảng1m)

Dưới tác dụng của các yếu tố trên, bùn cát trước chân đê biển bị tác động làm lỏng ra

và dịch chuyển đi theo dòng chảy Kết quả mặt cắt thay đổi như thế nào ở các vị tríkhác nhau trên mặt cắt là xói hay bồi tùy theo chênh lệch giữa lượng bùn cát được vậnchuyển đi và đến vị trí đó Trong điều kiện bão, mực nước biển dâng cao, sóng lớn tạo

ra nước dềnh do sóng cao, gió mạnh, sóng vỡ gần chân đê hơn dẫn đến dòng phản hồi

ra biển lớn hơn làm cho lượng vận chuyển bùn cát ra khỏi chân lớn hơn bùn cát đến(do phía trên mái đê được kiên cố hóa, lượng vận chuyển bùn cát không có trong khisuất chuyển cát lớn, nó sẽ đào bới lấy bùn cát đáy đem đi ra xa bờ) gây xói chân đê.Như vậy hố xói hình thành do sự mất cân bằng bùn cát cục bộ tại vị trí chân đê biểnkhi mái kè ngăn giữ không cho vật liệu thân đê bị sóng/dòng chảy mang đi Khi đếnngoài vùng sóng vỡ, dòng phản hồi nhỏ hơn, bùn cát lắng chìm xuống và tạo ra cáccồn cát ngầm Đối với các trận bão lớn, hình dạng mặt cắt thay đổi nhiều hơn, các bãicát ngầm vùng sóng vỡ có thể được hình thành cao hơn

1.2 Tổng quan nghiên cứu về dòng phản hồi

Theo nghiên cứu của Svendsen (1984), dòng phản hồi do sóng được phát hiện và mô

tả chi tiết từ cách đây hàng trăm năm bởi Johnson (1919) [21] Khi nghiên cứu sự hìnhthành bãi do sóng trong mô hình thí nghiệm máng sóng Bagnold (1940) đã thực hiệncác thí nghiệm quan trắc dòng ngang bờ và xác định dòng phản hồi [22] Nielsen vàSorensen (1970) đưa ra giải thích một cách định tính dựa trên phân tích lý thuyết vềhiện tượng dòng phản hồi [23] Đánh giá định lượng về ý tưởng của Nielsen vàSorensen (1970) dẫn đến các mô hình mô phỏng dòng phản hồi dựa trên phân bố theophương thẳng đứng của ứng suất bức xạ được thực hiện bởi Svendsen (1984) [8] [24]

Trong điều kiện cực trị, để có được bộ số liệu đo đạc hiện trường chi tiết đồng thời vớiphân bố dòng chảy trên mặt cắt ngang bờ và sóng để nghiên cứu dòng phản hồi là điềurất khó khăn trên toàn thế giới do hạn chế về thiết bị và điều kiện đo đạc Các đo đạcdòng chảy, sóng tương ứng với sự biến đổi địa hình trên mặt cắt ngang trong bão thuthập được trên các bài báo, báo cáo đã xuất bản chỉ có từ một số dự án nghiên cứu lớncủa các nước tiên tiến trên thế giới như Mỹ, Úc, Canada, Nhật… [25] [16] [26] [27][28]

[29] [30] [31] Chính vì vậy khi nghiên cứu dòng phản hồi theo các khía cạnh khácnhau nhằm dự tính cấu trúc của dòng phản hồi theo mô hình toán, các chuyên giathường sử dụng số liệu trong phòng thí nghiệm để phân tích tìm ra quy luật, đồng thời

sử dụng số liệu đo đạc trong mô hình vật lý đó để kiểm định mô hình [32] [33] [15][35] [37] [38]

[39] [40] [10] [41] [43]… Các thí nghiệm mô hình vật lý trong các nghiên cứu đã thựchiện chỉ giới hạn trong một số kịch bản nhất định về biên địa hình và biên thủy lực

Mô hình hóa được nghiên cứu phát triển, hoàn thiện trong nhiều năm và đã được ứngdụng rộng rãi [15] [44] [45] [46] [39]… Cách tiếp cận chung của mô hình được bắtnguồn từ các phương trình cơ bản mô phỏng dòng nước trong vùng sóng vỡ bao gồmphương trình cân bằng khối lượng (hay phương trình liên tục) và phương trình cânbằng động lượng Sự khác biệt giữa các mô hình hóa dòng phản hồi này là kỹ thuật môhình hóa để tính toán các khía cạnh của quá trình vật lý khi giải các phương trình nhưcác giả thiết, lựa chọn điều kiện biên, đặc trưng độ nhớt, xử lý lớp biên dưới đáy, loại

lý thuyết sóng [10] [47]

Trong phương trình chuyển động cơ bản, thành phần lực chủ đạo ở vùng sóng vỡ làthông lượng động lượng sóng vỡ (ứng suất bức xạ), chính là thành phần sinh ra dòngphản hồi Ma sát đáy là thành phần lực cản dòng chảy Phân bố ngang bờ của dòngphản hồi trở lại biển ở dưới mực nước chân sóng u(z) được xác định dựa trên phân bốcủa ứng suất tiếp rối theo phương thẳng đứng z theo công thức:

τ(z) = ρ (z) du

Trong đó  là khối lượng riêng của nước; t(z) là hệ số nhớt rối hay hệ số xáo trộn củachất lỏng

Giả thiết cơ bản đối với phương trình này là chỉ tính đến độ lệch của phân bố áp lựcthủy tĩnh do sóng gây ra Để giải u(z) thì phân bố độ nhớt theo phương thẳng đứng

t(z) được giả định và tích phân phương trình trên hai lần theo phương thẳng đứngcùng với các điều kiện biên phù hợp để xác định các hàm tích phân [48] [19]

Svendsen (1984) đã tính toán ứng suất tiếp trong trường hợp t = const (đồng nhất theo

độ sâu) và t biến đổi theo hàm mũ, trong khi Steetzel (1993) và Okayasu (1989) lại ápdụng quan hệ tuyến tính [19] [15], còn De Vriend (1987) chọn giả sử t có dạng phân

bố parabol [51] [52] Theo Steetzel (1993), sau khi tích phân theo phương thẳng đứng,

phương trình tính toán phân bố của dòng phản hồi u(z) trên mặt cắt ngang được viết

như sau [19]:

az u(z) = u0 + Klinz + Klogln [1 +

ε

Trong đó, u0 là vận tốc dòng chảy trung bình ngang bờ (m/s), Klin (s-1) và Klog (m/s) làhai hằng số liên quan đến thành phần tuyến tính và logarit, �0 là hệ số xáo trộn tại mựctham chiếu (m2/s), a là gradient xáo trộn theo phương đứng = /z

Nam (2013) khi nghiên cứu tính ảnh hưởng của độ lệch sóng đến dòng phản hồi đã sửdụng công thức phân bố dòng chảy theo phương thẳng đứng của Rattanapitikon vàShibayama (2000) cũng có dạng tương tự [39]:

u(z) = u0 + α1 (

ρ ) [α2 (d − 2) + α3 (ln d + 1)]

(1-3)

Trong đó �1, �2 , �3 là các hệ số, Dtotal là tổng năng lượng tiêu tán do sóng vỡ và

ma sát, dt là độ sâu dưới nước chân sóng

Svendsen (1984) đã tập trung vào nghiên cứu vai trò của chuyển động rối trong độnglực dòng phản hồi và tham số hóa yếu tố này trong hình thành của cuộn sóng (waveroller) Vì cuộn sóng trong sóng vỡ được chuyển từ mặt sóng tới (wave front) nên nóảnh hưởng tới sự dịch chuyển khối lượng Hơn nữa ứng suất bức xạ giữa cuộn sóng vàmặt sóng (wave surface) tác động lên phân bố dòng phản hồi theo phương thẳng đứng[8]

Trong phương trình liên tục, dòng phản hồi trung bình độ sâu nước được bù đắp bởithông lượng thể tích hướng vào bờ do sóng và cuộn sóng hoặc các sóng chưa vỡ.Bằng cách bỏ qua các thông lượng do dòng chảy sát đáy (streaming), thông lượng khốilượng tổng cộng (mt) do sóng vỡ bao gồm thông lượng khối lượng do chuyển độngsóng (mw) và thông lượng khối lượng do cuộn sóng trên mặt nước (mr) Các đại lượngnày được xác định bằng nhiều công thức khác nhau dựa trên giả thiết về phân bố vậntốc dựa trên mô hình của Fredsoe & Deigaard (1992) [44] [45] [19] [15].

Theo Yoshiaki và Toshiyuki (2000), có hai phương pháp để tính toán diện tích cuộnsóng trên mặt nước [45] Một phương pháp giả thiết rằng diện tích cuộn sóng tỉ lệ vớibình phương chiều cao sóng: Ar1= aH2 trong đó Svendsen (1984b) đưa ra a = 0,9;

Battjes (1985) [48] Buhr Hansen (1984) đề xuất a = 0,2 ÷ 0,4; Mặt khác Okayasu(1989) [15] cho thấy Ar có thể được tính bằng 0,06HL Phương pháp thứ hai giả thiết

rằng năng lượng của chuyển động sóng được chuyển thành năng lượng của cuộn sóngmặt ở trong sóng vỡ mà không bị tổn thất và sau đó năng lượng cuộn sóng mặt bị tiêután Có một điểm chung trong toàn bộ các nghiên cứu đó trong việc sử dụng định luậtbảo toàn khối lượng cục bộ theo phương thẳng đứng như một điều kiện biên [31].Trong kĩ thuật này, thông lượng khối lượng mt (kg/m/s) ở lớp trên mặt được cân bằngvới tổng lượng dòng phản hồi mà được tích phân từ đáy đến mực nước chân sóng:

Mô hình tính toán dòng phản hồi được áp dụng theo các bước sau: Từ số liệu đầu vàogồm chiều cao sóng (H0), chu kỳ sóng (T), điều kiện địa hình (zb) và mực nước, tiếnhành tính toán phân bố năng lượng sóng trên mặt cắt ngang, thông lượng khối lượng,ứng suất tại chân sóng và đáy, từ đó xác định được trắc diện dòng phản hồi Trong quátrình tính toán có thể lựa chọn các công thức phù hợp với đặc điểm của khu vực nghiêncứu, bởi sử dụng các công thức khác nhau sẽ cho các kết quả khác nhau như trongtrường hợp nghiên cứu của Hedegaard (1992) và Salmon (2016) Hầu hết các kết quảcủa mô hình được hiệu chỉnh và kiệm định phù hợp với số liệu đo đạc trong phòng thínghiệm như các nghiên cứu của Svendsen (1984), Stive & Wind (1986), Cox &Kobayashi (1996), Steetzel (1993), Hedegaard (1992), Svendsen & Hansen (1988),Okayasu (1989)

v.v Một số mô hình được thực hiện với cả số liệu hiện trường và số liệu trong phòngthí nghiệm như nghiên cứu của các tác giả Yoshiaki & Toshiyuki (2000), Brian (1998),Grasmeijer & Ruessink (2003) Các mô hình tính dòng phản hồi hiện tại cho kết quảkhá tốt với những số liệu đo đạc được đó Nhưng các số liệu đo đạc trong phòng thínghiệm phục vụ tính toán dòng phản hồi đó được thực hiện trong một máng sóng kínnghĩa là khối lượng nước được bảo toàn mà chưa tính đến quá trình tổn thất do sóngtràn đi và các kịch bản nếu có công trình thì thường là công trình tường đứng

Những nghiên cứu gần đây cho thấy nỗ lực cải thiện mô hình số về dòng phản hồi để

có kết quả sát với thực tế hơn, nhưng các quá trình thủy động lực học vùng sóng vỡ rấtphức tạp, đặc biệt là trong điều kiện bão, mực nước dâng cao, sóng lớn hơn vỡ sát bờ,sóng tràn xảy ra khá phổ biến đối với đê đỉnh thấp như ở Việt Nam thì cũng cần thiếtphải nghiên cứu xem xét tính toán đến hiện tượng sóng tràn và sóng phản xạ ảnhhưởng đến dòng phản hồi như thế nào

Sự thay đổi phân bố cấu trúc dòng phản hồi do sóng phản xạ của bãi biển đượcNeshaei (2009) mô phỏng dựa trên mô hình cho sóng đều của Okayasu (1989),Neshaei (2009) Thông lượng khối lượng tổng cộng của sóng vỡ và sóng phản xạ tínhtheo công thức mt

= (mw+mr)i – (mw)r với chiều cao sóng phản xạ bằng hệ số phản xạ nhân với chiều caosóng Tương tự như vậy thì ứng suất bức xạ Sxx do sóng cũng được tính bằng cáchcộng thêm ứng suất do sóng phản xạ: Sxx = (Sxx)i + (Sxx)r

Từ những vấn đề nêu trên có thể thấy sự cân bằng của thông lượng khối lượng do sónggây ra là nguyên tắc để xác định trắc diện dòng phản hồi trung bình, đại lượng có ýnghĩa quan trọng nhất của bất kỳ mô hình quá trình tính toán xói lở do bão gây ra vàhình thành các cồn cát ngầm ở phía trước đê biển Sự xuất hiện sóng tràn có xu hướnglàm giảm dòng phản hồi và do đó làm giảm xói trước chân đê Tuy nhiên, quan sát trênthực địa với một số loại cấu trúc mái dốc đê biển ở Việt Nam cho thấy sóng tràn có thểgây ra nhiều rối hơn ở chân phía biển của đê mái nghiêng, làm tăng cường độ dòngchảy và gây xói lở nhiều hơn Lập luận này nên được làm sáng tỏ bởi vì nhiều nghiêncứu về hiện tượng tràn sóng đã được phát triển, nhưng tác động của sóng tràn đếndòng phản hồi vẫn chưa được làm sáng tỏ.

Nghiên cứu của Neshaei (2009) và Goda & Suzuki (1976) đã chỉ ra rằng khi hệ sốphản xạ của bãi tăng thì độ lớn của dòng phản hồi giảm mà hệ số phản xạ giảm dần từmép nước ra ngoài khơi và mô hình Neshaei và cộng sự (2009) cũng chỉ mới kiểmđịnh với ba kịch bản sóng đều trong máng nhỏ, chưa tính đến sự tương tác của sóngvới công trình đê và đặc biệt đê biển có tường đỉnh, đồng thời sự suy giảm của hệ sốphản xạ theo quan hệ nào thì nghiên cứu cũng chưa có đề cập tới

1.3 Tổng quan mô hình vận chuyển bùn cát và xói lở ngang bờ

Trong những thập kỷ gần đây đã có rất nhiều nghiên cứu thuỷ lực, vận chuyển bùn cát

và xói lở ngang bờ từ đó phát triển mô hình toán mô phỏng từ đơn giản đến phức tạp.Các mô hình đơn giản như mô hình kinh nghiệm mặt cắt xói sau bão như nghiên cứuBruun (1954) với mặt cắt cân bằng dạng hàm mũ, Dean (1987, 1991) mô tả hệ số hìnhdạng mặt cắt liên hệ với mức độ tiêu hao năng lượng sóng trong vùng sóng đổ và tínhchất hạt bùn cát Đơn giản hơn là công thức tính chiều sâu hố xói tại vị trí chân côngtrình là hàm của đặc trưng sóng, độ sâu nước và độ dốc như nghiên cứu của Xie(1981), Hughes và Fowler (1992), Mc Dougal (1996) Các mô hình phổ biến được ápdụng hiện nay là các mô hình mô phỏng quá trình dựa trên các quá trình vật lý nhưsóng bất đối xứng, cân bằng thông lượng khối lượng, dòng phản hồi, dòng biên sátđáy… Hầu hết các mô hình này đều phân chia hai kiểu vận chuyển bùn cát đó là vậnchuyển bùn cát lơ lửng và vận chuyển bùn cát đáy Trong điều kiện sóng bão, sự khuấyđộng bùn cát diễn ra mạnh mẽ nên vận chuyển bùn cát lơ lửng do dòng phản hồi chiếm

ưu thế và

là nguyên nhân dẫn đến biến đổi hình thái mặt cắt ngang bờ nên Luận án tập trungnghiên cứu vận chuyển bùn cát lơ lửng.

Cấu trúc của các mô hình vận chuyển bùn cát và biến đổi hình thái thường bao gồm 4modul: (1) module sóng: tính toán truyền sóng ngang bờ (2) module dòng chảy: tínhphân bố dòng chảy do sóng (3) module vận chuyển bùn cát: tính toán tải lượng bùn cát

và (4) module hình thái: mô phỏng diễn biến mặt cắt ngang Mục 1.2 đã trình bày tổngquan về dòng chảy, phần dưới đây sẽ trình bày tổng quan về các module còn lại

1.3.1 Tổng quan mô hình sóng

Hiện nay có nhiều mô hình sóng đã được xây dựng để mô phỏng quá trình truyền sóng

từ vùng nước sâu đến vùng nước nông bằng cách sử dụng lý thuyết sóng tuyến tínhhay phi tuyến ứng với việc mô phỏng những con sóng đều và mô phỏng những consóng ngẫu nhiên đối với địa hình bãi biển bất kì hay thoải [54] Mỗi mô hình có giớihạn mô phỏng các quá trình vật lý của sóng gió trong vùng ven bờ bao gồm các quátrình hiệu ứng nước nông, khúc xạ, nhiễu xạ, sóng vỡ, sóng vỡ bạc đầu (whitecapping), tương tác sóng (triad and quadruplet), suy giảm năng lượng sóng do rừngngập mặn và ma sát đáy, công trình… mỗi mô hình có ưu, nhược điểm và sự phù hợpriêng, khi lựa chọn ứng dụng mô hình cần phải phân tích dựa trên bài toán cụ thể củavùng nghiên cứu [54] [55]

Mô hình cân bằng năng lượng sóng được phát triển sớm nhất và cập nhật các quá trìnhvật lý sóng qua các thế hệ khác nhau Dano Roelvink (2012) cho rằng thế hệ đầu tiênxem xét sự biến đổi của phổ tần với số liệu đầu vào là gió và hiện tượng sóng bạc đầu(while capping), thế hệ thứ hai có tính đến ảnh hưởng của sự tương tác giữa sóng phituyến, đồng thời thêm các quá trình như lan truyền theo hướng sóng, khúc xạ, hiệu ứngnước nông, ma sát đáy, tương tác sóng-dòng chảy, sóng vỡ do độ sâu nước hạn chế vớiđịa hình bất kỳ Thế hệ thứ ba đã giải đầy đủ sự biến đổi của mật độ năng lượng sóngtheo tần số và hướng sóng (chẳng hạn như mô hình WAM, WAVEWATCH, SWAN)…Nghiên cứu của Janssen (2006) được xem là thế hệ thứ tư, trong đó có bao gồmphương trình biến đổi góc pha [20] Mô hình sóng theo cách tiếp cận này phù hợp đểtính toán với các miền lớn bước tính rộng Cách tiếp cận khác là giải trực tiếp phươngtrình sóng phi tuyến ở vùng nước nông và vùng chuyển tiếp đó như các mô hìnhBoussinesq [40] [39], đây là mô hình liên tục pha theo thời gian nên cho kết quả tính

sóng khá tốt ở ven bờ Tuy nhiên mô hình không tính đến sự hình thành và phát triểnsóng từ năng lượng gió, mô hình này cần độ phân giải cao nên phù hợp cho vùng nhỏ,thời gian mô phỏng ngắn.

Mô hình cân bằng năng lượng sóng được sử dụng khá phổ biến trong các mô hình tínhvận chuyển bùn cát và xói bồi ven biển bởi mô hình có tính hiệu năng cao và độ tincậy phù hợp Các phương trình cơ bản tổng quát mô tả sự suy giảm chiều cao sóng cótính đến mực nước ngang bờ do sóng  (trung bình trong thời gian) gồm phương trìnhphương trình cân bằng năng lượng và động lượng ngang bờ như sau [1]:

d(cgE) + Dw + Df = 0

H rms: Chiều cao sóng quân phương; Dw: Tiêu tán năng lượng do sóng vỡ;

Df: Tiêu tán năng lượng do ma sát đáy; d:Độ sâu nước trung bình cục bộ.Young (1999) cho rằng trong số các quá trình vật lý được mô hình mô phỏng truyềnsóng vào bờ biển thì hiệu ứng nước nông và sóng vỡ do độ sâu nước hạn chế là quátrình chủ đạo; ma sát đáy, khúc xạ và tương tác sóng là các quá trình quan trọng;những quá trình khác là những quá trình ít quan trọng hơn [56] Jinhai (2008) cũng chorằng ở vùng nước nông, các quá trình sóng vỡ là quá trình giảm năng lượng chính sovới tất cả các quá trình khác, chi phối sự biến đổi theo không gian của chiều cao sóng[54]

Những tài liệu thu thập được cho thấy yếu tố khó khăn nhất trong việc mô hình hóaquá trình truyền sóng ven bờ là việc mô phỏng quá trình tiêu tán năng lượng sóng dohiện tượng sóng vỡ Điều này là do sự phức tạp của hiện tượng và hiếm có lời giải giảitích cho các trường sóng thực tế cũng như chưa hiểu biết rõ về hiện tượng này Kế đến

là sự tương tác của sóng với công trình cũng ảnh hưởng đáng kể đến chiều cao sóngtới, đặc biệt là sóng phản xạ Bất kể công trình ven biển nào cũng tạo ra sóng phản xạ,chính sóng phản xạ này sẽ làm ảnh hưởng đến sóng tới sóng leo và quá trình sóng vỡ

[58] Các mục dưới đây sẽ trình bày chi tiết hơn về các nghiên cứu tính toán sóng vỡ

và sóng phản xạ

Telewac (2017), Roelvink (2015) đã thống kê một số công thức tính toán sự tiêu tánnăng lượng Trong đó, công thức của Battjes (1978) được sử dụng rộng rãi để tính toán

sự suy giảm năng lượng sóng khi sóng đi qua vùng sóng vỡ trong một số mô hình phầnmềm thông dụng như Mike21SW, Litprof, Swan, Tomawac, Xbeach, Unibest-TC,Durosta, Wadibe… [60] Các công thức này đều được xây dựng dựa trên các thínghiệm trong mô hình vật lý, các hệ số được xác định dựa trên việc kiểm định ứng với

số các kịch bản thí nghiệm nhất định

Các nghiên cứu của Hedegaard (1992), Battjes và Jansen (1985, 2008), Salmon (2016)cũng đã tiến hành đánh giá các công thức sóng vỡ khác nhau trong các mô hình sóng[61] [48] [62] Kết quả cho thấy các công thức ứng với các tham số sóng vỡ khác nhau

có thể tạo ra sự khác biệt lớn về kết quả tính chiều cao sóng trong vùng sóng vỡ Từ sựkhác nhau đó có thể thấy việc lựa chọn công thức tính sóng vỡ phụ thuộc rất nhiều vào

số liệu quan trắc để hiệu chỉnh và kiểm định mô hình cho khu vực nghiên cứu

Trong vùng nước nông, chiều cao sóng tối đa của mỗi con sóng (Hmax) trong vùngsóng vỡ thường được biểu thị bằng một phần của độ sâu cục bộ: Hmax = γd, trong đó γđược gọi là hệ số sóng vỡ, d là độ sâu nước, độ lớn thực tế của nó thể thay đổi theo độdốc bãi biển và độ dốc của sóng tới Battjes (1974) cho rằng các giá trị của γ thườngnằm trong khoảng từ 0,7 đến 1,2 và chỉ số Iribarren  có thể được sử dụng để chỉ racác dạng sóng vỡ: sóng vỡ cuộn, vỡ đổ, sóng vỡ tràn, sóng vỡ dâng [64] Các dạngsóng vỡ này có thể ảnh hưởng đến quá trình vận chuyển bùn cát Andersen (2011) đưa

ra giới hạn của γs là từ 0,6 đến 0,8 đối với sóng tuyến tính và 0,5 đối với sóng phituyến (chiều cao sóng có ý nghĩa lớn nhất Hs) [13] Kris (2014) chỉ ra rằng hệ số nàydao động từ 0,3 - 0,6 và tăng theo độ dốc của bãi biển [65] Asbury (1985) cho thấy tỷ

lệ giữa chiều cao sóng quân phương với độ sâu nước địa phương (γrms) dao động trongkhoảng 0,29 đến 0,55 [66]

Hệ số sóng vỡ  là thông số quan trọng nhất trong các mô hình sóng bởi nó kiểm soátkhả năng sóng vỡ và mức năng lượng vùng nước nông nên cần thiết phải nghiên cứuhiệu chỉnh, kiểm định và lựa chọn thông số này Battjes (1985) dựa trên các số liệu thínghiệm và đo đạc hiện trường để kiểm định công thức Battjes và Jansen (1978) đã chỉ

ra rằng quá trình sóng vỡ ở vùng nước nông bị ảnh hưởng bởi độ dốc sóng ở vùngnước

sâu Sau đó, công thức Battjes &Janssen (2008) được sử dụng chuỗi số liệu khác cholại cho hệ số trong công thức [62] Ruessink (2003) tính toán với số liệu thực đo ở bãibiển Duck, Bắc Carolina được hệ số sóng vỡ trên mặt cắt ngang phụ thuộc và độ sâunước cục bộ [67] Nelson (1987) đưa ra công thức tính hệ số sóng vỡ là một hàm của

độ dốc đáy cục bộ Alex Apotsos (2008) đánh giá các mô hình trước đó và đưa ra côngthức tính  phụ thuộc vào chiều cao sóng nước sâu [68]và các hệ số được hiệu chỉnhtheo một bộ số liệu Salmon (2016) đã đề xuất một tỉ lệ kết hợp phụ thuộc vào cả độdốc đáy cục bộ và số sóng đặc trưng chuẩn hóa cục bộ (�� – với k là số sóng và d là)

độ sâu nước

Từ các nghiên cứu trên đây có thể thấy: quá trình sóng vỡ và tiêu tán năng lượng dosóng vỡ là quan trọng trong việc truyền động lượng từ sóng tới dòng chảy và là vấn đềphức tạp và vẫn cần phải tiếp tục nghiên cứu [62] Hầu hết các công thức tính toán đềuđược kiểm định bởi các số liệu có được từ số liệu đo đạc sóng trong phòng thí nghiệm(máng kín) Trên thực tế, khi áp dụng mô hình tính toán, việc kiểm định sự biến đổicủa chiều cao sóng hay phân bố của năng lượng trên mặt cắt ngang bờ trong vùng sóng

vỡ là cần thiết để đánh giá mô hình cũng như dự tính tin cậy dòng chảy trở lại bờ biểnphía dưới chân sóng

1.3.2 Mô hình vận chuyển bùn cát và hình thái

Với giả thiết rằng vận tốc chuyển động của hạt bùn cát tương ứng với vận tốc của hạtchất lỏng, phương trình cơ bản tính toán tốc độ vận chuyển bùn cát ngang bờ trên mặtthẳng đứng trên một đơn vị bề rộng được tính theo công thức [1]:

có gì biết được đối với sóng phi tuyến dọc theo đường bờ biển Hơn nữa trong thời kìbão ở vùng sóng vỡ, do dòng chảy đối ngược với dòng phản hồi, ảnh hưởng của sóngbất đối xứng đối với nồng độ bùn cát sát đáy lại càng khó hiểu biết một cách rõ ràng

Do đó cần phải

đơn giản hóa và sử dụng nồng độ trung bình theo thời gian và vận tốc của chất lỏng Tổng lượng vận chuyển bùn cát ngang bờ có thể được viết thành:

C(z) : nồng độ bùn cát trung bình theo thời gian

u(x, z ,t): vận tốc ngang bờ theo phương x [m/s]

u(z) : vận tốc trung bình theo thời gian tại một điểm [m/s]

z: vị trí trên phương ngang z=zb tại đáy [m]

Việc tính toán vận chuyển bùn cát theo cách tiếp cận đầu tiên ở trên thường có nhữngkhó khăn lớn khi định lượng hàm u và C Những hiểu biết về phản ứng của C dưới tácđộng của sóng đều rất hạn chế và gần như không có gì biết được đối với sóng phituyến ở khu vực ven biển Hơn nữa trong thời kì bão ở vùng sóng vỡ, do dòng chảyđối ngược với dòng phản hồi, ảnh hưởng của sóng bất đối xứng đối với nồng độ bùncát sát đáy lại

càng khó hiểu biết một cách rõ ràng Do đó vẫn cần phải đơn giản hóa nồng độ bùn cát

và vận tốc của chất lỏng, dẫn đến rất nhiều nghiên cứu với nhiều công thức tính khácnhau Mục 1.3 đã trình bày về u, phần dưới đây trình bày các nghiên cứu tính toánnồng độ (C)

Phương pháp phổ biến nhất để xác định nồng độ bùn cát (C) là dựa trên phương trìnhkhuếch tán như sau [19] [74]:

= ωs(z) + (ϑs(z) )

Trong đó ��là vận tốc lắng chìm của hạt cát và ��là hệ số khuếch tán

Trắc diện theo phương thẳng đứng của nồng độ bùn cát lơ lửng trên nền cát có thểđược mô tả bằng công thức: C(z) = C0fc(z), với C0 nồng độ bùn cát tham chiếu tại đáy;

fc(z) là hàm phân bố bùn cát không thứ nguyên và phụ thuộc vào phân bố của hệ sốxáo trộn bùn cát �� Cũng giống như phân bố hệ

số nhớt rối, hệ số xáo trộn bùn cát ��

được giả thiết với nhiều dạng phân bố

khác nhau như [1] đã trích dẫn tổng

quan nghiên cứu của Sistermans (2002)

trình bày trên Hình 1-3 Một số công

thức tính toán phân bố hệ số xáo trộn và

bùn cát đã được các tác giả đánh giá như Van Rijn (1993), Amiri (2009) [76] Trong số

đó, công

thức của Zanke (1977) được đánh giá là khá phổ biến đối với bờ biển cát có đườngkính hạt từ 0,1 đến 1mm.

Nồng độ tham chiếu tại đáy C0 thường được xác định tại biên trên cùng của vậnchuyển bùn cát đáy (z= δ) hay có thể nói đó là nồng độ bùn cát lơ lửng gần đáy Cónhiều công thức tính C0 được phát triển [52]: Theo Einstein (1950), C 0 được xác địnhtheo lưu lượng vận chuyển vận chuyển bùn cát đơn vị và vận tốc ma sát đáy Van Rijn(1984b) đã đưa ra mực tham chiếu (δ) tại chiều cao độ nhám tương đương hay một nửachiều cao hình dạng đáy và nồng độ tham chiếu phụ thuộc vào đường kính hạt trungbình, thông số kích thước hạt phi thứ nguyên và thông số ứng suất tiếp đáy không thứnguyên Zyserman và Fredsøe (1994) đề nghị tại mực tham chiếu tại z = 2d50, (d50 làđường kính hạt trung bình) và công thức tính nồng độ bùn cát tham

chiếu phụ thuộc tham số Shields và vận tốc

ma sát đáy Steetzel (1993) đưa ra công thức

tính nồng độ �0 phụ thuộc vào ứng suất tiếp

tới hạn, sự tiêu tán động năng, hàm không

thứ

nguyên liên quan đến đường kính hạt và đặc

biệt là sự có mặt của hàm mô tả ảnh hưởng

(1 − )� = −

Từ các công thức trên cho thấy có rất nhiều yếu tố ảnh hưởng đến vận chuyển bùn cát

và tiến thoái của mặt cắt địa hình Nhóm các yếu tố quan trọng đầu tiên liên quan đếnđặc điểm thủy động lực học gồm chiều cao, chu kì sóng, dòng chảy (dòng phản hồi),

độ sâu nước, độ nhớt, độ nhám đáy Nhóm yếu tố thứ hai đó là đặc trưng về bùn cátnhư đường kính, hình dạng, kích thước hạt, mật độ hạt Nhóm cuối cùng là đặc điểmhình học của cao độ đáy như độ dốc đáy Việc áp dụng các công thức khác nhau trongcác mô hình biến đổi mặt cắt ngang sẽ cho các kết quả khác nhau khá lớn Ví dụ điểnhình đối với sự khác nhau này đã được nghiên cứu bởi Hedegaard (1992) [61], sự sosánh kết quả của 6 mô hình được phát triển bởi các cơ quan khác nhau (1)DanishHydraulic Institute, (2)Delft Hydraulics, (3) HR Wallingford Ltd, (4)LaboratoireNational d'Hydraulique, (5) University of Liverpool và (6) Sogreah

Như vậy việc mô hình hóa tất cả các quá trình liên quan đến sự biến đổi hình thái mặtcắt ngang nên dựa trên sự sẵn có của số liệu tại vùng quan tâm để lựa chọn các phân

bố nhớt rồi và hệ số xáo trộn theo phương thẳng đứng đối với dòng chảy và nồng độbùn cát trong mô hình để có kết quả hợp lý Các số liệu đo đạc vận chuyển bùn cát vàxói lở ngang bờ ở ngoài thực tế là rất khó khăn trong điều kiện bão, chính vì thế đaphần số liệu của các nghiên cứu dựa trên các thí nghiệm mô hình vật lý để đánh giá vàhiệu chỉnh, kiểm định mô hình toán Các số liệu nhất định từ các mô hình vật lý chủyếu thí nghiệm với công trình tường đứng ứng với số kịch bản hạn chế nhất định [70][71] [72] [73]… Mặt khác các mô hình mô phỏng đều coi các biên công trình đê biểnnhư là biên cứng chưa xem xét đến ảnh hưởng của kết cấu mái đê có tường đỉnh đếndòng chảy và diễn biến chiều sâu hố xói như thế nào

1.4 Tổng quan nghiên cứu về sóng phản xạ

Hiện tượng phản xạ sóng được xem là có ảnh hưởng đến đặc trưng thủy động lực học

do sóng đến và các hoàn lưu ven bờ, làm gia tăng xói lở chân công trình Để đặc trưngcho hiện tượng này, các nhà khoa học sử dụng hệ số phản xạ Kr, được xác định là tỉ lệgiữa chiều cao sóng phản xạ và chiều cao sóng đến hoặc tỉ lệ giữa năng lượng/thônglượng sóng phản xạ và năng lượng/thông lượng sóng tới [12] [79] [80] Hệ số sóngphản xạ phụ thuộc vào đặc điểm bãi dốc và hình dạng kết cấu công trình trên đó, đốivới bãi