Tổng quan về luận án

Sự biến đổi sâu sắc của chế độ thủy văn - thủy lực trên hệ thống sông Hồng dưới tác động kép của biến đổi khí hậu toàn cầu và chuỗi hồ chứa thủy điện quy mô lớn ở thượng nguồn (như thủy điện Hòa Bình, Sơn La, Tuyên Quang và các hồ chứa trên lưu vực sông Đà phía Trung Quốc) đã làm suy giảm nghiêm trọng lượng bùn cát bồi đắp và gây hạ thấp mực nước mùa kiệt hạ lưu. Mực nước sông Hồng tại trạm Hà Nội liên tục suy giảm, ghi nhận các mốc lịch sử nghiêm trọng như năm 2010 chỉ đạt mức xấp xỉ $+0{,}1\text{ m}$, khiến hàng loạt công trình thủy lợi lấy nước tự chảy rơi vào tình trạng "treo đáy", mất hoàn toàn khả năng lấy nước trong các thời kỳ kiệt quệ.

Trong bối cảnh đó, sông Đáy – một phân lưu quan trọng dài $240\text{ km}$ của sông Hồng với diện tích lưu vực hơn $5.800\text{ km}^2$ trải qua 6 tỉnh thành (Hà Nội, Hòa Bình, Hà Nam, Nam Định, Ninh Bình) – đã bị cô lập thủy lực tự nhiên kể từ khi xây dựng đập Đáy năm 1937 để phục vụ mục tiêu phân lũ phòng thủ cho thủ đô. Hệ quả sinh thái và thủy văn là "sông Đáy bị chết dần, hầu như không còn dòng chảy trên đoạn 23km từ đập Đáy đến cầu Mai Lĩnh", gây suy thoái môi trường nước và khan hiếm nguồn nước tưới trầm trọng. Nhằm khôi phục hệ sinh thái dòng chảy và chuyển đổi chức năng phòng chống thiên tai theo Nghị định 04/2011/NĐ-CP (xóa bỏ khu phân chậm lũ, cấp nước thường xuyên mùa kiệt $30\div 100\text{ m}^3/\text{s}$, mùa lũ $500\div 800\text{ m}^3/\text{s}$ và phân lũ tối đa $2.500\text{ m}^3/\text{s}$), cụm công trình đầu mối Cẩm Đình – Hiệp Thuận đã được đầu tư xây dựng. Tuy nhiên, sau khi hoàn thành năm 2004, cụm công trình này bộc lộ khiếm khuyết kỹ thuật nghiêm trọng khi "lòng dẫn đã bị đất cát bồi lắng tương đối lớn trung bình so với thiết kế là 1,03m", làm tê liệt hoàn toàn chức năng cấp nước tự chảy mùa kiệt.

Khoảng trống học thuật (Research Gap) lớn nhất đặt ra là: Các nghiên cứu thủy lực và chỉnh trị sông truyền thống ở Việt Nam cũng như trên thế giới chủ yếu tiếp cận theo mô hình một chiều (1D) hoặc hai chiều phẳng (2D) tĩnh/động, hoàn toàn bỏ qua bản chất thủy động lực học ba chiều (3D) của dòng phân tầng mật độ (stratified/density flow) và quy luật phân bố dị trọng bùn cát theo phương đứng tại khu vực ngã ba phân lưu và cửa vào công trình lấy nước. Chưa có một cơ sở khoa học đa tiêu chí toàn diện nào xác lập được tương quan giữa hình thái khúc sông cong, trường vận tốc thứ sinh, góc lệch cửa lấy nước và cơ chế bồi lắng để xác định vị trí đặt công trình tối ưu phục vụ đồng thời hai mục tiêu mâu thuẫn: tối đa hóa lưu lượng lấy nước mùa kiệt và giảm thiểu xâm nhập bùn cát bồi lấp khi thoát lũ.

Luận án tiến sĩ kỹ thuật chuyên ngành Thủy văn học (Mã số: 9440224) của tác giả Trần Khắc Thạc, dưới sự hướng dẫn khoa học của GS.TS. Phạm Thị Hương Lan và PGS.TS. Tan Soon Keat tại Trường Đại học Thủy lợi, tập trung giải quyết 3 câu hỏi nghiên cứu cốt lõi:

  1. RQ1: Động lực học dòng chảy 3 chiều và cơ chế phân tầng dị trọng bùn cát tại khu vực cửa lấy nước biến đổi như thế nào dưới các kịch bản thủy văn cực đoan mùa kiệt và mùa lũ?
  2. RQ2: Những tiêu chí thủy lực, hình thái và bùn cát nào đóng vai trò quyết định trong việc định vị không gian cửa lấy nước trên đoạn sông cong phân lạch?
  3. RQ3: Đâu là vị trí và giải pháp công trình/phi công trình tối ưu trên đoạn sông Hồng qua cửa vào sông Đáy để đảm bảo khả năng tiếp nhận lưu lượng $30\div 100\text{ m}^3/\text{s}$ trong mùa kiệt và thích ứng với yêu cầu xả lũ đến $2.500\text{ m}^3/\text{s}$?

Nghiên cứu được thiết lập trên phạm vi địa lý $27\text{ km}$ dòng chính sông Hồng từ $\text{Km}30$ đến $\text{Km}47+500$ (tương ứng đoạn từ Sơn Tây - Cẩm Đình đến Hát Môn), sử dụng chuỗi số liệu thực đo thủy văn, địa hình chi tiết và các trận lũ/kiệt lịch sử (1996, 2001, 2002, 2003–2004). Kết quả nghiên cứu tạo ra bước đột phá định lượng, đặt nền móng lý thuyết và thực tiễn vững chắc cho công tác tái cấu trúc hạ tầng thủy lợi lưu vực sông Hồng - sông Đáy.

Literature Review và Positioning

Lịch sử nghiên cứu về công trình lấy nước (CTLN) trên các dòng sông bồi tích trên thế giới đã hình thành nhiều trường phái học thuật với những góc nhìn đa dạng:

  • Phân loại và hình học công trình lấy nước: Avery (1989) và GTZ (1989) đã hệ thống hóa CTLN thành dạng lấy nước bên hông (lateral intake) và lấy nước chính diện/lấy nước đáy (Tyrolean intake), chỉ ra rằng góc hợp dòng giữa trục cống và dòng chính quyết định tổn thất năng lượng thủy lực. Abbasi (2003) phân tích sâu cấu trúc đập dâng điều tiết mực nước thượng lưu nhưng chỉ giới hạn phạm vi áp dụng cho các sông miền núi có độ dốc đáy lớn ($1%\div 10%$).
  • Tranh luận học thuật về góc lệch dòng và vị trí tuyến công trình: Tồn tại sự bất đồng lớn giữa các kết quả nghiên cứu thực nghiệm. Ali Akbar Salehi (2008) và Alireza Masjedi (2011) dựa trên mô hình vật lý máng dẫn cho rằng góc lấy nước tối ưu nên dao động từ $45^\circ\div 90^\circ$ tại bờ lõm đoạn sông cong. Ngược lại, Mehdi Karami Moghadam (2010) khi nghiên cứu chế độ thủy động lực cống lấy nước trên sông Ohio (Hoa Kỳ) lại khẳng định góc lệch $30^\circ$ mang lại hiệu suất tháo lưu lượng lớn nhất và nồng độ bùn cát nhỏ nhất, dù nghiên cứu này bị giới hạn ở dải số Froude $Fr = 0{,}35\div 0{,}40$ và bỏ qua chế độ kiệt. Afrouza (2015) và Montaseri (2015) qua mô phỏng 3D dòng chảy quanh cống bằng mô hình Reynolds Stress (RSM) lại đề xuất góc tối ưu hẹp hơn rất nhiều, chỉ từ $17^\circ\div 20^\circ$ hoặc $15^\circ\div 30^\circ$ để triệt tiêu vùng xoáy phân ly và hạn chế bùn cát.
  • Nghiên cứu kiểm soát bùn cát và dòng thứ sinh: Azade Jamshidi và đồng sự (2016) chứng minh việc sử dụng các tấm lái dòng ngập nước (submerged vanes) đặt lệch góc $15^\circ\div 30^\circ$ có thể làm tăng lưu lượng vào cống từ $30{,}72%\div 36{,}34%$ thông qua can thiệp vào trường vận tốc xoắn đáy. Tại châu Âu, các tổng kết chỉnh trị sông Rhine (2006) sau các trận lũ lịch sử 1993 và 1995 nhấn mạnh sự kết hợp giữa đê phân lũ, khôi phục bãi ngập và kè mỏ hàn bảo vệ bờ sinh thái.
                  ┌────────────────────────────────────────────────────────┐
                  │              CÁC NHÁNH NGHIÊN CỨU QUỐC TẾ              │
                  └──────────────────────────┬─────────────────────────────┘
                                             │
         ┌───────────────────────────────────┼────────────────────────────────────┐
         ▼                                   ▼                                    ▼
┌──────────────────┐               ┌───────────────────┐                ┌──────────────────┐
│  Mô hình vật lý  │               │   Mô hình 2D      │                │  Mô hình số 3D   │
│  máng dẫn tĩnh   │               │   lưới thô        │                │  cục bộ vi mô    │
├──────────────────┤               ├───────────────────┤                ├──────────────────┤
│ Salehi (2008)    │               │ Elmoustafa (2011) │                │ Neary et al.(1999│
│ Masjedi (2011)   │               │ RMA2 (20x20m)     │                │ Montaseri (2015) │
│ Moghadam (2010)  │               │ Shahidan (2012)   │                │ Bỏ qua biến động │
│ Không xét dòng   │               │ Không mô tả dòng  │                │ lòng dẫn mùa kiệt│
│ kiệt sâu         │               │ phân tầng đứng    │                │ quy mô lưu vực   │
└────────┬─────────┘               └─────────┬─────────┘                └────────┬─────────┘
         │                                   │                                   │
         └───────────────────────────────────┼───────────────────────────────────┘
                                             ▼
                  ┌────────────────────────────────────────────────────────┐
                  │            KHOẢNG TRỐNG VÀ VỊ THẾ CỦA LUẬN ÁN          │
                  ├────────────────────────────────────────────────────────┤
                  │ • Tích hợp lý thuyết dòng phân tầng dị trọng 3D        │
                  │ • Mô hình thủy động lực - hình thái học MIKE 3 FM      │
                  │ • Thiết lập bộ tiêu chí liên hoàn: Thủy lực - Bùn cát  │
                  │ • Giải quyết đồng thời: Cấp nước mùa kiệt & Thoát lũ   │
                  └────────────────────────────────────────────────────────┘

Tại Việt Nam, các nghiên cứu trước đây tại cụm công trình đầu mối sông Hồng - sông Đáy còn bộc lộ nhiều khoảng trống phương pháp luận:

  1. Nghiên cứu của Viện Khoa học Thủy lợi Việt Nam (2001, 2005) và Viện Quy hoạch Thủy lợi (2009) chủ yếu sử dụng mô hình 1D (MIKE 11) hoặc mô hình vật lý chính thái đơn lẻ (Nguyễn Tuấn Anh, 2002), chỉ tập trung vào khả năng tháo lũ tổng thể mà chưa đánh giá được cơ chế bồi lấp cục bộ kênh dẫn vào mùa kiệt.
  2. Các nghiên cứu ứng dụng mô hình 2D (Phạm Đình, 2010; Phạm Thị Hương Lan, 2012) trên nền MIKE 21 FM-ST đã phân tích biến động lòng dẫn nhưng bị hạn chế bởi giả thiết vận tốc trung bình theo chiều đứng, không thể mô tả được sự phân lớp của nồng độ phù sa lơ lửng và bùn cát đáy khi dòng chảy rẽ nhánh.
  3. Nghiên cứu 3D của Phạm Đức Thắng (2010) tại cống Xuân Quan chỉ khảo sát giai đoạn mùa lũ khi cống đóng hoàn toàn, không phản ánh được tương tác động lực học trong quá trình mở cống lấy nước mùa kiệt.

So sánh trực tiếp với hai nghiên cứu quốc tế điển hình: Nghiên cứu của Elmoustafa (2011) trên sông Nile (Assuit, Ai Cập) sử dụng mô hình 2D RMA2 với lưới ô vuông thô $20\times 20\text{ m}$ nên làm mờ hoàn toàn cấu trúc hình học cửa cống; trong khi nghiên cứu của Shahidan và cs. (2012) trên sông Ijok (Malaysia) chỉ ghép nối lỏng lẻo HEC-RAS 1D với CCHE2D. Luận án của Trần Khắc Thạc đã vượt qua các giới hạn này bằng cách thiết lập mô hình 3D phi cấu trúc (MIKE 3 FM) có độ phân giải cao, tích hợp lý thuyết dòng phân tầng mật độ để lượng hóa đồng thời dòng chảy và bùn cát trên một đoạn sông tự nhiên phức tạp dài $27\text{ km}$.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng và làm sâu sắc thêm lý thuyết thủy lực sông ngòi (fluvial hydraulics) và động lực học lòng dẫn thông qua các đóng góp lý thuyết chuyên biệt:

  1. Làm sáng tỏ cơ chế phân tầng dị trọng tại cửa lấy nước bên sông: Tác giả đã tích hợp hệ phương trình Navier-Stokes trung bình hóa theo thời gian (RANS) với phương trình truyền tải - khuếch tán bùn cát 3 chiều: $$\frac{\partial \rho}{\partial t} + \frac{\partial (\rho u)}{\partial x} + \frac{\partial (\rho v)}{\partial y} + \frac{\partial (\rho w)}{\partial z} = 0$$ Nghiên cứu chứng minh rằng gradient mật độ $\partial \rho / \partial z$ sinh ra do sự phân tầng nồng độ phù sa lơ lửng theo phương đứng làm biến dạng cấu trúc dòng chảy xoắn thứ sinh tại bờ lõm, tạo nên một "nêm bùn cát" đậm đặc ở tầng đáy có xu hướng bị hút mạnh vào cửa lấy nước khi tỷ lệ phân lưu $Q_{\text{intake}}/Q_{\text{river}}$ thay đổi.
  2. Xác lập mô hình tương quan phi tuyến đa biến: Luận án đề xuất và kiểm chứng thực nghiệm mô hình tương quan giữa 5 tham số thủy động lực và hình thái: $$Q_c = f\left(H_{\text{SH}}, \alpha, L_{\text{cong}}, B_{\text{kênh}}, S_v(z)\right)$$ Trong đó $Q_c$ là lưu lượng lấy vào cống, $H_{\text{SH}}$ là mực nước sông Hồng tại cửa cống, $\alpha$ là góc lấy nước, $L_{\text{cong}}$ là chiều dài đoạn sông cong phát triển, $B_{\text{kênh}}$ là chiều rộng kênh dẫn, và $S_v(z)$ là hàm phân bố độ đục theo chiều sâu.
       BỜ LÕM SÔNG HỒNG                                    CỬA LẤY NƯỚC CẨM ĐÌNH
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~                     ┌────────────────────────────────
                               \   Dòng mặt         │
    Dòng chảy chính (Q_river)   \  (Nồng độ cát thấp)│  Tầng mặt: Lấy nước tưới
   ══════════════════════════>   \───────┐          │  (Chất lượng nước cao)
                                  \      │ α = 30°  │
                                   \     ▼          │
   -----------------------------    \ ════════════> │  Ngưỡng tràn trên (+9.5m)
    Dòng xoắn thứ sinh đáy          \               ├────────────────────────────────
    (Vận chuyển phù sa đậm đặc)      \   Dòng đáy   │  Tầng đáy: Bị bồi lắng
   ───────────────────────────>       \  (Dị trọng) │  Ngưỡng cống kiệt (+3.0m)
                                       \─────►      │  (Bồi lấp trung bình 1.03m)
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~            └────────────────────────────────

Khung phân tích độc đáo

Khung phân tích của luận án được xây dựng dựa trên sự hợp nhất giữa 3 trường phái lý thuyết: Lý thuyết dòng thứ sinh khúc sông cong (Rozovskii, 1957), Lý thuyết vận chuyển bùn cát không cân bằng (van Rijn, 1984), và Lý thuyết dòng dị trọng phân tầng (Turner, 1973). Khung phân tích này định hình nên Bộ tiêu chí 4 nhóm để lựa chọn vị trí cửa lấy nước:

  • Tiêu chí 1 (Thủy lực dòng dẫn): Vị trí đặt cống phải nằm ở vùng có trường vận tốc dòng chính ổn định, duy trì $v \ge v_{\text{không lắng}}$ ($0{,}6\div 0{,}8\text{ m/s}$) tại cửa vào trong suốt thời đoạn kiệt để hạn chế lắng đọng phù sa.
  • Tiêu chí 2 (Hình thái và động lực lòng sông): Nằm tại cung bờ lõm của đoạn sông cong tự nhiên ổn định, nơi dòng chảy mặt hướng vào bờ và dòng đáy mang bùn cát bị dòng xoắn đẩy ra ngoài tim sông.
  • Tiêu chí 3 (Góc hợp dòng $\alpha$): Tuyến lấy nước tạo với trục dòng chính một góc tiếp cận hợp lý ($\alpha \approx 30^\circ\div 45^\circ$) nhằm triệt tiêu vùng xoáy phân ly tĩnh tại mũi cống.
  • Tiêu chí 4 (Hàm lượng bùn cát xâm nhập): Giảm thiểu tối đa lưu lượng bùn cát lơ lửng và bùn cát đáy đi vào kênh dẫn trong cả hai chu kỳ lũ và kiệt.

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

Thiết kế nghiên cứu

Luận án vận dụng quan điểm triết học thực chứng (Positivism) kết hợp chủ nghĩa duy thực phê phán (Critical Realism) trong khoa học kỹ thuật thủy lợi. Thiết kế nghiên cứu tích hợp phương pháp tiếp cận hệ thống đa cấp độ (Multi-level Systemic Approach):

  • Cấp vĩ mô: Đánh giá thủy lực toàn hệ thống mạng lưới sông Hồng - sông Thái Bình bằng mô hình 1D (MIKE 11).
  • Cấp trung vi mô: Mô phỏng chi tiết trường dòng chảy 3 chiều và biến động xói bồi lòng dẫn trên miền tính $27\text{ km}$ sông Hồng qua cửa sông Đáy bằng mô hình 3D phần tử hữu hạn (MIKE 3 FM).

Quy trình nghiên cứu rigorous

Quy trình tính toán số trị tuân thủ nghiêm ngặt các chuẩn mực quốc tế về mô hình hóa động lực học thủy văn:

  1. Rời rạc hóa không gian: Sử dụng lưới tam giác phi cấu trúc (Flexible Mesh) linh hoạt với $24.850$ phần tử tam giác và tứ giác; độ phân giải lưới được làm mịn tối đa tại khu vực cửa cống Cẩm Đình, Hát Môn và các công trình điều tiết ($\Delta x \approx 2\div 5\text{ m}$), giãn dần ra ngoài dòng chính ($\Delta x \approx 50\div 100\text{ m}$). Theo phương đứng, trường dòng chảy được chia thành 5 lớp $\sigma$ (sigma-layers) tương thích với địa hình đáy thay đổi linh hoạt theo biến trình mực nước.
  2. Điều kiện biên và CSDL: Biên trên là chuỗi lưu lượng thực đo tại trạm thủy văn Sơn Tây; biên dưới là chuỗi mực nước tại Hà Nội và Thượng Cát. Nồng độ bùn cát biên được thiết lập từ đường quan hệ lưu lượng - nồng độ bùn cát thực đo nhiều năm.
  3. Hiệu chỉnh và kiểm định (Calibration & Validation): Mô hình được hiệu chỉnh với chuỗi số liệu lũ năm 1996 và kiệt năm 2001; kiểm định độc lập với trận lũ năm 2002 và mùa kiệt 2003–2004 tại các trạm Thượng Cát, Hà Nội và mặt cắt kiểm soát MC-1.
       QUY TRÌNH MÔ HÌNH HÓA THỦY ĐỘNG LỰC HỌC 3D (MIKE 3 FM)
 ┌───────────────────────────────────────────────────────────────────┐
 │                CSDL ĐỊA HÌNH & THỦY VĂN THỰC ĐO                   │
 │ (Địa hình 27km Km30-Km47+500, Mặt cắt cống, Trạm Sơn Tây, Hà Nội)  │
 └─────────────────────────────────┬─────────────────────────────────┘
                                   │
                                   ▼
 ┌───────────────────────────────────────────────────────────────────┐
 │                 THIẾT LẬP MÔ HÌNH MẠNG LƯỚI 1D                    │
 │         (MIKE 11 HD & ST toàn bộ hệ thống sông Hồng - Thái Bình)   │
 └─────────────────────────────────┬─────────────────────────────────┘
                                   │ (Cung cấp điều kiện biên)
                                   ▼
 ┌───────────────────────────────────────────────────────────────────┐
 │                 THIẾT LẬP MÔ HÌNH TOÁN 3 CHIỀU                    │
 │               (MIKE 3 FM: 24.850 phần tử, 5 lớp Sigma)            │
 └─────────────────────────────────┬─────────────────────────────────┘
                                   │
         ┌─────────────────────────┴─────────────────────────┐
         ▼                                                   ▼
 ┌──────────────────────────────────┐    ┌───────────────────────────────────┐
 │     GIAI ĐOẠN HIỆU CHỈNH         │    │       GIAI ĐOẠN KIỂM ĐỊNH         │
 │ • Mùa kiệt: 2001                 │    │ • Mùa kiệt: 2003 - 2004           │
 │ • Mùa lũ: Trận lũ 1996           │    │ • Mùa lũ: Trận lũ 2002            │
 │ • Hệ số Manning: n = 0.022-0.028 │    │ • Đánh giá tại trạm Thượng Cát,   │
 │ • Tiêu chuẩn Nash-Sutcliffe >0.88│    │   Hà Nội và mặt cắt MC-1          │
 └──────────────────────────────────┘    └───────────────────────────────────┘
                                   │
                                   ▼
 ┌───────────────────────────────────────────────────────────────────┐
 │           MÔ PHỎNG CÁC KỊCH BẢN VỊ TRÍ & QUY HOẠCH CTLN           │
 │ (Kịch bản 1: Cẩm Đình hiện trạng | Kịch bản 2: Vị trí mới hạ lưu  │
 │  600m | Kịch bản 3: Vị trí Vân Cốc | Kịch bản 4: Tuyến Hát Môn)   │
 └───────────────────────────────────────────────────────────────────┘

Data và phân tích

Độ tin cậy của mô hình thủy lực và bùn cát đạt mức rất cao:

  • Hệ số nhám Manning $n$ hiệu chỉnh trên đoạn sông dao động từ $0{,}022\div 0{,}028\text{ s/m}^{1/3}$; hệ số phân tán bùn cát nằm trong giới hạn vật lý $D_x = D_y = 1{,}5\div 3{,}0\text{ m}^2/\text{s}$.
  • Chỉ số đánh giá hiệu quả mô hình Nash-Sutcliffe ($NSE$) đối với mực nước tại trạm Hà Nội và Thượng Cát đạt $NSE \ge 0{,}91$; đối với lưu lượng đạt $NSE \ge 0{,}88$.
  • Sai số tương đối về lưu lượng bùn cát lơ lửng tại mặt cắt kiểm tra $\text{MC}-1$ qua các năm 1996 và 2002 đều $< 12%$, đảm bảo độ tin cậy tuyệt đối để chạy các kịch bản dự báo.

Phát hiện đột phá và implications

Những phát hiện then chốt

  1. Giải mã nguyên nhân bồi lấp cống Cẩm Đình hiện trạng: Kết quả mô phỏng 3D chỉ ra rằng cống Cẩm Đình hiện tại đặt ở vị trí thoái trào của khúc cong, góc lấy nước quá tù so với hướng dòng chính mùa lũ, tạo ra một vùng xoáy tĩnh (recirculation zone) rộng lớn ngay trước cửa vào cống. Vận tốc dòng chảy giảm đột ngột xuống $v < 0{,}15\text{ m/s}$, làm cho toàn bộ lượng bùn cát lơ lửng hạt mịn ($d_{50} = 0{,}029\text{ mm}$) bị lắng đọng với tốc độ bồi tụ thực tế đạt $1{,}03\text{ m}$ dọc theo đoạn kênh dẫn thượng lưu dài $700\text{ m}$, biến cống thành một "bẫy cát" khổng lồ.
  2. Khả năng tháo kiệt thực tế của cống hiện trạng: Trong suốt 5 tháng mùa kiệt (từ 01/12 đến 30/4 năm 2003–2004), mực nước sông Hồng hạ thấp liên tục khiến cống Cẩm Đình chỉ lấy được lưu lượng trung bình $Q \le 5\div 12\text{ m}^3/\text{s}$, hoàn toàn không đạt lưu lượng thiết kế $Q_{\text{TK}} = 36{,}24\text{ m}^3/\text{s}$ và không đáp ứng yêu cầu tối thiểu $30\div 100\text{ m}^3/\text{s}$ theo Nghị định 04/2011/NĐ-CP.
  3. Phát hiện vị trí vàng tại Kịch bản 2 (Cẩm Đình mới - dịch hạ lưu 600m): Khi dịch chuyển vị trí cửa lấy nước về phía hạ lưu $600\text{ m}$ so với cống cũ (tương ứng vị trí 2):
    • Đặt cửa cống đúng vào tâm động lực của dòng chủ lưu áp bờ lõm trong mùa kiệt.
    • Vận tốc trước cửa lấy nước luôn duy trì ở ngưỡng tối ưu $v = 0{,}65\div 0{,}85\text{ m/s}$, ngăn chặn hoàn toàn hiện tượng lắng đọng bùn cát trước cống.
    • Lưu lượng lấy nước mùa kiệt tăng vọt lên mức $35\div 48\text{ m}^3/\text{s}$ ở chế độ tự chảy tự nhiên và đạt $100\text{ m}^3/\text{s}$ khi kết hợp cải tạo nạo vét lòng dẫn sông Đáy.
  4. Quy luật vận chuyển bùn cát trong điều kiện phân lũ cực hạn ($Q = 450\div 2.500\text{ m}^3/\text{s}$): Tại vị trí mới, nhờ hình học dòng dẫn hợp lý ($\alpha = 30^\circ$), bùn cát tầng đáy bị lực ly tâm đẩy trôi về phía bãi bồi đối diện, hàm lượng bùn cát xâm nhập vào cống giảm $42{,}6%$ so với vị trí Cẩm Đình cũ.
Kịch bản vị trí công trình Vận tốc cửa vào mùa kiệt ($v$) Khả năng cấp nước kiệt ($Q_{\text{kiệt}}$) Mức độ bồi lắng kênh dẫn sau 1 chu kỳ năm Khả năng đáp ứng NĐ 04/2011/NĐ-CP
Vị trí 1: Cẩm Đình cũ (Hiện trạng) $0{,}10\div 0{,}18\text{ m/s}$ $5\div 12\text{ m}^3/\text{s}$ (Rất kém) Bồi tụ dày $+1{,}03\text{ m}$ Thất bại hoàn toàn
Vị trí 2: Cẩm Đình mới (Hạ lưu 600m) $\mathbf{0{,}65\div 0{,}85\text{ m/s}}$ $\mathbf{35\div 100\text{ m}^3/\text{s}}$ (Rất tốt) Tự xói sạch, bồi lắng $< 0{,}15\text{ m}$ Đạt 100% tiêu chí
Vị trí 3: Tuyến Vân Cốc $0{,}30\div 0{,}45\text{ m/s}$ $18\div 25\text{ m}^3/\text{s}$ (Trung bình) Bồi lắng cục bộ $+0{,}55\text{ m}$ Đạt một phần
Vị trí 4: Tuyến Hát Môn (Cửa Đáy cũ) $0{,}25\div 0{,}40\text{ m/s}$ $15\div 22\text{ m}^3/\text{s}$ (Thấp) Xói lở bờ đối diện, bồi lòng Không khả thi 
    SO SÁNH HIỆU QUẢ CỦA 2 VỊ TRÍ CÔNG TRÌNH (CẨM ĐÌNH CŨ VS CẨM ĐÌNH MỚI)
 ┌────────────────────────────────────────┬────────────────────────────────────────┐
 │       VỊ TRÍ 1: CẨM ĐÌNH CŨ            │       VỊ TRÍ 2: CẨM ĐÌNH MỚI           │
 │            (HIỆN TRẠNG)                │          (HẠ LƯU 600M)                 │
 ├────────────────────────────────────────┼────────────────────────────────────────┤
 │ • Vận tốc cửa vào: 0.10 - 0.18 m/s     │ • Vận tốc cửa vào: 0.65 - 0.85 m/s     │
 │ • Lưu lượng kiệt: 5 - 12 m3/s          │ • Lưu lượng kiệt: 35 - 100 m3/s        │
 │ • Bồi lắng kênh: +1.03 m (Tê liệt)     │ • Bồi lắng kênh: < 0.15 m (Tự xói sạch)│
 │ • Góc lấy nước: Tù, tạo xoáy tĩnh      │ • Góc lấy nước: 30 độ (Khớp dòng mặt)  │
 └────────────────────────────────────────┴────────────────────────────────────────┘

Implications đa chiều

  • Ý nghĩa lý thuyết: Xác lập khung lý thuyết hoàn chỉnh về tương tác 3 chiều giữa dòng chảy phân tầng mật độ và hình thái sông cong bồi tích trong việc định vị công trình phân lưu.
  • Đổi mới phương pháp luận: Chuẩn hóa quy trình ứng dụng mô hình số trị 3D phi cấu trúc (MIKE 3 FM) trong thẩm định, thiết kế quy hoạch hạ tầng đầu mối sông ngòi tại Việt Nam.
  • Ứng dụng thực tiễn: Cung cấp tọa độ và thông số hình học chuẩn xác để các cơ quan quản lý thực hiện nạo vét, chuyển dời hoặc xây dựng cống phân lũ mới thay thế cụm Cẩm Đình hiện nay.
  • Khuyến nghị chính sách: Kiến nghị Bộ Nông nghiệp & PTNT và UBND TP. Hà Nội điều chỉnh quy trình vận hành liên hồ chứa thượng nguồn, kết hợp quy hoạch chỉnh trị sông Đáy theo phương án nạo vét 2 cấp mặt cắt ($B = 150\text{ m}$, đáy $+1{,}0\text{ m}$ tại Đập Đáy đến $-2{,}5\text{ m}$ tại Ba Thá) để đảm bảo đồng bộ năng lực chuyển nước.

Limitations và Future Research

Mặc dù đạt được những bước tiến vượt bậc, luận án cũng thẳng thắn thừa nhận các giới hạn khoa học nội tại:

  1. Giới hạn trường dữ liệu địa hình đáy động: Mô hình toán 3D sử dụng dữ liệu đo đạc địa hình đáy tại các thời điểm rời rạc; chưa tích hợp module hình thái tự điều chỉnh tức thời (morphological dynamic feedback) ở bước thời gian siêu nhỏ trong các trận lũ cực đoan.
  2. Biên độ biến đổi khí hậu chưa bao quát kịch bản vỡ đập/xả lũ khẩn cấp: Nghiên cứu tập trung vào các kịch bản thủy văn lịch sử và tần suất thiết kế $P = 0{,}2%$, chưa mô phỏng kịch bản rủi ro cực đoan khi chuỗi hồ chứa thượng nguồn đồng loạt xả lũ khẩn cấp vượt thiết kế kết hợp bão biển dâng cao triều cường tại cửa Đáy (Như Tân).
  3. Thành phần bùn cát đáy: Số liệu bùn cát chủ yếu tập trung vào cấp hạt bùn cát lơ lửng ($d_{50} = 0{,}029\text{ mm}$), tỷ lệ bùn cát đáy di chuyển nhảy cóc (bed-load) được ước tính qua công thức bán thực nghiệm mà chưa có thiết bị đo sóng âm liên tục tại hiện trường.

Các hướng nghiên cứu tiếp nối trong giai đoạn 5–10 năm tới bao gồm:

  • Tích hợp trí tuệ nhân tạo (Machine Learning/Physics-Informed Neural Networks - PINNs) với mô hình MIKE 3 FM để dự báo theo thời gian thực (Real-time Forecasting) diễn biến bồi lấp cửa cống theo từng giờ xả nước của thủy điện Hòa Bình.
  • Nghiên cứu giải pháp công trình ngầm thông minh: Thiết kế hệ thống cánh lái dòng ngập nước tự động điều chỉnh góc nghiêng (Automated Smart Vanes) đặt tại đáy sông Hồng trước cửa lấy nước Cẩm Đình mới.
  • Mở rộng phạm vi đánh giá chất lượng nước tổng thể (WQI) và đa dạng sinh học thủy sinh khi dòng nước ngọt sông Hồng được chuyển tải thường xuyên vào lòng sông Đáy.

Tác động và ảnh hưởng

Nghiên cứu của luận án mang lại giá trị tác động sâu rộng trên nhiều bình diện:

  • Tác động học thuật: Mở ra tiền lệ ứng dụng mô hình động lực học 3D chi tiết cho các bài toán phân lưu bồi tích tại Việt Nam; tạo nguồn trích dẫn học thuật tin cậy cho các nghiên cứu về kỹ thuật sông ngòi, chỉnh trị và quản lý lưu vực sông quốc tế.
  • Chuyển đổi ngành thủy lợi - nông nghiệp: Trực tiếp giải quyết bài toán cấp nước tưới tự nguyện cho hơn $10.150\text{ ha}$ đất canh tác nông nghiệp công nghệ cao thuộc các huyện Phúc Thọ, Đan Phượng, Hoài Đức, Quốc Oai và lưu vực sông Tích - sông Đáy.
  • Hiệu quả kinh tế - xã hội: Tiết kiệm hàng chục tỷ đồng ngân sách nhà nước hàng năm chi phí nạo vét bùn cát khơi thông kênh dẫn Cẩm Đình – Hiệp Thuận; khôi phục cảnh quan môi trường sinh thái đô thị phía Tây Hà Nội, biến dòng sông Đáy "chết" thành trục không gian xanh sinh thái bền vững.
  • Tầm ảnh hưởng quốc tế: Phương pháp luận xác định vị trí cửa lấy nước dựa trên mô hình 3D dòng phân tầng có thể chuyển giao áp dụng trực tiếp cho các quốc gia đang phát triển có hệ thống sông bồi tích tương đồng như lưu vực sông Mê Kông, sông Ayeyarwady (Myanmar), sông Hằng (Ấn Độ).

Đối tượng hưởng lợi

  • Nghiên cứu sinh, Học viên cao học chuyên ngành Kỹ thuật tài nguyên nước, Thủy văn: Tiếp cận một khung phương pháp luận nghiên cứu mẫu mực, kết hợp nhuần nhuyễn giữa cơ sở lý thuyết RANS 3D, công cụ mô hình toán hiện đại (MIKE 3 FM) và bài toán thực tế phức tạp.
  • Các Giáo sư, Nhà khoa học đầu ngành: Nắm bắt các quy luật tương quan mới giữa trường dòng chảy xoắn thứ sinh, sự phân tầng dị trọng bùn cát và góc lệch cửa lấy nước trên sông cong.
  • Kỹ sư tư vấn thiết kế công trình thủy lợi (HEC, VIVASO...): Bộ tiêu chí kỹ thuật 4 nhóm và các công thức tương quan định lượng là cẩm nang trực tiếp để tính toán, định tuyến các cống lấy nước dưới đê, trạm bơm ven sông (như cống Liên Mạc, Xuân Quan, trạm bơm Phù Sa, Đan Hoài, Ấp Bắc).
  • Cơ quan quản lý nhà nước (Bộ Nông nghiệp & PTNT, Cục Quản lý Đê điều & PCTT, Viện Quy hoạch Thủy lợi): Sở hữu căn cứ khoa học chính xác để phê duyệt dự án đầu tư nâng cấp cụm đầu mối phân lũ sông Đáy và điều chỉnh quy hoạch phòng chống lũ lưu vực sông Hồng.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc tích hợp thành công lý thuyết dòng phân tầng mật độ/dị trọng (Stratified/Density Flow Theory của Turner) vào hệ phương trình thủy động lực học sông ngòi 3 chiều (RANS kết hợp mô hình rối $k$-$\varepsilon$ RNG). Luận án đã phá vỡ giả định truyền thống xem trường nồng độ bùn cát là đồng nhất theo phương đứng tại cửa lấy nước, chứng minh rằng sự biến thiên gradient mật độ $\partial \rho / \partial z$ ở khúc sông cong thúc đẩy hiện tượng bồi lắng chọn lọc tầng đáy, làm cơ sở lý giải vì sao các cống lấy nước đặt sai vị trí động lực đều bị bồi lấp nghiêm trọng.

2. Điểm đổi mới phương pháp luận khi so sánh với ít nhất 2 nghiên cứu quốc tế trước đây?

So với nghiên cứu của Elmoustafa (2011) trên sông Nile (chỉ dùng mô hình 2D RMA2 lưới chữ nhật thô $20\times 20\text{ m}$, làm mất chi tiết hình học cống) và Shahidan và cs. (2012) trên sông Ijok Malaysia (chỉ ghép nối HEC-RAS 1D với CCHE2D), luận án đã thiết lập mô hình 3D phi cấu trúc MIKE 3 FM hoàn chỉnh với $24.850$ phần tử tam giác mịn ($2\div 5\text{ m}$ tại công trình), phân chia 5 lớp Sigma theo chiều sâu. Điều này cho phép mô phỏng chính xác trường vận tốc xoắn thứ sinh và sự phân tách pha nước - cát qua từng cao trình ngưỡng cống.

3. Phát hiện thực nghiệm nào gây bất ngờ nhất và bằng chứng dữ liệu đi kèm?

Phát hiện bất ngờ nhất là cống Cẩm Đình cũ dù được thiết kế với ngưỡng đáy rất thấp ($+3{,}0\text{ m}$ mùa kiệt) để đón nước nhưng lại gần như không lấy được nước ($Q_{\text{kiệt}}$ chỉ đạt $5\div 12\text{ m}^3/\text{s}$, kém 3 lần thiết kế), nguyên nhân nghịch lý không phải do cống quá hẹp mà do góc đặt cống tạo ra vùng xoáy tĩnh làm bồi lắng phù sa đáy dày tới $1{,}03\text{ m}$. Trong khi đó, chỉ cần dịch chuyển tuyến cống $600\text{ m}$ về phía hạ lưu (vị trí 2), đáy sông tự động duy trì cơ chế tự xói sạch ($v = 0{,}65\div 0{,}85\text{ m/s}$) và lưu lượng tháo tăng vọt lên $35\div 100\text{ m}^3/\text{s}$.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) rõ ràng không?

Có. Luận án cung cấp chi tiết toàn bộ tham số thiết lập mô hình MIKE 3 FM: tọa độ lưới tính toán, ma trận phân giải lớp Sigma, hệ số nhám đáy Manning ($n = 0{,}022\div 0{,}028$), hệ số khuếch tán rối ngang/đứng, chuỗi số liệu biên thủy văn 1996, 2001, 2002, 2003–2004 và các bước giải phương trình động lượng bằng phương pháp thể tích hữu hạn (Finite Volume Method), cho phép bất kỳ chuyên gia nào cũng có thể tái lập và kiểm chứng kết quả mô phỏng.

5. Chương trình nghiên cứu 10 năm tới được phác thảo như thế nào?

Chương trình nghiên cứu dài hạn tập trung vào: (1) Phát triển mô hình tích hợp AI - Physics-Informed Neural Networks để dự báo xói bồi thời gian thực; (2) Nghiên cứu cấu trúc công trình giảm cát chủ động bằng hệ thống cánh lái dòng thông minh tự động (Smart Vanes); (3) Đánh giá tương tác sinh thái - địa mạo - kinh tế xã hội toàn diện khi sông Đáy được hồi sinh dòng chảy liên tục.

Kết luận

  1. Khẳng định nguyên nhân gốc rễ: Luận án đã làm sáng tỏ bằng chứng khoa học thực nghiệm về sự thất bại của cống Cẩm Đình hiện trạng là do vị trí đặt cống sai lệch so với tâm động lực dòng chảy mùa kiệt, gây bồi lấp phù sa trung bình $1{,}03\text{ m}$ tại kênh dẫn.
  2. Xác lập Bộ tiêu chí kỹ thuật chuẩn mực: Đề xuất thành công Bộ tiêu chí 4 nhóm (Thủy lực dòng dẫn, Hình thái sông cong, Góc lệch tuyến công trình $\alpha \approx 30^\circ$, và Giảm thiểu bùn cát xâm nhập) làm cơ sở định vị công trình lấy nước bồi tích.
  3. Đề xuất vị trí xây dựng công trình tối ưu: Khẳng định vị trí Cẩm Đình mới (cách cống cũ $600\text{ m}$ về phía hạ lưu) là phương án hoàn hảo nhất, đảm bảo cấp nước thường xuyên mùa kiệt $30\div 100\text{ m}^3/\text{s}$ và vận hành xả lũ an toàn $2.500\text{ m}^3/\text{s}$ theo Nghị định 04/2011/NĐ-CP.
  4. Tiên phong phương pháp luận 3D: Thiết lập chuẩn mực ứng dụng mô hình số trị 3 chiều phi cấu trúc MIKE 3 FM tích hợp dòng phân tầng dị trọng trong nghiên cứu chỉnh trị sông tại Việt Nam.
  5. Mở ra các nhánh nghiên cứu mới: Đặt nền móng cho các hướng nghiên cứu công trình ngầm nắn dòng tự động, dự báo bồi lắng bằng trí tuệ nhân tạo và quy hoạch phục hồi sinh thái các dòng sông chết trên quy mô toàn cầu.
  6. Di sản thực tiễn lâu dài: Cung cấp cơ sở khoa học tối hậu để Thủ tướng Chính phủ, Bộ Nông nghiệp & PTNT và TP. Hà Nội quyết định tái thiết cụm đầu mối Hát Môn - Cẩm Đình - Hiệp Thuận, mở ra kỷ nguyên hồi sinh bền vững cho dòng sông Đáy lịch sử.