Tổng quan nghiên cứu

Bài toán mô phỏng thủy lực và lan truyền chất ô nhiễm trên hệ thống kênh sông là một thách thức kỹ thuật lớn đối với công tác bảo vệ nguồn nước và quy hoạch môi trường bền vững. Tại các đô thị lớn thuộc vùng kinh tế trọng điểm phía Nam, quá trình đô thị hóa và công nghiệp hóa nhanh chóng đã gây áp lực nghiêm trọng lên mạng lưới đường thủy. Thống kê thực tế cho thấy tổng lượng nước thải xả vào lưu vực kênh sông Sài Gòn ghi nhận khoảng 449.144 m³/ngày-đêm, trong đó nước thải sinh hoạt chiếm hơn 80% tổng dung lượng và chủ yếu mang đặc tính ô nhiễm hữu cơ. Dự báo tải lượng này có thể gia tăng mạnh mẽ lên tới 1.587.880 m³/ngày-đêm khi quy mô dân số và sản xuất mở rộng.

Trước thực trạng đó, công trình nghiên cứu tập trung giải quyết bài toán mô phỏng dòng chảy một chiều có truyền tải vật chất thụ động và bảo toàn trên quy mô mạng lưới kênh sông rộng lớn. Mục tiêu cốt lõi của luận văn là xây dựng hệ thống mô hình toán học tích hợp, phát triển giải thuật tính toán song song trên môi trường phân tán nhằm tối ưu hóa hiệu năng xử lý, đồng thời thiết kế cơ sở dữ liệu quan hệ liên kết với công cụ mô phỏng đồ họa không gian. Phạm vi không gian của nghiên cứu bao phủ toàn bộ mạng lưới sông Sài Gòn - Đồng Nai với cấu trúc hình học phức tạp gồm 321 nhánh sông, 192 nút giao cắt hợp lưu và 919 mặt cắt định vị. Nghiên cứu mang ý nghĩa thực tiễn sâu sắc, giúp rút ngắn thời gian xử lý dữ liệu quy mô lớn, nâng cao độ chính xác dự báo chỉ số nhu cầu oxy sinh hóa (BOD) và lượng oxy hòa tan (DO), tạo tiền đề công nghệ phục vụ giám sát tài nguyên nước và cảnh báo ô nhiễm theo thời gian thực.

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

Khung lý thuyết áp dụng

Hệ thống lý thuyết của nghiên cứu được xây dựng trên nền tảng thủy lực sông ngòi và động học lan truyền vật chất trong môi trường chất lỏng một chiều:

  • Hệ phương trình Saint-Venant một chiều: Đóng vai trò mô tả quy luật động lực học của dòng chảy không ổn định trên từng nhánh kênh sông. Hệ phương trình bao gồm phương trình liên tục bảo toàn khối lượng nước và phương trình động lượng tính đến lực quán tính, trọng trường, áp suất thủy tĩnh và hệ số nhám đáy kênh Manning dao động trong khoảng 0.015 đến 0.035 tùy mặt cắt.
  • Phương trình truyền tải - khuếch tán đối lưu (Advection-Diffusion Equation): Mô hình hóa sự dịch chuyển của nồng độ chất ô nhiễm bảo toàn dưới tác động tổng hòa của hai cơ chế: truyền tải đối lưu vĩ mô do vận tốc dòng nước mang theo và khuếch tán rối dọc trục dòng chảy với hệ số phân tán E.
  • Lý thuyết động học chuyển hóa BOD - DO: Mô tả sự suy giảm nồng độ oxy hòa tan do quá trình oxy hóa chất hữu cơ và khả năng tự làm sạch của sông ngòi. Cơ chế này được định lượng qua hằng số phân hủy BOD (k1) và hằng số tái sục khí bề mặt (k2). Tỷ số fs = k2/k1 xác định hệ số tự làm sạch của dòng sông, thường đạt giá trị ngưỡng xấp xỉ 2.0 trong điều kiện tự nhiên.
  • Tiêu chuẩn số vô thứ nguyên Peclet (Pe): Đo lường mức độ áp đảo tương đối giữa thông lượng truyền tải đối lưu và thông lượng khuếch tán phân tử, từ đó làm cơ sở lựa chọn sơ đồ sai phân lũy thừa phù hợp nhằm triệt tiêu hiện tượng khuếch tán số học sai lệch.

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

Nghiên cứu kết hợp chặt chẽ giữa toán học ứng dụng, mô hình hóa số trị và kỹ thuật tính toán song song hiệu năng cao:

  • Nguồn dữ liệu và kích thước không gian: Toàn bộ mạng lưới nghiên cứu sông Sài Gòn - Đồng Nai được số hóa chi tiết từ hệ thống bản đồ địa hình UTM tỷ lệ 1:10.000 do Cục Bản đồ Bộ Quốc phòng cung cấp. Kích thước tập mẫu bao gồm 321 nhánh, 192 nút hợp lưu cùng 919 mặt cắt hình học ngang. Khoảng cách đo đạc giữa hai mặt cắt thực tế được thiết lập tối ưu từ 2 km đến 5 km nhằm đảm bảo tính phản ánh trung thực về độ dốc, bề rộng và cao trình đáy.
  • Phương pháp chọn mẫu và cấu trúc lưới: Không gian mặt bằng kênh sông được phân rã thành hệ thống lưới các phần tử tứ giác và tam giác vector. Dữ liệu địa hình và thủy lực được chuẩn hóa thành 9 bảng quan hệ thực thể, hỗ trợ ràng buộc toàn vẹn dữ liệu đa cấp.
  • Lý do lựa chọn phương pháp phân tích số: Luận văn sử dụng phương pháp Thể tích hữu hạn (Finite Volume Method - FVM) tích hợp trên các ô kiểm soát dọc tuyến sông, kết hợp sơ đồ sai phân ẩn 4 điểm Preissmann với trọng số thời gian theta = 0.6667 (tương đương 2/3). Lựa chọn FVM mang lại ưu thế tuyệt đối về khả năng bảo toàn thông lượng khối lượng vật chất tại mọi bề mặt phân chia. Hệ phương trình lan truyền chất sau rời rạc hóa được đưa về dạng hệ ba đường chéo và giải trực tiếp bằng thuật toán Thomas (TDMA) với độ phức tạp tính toán tuyến tính O(n), tối ưu hóa vượt bậc so với các thuật toán khử Gauss O(n³) truyền thống.
  • Kiến trúc tính toán song song: Triển khai giải thuật Master-Slave trên nền tảng máy ảo song song PVM (Parallel Virtual Machine) dưới hệ điều hành Linux. Quá trình truy đuổi kép khử ẩn tại các nhánh và giải ma trận thưa tại các nút hợp lưu được phân chia song song theo luồng dữ liệu, giúp khai thác triệt để năng lực của hệ thống nhiều vi xử lý.

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

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

Quá trình chạy thử nghiệm và mô phỏng số trên hệ thống máy tính phân tán đã mang lại các phát hiện khoa học và kỹ thuật quan trọng:

  • Đặc tính phân bố nguồn thải hữu cơ: Phân tích dữ liệu thực tế tại 5 lưu vực kênh rạch chính ở TP.HCM (gồm Tham Lương - Bến Cát, Nhiêu Lộc - Thị Nghè, Tân Hóa - Lò Gốm, Kênh Đôi - Tàu Hũ và Bến Nghé) với tổng lưu lượng nước thải 44.289 m³/ngày-đêm cho thấy chất thải sinh hoạt chiếm tỷ trọng vượt trội 90%, trong khi nước thải công nghiệp chỉ chiếm khoảng 10%. Kênh Tham Lương - Bến Cát ghi nhận tải lượng nước thải lớn nhất đạt 20.333 m³/ngày-đêm.
  • Quy luật lan truyền chất theo chế độ bán nhật triều: Vận tốc dòng chảy và khả năng tự làm sạch của mạng kênh sông phụ thuộc chặt chẽ vào độ chênh mức nước triều Biển Đông. Vào mùa khô (tháng 3 - tháng 4), khi lưu lượng thượng nguồn từ hồ Dầu Tiếng chỉ đạt khoảng 20 m³/s và hồ Trị An khoảng 220 m³/s, lưu lượng đỉnh triều tại trạm Nhà Bè dao động mạnh từ 15.000 m³/s đến 20.000 m³/s trong tháng 3, nhưng sụt giảm mạnh xuống xấp xỉ 10.000 m³/s vào tháng 4. Hiện tượng chân triều cao đã gây ứ đọng và tạo ra các "vùng bùng nhùng" ô nhiễm cục bộ tại các kênh rạch nội thành do nước bẩn bị đẩy ngược trở lại.
  • Tối ưu hóa tài nguyên tính toán qua cấu trúc ma trận thưa: Hệ phương trình đại số tại các nút hợp lưu được tổ chức lưu trữ thông qua 4 mảng vector chỉ số linh hoạt, giúp cắt giảm hơn 75% nhu cầu bộ nhớ RAM và giảm đáng kể lưu lượng thông điệp truyền tải qua mạng cục bộ giữa tiến trình Master và các tiến trình Slave.
  • Hiệu suất tăng tốc song song: Việc phân rã 321 nhánh sông theo tiêu chí cân bằng tổng số mặt cắt trên các máy trạm slave đã mang lại hiệu quả cân bằng tải đạt trên 92%, giúp thời gian hoàn thành bước tính thủy lực và lan truyền chất giảm tỷ lệ thuận với số lượng máy tính tham gia cụm mạng.

Thảo luận kết quả

Kết quả mô phỏng số đã khẳng định tính tương thích hoàn hảo giữa mô hình toán học và hiện tượng thủy lực thực tế. Diễn biến nồng độ oxy hòa tan DO và nhu cầu oxy sinh hóa BOD được biểu diễn trực quan qua hệ thống đồ thị diễn tiến nồng độ theo chu kỳ 24 giờ triều và bảng đối sánh trường nồng độ tại 919 mặt cắt. Các số liệu cho thấy khi vận tốc dòng chảy đạt cực đại trong pha triều rút, hệ số tái sục khí k2 tăng lên rõ rệt, thúc đẩy nồng độ DO phục hồi nhanh chóng. Ngược lại, vào thời điểm triều đứng hoặc pha triều dâng, nồng độ BOD tích tụ cao vượt ngưỡng quy chuẩn môi trường tại các ngã ba sông hợp lưu.

So với các nghiên cứu thủy lực truyền thống vốn chỉ tính toán tuần tự trên từng đoạn sông đơn lẻ, giải pháp tính toán song song kết hợp mô hình thể tích hữu hạn trong luận văn này đã giải quyết triệt để bài toán quy mô lớn trên toàn lưu vực. Việc xây dựng module đồ họa 3D tương tác hai chiều trên môi trường Windows liên kết trực tiếp với cơ sở dữ liệu quan hệ cho phép người dùng kích chuột vào bất kỳ phần tử nào trên bản đồ để tra cứu ngay các thông số lưu lượng, mực nước và hàm lượng ô nhiễm tương ứng.

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

Dựa trên các kết quả đạt được từ mô hình tính toán và mô phỏng thực nghiệm, luận văn đưa ra 4 nhóm giải pháp công nghệ và quản lý cụ thể:

  • Nâng cấp hạ tầng tính toán song song thời gian thực: Triển khai mở rộng cụm máy chủ xử lý phân tán từ 8 đến 16 nút mạng hiệu năng cao dưới hệ điều hành Linux. Giải pháp này hướng tới mục tiêu rút ngắn chu kỳ dự báo ô nhiễm xuống dưới 15 phút cho mỗi kịch bản vận hành xả tràn, phục vụ đắc lực cho Trung tâm Điều hành Chống ngập và Quản lý Môi trường đô thị trong giai đoạn 1 đến 2 năm tới.
  • Mở rộng mạng lưới quan trắc và số hóa mặt cắt định kỳ: Sở Tài nguyên và Môi trường cần phối hợp với các cơ quan chuyên môn tiến hành đo đạc cập nhật thông số hình học cho toàn bộ 919 mặt cắt định kỳ 2 năm/lần. Đồng thời, lắp đặt bổ sung các trạm quan trắc tự động đo DO, BOD và độ mặn theo chu kỳ 1 giờ/lần tại các vị trí xung yếu cách nhau từ 2 km đến 5 km trên sông Sài Gòn và sông Đồng Nai.
  • Quy hoạch kiểm soát nguồn thải sinh hoạt tại nguồn: Các cơ quan quản lý đô thị cần đẩy nhanh tiến độ xây dựng các nhà máy xử lý nước thải tập trung cho 5 lưu vực kênh rạch trọng điểm, áp dụng định mức phát thải chuẩn 40g BOD/người/ngày. Đặt mục tiêu xử lý tối thiểu 80% tổng lưu lượng nước thải sinh hoạt (tương đương hơn 350.000 m³/ngày-đêm) trước khi xả ra mạng kênh sông chính trong kế hoạch 5 năm.
  • Tích hợp mô hình số vào hệ thống thông tin địa lý (WebGIS) dùng chung: Các viện nghiên cứu và đơn vị phát triển phần mềm cần tiến hành đóng gói cấu trúc cơ sở dữ liệu quan hệ 9 phân hệ và module mô phỏng dòng chảy lên nền tảng đám mây. Dự án này cần được hoàn thành trong vòng 18 tháng để phục vụ công tác chia sẻ dữ liệu liên tỉnh giữa TP.HCM, Bình Dương và Đồng Nai.

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

Luận văn là tài liệu tham khảo học thuật và ứng dụng công nghệ có giá trị cao đối với nhiều nhóm chuyên môn:

  • Cơ quan quản lý nhà nước về môi trường và thoát nước đô thị: Cung cấp bức tranh tổng thể và công cụ định lượng khoa học để đánh giá sức tải môi trường của mạng kênh sông, hỗ trợ ra quyết định cấp phép xả thải và điều tiết nguồn nước xả từ các hồ chứa thượng nguồn như Trị An và Dầu Tiếng.
  • Kỹ sư tin học tính toán và phát triển phần mềm: Là tài liệu tham khảo mẫu mực về kỹ thuật lập trình song song Master-Slave thông qua PVM, phương pháp giải hệ ma trận thưa quy mô lớn và kỹ thuật liên kết dữ liệu đồ họa vector với cơ sở dữ liệu quan hệ.
  • Viện nghiên cứu thủy lợi và môi trường sinh thái: Cung cấp khung phương pháp luận hoàn chỉnh về việc áp dụng phương pháp thể tích hữu hạn (FVM) và thuật toán Thomas giải bài toán truyền tải - khuếch tán đối lưu phục vụ mô phỏng lan truyền mặn, phèn và phú dưỡng hóa.
  • Giảng viên, học viên cao học và sinh viên chuyên ngành Công nghệ thông tin / Kỹ thuật tài nguyên nước: Là giáo trình thực hành chuyên sâu kết hợp liên ngành giữa khoa học máy tính và kỹ thuật công trình thủy, mở ra hướng nghiên cứu tối ưu hóa thuật toán số trên các hệ thống tính toán phân tán.

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

  • Tại sao phương pháp thể tích hữu hạn (FVM) lại được ưu tiên lựa chọn thay vì phương pháp sai phân hữu hạn thông thường?
    Phương pháp thể tích hữu hạn đảm bảo tính bảo toàn khối lượng và thông lượng vật chất một cách tuyệt đối trên từng ô kiểm soát, ngăn ngừa hiện tượng mất mát số học khi lan truyền nồng độ qua các đoạn kênh có tiết diện thay đổi đột ngột. Kết hợp với sơ đồ lũy thừa theo số Peclet, FVM duy trì tính ổn định số trị ngay cả khi tỷ số đối lưu chiếm ưu thế tuyệt đối so với khuếch tán.

  • Kiến trúc Master-Slave trong PVM xử lý bài toán thủy lực mạng kênh sông như thế nào?
    Tiến trình Master chịu trách nhiệm phân phối dữ liệu các nhánh sông cho các tiến trình Slave tính toán hệ số truy đuổi. Sau đó, Master tập hợp dữ liệu để thiết lập và cùng các Slave giải hệ phương trình mực nước tại 192 nút hợp lưu bằng giải thuật khử Gauss song song theo dòng, rồi phân phối kết quả ngược lại để tính nồng độ chất trên từng nhánh.

  • Chế độ bán nhật triều ảnh hưởng như thế nào đến sự tích tụ ô nhiễm hữu cơ tại TP.HCM?
    Chế độ bán nhật triều gây ra hiện tượng nước lên xuống 2 lần mỗi ngày. Khi triều dâng, dòng nước đẩy ngược chất ô nhiễm trở lại thượng lưu tạo nên các vùng nước ứ đọng có vận tốc dòng chảy gần bằng 0 m/s. Vào các tháng kiệt nhất (tháng 4), lưu lượng triều tại Nhà Bè giảm xuống 10.000 m³/s khiến quá trình tự làm sạch bị suy giảm nghiêm trọng.

  • Cơ sở dữ liệu của phần mềm giải quyết mối liên kết giữa dữ liệu hình học và dữ liệu tính toán ra sao?
    Hệ thống sử dụng cấu trúc quan hệ chuẩn hóa với lược đồ liên kết GTLINK đóng vai trò cầu nối trung gian giữa bảng đồ họa lưới tứ giác (COORD, VER, EDG) và bảng thuộc tính thủy lực (nhánh, mặt cắt, nút hợp lưu). Cơ chế này cho phép truy xuất và hiển thị động trường nồng độ màu sắc trực tiếp trên bản đồ số.

  • Mô hình toán học trong luận văn có thể mở rộng cho các dạng ô nhiễm nào khác ngoài BOD và DO?
    Khung phương trình truyền tải tổng quát trong luận văn hoàn toàn áp dụng được cho bài toán xâm nhập mặn (đo theo nồng độ độ mặn S), mô hình hóa lan truyền phèn (dựa trên nồng độ Nhôm Al, Sunfat Su và giá trị pH), cũng như bài toán phú dưỡng hóa nguồn nước do các hợp chất Nitơ và Phốt pho gây ra hiện tượng thủy triều đỏ.

Kết luận

  • Đóng góp học thuật và công nghệ cốt lõi:
    • Xây dựng thành công mô hình toán học tích hợp hoàn chỉnh giải bài toán thủy lực dòng chảy một chiều và lan truyền chất ô nhiễm bảo toàn trên mạng lưới kênh sông phức hợp.
    • Ứng dụng xuất sắc phương pháp thể tích hữu hạn (FVM) đưa hệ phương trình đối lưu - khuếch tán về dạng ma trận ba đường chéo, tối ưu hóa thuật toán Thomas (TDMA) đạt hiệu năng tuyến tính O(n).
    • Thiết kế và triển khai giải thuật tính toán song song phân tán Master-Slave qua môi trường máy ảo PVM trên hệ điều hành Linux, giải quyết triệt để bài toán cân bằng tải trên mạng lưới 321 nhánh và 919 mặt cắt.
    • Xây dựng hệ cơ sở dữ liệu quan hệ hoàn chỉnh gồm 9 bảng thực thể liên kết chặt chẽ với module mô phỏng đồ họa không gian vector 3D tương tác trên Windows.
    • Kiểm chứng và định lượng thành công động học lan truyền ô nhiễm hữu cơ (BOD, DO) dưới tác động của chế độ bán nhật triều Biển Đông trên hệ thống sông Sài Gòn - Đồng Nai.
  • Định hướng phát triển tiếp theo: Tích hợp các thuật toán học máy dự báo biên dòng chảy tự động, nâng cấp mã nguồn song song lên kiến trúc điện toán đám mây và kết nối mạng lưới cảm biến IoT thời gian thực trên toàn lưu vực sông ngòi Nam Bộ.
  • Khuyến nghị hành động: Các đơn vị nghiên cứu và cơ quan quản lý tài nguyên môi trường cần nhanh chóng ứng dụng khung giải thuật và cơ sở dữ liệu của công trình vào thực tiễn quản lý lưu vực, góp phần kiểm soát ô nhiễm và thúc đẩy phát triển kinh tế - xã hội bền vững.