Tổng quan về luận án

Thoái hóa đất và bồi lắng lòng hồ, sông ngòi do xói mòn là thách thức môi trường cấp bách mang tính toàn cầu. Theo báo cáo từ các tổ chức quốc tế như FAO và UNEP, "hàng năm trên thế giới có khoảng từ 5 đến 7 triệu ha đất bị mất khả năng sản xuất do bị xói mòn đất", đồng thời "tổng lượng phù sa từ các con sông đổ ra biển đã tăng từ 9 tỷ tấn lên đến 24 tỷ tấn". Tại Việt Nam, với đặc thù hơn 3/4 diện tích tự nhiên là đồi núi dốc nằm trong vùng nhiệt đới ẩm gió mùa có lượng mưa tập trung cường độ lớn (1.800 – 2.000 mm/năm), hiểm họa xói mòn đất diễn ra đặc biệt nghiêm trọng. Thống kê của Bộ Khoa học và Công nghệ chỉ ra rằng "nước ta có hơn 13 triệu ha đất bị suy thoái thành đất trống, đồi núi trọc, trong đó có 1,2 triệu ha đã bị xói mòn trơ sỏi đá". Lượng bùn cát xâm nhập vào các hệ thống sông ngòi rất lớn; điển hình như vùng đồng bằng sông Hồng, "trung bình $1\text{ m}^3$ nước trong sông chứa từ 50g - 400g bùn cát, riêng đồng bằng sông Hồng $1000\text{ g/m}^3$ và có lúc đạt $2000\text{ g/m}^3$".

+-----------------------------------------------------------------------------------------+
|                                 BỐI CẢNH VÀ KHOẢNG TRỐNG KHOA HỌC                       |
+-----------------------------------------------------------------------------------------+
| • 13 triệu ha đất suy thoái tại VN | 1,2 triệu ha xói mòn trơ sỏi đá                    |
| • Mô hình kinh nghiệm (USLE, MUSLE): Giới hạn bước thời gian năm/tháng, thiếu bồi lắng  |
| • Mô hình quốc tế (WEPP, SWAT, ANSWERS): Mã nguồn đóng, chi phí cao, giả định dốc phẳng |
+-----------------------------------------------------------------------------------------+
+-----------------------------------------------------------------------------------------+
|                                  GIẢI PHÁP ĐỘT PHÁ CỦA LUẬN ÁN                          |
+-----------------------------------------------------------------------------------------+
| • Phát triển mô hình nhận thức vật lý phân bố (C++) giải hệ sóng động học 1D            |
| • Thuật toán sai phân hữu hạn Lax - Friedrichs cải tiến (trọng số không - thời gian)     |
| • Tích hợp GIS/DEM mô phỏng động học xói mòn liên rãnh, rãnh và truyền tải trong sông   |
| • Thử nghiệm thành công tại lưu vực Nậm Sập & Phiêng Hiềng (Sơn La)                     |
+-----------------------------------------------------------------------------------------+

Mặc dù có nhiều nỗ lực định lượng lượng đất mất bằng các mô hình kinh nghiệm truyền thống như USLE (Universal Soil Loss Equation) của Wischmeier & Smith (1978) hay MUSLE của Williams (1975), các công cụ này bộc lộ khoảng trống khoa học lớn: chỉ tính toán lượng đất xói mòn trung bình nhiều năm hoặc theo từng trận lũ đơn lẻ một cách tĩnh, không mô tả được diễn biến thủy động lực học theo thời gian thực ($t$), bỏ qua quá trình bồi lắng cục bộ và chuyển tải bùn cát trong mạng lưới sông suối, đồng thời thiếu khả năng phân rã động thái xói mòn theo từng bước thời gian ngắn (phút/giờ). Mặt khác, các mô hình nhận thức dựa trên quá trình vật lý tiên tiến của thế giới (như WEPP, EUROSEM, ANSWERS, SWAT) phần lớn là phần mềm thương mại mã nguồn đóng, đòi hỏi số liệu đầu vào cực kỳ chi tiết, nhiều tham số ngoại sinh được hiệu chỉnh trên điều kiện thổ nhưỡng ôn đới, áp dụng phương pháp giải tích giả định độ dốc sườn dốc đồng nhất nên khó thích ứng và gây sai số lớn khi áp dụng trực tiếp cho các lưu vực miền núi dốc đứng, phân cắt phức tạp ở Việt Nam.

Trước thực trạng đó, luận án tiến sĩ kỹ thuật thủy văn với đề tài "Nghiên cứu xây dựng mô hình toán mô phỏng dòng chảy và vận chuyển bùn cát trên lưu vực vừa và nhỏ" do nghiên cứu sinh Đào Tấn Quy thực hiện dưới sự hướng dẫn khoa học của PGS.TS Phạm Thị Hương Lan và PGS.TS Ngô Lê Long tại Trường Đại học Thủy Lợi đã giải quyết trọn vẹn khoảng trống trên. Luận án đặt ra hệ thống câu hỏi nghiên cứu và giả thuyết khoa học tường minh:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để xây dựng một khung toán học thủy động lực tích hợp mô tả chính xác quá trình hình thành dòng chảy tràn, xói mòn tách hạt (liên rãnh, rãnh) và vận chuyển bồi lắng bùn cát theo thời gian liên tục trên sườn dốc và mạng lưới sông suối?
  • Câu hỏi nghiên cứu 2 (RQ2): Lựa chọn và cải tiến sơ đồ số trị nào để giải hệ phương trình đạo hàm riêng hyperbolic phi tuyến mô phỏng dòng chảy và bùn cát nhằm triệt tiêu hiện tượng dao động số học và đảm bảo tính ổn định, hội tụ cao trên địa hình dốc phân cắt mạnh?
  • Câu hỏi nghiên cứu 3 (RQ3): Quy luật tương quan giữa cường độ mưa, độ dốc địa hình với động thái xói mòn liên rãnh và xói mòn rãnh tại các lưu vực vùng núi Tây Bắc Việt Nam diễn ra như thế nào?
  • Giả thuyết nghiên cứu 1 (H1): Việc kết hợp phương trình thấm Green-Ampt Mein-Larson, phương trình sóng động học 1D và phương trình liên tục bùn cát Bennett thông qua sơ đồ sai phân hữu hạn Lax-Friedrichs cải tiến bổ sung trọng số thời gian - không gian sẽ mô phỏng chuẩn xác đường quá trình lưu lượng $Q(t)$ và nồng độ bùn cát $C_s(t)$ cho từng trận mưa.
  • Giả thuyết nghiên cứu 2 (H2): Động lực xói mòn rãnh và liên rãnh trên nền đất dốc nhiệt đới có thể lượng hóa chính xác thông qua việc phân tích tương quan phi tuyến giữa năng lượng mưa, sức tải dòng chảy và các chỉ số địa hình trích xuất từ hệ thống thông tin địa lý (GIS).

Nghiên cứu tập trung vào phạm vi thực nghiệm tại hai lưu vực điển hình của tỉnh Sơn La: lưu vực Nậm Sập và lưu vực Phiêng Hiềng. Khung nghiên cứu đã làm chủ toàn bộ thuật toán, thiết kế mã nguồn phần mềm trên nền tảng C++, tạo ra công cụ tính toán thủy văn - môi trường độc lập, hiệu quả cao cho các lưu vực vừa và nhỏ tại Việt Nam.

Literature Review và Positioning

Nghiên cứu xói mòn và vận chuyển trầm tích trên thế giới đã trải qua quá trình tiến hóa từ các phương pháp thuần thực nghiệm sang các mô hình toán nhận thức vật lý phân bố (process-based distributed models).

TIẾN TRÌNH PHÁT TRIỂN CỦA CÁC LÝ THUYẾT VÀ MÔ HÌNH XÓI MÒN BÙN CÁT

[1940 - 1950: MÔ HÌNH KINH NGHIỆM BAN ĐẦU]
[1970 - 1980: KINH NGHIỆM PHỔ DỤNG & CẢI BIÊN]
[1980 - 1998: MÔ HÌNH NHẬN THỨC VẬT LÝ & THỦY ĐỘNG LỰC]
[ĐỊNH VỊ CỦA LUẬN ÁN - ĐÀO TẤN QUY (2017)]
      + RK4, thích ứng hoàn hảo cho địa hình dốc phi đồng nhất miền núi Việt Nam.

Trong lịch sử lý thuyết cơ chế xói mòn tồn tại hai trường phái học thuật đối trọng:

  1. Trường phái Ellison (1944, 1947): Khẳng định rằng sự phá hủy cấu trúc bề mặt và bong tách hạt đất chủ yếu do tác động động năng của hạt mưa rơi ($E = K \cdot V^{4.33}$); vì thế kiểm soát xói mòn bắt buộc phải tập trung vào lớp phủ thực vật và tán rừng nhằm triệt tiêu năng lượng giọt mưa.
  2. Trường phái Bennett (1939, 1974): Nhấn mạnh rằng sự hình thành rãnh xói và vận chuyển vật chất rắn diễn ra chủ yếu do động năng và sức cắt tiếp đáy của dòng chảy tràn tập trung; do đó kiểm soát xói mòn phải áp dụng các công trình cơ học nhằm phân tán dòng chảy, giảm độ dốc thủy lực mà không nhất thiết phải làm giảm tổng lượng dòng chảy mặt.

Về phương pháp mô hình hóa, các mô hình kinh nghiệm (Zingg, 1940; Musgrave, 1947; Wischmeier & Smith, 1978) đã đặt nền móng quản lý xói mòn đất nông nghiệp. Tuy nhiên, USLE chỉ áp dụng tính toán lượng mất đất trung bình nhiều năm (trên 20 năm) cho xói mòn bề mặt trên các sườn dốc thoải ($1 - 20%$), hoàn toàn mất tính chính xác trên các sườn dốc $>40%$ và không phản ánh được quá trình vận chuyển động học. Để khắc phục, các mô hình nhận thức như ANSWERS (Beasley et al., 1980), CREAMS (Knisel, 1980), AGNPS (Young et al., 1987), KINEROS (Woolhiser et al., 1990), WEPP (Flanagan & Nearing, 1995) và EUROSEM (Morgan et al., 1998) đã chuyển đổi sang tiếp cận hệ thống phương trình vi phân liên tục khối lượng và xung lượng.

So sánh vị thế học thuật của luận án với hai nghiên cứu quốc tế tiêu biểu:

  • So sánh với WEPP (USDA-ARS, Mỹ): WEPP sử dụng phương trình liên tục bùn cát của Foster & Meyer (1972) nhưng phần lớn áp dụng phương pháp giải tích trên các sườn dốc có độ dốc giả định đồng nhất, dẫn đến hạn chế lớn khi mô phỏng địa hình đồi núi thực tế có profile sườn dốc biến đổi phức tạp. Luận án của Đào Tấn Quy khắc phục điểm này bằng cách sử dụng phương pháp sai phân hữu hạn giải số trị trên mạng lưới ô lưới số hóa từ DEM, cho phép phản ánh chính xác sự phi đồng nhất của địa hình tự nhiên.
  • So sánh với EUROSEM (Morgan et al., Châu Âu): EUROSEM giải hệ phương trình sóng động học và bùn cát bằng sơ đồ sai phân hiện nhưng dễ gặp bất ổn định số học khi bước thời gian $\Delta t$ không thỏa mãn chặt chẽ điều kiện Courant-Friedrichs-Lewy (CFL) trên các triền dốc đứng nhiệt đới. Luận án đã giải quyết bài toán hội tụ bằng việc cải tiến sơ đồ sai phân Lax-Friedrichs có thêm các hệ số trọng số không gian và thời gian ($\theta_x, \theta_t$), tăng cường độ ổn định và giảm thiểu hiện tượng tiêu tán số (numerical diffusion).
  • So sánh với các nghiên cứu trong nước: Các nghiên cứu của Nguyễn Trọng Hà (1996), Phạm Ngọc Dũng (1996), Lại Vĩnh Cẩm (2002), Phạm Hữu Tỵ (2009), Nguyễn Văn Dũng (2011) chủ yếu dừng lại ở việc áp dụng USLE/RUSLE kết hợp GIS ở quy mô tĩnh; hay nghiên cứu của Lê Mạnh Hùng (2011) dùng SWAT đánh giá bùn cát lưu vực Mekong (đạt NSE = 0.86, PBIAS = 11.8% tại Kratie với tải lượng bùn cát 162 triệu tấn/năm) nhưng SWAT mang tính thương mại, tích hợp sẵn và sử dụng MUSLE cho từng tiểu lưu vực chứ không giải trực tiếp thủy động lực học sóng dòng chảy mặt theo từng phút. Công trình của Đào Tấn Quy là bước đột phá tự chủ hoàn toàn về mặt toán học và thuật toán mô phỏng thủy văn - bùn cát chi tiết tại Việt Nam.

Đó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à kế thừa nền tảng lý thuyết phương trình liên tục bùn cát một chiều của Bennett (1974), Foster & Meyer (1972) bằng cách thiết lập một khung toán học liên kết tương hỗ giữa quá trình thấm biến đổi, dòng chảy mặt sóng động học và động học vận chuyển bùn cát đa nguồn (liên rãnh, rãnh và dòng sông chính).

+-----------------------------------------------------------------------------------------+
|                    KHUNG LÝ THUYẾT VÀ HỆ PHƯƠNG TRÌNH LIÊN TỤC TÍCH HỢP                 |
+-----------------------------------------------------------------------------------------+
| 1. Thấm Green - Ampt Mein - Larson:                                                     |
|    f_i = K_e * [ 1 + ( (phi - theta_0) * psi_c ) / I ]                                  |
|                                                                                         |
| 2. Dòng chảy mặt (Sóng động học 1D):                                                    |
|    ∂h/∂t + ∂(alpha * h^m)/∂x = r(t)                                                     |
|    Biên & Ban đầu: h(0,t) = 0 ; h(x,0) = 0                                              |
|                                                                                         |
| 3. Dòng chảy trong kênh/sông:                                                           |
|    ∂A/∂t + ∂Q/∂x = q_A                                                                  |
|                                                                                         |
| 4. Vận chuyển bùn cát trên sườn dốc (Bennett, 1974):                                    |
|    ∂(C_s * h)/∂t + ∂q_s/∂x = theta(x, t)                                                |
|                                                                                         |
| 5. Giải thuật vi phân lưu trữ dòng bùn cát (Runge - Kutta bậc 4):                       |
|    dq(t)/dt = (1/k) * [ p(t) - q(t) ]  ==> f(t, q(t))                                   |
+-----------------------------------------------------------------------------------------+

Khung lý thuyết của luận án làm sáng tỏ hai cơ chế vận động vật chất:

  1. Cơ chế xói mòn liên rãnh (Inter-rill Erosion): Năng lượng hạt mưa rơi làm phân rã các liên kết hạt đất, tạo ra các phần tử bùn cát lơ lửng được lớp dòng chảy tràn mỏng chuyển tải vào hệ thống khe rãnh: $$D_i = K_i \cdot I_e^2 \cdot G_e \cdot C_e \cdot S_f$$ Trong đó, $K_i$ là tính xói mòn liên rãnh, $I_e$ là cường độ mưa hiệu quả, $G_e, C_e$ là hệ số điều chỉnh lớp phủ, và $S_f = 1.05 - 0.85\exp(-4\sin\theta)$ là nhân tố điều chỉnh ảnh hưởng của góc dốc $\theta$.
  2. Cơ chế xói mòn rãnh (Rill Erosion): Dòng chảy tràn khi tập trung vào các rãnh xói cục bộ sinh ra ứng suất tiếp xúc đáy $\tau = \gamma h S_0$. Khi $\tau > \tau_c$ (ứng suất cắt tới hạn), xói mòn rãnh được kích hoạt; tốc độ xói mòn phụ thuộc vào sức tải bùn cát cực đại $T_c$ của dòng nước so với nồng độ bùn cát hiện hữu $q_s$: $$D_r = K_r (\tau - \tau_c) \left(1 - \frac{q_s}{T_c}\right)$$

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

Mô hình phân tích của luận án tích hợp 4 cấu phần vật lý cốt lõi:

  • Mô đun thấm khí tượng: Giải phương trình Green-Ampt Mein-Larson (1973) tính toán cường độ thấm thực $f_i$ dựa trên độ dẫn thủy lực bão hòa $K_e$, độ ẩm ban đầu $\theta_0$, độ xốp $\phi$, thế áp lực mao dẫn sườn ẩm $\psi_c$ và độ sâu thấm lũy tích $I$: $$f_i = K_e \left[ 1 + \frac{(\phi - \theta_0)\psi_c}{I} \right]$$
  • Mô đun thủy lực dòng chảy mặt: Ứng dụng phương trình sóng động học kết hợp công thức ma sát Manning ($u = \alpha h^{m-1}$ với $\alpha = \frac{1}{n} S_0^{1/2}$ và $m = 5/3$).
  • Mô đun xói mòn - bồi lắng sườn dốc: Mô phỏng sự biến thiên nồng độ bùn cát $C_s(x,t)$ và lượng bùn cát bồi lắng/tách hạt $\theta(x,t)$ dọc theo phương dòng chảy $x$.
  • Mô đun diễn toán dòng chảy và bùn cát lòng kênh: Sử dụng phương trình Runge-Kutta bậc 4 (RK4) để giải phương trình vi phân lưu trữ dòng bùn cát: $$\frac{dq(t)}{dt} = \frac{1}{k}[p(t) - q(t)]$$ Đảm bảo tính liên tục của sóng bùn cát truyền từ toàn bộ các tiểu lưu vực về mặt cắt cửa ra.

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ủ thế giới quan thực chứng (positivism) và chủ nghĩa hiện thực phản biện (critical realism), xây dựng hệ thống mô phỏng tất định (deterministic process-based model) mô tả chi tiết các quá trình vật lý tự nhiên. Nghiên cứu triển khai thiết kế đa tầng bậc (multi-level design):

  • Tầng 1 (Ô lưới lưu vực - Mesh level): Lưu vực được số hóa thành lưới các ô nguyên tố chữ nhật ($i, j$) từ mô hình số độ cao DEM, đồng nhất các thông số vi mô về thổ nhưỡng, độ dốc, thảm phủ.
  • Tầng 2 (Tiểu lưu vực - Sub-catchment level): Tập hợp các dòng chảy tràn từ vùng liên rãnh dồn vào rãnh xói nhánh.
  • Tầng 3 (Mạng lưới kênh/sông - Channel network level): Tổng hợp lượng gia nhập bên ($q_A$), diễn toán truyền lũ và bùn cát về tuyến cửa ra khống chế.
SƠ ĐỒ PHÂN CẤP TÍNH TOÁN ĐA TẦNG BẬC CỦA MÔ HÌNH TOÁN

Quy trình nghiên cứu rigorous

Trọng tâm cải tiến phương pháp số học của luận án là việc phát triển sơ đồ sai phân hữu hạn Lax-Friedrichs có bổ sung trọng số không gian ($\theta_x$) và thời gian ($\theta_t$). Phương trình vi phân đạo hàm riêng của sóng động học: $$\frac{\partial h}{\partial t} + \frac{\partial (\alpha h^m)}{\partial x} = r$$

Được rời rạc hóa trên lưới sai phân $(j, n)$ theo công thức cải tiến: $$h_j^{n+1} = \theta_t \left[ \theta_x h_{j+1}^n + (1 - \theta_x) h_{j-1}^n \right] - \frac{\Delta t}{2\Delta x} \left[ \alpha (h_{j+1}^n)^m - \alpha (h_{j-1}^n)^m \right] + \Delta t \cdot r_j^n$$

SƠ ĐỒ SAI PHÂN LAX - FRIEDRICHS CẢI TIẾN TRONG KHÔNG - THỜI GIAN

      Thời gian (t)
                   j-1         j         j+1

Sơ đồ này khắc phục triệt để tính bất ổn định của sơ đồ sai phân hiện thuần túy và hiện tượng khuếch tán số học quá mức của sơ đồ Lax cổ điển khi áp dụng trên các lưu vực đồi núi có độ dốc thay đổi đột ngột. Toàn bộ thuật toán được đóng gói trong phần mềm lập trình hướng đối tượng bằng C++, gồm các module độc lập: Module nhập dữ liệu địa hình/thổ nhưỡng/mưa, Module giải thủy lực mặt đệm, Module xói mòn bùn cát sườn dốc, Module diễn toán thủy lực lòng dẫn và Module trích xuất biểu đồ quá trình $Q(t), C_s(t)$.

Data và phân tích

Nghiên cứu ứng dụng công nghệ viễn thám và GIS (ArcGIS, MapInfo) để trích xuất hệ thống dữ liệu không gian và thuộc tính vật lý tại hai lưu vực thử nghiệm:

  • Lưu vực Nậm Sập (Sơn La): Lưu vực sông miền núi điển hình có diện tích lưu vực vừa và nhỏ, độ dốc lớn, mạng lưới sông suối dốc đứng, rừng phòng hộ xen kẽ đất nương đồi canh tác ngô, sắn.
  • Lưu vực Phiêng Hiềng (Sơn La): Lưu vực có độ che phủ thực vật biến động mạnh, mức độ canh tác nông nghiệp trên đất dốc cao, nguy cơ xói mòn rửa trôi đặc biệt nghiêm trọng trong mùa lũ.
BẢNG THÔNG SỐ VẬT LÝ VÀ ĐIỀU KIỆN ĐẦU VÀO TẠI CÁC LƯU VỰC THỬ NGHIỆM

Mô hình được hiệu chỉnh và kiểm định thông qua tập số liệu các trận mưa lũ thực đo chi tiết. Đánh giá độ tin cậy mô hình thông qua đường quá trình lưu lượng $Q(t)$ và đường quá trình nồng độ bùn cát $C_s(t)$, kiểm chứng sai số bằng hệ số Nash-Sutcliffe Efficiency ($NSE > 0.80$) và độ nhạy của các thông số nhám Manning ($n$), độ dẫn thủy lực bão hòa ($K_e$), hệ số xói mòn đất ($K_i, K_r$).

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

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

  1. Thiết lập thành công phương trình tương quan định lượng giữa xói mòn và các nhân tố động lực: Luận án đã xây dựng phương trình hồi quy phi tuyến biểu diễn mối quan hệ chặt chẽ giữa lượng xói mòn liên rãnh ($E_{lr}$), xói mòn rãnh ($E_r$) với cường độ mưa ($I$) và độ dốc lưu vực ($S$) tại Nậm Sập và Phiêng Hiềng: $$E_{lr} = a \cdot I^b \cdot S^c$$ $$E_r = \alpha \cdot I^\beta \cdot S^\gamma$$ Kết quả chứng minh rằng cường độ mưa đóng vai trò lũy thừa bậc cao đối với sự bùng phát xói mòn rãnh, vượt trội hơn ảnh hưởng đơn thuần của tổng lượng mưa.
  2. Sự tương thích vượt trội của thuật toán sai phân Lax-Friedrichs cải tiến: Việc đưa trọng số không gian và thời gian ($\theta_x, \theta_t$) đã triệt tiêu hoàn toàn các sóng dao động số học giả tạo tại các điểm gián đoạn đáy và vùng sườn dốc chuyển tiếp, bảo toàn khối lượng nước và bùn cát với độ chính xác cao.
  3. Hiện tượng trễ pha (hysteresis) giữa đỉnh lũ và đỉnh bùn cát: Dữ liệu mô phỏng và thực nghiệm tại Nậm Sập và Phiêng Hiềng đã tái hiện chính xác hiện tượng nồng độ bùn cát $C_s(t)$ đạt cực đại trước hoặc đồng thời với đỉnh lưu lượng dòng chảy $Q(t)$ tùy thuộc vào vị trí nguồn cung cấp bùn cát trên sườn dốc gần hay xa mạng lưới sông suối.
  4. Định lượng vai trò kiểm soát xói mòn của thảm phủ thực vật: Phân tích độ nhạy chứng minh rằng khi độ che phủ tầng tán và thảm mục tăng từ $20%$ lên $70%$, cường độ xói mòn liên rãnh giảm từ $3.5$ đến $5.8$ lần, khẳng định vai trò của tầng tán rừng và thảm cỏ mục trong việc giảm động năng mưa rơi và tăng dung lượng thấm ngấm.

Implications đa chiều

SƠ ĐỒ HỆ THỐNG HÓA CÁC GIÁ TRỊ ỨNG DỤNG ĐA CHIỀU
  • Về mặt học thuật và phương pháp: Chuyển đổi phương thức nghiên cứu xói mòn tại Việt Nam từ các công thức kinh nghiệm bán định lượng sang mô hình số trị thủy động lực học phân bố theo thời gian thực. Cung cấp nền tảng mã nguồn C++ mở để các nhà khoa học tiếp tục tích hợp các mô đun chất lượng nước và ô nhiễm nông nghiệp.
  • Về mặt thực tiễn kỹ thuật: Cung cấp bộ công cụ tính toán chính xác bùn cát bồi lắng hồ chứa thủy lợi, thủy điện trên các lưu vực sông suối miền núi phía Bắc, tối ưu hóa dung tích chết của các công trình hồ chứa bậc thang tại Sơn La.
  • Về mặt chính sách và quản trị: Cung cấp cơ sở khoa học định lượng giúp các nhà hoạch định chính sách tại Sở Nông nghiệp & PTNT, Sở Tài nguyên & Môi trường tỉnh Sơn La phân vùng các khu vực đất dốc xung yếu, ban hành các quy định về khung canh tác bền vững (ruộng bậc thang, băng cây chắn cỏ, mô hình nông lâm kết hợp) nhằm hạn chế mất đất và suy thoái môi trường.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn chỉ ra các giới hạn nghiên cứu:

  • Giới hạn về dữ liệu đo đạc phân giải cao: Mô hình đòi hỏi chuỗi số liệu mưa tự ghi theo từng trận mưa (phút/giờ) và số liệu quan trắc nồng độ bùn cát liên tục. Tại nhiều lưu vực vùng sâu vùng xa của Việt Nam, mạng lưới trạm đo thủy văn còn thưa thớt, gây khó khăn cho việc thiết lập thông số đầu vào.
  • Mô hình hóa dòng chảy kênh 1D: Luận án tiếp cận dòng chảy và bùn cát trong kênh theo phương pháp sóng động học 1D. Điều này phù hợp với các sông suối dốc miền núi nhưng chưa phản ánh được hoàn toàn hiện tượng dòng thứ sinh, xói lở cục bộ hai bên bờ sông hoặc bồi lắng lòng dẫn 2D/3D ở các đoạn sông uốn khúc phức tạp.
  • Chưa tích hợp trượt lở khối (Mass Wasting): Mô hình tập trung vào xói mòn mặt, xói mòn liên rãnh và rãnh do nước chảy tràn, chưa liên kết với các quá trình tai biến trượt lở đất đá quy mô lớn do bão lũ kích hoạt.

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

  1. Nâng cấp mô đun thủy lực dòng dẫn từ sóng động học 1D lên hệ phương trình Saint-Venant 2D đầy đủ nhằm mô phỏng vùng ngập lụt ven sông và bồi lắng lòng hồ chi tiết.
  2. Tích hợp cơ chế liên kết giữa mô hình xói mòn thủy lực bề mặt với mô hình ổn định mái dốc địa kỹ thuật (geotechnical slope stability) để dự báo trượt lở đất và dòng bùn đá (debris flow).
  3. Tích hợp các kịch bản biến đổi khí hậu toàn cầu (IPCC) và kịch bản chuyển đổi mục đích sử dụng đất để dự báo nguy cơ mất đất trong tương lai dài hạn (2030 – 2050).
  4. Tự động hóa hiệu chỉnh thông số mô hình bằng các thuật toán trí tuệ nhân tạo (AI), giải thuật di truyền (Genetic Algorithm) và lọc hạt (Particle Filter).

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

  • Tác động học thuật (Academic Impact): Công trình là tài liệu tham khảo mẫu mực trong lĩnh vực mô hình hóa thủy văn và xói mòn bùn cát tại Việt Nam, mở ra hướng nghiên cứu làm chủ công nghệ mô hình số trị thay vì phụ thuộc hoàn toàn vào các phần mềm thương mại ngoại nhập.
  • Chuyển giao và ứng dụng thực tiễn: Cung cấp trực tiếp cho các cơ quan tư vấn thiết kế thủy lợi, thủy điện và quy hoạch tài nguyên nước công cụ giám sát, đánh giá bồi lắng và xói mòn lưu vực tại tỉnh Sơn La cũng như các tỉnh miền núi phía Bắc.
  • Bảo tồn sinh thái và xã hội: Góp phần giảm thiểu thiệt hại do lũ quét, lũ bùn đá cho cộng đồng dân cư sinh sống tại các thung lũng miền núi, bảo vệ độ phì nhiêu của đất sản xuất nông nghiệp, đảm bảo an ninh lương thực và phát triển sinh kế bền vững cho đồng bào vùng cao.

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

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 mở rộng hệ phương trình liên tục bùn cát của Bennett (1974) và hệ phương trình sóng động học 1D thông qua việc thiết lập sơ đồ liên kết thủy động lực liên tục giữa quá trình thấm biến đổi theo Green-Ampt Mein-Larson, quá trình tách hạt liên rãnh do mưa, xói mòn rãnh do lực cắt dòng chảy và truyền sóng bùn cát trong mạng lưới kênh dẫn giải bằng phương pháp Runge-Kutta bậc 4 (RK4).

2. Đột phá về mặt phương pháp nghiên cứu so với các công trình trước đây là gì?

Khác với WEPP hay SWAT (chủ yếu giải tích hoặc bán thực nghiệm trên sườn dốc đồng nhất) và EUROSEM (dễ phân kỳ trên độ dốc cao), luận án đã phát triển thành công thuật toán sai phân hữu hạn Lax-Friedrichs cải tiến có tích hợp các hệ số trọng số không gian ($\theta_x$) và thời gian ($\theta_t$). Đột phá này đảm bảo tính ổn định tuyệt đối của nghiệm số, triệt tiêu dao động số học và phản ánh trung thực địa hình dốc phân cắt phức tạp trên nền tảng ô lưới GIS.

3. Phát hiện thực nghiệm nào đáng ngạc nhiên nhất từ dữ liệu nghiên cứu?

Phát hiện thực nghiệm quan trọng nhất là mức độ phi tuyến tính cực cao của xói mòn rãnh đối với cường độ mưa thời đoạn ngắn ($I_{30}$) so với tổng lượng mưa. Tại lưu vực Nậm Sập và Phiêng Hiềng, xói mòn rãnh chỉ bùng phát khi cường độ mưa vượt qua một ngưỡng thấm bão hòa cục bộ xác định, tại đó ứng suất tiếp xúc dòng chảy vượt giá trị cắt tới hạn $\tau_c$, tạo nên sự gia tăng đột biến về nồng độ bùn cát cửa ra theo dạng hàm mũ.

4. Luận án có cung cấp đầy đủ giao thức để tái lặp (replication protocol) không?

Có. Luận án trình bày chi tiết toàn bộ hệ thống phương trình vi phân gốc, điều kiện biên ($h(0,t) = 0$), điều kiện ban đầu ($h(x,0) = 0$), thuật toán rời rạc hóa sai phân, sơ đồ khối logic chương trình, cấu trúc các module hàm C++ và bảng thông số thổ nhưỡng, thủy lực thực nghiệm tại Nậm Sập và Phiêng Hiềng, cho phép các nhà nghiên cứu độc lập lập trình và tái lập hoàn toàn kết quả mô phỏng.

5. Chương trình nghị sự nghiên cứu 10 năm tiếp theo được định hình như thế nào?

Lộ trình phát triển 10 năm bao gồm:

  • Mở rộng mô hình lên không gian 2 chiều (2D Hydrodynamics) cho các vùng bãi bồi và hồ chứa.
  • Tích hợp tương hỗ với mô hình ổn định mái dốc dự báo lũ quét, trượt lở đất và lũ bùn đá.
  • Tự động hóa hiệu chỉnh thông số qua nền tảng dữ liệu lớn (Big Data) và cảm biến viễn thám thời gian thực.
  • Ứng dụng mô phỏng chất lượng nước và vận chuyển hóa chất nông nghiệp (N, P, thuốc BVTV) bám dính theo hạt bùn cát.

Kết luận

Luận án tiến sĩ của tác giả Đào Tấn Quy là một công trình khoa học công phu, nghiêm túc và có hàm lượng học thuật cao. Tóm tắt 5 đóng góp cốt lõi của công trình:

  1. Hệ thống hóa toàn diện cơ sở khoa học: Tổng quan sâu sắc các trường phái lý thuyết và mô hình toán xói mòn - vận chuyển bùn cát trong và ngoài nước, chỉ rõ những tồn tại kỹ thuật của các mô hình kinh nghiệm và phần mềm ngoại nhập.
  2. Làm chủ khung toán học thủy động lực phân bố: Thiết lập thành công hệ phương trình liên tục dòng chảy và bùn cát tích hợp cơ chế thấm Green-Ampt Mein-Larson, sóng động học mặt phẳng và động lực xói mòn liên rãnh - rãnh.
  3. Đột phá giải thuật số trị: Đề xuất và ứng dụng thành công sơ đồ sai phân hữu hạn Lax-Friedrichs cải tiến có trọng số không - thời gian kết hợp giải thuật Runge-Kutta bậc 4 (RK4), giải quyết triệt để bài toán ổn định và hội tụ trên địa hình đồi núi dốc.
  4. Tự chủ công nghệ phần mềm: Thiết kế và xây dựng hoàn chỉnh chương trình tính toán bằng ngôn ngữ C++ với cấu trúc module hướng đối tượng, tích hợp dữ liệu không gian viễn thám và GIS.
  5. Xác lập quy luật thực nghiệm địa phương: Thiết lập hệ thống phương trình tương quan định lượng giữa cường độ mưa, độ dốc với xói mòn liên rãnh và rãnh cho các lưu vực Nậm Sập và Phiêng Hiềng (Sơn La), cung cấp công cụ dự báo và giải pháp thực tiễn quản lý đất dốc bền vững.

Công trình không chỉ nâng tầm năng lực nghiên cứu mô hình toán thủy văn tại Việt Nam mà còn là đóng góp khoa học quan trọng, mở ra hướng ứng dụng rộng rãi trong quy hoạch tài nguyên nước, bảo vệ môi trường và phòng chống thiên tai tại các vùng miền núi trên cả nước.