MÔ PHỎNG VẬN CHUYỂN BÙN CÁT VÀ BIẾN ĐỔI ĐỊA HÌNH ĐÁY KHU VỰC CỬA SÔNG THU BỒN

MÔ PHỎNG VẬN CHUYỂN BÙN CÁT VÀ BIẾN ĐỔI ĐỊA HÌNH ĐÁY KHU VỰC CỬA SƠNG THU BỒN Vũ Minh Cát1, Đặng Đình Đoan2, Tóm tắt: Bài báo trình bày kết mơ diễn biến địa hình đáy vùng cửa sông đáy biển trước cửa sông mơ hình MIKE21_FM Couple mơ địa hình đáy theo kịch khác Kết mô cho phép đánh giá định lượng thay đổi địa hình đáy sau trận lũ điển hình, trình xả lũ từ hồ chứa, bùn cát từ sông thay đổi trường mực nước trường sóng phía biển Kết mơ cho thấy phù hợp với số liệu khảo sát qui luật diễn biến cửa sơng Đó tượng xói lòng sau trận lũ bồi cửa nhanh vào mùa nước kiệt kết hợp với gió mùa tây nam Các kết nghiên cứu giúp nhà định lựa chọn giải pháp hợp lý nhằm ổn định cửa sơng phục vụ lũ, giao thơng thủy phòng tránh xói lở, bồi tụ cho khu vực nghiên cứu Từ khóa: Đáy biển, mơ phỏng, sóng nước sâu, vận chuyển bùn cát, MIKE21-FM Couple I MỞ ĐẦU1 Sơng Thu Bồn sơng với diện tích lưu vực khoảng 10,350 km2 thuộc lãnh thổ Quảng Nam thành phố Đà Nẵng Thu Bồn hai nhánh sơng chảy thành phố Hội An phía Nam có tên cửa Đại Với chế độ dòng chảy dồi dào, không năm thể mùa lũ kéo dài tháng với khoảng 80% lượng dòng chảy năm mùa cạn kéo dài tháng mà có 20% lượng nước Vào mùa lũ, lượng phù sa đáng kể mang từ lưu vực xuống vùng đồng hẹp dải ven biển gây tượng biến hình bờ, đáy sơng đường bờ biển [1,2] Do biến đổi dị thường thời tiết gây tương phản ngày khốc liệt mùa mưa mùa khô thêm vào hoạt động kinh tế xã hội thiếu quản lý chặt phá rừng, cày xới bề mặt lưu vực cho mục đích mưu sinh, phát triển hệ thống sở hạ tầng v.v… làm thay đổi chế độ dòng chảy, tăng đột ngột bùn cát làm vấn đề xói lở bồi lấp lòng dẫn, vùng cửa sơng sơng miền Trung nói chung Vu Gia - Thu Bồn nói riêng nghiêm trọng Ở phía biển, với tác động yếu tố thủy động lực sóng, gió, thủy triều, dòng chảy Trường ĐH Thủy lợi Viện Khoa học Thủy lợi Việt Nam 102 tổng hợp mang bùn cát chuyển động dọc bờ với mức độ khác gây xói lở, bồi tụ đường bờ lấp cửa sơng gây khó khăn cho vận tải thủy cửa Đại – cảng cửa sông với số lượng tàu thuyền vào lớn [3] Để xác định nguyên nhân, định lượng q trình vận chuyển trầm tích đánh giá thay đổi địa hình đáy, có nhiều phương pháp sử dụng Đó phương pháp đo đạc định kỳ theo tổ hợp sóng, triều, dòng chảy từ sơng, sau chồng ghép đồ ý đồ chun mơn khác định lượng thay đổi dịa hình đáy theo khơng thời gian Cũng dùng kỹ thuật viễn thám hệ thông tin địa lý để đánh giá biến động địa hình đáy Trong nghiên cứu này, tác giả sử dụng mơ hình số MIKE21 cho nghiên cứu [4] Modul MIKE 21/3 Coupled mô kết hợp chế độ dòng chảy triều sóng sinh kết hợp với modul vận chuyển bùn cát cho ta biết biến động địa hình đáy Mơ hình sử dụng để tính tốn tương tác sóng dòng chảy, kết hợp động lực học mơ đun dòng chảy mơ đun sóng Mơ hình bao gồm tính tốn kết hợp động lực học mô đun vận chuyển bùn, mơ đun dòng chảy sóng Do đó, tác động qua lại đầy đủ thay đổi độ sâu đến tính tốn sóng dòng chảy xem xét Hệ thống phương trình bao gồm: Phương trình liên tục, gọi phương KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ ĐẶC BIỆT (11/2013) trình bảo tồn khối lượng: u v w   =S x y z Phương trình bảo toàn động lượng theo phương ngang (x y), gọi phương trình chuyển động theo phương ngang: u u vu wu    = t x y z n  p g p  ỉ u fv - g -  dz Fn ỗ v ữ u z S x p x p z x z è z ø u u vu wu    = t x y z n - fu - g  - p - g  pdz y  Fv  p y p z y ổ v ỗ v ữ v zS z è z ø đó, t thời gian; x, y z tọa độ Đề các,  dao động mực nước, d độ sâu; h=  d độ sâu tổng cộng, u, v w thành phần vận tốc theo phương x, y z; f=2  sin  tham số Coriolis (  vận tốc góc  vĩ độ địa lý); g gia tốc trọng trường; r mật độ nước, vt nhớt rối thẳng đứng Pa áp suất khí quyển, r mật độ chuẩn; S lưu lượng nguồn nước đổ vào vùng nghiên cứu (us, vs) vận tốc dòng nước chảy từ nguồn vào miền tính Fu, Fv số hạng ứng suất theo phương ngang Phương trình bảo tồn vật chất hay phương trình vận chuyển bùn cát: C uC vC wC    = t x y z = FC  đó: C nồng độ đại lượng vơ hướng; kp tốc độ phân hủy đại lượng đó; Cs nồng độ đại lượng vô hướng điểm nguồn; Dv hệ số khuếch tán thẳng đứng; FC số hạng khuếch tán ngang MIKE 21 áp dụng sơ đồ sai phân hữu hạn phương pháp giải hiệu kỹ thuật ADI (Alternating Direction Imolicit) để giải phương trình bảo tồn khối lượng động lượng miền không gian thời gian II MÔ PHỎNG VẬN CHUYỂN BÙN CÁT VÀ BIẾN ĐỔI ĐỊA HÌNH ĐÁY KHU VỰC CỬA SƠNG THU BỒN Thiết lập miền tính, lưới tính Khu vực nghiên cứu thiết lập thành hai miền lưới tính: miền tính lớn miền tính nhỏ Miền tính lớn thiết lập để hiệu chỉnh kiểm định tạo biên mực nước sóng phục vụ cho tính tốn lưới mịn Miền tính nhỏ để mơ chi tiết biến đổi địa hình đáy khu vực cửa sơng theo kịch khác Miền tính lớn thiết lập từ số liệu đo đạc địa hình phía ngồi khơi lấy từ hải đồ tỷ lệ 1/25000 Hải quân Việt Nam đo đạc xuất năm 2009 Phía cửa sơng sơng lấy từ số liệu thực đo vào tháng 9/2009 Lưới tính bao gồm 9643 phần tử bao trùm toàn bán đảo Sơn Trà, phía Nam kéo qua mũi Tam Hải mở rộng khơi đến độ sâu 100m Miền tính có biên mực nước mực nước phía biển biên lưu lượng Miền tính nhỏ bao gồm cửa sơng cách cửa biển 7km phía thượng lưu, mở rộng phía biển khoảng km sang phía bờ biển khoảng 10 km vi 8093 phn t ổ C ỗ Dv ÷  k p C  CsS z è z ø KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ ĐẶC BIỆT (11/2013) 103 Hình 1: Số hóa địa hình miền tinh lớn HIỆU CHỈNH/KIỂM ĐỊNH MƠ HÌNH Hình 2: Địa hình chi tiết khu vực cửa Đại a Mơ hình triều Chuỗi tài liệu đo đạc liên tục 24/24 16 ngày (từ 15/8/2009 đến 23 ngày 31/8/2009) trạm Sơn Trà cửa Đại sử dụng để hiệu chỉnh mơ hình Kết kiểm định mực nước Sơn Trà cho thấy hệ số tương quan R2 = 0.93 số liệu thực đo mực nước mô trạm cửa Đại, hệ số tương quan R2 = 0.891 Bộ thông số mô triều trạm nói cao hồn tồn cho phép sử dụng mơ kịch Hình 3: Kết hiệu chỉnh mực nước Sơn Trà Hình 4: Kết hiệu chỉnh mực nước cửa Đại b Mơ hình sóng Chuỗi số liệu sóng thực đo 15 ngày từ 07 ngày 06/09/2009 đến 17 ngày 20/09/2009 sử dụng để hiệu chỉnh/kiểm định mơ hình Biên sóng lấy từ số liệu WAVEWATCH III khu vực khơi tọa độ 108°56'29.85"E 16°13'34.52"N Kết kiểm định cho thấy thời gian hai ngày đầu tính tốn kết thực đo mơ hình có sai khác tương đối lớn, ngày có tương đồng chiều cao sóng pha sóng thực đo mơ hình Kết hiệu chỉnh cho kết tương đối tốt, đảm bảo độ xác cho phép sử dụng thông số phục vụ mô kịch c Vận chuyển bùn cát biến đổi địa hình đáy Tại khu vực cửa Đại có số liệu thực đo 9/2009 10/2009 Đây chuỗi số liệu để phục vụ cho hiệu chỉnh để xác định bơ thơng số mơ hình vận chuyển bùn cát sử dụng để mô theo kịch Trên vùng nghiên cứu, chiết xuất mặt cắt để đánh giá kết mô với số liệu thực đo địa hình đáy vùng cửa sơng nghien cứu Các mặt cắt thể hình Kết mơ thể hình vẽ từ đến 10, theo thấy kết mơ cho sai khác khơng lớn dùng thông số để mô theo kịch 104 KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ ĐẶC BIỆT (11/2013) xây dựng phần So sánh địa hình đáy MC1 -1 10 15 20 25 30 Thực đo -2 Tính tốn (m) -3 -4 -5 -6 -7 -8 Hình 7: Địa hình cửa Đại sau mô -2 10 15 20 So sánh địa hình đáy MC3 25 30 Thực đo Tính tốn (m) -3 -4 -5 -6 -7 -8 Hình 9: Địa hình đáy thực đo tính tốn (MC2) KẾT QUẢ MÔ PHỎNG SỰ THAY ĐỔI ĐỊA HÌNH ĐÁY THEO CÁC KỊCH BẢN Sau hiệu chỉnh kiểm định trường mực nước, sóng, dòng chảy vận chuyển bùn cát, thấy thơng số mơ hình chiều sử dụng để mô cho tập kịch đề xuất Mục tiêu ứng dụng mơ hình nói chung để đánh giá diễn biến đường bờ đáy biển, cửa sông khu vực nghiên cứu sau thời khoảng khác mùa, năm nhiều năm Tuy nhiên, hạn chế điều kiện đầu vào thiếu tài liệu mực nước, sóng địa hình đồng khó khăn lớn cơng cụ tính tốn không cho phép mô thời gian dài thời gian chạy máy tính q lâu, gần tương đương với thời gian thực, mơ hình chiều sử dụng để mô biến đổi địa hình đáy khoảng thời gian ngắn, chẳng hạn cho trận lũ cho phép đánh giá biến động địa hình ứng với tổ hợp điều kiện biên xác định bao gồm lưu lượng sơng, mực nước, sóng v.v Trên sở đó, có nhận định diễn biến địa hình đáy Với ý nghĩa đó, tập kịch đề xuất để mô bao gồm trường hợp a Tập kịch mô thay đổi địa hình -1 10 15 20 25 30 Thực đo -2 Tính tốn -3 (m) -1 Hình 8: Địa hình đáy thực đo &tính tốn (MC1) So sánh địa hình đáy MC2 -4 -5 -6 -7 -8 Hình 10: Địa hình đáy thực đo tính tốn (MC3) đáy a.1 Mơ hình lũ 1999 có xét đến ảnh hưởng biến đổi khí hậu -Mô lũ 1999 không xét ảnh hưởng nước biển dâng (A1) -Mô lũ 1999 có xét tới ảnh hưởng nước biển dâng 50cm (A2) -Mô lũ 1999 xét tới ảnh hưởng nước biển dâng 100cm (A3) a.2 Mơ hình lũ 1999 có xét đến điều tiết hồ chứa ảnh hưởng biến đổi khí hậu -Lũ 1999 điều tiết không xét ảnh hưởng tượng nước biển dâng (B1) -Lũ 1999 điều tiết có xét tới ảnh hưởng nước biển dâng 50cm (B2) -Lũ 1999 điều tiết có xét tới ảnh hưởng nước biển dâng 100cm (B3) a.3 Mơ hình dòng chảy mùa cạn có xét đến ảnh hưởng biến đổi khí hậu -Dòng chảy từ đến 18/VII/1999 khơng xét ảnh hưởng nước biển dâng (C1) -Dòng chảy từ đến 18/VII/1999 có xét tới ảnh hưởng nước biển dâng 50cm C2) -Dòng chảy 01 đến 18/VII/1999, có xét tới ảnh hưởng nước biển dâng 100 cm (C3) b Kết mô vận chuyển bùn cát biến đổi địa hình đáy KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ ĐẶC BIỆT (11/2013) 105 Mục tiêu kịch phần đánh giá biến động địa hình đáy sơng địa hình đáy biển khu vực cửa tác động lũ sơng (chọn trận lũ điển hình từ 31/10 đến 11/11/1999), điều kiện bùn cát, mực nước và sóng truyền từ nước sâu vào vùng tính tốn thời gian nói Các liệu đầu vào để tính tốn biến đổi lòng dẫn bao gồm: (i) trình lũ thực đo từ 31/X đến 11/XI/1999 q trình lũ có xét tới q trình điều tiết hồ chứa thượng lưu cầu Câu Lâu tính từ mơ hình MIKE11, (ii) mực nước tính từ mơ hình triều tồn cầu sóng truyền từ khơi vào khu vực nghiên cứu thời gian với trận lũ nêu b.1 Nhóm kịch chọn lũ 1999 có xét đến ảnh hưởng biến đổi khí hậu Các kết trung gian mơ gồm q trình lũ cầu Câu Lâu, trường mực nước biên lỏng phía ngồi biển, trường sóng trường lưu tốc Khi có kết đó, tiến hành mơ biến đổi địa hình đáy cho kịch thuộc nhóm Dưới xin trích kết đánh giá biến đổi địa hình đáy mặt cắt chọn hình Hình 11: Biến đổi địa hình đáy mặt cắt Hình 12: Biến đổi địa hình đáy mặt cắt Hình 13: Biến đổi địa hình đáy mặt cắt Hình 14: Biến đổi địa hình đáy mặt cắt Hình 14: Biến đổi địa hình đáy mặt đáy MC5 Mặt cắt biểu diễn thay đổi cao độ đáy dọc sông vẽ với 5000m phía cửa sơng 2000m đáy biển phía ngồi cửa sơng Kết hình vẽ 14 cho thấy đoạn từ cửa sông vào khoảng 2000m đáy sông bị xói khoảng 106 1.5 đến 2.0m đoạn độ rộng sông thay đổi từ 700m đến 1100m từ cửa vào Khoảng 3000m nối tiếp lên thượng lưu đáy sông bồi khoảng 1m đoạn độ rộng sông phổ biến khoảng 2000m có chỗ lên tới 2500m Mặt cắt nằm phía cửa sơng, độ rộng mặt cắt lớn lần mặt cắt 2, đoạn phía cửa sơng bồi Tuy nhiên tượng bồi xảy từ trục động lực sang phía bờ huyện Duy Xun, phía bờ thuộc thành phố Hội An lại có tượng xói mạnh KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ ĐẶC BIỆT (11/2013) Mặt cắt bị xói tồn với lớp xói lên tới – 2m, hay nói cách khác độ rộng sơng khu vực gần cửa bị co hẹp đột ngột nên lưu tốc dòng chảy lũ có lúc lên tới 3m/s, nguyên nhân sinh xói sâu trên đáy sơng, tồn lượng bùn cát mang ngồi biển phía trước cửa sơng Từ hình vẽ 14 15 thấy mực nước biển dâng lên thêm 50cm hay 100cm đoạn sơng có độ rộng lớn, lớp bùn cát bồi lắng dày mực nước biển cao hơn, hay lý giải độ dốc mặt nước nhỏ lưu lượng biên không thay đổi Hiện tượng không rõ ràng đoạn sông sát cửa với mức co hẹp đáng kể, lưu tốc có giảm đi, lớn nhiều lưu tốc khơng lắng dẫn tới độ dày xói khơng khác nhiều kịch Đối với mặt cắt hình vẽ 13 15 cho kết đáy biển phía ngồi cửa sơng bồi sau trận lũ, mức độ bồi mặt cắt khác Sự khác phụ thuộc vào độ lớn hướng sóng mực nước biển, mức độ bồi đáy biển phía ngồi cửa khơng b.2 Mơ hình lũ 1999 có xét đến điều tiết hồ chứa ảnh hưởng biến đổi khí hậu Các kết trung gian mơ bao gồm q trình lũ xem xét tới điều tiết, cắt lũ hồ chứa thượng lưu Kết xác định lũ thiết kế điều tiết lấy từ nghiên cứu viện quy hoạch thủy lợi Việt nam thực [4] cầu Câu Lâu, trường mực nước biên lỏng phía ngồi biển, trường sóng trường lưu tốc Khi có kết đó, tiến hành mơ biến đổi địa hình đáy cho kịch thuộc nhóm Dưới xin trích số kết đánh giá biến đổi địa hình đáy mặt cắt điển hình chọn hình Biến đổi địa hình đáy MC1 -1 200 -2 400 600 800 1000 Hiện trạng -3 -5 -6 -7 -8 (m) Hình 16: Biến đổi địa hình đáy MC1 200 400 600 800 1000 1200 1400 400 500 Kịch B2 Kịch B3 (m) Biến đổi địa hình đáy MC4 1600 1800 2000 -1 -1 -2 -2 1000 2000 3000 4000 5000 6000 7000 -3 -3 (m) -4 -4 Hiện trạng -5 Kịch B1 Kịch B2 -6 Kịch B3 -5 -7 -8 300 Hiện trạng Kịch B1 (m) 200 Hình 17: Biến đổi địa hình đáy MC2 Biến đổi địa hình đáy MC3 100 (m) (m) Kịch B1 Kịch B2 Kịch B3 -4 Biến đổi địa hình đáy MC2 -1 -2 -3 -4 -5 -6 -7 -8 -9 -10 Hiện trạng Kịch B1 -6 Kịch B2 -7 Kịch B3 -8 (m) Hình 18: Biến đổi địa hình đáy MC3 -9 (m) Hình 19: Biến đổi địa hình đáy MC4 Về diễn biến đáy sông, đáy biển mặt cắt chọn giống với nhóm kịch A, lớp bồi mặt cắt 1, 3, và lớp xói mặt cắt nhỏ trường hợp kịch tương ứng nhóm kịch A Có thể lý giải sai khác đỉnh lũ sau điều KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ ĐẶC BIỆT (11/2013) 107 tiết qua hồ chứa nhỏ đỉnh lũ nguyên dạng tới 1000m3/s, điều kiện biên nguyên, dẫn tới vận tốc dòng chảy mặt cắt phía cửa sông nhỏ lũ chưa bị điều tiết, điều dẫn tới tình trạng xói đoạn gần cửa hay bồi đoạn nhỏ so với lũ nguyên dạng Đáy biển phía ngồi cửa sơng bồi bùn cát mang từ sông qui luật bồi phù hợp với điều kiện thủy động lực biển, mực nước biển dâng lên, độ lớn sóng tăng lên, hướng sóng khơng thay đổi dẫn tới lớp bùn cát lắng động tạo thành bar ngầm phía ngồi mỏng Cũng phải nhấn mạnh với kết mô cho trận lũ năm 1999 dù trường hợp tự nhiên hay điều tiết khó lý giải cách rõ ràng qui luật bồi xói vùng cửa sơng phía Tính phù hợp chúng trường hợp cụ thể lũ với điều kiện thủy động lực biển xét, thay đổi điều kiện biên dưới, chẳng hạn tổ hợp đỉnh lũ với chân triều khác với tổ hợp đỉnh lũ với đỉnh triều b.3 Mơ hình dòng chảy mùa cạn có xét đến ảnh hưởng biến đổi khí hậu Thời gian mơ từ đến 18/VII/1999 thời kỳ dòng chảy sông gần nhỏ nhất, số liệu mô MIKE11 Câu Lâu thời kỳ có trị số lưu lượng 50m3/s Chính vậy, vận tốc dòng chảy khu vực cửa nhỏ (hình vẽ 21) phổ biến 0.20 (m/s), lớn 0.4m/s Hình 21: Vận tốc dòng chảy tổng cộng kịch C1, C2, C3 108 KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ ĐẶC BIỆT (11/2013) Hình 22: Địa hình đáy – Kịch C1 Hình 23: Địa hình đáy – Kịch C2 Hình 24: Địa hình đáy – Kịch C3 Các hình vẽ 22 đến 24 biểu diễn địa hình đáy sau mô theo kịch C1, C2 C3 Từ kết mơ phỏng, thấy biến đổi địa hình đáy khơng lớn phổ biến khoảng cm lớn khoảng 12 cm Phía biển doi cát ngầm bị phá vỡ mực nước biển dâng lên lý sóng dòng chảy ven bờ tăng đáng kể tác động vào dải cát ngầm Khu vực phía cửa sơng tượng bồi lớn hơn, doi cát lên khỏi mặt nước mở rộng mực nước biển tăng lên Nếu tính thời gian mô khoảng 20 ngày, mà độ cao bồi gia tăng thêm – 10 cm nói số có ý nghĩa gây nên tình trạng “đóng cửa” bùn cát từ ngồi cửa dòng dọc bờ ngang bờ mang vào thời kỳ nước từ thượng lưu nhỏ chịu tác động chủ yếu gió mùa Tây Nam Kết luận So sánh kết mô mơ hình MIKE21FM số liệu đo đạc thực tế khẳng định mơ mơ hình tốn cho kết phù hợp với thực tế Tuy nhiên, với lực máy tính khơng thể mơ cho chu kỳ đại diện (ở nước ta chu kỳ đại diện năm trùng với chu kỳ khí hậu), nên thời gian mơ tốt cho điều kiện cực hạn, chẳng hạn trận lũ lớn, xả lũ từ hồ chứa kết hợp với điều kiện sóng, triều khác phía biển Việc mô để đánh giá xu diễn biến đáy sơng đáy biển dài hạn thực với trường sóng thủy triều đặc trưng KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ ĐẶC BIỆT (11/2013) 109 điều kiện dòng chảy, bùn cát từ sơng với tỉ lệ tương tự khác Cũng mơ biến đổi địa hình mơ hình hình thái Có thể khẳng định việc ứng dụng MIKE21FM hồn tồn khả thi kết mơ giúp định hướng giải pháp cơng trình nhằm mục đích đề ngăn chặn bùn cát gây bồi cửa sơng, giảm sóng gió tạo luồng giao thơng thủy đủ độ cơng trình lũ, đẩy cát bồi phía cửa ngồi xa ổn định cửa sông Tài liệu tham khảo [1] Nghiên cứu dự báo phòng chống sạt lở bờ sơng hệ thống sông miền Trung; Đề tài khoa học công nghệ cấp nhà nước, 2001 [2] Nghiên cứu giải pháp phòng tránh xói lở bồi lấp cửa sơng Thu Bồn – Vũ gia Đề tài khoa học công nghệ cấp Bộ, 2003 [3] Nghiên cứu đề xuất giải pháp ổn định cửa sông ven biển miền Trung; Đề tài khoa học công nghệ cấp nhà nước, 2009 [4] DHI - Mike11 Introduction and tutorial, User manual and Reference manual, 2008 Abstract SIMULATION OF SEDIMEN TRANSPORT AND BED EVOLUTION AT THU BON ESTUARY Paper presents the simulated results on river and sea bed evolution by using MIKE21_FM Couple and later simulate bed elevation variation under the different scenarios of typical flood events, flood release from reservoirs, sedimentation in the river as well as the different water level and wave spectrum in the sea The simulated results show that it is a good agreement between observed data and simulated results This is a river and sea bed erosion after big flood in the river and sedimentation in front of the river mouth during dry season and south west monsoon predominately Simulated results support decision makers to select appropriate solutions to keep river mouth stably for flood release, navigation as well as protection against erosion and sedimentation for study area Key words: sea bed, simulation, deep sea wave, sediment transport, MIKE21-FM Couple Người phản biện: TS Trần Thanh Tùng 110 BBT nhận bài: Phản biện xong: 25/10/2013 7/11/2013 KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ ĐẶC BIỆT (11/2013) ... mơ biến đổi địa hình đáy cho kịch thu c nhóm Dưới xin trích kết đánh giá biến đổi địa hình đáy mặt cắt chọn hình Hình 11: Biến đổi địa hình đáy mặt cắt Hình 12: Biến đổi địa hình đáy mặt cắt Hình. .. Kết mô vận chuyển bùn cát biến đổi địa hình đáy KHOA HỌC KỸ THU T THỦY LỢI VÀ MÔI TRƯỜNG - SỐ ĐẶC BIỆT (11/2013) 105 Mục tiêu kịch phần đánh giá biến động địa hình đáy sơng địa hình đáy biển khu. .. Hình 13: Biến đổi địa hình đáy mặt cắt Hình 14: Biến đổi địa hình đáy mặt cắt Hình 14: Biến đổi địa hình đáy mặt đáy MC5 Mặt cắt biểu diễn thay đổi cao độ đáy dọc sơng vẽ với 5000m phía cửa sơng