Tổng quan nghiên cứu

Việt Nam sở hữu mạng lưới thủy văn phong phú nhưng chịu sự phân hóa sâu sắc theo không gian và thời gian với lượng mưa trung bình cả nước đạt khoảng 1.940 mm/năm và địa hình đồi núi chiếm tới 75% diện tích lãnh thổ. Tại tỉnh Quảng Ninh – địa bàn trọng điểm kinh tế phía Bắc với tốc độ tăng trưởng kinh tế duy trì ở mức khoảng 12%/năm, áp lực từ quá trình công nghiệp hóa, đô thị hóa và khai thác khoáng sản than đã tạo nên những tác động tiêu cực rõ nét đến hệ thống sinh thái thủy vực. Tình trạng bồi lấp, chuyển đổi mặt nước và suy thoái nguồn nước tự nhiên đang diễn biến phức tạp trên hệ thống sông suối có mật độ từ 1,0 đến 2,4 km/km².

Đứng trước thực trạng trên, đề tài tập trung giải quyết bài toán định lượng và đánh giá quy luật biến động không gian của các vùng mặt nước có quy mô lớn trên 50 km² tại tỉnh Quảng Ninh. Mục tiêu cụ thể của công trình bao gồm: xây dựng bộ bản đồ hiện trạng mặt nước qua các giai đoạn 5 năm từ năm 1990 đến năm 2015, phân tích nguyên nhân tác động từ các hoạt động nhân sinh và yếu tố tự nhiên, dự báo xu thế biến động tài nguyên nước trong giai đoạn 10 năm tiếp theo đến năm 2025, đồng thời thiết lập hệ thống giải pháp kỹ thuật và quản lý bền vững.

Phạm vi nghiên cứu bao quát toàn bộ địa bàn tỉnh Quảng Ninh trong khung thời gian 25 năm (1990 - 2015), sử dụng chuỗi dữ liệu ảnh vệ tinh đa thời gian. Nghiên cứu mang ý nghĩa thực tiễn quan trọng, giúp tối ưu hóa độ chính xác giải đoán hiện trạng đạt trên 85%, giảm thiểu hơn 60% chi phí và thời gian so với phương pháp đo đạc truyền thống ngoài thực địa, đồng thời cung cấp luận cứ khoa học vững chắc phục vụ quy hoạch không gian kinh tế - xã hội gắn liền với bảo vệ tài nguyên nước.

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 dựa trên sự tích hợp giữa lý thuyết lan truyền bức xạ điện từ trong viễn thám quang học và lý thuyết phân tích không gian địa lý đa tầng của hệ thống thông tin địa lý. Theo nguyên lý vật lý viễn thám, các đối tượng tự nhiên trên bề mặt Trái Đất có đặc trưng phản xạ và bức xạ sóng điện từ hoàn toàn khác biệt tại từng bước sóng quang phổ, từ dải ánh sáng nhìn thấy (0,4 - 0,7 µm) đến dải hồng ngoại gần và hồng ngoại sóng ngắn. Mặt nước tự nhiên có đặc tính hấp thụ mạnh năng lượng bức xạ ở dải sóng hồng ngoại và phản xạ chọn lọc ở dải khả kiến, tạo nên cơ sở khoa học để thiết lập các chỉ số tách chiết thủy vực và xây dựng khóa giải đoán phổ.

Song song với đó, khung mô hình tích hợp dữ liệu địa không gian sử dụng các khái niệm nền tảng bao gồm: Hệ thống thông tin địa lý (GIS) với cấu trúc dữ liệu không gian tích hợp dữ liệu thuộc tính; Công nghệ viễn thám (Remote Sensing) khai thác dữ liệu vệ tinh đa phổ; Phương pháp phân loại ảnh có kiểm định (Supervised Classification) dựa trên thuật toán xác suất cực đại; và Khóa giải đoán ảnh (Interpretation Keys) gồm 8 chuẩn nhận biết trực quan về tông màu, hình dạng, kích thước, cấu trúc, hoa văn, độ bóng, độ đậm nhạt và mối liên hệ tương quan địa lý.

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

Nguồn dữ liệu sơ cấp và thứ cấp của đề tài bao gồm chuỗi ảnh vệ tinh Landsat 5 TM và Landsat 8 OLI/TIRS thu nhận tại 6 thời điểm then chốt: năm 1990, 1995, 2000, 2005, 2010 và 2015. Cỡ mẫu nghiên cứu bao gồm 6 cảnh ảnh vệ tinh toàn tỉnh với độ phân giải không gian gốc 30m, kết hợp dữ liệu đối chứng gồm 85 điểm khống chế mặt đất (GCPs) thu thập trên toàn bộ các khu vực cửa sông, hồ chứa lớn và vùng đầm phá ven biển. Phương pháp chọn mẫu là chọn mẫu phân tầng ngẫu nhiên kết hợp khảo sát thực địa nhằm đảm bảo tính đại diện cho cả vùng đồi núi cánh cung Đông Triều (chiếm 60,5% diện tích tự nhiên) và vùng trung du ven biển (chiếm 10% diện tích).

Quy trình xử lý dữ liệu được tiến hành nghiêm ngặt qua các giai đoạn kỹ thuật:

  • Hiệu chỉnh bức xạ và loại trừ ảnh hưởng khí quyển bằng các mô hình truyền xạ quang học, khắc phục hiện tượng mất dòng ảnh và nhiễu dải quét.
  • Hiệu chỉnh hình học phi hệ thống dựa trên đa thức bậc hai với các điểm khống chế mặt đất, đưa sai số vị trí hình học xuống dưới mức 1,0 pixel.
  • Tăng cường chất lượng ảnh bằng phương pháp PanSharpening, kết hợp kênh toàn sắc panchromatic độ phân giải cao để nâng độ chi tiết ảnh từ 30m lên 15m.
  • Phân loại có kiểm định bằng thuật toán Maximum Likelihood trên phần mềm chuyên dụng ENVI nhằm tối ưu hóa khả năng phân tách các lớp phủ mặt nước, đất trống, rừng và khu dân cư.
  • Phân tích không gian chồng lớp bằng các công cụ Intersect, Mosaic To New Raster và Extract by Mask trên phần mềm ArcGIS 10 để chiết xuất diện tích và lập bản đồ biến động.

Phương pháp phân loại có kiểm định Maximum Likelihood được lựa chọn nhờ khả năng tính toán hàm mật độ xác suất của từng pixel thuộc về một lớp đối tượng cụ thể, đem lại độ chính xác tin cậy cao nhất khi xử lý các vùng thủy vực có độ đục và hàm lượng phù sa biến động phức tạp ven bờ biển Quảng Ninh.

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

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

Phân tích chuỗi dữ liệu đa thời gian giai đoạn 1990 - 2015 đã làm rõ bức tranh biến động phức tạp của diện tích mặt nước tỉnh Quảng Ninh:

  • Biến động tổng diện tích mặt nước tự nhiên: Trong giai đoạn 1990 - 2005, diện tích mặt nước tự nhiên tại các vùng cửa sông và vũng vịnh ven bờ có xu hướng suy giảm trung bình khoảng 8,5% do áp lực san lấp mặt bằng, quai đê lấn biển phát triển đô thị tại thành phố Hạ Long, thành phố Cẩm Phả và thị xã Quảng Yên.
  • Sự gia tăng của diện tích mặt nước nhân tạo: Giai đoạn 2005 - 2015 ghi nhận sự mở rộng nhanh chóng của hệ thống mặt nước nhân tạo phục vụ nuôi trồng thủy sản và tích trữ nước ngọt, với mức tăng khoảng 22,4% tại khu vực ven biển Tiên Yên, Đầm Hà và Hải Hà.
  • Phân hóa không gian biến động mạnh mẽ: Vùng miền Tây (vùng núi Đông Triều, Uông Bí, Hoành Bồ) chứng kiến sự thay đổi chế độ dòng chảy mặt và bồi lắng lòng hồ do hoạt động khai thác than lộ thiên; trong khi vùng hải đảo (huyện Vân Đồn và Cô Tô với hơn 2.000 hòn đảo) diện tích mặt nước biến đổi cục bộ dưới 5% và chịu ảnh hưởng chủ yếu bởi biên độ triều ven bờ.
  • Độ tin cậy của mô hình viễn thám: Kết quả đánh giá độ chính xác tổng thể của bản đồ giải đoán năm 2015 đạt 88,6%, hệ số tương quan Kappa đạt 0,85, khẳng định chất lượng xuất sắc của quy trình xử lý dữ liệu viễn thám Landsat.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự chuyển đổi diện tích mặt nước bắt nguồn từ mối quan hệ tương tác giữa các yếu tố tự nhiên và phát triển kinh tế. Lượng mưa tại Quảng Ninh tập trung cao độ vào mùa hè (chiếm từ 75% đến 85% tổng lượng mưa cả năm từ 1.600 đến 2.700 mm), kết hợp độ dốc địa hình lớn khiến các con sông ngắn và dốc dễ xuất hiện lũ quét làm thay đổi lòng dẫn tự nhiên. Về mặt nhân sinh, quá trình đô thị hóa mạnh mẽ cùng tốc độ tăng trưởng kinh tế bình quân 12%/năm đã thúc đẩy các dự án chuyển đổi bãi bồi, đất ngập triều thành đất công nghiệp và khu dân cư tập trung. Hoạt động đổ thải đất đá từ các mỏ than tại Cẩm Phả và Uông Bí đã bồi lấp hàng loạt khe suối tự nhiên, làm suy giảm dung tích tích nước của các hồ chứa thượng nguồn.

So sánh với các nghiên cứu tương tự tại vùng châu thổ sông Hồng, tốc độ suy giảm mặt nước tự nhiên tại Quảng Ninh mang tính chất cục bộ gắn liền với công nghiệp khai khoáng và lấn biển, thay vì áp lực nông nghiệp thuần túy.

Trong báo cáo khoa học, dữ liệu nghiên cứu được trực quan hóa tối ưu thông qua:

  • Biểu đồ cột biểu diễn biến thiên diện tích các loại hình mặt nước qua từng chu kỳ 5 năm (1990, 1995, 2000, 2005, 2010, 2015).
  • Bảng ma trận chuyển đổi sử dụng đất phản ánh tỷ lệ diện tích chuyển hóa giữa 5 trạng thái: mặt nước, đất rừng, đất trống, khu dân cư và đất bãi triều.
  • Hệ thống bản đồ chuyên đề chồng lớp GIS tỷ lệ 1:100.000 thể hiện rõ các điểm nóng suy giảm mặt nước tại khu vực vịnh Cửa Lục, ven biển Quảng Yên và các hồ đập thủy lợi lớn.

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

Nhằm bảo tồn và nâng cao hiệu quả quản trị tài nguyên nước tỉnh Quảng Ninh trong bối cảnh biến đổi khí hậu, đề tài đưa ra 4 nhóm giải pháp trọng tâm:

  1. Thiết lập và vận hành hệ thống giám sát viễn thám định kỳ: Sở Tài nguyên và Môi trường tỉnh Quảng Ninh chủ trì triển khai cập nhật ảnh viễn thám quang học và radar với chu kỳ 6 tháng/lần, hoàn thành tích hợp cổng thông tin CSDL WebGIS quản lý tài nguyên nước mặt trước quý IV năm 2024, duy trì độ chính xác giám sát diện tích đạt trên 90%.
  2. Quy hoạch và cắm mốc hành lang bảo vệ nguồn nước mặt: Ủy ban nhân dân tỉnh Quảng Ninh chỉ đạo Sở Nông nghiệp và Phát triển Nông thôn tiến hành cắm mốc ranh giới bảo vệ cho 100% các hồ chứa nước ngọt có dung tích trên 50 triệu m³ (như hồ Yên Lập, hồ Cao Vân), phấn đấu phục hồi tối thiểu 15% diện tích mặt nước tự nhiên bị bồi lắng vào năm 2026.
  3. Kiểm soát dòng thải và phục hồi môi trường lưu vực khai thác khoáng sản: Tập đoàn Công nghiệp Than - Khoáng sản Việt Nam phối hợp chính quyền địa phương hoàn nguyên môi trường 100% các bãi thải mỏ lộ thiên xung quanh lưu vực sông suối, giảm thiểu 30% lượng bùn cát bồi lắng vào các thủy vực trước năm 2028.
  4. Tăng cường năng lực quản trị tài nguyên nước dựa vào cộng đồng: Ủy ban nhân dân các huyện, thị xã, thành phố (Hạ Long, Cẩm Phả, Đông Triều, Quảng Yên) tổ chức các chiến dịch truyền thông, tập huấn phổ biến Luật Tài nguyên nước, hướng tới mục tiêu hơn 80% người dân sinh sống ven các lưu vực sông lớn tham gia giám sát các hành vi xâm lấn mặt nước trái phép trước năm 2025.

Đố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 giá trị ứng dụng thực tiễn cao cho nhiều nhóm đối tượng:

  • Cơ quan quản lý nhà nước về tài nguyên và môi trường (Sở TN&MT, Sở NN&PTNT tỉnh Quảng Ninh và các tỉnh ven biển): Sử dụng làm tài liệu quy chuẩn để áp dụng quy trình công nghệ viễn thám - GIS vào công tác kiểm kê, thanh tra, theo dõi biến động diện tích mặt nước và cấp phép khai thác nguồn nước định kỳ.
  • Chuyên gia và đơn vị tư vấn quy hoạch không gian đô thị: Khai thác dữ liệu lịch sử biến động mặt nước 25 năm nhằm thiết lập các kịch bản thoát lũ, thiết kế hệ thống hồ điều hòa sinh thái và tối ưu hóa quy hoạch hạ tầng ven biển thích ứng với nước biển dâng.
  • Giảng viên, nhà nghiên cứu và học viên cao học chuyên ngành Quản lý Môi trường, Trắc địa - Viễn thám: Sử dụng công trình như một tài liệu tham khảo mẫu mực về phương pháp xử lý chuỗi ảnh Landsat (hiệu chỉnh khí quyển, tăng cường phân giải PanSharpening, phân loại Maximum Likelihood) và phân tích không gian trong ArcGIS.
  • Doanh nghiệp khai thác khoáng sản và nuôi trồng thủy sản: Vận dụng các đánh giá về tác động thủy văn để xây dựng báo cáo đánh giá tác động môi trường (ĐTM), thiết kế hệ thống hồ lắng xử lý nước thải mỏ đạt tiêu chuẩn và tối ưu hóa diện tích mặt nước đầm nuôi an toàn.

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

1. Tại sao công nghệ viễn thám và GIS là giải pháp tối ưu để giám sát diện tích mặt nước?
Công nghệ viễn thám cung cấp chuỗi dữ liệu ảnh vệ tinh chụp lặp lại theo chu kỳ với góc nhìn đồng bộ trên phạm vi rộng lớn, giúp tái hiện lịch sử biến động mặt nước 25 năm mà phương pháp đo đạc thực địa không thể thực hiện được. Kết hợp với GIS, các phép phân tích không gian chồng lớp giúp định lượng chính xác diện tích chuyển đổi, tiết kiệm hơn 60% chi phí điều tra.

2. Quy trình tiền xử lý ảnh Landsat cần lưu ý những khâu kỹ thuật then chốt nào?
Nghiên cứu cần thực hiện nghiêm ngặt 3 khâu: hiệu chỉnh bức xạ để loại trừ hiện tượng mất dòng ảnh và góc chiếu mặt trời; hiệu chỉnh khí quyển bằng hàm truyền bức xạ; và nắn chỉnh hình học phi hệ thống dựa trên đa thức bậc hai với các điểm khống chế mặt đất (GCPs) nhằm bảo đảm sai số vị trí dưới 1,0 pixel.

3. Hoạt động khai thác than tại Quảng Ninh đã ảnh hưởng cụ thể như thế nào đến tài nguyên nước mặt?
Khai thác than lộ thiên tạo ra các bãi thải đất đá lớn, khi gặp lượng mưa mùa hè chiếm 75 - 85% năm sẽ sinh ra dòng chảy tràn mang theo lượng bùn cát khổng lồ. Quá trình này gây bồi lấp nghiêm trọng lòng dẫn sông suối, làm suy giảm dung tích hữu ích của các hồ chứa tự nhiên và nhân tạo tại khu vực Uông Bí và Cẩm Phả.

4. Phương pháp nào được sử dụng để nâng cao độ chính xác giải đoán đối tượng mặt nước?
Nghiên cứu đã áp dụng kỹ thuật PanSharpening để kết hợp kênh toàn sắc Panchromatic và kênh đa phổ, nâng độ phân giải không gian từ 30m lên 15m. Đồng thời, việc ứng dụng thuật toán Maximum Likelihood kết hợp xây dựng khóa giải đoán 8 tiêu chuẩn đã giúp nâng độ chính xác toàn diện của bản đồ lên trên 88%.

5. Kết quả dự báo biến động mặt nước đến năm 2025 mang lại giá trị gì cho địa phương?
Dự báo xu thế 10 năm chỉ ra nguy cơ thu hẹp mặt nước tự nhiên tại các vùng đô thị hóa nhanh và cảnh báo thiếu hụt nước ngọt cục bộ vào mùa khô. Đây là căn cứ khoa học quyết định để tỉnh Quảng Ninh chủ động điều chỉnh quy hoạch sử dụng đất, bảo vệ nghiêm ngặt các vùng sinh thủy và xây dựng đập trữ nước ngọt chiến lược.

Kết luận

  • Luận văn đã thiết lập thành công quy trình ứng dụng công nghệ viễn thám đa phổ Landsat và GIS trong việc đánh giá định lượng biến động diện tích mặt nước tại tỉnh Quảng Ninh qua 6 mốc thời gian từ năm 1990 đến năm 2015.
  • Nghiên cứu đã số hóa và xây dựng bộ bản đồ hiện trạng diện tích mặt nước với độ chính xác phân loại đạt trên 88%, sai số hình học dưới 1,0 pixel, khẳng định tính khả thi vượt trội của phương pháp địa không gian.
  • Đề tài đã chỉ rõ quy luật chuyển dịch diện tích mặt nước, xác định áp lực từ hoạt động khai khoáng, san lấp phát triển đô thị và biến đổi lượng mưa là các nguyên nhân cốt lõi dẫn đến suy thoái thủy vực.
  • Hệ sinh thái 4 nhóm giải pháp kỹ thuật, quy hoạch, công nghệ và quản trị cộng đồng đã được đề xuất với lộ trình thực thi rõ ràng từ năm 2024 đến năm 2028 nhằm đảm bảo an ninh nguồn nước bền vững.
  • Đóng góp khoa học then chốt của công trình là cung cấp cơ sở dữ liệu số hóa hoàn chỉnh và phương pháp luận chuẩn xác phục vụ các cơ quan quản lý nhà nước ban hành chính sách bảo vệ tài nguyên nước mặt kịp thời.