Matlab 코드: 2차 원뿔 프로그래밍에 기반한 능동 분배 네트워크를 위한 최적의 전력 흐름 솔루션
참고: 활성 배전망의 다중 소스 협력 운영 최적화에 관한 연구_차오샨
초록: 최적의 전력 흐름 연구는 배전망 계획 및 운영에 있어 필수불가결하며, 특히 다수의 분산형 에너지원이 연결되는 능동형 배전망 환경에서 중요하다. 기존의 휴리스틱 알고리즘은 글로벌 최적 솔루션 및 솔루션 속도 측면에서 활성 배전망의 운영 요구 사항을 충족할 수 없으며, 선형화를 기반으로 한 최적 전력 흐름 방식도 고임피던스 배전망에서 취약합니다. 이를 기반으로 본 논문은 2차 콘 프로그래밍 기반의 동적 최적 전력 흐름 모델 프레임워크를 구축하여 빠른 솔루션을 위해 원래의 비선형 프로그래밍 모델의 완화를 SOPC로 변환하려고 시도합니다. 먼저 2차 콘 이완을 기반으로 한 배전망 동적 최적 전력 흐름의 기본 모델을 제시하고 이를 효율적으로 해결하기 위해 능동 배전망에 참여하는 중요 요소 각각에 대해 해당하는 선형 모델링 처리를 수행하며, 능동관리장비, 배전망 재구성, 수요반응 및 통합부하 등의 완화모델의 정확도와 근사등가물을 동시에 분석한다.
프로그램 플롯:
코드의 일부:
%Multi-period+SVC+CB+OLTC+DG SOCP_OPF Sbase=1MVA, Ubase=12.66KV
%목적 함수에 네트워크 손실만 있는 경우 OLTC는 항상 고속 기어에 있게 됩니다.전압이 높을수록 네트워크는 작아집니다. 메인 네트워크 전원 구매 또는 전압 밸런스와 같은 기능
%%
%부하 시 탭 체인저의 위치는 해당 노드에 있습니다.
%%
clear
clc
tic
warning off
%% 1. 매개변수 설정
mpc = IEEE33BW;
wind = mpc.wind;
pload = mpc.pload;
pload_prim = mpc.pload_prim/1000
; ;
a = 3.715; 모든 노드의 활성 용량 % 단일 기간, MW
b = 2.3, % 단일 기간에 있는 모든 노드의 무효 용량, MW
pload = pload/a, % 비례 계수 얻기
qload = pload/b, % 활성 부하 곡선이 다음과 같다고 가정합니다. 반응 부하 곡선
pload = pload_prim*pload;%는 33*24의 부하 값을 얻습니다, 각 기간에 각 노드의 부하 qload
= qload_prim*qload;
branch = mpc.branch;
branch(:,3) = branch(:,3)*1/(12.66^2);% 단위 임피던스 값 계산
R = real(branch(:,3));
X = imag( branch(:,3));
T = 24;% 기간은 24시간
nb = 33;% 노드 번호
nl = 32;% 분기 번호
nsvc = 3;%SVC 번호 Static Var 보상기
ncb = 2;%CB 번호 그룹 스위칭 커패시터 뱅크(커패시터뱅크, CB)
noltc = 1;%OLTC 부하시 탭 체인저 번호, OLTC 변압기
nwt = 2;%2 팬
ness = 2;% ESS 업스트림 번호
= 0(nb,nl);
dnstream = zeros(nb,nl);
for i = 1:nl
upstream(i,i)=1;
end
for i = [1:16,18:20,22: 23,25:31]
dnstream(i,i+1 )=1;
끝