Đang tải... (xem toàn văn)
Bằng cách tận dụng thế mạnh của mô hình toán mã nguồn mỡ Telemac, nghiên cứu đã thay đổi cấu trúc đáy phù hợp với đặc tính vật lý lòng dẫn thông qua áp đặt giá trị (ES1, ES2) = f(f1, f2, n, ES) để mô phỏng vận chuyển bùn cát theo đặc tính hỗn hợp (mixed sediment) trong sông và rời (non-cohesive sediment) phía ngoài biển áp dụng cho cửa sông Soài Rạp.
BÀI BÁO KHOA HỌC THIẾT LẬP MƠ HÌNH VẬN CHUYỂN BÙN CÁT HỖN HỢP THEO KHÔNG GIAN MÔ PHỎNG DIỄN BIẾN HÌNH THÁI LỊNG DẪN TẠI CỬA SƠNG SỒI RẠP BẰNG TELEMAC-2D Lê Ngọc Anh1,2, Hồng Trung Thống2, Ngơ Bảo Châu2 Tóm tắt: Thiết lập mơ hình vận chuyển bùn cát có đặc tính bùn cát biến đổi phức tạp theo không gian vấn đề cần đặt nghiên cứu động lực bùn cát hình thái vùng cửa sông Bằng cách tận dụng mạnh mô hình tốn mã nguồn mỡ Telemac, nghiên cứu thay đổi cấu trúc đáy phù hợp với đặc tính vật lý lịng dẫn thơng qua áp đặt giá trị (ES1, ES2) = f(f1, f2, n, ES) để mô vận chuyển bùn cát theo đặc tính hỗn hợp (mixed sediment) sơng rời (non-cohesive sediment) phía ngồi biển áp dụng cho cửa sơng Sồi Rạp Mơ hình thiết lập có độ tin cậy cao với sai số hàm lượng bùn cát lơ lửng MSE từ 0,001 - 0,013 Kết mô sau năm cho thấy, cửa Sồi Rạp, dịng chủ lưu xuất hố xói có chiều sâu xói từ 0,3 - 1,0m; lượng xói dịng chảy mang xa biển hình thành nên lưỡi cát có chiều cao bồi thêm từ 0,3 - 0,4m Từ khóa: Cửa sơng Sồi Rạp, Telemac, Sisyphe, mơ hình vận chuyển bùn cát hỗn hợp GIỚI THIỆU * Nghiên cứu q trình biến đổi hình thái cửa sơng thông tin cần thiết cho việc hoạch định chiến lược quản lý tài ngun nước Phương pháp mơ hình toán số sử dụng phổ biến nghiên cứu gần khả xây dựng kịch phong phú, chi phí thấp,thời gian mơ nhanh, phạm vi tính tốn rộng so với phương pháp mơ hình vật lý vốn khó đạt điều kiện tương tự bùn cát Nhiều nghiên cứu mô diễn biến lịng dẫn vùng HLSĐN mơ hình tốn số với đặc tính bùn cát rời (Nguyễn Đức Vượng, et al., 2005; Nguyễn Thế Biên, et al., 2012), đặc tính bùn dính (Nguyễn Thị Bảy, et al., 2006; Nguyễn Thị Bảy, et al., 2012; Lê Mạnh Hùng, et al., 2015; Trần Bá Hoằng, et al., 2014; Bùi Trọng Vinh, 2016) tiền đề cho nghiên cứu sau Tuy nhiên, đặc tính bùn cát khu vực thường không giống nhau, khu vực sông cửa sơng thể đặc tính bùn dính khu vực ngồi biển lại mang đặc tính rời Chính thế, việc mơ q trình vận chuyển bùn cát khu vực có đặc tính bùn cát biến đổi phức tạp theo không gian cần linh hoạt Như vậy, mục tiêu nghiên cứu thiết lập mơ hình vận chuyển bùn cát hỗn hợp thay đổi theo không gian dựa việc áp đặt chiều dày cát bùn theo quan hệ (ES1, ES2) = f(f1, f2, ES, n) mơ hình mã nguồn mở Telemac-2d ứng dụng khu vực cửa sông Sồi Rạp Cửa sơng Sồi Rạp nằm vùng hạ lưu hệ thống sơng Đồng Nai (HLSĐN) có diện tích lưu vực 40.700 km2 Sơng Sồi Rạp nhánh hệ thống HLSĐN, dài 45 km qua xã Phú Xuân (H.Nhà Bè), xã Bình Khánh (H.Cần Giờ) đổ biển Đơng cửa Sồi Rạp có chiều rộng ngang sơng khoảng 2.420 m, cao trình đáy sơng biến đổi từ -32,0m đến -8,0 m (Hình 1) PHƯƠNG PHÁP NGHIÊN CỨU 2.1 Lý thuyết vận chuyển bùn cát hỗn hợp 2.1.2 Phương trình vận chuyển bùn cát hai chiều bùn cát lơ lửng trung bình theo phương thẳng đứng C = C(x, y, t) NCS Trường ĐH Bách Khoa Tp.HCM Trường Đại học Tài nguyên & Môi trường Tp.HCM KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ 69 (6/2020) hCk hUCk hVCk Ck Ck k k hs hs E D t x y x x y y E k D k Zref s C eqk C refk Hình Cửa sơng Sồi Rạp 175 Hệ số k đại diện cho thành phần hạt k = 1, đại diện cho hạt cát bùn h = Zs – Zf ≈ Zs – Zref chiều sâu nước, giả thiết chiều dày lớp bùn cát đáy mỏng; (U,V) vận tốc trung bình theo phương x, y; E: suất xói lở (kg/m2/s); D: suất bồi tụ (kg/m2/s), (E – D) lượng trữ trầm tích lơ lửng; Ceq nồng độ bùn cát trạng thái cân sát đáy; Cref nồng độ bùn cát sát đáy Sự phân phối hàm lượng bùn cát lơ lửng tuân theo quy luật Rouse: R zh a C z CZref z ah với R ws số Rouse u* Sự thay đổi đáy lịng dẫn tính tốn dựa cân khối lượng lưu lượng bồi/xói z 1 b E D z Zref t Trong đó: hệ số độ rổng, zb cao trình đáy Ứng xử mô bùn cát hỗn hợp Tỷ lệ phần trăm khối lượng bùn lớp f j M s2j / M s1j M s2j sử dụng để xác định ứng suất đáy tới hạn trung bình lớp ce lưu lượng bồi/xói tương ứng j Ứng suất đáy tới hạn trung bình Với f j 30% (cát chiếm tỷ lệ lớn), vận chuyển bùn cát dựa đặc tính rời (non1 ce (với conhesive sediment), j ce j ce ứng suất tới hạn cát) Với f j 50% (bùn chiếm tỷ lệ lớn), vận chuyển bùn cát dựa đặc tính dính bùn ce (với (conhosive sediment), j ce j ce ứng suất tới hạn bùn) Ứng suất xói tới hạn xác định thí nghiệm Với 30% f2 j 50% thể đặc tính bùn hỗn hợp ứng suất tới hạn trung bình tính theo nội suy tuyến tính sau: f 0,3 ce2 ce1 ce ce j j j j j 0,5 0,3 Lượng xói trung bình lớp 176 Lưu lượng xói trung bình lớp E j xác định sau: Với f j 30% , lưu lượng xói tính: w Ceq f1 ;( b ce ) j j s Ej E 0;( b ce j ) Với f1j phần trăm thể tích cát chứa j lớp, b ứng suất tổng đáy tính theo b 0,5C f U V ; Cf hệ số ma sát tổng Với f j 50% , lưu lượng xói tính: M b 1 ; ( b ce ) j E j E 2j ce j 0; ( b ce j ) Với M(kg/m2/s) số xói KronePartheniades Với 30% f j 50% , lưu lượng xói tính dựa nguyên tắc nội suy tuyến tính tính sau: Ej E j f 2j 0, E 2j E 1j 0, 0, Lượng bồi lớp Đối với cát: D1 w s1.T2 Với T2 tỷ số hàm lượng bùn cát sát đáy hàm lượng bùn cát trung bình tính tính theo quy luật Rouse Đối với bùn: T 2 D w s 1 cr1 u*mud Với u*crmud (m/s) vận tốc tới hạn bồi sát đáy bùn T1 b / 2.2 Xây dựng công thức quan hệ (ES1, ES2) = f(f1, f2, n, ES) Gọi M1, M2 khối lượng tương ứng cát bùn khối đất; f1, f2 tỷ lệ khối lượng tương ứng cát bùn; n: độ rỗng hạt cát Do kích thước hạt bùn nhỏ nên bỏ qua độ rỗng hạt bùn (Hình 2) KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ 69 (6/2020) Hình Cấu tạo lớp bùn cát hỗn hợp Ta có: M 1 n ES1 ; M C.ES2 : khối lượng riêng cát (2.650 kg/m3); C: hàm lượng bùn (kg/m3) M1 M2 ; f2 M1 M M1 M M f C ES2 f M1 f1 1 n ES1 f1 f1 Cho trước chiều dày lớp đất ES ta chứng minh quan hệ sau: C.ES2 f C.ES f 1 n ES1 f1 1 n ES ES f1 ES2 2.3.1 Mơ hình Telemac Hệ thống mơ hình Telemac Phịng Thí Nghiệm Thủy Lực Mơi Trường Quốc Gia thuộc trung tâm quốc gia nghiên cứu Thủy lực Điện Lực Pháp (EDF) phát triển từ năm 1987 Đây mơ hình có khả mơ trình thủy động lực, vận chuyển bùn cát hình thái sơng/biển, chất lượng nước ứng dụng cho hồ chứa, sông, cửa sông, ven biển, đại dương Telemac mơ hình mã nguồn mở, hổ trợ tính tốn song song hồn tồn miễn phí Trong nghiên cứu này, mơ hình Telemac2d – Sisyphe – Tomawac kết nối để mơ q trình vận chuyển bùn cát; đó, mơ hình Telemac-2d đóng vai trị hạt nhân trao đổi thông tin thủy lực với mơ hình Sisyphe Tomawac 2.3.2 Sơ đồ mạng lưới thủy lực Miền tính bao gồm sơng Sồi Rạp sông Mekong xây dựng 154.459 phần tử tam giác, với 84.903 nút Kích thước lưới nhỏ từ 22 – 100m áp dụng cho sông vùng cửa sông; vùng tiếp giáp cửa sông biển Đơng có kích thước lưới từ 200 – 600m, ngồi biển Đơng có kích thước 800 – 12.700m Hình thể sơ đồ mạng lưới thủy lực f1 f2 1 n ES C f1 f1 f 1 n C f1 ES C f1 f1 f 1 n Hình 3: Minh họa thay đổi đặc tính bùn cát theo khơng gian ứng dụng cửa sơng Sồi Rạp ES1 ES ES2 Hình Sơ đồ mạng lưới thủy lực Hình Đặc tính bùn cát thay đổi theo khơng gian 2.3 Thiết lập mơ hình 2.3.3 Telemac-2d Lưu lượng biên thượng lưu lấy trạm Cần Thơ (sông Hậu), trạm Mỹ Thuận (sông Tiền), Hồ Trị An (sông Đồng Nai) Hồ Dầu Tiếng (sông Sài Gịn) Biên hạ lưu cách cửa Sồi Rạp 120km gồm mực nước (Z) thành phần vận tốc trung bình (U, V) theo hai phương X, Y áp đặt biên hở biển Các giá trị Z, U, V xem tổng nhiều thành phần sóng khác biểu diễn theo phương trình sau: KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ 69 (6/2020) 177 t Fi M, t fi t AFi M cos 2 Fi M ui0 vi t Ti h hi zb zmean ;U Ui ;V Vi i i i Trong đó: F M , t = độ sâu nước (h), thành phần vận tốc theo phương ngang U, V, Zb cao trình đáy, Zmean mực nước biển hiệu chỉnh, f i t hệ số suy biến, ui0 , vi t pha thời điểm t = t AFi M , Fi M số điều hòa đọc vào từ liệu triều (TPXO) Qui luận ma sát Manning áp dụng với hệ số ma sát từ 0,018 – 0,025 2.3.4 Sisyphe Thành phần hạt bùn cát sơng Sồi Rạp khơng đồng phức tạp bao gồm từ cuội sỏi đến bùn, sét có đường kính chủ yếu loại hạt có đường kính d = 0,01 – 10,0mm (Hoàng Văn Huân, et al., 2005) Trong nghiên cứu này, cấu trúc đáy thiết lập dạng bùn cát hổn hợp gồm hai thành phần hạt cát (d = 0,32.10-3m) bùn (d = 0,025.10-3m), chiều dày lớp cát (ES1) lớp bùn (ES2) lớp thay đổi theo không gian cách thay đổi lại hàm “INIT_COMPO_COH” dựa quan hệ (ES1, ES2) = f(f1, f2, n, ES) cho khu vực sông thể đặc tính bùn hỗn hợp (30% ≤ f2 ≤ 50%), cịn khu vực biển Đơng thể đặc tính bùn rời (f2 ≤ 30%) Hàm lượng bùn cát lớp C1 = 150; C2 = 250 (kg/m3) Ứng suất tới hạn bồi b 1000 Pa ; ứng suất tới hạn xói lớp ce1 0, 021( N / m ) ce2 0, 25( N / m ) Theo liệu đo đạc (McLachlan R L., et al., 2017) sông MeKong, vận tốc lắng biến đổi từ 0,01 – 1,0 mm/s phụ thuộc vào thời kỳ triều cao thấp dòng chảy theo mùa sông Trong nghiên cứu này, vận tốc lắng hạt chọn vs = 0,15 mm/s vs = 0,035 mm/s tương ứng với thành phần cát bùn Active layer thiết lập 0,05 m bước thời gian mô 45s 2.4 Hiệu chỉnh kiểm định mô hình trạm đo vận tốc dịng chảy hàm lượng bùn cát sơng Sồi Rạp dùng để hiệu chỉnh kiểm định mơ hình, sơ đồ vị trí Hình 178 Hình Sơ đồ tuyến đo vận tốc dòng chảy hàm lượng bùn cát lơ lửng sơng Sồi Rạp 2.4.1 Vận tốc dịng chảy Vận tốc dòng chảy yếu tố quan trọng ảnh hưởng trực tiếp đến kết mô q trình vận chuyển bùn cát sơng cửa sông Kết hiệu chỉnh cho kết tốt, sai số quân phương (MSE) từ 0,08 – 0,15, ME từ 0,05 – 0,22 MAE từ 0,15 – 0,32 Sai số vận tốc mộ số vị trí cịn chưa tốt trạm TV5 (sơng Đồng Nai xa cửa Soài Rạp), MAE 0,322 (xem Bảng 1) 2.4.2 Hàm lượng bùn cát (SSC) SSC hiệu chỉnh mùa kiệt 04/2014 kiểm định lại mùa lũ 11/2013 Kết hiệu chỉnh kiểm định trạm TV1, TV2, TV3, TV4, TV5, TV6 thể Bảng , Hình Kết hiệu chỉnh kiểm định cho thấy sai số quân phương nhỏ MSE từ 0,001 – 0,013, sai số trung bình MAE từ 0,014 – 0,128 Như vậy, kết mô tốt cho thấy mơ hình có độ tin cậy cao Bảng Kết hiệu chỉnh vận tốc trạm đo sơng Sồi Rạp Vị trí TV1 TV2 TV3 TV4 TV5 TV6 Sai số Max Min 0,078 0,001 0,046 0,011 0,083 0,043 0,343 0,010 0,409 0,032 0,091 0,001 ME 0,066 0,055 0,105 0,205 0,226 0,127 MAE 0,152 0,154 0,308 0,217 0,322 0,294 MSE 0,118 0,079 0,086 0,104 0,151 0,085 ME: sai số tuyệt đối trung bình, MAE: sai số trung bình, MSE: sai số quân phương KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ 69 (6/2020) Bảng Kết hiệu chỉnh kiểm định hàm lượng bùn cát Vị trí TV1 TV2 TV3 TV4 TV5 TV6 TV1 TV2 TV3 TV4 TV5 TV6 Sai số ME MAE Max Min Hiệu chỉnh (04/2014) 0,076 0,147 0,128 0,128 0,020 0,101 0,041 0,055 0,030 0,171 0,098 0,105 0,003 0,020 0,003 0,027 0,004 0,003 0,018 0,031 0,007 0,028 0,009 0,014 Kiểm định (11/2013) 0,046 0,373 0,068 0,099 0,010 0,085 0,040 0,044 0,012 0,148 0,011 0,059 0,047 0,099 0,047 0,054 0,010 0,011 0,019 0,033 0,019 0,135 0,005 0,018 MSE 0,003 0,002 0,005 0,001 0,001 0,000 0,013 0,001 0,005 0,001 0,001 0,001 tương tác động lực sông biển Dịng chảy sơng (nước ngọt) hình thành từ trình thủy văn phụ thuộc vào yếu tố mặt đệm lưu vực, chảy theo chiều biển gọi dòng dương Dòng chảy từ biển chủ yếu dòng triều (nước mặn) chảy theo hai chiều: (1) chảy từ biển vào hệ thống sông bên triều lên gọi dịng âm, (2) chảy từ sơng biển gọi dòng dương triều rút Khi triều rút, lượng dòng chảy dương tổng nước nước mặn tích sơng nên vận tốc dịng chảy theo chiều dương lớn so với chiều âm Hình cho thấy dịng chảy theo chiều dương (ra biển) có giá trị lớn tần suất xuất nhiều so với chiều ngược lại Tại cửa Soài Rạp (P1, P2, P3, P4), vận tốc dòng chảy lớn đạt 1,65 m/s, vận tốc trung bình đạt 0,52 m/s Hình Vận tốc dịng chảy cửa sơng Sồi Rạp Hình Hiệu chỉnh kiểm định hàm lượng bùn cát trạm TV1, TV2, TV3, TV4, TV5, TV6 KẾT QUẢ NGHIÊN CỨU Dòng chảy Dòng chảy vào, cửa sơng Sồi Rạp chịu Vận tốc dòng chảy kỳ triều cường lớn so với kỳ triều (Hình 8) Do tác động nghịch dịng chảy từ sơng với dịng triều nên triều lên, vận tốc dòng chảy mùa lũ nhỏ so với mùa kiệt; ngược lại triều xuống, vận tốc dòng chảy mùa lũ lớn so với mùa kiệt (Hình 8) Hàm lượng bùn cát SSC thời kỳ triều cường lớn so với kỳ triều Trong mùa kiệt, ảnh hưởng triều mạnh, SSC có giá trị lớn từ 0,2÷0,50 kg/m3 vị trí cửa sơng từ P3 đến P4; sang mùa lũ bùn cát dịng chảy sơng đẩy xa cửa sông SSC đạt giá trị lớn từ 0,05÷0,60 kg/m3 phạm vi từ P1 đến P3 (Hình 10) Trong kỳ triều cường, vị trí cửa sơng (P4) hàm lượng bùn cát đạt giá trị lớn 1/7 sườn triều xuống đạt giá trị nhỏ khoảng 1/2 sườn triều lên; Trong kỳ triều kém, hàm lượng bùn cát lớn chân triều Trong kỳ triều KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ 69 (6/2020) 179 kém, động lực triều suy yếu tương tác dịng chảy sơng biển khác cân bằng, giao động SSC không cho thấy phân biển rõ chân đỉnh (Hình 9) Hình thái lịng dẫn Hình 11 mơ tả thay đổi đáy khu vực cửa sơng Sồi Rạp năm Kết mơ cho thấy có hình thành lưỡi cát vùng cửa sơng theo mùa Mùa kiệt giai đoạn tích lũy bùn cát khu vực thượng lưu dòng chảy mang dần biển, tốc độ dịch chuyển chậm Sang đến mùa lũ, vận tốc dịng chảy sơng mạnh, kết hợp với dòng triều tạo thành dòng kết hợp với vận tốc lớn gây xói mạnh dịng chủ lưu đẩy khối bùn cát phía cửa sông biển Kết mô sau năm cho thấy, dòng chủ lưu xuất hố xói có chiều sâu xói từ 0,3 - 1,0m; lượng xói dịng chảy mang xa biển hình thành nên lưỡi cát có chiều cao bồi thêm từ 0,3 - 0,4m Kết luận thảo luận Mô vận chuyển bùn cát hỗn hợp theo hướng áp đặt linh hoạt đặc tính bùn cát từ sông biển cách tiếp cận phù hợp tốn có miền tính rộng với đặc tính bùn cát biến đổi phức tạp theo khơng gian Mơ hình hiệu chỉnh kiểm định đầy đủ yếu tố thủy lực bùn cát, có độ tin cậy cao thơng qua việc lựa chọn thơng số vật lý thơng số tốn học theo nguyên tắc: (1) thông số vật lý độ nhám thủy lực, kích thước hạt, vận tốc lắng, ứng suất đáy tới hạn… cần phải đảm bảo không sai khác lớn so với thực nghiệm; (2) thơng số tốn số cần đảm bảo cho ổn định sơ đồ toán số Kết cho thấy, vận chuyển bùn cát cửa sơng Sồi Rạp chịu chi phối triều Biển Đông chế độ dòng chảy theo mùa Trong mùa kiệt, lượng bùn cát với lưu lượng dòng chảy nhỏ nên phần lớn bùn cát tích trữ đây, chúng mang dần phía cửa sơng vào mùa lũ với tăng lên dòng chảy lũ Kết hình thành lưỡi cát phía cửa sông lan dần đến biển Giao động hàm lượng bùn cát vùng cửa sông kỳ triều theo mùa biến đổi hình thái cửa sơng Sồi Rạp năm cho thấy phù hợp với quy luật vận chuyển bùn cát chung cửa sông chịu ảnh hưởng thủy triều yếu tố mùa sâu sắc 180 Hình Vận tốc dịng chảy sơng Đồng Nai từ điểm P1 – P8 Hình SSC số vị trí cửa sông mùa kiệt lũ vào kỳ triều Hình 10 SSC kỳ triều cường vận tốc dịng chảy P4 (tại cửa Sồi Rạp) đạt giá trị lớn (a, b) vào mùa kiệt, (c, d) vào mùa lũ Hình 11 Quá trình biến đổi hình thái lòng dẫn sau năm KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ 69 (6/2020) TÀI LIỆU THAM KHẢO Bùi Trọng Vinh (2016) 'Xói lở bờ biển gị Cơng Đơng - Tiền Giang', Tạp chí phát triển KH&CN, 19(K1), pp 59-69 Hoàng Văn Huân, Lê Ngọc Bích, et al (2005), Nghiên cứu diễn biến lịng sơng, hình thái sơng loại lịng dẫn hạ du sơng Đồng Nai - Sài Gịn, Tp Hồ Chí Minh Lê Mạnh Hùng, Đinh Công Sản, et al (2015), Nghiên cứu khoa học liên quan đến dự án chỉnh trị luồng, đánh giá sa bồi sau nạo vét, Viện Khoa học Thủy lợi miền Nam Nguyễn Đức Vượng, Phạm Trung, et al (2005), Đề tài KC-08-29: "Nghiên cứu đề xuất giải pháp khoa học công nghệ để ổn định lịng dẫn hạ du hệ thống sơng Đồng Nai - Sài Gòn phục vụ cho phát triển kinh tế - xã hội vùng Đông Nam bộ", Viện Khoa Học Thủy Lợi Miền Nam Nguyễn Thế Biên, Hoàng Văn Huân, et al (2012), Nghiên cứu, điều tra, khảo sát, đánh giá ảnh hưởng đề giải pháp khắc phục, hạn chế sạt lở bờ khai thác cát địa bàn Thành phố Hồ Chí Minh, Viện Khoa học Thủy lợi miền Nam Nguyễn Thị Bảy, Mạnh Quỳnh Trang (2006) 'Mơ hình tốn tính chuyển tải bùn cát kết dính vùng ven biển', Science & Technology Development, 9(4), pp 31-39 Nguyễn Thị Bảy, Phan Văn Hoặc, et al (2012), Nghiên cứu đánh giá trình thay đổi luồng lạch diễn biến thuỷ văn nạo vét cửa sơng Sồi Rạp phục vụ giao thơng thủy, Viện Khí tượng thủy hải văn Mơi trường Trần Bá Hoằng, Nguyễn Duy Khang (2014), Nghiên cứu biến động chế độ thủy thạch động lực vùng cửa sông ven biển chịu tác động Dự án đê biển Vũng Tàu - Gị Cơng, Trung tâm NC Chỉnh trị sơng Phòng chống thiên tai - Viện Khoa học Thủy lợi Miền Nam McLachlan R L., Ogston A S., Allison M A (2017) 'Implications of tidally-varying bed stress and intermittent estuarine stratification on fine-sediment dynamics through the Mekong’s tidal river to estuarine reach', Continental Shelf Research, 147, pp 27-37 Abstract: ESTABLISHING A MIXED SEDIMENT TRANSPORT MODEL BASED ON SIMULATION SPACE OF MORPHOLOGICAL CHANGES OF THE MAINSTREAM AT SOAI RAP ESTUARY BY TELEMAC-2D Establishing a sediment transport model with spatial complex variation of characteristics of sediment is a matter for research in the sediment dynamics and estuary morphology By taking advantage of open source mathematical model Telemac, the study changed the riverbed structure suitable for physical characteristic by imposing the value (ES1, ES2) = f (f1, f2, n, ES) to transport simulation mixed sediment in the river and non-cohesive sediment outside the sea apply to Soai Rap estuary The model was established with high reliability and errors in the content of suspended sediment MSE were from 0.001 0.013 Simulation results after year show that, at the Soài Rạp estuary, in the middle of the mainstream there are erosion pits with the depth of erosion from 0.3 to 1.0 m; that erosion is carried far from the sea by the flow and forming sand dunes with an accretion height of 0.3 - 0.4m Keywords: Soai Rap estuary, Telemac, Sisyphe, mixed sediment transport model Ngày nhận bài: 06/5/2020 Ngày chấp nhận đăng: 30/6/2020 KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ 69 (6/2020) 181 ... Kết hình thành lưỡi cát phía cửa sông lan dần đến biển Giao động hàm lượng bùn cát vùng cửa sông kỳ triều theo mùa biến đổi hình thái cửa sơng Sồi Rạp năm cho thấy phù hợp với quy luật vận chuyển. .. biển hình thành nên lưỡi cát có chiều cao bồi thêm từ 0,3 - 0,4m Kết luận thảo luận Mô vận chuyển bùn cát hỗn hợp theo hướng áp đặt linh hoạt đặc tính bùn cát từ sông biển cách tiếp cận phù hợp. .. cân bằng, giao động SSC không cho thấy phân biển rõ chân đỉnh (Hình 9) Hình thái lịng dẫn Hình 11 mơ tả thay đổi đáy khu vực cửa sơng Sồi Rạp năm Kết mô cho thấy có hình thành lưỡi cát vùng cửa