Xác định giá điện thời thực

ĐẠI HỌC QUỐC GIA THÀNH PHỐ HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA CAO MINH THUẦN XÁC ĐỊNH GIÁ ĐIỆN THỜI THỰC CHUYÊN NGÀNH: THIẾT BỊ, MẠNG VÀ NHÀ MÁY ĐIỆN LUẬN VĂN THẠC SĨ TP HỒ CHÍ MINH, THÁNG 07 NĂM -i- CƠNG TRÌNH ĐƯỢC HỒN THÀNH TẠI TRƯỜNG ĐẠI HỌC BÁCH KHOA ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH Cán hướng dẫn khoa học : PGS.TS Phan Thị Thanh Bình Cán chấm nhận xét : Cán chấm nhận xét : Luận văn thạc sĩ bảo vệ HỘI ĐỒNG CHẤM BẢO VỆ LUẬN VĂN THẠC SĨ TRƯỜNG ĐẠI HỌC BÁCH KHOA, ngày tháng năm 2008 - ii - ĐẠI HỌC QUỐC GIA TP HCM TRƯỜNG ĐẠI HỌC BÁCH KHOA - CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc lập – Tự – Hạnh Phúc o0o Tp HCM, ngày tháng năm 2008 NHIỆM VỤ LUẬN VĂN THẠC SĨ X / Nữ Họ tên học viên:…… Cao Minh Thuần Giới tính: Nam Ngày, tháng, năm sinh:… 22 – 02 – 1975 Nơi sinh: Hậu giang Chuyên ngành:… Thiết bị, mạng nhà máy điện .MSHV: 01806507 Khoá (Năm trúng tuyển) : ….2006 I- TÊN ĐỀ TÀI: …… XÁC ĐỊNH GIÁ ĐIỆN THỜI THỰC II- NHIỆM VỤ LUẬN VĂN: Giới thiệu giá điện Nguyên tắc việc hình thành giá nút Phương pháp xác định chương trình mô giá điện nút Kết luận III- NGÀY GIAO NHIỆM VỤ : …….21/01/2008 IV- NGÀY HOÀN THÀNH NHIỆM VỤ:… 30/06/2008 V- CÁN BỘ HƯỚNG DẪN: PGS.TS PHAN THỊ THANH BÌNH Nội dung đề cương luận văn thạc sĩ Hội đồng chuyên ngành thông qua Ngày …… tháng …… năm 2008 CÁN BỘ HƯỚNG DẪN (Họ tên chữ ký) PGS.TS PHAN THỊ THANH BÌNH CHỦ NHIỆM BỘ MÔN QUẢN LÝ CHUYÊN NGHÀNH (Họ tên chữ ký) - iii - LỜI CẢM ƠN Trước hết, xin trân trọng bày tỏ lòng biết ơn chân thành đến thầy cô môn Hệ Thống Điện thuộc khoa Kỹ thuật điện – Điện tử trường đại học Bách khoa TP HCM, dạy dỗ suốt thời gian học cao học, đặc biệt bày tỏ lòng cảm ơn sâu sắc đến cô Phan Thị Thanh Bình, người tận tình hướng dẫn trực tiếp để hoàn thành luận văn Tôi xin gửi lời cám ơn bạn bè học lớp cao học khoá 2006, đồng nghiệp làm việc trường Cao Đẳng Điện Lực TP HCM, có nhiều ý kiến thiết thực giúp hoàn thành tốt đề tài Cho phép gửi lời cám ơn đến Ban giám hiệu trường Cao Đẳng Điện Lực TP HCM tạo điều kiện để tham gia lớp học Cuối xin cảm ơn đến người mẹ cố người cha vó đại cho sức mạnh tinh thần để vượt qua khó khăn công việc học tập TP HCM, tháng 07 năm 2008 Cao Minh Thuần - iv - TĨM TẮT LUẬN VĂN Mục đích đề tài để cung cấp mô tả chi tiết giá điện nút cách chia nhỏ giá nút thành nhiều phần khác dựa yếu tố liên quan như: nguồn phát, nghẽn mạch, giới hạn điện áp, ràng buộc khác Các thông tin giá nút dùng không để nâng cao hiệu hoạt động lưới điện quản lý, vận hành hệ thống, mà để tạo khung giá hợp lý cho hệ thống điện khuyến khích nhà đầu tư vào nguồn phát đường dây truyền tải điện Trên sở xây dựng chương trình mơ xác định giá điện nút Thiết kế giao diện chương trình Matlab, thân thiện với người sử dụng thông qua giao diện trực quan, dễ sử dụng -v- MỤC LỤC Giới thiệu đề tài 01 Chương I: Giới thiệu giá điện 06 I./ Giá cố định 06 1/ Quyết định giá doanh nghiệp có độc quyền định 06 2/ Quyết định giá có hai loại doanh nghiệp 07 II./ Giá nút 09 1/ Mục tiêu việc hình thành giá nút 11 2/ Sự phản hồi khách hàng 15 3/ Vận hành thị trường điện nước phát triển 18 4/ Vận hành thị trường điện nước phát triển 20 Chương II: Các nguyên tắc việc hình thành giá nút 23 I./ Công thức phương pháp thành lập 23 1/ Hàm mục tiêu 25 2/ Các ràng buộc cân 25 3/ Các ràng buộc không cân 26 II./ Giải vấn đề tối ưu phân bố công suất với (P, Q) 27 1/ Phương trình Newton 28 2/ Phương pháp qui hoạch tuyến tính 32 3/ Phương pháp hàm Lagrangian 36 Chương III: Xây dựng chương trình mô giá điện nút 39 I./ Xác định giá điện nút 39 1/ Hình thành giá điện nút 39 2/ Ví dụ 40 II./ Chương trình mơ giá điện nút 44 1/ Giới thiệu Matpower 44 - vi - 2/ Một số file M Matpower 45 III/ Xây dựng giao điện chương trình 46 1/ Giới thiệu giao diện 46 2/ Các kết 48 3/ Thay đổi số liệu từ giao diện 52 4/ Đánh giá kết chương trình 53 5/ Kết luận 53 IV./ Kết tính tốn phân tích kết 54 1./ Lưu đồ chương trình 54 2./ Giải thích giải thuật chương trình 55 3./ Sơ đồ lưới điện 56 4./ Số liệu đầu vào toán 58 5./ Kết chương trình 65 6./ Phân tích kết 68 Chương IV: Kết luận hướng phát triển đề tài 77 I/ Kết luận 77 II/ Hướng phát triển đề tài 78 Tài liệu tham khảo 80 Phụ lục PL-01 - vii - CÁC TỪ VIẾT TẮT Vm: biên độ điện áp Va: góc pha điện áp Vg: điện áp máy phát (cài đặt) Vmax: điện áp cực đại Vmin: điện áp cực tiểu Pg: công suất tác dụng máy phát Qg: công suất phản kháng máy phát Pmax: công suất tác dụng cực đại máy phát Qmax: công suất phản kháng cực đại máy phát Pmin: công suất tác dụng cực tiểu máy phát Qmin: công suất phản kháng cực tiểu máy phát r: điện trở dây dẫn x: điện kháng dây dẫn b=(y/2): dung dẫn đường dây Gs, điện kháng song song Bs, điện dung song song bus: nút gen: máy phát điện cost: chi phí branch: nhánh (từ nút đến nút…) gencost: hàm chi phí máy phát ratio: tỉ số máy biến áp angle: góc lệch pha máy biến áp status: tình trạng làm việc đường dây, máy biến áp fbus: từ nút… tbus: đến nút… -1- GIỚI THIỆU ĐỀ TÀI I Đặt vấn đề Ngành điện ngày trở thành sở hạ tầng kinh tế, mức tiêu thụ trở thành thước đo trình độ kinh tế quốc gia Để khắc phục tồn tại, đáp ứng thách thức đặt su hướng tồn cầu hố, thực thắng lợi cơng cơng nghiệp hố, đại hố đất nước Ngành điện cần phải tiếp tục xếp, đổi mới, phát triển nâng cao hiệu hoạt động sản xuất, kinh doanh, đảm bảo giá cả, độ tin cậy, chất lượng, quyền tự lựa chọn khách hàng đồng thời có khả cạnh tranh với công ty điện lực khu vực giới vấn đề thiết yếu Hiện thị trường điện hình thành phát triển nhiều nước giới Thị trường điện với chế mở, cạnh tranh hoạt động có hiệu nước cho thấy ưu điểm mặc kinh tế kỹ thuật hẳn hệ thống điện tập trung cấu theo chiều dọc truyền thống Thu nhận kết từ nước chuyển sang thị trường điện như: Anh, Mỹ, Canada, Brasil, Úc, Nhật, Trung Quốc, Thái Lan cho thấy hệ thống điện không ngừng phát triển không số lượng, chất lượng mà giá bán điện cho người sử dụng rẻ Thực tế giới, công ty điện lực trải qua quãng thời gian với nhiều thay đổi, đặc biệt cấu trúc thị trường điện sách kèm theo Với xuất nhà sản xuất điện độc lập (IPPs) thay đổi cấu trúc hệ thống cung cấp điện (Electricity Supply Industry – ESI), ngành điện bước vào thời kỳ phát triển cạnh tranh với yêu cầu phải chứng minh hiệu kinh tế độ tin cậy áp lực thị trường Lý thuyết giá nút (Nodal Price) hệ thống điện với cấu trúc phân quyền để đưa dấu hiệu kinh tế xác để sử dụng nguồn tài nguyên lưới nguồn hiệu Khác với cấu trúc ngành dọc truyền thống, giá nút thay đổi theo thời gian khơng gian Nó gồm chi phí cho việc sản xuất vận hành chi phí khác để trì độ tin cung cấp điện dịch vụ kèm theo Đến thời điểm có số phương pháp để tính tốn giá điện nút đưa khơng Thực hiện: Cao Minh Thuần ***** GVHD: PGS TS Phan Thị Thanh Bình -2- dựa vào chi phí vận hành ngắn hạn mà dịch vụ phụ khác (Ancillary Service) Hầu hết lý thuyết giá nút dựa vào hệ số phương trình Lagrangian để đánh giá giá trị yếu tố ràng buộc, độ tin cậy chất lượng điện Mặc dù hệ số Lagrangian cung cấp thêm thông tin để phản ánh ràng buộc (nghẽn mạch, giới hạn điện áp…), khơng thể cung cấp cách chi tiết mô tả giá nút Nói cách khác, khơng có rõ ràng ảnh hưởng giá nút nhà máy điện ràng buộc khác, không rõ ràng việc yếu tố ảnh hưởng đến giá nút Do để tìm hiểu kỹ thành phần cấu thành giá nút, có nhiều phương pháp đưa Trong luận án trình bày phương pháp xác định giá điện nút chia nhỏ thành thành phần độc lập tuyến tính, khơng máy nhà phát điện mà cịn có ràng buộc khác tồn hệ thống điện Mục đích việc khảo sát đưa phương pháp để mơ tả chi tiết giá nút Ví dụ chia nhỏ thành thành phần liên quan đến yếu tố nguồn phát, nghẽn mạch, giới hạn điện áp ràng buộc khác Sự phân tích thành phần nút đứng quan điểm kinh tế đạo hàm thành phần lấy dựa vào điều kiện biên Việc phân tích thông tin đầy đủ không giúp cho việc sử dụng tài nguyên lưới điện quản lý nghẽn mạch hiệu mà cịn có vai trị khuyến khích nhà đầu tư đầu tư vào ngành điện Xuất phát từ yêu cầu thực tế trên, đề tài tập trung nghiên cứu thị trường điện, cách thức hình thành giá điện nút hệ thống điện thơng qua việc giải tốn phân bố công suất hệ thống điện II Nhiệm vụ luận văn Dựa vào phần đặt vấn đề trên, nhiệm vụ luận văn đặt là: - Giới thiệu loại giá điện Thực hiện: Cao Minh Thuần ***** GVHD: PGS TS Phan Thị Thanh Bình PL-26 mpc = s; if ~isfield(mpc, 'version') mpc.version = '1'; end if strcmp(mpc.version, '1') [mpc.gen, mpc.branch] = mpc_1to2(mpc.gen, mpc.branch); mpc.version = '2'; end end end %% - Dinh nghia gia tr xuat if return_as_struct bus = info; elseif ~expect_opf_data areas = info; end - if info == %% Khong co loi if return_as_struct baseMVA = mpc; else baseMVA = mpc.baseMVA; bus = mpc.bus; gen = mpc.gen; branch = mpc.branch; if isfield(mpc, 'gencost') areas = mpc.areas; gencost = mpc.gencost; end end else %% Co loi cau truc if nargout == | nargout == | nargout == %% tra ve ma loi if return_as_struct baseMVA = struct([]); else baseMVA = []; bus = []; gen = []; branch = []; if expect_opf_data areas = []; gencost = []; end end else %% Thong bao loi switch info case 1, error('loadcase: Tham so vao se la cau truc hoac la chuoi filename'); case 2, error('loadcase: Xac dinh truong hop duong dan khon o MATLAB''s '); case 3, error('loadcase: Xac dinh file Mat khong ton tai'); case 4, error('loadcase: Xac dinh file M khong ton tai'); Thực hiện: Cao Minh Thuần ***** GVHD: PGS TS Phan Thị Thanh Bình PL-27 case 5, error(sprintf('loadcase: Loi cau truc or Khong xac dinh duoc ma tran so lieu file\n%s', err5)); otherwise, error('loadcase: Khong nhan dang duoc loi'); end end end return; %% - mpc_1to2 function [gen, branch] = mpc_1to2(gen, branch) %% Dinh nghia ten ma ran bus, gen, branch [GEN_BUS, PG, QG, QMAX, QMIN, VG, GEN_STATUS, PMAX, PMIN, MU_PMAX, MU_PMIN, MU_QMAX, MU_QMIN, PC1, PC2, QC1MIN, QC1MAX, QC2MIN, QC2MAX, RAMP_AGC, RAMP_10, RAMP_30, RAMP_Q, APF] = idxgen; [F_BUS, BR_R, BR_X, BR_B, RATE_A, RATE_B, RATE_C, TAP, SHIFT, BR_STATUS, PF, QF, PT, QT, MU_SF, MU_ST, ANGMIN, ANGMAX, MU_ANGMIN, MU_ANGMAX] = idxbrch; %% - gen %% shift = MU_PMAX - PMIN - 1; tmp = num2cell([MU_PMAX, MU_PMIN, MU_QMAX, MU_QMIN] - shift); [MU_PMAX, MU_PMIN, MU_QMAX, MU_QMIN] = drual(tmp{:}); %% Them cot phu vao may phat tmp = zeros(size(gen, 1), shift); if size(gen, 2) >= MU_QMIN gen = [ gen(:, 1:PMIN) tmp gen(:, MU_PMAX:MU_QMIN) ]; else gen = [ gen(:, 1:PMIN) tmp ]; end %% - Nhanh shift = PF - BR_STATUS - 1; tmp = num2cell([PF, QF, PT, QT, MU_SF, MU_ST] - shift); [PF, QF, PT, QT, MU_SF, MU_ST] = deal(tmp{:}); %% Them cot phu vao nhanh tmp = [ ones(size(branch, 1), 1) * [-360 360] ]; tmp2 = [ zeros(size(branch, 1), 2) ]; if size(branch, 2) >= MU_ST branch = [ branch(:, 1:BR_STATUS) tmp branch(:, PF:MU_ST) tmp2 ]; elseif size(branch, 2) >= QT branch = [ branch(:, 1:BR_STATUS) tmp branch(:, PF:QT) ]; else branch = [ branch(:, 1:BR_STATUS) tmp ]; end return; %% - FMINCON Thực hiện: Cao Minh Thuần ***** GVHD: PGS TS Phan Thị Thanh Bình PL-28 %% tim cuc tieu cac rang buoc cua cac ham so function [X,FVAL,EXITFLAG,OUTPUT,LAMBDA,GRAD,HESSIAN] = fmincon(FUN,X,A,B,Aeq,Beq,LB,UB,NONLCON,options,varargin) defaultopt = struct('Display','final','LargeScale','on', 'TolX',1e-6,'TolFun',1e-6,'TolCon',1e6,'DerivativeCheck','off', 'Diagnostics','off','FunValCheck','off', 'GradObj','off','GradConstr','off', 'HessMult',[], % HessMult [] mac dinh 'Hessian','off','HessPattern','sparse(ones(numberOfVariables))', 'MaxFunEvals','100*numberOfVariables', 'MaxSQPIter','10*max(numberOfVariables,numberOfInequalities+number OfBounds)', 'DiffMaxChange',1e-1,'DiffMinChange',1e-8, 'PrecondBandWidth',0,'TypicalX','ones(numberOfVariables,1)', 'MaxPCGIter','max(1,floor(numberOfVariables/2))', 'TolPCG',0.1,'MaxIter',400,'OutputFcn',[],'PlotFcns',[], 'RelLineSrchBnd',[],'RelLineSrchBndDuration',1,'NoStopIfFlatInfeas ','off', 'PhaseOneTotalScaling','off'); if nargin==1 && nargout computeLambda = 1; else computeLambda = 0; end caller='constr'; lenVarIn = length(varargin); XOUT=X(:); numberOfVariables=length(XOUT); %check for empty X if numberOfVariables == error('optim:fmincon:EmptyX','Can cung cap diem khoi dau khong rong.'); end switch optimget(options,'Display',defaultopt,'fast') case {'off','none'} verbosity = 0; case 'notify' verbosity = 1; case 'final' verbosity = 2; case 'iter' verbosity = 3; otherwise verbosity = 2; end % Set to column vectors B = B(:); Beq = Beq(:); % Tim giai thuat line_search = strcmp(optimget(options,'LargeScale',defaultopt,'fast'),'off'); [XOUT,l,u,msg] = checkbounds(XOUT,LB,UB,numberOfVariables); Thực hiện: Cao Minh Thuần ***** GVHD: PGS TS Phan Thị Thanh Bình PL-30 if ~isempty(msg) EXITFLAG = -2; [FVAL,LAMBDA,GRAD,HESSIAN] = deal([]); OUTPUT.iterations = 0; OUTPUT.funcCount = 0; OUTPUT.stepsize = []; if line_search OUTPUT.lssteplength = []; end if line_search OUTPUT.algorithm = 'medium-scale: SQP, Quasi-Newton, linesearch'; else OUTPUT.algorithm = 'large-scale: trust-region reflective Newton'; end OUTPUT.firstorderopt = []; if ~line_search OUTPUT.cgiterations = []; end OUTPUT.message = msg; X(:)=XOUT; if verbosity > disp(msg) end return end lFinite = l(~isinf(l)); uFinite = u(~isinf(u)); meritFunctionType = 0; mtxmpy = optimget(options,'HessMult',defaultopt,'fast'); if isequal(mtxmpy,'hmult') warning('optim:fmincon:HessMultNameClash', ['Ten ham dien ap dung voi ham tro giup Toolbox :\n', ' su dung ten canh ''hmult'' cho ham HessMult den\n', ' tranh loi hoac ket qua khong mong doi.']); end diagnostics = isequal(optimget(options,'Diagnostics',defaultopt,'fast'),'on'); funValCheck = strcmp(optimget(options,'FunValCheck',defaultopt,'fast'),'on'); gradflag = strcmp(optimget(options,'GradObj',defaultopt,'fast'),'on'); hessflag = strcmp(optimget(options,'Hessian',defaultopt,'fast'),'on'); if isempty(NONLCON) constflag = 0; else constflag = 1; Thực hiện: Cao Minh Thuần ***** GVHD: PGS TS Phan Thị Thanh Bình PL-31 end if ~isempty(FUN) % se kim tra chuoi rong, ma tran rong, mang rong [funfcn, msg] = optimfcnchk(FUN,'fmincon',length(varargin),funValCheck,gradflag,he ssflag); else error('optim:fmincon:InvalidFUN', ['FUN phai la ham can thiet;\n', ' hoac, FUN co the la mot mang ma rang buoc la can thiet.']); end if constflag % NONLCON la khong rong gradconstflag = strcmp(optimget(options,'GradConstr',defaultopt,'fast'),'on'); [confcn, msg] = optimfcnchk(NONLCON,'fmincon',length(varargin),funValCheck,gradcon stflag,false,1); else gradconstflag = false; confcn{1} = ''; end [rowAeq,colAeq]=size(Aeq); % Neu chi co l va u goi sfminbx if ~line_search && isempty(NONLCON) && isempty(A) && isempty(Aeq) && gradflag OUTPUT.algorithm = large; % Neu chi co Aeq beq va Aeq co so cot hon so hang, sfminle elseif ~line_search && isempty(NONLCON) && isempty(A) && isempty(lFinite) && isempty(uFinite) && gradflag && colAeq > rowAeq OUTPUT.algorithm = large; elseif ~line_search warning('optim:fmincon:SwitchingToMediumScale', ['phuong phap Large-scale khong su dung van de nay,\n' ' thay vi su dung medium-scale.']) if isequal(funfcn{1},'fungradhess') funfcn{1}='fungrad'; warning('optim:fmincon:HessianIgnored', ['Phuong pha Medium-scale la phuong phap Quasi-Newton va khong duoc su dung\n' 'phan tich Hessian Ky hieu Hessian lua chon se duoc bo qua.']) elseif isequal(funfcn{1},'fun_then_grad_then_hess') funfcn{1}='fun_then_grad'; warning('optim:fmincon:HessianIgnored', ['Phuong pha Medium-scale la phuong phap Quasi-Newton va khong duoc su dung\n' 'phan tich Hessian Ky hieu Hessian lua chon se duoc bo qua.']) end hessflag = 0; OUTPUT.algorithm = medium; Thực hiện: Cao Minh Thuần ***** GVHD: PGS TS Phan Thị Thanh Bình PL-32 elseif line_search OUTPUT.algorithm = medium; if issparse(Aeq) || issparse(A) warning('optim:fmincon:ConvertingToFull', 'Khong su dung ma tran vuong voi phuong phap mediumscale: Chuyen doi hoan toan.') end if line_search && hessflag % xung dot cac lua chon hessflag = 0; warning('optim:fmincon:HessianIgnored', ['Phuong phap Medium-scale la phuong phap Quasi-Newton kkhong dung phan tich Hessian.\n' 'Ky hieu Hessian lua chon se duoc bo qua (Ap dung Hessian se khong duoc su dung).']); if isequal(funfcn{1},'fungradhess') funfcn{1}='fungrad'; elseif isequal(funfcn{1},'fun_then_grad_then_hess') funfcn{1}='fun_then_grad'; end end % nguoc lai goi nlconst else error('optim:fmincon:InvalidOptions', 'Khong nhan su ket hop lua chon OPTIONS goi tuan tu.') end lenvlb=length(l); lenvub=length(u); if isequal(OUTPUT.algorithm,medium) % % Dam bao rang diem bat dau nam tren cac bien % i=1:lenvlb; lindex = XOUT(i)u(i); if any(uindex) XOUT(uindex)=u(uindex); end X(:) = XOUT; else % % Neu x ban dau khong o tren bien, thiet lap dan diem "boxcentered" % arg = (u >= 1e10); arg2 = (l % Phat hien thong tin msg = diagnose('fmincon',OUTPUT,gradflag,hessflag,constflag,gradconstfla g, line_search,options,defaultopt,XOUT,non_eq, non_ineq,lin_eq,lin_ineq,l,u,funfcn,confcn,f,GRAD,HESS,c,ceq,cGRAD ,ceqGRAD); end % Goi gia thuat if isequal(OUTPUT.algorithm,medium) [X,FVAL,lambda,EXITFLAG,OUTPUT,GRAD,HESSIAN]= nlconst(funfcn,X,l,u,full(A),B,full(Aeq),Beq,confcn,options,defaul topt, verbosity,gradflag,gradconstflag,hessflag,meritFunctionType, f,GRAD,HESS,c,ceq,cGRAD,ceqGRAD,varargin{:}); LAMBDA=lambda; else if (isequal(funfcn{1}, 'fun_then_grad_then_hess') || isequal(funfcn{1}, 'fungradhess')) Hstr=[]; elseif (isequal(funfcn{1}, 'fun_then_grad') || isequal(funfcn{1}, 'fungrad')) n = length(XOUT); Hstr = optimget(options,'HessPattern',defaultopt,'fast'); Thực hiện: Cao Minh Thuần ***** GVHD: PGS TS Phan Thị Thanh Bình PL-37 if ischar(Hstr) if isequal(lower(Hstr),'sparse(ones(numberofvariables))') Hstr = sparse(ones(n)); else error('optim:fmincon:InvalidHessPattern', 'Chon lua ''HessPattern'' phai la ma tran neu khong duoc mac dinh.') end end end if isempty(Aeq) [X,FVAL,LAMBDA,EXITFLAG,OUTPUT,GRAD,HESSIAN] = sfminbx(funfcn,X,l,u,verbosity,options,defaultopt,computeLambda,f, GRAD,HESS,Hstr,varargin{:}); else [X,FVAL,LAMBDA,EXITFLAG,OUTPUT,GRAD,HESSIAN] = sfminle(funfcn,X,sparse(Aeq),Beq,verbosity,options,defaultopt,comp uteLambda,f,GRAD,HESS,Hstr,varargin{:}); end end %% - MAKEYBUS - function [Ybus, Yf, Yt] = makeYbus(baseMVA, bus, branch) %% cac hang so j = sqrt(-1); nb = size(bus, 1); %% So nut nl = size(branch, 1); %% So nhanh %% Dinh nghia dua vao ma tran bus, branch [PQ, PV, REF, NONE, BUS_I, BUS_TYPE, PD, QD, GS, BS, BUS_AREA, VM, VA, BASE_KV, ZONE, VMAX, VMIN, LAM_P, LAM_Q, MU_VMAX, MU_VMIN] = idx_bus; [F_BUS, T_BUS, BR_R, BR_X, BR_B, RATE_A, RATE_B, RATE_C, TAP, SHIFT, BR_STATUS, PF, QF, PT, QT, MU_SF, MU_ST, ANGMIN, ANGMAX, MU_ANGMIN, MU_ANGMAX] = idx_brch; %% Kiem tra so nut co bang giua bang so lieu va nut hay khong (chi co so nut ma thoi) if any(bus(:, BUS_I) ~= [1:nb]') error('buses must appear in order by bus number') end %% moi nhanh,tinh phan cua nhanh ma tran tong dan %% %% | If | | Yff Yft | | Vf | %% | | = | | * | | %% | It | | Ytf Ytt | | Vt | %% stat = branch(:, BR_STATUS); %% ones at inservice branches Thực hiện: Cao Minh Thuần ***** GVHD: PGS TS Phan Thị Thanh Bình PL-38 Ys = stat / (branch(:, BR_R) + j * branch(:, BR_X)); %% tong dan noi tiep Bc = stat * branch(:, BR_B); %% Tong dan duong day tap = ones(nl, 1); %% Mac dinh ti so bien ap = i = find(branch(:, TAP)); %% Dinh nghia ti so MBA khac khong tap(i) = branch(i, TAP); %% Gan ti so MBA khac khong tap = tap * exp(j*pi/180 * branch(:, SHIFT)); %% Them su dich pha vao Ytt = Ys + j*Bc/2; Yff = Ytt / (tap * conj(tap)); Yft = - Ys / conj(tap); Ytf = - Ys / tap; %% Tinh tong dan song song %% Neu Psh la CS tac dung tieu hao noi song song o V = 1.0 p.u %% va Qsh la CS dua vao song song o V = 1.0 p.u %% Psh - j Qsh = V * conj(Ysh * V) = conj(Ysh) = Gs - j Bs, %% i.e Ysh = Psh + j Qsh, so Ysh = (bus(:, GS) + j * bus(:, BS)) / baseMVA; %% Vector cua ma tran tong dan %% Xay dung Ybus f = branch(:, F_BUS); %% Liet ke tu "from" nut t = branch(:, T_BUS); %% Liet ke tu "to" nut Cf = sparse(f, 1:nl, ones(nl, 1), nb, nl); %% Noi ma tran cho duong day & from nut Ct = sparse(t, 1:nl, ones(nl, 1), nb, nl); %% Noi ma tran cho duong day & to nut Ybus = spdiags(Ysh, 0, nb, nb) + %% Tong dan song song Cf * spdiags(Yff, 0, nl, nl) * Cf' + %% Yff dinh nghia tong dan nhanh Cf * spdiags(Yft, 0, nl, nl) * Ct' + %% Yft dinh nghia tong dan nhanh Ct * spdiags(Ytf, 0, nl, nl) * Cf' + %% Ytf dinh nghia tong dan nhanh Ct * spdiags(Ytt, 0, nl, nl) * Ct'; %% Ytt dinh nghia tong dan nhanh %% Xay dung Yf va Yt nhu la Yf * dua vao %% O moi nhanh "from" nut, va Yt if nargout > i = [[1:nl]'; [1:nl]']; Yf = sparse(i, [f; t], [Yff; Yt = sparse(i, [f; t], [Ytf; end return; Thực hiện: Cao Minh Thuần V la vector dong nhanh phuc duoc giong nhu "to" nut cuoi Yft], nl, nb); Ytt], nl, nb); ***** GVHD: PGS TS Phan Thị Thanh Bình PL-39 %% - VE DO THI -function Do_thi_ieee30 [MVAbase, bus, gen, gencost, branch, f, success, et ,nut,volM,volA,sumPG,sumQG,sumPL,sumQL,lamP,lamQ]=runopf('case_iee e30'); figure plot(nut,lamP) xlabel('Nut','fontsize',12,'fontweight','b','color','magenta'); ylabel('Gia dien','fontsize',12,'fontweight','b','color','magenta'); grid on end Thực hiện: Cao Minh Thuần ***** GVHD: PGS TS Phan Thị Thanh Bình PL-40 TĨM TẮT LÝ LỊCH TRÍCH NGANG Họ tên học viên:…Cao Minh Thuần Phái:……Nam Ngày, tháng, năm sinh:… 22 – 02 – 1975 Nơi sinh: …Hậu giang Địa liên hệ: Trường Cao Đẳng Điện Lực TP HCM, số 455 Hà Huy Giáp, P Thạnh Lộc, Quận 12, TP HCM Quá trình đào tạo: - Từ năm 1993 – 1998, học Khoa điện, trường Đại Học Sư Phạm Kỹ Thuật TP HCM - Từ năm 1996 – 1998, học cao học chuyên ngành Thiết Bị, Mạng Và Nhà Máy Điện, khoa điện điện tử trường Đại Học Bách Khoa TP HCM Quá trình công tác: - Từ năm 1998 – 2003: làm việc công ty Xây Lắp Điện - Từ năm 2003 đến nay: công tác khoa Lưới điện, trường Cao Đẳng Điện Lực Thành Phố Hồ Chí Minh Thực hiện: Cao Minh Thuần ***** GVHD: PGS TS Phan Thị Thanh Bình ... kim loại thực với nguyên tắc giá – lượng thấp (nghĩa xác xuất bị cắt điện cao giá thấp hơn) Trong thời điểm có khác biệt lớn giá thời điểm có xác xuất cắt điện thấp có xác xuất cắt điện cao,... Chương I: Giới thiệu giá điện I./ Giá cố định 1/ Quyết định giá doanh nghiệp có độc quyền định 2/ Quyết định giá có hai loại doanh nghiệp II./ Giá nút 1/ Mục tiêu việc hình thành giá nút 2/ Sự phản... trình mơ giá điện nút I./ Xác định giá điện nút 1/ Hình thành giá điện nút 2/ Ví dụ II./ Chương trình mơ giá điện nút 1/ Giới thiệu Matpower 2/ Một số file M Matpower III/ Xây dựng giao điện chương