一种基于量子行为粒子群算法的微电网调度方法
【技术领域】
[0001] 本发明属于微网技术领域，尤其针对风光储并网型微电网优化调度，涉及一种基 于量子行为粒子群算法的微电网调度方法。
【背景技术】
[0002] 微电网是由多种DG(distribution generation)、储能单元、负荷和控制保护系统 组成的，通过协调各DG，为一个小区或者海岛提供电力或热负荷的需求，微电网的技术特点 使其适合为一些边远地区供电。并网型风光储微电网作为一种典型的微电网在未来具有极 大的发展空间，并网运行时，由于大电网实施了峰谷电价，这使得并网型微电网可以在低电 价时从大电网购电，而在高电价时售电给大电网从而赚取差价，降低整个系统运行成本。
[0003] 微电网的调度的目标可以是单目标，也可以是多目标，目前单目标主要针对系统 某一阶段运行成本而言，比如日前调度则是考虑全天运行成本，而多目标主要考虑系统运 行成本，环境惩罚成本以及可靠性等。微电网优化调度问题是一个多维非线性优化问题，对 于传统的通过求解方程得到最优解方法难以应用于解决微网优化调度问题中。而随着智能 控制和智能算法的发展，智能算法对于解决高维非线性最优问题已经比较成熟，比较典型 的智能优化算法主要有进化遗传算法、鱼群算法、蚁群算法以及粒子群算法。通过对微网优 化调度问题的建模并结合智能算法，可以有效的求解此类高维非线性优化问题。
[0004] 现有微电网调度技术主要针对是经济性而言，其中包括机组运行维护成本、燃料 成本、折旧成本、启停机成本以及在并网状态下运行和大电网的交互成本，对于环境成本， 主要是针对柴油机、微型燃气轮机等有污染物排放机组而言的。如何降低微电网运行维护 成本和与大电网交互成本成为并网型风光储微网的关键问题，由于各机组出力与负荷以及 各机组的出力上下限均有关系，在随机类智能搜索算法中经常会出现变量越限或约束不满 足要求的情况，因此如何处理好约束也成为优化调度中的一个关键。此外，随机类智能算法 的容易早熟以及不收敛问题也给多维非线性问题的求解带来了一定的困难，且算法的全局 搜索能力和局部搜索能力对非线性问题的寻优也有一定的影响，而如何权衡两者使得算法 适合实际问题的需要也是目前需要研究的一个问题。
【发明内容】
[0005] 本发明主要针对并网型风光储微电网优化调度问题，提供了一种基于量子行为粒 子群算法的微电网优化调度方法，且通过对算法的参数控制，使得其适用于风光储微电网 优化调度问题，且能较好的处理各约束问题。本发明主要通过以下技术方案实现：
[0006] 为解决风光储并网型微电网优化调度问题，本发明主要采用的技术方案是基于量 子行为粒子群算法的并网型风光储微电网优化调度方法，所述方法主要步骤有：
[0007] A、基于日前调度计划安排各微电源出力，选定储能系统中的蓄电池出力为直接优 化变量，而微电网和大电网交换功率为间接优化变量;设定仿真的调度周期和时间维度，载 入光伏出力、风机出力和负荷功率预测数据以及微电网和大电网交换功率的实时电价，同 时设定算法的粒子个数、迭代次数、运行轮数以及粒子维数；
[0008] B、初始化蓄电池的荷电状态并计算蓄电池初始剩余电量;对每个粒子的每一维依 次进行循环，若蓄电池荷电状态满足所要求的约束范围，则对其在出力范围内进行随机初 始化，计算该时刻蓄电池荷电状态，然后对蓄电池荷电状态进行检测，如果蓄电池荷电状态 不在所要求的约束范围内，则对其进行边界化条件处理;接着计算微电网和大电网的交换 功率，并检测其是否超出传输线最大传输功率限制，如超出线路最大传输功率限制，则把传 输功率置为最大允许传输功率，同时由功率平衡再次计算蓄电池出力，由蓄电池出力再计 算其荷电状态，检测蓄电池荷电状态是否满足要求，如满足要求，则跳出该粒子循环，进入 下一个粒子的初始化；
[0009] C、当所有粒子初始化完毕，则对其初始适应度值进行计算，也即目标函数值f(x)， f(x)主要包括以下部分；
[0011] R(i)为微电网和大电网交换功率的实时电价，Pgrid(i)为微电网与大电网的交换 功率，i为调度的某个时间维度，T为调度周期，Pn为储能蓄电池的额定输出功率，Tn为储能蓄 电池年运行小时数，Cstart为储能系统的初始总投资成本，P b(i)为储能蓄电池在某调度时刻 的实际输出功率，K。·为储能系统的资本回收系数，其由下式计算可得：
[0013] s为储能蓄电池的折旧率，M为储能蓄电池的使用寿命；
[0014] D、将每个粒子初始适应度值设定为该粒子个体最优适应度值，将粒子初始位置设 定为粒子个体最优位置，比较所有粒子适应度值并找出全局最优粒子，并记录其位置和适 应度值；
[0015] E、对每个粒子进行更新迭代，迭代完毕后对蓄电池每一时刻荷电状态进行计算， 同时检测荷电状态是否满足约束条件，如不满足则对其进行边界化处理，并再次修正蓄电 池出力。当蓄电池荷电状态超出下边界SOC min时，则用(3)式修正蓄电池出力；当蓄电池荷电 状态超出上边界SOCmax时，则用(4)式修正蓄电池出力。
[0018] 其中Xt表示蓄电池在t时刻的输出功率，RL为蓄电池初始剩余容量，NL为蓄电池额 定容量。
[0019] F、由功率平衡计算微网和大电网交换功率，同时检测其交换功率是否超出约束范 围，若超出约束范围则进行边界化处理，并重新计算蓄电池出力和荷电状态，最后判断蓄电 池荷电状态是否满足条件，如满足则跳出该粒子循环，进入下一个粒子的循环；
[0020] G、所有粒子循环完毕后重新计算适应度值，并更新个体最优位置和个体最优适应 度值以及全局最优位置和全局最优适应度值，重新返回E进行循环；
[0021] H、当达到最大迭代次数时，终止迭代，输出最优出力、适应度值以及对应的蓄电池 荷电状态。
[0022] 本发明具有的优点及有益效果为：
[0023] 1、目前还未有将量子行为粒子群应用到并网型风光储微网优化调度中，本方法为 首次采用量子行为粒子群算法应用到此问题。与传统粒子群算法相比，量子粒子群算法采 用的是量子力学中粒子的波粒二象性的原理而引入的，通过对粒子建立势阱模型，把传统 粒子群算法通过速度来改变粒子位置变为通过吸引子来改变粒子位置，使粒子迭代更简 单，但寻优能力更强。同时在本发明中将粒子迭代过程中的收缩扩张系数以及吸引子模型 中与粒子个体最优位置和全局最优位置的两个固定参数均设置为跟随迭代次数动态变化 的，使算法在开始的时候全局搜索能力较强，而迭代结束时局部搜索能力较强，这样能使算 法较快的收敛到全局最优解，解决传统算法早熟和不收敛等问题。
[0024] 2、针对并网型风光储微电网功率平衡约束问题目前应用比较多的是在目标函数 中加入惩罚项，通过对偏离功率平衡的粒子相应目标函数值中加入惩罚项，从而使函数值 增大，这样通过算法寻优将这部分粒子过滤掉，但这样会造成粒子数量的减少，影响算法的 搜索性能。本发明通过将蓄电池作为直接优化变量和微网与大电网交换功率作为间接优化 变量的方式降低了粒子维度，并且对不满足荷电状态约束的粒子和与大电网交换功率约束 的粒子从新初始化和循环迭代，从而保证了粒子的搜索性能，并能相应提高收敛速度。
【附图说明】
[0025] 图1是本发明提供的基于量子行为粒子群算法的并网型风光储微电网优化调度流 程图；
[0026] 图2是本发明的一个具体实例优化调度最终方案图。
【具体实施方式】
[0027]下面结合附图对本发明的【具体实施方式】做进一步的详细介绍：
[0028] 本发明提出的基于量子行为粒子群算法的并网型风光储微电网优化调度方法流 程图如图1所示，主要包括如下步骤：
[0029] A、设置算法的各基本参数，所需设置参数主要有粒子个数、粒子维数、迭代次数， 粒子规模数取H = 20,粒子维数也即需要优化的蓄电池出力的时刻数，此处作为日前调度， 调度周期为24h，每Ih为一个时段，故粒子维数取T = 24,迭代次数取L = 500。
[0030] B、载入光伏、风机的预测出力和负荷的预测功率以及实时电价等数据。
[0031] C、设置所有粒子的初始荷电状态均一致，根据蓄电池总容量计算蓄电池剩余容量 RL
[0032] D、初始化粒子每一维值，也即在蓄电池出力上下限之间随机赋值，并计算蓄电池 核荷电状态荷电状态和微电网与大电网的交换功率Pgrid，其值由下式功率平衡约束决定：
[0033] Prgid = Pioad-Ppv-Pwt-Pb (5)
[0034] 其中Plciad为负荷的需求功率，Ppv为光伏发电的出力，Pwt为风机的出力，Pb为蓄电池 出力。
[0035] E、根据蓄电池荷电状态和微电网与大电网交换功率Pgrid约束对蓄电池和Pgrid进行 调整，蓄电池荷电状态约束用下述表达式表示：
[0036] SOCmin < SOC < SOCmax (6)
[0037] 其中SOCmin和SOCmax分别是蓄电池允许的荷电状态下限和上限；
[0038] Pgrid约束用下述表达式表示：
[0039] Pgrid min ^ Pgrid ^ Pgrid max (7)
[0040 ]其中Pgr i d_min和Pgr i d_max分别是风光储微电网和大电网联络线传输功率的最小值和 最大值；
[0041]蓄电池荷电状态小于其允许的最小值时，说明蓄电池处于过放状态，则让蓄电池 放到允许的最低荷电状态，并计算蓄电池所能允许的最高放电功率，把这个最高放电功率 作为蓄电池的实际出力；当蓄电池荷电状态大于其允许的最大值时，说明此时蓄电池处于 过充状态，则让蓄电池充到