Ứng dụng lý thuyết mờ để đánh giá rủi ro trong công trình ngầm

ĐẠI HỌC QUỐC GIA THÀNH PHỐ HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA NGUYỄN KIM CƯỜNG ỨNG DỤNG LÝ THUYẾT MỜ ĐỂ ĐÁNH GIÁ RỦI RO TRONG CƠNG TRÌNH NGẦM Chuyên ngành: XÂY DỰNG CẦU HẦM LUẬN VĂN THẠC SĨ TP HỒ CHÍ MINH THÁNG 06 NĂM 2012 CƠNG TRÌNH ĐƯỢC HỒ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 : PGS.TS BÙI CÔNG THÀNH Nhận xét ……………………………………………………………………………………… ……………………………………………………………………………………… ……………………………………………………………………………………… ……………………………………………………………………………………… ……………………………………………………………………………………… ……………………………………………………………………………………… ……………………………………………………………………………………… ………………………………………………………………… Cán chấm nhận xét :……………………………………………………………………………………… ……………………………………………………………………………………… ……………………………………………………………………………………… ……………………………………………………………………………………… ……………………………………………………………………………………… ……………………………………………………………………………………… ……………………………………………………………………………………… ………………………………………………………………… Cán chấm nhận xét :……………………………………………………………………………………… ……………………………………………………………………………………… ……………………………………………………………………………………… ……………………………………………………………………………………… ……………………………………………………………………………………… ……………………………………………………………………………………… …………………………………………………………………… Luận văn thạc sĩ bảo vệ Trường Đại học Bách Khoa, ĐHQG TP HCM Ngày tháng năm ĐẠI HỌC QUỐC GIA TP.HC CỘNG HOÀ XÃ HỘI CHỦ NGHIÃ VIỆT NAM TRƯỜNG ĐẠI HỌC BÁCH KHOA Độc Lập - Tự Do - Hạnh Phúc -oOo Tp HCM, ngày tháng năm NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên: NGUYỄN KIM CƯỜNG Ngày, tháng, năm sinh: Yên Phái: Nam 31-10-1982 Chuyên ngành: Xây dựng Cầu Hầm Nơi sinh: Phú MSHV: 03808504 1-TÊN ĐỀ TÀI: ỨNG DỤNG LÝ THUYẾT MỜ ĐỂ ĐÁNH GIÁ RỦI RO CƠNG TRÌNH NGẦM 2- NHIỆM VỤ VÀ NỘI DUNG: Tìm hiểu cần thiết nghiên cứu rủi ro cơng trình ngầm, cố cơng trình ngầm xảy q khứ Tìm hiểu tình hình nghiên cứu rủi ro ứng dụng lý thuyết mờ việc đánh giá rủi ro nước giới Tìm hiểu sở lý thuyết mờ cách nhận dạng, đánh giá rủi ro cơng trình ngầm Xây dựng chương trình tính tốn tính hấp dẫn dự án ngơn ngữ lập trình MatLab sau so sánh kết tính tốn rủi ro tác giả khác 3- NGÀY GIAO NHIỆM VỤ: 17/01/2011 4- NGÀY HOÀN THÀNH NHIỆM VỤ: 16/04/2012 5- HỌ VÀ TÊN CÁN BỘ HƯỚNG DẪN: PGS.TS BÙI CÔNG THÀNH CÁN BỘ HƯỚNG DẪN CHỦ NHIỆM BỘ MÔN KHOA QL CHUYÊN NGÀNH (Họ tên chữ ký) (Họ tên chữ ký) (Họ tên chữ ký) PGS.TS BÙI CÔNG THÀNH i LỜI CẢM ƠN Luận văn “ỨNG DỤNG LÝ THUYẾT MỜ ĐỂ ĐÁNH GIÁ RỦI RO TRONG CÔNG TRÌNH NGẦM” thực từ ngày 01/10/2011 đến ngày 16/04/2012 Mục đích đề tài đánh giá rủi ro Cơng Trình Ngầm, từ đánh giá tính hấp dẫn dự án xây dựng Cơng Trình Ngầm giúp chủ đầu tư cân nhắc đầu tư Em xin chân thành cảm ơn PGS.TS Bùi Cơng Thành có định hướng, dẫn tận tình để hồn thành luận văn Cảm ơn Công Ty TNHH Apave Việt Nam Đông Nam Á, nơi em cộng tác thực tập, nơi em thu tập nhiều thông tin từ dự án nước quốc tế, tiếp nhận nhiều kiến thức quí giá từ đồng nghiệp lãnh đạo Cảm ơn Hội Địa Kỹ Thuật Việt Nam, thông qua trao đổi thông tin liên lạc thành viên hội, em có thêm thơng tin cần thiết để viết thuật tốn chương trình Em xin chân thành cảm ơn giảng viên môn cầu đường khoa xây dựng, cảm ơn phòng đào tạo sau đại học Trường Đại Học Bách Khoa thành phố Hồ Chí Minh Xin cảm ơn trường cao đẳng xây dựng số (nay trường Đại Học Xây Dựng Miền Trung ) đồng nghiệp môn Cầu Đường tạo điều kiện thời gian học Cảm ơn giảng viên phản biện có nhận xét chân thành để luận văn hoàn thiện Vì thời gian làm luận văn có hạn hạn chế kiến thức thân, luận văn không tránh khỏi thiếu sót, mong bạn đồng nghiệp, q thầy đóng góp ý kiến để luận văn hồn chỉnh Tác giả luận văn Nguyễn Kim Cường ii TĨM TẮT Các cơng trình ngầm xây dựng phát sinh nhiều rủi ro chủ quan khách quan Một số rủi ro liệt kê sau: không chắn thông tin, số liệu thiếu sót, phối hợp thiếu chặt chẽ đơn vị khảo sát thiết kế, nhà thầu,…Các yếu tố rủi ro ảnh hưởng đáng kể đến việc hoàn thành dự án Do vậy, việc đánh giá cách chi tiết rủi ro quan trọng để có kế hoạch xử lý kiểm sốt phù hợp Luận văn ứng dụng lý thuyết mờ để phân tích rủi ro cơng trình ngầm Lý thuyết mờ đo lường độ xác phần tử tập hợp hàm liên thuộc, rủi ro nhận dạng chọn lọc qua nghiên cứu nhiều tác giả Số liệu đầu vào rủi ro thu thập thông qua phiếu đánh giá Một chương trình máy dựa sở lý thuyết mờ số liệu thu thập để đánh giá tính hấp dẫn cho dự án cơng trình ngầm viết ngơn ngữ lập trình MATLAB Thí dụ số áp dụng để đánh giá tính hấp dẫn Hầm Hải Vân số liệu khách quan thu thập thực có so sánh với nghiên cứu nhiều tác giả khác Kết số cho thấy lý thuyết mờ ứng dụng để đánh giá rủi ro cho dự án cơng trình ngầm iii ABSTRACT The construction of underground works as may arise quite risky both subjective and objective Some risks are listed as follows: uncertain information, data deficiencies, lack of close coordination between the survey and design units, contractors, these risk factors influence include the completion of the project Therefore, a detailed assessment of risks is important to plan and control process more consistent This paper applied fuzzy theory to analyze risks in underground work Fuzzy theory of measurement accuracy of the elements in the collection using a continuous function of, the risks are identified and selected through the study of many authors Data input of risk collected through the evaluation A computer program based on fuzzy theory and the data collected to assess the attractiveness for underground projects written in MATLAB programming language For example, some are applied to assess the attractiveness of the Hai Van tunnel and objective data gathered is done the research compared with many other authors The results indicate that fuzzy theory can be applied to risk assessment for an underground project iv MỤC LỤC LỜI CẢM ƠN……………………………………………………… ……………………i TĨM TẮT……………………………………………………………………………… ii DANH MỤC CÁC HÌNH VẼ…………………………………………………………Vii DANH MỤC CÁC BẢNG BIỂU………………………………….……………………ix CHƯƠNG MỞ ĐẦU……………….…………………… ………………………… 1.1 ĐẶT VẤN ĐỀ 1.2 MỤC TIÊU NGHIÊN CỨU ĐỀ TÀI 1.3 PHẠM VI VÀ ĐỐI TƯỢNG NGHIÊN CỨU 1.4 Ý NGHĨA KHOA HỌC ĐỀ TÀI CHƯƠNG TỔNG QUAN 2.1 GIỚI THIỆU 2.2 TÌNH HÌNH NGHIÊN CỨU TRONG NƯỚC 2.3 TÌNH HÌNH NGHIÊN CỨU TRÊN THẾ GIỚI 2.4 NHỮNG RỦI RO THƯỜNG GẶP TRONG CƠNG TRÌNH NGẦM……………7 2.5 KẾT LUẬN CHƯƠNG 11 CHƯƠNG CƠ SỞ LÝ THUYẾT 12 3.1 GIỚI THIỆU 12 3.2 CƠ SỞ LÝ THUYẾT MỜ 12 3.2.1 Khái niệm lý thuyết mờ 13 3.2.2 Các phép toán tập mờ 14 3.2.2.1 Phép hợp 14 3.2.2.2 Phép giao 14 3.2.2.3 Phép bù tập mờ 15 v 3.2.2.4 Quan hệ mờ 15 3.2.2.5 Hợp thành mờ 16 3.2.2.6 Số mờ 17 3.2.3 Biến ngôn ngữ 17 3.2.4 Logic mờ 17 3.2.4.1 Logic 18 3.2.4.2 Khái niệm Logic mờ 18 3.2.5 Suy diễn 18 3.2.5.1 Khái niệm 18 3.2.5.2 Phép suy diễn 18 3.2.5.3 Mệnh đề mờ 19 3.2.5.4 Hàm kéo theo mờ 21 3.2.5.5 Mệnh đề điều kiện 21 3.2.5.6 Suy diễn mờ 22 3.2.6 Khử mờ 25 3.2.6.1 Các phương pháp khử mờ phổ biến 25 3.2.6.2 Tiêu chuẩn chọn lựa phương pháp khử mờ 26 3.2.7 Giải toán hệ thống suy luận mờ 26 3.3 NHẬN DẠNG RỦI RO TRONG CƠNG TRÌNH NGẦM 38 3.3.1 Khái niệm rủi ro 39 3.3.2 Nguồn gốc rủi ro 40 3.3.3 Nhận dạng rủi ro cơng trình ngầm 41 CHƯƠNG CHƯƠNG TRÌNH ĐÁNH GIÁ TÍNH HẤP DẪN DỰ ÁN………… 51 4.1 SỐ LIỆU ĐẦU VÀO CHƯƠNG TRÌNH 52 4.2 ĐỊNH NGHĨA MỜ CHO NHẬN XÉT BẰNG NGÔN NGỮ 52 4.3 VIẾT THUẬT TỐN CHƯƠNG TRÌNH 53 4.3.1 Thuật tốn chương trình hàm hình thang nửa phải 54 4.3.2 Thuật tốn chương trình hàm hình thang nửa trái …55 vi 4.3.3 Thuật toán chương trình hàm hình tam giác 56 4.3.4 Thuật tốn chương trình đánh giá trễ tiến độ dự án 56 4.3.5 Thuật tốn chương trình đánh giá mức độ tác động đến dự án 58 4.3.6 Thuật tốn chương trình đánh giá tính hấp dẫn dự án 59 4.4 KIỂM TRA ĐỘ TIN CẬY CỦA CHƯƠNG TRÌNH 60 4.4.1 Thu thập số liệu cho biến số đầu vào phiếu thăm dò ý kiến chuyên gia… 60 4.4.2 So sánh kết với nghiên cứu Patrick.X.W.Zou [18] cộng A.Deviprasadh [23] 64 4.4.3 Phân tích kết nghiên cứu Shih-Tong Lu cộng [17] Sid-Ghod Jakkapan Jintanapakanont [5] với José Ramón Terés Flores [22] 67 4.5 ỨNG DỤNG CHƯƠNG TRÌNH NGƯỜI VIẾT ĐÁNH GIÁ TÍNH HẤP DẪN ĐƯỜNG HẦM ĐÈO HẢI VÂN 71 CHƯƠNG KẾT LUẬN 77 5.1 KẾT LUẬN 77 5.2 HƯỚNG NGHIÊN CỨU………………………….…… ……… ………… ….77 TÀI LIỆU THAM KHẢO .78 PHỤ LỤC 81 vii DANH MỤC CÁC HÌNH VẼ Hình 1.1 Xe buýt lọt hố tử thần đường Hầm Đức………………………… ………… Hình 1.2 Sập Hầm ảnh hưởng đến cơng trình xung quanh lún sụp………………… Hình 1.3 Sập tường chắn làm xuất hố sâu 30m đại lộ Nicoll kế bên khu vực làm Hầm………………………………………………………………………… ……….… …………2 Hình 2.1 Sạt lỡ gương đào……………………………………………………………………….10 Hình 3.1 Tập mờ tập rõ 13 Hình 3.2 dạng hàm thành viên Tập mờ 14 Hình 3.3 Đồ thị thể phép hợp hai tập mờ 14 Hình 3.4 Đồ thị thể phép giao hai tập mờ 15 Hình 3.5 Đồ thị quan hệ biến ngôn ngữ thuộc tính nhiệt độ 18 Hình 3.6 Đồ thị ví dụ thuộc tính cao nhiệt độ 19 Hình 3.7 Đồ thị hàm thành viên biến đinh tính 20 Hình 3.8 Đồ thị hàm liên thuộc miêu tả tính trẻ Tú 21 Hình 3.9 Sơ đồ khối mơ hình suy luận mờ 27 Hình 3.10 Đồ thị hàm liên thuộc biến đầu vào 29 Hình 3.11 Đồ thị hàm liên thuộc biến đầu vào mức độ tác động đến dự án 29 Hình 3.12 Đồ thị hàm liên thuộc biến hấp dẫn dự án 30 Hình 3.13 Hàm liên thuộc biến mờ trễ tiến độ dự án 32 Hình 3.15 Giao số thực với hàm thành viên tập mờ 34 Hình 3.16 Thể giá trị biến kích biến số đầu vào 35 Hình 3.17 Sơ đồ kết hợp qui tắc cho tập hợp biến mờ đầu 36 85 24 IF PD VH AND PI L THEN PA M 25 IF PD VH AND PI VL THEN PA M Phụ Lục D: Code chương trình đánh giá tính xác suất trễ tiến độ dự án % Chuong trinh tinh xac xuat rui ro tre tien du an cong trinh Ngam clear all; clc; x=(0:1:100); VL=Thangphai(x,12,23); L=Tamgiac(x,18,36,42); M=Tamgiac(x,39,50,63); H=Tamgiac(x,60,72,84); VH=Thangtrai(x,79,100); disp(' moi qui vi nhap xac suat xay cua cac bien so sau'); disp(' Nguyên nhân dan den rui ro quan ly du án'); EF=input('1:rui ro kinh te tai chinh nha thau EF= '); CP=input('2:rui ro hop dong phap ly CP='); SS=input('3:Rui ro lien quan den an toan xa hoiSS='); %Nhap cac bien co ban dan den rui ro khao sat dac GW=input('4:Hieu sai dia chat nuoc ngam GW='); CS=input('5:Nang luc cua doi khao sat CS='); CC=input('6:Sai lam tinh toan chep CC='); %Nhap cac bien co ban dan dan den rui ro thiet ke LI=input('7:Thieu hang muc cong viec LI='); CD=input('8:Do thiet ke thay doiCD= '); DD=input('9:Mau thuan tai lieu thiet ke DD= '); CT=input('10:Sai lam cong nghe, tinh toan phuc tap CT='); %Nhap cac bien co ban dan dan den rui ro thi cong lap dung QW=input('11:Chat luong Vat tu, Nhan cong, Thiet bi QW='); UNG=input('12:Hien tuong bat thuong dia chat khong xu ly kipUNG= '); RP=input('13:Ke hoach kiem soat rui ro kem RP= '); CM=input('14:Bien phap thi cong khong phu hop CM='); % Nhap cac bien co ban dan den rui ro bat kha khang RJ=input('15:Rui ro bat kha khangRJ='); %Tre tien ben thu ba ID=input('16:rui ro tre tien cac ben lien quan khacID= '); %Xac dinh ham lien thuoc cua bien rui ro quan ly du an RMVL=max([min([L(EF) L(CP) L(SS)]) VL(EF) VL(CP) VL(SS)]); RML=max([L(EF) L(CP) L(SS)]); RMM=max([M(EF) M(CP) M(SS)]); RMH=max([min([M(EF) M(CP) M(SS)]) H(EF) H(CP) H(SS)]); RMVH=max([min([H(EF) H(CP) H(SS)]) VH(EF) VH(CP) VH(SS)]); %Xac dinh ham lien thuoc cua bien rui ro thi cong lap dung RCVL=min([ VL(QW) VL(UNG) VL(RP) VL(CM)]); RCL=max([ L(QW) L(UNG) L(RP) L(CM)]); RCM=max([M(QW) M(UNG) M(RP) M(CM)]); RCH=max([min([M(QW) M(UNG) M(RP) M(CM)]) H(QW) H(UNG) H(RP) H(CM)]); 86 RCVH=max([min([H(QW) H(UNG) H(RP) H(CM)]) VH(QW) VH(UNG) VH(RP) VH(CM)]); %Xac dinh ham lien thuoc cua bien rui ro khao sat RSVL=max([VL(GW) VL(CS) VL(CC)]); RSL=max([L(GW) L(CS) L(CC)]); RSM=max([M(GW) M(CS) M(CC)]); RSH=max([min([M(GW) M(CS) M(CC)]) H(GW) H(CS) H(CC)]); RSVH=max([min([H(GW) H(CS) H(CC)]) VH(GW) VH(CS) VH(CC)]); %Xac dinh ham lien thuoc cua bien rui ro thiet ke RDVL=max([ VL(LI) VL(CD) VL(DD) VL(CT)]); RDL=max([L(LI) L(CD) L(DD) L(CT)]); RDM=max([M(LI) M(CD) M(DD) M(CT)]); RDH=max([min([M(LI) M(CD) M(DD) M(CT)]) H(LI) H(CD) H(DD) H(CT)]); RDVH=max([min([H(LI) H(CD) H(DD) H(CT)]) VH(LI) VH(CD) VH(DD) VH(CT)]); %Xac dinh cac tham so lien thuoc bien dau PDelayVL=max([min([RDVL RSVL]) VL(RJ) VL(ID) RMVL]); PDelayL=max([min([RCL RDL RSL]) L(RJ) L(ID) RML]); PDelayM=max([RCM RDM RSM L(RJ) RMM]); PDelayH=max([min([RCM RDM RSM]) M(RJ) RMH]); PDelayVH=max([min([RCH RDH RSH]) H(RJ) RMVH]); PDelay1=min(VL,PDelayVL); PDelay2=min(L,PDelayL); PDelay3=min(M,PDelayM); PDelay4=min(H,PDelayH); PDelay5=min(VH,PDelayVH); a=max(PDelay1,PDelay2); b=max(a,PDelay3); c=max(b,PDelay4); finished=max(c,PDelay5); figure u=(0:1:100); subplot(611), plot(u,PDelay1); subplot(612), plot(u,PDelay2); subplot(613), plot(u,PDelay3); subplot(614), plot(u,PDelay4); subplot(615), plot(u,PDelay5); subplot(616), plot(u,finished); centroid=(sum(x.*finished))/sum(finished); disp('gia tri xac suat rui ro tre tien la:'); disp(centroid); Phụ lục E: Đánh giá mức độ tác động đến dự án: % Chuong trinh tinh muc tac dong den du an cong trinh Ngam tre tien % clear all; clc; x=(0:1:100); VL=Thangphai(x,12,24); L=Tamgiac(x,19,28,40); M=Tamgiac(x,37,52,63); H=Tamgiac(x,58,69,78); VH=Thangtrai(x,76,93); disp(' Nhap cac bien so dau vao, tac dong den du an tre tien do'); 87 disp(' tac dong den du an cac yeu to cua quan ly'); EF=input('1:rui ro kinh te tai chinh nha thau EF= '); CP=input('2:rui ro hop dong phap ly CP='); SS=input('3:Rui ro lien quan den an toan xa hoiSS='); %Nhap cac bien co ban dan den tac dong den du an khao sat dac GW=input('4:Hieu sai dia chat nuoc ngam GW='); CS=input('5:Nang luc cua doi khao sat CS='); CC=input('6:Sai lam tinh toan chep CC='); %Nhap cac bien co ban dan den tac dong den du an rui ro thiet ke LI=input('7:Thieu hang muc cong viec LI='); CD=input('8:Do thiet ke thay doiCD= '); DD=input('9:Mau thuan tai lieu thiet ke DD= '); CT=input('10:Sai lam cong nghe, tinh toan phuc tap CT='); %Nhap cac bien co ban dan den tac dong den du an rui ro thi cong lap dung QW=input('11:Chat luong Vat tu, Nhan cong, Thiet bi QW='); UNG=input('12:Hien tuong bat thuong dia chat khong xu ly kipUNG= '); RP=input('13:Ke hoach kiem soat rui ro kem RP= '); CM=input('14:Bien phap thi cong khong phu hop CM='); % Nhap bien co ban dan tac dong den du an rui ro bat kha khang RJ=input('15:Rui ro bat kha khang RJ= '); %Tac dong den du an ben thu ba tre ID=input('16:rui ro tre tien cac ben lien quan khacID= '); %Xac dinh ham lien thuoc cua bien rui ro quan ly du an RMVL=max([VL(EF) VL(CP) VL(SS)]); RML=max([L(EF) M(CP) L(SS)]); RMM=max([M(EF) M(SS)]); RMH=max([min([M(EF) M(CP) M(SS)]) H(EF) H(SS)]); RMVH=max([min([H(EF) H(CP) H(SS)]) VH(CP) VH(EF)]); %Xac dinh ham lien thuoc cua bien rui ro thi cong lap dung RCVL=max([ VL(QW) VL(UNG) L(CM)]); RCL=max([ L(QW) VL(UNG) L(RP) L(CM)]); RCM=max([min([L(QW) L(UNG) M(CM)]) M(QW) M(UNG)]) ; RCH=max([min([M(QW) M(UNG) M(CM)]) H(QW) H(UNG) H(CM) H(RP)]); RCVH=max([min([H(QW) H(UNG) H(RP) H(CM)]) VH(QW) VH(UNG) VH(RP) VH(CM)]); %Xac dinh ham lien thuoc cua bien rui ro khao sat RSVL=max([VL(GW) VL(CS) L(CC)]); RSL=max([L(GW) L(CS) M(CC)]); RSM=max([M(GW) M(CS)]); RSH=max([min([M(GW) M(CS) M(CC)]) H(GW) H(CS) ]); RSVH=max([min([H(GW) H(CS) H(CC)]) VH(GW) VH(CS)]); %Xac dinh ham lien thuoc cua bien rui ro thiet ke RDVL=min([ VL(LI) VL(CD) VL(DD) VL(CT)]); RDL=min([L(LI) L(CD) L(DD) L(CT)]); RDM=max([M(LI) M(CD) M(DD) M(CT)]); RDH=max([min([M(LI) M(CD) M(DD) M(CT)]) H(LI) H(CD) H(DD) H(CT)]); RDVH=max([min([H(LI) H(CD) H(DD) H(CT)]) VH(LI) VH(CD) VH(DD)]); %Xac dinh cac tham so lien thuoc bien dau PImpactVL=max([min([RDVL RSVL RCVL]) VL(RJ) VL(ID) RMVL]); PImpactL=max([min([RCL RDL RSL VL(RJ)]) L(ID) RML]); PImpactM=max([RCM RDM L(RJ) RMM]); PImpactH=max([min([RCM RDM RSM]) M(RJ) RMH]); PImpactVH=max([min([RCH RDH RSH]) H(RJ) RMVH]); PDelay1=min(VL,PImpactVL); 88 PDelay2=min(L,PImpactL); PDelay3=min(M,PImpactM); PDelay4=min(H,PImpactH); PDelay5=min(VH,PImpactVH); a=max(PDelay1,PDelay2); b=max(a,PDelay3); c=max(b,PDelay4); finished=max(c,PDelay5); figure u=(0:1:100); subplot(611), plot(u,PDelay1); subplot(612), plot(u,PDelay2); subplot(613), plot(u,PDelay3); subplot(614), plot(u,PDelay4); subplot(615), plot(u,PDelay5); subplot(616), plot(u,finished); centroid=(sum(x.*finished))/sum(finished); disp('gia tri muc tac dong den du an :'); disp(centroid); PHỤ LỤC F: %Chương trình đánh giá trễ tính hấp dẫn dự án clear all; clc; x=(0:1:100); VL=Thangphai(x,10,22); L=Tamgiac(x,20,31,43); M=Tamgiac(x,38,49,59); H=Tamgiac(x,57,68,81); VH=Thangtrai(x,78,93); Pdelayprob=input('nhap vao xac suat tre tien du an:'); Pdelayprob=round(Pdelayprob); Pdelayimpact=input('nhap vao tac dong tre tien du an:'); Pdelayimpact=round(Pdelayimpact); c1=max(VL(Pdelayprob),VH(Pdelayimpact)); c2=max(VL(Pdelayprob),H(Pdelayimpact)); c3=max(VL(Pdelayprob),M(Pdelayimpact)); c4=max(VL(Pdelayprob),L(Pdelayimpact)); c5=max(VL(Pdelayprob),VL(Pdelayimpact)); c6=max(L(Pdelayprob),VH(Pdelayimpact)); c7=max(L(Pdelayprob),H(Pdelayimpact)); c8=max(L(Pdelayprob),M(Pdelayimpact)); c9=max(L(Pdelayprob),L(Pdelayimpact)); c10=max(L(Pdelayprob),VL(Pdelayimpact)); c11=max(M(Pdelayprob),VH(Pdelayimpact)); c12=max(M(Pdelayprob),H(Pdelayimpact)); c13=max(M(Pdelayprob),M(Pdelayimpact)); c14=max(M(Pdelayprob),L(Pdelayimpact)); c15=max(M(Pdelayprob),VL(Pdelayimpact)); c16=max(H(Pdelayprob),VH(Pdelayimpact)); c17=max(H(Pdelayprob),H(Pdelayimpact)); c18=max(H(Pdelayprob),M(Pdelayimpact)); c19=max(H(Pdelayprob),L(Pdelayimpact)); 89 c20=max(H(Pdelayprob),VL(Pdelayimpact)); c21=max(VH(Pdelayprob),VH(Pdelayimpact)); c22=max(VH(Pdelayprob),H(Pdelayimpact)); c23=max(VH(Pdelayprob),M(Pdelayimpact)); c24=max(VH(Pdelayprob),L(Pdelayimpact)); c25=max(VH(Pdelayprob),VL(Pdelayimpact)); PA1=min(H,c1); PA2=min(H,c2); PA3=min(VH,c3); PA4=min(VH,c4); PA5=min(VH,c5); PA6=min(H,c6); PA7=min(H,c7); PA8=min(VH,c8); PA9=min(VH,c9); PA10=min(H,c10); PA11=min(L,c11); PA12=min(M,c12); PA13=min(M,c13); PA14=min(H,c14); PA15=min(VH,c15); PA16=min(VL,c16); PA17=min(L,c17); PA18=min(L,c18); PA19=min(M,c19); PA20=min(H,c20); PA21=min(VL,c21); PA22=min(VL,c22); PA23=min(L,c23); PA24=min(M,c24); PA25=min(H,c25); a=max(PA1,PA2); b=max(a,PA3); c=max(b,PA4); d=max(c,PA5); e=max(d,PA6); f=max(e,PA7); g=max(f,PA8); h=max(g,PA9); i=max(h,PA10); j=max(i,PA11); k=max(j,PA12); l=max(k,PA13); m=max(l,PA14); n=max(m,PA15); o=max(n,PA16); p=max(o,PA17); q=max(p,PA18); r=max(q,PA19); s=max(r,PA20); t=max(s,PA21); u=max(t,PA22); v=max(u,PA23); 90 w=max(v,PA24); finished=max(w,PA25); u=(0:1:100); plot(u,finished); %gia tri khu mo dau bang phuong phap tri trung tam central=(sum(x.*finished))/(sum(finished)); disp('Gia tri dinh luong tinh hap dan du an la:'); disp(central); 91 PHỤ LỤC J PHIẾU THĂM DÒ Ý KIẾN CHUYÊN GIA Kính thưa quí vị! Để phục vụ cho việc nghiên cứu, xin q vị vài phút điền giúp bảng thơng tin đánh giá sau: Tôi xin cam kết, thông tin q vị đóng góp bảo mật, cảm ơn q vị bớt thời gian giúp tơi hồn thành chương trình nghiên cứu “Đánh giá tính hấp dẫn dự án xây dựng cơng trình Ngầm” Bảng 4.6 Đánh giá xác suất xảy trễ tiến độ dự án xây dựng Cơng Trình Ngầm (đơn vị: %) Rất ST TÊ T N 10 11 12 13 14 15 16 CÁC YẾU TỐ RỦI RO EF Kinh tế tài CP Hợp đồng pháp lý SS Rủi ro liên quan đến an toàn xã hội G Hiểu sai địa chất nước Ngầm W CS Năng lực đội khảo sát CC Sai lầm tính tốn chép LI Thiếu hạng mục công việc CD Do thay đổi thiết kế DD Mâu thuẫn tài liệu thiết kế CT Sai lầm cơng nghệ, tính tốn phức Q tạp Chất lượng Vật tư, Nhân Công, Thiết W UN bị Không xử lý kịp cố bất thường G RP địa hoạch chất phòng ngừa rủi ro Kế CM Biện pháp thi công không phù hợp RJ Bất khả kháng ID Do bên liên quan khác thấp Trung Thấp (0-23) (18-42) bình (39-63) Cao (6084) Rất Cao (79-100) 92 Đánh giá mức độ tác động đến dự án xây dựng Cơng Trình Ngầm (đơn vị: %) S T TÊ N Rất nhỏ CÁC YẾU TỐ RỦI RO (0-24) Trung Nhỏ (19-40) bình (37-63) Lớn Rất lớn (58-78) (76-100) T 1 1 1 EF Kinh tế tài CP Hợp đồng pháp lý SS Rủi ro liên quan đến an toàn xã hội G Hiểu sai địa chất nước Ngầm W CS Năng lực đội khảo sát CC Sai lầm tính tốn chép LI Thiếu hạng mục công việc CD Do thay đổi thiết kế DD Mâu thuẫn tài liệu thiết kế CT Sai lầm cơng nghệ, tính tốn phức Q tạp lượng Vật tư, Nhân Công, Thiết Chất W UN bị Không xử lý kịp cố bất thường G RP địa hoạch chất phòng ngừa rủi ro Kế CM Biện pháp thi công không phù hợp RJ Bất khả kháng ID Do bên liên quan khác Các yếu tố rủi ro khác: ……………………………………………………… Để gửi kết điều tra, chương trình đánh giá “Tính hấp dẫn dự án cơng trình ngầm” cho q vị, vui lịng cho biết thơng tin sau: (khơng bắt buộc) Họ tên: …………………………………………………………………… Số năm kinh nghiệm……………………………………………………… địa email: ……………………… , số điện thoại………………………… Xin chân thành cảm ơn đóng góp q vị! 93 LÝ LỊCH TRÍCH NGANG • Họ Và Tên: NGUYỄN KIM CƯỜNG • Ngày, tháng, năm sinh: 31-10-1982 Nơi sinh: Phú Yên • Địa liên lạc: F32, Khu tái định cư Trần Phú, Tuy Hồ, Phú n • Điện thoại: 0905824763 • Email: nkc311082@gmail.com; nguyenkimcuong@cuc.edu.vn Q TRÌNH ĐÀO TẠO • Từ năm 2000-2005: Học đại học trường đại học Giao Thông Vận Tải sở II Q9, TP.HCM • Từ năm 2008 đến nay: Học viên cao học trường đại học Bách Khoa TP HCM Q TRÌNH CƠNG TÁC • Từ tháng 05-2005 đến tháng 12-2005 làm việc công ty cổ phần đầu tư tư vấn thiết kế xây dựng giao thơng Phú n • Từ năm 2006 đến công tác trường Cao Đẳng Xây Dựng số (nay trường đại học Xây Dựng Miền Trung) • Từ năm 08/2008-12/2011 có hợp tác làm việc công ty TNHH Apave Việt Nam Đông Nam Á PHỤ LỤC G Bảng tính chi tiết số phiếu giá trị trung bình biến số thu từ kết điều tra dành cho PD Biến KH VL L M H VH VL L M H VH VL L M H VH VL L M H VH VL L M H VH VL L M H VH VL L M 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 20 37 31 40 30 49 65 60 26 40 61 70 85 11 62 35 30 31 32 30 42 62 40 40 70 37 38 39 40 41 42 43 44 45 46 47 48 21 42 70 65 60 44 49 70 22 40 65 60 82 85 22 17 25 35 45 20 10 20 39 50 60 65 70 40 17 70 47 50 40 57 52 70 50 72 35 70 49 40 65 32 42 65 71 50 50 67 50 64 35 39 42 47 52 45 47 45 48 70 27 19 41 60 41 47 37 57 50 40 55 39 35 33 45 36 53 65 10 36 65 58 20 43 39 15 42 33 70 21 40 55 39 57 65 71 10 35 18 40 42 57 15 63 ## 30 30 45 31 37 32 18 65 15 39 67 27 50 45 15 47 18 65 37 50 30 35 42 47 73 60 20 32 25 32 20 42 20 35 21 20 25 17 40 15 41 30 25 50 21 47 65 27 20 20 21 31 19 48 39 72 67 63 ## 21 33 49 69 67 65 73 20 20 45 60 53 65 10 40 33 60 47 32 21 49 48 50 65 20 35 55 65 35 32 49 67 95 85 30 50 20 23 36 25 35 42 ## 15 10 28 60 32 75 18 30 34 50 80 15 34 65 79 20 10 10 22 35 34 25 20 37 40 19 15 10 25 20 15 37 52 19 30 41 30 22 46 67 73 72 77 62 80 87 85 10 20 30 40 51 37 42 19 30 30 57 10 30 54 63 15 20 20 18 19 37 20 42 21 30 37 10 33 60 49 57 19 59 18 37 42 35 25 30 40 33 SP PD 17 31 11 45 14 66 89 18 10 31 22 46 67 0 13 12 34 23 49 69 0 16 20 27 12 46 10 67 100 15 17 32 49 67 88 14 14 15 31 50 70 86 13 20 28 12 51 10 11 12 13 14 H VH VL L M H VH VL L M H VH VL L M H VH VL L M H VH VL L M H VH VL L M H VH VL L M H 70 67 79 60 82 80 87 22 10 18 32 30 52 63 72 37 50 60 80 78 65 19 40 33 26 60 80 19 60 80 60 50 77 20 47 40 71 60 77 18 17 40 35 18 10 35 25 18 21 33 30 22 38 37 47 47 62 65 20 35 33 40 19 17 24 21 40 50 45 22 10 40 41 18 20 35 22 20 32 19 40 23 63 20 10 42 60 32 25 40 73 40 20 27 10 15 21 41 32 40 59 42 50 55 39 70 79 19 20 25 53 50 40 57 41 25 60 35 40 45 40 45 40 50 67 73 41 20 28 22 21 25 20 20 47 47 95 44 61 80 70 85 20 10 30 55 79 62 22 40 40 50 40 52 15 41 37 48 53 45 18 57 54 # 79 42 55 39 45 42 55 67 75 87 18 35 58 60 44 68 79 92 10 19 50 10 40 22 70 40 45 74 37 49 52 48 20 85 16 24 80 62 15 20 22 20 40 68 80 80 73 59 73 89 67 66 52 45 48 45 45 70 85 32 60 47 40 73 35 62 53 25 58 51 23 41 28 35 65 80 89 10 41 20 20 55 20 20 35 45 37 50 40 64 59 60 37 39 35 22 49 57 55 49 18 49 27 69 20 52 51 35 67 69 15 35 60 18 85 83 22 42 19 30 15 60 70 60 83 25 35 39 40 67 85 13 22 21 40 39 49 52 17 30 50 19 23 30 50 37 52 60 20 55 39 45 52 40 19 30 29 32 26 30 61 41 51 65 55 68 11 17 12 10 25 13 19 11 10 19 9 20 12 7 22 15 17 72 18 29 50 69 16 32 52 68 16 29 50 70 14 34 49 73 88 14 27 50 71 85 13 35 50 76 87 14 29 51 68 VH 79 90 90 79 VL 15 11 20 11 11 10 21 12 21 10 15 15 12 10 10 21 12 10 13 15 20 14 18 16 10 12 21 19 34 L 22 35 20 20 23 40 40 39 40 60 50 15 M H 84 60 70 VH VL 20 20 20 15 21 17 L 20 35 20 27 48 40 40 62 55 40 45 50 50 57 60 49 55 45 12 16 M H 69 72 61 60 84 76 84 60 75 82 65 80 72 71 67 62 79 67 18 VH 89 89 82 87 ghi chú: Trong đó: kí hiệu (KH) : Biến: 16 biến số rủi ro đầu vào Rất thấp: VL Thấp: L Trung bình: M Cao: H Rất cao: VH Số phiếu: SP Tính trung bình phiếu trọn làm trịn ta PD 85 13 29 47 71 15 32 51 71 PHỤ LỤC H Bảng tính chi tiết số phiếu giá trị trung bình biến số thu từ kết điều tra dành cho PI Biến KH RN N M L VL RN N M L VL RN N M L VL RN N M L VL RN N M L VL RN N M L VL RN N M 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 20 39 25 53 72 65 65 75 75 78 59 78 49 70 60 78 59 75 85 10 35 39 15 15 35 51 78 43 45 62 60 60 60 20 32 65 59 80 80 20 35 75 90 47 65 67 75 70 75 72 75 60 75 69 72 72 66 87 100 10 30 42 49 46 37 38 39 40 41 42 43 44 45 46 47 48 40 63 72 70 36 35 45 53 39 31 39 38 40 56 62 35 45 52 78 50 42 60 42 60 40 55 60 53 47 71 66 95 10 30 59 53 44 60 49 57 62 45 35 45 22 62 60 21 40 40 40 45 63 58 42 60 32 49 45 74 35 23 30 49 50 62 55 52 45 58 39 70 62 62 63 50 62 25 39 43 35 20 35 40 37 47 71 59 60 70 59 78 20 39 40 37 78 70 77 65 75 67 76 80 78 49 63 63 43 67 60 76 40 55 74 25 30 59 45 59 57 60 55 73 70 78 78 65 63 60 65 80 53 62 47 50 45 45 38 45 50 54 60 53 49 75 99 65 100 35 48 77 21 40 20 34 78 65 78 70 100 71 100 57 57 62 20 60 78 37 78 75 67 65 77 82 27 30 52 70 83 25 22 60 24 20 57 47 50 40 63 62 20 47 55 45 62 53 57 65 57 78 71 78 48 40 69 55 42 65 60 62 65 58 53 67 70 80 24 10 31 35 25 39 53 59 61 61 45 49 20 10 22 49 55 19 20 60 47 24 37 63 47 60 62 49 55 57 57 35 40 47 53 53 40 45 42 55 56 62 62 37 SP 31 10 26 4 31 11 24 25 21 18 29 PI 28 51 71 88 14 35 50 69 95 11 31 52 67 31 48 70 84 13 30 53 67 89 16 29 52 70 80 18 32 53 10 11 12 13 14 L VL RN N M L VL RN N M L VL RN N M L VL RN N M L VL RN N M L VL RN N M L VL RN N M L 70 72 78 72 10 25 27 49 47 59 40 40 52 62 22 37 45 41 70 32 60 47 60 55 60 59 44 20 60 68 15 30 35 33 37 32 39 49 49 39 48 78 60 60 60 47 38 39 55 67 78 33 32 40 35 25 50 45 39 40 65 10 35 47 32 27 50 49 57 64 24 35 30 62 49 40 25 49 72 45 52 41 75 80 20 41 45 55 42 39 72 75 57 40 60 78 63 73 75 18 25 40 64 78 10 27 36 61 70 65 76 75 76 37 72 78 75 71 75 62 70 65 78 80 24 22 20 20 49 62 67 60 78 75 90 80 72 70 70 65 78 72 78 78 62 60 78 33 60 50 57 62 59 63 85 74 77 35 25 78 45 78 78 78 65 78 50 40 63 60 60 70 73 78 75 75 72 78 67 66 78 78 82 33 21 59 35 42 20 32 47 10 24 22 30 37 55 58 40 50 70 85 35 80 80 80 93 25 72 68 75 19 40 44 60 81 60 92 10 75 78 77 60 65 67 22 17 82 15 40 40 62 48 60 60 49 55 78 62 43 42 49 45 71 82 70 80 20 30 48 63 59 73 48 63 40 70 40 45 49 40 73 25 50 20 25 63 64 35 40 53 65 78 65 70 60 59 72 38 67 69 72 65 65 65 37 42 45 49 48 49 29 60 35 39 17 40 53 38 57 40 22 4 45 57 62 13 78 72 65 24 1 10 75 78 71 77 63 79 29 49 60 65 65 73 67 27 7 32 39 35 35 13 45 42 19 49 15 75 70 70 75 78 23 73 19 30 50 66 10 33 47 68 80 13 34 51 72 80 24 24 58 71 87 15 37 52 72 83 20 33 50 70 79 20 29 50 69 VL RN N 63 15 M L VL 93 85 90 RN N 16 M 55 45 48 55 L VL ghi chú: Trong đó: kí hiệu (KH) : 90 80 80 21 10 30 30 40 37 50 40 60 20 15 85 22 30 37 75 88 78 79 88 41 58 89 87 13 45 60 40 79 93 40 42 65 37 25 89 85 40 40 10 50 78 85 95 85 21 22 15 20 55 40 45 44 78 Biến: 16 biến số rủi ro đầu vào Rất nhỏ: RN Nhỏ: N Trung bình: M Lớn: L Rất lớn: VL Số phiếu: SP Tính trung bình phiếu trọn làm trịn ta PI 25 44 42 49 58 70 90 92 85 85 90 90 85 85 85 82 40 45 62 78 52 48 65 72 53 52 53 45 72 6 25 25 83 15 34 50 81 87 14 34 49 70 ... tình hình nghiên cứu rủi ro ứng dụng lý thuyết mờ việc đánh giá rủi ro nước giới Tìm hiểu sở lý thuyết mờ cách nhận dạng, đánh giá rủi ro cơng trình ngầm Xây dựng chương trình tính tốn tính hấp... BÙI CÔNG THÀNH i LỜI CẢM ƠN Luận văn ? ?ỨNG DỤNG LÝ THUYẾT MỜ ĐỂ ĐÁNH GIÁ RỦI RO TRONG CƠNG TRÌNH NGẦM” thực từ ngày 01/10/2011 đến ngày 16/04/2012 Mục đích đề tài đánh giá rủi ro Cơng Trình Ngầm, ... 03808504 1-TÊN ĐỀ TÀI: ỨNG DỤNG LÝ THUYẾT MỜ ĐỂ ĐÁNH GIÁ RỦI RO CƠNG TRÌNH NGẦM 2- NHIỆM VỤ VÀ NỘI DUNG: Tìm hiểu cần thiết nghiên cứu rủi ro cơng trình ngầm, cố cơng trình ngầm xảy q khứ Tìm hiểu