Giới thiệu dự án
Xâm nhập mặn tại các vùng cửa sông ven biển hiện nay được xếp vào nhóm tai biến môi trường nghiêm trọng nhất dưới tác động của biến đổi khí hậu (BĐKH) và nước biển dâng. Theo thống kê của Tổ chức Lương thực và Nông nghiệp Liên Hợp Quốc (FAO), mỗi năm thế giới ghi nhận thêm khoảng 2 triệu hecta đất bị nhiễm mặn hóa, làm suy thoái tài nguyên đất, giảm năng suất cây trồng và đe dọa trực tiếp đến an ninh lương thực. Tại Việt Nam, vùng Đồng bằng sông Hồng (ĐBSH) với mạng lưới sông ngòi hạ lưu chịu ảnh hưởng sâu sắc của chế độ nhật triều vịnh Bắc Bộ đang đối mặt với nguy cơ xâm nhập mặn ngày càng phức tạp trong các tháng mùa kiệt (từ tháng 11 đến tháng 4 năm sau).
┌─────────────────────────────────────────────────────────────┐
│ KHÔNG GIAN NGHIÊN CỨU: HUYỆN GIAO THỦY (23.823,80 HA) │
└──────────────────────────────┬──────────────────────────────┘
│
┌───────────────────────┴───────────────────────┐
▼ ▼
┌───────────────────────────┐ ┌───────────────────────────┐
│ HỆ THỐNG THỦY NÔNG │ │ BIẾN ĐỘNG LỚP PHỦ / LULC │
│ - Cửa Ba Lạt & Hà Lạn │ │ - Dữ liệu: Landsat & SPOT │
│ - Xâm mặn: 18 - 37 km │ │ - Phân loại có kiểm định │
│ - Quan trắc: 1989 - 2012 │ │ - Chuỗi thời gian: 5 thời kỳ│
└─────────────┬─────────────┘ └─────────────┬─────────────┘
│ │
└───────────────────────┬───────────────────────┘
▼
┌─────────────────────────────────────────────────────────────┐
│ GIẢI PHÁP: CƠ SỞ DỮ LIỆU GIS & BẢN ĐỒ BIẾN ĐỘNG SỬ DỤNG ĐẤT │
│ TỐI ƯU LỊCH VẬN HÀNH CỐNG & QUY HOẠCH NÔNG - THỦY SẢN│
└─────────────────────────────────────────────────────────────┘
Huyện Giao Thủy, tỉnh Nam Định sở hữu vị trí địa lý đặc thù thuộc vùng hạ lưu sông Hồng, nằm kẹp giữa hai cửa sông lớn là cửa Ba Lạt (sông Hồng) và cửa Hà Lạn (sông Sò) với 32 km bờ biển. Toàn bộ hệ thống thủy nông phục vụ tưới tiêu nông nghiệp và nuôi trồng thủy sản (NTTS) của huyện phụ thuộc vào nguồn nước từ sông Hồng thông qua hệ thống cống đầu mối do Công ty TNHH MTV Khai thác Công trình Thủy lợi (CTTL) Xuân Thủy quản lý. Vấn đề cốt lõi đặt ra là sự xung đột sâu sắc giữa nhu cầu lấy nước ngọt tưới cho lúa vụ Đông Xuân (ĐX) và hiện tượng mặn xâm nhập sâu vào sông chính với nồng độ vượt ngưỡng chịu hạn của cây trồng ($> 1‰$), gây ra hàng loạt thiệt hại:
- Chết lúa diện rộng tại Xuân Thủy (vụ ĐX 1987-1988), Nghĩa Hưng (vụ ĐX 1998-1999).
- Trì hoãn thời vụ gieo cấy do không thể mở cống lấy nước ngọt (vụ ĐX 2002-2003).
- Đỉnh mặn kỷ lục xâm nhập sâu tới cống Hạ Miêu I ($7,2‰$, cách biển 26 km trên sông Hồng) và cống Múc 2 ($1,7‰$, cách biển 37 km trên sông Ninh Cơ) vào tháng 01/2006; nồng độ mặn tại cống Cồn Năm vượt ngưỡng $> 10‰$ vào năm 2010.
Nhằm giải quyết bài toán cấp thiết trên, đề tài khóa luận tập trung triển khai 4 mục tiêu cụ thể:
- Định lượng diễn biến biến đổi độ mặn trên hệ thống kênh tưới tiêu và các cống lấy nước chính thuộc huyện Giao Thủy trong chuỗi thời gian 23 năm (1989 – 2012).
- Xây dựng bộ cơ sở dữ liệu GIS (Geographic Information System) không gian chuẩn hóa về hệ thống công trình thủy lợi (kênh dẫn, cống tưới/tiêu) trên địa bàn nghiên cứu.
- Giải đoán ảnh vệ tinh đa thời gian (Landsat 5 TM, Landsat 7 ETM+, SPOT 5) nhằm thành lập bản đồ hiện trạng và bản đồ biến động sử dụng đất (LULC - Land Use / Land Cover) qua 5 mốc thời gian: 1989, 1995, 2003, 2007 và 2010.
- Đánh giá mối tương quan giữa quy luật xâm nhập mặn và xu thế chuyển dịch cơ cấu sử dụng đất nông nghiệp sang nuôi trồng thủy sản, cung cấp cơ sở khoa học cho quy hoạch không gian thích ứng BĐKH.
Phạm vi nghiên cứu giới hạn trên toàn bộ diện tích tự nhiên $23.823,80\text{ ha}$ của huyện Giao Thủy, tập trung vào tầng mặt nước tưới tiêu và các đơn vị hành chính ven đê, vùng đệm Vườn Quốc gia Xuân Thủy.
Phân tích và thiết kế giải pháp
Phân tích hiện trạng
Trước khi tích hợp công nghệ Viễn thám và GIS, công tác theo dõi độ mặn và quản lý biến động quỹ đất tại địa phương bộc lộ nhiều hạn chế về tính đồng bộ, chi phí và độ phủ không gian.
| Tiêu chí so sánh | Phương pháp quan trắc thủ công truyền thống | Mô phỏng thủy lực số (MIKE 11 HD/AD) | Phương pháp tích hợp Viễn thám & GIS (Đề tài áp dụng) |
|---|---|---|---|
| Độ phủ không gian | Rời rạc, chỉ cục bộ tại vị trí đặt trạm/cống lấy nước | Dọc theo các nút mạng sông 1D | Toàn diện trên diện rộng ($23.823,80\text{ ha}$), chi tiết đến từng pixel |
| Tính cập nhật lịch sử | Phụ thuộc vào sổ tay ghi chép, dễ thất lạc | Khó tái hiện quá khứ nếu thiếu số liệu biên | Truy xuất hồi cứu đa thời gian từ kho tư liệu ảnh vệ tinh 1989–2010 |
| Chi phí vận hành | Tốn kém nhân lực đo đạc trực tiếp định kỳ | Bản quyền phần mềm cao, yêu cầu máy chủ mạnh | Tối ưu chi phí nhờ nguồn ảnh Landsat miễn phí và phần mềm GIS chuyên dụng |
| Khả năng phân tích biến động | Không theo dõi được không gian chuyển đổi đất đai | Không tích hợp sẵn dữ liệu lớp phủ bề mặt | Ma trận chuyển dịch diện tích chéo (Cross-tabulation matrix) chính xác cao |
Phân tích yêu cầu hệ thống theo mô hình MoSCoW:
- Must have (Bắt buộc phải có): Phân loại chính xác 5 nhóm lớp phủ chính; trích xuất ma trận biến động đất đai 5 thời kỳ; xây dựng biểu đồ biến thiên độ mặn trung bình và cực đại tại các cống đầu mối (Cồn Nhất, Cồn Năm, Hạ Miêu).
- Should have (Nên có): Bản đồ số hóa mạng lưới kênh mương thủy lợi tỷ lệ 1/25.000; kiểm định độ chính xác phân loại bằng hệ số Kappa ($\kappa$).
- Could have (Có thể có): Tích hợp chỉ số thực vật NDVI và chỉ số nước LSWI để tách các vùng đất ngập nước/bãi bồi ven sông.
- Won't have (Chưa thực hiện đợt này): Mô phỏng 2 chiều (2D) dòng chảy mặn lan truyền tức thời theo thời gian thực (Real-time telemetry).
Thiết kế hệ thống
Kiến trúc xử lý dữ liệu được thiết kế theo quy trình module khép kín từ tiền xử lý quang học đến phân tích không gian GIS:
[Dữ liệu Vệ tinh (Landsat TM/ETM+, SPOT 5)] ──► [Tiền xử lý: Hiệu chỉnh hình học & bức xạ]
│
▼
[Mẫu ROI khảo sát thực địa] ──────────────────► [Phân loại có kiểm định Maximum Likelihood]
│
▼
[Dữ liệu đo mặn trạm thủy văn (1989-2012)] ────► [Đánh giá độ chính xác (Confusion Matrix / Kappa)]
│
▼
[Số liệu Thủy nông CTTL Xuân Thủy] ────────────► [Tích hợp GIS: Xây dựng Bản đồ Biến động LULC]
Ngăn xếp công nghệ sử dụng
- ENVI version 4.5 / 5.0: Xử lý ảnh số viễn thám, tính toán phản xạ phổ, chọn vùng mẫu (ROIs), thực hiện thuật toán phân loại Maximum Likelihood.
- MapInfo Professional version 10.0: Biên tập bản đồ chuyên đề, số hóa dữ liệu vector ranh giới, kênh mương, tính toán diện tích hình học (Polygon Area Calculation).
- ArcGIS Desktop version 10.2: Tích hợp cơ sở dữ liệu không gian Geodatabase, phân tích chồng lớp không gian (Spatial Overlay Analysis).
- Microsoft Excel / R Script: Xử lý thống kê chuỗi số liệu độ mặn thủy văn 1989–2012.
Cơ sở toán học và thuật toán cốt lõi
-
Công thức xác định độ mặn dựa trên hàm lượng Ion Clorua ($\text{Cl}^-$): Trong môi trường nước biển và cửa sông pha trộn, độ muối và độ mặn ($S‰$) được xác định thông qua hàm lượng Ion $\text{Cl}^-$ ($g/\text{kg}$ nước):
$$\text{Độ muối: } S(‰) = 0{,}030 + 1{,}8050 \times \text{Cl}^-(‰)$$
$$\text{Độ mặn thực tế trong sông: } S = 1{,}65 \times \text{Cl}^- \quad \Longleftrightarrow \quad \text{Cl}^- = 0{,}607 \times S$$
-
Chỉ số khác biệt thực vật chuẩn hóa (NDVI):
$$\text{NDVI} = \frac{\text{NIR} - \text{Red}}{\text{NIR} + \text{Red}} = \frac{\text{Band 4} - \text{Band 3}}{\text{Band 4} + \text{Band 3}} \quad (\text{đối với Landsat 5 TM})$$
-
Chỉ số nước bề mặt (LSWI):
$$\text{LSWI} = \frac{\text{NIR} - \text{SWIR}}{\text{NIR} + \text{SWIR}} = \frac{\text{Band 4} - \text{Band 5}}{\text{Band 4} + \text{Band 5}}$$
-
Đánh giá độ chính xác phân loại ảnh bằng hệ số Kappa ($\hat{K}$):
$$\hat{K} = \frac{N \sum_{i=1}^{k} x_{ii} - \sum_{i=1}^{k} (x_{i+} \cdot x_{+i})}{N^2 - \sum_{i=1}^{k} (x_{i+} \cdot x_{+i})}$$
Trong đó: $N$ là tổng số pixel mẫu kiểm tra; $x_{ii}$ là giá trị trên đường chéo chính của ma trận sai số (Confusion Matrix); $x_{i+}$ và $x_{+i}$ là tổng hàng $i$ và tổng cột $i$.
Methodology
Quy trình nghiên cứu áp dụng phương pháp luận kết hợp liên ngành (Interdisciplinary Methodology):
- Giai đoạn 1 - Thu thập và chuẩn hóa dữ liệu: Thu thập dữ liệu ảnh viễn thám đa phổ Landsat (1989, 1995, 2003, 2007, 2010), ảnh SPOT 5 độ phân giải cao ($10\text{ m}$); chuỗi số liệu quan trắc độ mặn thực tế tại các cống ven sông Hồng và sông Sò từ năm 1989 đến 2012.
- Giai đoạn 2 - Xử lý ảnh số & Kiểm định mẫu: Khử nhiễu khí quyển, nắn chỉnh tọa độ về lưới chiếu UTM WGS84 múi 48N. Lấy mẫu thực địa cho các đối tượng: Đất sản xuất nông nghiệp (SXN - lúa, hoa màu), Đất nuôi trồng thủy sản (NTS), Đất thổ cư/khu dân cư, Đất mặt nước (sông, kênh mương) và Đất bãi bồi/rừng ngập mặn.
- Giai đoạn 3 - Phân tích chồng lớp & Đánh giá thủy văn: Sử dụng ma trận biến động để bóc tách luồng chuyển đổi giữa các loại hình đất; kết nối với xu thế biến đổi mặn theo chu kỳ nhật triều và diễn biến đóng mở cống thủy nông.
Implementation và kết quả
Development process
Quá trình trích xuất thông tin không gian và chuyển đổi vector được tự động hóa thông qua pipeline xử lý ảnh số. Đoạn mã Python (sử dụng thư viện rasterio, numpy và scikit-learn) minh họa thuật toán phân loại và tính ma trận chuyển dịch đất đai:
import numpy as np
import rasterio
from sklearn.metrics import confusion_matrix, cohen_kappa_score
def compute_land_cover_transition(raster_t1_path: str, raster_t2_path: str):
"""
Tính toán ma trận chuyển dịch sử dụng đất giữa 2 thời kỳ (t1 -> t2).
"""
with rasterio.open(raster_t1_path) as src1, rasterio.open(raster_t2_path) as src2:
grid_t1 = src1.read(1)
grid_t2 = src2.read(1)
nodata = src1.nodata
# Lọc bỏ pixel NoData
valid_mask = (grid_t1 != nodata) & (grid_t2 != nodata)
t1_valid = grid_t1[valid_mask]
t2_valid = grid_t2[valid_mask]
# Tính toán ma trận biến động (Cross-tabulation matrix)
classes = np.unique(np.concatenate((t1_valid, t2_valid)))
transition_matrix = confusion_matrix(t1_valid, t2_valid, labels=classes)
return classes, transition_matrix
def calculate_accuracy_metrics(ground_truth: np.ndarray, classified: np.ndarray):
"""
Tính Overall Accuracy và Chỉ số Kappa
"""
oa = np.mean(ground_truth == classified) * 100.0
kappa = cohen_kappa_score(ground_truth, classified)
return oa, kappa
Testing và validation
Kết quả giải đoán ảnh viễn thám qua các năm được kiểm định thông qua ma trận sai số đối chiếu với số liệu điều tra thực địa và bản đồ kiểm kê đất đai huyện Giao Thủy:
- Độ chính xác tổng thể (Overall Accuracy): Đạt $91{,}20%$ trên bộ mẫu kiểm định năm 2010.
- Chỉ số Kappa ($\hat{K}$): Đạt giá trị dao động từ $0{,}82$ đến $0{,}89$ qua các cảnh ảnh Landsat 1989, 1995, 2003, 2007 và đạt $\sim 0{,}91$ trên ảnh SPOT 5.
+-------------------------------------------------------------------------+
| MA TRẬN KIỂM ĐỊNH ĐỘ CHÍNH XÁC (2010) |
| Overall Accuracy: 91.20% | Kappa Coefficient: 0.887 |
+--------------------------+-----------------------+----------------------+
| Loại hình lớp phủ | Producer Accuracy (%) | User Accuracy (%) |
+--------------------------+-----------------------+----------------------+
| Đất sản xuất lúa (SXN) | 92.45% | 90.15% |
| Đất NTTS nước lợ (NTS) | 94.10% | 93.80% |
| Đất khu dân cư / thổ cư | 88.30% | 87.50% |
| Đất rừng ngập mặn/Bãi bồi| 93.00% | 94.20% |
| Mặt nước sông ngòi | 98.20% | 97.60% |
+--------------------------+-----------------------+----------------------+
Kết quả đạt được
1. Quy luật biến đổi độ mặn trên hệ thống thủy nông (1989 - 2012)
- Tương quan triều - mặn: Vịnh Bắc Bộ có chế độ nhật triều thuần nhất với biên độ triều cực đại đạt tới $4{,}5\text{ m}$ tại Hòn Dấu; thời gian truyền sóng triều trên sông Hồng đạt tốc độ $15 - 20\text{ km/h}$. Đỉnh mặn xuất hiện chậm hơn đỉnh triều trung bình $1\text{ giờ}$, trong khi chân mặn xuất hiện đồng thời cùng chân triều.
- Biến động không gian và thời gian:
- Tốc độ truyền đỉnh mặn vùng cửa sông Hồng đạt $12{,}5\text{ km/h}$.
- Nồng độ mặn cao nhất trong năm luôn rơi vào các tháng mùa kiệt (tháng 1 đến tháng 3). Nồng độ mặn trung bình tại thời điểm lấy nước của các cống dọc sông Hồng có xu thế gia tăng rõ rệt từ năm 1989 đến 2012.
- Cống Cồn Nhất và Cồn Năm ghi nhận các đợt xâm nhập mặn đột biến $> 3‰$ đến $> 10‰$, buộc công ty thủy nông phải đóng cống hoàn toàn trong nhiều chu kỳ triều cường, gây thiếu nước ngọt nghiêm trọng cho $9.778{,}31\text{ ha}$ đất canh tác nông nghiệp.
Độ mặn (‰)
12 ┼ ▲ (2010: >10‰)
10 ┼ /
8 ┼ ▲ (2006: 7.2‰)
6 ┼ /
4 ┼ ▲ /
2 ┼ ▲ / \ /
0 ┼────────┴────────────┴───┴─────────────┴─────────────────
1989 1995 2003 2007 2012 (Năm)
── Xu thế độ mặn cực đại tại các cống đầu mối ven sông Hồng ──
2. Biến động sử dụng đất huyện Giao Thủy (1989 - 2010)
Tổng diện tích tự nhiên của huyện là $23.823{,}80\text{ ha}$. Quá trình phân tích ảnh vệ tinh đa thời gian đã chỉ ra sự tái cơ cấu không gian đất đai mạnh mẽ:
| Loại hình sử dụng đất | Mã | Năm 1989 (ha) | Năm 2003 (ha) | Năm 2010 (ha) | Cơ cấu 2010 (%) | Biến động 1989–2010 (%) |
|---|---|---|---|---|---|---|
| Đất sản xuất nông nghiệp | SXN | 13.250,40 | 11.025,10 | 9.778,31 | 41,04% | -26,20% |
| Đất nuôi trồng thủy sản | NTS | 840,50 | 2.450,20 | 3.875,34 | 16,27% | +361,08% |
| Đất làm muối | LMU | 780,20 | 620,15 | 515,46 | 2,16% | -33,93% |
| Đất sông suối, mặt nước | SMN | 1.980,00 | 1.890,50 | 1.852,39 | 7,78% | -6,44% |
| Đất khu dân cư / thổ cư | OCT | 2.150,00 | 2.480,30 | 2.719,65 | 11,42% | +26,50% |
| Đất bãi bồi / chưa sử dụng | CSD | 2.510,70 | 1.850,20 | 1.156,08 | 4,85% | -53,95% |
| Các loại đất khác | KHK | 2.312,00 | 3.507,35 | 3.926,57 | 16,48% | +69,83% |
| Tổng diện tích tự nhiên | TONG | 23.823,80 | 23.823,80 | 23.823,80 | 100,00% | 0,00% |
- Xu thế chuyển dịch nông nghiệp - thủy sản: Diện tích đất sản xuất nông nghiệp giảm liên tục $3.472{,}09\text{ ha}$ (giảm $26{,}20%$), trong khi đất nuôi trồng thủy sản tăng vọt từ $840{,}50\text{ ha}$ lên $3.875{,}34\text{ ha}$ (tăng $361{,}08%$, tương đương gấp 4,6 lần). Hiện tượng xâm nhập mặn khiến việc trồng lúa kém hiệu quả tại các xã ven biển (Giao Thiện, Giao An, Giao Lạc, Giao Xuân) đã thúc đẩy người dân chủ động chuyển đổi sang mô hình nuôi tôm sú, cua, cá vược nước lợ.
- Tốc độ bồi tụ phù sa ven biển: Vùng bãi bồi cửa Ba Lạt và Vườn Quốc gia Xuân Thủy được bồi lắng tự nhiên từ $150 - 200\text{ ha/năm}$, tạo ra các cồn cát mới (Cồn Lu, Cồn Ngạn) làm thay đổi đáng kể đường bờ biển và chế độ thủy động lực học.
Đổi mới và đóng góp
- Đổi mới phương pháp luận kết hợp: Khắc phục nhược điểm của các nghiên cứu đơn ngành (hoặc chỉ quan sát trạm đo mặn thủy văn, hoặc chỉ giải đoán ảnh viễn thám độc lập). Đề tài thiết lập mối liên kết nhân quả định lượng giữa động thái xâm nhập mặn trên hệ thống cống thủy lợi và tốc độ chuyển đổi cơ cấu lớp phủ đất đai.
- Chuỗi dữ liệu lịch sử kéo dài 23 năm: Cung cấp bộ bản đồ số hóa LULC qua 5 thời kỳ (1989, 1995, 2003, 2007, 2010) với độ phân giải nhất quán, làm sáng tỏ quá trình đô thị hóa và thủy sản hóa nông nghiệp tại vùng ĐBSH.
- Cơ sở khoa học cho vận hành công trình thủy lợi: Đưa ra dữ liệu thực nghiệm về độ trễ pha giữa đỉnh triều và đỉnh mặn ($1\text{ giờ}$), hỗ trợ các kỹ sư thủy nông lập quy trình mở cống đón đỉnh triều lấy nước ngọt chính xác, giảm thiểu rủi ro lấy nhầm nước mặn gây chết mùa màng.
Ứng dụng thực tế và triển khai
Kịch bản ứng dụng thực tế
- Quản lý vận hành cống thủy nông: Công ty TNHH MTV Khai thác CTTL Xuân Thủy ứng dụng biểu đồ dao động mặn để lập lịch vận hành đóng/mở tự động theo chu kỳ triều kém và triều cường.
- Quy hoạch vùng sản xuất thích ứng: Phòng Tài nguyên & Môi trường, Phòng Nông nghiệp huyện Giao Thủy sử dụng bản đồ biến động đất đai để quy hoạch phân ranh giới cố định giữa vùng chuyên lúa nước ngọt nội đồng và vùng chuyên canh NTTS nước lợ ven đê, tránh tình trạng tự phát đào ao nuôi tôm gây mặn hóa đồng ruộng.
QUY TRÌNH VẬN HÀNH ĐÓNG / MỞ CỐNG DỰA TRÊN QUY LUẬT TRIỀU - MẶN
[ Quan trắc đỉnh triều tại trạm ]
│
▼
[ Đỉnh triều xuất hiện (T_0) ] ──► [ Đỉnh mặn xuất hiện (T_0 + 1 giờ) ]
│ │
▼ ▼
[ Mở cống lấy nước tự chảy ] [ ĐÓNG CỐNG HOÀN TOÀN ]
(Thời điểm triều dâng mặn thấp) (Khi độ mặn đo vượt ngưỡng 1.0‰)
Lộ trình triển khai khuyến nghị
- Giai đoạn 1 (0 - 6 tháng): Tích hợp toàn bộ bản đồ chuyên đề định dạng Shapefile/MapInfo vào hệ thống thông tin địa lý cấp huyện.
- Giai đoạn 2 (6 - 18 tháng): Lắp đặt mạng lưới cảm biến đo độ dẫn điện/độ mặn (EC/Salinity Sensor) tự động truyền dữ liệu qua mạng viễn thông (GSM/IoT) tại 5 cống đầu mối chính.
- Giai đoạn 3 (18 - 36 tháng): Xây dựng mô hình ra quyết định tự động hỗ trợ điều tiết nước đa mục tiêu (Tưới lúa - Rửa mặn - Cấp nước nuôi trồng thủy sản).
Hạn chế và hướng phát triển
- Hạn chế kỹ thuật: Nguồn ảnh quang học Landsat trong quá khứ bị ảnh hưởng cục bộ bởi mây và sương mù mùa đông xuân ở miền Bắc; chưa tích hợp ảnh vệ tinh Radar giao thoa (SAR/Sentinel-1) để theo dõi độ ẩm đất và ngập lụt dưới mọi điều kiện thời tiết.
- Dữ liệu thủy văn: Số liệu đo mặn tại một số cống nhánh nội đồng chưa được đo đạc liên tục 24/24 theo từng giờ trong ngày mà chủ yếu đo vào thời điểm mở cống lấy nước.
- Hướng nghiên cứu tiếp theo:
- Ứng dụng thuật toán học máy (Random Forest, Support Vector Machine) trên nền tảng Google Earth Engine (GEE) để tự động hóa phân loại LULC chuỗi thời gian thực.
- Xây dựng mô hình thủy lực 2 chiều kết hợp (Coupled 2D Hydrodynamic - Salinity Transport Model) mô phỏng chi tiết sự lan truyền mặn trên toàn hệ thống sông Hồng - sông Đáy khi xả nước từ các hồ chứa thượng nguồn (Hòa Bình, Tuyên Quang).
Đối tượng hưởng lợi
- Học viên, sinh viên ngành Môi trường/Địa lý/GIS: Nguồn tài liệu tham khảo hoàn chỉnh về phương pháp luận xử lý ảnh Landsat/SPOT và thuật toán ma trận biến động.
- Kỹ sư GIS & Quản lý dữ liệu tài nguyên: Cung cấp quy trình làm sạch, số hóa bản đồ công trình thủy lợi và phương pháp kiểm định sai số phân loại viễn thám.
- Công ty Khai thác Công trình Thủy lợi: Nắm bắt chính xác quy luật truyền mặn theo pha triều để tối ưu hóa quy trình lấy nước ngọt, bảo vệ cây trồng vụ Đông Xuân.
- Các nhà hoạch định chính sách & Cơ quan quản lý nhà nước: Khung dữ liệu định lượng phục vụ xây dựng Kế hoạch hành động ứng phó với Biến đổi khí hậu tỉnh Nam Định tầm nhìn 2030 - 2050.
Câu hỏi thường gặp
1. Yêu cầu kỹ thuật phần cứng và phần mềm để tái lập nghiên cứu này là gì?
Hệ thống máy tính tối thiểu cần trang bị CPU 4 nhân (từ Intel Core i5 thế hệ 4 trở lên), 8GB RAM (khuyến nghị 16GB để xử lý ảnh đa băng phổ), ổ cứng lưu trữ tối thiểu 50GB trống. Về phần mềm: Cần cài đặt ENVI (phiên bản 4.5 hoặc mới hơn), MapInfo Professional (v10.0+) hoặc QGIS/ArcGIS (v10.2+) có hỗ trợ chuyển đổi hệ tọa độ VN-2000 và UTM WGS84.
2. Phương pháp phân loại Maximum Likelihood xử lý sự nhầm lẫn phổ giữa ao nuôi thủy sản và ruộng lúa ngập nước như thế nào?
Ruộng lúa chuẩn bị cấy vụ Đông Xuân và ao nuôi tôm đều có đặc trưng phản xạ phổ nước mặt (giá trị phổ hấp thụ mạnh ở dải hồng ngoại gần NIR). Để tách biệt, nghiên cứu sử dụng phương pháp kết hợp chỉ số thực vật thời gian đa vụ kết hợp kênh tỷ số $\text{NDVI} / \text{LSWI}$ và phân tích ngữ cảnh không gian (Spatial Context), đồng thời đối chiếu ảnh SPOT 5 độ phân giải cao $10\text{ m}$ để nhận diện bờ bao ngăn cách hình học đặc trưng của ao NTTS.
3. Làm thế nào để tích hợp cơ sở dữ liệu GIS này vào các hệ thống SCADA quản lý nước hiện đại?
Dữ liệu vector (Spatial Layers) của đề tài được lưu trữ dưới dạng chuẩn ESRI Shapefile và Geodatabase (chứa thuộc tính tọa độ, khẩu độ cống, cao trình đáy cống, lưu lượng thiết kế). Hệ thống SCADA có thể kết nối trực tiếp thông qua chuẩn dịch vụ Web Map Service (WMS) hoặc Web Feature Service (WFS) theo chuẩn OGC để hiển thị trạng thái đóng/mở và nồng độ mặn thời gian thực lên bản đồ nền GIS.
4. Chi phí duy tu và cập nhật cơ sở dữ liệu viễn thám định kỳ là bao nhiêu?
Chi phí cập nhật rất thấp nhờ việc tận dụng nguồn dữ liệu vệ tinh mở miễn phí chất lượng cao như Sentinel-2 ($10\text{ m}$) và Landsat 8/9 ($30\text{ m}$) của USGS/ESA. Chi phí vận hành định kỳ chủ yếu dành cho nhân sự kiểm định thực địa và bảo dưỡng các đầu mối quan trắc độ mặn (khoảng 20 - 50 triệu VNĐ/năm cho toàn huyện).
5. Khả năng nhân rộng mô hình này cho các huyện ven biển khác như thế nào?
Mô hình hoàn toàn có thể nhân rộng $100%$ cho các địa bàn ven biển lân cận có đặc điểm tương đồng tại miền Bắc như Hải Hậu, Nghĩa Hưng (Nam Định), Tiền Hải, Thái Thụy (Thái Bình) hoặc các vùng cửa sông Đồng bằng sông Cửu Long bằng cách thay đổi số liệu biên thủy văn và biên tập lại vùng ROI mẫu mặt đất tương ứng.
Kết luận
Luận văn đã giải quyết trọn vẹn bài toán thực tiễn về đánh giá định lượng biến đổi độ mặn và diễn biến sử dụng đất tại huyện Giao Thủy, tỉnh Nam Định giai đoạn 1989 – 2012. Việc kết hợp chặt chẽ giữa chuỗi quan trắc thủy văn $S(‰)$ dài hạn và công nghệ Viễn thám - GIS đã xác lập bức tranh tổng thể về mối liên hệ nhân quả giữa suy giảm nguồn nước ngọt do mặn xâm nhập sâu ($18 - 37\text{ km}$) và sự thu hẹp $26{,}20%$ diện tích lúa nước để chuyển dịch sang $3.875{,}34\text{ ha}$ nuôi trồng thủy sản nước lợ. Kết quả nghiên cứu không chỉ mang giá trị học thuật về phương pháp xử lý dữ liệu địa không gian mà còn là công cụ hỗ trợ ra quyết định trực tiếp cho ngành thủy lợi và chính quyền địa phương trong chiến lược quy hoạch kinh tế nông nghiệp bền vững, chủ động thích ứng với biến đổi khí hậu toàn cầu.