Tính toán quá trình lan truyền sóng trong rừng ngập mặn Nguyê ̃ n Viết Quỳnh Trường Đại học Khoa học Tự nhiên Luận văn ThS ngành: Hải dương học; Mã số: 60 44 97 Người hướng dẫn: PGS. TS. Nguyê ̃ n Minh Huấn Năm bảo vệ: 2012 Abstract: Phân tích quá trình tiêu tán năng lượng sóng trong rừng ngập mặn. Đưa ra Mô hình tính toán và dự báo sóng SWAN. Áp dụng mô hình SWAN tính toán suy giảm năng lượng sóng do rừng ngập mặn tại khu vực cửa Trà Lý – Thái Bình. Keywords: Hải dương học; Quá trình lan truyền sóng; Rừng ngập mặn; Tính toán Content Mở đầu Rừng ngập mặn (RNM) phát triển chủ yếu tồn tại và sinh trưởng tại khu vực bờ biển ngập nước. Rừng ngập mặn trên thế giới rải rác chủ yếu ở những khu vực nhiệt đới và cận nhiệt đới vì chúng không chịu được sự lạnh giá ( Taal, 1994). Trong quá khứ rừng ngập mặn được xem như là vô dụng và như một kết quả chúng bị biến mất nhanh chóng [3]. Chỉ gần đây rừng ngập mặn được nhận là một hệ sinh thái quan trọng, không chỉ vẻ đẹp mà còn vì tầm quan trọng đối với sự ổn định đường bờ biển và môi trường nuôi dưỡng nhiều sinh vật biển. Đặc biệt những khu rừng ngập mặn gần bờ có thể làm giảm độ cao sóng và thậm chí sóng thần. Tháng 7 năm 1996, cơn bão số 2 (Frankie) với sức gió 103 ÷ 117 km/giờ đổ bộ vào huyện Thái Thụy (Thái Bình) nhờ có dải RNM bảo vệ nên đê biển và nhiều bờ đầm không bị hư hỏng, trong lúc đó huyện Tiền Hải do phá phần lớn RNM nên các bờ đầm đều bị xói lở hoặc bị phá vỡ. Năm 2005, vùng ven biển huyện Thái Thụy tuy không nằm trong tâm bão số 7 (Damrey) nhưng sóng cao ở sông Trà Lý đã làm sạt lở hơn 650m đê nơi không có RNM ở thôn Tân Bồi, xã Thái Đô trong lúc phần lớn tuyến đê có RNM ở xã này không bị sạt lở vì thảm cây dày đặc đã làm giảm đáng kể cường độ sóng. Trận sóng thần khủng khiếp tại Ấn Độ Dương tháng 12 năm 2004, Kathiresan and Rajendran (2005) đã cho thấy tầm quan trọng của rừng ngập mặn trong việc làm suy giảm ảnh hưởng của sóng thần.Ví dụ, tại Indonesia tâm sóng thần rất gần với đảo Simeuleu, tuy nhiên số lượng người chết đặc biệt thấp bởi vì sự hiện diện của những khu rừng ngập mặn với mật độ dày đặc, phía đông nam của Ấn Độ, thiệt hại về kinh tế và con người ít tại những vùng có rừng ngập mặn rậm rạp. Do tầm quan trọng to lớn của chúng, rừng ngập mặn và hệ sinh thái rừng ngập mặn đã được chú trọng nghiên cứu trong thời gian gần đây nhưng sự hiểu biết về chúng còn rất hạn chế. Thiếu trầm trọng những nghiên cứu động lực và nghiên cứu về sự tương tác giữa trầm tích học và thủy động học. Đặc biệt, những nghiên cứu về quá trình lan truyền sóng trong rừng ngập mặn là không nhiều. Trong bối cảnh như vậy, đã lựa chọn “Tính toán quá trình lan truyền sóng trong rừng ngập mặn” là tên của đề tài nghiên cứu. Bố cục luận văn gồm 3 chương và phần kết luận: Chương 1: Quá trình tiêu tán năng lượng sóng trong rừng ngập mặn Chương 2: Mô hình tính toán và dự báo sóng SWAN Chương 3: Áp dụng mô hình SWAN tính toán suy giảm năng lượng sóng do rừng ngập mặn tại khu vực cửa Trà Lý – Thái Bình Kết luận và kiến nghị CHƢƠNG 1 – QUÁ TRÌNH TIÊU TÁN NĂNG LƢỢNG SÓNG TRONG RỪNG NGẬP MẶN Trong chƣơng này, tổng quan về rừng ngập mặn nói chung, phân bố rừng ngập mặn trên thế giới và ở Việt nam, các đặc điểm đƣợc giới thiệu trong phần đầu. Trong phần tiếp theo trình bày cơ sở khoa học của lý thuyết sóng tuyến tính, các phƣơng pháp xác định tiêu tán năng lƣợng sóng. 1.1 Rừng ngập mặn Rừng ngập mặn (RNM) là rừng có những loại cây đặc biệt, thường mọc ở ranh giới giữa những phần đất tiếp giáp bờ biển và biển, ở vùng nhiệt đới và cận nhiệt đới. RNM thường phát triển trong những vùng triều, giữa mực nước biển trung bình và mực nước triều lớn. Phía ngoài biển của rừng ngập mặn thường là vùng đất bùn, với độ dốc khoảng 1:1000. Phía sau RNM có thể là những đê biển, đầm muối và khu vực dân cư sinh sống. Phía trong RNM cây sẽ mọc song song với đường bờ biển với nhiều loài khác nhau. 1.2 Lý thuyết sóng tuyến tính và phổ sóng 1.2.1 Lý thuyết sóng tuyến tính Sóng được miêu tả một cách tổng quát bởi lý thuyết sóng tuyến tính hay lý thuyết sóng Airy (Airy, 1845). Trong đó: - Độ sâu nước là hằng số - Chuyển động của sóng là hai chiều - Sóng chuyển động trong dạng nhất định bất biến theo thời gian. - Chất lỏng không nén được - Nhớt, rối và sức căng bề mặt được bỏ qua - Chiều cao sóng (H) nhỏ so với chiều dài sóng (L) và độ sâu (d) Phương trình điều chỉnh vận tốc thế vị  (Young, (1999)),              zx wu  ,, theo hai hướng (x,z) 0 2 2 2 2       zx  (1.1) Điều kiện biên động lực trên bề mặt: 0           g t tz 0z Điều kiện biên đáy: 0   z  hz  Ở đây ɳ (m) là dao động mực nước, g (m/s 2 ) là gia tốc trọng trường,  thế vận tốc, u (m/s) là vận tốc theo hướng x, w (m/s) là vận tốc theo hướng z. Phân loại sóng theo độ sâu Sóng có thể được phân loại dựa theo mối quan hệ giữa chiều dài sóng với độ sâu hay tích số của số sóng k và độ sâu nước d (kd). Hầu hết tham số sóng sẽ thay đổi khi sóng lan truyền từ vùng nước sâu vào vùng nước nông. Chiều dài sóng, vận tốc nhóm sóng và chuyển động quỹ đạo của hạt nước sẽ thay đổi. 1.2.2 Năng lượng sóng và phổ năng lượng sóng Năng lượng sóng bao gồm động năng và thế năng (Battjes, (1998)). Năng lượng sóng tổng cộng bằng tổng động năng và thế năng (mật độ năng lượng) trên một đơn vị chiều rộng được tính bằng: 2 8 1 gHE   (1.2) Trong đó, là mật độ nước (kg/m), g là gia tốc trọng trường (m/s 2 ) và H là độ cao sóng (m). Phổ sóng Một phương pháp để nghiên cứu các quá trình dao động là xem xét chúng là một tổ hợp các sóng hình sin với độ lớn, tần số và pha khác nhau. Các hàm số thể hiện sự phân bố của biên độ và pha theo tần số được gọi là phổ tần số. Phổ có thể là liên tục hay rời rạc, nhưng phổ sóng là liên tục. Phổ mật độ năng lượng sóng đặc trưng được thể hiện trong hình 1.5 dưới đây. Trục ngang thể hiện tần số, trục thẳng đứng thể hiện mật độ năng lượng. 1.3 Các phƣơng pháp xác định tiêu tán năng lƣợng sóng 1.3.1 Phương pháp sử dụng hệ số ma sát đáy Sóng lan truyền ở vùng nước nông gần bờ chịu tác động ma sát do đáy sinh ra. Độ lớn của lực ma sát này phụ thuộc vào một loạt các tham số như vận tốc quỹ đạo sóng theo phương ngang, độ nhám của đáy và độ sâu nước. Trong nhiều tài liệu về hiện tượng suy giảm sóng trong rừng cây ngập mặn hiện tượng ma sát đáy rất ít được chú ý mặc dù các biểu thức của lực ma sát thường được sử dụng để mô tả các thành phần tiêu tán do rừng cây ngập mặn. Công thức của Van Rijn Theo Van Rijn (1989), tiêu tán năng lượng do ma sát đáy được thể hiện bằng đại lượng D f là công trung bình theo thời gian của lực ma sát   T bf dtU T D 0 1   (1.11) trong đó T là thời gian, τ b là ứng suất tiếp đáy, U δ vận tốc quỹ đạo sóng theo phương ngang. Công thức của Collins Ma sát đáy theo Collins (1972) dựa trên công thức thông thường đối với sóng có chu kỳ với các tham số phụ để phù hợp với trường sóng ngẫu nhiên. Cường độ tiêu tán năng lượng được xác định bằng công thức ma sát đáy thông thường: C bot = c f . g. U orb (1.15) trong đó c f là hệ số ma sát Collins c f = 0.015. 1.3.2 Phương pháp hình khối trụ Một phương pháp khác dựa trên công thức lý thuyết đối với tiêu tán năng lượng sóng do thực vật với các biểu thức riêng biệt trong đó xác định phần năng lượng bị tiêu tán do các quá trình quan trọng nhất Công thức của Morrison Morrison và nnk (1950) đã xác định được lực tác động lên một cọc mảnh đơn lẻ trong sóng bề mặt:   UUDC dt dU DCtf dm  2 1 4 1 2  (1.18) trong đó C m là hệ số quán tính, D là đường kính của cọc, C d hệ số tiêu tán và U là giá trị cực đại của vận tốc quỹ đạo theo phương ngang. Công thức của Dalrymple Dalrymple và nnk (1984) xác định một hệ số tiêu tán năng lượng dựa trên biểu thức của Morrison đối với các thực vật chìm trong nước ở đầm lầy. 3 3 3 3 cosh3 sinh3sinh 23 2 H khk hkhkkg NbC vvdv              (1.22) trong đó k là số sóng, σ là tần số góc, αh là độ cao trung bình của thực vật. Hiện tượng phản xạ sóng từ thực vật không được tính đến trong công thức này. Hai cơ chế cơ bản để mô tả hiện tượng tiêu tán năng lượng sóng do thực vật được mô tả ở các phần phía trên là phương pháp khối trụ và tăng cường tham số nhám đáy. Phương pháp tăng cường tham số nhám đáy dựa trên các kết quả của Collins (1968) được Vries và Roelvink (WL Delft Hydraulíc, 2004) nghiên cứu và áp dụng bổ sung vào mô hình SWAN, các kết quả này được hiệu chỉnh cho cỏ biển - Spartina Anglica. Phương pháp này có vẻ phù hợp khi hiệu chỉnh nhưng không tồn tại nguyên tắc vật lý để tham số hóa các dạng thực vật, các dạng thực vật phức tạp không dễ dàng đưa vào mô hình do thực vật chỉ được tham số hóa bằng một tham số. Mendez và Losada (2004) đã mở rộng công thức của Dalrymple (1984), kết quả này có thể áp dụng cho vùng đáy nghiêng và sóng vỡ. Các thí nghiệm vật lý trong máng chứng minh các kết quả rất triển vọng của phương pháp này, Mendez và Losada cho thấy rằng dạng mô hình này có thể dễ dàng áp dụng trong các mô hình tính toán lan truyền sóng chuẩn để tích hợp quá trình tiêu tán năng lượng sóng do thực vật và còn có thể dự tính dòng chảy sóng cũng như vận chuyển trầm tích trong vùng có thực vật. Kết luận là công thức của Dalrymple (1984) là phương pháp xấp xỉ quá trình tiêu tán năng lượng sóng do thực vật tốt nhất và thích hợp nhất để tích hợp vào mô hình SWAN. CHƢƠNG 2 - MÔ HÌNH TÍNH TOÁN VÀ DỰ BÁO SÓNG SWAN Giới thiệu mô hình tính toán sóng SWAN phiên bản mới nhất, đã trình bày tóm tắt cơ sở khoa học và thuật toán của mô hình SWAN. Phần phân tích thành phần tiêu tán đã tập trung vào thành phần tiêu tán năng lƣợng do cây để phục vụ thực hiện các nội dung của đề tài. 2.1. Giới thiệu mô hình SWAN SWAN là mô hình tính toán sóng thế hệ ba, tính toán phổ sóng hai chiều bằng cách giải phương trình cân bằng tác động sóng (trong trường hợp không có dòng chảy có thể dùng phương trình cân bằng năng lượng sóng) có tính tới sự lan truyền sóng từ vùng nước sâu vào vùng nước nông ven bờ, đồng thời trao đổi năng lượng với gió thông qua hàm nguồn cùng với sự tiêu tán năng lượng sóng. 2.1.1. Cơ sở lý thuyết của mô hình SWAN - Phương trình cân bằng tác động phổ: Trong SWAN sự tiến triển của phổ sóng được mô tả bằng phương trình cân bằng tác động phổ được viết trong hệ toạ độ Đề Các :   S NCNCNC y NC x N t yx                (2.1) Trong hệ toạ độ địa lý phương trình có thể viết dưới dạng:         S NCNCNCNCN t                 coscos 1 (2.2) Với  là kinh độ,  là vĩ độ. * Năng lượng gió truyền cho sóng: ),(),(  BEAS in  (2.3) Với A là hệ số tăng tuyến tính, B là hệ số tăng theo hàm mũ. Hệ số A và B phụ thuộc vào tần số và bước sóng đồng thời phụ thuộc vào vận tốc và hướng của gió. + Hệ số tăng tuyến tính A: Biểu thức của Cavaleri và Malanotile-Rizzli (1981) được sử dụng để loại bỏ sự tăng trưởng của sóng tại các tần số thấp hơn tần số Pierson-Moskovitz.    HU g A w 4 * 2 3 )cos(,0max 2 10.5,1     (2.4) Trong đó   là hướng gió, H là giá trị để lọc + Hệ số tăng theo hàm mũ B: Trong mô hình SWAN sử dụng biểu thức tính hệ số tăng theo hàm mũ B: Công thức của Kome và nnk (1984):                     1)cos(2825.0,0max * W phw a C U B (2.5) Với C ph là vận tốc pha,  a và  w là mật độ không khí và nước. - Mất mát năng lượng trong khi truyền sóng (S ds ): Mất mát năng lượng khi truyền sóng gây ra do ba quá trình: sự bạc đầu của sóng, ma sát đáy và sóng vỡ do ảnh hưởng của độ sâu địa hình. - Ma sát đáy. Mô hình ma sát đáy sử dụng trong SWAN bằng mô hình thực nghiệm của Jonswap, mô hình sức cản của Collins (1972) và mô hình nhớt rối của Madsen (1980). Công thức sử dụng trong các mô hình này là. ),( )(sinh ),( 22 2 ,    E kdg CS bottombds  (2.6) Ở đây C bottom là hệ số ma sát đáy, hệ số này phụ thuộc vào vận tốc quỹ đạo chuyển động sóng. - Sóng vỡ do ảnh hưởng của độ sâu địa hình. 2 _ 24 1 mbbjtot HQD               (2.7) Với  bj =1, Q b là phân số sóng vỡ, H m là chiều cao sóng cực đại H m = .d ở đây  là hệ số sóng vỡ, d là độ sâu nước tổng cộng, D tot là tỉ lệ tiêu tán năng lượng trung bình trên một đỉnh sóng. Công thức hàm nguồn năng lượng sóng vỡ có dạng ),(),( ,  E E D S tot tot brds  (2.8) * Tương tác phi tuyến giữa các sóng - Tương tác sóng bậc bốn Tương tác sóng bậc bốn được tính với các xấp xỉ tương tác riêng (DIA) của Hasselmann (1985). Hàm nguồn được cho dưới dạng ),(),(),( ** 14 * 1414  nnn SSS  (2.9) - Tương tác sóng bậc ba Theo xấp xỉ bậc ba Lumped (LTA) của Eldeberky và Battjes. Hàm nguồn có dạng. ),(),(),( 333    nlnlnl SSS (2.10) Với  eb là hệ số - Nước dâng do sóng. Trong mô hình 1D, nước dâng do sóng tính bằng cách tích phân phương trình cân bằng mô men, đó là sự cân bằng giữa tác động sóng (gradient ứng suất phát xạ) và gradient áp suất thuỷ tĩnh. 0    x gdF x  (2.11) Với d là độ sâu nước tổng cộng,  mực nước dâng so với mực nước trung bình. 2.1.2. Mô hình số của SWAN Trong mô hình SWAN, lời giải của phương trình cân bằng tác động sóng được triển khai bằng một số sơ đồ khác nhau trong cả năm chiều (thời gian, không gian địa lý, không gian phổ). - Sơ đồ quá trình truyền sóng Các sơ đồ được sử dụng trong SWAN dựa trên cơ sở có tính năng mạnh, sát thực tế và mang tính kinh tế. Do vậy, sơ đồ được sử dụng là dạng hàm ẩn trong cả không gian địa lý và không gian phổ, thêm vào đó là phép xấp xỉ trung tâm trong không gian phổ. SWAN bao gồm ba sơ đồ : * Sơ đồ BSBT (backward space - backward time): Sơ đồ này tính cho cả trường hợp sóng ổn định và sóng không ổn định. * Sơ đồ S&L Sơ đồ dạng này được sử dụng trong điều kiện sóng ổn định. Thời gian có thể bỏ đi như một biến, các bước lặp vẫn được tiến hành như trước. * Sơ đồ cấp độ gió bậc ba Sơ đồ này được áp dụng cho điều kiện sóng không ổn định. Hai thành phần đạo hàm theo hướng x và y Việc số hoá quá trình khuyếch tán sóng trong sơ đồ S & L được thực hiện được trong trường hợp sóng truyền trên khoảng cách lớn. Hiệu ứng này phụ thuộc vào phổ sóng. Giá trị của hiệu ứng phụ thuộc vào dạng phổ và thời gian sóng truyền. - Hàm nguồn Trong SWAN các hàm nguồn tạo sóng, tương tác sóng và tiêu tán sóng được viết dưới dạng hàm ẩn: )( 1 1 11            nn n nnn EE E S ES (2.12) Với ES   2.1.3. Điều kiện biên và điều kiện ban đầu Trong SWAN, điều kiện biên trong không gian địa lý và không gian phổ là biên hấp thụ đầy đủ với năng lượng sóng đi ra khỏi miền tính hoặc vào bờ. Trên các biên hở, năng lượng của sóng tới được cho bởi người sử dụng. Đối với vùng ven bờ điều kiện biên được cho dọc biên nước sâu và không cho dọc theo biên bên. Điều này có nghĩa là có những lỗi tại biên trong khi truyền vào vùng tính. Vùng ảnh hưởng là vùng tam giác với đỉnh tại góc giữa biên nước sâu và biên bên, phân tán tới bờ trong góc từ 30º tới 45º trên hướng truyền của sóng nước sâu. Để khắc phục vấn đề này thì biên bên nên lấy ở xa vùng cần tính để tránh những lỗi trong vùng tính. 2.2. Bổ sung tiêu tán năng lƣợng sóng do cây trong mô hình SWAN [6] Bổ sung sự tiêu tán năng lượng bởi cây thực vật nghĩa là phải đưa thêm một hạng mới vào phương trình. Với sự bổ sung này mô hình diễn tả quá trình tiêu tán năng lượng [...]... thay đổi mật độ rừng ngập mặn có một số kết luận sau: - Triển khai thành công mô hình SWAN có bổ sung thêm hiệu ứng cây, lưới lồng cho các kịch bản lan truyền sóng qua vùng có và không có sự tồn tại của rừng cây ngập mặn Khi lan truyền vào rừng ngập mặn độ cao sóng giảm và tốc độ suy giảm ở mỗi mực nước là khác nhau Với mực nước nhỏ rễ cây đóng vai trò quan trọng trong quá trình tiêu tán sóng Khi mực... ÁP DỤNG MÔ HÌNH SWAN TÍNH TOÁN SUY GIẢM SÓNG DO RỪNG NGẬP MẶN TẠI KHU VỰC CỬA TRÀ LÝ – THÁI BÌNH Trong chương này, trình bày việc áp dụng mô hình SWAN tính toán suy giảm sóng do rừng ngập mặn tại khu vực Trà Lý tỉnh Thái Bình Trình bày các đặc điểm về sóng, thủy triều, RNM của vùng nghiên cứu Các số liệu đã được xử lý để đưa vào mô hình tính và xây dựng các kịch bản tính toán Mô hình SWAN đã được kiểm... 1: 2000 Độ cao sóng Hs (m) 0.08 – 0.15 – 0.3 – 0.4 Chu kỳ sóng Tp (s) 6.5 3.2.2 Kết quả kiểm nghiệm mô hình Madza đã áp dụng công thức sau để tính suy giảm sóng khi lan truyền vào rừng ngập mặn: R = (Hs – HL)/Hs Hs – độ cao sóng trước RNM, HL – độ cao sóng ở khoảng cách L từ mép ngoài rừng Kết quả tính toán cho vị trí cách bìa rừng 100 m (L = 100 m) được thể hiện qua các kết quả tính toán (xem bảng... dụng, chu kỳ sóng và hướng sóng là sóng hướng Đông được dẫn về biên của lưới tính cho mô hình, một số đặc điểm cơ bản của các tham số sóng tại biên của lưới tính sóng được thể hiện trên bảng 3.11 dưới đây: Bảng 3.11 Đặc trưng đầu vào sóng cho mô hình tính toán Mực nước (m) 0.5- 1.86 – 3.5 Độ cao sóng Hs (m) 1,543 Chu kỳ sóng Tp (s) 5,155 Số liệu cây ngập mặn áp dụng cho vùng tính toán là cây Trang có chiều... quan trọng trong việc bảo vệ bờ biển Luận văn đã tiến hành thiết lập các điều kiện cho mô hình thủy động lực -sóng SWAN 40.85 có bổ sung tiêu tán năng lượng sóng do thực vật để mô phỏng lan truyền sóng trong RNM dựa trên các đặc điểm cây RNM như đường kính thân, chiều cao cây, mật độ cây, từ các kết quả phân tích so sánh trong các trườnng hợp không có rừng ngập mặn và có ảnh hưởng của rừng ngập mặn cũng... 3.17 dưới đây) Tuy nhiên, độ cao sóng giảm mạnh khi vào rừng ngập mặn có mật độ lớn hơn Với mật độ rừng là 5 cây/m2, trong mực nước 3.5 m độ cao sóng giảm từ 0.75 m ở ngoài bìa rừng còn 0.006 m khi vào bờ Khi số cây trong vùng nghiên cứu tăng lên, tức là mật độ tăng hoặc đường kính của thân cây ngập mặn tăng lên thì độ cao sóng giảm càng nhiều hơn Xu hướng biến đổi độ cao sóng khi vào RNM khi thay đổi... ràng, mật độ tăng thì tốc độ suy giảm sóng (r) sẽ tăng hay độ cao sóng sẽ giảm Với mục đích lựa chọn mật độ ổn định cho khu vực nghiên cứu, đảm bảo độ cao sóng là nhỏ nhất phía sau rừng ngập mặn và đảm bảo tính bền vững cho RNM Mật độ rừng nên lựa chọn là 3 cây/m2 và bề rộng rừng là khoảng từ 500 tới 1000 m cho khu vực nghiên cứu Kết luận và đề xuất Các dải rừng ngập mặn ở vùng ven không chỉ có ý nghĩa... nước 1.86 m độ cao sóng sẽ giảm mạnh nhất, khoảng 56%, hơn 18% so với mực nước 0.5 m tại vị trí 500 m cách bìa rừng, xem hình 3.9 Điều này được giải thích rằng sự đóng góp tiêu tán do vòm lá cho tiêu tán năng lượng sóng lớn hơn so với rễ cây Khi sóng lan truyền vào rừng ngập mặn thì độ cao sóng giảm và sự suy giảm sẽ khác nhau ứng với các mực nước khác nhau Độ cao sóng giảm mạnh nhất trong mực nước 1.86... ra sóng to, gió lớn, nước dâng làm phá vỡ đê kè, nhà cửa, biến dạng lòng dẫn và các val, bãi cát ngầm Tại khu vực nghiên cứu và lân cận có mật độ bão đổ bộ khá cao so với các tỉnh ven biển khác 3.1.5 Đặc điểm rừng ngập mặn Tại vùng ven biển Thái Bình hiện nay có 11.750 ha rừng ngập mặn, trong đó rừng nguyên sinh có hơn 500 ha, rừng bần và 3000 ha rừng sú, vẹt nguyên sinh, còn lại hơn 7000 ha là rừng. .. giá ảnh hưởng của RNM trong việc làm giảm độ cao sóng, sử dụng công thức [5]: R = (Hs – HL)/Hs Hs- độ cao sóng trước RNM, HL – độ cao sóng ở khoảng cách L từ mép ngoài rừng - Để đánh giá ảnh hưởng của chỉ riêng yếu tố RNM trong việc làm giảm độ cao sóng, một số tác giả [5] tính toán hệ số suy giảm độ cao sóng do RNM bằng công thức sau: R = (hkhông có RNM – hcó RNM) / hkhông có RNM trong đó: hkhông có . về quá trình lan truyền sóng trong rừng ngập mặn là không nhiều. Trong bối cảnh như vậy, đã lựa chọn Tính toán quá trình lan truyền sóng trong rừng ngập. năng lượng sóng do rừng ngập mặn tại khu vực cửa Trà Lý – Thái Bình. Keywords: Hải dương học; Quá trình lan truyền sóng; Rừng ngập mặn; Tính toán Content