Hình thái dòng sông khác gì so với động lực học dòng sông?

Hình thái dòng sông mô tả các đặc trưng hình học, kích thước và dạng mặt cắt của lòng dẫn tại các thời điểm quan sát. Trong khi đó, động lực học dòng sông tập trung nghiên cứu các quy luật vật lý, trường vận tốc và cơ chế thủy lực trực tiếp vận chuyển bùn cát và điều khiển quá trình xói lở hay bồi tụ.

Tại sao các dòng sông tự nhiên hiếm khi hoàn toàn thẳng?

Dòng chảy trong lòng dẫn tự nhiên luôn xuất hiện các cấu trúc dòng chảy rối và dòng chảy thứ sinh xoắn ốc gây ra sự phân bố không đối xứng của ứng suất cắt đáy. Sự mất đối xứng này làm khởi phát các bãi nông và hố sâu xen kẽ, khiến đường tim đáy sâu nhất liên tục uốn lượn và nhanh chóng biến lòng dẫn thẳng thành dạng uốn khúc.

Nguyên lý cân bằng Lane giải thích điều gì trong tiến hóa lòng sông?

Nguyên lý cân bằng Lane thể hiện mối tương quan tỷ lệ giữa tích của lưu lượng bùn cát với kích thước hạt trầm tích và tích của lưu lượng nước với độ dốc lòng dẫn. Khi một trong bốn yếu tố này bị thay đổi do đập thủy điện hoặc khai thác cát, lòng sông sẽ tự điều chỉnh độ dốc và kích thước hạt thông qua xói sâu hoặc bồi lắng để tái lập cân bằng mới.

Các hồ chứa lớn ảnh hưởng thế nào đến hình thái hạ lưu sông ngòi?

Các đập thủy điện giữ lại phần lớn bùn cát lơ lửng, tạo ra dòng nước xả thiếu hụt trầm tích có năng lượng bào mòn rất cao ở hạ du. Hiện tượng nước đói phù sa này gây xói sâu hạ thấp cao trình đáy sông, làm sụt giảm mực nước mùa kiệt, gia tăng nguy cơ sạt lở bờ và xâm nhập mặn tại các vùng đồng bằng châu thổ.

Hình thái dòng sông (river morphology) là gì?

Hình thái dòng sông là trạng thái hình học, kích thước và cấu trúc không gian của lòng dẫn nước tự nhiên, được định hình và biến đổi liên tục thông qua sự tương tác động lực giữa dòng chảy mặt, tải lượng bùn cát và điều kiện địa mạo lưu vực. Thuật ngữ này bao hàm toàn bộ các đặc trưng trắc diện dọc, trắc diện ngang và bình đồ mặt bằng của mạng lưới sông suối từ thượng nguồn vùng núi cao xuống đồng bằng châu thổ và cửa sông ven biển. Bài viết trình bày cơ sở lý thuyết, các quan hệ định lượng hình học thủy lực, hệ thống phân loại hình thái kinh điển và hiện đại, cùng những diễn biến thực tiễn dưới tác động nhân sinh tại các lưu vực sông lớn ở Việt Nam.

Khái niệm và bản chất của hình thái dòng sông

Trong địa mạo học sông ngòi và kỹ thuật thủy lực, hình thái dòng sông phản ánh sự cân bằng động giữa lực kéo của dòng nước chảy và sức kháng cản của vật liệu lòng bờ. Dòng sông không phải là một máng dẫn cứng nhắc cố định, mà là một hệ thống tự nhiên có khả năng tự tổ chức và tự thích ứng cao. Khi lưu lượng nước hoặc tải lượng trầm tích đưa vào hệ thống thay đổi, lòng dẫn sẽ tự điều chỉnh độ rộng, độ sâu, độ dốc hoặc độ uốn khúc để thiết lập trạng thái cân bằng mới phù hợp với điều kiện biên năng lượng.

Cần phân biệt rõ ràng giữa hình thái dòng sông và động lực học dòng sông nhằm tránh nhầm lẫn trong nghiên cứu và ứng dụng thực tiễn:

Hình thái dòng sông: Mô tả biểu hiện hình học bên ngoài và cấu trúc không gian của lòng dẫn tại một thời điểm hoặc diễn tiến biến đổi hình học qua các thời kỳ quan trắc, bao gồm bình đồ mặt bằng, hình dạng mặt cắt ngang và độ dốc đáy trắc diện dọc.
Động lực học dòng sông: Nghiên cứu các lực tác động, trường vận tốc thủy lực ba chiều, ứng suất cắt đáy và cơ chế vật lý trực tiếp làm dịch chuyển hạt trầm tích, gây xói mòn hoặc bồi tụ dẫn đến sự biến đổi hình thái.
Hình thái học sông ngòi: Phân ngành khoa học liên ngành nghiên cứu tổng hợp mối quan hệ biện chứng giữa quá trình động lực học thủy lực - trầm tích và kết quả biểu hiện hình thái của dòng sông theo không gian và thời gian.

Khái niệm cân bằng động giữ vai trò trung tâm trong phân tích hình thái dòng sông. Một đoạn sông ở trạng thái cân bằng động không có nghĩa là lòng sông hoàn toàn bất động. Ngược lại, các quá trình bồi tụ và xói lở vẫn diễn ra liên tục tại từng vị trí cục bộ, chẳng hạn như xói lở ở bờ lõm đi kèm bồi tụ bãi bồi ở bờ lồi của khúc uốn, nhưng kích thước mặt cắt ngang trung bình và độ dốc chung của toàn đoạn sông vẫn duy trì tương đối ổn định qua nhiều chu kỳ mùa lũ và mùa kiệt.

Các quy luật định lượng và quan hệ cân bằng hình thái

Sự hình thành và biến đổi của hình thái lòng sông được chi phối bởi các quy luật thủy lực và trầm tích có thể mô tả qua các công thức định lượng chuẩn xác.

Quan hệ hình học thủy lực Leopold và Maddock

Nghiên cứu kinh điển của Leopold và Maddock năm 1953 đã đặt nền móng định lượng cho việc nghiên cứu hình thái dòng sông thông qua khái niệm hình học thủy lực. Nhóm tác giả chứng minh rằng các biến số hình thái chính gồm bề rộng mặt nước, độ sâu dòng chảy trung bình và vận tốc dòng chảy trung bình đều biến thiên theo hàm lũy thừa đơn giản của lưu lượng nước:

w = a Q^b

d = c Q^f

v = k Q^m

Trong các phương trình trên, ký hiệu $w$ là bề rộng mặt nước của lòng dẫn, $d$ là độ sâu dòng chảy trung bình, $v$ là vận tốc dòng chảy trung bình của mặt cắt, và $Q$ là lưu lượng dòng chảy qua mặt cắt ngang. Các đại lượng $a$ , $c$ , và $k$ là các hệ số tỷ lệ thực nghiệm; trong khi $b$ , $f$ , và $m$ là các số mũ hình thái phản ánh tốc độ biến thiên tương đối của từng đại lượng.

Dựa trên phương trình bảo toàn lưu lượng liên tục $Q = w \cdot d \cdot v$ , tích các hệ số và tổng các số mũ bắt buộc phải thỏa mãn ràng buộc toán học:

a \cdot c \cdot k = 1

b + f + m = 1

Leopold và Maddock năm 1953 đã tổng hợp dữ liệu quan trắc thủy văn thực tế và xác định giá trị trung bình của các số mũ dọc theo hướng phát triển lưu vực từ thượng lưu về hạ lưu. Kết quả cho thấy số mũ bề rộng $b \approx 0.5$ , số mũ độ sâu $f \approx 0.4$ , và số mũ vận tốc $m \approx 0.1$ . Điều này phản ánh quy luật địa mạo rằng khi tiến dần về phía hạ lưu với lưu lượng dòng chảy tăng lên, lòng sông có xu hướng mở rộng bề ngang nhanh nhất, tiếp đến là sự gia tăng độ sâu, trong khi vận tốc dòng chảy trung bình chỉ tăng nhẹ hoặc thậm chí giữ giá trị xấp xỉ không đổi.

Quan hệ cân bằng định tính Lane

Để giải thích cơ chế phản ứng hình thái trước sự mất cân bằng giữa tải lượng bùn cát và dòng chảy, Lane năm 1955 đã đề xuất mối quan hệ cân bằng nổi tiếng:

Q_s \cdot D_{50} \propto Q \cdot S

Trong đó, đại lượng $Q_s$ là lưu lượng vận chuyển bùn cát của dòng sông, $D_{50}$ là đường kính hạt trầm tích trung vị của vật liệu đáy, $Q$ là lưu lượng nước của dòng chảy, và $S$ là độ dốc dọc của đáy lòng dẫn. Mối quan hệ Lane đóng vai trò như một chiếc cân đòn bẩy mô tả cân bằng trầm tích:

Khi tích số giữa dòng chảy và độ dốc lớn hơn khả năng vận chuyển bùn cát, năng lượng dư thừa của dòng nước sẽ bào xói vật liệu đáy và hai bên bờ, dẫn đến hiện tượng xói sâu hạ thấp đáy lòng dẫn.
Khi nguồn cung cấp bùn cát vượt quá khả năng vận chuyển của dòng nước, lượng trầm tích dư thừa sẽ bồi lắng lại, làm đáy sông được nâng cao dần và lòng dẫn có xu hướng mở rộng sang hai bên.

Hệ số uốn khúc và ngưỡng chuyển tiếp hình thái

Hệ số uốn khúc là một chỉ số hình học cơ bản dùng để định lượng mức độ ngoằn ngoèo của dòng sông trên bình đồ mặt bằng. Đại lượng này được xác định bằng tỷ số giữa chiều dài dòng dẫn thực tế và chiều dài thung lũng:

P = \frac{L_c}{L_v}

Trong công thức trên, ký hiệu $P$ là hệ số uốn khúc, $L_c$ là chiều dài lòng dẫn đo dọc theo tim dòng chảy sâu nhất, và $L_v$ là khoảng cách đường thẳng đo dọc theo trục thung lũng sông giữa hai điểm nghiên cứu. Theo tiêu chuẩn phân loại hình thái học kinh điển của Leopold và Wolman năm 1957:

Đoạn sông được xem là lòng dẫn thẳng khi hệ số uốn khúc $P < 1.05$ .
Đoạn sông có dạng uốn khúc chuyển tiếp hoặc uốn khúc nhẹ khi hệ số uốn khúc nằm trong khoảng từ 1.05 đến dưới 1.5.
Đoạn sông được phân loại là lòng dẫn uốn khúc điển hình khi hệ số uốn khúc đạt giá trị $P \ge 1.5$ .

Leopold và Wolman năm 1957 đã thiết lập phương trình ngưỡng thực nghiệm phân tách giữa kiểu lòng dẫn bện tóc và lòng dẫn uốn khúc dựa trên độ dốc và lưu lượng tạo lòng:

S = 0.06 Q^{-0.44}

Trong phương trình ngưỡng này của Leopold và Wolman năm 1957, đại lượng $S$ là độ dốc dọc lòng dẫn và $Q$ là lưu lượng tạo lòng tính theo đơn vị feet khối trên giây (cfs). Trong hệ đơn vị quốc tế (SI) với lưu lượng tính bằng mét khối trên giây, công thức ngưỡng tương đương có dạng xấp xỉ $S \approx 0.012 Q^{-0.44}$ . Với một lưu lượng nước nhất định, nếu độ dốc lòng dẫn vượt quá ngưỡng giới hạn, dòng chảy sẽ chuyển sang dạng bện tóc với nhiều nhánh rẽ đan xen; ngược lại, nếu độ dốc nhỏ hơn ngưỡng này, lòng sông sẽ ổn định dưới dạng lòng dẫn uốn khúc đơn.

Phân loại hình thái dòng sông

Trong lịch sử phát triển của hình thái học sông ngòi, nhiều hệ thống phân loại đã được xây dựng nhằm hệ thống hóa tính đa dạng của các lòng dẫn tự nhiên phục vụ quy hoạch, kỹ thuật chỉnh trị và bảo vệ môi trường sinh thái lưu vực.

Các kiểu hình thái lòng dẫn kinh điển

Dựa trên bình đồ mặt bằng và cơ chế bồi xói, các nhà địa mạo học truyền thống phân biệt bốn kiểu lòng dẫn cơ bản:

Kiểu hình thái	Đặc trưng bình đồ	Tải lượng trầm tích chủ đạo	Tính ổn định tương đối
Lòng dẫn thẳng	Lòng đơn ít uốn lượn, hệ số uốn khúc dưới 1.05	Hỗn hợp trầm tích lơ lửng và trầm tích đáy	Thấp đến trung bình do lạch sâu vẫn uốn lượn liên tục
Lòng dẫn uốn khúc	Lòng đơn phát triển các khúc uốn liên tục, hệ số uốn khúc từ 1.5 trở lên	Chủ yếu là hạt mịn và bùn cát lơ lửng	Trung bình, bờ sông dịch chuyển dần theo thời gian qua bồi lồi xói lõm
Lòng dẫn bện tóc	Mạng lưới nhiều lạch nông phân nhánh đan xen qua các cù lao cát	Tải lượng vật liệu đáy thô rất lớn và độ dốc cao	Rất kém ổn định, các lạch dòng chảy dịch chuyển liên tục sau mỗi mùa lũ
Lòng dẫn phân nhánh	Nhiều lòng dẫn sâu ổn định chia cắt bởi các đảo phù sa lớn có thực vật	Trầm tích hạt mịn kết hợp đất dính bền vững	Cao, hệ thống lòng dẫn duy trì vị trí ổn định qua nhiều thập kỷ

Hệ thống phân loại hình thái Rosgen

Hệ thống phân loại do Rosgen công bố năm 1994 là một trong những khung đánh giá định lượng được ứng dụng rộng rãi nhất trên thế giới. Phương pháp này phân loại lòng dẫn tự nhiên theo các bậc thứ bậc dựa trên bốn chỉ số hình học mặt cắt đo đạc thực tế:

Tỷ số giới hạn lòng dẫn: Tỷ số giữa bề rộng vùng ngập lũ đo ở cao trình bằng hai lần độ sâu tối đa thời điểm đầy bờ với bề rộng mặt nước đầy bờ. Tỷ số nhỏ hơn 1.4 đặc trưng cho lòng sông ăn sâu dạng hẻm vực hẹp, trong khi giá trị lớn hơn 2.2 biểu thị thung lũng thoải tiếp cận bãi bồi rộng.
Tỷ số bề rộng trên độ sâu: Tỷ số giữa bề rộng mặt nước đầy bờ và độ sâu trung bình đầy bờ. Các giá trị ngưỡng 12 và 40 được sử dụng để phân biệt lòng dẫn sâu hẹp với lòng dẫn nông rộng.
Hệ số uốn khúc: Đo mức độ uốn lượn mặt bằng với các ngưỡng phân tách 1.2 và 1.5.
Độ dốc lòng dẫn: Đo độ dốc đáy thực tế với các khoảng giá trị ranh giới 0.02 tức 2% và 0.04 tức 4%.

Dựa trên sự kết hợp của các chỉ số trên, Rosgen năm 1994 chia các dòng sông tự nhiên thành bảy kiểu cơ bản từ loại A đến loại G:

Kiểu lòng dẫn	Độ dốc đáy	Tỷ số giới hạn	Tỷ số rộng trên sâu	Hệ số uốn khúc	Mô tả hình thái
Loại A	Lớn hơn 0.04 (tức 4%)	Nhỏ hơn 1.4	Nhỏ hơn 12	Nhỏ hơn 1.2	Suối vùng núi cao dốc đứng, ăn sâu vào hẻm đá, dòng chảy xiết
Loại B	Từ 0.02 đến 0.04	Từ 1.4 đến 2.2	Lớn hơn 12	Lớn hơn 1.2	Đoạn sông trung du dốc vừa, lòng dẫn ổn định có bậc ghềnh và vũng sâu
Loại C	Nhỏ hơn 0.02 (dưới 2%)	Lớn hơn 2.2	Lớn hơn 12	Lớn hơn 1.2	Sông đồng bằng uốn khúc rộng rãi, bãi bồi thoải phát triển mạnh
Loại D	Nhỏ hơn 0.04	Không áp dụng	Lớn hơn 40	Rất nhỏ	Lòng dẫn bện tóc nhiều nhánh nông, biến động mạnh trên nền trầm tích thô
Loại E	Nhỏ hơn 0.02	Lớn hơn 2.2	Nhỏ hơn 12	Lớn hơn 1.5	Lòng dẫn uốn khúc rất mạnh, mặt cắt sâu hẹp, bờ sông ổn định nhờ thảm cỏ
Loại F	Nhỏ hơn 0.02	Nhỏ hơn 1.4	Lớn hơn 12	Lớn hơn 1.2	Lòng dẫn ăn sâu vào trầm tích phù sa, mất kết nối với bãi bồi lũ
Loại G	Từ 0.02 đến 0.04	Nhỏ hơn 1.4	Nhỏ hơn 12	Lớn hơn 1.2	Khe xói sâu dạng rãnh hẹp, thành đứng, độ bất ổn định và xói lở rất cao

Động lực biến đổi hình thái và tác động nhân sinh

Hình thái lòng dẫn không phải là một trạng thái tĩnh bất biến. Theo tổng quan khoa học của Church năm 2006, hình thái lòng sông phù sa là hệ quả vật lý trực tiếp của các quá trình vận chuyển bùn cát và bồi lắng liên tục. Sự cân bằng hình thái tự nhiên phụ thuộc chặt chẽ vào chu trình bổ cập trầm tích từ thượng nguồn. Giai đoạn từ năm 1950 đến nay, các hoạt động kinh tế - xã hội quy mô lớn của con người đã trở thành động lực chi phối mạnh mẽ nhất làm thay đổi căn bản hình thái các hệ thống sông ngòi trên toàn cầu cũng như tại Việt Nam.

Thiếu hụt bùn cát do hồ chứa và xói sâu hạ lưu sông Hồng

Hệ thống sông Hồng ở miền Bắc Việt Nam là minh chứng điển hình cho sự biến đổi hình thái dưới tác động của chuỗi hồ chứa thủy điện thượng lưu. Nghiên cứu thực nghiệm của Lu và cộng sự công bố năm 2015 đã chỉ ra rằng việc đưa vào vận hành các đập thủy điện lớn ở thượng nguồn sông Đà và sông Lô đã làm thay đổi hoàn toàn ngân sách bùn cát của vùng hạ du. Các hồ chứa đóng vai trò như những bẫy lắng nhân tạo khổng lồ, giữ lại khoảng 70% đến 80% tổng tải lượng bùn cát lơ lửng của dòng sông.

Dòng nước xả từ các đập thủy điện sau khi đã lắng sạch phù sa trở thành dòng nước thiếu trầm tích, có khả năng bào mòn mạnh mẽ khi di chuyển về hạ lưu. Hiện tượng này kết hợp với hoạt động khai thác cát xây dựng rầm rộ tại vùng lòng sông đã làm mất cân bằng trầm tích nghiêm trọng theo nguyên lý Lane năm 1955. Nghiên cứu của Lu và cộng sự năm 2015 đã chỉ ra rằng lòng dẫn đoạn hạ lưu sông Hồng chảy qua khu vực Hà Nội đã bị xói sâu đáy rõ rệt sau khi chuỗi đập thủy điện đi vào vận hành. Hiện tượng hạ thấp đáy lòng dẫn làm sụt giảm nghiêm trọng mực nước mùa kiệt, gây khó khăn cho việc lấy nước của các cống tưới tiêu nông nghiệp truyền thống và làm gia tăng nguy cơ sạt lở chân đê điều bảo vệ vùng đồng bằng.

Khai thác cát và sạt lở bờ sông tại đồng bằng sông Cửu Long

Tại miền Nam Việt Nam, hạ lưu sông Mê Kông chảy qua đồng bằng sông Cửu Long đang đối mặt với cuộc khủng hoảng hình thái lòng dẫn do áp lực kép từ các hồ chứa thượng nguồn xuyên biên giới và khai thác cát ồ ạt trong nội địa. Công trình nghiên cứu đột phá của Hackney và cộng sự năm 2020 trên tạp chí Nature Sustainability đã công bố những số liệu định lượng gây chấn động về mức độ mất cân bằng trầm tích:

Tổng lưu lượng cát tự nhiên bổ cập cho vùng hạ lưu châu thổ sông Mê Kông ước tính chỉ đạt khoảng 6.18 triệu tấn mỗi năm, bao gồm cả lượng cát đẩy đáy và cát lơ lửng.
Trong khi đó, tổng khối lượng cát bị khai thác thương mại ước tính lên tới khoảng 50 triệu tấn mỗi năm, vượt xa nguồn cung cấp tự nhiên từ thượng lưu.

Sự chênh lệch khổng lồ giữa nguồn cung và khối lượng khai thác đã tạo ra hố sụt trầm tích nghiêm trọng. Phân tích mô hình địa kỹ thuật của Hackney và cộng sự năm 2020 chỉ ra rằng tốc độ xói sâu lòng dẫn cục bộ có thể đạt 0.5 m mỗi năm, và mức hạ thấp đáy sông từ 1.0 m đến 3.0 m là ngưỡng kích hoạt trạng thái mất ổn định và sạt lở bờ sông nghiêm trọng ở đồng bằng sông Cửu Long. Đáy sông hạ thấp làm tăng chiều cao tương đối của bờ dốc và gia tăng ứng suất cắt trọng lực, dẫn đến hàng trăm vụ sạt lở bờ sông nghiêm trọng mỗi năm, cuốn trôi nhà cửa, đường giao thông và đất đai canh tác của người dân.

Ứng dụng trong kỹ thuật chỉnh trị và tiêu chuẩn quốc gia

Hiểu biết sâu sắc về diễn biến hình thái lòng sông là tiền đề bắt buộc cho công tác thiết kế công trình thủy lợi, cảng đường thủy và chỉnh trị bờ sông. Tại Việt Nam, Tiêu chuẩn Quốc gia TCVN 8419:2022 do cơ quan có thẩm quyền ban hành đã quy định chi tiết các yêu cầu kỹ thuật trong thiết kế công trình bảo vệ bờ sông chống lũ. Tiêu chuẩn này đòi hỏi kỹ sư phải tiến hành khảo sát toàn diện các quy luật diễn biến lòng dẫn, đánh giá chiều sâu xói lở cục bộ và tổng thể, cũng như dự báo độ dịch chuyển của luồng lạch chủ đạo trước khi quyết định lựa chọn hình thức kè lát mái, kè mỏ hàn nắn dòng hoặc kè mềm bảo vệ chân bờ.

Hạn chế của các mô hình hình thái và vấn đề khoa học còn mở

Mặc dù các lý thuyết hình thái dòng sông đã đạt được nhiều bước tiến quan trọng, việc áp dụng các mô hình tĩnh vào thực tế quản lý lưu vực vẫn bộc lộ nhiều điểm nghẽn khoa học cần tiếp tục hoàn thiện.

Hệ thống phân loại Rosgen dù rất phổ biến nhưng đã phải đối mặt với sự phê bình học thuật sâu sắc. Nhiều nhà địa mạo học và kỹ sư chỉnh trị sông ngòi đã chỉ ra rằng việc phân loại lòng dẫn dựa thuần túy trên các chỉ số hình thái tĩnh không thể phản ánh được các quá trình vật lý động lực học phức tạp như vận chuyển bùn cát không ổn định, tương tác dòng chảy rối quanh bãi bồi và các ngưỡng chuyển tiếp trạng thái. Việc áp dụng máy móc phương pháp thiết kế lòng dẫn tự nhiên dựa trên phân loại Rosgen trong các dự án phục hồi sông ngòi đã dẫn đến nhiều thất bại thực tế do hệ thống này không tính toán được phản ứng thích ứng lâu dài của lòng dẫn khi các điều kiện thủy văn lưu vực biến động.

Bên cạnh đó, nghiên cứu của Church năm 2006 cũng nhấn mạnh tính phức tạp của thời gian phản hồi và thời gian thư giãn trong hệ thống sông phù sa. Khi một tác động ngoại cảnh xuất hiện, chẳng hạn như xây đập hoặc biến đổi lượng mưa cực đoan, lòng sông không phản ứng tức thời mà trải qua một quá trình điều chỉnh phi tuyến tính kéo dài hàng thập kỷ hoặc hàng thế kỷ. Mức độ trễ này khiến việc dự báo chính xác diễn biến hình thái trong bối cảnh biến đổi khí hậu toàn cầu và các can thiệp nhân sinh phức tạp tiếp tục là một trong những thách thức khoa học mở hàng đầu của ngành khoa học địa cầu hiện đại.