Giải tích hệ thống phân phối hình tia với thiết bị bù tcvr dùng giải thuật lfb

ðại Học Quốc Gia Tp Hồ Chí Minh TRƯỜNG ðẠI HỌC BÁCH KHOA KHÚC THỊ KIM QUYÊN GIẢI TÍCH HỆ THỐNG PHÂN PHỐI HÌNH TIA VỚI THIẾT BỊ BÙ TCVR DÙNG GIẢI THUẬT LFB Chuyên ngành : Thiết bị, Mạng Nhà máy ñiện LUẬN VĂN THẠC SĨ TP HỒ CHÍ MINH, tháng 12 năm 2010 LỜI CẢM ƠN Em xin trân trọng bày tỏ lòng biết ơn sâu sắc lời cảm ơn chân thành ñến TS Hồ Văn Hiến, Thầy tận tình bảo, hướng dẫn tạo điều kiện tốt để em hồn thành luận văn thạc sĩ Xin chân thành cám ơn ñến tất quý Thầy, Cô Trường ðại Học Bách Khoa Tp.HCM ñã trang bị cho em khối lượng kiến thức bổ ích, đặc biệt xin chân thành biết ơn sâu sắc đến Thầy, Cơ Khoa ðiện ñã giúp ñỡ, tạo ñiều kiện thuận lợi hỗ trợ cho em nhiều thời gian học tập trường Cuối em xin gởi lời cám đến Gia đình Bạn bè lớp giúp ñỡ, ñộng viên, ủng hộ tạo cho em niềm tin để em nỗ lực cố gắng thực hoàn thành luận văn Do thời gian có hạn kiến thức cịn nhiều hạn chế nên luận văn chắn nhiều khiếm khuyết Em chân thành mong muốn nhận bảo, góp ý Thầy, Cơ bạn đọc quan tâm ñến nội dung luận văn Tp.Hồ Chí Minh, 06 / 12 / 2010 Khúc Thị Kim Quyên DANH MỤC CÁC HÌNH Hình Tên hình Trang Hình 1.1 Từ thơng mạch từ Hình 1.2 Từ thơng mạch từ Hình 1.3 Tiết đện dây dẫn hình trụ Hình 1.4 Từ thơng móc vịng bên ngồi Hình 1.5 Tiết đện dây dẫn hình trụ Hình 1.6 ðường dây pha hai dây 10 Hình 1.7 ðường dây ba pha bốn dây 11 Hình 1.8 ðường dây ba pha hốn vị đầy đủ 15 Hình 1.9 ðường dây ba pha lộ kép hốn vị ñầy ñủ 16 Hình 1.10 ðiện trường ñối với cực song song, điện tích điểm dây dẫn hình trụ 19 Hình 1.11 ðiện trường hai dây dẫn song song 21 Hình 1.12 ðiện dung đường dây ba pha đối xứng 23 Hình 1.13 ðiện dung ñường dây ba pha ñối xứng 24 Hình 1.14 ðiện dung pha đường dây tải điện 24 Hình 1.15 Ảnh hưởng ñất ñiện dung ñối với ñường dây ba pha 25 Hình 1.16 ðiện dung ñường dây lộ kép 26 Hình 1.17 Truyền tải ñiện ba pha 28 Hình 1.18 Truyền tải điện ba pha 28 Hình 1.19 Mạch tương đương đường dây ngắn 29 Hình 1.20 ðồ thị vectơ cho trường hợp tải có tính dung 29 Hình 2.1 Hệ thống phân phối hình tia 36 Hình 2.2 Mạch vịng thứ cấp 39 Hình 2.3 Phát tuyến nối vịng máy cắt thường đóng đường nối khẩn cấp 40 Hình 2.4 Phát tuyến nối mạch vịng 40 Hình 2.5 Những nét hệ thống mạng phân phối thứ cấp 42 Hình 2.6 Hệ thống mạng phân phối sơ cấp 44 Hình 2.7 Trạm biến áp phân phối đơn giản 46 Hình 2.8 Trạm biến áp phân phối ñơn giản với cách bố trí chiển mạch phía sơ cấp Hình 2.9 46 Trạm biến áp phân phối với máy cắt nối phân đoạn phía điện áp thấp 47 Hình 2.10 Trạm biến áp phân phối với góp đơi phía thứ cấp 48 Hình 3.1 Mơ hình chung thiết bị FACTS 50 Hình 3.2 Ứng dụng Facts thiết bị bù 51 Hình 3.3 Cấu trúc STATCOM 54 Hình 3.4 ðặc tính V-I STATCOM 54 Hình 3.5 Sơ đồ bù tĩnh SVC 55 Hình 3.6 ðặc tính V-I bù tĩnh SVC 56 Hình 3.7 GCSC 56 Hình 3.8 SSSC 58 Hình 3.9 TCSC 58 Hình 3.10 TSSC 59 Hình 3.11 TCSR 60 Hình 3.12 Mơ hình mạch điện TCVR với tải RL 60 Hình 3.13 Cấu trúc đầu phân áp tạo Thyristor 62 Hình 3.14 Cấu trúc ñầu phân áp ñược tạo Thyristor cho việc ñiều khiển liên tục ñiện áp ngõ 62 Hình 3.15 Hai điện áp đầu phân áp 64 Hình 3.16 UPEC 64 Hình 4.1 Mạng phân phối hình tia 13 nút, 12 dây 69 Hình 4.2 BFS đánh số hệ thống 13 nút 70 Hình 5.1 Bản vẽ quy hoạch cấp ñiện phân phối 22KV khu công nghiệp VSIPI 74 Hình 5.2 Sơ đồ ứng dụng mạng điện phân phối khu cơng nghiệp VSIPI 75 Hình 6.1 Sơ đồ mạng ñiện phân phối khu công nghiệp VSIPI sử dụng TCVR Hình 6.2 89 Biểu đồ điện áp trước sau gắn TCVR vào mạng lưới phân phối VSIPI 100 DANH MỤC CÁC TỪ VIẾT TẮT AC Alternating Current BFS Breadth – First – Search DC Direct Current FACTS Flexible Alternating Current Transmission System GCSC GTO Thyristor-Controller Series Capacitor GMR Bán kính trung bình hình học hay bán kính trung bình nhân GTO Gate Turn – Off LFB Line – Flow – Based SCR Semiconductor Controlled Rectifier SSSC Static Synchronous Series Compensator STATCOM Static Synchronous Compensator SVC Static Var Compensator TCR Thyristor Controlled Reator TCSC Thyristor Controlled Series Capacitor TCSR Thyristor Controlled Series Reactor TCVR Thyristor-Controlled Voltage Regulator TSC Thyristor Switched Capacitor TSR Thyristor Switched Reator TSSC Thyristor Switched Series Capacitor UPEC Unified Power Flow Controller VSIP I Vietnam Singapore Industrial Park I VSI Voltage Source Inverter MỤC LỤC MỞ ðẦU CHƯƠNG 1: LÝ THUYẾT VỀ ðƯỜNG DÂY TẢI ðIỆN 1.1 ðiện trở ñường dây tải ñiện 1.2 ðiện cảm ñường dây tải ñiện 1.2.1 Các hệ thức ñiện cảm 1.2.2 ðiện cảm bán kính trung bình nhân dây dẫn 1.2.3 ðiện cảm ñường dây pha 10 1.2.4 ðiện cảm ñường dây ba pha 11 1.2.5 ðường dây ñơn ba pha ñối xứng 14 1.2.6 ðường dây đơn ba pha hốn vị 14 1.2.7 ðường dây ba pha lộ kép 15 1.2.8 Dùng bảng ñể kiểm tra cảm kháng dây dẫn 17 1.3 ðiện dung ñường dây tải ñiện 19 1.3.1 ðiện trường cực song song, điện tích điểm dây dẫn hình trụ 19 1.3.2 ðiện dung ñường dây pha 21 1.3.3 ðiện dung ñường dây ba pha ñối xứng 23 1.3.4 ðiện dung ñường dây lộ kép 26 1.4 Mơ hình đường dây tải điện 27 1.4.1 Truyền tải ñiện ba pha 27 1.4.2 ðường dây truyền tải ngắn 29 1.4.3 Cách tính số mạch A , B , C , D ñường dây tải ñiện 30 CHƯƠNG 2: CẤU TRÚC CƠ BẢN CỦA HỆ THỐNG ðIỆN PHÂN PHỐI 34 2.1 Cấp phân phối hệ thống ñiện 34 2.2 Hệ thống phân phối hình tia 35 2.3 Hệ thống mạch vòng thứ cấp 38 2.4 Mạch vòng sơ cấp 39 2.5 Hệ thống phân phối mạng ñiện thứ cấp 41 2.6 Hệ thống mạng ñiện sơ cấp 44 2.7 Trạm biến áp phân phối 45 CHƯƠNG 3: CÁC THIẾT BỊ FACTS 50 3.1 Mô tả ðịnh Nghĩa 50 3.2 Cấu trúc 50 3.2.1 Thiết bị ñiều khiển nối tiếp 51 3.2.2 Thiết bị ñiều khiển song song 52 3.2.3 Thiết bị kết hợp ñiều khiển nối tiếp - nối tiếp 52 3.2.4 Thiết bị kết hợp ñiều khiển nối tiếp – song song 52 3.2.5 Các tính chất 53 3.3 Các ñiều khiển mắc song song 53 3.3.1 Static Synchronous Compensator (STATCOM) 53 3.3.2 Static Var Compensator (SVC) 54 3.4 Các ñiều khiển mắc nối tiếp 56 3.4.1 GTO Thyristor-Controller Series Capacitor (GCSC) 56 3.4.2 Static Synchronous Series Compensator (SSSC) 57 3.4.3 Thyristor Controlled Series Capacitor (TCSC) 58 3.4.4 Thyristor Switched Series Capacitor (TSSC) 59 3.4.5 Thyristor Controlled Series Reactor (TCSR) 60 3.4.6 Thyristor-Controlled Voltage Regulator (TCVR) 60 3.5 Thiết bị kết hợp ñiều khiển nối tiếp – song song UPEC 64 CHƯƠNG 4: THUẬT TOÁN LFB ( LINE – FLOW – BASED ) 66 4.1 Mơ hình dịng cơng suất phân phối hình tia LFB 66 4.1.1 Những phương trình cân cơng suất chung 66 4.1.2 Các phương trình điện áp nhánh 67 4.1.3 Những phương trình dịng cơng suất LFB 68 4.2 Bản chất Graph mạng lưới phân phối hình tia 69 4.3 Mơ hình LFB “phân lặp” hệ thống phân phối hình tia 71 CHƯƠNG 5: ỨNG DỤNG THUẬT TỐN LFB VÀO MẠNG LƯỚI PHÂN PHỐI HÌNH TIA THỰC TẾ 73 5.1 Ứng dụng thuật toán LFB vào mạng lưới phân phối 22KV hình tia thực tế 73 5.2 Ứng dụng thuật toán Newton Raphson vào mạng lưới phân phối hình tia 84 CHƯƠNG 6: ỨNG DỤNG THUẬT TOÁN LFB VÀO MẠNG LƯỚI PHÂN PHỐI HÌNH TIA THỰC TẾ CĨ SỬ DỤNG THIẾT BỊ BÙ TCVR87 6.1 Ứng dụng thuật toán LFB với thiết bị bù TCVR vào mạng lưới phân phối hình tia 87 6.2 Ứng dụng thuật toán LFB với thiết bị bù TCVR vào mạng lưới phân phối 22KV hình tia khu cơng nghiệp VSIPI 88 6.3 Ứng dụng thuật toán Newton Raphson vào mạng lưới phân phối hình tia có sử dụng thiết bị bù TCVR khu công nghiệp VSIPI 101 KẾT LUẬN 104 MỞ ðẦU ðặt vấn ñề Với phát triển điện tử cơng suất dùng thyristor tắt cổng (GTO – Gate turn – off ) ñã tạo tiền ñề cho ñời hệ thống truyền tải xoay chiều linh hoạt (Flexible Alternating Current transmission system – FACTS) Những lợi ích to lớn mà FACTS ñem lại mạng lưới truyền tải, nên ứng dụng phát triển nhiều nơi giới Trong xu hình thành phát triển thị trường điện giới, ngày nhiều cơng ty điện lực tiến hành đại hố lưới điện phân phối theo hướng áp dụng công nghệ FACTS Mục tiêu việc cài ñặt thiết bị FACTS vào mạng phân phối ñể cải thiện ñộ sụt áp, ñiều chỉnh hệ số cơng suất giảm tổn thất đường dây Với thiết bị FACTS tăng độ tin cậy hiệu hệ thống truyền tải phân phối Chúng cung cấp linh hoạt lớn điều khiển vận hành Mục đích luận văn phân tích mạng phân phối hình tia với thiết bị FACTS TCVR cài ñặt dùng giải thuật LFB ( Line – Flow – Based) Thuận lợi phương pháp ñộ lớn ñiện áp dịng cơng suất trên đường dây biến độc lập phương trình dịng cơng suất Vì dễ dàng điều khiển hoạt động thiết bị FACTS hệ thống phân phối cách trực tiếp dễ dàng thực thi Ngồi cịn sử dụng BFS (Breadth – First – Search), ma trận liên thuộc hệ thống phân phối hình tia, dẫn đến giảm sai số tính tốn Phạm vi nghiên cứu Luận văn trình bày: - Tổng quan lý thuyết ñường dây tải ñiện - Cấu trúc hệ thống phân phối - Mô tả số dạng FACTS mơ hình PHỤ LỤC A MATLAB MƠ PHỎNG THUẬT TỐN LFB ỨNG DỤNG VÀO MẠNG PHÂN PHỐI HÌNH TIA 22KV KHU CƠNG NGHIỆP VSIPI function [A,p,q,V] = LBF(data,P,Q,R,X) clc; data =[ 0.00409 0.0124 0.00218 0.00661 0.00491 0.01488 0.00273 0.00826 0.00164 0.00496 0.00205 0.0062 0.00436 0.01322 0.00341 0.01033 10 0.00191 0.00579 11 0.003 0.00909 12 0.00245 0.00744]; P = [-0.01 ; -0.006; -0.006; -0.002; -0.006; -0.008; -0.006; 0.0032; -0.01; -0.01; -0.016]; Q = [-0.0075; -0.0045; -0.0045; -0.0015; -0.0045; -0.006; -0.0045; 0.0024; -0.0075; -0.0075; -0.012]; M = size(data); G = M(1,1); R = zeros(G); X = zeros(G); l = zeros(G,1); m = zeros(G,1); k = zeros(G,1); RX = zeros(G); pq = zeros(G); A_c = zeros(G,1); A_c(1,1)=1; V_pv = zeros(G,1); V_pv(1,1) = 1; N = 0; for i = 1:G N = max(N,data(i,2)); end for i = 1:N for j = 1:(N - 1) A(i,j) = 0; end end for i = 1:G A(data(i,1),data(i,2)-1) = 1; A(data(i,2),data(i,2)-1) = -1; R(i,i) = data(i,3); X(i,i) = data(i,4); end A (1,:)= []; A_1 = A > A_2 = A < 0; A_2 = -1*A_2 B = ones (G,1); for i = 1: G RX(i,i) = R(i,i).^2 + X(i,i).^2; end p = A\(P - A_1*l); q = A\(Q - A_1*m); V_V = (A_1'+A_2')\(-A_c.*V_pv + k + 2*R*p + 2*X*q); V = V_V.^(1/2); for i = 2:K for i = 1:G pq(i,i) = p(i,1).^2 + q(i,1).^2; C(i,1) = V_V(i,1)\B(i,1); end l = pq*R*(-A_2')*C; m = pq*X*(-A_2')*C; k = pq*RX*(-A_2')*C; p = A\(P - A_1*l); q = A\(Q - A_1*m); V_V = (A_1'+A_2')\(-A_c.*V_pv + k + 2*R*p + 2*X*q); V = V_V.^(1/2); end V S = p+q*1i End PHỤ LỤC B MATLAB MÔ PHỎNG PHƯƠNG PHÁP NEWTON RAPHSON ỨNG DỤNG VÀO MẠNG PHÂN PHỐI HÌNH TIA 22KV KHU CƠNG NGHIỆP VSIPI function [ybus, power_temp,delta_power,ja,change_temp,u] newton_raphson (data,u,s,code) clc; data =[ 0.00409 + 0.0124i 0.00218 + 0.00661i = 0.00491 + 0.01488i 0.00273 + 0.00826i 0.00164 + 0.00496i 0.00205 + 0.0062i 0.00436 + 0.01322i 0.00341 + 0.01033i 10 0.00191 + 0.00579i 11 0.003 + 0.00909i 12 0.00245 + 0.00744i]; u = [0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 1; 1; 1; 1; 1; 1; 1; 1; 1; 1; 1; 1]; s = [0; -0.01 ; -0.006; -0.006; -0.002; -0.006; -0.008; -0.006; 0.0032; -0.01; -0.01; -0.016; 0; -0.0075; -0.0045; -0.0045; -0.0015; -0.0045; -0.006; -0.0045; -0.0024; -0.0075; -0.0075; -0.012]; attrib = [0; 2; 2; 2; 2; 2; 2; 2; 2; 2; 2; 2]; M = size(data); N = 0; for i = 1:M N = max(N,data(i,1)); end for i = 1:M N = max(N,data(i,2)); end ybus = zeros(N); for i = 1:M(1,1) if data(i,2) == ybus(data(i,1),data(i,1)) = ybus(data(i,1),data(i,1)) + 1/data(i,3); else ybus(data(i,1),data(i,1)) = ybus(data(i,1),data(i,1)) + 1/data(i,3); ybus(data(i,2),data(i,2)) = ybus(data(i,2),data(i,2)) + 1/data(i,3); ybus(data(i,1),data(i,2)) = -1/data(i,3); ybus(data(i,2),data(i,1)) = -1/data(i,3); end end N = length(ybus); fl = zeros(2*N,1); k = 5; for f =1:k power_temp = power_cal(u, ybus); delta_power = s - power_temp; ja = jacobi(u, ybus); for i=1:N if code(i,1) == fl(i,1) = 1; fl(i+N,1) = 1; elseif code(i,1) == fl(i+N,1) = 1; else continue; end end for i=1:2*N if fl(2*N+1-i,1) == ja(2*N+1-i,:)= []; ja(:,2*N+1-i)= []; delta_power(2*N+1-i,:)= []; end end change_temp = inv(ja)*delta_power; j = 1; for i=1:2*N if fl(i,1) == change(i,1) = 0; else change(i,1) = change_temp(j,1); j = j+1; end end u = u + change; end u m = length(u); for j= 1:m/2 a(j,1) = u((j+ m/2),1)*cos((u(j,1))); b(j,1) = u((j+ m/2),1)*sin((u(j,1)))*1i; U(j,1)= a(j,1)+ b(j,1); end for j=1:M(1,1) I(j,1)=(U(data(j,1),1)-U(data(j,2),1))/data(j,3); a(j,1)= real(I(j,1)); b(j,1)= imag(I(j,1))*1i; I_(j,1)= a(j,1) - b(j,1); S(j,1)=I_(j,1)*U(j+1,1); end S end function power = power_cal(u, ybus) N = length(ybus); for i=1:2*N power(i,1) = 0; end for i=1:N for n=1:N power(i,1) = power(i,1) + u(i+N,1)*u(n+N,1)*abs(ybus(i,n))*cos(angle(ybus(i,n))+ u(n,1)u(i,1)); power(i+N,1) = power(i+N,1) u(i+N,1)*u(n+N,1)*abs(ybus(i,n))*sin(angle(ybus(i,n))+u(n,1)-u(i,1)); end end end function Ja = jacobi(u, ybus) N = length(ybus); Ja = zeros(2*N); for i=1:N for n=1:N if n == i for j=1:N if j == i continue; else Ja(i,n) = Ja(i,n) + (u(i+N,1)*u(j+N,1)*abs(ybus(i,j))*sin(u(i,1)-u(j,1)angle(ybus(i,j)))); end end else Ja(i,n) = u(i+N,1)*u(n+N,1)*abs(ybus(i,n))*sin(u(i,1)-u(n,1)-angle(ybus(i,n))); end end end for i=1:N for n=(N+1):(2*N) if n == i+N for j=1:N if j == i Ja(i,n) = Ja(i,n) + 2*u(i+N,1)*abs(ybus(i,j))*cos(u(i,1)-u(j,1)-angle(ybus(i,j))); else Ja(i,n) = Ja(i,n) + u(j+N,1)*abs(ybus(i,j))*cos(u(i,1)-u(j,1)-angle(ybus(i,j))); end end else Ja(i,n) = u(i+N,1)*abs(ybus(i,n-N))*cos(u(i,1)u(n-N,1)-angle(ybus(i,n-N))); end end end for i=(N+1):(2*N) for n=1:N if n == i-N for j=1:N if j == i-N continue; else Ja(i,n) = Ja(i,n) + u(i,1)*u(j+N,1)*abs(ybus(i-N,j))*cos(u(i-N,1)-u(j,1)-angle(ybus(iN,j))); end end else Ja(i,n) = -u(i,1)*u(n+N,1)*abs(ybus(iN,n))*cos(u(i-N,1)-u(n,1)-angle(ybus(i-N,n))); end end end for i=(N+1):(2*N) for n=(N+1):(2*N) if n == i for j=1:N if j == i-N Ja(i,n) = Ja(i,n) + 2*u(i,1)*abs(ybus(iN,j))*sin(u(i-N,1)-u(j,1)-angle(ybus(i-N,j))); else Ja(i,n) = Ja(i,n) + u(j+N,1)*abs(ybus(iN,j))*sin(u(i-N,1)-u(j,1)-angle(ybus(i-N,j))); end end else Ja(i,n) = u(i,1)*abs(ybus(i-N,n-N))*sin(u(i-N,1)u(n-N,1)-angle(ybus(i-N,n-N))); end end end End PHỤ LỤC C MATLAB MÔ PHỎNG THUẬT TOÁN LFB ỨNG DỤNG VÀO MẠNG PHÂN PHỐI HÌNH TIA 22KV KHU CƠNG NGHIỆP VSIPI CĨ SỬ DỤNG THIẾT BỊ BÙ TCVR function [A,p,q,V] = LBF_TCVR(data,P,Q,R,X) clc; data =[ 0.00409 0.0124 0.00218 0.00661 0.00491 0.01488 0.00273 0.00826 0.00164 0.00496 0.00205 0.0062 0.00436 0.01322 0.00341 0.01033 10 0.00191 0.00579 11 0.003 0.00909 12 0.00245 0.00744]; P = [-0.01 ; -0.006; -0.006; -0.002; -0.006; -0.008; -0.006; 0.0032; -0.01; -0.01; -0.016]; Q = [-0.0075; -0.0045; -0.0045; -0.0015; -0.0045; -0.006; -0.0045; 0.0024; -0.0075; -0.0075; -0.012]; n = [2 4; 5; 6]; Vbu = [1; 1; 1]; K = 3; M = size(data); D = size(n); F = D(1,1); G = M(1,1); R = zeros(G); X = zeros(G); l = zeros(G,1); m = zeros(G,1); k = zeros(G,1); t = zeros(F,1); RX = zeros(G); pq = zeros(G); A_c = zeros(G,1); A_c(1,1)=1; V_pv = zeros(G,1); V_pv(1,1) = 1; N = 0; T = zeros(G); for i = 1:G T(i,i)=1; end for i = 1:F t(i,1)=1; end for i = 1:G N = max(N,data(i,2)); end for i = 1:N for j = 1:(N - 1) A(i,j) = 0; end end for i = 1:G A(data(i,1),data(i,2)-1) = 1; A(data(i,2),data(i,2)-1) = -1; R(i,i) = data(i,3); X(i,i) = data(i,4); end A (1,:)= []; A_1 = A > 0; A_2 = A < 0; A_2 = -1*A_2; B = ones (G,1); for i = 1: G RX(i,i) = R(i,i).^2 + X(i,i).^2; end p = A\(P - A_1*l); q = A\(Q - A_1*m); V_V = ((T*(A_1'))+A_2')\( k - T*A_c.*V_pv + 2*R*p + 2*X*q ); for i = 2:K for j = 1:F V_V(n(j,2)-1,1) = Vbu(j,1).^2; end for i = 1:G pq(i,i) = p(i,1).^2 + q(i,1).^2; C(i,1) = V_V(i,1)\B(i,1); end l = pq*R*(-A_2')*C; m = pq*X*(-A_2')*C; k = pq*RX*(-A_2')*C; for j = 1:F t(j,1) = ((V_V(n(j,2)-1,1)+2*(R(n(j,2)-1,n(j,2)-1)*p(n(j,2)1,1)+X(n(j,2)-1,n(j,2)-1)*q(n(j,2)-1,1))+k(n(j,2)1,1))\(V_V(n(j,1)-1,1))).^(1/2); end for j = 1:F T(n(j,2)-1,n(j,2)-1)= (t(j,1).^2)\1; end p = A\(P - A_1*l); q = A\(Q - A_1*m); V_V = ((T*(A_1'))+A_2')\( k - T*A_c.*V_pv + 2*R*p + 2*X*q ); for j = 1:F t(j,1) = ((V_V(n(j,2)-1,1)+2*(R(n(j,2)-1,n(j,2)-1)*p(n(j,2)1,1)+X(n(j,2)-1,n(j,2)-1)*q(n(j,2)-1,1))+k(n(j,2)1,1))\(V_V(n(j,1)-1,1))).^(1/2); end for j = 1:F T(n(j,2)-1,n(j,2)-1)= (t(j,1).^2)\1; end V = V_V.^(1/2); for i = 1:G pq(i,i) = p(i,1).^2 + q(i,1).^2; C(i,1) = V_V(i,1)\B(i,1); end end V t S = p+q*1i End PHỤ LỤC D MATLAB MÔ PHỎNG PHƯƠNG PHÁP NEWTON RAPHSON ỨNG DỤNG VÀO MẠNG PHÂN PHỐI HÌNH TIA 22KV KHU CƠNG NGHIỆP VSIPI CĨ SỬ DỤNG THIẾT BỊ BÙ TCVR basemva = 100; accuracy = busdata=[1 1 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 10 0.0 11 0.0 12 0.0 linedata=[1 0.00409 0.00218 0.00491 0.00273 0.00164 0.00205 0.00436 0.00341 10 0.00191 11 0.003 0.001; maxiter = 12; 0.0 0.0 0.0 0.0 1.0 0.75 0.0 0.0 0.6 0.45 0.0 0.0 0.6 0.45 0.0 0.0 0.2 0.15 0.0 0.0 0.6 0.45 0.0 0.0 0.8 0.60 0.0 0.0 0.6 0.45 0.0 0.0 0.32 0.24 0.0 0.0 1.0 0.75 0.0 0.0 1.0 0.75 0.0 0.0 1.6 1.20 0.0 0.0 0.0124 0.00661 0.01488 0.9984020001 0.00826 0.9983855792 0.00496 0.9984889369 0.0062 0.01322 0.01033 0.00579 0.00909 0 0 0 0 0 0 0 0 0 0 0 0]; 12 0.00245 0.00744 1]; j=sqrt(-1); i = sqrt(-1); nl = linedata(:,1); nr = linedata(:,2); R = linedata(:,3); X = linedata(:,4); Bc = j*linedata(:,5); a = linedata(:, 6); nbr=length(linedata(:,1)); nbus = max(max(nl), max(nr)); Z = R + j*X; y= ones(nbr,1)./Z; %branch admittance for n = 1:nbr if a(n)