Tổng quan nghiên cứu

Thành phố Hà Nội có lượng mưa trung bình năm dao động trong khoảng từ 1400 mm đến 2000 mm, trong đó lượng mưa mùa hè chiếm từ 80% đến 90% tổng lượng mưa cả năm. Trong bối cảnh biến đổi khí hậu toàn cầu, các hiện tượng thời tiết cực đoan tại khu vực nội thành diễn biến ngày càng phức tạp, điển hình là trận mưa lịch sử tháng 11 năm 1984 với lượng mưa ngày đạt 394 mm và đợt mưa lớn cuối tháng 10 đầu tháng 11 năm 2008 gây ngập lụt nghiêm trọng tại 63 điểm nội đô. Vấn đề nghiên cứu trọng tâm đặt ra là các tiêu chuẩn thiết kế công trình thoát nước hiện hữu xây dựng từ nhiều thập kỷ trước không còn phản ánh chính xác quy luật mưa cực đoan đang biến đổi.

Mục tiêu cụ thể của luận văn bao gồm: nghiên cứu, ứng dụng phương pháp hiệu chỉnh sai số thống kê phù hợp cho chuỗi số liệu mưa ngày mô phỏng từ các mô hình khí hậu toàn cầu; xây dựng đường cong quan hệ Cường độ - Thời gian - Tần suất của mưa cho trạm khí tượng Láng trong giai đoạn hiện tại và giai đoạn tương lai; định lượng mức độ biến đổi của cường độ mưa cực đoan dưới tác động của biến đổi khí hậu. Phạm vi không gian của đề tài tập trung tại khu vực nội thành Hà Nội thông qua số liệu quan trắc trạm Láng. Phạm vi thời gian nghiên cứu bao gồm chuỗi số liệu lịch sử 1961 - 2005 (trong đó giai đoạn 1976 - 2005 làm thời kỳ cơ sở, 1961 - 1975 làm thời kỳ đánh giá độc lập) và giai đoạn tương lai 2070 - 2099 theo hai kịch bản phát thải khí nhà kính.

Nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc khi giúp giảm sai số mô phỏng mưa cực trị từ mức âm 70% ở mô hình thô xuống còn dưới âm 17% sau hiệu chỉnh, cung cấp cơ sở dữ liệu thủy văn định lượng chính xác phục vụ thiết kế công trình và quy hoạch thoát nước bền vững cho thủ đô.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng dựa trên sự kết hợp chặt chẽ giữa lý thuyết giá trị cực trị khí hậu, lý thuyết phân tích tần suất thủy văn và lý thuyết quy mô thống kê. Khung lý thuyết chủ đạo sử dụng hàm phân bố cực trị Gumbel để mô tả xác suất xuất hiện của các giá trị lượng mưa lớn nhất năm tương ứng với các chu kỳ lặp lại khác nhau. Bên cạnh đó, nghiên cứu kế thừa lý thuyết quy mô đơn giản của Menabde để mô hình hóa mối quan hệ lũy thừa giữa các mômen thống kê của cường độ mưa theo các thời đoạn khác nhau.

Các khái niệm cốt lõi được vận dụng xuyên suốt bao gồm:

  • Đường cong Cường độ - Thời gian - Tần suất: Đồ thị thể hiện mối quan hệ giữa cường độ mưa trung bình, khoảng thời gian duy trì mưa và chu kỳ lặp lại của trận mưa.
  • Phân vị và hàm phân bố xác suất tích lũy: Công cụ thống kê xác định ngưỡng giá trị mà tại đó một tỷ lệ phần trăm số liệu quan sát nằm dưới mức đó.
  • Hiệu chỉnh sai số thống kê: Quá trình toán học điều chỉnh hàm phân bố lý thuyết của số liệu mô hình khí hậu tiệm cận với phân bố thực tế của chuỗi quan trắc.
  • Hạ quy mô thời gian: Kỹ thuật toán học cho phép suy giảm quy mô thời gian từ số liệu mưa ngày xuống các thời đoạn ngắn dưới ngày phục vụ tính toán thủy lực.

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

Nguồn dữ liệu sử dụng trong luận văn gồm 45 năm số liệu quan trắc mưa ngày tại trạm Láng giai đoạn 1961 - 2005 và chuỗi số liệu quan trắc mưa thời đoạn ngắn từ 15 phút đến 720 phút trong 30 năm giai đoạn 1976 - 2005 do Viện Khoa học Khí tượng Thủy văn và Biến đổi khí hậu cung cấp. Dữ liệu kịch bản tương lai được trích xuất từ 5 mô hình khí hậu toàn cầu thuộc dự án CMIP5 bao gồm ACCESS1-0, GFDL-CM3, GFDL-ESM2G, MRI-CGCM3 và NorESM1-M với độ phân giải không gian từ 1,1 độ đến 2,5 độ kinh vĩ.

Phương pháp chọn mẫu trong nghiên cứu là phương pháp lấy mẫu cực đại hàng năm đối với từng thời đoạn mưa để đảm bảo tính độc lập và đồng nhất của chuỗi số liệu theo tiêu chuẩn khí tượng thủy văn. Cỡ mẫu phân tích đường cong thời kỳ cơ sở gồm 30 giá trị cực đại hàng năm cho mỗi thời đoạn từ năm 1976 đến năm 2005.

Quy trình phân tích được thực hiện theo ba giai đoạn:

  1. Nội suy song tuyến tính đưa dữ liệu dạng lưới của mô hình khí hậu về tọa độ trạm Láng và áp dụng phương pháp ánh xạ phân vị dựa trên hàm phân bố Gamma hai tham số để khử sai số hệ thống về tần suất ngày mưa và cường độ mưa.
  2. Kiểm định kỹ năng hiệu chỉnh trên hai tập mẫu độc lập gồm thời kỳ phụ thuộc 1976 - 2005 (30 năm) và thời kỳ độc lập 1961 - 1975 (15 năm).
  3. Ứng dụng hàm phân bố Gumbel kết hợp kỹ thuật quy mô đơn giản để thiết lập đường cong quan hệ cường độ mưa thời đoạn ngắn cho giai đoạn cơ sở 1976 - 2005 và giai đoạn tương lai 2070 - 2099.

Lý do lựa chọn phương pháp ánh xạ phân vị và hàm Gumbel là vì tính ưu việt đã được kiểm chứng trên toàn cầu trong việc bảo toàn tín hiệu biến đổi dài hạn của khí hậu, đồng thời khắc phục triệt để hiện tượng mô hình khí hậu thường xuyên đánh giá thấp lượng mưa cực trị thời đoạn ngắn.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Thứ nhất, kết quả đánh giá kỹ năng phương pháp hiệu chỉnh ánh xạ phân vị cho thấy sự cải thiện vượt bậc trong mô phỏng lượng mưa cực trị. Đối với các giá trị mưa lớn hơn phân vị 95% trong thời kỳ 1976 - 2005, mô hình thô chưa hiệu chỉnh có sai số trung bình rất lớn, dao động từ âm 68,0% ở mô hình GFDL-ESM2G đến âm 73,9% ở mô hình GFDL-CM3. Sau khi áp dụng hiệu chỉnh, sai số trung bình của các mô hình giảm xuống rõ rệt, chỉ còn dao động trong khoảng từ âm 6,0% đến âm 16,9%. Trong giai đoạn thử nghiệm độc lập 1961 - 1975, sai số lượng mưa cực trị của mô hình thô khoảng âm 67,0% đến âm 72,4% cũng được đưa về mức dương 12,8% đến 23,5%, khẳng định độ tin cậy cao của phương pháp.

Thứ hai, đường cong quan hệ cường độ mưa tại trạm Láng thời kỳ cơ sở 1976 - 2005 xác lập các ngưỡng giá trị định lượng quan trọng. Ứng với chu kỳ lặp lại 100 năm, cường độ mưa 15 phút đạt 206,0 mm/h (tương đương lượng mưa tích lũy khoảng 51,5 mm trong 15 phút), cường độ mưa 60 phút đạt 148,6 mm/h và cường độ mưa 24 giờ đạt 18,5 mm/h. Ở chu kỳ lặp lại 25 năm, cường độ mưa 15 phút là 175,2 mm/h và 60 phút là 121,4 mm/h (từng xuất hiện 1 lần trong trận mưa năm 1984). Ở chu kỳ lặp lại 2 năm, cường độ mưa 15 phút là 113,6 mm/h (xuất hiện 18 lần trong 30 năm quan trắc) và cường độ mưa 60 phút là 66,6 mm/h (xuất hiện 16 lần).

Thứ ba, dưới tác động của biến đổi khí hậu trong giai đoạn 2070 - 2099, cường độ mưa cực đoan có xu thế gia tăng rõ rệt ở cả hai kịch bản phát thải trung bình và phát thải cao. Mức gia tăng cường độ mưa thể hiện rõ nét nhất ở các thời đoạn ngắn dưới 60 phút so với quy mô ngày 24 giờ, làm gia tăng nguy cơ quá tải tức thời cho mạng lưới thoát nước đô thị.

Thảo luận kết quả

Sự gia tăng mạnh mẽ của cường độ mưa cực đoan trong tương lai bắt nguồn từ sự gia tăng nhiệt độ bề mặt và độ ẩm trong khí quyển theo phương trình Clausius-Clapeyron, khiến đối lưu nhiệt đới diễn ra mạnh hơn và tạo ra các ổ mây giông có cường suất mưa rất lớn. Kết quả nghiên cứu tại Hà Nội hoàn toàn tương đồng với các công bố quốc tế tại vùng Địa Trung Hải, Bắc Mỹ và Đông Nam Á khi chỉ ra rằng tốc độ tăng cường độ mưa ở quy mô dưới ngày luôn cao hơn quy mô ngày khoảng từ 4% đến 9%.

Về mặt trình bày dữ liệu, các phát hiện nghiên cứu được hệ thống hóa trực quan thông qua đồ thị so sánh phân vị với đường chuẩn một - một thể hiện sự hội tụ của dữ liệu sau hiệu chỉnh sát với đường quan trắc thực tế. Đồng thời, biểu đồ kết hợp cột và đường thể hiện biến trình năm của lượng mưa và số ngày mưa minh chứng sự triệt tiêu các ngày mưa phùn nhân tạo của mô hình thô. Các đường cong quan hệ cường độ mưa biểu diễn trên hệ trục logarit kép cho thấy tính song song hoàn hảo giữa các tần suất lặp lại từ 2 năm đến 100 năm, khẳng định tính đúng đắn của phương pháp hạ quy mô đơn giản.

Đề xuất và khuyến nghị

Thứ nhất, rà soát và cập nhật tiêu chuẩn thiết kế hệ thống cấp thoát nước đô thị: Bộ Xây dựng và Ủy ban Nhân dân Thành phố Hà Nội cần tiến hành điều chỉnh hệ số cường độ mưa tính toán trong quy chuẩn kỹ thuật quốc gia về công trình hạ tầng kỹ thuật đô thị, nâng mức tính toán lưu lượng thiết kế thêm từ 15% đến 25% cho các cống chính và trạm bơm tiêu úng trong giai đoạn 2026 - 2030.

Thứ hai, tích hợp kịch bản mưa cực đoan tương lai vào đồ án quy hoạch tổng thể: Sở Quy hoạch - Kiến trúc phối hợp cùng Viện Quy hoạch Xây dựng Hà Nội đưa bộ thông số đường cong quan hệ cường độ mưa giai đoạn 2070 - 2099 vào 100% các đồ án quy hoạch phân khu đô thị và quy hoạch tiêu thoát nước lưu vực sông Tô Lịch, sông Nhuệ và sông Đáy trước năm 2028.

Thứ ba, nâng cấp mạng lưới quan trắc khí tượng tự động độ phân giải cao: Trung tâm Dự báo Khí tượng Thủy văn Quốc gia cần đầu tư lắp đặt bổ sung ít nhất 30 trạm đo mưa tự động với tần suất truyền tin từ 5 phút đến 10 phút tại các quận nội thành trong giai đoạn 2026 - 2028 nhằm liên tục cập nhật chuỗi số liệu thời đoạn ngắn phục vụ hiệu chỉnh mô hình.

Thứ tư, xây dựng hệ thống hạ tầng trữ nước ngầm thông minh: Ban Quản lý Dự án Đầu tư Xây dựng Công trình Hạ tầng Kỹ thuật Thành phố Hà Nội cần triển khai xây dựng các bể ngầm điều tiết ngập úng với tổng dung tích thiết kế trên 500.000 m3 tại các khu vực trũng thấp có mật độ xây dựng cao trong giai đoạn 2026 - 2035 nhằm cắt giảm đỉnh lũ tức thời do các trận mưa cực đoan gây ra.

Đối tượng nên tham khảo luận văn

Nhóm 1: Các kỹ sư thiết kế hạ tầng kỹ thuật và chuyên gia thủy văn đô thị. Luận văn cung cấp bảng tra cứu thông số cường độ mưa chính xác theo từng thời đoạn từ 15 phút đến 24 giờ và từng chu kỳ lặp lại từ 2 năm đến 100 năm, giúp tối ưu hóa kích thước hình học của cống rãnh, kênh dẫn và công suất trạm bơm tiêu thoát nước.

Nhóm 2: Các nhà quản lý quy hoạch đô thị và cơ quan hoạch định chính sách biến đổi khí hậu. Tài liệu cung cấp cơ sở định lượng về mức độ gia tăng rủi ro ngập lụt đô thị trong 50 đến 80 năm tới, hỗ trợ việc phân vùng rủi ro thiên tai và phân bổ nguồn vốn đầu tư công hiệu quả cho các dự án ứng phó biến đổi khí hậu.

Nhóm 3: Các nhà khoa học, giảng viên và học viên cao học ngành Khí tượng, Thủy văn và Biến đổi khí hậu. Luận văn là tài liệu tham khảo chuẩn mực về quy trình toán học chi tiết trong việc hiệu chỉnh sai số phân vị và kỹ thuật hạ quy mô thời gian từ mô hình khí hậu toàn cầu về điểm trạm địa phương.

Nhóm 4: Các tổ chức tư vấn môi trường và đơn vị quản lý vận hành thoát nước đô thị. Nghiên cứu giúp các đơn vị như Công ty Thoát nước Hà Nội xây dựng các kịch bản diễn tập ứng phó ngập lụt cục bộ tương ứng với các ngưỡng mưa cực đoan có chu kỳ xuất hiện 10 năm, 25 năm và 50 năm.

Câu hỏi thường gặp

Đường cong quan hệ cường độ mưa có vai trò gì trong quản lý thoát nước đô thị?

Đường cong quan hệ cường độ mưa cung cấp mối liên hệ định lượng giữa tốc độ mưa trung bình, thời gian duy trì và xác suất lặp lại. Đây là thông số đầu vào cốt lõi để các kỹ sư tính toán lưu lượng nước mưa cực đại chảy tràn trên bề mặt đô thị, từ đó quyết định chính xác đường kính cống thoát nước và công suất máy bơm chống ngập úng cục bộ.

Vì sao phải hiệu chỉnh sai số số liệu mưa từ mô hình khí hậu toàn cầu?

Các mô hình khí hậu toàn cầu có độ phân giải thô từ 100 km đến 250 km nên không mô tả được địa hình cục bộ và mây đối lưu quy mô nhỏ. Điều này dẫn đến sai số hệ thống nghiêm trọng, chẳng hạn như đánh giá thấp lượng mưa cực trị tới khoảng 70% tại trạm Láng, do đó bắt buộc phải hiệu chỉnh thống kê trước khi ứng dụng thực tế.

Phương pháp ánh xạ phân vị hoạt động theo nguyên lý nào?

Phương pháp ánh xạ phân vị điều chỉnh hàm phân bố xác suất tích lũy của số liệu mô hình về trùng khớp với hàm phân bố tích lũy của số liệu quan trắc thực tế. Kỹ thuật này gồm hai bước: loại bỏ các ngày mưa nhỏ giả định dưới ngưỡng 0,1 mm để khớp tần suất số ngày mưa, sau đó dùng hàm Gamma hiệu chỉnh quy mô và hình dạng của giá trị cường độ mưa.

Cường độ mưa cực đoan chu kỳ 100 năm tại trạm Láng lớn như thế nào?

Theo kết quả phân tích trong giai đoạn 1976 - 2005, trận mưa có chu kỳ lặp lại 100 năm tại trạm Láng có cường độ 15 phút lên tới 206,0 mm/h (khoảng 51,5 mm trong 15 phút) và cường độ 1 giờ đạt 148,6 mm/h. Đây là những giá trị cực trị lịch sử chưa từng được quan trắc trọn vẹn trong chuỗi dữ liệu 30 năm của thế kỷ trước.

Phương pháp trong luận văn có thể áp dụng cho các tỉnh thành khác không?

Quy trình hiệu chỉnh ánh xạ phân vị kết hợp mô hình quy mô đơn giản hoàn toàn có thể nhân rộng cho các đô thị khác như Thành phố Hồ Chí Minh, Đà Nẵng hay Hải Phòng. Điều kiện tiên quyết là địa phương cần có chuỗi số liệu quan trắc mưa ngày liên tục trên 30 năm và tối thiểu một chuỗi mưa tự động thời đoạn ngắn để xác định tham số quy mô.

Kết luận

  • Luận văn đã làm chủ và áp dụng thành công phương pháp hiệu chỉnh sai số ánh xạ phân vị trên nền hàm phân bố Gamma, giúp giảm sai số mô phỏng mưa cực trị từ âm 73,9% xuống dưới âm 17%.
  • Thiết lập thành công bộ đường cong quan hệ cường độ mưa chuẩn xác cho trạm Láng giai đoạn 1976 - 2005 với 9 cấp thời đoạn và 6 chu kỳ lặp lại từ 2 năm đến 100 năm.
  • Định lượng xu thế biến đổi của mưa cực đoan trong giai đoạn 2070 - 2099, chứng minh tốc độ gia tăng cường độ mưa thời đoạn ngắn diễn ra gay gắt hơn so với mưa quy mô ngày.
  • Cung cấp bộ công cụ toán học và cơ sở dữ liệu khoa học tin cậy giúp các cơ quan quản lý đô thị Hà Nội cập nhật quy chuẩn thiết kế công trình thoát nước thích ứng biến đổi khí hậu.
  • Lộ trình giai đoạn 2026 - 2030 cần tập trung mở rộng ứng dụng phương pháp cho toàn bộ mạng lưới trạm quan trắc tự động trên địa bàn thành phố và số hóa tích hợp vào nền tảng quản lý ngập lụt thông minh.

Các đơn vị quản lý quy hoạch, cơ quan nghiên cứu và doanh nghiệp hạ tầng kỹ thuật nên chủ động khai thác ngay các kết quả nghiên cứu định lượng này để nâng cao năng lực chống chịu thiên tai cho hạ tầng thủ đô Hà Nội.