Sử dụng số liệu mưa từ 1985–2014 tại 15 trạm: Vĩnh Yên, Sơn Tây, Láng, Hà Đông, Hưng Yên, Hà Nam, Nam Định, Văn Lý, Thái Bình, Ninh Bình, Nho Quan, Bắc Ninh, Hải Dương, Chí Linh, Hải Phòng.  

• Áp dụng phân phối Gumbel để tính tần suất lượng mưa lớn nhất năm theo các thời đoạn: 10 phút đến 24 giờ.  
• - Thiết lập quan hệ lượng mưa–thời gian–tần suất (DDF) bằng hàm số mũ, sau đó chuyển đổi thành công thức cường độ mưa:  
•   $$I_a(T) = a(T) \cdot d^{n(T)-1}$$  
• - Phân tích theo hai khoảng thời gian: d ≤ d* và d > d*, với điểm chuyển tiếp d* được xác định để đảm bảo tính liên tục.  
• - Tổng hợp công thức cường độ mưa cho 15 trạm ứng với T = 5, 10, 20 năm, với hệ số hồi quy R² > 0.99.  
• - So sánh với công thức TCVN 7957-2008 và Trần Việt Liễn cho thấy công thức mới có sai số thấp hơn rõ rệt (TB sai số: 2.6% so với 13.9% và 37.4%).  
• - Sử dụng ArcGIS 10 và phương pháp nội suy Spline để xây dựng bản đồ đẳng trị các tham số a, n của công thức IDF cho từng chu kỳ lặp lại.  
• - Bản đồ cho phép xác định cường độ mưa tại các vị trí không có trạm đo, phục vụ thiết kế hệ thống tiêu thoát nước đô thị, dân cư và công nghiệp.