Tổng quan về luận án

Nghiên cứu về lan truyền nước chua phèn tại vùng đồng bằng châu thổ sông Mê Kông là một trong những bài toán phức tạp và thách thức nhất của ngành Kỹ thuật Tài nguyên nước (Mã ngành: 9 58 02 12). Tại Đồng bằng sông Cửu Long (ĐBSCL), đất phèn (Acid Sulfate Soils - ASS) chiếm khoảng 1,5 triệu ha (tương đương 45% tổng diện tích tự nhiên), phân bố tập trung tại ba vùng trũng lớn: Đồng Tháp Mười (ĐTM), Bán đảo Cà Mau (BĐCM) và Tứ giác Long Xuyên (TGLX). Bản chất của đất phèn là quá trình hình thành axit sunfuric ($H_2SO_4$) tác động sâu sắc đến các đặc tính lý hóa của môi trường. Trong mùa khô, hiện tượng mao dẫn đưa các hợp chất chứa độc chất sinh phèn ($Fe^{2+}, Fe^{3+}, Al^{3+}, SO_4^{2-}, H^+$) lên tầng đất mặt; khi các trận mưa rào đầu mùa xuất hiện, quá trình thủy phân diễn ra mạnh mẽ, tạo ra các dòng chảy tràn mặt ruộng mang nồng độ axit cực cao đổ dồn vào mạng lưới kênh rạch, hình thành các "rốn phèn" làm suy thoái nghiêm trọng nguồn nước mặt, hủy hoại hệ sinh thái thủy sinh và đe dọa trực tiếp đến an ninh sản xuất nông - thủy sản.

Khoảng trống nghiên cứu (Research gap) then chốt mà luận án giải quyết xuất phát từ hạn chế của các mô hình toán truyền tải chất truyền thống: hầu hết các phần mềm thủy lực và chất lượng nước thương mại hiện nay (như MIKE 11 của DHI, SOBEK, QUAL2K) chủ yếu giải quyết bài toán lan truyền của các chất bảo toàn (conservative substances như độ mặn $NaCl$), hoặc chỉ mô phỏng các chất phân rã bậc nhất đơn giản. Bài toán lan truyền nước phèn là bài toán vận chuyển chất phi bảo toàn (non-conservative reactive transport) chịu sự chi phối đồng thời của động lực học dòng chảy một chiều vùng triều và các quy luật cân bằng hóa sinh phức tạp giữa các ion kim loại nặng và gốc axit. Luận án đã xác lập một hệ thống câu hỏi nghiên cứu và giả thuyết khoa học tường minh:

  • RQ1: Bản chất hóa sinh và các phương trình cân bằng chi phối sự hòa tan của phèn nhôm và phèn sắt tại vùng TGLX diễn ra theo quy luật nhiệt động học nào?
  • RQ2: Làm thế nào để đóng kín hệ phương trình vi phân mô tả quá trình lan truyền chất phi bảo toàn trong mạng lưới kênh hở chịu ảnh hưởng của triều và loại trừ hiện tượng khuếch tán số (numerical diffusion)?
  • RQ3: Mối tương quan động học giữa lượng mưa đầu mùa, quá trình rửa trôi trên ô ruộng và diễn biến không - thời gian của các trục tiêu phèn ra biển Tây diễn ra như thế nào?
  • H1: Sự hòa tan và thủy phân của phèn nhôm và phèn sắt trong môi trường nước TGLX tuân thủ nghiêm ngặt cân bằng nhiệt động Jurbanite ($AlOHSO_4$) và Jarosite ($KFe_3(SO_4)_2(OH)_6$).
  • H2: Việc tích hợp phương pháp phân rã toán tử (split-operator method) với phương pháp đường đặc trưng (method of characteristics) sẽ triệt tiêu sai số khuếch tán số trong bài toán hyperbolic bậc một.
  • H3: Việc vận hành hệ thống 30 cống ven biển Tây theo chu kỳ triều kết hợp giải pháp ém phèn mặt ruộng có thể kiểm soát hiệu quả sự lan truyền phèn ngay trong 30-45 ngày đầu mùa mưa.

Khung lý thuyết của công trình được xây dựng trên sự giao thoa giữa Thủy động lực học công trình thủy lợi (Saint-Venant 1D), Lý thuyết truyền tải - phân tán Advection-Dispersion và Hóa học đất ngập nước. Đóng góp đột phá của luận án là việc xây dựng thành công chương trình tính toán chuyên dụng ACID2020 (lập trình trên Digital Visual Fortran 64-bit), tích hợp mô-đun hóa sinh phi bảo toàn với lõi thủy lực mô hình DELTA (GS.TS Nguyễn Tất Đắc). Quy mô nghiên cứu bao phủ toàn bộ vùng TGLX với sơ đồ tính toán chi tiết gồm 1.061 mặt cắt tính toán, 424 nhánh sông kênh, 233 nút hợp lưu, 5 biên thủy lực chính (Châu Đốc, Vàm Nao, Cần Thơ, Rạch Giá, Hà Tiên) và tập trung phân tích sâu tiểu vùng Tứ giác Hà Tiên (kênh Hà Giang) với bước thời gian tính toán $\Delta t = 900$ giây trên chuỗi số liệu thực nghiệm năm 2016.

       +--------------------------------------------------------------------+
       |          HỆ THỐNG MẶT RUỘNG TRŨNG PHÈN TỨ GIÁC LONG XUYÊN          |
       |  Mùa khô: Mao dẫn ion Al3+, Fe3+, SO42- tích tụ tầng mặt ruộng     |
       +---------------------------------+----------------------------------+
                                         |
                                         v (Mưa đầu mùa rửa trôi & thủy phân)
       +--------------------------------------------------------------------+
       |              MÔ HÌNH Ô RUỘNG (Linked Field-Plot Model)             |
       |         V dC/dt = Qr Cr - Qs Cs - Qe Ce + F(Cg); Cg = Coi e^-a0t   |
       +---------------------------------+----------------------------------+
                                         |
                                         v (Dòng chảy tràn mang nồng độ pH ~ 2.0)
       +--------------------------------------------------------------------+
       |             MÔ HÌNH LAN TRUYỀN PHÈN SÔNG KÊNH (ACID2020)           |
       |  dS/dt + U dS/dx = E d2S/dx2 + Reaction Terms (Jurbanite/Jarosite) |
       |   + Lõi thủy lực: DELTA (1D Saint-Venant, 1061 mặt cắt, 424 nhánh) |
       |   + Giải phương trình tải: Phương pháp đường đặc trưng             |
       |   + Cân bằng nút hợp lưu: Tổng dòng ion vào = Tổng dòng ion ra    |
       +---------------------------------+----------------------------------+
                                         |
                                         v
       +--------------------------------------------------------------------+
       |      ĐIỀU TIẾT HỆ THỐNG CÔNG TRÌNH & TIÊU THOÁT BIỂN TÂY           |
       | 30 cống ven biển Tây (Đầm Chích, Hà Giang, Thần Nông) + Đập Tha La |
       +--------------------------------------------------------------------+

Literature Review và Positioning

Lịch sử nghiên cứu đất phèn thế giới ghi nhận từ năm 1735, Peelman đã phát hiện loại khoáng biến thành đất chua mang tên Agrilla Vitrolacea trên nền đá trầm tích kỷ Cambri vùng ven biển. Đến năm 1937, Doyne xác lập luận điểm nền tảng cho rằng nguồn gốc lưu huỳnh ($S$) trong đất phèn bắt nguồn từ nước biển và thảm thực vật sú vẹt ngập mặn (Rhizophora) chính là vật liệu hữu cơ tiền đề tạo nên pyrite ($FeS_2$). Tiến trình học thuật quốc tế về đất phèn được hệ thống hóa qua Hội nghị Đất phèn Quốc tế lần thứ nhất tại Đại học Wageningen (Hà Lan, 1972) nhằm chuẩn hóa các phương pháp phân loại (theo hệ thống USDA Soil Taxonomy và FAO-UNESCO). Đến Hội nghị Đất phèn Thế giới lần thứ tư tổ chức tại TP. Hồ Chí Minh (1992), cộng đồng khoa học đã ghi nhận bước tiến nhảy vọt về thiết bị đo đạc và phương pháp phân tích hóa sinh môi trường đất phèn nhiệt đới.

Trong nước, GS. Lê Huy Bá (2003) đã phân tích sâu sắc độc chất học môi trường đất phèn, khẳng định độc tính của các ion $Al^{3+}, Fe^{2+}, Fe^{3+}, SO_4^{2-}$ biến động phức tạp và gây tác hại gần như không thể đảo ngược đối với hệ sinh thái nông nghiệp và thủy sản. Về mặt hóa học nhiệt động, Nguyễn Thành Tín (1990) đã chứng minh cân bằng Jurbanite cho phèn nhôm tại Tân Thạnh (Đồng Tháp Mười) tuân thủ mối quan hệ tuyến tính giữa độ chua và hàm lượng sunfat:

"cân bằng Jurbanite cho phèn nhôm ở vùng Tân thạnh, Đồng Tháp Mười tuân thủ theo cân bằng: $pH = pSu + pAl + d$"

Đồng thời, Đặng Trung Thuận (2005) chỉ ra rằng trong môi trường TGLX, quá trình phong hóa khoáng sunfua tạo ra sản phẩm cuối cùng là oxit sắt không tan:

"khi pyrite bị phong hóa thì sản phẩm cuối cùng là oxit sắt, hợp chất này khó hòa tan và tích đọng dưới dạng mầu nâu"

Trong bức tranh tổng quan y văn, tồn tại hai luồng quan điểm học thuật đối lập về cơ chế kiểm soát chất lượng nước phèn:

  1. Trường phái Cân bằng Hóa sinh Cục bộ (Equilibrium-driven school): Đại diện bởi các nghiên cứu của Van Breemen (1973, 1976) và Dent (1986) tại Thái Lan và Hà Lan, cho rằng nồng độ các ion độc trong nước hoàn toàn bị chi phối bởi trạng thái cân bằng nhiệt động pha rắn - lỏng tại bề mặt tiếp xúc giữa đất phẫu diện và dung dịch đất.
  2. Trường phái Động lực học Thủy văn (Hydrodynamics-dominated school): Đại diện bởi các nghiên cứu của White et al. (1997) và Johnston et al. (2004, 2005) tại các vùng đất ngập nước trũng ven biển New South Wales (Úc), nhấn mạnh rằng sự phát tán axit trong hệ thống kênh rạch phụ thuộc chủ yếu vào cường độ mưa rào cực đoan, động lực tiêu thoát dòng tràn và chế độ truyền sóng triều hai chiều.

Luận án của Trần Ký định vị chính xác tại điểm giao thoa của hai trường phái trên: tích hợp mô hình phát thải nguồn phân tán từ ô ruộng (phản ứng hòa tan bề mặt) vào mạng lưới thủy lực động 1D phân tán (advection-dispersion). So sánh với các nghiên cứu quốc tế, mô hình của Johnston et al. (2005) tại Úc chủ yếu sử dụng mô hình hộp đơn giản (lumped box model) không phản ánh được hiện tượng đổi chiều dòng chảy do bán nhật triều không đều, trong khi mô hình ACID2020 của tác giả giải quyết tường minh sự lan truyền của sóng axít trên mạng lưới sông kênh đa nhánh với 233 nút hợp lưu phức tạp.

Đó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 lý thuyết truyền chất của Taylor (1954) và Fischer (1979) đối với các chất phản ứng không bảo toàn (reactive non-conservative solutes). Đóng góp lý thuyết cốt lõi là việc thiết lập các phương trình toán học bổ sung để đóng kín hệ phương trình vi phân chuyển tải:

  1. Đối với phèn nhôm (Cân bằng Jurbanite): $$\log[Al^{3+}] + \log[OH^-] + \log[SO_4^{2-}] = K_{sp}$$ Chuyển đổi sang phương trình quan hệ trạng thái: $$pH = pSu + pAl + d$$ Biến đổi vi phân theo không gian và thời gian: $$\frac{\partial H}{\partial t} = \frac{H}{Al} \frac{\partial Al}{\partial t} + \frac{H}{Su} \frac{\partial Su}{\partial t}$$

  2. Đối với phèn sắt (Cân bằng Jarosite và Oxy hóa Pyrite): $$FeS_2 + \frac{7}{2} O_2 + H_2O \rightarrow Fe^{2+} + 2SO_4^{2-} + 2H^+$$ Phương trình xác định pH theo nồng độ sắt và sunfat tức thời: $$pH = m \cdot pSu + \alpha \cdot pFe + \beta$$ $$\frac{\partial H}{\partial t} = \alpha \frac{H}{Fe} \frac{\partial Fe}{\partial t} + m \frac{H}{Su} \frac{\partial Su}{\partial t}$$

Sự chuyển dịch hệ hình (Paradigm shift) ở đây là việc chuyển từ việc coi pH như một đại lượng đo độc lập sang việc coi ion $H^+$ là một biến phụ thuộc động học, được liên kết trực tiếp với các đạo hàm riêng của $Al^{3+}, Fe^{3+}$ và $SO_4^{2-}$ thông qua các quy luật bảo toàn khối lượng và cân bằng hóa trị.

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

Khung phân tích của luận án tích hợp ba thành phần toán học chặt chẽ:

+-----------------------------------------------------------------------------------+
|                            KHUNG PHÂN TÍCH TOÁN HỌC ACID2020                      |
+-----------------------------------------------------------------------------------+
| 1. KHỐI Ô RUỘNG (Field Plot Dynamics):                                            |
|    - Cân bằng nước:  dV/dt = Qr - Qe - Qs                                         |
|    - Cân bằng chất:  d(VC)/dt = Qr Cr - Qe Ce - Qs Cs + F(Cg)                     |
|    - Hòa tan phèn:   Cg = Coi * exp(-a0 * t)                                      |
+-----------------------------------------------------------------------------------+
                                      |
                                      v (Trao đổi lưu lượng Qs & nồng độ Cs)
+-----------------------------------------------------------------------------------+
| 2. KHỐI KÊNH SÔNG MẠNG LƯỚI (1D Canal Network):                                   |
|    - Thủy lực dòng chảy: delta-system (Saint-Venant 1D)                           |
|    - Truyền tải chất:    dS/dt + U * dS/dx = E * d2S/dx2 + R(Si)                  |
|    - Tách toán tử:       (i) Tải thuần túy: dSt/dt + U * dSt/dx = 0 (Đặc trưng)   |
|                          (ii) Hiệu chỉnh khuếch tán và phản ứng hóa sinh          |
+-----------------------------------------------------------------------------------+
                                      |
                                      v
+-----------------------------------------------------------------------------------+
| 3. BẢO TOÀN NÚT HỢP LƯU (Confluence Boundary Condition):                          |
|    - Nồng độ nút: SN = (Tổng Qin * Sin) / (Tổng Qout)                             |
+-----------------------------------------------------------------------------------+

Điều kiện biên (Boundary conditions): Tại các biên sông lớn (Châu Đốc, Vàm Nao, Cần Thơ) và biên biển (Rạch Giá, Hà Tiên), nguồn nước được xem là không nhiễm phèn ($pH \ge 6.0$). Biên phát thải ô ruộng được thiết lập với diện tích chuẩn 10 ha tại vùng trũng kênh Hà Giang với nồng độ ban đầu cực đại tương ứng $pH = 2.0$.

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

Thiết kế nghiên cứu

Nghiên cứu đứng trên lập trường bản thể luận khách quan và nhận thức luận Chủ nghĩa thực chứng hậu kỳ (Post-positivism) kết hợp với phương pháp luận Cơ chế luận toán lý (Mechanistic physical-chemical modeling). Thiết kế nghiên cứu tuân thủ quy trình mô phỏng đa cấp độ (Multi-level design):

  • Cấp vi mô: Động học giải phóng chất chua trong ô ruộng cục bộ diện tích 10 ha.
  • Cấp trung mô: Động học lan truyền và pha loãng trên trục kênh cấp 1 và cấp 2 (kênh Hà Giang, kênh Trà Phô, kênh Quốc Phòng, kênh Nông Trường).
  • Cấp vĩ mô: Động lực học phân phối dòng axit trên toàn mạng lưới TGLX ra vùng biển Tây.
       +-------------------------------------------------------------------+
       |            QUY TRÌNH MÔ HÌNH HÓA VÀ GIẢI SỐ ACID2020              |
       +---------------------------------+---------------------------------+
                                         |
                                         v
       +-------------------------------------------------------------------+
       | Thu thập & Chuẩn hóa dữ liệu:                                     |
       | - Địa hình 2011 (1061 mặt cắt, 424 nhánh, 233 nút hợp lưu)        |
       | - Thủy văn 2016 (5 trạm biên mực nước, 4 trạm mưa)                |
       | - Hóa học: Dung dịch đất, phẫu diện đất, nồng độ ion Al, Fe, SO4  |
       +---------------------------------+---------------------------------+
                                         |
                                         v
       +-------------------------------------------------------------------+
       | Lõi Thủy lực DELTA (GS.TS Nguyễn Tất Đắc):                        |
       | - Giải hệ phương trình Saint-Venant 1 chiều                       |
       | - Xuất trường vận tốc U(x, t) và mực nước Z(x, t)                 |
       +---------------------------------+---------------------------------+
                                         |
                                         v
       +-------------------------------------------------------------------+
       | Thuật toán Phân rã Toán tử (Split-Operator Method):               |
       | + Bước 1: Giải truyền tải thuần túy bằng Phương pháp Đường đặc    |
       |           trưng (Method of Characteristics) -> Triệt tiêu sai số  |
       |           khuếch tán số.                                          |
       | + Bước 2: Hiệu chỉnh khuếch tán vật lý và phản ứng hóa học        |
       |           (Cân bằng Jurbanite/Jarosite).                          |
       +---------------------------------+---------------------------------+
                                         |
                                         v
       +-------------------------------------------------------------------+
       | Khối ô ruộng (Field Plot Model):                                  |
       | - Mô phỏng cơ chế rửa trôi phèn bề mặt: Cg = Coi * e^(-a0*t)      |
       | - Trao đổi hai chiều qua cống/mặt tràn với kênh trục              |
       +---------------------------------+---------------------------------+
                                         |
                                         v
       +-------------------------------------------------------------------+
       | Kiểm tra Điều kiện Ổn định & Biên dịch:                           |
       | - Bước thời gian dt = 900s; Thỏa mãn điều kiện Courant (CFL <= 1) |
       | - Biên dịch bằng Digital Visual Fortran 64-bit qua tệp Fvao.txt   |
       +-------------------------------------------------------------------+

Quy trình nghiên cứu rigorous

Luận án xử lý xuất sắc bài toán sai số số học. Khi giải phương trình truyền tải advection dạng hyperbolic: $$\frac{\partial S_t}{\partial t} + U \frac{\partial S_t}{\partial x} = 0$$ Các phương pháp sai phân hữu hạn thông thường (như Upwind, Lax-Wendroff) tạo ra hiện tượng khuếch tán số (numerical diffusion) nghiêm trọng, làm sai lệch biên độ đỉnh phèn và làm nhẵn nhân tạo độ dốc nồng độ. Tác giả đã áp dụng phương pháp đường đặc trưng (method of characteristics) để truy vết nghiệm chính xác dọc theo quỹ đạo dòng chảy $dx/dt = U$, sau đó mới thực hiện bước hiệu chỉnh khuếch tán phân tán và cân bằng hóa học.

Về quy chuẩn dữ liệu, tác giả thực hiện kiểm chứng chéo (Data triangulation) từ ba nguồn độc lập: (1) Số liệu dự án "Cải tạo đất chua phèn ĐTM" của Ủy ban sông Mê Kông; (2) Dữ liệu phẫu diện đất từ Phân viện Quy hoạch và Thiết kế Nông nghiệp; (3) Dữ liệu quan trắc chất lượng nước định kỳ của Viện Khoa học Thủy lợi Miền Nam và Viện Quy hoạch Thủy lợi Miền Nam.

Data và phân tích

  • Đặc trưng mạng lưới: 1.061 mặt cắt hình học, 424 đoạn kênh, 233 nút giao, tích hợp 2 đập ngăn lũ vùng biên (Tha La, Trà Sư) và 23 cống điều tiết triều ngăn mặn/xả phèn ven biển Kiên Giang.
  • Tham số khí tượng - thủy văn: Chuỗi dữ liệu đo mưa ngày tại Châu Đốc, Long Xuyên, Rạch Giá, Hà Tiên. Tổng lượng mưa vùng TGLX đạt từ $2.100 \div 2.200$ mm/năm, trong đó $90 \div 93%$ tập trung vào mùa mưa (từ tháng 5 đến tháng 12). Lượng mưa năm lớn nhất ghi nhận tại Rạch Giá (2.136 mm) và Hà Tiên (1.995 mm) với số ngày mưa trung bình lần lượt là 159 ngày và 131 ngày.
  • Thông số mô hình hóa: Bước thời gian tính toán $\Delta t = 900$ s (15 phút), bảo đảm hệ số Courant-Friedrichs-Lewy $CFL = U \frac{\Delta t}{\Delta x} \le 1$.
  • Tiêu chuẩn kiểm chứng: Đánh giá nồng độ sunfat ($SO_4^{2-}$) và tổng sắt ($Fe_{tot}$) theo Quy chuẩn Kỹ thuật Quốc gia về chất lượng nước tưới tiêu QCVN 39:2011/BTNMT (ngưỡng giới hạn sunfat là 600 mg/L).

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

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

+-----------------------------------------------------------------------------------+
|             BẢNG ĐỐI SOÁT ĐỘNG HỌC PHÈN THEO THỜI GIAN VÀ KHÔNG GIAN             |
+----------------------+-----------------------------+------------------------------+
| THÔNG SỐ QUAN TRẮC   | GIAI ĐOẠN ĐẦU MÙA (THÁNG 5) | GIAI ĐOẠN GIỮA MÙA (THÁNG 6) |
+----------------------+-----------------------------+------------------------------+
| pH Kênh Hà Giang     | Cực kỳ chua: 2.03 - 2.87    | Tăng mạnh: 6.12 - 6.15       |
| (Mặt cắt 28, 29A)    | (Nước màu đỏ gỉ/xanh đen)   | (Pha loãng do mưa & lũ về)   |
+----------------------+-----------------------------+------------------------------+
| pH Kênh Nông Trường  | Axit trung bình: 4.03       | Đạt ngưỡng ngọt hóa: 6.75    |
| (Mặt cắt 34)         | (Xa nguồn phát thải phèn)   | (Tiếp nhận nước Vĩnh Tế)     |
+----------------------+-----------------------------+------------------------------+
| Hàm lượng Sunfat     | Biến động: 0 - 600 mg/L     | Cơ bản dưới 600 mg/L         |
| (SO4 2-)             | (29A, 30A vượt chuẩn QCVN)  | (Đạt chuẩn QCVN 39:2011)     |
+----------------------+-----------------------------+------------------------------+
| Hàm lượng Sắt tổng   | 0.25 - 1.35 mg/L            | Giảm mạnh do kết tủa         |
| (Fe tot)             | (Đỉnh cao tại Nút 28 & 32)  | với hạt phù sa sông Hậu      |
+----------------------+-----------------------------+------------------------------+
| Động lực dòng chảy   | Mực nước kênh < Mặt ruộng   | Mực nước kênh dâng cao       |
|                      | (Dòng tràn tự do ra kênh)   | (Dòng triều đẩy ra biển Tây) |
+----------------------+-----------------------------+------------------------------+
  1. Cơ chế bùng phát rốn phèn đầu tháng 5: Sau các trận mưa rào đầu mùa, toàn bộ lượng phèn tích tụ do mao dẫn trong mùa khô bị thủy phân. Tại khu vực Bắc kênh Hà Giang, đường đồng mức $pH = 2.03$ bao trùm phần lớn diện tích. Tại mặt cắt 28, pH đạt 2.87; tại mặt cắt 29A và 30A (thượng và hạ lưu cống Hà Giang), pH tụt giảm xuống 2.62. Nước kênh có màu đỏ gỉ sắt hoặc xanh đen đặc trưng.
  2. Quy luật lan truyền không gian hướng biển Tây: Dòng nước chua từ kênh Hà Giang nhanh chóng xâm nhập vào các trục kênh giao cắt gồm kênh Trà Phô, kênh Quốc Phòng và hệ thống kênh nội đồng, tạo thành trục dịch chuyển axít hướng về kênh Rạch Giá - Hà Tiên trước khi xả ra vịnh Thái Lan qua các cửa cống ven biển Tây.
  3. Hiệu ứng pha loãng động lực trong tháng 6: Khi lượng mưa tích lũy gia tăng mạnh mẽ, mực nước trong hệ thống kênh dâng cao, nồng độ ion tự do được pha loãng nhanh chóng. Giá trị pH tại mặt cắt 28 và 29A tăng vọt lên mức $6.12 \div 6.15$. Tại mặt cắt 34 (kênh Nông Trường), do nằm xa nguồn phát thải và nhận nguồn nước ngọt bổ cập từ kênh Vĩnh Tế và sông Hậu, giá trị pH đạt tới 6.75.
  4. Hiện tượng tương tác mặn - phèn tại cống ven biển: Hàm lượng sunfat ($SO_4^{2-}$) trên phần lớn mạng lưới đều nằm trong giới hạn an toàn (< 600 mg/L theo QCVN 39:2011/BTNMT). Tuy nhiên, tại các mặt cắt 29A và 30A giáp biển, hàm lượng sunfat vượt ngưỡng chuẩn nhiều lần, chứng minh sự tồn tại của hiện tượng ô nhiễm kép: vừa nhiễm phèn do dòng tràn nội đồng, vừa nhiễm mặn do dòng triều biển Tây xâm nhập vào thời điểm chân triều.
  5. Động học lắng đọng sắt phù sa: Hàm lượng tổng sắt ($Fe$) dao động trong khoảng $0.25 \div 1.35$ mg/L, đạt cực đại tại các nút giao 28 và 32 trong tháng 5. Bước sang tháng 6, khi dòng nước ngọt mang phù sa từ sông Hậu tràn qua kênh Vĩnh Tế vào vùng TGLX, các ion sắt bị trung hòa và kết tủa cùng các hạt cặn lơ lửng, làm nước chuyển từ màu đỏ gỉ sắt sang màu xám rêu lắng cặn.

Implications đa chiều

  • Về mặt Lý thuyết: Luận án đã chứng minh thành công tính khả thi của việc liên kết phương trình cân bằng hóa sinh Jurbanite/Jarosite với hệ phương trình truyền tải thủy lực 1D trên mạng lưới kênh hở phức tạp.
  • Về mặt Phương pháp luận: Phần mềm ACID2020 cung cấp một công cụ mã nguồn mở, tối ưu hóa bộ nhớ trên nền Fortran 64-bit, cho phép các cơ quan nghiên cứu trong nước tự chủ hoàn toàn về công nghệ mô phỏng chất lượng nước mà không phụ thuộc vào các phần mềm thương mại đắt đỏ.
  • Về mặt Thực tiễn và Quản lý Thủy nông:
    • Đề xuất quy trình vận hành một chiều cho cụm 30 cống ven biển Tây (đặc biệt là cống Đầm Chích, cống Hà Giang, cống Thần Nông): Mở cống xả phèn triệt để khi triều rút trong tháng 5 và đóng cống ngăn mặn, giữ ngọt khi lũ sông Hậu bắt đầu về trong tháng 6-7.
    • Ứng dụng giải pháp "ém phèn": Duy trì một lớp nước ngập thường xuyên trên mặt ruộng trũng trong thời kỳ chuyển mùa khô - mưa nhằm hòa tan tại chỗ độc chất và tạo dòng thấm đứng đẩy các ion $Al^{3+}, Fe^{3+}$ xuống tầng nước ngầm, ngăn chặn sự bốc phèn mao dẫn.

Limitations và Future Research

Luận án thẳng thắn thừa nhận các giới hạn học thuật nội tại:

  1. Giả thiết hòa đều lý tưởng trong ô ruộng: Luận án đưa ra giả thiết:

"với giả thiết nước phèn trong từng ô ruộng được hòa đều thì các phương trình cân bằng trong từng ô ruộng trở nên đơn giản" Giả thiết này xem ô ruộng như một bể phản ứng khuấy trộn hoàn toàn (CSTR), chưa phản ánh được sự biến thiên nồng độ phèn theo phương thẳng đứng trong tầng đất mặt và gradient nồng độ từ giữa ruộng ra bờ kênh.

  1. Công thức thực nghiệm nguồn phát thải: Nồng độ sinh phèn $C_g$ được tính toán qua hàm mũ thực nghiệm $C_g = C_{oi} e^{-a_0 t}$ với hệ số $a_0$ phụ thuộc vào số liệu quan trắc, chưa mô phỏng trực tiếp quá trình phản ứng hóa sinh trong cột đất (Soil column reaction kinetics).
  2. Giới hạn tham số hóa oxy hóa khử: Mô hình chưa tích hợp phương trình cân bằng oxy hòa tan ($DO$) và thế oxy hóa khử ($Eh$), vốn là hai yếu tố chi phối tốc độ chuyển hóa giữa $Fe^{2+}$ và $Fe^{3+}$.
  3. Cập nhật dữ liệu địa hình: Cơ sở dữ liệu hình học mặt cắt kênh rạch chủ yếu được đo đạc năm 2011, cần được cập nhật bằng công nghệ quét LiDAR và đo sâu hồi âm đa chùm hiện đại.

Chương trình nghiên cứu tương lai (Future research agenda):

  • Xây dựng mô hình ghép nối hoàn chỉnh 3 chiều (Coupled 3D Soil-Porous Media and 1D Canal Network Model).
  • Xây dựng bản đồ số GIS-Atlas về các hằng số nhiệt động $a_0, d, m, \alpha, \beta$ cho toàn bộ các tiểu vùng sinh thái đất phèn ĐBSCL.
  • Phát triển hệ thống hỗ trợ ra quyết định thời gian thực (Real-time Decision Support System - DSS) tích hợp trí tuệ nhân tạo để điều khiển tự động 30 cống ven biển Tây.

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

  • Tác động Học thuật: Công trình là tài liệu tham khảo cốt lõi cho chương trình đào tạo sau đại học ngành Kỹ thuật Tài nguyên nước tại Viện Khoa học Thủy lợi Việt Nam, Viện KHTL Miền Nam và Trường Đại học Tài nguyên và Môi trường TP.HCM. Cung cấp nền tảng thuật toán cho các nghiên cứu mô phỏng ô nhiễm chất không bảo toàn.
  • Chuyển đổi Nông nghiệp & Thủy sản: Cung cấp cơ sở khoa học để tái cơ cấu lịch thời vụ lúa vụ Hè Thu tại vùng Tứ giác Hà Tiên; giúp nông dân xác định chính xác thời điểm "nước hết chua" (tháng 6) để tiến hành xuống giống, loại trừ nguy cơ mất trắng do nghẹt rễ lúa và sốc phèn thủy sản.
  • Hoạch định Chính sách & Quy hoạch: Đóng góp luận cứ khoa học định lượng phục vụ Quy hoạch vùng ĐBSCL thời kỳ 2021-2030, tầm nhìn đến năm 2050 của Thủ tướng Chính phủ trong bài toán chủ động thích ứng và kiểm soát tài nguyên nước vùng TGLX.

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

  • Nghiên cứu sinh & Nhà khoa học Thủy văn: Tiếp cận cấu trúc mã nguồn Fortran và thuật toán triệt tiêu khuếch tán số bằng phương pháp đường đặc trưng cho các chất phi bảo toàn.
  • Kỹ sư Quản lý Khai thác Công trình Thủy lợi: Sở Nông nghiệp & PTNT, Công ty TNHH MTV Khai thác Thủy lợi Kiên Giang và An Giang có được quy trình vận hành định lượng hệ thống cống đập (Tha La, Trà Sư, Đầm Chích, Hà Giang, Thần Nông).
  • Hợp tác xã & Nông dân vùng TGLX: Nắm bắt quy luật diễn biến chất lượng nước đầu mùa mưa để chủ động thau chua, ém phèn, tối ưu hóa năng suất lúa và tôm.

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 thiết lập và giải quyết thành công bài toán đóng kín hệ phương trình vi phân chuyển tải chất phi bảo toàn trong mạng lưới kênh sông vùng triều bằng cách tích hợp trực tiếp hai quy luật cân bằng nhiệt động khoáng vật: cân bằng Jurbanite ($pH = pSu + pAl + d$) cho phèn nhôm và cân bằng Jarosite/Pyrite ($pH = m \cdot pSu + \alpha \cdot pFe + \beta$) cho phèn sắt.

2. Đổi mới phương pháp luận so với các nghiên cứu quốc tế trước đây? So với phần mềm MIKE 11 (DHI - Đan Mạch) hay mô hình của Johnston et al. (2005 - Úc), mô hình ACID2020 áp dụng kỹ thuật phân rã toán tử kết hợp phương pháp đường đặc trưng cho bài toán hyperbolic thuần túy. Giải pháp này triệt tiêu hoàn toàn hiện tượng khuếch tán số (numerical diffusion) mà không đòi hỏi chi phí bản quyền thương mại đắt đỏ.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì? Sự biến động đảo chiều cực kỳ nhanh chóng của môi trường nước tại kênh Hà Giang: Từ trạng thái axít cực đoan ($pH = 2.03 \div 2.62$) trong nửa đầu tháng 5 đã chuyển hóa nhanh chóng sang trạng thái gần trung tính ($pH = 6.12 \div 6.15$) vào đầu tháng 6. Nguyên nhân do lượng mưa cực lớn ($2.100 \div 2.200$ mm) tập trung nhanh kết hợp với dòng nước ngọt mang phù sa sông Hậu chảy qua kênh Vĩnh Tế làm kết tủa ion sắt và trung hòa axit.

4. Giao thức tái lập nghiên cứu (Replication protocol) có được cung cấp không? Có. Luận án cung cấp cấu trúc chi tiết của tệp điều khiển đầu vào Fvao.txt, mã nguồn các chương trình con (subroutines) tính nhôm và sunfat trong môi trường Fortran, cùng toàn bộ thông số hình học của 1.061 mặt cắt và sơ đồ liên kết 233 nút hợp lưu.

5. Chương trình nghị sự nghiên cứu 10 năm tới được phác thảo thế nào? Tập trung vào 4 hướng chiến lược: (1) Xây dựng bản đồ GIS lý hóa đất phèn toàn ĐBSCL; (2) Tích hợp mô hình dòng thấm ngầm 2D/3D trong đất với mạng lưới kênh 1D; (3) Tự động hóa điều khiển cống xả phèn ven biển Tây; (4) Đánh giá tác động cộng hưởng của biến đổi khí hậu và nước biển dâng đến động lực tiêu phèn.

Kết luận

  1. Luận án đã giải quyết trọn vẹn bài toán lan truyền chất phi bảo toàn bằng việc xác lập bản chất phèn TGLX là phèn nhôm - sắt và xây dựng thành công hệ phương trình toán học chi phối dựa trên cân bằng Jurbanite và Jarosite.
  2. Xây dựng và hoàn thiện chương trình máy tính chuyên dụng ACID2020 trên nền Digital Visual Fortran 64-bit, ghép nối hoàn hảo với mô hình thủy lực DELTA để mô phỏng định lượng biến thiên các ion $Al^{3+}, Fe^{3+}, SO_4^{2-}$ và pH.
  3. Xác lập quy luật biến động không - thời gian của trục tiêu phèn kênh Hà Giang ra biển Tây trên mạng lưới 1.061 mặt cắt, chỉ rõ cơ chế bùng phát rốn phèn ($pH \approx 2.0$) vào tháng 5 và cơ chế tự pha loãng ($pH > 6.0$) trong tháng 6.
  4. Phát hiện hiện tượng ô nhiễm kép mặn - phèn tại các mặt cắt giáp biển (29A, 30A), cung cấp luận cứ khoa học để vận hành tối ưu hệ thống 30 cống ven biển Tây.
  5. Đề xuất hệ thống giải pháp công trình và phi công trình (rửa phèn một chiều, ém phèn mặt ruộng) có giá trị thực tiễn cao, đóng góp trực tiếp vào sự phát triển nông nghiệp bền vững của vùng Tứ giác Long Xuyên.
  6. Mở ra hướng nghiên cứu mới về mô hình hóa chất lượng nước liên kết đa môi trường (Đất - Nước mặt - Công trình điều tiết) cho toàn bộ vùng đồng bằng châu thổ sông Mê Kông.