Dòng chảy lũ khác gì so với hiện tượng ngập lụt?

Dòng chảy lũ là đại lượng động lực học thủy văn phản ánh khối lượng và tốc độ vận chuyển của nước trên lưu vực sau mưa bão. Ngược lại, ngập lụt là hiện tượng thủy lực xảy ra khi lưu lượng dòng chảy lũ vượt quá khả năng thoát nước của lòng sông hoặc hệ thống thoát nước, dẫn tới nước tràn bờ và nhấn chìm các vùng đất khô ráo xung quanh.

Những yếu tố nào quyết định lưu lượng đỉnh của dòng chảy lũ?

Lưu lượng đỉnh của dòng chảy lũ chịu chi phối bởi cường độ và thời gian mưa, diện tích và độ dốc lưu vực, dung lượng thấm của đất, và tỷ lệ che phủ thảm thực vật. Các lưu vực dốc, đất mỏng hoặc bề mặt bị bê tông hóa do đô thị hóa làm giảm khả năng thấm nước, khiến thời gian tập trung dòng chảy rút ngắn và lưu lượng đỉnh dâng cao đột ngột.

Phương pháp lý luận tính dòng chảy lũ áp dụng cho loại lưu vực nào?

Phương pháp lý luận chỉ được áp dụng tin cậy cho các lưu vực nhỏ có diện tích dưới 25 kilômét vuông, phổ biến trong thiết kế cống rãnh đô thị và thoát nước đường giao thông. Phương pháp này dựa trên giả thiết mưa phân bố đồng đều theo không gian và thời gian kéo dài ít nhất bằng thời gian tập trung dòng chảy của toàn bộ lưu vực.

Cơ chế dòng chảy tràn Hortonian khác gì so với dòng chảy tràn Dunne?

Dòng chảy tràn Hortonian xảy ra khi cường độ mưa tức thời vượt quá tốc độ thấm của bề mặt đất, thường gặp ở đất trơ trọi, vùng bán khô hạn hoặc khu vực đô thị hóa. Trong khi đó, dòng chảy tràn Dunne xảy ra khi toàn bộ tầng đất đã ngậm no nước từ đáy lên khiến mực nước ngầm dâng sát mặt đất, phổ biến tại các vùng thung lũng ẩm ướt và rừng nhiệt đới có tầng đất dày.

Dòng chảy lũ (flood runoff) là gì?

Dòng chảy lũ là dạng vận động tập trung quy mô lớn của khối lượng nước mặt và nước sát mặt trên lưu vực thủy văn trong và ngay sau các đợt mưa lớn hoặc tuyết tan, biểu hiện bằng sự gia tăng đột ngột của lưu lượng và mực nước sông suối vượt qua trạng thái dòng chảy kiệt cơ bản. Khái niệm này bao gồm toàn bộ quá trình động lực học từ khi các giọt mưa chuyển hóa thành dòng chảy tràn sườn dốc, tụ thủy vào mạng lưới sông nhánh, truyền sóng lũ dọc lòng dẫn chính và thoát ra cửa sông hoặc vùng hạ lưu. Bài viết trình bày toàn diện về bản chất vật lý, các cơ chế phát sinh dòng chảy, hệ thống phương pháp tính toán định lượng kinh điển và hiện đại, cùng những biến động sâu sắc của dòng chảy lũ dưới áp lực kép từ biến đổi khí hậu toàn cầu và hoạt động nhân sinh.

Bản chất vật lý và các đặc trưng cơ bản

Về mặt vật lý thủy văn, dòng chảy lũ đại diện cho pha động lực mãnh liệt nhất trong chu trình thủy văn của một lưu vực. Khi lượng mưa rơi xuống bề mặt vượt quá khả năng giữ nước của thảm phủ thực vật và dung lượng thấm của đất, lượng nước dư thừa sẽ liên kết lại thành dòng chảy mặt trực tiếp. Quá trình này diễn ra theo quy luật cơ học chất lưu với sự chuyển đổi thế năng địa hình thành động năng của dòng nước chảy xiết.

Để lượng hóa quy mô và cường độ của dòng chảy lũ, ngành thủy văn công trình sử dụng bốn nhóm thông số đặc trưng cốt lõi:

Lưu lượng đỉnh lũ: Ký hiệu là $Q_{\text{max}}$ hoặc $Q_P$ , là giá trị lưu lượng dòng chảy tức thời lớn nhất đo bằng mét khối trên giây xuất hiện tại một mặt cắt ngang khống chế trong suốt thời gian diễn ra trận lũ. Lưu lượng đỉnh lũ quyết định cao trình mực nước thiết kế cho các công trình đê điều, cống thoát nước và đập tràn hồ chứa.
Tổng lượng dòng chảy lũ: Ký hiệu là $W$ , biểu thị toàn bộ thể tích nước chảy qua mặt cắt lòng sông trong suốt trận lũ, tính bằng mét khối hoặc triệu mét khối. Về mặt toán học, tổng lượng dòng chảy lũ là tích phân của lưu lượng theo thời gian:
$W = \int_{t_0}^{t_k} Q(t) \, dt$
trong đó $W$ là tổng lượng dòng chảy lũ, $Q(t)$ là hàm lưu lượng biến thiên theo thời gian $t$ , $t_0$ là thời điểm bắt đầu xuất hiện lũ và $t_k$ là thời điểm trận lũ kết thúc hoàn toàn.
Đường quá trình lũ: Biểu đồ thể hiện sự biến thiên liên tục của lưu lượng nước theo thời gian. Đường quá trình lũ điển hình gồm hai pha rõ rệt: pha lũ lên tương ứng với thời gian lũ dâng do nước từ các sườn dốc tập trung dồn dập về lòng sông, và pha lũ rút tương ứng với thời gian lũ hạ khi nước trong hệ thống thoát dần ra ngoài kết hợp sự tháo cạn của các túi nước tạm thời.
Mô-đun đỉnh lũ: Lưu lượng đỉnh lũ tính trên một đơn vị diện tích lưu vực, phản ánh hiệu suất sinh lũ cục bộ và cho phép so sánh mức độ ác liệt của lũ giữa các lưu vực có quy mô không gian khác nhau.

Đặc trưng dòng chảy lũ	Ký hiệu quy ước	Đơn vị đo lường	Ý nghĩa trong thiết kế công trình
Lưu lượng đỉnh lũ lớn nhất	$Q_{\text{max}}$	Mét khối trên giây ( $\text{m}^3/\text{s}$ )	Xác định kích thước khẩu độ cống, cầu và chiều cao tràn xả lũ
Tổng lượng nước đợt lũ	$W$	Mét khối ( $\text{m}^3$ )	Tính toán dung tích phòng lũ của hồ chứa và khả năng cắt lũ hạ du
Thời gian duy trì lũ lên	$t_l$	Giờ ( $\text{h}$ )	Cung cấp thời gian dự báo cảnh báo sớm cho công tác ứng phó khẩn cấp
Thời gian duy trì lũ rút	$t_x$	Giờ ( $\text{h}$ )	Đánh giá thời gian ngập úng kéo dài của vùng đồng bằng trũng thấp

Các cơ chế hình thành dòng chảy lũ trên lưu vực

Sự hình thành dòng chảy lũ là kết quả tổng hòa của nhiều cơ chế vật lý tương tác phức tạp giữa khí quyển, thảm phủ thực vật, cấu trúc thổ nhưỡng và mạng lưới địa mạo. Khoa học thủy văn hiện đại phân định ba cơ chế phát sinh dòng chảy bão tố then chốt.

Dòng chảy tràn do mưa vượt cường độ thấm

Cơ chế dòng chảy tràn do mưa vượt cường độ thấm, thường gọi là dòng chảy tràn Hortonian, do nhà thủy văn học Robert E. Horton đề xuất trong công trình nghiên cứu công bố năm 1933. Cơ chế này xảy ra khi cường độ mưa tức thời vượt quá khả năng thấm nước của bề mặt đất tại một thời điểm xác định.

Horton (1933) đã làm sáng tỏ rằng khả năng thấm của đất không phải là một hằng số bất biến, mà suy giảm dần theo hàm mũ trong suốt đợt mưa từ tốc độ thấm ban đầu xuống tới tốc độ thấm ổn định. Nguyên nhân vật lý của hiện tượng này bao gồm sự lấp đầy các khoảng rỗng mao dẫn bởi các phân tử nước, sự trương nở của các hạt sét trong đất và tác động đập vỡ cơ học của hạt mưa làm nén chặt lớp bề mặt. Khi cường độ mưa lớn hơn tốc độ thấm, lượng nước thừa ban đầu lấp đầy các vùng trũng vi địa hình rồi liên kết lại tạo thành lớp màng nước chảy tràn liên tục trên sườn dốc hướng về mạng lưới lòng dẫn.

Dòng chảy tràn Hortonian là cơ chế chủ đạo tạo nên dòng chảy lũ tại các vùng đất trống đồi trọc, các sườn núi có lớp phủ đất mỏng, các vùng khí hậu bán khô hạn thiếu thảm thực vật hoặc các lưu vực đô thị hóa nơi bề mặt tự nhiên bị thay thế bởi bê tông và nhựa đường.

Dòng chảy tràn do bão hòa từ đáy lên

Trái ngược với quan điểm truyền thống của Horton, nghiên cứu thực nghiệm mang tính bước ngoặt của Thomas Dunne và Richard D. Black công bố năm 1970 đã chứng minh sự tồn tại của cơ chế dòng chảy tràn bão hòa, thường gọi là dòng chảy tràn Dunne. Cơ chế này xuất hiện khi toàn bộ tầng đất đã ngậm nước đến mức bão hòa hoàn toàn từ đáy lên đến mặt đất, khiến mực nước ngầm nông dâng cao trùng khớp với cao trình bề mặt.

Trong điều kiện đó, dù cường độ mưa rơi xuống rất nhỏ, đất cũng không thể tiếp nhận thêm bất kỳ lượng nước nào, và toàn bộ vũ lượng rơi trên vùng đất bão hòa lập tức biến thành dòng chảy tràn bề mặt. Dunne và Black (1970) phát hiện rằng dòng chảy lũ dạng này chỉ hình thành từ các diện tích đóng góp cục bộ chiếm từ 2% đến 10% diện tích lưu vực, chủ yếu tập trung tại các thung lũng đáy bằng, ven bờ sông suối hoặc chân các sườn dốc tích tụ nước. Cơ chế Dunne giữ vai trò quyết định tại các lưu vực nhiệt đới ẩm có lớp đất dày, thảm rừng rậm rạp và độ dẫn thủy lực cao.

Dòng chảy ngầm sát mặt và khái niệm diện tích đóng góp biến đổi

Bên cạnh hai dạng dòng chảy tràn bề mặt, dòng chảy ngầm sát mặt hình thành khi nước mưa thấm xuống tầng đất mặt nhưng bị chặn lại bởi một tầng đất sét hoặc nền đá ít thấm ở tầng dưới, buộc khối nước phải di chuyển ngang với vận tốc tương đối nhanh theo độ dốc địa hình để đổ vào lòng suối. Dòng chảy ngầm sát mặt có thể tái xuất hiện trên bề mặt tại các điểm chuyển tiếp địa hình, kết hợp với dòng chảy mặt làm gia tăng đáng kể lưu lượng đỉnh lũ.

Để mô phỏng động lực không gian của quá trình này, Keith J. Beven và Mike J. Kirkby công bố năm 1979 mô hình TOPMODEL dựa trên chỉ số địa hình: $\ln(a / \tan \beta)$ với $a$ là diện tích thu nước sườn dốc trên một đơn vị chiều rộng đường đồng mức và $\beta$ là góc dốc địa hình cục bộ. Beven và Kirkby (1979) đã phát triển khái niệm diện tích đóng góp biến đổi, chứng minh rằng diện tích sinh dòng chảy lũ không phải là một vùng cố định mà có tính động học rất cao: vùng sinh dòng chảy mở rộng nhanh chóng khi mưa kéo dài làm dâng mực nước ngầm, và co hẹp lại khi mưa ngớt.

Tiêu chí so sánh	Dòng chảy tràn Hortonian	Dòng chảy tràn Dunne
Nguyên nhân khởi phát	Cường độ mưa tức thời vượt quá tốc độ thấm của đất	Tầng đất ngậm no nước, mực nước ngầm dâng sát mặt
Thời điểm xuất hiện chính	Ngay trong các đợt mưa rào cường độ mạnh ban đầu	Sau thời gian mưa kéo dài khiến lưu vực tích ẩm tối đa
Phạm vi không gian điển hình	Phổ biến trên sườn dốc trơ trọi, vùng đất khô, khu đô thị	Khu trú ở dải đất trũng ven sông, đáy thung lũng ẩm ướt
Ảnh hưởng của thảm phủ	Rất nhạy cảm với lớp phủ thực vật và độ che phủ rừng	Chịu chi phối mạnh bởi dung tích trữ nước ngầm và địa mạo

Các phương pháp tính toán và mô hình hóa dòng chảy lũ

Tính toán dự báo các đại lượng dòng chảy lũ là yêu cầu kỹ thuật bắt buộc để thiết kế hệ thống công trình thủy lợi, cầu đường và kiểm soát rủi ro thiên tai. Trải qua hơn một thế kỷ phát triển, ngành thủy văn đã xây dựng hệ thống phương pháp từ bán thực nghiệm đơn giản đến các mô hình toán học tích hợp cao.

Phương pháp lý luận

Phương pháp lý luận có lịch sử phát triển từ đề xuất của Thomas Mulvaney năm 1851 và được Emil Kuichling hoàn thiện định lượng vào năm 1889. Đây là phương pháp bán thực nghiệm phổ biến nhất trên thế giới để xác định lưu lượng đỉnh lũ cho các công trình thoát nước quy mô nhỏ.

Trong hệ thống đơn vị đo lường quốc tế, phương pháp lý luận được biểu diễn qua công thức toán học:

Q = 0.278 \cdot C \cdot I \cdot A

trong đó

Q

là lưu lượng đỉnh lũ thiết kế tính bằng mét khối trên giây,

C

là hệ số dòng chảy lũ không thứ nguyên phản ánh tỷ lệ lượng mưa chuyển thành dòng chảy mặt,

I

là cường độ mưa tính toán tính bằng milimét trên giờ tương ứng với thời gian tập trung dòng chảy của lưu vực,

A

là diện tích lưu vực tính bằng kilômét vuông, và

0.278

là hệ số chuyển đổi đơn vị đo lường.

Phương pháp lý luận dựa trên giả thiết cơ bản rằng mưa phân bố hoàn toàn đồng đều trên toàn bộ lưu vực và lưu lượng đỉnh xuất hiện khi toàn bộ diện tích lưu vực đều tham gia đóng góp dòng chảy. Công trình nghiên cứu thực nghiệm của Young, McEnroe và Rome (2009) khi khảo sát các lưu vực có thời gian tập trung dao động từ 0,5 giờ đến 24 giờ đã chứng minh rằng hệ số dòng chảy $C$ phụ thuộc phi tuyến vào độ dốc, loại đất và chu kỳ lặp của trận mưa. Do các giả thiết lý tưởng hóa, Young, McEnroe và Rome (2009) nhấn mạnh rằng phương pháp lý luận chỉ được áp dụng an toàn cho các lưu vực nhỏ có diện tích dưới 25 km2 (tương đương 2500 ha).

Phương pháp Đường số cong dòng chảy

Phương pháp Đường số cong dòng chảy do Cơ quan Bảo tồn Đất đai Hoa Kỳ phát triển và được S. K. Mishra cùng V. P. Singh hệ thống hóa toàn diện vào năm 2003. Phương pháp này dựa trên nguyên lý cân bằng nước để tính toán lớp dòng chảy lũ tích lũy từ tổng lượng mưa của một trận bão.

Công thức cơ bản của phương pháp được viết dưới dạng:

P_e = \frac{(P - I_a)^2}{P - I_a + S} \quad \text{khi } P > I_a

trong đó

P_e

là lớp dòng chảy lũ tích lũy hiệu dụng tính bằng milimét,

P

là tổng lượng mưa của trận bão tính bằng milimét,

I_a

là lượng tổn thất ban đầu bao gồm tổn thất do thấm ban đầu, đọng trũng và chặn giữ bởi tán lá thực vật tính bằng milimét, và

S

là khả năng giữ nước tiềm năng tối đa của lưu vực tính bằng milimét.

Mishra và Singh (2003) phân tích rằng trong mô hình truyền thống, lượng tổn thất ban đầu thường được giả định bằng 0,2 $S$ (tương đương 20% khả năng giữ nước tối đa). Tuy nhiên, các nghiên cứu hiệu chỉnh hiện đại của Mishra và Singh (2003) chỉ ra rằng việc hạ tỷ lệ tổn thất ban đầu xuống 0,05 $S$ giúp mô hình phản ánh sát thực tế hơn đối với nhiều dạng địa hình sườn dốc. Đại lượng $S$ liên hệ trực tiếp với chỉ số Đường số cong $CN$ theo công thức hệ đơn vị mét:

S = \frac{25400}{CN} - 254

Chỉ số Đường số cong

CN

là một giá trị không thứ nguyên chạy trên thang đo từ 0 đến 100, được xác định dựa trên tổ hợp của bốn nhóm đất thủy văn, hiện trạng sử dụng đất, biện pháp canh tác và cấp độ ẩm tiền đề của lưu vực trước khi mưa xảy ra.

Phương pháp Đường quá trình đơn vị

Được Leroy K. Sherman đề xuất vào năm 1932, phương pháp đường quá trình đơn vị là một trong những bước tiến quan trọng nhất của thủy văn học hiện đại. Sherman (1932) định nghĩa đường quá trình đơn vị là đường quá trình dòng chảy mặt trực tiếp sinh ra bởi 1 cm hoặc 1 mm lớp nước mưa hiệu dụng phân bố đều trên khắp lưu vực trong một khoảng thời gian mưa đơn vị xác định.

Phương pháp này vận hành dựa trên hai nguyên lý cốt lõi của lý thuyết hệ thống tuyến tính dừng: nguyên lý tỷ lệ (lưu lượng tại mọi thời điểm tỷ lệ thuận với tổng lượng mưa hiệu dụng) và nguyên lý cộng dồn (đường quá trình tổng hợp của nhiều đợt mưa liên tiếp bằng tổng đại số các đường quá trình thành phần dịch chuyển theo trục thời gian). Mặc dù các trận lũ cực đoan thường bộc lộ tính chất phi tuyến do vận tốc truyền sóng lũ phụ thuộc vào độ sâu ngập, đường quá trình đơn vị vẫn là công cụ nền tảng trong các mô hình thủy văn lưu vực hiện hành.

Tính toán dòng chảy lũ theo tiêu chuẩn Việt Nam

Tại Việt Nam, công tác tính toán dòng chảy lũ thiết kế chịu sự chi phối chặt chẽ của các tiêu chuẩn quốc gia. Tiêu chuẩn Quốc gia TCVN 9845:2013 do Bộ Khoa học và Công nghệ công bố năm 2013 quy định phương pháp tính toán các đặc trưng dòng chảy lũ do mưa rào phục vụ thiết kế công trình thoát nước nhỏ trên các lưu vực sông có diện tích dưới 100 km2 không chịu ảnh hưởng của thủy triều và lũ bùn đá.

Theo quy định của TCVN 9845:2013, lưu lượng đỉnh lũ thiết kế ứng với tần suất tính toán được xác định theo công thức cường độ giới hạn:

Q_{P\%} = A_{P\%} \cdot \varphi \cdot H_{P\%} \cdot F \cdot \delta

trong đó các đại lượng được định nghĩa cụ thể:

$Q_{P\%}$ là lưu lượng đỉnh lũ thiết kế ứng với tần suất thiết kế $P\%$ , tính bằng mét khối trên giây;
$A_{P\%}$ là mô-đun tương đối đỉnh lũ ứng với tần suất thiết kế $P\%$ , phụ thuộc vào vùng mưa và thông số đặc trưng hình thái lưu vực;
$\varphi$ là hệ số dòng chảy lũ không thứ nguyên, được tra theo cấp đất cấu tạo lưu vực, diện tích hứng nước và lượng mưa ngày thiết kế;
$H_{P\%}$ là lượng mưa ngày lớn nhất ứng với tần suất thiết kế $P\%$ , tính bằng milimét;
$F$ là diện tích lưu vực tính toán, tính bằng kilômét vuông;
$\delta$ là hệ số chiết giảm đỉnh lũ do tác động điều tiết của ao hồ hoặc đầm lầy trên lưu vực; trường hợp lưu vực tự nhiên không có ao hồ đầm lầy điều tiết thì lấy hệ số $\delta = 1$ .

Tác động của biến đổi khí hậu và hoạt động nhân sinh

Chế độ dòng chảy lũ trên toàn cầu đang trải qua những biến động sâu sắc dưới tác động cộng hưởng của xu thế ấm lên toàn cầu và áp lực nhân sinh mãnh liệt lên bề mặt lưu vực.

Biến động dòng chảy lũ quy mô lục địa

Công trình nghiên cứu đột phá của Günter Blöschl và cộng sự công bố năm 2019 trên tạp chí Nature đã phân tích chuỗi dữ liệu quan trắc thực đo từ 3738 trạm thủy văn trải rộng khắp châu Âu trong giai đoạn kéo dài từ năm 1960 đến năm 2010. Kết quả nghiên cứu đã chứng minh rằng biến đổi khí hậu không làm tăng dòng chảy lũ một cách đồng nhất, mà tạo ra sự phân hóa không gian rõ rệt tùy thuộc vào cơ chế khí hậu chi phối.

Blöschl và cộng sự (2019) ghi nhận tại khu vực Trung và Tây Bắc Âu, lưu lượng dòng chảy lũ tăng tới 11,4% mỗi thập kỷ do lượng mưa mùa thu và mùa đông tăng cao khi độ ẩm khí quyển gia tăng. Ngược lại, tại khu vực Nam Âu, lưu lượng dòng chảy lũ lại sụt giảm tới 23,1% mỗi thập kỷ do tổng lượng mưa suy giảm và tốc độ bốc thoát hơi nước bề mặt gia tăng mạnh mẽ. Tại Đông Âu, dòng chảy lũ mùa xuân do tuyết tan giảm rõ rệt vì hiện tượng băng tuyết tan sớm hơn với lớp phủ tuyết mỏng hơn. Phát hiện này khẳng định tính bất định cao trong dự báo dòng chảy lũ tương lai.

Áp lực đô thị hóa và thay đổi sử dụng đất

Hoạt động nhân sinh, đặc biệt là quá trình đô thị hóa và chuyển đổi mục đích sử dụng đất, làm thay đổi căn bản cơ chế sinh dòng chảy lũ. Việc san ủi bề mặt và phủ kín các bề mặt thấm bằng bê tông, đường nhựa làm triệt tiêu khả năng thấm Hortonian tự nhiên của đất, giảm thiểu lượng trữ nước tại chỗ.

Hệ quả trực tiếp là hệ số dòng chảy lũ tăng vọt, thời gian tập trung dòng chảy bị rút ngắn nghiêm trọng do hệ thống thoát nước bê tông hóa vận chuyển nước nhanh hơn nhiều lần so với các lạch nước tự nhiên. Đỉnh lũ đô thị do đó xuất hiện sớm hơn, có biên độ dâng cao hơn và mức độ phá hủy khốc liệt hơn so với trạng thái tự nhiên ban đầu.

Nghiên cứu mô phỏng dòng chảy lũ tại Việt Nam

Tại Việt Nam, sự kết hợp giữa địa hình đồi núi dốc, chiều dài bờ biển hẹp và lượng mưa bão nhiệt đới tập trung cao khiến dòng chảy lũ luôn tiềm ẩn rủi ro thiên tai đặc biệt lớn. Nghiên cứu của Nguyen và cộng sự công bố năm 2023 trên tạp chí Vietnam Journal of Earth Sciences đã xây dựng khung mô hình thủy động lực học kết hợp học máy tiên tiến để dự báo độ sâu ngập lũ cho các tỉnh dải ven biển miền Trung Việt Nam.

Nguyen và cộng sự (2023) nhấn mạnh rằng các lưu vực dốc ven biển có đặc điểm địa hình hẹp ngang với thời gian tập trung dòng chảy lũ rất ngắn dưới 12 giờ. Khi các trận mưa lớn xảy ra đồng thời với hiện tượng triều cường dâng cao ngăn dòng chảy thoát ra biển, dòng chảy lũ bị nghẽn tắc cục bộ, gây ra ngập lụt diện rộng với vận tốc dòng chảy xiết và sức tàn phá dữ dội đối với cơ sở hạ tầng ven bờ.

Phạm vi áp dụng, hạn chế và các vấn đề mở

Mặc dù các lý thuyết và mô hình tính toán dòng chảy lũ đã đạt được nhiều thành tựu lớn, việc áp dụng vào thực tiễn kỹ thuật vẫn đối mặt với những rào cản và giới hạn khoa học quan trọng.

Trước hết là giới hạn về quy mô không gian và tính phi tuyến của lưu vực. Các công thức kinh nghiệm như phương pháp lý luận hay phương pháp cường độ giới hạn của TCVN 9845:2013 chỉ cho kết quả tin cậy trên các lưu vực quy mô nhỏ dưới 100 km2. Khi quy mô lưu vực mở rộng lên hàng nghìn kilômét vuông, sự không đồng đều của trường mưa theo không gian và thời gian trở nên lấn át, đồng thời hiệu ứng trữ nước và làm phẳng sóng lũ của lòng dẫn chính đòi hỏi phải tích hợp các phương trình vi phân thủy lực một chiều và hai chiều phức tạp thay vì các công thức giải tích đơn giản.

Thứ hai, cần phân biệt rạch ròi giữa dòng chảy lũ và hiện tượng ngập lụt. Dòng chảy lũ là một đại lượng động lực học thủy văn mô tả lượng nước và vận tốc nước chuyển động trên lưu vực; trong khi ngập lụt là hiện tượng ngập nước vùng đất khô xảy ra khi dòng chảy lũ vượt quá khả năng chuyển tải thủy lực của lòng dẫn tự nhiên hoặc nhân tạo. Một trận lũ lớn không nhất thiết gây ra thảm họa ngập lụt nếu lòng sông sâu, đê điều kiên cố và hành lang thoát lũ được bảo vệ nguyên vẹn.

Thách thức khoa học lớn nhất hiện nay là bài toán tính toán dòng chảy lũ tại các lưu vực không có số liệu quan trắc thủy văn đầy đủ. Trước thách thức này, các phương pháp chuyển dịch tham số vùng, ứng dụng ảnh viễn thám độ phân giải cao và các thuật toán học máy đang được tích cực phát triển để giải mã động lực dòng chảy lũ. Đồng thời, xu hướng quản trị rủi ro lũ lụt hiện đại đang chuyển dịch từ các công trình ngăn lũ cứng bằng bê tông sang các giải pháp dựa vào thiên nhiên, như phục hồi dải ngập lũ tự nhiên, bảo vệ rừng phòng hộ đầu nguồn và xây dựng mô hình đô thị bọt biển nhằm khôi phục khả năng giữ chậm dòng chảy lũ của lưu vực.