Chương trình mô phỏng tính toán thiết kế hệ thống điện vùng sâu, vùng xa, hải đảo

Đại Học Quốc Gia Tp Hồ Chí Minh TRƯỜNG ĐẠI HỌC BÁCH KHOA -! PHẠM HỒNG LĨNH CHƯƠNG TRÌNH MÔ PHỎNG TÍNH TOÁN THIẾT KẾ HỆ THỐNG CẤP ĐIỆN VÙNG SÂU, VÙNG XA, HẢI ĐẢO CHUYÊN NGÀNH : CÔNG NGHỆ NHIỆT MÃ SỐ : 60.52.80 LUẬN VĂN THẠC SĨ TP HỒ CHÍ MINH, THÁNG 10 NĂM 2004 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 NGUYỄN THẾ BẢO Cán chấm nhận xét 1: Cán chấm nhận xét 2: Luận văn thạc só bảo vệ tại: HỘI ĐỒNG CHẤM BẢO VỆ LUẬN VĂN THẠC SĨ TRƯỜNG ĐẠI HỌC BÁCH KHOA TP.HCM, ngày … tháng… năm 2004 Đại Học Quốc Gia Tp Hồ Chí Minh TRƯỜNG ĐẠI HỌC BÁCH KHOA CỘNG HOÀ XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc Lập – Tự Do – Hạnh Phúc NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên: PHẠM HỒNG LĨNH Phái: Nam Ngày, tháng, năm sinh: 14 - 08 - 1979 Nơi sinh: Đồng Nai Chuyên ngành: Công nghệ Nhiệt Mã số: 60.52.80 I TÊN ĐỀ TÀI: CHƯƠNG TRÌNH MÔ PHỎNG TÍNH TOÁN THIẾT KẾ HỆ THỐNG CẤP ĐIỆN VÙNG SÂU, VÙNG XA, HẢI ĐẢO II NHIỆM VỤ VÀ NỘI DUNG: Sơ lược nguồn lượng Cơ sở lý thuyết tính toán hệ thống cấp điện vùng sâu, vùng xa, hải đảo Phần mềm tính toán Đánh giá kết thu từ phần mềm tính toán Đánh giá kết luận III NGÀY GIAO NHIỆM VỤ: 09 - 02 - 2004 IV NGÀY HOÀN THÀNH NHIỆM VỤ: 30 -10 - 2004 V HỌ VÀ TÊN CÁN BỘ HƯỚNG DẪN: TS NGUYỄN THẾ BẢO VI HỌ VÀ TÊN CÁN BỘ CHẤM NHẬN XÉT 1: VII HỌ VÀ TÊN CÁN BỘ CHẤM NHẬN XÉT 2: CÁN BỘ HƯỚNG DẪN CHỦ NHIỆM NGÀNH BỘ MÔN QUẢN LÝ NGÀNH TS NGUYỄN THẾ BẢO PGS.TS LÊ CHÍ HIỆP PGS.TS LÊ CHÍ HIỆP 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 SĐH KHOA QUẢN LÝ NGÀNH LỜI CẢM ƠN Tác giả xin bày tỏ lòng biết ơn sâu sắc chân thành đến TS Nguyễn Thế Bảo tận tình hướng dẫn đóng góp nhiều ý kiến quan trọng quý giá cho nội dung luận văn Tác giả xin chân thành cảm ơn GS, PGS, TS cán Trường Đại Học Bách Khoa, Trường Đại Học Sư Phạm Kỹ Thuật tham gia giảng dạy có nhiều ý kiến đóng góp cho việc hoàn thành nội dung khoa học Luận văn Tác giả xin gửi lời cảm ơn đến cán Trung Tâm Tiết Kiệm Năng Lượng thuộc Sở Khoa Học Công Nghệ Môi Trường TPHCM, Ks Nguyễn Trần Phú công tác Công viên phần mềm Quang Trung, Ks Nguyễn Huy Hoàng công tác Công ty THHH TUICO nhiệt tình giúp đỡ, đóng góp vào thành công chung Luận văn Ngoài ra, tác giả chân thành cảm ơn đến người thân thương, bạn bè, đồng nghiệp hỗ trợ động viên tác giả suốt thời gian thực Luận văn Người thực PHẠM HỒNG LĨNH TÓM TẮT LUẬN VĂN Nội dung Luận văn trình bày số đánh giá nguồn lượng tái tạo gồm pin mặt trời turbine gió Trong phần tiếp theo, Luận văn đưa lý thuyết sử dụng việc tính toán mô tả tính hoạt động hệ thống cấp điện vùng sâu vùng xa thông qua việc sử dụng lượng tái tạo chủ yếu pin mặt trời, turbine gió kết hợp với việc sử dụng máy phát địên dùng xăng dầu kèm theo hệ thống lưu trữ điện số thiết bị phụ khác với việc tính toán tiêu kinh tế kỹ thuật Qua sở toán học ngôn ngữ lập trình Visual Basic, luận văn thể bật chương trình tính toán hệ thống cấp điện với phần mềm tính toán cho hệ thống cấp điện vùng sâu vùng xa Qua đó, phân tích kết rút từ chương trình tính toán để đánh giá đặc tính ưu điểm với khả ứng dụng tương lai hệ thống Các kết chương trình tính tham khảo làm sở cho việc thiết kế, tính toán chọn lựa thiết bị hệ thống để đạt tối ưu chi phí đầu tư ban đầu xem xét đánh giá việc ứng dụng hệ thống cấp điện dùng lượng tái tạo kết hợp sử dụng máy phát diesel Việt Nam MỤC LỤC MỞ ĐẦU CHƯƠNG SƠ LƯC VỀ CÁC NGUỒN NĂNG LƯNG 1.1 Tình hình sử dụng điện Việt Nam 1.2 Giơi thiệu nguồn lượng tái tạo 1.3 Năng lượng tái tạo 1.3.1 Thuỷ điện 1.3.2 Năng lượng mặt trời 1.3.3 Năng lượng gió .12 1.3.4 Năng lượng sinh khối .13 1.4 Các đánh giá ảnh hưởng tới môi trường 16 CHƯƠNG 18 CƠ SỞ LÝ THUYẾT 18 2.1 Xác định tải tiêu thụ 18 2.2 Hệ thống máy phát diesel, xăng 19 2.2.1 Hệ thống máy phát độc lập 19 2.2.2 Hệ thống máy phát sử dụng ắc quy 20 2.3 Hệ thống kết hợp máy phát - ắc quy lượng tái tạo 28 2.3.1 Sơ lược hệ thống pin mặt trời 28 2.3.2 Sơ lược hệ thống điện turbine gió 33 2.3.3 Tính toán hệ thống kết hợp 35 2.4 Tính toán kinh tế 47 2.4.1 Tính toán 47 2.4.2 Tính toán kinh tế 47 2.4.3 Tính chi phí máy phát 47 2.4.4 Chi phí nhiên liệu 48 2.4.5 Tính chi phí ắc quy .49 2.4.6 Tính chi phí pin mặt trời 49 2.4.7 Tính chi phí turbine gioù 50 2.4.8 Tính chi phí inverter, nạp aéc quy 50 2.4.9 Tính chi phí bảo trì vận hành 51 2.4.10.Chi phí sản xuất điện .55 2.5 Hệ thống điện dùng lượng tái tạo kết hợp ắc quy 55 CHƯƠNG 62 PHẦN MỀM TÍNH TOÁN 62 3.1 Khởi động chương trình 63 3.2 Nhập thông số đầu vào 64 3.2.1 Các thông số tải .64 3.2.2 Thông số máy phát .65 3.2.3 Thông số ắc quy .65 3.2.4 Thông số pin mặt trời 65 3.2.5 Thông số turbine gió 66 3.2.6 Thông số thiết bị khác .66 3.2.7 Thông số kinh tế 66 3.2.8 Thoâng số vận hành bảo dưỡng máy phát 68 3.3 Màn hình kết quaû 69 CHƯƠNG 73 ĐÁNH GIÁ CÁC KẾT QUẢ 73 4.1 Kiểm chứng độ tin cậy 73 4.2 Đánh giá tính khả thi số hệ thống 76 4.2.1 Đánh giá ảnh hưởng giá nhiên liệu 78 4.2.2 Đánh giá ảnh hưởng mức độ lạm phát 80 4.2.3 Đánh giá ảnh hưởng mức độ lạm phát giá nhiên liệu .81 4.2.4 Ảnh hưởng suất chiết khấu 82 4.2.5 Đánh giá chung tới chi phí điện .83 4.3 Sự ảnh hưởng lượng tái tạo hệ thống cấp điện kết hợp 83 4.3.1 Ảnh hưởng số lượng thiết bị hay công suất lắp đặt 83 4.3.2 Ảnh hưởng số lượng module tới chi phí điện 85 4.3.3 Ảnh hưởng góc nghiêng pin mặt trời 86 CHƯƠNG 87 ĐÁNH GIÁ - KẾT LUẬN 87 5.1 Khả ứng dụng Việt Nam 87 5.2 Ưu nhược điểm hệ thống điện dùng lượng tái tạo 90 5.2.1 Ưu điểm hệ thống điện dùng lượng tái tạo .91 5.2.2 Nhược điểm hệ thống điện dùng lượng tái tạo .92 5.3 Một số xu hướng phát triển hệ thống điện dùng lượng tái tạo 94 5.3.1 Các công nghệ chế tạo pin mặt trời .94 5.3.2 Công nghệ nano .95 5.3.3 Hệ thống điện mặt trời độc lập 96 5.3.4 Hệ thống điện mặt trời nối lưới 96 5.3.5 Hệ thống điện mặt trời phụ tải lưới điện 97 5.4 Kiến nghị kết luaän 97 KẾT LUẬN 100 PHUÏ LUÏC 103 TÀI LIỆU THAM KHAÛO 136 CÁC CHỮ VIẾT TẮT THƯỜNG DÙNG Ký hiệu Đơn vị Giải thích A Độ đồng chỗ tiếp xúc p-n APR Tỷ số vận tốc định mức turbine B Wh Điện ắc quy Tỷ số vận tốc cắt định mức CPR G W Công suất máy phát d % Suất chiết khấu DOD % Mức độ xả điện E eV Năng lượng hấp thụ bước sóng định F lít Lượng nhiên liệu tiêu hao h Js Hằng số Planck I W/m2 Cường độ xạ Iob W/m2 Bức xạ toàn phần tới bầu khí trái đất Isc Ampe Dòng điện ngắn mạch IT W/m2 Bức xạ tới mặt phẳng nghiêng i % Mức độ lạm phát giá if % Mức độ lạm phát giá nhiên liệu Hằng số Boltzmann K l Wh Tải tiêu thụ Tỷ số công suất turbine OPR P W Công suất thiết bị PW USD Giá trị SFC kWh/l Hệ số sử dụng nhiên liệu T o Nhiệt độ tuyệt đối tk năm Thời gian phân tích kinh tế TOE Tấn dầu tương đương v m/s Vận tốc vR m/s Tốc độ gió định mức turbine vW m/s Vận tốc gió bên V V Điện áp hở mạch Voc V điện áp hở mạch β độ Góc nghiêng thu γ độ Góc phương vị η % Hiệu suất thiết bị λ m Bước sóng ánh sáng φ độ Góc vó độ nơi tính toán, độ K Hệ số phản xạ mặt đất ρ CÁC CHỮ VIẾT TẮT ANN Mạng trí tuệ nhân tạo MPPT Bộ theo dõi công suất pin mặt trời PEFC Chất điện phân polymer PV Pin mặt trời RAPS Hệ thống cấp điện vùng sâu vùng xa SHS Hệ thống nhà mặt trời 124 OPR2 = 6.8772 * APR ^ - 18.97 * APR ^ + 16.743 * APR ^ 4.3451 * APR ^ + 0.3803 * APR - 0.0052 OPR = OPR1 + 1.8 * Abs(OPR2 - OPR1) * (CPR - 2) / Else OPR = 6.8772 * APR ^ - 18.97 * APR ^ + 16.743 * APR ^ - 4.3451 * APR ^ + 0.3803 * APR - 0.0052 End If 'CPR End If 'APR powertam = OPR * Ratedpower 'kWh windpower = powertam End Function CHƯƠNG TRÌNH CON TÍNH TOÁN CÔNG SUẤT PHÁT ĐIỆN TỪ NĂNG LƯNG TÁI TẠO Public Sub TinhSolar() rcst.MoveFirst rcst.MoveNext ngaythu = For thang = To 12 songay = sngthang(Val(thang)) For = To songay ngaythu = ngaythu + For gio = To 24 Windpwr(thang, ngay, gio) = 'Tinh buc xa mat troi Ibxmp = 0.0036 * Val(rcst.Fields("I")) If Ibxmp = Then It(thang, ngay, gio) = Else It(thang, ngay, gio) = Itilt(Val(phi), Val(beta), _ Val(gamma), Val(gio), _ Val(ngaythu), Val(pro), Val(Ibxmp)) / 0.0036 nhietdo = Val(rcst.Fields("TEMP")) hstheond = - (nhietdo - 5) * 0.38 / 100 'Tinh cho nhieu module It(thang, ngay, gio) = It(thang, ngay, gio) * hstheond * CsPhoto * Slmodule * Hsmodule * HsMPPT / 1000 End If rcst.MoveNext Next Next 125 Next thang End Sub Public Sub TinhWind() rcst.MoveFirst rcst.MoveNext For thang = To 12 songay = sngthang(Val(thang)) For = To songay For gio = To 24 It(thang, ngay, gio) = speed = Val(rcst.Fields("WV")) If speed cutoutspeed Then Windpwr(thang, ngay, gio) = Else speed = hieuchinh * speed Windpwr(thang, ngay,gio)= windpower(Val(cutoutspeed), Val(ratedspeed),Val(speed),Val(Ratedpower)) Windpwr(thang,ngay,gio)=Slturbine*Windpwr(thang,ngay ,gio) End If rcst.MoveNext Next Next Next thang End Sub Public Sub TinhSolarWind() rcst.MoveFirst rcst.MoveNext ngaythu = For thang = To 12 songay = sngthang(Val(thang)) For = To songay ngaythu = ngaythu + For gio = To 24 'Tinh buc xa mat troi Ibxmp = 0.0036 * Val(rcst.Fields("I")) If Ibxmp = Then It(thang, ngay, gio) = 126 Else It(thang, ngay, gio) = Itilt(Val(phi), Val(beta), _ Val(gamma), Val(gio), _ Val(ngaythu), Val(pro), Val(Ibxmp)) / 0.0036 nhietdo = Val(rcst.Fields("TEMP")) hstheond = - (nhietdo - 5) * 0.38 / 100 'Tinh cho nhieu module It(thang, ngay, gio) = It(thang, ngay, gio) * hstheond *CsPhoto * Slmodule * Hsmodule * HsMPPT / 1000 End If ' Tinh cong suat gio speed = Val(rcst.Fields("WV")) If speed cutoutspeed Then Windpwr(thang, ngay, gio) = Else speed = hieuchinh * speed Windpwr(thang, ngay, gio) = windpower(Val(cutoutspeed), Val(ratedspeed),Val(speed), Val(Ratedpower)) Windpwr(thang, ngay, gio) = Windpwr(thang, ngay, gio) * Slturbine End If rcst.MoveNext Next gio Next Next thang End Sub CHƯƠNG TRÌNH CON TÍNH TOÁN HỆ THỐNG KẾT HP MÁY PHÁT ẮC QUY VÀ NĂNG LƯNG TÁI TẠO Public Sub TinhBGPT() Dim Btamtruoc As Double Sumload = For gio = giomo To giongung Sumload = Sumload + Load(gio) Next 127 SlNoitiep = ltrlen(Val(DienapHT / DienapBat)) 'So luong Battery mac theo day Ahngay = (Diennangtt - Sumload) / (HsInverter * DienapHT * DOD) CsWhBat = CsAhBat * DienapBat 'Tinh cong suat dinh muc Watt hour If Sumload Diennangtt Then SlSongsong = ltrlen(Val(Ahngay / CsAhBat)) CsdmucBat = SlSongsong * DienapHT * CsAhBat DODtt = DOD * Ahngay * DienapHT / CsdmucBat ' Dien nang cua Battery bat dau nap BatEnergy = (1 - DOD) * CsdmucBat End If ' Xac dinh dien nap Battery vao thoi diem dau If giomo Then dem = giomo - Battery(1, 1, dem) = BatEnergy Do While dem > Battery(1, 1, dem - 1) = Battery(1, 1, dem) + Load(dem) / HsInverter _ - (It(1, 1, dem) * HsMPPT - Windpwr(1, 1, dem)) * HsBat dem = dem - Loop ElseIf giomo = And giongung 24 Then Battery(1, 1, 1) = BatEnergy Else SlSongsong = CsdmucBat = DienapHT * CsAhBat BatEnergy = (1 - DOD) * CsdmucBat Battery(1, 1, 1) = BatEnergy End If ' Xac dinh dien nap Battery vao thoi tiep theo Tai30Eng = 0.5 * CsdmucEng Tai90Eng = 0.9 * CsdmucEng Tai60Eng = 0.6 * CsdmucEng HsBBC = HsBat * HsBC Momay = False Ngung = False ' SolarUsed = ' WindUsed = SumUsed = demSFC = TotalSFC = 128 For thang = To 12 UsedMonthly(thang) = countPrg = / 12 ValuePrg = (countPrg / 25) * 100 + ValuePrg frmMainRaps.prgDangtinh.Value = ValuePrg songay = sngthang(Val(thang)) For = To songay UsedDaily(thang, ngay) = For gio = To 24 Renewable = It(thang, ngay, gio) + Windpwr(thang, ngay, gio) Lgio = Load(gio) RenLgio = Renewable - Lgio / HsInverter If (gio >= giomo) And (gio giomo) Then Btamtruoc = Battery(thang, ngay, gio - 1) If Lgio Then Btam = Btamtruoc + Tai90Eng * HsBBC + RenLgio * HsBat Call TestCSdmuc Etam = (Btam - Btamtruoc - RenLgio * HsBat) / HsBBC If Etam < Then Etam = End If Else Btam = Btamtruoc + (Tai90Eng + RenLgio * HsInverter) * HsBBC Call TestCSdmuc Etam = (Btam - Btamtruoc) / HsBBC - RenLgio * HsInverter End If Battery(thang, ngay, gio) = Btam LoadEng(thang, ngay, gio) = Etam loadtam = Etam / CsdmucEng SFCFuel(thang, ngay, gio) = SCF(Val(loadtam)) demSFC = demSFC + TotalSFC = TotalSFC + SFCFuel(thang, ngay, gio) Fuelcnsm(thang, ngay, gio) = Etam / (1000 * SFCFuel(thang, ngay, gio)) End If 'Kiem tra tai tieu thu so voi tai may phat ElseIf thang = And = And (gio = Or gio < giomo) Then Battery(thang, ngay, gio) = Battery(thang, ngay, gio) Else Btam = TtaiBat(Val(thang), Val(ngay), Val(gio)) If Btam > CsdmucBat Then Btam = CsdmucBat ElseIf Btam < Then Ngung = True 130 Exit Sub End If Battery(thang, ngay, gio) = Btam End If ' Tinh dien nang su dung tu nang luong tai tao Etam = LoadEng(thang, ngay, gio) If Btam < CsdmucBat Or (Btam = CsdmucBat And Etam 0) Then Usedtam = Renewable Else If gio Then Btamtruoc = Battery(thang, ngay, gio - 1) Usedtam = (Btam - Btamtruoc) / HsBat + Load(gio) / HsInverter ElseIf gio = And Then Btamtruoc = Battery(thang, - 1, 24) Usedtam = (Btam - Btamtruoc) / HsBat + Load(gio) / HsInverter ElseIf gio = And = And thang Then Btamtruoc = Battery(thang - 1, Sngthangtruoc, 24) Usedtam = (Btam - Btamtruoc) / HsBat + Load(gio) / HsInverter Else Usedtam = End If End If Used(thang, ngay, gio) = Usedtam UsedDaily(thang, ngay) = UsedDaily(thang, ngay) + Used(thang, ngay, gio) Next gio UsedMonthly(thang) = UsedMonthly(thang) + UsedDaily(thang, ngay) Next Sngthangtruoc = songay SumUsed = SumUsed + UsedMonthly(thang) Next thang End Sub CHƯƠNG TRÌNH CON TÍNH TOÁN KINH TẾ Public Sub GanKinhte() timeanalys = Val(frmMainRaps.txtKTime) interest = Val(frmMainRaps.txtKtlaisuat) / 100 inflation = Val(frmMainRaps.txtKtlamphat) / 100 fuelinflation = Val(frmMainRaps.txtKLamphatNL) / 100 End Sub Public Sub Tinhkinhte() Dim TestVH As Boolean Call GanKinhte 131 TestVH = (frmMainRaps.chkVanhanh.Value) ValuePrg = (1 / 25) * 100 + ValuePrg frmMainRaps.prgDangtinh.Value = ValuePrg If chkGen And TestVH Then Call TinhCostGen Call TinhCostVH ElseIf chkGen Then Call TinhCostGen TotalPWVH = PWVC = InicostVH = InicostVC = Else TotalPWVH = PWVC = IniCostGen = InicostVH = InicostVC = PWGen = PWGen = PWFuel = End If ' Tinh Battery Call TinhCostBat ' Tinh Solar module Call TinhCostPho ' Tinh wind turbine TinhCostTurb ' Tong chi phi dau tu ban dau TotalIniCost = IniCostGen + IniCostMod + IniCostTurb _ + IniCostBC + IniCostBat + IniCostInv + InicostVH ' Tinh tong gia tri hien tai cua chi phi dau tu TotalCapCost = Round(PWGen + PWBat + PWMod + PWTurb + PWInv _ + PWBC) ' Tong gia tri hien tai TotalPreCost = Round(TotalCapCost + PWFuel + PWVC + TotalPWVH) ' Chi phi kWh dien PrElecgenerate = Round(TotalPreCost * 1000 / (timeanalys Diennangtt * 365), 4) End Sub * Public Sub TinhCostVH() Call GanCPVanhanh InicostVH = costPhukien + costDaydien + costDK + costLapdat + costKhac PWVC = For time = To timeanalys PWVC = PWVC + costLaodong * PreWorth(Val(inflation), Val(interest), Val(rdown(Val(time)))) Next time SolanDtu = lifeGen / timeDT SolanDtu = ltrlen(Val(SolanDtu)) - SolanTtu = timeDT / timeTT SolanTtu = ltrlen(Val(SolanTtu)) - SolanBD = timeTT / timeBD SolanBD = ltrlen(Val(SolanBD)) - yearDT = yearTT = yearBD = PWDT = PWTT = PWBD = PWLD = For time = To noGenreplace 132 year = time * lifeGen / gioHD If time < noGenreplace Then Call TTheGen Else Call TTheLancuoi End If Next time TotalPWVH = PWDT + PWTT + PWBD + PWLD End Sub Public Sub TTheGen() For lanDT = To SolanDtu yearDT = year + timeDT * lanDT / gioHD If lanDT < SolanDtu Then Call BatdauDT Else DTtam = lifeGen - timeDT * lanDT Call KetthucDT Exit Sub End If Next lanDT End Sub Public Sub TTheLancuoi() Reptam = timeanalys * gioHD - time * lifeGen For lanDT = To SolanDtu If lanDT * timeDT < Reptam Then yearDT = year + timeDT * lanDT / gioHD Call BatdauDT Else DTtam = Reptam - timeDT * lanDT Call KetthucDT End If Next lanDT End Sub Public Sub BatdauDT() If lanDT > Then PWDT = PWDT + costDT * PreWorth(Val(inflation), Val(interest), Val(rdown(Val(yearDT)))) PWLD = PWLD + costLaodong * PreWorth(Val(inflation), Val(interest), Val(rdown(Val(yearDT)))) End If For lanTT = To SolanTtu yearTT = yearDT + timeTT * lanTT / gioHD If lanTT < SolanTtu Then 'Trung tu va bao duong If lanTT > Then PWTT = PWTT + costTT * PreWorth(Val(inflation), Val(interest), Val(rdown(Val(yearTT)))) End If For lanBD = To SolanBD yearBD = yearTT + timeBD * lanBD / gioHD PWBD = PWBD + costBD * PreWorth(Val(inflation), Val(interest), Val(rdown(Val(yearBD)))) Next lanBD Else TTtam = timeDT - lanTT * timeTT If TTtam > timeBD Then ' Trung tu PWTT = PWTT + costTT * PreWorth(Val(inflation), Val(interest), Val(rdown(Val(yearTT)))) For lanBD = To SolanBD If lanBD * timeBD < TTtam Then yearBD = yearTT + timeBD * lanBD / gioHD 133 PWBD = PWBD + costBD * PreWorth(Val(inflation), Val(interest), Val(rdown(Val(yearBD)))) Else: Exit Sub End If Next lanBD ElseIf TTtam < timeBD And TTtam > Then yearBD = yearTT 'Chi bao duong End If End If Next lanTT End Sub Public Sub KetthucDT() If DTtam timeBD Then 'Trung tu If lanTT > Then PWTT = PWTT + costTT * PreWorth(Val(inflation), Val(interest), Val(rdown(Val(yearTT)))) End If For lanBD = To SolanBD If lanBD * timeBD < TTtam Then yearBD = yearTT + timeBD * lanBD / gioHD PWBD = PWBD + costBD * PreWorth(Val(inflation), Val(interest), Val(rdown(Val(yearBD)))) Else: Exit For End If Next lanBD ElseIf TTtam > And TTtam < timeBD Then yearBD = yearTT 'Chi bao duong PWBD = PWBD + costBD * PreWorth(Val(inflation), Val(interest), Val(rdown(Val(yearBD)))) Exit For End If Next lanTT End If End Sub Public Sub TinhCostGen() sumWhGen = CsdmucEng * 0.6 * lifeGen gioHD = For thang = To 12 songay = sngthang(Val(thang)) 134 For = To songay For gio = To 24 If LoadEng(thang, ngay, gio) Then gioHD = gioHD + End If Next gio Next Next thang ' ltimeGen = ltrlen(Val(sumWhGen / SumElGen)) ltimeGen = ltrlen(Val(lifeGen / gioHD)) noGenreplace = ltrlen(Val(timeanalys / ltimeGen)) - IniCostGen = Round(costGen * CsWEngine) PWGen = IniCostGen For noReplace = To noGenreplace PWGen = PWGen + IniCostGen * PreWorth(Val(inflation), Val(interest), Val(ltimeGen * noReplace)) Next noReplace ' Tinh Fuel sumcostFuel = Round(SumFuel * costFuel) PWFuel = sumcostFuel * PWF(Val(fuelinflation), Val(interest), Val(timeanalys)) End Sub Public Sub TinhCostBat() If TestBat Then IniCostBat = Round(CsdmucBat * costBat) ' yearBat = (-3457.9439 * DOD + 4219.6) / 365 ' yearBat = ltrlen(Val(yearBat)) noBatreplace = ltrlen(Val(timeanalys / lifeBat)) - PWBat = IniCostBat For noReplace = To noBatreplace PWBat = PWBat + IniCostBat * PreWorth(Val(inflation), Val(interest), Val(yearBat * noReplace)) Next noReplace ' Tinh Battery charger IniCostBC = Round(CsdmucEng * costBC) noBCreplace = ltrlen(Val(timeanalys / lifeBC)) - PWBC = IniCostBC For noReplace = To noBCreplace PWBC = PWBC + IniCostBC * PreWorth(Val(inflation), Val(interest), Val(lifeBC * noReplace)) Next noReplace ' Tinh Inverter IniCostInv = Round(CsInverter * costInv) noInvreplace = ltrlen(Val(timeanalys / lifeInv)) - PWInv = IniCostInv For noReplace = To noInvreplace PWInv = PWInv + IniCostInv * PreWorth(Val(inflation), Val(interest), Val(lifeInv * noReplace)) Next noReplace Else IniCostBat = IniCostBC = IniCostInv = PWBat = PWInv = PWBC = End If End Sub Public Sub TinhCostPho() If TestPho Then IniCostMod = Round(CsPhoto * costMod * Slmodule) noModreplace = ltrlen(Val(timeanalys / lifeMod)) - PWMod = IniCostMod For noReplace = To noModreplace 135 PWMod = PWMod + IniCostMod * PreWorth(Val(inflation), Val(interest), Val(lifeMod * noReplace)) Next noReplace Else IniCostMod = PWMod = End If End Sub Public Sub TinhCostTurb() If TestTurb Then IniCostTurb = Round(Ratedpower * costTurb * Slturbine) noTurbreplace = ltrlen(Val(timeanalys / lifeTurb)) - PWTurb = IniCostTurb For noReplace = To noModreplace PWTurb = PWTurb + IniCostTurb * PreWorth(Val(inflation), Val(interest), Val(lifeTurb * noReplace)) Next noReplace Else IniCostTurb = PWTurb = End If End Sub 136 TÀI LIỆU THAM KHẢO A Akbarzadeh, “Fundamentals of Remote Area Power Supply System”, Energy Victoria 1992 HP Garg vaø J Prakash, “Solar Energy”, 1996 J.A Duffie and W.A Beckman, “Solar Energy Thermal Processes”, Wiley Interscience 1974 Clark, D.R, S.A Klein and W.A Beckman, “Solar Energy”, 1984 R.J Van Overstraeten, “Advanced in Silicon Solar Cell Processing”, D Riedel Publishing Company 1980 Evans et al, “Simplified design guide for estimating photovoltaic flat array and system performance”, Sandia National Laboratories 1980 J Gordon, “Optimal sizing of stand – alone photovoltaic power systems”, 1987 “Solar 96 Conference Proceedings”, Australian and New Zealand solar energy Society 1996 Dr Trevor Pryor, “Energy Systems M390”, Murdoch University 1994 10 F.R Eldrige, “Wind machines”, Nostrand Co 1980 11 J Jensen and B Sorensen, ”Fundamentals of Energy Store”, New York 1980 12 C.L Mantell, “ Batteries and Energy Systems”, McGraw Hill 1983 13 Shugar, D.S., “Photovoltaic in the Utility Distribution System: The Evaluation of System and Distributed” 14 Hanazawa, K., et al., “Purification of Metallurgical Grade Silicon by Electron Beam Melting,” Proceedings of the Intern Photo Science and Engineering Conference, Miyazaki, Japan, (November 1996) 15 Shugar, D.S., M.G Real, and P Aschenbrenner, “Comparison of Selected Economics Factors for Large Ground-Mounted Photovoltaic Systems with Roof-Mounted Photovoltaic Systems in Switzerland and the USA,” Proceedings of the 11th E.C Photo Solar Energy Conf (1992) 137 16 Osborn, D.E., and D.E Collier, “Utility Grid-connected Photovoltaic Distributed Power Systems,” Proceedings of the National Solar Energy Conf., ASES Solar 96, Asheville, NC (April 1996) 17 Strong, S.J., “Power Windows: Building-integrated Photovoltaics,” IEEE Spectrum, October 1996 18 Lamarre, L., "Utility Customers Go for Green," EPRI Journal, March/April 1997 19 “Renewable energy technology in Asia”,Viện lượng 2003 20 “Hội thảo sử dụng hiệu lượng bảo vệ môi trường”, TTTKNL 2003 21 Một số seminar sử dụng hiệu lượng 22 Một số website http://www.energy.ca.gov/reports http://www.solarlectrix.com http://www.acourt.co.nz http://www.pege.org http://www.fsec.ucf.edu http://www.bpsolar.com http://www.aseanenergy.org http://www.self.org http://www.wordbank.org.vn TÓM TẮT LÝ LỊCH TRÍCH NGANG Họ tên: PHẠM HỒNG LĨNH Ngày, tháng, năm sinh: 14 - 08 - 1979 Nơi sinh: Đồng Nai Địa liên lạc: Trung tâm Kiểm định Kỹ thuật an toàn Khu vực 193A, Lý Chính Thắng Q.3, Tp.Hồ Chí Minh Trạm Kiểm định Linh Trung Lô 108, Khu chế xuất Linh Trung, Thủ Đức, Tp.Hồ Chí Minh QUÁ TRÌNH ĐÀO TẠO: 1997-2002: Sinh viên trường Đại học Bách Khoa Thành phố Hồ Chí Minh 2002-2004: Học viên Cao học khoá 13 Trường Đại học Bách Khoa Thành phố Hồ Chí Minh QUÁ TRÌNH CÔNG TÁC: 2002: Công tác Công ty Cơ nhiệt lạnh Searefico, Q.5, Tp HCM 2003: Công tác Công ty TNHH Cơ nhiệt TST, Q Tân Bình, Tp.HCM 2004 - đến nay: Trung tâm Kiểm định Kỹ thuật an toàn Khu vực ... Phần mềm tính toán thiết kế hệ thống điện RAPS gồm số hệ thống cấp điện đơn giản hệ thống máy phát điện diesel, hệ thống điện dùng lượng mặt trời, hệ thống điện dùng lượng gió số hệ thống kết hợp... phát triển hệ thống điện vùng sâu, vùng xa (RAPS), việc tính toán thiết kế với hỗ trợ máy tính thông qua phần mềm tính toán đơn giản giúp người thiết kế tính toán, mô hệ thống cấp điện phương... văn thể bật chương trình tính toán hệ thống cấp điện với phần mềm tính toán cho hệ thống cấp điện vùng sâu vùng xa Qua đó, phân tích kết rút từ chương trình tính toán để đánh giá đặc tính ưu điểm