The scale of the central beach change processes in the west coast of Nha Trang bay

The central beach in west coast of Nha Trang bay has experienced erosion in recent years. The determination of hydrodynamic regime and causes of this beach fluctuation is still an open problem and is concerned by scientists and managers. This study shows the causes and scales of the shoreline change processes based on the results of shoreline change model and the analysis of remote sensing images of shoreline location.

Vietnam Journal of Marine Science and Technology; Vol 20, No 1; 2020: 25–38 DOI: https://doi.org/10.15625/1859-3097/20/1/15039 http://www.vjs.ac.vn/index.php/jmst The scale of the central beach change processes in the west coast of Nha Trang bay Vu Cong Huu*, Dinh Van Uu VNU University of Science, Hanoi, Vietnam * E-mail: vuconghuu80@gmail.com Received: January 2019; Accepted: July 2019 ©2020 Vietnam Academy of Science and Technology (VAST) Abstract The central beach in west coast of Nha Trang bay has experienced erosion in recent years The determination of hydrodynamic regime and causes of this beach fluctuation is still an open problem and is concerned by scientists and managers This study shows the causes and scales of the shoreline change processes based on the results of shoreline change model and the analysis of remote sensing images of shoreline location Keywords: Nha Trang beach, EBED wave model, shoreline change model, camera image Citation: Vu Cong Huu, Dinh Van Uu, 2020 The scale of the central beach change processes in the west coast of Nha Trang bay Vietnam Journal of Marine Science and Technology, 20(1), 25–38 25 Tạp chí Khoa học Cơng nghệ Biển, Tập 20, Số 1; 2020: 25–38 DOI: https://doi.org/10.15625/1859-3097/20/1/15039 http://www.vjs.ac.vn/index.php/jmst Quy mô nguyên nhân trình biến động bãi biển trung tâm thuộc bờ tây vịnh Nha Trang Vũ Công Hữu*, Đinh Văn Ưu Trường Đại học Khoa học Tự nhiên, Đại học Quốc gia Hà Nội, Hà Nội, Việt Nam * E-mail: vuconghuu80@gmail.com Nhận bài: 2-1-2019; Chấp nhận đăng: 8-7-2019 Tóm tắt Bãi biển trung tâm thuộc bờ tây vịnh Nha Trang xảy q trình xói năm gần Việc xác định chế độ thủy động lực nguyên nhân gây biến động bãi biển toán mở nhà khoa học, quản lý quan tâm Nghiên cứu cho thấy nguyên nhân quy mô q trình gây biến động bãi biển thơng qua kết mơ hình tốn biến động bãi biển kết phân tích ảnh viễn thám (camera giám sát) vị trí đường bờ Các kết có ý nghĩa tham khảo cho nghiên cứu tính tốn q trình ven bờ vịnh Nha Trang Từ khóa: Bãi biển Nha Trang, mơ hình sóng EBED, mơ hình biến đổi đường bờ, ảnh camera MỞ ĐẦU Xói lở bãi biến đổi đường bờ diễn nghiêm trọng nhiều quốc gia nước ta Với 3.200 km chiều dài bờ biển, chưa kể đảo, Việt Nam chịu ảnh hưởng nặng nề q trình xói lở bờ biển diễn nhiều nơi, chẳng hạn vùng biển Hải Hậu (Nam Định), đầm phá Tam Giang (Huế), Lý Hồ (Quảng Bình), bi biển bắc cửa Tùng (Quảng Trị), bờ biển Bạc Liêu Mũi Cà Mau, bờ biển Hội An–Cửa Đại, bãi biển Nha Trang,… Hình Vịnh Nha Trang (phải) bãi biển trung tâm (trái) 26 The scale of the central beach change processes Vịnh Nha Trang 29 vịnh đẹp giới, có chiều dài khoảng 16 km rộng khoảng 13 km Vịnh thơng với biển ngồi qua hai cửa: Cửa phía đơng bắc, cửa nhỏ phía đơng nam (hình 1, bên phải) [1] Bãi biển bị xói mùa gió Đơng Bắc (ảnh T03/2013) Bên cạnh mạnh du lịch, bãi biển trung tâm thuộc bờ tây vịnh Nha Trang (hình 1, bên trái) phải đối mặt với biến đổi bãi biển trầm trọng Cụ thể biến đổi đường bờ bãi biển bị thu hẹp (hình 2) [1–3] Bãi biển bồi mùa gió Tây Nam (ảnh T07/2013) Hình Biến động bãi biển trung tâm vịnh Nha Trang năm 2013 Cho đến nay, bãi biển tồn số vấn đề sau: a) Bãi biển chịu tác động sóng lớn mùa gió Đơng Bắc có bão [1– 3]; b) Bãi biển bị tác động bồi xói, biến đổi mạnh theo mùa [1–3]; c) Bãi biển hẹp, có độ dốc lớn sâu, gây bất lợi cho việc tắm biển vào thời kỳ mùa gió Đơng Bắc; d) Sự phát triển cồn ngầm cửa sông Cái tác động tới bãi tắm lân cận cửa [3]; c) Diễn biến bãi biển chịu chi phối sóng mực nước [2, 3] Các nguyên nhân gây biến động bãi biển chưa xác định rõ ràng Kết nghiên cứu [1] luận chứng cần thiết cấp bách việc cải tạo tôn tạo bãi tắm biển thành phố Nha Trang, đưa định hướng khoa học công nghệ dựa kết phân tích liệu khảo sát địa hình, liệu giám sát camera kết mô hình số Tiếp đó, nghiên cứu [4] khai thác đầy đủ số liệu thực đo, từ camera sử dụng kỹ thuật phân tích hàm trực giao để xác định biến đổi thành phần thu phương trình dự báo vị trí đường bờ Các kết cho thấy thành phần có mối quan hệ với diễn biến độ cao sóng nước sâu theo mùa Tuy nhiên, phương pháp không quy mô thời gian thành phần gây biến đổi đường bờ bao quát hết quy mô trình gây biến động bãi biển Nghiên cứu tính tốn đặc trưng trường sóng ven bờ dịng vận chuyển trầm tích dọc bờ [5] cho thấy dịng vận chuyển trầm tích có xu từ Bắc xuống Nam cho thấy tồn trình vận chuyển ngang bờ làm biến đổi bãi biển thiếu hụt nguồn trầm tích cung cấp từ phía cửa sơng Cái Phương pháp mơ hình tốn đề cập đến cơng trình nghiên cứu [2, 3] tính tốn đặc trưng trường sóng chi tiết diễn biến đường bờ kết hợp sóng mực nước Kết mơ hình cho thấy diễn biến vị trí đường bờ có quan hệ mật thiết với diện kiện gió mùa, bão nước dâng Nguyên nhân mối quan hệ sóng truyền từ ngồi khơi vào kết hợp với dao động mực nước làm biến động bãi biển Nha Trang Nghiên cứu trình bày kết hệ thống mơ hình số biến đổi đường bờ cơng cụ phân tích, xử lý liệu đường bờ Từ đó, ngun nhân quy mơ q trình gây biến động bãi biển HỆ THỐNG CÁC MƠ HÌNH VÀ CƠNG CỤ PHÂN TÍCH Mơ hình biến động đường bờ: Biến đổi đường bờ trình phức tạp với nhiều quy mô khác mối tương tác q trình vận chuyển trầm tích ngang bờ dọc bờ biển Nghiên cứu áp dụng phương pháp mô biến đổi đường bờ 27 Vu Cong Huu, Dinh Van Uu rút từ kết nghiên cứu thí nghiệm số phịng thí nghiệm biến đổi đường bờ với quy mơ lớn nhỏ Cả hai thí nghiệm quy mô lớn nhỏ cho thấy điều kiện tác động sóng mực nước bãi biển tồn xu cân theo quy mô thời gian xấp xỉ dạng hàm mũ Các mô số [6, 7] cho thấy biến đổi đường bờ mơ hình hóa theo phương trình có dạng: dy  t   k  yeq  t   y  t   dt (1) Với y(t) yeq(t) tương ứng vị trí đường bờ vị trí đường bờ cân thời điểm t; k hệ số biến đổi đường bờ so với đường bờ cân Phương trình (1) dạng cân cải tiến từ phương trình cân dạng kinh điển cho thấy biến đổi đường bờ tỷ lệ với mức độ cân đường bờ Sơ đồ sai phân bán Crank-Nicholson áp dụng để giải số phương trình (1) Từ đó, thu cơng thức tính có dạng [2, 3, 6, 7]:   vx S  x       v S    vy S y   Mơ hình xử lý ảnh camera [1]: Cơng nghệ giải đốn đường bờ đặc trưng sóng, dịng chảy ven bờ Video-Camera phương pháp quan trắc mới, đại áp dụng để phân tích diễn biến bãi biển, đặc trưng thủy động lực ven bờ Phương pháp cho phép quan trắc liên tục với khoảng thời gian đo đạc rộng, từ vài giây đến nhiều năm quy mô không gian từ mét đến km Việc phân tích, xử lý số liệu từ camera dựa vào nguyên tắc  y n  A  yeqn 1  yeqn   y n 1 A  , A  k t (2) Với n số theo thời gian Hệ số biến đổi đường bờ tham số hóa theo đặc trưng mực nước, sóng, trầm tích Trong nghiên cứu này, xây dựng chương trình tính tốn ngơn lập trình Fortran sử dụng thuật tốn tự hiệu chỉnh mơ hình nhóm tác giả Jon K Miller Robert G Dean (2004) nhằm lựa chọn tham số thích hợp địa điểm trong thực tế [6, 7] Dữ liệu đầu vào cho mơ hình gồm có kích thước hạt trầm tích, tham số hình dạng mặt cắt ngang cân bằng, tham số sóng vỡ mực nước Thuật tốn tự hiệu chỉnh mơ hình cho phép thu hệ số tóc độ biến đổi k, độ lệch vị trí đường bờ so với trạng thái cân biến đổi chúng Mơ hình kết hợp với mơ hình tính tốn sóng (tham số đầu vào) xử lý ảnh camera (dữ liệu phục vụ hiệu chỉnh) Mơ hình phổ sóng dừng EBED [5, 8, 9] Phương trình phổ sóng dừng có dạng:    CCg cos  S y 2  Trong đó: S mật độ phổ tần số góc hướng; (x, y) toạ độ phương ngang; θ hướng sóng tới tính từ trục x ngược chiều kim đồng hồ; ω tần số; C tốc độ sóng; Cg- tốc độ nhóm; h mực nước tĩnh; Ƙ, K tham số tự để tối ưu hóa tác động trình nhiễu xạ, khúc xạ tiêu tán lượng; Sstab mật độ phổ sóng ổn định Với vx, vy, v vận tốc lan truyền theo hướng tọa độ tương ứng 28 y n 1   y  K  CCg cos2  S yy   Cg  S  S stab  (3)  h hình học ảnh Xét hệ tọa độ xyz, với trục x vng góc bờ biển dương hướng xa bờ, trục y vng góc với trục x, trục z hướng thẳng đứng lên phía với mực chuẩn tham chiếu (z = 0), thường đặt trùng với mực nước triều trung bình mực chuẩn quốc gia [1, 4] Hình Sơ đồ quan hệ hình học tâm camera (X0, Y0, Z0), tọa độ ảnh (u, v) tọa độ thực (X, Y, Z) The scale of the central beach change processes Hình Sơ đồ hệ thống mơ hình cơng cụ Hệ phương trình quang trắc gồm phương trình liên kết tọa độ mặt phẳng cảm biến (2D) hay mặt phẳng ảnh (u, v) với tọa độ vật thể (3D) hay tọa độ thực (X, Y, Z) Các phương trình thiết lập từ phép chiếu tâm từ điểm vật thể qua tâm quang camera đến ảnh mặt phẳng ảnh Hệ phương trình sau:  m ( x  xc )  m12 ( y  yc )  m13 ( z  zc )  u  u0   f / u  11   m31 ( x  xc )  m32 ( y  yc )  m33 ( z  zc )   m ( x  xc )  m22 ( y  yc )  m23 ( z  zc )  v  v0   f / v  21   m31 ( x  xc )  m32 ( y  yc )  m33 ( z  zc )  Với mij ma trận 3×3 góc nghiêng (τ),  cos( ) sin(  )    M   sin(  ) cos( )  cos( )  0  sin(  )  Hệ phương trình (4) bao gồm 11 số chưa biết: Góc nghiêng (τ), góc phương vị (φ), góc quay (σ), tọa độ thực tâm camera (xc, yc, zc), tọa độ tâm ảnh (u0; v0), chiều dài tiêu cự f, hệ số tỷ lệ λu, λv Kết nối mơ hình: Mơ hình biến đổi vị trí đường bờ kết hợp với mơ hình tính sóng mơ hình xử lý ảnh camera tạo thành hệ thống mơ trình bày hình Cơng cụ phân tích hàm điều hịa: Trên giới áp dụng phương pháp phân tích hàm điều hịa để phân tích diễn biến đường bờ (4) phương vị (φ), góc quay (σ):   cos( )   sin(  )   sin(  ) cos( )   sin(  )   cos( )   Từ đó, xác định đươc quy mơ thời gian độ lớn thành phần điều hòa khác đóng góp vào diễn biến vị trí đường bờ Có thể kể đến tác Fi-John Chang Horng-Cherng Lai (2014) áp dụng phân tích diễn biến đường bờ Đài Loan [10] Bản chất phương pháp triển khai vị trí đường bờ theo chuỗi Fourier hay hàm tuần hoàn tách thành tổng hàm dao động đơn lẻ sin cosin: Giả thiết Yˆt giá trị nhận từ phân tích điều hịa Yt: N 2 2   Yˆt  Yt   t      Ai sin it  Bi cos it , t  1, 2, N P P  i 1  (5) 29 Vu Cong Huu, Dinh Van Uu Với: µ giá trị trung bình, Ai Bi hệ số Fourier cho ngày năm, P chu kỳ đầy đủ hay chu kỳ hàm tuần hoàn, P có giá trị N khoảng gián đoạn thời gian 1, số i số hiệu hài điều hịa thứ i Dựa phương pháp bình phương tối thiểu biến µ, Ai Bi tính theo cơng thức: P 2 2   Yˆt  Yt   t      Ai sin it  Bi cos it , t  1, 2, P P P  i 1   P t 1 P P P  2 2  Yˆt  Yt    Yt     Ai sin it   Bi cos it  P P  t 1  i 1 i 1   P   Y cos  t  Nt  i 1 P P  Trong công thức trên, số hạng chứa Ai sin i  (7) (9) sin cos viết gộp lại: 2 2 2   it  Bi cos it  Ci cos i  t   i  P P P   Ai2  Bi2 gọi biên độ hài thứ i, A P arc  i  pha hài điều hòa thứ i 2 i  Ai  Với P chu kỳ bản, hay chu kỳ đầy đủ hàm tuần hoàn, đơn vị đo đơn vị thời gian P không luôn N N trị số khoảng gián đoạn quan trắc Δt = Đại lượng i gọi số hiệu hài điều hòa số nguyên 2/N Đơn vị đo t P phải Hai số hạng dấu tổng biến thiên chu trình đầy đủ chu kỳ Các số hạng thứ biến thiên nhanh gấp đơi, hồn thành chu trình đầy đủ thời gian nửa chu kỳ Số hạng cuối có chu kỳ biến thiên N/2P Dao động với số hiệu i đóng góp vào phương sai (Dy) chung hàm Y lượng Ci2 , ngoại trừ hài cuối cùng, Ci2 Tỷ lệ % đóng góp hài thứ i Ci2 100% Dx 30 P P P P  2   2  Y , A  Y sin it , B  Y cos it , i  1, 2,  (8)    t i t i t     P i 1 P i 1 P i 1  P   P  AP  0, BP  Với Ci  (6) (10) ÁP DỤNG MÔ PHỎNG BIẾN ĐỘNG BÃI BIỂN TRUNG TÂM THUỘC BỜ TÂY VỊNH NHA TRANG Hiệu chỉnh kiểm định mơ hình tính sóng: Mơ hình sóng hiệu chỉnh kiểm định dựa số liệu sóng thực đo thời đoạn mùa gió Đơng Bắc gió mùa Tây Nam Sau hiệu chỉnh kiểm định, mơ hình sử dụng để khơi phục trường sóng phục vụ làm đầu vào cho mơ hình biến đổi đường bờ [5] Mực nước phục vụ tính toán số liệu thực đo trạm hải văn Cầu Đá, Nha Trang Hiệu chỉnh mơ hình xử lý ảnh camera: Mơ hình xử lý ảnh camera cho kết vị trí đường bờ (đường mép nước), đặc trưng hình thái bãi đặc trưng sóng ven bờ Trong giới hạn nghiên cứu này, xét đến vị trí đường bờ Số liệu phục vụ cho việc hiệu chỉnh mơ hình xử lý ảnh camera gồm tọa độ điểm khống chế (GCP-ground control points) tọa độ vị trí đường mép nước (đường bờ) 20 điểm đo máy toàn đạc [1] The scale of the central beach change processes Hình So sánh tọa độ thực đo tính tốn Sai số đánh giá thông qua công thức sai số bình phương quân phương (NMSE) Kết thu RMSE = 3,83 Theo tác giả Dean, khoảng sai số đánh giá dự báo sau: < NMSE < 0,3 (rất tốt), 0,3 < NMSE < 0,6 (tốt), 0,6 < NMSE < 0,8 (hợp lý) 0,8 < NMSE 0,55) hệ số tương quan mực nước so với vị trí đường bờ giảm xuống (0,03 < 0,35) Như vậy, vai trị sóng tháng mùa gió Đông Bắc chiếm ưu nhiều so với vai trị mực nước đến q trình diễn biến đường bờ Trong mùa gió Tây Nam, độ cao sóng nhỏ (đa phần nhỏ 0,5 m) vai trị sóng lại nhỏ so với vai trị mực nước đến q trình diễn biến đường bờ QUY MƠ THỜI GIAN CỦA CÁC Q TRÌNH Kết phân tích điều hịa cho thấy, biến đổi vị trí đường bờ chủ yếu thành phần định: Thành phần thứ The scale of the central beach change processes ứng với số hiệu i = 1, có chu kỳ 1.097 ngày (3 năm), biên độ 0,81 m, pha ban đầu 621 ngày, dao động gây nên 1,5% phương sai chung dao động năm Thành phần thứ ứng với số hiệu i = 3, có chu kỳ 365 ngày, biên độ m, pha ban đầu 86 ngày, góp 80,3% vào phương sai chung dao động Thành phần thứ ứng với số hiệu i = 6, có chu kỳ 182 ngày, biên độ 2,3 m, pha ban đầu 101 ngày, góp 12% vào phương sai chung dao động Còn lại thành phần có chu kỳ dao động biến thiên từ quy mơ kiện đến quy mô tháng chiếm tỷ trọng 7,7% lại gây biến động vị trí đường bờ cục làm biến động bãi biển (hình 15) Hình 15 So sánh diến biến đường bờ sau lọc chu kỳ nửa năm, năm năm Hình 16 Diễn biến đường bờ theo phương pháp vị trí 200 m Tổng hợp hai thành phần thứ thứ đóng góp 92,3% phương sai vào dao động năm Kết hợp hai dao động tạo nên biến trình điển hình vị trí đường bờ gồm cực đại vào thời điểm tháng (trong mùa gió Tây Nam) cực tiểu vào tháng (trong mùa gió Đơng Bắc) đặc trưng cho bãi biển Nha Trang Kết hợp hai dao động tạo nên biến trình điển hình gồm cực đại vị trí đường bờ vào thời điểm tháng (trong mùa gió Tây Nam) cực tiểu vào tháng (trong mùa gió Đơng Bắc) đặc trưng cho bãi biển Nha Trang Như vậy, qua việc phân tích xác định thành phần gây biến động vị trí đường bờ với chu kỳ, biên độ pha khác Từ đó, việc dự báo diễn biến đường bờ theo công thức (6) Số liệu thực đo lưu lượng trạm thủy văn Đồng Trăng sông Cái cho thấy suất lũ sớm vào khoảng tháng đến lũ vụ vào khoảng tháng đến tháng 12 So sánh quy mô thời gian với kết mơ hình, phân tích ảnh camera phân tích điều hịa cho thấy có khớp pha hay trình diễn biến bãi 35 Vu Cong Huu, Dinh Van Uu biển quan tâm diễn biến theo chu kỳ lũ sông Cái Như vậy, nguồn trầm tích từ sơng yếu tố có ảnh hưởng đến qua trình bồi, xói bãi biển quan tâm DIỄN BIẾN ĐƯỜNG BỜ BÃI BIỂN THƠNG QUA PHÂN TÍCH ẢNH VỆ TINH LANDSAT Cửa sơng nguyên nhân gây trình bồi xói bãi biển trung tân thuộc bờ tây vịnh Nha Trang Chính vậy, việc phân tích diễn biến cửa sơng góp phần bổ sung thêm thơng tin nguyên nhân gây bồi xói khu vực bãi biển Cửa sông Cái bãi biển quan tâm phân tích qua ảnh vệ tinh nhằm xác định xu diễn biến qua thời kỳ đường bờ Kết phân tích cho thấy, thời gian từ năm 1988 đến năm 2014 khu vực cửa sơng Cái có biến động mạnh nhất, đặc biệt mũi cát (spit) phía nam cửa sơng Diễn biến cửa sơng Cái thời kì 1988–1995 Diễn biến cửa sơng Cái thời kì 1988–1999–2014 Diễn biến cửa sơng Cái thời kì tháng 6/1988 tháng 6/1996 Diễn biến cửa sông Cái thời kì tháng 2, 6, năm 1996 tháng 3/1997 Hình 17 Diễn biến cửa sơng bãi biển qua thời kỳ Từ năm 1988 đến năm 1995 khơng có thay đổi, độ rộng cửa sông dao động khoảng từ 80–120 m Đến năm 1997, có thay đổi hình dạng mũi cát phía nam Đến năm 2014 khu vực bãi giáp cửa có biến động mạnh bao gồm tượng bãi cát xu biển tiến vào thời kì 2013, 2014, khoảng 36 cách đo khoảng 10 m Vào tháng năm 1996 khu vực bãi biển phía trước UBND tỉnh bị xói sâu vào (khoảng 14 m) khu vực lại tái tạo vào tháng năm sau (1996), khả tái tạo bãi gần hồn tồn so với trước bị xói The scale of the central beach change processes Hình 18 Diễn biến vị trí đường bờ qua năm CÁC NGUYÊN NHÂN GÂY BIẾN ĐỘNG BÃI BIỂN Kết tính tốn mơ hình mơ biến đổi vị trí đường bờ phân tích thành thành phần với biên độ chu kỳ khác cho thấy: Các thành phần quy mô ngắn hạn, kiện làm biến đổi bãi biển cục bãi biển khơi phục lại sau kiện Các thành phần có chu kỳ biến động theo theo mùa, gắn với quy mơ gió mùa Đơng Bắc Tây Nam nguyên nhân tác động theo mùa Các hệ số tương quan cho thấy vai trị sóng chủ đạo Các thành phần chu kỳ dao động nửa năm năm trùng với chu kỳ lũ lũ sớm hay chu kỳ cung cấp nguồn trầm tích sơng Cái Các kết phân tích ảnh vệ tinh cho thấy, q trình xói bãi biển xuất đồng thời với biến roi cát cửa sông Như vậy, nguyên nhân thiếu hụt nguồn trầm tích Kết tính tốn dịng vận chuyển dọc bờ cho thấy bãi biển trung tâm thuộc bờ tây vịnh có xu bồi thực tế bãi biển xói Do đó, tồn q trình vận chuyển trầm tích ngag bờ làm xói bãi biển thiếu hụt nguồn trầm tích từ cửa sơng Cái Tổng hợp nguyên nhân cho thấy, thiếu hụt nguồn trầm tích từ sơng kết hợp với điều kiện sóng gió mùa Đơng Bắc ngun nhân gây biến đổi bãi biển trung tâm thuộc bờ tây vịnh Nha Trang Hình 19 Sơ đồ nguyên lý biến đổi bãi biển KẾT LUẬN Nghiên cứu lựa chọn ứng dụng thành công mơ hình lai ghép thống kê số trị mô quy luật biến động nhân tố thủy động lực hình thái bãi biển trung tâm thuộc bờ tây vịnh Nha Trang, đáp ứng nhu cầu mô quy mô từ ngấn hạn đến dài hạn (Quy mơ cơng trình), tích hợp với mơ hình 37 Vu Cong Huu, Dinh Van Uu sóng, mơ hình xử lý ảnh camera cơng cụ xử lý, phân tích kết Những kết thu góp phần bổ sung tạo tiền đề cho nghiên cứu mô biến đổi đường bờ bãi biển nước ta để nghiên cứu ngày phát triển hoàn thiện tương lai Các kết từ mơ hình phân tích xử lý theo phương pháp phân tích điều hịa phân tích tương quan xác định quy mô thời gian trình chưa làm rõ nghiên cứu trước Từ đó, xác định rõ quy mơ nguyên nhân trình làm biến động bãi biển Lời cảm ơn: Nhóm tác giả xin cảm ơn ban chủ nhiệm đề tài cấp nhà nước KC.09.14/16–20 GS.TS Đinh Văn Ưu làm chủ nhiệm tạo điều kiện kinh phí sử dụng số liệu, tài liệu phần mềm liên quan đến thủy thạch động lực vùng bờ TÀI LIỆU THAM KHẢO [1] Nguyen Trung Viet, 2014 Project of protocol “Study on hydrodynamic regime and sediment transport in estuarine and coastal zones of Nha Trang bay, Khanh Hoa province, 2013– 2014” (in Vietnamese) [2] Vu Cong Huu, Dinh Van Uu, Nguyen Kim Cuong, Le Xuan Hoan, Duong Cong Dien, Duong Hai Thuan, 2014 Toward a prediction and warning system of shoreline change due to water level and wave conditions Vietnam association for fluid mechanics ISBN: 1859-4182, pp 287–295 (in Vietnamese) [3] Vu Cong Huu, Nguyen Kim Cuong, Dinh Van Uu, Nguyen Minh Huan, Nguyen 38 Trung Viet, 2015 Characteristics of wave field and shoreline change in the Nha Trang bay VNU Journal of Science: Natural Sciences and Technology, 31(3S), 179–185 (in Vietnamese) [4] Le Thanh Binh, 2017 Study on shoreline evolution and structural measures for beach protection in Nha Trang city Doctoral dissertation, Code No: 62-58-02-02, Thuyloi University, Hanoi (in Vietnamese) [5] Vu Cong Huu, Dinh Van Uu, 2016 Calculation of wave regime and longshore sediment transport in Nha Trang bay, Khanh Hoa province VNU Journal of Science: Earth and Environmental Sciences, 32(3S), 122– 129 (in Vietnamese) [6] Miller, J K., and Dean, R G., 2004 A simple new shoreline change model Coastal Engineering, 51(7), 531–556 [7] Miller, J K., and Dean, R G., 2005 A simple new shoreline evolution model In Coastal Engineering 2004: (In Volumes) (pp 2009–2021) [8] Mase, H., 2001 Multi-directional random wave transformation model based on energy balance equation Coastal Engineering Journal, 43(04), 317–337 [9] Mase, H., Oki, K., Hedges, T S., and Li, H J., 2005 Extended energy-balanceequation wave model for multidirectional random wave transformation Ocean Engineering, 32(8–9), 961–985 [10] Chang, F J., and Lai, H C., 2014 Adaptive neuro-fuzzy inference system for the prediction of monthly shoreline changes in northeastern Taiwan.Ocean engineering, 84, 145–156  ... bờ biển Hội An–Cửa Đại, bãi biển Nha Trang, … Hình Vịnh Nha Trang (phải) bãi biển trung tâm (trái) 26 The scale of the central beach change processes Vịnh Nha Trang 29 vịnh đẹp giới, có chiều... Project of protocol “Study on hydrodynamic regime and sediment transport in estuarine and coastal zones of Nha Trang bay, Khanh Hoa province, 2013– 2014” (in Vietnamese) [2] Vu Cong Huu, Dinh Van... Nguyen Minh Huan, Nguyen 38 Trung Viet, 2015 Characteristics of wave field and shoreline change in the Nha Trang bay VNU Journal of Science: Natural Sciences and Technology, 31(3S), 179–185 (in Vietnamese)