Tổng quan về luận án

Nghiên cứu của Tiến sĩ Kỹ thuật Trần Khắc Thạc (Trường Đại học Thủy lợi, năm 2019) với đề tài "Nghiên cứu cơ sở khoa học và thực tiễn đề xuất vị trí hợp lý của cửa lấy nước vào sông Đáy đảm bảo yêu cầu cấp nước và tạo dòng chảy thường xuyên cho sông Đáy" (Chuyên ngành Thủy văn học, Mã số: 9440224) giải quyết bài toán cấp bách về an ninh nguồn nước và động lực học dòng chảy hạ du sông Hồng. Bối cảnh nghiên cứu xuất phát từ tác động kép của biến đổi khí hậu toàn cầu và hệ quả điều tiết của chuỗi hồ chứa thủy điện thượng nguồn (điển hình là các hồ chứa phía Trung Quốc tích trữ từ 10% đến hơn 30% lượng nước giai đoạn 2007–2009), khiến chế độ thủy văn hạ du biến đổi sâu sắc. Minh chứng cụ thể từ số liệu thực đo: "mực nước sông Hồng vào mùa kiệt thường xuyên hạ thấp, nhất là năm 2010 mực nước tại Hà Nội còn +0,1m", đe dọa nghiêm trọng khả năng vận hành của các công trình lấy nước trọng điểm.

Khoảng trống nghiên cứu (research gap) cốt lõi được định vị: Sự suy thoái nghiêm trọng của sông Đáy (chiều dài 240 km, diện tích lưu vực $5.910 \text{ km}^2$) kể từ sau khi xây dựng Đập Đáy năm 1937, làm đoạn sông dài 23 km từ Đập Đáy đến Cầu Mai Lĩnh gần như biến thành "dòng sông chết". Dù cụm công trình cống Cẩm Đình – Hiệp Thuận được xây dựng giai đoạn 2002–2004 với thiết kế lưu lượng $Q_{\text{kiệt}} = 36,24 \text{ m}^3/\text{s}$ và $Q_{\text{lũ}} = 70 \text{ m}^3/\text{s}$, thực tế vận hành hoàn toàn tê liệt vào mùa kiệt do hiện tượng bồi lắng lòng dẫn nghiêm trọng. Luận án ghi nhận: "Kênh dẫn thượng lưu cống Cẩm Đình dài 700m từ cửa cống ra sông Hồng có bề rộng đáy kênh là 22m. Theo điều tra của nghiên cứu sinh năm 2014 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."

Để hiện thực hóa chủ trương chiến lược: "Theo nghị định 04/2011/NĐ-CP của Chính phủ, với mục tiêu làm sống lại sông Đáy quy định yêu cầu cấp nước thường xuyên: Mùa kiệt từ 30 - 100m3/s; Mùa lũ là 500 - 800m3/s; Và đảm nhiệm phân lũ với lưu lượng 2.500m3/s", luận án đặt ra hai câu hỏi nghiên cứu then chốt:

  1. RQ1: Cơ chế động lực học phân tầng và tương tác dòng chảy - bùn cát tại khu vực cửa lấy nước bên sông biến đổi ra sao dưới điều kiện mực nước kiệt lịch sử và vận hành hồ chứa thượng nguồn?
  2. RQ2: Vị trí, góc lấy nước và quy mô hình học tối ưu nào của cửa lấy nước cho phép tối đa hóa lưu lượng dẫn vào mùa kiệt ($30 - 100 \text{ m}^3/\text{s}$), đáp ứng thoát lũ thường xuyên ($500 - 800 \text{ m}^3/\text{s}$) và cắt lũ khẩn cấp ($2.500 \text{ m}^3/\text{s}$), đồng thời kiểm soát bồi lắng bền vững?

Hai giả thuyết nghiên cứu tương ứng:

  • H1: Hiện tượng bồi lấp cống Cẩm Đình hiện trạng là do vị trí đặt cửa lấy nước chưa tối ưu theo quy luật dòng dị trọng và độ cong luồng lạch, dẫn đến vận tốc dòng chảy trước cửa cống giảm sâu dưới ngưỡng vận chuyển bùn cát lơ lửng ($d_{50} = 0,029 \text{ mm}$).
  • H2: Việc dịch chuyển tuyến công trình về phía hạ lưu kết hợp điều chỉnh góc lấy nước hợp lý sẽ tạo ra trường thủy lực thuận lợi, giảm thiểu nồng độ bùn cát đáy xâm nhập và duy trì ổn định dòng chảy tự chảy mùa kiệt.

Phạm vi không gian nghiên cứu bao quát đoạn sông Hồng dài 27 km từ Km30 đến Km47+500 (khu vực cửa vào sông Đáy), ứng dụng khung lý thuyết dòng chảy phân tầng/dị trọng kết hợp mô hình số trị thủy động lực học 3D (MIKE 3 FM) liên kết mạng thủy lực 1D (MIKE 11). Nghiên cứu tạo bước đột phá trong việc thiết lập cơ sở khoa học để tái định vị công trình lấy nước đầu mối, phục vụ trực tiếp cho quy hoạch Thủ đô Hà Nội tầm nhìn 2030–2050.

Literature Review và Positioning

Tổng quan y văn quốc tế và trong nước về công trình lấy nước (CTLN) và động lực học bồi lắng lòng dẫn bộc lộ rõ các luồng quan điểm và hạn chế phương pháp luận:

                      TIẾN TRÌNH NGHIÊN CỨU THỦY LỰC CỬA LẤY NƯỚC
┌─────────────────────────────────────────────────────────────────────────────────┐
│ Mô hình 1D/2D Thô (HEC-RAS, RMA2)                                              │
│ • Nguyễn Thọ Sáo (2005), Elmoustafa (2011)                                     │
│ • Nhược điểm: Bỏ qua biến dạng trường vận tốc theo chiều sâu (Stratified Flow)  │
└──────────────────────────────────────┬──────────────────────────────────────────┘
                                       │
                                       ▼
┌─────────────────────────────────────────────────────────────────────────────────┐
│ Mô hình Vật lý & Bán Thực nghiệm (Froude 0.35 - 0.40)                           │
│ • Karami Moghadam (2010), Masjedi (2011), N.F. Danhelia (Bán kính cong)        │
│ • Nhược điểm: Hiệu ứng tỷ lệ bùn cát, góc lấy nước dao động rộng (15° - 90°)    │
└──────────────────────────────────────┬──────────────────────────────────────────┘
                                       │
                                       ▼
┌─────────────────────────────────────────────────────────────────────────────────┐
│ Luận án: Khung Mô hình Hóa 3D Đa quy mô (MIKE 11 HD/ST + MIKE 3 FM Lưới Tam giác)│
│ • Tích hợp Cơ học Dòng Dị trọng/Phân tầng + Mô hình Rối k-ε/RNG                 │
│ • Giải quyết đồng thời: Lấy nước mùa kiệt (30-100 m³/s) & Thoát lũ (2.500 m³/s) │
└─────────────────────────────────────────────────────────────────────────────────┘

Trong dòng nghiên cứu quốc tế:

  • Phân loại và hình học công trình: Avery (1989) và GTZ (1989) phân loại CTLN bên cạnh (lateral intake, góc lấy $\approx 90^\circ$), CTLN chính diện và CTLN đáy dạng Tyrolean. Abbasi (2003) nhấn mạnh vai trò của đập dâng kiểm soát mực nước tối thiểu $Q_r$.
  • Góc lấy nước và bồi lắng: Tồn tại tranh luận học thuật sâu sắc về góc lấy nước tối ưu. Mehdi Karami Moghadam (2010) nghiên cứu thực nghiệm trên sông Ohio kết luận góc $30^\circ$ cho lưu lượng lớn nhất và bùn cát nhỏ nhất, nhưng chỉ giới hạn trong dải số Froude $Fr = 0,35 - 0,40$ và hoàn toàn bỏ qua dòng chảy kiệt. Alireza Masjedi (2011) đề xuất góc $45^\circ$, trong khi Afrouzi (2015) khuyến nghị $30^\circ - 45^\circ$, và Montaseri (2015) dùng mô hình số 3D (phần mềm Fluent với mô hình Reynolds Stress) khẳng định góc tối ưu là $17^\circ - 20^\circ$ để giảm phù sa.
  • Mô hình hóa số trị: Neary et al. (1999) tiên phong giải phương trình Navier-Stokes 3D với mô hình rối Wilcox cho kênh chữ T. Ngược lại, Ashraf Elmoustafa (2011) khi đánh giá CTLN bờ Tây sông Nile (Assiut, Ai Cập) lại dùng mô hình 2D RMA2 với lưới chữ nhật thô $20 \times 20 \text{ m}$, không phản ánh được cấu trúc dòng chảy sát đáy và hình học công trình chi tiết. Noor et al. (2012) kết hợp HEC-RAS và CCHE2D trên sông Ijok (Malaysia), xác nhận vận tốc nhỏ trước cửa lấy nước gây bồi tích cục bộ làm suy giảm dung lượng lấy nước. Về mặt chỉnh trị hình thái, Margriet (2006) phân tích các đợt lũ lịch sử 1993, 1995 trên sông Rhine ($3.200 \text{ km}$, lưu vực $185.000 \text{ km}^2$), chứng minh việc lấn chiếm bãi sông làm tê liệt khả năng cấp nước mùa kiệt và phá hủy lòng dẫn ổn định.

Tại Việt Nam, các nghiên cứu dọc tuyến sông Hồng đã phát triển qua nhiều giai đoạn:

  • Viện Khoa học Thủy lợi Việt Nam (2001, 2003, 2005, 2008–2010) và Hà Văn Khôi et al. (2009) tập trung chủ yếu vào các phương án công trình bảo vệ bờ (Hải Bối, Phú Gia - Tứ Liên, Bát Tràng), xóa bỏ vùng chậm lũ và cải tạo thoát lũ sông Đáy ở quy mô 1D.
  • Nguyễn Thọ Sáo (2005) dùng mô hình 1D phân tích thủy động lực học nhưng chưa thể hiện biến dạng lòng dẫn.
  • Phạm Đức Thắng (2010) ứng dụng mô hình số 3D nghiên cứu bồi lắng cống Xuân Quan nhưng chỉ mô phỏng trong điều kiện mùa lũ khi đóng cống không lấy nước, bỏ qua tương tác động lực dòng chảy lấy nước mùa kiệt.
  • Phạm Thị Hương Lan (2012) sử dụng MIKE 11 ST và MIKE 21 FM mô phỏng hạ lưu cửa vào sông Đáy theo Nghị định 04/2011/NĐ-CP, song chưa định lượng được cơ chế bồi lấp cống Cẩm Đình cũng như tiêu chí xác định vị trí tuyến công trình tối ưu.

Định vị của luận án: Luận án của Trần Khắc Thạc đã lấp đầy khoảng trống phương pháp luận bằng cách tích hợp lý thuyết dòng phân tầng/dị trọng với mô hình thủy động lực học 3D lưới phi cấu trúc (MIKE 3 FM), giải quyết bài toán đa mục tiêu: vừa cấp nước kiệt tự chảy vừa đảm bảo năng lực xả/thoát lũ khẩn cấp $2.500 \text{ m}^3/\text{s}$, đưa ra bộ tiêu chí định lượng chuẩn xác thay thế các công thức kinh nghiệm bán kính cong đơn thuần của N.F. Danhelia vốn không còn thích ứng với biến đổi hình thái sông hiện đại.

Đó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 các lý thuyết nền tảng trong thủy lực học công trình và động lực học sông ngòi:

                            KHUNG PHÂN TÍCH LÝ THUYẾT DÒNG DỊ TRỌNG
┌────────────────────────────────────────────────────────────────────────────────────────┐
│ PHƯƠNG TRÌNH NAVIER-STOKES 3D (DÒNG PHÂN TẦNG LIÊN TỤC)                                │
│                                                                                        │
│   ∂u/∂t + u(∂u/∂x) + v(∂u/∂y) + w(∂u/∂z) = - (1/ρ)(∂p/∂x) + ν ∇²u + f_cor v           │
│   ∂v/∂t + u(∂v/∂x) + v(∂v/∂y) + w(∂v/∂z) = - (1/ρ)(∂p/∂y) + ν ∇²v - f_cor u           │
│   ∂w/∂t + u(∂w/∂x) + v(∂w/∂y) + w(∂w/∂z) = - (1/ρ)(∂p/∂z) + ν ∇²w - g [1 + β_s (S - S₀)]│
│                                                                      ▲                 │
│                                              Thành phần Gradient Khối lượng riêng      │
└──────────────────────────────────────────────────┬─────────────────────────────────────┘
                                                   │ Tích hợp vận chuyển phù sa
                                                   ▼
┌────────────────────────────────────────────────────────────────────────────────────────┐
│ KHUẾCH TÁN BÙN CÁT 3 CHIỀU & CÂN BẰNG NĂNG LƯỢNG RỐI                                   │
│   ∂S/∂t + u(∂S/∂x) + v(∂S/∂y) + (w - ω_s)(∂S/∂z) = ∂/∂x(ε_x ∂S/∂x) + ∂/∂y(ε_y ∂S/∂y)  │
│                                                   + ∂/∂z(ε_z ∂S/∂z)                    │
│                                                                                        │
│ • Đóng góp: Xác lập cơ chế phân tách tầng nước ngọt mặt và tầng bùn cát đáy tại góc rẽ│
│ • Phá vỡ giả định trường vận tốc đồng nhất (Uniform Velocity Profile) của mô hình 2D/1D│
└────────────────────────────────────────────────────────────────────────────────────────┘
  1. Lý thuyết Dòng dị trọng/Dòng phân tầng (Density-Stratified Flow): Ứng dụng hệ phương trình liên tục và phương trình động lượng Navier-Stokes 3 chiều có xét tới độ dốc mật độ chất lưu do nồng độ bùn cát lơ lửng $S(x,y,z,t)$: $$\rho = \rho_0 [1 + \beta_s (S - S_0)]$$ Luận án chứng minh rằng tại khu vực cửa lấy nước ven sông, gradient mật độ thẳng đứng kết hợp với lực ly tâm của dòng chảy vòng tại khúc sông cong tạo ra hoàn lưu thứ cấp (secondary circulation). Hoàn lưu này đẩy bùn cát đáy tích tụ mạnh về phía cửa lấy nước nếu góc rẽ không tương thích với trường động năng dòng chính.
  2. Mô hình Vận chuyển Bùn cát Đa chiều: Bổ sung cơ chế tương tác giữa vận tốc lắng chìm của hạt $\omega_s$ (đối với đường kính $d_{50} = 0,029 \text{ mm}$) và ten-sơ ứng suất rối Reynolds trong điều kiện dòng chảy kiệt biến đổi chậm, làm sáng tỏ hiện tượng suy giảm sức tải cát khi dòng nước chuyển hướng đột ngột vào kênh dẫn $B = 22 \text{ m}$.
  3. Mô hình Tiêu chí Đa thông số Hình học - Thủy lực: Chuyển dịch từ việc áp dụng bán kính cong hình học tĩnh $R$ (theo Danhelia) sang hàm tương quan động lực đa biến: $$Q_{\text{intake}} = f(H_{\text{sông}}, \theta_{\text{góc}}, L_{\text{đoạn cong}}, L_{\text{kênh dẫn}}, S_{\text{bùn cát}})$$

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

Khung phân tích của luận án tích hợp liên hoàn 3 trụ cột lý thuyết: Cơ học chất lỏng thủy lực sông (River Hydraulics), Động lực học bùn cát (Sediment Transport Dynamics) và Lý thuyết Điều khiển Tối ưu Hệ thống Thủy công.

┌─────────────────────────────────────────────────────────────────────────────────────────┐
│                    BỘ TIÊU CHÍ KỸ THUẬT LỰA CHỌN VỊ TRÍ CỬA LẤY NƯỚC                     │
└────────────────────────────────────────────┬────────────────────────────────────────────┘
                                             │
      ┌──────────────────────────────────────┼─────────────────────────────────────┐
      │                                      │                                     │
      ▼                                      ▼                                     ▼
┌───────────────────────────┐ ┌───────────────────────────┐ ┌───────────────────────────┐
│   TIÊU CHÍ ĐỘNG LỰC HỌC   │ │    TIÊU CHÍ HÌNH THÁI     │ │    TIÊU CHÍ VẬN HÀNH      │
│         DÒNG CHẢY         │ │       LÒNG DẪN            │ │        CÔNG TRÌNH         │
├───────────────────────────┤ ├───────────────────────────┤ ├───────────────────────────┤
│• Vận tốc tiếp cận v ≥ v_k │ │• Nằm ở bờ lõm đoạn sông   │ │• Q_kiệt = 30 - 100 m³/s   │
│  (vận tốc không bồi lắng) │ │  cong ổn định động lực    │ │  (tự chảy ứng H_kiệt)     │
│• Góc lấy nước θ tối ưu:   │ │• Chiều dài đoạn cong L_c  │ │• Năng lực xả lũ kiểm soát:│
│  30° ≤ θ ≤ 45°            │ │  đủ duy trì tim chủ lưu   │ │  Q_lũ = 500 - 800 m³/s    │
│• Kiểm soát hoàn lưu đáy   │ │• Chiều dài kênh dẫn L_k   │ │• Khả năng phân lũ sự cố:  │
│  ngăn bùn cát thô xâm nhập│ │  ngắn nhất (giảm tổn thất)│ │  Q_max = 2.500 m³/s       │
└───────────────────────────┘ └───────────────────────────┘ └───────────────────────────┘
  • Điều kiện biên giới hạn (Boundary Conditions): Khung phân tích áp dụng cho các hệ thống phân lưu vùng đồng bằng châu thổ chịu ảnh hưởng điều tiết hồ chứa thượng nguồn, biên độ dao động mực nước giữa lũ và kiệt lớn ($H_{\text{kiệt}} \approx +0,1 \text{ m}$ đến $H_{\text{lũ}} > +13 \text{ m}$ tại trạm Hà Nội), và lòng dẫn cấu tạo bởi trầm tích cát mịn hạt rời.

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

Thiết kế nghiên cứu

Luận án tuân thủ nghiêm ngặt quan điểm thực chứng (positivism) kết hợp mô hình hóa số trị định lượng đa cấp (multi-level computational design):

  • Cấp vĩ mô (Hệ thống lưu vực 1D): Toàn bộ mạng lưới sông Hồng – sông Thái Bình nhằm thiết lập chính xác các điều kiện biên thủy lực thượng lưu (Sơn Tây, Hòa Bình, Thác Bà) và biên triều hạ lưu.
  • Cấp vi mô (Đoạn sông chi tiết 3D): Đoạn sông Hồng dài 27 km (Km30 đến Km47+500) kết hợp toàn tuyến kênh dẫn Cẩm Đình – Hiệp Thuận ($11.295 \text{ m}$) và khu vực cửa vào sông Đáy.

Quy trình nghiên cứu rigorous

┌────────────────────────────────────────────────────────────────────────────────┐
│                       QUY TRÌNH MÔ HÌNH HÓA VÀ KIỂM ĐỊNH                       │
└───────────────────────────────────────┬────────────────────────────────────────┘
                                        │
      ┌─────────────────────────────────┴─────────────────────────────────┐
      │                                                                   │
      ▼                                                                   ▼
┌───────────────────────────────────────────────┐ ┌───────────────────────────────────────────────┐
│     MÔ HÌNH 1D TOÀN HỆ THỐNG (MIKE 11)        │ │        MÔ HÌNH 3D CHI TIẾT (MIKE 3 FM)        │
├───────────────────────────────────────────────┤ ├───────────────────────────────────────────────┤
│• Mạng sông: Hệ thống Hồng - Thái Bình         │ │• Phạm vi: 27 km sông Hồng (Km30 - Km47+500)   │
│• Biên thượng lưu: Thủy văn Sơn Tây, Hòa Bình  │ │• Lưới phi cấu trúc (Flexible Mesh):           │
│• Biên hạ lưu: Dao động triều ven biển         │ │  Tam giác mịn hóa tại công trình cống         │
│• Nhiệm vụ: Cung cấp chuỗi Q(t), H(t) biên cho │ │• Phân tầng đứng: Sigma-coordinates (5-10 lớp) │
│  mô hình MIKE 3 FM mùa lũ và mùa kiệt         │ │• Giải phương trình Navier-Stokes 3D & Mud/ST  │
└───────────────────────┬───────────────────────┘ └───────────────────────┬───────────────────────┘
                        │                                                 │
                        └────────────────────────┬────────────────────────┘
                                                 │
                                                 ▼
┌────────────────────────────────────────────────────────────────────────────────────────┐
│                        HIỆU CHỈNH VÀ KIỂM ĐỊNH MÔ HÌNH ĐỘC LẬP                         │
├────────────────────────────────────────────────────────────────────────────────────────┤
│• Chuỗi số liệu Lũ: Trận lũ lớn năm 1996 (Hiệu chỉnh) và Trận lũ năm 2002 (Kiểm định)   │
│• Chuỗi số liệu Kiệt: Mùa kiệt 2001 (Hiệu chỉnh) và Mùa kiệt 2003, 2003-2004 (Kiểm định)│
│• Vị trí kiểm tra: Trạm Thượng Cát, Trạm Hà Nội và Trạm Mặt cắt chi tiết MC-1            │
│• Tiêu chuẩn khớp: Đường tương quan Q_nước và Q_bùn cát đạt R² > 0.90                   │
└────────────────────────────────────────────────────────────────────────────────────────┘

Quy trình bảo đảm tính hợp lệ cấu trúc (construct validity) qua việc đồng bộ hóa dữ liệu địa hình lòng sông từ các bình đồ đo đạc mới nhất kết hợp ảnh viễn thám và khảo sát mặt cắt ngang thực tế năm 2014 của Nghiên cứu sinh.

Data và phân tích

  • Môi trường tính toán: Bộ phần mềm thủy lực hàng đầu thế giới của Viện Thủy động lực Đan Mạch (DHI): MIKE 11 (HD/ST) liên kết MIKE 3 FM (Hydrodynamic & Sand Transport Modules).
  • Hệ thống lưới tính: Lưới tam giác phi cấu trúc (Flexible Mesh) linh hoạt, tăng dày mật độ tại khu vực cửa cống lấy nước và kênh dẫn để bắt trọn các dòng xoáy cục bộ và hiện tượng tách dòng. Phân chia theo phương đứng sử dụng hệ tọa độ $\sigma$ (sigma-coordinates) đảm bảo mô phỏng chính xác ứng suất đáy và phân bố nồng độ bùn cát theo phương thẳng đứng.
  • Dữ liệu đầu vào và kiểm định:
    • Mùa kiệt: Số liệu quan trắc liên tục các mùa kiệt 2001 (hiệu chỉnh), 2003 và chuỗi 5 tháng mùa kiệt 2003–2004 (từ 01/12 đến 30/04 năm sau).
    • Mùa lũ: Trận lũ lịch sử năm 1996 với lưu lượng đỉnh lớn (hiệu chỉnh) và trận lũ năm 2002 (kiểm định).
    • Đặc trưng trầm tích: Phù sa lơ lửng lòng sông Hồng đoạn qua Sơn Tây có đường kính hạt trung bình $d_{50} = 0,029 \text{ mm}$, hệ số thấm tầng chứa nước cát cuội sỏi thềm sông $K = 5 \times 10^{-3} \text{ cm/s}$.

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

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

┌───────────────────────────────────────────────────────────────────────────────────────────┐
│                 SO SÁNH CÁC PHƯƠNG ÁN VỊ TRÍ CỬA LẤY NƯỚC VÀO SÔNG ĐÁY                    │
├──────────────────────────┬─────────────────────────────┬──────────────────────────────────┤
│ THÔNG SỐ ĐÁNH GIÁ        │ PHƯƠNG ÁN HIỆN TRẠNG        │ PHƯƠNG ÁN ĐỀ XUẤT MỚI            │
│                          │ (Cống Cẩm Đình Cũ - Vị trí 1)│ (Hạ lưu 600m - Vị trí 2)         │
├──────────────────────────┼─────────────────────────────┼──────────────────────────────────┤
│ Vị trí tuyến             │ Tuyến cống hiện có tại K0   │ Dịch chuyển 600m về phía hạ lưu  │
│ Góc lấy nước (θ)         │ Góc lấy nước tù/bất lợi     │ Tối ưu hóa 30° - 45°             │
│ Vận tốc dòng dẫn kiệt    │ v < 0,2 m/s (Bồi lắng mạnh) │ v = 0,45 - 0,75 m/s              │
│ Cao độ bồi lắng trung bình│ +1,03 m (Lấp tắc kênh dẫn)  │ Duy trì ổn định xói rửa tự nhiên │
│ Khả năng lấy nước kiệt   │ Gần như bằng 0 (Tê liệt)    │ Đảm bảo 30 - 100 m³/s            │
│ Năng lực chuyển lũ       │ Không đạt Q_tk = 70 m³/s    │ Đáp ứng 500 - 800 m³/s (cắt lũ)  │
└──────────────────────────┴─────────────────────────────┴──────────────────────────────────┘
  1. Giải mã nguyên nhân bồi lấp cống Cẩm Đình hiện trạng (Vị trí 1): Mô phỏng 3D chỉ rõ vận tốc dòng chảy trước cửa lấy nước Cẩm Đình trong mùa kiệt giảm xuống rất thấp ($v < 0,2 \text{ m/s}$), thấp hơn nhiều so với ngưỡng vận tốc không bồi lắng đối với hạt phù sa $d_{50} = 0,029 \text{ mm}$. Kết quả là bùn cát lơ lửng chuyển thành bồi lắng nhanh chóng, làm lòng dẫn bị nâng cao trung bình $1,03 \text{ m}$ (có vị trí bồi cục bộ tới $2 \text{ m}$ sau 5 năm), triệt tiêu hoàn toàn khả năng lấy nước tự chảy thiết kế $Q = 36,24 \text{ m}^3/\text{s}$.
  2. Khám phá vị trí lấy nước tối ưu (Vị trí 2 - Cách cống cũ 600m về hạ lưu): Kết quả chạy các kịch bản thủy lực 3D khẳng định việc dịch chuyển cửa lấy nước về phía hạ lưu $600 \text{ m}$ (nằm tại vị trí bờ lõm của cung cong có bán kính động lực học ổn định) giúp đưa dòng chủ lưu áp sát cửa lấy nước. Tại đây, trường vận tốc đạt giá trị thuận lợi ($0,45 - 0,75 \text{ m/s}$ vào mùa kiệt), hạn chế tối đa góc xoáy phân lưu chết và giảm trên $65%$ lượng bùn cát hạt thô xâm nhập vào kênh dẫn.
  3. Mối quan hệ định lượng giữa Hình học - Lưu lượng - Mực nước: Xác lập biểu đồ tương quan 3 biến giữa Lưu lượng nước lấy vào ($Q_{\text{intake}}$) – Mực nước sông Hồng ($H_{\text{sông}}$) – Góc lấy nước ($\theta$). Góc lấy nước $\theta = 30^\circ - 45^\circ$ kết hợp chiều dài kênh dẫn thu ngắn mang lại hiệu suất thủy lực cao nhất, bảo đảm đưa đủ $30 - 100 \text{ m}^3/\text{s}$ nước kiệt vào sông Đáy ngay cả khi mực nước sông Hồng hạ thấp cực đoan.
  4. Hiệu ứng phân phối lũ và hồi sinh sinh thái: Khi thực hiện phương án cải tạo lòng dẫn sông Đáy kết hợp công trình phân lũ mới tại vị trí đề xuất, khả năng thoát lũ định kỳ $500 - 800 \text{ m}^3/\text{s}$ và phân lũ khẩn cấp $2.500 \text{ m}^3/\text{s}$ được đảm bảo an toàn tuyệt đối mà không gây ngập úng bãi sông vùng Chương Mỹ, Mỹ Đức, đồng thời khống chế mực nước lũ tại trạm Hà Nội dưới ngưỡng an toàn $+13,40 \text{ m}$.

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp bằng chứng thực nghiệm và số trị khẳng định giới hạn của lý thuyết Danhelia; xác lập phương pháp luận tính toán CTLN dựa trên động lực dòng dị trọng 3D thay vì thủy lực 1D/2D truyền thống.
  • Về mặt phương pháp: Chuẩn hóa quy trình liên kết động mô hình 1D quy mô lưu vực (MIKE 11) với mô hình 3D lòng động chi tiết (MIKE 3 FM) trong điều kiện dữ liệu địa hình sông ngòi nhiệt đới biến động mạnh.
  • Về mặt thực tiễn và chính sách: Đưa ra giải pháp công trình và phi công trình hoàn chỉnh: Đề xuất dịch chuyển tuyến công trình lấy nước Cẩm Đình mới; cải tạo kênh dẫn Cẩm Đình – Hiệp Thuận với mặt cắt đáy rộng $B = 150 \text{ m}$; xây dựng quy trình vận hành liên hồ chứa thượng nguồn kết hợp xả lũ/lấy nước chủ động theo mùa, làm cơ sở khoa học trực tiếp cho Bộ Nông nghiệp và PTNT, UBND TP. Hà Nội điều chỉnh Quy hoạch phòng chống lũ và quy hoạch thủy lợi sông Đáy.

Limitations và Future Research

Luận án thẳng thắn chỉ rõ các giới hạn học thuật nội tại:

  1. Dữ liệu bùn cát đáy: Do điều kiện đo đạc mùa lũ nguy hiểm, số liệu quan trắc bùn cát đáy (bed load) còn hạn chế; mô hình chủ yếu hiệu chỉnh theo bùn cát lơ lửng (suspended load).
  2. Biến động hình thái dài hạn: Mô hình 3D chưa tích hợp đầy đủ quá trình xói mòn sinh học và biến đổi dòng chảy dưới tác động của toàn bộ chuỗi hồ chứa bậc thang sông Đà (Sơn La, Lai Châu) trong chu kỳ 30–50 năm.
  3. Tính chất cơ lý lòng dẫn: Giả định nền đáy sông có độ nhám đồng nhất cục bộ theo vùng, chưa phản ánh hết sự phân lớp phức tạp của tầng cát pha cuội sỏi tại các hố xói sâu.

Định hướng nghiên cứu tiếp nối:

  • Mở rộng mô hình hình thái sông 3D liên tục (morphodynamic 3D) mô phỏng dài hạn từ 10 đến 20 năm cho toàn dải phân lưu sông Hồng – sông Đáy.
  • Nghiên cứu giải pháp tấm lái dòng ngập (submerged vanes) tự động điều chỉnh theo mực nước để nắn dòng chủ lưu đẩy cát xa cửa cống.
  • Tích hợp mô hình chất lượng nước (MIKE 11 WQ/ECO Lab) đánh giá khả năng tự làm sạch của sông Đáy khi được bổ cập lưu lượng thường xuyên.

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

                                  MA TRẬN TÁC ĐỘNG VÀ ẢNH HƯỞNG
┌────────────────────────┬───────────────────────────────────────────────────────────────────────┐
│ TRỤ CỘT TÁC ĐỘNG       │ MINH CHỨNG VÀ CHỈ SỐ LƯỢNG HÓA                                        │
├────────────────────────┼───────────────────────────────────────────────────────────────────────┤
│ 1. Khoa học & Học thuật│ • Khung phương pháp luận chuẩn xác cho thiết kế CTLN vùng đồng bằng   │
│                        │ • Dự báo tạo nguồn trích dẫn cao trong ngành Thủy văn & Kỹ thuật Sông  │
├────────────────────────┼───────────────────────────────────────────────────────────────────────┤
│ 2. Nông nghiệp & Dân sinh│ • Đảm bảo cấp nước tưới tự chảy trực tiếp cho >10.000 ha đất canh tác │
│                        │ • Cắt giảm 100% chi phí điện năng bơm dã chiến và nạo vét sự cố      │
├────────────────────────┼───────────────────────────────────────────────────────────────────────┤
│ 3. Môi trường Sinh thái│ • Hồi sinh 23 km "dòng sông chết" từ Đập Đáy đến Mai Lĩnh             │
│                        │ • Khôi phục hệ sinh thái thủy sinh và cảnh quan đô thị sinh thái      │
├────────────────────────┼───────────────────────────────────────────────────────────────────────┤
│ 4. Thể chế & Quy hoạch │ • Cung cấp luận cứ cốt lõi thực thi Nghị định 04/2011/NĐ-CP           │
│                        │ • Tái định hình quy hoạch phân lũ sông Đáy trong bối cảnh hồ Sơn La   │
└────────────────────────┴───────────────────────────────────────────────────────────────────────┘
  • Tác động học thuật: Đặt nền móng chuyển đổi phương pháp luận thiết kế công trình cửa lấy nước bên sông tại Việt Nam từ kinh nghiệm thực nghiệm sang mô phỏng số 3D dòng phân tầng chính xác cao.
  • Tác động kinh tế – xã hội: Khôi phục lưu thông dòng chảy giúp tái sinh hệ thống thủy nông Đan Hoài, Phù Sa, sông Tích; bảo đảm nguồn nước sinh hoạt, công nghiệp và du lịch cảnh quan cho vùng Tây Nam Hà Nội và các tỉnh hạ du (Hà Nam, Ninh Bình, Nam Định).
  • Tác động phòng chống thiên tai: Tối ưu hóa dung tích phân lũ $2.500 \text{ m}^3/\text{s}$, bảo vệ an toàn tuyệt đối cho trung tâm Thủ đô Hà Nội trước các đợt lũ chu kỳ 500 năm.

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

                                     BẢN ĐỒ ĐỐI TƯỢNG HƯỞNG LỢI
┌────────────────────────┬───────────────────────────────────────────────────────────────────────┐
│ NHÓM ĐỐI TƯỢNG         │ GIÁ TRỊ VÀ LỢI ÍCH TIẾP NHẬN CỤ THỂ                                   │
├────────────────────────┼───────────────────────────────────────────────────────────────────────┤
│ Nghiên cứu sinh &      │ Khung phương pháp luận tích hợp 1D-3D MIKE FM; bộ tiêu chí định lượng │
│ Học giả ngành Nước     │ thiết kế cửa lấy nước và cơ chế dòng phân tầng / dị trọng.            │
├────────────────────────┼───────────────────────────────────────────────────────────────────────┤
│ Đơn vị Tư vấn Thiết kế │ Thông số kỹ thuật chuẩn xác về góc lấy nước (30°-45°), cao trình      │
│ Thủy lợi (HEC, VNKHTL) │ ngưỡng cống, mặt cắt kênh dẫn và vị trí tuyến công trình tối ưu.     │
├────────────────────────┼───────────────────────────────────────────────────────────────────────┤
│ Cơ quan Quản lý        │ Cơ sở khoa học vững chắc phục vụ rà soát quy hoạch thoát lũ sông Đáy, │
│ Nhà nước & Quy hoạch   │ điều hành liên hồ chứa và xóa bỏ vùng phân chậm lũ theo Nghị định 04. │
├────────────────────────┼───────────────────────────────────────────────────────────────────────┤
│ Doanh nghiệp Khai thác │ Giảm thiểu tần suất nạo vét kênh dẫn Cẩm Đình - Hiệp Thuận; vận hành │
│ Công trình Thủy lợi    │ ổn định công trình lấy nước tự chảy mùa kiệt, tiết kiệm điện nă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 xác lập cơ chế phân tầng và tách dòng động lực học (Dynamic Stratified Separation Mechanism) của dòng dị trọng chứa bùn cát tại cửa lấy nước bên sông. Luận án đã mở rộng lý thuyết dòng phân tầng kinh điển (Stratified Flow Theory) và khắc phục khiếm khuyết của quy tắc bán kính cong hình học N.F. Danhelia bằng cách chứng minh: khả năng lấy nước và bồi lắng không chỉ phụ thuộc vào bán kính cong $R$ của đoạn sông, mà bị chi phối trực tiếp bởi tương tác giữa gradient nồng độ bùn cát thẳng đứng $\partial S/\partial z$, góc lấy nước $\theta$ và trường vận tốc 3 chiều trước cửa công trình.

2. Điểm cải tiến phương pháp luận so với các nghiên cứu quốc tế và trong nước trước đây?

So với nghiên cứu của Ashraf 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}$ không thể hiện được biến dạng dòng chảy sát đáy, và nghiên cứu của Phạm Đức Thắng (2010) tại cống Xuân Quan chỉ mô phỏng 3D trong điều kiện mùa lũ đóng cống không lấy nước, luận án đã:

  1. Thiết lập mô hình 3D MIKE 3 FM với lưới phi cấu trúc tam giác mịn hóa cao độ tại công trình, giải quyết đồng thời bài toán dòng chảy 3 chiều trong cả mùa kiệt sâu lẫn mùa lũ lớn khi cống mở lấy nước.
  2. Liên kết động mô hình 1D toàn lưu vực (MIKE 11) với mô hình 3D cục bộ để truyền biên thủy lực chính xác tuyệt đối.

3. Phát hiện bất ngờ nhất (counter-intuitive finding) có số liệu chứng minh là gì?

Phát hiện bất ngờ nhất là cống Cẩm Đình hiện trạng dù được đầu tư xây dựng quy mô rất lớn với 2 tầng cống hở nhưng hoàn toàn "bất động" trong mùa kiệt không phải do thiếu nguồn nước từ sông Hồng, mà do chính sai lệch về vị trí đặt tuyến công trình. Mặc dù lưu lượng sông Hồng mùa kiệt vẫn duy trì hàng nghìn $\text{m}^3/\text{s}$, nhưng việc đặt cống tại vị trí tách dòng tạo nên vùng "nước chết" với vận tốc $v < 0,2 \text{ m/s}$, khiến $1,03 \text{ m}$ bùn cát bồi lắng lấp kín đáy kênh $B = 22 \text{ m}$, biến công trình cấp nước trọng điểm thành bẫy giữ bùn cát phù sa.

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

Có. Luận án cung cấp chi tiết toàn bộ quy trình:

  • Cơ sở dữ liệu mặt cắt địa hình sông Hồng Km30–Km47+500 và kênh dẫn.
  • Thiết lập bước lưới tam giác phi cấu trúc và thông số phân tầng $\sigma$ trong MIKE 3 FM.
  • Bộ thông số nhám Manning ($n$), hệ số khuếch tán rối ngang/đứng ($k-\epsilon$).
  • Chuỗi số liệu kiểm định mực nước, lưu lượng và hàm lượng bùn cát thực đo các năm 1996, 2001, 2002, 2003 tại các trạm Thượng Cát, Hà Nội và MC-1.

5. Chương trình nghị sự nghiên cứu 10 năm (10-Year Research Agenda) tiếp theo gồm những gì?

  1. Giai đoạn 1–3 năm: Nghiên cứu thực nghiệm lắp đặt tấm lái dòng ngập (submerged vanes) tại cửa vào cống Cẩm Đình mới để kiểm soát bồi lắng chủ động.
  2. Giai đoạn 3–6 năm: Phát triển mô hình hình thái 3D tích hợp AI dự báo bồi xói theo thời gian thực dưới tác động của các kịch bản xả nước gia tăng từ hồ thủy điện Trung Quốc và Việt Nam.
  3. Giai đoạn 6–10 năm: Thiết lập hệ thống hỗ trợ ra quyết định (DSS) vận hành liên động toàn bộ cụm đầu mối sông Đáy (Cẩm Đình – Hiệp Thuận – Đập Đáy – Vân Cốc) thích ứng với biến đổi khí hậu toàn cầu.

Kết luận

  1. Xác lập cơ sở khoa học hoàn chỉnh: Luận án đã chứng minh một cách khoa học nguyên nhân bồi lấp làm tê liệt cống Cẩm Đình hiện trạng là do vị trí đặt tuyến chưa hợp lý, gây suy giảm vận tốc tiếp cận ($v < 0,2 \text{ m/s}$) và bồi lắng trung bình $1,03 \text{ m}$ bùn cát mịn $d_{50} = 0,029 \text{ mm}$.
  2. Đề xuất vị trí công trình mới đột phá: Khẳng định vị trí xây dựng cửa lấy nước hợp lý là dịch chuyển $600 \text{ m}$ về phía hạ lưu cống Cẩm Đình cũ, bố trí ở bờ lõm đoạn cong ổn định với góc lấy nước tối ưu $\theta = 30^\circ - 45^\circ$.
  3. Giải quyết hài hòa bài toán đa mục tiêu: Đảm bảo cấp nước tự chảy thường xuyên mùa kiệt đạt $30 - 100 \text{ m}^3/\text{s}$, thoát lũ mùa mưa $500 - 800 \text{ m}^3/\text{s}$ và tham gia phân lũ khẩn cấp an toàn $2.500 \text{ m}^3/\text{s}$ theo đúng Nghị định 04/2011/NĐ-CP.
  4. Tiên phong phương pháp luận mô hình 3D: Ứng dụng thành công chuỗi mô hình MIKE 11 liên kết MIKE 3 FM lưới tam giác phi cấu trúc giải phương trình Navier-Stokes 3 chiều có xét đến dòng dị trọng trong thủy văn công trình tại Việt Nam.
  5. Đóng góp giải pháp quản lý - công trình đồng bộ: Đề xuất cải tạo mở rộng kênh dẫn Cẩm Đình – Hiệp Thuận ($B = 150 \text{ m}$) kết hợp quy trình điều tiết mực nước liên hồ chứa, tạo tiền đề vững chắc hồi sinh toàn diện môi trường sinh thái sông Đáy và phát triển bền vững Thủ đô Hà Nội.