Tổng quan nghiên cứu

Thành phố Nha Trang là trung tâm kinh tế - du lịch trọng điểm của khu vực miền Trung với tổng công suất cấp nước sạch đạt khoảng 125.000 m³/ngày đêm. Theo số liệu điều tra dân số, toàn thành phố có 396.066 người, trong đó mật độ dân cư khu vực nội thành lên đến gần 30.000 người/km², trong khi các xã ngoại thành có địa hình đồi núi dốc như xã Vĩnh Lương chỉ đạt mật độ khoảng 320 - 370 người/km². Sự phân bố dân cư không đồng đều cùng đặc thù địa hình phân tầng gây ra nhiều thách thức lớn trong công tác vận hành mạng lưới cấp nước và kiểm soát chất lượng nước sinh hoạt.

Vấn đề nghiên cứu cốt lõi xuất phát từ thực trạng suy giảm nồng độ hóa chất khử trùng trên mạng lưới phân phối. Quy chuẩn kỹ thuật quốc gia quy định hàm lượng Clo dư trong nước sinh hoạt phải duy trì trong ngưỡng từ 0,3 mg/L đến 0,5 mg/L. Tuy nhiên, tại khu dân cư Vĩnh Lương, số liệu quan trắc thực tế cho thấy nồng độ Clo dư ở đầu mạng lưới thường vượt mức 0,55 mg/L, nhưng khi vận chuyển đến các điểm cuối mạng lưới lại suy giảm xuống dưới 0,25 mg/L, tiềm ẩn nguy cơ tái nhiễm vi sinh vật gây hại cho sức khỏe người tiêu dùng.

Luận văn thạc sĩ chuyên ngành Kỹ thuật Tài nguyên nước tập trung vào mục tiêu xây dựng, cân chỉnh và kiểm định mô hình thủy lực kết hợp mô hình chất lượng nước cho hệ thống mạng lưới gồm 423 nút và 413 đoạn ống tại xã Vĩnh Lương, thành phố Nha Trang. Thông qua phần mềm WaterGEMS và công cụ tối ưu hóa tự động, nghiên cứu mô phỏng chính xác sự biến thiên nồng độ Clo dư theo thời gian mở rộng, đồng thời thiết lập các kịch bản dự báo khi xuất hiện sự cố rò rỉ hoặc thay đổi liều lượng châm Clo đầu nguồn. Kết quả nghiên cứu mang ý nghĩa thực tiễn cao, giúp đơn vị cấp nước tối ưu hóa quy trình khử trùng, giảm thiểu chi phí lắp đặt thiết bị đo tự động và đảm bảo 100% người dân ở các điểm bất lợi cuối nguồn đều được tiếp cận nguồn nước đạt chuẩn 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 thủy lực đường ống và động học chất lượng nước trong hệ thống phân phối:

  • Phương trình liên tục và bảo toàn năng lượng Bernoulli: Dòng chảy ổn định đều trong mạng lưới đường ống kín tuân theo định luật bảo toàn khối lượng và phương trình Bernoulli để xác định phân bố lưu lượng, cột áp thủy tĩnh và cột áp toàn phần tại các nút. Tổn thất dọc đường do ma sát được xác định qua công thức Hazen-Williams với hệ số nhám C, trong khi tổn thất cục bộ tại van và mối nối được tính toán qua hệ số trở lực riêng biệt.
  • Mô hình động học truyền tải và phân rã Clo dư: Sự biến đổi nồng độ Clo dư được mô tả thông qua phương trình truyền tải một chiều kết hợp phương trình cân bằng pha trộn hoàn toàn tại các nút mạng. Động học phân rã Clo tuân theo quy luật phản ứng bậc một, tích hợp đồng thời hai cơ chế: phản ứng phân hủy khối trong lòng dòng chảy ($k_b = -0,300\text{ ngày}^{-1}$) và phản ứng tương tác tại thành ống ($k_w$). Tốc độ phản ứng thành ống được liên kết chặt chẽ với hệ số nhám $C$, hệ số chuyển khối chất lỏng qua số Reynolds ($Re > 2.300$) và số Sherwood ($Sh \ge 2$).
  • Giải thuật di truyền (Genetic Algorithm - GA): Ứng dụng giải thuật tiến hóa Darwin trong module Darwin Calibrator để tự động giải bài toán nghịch đảo thủy lực, tối ưu hóa hệ số nhám ống và hệ số tiêu thụ nước nhằm cực tiểu hóa sai số giữa kết quả mô phỏng và số liệu quan trắc thực địa.

Các khái niệm chuyên ngành then chốt bao gồm: nồng độ Clo dư tự do, vùng phân vùng tách mạng (DMA), mô phỏng thời gian mở rộng (EPS) và bộ ghi nhận dữ liệu áp lực tự động (Datalogger).

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

  • Nguồn dữ liệu và tích hợp GIS: Tích hợp dữ liệu không gian và thuộc tính của 413 đoạn ống (vật liệu HDPE D33 - D300 PN16, tổng chiều dài hơn 26.000 m) từ bản đồ địa chính tỷ lệ 1:2.000 theo hệ tọa độ VN-2000. Dữ liệu nhu cầu dùng nước được gán tự động vào 423 nút mạng thông qua công cụ LoadBuilder từ dữ liệu tiêu thụ thực tế của 100% đồng hồ khách hàng (90% hộ dân sinh hoạt D15 và 10% cơ sở dịch vụ, hành chính).
  • Phương pháp chọn mẫu và quan trắc hiện trường: Thu thập dữ liệu áp lực và lưu lượng bằng phương pháp chọn mẫu phân tầng tại các điểm xung yếu. Bố trí 4 thiết bị datalogger đo áp lực liên tục 24/24 giờ tại các nút đại diện (nút J-2200, J-6228, J-8841, J-8103, J-1867) vào các ngày 15/4/2017 và 16/4/2017; kết hợp 1 đồng hồ tổng DMA và 1 trụ cứu hỏa Ø100 mm (TCH 68 đường Giáp Văn Cương) phục vụ xả nước thí nghiệm. Mẫu nồng độ Clo dư được thu thập định kỳ từ ngày 07/5/2017 đến 11/5/2017 tại bể chứa Vĩnh Lương (dung tích 1.000 m³) và các vị trí cuối nguồn bất lợi (nút J-13310 và J-11363).
  • Lý do lựa chọn phương pháp phân tích: Sử dụng phần mềm Bentley WaterGEMS kết hợp công cụ Constituent Analysis trong thời gian mở rộng 144 giờ giúp xử lý đồng thời tính phi tuyến phức tạp của mạng lưới cấp nước đồi núi và tái hiện trực quan động học chất lượng nước với độ chính xác cao.

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

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

Quá trình mô phỏng, cân chỉnh và kiểm định mô hình trên phần mềm WaterGEMS đã mang lại những phát hiện quan trọng:

  1. Độ chính xác cao trong cân chỉnh thủy lực: Sau khi tối ưu hóa hệ số nhám Hazen-Williams cho từng nhóm ống HDPE (hệ số $C$ dao động từ 130 đến 145), sai số bình phương trung bình (RMSE) giữa cột áp mô phỏng và số liệu quan trắc thực tế tại 5 vị trí datalogger giảm xuống dưới 0,5 - 1,0 m cột nước. Đồ thị áp lực ngày 15/4/2017 và ngày kiểm định 16/4/2017 đạt độ tương quan tuyến tính trên 95%.
  2. Xác định mức độ suy giảm Clo dư theo chiều dài mạng lưới: Nồng độ Clo dư giảm mạnh từ mức bình quân 0,58 mg/L tại bể chứa Vĩnh Lương xuống còn 0,22 - 0,26 mg/L tại nút cuối mạng J-13310 (giảm hơn 55%). Mô hình sau khi hiệu chỉnh hệ số phản ứng thành ống $k_w$ khớp hoàn toàn với số liệu kiểm nghiệm thực tế tại nút J-11363 với sai số nồng độ dưới 0,03 mg/L.
  3. Ảnh hưởng rõ rệt của sự cố rò rỉ đối với chất lượng nước: Khi mô phỏng kịch bản xuất hiện điểm bể gây rò rỉ lưu lượng qua hệ số emitter tại nút J-1183, vận tốc dòng chảy thay đổi làm xáo trộn thời gian lưu nước, khiến nồng độ Clo dư tại các nút hạ lưu lân cận giảm thêm 15% đến 22% so với điều kiện vận hành thông thường.
  4. Xác lập liều lượng châm Clo tối ưu: Kịch bản giả định điều chỉnh nồng độ Clo châm đầu nguồn tại bể chứa Vĩnh Lương lên mức 0,65 - 0,70 mg/L chứng minh 100% các nút trên toàn mạng lưới duy trì được hàm lượng Clo dư ổn định từ 0,30 đến 0,38 mg/L, triệt tiêu hoàn toàn vùng nước có nguy cơ mất an toàn sinh học.

Thảo luận kết quả

Sự suy giảm nồng độ Clo dư tại Vĩnh Lương chịu tác động chi phối bởi thời gian lưu nước (Water Age) kéo dài trên tuyến ống truyền tải chính HDPE D300 dài 5.000 m từ đèo Rù Rì và mạng phân phối nhánh D50, D33 có chiều dài hơn 14.000 m. Do nhu cầu sử dụng nước ban đêm giảm sâu, nước bị lưu cữu lâu trong đường ống, thúc đẩy phản ứng oxy hóa chất hữu cơ trong dòng chảy và tiếp xúc bề mặt thành ống.

Dữ liệu nghiên cứu được biểu diễn trực quan thông qua các đường cong biến thiên áp lực 24 giờ, đồ thị tương quan hồi quy giữa giá trị tính toán và thực đo, cùng bảng tổng hợp phân tích độ nhạy của hệ số phản ứng thành ống. Kết quả này tương đồng với các công bố quốc tế của J. Bohórquez (2016) tại Colombia và Bassam Halabi tại New Zealand, khẳng định phản ứng trên thành ống đóng vai trò then chốt trong quá trình phân rã chất khử trùng đối với mạng lưới cấp nước đô thị quy mô vừa và nhỏ.

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

Dựa trên kết quả mô phỏng và phân tích thực nghiệm, tác giả đề xuất 4 nhóm giải pháp cụ thể nhằm tối ưu hóa công tác quản lý chất lượng nước:

  1. Thiết lập quy trình châm Clo tự động theo lưu lượng tại bể chứa Vĩnh Lương: Điều chỉnh nồng độ Clo dư đầu nguồn dao động ổn định trong khoảng 0,65 - 0,70 mg/L tương ứng với từng khung giờ tiêu thụ cao điểm và thấp điểm. Giải pháp này giúp kiểm soát nồng độ Clo cuối nguồn luôn đạt trên 0,30 mg/L, triển khai trong vòng 3 tháng do Phòng Kỹ thuật và Vận hành Công ty Cổ phần Cấp thoát nước Khánh Hòa chủ trì thực hiện.
  2. Lắp đặt trạm châm Clo bổ sung (Booster Chlorination) tại phân vùng thôn Võ Tánh: Xây dựng 1 điểm châm bổ sung cục bộ tại khu vực có mật độ dân cư thưa thớt cuối tuyến để nâng nồng độ Clo dư từ mức <0,25 mg/L lên ngưỡng mục tiêu 0,32 - 0,35 mg/L. Timeline thực hiện trong 6 tháng do Ban Quản lý Dự án phụ trách.
  3. Định kỳ súc xả mạng lưới tại các tuyến ống cụt D33 và D50: Tận dụng trụ cứu hỏa TCH 68 đường Giáp Văn Cương và các van xả cặn cuối tuyến để thực hiện súc xả định kỳ 2 tuần/lần, giảm thời gian lưu nước quá 48 giờ và hạn chế tối đa sự hình thành màng sinh học bám thành ống. Hoạt động do Đội Quản lý Mạng lưới trực tiếp vận hành.
  4. Nâng cấp tích hợp mô hình WaterGEMS với hệ thống SCADA và GIS: Tự động hóa quá trình đồng bộ số liệu áp lực, lưu lượng từ 4 trạm datalogger về trung tâm điều hành để phát hiện sớm các điểm rò rỉ nước ngầm, hoàn thành trong lộ trình 12 tháng.

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

Công trình nghiên cứu mang tính ứng dụng thực tiễn cao, là tài liệu giá trị cho 4 nhóm đối tượng:

  1. Kỹ sư vận hành và quản lý mạng lưới cấp nước: Cung cấp quy trình từng bước ứng dụng module Darwin Calibrator và Constituent trong WaterGEMS để cân chỉnh mạng lưới cấp nước thực tế, giúp giảm 5% - 10% tỷ lệ thất thoát nước và tối ưu hóa hóa chất khử trùng.
  2. Chuyên gia quy hoạch và tư vấn thiết kế hạ tầng kỹ thuật: Tham khảo phương pháp phân vùng tách mạng DMA, tính toán thủy lực vùng địa hình đồi núi chênh cao áp lực lớn (4 - 6 PN) và thiết kế đường kính ống HDPE hợp lý nhằm hạn chế vùng nước chết.
  3. Học viên cao học và nghiên cứu sinh ngành Kỹ thuật Tài nguyên nước, Môi trường: Làm tài liệu tham khảo mẫu mực về phương pháp kết hợp giữa lý thuyết động học phân rã chất khử trùng, giải thuật di truyền và dữ liệu không gian GIS trong bài toán quản lý tài nguyên nước.
  4. Cơ quan quản lý nhà nước về y tế dự phòng và cấp nước an toàn: Cung cấp cơ sở khoa học định lượng để xây dựng quy chuẩn giám sát chất lượng nước sạch trực tuyến theo quy chuẩn QCVN 01-1:2018/BYT.

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

  • Tại sao cần cân chỉnh mô hình thủy lực trước khi xây dựng mô hình chất lượng nước? Mô hình chất lượng nước sử dụng trực tiếp trường vận tốc, lưu lượng và hướng dòng chảy từ mô hình thủy lực để tính toán quá trình vận chuyển và pha trộn chất tan. Nếu mô hình thủy lực chưa khớp với áp lực thực tế (sai số cột áp vượt quá 1,0 m), thời gian lưu nước tính toán sẽ sai lệch, dẫn đến dự báo nồng độ Clo dư không chính xác.

  • Giải thuật di truyền (GA) hỗ trợ cân chỉnh hệ số nhám Hazen-Williams như thế nào? Giải thuật di truyền trong Darwin Calibrator tự động tạo lập quần thể các giá trị hệ số nhám C giả định, thực hiện quá trình lai ghép và đột biến để tìm ra bộ thông số tối ưu nhất. Phương pháp này giúp cực tiểu hóa hàm mục tiêu sai số bình phương (RMSE) giữa áp lực đo thực tế tại 4 datalogger và áp lực tính toán nhanh hơn 80% so với phương pháp thử dần thủ công.

  • Những nguyên nhân chính gây suy giảm nồng độ Clo dư tại khu dân cư Vĩnh Lương là gì? Sự suy giảm nồng độ Clo dư (từ >0,55 mg/L xuống <0,25 mg/L) bắt nguồn từ thời gian lưu nước kéo dài trên mạng lưới đường ống dài hơn 26.000 m, địa hình đồi núi dốc khiến chế độ thủy lực biến động và phản ứng oxy hóa Clo diễn ra đồng thời trong khối chất lỏng lẫn bề mặt thành ống HDPE.

  • Phần mềm WaterGEMS có ưu điểm gì vượt trội so với EPANET trong nghiên cứu này? WaterGEMS cho phép liên kết trực tiếp với cơ sở dữ liệu GIS bản đồ địa chính tỷ lệ 1:2.000, tự động phân bổ lưu lượng tiêu thụ từ 100% đồng hồ khách hàng bằng công cụ LoadBuilder và tích hợp sẵn giải thuật di truyền tự động cân chỉnh tham số nhám mà EPANET tiêu chuẩn không hỗ trợ.

  • Mô phỏng kịch bản điểm bể đường ống mang lại lợi ích gì trong quản lý cấp nước? Mô phỏng điểm vỡ đường ống bằng hệ số rò rỉ emitter giúp xác định chính xác khu vực bị tụt áp dưới 2 PN và mức độ suy giảm Clo dư (giảm 15% - 22% tại hạ lưu). Nhờ đó, đơn vị vận hành có thể chủ động cô lập sự cố, xả cặn và châm tăng cường Clo kịp thời khi khắc phục sửa chữa ống.

Kết luận

  • Luận văn đã hoàn thành xuất sắc việc xây dựng và cân chỉnh đồng thời mô hình thủy lực và mô hình chất lượng nước cho mạng lưới cấp nước 423 nút tại xã Vĩnh Lương, đạt độ tương quan áp lực trên 95% và sai số nồng độ Clo dưới 0,03 mg/L.
  • Xác định rõ cơ chế phân rã Clo dư trên hệ thống đường ống HDPE chịu ảnh hưởng đồng thời bởi phản ứng khối ($k_b = -0,300\text{ ngày}^{-1}$) và phản ứng thành ống $k_w$, lý giải tình trạng suy giảm Clo xuống dưới 0,25 mg/L tại cuối nguồn.
  • Đề xuất thành công kịch bản điều chỉnh nồng độ Clo châm đầu nguồn tối ưu ở mức 0,65 - 0,70 mg/L tại bể chứa 1.000 m³, đảm bảo 100% người dân trong khu vực nhận được nguồn nước đạt quy chuẩn 0,3 - 0,5 mg/L.
  • Lộ trình triển khai các giải pháp châm Clo tự động, bố trí trạm booster và súc xả định kỳ qua trụ cứu hỏa TCH 68 được xác lập rõ ràng từ 3 đến 12 tháng.
  • Đóng góp quy trình công nghệ chuẩn hóa, có khả năng nhân rộng hiệu quả cho các mạng lưới cấp nước đô thị có điều kiện địa hình đồi núi tương tự trên toàn quốc.