-1- Đại Học Quốc Gia Tp Hồ Chí Minh TRƯỜNG ĐẠI HỌC BÁCH KHOA NGUYỄN VŨ HUY TỐI ƯU QUẢN LÍ NGUỒN NƯỚC TRÊN LƯU VỰC SÔNG ĐỒNG NAI Chuyên ngành : XÂY DỰNG CÔNG TRÌNH THỦY Mã số ngành : 60 58 40 LUẬN VĂN THẠC SĨ TP HỒ CHÍ MINH, tháng 10 năm 2004 -2- CÔNG TRÌNH ĐƯC HOÀN THÀNH TẠI TRƯỜNG ĐẠI HỌC BÁCH KHOA ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH Cán hướng dẫn khoa học: TS LÊ ĐÌNH HỒNG Cán chấm nhận xét : TS TRƯƠNG CHÍ HIỀN Cán chấm nhận xét : TS NGUYỄN THỐNG Luận văn thạc só bảo vệ HỘI ĐỒNG CHẤM BẢO VỆ LUẬN VĂN THẠC SĨ TRƯỜNG ĐẠI HỌC BÁCH KHOA, ngày 29 tháng 12 năm 2004 -3- TRƯỜNG ĐẠI HỌC BÁCH KHOA CỘNG HOÀ XÃ HỘI CHỦ NGHĨA VIỆT NAM PHÒNG ĐÀO TẠO SĐH ĐỘC LẬP – TỰ DO – HẠNH PHÚC Tp Hồ Chí Minh, ngày … tháng … năm 2003 NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên : NGUYỄN VŨ HUY Phái: Nam Ngày, tháng, năm sinh : 25 –094 -1977 Nơi sinh : NAM ĐỊNH Chuyên ngành : Xây dựng công trình thủy MSHV : CTTH 13.002 I TÊN ĐỀ TÀI : Tối ưu quản lí nguồn nước lưu vực sông Đồng Nai II NHIỆM VỤ VÀ NỘI DUNG : Trước vấn đề quản lí nguồn nước đặt cho LVSĐN, luận văn lần tập trung nghiên cứu nhằm tìm qui trình vận hành hệ thống phân phối sử dụng nguồn nước lưu vực sông cách tối ưu mặt kinh tế, đáp ứng nhu cầu sử dụng nước cho phát triển kinh tế đáp ứng yêu cầu bảo vệ môi trường, bảo vệ nguồn nước III NGÀY GIAO NHIỆM VỤ: 12 - 2003 IV NGÀY HOÀN THÀNH NHIỆM VỤ: 31 - 10 - 2004 V HỌ VÀ TÊN CÁN BỘ HƯỚNG DẪN: Tiến só Lê Đình Hồng - Trường Đại học Bách Khoa thành phố Hồ Chí Minh CÁN BỘ HƯỚNG DẪN TS LÊ ĐÌNH HỒNG CHỦ NHIỆM NGÀNH TS NGUYỄN VĂN ĐĂNG BM QUẢN LÝ CHUYÊN NGÀNH TS HUỲNH THANH SƠN Nội dung đề cương luận văn thạc só Hội Đồng Chuyên Ngành thông qua Ngày … tháng … năm 2004 PHÒNG ĐÀO TẠO SAU ĐẠI HỌC KHOA QUẢN LÝ NGÀNH -4- Lời cam đoan Tôi xin cam đoan tất kết nghiên cứu luận văn cá nhân thực khoá học cao học chịu hoàn toàn trách nhiệm lời cam đoan -5- Lời cảm ơn Luận văn thực khoá 13 cao học ngành Xây dựng công trình thủy trường Đại học Bách khoa thành phố Hồ Chí Minh từ năm 2002 đến 2004 Trong suốt trình học tập nghiên cứu luận văn này, nhận giúp đỡ chân thành nhiệt tình q thầy cô, bạn bè đồng nghiệp Tôi xin bày tỏ lòng cảm ơn sâu sắc đến thầy TS Lê Đình Hồng, môn Kó thuật tài nguyên nước, người trực tiếp hướng dẫn tận tình thực luận văn Tôi xin chân thành cảm ơn TS Claudia Ringler, Viện nghiên cứu Chính sách Lương thực Quốc tế (IFPRI), giúp đỡ nhiệt tình việc phát triển mô hình tối ưu Tôi xin chân thành gửi lời cảm ơn đến thầy cô môn Kó thuật Tài nguyên nước tận tình giảng dạy giúp đỡ thời gian học tập trường Tôi xin chân thành cảm ơn bạn bè lớp cao học XDCTT K13 giúp đỡ chia xẻ kinh nghiệm học tập Tôi xin chân thành cảm ơn Phân viện Khảo sát Qui hoạch Thủy lợi Nam Bộ giúp đỡ tạo điều kiện cho tham gia khoá học -6- TÓM TẮT Nguồn tài nguyên nước lưu vực sông Đồng Nai đóng vai trò quan trọng trình tồn phát triển kinh tế - xã hội tỉnh miền Đông Nam Bộ nằm lưu vực, đặc biệt thành phố Hồ Chí Minh trung tâm kinh tế lớn nước Mâu thuẫn yêu cầu phát triển nhanh kinh tế bảo vệ bền vững nguồn nước đặt công tác quản lí tổng hợp nguồn nước lưu vực sông Đồng Nai đứng trước khó khăn thách thức Giải quiết vấn đề quản lí nguồn nước lưu vực sông Đồng Nai lưu vực khác nước vấn đề cần quan tâm nghiên cứu nhà khoa học Luận văn thạc só “Quản lí tối ưu nguồn nước lưu vực sông Đồng Nai” thực khoá 13 cao học trường Đại học Bách khoa thành phố Hồ Chí Minh Luận văn thực phương pháp nghiên cứu lí thuyết ứng dụng mô hình tối ưu toán học cho lưu vực sông Luận văn tiếp cận xem xét vấn đề phân phối sử dụng nước bối cảnh lưu vực sông chịu tác động điều kiện tự nhiên, kinh tế xã hội Bằng phương pháp phân tích tổng hợp luận văn xác định mối quan hệ đối tượng sử dụng nước với qui luật vật lí qui luật kinh tế Kết từ mô hình tối ưu đưa qui trình vận hành hệ thống phân phối sử dụng nước lưu vực sông Đồng Nai quan điểm tối ưu kinh tế có xét đến bảo vệ môi trường Từ kết nghiên cứu luận văn đưa kiến nghị quản lí tổng hợp nguồn nước cho lưu vực sông Đồng Nai phương pháp tiếp cận quản lí nguồn nước cho lưu vực sông khác -7- ABSTRACT Water resources in Dong Nai river basin play an important role on the existing and development of socio-economic process of provinces in the South-East region of Viet Nam, particularly of Ho Chi Minh city - a biggest economic center of country The conflict between rapidly developing of economy and sustainable maintaining the water resources is putting the integrated water resources management in facing difficulties and challenges Solving the water resources management problems in Dong Nai river basin are the new study problem, which has drawn interest from various related scientists The thesis namely “Water Resources Management Optimization in Dong Nai River Basin” has been carried out during the master course at Ho Chi Minh City University of Technology The thesis has been theoretically studied and applied an optimization model for a river basin The thesis has approached and considered the water distribution and uses of a river basin under influence of natural factor, economic factor and social factor The model result has proposed an operational regulation for distribution and uses of water in Dong Nai river basin under economic viewpoint and environmental protection The thesis has proposed recommendations related to integrated water resources management for Dong Nai river basin as well as an approach of water management for other river basin -8- MỤC LỤC TÓM TẮT ABSTRACT MUÏC LUÏC HÌNH VẼ BAÛNG 10 MỞ ĐẦU 11 LÍ DO CỦA NGHIÊN CỨU 11 MUÏC ĐÍCH CỦA NGHIÊN CỨU 13 ĐỐI TƯNG VÀ PHẠM VI NGHIÊN CỨU 13 Ý NGHĨA KHOA HỌC VÀ THỰC TIỄN CỦA NGHIÊN CỨU 14 CHƯƠNG 1: TỔNG QUAN 16 1.1 CÁC NGHIÊN CỨU TRÊN THẾ GIỚI 16 1.2 CÁC NGHIÊN CỨU VỀ LƯU VỰC SÔNG ĐỒNG NAI 19 CHƯƠNG 2: NỘI DUNG VÀ PHƯƠNG PHÁP NGHIÊN CỨU 21 2.1 NỘI DUNG NGHIÊN CỨU 21 2.2 PHƯƠNG PHÁP NGHIÊN CỨU 21 2.1.1 Phương pháp luận 21 2.1.2 Phương pháp cụ thể 22 CHƯƠNG 3: KẾT QUẢ VÀ THẢO LUẬN 24 3.1 ĐẶC ĐIỂM THỦY VĂN LƯU VỰC SÔNG ĐỒNG NAI 24 3.1.1 Hệ thống sông ngòi LVSĐN 24 3.1.2 Đặc điểm dòng chảy LVSĐN 28 3.2 HIỆN TRẠNG SỬ DỤNG NƯỚC TRÊN LƯU VỰC SÔNG ĐỒNG NAI 33 3.2.1 Các công trình thủy điện dòng 33 3.2.2 Nước sử dụng cho nông nghiệp 37 3.2.3 Nước sử dụng cho sinh hoạt SXCN 38 3.2.4 Nước sử dụng cho bảo vệ môi trường đẩy mặn 40 3.3 MÔ HÌNH TỐI ƯU PHÂN PHỐI NGUỒN NƯỚC 40 3.3.1 Mô hình tối ưu phân phối nguồn nước LVSĐN 40 3.3.2 Kết mô hình tối ưu 56 CHƯƠNG 4: KẾT LUẬN VÀ KIẾN NGHỊ 75 4.1 KẾT LUẬN 75 4.2 KIẾN NGHỊ 75 DANH MỤC TÀI LIỆU THAM KHẢO 76 PHUÏ LUÏC 77 -9- HÌNH VẼ Hình 3.1: Bản đồ hành LVSĐN vùng phụ cận 25 Hình 3.2: Sơ đồ mô hình tổng hợp kinh tế – thủy văn LVSĐN 43 Hình 3.3: Đường quan hệ lượng nước tưới suất lúa Hè Thu 54 Hình 3.4: Sơ đồ LVSĐN mô hình tối öu 57 Hình 3.5: Đường mực nước hồ Trị An 71 Hình 3.6: Điện lượng NMTĐ Trị An 71 Hình 3.7: Đường mực nước hồ Thác Mơ 72 Hình 3.8: Điện lượng NMTĐ Thác Mơ 72 Hình 3.9: Đường mực nước hồ Đa Nhim 73 Hình 3.10: Điện lượng NMTĐ Đa Nhim 73 Hình 3.11: Đường mực nước hồ Dầu Tiếng 74 Hình 3.12: Dòng chảy trường hợp cửa sông Đồng Nai 74 - 10 - BẢNG Bảng 3.1: Các tiêu trung bình năm trường hợp 69 Bảng 3.2: Mực nước hồ chưa trường hợp 69 Bảng 3.3: Diện tích tưới loại trồng trường hợp 70 - 117 - be1.26 be1.27 be1.28 be1.29 be1.30 be1.31 95.51 30.87 65.54 32.08 64.23 43.18 65.70 33.05 178.70 94.69 135.39 80.51 20.33 19.85 24.69 21.78 76.71 81.37 14.41 14.65 50.31 205.35 474.85 555.39 597.15 521.61 219.58 26.78 153.26 306.66 417.03 553.60 553.94 582.69 509.98 212.56 19.03 48.24 128.62 401.84 779.47 712.02 687.17 604.27 248.29 15.16 14.44 23.02 125.00 336.03 416.12 447.37 380.35 223.06 191.83 297.36 535.51 686.61 921.45 593.33 622.54 530.63 249.98 117.70 65.19 189.01 362.12 417.94 556.27 991.68 520.47 194.64 be2.1 be2.2 be2.3 be2.4 be2.5 be2.6 be2.7 be2.8 be2.9 be2.10 be2.11 be2.12 be2.13 be2.14 be2.15 be2.16 be2.17 be2.18 be2.19 be2.20 be2.21 be2.22 be2.23 be2.24 be2.25 be2.26 be2.27 be2.28 be2.29 be2.30 be2.31 85.55 44.01 24.00 26.41 65.25 160.76 269.39 377.71 615.37 335.52 282.42 157.54 66.45 34.76 22.36 25.95 65.00 169.13 277.96 535.92 485.20 394.64 249.14 145.17 79.44 42.72 21.45 23.33 57.51 135.51 239.74 467.41 411.97 369.99 229.55 133.25 52.50 28.55 16.55 17.70 38.33 80.27 335.39 397.50 567.49 580.84 221.51 149.80 72.96 41.47 31.28 33.13 60.56 151.14 265.03 690.71 433.51 445.55 409.25 152.64 66.96 39.00 30.67 31.78 33.96 279.78 508.09 487.87 550.62 510.93 243.70 110.22 51.00 29.64 18.21 17.03 48.61 197.07 273.20 533.59 760.86 505.57 268.38 116.83 58.04 24.19 8.38 1.14 13.77 22.11 128.67 264.09 496.13 324.49 201.66 96.18 47.78 20.13 9.99 13.17 36.59 98.99 188.83 434.22 934.70 598.84 251.37 137.67 62.81 28.28 15.75 12.29 63.93 142.61 336.38 470.09 330.69 452.97 175.63 87.72 55.34 26.54 13.20 8.42 41.65 159.51 212.61 346.85 675.55 526.04 328.69 182.27 102.93 57.50 28.28 21.80 58.90 212.80 268.59 638.80 448.65 438.27 295.82 147.95 79.41 33.99 16.74 23.59 42.21 159.20 382.37 536.54 930.76 567.20 396.42 190.76 80.94 42.09 27.27 12.80 12.53 53.55 158.80 270.52 296.89 408.54 263.84 95.27 68.51 24.68 16.53 24.31 36.77 74.55 132.90 366.46 365.37 476.62 173.48 89.08 67.04 33.22 29.30 28.87 100.63 161.51 260.80 389.87 445.62 518.86 323.53 179.93 83.22 40.09 22.39 18.92 50.27 35.77 172.62 969.96 793.07 670.83 379.88 211.73 103.92 44.05 24.91 16.10 18.96 104.43 289.56 296.15 464.28 493.52 328.98 176.35 95.97 46.96 28.02 34.45 102.23 192.61 228.79 366.54 417.75 674.88 300.65 147.18 69.10 29.27 33.91 30.92 64.17 129.73 391.74 609.18 771.12 634.97 305.47 152.53 74.83 40.86 34.23 20.68 40.79 298.39 254.82 651.60 738.59 561.10 436.23 206.56 107.64 47.95 30.96 27.16 35.81 130.61 513.80 755.50 890.33 1008.20 361.17 147.53 89.19 44.90 29.97 29.52 38.81 257.64 320.98 1066.27 794.29 725.85 310.68 144.02 78.50 37.64 48.43 33.33 123.23 177.19 421.90 672.57 780.17 1017.68 444.42 231.36 68.73 31.09 21.35 1.48 10.77 157.72 302.34 539.67 419.39 310.45 327.42 93.99 36.16 28.52 17.49 13.12 12.11 18.66 219.76 539.75 471.90 388.29 379.60 145.62 69.40 38.32 21.67 19.88 105.42 228.69 332.26 560.80 432.73 435.59 507.23 188.08 83.94 40.18 26.09 18.20 38.03 72.86 367.72 741.49 544.29 495.02 492.22 155.35 76.98 41.20 21.19 16.20 15.03 7.72 95.97 305.07 219.59 230.13 298.21 188.51 175.22 92.82 75.21 188.10 291.54 525.04 673.19 903.42 581.75 610.38 520.24 245.07 132.77 78.94 79.79 115.40 63.91 185.30 355.02 409.74 545.38 972.29 510.29 190.84 sg1.1 sg1.2 sg1.3 sg1.4 sg1.5 sg1.6 sg1.7 sg1.8 sg1.9 sg1.10 sg1.11 sg1.12 sg1.13 sg1.14 sg1.15 sg1.16 70.31 53.54 52.01 48.08 52.15 53.60 112.95 191.91 288.52 373.93 282.97 107.00 81.74 57.50 53.86 49.71 74.49 148.81 185.35 408.24 465.70 335.98 255.96 116.19 85.68 59.95 54.85 52.07 57.18 82.56 139.92 134.91 365.78 410.06 211.64 102.29 72.61 54.24 51.69 48.83 52.07 62.42 135.07 241.38 333.95 279.57 156.92 90.74 75.13 56.25 53.46 52.51 70.82 131.03 213.92 220.03 230.48 380.04 279.88 107.19 62.46 45.72 45.45 41.81 56.35 147.56 204.60 165.15 351.79 561.45 310.86 105.61 76.87 63.38 56.78 47.95 56.51 90.98 100.98 197.40 430.27 348.19 215.91 98.46 76.87 63.38 56.78 47.95 45.26 45.88 146.51 249.63 349.92 273.20 187.14 86.51 93.42 56.39 48.02 49.30 187.62 188.65 179.51 261.95 419.10 350.36 174.16 126.23 99.98 71.80 45.59 60.55 161.13 155.49 197.96 192.28 306.06 315.06 145.80 103.92 80.89 52.28 39.59 21.51 37.34 145.70 194.99 267.36 310.00 313.99 249.12 131.00 104.32 75.43 45.24 33.02 42.18 97.64 184.38 233.13 251.35 265.43 172.96 104.30 92.08 71.87 47.97 62.57 123.74 106.53 162.93 369.35 452.49 313.56 229.13 139.20 114.80 86.03 55.84 40.59 19.10 17.37 68.35 217.97 363.89 340.32 232.58 117.42 105.15 77.51 50.76 39.42 54.85 87.22 143.91 192.47 285.20 349.93 204.09 122.80 100.01 72.70 55.42 37.01 54.26 77.68 79.17 123.80 212.00 305.55 256.66 126.58 - 118 - sg1.17 sg1.18 sg1.19 sg1.20 sg1.21 sg1.22 sg1.23 sg1.24 sg1.25 sg1.26 sg1.27 sg1.28 sg1.29 sg1.30 sg1.31 107.46 84.28 54.59 42.74 33.64 28.20 59.54 326.66 511.38 302.93 132.24 122.83 104.83 79.81 51.75 41.45 33.53 31.83 43.44 111.31 315.68 402.56 267.08 127.71 110.51 83.51 59.84 78.15 154.84 147.67 167.94 194.13 317.23 320.68 185.20 123.50 101.03 73.69 57.61 89.55 106.41 145.54 286.75 373.26 473.30 458.86 273.84 151.33 122.54 84.48 57.77 44.01 40.66 107.93 204.50 361.24 474.83 404.25 250.93 144.74 119.19 86.75 59.49 46.22 41.92 71.98 230.88 310.91 381.10 412.21 218.09 127.41 99.48 65.39 48.85 46.50 74.59 88.96 161.08 313.51 428.07 353.41 219.05 130.12 107.43 76.71 80.49 61.90 46.68 94.45 192.15 210.39 310.31 365.71 230.82 121.89 99.85 72.14 59.81 78.41 78.40 184.27 203.99 173.88 388.00 376.32 153.71 108.74 83.33 57.89 44.38 26.39 28.93 57.85 127.09 209.77 335.59 329.82 195.05 122.08 103.17 72.04 49.10 52.41 123.53 123.38 146.40 242.42 372.60 491.86 436.93 196.54 140.35 104.34 71.46 93.91 155.78 81.23 154.62 249.49 242.35 318.06 256.04 124.52 101.30 73.54 44.60 38.98 46.58 56.89 87.10 280.51 455.88 406.02 347.66 253.43 124.41 69.19 54.56 104.04 178.54 276.72 311.85 434.49 444.79 546.90 467.88 195.47 91.90 59.95 58.39 84.50 71.30 212.00 393.00 272.82 235.64 814.82 289.45 149.40 sg2.1 sg2.2 sg2.3 sg2.4 sg2.5 sg2.6 sg2.7 sg2.8 sg2.9 sg2.10 sg2.11 sg2.12 sg2.13 sg2.14 sg2.15 sg2.16 sg2.17 sg2.18 sg2.19 sg2.20 sg2.21 sg2.22 sg2.23 sg2.24 sg2.25 sg2.26 sg2.27 sg2.28 sg2.29 sg2.30 sg2.31 30.48 23.20 35.44 24.92 37.12 25.98 31.47 23.51 32.57 24.36 27.08 19.81 33.32 27.48 33.32 27.48 40.47 24.43 43.31 31.11 35.03 22.64 45.18 32.68 39.88 31.14 49.74 37.28 45.56 33.58 43.34 31.50 46.55 36.51 45.43 34.57 47.89 36.19 43.77 31.93 53.09 36.60 51.64 37.59 43.10 28.33 46.55 33.24 43.26 31.26 36.10 25.09 44.70 31.21 60.80 45.21 43.87 31.86 53.92 29.97 39.83 25.98 sg3.1 sg3.2 sg3.3 sg3.4 sg3.5 sg3.6 sg3.7 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 22.55 23.36 23.78 22.42 23.17 19.69 24.61 24.61 20.81 19.74 17.14 19.61 20.78 24.19 21.99 24.00 23.65 22.42 25.93 24.96 25.02 25.77 21.16 34.87 25.90 19.23 21.27 30.96 19.31 23.65 25.31 0.00 0.00 0.00 0.00 0.00 0.00 0.00 20.84 21.54 22.58 21.18 22.76 18.12 20.79 20.79 21.36 26.23 9.33 14.31 27.11 17.57 17.08 16.04 18.51 17.96 33.85 38.80 19.08 20.01 20.14 26.83 33.98 11.43 22.71 40.69 16.90 45.07 36.62 0.00 0.00 0.00 0.00 0.00 0.00 0.00 22.61 32.27 24.78 22.58 30.69 24.43 24.48 19.61 81.29 69.80 16.18 18.27 53.62 8.28 23.76 23.52 14.57 14.52 67.09 46.10 17.62 18.16 32.33 20.22 33.96 12.53 53.51 67.50 20.17 77.38 30.91 0.00 0.00 0.00 0.00 0.00 0.00 0.00 23.22 48.96 83.16 125.04 162.07 122.63 46.36 64.49 80.33 176.94 201.84 145.62 110.94 50.35 35.77 60.64 58.47 158.53 177.71 91.73 44.33 27.03 58.55 104.62 144.74 121.17 68.01 39.32 56.79 92.73 95.35 99.90 164.72 121.31 46.44 63.94 88.68 71.57 152.46 243.33 134.73 45.77 39.42 43.77 85.55 186.47 150.90 93.57 42.67 19.88 63.50 108.18 151.66 118.41 81.10 37.50 81.73 77.78 113.48 181.57 151.78 75.45 54.69 67.37 85.76 83.30 132.61 136.49 63.17 45.02 63.12 84.48 115.84 134.29 136.04 107.93 56.76 42.30 79.87 101.00 108.89 114.98 74.93 45.18 46.16 70.58 160.01 196.03 135.85 99.27 60.32 7.52 29.62 94.44 157.65 147.45 100.75 50.86 37.79 62.35 83.38 123.56 151.60 88.41 53.19 33.64 34.31 53.62 91.83 132.37 111.20 54.83 12.21 25.79 141.53 221.54 131.24 57.28 53.22 13.79 18.83 48.24 136.75 174.39 115.71 55.34 63.97 72.75 84.10 137.43 138.93 80.22 53.51 63.06 124.22 161.72 205.05 198.79 118.64 65.57 46.76 88.60 156.50 205.70 175.14 108.71 62.70 31.18 100.01 134.70 165.11 178.57 94.48 55.20 38.54 69.77 135.82 185.46 153.10 94.89 56.38 40.93 83.24 91.15 134.45 158.43 100.00 52.82 79.83 88.36 75.34 168.09 163.03 66.59 47.11 25.06 55.07 90.88 145.39 142.89 84.50 52.90 53.45 63.42 105.02 161.43 213.09 189.29 85.15 35.20 66.99 108.10 105.00 137.80 110.91 53.94 24.65 37.74 121.52 197.51 175.89 150.62 109.79 119.91 135.13 188.26 192.74 236.98 202.75 84.69 91.86 170.29 118.22 102.10 353.09 125.43 64.74 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 - 119 - sg3.8 sg3.9 sg3.10 sg3.11 sg3.12 sg3.13 sg3.14 sg3.15 sg3.16 sg3.17 sg3.18 sg3.19 sg3.20 sg3.21 sg3.22 sg3.23 sg3.24 sg3.25 sg3.26 sg3.27 sg3.28 sg3.29 sg3.30 sg3.31 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 vd1.1 vd1.2 vd1.3 vd1.4 vd1.5 vd1.6 vd1.7 vd1.8 vd1.9 vd1.10 vd1.11 vd1.12 vd1.13 vd1.14 vd1.15 vd1.16 vd1.17 vd1.18 vd1.19 vd1.20 vd1.21 vd1.22 vd1.23 vd1.24 vd1.25 vd1.26 vd1.27 vd1.28 vd1.29 vd1.30 34.63 23.76 19.55 15.29 14.81 12.10 68.57 181.81 272.89 497.73 273.46 81.77 59.65 32.47 23.86 19.21 61.36 164.00 173.91 441.37 533.82 431.33 242.02 100.33 68.27 37.79 26.17 24.91 25.28 58.27 107.81 113.40 386.67 560.96 190.41 72.29 39.67 25.30 18.83 17.11 14.65 26.15 100.76 241.16 339.81 332.58 126.72 49.01 51.32 27.92 19.98 21.15 44.70 103.34 183.98 148.04 134.37 342.06 262.88 95.30 70.42 36.36 23.01 20.87 35.94 148.83 224.21 140.86 306.01 724.37 413.09 97.73 59.94 30.39 20.97 24.91 32.33 70.11 74.73 165.28 534.11 549.05 250.36 82.79 52.18 24.02 17.54 12.23 11.65 26.15 175.35 240.01 265.65 321.94 197.51 67.12 30.51 18.22 15.08 13.50 38.60 117.13 200.37 250.99 248.39 317.85 199.14 48.32 29.19 19.64 14.38 13.58 27.80 80.95 177.02 176.29 168.92 284.34 102.64 33.24 22.31 14.20 10.18 23.59 58.26 117.34 197.80 361.93 343.54 440.33 209.98 50.92 30.88 22.06 17.33 13.66 29.46 51.14 188.08 261.49 287.61 356.25 225.04 53.17 34.31 23.93 19.87 23.72 24.00 64.02 114.64 146.40 343.36 435.05 233.28 50.35 33.35 22.23 17.57 13.14 54.24 76.62 57.88 155.83 289.89 547.12 251.76 56.78 35.03 23.42 19.79 25.48 96.96 142.17 144.96 178.09 357.20 605.91 238.72 56.09 36.94 25.74 21.91 42.28 128.83 96.60 145.17 157.01 330.38 549.71 256.53 67.07 40.02 28.79 21.56 18.01 19.69 43.99 111.85 326.93 539.14 623.16 189.03 58.42 36.32 23.15 19.95 15.06 11.46 24.60 37.23 86.46 291.57 512.06 276.83 55.36 32.70 20.61 16.69 15.50 39.24 101.22 167.75 213.79 385.30 424.69 143.93 51.72 33.00 23.47 38.19 59.69 25.44 110.45 327.41 352.18 488.23 702.62 278.93 70.31 46.20 30.36 24.56 17.39 19.31 58.61 159.58 460.55 505.44 466.12 283.36 64.79 42.16 28.28 23.38 19.67 23.81 71.93 250.83 291.33 244.66 533.32 159.69 52.31 38.81 25.35 19.79 16.28 12.40 37.58 138.26 338.92 285.69 442.79 161.04 47.73 35.06 25.33 27.59 17.65 17.06 90.72 182.27 195.87 313.14 573.55 209.77 63.80 32.01 22.57 16.37 22.84 63.53 160.24 156.39 73.01 403.78 499.71 112.16 29.62 30.67 20.64 16.45 12.03 25.39 49.84 91.71 163.68 530.82 674.66 152.95 49.23 39.67 24.19 19.34 18.20 45.43 113.87 181.14 177.02 424.13 796.61 459.46 102.61 57.29 39.53 31.52 29.83 83.62 97.43 188.43 217.67 253.89 482.86 273.40 57.51 37.82 26.01 20.95 17.18 16.04 26.33 91.28 190.09 388.10 498.45 271.85 137.46 148.04 84.12 65.54 107.65 196.92 276.80 280.27 403.88 546.63 712.88 610.60 244.32 - 120 - vd1.31 101.06 68.71 65.86 95.33 92.59 277.03 509.38 294.33 184.68 970.36 284.83 173.37 vd2.1 vd2.2 vd2.3 vd2.4 vd2.5 vd2.6 vd2.7 vd2.8 vd2.9 vd2.10 vd2.11 vd2.12 vd2.13 vd2.14 vd2.15 vd2.16 vd2.17 vd2.18 vd2.19 vd2.20 vd2.21 vd2.22 vd2.23 vd2.24 vd2.25 vd2.26 vd2.27 vd2.28 vd2.29 vd2.30 vd2.31 20.03 13.74 34.50 18.77 39.48 21.85 22.95 14.64 29.68 16.14 40.74 21.02 34.66 17.56 30.19 13.89 19.04 10.91 22.53 12.68 16.47 8.49 23.19 14.81 29.06 18.29 26.81 15.77 28.69 16.72 26.86 17.44 30.53 20.32 25.85 15.89 23.92 14.08 24.43 15.51 32.92 20.42 29.09 18.46 33.16 20.32 31.44 20.59 24.03 15.63 25.95 15.48 27.91 16.02 47.78 32.42 29.19 18.80 85.63 48.65 56.62 38.49 vd3.1 vd3.2 vd3.3 vd3.4 vd3.5 vd3.6 vd3.7 vd3.8 vd3.9 vd3.10 vd3.11 vd3.12 vd3.13 vd3.14 vd3.15 vd3.16 vd3.17 vd3.18 vd3.19 vd3.20 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.34 0.34 0.34 0.34 0.34 0.34 0.34 0.34 0.34 0.34 0.34 0.34 0.34 0.34 0.34 0.34 0.34 0.34 0.34 0.34 11.30 8.84 8.57 7.00 39.67 105.15 157.85 287.90 158.16 47.30 13.79 11.09 35.49 94.84 100.60 255.31 308.78 249.49 139.99 58.01 15.13 14.41 14.62 33.72 62.35 65.59 223.66 324.46 110.13 41.81 10.90 9.90 8.49 15.14 58.28 139.49 196.55 192.36 73.30 28.34 11.54 12.23 25.85 59.77 106.41 85.63 77.71 197.85 152.05 55.12 13.31 12.08 20.78 86.08 129.69 81.48 176.98 418.98 238.93 56.54 12.13 14.41 18.70 40.54 43.23 95.59 308.94 317.58 144.82 47.89 10.15 7.08 6.72 15.14 101.43 138.82 153.65 186.23 114.23 38.81 8.14 8.76 35.14 55.73 129.47 113.48 162.36 209.64 107.15 42.13 8.76 10.50 23.41 46.40 83.57 69.34 61.77 132.42 55.11 26.73 5.65 7.83 29.01 94.56 149.16 255.25 183.49 282.20 126.75 42.16 11.01 9.18 32.03 51.53 103.20 113.64 134.27 179.16 170.81 47.43 15.35 22.89 25.87 66.20 72.16 97.95 252.41 206.18 118.30 43.10 11.33 7.70 16.07 30.30 63.88 159.18 160.19 346.88 165.08 50.38 14.14 22.71 40.85 65.16 88.33 104.89 220.53 374.92 113.30 44.41 14.49 20.61 67.58 45.75 75.66 88.09 186.96 301.29 146.86 56.86 14.89 11.90 16.55 26.67 66.56 237.87 257.15 295.24 89.37 44.70 14.22 10.73 8.92 27.60 31.66 35.52 105.70 191.75 120.94 42.83 10.77 11.64 24.08 35.10 75.08 109.39 163.11 191.32 97.93 41.38 22.66 19.96 21.29 56.48 177.87 208.08 249.01 357.08 145.72 52.26 18.16 12.05 15.96 43.00 47.68 313.00 274.54 197.75 184.45 49.79 15.43 12.91 14.14 43.26 123.47 129.18 153.91 443.46 94.79 45.51 15.80 13.27 11.57 45.54 79.84 191.69 172.78 367.29 88.67 43.28 25.50 15.55 16.53 45.10 63.50 67.20 128.38 369.65 122.19 56.43 10.63 21.38 65.49 120.24 122.14 59.14 292.84 366.81 85.17 21.48 13.02 10.29 12.11 18.20 40.71 51.85 221.36 277.62 71.88 30.94 12.45 13.25 48.10 87.87 155.32 112.41 192.51 477.26 329.73 95.38 25.58 25.58 50.11 62.88 166.25 108.69 155.08 274.62 141.60 47.27 15.35 13.69 13.98 28.90 75.69 128.19 227.29 384.91 262.18 138.02 37.90 62.26 113.89 160.11 162.10 233.61 316.17 412.34 353.19 141.34 36.91 53.42 51.88 155.23 285.41 164.91 103.47 543.72 159.59 97.15 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 - 121 - vd3.21 vd3.22 vd3.23 vd3.24 vd3.25 vd3.26 vd3.27 vd3.28 vd3.29 vd3.30 vd3.31 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.34 0.34 0.34 0.34 0.34 0.34 0.34 0.34 0.34 0.34 0.34 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.36 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37 vd4.1 vd4.2 vd4.3 vd4.4 vd4.5 vd4.6 vd4.7 vd4.8 vd4.9 vd4.10 vd4.11 vd4.12 vd4.13 vd4.14 vd4.15 vd4.16 vd4.17 vd4.18 vd4.19 vd4.20 vd4.21 vd4.22 vd4.23 vd4.24 vd4.25 vd4.26 vd4.27 vd4.28 vd4.29 vd4.30 vd4.31 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ca1.1 ca1.2 ca1.3 ca1.4 ca1.5 ca1.6 ca1.7 ca1.8 ca1.9 ca1.10 ca1.11 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 - 122 - ca1.12 ca1.13 ca1.14 ca1.15 ca1.16 ca1.17 ca1.18 ca1.19 ca1.20 ca1.21 ca1.22 ca1.23 ca1.24 ca1.25 ca1.26 ca1.27 ca1.28 ca1.29 ca1.30 ca1.31 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ca2.1 0.00 0.00 0.00 0.00 0.00 ca2.2 0.00 0.00 0.00 0.00 0.00 ca2.3 0.00 0.00 0.00 0.00 0.00 ca2.4 0.00 0.00 0.00 0.00 0.00 ca2.5 0.00 0.00 0.00 0.00 0.00 ca2.6 0.00 0.00 0.00 0.00 0.00 ca2.7 0.00 0.00 0.00 0.00 0.00 ca2.8 0.00 0.00 0.00 0.00 0.00 ca2.9 0.00 0.00 0.00 0.00 0.00 ca2.10 0.00 0.00 0.00 0.00 0.00 ca2.11 0.00 0.00 0.00 0.00 0.00 ca2.12 0.00 0.00 0.00 0.00 0.00 ca2.13 0.00 0.00 0.00 0.00 0.00 ca2.14 0.00 0.00 0.00 0.00 0.00 ca2.15 0.00 0.00 0.00 0.00 0.00 ca2.16 0.00 0.00 0.00 0.00 0.00 ca2.17 0.00 0.00 0.00 0.00 0.00 ca2.18 0.00 0.00 0.00 0.00 0.00 ca2.19 0.00 0.00 0.00 0.00 0.00 ca2.20 0.00 0.00 0.00 0.00 0.00 ca2.21 0.00 0.00 0.00 0.00 0.00 ca2.22 0.00 0.00 0.00 0.00 0.00 ca2.23 0.00 0.00 0.00 0.00 0.00 ca2.24 0.00 0.00 0.00 0.00 0.00 ca2.25 0.00 0.00 0.00 0.00 0.00 ca2.26 0.00 0.00 0.00 0.00 0.00 ca2.27 0.00 0.00 0.00 0.00 0.00 ca2.28 0.00 0.00 0.00 0.00 0.00 ca2.29 0.00 0.00 0.00 0.00 0.00 ca2.30 0.00 0.00 0.00 0.00 0.00 ca2.31 0.00 0.00 0.00 0.00 0.00 ; EQUATIONS PHUONG TRINH 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 - 123 - OBJFUN PT Ham muc tieu N_BAL(rn,yr,pd) PT CBN nut song R_BAL(rev,yr,pd) PT CBN ho chua VOL_ELEV(rev,yr,pd) PT quan he Z-V AR_VOL(rev,yr,pd) PT quan he F-V dm_tfla(agdm,yr,pd) PT tong luong nuoc cap cho khu tuoi dm_tflm(mundm,yr,pd) PT tong luong nuoc cap cho khu SH dm_tfli(indm,yr,pd) PT tong luong nuoc cap cho khu SXCN hydr_rev(pwst,yr,pd) PT phat dien w_bala(agdm,yr,pd) PT CBN mat ruong cal_ylda(agdm,crop,yr) PT nang suat cay du nuoc cal_yldst(agdm,crop) PT nang suat cay thieu nuoc cal_muval(mundm) PT tinh loi nhuan tu cap nuoc SH cal_muvalu(mundm) PT tinh loi nhuan tu cap nuoc SH cho thi cal_muvalr(mundm) PT tinh loi nhuan tu cap nuoc SH cho nong thon cal_arval(agdm) PT tinh loi nhuan tu cap nuoc tuoi cal_powval(pwst) PT tinh loi nhuan tu phat dien cal_inval(indm) PT tinh loi nhuan twf cap nuoc SXCN def_dft(agdm,crop,yr,pd) PT xac dinh luong nuoc thieu tung thoi doan def_mdft1(agdm,crop,yr,pd) PT xac dinh luong nuoc thieu lon nhat area_lmt1(agdm,yr,pd) PT gioi han dien tich cay eta_lmt(agdm,crop,yr,pd) PT xac dinh ETm tung thoi doan cal_dp(agdm,crop,yr) PT tinh luong nuoc cay su dung cua ca mua vu cal_dp_stg(agdm,crop,yr,pd) PT tinh luong nuoc cay su dung tung thoi doan cal_eta_se(agdm,crop,yr) PT tinh ETa cua ca mua vu cal_et_st(agdm,crop,yr,pd) PT tinh ETa cua tung thoi doan cal_twacp(agdm,crop,yr) PT tinh tong luong nuoc cap cho cay cua ca mua vu cal_wacp(agdm,crop,yr,pd) PT tinh luong nuoc cap cho cay tung thoi doan cal_agwto(agdm,yr) PT tinh tong luong nuoc lay tu song cua nam cal_rf_m(mundm,yr,pd) PT tinh luong nuoc hoi qui tu khu su dung nuoc SH cal_rf_i(indm,yr,pd) PT tinh luong nuoc hoi qui tu khu su dung nuoc CN cal_rf_a(agdm,yr,pd) PT tinh luong nuoc hoi qui tu khu su dung nuoc tuoi cal_muwat(mundm,yr,pd) PT tinh luong nuoc cap cho khu su dung nuoc SH cal_muwatsf(mundm, yr,pd) PT tinh nhu cau nuoc SH tmuwat_lmt(mundm,yr,pd) PT tinh nhu cau nuoc SH gioi han yldm_lmt(agdm,crop,yr) PT tinh nhu cau nuoc tuoi gioi han cal_inwat(indm, yr,pd) PT tinh nhu cau nuoc SXCN miflow_lt(rn, rn_lo, yr,pd) PT xac dinh dong chay toi thieu luu cal_effy(agdm,crop,yr) PT tinh hieu suat tuoi eta_se_st(agdm,crop,yr) PT tinh luong nuoc boc hoi cua ca vu thieu nuoc ch_urb(mundm,yr,pd) PT tinh luong nuoc cap cho khu vuc thi ch_rur(mundm,yr,pd) PT tinh luong nuoc cap cho khu vuc nong thon ; VARIABLES BIEN QUIET DINH OBJ Bien muc tieu arval(agdm) Loi nhuan tu nong nghiep muval(mundm) Loi nhuan tu cap nuoc SH inval(indm) Loi nhuan tu cap nuoc cong nghiep powval(pwst) Loi nhuan tu phat dien muvalu(mundm) Loi nhuan tu cap nuoc SH cho khu vuc thi muvalr(mundm) Loi nhuan tu cap nuoc SH cho khu vuc nong thon ; - 124 - POSITIVE VARIABLES BIEN KHONG AM FLOW(rn,rn_lo,yr,pd) Dong chay giua cac nut song RES_ST(rev,yr,pd) Dung tich cac ho chua RES_AREA(rev,yr,pd) Dien tich mat ho cac ho chua RES_Da(rev,agdm,yr,pd) Luong nuoc lay tu ho cap cho khu tuoi RES_Dm (rev,mundm,yr,pd) Luong nuoc lay tu ho cap cho SH RES_Di (rev,indm,yr,pd) Luong nuoc lay tu ho cap cho KCN RES_N_T(rev,rn,yr,pd) Dong chay tu ho chua qua tua-bin RES_N_SP(rev,rn,yr,pd) Dong chay tu ho cha qua tran xa lu RES_RES(rev_up,rev,yr,pd) Dong chay tu ho TL qua ho HL h(rev,yr,pd) Muc nuoc cac ho chua N_RES(rn,rev,yr,pd) Dong chay tu nut song vao ho chua N_Da(rn,agdm,yr,pd) Luong nuoc lay tu song cap cho khu tuoi N_Dm(rn,mundm,yr,pd) Luong nuoc lay tu song cap cho khu SH N_Di(rn,indm,yr,pd) Luong nuoc lay tu song cap cho khu CN to_infa(agdm,yr,pd) Tong luong nuoc cap cho khu tuoi nam to_infi(indm,yr,pd) Tong luong nuoc cap cho khu SH nam to_infm(mundm,yr,pd) Tong luong nuoc cap cho khu CN nam power(pwst,yr,pd) Dien luong phat hang thang wacp(agdm,crop,yr,pd) Luong nuoc co san cho cay theo thang twacp(agdm,crop,yr) Luong nuoc co san cho cay theo nam tinwat(indm,yr,pd) Luong nuoc cap cho KCN tmuwat(mundm,yr,pd) Luong nuoc cap cho khu SH tmuwatr(mundm,yr,pd) Luong nuoc cap cho khu SH thi tmuwatu(mundm,yr,pd) Luong nuoc cap cho khu SH nong thon area(agdm,crop,yr) Dien tich cua tung loai cay ylda(agdm,crop,yr) Nang suat cay cap du nuoc yldst(agdm,crop) Nang suat cay khong du nuoc agwto(agdm,yr) Tong luong nuoc cap cho tuoi nam mdft(agdm,crop,yr) Luong nuoc thieu hut cho tuoi tung nam dft(agdm,crop,yr,pd) Luong nuoc thieu hut tung thang eta_st(agdm,crop,yr,pd) Boc thoat hoi thuc te theo thang eta_se(agdm,crop,yr) Boc thoat hoi thuc te theo vu mua dp(agdm,crop,yr) Luong nuoc cay su dung theo vu mua dp_stg(agdm,crop,yr,pd) Luong nuoc cay su dung theo thang retn_fla(agdm,yr,pd) Dong chay hoi qui tu khu tuoi retn_flm(mundm,yr,pd) Dong chay hoi qui tu khu SH retn_fli(indm,yr,pd) Dong chay hoi qui tu KCN inpts(agdm,crop,IN) Chi phi dau vao cho SXNN ; *CONSTRAINS CAC RANG BUOC KHAC inpts.fx(agdm,crop, "WATER")$dem_cp(agdm,crop)=maxwat(crop); inpts.fx(agdm,crop,IN)$(INS(IN)and dem_cp(agdm,crop))= inpts0(crop, IN); area.lo(agdm,crop,yr) = 0.5*areabnd(agdm,crop); area.up(agdm,crop,yr) = 1.2*areabnd(agdm,crop); tmuwat.lo(mundm,yr,pd) = 0.5*(mwd(mundm,pd)*pop(mundm)/1000000); tmuwat.up(mundm,yr,pd) = 1.5*(mwd(mundm,pd)*pop(mundm)/1000000); tmuwatr.lo(mundm,yr,pd)= 0.5*(mwdr(mundm,pd)*popr(mundm)/1000000); tmuwatr.up(mundm,yr,pd)= 1.5*(mwdr(mundm,pd)*popr(mundm)/1000000); tmuwatu.lo(mundm,yr,pd)= 0.5*(mwdu(mundm,pd)*popu(mundm)/1000000); tmuwatu.up(mundm,yr,pd)= 1.5*(mwdu(mundm,pd)*popu(mundm)/1000000); muvalr.fx('MTDHC')=0; tinwat.lo(indm,yr,pd) = 0.5*iwd(indm,pd); - 125 - tinwat.up(indm,yr,pd) = 1.5*iwd(indm,pd); flow.lo('dn4','dn5',yr,pd)=down_fldn(pd); flow.lo('sg2','sg3',yr,pd)=down_flsg(pd); RES_ST.up(rev,yr,pd)=R_up(rev); RES_ST.lo(rev,yr,pd)=stg0(rev); RES_N_T.up(rev,rn,yr,pd)=sum(pwst$RPLINK(rev,pwst),turb_max(pwst)); RES_N_T.fx('revDN','dn2',yr,pd)=0; RES_N_T.fx('revHT','ln3',yr,pd)=0; eta_st.fx(agdm,crop,yr,pd)$(NOT(dem_cp(agdm,crop) and cp_prd(agdm,crop,pd)))=0.0; dp_stg.fx(agdm,crop,yr,pd)$(NOT(dem_cp(agdm,crop) and cp_prd(agdm,crop,pd)))=0.0; wsfacp.fx(agdm,crop,yr,pd)$(NOT(dem_cp(agdm,crop) and cp_prd(agdm,crop,pd)))=0.0; h.up('revTA',yr,pd)= 62.0; h.up('revDT',yr,pd)= 24.4; h.up('revDN',yr,pd)= 1042.0; h.up('revHT',yr,pd)= 605.0; h.up('revDMi',yr,pd)= 325.0; h.up('revTM',yr,pd)= 218.0; h.lo('revTA',yr,pd)= 50.0; h.lo('revDT',yr,pd)= 17.0; h.lo('revDN',yr,pd)= 1018.0; h.lo('revHT',yr,pd)= 575.0; h.lo('revDMi',yr,pd)= 323.0; h.lo('revTM',yr,pd)= 198.0; *** *** objfun obj =E= SUM(agdm,arval(agdm)) + SUM(mundm,muval(mundm)) + SUM(pwst,powval(pwst)) + SUM(indm,inval(indm)) - pwgt1*sum(yr, (SUM(dem_cp(agdm,crop),yldst(agdm,crop)*area(agdm,crop,yr)*price(crop)* (mdft(agdm,crop,yr)- SUM(pd$cp_prd(agdm,crop,pd), dft(agdm,crop,yr,pd))/ (SUM(pd$cp_prd(agdm,crop,pd),1)+0.0000000001))))); N_BAL(rn,yr,pd)$((ord(rn)ne 7)and (ord(rn) ne 22)) SUM(rn_lo$RVLINK(rn,rn_lo), FLOW(rn,rn_lo,yr,pd))+ SUM(rev$NRLINK(rn,rev), N_RES(rn,rev,yr,pd)) + SUM(agdm$NDLINKa(rn,agdm), N_Da(rn,agdm,yr,pd)) + SUM(mundm$NDLINKm(rn,mundm), N_Dm(rn,mundm,yr,pd))+ SUM(indm$NDLINKi(rn,indm), N_Di(rn,indm,yr,pd)) =e= SUM(rn_up$RVLINK(rn_up,rn),FLOW(rn_up,rn,yr,pd)) + source(rn,yr,pd) + SUM(rev$RNLINK(rev,rn), RES_N_T(rev,rn,yr,pd)) + SUM(rev$RNLINK(rev,rn), RES_N_SP(rev,rn,yr,pd)) + SUM(dnlinka(agdm,rn), RETN_FLa(agdm,yr,pd)) + SUM(dnlinki(indm,rn), RETN_FLi(indm,yr,pd)) + SUM(dnlinkm(mundm,rn),RETN_FLm(mundm,yr,pd)); miflow_lt(rn,rn_lo,yr,pd)$(rvlink(rn,rn_lo) and (ord(rn) ne 1)) flow(rn,rn_lo,yr,pd) =g= 0.1*source(rn,yr,pd); - 126 - R_BAL(rev,yr,pd) RES_ST(rev,yr,pd-1)$(ord(pd) ne 1) +res_st0(rev)$((ord(pd) eq 1)and (ord(yr) eq 1)) +RES_ST(rev,yr-1,'12')$((ord(yr) ne 1) and (ord(pd)eq 1)) +SUM(rn$NRLINK(rn,rev), N_RES(rn,rev,yr,pd)) +SUM(rev_up$RRLINK(rev_up, rev), RES_RES(rev_up,rev,yr,pd)) =e= RES_ST(rev,yr,pd) +SUM(rn$RNLINK(rev,rn), RES_N_T(rev,rn,yr,pd)) +SUM(rn$RNLINK(rev,rn), RES_N_SP(rev,rn,yr,pd)) +SUM(indm$rdlinki(rev,indm), RES_Di(rev,indm,yr,pd)) +SUM(agdm$rdlinka(rev, agdm), RES_Da(rev,agdm,yr,pd)) +SUM(mundm$rdlinkm(rev,mundm), RES_Dm(rev,mundm,yr,pd)) +SUM(rev_lo$RRLINK(rev, rev_lo), RES_RES(rev, rev_lo,yr,pd)) +r_evap(rev,pd)*RES_AREA(rev,yr,pd)* 0.001; VOL_ELEV(rev,yr,pd) RES_ST(rev,yr,pd) =e= a(rev) * h(rev,yr,pd)**3 + b(rev)*h(rev,yr,pd)**2 + c(rev) *h(rev,yr,pd) + d(rev); AR_VOL(rev,yr,pd) RES_AREA(rev,yr,pd)=e= (e(rev)* (RES_ST(rev,yr,pd))**3 + f(rev)*(RES_ST(rev,yr,pd))**2 + g(rev)*(RES_ST(rev,yr,pd)) + i(rev)); Hydr_rev(pwst,yr,pd)$for_dam(pwst) SUM(rev$RPLINK(rev,pwst), SUM(rn$RNLINK(rev,rn),RES_N_T(rev,rn,yr,pd)* (h(rev,yr,pd)-h0(pwst))* pw_effy(pwst)*2.725) + SUM(rev_lo$RRLINK(rev,rev_lo),RES_RES(rev,rev_lo,yr,pd)* (h(rev,yr,pd)-h0(pwst))* pw_effy(pwst)*2.725)) =e= power(pwst,yr,pd); dm_tfla(agdm,yr,pd) to_infa(agdm,yr,pd) =e= (SUM(rn$NDLINKa(rn,agdm),N_Da(rn,agdm,yr,pd)) + SUM(rev$rdlinka(rev,agdm),RES_Da(rev,agdm,yr,pd)))*(1-lossa(agdm)); dm_tflm(mundm,yr,pd) to_infm(mundm,yr,pd) =e= (SUM(rn$NDLINKm(rn,mundm),N_Dm(rn,mundm,yr,pd)) + SUM(rev$rdlinkm(rev,mundm),RES_Dm(rev,mundm,yr,pd)))*(1-lossm(mundm)); dm_tfli(indm,yr,pd) to_infi(indm,yr,pd) =e= (SUM(rn$NDLINKi(rn,indm),N_Di(rn,indm,yr,pd)) + SUM(rev$rdlinki(rev,indm),RES_Di(rev,indm,yr,pd)))*(1-lossi(indm)); w_bala(agdm,yr,pd) SUM(crop$(cp_prd(agdm,crop,pd) and dem_cp(agdm,crop)), wsfacp(agdm,crop,yr,pd))=e= to_infa(agdm,yr,pd); cal_wacp(agdm,crop,yr,pd)$(dem_cp(agdm,crop) and cp_prd(agdm,crop,pd)) wacp(agdm,crop,yr,pd) =e= wsfacp(agdm,crop,yr,pd); - 127 - cal_agwto(agdm,yr) agwto(agdm,yr) =e= sum(rn,sum(pd,N_Da(rn,agdm,yr,pd))) + sum(rev,sum(pd,Res_Da(rev,agdm,yr,pd))); cal_twacp(agdm,crop,yr)$dem_cp(agdm,crop) twacp(agdm,crop,yr) =e= SUM(pd$cp_prd(agdm,crop,pd),wacp(agdm,crop,yr,pd) + area(agdm,crop,yr)*er(crop,pd)*0.00001); area_lmt1(agdm,yr,pd) SUM(crop$cp_prd(agdm,crop,pd),area(agdm,crop,yr)) =l= 1.2*tarea(agdm); cal_dp_stg(agdm,crop,yr,pd)$dem_cp(agdm,crop) dp_stg(agdm,crop,yr,pd) =e= wacp(agdm,crop,yr,pd)*(1 - effy(agdm,crop)); cal_dp(agdm,crop,yr)$dem_cp(agdm,crop) dp(agdm,crop,yr) =e= sum(pd,dp_stg(agdm,crop,yr,pd)); cal_eta_se(agdm,crop,yr) eta_se(agdm,crop,yr) =l= twacp(agdm,crop,yr) - dp(agdm,crop,yr); cal_effy(agdm,crop,yr) effy(agdm,crop)*SUM(pd$cp_prd(agdm,crop,pd),wacp(agdm,crop,yr,pd)) =e= eta_se(agdm,crop,yr) - ter(agdm,crop)*area(agdm,crop,yr)*0.00001; cal_et_st(agdm,crop,yr,pd)$(dem_cp(agdm,crop) and cp_prd(agdm,crop,pd)) eta_st(agdm,crop,yr,pd) =e= wacp(agdm,crop,yr,pd)*effy(agdm,crop)+ area(agdm,crop,yr)*er(crop,pd)*0.00001; eta_lmt(agdm,crop,yr,pd)$cp_prd(agdm,crop,pd) eta_st(agdm,crop,yr,pd) =l= etm(crop,pd)*area(agdm,crop,yr)*0.00001; eta_se_st(agdm,crop,yr) eta_se(agdm,crop,yr) =e= sum(pd$cp_prd(agdm,crop,pd),eta_st(agdm,crop,yr,pd)); cal_rf_a(agdm,yr,pd) retn_fla(agdm,yr,pd) =e= SUM(crop$dem_cp(agdm,crop),dp_stg(agdm,crop,yr,pd)*drna(agdm)) + (SUM(rn$NDLINKa(rn,agdm),N_Da(rn,agdm,yr,pd)) +SUM(rev$rdlinka(rev,agdm),RES_Da(rev,agdm,yr,pd)))*lossa(agdm)*0.40; cal_rf_m(mundm, yr,pd) retn_flm(mundm, yr,pd) =e= ( SUM(rn$NDLINKm(rn,mundm),N_Dm(rn,mundm,yr,pd)) + SUM(rev$RDLINKm(rev,mundm),RES_Dm(rev,mundm,yr,pd)))*lossm(mundm)*0.50 + (to_infm(mundm,yr,pd) *drnm(mundm)*0.5) + (to_infm(mundm,yr,pd)*0.9*(1-consfrm(mundm))*0.8); cal_rf_i(indm,yr,pd) retn_fli(indm,yr,pd) =e= - 128 - (SUM(rn$NDLINKi(rn,indm), N_Di(rn,indm,yr,pd)) + SUM(rev$RDLINKi(rev, indm), RES_Di(rev,indm,yr,pd)))*lossi(indm)*0.50 + to_infi(indm,yr,pd)*drni(indm)*0.5 + to_infi(indm,yr,pd)*(1- drni(indm))*(1-consfri(indm))*0.8; def_dft(agdm,crop,yr,pd)$cp_prd(agdm,crop,pd) dft(agdm,crop,yr,pd)*(area(agdm,crop,yr)*0.00001*etm(crop,pd)) =e= kym(crop,pd)*((area(agdm,crop,yr)*0.00001*etm(crop,pd)) - eta_st(agdm,crop,yr,pd)); def_mdft1(agdm,crop,yr,pd)$cp_prd(agdm,crop,pd) mdft(agdm,crop,yr) =g= dft(agdm,crop,yr,pd); yldm_lmt(agdm,crop,yr) ylda(agdm,crop,yr) =g= 0.5*yldst(agdm,crop); cal_yldst(agdm,crop)$dem_cp(agdm,crop) yldst(agdm,crop) =e= SUM(IN,ALPHA(crop,IN)*INPTS(agdm,crop,IN)) -SUM(IN,INPTS(agdm,crop,IN)*SUM(IIN,ZETA(crop,IN,IIN)*INPTS(agdm,crop,IIN))); cal_ylda(agdm,crop,yr)$dem_cp(agdm,crop) ylda(agdm,crop,yr)=e= yldst(agdm,crop)*(1-ky(crop)*(1-eta_se(agdm,crop,yr) /SUM(pd$cp_prd(agdm,crop,pd),area(agdm,crop,yr)*0.00001*etm(crop,pd)))); cal_muwatsf(mundm,yr,pd) muwatsf(mundm,yr,pd) =e= to_infm(mundm,yr,pd); cal_muwat(mundm,yr,pd) tmuwat(mundm,yr,pd) =e= (to_infm(mundm,yr,pd) *(1-drnm(mundm))) ; ch_rur(mundm,yr,pd) MGWR(mundm,pd)*tmuwatr(mundm,yr,pd)+ MSFR(mundm,pd)*tmuwatr(mundm,yr,pd) =e=tmuwatr(mundm,yr,pd); ch_urb(mundm,yr,pd) MGWU(mundm,pd)*tmuwatu(mundm,yr,pd)+ MSFU(mundm,pd)*tmuwatu(mundm,yr,pd) =e=tmuwatu(mundm,yr,pd); tmuwat_lmt(mundm,yr,pd) tmuwat(mundm,yr,pd) =e= (tmuwatu(mundm,yr,pd)*MGWU(mundm,pd) +tmuwatu(mundm,yr,pd)*MSFU(mundm,pd)) +(tmuwatr(mundm,yr,pd)*MGWR(mundm,pd) +tmuwatr(mundm,yr,pd)*MSFR(mundm,pd)); cal_inwatsf(indm,yr,pd) inwatsf(indm,yr,pd) =e= to_infi(indm,yr,pd); cal_inwat(indm,yr,pd) tinwat(indm,yr,pd) =e= to_infi(indm,yr,pd)*(1-drni(indm)); cal_muvalu(mundm)$urban(mundm) muvalu(mundm)$urban(mundm) =e= sum(yr,sum(pd, - 129 - exp(intp1/elas_pr1)* incu(mundm)**(elas_in1/elas_pr1)* 1/(1-1/elas_pr1)* (((tmuwatu(mundm,yr,pd)+0.0001)*1000000 /popu(mundm))**(-1/elas_pr1+1) -(0.9*mwdu(mundm,pd))**(-1/elas_pr1+1))*popu(mundm)/1000000 to_infm(mundm,yr,pd)$urban(mundm)*(1-drnm(mundm))*MSFU(mundm,pd)*w_cmu(mundm))); cal_muvalr(mundm)$rural(mundm) muvalr(mundm)$rural(mundm) =e=sum(yr,sum(pd, exp(intp2/elas_pr2)* incr(mundm)**(elas_in2/elas_pr2)* 1/(1-1/elas_pr2)* (((tmuwatr(mundm,yr,pd)+0.0001)*1000000 /popr(mundm))**(-1/elas_pr2+1) -(0.9*mwdr(mundm,pd))**(-1/elas_pr2+1)) * popr(mundm)/1000000 -to_infm(mundm,yr,pd)*(1-drnm(mundm))*MSFR(mundm,pd)* w_cmr(mundm))); cal_muval(mundm) muval(mundm) =e=muvalu(mundm)$urban(mundm) + muvalr(mundm); cal_arval(agdm) arval(agdm) =e= sum(yr,sum(crop$dem_cp(agdm,crop), area(agdm,crop,yr)*(price(crop)*ylda(agdm,crop,yr) - (INPTS(agdm,crop,'labor')*inpcos(crop,'labc') + INPTS(agdm,crop,'fertn')*inpcos(crop,'fertnc')+ INPTS(agdm,crop,'fertp')*inpcos(crop,'fertpc')+ INPTS(agdm,crop,'fertk')*inpcos(crop,'fertkc')+ INPTS(agdm,crop,'mach') + INPTS(agdm,crop,'pesti')))/1000000 -(0.7* w_ca(crop)*SUM(pd$cp_prd(agdm,crop,pd),wsfacp(agdm,crop,yr,pd))))); cal_inval(indm) inval(indm) =e= sum(yr,sum(pd, (exp(alfa(indm)/elas_ind)/(1-1/elas_ind)* ((tinwat(indm,yr,pd)*1000000)**(-1/elas_ind+1) - (iwd0(indm))**(-1/elas_ind+1))/1000000) - to_infi(indm,yr,pd) *w_ci(indm)*(1-drni(indm)))) ; cal_powval(pwst) powval(pwst) =e= (sum(yr,sum(pd,(pben(pwst)-pcost(pwst))*power(pwst,yr,pd))))/1000; *** MODEL DONGNAI / objfun N_BAL R_BAL hydr_rev - 130 - w_bala dm_tfla dm_tflm dm_tfli cal_ylda cal_yldst cal_muval cal_muvalu cal_muvalr cal_inval cal_arval cal_powval miflow_lt tmuwat_lmt yldm_lmt def_dft def_mdft1 cal_twacp cal_wacp cal_et_st cal_dp cal_dp_stg cal_agwto cal_rf_a cal_rf_m cal_rf_i eta_lmt area_lmt1 cal_eta_se cal_muwatsf cal_muwat cal_inwat cal_inwatsf ar_vol vol_elev eta_se_st ch_urb ch_rur / ; {select a solver} option nlp=conopt3; dongnai.optfile=1; SOLVE DONGNAI USING NLP MAXIMIZING OBJ; *** *** - - 131 - LÝ LỊCH TRÍCH NGANG Họ tên: NGUYỄN VŨ HUY Ngày, tháng, năm sinh: 25-09-1977 Nơi sinh: NAM ĐỊNH Địa liên lạc: Số 17 đường 15, Phường Bình An, Quận 2, TPHCM Điện thoại : 8.875.294 – 0913.985.041 E-mail : huy77nguyen@yahoo.com Q TRÌNH ĐÀO TẠO 1995-2000: Học đại học hệ qui tập trung trường Đại học Bách khoa TP Hồ Chí Minh, tốt nghiệp văn Kĩ sư Xây dựng năm 2000 2002-2004 : Học cao học ngành Xây dựng Cơng trình thủy Đại học Bách khoa TP Hồ Chí Minh Q TRÌNH CƠNG TÁC Từ năm 2000 đến : công tác Phân viện Khảo sát Qui hoạch Thủy lợi Nam Bộ, chức vụ kĩ sư qui hoạch thủy lợi thuộc phòng Qui hoạch thủy lợi lưu vực sông Đồng Nai Địa quan: 271/3 An Dương Vương, Phường 3, Quận 5, TPHCM ... quiết vấn đề quản lí nguồn nước lưu vực sông Đồng Nai lưu vực khác nước vấn đề cần quan tâm nghiên cứu nhà khoa học Luận văn thạc só ? ?Quản lí tối ưu nguồn nước lưu vực sông Đồng Nai? ?? thực khoá 13... văn đưa kiến nghị quản lí tổng hợp nguồn nước cho lưu vực sông Đồng Nai phương pháp tiếp cận quản lí nguồn nước cho lưu vực sông khác -7- ABSTRACT Water resources in Dong Nai river basin play... nguồn nước giữ vai trò quan trọng để đảm bảo đạt mục tiêu phát triển kinh tế bảo vệ môi trường chống suy thoái nguồn nước Lưu vực sông Đồng Nai (LVSĐN) lưu vực sông lớn thứ ba Việt Nam lại lưu