Hong, “Reconfiguration of electric distri- bution networks for resistive line losses reduction,” IEEE Trans.. Lee, “Distribution feederreconfiguration for loss reduction,” IEEE Trans.. Pow
Trang 11
4
5
6
7
9
P 9
9
9
13
13
15
15
16
17
1.2.3.1 17
18
19
19
19
2.2 Cá 19
20
20
22
24
25
25
1708177943651a52c5a05-9b3b-4bf3-8979-dca038a0fead
Trang 2 27
27
27
27
29
-Raphson 31
31
31
32
34
4 39
39
39
39
40
4.1.3 Cá 41
44
45
50
54
54
54
58
60
67
72
76
78
79
Trang 5Hình 2.1: Gii thu t c a Merlin và Back 21
Hình 2.2: 23
Hình 3.1: Newton Raphson 32
Hình 3.2: 34
Hình 3.2: 38
Hình 4.1: 42
Hình 4.2: 2 .43
Hình 4.3: 44
Hình 4.4: 48
Hình 4.5: 49
Hình 4.6: 50
Hình 4.7: 51
Hình 4.8: 52
Hình 4.9: 52
Hình 5.1: 54
Hình 5.2: 55
Hình 5.3: 56
Hình 5.4: 472 E3.4, 475 E3.7 và 471 E3.14 69
Hình 5.5: 75
Hình 5.6: 77
Trang 6B ng 2.1: Phm vi ng d ng c a các bài toán tái c 20
B ng 2.2: 26
B ng 4.1 : 44
B ng 4.2 : 45
B ng 4.3 : 45
B ng 4.4 : 45
B ng 4.5 : 46
B ng 4.6 : 46
B ng 4.7 : 46
B ng 4.8 : 47
B ng 4.9 : 47
B ng 4.10 : 47
B ng 4.11 : 48
B ng 4.12: 48
B ng 4.13 : 49
B ng 4.14 : 49
B ng 4.15 Danh sách : 50
B ng 5.1 : 70
B ng 5.2 : 71
B ng 5.3 : 72
B ng 5.4 : 72
B ng 5.5 : 73
B ng 5.5 : 74
B ng 5.7 : 74
B ng 5.8 : 74
Trang 8Trong quá t
n thân, tôi thành
Trang 10
hàng
Trang 11
Trang 12 ,
c
Trang 13 NL/KHKT ngày 24 tháng 3
2
- : TBA g MBA (MBA): MBA T1
n áp 110/35/22 kV; MBA 110/35/22k
Trang 14MBA
- TBA có 2 MBA T1 25 MVA,
110/35/10 kV và MBA T2 25 MVA 110/35/22 kV TBA có
MBA
Trang 15- TBA có MBA 2 TBA có
Trang 17MBA MBA có công
, 400 kVA, 560kVA và 630kVA :
-
-E3.max dây 472- E3 max -E3.max là 8,118 MW
-
- max là 4,39
Trang 21* Theo D Shirmohammadi and H.W Hong, “Reconfiguration of electric distri-
bution networks for resistive line losses reduction,” IEEE Trans Power
Del., vol 4, no 2, pp 1492 1498, Apr 1989–
Trang 22* Theo S Civanlar, J J Grainger, H Yin, and S S H Lee, “Distribution feeder
reconfiguration for loss reduction,” IEEE Trans Power Del., vol 3, no.
3, pp 1217 1223, Jul 1988.–
Trang 24* Theo T E McDermott, I Drezga, and R P Broadwater, “A heuristic non-
linear constructive method for distribution system reconfiguration,”
IEEE Trans Power Syst., vol 1 no 2, pp 478 483, May 1999 4, –
Trang 26 tính toán hình
Trang 28 ,
- = ,
(2.6) (2.7)
Trang 32x2 = - + x1 (2.23)
k+1 k: Xk+1 - = + xk (2.24)
kk - f(x= k) (2.25) Tk = xk+1 - xk
kk theo xk+1 k
Hình 3
3.2.1.2
1 x2n]T
Trang 33F(x) = (2.26)
i hóa), ta có:
Trang 34= Ui = Ui(cos + jsinCông sut ph c liên h p t i nút i:
Pi jQi= = (2.30) Khai tring ph n th c và o ta có:
Pi = (2.31)
Qi - = (2.32)
u Gán i = 0
Trang 35 Pi,d và Qi,d là công sut tác d ng và ph nút i
Pi,t vàQi,t là công su t tác d ng và ph c t i nút i theo công thc (10) và (11)
Trong n+1
0 và n nút: nc nút PQ và ng nút PV
có bi n tr
i ta có:
i = (i=1,nc) (2.39)
Trang 36(2.42)
1
Ui{2Ui.Gii + (2.43)Các ph n t ng chéo c a ma tr n J 2c xnh b ng:
Trang 37nh theo các công th c sau:
m l n là t h i t nhanh n u sai s u nh ( chênh l ch gi a nghi m
u và nghi m th c) Trái l i n u sai s u l n, quá trình l p có th phân kì
Nc x p x u th a mãn, ta có th chn giá tr n áp các nút PQ bnh mn áp các nút PV theo giá tr a các nút ly b ng 0 (b ng nút cân b ng)
Tuy th i gian tính m c l p c -Raphson l n,
c l p yêu cng nh , h thu c vào s nút c a
d i cho tính toán h
thn có s nút b t k
Trang 38
S
Raphson trong tính
Trang 42g Sai
Trang 44Nút 4: nút PQ, cho S4 -400 50i (kVA) =
Nút 5: nút PQ, cho S5 -600 100i (kVA) =
Trang 48Hình 4.4:
các khóasách các khóa Khóa
1
Trang 491
Trang 50các khóa => danh sách các khóa
Danh sách các khóa tia, các khóa
1
Trang 511
Trang 521
Trang 545 TÍNH TOÁN
Trang 55
Nút Bus data
Hình 5.2:
Nút Line data
Trang 57% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
busdatas = xlsread('busdatas.xlsx');
% hObject handle to Nhap_nhanh (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
linedatas = xlsread('linedatas.xlsx');
% - Executes on button press in Chay
% hObject handle to Chay (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
Trang 58% handles structure with handles and user data (see GUIDATA)
xlswrite('khoa_mo.xlsx', khoamo2);
set(handles.delta_p,'String',num2str(deltaP)); winopen('khoa_mo.xlsx');
Trang 59else if nhanh(m,3) == kiemtra(n)
Trang 61if isempty(s) && khoamo1(l,t) ~= 0
else if nhanh(m,3) == kiemtra(n)
Trang 62 Newton-
% Program for Newton-Raphson Load Flow Analysis
% Calling ybusppg.m to get
Y = ybusppg(lined);
Y-Bus Matrix
% Base kVA
BMva = 100000;
bus = busd(:,1); % Bus Number
type = busd(:,2); % Type of Bus 1-Slack,
2-PV, 3-PQ
V = busd(:,3); % Specified Voltage
del = busd(:,4); % Voltage Angle
G = real(Y); % Conductance matrix
B = imag(Y); % Susceptance matrix
pv = find(type == 2 | type == 1); % PV Buses
Q(i) = Q(i) + V(i)*
V(k)*(G(i,k)*sin(del(i)-del(k)) - del(k)));
end
Trang 63end
if Iter <= 7 && Iter > 2 % Only checked up to 7th iterations
V(n)*(-del(n)));
end
J1(i,k) = J1(i,k) - V(m)^2*B(m,m); else
J1(i,k) = V(m)*
V(n)*(G(m,n)*sin(del(m)-del(n)) - del(n)));
Trang 64end
J2(i,k) = J2(i,k) + V(m)*G(m,m); else
J2(i,k) = del(n)) + B(m,n)*sin(del(m)-del(n)));
end
J3(i,k) = J3(i,k) - V(m)^2*G(m,m); else
end
Trang 65J4(i,k) = J4(i,k) - V(m)*B(m,m); else
J4(i,k) = del(n)) - B(m,n)*cos(del(m)-del(n)));
end
end
end
J = [J1 J2; J3 J4]; % Jacobian Matrix
X = inv(J)*M; % Correction Vector
dTh = X(1:nbus-1); % Change in Voltage
Angle
dV = X(nbus:end); % Change in Voltage
Magnitude
del(2:nbus) = dTh + del(2:nbus); % Voltage Angle
fb = lined(:,2); % From bus number
tb = lined(:,3); % To bus number
nl = length(fb); % No of Branches
Trang 66fb = lined(:,2); % From bus number
tb = lined(:,3); % To bus number
b = 1i*b; % Make B imaginary
Trang 67nb = max(max(fb),max(tb)); % No of buses
nl = length(fb); % No of branches
Y = zeros(nb,nb); % Initialise YBus
% Formation of the Off Diagonal Elements
%Y; % Bus Admittance Matrix
%Z = inv(Y); % Bus Impedance Matrix
-g
-% Polar to Rectangular Conversion
% [RECT] = RECT2POL(RHO, THETA)
% RECT - Complex matrix or number, RECT = A + jB, A = Real, B = Imaginary
% RHO - Magnitude
% THETA - Angle in radians
rect = rho.*cos(theta) + j*rho.*sin(theta);