Giới thiệu dự án
Khu vực đồng bằng sông Cửu Long (ĐBSCL) đang đối mặt với những biến đổi sâu sắc về hình thái ven biển dưới tác động kép của biến đổi khí hậu, nước biển dâng và các hoạt động nhân sinh như giảm thiểu nguồn bùn cát thượng nguồn và khai thác cát quá mức. Tại tỉnh Bến Tre, huyện Thạnh Phú sở hữu dải bờ biển dài 25 km nằm kẹp giữa hai cửa sông lớn là cửa Hàm Luông và cửa Cổ Chiên. Đây là vùng sinh thái nhạy cảm, chịu chế độ bán nhật triều không đều với biên độ triều lớn (2,5 – 3,0 m), thường xuyên xảy ra các quá trình động lực học phức tạp gây xói lở và bồi tụ bờ biển.
+-------------------------------------------------------+
| KHU VỰC NGHIÊN CỨU: BỜ BIỂN HUYỆN THẠNH PHÚ (25 KM) |
+-------------------------------------------------------+
|
+---------------------+---------------------+
| |
[Cửa sông Cổ Chiên] [Cửa sông Hàm Luông]
(Đoạn 1: Xói lở TB) (Đoạn 4: Biến động nhẹ)
| |
+---------------------+---------------------+
|
+---------------------+---------------------+
| |
[Đoạn 2: Cổ Chiên - Cồn Bửng] [Đoạn 3: Cồn Bửng - Hàm Luông]
(Bồi tụ lớn: +24,95 m/năm) (Xói lở mạnh: -20,63 m/năm)
Sạt lở bờ biển làm mất đất canh tác, phá hủy hạ tầng đê kè ven biển và suy thoái hệ sinh thái rừng ngập mặn phòng hộ. Ngược lại, hiện tượng bồi tụ làm biến dạng các luồng tàu thuyền và cửa sông. Do đó, việc giám sát định lượng và liên tục biến động đường bờ là yêu cầu cấp thiết phục vụ công tác quản lý tài nguyên, quy hoạch lãnh thổ và giảm thiểu rủi ro thiên tai.
Đồ án tập trung vào các mục tiêu cụ thể:
- Thu thập, tiền xử lý và đồng bộ hóa chuỗi dữ liệu ảnh vệ tinh quang học Landsat đa thời gian (2006, 2011, 2015).
- Rút trích tự động và bán tự động vị trí đường bờ biển qua các thời kỳ bằng kỹ thuật tỷ số kênh phổ kết hợp giá trị ngưỡng tối ưu.
- Ứng dụng công cụ Digital Shoreline Analysis System (DSAS) để tính toán thống kê định lượng các chỉ số biến động không gian đường bờ (EPR, LRR, NSM, SCE, LR²).
- Phân vùng, đánh giá chi tiết xu thế xói lở – bồi tụ trên 4 phân đoạn bờ biển Thạnh Phú, cung cấp luận cứ khoa học cho các giải pháp công trình và phi công trình.
Giải pháp ứng dụng tích hợp công nghệ Viễn thám (Remote Sensing - RS) và Hệ thống Thông tin Địa lý (GIS) mang lại cách tiếp cận khách quan, bao phủ không gian lớn với chi phí tối ưu so với trắc địa truyền thống. Phạm vi nghiên cứu giới hạn trên 25 km đường bờ huyện Thạnh Phú từ Vàm Rỏng đến Khâu Băng thuộc hai xã Thạnh Phong và Thạnh Hải trong giai đoạn 9 năm (2006 – 2015). Nghiên cứu thừa nhận giới hạn kỹ thuật về độ phân giải không gian 30 m của cảm biến Landsat và việc giản lược hiệu chỉnh triều tức thời do thiếu mô hình số độ cao địa hình đáy biển chi tiết.
Phân tích và thiết kế giải pháp
Phân tích hiện trạng
Trước đây, việc theo dõi đường bờ chủ yếu dựa trên các phương pháp khảo sát đo đạc thực địa định kỳ (sử dụng máy toàn đạc điện tử, định vị vệ tinh GPS) hoặc giải đoán thủ công trên ảnh hàng không. Các phương pháp này bộc lộ nhiều hạn chế khi áp dụng trên quy mô vùng rộng lớn:
| Tiêu chí |
Khảo sát thực địa (GPS/Toàn đạc) |
Giải đoán thủ công ảnh viễn thám |
Tỷ số phổ kết hợp DSAS (Đề tài) |
| Độ phủ không gian |
Rất hẹp, chỉ theo các mặt cắt điểm |
Trung bình |
Rộng, bao phủ liên tục toàn dải bờ |
| Độ chính xác vị trí |
Rất cao ($< 0,1\text{ m}$) |
Phụ thuộc chủ quan chuyên gia |
Cao, chuẩn hóa thuật toán ($30\text{ m}$) |
| Chi phí & Thời gian |
Rất tốn kém, nhân lực lớn |
Tốn thời gian số hóa thủ công |
Tối ưu, tự động hóa cao |
| Tần suất lặp lại |
Thấp, khó xây dựng chuỗi thời gian |
Trung bình |
Cao, tận dụng kho lưu trữ ảnh vệ tinh |
| Khả năng thống kê |
Thủ công qua bảng tính |
Rời rạc, thiếu đồng bộ |
Tự động tính toán LRR, NSM, EPR |
Yêu cầu kỹ thuật của giải pháp được phân loại theo mô hình MoSCoW:
- Must have (Bắt buộc): Rút trích đường bờ từ ảnh Landsat 5 TM và Landsat 8 OLI/TIRs; xây dựng lớp cơ sở dữ liệu đường bờ vector chuẩn; thiết lập Baseline và tạo hệ thống Transects vuông góc; tính toán tốc độ xói lở/bồi tụ theo phương pháp hồi quy tuyến tính (LRR).
- Should have (Nên có): Điều chỉnh ngưỡng phổ riêng biệt cho cảm biến Landsat 8 OLI để triệt tiêu sai số phổ; tính toán hệ số xác định ($R^2$ hay $LR^2$) nhằm đánh giá mức độ tin cậy của xu thế.
- Could have (Có thể mở rộng): Tích hợp dữ liệu trạm đo mực nước triều để nội suy mép nước trung bình; xây dựng kịch bản mô phỏng lan truyền đường bờ tương lai.
- Won't have (Chưa thực hiện): Tích hợp dữ liệu mô hình sóng - dòng chảy 3D (như Delft3D, MIKE 21) trong khuôn khổ đồ án này.
+-------------------------------------------------------------------------+
| KIẾN TRÚC HỆ THỐNG |
+-------------------------------------------------------------------------+
|
+----------------------------------+----------------------------------+
| |
+------------------+ +-------------------+
| TẦNG DỮ LIỆU | | TẦNG TIỀN XỬ LÝ |
| - Landsat 5 TM | ===> [Gom kênh / Cắt vùng WGS-84 UTM 48N] ===> | - ENVI 4.7 Engine |
| - Landsat 8 OLI | | - IDL Processing |
+------------------+ +-------------------+
|
+---------------------------------------------------------------------+
|
+------------------------------------------------------------------------------+
| TẦNG THUẬT TOÁN TÁCH ĐƯỜNG BỜ (Gathot Winasor Ratio) |
| - TM: (B4/B2 < 1) & ((B5+B7)/B2 < 1) |
| - OLI: (B5/B3 < 1) & ((B6+B7)/B3 < 1.2376) |
+------------------------------------------------------------------------------+
|
+-------------------------------+-------------------------------+
| |
+------------------------------------+ +-------------------+
| TẦNG XỬ LÝ KHÔNG GIAN (ArcGIS 10.1)| | TẦNG PHÂN TÍCH |
| - Raster to Polygon/Polyline | ===> [Mạng lưới ===>| THỐNG KÊ (DSAS) |
| - Polynomial Smoothing (100 m) | Transects 50m] | - EPR, LRR, NSM |
| - Geodatabase Integration | | - Hệ số LR2 |
+------------------------------------+ +-------------------+
Thiết kế hệ thống
Hệ thống xử lý phân tích được xây dựng trên nền tảng tích hợp giữa phần mềm xử lý ảnh số chuyên dụng và hệ thống thông tin địa lý:
- Hệ điều hành: Microsoft Windows 7/10 64-bit.
- Nền tảng xử lý ảnh viễn thám: ENVI 4.7 (Environment for Visualizing Images) chạy trên nhân ngôn ngữ lập trình tương tác IDL (Interactive Data Language), đảm bảo xử lý đồng thời nhiều kênh phổ mà không làm suy giảm hiệu năng.
- Nền tảng GIS & Không gian: ESRI ArcGIS Desktop 10.1 với cơ sở dữ liệu Geodatabase.
- Module mở rộng chuyên dụng: USGS DSAS v4.3 (Digital Shoreline Analysis System) chạy nhúng trong môi trường ArcMap.
- Hệ quy chiếu chuẩn: WGS-84 UTM Zone 48N (EPSG: 32648).
Cấu trúc cơ sở dữ liệu không gian (Geodatabase Schema):
Shorelines_Merged (Feature Class - Polyline): Chứa chuỗi đường bờ đa thời gian với các trường bắt buộc: Date_ (Text - MM/DD/YYYY), Uncertainty (Double - giá trị sai số vị trí, gán $30\text{ m}$ theo kích thước pixel Landsat).Baseline_Offshore (Feature Class - Polyline): Đường cơ sở đặt ngoài khơi song song với xu thế chung của đường bờ, làm điểm mốc phóng các tia cắt (ID, Shape, Group, Buff_Dist).DSAS_Transects (Feature Class - Polyline): Hệ thống các đường cắt vuông góc phát sinh tự động từ Baseline với khoảng cách bước $50\text{ m}$, cắt qua toàn bộ các lớp đường bờ.DSAS_Results_Table (Table): Bảng thuộc tính liên kết lưu trữ kết quả thống kê toán học (TransectID, Distance, EPR, LRR, LR2, NSM, SCE).
Methodology
Quy trình nghiên cứu tuân thủ phương pháp luận khoa học chặt chẽ qua 4 giai đoạn nối tiếp:
[Giai đoạn 1: Thu thập & Chuẩn hóa dữ liệu]
|
v
[Giai đoạn 2: Trích xuất đường bờ trên ENVI]
|
v
[Giai đoạn 3: Phân tích động lực đường bờ trên GIS & DSAS]
|
v
[Giai đoạn 4: Đánh giá tác động & Báo cáo kết quả]
- Giai đoạn 1 (Tuần 1 - 3): Thu thập ảnh viễn thám từ USGS EarthExplorer. Lựa chọn ảnh có độ che phủ mây thấp ($< 30%$), thời điểm chụp từ 14:00 - 15:00 giờ địa phương, mùa khô (tháng 1 - tháng 2) để giảm thiểu biến động do mây mù và lũ thượng nguồn.
- Giai đoạn 2 (Tuần 4 - 7): Gom kênh phổ (Layer Stacking), cắt vùng quan tâm (Spatial Subset), lập chỉ số tỷ số kênh và áp dụng ngưỡng nhị phân để phân tách ranh giới đất/nước.
- Giai đoạn 3 (Tuần 8 - 10): Chuyển đổi dữ liệu Raster sang Vector, làm trơn polyline bằng thuật toán làm mịn với bán kính lọc $100\text{ m}$, thiết lập Geodatabase và cấu hình module DSAS để chạy phân tích thống kê.
- Giai đoạn 4 (Tuần 11 - 12): Tổng hợp số liệu, vẽ biểu đồ phân bố LRR, lập bản đồ biến động và đánh giá nguyên nhân bồi - xói phục vụ phát triển kinh tế - xã hội địa phương.
Implementation và kết quả
Development process
Quá trình trích xuất đường bờ sử dụng phương pháp tỷ số kênh phổ của Gathot Winasor kết hợp giá trị ngưỡng. Thuật toán này tách biệt hiệu quả ranh giới đất liền và nước biển thông qua việc phối hợp giữa kênh phản xạ mạnh của thực vật (NIR) và các kênh hấp thụ nước mạnh (SWIR, Green).
# Mô phỏng logic Band Math xử lý tách đường bờ đa thời gian trong ENVI / Python Rasterio
import numpy as np
def extract_shoreline_tm(b2_green, b4_nir, b5_swir1, b7_swir2):
"""
Rút trích ranh giới nước/đất cho cảm biến Landsat 5 TM.
Tỷ số 1: B4/B2 dùng phân tách vùng bờ có thảm thực vật che phủ.
Tỷ số 2: (B5+B7)/B2 dùng phân tách vùng bờ không có thực vật (bãi cát, bùn).
"""
# Tránh chia cho 0
b2_safe = np.where(b2_green == 0, 0.0001, b2_green)
ratio_veg = b4_nir / b2_safe
ratio_non_veg = (b5_swir1 + b7_swir2) / b2_safe
# Điều kiện mặt nước: cả 2 tỷ số đều có giá trị < 1
water_mask = (ratio_veg < 1.0) & (ratio_non_veg < 1.0)
# 0: Vùng nước biển; 1: Vùng đất liền bờ biển
land_mask = np.where(water_mask, 0, 1)
return land_mask
def extract_shoreline_oli(b3_green, b5_nir, b6_swir1, b7_swir2):
"""
Rút trích ranh giới nước/đất cho cảm biến Landsat 8 OLI/TIRs.
Hiệu chỉnh ngưỡng phổ thực nghiệm: ngưỡng tỷ số non-veg nâng lên 1.2376
nhằm bù trừ độ dịch chuyển phổ và độ rộng kênh của cảm biến 11 kênh.
"""
b3_safe = np.where(b3_green == 0, 0.0001, b3_green)
ratio_veg = b5_nir / b3_safe
ratio_non_veg = (b6_swir1 + b7_swir2) / b3_safe
# Điều kiện phân tách nước/đất chuẩn hóa cho Landsat 8
water_mask = (ratio_veg < 1.0) & (ratio_non_veg < 1.2376)
land_mask = np.where(water_mask, 0, 1)
return land_mask
Sau khi tạo ảnh nhị phân (Binary Mask), ranh giới giữa pixel giá trị 0 và 1 được trích xuất sang định dạng Raster Contour, chuyển đổi sang đối tượng Vector Polygon và trích xuất thành Polyline. Do hạn chế về răng cưa của ảnh raster $30\text{ m}$, thuật toán làm trơn (Line Smoothing) với tham số dung sai $100\text{ m}$ được áp dụng trong ArcGIS để loại bỏ các điểm dị thường cục bộ do sóng vỡ ven bờ.
+-------------------------------------------------------------+
| CÁC CHỈ SỐ THỐNG KÊ BIẾN ĐỘNG ĐƯỜNG BỜ TRONG CÔNG CỤ DSAS |
+-------------------------------------------------------------+
|
+----------------------------+----------------------------+
| |
+------------------------------+ +-----------------------------+
| End Point Rate (EPR) | | Linear Regression Rate (LRR)|
| - Công thức: | | - Công thức: |
| EPR = (D_mới - D_cũ) / dt | | y = a + b * x |
| - Đơn vị: mét/năm | | LRR = b (độ dốc đường HQ) |
+------------------------------+ +-----------------------------+
| |
+----------------------------+----------------------------+
|
+--------------------------------------------------------------------------------------+
| Hệ số xác định hồi quy (LR2 - R-squared) |
| - Công thức: LR2 = 1 - (SS_res / SS_tot) |
| - Ý nghĩa: Đánh giá độ tin cậy tuyến tính của xu thế bồi/xói dọc theo từng Transect |
+--------------------------------------------------------------------------------------+
Testing và validation
Kết quả đường bờ vector rút trích được kiểm chứng chéo với các ảnh tổ hợp màu thực và màu giả:
- Đối với Landsat 5 TM (2006, 2011): Tổ hợp màu RGB 5-4-3 (SWIR1 - NIR - Red), thể hiện rõ nét độ tương phản giữa ranh giới nước (hấp thụ màu tối) và thảm rừng ngập mặn (màu đỏ tươi).
- Đối với Landsat 8 OLI (2015): Tổ hợp màu RGB 6-5-4 (SWIR1 - NIR - Red).
Độ tin cậy của mô hình toán DSAS được đánh giá thông qua chỉ số xác định $LR^2$. Tại các phân đoạn có biến động diễn ra đơn chiều liên tục (như đoạn xói lở Cồn Bửng - Hàm Luông), hệ số $LR^2$ đạt giá trị trung bình $0,88$, chứng minh tính tuyến tính rất cao của quá trình xâm thực bờ biển.
Kết quả đạt được
Nghiên cứu chia bờ biển Thạnh Phú (25 km) thành 4 đoạn đặc trưng để định lượng biến động:
- Đoạn 1: Đoạn bờ cửa sông Cổ Chiên.
- Đoạn 2: Đoạn bờ biển từ cửa sông Cổ Chiên đến Cồn Bửng.
- Đoạn 3: Đoạn bờ biển từ Cồn Bửng đến cửa sông Hàm Luông.
- Đoạn 4: Đoạn bờ cửa sông Hàm Luông.
1. Giai đoạn 2006 – 2011: Xu thế xói lở chiếm ưu thế
- Tổng diện tích biến động: Diện tích xói lở đạt $211,87\text{ ha}$, gấp hơn 2 lần diện tích bồi tụ ($99,15\text{ ha}$).
- Đoạn 1: Xói lở trung bình $-72,61\text{ m}$, điểm xói mạnh nhất lên tới $308,07\text{ m}$ (mất $61,55\text{ ha}$). Bồi tụ không đáng kể ($0,08\text{ ha}$).
- Đoạn 2: Bồi tụ chiếm ưu thế với bề rộng bồi trung bình $+55,49\text{ m}$, điểm bồi lớn nhất đạt $609,89\text{ m}$ (tăng $99,07\text{ ha}$).
- Đoạn 3: Xói lở toàn diện, bề rộng xói trung bình $-77,27\text{ m}$, điểm xói lớn nhất $298,06\text{ m}$ (mất $97,72\text{ ha}$). Không ghi nhận bồi tụ.
- Đoạn 4: Xói lở trung bình $-18,46\text{ m}$, xói cực đại $55,38\text{ m}$ (mất $43,05\text{ ha}$).
2. Giai đoạn 2011 – 2015: Xu thế bồi tụ phục hồi mạnh mẽ
- Tổng diện tích biến động: Diện tích bồi tụ đạt $154,78\text{ ha}$, gấp gần 6 lần diện tích xói lở ($26,00\text{ ha}$).
- Đoạn 1: Chuyển dịch hoàn toàn sang bồi tụ, bề rộng bồi trung bình $+129,99\text{ m}$, cực đại $+367,98\text{ m}$ (tăng $62,22\text{ ha}$ đất).
- Đoạn 2: Tiếp tục bồi tụ mạnh với diện tích $77,75\text{ ha}$, bồi trung bình $+115,22\text{ m}$, cực đại $+336,87\text{ m}$.
- Đoạn 3: Giảm tốc độ xói nhưng vẫn mất đất, xói trung bình $-22,79\text{ m}$, diện tích xói $24,10\text{ ha}$.
- Đoạn 4: Xu thế bồi nhẹ chiếm ưu thế với $+10,52\text{ ha}$, bồi trung bình $+26,81\text{ m}$.
3. Toàn giai đoạn 2006 – 2015 và Tốc độ biến động tuyến tính (LRR)
Tổng hợp 9 năm, tổng diện tích bồi tụ toàn huyện đạt $235,84\text{ ha}$, vượt qua tổng diện tích xói lở ($186,19\text{ ha}$).
| Phân đoạn bờ biển |
Thay đổi diện tích Bồi (ha) |
Thay đổi diện tích Xói (ha) |
Tốc độ biến động LRR lớn nhất (m/năm) |
Tốc độ biến động LRR nhỏ nhất (m/năm) |
Tốc độ LRR trung bình (m/năm) |
Hệ số tin cậy $LR^2$ trung bình |
Đánh giá mức độ |
| Đoạn 1 (Cửa Cổ Chiên) |
19,03 |
33,10 |
+39,70 (Bồi) |
-15,07 (Xói) |
-2,76 |
0,34 |
Xói lở trung bình (bồi xói đan xen) |
| Đoạn 2 (Cổ Chiên - Cồn Bửng) |
215,69 |
0,00 |
+9,19 (Bồi) |
-3,91 (Xói) |
+24,95 |
0,76 |
Bồi tụ rất lớn (tăng trưởng ổn định) |
| Đoạn 3 (Cồn Bửng - Hàm Luông) |
0,00 |
120,22 |
0,00 |
-41,18 (Xói) |
-20,63 |
0,88 |
Xói lở nghiêm trọng (tính tuyến tính cao) |
| Đoạn 4 (Cửa Hàm Luông) |
1,13 |
32,87 |
+1,99 (Bồi) |
-4,58 (Xói) |
-1,37 |
0,33 |
Xói lở trung bình (biến động nhẹ) |
| Toàn dải bờ Thạnh Phú |
235,84 |
186,19 |
+39,70 |
-41,18 |
+0,19 net |
0,58 |
Bồi tụ chiếm ưu thế tổng thể |
Đổi mới và đóng góp
- Cải tiến giá trị ngưỡng phổ cho cảm biến thế hệ mới: Phát hiện và xử lý sự khác biệt cấu trúc phổ giữa Landsat 5 TM (7 kênh phổ) và Landsat 8 OLI (11 kênh phổ). Đề tài đã hiệu chỉnh ngưỡng tỷ số $(B6+B7)/B3$ từ mức mặc định $1,0$ lên ngưỡng tối ưu $1,2376$. Cải tiến này giúp loại bỏ triệt để nhiễu phản xạ bọt sóng ven bờ và sai số phân loại giữa bãi bùn ẩm ướt với vùng nước nông ven bờ.
- Chuẩn hóa quy trình giám sát bán tự động: So với phương pháp chỉ số khác biệt thực vật (NDVI) hoặc chỉ số nước khác biệt chuẩn hóa (NDWI) đơn lẻ, thuật toán tỷ số kép Gathot Winasor kết hợp DSAS nâng cao độ chính xác trích xuất ranh giới bờ biển thảm thực vật rừng ngập mặn lên trên $92%$, giảm thời gian xử lý dữ liệu vector đến $85%$ so với số hóa thủ công.
- Cung cấp bức tranh động lực học bờ biển độ phân giải cao: Lần đầu tiên xây dựng được chuỗi dữ liệu định lượng chi tiết theo từng bước cắt $50\text{ m}$ (transect) cho toàn bộ 25 km bờ biển Thạnh Phú, chỉ rõ điểm nóng xói lở nghiêm trọng nhất tại Đoạn 3 (xói lở lên tới $41,18\text{ m/năm}$) và khu vực bồi tụ bãi bồi Cồn Bửng ($24,95\text{ m/năm}$).
Ứng dụng thực tế và triển khai
+-------------------------------------------------------------------------------+
| KHUNG ỨNG DỤNG THỰC TẾ & QUẢN LÝ |
+-------------------------------------------------------------------------------+
|
+---------------------------------+---------------------------------+
| |
+-----------------------------------+ +-----------------------------------+
| QUẢN LÝ LÂM NGHIỆP & PHÒNG HỘ | | QUY HOẠCH KINH TẾ & DU LỊCH |
| - Quy hoạch trồng rừng ngập mặn | | - Bảo vệ khu du lịch sinh thái |
| (Bần chua, Mắm tại vùng xói lở) | | Cồn Bửng (xã Thạnh Hải) |
| - Cố định giồng cát ven biển | | - Kiểm soát hoạt động khai thác |
| (Phi lao, Keo gai tại bãi bồi) | | cát lòng sông Cổ Chiên/Hàm Luông|
+-----------------------------------+ +-----------------------------------+
- Ứng dụng trong lâm nghiệp và phòng hộ đê biển: Kết quả nghiên cứu cung cấp cơ sở không gian chuẩn xác để Sở NN&PTNT tỉnh Bến Tre khoanh vùng phục hồi rừng ngập mặn. Tại Đoạn 3 (tốc độ sạt lở $> 20\text{ m/năm}$), đề xuất ưu tiên trồng cây con cỡ lớn trong bầu có giá đỡ bằng loài Bần chua (Sonneratia caseolaris) hoặc Mắm (Avicennia). Tại các dải cồn cát mới bồi ở Đoạn 2, triển khai trồng Phi lao (Casuarina equisetifolia) để chắn gió bão và cố định đụn cát.
- Quản lý tài nguyên cát và chỉnh trị dòng hàn: Định lượng mức độ sạt lở giúp cơ quan quản lý môi trường siết chặt hoạt động khai thác cát không phép tại hai cửa sông Cổ Chiên và Hàm Luông – một trong những nguyên nhân nhân sinh làm hạ thấp thềm bãi và gia tăng năng lượng sóng phá hủy chân bờ.
- Hiệu quả kinh tế (ROI) và khả năng mở rộng: Chi phí thực hiện một chu kỳ quan trắc bờ biển bằng công nghệ Viễn thám - GIS giảm hơn $90%$ kinh phí so với phương pháp thuê tàu thuyền đo trắc địa hồi âm ven bờ. Giải pháp sẵn sàng mở rộng và nhân rộng cho toàn bộ 65 km bờ biển tỉnh Bến Tre (Bình Đại, Ba Tri, Thạnh Phú) và toàn dải ven biển ĐBSCL.
Hạn chế và hướng phát triển
- Hạn chế kỹ thuật:
- Kích thước pixel ảnh Landsat ($30\text{ m} \times 30\text{ m}$) tạo ra biên độ sai số vị trí $\pm 15\text{ m}$ tại các vùng bờ có độ dốc thấp.
- Chưa tích hợp mô hình hiệu chỉnh mực nước triều tức thời do thiếu số liệu đo đạc thực tế địa hình đáy biển tỷ lệ lớn tại thời điểm vệ tinh quét qua.
- Dữ liệu quan sát phụ thuộc vào mùa khô (tháng 1 - tháng 2) do mùa mưa bị ảnh hưởng lớn bởi mây che phủ.
- Hướng phát triển:
- Ứng dụng nguồn ảnh vệ tinh quang học độ phân giải cao (Sentinel-2A/2B $10\text{ m}$, PlanetScope $3\text{ m}$) kết hợp ảnh Radar khẩu độ tổng hợp (SAR Sentinel-1 C-band) để quan trắc xuyên mây trong cả mùa bão lũ.
- Xây dựng pipeline tự động hóa hoàn toàn trên nền tảng điện toán đám mây Google Earth Engine (GEE) bằng ngôn ngữ Python/JavaScript.
- Kết hợp mô hình thủy động lực học sóng - dòng chảy (SWAN, Delft3D) và mô hình số độ cao hồi âm đáy biển để dự báo vị trí đường bờ giai đoạn 2030 – 2050 theo các kịch bản biến đổi khí hậu (RCP 4.5 và RCP 8.5).
Đối tượng hưởng lợi
- Cơ quan quản lý nhà nước (Sở TN&MT, Chi cục Thủy lợi Bến Tre): Nắm bắt dữ liệu định lượng tin cậy để cấp phép sử dụng bãi bồi ven biển, hoạch định ngân sách đầu tư đê kè chắn sóng và cắm mốc hành lang bảo vệ bờ biển.
- Kỹ sư GIS & Viễn thám: Tiếp cận bộ quy trình kỹ thuật hoàn chỉnh từ xử lý ảnh Landsat TM/OLI, Band Math hiệu chỉnh ngưỡng đến phân tích không gian với DSAS.
- Nhà nghiên cứu & Sinh viên chuyên ngành: Tài liệu tham khảo ứng dụng thực tiễn về động lực học vùng bờ ĐBSCL, cung cấp phương pháp luận kết hợp giữa viễn thám đa phổ và thống kê địa lý.
- Cộng đồng dân cư địa phương: Hưởng lợi gián tiếp từ các công trình đê kè phòng hộ, bảo vệ diện tích nuôi trồng thủy sản và giảm thiểu rủi ro thiên tai sạt lở cuốn trôi nhà cửa.
Câu hỏi thường gặp
1. Yêu cầu hệ thống phần cứng và phần mềm tối thiểu để triển khai giải pháp này là gì?
Hệ thống yêu cầu máy trạm tối thiểu CPU Intel Core i5/AMD Ryzen 5 thế hệ 4 trở lên, 8GB RAM (khuyến nghị 16GB), ổ cứng SSD trống tối thiểu 50GB. Về phần mềm: ArcGIS Desktop 10.x cài đặt kèm phần mở rộng DSAS v4.3, ENVI 4.7 trở lên, và Microsoft Access Database Engine (để DSAS ghi kết quả vào Personal/File Geodatabase).
2. Làm thế nào để loại bỏ sai số do sự khác biệt mực nước thủy triều khi ảnh vệ tinh được chụp ở các ngày khác nhau?
Để hiệu chỉnh triều triệt để, cần thu thập số liệu mực nước giờ từ trạm hải văn gần nhất (như trạm Bình Đại hoặc Vũng Tàu) tại thời điểm chụp ảnh, kết hợp với mô hình số độ cao địa hình bãi triều (DEM). Đường bờ sẽ được hiệu chỉnh dời vị trí theo công thức: $\Delta X = \frac{\Delta H}{\tan(\beta)}$, trong đó $\Delta H$ là chênh lệch mực triều và $\beta$ là góc dốc bãi biển. Khi thiếu DEM, phương pháp chọn ảnh tại cùng một thời điểm mùa khô và giờ chụp tương đồng (14:00 - 15:00) giúp hạn chế tối đa độ lệch này.
3. Tại sao ngưỡng tỷ số (B6+B7)/B3 của Landsat 8 OLI phải đặt là 1,2376 thay vì 1,0 như Landsat 5 TM?
Cảm biến OLI trên Landsat 8 có bước sóng và độ nhạy phổ được thu hẹp, cải tiến so với cảm biến TM để tránh các dải hấp thụ hơi nước khí quyển (Kênh 6 OLI: 1,56 – 1,66 $\mu\text{m}$; Kênh 7 OLI: 2,10 – 2,30 $\mu\text{m}$; Kênh 3 OLI: 0,525 – 0,600 $\mu\text{m}$). Việc nâng ngưỡng lên 1,2376 qua hiệu chuẩn thực nghiệm bảo đảm giá trị phân loại nhị phân đất/nước của Landsat 8 hoàn toàn đồng nhất với không gian phổ của Landsat 5 TM.
4. Chỉ số LRR (Linear Regression Rate) có ưu điểm gì vượt trội so với EPR (End Point Rate)?
Chỉ số EPR chỉ tính toán tốc độ thay đổi dựa trên điểm đầu và điểm cuối thời gian ($EPR = \frac{D_{cuối} - D_{đầu}}{\Delta t}$), do đó hoàn toàn bỏ qua các biến động trung gian và dễ bị bóp méo bởi các sự kiện dị thường cục bộ. Ngược lại, LRR sử dụng phương pháp bình phương tối thiểu để khớp một đường hồi quy tuyến tính qua tất cả các mốc thời gian quan sát, giảm thiểu sai số ngẫu nhiên và cung cấp kèm hệ số $R^2$ để định lượng độ tin cậy của xu thế.
5. Chi phí thu thập dữ liệu ảnh viễn thám thực hiện dự án này là bao nhiêu?
Toàn bộ kho dữ liệu ảnh Landsat 5 TM và Landsat 8 OLI/TIRs được cung cấp hoàn toàn miễn phí bởi Cơ quan Địa chất Hoa Kỳ (USGS) thông qua cổng EarthExplorer hoặc GloVis. Chi phí đầu tư dự án chỉ tập trung vào nhân lực xử lý và bản quyền phần mềm phân tích chuyên dụng.
Kết luận
Đồ án tốt nghiệp đã giải quyết trọn vẹn bài toán quan trắc và phân tích động lực học dải bờ biển huyện Thạnh Phú, tỉnh Bến Tre trong giai đoạn 2006 – 2015. Bằng việc kết hợp thành công thuật toán tỷ số kênh phổ Gathot Winasor cải tiến trên ENVI và công cụ thống kê không gian DSAS trên ArcGIS, đề tài đã trích xuất chuẩn xác vị trí đường bờ và định lượng chi tiết tốc độ bồi - xói dọc theo 25 km bờ biển.
Kết quả nghiên cứu khẳng định: Mặc dù toàn dải bờ Thạnh Phú có xu thế bồi tụ nhẹ về mặt diện tích tổng thể ($+235,84\text{ ha}$ bồi so với $-186,19\text{ ha}$ xói), quá trình sạt lở diễn ra đặc biệt khốc liệt và mang tính tuyến tính cao tại khu vực từ Cồn Bửng đến cửa sông Hàm Luông (tốc độ sạt lở trung bình $-20,63\text{ m/năm}$, cực đại $-41,18\text{ m/năm}$). Những phát hiện này là cơ sở khoa học cốt lõi phục vụ công tác quy hoạch trồng rừng ngập mặn phòng hộ, quản lý khai thác tài nguyên và bảo vệ an toàn sinh kế bền vững cho cộng đồng dân cư ven biển Bến Tre trước thách thức biến đổi khí hậu toàn cầu.
