分布式电源接入对配电网的影响matlab程序（IEEE9节点系统算例）

分布式电源接入对配电网的影响matlab程序（IEEE9节点系统算例）

摘 要：分布式电源的接入使得配电系统从放射状无源网络变为分布有中小型电源的有源网络。带来了使单向流动的电流方向具有了不确定性等等问题，使得配电系统的控制和管理变得更加复杂。但同时，分布式电源又具有提高电网可靠性，绿色节能等优点，所以为更好的利用分布式电源为人类造福，我们必须对其进行研究与分析。本文利用仿真软件Matlab编写计算潮流程序模拟分布式电源接入配电网的模型进行潮流计算的方法对分布式电源的稳态影响进行探索与分析。选取了9节点的配电网网络模型，通过对单个分布式电源的接入位置以及容量的不同情况对9节点配电网的网损以及节点电压状况进行了分析。
关键字：分布式电源，配电网，牛顿拉夫逊法

1 配电网潮流计算方法
研究其影响，就是要研究其电能质量等的各项指标，所以采取何种潮流计算方法是第一步要做的工作。自60年代稀疏矩阵技术应用于牛顿法以来，经过几十年发展，已经成为求解电力系统潮流问题的最广泛的一种方法。当以节点功率为注入量时，潮流方程为一组线性方程，牛顿法为求解非线性组最有效的方法之一。牛顿法的极坐标方程为：

对公式进行泰勒展开，取一次项，即可得到牛顿-拉夫逊潮流算法的修正方程组。

式中：，为潮流方程的残差向量，,为母线的电压修正量，J为雅克比矩阵。
.。。。。。。。。潮流计算理论省略
详细 参考：https://blog.csdn.net/weixin_47365903/article/details/120895286?spm=1001.2014.3001.5502

2 IEEE9节点系统算例
在4节点接入有功60Mw，无功60Mvar分布式电源
2.1 系统结构

2.2 系统节点参数

2.3 系统变压器信息

2.4 支路参数信息

3 程序运行结果
1）潮流计算结果
接入前

接入后

2）电压分分布结果

4 matlab程序

clc
clear
close all
%%运行潮流
t0 = clock;
%% 列索引
[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;
[GEN_BUS, PG, QG, QMAX, QMIN, VG, MBASE, 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] = idx_gen;


casename = 'case9'; 
mpopt = mpoption;       %% use default options
%% 打印选项
verbose = mpopt(31);
%% 1.初始潮流
mpc = loadcase(casename);
mpc.branch
[baseMVA, bus, gen, branch] = deal(mpc.baseMVA, mpc.bus, mpc.gen, mpc.branch);
[ref, pv, pq] = bustypes(bus, gen);
[Ybus, Yf, Yt] = makeYbus(baseMVA, bus, branch);
Sbus = makeSbus(baseMVA, bus, gen);
V0    = ones(size(bus, 1), 1);   
[V, success, iterations] = newtonpf(Ybus, Sbus, V0, ref, pv, pq, mpopt);
V1=sqrt(abs(V).^2);
figure
plot(V1,'-*')
[bus, gen, branch] = pfsoln(baseMVA, bus, gen, branch, Ybus, Yf, Yt, V, ref, pv, pq);
mpc.et = etime(clock, t0);
mpc.success = success;
%% -----  输出结果  -----
[mpc.bus, mpc.gen, mpc.branch] = deal(bus, gen, branch);
printpf(mpc, 1, mpopt);



%% 2.加入分布式电源后潮流
% 读数据
fprintf('请输入分布式电源的节点位置，有功与无功大小    bus  Pg   Qg  ') 
DGs2=[0  0  1  0  1  0  0  0  0  0  0  0  0  0  0  0  0  0]
DGbus=input('请输入DG节点位置的值：');%4
DGPg=input('请输入DG有功大小的值：');%30
DGQg=input('请输入DG无功大小的值：');%30
DGs1=[DGbus DGPg DGQg] 
DGs=[DGs1,DGs2]
mpc.gen=[mpc.gen;DGs]
。。。。。。。略
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47