Luận án tiến sĩ nghiên cứu xác định tiến độ thi công công trình có tính đến yếu tố bất định áp dụng đối với một số

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI TRẦN HỮU LÂN NGHIÊN CỨU XÁC ĐỊNH TIẾN ĐỘ THI CƠNG CƠNG TRÌNH CĨ TÍNH ĐẾN YẾU TỐ BẤT ĐỊNH ÁP DỤNG ĐỐI VỚI MỘT SỐ CÔNG TRÌNH Ở VIỆT NAM Chuyên ngành : Kinh tế học Mã số : 62.31.03.01 LUẬN ÁN TIẾN SĨ KINH TẾ NGƯỜI HƯỚNG DẪN KHOA HỌC PGS TS NGUYỄN MINH DUỆ Hà Nội – 2012 LỜI CAM ĐOAN Tôi tác giả luận văn xin cam đoan: Đây cơng trình nghiên cứu riêng hướng dẫn khoa học PGS.TS Nguyễn Minh Duệ Các tài liệu, số liệu sử dụng luận án có nguồn gốc rõ ràng, công bố theo quy định Các kết nghiên cứu luận án tác giả tự nghiên cứu độc lập, trung thực, không chép Hà Nội, ngày 26 tháng 06 năm 2012 Người cam đoan Trần Hữu Lân LỜI CẢM ƠN Tôi xin bày tỏ lòng biết ơn sâu sắc PGS TS Nguyễn Minh Duệ, người thầy người hướng dẫn khoa học tận tình hướng dẫn, động viên, khích lệ giúp đỡ suốt thời gian thực luận án Tôi xin chân thành cảm ơn cán thuộc Viện đào tạo Bồi dưỡng sau đại học, Trường Đại học Bách khoa Hà Nội tận tình, tạo điều kiện thuận lợi nghiên cứu thực thủ tục quy định cho suốt thời gian nghiên cứu học tập vừa qua Tơi xin bày tỏ lịng biết ơn thầy cô giáo thuộc Viện Kinh tế Quản lý, Trường Đại học Bách khoa Hà Nội; viện nghiên cứu, ban quản lý dự án, doanh nghiệp, đặc biệt Tạp chí Kinh tế Dự báo; Cơng ty TNHH Minh Tuấn, đặc biệt Giám đốc Nguyễn Đình Bích; Công ty TNHH Năng lượng mới, Giám đốc Nguyễn Hùng Minh; Ban quản lý dự án xây dựng Đà Nẵng sẵn lòng giúp đỡ hợp tác để vận dụng kết nghiên cứu vào thực tiễn thời gian qua tương lai Tôi đồng thời xin bày tỏ cảm ơn sâu sắc tới KS Phạm Thành Long, người anh, người bạn TS KTS Nguyễn Đức Anh, thuộc Viện khoa học Công nghệ quốc gia Đài Loan ý kiến gợi mở ban đầu cho việc hình thành đề tài Tôi bày tỏ cảm ơn sâu sắc đến lời động viên, giúp đỡ anh Lê Anh Thư, Nguyễn Đức Hùng, Chu Tấn Phát Tôi xin bày tỏ lòng biết ơn đến vợ, con, cha mẹ hai bên nội ngoại tạo điều kiện thuận lợi tốt vật chất tinh thần suốt trình thực luận án Đà Nẵng Thủ đô Hà Nội Cuối cùng, vô biết ơn Ban Giám đốc, Phòng Kinh tế Ngành, Phòng Tổ chức, Văn phòng, Phòng Xây dựng – Thẩm định đồng nghiệp thuộc Sở Kế hoạch Đầu tư thành phố Đà Nẵng tạo điều kiện quý báu cho mặt pháp lý thời gian mà thiếu tơi khơng thể hồn thành luận án này./ Trần Hữu Lân MỤC LỤC Trang phụ bìa Trang Lời cam đoan i Lời cảm ơn ii Mục lục iii Danh mục kí hiệu chữ viết tắt iv Danh mục bảng V Danh mục sơ đồ, hình vẽ đồ thị Vi PHẦN MỞ ĐẦU 17 Sự cần thiết đề tài 17 Mục tiêu nghiên cứu 18 Đối tượng phạm vi nghiên cứu 18 3.1 Đối tượng nghiên cứu 18 3.2 Phạm vi nghiên cứu 19 Phương pháp nghiên cứu 19 Những kết đạt được, ý nghĩa khoa học, ý nghĩa thực tiễn, điểm luận án 19 5.1 Những kết đạt điểm luận án 19 5.2 Ý nghĩa khoa học 20 5.3 Ý nghĩa thực tiễn 20 Các nghiên cứu trước có quan 21 Bố cục luận án 22 Chương MỘT SỐ VẤN ĐỀ LÝ LUẬN CƠ BẢN VỀ XÁC ĐỊNH TIẾN ĐỘ THI CÔNG CÔNG TRÌNH CĨ TÍNH ĐẾN YẾU TỐ BẤT ĐỊNH 23 1.1 Tiến độ thi công cần thiết xác định tiến độ thi công 23 1.1.1.Tiến độ thi cơng 23 1.1.2 Sự cần thiết xác định tiến độ thi công 24 1.1.3 Yêu cầu xác định tiến độ thi công 24 1.2 Nội dung chủ yếu xác định tiến độ thi công 25 1.3 Các phương pháp xác định tiến độ thi công 29 1.3.1 Nhóm phương pháp tiến độ theo đường thẳng 30 1.3.1.1 Phương pháp tiến độ ngang (Gantt Charts Method) 30 1.3.1.2 Phương pháp đường tiến độ xiên 31 1.3.2 Nhóm phương pháp tiến độ theo sơ đồ mạng 33 1.3.2.1 Giới thiệu chung 33 1.3.2.2 Trình tự lập sơ đồ mạng 34 1.3.2.3 Phương pháp sơ đồ mạng CPM (Critical Path Method) 36 1.3.2.4 Phương pháp sơ đồ mạng PERT (Program Evaluation Review Technique) 40 1.3.3 Phương pháp quản trị giá trị thu EVM (Earned Value Management) 42 1.3.4 Phương pháp Mô Monte Carlo (Monte Carlo Simulation) 44 1.3.5 Đánh giá phương pháp xác định tiến độ thi công 45 1.3.5.1 Đối với nhóm phương pháp tiến độ theo đường thẳng 45 1.3.5.2 Đối với nhóm phương pháp tiến độ sơ đồ mạng 46 1.3.5.3 Đối với phương pháp quản trị giá trị thu (EVM) 47 1.4 Áp dụng xác định tiến độ thi công nước ta 48 1.4.1 Quan điểm tiến độ thi công chủ đầu tư 48 1.4.2 Quan điểm tiến độ thi công nhà thầu 48 1.4.3 Mối liên hệ quan điểm xác định tiến độ 49 1.5 Một số nhận xét xác định tiến độ thi công nước ta 50 1.5.1 Một số nhận xét 50 1.5.2 Kết luận hướng nghiên cứu 51 Tóm tắt chương 53 Chương XÁC ĐỊNH TIẾN ĐỘ THI CƠNG CƠNG TRÌNH CĨ TÍNH ĐẾN YẾU TỐ BẤT ĐỊNH THEO PHƯƠNG PHÁP DỰ BÁO XÁC SUẤT 54 2.1 Yếu tố bất định tác động đến thực tiến độ thi công 54 2.1.1 Thế yếu tố bất định 54 2.1.1.1 Khái niệm 54 2.1.1.2 Nhận dạng yếu tố bất định theo mức độ tác động 54 2.1.1.3 Nhận dạng yếu tố bất định theo nguyên nhân 55 2.1.2 Tác động yếu tố bất định đến tiến độ thi công 55 2.1.2.1 Tại phải xem xét yếu tố bất định 55 2.1.2.2 Các minh chứng điển hình tác động yếu tố bất định 55 2.1.3 Các yếu tố bất định liên quan đến tiến độ thi công 56 2.1.3.1 Nhóm yếu tố bất định bên ngồi (ký hiệu βni) 56 2.1.3.2 Nhóm yếu tố bất định bên (ký hiệu βti) 56 2.1.3.3 Xác định thời gian yếu tố bất định tác động lên tiến độ 57 2.2 Mơ hình đường cong tiến độ khối lượng thi cơng biến thiên có tác động yếu tố bất định 58 2.2.1 Mơ hình đường cong tiến độ 58 2.2.2 Sự biến thiên đường cong tiến độ có tác động yếu tố bất định 58 2.2.3 Mục đích việc xác định mơ hình 60 2.3 Nội dung phương pháp dự báo tiến độ thi công có xét đến yếu tố bất định Phương pháp dự báo xác suất KALMAN 61 2.3.1 Tổng quan Dự báo Kalman (KMDB) 61 2.3.2 Nội dung Dự báo Kalman 62 2.4 Phương pháp dự báo xác suất Kalman việc áp dụng dự báo tiến độ thi công công trình có tính đến yếu tố bất định 67 2.4.1 Cơ sở tính khả thi việc áp dụng 68 2.4.1.1 Cơ sở áp dụng 68 2.4.1.2 Tính khả thi việc áp dụng 68 2.4.2 Phương pháp dự báo xác suất Kalman dự báo tiến độ thi công 69 2.4.2.1 Điều kiện dự án để áp dụng Dự báo Kalman 69 2.4.2.2 Mô hình dự báo tiến độ 69 2.4.2.3 Nội dung dự báo tiến độ 70 2.4.2.4 Xác định phân phối xác suất thời điểm hồn thành ρ(T) 77 2.5 Mơ hình chi tiết liệu đầu vào đầu cho dự báo tiến độ áp dụng Phương pháp dự báo xác suất Kalman 80 2.5.1 Mơ hình chi tiết liệu đầu vào 80 2.5.2 Mô hình chi tiết liệu đầu 82 2.6 Xây dựng thuật toán ứng dụng máy tính 83 2.6.1 Ngơn ngữ viết lập trình phần mềm quản lý dự án thông dụng 83 2.6.2 Xây dựng sơ đồ hệ thống 84 2.6.3 Xây dựng sơ đồ khối 85 2.6.4 Chú giải sơ đồ hệ thống sơ đồ khối 86 Tóm tắt chương 89 Chương ÁP DỤNG PHƯƠNG PHÁP DỰ BÁO XÁC SUẤT KALMAN XÁC ĐỊNH TIẾN ĐỘ MỘT SỐ CƠNG TRÌNH THỰC TẾ Ở VIỆT NAM 90 3.1 Giới thiệu 90 3.2 Áp dụng dự báo tiến độ thi công cho cơng trình lúc nhiều cơng trình thực tế Việt Nam 90 3.2.1 Áp dụng thứ nhất, dự báo tiến độ cho cơng trình độc lập 90 3.2.1.1 Các thông tin đầu vào chủ yếu 91 3.2.1.2 Xác định tiến độ dự án 91 3.2.1.3 Giải thích kết dự báo giải cách tính 96 3.2.1.4 Nhận xét 98 3.2.2 Áp dụng thứ hai, xác định tiến độ thi cơng đồng thời ba dự án có thời gian khởi công khác 99 3.2.2.1 Các thông tin đầu vào chủ yếu 99 3.2.2.2 Xác định tiến độ dự án 103 3.2.2.3 Giải thích kết dự báo giải cách tính 111 3.2.2.4 Nhận xét áp dụng thứ hai 112 3.2.3 Áp dụng thứ ba, xác định tiến độ thi công cơng trình nhiều phương pháp khác 113 3.2.3.1 Các thông tin đầu vào chủ yếu 113 3.2.3.2 Xác định tiến độ 114 3.2.3.3 Giải thích kết dự báo giải cách tính 120 3.2.3.4 Nhận xét áp dụng thứ ba 122 3.3 So sánh đánh giá phương pháp dự báo tiến độ thi cơng có xét đến yếu tố bất định với phương pháp truyền thống 123 3.3.1 So sánh phương pháp dự báo tiến độ 123 3.3.2 Đánh giá ưu việt khả dự báo tiến độ áp dụng Phương pháp dự báo xác suất Kalman 125 Tóm tắt chương 128 KẾT LUẬN 130 DANH MỤC CÁC CƠNG TRÌNH ĐÃ CÔNG BỐ CỦA LUẬN ÁN 133 DANH MỤC TÀI LIỆU THAM KHẢO 134 Tài liệu Tiếng Việt 134 Tài liệu tiếng nước 137 Tài liệu tham khảo Website 140 PHẦN PHỤ LỤC 142 DANH MỤC CÁC KÍ HIỆU VÀ CHỮ VIẾT TẮT a Hệ số Góc phương trình Gompertz A Ma trận chuyển trạng thái, có kích thước (n x n) AC Chi phí thực tế phương pháp EVM (Actual Cost) B Ma trận điều khiển, có kích thước (n x l) BAC Tổng chi phí theo kế hoạch thời điểm hoàn thành phương pháp EVM (Budgeted At completion) CPM Phương pháp sơ đồ mạng (Critical path method) di Định mức thời gian cơng việc i E(tij) Kỳ vọng tốn học phương pháp sơ đồ PERT ES(t) Thời gian theo tiến độ kế hoạch phương pháp EVM (Earned schedule) EV Giá trị thu được, phương pháp EVM (Earned Value) EVM Phương pháp quản trị giá trị thu (Earned Value Management) GANTT Phương pháp sơ đồ Gantt (Gantt Chart Method) H Ma trận quan sát, có kích thước (m x n) KMDB Phương pháp dự báo xác suất Kalman Kk Gia số Kalman m Hệ số chặn phương trình Gompertz N(ci ) Số nhân lực thực công việc i P0 Hiệp phương sai sai số giá trị x0 ρt(k) Phân phối xác suất trước biến thời gian k ρs(k) Phân phối xác suất sau biến thời gian k ρ(T) Phân phối xác suất thời điểm hồn thành PV Chi phí kế hoạch phương pháp EVM (Planned Value) Q Ma trận hiệp phương sai nhiễu trình Qi Khối lượng thực công việc i R Ma trận hiệp phương sai nhiễu phép đo ST Giá trị hàm w(t) thời điểm T Phương trình Gompertz tk Biến thời gian dự báo thời điểm hồn thành cơng trình Phương pháp dự báo xác suất Kalman tk(t) Giá trị thời gian dự báo thực hồn thành cơng trình, phương pháp EVM tm Thời gian thường gặp thực hồn thành cơng việc i Thời gian bi quan thực hồn thành cơng việc i tei Thời gian dự tính thực hồn thành cơng việc i Tcp Tiến độ kế hoạch chủ đầu tư TDB Thời gian hồn thành cơng trình theo dự báo TĐHT Thời gian hồn thành cơng trình theo kế hoạch TĐTC Tiến độ thi cơng TĐTX Điểm 100% tiến độ hồn thành Ti Thời gian thi công công việc i Tj Thời gian thi công công việc j Ttt Thời gian thực tế thi công đến thời điểm dự báo TVk Vector biến dự báo tiến độ t kts ij Thời gian kết thúc sớm sơ đồ mạng CPM t bds ij Thời gian bắt đầu sớm sơ đồ mạng CPM 10 12 54.260 45.598 0.24 0.00 2% 2.54 0.000 32.7 13 57.790 49.016 0.36 0.00 1% 2.40 0.000 32.8 14 57.790 52.421 0.50 0.01 0% 2.41 0.000 33.8 15 61.310 0.64 0.02 16 68.340 0.76 0.06 18 75.380 0.86 0.12 19 78.890 0.92 0.23 20 83.860 0.96 0.37 21 84.500 0.98 0.54 22 87.600 0.99 0.69 23 88.987 1.00 0.82 24 90.231 1.00 0.91 25 92.456 1.00 0.96 26 98.960 1.00 0.98 27 99.480 1.00 1.00 28 100.000 1.00 1.00 145 Phụ lục 2.2 Bảng dự báo tiến độ theo Phương pháp dự báo xác suất Kalman Áp dụng thứ 2, Dự án 3: Cầu qua Vịnh Đà Nẵng Tiến độ (%) Tháng Phân phối xác suất (100%) Trước Sau ρt(k) ρs(k) D E 0.00 0.00 Xác suất hoàn thành ρ(T) Độ lệch chuẩn F 50% G 3.00 H 0.000 I 30.0 Tác động Giá trị TĐTX (tháng) bất định ∆Ti (tháng) Kế hoạch Thực tế A B 0.000 C 0.000 0.000 2.736 0.00 0.00 50% 3.00 0.000 30.0 5.830 7.034 0.00 0.00 50% 2.85 0.000 30.0 11.270 11.141 0.01 0.00 50% 2.85 0.500 30.0 18.510 15.299 0.02 0.00 48% 2.98 0.300 30.2 25.800 18.967 0.04 0.00 36% 3.08 0.600 31.1 29.590 22.215 0.08 0.00 17% 3.31 1.200 33.1 32.950 25.670 0.15 0.00 5% 3.50 0.900 35.9 35.510 29.513 0.24 0.00 2% 3.83 0.000 38.2 43.480 32.958 0.36 0.00 1% 3.83 0.000 39.5 10 47.040 36.257 0.50 0.00 0% 3.68 0.000 39.7 11 52.392 0.64 0.01 12 59.400 0.76 0.02 146 Ghi K 13 64.408 0.85 0.04 14 69.416 0.92 0.06 15 71.920 0.96 0.10 16 76.912 0.98 0.16 17 79.432 0.99 0.24 18 80.920 1.00 0.33 19 82.650 1.00 0.43 20 84.560 1.00 0.54 21 86.980 1.00 0.64 22 87.340 1.00 0.74 23 88.810 1.00 0.82 24 89.980 1.00 0.88 25 90.448 1.00 0.93 26 91.952 1.00 0.96 27 94.456 1.00 0.98 28 96.960 1.00 0.99 29 99.480 1.00 0.99 30 100.000 1.00 1.00 147 Phụ lục Bảng đối chiếu thời gian dự báo lúc cho cơng trình đồ thị 3.2, áp dụng Thứ tự mốc thời gian (tháng) 10 11 12 13 14 15 16 17 18 19 20 Dự án 1 10 11 12 13 14 15 16 17 18 19 20 Tiến độ kế hoạch 0.0 0.3 1.0 3.6 8.2 12.3 19.0 23.1 28.3 32.4 37.5 40.8 45.2 52.5 57.1 65.7 73.7 80.1 85.4 90.8 94.0 Tiến độ thực tế 0.0 0.0 0.0 4.1 8.0 11.6 16.6 21.6 25.9 28.7 33.4 36.5 37.7 43.1 47.9 57.1 60.8 65.5 68.8 Dự án Dự án khởi công 10 11 12 13 14 15 16 Tiến độ kế hoạch 0.0 0.0 6.3 14.2 17.3 25.1 29.9 33.5 36.1 40.0 43.6 47.2 54.3 57.8 57.8 61.3 68.3 Tiến độ thực tế 0.0 3.2 7.9 12.7 16.7 20.5 24.3 27.9 31.4 34.9 38.6 42.2 45.6 49.0 52.4 Dự án Dự án khởi công 10 11 Tiến độ kế hoạch 0.0 0.0 5.8 11.3 18.5 25.8 29.6 33.0 35.5 43.5 47.0 52.4 59.4 Tiến độ thực tế 0.0 2.7 7.0 11.1 15.3 19.0 22.2 25.7 29.5 33.0 36.3 Thời điểm dự báo 03 dự án 148 12 Phụ lục 4.1: Bảng tổng tiến độ Trung tâm giáo dục thường xun huyện Hịa Vang, phần Máy thi cơng BiỂU ĐỒ MÁY THI CƠNG CƠNG TRÌNH: TRUNG TÂM GIÁO DỤC THƯỜNG XUYÊN HUYỆN HÒA VANG ẠI THIẾT Đ LO PH MÁY ĐÀO ĐẤT 10 11 B 10 c Ự ĐỔ 1 2 2 3 2 Cái Cái Cái Cái Cái Cái Cái ẦN THÂN MÁY ĐẦM DÙI ỘN VỮA MÁY TR MÁY T ỜI Cái Cái Cái 2 2 3 2 ẬN THĂNG MÁY V ỘN VỮA MÁY TR ẮT GẠCH MÁY C MÁY TOÀN ĐẠC MÁY T ỜI MÁY K HOAN MÁY HÀN ĐIỆN Cái Cái Cái Cái Cái Cái Cái 5 MÁY B ẮT UỐN THÉP MÁY C MÁY HÀN ĐIỆN MÁY TOÀN ĐẠC MÁY ĐẦM B ÀN PH 10 ẬN THĂNG ẠM TRỘN BÊ TÔ 20M3/H ƠM BÊ TÔNG MÁY V Cái ẠM TRỘN BÊ TÔ TR 20M3/H ẬN CHUY Ô TÔ V MÁY ĐẦM ĐẤT TACOM ẮT UỐN THÉP MÁY C MÁY HÀN ĐIỆN MÁY TOÀN ĐẠC MÁY ĐẦM B ÀN MÁY ĐẦM DÙI ỘN VỮA MÁY TR TR THÁNG 12/2007 Cái Cái Cái Cái Cái Cái Cái Cỏi Cỏi Cỏi PH Năm 2008 S l 15 20 THÁNG 01/2008 25 31 10 15 1 4 3 1 2 1 2 THÁNG 02/2008 THÁNG 03/2008 20 25 31 10 15 20 25 1 2 1 2 3 1 2 3 1 2 3 2 2 2 1 2 3 2 1 2 3 2 1 2 3 2 3 2 28 THÁNG 04/2008 THÁNG 5/2008 10 15 20 25 31 10 15 20 25 30 5 2 1 1 1 5 1 5 1 5 2 5 2 5 1 5 1 5 1 5 1 5 ẦN MĨNG A Ơ TƠ T ơn vị ẦN HỒN TH Số lượng máy móc, thiết bị (cái) stt 2 2 2 2 2 2 149 2 2 10 Phụ lục 4.2: Bảng tổng tiến độ Trung tâm giáo dục thường xuyên Huyện Hòa Vang 150 Biểu đồ nhân lực Nhân công Thời gian 151 Phụ lục 4.3 Bảng dự báo tiến độ theo Phương pháp dự báo xác suất Kalman Dự án: Trung tâm giáo dục thường xuyên huyện Hòa Vang Tiến độ (%) Tháng Phân phối xác suất (100%) Trước Sau ρt(k) ρs(k) D E 0.00 0.00 Xác suất hoàn thành ρ(T) Độ lệch chuẩn F 50% G 3.30 H 0.000 I 33.0 Tác động bất Giá trị TĐTX (tháng) định ∆Ti (tháng) Kế hoạch Thực tế A B 0.000 C 0.000 2.990 5.784 0.00 0.00 47% 3.15 0.000 33.0 9.260 8.987 0.00 0.00 51% 3.30 0.000 33.2 9.260 12.204 0.00 0.00 58% 3.15 0.000 32.9 12.570 15.346 0.01 0.00 61% 3.13 0.000 32.4 15.800 18.197 0.01 0.00 48% 3.09 0.500 32.2 22.863 21.084 0.03 0.00 46% 3.17 0.600 33.1 24.737 23.860 0.06 0.00 34% 3.07 1.000 33.3 30.713 26.917 0.10 0.00 18% 3.09 0.500 34.3 35.627 29.893 0.17 0.00 7% 3.09 0.500 35.8 10 44.889 32.646 0.26 0.00 2% 3.24 1.000 37.8 11 47.583 35.904 0.38 0.00 1% 3.38 0.000 39.7 12 49.962 38.938 0.50 0.01 0% 3.69 0.000 42.4 13 52.150 0.62 0.01 152 Ghi K 14 62.580 0.74 0.02 15 64.110 0.83 0.04 16 67.300 0.90 0.07 17 68.100 0.94 0.12 18 70.090 0.97 0.18 19 72.020 0.99 0.26 20 73.080 0.99 0.35 21 75.070 1.00 0.46 22 76.080 1.00 0.57 23 77.121 1.00 0.67 24 78.123 1.00 0.76 25 79.547 1.00 0.84 26 80.500 1.00 0.89 27 82.060 1.00 0.94 28 90.800 1.00 0.96 29 94.000 1.00 0.98 30 97.600 1.00 0.99 31 98.500 1.00 1.00 32 99.900 1.00 1.00 33 100.000 1.00 1.00 153 Phụ lục Diễn giải chi tiết cách tính tốn theo phần Code lập trình phần mềm dùng cho việc trợ giúp máy tính + Khai báo nhập giá trị ma trận A, H, Q, R: + Ma trận A xác định theo công thức (2.22), trường hợp kỳ báo cáo đơn vị thời gian, ma trận A=1 + Ma trận H xác định theo công thức (2.25), trình quan sát đồng cho trình dự báo ma trận H=1 + Ma trận Q xác định theo công thức (2.43), trường hợp kỳ báo cáo đơn vị thời gian, ma trận Q=1 + Ma trận R xác định theo công thức (2.26) Trong trường hợp này, sai số phương pháp đo dụng cụ đo không, giá trị Rk phụ thuộc vào tác động yếu tố bất định bước k + Xác định tác động bất định theo công thức (2.39) (2.40) + Khai báo biến : dtResults = new DataTable(); dtResults.Columns.Add("Kỳ", typeof(int)); dtResults.Columns.Add("Trị trung bình", typeof(double)); dtResults.Columns.Add("Độ lệch chuẩn", typeof(double)); dtResults.Columns.Add("Xác suất P(t)", typeof(string)); dtXS = new DataTable(); dtXS.Columns.Add("Kỳ", typeof(int)); dtXS.Columns.Add("Xác suất trước", typeof(double)); dtXS.Columns.Add("Xác suất sau", typeof(double)); + Xác định giá trị tiền nghiệm (XT), hiệp phương sai sai số trước (PT) theo công thức (2.27), (2.31) 154 + Xác định giá trị tiền nghiệm theo công thức (2.31) XP = A*X0 MATMUL(ref TBUOC[k].XM, TRANGTHAI, X0); + Xác định hiệp phương sai sai số trước theo công thức (2.27) PT = A*P0*A' + Q MATQUA(ref TBUOC[k].PM, TRANGTHAI, P0); MATADD(ref TBUOC[k].PM, TBUOC[k].PM, PROC_BATDINH); + Xác định tính tốn lại theo bước thời gian k: XP = A*X MATMUL(ref TBUOC[k].XM, TRANGTHAI, TBUOC[k - 1].X); PT = A*P*A' + Q MATQUA(ref TBUOC[k].PM, TRANGTHAI, TBUOC[k - 1].P); MATADD(ref TBUOC[k].PM, TBUOC[k].PM, PROC_ BATDINH); + Xác định gia số Kalman theo công thức (2.36) Cập nhật gia số Kalman, K = (PT*H')(H*PT*H' + R)^(-1) + Xác ma trận H theo công thức (2.25) MATTRA(ref HT, QUANSAT); + Xác định tính tốn lại theo bước thời gian k, gia số Kalman: K = PT*HT*(H*PT*HT + R)^(-1) TBUOC[k].KG[i, 0] = PMHT[i, 0] / HPMHTR[0, 0]; + Cập nhật xác định giá trị hậu nghiệm X P TBUOC[k].X[0, i] = TBUOC[k].XM[0, i] + TBUOC[k].KG[i,0]*TBUOC[k].RES[0, 0]; P = [I-K*H]*PT MATMUL(ref KH, TBUOC[k].KG, QUANSAT); MATSUB(ref IKH, I44, KH); 155 MATMUL(ref TBUOC[k].P, IKH, TBUOC[k].PM); + Tính tốn phần sai số (Phần dư) MATMUL(ref HXM, QUANSAT, TBUOC[k].XM); + Xác định vector Zk “phần dư” theo công thức (2.24) TBUOC[k].Z[0, 0] = TV[k]; TBUOC[k].PHANDU[0, 0] = TBUOC[k].Z[0, 0] - HXM[0, 0]; + Cập nhật xác định giá trị hậu nghiệm theo công thức (2.35) TV_LAN = TBUOC[k].X[0, 0]; PT_LAN = RT[k] + TV_LAN; + ' Xác định hiệp phương sai sai số sau theo công thức (2.28): P = [I-K*H]*PT Call MATMUL(KH, KG, QUANSAT) + Xác định tính tốn TĐTX theo công thức (2.45), (2.46), (2.47) If PT_LAN >= TDKH Then TĐTX = RT[k - 1] + (TDKH - PT_LAN_CU) / (PT_LAN – PT_LAN_CU); StDev = Math.Sqrt(TBUOC[k].P[0,0]); PR = NORMDIST(TDKH, TĐTX, StDev, true); PT_LAN_CU = Math.Round(PT_LAN,14); TĐTX = RT[k - 1] + (TDKH - PT_TRUOC) / (PT_SAU - PT_TRUOC); + Xác định tính tốn xác suất thời điểm hồn thành theo (2.40a), (2.40b) PR = NORMDIST(TDKH, TDKX, StDev, true); + Xác định Trị trung bình độ lệch chuẩn theo công thức (E(tij) = (tm + 4tei + tp)/6 ; Độ lệch chuẩn: σij = (tp - tm)/6; Phương sai: : σij2 = (tp - tm)2/36) Câu lệnh C#: StDev = Math.Sqrt(TBUOC[k].P[0,0]); 156 TDKH_SD = (TDKH_P - TDKH_O) / 6; TDKH_TBINH = (TDKH_O + * (double)TDKH + TDKH_P) / 6; + Xác định tính tốn lại độ lệch chuẩn double fact = standard_dev * Math.Sqrt(2.0 * Math.PI); double expo = (x - mean) * (x - mean) / (2.0 * standard_dev * standard_dev); return Math.Exp(-expo) / fact; pi = 3.14159265358979 PRs(ITDDB) = Math.Round(normdist(TDKH, TĐTX, Sqr(.) * 100, + Xác định phân phối xác suất trước sau theo công thức (3.1) (3.2) normdist = normdistx(x, mean, standard_dev) dr[1] = Math.Round(NORMDIST(kXS, TDKH, TDKH_SD, true), 5); dr[2] = Math.Round(NORMDIST(kXS, Maxst, true), 5); 157 Phụ lục Giới thiệu tổng quan dự án áp dụng A Dự án thứ nhất: Khu học xá quốc tế Đà Nẵng Chủ đầu tư: UBND thành phố Đà Nẵng Đơn vị quản lý, điều hành dự án: Ban sở Hạ tầng ưu tiên Địa điểm xây dựng dự án: Phường Hòa Hải, quận Ngũ Hành Sơn, thành phố Đà Nẵng; Quy mô dự án : - Xây dựng 02 tòa nhà ký túc xá tầng, 02 tòa nhà tầng, khối nhà ăn, bệnh xá, nhà để xe, nhà làm việc ban quản lý - Tiến độ thi công theo Hồ sơ dự thầu: 25 tháng Tổng vốn đầu tư : 180 tỷ đồng B Dự án thứ hai: Đường vành đai Sơn Trà Chủ đầu tư: Sở giao thông vận tải thành phố Đà Nẵng Đơn vị quản lý, điều hành dự án: Ban cơng trình Trọng điểm Địa điểm xây dựng dự án: Phường Thọ Quang, quận Sơn Trà, thành phố Đà Nẵng; Quy mô dự án : - Xây dựng đường giao thơng với tổng chiều dài: L = 7.182,00m, có mặt cắt ngang Bn = 22,5m(3 + 10,5 + 3,0)m; độ dốc ngang 2%, mái - Xây dựng hệ thống Bó vỉa- vỉa hè: Bó vỉa, kết cấu thân bê tông M250 đá (1x2)cm, đúc sẵn lắp ghép Phần chân đế bê tông M200 đá (1x2)cm đổ chỗ, lớp đệm đá dăm dày 10cm Chiều cao bó vỉa 15cm, 158 chiều rộng rãnh biên 20-30cm, mặt vát 2/1; Vỉa hè Lót gạch Block M200 màu ghi dày 6cm, lớp cát đầm chặt K = 0, 90 dày 10cm - Hệ thống thoát nước: Tổng chiều dài L = 8.233,0m Sử dụng mương bê tông ly tâm (BTLT) vỉa hè, độ D = 60cm - Tiến độ thi công theo Hồ sơ dự thầu: 28 tháng Tổng vốn đầu tư : 120 tỷ đồng C Dự án thứ ba: Cầu qua Vịnh Đà Nẵng Chủ đầu tư: UBND thành phố Đà Nẵng Đơn vị quản lý, điều hành dự án: Ban cơng trình Trọng điểm Địa điểm xây dựng: Phường Nại Hiên Đông, quận Sơn Trà, Đà Nẵng; Quy mô dự án: Xây dựng cầu qua vịnh Đà Nẵng, tổng chiều dài toàn cầu 1.635 mét, rộng 15 mét, ba tháp cao 40 mét, khoảng cách ba trụ tháp 205 mét, độ tĩnh không thông thuyền 20 mét Tổng vốn đầu tư : 700 tỷ đồng D Dự án thứ tư: Trung Tâm giáo dục Thường xuyên Hòa Vang Chủ đầu tư: Sở Giáo dục thành phố Đà Nẵng Đơn vị quản lý, điều hành dự án: UBND Huyện Hòa Vang Địa điểm xây dựng: Huyện Hịa Vang, thành phố Đà Nẵng; Quy mơ dự án: Khối phòng học nhà tầng, khối hiệu bộ, thí nghiệm 02 nhà tầng, nhà để xe, hệ thống thử nghiệm; - Tiến độ thi công theo Hồ sơ dự thầu: 7.7 tháng Tổng vốn đầu tư : 70 tỷ đồng 159 ... 1: Một số vấn đề lý luận xác định tiến độ thi cơng cơng trình có tính đến yếu tố bất định Chương 2: Xác định tiến độ thi cơng cơng trình có tính đến yếu tố bất định theo phương pháp dự báo xác. .. Chương XÁC ĐỊNH TIẾN ĐỘ THI CƠNG CƠNG TRÌNH CĨ TÍNH ĐẾN YẾU TỐ BẤT ĐỊNH THEO PHƯƠNG PHÁP DỰ BÁO XÁC SUẤT 54 2.1 Yếu tố bất định tác động đến thực tiến độ thi công 54 2.1.1 Thế yếu tố bất định. .. bất định lĩnh vực thi công cơng trình Qua đó, đánh giá việc tác động yếu tố đến trình thực tiến độ - Đề xuất phương pháp xác định tiến độ thi cơng có tính đến tác động yếu tố bất định Phương pháp