Lựa chọn thông số và vị trí đặt các bộ lọc sóng hài thụ động

MỤC LỤC LỜI CAM ĐOAN LỜI CẢM ƠN TÓM TẮT NỘI DUNG DANH MỤC BẢNG BIỂU DANH MỤC HÌNH VẼ CHƢƠNG GIỚI THIỆU VỀ SÓNG HÀI TRONG HỆ THỐNG ĐIỆN 1.1 Khái niệm sóng hài 1.2 Các số đánh giá sóng hài hệ thống điện 10 1.3 Nguồn phát sinh sóng hài hệ thống điện 13 1.4 Ảnh hƣởng sóng hài tới thiết bị hệ thống điện 18 CHƢƠNG GIỚI THIỆU VỀ BỘ LỌC SÓNG HÀI THỤ ĐỘNG 22 2.1 Tổng quan lọc sóng hài thụ động 22 2.2 Phân tích lọc cộng hƣởng đơn 26 2.2.1 Đặc tính tổng trở lọc cộng hƣởng đơn 26 2.2.2 Phân tích hệ số chất lƣợng Q cho lọc cộng hƣởng đơn 2827 2.2.3 Phân tích tần số cộng hƣởng cho lọc cộng hƣởng đơn 29 a Lý không nên chọn tần số cộng hƣởng tần số sóng hài 29 b Các yếu tố gây tƣợng dịch chuyển tần số cộng hƣởng 31 2.3 Hƣớng nghiên cứu luận văn 32 CHƢƠNG NGHIÊN CỨU TÍNH TOÁN LỰA CHỌN THÔNG SỐ CHO CÁC BỘ LỌC CỘNG HƢỞNG ĐƠN VÀ VỊ TRÍ LẮP ĐẶT 3433 3.1 Xây dựng vấn đề nghiên cứu 3433 3.2 Qui trình thiết kế nhánh lọc theo ràng buộc bù công suất phản kháng4038 3.3 Nghiên cứu bBài toán phân bổ công suất phản kháng tối ƣu cho nhánh lọc nhóm lọc 4542 3.3.1 Đặt vấn đề 4642 3.3.2 Xác định chi phí cho tụ kháng nhánh lọc 4743 3.3.3 Xây dựng hàm mục tiêu ràng buộc 5147 3.3.4 3.4 Phƣơng pháp giải toán tối ƣu 5349 Nghiên cứu ràng buộc lựa chọn vị trí đặt nhóm lọc 5450 CHƢƠNG TÍNH TOÁN KIỂM NGHIỆM VÀ KẾT LUẬN 5852 4.1 Giới thiệu phần mềm PSCAD mô hình nghiên cứu 5852 4.1.1 Tính chung PSCAD 5852 4.1.2 Mô hình, thông số đầu vào phục vụ tính toán 6054 4.2 Giới thiệu phần mềm Matlab công cụ tính toán tối ƣu sẵn có 70595960 4.2.1 Giới thiệu phần mềm Matlab 70595960 4.2.2 Công cụ tính toán tối ƣu sẵn có Matlab 73626263 4.2.3 Kết tính toán dựa Matlab nhận xét 7363 KẾT LUẬN: 67676768 CHƢƠNG KẾT LUẬN VÀ HƢỚNG NGHIÊN CỨU TRONG TƢƠNG LAI 70696969 5.1 Kết luận 70696969 5.2 Hƣớng nghiên cứu tƣơng lai 71707070 TÀI LIỆU THAM KHẢO 73717171 PHỤ LỤC 75727272 LỜI CAM ĐOAN Tôi xin cam đoan luận văn kết nghiên cứu riêng tôi, không chép Các số liệu, kết nêu luận văn trung thực chƣa đƣợc công bố công trình khác Nội dung luận văn có tham khảo sử dụng tài liệu, thông tin đƣợc đăng tải tác phẩm, tạp chí, báo trang web theo danh mục tài liệu tham khảo luận văn Tác giả Hà Đức Nhẹn Formatted: Indent: First line: 0.5" LỜI CẢM ƠN Tôi xin chân thành cảm ơn hƣớng dẫn nhiệt tình TS Nguyễn Xuân Tùng, giảng viên Bộ môn Hệ thống điện, Viện Điện, Trƣờng đại học Bách khoa Hà Nội hƣớng dẫn hoàn thành luận văn Tôi xin gửi lời cảm ơn sâu sắc tới thầy cô tham gia giảng dạy khóa học, thầy cô Viện Điện, tạo điều kiện thuận lợi giúp hoàn thành khóa học Tôi xin gửi lời cảm ơn tới cán hành Viện Điện Viện Đào tạo Sau đại học giúp đỡ trình học tập trƣờng Lời cuối cùng, chân thành cảm ơn động viên gia đình, bạn bè, ngƣời tạo điều kiện nhiều cho suốt chặng đƣờng học tập qua Hà Đức Nhẹn TÓM TẮT NỘI DUNG Một vấn đề chất lƣợng điện phổ biến hệ thống điện sóng hài Sóng hài đƣợc sinh lan truyền hệ thống việc sử dụng tải phi tuyến, thiết bị điện tử công suất, lò hồ quang lò cảm ứng… Sóng hài hệ thống điện gây nhiều vấn đề thiết bị có lõi từ nhƣ gây phát nóng mức, gây rung động thiết bị quay, làm tải dây trung tính nhiều vấn đề khác cho thiết bị điều khiển Giải pháp loại trừ sóng hài hệ thống điện chia ba nhánh chính: sử dụng lọc thụ động; sử dung lọc tích cực lọc lai ghép hai dạng Bộ lọc sóng hài thụ động đƣợc lựa chọn đối tƣợng tính toán nghiên cứu luận văn có chi phí đầu tƣ bảo dƣỡng thấp, phù hợp với công nghiệp Trong luận văn tìm hiểu qui trình tính toán thiết kế lọc thụ động áp dụng nghiên cứu lựa chọn vị trí đặt tối ƣu cho lọc thụ động lƣới điện…Mô hình áp dụng luận văn lƣới điện phân phối Công ty TNHH Điện Stanley Việt Nam (lƣới điện cấp cho nhà máy II) Do trình toán toán phức tạp với nhiều số liệu nên luận văn sử dụng phần mềm MATLAB để hỗ trợ tính toán Về mặt cấu trúc luận văn đƣợc chia thành 54 chƣơng  Chương 1: Giới thiệu chung sóng hài hệ thống điện Các nguyên nhân gây sóng hài ảnh hƣởng sóng hài tới hệ thống thiết bị Các tiêu chuẩn qui định mức độ méo sóng hài cho phép  Chương 2: Giới thiệu lọc sóng hài thụ động Phân tích cấu trúc, ƣu nhƣợc điểm giải pháp lọc sóng hài sử dụng thiết bị lọc thụ động Đồng thời đặt vấn đề nghiên cứu chuyên sâu luận văn  Chương 3: Xây dựng qui trình tính toán thiết kế, kiểm tra thông số lọc thụ động Xây dựng ràng buộc sử dụng nhiều lọc thụ động có nhiều vị trí lựa chọn để đặt lọc Formatted: Font: Not Italic  Chương 4: Áp dụng tính toán với mô hình lƣới điện phân phối nhà máy II Công ty TNHH Stanley Việt Nam đánh giá kết đạt đƣợc  Chương 5: Kết luận đề xuất hƣớng nghiên cứu tƣơng lai DANH MỤC BẢNG BIỂU Bảng 1.2.1 Tiêu chuẩn điện áp theo Thông tƣ 12 32 11 Bảng 1.2.2 Giới hạn tổng độ méo nhu cầu dòng điện hệ thống điện phân phối có điện áp từ 120V đến 69kV 12 Bảng 1.2.3 Tổng méo dạng nhu cầu dòng điện hệ thống điện truyền tải 69kV đến 161kV 12 Bảng 1.2.4 Tổng méo dạng nhu cầu dòng điện hệ thống điện truyền tải nguồn phân tán 13 Bảng 1.3.1 Dạng sóng dòng điện, phổ tần tổng độ méo sóng hài số tải phi tuyến khác 18 Bảng 4.1.1 Dòng điện dòng điện hài vị trí đặt lọc (tính theo Ampe)………………………………………………………………………………57 Bảng 4.1.2 CSPK yêu cầu vị trí lắp đặt lọc………………………… 578 DANH MỤC HÌNH VẼ Hình 1.1.1 Phân tích Fourer sóng bị méo dạng .9 Hình 1.1.2 Sóng méo dạng phân tích Furier tƣơng ứng .10 Hình 1.3.1 Sóng hài sinh tƣợng bão hòa mạch từ máy biến áp .14 Hình 1.3.2 Dạng sóng phổ dòng pha A máy biến áp hoạt động điều kiện áp 10% điện áp định mức 14 Hình 1.3.3 Dòng điện máy lạnh với THD=6,3% .14 Hình 1.3.4 Dòng điện điều hòa với THD=10,5% .15 Hình 1.3.5 Dạng sóng điện áp & dòng pha A chỉnh lƣu .15 Hình 1.3.6 Nguồn chuyển mạch pha dạng sóng dòng điện 16 Hình 1.3.7 Sóng dòng điện phổ tần đèn huỳnh quang 17 Hình 2.1.1 Cấu hình loại lọc thụ động phổ biến .22 Hình 2.1.2 Mô hình hệ thống điện có tải phi tuyến 23 Hình 2.1.3 Sơ đồ tƣơng đƣơng mô hình hệ thống có tải phi tuyến 24 Hình 2.1.4 Đặc tính tổng trở độ rộng băng thông lọc cộng hƣởng đơn 25 Hình 2.2.1 Đặc tính tổng trở theo tần số lọc cộng hƣởng đơn 27 Hình 2.2.2 Bộ lọc cộng hƣởng đơn đặc tính tổng trở theo tần số 28 Hình 2.2.3 Đặc tính tổng trở lọc với hệ số chất lƣợng khác 29 Hình 2.2.4 Sơ đồ nối lọc hệ thống điện có tải phi tuyến 30 Hình 2.2.5 Các tần số cộng hƣởng xuất có lọc hệ thống 31 Hình 3.1.1 Mô tả nhánh & nhóm lọc .3534 Hình 3.1.2 Minh họa phân bổ công suất phản kháng cho nhánh lọc 3635 Hình 3.1.3 Vị trí đặt lọc 3836 Hình 3.1.4 Vấn đề nghiên cứu luận văn 4037 Hình 3.2.1 Sơ đồ khối trình tính toán lọc cộng hƣởng đơn .4541 Hình 3.4.1 Các vị trí lắp đặt lọc sóng hài 5650 Hình 4.1.1 Giao diện làm việc PSCAD……………………………………… 52 Hình 4.1.2 Thƣ viện linh kiện PSCAD……………………………………54 Hình 4.1.3 Sơ đồ sợi lƣới điện nhà máy II công ty TNHH Điện Stanley Việt Nam…………………………………………………………… ……………55 Hình 4.1.4 Sơ đồ lƣới điện nhà máy II công ty TNHH Điện Stanley Việt Nam mô phần mềm PSCAD………………………………………………….….56 Hình 4.2.1 Giao diện MATLAB………………………………….…… 612 Hình 4.2.2 Giao diện cửa sổ soạn thảo lệnh……………………… 62 Hình 4.2.3 Sơ đồ sợi lƣới điện sau lắp đặt lọc vị trí thứ cấp máy biến Formatted: Font: Not Bold áp………………………………………………………………… ………………68 Formatted: Left Formatted: Font: Not Italic Formatted: Font: Not Italic CHƢƠNG GIỚI THIỆU VỀ SÓNG HÀI TRONG HỆ THỐNG ĐIỆN 1.1 Khái niệm sóng hài Sóng hài dạng nhiễu không mong muốn xuất dƣới dạng dòng điện hay điện áp có tần số số nguyên lần tần số nguồn cung cấp (thƣờng đƣợc gọi tần số sóng bản) Các dòng điện, điện áp bị méo đƣợc phân tích thành tổng sóng có tần số thành phần sóng hài Các thành phần sóng hài tải phi tuyến sinh Công cụ toán học để phân tích mức độ méo dạng sóng dòng điện có chu kỳ phân tích Fourier Phƣơng pháp dựa nguyên lý dạng sóng không sin, có chu kỳ phân tích đƣợc thành tổng dạng sóng điều hòa hình sin, bao gồm:  Sóng hình sin với tần số  Các sóng hình sin khác với tần số hài cao hơn, bội tần số Hình 1.1.1 Phân tích Fourer sóng bị méo dạng Trong trƣờng hợp lý tƣởng, tất sóng điện áp dòng điện hệ thống điện có dạng hình sin với tần số tần số bản.Tuy nhiên, điện áp dòng điện thực tế hệ thống điện không túy hình sin Khi đó, sóng điện áp dòng điện tổng sóng điều hòa sóng điều hòa có bậc bội số sóng K7TU=(MC+ML)+UC*3*Anpha7*V1*sqrt(I7^2*(1+1/h7^2+((3*I7*V1)/(h7*x(2)))^2)+ (x(2)/(3*V1))^2)+ 3*UL*Anpha7*V1/h7^2*sqrt(I7^2*(1+h7^2+((3*h7*I7*V1)/(x(2)))^2)+(x(2)/(3*V 1))^2)+ (Pv*8760*kc*x(2)*c)+(Pv*8760*c*((x(2)^2/(3*V1^2))+I7^2)*h7*V1^2/((h7^21)*Q*x(2))) K11TU= (MC+ML)+UC*3*Anpha11*V1*sqrt(I11^2*(1+1/h11^2+((3*I11*V1)/(h11*x(3)))^2)+ (x(3)/(3*V1))^2)+ 3*UL*Anpha11*V1/h11^2*sqrt(I11^2*(1+h11^2+((3*h11*I11*V1)/(x(3)))^2)+(x(3 )/(3*V1))^2)+ (Pv*8760*kc*x(3)*c)+(Pv*8760*c*((x(3)^2/(3*V1^2))+I11^2)*h11*V1^2/((h11^2 -1)*Q*x(3))) K13TU=(MC+ML)+UC*3*Anpha13*V1*sqrt(I13^2*(1+1/h13^2+((3*I13*V1)/(h13*x(4) ))^2)+(x(4)/(3*V1))^2)+ 3*UL*Anpha13*V1/h13^2*sqrt(I13^2*(1+h13^2+((3*h13*I13*V1)/(x(4)))^2)+(x(4 )/(3*V1))^2)+ (Pv*8760*kc*x(4)*c)+(Pv*8760*c*((x(4)^2/(3*V1^2))+I13^2)*h13*V1^2/((h13^2 -1)*Q*x(4))) Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt % Tinh toan voi truong hop luong CSPK duoc chia deu cho cac bo loc Formatted: Font: 10 pt x(1)=Qyeucau/4; x(2)=x(1); x(3)=x(1); x(4)=x(1); Formatted: Font: 10 pt Formatted: Font: 10 pt disp('Chi phi tinh toan cho tung bo loc theo phan bo deu') Formatted: Font: 10 pt K5PBD=(MC+ML)+UC*3*Anpha5*V1*sqrt(I5^2*(1+1/h5^2+((3*I5*V1)/(h5*x(1)))^2) +(x(1)/(3*V1))^2)+ 3*UL*Anpha5*V1/h5^2*sqrt(I5^2*(1+h5^2+((3*h5*I5*V1)/(x(1)))^2)+(x(1)/(3*V 1))^2)+ (Pv*8760*kc*x(1)*c)+(Pv*8760*c*((x(1)^2/(3*V1^2))+I5^2)*h5*V1^2/((h5^21)*Q*x(1))) Formatted: Font: 10 pt K7PBD=(MC+ML)+UC*3*Anpha7*V1*sqrt(I7^2*(1+1/h7^2+((3*I7*V1)/(h7*x(2)))^2) +(x(2)/(3*V1))^2)+ 3*UL*Anpha7*V1/h7^2*sqrt(I7^2*(1+h7^2+((3*h7*I7*V1)/(x(2)))^2)+(x(2)/(3*V 1))^2)+ (Pv*8760*kc*x(2)*c)+(Pv*8760*c*((x(2)^2/(3*V1^2))+I7^2)*h7*V1^2/((h7^21)*Q*x(2))) Formatted: Font: 10 pt K11PBD= (MC+ML)+UC*3*Anpha11*V1*sqrt(I11^2*(1+1/h11^2+((3*I11*V1)/(h11*x(3)))^2)+ (x(3)/(3*V1))^2)+ 3*UL*Anpha11*V1/h11^2*sqrt(I11^2*(1+h11^2+((3*h11*I11*V1)/(x(3)))^2)+(x(3 )/(3*V1))^2)+ (Pv*8760*kc*x(3)*c)+(Pv*8760*c*((x(3)^2/(3*V1^2))+I11^2)*h11*V1^2/((h11^2 -1)*Q*x(3))) K13PBD=(MC+ML)+UC*3*Anpha13*V1*sqrt(I13^2*(1+1/h13^2+((3*I13*V1)/(h13*x(4 )))^2)+(x(4)/(3*V1))^2)+ 3*UL*Anpha13*V1/h13^2*sqrt(I13^2*(1+h13^2+((3*h13*I13*V1)/(x(4)))^2)+(x(4 )/(3*V1))^2)+ (Pv*8760*kc*x(4)*c)+(Pv*8760*c*((x(4)^2/(3*V1^2))+I13^2)*h13*V1^2/((h13^2 -1)*Q*x(4))) KQ_phan_bo_deu =(K5PBD+K7PBD+K11PBD+K13PBD); Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt 85 disp('Ham chi phi tinh toan theo phan bo deu'); disp(KQ_phan_bo_deu) Hieu_qua=100*(KQ_phan_bo_deu-fval)/KQ_phan_bo_deu Formatted: Font: 10 pt % Tham so cua cac phan tu cua cac bo loc (C: F; L; H) Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt disp('Tham so cua cac phan tu voi phan bo toi uu') Xc5=(h5^2*V1^2/(h5^2-1))/Q5; C5=1/(2*pi*50*Xc5*1000) XL5=Xc5/h5^2; L5=XL5*1000/(2*pi*50) Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Xc7=(h7^2*V1^2/(h7^2-1))/Q7; C7=1/(2*pi*50*Xc7*1000) XL7=Xc7/h7^2; L7=XL7*1000/(2*pi*50) Formatted: Font: 10 pt Formatted: Font: 10 pt Xc11=(h11^2*V1^2/(h11^2-1))/Q11; C11=1/(2*pi*50*Xc11*1000) XL11=Xc11/h11^2; L11=XL11*1000/(2*pi*50) Formatted: Font: 10 pt Xc13=(h13^2*V1^2/(h13^2-1))/Q13; C13=1/(2*pi*50*Xc13*1000) XL13=Xc13/h13^2; L13=XL13*1000/(2*pi*50) Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt disp('Tham so cua cac phan tu voi phan bo deu') Xc5=(h5^2*V1^2/(h5^2-1))/x(1); C5=1/(2*pi*50*Xc5*1000) XL5=Xc5/h5^2; L5=XL5*1000/(2*pi*50) Formatted: Font: 10 pt Xc7=(h7^2*V1^2/(h7^2-1))/x(2); C7=1/(2*pi*50*Xc7*1000) XL7=Xc7/h7^2; L7=XL7*1000/(2*pi*50) Formatted: Font: 10 pt Xc11=(h11^2*V1^2/(h11^2-1))/x(3); C11=1/(2*pi*50*Xc11*1000) XL11=Xc11/h11^2; L11=XL11*1000/(2*pi*50) Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Xc13=(h13^2*V1^2/(h13^2-1))/x(4); C13=1/(2*pi*50*Xc13*1000) XL13=Xc13/h13^2; L13=XL13*1000/(2*pi*50) clc clear all close all format short Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt % Truong hop bo loc dat tai vi tri thu cap may bien ap phan phoi Formatted: Font: 10 pt % Xet tung nhanh thu cap may bien ap - Nhap gia tri tuong ung tung nhanh; Formatted: Font: 10 pt Formatted: Font: 10 pt % Dien ap dinh muc thu cap may bien ap; Formatted: Font: 10 pt V1 = input('Nhap vao dien ap thu cap may bien ap phan phoi V1(kV): '); % Bien cac phan song hai cua nhanh lap dat bo loc I5 = input('Nhap vao dong dien hai bac cua nhanh thu cap I5(A): '); I7 = input('Nhap vao dong dien hai bac cua nhanh thu cap I7(A): '); I11 = input('Nhap vao dong dien hai bac 11 cua nhanh thu cap I11(A): '); I13 = input('Nhap vao dong dien hai bac 13 cua nhanh thu cap I13(A): '); % Tong luong CSPK yeu cau cua nhanh thu cap lap dat bo loc Qyeucau = input('Nhap vao CSPK yeu cau cua nhanh thu cap Qyeucau(VAr): '); 86 % -Bac song hai can loai tru h5=4.75; h7=6.65; h11=10.45; h13=12.35; Formatted: Font: 10 pt, English (U.S.) % -Chi phi cho tu va khang MC=1; % Chi phi co dinh cua tu dien UC=1; % He so chi phi phu thuoc ML=17.5;% Chi phi co dinh cua cuon khang UL=17.5; % He so chi phi phu thuoc cua cuon khang Formatted: Font: 10 pt % -He so anpha cua tung bac song hai Anpha5 = 1.046; Anpha7 = 1.023; Anpha11 = 1.009; Anpha13 = 1.007; % -He so chat luong cua bo loc Q=40; % -He so ton hao cong suat bo tu kc=0.0004; % He so quy doi ve thoi diem hien tai Pv=12.46 % Tinh vong doi bo loc la 20 nam voi lai suat la 0.05 Formatted: Font: 10 pt, English (U.S.) % -Don gia dien nang chi phi ton hao tu C va dien tro R bo loc c=0.1 Formatted: Font: 10 pt %===========================================================% Thiet lap ham fun=@(x)(MC+ML)+UC*3*Anpha5*V1*sqrt(I5^2*(1+1/h5^2+((3*I5*V1)/(h5*x(1)))^ 2)+(x(1)/(3*V1))^2)+ 3*UL*Anpha5*V1/h5^2*sqrt(I5^2*(1+h5^2+((3*h5*I5*V1)/(x(1)))^2)+(x(1)/(3*V 1))^2)+ (Pv*8760*kc*x(1)*c)+(Pv*8760*c*((x(1)^2/(3*V1^2))+I5^2)*h5*V1^2/((h5^21)*Q*x(1)))+ (MC+ML)+UC*3*Anpha7*V1*sqrt(I7^2*(1+1/h7^2+((3*I7*V1)/(h7*x(2)))^2)+(x(2) /(3*V1))^2)+ 3*UL*Anpha7*V1/h7^2*sqrt(I7^2*(1+h7^2+((3*h7*I7*V1)/(x(2)))^2)+(x(2)/(3*V 1))^2)+ (Pv*8760*kc*x(2)*c)+(Pv*8760*c*((x(2)^2/(3*V1^2))+I7^2)*h7*V1^2/((h7^21)*Q*x(2)))+ (MC+ML)+UC*3*Anpha11*V1*sqrt(I11^2*(1+1/h11^2+((3*I11*V1)/(h11*x(3)))^2)+ (x(3)/(3*V1))^2)+ 3*UL*Anpha11*V1/h11^2*sqrt(I11^2*(1+h11^2+((3*h11*I11*V1)/(x(3)))^2)+(x(3 )/(3*V1))^2)+ (Pv*8760*kc*x(3)*c)+(Pv*8760*c*((x(3)^2/(3*V1^2))+I11^2)*h11*V1^2/((h11^2 -1)*Q*x(3)))+ (MC+ML)+UC*3*Anpha13*V1*sqrt(I13^2*(1+1/h13^2+((3*I13*V1)/(h13*x(4)))^2)+ (x(4)/(3*V1))^2)+ 3*UL*Anpha13*V1/h13^2*sqrt(I13^2*(1+h13^2+((3*h13*I13*V1)/(x(4)))^2)+(x(4 )/(3*V1))^2)+ (Pv*8760*kc*x(4)*c)+(Pv*8760*c*((x(4)^2/(3*V1^2))+I13^2)*h13*V1^2/((h13^2 -1)*Q*x(4))) 87 % Cac rang buoc constr=@(x)sum(x); Aeq = [1 1 1]; beq = Qyeucau; lb = zeros(1,4); Formatted: Font: 10 pt, English (U.S.) % Nghiem khoi dau va cac lua chon tinh x0 = [10;10;10;10]; options = optimset('Algorithm','interior-point'); Formatted: Font: 10 pt % Lenh tinh toan tim cuc tieu [x, fval] = fmincon(fun,x0,[],[],Aeq,beq,lb,[],[],options); Formatted: Font: 10 pt, English (U.S.) % Hien thi ket qua disp('Cong suat cua moi bo loc theo tinh toan phan bo toi uu (kVar)') disp(['x = ']);disp(x) % Don vi cua cong suat la kVAr disp('Ham chi phi tinh toan theo phan bo toi uu'); disp(fval) Formatted: Font: 10 pt Q5=x(1); Q7=x(2); Q11=x(3); Q13=x(4); disp('Chi phi tinh toan cho tung bo loc theo phan bo toi uu') K5TU=(MC+ML)+UC*3*Anpha5*V1*sqrt(I5^2*(1+1/h5^2+((3*I5*V1)/(h5*x(1)))^2)+ (x(1)/(3*V1))^2)+ 3*UL*Anpha5*V1/h5^2*sqrt(I5^2*(1+h5^2+((3*h5*I5*V1)/(x(1)))^2)+(x(1)/(3*V 1))^2)+ (Pv*8760*kc*x(1)*c)+(Pv*8760*c*((x(1)^2/(3*V1^2))+I5^2)*h5*V1^2/((h5^21)*Q*x(1))) K7TU=(MC+ML)+UC*3*Anpha7*V1*sqrt(I7^2*(1+1/h7^2+((3*I7*V1)/(h7*x(2)))^2)+ (x(2)/(3*V1))^2)+ 3*UL*Anpha7*V1/h7^2*sqrt(I7^2*(1+h7^2+((3*h7*I7*V1)/(x(2)))^2)+(x(2)/(3*V 1))^2)+ (Pv*8760*kc*x(2)*c)+(Pv*8760*c*((x(2)^2/(3*V1^2))+I7^2)*h7*V1^2/((h7^21)*Q*x(2))) K11TU= (MC+ML)+UC*3*Anpha11*V1*sqrt(I11^2*(1+1/h11^2+((3*I11*V1)/(h11*x(3)))^2)+ (x(3)/(3*V1))^2)+ 3*UL*Anpha11*V1/h11^2*sqrt(I11^2*(1+h11^2+((3*h11*I11*V1)/(x(3)))^2)+(x(3 )/(3*V1))^2)+ (Pv*8760*kc*x(3)*c)+(Pv*8760*c*((x(3)^2/(3*V1^2))+I11^2)*h11*V1^2/((h11^2 -1)*Q*x(3))) K13TU=(MC+ML)+UC*3*Anpha13*V1*sqrt(I13^2*(1+1/h13^2+((3*I13*V1)/(h13*x(4) ))^2)+(x(4)/(3*V1))^2)+ 3*UL*Anpha13*V1/h13^2*sqrt(I13^2*(1+h13^2+((3*h13*I13*V1)/(x(4)))^2)+(x(4 )/(3*V1))^2)+ (Pv*8760*kc*x(4)*c)+(Pv*8760*c*((x(4)^2/(3*V1^2))+I13^2)*h13*V1^2/((h13^2 -1)*Q*x(4))) % Tinh toan voi truong hop luong CSPK duoc chia deu cho cac bo loc 88 x(1)=Qyeucau/4; x(2)=x(1); x(3)=x(1); x(4)=x(1); Formatted: Font: 10 pt, English (U.S.) disp('Chi phi tinh toan cho tung bo loc theo phan bo deu') Formatted: Font: 10 pt K5PBD=(MC+ML)+UC*3*Anpha5*V1*sqrt(I5^2*(1+1/h5^2+((3*I5*V1)/(h5*x(1)))^2) +(x(1)/(3*V1))^2)+ 3*UL*Anpha5*V1/h5^2*sqrt(I5^2*(1+h5^2+((3*h5*I5*V1)/(x(1)))^2)+(x(1)/(3*V 1))^2)+ (Pv*8760*kc*x(1)*c)+(Pv*8760*c*((x(1)^2/(3*V1^2))+I5^2)*h5*V1^2/((h5^21)*Q*x(1))) K7PBD=(MC+ML)+UC*3*Anpha7*V1*sqrt(I7^2*(1+1/h7^2+((3*I7*V1)/(h7*x(2)))^2) +(x(2)/(3*V1))^2)+ 3*UL*Anpha7*V1/h7^2*sqrt(I7^2*(1+h7^2+((3*h7*I7*V1)/(x(2)))^2)+(x(2)/(3*V 1))^2)+ (Pv*8760*kc*x(2)*c)+(Pv*8760*c*((x(2)^2/(3*V1^2))+I7^2)*h7*V1^2/((h7^21)*Q*x(2))) K11PBD= (MC+ML)+UC*3*Anpha11*V1*sqrt(I11^2*(1+1/h11^2+((3*I11*V1)/(h11*x(3)))^2)+ (x(3)/(3*V1))^2)+ 3*UL*Anpha11*V1/h11^2*sqrt(I11^2*(1+h11^2+((3*h11*I11*V1)/(x(3)))^2)+(x(3 )/(3*V1))^2)+ (Pv*8760*kc*x(3)*c)+(Pv*8760*c*((x(3)^2/(3*V1^2))+I11^2)*h11*V1^2/((h11^2 -1)*Q*x(3))) K13PBD=(MC+ML)+UC*3*Anpha13*V1*sqrt(I13^2*(1+1/h13^2+((3*I13*V1)/(h13*x(4 )))^2)+(x(4)/(3*V1))^2)+ 3*UL*Anpha13*V1/h13^2*sqrt(I13^2*(1+h13^2+((3*h13*I13*V1)/(x(4)))^2)+(x(4 )/(3*V1))^2)+ (Pv*8760*kc*x(4)*c)+(Pv*8760*c*((x(4)^2/(3*V1^2))+I13^2)*h13*V1^2/((h13^2 -1)*Q*x(4))) KQ_phan_bo_deu =(K5PBD+K7PBD+K11PBD+K13PBD); disp('Ham chi phi tinh toan theo phan bo deu'); disp(KQ_phan_bo_deu) Hieu_qua=100*(KQ_phan_bo_deu-fval)/KQ_phan_bo_deu % Tham so cua cac phan tu cua cac bo loc (C: F; L; H) disp('Tham so cua cac phan tu voi phan bo toi uu') Xc5=(h5^2*V1^2/(h5^2-1))/Q5; C5=1/(2*pi*50*Xc5*1000) XL5=Xc5/h5^2; L5=XL5*1000/(2*pi*50) Xc7=(h7^2*V1^2/(h7^2-1))/Q7; C7=1/(2*pi*50*Xc7*1000) XL7=Xc7/h7^2; L7=XL7*1000/(2*pi*50) Xc11=(h11^2*V1^2/(h11^2-1))/Q11; C11=1/(2*pi*50*Xc11*1000) XL11=Xc11/h11^2; L11=XL11*1000/(2*pi*50) Xc13=(h13^2*V1^2/(h13^2-1))/Q13; C13=1/(2*pi*50*Xc13*1000) XL13=Xc13/h13^2; L13=XL13*1000/(2*pi*50) disp('Tham so cua cac phan tu voi phan bo deu') Xc5=(h5^2*V1^2/(h5^2-1))/x(1); C5=1/(2*pi*50*Xc5*1000) XL5=Xc5/h5^2; L5=XL5*1000/(2*pi*50) Xc7=(h7^2*V1^2/(h7^2-1))/x(2); C7=1/(2*pi*50*Xc7*1000) 89 XL7=Xc7/h7^2; L7=XL7*1000/(2*pi*50) Xc11=(h11^2*V1^2/(h11^2-1))/x(3); C11=1/(2*pi*50*Xc11*1000) XL11=Xc11/h11^2; L11=XL11*1000/(2*pi*50) Xc13=(h13^2*V1^2/(h13^2-1))/x(4); C13=1/(2*pi*50*Xc13*1000) XL13=Xc13/h13^2; L13=XL13*1000/(2*pi*50) Formatted: Font: 10 pt Code 3Phụ lục 3: Các câu lệnh Matlab 2011a tính chi phí, thông số cho nhóm lọc sóng hài thụ Formatted: Font: 10 pt, Italic động đặt sơ cấp MBA sơ đồ lưới điện nhà máy II công ty TNHH Điện Stanley Việt Nam Formatted: Font: 10 pt Formatted: Font: 10 pt clc clear all close all format short Formatted: Font: 10 pt % Truong hop bo loc dat tai so cap may bien ap phan phoi; Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt 90 Formatted Formatted Formatted Formatted Formatted Formatted Formatted % Xet tung bo loc dat tai so cap may bien ap - Nhap gia tri tuong ung tung nhanh; Formatted Formatted % Cong suat dinh muc cua may bien ap phan phoi; Formatted Sr = input('Nhap cong suat dinh muc cua may bien ap Sr (kVA): '); Formatted Formatted % Cong suat thuc cua may bien ap co bo loc dat o so cap; Formatted S = input('Nhap cong suat thuc cua may bien ap S(kVA): '); Formatted Formatted % He so ton hao phu o tan so co ban Formatted X = 0.15; % Lay bang 0.15 Formatted Formatted % He so tai song hai Formatted LF = 0.33; % Gia tri la gia tri gia su Formatted % Don gia dien nang cong thuc tinh ton hao dien nang quan may bien ap Formatted Formatted Formatted UE = 0.1; Formatted Formatted % Dien ap dinh muc cua nhanh so cap may bien ap; Formatted V1 = input('Nhap vao dien ap cao ap cua may bien ap lap bo loc V1(kV): '); Formatted Formatted % Dong dien so cap co ban cua may bien ap lap dat bo loc Formatted I1 = input('Nhap dong dien so cap co ban cua may bien ap I1(A): '); Formatted Formatted % Bien cac phan song hai cua nhanh lap dat bo loc Formatted I5 = input('Nhap vao dong dien hai bac cua nhanh I5(A): '); Formatted Formatted I7 = input('Nhap vao dong dien hai bac cua nhanh I7(A): '); Formatted I11 = input('Nhap vao dong dien hai bac 11 cua nhanh I11(A): '); Formatted Formatted I13 = input('Nhap vao dong dien hai bac 13 cua nhanh I13(A): '); Formatted Formatted % Tong luong CSPK yeu cau cua nhanh lap dat bo loc Formatted Qyeucau = input('Nhap vao CSPK yeu cau cua nhanh Qyeucau(kVAr): '); Formatted Formatted % -Bac song hai can loai tru h5=4.75; h7=6.65; h11=10.45; h13=12.35; Formatted Formatted Formatted Formatted % -Chi phi cho tu va khang MC=1; % Chi phi co dinh cua tu dien UC=1; % He so chi phi phu thuoc ML=45; % Chi phi co dinh cua cuon khang UL=45; % He so chi phi phu thuoc cua cuon khang Formatted Formatted Formatted Formatted 91 Formatted Formatted Formatted Formatted Formatted Formatted Formatted Formatted Formatted Formatted % -He so anpha cua tung bac song hai - Formatted Anpha5 = 1.046; Anpha7 = 1.023; Anpha11 = 1.009; Anpha13 = 1.007; Formatted Formatted Formatted Formatted % -He so chat luong cua bo loc Q=40; Formatted Formatted % -He so ton hao cong suat bo tu kc=0.0004; Formatted Formatted % He so quy doi ve thoi diem hien tai Pv=12.46 % Tinh vong doi bo loc la 20 nam voi lai suat la 0.05 Formatted % -Don gia chi phi ton hao qua tu C va dien tro R cua bo loc Formatted Formatted Formatted c=0.1 Formatted Formatted % Ton hao cong suat cuon day dinh muc Formatted PR = 0.003*Sr Formatted Formatted % He so K - factor theo cong thuc luan van Formatted Kfactor = (h5*I5^2+h7*I7+h11*I11^2+h13*I13^2)/(I5^2+I7^2+I11^2+I13^2) Formatted Formatted % Ton hao cong suat cuon day Formatted PT=PR*S/Sr*(1+X*(Kfactor))/(1+X)*((I5+I7+I11+I13)/I1)^2 Formatted Formatted % Chi phi ton hao song hai bo loc dat o so cap may bien ap KL=8760*Pv*PT*LF*UE Formatted % Chi phi ton that qua may bien ap phan phoi Formatted %===========================================================% Thiet lap ham fun=@(x)(MC+ML)+UC*3*Anpha5*V1*sqrt(I5^2*(1+1/h5^2+((3*I5*V1)/(h5*x(1)))^ 2)+(x(1)/(3*V1))^2)+ 3*UL*Anpha5*V1/h5^2*sqrt(I5^2*(1+h5^2+((3*h5*I5*V1)/(x(1)))^2)+(x(1)/(3*V 1))^2)+ (Pv*8760*kc*x(1)*c)+(Pv*8760*c*((x(1)^2/(3*V1^2))+I5^2)*h5*V1^2/((h5^21)*Q*x(1)))+ (MC+ML)+UC*3*Anpha7*V1*sqrt(I7^2*(1+1/h7^2+((3*I7*V1)/(h7*x(2)))^2)+(x(2) /(3*V1))^2)+ 3*UL*Anpha7*V1/h7^2*sqrt(I7^2*(1+h7^2+((3*h7*I7*V1)/(x(2)))^2)+(x(2)/(3*V 1))^2)+ (Pv*8760*kc*x(2)*c)+(Pv*8760*c*((x(2)^2/(3*V1^2))+I7^2)*h7*V1^2/((h7^21)*Q*x(2)))+ (MC+ML)+UC*3*Anpha11*V1*sqrt(I11^2*(1+1/h11^2+((3*I11*V1)/(h11*x(3)))^2)+ (x(3)/(3*V1))^2)+ 3*UL*Anpha11*V1/h11^2*sqrt(I11^2*(1+h11^2+((3*h11*I11*V1)/(x(3)))^2)+(x(3 )/(3*V1))^2)+ (Pv*8760*kc*x(3)*c)+(Pv*8760*c*((x(3)^2/(3*V1^2))+I11^2)*h11*V1^2/((h11^2 -1)*Q*x(3)))+ Formatted Formatted Formatted Formatted Formatted Formatted Formatted Formatted Formatted Formatted Formatted Formatted Formatted Formatted Formatted Formatted 92 (MC+ML)+UC*3*Anpha13*V1*sqrt(I13^2*(1+1/h13^2+((3*I13*V1)/(h13*x(4)))^2)+ (x(4)/(3*V1))^2)+ 3*UL*Anpha13*V1/h13^2*sqrt(I13^2*(1+h13^2+((3*h13*I13*V1)/(x(4)))^2)+(x(4 )/(3*V1))^2)+ (Pv*8760*kc*x(4)*c)+(Pv*8760*c*((x(4)^2/(3*V1^2))+I13^2)*h13*V1^2/((h13^2 -1)*Q*x(4)))+KL Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt % Cac rang buoc constr=@(x)sum(x); Aeq = [1 1 1]; beq = Qyeucau; lb = zeros(1,4); Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt % Nghiem khoi dau va cac lua chon tinh x0 = [10;10;10;10]; options = optimset('Algorithm','interior-point'); Formatted: Font: 10 pt % Lenh tinh toan tim cuc tieu Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt [x, fval] = fmincon(fun,x0,[],[],Aeq,beq,lb,[],[],options); Formatted: Font: 10 pt Formatted: Font: 10 pt % Hien thi ket qua Formatted: Font: 10 pt disp('Cong suat cua moi bo loc theo tinh toan phan bo toi uu (kVar)') disp(['x = ']);disp(x) Formatted: Font: 10 pt Formatted: Font: 10 pt % Don vi cua cong suat la kVAr Formatted: Font: 10 pt disp('Ham chi phi tinh toan theo phan bo toi uu'); disp(fval) Formatted: Font: 10 pt Formatted: Font: 10 pt Q5=x(1); Q7=x(2); Q11=x(3); Q13=x(4); Formatted: Font: 10 pt disp('Chi phi tinh toan cho tung bo loc theo phan bo toi uu') Formatted: Font: 10 pt K5TU=(MC+ML)+UC*3*Anpha5*V1*sqrt(I5^2*(1+1/h5^2+((3*I5*V1)/(h5*x(1)))^2)+ (x(1)/(3*V1))^2)+ 3*UL*Anpha5*V1/h5^2*sqrt(I5^2*(1+h5^2+((3*h5*I5*V1)/(x(1)))^2)+(x(1)/(3*V 1))^2)+ (Pv*8760*kc*x(1)*c)+(Pv*8760*c*((x(1)^2/(3*V1^2))+I5^2)*h5*V1^2/((h5^21)*Q*x(1))) Formatted: Font: 10 pt K7TU=(MC+ML)+UC*3*Anpha7*V1*sqrt(I7^2*(1+1/h7^2+((3*I7*V1)/(h7*x(2)))^2)+ (x(2)/(3*V1))^2)+ 3*UL*Anpha7*V1/h7^2*sqrt(I7^2*(1+h7^2+((3*h7*I7*V1)/(x(2)))^2)+(x(2)/(3*V 1))^2)+ (Pv*8760*kc*x(2)*c)+(Pv*8760*c*((x(2)^2/(3*V1^2))+I7^2)*h7*V1^2/((h7^21)*Q*x(2))) Formatted: Font: 10 pt K11TU= (MC+ML)+UC*3*Anpha11*V1*sqrt(I11^2*(1+1/h11^2+((3*I11*V1)/(h11*x(3)))^2)+ (x(3)/(3*V1))^2)+ 3*UL*Anpha11*V1/h11^2*sqrt(I11^2*(1+h11^2+((3*h11*I11*V1)/(x(3)))^2)+(x(3 )/(3*V1))^2)+ (Pv*8760*kc*x(3)*c)+(Pv*8760*c*((x(3)^2/(3*V1^2))+I11^2)*h11*V1^2/((h11^2 -1)*Q*x(3))) K13TU=(MC+ML)+UC*3*Anpha13*V1*sqrt(I13^2*(1+1/h13^2+((3*I13*V1)/(h13*x(4) ))^2)+(x(4)/(3*V1))^2)+ Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt 93 3*UL*Anpha13*V1/h13^2*sqrt(I13^2*(1+h13^2+((3*h13*I13*V1)/(x(4)))^2)+(x(4 )/(3*V1))^2)+ (Pv*8760*kc*x(4)*c)+(Pv*8760*c*((x(4)^2/(3*V1^2))+I13^2)*h13*V1^2/((h13^2 -1)*Q*x(4))) Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt % Tinh toan voi truong hop luong CSPK duoc chia deu cho cac bo loc Formatted: Font: 10 pt x(1)=Qyeucau/4; x(2)=x(1); x(3)=x(1); x(4)=x(1); Formatted: Font: 10 pt Formatted: Font: 10 pt disp('Chi phi tinh toan cho tung bo loc theo phan bo deu') K5PBD=(MC+ML)+UC*3*Anpha5*V1*sqrt(I5^2*(1+1/h5^2+((3*I5*V1)/(h5*x(1)))^2) +(x(1)/(3*V1))^2)+ 3*UL*Anpha5*V1/h5^2*sqrt(I5^2*(1+h5^2+((3*h5*I5*V1)/(x(1)))^2)+(x(1)/(3*V 1))^2)+ (Pv*8760*kc*x(1)*c)+(Pv*8760*c*((x(1)^2/(3*V1^2))+I5^2)*h5*V1^2/((h5^21)*Q*x(1))) Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt K7PBD=(MC+ML)+UC*3*Anpha7*V1*sqrt(I7^2*(1+1/h7^2+((3*I7*V1)/(h7*x(2)))^2) +(x(2)/(3*V1))^2)+ 3*UL*Anpha7*V1/h7^2*sqrt(I7^2*(1+h7^2+((3*h7*I7*V1)/(x(2)))^2)+(x(2)/(3*V 1))^2)+ (Pv*8760*kc*x(2)*c)+(Pv*8760*c*((x(2)^2/(3*V1^2))+I7^2)*h7*V1^2/((h7^21)*Q*x(2))) Formatted: Font: 10 pt K11PBD= (MC+ML)+UC*3*Anpha11*V1*sqrt(I11^2*(1+1/h11^2+((3*I11*V1)/(h11*x(3)))^2)+ (x(3)/(3*V1))^2)+ 3*UL*Anpha11*V1/h11^2*sqrt(I11^2*(1+h11^2+((3*h11*I11*V1)/(x(3)))^2)+(x(3 )/(3*V1))^2)+ (Pv*8760*kc*x(3)*c)+(Pv*8760*c*((x(3)^2/(3*V1^2))+I11^2)*h11*V1^2/((h11^2 -1)*Q*x(3))) Formatted: Font: 10 pt K13PBD=(MC+ML)+UC*3*Anpha13*V1*sqrt(I13^2*(1+1/h13^2+((3*I13*V1)/(h13*x(4 )))^2)+(x(4)/(3*V1))^2)+ 3*UL*Anpha13*V1/h13^2*sqrt(I13^2*(1+h13^2+((3*h13*I13*V1)/(x(4)))^2)+(x(4 )/(3*V1))^2)+ (Pv*8760*kc*x(4)*c)+(Pv*8760*c*((x(4)^2/(3*V1^2))+I13^2)*h13*V1^2/((h13^2 -1)*Q*x(4))) Formatted: Font: 10 pt KQ_phan_bo_deu =(K5PBD+K7PBD+K11PBD+K13PBD); disp('Ham chi phi tinh toan theo phan bo deu'); disp(KQ_phan_bo_deu) Hieu_qua=100*(KQ_phan_bo_deu-fval)/KQ_phan_bo_deu Formatted: Font: 10 pt % Tham so cua cac phan tu cua cac bo loc (C: F; L; H) Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt disp('Tham so cua cac phan tu voi phan bo toi uu') Xc5=(h5^2*V1^2/(h5^2-1))/Q5; C5=1/(2*pi*50*Xc5*1000) XL5=Xc5/h5^2; L5=XL5*1000/(2*pi*50) Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Xc7=(h7^2*V1^2/(h7^2-1))/Q7; C7=1/(2*pi*50*Xc7*1000) XL7=Xc7/h7^2; L7=XL7*1000/(2*pi*50) Formatted: Font: 10 pt Formatted: Font: 10 pt Xc11=(h11^2*V1^2/(h11^2-1))/Q11; C11=1/(2*pi*50*Xc11*1000) XL11=Xc11/h11^2; L11=XL11*1000/(2*pi*50) Formatted: Font: 10 pt Xc13=(h13^2*V1^2/(h13^2-1))/Q13; C13=1/(2*pi*50*Xc13*1000) XL13=Xc13/h13^2; L13=XL13*1000/(2*pi*50) Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt 94 Formatted: Font: 10 pt Formatted: Font: 10 pt disp('Tham so cua cac phan tu voi phan bo deu') Xc5=(h5^2*V1^2/(h5^2-1))/x(1); C5=1/(2*pi*50*Xc5*1000) XL5=Xc5/h5^2; L5=XL5*1000/(2*pi*50) Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Xc7=(h7^2*V1^2/(h7^2-1))/x(2); C7=1/(2*pi*50*Xc7*1000) XL7=Xc7/h7^2; L7=XL7*1000/(2*pi*50) Formatted: Font: 10 pt Formatted: Font: 10 pt Xc11=(h11^2*V1^2/(h11^2-1))/x(3); C11=1/(2*pi*50*Xc11*1000) XL11=Xc11/h11^2; L11=XL11*1000/(2*pi*50) Formatted: Font: 10 pt Formatted: Font: 10 pt Xc13=(h13^2*V1^2/(h13^2-1))/x(4); C13=1/(2*pi*50*Xc13*1000) XL13=Xc13/h13^2; L13=XL13*1000/(2*pi*50) Formatted: Font: 10 pt disp('Chi phi ton hao song hai qua may bien ap') disp(KL) Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt clc clear all close all format short Formatted: Font: 10 pt Formatted: Font: 10 pt Formatted: Font: 10 pt % Truong hop bo loc dat tai so cap may bien ap phan phoi; % Xet tung bo loc dat tai so cap may bien ap - Nhap gia tri tuong ung tung nhanh; % Cong suat dinh muc cua may bien ap phan phoi; Sr = input('Nhap cong suat dinh muc cua may bien ap Sr (kVA): '); % Cong suat thuc cua may bien ap co bo loc dat o so cap; S = input('Nhap cong suat thuc cua may bien ap S(kVA): '); % He so ton hao phu o tan so co ban X = 0.15; % Lay bang 0.15 % He so tai song hai LF = 0.33; % Gia tri la gia tri gia su % Don gia dien nang cong thuc tinh ton hao dien nang quan may bien ap UE = 0.1; % Dien ap dinh muc cua nhanh so cap may bien ap; V1 = input('Nhap vao dien ap cao ap cua may bien ap lap bo loc V1(kV): '); % Dong dien so cap co ban cua may bien ap lap dat bo loc I1 = input('Nhap dong dien so cap co ban cua may bien ap I1(A): '); 95 % Bien cac phan song hai cua nhanh lap dat bo loc I5 = input('Nhap vao dong dien hai bac cua nhanh I5(A): '); I7 = input('Nhap vao dong dien hai bac cua nhanh I7(A): '); I11 = input('Nhap vao dong dien hai bac 11 cua nhanh I11(A): '); I13 = input('Nhap vao dong dien hai bac 13 cua nhanh I13(A): '); % Tong luong CSPK yeu cau cua nhanh lap dat bo loc Qyeucau = input('Nhap vao CSPK yeu cau cua nhanh Qyeucau(kVAr): '); % -Bac song hai can loai tru h5=4.75; h7=6.65; h11=10.45; h13=12.35; % -Chi phi cho tu va khang MC=1; % Chi phi co dinh cua tu dien UC=1; % He so chi phi phu thuoc ML=45; % Chi phi co dinh cua cuon khang UL=45; % He so chi phi phu thuoc cua cuon khang % -He so anpha cua tung bac song hai Anpha5 = 1.046; Anpha7 = 1.023; Anpha11 = 1.009; Anpha13 = 1.007; % -He so chat luong cua bo loc Q=40; % -He so ton hao cong suat bo tu kc=0.0004; % He so quy doi ve thoi diem hien tai Pv=12.46 % Tinh vong doi bo loc la 20 nam voi lai suat la 0.05 % -Don gia chi phi ton hao qua tu C va dien tro R cua bo loc c=0.1 % Ton hao cong suat cuon day dinh muc PR = 0.003*Sr % He so K - factor theo cong thuc luan van Kfactor = (h5*I5^2+h7*I7+h11*I11^2+h13*I13^2)/(I5^2+I7^2+I11^2+I13^2) % Ton hao cong suat cuon day PT=PR*S/Sr*(1+X*(Kfactor))/(1+X)*((I5+I7+I11+I13)/I1)^2 96 % Chi phi ton hao song hai bo loc dat o so cap may bien ap KL=8760*Pv*PT*LF*UE % Chi phi ton that qua may bien ap phan phoi %===========================================================% Thiet lap ham fun=@(x)(MC+ML)+UC*3*Anpha5*V1*sqrt(I5^2*(1+1/h5^2+((3*I5*V1)/(h5*x(1)))^ 2)+(x(1)/(3*V1))^2)+ 3*UL*Anpha5*V1/h5^2*sqrt(I5^2*(1+h5^2+((3*h5*I5*V1)/(x(1)))^2)+(x(1)/(3*V 1))^2)+ (Pv*8760*kc*x(1)*c)+(Pv*8760*c*((x(1)^2/(3*V1^2))+I5^2)*h5*V1^2/((h5^21)*Q*x(1)))+ (MC+ML)+UC*3*Anpha7*V1*sqrt(I7^2*(1+1/h7^2+((3*I7*V1)/(h7*x(2)))^2)+(x(2) /(3*V1))^2)+ 3*UL*Anpha7*V1/h7^2*sqrt(I7^2*(1+h7^2+((3*h7*I7*V1)/(x(2)))^2)+(x(2)/(3*V 1))^2)+ (Pv*8760*kc*x(2)*c)+(Pv*8760*c*((x(2)^2/(3*V1^2))+I7^2)*h7*V1^2/((h7^21)*Q*x(2)))+ (MC+ML)+UC*3*Anpha11*V1*sqrt(I11^2*(1+1/h11^2+((3*I11*V1)/(h11*x(3)))^2)+ (x(3)/(3*V1))^2)+ 3*UL*Anpha11*V1/h11^2*sqrt(I11^2*(1+h11^2+((3*h11*I11*V1)/(x(3)))^2)+(x(3 )/(3*V1))^2)+ (Pv*8760*kc*x(3)*c)+(Pv*8760*c*((x(3)^2/(3*V1^2))+I11^2)*h11*V1^2/((h11^2 -1)*Q*x(3)))+ (MC+ML)+UC*3*Anpha13*V1*sqrt(I13^2*(1+1/h13^2+((3*I13*V1)/(h13*x(4)))^2)+ (x(4)/(3*V1))^2)+ 3*UL*Anpha13*V1/h13^2*sqrt(I13^2*(1+h13^2+((3*h13*I13*V1)/(x(4)))^2)+(x(4 )/(3*V1))^2)+ (Pv*8760*kc*x(4)*c)+(Pv*8760*c*((x(4)^2/(3*V1^2))+I13^2)*h13*V1^2/((h13^2 -1)*Q*x(4)))+KL % Cac rang buoc constr=@(x)sum(x); Aeq = [1 1 1]; beq = Qyeucau; lb = zeros(1,4); % Nghiem khoi dau va cac lua chon tinh x0 = [10;10;10;10]; options = optimset('Algorithm','interior-point'); % Lenh tinh toan tim cuc tieu [x, fval] = fmincon(fun,x0,[],[],Aeq,beq,lb,[],[],options); % Hien thi ket qua disp('Cong suat cua moi bo loc theo tinh toan phan bo toi uu (kVar)') disp(['x = ']);disp(x) % Don vi cua cong suat la kVAr disp('Ham chi phi tinh toan theo phan bo toi uu'); disp(fval) Q5=x(1); Q7=x(2); Q11=x(3); Q13=x(4); disp('Chi phi tinh toan cho tung bo loc theo phan bo toi uu') 97 K5TU=(MC+ML)+UC*3*Anpha5*V1*sqrt(I5^2*(1+1/h5^2+((3*I5*V1)/(h5*x(1)))^2)+ (x(1)/(3*V1))^2)+ 3*UL*Anpha5*V1/h5^2*sqrt(I5^2*(1+h5^2+((3*h5*I5*V1)/(x(1)))^2)+(x(1)/(3*V 1))^2)+ (Pv*8760*kc*x(1)*c)+(Pv*8760*c*((x(1)^2/(3*V1^2))+I5^2)*h5*V1^2/((h5^21)*Q*x(1))) K7TU=(MC+ML)+UC*3*Anpha7*V1*sqrt(I7^2*(1+1/h7^2+((3*I7*V1)/(h7*x(2)))^2)+ (x(2)/(3*V1))^2)+ 3*UL*Anpha7*V1/h7^2*sqrt(I7^2*(1+h7^2+((3*h7*I7*V1)/(x(2)))^2)+(x(2)/(3*V 1))^2)+ (Pv*8760*kc*x(2)*c)+(Pv*8760*c*((x(2)^2/(3*V1^2))+I7^2)*h7*V1^2/((h7^21)*Q*x(2))) K11TU= (MC+ML)+UC*3*Anpha11*V1*sqrt(I11^2*(1+1/h11^2+((3*I11*V1)/(h11*x(3)))^2)+ (x(3)/(3*V1))^2)+ 3*UL*Anpha11*V1/h11^2*sqrt(I11^2*(1+h11^2+((3*h11*I11*V1)/(x(3)))^2)+(x(3 )/(3*V1))^2)+ (Pv*8760*kc*x(3)*c)+(Pv*8760*c*((x(3)^2/(3*V1^2))+I11^2)*h11*V1^2/((h11^2 -1)*Q*x(3))) K13TU=(MC+ML)+UC*3*Anpha13*V1*sqrt(I13^2*(1+1/h13^2+((3*I13*V1)/(h13*x(4) ))^2)+(x(4)/(3*V1))^2)+ 3*UL*Anpha13*V1/h13^2*sqrt(I13^2*(1+h13^2+((3*h13*I13*V1)/(x(4)))^2)+(x(4 )/(3*V1))^2)+ (Pv*8760*kc*x(4)*c)+(Pv*8760*c*((x(4)^2/(3*V1^2))+I13^2)*h13*V1^2/((h13^2 -1)*Q*x(4))) % Tinh toan voi truong hop luong CSPK duoc chia deu cho cac bo loc x(1)=Qyeucau/4; x(2)=x(1); x(3)=x(1); x(4)=x(1); disp('Chi phi tinh toan cho tung bo loc theo phan bo deu') K5PBD=(MC+ML)+UC*3*Anpha5*V1*sqrt(I5^2*(1+1/h5^2+((3*I5*V1)/(h5*x(1)))^2) +(x(1)/(3*V1))^2)+ 3*UL*Anpha5*V1/h5^2*sqrt(I5^2*(1+h5^2+((3*h5*I5*V1)/(x(1)))^2)+(x(1)/(3*V 1))^2)+ (Pv*8760*kc*x(1)*c)+(Pv*8760*c*((x(1)^2/(3*V1^2))+I5^2)*h5*V1^2/((h5^21)*Q*x(1))) K7PBD=(MC+ML)+UC*3*Anpha7*V1*sqrt(I7^2*(1+1/h7^2+((3*I7*V1)/(h7*x(2)))^2) +(x(2)/(3*V1))^2)+ 3*UL*Anpha7*V1/h7^2*sqrt(I7^2*(1+h7^2+((3*h7*I7*V1)/(x(2)))^2)+(x(2)/(3*V 1))^2)+ (Pv*8760*kc*x(2)*c)+(Pv*8760*c*((x(2)^2/(3*V1^2))+I7^2)*h7*V1^2/((h7^21)*Q*x(2))) K11PBD= (MC+ML)+UC*3*Anpha11*V1*sqrt(I11^2*(1+1/h11^2+((3*I11*V1)/(h11*x(3)))^2)+ (x(3)/(3*V1))^2)+ 3*UL*Anpha11*V1/h11^2*sqrt(I11^2*(1+h11^2+((3*h11*I11*V1)/(x(3)))^2)+(x(3 )/(3*V1))^2)+ (Pv*8760*kc*x(3)*c)+(Pv*8760*c*((x(3)^2/(3*V1^2))+I11^2)*h11*V1^2/((h11^2 -1)*Q*x(3))) 98 K13PBD=(MC+ML)+UC*3*Anpha13*V1*sqrt(I13^2*(1+1/h13^2+((3*I13*V1)/(h13*x(4 )))^2)+(x(4)/(3*V1))^2)+ 3*UL*Anpha13*V1/h13^2*sqrt(I13^2*(1+h13^2+((3*h13*I13*V1)/(x(4)))^2)+(x(4 )/(3*V1))^2)+ (Pv*8760*kc*x(4)*c)+(Pv*8760*c*((x(4)^2/(3*V1^2))+I13^2)*h13*V1^2/((h13^2 -1)*Q*x(4))) KQ_phan_bo_deu =(K5PBD+K7PBD+K11PBD+K13PBD); disp('Ham chi phi tinh toan theo phan bo deu'); disp(KQ_phan_bo_deu) Hieu_qua=100*(KQ_phan_bo_deu-fval)/KQ_phan_bo_deu % Tham so cua cac phan tu cua cac bo loc (C: F; L; H) disp('Tham so cua cac phan tu voi phan bo toi uu') Xc5=(h5^2*V1^2/(h5^2-1))/Q5; C5=1/(2*pi*50*Xc5*1000) XL5=Xc5/h5^2; L5=XL5*1000/(2*pi*50) Xc7=(h7^2*V1^2/(h7^2-1))/Q7; C7=1/(2*pi*50*Xc7*1000) XL7=Xc7/h7^2; L7=XL7*1000/(2*pi*50) Xc11=(h11^2*V1^2/(h11^2-1))/Q11; C11=1/(2*pi*50*Xc11*1000) XL11=Xc11/h11^2; L11=XL11*1000/(2*pi*50) Xc13=(h13^2*V1^2/(h13^2-1))/Q13; C13=1/(2*pi*50*Xc13*1000) XL13=Xc13/h13^2; L13=XL13*1000/(2*pi*50) disp('Tham so cua cac phan tu voi phan bo deu') Xc5=(h5^2*V1^2/(h5^2-1))/x(1); C5=1/(2*pi*50*Xc5*1000) XL5=Xc5/h5^2; L5=XL5*1000/(2*pi*50) Xc7=(h7^2*V1^2/(h7^2-1))/x(2); C7=1/(2*pi*50*Xc7*1000) XL7=Xc7/h7^2; L7=XL7*1000/(2*pi*50) Xc11=(h11^2*V1^2/(h11^2-1))/x(3); C11=1/(2*pi*50*Xc11*1000) XL11=Xc11/h11^2; L11=XL11*1000/(2*pi*50) Xc13=(h13^2*V1^2/(h13^2-1))/x(4); C13=1/(2*pi*50*Xc13*1000) XL13=Xc13/h13^2; L13=XL13*1000/(2*pi*50) disp('Chi phi ton hao song hai qua may bien ap') disp(KL) 99 ... tần số cộng hưởng cho lọc cộng hưởng đơn a Lý không nên chọn tần số cộng hưởng tần số sóng hài Bộ lọc sóng hài thụ động đƣợc điều chỉnh để hút sóng hài tần số đƣợc tính toán, tần số sóng hài. .. Tổng quan lọc sóng hài thụ động Bộ lọc sóng hài đƣợc sử dụng phổ biến công nghiệp giá thành rẻ dễ bảo dƣỡng, sửa chữa Các lọc sóng hài thụ động đƣợc cấu thành từ tụ kháng với trị số thích hợp... c Bộ lọc thông cao bậc d Bộ lọc thông cao bậc e Bộ lọc thông cao kiểu C Nguyên lý làm việc lọc cộng hƣởng đơn: Bộ lọc cộng hƣởng đơn lọc có tụ điện kháng điện đấu nối tiếp với nhau, t Thông số