一种基于改进PSO的主动配电网可再生DG规划方法与流程

一种基于改进pso的主动配电网可再生dg规划方法
技术领域
1.本发明涉及电力系统分布式电源规划领域，尤其是涉及一种基于改进pso的主动配电网可再生dg规划方法。


背景技术：

2.能源对人类的生存与社会的发展来说是不可或缺的重要部分，随着电力需求的持续增长、传统能源紧缺形式的不断加深以及环境问题的日益突出，分布式电源(dg，distributed generation)尤其是可再生能源发电技术的发展获得了广泛支持。分布式电源具有能源利用效率高、清洁环保、安装地点灵活等多方面的优点，可以有效地解决目前电网中成本高、运行难度大等许多潜在问题。主动配电网(active distribution network，adn)以其灵活、兼容以及优化等特性，是未来智能配电网实现对大量接入的分布式电源进行主动管理的有效解决方案。目前，国内外学者对分布式电源接入主动配电网进行规划已有一定成果。
3.中国专利cn201910086985.7公开了一种含分布式电源的配电网分层规划方法，利用psco优化算法得到最终网架结构和dg接入位置与容量配置，提高了供电可靠性。但是，由于可再生dg的输出功率有波动性，目前没有合适的模型或方式能完整描述可再生dg输出功率的不确定性变化，直接将可再生dg简单的当成出力稳定的分布式电源，会忽略其对系统的具体影响，导致与实际的电网运行相差较大。可再生dg参与配网调度运行的需求以及可再生dg出力的间歇性和不确定性给主动配电网的可再生dg协调控制带来了极大的挑战。
4.此外，求解主动配电网的可再生dg规划时，一般通过建立目标函数和约束条件描述主动配电网规划模型，利用数学方法对目标函数进行求解，会导致计算量过大且甚至可能导致不收敛；在利用单一智能算法对目标函数进行求解时，如使用粒子群优化算法(particle swarm optimization，pso)，极有可能陷入局部最优值的问题，且收敛速度和收敛性较差。


技术实现要素：

5.本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种基于改进pso的主动配电网可再生dg规划方法，先确定可再生dg接入主动配电网的较优位置及较优接入容量，再建立主动配电网双层规划模型，并利用改进pso算法结合跟踪中心轨迹内点法对双层规划模型进行求解，得到可再生dg接入主动配电网的最佳接入位置和容量。
6.本发明的目的可以通过以下技术方案来实现：
7.一种基于改进pso的主动配电网可再生dg规划方法，包括以下步骤：
8.s1：对可再生dg的输出功率和主动配电网的负荷进行数学建模，所述可再生dg的输出功率包括风力发电机输出功率和光伏太阳能板输出功率；
9.s2：根据规划区域内的可再生资源总量确定可再生dg的最佳装机容量；
10.s3：根据可再生dg接入主动配电网对系统潮流的影响，通过潮流计算确定可再生
dg接入主动配电网的较优位置及较优接入容量；
11.s4：建立主动配电网上层模型和下层模型，上层模型的目标函数为年综合成本最小，下层模型的目标函数为网络损耗最小，结合负荷的数学模型、可再生dg的输出功率数学模型和可再生dg的最佳装机容量，构建上层模型和下层模型的约束条件；
12.s5：将pso算法与遗传算法相结合并进行自适应改进以求解上层模型，得到可再生dg接入位置和容量信息，并将可再生dg接入位置和容量信息传递到下层模型；
13.s6：通过跟踪中心轨迹内点法求解下层模型，得到可再生dg接入主动配电网的网络损耗并将结果返回给上层模型；
14.s7：计算上层模型是否满足预设置的优化标准，若满足，则输出可再生dg接入位置和容量信息，否则，重复步骤s5。
15.进一步的，所述步骤s1包括以下步骤：
16.s11：对风力发电机的输出功率进行数学建模，如下式所示：
[0017][0018]
其中，v
ci
表示切入风速，单位为m/s；v
r
表示额定风速，单位为m/s；v
co
表示切出风速，单位为m/s；p
r
表示风机额定输出功率，单位为kw；风速数据根据风速的双参数weibull分布确定；
[0019]
s12：对光伏太阳能板的输出功率进行数学建模，如下所示：
[0020][0021]
其中，p
s
_
rated
表示额定输出功率，单位为w；s
r
表示额定光照强度，单位为w/m2；光照强度数据根据光照强度的beta分布确定；
[0022]
s13：对主动配电网的负荷进行数学建模，采用正态分布模型来描述负荷概率曲线，如下所示：
[0023][0024][0025]
其中，p
l,i
、q
l,i
表示节点i有功、无功负荷随机量；μ
p,i
、μ
q,i
表示节点i有功、无功负荷期望；σ
p,i
、σ
q,i
表示节点i有功、无功负荷方差。
[0026]
进一步的，所述步骤s2中，利用区域可再生资源评估方法确定规划区域内可利用的太阳能资源总量，具体为：根据规划区域的经纬度、日照百分率计算可利用的太阳能资源的总量；获取规划区域内太阳能光伏电板的可安装面积；根据太阳能资源的总量和太阳能光伏电板的可安装面积确定最佳光伏装机容量。
[0027]
进一步的，所述步骤s2中，利用区域可再生资源评估方法确定规划区域内可利用的风能资源总量，具体为：根据规划区域的风速数据计算风能资源量化指标；获取规划区域内风机的可安装面积；根据风能资源量化指标和风机的可安装面积确定最佳风机装机容量。
[0028]
进一步的，所述步骤s3包括以下步骤：
[0029]
s31：将相同容量的可再生dg分别接入主动配电网的各个节点，采用牛顿
‑
拉夫逊算法分别计算可再生dg接入各个节点后的潮流参数，根据可再生dg接入主动配电网后引起的网络损耗变化率与所接入节点的负荷占比，确定可再生dg接入主动配电网的较优位置；
[0030]
s32：选取一个节点，分别接入不同容量的可再生dg，采用牛顿
‑
拉夫逊算法分别计算不同容量的可再生dg接入该节点后的潮流参数，根据该节点接入不同容量的可再生dg后的网络损耗值和节点电压，确定该节点的较优接入容量；
[0031]
s33：重复步骤s32，直至得到各个节点的较优接入容量。
[0032]
更进一步的，所述步骤s31中，根据可再生dg接入主动配电网后引起的网络损耗变化率与所接入节点的负荷占比计算综合比率，综合比率较大的节点为较优位置，综合比率的计算公式具体为：
[0033]
z(％)＝k1w+k2f
[0034]
其中，z表示综合比率，w表示可再生dg接入主动配电网后的网络损耗变化率，f表示可再生dg所接入的节点负荷占总负荷的比率，k1、k2表示预设置的权重系数。
[0035]
进一步的，所述步骤s4包括以下步骤：
[0036]
s41：采用蒙特卡洛模拟法根据可再生dg的输出功率数学模型和主动配电网的负荷数学模型生成可再生dg的输出功率数据和主动配电网的负荷数据，获取可再生dg的最佳装机容量；
[0037]
s42：以经济性为目标建立主动配电网建设及运行的上层模型；
[0038]
目标函数为年综合成本最小，如下式所示：
[0039]
c
min
＝c
dg
+c
w
+c
g
+c
l
‑
c
h
[0040]
其中，c
min
表示年综合成本，c
dg
表示可再生dg的投资年费用，c
w
表示可再生dg运行维护费用，c
g
表示向上级电网购电成本，c
l
表示网络损耗费用，c
h
表示政府补贴费用，计算公式分别为：
[0041][0042][0043][0044]
c
l
＝c
k
p
loss
τ
max
[0045][0046]
其中，r表示贴现率，k表示可再生dg使用年限，n表示主动配电网节点数，c
pv,i
表示
第i个节点的单位光伏投资费用，p
pv,i
表示第i个节点的光伏安装量，c
wt,i
表示第i个节点的单位风机投资费用，p
wt,i
表示第i个节点的风机安装量，单位为mw；c
p,i
表示第i个节点的光伏运行维护费用，e
pv,i
表示第i个节点的光伏发电量，c
w,i
表示第i个节点的风机运行维护费用，e
wt,i
表示第i个节点的风机发电量；c
e
表示上级电网购电电价，p
y
表示主动配电网总有功负荷量，t
max
表示年最大负荷利用小时数；c
k
表示网络损耗电价，c
k
表示网络损耗，τ
max
表示年最大负荷损耗小时数；c
b
表示环境补贴费用；
[0047]
约束条件包括：系统节点容量约束、可再生dg安装总容量约束和可再生dg渗透率约束，公式表示依次为：
[0048][0049][0050][0051]
其中，p
dg,i
表示第i个节点可再生dg安装量，k
i
表示主动配电网节点负荷值；p
dg,max
表示主动配电网最大准入容量；
[0052]
s43：以网络损耗最小为目标建立主动配电网建设及运行的下层模型；
[0053]
目标函数为网络损耗最小，如下式所示：
[0054]
f＝min(p
loss
)
[0055]
约束条件包括：节点功率平衡约束、节点电压约束、支路传输功率约束、无功补偿投入装置容量约束、可再生dg出力约束和有载调压变压器分接头调节约束，公式表示依次为：
[0056]
节点功率平衡约束，如下式所示：
[0057][0058]
其中，p
is
表示节点i的有功注入，q
is
表示节点i的无功注入，u
i
、u
j
表示节点i、j电压向量幅值，g
ij
和b
ij
表示电导和电纳；θ
ij
表示节点i和节点j电压相角差；
[0059]
节点电压约束，如下式所示：
[0060][0061]
其中，u
i
表示节点i的电压幅值，n
j
表示系统节点数；u
i,min
和u
i,max
表示节点i电压幅值的下限和上限；
[0062]
支路传输功率约束，如下式所示：
[0063][0064]
其中，s
n
表示支路s
n,max
的传输功率；n
l
表示系统支路数；s
n,max
表示支路s
n,max
传输功
率上限；
[0065]
无功补偿投入装置容量约束，如下式所示：
[0066][0067]
其中，q
ci
表示第i个可再生dg安装节点无功补偿设备投入量，q
ci,min
和q
ci,max
表示无功补偿设备投入量下限和上限；
[0068]
可再生dg出力约束，如下式所示：
[0069][0070]
其中，p
i
表示第i个节点可再生dg有功出力量，p
i,min
和p
i,max
表示可再生dg有功出力量的下限和上限；
[0071]
有载调压变压器分接头调节约束，如下式所示：
[0072]
k
p,min
≤k
p
≤k
p,max
[0073]
其中，k
p
表示有载调压变压器抽头位置，k
p,min
和k
p,max
表示有载调压变压器抽头调节范围的下限与上限。
[0074]
进一步的，所述步骤s5包括以下步骤：
[0075]
s51：根据可再生dg的接入位置生成编码向量作为个体，个体的粒子速度表示可再生dg的容量信息，生成初始种群；
[0076]
s52：设置pso算法的参数值和遗传算法的参数值，对pso算法的惯性权重ω和学习因子c1、c2进行自适应改进，对遗传算法的交叉算子p
cro
和变异算子p
mut
进行自适应改进：
[0077][0078]
其中，ω
max
与ω
min
为最大惯性权重和最小惯性权重，maxiter为最大迭代次数，iter为当前迭代次数；
[0079][0080]
其中，c
max
与c
min
为学习因子的最大和最小值；
[0081][0082]
其中，p
cro
、p
cro1
、p
cro2
表示交叉因子，g1表示常数，f(k+1)表示待交叉个体的适应度函数值，f(avg)表示群体所有个体平均适应度函数值，f(max)表示群体适应度最大个体的适应度函数值；
[0083][0084]
其中，p
mut
、p
mut1
、p
mut2
表示交叉因子，g2表示常数；
[0085]
s53：计算种群中每个个体的适应度，所述适应度是根据上层模型的目标函数值计算得到的，记录该代种群中的最优个体，若种群收敛，则执行步骤s55，否则，执行步骤s54；
[0086]
s54：对种群进行选择操作、交叉操作和变异操作，更新pso算法的惯性权重和学习因子，计算每个个体的粒子速度，并更新每个个体的粒子位置，得到新一代种群，重复步骤s53；
[0087]
s54：输出最优个体所对应的可再生dg接入位置和容量信息，并将其传递到下层模型。
[0088]
更进一步的，所述步骤s51中，生成初始种群包括以下步骤：
[0089]
a1：随机生成编码向量和粒子速度；
[0090]
a2：判断该编码向量和粒子速度对应的可再生dg接入位置和容量信息是否满足上层模型的约束条件，若满足，则将该编码向量和粒子速度作为初始种群中的个体，否则，舍弃该编码向量和粒子速度；
[0091]
a3：重复步骤a1，直至初始种群中的个体数量满足预设置的种群规模。
[0092]
更进一步的，所步骤s53中，种群收敛具体为：最优个体的适应度满足预设置的收敛精度、或在kp(kp>1)次迭代中种群的最优个体的适应度的变化率小于预设置稳定阈值、或种群代数等于预设置的最大迭代次数，kp为预设置的稳定次数。
[0093]
与现有技术相比，本发明具有以下有益效果：
[0094]
(1)先确定可再生dg接入主动配电网的较优位置及较优接入容量，再建立主动配电网双层规划模型，并利用改进pso算法结合跟踪中心轨迹内点法对双层规划模型进行求解，得到可再生dg接入主动配电网的最佳接入位置和容量。
[0095]
(2)将pso算法结合遗传算法求解上层模型，与现有技术中使用单一智能算法求解相对比，算法的可行性更高，能避免陷入局部最优解，提高了算法的搜索性能，并通过对参数进行自适应改进以达到更好的收敛效果。
[0096]
(3)利用区域可再生资源评估方法，综合考虑规划区域的环境因素及地理位置因素确定规划区域内的可再生资源总量，得到的数据能更好地用于确定可再生dg的最佳装机容量，方法合理科学。
[0097]
(4)对可再生dg出力和主动配电网的负荷的不确定性进行数学建模，将可再生dg的输出功率及负荷视为随机变量，并采用蒙特卡洛模拟法对其模拟并描述，降低了可再生dg输出功率不确定性对主动配电网的管理和规划的不利影响。
[0098]
(5)利用牛顿
‑
拉夫逊潮流计算可再生dg的接入对主动配电网潮流的影响，综合考虑网络损耗变化率及所接入节点的负荷占比，确定可再生dg接入主动配电网的较佳位置；再比较将不同容量的可再生dg接入节点对网络损耗及节点电压分布的影响，确定可再生dg接入主动配电网的较佳容量，初步确定了可再生dg接入位置和容量信息，为后续的可再生
dg规划方案求解奠定了基础。
[0099]
(6)考虑可再生dg接入主动配电网的经济性与系统运行稳定性问题，建立以年综合成本最小为目标函数的上层模型，以网络损耗最小为目标函数的下层规划模型，并通过上层模型将初始信息传递给下层模型、下层模型得到的结果反馈给上层模型的方式对双层模型进行求解，得到最优的可再生dg规划方案。
附图说明
[0100]
图1为本发明的流程图；
[0101]
图2为实施例中风机输出功率与风速的关系图；
[0102]
图3为实施例中ieee33节点的配电网图；
[0103]
图4为实施例中可再生dg接入各个节点的网络损耗图；
[0104]
图5为实施例中蒙特卡洛风速概率密度模拟图；
[0105]
图6为实施例中双层模型求解模式示意图；
[0106]
图7为实施例中传统pso算法与改进pso算法的性能比较示意图。
具体实施方式
[0107]
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。
[0108]
实施例1：
[0109]
一种基于改进pso的主动配电网可再生dg规划方法，如图1所示，包括以下步骤：
[0110]
s1：对可再生dg的输出功率和主动配电网的负荷进行数学建模，所述可再生dg的输出功率包括风力发电机输出功率和光伏太阳能板输出功率；
[0111]
s2：根据规划区域内的可再生资源总量确定可再生dg的最佳装机容量；
[0112]
s3：根据可再生dg接入主动配电网对系统潮流的影响，通过潮流计算确定可再生dg接入主动配电网的较优位置及较优接入容量；
[0113]
s4：建立主动配电网上层模型和下层模型，上层模型的目标函数为年综合成本最小，下层模型的目标函数为网络损耗最小，结合负荷的数学模型、可再生dg的输出功率数学模型和可再生dg的最佳装机容量，构建上层模型和下层模型的约束条件；
[0114]
s5：将pso算法与遗传算法相结合并进行自适应改进以求解上层模型，得到可再生dg接入位置和容量信息，并将可再生dg接入位置和容量信息传递到下层模型；
[0115]
s6：通过跟踪中心轨迹内点法求解下层模型，得到可再生dg接入主动配电网的网络损耗并将结果返回给上层模型；
[0116]
s7：计算上层模型是否满足预设置的优化阈值，若满足，则输出可再生dg接入位置和容量信息，否则，重复步骤s5。
[0117]
具体的，可再生dg中，太阳能和风能的开发与利用最为广泛，本实施例中对风力发电机的输出功率、光伏太阳能板的输出功率进行数学建模，步骤s1包括以下步骤：
[0118]
s11：对风力发电机的输出功率进行数学建模：
[0119]
风机发电的输出功率大小与风速有关，风速有很强的不确定性导致风机发电的出
力也有波动性，大量实验结果表明，风速概率密度函数满足双参数weibull分布，其表达式如下所示：
[0120][0121]
其中，v表示风速，单位为m/s；k代表威布尔分布形状参数(无量纲)；c代表威布尔分布尺度参数。k、c可根据风速的均值和方差计算得到，计算公式下：
[0122][0123]
其中，μ表示区域风速统计值的平均值。
[0124][0125]
其中，σ表示区域风速统计值的均方差。
[0126]
使用分段函数表示风机输出有功功率：
[0127][0128]
其中，v
ci
表示切入风速，单位为m/s；v
r
表示额定风速，单位为m/s；v
co
表示切出风速，单位为m/s；p
r
表示风机额定输出功率，单位为kw。
[0129]
风机输出功率与风速的函数关系如图2所示。
[0130]
s12：对光伏太阳能板的输出功率进行数学建模：
[0131]
光伏发电输出的功率会随着光照强度的变化而改变，由于光照强度的不确定性导致光伏太阳能发电输出的功率波动较大，光照强度具有随机性的特点，根据大量实验统计，光照强度可近似认为满足beta分布，其概率密度函数如下所示：
[0132][0133]
其中，s表示光照强度，单位为w/m2；s
max
表示光照强度的最大值，单位为w/m2；α和β表示beta分布的两个形状参数；γ表示伽马函数。
[0134]
beta分布的形状参数由光照强度的平均值μ及其方差σ计算得出：
[0135][0136]
光伏发电输出功率的函数表达式如下所示：
[0137][0138]
其中，p
s
_
rated
表示额定输出功率，单位为w；s
r
表示额定光照强度，单位为w/m2。
[0139]
本实施例中，光伏发电额定光照强度s
r
的大小为1000w/m2，beta分布的形状参数α＝0.89，β＝0.83。
[0140]
s13：对主动配电网的负荷进行数学建模：
[0141]
负荷的波动性会使预测的负荷值与实际负荷值间产生一定误差，主动配电网的负荷如果用确定性方法进行预测描述会造成规划结果与实际存在偏差，本申请采用正态分布模型来描述负荷概率曲线，如下所示：
[0142][0143][0144]
其中，p
l,i
、q
l,i
表示节点i有功、无功负荷随机量；μ
p,i
、μ
q,i
表示节点i有功、无功负荷期望；σ
p,i
、σ
q,i
表示节点i有功、无功负荷方差。
[0145]
对任何一个规划区域来说，区域内可利用的太阳能和风能资源是有限的，如果装机容量超出该区域可利用的太阳能和风能的总量，会减少其经济性。但如果装机容量过小，则导致太阳能和风能资源无法被充分利用，又会造成能源上的浪费。因此，本申请通过区域可再生能源发电资源评估方法对规划区域内太阳能和风能的最大可利用量分析评估，确定与区域资源总量相符合的光伏和风机装机容量，并用于之后上层模型和下层模型中约束条件的构建。
[0146]
步骤s2中，利用区域可再生资源评估方法确定规划区域内可利用的太阳能资源总量，具体为：根据规划区域的经纬度、日照百分率计算可利用的太阳能资源的总量；获取规划区域内太阳能光伏电板的可安装面积；根据太阳能资源的总量和太阳能光伏电板的可安装面积确定最佳光伏装机容量。
[0147]
步骤s2中，利用区域可再生资源评估方法确定规划区域内可利用的风能资源总量，具体为：根据规划区域的风速数据计算风能资源量化指标；获取规划区域内风机的可安装面积；根据风能资源量化指标和风机的可安装面积确定最佳风机装机容量。
[0148]
太阳能光伏板接收的太阳辐射量直接决定了光伏发电量，光照强度与光伏板的发电量具有强相关性，通过matlab仿真平台对太阳能光伏发电搭建simulink仿真模型，利用搭建的光伏模型进行光伏发电仿真。以a市为例，获取其经纬度，规定日照散射比为70％，将0～100
°
范围内的所有角度分别作为光伏太阳能板的倾角，计算出所有角度下的太阳能辐射量，并取年太阳辐射量最大值所在的角度作为光伏最佳安装倾角，再确定规划区域内的光伏安装面积，进而确定光伏板的太阳辐射接收量和发电量，确定最佳光伏装机容量。
[0149]
同理，可以确定最佳风机装机容量。
[0150]
步骤s3中确定可再生dg接入主动配电网的较优位置和较优接入容量。
[0151]
当可再生dg接入主动配电网系统后，不论可再生dg接入哪个节点，系统的网络损耗与电压都发生变化，这是由于当主动配电网接入可再生dg后，潮流流动模式会发生改变，从单向流动转变为双向流动，进而改变系统网络损耗以及节点电压大小，提高主动配电网系统经济性运行。由于可再生dg接入主动配电网后，可以减少主动配电网系统的有功传输，
从而使系统的网络损耗减少，所以接入可再生dg后主动配电网系统的网络损耗会比未接入可再生dg时小，且节点在主动配电网的位置、负荷不同，相同容量的可再生dg接入不同节点带来的网络损耗也不同。
[0152]
主动配电网系统的网络损耗和节点电压与接入可再生dg的容量有关，随着接入容量的改变，系统潮流也发生相应的变化，因为当可再生dg接入主动配电网系统时其发出的电能越多，主动配电网系统传输的电能就越少，因此系统网络损耗也会相应地减少，大大减少了主动配电网系统上电能的损耗。可再生dg接入主动配电网会导致系统负荷减小，流过主动配电网线路的电流会减少，可以减少系统网络损耗；但当可再生dg接入容量大于所接入节点负荷值时，会造成系统反向潮流，系统减小的电流量趋于饱和，甚至增加了系统的电流量，进而增大网络损耗。
[0153]
本实施例中采用ieee33节点配电网为仿真算例，系统电压的等级为12.66kv，系统总有功负荷为3.715mw，系统总无功负荷为2.30mvar。ieee33节点的配电网图如图3所示。通过牛顿
‑
拉夫逊算法进行潮流计算，进而得到可再生dg的较佳接入位置和较佳接入容量，在其他实施方式中，也可以使用其他潮流计算方法，如p
‑
q分解法。
[0154]
首先计算综合比率以确定可再生dg的较佳接入系统位置，综合比率综合考虑了网络损耗变化率与接入节点的负荷占比，综合比率越大，则其对系统潮流影响越大，即为较佳接入点。综合比率计算公式如下：
[0155]
z(％)＝k1w+k2f
[0156]
其中，z表示综合比率，w表示可再生dg接入后系统网络损耗变化率，f表示可再生dg所接入的节点负荷占总负荷的比率，k1、k2表示预设置的权重系数，为常量。
[0157]
其中可再生dg接入后系统网络损耗变化率w与可再生dg所接入的节点负荷占总负荷的比率f的公式如下所示：
[0158][0159][0160]
其中，w
y
表示未接可再生dg时系统网络损耗，单位为kw；w
x
表示可再生dg接入后系统网络损耗，单位为kw；f
j
表示可再生dg所接入的节点负荷，单位为kva；f
z
表示系统负荷总量，单位为kva。
[0161]
本实施例中，将容量为600kw、功率因数为0.9且视为pq节点的可再生dg分别接入到ieee33配电网各节点中，利用牛顿
‑
拉夫逊潮流计算的方法计算系统潮流，接入各个节点带来的网络损耗如图4所示。
[0162]
利用上述综合比率算法求解综合比率，k1与k2取常数0.5，节点30、32、31、14、33、29、13、15节点的综合比率较大，因此取这些节点为较优接入位置。
[0163]
选取一个节点，如15号节点，将功率因数为0.9且视为pq节点的不同容量可再生dg接入，所得结果如表1所示。
[0164]
表1：不同容量可再生dg接入15节点潮流参数比较
[0165][0166]
同理，对其他节点也接入不同容量的可再生dg，得到各个节点的较优接入容量。
[0167]
s4：建立主动配电网上层模型和下层模型，上层模型的目标函数为年综合成本最小，下层模型的目标函数为网络损耗最小，结合负荷的数学模型、可再生dg的输出功率数学模型和可再生dg的最佳装机容量，构建上层模型和下层模型的约束条件。步骤s4包括以下步骤：
[0168]
s41：采用蒙特卡洛模拟法根据可再生dg的输出功率数学模型和主动配电网的负荷数学模型生成可再生dg的输出功率数据和主动配电网的负荷数据；
[0169]
以风机发电为例，应用蒙特卡洛对风机出力不确定性进行模拟，通过homer软件取得某地区的风速数据，根据双参数威布尔分布相关公式求解出weibull分布参数k和c分别为2.28和8.87，利用风速概率密度函数，通过10000次蒙特卡洛仿真模拟，得到蒙特卡洛风速概率密度模拟图如图5所示。
[0170]
结合风力发电概率密度函数，对风机发电出力不确定性进行描述，模拟得到其输出功率。同理，利用光伏概率密度函数及输出功率曲线，模拟光伏的有功出力。
[0171]
s42：以经济性为目标建立主动配电网建设及运行的上层模型；
[0172]
在上层可再生dg规划模型中，为了最大程度减少可再生dg接入主动配电网所产生的费用，充分考虑可再生dg接入主动配电网的经济性，目标函数为年综合成本最小，如下式所示：
[0173]
c
min
＝c
dg
+c
w
+c
g
+c
l
‑
c
h
[0174]
其中，c
min
表示年综合成本，c
dg
表示系统可再生dg的投资年费用，c
w
表示可再生dg运行维护费用，c
g
表示向上级电网购电成本，c
l
表示网络损耗费用，c
h
表示政府补贴费用。
[0175]
系统可再生分布式电源的投资年费用c
dg
的计算公式如下：
[0176][0177]
其中，r表示贴现率；k表示可再生dg使用年限，单位为年；n表示主动配电网节点数；c
pv,i
表示第i个节点的单位光伏投资费用，单位为万元/mw；p
pv,i
表示第i个节点的光伏安装量，单位为mw；c
wt,i
表示第i个节点的单位风机投资费用，单位为万元/mw；p
wt,i
表示第i个节点的风机安装量，单位为mw。本实施例中，可再生dg的经济使用年限为20年，折现率为10％。
[0178]
可再生分布式电源运行维护费用c
w
的计算公式如下：
[0179][0180]
其中，c
p,i
表示第i个节点的光伏运行维护费用，单位为万元/mwh；e
pv,i
表示第i个节点的光伏发电量，单位为mwh；c
w,i
表示第i个节点的风机运行维护费用，单位为万元/mwh；
e
wt,i
表示第i个节点的风机发电量，单位为mwh。
[0181]
向上级电网购电成本c
g
的计算公式如下：
[0182][0183]
其中，c
e
表示上级电网购电电价，单位为万元/mwh；p
y
表示主动配电网总有功负荷量，单位为mw；t
max
表示年最大负荷利用小时数，单位为h。本实施例中，电价为0.5元/kwh。
[0184]
网络损耗费用c
l
的计算公式如下：
[0185]
c
l
＝c
k
p
loss
τ
max
[0186]
其中，c
k
表示网络损耗电价，单位为万元/mwh；c
k
表示网络损耗，单位为mw；τ
max
表示年最大负荷损耗小时数，单位为h。
[0187]
政府补贴费用c
h
的计算公式如下：
[0188][0189]
其中，c
b
表示环境补贴费用，万元/mwh。
[0190]
为了保证系统稳定运行且考虑系统对可再生分布式电源的接纳能力，上层规划模型的约束条件包括系统节点容量约束、可再生dg安装总容量约束和可再生dg渗透率约束。
[0191]
系统节点容量约束，如下式所示：
[0192][0193]
其中，p
dg,i
表示第i个节点可再生dg安装量，单位为mw；k
i
表示主动配电网节点负荷值，单位为mw。
[0194]
本实施例中，可再生dg总装机容量不超过网络总负荷的25％，且考虑步骤s2中得到的最佳光伏装机容量和最佳风机装机容量。
[0195]
可再生dg安装总容量约束，如下式所示：
[0196][0197]
其中，p
dg,max
表示主动配电网最大准入容量。
[0198]
可再生dg渗透率约束，如下式所示：
[0199][0200]
s43：以网络损耗最小为目标建立主动配电网建设及运行的下层模型；
[0201]
在下层可再生dg规划模型中，以网络损耗最小为目标函数，其表达式如下式所示：
[0202]
f＝min(p
loss
)
[0203]
为了满足系统可靠运行，并考虑主动管理模式，下层规划模型的约束条件包括节点功率平衡约束、节点电压约束、支路传输功率约束、无功补偿投入装置容量约束、可再生dg出力约束
[0204]
节点功率平衡约束，如下式所示：
[0205][0206]
其中，p
is
表示节点i的有功注入，单位为mw；q
is
表示节点i的无功注入，单位为mvar；u
i
、u
j
表示节点i、j电压向量幅值，单位为mv；g
ij
和b
ij
表示电导和电纳；θ
ij
表示节点i和节点j电压相角差，单位为
°
。
[0207]
节点电压约束，如下式所示：
[0208][0209]
其中，u
i
表示节点i的电压幅值，单位为mv；n
j
表示系统节点数；u
i,min
和u
i,max
表示节点i电压幅值的下限和上限，单位为mv。本实施例中，设定节点电压范围为0.95～1.05p.u.。
[0210]
支路传输功率约束，如下式所示：
[0211][0212]
其中，s
n
表示支路s
n,max
的传输功率，单位为mw；n
l
表示系统支路数；s
n,max
表示支路s
n,max
传输功率上限，单位为mw。本实施例中，支路传输功率上限为4.5mw。
[0213]
无功补偿投入装置容量约束，如下式所示：
[0214][0215]
其中，q
ci
表示第i个可再生dg安装节点无功补偿设备投入量，单位为mvar；q
ci,min
和q
ci,max
表示无功补偿设备投入量下限和上限，单位为mvar。
[0216]
可再生dg出力约束，如下式所示：
[0217][0218]
其中，p
i
表示第i个节点可再生dg有功出力量，单位为mw；p
i,min
和p
i,max
表示可再生dg有功出力量的下限和上限，单位为mw。
[0219]
有载调压变压器分接头调节约束，如下式所示：
[0220]
k
p,min
≤k
p
≤k
p,max
[0221]
其中，k
p
表示有载调压变压器抽头位置；k
p,min
和k
p,max
表示有载调压变压器抽头调节范围的下限与上限。
[0222]
建立双层模型后，使用改进pso算法和跟踪中心内点轨迹法求解，整体求解模式如图6所示。
[0223]
步骤s5包括以下步骤：
[0224]
s51：根据可再生dg的接入位置生成编码向量作为个体，个体的粒子速度表示可再生dg的容量信息，生成初始种群；
[0225]
生成初始种群包括以下步骤：
[0226]
a1：随机生成编码向量和粒子速度；
[0227]
a2：判断该编码向量和粒子速度对应的可再生dg接入位置和容量信息是否满足上层模型的约束条件，若满足，则将该编码向量和粒子速度作为初始种群中的个体，否则，舍弃该编码向量和粒子速度；具体的，本实施例采用二进制编码。
[0228]
a3：重复步骤a1，直至初始种群中的个体数量满足预设置的种群规模。本实施例中，种群规模为100。
[0229]
s52：设置pso算法的参数值和遗传算法的参数值，对pso算法的惯性权重ω和学习因子c1、c2进行自适应改进，对遗传算法的交叉算子p
cro
和变异算子p
mut
进行自适应改进：
[0230][0231]
其中，ω
max
与ω
min
为最大惯性权重和最小惯性权重，maxiter为最大迭代次数，iter为当前迭代次数；本实施例中，ω
max
＝0.9，ω
min
＝0.5，最大迭代次数＝300。
[0232][0233]
其中，c
max
与c
min
为学习因子的最大和最小值；本实施例中，c
max
＝2，c
min
＝1。
[0234]
增大ω值能够提高pso全局搜索性能，减小ω值能够使得pso的局部搜索性能得到提高。合理确定惯性权重ω的值，能避免陷入局部最优，提高搜索性能。学习因子对收敛性以及算法寻优也有一定的影响，因此，对pso算法的惯性权重ω和学习因子c1、c2进行了自适应改进。
[0235]
为验证改进pso算法的性能，构建一个目标函数，使用传统pso算法对其进行求解；再使用改进pso算法，将惯性权重ω和学习因子c1、c2进行自适应改进，对目标函数进行求解，传统pso算法和改进pso算法求解目标函数迭代曲线如图7所示，传统pso算法在第238代收敛，而改进传统pso算法在63代收敛，不仅具有较快的收敛性，且目标函数值也比传统粒子群算法求得的值更优，由此说明改进pso算法的优越性。
[0236][0237]
其中，p
cro
、p
cro1
、p
cro2
表示交叉因子，g1表示常数，f(k+1)表示待交叉个体的适应度值，f(avg)表示群体所有个体平均适应度值，f(max)表示群体适应度最大个体的适应度值；
[0238]
本实施例中，p
cro1
＝0.9，p
cro2
＝0.6。
[0239]
[0240]
其中，p
mut
、p
mut1
、p
mut2
表示交叉因子，g2表示常数；
[0241]
本实施例中，p
mut1
＝0.1，p
mut2
＝0.06。
[0242]
s53：计算种群中每个个体的适应度，适应度是根据上层模型的目标函数值计算得到的，记录该代种群中的最优个体，若种群收敛，则执行步骤s55，否则，执行步骤s54；
[0243]
种群收敛具体为：最优个体的适应度满足预设置的收敛精度、或在kp(kp>1)次迭代中种群的最优个体的适应度的变化率小于预设置稳定阈值、或种群代数等于预设置的最大迭代次数，kp为预设置的稳定次数。
[0244]
s54：对种群进行选择操作、交叉操作和变异操作，更新pso算法的惯性权重和学习因子，计算每个个体的粒子速度，并更新每个个体的粒子位置，得到新一代种群，重复步骤s53。
[0245]
选择是对种群的优质个体进行保留并舍弃劣质个体，使优质个体传递给下一代种群。本实施例中，选用轮盘赌的方法，个体适应度对整体适应度总和的占比为被选概率，个体适应度值越大，则被选择概率越高，反之亦然。
[0246]
本实施例中，采用双切点交叉操作双切点交叉指在个体中随机设置一个交叉点，利用自适应公式交换交叉点前后两个体的部分染色体，形成两个新个体。对交叉算子进行了自适应改进，由于交叉算子与全局搜索能力相关，对其进行自适应改进能够提高遗传算法的性能。
[0247]
变异操作是先根据变异概率对群体中所有个体进行判断，如需要进行变异操作则随机选择个体进行变异。对变异算子进行了自适应改进，变异算子与局部搜索能力有关，变异概率对遗传算法全局最优解的求取也有重要作用，对其进行自适应改进可以加快对最优解求取的速度，防止出现“早熟现象”。
[0248]
s54：输出最优个体所对应的可再生dg接入位置和容量信息，并将其传递到下层模型。
[0249]
s6：通过跟踪中心轨迹内点法求解下层模型，得到可再生dg接入主动配电网的网络损耗并将结果返回给上层模型；
[0250]
s7：计算上层模型是否满足预设置的优化标准，如年综合成本是否满足预设置的成本指标，若满足，则输出可再生dg接入位置和容量信息，否则，重复步骤s5。
[0251]
以上详细描述了本发明的较佳具体实施例。应当理解，本领域的普通技术人员无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此，凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案，皆应在由权利要求书所确定的保护范围内。