Bài báo trình bày kết quả tính toán phân bố dòng chảy dọc bờ do sóng đổ nhào phía ngoài đầm Ô Loan thuộc vùng biển Tuy An, Phú Yên ứng với các đợt gió mùa điển hình cấp 6 (Vnk = 13 m/s, hướng Đông Bắc (NE), Đông (E) và Đông Nam (SE).
Tuyển Tập Nghiên Cứu Biển, 2013, tập 19: 17-26 PHÂN BỐ DỊNG CHẢY DỌC BỜ DO SĨNG ĐỔ NHÀO TRONG CÁC TRƯỜNG GIĨ ĐIỂN HÌNH TẠI VÙNG BIỂN TUY AN, PHÚ YÊN Đỗ Như Kiều1, Lê Đình Mầu2 Trường Đại học Khoa học Tự nhiên, Đại học Quốc gia TP HCM Viện Hải dương học, Viện Hàn lâm Khoa học & Cơng nghệ Việt Nam Tóm tắt Bài báo trình bày kết tính tốn phân bố dịng chảy dọc bờ sóng đổ nhào phía ngồi đầm Ô Loan thuộc vùng biển Tuy An, Phú Yên ứng với đợt gió mùa điển hình cấp (Vnk = 13 m/s, hướng Đông Bắc (NE), Đông (E) Đơng Nam (SE) So sánh độ xác cơng thức thực nghiệm tính tốc độ dịng chảy dọc bờ sóng đổ nhào theo qui phạm bảo vệ bờ biển Hải quân Mỹ (SPM, 1984) Rattanapitikon – Shibayama (2006) với số liệu thực đo Kết tính tốn cho thấy cơng thức Rattanapitikon – Shibayama (2006) phù hợp Với sóng hướng NE: dịng có hướng từ bắc xuống nam với vận tốc đạt giá trị cực đại 0,7 m/s cửa An Hải vị trí cách cửa An Hải khoảng km phía nam Với sóng hướng E, nhìn chung dịng có hướng từ nam lên bắc, xuất điểm hội tụ dịng có vận tốc nhỏ 0,1 m/s Vận tốc dòng đạt giá trị cực đại 0,91 m/s cửa An Hải Tại cửa Lễ Thịnh dịng có vận tốc từ 0,25 đến 0,5 m/s Với sóng hướng SE, dịng có hướng từ nam lên bắc vận tốc dịng lớn so với trường hợp sóng hướng NE E Vận tốc dòng cực đại 1,02 m/s cửa An Hải Khu vực hội tụ sóng có vận tốc từ 0,5 đến 0,7 m/s DISTRIBUTION OF LONGSHORE CURRENT INDUCED BY BREAKING WAVE CORRESPONDING TO TYPICAL WIND CONDITIONS ALONG TUY AN COAST, PHU YEN PROVINCE Do Nhu Kieu1, Le Dinh Mau2 University of Science, Vietnam National University – Ho Chi Minh City Institute of Oceanography, Vietnam Academy of Science & Technology Abstract This paper presents the calculated results of wave and longshore current induced by breaking wave along Tuy An coast, Phu Yen province corresponding to typical wind conditions (Vnk = 13 m/s in NE, E and SE directions) Comparing the accuracy degree between two experimental formulae of SPM (1984) and Rattanapitikon–Shibayama (2006) for calculation of longshore current induced by breaking wave and the practical data Study results show that Rattanapitikon–Shibayama formula is more suitable for calculation of longshore current induced by breaking wave in Tuy An area For NE wave direction, the current flows from north to south The maximum current velocity is 0.7 m/s at An Hai inlet and the position from An Hai inlet about km towards south For E wave direction, the current flows from south to north and appear the convergence points with 17 small velocity, less than 0.1 m/s The maximum current velocity is 0.91 m/s at An Hai inlet The current velocity is from 0.25 to 0.5 m/s around Le Thinh inlet For SE wave direction, the current flows from south to north The current velocity attains the highest value in three calculated cases The maximum current velocity is 1.02 m/s at An Hai inlet The velocity in the area where occurs convergent phenomenon is from 0.5 to 0.7 m/s I MỞ ĐẦU Shibayama (2006) Các công thức thể Bảng Ngoài ra, số mơ hình động lực khác tích hợp thêm module tính dịng dọc bờ sóng đổ nhào như: mơ hình hồn lưu ven bờ SYMPHONIE, mơ hình động lực ven bờ MORPHODYN tính tốn dựa phương trình cân lượng phương trình liên tục Việc nghiên cứu chế hình thành tính tốn phân bố dịng chảy dọc bờ sóng đổ nhào thực nhiều tác giả khác Đồng thời, tác giả đưa số công thức khác cho việc tính tốn hầu hết cơng thức thực nghiệm như: Longuet– Higgins (1970), Komar Inman (1970), Galvin (1987), Rattanapitikon Bảng Một số công thức thực nghiệm tính dịng dọc bờ sóng đổ nhào Table Some experimental formulae for calculation of longshore current induced by breaking wave Năm Tác giả 1970 Longuet – Higgins Công thức thực nghiệm vL = 5π ς S Cf Ghi ghb sin αb cos αb Hb: độ cao sóng vỡ (m) Hbs: độ cao sóng vỡ có vL = 2.7um sin α b cos α b 1970 nghĩa (m) um = ghb = gH b / γ Komar Inman S: độ dốc vL = 1.0 gH bs sin α b cos α b 1984 Longuet –Higgins (trong SPM, 1984) vL = 20.7 S gH b sin(2α b ) 1987 Galvin vL = gST sin 2α b 1999 Pilirczyk (trong Kamran cs., 2012) vL = K L gH b sin 2α b 2006 Rattanapitikon – Shibayama K L = 0.3 ÷ 0.6 vL = 2.7uˆb sin α b cos α b uˆb = ( −0.57S ) + 0.31S + 0.58 π Cb H ( kb hb ) L0 0.83 hb: độ sâu sóng vỡ (m) γ : tiêu sóng vỡ Cf: hệ số ma sát kéo α b : góc sóng vỡ (độ) H0: độ cao sóng nước sâu (m) L0: chiều dài sóng nước sâu (m) kb: số sóng vùng sóng vỡ T: chu kì sóng (s) thực đo như: “Đặc điểm phân bố đặc trưng sóng vịnh Nha Trang trường gió mùa điển hình” (Lê Đình Mầu cs., 2010) “Đặc điểm phân bố đặc trưng sóng vùng biển Lagi (Bình Thuận) tác động chúng đến q trình xói lở – bồi tụ” (Lê Đình Mầu, 2010) Trong báo này, phân bố dịng chảy dọc bờ sóng Tại Việt Nam, có nghiên cứu tính tốn dịng chảy dọc bờ sóng đổ nhào Các tính tốn hầu hết dựa nghiên cứu, giả thiết Longuet–Higgins (trong SPM, 1984) phần lớn mang tính chất định tính hướng, riêng độ lớn vận tốc dịng có khác biệt tương đối lớn so với số liệu 18 trình có liên quan Thường mơ hình Longuet–Higgins phù hợp cho vùng biển thống, điều kiện sóng mạnh (bên bờ đại dương) nên cho tốc độ dịng lớn thực tế áp dụng Việt Nam Trong nghiên cứu này, hai mơ hình tiến hành tính tốn so sánh với số liệu thực đo Từ đó, chọn cơng thức phù hợp để tính tốn cho khu vực nghiên cứu Phạm vi tính tốn thể Hình Địa hình khu vực nghiên cứu vùng biển nơng ven bờ, có địa hình đáy phức tạp Độ dốc giảm dần từ bắc xuống nam Tuy nhiên, đường bờ có dạng vịng cung che chắn mỏm đá phía bắc cửa Lễ Thịnh đảo Mái Nhà nên lượng sóng truyền vào vùng ven bờ phân bố khơng đồng đều, dẫn đến phân bố không đồng vận tốc dòng chảy dọc bờ sinh sóng đổ nhào Vị trí khu vực tập trung lượng sóng tác động vào bờ phụ thuộc vào hướng sóng tới Do khn khổ báo, chúng tơi trình bày phân bố dịng chảy dọc bờ sóng đổ nhào theo hướng sóng tác động chính: Đơng Bắc (NE), Đông (E) Đông Nam (SE) đổ nhào điều kiện gió mùa điển hình tính tốn theo quy trình sau: - Xác định đặc trưng sóng ngồi khơi vùng biển Tuy An từ đặc trưng gió thống kê trạm Tuy Hịa: 1988 – 2007 (Hình 2) mơ hình Dolphin (Mandal Holthuijsen, 1985) - Xác định đặc trưng sóng ven bờ mơ hình SWAN (Holthuijsen cs., 2003) với điều kiện biên sóng ngồi khơi, gió địa phương phân bố độ sâu vùng nghiên cứu - Tính vận tốc dịng dọc bờ sóng đổ nhào mơ hình thực nghiệm: Longuet–Higgins (trong SPM, 1984) Rattanapitikon–Shibayama (2006) Trong đó, mơ hình Longuet–Higgins mơ hình đơn giản nhất, phát triển sớm cơng thức chuẩn để tính tốc độ dịng dọc bờ sóng đổ nhào Hải quân Mỹ Tuy nhiên, hệ số thực nghiệm dao động lớn, điều yêu cầu phải có nhiều số liệu đo đạc thực địa đồng nhằm xác định xác hệ số thực nghiệm Mơ hình Rattanapitikon–Shibayama hơn, đại hơn, có xét đến nhiều Hình Đặc điểm khu vực nghiên cứu Figure Feature of study area 19 Hình Hoa gió trạm Tuy Hòa (1988 – 2007) Figure Wind rose diagram at Tuy Hoa station (1988 – 2007) II TÀI LIỆU VÀ PHƯƠNG PHÁP triển bền vững kinh tế xã hội, bảo vệ môi trường dải ven biển Nam Trung Bộ (Đà Nẵng – Bình Thuận)” Viện Hải dương học chủ trì (2010-2014) Tài liệu Số liệu gió thu thập trạm Tuy Hịa, Phú n (1988 – 2007) phân tích, thống kê từ số liệu lưu trữ đề tài độc lập cấp Viện Hàn lâm Khoa học Công nghệ Việt Nam (VAST): “Đánh giá tác động trường sóng gió mùa đến dải ven biển Nam Trung Bộ từ Phú Yên đến Bình Thuận đề xuất giải pháp giảm nhẹ thiệt hại phục vụ phát triển bền vững” (2007-2009) Độ sâu khu vực nghiên cứu lấy từ hải đồ tỉ lệ 1/50.000 xuất năm 2005 Bộ Tài nguyên Môi trường Các đặc trưng sóng đo máy ALEC AWH – 16M điểm 109017’54”E, 13017’28” N (bờ phía nam cửa An Hải) sử dụng để kiểm chứng mơ hình SWAN, dịng chảy đo phương pháp thả phao trơi gió đo máy đo gió cầm tay ANEMOMETER – AVM 01 Các liệu lấy từ đề tài bảo vệ mơi trường: “Áp dụng mơ hình đại nhằm đánh giá, phòng tránh giảm thiểu thiệt hại tác động mơi trường tượng đóng/mở cửa sông, đầm phá phục vụ chiến lược phát Phương pháp 2.1 Tính tốn đặc trưng sóng biển sâu Mơ hình số trị tính sóng biển khơi Dolphin (Mandal Holthuijsen, 1985) tính tốn kết hợp sóng gió sóng lừng Phương trình cân lượng phổ sóng hai chiều: - Trường hợp sóng gió hướng θ: ∞ dE1 (θ ) = S1 (θ ) + ∫ S2 ( f ,θ )df (1) dt Trong đó: E1(θ) mật độ lượng phổ sóng gió theo hướng θ S1(θ): tốc độ biến đổi E1(θ) S2(f, θ) tốc độ chuyển giao lượng sóng lừng chiều E2 sang sóng gió có tần số f - Trường hợp sóng lừng từ hướng θ: dE2 = − BS11 ( f ,θ ) − S ( f ,θ ) − S3 ( f ,θ ) (2) dt S1 < B =1 S1 > B = 20 Các đặc trưng sóng đổ nhào (độ cao, hướng, vị trí) lấy từ kết tính sóng mơ hình SWAN theo tiêu chuẩn γ b = H b hb với γ b lấy ≥ 0.6 độ cao sóng sử dụng để tính tốn độ cao sóng có nghĩa sóng ngẫu nhiên (SPM, 1984) Vận tốc dịng chảy dọc bờ sóng đổ nhào gây tính theo cơng thức Longuet–Higgins (trong SPM, 1984) Rattanapitikon – Shibayama (2006) So sánh kết tính tốn với số liệu đo đạc Chọn cơng thức tính phù hợp cho khu vực nghiên cứu + Longuet – Higgins (trong SPM, 1984): Trong đó: E2(f,θ) mật độ lượng phổ sóng lừng với tần số f theo hướngθ S11 tốc độ biến đổi E2 S2 tốc độ chuyển giao lượng sóng lừng E2 sang sóng gió có tần số S3 tốc độ biến đổi lượng sóng lừng gây hiệu ứng nước nông Số liệu đầu vào mơ hình: hướng tốc độ gió, thời gian gió tác động chọn từ việc thống kê số liệu gió nhiều năm Số liệu đầu mơ hình: độ cao sóng có nghĩa Hs, chu kỳ đỉnh phổ sóng Tp, độ phân tán phổ lượng DSPR hướng sóng θ 2.2 Tính tốn đặc trưng sóng nước nơng ven bờ Các đặc trưng sóng vùng biển ven bờ xác định mơ hình SWAN (Holthuijsen cs., 2003) sở phương trình cân tác động phổ ∂ ∂ ∂ ∂ ∂ S N + Cx N + C y N + Cσ N + Cθ N = (3) ∂t ∂x ∂y ∂σ ∂θ vL = 20.7 S gH b sin(2α b ) (4) + Rattanapitikon – Shibayama (2006): vL = 2.7uˆb sin α b cos α b (5) Với: uˆb ( −0.57S = ) + 0.31S + 0.58 π Cb H ( kb hb ) L0 0.83 (6) Trong đó: vL vận tốc dịng chảy dọc bờ sóng đổ nhào gây (m/s) ûb vận tốc truyền sóng vị trí sóng vỡ S độ dốc bãi biển kb số sóng vị trí sóng vỡ hb độ sâu vị trí sóng vỡ H0, L0 độ cao sóng chiều dài sóng nước sâu αb hướng sóng vị trí sóng vỡ σ Trong đó: N mật độ phổ tác động (action density); Cx, Cy tốc độ lan truyền N theo không gian x y; Cσ, Cθ tốc độ biến đổi N theo tần số tương đối σ hướng θ; S = S(σ, θ) hàm nguồn Mơ hình SWAN tính hiệu ứng nước nơng, khúc xạ, tạo sóng gió địa phương, sóng bạc đầu, tương tác sóng – sóng, ma sát đáy, sóng đổ nhào,…Phạm vi áp dụng mơ hình khu vực vùng biển bên ngồi đầm Ơ Loan với diện tích khoảng 56 km2 Kích thước lưới 48 m x 48 m Số liệu đầu vào mô hình: - Các tham số sóng ngồi khơi: Hs, Tp, DSPR θ lấy từ kết mô hình Dolphin - Trường gió địa phương ổn định theo thời gian (V = 10 m/s, hướng NE, E SE) - Phân bố độ sâu vùng nghiên cứu III KẾT QUẢ Các đặc trưng sóng ngồi khơi vùng biển Tuy An Kết tính tốn từ mơ hình Dolphin cho đặc trưng sóng ngồi khơi vùng biển nghiên cứu ứng với trường gió ngồi khơi điển hình Vnk = 13 m/s theo hướng NE, E SE trình bày Bảng 2 Các đặc trưng sóng ven bờ dịng chảy dọc bờ sóng đổ nhào ứng với trường sóng tác động khu vực nghiên cứu 2.1 So sánh kết tính tốn dịng chảy dọc bờ sóng đổ nhào theo Longuet – Higgins (trong SPM, 1984) Rattanapitikon – Shibayama (2006) Tính tốn dịng chảy dọc bờ sóng đổ nhào gây 21 Kết so sánh cho thấy giá trị vận tốc dòng dọc bờ sóng đổ nhào tính từ cơng thức Longuet–Higgins (trong SPM, 1984) lớn so với kết khảo sát thực tế sơ Ngược lại, kết tính tốn sử dụng cơng thức Rattanapitikon– Shibayama (2006) lại cho giá trị vận tốc dòng nhỏ giá trị thực đo với sai số nhỏ nhiều so với công thức Longuet–Higgins (trong SPM, 1984) Do đó, cơng thức Rattanapitikon–Shibayama (2006) chọn để tính tốn dịng chảy dọc bờ sóng đổ nhào khu vực nghiên cứu Những liệu tính tốn, so sánh trình bày Bảng Bảng Các đặc trưng sóng ngồi khơi vùng biển Tuy An Table Offshore wave characteristics in Tuy An waters Các đặc trưng sóng vùng khơi Độ cao sóng có nghĩa (Hs) Chu kỳ sóng đỉnh phổ (Tp) Hệ số phân tán phổ sóng (DSPR) Độ lớn 2,33 m 6,95s 27,080 Bảng So sánh vận tốc dòng dọc bờ sóng đổ nhào đo đạc với kết tính toán theo Longuet– Higgins (trong SPM, 1984) Rattanapitikon – Shibayama (2006) Table Comparison between calculated results of longshore current induced by breaking wave corresponding to Longuet–Higgins (in SPM, 1984) Rattanapitikon–Shibayama (2006) Tính tốn Tọa độ Thực đo Longuet Rattanapitikon – Higgins – Shibayama Điểm Hướng x y Vận tốc (m/s) Hướng Vận tốc (m/s) Vận tốc (m/s) 109,295 13,2916 0,2 Bắc xuống Nam 0,41 0,17 Bắc xuống Nam 109,2929 13,2961 0,26 Nam lên Bắc 1,03 0,22 Nam lên Bắc 109,2896 13,3023 0,22 Nam lên Bắc 0,42 0,18 Nam lên Bắc 109,2964 13,2902 0,28 Bắc xuống Nam 0,44 0,24 Bắc xuống Nam 109,286 13,3342 0,18 0,03 Bắc xuống Nam Hình So sánh vận tốc dịng dọc bờ sóng đổ nhào thực đo với kết tính tốn theo Longuet – Higgins (trong SPM, 1984) Rattanapitikon – Shibayama (2006) Figure Comparison between calculated results of longshore current induced by breaking wave corresponding to Longuet –Higgins (in SPM, 1984) Rattanapitikon – Shibayama (2006) 22 2.2 Dòng chảy dọc bờ sóng đổ nhào + Trường hợp sóng hướng Đơng Bắc (NE): Kết mơ trường sóng cho thấy, phần hứng sóng đảo Mái Nhà mỏm đá phía bắc cửa Lễ Thịnh khu vực chịu tác động mạnh sóng Tuy nhiên, đoạn bờ từ cửa An Hải dài khoảng 1,5 km phía nam khu vực đáng ý với đường đẳng độ cao sóng m tiến sát vào bờ Đoạn bờ đến mỏm đá phía nam từ cửa An Hải phía bắc gần km, sóng vào bờ có độ cao khoảng 1,5 m Đoạn bờ dài khoảng 1,5 km gần cửa Lễ Thịnh có độ cao sóng giảm dần từ 1,0 đến 0,5 m cửa (Hình 4a) Kết tính tốn dịng chảy dọc bờ sóng đổ nhào cho thấy, vận tốc phân bố khơng tồn khu vực tính tốn Các vị trí có độ dốc thay đổi nhiều thường có vận tốc dịng lớn Vùng bờ gần cửa Lễ Thịnh có vận tốc dịng nhỏ nằm vùng che chắn Vận tốc dịng trung bình tổng số điểm sóng vỡ khoảng 0,14 m/s Vận tốc cực đại vị trí phía bắc gần cửa An Hải vị trí phía nam cách cửa An Hải khoảng km với giá trị 0,7 m/s Tại cửa Lễ Thịnh vận tốc dòng đạt giá trị cực tiểu Các khu vực phía bắc cửa An Hải có vận tốc dịng từ 0,05 đến 0,5 m/s Nhìn chung, dịng có hướng từ bắc xuống nam (Hình 4b) (a) (b) Hình (a) Trường độ cao sóng hữu hiệu; (b): Trường dịng chảy dọc bờ sóng đổ nhào Gió ngồi khơi: Vnk = 13 m/s; Gió ven bờ: Vvb = 10 m/s; Hướng: NE Figure (a) Pattern of significant wave height; (b): Pattern of longshore current induced by breaking wave Offshore wind: Vnk = 13 m/s; Nearshore wind: Vvb = 10 m/s; Direction = NE + Trường hợp sóng hướng Đơng (E): Trong trường hợp này, gần toàn khu vực bờ vùng nghiên cứu chịu tác động mạnh sóng với độ cao sóng từ 1,5 đến 2,0 m trừ cửa Lễ Thịnh đoạn bờ 500 m gần mỏm đá phía nam có độ cao sóng vào bờ khoảng 1,0 m Vùng khuất sóng gần đảo Mái Nhà có độ cao sóng nằm khoảng 0,5 đến 1,0 m tăng dần xa đảo Vùng tập trung lượng sóng lớn có chiều dài khoảng 1,5 km, cách cửa An Hải khoảng 500 m phía bắc (Hình 5a) Nhìn chung, dịng có hướng từ nam lên bắc đặc trưng địa hình khu vực nghiên cứu có xu hướng dốc phía bắc vận tốc dịng phân bố khơng đồng Tuy nhiên, xuất số vị trí hội tụ dịng, điển hình khu vực gần cửa Lễ Thịnh 23 vị trí cách cửa An Hải khoảng km phía nam (được đánh dấu hình 5b) Vận tốc dịng có xu hướng giảm nhanh vị trí dịng hội tụ có giá trị nhỏ 0,1 m/s Vận tốc dòng mạnh đạt giá trị 0,91 m/s gần An Hải Tại khu vực cửa Lễ Thịnh độ lớn vận tốc dòng từ 0,25 đến 0,5 m/s Vận tốc dịng trung bình khoảng 0,142 m/s (Hình 5b) (a) (b) Hình (a) Trường độ cao sóng hữu hiệu; (b): Trường dịng chảy dọc bờ sóng đổ nhào Gió ngồi khơi: Vnk = 13 m/s; Gió ven bờ: Vvb = 10 m/s; Hướng: E Figure (a) Pattern of significant wave height; (b): Pattern of longshore current induced by breaking wave Offshore wind: Vnk = 13 m/s; Nearshore wind: Vvb = 10 m/s; Direction: E + Trường hợp sóng hướng Đơng Nam (SE): Sóng tác động vào bờ khu vực nghiên cứu trường hợp nhỏ nhiều so với trường hợp sóng hướng NE E Khu vực hứng sóng đảo Mái Nhà chịu tác động sóng với độ cao m Mỏm đá phía bắc cửa Lễ Thịnh khu vực từ mỏm đá phía nam phía nam chịu tác động sóng có độ cao 1,5 m Toàn khu vực đường bờ nghiên cứu chịu tác động sóng có độ cao khoảng 1,0 m Riêng vị trí phía bắc cửa An Hải, cách bờ khoảng 200 m xảy tượng hội tụ sóng Vùng khuất sóng gần đảo Mái Nhà có độ cao sóng khoảng 1,0 m tăng dần xa đảo (Hình 6a) Kết tính tốn dịng dọc bờ sóng đổ nhào cho thấy, khu vực xuất hội tụ sóng có vận tốc dịng từ 0,5 đến 0,7 m/s, lớn nhiều so với hầu hết khu vực khác Mặc dù, vận tốc dòng phân bố khơng đồng nhìn chung vùng bờ phía nam cửa An Hải che chắn đảo Mái Nhà mỏm đá phía nam, vận tốc dịng có xu hướng giảm dần phía cửa An Hải với giá trị từ 0,6 đến 0,3 m/s Vận tốc dòng cực đại lớn, đạt giá trị 1,02 m/s gần cửa An Hải (bờ phía bắc) Vận tốc trung bình trường hợp 0,213 m/s Tại cửa Lễ Thịnh, vận tốc dịng có giá trị nhỏ 0,3 m/s Vận tốc nhỏ khu vực sát bờ đoạn bờ uốn cong (Hình 6b) 24 (a) (b) Hình (a) Trường độ cao sóng hữu hiệu; (b): Trường dịng chảy dọc bờ sóng đổ nhào Gió ngồi khơi: Vnk = 13 m/s; Gió ven bờ: Vvb = 10 m/s; Hướng: SE Figure (a) Pattern of significant wave height; (b): Pattern of longshore current induced by breaking wave Offshore wind: Vnk = 13 m/s; Nearshore wind: Vvb = 10 m/s; Direction: SE IV THẢO LUẬN Với sóng hướng E: dịng có hướng từ nam lên bắc, nhiên xuất số điểm hội tụ dịng ảnh hưởng địa hình hướng sóng tác động Vận tốc đạt giá trị cực đại 0,91 m/s cửa An Hải Vận tốc dòng gần cửa Lễ Thịnh có giá trị từ 0,25 đến 0,5 m/s Tại điểm hội tụ, vận tốc đạt giá trị nhỏ 0,1 m/s Với sóng hướng SE: dịng có hướng từ nam lên bắc vận tốc dòng lớn so với trường hợp sóng hướng NE E địa hình đáy dốc phía bắc Vận tốc dòng cực đại 1,02 m/s cửa An Hải Khu vực hội tụ sóng có vận tốc dịng từ 0,5 đến 0,7 m/s Kết tính tốn hướng dịng chảy dọc bờ sóng đổ nhào cho thấy phù hợp với thực tế, thông số quan trọng phục vụ thiết kế công trình bảo vệ bờ biển, cửa sơng Tuy nhiên, số liệu đo đạc vận tốc dòng chảy để kiểm chứng kết tính tốn cịn hạn chế số lượng kỹ thuật đo đạc nên độ tin cậy chưa cao Do đó, kết tính tốn tốc độ dòng chảy Các vùng khác khu vực nghiên cứu chịu tác động sóng mang lượng khác phụ thuộc vào hướng sóng tới bị chi phối nhiều điều kiện địa lý hình dạng đường bờ địa hình đáy ven bờ Dịng dọc bờ sóng đổ nhào sinh phụ thuộc vào hướng sóng tác động độ dốc khu vực ven bờ đới sóng vỡ Nhìn chung, vận tốc dịng phân bố khơng tỉ lệ thuận với độ cao sóng tác động, trừ vị trí có hội tụ phân kì dịng Đới sóng vỡ mở rộng phía biển đến độ sâu từ đến m tùy thuộc vào đặc trưng sóng tới Với sóng hướng NE: dịng có hướng từ bắc xuống nam Vận tốc dòng nhỏ so với trường hợp sóng hướng E SE Vận tốc dòng đạt giá trị cực đại cửa An Hải vị trí cách cửa An Hải khoảng km phía nam Các khu vực phía bắc cửa An Hải có giá trị thấp thấp cửa Lễ Thịnh 25 dọc bờ sóng đổ nhào nghiên cứu có tính chất tham khảo applied scientific research, 2(5): 4556 4562 Komar, P.D., and Inman, D.L., 1970 Longshore sand transport on beaches J Geophysical Research, 75:5914–5927 Lê Đình Mầu, 2010 Đặc điểm phân bố đặc trưng sóng vùng biển cửa La Gi (Bình Thuận) tác động chúng đến q trình xói lở - bồi tụ Kỷ yếu Hội nghị khoa học kỷ niệm 35 năm Viện KH&CN VN (1975-2010), Tiểu ban KHCN Biển, Hà Nội: 10/2010, tr 211-216 Lê Đình Mầu, Nguyễn Văn Tuân, Phạm Thị Phương Thảo, 2010 Đặc điểm phân bố đặc trưng sóng vịnh Nha Trang trường gió mùa điển hình Tuyển tập Nghiên cứu biển NXB KH&KT, ISSN:1859-2120, XVII: - 17 Longuet-Higgins, 1970 Longshore currents generated by obliquely incident sea waves J Geoph Res., 75:6778-6789 Mandal S and L.H Holthuijsen, 1985 A numerical wave prediction model DOLPHIN: Theory and test results Report No 3-85, Delft University of Technology, Department of Civil Engineering, Group of Fluid Mechanics, p - 70 Rattanapitikon W., and Shibayama T., 2006 Breaking wave formulas for breaking depth and orbital to phase velocity ratio Coast Eng J., 48(4):395416 SPM, 1984 Shore Protection Manual U.S Army Coastal Engineering Research Centre, Department of the Army Corps of Engineers, Washington DC USA Vol Lời cám ơn: Tập thể tác giả chân thành cám ơn ban chủ nhiệm đề tài: “Áp dụng mơ hình đại nhằm đánh giá, phịng tránh giảm thiểu thiệt hại tác động môi trường tượng đóng/mở cửa sơng, đầm phá phục vụ chiến lược phát triển bền vững kinh tế xã hội, bảo vệ môi trường dải ven biển Nam Trung Bộ (Đà Nẵng – Bình Thuận)” tập thể cán phòng Vật lý biển, Viện Hải dương học tạo điều kiện việc thu thập số liệu lưu trữ khảo sát thực địa động viên khích lệ q trình hồn thiện báo TÀI LIỆU THAM KHẢO Chistos A., 2013 New modified equation of longshore current velocity at the breaking point (for mixed and gravel beaches) Journal of Coastal Development, 16(2): 121 - 134 Galvin C J., 1987 Vertical profile of littoral sand tracers from a distribution of waiting time In: Proceedings of Coastal Sediments '87 (New Orleans, Louisiana, ASCE), pp 436–451 Holthuijsen L H., N Booij, R C Ris, I J G Haagsma, A T M M Kieftenburg, E E Kriezi, and M Zijlema, 2003 SWAN cycle III version 40.20 User Manual Delft University of Technology, the Netherlands Kamran L., S Kheiri, A Karami, M.Torabi Azad, M Abrehdary, 2012 Field study of longshore current along the Anzali coast in Caspian Sea Journal of basic and 26 ... Bảng 2 Các đặc trưng sóng ven bờ dịng chảy dọc bờ sóng đổ nhào ứng với trường sóng tác động khu vực nghiên cứu 2.1 So sánh kết tính tốn dịng chảy dọc bờ sóng đổ nhào theo Longuet – Higgins (trong. .. Mầu, 2010) Trong báo này, phân bố dòng chảy dọc bờ sóng Tại Việt Nam, có nghiên cứu tính tốn dịng chảy dọc bờ sóng đổ nhào Các tính tốn hầu hết dựa nghiên cứu, giả thiết Longuet–Higgins (trong SPM,... hướng sóng tới Do khn khổ báo, chúng tơi trình bày phân bố dịng chảy dọc bờ sóng đổ nhào theo hướng sóng tác động chính: Đơng Bắc (NE), Đơng (E) Đông Nam (SE) đổ nhào điều kiện gió mùa điển hình