Tổng quan nghiên cứu

Quá trình công nghiệp hóa và đô thị hóa nhanh chóng tại tỉnh Bắc Ninh đã tạo ra áp lực nặng nề lên môi trường tự nhiên, đặc biệt là nguồn tài nguyên nước ngầm. Tỉnh hiện có 62 làng nghề cùng 15 khu công nghiệp tập trung, trong đó 10 khu công nghiệp đã đi vào vận hành. Đi kèm với sự tăng trưởng kinh tế là tình trạng khai thác và xả thải thiếu kiểm soát: kết quả thanh tra cho thấy có tới 80% đơn vị (8/10 doanh nghiệp được kiểm tra) không thực hiện đầy đủ việc quan trắc, giám sát mực nước và chất lượng nước ngầm theo Quyết định số 15/2008/QĐ-BTNMT, trong khi 60% đơn vị chưa thực hiện báo cáo định kỳ. Bên cạnh đó, mỗi ngày toàn tỉnh phát sinh khoảng 500 tấn chất thải sinh hoạt và 200 tấn chất thải công nghiệp, khiến nguy cơ suy thoái và ô nhiễm tầng chứa nước ngày càng gia tăng.

Trước thực trạng quản lý hồ sơ dạng giấy tờ phân tán và thiếu công cụ trực quan hóa không gian, luận văn tập trung giải quyết bài toán: ứng dụng công nghệ Hệ thống Thông tin Địa lý (GIS) nhằm xây dựng hệ thống quản lý, phân tích và cảnh báo chất lượng nước ngầm trên địa bàn tỉnh Bắc Ninh. Mục tiêu cụ thể của nghiên cứu là thiết lập cơ sở dữ liệu không gian tích hợp, ứng dụng các thuật toán nội suy không gian để biểu diễn hiện trạng phân bố các chỉ số ô nhiễm (như hàm lượng Sắt - Fe, Mangan - Mn, pH) và xây dựng phần mềm hỗ trợ quản trị trực quan. Phạm vi dữ liệu được thu thập và phân tích từ 20 trạm quan trắc nước ngầm phân bố trên địa bàn tỉnh trong năm 2014. Nghiên cứu mang ý nghĩa thực tiễn to lớn khi giúp các nhà quản trị môi trường rút ngắn trên 60% thời gian tra cứu hồ sơ địa chất thủy văn, đồng thời nâng cao độ chính xác trong công tác khoanh vùng ô nhiễm và quy hoạch khai thác bền vững.

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 lý thuyết Hệ thống Thông tin Địa lý (GIS) kết hợp với các mô hình toán địa thống kê (Geostatistics) trong quản lý môi trường. GIS là hệ thống máy tính tích hợp phần cứng, phần mềm, cơ sở dữ liệu và tri thức chuyên gia, thực hiện 5 nhóm chức năng cốt lõi: thu thập, xử lý sơ bộ, lưu trữ - truy vấn, tìm kiếm - phân tích không gian và hiển thị đồ họa.

Trong cấu trúc dữ liệu không gian, nghiên cứu ứng dụng 3 mô hình cơ bản:

  • Mô hình Vector: Biểu diễn các đối tượng địa lý rời rạc gồm điểm (Point - vị trí 20 trạm quan trắc), đường (Line - hệ thống sông ngòi, mạng lưới giao thông) và vùng (Polygon - ranh giới hành chính các huyện, thị xã).
  • Mô hình Raster: Biểu diễn các bề mặt liên tục thông qua ma trận các ô lưới (grid cells), được sử dụng để hiển thị các bản đồ chuyên đề về nồng độ ô nhiễm phân tán sau khi nội suy.
  • Mô hình Mạng lưới tam giác bất quy tắc (TIN): Hỗ trợ phân tích địa hình và mô hình hóa bề mặt 3 chiều.

Bên cạnh đó, các phép phân tích không gian chuyên sâu như vùng đệm (Buffer), phân tích mạng (Network Analysis) và kỹ thuật chồng lớp bản đồ (Overlay với các toán tử Union, Intersect, Identity) được vận dụng để liên kết dữ liệu thuộc tính chất lượng nước với vị trí địa lý. Tiêu chuẩn đánh giá chất lượng nước ngầm trong nghiên cứu dựa trên Quy chuẩn kỹ thuật quốc gia QCVN 09:2008/BTNMT.

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

Nguồn dữ liệu thực nghiệm của luận văn bao gồm tọa độ không gian và kết quả phân tích hóa lý của 20 trạm quan trắc nước ngầm (ký hiệu từ G1 đến G20) trên địa bàn tỉnh Bắc Ninh qua 2 đợt đo đạc trong năm 2014. Phương pháp chọn mẫu là chọn mẫu mạng lưới không gian có chủ đích, bao phủ các khu vực trọng điểm gồm đô thị (Thành phố Bắc Ninh, Thị xã Từ Sơn), khu công nghiệp tập trung (Yên Phong, Tiên Sơn, Quế Võ) và các vùng phụ cận sông Cầu, sông Đuống nhằm đảm bảo tính đại diện địa chất thủy văn.

Để xây dựng bề mặt phân bố ô nhiễm từ 20 điểm đo thực tế, nghiên cứu tiến hành khảo sát và so sánh 3 thuật toán nội suy không gian phổ biến:

  • Thuật toán nội suy nghịch đảo khoảng cách có trọng số (Inverse Distance Weighting - IDW): Ước tính giá trị điểm chưa biết bằng cách tính trung bình trọng số nghịch đảo khoảng cách từ các điểm đã biết.
  • Thuật toán Spline: Sử dụng các đa thức bậc thấp để làm mịn bề mặt uốn cong qua các điểm mẫu với tổng độ cong tối thiểu.
  • Thuật toán Kriging: Phương pháp địa thống kê dựa trên biểu đồ bán phương sai (Semi-variogram) kết hợp mối tương quan không gian giữa các mẫu.

Độ chính xác và độ tin cậy của thuật toán được kiểm định thông qua hệ số xác định $R^2$ (biến thiên từ 0 đến 1) và chỉ số hiệu quả Nash–Sutcliffe (NSI, chạy từ âm vô cực đến 1). Sau quá trình đánh giá đối sánh, phương pháp IDW được lựa chọn làm giải pháp chính nhờ tính ổn định cao, thời gian xử lý nhanh và khả năng hạn chế hiện tượng dao động biên cực trị khi mật độ điểm mẫu ở mức 20 trạm. Toàn bộ cơ sở dữ liệu và thuật toán được lập trình, cấu hình trên bộ công cụ ArcGIS Desktop (ArcMap, ArcCatalog, ArcToolbox) kết hợp phần mềm tính toán ma trận Matlab R2015a trong giai đoạn 2014-2017.

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

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

Phân tích dữ liệu hiện trạng môi trường và kết quả mô phỏng GIS trên địa bàn tỉnh Bắc Ninh đã ghi nhận những phát hiện định lượng nổi bật:

Thứ nhất, mức độ ô nhiễm tại các làng nghề sản xuất ở mức báo động. Khảo sát chất lượng nước tại các khu vực làng nghề truyền thống và mới cho thấy các chỉ số nhu cầu oxy hóa học (COD), nhu cầu oxy sinh hóa (BOD5), oxy hòa tan (DO) và pH đều vượt giới hạn cho phép từ 2,5 đến 5,3 lần so với quy chuẩn môi trường hiện hành.

Thứ hai, kết quả phân tích nồng độ kim loại nặng tại 20 trạm quan trắc (G1 - G20) trong 2 đợt năm 2014 chỉ ra sự tích tụ Sắt (Fe) và Mangan (Mn) vượt ngưỡng giới hạn theo QCVN 09:2008/BTNMT tại nhiều khu vực dân cư và công nghiệp. Việc áp dụng thuật toán IDW đã tạo ra các bản đồ chuyên đề phân lớp (Classify) phản ánh rõ nét ranh giới vùng ô nhiễm tập trung cao tại các khu vực phụ cận Yên Phong và Quế Võ.

Thứ ba, công tác chấp hành quy định pháp luật về tài nguyên nước của các cơ sở sản xuất còn rất yếu kém. Kết quả thanh tra thực tế 10 cơ sở trên địa bàn ghi nhận: 10% cơ sở không có giấy phép xả nước thải vào nguồn nước, 10% không thông báo chuyển đổi tên doanh nghiệp, 60% không nộp báo cáo định kỳ tình hình khai thác và đặc biệt có tới 80% cơ sở (8/10 đơn vị) hoàn toàn bỏ qua hoặc thực hiện không đầy đủ quy định quan trắc giám sát lưu lượng và chất lượng nước dưới đất.

Thứ tư, lượng phát thải chất thải rắn và nước thải không ngừng gia tăng gây áp lực thẩm thấu trực tiếp xuống tầng nước ngầm. Mỗi ngày toàn tỉnh phát sinh 500 tấn rác thải sinh hoạt (riêng TP Bắc Ninh và TX Từ Sơn chiếm trên 200 tấn, tương đương 40%), 200 tấn chất thải công nghiệp và 2 tấn rác thải y tế. Tốc độ gia tăng chất thải rắn công nghiệp đạt mức 15%/năm, sinh hoạt tăng 10%/năm và y tế tăng 8%/năm, trong khi chỉ có 5/10 khu công nghiệp đang hoạt động có hệ thống xử lý nước thải tập trung hoàn thiện.

Thảo luận kết quả

Nguyên nhân cốt lõi dẫn đến tình trạng ô nhiễm nước ngầm cục bộ tại Bắc Ninh bắt nguồn từ sự phát triển công nghiệp nóng nhưng thiếu đồng bộ về hạ tầng môi trường. Việc 50% khu công nghiệp đang hoạt động chưa hoàn thiện trạm xử lý nước thải tập trung, kết hợp với hơn 62 làng nghề xả thải trực tiếp ra môi trường xung quanh, đã biến nước mặt và đất thành nguồn phát tán ô nhiễm ngấm sâu vào các tầng chứa nước ngầm nông. Bên cạnh đó, bộ máy quản lý nhà nước cấp huyện và xã còn mỏng; đa số cấp xã chỉ có 1 cán bộ địa chính kiêm nhiệm nhiều lĩnh vực, gần như không có chuyên môn sâu về tài nguyên nước, khiến các hành vi vi phạm không được ngăn chặn kịp thời.

Về mặt kỹ thuật hiển thị, dữ liệu nghiên cứu được biểu diễn sinh động qua hệ thống bản đồ Raster đa lớp kết hợp với các bảng thuộc tính chi tiết trong ArcMap. Người quản lý có thể thực hiện truy vấn tức thời một trạm quan trắc cụ thể (chẳng hạn như trạm G11) để hiển thị đồng thời tọa độ không gian, biểu đồ cột so sánh nồng độ Fe, Mn giữa Đợt 1 và Đợt 2 so với ngưỡng chuẩn QCVN 09:2008/BTNMT. Phương pháp hiển thị dải màu phân cấp (Classify) giúp phân biệt trực quan các vùng có nồng độ kim loại an toàn (màu xanh lá) và vùng vượt ngưỡng nguy hiểm (màu đỏ cam).

Khi so sánh với các nghiên cứu mô hình hóa khác, việc sử dụng thuật toán IDW trong phạm vi 20 điểm quan trắc tại Bắc Ninh chứng minh tính phù hợp vượt trội so với Spline và Kriging. Thuật toán Spline có xu hướng tạo ra các giá trị ngoại lai cực đại hoặc cực tiểu vượt quá giới hạn đo thực tế tại các vùng chuyển tiếp địa hình, trong khi Kriging đòi hỏi xây dựng hàm biến thiên bán phương sai phức tạp và tốn nhiều tài nguyên tính toán nhưng không mang lại sự khác biệt đáng kể về độ chính xác trên một diện tích tự nhiên tương đối bằng phẳng như tỉnh Bắc Ninh. Do đó, mô hình IDW tích hợp trong ArcGIS Desktop là lựa chọn cân bằng tối ưu giữa độ chính xác và hiệu quả ứng dụng thực tiễn.

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

Nhằm tăng cường hiệu lực quản lý nhà nước và bảo vệ tài nguyên nước ngầm tỉnh Bắc Ninh, nghiên cứu đề xuất 4 nhóm giải pháp trọng tâm:

  1. Chuẩn hóa và nâng cấp mạng lưới quan trắc GIS tự động: Sở Tài nguyên và Môi trường tỉnh Bắc Ninh cần phê duyệt đề án mở rộng mạng lưới quan trắc từ 20 trạm hiện tại lên tối thiểu 50 trạm quan trắc tự động truyền dữ liệu theo thời gian thực (IoT kết hợp GIS) trước năm 2028. Hệ thống cần tích hợp cơ sở dữ liệu không gian Geodatabase để theo dõi liên tục mực nước và nồng độ chất ô nhiễm tại 100% các khu, cụm công nghiệp.

  2. Thắt chặt kiểm tra và xử phạt vi phạm xả thải: Chi cục Bảo vệ Môi trường và Thanh tra liên ngành cần triển khai kiểm toán môi trường định kỳ tại 15 khu công nghiệp và 62 làng nghề. Mục tiêu trong vòng 24 tháng tới là kéo giảm tỷ lệ cơ sở vi phạm nghĩa vụ quan trắc từ mức 80% hiện nay xuống dưới 15%, đồng thời xử lý triệt để 100% các cơ sở chưa có giấy phép xả thải.

  3. Hoàn thiện hạ tầng xử lý chất thải và nước thải tập trung: Ban Quản lý các Khu công nghiệp tỉnh Bắc Ninh cần yêu cầu 5 khu công nghiệp còn lại khẩn trương hoàn thành trạm xử lý nước thải tập trung đạt chuẩn trước quý IV năm 2027. Đối với các làng nghề tại Yên Phong, Từ Sơn, cần bố trí vốn ngân sách xây dựng khu xử lý nước thải tập trung với công suất đáp ứng xử lý toàn bộ 200 tấn chất thải công nghiệp và nước thải phát sinh mỗi ngày.

  4. Đào tạo nâng cao năng lực ứng dụng công nghệ cho cán bộ cơ sở: Ủy ban nhân dân tỉnh Bắc Ninh cần phối hợp với các cơ sở đào tạo tổ chức các khóa tập huấn chuyên sâu về phần mềm ArcGIS và hệ thống thông tin địa lý cho 100% cán bộ phòng Tài nguyên và Môi trường cấp huyện và công chức địa chính - môi trường cấp xã trước năm 2027, đảm bảo vận hành thành thạo các chức năng tra cứu, cập nhật dữ liệu không gian.

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

Nội dung và kết quả nghiên cứu của luận văn mang lại giá trị học thuật và ứng dụng thiết thực cho các nhóm đối tượng sau:

  • Cơ quan quản lý nhà nước về môi trường: Lãnh đạo và chuyên viên Sở Tài nguyên và Môi trường, Chi cục Bảo vệ Môi trường tỉnh Bắc Ninh và các địa phương lân cận có thể sử dụng khung cơ sở dữ liệu GIS của đề tài để số hóa hồ sơ quản lý tài nguyên nước, hỗ trợ ra quyết định cấp phép khai thác và khoanh định hành lang bảo vệ nguồn nước.
  • Ban quản lý khu công nghiệp và doanh nghiệp hạ tầng: Các đơn vị quản lý 15 khu công nghiệp tại Bắc Ninh có thể áp dụng mô hình phân tích không gian để xác định vị trí tối ưu cho các giếng quan trắc nội bộ, theo dõi biến động tầng chứa nước và lập phương án phòng ngừa sự cố ô nhiễm.
  • Nhà nghiên cứu, giảng viên và sinh viên ngành Khoa học máy tính, Địa tin học (GIS), Môi trường: Tài liệu cung cấp quy trình chi tiết từ thu thập, chuẩn hóa tọa độ trạm quan trắc, so sánh thuật toán nội suy (IDW, Spline, Kriging) đến xây dựng giao diện hiển thị trên ArcGIS và Matlab, đóng vai trò là case study thực nghiệm mẫu mực cho các đề tài tốt nghiệp và dự án nghiên cứu.
  • Kỹ sư môi trường và chuyên gia tư vấn đánh giá tác động môi trường (ĐTM): Tham khảo phương pháp nội suy không gian IDW để mô phỏng phát tán chất ô nhiễm (Fe, Mn, BOD5, COD), phục vụ lập báo cáo hiện trạng môi trường và thiết kế mạng lưới quan trắc theo đúng quy chuẩn kỹ thuật.

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

Tại sao thuật toán nội suy IDW lại được ưu tiên lựa chọn hơn Kriging và Spline trong nghiên cứu này?
Phương pháp IDW được chọn vì tính toán trực tiếp dựa trên trung bình trọng số khoảng cách, giúp đơn giản hóa cấu trúc mô hình bản đồ mà vẫn đảm bảo độ chính xác cao khi xử lý tập dữ liệu 20 điểm quan trắc. Ngược lại, Spline dễ tạo ra các giá trị biến động cực đoan vượt quá ngưỡng thực tế tại các vùng chuyển tiếp, trong khi Kriging đòi hỏi phân tích bán phương sai phức tạp và phụ thuộc nhiều vào giả định phân phối chuẩn của số liệu.

Nghiên cứu đã sử dụng bộ dữ liệu quan trắc cụ thể như thế nào để thử nghiệm hệ thống?
Luận văn đã thu thập và số hóa dữ liệu tọa độ không gian cùng các chỉ số phân tích hóa lý từ 20 trạm quan trắc nước ngầm (ký hiệu từ G1 đến G20) trên địa bàn tỉnh Bắc Ninh. Dữ liệu được đo đạc qua 2 đợt trong năm 2014, bao gồm mực nước tĩnh, hàm lượng kim loại nặng (Fe, Mn) và các thông số chất lượng nước đối chiếu theo chuẩn QCVN 09:2008/BTNMT.

Thực trạng tuân thủ quy định bảo vệ tài nguyên nước của các doanh nghiệp tại Bắc Ninh được ghi nhận ra sao?
Báo cáo thanh tra trong luận văn chỉ ra rằng 100% đơn vị được thanh tra đều có hành vi chưa tuân thủ quy định pháp luật. Đáng chú ý, có tới 80% cơ sở (8/10 đơn vị) không thực hiện hoặc thực hiện không đầy đủ việc quan trắc, giám sát lưu lượng và chất lượng nước ngầm theo Quyết định 15/2008/QĐ-BTNMT, và 60% đơn vị không nộp báo cáo định kỳ.

Hệ thống GIS xây dựng trong luận văn cung cấp những tính năng nghiệp vụ chính nào?
Hệ thống cho phép thực hiện 4 nhóm chức năng cốt lõi: lưu trữ và số hóa dữ liệu không gian các trạm quan trắc; truy vấn thuộc tính chi tiết của từng giếng khoan; tự động chạy thuật toán nội suy IDW để tạo bản đồ phân vùng nồng độ ô nhiễm kim loại (Fe, Mn); và thực hiện phép chồng lớp bản đồ (Overlay) phục vụ đánh giá tổng hợp chất lượng môi trường nước.

Làm thế nào để ứng dụng mở rộng mô hình nghiên cứu này cho các địa phương khác?
Hệ thống được thiết kế theo cấu trúc mô đun hóa chuẩn của ArcGIS Desktop, do đó các địa phương khác hoàn toàn có thể tái sử dụng bằng cách: chuẩn hóa lại lớp bản đồ ranh giới hành chính, cập nhật bảng tọa độ các giếng quan trắc thực tế, và nhập các chỉ số phân tích môi trường mới vào cơ sở dữ liệu để tự động xuất bản đồ chuyên đề phân vùng ô nhiễm.

Kết luận

Luận văn đã giải quyết thành công bài toán tích hợp công nghệ Hệ thống Thông tin Địa lý vào công tác quản trị môi trường thông qua 5 kết quả cốt lõi:

  • Hệ thống hóa toàn diện cơ sở lý thuyết GIS, các mô hình dữ liệu không gian (Vector, Raster, TIN) và so sánh định lượng 3 thuật toán nội suy không gian (IDW, Spline, Kriging).
  • Phân tích hiện trạng chất lượng nước ngầm tại 20 trạm quan trắc tỉnh Bắc Ninh năm 2014, chỉ rõ các khu vực có hàm lượng Fe, Mn vượt ngưỡng QCVN 09:2008/BTNMT và thực trạng 80% cơ sở vi phạm nghĩa vụ giám sát môi trường.
  • Ứng dụng thành công thuật toán IDW trên ArcGIS Desktop để xây dựng các bản đồ chuyên đề phân vùng ô nhiễm nước ngầm có độ trực quan và chính xác cao.
  • Thiết lập giải pháp quản lý tích hợp cho phép tra cứu, hiển thị thuộc tính và chồng phủ các lớp dữ liệu không gian, giúp nâng cao hiệu suất làm việc của cán bộ quản lý môi trường.
  • Đề xuất hệ thống giải pháp đồng bộ về chính sách, hạ tầng xử lý nước thải tại 15 khu công nghiệp và đào tạo nhân lực kỹ thuật cho tỉnh Bắc Ninh giai đoạn đến năm 2030.

Nghiên cứu là tài liệu tham khảo khoa học và ứng dụng thực tiễn giá trị cao. Các cơ quan quản lý và nhà nghiên cứu quan tâm có thể khai thác mô hình dữ liệu này để tiếp tục phát triển các hệ thống giám sát môi trường WebGIS thời gian thực phục vụ phát triển kinh tế - xã hội bền vững.