Luận văn thạc sĩ hus ứng dụng mô hình vnu mdec tính toán chế độ thủy động lực và vận chuyển trầm tích vùng cửa sông ven biển hải phòng

Luận văn thạc sĩ HUS nghiên cứu ứng dụng mô hình VNU MDEC trong tính toán chế độ thủy động lực và vận chuyển trầm tích tại Hải Phòng.

• Tác giả: Phạm Văn Tiến
• Người hướng dẫn: GS. Đinh Văn Ưu
• Trường học: Trường Đại học Khoa học Tự nhiên
• Chuyên ngành: Hải dương học
• Đề tài: Ứng dụng mô hình (VNU/MDEC) tính toán chế độ thủy động lực và vận chuyển trầm tích khu vực cửa sông ven biển Hải Phòng
• Loại tài liệu: Luận văn thạc sỹ
• Năm xuất bản: 2012
• Địa điểm: Hà Nội

Mục lục chi tiết

• 1. CHƯƠNG 1: TỔNG QUAN VỀ MÔ HÌNH THỦY ĐỘNG LỰC VÀ VẬN CHUYỂN TRẦM TÍCH
• 1.1. Các nghiên của nước ngoài
• 2. CHƯƠNG 2: MÔ HÌNH VNU/MDEC
• 3. CHƯƠNG 3: KẾT QUẢ NGHIÊN CỨU
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Tổng quan về mô hình VNU MDEC trong tính toán thủy động lực

Mô hình VNU/MDEC là một công cụ mạnh mẽ trong việc tính toán chế độ thủy động lực và vận chuyển trầm tích tại Hải Phòng. Mô hình này được phát triển nhằm đáp ứng nhu cầu nghiên cứu và quản lý tài nguyên nước, đặc biệt là trong các khu vực ven biển. Việc áp dụng mô hình này giúp hiểu rõ hơn về các quá trình động lực học diễn ra trong môi trường nước, từ đó đưa ra các giải pháp quản lý hiệu quả hơn.

1.1. Đặc điểm của mô hình VNU MDEC trong nghiên cứu thủy động lực

Mô hình VNU/MDEC được thiết kế để mô phỏng các quá trình thủy động lực phức tạp, bao gồm cả dòng chảy và vận chuyển trầm tích. Mô hình này sử dụng các phương trình động lực học cơ bản để mô tả sự tương tác giữa nước và trầm tích, từ đó cung cấp cái nhìn sâu sắc về các yếu tố ảnh hưởng đến chế độ thủy động lực.

1.2. Lợi ích của việc sử dụng mô hình VNU MDEC tại Hải Phòng

Việc áp dụng mô hình VNU/MDEC tại Hải Phòng mang lại nhiều lợi ích, bao gồm khả năng dự đoán chính xác các hiện tượng thủy động lực, hỗ trợ trong việc quy hoạch và quản lý tài nguyên nước, cũng như giảm thiểu tác động tiêu cực đến môi trường.

II. Vấn đề và thách thức trong tính toán thủy động lực tại Hải Phòng

Hải Phòng là một trong những khu vực có chế độ thủy động lực phức tạp do sự tương tác giữa sông và biển. Các vấn đề như xói mòn, lắng đọng trầm tích và biến đổi môi trường đang đặt ra nhiều thách thức cho công tác quản lý. Việc hiểu rõ các vấn đề này là rất cần thiết để có thể áp dụng các giải pháp hiệu quả.

2.1. Các vấn đề chính trong vận chuyển trầm tích tại Hải Phòng

Vận chuyển trầm tích tại Hải Phòng gặp phải nhiều vấn đề như sự thay đổi dòng chảy do triều, ảnh hưởng của gió và các hoạt động con người. Những yếu tố này làm cho việc dự đoán và quản lý trầm tích trở nên khó khăn hơn.

2.2. Thách thức trong việc áp dụng mô hình VNU MDEC

Mặc dù mô hình VNU/MDEC có nhiều ưu điểm, nhưng việc áp dụng nó cũng gặp phải một số thách thức như yêu cầu về dữ liệu đầu vào chính xác, khả năng tính toán và hiệu chỉnh mô hình để phù hợp với điều kiện thực tế tại Hải Phòng.

III. Phương pháp chính trong mô hình VNU MDEC cho tính toán thủy động lực

Mô hình VNU/MDEC sử dụng nhiều phương pháp khác nhau để tính toán chế độ thủy động lực và vận chuyển trầm tích. Các phương pháp này bao gồm việc sử dụng các phương trình động lực học, điều kiện biên và các kỹ thuật số để mô phỏng các quá trình tự nhiên.

3.1. Phương pháp biến đổi tọa độ cong σ trong mô hình

Phương pháp biến đổi tọa độ cong σ là một trong những kỹ thuật quan trọng trong mô hình VNU/MDEC. Phương pháp này giúp mô phỏng chính xác các dòng chảy trong các khu vực có địa hình phức tạp, từ đó cải thiện độ chính xác của kết quả tính toán.

3.2. Điều kiện biên trong mô hình VNU MDEC

Điều kiện biên là yếu tố quan trọng trong mô hình VNU/MDEC. Việc xác định đúng các điều kiện biên giúp mô hình phản ánh chính xác các quá trình thủy động lực và vận chuyển trầm tích, từ đó cung cấp thông tin hữu ích cho công tác quản lý.

IV. Kết quả nghiên cứu từ mô hình VNU MDEC tại Hải Phòng

Kết quả nghiên cứu từ mô hình VNU/MDEC cho thấy rõ ràng các đặc điểm của chế độ thủy động lực và vận chuyển trầm tích tại Hải Phòng. Những thông tin này không chỉ giúp hiểu rõ hơn về các quá trình tự nhiên mà còn hỗ trợ trong việc đưa ra các quyết định quản lý hiệu quả.

4.1. Trường dòng chảy và mực nước triều tại Hải Phòng

Mô hình VNU/MDEC đã cung cấp các thông tin chi tiết về trường dòng chảy và mực nước triều tại Hải Phòng. Những thông tin này rất quan trọng cho việc quy hoạch và quản lý các hoạt động hàng hải trong khu vực.

4.2. Kết quả tính toán vận chuyển trầm tích lơ lửng

Kết quả tính toán cho thấy sự vận chuyển trầm tích lơ lửng dưới tác động của thủy triều và dòng chảy tổng hợp. Những thông tin này giúp đánh giá tác động của các yếu tố tự nhiên đến quá trình vận chuyển trầm tích trong khu vực.

V. Kết luận và tương lai của nghiên cứu mô hình VNU MDEC

Nghiên cứu mô hình VNU/MDEC đã mở ra nhiều cơ hội mới trong việc hiểu và quản lý chế độ thủy động lực và vận chuyển trầm tích tại Hải Phòng. Tương lai của nghiên cứu này hứa hẹn sẽ mang lại nhiều giá trị cho công tác quản lý tài nguyên nước và bảo vệ môi trường.

5.1. Tầm quan trọng của mô hình VNU MDEC trong nghiên cứu thủy động lực

Mô hình VNU/MDEC đóng vai trò quan trọng trong việc nghiên cứu và quản lý tài nguyên nước tại Hải Phòng. Việc áp dụng mô hình này giúp nâng cao hiệu quả trong công tác quản lý và bảo vệ môi trường.

5.2. Hướng phát triển tương lai cho mô hình VNU MDEC

Trong tương lai, mô hình VNU/MDEC có thể được cải tiến và phát triển thêm để đáp ứng tốt hơn các yêu cầu nghiên cứu và quản lý. Việc tích hợp công nghệ mới và dữ liệu lớn sẽ giúp nâng cao độ chính xác và hiệu quả của mô hình.

Nội dung chính

Tổng quan nghiên cứu

Việt Nam sở hữu hơn 3.260 km bờ biển cùng vùng đặc quyền kinh tế rộng trên 1.000.000 km2, đóng vai trò huyết mạch trong chiến lược phát triển kinh tế biển. Tại khu vực duyên hải Đông Bắc Bộ, thành phố Hải Phòng là đầu mối hàng hải trọng điểm với cụm cảng Hải Phòng và cảng Đình Vũ. Tuy nhiên, vùng cửa sông ven biển nơi đây chịu sự chi phối phức tạp của 5 cửa sông lớn gồm Bạch Đằng, Cấm, Lạch Tray, Văn Úc và Thái Bình, kết hợp cùng chế độ thủy động lực biển vịnh Bắc Bộ. Hiện tượng bồi lắng luồng tàu và xói lở bờ biển đang đặt ra bài toán cấp thiết cho công tác quản trị cảng biển, khi trung bình 20 km đường bờ biển lại xuất hiện một cửa sông có biến động địa hình đáy dữ dội.

Luận văn thạc sĩ chuyên ngành Hải dương học của tác giả Phạm Văn Tiến, dưới sự hướng dẫn khoa học của Giáo sư Đinh Văn Ưu tại Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội, tập trung giải quyết bài toán thủy thạch động lực học tại vùng biển Hải Phòng. Mục tiêu cốt lõi của nghiên cứu là ứng dụng mô hình số trị ba chiều VNU/MDEC để mô phỏng chi tiết trường thủy động lực, hoàn lưu dòng chảy và quy luật lan truyền, vận chuyển trầm tích lơ lửng. Phạm vi không gian nghiên cứu bao phủ từ vùng biển ven bờ được che chắn bởi đảo Cát Bà, Cát Hải đến bán đảo Đồ Sơn và Đình Vũ. Nghiên cứu thực hiện mô phỏng chuỗi 15 ngày liên tục theo thời gian Mặt Trời trung bình, cung cấp cơ sở khoa học để kiểm soát sa bồi luồng hàng hải, định hướng quy hoạch đê chắn cát và tối ưu hóa chi phí nạo vét luồng hàng năm.

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

Khung lý thuyết áp dụng

Nghiên cứu ứng dụng khung lý thuyết thủy thạch động lực học biển hiện đại dựa trên nền tảng hệ phương trình động lực biển nguyên thủy ba chiều kết hợp với hệ phương trình bình lưu - khuếch tán vật chất. Mô hình giải quyết các tương tác phi tuyến giữa ba yếu tố động lực chủ đạo gồm thủy triều, lưu lượng xả từ sông và trường gió mùa. Về lý thuyết vận chuyển bùn cát, luận văn kế thừa các nguyên lý kinh điển về ứng suất trượt đáy tới hạn của Shields, công thức dòng vận chuyển di đáy của Van Rijn và mô hình trao đổi bùn cát ranh giới đáy của Meyer.

Các khái niệm chính được làm rõ bao gồm:

• Hoàn lưu tà áp và chính áp trong đầm phá cửa sông chịu tác động phân tầng mật độ;
• Động năng rối và quãng đường xáo trộn rối theo lý thuyết Kolmogorov kết hợp số Richardson thông lượng;
• Cơ chế bứt xói và lắng đọng của hạt trầm tích lơ lửng dựa trên kích thước hạt d50 dao động từ 0,001 mm đến 1,0 m;
• Tọa độ cong sigma kép giúp chia cột nước thành 5 tầng thẳng đứng với ranh giới phân tầng HLIM bằng 1,86 m, loại bỏ sai số tính toán tại các khu vực sườn dốc ngầm.

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

Nghiên cứu sử dụng phương pháp mô hình hóa số trị 3D thông qua bộ mô hình VNU/MDEC, vốn được phát triển từ mô hình GHER của Đại học Liege. Nguồn dữ liệu đầu vào bao gồm bình đồ độ sâu số hóa từ Bộ Tư lệnh Hải quân, mô hình số hóa độ cao địa hình DEM 90 m và chuỗi số liệu khí tượng - hải văn quan trắc liên tục 42 năm từ 1960 đến 2002 tại trạm hải văn Hòn Dáu. Số liệu thủy văn lưu lượng sông và nồng độ phù sa được thu thập từ mạng lưới quan trắc của Viện Tài nguyên và Môi trường biển Hải Phòng kết hợp cùng các trạm thủy văn thượng lưu.

Lý do lựa chọn phương pháp thể tích hữu hạn trên lưới so le Arakawa C là khả năng bảo toàn khối lượng và xung lượng tối ưu khi giải bài toán cửa sông có biên địa hình phức tạp. Kỹ thuật tách mode cho phép giải tách biệt phương trình 2D chính áp theo bước thời gian nhỏ và phương trình 3D tà áp theo bước thời gian lớn hơn, đảm bảo tính ổn định số trị và tối ưu tài nguyên tính toán. Luận văn thiết lập 19 phương án mô phỏng thực nghiệm số trị từ HP01 đến HP19 với thời gian chạy mô hình chuẩn hóa 15 ngày, bao quát đầy đủ các kịch bản tương tác mùa lũ, mùa kiệt, chế độ nhật triều đều và các hướng gió thịnh hành.

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

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

Mô phỏng từ 19 kịch bản số trị đã chỉ ra những đặc tính cốt lõi của chế độ thủy thạch động lực học vùng cửa sông Hải Phòng:

• Thứ nhất, trường dòng chảy tổng hợp chịu sự chi phối tuyệt đối của dòng triều với chế độ nhật triều đều thuần nhất. Vận tốc dòng chảy tổng hợp đạt giá trị trung bình từ 0,4 m/s đến 1,0 m/s. Vào mùa đông, dòng chảy tầng mặt duy trì ổn định theo hướng tây nam với vận tốc từ 20 cm/s đến 25 cm/s; trong khi mùa hè chuyển dịch sang hướng đông bắc với vận tốc từ 15 cm/s đến 20 cm/s.
• Thứ hai, lưu lượng dòng chảy sông mùa lũ chiếm tới 75% đến 85% tổng lượng nước cả năm, tập trung cao điểm vào các tháng 7, 8 và 9 với tỷ trọng 50% đến 70%. Dòng xả lũ cực đại tạo nên gradient độ muối lớn, hình thành lưỡi nước ngọt mở rộng về phía đảo Trà Bản và làm gia tăng mực nước cửa sông từ 0,001 m đến 0,003 m tại cửa Nam Triệu.
• Thứ ba, phân bố nồng độ trầm tích lơ lửng có sự tương phản sâu sắc theo mùa. Mùa lũ ghi nhận nồng độ phù sa từ 53 mg/l đến 215 mg/l, đạt cực đại từ 700 mg/l đến 964 mg/l tại luồng Cửa Cấm. Ngược lại, mùa kiệt có nồng độ nền thấp từ 42 mg/l đến 94 mg/l, nhưng tại các vùng cạn ven bờ ngoài cửa sông, tác động khuấy đục đáy của sóng và dòng triều làm nồng độ trầm tích cục bộ tăng vọt lên 252 mg/l đến 860 mg/l.
• Thứ tư, lượng bùn cát mịn vận chuyển ra biển trong mùa lũ chiếm hơn 90% tổng lượng bồi tích năm, trong đó tháng 8 chiếm tới 35% đến 40%. Sông Cấm và Lạch Tray là hai nguồn phát tán trầm tích mạnh nhất, gây bồi lắng nghiêm trọng tại khu vực luồng tàu cửa Nam Triệu và Lạch Huyện.

Thảo luận kết quả

Cơ chế bồi tụ và xói lở tại cửa sông ven biển Hải Phòng bắt nguồn từ sự giao thoa giữa sóng triều có biên độ cao lên tới 2,0 m và quá trình phân tầng mật độ nước sông - biển. Khi triều dâng trong khoảng thời gian từ 11 đến 12 giờ, dòng nước mặn ép sâu vào lòng sông làm chậm vận tốc dòng chảy đáy, kích thích quá trình keo tụ bùn sét có kích thước hạt từ 0,001 mm đến 0,007 mm. Khi triều rút kéo dài từ 13 đến 14 giờ, kết hợp cùng lưu lượng sông mùa lũ, vận tốc dòng chảy tầng mặt tăng mạnh đẩy nồng độ phù sa lơ lửng phát tán xa ra vịnh ngoài.

Dữ liệu mô phỏng khi được trực quan hóa qua biểu đồ không gian và các bảng mặt cắt dọc cho thấy sự hình thành rõ rệt của nêm mặn và lớp biến đổi nồng độ trầm tích sát đáy. Các trường dòng chảy và đường đẳng trị nồng độ trầm tích lơ lửng được thể hiện chi tiết tại các mốc thời gian 35 giờ, 102 giờ và 141 giờ, phản ánh chính xác các pha triều lớn và triều ròng. So với các nghiên cứu trước đây tại khu vực vịnh Bắc Bộ sử dụng mô hình 2D tích phân độ sâu, mô hình VNU/MDEC 3D đã tái hiện chuẩn xác hiện tượng dòng chảy ngược tầng đáy, giải thích thỏa đáng vì sao khu vực cửa Nam Triệu luôn bị bồi lấp nhanh hơn dự kiến ngay cả trong mùa kiệt.

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

Dựa trên các kết quả tính toán thủy thạch động lực học, nghiên cứu đưa ra 4 nhóm giải pháp công trình và phi công trình trọng tâm:

1. Điều chỉnh lịch trình nạo vét duy tu luồng hàng hải Hải Phòng: Cục Hàng hải Việt Nam và Tổng công ty Bảo đảm an toàn hàng hải miền Bắc cần tổ chức nạo vét tập trung vào cuối mùa kiệt từ tháng 3 đến tháng 4 hàng năm. Giải pháp này giúp giảm thiểu từ 15% đến 20% khối lượng bồi lắng lặp lại do lũ mùa hè mang lại, tiết kiệm hàng chục tỷ đồng ngân sách bảo trì luồng tàu.
2. Xây dựng và kéo dài đê chắn cát, đê ngăn triều tại cửa Nam Triệu và Lạch Huyện: Sở Giao thông Vận tải phối hợp cùng Ban Quản lý Dự án Hàng hải thực hiện xây dựng tuyến đê dẫn dòng định hướng trong giai đoạn 2026 - 2028. Mục tiêu kỹ thuật là duy trì vận tốc dòng chảy đáy qua luồng luôn vượt ngưỡng vận tốc lắng đọng 0,008 m/s, hạn chế tình trạng sa bồi tự nhiên.
3. Thiết lập hệ thống trạm quan trắc tự động đo độ đục và dòng chảy thời gian thực: Ủy ban nhân dân thành phố Hải Phòng và Sở Tài nguyên và Môi trường cần đầu tư lắp đặt mạng lưới phao quan trắc tự động tại 5 cửa sông chính trong vòng 24 tháng tới. Hệ thống cung cấp dữ liệu liên tục 24/7 với chu kỳ đo 10 phút một lần nhằm cảnh báo sớm hiện tượng bồi cạn luồng tàu sau các đợt bão lớn.
4. Tích hợp mô hình VNU/MDEC vào hệ thống quy hoạch không gian biển và cảng biển: Các viện nghiên cứu và cơ quan quản lý đô thị cần chuẩn hóa quy trình mô phỏng thủy động lực 3D cho 100% dự án lấn biển, xây dựng cầu cảng mới tại khu vực Đình Vũ và Cát Hải, đảm bảo không làm cản trở thoát lũ và không gây xói lở thứ cấp tại bờ biển Đồ Sơn.

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

Nội dung và phương pháp luận của công trình nghiên cứu mang lại giá trị thực tiễn cao cho 4 nhóm đối tượng chính:

• Cơ quan quản lý nhà nước và quản lý cảng biển: Cục Hàng hải Việt Nam, Cảng vụ Hàng hải Hải Phòng và các doanh nghiệp khai thác cảng tại Đình Vũ có thể ứng dụng kết quả để xây dựng kế hoạch duy tu luồng tàu, tối ưu hóa chi phí nạo vét định kỳ từ 15% đến 25%.
• Kỹ sư thiết kế công trình thủy công và ven biển: Các kỹ sư thuộc các viện thiết kế giao thông thủy sử dụng bộ thông số về trường vận tốc dòng chảy, độ cao sóng từ 0,64 m đến 5,6 m và ứng suất trượt đáy để tính toán kết cấu đê kè, cầu cảng và luồng tàu đạt chuẩn an toàn.
• Nhà khoa học và nghiên cứu sinh chuyên ngành Hải dương học: Các nhà nghiên cứu về động lực học vùng ven bờ có thể kế thừa thuật toán tọa độ sigma kép, sơ đồ lưới Arakawa C và kỹ thuật xử lý biên hở cửa sông có triều áp đảo để phát triển các mô hình thủy nhiệt động lực 3D cho các vùng biển khác.
• Giảng viên và học viên cao học ngành Khí tượng, Thủy văn, Môi trường biển: Tài liệu là tài liệu tham khảo mẫu mực về quy trình xây dựng kịch bản mô phỏng, phương pháp hiệu chỉnh mô hình số trị và phân tích tương tác sông - biển phức tạp trong các đồ án tốt nghiệp và luận văn chuyên ngành.

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

Tại sao khu vực cửa sông ven biển Hải Phòng lại có chế độ bồi xói phức tạp bậc nhất miền Bắc? Vùng biển Hải Phòng chịu sự tương tác đồng thời của 5 nhánh sông thuộc hệ thống sông Thái Bình và chế độ nhật triều đều có biên độ lớn tới 2,0 m. Mùa lũ tập trung tới 85% lượng nước và 90% lượng bùn cát cả năm, tạo nên sự tranh chấp động lực gay gắt giữa dòng xả sông và sóng triều biển.

Mô hình VNU/MDEC có ưu điểm gì vượt trội so với các mô hình 2D truyền thống? Mô hình VNU/MDEC sử dụng hệ tọa độ sigma kép với 5 tầng thẳng đứng và sơ đồ Arakawa C, giúp tính toán chính xác cấu trúc dòng chảy ba chiều và hiện tượng phân tầng mật độ tà áp. Điều này loại bỏ hoàn toàn các sai số biên tại vùng địa hình đáy dốc đứng và bãi triều lộ thiên.

Yếu tố nào đóng vai trò chủ đạo chi phối trường dòng chảy tại khu vực nghiên cứu? Dòng triều toàn nhật là yếu tố áp đảo lớn gấp 5 đến 10 lần dòng bán nhật, quy định hướng và độ lớn của dòng chảy tổng hợp với vận tốc từ 0,4 m/s đến 1,0 m/s. Dòng chảy sông chỉ đóng vai trò thứ cấp chi phối mạnh tại vùng cận cửa trong các đợt lũ lớn tháng 7 và 8.

Vận tốc lắng đọng và bứt xói của bùn cát đáy được xác định như thế nào trong mô hình? Mô hình áp dụng công thức của Van Rijn và Pohlmann với ngưỡng vận tốc lắng đọng tới hạn là 0,008 m/s và vận tốc động lực bứt xói tới hạn là 0,028 m/s. Khi ứng suất đáy vượt ngưỡng tới hạn, trầm tích đáy sẽ bị cuốn hút trở lại cột nước lơ lửng.

Làm thế nào để ứng dụng 19 phương án tính toán của luận văn vào thực tế vận hành cảng biển? Các phương án từ HP01 đến HP19 bao quát toàn bộ các tổ hợp khí tượng thủy văn bất lợi nhất trong năm. Đơn vị quản lý luồng tàu có thể tra cứu nhanh trường bồi lắng ứng với từng hướng gió mùa đông bắc hoặc đông nam để chủ động phân luồng tàu thuyền an toàn.

Kết luận

Luận văn thạc sĩ của tác giả Phạm Văn Tiến đã giải quyết trọn vẹn bài toán mô phỏng thủy thạch động lực học ba chiều tại vùng cửa sông ven biển Hải Phòng thông qua 5 đóng góp khoa học nổi bật:

• Ứng dụng thành công mô hình 3D VNU/MDEC với tọa độ sigma kép và sơ đồ phân tách mode hiện đại, mô phỏng ổn định trên miền tính có địa hình phức tạp;
• Xây dựng và thực thi hoàn chỉnh 19 kịch bản tính toán chi tiết, bao quát đầy đủ sự tương tác đa quy mô giữa thủy triều, dòng chảy lũ và gió mùa;
• Định lượng chính xác trường vận tốc dòng chảy tổng hợp từ 0,4 m/s đến 1,0 m/s cùng bức tranh phân bố trầm tích lơ lửng theo không gian và thời gian;
• Làm rõ cơ chế sa bồi cục bộ tại luồng Cửa Cấm và cửa Nam Triệu dưới tác động kết hợp của dòng xả sông mùa lũ và sóng triều khuấy đục đáy;
• Đề xuất hệ thống giải pháp công trình và phi công trình có tính ứng dụng thực tiễn cao cho ngành hàng hải và quy hoạch kinh tế biển địa phương.

Để triển khai các kết quả nghiên cứu vào thực tế, lộ trình từ 12 đến 36 tháng tới cần tập trung hoàn thiện mô hình dự báo trực tuyến và số hóa bản đồ sa bồi luồng tàu phục vụ trực tiếp cho các doanh nghiệp cảng biển. Các nhà quản lý cảng, kỹ sư công trình biển và các nhà nghiên cứu hãy tham khảo chi tiết toàn văn luận văn thạc sĩ này để làm chủ công nghệ mô phỏng số trị 3D và tối ưu hóa các giải pháp bảo vệ công trình ven bờ bền vững.

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

Chương 1. Tổng quan về mô hình thủy động lực và vận chuyển trầm tích 1. Các nghiên của nước ngoài Vận chuyển trầm tích được nghiên cứu từ rất sớm như ở Trung Quốc cổ đại, Lương Hà, Hy Lạp và Đế quốc La Mã. Nghiên cứu bằng phương pháp lý thuyết và thực nghiệm sớm nhất được thực hiện bởi nhà khoa học DuBuat (1738-1809) người Pháp.

Ông xác định vận tốc dòng chảy gây ra xói mòn đáy, trong đó có xem xét đến sự khác nhau của vật liệu đáy. DuBuat đã phát triển khái niệm ma sát trượt. Hagen (1797-1884) người Đức và Dupuit (1804-1866) người Pháp mô tả về chuyển động dọc theo đáy và chuyển động lơ lửng của trầm tích. Brahms (1753) đề xuất một công thức tính vận tốc tới hạn trên đáy với vật liệu là đá.

Công thức vận tải đáy đầu tiên dựa vào độ dốc và độ sâu được DuBoys (1847-1924) người Pháp đề xuất, Ông khái quát quá trình vận chuyển như chuyển động của các hạt trầm tích trong một loạt các lớp. Đến khoảng năm 1900, mô hình biến đổi đáy đầu tiên được Fargue (1827- 1910) người Pháp và Reynolds (1892-1912) người Anh xây dựng. Cơ sở nghiên cứu vận chuyển trầm tích trong các máng thí nghiệm được bắt đầu bởi Engels (1854- 1945) người Đức và Gilbert (1843-1918) người Mỹ. Lý thuyết vận chuyển trầm tích được viết bởi Forchheimer (1852-1933) và Schoklitsch (1888-1969) người Đức.

Đến năm 1914, phát triển phương trình tích ứng suất trượt đáy tới hạn (bắt đầu chuyển động của một hạt) theo chiều dọc của đáy dốc. Phương trình tương tự cho một hạt dừng chuyển động theo chiều ngang một đáy dốc được Leiner đề xuất năm 1912. Năm 1936, Shields có một đóng góp quan trọng liên quan đến ứng suất trượt đáy tới hạn cho sự khởi đầu chuyển động của các hạt trầm tích. Các đường cong được đề xuất gọi là đường cong “Shislds”.

Các nghiên cứu đầu tiên liên quan đến động lực học chất lỏng và vận chuyển bùn cát được thực hiện bởi Bagnold năm 1936, 1937. Đến năm 1950, Einstein và 9 LUAN VAN CHAT LUONG download : add [e-mail] các cộng sự nhờ vào sự phát triển của năng lực tính toán, biến các mô hình toán vận chuyển bùn cát thành một chủ đề quan trọng trong lĩnh vực khoa học ven biển. Apmann và Ralph R. Rumer nghiên cứu quá trình phát tán các hạt trầm tính do khuếch tán rối trong dòng chảy bất đồng nhất dựa trên mô hình toán.

Thí nghiệm được tiến hành trong một máng dài với 3 lớp trầm tích. Hệ số khuếch tán được xác định là một hàm của đặc trưng trầm tích và vận tốc dòng chảy. Một trong những nghiên cứu đầu tiên liên quan đến điều kiện bùn lỏng được thực hiện bởi Einstein và Chien năm 1955, hai quá trình kết bông và cố kết đáy đã được nghiên cứu. Tác giả đã nhận định rằng độ mặn tối thiểu 1‰ là giới hạn khởi đầu cho quá trình kết bông.

Odd và Owen, 1972 sử dụng mô hình 1D xem xét tốc độ xói mòn và lắng đọng dựa trên công thức đề xuất của Krone 1962 và Partheniades 1965. Smith và Kirby, 1989 đã ứng dụng các mô hình 1D để mô phỏng vận chuyển bùn cát và thay đổi hình thái quy mô lớn ở các sông De Vries, trong kênh thủy triều Dyer và Evans, mô phỏng quá trình hình thành “lutocline” ở các cửa sông Ross và Mehta. Năm 1971, O'Connor trình bày mô hình 2D tích phân theo độ sâu. Ariathurai và Krone, 1976 đã trình bày một mô hình phần tử hữu hạn áp dụng các yếu tố hình tam giác với một xấp xỉ bậc hai cho nồng độ và phương pháp trọng số thặng dư Galerkian.

Mô hình sử dụng các công thức cổ điển xác định quá trình xói mòn và lắng đọng trầm tích. Quá trình keo tụ được tính toán bằng cách xác định vận tốc chìm lắng trên mỗi phần tử lưới là một hàm của thời gian. Mulder và Udink 1991 áp dụng mô hình 2D cho cửa sông Western Scheldt có tính đến thủy triều và sóng gió. Mô hình giải một phương trình cân bằng tác động phổ, nội suy độ cao và chu kỳ sóng tính toán theo các thời kỳ triều khác nhau để xác định vận tốc quỹ đạo và thành phần ứng suất trượt đáy do sóng.

Sử dụng các công thức thực nghiệm để tính toán xói mòn và lắng đọng trầm tích và sử dụng các giá trị đồng nhất cho ứng suất trượt tới hạn của quá trình xói mòn, lắng đọng và vận tốc chìm lắng. 10 LUAN VAN CHAT LUONG download : add [e-mail] Li và cộng sự, 1994 phát triển mô hình 2DV tích hợp giữa mô hình thủy động lực học và mô hình vận chuyển bùn cát cho cửa sông Gironde nước Pháp, trong đó có sử dụng mô hình khép kín rối để tính hệ số nhớt rối và hệ số khuếch tán, mô hình có tính đến quá trình trao đổi trầm tính đáy. Năm 2002, Wen-Cheng Liu, Ming-Hsi Hsu và Albert Y. Kuo áp dụng mô hình hai chiều trung bình độ sâu nghiên cứu đặc điểm thủy động lực và vận chuyển bùn cát lơ lửng trong cửa sông của hệ thống sông Tanshui Rivers, Đài Loan.

Beckers, 1991, trong một nghiên cứu dòng chảy tổng hợp vùng biển Tây Địa Trung Hải trong điều kiện mùa đông điển hình bằng mô hình GHER-3D, cho rằng mô hình có thể khôi phục các quá trình vật lý và xu hướng chính của dòng chảy tổng hợp trong khu vực. Năm 1994, Beckers và cộng sự nghiên cứu thủy động lực học vùng biển Tây Địa Trung Hải bằng mô hình 3D. Trong nghiên cứu này, các tác giả đã sử dụng 2 mô hình: mô hình “metagnostic” (định hướng hệ thống) và mô hình chuẩn đoán (định hướng quá trình), được chạy đồng thời và có tính đến tương tác. Nghiên cứu chỉ ra quá cấu trúc và sự bất ổn định của dòng Algeria.

O'Connor và Nicholson, 1988 cung cấp một mô hình 3D đầy đủ, bao gồm một mô hình vận chuyển bùn lỏng, có tính đến sự kết bông và cố kết. Katopodi và Ribberink 1992 đã phát triển một mô hình tựa 3D cho vận chuyển bùn cát lơ lửng trên cơ sở của phương trình bình lưu khuếch tán cho dòng chảy và sóng, phân tích độ nhạy của các tham số sóng và dòng chảy. Năm 1994, Leonor Cancino và Ramiro Neves mô tả và ứng dụng hệ thống mô hình thuỷ động lực và vận chuyển trầm tích 3D (dạng nghiêng áp, sử dụng phương pháp sai phân hữu hạn). Mô hình thủy động lực dựa trên xấp xỉ thuỷ tĩnh và xấp xỉ Boussinesq, sử dụng tọa độ sigma kép cho chiều thẳng đứng với lưới so le và sơ đồ bán ẩn bậc hai.

Ngoài phương trình động lượng và phương trình liên tục, mô hình giải hai phương trình vận chuyển nhiệt độ, độ muối và một phương trình trạng 11 LUAN VAN CHAT LUONG download : add [e-mail] thái có tính đến hiệu ứng nghiêng áp. Mô phỏng quá trình vận chuyển trầm tích gắn kết được thực hiện bằng cách giải các phương trình bảo toàn, bình lưu - khuếch tán 3D, trong cùng một lưới sử dụng trong mô hình thủy động lực. Qúa trình cố kết, xói mòn và lắng đọng của trầm tích được biểu diễn bằng các công thức thực nghiệm. Các mô hình đã được thử nghiệm và hiệu chỉnh bằng cách mô phỏng dòng triều và vận chuyển bùn cát lơ lửng ở các cửa sông.

Hai ứng dụng ở cửa sông Western Scheldt (Hà Lan) và Gironde (Pháp) cho thấy sự phù hợp tốt giữa kết quả tính toán và đo đạc thực địa. Năm 2003, Changsheng Chen và Hedong Liu phát triển mô hình 3D tính hoàn lưu khu vực ven biển và cửa sông. Mô hình dựa trên hệ phương trình nguyên thủy 3 chiều gồm các phương trình động lượng, liên tục, nhiệt, muối, mật độ và sử dụng mô hình khép kín rối bậc 2,5 của Mellor và Yamada. Mô hình sử dụng hệ tọa độ chuyển đổi sigma cho phương thẳng đứng, phương ngang sử dụng lưới cấu trúc hình tam giác.

Mô hình toán sử dụng phương pháp sai phân hữu hạn, thể tích hữu hạn và phần tử hữu hạn. Mô hình đã được áp dụng cho biển Bột Hải, cửa sông Satilla River. Pandoe và Billy L. Edge ứng dụng mô hình ADCIRC-3D tính toán dòng chảy và vận chuyển bùn cát dọc bờ biển vịnh Mexico và dọc bờ biển Texas, kết quả cho thấy mô hình cho kết quả tốt khi áp dụng cho các khu vực cửa sông có độ dốc nhỏ.

Hu phát triển mô hình tính toán vận chuyển trầm tích cho vùng cửa sông Pearl River (vịnh Lingding). Mô hình sử dụng kỹ thuật tách để giải các phương trình chủ đạo: giải các số hạng bình lưu bằng phương pháp Eulerian-Lagrangian, sử dụng phương pháp phần tử hữu hạn cho các số hạng khuếch tán theo phương ngang và phương pháp sai phân hữu hạn cho số hạng khuếch tán theo phương thẳng đứng. Sơ đồ khép kín rối bậc 2,5 của Mellor- Yamada được sử dụng kết hợp để xác định tham số nhớt rối thẳng đứng. 12 LUAN VAN CHAT LUONG download : add [e-mail] Guy Simpsona, Sebastien Castelltort, 2005 trình bày mô hình coupled giữa mô hình dòng chảy mặt, vận chuyển trầm tích và diễn biến hình thái.

Mô hình sử dụng các phương trình nước nông cho dòng chảy, bảo toàn nồng độ trầm tích, hàm thực nghiệm cho ma sát đáy, xói mòn và lắng đọng. Quá trình xói mòn và lắng đọng được xử lý độc lập và tác động đến thông lượng trầm tích thông qua trao đổi vuông góc với biên đáy của dòng chảy. Harris và Hernan G. Arangoc phát triển mô hình 3D couple sóng, dòng chảy và vận chuyển bùn cát bằng công cụ MCT (Model Coupling Toolkit) và áp dụng tính toán cho vịnh Massachusetts.

Mô hình là sự kết hợp giữa mô hình hoàn lưu ven biển ROM v3.0 và mô hình tính sóng vùng nước nông SWAN. Ứng suất sóng 2 chiều được đưa vào phương trình động lượng, cùng với hiệu ứng của sóng mặt. Vận chuyển trầm tích được xem xét trong nhiều lớp, mỗi lớp có các đặc điểm riêng như đường kính hạt, mật độ, vận tốc lắng đọng, ứng suất tới hạn cho quá trình xói mòn. Vận chuyển trầm tích lơ lửng trong cột nước được tính giống thuật toán bình lưu khếch tán và bổ sung thuận toán giải theo chiều thẳng đứng mà không phụ thuộc vào tiêu chuển CFL.

Ngoài ra, còn có mô hình lớp biên đáy tính toán tương tác sóng - dòng chảy, làm tăng ứng suất đáy, tạo điều kiện thuận lợi cho vận chuyển trầm tích và làm tăng ma sát đáy, tạo ra tác động ngược trở lại dòng chảy. Năm 2008, Idris Mandang và Tetsuo Yanagi áp dụng mô hình 3D ECOMSED được phát triển bởi HydroQual (2002) vào tính toán vận chuyển trầm tích khu vực cửa sông Mahakam, phía Đông Kalimantan, Indonesia. Mô hình có sử dụng phép xấp xỉ Bousinesq và xấp xỉ thủy tĩnh.