Chào mừng các bạn tới với Chặng 3 trong Series "Lộ trình làm chủ GIS từ số 0"!
Nếu là người Việt Nam, có lẽ ai trong chúng ta cũng đã nhiều lần chứng kiến những thiệt hại nặng nề về người và tài sản sau mỗi mùa mưa bão. Và rồi không ít lần tự hỏi: Tại sao cùng một lượng mưa trút xuống, nhưng ngã tư này lại nhanh chóng biến thành "dòng sông", trong khi con phố chỉ cách đó vài trăm mét vẫn gần như khô ráo? Hay vì sao những trận lũ quét ở vùng núi lại dồn toàn bộ sức mạnh vào đúng một thung lũng, gây nên sức tàn phá khủng khiếp?
Câu trả lời nằm ở một yếu tố tưởng chừng đơn giản nhưng lại quyết định gần như toàn bộ hành trình của nước trên bề mặt Trái Đất: hình thái địa hình.
Trong hai bài trước của series Phân tích bề mặt địa hình trong GIS, chúng ta đã tìm hiểu cách mô hình số độ cao (DEM – Digital Elevation Model) biểu diễn địa hình và cách khai thác các đặc trưng như độ dốc, hướng dốc hay địa hình bóng đổ. Tuy nhiên, DEM không chỉ là một mô hình mô tả địa hình. Khi được xử lý bằng các thuật toán thủy văn, nó còn cho phép chúng ta mô phỏng cách nước di chuyển trên bề mặt, từ những giọt mưa đầu tiên cho đến khi chúng hội tụ thành các con suối, dòng sông và cuối cùng thoát ra khỏi lưu vực.
Đó chính là phân tích thủy văn từ DEM (Surface Hydrology Analysis) – tập hợp các kỹ thuật giúp xác định hướng dòng chảy, tính toán mức độ tích tụ dòng chảy, xây dựng mạng lưới sông suối và phân chia lưu vực. Đây là nền tảng của nhiều ứng dụng quan trọng như dự báo ngập lụt, đánh giá nguy cơ lũ quét, quy hoạch hệ thống thoát nước đô thị, quản lý tài nguyên nước và quy hoạch hạ tầng.
Trong bài viết thứ ba của series này, chúng ta sẽ cùng khám phá toàn bộ quy trình phân tích thủy văn từ DEM, từ việc xử lý các vùng trũng (Sink Filling), xác định hướng dòng chảy (Flow Direction), tính tích tụ dòng chảy (Flow Accumulation) cho đến phân chia lưu vực (Watershed Delineation). Sau khi hoàn thành, bạn sẽ hiểu cách một mô hình độ cao tưởng chừng chỉ chứa các giá trị cao độ lại có thể được chuyển hóa thành mô hình mô phỏng toàn bộ quá trình nước vận động trên bề mặt địa hình.
 |
| Nguồn: Victor Nguyễn |
DEM và phân tích thủy văn: Nước có "ngôn ngữ" riêng của nó
Có một quy luật gần như bất biến trong tự nhiên: nước luôn chảy từ nơi có thế năng cao xuống nơi có thế năng thấp. Chính quy luật đơn giản này là nền tảng của hầu hết các thuật toán phân tích thủy văn trong GIS.
Về bản chất, mô hình số độ cao (DEM – Digital Elevation Model) là một lưới raster, trong đó mỗi ô (cell) lưu trữ một giá trị cao độ (z). Khi xem xét mối quan hệ giữa một ô và các ô lân cận, hệ thống GIS có thể tính toán chênh lệch độ cao để xác định hướng mà nước có xu hướng chảy.
Từ những phép tính tưởng chừng đơn giản đó, hàng loạt thông tin thủy văn quan trọng có thể được suy ra:
- Dòng chảy sẽ di chuyển theo hướng nào (Flow Direction)?
- Nước sẽ hội tụ ở những vị trí nào (Flow Accumulation)?
- Mạng lưới sông, suối tự nhiên hình thành ra sao (Stream Network)?
- Ranh giới giữa các lưu vực được xác định ở đâu (Watershed Delineation)?
Nói cách khác, DEM chính là "ngôn ngữ" mà GIS sử dụng để mô phỏng sự vận động của nước trên bề mặt địa hình. Chỉ cần mô hình độ cao đủ chính xác và được tiền xử lý phù hợp, chúng ta có thể tái hiện khá sát cách nước di chuyển trong thực tế.
Đó cũng là lý do vì sao chất lượng của DEM đóng vai trò quyết định đối với kết quả phân tích thủy văn. Một DEM được hiệu chỉnh tốt sẽ là nền tảng cho các bài toán dự báo ngập lụt, đánh giá nguy cơ lũ quét, thiết kế hệ thống thoát nước và quản lý tài nguyên nước.
Quy trình 4 bước cốt lõi trong phân tích thủy văn từ DEM
Để chuyển từ một mô hình số độ cao (DEM) thành các sản phẩm như mạng lưới sông suối hay bản đồ lưu vực, dữ liệu thường phải trải qua một quy trình gồm bốn bước chính. Mỗi bước đều sử dụng kết quả của bước trước, tạo thành một pipeline phân tích liên hoàn.
DEM → Find Sink → Flow Direction → Flow Accumulation → Watershed
Bước 1: Xác định và xử lý vùng trũng (Find Sink)
Trước khi mô phỏng dòng chảy, DEM cần được kiểm tra và xác định các vùng trũng (Sink hoặc Depression).
Sink là những ô raster có cao độ thấp hơn các ô xung quanh, khiến nước mô phỏng bị "mắc kẹt" và không thể tiếp tục chảy xuống hạ lưu.
Sink có thể xuất hiện vì hai nguyên nhân:
- Địa hình thực tế, chẳng hạn hồ tự nhiên, hố sụt hoặc các vùng trũng khép kín.
- Sai số dữ liệu, phát sinh trong quá trình đo đạc LiDAR, nội suy DEM hoặc xử lý ảnh viễn thám. Đây là trường hợp phổ biến hơn.
Trong hầu hết các bài toán phân tích dòng chảy bề mặt, mục tiêu là loại bỏ các sink do sai số. Công cụ Fill sẽ nâng cao độ của các ô trong vùng trũng lên mức vừa đủ để tạo ra một đường thoát nước liên tục.
Sau bước này, DEM trở thành một bề mặt có khả năng thoát nước hợp lý, sẵn sàng cho các phép phân tích tiếp theo.
 |
| Nguồn: Victor Nguyễn |
Bước 2: Xác định hướng dòng chảy (Flow Direction)
Sau khi DEM đã được hiệu chỉnh, bước tiếp theo là xác định nước sẽ chảy từ mỗi ô raster sang ô lân cận nào.
Thuật toán được sử dụng phổ biến nhất là D8 (Eight-direction). Với mỗi ô raster, thuật toán sẽ tính độ dốc đến tám ô xung quanh và giả định dòng chảy sẽ đi theo hướng có độ dốc lớn nhất.
Kết quả là một raster Flow Direction, trong đó mỗi ô lưu trữ một mã hướng dòng chảy (ví dụ: Đông, Đông Nam, Nam...). Raster này đóng vai trò là "bản đồ dẫn đường" cho toàn bộ các phép phân tích thủy văn phía sau.
 |
| Nguồn: Victor Nguyễn |
Bước 3: Tính tích tụ dòng chảy (Flow Accumulation)
Nếu Flow Direction cho biết nước chảy đi đâu, thì Flow Accumulation trả lời câu hỏi có bao nhiêu dòng chảy hội tụ tại mỗi vị trí.
Thuật toán sẽ tính số lượng ô raster ở thượng nguồn có dòng chảy đi qua từng ô trên DEM.
Thông thường:
- Các ô nằm trên đỉnh hoặc sườn núi có giá trị tích tụ rất nhỏ.
- Các ô nằm trong khe tụ nước hoặc dọc theo lòng suối có giá trị tích tụ ngày càng lớn vì nhận nước từ nhiều khu vực phía thượng nguồn.
Giá trị tích tụ càng lớn thì vị trí đó càng có xu hướng trở thành dòng chảy tự nhiên.
Bằng cách lựa chọn một ngưỡng (Threshold), chẳng hạn chỉ giữ lại các ô có giá trị tích tụ lớn hơn một giá trị nhất định, chúng ta có thể trích xuất mạng lưới sông suối từ DEM. Việc lựa chọn ngưỡng phụ thuộc vào độ phân giải DEM và quy mô lưu vực nghiên cứu.
Bước 4: Phân chia lưu vực (Watershed Delineation)
Bước cuối cùng là xác định phạm vi của từng lưu vực (Watershed hoặc Basin).
Lưu vực là toàn bộ diện tích mà nước mưa rơi xuống sẽ tập trung và thoát về cùng một điểm thoát nước (Outlet hoặc Pour Point), chẳng hạn cửa sông, vị trí hồ chứa, đập thủy điện hoặc trạm thủy văn.
Dựa trên raster Flow Direction và vị trí của Outlet, GIS sẽ truy vết ngược dòng để xác định toàn bộ khu vực đóng góp dòng chảy cho điểm đó, từ đó tạo ra ranh giới lưu vực.
Thông tin này là nền tảng cho nhiều ứng dụng như quản lý tài nguyên nước, thiết kế hồ chứa, đánh giá nguy cơ lũ quét, quy hoạch thoát nước đô thị và mô hình hóa thủy văn.
 |
| Nguồn: Victor Nguyễn |
Ứng dụng thực tế: Từ màn hình máy tính đến bài toán thực tiễn
Phân tích thủy văn từ DEM không chỉ là một quy trình xử lý dữ liệu trong GIS mà còn là nền tảng cho nhiều bài toán quan trọng trong quản lý tài nguyên, quy hoạch và giảm nhẹ thiên tai.
Phòng chống thiên tai
Thông qua việc xác định hướng dòng chảy, vùng hội tụ nước và ranh giới lưu vực, các nhà chuyên môn có thể nhận diện những khu vực có nguy cơ tập trung dòng chảy lớn sau mưa, hỗ trợ đánh giá nguy cơ ngập lụt, lũ quét và sạt lở đất. Đây cũng là dữ liệu đầu vào quan trọng cho các mô hình dự báo thủy văn và cảnh báo thiên tai.
Quy hoạch đô thị và thiết kế san nền
Trong quá trình quy hoạch khu đô thị, khu công nghiệp hoặc khu dân cư mới, phân tích thủy văn từ DEM giúp các kỹ sư hiểu được hướng thoát nước tự nhiên trước khi tiến hành san nền. Nhờ đó, phương án thiết kế có thể hạn chế việc cắt đứt dòng chảy hoặc tạo ra các vùng trũng nhân tạo gây ngập úng sau mưa.
Kết quả phân tích còn là cơ sở để xác định cao độ thiết kế, bố trí hồ điều hòa, kênh thoát nước và hệ thống cống, giúp cân bằng giữa khối lượng đào đắp và khả năng tiêu thoát nước của toàn khu vực.
Quản lý tài nguyên nước và thủy lợi
Kết quả phân chia lưu vực và xác định mạng lưới sông suối hỗ trợ tính toán diện tích góp nước cho hồ chứa, quy hoạch công trình thủy lợi, quản lý nguồn nước theo lưu vực sông và cung cấp dữ liệu đầu vào cho các mô hình mô phỏng dòng chảy phục vụ vận hành hồ chứa.
Thiết kế hạ tầng giao thông
Trong quy hoạch và thiết kế đường bộ, đường sắt hay cầu cống, phân tích thủy văn từ DEM giúp xác định các tuyến thoát nước tự nhiên và khu vực tập trung dòng chảy. Những thông tin này là cơ sở để lựa chọn vị trí công trình, bố trí cống thoát nước và kết hợp với các tính toán thủy văn nhằm xác định kích thước công trình phù hợp, góp phần giảm nguy cơ hư hỏng do mưa lũ.
Khai thác mỏ và các dự án hạ tầng quy mô lớn
Đối với các dự án khai thác mỏ, bãi thải, sân bay, cảng biển hay khu công nghiệp, việc thay đổi địa hình có thể làm biến đổi hoàn toàn chế độ dòng chảy tự nhiên. Phân tích thủy văn từ DEM giúp đánh giá tác động của các phương án san gạt địa hình, thiết kế hệ thống thu gom và thoát nước mặt, đồng thời giảm nguy cơ xói mòn, ngập úng và ảnh hưởng đến các khu vực xung quanh sau khi dự án đi vào hoạt động.
Lời kết
Qua ba bài viết trong series Phân tích bề mặt địa hình trong GIS, chúng ta đã cùng đi từ những khái niệm nền tảng về mô hình số độ cao (DEM), khai thác các đặc trưng địa hình như độ dốc, hướng dốc, địa hình bóng đổ, cho đến quy trình phân tích thủy văn với các bước Fill Sink→ Flow Direction → Flow Accumulation → Watershed.
Có thể thấy, DEM không đơn thuần là một tập hợp các giá trị cao độ. Khi được kết hợp với các thuật toán phân tích không gian, nó trở thành nguồn dữ liệu quan trọng giúp chúng ta hiểu rõ hình thái địa hình, mô phỏng sự vận động của dòng chảy và hỗ trợ ra quyết định trong nhiều lĩnh vực như quy hoạch đô thị, thiết kế san nền, quản lý tài nguyên nước, phát triển hạ tầng và phòng chống thiên tai.
Hy vọng series này đã giúp bạn có được một cái nhìn hệ thống về phân tích bề mặt địa hình trong GIS, từ những khái niệm cơ bản đến các quy trình xử lý thường gặp trong thực tế. Đây cũng là nền tảng để tiếp cận các chủ đề chuyên sâu hơn như trích xuất mạng lưới sông suối, phân cấp thứ bậc sông ngòi, mô hình thủy văn, phân tích lưu vực hay các ứng dụng của GIS trong quy hoạch và quản lý tài nguyên.
Cảm ơn bạn đã đồng hành cùng mình trong suốt series. Hẹn gặp lại ở những bài viết tiếp theo, nơi chúng ta sẽ tiếp tục khám phá những kỹ thuật và ứng dụng thú vị khác của GIS trong việc giải quyết các bài toán của thế giới thực.
0 Nhận xét