1. Mức 100 năm không phải lịch hẹn
Trong hồ sơ khí tượng hải văn, cụm từ “mức 100 năm” thường xuất hiện cạnh chiều cao sóng, tốc độ gió hoặc mực nước. Nếu đại lượng đang xét là cực đại từng năm và giả thiết thống kê không đổi, xác suất vượt mức đó trong một năm là p = 1/100. USGS giải thích cùng cách diễn đạt đối với lũ: chu kỳ lặp trung bình không có nghĩa các sự kiện cách nhau đúng số năm ấy. [1]
Cần ghi rõ biến được phân tích. Cực đại năm của chiều cao sóng có nghĩa Hₛ khác với chiều cao một con sóng riêng lẻ; mực nước tổng cộng khác phần nước dâng do bão. Một con số có nhãn 100 năm nhưng thiếu định nghĩa biến, khoảng lấy trung bình, mốc cao độ hoặc vị trí sẽ không đủ để đưa trực tiếp vào mô hình tải.
Bài này dùng toán xác suất và một phân bố giả định để minh họa thao tác. Không có dữ liệu đo của vùng biển Việt Nam trong ví dụ; các kết quả không phải thông số thiết kế cho một mỏ hay dự án cụ thể.
2. Xác suất ít nhất một lần vượt trong tuổi thọ
Giả sử mỗi năm độc lập và cùng xác suất vượt p. Xác suất không vượt trong L năm bằng (1 − p)ᴸ. Vì “có ít nhất một lần” là biến cố bổ sung:
Pᴸ = 1 − (1 − p)ᴸ
p = 1/T
T tính bằng năm theo định nghĩa từ chuỗi cực đại năm; L là số năm xét. Với T = 100 và L = 25, Pᴸ = 22,2179%. Nếu tuổi thọ chỉ 20 năm, giá trị là 18,2093%. Thời gian sử dụng ngắn hơn chu kỳ lặp không đồng nghĩa xác suất vượt gần bằng không.
Với cùng L = 25, mức T = 50 năm có xác suất vượt 39,6535%; mức T = 200 năm có xác suất 11,7780%. Phép gần đúng L/T cho lần lượt 50%, 25% và 12,5%, cao hơn kết quả chính xác. Khi Lp không nhỏ, không nên dùng phép gần đúng này làm kết quả cuối.
Nếu muốn xác suất ít nhất một lần vượt một ngưỡng chỉ bằng 10% trong 25 năm, phương trình ngược cho:
p = 1 − (1 − 0,10)^(1/25)
Ta được p = 0,00420555, tương ứng T = 237,781 năm. Đây là chuyển đổi xác suất, chưa phải yêu cầu phải thiết kế mọi công trình theo mức 238 năm. Mức độ chấp nhận, trạng thái giới hạn và hệ quả hư hỏng cần được xác định trong cơ sở thiết kế riêng.
3. Ví dụ giả định với phân bố Gumbel cực đại
Chọn X là cực đại năm của Hₛ tại một vị trí giả định. Đặt tham số vị trí μ = 6 m và tỷ lệ β = 1,2 m. Hàm phân phối tích lũy Gumbel cực đại có dạng:
F(x) = exp{−exp[−(x − μ)/β]}
NIST trình bày riêng phân bố Gumbel cực tiểu và cực đại; hai dạng khác dấu nên cần chọn đúng trước khi sử dụng hàm nghịch đảo. [2] Với cực đại năm, F(xT) = 1 − 1/T, suy ra:
xT = μ − β ln[−ln(1 − 1/T)]
Các logarit là logarit tự nhiên, đối số không thứ nguyên, β phải dương. Tại 50 năm, biến giảm zT = −ln[−ln(1 − 1/T)] bằng 3,901939 và xT = 10,6823 m. Tại 100 năm, zT = 4,600149, xT = 11,5202 m. Tại 200 năm, zT = 5,295812, xT = 12,3550 m.
Tăng chu kỳ lặp gấp đôi không làm chiều cao gấp đôi. Trong ví dụ, từ 100 lên 200 năm chỉ tăng 0,8348 m. Ngược lại, xác suất trong tuổi thọ thay đổi đáng kể như phần trước. Nên trình bày cả mức biến vật lý và xác suất, thay vì chỉ so các nhãn 50–100–200 năm.
Với mục tiêu xác suất 10% trong 25 năm vừa giải, xT = 12,5631 m. Có thể kiểm ngược bằng F(xT), sau đó nâng F(xT) lên lũy thừa 25; kết quả không vượt cả giai đoạn phải bằng 0,90. Phép kiểm khép kín này bắt được lỗi dùng p thay cho 1 − p trong hàm phân vị.
4. Độ nhạy tham số và độ không chắc chắn
Giữ β = 1,2 m nhưng tăng μ thêm 0,2 m làm tất cả xT tăng đúng 0,2 m. Giữ μ = 6 m, tăng β từ 1,2 lên 1,4 m làm x100 tăng từ 11,5202 lên 12,4402 m, chênh 0,9200 m. Độ nhạy với β là zT; tại chu kỳ lặp dài hơn, cùng một thay đổi β tác động lớn hơn.
Hai kịch bản này chỉ là khảo sát tham số, không phải khoảng tin cậy 95% hay 90%. Muốn có khoảng tin cậy phải xét phương pháp ước lượng, độ dài chuỗi, phân tán và sự phụ thuộc của các tham số. Không nên gọi hai đường “thấp–cao” tự chọn là biên thống kê đã được chứng minh.
Sản phẩm Sea Level Trends and Extremes của NOAA trình bày các mức xác suất vượt mực nước cùng khoảng tin cậy; phần giải thích cho biết sự kiện hiếm có ít điểm dữ liệu hơn và khoảng tin cậy rộng hơn. [3] Đây là ví dụ về cách báo cáo độ không chắc chắn. Nó không xác nhận phân bố Gumbel hai tham số đang chọn sẽ phù hợp với Hₛ tại mọi địa điểm.
Nếu chỉ có 30 cực đại năm, tăng chu kỳ lặp tính toán lên hàng trăm năm là ngoại suy đuôi phân bố. Cần khảo sát mức độ phù hợp của mô hình và ảnh hưởng của những năm lớn nhất, thay vì đánh giá chỉ bằng việc đường cong đi qua phần giữa dữ liệu. Thêm chữ số thập phân không làm phần ngoại suy chắc chắn hơn.
5. Khi xác suất thay đổi theo thời gian
Nếu các năm vẫn được giả thiết độc lập nhưng xác suất năm j là pⱼ khác nhau, công thức trở thành:
Pᴸ = 1 − ∏(1 − pⱼ)
Ví dụ giả định: 10 năm đầu p = 0,01, 15 năm sau p = 0,02 cho Pᴸ = 33,2051%. Giữ p = 0,01 cho toàn bộ 25 năm chỉ cho 22,2179%. Đây là minh họa ảnh hưởng tính không dừng, không phải dự báo xu hướng biển tại một khu vực.
Nếu có phụ thuộc giữa các năm, tích các xác suất riêng cũng không còn đủ. Tương tự, không thể lấy công thức xác suất năm để suy trực tiếp xác suất của một chiến dịch thi công sáu giờ bằng cách chia 8.760, nếu không có mô hình mùa, thời gian tồn tại sự kiện và tương quan phù hợp. Tải cực trị dài hạn và cửa sổ thời tiết ngắn hạn là hai câu hỏi thống kê khác nhau.
6. Không ghép máy móc các mức cực trị riêng
Giả sử gió và sóng mỗi biến có xác suất vượt ngưỡng 1% trong năm. Nếu tạm giả thiết hai biến độc lập, xác suất cả hai vượt trong cùng năm bằng 0,01%, còn ít nhất một biến vượt bằng 1,99%. Ngay cả “cùng năm” cũng không bảo đảm cùng giờ xảy ra. Nếu chúng phụ thuộc do cùng bão, các phép nhân độc lập này không đại diện tình huống thực.
Vì vậy tổ hợp sóng 100 năm, gió 100 năm và dòng 100 năm không tự có chu kỳ lặp 100 năm. Cần mô hình đồng thời của các biến và đáp ứng của công trình. Hồ sơ kết quả nên ghi biến, định nghĩa thời gian, dữ liệu và thời kỳ quan trắc, mô hình phân bố, mức vượt hằng năm, khoảng tin cậy và giả thiết tổ hợp. Nhờ những thông tin đó, người tính kết cấu mới biết con số cực trị đang đại diện điều gì.
Nguồn tham khảo
[1] Floods and Recurrence Intervals — U.S. Geological Survey. Truy cập ngày 09/11/2026.
[2] Engineering Statistics Handbook — Extreme Value Type I Distribution — NIST/SEMATECH. Truy cập ngày 09/11/2026.
[3] Sea Level Trends and Extremes — Main Product, Exceedance Probability Curves — NOAA CO-OPS. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Trong hồ sơ khí tượng hải văn, cụm từ “mức 100 năm” thường xuất hiện cạnh chiều cao sóng, tốc độ gió hoặc mực nước. Nếu đại lượng đang xét là cực đại từng năm và giả thiết thống kê không đổi, xác suất vượt mức đó trong một năm là p = 1/100. USGS giải thích cùng cách diễn đạt đối với lũ: chu kỳ lặp trung bình không có nghĩa các sự kiện cách nhau đúng số năm ấy. [1]
Cần ghi rõ biến được phân tích. Cực đại năm của chiều cao sóng có nghĩa Hₛ khác với chiều cao một con sóng riêng lẻ; mực nước tổng cộng khác phần nước dâng do bão. Một con số có nhãn 100 năm nhưng thiếu định nghĩa biến, khoảng lấy trung bình, mốc cao độ hoặc vị trí sẽ không đủ để đưa trực tiếp vào mô hình tải.
Bài này dùng toán xác suất và một phân bố giả định để minh họa thao tác. Không có dữ liệu đo của vùng biển Việt Nam trong ví dụ; các kết quả không phải thông số thiết kế cho một mỏ hay dự án cụ thể.
2. Xác suất ít nhất một lần vượt trong tuổi thọ
Giả sử mỗi năm độc lập và cùng xác suất vượt p. Xác suất không vượt trong L năm bằng (1 − p)ᴸ. Vì “có ít nhất một lần” là biến cố bổ sung:
Pᴸ = 1 − (1 − p)ᴸ
p = 1/T
T tính bằng năm theo định nghĩa từ chuỗi cực đại năm; L là số năm xét. Với T = 100 và L = 25, Pᴸ = 22,2179%. Nếu tuổi thọ chỉ 20 năm, giá trị là 18,2093%. Thời gian sử dụng ngắn hơn chu kỳ lặp không đồng nghĩa xác suất vượt gần bằng không.
Với cùng L = 25, mức T = 50 năm có xác suất vượt 39,6535%; mức T = 200 năm có xác suất 11,7780%. Phép gần đúng L/T cho lần lượt 50%, 25% và 12,5%, cao hơn kết quả chính xác. Khi Lp không nhỏ, không nên dùng phép gần đúng này làm kết quả cuối.
Nếu muốn xác suất ít nhất một lần vượt một ngưỡng chỉ bằng 10% trong 25 năm, phương trình ngược cho:
p = 1 − (1 − 0,10)^(1/25)
Ta được p = 0,00420555, tương ứng T = 237,781 năm. Đây là chuyển đổi xác suất, chưa phải yêu cầu phải thiết kế mọi công trình theo mức 238 năm. Mức độ chấp nhận, trạng thái giới hạn và hệ quả hư hỏng cần được xác định trong cơ sở thiết kế riêng.
3. Ví dụ giả định với phân bố Gumbel cực đại
Chọn X là cực đại năm của Hₛ tại một vị trí giả định. Đặt tham số vị trí μ = 6 m và tỷ lệ β = 1,2 m. Hàm phân phối tích lũy Gumbel cực đại có dạng:
F(x) = exp{−exp[−(x − μ)/β]}
NIST trình bày riêng phân bố Gumbel cực tiểu và cực đại; hai dạng khác dấu nên cần chọn đúng trước khi sử dụng hàm nghịch đảo. [2] Với cực đại năm, F(xT) = 1 − 1/T, suy ra:
xT = μ − β ln[−ln(1 − 1/T)]
Các logarit là logarit tự nhiên, đối số không thứ nguyên, β phải dương. Tại 50 năm, biến giảm zT = −ln[−ln(1 − 1/T)] bằng 3,901939 và xT = 10,6823 m. Tại 100 năm, zT = 4,600149, xT = 11,5202 m. Tại 200 năm, zT = 5,295812, xT = 12,3550 m.
Tăng chu kỳ lặp gấp đôi không làm chiều cao gấp đôi. Trong ví dụ, từ 100 lên 200 năm chỉ tăng 0,8348 m. Ngược lại, xác suất trong tuổi thọ thay đổi đáng kể như phần trước. Nên trình bày cả mức biến vật lý và xác suất, thay vì chỉ so các nhãn 50–100–200 năm.
Với mục tiêu xác suất 10% trong 25 năm vừa giải, xT = 12,5631 m. Có thể kiểm ngược bằng F(xT), sau đó nâng F(xT) lên lũy thừa 25; kết quả không vượt cả giai đoạn phải bằng 0,90. Phép kiểm khép kín này bắt được lỗi dùng p thay cho 1 − p trong hàm phân vị.
4. Độ nhạy tham số và độ không chắc chắn
Giữ β = 1,2 m nhưng tăng μ thêm 0,2 m làm tất cả xT tăng đúng 0,2 m. Giữ μ = 6 m, tăng β từ 1,2 lên 1,4 m làm x100 tăng từ 11,5202 lên 12,4402 m, chênh 0,9200 m. Độ nhạy với β là zT; tại chu kỳ lặp dài hơn, cùng một thay đổi β tác động lớn hơn.
Hai kịch bản này chỉ là khảo sát tham số, không phải khoảng tin cậy 95% hay 90%. Muốn có khoảng tin cậy phải xét phương pháp ước lượng, độ dài chuỗi, phân tán và sự phụ thuộc của các tham số. Không nên gọi hai đường “thấp–cao” tự chọn là biên thống kê đã được chứng minh.
Sản phẩm Sea Level Trends and Extremes của NOAA trình bày các mức xác suất vượt mực nước cùng khoảng tin cậy; phần giải thích cho biết sự kiện hiếm có ít điểm dữ liệu hơn và khoảng tin cậy rộng hơn. [3] Đây là ví dụ về cách báo cáo độ không chắc chắn. Nó không xác nhận phân bố Gumbel hai tham số đang chọn sẽ phù hợp với Hₛ tại mọi địa điểm.
Nếu chỉ có 30 cực đại năm, tăng chu kỳ lặp tính toán lên hàng trăm năm là ngoại suy đuôi phân bố. Cần khảo sát mức độ phù hợp của mô hình và ảnh hưởng của những năm lớn nhất, thay vì đánh giá chỉ bằng việc đường cong đi qua phần giữa dữ liệu. Thêm chữ số thập phân không làm phần ngoại suy chắc chắn hơn.
5. Khi xác suất thay đổi theo thời gian
Nếu các năm vẫn được giả thiết độc lập nhưng xác suất năm j là pⱼ khác nhau, công thức trở thành:
Pᴸ = 1 − ∏(1 − pⱼ)
Ví dụ giả định: 10 năm đầu p = 0,01, 15 năm sau p = 0,02 cho Pᴸ = 33,2051%. Giữ p = 0,01 cho toàn bộ 25 năm chỉ cho 22,2179%. Đây là minh họa ảnh hưởng tính không dừng, không phải dự báo xu hướng biển tại một khu vực.
Nếu có phụ thuộc giữa các năm, tích các xác suất riêng cũng không còn đủ. Tương tự, không thể lấy công thức xác suất năm để suy trực tiếp xác suất của một chiến dịch thi công sáu giờ bằng cách chia 8.760, nếu không có mô hình mùa, thời gian tồn tại sự kiện và tương quan phù hợp. Tải cực trị dài hạn và cửa sổ thời tiết ngắn hạn là hai câu hỏi thống kê khác nhau.
6. Không ghép máy móc các mức cực trị riêng
Giả sử gió và sóng mỗi biến có xác suất vượt ngưỡng 1% trong năm. Nếu tạm giả thiết hai biến độc lập, xác suất cả hai vượt trong cùng năm bằng 0,01%, còn ít nhất một biến vượt bằng 1,99%. Ngay cả “cùng năm” cũng không bảo đảm cùng giờ xảy ra. Nếu chúng phụ thuộc do cùng bão, các phép nhân độc lập này không đại diện tình huống thực.
Vì vậy tổ hợp sóng 100 năm, gió 100 năm và dòng 100 năm không tự có chu kỳ lặp 100 năm. Cần mô hình đồng thời của các biến và đáp ứng của công trình. Hồ sơ kết quả nên ghi biến, định nghĩa thời gian, dữ liệu và thời kỳ quan trắc, mô hình phân bố, mức vượt hằng năm, khoảng tin cậy và giả thiết tổ hợp. Nhờ những thông tin đó, người tính kết cấu mới biết con số cực trị đang đại diện điều gì.
Nguồn tham khảo
[1] Floods and Recurrence Intervals — U.S. Geological Survey. Truy cập ngày 09/11/2026.
[2] Engineering Statistics Handbook — Extreme Value Type I Distribution — NIST/SEMATECH. Truy cập ngày 09/11/2026.
[3] Sea Level Trends and Extremes — Main Product, Exceedance Probability Curves — NOAA CO-OPS. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
