Mô phỏng động lực học và vận chuyển trầm tích trong vịnh Delaware - Luận án tiến sĩ

Delaware Bay nằm giữa bang Delaware và New Jersey. Cửa sông có chiều dài khoảng 85 km từ cửa ra biển đến đầu nguồn. Độ sâu trung bình dao động từ 8-15 mét tại kênh chính. Dòng chảy triều chi phối chế độ thủy động lực với biên độ khoảng 1.5-2 mét. Lưu lượng nước ngọt từ sông Delaware đổ vào khoảng 340 m³/s trung bình. Địa hình đáy phức tạp với các kênh sâu và bãi cạn xen kẽ. Chế độ triều bán nhật triều với hai lần triều lên xuống mỗi ngày. Vùng nước mặn xâm nhập sâu vào nội địa theo mùa.

Mô hình số trị giúp dự báo chính xác dòng chảy và vận chuyển trầm tích. Các công cụ này hỗ trợ quy hoạch phát triển bền vững vùng ven biển. Nghiên cứu cung cấp thông tin cho quản lý chất lượng nước và sinh thái. Mô hình ECOMSED và mô hình Delft3D là các công cụ phổ biến trong lĩnh vực này. Hiểu biết về xói mòn đáy và lắng đọng giúp bảo vệ hệ sinh thái. Dữ liệu mô phỏng hỗ trợ thiết kế công trình thủy lợi và hàng hải. Độ đục nước ảnh hưởng trực tiếp đến sinh vật thủy sinh và chất lượng môi trường.

Nghiên cứu kết hợp quan trắc thực địa và mô phỏng số. Dữ liệu từ các trạm đo USGS được thu thập liên tục. Phân tích điều hòa triều xác định các thành phần dao động chính. Lưới tính toán được xây dựng bằng công cụ JANET chuyên dụng. Mô hình ba chiều giải quyết phương trình Navier-Stokes. Điều kiện biên được thiết lập từ số liệu quan trắc thực tế. Quá trình hiệu chỉnh đảm bảo độ tin cậy của kết quả mô phỏng.

Lưới trực giao giảm thiểu sai số số học trong tính toán. Tính trực giao đảm bảo độ chính xác cao cho gradient áp suất. Lưới không cấu trúc thích ứng tốt với địa hình tự nhiên. Chi phí tính toán giảm so với lưới cấu trúc có cùng độ chính xác. Khả năng làm mịn cục bộ giúp tập trung vào vùng quan trọng. Phương pháp phù hợp cho mô phỏng dòng chảy phức tạp ở cửa sông. Độ ổn định số học được cải thiện đáng kể.

GLS model thống nhất các mô hình k-ε, k-ω trong một khung toán học. Hai phương trình vận chuyển giải cho năng lượng rối và tỷ lệ tiêu tán. Các hằng số mô hình được điều chỉnh theo điều kiện cửa sông. Mô hình tái tạo tốt cấu trúc rối trong lớp biên đáy. Ứng suất Reynolds được tính từ giả thuyết độ nhớt rối. Phương pháp ổn định và cho kết quả tin cậy trong mô phỏng dài hạn. Khả năng mô phỏng phân tầng mật độ là điểm mạnh quan trọng.

JANET là phần mềm chuyên dụng tạo lưới không cấu trúc. Công cụ cho phép nhập dữ liệu địa hình từ nhiều nguồn khác nhau. Thuật toán tự động tối ưu hóa chất lượng lưới. Người dùng kiểm soát độ phân giải theo vùng quan tâm. Giao diện trực quan giúp điều chỉnh lưới dễ dàng. Lưới được xuất ra định dạng tương thích với các mô hình số. Quá trình tạo lưới nhanh chóng và hiệu quả cho vùng nghiên cứu lớn.

USGS vận hành mạng lưới trạm đo mực nước và dòng chảy. Trạm Cape May nằm ở cửa vào phía Nam Delaware Bay. Trạm Lewes đặt ở bờ phía Bắc gần cửa biển. Brandywine Shoal Light và Ship John Shoal Light nằm trong cửa sông. Dữ liệu được ghi liên tục với tần suất cao. Các trạm đo mực nước, vận tốc dòng chảy và độ mặn. Thông tin này quan trọng cho hiệu chỉnh và kiểm định mô hình.

Phân tích điều hòa phân tách tín hiệu triều thành các thành phần. Phương pháp bình phương tối thiểu xác định biên độ và pha. Các thành phần chính bao gồm M2, S2, N2, K1, O1. Tần số mỗi thành phần liên quan đến chu kỳ thiên văn. Nodal factor điều chỉnh biên độ theo vị trí Mặt Trăng. Equilibrium argument xác định pha ban đầu tại thời điểm chuẩn. Kết quả phân tích cho phép dự báo triều chính xác.

Dòng chảy dư là thành phần không dao động của dòng chảy. Nó được tính bằng cách lọc các thành phần triều từ dữ liệu. Dòng dư quan trọng cho vận chuyển vật chất và sinh vật. Phân bố không gian của dòng dư ảnh hưởng đến tuần hoàn nước. Vùng có dòng dư mạnh thường có độ đục nước cao. Thông tin này hỗ trợ nghiên cứu sinh thái và chất lượng nước. Mô hình tái tạo tốt cấu trúc dòng dư quan trắc được.

Xói mòn đáy xảy ra khi ứng suất đáy vượt ngưỡng tới hạn. Ứng suất này phụ thuộc vào vận tốc ma sát đáy. Công thức Shields xác định điều kiện khởi động chuyển động hạt. Kích thước hạt và mật độ ảnh hưởng đến ngưỡng xói. Trầm tích dính kết có ngưỡng xói cao hơn cát rời. Mô hình tính toán ứng suất đáy từ trường vận tốc. Vùng có dòng chảy mạnh thường xảy ra xói mòn tích cực.

Lắng đọng trầm tích xảy ra khi vận tốc dòng chảy giảm. Vận tốc lắng của hạt phụ thuộc vào kích thước và hình dạng. Công thức Stokes áp dụng cho hạt mịn trong điều kiện tầng. Độ đục nước đo lường nồng độ trầm tích lơ lửng. Vùng có độ đục cao ảnh hưởng đến sinh thái thủy sinh. Ánh sáng không xuyên qua được nước đục, hạn chế quang hợp. Mô hình dự báo phân bố độ đục theo không gian và thời gian.

Kết quả mô hình được kiểm định với dữ liệu thực địa. Phân bố trầm tích đáy từ mô hình phù hợp với khảo sát. Vùng lắng đọng chính tập trung ở kênh sâu và vùng nước tĩnh. Xói mòn mạnh xảy ra ở các bãi cạn có dòng chảy nhanh. Độ đục nước mô phỏng tương quan tốt với đo đạc. Sai số trung bình nằm trong giới hạn chấp nhận được. Mô hình có thể sử dụng cho dự báo và quy hoạch.

ECOMSED viết tắt của Estuarine Coastal Ocean Model with Sediment. Mô hình phát triển từ POM (Princeton Ocean Model) nổi tiếng. Lưới sigma theo địa hình thích hợp cho vùng nước nông. Phương trình liên tục và động lượng được giải bằng phương pháp mode splitting. Mô đun trầm tích tính toán xói, vận chuyển và lắng đọng. Nhiều lớp trầm tích đáy được mô phỏng với độ cứng khác nhau. Mô hình đã được áp dụng rộng rãi cho các cửa sông Bắc Mỹ.

Delft3D bao gồm nhiều mô đun cho các quá trình khác nhau. Mô đun FLOW giải quyết thủy động lực học hai và ba chiều. Mô đun SED mô phỏng vận chuyển trầm tích không dính kết. Mô đun MUD xử lý trầm tích dính kết và bùn. WAQ mô đun chất lượng nước với nhiều thông số sinh địa hóa. WAVE tính toán sóng và tương tác sóng-dòng chảy. Các mô đun tích hợp chặt chẽ qua giao diện chuẩn.

Việc lựa chọn mô hình phụ thuộc vào mục tiêu nghiên cứu. ECOMSED phù hợp cho nghiên cứu tổng hợp thủy động và trầm tích. Delft3D linh hoạt hơn với nhiều mô đun mở rộng. Cả hai mô hình đều có khả năng mô phỏng Delaware Bay. Dữ liệu đầu vào và điều kiện biên cần chuẩn bị kỹ lưỡng. Quá trình hiệu chỉnh đòi hỏi thời gian và dữ liệu quan trắc. Kết quả từ cả hai mô hình có thể so sánh để tăng độ tin cậy.

Luận án phát triển mô hình ba chiều cho Delaware Bay. Lưới tính toán được tối ưu hóa cho vùng nghiên cứu rộng lớn. Mô hình tái tạo chính xác dòng chảy triều và dòng dư. Quá trình vận chuyển trầm tích được mô phỏng chi tiết. Kết quả được kiểm định với dữ liệu từ nhiều nguồn. Công cụ này hỗ trợ hiểu biết sâu về động lực học cửa sông. Phương pháp có thể áp dụng cho các cửa sông khác.

Mô hình chưa tích hợp đầy đủ quá trình sinh địa hóa. Dữ liệu trầm tích đáy còn hạn chế ở một số vùng. Hiệu chỉnh các thông số mô hình đòi hỏi nhiều thời gian. Độ phân giải lưới cần tăng để mô phỏng chi tiết hơn. Tương tác sóng-dòng chảy chưa được xem xét đầy đủ. Ảnh hưởng của thực vật thủy sinh cần nghiên cứu thêm. Chi phí tính toán cao cho mô phỏng dài hạn.

Tích hợp mô đun chất lượng nước để đánh giá ô nhiễm. Nghiên cứu tác động của biến đổi khí hậu lên động lực cửa sông. Phát triển hệ thống dự báo thời gian thực cho quản lý. Áp dụng mô hình cho kịch bản quy hoạch và phát triển. Nghiên cứu vận chuyển ô nhiễm và vi sinh vật. Đánh giá tác động của nạo vét lên môi trường trầm tích. Kết hợp dữ liệu viễn thám để cải thiện độ chính xác mô hình.