Ứng dụng phương pháp nhánh và cận để quy hoạch mở rộng lưới điện khi xét đến mạng điện có sẵn

B GIO DC V O TO TRNG I HC BCH KHOA H NI VN SNG NG DNG PHNG PHP NHNH V CN QUY HOCH M RNG LI IN KHI XẫT N MNG IN Cể SN Chuyờn ngnh: H Thng in LUN VN THC S KHOA HC NGI HNG DN KHOA HC PGS.TS NGUYN LN TRNG H Ni 2010 LI CAM OAN Tụi xin cam oan lun ny l cụng trỡnh nghiờn cu ca riờng tụi, c hp t nhiu ngun ti liu, s liu khỏc T ú dng nhng kin thc ó c hc v tip thu t thc t hon thnh cụng trỡnh ny, khụng chộp ca bt k lun no trc ú Tụi xin cam hon ton chu trỏch nhim v ni dung ca lun ny H Ni, ngy 29 thỏng 10 nm 2010 VN SNG Khúa: CH 2008 - 2010 Lời cảm ơn Tác giả xin chân thành cảm ơn PGS.TS Nguyễn Lân Tráng có gợi mở dẫn dắt tận tình, cung cấp tài liệu quý giá để tác giả hoàn thành luận văn Tác giả xin chân thành cảm ơn thầy cô giáo môn Hệ thống điện - Đại học Bách Khoa Hà Nội ngời giúp đỡ hoàn thành luận văn Do kiến thức hạn chế nên luận văn khó tránh khỏi sai sót, tác giả mong nhận đợc bảo, góp ý thầy cô giáo môn Hệ thống điện ngời quan tâm MC LC M U Chng TNG QUAN V QUY HOCH H THNG IN 1.1 Mụ hỡnh tuyn tớnh tnh 1.2 Mụ hỡnh tuyn tớnh ng .11 1.3 Mụ hỡnh phc 16 1.4 Tng kt chng 18 Chng PHNG PHP NHNH V CN XC NH CU TRC TI U CA LI IN 19 2.1 t bi toỏn 19 2.2 Phng phỏp gii bi toỏn bng thut toỏn nhỏnh cn 20 2.3 Chng trỡnh nhỏnh v cn xỏc nh cu trỳc ti u ca li in 34 2.4 Tng kt chng 41 Chng P DNG PHNG PHP NHNH V CN TRONG TRNG HP TN TI MT MNG IN Cể SN 42 3.1 Hn ch ca phng phỏp nhỏnh v cn .42 3.2 Phng phỏp ỏp dng thut toỏn nhỏnh v cn vo xõy dng phng ỏn quy hoch li in ti u cú tớnh n mng in cú sn 42 3.3 Chng trỡnh nhỏnh v cn ci biờn xỏc nh cu trỳc ti u ca li in xột n mng in cú sn 48 3.4 Tng kt chng 57 Chng - NG DNG PHNG PHP NHNH V CN QUY HOCH M RNG LI IN 220KV THNH PH H NI KHI Cể XẫT N MNG IN Cể SN 58 4.1 Phõn tớch hin trng ngun, li in Thnh ph H Ni 58 4.2 c im chung v phng hng phỏt trin kinh t-xó hi Thnh ph H Ni giai on 2010-2015-2020 .68 4.3 D bỏo nhu cu ph ti Thnh ph H Ni n nm 2015, 2020 71 4.4 Phng ỏn quy hoch li in Thnh ph H Ni xột n mng in cú sn 77 4.5 Tng kt chng 95 KT LUN 96 TI LIU THAM KHO PH LC NG DNG PHNG PHP NHNH CN QUY HOCH M RNG LI IN KHI XẫT N MNG IN Cể SN PHN M U Lý chn ti Trong nhng nm qua cựng vi s phỏt trin nhanh chúng ca nn kinh t xó hi nhu cu tiờu th in cng khụng ngng tng trng vi tc cao Do vy m nhiu nh mỏy in v khu cụng nghờp ang c xõy dng v sp sa i Vic xõy dng cỏc ng dõy in a in nng t cỏc nh mỏy cung cp cho ph ti ũi hi phi cú u t rt ln Do vy t l phi cú mt phng phỏp thớch ỏng xỏc nh cu hỡnh mng li in da trờn mng in ó tn ti cho chi phớ tớnh toỏn l nh nht Hin nay, thụng thng quy hoch m rng li in, phng phỏp lit kờ l phng phỏp thng c s dng Phng phỏp ny da da trờn c s phõn tớch trc quan lit kờ mt s phng ỏn kh thi, nú cú quan h cht ch vi trỡnh v kinh nghim ca cỏc chuyờn gia a cỏc s kh thi Sau ú cn so sỏnh kinh t k thut tỡm phng ỏn ti u Tuy vy nu h thng thit k l bao gm rt nhiu nỳt thỡ s phng ỏn m chuyờn gia a phõn tớch l rt ln Mt khỏc nhng phng ỏn a phõn tớch hon ton l trc quan vy cú th b qua nhng phng ỏn hp lý khc phc nhc im ú, ngi ta s dng cỏc phng phỏp quy hoch toỏn hc v phng phỏp c s dng hin l phng phỏp nhỏnh v cn Tuy nhiờn, phng phỏp ny li cú nhc im ln l ch xột cho mt h thng mi hon ton v coi kh nng ti ca cỏc ng dõy l vụ tn iu ú ụi khụng thớch hp cho vic m rng li in thc t Phng phỏp c gii thiu sau õy ó khc phc c nhc im ú ca phng phỏp nhỏnh cn Nú cho phộp xỏc nh cu trỳc ti u ca li in c cú s phỏt trin ca mt li in ó tn ti v cho kt qu khỏ tin cy Mc ớch, i tng v phm vi nghiờn cu + Mc ớch nghiờn cu ng dng mt phng phỏp quy hoch toỏn hc cht ch vo thc tin, c th l phng phỏp ỏp dng k thut nhỏnh cn, cn c trờn hin trng thc t, nh Ngi thc hin: Vn Sỏng NG DNG PHNG PHP NHNH CN QUY HOCH M RNG LI IN KHI XẫT N MNG IN Cể SN hng phỏt trin kinh t xó hi, a mt mụ hỡnh cu trỳc ti u cho li in phự hp vi nhu cu ph ti giai on tng lai + i tng v phm vi nghiờn cu i tng nghiờn cu: i tng nghiờn cu chớnh ca lun l h thng truyn ti in cp in ỏp 220 kV cung cp in cho cỏc trm bin ỏp 110/35/22 kV ti cỏc trung tõm ph ti Phm vi nghiờn cu: Trong khuụn kh ca lun vn, chỳng tụi ch cp n mng in khu vc, c th l mng in cung cp cho Thnh ph H Ni Kho sỏt hin trng li in ang hot ng, thu thp s liu v nh hng phỏt trin ca tnh a nhu cu phỏt trin ph ti in thi gian ti p dng thut toỏn nhỏnh v cn xõy dng phng phỏp quy hoch m rng li in cú xột n mng in hin cú trờn c s phng phỏp nhỏnh v cn ó c nhiu nh khoa hc nghiờn cu Cỏc lun im c bn v úng gúp ca ti - Khc phc c nhc im ca ca phng phỏp nhỏnh cn ú l coi mng in l mi hon ton cú th ng dng vo quy hoch li in thc t - Xõy dng chng trỡnh mu cú th chy trờn mỏy tớnh lp quy hoch li in trờn c s lý thuyt mụ hỡnh toỏn quy hoch ti u Phng phỏp nghiờn cu Vic la chn phng phỏp nghiờn cu l mc ớch v i tng nghiờn cu quyt nh lun ny, chỳng tụi s dng hai phng phỏp chung theo s phõn loi ca khoa hc lun: phng phỏp nghiờn cu lý thuyt v phng phỏp nghiờn cu phi thc nghim Trong cỏc phng phỏp chung ú, nhng phng phỏp nghiờn cu c th c phi hp s dng: phng phỏp thng kờ, phng phỏp phõn tớch tng hp, phng phỏp d bỏo nhu cu, phng phỏp mụ hỡnh toỏn hc ti u, phng phỏp cn v nhỏnh Trong ú, phng phỏp cn v nhỏnh c s Ngi thc hin: Vn Sỏng NG DNG PHNG PHP NHNH CN QUY HOCH M RNG LI IN KHI XẫT N MNG IN Cể SN dng kt hp vi phng phỏp phõn tớch tng hp giỳp cho nhng nhn xột, kt lun cú c s khoa hc v tớnh thuyt phc cao B cc ca lun ti gm 96 trang Ngoi phn M u (03 trang), Kt lun (01 trang), ni dung c th hin chng (92 trang) Chng I Tng quan v quy hoch li in Chng II Phng phỏp nhỏnh v cn xỏc nh cu trỳc ti u ca li in Chng III p dng phng phỏp nhỏnh v cn trng hp ó tn ti mt mng in cú sn Chng IV ng dng phng phỏp nhỏnh v cn quy hoch m rng li in 220 kV Thnh ph H Ni xột n mng in ó cú sn Cui cựng l Phn ph lc v Danh mc ti liu tham kho Ngi thc hin: Vn Sỏng NG DNG PHNG PHP NHNH CN Chơng Tổng quan quy hoạch lới điện Bi toỏn quy hoch ti u h thng in bao gm tớnh toỏn nhm tho mc tiờu cc tiu chi phớ t trang thit b v cc tiu chi phớ hnh cỏc phn t li in, ng thi m bo tuõn th cỏc nh mc k thut i vi cht lng v tin cy cung cp in cỏc giai on phỏt trin s h thng in iu kin ti u c thit lp dng cc tiu tng cỏc chi phớ tớnh toỏn xột trờn tt c cỏc phn t li in cú xột n s bin i cỏc chi phớ ny theo tiờu chun kinh t k thut Khi so sỏnh cỏc phng ỏn khỏc xõy dng s li in thỡ cú th gi thit rng thnh phn chi phớ c nh i vi cỏc cụng trỡnh hin cú l bng i vi cỏc phng ỏn, v cú th loi b hm mc tiờu Thc hin ti u húa v tit din dõy dn truyn ti in sau chn c phng ỏn s phỏt trin li in Tựy theo cụng sut ti trờn ng dõy bin i vựng kinh t - k thut no m chn c s l ti u, hm chi phớ ti u tựy thuc vo s mch song song ca ng dõy Cỏc s liu ban u xut phỏt ban u ca bi toỏn quy hoch m rng li in thng khụng y v s lng cng nh cht lng Bi vy, cú th chp nhn cỏc phng ỏn gn ti u theo cỏc cỏch tip cn gii quyt bi toỏn nh sau: Bin i bi toỏn ti u húa h thng in v dng cho phộp ỏp dng c mt phng phỏp toỏn ti u cú gii thut cht ch, cho phộp dn n mt dng ỏp s ti u; cú th b qua mt s iu kin rng buc, nh th cho phộp gii c bi toỏn mt cỏch cht ch v mt toỏn hc, nhiờn nh vy ta ch nhn c li gii gn ti u v mt thc t Cũn cỏch gii quyt bi toỏn quy hoch phỏt trin m rng h thng in theo tỡnh ng ca cỏc giai on phỏt trin kinh t, xột theo mt quỏ trỡnh nhiu giai on tớnh toỏn liờn tip lp li thỡ s thu c mt gii phỏp tt nht iu kin cho phộp phỏt trin h thng theo tng giai on phi c tha trờn mi giai on phỏt trin h thng in, phi m bo cỏc quy phm k thut v l mt h thng m cú th phỏt trin chuyn bin t mt trng thỏi s ny Ngi thc hin: Vn Sỏng NG DNG PHNG PHP NHNH CN sang mt trng thỏi s khỏc Trong trng hp u t thc hin cho phng ỏn phỏt trin cỏc thi im khỏc thỡ chi phớ bin i theo thi gian, so sỏnh cỏc phng ỏn khỏc thỡ cỏc chi phớ mi thi im khỏc phi c quy dn v cựng mt thi im Bài toán quy hoạch phát triển mở rộng hệ thống điện thờng đợc đặt nh sau: Đã biết cấu trúc số liệu trạng lới điện vận hành, dự báo tốc độ phát triển phụ tải nút hệ thống, phụ tải phát triển khu vực suốt thời gian tính toán, cho trớc đờng dây, trạm biến áp đợc xây dựng tiêu kinh tế kỹ thuật có Yêu cầu tìm phơng án kết cấu lới tối u để xây dựng đa đối tợng vào vận hành, đảm bảo đạt đợc yêu cầu hàm chi phí tính toán cực tiểu giai đoạn thời gian T, đồng thời thoả mãn ràng buộc kinh tế kỹ thuật Nhìn chung, toán quy hoạch phát triển mở rộng lới điện toán tối u đa mục tiêu, nhiều cực trị với biến số liên tục rời rạc, ràng buộc toán có dạng đẳng thức bất đẳng thức Tuỳ theo lợng thông tin ban đầu có đợc, ngời ta thờng chia toán quy hoạch thành loại: tất định, ngẫu nghiên xác định bất định Đối với loại toán ngời ta có mô hình toán học phơng pháp giải khác Đối với toán tất định, thông tin ban đầu đợc xác định Căn vào yếu tố thời gian, lại đợc phân thành toán tĩnh, không xét đến yếu tố thay đổi theo thời gian toán động có xét đến thay đổi theo thời gian tham số Thông thờng, toán đợc chia thành dạng: mô hình tĩnh mô hình động 1.1 Bài toán tuyến tính tĩnh Hàm mục tiêu toán đợc biểu diễn nh sau: Z= m n n n i =1 j=1 k =1 j=1 Ci j x i j + C kj x kj Ngi thc hin: Vn Sỏng NG DNG PHNG PHP NHNH CN ij kj Trong : Cij l sut chi phớ tớnh toỏn ca ng dõy truyn ti t nỳt i n nỳt j (/MW) x ij l lng cụng sut truyn ti trờn ng dõy i j m: số nút nguồn cung cấp n: số nút phụ tải xkj l lng cụng sut truyn ti trờn t nỳt k ti nỳt j Vi cỏc rng buc v cõn bng cụng sut ch ph ti cỏc nỳt: x i j ij x kj = Pj j= 1,n k j Trong ú Pj l cụng sut ca nỳt ph ti th j (MW) Mụ hỡnh trờn l mụ hỡnh n gin ch xột ti rng buc v cụng sut Song khú khn õy l xỏc nh cỏc h s sut chi phớ tớnh toỏn hm mc tiờu vỡ quan h gia chi phớ tớnh toỏn v cụng sut truyn ti trờn ng dõy l quan h phi tuyn ph thuc vo tit din dõy dn a c bi toỏn v mụ hỡnh tuyn tuyn tớnh tnh ta phi tuyn tớnh hoỏ mi quan h ny Cho n cú rt nhiu cụng trỡnh nghiờn cu ca cỏc tỏc gi ó t bi toỏn quy hoch phỏt trin li in v xột thờm rt nhiu rng buc hn v kinh t k thut Trong s cỏc mụ hỡnh tnh ỏng lu ý nht l cỏc phng phỏp trỡnh by sau õy 1.1.2 Phng phỏp xõy dng v la chn cỏc phng ỏn cựng h Cú th thnh lp mt s ni dõy ban u no ú lm cn c xut phỏt, xõy dng cỏc phng ỏn cựng h vi nú S ni dõy ban u ny gm cú cỏc ng dõy vi chiu di ó bit, ni cỏc nỳt vi Ngoi cn phi d kin Ngi thc hin: Vn Sỏng {Tinh z2 cua phuong an minz1} assignfile(f3,tf6);reset(f3); read(f3,m,m1,m2,n,gz,gh); for i:=1 to (m+1)*(n+1) read(f3,s[i]); closefile(f3); (* goi thu tuc don hinh *) dhct(m,n,n1,gz,s,ib,x,y,val,it); Zm:=W1+val*Cp2; Z2tu:=val*Cp2; Ztu:=Zm; { } {Tinh W2} assignfile(f3,tf7);reset(f3); read(f3,m,m1,m2,n,gz,gh); for i:=1 to (m+1)*(n+1) read(f3,s[i]); closefile(f3); (* goi thu tuc don hinh *) dhct(m,n,n1,gz,s,ib,x,y,val,it); W2:=val*Cp2; {=================================================} {Thuat toan nhanh can tim phuong an toi uu} {=================================================} m:=nut;Tc:=m; for i:=1 to m for j:=1 to m C2[i,j]:=C[i,j];Sc:=0; { -} if checkbox1.checked=true then begin button1.visible:=false;checkbox1.visible:=false;button2.visible:=true;label5.visible:=true; label6.visible:=true; tf3:='ketquacuoi.txt'; assignfile(f3,tf3);rewrite(f3); writeln(f3,'Phuong an co so la phuong an bao gom cac canh: '); for i:=1 to m-1 for j:=i+1 to m if (Ctu[i,j]0) then writeln(f3,' Canh ',i,' - ',j,' dai ',Ctu[i,j]:5:2,' km'); writeln(f3,'Gia tri ham muc tieu: Z = ',Ztu:2:2); writeln(f3,' Z1 = ',Z1tu:2:2); writeln(f3,' Z2 = ',Z2tu:2:2); closefile(f3); memo1.lines.loadfromfile('ketquacuoi.txt'); memo1.visible:=true;buoc:=1;goto label2; end; { } {tim cay bao trum moi} for i:=1 to m for j:=1 to m begin Ak[i]:=0;At[i]:=i;C[i,j]:=0; end; r:=ln[1,2]; for i:=1 to m-1 for j:=i+1 to m if (r>ln[i,j])and(ln[i,j]0) then r:=ln[i,j]; for i:=1 to m for j:=1 to m if (r=ln[i,j])and(r0) then begin Ak[i]:=i;Ak[j]:=j;At[i]:=0;At[j]:=0;C[i,j]:=ln[i,j]; end; for k:=2 to m-1 begin for i:=1 to m-1 for j:=i+1 to m h[i,j]:=0; for i:=1 to m for j:=1 to m begin if (Ak[i]>0) and (At[j]>0) then begin if At[j]>Ak[i] then h[Ak[i],At[j]]:=ln[Ak[i],At[j]] else h[At[j],Ak[i]]:=ln[At[j],Ak[i]]; end;end; for i:=1 to m-1 for j:=i+1 to m if h[i,j]>0 then begin r:=h[i,j];break; end; for i:=1 to m-1 for j:=i+1 to m if (r>h[i,j]) and (h[i,j]>0) then r:=h[i,j]; for i:=1 to m for j:=1 to m if (r=ln[i,j])and(r0) then begin Ak[i]:=i;Ak[j]:=j;At[i]:=0;At[j]:=0;C[i,j]:=ln[i,j];end; end; Z1:=0; for i:=1 to m-1 for j:= i+1 to m Z1:=Z1+C[i,j];Z1:=Z1*Cp1; if Z1+W2Zm*) end; if Sc=nut then begin tf3:='ketquacuoi.txt'; assignfile(f3,tf3);rewrite(f3); writeln(f3,'Phuong an toi uu la phuong an bao gom cac canh: '); for i:=1 to nut-1 for j:=i+1 to nut if (Ctu[i,j]0) then writeln(f3,' Canh ',i,' - ',j,' dai ',Ctu[i,j]:5:2,' km'); writeln(f3,'Gia tri ham muc tieu: Z = ',Ztu:2:2); writeln(f3,' Z1 = ',Z1tu:2:2); writeln(f3,' Z2 = ',Z2tu:2:2); writeln(f3);closefile(f3); end; memo1.lines.loadfromfile('ketquacuoi.txt'); memo1.visible:=true; label2:end; end Ph lc Mó ngun chng trỡnh nhỏnh v cn ci biờn unit QHLD; uses Windows, Messages, SysUtils, Classes, Graphics, Controls, Forms, Dialogs,StdCtrls; const mmax=200;nmax=500; type mang1=array[1 (mmax+1)*(nmax+1)]of real;mang2=array[1 nmax+mmax]of real; mang3=array[1 mmax]of integer;mang4=array[1 mmax,1 mmax]of real; mang5=array[1 mmax,1 mmax]of integer; TForm1 = class(TForm) GroupBox1: TGroupBox;Label1: TLabel;Label2: TLabel;Label5: TLabel;Edit1: TEdit;Edit2: TEdit; Edit3: TEdit;Button1: TButton;Label3: TLabel;Edit4: TEdit;Edit5: TEdit;Label4: TLabel; Memo1: TMemo;Label6: TLabel;Label7: TLabel;Button2: TButton;Label8: TLabel; CheckBox1: TCheckBox; procedure Button1Click(Sender: TObject); procedure Button2Click(Sender: TObject); private { Private declarations } public { Public declarations } end; var s:mang1;st,l,h,lm,ln,lc,T1,T2,ls:mang4;x,y,id:mang2;ib:mang3; m,mc,m1,m2,n,n1,i,j,it,u,gh,k,Sc,Tc,nut,v,v1,tt,buoc,dem:integer; g,gz,val,r,W1,W2,Z1,Z2,Zm,Z1tu,Z2tu,Ztu,Z,dienap:real; f1,f2,f3,f4,f5:text;tf1,tf2,tf3,tf4,tf5,tf6,tf7:string[30]; At,Ak:mang3;C,C1,C2,s1,Ctu:mang4; Pg,Pt,dongtai,cos:mang2;phu1,phu2,phu3:string[100]; phu4:string[7];phu5:string[20]; Form1: TForm1; implementation {$R *.DFM} procedure dhct(m,n,n1:integer;gz:real; s:mang1;ib:mang3;var x,y:mang2; var val:real; var it:integer); label l1,l2; var dk,dk1,te,r:real; i,j,mf,nf,fj,ff,ij,k,ik:integer; fi,lk,l,lj,kj:integer; begin it:=0;mf:=m+1;nf:=n+1;ff:=mf*nf; for i:=1 to m ib[i]:=n+i;fj:=ff-nf; for j:=1 to nf begin fj:=fj+1;r:=0;ij:=j; for i:=1 to m begin r:=r+s[ij];ij:=ij+nf; end; s[fj]:=gz*r-s[fj]; end; (* Bat dau buoc lap moi -*) l1: it:=it+1;dk:=1.0e+30; (* Tim cot xoay va kiem tra toi uu *) l2: k:=1;fj:=ff-n;r:=s[fj]; for j:=2 to n begin fj:=fj+1; if (s[fj]>r)and(s[fj]>0.00001)and(s[fj]0) and (At[j]>0) then begin if At[j]>Ak[i] then h[Ak[i],At[j]]:=l[Ak[i],At[j]] else h[At[j],Ak[i]]:=l[At[j],Ak[i]];end;end; for i:=1 to m-1 for j:=i+1 to m if h[i,j]>0 then begin r:=h[i,j];break; end; for i:=1 to m-1 for j:=i+1 to m if (r>h[i,j]) and (h[i,j]>0) then r:=h[i,j]; for i:=1 to m for j:=1 to m if (r=l[i,j])and(r0) then begin Ak[i]:=i;Ak[j]:=j;At[i]:=0;At[j]:=0;C[i,j]:=l[i,j];end;end; W1:=0; for i:=1 to m-1 for j:= i+1 to m W1:=W1+C[i,j]; for i:=1 to m for j:=1 to m Ctu[i,j]:=C[i,j];Z1tu:=W1; {=========================================} {tim z2 cua phuong an co z1 va phuong an co z2} {==========================================} { } {du lieu tinh z2 cua phuong an minz1} tf6:='Z2.SLI'; assignfile(f3,tf6);rewrite(f3);k:=1; for i:=1 to m-1 for j:=i+1 to m if C[i,j]0 then begin C1[i,k]:=-1;C1[j,k]:=1;C1[i,k+m-1]:=1;C1[j,k+m-1]:=-1;C1[m,k]:=C[i,j];C1[m,k+m-1]:=C[i,j];k:=k+1; end; write(f3,m-1:10);write(f3,0:10);write(f3,0:10);write(f3,(m-1)*2:10);write(f3,1E+7:10);write(f3,1:10); writeln(f3); for i:=1 to m begin for j:=1 to (m-1)*2 write(f3,C1[i,j]:10:2);write(f3,dongtai[i]:10:2);writeln(f3);end; closefile(f3); { } { Du lieu tinh minz2} for i:=1 to (m-1) for j:=1 to m*(m-1) s1[i,j]:=0; for i:=1 to (m-1) for j:=((i-1)*(m-1)+1) to ((i-1)*(m-1)+(m-1)) s1[i,j]:=-1; for i:=1 to (m-1) for j:=1 to m begin if (s1[i,j*(m-1)+i]=0)and(j*(m-1)+i>((i-1)*(m-1)+(m-1))) then s1[i,j*(m-1)+i]:=1; if (i>1)and(s1[i,(j-1)*(m-1)+i-1]=0)and (((j-1)*(m-1)+i-1)0 then write(f5,l[i,j]:10:2);write(f5,0:10);writeln(f5); closefile(f5); { } {Tinh z2 cua phuong an minz1} assignfile(f3,tf6);reset(f3); read(f3,m,m1,m2,n,gz,gh); for i:=1 to (m+1)*(n+1) read(f3,s[i]); closefile(f3); (* goi thu tuc don hinh *) dhct(m,n,n1,gz,s,ib,x,y,val,it); Zm:=W1+val; Z2tu:=val; Ztu:=Zm; { } {Tinh W2} assignfile(f3,tf7);reset(f3); read(f3,m,m1,m2,n,gz,gh); for i:=1 to (m+1)*(n+1) read(f3,s[i]); closefile(f3); (* goi thu tuc don hinh *) dhct(m,n,n1,gz,s,ib,x,y,val,it); W2:=val; {=================================================} {Thuat toan nhanh can tim phuong an toi uu} {=================================================} { } {tim cay bao trum moi} for i:=1 to m for j:=1 to m begin Ak[i]:=0;At[i]:=i;C[i,j]:=0; end; r:=1000; for i:=1 to m-1 for j:=i+1 to m if (r>ln[i,j])and(ln[i,j]0) then r:=ln[i,j]; for i:=1 to m for j:=1 to m if (r=ln[i,j])and(r0) then begin Ak[i]:=i;Ak[j]:=j;At[i]:=0;At[j]:=0;C[i,j]:=ln[i,j]; end; for k:=2 to m-1 begin for i:=1 to m-1 for j:=i+1 to m h[i,j]:=0; for i:=1 to m for j:=1 to m begin if (Ak[i]>0) and (At[j]>0) then begin if At[j]>Ak[i] then h[Ak[i],At[j]]:=ln[Ak[i],At[j]] else h[At[j],Ak[i]]:=ln[At[j],Ak[i]]; end; end; for i:=1 to m-1 for j:=i+1 to m if h[i,j]>0 then begin r:=h[i,j]; break; end; for i:=1 to m-1 for j:=i+1 to m if (r>h[i,j]) and (h[i,j]>0) then r:=h[i,j]; for i:=1 to m for j:=1 to m if (r=ln[i,j])and(r0) then begin Ak[i]:=i;Ak[j]:=j;At[i]:=0;At[j]:=0;C[i,j]:=ln[i,j]; end; end; tf3:='NC.SLI'; assignfile(f3,tf3);rewrite(f3); for i:=1 to m begin for j:=1 to m write(f3,C[i,j]:7:2); writeln(f3); end; writeln(f3); Z1:=0; for i:=1 to m-1 for j:= i+1 to m Z1:=Z1+C[i,j]; writeln(f3,'Z1 = ',Z1:5:5);writeln(f3,Z1+w2:5:5);writeln(f3,sc); closefile(f3); if Z1+W2Zm*) end; if Sc=nut then begin tf3:='ketqua1.txt'; assignfile(f1,tf3);rewrite(f1); writeln(f1,'Phuong an toi uu chua xet den mang dien co san: '); for i:=1 to nut-1 for j:=i+1 to nut if (Ctu[i,j]0) then writeln(f1,' Canh ',i,' - ',j,' dai ',Ctu[i,j]:5:2,' km'); writeln(f1,'Gia tri ham muc tieu: Z = ',Ztu:1:5); writeln(f1,' Z1 = ',Z1tu:1:5); writeln(f1,' Z2 = ',Z2tu:1:5); writeln(f1); {closefile(f1); } end; {memo1.lines.loadfromfile('ketquacuoi.txt'); memo1.visible:=true; } {=============================================} {Nhanh can cai bien} {=============================================} {Xac dinh tap T1, T2} m:=nut;v:=0;v1:=0; for i:=1 to m-1 for j:=i+1 to m if (lc[i,j]0) then v:=v+1; for i:=1 to m for j:=1 to m begin if (Ctu[i,j]0) and (lc[i,j]0) then begin T1[i,j]:=lc[i,j]; end; if (Ctu[i,j]=0) and (lc[i,j]0) then begin T2[i,j]:=lc[i,j];end;end; for i:=1 to m-1 for j:=i+1 to m if (T2[i,j]0) then v1:=v1+1; {Tim cay bao trum nho nhat chua luoi dien cu} for i:=1 to m for j:=1 to m begin Ak[i]:=0;At[i]:=i; C[i,j]:=0;end; for i:=1 to m for j:=1 to m if lc[i,j]0 then begin Ak[i]:=i;Ak[j]:=j;At[i]:=0;At[j]:=0;C[i,j]:=l[i,j]; end; for k:=mc to m-1 begin for i:=1 to m-1 for j:=i+1 to m h[i,j]:=0; for i:=1 to m for j:=1 to m begin if (Ak[i]>0) and (At[j]>0) then begin if At[j]>Ak[i] then h[Ak[i],At[j]]:=l[Ak[i],At[j]] else h[At[j],Ak[i]]:=l[At[j],Ak[i]]; end; end; for i:=1 to m-1 for j:=i+1 to m if h[i,j]>0 then begin r:=h[i,j]; break; end; for i:=1 to m-1 for j:=i+1 to m if (r>h[i,j]) and (h[i,j]>0) then r:=h[i,j]; for i:=1 to m for j:=1 to m if (r=l[i,j])and(r0) then begin Ak[i]:=i;Ak[j]:=j;At[i]:=0;At[j]:=0;C[i,j]:=l[i,j];end;end; W1:=0; for i:=1 to m-1 for j:= i+1 to m if (C[i,j]0) and (lc[i,j]=0) then W1:=W1+C[i,j]; for i:=1 to m-1 for j:= i+1 to m if (C[i,j]0) and (T1[i,j]>0) then W1:=W1+C[i,j]; { } {Tinh Z2' cua phuong an D1'} (* du lieu tim Z2 cua cay bao trum moi*) tf3:='TG.SLI'; assignfile(f3,tf3);rewrite(f3); for i:=1 to m for j:=1 to (m-1)*2 C1[i,j]:=0;k:=1; for i:=1 to m-1 for j:=i+1 to m if C[i,j]0 then begin C1[i,k]:=-1;C1[j,k]:=1;C1[i,k+m-1]:=1;C1[j,k+m-1]:=-1;C1[m,k]:=C[i,j];C1[m,k+m-1]:=C[i,j];k:=k+1; end; write(f3,m-1:10);write(f3,0:10);write(f3,0:10);write(f3,(m-1)*2:10);write(f3,1E+7:10);write(f3,1:10); writeln(f3); for i:=1 to m begin for j:=1 to (m-1)*2 write(f3,C1[i,j]:10:2);write(f3,dongtai[i]:10:2);writeln(f3);end; closefile(f3); { } {Tinh z2 cua phuong an cay bao trum moi} assignfile(f3,tf3);reset(f3); read(f3,m,m1,m2,n,gz,gh); for i:=1 to (m+1)*(n+1) read(f3,s[i]); closefile(f3); (* goi thu tuc don hinh *) dhct(m,n,n1,gz,s,ib,x,y,val,it); Z2:=val;Zm:=W1+val; writeln(f1,'Phuong an cay bao trum nho nhat chua mang dien co san: '); for i:=1 to nut-1 for j:=i+1 to nut if (C[i,j]0) then writeln(f1,' Canh ',i,' - ',j,' dai ',C[i,j]:5:2,' km'); writeln(f1,'Gia tri ham muc tieu: Z = ',Zm:1:5); writeln(f1,' Z1 = ',W1:1:5); writeln(f1,' Z2 = ',val:1:5); writeln(f1); closefile(f1); if Ztu>Zm then begin Ztu:=Zm;Z1tu:=W1;Z2tu:=val; for i:=1 to nut for j:=1 to nut Ctu[i,j]:=C[i,j]; end else Zm:=Ztu; {=================================================} {Thuat toan nhanh can cai bien tim phuong an toi uu} {=================================================} m:=nut;Tc:=m; for i:=1 to m for j:=1 to m C2[i,j]:=C[i,j]; for i:=1 to nut for j:=1 to nut begin ls[i,j]:=l[i,j];lm[i,j]:=ls[i,j]; end; Sc:=0;buoc:=0; if checkbox1.checked=true then begin button1.visible:=false;checkbox1.visible:=false;button2.visible:=true; label6.visible:=true;label7.visible:=true; memo1.lines.loadfromfile('ketqua1.txt');memo1.visible:=true;goto laben3; end; label2:Sc:=Sc+1; if (Scr)and(lc[i,j]=0) then r:=C2[i,j]; for i:=1 to m for j:=1 to m if ln[i,j]=r then begin ln[i,j]:=0;C2[i,j]:=0; end; { } {tim cay bao trum moi} for i:=1 to m for j:=1 to m begin Ak[i]:=0;At[i]:=i;C[i,j]:=0; end; for i:=1 to m for j:=1 to m if lc[i,j]0 then begin Ak[i]:=i;Ak[j]:=j;At[i]:=0;At[j]:=0;C[i,j]:=l[i,j]; end; for k:=mc to m-1 begin for i:=1 to m-1 for j:=i+1 to m h[i,j]:=0; for i:=1 to m for j:=1 to m begin if (Ak[i]>0) and (At[j]>0) then begin if At[j]>Ak[i] then h[Ak[i],At[j]]:=ln[Ak[i],At[j]] else h[At[j],Ak[i]]:=ln[At[j],Ak[i]]; end; end; for i:=1 to m-1 for j:=i+1 to m if h[i,j]>0 then begin r:=h[i,j];break;end; for i:=1 to m-1 for j:=i+1 to m if (r>h[i,j]) and (h[i,j]>0) then r:=h[i,j]; for i:=1 to m for j:=1 to m if (r=ln[i,j])and(r0) then begin Ak[i]:=i;Ak[j]:=j;At[i]:=0;At[j]:=0;C[i,j]:=ln[i,j];end;end; Z1:=0; for i:=1 to m-1 for j:= i+1 to m if (C[i,j]>0)and(lc[i,j]=0) then Z1:=Z1+C[i,j]; for i:=1 to m-1 for j:= i+1 to m if (C[i,j]0) and (T1[i,j]>0) then Z1:=Z1+C[i,j]; (* -*) (* du lieu tim Z2 cua cay bao trum moi*) tf6:='Z2.SLI'; assignfile(f3,tf6);rewrite(f3); for i:=1 to m for j:=1 to (m-1)*2 C1[i,j]:=0;k:=1; for i:=1 to m-1 for j:=i+1 to m if C[i,j]0 then begin C1[i,k]:=-1;C1[j,k]:=1;C1[i,k+m-1]:=1;C1[j,k+m-1]:=-1;C1[m,k]:=C[i,j];C1[m,k+m-1]:=C[i,j];k:=k+1; end; write(f3,m-1:10);write(f3,0:10);write(f3,0:10);write(f3,(m-1)*2:10);write(f3,1E+7:10);write(f3,1:10); writeln(f3); for i:=1 to m begin for j:=1 to (m-1)*2 write(f3,C1[i,j]:10:2);write(f3,dongtai[i]:10:2);writeln(f3); end; closefile(f3); { } {Tinh z2 cua phuong an cay bao trum moi} assignfile(f3,tf6);reset(f3); read(f3,m,m1,m2,n,gz,gh); for i:=1 to (m+1)*(n+1) read(f3,s[i]); closefile(f3); (* goi thu tuc don hinh *) dhct(m,n,n1,gz,s,ib,x,y,val,it); Z2:=val; if Z1+Z2 Zm*) {end; (*ket thuc dieu kien Z1+W2=Zm then begin for i:=1 to nut for j:=1 to nut lm[i,j]:=ls[i,j];goto label1; end; } (*ket thuc dieu kien Z1+W2>Zm*) end; if (Sc=nut-1) then begin tf3:='ketqua1.txt'; assignfile(f3,tf3);rewrite(f3); writeln(f3,'Phuong an toi uu la phuong an bao gom cac canh: '); for i:=1 to nut-1 for j:=i+1 to nut if (Ctu[i,j]0) then writeln(f3,' Canh ',i,' - ',j,' dai ',Ctu[i,j]:5:2,' km'); writeln(f3,'Gia tri ham muc tieu: Z = ',Ztu:1:5); writeln(f3,' Z1 = ',Z1tu:1:5); writeln(f3,' Z2 = ',Z2tu:1:5); writeln(f3); closefile(f3); memo1.lines.loadfromfile('ketqua1.txt'); memo1.visible:=true;button1.visible:=false;label6.visible:=false;label7.visible:=false; checkbox1.visible:=false; end;laben3:end; end ... dụng phơng pháp đợc nhiều nhà khoa học nghiên cứu, phơng pháp nhánh cận Đây phơng pháp quy hoạch tối u có thuật toán rõ ràng, bớc chặt chẽ, lập trình để tạo phần mềm chuyên dụng 2.2 Phơng pháp giải... NHNH CN Chơng phơng pháp nhánh cận xác định cấu trúc tối u lới điện 2.1 Đặt toán Quy hoạch lới điện phần quan trọng quy hoạch hệ thống lợng Nhiệm vụ xác định cấu trúc tối u lới điện theo tăng trởng... sơ đồ quy hoạch nguồn thời gian quy hoạch ứng với yêu cầu phân phối điện cách an toàn kinh tế Nguyên lý quy hoạch cực tiểu cấu trúc lới chi phí vận hành nhằm thoả mãn yêu cầu phân phối điện an