Tổng quan về giáo trình

Thủy lực sông ngòi – Tập I: Các nội dung lý thuyết do GS.TS. Hà Văn Khối, GS.TSKH. Nguyễn Ân Niên và PGS.TS. Đỗ Tất Túc biên soạn, được xuất bản năm 2007 trong khuôn khổ dự án DANIDA hỗ trợ tăng cường năng lực cho Trường Đại học Thủy lợi. Tài liệu được thiết kế phục vụ chương trình đào tạo trình độ cao học và nghiên cứu sinh thuộc các chuyên ngành: Thủy văn học, Chỉnh trị sông và bờ biển, và Phát triển nguồn nước.

Mục tiêu học tập của giáo trình là trang bị hệ thống lý thuyết thủy lực kênh hở nâng cao, phương pháp thiết lập bài toán thủy động lực học dòng chảy sông ngòi và các thuật toán số ứng dụng trong diễn toán thủy văn – thủy lực. Sau khi hoàn thành nội dung, người học có khả năng phân tích diễn biến hình thái sông bồi tích, thiết lập mô hình tính toán một chiều (1D), hai chiều (2D), xác định quy luật vận chuyển bùn cát, tính toán xói lở – bồi lắng và mô phỏng quá trình truyền chất, chất lượng nước cũng như xâm nhập mặn trong hệ thống sông.

Cấu trúc giáo trình gồm 4 phần chính với 16 chương lý thuyết:

  • Phần 1: Diễn toán dòng chảy trong sông (Chương 1 đến Chương 5, do GS.TS. Hà Văn Khối biên soạn)
  • Phần 2: Thủy động lực học sông (Chương 6 đến Chương 10, do GS.TSKH. Nguyễn Ân Niên biên soạn)
  • Phần 3: Vận chuyển bùn cát và chỉnh trị sông (Chương 11 đến Chương 14, do PGS.TS. Đỗ Tất Túc biên soạn)
  • Phần 4: Chất lượng nước (Chương 15 và Chương 16, do GS.TSKH. Nguyễn Ân Niên biên soạn)

Cách tiếp cận của giáo trình kết hợp chặt chẽ giữa cơ học chất lỏng lý thuyết, thực nghiệm thủy lực và kỹ thuật mô phỏng số hiện đại. Điểm đặc sắc là việc tích hợp các mô hình toán kinh điển của thế giới (như hệ phương trình Saint-Venant, mô hình sóng động học, sóng khuếch tán, các mô hình chảy rối NKE và RSM) với thực tiễn ứng dụng tại hệ thống sông Việt Nam (sông Hồng, sông Cửu Long) và các mô hình do các nhà khoa học trong nước phát triển như VRSAP, KOD01, KOD02.

Nội dung kiến thức cốt lõi

Các chương/chủ đề chính

Cấu trúc nội dung giáo trình phát triển theo trình tự logic chặt chẽ từ cơ sở động lực học dòng chảy hở, các mô hình diễn toán thủy văn - thủy lực, đến các quá trình tương tác đáy lòng dẫn (vận chuyển bùn cát, biến hình lòng sông) và lan truyền vật chất trong môi trường nước:

  • Chương 1: Mở đầu: Thiết lập hệ thống khái niệm về hệ thống sông 3 vùng theo phân loại của Schumm (1977) gồm vùng sinh thủy (vùng 1), vùng chuyển tiếp ổn định tương đối (vùng 2) và vùng cửa sông chịu ảnh hưởng triều (vùng 3). Phân loại các biến độc lập và phụ thuộc trong sông bồi tích theo quy mô thời gian (ngắn hạn, dần dần theo thời gian, thời gian địa chất). Giới thiệu khái niệm trạng thái bán cân bằng động lực (Mackin), công thức lưu lượng tạo lòng của William: $$Q = 4{,}28 \cdot A_f \cdot S^{0{,}5}$$ quy luật mặt cắt dọc sông theo Shulits (1941) $S = S_0 e^{-\alpha x}$, định luật giảm cấp phối hạt $d = d_0 e^{-\beta x}$, các tham số không thứ nguyên như số Froude ($F = \frac{V}{\sqrt{gh}}$), số Reynolds và phân tích các mô hình dòng chảy rối phi tuyến NKE (Pope), mô hình vận chuyển ứng suất Reynolds (RSM).
  • Chương 2: Diễn toán dòng chảy trong sông theo mô hình thông số tập trung: Thiết lập phương trình cân bằng nước $\frac{dS}{dt} = I(t) - Q(t)$. Phân tích đặc tính hồ chứa tĩnh (mặt nước nằm ngang, quan hệ $S = f(Q)$ đơn trị) và hồ chứa động (ảnh hưởng nước vật, quan hệ $S - Q$ dạng vòng dây). Trình bày thuật giải bán đồ giải V.T. Chow, phương pháp Runge-Kutta bậc 3, phương pháp lặp trực tiếp và phương pháp Muskingum với phương trình lượng trữ tuyến tính $S = K[xI + (1-x)Q]$.
  • Chương 3, 4, 5: Mô hình thông số phân bố, sóng động học và sóng khuếch tán: Xây dựng hệ phương trình Saint-Venant 1D đầy đủ gồm phương trình liên tục và phương trình động lượng có kể đến lưu lượng gia nhập bên ($q$), độ dốc ma sát ($J$), tổn thất cục bộ ($i_e$), lực ma sát gió ($W_f$). Phân loại các cấp độ mô hình giản lược: mô hình sóng động học giải bằng sai phân hữu hạn; mô hình sóng khuếch tán giải bằng biến đổi Laplace cho kênh tuyến tính không có/có nhập lưu bên và bài toán phi tuyến.
  • Chương 7, 8, 9, 10: Thiết lập bài toán thủy lực, dòng chảy ổn định, không ổn định và dòng trong sông cong: Phương pháp giải số hệ Saint-Venant cho dòng 1D, sơ đồ mạng lưới, điều kiện biên, số liệu đầu vào, hiệu chỉnh và kiểm định mô hình. Phân tích cấu trúc dòng chảy xoắn thứ cấp (secondary currents) trong đoạn sông uốn khúc và hệ phương trình cơ bản cho dòng chảy cong.
  • Chương 11, 12, 13, 14: Dòng chảy bùn cát, hình thái và chỉnh trị sông: Xác định trạng thái chuyển động ban đầu của bùn cát không cố kết, thiết kế mặt cắt kênh ổn định, tính toán độ nhám lòng dẫn. So sánh các phương pháp xác định quan hệ hình thái - lưu lượng gồm phương pháp Engelund, Van Rijn, Karim-Kennedy. Công thức tính suất chuyển cát đáy, sức tải cát lơ lửng, tổng lượng bùn cát và các giải pháp công trình chỉnh trị, bảo vệ bờ sông.
  • Chương 15, 16: Mô hình chất lượng nước: Phương trình vi phân cơ bản truyền chất một chiều mô tả quá trình đối lưu (advection) và phân tán (dispersion), điều kiện biên, thuật toán số và mô phỏng quá trình xâm nhập mặn ở vùng cửa sông.
[Hệ thống sông & Dòng chảy hở (Ch. 1)]

Kiến thức nền tảng được xây dựng

  • Các lý thuyết nền tảng (Fundamental theories): Định lý truyền tải Reynolds áp dụng cho thể tích kiểm tra (hệ thống bảo toàn khối lượng, động lượng, năng lượng); lý thuyết sóng dài trong lòng dẫn hở; lý thuyết dòng chảy rối (mô hình đẳng hướng Boussinesq, mô hình nhớt xoáy phi tuyến NKE, mô hình ứng suất Reynolds RSM); lý thuyết chuyển động bùn cát đáy và lơ lửng.
  • Các nguyên lý cốt lõi (Core principles): Nguyên lý cân bằng động lực trong sông bồi tích; nguyên lý bảo toàn vật chất trong phương trình liên tục; nguyên lý tương quan giữa lực quán tính và trọng lực (biểu thị qua số Froude); nguyên lý tương tác biên giữa dòng chảy và lòng dẫn di động.
  • Các khung phân tích chủ đạo (Essential frameworks): Khung phân loại hệ thống sông Schumm; hệ phương trình Saint-Venant 1D hoàn chỉnh; khung phân tích lượng trữ hình lăng trụ và hình nêm Muskingum; các mô hình quan hệ hình thái - lưu lượng (Engelund, Van Rijn, Karim-Kennedy).

Kỹ năng phát triển

  • Kỹ năng kỹ thuật (Technical skills): Thiết lập và giải số hệ phương trình đạo hàm riêng mô tả dòng chảy không ổn định bằng phương pháp sai phân hữu hạn; áp dụng thuật toán Runge-Kutta bậc 3 và thuật toán lặp trực tiếp để giải bài toán điều tiết lũ qua hồ chứa; tính toán lưu lượng tháo qua đập tràn tự do, cống xả tự do và cống xả ngập.
  • Kỹ năng phân tích (Analytical skills): Phân tích độ nhạy của các biến hình thái sông bồi tích theo quy mô thời gian; nhận diện hiện tượng vòng dây trong quan hệ mực nước - lưu lượng ($Z \sim Q$) do ảnh hưởng của nước vật; đánh giá cơ chế phân bố ứng suất và năng lượng rối trong mặt cắt ngang kênh dẫn.
  • Năng lực ứng dụng thực tiễn (Practical competencies): Xác định bộ thông số mô hình Muskingum ($K, x$) cho mạng lưới sông phân nhánh phức tạp dựa trên chuỗi số liệu thủy văn thực đo; tính toán thiết kế mặt cắt kênh ổn định không xói, không bồi; lập sơ đồ công trình chỉnh trị và đánh giá nguy cơ xâm nhập mặn vùng hạ lưu sông.

Phương pháp giảng dạy và học tập

Giáo trình định hướng phương pháp tiếp cận sư phạm theo mô hình diễn dịch kỹ thuật: đi từ việc thiết lập các phương trình vi phân chủ đạo dựa trên định luật bảo toàn vật lý, phân tích các giả thiết đơn giản hóa, chuyển đổi sang sơ đồ rời rạc hóa số học, đến triển khai thuật toán tính toán cụ thể cho các bài toán công trình thủy lợi, thủy điện và chỉnh trị sông ngòi.

Định luật bảo toàn vật lý (Khối lượng, Động lượng, Năng lượng)
      Phương trình vi phân liên tục (Saint-Venant)
       Rời rạc hóa số học (Sai phân hữu hạn)
Thuật toán giải số & Kiểm định dữ liệu thực địa (Sông Hồng, ĐBSCL)

Nội dung giáo trình được bổ trợ bằng các ví dụ mẫu có số liệu định lượng chi tiết:

  • Bài toán diễn toán lũ qua hồ chứa tĩnh: Minh họa thuật giải bán đồ giải của V.T. Chow thông qua bảng quan hệ phụ trợ $[\frac{2S}{\Delta t} + Q] \sim Q$, chia bước thời gian $\Delta t$, sử dụng công thức nội suy tuyến tính: $$y = y_1 + \frac{y_2 - y_1}{x_2 - x_1}(x - x_1)$$
  • Bài toán diễn toán lũ sông đơn bằng phương pháp Muskingum: Ví dụ mẫu với chuỗi số liệu lưu lượng gồm 20 thời đoạn, lưu lượng đỉnh đến $I_{\max} = 691\text{ m}^3/\text{s}$, các tham số $K = 2{,}3$, $x = 0{,}15$, $\Delta t = 1\text{h}$, kiểm tra điều kiện tổng hệ số $C_1 + C_2 + C_3 = 1{,}0$ và tính toán chi tiết từng bước thời gian.
  • Bài toán mạng lưới sông phân nhánh: Áp dụng cho hệ thống sông Hồng gồm 8 nút sông và 7 đoạn sông (từ Hòa Bình, Yên Bái, Hàm Yên, Chiêm Hóa, Tuyên Quang, Thác Bà, Vụ Quang, Việt Trì đến Sơn Tây), sử dụng chuỗi số liệu trận lũ thực đo từ ngày 12/VIII đến 26/VIII/1969 với thời khoảng 6 giờ để xác định bộ thông số $K$ và $x$ bằng phương pháp thử sai.

Hệ thống câu hỏi ôn tập được bố trí ở cuối mỗi chương tập trung vào việc kiểm tra bản chất vật lý của các phương trình, sự khác biệt giữa dòng chảy kênh nhân tạo và sông thiên nhiên, cũng như phạm vi áp dụng của từng phương pháp giải.

Phương pháp tự học khuyến nghị người học:

  1. Nắm vững việc thiết lập các phương trình vi phân cơ bản từ thể tích kiểm tra $\Delta x$.
  2. Thực hành lập bảng tính Excel hoặc lập trình giải các sơ đồ tính toán lặp, thuật toán Muskingum và sai phân Saint-Venant.
  3. Liên hệ đối chiếu các giả thiết tính toán với điều kiện địa hình, thủy lực thực tế của từng lưu vực sông.

Điểm nổi bật và cập nhật

Giáo trình phản ánh giai đoạn chuyển giao quan trọng trong nghiên cứu thủy lực sông ngòi tại Việt Nam, từ các phương pháp kinh nghiệm hoặc đồ giải truyền thống (như phương pháp đồ giải Pô-ta-pôp) sang các mô hình toán học giải tích và mô phỏng số hiện đại.

Các cập nhật lý thuyết và công nghệ nổi bật trong giáo trình bao gồm:

  • Cơ chế mô hình hóa dòng chảy rối tiên tiến: Giới thiệu các hạn chế của giả thiết độ nhớt xoáy tuyến tính Boussinesq trong dòng kênh hở (không mô tả được dòng thứ cấp, không nhạy với độ cong dọc, thiếu cơ chế ảnh hưởng của mặt thoáng tự do). Phân tích chi tiết mô hình nhớt xoáy phi tuyến (NKE) của Pope sử dụng đạo hàm Oldroyd và mô hình vận chuyển ứng suất Reynolds (RSM) gồm 5 thành phần khuếch tán và tái phân phối năng lượng rối.
  • Hệ thống hóa các công cụ mô hình toán: Tổng hợp và phân tích nguyên lý của các phần mềm mô phỏng thủy lực, thủy văn trong và ngoài nước:
    • Mô hình trong nước: Mô hình mạng lưới sông VRSAP (1978) của GS. Nguyễn Như Khuê; các mô hình KOD01 (1974) và KOD02 (1985) của GS. Nguyễn Ân Niên.
    • Mô hình quốc tế: SOGREAH (1964) ứng dụng cho lũ tràn đồng Đồng bằng sông Cửu Long; các bộ phần mềm MIKE 11 (1996), MIKE Flood (2000), MIKE 21C, FLO-2D, HEC-1, HEC-HMS, HEC-ResSim.
  • Gắn kết thực tiễn thủy lợi và phòng chống thiên tai: Giáo trình giải quyết trực tiếp các bài toán kỹ thuật cấp bách tại Việt Nam như kiểm soát lũ hệ thống sông Hồng và sông Cửu Long, tính toán diễn biến xói lở - bồi lắng lòng dẫn bồi tích, ổn định kênh dẫn và đánh giá mức độ lan truyền ô nhiễm, xâm nhập mặn ở các vùng hạ lưu chịu ảnh hưởng của chế độ bán nhật triều/nhật triều.

Đối tượng sử dụng giáo trình

Giáo trình được biên soạn chuyên biệt cho các nhóm đối tượng sau trong hệ thống đào tạo đại học và sau đại học ngành kỹ thuật tài nguyên nước:

Nhóm đối tượng Mục đích & Phạm vi sử dụng Yêu cầu kiến thức tiên quyết (Prerequisites)
Học viên cao học & Nghiên cứu sinh Tài liệu học tập chính cho các học phần Thủy lực sông ngòi nâng cao, Thủy lực công trình, Chỉnh trị sông; tài liệu tham khảo để thực hiện luận văn, luận án về mô hình hóa dòng chảy và bùn cát. Cơ học chất lỏng/Thủy lực cơ sở, Thủy văn công trình, Toán giải tích, Phương pháp số và Lập trình tính toán.
Sinh viên đại học năm cuối Tài liệu tham khảo chuyên đề cho sinh viên năm 4 - năm 5 các ngành Kỹ thuật tài nguyên nước, Kỹ thuật công trình thủy, Thủy văn học, Quản lý tài nguyên nước và Môi trường. Thủy lực đại cương, Toán cao cấp, Thủy văn đại cương.
Giảng viên & Cán bộ nghiên cứu Khung tham chiếu học thuật để xây dựng đề cương bài giảng, bài tập lớn, đồ án môn học và thiết lập các mô hình số phục vụ nghiên cứu khoa học. Năng lực chuyên sâu về cơ học dòng chảy hở và mô hình hóa thủy động lực học.
Kỹ sư tư vấn & Quản lý dự án Sổ tay lý thuyết phục vụ việc lựa chọn thông số đầu vào, kiểm tra điều kiện biên và đánh giá độ tin cậy của các kết quả tính toán thủy lực từ các phần mềm thương mại. Kiến thức thực tế về quy hoạch hệ thống thủy lợi, chỉnh trị sông và phòng chống thiên tai.

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

1. Giáo trình này phù hợp với đối tượng người học nào?

Giáo trình được biên soạn phục vụ học viên cao học, nghiên cứu sinh và kỹ sư thuộc các chuyên ngành Thủy văn học, Chỉnh trị sông và bờ biển, Phát triển nguồn nước, Kỹ thuật tài nguyên nước và Công trình thủy. Sinh viên đại học các năm cuối cũng có thể sử dụng làm tài liệu tham khảo cho đồ án tốt nghiệp và các chuyên đề nâng cao.

2. Người học cần chuẩn bị kiến thức nền tảng nào trước khi tiếp cận giáo trình?

Người học cần hoàn thành các học phần tiên quyết bao gồm: Cơ học chất lỏng/Thủy lực cơ sở (các phương trình Bernoulli, Euler, phương trình liên tục), Thủy văn công trình, Giải tích toán học (phương trình vi phân đạo hàm riêng, phép biến đổi Laplace) và Phương pháp số (sai phân hữu hạn, các thuật toán giải hệ phương trình phi tuyến).

3. Điểm khác biệt căn bản giữa giáo trình này và các tài liệu Thủy lực đại cương là gì?

Khác với Thủy lực đại cương vốn tập trung chủ yếu vào dòng chảy dừng, đều trong kênh nhân tạo có mặt cắt lăng trụ cố định, giáo trình Thủy lực sông ngòi tập trung giải quyết bài toán dòng không ổn định, không đều trong sông thiên nhiên với lòng dẫn bồi tích biến đổi, có xét đến các hiệu ứng phức tạp như chảy rối phi tuyến, dòng thứ cấp tại đoạn sông cong, vận chuyển bùn cát và lan truyền chất hòa tan.

4. Làm thế nào để tự học và làm chủ nội dung giáo trình hiệu quả?

Người học nên kết hợp việc đọc hiểu các giả thiết vật lý với việc tự lập trình (hoặc sử dụng bảng tính) để thực thi lại các thuật toán mẫu được trình bày trong sách, chẳng hạn như sơ đồ Runge-Kutta bậc 3 cho hồ chứa, phương pháp Muskingum cho đoạn sông và mạng sông Hồng trận lũ năm 1969. Sau đó, mở rộng đối chiếu lý thuyết với việc ứng dụng các phần mềm chuyên dụng như HEC-HMS, HEC-RAS hay MIKE 11.

5. Giáo trình có các tài liệu tham khảo bổ trợ tiêu biểu nào?

Giáo trình trích dẫn và liên kết với nhiều tài liệu kinh điển trong nước và quốc tế, bao gồm:

  • Applied Hydrology (Chow, V.T., Maidment, D.R., Mays, L.W., 1988/1998)
  • Kinematic Wave Modelling in Water Resources (Singh, V.P., 1996)
  • Open Channel Hydraulics (Sturm, T., 2001)
  • River Mechanics (Julien, P.Y., 2002)
  • Fluvial Processes in River Engineering (Chang, H.H., 2003)
  • Floodplain Processes (Anderson, M.G. et al., 1996)
  • River Training Techniques (Przedwojski, B., Blazejewski, R. et al., 1995)

Kết luận

Giáo trình Thủy lực sông ngòi – Tập I: Các nội dung lý thuyết cung cấp hệ thống tri thức hoàn chỉnh và chuẩn mực về cơ sở lý thuyết thủy động lực học dòng chảy trong sông, các phương pháp diễn toán dòng không ổn định từ mô hình thông số tập trung đến mô hình thông số phân bố, cơ chế vận chuyển bùn cát và mô hình hóa chất lượng nước.

Lộ trình học tập đề xuất cho người học là nghiên cứu tuần tự từ việc làm chủ các phương trình động lực học dòng chảy hở (Phần I, II), tiến đến tính toán tương tác lòng dẫn - bùn cát (Phần III) và hoàn thiện với bài toán lan truyền chất, xâm nhập mặn (Phần IV). Để tối ưu hóa quá trình tiếp thu, người học nên khai thác kết hợp các bộ phần mềm mô phỏng (HEC-HMS, FLO-2D, MIKE 11) cùng hệ thống tài liệu tham khảo quốc tế được chỉ dẫn chi tiết trong từng chương mục của giáo trình.