基于改进多目标粒子群的电动汽车充电优化方法和系统与流程

[0001]
本发明属于电力企业管理领域，尤其是涉及一种基于改进多目标粒子群的电动汽车充电优化方法和系统。


背景技术：

[0002]
随着供给侧结构性和电力体制的变化的逐步深入，加强电力需求侧管理已成为完善电力供需平衡机制、优化能源消费结构、促进可再生能源消纳和提升智能用电水平的重要支撑。
[0003]
然而，大规模电动汽车在接入电网进行无序充电时也极易对电力系统的安全稳定运行造成影响。相关研究已经表明，电动汽车在没有相关组织机构进行规范的情况下，其充放电行为容易增大电网的负荷峰谷差和损耗，同时也会缩短变压器等电力设备的使用周期，从而不利于电网的稳定运行。因此，研究电动汽车最优的充电策略，正确引导电动汽车有序充放电，既可以大幅度地降低电动汽车用户充电费用，同时又能够有效地平滑电网负荷曲线，对电网安全稳定、经济运行具有重要的现实意义。
[0004]
多目标粒子群优化算最初是在对动物集群活动行为观察基础上，利用群体中的个体对信息的共享使整个群体的运动在问题求解空间中产生从无序到有序的演化过程，从而获得最优解。对于现有的多目标粒子算法来说虽然它搜索速度快、效率高、算法简单，但是对于离散的优化问题处理不佳，容易陷入局部最优等一些问题导致他寻优速度和准备性方面受到一定的影响。


技术实现要素：

[0005]
为克服上述现有技术的不足，本发明提出一种基于改进多目标粒子群的电动汽车充电优化方法，包括：
[0006]
获取电动汽车充电信息和电网负荷信息的数据输入到预先建立的电动汽车充电策略模型；
[0007]
采用多目标粒子群算法和小生境算法对电动汽车充电策略模型求解，得到各个时段电动汽车充电最优策略；
[0008]
基于各个时段电动汽车充电最优策略得到全时段电动汽车的充电最优充电策略；
[0009]
其中，所述电动汽车充电模型是基于充电费用最小和电网负荷方差最小建立的,所述电动汽车充电模型包括多个分别与电动汽车在各时刻内的充电功率对应的粒子。
[0010]
优选的，所述采用多目标粒子群算法和小生境算法对电动汽车充电策略模型求解，得到各个时段电动汽车充电最优策略，包括：
[0011]
建立各电动汽车在各时刻内的充电功率与改进多目标粒子群算法模型预设时段的粒子位置变量间的对应关系；
[0012]
根据多目标粒子群算法和小生境算法对电动汽车充电策略模型中的充电功率进行粒子位置更新；
[0013]
根据所述粒子位置更新的结果确定更新全局的最优位置；
[0014]
利用小生境算法对所述全局的最优位置的粒子进行精英解集的更新；
[0015]
判断更新的精英解集是否能够生成pareto最优解集，若是则在pareto最优解集中采用topsis的方法确定电动汽车充电策略模型中指定时段的充电功率的最优解，若不是则重新更新个体位置，直到能够生成pareto最优解集。
[0016]
优选的，所述根据多目标粒子群算法和小生境算法对电动汽车充电策略模型中的充电功率进行粒子位置更新，包括：
[0017]
在最优位置的集合里根据小生境算法确定电动汽车充电策略模型中孤立性最好的粒子；
[0018]
对所述电动汽车充电策略模型中孤立性最好的粒子和最优位置集合中存在的pareto改进的粒子进行变异操作，得到新的粒子；
[0019]
根据pareto优劣性，将所述新的粒子与其之前所经过的最优位置的粒子进行选择，得到该新的粒子的最新位置。
[0020]
优选的，所述在最优位置的集合里根据小生境算法确定电动汽车充电策略模型中孤立性最好的粒子，包括：
[0021]
根据当前迭代过程中电动汽车充电策略模型中每个粒子分别经过的最优位置的集合计算两个粒子之间的当量距离；
[0022]
根据两个粒子之间的当量距离计算粒子之间的共享值；
[0023]
依据粒子之间的当量距离对粒子的小生境半径进行动态更新并根据所述粒子之间的共享值计算粒子的小生境数；
[0024]
根据小生境数中最小的粒子确定独立性最好的粒子。
[0025]
优选的，所述根据所述粒子位置更新的结果确定更新全局的最优位置，包括：
[0026]
计算所有粒子的pareto支配数，选取支配数值最大所对应粒子的位置作为预设次数迭代的全局最优位置；
[0027]
当多个粒子的支配数值同时为最大值，任意选取其中一个粒子的位置作为预设次数迭代的全局最优位置。
[0028]
优选的，所述利用小生境算法对所述全局的最优位置的粒子进行精英解集的更新，包括：
[0029]
将每次寻优迭代中的全局最优粒子保存在种群外部的预设规模的精英解集中，并依据小生境数对该解集进行更新；
[0030]
如果精英解集没有饱和，则将全局最优粒子直接放入至精英解集；
[0031]
如果精英解集已经饱和，则根据按照小生境数去除数值较大的粒子。
[0032]
优选的，所述电动汽车充电策略模型的建立，包括：
[0033]
根据电动汽车充电信息和电网负荷信息以充电费用最小和电网负荷方差最小作为电动汽车充电策略的目标函数；
[0034]
以充电功率上限约束、充电功率波动约束和充电量约束作为电动汽车充电优化的约束条件，建立电动汽车充电策略模型。
[0035]
优选的，所述根据电动汽车充电信息和电网负荷信息以充电费用最小和电网负荷方差最小作为电动汽车充电策略的目标函数，包括：
[0036]
所述的充电费用最小的计算公式如下：
[0037][0038]
所述的电网负荷方差最小的计算公式如下：
[0039][0040]
式中minf
1
表示电动汽车的充电费用最小；minf
2
表示电网负荷方差最小；c
b
表示第b时段内的充电单价；(p
a,t
δtp)表示第a辆电动汽车在第b时段内的充电电量；p
a,b
表示第a辆电动汽车在第b时段内的充电功率；δtp表示每个时段的时间间隔；b表示时段总数；n表示电动汽车总数；p
load,b
表示第b时段内区域的普通负荷；表示第a辆电动汽车在第b时段内n个电动汽车的总负荷。
[0041]
优选的，所述建立各电动汽车在各时刻内的充电功率与改进多目标粒子群算法模型预设时段的粒子位置变量间的对应关系之后，且根据多目标粒子群算法和小生境算法对电动汽车充电策略模型中的充电功率进行粒子位置更新之前，还包括：
[0042]
确定粒子的初始位置和速度及其变化范围；
[0043]
利用粒子初始化的位置对个体最优位置和全局最优粒子位置进行初始化设置；
[0044]
设定精英解集为空集，并设置迭代次数的初值；
[0045]
对所有粒子的目标函数值进行计算。
[0046]
基于同一发明构思，本申请还提供了一种基于改进多目标粒子群的电动汽车充电优化系统，包括：获取数据模块、各时段充电最优模块和全时段充电最优模块；
[0047]
其中获取数据模块，用于获取电动汽车充电信息和电网负荷信息的数据输入到预先建立的电动汽车充电策略模型；
[0048]
其中各时段充电最优模块，用于采用多目标粒子群算法和小生境算法对电动汽车充电策略模型求解，得到各个时段电动汽车充电最优策略；
[0049]
其中全时段充电最优模块，用于基于各个时段电动汽车充电最优策略得到全时段电动汽车的充电最优充电策略；
[0050]
其中，所述电动汽车充电模型是基于充电费用最小和电网负荷方差最小建立的,所述电动汽车充电模型包括多个分别与电动汽车在各时刻内的充电功率对应的粒子。
[0051]
优选的，所述各时段充电最优模块，包括：对应关系子模块、粒子位置更新子模块、更新全局位置子模块、精英解集子模块和判断生成子模块；
[0052]
其中对应关系子模块，用于建立各电动汽车在各时刻内的充电功率与改进多目标粒子群算法模型预设时段的粒子位置变量间的对应关系；
[0053]
其中粒子位置更新子模块，用于根据多目标粒子群算法和小生境算法对电动汽车充电策略模型中的充电功率进行粒子位置更新；
[0054]
其中更新全局位置子模块，用于根据所述粒子位置更新的结果确定更新全局的最优位置；
[0055]
其中精英解集子模块，用于利用小生境算法对所述全局的最优位置的粒子进行精
英解集的更新；
[0056]
其中判断生成子模块，用于判断更新的精英解集是否能够生成pareto最优解集，若是则在pareto最优解集中采用topsis的方法确定电动汽车充电策略模型中指定时段的充电功率的最优解，若不是则重新更新个体位置，直到能够生成pareto最优解集。
[0057]
优选的，所述粒子位置更新子模块，包括：确定孤立性粒子单元、得到新粒子单元和最新位置单元；
[0058]
其中确定孤立性粒子单元，用于在最优位置的集合里根据小生境算法确定电动汽车充电策略模型中孤立性最好的粒子；
[0059]
其中得到新粒子单元，用于对所述电动汽车充电策略模型中孤立性最好的粒子和最优位置集合中存在的pareto改进的粒子进行变异操作，得到新的粒子；
[0060]
其中最新位置单元，用于根据pareto优劣性，将所述新的粒子与其之前所经过的最优位置的粒子进行选择，得到该新的粒子的最新位置。
[0061]
优选的，所述确定孤立性粒子单元，包括：确定距离子单元、共享值子单元、计算小生境数子单元和独立性最好子单元；
[0062]
其中确定距离子单元，用于根据当前迭代过程中电动汽车充电策略模型中每个粒子分别经过的最优位置的集合计算两个粒子之间的当量距离；
[0063]
其中共享值子单元，用于根据两个粒子之间的当量距离计算粒子之间的共享值；
[0064]
其中计算小生境数子单元，用于依据粒子之间的当量距离对粒子的小生境半径进行动态更新并根据所述粒子之间的共享值计算粒子的小生境数；
[0065]
其中独立性最好子单元，用于根据小生境数中最小的粒子确定独立性最好的粒子。
[0066]
优选的，所述更新全局位置子模块，包括：选择位置单元和选取一个位置单元；
[0067]
其中选择位置单元，用于计算所有粒子的pareto支配数，选取支配数值最大所对应粒子的位置作为预设次数迭代的全局最优位置；
[0068]
其中选取一个位置单元，用于当多个粒子的支配数值同时为最大值，任意选取其中一个粒子的位置作为预设次数迭代的全局最优位置。
[0069]
与最接近的现有技术相比，本发明具有的有益效果如下：
[0070]
1、本发明实现了一种基于改进多目标粒子群的电动汽车充电优化方法和系统，所述方法包括：获取电动汽车充电信息和电网负荷信息的数据输入到预先建立的电动汽车充电策略模型；采用多目标粒子群算法和小生境算法对电动汽车充电策略模型求解，得到各个时段电动汽车充电最优策略；基于各个时段电动汽车充电最优策略得到全时段电动汽车的充电最优充电策略；其中，所述电动汽车充电模型是基于充电费用最小和电网负荷方差最小建立的,所述电动汽车充电模型包括多个分别与电动汽车在各时刻内的充电功率对应的粒子，优化后可以大幅度降低充电费用和平滑电网负荷曲线，当出现大规模电动汽车接入电网时可以进行有序充电，提高了电动汽车充电的稳定性和准确性；
[0071]
2、本发明同时利用多目标粒子群算法和小生境算法提高了在粒子寻优过程中寻优的速度和准确性。
附图说明
[0072]
图1为本发明提供的基于改进多目标粒子群的电动汽车充电优化方法流程示意图；
[0073]
图2为本发明提供的基于改进多目标粒子群算法的电动汽车充电优化流程示意图；
[0074]
图3为本发明实施例中的商务办公区域日常电网负荷；
[0075]
图4为本发明实施例中的pareto前沿；
[0076]
图5为本发明实施例中的电网负荷比较；
[0077]
图6为本发明提供的基于改进多目标粒子群的电动汽车充电优化系统结构示意图。
具体实施方式
[0078]
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
[0079]
实施例1：
[0080]
下面结合附图及具体实施例对本发明的应用原理作进一步描述。
[0081]
如图1所示,本发明实施例的基于改进多目标粒子群的电动汽车充电优化方法包括：
[0082]
步骤1:获取电动汽车充电信息和电网负荷信息的数据输入到预先建立的电动汽车充电策略模型；
[0083]
步骤2:采用多目标粒子群算法和小生境算法对电动汽车充电策略模型求解，得到各个时段电动汽车充电最优策略；
[0084]
步骤3:基于各个时段电动汽车充电最优策略得到全时段电动汽车的充电最优充电策略；
[0085]
其中，所述电动汽车充电模型是基于充电费用最小和电网负荷方差最小建立的,所述电动汽车充电模型包括多个分别与电动汽车在各时刻内的充电功率对应的粒子。
[0086]
为了实现对电动汽车充电策略的优化，提出了基于改进多目标粒子群算法的电动汽车充电策略优化方法，选取充电费用最小和电网负荷方差最小作为目标函数，运用多目标粒子群算法进行寻优，同时加入小生境算法以提升寻优速度和准确性。并通过商务办公区的电动汽车充电策略优化进行仿真分析，验证了优化方法的效果，优化后的充电策略够较大幅度地降低充电费用，同时能够有效地平滑电网负荷曲线。
[0087]
步骤1:获取电动汽车充电信息和电网负荷信息的数据输入到预先建立的电动汽车充电策略模型；其中，所述电动汽车充电模型是基于充电费用最小和电网负荷方差最小建立的,所述电动汽车充电模型包括多个分别与电动汽车在各时刻内的充电功率对应的粒子。具体如下:
[0088]
1.目标函数的建立
[0089]
综合考虑用户侧的利益和电网侧的安全性，以充电费用最少、电网负荷方差最小作为优化目标，如下所示：
[0090]
(1)充电费用最少
[0091][0092]
式中c
b
——表示第b时段内的充电单价；
[0093]
(p
a,t
δtp)——表示第a辆电动汽车在第b时段内的充电电量；
[0094]
p
a,b
——表示第a辆电动汽车在第b时段内的充电功率；
[0095]
δtp——表示每个时段的时间间隔；
[0096]
b——表示时段总数；
[0097]
n——表示电动汽车总数。
[0098]
(2)电网负荷方差最小
[0099][0100]
式中p
load,b
——表示第b时段内区域的普通负荷；
[0101]
b——表示时段总数；
[0102]
n——表示电动汽车总数；
[0103]
——表示第a辆电动汽车在第b时段内n个电动汽车的总负荷。
[0104]
2.约束条件的建立
[0105]
电动汽车充电优化的约束条件如下所示：
[0106]
(1)充电功率上限约束
[0107]
电动汽车在充电时充电功率不能超过充电站的充电功率上限，如式(3)所示。
[0108]
0≤p
a,b
≤p
max
ꢀꢀ
(3)
[0109]
式中p
a,b
——表示第a辆电动汽车在第b时段内的充电功率；
[0110]
p
max
——表示电动汽车充电功率的上限。
[0111]
(2)充电功率波动约束
[0112]
对电动汽车充电时的充电功率波动范围进行限制，用以防止过大的充电功率波动对电网负荷造成较大的影响，如式(4)所示。
[0113][0114]
式中p
a,b
——表示第a辆电动汽车在第b时段内的充电功率；
[0115]
δtp——表示每个时段的时间间隔；
[0116]
n——表示电动汽车总数。
[0117]
p
a,b+1
——表示第a辆电动汽车在第b+1时段内的充电功率；
[0118]
p
f
——表示充电过程中两个相邻时段的充电功率波动最大值。
[0119]
(3)充电量约束
[0120]
电动汽车充电结束时的电量需小于等于电动汽车电池的最大容量，如式(5)所示。
[0121]
[0122]
式中——表示第a辆电动汽车预期的充电电量；
[0123]
——表示第a辆电动汽车充电前的已有电量；
[0124]
——表示第a辆电动汽车电池最大容量；
[0125]
p
a,b
——表示第a辆电动汽车在第b时段内的充电功率；
[0126]
δtp——表示每个时段的时间间隔；
[0127]
b——表示时段总数；。
[0128]
步骤2:采用多目标粒子群算法和小生境算法对电动汽车充电策略模型求解，得到各个时段电动汽车充电最优策略；具体如下：
[0129]
3.基于改进多目标粒子群的电动汽车充电优化
[0130]
(1)多目标优化
[0131]
多目标优化问题表示的是在实际问题中绝大多数都有2个及以上目标需要能够同时满足，即对这些目标进行同时优化处理。但是这些目标往往是相互冲突的，当对其中一个目标的性能进行提升时通常会导致其他目标的性能下降。为了解决这一问题，通常对各目标赋予一定的权重系数，将其转化为单目标优化问题。但是，采用这种方法得到的解很有可能是局部最优解，即仍存在一个或多个别的解能同时提升多目标的性能，此外，由于目标之间的量纲存在差异，所以进行组合时很难保证与其实际意义保持一致。鉴于此，本文采用pareto最优解集的相关理论方法来求解多目标优化问题。
[0132]
定义1：pareto支配：称一个向量u＝(u
1
,u
2
,
…
,u
m
)支配(或非劣于)向量v＝(v
1
,v
2
,
…
,v
m
)，当且仅当对于使u
i
＜v
i
，记为u＜v。
[0133]
定义2：pareto最优：如果x
*
∈s，并且在s中不存在比x
*
更优的解x，那么称x
*
是可行解集合s的pareto最优解。
[0134]
定义3：pareto最优解集：由满足定义pareto最优解条件的所有x
*
构成的集合，记为：
[0135][0136]
定义4：pareto前沿：所有pareto最优解对应的目标函数值所构成的区域称为pareto前沿，记为：
[0137]
p
f
＝{f(x)＝(f
1
(x),f
1
(x),
…
,f
m
(x))|x∈p
*
}
ꢀꢀ
(7)
[0138]
(2)粒子群优化
[0139]
粒子群优化算法(particle swarm optimization，pso)是一种模拟鸟群觅食的智能优化算法。种群中的每个个体被看成一个可能的最优解，设定种群在n维空间中进行搜索，其规模为n，第i个粒子的速度和位置记为v
i
＝(v
i1
,v
i2
,
…
,v
in
)和x
i
＝(x
i1
,x
i2
,
…
,x
in
)，记录该粒子在搜索过程中经过的最优位置p
i
＝(p
i1
,p
i2
,
…
,p
in
)以及整个群体在搜索过程中经过的最优位置p
g
。种群个体将自身所经过的最优位置与群体最优位置进行比较来不断更新自身的速度和位置，使其离群体最优位置的距离缩短，从而实现种群的进化。第i个粒子的第z维速度和位置更新方式见式(8)和式(9)。
[0140]
[0141][0142]
式中t——表示种群当前的迭代次数；
[0143]
c
1
，c
2
——表示加速常数；
[0144]
r
1
，r
2
——表示属于[0,1]的随机数；
[0145]
——表示在种群当前的第t次迭代次数第i个粒子的第z维速度；
[0146]
——表示在种群当前的第t次迭代次数第i个粒子的第z维速度在搜索过程中经过的最优位置；
[0147]
——表示在种群当前的第t次迭代次数第i个粒子的第z维速度的位置；
[0148]
——表示在种群当前的第t次迭代次数第g个粒子的第z维速度在搜索过程中经过的最优位置；
[0149]
——表示在种群当前的第i+1次迭代次数第i个粒子的第z维速度；
[0150]
w——表示惯性权重系数。
[0151]
(3)改进的多目标粒子群优化
[0152]
标准的粒子群优化算法由于其加速参数通常为固定值，所以容易限制粒子个体的更新，不利于其寻找到全局最优解。鉴于此，本文将小生境引入到粒子群优化算法中以提高其寻优的效率，主要操作步骤如下所示：
[0153]
1)更新粒子所经过的最优位置
[0154]
更新粒子所经过的最优位置包括确定孤立粒子、对存在pareto改进的粒子进行变异操作和选择操作。
[0155]
step1：根据小生境数确定孤立性最好的粒子。令当前迭代过程中所有个体所经过的最优位置的集合为p＝{p
i
＝(p
i1
,p
i2
,
…
,p
in
)|i＝1,2,
…
,n}，对于p中任意两个粒子i、j之间的当量距离可以用式(10)来表示。
[0156][0157]
式中d
ij
——表示粒子i和粒子j之间的当量距离；
[0158]
——表示粒子i的第s个目标；
[0159]
——表示粒子j第s个目标；
[0160]
m
——表示粒子数量；
[0161]
——表示第s个目标的最大值；
[0162]
——表示第s个目标的最小值。
[0163]
采用式(11)来计算粒子之间的共享值。
[0164][0165]
式中d
ij
——表示粒子i和粒子j之间的当量距离；
[0166]
σ
share
——表示粒子的小生境半径。
[0167]
依据粒子之间的当量距离对粒子的小生境半径进行动态更新，如式(12)所示。
[0168][0169]
则粒子i的小生境数记为：
[0170][0171]
那么小生境数最小的粒子为独立性最好的粒子，记为：q＝(q
1
,q
2
,
…
,q
n
)。
[0172]
step2：将p中相比较粒子q存在pareto改进的粒子与q进行变异操作，得到新的粒子且
[0173]
step3：根据pareto优劣性，将新得到的粒子与其之前所经过的最优位置的粒子进行比较和选择，如式(14)所示。
[0174][0175]
将小生境数变异操作和选择操作引入到多目标粒子群优化算法中能够提高粒子的多样性，进而有利于粒子全局寻优效率。
[0176]
2)更新全局最优位置
[0177]
对全局最优位置的更新包括如下步骤：
[0178]
step1：在寻优迭代到第t次时，根据第1)步的操作，完成对粒子所经过的最优位置的更新后，计算所有粒子的pareto支配数rand(u)，选取支配数值最大所对应粒子的位置作为第t次迭代的全局最优位置。
[0179]
step2：如果有多个粒子的支配数值同时为最大值，由于这些粒子之间并不存在支配关系并且均是pareto最优解，所以可以任意选取其中一个粒子的位置作为第t次迭代的全局最优位置，如式(15)所示。
[0180][0181]
依据pareto支配数对迭代过程中的全局最优位置进行更新能够让第t次迭代中的具有全局最优位置的粒子加入到第t+1次的迭代中，促使群体朝向全局最优方向运动，提高了算法的收敛效率。
[0182]
3)更新精英解集
[0183]
将每次寻优迭代中的全局最优粒子保存在种群外部的一个规模为m的精英解集中，并依据小生境数对该解集进行更新。
[0184]
如果精英解集没有饱和，则可以将全局最优粒子直接放入至精英解集。
[0185]
如果精英解集已经饱和，则根据按照小生境数去除数值较大的粒子。
[0186]
利用小生境数对精英解集进行更新能够有效去除位于密集区域的粒子，并将分散区域的粒子保留下来，从而使pareto前沿的分布变得更加均匀。
[0187]
步骤3:基于各个时段电动汽车充电最优策略得到全时段电动汽车的充电最优充电策略；具体包括：根据上述步骤得到了各个时段的电动汽车充电优化策略，对每个时段的
策略进行汇总，得到全时段电动汽车的充电最优策略。
[0188]
4.基于改进多目标粒子群的电动汽车充电优化流程，如图2所示。
[0189]
运用前文提出了改进多目标粒子群算法对电动汽车的充电策略进行优化，主要包括如下步骤：
[0190]
(1)设定第i辆电动汽车在第t时段内的充电功率p
i,t
为待寻优的变量。
[0191]
(2)对粒子群进行初始化，主要包括粒子的初始位置和速度及其变化范围。利用初始化粒子的位置对个体最优位置和全局最优粒子位置进行初始化设置。设定精英解集为空集，迭代次数t＝0。
[0192]
(3)对所有粒子的目标函数值f
1
,f
2
,
…
,f
m
进行计算。
[0193]
(4)利用式(8)和式(9)对粒子的速度和位置进行更新。
[0194]
(5)利用式(10)至式(14)对孤立粒子进行确定，运用变异操作和选择操作对粒子所经过的最优位置集合进行更新。
[0195]
(6)利用式(15)对群体的全局最优位置进行更新。
[0196]
(7)判定约束条件是否满足，如果满足，则输出一组pareto最优解集，否则t＝t+1并重复步骤(2)至步骤(7)，继续寻优迭代。
[0197]
(8)采用topsis法在对pareto最优解集中对充电策略的最优解进行选取。
[0198]
实施例2
[0199]
本专利选取某商务办公区域的电动汽车充电站的用户为例，主要是该区域的各企业员工，即企业员工上班时将电动汽车放在充电站进行充电，下班时将电动汽车开走。所以商务办公区的电动汽车充电负荷主要集中的白天，这里对8:00-18:00时段的电动汽车充电策略进行仿真优化。
[0200]
1.参数设置
[0201]
商务办公区充电站的充电价格采用某市区一般工商业销售电价，具体划分时段如表1所示。
[0202]
表1分时电价设置
[0203][0204]
对约束条件及改进粒子群算法中的参数进行设置，见表2。
[0205]
表2商务办公区域场景的参数设置
[0206][0207][0208]
该商务办公区域日常的电网负荷特性如图3所示，可以看出，商务办公区域的日常负荷高峰出现在11:00和15:00，负荷低谷出现在13:00。
[0209]
该商务办公区域日常的电网负荷特性如图3所示，可以看出，商务办公区域的日常负荷高峰出现在11:00和15:00，负荷低谷出现在13:00。
[0210]
2.结果分析
[0211]
运用改进多目标粒子群算法对商务办公区域电动汽车充电策略进行优化，图4展示的是改进多目标粒子群算法在进行优化计算时所得到的pareto前沿。
[0212]
从图4可以看出运用改进多目标粒子群算法得到的pareto前沿分布比较均匀和稳定，说明所提算法能够有效求解对求解商务办公区域电动汽车充电策略多目标优化问题。随着总负荷标准差的增大，充电费用将会逐渐降低，这是由于安排在低价时段的电动汽车充电负荷越大，总的充电费用是越低的。
[0213]
为了从pareto最优解集中选出更合理的电动汽车优化充电策略，使用topsis获取最终的充电策略，结果如表3、表4和图5所示。
[0214]
表3商务办公区域的充电功率
[0215][0216]
表4商务办公区域充电优化策略
[0217][0218]
从表3和表4可以看出，经过优化的充电策略在高峰电价时段的充电功率维持在一个较低的水平，而在平段电价时段的充电功率较大，总体充电费用为1102.43元，相比较无序充电的费用1422.44元降低了22.50％。从图5可以看出，电动汽车无序充电会增大电网负荷方差，即增大电网负荷的波动，对电网负荷产生不利影响；采取优化充电策略后，电网负荷方差明显小于无序充电负荷以及日常负荷，从19080.74kw
2
下降到5055.96kw
2
，说明该充电策略对平滑电网负荷曲线的效果较为明显。综上，对商务办公区电动汽车充电策略进行优化能够较大幅度地降低充电费用，同时能够有效地平滑电网负荷曲线。
[0219]
实施例3：
[0220]
基于同一发明构思，本发明还提供了基于改进多目标粒子群的电动汽车充电优化系统，如图6所示，包括：获取数据模块、各时段充电最优模块和全时段充电最优模块；
[0221]
其中获取数据模块，用于获取电动汽车充电信息和电网负荷信息的数据输入到预先建立的电动汽车充电策略模型；
[0222]
其中各时段充电最优模块，用于采用多目标粒子群算法和小生境算法对电动汽车充电策略模型求解，得到各个时段电动汽车充电最优策略；
[0223]
其中全时段充电最优模块，用于基于各个时段电动汽车充电最优策略得到全时段电动汽车的充电最优充电策略；
[0224]
其中，所述电动汽车充电模型是基于充电费用最小和电网负荷方差最小建立的,所述电动汽车充电模型包括多个分别与电动汽车在各时刻内的充电功率对应的粒子。
[0225]
所述各时段充电最优模块，包括：对应关系子模块、粒子位置更新子模块、更新全局位置子模块、精英解集子模块和判断生成子模块；
[0226]
其中对应关系子模块，用于建立各电动汽车在各时刻内的充电功率与改进多目标粒子群算法模型预设时段的粒子位置变量间的对应关系；
[0227]
其中粒子位置更新子模块，用于根据多目标粒子群算法和小生境算法对电动汽车充电策略模型中的充电功率进行粒子位置更新；
[0228]
其中更新全局位置子模块，用于根据所述粒子位置更新的结果确定更新全局的最优位置；
[0229]
其中精英解集子模块，用于利用小生境算法对所述全局的最优位置的粒子进行精英解集的更新；
[0230]
其中判断生成子模块，用于判断更新的精英解集是否能够生成pareto最优解集，若是则在pareto最优解集中采用topsis的方法确定电动汽车充电策略模型中指定时段的充电功率的最优解，若不是则重新更新个体位置，直到能够生成pareto最优解集。
[0231]
所述粒子位置更新子模块，包括：确定孤立性粒子单元、得到新粒子单元和最新位置单元；
[0232]
其中确定孤立性粒子单元，用于在最优位置的集合里根据小生境算法确定电动汽车充电策略模型中孤立性最好的粒子；
[0233]
其中得到新粒子单元，用于对所述电动汽车充电策略模型中孤立性最好的粒子和最优位置集合中存在的pareto改进的粒子进行变异操作，得到新的粒子；
[0234]
其中最新位置单元，用于根据pareto优劣性，将所述新的粒子与其之前所经过的最优位置的粒子进行选择，得到该新的粒子的最新位置。
[0235]
所述确定孤立性粒子单元，包括：确定距离子单元、共享值子单元、计算小生境数子单元和独立性最好子单元；
[0236]
其中确定距离子单元，用于根据当前迭代过程中电动汽车充电策略模型中每个粒子分别经过的最优位置的集合计算两个粒子之间的当量距离；
[0237]
其中共享值子单元，用于根据两个粒子之间的当量距离计算粒子之间的共享值；
[0238]
其中计算小生境数子单元，用于依据粒子之间的当量距离对粒子的小生境半径进行动态更新并根据所述粒子之间的共享值计算粒子的小生境数；
[0239]
其中独立性最好子单元，用于根据小生境数中最小的粒子确定独立性最好的粒子。
[0240]
所述更新全局位置子模块，包括：选择位置单元和选取一个位置单元；
[0241]
其中选择位置单元，用于计算所有粒子的pareto支配数，选取支配数值最大所对应粒子的位置作为预设次数迭代的全局最优位置；
[0242]
其中选取一个位置单元，用于当多个粒子的支配数值同时为最大值，任意选取其中一个粒子的位置作为预设次数迭代的全局最优位置。
[0243]
所述精英解集子模块，包括：更新单元、饱和单元和不饱和单元；
[0244]
其中更新单元，用于将每次寻优迭代中的全局最优粒子保存在种群外部的预设规模的精英解集中，并依据小生境数对该解集进行更新；
[0245]
其中饱和单元，用于如果精英解集没有饱和，则将全局最优粒子直接放入至精英解集；
[0246]
其中不饱和单元，用于如果精英解集已经饱和，则根据按照小生境数去除数值较大的粒子。
[0247]
该系统还包括：模型建立模块；
[0248]
其中模型建立模块，包括目标函数子模块和约束子模块；
[0249]
其中目标子模块，用于根据电动汽车充电信息和电网负荷信息以充电费用最小和电网负荷方差最小作为电动汽车充电策略的目标函数；
[0250]
其中约束子模块，用于以充电功率上限约束、充电功率波动约束和充电量约束作为电动汽车充电优化的约束条件，建立电动汽车充电策略模型。
[0251]
目标函数中，充电费用最小的计算公式如下：
[0252][0253]
电网负荷方差最小的计算公式如下：
[0254][0255]
式中min f
1
表示电动汽车的充电费用最小；min f
2
表示电网负荷方差最小；c
b
表示
第b时段内的充电单价；(p
a,t
δtp)表示第a辆电动汽车在第b时段内的充电电量；p
a,b
表示第a辆电动汽车在第b时段内的充电功率；δtp表示每个时段的时间间隔；b表示时段总数；n表示电动汽车总数；p
load,b
表示第b时段内区域的普通负荷；表示第a辆电动汽车在第b时段内n个电动汽车的总负荷。
[0256]
该系统还包括：确定速度范围模块、初始化设置模块、设置迭代初值模块和目标计算模块；
[0257]
其中确定速度范围模块，用于确定粒子的初始位置和速度及其变化范围；
[0258]
其中初始化设置模块，用于利用粒子初始化的位置对个体最优位置和全局最优粒子位置进行初始化设置；
[0259]
其中设置迭代初值模块，用于设定精英解集为空集，并设置迭代次数的初值；
[0260]
其中目标计算模块，用于对所有粒子的目标函数值进行计算。本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
[0261]
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0262]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0263]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0264]
最后应当说明的是:以上实施例仅用于说明本申请的技术方案而非对其保护范围的限制,尽管参照上述实施例对本申请进行了详细的说明,所属领域的普通技术人员应当理解:本领域技术人员阅读本申请后依然可对申请的具体实施方式进行种种变更、修改或者等同替换，但这些变更、修改或者等同替换，均在申请待批的权利要求保护范围之内。