(Luận văn thạc sĩ) tính toán hiệu quả của hệ thống phát điện kết hợp từ thủy động, địa nhiệt và nhiệt mặt trời

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH LUẬN VĂN THẠC SĨ TRẦN THÀNH NIỆM TÍNH TỐN HIỆU QUẢ CỦA HỆ THỐNG PHÁT ÐIỆN KẾT HỢP TỪ THỦY ÐỘNG, ÐỊA NHIỆT VÀ NHIỆT MẶT TRỜI NGÀNH: KỸ THUẬT ÐIỆN - 60520202 SKC005809 Tp Hồ Chí Minh, tháng 10/2018 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƢỜNG ĐẠI HỌC SƢ PHẠM KỸ  THUẬT THÀNH PHỐ HỒ CHÍ MINH - - LUẬN VĂN THẠC SĨ TRẦN THÀNH NIỆM TÍNH TỐN HIỆU QUẢ CỦA HỆ THỐNG PHÁT ĐIỆN KẾT HỢP TỪ THỦY ĐỘNG, ĐỊA NHIỆT VÀ NHIỆT MẶT TRỜI NGÀNH: KỸ THUẬT ĐIỆN – 60520202 Tp Hồ Chí Minh, tháng 10/2018 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƢỜNG ĐẠI HỌC SƢ PHẠM KỸ  THUẬT THÀNH PHỐ HỒ CHÍ MINH - - LUẬN VĂN THẠC SĨ TRẦN THÀNH NIỆM TÍNH TỐN HIỆU QUẢ CỦA HỆ THỐNG PHÁT ĐIỆN KẾT HỢP TỪ THỦY ĐỘNG, ĐỊA NHIỆT VÀ NHIỆT MẶT TRỜI NGÀNH: KỸ THUẬT ĐIỆN – 60520202 Hƣớng dẫn khoa học: PGS.TS LÊ CHÍ KIÊN Tp Hồ Chí Minh, tháng 10/2018 Luận văn LÝ LỊCH KHOA HỌC I LÝ LỊCH SƠ LƢỢC: Họ tên: Trần Thành Niệm Giới tính: Nam Ngày, tháng, năm sinh: 26 – 04 – 1986 Nơi sinh: Tiền Giang Quê quán: Tiền Giang Dân tộc: Kinh Địa liên lạc: Tổ 18, Ấp 3, Phú An, Cai Lậy, Tiền Giang Điện thoại riêng: 0906 995 034 E-mail: niempctg@gmail.com II QUÁ TRÌNH ĐÀO TẠO: Đại học: Hệ đào tạo: Khơng quy Thời gian đào tạo từ 10/2005 đến 10/2010 Nơi học (trường, thành phố): Trường đại học Bách Khoa thành phố Hồ Chí Minh Ngành học: Hệ Thống Điện Tên đồ án, luận án môn thi tốt nghiệp: Thiết kế phần điện – Bảo vệ chống sét cho nhà máy nhiệt điện trạm biến áp Ngày & nơi bảo vệ đồ án, luận án thi tốt nghiệp: Năm 2010, trường đại học Bách Khoa thành phố Hồ Chí Minh Người hướng dẫn: Ts Đinh Hoàng Bách Thạc sĩ: Hệ đào tạo: Chính quy Thời gian đào tạo từ 10/2016 đến 10/2018 Nơi học (trường, thành phố): Trường Đại học Sư phạm kỹ thuật thành phố Hồ Chí Minh Ngành học: Kỹ thuật điện Tên luận văn: Tính tốn hiệu hệ thống phát điện kết hợp Từ thủy động, địa nhiệt nhiệt mặt trời HVTH: Trần Thành Niệm i GVHD: PGS.TS Lê Chí Kiên Luận văn Ngày & nơi bảo vệ luận văn: Tháng 10/2018 Trường đại học Sư phạm kỹ thuật thành phố Hồ Chí Minh Người hướng dẫn: PGS TS Lê Chí Kiên BI QUÁ TRÌNH CƠNG TÁC CHUN MƠN KỂ TỪ KHI TỐT NGHIỆP ĐẠI HỌC: Thời gian 04/2011 đến HVTH: Trần Thành Niệm ii GVHD: PGS.TS Lê Chí Kiên Luận văn LỜI CAM ĐOAN Tơi cam đoan cơng trình nghiên cứu Các số liệu, kết nêu luận văn trung thực chưa cơng bố cơng trình khác Tp Hồ Chí Minh, ngày 01 tháng 10 năm 2018 (Ký tên ghi rõ họ tên) HVTH: Trần Thành Niệm iii GVHD: PGS.TS Lê Chí Kiên Luận văn CẢM TẠ Sau thời gian học tập nghiên cứu hướng dẫn tận tình thầy PGS.TS Lê Chí Kiên, học viên hoàn thành luận văn thời hạn mục tiêu đề ban đầu Học viên xin trân trọng bày tỏ lòng biết ơn sâu sắc đến thầy PGS.TS Lê Chí Kiên kiến thức quý báu phương pháp nghiên cứu mà thầy truyền đạt Học viên xin gửi lòng tri ân đến quý thầy cô trường Đại Học Sư Phạm Kỹ Thuật Tp.Hồ Chí Minh trang bị cho học viên lượng kiến thức bổ ích, đặc biệt xin chân thành cảm ơn q thầy Khoa Điện – Điện Tử tạo điều kiện thuận lợi hỗ trợ cho học viên nhiều trình học tập thời gian làm luận văn Cuối học viên xin gửi lời cảm ơn đến đồng nghiệp, gia đình, bạn bè hỗ trợ học viên mặt tinh thần để học viên hoàn thành luận văn HVTH: Trần Thành Niệm iv GVHD: PGS.TS Lê Chí Kiên Luận văn TĨM TẮT Hiện nay, ngành điện nước ta đối mặt với nhiệm vụ khó khăn cung ứng nhu cầu ngày cao khách hàng Bên cạnh vấn đề nhiễm môi trường, khủng hoảng cung cầu lượng diễn gay gắt Việt Nam mà xuất nhiều quốc gia giới Nhận thức điều quốc gia giới nghiên cứu đưa vào nguồn lượng mới, lượng tái tạo Bao gồm lượng mặt trời, lượng gió, địa nhiệt, khí sinh học… Tuy khơng phổ biến nguồn lượng tái tạo khác lượng địa nhiệt lượng mặt trời nguồn lượng sạch, thân thiện gần vơ tận, đáp ứng khát lượng Nguồn lượng khơng gây nhiễm khí hậu hay mơi trường Gần đây, công nghệ phát điện máy phát MHD có bước phát triển Nó ứng dụng rộng rãi chu trình phát điện hỗn hợp để nâng cao hiệu suất nhà máy điện Đề tài: "Tính tốn hiệu hệ thống phát điện kết hợp Từ thủy động, địa nhiệt nhiệt mặt trời" nghiên cứu với mục tiêu tính toán so sánh hiệu suất phát điện hiệu kinh tế hai mơ hình phát điện kết hợp máy phát điện Từ thủy động với địa nhiệt nhiệt mặt trời Đề tài định hướng nghiên cứu để xây dựng nhà máy địa nhiệt điện điện mặt trời hiệu suất cao phù hợp với kinh tế kỹ thuật nước ta HVTH: Trần Thành Niệm v GVHD: PGS.TS Lê Chí Kiên Luận văn ABSTRACT At present, our country's power sector is facing a very difficult task that is supplying the increasing demand of our customers Besides the environmental pollution problem, a crisis of demand and supply of energy is happening very harshly not only in Vietnam but also appeared in many countries around the world Recognizing that countries around the world are included in the research of new energy sources, renewable energy Including solar, wind, geothermal, biogas Although not as common as other renewable energy sources but geothermal energy and solar energy are a clean energy sources, friendly and virtually inexhaustible, could meet today's energy thirst This energy does not pollute the environment or climate Recently, power generation technology of MHD generators with the new developments It has been widely used in the power generation mix cycle to improve the efficiency of power plants Topic: “Miscalculate outcome of combinary systems MHD generators, geothermal and solar energy” The purpose of the study was to compare the differences in economic efficiency of MHD with geothermal energy and MHD with solar energy The themeoriented research also recommended the construction of geothermal power or solar power plants suitable for high performance and technical economy of the country HVTH: Trần Thành Niệm vi GVHD: PGS.TS Lê Chí Kiên Luận văn Q2B=G2*Cp*T2B Q3B=G2*Cp*T3B Q4B=G2*Cp*T4B %dien nang may phat MHD PMHD=Q3M*nuyEE %cong suat may nen N1=1.5;%thong so may nen N2=1.75;% thong so may nen Pnen1=N1*(T2M-T1M)*G1*Cp Pnen2=N2*(T2B-T1B)*G2*Cp %cong suat ion hoa MHD nion=0.5; MassNo=0.004; Seed=0.0001; Wion=G1/(MassNo*Seed*60220*1620*12.13) Pion=Wion/nion %cong suat tuabin Pkhi=Q3B-Q4B %cong suat tuabin hoi nuyHOI=0.15;%hieu suat tuabin hoi cai dat PHOI=((Q5M-Q1M)+(Q4B-Q1B))*nuyHOI %hieu suat chu trinh Ptongphat=PMHD+Pkhi+PHOI Ptongthu=Pion+Pnen1+Pnen2 nuytong=(Ptongphat-Ptongthu)/Qin1 %tinh entropy Tref=298; Pref=1.03e+005; %chu trinh MHD SM=(Cp*log(T3M/Tref)-Cp*((y-1)/y)*log(P3M/Pref))*G1 S3M=SM+(Cp*log(T3M/T2M)-Cp*((y-1)/y)*log(P3M/P2M))*G1 S4M=S3M+(Cp*log(T4M/T3M)-Cp*((y-1)/y)*log(P4M/P3M))*G1 S5M=S4M+(Cp*log(T5M/T4M)-Cp*((y-1)/y)*log(P5M/P4M))*G1 S1M=S5M+(Cp*log(T1M/T5M)-Cp*((y-1)/y)*log(P1M/P5M))*G1 S2M=S1M+(Cp*log(T2M/T1M)-Cp*((y-1)/y)*log(P2M/P1M))*G1 %chu trinh brayton SB=(Cp*log(T3B/Tref)-Cp*((y-1)/y)*log(P3B/Pref))*G2 S3B=SB+(Cp*log(T3B/T2B)-Cp*((y-1)/y)*log(P3B/P2B))*G2 S4B=S3B+(Cp*log(T4B/T3B)-Cp*((y-1)/y)*log(P4B/P3B))*G2 S1B=S4B+(Cp*log(T1B/T4B)-Cp*((y-1)/y)*log(P1B/P4B))*G2 HVTH: Trần Thành Niệm 79 GVHD: PGS.TS Lê Chí Kiên Luận văn S1B=S1M;% thong so dau noi hoi S2B=S1B+ (Cp*log(T2B/T1B)-Cp*((y-1)/y)*log(P2B/P1B))*G2 x1=[S3M S4M S5M S1M S2M S3M]; y1=[T3M T4M T5M T1M T2M T3M]; x2=[S3B S4B S1M S2B S3B]; y2=[T3B T4B T1M T2B T3B]; h=plot(x1,y1,'b-*',x2,y2,'b-+'); grid xlabel('Truc Entropy') ylabel('Truc Nhiet Do (K)') - Chƣơng trình tính tốn mơ thơng số chu trình Brayton dùng NLMT %phan tich chu trinh NLMT hon hop truyen thong %thong so ban dau nuyquanghoc=0.65 Tthuc= 3310; Tthu=nuyquanghoc*Tthuc nuythu=0.93; Tbothu=Tthu*nuythu T3B=2000; P3B=4e+005; y=1.6667; ns=0.93; et=0.87; Qin=1.00e+007;% nang luong NLMT thu 10Kmvuong nuyS=0.3; Pbom=3.00e+004; %phan tich cac nut %nut so %tuabin hoi tan N=18:1 N=(1/18); P4B=N*P3B totuabin=N^((et*(y-1)/y)) T4B=T3B*totuabin %nut so deltaHN=0.01; nuyHN=1; P1B=P4B*(1-deltaHN*nuyHN) T2T=600; HVTH: Trần Thành Niệm 80 GVHD: PGS.TS Lê Chí Kiên Luận văn T1T=281.58; T1B=T4B-((T2T-T1T)/(1-deltaHN)) %nut so tysonen=6; P2B=P1B*tysonen tonen=1+(1/ns) T2B=T1B*tonen %phan tich nhiet luong Cp=5196.5; G=Qin/(Cp*(T3B-T2B)) %nhiet luong tai tung nut chu trinh Q1B=G*Cp*T1B Q2B=G*Cp*T2B Q3B=G*Cp*T3B Q4B=G*Cp*T4B %cong suat tuabin P1=Q3B-Q4B %cong suat tubin hoi QinS=Q4B-Q1B P2=nuyS*QinS %cong suat may nen hsn=2 Pnen=hsn*(T2B-T1B)*G*Cp %hieu suat chu trinh nuytong=(P1+P2-Pnen)/Qin %tinh entropy Tref=298; Pref=1.03e+005; SB=(Cp*log(T3B/Tref)-Cp*((y-1)/y)*log(P3B/Pref))*G S3B=SB+(Cp*log(T3B/T2B)-Cp*((y-1)/y)*log(P3B/P2B))*G S4B=S3B+(Cp*log(T4B/T3B)-Cp*((y-1)/y)*log(P4B/P3B))*G S1B=S4B+(Cp*log(T1B/T4B)-Cp*((y-1)/y)*log(P1B/P4B))*G S2B=S1B+(Cp*log(T2B/T1B)-Cp*((y-1)/y)*log(P2B/P1B))*G x=[S3B S4B S1B S2B S3B]; y=[T3B T4B T1B T2B T3B]; h=plot(x,y,'b-*'); grid xlabel('Truc Entropy') ylabel('Truc Nhiet Do (K)') HVTH: Trần Thành Niệm 81 GVHD: PGS.TS Lê Chí Kiên Luận văn - Chƣơng trình tính tốn mơ thơng số chu trình LMMHD dùng NLMT: % chu trinh LM MHD NLMT % phan tich chu trinh y=110;% he so dan nhiet cua kim loai long Tbothu=1000; Qin=1.0e+005;% nhiet luong NLMT 1Kmvuong P2=8.0e+004; T2=890; Ptron=0.2e+006; nEE=0.3; nMHD=0.8; deltaQMHD=0.005; ntach=0.9; deltaQTDN=0.01; nTDN=1; T1s=298;%thong so chu trinh hoi T2s=598;%thong so chu trinh hoi tonen=2; to=1.25; Cp=982; %nhiet dung rieng NaK deltaTtron=5; deltaTtach=5; nuyS=0.3; % nut T3=T2deltaTtron P3=Ptron %nut so T4=T3*(1-nEEdeltaQMHD) P4=P3*((1(nEE/nMHD))^(y/(y-1))) %nut so P5=P4/2 P5a=P5 T5=T4-deltaTtach %nut so P1=P5*(1-deltaQTDN*nTDN) T1=T5-((T2s-T1s)/(1-deltaQTDN)) %nut so P2=P1*tonen T2=T1*to HVTH: Trần Thành Niệm 82 GVHD: PGS.TS Lê Chí Kiên Luận văn %chu trinh brayton G=Qin/(Cp*(T3-T2)) %nhiet luong tai tung nut chu trinh Q1=G*Cp*T1 Q2=G*Cp*T2 Q3=G*Cp*T3 Q4=G*Cp*T4 Q5=G*Cp*T5 % cong suat may phat MHD PMHD=Q3*nEE % cong suat tuabin hoi Ptuabin=(Q5-Q1)*nuyS % cac cong suat tieu thu chu trinh hesobom=1.25;%thong so tra cuu Ptron=6.0e+003;%thong so tra cuu Ptach=4.0e+003;%thong so tra cuu Pbom=hesobom*(T2-T1)*G*Cp Pbomdt=250; %hieu suat he thong LMMHD nuytong=(PMHD+PtuabinPtron-Ptach-Pbom-Pbomdt)/Qin %tinh entropy Tref=298; Pref=1.03e+005; %chu trinh MHD S=(Cp*log(T3/Tref)-Cp*((y-1)/y)*log(P3/Pref))*G S3=S+(Cp*log(T3/T2)-Cp*((y-1)/y)*log(P3/P2))*G S4=S3+(Cp*log(T4/T3)-Cp*((y-1)/y)*log(P4/P3))*G S5=S4+(Cp*log(T5/T4)-Cp*((y-1)/y)*log(P5/P4))*G S1=S5+(Cp*log(T1/T5)-Cp*((y-1)/y)*log(P1/P5))*G S2=S1+(Cp*log(T2/T1)-Cp*((y-1)/y)*log(P2/P1))*G x=[S3 S4 S5 S1 S2 S3]; y=[T3 T4 T5 T1 T2 T3]; h=plot(x,y,'b-*'); grid xlabel('Truc Entropy') ylabel('Truc Nhiet Do (K)') - Chƣơng trình tính tốn mơ thơng số chu trình LMMHD dùng nhiệt thải: % chu trinh LM MHD Nhiet thai HVTH: Trần Thành Niệm 83 GVHD: PGS.TS Lê Chí Kiên Luận văn % phan tich chu trinh y=110;% he so dan nhiet cua kim loai long Qin=1.0e+005; P2=8.0e+004; T2=890; Ptron=0.2e+006; ntach=0.9; nEE=0.3; nMHD=0.8; deltaQMHD=0.005; deltaQngung=0.01; nuyTngung=1; Tin=300;%thong so vao TDN Tout=470;%thong so vao TDN tonen=2; to=1.25; Cp=982; %nhiet dung rieng NaK deltaTtron=5; deltaTtach=5; % nut T3=T2-deltaTtron P3=Ptron %nut so T4=T3*(1-nEE-deltaQMHD) P4=P3*((1-(nEE/nMHD))^(y/(y-1))) %nut so P5=P4/2 P5a=P5 T5=T4-deltaTtach %nut so P1=P5*(1-deltaQngung*nuyTngung) T1=T5-((Tout-Tin)/(1-deltaQngung)) %nut so P2=P1*tonen T2=T1*to %chu trinh brayton G=Qin/(Cp*(T3-T2)) %nhiet luong tai tung nut chu trinh Q1=G*Cp*T1 Q2=G*Cp*T2 HVTH: Trần Thành Niệm 84 GVHD: PGS.TS Lê Chí Kiên Luận văn Q3=G*Cp*T3 Q4=G*Cp*T4 Q5=G*Cp*T5 % cong suat may phat MHD PMHD=Q3*nEE % cac cong suat tieu thu chu trinh hesobom=1.25;%thong so cai dat chu trinh Ptron=4000;%thong so cai dat chu trinh Ptach=2000;%thong so cai dat chu trinh Pbom=hesobom*(T2-T1)*G*Cp Pbomdt=250;% tra cuu %hieu suat he thong LMMHD nuytong=(PMHDPtron-Ptach-Pbom-Pbomdt)/Qin %tinh entropy Tref=298; Pref=1.03e+005; %chu trinh MHD S=(Cp*log(T3/Tref)-Cp*((y-1)/y)*log(P3/Pref))*G S3=S+(Cp*log(T3/T2)-Cp*((y-1)/y)*log(P3/P2))*G S4=S3+(Cp*log(T4/T3)-Cp*((y-1)/y)*log(P4/P3))*G S5=S4+(Cp*log(T5/T4)-Cp*((y-1)/y)*log(P5/P4))*G S1=S5+(Cp*log(T1/T5)-Cp*((y-1)/y)*log(P1/P5))*G S2=S1+(Cp*log(T2/T1)-Cp*((y-1)/y)*log(P2/P1))*G x=[S3 S4 S5 S1 S2 S3]; y=[T3 T4 T5 T1 T2 T3]; h=plot(x,y,'b-*'); grid xlabel('Truc Entropy') ylabel('Truc Nhiet Do (K)') - Chƣơng trình tính tốn dịng điện, hiệu suất máy phá AC MHD theo trƣờng hợp clc clear all %thong so cai dat ban dau J=31.8*1.00e+006; S=3.14*1.00e-006; L=0.154; muy=4*pi*1.00e-007; sigma1=1.00e+007; HVTH: Trần Thành Niệm 85 GVHD: PGS.TS Lê Chí Kiên Luận văn sigma2=1.00e+007; sigma3=1.00e+006; f=50; s1=-2; s2=-4; s3=-4.5; lamda=0.154; c=0.004; U=7.7; %thong so gian tiep I=J*S H=I/L Rmf1=muy*sigma1*U*lamda/(2*pi) Rmf2=muy*sigma2*U*lamda/(2*pi) Rmf3=muy*sigma3*U*lamda/(2*pi) omega=2*pi*f B=muy*H %ket qua ngo I1=B*c*lamda*sigma1*s1*U/2 I2=B*c*lamda*sigma2*s2*U/2 I3=B*c*lamda*sigma3*s3*U/2 nuy1=1/(1-s1) nuy2=1/(1-s2) nuy3=1/(1-s3) x=linspace(0,0.08,500) y1=I1*cos(omega*x)+(I/s1^2)*cos(s1*omega*x) y2=I2*cos(omega*x)+(I/s2^2)*cos(s2*omega*x) y3=I3*cos(omega*x)+(I/s3^2)*cos(s3*omega*x) plot(x,y1,'r',x,y2,'b',x,y3,'g') grid xlabel('Thoi gian (s)') ylabel('Dong dien (A)') - Chƣơng trình tính tốn hiệu suất so sánh dạng sóng dịng điện ngõ máy phát AC MHD đƣợc chọn với trƣờng hợp clc clear all %thong so cai dat ban dau J=31.8*1.00e+006; S=3.14*1.00e-006; HVTH: Trần Thành Niệm 86 GVHD: PGS.TS Lê Chí Kiên Luận văn L=0.154; muy=4*pi*1.00e-007; sigma1=0.5*1.00e+007; sigma2=1.00e+007; sigma3=1.00e+006; sigma4=1.00e+007; f=50; s1=-2; s2=-4; s3=-4.5; s4=-0.6; lamda=0.154; c=0.004; U=7.7; %thong so gian tiep I=J*S H=I/L Rmf1=muy*sigma1*U*lamda/(2*pi) Rmf2=muy*sigma2*U*lamda/(2*pi) Rmf3=muy*sigma3*U*lamda/(2*pi) Rmf4=muy*sigma4*U*lamda/(2*pi) omega=2*pi*f B=muy*H %ket qua ngo I1=B*c*lamda*sigma1*s1*U/2 I2=B*c*lamda*sigma2*s2*U/2 I3=B*c*lamda*sigma3*s3*U/2 I4=B*c*lamda*sigma4*s4*U/2 nuy1=1/(1-s1) nuy2=1/(1-s2) nuy3=1/(1-s3) nuy4=1/(1-s4) x=linspace(0,0.08,500) y1=I1*cos(omega*x)+(I/s1^2)*cos(s1*omega*x) y2=I2*cos(omega*x)+(I/s2^2)*cos(s2*omega*x) y3=I3*cos(omega*x)+(I/s3^2)*cos(s3*omega*x) y4=I4*cos(omega*x)+(I/s4^2)*cos(s4*omega*x) plot(x,y1,'r',x,y2,'b',x,y3,'g',x,y4,'y') grid xlabel('Thoi gian (s)') HVTH: Trần Thành Niệm 87 GVHD: PGS.TS Lê Chí Kiên Luận văn ylabel('Dong dien (A)') - Chƣơng trình tính tốn hiệu suất chu trình địa nhiệt truyền thống clc clear all %thong so cai dat ban dau P1=10;%KPa P2=400;%KPa P3=P2; h3=3273.41;%KJ/Kg s3=7.8984;%KJ/Kg.K h1=191.8;%KJ/Kg v1=0.00101;%m^3/Kg Wp=v1*(P2-P1) h2=h1+Wp Qh=h3-h2 gamaSteam=0.6493; deltaTuabin=7.5009; x4=(s3-gamaSteam)/deltaTuabin; gamaCondenser=2392.8; h4=h1+x4*gamaCondenser %nang luong sinh WT=h3-h4 Wnet=WT-Wp %hieu suat he thong nuyHT=Wnet/Qh %nhiet cac nut T3=673; toTuabin=0.68; T4=T3*toTuabin Nbom=1.07; T1=333; T2=T1*Nbom %do thi Ts x=[h3 h4 h1 h2 h3] y=[T3 T4 T1 T2 T3] h=plot(x,y,'b-*') grid xlabel('Truc Entropy') HVTH: Trần Thành Niệm 88 GVHD: PGS.TS Lê Chí Kiên Luận văn ylabel('Truc nhiet (K)') - Chƣơng trình tính tốn hiệu suất chu trình địa nhiệt sủ dụng tuabin cấp clc clear all %thong so cai dat ban dau P1=10;%KPa P6=P1 P2=400;%KPa P3=P2; P4=200;%KPa P5=P4 h3=3273.41;%KJ/Kg s3=7.8984;%KJ/Kg.K h1=191.8;%KJ/Kg v1=0.00101;%m^3/Kg h4s=504.68; gamaSteam1=1.7766; deltaTuabin1=5.1193; gamaCondenser1=2133.8; h5=3276.55;%KJ/Kg s5=8.2217;%KJ/Kg.K x4=(s3gamaSteam1)/deltaTuabin1 h4=h4s+x4*gamaCondenser1 Whp=h3-h4 gamaSteam2=0.6493; deltaTuabin2=7.5009; x6=(s5gamaSteam2)/deltaTuabin2 gamaCondenser2=2392.8; h6=h1+x6*gamaCondenser2 Wlp=h5-h6 Wt=(h3-h4)+(h5-h6) Wp=v1*(P2-P1) h2=h1+Wp Qh=(h3h2)+(h5-h4) %nang luong sinh Wnet=Wt-Wp %hieu suat he thong nuyHT=Wnet/Qh HVTH: Trần Thành Niệm 89 GVHD: PGS.TS Lê Chí Kiên Luận văn %nhiet cac nut T3=673; toTuabin=0.68; T4=T3*toTuabin T5=T3 Nbom=1.07; toCondenser=0.56; T1=333; T6=T5*toCondenser T2=T1*Nbom %do thi Ts x=[h3 h4 h5 h6 h1 h2 h3] y=[T3 T4 T5 T6 T1 T2 T3] h=plot(x,y,'b-*') grid xlabel('Truc Entropy') ylabel('Truc nhiet (K)') - Chƣơng trình tính tốn hiệu suất chu trình địa nhiệt – AC MHD clc clear all y=110;%he so dan nhiet cua kim loai long Qin=3.50e+003;%nhiet luong he thong dia nhiet cung cap P2=4.0e+005; Ttraodoinhiet=673;%nhiet bo trao doi nhiet T2=Ttraodoinhiet Ptron=0.2e+004; ntach=0.9; nEE=0.625; nACMHD=0.8; deltaQACMHD=0.005; deltaQngung=0.01; nuyTngung=1; tonen=2; to=0.75; Cp=0.27;%nhiet dung rieng NaK deltaTtron=5; deltaTtach=5; %nut T3=T2deltaTtron HVTH: Trần Thành Niệm 90 GVHD: PGS.TS Lê Chí Kiên Luận văn P3=Ptron %nut T4=T3*(1-nEE-deltaQACMHD) P4=P3*((1-(nEE/nACMHD))^(y/(y-1))) %nut P5=P4/2 P5a=P5 T5=T4-deltaTtach %nut Tout=680; Tin=630; P1=P5*(1-deltaQngung*nuyTngung) T1=T5-((Tout-Tin)/(1-deltaQngung)) %nut P2=P1*tonen T2=T1*to %luu luong kim loai long G=Qin/(Cp*(T3-T2)) %nhiet luong tai tung nut cua chu trinh Q1=G*Cp*T1 Q2=G*Cp*T2 Q3=G*Cp*T3 Q4=G*Cp*T4 Q5=G*Cp*T5 %cong suat may phat AC MHD WACMHD=Q3*nEE %cac cong suat tieu thu chu trinh hesobom=1.25;%thong so cai dat chu trinh Wtron=1.0e+003;%thong so cai dat chu trinh Wtach=0.5e+003;%thong so cai dat chu trinh Wbom=hesobom*0.02*G*Cp Wbomdt=150;%thong so tra cuu %hieu suat he thong ACMHD Wnet2=WACMHD-Wtron-Wtach-Wbom-Wbomdt nuytong=Wnet2/Qin %tinh Entropy Tref=298; Pref=1.03e+005; %chu trinh ACMHD h=(Cp*log(T3/Tref)-Cp*((y-1)/y)*log(P3/Pref))*G HVTH: Trần Thành Niệm 91 GVHD: PGS.TS Lê Chí Kiên Luận văn h3=h+(Cp*log(T3/T2)-Cp*((y-1)/y)*log(P3/P2))*G h4=h3+(Cp*log(T4/T3)-Cp*((y-1)/y)*log(P4/P3))*G h5=h4+(Cp*log(T5/T4)-Cp*((y-1)/y)*log(P5/P4))*G h1=h5+(Cp*log(T1/T5)-Cp*((y-1)/y)*log(P1/P5))*G h2=h1+(Cp*log(T2/T1)-Cp*((y-1)/y)*log(P2/P1))*G x=[h3 h4 h5 h1 h2 h3]; y=[T3 T4 T5 T1 T2 T3]; h=plot(x,y,'b-*'); grid xlabel('Truc Entropy') ylabel('Truc nhiet (K)') HVTH: Trần Thành Niệm 92 GVHD: PGS.TS Lê Chí Kiên ... động, địa nhiệt nhiệt mặt trời" nghiên cứu với mục tiêu tính tốn so sánh hiệu suất phát điện hiệu kinh tế hai mô hình phát điện kết hợp máy phát điện Từ thủy động với địa nhiệt nhiệt mặt trời. .. nghệ phát điện máy phát MHD có bước phát triển Nó ứng dụng rộng rãi chu trình phát điện hỗn hợp để nâng cao hiệu suất nhà máy điện Đề tài: "Tính tốn hiệu hệ thống phát điện kết hợp Từ thủy động,. .. DỤC VÀ ĐÀO TẠO TRƢỜNG ĐẠI HỌC SƢ PHẠM KỸ  THUẬT THÀNH PHỐ HỒ CHÍ MINH - - LUẬN VĂN THẠC SĨ TRẦN THÀNH NIỆM TÍNH TỐN HIỆU QUẢ CỦA HỆ THỐNG PHÁT ĐIỆN KẾT HỢP TỪ THỦY ĐỘNG, ĐỊA NHIỆT VÀ NHIỆT MẶT