Tổng quan về luận án
Tài nguyên nước mặt giữ vai trò quyết định đối với cân bằng sinh thái, an ninh lương thực và sự phát triển kinh tế – xã hội của mỗi quốc gia. Tại Việt Nam, tổng lượng dòng chảy mặt hàng năm đạt khoảng 830 – 840 tỷ m³/năm, tuy nhiên tài nguyên này phân bố không đồng đều theo không gian và thời gian; trong đó nguồn nước nội sinh chỉ chiếm 310 – 315 tỷ m³, còn lại 520 – 525 tỷ m³ (chiếm gần 63%) bắt nguồn từ lãnh thổ các quốc gia thượng nguồn. Lượng nước bình quân đầu người tại Việt Nam năm 2012 đạt khoảng 9.560 m³/người/năm, thấp hơn ngưỡng chuẩn 10.000 m³/người/năm đối với quốc gia có tài nguyên nước trung bình theo tiêu chí của Hiệp hội Nước Quốc tế (IWRA). Cùng với tốc độ đô thị hóa và công nghiệp hóa nhanh chóng, hệ thống nước mặt nội địa, đặc biệt tại các đô thị lớn như Thủ đô Hà Nội – nơi tập trung 22% số làng nghề cả nước (1.350 làng có nghề, 286 làng nghề truyền thống) với thải lượng nước thải chế biến lương thực lên tới 1.000 – 7.000 m³/ngày và lượng nước thải khu công nghiệp đạt 36.926 m³/ngày – đang đứng trước nguy cơ ô nhiễm hữu cơ và phú dưỡng nghiêm trọng.
┌──────────────────────────────────────────────┐
│ HỆ THỐNG VỆ TINH VNREDSAT-1A │
│ (B1: Blue, B2: Green, B3: Red, B4: NIR) │
└──────────────────────┬───────────────────────┘
│
Hiệu chỉnh bức xạ (Gain, Bias)
& Hiệu chỉnh khí quyển (DOS/FLAASH)
│
▼
┌─────────────────────────┐ ┌──────────────────┐ ┌─────────────────────────┐
│ QUANG TRẮC HIỆN TRƯỜNG │ │ PHẢN XẠ BỀ MẶT │ │ HÓA NGHIỆM TIÊU CHUẨN │
│ Máy GER 1500 (350-1050nm)├───► │ BOA: ρ(λ) │ ◄────┤ QCVN 08-MT:2015/BTNMT │
│ Tọa độ GPS trạm đo │ │ 4 kênh phổ (10m) │ │ TSS, Turbidity, COD, BOD│
└─────────────────────────┘ └────────┬─────────┘ └─────────────────────────┘
│
Phân tích tương quan &
Mô hình hóa hồi quy đa biến
│
▼
┌──────────────────────────────────────────────┐
│ BẢN ĐỒ CHẤT LƯỢNG NƯỚC MẶT HÀ NỘI │
│ Phân hóa: Sông Hồng & Hồ đầm nội thành │
└──────────────────────────────────────────────┘
Phương pháp quan trắc truyền thống dựa trên lấy mẫu điểm tại hiện trường và phân tích hóa nghiệm tuy mang lại độ tin cậy tại điểm đo nhưng bộc lộ hạn chế lớn về tính đại diện không gian, độ trễ thời gian và chi phí vận hành cao, không thể bao quát toàn diện sự biến động môi trường nước trên diện rộng. Công nghệ viễn thám quang học độ phân giải cao mở ra khả năng giám sát liên tục, đồng bộ và diện rộng. Luận án tiến sĩ kỹ thuật của tác giả Đinh Thị Thu Hiền (2018) với đề tài "Nghiên cứu xây dựng quy trình công nghệ giám sát môi trường nước mặt khu vực Hà Nội từ dữ liệu ảnh vệ tinh VNREDSat-1A" do PGS.TS. Trần Xuân Trường và TS. Lê Quốc Hưng hướng dẫn tại Trường Đại học Mỏ – Địa chất và Cục Viễn thám Quốc gia, là công trình tiên phong khai thác dữ liệu từ vệ tinh quan sát tài nguyên đầu tiên của Việt Nam (VNREDSat-1A, phóng lên quỹ đạo tháng 5/2013) nhằm giải quyết bài toán giám sát chất lượng nước mặt đô thị và sông ngòi phù sa nhiệt đới.
Khoảng trống nghiên cứu (research gap) trọng tâm mà luận án giải quyết bao gồm:
- Thiếu hụt quy trình công nghệ chuẩn hóa tích hợp đồng bộ từ tiền xử lý quang học bức xạ, hiệu chỉnh khí quyển đến mô hình hóa bán thực nghiệm phù hợp với đặc tính dải phổ và độ phân giải không gian 10m của bộ cảm đa phổ VNREDSat-1A.
- Thiếu các mô hình hồi quy đa kênh chuyên biệt phân tách rõ rệt giữa hai thủy vực có bản chất quang học tương phản: hệ thống sông phù sa tải lượng khoáng cao (sông Hồng) và hệ thống hồ đầm đô thị chịu hiện tượng phú dưỡng và ô nhiễm hữu cơ nhân sinh.
- Hạn chế trong việc chuyển đổi và mở rộng khả năng viễn thám từ các thông số hoạt tính quang học (optically active parameters) như tổng chất rắn lơ lửng (TSS), độ đục sang các thông số chất lượng nước phi hoạt tính quang học (non-optically active parameters) như nhu cầu oxy hóa học (COD) và nhu cầu oxy sinh hóa (BOD5).
Hệ thống câu hỏi nghiên cứu và giả thuyết khoa học được thiết lập chặt chẽ:
- RQ1: Dữ liệu phổ phản xạ bề mặt thu nhận từ bộ cảm đa phổ VNREDSat-1A sau khi hiệu chỉnh khí quyển có đáp ứng độ chính xác tương thích với phổ đo bức xạ hiện trường hay không?
- RQ2: Sự kết hợp giữa các kênh phổ đơn lẻ, tỷ số kênh hay tổ hợp đa kênh tuyến tính mang lại độ chính xác tối ưu trong ước lượng nồng độ TSS, độ đục, COD và BOD5?
- RQ3: Có tồn tại sự khác biệt mang tính quy luật về phương trình tương quan phổ – chất lượng nước giữa thủy vực hở dòng chảy mạnh (sông Hồng) và thủy vực kín nước đứng (hồ đô thị Hà Nội)?
- H1 (Giả thuyết 1): Tồn tại mối tương quan thống kê cao và có ý nghĩa ($R^2 > 0.80$) giữa phổ phản xạ bề mặt từ ảnh vệ tinh VNREDSat-1A đã hiệu chỉnh khí quyển với hàm lượng các thông số chất lượng nước mặt.
- H2 (Giả thuyết 2): Quy trình công nghệ tích hợp đo phổ hiện trường GER 1500, hiệu chỉnh hình học – khí quyển chính xác và mô hình hồi quy đa biến tối ưu sẽ cho phép thiết lập bản đồ phân bố định lượng chất lượng nước mặt đạt độ tin cậy tương đương số liệu quan trắc tiêu chuẩn theo QCVN 08-MT:2015/BTNMT.
Khung lý thuyết của công trình dựa trên Lý thuyết truyền xạ quang học thủy vực (Radiative Transfer Theory - Gordon, 1988), Mô hình quang học sinh học (Bio-optical Model - Morel & Prieur, 1977) phân loại nước Case-1 và Case-2, cùng Lý thuyết hiệu chỉnh bức xạ khí quyển Dark Object Subtraction (DOS - Chavez, 1988). Phạm vi nghiên cứu thực nghiệm được tiến hành trên địa bàn Hà Nội qua 02 cảnh ảnh VNREDSat-1A chụp vào ngày 20/10/2016 (mùa kiệt chuyển tiếp) và ngày 21/12/2017 (chính mùa đông khô), đồng bộ với hai đợt lấy mẫu hóa nghiệm và đo phổ phản xạ mặt nước bằng máy quang phổ kế GER 1500 trên hệ thống sông Hồng và các hồ đại diện (Hồ Tây, Hoàn Kiếm, Yên Sở, Linh Đàm).
Literature Review và Positioning
Nghiên cứu ứng dụng viễn thám trong giám sát chất lượng nước mặt đã trải qua các giai đoạn phát triển rõ rệt trên phạm vi quốc tế và trong nước. Dòng nghiên cứu nền tảng khởi đầu từ việc khám phá mối tương quan tuyến tính giữa phổ phản xạ kênh đơn ở dải sóng nhìn thấy và cận hồng ngoại (0.7 – 0.8 µm) với hàm lượng chất rắn lơ lửng (Ritchie et al., 1976; Novo et al., 1991). Khi nồng độ hạt lơ lửng tăng, năng lượng tán xạ ngược của khối nước trong dải sóng đỏ và cận hồng ngoại tăng tuyến tính trước khi đạt mức bão hòa quang học.
Dòng nghiên cứu thứ hai tập trung vào phát triển các chỉ số tỷ số kênh (band ratios) nhằm triệt tiêu hiệu ứng phông nhiễu từ khí quyển và góc chiếu mặt trời. Doxaran et al. (2002, 2006, 2009) khi nghiên cứu vùng cửa sông Gironde (Pháp) trên dữ liệu ảnh SPOT HRV và Landsat ETM+ đã chứng minh tỷ số giữa kênh cận hồng ngoại (NIR) và kênh xanh lục (Green), hoặc giữa kênh đỏ (Red) và xanh lục, cho phép ước tính hàm lượng chất lơ lửng (SPM) với hệ số tương quan vượt trội, khắc phục triệt để hiện tượng bất đồng nhất của dòng chảy phù sa. Ouillon et al. (2008) phân tích dữ liệu đa cảm biến (MODIS, MERIS, TM, SeaWiFS) tại các thủy vực ven biển Cuba, Fiji và New Caledonia, xác lập tỷ số phản xạ viễn thám $R_{rs}(510)/R_{rs}(681)$ và bước sóng 0.620 µm là tối ưu cho việc xây dựng bản đồ độ đục (Turbidity).
Dòng nghiên cứu thứ ba khai thác ảnh vệ tinh độ phân giải không gian cao (IKONOS, QuickBird, WorldView) kết hợp phương pháp hồi quy đa biến (MLR) và trí tuệ nhân tạo để giám sát đồng thời các thông số dinh dưỡng và hữu cơ. Liu et al. (2015) ứng dụng ảnh IKONOS và mô hình mạng nơ-ron tại hồ Cihu (Trung Quốc), xác định hàm lượng tổng nitơ (TN) và tổng photpho (TP) đạt hệ số $R^2 > 0.86$ và sai số $RMSE < 0.89$. Ustun (2011) sử dụng 4 kênh phổ IKONOS để mô hình hóa hàm lượng chất lơ lửng và nitơ tại hồ Buyukcekmece (Thổ Nhĩ Kỳ) với $R^2 > 0.90$. Somvanshi et al. (2011) phát triển chỉ số độ đục chuẩn hóa NDTI (Normalized Difference Turbidity Index) từ ảnh QuickBird trên sông Gomti Lucknow (Ấn Độ). Weipi He et al. (2008) ứng dụng ảnh Landsat TM và 76 điểm quan trắc tại các hồ Bắc Kinh để tính toán chlorophyll-a, COD, TN và $NH_3\text{-}N$.
┌─────────────────────────────────────────────────────────────────────────────┐
│ TIẾN TRÌNH LITERATURE REVIEW │
├─────────────────────────────────────────────────────────────────────────────┤
│ 1. Kênh đơn tuyến tính: Ritchie et al. (1976), Novo et al. (1991) │
│ └─► Ước tính TSS nồng độ thấp đến trung bình (0.7 - 0.8 µm) │
│ │
│ 2. Tỷ số kênh (Band Ratios): Doxaran et al. (2002), Ouillon et al. (2008) │
│ └─► Triệt tiêu hiệu ứng phông khí quyển (NIR/Green, Red/Green, NDTI) │
│ │
│ 3. Đa biến độ phân giải cao: Liu et al. (2015), Ustun (2011), He (2008) │
│ └─► Ước tính gián tiếp COD, BOD5, TN, TP qua tương quan đa phổ │
│ │
│ 4. Định vị của Luận án (Đinh Thị Thu Hiền, 2018): │
│ └─► Chuẩn hóa trên VNREDSat-1A (10m), giải mã Case-2 thủy vực nhiệt đới │
└─────────────────────────────────────────────────────────────────────────────┘
Trong bức tranh học thuật đó, tồn tại hai luồng tranh luận học thuật lớn:
- Tranh luận về mô hình kênh đơn so với mô hình tỷ số / tổ hợp đa kênh: Trường phái truyền thống (Ritchie, Novo) lập luận rằng kênh phổ đơn dải đỏ/cận hồng ngoại phản ánh trực tiếp vật lý tán xạ hạt lơ lửng. Ngược lại, trường phái quang học hiện đại (Doxaran, Ouillon) khẳng định môi trường nước lục địa có tính chất quang học phức tạp (Case-2 Waters), sự biến thiên của chất hữu cơ hòa tan màu (CDOM) và diệp lục sẽ gây nhiễu tín hiệu kênh đơn, do đó bắt buộc phải sử dụng tỷ số kênh hoặc hồi quy đa biến đa kênh để hiệu chỉnh chéo.
- Tranh luận về khả năng viễn thám các thông số phi hoạt tính quang học (BOD5, COD): Các nhà quang học thuần túy (Mobley, Morel) giữ quan điểm hoài nghi vì các phân tử hữu cơ hòa tan không có cấu trúc hấp thụ/tán xạ trực tiếp đặc trưng trong dải phổ khả kiến. Ngược lại, các nhà nghiên cứu môi trường thực nghiệm (Weipi He et al., 2008; Oki et al., 2012) chứng minh rằng trong điều kiện ô nhiễm xả thải đô thị tập trung, các chất ô nhiễm hữu cơ và vi sinh gắn kết chặt chẽ với các hạt lơ lửng và sinh khối tảo, cho phép viễn thám gián tiếp thông qua tương quan đồng biến (co-variation).
Tại Việt Nam, các nghiên cứu trước đây (Lê Trình, 2008; Tôn Thất Lãng, 2010; Nguyễn Duy Phú, 2012; Phạm Thế Anh & Nguyễn Văn Huy, 2013) chủ yếu áp dụng chỉ số chất lượng nước WQI trên số liệu quan trắc trạm đo hoặc sử dụng ảnh vệ tinh MODIS, Landsat có độ phân giải không gian từ 30m đến 1000m, vốn không phù hợp để quan sát các đoạn sông hẹp và hồ nhỏ trong đô thị. Luận án của Đinh Thị Thu Hiền định vị chính xác vào khoảng trống này: thiết lập một hệ phương pháp chuẩn mực khai thác dữ liệu ảnh độ phân giải cao 10m từ vệ tinh VNREDSat-1A của Việt Nam, xây dựng các thuật toán hồi quy đa kênh thích ứng với tính chất quang học phức tạp của sông hồ nhiệt đới Hà Nội, đặt nền móng độc lập tự chủ công nghệ viễn thám quốc gia.
Đóng góp lý thuyết và khung phân tích
Đóng góp cho lý thuyết
Luận án mở rộng và cụ thể hóa Lý thuyết truyền xạ quang học thủy vực của Gordon (1988) và Mô hình quang học sinh học của Morel & Prieur (1977) đối với loại nước Case-2 chịu tác động nhân sinh cực đoan tại vùng khí hậu nhiệt đới gió mùa. Công trình chỉ ra giới hạn của Định luật suy giảm Beer-Lambert-Bouguer khi áp dụng trực tiếp cho nước mặt đô thị có sự pha trộn phức tạp giữa hạt lơ lửng khoáng vô cơ (phù sa sông Hồng) và tải lượng hữu cơ phân hủy sinh học cao từ các làng nghề và khu công nghiệp.
┌────────────────────────────────────────────────────────┐
│ KHUNG LÝ THUYẾT TÍCH HỢP TỔNG THỂ │
└───────────────────────────┬────────────────────────────┘
│
┌────────────────────────────────────┼────────────────────────────────────┐
│ │ │
▼ ▼ ▼
┌───────────────────────┐ ┌───────────────────────┐ ┌───────────────────────┐
│ Lý thuyết Truyền xạ │ │ Mô hình Quang sinh học│ │ Lý thuyết Hiệu chỉnh │
│ Gordon (1988) │ │ Morel & Prieur (1977) │ │ Khí quyển (Chavez) │
│ Tán xạ ngược / Hấp thụ│ │ Nước Case-2 phức hợp │ │ Loại bỏ tán xạ Rayleigh│
└───────────┬───────────┘ └───────────┬───────────┘ └───────────┬───────────┘
│ │ │
└────────────────────────────────────┼────────────────────────────────────┘
│
▼
┌────────────────────────────────────────────────────────┐
│ MÔ HÌNH HỒI QUY ĐA KÊNH VNREDSAT-1A │
│ Sông Hồng: f(B1, B2, B3, B4) -> TSS, Độ đục │
│ Hồ đầm: f(B1, B2, B3, B4) -> TSS, COD, BOD5 │
└────────────────────────────────────────────────────────┘
Mô hình lý thuyết của luận án thiết lập ba mệnh đề khoa học (Propositions):
- Mệnh đề 1 ($P_1$): Tín hiệu bức xạ phản xạ rời mặt nước ($R_{rs}$) trong các dải phổ xanh lam (B1: 0.45 – 0.52 µm), xanh lục (B2: 0.53 – 0.59 µm), đỏ (B3: 0.63 – 0.69 µm) và cận hồng ngoại (B4: 0.76 – 0.89 µm) của VNREDSat-1A phản ánh sự đóng góp tổng hợp của các đặc tính quang học nội tại (IOPs - Inherent Optical Properties: hệ số hấp thụ $a(\lambda)$ và hệ số tán xạ ngược $b_b(\lambda)$).
- Mệnh đề 2 ($P_2$): Sự kết hợp tương hỗ giữa phổ phản xạ dải đỏ (B3 - nhạy cảm với tán xạ hạt lơ lửng) và dải cận hồng ngoại (B4 - hấp thụ mạnh bởi phân tử nước thuần khiết, nhạy cảm với nồng độ hạt đậm đặc) cho phép cô lập và lượng hóa chính xác nồng độ TSS và độ đục trong môi trường nước phù sa sông Hồng.
- Mệnh đề 3 ($P_3$): Các thông số ô nhiễm hữu cơ (COD, BOD5) thể hiện tính chất đồng biến quang học với phổ phản xạ đa kênh thông qua mối liên kết vật lý – sinh hóa giữa hàm lượng chất rắn lơ lửng hữu cơ, sinh khối vi sinh vật và các hợp chất tiêu hao oxy hòa tan.
Công trình đánh dấu bước chuyển dịch mô hình (paradigm shift) từ phương pháp tiếp cận nội suy không gian dựa trên các điểm đo tĩnh rời rạc (point-based deterministic interpolation) sang mô hình trường phân bố động lực học quang học viễn thám liên tục (continuous bio-optical remote sensing field).
Khung phân tích độc đáo
Khung phân tích của luận án tích hợp bốn cấu phần lý thuyết và kỹ thuật:
- Lý thuyết truyền xạ quang học thủy vực (Gordon, 1988; Mobley, 1994).
- Mô hình tính chất quang học biểu kiến và nội tại (Morel & Prieur, 1977).
- Lý thuyết hiệu chỉnh bức xạ khí quyển vật lý thực nghiệm (Dark Object Subtraction - Chavez, 1988).
- Phương pháp phân tích hồi quy toán học đa biến từng bước (Stepwise Multivariate Linear Regression).
Các khái niệm then chốt được định nghĩa tường minh:
- Phổ phản xạ viễn thám ($R_{rs}(\lambda)$): Tỷ số giữa độ chói bức xạ rời mặt nước $L_w(\lambda)$ và độ rọi bức xạ chiếu xuống $E_d(\lambda)$ tại bước sóng $\lambda$ ($sr^{-1}$).
- Phản xạ bề mặt ($\rho(\lambda)$): Giá trị phản xạ đáy khí quyển (BOA - Bottom-of-Atmosphere) đã loại bỏ ảnh hưởng tán xạ phân tử Rayleigh và tán xạ sol khí (Aerosol scattering).
- Độ đục và TSS: Đại lượng đại diện cho mức độ suy giảm độ truyền quang và tổng khối lượng chất lơ lửng dạng hạt tính trên một đơn vị thể tích nước ($mg/l$ hoặc NTU).
Điều kiện biên (boundary conditions) của mô hình được xác lập: áp dụng tối ưu cho các thủy vực nước ngọt nội địa vùng nhiệt đới, độ sâu tầng nước $> 1.5m$ để triệt tiêu hiệu ứng phản xạ đáy sông hồ; điều kiện khí quyển quang đãng với độ che phủ mây $< 10%$; góc thiên đỉnh mặt trời tại thời điểm chụp đảm bảo không xuất hiện hiện tượng phản xạ gương (sun glint).
Phương pháp nghiên cứu tiên tiến
Thiết kế nghiên cứu
Nghiên cứu được thiết kế theo trường phái thực chứng luận (Positivism), kết hợp chặt chẽ giữa đo đạc thực địa hiện trường đồng bộ, phân tích hóa nghiệm tiêu chuẩn phòng thí nghiệm và xử lý ảnh số viễn thám đa phổ độ phân giải cao. Thiết kế nghiên cứu đa tầng bậc (Multi-level design) được triển khai trên ba cấp độ:
- Cấp độ vi mô: Phổ bức xạ siêu phổ mặt nước đo bằng máy GER 1500 (dải phổ 350 – 1050 nm, 512 kênh phổ).
- Cấp độ điểm trắc vi mô: Các trạm lấy mẫu chất lượng nước mặt phân tích trong phòng thí nghiệm.
- Cấp độ vĩ mô: Không gian toàn cảnh đô thị Hà Nội từ ảnh vệ tinh VNREDSat-1A (kích thước khung chụp 17.5 km × 17.5 km, độ phân giải không gian 10m ở chế độ đa phổ và 2.5m ở chế độ toàn sắc).
┌─────────────────────────────────────────────────────────────────────────────┐
│ QUY TRÌNH NGHIÊN CỨU RIGOROUS │
├─────────────────────────────────────────────────────────────────────────────┤
│ 1. DỮ LIỆU ĐẦU VÀO: │
│ ├─ Ảnh VNREDSat-1A (20/10/2016 & 21/12/2017) │
│ ├─ Phổ GER 1500 (512 kênh: 350 - 1050 nm) │
│ └─ Mẫu nước thực địa (TSS, Turbidity, COD, BOD5) │
│ │
│ 2. TIỀN XỬ LÝ ẢNH CHUYÊN SÂU: │
│ ├─ Hiệu chuẩn bức xạ: L(λ) = Gain × DN + Bias │
│ ├─ Hiệu chỉnh khí quyển: DOS / FLAASH -> Phản xạ bề mặt ρ(λ) │
│ └─ Hiệu chỉnh hình học: Nắn chỉnh trực giao DEM + GCPs (RMSE < 0.5 pixel)│
│ │
│ 3. MÔ HÌNH HÓA THỰC NGHIỆM: │
│ ├─ Hồi quy đơn kênh: f(Bi) │
│ ├─ Hồi quy tỷ số kênh: f(Bi/Bj) │
│ ├─ Hồi quy tổ hợp 3 kênh: f(Bi, Bj, Bk) │
│ └─ Hồi quy toàn kênh 4 phổ: f(B1, B2, B3, B4) -> Tối ưu R2, RMSE │
│ │
│ 4. THÀNH LẬP BẢN ĐỒ CHẤT LƯỢNG NƯỚC (GIS): │
│ └─ Phân hạng theo QCVN 08-MT:2015/BTNMT (A1, A2, B1, B2) │
└─────────────────────────────────────────────────────────────────────────────┘
Quy trình nghiên cứu rigorous
Quy trình thực nghiệm được chuẩn hóa nghiêm ngặt qua 4 bước:
- Lấy mẫu và đo đạc thực địa: Các đợt khảo sát hiện trường được tổ chức trùng khớp với thời gian vệ tinh VNREDSat-1A bay qua chụp ảnh khu vực Hà Nội. Mẫu nước mặt được thu thập tại tầng mặt ($0 - 0.5m$) theo tiêu chuẩn TCVN, bảo quản lạnh ở nhiệt độ $4^\circ C$ và chuyển ngay về phòng thí nghiệm trong vòng 24 giờ để xác định các chỉ tiêu: Độ đục (máy đo độ đục), TSS (phương pháp sấy khô ở $105^\circ C$), COD (phương pháp chuẩn độ Bicromat), BOD5 (phương pháp ủ 5 ngày ở $20^\circ C$).
- Đo phổ phản xạ mặt nước: Sử dụng máy quang phổ kế bức xạ GER 1500, đo đạc phổ phản xạ của khối nước và tấm chuẩn trắng phản xạ tuyệt đối (Spectralon panel) để tính toán phổ phản xạ thực địa chuẩn xác.
- Tiền xử lý ảnh vệ tinh:
- Hiệu chuẩn bức xạ: Chuyển đổi giá trị số số học (Digital Number - DN) sang giá trị bức xạ phổ tại đỉnh khí quyển (TOA Radiance) thông qua các hệ số Gain và Bias do Viện Hàn lâm Khoa học và Công nghệ Việt Nam (VAST) cung cấp.
- Hiệu chỉnh khí quyển: Áp dụng phương pháp vật thể tối DOS và mô hình truyền xạ nhằm loại bỏ hoàn toàn tán xạ Rayleigh và hấp thụ của sol khí, chuyển đổi TOA Radiance thành phổ phản xạ bề mặt (Surface Reflectance).
- Hiệu chỉnh hình học: Sử dụng các điểm khống chế mặt đất (GCPs) từ bản đồ địa hình tỷ lệ lớn kết hợp mô hình số độ cao (DEM) để nắn chỉnh trực giao, đạt sai số vị trí hình học $RMSE < 0.5$ pixel (tương đương sai số mặt đất $< 5m$).
- Tam giác đạc (Triangulation): Kiểm định chéo độ chính xác giữa phổ phản xạ ảnh vệ tinh sau hiệu chỉnh khí quyển với phổ đo hiện trường GER 1500, đảm bảo độ tương đồng phổ trước khi đưa vào mô hình hồi quy.
Data và phân tích
Cơ sở dữ liệu thực nghiệm bao gồm hai cảnh ảnh VNREDSat-1A (20/10/2016 và 21/12/2017) với 4 kênh phổ:
- Kênh 1 (Blue): $0.45 - 0.52,\mu m$
- Kênh 2 (Green): $0.53 - 0.59,\mu m$
- Kênh 3 (Red): $0.63 - 0.69,\mu m$
- Kênh 4 (NIR): $0.76 - 0.89,\mu m$
Phương pháp phân tích hồi quy toán học được triển khai toàn diện qua 4 cấp độ cấu hình biến:
- Mô hình hồi quy đơn kênh: $Y = f(B_i)$ với $i \in {1, 2, 3, 4}$.
- Mô hình hồi quy tỷ số 2 kênh: $Y = f(B_i / B_j)$ (các tổ hợp $B1/B2, B1/B3, B1/B4, B2/B3, B2/B4, B3/B4$).
- Mô hình hồi quy tổ hợp 3 kênh: $Y = f(B_i, B_j, B_k)$ (các tổ hợp $B1\text{-}B2\text{-}B3, B1\text{-}B2\text{-}B4, B1\text{-}B3\text{-}B4, B2\text{-}B3\text{-}B4$).
- Mô hình hồi quy toàn kênh 4 phổ: $Y = f(B1, B2, B3, B4)$.
Phần mềm chuyên dụng được sử dụng bao gồm ENVI (xử lý ảnh số quang học), ArcGIS (phân tích không gian và biên tập bản đồ chuyên đề), và SPSS/R (phân tích hồi quy đa biến và kiểm định thống kê). Các chỉ số đánh giá độ chính xác bao gồm: Hệ số tương quan Pearson ($R$), Hệ số xác định ($R^2$), Sai số trung phương ($RMSE$) và Sai số tương đối trung bình ($MRE%$).
Phát hiện đột phá và implications
Những phát hiện then chốt
| Thông số |
Thủy vực |
Cấu hình mô hình tối ưu |
Kênh phổ sử dụng |
Hệ số $R^2$ |
Đánh giá sai số (RMSE) |
| Độ đục |
Sông Hồng |
Hồi quy đa biến 4 kênh |
B1, B2, B3, B4 |
$0.88 - 0.93$ |
Rất thấp ($p < 0.001$) |
| TSS |
Sông Hồng |
Hồi quy đa biến 4 kênh / Tỷ số |
B3, B4 / B1, B2, B3, B4 |
$0.86 - 0.91$ |
Sai số $< 10%$ |
| Độ đục |
Hồ, đầm nội thành |
Hồi quy đa biến 4 kênh |
B1, B2, B3, B4 |
$0.84 - 0.89$ |
Sai số $< 12%$ |
| TSS |
Hồ, đầm nội thành |
Hồi quy đa biến 4 kênh |
B1, B2, B3, B4 |
$0.82 - 0.87$ |
Sai số $< 12%$ |
| COD |
Hồ, đầm nội thành |
Hồi quy tổ hợp đa kênh |
B1, B2, B3, B4 |
$0.78 - 0.83$ |
Đạt độ tin cậy ($p < 0.01$) |
| BOD5 |
Hồ, đầm nội thành |
Hồi quy tổ hợp đa kênh |
B1, B2, B3, B4 |
$0.75 - 0.81$ |
Đạt độ tin cậy ($p < 0.01$) |
- Hiệu năng vượt trội của mô hình tổ hợp toàn kênh (4 kênh phổ): Kết quả thực nghiệm khẳng định việc sử dụng đồng thời 4 kênh phổ của VNREDSat-1A cho hệ số xác định $R^2$ cao nhất ($R^2$ đạt $0.82 - 0.93$) và sai số RMSE nhỏ nhất đối với cả độ đục và TSS, vượt trội rõ rệt so với mô hình hồi quy đơn kênh ($R^2$ dao động từ $0.45 - 0.68$) hoặc mô hình tỷ số hai kênh ($R^2$ dao động từ $0.60 - 0.76$).
- Vai trò quyết định của hiệu chỉnh khí quyển: So sánh đối chứng phương trình hồi quy trước và sau hiệu chỉnh khí quyển cho thấy, việc loại bỏ bức xạ tán xạ khí quyển giúp cải thiện hệ số xác định $R^2$ thêm 15% – 30%, biến đổi phổ TOA sai lệch thành phổ phản xạ bề mặt tương thích hoàn toàn với phổ đo thực địa GER 1500.
- Phân hóa quang phổ sâu sắc giữa sông Hồng và hồ đầm đô thị: Luận án phát hiện quy luật phản xạ phổ phân kỳ rõ rệt giữa hai môi trường thủy vực. Sông Hồng có đỉnh phản xạ dịch chuyển mạnh về dải đỏ và cận hồng ngoại do nồng độ hạt phù sa vô cơ cao; trong khi các hồ nội thành có đỉnh phản xạ ở dải xanh lục và hấp thụ mạnh ở dải đỏ/cận hồng ngoại do nồng độ chlorophyll và chất hữu cơ cao. Do đó, việc áp dụng chung một phương trình hồi quy là phi khoa học; luận án đã xây dựng hai bộ phương trình đặc thù riêng biệt cho từng loại thủy vực.
- Khả năng ước lượng gián tiếp các thông số ô nhiễm hữu cơ (COD, BOD5): Dù COD và BOD5 là các thông số phi hoạt tính quang học, mô hình tổ hợp 4 kênh phổ trên các hồ đầm Hà Nội vẫn đạt $R^2 > 0.75$, chứng minh sự tồn tại của mối liên kết sinh địa hóa gián tiếp giữa chất hữu cơ hòa tan/lơ lửng với tín hiệu tán xạ đa phổ.
- Thực trạng ô nhiễm nước mặt Hà Nội được định lượng chi tiết: Bản đồ chất lượng nước thành lập từ ảnh năm 2016 và 2017 phản ánh sự gia tăng độ đục và TSS tại các đoạn sông Hồng chảy qua khu vực ngoại thành tiếp nhận nước thải làng nghề (Hoài Đức, Đan Phượng) và sự suy giảm chất lượng nước nghiêm trọng tại các hồ đô thị phía Nam (Yên Sở, Linh Đàm) với nồng độ COD và BOD5 thường xuyên vượt giới hạn B1, chạm ngưỡng B2 theo QCVN 08-MT:2015/BTNMT.
Implications đa chiều
- Ý nghĩa lý thuyết: Xác lập cơ sở khoa học quang học thủy văn cho dữ liệu viễn thám độ phân giải cao của Việt Nam; hoàn thiện khung lý thuyết đánh giá nước Case-2 phức hợp tại vùng nhiệt đới châu Á.
- Đổi mới phương pháp luận: Chuẩn hóa quy trình khép kín từ tiền xử lý quang học đến mô hình hồi quy đa phổ thích ứng, có khả năng chuyển giao áp dụng cho các hệ thống vệ tinh quan sát trái đất hiện tại và tương lai của Việt Nam.
- Ứng dụng thực tiễn: Cung cấp công cụ công nghệ cho phép các cơ quan quản lý môi trường thành lập bản đồ chất lượng nước diện rộng định kỳ với chi phí chỉ bằng một phần nhỏ so với quan trắc truyền thống, phát hiện sớm các điểm nóng xả thải ô nhiễm.
- Khuyến nghị chính sách: Đề xuất Sở Tài nguyên và Môi trường Hà Nội cùng Bộ Tài nguyên và Môi trường tích hợp quy trình viễn thám VNREDSat-1A vào Hệ thống Quan trắc Môi trường Tự động Quốc gia; xây dựng cơ chế chia sẻ dữ liệu viễn thám định kỳ phục vụ kiểm soát xả thải làng nghề và quy hoạch không gian mặt nước Thủ đô.
Limitations và Future Research
Công trình nghiên cứu thừa nhận các giới hạn khoa học cụ thể:
- Giới hạn điều kiện thời tiết quang học: Vệ tinh quang học VNREDSat-1A hoàn toàn phụ thuộc vào độ che phủ mây và điều kiện chiếu sáng mặt trời. Tại miền Bắc Việt Nam, mùa mưa ẩm có mật độ mây cao làm giảm tần suất thu nhận ảnh sử dụng được.
- Độ phân giải phổ của bộ cảm đa phổ: Cảm biến VNREDSat-1A chỉ có 4 kênh phổ băng rộng, chưa thể phân tách chi tiết các phân lớp diệp lục ($Chl\text{-}a$, $Chl\text{-}b$, Phycocyanin) và chất hữu cơ hòa tan màu vàng (CDOM) như các bộ cảm siêu phổ (Hyperspectral) với hàng trăm kênh.
- Độ trễ dòng chảy sông ngòi: Tốc độ dòng chảy sông Hồng biến động nhanh theo chế độ xả lũ của các hồ chứa thượng nguồn (Hòa Bình, Tuyên Quang, Sơn La), tạo ra sự chênh lệch nhỏ về động lực học thủy văn giữa thời điểm lấy mẫu thực địa và thời điểm vệ tinh quét qua.
Chương trình nghiên cứu tương lai (Future Research Agenda):
- Phát triển thuật toán học máy (Machine Learning / Deep Learning như Random Forest, XGBoost, CNN) trên chuỗi thời gian ảnh viễn thám nhiều năm để nâng cao độ chính xác dự báo.
- Tích hợp dữ liệu viễn thám siêu phổ (Hyperspectral) từ các cảm biến hàng không (UAV/Máy bay) hoặc vệ tinh thế hệ mới để bóc tách riêng biệt từng hợp chất hữu cơ độc hại.
- Kết hợp công nghệ Radar độ mở tổng hợp (SAR) và viễn thám nhiệt (TIR) để khắc phục nhược điểm mây mù và giám sát ô nhiễm nhiệt độ bề mặt nước.
- Mở rộng quy trình công nghệ cho toàn bộ lưu vực sông Hồng – Thái Bình và lưu vực sông Đồng Nai – Cửu Long.
Tác động và ảnh hưởng
- Tác động học thuật: Khẳng định giá trị khoa học của nguồn dữ liệu viễn thám bản địa VNREDSat-1A trên trường quốc tế; cung cấp hệ số thực nghiệm và phương pháp luận làm tài liệu tham khảo chuẩn cho các nghiên cứu đo ảnh, viễn thám và quản lý tài nguyên nước tại các trường đại học, viện nghiên cứu chuyên ngành.
- Chuyển đổi ngành công nghiệp trắc địa - môi trường: Cung cấp quy trình công nghệ sẵn sàng chuyển giao cho các trung tâm viễn thám, viện quy hoạch môi trường và các doanh nghiệp tư vấn công nghệ địa không gian (Geospatial).
- Tác động chính sách công: Cung cấp cơ sở dữ liệu không gian trực quan, khách quan phục vụ việc thanh tra, giám sát môi trường và đánh giá hiệu quả thực thi Luật Bảo vệ Môi trường và các quy chuẩn kỹ thuật quốc gia (QCVN).
- Lợi ích xã hội và môi trường: Hỗ trợ bảo vệ nguồn nước sinh hoạt cho hơn 8 triệu dân cư Thủ đô Hà Nội; cảnh báo sớm nguy cơ suy thoái hệ sinh thái thủy sinh và hỗ trợ phục hồi các hồ đô thị bị ô nhiễm.
Đối tượng hưởng lợi
- Nghiên cứu sinh và học viên sau đại học: Tiếp cận quy trình nghiên cứu viễn thám môi trường nước chuẩn mực, từ đo phổ thực địa GER 1500, hiệu chỉnh khí quyển DOS đến mô hình hóa hồi quy đa biến.
- Các nhà khoa học viễn thám và môi trường: Kế thừa hệ thống phương trình tương quan phổ và phương pháp luận phân hóa thủy vực sông ngòi - hồ đầm trong môi trường nhiệt đới gió mùa.
- Kỹ sư và chuyên gia R&D môi trường: Ứng dụng quy trình công nghệ tự động hóa xử lý ảnh viễn thám để xây dựng hệ thống phần mềm GIS giám sát môi trường nước theo thời gian thực.
- Các nhà quản lý chính sách và tài nguyên nước: Sử dụng các bản đồ hiện trạng chất lượng nước để đưa ra quyết định quản lý lưu vực, phê duyệt quy hoạch đô thị và kiểm soát các nguồn xả thải công nghiệp, làng nghề.
Câu hỏi chuyên sâu
1. Đóng góp lý thuyết độc đáo nhất của luận án là gì?
Luận án đã mở rộng Mô hình quang học thủy vực Case-2 của Morel & Prieur (1977) và Lý thuyết truyền xạ Gordon (1988) bằng việc chứng minh và lượng hóa hiện tượng phân hóa quang phổ tương phản giữa thủy vực phù sa dòng chảy xiết (sông Hồng) và thủy vực phú dưỡng nước đứng (hồ đầm Hà Nội), từ đó chứng minh tính khả thi của việc viễn thám gián tiếp các thông số ô nhiễm hữu cơ phi hoạt tính quang học (COD, BOD5) qua tổ hợp hồi quy 4 kênh đa phổ VNREDSat-1A.
2. Điểm mới về phương pháp luận so với các nghiên cứu quốc tế là gì?
So với nghiên cứu của Doxaran et al. (2002, 2006) tại cửa sông Gironde (chỉ dùng tỷ số 2 kênh trên ảnh SPOT/Landsat) và Liu et al. (2015) tại hồ Cihu (chỉ nghiên cứu hồ kín trên ảnh IKONOS), luận án đã phát triển một quy trình tổng thể đa tầng bậc: tích hợp hiệu chuẩn bức xạ tuyệt đối, hiệu chỉnh khí quyển DOS/FLAASH, đồng bộ phổ thực địa GER 1500 và xây dựng hai hệ thống mô hình hồi quy toàn kênh 4 phổ ($B1\text{-}B2\text{-}B3\text{-}B4$) tối ưu hóa riêng biệt cho cả hai loại hình thủy vực sông ngòi và hồ đầm đô thị nhiệt đới.
3. Phát hiện bất ngờ nhất từ số liệu thực nghiệm là gì?
Phát hiện bất ngờ nhất là việc các mô hình hồi quy đa kênh tuyến tính kết hợp đầy đủ 4 kênh phổ ($B1, B2, B3, B4$) mang lại độ chính xác ước lượng TSS và độ đục vượt trội hoàn toàn ($R^2 > 0.85 - 0.93$) so với các mô hình tỷ số kênh truyền thống (như $B4/B2$ hay $B3/B2$ vốn chỉ đạt $R^2 \approx 0.60 - 0.76$). Điều này chứng minh rằng việc bổ sung kênh Blue (B1) và Green (B2) giúp hiệu chỉnh khử nhiễu nền hữu cơ và chất màu hòa tan hiệu quả hơn trong môi trường nước phù sa nhiệt đới.
4. Luận án có cung cấp quy trình tái lặp (Replication Protocol) không?
Có. Luận án cung cấp chi tiết từng bước quy trình công nghệ khép kín bao gồm: thông số Gain/Bias cho từng kênh phổ VNREDSat-1A, thuật toán hiệu chỉnh khí quyển vật thể tối, quy trình hiệu chuẩn tấm trắng Spectralon và vận hành máy GER 1500 ngoài hiện trường, phương pháp lấy mẫu bảo quản theo TCVN và bảng cấu trúc các phương trình hồi quy toán học cụ thể.
5. Định hướng chương trình nghiên cứu 10 năm tiếp theo được phác thảo như thế nào?
Định hướng 10 năm tập trung vào việc tự động hóa toàn diện quy trình giám sát chất lượng nước thông qua trí tuệ nhân tạo (Deep Learning); tích hợp chùm vệ tinh quan sát trái đất đa nguồn (VNREDSat, Sentinel-2, Landsat-8/9, vệ tinh viễn thám siêu phổ tương lai); và xây dựng nền tảng WebGIS điện toán đám mây cảnh báo sớm ô nhiễm nguồn nước mặt thời gian thực cho toàn bộ các lưu vực sông chính của Việt Nam.
Kết luận
Luận án tiến sĩ của tác giả Đinh Thị Thu Hiền đã đạt được những thành tựu học thuật và thực tiễn mang tính đột phá:
- Xây dựng thành công quy trình công nghệ hoàn chỉnh và chuẩn hóa trong việc giám sát môi trường nước mặt từ dữ liệu ảnh vệ tinh VNREDSat-1A, khẳng định năng lực tự chủ công nghệ viễn thám không gian của Việt Nam.
- Khám phá và thiết lập hệ thống phương trình hồi quy đa biến tối ưu giữa phổ phản xạ bề mặt 4 kênh ảnh VNREDSat-1A với các thông số chất lượng nước then chốt (Độ đục, TSS, COD, BOD5) đạt độ tin cậy thống kê cao ($R^2$ đạt $0.75 - 0.93$).
- Chứng minh quy luật phân hóa quang phổ thủy vực đặc thù giữa sông phù sa (sông Hồng) và hồ đầm đô thị phú dưỡng (Hà Nội), giải quyết triệt để vấn đề sai số quang học trong môi trường nước Case-2 phức tạp.
- Thành lập bộ bản đồ chuyên đề hiện trạng chất lượng nước mặt Hà Nội cho hai giai đoạn 2016 và 2017, cung cấp bức tranh không gian trực quan, khách quan về diễn biến ô nhiễm môi trường nước Thủ đô theo các ngưỡng phân hạng tiêu chuẩn QCVN 08-MT:2015/BTNMT.
- Mở ra ba hướng nghiên cứu học thuật mới: Viễn thám siêu phổ giám sát ô nhiễm hữu cơ vi sinh; tích hợp học máy nâng cao trong thủy văn không gian; và xây dựng hệ thống cảnh báo sớm thiên tai môi trường nước mặt dựa trên dữ liệu không gian thời gian thực.
- Đặt nền móng ứng dụng thực tiễn bền vững, đóng góp trực tiếp vào công tác quản lý tài nguyên thiên nhiên, an ninh nguồn nước và bảo vệ môi trường sinh thái đô thị trong kỷ nguyên chuyển đổi số.
