Tổng quan luận án
Công trình khoa học "Nghiên cứu xác định một số tham số về mưa góp phần hoàn thiện công thức tính lưu lượng thiết kế công trình thoát nước nhỏ trên đường trong điều kiện khí hậu Việt Nam" được tác giả Nguyễn Anh Tuấn thực hiện tại Trường Đại học Giao thông Vận tải dưới sự hướng dẫn khoa học của PGS. Bùi Xuân Cậy và GS.TSKH Nguyễn Xuân Trục. Luận án thuộc chuyên ngành Xây dựng đường ôtô và đường thành phố (Mã số: 62.01), bảo vệ năm 2014.
Tính cấp thiết và khoảng trống nghiên cứu
Hệ thống công trình thoát nước nhỏ trên đường (bao gồm cầu nhỏ, cống và rãnh thoát nước mặt) chiếm tỷ trọng kinh phí xây dựng không lớn so với nền và mặt đường, nhưng năng lực thoát lũ lại mang tính chất quyết định đối với độ bền vững, tuổi thọ khai thác và an toàn giao thông. Những hạn chế trong khâu tính toán lưu lượng thiết kế ($Q_p$) có thể dẫn đến các sự cố nghiêm trọng:
- Hiện tượng khẩu độ công trình không đủ tiêu thoát gây ngập úng thượng lưu, phá hủy nền mặt đường, làm gián đoạn giao thông và gây tổn thất kinh tế - xã hội (điển hình như các đợt ngập úng tại miền Trung hoặc ngập lụt đô thị).
- Hiện tượng tạo dòng chảy xiết và lũ quét nguy hiểm ở hạ lưu do thế năng chuyển thành động năng (điển hình như trận lũ quét lịch sử tại Sơn Diệm năm 2002 trên tuyến đường Hồ Chí Minh).
- Sự biến đổi khí hậu toàn cầu và hiện tượng nước biển dâng làm gia tăng các hình thế thời tiết cực đoan, thay đổi chế độ mưa tại Việt Nam. Trong khi đó, các cơ sở dữ liệu và tham số khí tượng - thủy văn phục vụ thiết kế công trình thoát nước trên đường ở nước ta chủ yếu được xác lập từ năm 1987 trở về trước (như trong tiêu chuẩn TCVN 9845:2013), khiến độ tin cậy của kết quả tính toán theo thời gian bị suy giảm.
Mục tiêu nghiên cứu
Luận án xác định hai mục tiêu chính:
- Mục tiêu lý thuyết: Nghiên cứu bản chất thủy văn, đề xuất nguyên tắc, phương pháp và hoàn thiện hệ thống công thức xác định các tham số về mưa (lượng mưa ngày tính toán, hệ số đặc trưng hình dạng cơn mưa, cường độ mưa tính toán và các thông số liên quan) phục vụ tính toán lưu lượng thiết kế công trình thoát nước nhỏ trên đường trong điều kiện khí hậu Việt Nam.
- Mục tiêu thực tiễn: Ứng dụng chuỗi số liệu quan trắc thực tế giai đoạn 1960 - 2010 để xác lập giá trị cụ thể của các tham số mưa tại 12 trạm khí tượng đại diện, cung cấp số liệu tham khảo cho công tác tư vấn thiết kế và quản lý công trình hạ tầng giao thông.
Đối tượng và phạm vi nghiên cứu
- Đối tượng nghiên cứu: Công trình thoát nước nhỏ (cầu nhỏ, cống thoát nước ngang, rãnh thoát nước mặt) trên hệ thống đường bộ, đường sắt, đường đô thị và sân bay.
- Phạm vi nghiên cứu: Áp dụng cho các lưu vực nhỏ (thường có diện tích $F \le 100\text{ km}^2$, đặc biệt là các lưu vực từ dưới $1\text{ km}^2$ đến một vài $\text{km}^2$, hoặc $F \le 30\text{ km}^2$) trên phạm vi lãnh thổ Việt Nam, khảo sát chuyên sâu tại 12 trạm khí tượng điển hình thuộc các vùng khí hậu khác nhau.
Tổng quan tài liệu và vị trí của luận án
Các hướng nghiên cứu và công thức tính toán trên thế giới
Lý thuyết nền tảng cho việc xác định lưu lượng đỉnh lũ thiết kế lưu vực nhỏ bắt nguồn từ Lý thuyết tập trung nước từ lưu vực do N. Prôtôđiakônov (Liên Xô cũ) khởi xướng, làm cơ sở hình thành công thức căn nguyên dòng chảy và công thức cường độ giới hạn. Luận án đã hệ thống hóa các mô hình tính toán quốc tế tiêu biểu:
- Tại Anh: Sử dụng công thức cường độ giới hạn $Q = K \cdot C \cdot I \cdot A$, trong đó lưu lượng phụ thuộc hệ số dòng chảy $C$, cường độ mưa trung bình $I$ ứng với thời gian tập trung nước và diện tích $A$.
- Tại Mỹ: Áp dụng công thức $Q_p = K \cdot \alpha \cdot C \cdot I \cdot A = \alpha \cdot C \cdot q \cdot A$, bổ sung hệ số phân bố mưa rào $\alpha$ dạng đồ thị tra cứu.
- Tại Pháp: Mô hình Albert Caquot dựa trên nguyên lý hiệu quả năng lực mạng lưới và cân bằng thể tích: $Q = \alpha \cdot C \cdot H / t_c^\beta \dots$ chuyển đổi trực tiếp lượng mưa thời đoạn thành cường độ mưa tính toán $a_{t_c}$.
- Tại Nga và các nước SNG: Công thức của Viện Thiết kế đường Liên Xô cũ (B.I. Rôtpurt) xét tới hệ số hình dạng lưu vực $K_d$, độ dốc $K_i$ và hệ số triết giảm diện tích $\delta$; công thức theo quy trình CH435-72 của Viện Thủy văn Nhà nước Liên Xô tính toán qua mô đuyn đỉnh lũ tương đối $A_{1%}$ và lượng mưa ngày $H_{n,1%}$.
- Tại Trung Quốc: Viện Nghiên cứu Thủy văn đưa ra công thức suy giảm lượng mưa thời đoạn ngắn $H_{T,p} = H_{1h,p} \cdot T^{1-m}$ dựa trên chỉ số triết giảm mưa rào $m$.
Các công trình và tiêu chuẩn tính toán tại Việt Nam
- TCVN 9845:2013 (Tính toán các đặc trưng dòng chảy lũ): Sử dụng công thức cường độ giới hạn chuyển đổi từ lượng mưa ngày tính toán $Q_p = A_p \cdot \alpha \cdot H_{n,p} \cdot \delta_1 \cdot F$, giải lặp thông qua mô đuyn dòng chảy tương đối $A_p = \frac{16.67 \cdot \phi_\tau}{\tau}$ và các đặc trưng địa mạo thủy văn sườn dốc ($\Phi_{sd}$), lòng sông ($\Phi_{ls}$).
- Công thức Đại học Xây dựng Hà Nội (GS.TSKH Nguyễn Xuân Trục đề xuất từ thập niên 1980 cho $F \le 30\text{ km}^2$, đưa vào Sổ tay tính toán thủy văn của Bộ GTVT): Tính trực tiếp qua cường độ mưa tính toán $Q_p = 16.67 \cdot \alpha \cdot a_{\tau,p} \cdot \delta \cdot \delta_1 \cdot F$.
- TCVN 7957:2008 (Thoát nước - Mạng lưới và công trình bên ngoài): Dùng mô đuyn dòng chảy mưa $q_{\psi,p}$ cho mạng lưới thoát nước đô thị.
- Công thức Sôkôlôpsky (Quy trình Khảo sát và thiết kế cầu vượt sông 1972) và Phương pháp cân bằng lượng nước (Quy phạm BCH 63-61 Liên Xô): Dùng cho lưu vực vừa và lớn ($F > 100\text{ km}^2$) hoặc bài toán xác định đường quá trình lũ chi tiết.
Lịch sử nghiên cứu phân vùng mưa tại Việt Nam
- Quy phạm thủy lợi QP.C-6-77 (năm 1977): Phân chia miền Bắc thành 10 khu mưa rào.
- Nghiên cứu của Hoàng Minh Tuyển (1980): Sử dụng số liệu 60 trạm khí tượng (chuỗi 10 - 20 năm), chia cả nước thành 15 vùng mưa.
- Nghiên cứu của Hoàng Niêm và Đỗ Đình Khôi (1991): Phân chia 18 vùng mưa với 18 đường cong $\phi_T \sim T$.
- Nghiên cứu của TS. Trịnh Nhân Sâm (1993): Hiệu chỉnh 18 vùng mưa và đường cong $\phi_T$, làm cơ sở cho quy chuẩn ngành và TCVN 9845:2013 (chuỗi số liệu kết thúc năm 1987).
Khoảng trống nghiên cứu luận án tập trung giải quyết
- Chuỗi số liệu mưa trước đây đã quá cũ (hơn 25 năm chưa cập nhật), chưa phản ánh được tác động của biến đổi khí hậu và các trận mưa cực đoan sau năm 1987.
- Việc phân chia 18 vùng mưa hiện hành còn quá rộng đối với các lưu vực nhỏ.
- Số lượng trạm có máy đo mưa tự ghi còn hạn chế (khoảng 20% tổng số điểm đo), gây khó khăn cho việc xác định trực tiếp cường độ mưa tính toán $a_{T,p}$, đòi hỏi phải hoàn thiện các công thức thực nghiệm gián tiếp có độ chính xác cao.
Cơ sở lý thuyết và phương pháp nghiên cứu
Cơ sở lý thuyết
- Lý thuyết tập trung nước từ lưu vực: Dựa trên các giả thiết của M. Prôtôđiakônov về sự phân bố đồng đều của mưa và tổn thất, quan niệm lưu lượng đỉnh lũ đạt cực đại ($Q_{max}$) khi giọt nước xa nhất trên lưu vực chảy về tới mặt cắt công trình ứng với thời gian tập trung nước $\tau$.
- Phương trình cân bằng nước: Áp dụng định luật bảo toàn vật chất $W = W_{sd} + W_{ls} + W_Q$ nhằm liên kết các thành phần dòng chảy sườn dốc, dòng chảy lòng dẫn và lưu lượng xả qua công trình.
- Lý thuyết phân tích thống kê xác suất thủy văn: Xử lý chuỗi mẫu cực trị hàng năm ($H_{ngàymax}$, $a_{Tmax}$), giải bài toán phân phối xác suất và đường tần suất lý luận (Pearson III hoặc Kritski-Menkel).
[Tài liệu mưa quan trắc thực tế (1960 - 2010)]
|
+-----------------------+-----------------------+
| |
[Đo lượng mưa ngày] [Đo mưa tự ghi]
| |
(Kiểm định Smirnov-Kolmogorov) (Đường cong a - T - p trực tiếp)
| |
Tham số H_{n,p} Tham số a_{T,p}
\ /
\ /
+---> [7 Công thức thực nghiệm tính a_{T,p}] <---+
- Sức mưa Sp, hệ số hình dạng m
- Hệ số khí hậu A, B, hệ số beta
- Cường độ chuẩn a_{T0,p}, trạm tựa
|
[Lưu lượng thiết kế công trình nhỏ Q_p]
Phương pháp nghiên cứu thực tế
- Phương pháp thống kê xác suất: Thu thập và chỉnh lý chuỗi số liệu đo mưa ngày và đo mưa tự ghi liên tục 51 năm (1960 - 2010). Kiểm định tính độc lập, tính đại diện và tính đồng nhất của mẫu số liệu; kiểm định sự phù hợp của đường tần suất lý luận bằng tiêu chuẩn Smirnov - Kolmogorov ở mức ý nghĩa $\alpha = 5%$.
- Phương pháp hồi quy phi tuyến và đa biến: Thiết lập tương quan giữa các đặc trưng mưa (sức mưa $S_p$, hệ số hình dạng $m$, hệ số khí hậu $A, B$, hệ số hồi quy $\beta$) với hệ số xác định $R^2$.
- Nguồn dữ liệu thực nghiệm: Chuỗi số liệu quan trắc khí tượng thủy văn tại 12 trạm đại diện do Trung tâm Khí tượng Thủy văn Quốc gia quản lý, cùng số liệu rà soát năm 2008 của 54 trạm khí tượng trên toàn quốc.
Nội dung chính theo từng chương
Danh mục bảng biểu và cấu trúc dữ liệu chính trong luận án
| Số hiệu bảng |
Nội dung và tiêu đề bảng theo văn bản |
Mục tiêu khoa học |
| Bảng 2.1 |
Thông tin số liệu đo mưa ở 12 trạm khí tượng chọn nghiên cứu (1960 - 2010) |
Xác định cơ sở dữ liệu quan trắc |
| Bảng 2.2 |
Tổng hợp tháng mưa nhiều ngày và so sánh với tháng mùa mưa |
Đánh giá quy luật phân bố mùa |
| Bảng 2.3 |
Xu hướng biến thiên lượng mưa năm và số ngày mưa trong năm |
Nhận diện tác động của biến đổi khí hậu |
| Bảng 2.4 |
Xu hướng biến thiên của $H_{ngàymax}$ và $a_{Tmax}$ ($T = 5\text{ ph} \div 1440\text{ ph}$) |
Phân tích biến động cực đoan thời đoạn |
| Bảng 2.5 |
So sánh giá trị $(H_{ngàymax})^$, $(a_{Tmax})^$ đột biến với $H_{n,p}, a_{T,p}$ ($p = 4%, 1%$) |
Xác định tính bất thường của mưa cực hạn |
| Bảng 2.6 - 2.7 |
Giá trị $C_v, C_s$ của $H_{ngàymax}$ và $a_{Tmax}$ tại 12 trạm nghiên cứu |
Phân tích tham số thống kê thủy văn |
| Bảng 2.8 |
Chu kỳ biến đổi lớn - nhỏ - trung bình của $H_{ngàymax}$ và $a_{Tmax}$ |
Xác định chu kỳ lặp thời tiết |
| Bảng 2.9 - 2.10 |
Mức độ trùng lặp thời điểm xuất hiện $H_{ngàymax}$ và $a_{Tmax}$ ($5\text{ ph} \div 1440\text{ ph}$) |
Kiểm tra quan hệ đồng thời giữa lượng mưa ngày và cường độ mưa |
| Bảng 3.1 - 3.3 |
Số năm quan trắc cần thiết, lượng mưa đột biến và kiểm định Smirnov - Kolmogorov |
Đánh giá độ tin cậy của chuỗi $H_{n,p}$ |
| Bảng 3.4 - 3.6 |
Quá trình xác định và đánh giá sai số hệ số đặc trưng hình dạng cơn mưa $\phi_T$ |
Chuẩn hóa đường cong $\phi_T \sim T$ |
| Bảng 4.1 - 4.8 |
Trình tự hồi quy và hệ số tương quan $R^2$ xác định $m, S_p, A, B, \beta$ |
Thiết lập các công thức thực nghiệm $a_{T,p}$ |
| Bảng 4.9 |
Đánh giá và so sánh mức độ sai số của các công thức tính $a_{T,p}$ |
Kiểm chứng độ chính xác mô hình |
Chương 1: Tổng quan về vấn đề nghiên cứu
Chương này tập trung phân tích cơ chế hình thành dòng chảy lũ do mưa trên lưu vực nhỏ, chỉ ra ba giai đoạn chính: giai đoạn tổn thất hoàn toàn ($a_t \le g_t$), giai đoạn sinh dòng chảy ($a_t > g_t$) và giai đoạn sau khi kết thúc sinh dòng chảy. Luận án phân tích chi tiết bản chất của cường độ mưa hiệu quả $a_{hqt} = a_t - g_t$ trong thời gian cung cấp nước $T_{cn}$.
Cường độ
^
| /-----\ Đường quá trình mưa a(t)
| / | \
| ...../....|....\..... Đường tổn thất g(t)
| : / | \ :
| : / a_hqt \ :
| : / | \ :
--+--+-/--------+--------\-+---> Thời gian (t)
0 | <------ Tcn ------> |
|<------- T_mưa ----->|
Phân tích thứ tự ảnh hưởng của các nhân tố đến dòng chảy lũ lưu vực nhỏ:
- Nhân tố mưa: Mang tính chất quyết định và duy nhất cung cấp nguồn dòng chảy tại Việt Nam; lượng mưa mùa mưa chiếm khoảng 80% cả năm, tương ứng dòng chảy mùa lũ chiếm khoảng 70%.
- Nhân tố mặt đệm: Địa hình (độ cao, hướng sườn đón gió ẩm), địa mạo (lớp phủ, độ nhám sườn dốc $n_{sd}$, độ nhám lòng suối $n_{ls}$, tỷ lệ ao hồ đầm lầy) quyết định thời gian tập trung nước $\tau$ và tổn thất thấm $i$.
- Tần suất thiết kế ($p$): Quy định mức độ an toàn công trình theo cấp đường ($p = 1% \div 4%$ theo TCVN 4054:2005 đối với cầu nhỏ, cống, rãnh đường ô tô; chu kỳ $N = 0.33 \div 10\text{ năm}$ theo TCVN 7957:2008 đối với thoát nước đô thị).
Chương 1 khẳng định tham số cường độ mưa tính toán $a_{\tau,p} = H_{\tau,p}/\tau$ là tham số mang tính chất tổng hợp cao nhất, phản ánh đồng thời chế độ mưa, hình dạng cơn mưa, đặc trưng mặt đệm và tần suất thiết kế.
Chương 2: Nghiên cứu đặc điểm mưa chịu tác động của hiện tượng biến đổi khí hậu trong tính toán lưu lượng thiết kế công trình thoát nước nhỏ trên đường
Chương 2 tổng kết đặc trưng khí hậu nhiệt đới gió mùa ẩm của Việt Nam:
- Độ ẩm không khí trung bình năm trên 80% (77% - 87%).
- Lượng mưa hàng năm dao động từ 1100 mm đến 8000 mm; số ngày mưa từ 67 đến 223 ngày/năm.
- Sự phân hóa rõ rệt giữa hai miền: Miền Bắc (nhiệt độ trung bình $< 24^\circ\text{C}$, giờ nắng $< 2000\text{ h}$, tổng xạ $< 586\text{ KJ/cm}^2$), Miền Nam (nhiệt độ trung bình $24 - 28^\circ\text{C}$, giờ nắng $> 2000\text{ h}$, tổng xạ $> 586\text{ KJ/cm}^2$).
- Các trung tâm mưa lớn (Đông Triều, Bắc Quang, Hoàng Liên Sơn, Tam Đảo, Bà Nà, Trà Mi, Bạch Mã với lượng mưa trên 8000 mm/năm; dải "vĩ tuyến nước" $15^\circ\text{B} - 16^\circ\text{B}$) và trung tâm mưa nhỏ (An Châu, Sơn La, Mường Xén, Ninh Thuận khoảng 800 - 1000 mm/năm; dải khô hạn duyên hải cực Nam Trung Bộ $10^\circ\text{B} - 12^\circ\text{B}$).
Chương này đánh giá mạng lưới trạm quan trắc (toàn quốc có khoảng 172 trạm khí tượng, 231 trạm thủy văn, trên 400 điểm đo mưa nhân dân, đạt mật độ bình quân $433\text{ km}^2/\text{điểm}$, trong đó khoảng 20% số điểm có máy đo tự ghi) và khảo sát chuỗi số liệu 1960 - 2010 tại 12 trạm nghiên cứu:
- Danh sách 12 trạm khí tượng: TX. Điện Biên, TP. Tuyên Quang, TP. Lạng Sơn, Trạm Láng (Hà Nội), Trạm Hà Đông (Hà Nội), TX. Sơn Tây (Hà Nội), TP. Nghệ An, TP. Quảng Bình, TP. Khánh Hòa, TP. Buôn Ma Thuột (Đắk Lắk), TP. Cần Thơ.
- Kết quả phân tích:
- Nhận diện xu thế biến thiên của lượng mưa năm, số ngày mưa, lượng mưa ngày lớn nhất năm ($H_{ngàymax}$) và cường độ mưa lớn nhất năm ($a_{Tmax}$) ở các thời đoạn ngắn ($T = 5\text{ ph} \div 1440\text{ ph}$).
- Tính toán các đặc trưng thống kê: Giá trị trung bình nhiều năm $\bar{X}$, hệ số biến sai $C_v$, hệ số thiên lệch $C_s$.
- Phân tích hiện tượng xuất hiện các trận mưa lớn đột biến $(H_{ngàymax})^, (a_{Tmax})^$ vượt qua các mức tần suất thiết kế thông thường ($p = 1%, 4%$).
- Đánh giá mức độ trùng lặp thời điểm xuất hiện giữa $H_{ngàymax}$ và $a_{Tmax}$, chỉ ra tính độc lập tương đối giữa cường độ mưa thời đoạn ngắn và tổng lượng mưa ngày tại các lưu vực nhỏ.
Chương 3: Xác định lượng mưa ngày tính toán và nghiên cứu xác định hệ số đặc trưng hình dạng cơn mưa
Nội dung chương 3 giải quyết hai bài toán then chốt:
-
Xác định lượng mưa ngày tính toán ($H_{n,p}$):
- Thực hiện lấy mẫu thống kê và kiểm định số năm quan trắc cần thiết để đảm bảo sai số lấy mẫu cho phép.
- Xử lý các trường hợp mẫu có một hoặc nhiều cơn mưa đặc biệt lớn (ĐBL).
- Kiểm định sự phù hợp của đường tần suất lý luận $H_{n,p}$ với tài liệu thực đo tại 12 trạm bằng tiêu chuẩn Smirnov - Kolmogorov ở mức ý nghĩa $\alpha = 5%$.
- So sánh giá trị $H_{n,p}$ giai đoạn 1960 - 2010 với số liệu trong TCVN 9845:2013 (tính đến năm 1987), chỉ ra sự thay đổi rõ rệt của giá trị mưa tính toán dưới tác động của biến đổi khí hậu.
-
Xác định hệ số đặc trưng hình dạng cơn mưa ($\phi_T$) và tiêu chí phân vùng:
- Xây dựng hàm quan hệ $\phi_T \sim T$ cho 12 trạm nghiên cứu từ chuỗi số liệu 1960 - 2010.
- Đánh giá sai số của đường cong $\phi_T$ trung bình vùng so với các đường $\phi_{T,p}$ ở các mức tần suất khác nhau.
- Khảo sát thực tế tại 3 trạm lân cận thuộc khu vực Hà Nội (Láng, Hà Đông, Sơn Tây) để chứng minh tính đồng nhất cục bộ và sự khác biệt liên vùng.
- Đề xuất phương pháp và tiêu chí sử dụng hệ số đặc trưng hình dạng cơn mưa $\phi_T$ làm công cụ phân vùng mưa lũ tỷ lệ lớn phục vụ tính toán lưu vực nhỏ.
Chương 4: Nghiên cứu xác định tham số cường độ mưa trong tính toán lưu lượng thiết kế công trình thoát nước nhỏ trên đường ở Việt Nam
Chương 4 hoàn thiện phương pháp luận và công thức thực nghiệm xác định cường độ mưa tính toán $a_{T,p}$:
- Phương pháp trực tiếp:
- Xây dựng hệ đường cong $a - T - p$ trực tiếp từ số liệu máy đo mưa tự ghi (xử lý cho cả trường hợp chuỗi số liệu liên tục và trường hợp gián đoạn quan trắc).
- Hệ thống 7 công thức thực nghiệm tính gián tiếp $a_{T,p}$:
- Dạng 1: Tính $a_{T,p}$ dựa vào lượng mưa ngày tính toán $H_{n,p}$ và hệ số hình dạng cơn mưa $\phi_T$.
- Dạng 2: Tính $a_{T,p}$ theo đặc trưng sức mưa $S_p$ và hệ số hình dạng cơn mưa $m$: $a_{T,p} = f(S_p, m, T)$.
- Dạng 3: Tính $a_{T,p}$ theo hệ số vùng khí hậu $A, B$ và hệ số hình dạng cơn mưa $m$.
- Dạng 4: Tính $a_{T,p}$ theo hệ số hồi quy vùng khí hậu $\beta$, hệ số $m$ và lượng mưa ngày tính toán $H_{n,p}$.
- Dạng 5: Tính $a_{T,p}$ theo cường độ mưa chuẩn $a_{T0,p}$.
- Dạng 6: Phương pháp sử dụng trạm tựa nội suy theo lượng mưa ngày tính toán $H_{n,p}$.
- Dạng 7: Phương pháp sử dụng trạm tựa nội suy theo đặc trưng sức mưa $S_p$.
- Đánh giá sai số và điều kiện áp dụng:
- Thực hiện hồi quy xác định các tham số $m, S_p, A, B, \beta$ với hệ số tương quan $R^2$ cao tại 12 trạm nghiên cứu.
- So sánh mức độ sai số của từng công thức thực nghiệm so với đường cong thực đo trực tiếp $a - T - p$, phân loại điều kiện áp dụng phù hợp với từng mức độ sẵn có của tài liệu khí tượng tại địa phương thiết kế.
Kết quả và những đóng góp mới
Những đóng góp mới về mặt khoa học và thực tiễn
- Cập nhật chuỗi số liệu lượng mưa ngày tính toán ($H_{n,p}$): Xác lập bảng giá trị $H_{n,p}$ ở các cấp tần suất thiết kế cho 12 trạm khí tượng đại diện với chuỗi số liệu quan trắc 51 năm (1960 - 2010), phản ánh đầy đủ diễn biến thời tiết trong giai đoạn có biến đổi khí hậu, nâng cao độ tin cậy khi tính toán theo TCVN 9845:2013 hoặc công thức Sôkôlôpsky.
- Chuẩn hóa hệ số đặc trưng hình dạng cơn mưa ($\phi_T$): Xác lập đường cong $\phi_T \sim T$ cập nhật cho 12 trạm nghiên cứu; cung cấp công cụ chuyển đổi chính xác từ lượng mưa ngày $H_{n,p}$ sang lượng mưa thời đoạn ngắn $H_{T,p}$ phục vụ tính toán dòng chảy bằng mô hình thủy văn NAM - MIKE và làm tiêu chuẩn phân vùng mưa chi tiết.
- Phát triển và hoàn thiện 7 công thức thực nghiệm xác định $a_{T,p}$: Tổng kết, cải tiến và xây dựng mới 7 dạng công thức tính cường độ mưa (trong đó có 3 dạng công thức cải tiến mới), đồng thời xác định toàn bộ hệ số thực nghiệm ($m, S_p, A, B, \beta, a_{T0,p}$) cho 12 trạm nghiên cứu, giải quyết bài toán thiếu hụt thiết bị đo mưa tự ghi tại Việt Nam.
Kiến nghị và giải pháp đề xuất
- Cần hiệu chỉnh, cập nhật cơ sở dữ liệu khí tượng thủy văn trong các tiêu chuẩn thiết kế hiện hành (thay thế dần dữ liệu đo mưa trước năm 1987).
- Khi thiết kế công trình thoát nước nhỏ ở tần suất lũ cao ($p$ nhỏ như $1%$), cần cân nhắc giải pháp kết hợp giữa tăng khẩu độ công trình và tăng số lượng công trình thoát nước trên tuyến để hạn chế việc phải đắp nền đường quá cao, giảm chi phí xây dựng và nguy cơ ngập úng thượng lưu.
- Áp dụng các công thức thực nghiệm gián tiếp có kiểm chứng sai số để xác định $a_{\tau,p}$ cho các vùng chưa có trạm đo mưa tự ghi.
Hạn chế và hướng nghiên cứu tiếp
- Hạn chế về phạm vi số liệu: Do hạn chế về số lượng trạm có thiết bị đo mưa tự ghi đủ dài và liên tục, luận án mới chỉ thực hiện tính toán chi tiết và kiểm chứng cho 12 trạm khí tượng đại diện.
- Hướng nghiên cứu tiếp theo:
- Mở rộng phạm vi thu thập và xử lý số liệu đo mưa tự ghi cho toàn bộ hệ thống các trạm khí tượng trên phạm vi cả nước.
- Xây dựng bản đồ số hóa phân vùng mưa lũ tỷ lệ lớn cho toàn lãnh thổ Việt Nam dựa trên tiêu chí hệ số đặc trưng hình dạng cơn mưa $\phi_T$ và các hệ số vùng khí hậu.
- Tiếp tục theo dõi, cập nhật chuỗi số liệu quan trắc sau năm 2010 nhằm đánh giá sâu hơn các tác động dài hạn của biến đổi khí hậu cực đoan.
Giá trị tham khảo
Luận án là tài liệu chuyên khảo có giá trị ứng dụng và tham khảo cho nhiều đối tượng:
- Các kỹ sư, chuyên gia tư vấn thiết kế giao thông và thủy lợi: Sử dụng hệ thống 7 công thức thực nghiệm và các bảng tra tham số mưa ($H_{n,p}, \phi_T, a_{T,p}, S_p, m, A, B, \beta$) tại phụ lục của luận án để tính toán lưu lượng đỉnh lũ thiết kế $Q_p$ cho các dự án xây dựng đường bộ, đường sắt, đô thị.
- Các cơ quan quản lý nhà nước và viện nghiên cứu: Cung cấp cơ sở khoa học và số liệu thực nghiệm phục vụ công tác soát xét, sửa đổi, bổ sung các tiêu chuẩn kỹ thuật quốc gia như TCVN 9845:2013 và TCVN 7957:2008.
- Nghiên cứu sinh, học viên cao học và giảng viên: Phục vụ công tác giảng dạy, nghiên cứu chuyên sâu trong các bộ môn Đường bộ, Công trình Thủy văn, Kỹ thuật Hạ tầng Giao thông tại các trường đại học khối kỹ thuật.
Câu hỏi thường gặp
1. Tại sao đối với lưu vực nhỏ của công trình thoát nước trên đường, việc tính lưu lượng theo cường độ mưa tính toán ($a_{\tau,p}$) lại chính xác hơn lượng mưa ngày ($H_{n,p}$)?
Lưu vực của công trình thoát nước nhỏ trên đường thường có diện tích bé (từ dưới $1\text{ km}^2$ đến vài $\text{km}^2$), dẫn đến thời gian tập trung nước $\tau$ rất ngắn (từ vài phút đến vài giờ). Cường độ mưa tính toán $a_{\tau,p}$ phản ánh chính xác lượng mưa lớn nhất rơi trong chính khoảng thời gian tập trung nước $\tau$ bất lợi nhất và thể hiện được hình dạng cơn mưa. Ngược lại, lượng mưa ngày $H_{n,p}$ phân bố trong 24 giờ (1440 phút) không phản ánh được cường độ tập trung mưa cực hạn trong các thời đoạn ngắn.
2. Luận án đã khảo sát chuỗi số liệu đo mưa tại những trạm khí tượng cụ thể nào?
Luận án thu thập và phân tích chuỗi số liệu quan trắc liên tục 51 năm (1960 - 2010) tại 12 trạm khí tượng đại diện: TX. Điện Biên, TP. Tuyên Quang, TP. Lạng Sơn, Trạm Láng (Hà Nội), Trạm Hà Đông (Hà Nội), TX. Sơn Tây (Hà Nội), TP. Nghệ An, TP. Quảng Bình, TP. Khánh Hòa, TP. Buôn Ma Thuột (Đắk Lắk), TP. Cần Thơ.
3. Hệ số đặc trưng hình dạng cơn mưa ($\phi_T$) có vai trò gì trong tính toán thủy văn công trình?
Hệ số $\phi_T$ biểu thị tỷ lệ giữa lượng mưa thời đoạn $T$ so với lượng mưa ngày tính toán. Tham số này có vai trò:
- Xác định mô đuyn đỉnh lũ $A_p$ trong công thức tiêu chuẩn TCVN 9845:2013.
- Chuyển đổi lượng mưa ngày tính toán $H_{n,p}$ thành lượng mưa các thời đoạn ngắn $H_{T,p}$ phục vụ tính toán dòng chảy lũ trong mô hình Sôkôlôpsky và mô hình mô phỏng mưa - dòng chảy NAM - MIKE.
- Làm tiêu chí khoa học để phân vùng mưa rào phục vụ quy hoạch thoát nước.
4. Luận án đã đề xuất bao nhiêu công thức thực nghiệm xác định cường độ mưa tính toán ($a_{T,p}$)?
Luận án đã tổng kết, cải tiến và xây dựng 7 dạng công thức thực nghiệm xác định $a_{T,p}$ ứng với các điều kiện dữ liệu khác nhau của vùng thiết kế (tính theo $H_{n,p}$ và $\phi_T$; theo sức mưa $S_p$ và $m$; theo hệ số khí hậu $A, B$ và $m$; theo hệ số hồi quy $\beta$; theo cường độ mưa chuẩn $a_{T0,p}$; và 2 phương pháp trạm tựa nội suy theo $H_{n,p}$ hoặc $S_p$).
5. Tiêu chuẩn thống kê nào được sử dụng để kiểm định đường tần suất lý luận lượng mưa ngày?
Luận án sử dụng tiêu chuẩn kiểm định thống kê phi tham số Smirnov - Kolmogorov để đánh giá mức độ phù hợp giữa đường tần suất lý luận $H_{n,p}$ với chuỗi số liệu thực tế tại 12 trạm khí tượng ở mức ý nghĩa cho phép $\alpha = 5%$.
Kết luận
Luận án tiến sĩ của tác giả Nguyễn Anh Tuấn đã giải quyết một vấn đề kỹ thuật quan trọng trong tính toán thủy văn công trình giao thông tại Việt Nam: chuẩn hóa và lượng hóa các tham số mưa phục vụ xác định lưu lượng thiết kế công trình thoát nước nhỏ. Dựa trên chuỗi số liệu thực đo 51 năm (1960 - 2010) tại 12 trạm khí tượng điển hình, công trình đã cập nhật các giá trị lượng mưa ngày tính toán $H_{n,p}$, hệ số đặc trưng hình dạng cơn mưa $\phi_T$ và hoàn thiện hệ thống 7 công thức thực nghiệm tính cường độ mưa $a_{T,p}$. Kết quả nghiên cứu cung cấp cơ sở khoa học tin cậy để nâng cao độ chính xác trong thiết kế công trình hạ tầng giao thông và đóng góp tài liệu phục vụ công tác rà soát, hoàn thiện các tiêu chuẩn kỹ thuật quốc gia.
