一种孤岛混合可再生能源系统的多目标组合优化配置方法
一种孤岛混合可再生能源系统的多目标组合优化配置方法
【技术领域】
[0001] 本发明涉及一种孤岛混合可再生能源系统的多目标组合优化配置方法，具体而言 是考虑能源系统中的各单元的类型和数量的配置，最优化系统的各个目标，属于计算机领 域。
【背景技术】
[0002] 能源问题始终伴随着人类文明发展的整个进程，特别是进入21世纪后，工业经济 的急速发展对能源需求提出了更高的要求。然而煤、石油之类的化石能源日益枯竭，并且使 用化石燃料带来的环境问题日益凸显，能源安全和环境污染已经成为当前世界各国面临的 巨大挑战。为了应对能源和环境的双重危机，大力发展和利用可再生能源，如太阳能、风能 等，成为当务之急。
[0003]混合可再生能源系统（Hybrid Renewable Energy Systems,HRES)顾名思义，是一 种混合使用化石能源和可再生能源的系统。具体而言是包括风力发电、光伏发电、柴油发电 以及储能装置的能源系统。孤岛混合可再生能源系统指HRES没有接入系统电网。孤岛混合 可再生能源系统通常用于海岛之类的偏远地区。用户相对较少，用电量需求较低，引入大电 网比较困难且经济成本较高。
[0004] 孤岛混合可再生能源系统规划需要解决的问题是：如何合理配置风力机、光伏板、 柴油机以及储能设备以使得整个能源系统满足区域用能需求，同时经济成本最低、对环境 危害最小。
[0005] 目前关于孤岛混合可再生能源系统的组合优化配置研究主要集中于单目标的优 化模型，如在满足供需平衡下，最小化HRES全寿命周期成本或最小化HRES的温室气体排放 量等。然而，从实际角度出发，HRES的组合优化配置需要同时考虑全寿命周期成本、温室气 体排放量等评估指标，是一个多目标优化问题。同时，HRES的组合优化配置中涉及到风光柴 储设备的选择，柴油使用量等变量的规划，是一个多类型变量的（连续、离散变量)优化问 题。此外，HRES组合优化配置还需满足供需平衡约束。综上所述，HRES的组合优化配置是一 个多变量、多目标、带约束的复杂优化问题，目前关于HRES优化配置的相关研究(包括模型 及求解算法)远不能满足实际需求。
[0006] 所述的多目标优化问题是指：同时对多个目标进行优化，由于各个目标之间通常 是耦合在一起且互相制约、互相竞争，即某个目标的改善可能引起其他目标性能的降低，很 难找到一个真正意义上的最优解使得各个目标同时达到最优，因此多目标优化问题的最优 解通常不是单一的，而是一组多个互有利弊的非支配解，是一个非劣解的集合，即帕累托 (Pareto)最优解集。求解多目标优化问题的核心是找到一组分布均勾的Pareto最优解。
[0007] 所述的智能多目标优化算法是指:传统的处理多目标优化问题的方法，如加权法、 约束法、目标规划法等，通过构建一个评价函数，将多目标问题转化为单目标优化问题，然 后利用一般的求解方法计算得到问题的一个解。智能多目标优化算法是通过模拟某些自然 过程发展而来的基于种群的优化算法，其思想和内容涉及数学、生物学和计算机学科等。该 类算法不依赖于梯度信息，一次运行能够找到一组Pareto最优解，具有全局、并行、高效、鲁 棒和通用性强等特点。是求解复杂非线性多目标优化问题的有效方法。
【发明内容】
[0008] 本发明要解决的技术问题是如何合理配置风力机、光伏板、柴油机以及储能设备 以使得整个孤岛混合可再生能源系统满足区域用能需求，且经济成本最低、环境危害最小。
[0009] 为解决该问题，本发明所采取的技术方案是:建立HRES的组合优化配置模型，利用 智能多目标优化算法求解，得到一组分布均勾的Pareto最优解，从这组Pareto最优解中，选 择一个作为最终配置方案来对HRES进行配置。
[0010]步骤1:以HRES全寿命周期成本最低，HRES对环境危害最小为目标，以满足区域能 源需求为约束条件，建立HRES的多目标组合优化配置模型：
[0011]
[0012]
[0013]
[0014] 式（1)表示以最小化HRES的全寿命周期成本和系统的C〇2排放量为目标，以系统供 能满足用户用电需求为约束条件，构造出带约束的两目标优化模型；
[0015] 其中Fccist(x)表示HRES配置方案为X下的全寿命周期（T年）成本，包括孤岛风光柴 储能源系统中各类设备的初始安装费用和使用过程中的维护费用；
[0016] &〇2(Χ)表示HRES配置方案为X下系统的 C〇2排放量，包括HRES中风光柴储设备在生 产制造过程的c〇2排放量和柴油机使用过程中的c〇2排放量，其中风力机、光伏板以及储能设 备在使用过程中的c〇2排放量不予考虑；
[0017] Fsupply(x)表示HRES配置方案为X下的供电量；
[0018] Fdemand表示某一地区的用电需求量；
[0019] 所述配置方案X是指HRES系统中光伏组件、风力机、储能设备、柴油机四类组件的 配置数量以及柴油使用量，配置方案X表示如下：
[0020]
[0021] xi表示HRES中第i类型的光伏板个数，ie{l，2，. . .，NPV}; Vi表示HRES中第j类型 的风力机个数，je{l，2, . . .，Nwt};xvv^表示HRES中第k类型的储能设备个数，ke{l， 2，. . .，Nba}; 表示HRES中第 1 类型的柴油机个数，1 e {1，2，. . .，Nfg}; 示HRES中第1类型的柴油机的柴油使用量，单位为升，取值范围为[0，C]，C表示柴油使用量 的上限值，对于柴油发电机，若未使用某一类型的柴油机，那么其对应的柴油使用量设定为 〇;配置方案X的前Npv+N wt+Nba+Nfg个变量取值范围为(0，Nl]，Nl表示设备个数的上限值;N PV表 示光伏板的类型数，Nwt表示风力机的类型数，他3表示备选储能设备的类型数，心 8表示柴油 机的类型数；
[0022] Fccist(x)表示为：
[0023]
[0024] Fc〇stPV(x)表不腿S配置方案为x下的光伏板全寿命周期费用，匕燃(》) = 1:'-; <?", + C2n， Fcos· (x)表不HRE S配置方案为x下的风力机全寿命周期费用，4tWT W = Σ ,1卜n?β+ ""r) * 心他⑴表示TOS配置方案为x下的备选储能设备全寿命周期费用，，; FmstreW^^^RES配置方案为x下的柴油机周期费用， Clpvi表不第i类型光伏板的初装费用，C2pvi表不第i类型光伏板的维护费用；ClwTj表不第j类 型风力机的初装费用，C2叩表示第j类型风力机的维护费用;ClBAk表示第k类型备选储能设 备的初装费用，C2 BAk表示第k类型备选储能设备的维护费用;C1FC1表示第1类型柴油机的初 装费用，C2?表示第1类型柴油机的维护费用;C FG表示柴油价格；
[0025]
[0026]
[0027] Fc〇2PV(x)ifcrii?Sl￡Kir案为 X 下的光伏 表不HRES配置方案为X下的风力机全寿命周期碳排放量，= ; .&α2ΒΑ(χ_.)_ 表示HRES配置方案为X下的备选储能设备全寿命周期碳排放量，&2βα> = ΣΓ:;W
表不HRES配置方案χ下的柴油设备全寿命周期碳排放量，+~％.??.*+?.，*<)* Pm表示第i类型光伏板的额定功率;〇PVl表示制造单位面积的第i类型光伏板产生的碳排放 量;Own表示制造第j类型风力机产生的碳排放量;〇BAk表示制造第k类型备选储能设备产生 的碳排放量;〇 FC1表示第1类型柴油发电设备产生的碳排放量;表示燃烧每升柴油的碳排 放量；
[0028]
[0029]
[0030] Fsuppiypv(x)表不HRES配置方案为X下的光伏板发电量，.FsupplyPV(x) = i：、； L为HRES建设地区年平均光照强度，单位为勒克斯，相关数据可以通过查询当地历史天气数 据分析得到；FsupplyWT(X)表不HRES配置方案为X下的风力机发电量， AUpplyWT(X) = * Γ * ; V为年平均风速，单位为米/秒，相关数据可以通过查询当 地历史天气数据分析得到；FsupplyBA(x)表示HRES配置方案为X下的备选储能设备供电量， 巧upplyBA(X) = Σ?? * ' J :; FsupplyFG ( X )表不HRES配置方案为X下的柴油设备发电直， W / V ) - \ ^ ，Γ ρ '1·· γ 承：y .? Siippl5Hrv ~ Lui-^X · . -十/ ?卡1+"
[0031] PPVl表示第i类型光伏板的额定功率;PWTj表示第j类型风力机的额定功率，P BAk表示 第k类型备选储能设备的额定功率，ρ?表示第1类型柴油机的额定功率；
[0032] Fdemand表不为：
[0033] Fd_nd = Q，其中Q为某地区年平均用电量为Q(瓦）；
[0034] 步骤2:利用智能多目标优化算法NSGA-II求解HRES组合优化模型，具体过程如下：
[0035] 步骤2.1:设置算法参数:包括种群规模和终止条件，种群规模为N，取值范围100至 500，终止条件采用最大运行代数maxgen，取值范围50至500;
[0036] 步骤2.2:初始化种群：随机生成种群规模N= 100父代种群S;种群中每个个体X包 含η个编码值，其中n = Npv+Nwt+Nba+2*Nf g，令当前运行代数gen = 1;
[0037] 步骤2.3:基于当前种群S，通过交叉变异操作产生子代种群Sc ，种群规模也为N。所 述交叉变异操作实施方法如下：
[0038]步骤2.3.1:交叉操作；
[0039]步骤2.3.1.1:令q = 1，并初始化集合Sc '为空集；
[0040] 步骤2.3.1.2:选取个体12。-1和12。^2。-1三3，12。三3，，首先随机生成一个位于区间
[l，n]之间的整数，作为交叉位置，记为i'，然后用X2H的前半段第1到第i'个编码和^。的后 半段即第i'+l到最后一个编码组成子个体的后半段即第i'+l到最后一个编码和 x2q的前半段即第1到第i'个编码组成子个体x'2q;将产生的个体^2< 1-1和1'2(1存入集合3(3'， 艮P Sc ' = Sc ' U X ' 2q-1 U X ' 2q ;
[0041 ] 步骤2.3.1.3:令9 = 9+1;若9>~/2，转向步骤2.3.2，否则转向步骤2.3.1.2;
[0042] 步骤2.3.2:变异操作，设置变异操作的概率为σ，σ e (〇，1);
[0043] 步骤2 · 3 · 2 · 1:令q = 1，初始化Sc为空集；
[0044] 步骤2 · 3 · 2 · 2 :针对个体xqe Sc '，随机生成一个位于(0，1)区间的实数，记为rand。 若rand〈〇,则随机选取一个变异位置v，ve [1，n]，若ve [1 ,η-Nfg]，则随机选取[1，Nl]区间 的一个正整数代替个体Xq的第V个变量值，得到变异个体X'q,并存入集合SC，即SC = SC U X 'q;若ve[n-Nfg+l，n]，则随机选取位于(0，C]区间的实数代替个体 Xq的第v个变量值，得到 变异个体X ' q，并存入集合Sc，即Sc = Sc U X ' q;若rand 2 〇，将Xq存入集合Sc，即Sc = Sc U Xq; [0045] 步骤2.3.2.3:令q = q+l;若q>N，转向步骤2.4,否则转向步骤2.3.2.2;
[0046]步骤2.4:将父代种群5与子代种群5(：合并，得到规模为2_勺合种群5311 = 51^(3，对 合种群Sall中的个体进行非劣分层，然后计算每一个非劣层的个体局部拥挤距离，最后依据 个体所处的非劣层及拥挤距离，对所有个体进行排序；
[0047] 步骤2.4.1:对种群Saii中的个体进行非劣分层。下面介绍非劣分层的方法；
[0048] 步骤 2.4.1.1:令h = l;
[0049] 步骤2.4.1.2:找出种群5311中不被任何个体Pareto约束支配的个体，并保存在集 合Ah中，即集合Ah中的个体均属于第h非劣层;所述个体x Pareto约束支配个体y，是指当下 列其中一个条件满足时：（i)个体X和y均满足约束条件且κ y; (i i)个体X满足约束条件，而 y不满足约束条件；（i ii)个体X和y均不满足约束条件，个体X违反约束条件的程度小于个体 y违反约束条件的程度，则称个体X Pareto约束支配个体y;所述的关系式x *<y表示个体X Pareto支配个体y，是在所有目标函数值上，个体X的目标函数值都不大于个体y的目标函数 值即VX，f m(xHfm(y)，并且至少存在一个目标函数上个体X的目标函数值小于个体y的目 标函数值，即3九.上令)〈匕，（7)，所述匕，111=1，2为步骤1的公式（1)中的两个目标函数4 1表 示目标函数F?st(x)的值，5表示目标函数巧。2(幻的值；
[0050] 步骤2.4.1.3 : WSaii中除去所有在集合Ah中的个体，判断当前Saii是否为空集，若 是，则分层完毕，转入步骤2.4.2;否则转入步骤2.4.1.2，令h = h+l;
[0051] 步骤2.4.2:计算种群中个体的拥挤距离，拥挤距离越小，个体周围越稠密;拥挤距 离具体计算方法如下：
[0052] 步骤2.4.2.1:计算种群Sall内每个个体的目标函数值，并对其进行排序;所述目标 函数值指在步骤1的公式(1)中的目标函数F? st(x)和匕):(x)的值；
[0053] 步骤2.4.2.2:计算拥挤距离;首先定义在种群中目标函数心或5最小的个体为边 界个体，边界个体的拥挤距离dist定义为无穷大；其次计算除边界外的个体^拥挤距离 dist(q)为
[0054]
[0055] 其中./Γχ和./Γ分别表示当前种群中目标函数匕的最大和最小值；/T1和表 示个体X(rl和个体Xq+1的目标函数fVf直，fm=1表示目标函数Fc^stU)的值，fm=2表示目标函数 &>)的值；
[0056]步骤2.4.3:基于个体所处的非劣层及拥挤距离，对所有个体进行排序，排序准则 为:处于第h非劣层的个体排序在第h '非劣层的个体之前，其中h ' >h;若两个个体处于同一 非劣层，那么拥挤距离大的个体排序在拥挤距离小的个体之前；
[0057]步骤2.5:在排序之后的个体种群中选取第1到N个个体作为新的父代种群S，令gen = gen+l ；
[0058]步骤2 · 6:判断gen是否大于maxgen，若大于，则输出种群S中的处于第1非劣层中的 个体作为步骤1中所确定的多目标优化模型的Pareto最优解集，转步骤3;否则转步骤2.3。 [0059]步骤3:从Pareto最优解集中选择最优HRES方案；
[0060] 步骤3.1当要求整个系统的全寿命周期成本小于等于给定值CostA时，选择Pareto 最优解集中目标函数值6:+00最小的解，即使系统C02排放量最小的个体X，作为最终的配 置方案X，所述CostA指HRES配置过程所要求的最大成本额度；
[0061 ]步骤3.2当要求整个系统的⑶2排放量小于等于给定值⑶2A时，选择Pareto最优解 集中目标函数值F_t(X)最小的解，即使系统全寿命周期成本最小的个体X，作为最终的配 置方案X，所述C0 2#BHRES配置过程所要求的最大排放量；
[0062]本发明的优点在于：
[0063] (1)本发明针对HRES中设备的组合配置构建了带约束的多目标组合优化模型，更 符合实际情况、方案可行性更强。
[0064] (2)本发明采用多目标智能优化算法求解问题，能够同时找到一组Pareto最优解 集，满足不同情况(即，对成本和环境污染的不同要求)下的HRES方案配置。
【附图说明】
[0065] 图1典型孤岛混合可再生能源系统示意图
[0066]图2是利用智能多目标优化算法NSGA-II求解HRES组合优化配置问题的流程图
[0067] 图3 HRES方案设计编码
[0068]图4单点交叉和均匀变异算子示意图
[0069]图5个体非劣分层及拥挤距离计算示意图
[0070]图6 HRES规划的Pareto最优解集
[0071] 图7 HRES规划的决策过程示意图
【具体实施方式】
[0072] 下面结合附图和实例对本发明作进一步的详细说明。
[0073] 本发明提供的孤岛混合可再生能源系统的多目标组合优化配置方法，具体步骤如 下：
[0074]表1至4给出了HRES系统（如图1所示）中各设备的相关数据，同时已知在HRES建设 地区年平均光照强度L = 20(单位:勒克斯）、年平均风速V = 5(单位:米/秒）以及年平均用电 量Q=1000(单位:千瓦）；
[0075] 表1备选光伏板相关参数(假设每块光伏板面积为1米2)
[0076]
[0077]表2备选风力机相关参数 [0078]
[0079] 表3备选储能设备相关参数
[0080]
[0081 ]表4备选柴油机相关参数 [0082]
[0083] 步骤1:考虑HRES的全寿命周期T = 20年，以HRES全寿命周期成本 [0084]最低，系统排放C02量最小为目标，以满足区域能源需求为约束条件，
[0085]建立HRES的多目标组合优化配置模型；
[0086] mLu(Fc〇st(x),FCOi(x)!
[0087] S.t.Fsupply(x) ^ Fdemand
[0088] 步骤2:利用智能多目标优化算法NSGA-II求解HRES组合优化模型，算法流程见图 2,具体过程描述如下：
[0089] 步骤2.1:设置算法参数:包括种群规模和终止条件，种群规模N设置为100，终止条 件采用最大运行代数maxgen，设置为maxgen = 100。
[0090] 步骤2.2:初始化种群:令当前运行代数8的=1。随机生成~=100父代种群5。种群 中每个个体X包含η个编码值，其中n = Npv+Nwt+Nba+2*Nfg，且N L取5，C取100;由表1 -4数据得 到，本实例中所考虑的光伏板类型数Npv = 5,风力机类型数Nwt = 5,储能设备类型数Nba = 3， 柴油发电机类型数Nfg = 5,个体编码长度n = 5+5+3+2*5 = 23,即每个个体包括23个编码值。 图3给出了针对此算例的一个个体编码。
[0091]步骤2.3:基于当前种群S，通过交叉变异操作（如图4所示)产生子代种群Sc，种群 规模也为N。所述交叉变异操作实施方法如下：
[0092]步骤2.3.1:交叉操作；
[0093]步骤2 · 3 · 1 · 1:令q = 1，并初始化集合Sc '为空集；
[0094] 步骤2.3.1.2:选取个体12。-1和12。^2。-1三3，12。三3，，首先随机生成一个位于区间 [l，n]之间的整数，作为交叉位置，记为i'，然后用X2 H的前半段第1到第i'个编码和^。的后 半段即第i'+l到最后一个编码组成子个体的后半段即第i'+l到最后一个编码和 x2q的前半段即第1到第i'个编码组成子个体x'2q;将产生的个体x'kl·和X' 2q#入集合Sc_l， 艮P Sc ' = Sc ' U X ' 2q-1 U X ' 2q ;
[0095] 步骤2.3.1.3:令9 = 9+1;若9>~/2，转向步骤2.3.2，否则转向步骤2.3.1.2;
[0096] 步骤2.3.2:变异操作，设置变异操作的概率为σ，〇6(0，1)，取〇 = 0.01;
[0097] 步骤2 · 3 · 2 · 1:令q = 1，初始化Sc为空集；
[0098] 步骤2 · 3 · 2 · 2 :针对个体xqe Sc '，随机生成一个位于(0，1)区间的实数，记为rand。 若rand〈〇,则随机选取一个变异位置v，ve [1，n]，若ve [1 ,η-Nfg]，则随机选取[1，Nl]区间 的一个正整数代替个体Xq的第V个变量值，得到变异个体X'q,并存入集合SC，即SC = SC U X 'q;若ve[n-Nfg+l，n]，则随机选取位于(0，C]区间的实数代替个体 Xq的第v个变量值，得到 变异个体X ' q，并存入集合Sc，即Sc = Sc U X ' q;若rand 2 〇，将Xq存入集合Sc，即Sc = Sc U Xq;
[0099] 步骤2 · 3 · 2 · 3:令q = q+1;若q>N，转向步骤2 · 4，否则
[0100] 转向步骤2.3.2.2;
[0101]步骤2.4:将父代种群3与子代种群5(：合并，得到规模为2_勺合种群3311 = 31^(3，对 合种群Sall中的个体进行非劣分层，然后计算每一个非劣层的个体局部拥挤距离，最后依据 个体所处的非劣层及拥挤距离，对所有个体进行排序；
[0102] 步骤2.4.1:对种群Sall中的个体进行非劣分层。下面介绍非劣分层的方法；
[0103] 步骤 2.4.1.1:令h = l;
[0104] 步骤2.4.1.2:找出种群5311中不被任何个体Pareto约束支配的个体，并保存在集 合Ah中，即集合Ah中的个体均属于第h非劣层;所述个体X Pareto约束支配个体y，是指当下 列其中一个条件满足时：（i)个体X和y均满足约束条件且个体X满足约束条件，而 y不满足约束条件；（i ii)个体X和y均不满足约束条件，个体X违反约束条件的程度小于个体 y违反约束条件的程度，则称个体X Pareto约束支配个体y;所述的关系式χ ?表示个体X Pareto支配个体y，是在所有目标函数值上，个体X的目标函数值都不大于个体y的目标函数 值即f m(y)，并且至少存在一个目标函数上个体X的目标函数值小于个体y的目 标函数值，即3人^七)仏，（7)，所述匕，111=1，2为步骤1的公式(1)中的两个目标函数4 1表 示目标函数F?st(x)的值，f2表示目标函数的值；
[0105] 步骤2.4.1.3 : Wsaii中除去所有在集合Ah中的个体，判断当前Sail是否为空集，若 是，则分层完毕，转入步骤2.4.2;否则转入步骤2.4.1.2，令h = h+l;
[0106] 步骤2.4.2:计算种群中个体的拥挤距离，拥挤距离越小，个体周围越稠密;拥挤距 离具体计算方法如下：
[0107] 步骤2.4.2.1:计算种群Sall内每个个体的目标函数值，并对其进行排序;所述目标 函数值指在步骤1中的目标函数F? st(x)和/*& (>0的值；
[0108] 步骤2.4.2.2:计算拥挤距离;首先定义在种群中目标函数心或5最小的个体为边 界个体，边界个体的拥挤距离dist定义为无穷大；其次计算除边界外的个体^拥挤距离 dist(q)为
[0109]
[0110] 六r T Jrn yj =| HU TT群中目标函数fm的最大和最小值；尤丨__和/0丨：表 示个体X(rl和个体Xq+1的目标函数fVf直，fm=1表示目标函数Fc^stU)的值，fm=2表示目标函数 匕K(X)的值；
[0111] 步骤2.4.3:基于个体所处的非劣层及拥挤距离，对所有个体进行排序，排序准则 为:处于第h非劣层的个体排序在第h '非劣层的个体之前，其中h ' >h;若两个个体处于同一 非劣层，那么拥挤距离大的个体排序在拥挤距离小的个体之前；
[0112] 步骤2.5:在排序之后的个体种群中选取第1到N个个体作为新的父代种群S，令gen = gen+l ；
[0113]步骤2 · 6:判断gen是否大于maxgen，若大于，则输出种群S中的处于Αι集合（即第1 非劣层中）的个体作为步骤1中所确定的多目标优化模型的Pareto最优解集，转步骤3;否则 转步骤2.3。
[0114]步骤3:从Pareto最优解集中选择最优HRES方案：依据本实例中建立的模型，利用 智能多目标优化算法NSGA-II求解模型，得到如图6所示的一组Pareto最优解集。
[0115]步骤3.1:当整个系统的全寿命周期成本小于等于给定值CostA时，选择Pareto最 优解集中目标函数值^:(x)最小的解，即使系统C02排放量最小的个体X，作为最终的配置 方案X，如图7所示，所述CostA指HRES配置过程所要求的最大成本额度。当要求HRES系统的 全寿命周期成本小于5.26*104千元，即CostA=5.26*104，HRES系统配置为使用36个#1类型 光伏板，19个#2类型光伏板，不使用#3、#4和#5类型的光伏板;使用41个#1类型风力机，22 个#3类型风力机，14个#4类型风力机，不使用#2、#5类型的风力机;使用26个#1类型储能设 备，不使用#2和#3类型的储能设备;使用4台#2类型柴油机，每台使用柴油183升，不使用其 他类型的柴油机；
[0116] 步骤3.2:当整个系统的C02排放量小于等于给定值C02A时，选择Pareto最优解集中 目标函数值F_ t(x)最小的解，即使系统全寿命周期成本最小的个体X，作为最终的配置方 案X，如图7所示。所述C02#gHRES配置过程所要求的最大排放量。当要求HRES系统的C0 2排放 量小于0.9*104克，即C02A=0.9*104，HRES系统配置为使用69个#1类型光伏板，不使用#2、# 3、#4和#5类型的光伏板;使用72个#1类型风力机，31个#2类型风力机，不使用#3、#4、#5类型 的风力机;使用52个#1类型储能设备，12个#2类型储能设备，不使用#3类型的储能设备;使 用1台#1类型柴油机，使用柴油11升，不使用其他类型的柴油机；
[0117] 本发明不仅限于上述【具体实施方式】，本领域一般技术人员根据本发明公开的内 容，可以采用其它多种具体方式实施本发明，如采用其他类似的多目标智能优化算法（如 MOGA、MOEA/D)。凡是采用本发明的技术思路，做一些简单的变化或改善的设计，如修改算法 参数等，都属于本发明保护的范围。
【主权项】
1. 一种孤岛混合可再生能源系统的多目标组合优化配置方法，其特征在于所述方法包 括以下具体步骤： 步骤1:建立HRES的多目标组合优化配置模型； 以HRES全寿命周期成本最低，HRES对环境危害最小为目标，以满足区域能源需求为约 束条件； 步骤2:利用智能多目标优化算法NSGA-II求解HRES组合优化模型，输出Pareto最优解 集； 步骤3:从Pareto最优解集中选择最优HRES方案； 步骤3.1当要求整个系统的全寿命周期成本小于等于给定值CostA时，选择Pareto最优 解集中目标函数值己,Jx)最小的解，即使系统CO2排放量最小的个体X，作为最终配置方案， 所述CostA指HRES配置过程所要求的最大成本额度； 步骤3.2当要求整个系统的C〇2排放量小于等于给定值C〇2A时，选择Pareto最优解集中 目标函数值Frast(X)最小的解，即使系统全寿命周期成本最小的个体X，作为最终配置方案， 所述CO 2WBHREs配置过程所要求的最大排放量。2. 根据权利要求1所述的一种孤岛混合可再生能源系统的多目标组合优化配置方法， 其特征在于所述的组合优化配置模型为：式（1)表示以最小化HRES的全寿命周期成本和系统的CO2排放量为目标，以系统供能满 足用户用电需求为约束条件，构造出带约束的两目标优化模型； 其中Frast(X)表示HRES配置方案为X下全寿命周期T年的成本，包括孤岛风光柴储能源系 统中各类设备的初始安装费用和使用过程中的维护费用； 々·〇_.(X)表示HRES配置方案为X下系统的CO2排放量，包括HRES中风光柴储设备在生产制 造过程的CO2排放量和柴油机使用过程中的CO2排放量，其中风力机、光伏板以及储能设备在 使用过程中的CO2排放量不予考虑； Fsupply(X)表示HRES配置方案为X下的供电量； Fd_nd表示某一地区的用电需求量； 所述配置方案X是指HRES系统中光伏组件、风力机、储能设备、柴油机四类组件的配置 数量以及柴油使用量，配置方案X表示如下：Xi表示HRES中第i类型的光伏板个数，ie{l，2，. ..，NPV}; V,表示HRES中第j类型的风 力机个数，je{l，2,…表示HRES中第k类型的储能设备个数，ke{l，2,…， Nba}; ' · w?表示HRES中第1类型的柴油机个数，I e {1，2，. . .，Nfg}; 表示HRES 中第1类型的柴油机的柴油使用量，单位为升，取值范围为[〇，C]，C表示柴油使用量的上限 值，对于柴油发电机，若未使用某一类型的柴油机，那么其对应的柴油使用量设定为0;配置 方案X的前Npv+Nwt+N ba+Nfg个变量取值范围为(0，Nl]，Nl表示设备个数的上限值;N pv表示光伏 板的类型数，Nwt表示风力机的类型数，他3表示备选储能设备的类型数，阶 8表示柴油机的类 型数； Fcost ( X )表为：Fcdsir<x)表不HRES配置方案为X下的光伏板全寿命周期费用，FccsfflT(X)表示HRES配置方案为X下的风力机全寿命周期费用，FcosM(X)表示HRES配置方案为X下的备选储能设备全寿命周期费用，-00????????%下的柴槲/1^#^周期费甩Clm表示第i类型光伏板的初装费用，C2PVl表示第i类型光伏板的维护费用;Ch^表示第j类 型风力机的初装费用，C2叩表示第j类型风力机的维护费用;ClBAk表示第k类型备选储能设 备的初装费用，C2 BAk表示第k类型备选储能设备的维护费用;CIfc1表示第1类型柴油机的初 装费用，C2ra表示第1类型柴油机的维护费用;C fg表示柴油价格； Go2 表不为表不HRES配置方案为X下的光伏板全寿命周期碳排放量，Ao2WT(X)表不HRES配置方案为X下的风力机全寿命周期碳排放量，Α〇2βα(χ)表示HRES配置方案为X下的备选储能设备全寿命周期碳排放量，砰㈤表示HRES配置方案X下的柴油设备全寿命周期碳排放量，Pm表示第i类型光伏板的额定功率;Om表示制造单位面积的第i类型光伏板产生的碳 排放量;0叩表示制造第j类型风力机产生的碳排放量;〇BAk表示制造第k类型备选储能设备 产生的碳排放量;Ofc1表示第1类型柴油发电设备产生的碳排放量;B wcl表示燃烧每升柴油的 碳排放量； Fsu卯Iy(X)表不为：Fsuppiypv(x)表不HRES配置方案为X下的光伏板发电量，；:L为 HRES建设地区年平均光照强度，单位为勒克斯，相关数据可以通过查询当地历史天气数据 分析得到;FsupplyWT(x)表示HRES配置方案为X下的风力机发电量，V为年平均风速，单位为米/秒，相关数据可以通过查询当地历史天气数据分析得到;FsupplyBA (X)表示HRES配置方案为X下的备选储能设备供电量，FsupplyFC (X)表示HRES配置方案为X下的柴油设备发电量，Pm表示第i类型光伏板的额定功率;Pm表示第j类型风力机的额定功率，PBAk表示第k 类型备选储能设备的额定功率，Pfw表示第1类型柴油机的额定功率； Fdemand表为： Fd_nd = Q，其中Q为某地区年平均用电量为Q(瓦）。3. 根据权利要求1所述的一种孤岛混合可再生能源系统的多目标组合优化配置方法， 其特征在于所述的智能多目标优化算法NSGA-II的具体过程为： 步骤2.1:设置算法参数:包括种群规模和终止条件，种群规模为N，取值范围100至500， 终止条件采用最大运行代数maxgen，取值范围50至500; 步骤2.2:初始化种群:随机生成种群规模N= 100父代种群S;种群中每个个体X包含η个 编码值，其中n = Npv+Nwt+Nba+2*Nf g，令当前运行代数gen = 1; 步骤2.3:基于当前种群S，通过交叉变异操作产生子代种群Sc，种群规模也为N; 步骤2.3.1:交叉操作； 步骤2.3.2:变异操作，设置变异操作的概率为σ，σ e (〇，1); 步骤2.4:将父代种群S与子代种群Sc合并，得到规模为2N的合种群5311 = S U Sc，对合种 群Sall中的个体进行非劣分层，然后计算每一个非劣层的个体局部拥挤距离，最后依据个体 所处的非劣层及拥挤距离，对所有个体进行排序； 步骤2.4.1:对种群Sall中的个体进行非劣分层； 步骤2.4.2:计算种群中个体的拥挤距离，拥挤距离越小，个体周围越稠密； 步骤2.4.3:基于个体所处的非劣层及拥挤距离，对所有个体进行排序； 步骤2.5:在排序之后的个体种群中选取第1到N个个体作为新的父代种群S，令gen = gen+1 ； 步骤2.6:判断gen是否大于maxgen，若大于，则输出种群S中的处于第1非劣层中的个体 作为步骤1中所确定的多目标优化模型的Pareto最优解集，转步骤3;否则转步骤2.3。4. 根据权利要求3所述的一种孤岛混合可再生能源系统的多目标组合优化配置方法， 其特征在于所述步骤2.1中的种群规模取值为100，最大运行代数取值为100。5. 根据权利要求3所述的一种孤岛混合可再生能源系统的多目标组合优化配置方法， 其特征在于所述交叉操作的具体方法为： 步骤2.3.1.1:令q= I，并初始化集合Sc '为空集； 步骤2.3.1.2:选取个体12。-1和12。，12。-1已3，12。已3，，首先随机生成一个位于区间[1，11] 之间的整数，作为交叉位置，记为i'，然后用X2trl的前半段第1到第i'个编码和x2q的后半段 即第i'+l到最后一个编码组成子个体的后半段即第i'+l到最后一个编码和 X2q 的前半段即第1到第i'个编码组成子个体x'2q;将产生的个体^2<1- 1和1'2(1存入集合3(3'，即 Scj=Sc' Ux' 2q-i Ux^q; 步骤2.3.1.3:令q = q+1;若q>N/2，转向步骤2.3.2，否则转向步骤2.3.1.2。6. 根据权利要求3所述的一种孤岛混合可再生能源系统的多目标组合优化配置方法， 其特征在于所述变异操作的具体方法为： 步骤2.3.2.1:令q= 1，初始化Sc为空集； 步骤2.3.2.2:针对个体Xq e Sc '，随机生成一个位于（O，1)区间的实数，记为rand;若 rand〈〇,则随机选取一个变异位置v，ve [I，n]，若ve [I，n-Nfg]，则随机选取[I，Nl]区间的 一个正整数代替个体Xq的第V个变量值，得到变异个体X'q,并存入集合Sc，即Sc = Sc U X'q; 若ve [n-Nfg+l，n]，则随机选取位于(0，C]区间的实数代替个体Xq的第V个变量值，得到变异 个体X ' q，并存入集合Sc，即Sc = Sc U X ' q;若rand 2 〇，将Xq存入集合Sc，即Sc = Sc U Xq; 步骤2.3.2.3:令q = q+l;若q>N，转向步骤2.4,否则转向步骤2.3.2.2。7. 根据权利要求3所述的一种孤岛混合可再生能源系统的多目标组合优化配置方法， 其特征在于所述步骤2.3.2中〇 = 〇.01。8. 根据权利要求3所述的一种孤岛混合可再生能源系统的多目标组合优化配置方法， 其特征在于非劣分层方法为： 步骤 2.4.1.1:令11=1; 步骤2.4.1.2:找出种群5311中不被任何个体Pareto约束支配的个体，并保存在集合Ah 中，即集合Ah中的个体均属于第h非劣层； 所述个体X Pareto约束支配个体y，是指当下列其中一个条件满足时：（i)个体X和y均 满足约束条件且1<7; (ii)个体X满足约束条件，而y不满足约束条件；（iii)个体X和y均不 满足约束条件，个体X违反约束条件的程度小于个体y违反约束条件的程度，则称个体X Pareto约束支配个体y;所述的关系式1<又表示个体X Pareto支配个体y，是在所有目标函 数值上，个体X的目标函数值都不大于个体y的目标函数值即V/；?，厶㈨s/Jy),并且至少 存在一个目标函数上个体X的目标函数值小于个体y的目标函数值，即3./；.,fm，（x)〈f m，（y)， 所述fm，m=l，2为步骤1的公式（1)中的两个目标函数，表示目标函数Frast(X)的值，f 2表示 目标函数^：。:(幻的值； 步骤2.4.1.3: WSaIi中除去所有在集合Ah中的个体，判断当前Saii是否为空集，若是，则 分层完毕，转入步骤2.4.2;否则转入步骤2.4.1.2，令h = h+l。9. 根据权利要求3所述的一种孤岛混合可再生能源系统的多目标组合优化配置方法， 其特征在于所述拥挤距离计算方法为： 步骤2.4.2.1:计算种群Sall内每个个体的目标函数值，并对其进行排序;所述目标函数 值指在步骤1中的目标函数Frast(X)和/(X)的值； 步骤2.4.2.2:计算拥挤距离;首先定义在种群中目标函数ASf2最小的个体为边界个 体，边界个体的拥挤距离dist定义为无穷大;其次计算除边界外的个体XJi挤距离dist(q) 为其中/Tx和分别表示当前种群中目标函数的最大和最小值； /Γ1和/?Γ1表示个体Xq-I和个体Xq+ι的目标函数值。10.根据权利要求3所述的一种孤岛混合可再生能源系统的多目标组合优化配置方法， 其特征在于所述排序准则为:处于第h非劣层的个体排序在第h '非劣层的个体之前，其中h ' >h;若两个个体处于同一非劣层，那么拥挤距离大的个体排序在拥挤距离小的个体之前。
【专利摘要】本发明公开了一种孤岛混合可再生能源系统的多目标组合优化配置方法，解决的技术问题是如何合理配置风力机、光伏板、柴油机以及储能设备以使得整个孤岛混合可再生能源系统满足区域用能需求，且经济成本最低、环境危害最小。采用的技术方案是建立HRES的组合优化配置模型，利用智能多目标优化算法求解，得到一组分布均匀的Pareto最优解集，从这组最优解集中，选择一个作为最终配置方案来进行配置。取得的有益效果一是针对HRES中设备的组合配置构建了带约束的多目标组合优化模型，更符合实际情况、方案可行性更强；二是采用多目标智能优化算法求解，能够同时找到一组Pareto最优解集，满足不同情况要求下的HRES方案配置。
【IPC分类】G06N3/00, G06Q50/06, G06Q10/04
【公开号】CN105488584
【申请号】CN201510816203
【发明人】王锐, 雷洪涛, 张涛, 查亚兵, 刘亚杰, 史志超
【申请人】中国人民解放军国防科学技术大学
【公开日】2016年4月13日
【申请日】2015年11月20日
【技术领域】
[0001] 本发明涉及一种孤岛混合可再生能源系统的多目标组合优化配置方法，具体而言 是考虑能源系统中的各单元的类型和数量的配置，最优化系统的各个目标，属于计算机领 域。
【背景技术】
[0002] 能源问题始终伴随着人类文明发展的整个进程，特别是进入21世纪后，工业经济 的急速发展对能源需求提出了更高的要求。然而煤、石油之类的化石能源日益枯竭，并且使 用化石燃料带来的环境问题日益凸显，能源安全和环境污染已经成为当前世界各国面临的 巨大挑战。为了应对能源和环境的双重危机，大力发展和利用可再生能源，如太阳能、风能 等，成为当务之急。
[0003]混合可再生能源系统（Hybrid Renewable Energy Systems,HRES)顾名思义，是一 种混合使用化石能源和可再生能源的系统。具体而言是包括风力发电、光伏发电、柴油发电 以及储能装置的能源系统。孤岛混合可再生能源系统指HRES没有接入系统电网。孤岛混合 可再生能源系统通常用于海岛之类的偏远地区。用户相对较少，用电量需求较低，引入大电 网比较困难且经济成本较高。
[0004] 孤岛混合可再生能源系统规划需要解决的问题是：如何合理配置风力机、光伏板、 柴油机以及储能设备以使得整个能源系统满足区域用能需求，同时经济成本最低、对环境 危害最小。
[0005] 目前关于孤岛混合可再生能源系统的组合优化配置研究主要集中于单目标的优 化模型，如在满足供需平衡下，最小化HRES全寿命周期成本或最小化HRES的温室气体排放 量等。然而，从实际角度出发，HRES的组合优化配置需要同时考虑全寿命周期成本、温室气 体排放量等评估指标，是一个多目标优化问题。同时，HRES的组合优化配置中涉及到风光柴 储设备的选择，柴油使用量等变量的规划，是一个多类型变量的（连续、离散变量)优化问 题。此外，HRES组合优化配置还需满足供需平衡约束。综上所述，HRES的组合优化配置是一 个多变量、多目标、带约束的复杂优化问题，目前关于HRES优化配置的相关研究(包括模型 及求解算法)远不能满足实际需求。
[0006] 所述的多目标优化问题是指：同时对多个目标进行优化，由于各个目标之间通常 是耦合在一起且互相制约、互相竞争，即某个目标的改善可能引起其他目标性能的降低，很 难找到一个真正意义上的最优解使得各个目标同时达到最优，因此多目标优化问题的最优 解通常不是单一的，而是一组多个互有利弊的非支配解，是一个非劣解的集合，即帕累托 (Pareto)最优解集。求解多目标优化问题的核心是找到一组分布均勾的Pareto最优解。
[0007] 所述的智能多目标优化算法是指:传统的处理多目标优化问题的方法，如加权法、 约束法、目标规划法等，通过构建一个评价函数，将多目标问题转化为单目标优化问题，然 后利用一般的求解方法计算得到问题的一个解。智能多目标优化算法是通过模拟某些自然 过程发展而来的基于种群的优化算法，其思想和内容涉及数学、生物学和计算机学科等。该 类算法不依赖于梯度信息，一次运行能够找到一组Pareto最优解，具有全局、并行、高效、鲁 棒和通用性强等特点。是求解复杂非线性多目标优化问题的有效方法。
【发明内容】
[0008] 本发明要解决的技术问题是如何合理配置风力机、光伏板、柴油机以及储能设备 以使得整个孤岛混合可再生能源系统满足区域用能需求，且经济成本最低、环境危害最小。
[0009] 为解决该问题，本发明所采取的技术方案是:建立HRES的组合优化配置模型，利用 智能多目标优化算法求解，得到一组分布均勾的Pareto最优解，从这组Pareto最优解中，选 择一个作为最终配置方案来对HRES进行配置。
[0010]步骤1:以HRES全寿命周期成本最低，HRES对环境危害最小为目标，以满足区域能 源需求为约束条件，建立HRES的多目标组合优化配置模型：
[0011]
[0012]
[0013]
[0014] 式（1)表示以最小化HRES的全寿命周期成本和系统的C〇2排放量为目标，以系统供 能满足用户用电需求为约束条件，构造出带约束的两目标优化模型；
[0015] 其中Fccist(x)表示HRES配置方案为X下的全寿命周期（T年）成本，包括孤岛风光柴 储能源系统中各类设备的初始安装费用和使用过程中的维护费用；
[0016] &〇2(Χ)表示HRES配置方案为X下系统的 C〇2排放量，包括HRES中风光柴储设备在生 产制造过程的c〇2排放量和柴油机使用过程中的c〇2排放量，其中风力机、光伏板以及储能设 备在使用过程中的c〇2排放量不予考虑；
[0017] Fsupply(x)表示HRES配置方案为X下的供电量；
[0018] Fdemand表示某一地区的用电需求量；
[0019] 所述配置方案X是指HRES系统中光伏组件、风力机、储能设备、柴油机四类组件的 配置数量以及柴油使用量，配置方案X表示如下：
[0020]
[0021] xi表示HRES中第i类型的光伏板个数，ie{l，2，. . .，NPV}; Vi表示HRES中第j类型 的风力机个数，je{l，2, . . .，Nwt};xvv^表示HRES中第k类型的储能设备个数，ke{l， 2，. . .，Nba}; 表示HRES中第 1 类型的柴油机个数，1 e {1，2，. . .，Nfg}; 示HRES中第1类型的柴油机的柴油使用量，单位为升，取值范围为[0，C]，C表示柴油使用量 的上限值，对于柴油发电机，若未使用某一类型的柴油机，那么其对应的柴油使用量设定为 〇;配置方案X的前Npv+N wt+Nba+Nfg个变量取值范围为(0，Nl]，Nl表示设备个数的上限值;N PV表 示光伏板的类型数，Nwt表示风力机的类型数，他3表示备选储能设备的类型数，心 8表示柴油 机的类型数；
[0022] Fccist(x)表示为：
[0023]
[0024] Fc〇stPV(x)表不腿S配置方案为x下的光伏板全寿命周期费用，匕燃(》) = 1:'-; <?", + C2n， Fcos· (x)表不HRE S配置方案为x下的风力机全寿命周期费用，4tWT W = Σ ,1卜n?β+ ""r) * 心他⑴表示TOS配置方案为x下的备选储能设备全寿命周期费用，，; FmstreW^^^RES配置方案为x下的柴油机周期费用， Clpvi表不第i类型光伏板的初装费用，C2pvi表不第i类型光伏板的维护费用；ClwTj表不第j类 型风力机的初装费用，C2叩表示第j类型风力机的维护费用;ClBAk表示第k类型备选储能设 备的初装费用，C2 BAk表示第k类型备选储能设备的维护费用;C1FC1表示第1类型柴油机的初 装费用，C2?表示第1类型柴油机的维护费用;C FG表示柴油价格；
[0025]
[0026]
[0027] Fc〇2PV(x)ifcrii?Sl￡Kir案为 X 下的光伏 表不HRES配置方案为X下的风力机全寿命周期碳排放量，= ; .&α2ΒΑ(χ_.)_ 表示HRES配置方案为X下的备选储能设备全寿命周期碳排放量，&2βα> = ΣΓ:;W
表不HRES配置方案χ下的柴油设备全寿命周期碳排放量，+~％.??.*+?.，*<)* Pm表示第i类型光伏板的额定功率;〇PVl表示制造单位面积的第i类型光伏板产生的碳排放 量;Own表示制造第j类型风力机产生的碳排放量;〇BAk表示制造第k类型备选储能设备产生 的碳排放量;〇 FC1表示第1类型柴油发电设备产生的碳排放量;表示燃烧每升柴油的碳排 放量；
[0028]
[0029]
[0030] Fsuppiypv(x)表不HRES配置方案为X下的光伏板发电量，.FsupplyPV(x) = i：、； L为HRES建设地区年平均光照强度，单位为勒克斯，相关数据可以通过查询当地历史天气数 据分析得到；FsupplyWT(X)表不HRES配置方案为X下的风力机发电量， AUpplyWT(X) = * Γ * ; V为年平均风速，单位为米/秒，相关数据可以通过查询当 地历史天气数据分析得到；FsupplyBA(x)表示HRES配置方案为X下的备选储能设备供电量， 巧upplyBA(X) = Σ?? * ' J :; FsupplyFG ( X )表不HRES配置方案为X下的柴油设备发电直， W / V ) - \ ^ ，Γ ρ '1·· γ 承：y .? Siippl5Hrv ~ Lui-^X · . -十/ ?卡1+"
[0031] PPVl表示第i类型光伏板的额定功率;PWTj表示第j类型风力机的额定功率，P BAk表示 第k类型备选储能设备的额定功率，ρ?表示第1类型柴油机的额定功率；
[0032] Fdemand表不为：
[0033] Fd_nd = Q，其中Q为某地区年平均用电量为Q(瓦）；
[0034] 步骤2:利用智能多目标优化算法NSGA-II求解HRES组合优化模型，具体过程如下：
[0035] 步骤2.1:设置算法参数:包括种群规模和终止条件，种群规模为N，取值范围100至 500，终止条件采用最大运行代数maxgen，取值范围50至500;
[0036] 步骤2.2:初始化种群：随机生成种群规模N= 100父代种群S;种群中每个个体X包 含η个编码值，其中n = Npv+Nwt+Nba+2*Nf g，令当前运行代数gen = 1;
[0037] 步骤2.3:基于当前种群S，通过交叉变异操作产生子代种群Sc ，种群规模也为N。所 述交叉变异操作实施方法如下：
[0038]步骤2.3.1:交叉操作；
[0039]步骤2.3.1.1:令q = 1，并初始化集合Sc '为空集；
[0040] 步骤2.3.1.2:选取个体12。-1和12。^2。-1三3，12。三3，，首先随机生成一个位于区间
[l，n]之间的整数，作为交叉位置，记为i'，然后用X2H的前半段第1到第i'个编码和^。的后 半段即第i'+l到最后一个编码组成子个体的后半段即第i'+l到最后一个编码和 x2q的前半段即第1到第i'个编码组成子个体x'2q;将产生的个体^2< 1-1和1'2(1存入集合3(3'， 艮P Sc ' = Sc ' U X ' 2q-1 U X ' 2q ;
[0041 ] 步骤2.3.1.3:令9 = 9+1;若9>~/2，转向步骤2.3.2，否则转向步骤2.3.1.2;
[0042] 步骤2.3.2:变异操作，设置变异操作的概率为σ，σ e (〇，1);
[0043] 步骤2 · 3 · 2 · 1:令q = 1，初始化Sc为空集；
[0044] 步骤2 · 3 · 2 · 2 :针对个体xqe Sc '，随机生成一个位于(0，1)区间的实数，记为rand。 若rand〈〇,则随机选取一个变异位置v，ve [1，n]，若ve [1 ,η-Nfg]，则随机选取[1，Nl]区间 的一个正整数代替个体Xq的第V个变量值，得到变异个体X'q,并存入集合SC，即SC = SC U X 'q;若ve[n-Nfg+l，n]，则随机选取位于(0，C]区间的实数代替个体 Xq的第v个变量值，得到 变异个体X ' q，并存入集合Sc，即Sc = Sc U X ' q;若rand 2 〇，将Xq存入集合Sc，即Sc = Sc U Xq; [0045] 步骤2.3.2.3:令q = q+l;若q>N，转向步骤2.4,否则转向步骤2.3.2.2;
[0046]步骤2.4:将父代种群5与子代种群5(：合并，得到规模为2_勺合种群5311 = 51^(3，对 合种群Sall中的个体进行非劣分层，然后计算每一个非劣层的个体局部拥挤距离，最后依据 个体所处的非劣层及拥挤距离，对所有个体进行排序；
[0047] 步骤2.4.1:对种群Saii中的个体进行非劣分层。下面介绍非劣分层的方法；
[0048] 步骤 2.4.1.1:令h = l;
[0049] 步骤2.4.1.2:找出种群5311中不被任何个体Pareto约束支配的个体，并保存在集 合Ah中，即集合Ah中的个体均属于第h非劣层;所述个体x Pareto约束支配个体y，是指当下 列其中一个条件满足时：（i)个体X和y均满足约束条件且κ y; (i i)个体X满足约束条件，而 y不满足约束条件；（i ii)个体X和y均不满足约束条件，个体X违反约束条件的程度小于个体 y违反约束条件的程度，则称个体X Pareto约束支配个体y;所述的关系式x *<y表示个体X Pareto支配个体y，是在所有目标函数值上，个体X的目标函数值都不大于个体y的目标函数 值即VX，f m(xHfm(y)，并且至少存在一个目标函数上个体X的目标函数值小于个体y的目 标函数值，即3九.上令)〈匕，（7)，所述匕，111=1，2为步骤1的公式（1)中的两个目标函数4 1表 示目标函数F?st(x)的值，5表示目标函数巧。2(幻的值；
[0050] 步骤2.4.1.3 : WSaii中除去所有在集合Ah中的个体，判断当前Saii是否为空集，若 是，则分层完毕，转入步骤2.4.2;否则转入步骤2.4.1.2，令h = h+l;
[0051] 步骤2.4.2:计算种群中个体的拥挤距离，拥挤距离越小，个体周围越稠密;拥挤距 离具体计算方法如下：
[0052] 步骤2.4.2.1:计算种群Sall内每个个体的目标函数值，并对其进行排序;所述目标 函数值指在步骤1的公式(1)中的目标函数F? st(x)和匕):(x)的值；
[0053] 步骤2.4.2.2:计算拥挤距离;首先定义在种群中目标函数心或5最小的个体为边 界个体，边界个体的拥挤距离dist定义为无穷大；其次计算除边界外的个体^拥挤距离 dist(q)为
[0054]
[0055] 其中./Γχ和./Γ分别表示当前种群中目标函数匕的最大和最小值；/T1和表 示个体X(rl和个体Xq+1的目标函数fVf直，fm=1表示目标函数Fc^stU)的值，fm=2表示目标函数 &>)的值；
[0056]步骤2.4.3:基于个体所处的非劣层及拥挤距离，对所有个体进行排序，排序准则 为:处于第h非劣层的个体排序在第h '非劣层的个体之前，其中h ' >h;若两个个体处于同一 非劣层，那么拥挤距离大的个体排序在拥挤距离小的个体之前；
[0057]步骤2.5:在排序之后的个体种群中选取第1到N个个体作为新的父代种群S，令gen = gen+l ；
[0058]步骤2 · 6:判断gen是否大于maxgen，若大于，则输出种群S中的处于第1非劣层中的 个体作为步骤1中所确定的多目标优化模型的Pareto最优解集，转步骤3;否则转步骤2.3。 [0059]步骤3:从Pareto最优解集中选择最优HRES方案；
[0060] 步骤3.1当要求整个系统的全寿命周期成本小于等于给定值CostA时，选择Pareto 最优解集中目标函数值6:+00最小的解，即使系统C02排放量最小的个体X，作为最终的配 置方案X，所述CostA指HRES配置过程所要求的最大成本额度；
[0061 ]步骤3.2当要求整个系统的⑶2排放量小于等于给定值⑶2A时，选择Pareto最优解 集中目标函数值F_t(X)最小的解，即使系统全寿命周期成本最小的个体X，作为最终的配 置方案X，所述C0 2#BHRES配置过程所要求的最大排放量；
[0062]本发明的优点在于：
[0063] (1)本发明针对HRES中设备的组合配置构建了带约束的多目标组合优化模型，更 符合实际情况、方案可行性更强。
[0064] (2)本发明采用多目标智能优化算法求解问题，能够同时找到一组Pareto最优解 集，满足不同情况(即，对成本和环境污染的不同要求)下的HRES方案配置。
【附图说明】
[0065] 图1典型孤岛混合可再生能源系统示意图
[0066]图2是利用智能多目标优化算法NSGA-II求解HRES组合优化配置问题的流程图
[0067] 图3 HRES方案设计编码
[0068]图4单点交叉和均匀变异算子示意图
[0069]图5个体非劣分层及拥挤距离计算示意图
[0070]图6 HRES规划的Pareto最优解集
[0071] 图7 HRES规划的决策过程示意图
【具体实施方式】
[0072] 下面结合附图和实例对本发明作进一步的详细说明。
[0073] 本发明提供的孤岛混合可再生能源系统的多目标组合优化配置方法，具体步骤如 下：
[0074]表1至4给出了HRES系统（如图1所示）中各设备的相关数据，同时已知在HRES建设 地区年平均光照强度L = 20(单位:勒克斯）、年平均风速V = 5(单位:米/秒）以及年平均用电 量Q=1000(单位:千瓦）；
[0075] 表1备选光伏板相关参数(假设每块光伏板面积为1米2)
[0076]
[0077]表2备选风力机相关参数 [0078]
[0079] 表3备选储能设备相关参数
[0080]
[0081 ]表4备选柴油机相关参数 [0082]
[0083] 步骤1:考虑HRES的全寿命周期T = 20年，以HRES全寿命周期成本 [0084]最低，系统排放C02量最小为目标，以满足区域能源需求为约束条件，
[0085]建立HRES的多目标组合优化配置模型；
[0086] mLu(Fc〇st(x),FCOi(x)!
[0087] S.t.Fsupply(x) ^ Fdemand
[0088] 步骤2:利用智能多目标优化算法NSGA-II求解HRES组合优化模型，算法流程见图 2,具体过程描述如下：
[0089] 步骤2.1:设置算法参数:包括种群规模和终止条件，种群规模N设置为100，终止条 件采用最大运行代数maxgen，设置为maxgen = 100。
[0090] 步骤2.2:初始化种群:令当前运行代数8的=1。随机生成~=100父代种群5。种群 中每个个体X包含η个编码值，其中n = Npv+Nwt+Nba+2*Nfg，且N L取5，C取100;由表1 -4数据得 到，本实例中所考虑的光伏板类型数Npv = 5,风力机类型数Nwt = 5,储能设备类型数Nba = 3， 柴油发电机类型数Nfg = 5,个体编码长度n = 5+5+3+2*5 = 23,即每个个体包括23个编码值。 图3给出了针对此算例的一个个体编码。
[0091]步骤2.3:基于当前种群S，通过交叉变异操作（如图4所示)产生子代种群Sc，种群 规模也为N。所述交叉变异操作实施方法如下：
[0092]步骤2.3.1:交叉操作；
[0093]步骤2 · 3 · 1 · 1:令q = 1，并初始化集合Sc '为空集；
[0094] 步骤2.3.1.2:选取个体12。-1和12。^2。-1三3，12。三3，，首先随机生成一个位于区间 [l，n]之间的整数，作为交叉位置，记为i'，然后用X2 H的前半段第1到第i'个编码和^。的后 半段即第i'+l到最后一个编码组成子个体的后半段即第i'+l到最后一个编码和 x2q的前半段即第1到第i'个编码组成子个体x'2q;将产生的个体x'kl·和X' 2q#入集合Sc_l， 艮P Sc ' = Sc ' U X ' 2q-1 U X ' 2q ;
[0095] 步骤2.3.1.3:令9 = 9+1;若9>~/2，转向步骤2.3.2，否则转向步骤2.3.1.2;
[0096] 步骤2.3.2:变异操作，设置变异操作的概率为σ，〇6(0，1)，取〇 = 0.01;
[0097] 步骤2 · 3 · 2 · 1:令q = 1，初始化Sc为空集；
[0098] 步骤2 · 3 · 2 · 2 :针对个体xqe Sc '，随机生成一个位于(0，1)区间的实数，记为rand。 若rand〈〇,则随机选取一个变异位置v，ve [1，n]，若ve [1 ,η-Nfg]，则随机选取[1，Nl]区间 的一个正整数代替个体Xq的第V个变量值，得到变异个体X'q,并存入集合SC，即SC = SC U X 'q;若ve[n-Nfg+l，n]，则随机选取位于(0，C]区间的实数代替个体 Xq的第v个变量值，得到 变异个体X ' q，并存入集合Sc，即Sc = Sc U X ' q;若rand 2 〇，将Xq存入集合Sc，即Sc = Sc U Xq;
[0099] 步骤2 · 3 · 2 · 3:令q = q+1;若q>N，转向步骤2 · 4，否则
[0100] 转向步骤2.3.2.2;
[0101]步骤2.4:将父代种群3与子代种群5(：合并，得到规模为2_勺合种群3311 = 31^(3，对 合种群Sall中的个体进行非劣分层，然后计算每一个非劣层的个体局部拥挤距离，最后依据 个体所处的非劣层及拥挤距离，对所有个体进行排序；
[0102] 步骤2.4.1:对种群Sall中的个体进行非劣分层。下面介绍非劣分层的方法；
[0103] 步骤 2.4.1.1:令h = l;
[0104] 步骤2.4.1.2:找出种群5311中不被任何个体Pareto约束支配的个体，并保存在集 合Ah中，即集合Ah中的个体均属于第h非劣层;所述个体X Pareto约束支配个体y，是指当下 列其中一个条件满足时：（i)个体X和y均满足约束条件且个体X满足约束条件，而 y不满足约束条件；（i ii)个体X和y均不满足约束条件，个体X违反约束条件的程度小于个体 y违反约束条件的程度，则称个体X Pareto约束支配个体y;所述的关系式χ ?表示个体X Pareto支配个体y，是在所有目标函数值上，个体X的目标函数值都不大于个体y的目标函数 值即f m(y)，并且至少存在一个目标函数上个体X的目标函数值小于个体y的目 标函数值，即3人^七)仏，（7)，所述匕，111=1，2为步骤1的公式(1)中的两个目标函数4 1表 示目标函数F?st(x)的值，f2表示目标函数的值；
[0105] 步骤2.4.1.3 : Wsaii中除去所有在集合Ah中的个体，判断当前Sail是否为空集，若 是，则分层完毕，转入步骤2.4.2;否则转入步骤2.4.1.2，令h = h+l;
[0106] 步骤2.4.2:计算种群中个体的拥挤距离，拥挤距离越小，个体周围越稠密;拥挤距 离具体计算方法如下：
[0107] 步骤2.4.2.1:计算种群Sall内每个个体的目标函数值，并对其进行排序;所述目标 函数值指在步骤1中的目标函数F? st(x)和/*& (>0的值；
[0108] 步骤2.4.2.2:计算拥挤距离;首先定义在种群中目标函数心或5最小的个体为边 界个体，边界个体的拥挤距离dist定义为无穷大；其次计算除边界外的个体^拥挤距离 dist(q)为
[0109]
[0110] 六r T Jrn yj =| HU TT群中目标函数fm的最大和最小值；尤丨__和/0丨：表 示个体X(rl和个体Xq+1的目标函数fVf直，fm=1表示目标函数Fc^stU)的值，fm=2表示目标函数 匕K(X)的值；
[0111] 步骤2.4.3:基于个体所处的非劣层及拥挤距离，对所有个体进行排序，排序准则 为:处于第h非劣层的个体排序在第h '非劣层的个体之前，其中h ' >h;若两个个体处于同一 非劣层，那么拥挤距离大的个体排序在拥挤距离小的个体之前；
[0112] 步骤2.5:在排序之后的个体种群中选取第1到N个个体作为新的父代种群S，令gen = gen+l ；
[0113]步骤2 · 6:判断gen是否大于maxgen，若大于，则输出种群S中的处于Αι集合（即第1 非劣层中）的个体作为步骤1中所确定的多目标优化模型的Pareto最优解集，转步骤3;否则 转步骤2.3。
[0114]步骤3:从Pareto最优解集中选择最优HRES方案：依据本实例中建立的模型，利用 智能多目标优化算法NSGA-II求解模型，得到如图6所示的一组Pareto最优解集。
[0115]步骤3.1:当整个系统的全寿命周期成本小于等于给定值CostA时，选择Pareto最 优解集中目标函数值^:(x)最小的解，即使系统C02排放量最小的个体X，作为最终的配置 方案X，如图7所示，所述CostA指HRES配置过程所要求的最大成本额度。当要求HRES系统的 全寿命周期成本小于5.26*104千元，即CostA=5.26*104，HRES系统配置为使用36个#1类型 光伏板，19个#2类型光伏板，不使用#3、#4和#5类型的光伏板;使用41个#1类型风力机，22 个#3类型风力机，14个#4类型风力机，不使用#2、#5类型的风力机;使用26个#1类型储能设 备，不使用#2和#3类型的储能设备;使用4台#2类型柴油机，每台使用柴油183升，不使用其 他类型的柴油机；
[0116] 步骤3.2:当整个系统的C02排放量小于等于给定值C02A时，选择Pareto最优解集中 目标函数值F_ t(x)最小的解，即使系统全寿命周期成本最小的个体X，作为最终的配置方 案X，如图7所示。所述C02#gHRES配置过程所要求的最大排放量。当要求HRES系统的C0 2排放 量小于0.9*104克，即C02A=0.9*104，HRES系统配置为使用69个#1类型光伏板，不使用#2、# 3、#4和#5类型的光伏板;使用72个#1类型风力机，31个#2类型风力机，不使用#3、#4、#5类型 的风力机;使用52个#1类型储能设备，12个#2类型储能设备，不使用#3类型的储能设备;使 用1台#1类型柴油机，使用柴油11升，不使用其他类型的柴油机；
[0117] 本发明不仅限于上述【具体实施方式】，本领域一般技术人员根据本发明公开的内 容，可以采用其它多种具体方式实施本发明，如采用其他类似的多目标智能优化算法（如 MOGA、MOEA/D)。凡是采用本发明的技术思路，做一些简单的变化或改善的设计，如修改算法 参数等，都属于本发明保护的范围。
【主权项】
1. 一种孤岛混合可再生能源系统的多目标组合优化配置方法，其特征在于所述方法包 括以下具体步骤： 步骤1:建立HRES的多目标组合优化配置模型； 以HRES全寿命周期成本最低，HRES对环境危害最小为目标，以满足区域能源需求为约 束条件； 步骤2:利用智能多目标优化算法NSGA-II求解HRES组合优化模型，输出Pareto最优解 集； 步骤3:从Pareto最优解集中选择最优HRES方案； 步骤3.1当要求整个系统的全寿命周期成本小于等于给定值CostA时，选择Pareto最优 解集中目标函数值己,Jx)最小的解，即使系统CO2排放量最小的个体X，作为最终配置方案， 所述CostA指HRES配置过程所要求的最大成本额度； 步骤3.2当要求整个系统的C〇2排放量小于等于给定值C〇2A时，选择Pareto最优解集中 目标函数值Frast(X)最小的解，即使系统全寿命周期成本最小的个体X，作为最终配置方案， 所述CO 2WBHREs配置过程所要求的最大排放量。2. 根据权利要求1所述的一种孤岛混合可再生能源系统的多目标组合优化配置方法， 其特征在于所述的组合优化配置模型为：式（1)表示以最小化HRES的全寿命周期成本和系统的CO2排放量为目标，以系统供能满 足用户用电需求为约束条件，构造出带约束的两目标优化模型； 其中Frast(X)表示HRES配置方案为X下全寿命周期T年的成本，包括孤岛风光柴储能源系 统中各类设备的初始安装费用和使用过程中的维护费用； 々·〇_.(X)表示HRES配置方案为X下系统的CO2排放量，包括HRES中风光柴储设备在生产制 造过程的CO2排放量和柴油机使用过程中的CO2排放量，其中风力机、光伏板以及储能设备在 使用过程中的CO2排放量不予考虑； Fsupply(X)表示HRES配置方案为X下的供电量； Fd_nd表示某一地区的用电需求量； 所述配置方案X是指HRES系统中光伏组件、风力机、储能设备、柴油机四类组件的配置 数量以及柴油使用量，配置方案X表示如下：Xi表示HRES中第i类型的光伏板个数，ie{l，2，. ..，NPV}; V,表示HRES中第j类型的风 力机个数，je{l，2,…表示HRES中第k类型的储能设备个数，ke{l，2,…， Nba}; ' · w?表示HRES中第1类型的柴油机个数，I e {1，2，. . .，Nfg}; 表示HRES 中第1类型的柴油机的柴油使用量，单位为升，取值范围为[〇，C]，C表示柴油使用量的上限 值，对于柴油发电机，若未使用某一类型的柴油机，那么其对应的柴油使用量设定为0;配置 方案X的前Npv+Nwt+N ba+Nfg个变量取值范围为(0，Nl]，Nl表示设备个数的上限值;N pv表示光伏 板的类型数，Nwt表示风力机的类型数，他3表示备选储能设备的类型数，阶 8表示柴油机的类 型数； Fcost ( X )表为：Fcdsir<x)表不HRES配置方案为X下的光伏板全寿命周期费用，FccsfflT(X)表示HRES配置方案为X下的风力机全寿命周期费用，FcosM(X)表示HRES配置方案为X下的备选储能设备全寿命周期费用，-00????????%下的柴槲/1^#^周期费甩Clm表示第i类型光伏板的初装费用，C2PVl表示第i类型光伏板的维护费用;Ch^表示第j类 型风力机的初装费用，C2叩表示第j类型风力机的维护费用;ClBAk表示第k类型备选储能设 备的初装费用，C2 BAk表示第k类型备选储能设备的维护费用;CIfc1表示第1类型柴油机的初 装费用，C2ra表示第1类型柴油机的维护费用;C fg表示柴油价格； Go2 表不为表不HRES配置方案为X下的光伏板全寿命周期碳排放量，Ao2WT(X)表不HRES配置方案为X下的风力机全寿命周期碳排放量，Α〇2βα(χ)表示HRES配置方案为X下的备选储能设备全寿命周期碳排放量，砰㈤表示HRES配置方案X下的柴油设备全寿命周期碳排放量，Pm表示第i类型光伏板的额定功率;Om表示制造单位面积的第i类型光伏板产生的碳 排放量;0叩表示制造第j类型风力机产生的碳排放量;〇BAk表示制造第k类型备选储能设备 产生的碳排放量;Ofc1表示第1类型柴油发电设备产生的碳排放量;B wcl表示燃烧每升柴油的 碳排放量； Fsu卯Iy(X)表不为：Fsuppiypv(x)表不HRES配置方案为X下的光伏板发电量，；:L为 HRES建设地区年平均光照强度，单位为勒克斯，相关数据可以通过查询当地历史天气数据 分析得到;FsupplyWT(x)表示HRES配置方案为X下的风力机发电量，V为年平均风速，单位为米/秒，相关数据可以通过查询当地历史天气数据分析得到;FsupplyBA (X)表示HRES配置方案为X下的备选储能设备供电量，FsupplyFC (X)表示HRES配置方案为X下的柴油设备发电量，Pm表示第i类型光伏板的额定功率;Pm表示第j类型风力机的额定功率，PBAk表示第k 类型备选储能设备的额定功率，Pfw表示第1类型柴油机的额定功率； Fdemand表为： Fd_nd = Q，其中Q为某地区年平均用电量为Q(瓦）。3. 根据权利要求1所述的一种孤岛混合可再生能源系统的多目标组合优化配置方法， 其特征在于所述的智能多目标优化算法NSGA-II的具体过程为： 步骤2.1:设置算法参数:包括种群规模和终止条件，种群规模为N，取值范围100至500， 终止条件采用最大运行代数maxgen，取值范围50至500; 步骤2.2:初始化种群:随机生成种群规模N= 100父代种群S;种群中每个个体X包含η个 编码值，其中n = Npv+Nwt+Nba+2*Nf g，令当前运行代数gen = 1; 步骤2.3:基于当前种群S，通过交叉变异操作产生子代种群Sc，种群规模也为N; 步骤2.3.1:交叉操作； 步骤2.3.2:变异操作，设置变异操作的概率为σ，σ e (〇，1); 步骤2.4:将父代种群S与子代种群Sc合并，得到规模为2N的合种群5311 = S U Sc，对合种 群Sall中的个体进行非劣分层，然后计算每一个非劣层的个体局部拥挤距离，最后依据个体 所处的非劣层及拥挤距离，对所有个体进行排序； 步骤2.4.1:对种群Sall中的个体进行非劣分层； 步骤2.4.2:计算种群中个体的拥挤距离，拥挤距离越小，个体周围越稠密； 步骤2.4.3:基于个体所处的非劣层及拥挤距离，对所有个体进行排序； 步骤2.5:在排序之后的个体种群中选取第1到N个个体作为新的父代种群S，令gen = gen+1 ； 步骤2.6:判断gen是否大于maxgen，若大于，则输出种群S中的处于第1非劣层中的个体 作为步骤1中所确定的多目标优化模型的Pareto最优解集，转步骤3;否则转步骤2.3。4. 根据权利要求3所述的一种孤岛混合可再生能源系统的多目标组合优化配置方法， 其特征在于所述步骤2.1中的种群规模取值为100，最大运行代数取值为100。5. 根据权利要求3所述的一种孤岛混合可再生能源系统的多目标组合优化配置方法， 其特征在于所述交叉操作的具体方法为： 步骤2.3.1.1:令q= I，并初始化集合Sc '为空集； 步骤2.3.1.2:选取个体12。-1和12。，12。-1已3，12。已3，，首先随机生成一个位于区间[1，11] 之间的整数，作为交叉位置，记为i'，然后用X2trl的前半段第1到第i'个编码和x2q的后半段 即第i'+l到最后一个编码组成子个体的后半段即第i'+l到最后一个编码和 X2q 的前半段即第1到第i'个编码组成子个体x'2q;将产生的个体^2<1- 1和1'2(1存入集合3(3'，即 Scj=Sc' Ux' 2q-i Ux^q; 步骤2.3.1.3:令q = q+1;若q>N/2，转向步骤2.3.2，否则转向步骤2.3.1.2。6. 根据权利要求3所述的一种孤岛混合可再生能源系统的多目标组合优化配置方法， 其特征在于所述变异操作的具体方法为： 步骤2.3.2.1:令q= 1，初始化Sc为空集； 步骤2.3.2.2:针对个体Xq e Sc '，随机生成一个位于（O，1)区间的实数，记为rand;若 rand〈〇,则随机选取一个变异位置v，ve [I，n]，若ve [I，n-Nfg]，则随机选取[I，Nl]区间的 一个正整数代替个体Xq的第V个变量值，得到变异个体X'q,并存入集合Sc，即Sc = Sc U X'q; 若ve [n-Nfg+l，n]，则随机选取位于(0，C]区间的实数代替个体Xq的第V个变量值，得到变异 个体X ' q，并存入集合Sc，即Sc = Sc U X ' q;若rand 2 〇，将Xq存入集合Sc，即Sc = Sc U Xq; 步骤2.3.2.3:令q = q+l;若q>N，转向步骤2.4,否则转向步骤2.3.2.2。7. 根据权利要求3所述的一种孤岛混合可再生能源系统的多目标组合优化配置方法， 其特征在于所述步骤2.3.2中〇 = 〇.01。8. 根据权利要求3所述的一种孤岛混合可再生能源系统的多目标组合优化配置方法， 其特征在于非劣分层方法为： 步骤 2.4.1.1:令11=1; 步骤2.4.1.2:找出种群5311中不被任何个体Pareto约束支配的个体，并保存在集合Ah 中，即集合Ah中的个体均属于第h非劣层； 所述个体X Pareto约束支配个体y，是指当下列其中一个条件满足时：（i)个体X和y均 满足约束条件且1<7; (ii)个体X满足约束条件，而y不满足约束条件；（iii)个体X和y均不 满足约束条件，个体X违反约束条件的程度小于个体y违反约束条件的程度，则称个体X Pareto约束支配个体y;所述的关系式1<又表示个体X Pareto支配个体y，是在所有目标函 数值上，个体X的目标函数值都不大于个体y的目标函数值即V/；?，厶㈨s/Jy),并且至少 存在一个目标函数上个体X的目标函数值小于个体y的目标函数值，即3./；.,fm，（x)〈f m，（y)， 所述fm，m=l，2为步骤1的公式（1)中的两个目标函数，表示目标函数Frast(X)的值，f 2表示 目标函数^：。:(幻的值； 步骤2.4.1.3: WSaIi中除去所有在集合Ah中的个体，判断当前Saii是否为空集，若是，则 分层完毕，转入步骤2.4.2;否则转入步骤2.4.1.2，令h = h+l。9. 根据权利要求3所述的一种孤岛混合可再生能源系统的多目标组合优化配置方法， 其特征在于所述拥挤距离计算方法为： 步骤2.4.2.1:计算种群Sall内每个个体的目标函数值，并对其进行排序;所述目标函数 值指在步骤1中的目标函数Frast(X)和/(X)的值； 步骤2.4.2.2:计算拥挤距离;首先定义在种群中目标函数ASf2最小的个体为边界个 体，边界个体的拥挤距离dist定义为无穷大;其次计算除边界外的个体XJi挤距离dist(q) 为其中/Tx和分别表示当前种群中目标函数的最大和最小值； /Γ1和/?Γ1表示个体Xq-I和个体Xq+ι的目标函数值。10.根据权利要求3所述的一种孤岛混合可再生能源系统的多目标组合优化配置方法， 其特征在于所述排序准则为:处于第h非劣层的个体排序在第h '非劣层的个体之前，其中h ' >h;若两个个体处于同一非劣层，那么拥挤距离大的个体排序在拥挤距离小的个体之前。
【专利摘要】本发明公开了一种孤岛混合可再生能源系统的多目标组合优化配置方法，解决的技术问题是如何合理配置风力机、光伏板、柴油机以及储能设备以使得整个孤岛混合可再生能源系统满足区域用能需求，且经济成本最低、环境危害最小。采用的技术方案是建立HRES的组合优化配置模型，利用智能多目标优化算法求解，得到一组分布均匀的Pareto最优解集，从这组最优解集中，选择一个作为最终配置方案来进行配置。取得的有益效果一是针对HRES中设备的组合配置构建了带约束的多目标组合优化模型，更符合实际情况、方案可行性更强；二是采用多目标智能优化算法求解，能够同时找到一组Pareto最优解集，满足不同情况要求下的HRES方案配置。
【IPC分类】G06N3/00, G06Q50/06, G06Q10/04
【公开号】CN105488584
【申请号】CN201510816203
【发明人】王锐, 雷洪涛, 张涛, 查亚兵, 刘亚杰, 史志超
【申请人】中国人民解放军国防科学技术大学
【公开日】2016年4月13日
【申请日】2015年11月20日
最新回复(0)