Tính toán phân bố dòng có điều khiển với việc sử dụng bộ điều chỉnh góc pha áp dụng vào hệ thống điện miền bắc

Cao hc 2006 -2008 giáo dục đào tạo trờng đại học bách khoa hà nội - luận văn thạc sĩ khoa học tính toán phân bố dòng có điều khiển với việc sử dụng điều chỉnh góc pha áp dụng vào hệ thống điện miền bắc ngành : hệ thống điện m số:23.04.3898 vũ ngọc trung Ngời hớng dẫn khoa học : TS đỗ xuân khôi Hà Nội 2008 _ Hc viờn: V Ngc Trung Cao hc 2006 -2008 MC LC DANH MC CC CH VIT TT: HT - h thng in CXL - ch xỏc lp PST - Mỏy bin ỏp iu chnh PS - B iu chnh gúc pha FACTS - Flexible AC Transmission Systems TCSC - Thyristor Controlled Series Capacitor TCPST - Thyristor-Controlled Phase Shifting Transformer SSSC - Static Synchronous Series Compentation UPFC - Unified Power Flow Controller DANH MC BNG BIU: Bng 1: Cỏc dng thit b FACTS 13 Bng 3.1 phn ph lc Bng 3.2 phn ph lc Bng 3.3: kt qu tớnh toỏn in ỏp nỳt (mụun, gúc pha) ca HT 220-500kV Bc nm 2009 66 Bng 3.4: kt qu tớnh toỏn dũng in v cụng sut trờn cỏc nhỏnh ng dõy 67 Bng 3.5: kt qu tớnh toỏn in ỏp nỳt (mụun, gúc pha) 70 Bng 3.6: l kt qu tớnh toỏn dũng in v cụng sut trờn cỏc nhỏnh ng dõy 71 Bng 5.1: Kt qu tớnh toỏn lng cụng sut tỏc dng trờn ng dõy t PS thay i t 200 (MW) ữ 450 (MW) .92 Bng 5.2: Kt qu tớnh toỏn in ỏp (moun v gúc pha) t PS trờn ng dõy Thng Tớn H ụng v cụng sut tỏc dng truyn vo ng dõy ny l 450 (MW) .94 Bng 5.3: Kt qu tớnh toỏn phõn b dũng v cụng sut t PS trờn ng dõy Thng Tớn H ụng v cụng sut tỏc dng truyn vo ng dõy ny l 450 (MW) .95 Bng 5.4: Kt qu tớnh toỏn lng cụng sut tỏc dng trờn ng dõy t PS thay i t 200 (MW) ữ 450 (MW) .98 Bng 5.5: Kt qu tớnh toỏn in ỏp (moun v gúc pha) t PS trờn ng dõy Thng Tớn H ụng v cụng sut tỏc dng truyn vo ng dõy ny l 450 (MW) .100 _ Hc viờn: V Ngc Trung Cao hc 2006 -2008 Bng 5.6: Kt qu tớnh toỏn phõn b dũng, v cụng sut t PS trờn ng dõy Thng Tớn H ụng v cụng sut tỏc dng truyn vo ng dõy ny l 450 (MW) .100 Bng 5.7: Kt qu tớnh toỏn in ỏp (moun v gúc pha) t PS trờn ng dõy Ho Bỡnh Nho Quan v cụng sut tỏc dng truyn vo ng dõy ny l 900 (MW) .103 Bng 5.8: Kt qu tớnh toỏn phõn b dũng, v cụng sut t PS trờn ng dõy Ho Bỡnh Nho Quan v cụng sut tỏc dng truyn vo ng dõy ny l 900 (MW) .104 Bng 5.9: Kt qu tớnh toỏn in ỏp (moun v gúc pha) t PS trờn ng dõy Hũa Bỡnh Nho Quan v cụng sut tỏc dng truyn vo ng dõy ny l 900 (MW) .106 Bng 5.10: Kt qu tớnh toỏn phõn b dũng, v cụng sut t PS trờn ng dõy Hũa Bỡnh Nho Quan v cụng sut tỏc dng truyn vo ng dõy ny l 900 (MW) .107 DANH MC HèNH V: Hỡnh 1.1: iu khin cụng sut trờn ng dõy truyn ti in xoay chiu 10 Hỡnh 1.2: Cỏc dng thit b FACTS iu khin dũng cụng sut .13 Hỡnh 1.3: Cỏc dng thit b FACTS iu khin dũng cụng sut tỏc dng 14 Hỡnh 1.4: Thit b FACTS iu khin c dũng cụng sut tỏc dng v .16 Hỡnh 1.5: S nguyờn lý ca B iu chnh gúc pha 19 Hỡnh 1.6: B iu chnh gúc pha ni gia h thng in .21 Hỡnh 1.7: th vộc-t v c tuyn cụng sut ph thuc vo 22 Hỡnh 1.8: S hai ngun in ni vi 23 Hỡnh 1.9: Hai h thng in ni vi ng dõy song song 24 Hỡnh 1.10: th vộc t in ỏp khụng cn bng 24 Hỡnh 1.11: c tớnh n nh quỏ ca h thng ng c n gin, 26 Hỡnh 1.12: Gim dao ng cụng sut bi B iu chnh gúc pha 28 Hỡnh 2.1: S thut toỏn Newton-Raphson tớnh toỏn ch xỏc lp 39 Hỡnh 2.2: S hm MaiNewto-Raphson 41 Hỡnh 2.3: S thut toỏn tớnh ma trn tng dn .46 Hỡnh 2.4: S thut toỏn b xung nỳt cú Shunt vo ma trn tng dn 46 Hỡnh 2.5: S thut toỏn kim tra gii hn 49 Hỡnh 2.6: S tớnh toỏn dũng cụng sut gi, nhn, 51 Hỡnh 3.1: Bn HT 220-500kV Bc 52 Hỡnh 3.2: S nguyờn lý HT 220-500kV Bc 53 Hỡnh 3.3: S thay th ng dõy 60 Hỡnh 3.4: S thay th mỏy bin ỏp t ngu 61 Hỡnh 3.5: S thay th Khỏng bự ngang v T bự dc 62 Hỡnh 3.6: Chiu cụng sut v t l truyn ti theo t l phỏt núng ca ng dõy 69 Hỡnh 4.1: S mỏy bin bin ỏp iu chnh 77 Hỡnh 4.2: S thay th mỏy bin ỏp iu chnh .78 _ Hc viờn: V Ngc Trung Cao hc 2006 -2008 Hỡnh 4.3: Mụ hỡnh nỳt ca nhỏnh mỏy bin ỏp iu chnh .79 LI NểI U CHNG 1:PHNG PHP IU KHIN DềNG CễNG SUT GII THIU V MT S THIT B FACTS IU KHIN DềNG CễNG SUT 1.1 Phng phỏp iu khin dũng cụng sut HT: 1.1.1 t : 1.1.2 Mi quan h gia cỏc thụng s HT c iu khin trờn ng dõy ti in: 1.2 Gii thiu v mt s thit b FACTS iu chnh dũng cụng sut dũng cụng sut HT: 13 1.3 Cỏc gii hn v kh nng truyn ti cụng sut HT: .17 1.4 Nguyờn lý hot ng v tỏc dng ca B iu chnh gúc pha: 18 1.4.1 S nguyờn lý hot ng ca B iu chnh gúc pha: 18 1.4.2 Tỏc dng ca B iu chnh gúc pha: .19 1.4.3 Túm tt chc nng yờu cu ca B iu chnh gúc pha: 28 CHNG 2:H PHNG TRèNH TNH TON CXL V PHNG PHP NEWTON-RAPHSON 30 2.1 H phng trỡnh tớnh toỏn CXL 30 2.2 Phng phỏp Newton-Raphson 33 2.2.1 Phng trỡnh lp Newton-Raphson 33 2.2.2 Jacobian ca cỏc phng trỡnh ch xỏc lp 34 2.3 Tỡm hiu v xõy dng chng trỡnh tớnh toỏn phõn b CXL bng phng phỏp lp Newton-Raphson 36 2.3.1 c im phng phỏp Newton-Raphson chun v ỏp dng cho bi toỏn gii tớch li in .36 2.3.2 Tỡm hiu v phỏt trin chng trỡnh tớnh toỏn phõn b dũng bng phng phỏp N-R MatLab 40 CHNG 3: TNH TON PHN B DềNG CXL CA H THNG IN 220-500KV MIN BC 53 3.1 c im ca ch xỏc lp ca HT: .53 3.2 S hnh v nhng c im ca h thng in 220-500kV Bc hin ti .53 3.2.1 Hin trng ngun in Bc: .54 3.2.2 Hin trng li in 220-500kV: 57 3.3 Mụ hỡnh hoỏ cỏc phn t HT li in 220-500kV Bc 60 3.3.1 Mụ hỡnh hoỏ ng dõy: .60 3.3.2 Mụ hỡnh hoỏ mỏy bin ỏp t ngu: .60 3.3.3 Mụ hỡnh hoỏ khỏng bự ngang v T bự doc: 61 _ Hc viờn: V Ngc Trung Cao hc 2006 -2008 3.4 p dng chng trỡnh Newton Raphson chun tớnh toỏn phõn b dũng CXL cho li in Bc: 63 3.4.1 Tớnh toỏn h n v tng i: 63 3.4.2 S liu nỳt v nhỏnh: 64 3.4.3 Tớnh toỏn phõn b dũng ch hnh bỡnh thng (ch xỏc lp) ca h thng li in 220-500kV Bc nm 2009: 66 3.4.4 Tớnh toỏn phõn b dũng ch hnh ct mch 220kV Ho BỡnhH ụng ca h thng li in 220-500kV Bc: 70 CHNG 4: Mễ HèNH HểA B IU CHNH GểC PHA TRONG LI IN 77 4.1 Mỏy bin ỏp iu chnh (PST): 77 4.1.1 Phng trỡnh mụ t PST h thng li in: 78 4.2 Chng trỡnh tớnh toỏn phõn b dũng bng phng phỏp Newton-Raphson cú B iu chnh gúc pha (PS) Matlab: 85 4.2.1 Hm MaiPSNewton-Raphson (Chng trỡnh chớnh): .85 4.2.2 Hm PSNewtonRaphson (thc hin vũng lp): 88 4.2.3 Hm PSCalculatedPowers: 88 4.2.4 Hm PSPowerMismatches: .88 4.2.5 Hm PSNewtonRaphsonJacobian: .89 4.2.6 Hm PSUpdates: 89 4.2.7 PSLimits: 89 4.2.8 Hm PSPQflows: 89 CHNG 5: HIU QU CA B IU CHNH GểC PHA NNG CAO KH NNG TRUYN TI CễNG SUT TRONGH THNG IN 220500KV MIN BC 91 5.1 t : 91 5.2 Bin phỏp nõng cao kh nng truyn ti cụng sut trờn HT 220-500kV Bc bng cỏch t B iu chnh gúc pha (PS): 91 5.3 Phng ỏn t PS trờn ng dõy 220kV Thng Tớn H ụng c ch : ch hnh bỡnh thng v ch ct mch ng dõy Ho bỡnh H ụng: 92 5.3.1 Ch hnh bỡnh thng: 92 5.3.2 Ch ct mch ng dõy Ho Bỡnh H ụng: 98 5.4 Phng ỏn t PS trờn ng dõy 220kV Ho Bỡnh Nho Quan ch : ch hnh bỡnh thng v ch ct mch ng dõy Ho bỡnh H ụng: 103 5.4.1 Ch hnh bỡnh thng: 103 5.4.2 Ch ct mch ng dõy Ho Bỡnh H ụng: 106 5.5 Kt lun: 109 KT LUN V KIN NGH 111 TI LIU THAM KHO 113 _ Hc viờn: V Ngc Trung Cao hc 2006 -2008 PH LC LI NểI U Nhu cu nng lng ton cu núi chung v Vit Nam núi riờng ngy cng mt gia tng ú nng lng in úng vai trũ quan trng ỏp ng nhu cu trờn h thng in (HT) cng ngy cng phỏt trin v m rng, nhiu ng dõy truyn ti in di in ỏp cao v siờu cao c hỡnh thnh liờn kt cỏc h thng in khu vc vi hỡnh thnh HT ln, vi mc ớch truyn ti trao i cụng sut gia cỏc khu vc nhm m bo hnh kinh t v nõng cao tin cy cho HT Theo quy hoch ca h thng in Vit Nam n nm 2010 thỡ li in quc gia ó ph n hu ht cỏc vựng c nc v ang hỡnh thnh liờn kt h thng li in Quc gia vi h thng li in cỏc nc khu vc Nn kinh t nc ang phỏt trin n nh t tc cao vi vic ang hỡnh thnh cỏc khu, cm cụng nghip, vựng kinh t trung iu ú ó lm cho s gia tng ph ti ngy cng cao v nú ó tr thnh ỏp lc i vi ngnh in, cn phi xõy dng thờm nhiu nh mỏy v ng dõy ti in m ngun u t cho cỏc cụng trỡnh mi ca ngnh in cũn rt hn ch, ú hnh HT ngy cng phi m bo cỏc yờu cu v cung cp in liờn tc, an ton, tin cy v n nh ũi hi rt cao phỏt trin nn kinh t Mt khỏc h thng truyn ti in hin cng ó hỡnh thnh nhng mch vũng liờn kt gia cỏc nh mỏy in v ph ti, nhng quỏ trỡnh hnh thng xuyờn xut hin mt s ng dõy ti in b quỏ ti ú mt s ng dõy ti in khỏc luụn tỡnh trng non ti Cựng vi s phỏt trin vt bc lnh vc iu khin t ng, c bit lnh vc in t cụng sut vi vic ch to cỏc Thyristor cú cụng sut ln v k thut o lng iu khin h thng in (SCADA, DCS,) iu ú ó cho phộp to cỏc thit b iu khin linh hot (cỏc thit b FACTS) cỏc thụng s hnh h thng in xoay chiu _ Hc viờn: V Ngc Trung Cao hc 2006 -2008 Xut phỏt t nhng yờu cu trờn ti Tớnh toỏn phõn b dũng cú iu khin vi vic s dng B iu chnh gúc pha p dng vo h thng in Bc khụng ngoi mc ớch hnh HT hiu qu v kinh t - iu khin mt cỏch hp lý dũng cụng sut trờn HT cú th trỏnh c tỡnh trng mt s ng dõy truyn ti in hnh thng xuyờn b quỏ ti mt s ng dõy truyn ti in khỏc luụn tỡnh trng non ti, t ú gim c cỏc s c trờn cỏc ng dõy truyn ti, gim dao ng ca h thng, tng tỏc dng liờn kt ca li truyn ti - Nõng cao n nh in ỏp, n nh tnh cng nh n ng ca ton h thng - Gim lng cụng sut d tr ton h thng, h thng ln liờn kt vi cú kh nng huy ng ti u cụng sut t nhiu ngun phỏt T ú gim c chi phớ u t vo cỏc cụng trỡnh ngun - Gim chi phớ hnh, nõng cao tin cy cung cp in v kh hnh linh hot ca h thng in cụng sut chung ca h thng in ln - Tng tớnh kinh t chung v hiu qu hnh ca h thng in tn dng c cỏc ngun phỏt cú giỏ thnh xn sut in nng thp Ngoi phn m u, kt lun, ni dung chớnh ca lun c trỡnh by chng Chng 1: Phng phỏp iu khin dũng cụng sut - Gii thiu v mt s thit b FACTS iu khin dũng cụng sut Chng 2: H phng trỡnh ch xỏc lp v phng phỏp NewtonRaphson Chng 3: Tớnh toỏn phõn b dũng ch xỏc lp ca h thng in 220500kV Bc Chng 4: Mụ hỡnh hoỏ B iu chnh gúc pha li in Chng 5: Hiu qu ca B iu chnh gúc pha nõng cao kh nng truyn ti cụng sut h thng in 220-500kV Bc _ Hc viờn: V Ngc Trung Cao hc 2006 -2008 CHNG PHNG PHP IU KHIN DềNG CễNG SUT GII THIU V MT S THIT B FACTS IU KHIN DềNG CễNG SUT 1.1 Phng phỏp iu khin dũng cụng sut HT: 1.1.1 t : Cựng vi s phỏt trin rt nhanh ca HT, ũi hi nhng cụng ngh mi c ỏp dng khai thỏc trit cỏc kh nng sn cú ca h thng m khụng nh hng n an ton truyn ti v phõn phi in Cỏc nghiờn cu v h thng truyn ti in xoay chiu linh hot FACTS (Flexible AC Transmission Systems) ó c bt u vo cui nm 1980 nhm gii quyt chớnh: tng kh nng truyn ti ca cỏc h thng in v m bo cụng sut truyn ti khong gii hn thc hin cỏc yờu cu trờn ũi hi phi cú cỏc phng phỏp iu khin hiu qu dũng cụng sut Cụng ngh FACTS l nht cú kh nng iu khin nhanh mt cỏch linh hot dũng cụng sut tỏc dng v phn khỏng 1.1.2 Mi quan h gia cỏc thụng s HT c iu khin trờn ng dõy ti in: xem xột n s thay i ca cỏc thụng s ch ca HT liờn quan n iu khin dũng cụng sut, ta nghiờn cu vớ d n gin v ng dõy truyn ti trờn hỡnh 1.1(a) Nỳt v l cỏi ca cỏc nỳt ngun ln v n gin c gi thit l cỏc nỳt cú cụng sut vụ cựng ln Hai nỳt c ni vi thụng qua ng dõy truyn ti ch cú in khỏng X (b qua in tr v dung dn ca ng dõy) E1 v E2 l biờn in ỏp ti cỏc nỳt tng ng v gúc lch in ỏp gia chỳng l EL l vộc t in ỏp giỏng trờn in khỏng X ca ng dõy Biờn dũng in trờn ng dõy nhn c theo cụng thc: I = EL/ X, v lch vi in ỏp EL gúc 900 ( Hỡnh 1.1(b)) Dũng cụng sut chy trờn ng dõy cú th c iu khin bng vic iu khin in ỏp E1, E2, EL, X hoc Tht vy: _ Hc viờn: V Ngc Trung Cao hc 2006 -2008 10 Hỡnh 1.1(b), 1.1(c), ch mi quan h theo th vộc-t gia cỏc dũng in tỏc dng v phn khỏng vi in ỏp ti hai im nỳt P, Q E1 a) EL E1 E2 EL b) c) I (E1 - E2 cos) E1 E1 sin E2 cos E1 sin Ip1 = E2 sin/X E2 sin Ip2 = E1 sin/X E1 cos (E2 - E1cos) Iq1= (E1 - E2 cos)/X Iq2= (E2 - E1 cos)/X E2 Pmax E1E2 P = sin X d) E1 E1 - E2 I e) 90 180 Điện áp bơm I E2 Điện áp bơm E1 f) I E1 I E1 - E2 E2 E1 - E2 I g) I E2 Hỡnh 1.1: iu khin cụng sut trờn ng dõy truyn ti in xoay chiu (a) S h thng in n gin (b) th vộc t dũng in vuụng gúc vi in ỏp ri trờn ng dõy (c) th vộc t liờn quan n cụng sut tỏc dng v phn khỏng (d) Cỏc ng cong c tớnh ca cụng sut ph thuc vo gúc lch i vi in khỏng X khỏc (e) iu chnh biờn in ỏp lm thay i cụng sut phn khỏng (f) in ỏp bm vo vuụng gúc vi dũng in ng dõy lm thay i cụng sut tỏc dng (g) Ngun in ỏp bm vo dc theo in ỏp nỳt ca ng dõy Ti nỳt 1: - Thnh phn tỏc dng ca dũng in l: Ip1 = (E2 sin)/X _ Hc viờn: V Ngc Trung Cao hc 2006 -2008 11 ú P1 = E1Ip1 = E1(E2 sin)/X - Thnh phn phn khỏng ca dũng in l: Iq1 = (E1 E2 cos)/X ú Q1 = E1Iq1 = E1(E1 E2 cos)/X (1.1) Ti nỳt 2: - Thnh phn tỏc dng ca dũng in l: Ip2 = (E1 sin)/X ú P2 = E2Ip2 = E2(E2 sin)/X - Thnh phn phn khỏng ca dũng in l: Iq2 = (E2 E1 cos)/X ú Q2 = E2Iq2 = E2(E2 E1 cos)/X (1.2) Vỡ ó gi thit b qua in tr v dung dn ca ng dõy cho nờn khụng cú tn tht cụng sut tỏc dng trờn ng dõy Do ú: P1 = P2 = P = E1(E2 sin)/X (1.3) Theo cỏc biu thc (1.1), (1.2), (1.3) thay i cỏc giỏ tr E1, E2, EL = E1 E2, X hoc cỏc giỏ tr P, Q1, Q2 s thay i theo thy mc thay i ca P, Q1, Q2 thay i tng cỏc thụng s trờn, ta xem xột cỏc sau õy: - Gi thit rng, E1 v E2 l sut in ng bờn ca ngun quy i mụ t h thng v in khỏng X bao gm c in khỏng ca ngun quy i, hỡnh 1.1(d) ch mt na ng cong hỡnh sin cụng sut tỏc dng tng ti cc i gúc tng ti 90o Cụng sut sau ú gim tng ng vi gúc tng v cụng sut gim n = 180o D dng nhn thy cụng sut truyn ti trờn ng dõy t giỏ tr ln nht vi gúc 90o - Hỡnh 1.1(d) cng ch rng ta thay i giỏ tr X thỡ biờn ca ng cong cụng sut cng thay i, t ú nõng cao c gii hn cụng sut truyn ti trờn ng dõy - Cụng sut (dũng in) cú th cng c iu khin bi vic iu chnh in ỏp E1 hoc E2 Tuy nhiờn, trờn hỡnh 1.1(e) cho thy vi s thay i ca E1 thỡ hiu in ỏp E1 E2 cng khụng thay i nhiu, nhng gúc pha ca nú thỡ thay i nhiu iu ú cú ngha rng, vic thay i din ỏp E1 hoc E2 cú nh hng nhiu n cụng sut phn khỏng hn l _ Hc viờn: V Ngc Trung Cao hc 2006 -2008 116 PH LC 2: Bng 3.2: S liu ngun v ph ti li in 220-500kV Bc Nỳt nbb 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 Tờn nỳt hoabinh220 xuanmai hadong chem nhoquan220 ninhbinh bache namdinh thaibinh nghison hungdong hatinh220 donghoi maidong thuongtin220 phonoi phailai donghoa vatcach trangbach hoanbo uongbi sondong sosson tuyenquang bacgiang bacninh thainguyen tdtuyenquang ndhaiphong dinhvu tdbanla doluong ndcampha quangninh220 Um Pphỏt (kV) (MW) 230 230 330 230 230 230 230 230 230 230 230 230 230 230 230 230 230 230 230 230 230 230 230 230 230 230 230 230 242 230 230 242 230 230 230 1920 0 0 0 0 0 0 0 830 0 0 320 110 0 0 228 300 300 0 300 Qphỏt Pti Qti Pmin Qmax (MVAr) (MW) (MVAr) (MW) (MVAr) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 84 249 445 436 173 160 184 110 105 160 80 94 425 230 234 236 170 115 137 82 135 125 100 143 197 70 65 50 110 90 30 92 141 143 57 54 60 35 34 51 26 28 127 76 74 78 58 37 45 27 44 41 32 47 64 23 21 16 36 29 -1040 0 0 0 0 0 0 0 -140 0 0 -40 -15 0 0 -120 -80 -120 0 -40 1280 0 0 0 0 0 0 0 560 0 0 160 80 0 0 40 150 240 180 0 120 _ Hc viờn: V Ngc Trung Cao hc 2006 -2008 36 37 38 39 40 41 42 43 44 45 hoabinh500 thuongtin500 quangninh500 nhoquan500 hatinh500 tuhatinh500 tunhoquan500 nhoquanhatinhdoan1 nhoquanhatinhdoan2 nhoquanhatinhdoan3 117 500 500 515 500 500 500 500 500 500 500 0 600 0 0 0 0 0 0 0 0 0.25 0.125 0.125 0.125 0.125 0 0 4.008 2.004 2.004 2.004 2.004 0 0 0 -80 0 0 0 0 240 0 0 0 _ Hc viờn: V Ngc Trung Cao hc 2006 -2008 PH LC 3: 118 Cỏc hm chng trỡnh Newton-Raphson Hm MaiNewtonRaphson (Chng trỡnh chớnh) %Chuong trinh chinh NewtonRaphson chuan stop=0; while ~stop choice = menu(' Choose onething for you ', 'Nhap du lieu cho chuong trinh', 'Thanh lap ma tran tong dan cua luoi', 'Thuc hien vong lap de tim ket qua', 'Tinh cong suat nhanh', 'Tinh dong dien nhanh', 'Ket qua cuoi cung', 'Thoat khoi chuong trinh '); switch choice case delay = input('An Enter de nhap so lieu cho luoi'); PowerFlowsData5b1 delay = input('An mot phim bat ky de chay tiep '); case delay = input('An Enter de tinh ma tran tong dan cho luoi'); [YR,YI]=YBus(tlsend,tlrec,tlresis,tlreac,tlsuscep,tlcond,shbus,shresis,sh reac,ntl,nbb,nsh) delay = input('An mot phim bat ky de chay tiep '); case delay = input('An Enter de tinh dien ap cac nut cho luoi'); [VM,VA,it]=Newtonraphson(nmax,tol,itmax,ngn,nld,nbb,bustype,genbus,loadbu s,PGEN,QGEN,QMAX,QMIN,PLOAD,QLOAD,YR,YI,VM,VA) delay = input('An mot phim bat ky de chay tiep '); case delay = input('An Enter de tinh cong suat chay tren cac nhanh cho luoi'); [PQsend,PQrec,PQloss,PQbus,PQGEN]=PQflows(nbb,ngn,ntl,nld,genbus,loadbus, tlsend,tlrec,tlresis,tlreac,tlcond,tlsuscep,PLOAD,QLOAD,VM,VA,nsh,shbus,s hreac,shresis) delay = input('An mot phim bat ky de chay tiep '); case delay = input('An Enter de tinh cong dien chay tren cac nhanh cho luoi'); [CURRENTsend,CURRENTrec]=CURRENTflows(nbb,ngn,ntl,nld,genbus,loadbus,tlse nd,tlrec,tlresis,tlreac,tlcond,tlsuscep,PLOAD,QLOAD,VM,VA,nsh,shbus,shrea c,shresis) delay = input('An mot phim bat ky de chay tiep '); case delay = input('An Enter de tinh ket qua cuoi cung'); it VM VA=VA*180/pi PQsend=PQsend*100 PQrec=PQrec*100 CURRENTsend=(CURRENTsend*100)/(sqrt(3)*230) _ Hc viờn: V Ngc Trung Cao hc 2006 -2008 119 CURRENTrec=(CURRENTrec*100)/(sqrt(3)*230) delay = input('An mot phim bat ky de chay tiep '); case stop=1; msgbox(' Thanks for using my program '); end end -2 Hm YBus (Thnh lp ma trn tng dn) %xay dung ma tran tong dan Y _ function [YR,YI]=YBus(tlsend,tlrec,tlresis,tlreac,tlsuscep,tlcond,shbus,shresis,sh reac,ntl,nbb,nsh) YR = zeros(nbb,nbb); YI = zeros(nbb,nbb); for kk=1:ntl ii=tlsend(kk); jj=tlrec(kk); denom=tlresis(kk)^2+tlreac(kk)^2; YR(ii,ii)=YR(ii,ii) + tlresis(kk)/denom + 0.5*tlcond(kk); YI(ii,ii)=YI(ii,ii) - tlreac(kk)/denom + 0.5*tlsuscep(kk); YR(ii,jj)=YR(ii,jj) - tlresis(kk)/denom; YI(ii,jj)=YI(ii,jj) + tlreac(kk)/denom; YR(jj,ii)=YR(jj,ii) - tlresis(kk)/denom; YI(jj,ii)=YI(jj,ii) + tlreac(kk)/denom; YR(jj,jj)=YR(jj,jj) + tlresis(kk)/denom + 0.5*tlcond(kk); YI(jj,jj)=YI(jj,jj) - tlreac(kk)/denom + 0.5*tlsuscep(kk); end for kk=1:nsh ii=shbus(kk); denom=shresis(kk)^2+shreac(kk)^2; YR(ii,ii)=YR(ii,ii) + shresis(kk)/denom; YI(ii,ii)=YI(ii,ii) - shreac(kk)/denom; end %ket thuc viec lap ma tran tong dan Y! -3 Hm NewtonRaphson (Thc hin vũng lp) %VONG LAP CUA PHUONG PHAP NEWTONRAPHSON function [VM,VA,it]=Newtonraphson(nmax,tol,itmax,ngn,nld,nbb,bustype,genbus,loadbu s,PGEN,QGEN,QMAX,QMIN,PLOAD,QLOAD,YR,YI,VM,VA) D=zeros(1,nmax); flag=0; it=1; %tinh cong suat dat tai cac nut [PNET,QNET]=NetPowers(nbb,ngn,nld,genbus,loadbus,PGEN,QGEN,PLOAD,QLOAD); while (it