Luận án Tiến sĩ Kỹ thuật Nghiên cứu xây dựng cơ sở khoa học đề xuất giải pháp nhằm nâng cao năng lực và hiệu quả của hệ thống tiêu Bắc Nam Hà trong điều kiện biến đổi khí hậu n

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 LA ĐỨC DŨNG NGHIÊN CỨU XÂY DỰNG CƠ SỞ KHOA HỌC ĐỀ XUẤT GIẢI PHÁP NHẰM NÂNG CAO NĂNG LỰC VÀ HIỆU QUẢ CỦA HỆ THỐNG TIÊU BẮC NAM HÀ TRONG ĐIỀU KIỆN BIẾN ĐỔI KHÍ HẬU NƯỚC BIỂN DÂNG LUẬN ÁN TIẾN SĨ KỸ THUẬT HÀ NỘI, NĂM 2017 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 LA ĐỨC DŨNG NGHIÊN CỨU XÂY DỰNG CƠ SỞ KHOA HỌC ĐỀ XUẤT GIẢI PHÁP NHẰM NÂNG CAO NĂNG LỰC VÀ HIỆU QUẢ CỦA HỆ THỐNG TIÊU BẮC NAM HÀ TRONG ĐIỀU KIỆN BIẾN ĐỔI KHÍ HẬU NƯỚC BIỂN DÂNG Chuyên ngành: Mã số: Quy hoạch quản lý tài nguyên nƣớc 62-62-30-01 NGƢỜI HƢỚNG DẪN KHOA HỌC: 1: PGS.TS Nguyễn Tuấn Anh 2: GS TS Đào Xuân Học HÀ NỘI, NĂM 2017 LỜI CAM ĐOAN Tác giả xin cam đoan cơng trình nghiên cứu tác giả, kết nghiên cứu luận án trung thực chƣa đƣợc công bố cơng trình khác Việc tham khảo nguồn tài liệu (nếu có) đƣợc trích dẫn theo quy định Tác giả luận án La Đức Dũng i LỜI CẢM ƠN Sau thời gian miệt mài nghiên cứu, luận án hoàn thành dƣới hƣớng dẫn tận tình GS TS Đào Xuân Học PGS TS Nguyễn Tuấn Anh Tác giả xin bày tỏ lòng biết ơn sâu sắc tới Thầy trực tiếp hƣớng dẫn, bảo tận tình tạo điều kiện thuận lợi suốt trình nghiên cứu để tác giả hoàn thành luận án Tác giả xin chân thành cảm ơn đơn vị, cá nhân tạo điều kiện cung cấp tài liệu tham gia góp ý cho tác giả suốt trình nghiên cứu Xin cảm ơn tác giả cơng trình nghiên cứu cung cấp nguồn tƣ liệu kiến thức liên quan quý báu để tác giả sử dụng q trình nghiên cứu trích dẫn luận án Tác giả xin chân thành cảm ơn Bộ môn Kỹ thuật hạ tầng PTNT; Khoa Kỹ thuật tài nguyên nƣớc; Phòng Đào tạo đại học sau đại học; Trƣờng Đại học Thủy lợi; Tổng cục Thủy lợi; Viện Quy hoạch thủy lợi; Viện Khoa học Thủy lợi Việt Nam; Viện kỹ thuật tài nguyên nƣớc; Công ty TNHH thành viên khai thác công trình thủy lợi Bắc Nam Hà; Trung tâm Khí tƣợng Thủy văn quốc gia; Viện Khoa học Khí tƣợng thủy văn & Biến đổi khí hậu; Cục Khí tƣợng thủy văn & Biến đổi khí hậu tồn thể thầy, cô giáo; bạn bè; đồng nghiệp; quan gia đình tạo điều kiện, chia sẻ khó khăn, tham gia ý kiến giúp đỡ tác giả suốt q trình học tập hồn thành luận án ii MỤC LỤC MỞ ĐẦU 1 Tính cấp thiết đề tài nghiên cứu Mục tiêu nghiên cứu .2 Đối tƣợng phạm vi nghiên cứu Nội dung nghiên cứu Phƣơng pháp nghiên cứu Ý nghĩa khoa học thực tiễn nghiên cứu Bố cục luận án CHƢƠNG TỔNG QUAN 1.1 Tổng quan hệ thống tiêu Bắc Nam Hà 1.1.1 Vị trí địa lý đặc điểm địa hình 1.1.2 Đặc điểm mƣa, mạng lƣới sơng ngịi chế độ thủy văn 1.1.3 Hiện trạng hệ thống tiêu Bắc Nam Hà 10 1.2 Tổng quan tình hình nghiên cứu nƣớc tác động BĐKH, NBD đến hệ thống tiêu nƣớc .14 1.2.1 Tổng quan nghiên cứu giới tác động BĐKH, NBD đến hệ thống tiêu nƣớc 14 1.2.2 Tổng quan nghiên cứu tác động BĐKH, NBD đến hệ thống tiêu nƣớc Việt Nam 16 1.3 Tổng quan tình hình nghiên cứu phƣơng pháp xác định mực nƣớc sông thiết kế cho trạm bơm tiêu 21 1.4 Tổng quan tình hình nghiên cứu quản lý vận hành hệ thống tiêu Bắc Nam Hà 25 1.5 Kết luận chƣơng 27 CHƢƠNG ĐÁNH GIÁ KHẢ NĂNG TIÊU ÚNG CỦA HỆ THỐNG TIÊU BẮC NAM HÀ TRONG ĐIỀU KIỆN BIẾN ĐỔI KHÍ HẬU, NƢỚC BIỂN DÂNG 29 2.1 Đánh giá khả tiêu úng hệ thống tiêu Bắc Nam Hà giai đoạn 29 2.1.1 Xác định nhu cầu tiêu hệ thống 29 2.1.2 Mô hệ thống để kiểm tra lực tiêu 43 iii 2.2 Đánh giá khả tiêu úng hệ thống tiêu Bắc Nam Hà giai đoạn 2030 2050 70 2.2.1 Xác định kịch BĐKH, NBD cho vùng nghiên cứu 70 2.2.2 Lựa chọn tài liệu mƣa để tính tốn cho vùng nghiên cứu 72 2.2.3 Xác định mực nƣớc sông tƣơng lai ứng với kịch phát thải trung bình cho vùng nghiên cứu 73 2.2.4 Xây dựng mơ hình mƣa tiêu điều kiện biến đổi khí hậu 74 2.2.5 Xây dựng kịch sử dụng đất giai đoạn 2030, 2050 81 2.2.6 Xác định hệ số tiêu thiết kế hệ thống Bắc Nam Hà điều kiện BĐKH NBD 81 2.2.7 Mô hệ thống tiêu ứng với trận mƣa thiết kế trạng hệ thống tƣơng ứng với giai đoạn 2030, 2050 82 2.3 Đánh giá khả tiêu úng hệ thống BNH điều kiện BĐKH, NBD 84 2.4 Kết luận chƣơng 84 CHƢƠNG XÂY DỰNG CƠ SỞ KHOA HỌC CHO VIỆC ĐỀ XUẤT GIẢI PHÁP NHẰM NÂNG CAO NĂNG LỰC VÀ HIỆU QUẢ CỦA HỆ THỐNG TIÊU BNH TRONG ĐIỀU KIỆN BĐKH, NBD 86 3.1 Tính tốn cân nƣớc cho hệ thống giai đoạn tại, 2030 2050 86 3.1.1 Phƣơng trình cân nƣớc lƣu vực tiêu 86 3.1.2 Tính tốn cân nƣớc cho hệ thống qua giai đoạn tại, 2030 2050 .87 3.2 Xây dựng quan hệ hệ số tiêu với tỷ lệ diện tích hồ điều hòa, chiều sâu trữ định hƣớng quy hoạch hồ điều hòa cho hệ thống BNH điều kiện BĐKH, NBD 92 3.2.1 Xây dựng quan hệ hệ số tiêu với tỷ lệ diện tích hồ điều hịa chiều sâu trữ cho hệ thống Bắc Nam Hà điều kiện BĐKH NBD 92 3.2.2 Định hƣớng công tác quy hoạch hồ điều hòa cho hệ thống Bắc Nam Hà điều kiện BĐKH NBD 97 3.3 Đề xuất phƣơng pháp xác định tần suất mực nƣớc sông thiết kế trạm bơm tiêu có tổng chi phí xây dựng quản lý vận hành nhỏ 101 3.3.1 Đánh giá mức độ lãng phí điện bơm tiêu thiết kế trạm bơm với tần suất mực nƣớc sông thiết kế 10% 104 3.3.2 Phân tích quan hệ mực nƣớc bể tháo mực nƣớc sông tiêu 111 iv 3.3.3 Ứng dụng lý thuyết phân tích hệ thống để xây dựng phƣơng pháp xác định mực nƣớc sơng thiết kế hợp lý có tổng chi phí xây dựng vận hành trạm bơm nhỏ .113 3.3.4 Thuật toán chƣơng trình tính 119 3.3.5 Ứng dụng phƣơng pháp xác định mực nƣớc sông thiết kế có tổng chi phí xây dựng vận hành quy năm đầu nhỏ cho trạm bơm Hữu Bị 123 3.4 Kết luận chƣơng 138 KẾT LUẬN VÀ KIẾN NGHỊ .139 Kết đạt đƣợc luận án 139 Những đóng góp luận án 141 Những tồn hƣớng nghiên cứu tiếp 141 Kiến nghị 142 DANH MỤC CƠNG TRÌNH ĐÃ CƠNG BỐ 143 TÀI LIỆU THAM KHẢO 144 PHỤ LỤC 149 v DANH MỤC HÌNH ẢNH Hình 1.1 Sơ họa địa giới hành Hệ thống Bắc Nam Hà Hình 2.1 Sơ đồ xác định loại cột nƣớc tiêu qua đỉnh đƣờng tràn 35 Hình 2.2 Sơ đồ tiêu nƣớc qua đƣờng tràn có chế độ chảy tự 36 Hình 2.3 Sơ đồ tiêu nƣớc qua đƣờng tràn có chế độ chảy ngập 37 Hình 2.4 Sơ đồ loại mực nƣớc hồ điều hồ 40 Hình 2.5 Bản đồ phân vùng tiêu hệ thống Bắc Nam Hà 54 Hình 2.6 Các thành phần hệ thống mô SWMM 5.0 56 Hình 2.7 Sơ đồ hệ thống tiêu Bắc Nam Hà SWMM 5.0 58 Hình 2.8 Sơ đồ khối tính tốn dòng chảy cho hệ thống BNH 59 Hình 2.9 Diễn biếnchiều sâu mực nƣớc điểm đo CG19 Kết kiểm tra: Ngày 24-30/09/2009 62 Hình 2.10 Diễn biến chiều sâu mực nƣớc điểm đo NTR Kết kiểm tra: Ngày 24-30/09/2009 63 Hình 2.11 Đƣờng trình mực nƣớc điểm đo NTR- Sơng Nhâm Tràng 66 Hình 2.12 Đƣờng trình mực nƣớc điểm CGI-Sông T5 67 Hình 2.13 Các nút ngập từ kết chạy mơ hình 68 Hình 2.14 Mô đoạn kênh tràn bờ sông Long Xuyên 68 Hình 2.15 Các đoạn kênh tràn bờ từ kết chạy mơ hình 69 Hình 2.16 Quan hệ mƣa năm mƣa ngày max trạm Hƣng Yên (1980-1999) 72 Hình 2.17 Quan hệ mƣa mùa mƣa ngày max trạm Hƣng Yên (1980-1999) 73 Hình 2.18 Xu mƣa 3, 5, ngày Max Trạm Phủ Lý (Giai đoạn 1980-2059) 79 Hình 2.19 Xu mƣa 3, 5, ngày Max Trạm Nam Định (Giai đoạn 1980-2059) 79 Hình 2.20 Xu thời gian ngập qua thời kỳ 83 Hình 3.1 Năng lực nhu cầu tiêu giai đoạn 89 Hình 3.2 Nhu cầu tăng lực bơm tiêu hệ thống qua giai đoạn 89 Hình 3.3 Biểu đồ quan hệ q ~ α ~ Htrữ giai đoạn 2030 94 Hình 3.4 Biểu đồ quan hệ q ~ α ~ Htrữ giai đoạn 2050 96 Hình 3.5 Ảnh phân bố ao hồ tháng 9/2007 – Làng Chều, xã Nguyên Lý, huyện Lý Nhân, tỉnh Hà Nam 99 Hình 3.6 So sánh diện tích ao bị lấp – Làng Chều, xã Nguyên Lý, huyện Lý Nhân, tỉnh Hà Nam 99 vi Hình 3.7 Ảnh phân bố ao tháng 9/2007 – Xóm Già, xã Kim Thái, huyện Vụ Bản, tỉnh Nam Định 100 Hình 3.8 So sánh diện tích ao bị lấp – Xóm Già, xã Kim Thái, huyện Vụ Bản, tỉnh Nam Định 100 Hình 3.9 So sánh tần suất mƣa đồng mực nƣớc sơng trạm Nhƣ Trác 103 Hình 3.10 Cắt dọc bể tháo trạm bơm 108 Hình 3.11: Trạm bơm tiêu làm việc trƣờng hợp mực nƣớc sông tiêu nhỏ mực nƣớc bể tháo trạm bơm Phù Đổng 112 Hình 3.12: Sơ đồ tính cột nƣớc địa hình dƣ thừa 112 Hình 3.13: Sơ đồ khối xác định mực nƣớc bể tháo tối ƣu 119 Hình 3.14: Sơ đồ khối tính tốn điện tiêu thụ trạm bơm tiêu 120 Hình 3.15 Giao diện chƣơng trình tính tốn xác định phƣơng án tối ƣu 122 Hình 3.16 Xu chi phí xây dựng theo mực nƣớc sông thiết kế 128 Hình 3.17 Xu Cđnj ~ ZsTK theo phƣơng án I 131 Hình 3.18: Xu Cđnj ~ ZsTK theo phƣơng án II 131 Hình 3.19: Diễn biến Ctổng theo mực nƣớc sông thiết kế (Phƣơng án I) 133 Hình 3.20: Diễn biến Ctổng theo mực nƣớc sông thiết kế (Phƣơng án II) 133 vii DANH MỤC BẢNG BIỂU Bảng 1.1 Phân bố cao độ ruộng đất khu vực Bắc Nam Hà Bảng 1.2 Các mực nƣớc tháng sông Hồng, sông Đáy, sông Đào Bảng 1.3 Mức báo động số vị trí sơng Bảng 1.4 Thông số kỹ thuật số trạm bơm 10 Bảng 1.5 Quy mô số trạm bơm hệ thống 11 Bảng 1.6 Hệ thống kênh trục cống điều tiết trạm bơm lớn 12 Bảng 2.1 Thời điểm diễn trận mƣa mƣa 3, 5, ngày max tỷ trọng mƣa trạm Ninh Bình giai đoạn 30 Bảng 2.2 Sự phụ thuộc trận mƣa thời đoạn ngắn vào trận mƣa thời đoạn dài trạm Ninh Bình 32 Bảng 2.3 Tổng hợp phụ thuộc trận mƣa thời đoạn ngắn vào trận mƣa thời đoạn dài tỷ trọng trận mƣa giai đoạn 32 Bảng 2.4 Mơ hình mƣa tiêu thiết kế (10%) thời kỳ - Trạm Hƣng Yên 33 Bảng 2.5 Mơ hình mƣa tiêu thiết kế (10%) thời kỳ – Trạm Nam Định 34 Bảng 2.6 Mơ hình mƣa tiêu thiết kế (10%) thời kỳ – Trạm Ninh Bình 34 Bảng 2.7 Mơ hình mƣa tiêu thiết kế (10%) thời kỳ – Trạm Phủ Lý 34 Bảng 2.8 Hệ số dòng chảy C cho đối tƣợng tiêu nƣớc hệ thống thủy lợi theo TCVN 10406: 2015 41 Bảng 2.9 Hệ số tiêu thiết kế theo trạm đo mƣa khống chế giai đoạn 42 Bảng 2.10 Cơng trình đầu mối cống điều tiết khu tiêu 52 Bảng 2.11 Bảng đánh giá sai số q trình tính tốn thực đo node CG19-Mỹ Hà 64 Bảng 2.12 Bảng đánh giá sai số trình tính tốn thực đo NTR - Nhâm Tràng 64 Bảng 2.13 Bảng đánh giá sai số q trình tính tốn thực đo vị trí sơng Nhƣ Trác, Sắt, Châu Giang, Nhâm Tràng 66 Bảng 2.14 Bảng đánh giá sai số trình tính tốn thực đo vị trí sông Sắt, sông Châu Giang, sông T5 66 Bảng 2.15 So sánh mực NBD kịch B2 RCP4.5 giai đoạn 2030 2050 74 Bảng 2.16 Trận mƣa 3, 5,7 ngày max trạm Ninh Bình thời kỳ 2030 75 viii End Sub (*Chương trình con: Tính điện năng*) Private Sub cmdChay_CT_Click() Dim Tenfile1$, Tenfile2$, Tenfile3$, Tenfile4$, Tenfile5$, Tenfile$, Tenthumuc$, st$, tentran$ Dim i%, j%, jj%, k%, kk%, n1%, n2%, n3%, N%, sonam% Dim x!, delt!, Tmua!, Tdongchay!, Tsaumua!, qr_max! Dim qr1!, Hp!, qmax!, qmax1!, Qtongmax! Dim Hc1!, Hc2! Dim Ebtran!(100, 50) Dim nam$(100), sotran%(100), Ebnam!(100) 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! Dim dientichlua!, dientichmau!, dientichthocu!, dientichctcc!, dientichao_kodt! Dim dientichao_codt!, Cmau!, Cthocu!, Cctcc!, htru!, Wtru!, ik!, nk! Dim ndiem%, somay%, somaychay% Dim Hdactinh!(20), Qdactinh!(20), hsdactinh!(20), hdhdactinh!(20) Dim S_ong!, hsdc! Dim h_kenhhut!, b_kenhhut!, m_kenhhut!, n_kenhhut!, i_kenhhut!, Zdaykenhhut! Dim b_kenhthao!, h_kenhthao!, m_kenhthao!, n_kenhthao!, i_kenhthao!, htt_kenhthao!, Zdaykenhthao! Dim chuvi_uot!, dientich_uot!, Eb!, Ebtong! Const ep = 0.00001 Const beta = 1.03 Const hs = 0.1 Const mt = 0.34 Const hsm = 0.9 Print "Dang chay chuong trinh " ' Nhap du lieu ve tram bom, luu vuc tieu Tenthumuc = "D:\TDN" Tenfile = "SL_TB" Tenfile1 = Tenthumuc + "\" + Tenfile + ".txt" Open Tenfile1 For Input As 'Open "D:\SL_TB.txt" For Input As Input #1, st$ Input #1, st$ Input #1, ndiem Input #1, st$ For i% = To ndiem Input #1, Hdactinh(i) Next Input #1, st$ For i% = To ndiem Input #1, Qdactinh(i) Next Input #1, st$ For i% = To ndiem Input #1, hsdactinh(i) Next Input #1, st$ PL2.1-4 Input #1, somay Input #1, st$ Input #1, S_ong 'He so ma sat duong ong Input #1, st$ Input #1, hsdc 'Hieu suat dong co Input #1, st$ 'So lieu dien tich luu vuc (ha) Input #1, st$ Input #1, dientichlua Input #1, st$ Input #1, dientichmau Input #1, st$ Input #1, dientichthocu Input #1, st$ Input #1, dientichctcc Input #1, st$ Input #1, dientichao_kodt Input #1, st$ Input #1, dientichao_codt Input #1, st$ Input #1, b 'chieu rong tran mat ruong Input #1, st$ Input #1, Cmau Input #1, st$ Input #1, Cthocu Input #1, st$ Input #1, Cctcc Input #1, st$ Input #1, htru 'met Input #1, st$ 'So lieu kenh hut, thao Input #1, b_kenhhut Input #1, st$ Input #1, m_kenhhut Input #1, st$ Input #1, n_kenhhut Input #1, st$ Input #1, i_kenhhut Input #1, st$ Input #1, Zdaykenhhut Input #1, st$ Input #1, b_kenhthao Input #1, st$ Input #1, m_kenhthao Input #1, st$ Input #1, n_kenhthao Input #1, st$ Input #1, i_kenhthao Input #1, st$ Input #1, Zdaykenhthao Input #1, st$ Input #1, httkthao Close Cls ' Nhap so lieu So tran mua cac nam PL2.1-5 'Tenfile2 = InputBox("Nhap ten file du lieu so nam, so tran mua(VD: D:\TDN\Sotran.txt): ") Tenfile2 = Tenthumuc + "\" + "Sotran.txt" Open Tenfile2 For Input As Input #2, st$ Input #2, st$ Input #2, sonam Input #2, st$ For j = To sonam Input #2, nam(j) Next Input #2, st$ For j = To sonam Input #2, sotran(j) Next Close Tenfile = InputBox("Nhap ten file ket qua (VD: KQ_Zdk25") Tenfile5 = Tenthumuc + "\" + Tenfile + ".txt" Open Tenfile5 For Append As ' Bat dau doc SL va tinh cho tung tran mua jj = Do Until jj > sonam kk = Do Until kk > sotran(jj) ' -tentran = LTrim(Str(kk)) Tenfile3 = Tenthumuc + "\" + nam(jj) + "\" + tentran + ".txt" Open Tenfile3 For Input As Input #3, st$ Input #3, st$ Input #3, delt 'gio Input #3, st$ Input #3, Tsaumua Input #3, st$ Input #3, n1 'so thoi doan mua Input #3, st$ For i = To n1 Input #3, P(i) 'mm/thoi doan Next Input #3, st$ Input #3, n2 'so thoi doan muc nuoc Input #3, st$ For i = To n2 Input #3, Zsong(i) Next Close 'Print "Da doc xong" Tenfile4 = Tenthumuc + "\" + nam(jj) + "\T" + tentran + ".txt" Open Tenfile4 For Input As Input #4, st$ Input #4, st$ Input #4, n3 'so thoi doan tieu Input #4, st$ PL2.1-6 For i = To n3 Input #4, qr(i) 'he so tieu Next Close '(* *) Tmua = n1 * delt 'Thoi gian mua, gio Tdongchay = Tmua + Tsaumua 'Thoi gian mo phong, gio N = Round(Tdongchay / delt) For i = To n1 Pr(i) = P(i) Next For i = (n1 + 1) To N Pr(i) = Next ' -Tinh luu luong tieu Qt cua ca luu vuc Qtongmax = For i = To N Qlua(i) = qr(i) * dientichlua / 1000 Qmau(i) = Cmau * (Pr(i) / (delt * 360)) * dientichmau Qthocu(i) = Cthocu * (Pr(i) / (delt * 360)) * dientichthocu Qctcc(i) = Cctcc * (Pr(i) / (delt * 360)) * dientichctcc Qao_kodt(i) = (Pr(i) / (delt * 360)) * dientichao_kodt Qao_codt(i) = (Pr(i) / (delt * 360)) * dientichao_codt Qtong(i) = Qlua(i) + Qmau(i) + Qthocu(i) + Qctcc(i) + Qao_kodt(i) + Qao_codt(i) If Qtong(i) > Qtongmax Then Qtongmax = Qtong(i) Tmax = i End If Next Wtru = dientichao_codt * 10000 * htru Qtong(Tmax) = Qtongmax - Wtru / (delt * 3600) Ebtong = ' -Tinh dien nang cho tung thoi doan i For i = To N ' Tinh hkenh hut va muc nuoc be hut H = 0.02 b = b_kenhhut m = m_kenhhut nk = n_kenhhut ik = i_kenhhut dientich_uot = (b + m * H) * H chuvi_uot = * H * Sqr(1 + m * m) + b Q = dientich_uot * (1 / nk) * Sqr(ik) * (dientich_uot / chuvi_uot) ^ (2 / 3) Do Until Q > Qtong(i) H = H + 0.01 dientich_uot = (b + m * H) * H chuvi_uot = * H * Sqr(1 + m * m) + b Q = dientich_uot * (1 / nk) * Sqr(ik) * (dientich_uot / chuvi_uot) ^ (2 / 3) Loop h_kenhhut = H Zbhut = Zdaykenhhut + h_kenhhut ' -Tinh hkenh thao va muc nuoc be thao PL2.1-7 H = 0.02 b = b_kenhthao m = m_kenhthao nk = n_kenhthao ik = i_kenhthao dientich_uot = (b + m * H) * H chuvi_uot = * H * Sqr(1 + m * m) + b Q = dientich_uot * (1 / nk) * Sqr(ik) * (dientich_uot / chuvi_uot) ^ (2 / 3) Do Until Q > Qtong(i) H = H + 0.01 dientich_uot = (b + m * H) * H chuvi_uot = * H * Sqr(1 + m * m) + b Q = dientich_uot * (1 / nk) * Sqr(ik) * (dientich_uot / chuvi_uot) ^ (2 / 3) Loop h_kenhthao = H Zbthao1 = Zdaykenhthao + h_kenhthao Zbthao2 = Zsong(i) + httkenhthao If Zbthao1 > Zbthao2 Then Zbthao = Zbthao1 Else Zbthao = Zbthao2 End If ' -Tinh hdh, Hbom, Q1m, somaychay, N, Ebom Hdh = Zbthao - Zbhut For k = To ndiem hdhdactinh(k) = Hdactinh(k) - S_ong * Qdactinh(k) * Qdactinh(k) Next 'Tu Hdh tra duong dac tinh Hdh-Q de tim Q1m If Hdh > hdhdactinh(1) Then Q1m = Else If Hdh < hdhdactinh(ndiem) Then Q1m = Else For k = To ndiem If (Hdh >= hdhdactinh(k)) And (Hdh = Qdactinh(k - 1)) And (Q1m Ctong(i) Then Cmin = Ctong(i) Pan = i PL2.1-10 End If Next Zstoiuu = Zstk(Pan) Print " KET QUA TINH TOAN " Print "Muc nuoc song thiet ke toi uu, Zs toi uu = ", Format(Zstoiuu, "###.##"), "m" Print "Tong chi phi nho nhat tuong ung Zs toi uu, Cmin = ", Format(Cmin, "#####.###"), "ty dong" End Sub ‘ -Hết chương trình PL2.1-11 Bảng PL 2.3 Kinh phí xây dựng cơng trình tháo theo phương án mực nước STT 4 Công việc Đơn vị Bể xả Đào đất máy Đắp đất cơng trình Bê tông M200 Thép loại Kênh xả cống Đào đất máy Mua đất để đắp (phần thiếu) Đắp đất cơng trình Bê tơng loại Thép loại Cống qua đê Đào đất máy Đắp đất công trình Bê tơng loại Thép loại Kênh xả ngồi cống Đào đất máy Đắp đất cơng trình Cộng m3 m3 m3 kg Khối lượng MN+5,79 MN+5,29 MN+4,79 MN+4,29 MN+3,79 MN+3,29 MN+2,79 12160 13631 15102 16572 18043 19513 7025 7878 8787 10192 10876 11937 2233 2288 2342 2396 2451 2505 167498 171563 175613 179663 183826 187876 MN+5,79 MN+5,29 MN+4,79 Thành tiền MN+4,29 MN+3,79 MN+3,29 MN+2,79 20984 17.000 13046 13.500 2559 1.700.000 191942 22.500 206.724.950 94.841.010 3.796.611.700 3.768.695.438 231.725.660 106.353.810 3.888.768.700 3.860.174.813 256.726.370 118.621.260 3.980.568.700 3.951.299.813 281.727.080 137.592.810 4.072.368.700 4.042.424.813 306.727.790 146.824.110 4.166.718.700 4.136.081.063 331.728.500 161.144.910 4.258.518.700 4.227.206.063 356.729.210 176.115.060 4.350.675.700 4.318.685.438 18464 17.000 11757 50.000 26744 13.500 1538 1.700.000 115350 22.500 90.712.000 1.527.626.000 353.646.000 1.931.200.000 1.917.000.000 119.068.000 1.449.158.000 354.618.000 2.046.800.000 2.031.750.000 150.178.000 1.360.398.000 355.158.000 2.159.000.000 2.143.125.000 186.150.000 1.262.818.000 356.778.000 2.274.600.000 2.257.875.000 226.372.000 837.100.000 359.100.000 2.386.800.000 2.369.250.000 268.736.000 718.376.000 360.504.000 2.502.400.000 2.484.000.000 313.888.000 587.836.000 361.044.000 2.614.600.000 2.595.375.000 60.143.076 54.200.664 1.581.081.600 1.569.456.000 67.905.684 66.359.412 1.581.081.600 1.569.456.000 76.015.296 79.035.264 1.581.081.600 1.569.456.000 84.471.912 104.339.340 1.581.081.600 1.569.456.000 93.271.248 105.911.064 1.581.081.600 1.569.456.000 102.426.156 120.111.012 1.581.081.600 1.569.456.000 111.906.648 134.828.064 1.581.081.600 1.569.456.000 Đơn giá m3 m3 m3 m3 kg 5336 30553 26196 1136 85200 7004 28983 26268 1204 90300 8834 10950 13316 15808 27208 25256 16742 14368 26308 26428 26600 26704 1270 1338 1404 1472 95250 100350 105300 110400 m3 m3 m3 kg 3538 4015 930 69754 3994 4916 930 69754 4471 5854 930 69754 4969 7729 930 69754 5487 7845 930 69754 6025 8897 930 69754 6583 17.000 9987 13.500 930 1.700.000 69754 22.500 m3 m3 8796 11052 11688 11052 14877 11052 18366 11052 22158 11052 26247 11052 30633 11052 17.000 13.500 149.532.000 198.696.000 252.909.000 312.222.000 376.686.000 446.199.000 520.761.000 149.202.000 149.202.000 149.202.000 149.202.000 149.202.000 149.202.000 149.202.000 17.250.672.438 17.721.117.679 18.182.774.303 18.673.107.255 18.830.581.575 19.281.089.941 19.742.183.720 Pl 2.3-1 Bảng PL 2.4 Kinh phí xây dựng cơng trình tháo theo phương án mực nước STT 4 Công việc Bể xả Đào đất máy Đắp đất cơng trình Bê tơng M200 Thép loại Kênh xả cống Đào đất máy Mua đất để đắp (phần thiếu) Đắp đất cơng trình Bê tơng loại Thép loại Cống qua đê Đào đất máy Đắp đất cơng trình Bê tơng loại Thép loại Kênh xả cống Đào đất máy Đắp đất cơng trình Cộng Đơn vị Khối lượng MN+5,54 MN+5,04 MN+4,54 MN+4,04 MN+3,54 MN+3,04 Đơn giá MN+5,54 MN+5,04 Thành tiền MN+4,54 MN+4,04 MN+3,54 MN+3,04 m3 m3 m3 kg 12896 7452 2260 169530 14366 8332 2315 173588 15837 9489 2369 177638 17307 10534 2423 181745 18778 11406 2478 185851 20249 12491 2532 189909 17.000 13.500 1.700.000 22.500 219.225.305 100.597.410 3.842.690.200 3.814.435.125 244.226.015 112.487.535 3.934.668.700 3.905.737.313 269.226.725 128.107.035 4.026.468.700 3.996.862.313 294.227.435 142.208.460 4.119.543.700 4.089.252.938 319.228.145 153.984.510 4.212.618.700 4.181.643.563 344.228.855 168.629.985 4.304.597.200 4.272.945.750 m3 m3 m3 m3 kg 6158 29807 26252 1170 87750 7894 28121 26288 1236 92700 9824 26289 26360 1304 97800 12118 22358 26520 1370 102750 14542 15579 26656 1438 107850 17098 13100 26724 1506 112950 17.000 50.000 13.500 1.700.000 22.500 104.686.000 1.490.362.000 354.402.000 1.989.000.000 1.974.375.000 134.198.000 1.406.028.000 354.888.000 2.101.200.000 2.085.750.000 167.008.000 1.314.460.000 355.860.000 2.216.800.000 2.200.500.000 206.006.000 1.117.900.000 358.020.000 2.329.000.000 2.311.875.000 247.214.000 778.964.000 359.856.000 2.444.600.000 2.426.625.000 290.666.000 655.006.000 360.774.000 2.560.200.000 2.541.375.000 m3 m3 m3 kg 3538 4015 930 69754 3994 4916 930 69754 4471 5854 930 69754 4969 7729 930 69754 5487 7845 930 69754 6583 9987 930 69754 17.000 13.500 1.700.000 22.500 60.143.076 54.200.664 1.581.081.600 1.569.456.000 67.905.684 66.359.412 1.581.081.600 1.569.456.000 76.015.296 79.035.264 1.581.081.600 1.569.456.000 84.471.912 104.339.340 1.581.081.600 1.569.456.000 93.271.248 105.911.064 1.581.081.600 1.569.456.000 111.906.648 134.828.064 1.581.081.600 1.569.456.000 m3 m3 8796 11052 11688 11052 14877 11052 18366 11052 22158 11052 30633 11052 17.000 13.500 149.532.000 149.202.000 17.453.388.380 198.696.000 149.202.000 17.911.884.259 252.909.000 149.202.000 18.382.991.933 312.222.000 149.202.000 18.768.806.385 376.686.000 149.202.000 19.000.341.830 520.761.000 149.202.000 19.565.658.102 Pl 2.4-1