Tổng quan nghiên cứu

Trong bối cảnh đô thị hóa và phát triển công nghiệp nhanh chóng, nhu cầu khai thác nước dưới đất phục vụ sinh hoạt và sản xuất tại các đô thị ước tính chiếm từ 50% đến 60% tổng công suất cấp nước. Tuy nhiên, chất lượng nguồn nước ngầm đang đối mặt với hiểm họa ô nhiễm nghiêm trọng do sự xâm nhập của các dòng nước thải chưa qua xử lý từ hệ thống sông ngòi và kênh rạch, nơi ghi nhận mức độ ô nhiễm gia tăng khoảng 15% đến 20% mỗi năm. Sự lan truyền các chất bẩn hòa tan từ đáy kênh vào tầng chứa nước không chỉ làm suy thoái trữ lượng nước sạch mà còn đe dọa trực tiếp đến an toàn sức khỏe cộng đồng.

Vấn đề cốt lõi đặt ra là việc kiểm soát và dự báo chính xác động thái lan truyền chất ô nhiễm trong môi trường xốp chịu tác động phức tạp của các yếu tố địa chất và thủy lực tự nhiên. Nghiên cứu tập trung giải quyết bài toán mô phỏng quá trình khuếch tán chất bẩn hòa tan hai chiều theo phương đứng từ kênh ô nhiễm vào tầng nước ngầm trong giai đoạn quá độ. Mục tiêu cụ thể là thiết lập mô hình toán số chính xác, đánh giá mức độ ảnh hưởng của độ rỗng môi trường thấm và chế độ dòng chảy ngầm ngang đến trường vận tốc, cột áp và nồng độ chất bẩn theo thời gian.

Phạm vi không gian nghiên cứu được xác lập trên miền thấm hai chiều có kích thước hình học 100m chiều rộng và 50m chiều sâu, với chu kỳ mô phỏng động học kéo dài 365 ngày liên tục. Nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp luận chứng tin cậy giúp các cơ quan quản lý giảm thiểu 25% đến 30% chi phí xử lý sự cố môi trường nước ngầm và tối ưu hóa quy hoạch mạng lưới cấp nước an toàn.

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 được xây dựng trên nền tảng cơ học chất lưu trong môi trường rỗng và lý thuyết truyền chất, tích hợp hai phương trình vi phân đạo hàm riêng phi tuyến:

  • Phương trình dòng thấm Richard: Áp dụng cho môi trường rỗng hai chiều theo phương thẳng đứng (y, z) ở trạng thái không ổn định. Phương trình mô tả sự biến thiên của cột áp toàn phần theo dung tích ẩm, độ bão hòa của đất và hệ số quan hệ thấm giữa trạng thái tự nhiên và bão hòa. Đặc tính thủy lực của đất được tích hợp qua mô hình phi tuyến Genuchten, xác định các thông số độ chứa nước bão hòa, độ chứa nước sau thoát nước và hệ số thấm bão hòa.
  • Phương trình lan truyền chất Advection - Dispersion: Mô tả quá trình vận chuyển khối lượng hai chiều của chất ô nhiễm hòa tan với nồng độ ban đầu $1.0\text{ kg/m}^3$. Phương trình tích hợp đồng thời các hiện tượng bình lưu do dòng thấm, khuếch tán phân tử, phân tán thủy động lực học qua tensor khuếch tán, hiện tượng hấp phụ tuyến tính qua hệ số phân bố và sự phân rã sinh - hóa trong đất và nước.

Các khái niệm then chốt bao gồm: trường vận tốc thấm Darcy hai chiều, thông lượng thấm biên, tensor khuếch tán phân tán thủy lực và hệ số cản làm chậm lan truyền chất trong lỗ rỗng.

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

Nghiên cứu ứng dụng phương pháp mô hình hóa toán số kết hợp phần tử hữu hạn thông qua phần mềm chuyên dụng COMSOL Multiphysics, đối chiếu kiểm chứng với nghiệm giải tích cổ điển:

  • Cỡ mẫu và phân chia ô lưới: Miền tính toán $100\text{m} \times 50\text{m}$ được rời rạc hóa thành 11.312 phần tử tam giác đều. Để đảm bảo độ chính xác cao tại các vị trí biến thiên gradient thủy lực lớn, kích thước ô lưới tại đáy kênh được làm mịn xuống 0,25m và tại các điểm góc kênh được chia nhỏ hơn 0,1m. Cỡ mẫu ô lưới lớn này đảm bảo triệt tiêu sai số làm tròn và duy trì tính hội tụ của bài toán phi tuyến.
  • Phương pháp chọn mẫu và điều kiện biên: Dữ liệu tính toán sử dụng thông lượng thấm đơn vị đáy kênh $0,00454\text{ m}^3/\text{s}$, mực nước ngập trong kênh cố định 3m, góc mái dốc kênh 34 độ, nồng độ chất ô nhiễm nguồn tại kênh $c_0 = 1,0\text{ kg/m}^3$. Điều kiện biên Dirichlet được gán cho cột áp và nồng độ đáy kênh; điều kiện biên Neumann không thấm và không thông lượng khuếch tán được gán cho các biên đáy và biên bên ngoài.
  • Lý do lựa chọn phương pháp: Phương pháp phần tử hữu hạn trong COMSOL cho phép giải đồng thời hệ phương trình liên kết phi tuyến dòng chảy và truyền chất trong môi trường dị tính không bão hòa, vượt trội hoàn toàn so với các phương pháp giải tích truyền thống vốn bị giới hạn ở các hình học đơn giản.
  • Timeline nghiên cứu: Mô hình mô phỏng quá trình quá độ từ trạng thái ban đầu chưa ô nhiễm tại thời điểm $t = 0$, theo dõi liên tục và trích xuất dữ liệu tại các mốc thời gian 120, 180, 240, 300 và 365 ngày.

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

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

Nghiên cứu đã định lượng hóa toàn diện quy luật biến đổi của trường thủy lực và trường nồng độ chất ô nhiễm:

  • Đặc tính trường vận tốc thấm: Vận tốc thấm đạt giá trị cực đại tại các góc đáy kênh hình thang với $0,472\text{ m/s}$ tại góc trái và $0,489\text{ m/s}$ tại góc phải, cao hơn khoảng 56% đến 62% so với vận tốc tại tim đáy kênh ($0,301\text{ m/s}$). Sai số giữa kết quả số và lời giải giải tích theo công trình của V. Vedernikov chỉ là 0,22%, khẳng định độ tin cậy tuyệt đối của mô hình.
  • Diễn biến cột áp theo không gian và thời gian: Cột áp đo áp toàn phần tăng đơn điệu theo chiều sâu và tăng dần theo thời gian mô phỏng. Tại các điểm quan trắc ở độ sâu -10m và -30m, cột áp ghi nhận mức tăng trưởng ổn định từ ngày thứ 120 đến ngày thứ 365, thể hiện sự hình thành thấu kính nước dâng dưới đáy kênh.
  • Động thái nồng độ chất ô nhiễm: Tại vị trí tim kênh ở độ sâu -10m, nồng độ chất ô nhiễm tăng nhanh và đạt đỉnh $0,196\text{ kg/m}^3$ sau 180 ngày trước khi giảm nhẹ do quá trình phân tán ngang. Ở tầng sâu -30m, đỉnh nồng độ xuất hiện muộn hơn vào khoảng ngày thứ 240 với giá trị $5,5 \times 10^{-4}\text{ kg/m}^3$, cho thấy độ trễ truyền chất theo phương đứng.
  • Hiệu ứng tương tác đa kênh: Trong kịch bản hai kênh ô nhiễm song song cách nhau 50m, hiện tượng giao thoa trường thấm làm gia tăng đột biến nồng độ chất bẩn. Tại điểm quan trắc biên (25, -10), nồng độ đỉnh đạt tới $0,225\text{ kg/m}^3$ tại ngày thứ 240, gấp hơn 10.000 lần so với trường hợp đơn kênh ($1,83 \times 10^{-5}\text{ kg/m}^3$).

Thảo luận kết quả

Cơ chế lan truyền chất bẩn chịu sự chi phối áp đảo của lực đẩy thủy lực tại giai đoạn đầu và quá trình khuếch tán phân tán ở giai đoạn sau. Sự tập trung vận tốc lớn tại góc kênh bắt nguồn từ gradient áp lực cực đại tại điểm uốn hình học của mặt cắt ướt hình thang. Khi độ rỗng môi trường địa chất tăng từ 0,2 lên 0,4, tốc độ di chuyển của ranh giới ô nhiễm tăng khoảng 35% đến 45%, đẩy nhanh thời gian chất độc hại thâm nhập vào các tầng chứa nước sâu.

Dữ liệu mô phỏng được biểu diễn trực quan thông qua các bản đồ đường đồng mức nồng độ hai chiều và các biểu đồ phân bố dọc theo các trục thẳng đứng 25-25, 50-50, 75-75 qua các mốc thời gian 120, 180, 240, 300, 365 ngày. So sánh với các công bố của Nguyễn Thống (2010) trên mô hình dòng Darcy ổn định và kết quả nghiên cứu của Shaul Sorek (1985), việc giải phương trình Richard phi tuyến trong điều kiện quá độ phản ánh sát thực hơn động thái hình thành vệt ô nhiễm dạng hình chuông mở rộng dưới đáy kênh. Khi xuất hiện dòng chảy ngầm tự nhiên theo phương ngang, trục đối xứng của đám mây ô nhiễm bị lệch góc từ 15 đến 30 độ về phía hạ lưu dòng ngầm, làm tăng diện tích tầng nước ngầm bị ảnh hưởng lên hơn 40%.

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

Dựa trên các quy luật lan truyền chất ô nhiễm đã được chứng minh, bốn nhóm giải pháp hành động cụ thể được kiến nghị thực hiện:

  • Thiết lập hành lang cách ly vệ sinh an toàn: Ban hành quy chuẩn bắt buộc cấm khai thác nước ngầm sinh hoạt trong phạm vi bán kính tối thiểu 75m tính từ mép các kênh mương hở dẫn nước thải. Cơ quan quản lý tài nguyên môi trường cần hoàn thành việc cắm mốc cảnh báo trong vòng 6 đến 12 tháng, hướng tới mục tiêu giảm 95% nguy cơ nhiễm độc nguồn nước sinh hoạt tại các vùng ven đô.
  • Gia cố lớp chống thấm nhân tạo tại đáy kênh: Triển khai thi công lớp lót địa kỹ thuật HDPE hoặc đệm sét nén dày tối thiểu 0,5m với hệ số thấm $K < 10^{-8}\text{ m/s}$ tại các đoạn kênh đi qua khu công nghiệp. Ban quản lý dự án công trình thủy lợi cần triển khai nâng cấp hệ thống kênh trong giai đoạn 12 đến 24 tháng tới nhằm cắt giảm hơn 85% lưu lượng thấm ô nhiễm vào lòng đất.
  • Thiết lập mạng lưới giếng quan trắc tự động tầng nông: Xây dựng hệ thống 5 đến 7 cụm giếng quan trắc đa tầng tại các độ sâu chiến lược -10m và -30m dọc theo hành lang các tuyến kênh trọng yếu. Trung tâm quan trắc môi trường thực hiện đo đạc tự động với tần suất 15 ngày/lần, đảm bảo phát hiện sớm sự biến động nồng độ vượt ngưỡng $0,05\text{ kg/m}^3$.
  • Ứng dụng mô hình số trong thẩm định quy hoạch đô thị: Tích hợp quy trình mô phỏng phần tử hữu hạn hai chiều trên COMSOL vào công tác đánh giá tác động môi trường của các dự án thoát nước và khu đô thị mới. Các viện quy hoạch xây dựng cần áp dụng bộ công cụ này trong lộ trình 2026-2030 để tối ưu hóa 100% vị trí bãi giếng khoan cấp nước tập trung.

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

Nội dung luận văn mang giá trị học thuật và ứng dụng thực tiễn cao cho nhiều nhóm chuyên gia:

  • Kỹ sư công trình thủy và kỹ thuật tài nguyên nước: Nắm vững phương pháp tính toán trường vận tốc thấm không ổn định qua mặt cắt kênh hình thang và áp dụng quy trình mô hình hóa số để thiết kế các công trình kênh dẫn, hồ chứa có biện pháp kiểm soát rò rỉ tối ưu.
  • Chuyên gia và cán bộ quản lý môi trường: Sử dụng dữ liệu lan truyền chất ô nhiễm và khoảng cách an toàn được định lượng trong luận văn làm căn cứ khoa học để xây dựng chính sách cấp phép xả thải và quy hoạch mạng lưới quan trắc môi trường ngầm cấp địa phương.
  • Nhà quy hoạch đô thị và hạ tầng kỹ thuật: Khai thác các kết quả tương tác thủy lực giữa kênh ô nhiễm và tầng nước ngầm để bố trí hợp lý khoảng cách giữa trục giao thông thủy, kênh thoát nước và các khu dân cư tập trung.
  • Học viên cao học, nghiên cứu sinh và giảng viên chuyên ngành: Tiếp cận tài liệu tham khảo chi tiết về thuật toán giải phương trình Richard phi tuyến kết hợp phương trình lan truyền chất trên phần mềm COMSOL Multiphysics, mở rộng cho các bài toán phân tán chất ô nhiễm phức tạp.

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

Mô hình 2D theo phương đứng có phản ánh đầy đủ thực tế hiện trường không?

Mô hình hai chiều theo phương thẳng đứng hoàn toàn phù hợp với các tuyến kênh dài có mặt cắt ướt và điều kiện địa chất đồng nhất dọc theo trục kênh. Sai số tính toán so với nghiệm giải tích dưới 0,22% chứng minh độ chính xác cao trong việc dự báo sự dịch chuyển chất bẩn theo phương ngang và phương thẳng đứng.

Tại sao vận tốc dòng thấm tại góc đáy kênh lại lớn hơn trục giữa kênh?

Do hình học góc vát đáy kênh 34 độ tạo ra sự co hẹp đường dòng và tập trung gradient thế năng áp lực lớn nhất tại các điểm gãy khúc. Vận tốc tại góc đạt $0,489\text{ m/s}$, cao hơn đáng kể so với mức $0,301\text{ m/s}$ tại đáy giữa kênh, khiến vị trí góc trở thành điểm xâm nhập ô nhiễm mạnh nhất.

Khoảng cách an toàn tối thiểu từ kênh ô nhiễm đến giếng khoan là bao nhiêu?

Theo kết quả mô phỏng trong vòng 365 ngày, nồng độ chất ô nhiễm suy giảm mạnh ngoài phạm vi cách tâm kênh 50m. Để đảm bảo an toàn dự phòng trước sự dao động của mực nước ngầm theo mùa, khoảng cách ly tối thiểu được khuyến nghị là từ 75m đến 100m.

Yếu tố địa chất nào ảnh hưởng mạnh nhất đến quá trình khuếch tán chất bẩn?

Độ rỗng hiệu dụng và hệ số thấm bão hòa của đất là hai yếu tố then chốt nhất. Độ rỗng môi trường càng lớn thì trở lực thấm càng giảm, làm gia tăng vận tốc bình lưu và thúc đẩy nồng độ chất hòa tan đạt đỉnh ở tầng sâu -30m nhanh hơn từ 30 đến 60 ngày so với nền đất sét chặt.

Phương trình Richard có ưu điểm gì vượt trội so với định luật Darcy truyền thống?

Định luật Darcy tuyến tính chỉ mô tả chính xác dòng thấm ổn định trong môi trường hoàn toàn bão hòa. Phương trình Richard phi tuyến giải quyết được bài toán dòng thấm biến đổi không ổn định trong đới bão hòa cục bộ và không bão hòa phía trên mực nước ngầm, giúp mô tả chính xác giai đoạn thấm quá độ từ đáy kênh.

Kết luận

  • Xây dựng thành công mô hình toán số hai chiều mô phỏng liên kết giữa phương trình thấm Richard phi tuyến và phương trình lan truyền chất Advection - Dispersion trong môi trường xốp quá độ.
  • Định lượng hóa chính xác trường vận tốc thấm với độ tập trung cực đại tại góc đáy kênh ($0,489\text{ m/s}$), kiểm chứng độ tin cậy qua sai số giải tích chỉ 0,22%.
  • Chứng minh hiện tượng tương tác giữa hai kênh song song làm nồng độ ô nhiễm tích tụ tăng vọt lên $0,225\text{ kg/m}^3$, đòi hỏi các tiêu chuẩn cách ly nghiêm ngặt hơn so với kênh đơn.
  • Xác lập quy luật phân bố nồng độ chất bẩn theo độ sâu và thời gian 365 ngày, khẳng định vai trò quyết định của độ rỗng địa chất và gradient dòng ngầm ngang.
  • Cung cấp bộ thông số đầu vào và cơ sở khoa học vững chắc phục vụ công tác quy hoạch, thiết kế lớp chống thấm đáy kênh và bảo vệ nguồn tài nguyên nước dưới đất.

Đóng góp lớn nhất của luận văn là hoàn thiện phương pháp mô phỏng số chi tiết hiện tượng lan truyền chất bẩn hòa tan từ công trình thủy lợi vào tầng ngầm, mở ra hướng ứng dụng hiệu quả cho ngành tài nguyên nước. Trong 12 đến 18 tháng tới, các nghiên cứu tiếp theo cần mở rộng mô hình sang không gian ba chiều và tích hợp các phản ứng phân hủy sinh học hiếu khí - kỵ khí. Các cơ quan quản lý và đơn vị tư vấn hạ tầng cần sớm áp dụng khung mô hình này vào thực tiễn để bảo vệ bền vững an ninh nguồn nước đô thị.