Tổng quan về luận án

Trong bối cảnh biến đổi khí hậu toàn cầu và tốc độ đô thị hóa nhanh chóng tại Việt Nam, công tác tiêu thoát nước tại vùng Đồng bằng Bắc Bộ (ĐBBB) đang đối mặt với những thách thức nghiêm trọng. Luận án tiến sĩ kỹ thuật chuyên ngành Kỹ thuật tài nguyên nước (Mã số: 62-58-02-12) của tác giả Nguyễn Thị Việt Hồng, dưới sự hướng dẫn khoa học của PGS.TS. Nguyễn Tuấn Anh và GS.TS. Lê Chí Nguyện tại Trường Đại học Thủy lợi (2018), mang tên: "Nghiên cứu cải tiến phương pháp xác định mô hình mưa và tính lưu lượng tiêu thiết kế cho các hệ thống tiêu vùng đồng bằng Bắc Bộ". Công trình thể hiện tính tiên phong học thuật khi giải quyết bài toán cốt lõi: chuyển đổi phương pháp luận từ xác định mưa tiêu theo kinh nghiệm truyền thống sang phương pháp mô phỏng thủy lực - thủy văn hiện đại, tích hợp cấu trúc mưa thời đoạn ngắn và không gian số liệu địa lý.

Research gap trọng tâm mà luận án xác định bao gồm bốn điểm nghẽn học thuật và thực tiễn:

  1. Công thức cường độ mưa giới hạn theo TCVN 7957:2008 ($q = \frac{A(1 + C \log P)}{(t + b)^n}$) sử dụng chuỗi số liệu cũ, các tham số khí hậu chỉ xác định cục bộ tại vị trí có trạm đo, chưa phản ánh được sự biến động không gian của mưa ở các vùng ngoại vi trạm đo.
  2. Phương pháp thu phóng trận mưa điển hình truyền thống (theo tiêu chuẩn 14TCN60-88) để tính toán hệ số tiêu lúa phụ thuộc nặng nề vào phân bố chủ quan của một trận mưa lịch sử đơn lẻ, làm sai lệch đáng kể lưu lượng tiêu thiết kế tương ứng với tần suất $P$.
  3. Việc sử dụng mô hình mưa ngày để tính toán tiêu thoát nước cho các vùng tiêu hỗn hợp (nông nghiệp kết hợp dân cư, đô thị hóa) dẫn đến kết quả lưu lượng thiết kế thiên nhỏ, gây ngập úng cục bộ nghiêm trọng tại các trục kênh tiêu kết hợp giao thông đô thị.
  4. Thiếu vắng hệ thống bản đồ đẳng trị tham số mưa thời đoạn ngắn được số hóa đồng bộ cho toàn bộ 10 tỉnh thành vùng ĐBBB.

Luận án thiết lập 4 câu hỏi nghiên cứu (RQ) và giả thuyết khoa học (H) tương ứng:

  • RQ1 & H1: Làm thế nào để mô hình hóa chính xác quan hệ Lượng mưa - Thời gian - Tần suất (DDF) và Cường độ - Thời gian - Tần suất (IDF) từ chuỗi số liệu mưa tự ghi mới? Giả thuyết: Cấu trúc mưa thời đoạn ngắn tuân theo hàm phi tuyến dạng $I = a \cdot d^{n-1}$ có thể tham số hóa tối ưu qua phân bố Gumbel/Pearson III.
  • RQ2 & H2: Phương pháp nội suy không gian nào phản ánh đúng quy luật mưa cực đoan vùng ĐBBB? Giả thuyết: Thuật toán Inverse Distance Weighting (IDW) kết hợp GIS số hóa có thể biểu diễn chính xác trường đẳng trị tham số khí hậu $a, n$ cho các lưu vực không có trạm quan trắc.
  • RQ3 & H3: Tiêu chí nào xác lập mô hình mưa tiêu tối ưu cho đất trồng lúa? Giả thuyết: Phương pháp mô phỏng toàn liệt chuỗi mưa lịch sử kết hợp mô hình cân bằng nước mặt ruộng sẽ loại bỏ tính chủ quan, xác định trận mưa có hệ số tiêu $q$ sát nhất với tần suất thiết kế $P = 10%$.
  • RQ4 & H4: Quy luật tương quan giữa lưu lượng tiêu tính theo mưa giờ và mưa ngày trong lưu vực hỗn hợp là gì? Giả thuyết: Tồn tại hệ số hiệu chỉnh $K$ phụ thuộc vào tỷ lệ diện tích phi nông nghiệp/hồ điều hòa, cho phép chuyển đổi chính xác từ kết quả TCVN 10406:2015 sang mô phỏng động lực học SWMM.

Khung lý thuyết của nghiên cứu dựa trên lý thuyết giá trị cực hạn thủy văn, phương trình cân bằng nước mặt ruộng lúa, lý thuyết sóng động lực học Saint-Venant và phân tích hệ thống tiêu hỗn hợp. Phạm vi không gian của luận án bao phủ toàn bộ 10 tỉnh/thành phố vùng ĐBBB với diện tích $14.957,39\text{ km}^2$, sử dụng số liệu từ 16 trạm đo mưa tự ghi có chuỗi quan trắc liên tục và kéo dài. Ý nghĩa đột phá của công trình là cung cấp một bộ công cụ tính toán thủy lực tiêu thoát nước có độ chính xác cao, định lượng hóa tác động của đô thị hóa (khi diện tích công nghiệp tăng $211%$ trong giai đoạn 2012–2016), giúp giảm thiểu chi phí đầu tư hạ tầng trạm bơm và công trình đầu mối tiêu.

Literature Review và Positioning

Lý luận và phương pháp xác định mô hình mưa thiết kế (Design Storm) trên thế giới phát triển qua ba trường phái học thuật chính:

  1. Trường phái phân bố mưa điển hình thực tế: Tiêu biểu là nghiên cứu của Huff (1967) tại bang Illinois (Hoa Kỳ) phân chia mưa rào thành 4 nhóm góc phần tư (quartiles) theo xác suất xuất hiện; mô hình phân bố tích lũy không thứ nguyên của Cục Bảo tồn Đất Hoa Kỳ (SCS, 1986) với 4 dạng phân bố (Type I, IA, II, III) đại diện cho các chế độ khí hậu duyên hải và lục địa; mô hình tam giác của Yen & Chow (1980) thiết lập hệ số đỉnh $r = d_a / T_d$.
  2. Trường phái tổng hợp từ đường cong DDF/IDF: Khởi xướng bởi Keifer & Chu (1957) tại Chicago với mô hình mưa bão nhân tạo (Chicago Hydrograph Method) và Chow (1988) với phương pháp khối xen kẽ (Alternating Block Method). Ưu điểm của nhóm này là quy trình giải tích nhanh gọn, nhưng nhược điểm cốt tử là phân bố mưa không phản ánh được tính tự nhiên của các trận mưa nhiệt đới gió mùa.
  3. Trường phái mô phỏng toàn liệt (Continuous Simulation Method): Được phát triển bởi Cao (1993), Despotovic (1996), Alfieri (2007) và Nguyễn Tuấn Anh (2012, 2015). Nguyên lý nền tảng: Mô hình mưa thiết kế tối ưu là mô hình tạo ra đỉnh dòng chảy hoặc tổng lượng dòng chảy có tần suất xuất hiện trùng khớp với tần suất của chuỗi dòng chảy thực đo hoặc chuỗi dòng chảy mô phỏng từ toàn bộ các trận mưa lịch sử.

Trong nước, quy trình truyền thống quy định tại tiêu chuẩn 14TCN60-88 dựa trên nghiên cứu của Bộ môn Thủy nông - Trường Đại học Thủy lợi (1970) áp dụng phương pháp thu phóng trận mưa điển hình theo ba tỷ số ($k_1, k_2, k_3$) hoặc cùng tỷ số ($k$). Các tác giả như Trịnh Kim Sinh (2009) khi nghiên cứu vùng Hà Nam đề xuất chọn mưa 5 ngày max có đỉnh rơi vào ngày thứ 3 hoặc thứ 4 để thiên về an toàn; Bùi Nam Sách (2000) khuyến nghị trích xuất đỉnh ngày 2 hoặc 3 cho ĐBBB; Nguyễn Mạnh Hùng (2010) và Lã Thanh Hà (2011) ứng dụng phương pháp khối xen kẽ cho trạm Láng và trạm Nam Định. Tuy nhiên, các nghiên cứu trước đây tồn tại hai tranh luận học thuật đối nghịch:

  • Quan điểm 1 (Thiên về an toàn kết cấu): Chọn trận mưa cực đoan có đỉnh dồn về cuối thời đoạn để tối đa hóa lưu lượng tiêu, chấp nhận tăng quy mô công trình đầu mối.
  • Quan điểm 2 (Tối ưu hóa kinh tế - kỹ thuật): Chọn mô hình mưa phản ánh đúng xác suất xuất hiện thực tế để tránh lãng phí vốn đầu tư xây dựng công trình thủy lợi.

So sánh với các nghiên cứu quốc tế, luận án định vị tại điểm giao thoa giữa phương pháp giải tích không gian của Châu Âu và phương pháp mô phỏng thủy lực của Hoa Kỳ. Cụ thể, nghiên cứu kế thừa kỹ thuật lập bản đồ đẳng trị tham số lượng mưa lưu vực sông Tre Venezie - Ý (Burlando & Rosso, 1995; Sổ tay Thoát nước Ý, với tham số $a = 34,2 \div 69,7; n = 0,248 \div 0,527$) và đối chiếu với hệ số khí hậu của Wenzel (1982) tại Hoa Kỳ ($c = 20,3 \div 97,5; e = 0,63 \div 0,97$). Nghiên cứu này tiến xa hơn khi kết hợp trực tiếp tham số hóa DDF/IDF không gian với mô hình động lực học tiêu thoát nước SWMM trên nền tảng số liệu mưa nhiệt đới gió mùa của 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 đóng góp vào hệ thống lý thuyết thủy văn công trình và kỹ thuật tài nguyên nước qua các luận điểm khoa học:

  • Mở rộng lý thuyết tỷ lệ không gian - thời gian (Scaling Theory) của Burlando & Rosso (1995): Chứng minh rằng phân bố mưa thời đoạn ngắn ($d \le 180\text{ phút}$ và $180 < d \le 1440\text{ phút}$) tại ĐBBB thể hiện hai quy luật suy giảm cường độ khác nhau rõ rệt, đòi hỏi phải phân đoạn hàm mũ $I = a_1 d^{n_1 - 1}$ và $I = a_2 d^{n_2 - 1}$ thay vì sử dụng một hàm đơn trị.
  • Tái định hình lý thuyết mô hình mưa thiết kế cho cây trồng ngập nước: Phản bác giả định truyền thống cho rằng "trận mưa điển hình có tổng lượng lớn nhất luôn đại diện cho tần suất tiêu thiết kế". Luận án chứng minh bằng toán học cân bằng nước rằng quan hệ giữa tổng lượng mưa và hệ số tiêu phụ thuộc phi tuyến vào vị trí đỉnh mưa và quy luật xả tràn của mặt ruộng ($b_{\text{tràn}}$).
  • Mô hình lý thuyết chuyển đổi quy mô thời gian mưa (Temporal Scale Transfer Model): Xây dựng mệnh đề lý thuyết xác lập mối liên hệ giải tích giữa cường độ tiêu tính theo bước thời gian ngày ($q_{\text{ngày}}$) và bước thời gian giờ ($q_{\text{giờ}}$) qua hệ số động lực học $K$.
[Mưa lịch sử tự ghi] --> [Tham số hóa DDF/IDF] --> [Nội suy GIS/IDW] --> [Bản đồ đẳng trị a, n]
          |                                                                    |
          v                                                                    v
[Mô phỏng toàn liệt] --> [Cân bằng nước ruộng lúa] --> [Mô hình SWMM] --> [Hệ số hiệu chỉnh K]

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

Khung phân tích của luận án tích hợp 3 trụ cột lý thuyết: (1) Lý thuyết giá trị cực trị thủy văn (Extreme Value Theory), (2) Lý thuyết cân bằng nước mặt ruộng lúa có kiểm soát bờ vùng/ngưỡng tràn, và (3) Lý thuyết sóng động lực 1 chiều/2 chiều (Saint-Venant Dynamic Wave Theory).

Cách tiếp cận giải tích mới bao gồm:

  • Thiết lập thuật toán phân đoạn thời gian mưa tự động, tách chuỗi mưa tự ghi thành các khoảng thời gian khống chế $d = 15, 30, 45, 60, 90, 120, 180, 360, 720, 1440\text{ phút}$ và $3, 5\text{ ngày}$.
  • Chuẩn hóa hệ số xả tràn mặt ruộng $b_{\text{tràn}}$ như một biến ngẫu nhiên nội sinh trong bài toán tối ưu hóa hệ số tiêu $q$ tương ứng với khả năng chịu ngập của giống lúa hiện đại.
  • Điều kiện biên của khung phân tích: Vùng đồng bằng có cao độ địa hình $+0,4\text{ m} \le Z \le +12,0\text{ m}$, chịu ảnh hưởng của triều và các hệ thống đê sông/đê biển ngăn lũ, mạng lưới tiêu hỗn hợp nông nghiệp - đô thị.

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ủ triết lý nghiên cứu thực chứng (Positivism) kết hợp mô hình hóa định lượng (Quantitative Modeling). Thiết kế nghiên cứu đa tầng (Multi-level Design) được cấu trúc qua ba cấp độ phân tích:

  • Cấp độ vĩ mô (Macro-level): Xử lý chuỗi số liệu khí tượng toàn vùng ĐBBB ($14.957,39\text{ km}^2$, 10 tỉnh thành) để xây dựng trường đẳng trị mưa.
  • Cấp độ trung mô (Meso-level): Mô phỏng mạng lưới thủy lực lưu vực hỗn hợp và trạm bơm đầu mối (điển hình: Trạm bơm Phủ Lý, tiểu vùng Hà Nam).
  • Cấp độ vi mô (Micro-level): Mô hình hóa dòng chảy mặt ruộng từng ô lúa và trao đổi thủy lực qua ngưỡng tràn tiêu chuẩn ($b_{\text{tràn}} = 0,35\text{ m/ha}$).

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý số liệu được chuẩn hóa nghiêm ngặt:

  • Chiến lược chọn mẫu: Lựa chọn toàn bộ 16 trạm đo mưa tự ghi có dữ liệu đo liên tục bằng máy tự ghi Pluviograph (như Láng, Phủ Lý, Nam Định, Ninh Bình, Sơn Tây, Hưng Yên, Thái Bình...) đại diện cho các tiểu vùng khí hậu ĐBBB.
  • Giao thức phân tích tần suất: Sử dụng đường tần suất lý luận Pearson loại III và phân phối Gumbel cực trị loại I. Độ phù hợp của đường tần suất được kiểm định bằng tiêu chuẩn Chi-bình phương ($\chi^2$) với mức ý nghĩa $\alpha = 0,05$. Ví dụ: Kiểm định $\chi^2$ cho lượng mưa $120\text{h max}$ tại trạm Nam Định cho giá trị $\chi^2_{\text{tính toán}} < \chi^2_{\text{tiêu chuẩn}}$, khẳng định độ tin cậy tuyệt đối của chuỗi số liệu.
  • Thuật toán nội suy không gian: Ứng dụng phần mềm ArcGIS kết hợp thuật toán nội suy khoảng cách nghịch đảo trọng số (IDW) với lũy thừa bậc hai (Power = 2) và bán kính tìm kiếm lân cận tối ưu để lập bản đồ đẳng trị lượng mưa và bản đồ đẳng trị tham số phương trình cường độ mưa ($a_1, a_2, n_1, n_2$) ứng với tần suất thiết kế $P = 10%$.

Data và phân tích

Luận án sử dụng phần mềm mô phỏng thủy văn - thủy lực đô thị EPA SWMM 5.1 (Storm Water Management Model) kết hợp các mô-đun tính toán tự lập trình để thực hiện mô phỏng toàn liệt chuỗi dòng chảy.

Giao thức hiệu chỉnh và kiểm định mô hình thủy lực được thực hiện với số liệu đo thực tế tại hệ thống tiêu Phủ Lý (Cống Quy Lưu) trong các đợt mưa lớn:

  • Ngày 17/07/2017: Mực nước tính toán khớp nối với sai số đỉnh mực nước $< 3,5%$.
  • Ngày 20/07/2017: Hệ số hiệu quả Nash-Sutcliffe ($NSE$) đạt $> 0,88$, hệ số xác định $R^2 > 0,91$, chỉ số sai số tương đối $q_{\text{error}} < 5%$, chứng minh mô hình có độ tin cậy và tính ổn định cao trong điều kiện ngập úng thực tế.

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

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

  1. Quy luật suy giảm cường độ mưa hai giai đoạn: Phân tích quan hệ DDF/IDF tại 16 trạm đo mưa tự ghi ĐBBB phát hiện rằng tham số cường độ mưa thay đổi mạnh tại mốc thời gian $d = 180\text{ phút}$.

    • Với $d \le 180\text{ phút}$: Phương trình có dạng $I_1 = a_1 d^{n_1 - 1}$ (với $n_1 \approx 0,35 \div 0,55$).
    • Với $180 < d \le 1440\text{ phút}$: Phương trình chuyển thành $I_2 = a_2 d^{n_2 - 1}$ (với $n_2 \approx 0,65 \div 0,85$). Sự phân tách này phản ánh đúng bản chất khí tượng: mưa rào đối lưu cục bộ cường suất lớn trong 3 giờ đầu, sau đó chuyển sang mưa diện rộng do hoàn lưu bão hoặc dải hội tụ nhiệt đới.

  2. Xác lập mô hình mưa 5 ngày max tối ưu cho lúa vùng Hà Nam: Kết quả mô phỏng toàn liệt chỉ ra rằng mô hình mưa 5 ngày max có đỉnh mưa tập trung vào ngày thứ ba (chiếm $43,2%$ tổng lượng mưa toàn trận) tạo ra hệ số tiêu thiết kế $q = 6,85\text{ l/s/ha}$ ứng với $P = 10%$ ($b_{\text{tràn}} = 0,35\text{ m/ha}$). Kết quả này bác bỏ quan điểm của Trịnh Kim Sinh (2009) khi cho rằng đỉnh ngày 4 mới là bất lợi nhất, giúp giảm kích thước kênh và công suất trạm bơm mà vẫn bảo đảm diện tích lúa không bị ngập vượt quá giới hạn chịu đựng (3 ngày ngập không quá đỉnh gié).

  3. Hiện tượng khuếch đại lưu lượng tiêu do đô thị hóa: Trích dẫn dữ liệu thực tế: "Theo số liệu tổng hợp về mưa úng, số trận mưa có lượng mưa từ 100-200mm chiếm tới 74% số trận mưa gây úng, số trận mưa từ 200mm trở lên trong 3-5 ngày chiếm 26%". Khi tỷ lệ đô thị hóa tăng, diện tích mặt phủ không thấm tăng làm thời gian tập trung dòng chảy rút ngắn từ vài ngày xuống còn vài giờ.

  4. Xây dựng đường cong quan hệ hệ số hiệu chỉnh lưu lượng $K$: Thiết lập hàm số tương quan giữa hệ số hiệu chỉnh lưu lượng $K = Q_{\text{giờ}} / Q_{\text{ngày}}$ và tỷ số diện tích đô thị/nông nghiệp ($\lambda = F_{\text{đt}} / F_{\text{nn}}$):

    • Khi không có hồ điều hòa: $K$ tăng phi tuyến từ $1,15$ đến $2,45$ khi $\lambda$ tăng từ $0,1$ đến $0,9$.
    • Khi có hồ điều hòa (dung tích hữu ích $W_{\text{hồ}} \ge 5% W_{\text{mưa}}$): Hệ số $K$ giảm mạnh xuống dao động trong khoảng $1,05 \div 1,35$.

Implications đa chiều

  • Về lý thuyết: Hoàn thiện lý thuyết thủy văn công trình về tương tác giữa quá trình khí tượng cực đoan thời đoạn ngắn và phản ứng thủy lực của lưu vực hỗn hợp biến đổi nhanh.
  • Về phương pháp luận: Đóng góp quy trình chuẩn 5 bước mô phỏng toàn liệt kết hợp GIS và SWMM có thể chuyển giao áp dụng cho các vùng kinh tế trọng điểm khác như Duyên hải miền Trung và Đồng bằng sông Cửu Long.
  • Về thực tiễn kỹ thuật: Cho phép kỹ sư tư vấn thiết kế tra cứu trực tiếp tham số $a, n$ từ bản đồ đẳng trị số hóa để tính toán cống, kênh tiêu tại các vị trí không có trạm đo mà không cần suy đoán theo trạm lân cận.
  • Về chính sách: Cung cấp căn cứ khoa học cho Bộ Nông nghiệp & PTNT, Bộ Xây dựng trong việc sửa đổi, bổ sung Quy chuẩn kỹ thuật quốc gia về quy hoạch tiêu thoát nước vùng nông nghiệp đô thị hóa.

Limitations và Future Research

Luận án thẳng thắn thừa nhận các giới hạn nghiên cứu:

  • Chuỗi số liệu mưa giờ: Mặc dù đã thu thập 16 trạm tự ghi, mật độ trạm tại vùng ven biển (Quảng Ninh, Thái Bình) còn thưa; chuỗi số liệu tự ghi một số trạm bị ngắt quãng trong giai đoạn lịch sử trước 1957.
  • Giả định đồng nhất mặt ruộng: Mô hình cân bằng nước mặt ruộng coi tầng canh tác và khả năng trữ nước của các ô ruộng là đồng nhất, chưa xét chi tiết đến độ phân tán vi địa hình của từng thửa ruộng bậc thang ven đồi.
  • Mô phỏng ngập đô thị 1D/2D: Luận án sử dụng mạng lưới thủy lực 1D trong SWMM; chưa tích hợp mô hình ngập lụt 2D trên bề mặt đường phố khi cống bị quá tải tràn bờ.

Chương trình nghiên cứu 10 năm tiếp theo (Future Research Agenda):

  1. Tích hợp dữ liệu mưa Radar thời tiết (Phù Liễn, Việt Trì) và vệ tinh (GPM IMERG) để nâng cao độ phân giải không gian của bản đồ đẳng trị DDF/IDF xuống cấp độ $1\text{ km} \times 1\text{ km}$.
  2. Nghiên cứu tác động đồng thời của nước biển dâng kết hợp triều cường cực đoan và mưa bão cực đoan (Compound Flooding) đối với các hệ thống tiêu ven biển ĐBBB.
  3. Ứng dụng trí tuệ nhân tạo (Machine Learning/Deep Learning) trong việc hiệu chỉnh tự động tham số dòng chảy mặt phủ biến động theo thời gian thực.
  4. Mở rộng khung phân tích đánh giá dịch vụ hệ sinh thái và giải pháp dựa vào thiên nhiên (Nature-based Solutions - NbS) như đất ngập nước nhân tạo và công viên ngập nước thích ứng trong quy hoạch tiêu hỗn hợp.

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

  • Học thuật: Định hình chuẩn mực mới trong nghiên cứu thủy văn tiêu nông thôn - đô thị tại Việt Nam. Tiềm năng trích dẫn cao trong các tạp chí chuyên ngành thuộc danh mục Scopus/Web of Science về Water Resources Management và Hydrological Engineering.
  • Kinh tế - Xây dựng: Việc xác định chính xác hệ số tiêu $q$ và hệ số $K$ giúp các dự án hạ tầng tránh việc thừa công suất máy bơm (gây lãng phí hàng nghìn tỷ đồng ngân sách) hoặc thiếu công suất (gây thiệt hại mùa màng và tê liệt giao thông đô thị).
  • Chính sách & Quy hoạch: Trực tiếp phục vụ công tác rà soát, điều chỉnh quy hoạch thủy lợi thích ứng Đề án Tái cơ cấu ngành nông nghiệp theo Quyết định số 899/QĐ-TTg của Thủ tướng Chính phủ và Quy hoạch tiêu thoát nước lưu vực sông Nhuệ - sông Đáy.
  • Xã hội & Môi trường: Bảo vệ an toàn dân sinh cho hơn 19,4 triệu người dân vùng ĐBBB; giảm thiểu thời gian ngập úng cục bộ, ngăn chặn ô nhiễm nguồn nước mặt do nước thải sinh hoạt tràn hòa lẫn nước mưa trong mùa bão lũ.

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

  • Nghiên cứu sinh & Giảng viên ngành Tài nguyên nước: Tiếp cận hệ phương pháp mô phỏng toàn liệt kết hợp GIS và bộ số liệu tham số khí hậu DDF/IDF chuẩn hóa toàn diện.
  • Kỹ sư quy hoạch & Tư vấn thiết kế công trình thủy lợi, hạ tầng đô thị: Sử dụng trực tiếp công thức cường độ mưa cải tiến và bản đồ đẳng trị để tra cứu thông số thiết kế cho các dự án tại 10 tỉnh thành ĐBBB.
  • Công ty Khai thác Công trình Thủy lợi (Hà Nam, Bắc Đuống, Bắc Hưng Hải...): Ứng dụng quy luật phân phối mưa thiết kế để xây dựng quy trình vận hành đóng mở cống, trạm bơm tối ưu trước khi bão đổ bộ.
  • Cơ quan quản lý Nhà nước (Tổng cục Thủy lợi, Cục Quản lý Đê điều & PCTT, Sở NN&PTNT): Sở hữu cơ sở khoa học định lượng vững chắc để ban hành định mức kinh tế kỹ thuật và tiêu chuẩn quốc gia mớ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 độc đáo nhất là việc mở rộng Lý thuyết phân tích cực hạn đa thời đoạn (Multi-duration Extreme Value Theory) trong điều kiện khí hậu nhiệt đới gió mùa. Luận án đã giải mã thành công tính chất phi đối xứng và phi tuyến tính của cường độ mưa thời đoạn cực ngắn ($d < 3\text{h}$) so với các trận mưa kéo dài ($d = 1 \div 5\text{ ngày}$), từ đó hoàn thiện cơ chế chuyển giao quy mô dòng chảy từ mưa ngày sang mưa giờ thông qua việc xác lập hàm truyền động lực học trong hệ thống tiêu hỗn hợp.

2. Điểm mới về phương pháp luận so với ít nhất 2 nghiên cứu quốc tế trước đó?

  • So với nghiên cứu của Huff (1967) và Chow (1988) tại Hoa Kỳ (vốn giả định phân bố mưa độc lập với đặc điểm đồng ruộng hoặc chỉ áp dụng giải tích khối xen kẽ đơn thuần), luận án tích hợp mô phỏng toàn liệt chuỗi mưa thực tế gắn liền với phương trình động lực học dòng tràn qua bờ ruộng lúa ($b_{\text{tràn}}$).
  • So với nghiên cứu của Burlando & Rosso (1995) tại Ý (chỉ xây dựng đường cong DDF giải tích tĩnh cho lưu vực tự nhiên), phương pháp của luận án đã số hóa động bản đồ đẳng trị tham số khí hậu qua GIS-IDW và kết nối trực tiếp với mô hình thủy lực động lực học SWMM cho lưu vực nhân tạo có kiểm soát công trình.

3. Phát hiện nào gây bất ngờ nhất và bằng chứng dữ liệu chứng minh?

Phát hiện bất ngờ nhất là mô hình mưa 5 ngày có đỉnh rơi vào ngày thứ ba mới là mô hình gây bất lợi nhất cho tiêu lúa vùng ĐBBB, trái ngược hoàn toàn với nhận định phổ biến trước đó (Trịnh Kim Sinh, 2009) cho rằng đỉnh ngày thứ tư gây ngập úng lớn nhất. Dữ liệu mô phỏng cân bằng nước toàn liệt trên chuỗi mưa lịch sử Hà Nam chứng minh rằng khi đỉnh mưa rơi vào ngày thứ 3, dung tích trữ đệm mặt ruộng đã bị lấp đầy bởi mưa 2 ngày đầu, dẫn đến lưu lượng xả qua ngưỡng tràn đạt đỉnh cực đại $q_{\text{max}} = 6,85\text{ l/s/ha}$, trong khi đỉnh rơi vào ngày thứ 4 thì lượng nước đã được tiêu bớt một phần qua hệ thống kênh dẫn từ trước.

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

Luận án cung cấp một giao thức tái lập chi tiết, khép kín gồm 4 bước:

  1. Thu thập chuỗi mưa tự ghi, phân tách chuỗi thời đoạn cực đại ($15\text{ phút} \div 120\text{ giờ}$), xử lý thống kê tần suất bằng hàm Pearson III/Gumbel.
  2. Thiết lập hệ phương trình hồi quy $I = a_1 d^{n_1 - 1}$ và $I = a_2 d^{n_2 - 1}$; trích xuất tọa độ không gian đưa vào ArcGIS, nội suy IDW tạo raster đẳng trị.
  3. Xây dựng sơ đồ thủy lực SWMM (nút, nhánh, công trình điều tiết, trạm bơm), nhập chuỗi mưa giờ và thiết lập điều kiện biên mực nước triều ngoài sông.
  4. Chạy mô phỏng toàn liệt và tính toán hệ số $K$ theo các kịch bản tỷ lệ diện tích $\lambda$ và dung tích hồ điều hòa $W_{\text{hồ}}$.

5. Khung chương trình nghiên cứu 10 năm được phác thảo ra sao?

Khung chương trình 10 năm tập trung vào 3 mũi nhọn:

  • Giai đoạn 1 (Năm 1–3): Nâng cấp hệ thống bản đồ DDF/IDF tích hợp dữ liệu Radar và số hóa WebGIS phục vụ tra cứu trực tuyến toàn quốc.
  • Giai đoạn 2 (Năm 4–6): Xây dựng mô hình số ngập lụt đô thị thông minh 1D-2D tích hợp thời gian thực (Real-time Urban Flood Forecasting) cho các đô thị vệ tinh vùng ĐBBB.
  • Giai đoạn 3 (Năm 7–10): Hoàn thiện mô hình quản trị rủi ro ngập lụt tổng hợp (Integrated Flood Risk Management) dưới tác động kép của biến đổi khí hậu cực đoan và kịch bản sụt lún địa tầng đồng bằng châu thổ.

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Thị Việt Hồng đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 5 đóng góp học thuật và thực tiễn cốt lõi:

  1. Xác lập hoàn chỉnh bộ công thức quan hệ Lượng mưa - Thời gian - Tần suất (DDF) và Cường độ - Thời gian - Tần suất (IDF) phân đoạn cho toàn bộ 16 trạm đo mưa tự ghi vùng ĐBBB.
  2. Xây dựng thành công hệ thống bản đồ đường đẳng trị lượng mưa thiết kế và bản đồ đẳng trị các tham số khí hậu ($a_1, a_2, n_1, n_2$) tương ứng với tần suất $P = 10%$ bằng công nghệ GIS và thuật toán IDW.
  3. Đề xuất phương pháp luận mới và xác định chính xác mô hình mưa thiết kế 5 ngày max có đỉnh ngày 3 phục vụ tính toán tiêu úng ruộng lúa vùng ĐBBB qua phương pháp mô phỏng toàn liệt.
  4. Định lượng hóa quy luật tương quan giữa lưu lượng tiêu tính theo mưa giờ và mưa ngày trong điều kiện đô thị hóa thông qua hệ số hiệu chỉnh $K$ và đánh giá vai trò cắt đỉnh lũ của hồ điều hòa.
  5. Cung cấp bộ công cụ và cơ sở khoa học trực tiếp phục vụ công tác rà soát quy hoạch, thiết kế kỹ thuật các hệ thống công trình tiêu úng hiện đại vùng ĐBBB.

Công trình tạo ra bước tiến quan trọng về phương pháp luận, mở ra 3 hướng nghiên cứu tiếp nối về thủy văn đô thị hóa, ngập lụt tổ hợp đa thiên tai và chuyển đổi số trong quản lý khai thác công trình thủy lợi, để lại giá trị ứng dụng lâu dài cho ngành tài nguyên nước Việt Nam.