Tổng quan về luận án

Nghiên cứu cơ chế thủy động lực học và quá trình lan truyền mặn tại vùng cửa sông ven biển giữ vị trí chiến lược sống còn đối với sự phát triển kinh tế - xã hội và bảo đảm an ninh lương thực quốc gia tại Đồng bằng sông Cửu Long (ĐBSCL). Với diện tích tự nhiên khoảng 4,0 triệu ha, ĐBSCL đóng góp 50% sản lượng lúa, 65% sản lượng thủy sản nuôi trồng, 70% sản lượng trái cây và 95% lượng gạo xuất khẩu của cả nước. Tuy nhiên, toàn vùng có tới 1,4 đến 2,1 triệu ha đất tự nhiên chịu ảnh hưởng trực tiếp của xâm nhập mặn trong mùa kiệt do mạng lưới sông ngòi tiếp giáp với bờ biển dài hơn 700 km chịu tác động đồng thời của chế độ bán nhật triều không đều Biển Đông và nhật triều Biển Tây. Trong các năm hạn mặn cực đoan lịch sử 2015-2016 và 2019-2020: "Trên sông Tiền và sông Hậu, độ mặn là trên 4,5‰, xâm nhập sâu tới 70 km tính từ cửa sông, thậm chí có nơi lên đến 85 km" (Viện Khoa học Thủy lợi miền Nam, 2016), gây tê liệt hệ thống cấp nước sinh hoạt và sản xuất nông nghiệp.

Khoảng trống nghiên cứu then chốt (research gap) nằm ở chỗ hầu hết các công trình nghiên cứu và quy hoạch kiểm soát mặn tại ĐBSCL từ trước đến nay đều dựa trên các chuỗi số liệu quan trắc tầng mặt đơn điểm kết hợp với mô hình thủy lực - truyền chất một chiều (1D), vốn giả định nồng độ mặn được xáo trộn đồng nhất và đại diện bằng giá trị trung bình trên toàn bộ mặt cắt ướt. Cách tiếp cận truyền thống này bỏ qua hoàn toàn cấu trúc phân tầng mật độ ba chiều (3D density stratification), sự hình thành nêm mặn (salt wedge) và các động lực vi mô tại mặt phân cách nước ngọt - mặn. Luận án tiến sĩ kỹ thuật của NCS Đỗ Đắc Hải thuộc chuyên ngành Kỹ thuật Tài nguyên Nước (Mã ngành: 9 58 02 12) thực hiện dưới sự hướng dẫn của TS. Huỳnh Thanh Sơn và TSKH. Nguyễn Ân Niên tại Viện Khoa học Thủy lợi Việt Nam đã tiên phong giải quyết khoảng trống này thông qua ba câu hỏi nghiên cứu và giả thuyết tương ứng:

  1. Câu hỏi nghiên cứu 1 (RQ1): Cơ chế thủy động lực học nào chi phối sự hình thành, cấu trúc không gian và sự vận động của nêm mặn dọc theo các phân lưu cửa sông Cửu Long trong chu kỳ triều mùa kiệt? Giả thuyết 1 (H1): Chênh lệch tỷ trọng giữa khối nước ngọt thượng nguồn và khối nước biển dẫn đến sự phân tầng mật độ mạnh hoặc phân tầng một phần tại vùng cửa sông, tạo ra chênh lệch nồng độ mặn thẳng đứng $\Delta S = S_{bot} - S_{surf} > 0$ biến thiên theo pha triều.
  2. Câu hỏi nghiên cứu 2 (RQ2): Ngưỡng phân tầng thực nghiệm và định lượng động lực nêm mặn trên 7 nhánh sông chính của ĐBSCL được xác lập như thế nào? Giả thuyết 2 (H2): Hệ số phân tầng Simmons cải tiến $\eta$ và chỉ số Pritchard $n$ phản ánh chính xác ranh giới chuyển tiếp giữa trạng thái phân tầng mạnh, xáo trộn từng phần và đồng nhất trên toàn hệ thống sông.
  3. Câu hỏi nghiên cứu 3 (RQ3): Cấu trúc phân tầng mặn 3D cung cấp cơ sở kỹ thuật nào để tái cấu trúc công nghệ công trình và lịch trình khai thác nước ngọt tối ưu? Giả thuyết 3 (H3): Khai thác nước tầng mặt theo cửa van phân tầng và trạm bơm di động trong các khoảng thời gian chân triều cho phép thu hồi nước ngọt chất lượng cao ($S < 4$ g/l) ngay tại các vùng không gian vốn bị coi là nhiễm mặn hoàn toàn theo tính toán 1D.

Phạm vi nghiên cứu bao quát 7 cửa sông chính vùng ĐBSCL (Cửa Tiểu, Cửa Đại, Hàm Luông, Cổ Chiên, Cung Hầu, Định An, Trần Đề), tập trung chuyên sâu vào hạ lưu sông Hậu qua hai nhánh Định An và Trần Đề trên chiều dài hơn 80 km. Nghiên cứu tích hợp chuỗi số liệu thực đo chuyên sâu mùa khô 2015-2016 với mô phỏng thủy lực 3 chiều phi thủy tĩnh, tạo bước đột phá định lượng hóa tài nguyên nước ngọt tầng mặt phục vụ ứng phó biến đổi khí hậu và nước biển dâng.

Literature Review và Positioning

Khung lý thuyết về thủy động lực học cửa sông khởi nguồn từ định nghĩa kinh điển của Pritchard (1955): "Estuary là một thuỷ vực nửa kín ven bờ thông với biển khơi, trong đó có sự hoà trộn nhất định giữa nước biển mặn và nước ngọt chảy ra từ lục địa". Pritchard (1955) và sau đó là Cameron và Pritchard (1963) đã đặt nền móng phân loại cửa sông dựa trên cấu trúc phân tầng độ mặn thẳng đứng và chế độ tuần hoàn trọng lực thành ba dạng cơ bản: cửa sông phân tầng mạnh (highly stratified / salt wedge), cửa sông xáo trộn một phần (partially mixed) và cửa sông xáo trộn hoàn toàn (well-mixed / homogeneous). Leeder (1984) phát triển hệ thống phân loại này thành 4 nhóm hình thái động lực tương ứng với cân bằng giữa năng lượng triều, dòng chảy sông và quá trình bồi xói trầm tích.

Lịch sử nghiên cứu thế giới ghi nhận hai trường phái đối thoại học thuật sâu sắc về động lực học hòa trộn cửa sông:

  • Trường phái dòng chảy hai lớp và nêm mặn giải tích: Schijf và Schönfeld (1953), Farmer và Morgan (1953), Keulegan (1957) tiếp cận hiện tượng nêm mặn thông qua hệ phương trình vi phân liên tục và động lượng cho hệ thống hai lớp chất lỏng có tỷ trọng khác nhau, xác định chiều dài nêm mặn $L_0/H_0$ phụ thuộc vào số Reynolds ($R$) và số Froude mật độ ($Fr_0$). Trường phái này nhấn mạnh tính bền vững của mặt phân cách khi gradient tỷ trọng triệt tiêu sự rối thẳng đứng.
  • Trường phái bất ổn định thủy động và xáo trộn mãnh liệt (IMP): West và Mangat (1986), Thorpe (1987), Geyer và Farmer (1989), Lewis (1996) chứng minh rằng mặt phân cách hai lớp không cố định mà liên tục chịu tác động của sóng nội và hiện tượng trượt vận tốc. Khi số Richardson gradient $Ri_g = -\frac{g}{\rho}\frac{\partial \rho / \partial z}{(\partial u / \partial z)^2} < 0.25$, bất ổn định Kelvin-Helmholtz (K-H) bùng phát, phá vỡ cấu trúc hai lớp và tạo ra các chu kỳ xáo trộn mạnh (Intense Mixing Periods - IMP) khiến cột nước chuyển hóa nhanh chóng từ phân tầng sang đồng nhất trong các pha triều rút.

So sánh với các nghiên cứu quốc tế điển hình:

  1. Hệ thống sông Mississippi (Hoa Kỳ): Nghiên cứu của Harleman (1961) và US Army Corps of Engineers (1988) xác định Mississippi là cửa sông phân tầng mạnh điển hình vào mùa kiệt khi lưu lượng giảm xuống dưới 2.831 m³/s, nêm mặn xâm nhập sâu 160 km tới New Orleans với tốc độ 3-5 km/ngày. Giải pháp đột phá được áp dụng thành công là xây dựng đập ngầm ngăn mặn (underwater sill) dài 500 m, cao trình đỉnh đập cách mặt nước 14,1 m tại vị trí cách biển 100 km, chặn dòng nêm mặn tầng đáy mà vẫn đảm bảo mớn nước cho tàu biển trọng tải lớn lưu thông trên tầng mặt (Hình 1-19).
  2. Cửa sông Fraser (Canada): Geyer và Farmer (1989) quan trắc thực địa xác nhận chênh lệch độ mặn giữa tầng mặt và đáy lên tới 26,0 g/l trong pha triều rút, khẳng định cơ chế bất ổn định trượt K-H chi phối sự tiêu tán năng lượng và gia tăng hệ số khuếch tán dọc $K_{ze}$.
  3. Cửa sông Tees và Tamar (Anh): New và cộng sự (1986), Lewis (1996) ghi nhận các đợt IMP xuất hiện từ 1,5 đến 5,0 giờ sau thời điểm nước ngưng đỉnh triều (HWS), làm phẳng hoàn toàn gradien độ mặn phương đứng.

Luận án của NCS Đỗ Đắc Hải định vị chính xác điểm giao thoa giữa các lý thuyết thủy động lực học nêm mặn quốc tế với bối cảnh đặc thù của vùng đồng bằng châu thổ nhiệt đới gió mùa hạ lưu sông Mê Công, nơi có mạng lưới kênh rạch chằng chịt, biên độ triều lớn (trên 3,5 m tại Biển Đông) và dòng chảy thượng lưu bị biến đổi phức tạp bởi các bậc thang thủy điện xuyên biên giới sau năm 2012.

Đó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à kiểm chứng các lý thuyết thủy động lực học phân tầng mật độ trong điều kiện cửa sông đồng bằng thấp:

  1. Mở rộng lý thuyết phân loại cửa sông của D. Pritchard (1955) thông qua chỉ số phân tầng không thứ nguyên $n$: $$n = \frac{S_{bot} - S_{surf}}{S_m}$$ Trong đó $S_{bot}$ là độ mặn đáy, $S_{surf}$ là độ mặn mặt, và $S_m$ là độ mặn trung bình mặt cắt đứng. Luận án thiết lập miền giá trị giới hạn thực tế cho các nhánh sông ĐBSCL, chứng minh trạng thái chuyển đổi linh hoạt giữa phân tầng một phần ($0.1 \le n \le 0.5$) và phân tầng mạnh ($n > 0.5$) theo từng pha triều.
  2. Ứng dụng và chuẩn hóa lý thuyết tham số phân tầng Simmons cải tiến của GS.TSKH Nguyễn Ân Niên: $$\eta = \frac{W^+}{W^-}$$ Trong đó $W^+$ là tổng lượng dòng chảy sông trong một pha triều xuống và $W^-$ là lượng dòng chảy ngược (dòng triều) trong pha triều lên. Luận án chứng minh hệ số $\eta$ phản ánh trung thực tương quan động lực giữa lực đẩy nước ngọt và lực đẩy triều mặn dọc theo 7 nhánh sông Cửu Long hơn hẳn tham số thể tích lăng trụ triều cổ điển $\alpha = Q_m \cdot T_s / P_t$.
  3. Xác lập luận điểm khoa học mới: "Trên các nhánh sông Cửu Long có xuất hiện sự phân tầng mặn vào một số thời điểm trong mùa kiệt" (Đỗ Đắc Hải, 2022). Phát hiện này bác bỏ quan điểm truyền thống cho rằng các cửa sông ĐBSCL luôn ở trạng thái xáo trộn hoàn toàn do biên độ triều cao, từ đó tạo ra bước chuyển đổi hệ hình (paradigm shift) từ thủy lực truyền chất đơn tầng (1D) sang thủy động lực học cấu trúc mật độ đa tầng (3D).

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

Khung phân tích của luận án tích hợp ba trụ cột lý thuyết: (1) Động lực học chất lỏng hai lớp Schijf - Schönfeld, (2) Lý thuyết bất ổn định ranh giới mặt phân cách Kelvin-Helmholtz, và (3) Mô hình trường phân bố độ mặn ba chiều (3D Baroclinic Salinity Dispersion Model).

                                [Nguồn động lực biên]
               +--------------------------+--------------------------+
               | Thượng nguồn: Lưu lượng Q | Hạ lưu: Dao động triều H |
               +--------------------------+--------------------------+
                                          |
                                          v
                [Khung phân tích thủy động lực phân tầng mật độ 3D]
               +-----------------------------------------------------+
               | 1. Tham số kinh nghiệm: Simmons cải tiến (η) &      |
               |    Pritchard (n) xác lập trạng thái phân tầng       |
               | 2. Phương trình trạng thái mật độ ρ(S, T, Css)      |
               | 3. Cân bằng động lượng & bất ổn định mặt phân cách  |
               |    theo số Richardson (Ri) và số Froude mật độ (Fr) |
               +-----------------------------------------------------+
                                          |
                                          v
             [Mô phỏng số 3 chiều MIKE 3 FM (Non-hydrostatic/Baroclinic)]
               +-----------------------------------------------------+
               | - Tái hiện hình thái nêm mặn, độ dốc mặt phân cách  |
               | - Phân bố nồng độ mặn theo 3 trục (x, y, z) & pha t |
               +-----------------------------------------------------+
                                          |
                                          v
               [Định hướng giải pháp công trình & vận hành hớt ngọt]
               +-----------------------------------------------------+
               | - Khai thác nước ngọt tầng mặt (S < 4 g/l)          |
               | - Cửa van nhiều lớp / Trạm bơm nổi cơ động          |
               | - Đập ngầm ngăn mặn tầng đáy theo pha triều         |
               +-----------------------------------------------------+

Khung phân tích xác định rõ các điều kiện biên giới hạn (boundary conditions): áp dụng cho vùng cửa sông hình phễu ven biển chịu ảnh hưởng bán nhật triều biên độ 2,5 - 4,0 m, địa hình đáy biến đổi từ 0 m đến -25 m, lưu lượng kiệt thượng nguồn sông Hậu dao động từ 1.500 m³/s đến 3.500 m³/s, nồng độ muối nước biển dao động từ 30‰ đến 34‰.

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ chặt chẽ thế giới quan thực chứng luận (positivism) kết hợp mô hình hóa định lượng và đối sánh thực nghiệm hiện trường đa điểm. Thiết kế nghiên cứu đa tầng bao gồm:

  • Cấp độ vĩ mô (Macro-scale): Đánh giá phân tầng sơ bộ trên 7 phân lưu chính của sông Cửu Long (Cửa Tiểu, Cửa Đại, Hàm Luông, Cổ Chiên, Cung Hầu, Định An, Trần Đề) bằng giải tích chuỗi số liệu thủy văn và tham số thực nghiệm $\eta$.
  • Cấp độ vi mô (Micro-scale): Khảo sát đo đạc chuyên sâu và thiết lập mô hình toán 3 chiều chi tiết cho hệ thống sông Hậu (nhánh Định An và Trần Đề) từ trạm Đại Ngãi ra đến mép biển Đông ngoài đường đẳng sâu 20 m.

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu hiện trường được chuẩn hóa theo quy chuẩn quan trắc hải dương học và thủy văn quốc tế:

  1. Chiến dịch khảo sát hiện trường đồng bộ: Đo đạc cấu trúc độ mặn theo chiều sâu cột nước tại 6 vị trí trắc đoạn trọng điểm (T1 đến T6) dọc theo hai nhánh Định An và Trần Đề trong thời kỳ triều cường và triều kém mùa khô năm 2016 (từ 31/3/2016 đến 01/4/2016).
  2. Thiết bị và giao thức quan trắc: Sử dụng máy đo độ mặn tự ghi đa chỉ tiêu đa tầng kết hợp thiết bị đo dòng chảy biên dạng âm thanh ADCP (Acoustic Doppler Current Profiler) để đo đồng thời trường vận tốc dòng chảy 3 chiều và nồng độ mặn tại các thời điểm chuẩn: Nước ngưng đỉnh triều (High Water Slack - HWS) và nước ngưng chân triều (Low Water Slack - LWS).
  3. Quy trình tam giác giác hóa (Triangulation): Kiểm chứng chéo kết quả giữa (1) số liệu đo đạc thực tế tầng mặt - tầng đáy, (2) kết quả tính toán giải tích tham số $n$ và $\eta$, (3) kết quả mô phỏng số 3 chiều MIKE 3.

Data và phân tích

Luận án ứng dụng bộ công cụ phần mềm MIKE 3 Flow Model FM (Flexible Mesh) do DHI phát triển, sử dụng hệ tọa độ sigma chia cột nước thành nhiều lớp linh hoạt theo chiều sâu:

  • Dữ liệu địa hình và lưới tính: Lưới phi cấu trúc tam giác/tứ giác bao phủ toàn bộ vùng cửa sông Hậu và dải ven biển lân cận, độ phân giải lưới mịn tại khu vực lòng dẫn và thắt hẹp ($20 - 50$ m), thưa dần về phía biển khơi ($200 - 500$ m).
  • Điều kiện biên khí tượng - thủy văn:
    • Biên thượng lưu: Chuỗi lưu lượng thực đo theo thời gian tại trạm Châu Đốc và Cần Thơ.
    • Biên biển: Dao động mực nước điều hòa từ 9 sóng triều chính trích xuất từ mô hình toàn cầu kết hợp trạm hải văn Côn Đảo và Vũng Tàu; nồng độ mặn biên biển gán giá trị ổn định 32 - 34‰.
    • Biên khí quyển: Trường gió bề mặt 2 chiều (U, V) với bước thời gian 1 giờ được trích xuất từ dữ liệu tái phân tích khí hậu toàn cầu CFSR (Climate Forecast System Reanalysis) của NOAA.
  • Kiểm định và đánh giá độ tin cậy mô hình: Sử dụng hệ số Nash-Sutcliffe ($NSE$) và hệ số tương quan tuyến tính ($R^2$): $$NSE = 1 - \frac{\sum_{i=1}^N (O_i - P_i)^2}{\sum_{i=1}^N (O_i - \bar{O})^2}, \quad R^2 = \frac{\left[\sum_{i=1}^N (O_i - \bar{O})(P_i - \bar{P})\right]^2}{\sum_{i=1}^N (O_i - \bar{O})^2 \sum_{i=1}^N (P_i - \bar{P})^2}$$ Kết quả hiệu chỉnh và kiểm định tại các trạm kiểm tra (Mỹ Thanh, Đại Ngãi, T1, T2, T3, T4, T5, T6) đạt độ chính xác rất cao: $R^2$ dao động từ 0,86 đến 0,95 đối với mực nước và từ 0,81 đến 0,91 đối với độ mặn; $NSE$ đạt từ 0,78 đến 0,89, khẳng định độ tin cậy tuyệt đối của công cụ mô phỏng.

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

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

  1. Xác nhận sự tồn tại của nêm mặn và hiện tượng phân tầng mật độ: Kết quả đo đạc ngày 31/3/2016 và 01/4/2016 chứng minh rõ nét sự phân tầng trên cả hai nhánh sông Hậu. Tại nhánh Định An thời điểm chân triều (LWS), chỉ số Pritchard $n$ đạt các giá trị từ 0,42 đến 0,85; chênh lệch nồng độ mặn giữa tầng đáy và tầng mặt ($\Delta S = S_{bot} - S_{surf}$) dao động từ 4,0 g/l đến 11,5 g/l tại các vị trí cách cửa biển từ 10 km đến 25 km.
  2. Quy luật động lực nêm mặn theo chu kỳ triều:
    • Pha triều lên: Dòng triều mang nêm mặn có tỷ trọng lớn hơn luồn sâu dưới đáy sông, đẩy khối nước ngọt nổi lên trên và trôi ngược về phía thượng lưu. Độ dốc mặt phân cách nêm mặn tăng cao nhất tại thời điểm nước ngưng đỉnh triều (HWS).
    • Pha triều xuống: Khi gradien mực nước hướng ra biển gia tăng, dòng nước ngọt thượng nguồn đổ dồn với vận tốc lớn ở lớp mặt ($0,8 - 1,4$ m/s), trượt trên nêm mặn tầng đáy. Hiện tượng trượt gây mất ổn định ranh giới phân cách, hình thành các sóng xoáy cuộn rối K-H làm bề dày nêm mặn mỏng dần và bị đẩy lùi nhanh chóng ra phía biển.
  3. Sự bất đối xứng không gian giữa nhánh Định An và Trần Đề: Nhánh Định An có lòng dẫn rộng và sâu hơn nên nêm mặn xâm nhập xa hơn (vượt qua cồn Dung vào sâu hơn nhánh Trần Đề từ 3 đến 5 km trong cùng pha triều). Trên mặt cắt ngang, nêm mặn tập trung sâu nhất ở luồng lạch chính (thalweg), trong khi khu vực bãi nông ven bờ có độ mặn thấp hơn và mức độ xáo trộn đồng nhất diễn ra sớm hơn.
  4. Tác động phi tuyến của nước biển dâng (NBD): Kịch bản biến đổi khí hậu và NBD làm gia tăng chiều dài nêm mặn thêm 4,5 - 7,2 km về phía thượng lưu sông Hậu, nhưng đồng thời kéo dài thời gian duy trì cấu trúc phân tầng tại các hố sâu địa hình đáy tự nhiên.

Implications đa chiều

  • Về mặt lý thuyết học thuật: Chứng minh tính quy luật của dòng chảy mật độ hai pha trong các hệ thống sông nhiệt đới chịu triều cường; cung cấp cơ sở dữ liệu thực nghiệm quý giá cho mạng lưới nghiên cứu thủy văn quốc tế về lưu vực sông Mê Công.
  • Đổi mới phương pháp luận kỹ thuật: Khẳng định sự cần thiết phải chuyển dịch từ mô hình toán 1D sang mô hình 3D trong tính toán dự báo xâm nhập mặn phục vụ thiết kế công trình cửa sông. Mô hình 1D luôn ước tính nồng độ mặn mặt cắt cao hơn thực tế nồng độ mặn tầng mặt, dẫn đến việc đóng cống ngăn mặn quá sớm và lãng phí nguồn nước ngọt quý giá.
  • Ứng dụng công trình thực tiễn:
    • Ý tưởng cửa van cống phân tầng nhiều lớp / Cửa van mẹ bồng con: Thiết kế công trình cống vùng ven biển có nhiều tầng cửa van cho phép linh hoạt hạ cánh van đáy để chặn nêm mặn tỷ trọng cao xâm nhập, đồng thời mở cánh van mặt để lấy trọn vẹn lớp nước ngọt ($S < 1 - 2$ g/l) vào nội đồng trong pha triều xuống (Hình 4-5).
    • Ứng dụng trạm bơm nổi di động: Điển hình như mô hình trạm bơm nổi Xuân Hòa trên sông Tiền (Hình 4-6), định vị ống hút nước ở độ sâu $0,5 - 1,5$ m tính từ mặt nước để bơm cấp nước ngọt trong các khung giờ vàng phân tầng, kéo dài thời gian lấy nước thêm $2,5 - 4,0$ giờ mỗi ngày so với trạm bơm cố định đáy sâu.
  • Chính sách quản lý tài nguyên nước: Cung cấp cơ sở khoa học để Bộ Nông nghiệp & PTNT cùng các tỉnh Tây Nam Bộ ban hành quy trình vận hành liên hồ chứa, liên cống kiểm soát mặn theo thời gian thực dựa trên cấu trúc phân tầng độ mặn.

Limitations và Future Research

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

  1. Phạm vi không gian mô phỏng 3D: Do khối lượng tính toán của mô hình phi thủy tĩnh 3D rất lớn, nghiên cứu tập trung chi tiết vào hạ lưu sông Hậu (Định An và Trần Đề), chưa thiết lập mô hình 3D đồng thời cho toàn bộ 6 cửa sông còn lại của sông Tiền.
  2. Thời gian đo đạc khảo sát thực địa: Chuỗi số liệu đo mặt cắt đứng chi tiết tập trung chủ yếu vào đợt khô hạn khốc liệt năm 2016, chưa bao quát liên tục nhiều năm thủy văn ẩm và trung bình.
  3. Biến động hình thái lòng dẫn: Mô hình chưa tích hợp động lực bồi xói bùn cát và sự thay đổi hình thái lòng dẫn theo thời gian thực do hoạt động khai thác cát và nạo vét luồng hàng hải Định An.
  4. Tương tác sóng - dòng chảy tầng mặt: Ảnh hưởng của sóng biển vỡ ở cửa sông nông chưa được ghép nối hoàn toàn (fully coupled) với mô hình truyền chất 3D.

Chương trình nghiên cứu tiếp nối trong tương lai bao gồm:

  • Mở rộng thiết lập mô hình 3D toàn diện cho toàn bộ hệ thống các cửa sông Cửu Long và hệ thống sông Vàm Cỏ.
  • Nghiên cứu mô hình liên kết thủy lực - bùn cát - độ mặn 3D để đánh giá vùng lắng đọng phù sa cực đại (Estuarine Turbidity Maximum - ETM) tại mặt phân cách nêm mặn.
  • Xây dựng hệ thống cảm biến đo mặn đa tầng tự động truyền dữ liệu IoT theo thời gian thực kết hợp thuật toán trí tuệ nhân tạo (AI/Machine Learning) dự báo nêm mặn phục vụ điều khiển tự động hệ thống cống thông minh.

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

  • Tác động học thuật: Luận án là công trình tiên phong tại Việt Nam nghiên cứu toàn diện cơ chế phân tầng nêm mặn 3D tại ĐBSCL, mở ra hướng nghiên cứu thủy lực tài nguyên nước chuyên sâu với tiềm năng trích dẫn cao trong các tạp chí quốc tế thuộc danh mục ISI/Scopus (Q1/Q2) về cửa sông và thủy văn duyên hải.
  • Tác động chuyển đổi ngành thủy lợi: Cung cấp luận cứ khoa học trực tiếp thay đổi tiêu chuẩn thiết kế công trình cống và trạm bơm vùng ven biển Việt Nam từ phương pháp lấy nước đơn tầng truyền thống sang công nghệ lấy nước phân tầng đa chiều.
  • Tác động kinh tế - xã hội: Việc gia tăng khả năng lấy nước ngọt thêm $15 - 25%$ khối lượng trong các tháng mùa kiệt giúp bảo vệ hàng trăm nghìn hecta diện tích lúa Đông Xuân, vườn cây ăn trái đặc sản tại các tỉnh Bến Tre, Trà Vinh, Sóc Trăng, Tiền Giang, tiết kiệm hàng trăm tỷ đồng chi phí đắp đập tạm ngăn mặn hàng năm.

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

  • Nghiên cứu sinh và học giả chuyên ngành: Tiếp cận khung lý thuyết chuẩn hóa, hệ phương trình trạng thái mật độ và bộ tham số thực nghiệm $\eta, n$ làm nền tảng phát triển các đề tài nghiên cứu chuyên sâu về động lực học vùng biển ven bờ.
  • Kỹ sư thiết kế và tư vấn thủy lợi: Ứng dụng quy trình mô phỏng 3D MIKE 3 và giải pháp công nghệ cửa van mẹ bồng con, trạm bơm nổi hớt ngọt vào các dự án công trình kiểm soát nguồn nước thực tế.
  • Cơ quan quản lý nhà nước (Bộ NN&PTNT, Viện Quy hoạch Thủy lợi, Chi cục Thủy lợi các tỉnh): Sử dụng các biểu đồ phân tầng và ngưỡng xâm nhập nêm mặn để lập lịch vận hành công trình ngăn mặn, trữ ngọt chính xác và linh hoạt.
  • Cộng đồng nông dân và hợp tác xã nông nghiệp ven biển: Hưởng lợi trực tiếp từ nguồn nước ngọt bổ sung ổn định trong mùa kiệt, giảm thiểu rủi ro thiệt hại mùa màng do thiên tai hạn mặn.

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ì? Đó là việc chứng minh và định lượng hóa thành công sự tồn tại của cấu trúc nêm mặn phân tầng hai lớp trên các nhánh sông hạ lưu ĐBSCL trong mùa kiệt. Luận án đã mở rộng lý thuyết phân loại của D. Pritchard (1955) và chuẩn hóa tham số Simmons cải tiến $\eta = W^+/W^-$ của TSKH. Nguyễn Ân Niên, xác lập ranh giới định lượng giữa các trạng thái phân tầng động lực học cửa sông nhiệt đới.

  2. Điểm mới về phương pháp luận nghiên cứu so với các công trình trước đây? Trước đây, các nghiên cứu trong nước chủ yếu áp dụng mô hình 1D (như MIKE 11, ISIS) chỉ tính được nồng độ mặn trung bình mặt cắt ngang. Luận án đã đột phá ứng dụng mô hình phi thủy tĩnh 3 chiều MIKE 3 FM kết hợp lưới phi cấu trúc và tọa độ sigma, hiệu chỉnh bằng chuỗi dữ liệu thực đo phân tầng đa điểm theo thời gian thực (HWS, LWS) với các chỉ số độ chính xác cao ($R^2 > 0,85$; $NSE > 0,78$).

  3. Phát hiện nào bất ngờ nhất có số liệu thực nghiệm chứng minh? Đó là hiện tượng nồng độ mặn tầng mặt tại một số thời điểm chân triều ở nhánh Định An và Trần Đề vẫn duy trì ở ngưỡng an toàn cho nông nghiệp ($S < 2 - 4$ g/l) trong khi nồng độ mặn tầng đáy đã vượt ngưỡng $10 - 15$ g/l. Số liệu đo đạc ngày 31/3/2016 tại trạm cách biển 15 km ghi nhận chênh lệch mặn thẳng đứng $\Delta S$ lên tới 11,5 g/l, tạo cơ hội lấy nước ngọt tầng mặt ngay trong giai đoạn đỉnh điểm của đợt hạn mặn lịch sử.

  4. Luận án có cung cấp quy trình tái lập nghiên cứu (replication protocol) không? Có. Luận án cung cấp đầy đủ và chi tiết toàn bộ thông số thiết lập mô hình (Bảng 3-1), hệ số nhám đáy Manning, hệ số khuếch tán rối ngang/đứng, điều kiện biên dòng chảy Châu Đốc - Cần Thơ, biên sóng triều hải văn, trường gió CFSR và tọa độ chính xác các mặt cắt trích xuất kết quả (Hình 3-1, 3-3, 3-31), cho phép các nhà khoa học độc lập tái lập và kiểm chứng mô hình.

  5. Chương trình nghiên cứu 10 năm tới định hình ra sao? Tập trung tích hợp mô hình thủy động lực học nêm mặn 3D vào hệ thống giám sát tự động viễn thám kết hợp mạng lưới quan trắc thông minh IoT, xây dựng nền tảng bản sao số (Digital Twin) lưu vực sông Cửu Long nhằm tự động hóa điều khiển thời gian thực các cụm công trình thủy lợi liên vùng thích ứng với biến đổi khí hậu toàn cầu.

Kết luận

  1. Luận án đã làm sáng tỏ một cách toàn diện và khoa học cơ chế xâm nhập mặn, sự hình thành, cấu trúc không gian và quy luật vận động của nêm mặn phân tầng tại vùng cửa sông ven biển ĐBSCL.
  2. Ứng dụng thành công hệ thống tham số kinh nghiệm ($\eta, n$) để định lượng và phân loại mức độ phân tầng cho 7 nhánh sông Cửu Long chính, xác nhận tính chất phân tầng một phần và phân tầng mạnh theo chu kỳ triều.
  3. Thiết lập thành công mô hình toán 3 chiều MIKE 3 FM mô phỏng chính xác trường độ mặn 3 chiều và động lực học nêm mặn tại hạ lưu sông Hậu với độ tin cậy được kiểm chứng qua hệ số $NSE > 0,78$ và $R^2 > 0,85$.
  4. Đề xuất hệ thống giải pháp công trình và phi công trình đột phá: giải pháp khai thác nước ngọt phân tầng theo không gian và thời gian pha triều, ý tưởng thiết kế cửa van nhiều lớp (cửa van mẹ bồng con) và ứng dụng trạm bơm nổi cơ động.
  5. Tạo bước chuyển biến mang tính hệ hình (paradigm shift) từ tư duy kiểm soát mặn 1D thụ động sang tư duy khai thác tài nguyên nước ngọt 3D chủ động thích ứng, mở ra hướng đi chiến lược cho sự phát triển bền vững của Đồng bằng sông Cửu Long trước thách thức kép của biến đổi khí hậu và tác động thượng nguồn Mê Công.