Tổng quan về luận án
Nghiên cứu đánh giá rủi ro xâm nhập mặn (XNM) tại khu vực ven biển Thái Bình – Nam Định ứng với các kịch bản mực nước triều và nước biển dâng (NBD) là công trình luận án tiến sĩ kỹ thuật chuyên ngành Thủy văn học (mã số: 9440224) của nghiên cứu sinh Nguyễn Văn Đào, được thực hiện dưới sự hướng dẫn khoa học của PGS. Nguyễn Mai Đăng và GS. Trần Hồng Thái tại Trường Đại học Thủy Lợi (2023). Đặt trong bối cảnh biến đổi khí hậu (BĐKH), suy giảm dòng chảy thượng nguồn và áp lực khai thác nước ngầm gia tăng, công trình giải quyết một khoảng trống nghiên cứu cốt lõi: sự thiếu hụt phương pháp luận chuẩn xác để mô hình hóa và bản đồ hóa rủi ro XNM vào tận hệ thống thủy lợi nội đồng cấp xã khi không có đầy đủ số liệu mặt cắt hình học và quan trắc độ mặn chi tiết trên toàn mạng lưới kênh mương.
Luận án tập trung giải quyết 3 câu hỏi nghiên cứu then chốt:
- Làm thế nào để xác định nguy cơ và lập bản đồ hiểm họa XNM trong nội đồng khi việc lấy nước từ sông chính vào phục vụ nông nghiệp phụ thuộc chặt chẽ vào quy trình vận hành và thời gian đóng mở của các cống đầu mối?
- Bộ tiêu chí nào phản ánh toàn diện tính dễ bị tổn thương (TDBTT) do XNM cho các xã ven biển thuần nông, và trọng số của từng tiêu chí thành phần được lượng hóa như thế nào?
- Quy trình tích hợp không gian đa tiêu chí nào cho phép tính toán định lượng chỉ số rủi ro XNM chi tiết đến 243 xã ven biển phục vụ công tác quy hoạch và ứng phó thiên tai?
Luận án vận dụng khung lý thuyết tích hợp đánh giá rủi ro thiên tai của Ủy ban Liên chính phủ về Biến đổi Khí hậu (IPCC, 2007; IPCC, 2012), kết hợp phương pháp luận động lực dòng chảy 1D (mô hình MIKE 11 HD/AD), khung phân tích quan hệ nhân quả DPSIR, nguyên tắc lựa chọn chỉ thị SMART và phương pháp phân tích thứ bậc AHP của Thomas L. Saaty. Phạm vi nghiên cứu bao quát toàn bộ 243 xã thuộc 8 huyện ven biển trọng điểm: Nghĩa Hưng, Hải Hậu, Giao Thủy, Xuân Trường, Trực Ninh (tỉnh Nam Định) và Kiến Xương, Tiền Hải, Thái Thụy (tỉnh Thái Bình) trong chuỗi thời gian mùa kiệt từ tháng I đến tháng III, xét đến các mốc kịch bản NBD 2030, 2040 và 2050 (theo kịch bản RCP4.5 và RCP8.5 của Bộ Tài nguyên và Môi trường). Kết quả nghiên cứu mang lại tác động định lượng đột phá: cung cấp cơ sở khoa học tối ưu hóa quy trình xả nước gia tăng từ các hồ chứa thượng nguồn (Hòa Bình, Tuyên Quang, Thác Bà) với việc chứng minh duy trì mực nước Hà Nội ở mức +1,8 m thay vì +2,2 m trong đợt xả đầu vẫn đảm bảo vận hành lấy nước thau rửa hệ thống cho các tỉnh ven biển, giúp tiết kiệm tới 302 triệu m³ nước xả (tương đương lưu lượng giảm 2.972 m³/s trong 4 ngày xả).
Literature Review và Positioning
Tổng quan nghiên cứu quốc tế và trong nước phản ánh sự phát triển từ việc tiếp cận thuần túy thủy động lực học sang tích hợp quản trị rủi ro đa chiều. Trong dòng nghiên cứu mô phỏng truyền mặn cửa sông:
- Sangman Jeong và cộng sự (2010) sử dụng mô hình số trị 3D EFDC tại hạ lưu sông Geum (Hàn Quốc), chứng minh phạm vi XNM đạt 50,72 km trong điều kiện hạn hán gay gắt và co hẹp về 42,10 km khi có dòng chảy lũ.
- Rice và cộng sự (2012) áp dụng mô hình thủy động lực – phú dưỡng 3D HEM-3D trên sông James và Chickahominy thuộc Vịnh Chesapeake (Hoa Kỳ), lượng hóa rằng mức NBD 30 cm làm số ngày có độ mặn vượt ngưỡng 0,1‰ tăng hơn 100 ngày tại các điểm lấy nước sinh hoạt.
- Wei Chen và cộng sự (2016) trên cửa sông Dương Tử (Trung Quốc) thiết lập mô hình 3D lưới phi cấu trúc, chỉ ra đường đẳng mặn dịch chuyển ngược dòng phi tuyến tính theo các mức NBD từ 0,5 m đến 2 m.
- Lucy M. Bricheno và cộng sự (2021) trên đồng bằng sông Hằng – Brahmaputra – Meghna đã chỉ ra sự khuếch đại phân hóa độ mặn giữa mùa khô và mùa mưa dưới tác động của BĐKH toàn cầu.
Tại Việt Nam, các nghiên cứu của Hoàng Văn Đại và cộng sự (2014) trên sông Mã, Nguyễn Kỳ Phùng và cộng sự (2017) trên sông Đồng Nai, hay Tô Văn Trường và cộng sự (2017), Đỗ Đức Thắng và cộng sự (2019) tại đồng bằng sông Hồng đã làm sáng tỏ quy luật truyền triều và nêm mặn trên các trục sông chính. Tuy nhiên, tồn tại hai luồng quan điểm và hạn chế phương pháp luận nổi cộm:
- Tranh luận về phương pháp xác lập ranh giới mặn: Luồng quan điểm truyền thống chủ yếu dựa trên việc nội suy đường đẳng mặn (1‰, 4‰) dọc theo dòng chính hoặc dùng công thức triết giảm khoảng cách (Lưu Đức Dũng và ccs, 2014; Hoàng Thanh Sơn và ccs, 2018). Luồng quan điểm này bị phê phán vì mang tính chủ quan và định tính khi gán trực tiếp mức độ nhiễm mặn của sông chính cho toàn bộ vùng nội đồng lân cận, bỏ qua hoàn toàn cơ chế lấy nước qua hệ thống công trình thủy nông.
- Tranh luận về khung đánh giá thiệt hại: Nhiều nghiên cứu trước đây chỉ thống kê thiệt hại quá khứ (ex-post) mà chưa xây dựng được khung đánh giá rủi ro đón đầu (ex-ante) kết hợp giữa hiểm họa vật lý và tính dễ tổn thương kinh tế - xã hội theo phân cấp hành chính cấp xã.
Luận án của NCS. Nguyễn Văn Đào đã định vị chính xác khoảng trống này: thiết lập một phương pháp luận bắc cầu (bridging methodology), tích hợp mô hình thủy lực 1D trên sông chính với thuật toán đánh giá hiểm họa thông qua khả năng vận hành lấy nước thực tế của các cống tưới tiêu (tổ hợp giữa độ mặn cực đại $H_S$ và thời gian lấy nước an toàn $H_T$) để truyền dẫn hiểm họa vào nội đồng 243 xã ven biển mà không đòi hỏi chi phí đo đạc mặt cắt kênh mương nội đồng không khả thi.
Đó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 khung lý thuyết đánh giá rủi ro thiên tai của IPCC (2012) và lý thuyết tính dễ bị tổn thương xã hội của Dow & Downing (1995), Weichselgartner & Bertens (2000) trong phân ngành Thủy văn học và Quản lý tài nguyên nước. Cụ thể:
- Tái định nghĩa cấu trúc Hiểm họa Xâm nhập mặn ($H$): Không đồng nhất hiểm họa mặn với nồng độ mặn cực đại ($S_{max}$), luận án đề xuất mô hình hàm tổ hợp $H = f(H_S, H_T)$, trong đó $H_S$ phản ánh cường độ độ mặn đỉnh triều và $H_T$ phản ánh thời lượng lấy nước hữu dụng (giờ lấy nước khả dụng khi nồng độ mặn $S \le S_{cp}$ và mực nước sông $H_{river} \ge H_{gate}$). Đây là bước chuyển dịch từ đánh giá trạng thái tĩnh sang đánh giá tương tác động lực - công trình.
- Khung khái niệm Tính dễ bị tổn thương ($V$): Được cụ thể hóa thành hệ thống 3 thành phần liên kết hữu cơ theo công thức chuẩn hóa:
$$V = f(S, E, A) = \alpha \cdot S + \beta \cdot E - \gamma \cdot A$$
(hoặc tích hợp đa biến chuẩn hóa), trong đó $S$ là Tính nhạy (Sensitivity - cấu trúc nhân khẩu học, tỷ trọng đất lúa, hạ tầng nông nghiệp), $E$ là Mức độ phơi nhiễm (Exposure - diện tích canh tác, tài sản nằm trong vùng ảnh hưởng), và $A$ là Khả năng thích ứng (Adaptive Capacity - thu nhập, tỷ lệ đê kiên cố, diện tích rừng ngập mặn, năng lực thể chế PCTT).
Khung phân tích độc đáo
Khung phân tích của luận án là sự hợp nhất của 4 nền tảng phương pháp luận:
- Khung DPSIR (Động lực - Áp lực - Trạng thái - Tác động - Phản hồi): Đóng vai trò cấu trúc hóa các mối quan hệ nguyên nhân - kết quả giữa quá trình phát triển kinh tế, BĐKH, chế độ thủy văn cửa sông và tổn thất sinh kế.
- Bộ nguyên tắc SMART: Chuẩn hóa việc tuyển chọn các chỉ thị đầu vào, đảm bảo tính cụ thể, đo lường được, khả thi, tương thích và gắn với chu kỳ thời gian mùa kiệt.
- Phương pháp Phân tích thứ bậc AHP: Thiết lập ma trận so sánh cặp, kiểm soát tính nhất quán logic thông qua Chỉ số nhất quán ngẫu nhiên ($RI$) và Tỷ số nhất quán ($CR < 0,1$) để lượng hóa khách quan trọng số các biến thành phần.
- Điều kiện biên giới hạn (Boundary Conditions): Khung phân tích được xác lập trong giới hạn các hệ thống thủy nông có kiểm soát cống tự chảy và bơm tát vùng hạ du sông Hồng - Thái Bình, lấy ngưỡng chịu mặn $S_{cp} \le 1‰$ (tưới lúa và sinh hoạt) và $S_{cp} \le 4‰$ (ngưỡng thiên tai theo Quyết định 46/2014/QĐ-TTg) làm tiêu chuẩn vận hành.
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 nhận thức luận thực chứng (Positivism) kết hợp chủ nghĩa hiện thực phản biện (Critical Realism), sử dụng thiết kế nghiên cứu đa phương pháp định lượng chuyên sâu (Rigorous Quantitative Multi-method Design) bao gồm 3 cấp độ phân tích không gian:
- Cấp vĩ mô: Toàn bộ mạng lưới sông đồng bằng sông Hồng - sông Thái Bình kết nối 9 cửa sông đổ ra biển Đông.
- Cấp trung gian: Hệ thống công trình thủy lợi, cống tưới tiêu dọc các triền sông Đáy, sông Ninh Cơ, sông Hồng, sông Trà Lý, sông Thái Bình.
- Cấp vi mô: 243 đơn vị hành chính cấp xã thuộc 8 huyện ven biển của 2 tỉnh Thái Bình và Nam Định.
Quy trình nghiên cứu rigorous
Quy trình nghiên cứu tuân thủ các giao thức kiểm định nghiêm ngặt:
- Giao thức thu thập và chuẩn hóa dữ liệu: Số liệu địa hình mặt cắt sông chính, chuỗi thủy văn - mặn thực đo mùa kiệt (tháng I - III) từ hệ thống trạm quốc gia (Sơn Tây, Hà Nội, Thượng Cát, Ba Lạt, Hòn Dáu...). Dữ liệu kinh tế - xã hội, hiện trạng đê điều, diện tích rừng ngập mặn, cơ cấu cây trồng được điều tra tại các Công ty Khai thác Công trình Thủy lợi (Xuân Thủy, Hải Hậu, Nghĩa Hưng, Nam Ninh...) và Chi cục Thủy sản, Nông nghiệp 2 tỉnh.
- Tam giác đạc dữ liệu (Data Triangulation): Đối chiếu chéo giữa số liệu mô phỏng mô hình toán với tài liệu đo đạc thực tế tại các cống đầu mối (cống Ngô Đồng, Cồn Nhất, cống Số 7, cống Trệ, cống Tam Tòa, cống Thái Phúc, cống Bùi...) và kết quả phỏng vấn chuyên gia quản lý khai thác thủy nông.
- Kiểm soát độ tin cậy AHP: Ma trận so sánh cặp so sánh mức độ ưu tiên của các biến tính nhạy ($S$), phơi nhiễm ($E$) và khả năng ứng phó ($A$) thỏa mãn nghiêm ngặt điều kiện tỷ số nhất quán $CR = \frac{CI}{RI} < 0,10$.
+-------------------------------------------------------------------------------+
| QUY TRÌNH ĐÁNH GIÁ RỦI RO XÂM NHẬP MẶN |
+-------------------------------------------------------------------------------+
|
+--------------------------------+-------------------------------+
| |
v v
[MÔ PHỎNG THỦY LỰC - TRUYỀN MẶN] [ĐÁNH GIÁ TÍNH DỄ TỔN THƯƠNG]
- MIKE 11 HD & AD trên mạng sông chính - Khung DPSIR & Chỉ thị SMART
- Kịch bản Triều (P = 1% ÷ 50%) - Trọng số hóa bằng AHP (CR < 0.1)
- Kịch bản NBD (RCP4.5, RCP8.5: 2030-2050) - 3 Trụ cột: Nhạy (S), Phơi nhiễm (E),
| Ứng phó (A)
v |
[HIỂM HỌA XÂM NHẬP MẶN (H)] |
- Tổ hợp: Độ mặn cực đại (Hs) v
& Thời gian lấy nước an toàn (Ht) [CHỈ SỐ TÍNH DỄ TỔN THƯƠNG (V)]
| |
+--------------------------------+-------------------------------+
|
v
[BẢN ĐỒ RỦI RO XÂM NHẬP MẶN (R)]
- Tích hợp không gian GIS (R = H x V)
- Phân cấp chi tiết 243 xã ven biển
Data và phân tích
- Thiết lập mô hình MIKE 11:
- Mô-đun Thủy lực (HD) giải hệ phương trình Saint-Venant 1 chiều kết hợp mô-đun Truyền chất đối lưu khuếch tán (AD).
- Biên trên: Lưu lượng dòng chảy tại các trạm Hòa Bình, Yên Bái, Thác Bà, Hàm Yên, Na Hang, Liễu Sơn, Phú Cường, Thác Bưởi, Cầu Sơn, Chũ, Sơn Tây.
- Biên dưới: Mực nước thủy triều 1 giờ/lần tại 9 cửa sông chính (Đáy, Ninh Cơ, Ba Lạt, Trà Lý, Thái Bình, Văn Úc, Lạch Tray, Cấm, Bạch Đằng) kết hợp các tần suất triều $P = 1%, 5%, 10%, 20%, 25%, 50%$ và các kịch bản NBD trích xuất từ Kịch bản BĐKH của Bộ TNMT cho dải ven biển Hòn Dáu – Đèo Ngang (năm 2030: 12-13 cm; 2040: 17-18 cm; 2050: 22-25 cm).
- Hiệu chỉnh và kiểm định mô hình:
- Đánh giá thông qua hệ số tương quan ($R^2$) và chỉ số hiệu quả mô hình Nash-Sutcliffe ($NSE$). Giai đoạn hiệu chỉnh lưu lượng tại Sơn Tây và Hà Nội (1/I - 31/III/2012) và kiểm định mực nước tại Sơn Tây, Thượng Cát (1/I - 31/III/2011) cho thấy đường quá trình tính toán khớp chặt chẽ với thực đo.
- Mô-đun AD truyền mặn được hiệu chỉnh và kiểm định với số liệu độ mặn thực đo tại các trạm kiểm tra vùng cửa sông, đạt chỉ số $NSE > 0,75$, xác nhận mô hình có độ tin cậy khoa học cao trong tái hiện nêm mặn mùa cạn.
Phát hiện đột phá và implications
Những phát hiện then chốt
- Quy luật lan truyền mặn trên các trục sông chính:
- Trong điều kiện mùa kiệt trung bình, chiều dài xâm nhập của ranh mặn 1‰ trên các sông là: sông Hồng từ 14 – 33 km, sông Đáy 41 km, sông Ninh Cơ 36 km, sông Trà Lý 51 km và sông Thái Bình từ 13 – 49 km.
- Khi xuất hiện tổ hợp triều kiệt cực đoan (kịch bản $P = 1%$), ranh mặn xâm nhập sâu bất thường: tại cống Số 7 trên triền sông Hồng (cách biển 28 km), độ mặn cực đại ghi nhận đạt 3,75‰ (vượt ngưỡng cho phép lấy nước là 2,75‰); cống Trệ trên sông Ninh Cơ (cách biển 29 km) đạt 3,0‰ (ngưỡng cho phép 2,0‰); cống Tam Tòa trên sông Đáy (cách biển 35 km) đạt 2,0‰ (ngưỡng cho phép 1,0‰).
- Sự suy giảm nghiêm trọng của thời lượng lấy nước an toàn ($H_T$):
- Hiện tượng XNM làm tê liệt khả năng lấy nước tự chảy. Điển hình tại cống Ngô Đồng (Nam Định), trong chu kỳ triều tháng III, tổng thời gian có độ mặn $> 1‰$ lên tới gần 320 giờ, trong khi khoảng thời gian nước ngọt khả dụng ($S \le 1‰$) chỉ vỏn vẹn 40 giờ (tương tự tại cống Cồn Nhất chỉ đạt 43,9 giờ). Nhiều cống ven biển mỗi ngày chỉ tranh thủ mở được 15 – 30 phút vào chân triều thấp, buộc các trạm bơm nội đồng phải vận hành cưỡng bức.
- Phân hóa không gian tính dễ tổn thương và rủi ro:
- Kết quả đánh giá trên 243 xã cho thấy: mức độ rủi ro cao và rất cao tập trung đậm đặc tại các huyện trực tiếp giáp biển như Tiền Hải, Thái Thụy (Thái Bình) và Giao Thủy, Hải Hậu, Nghĩa Hưng (Nam Định).
- Tại huyện Thái Thụy, 91% diện tích chịu nhiễm mặn $\ge 1‰$; huyện Tiền Hải có 100% diện tích nằm trong vùng nguy cơ. Ngược lại, các xã nội địa thuộc Kiến Xương, Xuân Trường, Trực Ninh chịu mức độ rủi ro trung bình và thấp hơn nhờ khoảng cách tiêu tán mặn và năng lực công trình ngăn mặn tốt hơn.
- Tác động của nước biển dâng theo kịch bản RCP4.5 và RCP8.5:
- Đến năm 2050, khi mực nước biển dâng thêm 22 – 25 cm, phạm vi ảnh hưởng của hiểm họa mặn cấp cao mở rộng sâu vào nội địa thêm từ 3,5 – 6,0 km tùy phân lưu, làm gia tăng số lượng xã rơi vào nhóm rủi ro rất cao từ 18% lên hơn 32% tổng số xã nghiên cứu.
| Chỉ số / Thông số đặc trưng |
Giá trị định lượng |
Địa bàn / Trạm quan trắc |
Ghi chú & Nguồn dẫn chứng |
| Phạm vi xâm nhập ranh mặn 1‰ |
14 – 51 km |
Các cửa sông chính ĐBSH |
Trà Lý xa nhất (51 km), Đáy (41 km) |
| Độ mặn cực đại đo đạc mùa kiệt |
3,75‰; 3,0‰; 2,0‰ |
Cống Số 7, cống Trệ, cống Tam Tòa |
Vượt xa ngưỡng lấy nước cho phép |
| Thời gian lấy nước an toàn ($S \le 1‰$) |
40,0 – 43,9 giờ |
Cống Ngô Đồng, Cồn Nhất |
Ghi nhận trong đợt triều kiệt tháng 3 |
| Diện tích chịu ảnh hưởng mặn 1‰ |
91% và 100% |
Huyện Thái Thụy, Tiền Hải |
Số liệu thực địa vùng ven biển Thái Bình |
| Tiết kiệm nguồn nước xả hồ chứa |
302 triệu m³ |
Hệ thống hồ Hòa Bình - Tuyên Quang |
Giảm mực nước Hà Nội về +1,8 m |
| Đê biển đã kiên cố hóa chống bão |
56,8 / 76,6 km |
Tỉnh Nam Định |
Đáp ứng bão cấp 10, tần suất 5% |
| Diện tích rừng ngập mặn phòng hộ |
3.600 ha |
Vùng bờ biển Nam Định |
Đóng góp trực tiếp vào chỉ số $A$ |
Implications đa chiều
- Về mặt lý thuyết và phương pháp: Chứng minh tính ưu việt của việc kết hợp chỉ số thời gian vận hành công trình ($H_T$) vào bản đồ hiểm họa XNM, cung cấp phương pháp luận tổng quát có thể nhân rộng cho các đồng bằng châu thổ ven biển có hệ thống đê cống phức tạp như Đồng bằng sông Cửu Long, đồng bằng sông Mã, sông Cả.
- Về mặt thực tiễn và chính sách:
- Cung cấp cơ sở khoa học để Bộ Nông nghiệp & PTNT cùng Tập đoàn Điện lực Việt Nam (EVN) tái cơ cấu quy trình vận hành liên hồ chứa lưu vực sông Hồng trong các đợt điều tiết nước đổ ải vụ Đông Xuân.
- Định hướng các địa phương ven biển chuyển đổi mô hình sinh kế và cơ cấu mùa vụ theo chuỗi thích ứng không gian: ưu tiên mô hình rừng ngập mặn kết hợp nuôi trồng thủy sản ở vùng bãi triều; chuyển đổi đất lúa kém hiệu quả thường xuyên nhiễm mặn sang cây màu chịu hạn mặn hoặc thủy sản nước lợ; bố trí diện tích lúa gieo thẳng tập trung (như 6.000 ha tại Thái Thụy) để né mặn đầu vụ.
Limitations và Future Research
Luận án thẳng thắn chỉ rõ 3 giới hạn nghiên cứu chính:
- Giới hạn chiều không gian mô hình thủy lực: Sử dụng mô hình 1D MIKE 11 cho toàn hệ thống sông Hồng - Thái Bình phản ánh tốt diễn biến dọc trục sông nhưng chưa thể hiện được cấu trúc phân tầng độ mặn 3D (nêm mặn, hoàn lưu trọng lực) cục bộ ngay tại vùng cửa sông tiếp giáp biển như các mô hình 3D (MIKE 3, FVCOM, EFDC).
- Số liệu mặt cắt kênh mương nội đồng: Do hạn chế kinh phí điều tra đo đạc trên diện rộng 243 xã, mô hình lan truyền mặn trong nội đồng vẫn phải giả định thông qua cơ chế phân phối nước của các cống đầu mối chính thay vì mô phỏng thủy lực dòng chảy 2D/1D chi tiết trên toàn bộ hệ thống kênh cấp 1, cấp 2.
- Đặc tính chuỗi số liệu mặn: Mạng lưới quan trắc mặn vùng ven biển đồng bằng Bắc Bộ còn thưa (chỉ khoảng 18 trạm đo), chuỗi số liệu đo mặn nhiều trạm bị gián đoạn và không đồng bộ, tạo ra những thách thức nhất định trong việc phân tích thống kê chuỗi cực trị đa năm.
Chương trình nghiên cứu tiếp nối trong tương lai cần tập trung vào:
- Tích hợp mô hình ghép 1D - 2D/3D (couplings) giữa dòng chính và hệ thống ô đồng nội đồng để mô phỏng tương tác xâm nhập mặn mặt ruộng chi tiết.
- Lồng ghép các tác động của sụt lún địa tầng do khai thác nước ngầm và biến động bùn cát đáy sông (xâm thực lòng dẫn) vào quá trình truyền mặn.
- Xây dựng hệ thống hỗ trợ ra quyết định thời gian thực (Real-time Decision Support System - DSS) trên nền tảng WebGIS và AI để cảnh báo sớm rủi ro mặn cho từng cống thủy lợi.
Tác động và ảnh hưởng
- Giá trị học thuật: Công trình đóng góp hệ phương pháp luận hoàn chỉnh và dữ liệu thực nghiệm chuẩn mực cho chuyên ngành Thủy văn học và Biến đổi khí hậu tại Việt Nam, mở ra hướng trích dẫn khoa học cao trong các nghiên cứu đánh giá rủi ro thiên tai ven biển.
- Chuyển đổi ngành Thủy lợi và Nông nghiệp: Cung cấp công cụ tra cứu định lượng cho 8 Công ty Khai thác Công trình Thủy lợi tại Nam Định và Thái Bình trong việc lập lịch đóng mở cống chính xác theo từng con triều, giảm thiểu tối đa rủi ro đưa nước mặn xâm nhập vào đồng ruộng gây suy giảm năng suất lúa.
- Hiệu quả kinh tế - xã hội: Khuyến nghị tối ưu hóa quy trình xả hồ chứa thượng nguồn giúp tiết kiệm hàng trăm triệu mét khối nước ngọt cho phát điện vào cuối mùa khô, đồng thời bảo vệ an ninh nguồn nước cho hơn 38.000 ha đất canh tác ven biển Nam Định và hơn 90% diện tích đất nông nghiệp các huyện ven biển Thái Bình.
Đối tượng hưởng lợi
- Nghiên cứu sinh và Giảng viên ngành Khoa học Trái đất: Tiếp cận khung phương pháp tích hợp AHP - DPSIR - MIKE 11 - GIS để giải quyết bài toán rủi ro thiên tai phức hợp cấp cơ sở.
- Cán bộ quản lý Thủy nông và PCTT cấp tỉnh/huyện: Sử dụng bộ bản đồ rủi ro chi tiết của 243 xã để phân bổ nguồn vốn đầu tư công trình (nâng cấp 76,6 km đê biển Nam Định, nạo vét sông trục, cứng hóa kênh mương) vào đúng các "điểm nóng" có rủi ro cao nhất.
- Chính quyền cấp xã và Nông dân ven biển: Chủ động nắm bắt lịch thời vụ, lịch lấy nước đổ ải vụ Đông Xuân và chuyển đổi cơ cấu cây trồng, vật nuôi phù hợp với ngưỡng chịu mặn của từng tiểu vùng sinh thá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ì?
Luận án đã mở rộng khung lý thuyết đánh giá rủi ro của IPCC (2012) bằng việc định nghĩa lại cấu trúc Hiểm họa Xâm nhập mặn ($H$) thông qua hàm tương tác hai biến $H = f(H_S, H_T)$. Thay vì chỉ căn cứ vào nồng độ mặn cực đại ($H_S$), luận án đưa chỉ số thời gian lấy nước an toàn ($H_T$) qua công trình vào cấu trúc hiểm họa, giải quyết triệt để hạn chế của việc nội suy đường đẳng mặn hình học truyền thống.
2. Tính đổi mới về phương pháp luận so với các nghiên cứu tiền nhiệm?
So với các nghiên cứu của Hoàng Đức Nghĩa (2011) và Nguyễn Thanh Hùng (2013) vốn chỉ dừng lại ở việc lập bản đồ nồng độ mặn dọc sông chính bằng MIKE 11, luận án đã tích hợp đồng thời: mô phỏng thủy lực AD với phương pháp AHP (có kiểm soát nghiêm ngặt chỉ số $CR < 0,1$), khung DPSIR và công nghệ GIS để nội suy hiểm họa và tính dễ tổn thương vào không gian nội đồng của 243 xã một cách có cơ sở khoa học chặt chẽ.
3. Phát hiện bất ngờ và có giá trị thực tiễn cao nhất?
Phát hiện về hiệu quả điều tiết xả nước từ hồ chứa: Luận án chứng minh rằng trong đợt xả ải đầu tiên, việc giảm mực nước khống chế tại Hà Nội từ +2,2 m xuống +1,8 m vẫn hoàn toàn đáp ứng nhu cầu lấy nước thau rửa hệ thống của các tỉnh ven biển cuối nguồn (Thái Bình, Nam Định), giúp tiết kiệm tới 302 triệu m³ nước ngọt xả từ các hồ chứa thủy điện Hòa Bình, Tuyên Quang, Thác Bà.
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 ma trận so sánh cặp AHP, bảng trọng số các biến ($S, E, A$), hệ số nhám Manning trên từng đoạn sông, điều kiện biên thủy lực 9 cửa sông và quy trình chuẩn hóa số liệu cấp xã, cho phép tái lập quy trình tính toán cho các vùng đồng bằng duyên hải khác.
5. Định hướng nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?
Xây dựng mô hình số trị tích hợp đa quy mô 1D-2D-3D kết nối hoàn lưu biển ven bờ vịnh Bắc Bộ với mạng lưới kênh mương nội đồng, tích hợp dữ liệu viễn thám độ ẩm đất và thuật toán học máy (Machine Learning) để dự báo xâm nhập mặn thời gian thực phục vụ nông nghiệp thông minh thích ứng biến đổi khí hậu.
Kết luận
Luận án tiến sĩ của NCS. Nguyễn Văn Đào là công trình khoa học công phu, giải quyết trọn vẹn và xuất sắc bài toán đánh giá rủi ro xâm nhập mặn vùng ven biển Thái Bình – Nam Định với 5 đóng góp cốt lõi:
- Xác lập thành công phương pháp luận mới đánh giá hiểm họa XNM nội đồng thông qua tổ hợp nồng độ mặn và thời gian vận hành lấy nước qua công trình thủy lợi.
- Xây dựng bộ tiêu chí và lượng hóa chính xác trọng số đánh giá Tính dễ bị tổn thương ($V$) do XNM cho 243 xã thuần nông ven biển bằng phương pháp AHP.
- Thiết lập hệ thống bản đồ hiểm họa, bản đồ tổn thương và bản đồ rủi ro XNM chi tiết cấp xã ứng với các kịch bản tần suất triều ($P = 1% \div 50%$) và kịch bản nước biển dâng (RCP4.5, RCP8.5 đến năm 2050).
- Cung cấp bằng chứng định lượng thuyết phục giúp tối ưu hóa quy trình xả nước liên hồ chứa thượng nguồn sông Hồng, tiết kiệm 302 triệu m³ nước xả trong thời kỳ đổ ải vụ Đông Xuân.
- Đề xuất hệ thống giải pháp công trình và phi công trình thích ứng mặn có tính khả thi cao, mở ra 3 hướng nghiên cứu mới về mô hình hóa tích hợp đa quy mô, kinh tế học thích ứng BĐKH và quản trị rủi ro thiên tai dựa vào cộng đồng cho các dải ven biển châu thổ trên toàn cầu.