基于改进灰狼算法的并网型微电网容量优化研究.pdf

收稿日期：2022-05-06作者简介：田雨明（1997-），男，辽宁兴城人，硕士研究生。基于改进灰狼算法的并网型微电网容量优化研究田雨明a，尹常永b（沈阳工程学院 a.电力学院；b.新能源学院；辽宁 沈阳 110136）摘要：为了合理配置并网型微电网中各分布式电源的容量，降低微电网的经济成本，对并网型微电网的容量优化配置方法进行了研究。首先，建立了并网型微电网中各单元的数学模型，并提出了微电网的运行策略；其次，综合考虑系统的年均投资成本、年运行维护成本、置换成本、年停电补偿费、年购电费用、年运行收益和设备报废回收效益等因素建立目标函数，并结合微电网的功率平衡约束、设备的运行约束及与主网的交换功率约束来构建并网型微电网的容量优化配置模型；再次，采用改进的灰狼算法对优化模型进行求解；最后，通过仿真分析验证了所提出模型的有效性。关键词：微电网；并网型；容量配置；改进灰狼优化算法中图分类号：TM715文献标识码：A文章编号：1673-1603（2023）04-0056-08为实现“双碳”目标，以风能、太阳能等为主导的绿色能源被大规模应用于电力生产当中，有力推动了分布式发电形式的快速发展。与集中式发电相比，分布式发电形式具备安装位置灵活、建设周期短和污染小等诸多优点，同时还是解决能源短缺问题的重要手段1。但是，随着电网中分布式发电的占比不断扩大，其自身的波动性与不确定性会给电网安全运行造成严重的影响。目前，微电网技术可以很好地缓解这类问题。微电网内集成各类分布式电源与负荷，能够对电源进行合理有效的管理，已经得到了广泛的应用2-3。微电网容量优化作为建设微电网前的重要任务，合理地配置微网中各类分布式电源的容量不仅能够提高微电网运行的可靠性，还可以降低系统的投资运行成本。国内外已有众多学者对离网型和并网型微电网的容量优化配置方法展开研究。文献 4 同时考虑了经济性目标和可靠性目标建立离网型风/光/氢/储微电网容量优化配置模型，并采用多目标天牛须探索算法进行求解。文献 5 根据偏远地区的实时气象数据对独立型的风/光/储/柴智能微电网进行容量配置，以度电成本与负荷缺电率为优化目标，采用鲸鱼算法求得最优容量配置方案。相对于离网型微电网，在配置并网型微电网容量时还需要考虑微电网与大电网之间的功率交互，所以优化模型中应引入微电网的购售电成本和交换功率约束。文献 6针对并网型海岛微电网，以各设备安装成本、运行维护成本、污染治理成本、余电上网收益及购电成本来构建目标函数，采用遗传算法对容量进行配置。文献 7 从系统运行的净收益和水风互补度的角度来构建微电网容量优化配置模型，对引入需求侧响应下的微电网配置结果进行分析，证明了引入需求侧响应可以提高系统的经济性。但在以上文DOI：10.13888/ki.jsie（ns）.2023.04.012第 19 卷第 4 期2 0 2 3 年 10 月Vol.19 No.4Oct.2023沈阳工程学院学报（自然科学版）Journal of Shenyang Institute of Engineering（Natural Science）献中，并网型微电网容量优化配置模型均未考虑储能电池的置换成本，由于现阶段电池储能寿命周期年限较小，有必要在模型中加入储能电池的置换成本。针对上述并网型微电网容量优化配置模型存在的问题，本文对并网型微电网容量优化配置方法展开研究，以北方某地区作为研究对象，考虑储能的置换成本，建立以系统年综合经济成本为优化目标的并网型风光储微电网容量优化配置模型。为提升灰狼算法（Grey Wolf Optimization，GWO）的全局寻优性能，在 GWO 中引入了非线性收敛因子。最后通过MATLAB软件编程进行仿真，分析不同停电惩罚系数对容量配置结果的影响，验证微电网优化配置模型的有效性。1分布式电源模型1.1风机模型风力发电机输出功率主要受风速影响，具有不确定性，两者之间的关系可由如下关系式表示8：Pwt()v=0,v viP()v,vi v vNPN,vN v vo0,v vo（1）P()v=PNv-vivN-vi（2）式中，PN为风机额定输出功率；v为风机轮毂处实际风速；vi为切入风速；vN为额定风速；vo为切出风速；P（v）为实际风速的函数。当风速处于vi和vN之间时，风机的输出功率可由式（2）求出。由于实际风速与风机轮毂高度有关，需要对风速进行修正，风机轮毂处的实际风速可表示为v=vm()hhmb（3）式中，vm为气象站实际测量风速；hm为测量高度；h为风机轮毂高度；b为地面粗糙度因子，开阔、平坦地区一般为0.1439。1.2光伏阵列模型光伏发电系统输出功率主要受太阳辐射强度和及环境温度影响，其输出功率10表达式为Ppv=PNfGGSTC1+k()Ta-TSTC（4）式中，PN为标准测试条件下光伏的额定输出功率；f为光伏降额因数；G为光伏倾斜面的实际辐射强度；GSTC、TSTC为标准测试条件下的辐照强度和光伏电池温度；k为功率温度系数；Ta为光伏表面的实际温度。Ta表达式如下所示：Ta=Tr+301 000G（5）式中，Tr为实际环境温度。1.3储能模型本文选择锂电池储能。在并网型微电网中，锂电池可以将新能源出力高峰时的电量进行储存，当新能源的出力由于天气原因不满足负荷需求时，锂电池进行放电，从而提高了系统供电的可靠性。微网中锂电池充放电过程的数学模型11可表示为SOC()t=()1-SOC()t-1+Pchess()t tEssnc（6）SOC()t=()1-SOC()t-1-Pdisess()t tEssnd（7）式中，SOC()t、SOC()t-1为锂电池在t、t-1时刻的荷电状态；为电池自放电率；Pchess()t、Pdisess()t、c、d为锂电池储能在t时刻的充放电功率和充放电效率；Essn为电池储能的额定容量；t为t时间段长度，设定为1 h。2并网型微电网容量优化配置模型2.1目标函数风光储微电网目标函数为系统的年综合经济成本最小，其数学表达式为C=min()C1-C2（8）C1=Cinv+Crep+Com+Cp+Cbuy（9）第 4 期田雨明，等：基于改进灰狼算法的并网型微电网容量优化研究57第 19 卷沈阳工程学院学报（自然科学版）C2=Csell+Cd（10）式中，Cinv为微电网的年均投资成本；Crep为各设备的置换成本；Com为设备的年运行维护成本；Cp为年停电惩罚成本；Cbuy为年购电成本；Csell为年售电收益；Cd为设备的报废收益。2.1.1年均投资成本与年运行维护成本微电网的年均投资费用为Cinv=r()1+rY()1+rY-1i=1mNiCinv,iSi（11）式中，r为贴现率；Y为项目寿命周期；m为设备类型数目；Ni为第i类设备的数量；Cinv,i为第i类设备的单位容量投资成本；Si为第i类设备的单机容量。微电网的年运行维护成本为Com=i=1mNikiSi（12）式中，ki为第i类设备的单位容量年运维成本。2.1.2年置换成本微电网项目周期一般为20年，而风机、光伏的寿命均为20年，因此不需要考虑其置换成本。电池储能寿命一般不超过20年，所以本文按照标准工况下锂电池储能寿命为10年，考虑对其进行一次置换：Crep=NessCinv,essSessRr()1+rY()1+rY-1（13）R=1()1+rYess（14）式中，Ness为储能的数量；Cinv,ess为储能单位容量投资成本；Cinv,ess为单个储能电池的容量；R为置换系数；Yess为储能的寿命。2.1.3年购电成本与年停电惩罚成本当新能源出力不满足负荷需求时，需向主网购电，微电网的年购电成本为Cbuy=t=18 760cb()t Pgrid,b()t（15）式中，cb()t为t时段从主网购电的电价；Pgrid,b()t为t时段系统从主网的购电功率。虽然主网可以对微电网提供功率支撑，但某些时刻微电网内仍可能会出现供电不足的现象，需对其进行惩罚，惩罚成本为Cp=t=18 760lPloss()t（16）式中，l为停电惩罚系数；Ploss()t为t时刻不满足负荷的缺额功率。2.1.4年售电收益微电网的年售电收益包含上网售电收益和向负荷供电收益：Csell=Cin+Cin_v（17）Cin=t=18 760cs()tPload()t-Ploss()t（18）Cin_v=t=18 760cs()t Pgrid,s()t（19）式中，Cin和Cin_v分别为上网售电收益和向负荷供电收益；cs()t为t时刻微电网向主网和负荷售电的电价；Pload()t和Pgrid,s()t为t时刻负荷的需求和上网售电功率。2.1.5报废回收效益微电网各类设备的报废收益为Cd=()Cpd+Cwd+Cessdr()1+rY-1（20）Cpd=Cinv,pvNpvSpv（21）Cwd=Cinv,wtNwtSwt（22）Cessd=Ness()Cinv,essSess+RCinv,essSess（23）式中，Cpd、Cessd、Cwd分别为光伏、储能、风机报废时获得的收益；为设备回收系数，一般取 8%12；Cinv，pv、Cinv，wt、Cinv，ess分别为光伏、风机、储能的单位容量投资成本；Npv、Nwt、Spv、Swt分别为光伏、风机的数量和单机容量。锂电池储能需要进行一次替换，所以在模型中应考虑储能设备两次报废所获得的收益。2.2约束条件2.2.1微电网内各类设备的安装数量约束0 Npv Npv,max（24）0 Nwt Nwt,max（25）0 Ness Ness,max（26）58式中，Npv，max、Nw，max、Ness，max分别为光伏、风机、储能的最大装机数量。2.2.2锂电池储能约束锂电池储能在运行时应满足以下约束条件：SOC,min SOC()t SOC,max（27）0 Pchess()t Pchess,max（28）0 Pdisess()t Pdisess,max（29）式中，SOC,min、SOC,max分别为锂电池储能的最小荷电状态和最大荷电状态；Pchess,max、Pdisess,max分别为最大充电功率和最大放电功率。2.2.3并网点交换功率约束当微电网与主网进行能量交互时，应设置允许交换功率的限值。Pgrid,min Pgrid()t Pgrid,max（30）式中，Pgrid，min、Pgrid，max分别为微电网与主网进行功率交互的下限和上限。2.2.4微电网功率平衡约束为保证系统安全稳定运行，微电网内需要实时保持有功功率平衡。其约束表达式为Pwt()t+Ppv()t+Pdisess()t+Pgrid,b()t=Pload()t+Pchess+Pgrid,s()t（31）在系统运行时，锂电池不能同时充电和放电，微电网内也不能同时出现馈入主网的功率和主网提供的功率，即Pdisess()t Pchess()t=0（32）Pgrid,s()t Pgrid,b()t=0（33）2.3并网型微电网运行策略为充分利用微网内的清洁能源，并网型微电网在运行时，新能源发电设备应优先为区域内负荷提供所需电能。但风力发电和光伏发电存在不确定性和随机性，微网内将会出现不平衡功率，所以首先需要判断微电网在t时刻产生的不平衡功率：P=Pload()t-Pwt()t-Ppv()t（34）如果P 0，并且t时刻储能的荷电状态在规定范围之内，则优先通过锂电池储能放电为负荷供电；当光伏、风电和储能出力仍小于负荷需求时，则向大电网购电，即Pgrid,b()t=P-Pdisess()t（36）若大电网以最大联络线功率进行供电时，仍不满足负荷需求，则需要切掉非重要负荷以保证微网内的功率平衡，同时计算出不满足负荷的缺额功率并进行惩罚。不满足负荷的缺额功率为Ploss()t=P-Pdisess()t-Pgrid,max（37）3灰狼优化算法3.1灰狼优化算法原理及改进灰狼优化算法是一种群体智能优化算法，目前主要应用于参数优化、复杂函数优化、经济调度等领域13。该算法采用精英群体领导策略，在每次迭代中根据适应度值排名依次选取出前三名适应度最优的精英个体，即狼、狼、狼，其余个体则被定义为狼。在优化过程中，根据前三头狼的位置来判断猎物（全局最优解）的潜在位置，其余狼则必须服从前三头狼的引导，并不断更新自身位置向猎物靠近。灰狼更新自身位置的机制如图1所示。图1灰狼个体位置更新机制图1中，D、D、D分别为3只领头狼与其余灰狼个体的距离，可以表示为14第 4 期田雨明，等：基于改进灰狼算法的并网型微电网容量优化研究59第 19 卷沈阳工程学院学报（自然科学版）D=|C1X()n-X()n（38）D=|C2X()n-X()n（39）D=|C3X()n-X()n（40）Ck=2r1()k=1,2,3（41）式中，X()n、X()n、X()n为 3 只领头狼的当前位置；X()n为其余灰狼的当前位置；r1为 0，1 之间的随机数；Ck为系数向量，表示个体所在的位置对猎物影响的随机权重，有助于提升算法寻找全局最优解的能力。在确定领头狼与其余灰狼个体的距离之后，狼群将开始逐渐逼近包围猎物15。其余灰狼个体朝向3只领头狼前进的步长为X1=X()n-A1D（42）X2=X()n-A2D（43）X3=X()n-A3D（44）Ak=a()2r2-1(k=1,2,3)（45）a=2-2nnmax（46）其余个体的最终位置为X()n+1=()X1+X2+X33（47）式中，n为算法的迭代次数；nmax为算法最大迭代次数；a为收敛因子；Ak为系数向量；r2为 0，1 之间的随机数。从上述式子可以看出：随着迭代次数的增加，收敛因子线性递减至0，导致Ak的值也在对应区间内发生变化。当|Ak|1 时，灰狼将与局部最优分离，并扩大范围去寻找更优异的解；当|Ak|1时，狼群则向猎物发起攻击，搜索局部区域的最优解。虽然自适应的a与Ak保证了全局优化和局部优化，但由于寻优过程复杂，收敛因子线性递减不符合实际优化过程对Ak值的要求。本文在算法中采用了文献 16 提出的非线性收敛因子，其形式如下：a=amax-()amax-amin()n-1N-1cos()nN（48）式中，amax、amin分别为非线性收敛因子的上限和下限；N为迭代阀值。该非线性收敛因子在迭代次数小于迭代阀值时，收敛速度缓慢；而当大于迭代阀值时，非线性收敛因子为常数，既有助于增强算法的全局搜索能力，又能够进行细致的局部搜索，可以较好地均衡算法全局搜索和局部搜索性能。3.2微电网容量优化配置流程针对并网型风光储微电网容量优化模型，将风机、光伏组件和储能电池的数量设置为改进灰狼算法的决策变量，种群数目取 40，最大迭代次数取100，非线性收敛因子的上、下限取2和0.01，迭代阀值取40。具体优化流程如下：1）输入全年8 760 h风速、光伏倾斜面的实际辐射强度、环境温度、历史负荷及各类设备参数等数据。2）依据风光储容量的上下限约束，在搜索空间中随机产生初始种群，并将目标函数转换为适应度函数。3）在运行策略的作用下，根据气象数据和分布式电源、储能的出力模型计算各单元的实时出力，并通过计算适应度函数求出每只灰狼个体的适应度值。4）根据每个灰狼个体的适应度值排名，获取适应度值最优的、狼。5）根据式（48）计算非线性收敛因子，并更新其他狼的位置。6）判断是否达到最大迭代次数nmax，若未达到，则重复步骤3）步骤5）；若达到，则输出最优灰狼个体。4算例分析以北纬4045、东经12051 地区的某区域作为研究对象，利用Meteonorm气象软件获得该地区的全年气象数据。图2为利用该软件计算得到的光伏倾斜面的辐射强度，在设置36倾斜角后，春季60和冬季中光伏所接收的辐射量可以得到显著提升。图3为该地区温度变化数据，图4为风机轮毂高度处的时序风速，是根据式（3）对原始数据修正所得。图5为全年用电负荷数据，从图中可以看出该地区在夏季的用电需求最大。各类设备的成本数据如表1所示17。分时电价信息如表2所示，微电网向负荷售电的电价和上网电价相同。最后采用文中方法对微电网进行容量优化，仿真时间为8 760 h，步长设置为1 h。锂电池储能装置SOC的上、下限分别为90%和10%，初始状态为50%。图2光伏阵列倾斜面上的辐照度图3全年风速数据图4全年温度数据图5全年负荷数据表1各类设备成本数据元/kW类型风机光伏储能单机容量10 kW1 kW2 kWh初始投资成本6 0003 7401 500年运维成本405118表2微电网的售电电价和购电电价元/（kW h）时段0：007：007：0010：0015：0018：0021：0024：0010：0015：0018：0021：00售电电价0.2200.4200.650购电电价0.2640.5040.7804.1不同优化算法下优化结果对比分别采用灰狼算法和改进的灰狼算法（Improved Grey Wolf Optimization，IGWO）对微电网容量优化配置模型进行求解，两种算法的种群规模均取40，最大迭代次数为100。两种算法的迭代寻优曲线如图 6 所示，从图中可以看出：虽然 IGWO算法前期收敛速度与 GWO 算法差别不大，但IGWO算法的优化精度高于GWO算法，IGWO算法求出的最优目标函数值为-9 406.86，而 GWO算法得到的最优目标函数值为-8 475.23。这是由于IGWO算法采用的非线性收敛因子增强了算法前期的全局搜索性能，较好地提升了算法的寻优精度。第 4 期田雨明，等：基于改进灰狼算法的并网型微电网容量优化研究61第 19 卷沈阳工程学院学报（自然科学版）图6两种算法的收敛曲线4.2不同停电惩罚系数下的优化结果分析为进一步验证本文所提微电网容量优化配置模型的准确性，分别设置停电惩罚系数为2、3、4，并采用IGWO算法进行仿真求解，得到的配置结果如表3所示，优化运行结果如表4所示。从表3和表4可以看出：当惩罚系数设置为2时，微电网的年综合经济成本为负值，说明此时微电网运营商可以通过向负荷供电、余电上网两种方式获得一定的投资收益；但随着停电惩罚系数的逐渐增大，年综合经济成本均变为正值，这是因为停电惩罚系数增大所造成的经济损失急剧增加，为保证系统运行的经济性，需要配置更多的分布式电源来弥补系统中的功率缺额，以减少停电惩罚费用。由表3可知，由于风机的成本较高，导致风机的容量也相对较少，而储能与光伏的安装容量较大。由表4可以看出：惩罚系数设为3、4时，微电网的停电惩罚成本、购电成本均有所减少，售电收益有所增加，微电网的供电可靠性得到提升，但同时也导致微电网的初始投资成本、储能系统的置换成本以较大幅度增加，所以微电网所能产生的经济效益将会逐渐降低。表3不同停电惩罚系数的优化配置结果惩罚系数234风机/kW390360340光伏/kW7449131 005储能/(kWh)7881 3841 656年综合经济成本/万元-0.946.2011.13表4不同停电惩罚系数下的优化运行结果万元惩罚系数234初始投资成本46.3956.3061.26运行维护成本6.778.599.48置换成本5.8810.3212.32停电惩罚成本18.8516.9817.02报废收益1.912.432.69售电收益91.6292.9993.78购电成本14.709.437.525结论研究并网型微电网容量优化对减小系统成本和提升系统运行的稳定性具有重要意义。为此，从微电网初始投资成本、运行维护成本、停电惩罚成本、替换成本、报废收益及售购电成本等经济性角度提出了一种并网型微电网容量优化配置方法。在分析各分布式电源出力特性的基础上，建立以年综合经济成本最小为目标函数，以功率平衡约束、储能系统约束和联络线交换功率约束为条件的微电网容量优化配置模型，采用改进的灰狼算法对该模型进行求解。通过对比分析不同停电惩罚系数对容量配置结果的影响可知：适当减少停电惩罚系数可降低微电网的年综合经济成本，验证了所提出的优化配置模型的有效性。但由于模型中未考虑可靠性指标，后续工作将进一步研究并完善满足可靠运行指标的并网型微电网容量优化配置模型。参考文献1孟明，陈世超，赵树军，等.新能源微电网研究综述 J.现代电力，2017，34（1）：1-7.2 HATZIARGYRIOU N，ASANO H，IRAVANI R，et a1.Microgrids：an overview of ongoing research，development and demonstration projectsJ.IEEE Power andEnergy Magazine，2007，5（4）：78-94.3王彦宇，郭权利.微电网示范工程综述 J.沈阳工程学院学报（自然科学版），2019，15（1）：82-87.4于凯旋，李凤霞，李双拾.基于改进BASDE算法的孤立微电网多目标容量优化配置仿真 J.电网与清洁能62源，2021，37（8）：109-117.5 DIAB AA Z，SULTAN H M，MOHAMED I S，et al.Application of different optimization algorithms for optimalsizing of PV/wind/diesel/battery storage stand-alone hybrid microgrid J.IEEE Access，2019，7：119 223-119245.6任岷.并网海岛微电网容量优化配置研究 D.大连：大连理工大学，2021.7崔明勇，杨林林，杨少华，等.基于需求侧响应的水风互补微电网容量优化配置 J.电工电能新技术，2022，41（2）：21-34.8张节谭，程浩忠，胡泽春，等.含风电场的电力系统随机生产模拟 J.中国电机工程学报，2009，29（28）：34-39.9茆美琴，孙树娟，苏建徽.包含电动汽车的风/光/储微电网经济性分析 J.电力系统自动化，2011，35（14）：30-35.10王成山，洪博文，郭力.不同场景下的光蓄微电网调度策略 J.电网技术，2011，37（7）：1 775-1 782.11 丁明，王波，赵波，等.独立风光柴储微网系统容量优化配置 J.电网技术，2013，37（3）：575-581.12 田梁玉，唐忠，焦婷，等.基于多目标评判法的微电网经济性指标研究 J.电力系统保护与控制，2016，44（19）：122-127.13 李春雷，王洪江，尹常永，等.基于改进灰狼算法的微电网优化运行 J.沈阳工程学院学报（自然科学版），2023，19（2）：50-56.14 张晓凤，王秀英.灰狼优化算法研究综述 J.计算机科学，2019，46（3）：30-38.15张涛，张东方，王凌云.基于灰狼算法的分布式电源优化配置 J.水电能源科学，2018，36（4）：204-207.16 蔡国伟，刘旭，张旺，等.基于改进灰狼优化算法的分布式电源优化配置 J.太阳能学报，2019，40（1）：134-141.17 王晓毅，唐忠.考虑供需自平衡和独立运行能力的并网型微电网容量优化配置 J.太阳能学报，2021，42（5）：74-82.Research on Capacity Optimization of Grid-connected MicrogridBased on Improved Gray Wolf AlgorithmTIAN Yuminga,YIN Changyongb（a.School of Electric Power;b.School of Renewable Energy,Shenyang Institute of Engineering,Shenyang 110136,Liaoning Province）Abstract：In order to rationally allocate the capacity of distributed power sources in grid-connected microgridand reduce the economic cost of microgrid,the optimal capacity allocation method of grid-connected microgridis studied.The mathematical model of each unit in grid-connected microgrid is established,and the operationstrategy of microgrid is put forward.Secondly,the objective function is established by comprehensivelyconsidering the systems average annual investment cost,annual operation and maintenance cost,replacementcost,annual power outage compensation fee,annual power purchase cost,annual operation income andequipment scrap recovery benefit,and combining the power balance constraints of the microgrid,the operationconstraints of the equipment and the exchange power constraints with the main network to build the capacityoptimization model of the grid-connected microgrid.Then the improved gray Wolf algorithm was used to solvethe optimization model.Finally,the effectiveness of the proposed model is verified by simulation analysis.Keywords：Microgrid;Grid-connected;Capacity configuration;Improved grey Wolf optimization algorithm第 4 期田雨明，等：基于改进灰狼算法的并网型微电网容量优化研究63