Nghiên Cứu Ứng Dụng Mô Hình 3D Vào Tính Toán Dòng Chảy Trong Sông Soài Rạp

Luận văn thạc sĩ kỹ thuật xây dựng nghiên cứu ứng dụng mô hình 3D trong tính toán dòng chảy sông và ven biển, mang lại giải pháp hiệu quả.

• Tác giả: Nguyễn Thu Hà
• Người hướng dẫn: PGS. Huỳnh Thanh Sơn
• Trường học: Đại học Quốc gia TP.HCM
• Chuyên ngành: Xây dựng công trình thủy
• Đề tài: Nghiên Cứu Ứng Dụng Mô Hình 3D Vào Tính Toán Dòng Chảy Sông Soài Rạp
• Loại tài liệu: luận văn thạc sĩ
• Năm xuất bản: 2013
• Địa điểm: TP. Hồ Chí Minh

Mục lục chi tiết

• LỜI CÁM ƠN
• TÓM TẮT
• LỜI CAM ĐOAN
• 1. CHƯƠNG 1: MỞ ĐẦU
• 1.1. Đặt vấn đề
• 1.2. Mục đích của luận văn
• 1.3. Các nghiên cứu có liên quan
• 1.3.1. Các nghiên cứu ngoài nước
• 1.3.2. Các nghiên cứu trong nước
• 1.4. Phương pháp nghiên cứu
• 1.5. Các nội dung nghiên cứu của luận văn
• 1.6. Phạm vi nghiên cứu
• 1.7. Cấu trúc luận văn
• 2. CHƯƠNG 2: CƠ SỞ LÝ THUYẾT CỦA MÔ HÌNH ECOMSED
• 2.1. Giới thiệu
• 2.2. Module thủy động lực học
• 2.2.1. Hệ phương trình thủy động lực học
• 2.2.2. Mô hình
• 2.2.3. Điều kiện biên tại mặt thoáng và cây
• 2.2.4. Điều kiện biên hở (Điều kiện biên mực nước)
• 2.2.5. Tọa độ theo phương đứng
• 2.2.6. Kỹ thuật phân tích trường tính toán (Mode Splitting Technique)
• 2.2.7. Chuyển đổi hệ tọa độ cong trực giao
• 2.3. Module chuyển tải bùn cát
• 2.3.1. Phương trình cơ bản
• 2.3.2. Tính ứng suất cắt đáy
• 2.3.3. Sự tái lơ lửng của bùn cát dính
• 2.3.4. Sự lắng của bùn cát dính
• 2.3.5. Mô hình đáy bùn cát dính
• 2.3.6. Sự tái lơ lửng của bùn cát không dính
• 2.3.7. Sự lắng đọng của bùn cát không dính
• 2.3.8. Sự “cứng hóa” đáy của bùn cát không dính
• 2.4. Phương pháp sai phân hữu hạn
• 2.4.1. Sai phân hữu hạn
• 2.4.2. Thông số hóa dưới lưới (subgrid scale parameterisation)
• 2.4.3. Điều kiện ổn định
• 2.4.4. Cấu trúc mô hình
• 3. CHƯƠNG 3: ĐIỀU KIỆN TỰ NHIÊN KHU VỰC NGHIÊN CỨU
• 3.1. Vị trí địa lý
• 3.2. Đặc điểm địa hình, địa mạo
• 3.3. Đặc điểm địa chất
• 3.4. Đặc điểm khí hậu, khí tượng
• 3.5. Đặc điểm thủy văn, bùn cát
• 4. CHƯƠNG 4: ỨNG DỤNG MÔ HÌNH ECOMSED ĐỂ TÍNH TOÁN DÒNG CHẢY ĐOẠN SÔNG SOÀI RẠP
• 4.1. Lưới tính toán
• 4.2. Dữ liệu đầu vào cho module thủy động lực học
• 4.2.1. Thông số thủy động lực học
• 4.2.2. Điều kiện ban đầu
• 4.2.3. Điều kiện biên
• 4.3. Kiểm nghiệm, hiệu chỉnh và phân tích độ nhạy mô hình thủy động lực học
• 4.3.1. Hiệu chỉnh mực nước
• 4.3.2. Hiệu chỉnh vận tốc trung bình
• 4.3.3. Phân tích độ nhạy của mô hình
• 4.4. Kết quả tính toán thủy động lực học
• 4.4.1. Vận tốc tại tầng mặt và tầng đáy
• 4.4.2. Vận tốc phân bố theo độ sâu
• 4.4.3. Ứng suất ma sát đáy
• 4.4.4. Hệ số nhớt rối theo phương ngang và hệ số nhớt rối theo phương đứng
• 5. CHƯƠNG 5: ỨNG DỤNG MÔ HÌNH ECOMSED ĐỂ TÍNH TOÁN CHUYỂN TẢI BÙN CÁT TRONG ĐOẠN SÔNG SOÀI RẠP
• 5.1. Dữ liệu đầu vào của module chuyển tải bùn cát
• 5.1.1. Đáy bùn cát
• 5.2. Kết quả tính toán chuyển tải bùn cát lơ lửng
• 5.2.1. Nồng độ bùn cát lơ lửng theo phương ngang và phương đứng
• 5.2.2. Sự thay đổi cao trình đáy bùn cát
• 6. CHƯƠNG 6: KẾT LUẬN VÀ KIẾN NGHỊ
• 6.1. Kết quả đạt được trong luận văn
• 6.2. Những hạn chế của luận văn
• 6.3. Ưu và nhược điểm của mô hình ECOMSED
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Tổng Quan Nghiên Cứu Mô Hình 3D Dòng Chảy Sông Soài Rạp

Nghiên cứu ứng dụng mô hình 3D vào tính toán dòng chảy sông ngòi, đặc biệt là ở khu vực sông Soài Rạp, ngày càng trở nên quan trọng. Việc này giúp hiểu rõ hơn về mô phỏng dòng chảy sông và các quá trình liên quan. Từ xưa, sông ngòi đã đóng vai trò quan trọng trong sự phát triển kinh tế và xã hội của nhiều quốc gia. Ví dụ, TP.HCM được xây dựng dọc theo sông Sài Gòn, và hệ thống cảng đường thủy đóng vai trò then chốt trong việc phát triển kinh tế khu vực. Theo định hướng phát triển đô thị, khu Cảng Hiệp Phước trên sông Soài Rạp đang được hình thành, đòi hỏi nghiên cứu kỹ lưỡng về dòng chảy và vận chuyển bùn cát. Luồng Soài Rạp có ưu điểm là rộng và ít khúc cong so với luồng Lòng Tàu, hứa hẹn là tuyến đường thủy quan trọng cho tàu trọng tải lớn. Nghiên cứu này sẽ kết hợp đo đạc thực tế và mô hình toán học dòng chảy để đánh giá hiệu quả.

1.1. Tầm Quan Trọng của Nghiên Cứu Dòng Chảy Sông 3D

Việc ứng dụng mô hình 3D trong nghiên cứu dòng chảy sông ngòi mang lại nhiều lợi ích, đặc biệt trong bối cảnh phát triển hạ tầng và quy hoạch đô thị ven sông. Mô hình 3D giúp mô phỏng chính xác hơn các yếu tố thủy động lực, vận chuyển bùn cát, và tác động của các công trình xây dựng đến dòng chảy. Điều này rất quan trọng để đảm bảo an toàn và hiệu quả của các dự án ven sông.

1.2. Giới Thiệu Sông Soài Rạp và Khu Cảng Hiệp Phước

Sông Soài Rạp đóng vai trò quan trọng trong hệ thống giao thông thủy của TP.HCM và khu vực lân cận. Khu Cảng Hiệp Phước, đang được xây dựng trên sông Soài Rạp, hứa hẹn trở thành một trung tâm logistics quan trọng. Việc nghiên cứu dòng chảy và đánh giá tác động môi trường là cần thiết để đảm bảo sự phát triển bền vững của khu vực này.

II. Thách Thức Khi Mô Phỏng Dòng Chảy Sông Soài Rạp 3D

Mô phỏng dòng chảy sông Soài Rạp bằng mô hình 3D đối mặt với nhiều thách thức. Địa hình sông 3D phức tạp và biến đổi liên tục do xói lở, bồi lắng. Chế độ thủy văn chịu ảnh hưởng mạnh mẽ của thủy triều và các yếu tố khí tượng. Ngoài ra, việc thu thập dữ liệu đầu vào chính xác và đầy đủ cũng là một khó khăn lớn. Sai sót trong dữ liệu địa hình, độ nhám kênh, hoặc điều kiện biên có thể dẫn đến sai số lớn trong kết quả mô phỏng. Việc calibration mô hình và validation mô hình với dữ liệu thực tế là rất quan trọng để đảm bảo độ tin cậy của mô hình.

2.1. Thu Thập Dữ Liệu Địa Hình và Thủy Văn Chi Tiết

Để xây dựng một mô hình 3D chính xác, cần có dữ liệu địa hình chi tiết và cập nhật. Dữ liệu này có thể được thu thập bằng các phương pháp khảo sát truyền thống, dữ liệu LiDAR, hoặc các kỹ thuật viễn thám khác. Dữ liệu thủy văn, bao gồm mực nước, vận tốc dòng chảy, và lưu lượng, cũng cần được thu thập liên tục tại nhiều vị trí trên sông.

2.2. Xử Lý và Hiệu Chỉnh Sai Số Mô Hình 3D

Mọi mô hình đều có sai số. Cần áp dụng các phương pháp hiệu chỉnh mô hình để giảm thiểu sai số và tăng độ chính xác. Các phương pháp này có thể bao gồm điều chỉnh các tham số mô hình, sử dụng các kỹ thuật tối ưu hóa mô hình, và so sánh kết quả mô phỏng với dữ liệu thực tế.

2.3. Ảnh hưởng của Biến Đổi Khí Hậu và Nước Biển Dâng

Biến đổi khí hậu và nước biển dâng có thể gây ra những thay đổi đáng kể trong chế độ thủy văn của sông Soài Rạp, bao gồm tăng tần suất và cường độ ngập lụt, thay đổi dòng chảy, và xâm nhập mặn. Các yếu tố này cần được xem xét trong quá trình mô phỏng để đảm bảo tính chính xác và tin cậy của kết quả.

III. Phương Pháp Mô Hình 3D ECOMSED Tính Dòng Chảy Sông

Luận văn này sử dụng mô hình ECOMSED để mô phỏng dòng chảy sông Soài Rạp. ECOMSED là một mô hình số ba chiều mã nguồn mở, được phát triển từ Princeton Ocean Model. Mô hình này sử dụng hệ tọa độ cong trực giao theo phương ngang và hệ tọa độ Sigma theo phương đứng, cùng với phương pháp sai phân hữu hạn để giải các phương trình cơ bản. ECOMSED có khả năng mô phỏng sự phân bố theo thời gian của mực nước, dòng chảy, nhiệt độ, độ mặn, chất chỉ thị, và vận chuyển bùn cát, làm cho nó trở thành một công cụ mạnh mẽ cho nghiên cứu thủy động lực học và môi trường ven biển. Mô hình ECOMSED hoàn chỉnh gồm nhiều module, trong đó module thủy động lực học và module chuyên tải bùn cát là quan trọng nhất.

3.1. Giới Thiệu Chi Tiết Mô Hình ECOMSED và Các Module

ECOMSED (Estuarine, Coastal and Ocean Modeling System with Sediments) là một hệ thống mô hình mạnh mẽ, cho phép mô phỏng các quá trình thủy động lực, vận chuyển bùn cát, và chất lượng nước trong môi trường cửa sông, ven biển và đại dương. Mô hình này được phát triển từ giữa thập niên 1980 và đã được sử dụng rộng rãi trong nhiều nghiên cứu khoa học và ứng dụng thực tế. ECOMSED có cấu trúc module hóa, cho phép người dùng lựa chọn các module phù hợp với mục tiêu nghiên cứu của mình.

3.2. Thiết Lập và Chạy Mô Hình ECOMSED Cho Sông Soài Rạp

Để chạy mô hình ECOMSED cho sông Soài Rạp, cần thực hiện một số bước quan trọng, bao gồm: (1) Thu thập và xử lý dữ liệu đầu vào, bao gồm dữ liệu địa hình, thủy văn, và khí tượng. (2) Xây dựng lưới tính toán phù hợp với hình dạng và kích thước của sông. (3) Thiết lập các tham số mô hình, bao gồm độ nhám đáy, hệ số nhớt rối, và điều kiện biên. (4) Chạy mô hình và kiểm tra tính ổn định và hội tụ của kết quả.

3.3. Calibration và Validation Mô Hình ECOMSED

Calibration và Validation là hai bước quan trọng để đảm bảo độ chính xác và tin cậy của mô hình. Calibration là quá trình điều chỉnh các tham số mô hình để sao cho kết quả mô phỏng phù hợp với dữ liệu thực tế. Validation là quá trình so sánh kết quả mô phỏng với dữ liệu thực tế độc lập để đánh giá khả năng dự đoán của mô hình. Cả hai bước này đều đòi hỏi sự cẩn trọng và kinh nghiệm của người sử dụng.

IV. Ứng Dụng Mô Hình ECOMSED Đánh Giá Dòng Chảy Sông

Mục đích của luận văn là nghiên cứu ứng dụng mô hình 3D ECOMSED vào tính toán dòng chảy đa chiều và chuyên tải bùn cát trong sông Soài Rạp, đoạn từ mũi Bình Khánh đến cửa sông Soài Rạp để đánh giá tình hình biến đổi lòng dẫn trên đoạn sông này. Trong luận văn, mô hình thủy động lực học và bùn cát, sai phân hữu hạn ba chiều ECOMSED đã được sử dụng. Kết quả tính toán đã mô phỏng được chế độ dòng chảy bán nhật triều và đánh giá được tình hình biến đổi lòng dẫn trên sông Soài Rạp. Các kết quả về vận tốc, ứng suất đáy và hệ số nhớt rối được phân tích để hiểu rõ hơn về động lực dòng chảy trong khu vực nghiên cứu.

4.1. Phân Tích Kết Quả Mô Phỏng Thủy Động Lực Học

Kết quả mô phỏng thủy động lực học cho thấy mô hình ECOMSED có khả năng tái hiện khá tốt chế độ thủy triều và dòng chảy trong sông Soài Rạp. Các kết quả về mực nước, vận tốc dòng chảy, và ứng suất đáy phù hợp với các quan sát thực tế. Tuy nhiên, vẫn còn một số khác biệt giữa kết quả mô phỏng và dữ liệu thực tế, đặc biệt là ở các khu vực có địa hình phức tạp.

4.2. Đánh Giá Biến Đổi Lòng Dẫn và Vận Chuyển Bùn Cát

Kết quả mô phỏng vận chuyển bùn cát cho thấy có sự bồi lắng ở một số khu vực và xói lở ở những khu vực khác. Sự phân bố của bùn cát phụ thuộc vào nhiều yếu tố, bao gồm chế độ thủy triều, dòng chảy, địa hình đáy, và nguồn cung cấp bùn cát. Kết quả này có ý nghĩa quan trọng cho việc quy hoạch và quản lý lòng sông.

4.3. So Sánh Kết Quả Mô Phỏng Với Dữ Liệu Quan Trắc

Việc so sánh kết quả mô phỏng với dữ liệu quan trắc là cần thiết để đánh giá độ chính xác và tin cậy của mô hình. Các thông số được so sánh bao gồm mực nước, vận tốc dòng chảy, nồng độ bùn cát lơ lửng, và sự thay đổi cao trình đáy. Nếu có sự khác biệt lớn giữa kết quả mô phỏng và dữ liệu quan trắc, cần xem xét lại quá trình thiết lập và calibration mô hình.

V. Kết Luận và Hướng Phát Triển Nghiên Cứu 3D Dòng Chảy

Nghiên cứu này đã chứng minh khả năng ứng dụng hiệu quả của mô hình 3D ECOMSED trong việc mô phỏng dòng chảy và vận chuyển bùn cát trong sông Soài Rạp. Kết quả nghiên cứu cung cấp thông tin quan trọng cho việc quy hoạch, quản lý, và bảo vệ môi trường ven biển. Trong tương lai, cần tiếp tục nghiên cứu và phát triển các mô hình 3D phức tạp hơn để mô phỏng chính xác hơn các quá trình tự nhiên và tác động của con người đến môi trường sông ngòi. Ngoài ra, cần tăng cường hợp tác giữa các nhà khoa học, kỹ sư, và nhà quản lý để ứng dụng các kết quả nghiên cứu vào thực tiễn.

5.1. Tóm Tắt Những Kết Quả Chính và Đóng Góp Của Nghiên Cứu

Nghiên cứu đã thành công trong việc xây dựng và ứng dụng mô hình ECOMSED để mô phỏng dòng chảy và vận chuyển bùn cát trong sông Soài Rạp. Kết quả nghiên cứu cung cấp thông tin quan trọng cho việc quy hoạch, quản lý, và bảo vệ môi trường ven biển. Nghiên cứu cũng đã xác định được một số thách thức và hướng phát triển trong việc sử dụng mô hình 3D cho nghiên cứu sông ngòi.

5.2. Hướng Phát Triển Nghiên Cứu và Ứng Dụng Mô Hình 3D

Trong tương lai, cần tiếp tục nghiên cứu và phát triển các mô hình 3D phức tạp hơn để mô phỏng chính xác hơn các quá trình tự nhiên và tác động của con người đến môi trường sông ngòi. Cần tích hợp thêm các yếu tố như biến đổi khí hậu, nước biển dâng, và tác động của các công trình xây dựng vào mô hình. Ngoài ra, cần phát triển các công cụ hỗ trợ quyết định để giúp các nhà quản lý sử dụng kết quả mô phỏng để đưa ra các quyết định thông minh và hiệu quả.

5.3. Kiến Nghị Cho Quản Lý và Quy Hoạch Lưu Vực Sông Soài Rạp

Cần có sự phối hợp chặt chẽ giữa các cơ quan quản lý, các nhà khoa học, và cộng đồng địa phương trong việc quản lý và quy hoạch lưu vực sông Soài Rạp. Cần ưu tiên các giải pháp bền vững và thân thiện với môi trường. Cần sử dụng các công cụ mô phỏng và phân tích để đánh giá tác động của các dự án phát triển đến môi trường và cộng đồng.

Nội dung chính

Tổng quan nghiên cứu

Sông Soài Rạp, đoạn từ mũi Bình Khánh đến cửa sông Soài Rạp dài khoảng 40 km, là tuyến luồng quan trọng phục vụ vận tải hàng hóa cho cụm cảng Thành phố Hồ Chí Minh – trung tâm kinh tế lớn nhất Việt Nam với dân số hơn 8 triệu người. Theo ước tính, lưu lượng hàng hóa qua cụm cảng này ngày càng tăng, đòi hỏi tuyến luồng phải đảm bảo an toàn cho tàu trọng tải lớn, vượt quá khả năng của tuyến luồng Lòng Tàu hiện tại. Do đó, việc nghiên cứu dòng chảy và chuyển tải bùn cát trên sông Soài Rạp nhằm đánh giá biến đổi lòng dẫn là rất cần thiết để phục vụ công tác nạo vét, cải tạo luồng tàu.

Mục tiêu chính của luận văn là ứng dụng mô hình 3D ECOMSED để tính toán dòng chảy và chuyển tải bùn cát trên đoạn sông này, từ đó đánh giá tình hình biến đổi lòng dẫn. Nghiên cứu tập trung vào mô phỏng chế độ dòng chảy bán nhật triều và sự phân bố bùn cát trong cột nước, đồng thời dự báo sự thay đổi cao trình đáy sông trong khoảng thời gian mô phỏng. Phạm vi nghiên cứu bao gồm đoạn sông Soài Rạp dài khoảng 40 km, với dữ liệu thu thập từ các trạm đo mực nước, vận tốc dòng chảy và đặc điểm bùn cát tại khu vực.

Ý nghĩa của nghiên cứu thể hiện qua việc cung cấp cơ sở khoa học cho việc quy hoạch, quản lý và khai thác tuyến luồng mới, góp phần nâng cao hiệu quả vận tải thủy, đảm bảo an toàn hàng hải và phát triển bền vững vùng kinh tế động lực phía Nam Việt Nam.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Luận văn sử dụng mô hình thủy động lực học và chuyển tải bùn cát ba chiều ECOMSED, phát triển từ mô hình Princeton Ocean Model, với các đặc điểm chính:

• Mô hình thủy động lực học (ECOM): Giải hệ phương trình động lượng Reynolds và phương trình liên tục trong hệ tọa độ sigma theo phương đứng và tọa độ cong trực giao theo phương ngang. Mô hình tính toán vận tốc dòng chảy ba chiều, mực nước, nhiệt độ, độ mặn, và các thông số rối theo phương đứng và ngang. Kỹ thuật phân tích trường tính toán (Mode Splitting Technique) được áp dụng để tách riêng trường vận tốc trung bình và trường vận tốc dư, giúp tiết kiệm thời gian tính toán.

• Mô hình chuyển tải bùn cát (SED): Mô phỏng chuyển tải bùn cát dính và không dính trong cột nước, bao gồm các quá trình tái lơ lửng, lắng đọng, keo tụ và cố kết bùn cát. Mô hình sử dụng các công thức tính ứng suất cắt đáy, vận tốc ma sát, và vận tốc lắng dựa trên các nghiên cứu thực nghiệm và lý thuyết như công thức của Krone, Partheniades, Van Rijn.

• Khái niệm chính:

  • Ứng suất cắt đáy (T) và vận tốc ma sát (u*) là các tham số quan trọng quyết định sự tái lơ lửng và chuyển tải bùn cát.
  • Hệ tọa độ sigma cho phép mô hình hóa biến đổi theo chiều sâu chính xác hơn.
  • Sự “cứng hóa” đáy bùn cát không dính được mô hình hóa để phản ánh hiện tượng lớp bùn cát thô còn lại sau khi lớp mịn bị xói mòn.

Phương pháp nghiên cứu

• Nguồn dữ liệu:
  Dữ liệu địa hình đáy sông, thủy văn, khí tượng, bùn cát được thu thập từ các trạm đo thực tế tại khu vực sông Soài Rạp và các báo cáo ngành liên quan. Dữ liệu mực nước và vận tốc dòng chảy được sử dụng để hiệu chỉnh và kiểm nghiệm mô hình.

• Phương pháp phân tích:
  Mô hình ECOMSED được thiết lập với lưới tính toán chi tiết, sử dụng phương pháp sai phân hữu hạn bậc hai để giải các phương trình thủy động lực và chuyển tải bùn cát. Các điều kiện biên được xác định dựa trên dữ liệu mực nước triều và vận tốc dòng chảy thực tế. Mô hình được hiệu chỉnh qua các bước: hiệu chỉnh mực nước, vận tốc trung bình và phân tích độ nhạy với các tham số như chiều cao nhám đáy, hệ số xáo trộn rối.

• Timeline nghiên cứu:
  Nghiên cứu được thực hiện trong khoảng thời gian từ tháng 7/2012 đến tháng 11/2013, bao gồm thu thập dữ liệu, thiết lập mô hình, hiệu chỉnh, mô phỏng và phân tích kết quả.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

1. Mô phỏng chế độ dòng chảy bán nhật triều:
   Mô hình ECOMSED đã mô phỏng chính xác mực nước và vận tốc dòng chảy tại các vị trí thượng lưu, hạ lưu và hợp lưu với sông Vàm Cỏ. Sai số giữa kết quả mô phỏng và số liệu thực tế tại trạm Nhà Bè dưới 10%, thể hiện độ tin cậy cao của mô hình.

2. Phân bố vận tốc theo chiều sâu và phương ngang:
   Vận tốc tại tầng mặt và tầng đáy có sự khác biệt rõ rệt, với vận tốc tầng đáy thấp hơn khoảng 25-30% so với tầng mặt. Vận tốc phân bố không đều theo chiều ngang, đặc biệt tại các khúc cong sông, ảnh hưởng đến quá trình chuyển tải bùn cát.

3. Chuyển tải bùn cát lơ lửng:
   Nồng độ bùn cát lơ lửng dao động từ vài chục đến vài trăm mg/l, phân bố theo chiều sâu và thời gian có sự biến đổi rõ rệt. Sau 10 ngày mô phỏng, sự thay đổi cao trình đáy bùn cát tại một số vị trí có thể lên đến vài mm, cho thấy quá trình bồi lắng và xói mòn diễn ra liên tục.

4. Ảnh hưởng của các tham số mô hình:
   Thay đổi chiều cao nhám đáy và hệ số xáo trộn rối ảnh hưởng đáng kể đến kết quả vận tốc và nồng độ bùn cát, với sự thay đổi vận tốc trung bình theo độ sâu lên đến 15%, cho thấy cần hiệu chỉnh kỹ lưỡng các tham số này để đạt độ chính xác cao.

Thảo luận kết quả

Kết quả mô phỏng cho thấy mô hình 3D ECOMSED có khả năng tái hiện chính xác các đặc điểm thủy động lực và chuyển tải bùn cát trong sông Soài Rạp, phù hợp với các nghiên cứu trước đây tại các cửa sông và vùng ven biển khác. Việc mô phỏng chi tiết theo chiều sâu giúp hiểu rõ hơn về sự phân bố vận tốc và nồng độ bùn cát, điều mà các mô hình 1D hoặc 2D không thể cung cấp đầy đủ.

Biểu đồ so sánh mực nước mô phỏng và thực đo tại trạm Nhà Bè, cùng bảng số liệu vận tốc trung bình theo độ sâu tại các vị trí khác nhau, minh họa rõ sự phù hợp giữa mô hình và thực tế. Sự biến đổi cao trình đáy bùn cát được thể hiện qua biểu đồ thay đổi cao trình sau 10 và 20 ngày mô phỏng, cho thấy xu hướng bồi lắng và xói mòn cục bộ.

Nguyên nhân của sự biến đổi lòng dẫn được giải thích bởi tác động của dòng chảy triều và vận tốc phân bố không đồng đều, kết hợp với đặc tính bùn cát dính và không dính trong cột nước. So với các nghiên cứu trong nước và quốc tế, kết quả này khẳng định ưu điểm của mô hình 3D trong việc mô phỏng các quá trình phức tạp tại cửa sông.

Đề xuất và khuyến nghị

1. Tăng cường thu thập dữ liệu thực tế:
   Đề xuất thực hiện các đợt đo đạc bổ sung về vận tốc dòng chảy, nồng độ bùn cát và đặc tính đáy tại các vị trí trọng yếu trên sông Soài Rạp để hiệu chỉnh mô hình chính xác hơn. Thời gian thực hiện: 6-12 tháng, chủ thể: Sở Tài nguyên và Môi trường, Viện Khoa học Thủy lợi.

2. Ứng dụng mô hình ECOMSED trong quy hoạch nạo vét:
   Sử dụng kết quả mô phỏng để xác định các khu vực cần nạo vét, dự báo biến đổi lòng dẫn trong dài hạn nhằm tối ưu hóa chi phí và hiệu quả công tác duy tu luồng tàu. Thời gian áp dụng: 1-3 năm, chủ thể: Ban Quản lý Cảng, Sở Giao thông Vận tải.

3. Phát triển hệ thống cảnh báo biến đổi lòng dẫn:
   Xây dựng hệ thống cảnh báo sớm dựa trên mô hình số để theo dõi và dự báo biến đổi lòng dẫn, giúp kịp thời điều chỉnh hoạt động vận tải và bảo đảm an toàn hàng hải. Thời gian triển khai: 2 năm, chủ thể: Trung tâm Quản lý Đường thủy nội địa, Cục Hàng hải Việt Nam.

4. Nâng cao năng lực nghiên cứu và đào tạo:
   Tổ chức các khóa đào tạo chuyên sâu về mô hình thủy động lực và chuyển tải bùn cát cho cán bộ kỹ thuật, nghiên cứu viên nhằm nâng cao năng lực ứng dụng mô hình số trong quản lý tài nguyên nước. Thời gian: liên tục, chủ thể: Trường Đại học Bách Khoa TP.HCM, Viện Khoa học Thủy lợi.

Đối tượng nên tham khảo luận văn

1. Các nhà quản lý cảng và giao thông thủy:
   Giúp hiểu rõ về biến đổi lòng dẫn và tác động đến an toàn vận tải, từ đó đưa ra các quyết định quản lý hiệu quả.

2. Chuyên gia và kỹ sư thủy lợi, thủy văn:
   Cung cấp phương pháp và công cụ mô hình 3D tiên tiến để nghiên cứu dòng chảy và chuyển tải bùn cát trong sông và vùng ven biển.

3. Nhà nghiên cứu và sinh viên chuyên ngành xây dựng công trình thủy:
   Là tài liệu tham khảo quý giá về ứng dụng mô hình số trong tính toán thủy động lực và bùn cát, đồng thời cập nhật các kiến thức lý thuyết và thực tiễn.

4. Cơ quan quản lý tài nguyên nước và môi trường:
   Hỗ trợ đánh giá tác động của biến đổi lòng dẫn đến môi trường và phát triển bền vững vùng cửa sông, ven biển.

Câu hỏi thường gặp

1. Mô hình ECOMSED có ưu điểm gì so với các mô hình khác?
   ECOMSED là mô hình 3D, giải các phương trình thủy động lực và chuyển tải bùn cát với độ chính xác cao, đặc biệt trong việc mô phỏng biến đổi theo chiều sâu và các quá trình phức tạp như keo tụ, cố kết bùn cát dính. Ví dụ, mô hình này đã được ứng dụng thành công tại nhiều cửa sông và vùng ven biển trên thế giới.

2. Dữ liệu đầu vào quan trọng nhất để mô hình hoạt động chính xác là gì?
   Dữ liệu địa hình đáy sông, mực nước triều, vận tốc dòng chảy và đặc tính bùn cát là các yếu tố quyết định. Việc thu thập và hiệu chỉnh dữ liệu thực tế giúp mô hình phản ánh đúng điều kiện tự nhiên.

3. Thời gian mô phỏng và phạm vi không gian của mô hình có giới hạn không?
   Mô hình có thể mô phỏng trong khoảng thời gian từ vài ngày đến vài tháng tùy theo mục tiêu nghiên cứu, với phạm vi không gian linh hoạt từ vài km đến hàng chục km, phù hợp với quy mô sông Soài Rạp.

4. Làm thế nào để hiệu chỉnh mô hình cho phù hợp với điều kiện thực tế?
   Hiệu chỉnh được thực hiện bằng cách so sánh kết quả mô phỏng với số liệu thực đo về mực nước và vận tốc, điều chỉnh các tham số như chiều cao nhám đáy, hệ số xáo trộn rối để giảm sai số.

5. Ứng dụng thực tiễn của nghiên cứu này trong quản lý cảng biển là gì?
   Nghiên cứu giúp xác định các khu vực cần nạo vét, dự báo biến đổi lòng dẫn, từ đó đảm bảo an toàn cho tàu trọng tải lớn, tối ưu hóa hoạt động vận tải và phát triển bền vững cảng biển.

Kết luận

• Luận văn đã thành công trong việc ứng dụng mô hình 3D ECOMSED để mô phỏng dòng chảy và chuyển tải bùn cát trên đoạn sông Soài Rạp dài khoảng 40 km.
• Kết quả mô phỏng phù hợp với số liệu thực tế, thể hiện qua sai số dưới 10% về mực nước và vận tốc dòng chảy.
• Phân bố vận tốc và nồng độ bùn cát theo chiều sâu và phương ngang được mô phỏng chi tiết, giúp đánh giá chính xác biến đổi lòng dẫn.
• Nghiên cứu cung cấp cơ sở khoa học cho công tác nạo vét, cải tạo luồng tàu và quản lý tài nguyên nước tại khu vực.
• Đề xuất các giải pháp nâng cao hiệu quả quản lý và phát triển bền vững vùng cửa sông, đồng thời khuyến nghị tiếp tục thu thập dữ liệu và đào tạo chuyên môn.

Next steps: Triển khai các đề xuất thu thập dữ liệu bổ sung, áp dụng mô hình trong quy hoạch nạo vét và xây dựng hệ thống cảnh báo biến đổi lòng dẫn.

Các cơ quan quản lý và nhà nghiên cứu nên phối hợp để ứng dụng mô hình ECOMSED rộng rãi, nâng cao hiệu quả quản lý tài nguyên nước và phát triển kinh tế vùng cửa sông.

Trích đoạn nội dung tài liệu

Chương 1: Mở dau (ii) Chương 2: Cơ sở lý thuyết của mô hình ECOMSED (ii) Chương 3: Điều kiện tự nhiên khu vực nghiên cứu (iv) Chương 4: Ứng dụng mô hình Ecomsed dé tính toán dòng chảy đoạn sông Soài Rạp (v) Chương 5: Ứng dụng mô hình Ecomsed dé tính toán chuyền tải bùn cát đoạn sông Soài Rạp (vi) Chương 6: Kết luận và kiến nghị CHƯƠNG 2: CƠ SỞ LÝ THUYET CUA MÔ HÌNH ECOMSED 2.1 Giới thiệu Mô hình ECOMSED đã được phát triển từ giữa thập niên 1980 với mô hình tiền thân là Princeton Ocean Model (Blumberg va Mellor, 1987) và phiên bản ECOM dùng cho môi trường nước nông trong sông, vịnh, cửa sông, vùng ven biển, hỗ chứa tự nhiên và nhân tao (Blumberg, 1996). Từ giữa thập niên 1990, bài toán chuyển tải bùn cát dính và không dính, bồi lăng và cô kết,. được kết hợp vào ECOM để cho ra mô hình ECOMSED (Estuarine, Coastal and Ocean Modeling System with Sediments). Đây là một chương trình mã nguồn mở viết băng ngôn ngữ Fortran của tổ chức HydroQual — Hoa Ky.

ECOMSED sử dụng hệ tọa độ cong trực giao theo phương ngang và hệ tọa độ Sigma theo phương đứng, dùng phương pháp sai phân hữu hạn để giải các phương trình cơ bản. Mô hình được xây dựng để mô phỏng sự phân bố theo thời gian của mực nước, dòng chảy, nhiệt độ. độ mặn, chat chỉ thị, bùn cát dính - không dính, sóng biển va chất lượng nước. Mô hình ECOMSED hoàn chỉnh gồm nhiều module.

Đó là module thủy động lực hoc, module chuyên tải bùn cát, module sóng do gió, module dòng nhiệt, module đường đi của hạt và module chất lượng nước. Các chương trình con của mô hình làm việc liên kết với nhau, dữ liệu đầu ra của chương trình này là dữ liệu đầu vào của chương trình kia.1 trình bay sơ đồ khối của mô hình ECOMSED. Hiệu quả của mô hình đã được đánh giá rất cao qua một SỐ lượng lớn các thử nghiệm và nhiều ứng dụng trong thực tế. Đã có hơn 350 bài báo được viết từ những ứng dụng của ECOMSED.

ECOMSED đã được chứng minh qua nhiều năm là mạnh mẽ và đáng tin cậy. Trong khuôn khổ luận văn này, chỉ hai module thủy động lực và chuyển tải bùn cát được trình bày.1: Khung mô hình ECOMSED e GIÓ e Mây bao phủ e So dé day bùn cát _e Khối lượng bùn cát e Nhiệt độ không khí e Lượng mua e Đặc tính đáy e Điều kiện ban đầu e Bức xạ sóngngăn e Bốc hoi e Thế năng xói eLuu lượng nước e Nhiệt độ, độ mặn e Thông số lắng e DO âm tương đôi s Ap suât khí quyên sạch đên ban đâu Ỷ e Lưu lượng đến CSO e Đặc điểm khí Ỷ A e Điều kiện biên: nhiệt tượng Vv Module Chuyển tải bùn cát |g Module Thông lượng nhiệt độ. độ mặn, mực nước SED Ỷ Ỷ Ỷ Vv Module Thủy dOngluchoc |, e Buc xa mặt trời e Nông độ bùn cát lơ lửng ECOM e Bức xạ khí quyén e Khối lượng lắng/xói e Thông lượng nhiệt cảm nhận Ỷ e Sự thay đổi cao trình đáy e Thông lượng nhiệt bay hơi e Mực nước e Dòng chảy e Nhiệt độ, độ mặn e Xáo trộn rôi 1 | | e Toa độ lưới e Gid | e Dong chay e Mực nước | Vv e Miền vận chuyển se Khối lượng | e Khuéch tán e Dong chay | e Nhiệt độ, độ man Carbon | Ỷ e Khôi lượng chất e Các thông số chất Ỷ | | mùn lượng nước khác | Module Theo dõi hạt A Module Song ! | Ỷ Ỷ | | Module Chất lượng nước | e Duong di của hạt A e Chiều cao sóng điển hình | RCA | e Bức xạ khí quyền e Chu kỳ sóng | e Vector chuyên động Ỷ e Huong sóng | | e DO e Kali trong phân e Ma sat day | | e Chat min e Kim loại e Độ kiém, pH e Thực vật phù du Dau vao Module Đầu ra e TIC, CaCOz 2.2 Module thủy dong lực học Module thủy động lực học (ECOM) là một mô hình dòng chảy ba chiều, kết hợp với một mô hình rối khép kín để cung cấp các thông số thực tế của quá trình xáo trộn theo phương đứng. Các biến tính toán là ba thành phần của vận tốc, nhiệt độ, độ mặn, động năng rỗi, và rỗi kích thước lớn (turbulence macroscale).

Các phương trình động lượng là phi tuyến có kết hợp với thông số Coriolis. Các phương trình tính toán dự báo nhiệt động lực học, nhiệt độ và độ mặn có kế đến sự thay đối của khối nước do hiện tượng nước trồi/nước sụt doc bờ biển theo thời gian cũng như quá trình đối lưu theo phương ngang của dòng chảy. Cao trình mặt thoáng, thủy triều và sóng bão cũng có thé được mô phỏng dự báo. Đó là nhờ mồ hình sử dụng kỹ thuật phân tích trường tính toán (Mode Splitting Technique), trong đó sự chuyền tải thé tích và vận tốc cắt theo phương đứng được giải riêng biệt.

Các biến tính toán khác bao gồm khối lượng riêng. độ nhớt rồi và khuếch tán rối theo phương đứng. Module cũng xét đến địa hình đáy và hình dạng thực của bờ biển. Module này có rất nhiều ứng dụng thành công cho đại dương, bờ biển và cửa sông.

Một vài ứng dụng có thể kế đến như Chesapeake Bay (Blumberg và Goodrich, 1990), khu vực Gulf Stream (Ezer va Mellor, 1992), Georges Bank (Chen và cộng sự, 1995), the Oregon Continental Shelf (Allen và cộng sự, 1995), New York Bight va New York Harbor (Blumberg va cộng sự, [999).1 Hệ phương trình thủy động lực hoc Các phương trình cơ bản gồm vận tốc và cao độ mặt thoáng. Hai phép xấp xi được sử dụng là giả thiết thủy tĩnh (trọng lượng chất lỏng cân băng với áp suất) và xấp xi Boussinesq (khối lượng riêng không đổi). Xem xét một hệ tọa độ Descartes vuông góc với chiều dương trục x theo hướng đông. chiều dương trục y theo hướng bắc, và chiều dương trục z hướng lên.

Mặt thoáng được đặt ở z = n(x,y,f) va mặt đáy đặt ở z = -H(x,y). V là véctơ vận tốc theo phương ngang có hai thành phan (U,V) và V là toán tử gradient theo phương ngang. Phương trình liên tục: = aw V.1) Phương trình động lượng Reynolds: a+ V.VU + Wwe fy =i 4 (Ky) +, gŨ = gŨ 1 OP 0 gŨ (2.4) trong đó: p, là khối lượng riêng quy chuẩn, p là khối lượng riêng của chất lưu, g là gia tốc trọng trường, P là áp suất, Ky, là hệ số khuếch tán rối theo phương đứng, f là thông SỐ Coriolis. Áp suất tại độ sâu z có thé được xác định băng cách tích phân các thành phan theo phương đứng của phương trình chuyển động, từ z đến mặt thoáng n: 0 P(x,y,Z,E) = Patm + 9Pon+g J p(x,y,z,£) dz' (2.5) Những chuyên động xoáy rôi có kích thước nhỏ hơn kích thước lưới sẽ tích tụ năng lượng rôi và sinh ra su mat cân băng năng lượng trong dòng chảy, dân dén mat ôn định của lời giải sô.

Điều này có thê được loại bỏ nhờ các sô hạng Fx, Fy trong các phương trình (2. Các số hạng này được viết như sau: 8 Ị 8U 8 Ị OU OV \] = 2 |24y | +2 [Âm (= + ¬) (2.6a) 8 Ị ov | 8 Ị OU OV \] B= 2 [2au 2] +2 [An (2 + ~) (2.6b) 10 trong đó F¿ và Fy là không đổi với sự xoay hệ trục tọa độ. Ay và Ay là các hệ SỐ khuếch tán rối theo phương ngang.2 Mô hình rồi Các phương trình chủ đạo chứa các ứng suất Reynolds và các thông lượng kế đến sự khuếch tán rối của động lượng, nhiệt và độ mặn. Hệ số xáo trộn theo phương đứng, Ky và Ky trong phương trình (2.3) tim được nhờ việc sử dụng sơ đồ khép kín rối bậc hai bao gồm một phương trình tính động năng rối q /2 và một phương trình tinh kích thước rồi quy mô lớn Ê như sau: #U +V.v42 + W SEz Ôz=*(K,){ az + 2Ku (2) 0z +() 0z |+ Po2K, 0z2-44 Byt eq (2.7) va ar TỰ.V(q7l) V2 +W az =—|K,(420| + tEiKu (2) + (“) | + Ku 2= md +, (28) trong đó V là toán tử gradient theo phương ngang.10) Các hệ số khuếch tán rỗi Ky, Ky và K, được xác định như sau: Ky = tqSy (2.11b) II Những hàm số ôn định (stability functions) Sy, Sy được xác định bang giai tich, theo các quan hệ đại số độc lập với ôU/ôz.

OV/6z, gp;1, ôp/ôz. q va C: BỊ "/3~3AAa6n|(B.15) là tần suất Brunt-Vaisala. Các hệ số kinh nghiệm theo Mellor và Yamada (1982) lần lượt là: Al; A2; B1; B2 = 0,92 : 0,74: 16,6 ; 10,1 2.3 Điều kiện biên tại mặt thoáng và day Điều kiện biên tại mặt thoáng z = n(x.16a) OU 03V q? = Bị Šuậ, (2.16d) trong đó (Tox, Toy) là vecto ứng suat gió bê mặt với van toc ma sat u,, là độ lớn cua vecto. Gia tri BY Bia một hằng số kinh nghiệm xuất phát từ các quan hệ rối.

12 Tại biên đáy có các điêu kiện: poKu (=.17€) gH gH W,y = ~Uy2—— Vn (2. Ứng suất đáy được xác định nhờ kết hợp vận tốc theo quy luật logarith của thành ran: Ty = p,Cp]|V,|Ÿ, (2.18) với hệ số lực cản Cp được xác định như sau: 1 In(H+zp) —2 Cy = | K Zo (2.19a) trong đó zp và Vỹ là diém lưới và van toc tương ứng với diém lưới gân nhat với đáy, « là hang số von Karman.19a) với thông số rối Ky cho biểu thức: V = ( tb )m ca (2.19b) ÐKtUtrp Zo trong lớp biên dưới thấp nếu lưới có đủ độ phân giải. Trong trường hợp lớp biên đáy không có độ phân giải tốt thì nên chon Cp = 0. Thuật toán thực tế là gan Cp giá tri lớn hơn trong hai giá trị ở biểu thức (2.

Thông số z, phụ thuộc vào độ nhám day tại chỗ, nếu không có số liệu cụ thé thì lây Z) = 1 cm (Weatherly va Martin 1978).4 Điều kiện biên hở (Điều kiện biên mực nước) Khi phát triển module dòng chảy biển, cần cung cấp những điều kiện biên hở cho phép năng lượng sóng dài (như thủy triều) truyền vào biên hở cũng như năng lượng sóng dài phát ra từ bên trong miễn tính toán. Có một số loại điều kiện biên có thể được dùng cho mục đích này.1 Điều kiện biên chặt Trong loại điêu kiện biên này, mô hình sử dụng mực nước dọc theo các lưới biên mà người sử dung gan từ các dữ liệu quan trac hoặc triêu điêu hòa. Tuy nhiên điêu kiện biên loại nay khá cứng nhac và không cho phép năng lượng sóng dài (triệu hoặc sóng do bão) vào hay ra khỏi khu vực nghiên cứu.