Luận án Nghiên cứu cải tiến phương pháp xác định mô hình mua và lưu lượng tiêu thiết kế cho các hệ thống tiêu vùng đồng bằng Bắc Bộ

BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT 
TRƯỜNG ĐẠI HỌC THỦY LỢI 
NGUYỄN TH VIỆT HỒNG 
NGHIÊN CỨU CẢI TIẾN PHƯ NG PHÁP ÁC Đ NH Ô H NH Ư 
VÀ LƯU LƯỢNG TIÊU THIẾT Ế CHO CÁC HỆ THỐNG TIÊU 
V NG ĐỒNG BẰNG BẮC BỘ 
LUẬN ÁN TIẾN SĨ Ỹ THUẬT 
HÀ NỘI, NĂM 2018 
BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT 
TRƯỜNG ĐẠI HỌC THỦY LỢI 
NGUYỄN TH VIỆT HỒNG 
NGHIÊN CỨU CẢI TIẾN PHƯ NG PHÁP ÁC Đ NH Ô H NH Ư 
VÀ LƯU LƯỢNG TIÊU THIẾT Ế CHO CÁC HỆ THỐNG TIÊU 
V NG ĐỒNG BẰNG BẮC BỘ 
Chuyên ngành: ỹ thuật tài nguyên nước 
 ã số: 62-58-02-12 
NGƯỜI HƯỚNG DẪN HO HỌC: 1. PGS.TS Nguyễn Tuấn nh 
2. GS.TS Lê Chí Nguyện 
HÀ NỘI, NĂM 2018
i 
LỜI C ĐO N 
Tác giả xin cam đoan đây là công trình nghiên cứu của bản thân tác giả Các kết quả 
nghiên cứu và các kết lu n trong lu n án là trung th c, không sao ch p t b t k một 
nguồn nào và dưới b t k hình thức nào Việc tham khảo các nguồn tài liệu (nếu có) đã 
được th c hiện trích dẫn và ghi nguồn tài liệu tham khảo đúng quy định 
 Tác giả lu n án 
Chữ ký 
Nguyễn Th Việt Hồng 
ii 
LỜI CẢ N 
Trước hết tác giả xin được trân trọng cảm ơn sâu sắc tới PGS TS Nguyễn Tu n Anh, 
GS TS Lê Chí Nguyện đã định hướng về mặt khoa học và t n tình hướng dẫn để tác 
giả có thể hoàn thành lu n án này 
Tác giả xin được trân trọng cảm ơn đến các Thầy trong các Hội đồng ch m chuyên đề, 
Hội đồng ch m lu n án của tác giả, đã đóng góp các ý kiến quý báu để tác giả có thể 
hoàn thành lu n án Tác giả xin trân trọng cảm ơn Ban giám hiệu, Phòng Đào tạo Đại 
học và Sau đại học, Khoa Kỹ thu t tài nguyên nước; Bộ môn Kỹ thu t hạ tầng và Phát 
triển nông thôn Trường Đại học Thủy lợi đã tạo mọi điều kiện thu n lợi để tác giả 
hoàn thành lu n án 
Tác giả xin chân thành cảm ơn Trường Cao đẳng Thủy lợi Bắc bộ, nơi tác giả đang 
công tác, đã tạo điều kiện, giúp đỡ, động viên tác giả trong suốt thời gian th c hiện 
lu n án Xin được cảm ơn Sở Nông nghiệp và Phát triển nông thôn Hà Nam, Sở Xây 
d ng Hà Nam, Công ty cổ phần môi trường và công trình đô thị Hà Nam đã trao đổi, 
cung c p tài liệu cũng như kiến thức th c tế quý báu giúp tác giả hoàn thành các nội 
dung nghiên cứu của lu n án 
Cuối cùng tác giả xin bày tỏ lòng biết ơn đến gia đình, bạn bè luôn quan tâm, động 
viên, giúp đỡ tác giả vượt qua mọi khó khăn để tác giả có thể hoàn thành được lu n án 
của mình 
 Tác giả lu n án 
Chữ ký 
Nguyễn Th Việt Hồng 
iii 
 ỤC LỤC 
 ỤC LỤC iii 
D NH ỤC BẢNG BIỂU ..........................................................................................ix 
D NH ỤC CÁC TỪ VIẾT TẮT .............................................................................xi 
 Ở Đ U...................... ................................................................................................... 1 
1 Tính c p thiết của đề tài ............................................................................................... 1 
2 M c tiêu nghiên cứu .................................................................................................... 2 
3 Đối tượng và phạm vi nghiên cứu ............................................................................... 3 
4 Nội dung nghiên cứu ................................................................................................... 3 
5 Phương pháp nghiên cứu ............................................................................................. 3 
6 Ý nghĩa khoa học và th c tiễn ..................................................................................... 4 
7 C u trúc lu n án ........................................................................................................... 4 
CHƯ NG 1 TỔNG QU N ....................................................................................... 5 
1.1 Tổng quan về mô hình mưa tiêu thiết kế .......................................................... 5 
1.1.1 Khái niệm về mô hình mưa tiêu thiết kế ................................................... 5 
1.1.2 Các nghiên cứu về xác định mô hình mưa tiêu thiết kế ............................ 6 
1.1.3 Các nghiên cứu xác định mô hình mưa tiêu thiết kế ở Việt Nam ........... 11 
1.2 Tổng quan các nghiên cứu về quan hệ lượng mưa – thời gian mưa – tần su t 
(DDF) và quan hệ cường độ mưa-thời gian mưa-tần su t (IDF) .............................. 16 
1.2.1 Tổng quan các nghiên cứu về quan hệ lượng mưa – thời gian mưa – tần 
su t...................................................................................................................... 16 
1.2.2 Tổng quan các nghiên cứu về quan hệ cường độ mưa – thời gian – tần 
su t.............. ....................................................................................................... 17 
1.3 Tổng quan các nghiên cứu về xây d ng bản đồ đẳng trị mưa ........................ 22 
1.4 Tổng quan các nghiên cứu về tính lưu lượng tiêu thiết kế và hệ số hiệu chỉnh 
lưu lượng theo bước thời gian mưa ........................................................................... 28 
1.4.1 Các nghiên cứu về tính lưu lượng tiêu thiết kế cho các khu v c dân cư, 
đô thị.... ............................................................................................................... 28 
1.4.2 Các nghiên cứu về tính lưu lượng tiêu thiết kế cho các vùng nông nghiệp 
và hỗn hợp dân cư, nông nghiệp ....................................................................... 29 
1.5 Giới thiệu về đồng bằng Bắc Bộ ..................................................................... 33 
iv 
1.5.1 Điều kiện t nhiên vùng đồng bằng Bắc Bộ ........................................... 33 
1.5.2 Phát triển kinh tế xã hội trong vùng ........................................................ 36 
1.5.3 Tình hình đo mưa trong khu v c ............................................................. 37 
1.5.4 Tình hình tiêu nước trong khu v c .......................................................... 39 
1.6 Kết lu n chương 1 ........................................................................................... 42 
CHƯƠNG 2 XÂY DỰNG PHƯƠNG PHÁP XÁC ĐỊNH MÔ HÌNH MƯA VÀ 
LƯU LƯỢNG THIẾT KẾ CHO HỆ THỐNG TIÊU .............................................. 45 
2.1 Phương pháp xây d ng quan hệ lượng mưa - thời gian - tần su t và (DDF) 
quan hệ cường độ mưa - thời gian - tần su t (IDF) thời đoạn ngắn, t đó thiết l p 
công thức tính cường độ mưa cho các trạm đo mưa t ghi ...................................... 45 
2.1.1 Phương pháp xây d ng quan hệ lượng mưa – thời gian – tần su t ......... 46 
2.1.2 Thiết l p công thức tính cường độ mưa (quan hệ cường độ mưa – thời 
gian – tần su t) cho các trạm đo mưa t ghi ...................................................... 51 
2.2 Phương pháp xây d ng các bản đồ đẳng trị lượng mưa thiết kế tương ứng với 
thời gian mưa và tần su t khác nhau cho vùng Đồng bằng Bắc bộ........................... 51 
2.2.1 L a chọn phần mềm vẽ bản đồ ............................................................... 52 
2.2.2 L a chọn phương pháp nội suy trong phần mềm .................................... 53 
2.3 Phương pháp mới đề xu t về l a chọn mô hình mưa tiêu thiết kế phù hợp 
nh t cho các khu v c trồng lúa vùng đồng bằng Bắc Bộ bằng mô phỏng toàn liệt các 
tr n mưa đã xảy ra trong th c tế ............................................................................... 55 
2.3.1 Giới thiệu phương pháp .......................................................................... 56 
2.3.2 Xây d ng các mô hình mưa điển hình .................................................... 58 
2.3.3 Phương pháp tính toán hệ số tiêu cho lúa ............................................... 59 
2.4 Phương pháp xây d ng hệ số hiệu chỉnh lưu lượng tiêu thiết kế t mô hình 
mưa giờ và mô hình mưa ngày cho vùng tiêu hỗn hợp nông nghiệp và dân cư đô thị 
 ........................................................................................................................64 
2.4.1 Tính toán lưu lượng tiêu của hệ thống theo mô hình SWMM ................ 65 
2.4.2 Tính toán lưu lượng tiêu của hệ thống theo TCVN 10406:2015 ............ 71 
2.5 Kết lu n chương 2 ........................................................................................... 71 
CHƯ NG 3 ẾT QUẢ NGHIÊN CỨU VÀ THẢO LUẬN ................................ 73 
3.1 Kết quả xây d ng quan hệ lượng mưa - thời gian - tần su t (DDF) và quan hệ 
cường độ mưa - thời gian - tần su t (IDF) ................................................................ 73 
v 
3.1.1 Kết quả xây d ng quan hệ DDF .............................................................. 73 
3.1.2 Đánh giá khả năng ứng d ng của công thức tính cường độ mưa đã thiết 
l p.............. ......................................................................................................... 84 
3.2 Kết quả xây d ng các bản đồ đẳng trị mưa thời đoạn ngắn ứng với tần su t 
thiết kế P = 10% cho vùng đồng bằng Bắc bộ........................................................... 97 
3.2.1 Kết quả xây d ng các bản đồ đẳng trị lượng mưa ứng với tần su t P = 
10%.............. ....................................................................................................... 97 
3.2.2 Kết quả xây d ng các bản đồ đẳng trị tham số của phương trình cường 
độ mưa ứng với tần su t P = 10% ................................................................... 104 
3.3 Kết quả xác định mô hình mưa thiết kế hợp lý cho vùng trồng lúa ứng d ng 
cho khu v c Hà Nam ............................................................................................... 107 
3.3.1 Tính toán và vẽ đường tần su t ............................................................. 107 
3.3.2 Kết quả l a chọn các tr n mưa điển hình .............................................. 108 
3.3.3 Kết quả l a chọn tr n mưa thiết kế ....................................................... 114 
3.4 Kết quả xây d ng hệ số hiệu chỉnh lưu lượng cho vùng tiêu hỗn hợp ......... 122 
3.4.1 Thiết l p các thông số giả định của hệ thống ........................................ 123 
3.4.2 Xây d ng hệ số hiệu chỉnh lưu lượng trong trường hợp chưa x t đến ảnh 
hưởng của hồ điều hòa trong hệ thống ............................................................ 124 
3.4.3 Xây d ng hệ số hiệu chỉnh lưu lượng có x t đến ảnh hưởng của hồ điều 
hòa...... .............................................................................................................. 128 
3.5 Kết lu n chương 3 ......................................................................................... 134 
 ẾT LUẬN VÀ IẾN NGH ................................................................................... 135 
D NH ỤC CÁC BÀI BÁO ĐÃ CÔNG BỐ.........................................................138 
TÀI LIỆU TH HẢO.........................................................................................139 
PHỤ LỤC........ ........................................................................................................... 147 
vi 
D NH ỤC CÁC H NH ẢNH 
Hình 1 1 Phân bố xác su t của các tr n mưa nhóm thứ nh t tại Illinois[1] .................... 7 
Hình 1 2 Biểu đồ xác su t 10% các tr n mưa nhóm thứ nh t tại Illinois ........................ 8 
Hình 1 3 Biểu đồ xác su t 50% các tr n mưa nhóm thứ nh t tại Illinois ........................ 8 
Hình 1 4 Biểu đồ xác su t 90% các tr n mưa nhóm thứ nh t tại Illinois ........................ 8 
Hình 1 5 Các đường quá trình mưa lũy tích mưa trong 24 giờ của SCS Hoa K ........... 9 
Hình 1 6 Phạm vi áp d ng các loại đường quá trình mưa của SCS trên lãnh thổ Hoa K 
 ......................................................................................................................................... 9 
Hình 1 7 Biểu đồ mưa thiết kế xây d ng bằng phương pháp khối xen kẽ của Chow [3]
 ....................................................................................................................................... 10 
Hình 1 8 Biểu đồ quá trình mưa thiết kế hình tam giác ................................................ 10 
Hình 1.9 Các đường cong DDF của mưa lớn nh t thời đoạn ngắn tại trạm Phủ lý .... 187 
Hình 1.10 Các đường cong IDF của mưa lớn nh t tại Chicago, Hoa K [3] ................ 18 
Hình 1.11 Phân bố mưa (in) thời k xu t hiện lại 2 năm, thời gian mưa 60 phút trên 
lãnh thổ Hoa K [38] ..................................................................................................... 24 
Hình 1.12 Bản đồ đẳng trị tham số của phương trình lượng mưa lưu v c sông Tre 
Venezie (Ý) .................................................................................................................... 24 
Hình 1.13 Bản đồ xác định khu v c áp d ng các công thức tính lượng mưa được thiết 
l p thuộc lưu v c sông Tre Venezie (Ý) ....................................................................... 25 
Hình 1.14 Một ví d về đồ thị quan hệ độ sâu mưa và diện tích mưa để tính các giá trị 
trung bình của mưa diện (Tổ chức khí tượng thế giới, 1983) ....................................... 27 
Hình 1.15 Bản đồ vị trí vùng Đồng bằng Bắc bộ .......................................................... 33 
Hình 1.16 Bản đồ mạng lưới trạm khí tượng thủy văn vùng ĐBBB ............................ 38 
Hình 1.17 Bản đồ hệ thống thủy nông đồng bằng sông Hồng [51] ............................... 40 
Hình 1 18 Sơ đồ nội dung nghiên cứu và cách tiếp c n...............................................44 
Hình 2 1 Sơ đồ tiến trình xây d ng quan hệ DDF và IDF ............................................ 45 
Hình 2.2 Quan hệ giữa lượng mưa và thời gian mưa ................................................... 49 
Hình 2 3 Đường quan hệ Hd được mô tả bởi phương trình (2-14) và (2-15) và những 
điểm th c nghiệm .......................................................................................................... 50 
Hình 2 4 Bản đồ các trạm đo mưa t ghi vùng ĐBBB .................................................. 52 
Hình 2 5 Phương pháp nội suy IDW ............................................................................. 54 
vii 
Hình 2 6 Sơ đồ khối tính toán l a chọn MHM tiêu cho lúa phù hợp nh t bằng phương 
pháp so sánh mô phỏng toàn liệt ................................................................................... 57 
Hình 2 7 Dòng chảy t ô ruộng vào kênh .................................................................... 59 
Hình 2 8 Sơ đồ khối tính hệ số tiêu thiết kế .................................................................. 62 
Hình 2 9 Sơ đồ khối mô phỏng hệ số tiêu q~t toàn liệt ................................................. 63 
Hình 2 10 Sơ đồ một hệ thống tiêu hỗn hợp dân cư và nông nghiệp (TB: trạm bơm) . 64 
Hình 2 11 Sơ đồ khối tính toán dòng chảy cho lưu v c nghiên cứu ............................. 70 
Hình 3 1 Quan hệ giữa lượng mưa và thời gian mưa trạm Phủ Lý .............................. 74 
Hình 3 2 Đường cong H  d ứng với  ... * (x1 + H1) * (x1 + H1) / 8) 
 F2 = W - x2 * S * 10000 - beta * 3600 * delt * mt * sig * Sqr(2 * 9.81) * b * Sqr((x2 
+ H1) * (x2 + H1) * (x2 + H1) / 8) 
Do Until (F1 * F2) < 0 
 x2 = x2 + 5 
 F1 = W - x1 * S * 10000 - beta * 3600 * delt * mt * sig * Sqr(2 * 9.81) * b * Sqr((x1 
+ H1) * (x1 + H1) * (x1 + H1) / 8) 
 F2 = W - x2 * S * 10000 - beta * 3600 * delt * mt * sig * Sqr(2 * 9.81) * b * Sqr((x2 
+ H1) * (x2 + H1) * (x2 + H1) / 8) 
 Loop 
 Do 
 d = (x1 + x2) / 2 
 F1 = W - x1 * S * 10000 - beta * 3600 * delt * mt * sig * Sqr(2 * 9.81) * b * Sqr((x1 
+ H1) * (x1 + H1) * (x1 + H1) / 8) 
176 
 Fd = W - d * S * 10000 - beta * 3600 * delt * mt * sig * Sqr(2 * 9.81) * b * Sqr((d + 
H1) * (d + H1) * (d + H1) / 8) 
 F2 = W - x2 * S * 10000 - beta * 3600 * delt * mt * sig * Sqr(2 * 9.81) * b * Sqr((x2 
+ H1) * (x2 + H1) * (x2 + H1) / 8) 
 If Fd * F1 < 0 Then 
 x1 = x1 
 x2 = d 
 Else 
 If Fd * F1 > 0 Then 
 x1 = d 
 x2 = x2 
 Else 
 If Fd * F1 = 0 Then 
 Fd = 0 
 End If 
 End If 
 End If 
 Loop While Abs(x2 - x1) > ep 
 H2 = d 
 Hm = (H2 + H1) / 2 '(*------H mean-----*) 
 Q = (W - H2 * S * 10000) / (beta * delt * 3600) '(* m3/s *) 
 If Q < 0 Then 
 Q = 0 
 End If 
End If 
177 
 a = 1000 * (hs + Hm) '(* Chieu sau lop nuoc: mm *) 
 If a >= a1 Then 
 tn1 = tn1 + 1 
 End If 
 If a >= a2 Then 
 tn2 = tn2 + 1 
 End If 
 If a >= a3 Then 
 tn3 = tn3 + 1 
 End If 
 If a >= a4 Then 
 tn4 = tn4 + 1 
 End If 
 If a >= a5 Then 
 tn5 = tn5 + 1 
 End If 
 H1 = H2 
 qr(i) = 1000 * Q 
 If qr_max < qr(i) Then 
 qr_max = qr(i) 
 End If 
 Next 'Quay lai vong lap for 
Loop Until (tn1 <= 1) And (tn2 <= 2) And (tn3 <= 3) And (tn4 <= 4) And (tn5 <= 5) 
And (a <= 110) 
Print 
178 
Print " KET QUA TINH TOAN " 
Print " Thoi doan Mua (mm) He so tieu (l/s/ha)" 
Write #2, " KET QUA TINH TOAN" 
Write #2, " Thoi doan He so tieu(l/s/ha)" 
For i = 1 To N 
Print "", i, Pr(i), Format(qr(i), "##.##") 
Write #2, " ", i, " ", Format(qr(i), "##.##") 
Next 
Write #2, "Chieu rong tran mat ruong, b = ", Format(b, "##.###"), "m/ha" 
Write #2, "He so tieu cho lua max, qmax = ", Format(qr_max, "##.##"), "l/s/ha" 
Close 2 
Print "He so tieu cho lua max, qmax = ", Format(qr_max, "##.##"), "l/s/ha" 
Print 
Print "Chieu rong tran mat ruong, b = ", Format(b, "##.###"), "m/ha" 
Print 
Print "Ghi chu: Ket qua da duoc luu tai file ket qua." 
End Sub 
Private Sub cmdXoa_Click() 
Cls 
End Sub 
Private Sub cmdKet_Click() 
Beep 
End 
End Sub 
PL 3.4: Chương trình tính hệ số tiêu vận hành 
Private Sub cmdTinh_HSTVH_Click() 
179 
Dim Tenfile1$, Tenfile2$, Tenfile3$, st$, Tenthumuc$, Tennam$, tentran$ 
Dim i%, k%, n1%, N% 
Dim x!, delt!, Tmua!, Tdongchay!, Tsaumua!, qr_max! 
Dim qr1!, Hp!, qmax!, qmax1!, Qtongmax! 
Dim Hc1!, Hc2! 
Dim P!(2000), Pr!(2000), Zsong!(2000), qr!(2000), Qtong!(2000), Qlua!(2000), 
Qmau!(2000) 
Dim Qthocu!(2000), Qctcc!(2000), Qao_kodt!(2000), Qao_codt!(2000) 
Dim b!, Pt!, inf!, Ho!, H!, sig!, a!, S!, H1!, H2!, Q!, Zbhut!, Zbthao!, Zbthao1!, 
Zbthao2!, Hdh! 
Dim W!, x1!, x2!, Fd!, F1!, F2!, d!, Hm! 
Const ep = 0.00001 
Const beta = 1.03 
Const hs = 0.1 
Const mt = 0.34 
Const hsm = 0.9 
'----------Nhap du lieu ve tram bom, luu vuc tieu 
Tenthumuc = "D:\HSTL" 
Tenfile1 = Tenthumuc + "\" + "b_tran.txt" 
Open Tenfile1 For Input As 1 
 Input #1, st$ 
 Input #1, b 
 Close 1 
 '------Nhap so lieu tran mua, muc nuoc song 
Tennam = InputBox("Nhap ten nam tinh toan (VD: 1990):") 
180 
tentran = InputBox("Nhap ten tran mua tinh toan (VD: 1):") 
Tenfile2 = Tenthumuc + "\" + Tennam + "\" + tentran + ".txt" 
Open Tenfile2 For Input As 2 
 Input #2, st$ 
 Input #2, st$ 
 Input #2, delt 'gio 
 Input #2, st$ 
 Input #2, Tsaumua 
 Input #2, st$ 
 Input #2, n1 'so thoi doan mua 
 Input #2, st$ 
For i = 1 To n1 
 Input #2, P(i) 'mm/thoi doan 
Next 
Close 2 
Tenfile3 = Tenthumuc + "\" + Tennam + "\TL" + tentran + ".txt" 
Open Tenfile3 For Append As 3 
 Write #3, "KET QUA TINH TOAN " 
 Print 
 Print 
 Print " KET QUA TINH TOAN " 
 Print " Thoi doan Mua (mm) He so tieu (l/s/ha) " 
 Write #3, "Thoi doan - He so tieu (l/s/ha) " 
 '(*--------BEGIN COMPUTING------*) 
181 
 Tmua = n1 * delt 'Thoi gian mua, gio 
 Tdongchay = Tmua + Tsaumua 'Thoi gian mo phong, gio 
 N = Round(Tdongchay / delt) 
 For i = 1 To n1 
 Pr(i) = P(i) 
 Next 
 For i = (n1 + 1) To N 
 Pr(i) = 0 
 Next 
'(*---------Tinh he so tieu vung lua------*) 
 sig = 1 
 H1 = 0 
 qr_max = 0 
 S = 1 'ha 
 inf = 0.23 'mm/h = 5.5 mm/ngay 
 H = inf 
'(*----P: mm/delt, delt: gio, S: ha, q:m3/s, a:mm, H1: m, H2:m---*) 
 For i = 1 To N 
 Pt = Pr(i) 
W = (Pt / 1000) * beta * S * 10000 - (H / 1000) * delt * beta * S * 10000 + H1 * S * 
10000 
 If W <= 0 Then 
 Hm = 0 
 Q = 0 
 H2 = 0 
182 
 Else 
 '(*--------Tim H2-----------*) 
 x1 = 0 
 x2 = 0 
 F1 = W - x1 * S * 10000 - beta * 3600 * delt * mt * sig * Sqr(2 * 9.81) * b * Sqr((x1 + 
H1) * (x1 + H1) * (x1 + H1) / 8) 
 F2 = W - x2 * S * 10000 - beta * 3600 * delt * mt * sig * Sqr(2 * 9.81) * b * Sqr((x2 + 
H1) * (x2 + H1) * (x2 + H1) / 8) 
Do Until (F1 * F2) < 0 
 x2 = x2 + 5 
 F1 = W - x1 * S * 10000 - beta * 3600 * delt * mt * sig * Sqr(2 * 9.81) * b * Sqr((x1 + 
H1) * (x1 + H1) * (x1 + H1) / 8) 
 F2 = W - x2 * S * 10000 - beta * 3600 * delt * mt * sig * Sqr(2 * 9.81) * b * Sqr((x2 + 
H1) * (x2 + H1) * (x2 + H1) / 8) 
 Loop 
 Do 
 d = (x1 + x2) / 2 
 F1 = W - x1 * S * 10000 - beta * 3600 * delt * mt * sig * Sqr(2 * 9.81) * b * Sqr((x1 + 
H1) * (x1 + H1) * (x1 + H1) / 8) 
 Fd = W - d * S * 10000 - beta * 3600 * delt * mt * sig * Sqr(2 * 9.81) * b * Sqr((d + 
H1) * (d + H1) * (d + H1) / 8) 
 F2 = W - x2 * S * 10000 - beta * 3600 * delt * mt * sig * Sqr(2 * 9.81) * b * Sqr((x2 + 
H1) * (x2 + H1) * (x2 + H1) / 8) 
 If Fd * F1 < 0 Then 
 x1 = x1 
 x2 = d 
 Else 
 If Fd * F1 > 0 Then 
 x1 = d 
 x2 = x2 
 Else 
 If Fd * F1 = 0 Then 
 Fd = 0 
 End If 
 End If 
 End If 
 Loop While Abs(x2 - x1) > ep 
 H2 = d 
183 
 Hm = (H2 + H1) / 2 '(*------H mean (m)-----*) 
 Q = (W - H2 * S * 10000) / (beta * delt * 3600) '(* m3/s *) 
 If Q < 0 Then 
 Q = 0 
 End If 
End If 
 a = 1000 * (hs + Hm) '(* water layer depth: mm *) 
 '--------------------------------------- 
 H1 = H2 
 qr(i) = 1000 * Q 
 If qr_max < qr(i) Then 
 qr_max = qr(i) 
 End If 
 Print "", i, Pr(i), Format(qr(i), "##.##") 
 Write #3, " ", i, " ", Format(qr(i), "##.##") 
 Next 
Write #3, "He so tieu cho lua max, qmax = ", Format(qr_max, "##.##"), "l/s/ha" 
Close 3 
Print 
Print "He so tieu cho lua max, qmax =", Format(qr_max, "##.##"), "l/s/ha" 
Print 
Print "Ghi chu: Ket qua da duoc luu tai file ket qua." 
End Sub 
184 
 Phụ lục 3.5: Kết quả tính hệ số tiêu toàn liệt lớn nhất ứng với btran
max 
STT Năm 
Hệ số tiêu vận hành lớn nhất củ năm qmaxmax (l/s/ha) 
btràn
1
 = 0,33 
(m/ha) 
btràn
2
 = 0,29 
(m/ha) 
btràn
3
 = 0,35 
(m/ha) 
btràn
4
 = 0,42 
(m/ha) 
btràn
5 = 0,43 
(m/ha) 
1 1985 15.97 15.87 16.47 16.84 16.95 
2 1986 15.09 15.05 15.16 15.68 15.99 
3 1987 3.87 3.82 4.14 4.35 4.4 
4 1988 4.15 4.11 4.36 4.5 4.66 
5 1989 9.95 9.89 10.28 10.49 10.54 
6 1990 12.67 12.57 13.1 13.36 13.42 
7 1991 4.47 4.40 4.86 5.17 5.26 
8 1992 11.43 11.33 11.96 12.36 12.48 
9 1993 16.03 15.88 16.79 17.33 17.48 
10 1994 21.77 21.57 22.75 23.45 23.64 
11 1995 8.89 8.83 9.16 9.31 9.34 
12 1996 21.02 20.73 22.51 23.72 24.07 
13 1997 10.76 10.59 11.68 12.43 12.66 
14 1998 4.96 4.91 5.21 5.38 5.43 
15 1999 5.43 5.38 5.68 5.83 5.87 
16 2000 8.77 8.73 8.97 9.06 9.08 
17 2001 8.37 8.31 8.7 9.0 9.1 
18 2002 10.64 10.52 11.25 11.7 11.82 
19 2003 10.82 10.73 11.29 11.63 11.73 
20 2004 5.95 5.89 6.28 6.51 6.57 
21 2005 8.91 8.83 9.32 9.6 9.67 
22 2006 7.34 7.28 7.62 7.79 7.84 
23 2007 8.63 8.55 9.03 9.30 9.37 
24 2008 16.28 16.21 16.6 16.75 16.77 
25 2009 7.61 7.52 8.05 8.37 8.46 
26 2010 7.7 7.66 7.91 8.04 8.08 
27 2011 6.08 6.03 6.34 6.5 6.54 
28 2012 7.08 7.0 7.51 7.84 7.93 
29 2013 5.84 5.79 6.14 6.35 6.4 
30 2014 11.34 11.27 11.62 11.75 11.77 
185 
PL 3.6: Kết quả tính hệ số tiêu toàn liệt lớn nhất ứng với btran
max 
STT Năm 
Hệ số tiêu vận hành lớn nhất củ năm qmaxmax (l/s/ha) 
btràn
6
 = 0,4 
(m/ha) 
btràn
7
 = 0, 32 
(m/ha) 
btràn
8
 = 0,42 
(m/ha) 
1 1985 16.79 16.24 16.9 
2 1986 15.53 15.15 15.84 
3 1987 4.32 4.02 4.38 
4 1988 4.48 4.26 4.51 
5 1989 10.46 10.13 10.52 
6 1990 13.32 12.91 13.39 
7 1991 5.12 4.67 5.21 
8 1992 12.3 11.72 12.42 
9 1993 17.25 16.44 17.41 
10 1994 23.35 22.3 23.55 
11 1995 9.28 9.04 9.33 
12 1996 23.53 21.8 23.9 
13 1997 12.32 11.24 12.55 
14 1998 5.36 5.1 5.41 
15 1999 5.81 5.57 5.86 
16 2000 9.03 8.89 9.07 
17 2001 8.95 8.54 9.05 
18 2002 11.63 10.97 11.76 
19 2003 11.58 11.07 11.68 
20 2004 6.48 6.15 6.54 
21 2005 9.56 9.13 9.63 
22 2006 7.77 7.5 7.81 
23 2007 9.26 8.85 9.34 
24 2008 16.73 16.47 16.76 
25 2009 8.32 7.84 8.41 
26 2010 8.02 7.84 8.06 
27 2011 6.48 6.23 6.52 
28 2012 7.79 7.31 7.89 
29 2013 6.32 6.0 6.38 
30 2014 11.73 11.5 11.76 
186 
PL 3.7: Kết quả vẽ ường t n suất hệ số tiêu toàn liệt 
Hình 1: Đường tần su t hệ số tiêu toàn liệt ứng với btràn
1
 = 0,33m/ha 
Với P = 10% tra được qTL
6
 = 15,57 (l/s.ha) 
Hình 2: Đường tần su t hệ số tiêu toàn liệt ứng với btràn
2
 = 0,29 m/ha 
Với P = 10% tra được qTL
6
 = 14,97 (l/s.ha) 
187 
Hình 3: Đường tần su t hệ số tiêu toàn liệt ứng với btràn
3
 = 0,35m/ha 
Với P = 10% tra được qTL
6
 = 16,44 (l/s.ha) 
Hình 4: Đường tần su t hệ số tiêu toàn liệt ứng với btràn
4
 = 0,38 m/ha 
Với P = 10% tra được qTL
6
 = 17,03 (l/s.ha) 
188 
Hình 5: Đường tần su t hệ số tiêu toàn liệt ứng với btràn
5
 = 0,43 m/ha 
Với P = 10% tra được qTL
6
 = 17,11 (l/s.ha) 
Hình 6: Đường tần su t hệ số tiêu toàn liệt ứng với btràn
6
 = 0,41 m/ha 
Với P = 10% tra được qTL
6
 = 16,87 (l/s.ha) 
189 
Hình 7: Đường tần su t hệ số tiêu toàn liệt ứng với btràn
7
 = 0,32 m/ha 
Với P = 10% tra được qTL
7
 = 16,24 (l/s.ha) 
Hình 8: Đường tần su t hệ số tiêu toàn liệt ứng với btràn
8
 = 0,42 m/ha 
Với P = 10% tra được qTL
8
 = 17,03 (l/s.ha) 
190 
PHỤ LỤC 4 
 ết quả tính toán hệ số tiêu sơ bộ củ hệ thống theo mô hình mư ngày 
 (Theo TCVN 10604:2015) 
PL 4.1: Kết quả tính toán hệ số tiêu sơ bộ của hệ thống theo cơ c u sử d ng đ t: 
 l/ dc = 0,8/0,2 
Ngày 
Pi 
(mm) 
 l ql l.qli dc 
qdc 
(l/s/ha) 
 l.qdci 
qi 
(l/s/ha) 
 
(ha) 
Q 
(m
3
/s) 
 
(ha) 
Q 
(m3/s) 
 
(ha) 
Q 
(m
3
/s) 
1 11.33 0.8 0.07 0.06 0.20 0.85 0.17 0.23 300 0.07 500 0.11 1000 0.23 
2 40.68 0.8 1.26 1.01 0.20 3.06 0.61 1.62 300 0.49 500 0.81 1000 1.62 
3 220.71 0.8 12.84 10.27 0.20 16.60 3.32 13.59 300 4.08 500 6.80 1000 13.59 
4 78.63 0.8 16.81 13.45 0.20 5.91 1.18 14.63 300 4.39 500 7.32 1000 14.63 
5 44.63 0.8 6.63 5.30 0.20 3.36 0.67 5.98 300 1.79 500 2.99 1000 5.98 
6 
0.8 2.71 2.17 0.20 0.00 0.00 2.17 300 0.65 500 1.08 1000 2.17 
7 
0.8 0.51 0.41 0.20 0.00 0.00 0.41 300 0.12 500 0.20 1000 0.41 
PL 4.2: Kết quả tính toán hệ số tiêu sơ bộ của hệ thống theo cơ c u sử d ng đ t: 
 l/ dc = 0,7/0,3 
Ngày 
Pi 
(mm) 
 l ql l.qli dc 
qdc 
(l/s/ha) 
 l.qdci 
qi 
(l/s/ha) 
 
(ha) 
Q 
(m
3
/s) 
 
(ha) 
Q 
(m3/s) 
 
(ha) 
Q 
(m
3
/s) 
1 11.33 0.7 0.07 0.05 0.30 0.85 0.26 0.30 300 0.09 500 0.15 1000 0.30 
2 40.68 0.7 1.26 0.88 0.30 3.06 0.92 1.80 300 0.54 500 0.90 1000 1.80 
3 220.7 0.7 12.84 8.99 0.30 16.60 4.98 13.97 300 4.19 500 6.98 1000 13.97 
4 78.63 0.7 16.81 11.77 0.30 5.91 1.77 13.54 300 4.06 500 6.77 1000 13.54 
5 44.63 0.7 6.63 4.64 0.30 3.36 1.01 5.65 300 1.69 500 2.82 1000 5.65 
6 
0.7 2.71 1.90 0.30 0.00 0.00 1.90 300 0.57 500 0.95 1000 1.90 
7 
0.7 0.51 0.36 0.30 0.00 0.00 0.36 300 0.11 500 0.18 1000 0.36 
191 
PL 4.3: Kết quả tính toán hệ số tiêu sơ bộ của hệ thống theo cơ c u sử d ng đ t: 
 l/ dc = 0,6/0,4 
Ngày 
Pi 
(mm) 
 l ql l.qli dc 
qdc 
(l/s/ha) 
 l.qdci 
qi 
(l/s/ha) 
 
(ha) 
Q 
(m
3
/s) 
 
(ha) 
Q 
(m3/s) 
 
(ha) 
Q 
(m
3
/s) 
1 11.33 0.6 0.07 0.04 0.40 0.85 0.34 0.38 300 0.11 500 0.19 1000 0.38 
2 40.68 0.6 1.26 0.76 0.40 3.06 1.22 1.98 300 0.59 500 0.99 1000 1.98 
3 220.7 0.6 12.84 7.70 0.40 16.60 6.64 14.34 300 4.30 500 7.17 1000 14.34 
4 78.63 0.6 16.81 10.09 0.40 5.91 2.37 12.45 300 3.74 500 6.23 1000 12.45 
5 44.63 0.6 6.63 3.98 0.40 3.36 1.34 5.32 300 1.60 500 2.66 1000 5.32 
6 
0.6 2.71 1.63 0.40 0.00 0.00 1.63 300 0.49 500 0.81 1000 1.63 
7 
0.6 0.51 0.31 0.40 0.00 0.00 0.31 300 0.09 500 0.15 1000 0.31 
PL 4.4: Kết quả tính toán hệ số tiêu sơ bộ của hệ thống theo cơ c u sử d ng đ t: 
 l/ dc = 0,5/0,5 
Ngày 
Pi 
(mm) 
 l ql l.qli dc 
qdc 
(l/s/ha) 
 l.qdci 
qi 
(l/s/ha) 
 
(ha) 
Q 
(m
3
/s) 
 
(ha) 
Q 
(m3/s) 
 
(ha) 
Q 
(m
3
/s) 
1 11.33 0.5 0.07  0.50   0.46  0.14  0.23  0.46 
2 40.68 0.5 1.26 0.63 0.50 3.06 1.53 2.16 300 0.65 500 1.08 1000 2.16 
3 220.7 0.5 12.84 6.42 0.50 16.60 8.30 14.72 300 4.42 500 7.36 1000 14.72 
4 78.63 0.5 16.81 8.41 0.50 5.91 2.96 11.36 300 3.41 500 5.68 1000 11.36 
5 44.63 0.5 6.63 3.32 0.50 3.36 1.68 4.99 300 1.50 500 2.50 1000 4.99 
6 
0.5 2.71 1.36 0.50 0.00 0.00 1.36 300 0.41 500 0.68 1000 1.36 
7 
0.5 0.51 0.26 0.50 0.00 0.00 0.26 300 0.08 500 0.13 1000 0.26 
192 
PL 4.5: Kết quả tính toán hệ số tiêu sơ bộ của hệ thống theo cơ c u sử d ng đ t: 
 l/ dc = 0,6/0,4 
Ngày 
Pi 
(mm) 
 l ql l.qli dc 
qdc 
(l/s/ha) 
 l.qdci 
qi 
(l/s/ha) 
 
(ha) 
Q 
(m
3
/s) 
 
(ha) 
Q 
(m3/s) 
 
(ha) 
Q 
(m
3
/s) 
1 11.33 0.4 0.07 0.03 0.40 0.85 0.34 0.37 300 0.11 500 0.18 1000 0.37 
2 40.68 0.4 1.26 0.50 0.40 3.06 1.22 1.73 300 0.52 500 0.86 1000 1.73 
3 220.7 0.4 12.84 5.14 0.40 16.60 6.64 11.77 300 3.53 500 5.89 1000 11.77 
4 78.63 0.4 16.81 6.72 0.40 5.91 2.37 9.09 300 2.73 500 4.54 1000 9.09 
5 44.63 0.4 6.63 2.65 0.40 3.36 1.34 3.99 300 1.20 500 2.00 1000 3.99 
6 
0.4 2.71 1.08 0.40 0.00 0.00 1.08 300 0.33 500 0.54 1000 1.08 
7 
0.4 0.51 0.20 0.40 0.00 0.00 0.20 300 0.06 500 0.10 1000 0.20 
Bảng 4.6: Kết quả tính toán hệ số tiêu sơ bộ của hệ thống theo cơ c u sử d ng đ t: 
 l/ dc = 0,7/0,3 
Ngày 
Pi 
(mm) 
 l ql l.qli dc 
qdc 
(l/s/ha) 
 l.qdci 
qi 
(l/s/ha) 
 
(ha) 
Q 
(m
3
/s) 
 
(ha) 
Q 
(m3/s) 
 
(ha) 
Q 
(m
3
/s) 
1 11.33 0.3 0.07 0.02 0.40 0.85 0.34 0.36 300 0.11 500 0.18 1000 0.36 
2 40.68 0.3 1.26 0.38 0.40 3.06 1.22 1.60 300 0.48 500 0.80 1000 1.60 
3 220.7 0.3 12.84 3.85 0.40 16.60 6.64 10.49 300 3.15 500 5.25 1000 10.49 
4 78.63 0.3 16.81 5.04 0.40 5.91 2.37 7.41 300 2.22 500 3.70 1000 7.41 
5 44.63 0.3 6.63 1.99 0.40 3.36 1.34 3.33 300 1.00 500 1.67 1000 3.33 
6 
0.3 2.71 0.81 0.00 0.00 0.81 300 0.24 500 0.41 1000 0.81 
7 
0.3 0.51 0.15 0.00 0.00 0.15 300 0.05 500 0.08 1000 0.15