Tổng quan về luận án

Biến đổi khí hậu (BĐKH) toàn cầu đang làm gia tăng tần suất và cường độ của các hiện tượng thời tiết cực đoan, đặt ra thách thức nghiêm trọng đối với an toàn công trình hạ tầng và an ninh nguồn nước. Tại Việt Nam, dải duyên hải Nam Trung Bộ—kéo dài qua 8 tỉnh, thành phố từ Đà Nẵng đến Bình Thuận—là khu vực có tính chất nhạy cảm sinh thái và thủy văn đặc thù bậc nhất. Địa hình sườn đông dãy Trường Sơn Nam có chiều ngang hẹp (trung bình 40–50 km), phân cắt mạnh với hệ thống sông ngắn dốc, khiến thời gian tập trung dòng chảy lũ diễn ra nhanh, tạo nên những đợt lũ kép, lũ quét khốc liệt gây tổn thất lớn về kinh tế và nhân mạng.

flowchart TD
    subgraph Data_Input ["1. Dữ liệu đầu vào"]
        A["11 Mô hình Khí hậu Toàn cầu GCMs (CMIP5)"]
        B["Chuỗi mưa thực đo tại 93 Trạm Khí tượng Thủy văn"]
    end

    subgraph Climate_Downscaling ["2. Xử lý Khí hậu & Không gian"]
        C["Chi tiết hóa thống kê & Hiệu chỉnh sai số (Quantile Mapping)"]
        D["Nội suy không gian GIS (IDW)"]
        A & B --> C --> D
        E["Lượng mưa 1 ngày lớn nhất Rx1day (RCP4.5 & RCP8.5)"]
        D --> E
    end

    subgraph Hydrological_Modeling ["3. Mô hình hóa Thủy văn Đa quy mô"]
        F["Lưu vực lớn (Vu Gia - Thu Bồn, Sông Ba, Sông Kôn)"]
        G["Lưu vực vừa và nhỏ (Tiểu lưu vực sườn dốc)"]
        E --> F & G
        H["Mô hình bán phân bố NAM + Tương quan Qmax - Qngày max"]
        I["Công thức Sokolovsky & Cường độ giới hạn"]
        F --> H
        G --> I
    end

    subgraph Output_Applied ["4. Kết quả & Ứng dụng Thiết kế"]
        J["Lưu lượng đỉnh lũ thiết kế Qmax tương lai"]
        K["Bản đồ phân vùng biến động dòng chảy lũ 2040-2069 & 2070-2099"]
        H & I --> J --> K
    end

Khoảng trống nghiên cứu then chốt trong thủy văn học công trình tại Việt Nam nằm ở việc thiếu vắng phương pháp luận chi tiết hóa lượng mưa cực trị từ quy mô lưới mô hình toàn cầu về quy mô điểm trạm cho các lưu vực vừa và nhỏ. Các kịch bản BĐKH do Bộ Tài nguyên và Môi trường công bố trước đây (như phiên bản 2012 dựa trên phát thải SRES A1B, B2, A2) chủ yếu sử dụng giá trị trung bình hóa trên không gian rộng và thời gian dài (tháng, mùa), dẫn đến sai số rất lớn khi đưa vào các mô hình toán thủy văn - thủy lực tính toán lũ thiết kế. Luận án giải quyết triệt để khoảng trống này thông qua các câu hỏi và giả thuyết khoa học:

  • Câu hỏi nghiên cứu 1 ($Q_1$): Làm thế nào để loại bỏ sai số hệ thống và chi tiết hóa dữ liệu từ tổ hợp đa mô hình khí hậu toàn cầu về 93 trạm quan trắc địa phương tại Nam Trung Bộ?
    Giả thuyết 1 ($H_1$): Phương pháp hiệu chỉnh sai số phân phối xác suất phi tham số cho phép tái hiện chính xác chế độ mưa cực trị lịch sử và bảo toàn tính đặc thù địa hình ven biển.
  • Câu hỏi nghiên cứu 2 ($Q_2$): Xu thế biến đổi của lượng mưa 1 ngày lớn nhất ($R_{x1day}$) theo các kịch bản nồng độ phát thải đại diện (RCP4.5 và RCP8.5) trong các giai đoạn 2040–2069 và 2070–2099 diễn ra như thế nào trên không gian toàn vùng?
    Giả thuyết 2 ($H_2$): Lượng mưa cực trị 1 ngày gia tăng không đồng nhất, tập trung cao điểm tại khu vực đón gió phía Bắc (Quảng Nam - Đà Nẵng) và phân hóa phức tạp ở phía Nam (Ninh Thuận, Bình Thuận).
  • Câu hỏi nghiên cứu 3 ($Q_3$): Phương thức nào tối ưu hóa việc chuyển đổi lượng mưa tương lai sang lưu lượng đỉnh lũ thiết kế ($Q_{\text{max}}$) cho cả lưu vực lớn và lưu vực vừa/nhỏ?
    Giả thuyết 3 ($H_3$): Tích hợp mô hình toán bán phân bố (NAM) cho lưu vực lớn kết hợp công thức kinh nghiệm giải tích (Sokolovsky, cường độ giới hạn) cho lưu vực nhỏ cung cấp cơ sở tin cậy để lập bản đồ phân vùng biến động dòng chảy lũ.

Khung lý thuyết của công trình tích hợp cơ sở động lực khí hậu của Dự án Đối chứng Mô hình Khí hậu Giai đoạn 5 (CMIP5/IPCC), lý thuyết phân phối xác suất thống kê thủy văn, và lý thuyết hình thành dòng chảy lưu vực (Chow & Maidment, 1988). Với quy mô nghiên cứu bao phủ 93 trạm mưa thực đo, 11 mô hình GCM, phân tích trên 2 giai đoạn tương lai (2040–2069 và 2070–2099), công trình mang lại giá trị định lượng trực tiếp phục vụ quy hoạch đa ngành, an toàn hồ đập và tiêu chuẩn thiết kế công trình hạ tầng.

Literature Review và Positioning

Nghiên cứu về biến đổi lượng mưa và dòng chảy lũ cực đoan đã được phát triển sâu rộng trên thế giới. Frich và cộng sự (2002) cùng Alexander và cộng sự (2006) khi phân tích bộ dữ liệu toàn cầu thế kỷ 20 đã chỉ ra xu thế gia tăng của mưa cực trị tại phần lớn các lục địa, mặc dù có sự suy giảm cục bộ ở Đông Á và Siberia. Ở cấp độ khu vực, Re & Barros (2009) tại Đông Nam Nam Mỹ và Penalba & Robledo (2009) tại lưu vực La Plata ghi nhận sự gia tăng tần suất và cường độ mưa lớn trong các mùa xuân - hè. Về mặt phương pháp luận thủy văn, nền tảng tính toán dòng chảy lũ thiết kế được định hình bởi Chow & Maidment (1988), Shaw (1964), Vijay (2002) và Raghunath (2006) thông qua các mô hình tập trung dòng chảy, đường lũ đơn vị (SCS, Clark, Nash) và đường cong cường độ - thời gian - tần suất (IDF).

Trong y văn tồn tại hai cuộc tranh luận học thuật lớn:

  1. Động lực học khu vực (Dynamical Downscaling) đối đầu Thống kê chi tiết hóa (Statistical Downscaling): Các mô hình RCM độ phân giải cao (như HIRHAM, PRECIS, RegCM) tuy mô phỏng tốt các quá trình vật lý khí quyển nhưng đòi hỏi tài nguyên tính toán khổng lồ và vẫn tồn tại sai số hệ thống quy mô địa phương (Dankers & Feyen, 2008). Ngược lại, phương pháp thống kê hiệu chỉnh sai số (bias correction) trực tiếp theo phân bố xác suất tại trạm mặt đất đem lại hiệu quả vượt trội về thời gian và độ chính xác điểm trạm (Gudmundsson et al.).
  2. Tính bất định của biến đổi dòng chảy lũ cực trị: Kết quả quốc tế cho thấy tác động của BĐKH lên dòng chảy lũ không đồng nhất. Tại Na Uy, cơ quan quản lý yêu cầu cập nhật chu kỳ lũ thiết kế mỗi 15–20 năm với biên độ điều chỉnh tăng từ 0% đến 40% tùy lưu vực. Tại Anh, Kay và cộng sự (2006) chỉ ra lũ đỉnh 50 năm tăng tới 50% ở phía Tây/Bắc nhưng lại có xu hướng giảm ở phía Đông/Nam. Tại Phần Lan, Veijalainen & Vehvilainen (2008) nhận thấy BĐKH ít tác động đến lũ thiết kế hồ đập; trong khi tại Hàn Quốc, Hyun-Han Kwon và cộng sự (2011) áp dụng mô hình trạng thái thời tiết ngẫu nhiên kết hợp chuỗi Markov-Bayesian cho thấy sự gia tăng rõ rệt của lũ thiết kế tại đập Soyang.

Tại Việt Nam, các nghiên cứu của Lê Đình Thành (1997), Phạm Ngọc Quý và cộng sự (2005), Vũ Thị Thu Lan & Hoàng Thanh Sơn (2012), Ngô Lê An và cộng sự (2015) trên hồ A Vương, hay Đoàn Thị Nội (2016) trên vùng núi Đông Bắc đã từng bước áp dụng lý thuyết mưa lũ thiết kế. Tuy nhiên, phần lớn các công trình chỉ sử dụng đơn lẻ một mô hình khí hậu hoặc dựa trên kịch bản phát thải cũ (SRES). Bằng cách khai thác đồng thời 11 mô hình GCM thế hệ CMIP5 dưới các kịch bản RCP4.5 và RCP8.5, luận án định vị vững chắc vai trò tiên phong trong việc cung cấp bức tranh định lượng chi tiết hóa đa mô hình cho toàn bộ khu vực Nam Trung Bộ.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Công trình mở rộng lý thuyết tính toán thủy văn truyền thống trong điều kiện chuỗi thủy văn phi dừng (non-stationarity) dưới tác động của biến đổi khí hậu. Luận án chứng minh rằng các tham số thống kê lịch sử ($C_v, C_s, \bar{X}$) không còn cố định mà biến động theo mức độ gia tăng bức xạ tác động ($4.5\text{ W/m}^2$ và $8.5\text{ W/m}^2$). Cụ thể:

  • Bổ sung cơ sở khoa học chuyển đổi từ kịch bản nồng độ khí nhà kính toàn cầu sang phân bố xác suất mưa cực trị địa phương tại các tiểu vùng có địa hình chắn gió phức tạp.
  • Phát triển mệnh đề khoa học: Tác động của BĐKH lên đỉnh lũ thiết kế ($Q_{\text{max}}$) có tính phi tuyến và khuếch đại mạnh hơn so với tỷ lệ biến động của lượng mưa ngày lớn nhất ($R_{x1day}$), do hiệu ứng bão hòa độ ẩm lưu vực và rút ngắn thời gian tập trung dòng chảy lũ.

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 3 trụ cột lý thuyết:

  1. Lý thuyết hoàn lưu khí quyển toàn cầu: Mô phỏng phản ứng bức xạ nhiệt của Trái Đất thông qua 11 mô hình GCM (Access1-3, CanESM2, CMCC-CMS, CNRM-CM5,...).
  2. Lý thuyết thống kê hiệu chỉnh phi tuyến: Áp dụng phương pháp khớp hàm phân phối xác suất tích lũy (Quantile Mapping) theo tiếp cận của Gudmundsson và cộng sự để loại bỏ sai số hệ thống giữa mô phỏng và thực tế đo đạc.
  3. Lý thuyết cân bằng ẩm và động lực dòng chảy lưu vực: Vận hành mô hình bể chứa bán phân bố NAM (Nielsen & Hansen, 1973) kết hợp phương pháp thủy lực chuyển lũ trên sông cho các lưu vực lớn (Vu Gia - Thu Bồn, Sông Kôn, Sông Ba) và phương pháp giải tích Sokolovsky/cường độ giới hạn cho các lưu vực nhỏ.

Điều kiện biên của khung phân tích được xác lập chặt chẽ: Áp dụng cho các lưu vực sông ngắn dốc thuộc vùng khí hậu nhiệt đới gió mùa Đông Trường Sơn, với chuỗi quan trắc khí tượng thủy văn đồng bộ sau năm 1975.

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

Thiết kế nghiên cứu

Nghiên cứu vận hành trên lập trường Thực chứng luận (Positivism), sử dụng phương pháp luận định lượng đa cấp độ (Multi-level design) kết hợp giữa mô hình hóa khí hậu, thống kê toán học và mô phỏng thủy văn:

  • Cấp độ toàn cầu/khu vực: 11 mô hình GCM thuộc CMIP5 với độ phân giải lưới khác nhau.
  • Cấp độ điểm trạm: 93 trạm quan trắc khí tượng phủ kín 8 tỉnh Nam Trung Bộ.
  • Cấp độ lưu vực: Phân cấp thành các lưu vực lớn ($A > 1000\text{ km}^2$) và các lưu vực vừa/nhỏ với đặc tính hình thái sông suối riêng biệt.

Quy trình nghiên cứu rigorous

Quy trình thực nghiệm được chuẩn hóa qua 5 bước nghiêm ngặt:

  1. Thu thập và kiểm soát chất lượng dữ liệu: Đồng bộ hóa chuỗi dữ liệu mưa ngày thực đo tại 93 trạm khí tượng thủy văn trong quá khứ (tập trung giai đoạn 1977–2010).
  2. Chi tiết hóa thống kê và hiệu chỉnh sai số: Trích xuất dữ liệu mưa ngày từ 11 GCMs cho thời kỳ nền và 2 giai đoạn tương lai (2040–2069 và 2070–2099) theo 2 kịch bản RCP4.5 và RCP8.5. Sai số trung bình và độ lệch chuẩn giữa tính toán mô hình và thực đo được hiệu chỉnh triệt để thông qua hàm phân phối xác suất thực nghiệm.
  3. Nội suy không gian: Áp dụng phương pháp nghịch đảo khoảng cách (IDW) trên nền tảng Hệ thống Thông tin Địa lý (GIS) để thành lập bản đồ biến động mưa $R_{x1day}$.
  4. Mô phỏng thủy văn lưu vực lớn: Thiết lập mô hình NAM cho các hệ thống lưu vực Nông Sơn, Thành Mỹ (Thu Bồn - Vu Gia), Bình Tường (Sông Kôn) và Củng Sơn (Sông Ba). Tiến hành hiệu chỉnh (calibration) và kiểm định (validation) thông số mô hình dựa trên chuỗi dòng chảy thực đo, tối ưu hóa hệ số tương quan Nash-Sutcliffe và sai số lưu lượng.
  5. Xây dựng tương quan $Q_{\text{max}} - Q_{\text{ngày max}}$ và công thức giải tích: Thiết lập quan hệ thực nghiệm giữa lưu lượng đỉnh lũ tức thời ($Q_{\text{max}}$) và lưu lượng bình quân ngày lớn nhất ($Q_{\text{ngày max}}$) tại các trạm kiểm soát, đồng thời áp dụng công thức Sokolovsky và cường độ giới hạn (TCVN 9845:2013, QP.TL C-6-77) cho lưu vực nhỏ.

Data và phân tích

Bộ số liệu thực nghiệm bao gồm 93 trạm đo mưa với chuỗi quan trắc liên tục từ 25 đến 36 năm. Phân tích thống kê xu thế lịch sử bằng phương trình hồi quy tuyến tính cho thấy sự phân hóa sâu sắc: Trong khi các trạm thuộc lưu vực Vu Gia - Thu Bồn có xu thế tăng rõ rệt (Tam Kỳ tăng $+2.65\text{ mm/năm}$, Giao Thủy $+1.75\text{ mm/năm}$, Nông Sơn $+1.46\text{ mm/năm}$, Cửa Lục $+1.27\text{ mm/năm}$), thì một số trạm phía Nam lại giảm (Hoài An $-2.53\text{ mm/năm}$, Vân Canh $-1.78\text{ mm/năm}$, Ayun Pa $-1.39\text{ mm/năm}$). Các công cụ phần mềm chuyên sâu được khai thác gồm gói mô hình MIKE 11 (module NAM), công cụ phân tích không gian ArcGIS và các thuật toán chi tiết hóa thống kê chuyên dụng.

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

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

Phân tích kết quả từ 11 mô hình GCM và mô phỏng thủy văn đã mang lại các phát hiện khoa học mang tính đột phá:

Hạng mục / Khu vực Giai đoạn 2040–2069 (RCP4.5 & RCP8.5) Giai đoạn 2070–2099 (RCP4.5 & RCP8.5) Bằng chứng thực nghiệm & Cơ chế tác động
Lưu vực Vu Gia - Thu Bồn (Trạm Nông Sơn, Thành Mỹ) Lượng mưa 1 ngày lớn nhất ($R_{x1day}$) tăng từ 10.5% đến 22.4% so với thời kỳ nền. $R_{x1day}$ tăng đột biến, đạt mức tăng từ 18.2% đến trên 35.6% dưới kịch bản RCP8.5. Đón gió mùa Đông Bắc và xoáy thuận nhiệt đới (XTNĐ) kết hợp địa hình núi cao (Ngọc Linh, Bà Nà).
Lưu vực Sông Kôn - Hà Thanh (Bình Định) $R_{x1day}$ tại vùng đồng bằng ven biển tăng 8.0% - 15.3%; vùng núi biến động nhẹ. $R_{x1day}$ toàn lưu vực tăng từ 12.1% đến 24.5%; tần suất mưa cực đoan gia tăng. Tác động trực tiếp từ bão đổ bộ; chuyển tiếp địa hình từ dãy Trường Sơn ra vịnh Quy Nhơn.
Lưu vực Sông Ba (Gia Lai, Phú Yên) Phân hóa mạnh: Thượng lưu (Kon Tum, Pleiku) tăng 5.2% - 12.0%; Hạ lưu (Củng Sơn) tăng thấp hoặc cục bộ giảm. Xu thế tăng mở rộng trên toàn lưu vực, trung bình đạt 10.4% - 21.8% ở giai đoạn cuối thế kỷ. Hiệu ứng phơn khuất gió và diện tích lưu vực rộng lớn ($13,900\text{ km}^2$) tạo độ trễ dòng chảy.
Đỉnh lũ thiết kế ($Q_{\text{max}}$) trạm Thành Mỹ Lưu lượng đỉnh lũ các tần suất thiết kế ($P = 0.1%, 1%, 5%$) tăng từ 12.5% đến 23.8%. Đỉnh lũ thiết kế $Q_{\text{max}}$ tăng vọt từ 22.0% đến 34.2% tương ứng với các kịch bản phát thải cao. Khuếch đại phi tuyến giữa mưa cực trị và dòng chảy mặt trên nền địa hình dốc sườn Đông.
Phân vùng nguy cơ dòng chảy lũ Xác lập vùng nguy cơ cao tập trung tại các tỉnh phía Bắc khu vực nghiên cứu (Đà Nẵng, Quảng Nam, Quảng Ngãi). Vùng nguy cơ rất cao mở rộng về phía Nam, bao trùm thêm phần lớn tỉnh Bình Định và Phú Yên. Bản đồ tích hợp GIS phân loại các cấp độ biến động dòng chảy lũ phục vụ an toàn công trình.

Các phân tích định lượng làm nổi bật những đặc trưng điển hình:

  1. Sự gia tăng áp đảo của mưa 1 ngày cực trị: Dưới kịch bản RCP8.5 giai đoạn 2070–2099, mức tăng lượng mưa 1 ngày lớn nhất đạt giá trị cực đại tại các tâm mưa vùng núi cao như Trà My (lượng mưa nhiều năm vốn đạt $3000 - 4000\text{ mm}$), làm thay đổi căn bản đường quá trình lũ thiết kế.
  2. Khuếch đại lưu lượng đỉnh lũ: Mức độ gia tăng của $Q_{\text{max}}$ cao hơn tốc độ gia tăng của lượng mưa từ 1.2 đến 1.5 lần, chứng minh phản ứng nhạy cảm của các lưu vực sông ngắn dốc miền Trung đối với các trận mưa cực đoan.
  3. Mâu thuẫn không gian sâu sắc: Nghiên cứu phát hiện hiện tượng trái quy luật tại một số trạm nội địa sườn Tây lưu vực sông Ba trong giai đoạn ngắn do hiệu ứng che chắn địa hình, khẳng định tính tất yếu của việc không sử dụng giá trị trung bình hóa toàn quốc.

Implications đa chiều

  • Đóng góp lý thuyết: Khẳng định tính phi dừng của các chuỗi số liệu khí tượng thủy văn dưới tác động của BĐKH, cung cấp cơ sở để tái cấu trúc lý thuyết tính toán an toàn công trình thủy lợi, giao thông.
  • Đổi mới phương pháp luận: Thiết lập quy trình chuẩn hóa kết hợp giữa chi tiết hóa thống kê 11 mô hình GCM với mô hình toán bán phân bố NAM và công thức kinh nghiệm, có khả năng chuyển giao áp dụng cho các vùng kinh tế sinh thái khác như Bắc Trung Bộ và Tây Nguyên.
  • Khuyến nghị thực tiễn & Chính sách:
    • Bộ Nông nghiệp & PTNT và Bộ Giao thông Vận tải: Cập nhật tiêu chuẩn TCVN 9845:2013 và quy phạm QP.TL C-6-77, bắt buộc bổ sung hệ số gia tăng lũ thiết kế từ 15% đến 35% khi tính toán khẩu độ cầu cống và dung tích phòng lũ hồ chứa tại Nam Trung Bộ.
    • Ban Chỉ huy PCTT & TKCN các tỉnh: Ứng dụng bản đồ phân vùng biến động dòng chảy lũ để rà soát quy trình vận hành liên hồ chứa trên lưu vực sông Vu Gia - Thu Bồn và sông Ba, chủ động phương án ứng phó sự cố vỡ đập tương tự bài học đập Khe Mơ (2010) hay đập Delhi (Hoa Kỳ, 2010).

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận các giới hạn học thuật:

  1. Giới hạn kịch bản sử dụng đất: Mô hình NAM trong nghiên cứu giả định thảm phủ thực vật và điều kiện sử dụng đất không thay đổi theo thời gian, chưa tích hợp quá trình đô thị hóa và suy thoái rừng đầu nguồn vào mô phỏng dòng chảy.
  2. Độ phân giải chuỗi số liệu thời gian: Dữ liệu chi tiết hóa từ GCMs dừng lại ở thời đoạn ngày ($R_{x1day}$), chưa đạt tới thời đoạn ngắn (mưa giờ, mưa 15 phút) để xây dựng đường cong IDF chi tiết cho đô thị hóa nhanh.
  3. Thế hệ mô hình khí hậu: Nghiên cứu sử dụng kịch bản RCPs thuộc CMIP5; các kịch bản SSPs (Shared Socioeconomic Pathways) mới nhất thuộc CMIP6 chưa được tích hợp tại thời điểm thực hiện luận án.

Chương trình nghiên cứu 10 năm tới cần tập trung vào 4 hướng chiến lược:

  1. Ứng dụng thế hệ mô hình CMIP6 kết hợp các kịch bản kinh tế - xã hội SSP1-2.6 đến SSP5-8.5 để tinh chỉnh dự báo cực trị khí hậu.
  2. Phát triển mô hình thủy văn - thủy lực liên kết 2 chiều (2D hydrodynamic) mô phỏng ngập lụt chi tiết vùng hạ du dưới tác động tổng hợp của lũ thượng nguồn, triều cường và nước biển dâng.
  3. Ứng dụng trí tuệ nhân tạo (Deep Learning/LSTM) để tự động hóa quá trình hiệu chỉnh tham số mô hình thủy văn phi tuyến.
  4. Đánh giá rủi ro vỡ đập dây chuyền (cascade failure) hệ thống bậc thang thủy điện trên lưu vực sông Vu Gia - Thu Bồn và sông Ba trong điều kiện lũ vượt thiết kế.

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

  • Học thuật: Mở ra hướng tiếp cận chuẩn mực trong việc tích hợp đa mô hình GCM vào tính toán thủy văn công trình tại Việt Nam, là tài liệu tham khảo nền tảng với tiềm năng trích dẫn cao trong các nghiên cứu biến đổi khí hậu khu vực Đông Nam Á.
  • Công nghiệp & Năng lượng: Cung cấp cơ sở khoa học trực tiếp cho Tập đoàn Điện lực Việt Nam (EVN) và các chủ đầu tư thủy điện nhằm đánh giá lại an toàn đập tràn xả lũ, quy hoạch dung tích đón lũ và tối ưu hóa hiệu quả phát điện.
  • Chính sách & Xã hội: Cung cấp công cụ hỗ trợ ra quyết định dựa trên bằng chứng (evidence-based policy) cho Ủy ban Quốc gia về Ứng phó Sự cố Thiên tai và Tìm kiếm Cứu nạn, góp phần bảo vệ tính mạng và tài sản cho hơn 9 triệu dân cư vùng duyên hải Nam Trung Bộ.

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

  • Nghiên cứu sinh & Giới học thuật: Tiếp cận phương pháp luận chi tiết hóa thống kê đa mô hình và quy trình hiệu chỉnh mô hình thủy văn bán phân bố.
  • Kỹ sư tư vấn thiết kế công trình Thủy lợi - Thủy điện - Giao thông: Sử dụng các thông số mưa, đỉnh lũ thiết kế cập nhật theo kịch bản BĐKH để xác định kích thước công trình an toàn, tránh nguy cơ xói lở, vỡ đập hoặc ngập lụt cục bộ.
  • Nhà hoạch định chính sách & Cơ quan quản lý tài nguyên nước: Khai thác bản đồ phân vùng biến động dòng chảy lũ để phân bổ nguồn vốn đầu tư công trình phòng chống thiên tai và quy hoạch không gian thoát lũ bền vững.

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 lý thuyết tính toán thủy văn truyền thống (Chow & Maidment, 1988) bằng việc thiết lập cơ chế chuyển đổi phi tuyến từ bức xạ tác động của khí nhà kính (RCP4.5, RCP8.5) sang lưu lượng đỉnh lũ thiết kế $Q_{\text{max}}$ thông qua việc phá vỡ giả định chuỗi dừng (stationarity), chứng minh hệ số khuếch đại dòng chảy đỉnh lũ tại các lưu vực ngắn dốc ven biển.

2. Đổi mới phương pháp luận thể hiện như thế nào khi so sánh với các nghiên cứu quốc tế?
So với nghiên cứu của Dankers & Feyen (2008) tại Châu Âu (dùng đơn lẻ HIRHAM-LISFLOOD) hay Hyun-Han Kwon và cộng sự (2011) tại Hàn Quốc (dùng mô hình ngẫu nhiên cho một hồ chứa đơn lẻ), luận án đã kết hợp đồng thời 11 mô hình GCMs với phương pháp chi tiết hóa thống kê hiệu chỉnh sai số phi tuyến trên mạng lưới 93 trạm dày đặc, tích hợp đa dạng phương pháp thủy văn (mô hình NAM cho lưu vực lớn và Sokolovsky cho lưu vực nhỏ) trên một không gian lãnh thổ rộng lớn 8 tỉnh.

3. Phát hiện bất ngờ nhất về mặt số liệu là gì?
Mặc dù xu thế chung là gia tăng mưa lớn, nghiên cứu phát hiện sự suy giảm cục bộ và phân hóa ngược chiều của mưa 1 ngày lớn nhất tại một số trạm nội địa lưu vực sông Ba (như Ayun Pa giảm $-1.39\text{ mm/năm}$, Hoài An $-2.53\text{ mm/năm}$ trong quá khứ) và mức tăng tương lai rất thấp so với lưu vực Vu Gia - Thu Bồn, chứng minh vai trò chi phối áp đảo của hiệu ứng phơn và cấu trúc sườn núi đối với phân bố mưa cực trị.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) không?
Quy trình nghiên cứu được mô tả chi tiết từng bước: từ thu thập số liệu trạm, trích xuất dữ liệu lưới GCMs, công thức hàm phân phối tích lũy hiệu chỉnh sai số, thiết lập thông số mô hình NAM, xây dựng đường tương quan $Q_{\text{max}} - Q_{\text{ngày max}}$ đến thuật toán nội suy IDW trên GIS, đảm bảo khả năng tái lập hoàn toàn trên các lưu vực khác.

5. Chương trình nghiên cứu 10 năm được định hình ra sao?
Chương trình nghiên cứu mở rộng tập trung vào: Tích hợp kịch bản phát thải thế hệ mới CMIP6 (SSPs); xây dựng mô hình liên kết thủy lực ngập lụt 2 chiều thời gian thực; và đánh giá rủi ro vỡ đập bậc thang dưới tác động kép của BĐKH và biến động thảm phủ lưu vực.

Kết luận

  1. Luận án đã giải quyết thành công bài toán chi tiết hóa lượng mưa cực trị từ 11 mô hình khí hậu toàn cầu (GCMs) thế hệ CMIP5 về 93 trạm quan trắc mặt đất trên toàn khu vực Nam Trung Bộ.
  2. Định lượng hóa xu thế gia tăng của lượng mưa 1 ngày lớn nhất ($R_{x1day}$) theo kịch bản RCP4.5 và RCP8.5, chỉ ra mức tăng cực đại vượt trên 35% tại các lưu vực đón gió phía Bắc trong giai đoạn 2070–2099.
  3. Chứng minh hiện tượng khuếch đại dòng chảy đỉnh lũ thiết kế ($Q_{\text{max}}$ tăng từ 12% đến trên 34% tại trạm Thành Mỹ), khẳng định sự cần thiết phải thay đổi tiêu chuẩn thiết kế công trình thoát nước và hồ chứa.
  4. Xây dựng thành công bộ bản đồ phân vùng biến động dòng chảy lũ thiết kế giai đoạn 2040–2069 và 2070–2099 trên nền tảng GIS, cung cấp công cụ trực quan phục vụ quy hoạch lãnh thổ.
  5. Tạo bước tiến đột phá về mặt phương pháp luận tính toán thủy văn có xét đến BĐKH, kết nối hoàn chỉnh chuỗi mô hình từ khí hậu toàn cầu $\rightarrow$ chi tiết hóa thống kê $\rightarrow$ mô hình toán thủy văn $\rightarrow$ tiêu chuẩn kỹ thuật công trình.
  6. Đặt nền móng vững chắc cho các nghiên cứu tiếp theo về an toàn đập, quản trị rủi ro thiên tai và thích ứng biến đổi khí hậu tại dải duyên hải miền Trung Việt Nam.