Luận án nghiên cứu xây dựng mô hình hiện đại hóa quản lý vận hành hệ thống tưới

LỜI CAM ĐOAN Tác giả luận án xin cam đoan cơng trình nghiên cứu thân tác giả Các kết nghiên cứu kết luận luận án trung thực, không chép từ nguồn dƣới hình thức Việc tham khảo, trích dẫn nguồn tài liệu đƣợc ghi rõ nguồn tài liệu tham khảo quy định Tác giả luận án Ngô Đăng Hải i LỜI CẢM ƠN Tác giả luận án xin trân trọng cảm ơn GS.TS Bùi Hiếu PGS.TS Lê Văn Ƣớc tận tình hƣớng dẫn tác giả nghiên cứu hoàn thành Luận án tiến sĩ Tác giả xin chân thành cảm ơn Ban giám hiệu Trƣờng Đại học Thủy lợi, Phòng Đào tạo Đại học Sau Đại học, Khoa Kỹ thuật tài nguyên nƣớc, Bộ môn Kỹ thuật tài nguyên nƣớc tạo điều kiện thuận lợi, nhiệt tình giúp đỡ động viên tác giả trình nghiên cứu, thực luận án Tác giả luận án xin chân thành cảm ơn đồng nghiệp, bạn bè, gia đình hỗ trợ, động viên, chia sẻ khó khăn suốt thời gian nghiên cứu, làm luận án ii MỤC LỤC LỜI CAM ĐOAN .i LỜI CẢM ƠN ii MỤC LỤC iii DANH MỤC CÁC HÌNH VẼ, ĐỒ THỊ vii DANH MỤC CÁC BẢNG BIỂU ix DANH MỤC CÁC TỪ VIẾT TẮT x MỞ ĐẦU 1 Tính cấp thiết đề tài luận án Mục tiêu nghiên cứu Đối tƣợng phạm vi nghiên cứu .4 Ý nghĩa khoa học thực tiễn .5 Những đóng góp luận án .6 Cấu trúc luận án CHƢƠNG 1.1 TỔNG QUAN VỀ VẤN ĐỀ NGHIÊN CỨU .8 Tổng quan tình hình nghiên cứu giới 1.1.1 Các cơng trình nghiên cứu mơ hình lập kế hoạch quản lý vận hành giám sát, điều khiển bán tự động .9 1.1.1.1 1.1.1.2 1.1.1.3 1.1.1.4 Khái quát chung Ưu điểm 13 Hạn chế 13 Phạm vi áp dụng khả ứng dụng vào Việt Nam 14  Nhận xét đánh giá 15 1.1.2 Các công trình nghiên cứu mơ hình giám sát điều khiển tự động từ xa (SCADA) 16 1.1.2.1 1.1.2.2 1.1.2.3 1.1.2.4 Khái quát chung 16 Ưu điểm 18 Hạn chế 19 Phạm vi áp dụng khả ứng dụng vào Việt Nam 19  Nhận xét đánh giá 20 1.2 Tổng quan tình hình nghiên cứu Việt Nam 21 1.2.1 Các công trình nghiên cứu mơ hình lập kế hoạch quản lý vận hành giám sát, điều khiển bán tự động……… .21 1.2.1.1 1.2.1.2 1.2.1.3 Quản lý điều hành HTTN Đan Hoài (QLĐH Đan Hoài) 21 Hệ thống bán tự động quản lý điều hành HTTN Thạch Nham 22 Quản lý điều hành HTTN Bắc Hưng Hải (QLĐH Bắc Hưng Hải) 23 iii 1.2.1.4 Nghiên cứu cải tiến mơ hình IMSOP để điều hành phân phối nước hợp lý, nâng cao hiệu sử dụng nước hệ thống tưới bơm đồng sơng Hồng (Mơ hình IMSOP cải tiến) 24 1.2.2 xa Các cơng trình nghiên cứu mơ hình giám sát điều khiển tự động từ 24 1.2.2.1 Nghiên cứu áp dụng khoa học công nghệ tiên tiến bước đại hố cơng tác quản lý điều hành hệ thống thuỷ nơng Phù Sa, Hà Tây (Mơ hình SCADA/MAC) 24 1.2.2.2 Nghiên cứu áp dụng khoa học công nghệ tiên tiến quản lý vận hành hệ thống thuỷ nơng Ấp Bắc - Nam Hồng (Mơ hình VKHTLMB SCADA) 26 1.2.2.3 Nghiên cứu xây dựng hệ thống thông tin quản lý điều hành hệ thống thuỷ lợi Gị Cơng (Mơ hình VKHTLMN SCADA) 27 1.2.2.4 Báo cáo thiết kế chi tiết SCADA n Lập (Mơ hình BCEOM1 SCADA) 28 1.2.2.5 Báo cáo thiết kế chi tiết SCADA Cầu Sơn - Cấm Sơn Kẻ Gỗ (Mơ hình BCEOM2 SCADA) 29 1.2.2.6 Hiện đại hóa hệ thống tưới Phú Ninh Đá Bàn (Mơ hình HASKONING SCADA) 30 1.2.2.7 Hệ thống SCADA cho kênh Dầu Tiếng hệ thống tưới Củ Chi (Mơ hình BRLI SCADA) 30 * Nhận xét đánh giá: 31 1.3 Định hƣớng nghiên cứu xây dựng, phát triển mơ hình đại hóa quản lý vận hành hệ thống tƣới 32 Kết luận chƣơng 35 CHƢƠNG NỘI DUNG VÀ PHƢƠNG PHÁP NGHIÊN CỨU 37 2.1 Nội dung nghiên cứu .37 2.2 Phƣơng pháp nghiên cứu 45 2.2.1 Phƣơng pháp mơ hình hóa .45 2.2.1.1 2.2.1.2 thực tế Mơ hình hóa tốn quản lý khai thác “tối ưu“ hệ thống tưới 45 Mơ hình hóa toán quản lý vận hành hệ thống tưới theo số liệu quan trắc 53 2.2.2 Phƣơng pháp dự báo thích nghi yếu tố khí tƣợng, thuỷ văn .57 2.2.3 Phƣơng pháp lập trình .60 2.2.4 Phƣơng pháp thí nghiệm, thực nghiệm 60 Kết luận chƣơng 61 CHƢƠNG KẾT QUẢ NGHIÊN CỨU VÀ THẢO LUẬN 62 3.1 Kết nghiên cứu xây dựng mơ hình đại hóa quản lý vận hành hệ thống tƣới 62 3.1.1 tƣới Mơ hình hóa tốn lập kế hoạch quản lý khai thác ”tối ƣu” hệ thống 62 3.1.1.1 3.1.1.2 3.1.1.3 3.1.1.4 Mơ hình hóa hàm mục tiêu 62 Các điều kiện ràng buộc 66 Các biến (tham số) định 68 Các kết đầu 68 iv 3.1.1.5 Kết nghiên cứu phát triển thuật toán lập kế hoạch quản lý khai thác “tối ưu“ hệ thống tưới 69 3.1.1.6 Kết xây dựng phần mềm lập kế hoạch quản lý khai thác “tối ưu“ hệ thống tưới 73 * Nhận xét đánh giá: 75 3.1.2 Kết nghiên cứu xây dựng mơ hình quản lý vận hành hệ thống tƣới theo số liệu quan trắc (thời gian thực) 75 3.1.2.1 thực tế 3.1.2.2 xa 3.1.2.3 trắc Kết nghiên cứu xây dựng công cụ giám sát, cập nhật số liệu quan trắc 75 Kết nghiên cứu xây dựng công cụ điều khiển hệ thống tưới tự động từ 79 Kết xây dựng phần mềm quản lý vận hành hệ thống tưới theo số liệu quan 81 3.1.3 Kết nghiên cứu đại hóa quản lý vận hành giám sát, điều khiển hệ thống tƣới thông qua mạng Internet hệ thống viễn thơng di động tồn cầu 82 3.1.3.1 Trang web đại hóa quản lý vận hành hệ thống tưới 82 3.1.3.2 Cơ sở liệu trực tuyến 84 3.1.3.3 Kết xây dựng giao diện giám sát, điều khiển phần mềm đại hóa quản lý vận hành hệ thống tưới 86 3.1.3.4 Kết nghiên cứu bước đại hóa giám sát, điều khiển cho HTT quản lý vận hành thủ công giai đoạn độ tiến tới tự động hóa 89 3.1.4 Kết nghiên cứu ứng dụng, phát triển công nghệ giám sát điều khiển hệ thống tƣới 91 3.1.4.1 Kết phát triển, ứng dụng công nghệ giám sát điều khiển HTT theo hướng truyền thông đại lựa chọn thiết bị GS&ĐK kèm hợp lý 91 3.1.4.2 Kết nghiên cứu xây dựng cơng cụ thực hiện đại hóa quản lý vận hành miễn phí cho cơng ty QLKT CTTL 93 3.1.5 Kết nghiên cứu xây dựng sở liệu trực tuyến phục vụ công tác quản lý công trình hệ thống tƣới 94 3.1.6 Kết nghiên cứu cải tiến, phát triển công cụ đánh giá hiệu quản lý khai thác hỗ trợ quy trình nâng cao hiệu hệ thống tƣới 96 3.1.6.1 Cải tiến, phát triển số đánh giá hiệu quản lý khai thác HTT 96 3.1.6.2 Xây dựng công cụ đánh giá số hiệu hỗ trợ nâng cao hiệu quản lý khai thác HTT 100 3.1.7 Kết tổng hợp liên kết thành phần mô hình đại hóa quản lý vận hành hệ thống tƣới 102 3.2 Kết ứng dụng mơ hình đại hóa quản lý vận hành hệ thống tƣới 104 3.2.1 Địa điểm ứng dụng mơ hình 104 3.2.2 Kết lập kế hoạch quản lý khai thác ”tối ƣu” hệ thống tƣới Phù Sa 108 3.2.2.1 Các số liệu đầu vào 108 3.2.2.2 Kết ứng dụng mơ hình lập kế hoạch quản lý khai thác "tối ưu“ HTT Phù Sa 115 v 3.2.3 Kết ứng dụng mơ hình quản lý vận hành hệ thống tƣới theo số liệu quan trắc thực tế 125 3.2.3.1 Kết tính tốn kiểm định mơ hình quản lý vận hành hệ thống tưới số liệu quan trắc năm 2011 125 3.2.3.2 Kết ứng dụng mơ hình quản lý vận hành hệ thống tưới theo số liệu quan trắc thực tế 126 3.2.4 Kết ứng dụng cơng nghệ đại hóa quản lý vận hành, GS&ĐK hệ thống tƣới thông qua mạng Internet hệ thống viễn thơng di động tồn cầu 129 3.2.4.1 Kết ứng dụng website đại hóa quản lý vận hành hệ thống tưới http://www.hiendaihoaqlvh.com 129 3.2.4.2 Kết ứng dụng giao diện giám sát, điều khiển phần mềm đại hóa quản lý vận hành hệ thống tưới 131 3.2.5 tƣới Kết ứng dụng, phát triển công nghệ giám sát điều khiển hệ thống .132 3.2.6 Kết xây dựng sở liệu trực tuyến phục vụ công tác quản lý cơng trình hệ thống tƣới Phù Sa 135 3.2.7 Kết ứng dụng công cụ đánh giá số hiệu hỗ trợ nâng cao hiệu quản lý khai thác HTT Phù Sa .136 Kết luận chƣơng .137 KẾT LUẬN VÀ KIẾN NGHỊ .138 Kết luận 138 Kiến nghị 140 DANH MỤC CÁC CƠNG TRÌNH ĐÃ CƠNG BỐ KẾT QUẢ NGHIÊN CỨU CỦA ĐỀ TÀI LUẬN ÁN .142 TÀI LIỆU THAM KHẢO 143 PHỤ LỤC 148 Phụ lục 1: Thông tin chi tiết mơ hình đại hóa QLVH HTT .148 Phụ lục 2: Các số liệu khí tƣợng, thủy văn; chế độ tƣới tài liệu liên quan đến quản lý khai thác Công ty TNHH MTV Thủy lợi Sơng Tích 156 Phụ lục 3: Kết tính tốn dự báo khí tƣợng, thủy văn theo mơ hình ARIMA .175 Phụ lục 4: Các kết tính tốn chi tiết lập kế hoạch khai thác “tối ƣu” HTT .180 Phụ lục 5: Những sơ đồ thuật toán chƣơng trình phần mềm 194 vi DANH MỤC CÁC HÌNH VẼ, ĐỒ THỊ Hình 1.1: Diễn biến lớp nƣớc mặt ruộng độ ẩm đất đƣợc quan trắc từ vệ tinh 17 Hình 1.2: FarmConnect giám sát thông báo độ ẩm đất từ vệ tinh 18 Hình 1.3: Cấu trúc tổng thể mơ hình đại hóa quản lý vận hành HTT 34 Hình 2.1: HTT với cấu trúc phân cấp, đa mục tiêu nhiều ngƣời định 48 Hình 2.2: Mơ hình cấu trúc tổng qt toán quản lý khai thác HTT 49 Hình 2.3: Cấu trúc HTT mối quan hệ hữu thành phần 50 Hình 2.4: Sơ đồ khối xây dựng dự báo theo mơ hình ARIMA(p, d, q) 59 Hình 3.1: Tổng mức tƣới thay đổi dịch chuyển thời vụ gieo trồng .65 Hình 3.2: Thuật tốn tiến hóa sai phân đa mục tiêu mở rộng .72 Hình 3.3: Giao diện khai báo số biến định toán .73 Hình 3.4: Giao diện thiết lập tốn “tối ƣu“ nhập số liệu đầu vào 73 Hình 3.5: Giao diện tính tốn biến định T0_c,v, Sc,v, Mc,v chế độ tƣới .74 Hình 3.6: Giao diện tính tốn phƣơng án tối ƣu Pareto 74 Hình 3.7: Sơ đồ thuật toán cập nhập số liệu quan trắc 76 Hình 3.8: Mơ hình quản lý vận hành hệ thống tƣới theo số liệu quan trắc 77 Hình 3.9: Giao diện website dùng cho cập nhật liệu quan trắc thủ cơng 78 Hình 3.10: Sơ đồ thuật tốn tính toán quản lý vận hành HTT theo số liệu quan trắc 80 Hình 3.11: Giao diện trang web http://www.hiendaihoaqlvh.com 82 Hình 3.12: Giao diện trình đơn trang “Mơ hình HĐH quản lý vận hành” 83 Hình 3.13: Trang cập nhật yếu tố khí tƣợng từ chuyên mục Dự báo thời tiết 84 Hình 3.14: Cấu trúc sở liệu đại hóa quản lý vận hành HTT 84 Hình 3.15: Cơ sở liệu OnMs vận hành cơng trình .85 Hình 3.16: Cập nhật, lƣu trữ truy xuất liệu trang chủ 85 Hình 3.17: Giao diện khai báo số số cần giám sát số tham số cần điều khiển 86 Hình 3.18: Giao diện nhập tên số cần giám sát tham số cần điều khiển 87 Hình 3.19: Sơ đồ thuật giải phần mềm HĐHQLVH .88 Hình 3.20: Mơ hình quản lý vận hành hệ thống tƣới theo số liệu quan trắc 88 Hình 3.21: Giao diện trang web “Giám sát Điều khiển“ 90 Hình 3.22: Phƣơng thức cập nhật số liệu khí tƣợng, thủy văn qua mạng Internet 90 Hình 3.23: Các chức vi xử lý modem IP F2164 3G RTU 91 Hình 3.24: Cấu trúc mạng GS&ĐK sử dụng modem F2X64 F2X14 .91 Hình 3.25: Thiết lập hệ thống giám sát điều khiển 93 Hình 3.26: Cấu trúc sở liệu thông số quản lý cơng trình 95 Hình 3.27: Tính tốn số hiệu so sánh với tiêu phấn đấu 101 Hình 3.28: Mơ hình đại hóa quản lý vận hành hệ thống tƣới 103 Hình 3.29: Bản đồ hệ thống tƣới Phù Sa .105 Hình 3.30: Giao diện kết mơ hình hóa, thiết lập tốn lập kế hoạch .115 Hình 3.31: Tính tốn chế độ tƣới cho loại trồng 117 vii Hình 3.32: Kết tính tốn biến định T0_c,v, Sc,v, Mc,v , chế độ tƣới .118 Hình 3.33: Kết tính tốn xác định phƣơng án tối ƣu Pareto 119 Hình 3.34: Kết tính toán xác định biến định “tối ƣu” 121 Hình 3.35: Cảm biến đo mực nƣớc sông Hồng trƣớc cống Phù Sa 126 Hình 3.36: Cảm biến đo mực nƣớc bể xả trạm bơm Phù Sa .127 Hình 3.37: Thiết bị điều khiển IP Modem - RTU F2164 trạm bơm Phù Sa 127 Hình 3.38: Đồ thị cập nhật trình lớp nƣớc mặt ruộng độ ẩm đất HTT Phù Sa .128 Hình 3.39: Website http://www.hiendaihoaqlvh.com cung cấp, phổ biến thông tin, tài liệu phục vụ công tác đại hóa tƣới tồn quốc 130 Hình 3.40: Cơng cụ hỗ trợ triển khai thực cơng tác đại hóa 130 Hình 3.41: Giao diện hƣớng dẫn triển khai, thực cơng tác đại hóa QLVH 131 Hình 3.42: Giao diện QLVH giám sát, điều khiển HTT Website đại hóa http://www.hiendaihoaqlvh.com 132 Hình 3.43: Giao diện quản lý cơng trình trực tuyến (Bản đồ số hóa sở liệu) .135 Hình 3.44: Cơng cụ theo dõi lịch tu, sửa chữa định kỳ nhắc việc tự động 135 Hình 3.45: Kết ứng dụng cơng cụ đánh giá số hiệu .136 viii DANH MỤC CÁC BẢNG BIỂU Bảng 1.1: Phân loại tính mơ hình đại hóa quản lý vận hành HTT đƣợc áp dụng phổ biến giới 10 Bảng 1.2: Tổng hợp ƣu, nhƣợc điểm mơ hình đại hóa QLVH hƣớng phát huy, khắc phục 33 Bảng 3.1: 20 số hiệu đƣợc lựa chọn sửa đổi từ số “BenchMarking in Irrigation and Drainage Sector“ 97 Bảng 3.2: Bộ số (mở rộng) đánh giá hiệu quản lý khai thác HTT 99 Bảng 3.3: Nhiệt độ khơng khí khu vực Phù Sa (oC) 106 Bảng 3.4: Độ ẩm khơng khí khu vực Phù Sa (%) .106 Bảng 3.5: Lƣợng bốc đo ống Piche (mm) 106 Bảng 3.6: Thống kê yếu tố khí tƣợng mực nƣớc theo năm 109 Bảng 3.7: Kết xác định tham số p, q thơng số cho mơ hình .110 Bảng 3.8: Lịch lấy nƣớc (xả nƣớc) phục vụ gieo cấy vụ Đông - Xuân 2013 – 2014 113 Bảng 3.9: Các tài liệu nông nghiệp 113 Bảng 3.10: Kết tính tốn biến định “tối ƣu” thời vụ gieo trồng, .120 Bảng 3.11: Các hàm mục tiêu theo phƣơng án tối ƣu Pareto 121 Bảng 3.12: So sánh, đánh giá kết tính tốn số liệu QLVH 122 Bảng 3.13: Quá trình lƣu lƣợng yêu cầu cho vụ chiêm 2014 đầu HTT Phù Sa .123 Bảng 3.14: Kế hoạch vận hành HTT Phù Sa vụ chiêm 2014 124 Bảng 3.15: Thời gian sinh trƣởng công thức tƣới tăng sản 125 Bảng 3.16: Các tiêu lý đất khu vực Phù Sa 125 Bảng 3.17: Lƣợng nƣớc tƣới thiếu ƣớc tính hệ số suất tƣơng đối 129 Bảng 3.18: Chi phí xây dựng hệ thống SCADA Yên Lập theo dự toán VWRAP (Bản gốc tiếng Anh) 133 Bảng 3.19: Chi phí xây dựng hệ thống SCADA Yên Lập theo dự toán VWRAP (Dịch từ gốc sang tiếng Việt) 134 Bảng 3.20: Chi phí trực tiếp xây dựng hệ thống SCADA Yên Lập theo công nghệ SCADA (IP Modem - RTU) đƣợc nghiên cứu đề xuất 134 ix DANH MỤC CÁC TỪ VIẾT TẮT ADB The Asian Development Bank (Ngân hàng phát triển châu Á) CSDL Cơ sở liệu CSDLQH Cơ sở liệu quan hệ ĐHTL Đại học Thuỷ lợi EAO Exploitation Assistée Par Ordinateur (Quản lý khai thác với trợ giúp máy tính) FAO The Food and Agriculture Organization of the United Nations (Tổ chức Lƣơng thực Nông nghiệp Liên hợp quốc) FQM Flow and Quality Management (Quản lý dòng chảy chất lƣợng nƣớc) GS&ĐK Giám sát Điều khiển HĐH Hiện đại hóa HT Hệ thống HTC Hệ thống HTT Hệ thống tƣới HTTL Hệ thống thuỷ lợi HTTN Hệ thống thuỷ nông HTTT Hệ thống tƣới tiêu HTX Hợp tác xã ICID The International Commission on Irrigation and Drainage (Uỷ ban tƣới tiêu Quốc tế) IIASA The International Institute for Applied Systems Analysis (Viện phân tích hệ thống ứng dụng Quốc tế) INCA Irrigation Network Control and Analysis (Điều khiển phân tích hệ thống tƣới) KHDN Kế hoạch dùng nƣớc KH&KT Khoa học kỹ thuật KHTL Khoa học thuỷ lợi KTTV Khí tƣợng thuỷ văn MH Mơ hình MP Mơ x cao Khoai/ Lạc 102 vùng đất cao Rau màu vùng 103 đất cao Lúa vùng đất 104 69.98 06/02 1053.50 451.50 602.00 139.94 03/02 806.25 301.00 505.25 3193.54 19/02 4647.06 2147.06 2500.00 112.99 21/02 940.63 564.38 376.25 98.33 22/02 967.50 451.50 516.00 205.91 19/02 978.25 301.00 677.25 392.17 19/01 5775.00 1775.00 4000.00 24.66 23/01 967.50 451.50 516.00 429.33 14/01 1284.63 564.38 720.25 56.31 28/02 1053.50 451.50 602.00 127.18 21/02 817.00 301.00 516.00 3383.07 13/02 5147.06 2147.06 3000.00 143.00 17/02 1284.63 564.38 720.25 Ngơ vùng đất 105 Khoai/Lạc 106 vùng đất Rau màu vùng 107 đất Lúa vùng đất 108 thấp, trũng Khoai/Lạc 109 vùng đất thấp 110 111 Phƣơng án 12 Ngô vùng đất cao 112 113 114 Khoai/ Lạc vùng đất cao Rau màu vùng đất cao Lúa vùng đất 115 Ngơ vùng đất 187 Khoai/Lạc 116 vùng đất Rau màu vùng 117 đất Lúa vùng đất 118 thấp, trũng Khoai/Lạc 119 vùng đất thấp 97.34 21/02 967.50 451.50 516.00 106.51 19/02 978.25 301.00 677.25 595.96 21/01 5275.00 1775.00 3500.00 63.11 27/01 1053.50 451.50 602.00 496.56 02/02 940.63 564.38 376.25 82.64 07/02 1053.50 451.50 602.00 162.83 09/02 806.25 301.00 505.25 3078.69 12/02 5147.06 1647.06 3500.00 177.80 24/02 940.63 564.38 376.25 101.50 17/02 1311.50 451.50 860.00 155.21 22/02 817.00 301.00 516.00 609.58 22/01 5275.00 1775.00 3500.00 55.12 17/01 1311.50 451.50 860.00 120 121 Phƣơng án 13 Ngô vùng đất cao 122 123 124 Khoai/ Lạc vùng đất cao Rau màu vùng đất cao Lúa vùng đất 125 Ngơ vùng đất 126 127 128 129 Khoai/Lạc vùng đất Rau màu vùng đất Lúa vùng đất thấp, trũng Khoai/Lạc vùng đất thấp 130 188 131 Phƣơng án 14 Ngô vùng đất 269.22 04/02 940.63 564.38 376.25 37.24 30/01 1053.50 451.50 602.00 149.22 07/02 806.25 301.00 505.25 3229.34 22/02 4500.00 2000.00 2500.00 292.07 22/02 940.63 564.38 376.25 39.95 25/02 967.50 451.50 516.00 198.78 25/02 817.00 301.00 516.00 630.23 15/01 5850.00 1850.00 4000.00 21.77 21/01 1311.50 451.50 860.00 410.95 02/02 940.63 564.38 376.25 37.73 13/02 1053.50 451.50 602.00 133.30 04/02 806.25 301.00 505.25 3038.12 18/02 5147.06 2147.06 3000.00 228.65 16/02 1284.63 564.38 720.25 cao Khoai/ Lạc 132 vùng đất cao Rau màu vùng 133 đất cao Lúa vùng đất 134 Ngơ vùng đất 135 Khoai/Lạc 136 vùng đất Rau màu vùng 137 đất Lúa vùng đất 138 thấp, trũng Khoai/Lạc 139 vùng đất thấp 140 141 Phƣơng án 15 Ngô vùng đất cao 142 143 144 Khoai/ Lạc vùng đất cao Rau màu vùng đất cao Lúa vùng đất 145 Ngơ vùng đất 189 Khoai/Lạc 146 vùng đất Rau màu vùng 147 đất Lúa vùng đất 148 thấp, trũng Khoai/Lạc 149 vùng đất thấp 45.30 20/02 1311.50 451.50 860.00 105.81 21/02 817.00 301.00 516.00 515.52 21/01 5275.00 1775.00 3500.00 28.00 19/01 1311.50 451.50 860.00 204.38 30/01 940.63 564.38 376.25 33.00 13/02 1053.50 451.50 602.00 135.37 25/01 806.25 301.00 505.25 3186.50 20/02 4647.06 2147.06 2500.00 276.84 24/02 940.63 564.38 376.25 49.79 22/02 967.50 451.50 516.00 163.71 25/02 817.00 301.00 516.00 327.58 21/01 5275.00 1775.00 3500.00 42.26 20/01 1311.50 451.50 860.00 150 151 Phƣơng án 16 Ngô vùng đất cao 152 153 154 Khoai/ Lạc vùng đất cao Rau màu vùng đất cao Lúa vùng đất 155 Ngơ vùng đất 156 157 158 159 Khoai/Lạc vùng đất Rau màu vùng đất Lúa vùng đất thấp, trũng Khoai/Lạc vùng đất thấp 190 160 161 Phƣơng án 17 Ngô vùng đất 463.22 25/02 940.63 564.38 376.25 16.96 22/02 967.50 451.50 516.00 163.71 06/02 806.25 301.00 505.25 2977.03 28/02 4029.41 1529.41 2500.00 311.78 21/02 940.63 564.38 376.25 15.23 26/02 1311.50 451.50 860.00 201.62 23/02 817.00 301.00 516.00 471.53 18/01 5775.00 1775.00 4000.00 24.93 27/01 1053.50 451.50 602.00 137.68 28/01 940.63 564.38 376.25 24.02 18/02 1311.50 451.50 860.00 156.48 26/01 806.25 301.00 505.25 3431.87 21/02 4617.65 2117.65 2500.00 cao Khoai/ Lạc 162 vùng đất cao Rau màu vùng 163 đất cao Lúa vùng đất 164 Ngơ vùng đất 165 Khoai/Lạc 166 vùng đất Rau màu vùng 167 đất Lúa vùng đất 168 thấp, trũng Khoai/Lạc 169 vùng đất thấp 170 171 Phƣơng án 18 Ngô vùng đất cao 172 173 174 Khoai/ Lạc vùng đất cao Rau màu vùng đất cao Lúa vùng đất 191 Ngơ vùng đất 175 369.29 26/02 940.63 564.38 376.25 59.81 19/02 1311.50 451.50 860.00 183.46 18/02 817.00 301.00 516.00 536.66 19/01 5775.00 1775.00 4000.00 68.02 24/01 967.50 451.50 516.00 271.72 23/01 1284.63 564.38 720.25 82.45 30/01 1053.50 451.50 602.00 147.36 02/02 806.25 301.00 505.25 2979.72 08/02 5205.88 1705.88 3500.00 320.29 18/02 1284.63 564.38 720.25 25.54 17/02 1311.50 451.50 860.00 110.72 25/02 817.00 301.00 516.00 526.58 21/01 5275.00 1775.00 3500.00 16.93 22/01 1311.50 451.50 860.00 Khoai/Lạc 176 vùng đất Rau màu vùng 177 đất Lúa vùng đất 178 thấp, trũng Khoai/Lạc 179 vùng đất thấp 180 181 Phƣơng án 19 Ngô vùng đất cao 182 183 184 Khoai/ Lạc vùng đất cao Rau màu vùng đất cao Lúa vùng đất 185 Ngơ vùng đất 186 187 188 189 Khoai/Lạc vùng đất Rau màu vùng đất Lúa vùng đất thấp, trũng Khoai/Lạc 192 vùng đất thấp 190 191 Phƣơng án 20 Ngô vùng đất 368.01 25/01 940.63 564.38 376.25 76.83 23/01 967.50 451.50 516.00 126.04 28/02 806.25 301.00 505.25 3077.31 27/02 4500.00 2000.00 2500.00 276.29 28/02 940.63 564.38 376.25 14.46 27/02 1311.50 451.50 860.00 113.89 20/02 978.25 301.00 677.25 519.63 14/01 5850.00 1850.00 4000.00 51.69 21/01 1311.50 451.50 860.00 cao 192 193 194 Khoai/ Lạc vùng đất cao Rau màu vùng đất cao Lúa vùng đất 195 Ngơ vùng đất 196 197 198 199 Khoai/Lạc vùng đất Rau màu vùng đất Lúa vùng đất thấp, trũng Khoai/Lạc vùng đất thấp 193 Phụ lục 5: Những sơ đồ thuật toán chƣơng trình phần mềm * Chương trình phần mềm (PL5): Thiết lập giao diện toán QLKT “tối ƣu“ HTT protected void Creating_Click(object sender, EventArgs e) { GridView1.Visible = false; UpdPnlInput.Visible = true; CmdInitializeSolutions.Visible = true; CmdObjFunctions.Visible = true; LblCropCodes.Visible = true; TxtS0C1.Visible = true; TxtS0C2.Visible = true; TxtS0C3.Visible = true; TxtS0C4.Visible = true; LblY5.Visible = true; if (TxtNoCrops.Text == "1") { TxtCrop1.Visible = true; } if (TxtNoCrops.Text == "2") { TxtCrop1.Visible = true; TxtCrop2.Visible = true; } if (TxtNoCrops.Text == "3") { TxtCrop1.Visible = true; TxtCrop2.Visible = true; TxtCrop3.Visible = true; } if (TxtNoCrops.Text == "4") { TxtCrop1.Visible = true; TxtCrop2.Visible = true; TxtCrop3.Visible = true; TxtCrop4.Visible = true; } LblCrops.Visible = true; LblSoils.Visible = true; LblCropCodes.Visible = true; if (TxtNoSoils.Text == "1") { TxtSoil1.Visible = true; TxtS1C1.Visible = true; TxtS1C2.Visible = true; TxtS1C3.Visible = true; TxtS1C4.Visible = true; TxtS1Y1.Visible = true; TxtS1Y2.Visible = true; TxtS1Y3.Visible = true; TxtS1Y4.Visible = true; } if (TxtNoSoils.Text == "2") { TxtSoil1.Visible = true; TxtS1C1.Visible = true; TxtS1C2.Visible = true; TxtS1C3.Visible = true; TxtS1C4.Visible = true; TxtSoil2.Visible = true; 194 TxtS2C1.Visible = true; TxtS2C2.Visible = true; TxtS2C3.Visible = true; TxtS2C4.Visible = true; TxtS1Y1.Visible = true; TxtS1Y2.Visible = true; TxtS1Y3.Visible = true; TxtS1Y4.Visible = true; LblY2.Visible = true; TxtS2Y1.Visible = true; TxtS2Y2.Visible = true; TxtS2Y3.Visible = true; TxtS2Y4.Visible = true; LblLowAmax.Visible = false; TxtAmax31.Visible = false; TxtAmax32.Visible = false; TxtAmax33.Visible = false; TxtAmax34.Visible = false; LblLowAmin.Visible = false; TxtAmin31.Visible = false; TxtAmin32.Visible = false; TxtAmin33.Visible = false; TxtAmin34.Visible = false; LblLowTpmax.Visible = false; TxtTpmax31.Visible = false; TxtTpmax32.Visible = false; TxtTpmax33.Visible = false; TxtTpmax34.Visible = false; LblLowTpmin.Visible = false; TxtTpmin31.Visible = false; TxtTpmin32.Visible = false; TxtTpmin33.Visible = false; TxtTpmin34.Visible = false; } if (TxtNoSoils.Text == "3") { TxtSoil1.Visible = true; LblCropCodes1.Visible = true; LblCropCodes2.Visible = true; LblCropCodes3.Visible = true; TxtS1C1.Visible = true; TxtS1C2.Visible = true; TxtS1C3.Visible = true; TxtS1C4.Visible = true; TxtSoil2.Visible = true; TxtS2C1.Visible = true; TxtS2C2.Visible = true; TxtS2C3.Visible = true; TxtS2C4.Visible = true; TxtSoil3.Visible = true; TxtS3C1.Visible = true; TxtS3C2.Visible = true; TxtS3C3.Visible = true; TxtS3C4.Visible = true; TxtS1Y1.Visible = true; TxtS1Y2.Visible = true; TxtS1Y3.Visible = true; TxtS1Y4.Visible = true; LblY2.Visible = true; TxtS2Y1.Visible = true; TxtS2Y2.Visible = true; TxtS2Y3.Visible = true; 195 TxtS2Y4.Visible LblY3.Visible = TxtS3Y1.Visible TxtS3Y2.Visible TxtS3Y3.Visible TxtS3Y4.Visible = true; true; = true; = true; = true; = true; } if (TxtNoSoils.Text == "4") { TxtSoil1.Visible = true; TxtS1C1.Visible = true; TxtS1C2.Visible = true; TxtS1C3.Visible = true; TxtS1C4.Visible = true; TxtSoil2.Visible = true; TxtS2C1.Visible = true; TxtS2C2.Visible = true; TxtS2C3.Visible = true; TxtS2C4.Visible = true; TxtSoil3.Visible = true; TxtS3C1.Visible = true; TxtS3C2.Visible = true; TxtS3C3.Visible = true; TxtS3C4.Visible = true; TxtS1Y1.Visible = true; TxtS1Y2.Visible = true; TxtS1Y3.Visible = true; TxtS1Y4.Visible = true; LblY2.Visible = true; TxtS2Y1.Visible = true; TxtS2Y2.Visible = true; TxtS2Y3.Visible = true; TxtS2Y4.Visible = true; LblY3.Visible = true; TxtS3Y1.Visible = true; TxtS3Y2.Visible = true; TxtS3Y3.Visible = true; TxtS3Y4.Visible = true; } PnlLatitude.Visible = true; } * Chương trình phần mềm (PL5): Tính tốn biến định To_c,v, Sc,v, Mc,v chế độ tƣới protected void CmdSolutions_Click(object sender, EventArgs e) { Random Rand = new Random(); int Tg = 1, Tsk = 1, t = 1, Tini = 1, Tdev = 1, Tmid = 1, Tlate = 1, TsOrder = 1, TsLastOrder = 1; float IniHmin = 0, IniHmax = 0, DevHmin = 0, DevHmax = 0, MidHmin = 0, MidHmax = 0, LateHmin = 0, LateHmax = 0; float IniHmin0 = 0, IniHmax0 = 0, DevHmin0 = 0, DevHmax0 = 0, MidHmin0 = 0, MidHmax0 = 0, LateHmin0 = 0, LateHmax0 = 0; float Hini = 0, Hdev = 0, Hmid = 0, Hlate = 0; float Kc_ini = 0, Kc_dev = 0, Kc_mid = 0, Kc_late = 0; Tg = Convert.ToInt16(TxtTg.Text); Tsk = Convert.ToInt16(TxtSkTime.Text); Tini = Convert.ToInt16(TxtIniTime.Text); Tdev = Convert.ToInt16(TxtDevTime.Text); Tmid = Convert.ToInt16(TxtMidTime.Text); 196 Tlate = Convert.ToInt16(TxtLateTime.Text); Kc_ini = Convert.ToSingle(TxtKc_ini.Text); Kc_mid = Convert.ToSingle(TxtKc_mid.Text); Kc_late = Convert.ToSingle(TxtKc_end.Text); Tstrt = Convert.ToSingle(Math.Pow((A * Hstrt * (1 - MstAtTo) / (K1 / (1 Anfa))), (1 / (1 - Anfa)))); Tstrt = Math.Round(Tstrt, 0); KstrtMean = Convert.ToSingle(A * Hstrt * (1 - MstAtTo) / Tstrt); for (s = 1; s < Ns; s++) { Urand[n, s] = Convert.ToSingle(Rand.NextDouble()); D[n, s, g, obj] = Amin[n] + (Amax[n] - Amin[n]) * Convert.ToSingle(Urand[n, s]); } foreach (OnM OnM in OnMs) k = Order1stDate[t]; if (OnM.ProductID > k - 1) { if (OnM.ProductID == k) { h[k - 1] = 10 * Hini * A * MstAtTo; h[i] = Convert.ToSingle((h[k - 1] + 10 * P[i]) + dH[i] - 10 * (EToIni[i] + EToDev[i] + EToMid[i] + EToLate[i])); if (h[i] hmax[i]) h[i] = Convert.ToSingle(hmax[i]); if (OnM.ProductID > (k - 1) + && OnM.ProductID 0) { Txt10.Text = Convert.ToString(QTBcc[t]); PpStAlarm.Visible = true; TxtQTBcc.Text = Convert.ToString(QTBcc[t]); LblPpStONnOFF.Text = "Mưa to - Hãy dừng máy bơm"; } if (QCcc[t] > 0) { PnlSluice.Visible = true; LblSluiceOnC.Text = "Mưa to - Hãy đóng cống"; } } else { if (Convert.ToSingle(TxtFieldWL.Text)