CN103854070A - 独立微网系统的随机优化规划方法 - Google Patents

独立微网系统的随机优化规划方法 Download PDF

Info

Publication number
CN103854070A
CN103854070A CN201410083264.8A CN201410083264A CN103854070A CN 103854070 A CN103854070 A CN 103854070A CN 201410083264 A CN201410083264 A CN 201410083264A CN 103854070 A CN103854070 A CN 103854070A
Authority
CN
China
Prior art keywords
capital
electricity price
independent micro
total
grid system
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201410083264.8A
Other languages
English (en)
Inventor
郭力
王成山
余舟子
赵波
张雪松
周丹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
Electric Power Research Institute of State Grid Zhejiang Electric Power Co Ltd
Original Assignee
Tianjin University
Electric Power Research Institute of State Grid Zhejiang Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tianjin University, Electric Power Research Institute of State Grid Zhejiang Electric Power Co Ltd filed Critical Tianjin University
Priority to CN201410083264.8A priority Critical patent/CN103854070A/zh
Publication of CN103854070A publication Critical patent/CN103854070A/zh
Pending legal-status Critical Current

Links

Images

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种独立微网系统的随机优化规划方法,提出了计及不确定因素的多目标随机机会约束规划模型,随机数据采用马尔科夫状态转移矩阵实现模拟;控制策略采用准稳态仿真策略进行模拟;优化指标为全寿命周期内的资本金现金流收益净现值和污染物排放量,所述全寿命周期内的资本金现金流收益净现值为所述独立微网系统在末年的资本金累计现金流净现值;所述污染物排放量为全寿命周期内各类污染物的排放总量;选择负荷容量缺失率作为概率约束条件;选择工程动态经济分析方法进行经济评估,提供了还本付息参考电价和期望收益参考电价。本发明可实现独立微网系统设备类型与设备容量的组合优化。

Description

独立微网系统的随机优化规划方法
技术领域
本发明属于微网优化规划设计领域,提出了包含微网全寿命周期内资本金现金流收益净现值和污染排放水平在内的风光柴储微网系统的多目标随机机会约束优化规划方法。
背景技术
随着能源需求的不断增长,温室效应愈加明显,风机、光伏等分布式发电单元,在微网经济性、环境性等方面起到的作用更加重要。独立微网系统可将多种类型的分布式发电单元组合在一起,有效发挥单一能源系统的优点,实现多种能源互补,提高整个独立微网系统的效率、能源利用率和供能可靠性。在偏远和海岛地区,相比线路搭建和线路损耗,独立微网系统具有更低的成本,是解决其供电的有效手段之一。
开展独立微网系统的优化规划设计,需要充分考虑独立微网系统内分布式电源的组合方案和运行控制策略。在确定性优化规划建模上,国内外的研究多采用基于气象数据和负荷数据的准稳态逐时仿真优化设计方法,根据已获得系统全寿命周期内的风速、光照强度、温度和负荷的延时数据(通常以小时为单位),利用准稳态仿真程序,计算不同类型可再生能源组合方案下系统的各项指标。
然而,采用确定性的规划方法,需要获取全寿命周期内小时级的气象和负荷数据。在实际规划设计阶段获取相关数据存在着一定的困难。同时,大量的不确定性因素,如风速、光照强度,负荷水平,电价、燃料价格等,将对独立微网系统的优化规划结果产生影响。在独立微网系统中确定性规划中引入随机变量,构建微网的随机优化规划模型,能够更加符合实际情况,从而使决策更加合理。
现阶段,关于计及不确定性条件下的随机优化规划研究主要将其转化为确定性等价类,采用确定性模型和算法对其优化求解。在优化规划设计上,现有规划设计方法中并未考虑系统内常规电源,如柴油发电的单机容量和设备台数的组合优化、多机设备的组合开机方式、以及系统备用容量和设备选型等问题。在计及不确定性因素的随机规划过程中,现有方法一般采用蒙特卡罗方法或类似方法仿真模拟不确定性因素,并通过大规模的统计和实验获得优化结果,这种方法计算量大,比较费时,而且不能体现随机变量的时序性。另外,现有多以系统总成本作为经济目标,既没有考虑项目建设期和银行贷款利息等实际工程经济因素,也没有反映出投资者关心的项目现金流量情况。
发明内容
本发明目的在于提出一种适用于独立风光柴储微网系统的多目标随机机会约束优化规划方法,有效计及独立微网系统规划过程中的各种不确定因素,以实现独立微网系统随机优化规划设计运行更高的经济性、可靠性和资源利用率。本发明的技术方案如下:
一种适用于独立风光柴储微网系统的随机优化规划方法,其特征在于,包括:
获取所述风力发电机、所述光伏电池、所述柴油发电机和所述储能电池的设备参数;
基于风速、大气清洁指数、负荷随机数据的概率模型,采用马尔可夫状态转换概率矩阵模拟上述随机数据;
根据所述设备参数和所述随机数据,利用NSGA-II多目标遗传算法进行优化求解,在求解过程中采用随机机会约束规划,获得所述独立微网系统的优化指标,最终通过所述多目标遗传算法获得优化结果,并依据经济评价环节获取优化结果的具体经济性信息;
所述随机机会约束规划采用准稳态仿真策略进行模拟,以负荷容量缺失率为概率约束条件;所述负荷容量缺失率指将系统备用容量计入净负荷中,所得系统未满足的负荷容量与整个负荷需求容量的比值;
所述独立微网系统的优化指标为全寿命周期内的资本金现金流收益净现值和污染物排放量,所述全寿命周期内的资本金现金流收益净现值为所述独立微网系统在末年的资本金累计现金流净现值;所述污染物排放量为全寿命周期内各类污染物的排放总量;
所述经济评价环节采用工程经济学中的动态分析方法,所述动态分析方法计入项目建设期、资本金比例、银行贷款利息工程财务评估要素,获得优化结果对应的还本付息参考电价和期望收益参考电价。
作为优选实施方式,利用NSGA-II多目标遗传算法进行优化求解的过程中,选取所述独立微网系统中风力发电机的类型、风力发电机的台数、柴油发电机的类型、柴油发电机的台数,光伏电池的容量,储能电池的类型、并联支路数、电池用双向变流器的容量为优化变量。
利用NSGA-II多目标遗传算法进行优化求解的过程中,建立了随机机会约束规划模型,采用准稳态仿真策略进行模拟,以负荷容量缺失率为概率约束条件,具体为:
以小时为仿真步长,根据马尔可夫状态转换概率矩阵生成的风力资源、光照资源和负荷需求情况,对所述独立微网系统内的全部设备进行全寿命周期内的准稳态仿真,根据硬充电策略,确定每个仿真步长内,所述柴油发电机的启停状态和燃料消耗、所述储能电池的充放电功率和剩余容量,以及所述仿真步长内未满足的负荷量,计算得到全寿命周期内的资本金现金流收益净现值、污染物排放量和负荷容量缺失率,判断负荷容量缺失率是否满足预设的约束条件;
上述为一次完整的随机机会规划迭代过程,根据预设的随机机会规划迭代次数反复执行所述随机机会规划迭代过程,终止时可获得整个随机机会规划过程中的负荷容量缺失率满足约束条件的概率;若所述概率不低于预设的置信水平,认为该决策是可行的,以所述迭代过程中的平均资本金现金流收益净现值和平均污染物排放量作为优化指标;若所述概率低于预设的置信水平,认为该决策不可行,将不满足约束条件个体的相关约束的绝对值之和作为惩罚项增加到目标函数中。
获得所述独立微网系统全寿命周期内的资本金现金流收益净现值的过程中,以资本金现金流量表作为分析依据,具体为:
获取所述独立微网系统及其设备的经济参数,根据所述经济参数,计算初始投资中的资本金金额和长期银行贷款金额,根据长期银行贷款金额和系统建设期确定建设期内利息,得到含建设期利息的银行欠款总额;
对系统寿命周期中的每一年,根据售电量和售电价格确定售电收益,根据设备剩余寿命确定残值收益;根据银行欠款额、长期贷款利率和银行还款年限确定还本付息费用,根据设备替换情况确定更新费用,根据设备开启时间确定维护费用,根据柴油发电机燃料消耗和燃料价格确定燃料费用;根据收益项计算该年的现金流入,根据费用项计算该年的现金流出;
依据所述现金流入和现金流出,计算每一年的净收益,所述净收益在每年年初修正银行欠款总额,得到该年的资本金累计净现金流量,并根据折现率计算其净现值,将系统寿命周期最后一年的资本金累计净现金流量净现值作为优化中的经济指标。
所述独立微网系统优化结果的经济评价环节中,获取还本付息参考电价和期望收益参考电价的方法为:
计算银行还款年限内的系统总成本现值,所述总成本现值为资本金金额与还款年限内每年的现金流出之和;计算年售电总量,所述售电总量为每个步长能得到满足的负荷之和;根据还款期限内的总收益现值等于该期限内的总成本总现值,计算还本付息参考电价,所述还本付息参考电价反映了在还款期限内还清银行贷款的最低电价;
计算系统寿命周期内的系统总成本现值,所述总成本现值为资本金金额、寿命周期内每年的现金流出与最后一年系统残值收益之和;计算年售电总量;根据预设的期望收益值,令系统在整个寿命周期内的总盈利等于预设的期望收益值,计算期望收益参考电价,所述期望收益参考电价反映了系统达到期望盈利要求的最低电价。
本发明的技术效果包括:
1.所述利用NSGA-II多目标遗传算法进行优化求解的过程中,选取所述独立微网系统中风力发电机的类型、风力发电机的台数、柴油发电机的类型、柴油发电机的台数,光伏电池的容量,储能电池的类型、并联支路数、电池用双向变流器的容量为优化变量,从而实现了设备类型与设备容量的组合优化,大大节约了进行独立微网随机优化规划设计的时间。
2.所述利用NSGA-II多目标遗传算法进行优化求解的过程中,建立了随机机会约束规划模型,采用准稳态仿真策略进行模拟,以负荷容量缺失率为概率约束条件。利用准稳态仿真策略,实现了计及多台柴油机组合运行的柴油机/蓄电池协调运行模拟;以负荷容量缺失率为概率约束条件,有效计入了风、光等自然资源的随机性影响,提高了独立微网系统随机优化规划设计运行的可靠性。
3.获得所述独立微网系统全寿命周期内的资本金现金流收益净现值的过程中,以资本金现金流量表作为分析依据。此方法考虑了工程建设期、项目资本金、长期银行贷款利息、项目投产后还本付息等实际工程建设中的经济要素,可得到项目寿命周期中每一年内的资金流入和流出情况,充分反映出独立微网系统规划者所关心的经济细节,提高了独立微网系统随机优化规划设计中经济模型的精确性和实用性。
4.所述独立微网系统优化结果的经济评价环节中,提出了还本付息参考电价和期望收益参考电价。还本付息参考电价是在还款期限内还清银行贷款的最低电价,期望收益参考电价是令方案在寿命周期末达到盈利要求的最低电价。所述两种参考电价可为独立微网系统随机优化规划设计提供合理的收益区间,从而为独立微网系统的电价设置提供依据。
附图说明
图1为本发明独立微网系统结构拓扑图;
图2为本发明优化计算流程图。
具体实施方式
下面根据说明书附图,详细地介绍一下本发明的技术方案。
独立微网系统结构如图1所示。风力发电机、光伏发电和储能系统等通过各自的变流器接入交流微网系统。柴油发电机采用同步发电机发电,直接并入交流微网。在优化规划设计中,忽略线路阻抗引起的损耗。
本发明建立了多目标随机机会约束规划模型。机会约束规划采用置信度水平描述目标函数和约束条件,提供了一种显示表示风险程度的手段。机会约束规划是在一定概率意义下达到最优,主要用于解决约束条件中含有随机变量且必须在观测到随机变量实现之前做出的决策优化问题。随机优化的表现形式:
minf(x)
                                        (1)
s.t Pr{g(x,ξ)≤0,j=1,2,...k}>β
其中,x是n维决策变量,ξ为已知概率密度函数的随机向量,f(x)为目标函数,g(x,ξ)为随机约束函数,Pr{g(x,ξ)≤0,j=1,2,…k}>β表示g(x,ξ)中的事件成立的概率,为给定约束条件的置信水平。
选择负荷容量缺失率(LOC,Loss of Capacity)作为概率约束条件,LOC是将系统备用容量计入净负荷中,所得系统未满足的负荷容量与整个负荷需求容量的比值。计入备用容量的净负荷可表示如下:
P2(t)=Pres(t)+P1(t)=rloadPprim(t)+rpeakloadPpeakprim
                                                             (2)
+rwindPwind(t)+rsolarPpv(t)+P1(t)
其中:P1(t)表示t时刻系统净负荷;P2(t)表示t时刻系统计入备用容量的净负荷;Pres(t)为t时刻需要考虑的系统备用容量;rload为负荷的备用容量系数;Pprim(t)为t时刻的负荷有功需求;rpeakload为全年范围内峰值负荷的备用容量系数;Ppeakprim为全年范围内的峰值负荷;rwind为风力发电机的备用容量系数;Pwind(t)为t时刻风力发电机的有功出力;rsolar为光伏发电的备用容量系数;Ppv(t)为t时刻光伏发电的有功功率。
LOC作为概率约束条件表示为:
Pr ( LOC = E CS E tot ≤ 0.01 ) ≥ 80 % - - - ( 3 )
其中:ECS为总的未满足负荷容量;Etot为总的电负荷需求容量;LOC表示未满足需求的负荷容量与整个负荷需求容量的比值。LOC的取值在0~1之间,数值越小,供电可靠性越高。
优化设计重点考察两个指标:全寿命周期内的资本金现金流收益净现值和污染物排放。通常情况下,收益净现值与污染物排放存在着一定的矛盾,如果污染物排放较小,则需要的可再生能源发电单元和储能电池则越多,经济性相对较差。因此该问题为一个多目标优化问题,目标函数表示为:
minfi(X)i=1,2    (4)
全寿命周期内的资本金现金流收益净现值需通过资本金现金流分析获得。第k年累计净现金流量(收益为正)的净现值可表示为:
CNPV ( k ) = CNPV ( k - 1 ) - ( C in ( k ) - c out ( k ) ) ( 1 + r ) k k = Y c + 1 , Y c + 2 , · · · , Y c + Y l - - - ( 5 )
其中,CPNV(k)为第k年累计净现金流量的净现值,CPNV(Yc)代表项目的初始资本金投入;r为折现率;Yc为项目建设期;Yl为项目寿命周期;Cin(k)、Cout(k)分别为第k年的现金流入和现金流出,可表示为:
Cin(k)=CE(k)+CS(k)
                                                                  (6)
Cout(k)=CD(k)+CR(k)+CM(k)+CF(k)
其中,CE(k)和CS(k)分别为第k年的售电收益和残值收益;CD(k),CR(k),CM(k),CF(k)分别为第k年的还本付息、更新、维护、燃料费用。
由于资本金现金流表得到的是项目的收益现值而非成本现值,本方法选择项目寿命周期最后一年的累计净现金流量净现值的负值作为目标函数,即:
f1(X)=-CNPV(Yc+Yl)    (7)
柴油发电机以柴油作为主要燃料,其污染物包括CO、CO2、燃料中未燃烧尽的碳氢化合物、硫化物和一氧化氮。污染物的排放量与燃料消耗量直接相关。微网年排放水平表示如下:
f 2 ( X ) = Σ k = 1 K ( σ CO 2 + σ CO + σ HC + σ NO + σ S ) v fule ( k ) - - - ( 8 )
其中:σCO2、σCO、σHC、σNO、σS分别代表各污染物的排放系数[kg/L];vfule(k)代表第k年柴油发电机的年消耗量,单位是升。
选择独立微网中的风机类型Windtype、风机台数WindNum、柴油发电机的类型Diseltype、柴油发电机的台数Diselnum,光伏电池的容量PVcap,储能电池的类型Battype、并联支路数Batparrel(此处假定电池的串联个数根据逆变器的直流工作电压已确定)、电池用双向变流器的容量Concap为优化变量。定义优化变量为:
X=[Windtype,WindNum,Diseltype,Diselnum,PVcap,
                                             (9)
Battype,BatParrel,Concap]
柴油发电机与储能系统共同使用时,采用柴油发电机与蓄电池可轮流做主电源满足净负荷需求的运行策略。本发明采用硬充电策略作为独立微网系统的控制策略。该策略的主要目的是尽量减少柴油机组的运行时间,一般适用于因环境或能源约束对柴油机组运行时间有限制的场合,以及可再生能源资源充沛的地区。该策略下,柴油发电机与蓄电池可轮流做主电源满足净负荷(本发明中净负荷表示为可再生能源出力减去负荷)需求,允许柴油机给蓄电池充电。
本发明在经济评价环节中提供了还本付息参考电价和期望收益参考电价,为独立微网系统的电价设置提供依据,具体方法如下。
还本付息参考电价是在还款期限内还清银行贷款的最低电价,考虑到资金的时间价值(即折现率)的影响,还本付息参考电价应使还款期限内的总收益现值等于该期限内的总成本总现值:
Σ k = Y c + 1 Y c + Y r P r · E ( k ) ( 1 + r ) k = CNPV ( Y c ) + Σ k = Y c + 1 Y c + Y r C out ( k ) ( 1 + r ) k - - - ( 10 )
式中,Pr为还本付息参考电价;E(k)为第k年的总售电量;Yr为银行规定的还款期限。
期望收益参考电价是令方案在寿命周期末达到盈利要求的最低电价,它应使项目整个寿命期限内的总盈利应等于投资者的期望收益值:
P b = CNPV ( Y c ) + Σ k = Y c + 1 Y c + Y l ( C out ( k ) - C s ( k ) ) + R expect Σ k = Y c + 1 Y c + Y l E ( k ) - - - ( 11 )
式中,Pb为项目期望收益参考电价;Rexpect为还本付息期结束后的期望收益值。
本发明基于随机数据(风速、清洁指数、负荷)的概率密度函数,采用马尔科夫状态转移概率矩阵实现数据的随机模拟,具体方法如下。
风速数据采用韦布尔概率密度函数生成,表达式为:
p ( W ) = ( k c ) ( W c ) k - 1 exp ( - ( W c ) k ) - - - ( 12 )
k和c的的计算公式如下:
k = ( σ W ‾ ) - 1.086 - - - ( 13 )
c = W ‾ Γ ( 1 + 1 / k ) - - - ( 14 )
其中:
Figure BDA0000474400350000065
表示风速平均值;σ是标准差;Γ(·)表示伽玛函数。
负荷采用原始数据乘以扰动因子得到。输入数据为一年不同月份每一天24小时负荷值,共12×24个数值。扰动系数公式如下:
α=1+δdh    (15)
其中:δd为日扰动因子,δh为时扰动因子。日扰动因子服从正态分布,平均值为零。时扰动因子也服从正态分布,平均值为零。
由于光照强度具有非常强烈的昼夜模式,不能直接由概率密度函数建模。本发明采用清洁指数代表天气状况的不确定性,进而表示光照强度的随机特性。清洁指数指地表的水平光照总量I与地球大气顶层的平均小时光照量I0之比:
k T = I I o - - - ( 16 )
地表的水平光照总量I公式如下:
I=Ib+Id    (17)
其中,其中Ib表示束辐射,Id漫辐射,其表达如下:
I d I = 1.0 - 0.09 k T for k T &le; 0.22 0.9511 - 0.1604 k T + 4.388 k T 2 - 16.638 k T 3 + 12.336 k T 4 for 0.22 < k T < 0.80 0.165 for k T > 0.80 - - - ( 18 )
I T = ( I b + I d I b I o ) cos &theta; cos &theta; z + I d ( 1 - I b I o ) ( 1 + cos &beta; 2 ) [ 1 + I b I si n 3 ( &beta; 2 ) ] + I&rho; g ( 1 - cos &beta; 2 ) - - - ( 19 )
其中,θ太阳光于任意时刻到达任意朝向的光伏板表面的入射角,θz太阳光和地表水平面的入射角,β为表面倾斜角,ρg为地面反射系数。
大气清洁指数的概率密度分布函数表示为:
P ( k t ) = C ( k tu - k t ) k tu exp ( &lambda; k t ) - - - ( 20 )
其中:
Figure BDA0000474400350000073
表示平均值;最大清洁指数ktu理论计算值为0.864。系数C的表达式为:
C = &lambda; 2 k tu / ( e &lambda; k tu - 1 - &lambda; k tu ) - - - ( 21 )
系数λ的表达式为:
λ=(2τ-17.519exp(-1.3118τ)
                                                            (22)
-1062exp(-5.0426τ))/ktu
τ的表达式为:
&tau; = k tu / ( k tu - k t &OverBar; ) - - - ( 23 )
由于风速,负荷,清洁指数存在着时序性,即下一时刻的数据与上一时刻的数据有一定的关联,故本发明采用马尔可夫状态转移概率矩阵的方法生成随机数据。
马尔可夫状态转移概率矩阵[T]由三个矩阵的乘积构成,从左至右分别是标准化矩阵[N],衰减矩阵[G]和初始概率矩阵[P]。其中,[N]、[P]是对角元素矩阵。若某时刻状态为i,下一刻状态转化为j的概率即为Ti,j。由于状态i在下一时刻一定会转移到[T]矩阵中i行的某一个状态中,故转移矩阵[T]中每行元素之和为1。
衰减矩阵[G]中的元素Gi,j如下所示:
Gi,j=gbase∧(i-j)    (24)
其中:gbase表示初始基准值。
标准化矩阵[N]为对角矩阵,其作用是使得[T]中每一行元素之和为1。[N]中元素Ni,i表达式如下:
N i , i = 1 &Sigma; j = 1 K p i , i g i , j - - - ( 25 )
先根据概率模型得到稳态分布向量{R},{R}中的元素为不同数据下的概率值,然后采用如下方法得到初始概率矩阵[P]:
首先假设初始概率向量{P}等于稳态分布向量{R},过渡稳态向量{Rt}表示为:
{Rt}=[P][G]{P}    (26)
其中:{P}中元素与[P]中的对角线元素相等。
其次迭代修正初始概率向量{P},即:
{P}={P}+0.9×({R}-{Rt})    (27)
当稳态分布向量与过渡稳态向量之间差值小于某一个阀值时,停止迭代。此时的向量{P}即为所需的初始概率向量,即[P]上的对角线元素。然后便可利用此初始概率矩阵[P]生成所需的转移矩阵[T]:
[T]=[N][G][P]    (28)
根据[T]生成累积概率矩阵,从而生成所需的随机数据。假设累积概率矩阵为[C],其元素为Ci,j,则:
Ci,1=Ti,1,Ci,j=Ci,j-1+Ti,j    (29)
随机生成一个在(0,1)上均匀分布的数x,假设此时状态为i,若x≥Ci,j且x<Ci,j+1,则下一时刻状态由i转移到j。
自相关系数用于判断某一时间间隔前后数据的线性相关性,公式如下:
R ( r&Delta;t ) = 1 &sigma; x ( N - r ) &Sigma; i = 1 N - r ( x i - x &OverBar; ) ( x i + r - x &OverBar; ) - - - ( 30 )
其中:N为生成的总数据数,r是所选的时间间隔,σx为生成数据的方差,xi、xi+r分别为生成的随机数据,为生成随机数据平均值。若得到的自相关系数过小,则可通过增大gbase,从而使得某时刻的状态i在下一时刻转移到i附近的概率增大,即增大了两者之间的线性相关性,反之,减小gbase
本发明采用了基于NSGA-II的多目标遗传算法进行求解。首先,对当前种群P个体进行Pareto分层排序,即按个体支配关系进行分层并计算聚集距离,排序时优先考虑层级关系,同层则按聚集距离大小排序,层级低、聚集距离大的个体优先。然后对种群个体进行遗传算法的选择、交叉、变异和重插操作,得到子代种群Q。其中选择算子采用锦标赛法,比赛规模为2个,比赛中优先选取Pareto分层排序较高的个体;交叉算子采用单点交叉方式;变异算子采用均匀变异方式。最后,将当前种群P与子代种群Q合并,重新进行Pareto分层排序操作,按排序次序选择满足群体规模的个体组成新一代种群P。重复上述操作,直至满足终止条件。本发明将进化代数作为进化算法终止判据。最后一代种群的所有非支配解构成的集合,即所研究问题的Pareto最优解集;各非支配解对应的目标函数值构成的解集,即为该问题的Pareto最优边界,也即本发明多目标问题的最优解集。
多目标优化计算的求解流程如图2所示,描述如下:
1)系统初始化。读取系统中蓄电池、风机、光伏、柴油发电机、变流器等设备和遗传算法的参数。
2)初始化种群P,通过随机函数产生N个个体,作为初始种群P。
3)计算种群P适应度函数值:
首先进行随机机会约束规划,在一定的迭代次数内,利用马尔可夫状态转移概率矩阵生成全年的逐时数据(风速、清洁指数和负荷),对种群中的每个个体调用控制策略仿真计算,记录每个个体的目标函数和LOC满足次数。满足迭代次数后对目标函数求取平均值,计算约束值。然后按照式(26)计算各个体适应度函数值。
其中:f1,max(X)表示所有个体的第1个目标函数值中的最大值,f2,max(X)表示所有个体的第2个目标函数值中的最大值,Δ为不满足约束条件个体的相关约束的绝对值之和。
4)对初始种群P进行Pareto排序操作
a.按个体目标函数值计算各个体之间的支配关系。
b.计算各个体间聚集距离。
c.按个体间支配关系和聚集距离执行Pareto排序操作。
5)选择、交叉、遗传操作
a.通过轮盘赌转法从当前(父代)种群P中选择交叉、变异种群P2。
b.通过交叉、变异算子更新种群P2的子代种群P3。
6)计算种群P3适应度函数值。
7)对父代种群和子代种群P∪P3进行Pareto排序。
8)更新种群P。按Pareto排序结果,从种群P∪P3选择N个个体形成新的种群P。
9)终止条件:本发明将进化代数作为进化算法终止判据。最后一代种群的所有非支配解构成的集合,即所研究问题的Pareto最优解集;各非支配解对应的目标函数值构成的解集,即为该问题的Pareto最优边界。若满足,则对最后一代种群调用准稳态仿真策略后输出最终的优化结果,否则返回5)。

Claims (5)

1.一种独立微网系统的随机优化规划方法,其特征在于,包括:
获取风力发电机、光伏电池、柴油发电机和所述储能电池的设备参数;
基于风速、大气清洁指数、负荷随机数据的概率模型,采用马尔可夫状态转换概率矩阵模拟上述随机数据;
根据所述设备参数和所述随机数据,利用NSGA-II多目标遗传算法进行优化求解,在求解过程中采用随机机会约束规划,获得所述独立微网系统的优化指标,最终通过所述多目标遗传算法获得优化结果,并依据经济评价环节获取优化结果的具体经济性信息;
所述随机机会约束规划采用准稳态仿真策略进行模拟,以负荷容量缺失率为概率约束条件;所述负荷容量缺失率指将系统备用容量计入净负荷中,所得系统未满足的负荷容量与整个负荷需求容量的比值;
所述独立微网系统的优化指标为全寿命周期内的资本金现金流收益净现值和污染物排放量,所述全寿命周期内的资本金现金流收益净现值为所述独立微网系统在末年的资本金累计现金流净现值;所述污染物排放量为全寿命周期内各类污染物的排放总量;
所述经济评价环节采用工程经济学中的动态分析方法,所述动态分析方法计入项目建设期、资本金比例、银行贷款利息工程财务评估要素,获得优化结果对应的还本付息参考电价和期望收益参考电价。
2.根据权利要求1所述的方法,其特征在于,利用NSGA-II多目标遗传算法进行优化求解的过程中,选取所述独立微网系统中风力发电机的类型、风力发电机的台数、柴油发电机的类型、柴油发电机的台数,光伏电池的容量,储能电池的类型、并联支路数、电池用双向变流器的容量为优化变量。
3.根据权利要求1所述的方法,其特征在于,利用NSGA-II多目标遗传算法进行优化求解的过程中,建立了随机机会约束规划模型,采用准稳态仿真策略进行模拟,以负荷容量缺失率为概率约束条件,具体为:
以小时为仿真步长,根据马尔可夫状态转换概率矩阵生成的风力资源、光照资源和负荷需求情况,对所述独立微网系统内的全部设备进行全寿命周期内的准稳态仿真,根据硬充电策略,确定每个仿真步长内,所述柴油发电机的启停状态和燃料消耗、所述储能电池的充放电功率和剩余容量,以及所述仿真步长内未满足的负荷量,计算得到全寿命周期内的资本金现金流收益净现值、污染物排放量和负荷容量缺失率,判断负荷容量缺失率是否满足预设的约束条件;
上述为一次完整的随机机会规划迭代过程,根据预设的随机机会规划迭代次数反复执行所述随机机会规划迭代过程,终止时可获得整个随机机会规划过程中的负荷容量缺失率满足约束条件的概率;若所述概率不低于预设的置信水平,认为该决策是可行的,以所述迭代过程中的平均资本金现金流收益净现值和平均污染物排放量作为优化指标;若所述概率低于预设的置信水平,认为该决策不可行,将不满足约束条件个体的相关约束的绝对值之和作为惩罚项增加到目标函数中。
4.根据权利要求1所述的方法,其特征在于,获得所述独立微网系统全寿命周期内的资本金现金流收益净现值的过程中,以资本金现金流量表作为分析依据,具体为:
获取所述独立微网系统及其设备的经济参数,根据所述经济参数,计算初始投资中的资本金金额和长期银行贷款金额,根据长期银行贷款金额和系统建设期确定建设期内利息,得到含建设期利息的银行欠款总额;
对系统寿命周期中的每一年,根据售电量和售电价格确定售电收益,根据设备剩余寿命确定残值收益;根据银行欠款额、长期贷款利率和银行还款年限确定还本付息费用,根据设备替换情况确定更新费用,根据设备开启时间确定维护费用,根据柴油发电机燃料消耗和燃料价格确定燃料费用;根据收益项计算该年的现金流入,根据费用项计算该年的现金流出;
依据所述现金流入和现金流出,计算每一年的净收益,所述净收益在每年年初修正银行欠款总额,得到该年的资本金累计净现金流量,并根据折现率计算其净现值,将系统寿命周期最后一年的资本金累计净现金流量净现值作为优化中的经济指标。
5.根据权利要求1所述的方法,其特征在于,所述独立微网系统优化结果的经济评价环节中,获取还本付息参考电价和期望收益参考电价的方法为:
计算银行还款年限内的系统总成本现值,所述总成本现值为资本金金额与还款年限内每年的现金流出之和;计算年售电总量,所述售电总量为每个步长能得到满足的负荷之和;根据还款期限内的总收益现值等于该期限内的总成本总现值,计算还本付息参考电价,所述还本付息参考电价反映了在还款期限内还清银行贷款的最低电价;
计算系统寿命周期内的系统总成本现值,所述总成本现值为资本金金额、寿命周期内每年的现金流出与最后一年系统残值收益之和;计算年售电总量;根据预设的期望收益值,令系统在整个寿命周期内的总盈利等于预设的期望收益值,计算期望收益参考电价,所述期望收益参考电价反映了系统达到期望盈利要求的最低电价。
CN201410083264.8A 2014-03-07 2014-03-07 独立微网系统的随机优化规划方法 Pending CN103854070A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410083264.8A CN103854070A (zh) 2014-03-07 2014-03-07 独立微网系统的随机优化规划方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410083264.8A CN103854070A (zh) 2014-03-07 2014-03-07 独立微网系统的随机优化规划方法

Publications (1)

Publication Number Publication Date
CN103854070A true CN103854070A (zh) 2014-06-11

Family

ID=50861700

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410083264.8A Pending CN103854070A (zh) 2014-03-07 2014-03-07 独立微网系统的随机优化规划方法

Country Status (1)

Country Link
CN (1) CN103854070A (zh)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104683144A (zh) * 2014-12-03 2015-06-03 北京邮电大学 一种基于马尔可夫模型的物联网设备依赖度评价方法
CN105182749A (zh) * 2015-09-02 2015-12-23 湘潭大学 基于马尔科夫决策过程的变流器统一控制策略的方法
CN105406474A (zh) * 2015-12-25 2016-03-16 广东电网有限责任公司中山供电局 含分布式电源的配电网设备的利用率检测方法及其系统
CN105406507A (zh) * 2015-12-07 2016-03-16 浙江工业大学 一种光伏微网微源动态切换方法
CN105976108A (zh) * 2016-05-05 2016-09-28 国网浙江省电力公司电力科学研究院 一种配电网分布式储能规划方法
CN106100002A (zh) * 2016-07-28 2016-11-09 华北电力大学(保定) 一种交直流混合微网的优化运行方法
CN105226688B (zh) * 2015-10-12 2018-11-27 中国电力科学研究院 基于机会约束模型的多类型储能系统容量优化配置方法
CN109102185A (zh) * 2018-08-07 2018-12-28 中国电力科学研究院有限公司 一种评估光储充一体化电站全寿命周期经济性的方法和系统
CN109767025A (zh) * 2017-11-10 2019-05-17 通用电气公司 用以生成优化的工作范围的设备、存储介质和方法
CN109783909A (zh) * 2018-12-29 2019-05-21 天津大学 考虑多种储能设备的综合能源站设备选型与容量规划方法
CN111709695A (zh) * 2020-04-16 2020-09-25 青海绿能数据有限公司 一种光伏运营管理系统
CN112036655A (zh) * 2020-09-07 2020-12-04 南通大学 基于机会约束的光伏电站与电动汽车充电网络规划方法
CN113315156A (zh) * 2021-05-31 2021-08-27 广东电网有限责任公司 动力电池的入网配置方法、装置、计算机设备和存储介质
CN114944698A (zh) * 2022-04-29 2022-08-26 南通电博士自动化设备有限公司 一种智能的柴油发电机群组的控制方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101045326B1 (ko) * 2009-09-29 2011-06-30 한국전력공사 마이크로그리드 운영 시스템 및 방법
US20130024014A1 (en) * 2011-07-20 2013-01-24 Nec Laboratories America, Inc. Optimal energy management of a rural microgrid system using multi-objective optimization
CN103151798A (zh) * 2013-03-27 2013-06-12 浙江省电力公司电力科学研究院 独立微网系统的优化方法
CN103166248A (zh) * 2013-04-01 2013-06-19 国电联合动力技术有限公司 一种独立风-柴-储微电网系统容量的工程配置方法
CN103544372A (zh) * 2013-08-05 2014-01-29 中国科学院电工研究所 一种微型燃气轮机热电联供系统的机组选型方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101045326B1 (ko) * 2009-09-29 2011-06-30 한국전력공사 마이크로그리드 운영 시스템 및 방법
US20130024014A1 (en) * 2011-07-20 2013-01-24 Nec Laboratories America, Inc. Optimal energy management of a rural microgrid system using multi-objective optimization
CN103151798A (zh) * 2013-03-27 2013-06-12 浙江省电力公司电力科学研究院 独立微网系统的优化方法
CN103166248A (zh) * 2013-04-01 2013-06-19 国电联合动力技术有限公司 一种独立风-柴-储微电网系统容量的工程配置方法
CN103544372A (zh) * 2013-08-05 2014-01-29 中国科学院电工研究所 一种微型燃气轮机热电联供系统的机组选型方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘超: ""不同运营模式的微网综合效益评价研究"", 《中国优秀硕士学位论文全文数据库 经济与管理科学辑》 *
郭力等: ""独立微网系统的多目标优化规划设计方法"", 《中国电机工程学报》 *

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104683144A (zh) * 2014-12-03 2015-06-03 北京邮电大学 一种基于马尔可夫模型的物联网设备依赖度评价方法
CN104683144B (zh) * 2014-12-03 2018-01-09 北京邮电大学 一种基于马尔可夫模型的物联网设备依赖度评价方法
CN105182749A (zh) * 2015-09-02 2015-12-23 湘潭大学 基于马尔科夫决策过程的变流器统一控制策略的方法
CN105226688B (zh) * 2015-10-12 2018-11-27 中国电力科学研究院 基于机会约束模型的多类型储能系统容量优化配置方法
CN105406507B (zh) * 2015-12-07 2018-01-12 浙江工业大学 一种光伏微网微源动态切换方法
CN105406507A (zh) * 2015-12-07 2016-03-16 浙江工业大学 一种光伏微网微源动态切换方法
CN105406474A (zh) * 2015-12-25 2016-03-16 广东电网有限责任公司中山供电局 含分布式电源的配电网设备的利用率检测方法及其系统
CN105976108A (zh) * 2016-05-05 2016-09-28 国网浙江省电力公司电力科学研究院 一种配电网分布式储能规划方法
CN106100002A (zh) * 2016-07-28 2016-11-09 华北电力大学(保定) 一种交直流混合微网的优化运行方法
CN109767025A (zh) * 2017-11-10 2019-05-17 通用电气公司 用以生成优化的工作范围的设备、存储介质和方法
CN109767025B (zh) * 2017-11-10 2023-09-12 通用电气公司 用以生成优化的工作范围的设备、存储介质和方法
CN109102185A (zh) * 2018-08-07 2018-12-28 中国电力科学研究院有限公司 一种评估光储充一体化电站全寿命周期经济性的方法和系统
CN109102185B (zh) * 2018-08-07 2022-04-15 中国电力科学研究院有限公司 一种评估光储充一体化电站全寿命周期经济性的方法和系统
CN109783909A (zh) * 2018-12-29 2019-05-21 天津大学 考虑多种储能设备的综合能源站设备选型与容量规划方法
CN111709695A (zh) * 2020-04-16 2020-09-25 青海绿能数据有限公司 一种光伏运营管理系统
CN112036655A (zh) * 2020-09-07 2020-12-04 南通大学 基于机会约束的光伏电站与电动汽车充电网络规划方法
CN113315156A (zh) * 2021-05-31 2021-08-27 广东电网有限责任公司 动力电池的入网配置方法、装置、计算机设备和存储介质
CN114944698A (zh) * 2022-04-29 2022-08-26 南通电博士自动化设备有限公司 一种智能的柴油发电机群组的控制方法及系统
CN114944698B (zh) * 2022-04-29 2024-05-03 湛江伟力机电设备有限公司 一种智能的柴油发电机群组的控制方法及系统

Similar Documents

Publication Publication Date Title
CN103854070A (zh) 独立微网系统的随机优化规划方法
Wang et al. Optimal design of integrated energy system considering economics, autonomy and carbon emissions
Pavičević et al. The potential of sector coupling in future European energy systems: Soft linking between the Dispa-SET and JRC-EU-TIMES models
Al-Maamary et al. The impact of oil price fluctuations on common renewable energies in GCC countries
Zhang et al. Capacity configuration optimization of multi-energy system integrating wind turbine/photovoltaic/hydrogen/battery
Raza et al. Holistic and scientific approach to the development of sustainable energy policy framework for energy security in Pakistan
Wen et al. Optimizing the sizes of wind and photovoltaic plants complementarily operating with cascade hydropower stations: Balancing risk and benefit
Elsied et al. An advanced energy management of microgrid system based on genetic algorithm
CN107464010A (zh) 一种虚拟电厂容量优化配置方法
CN105977991A (zh) 一种考虑价格型需求响应的独立型微网优化配置方法
CN104392394B (zh) 一种微电网储能裕度的检测方法
CN103440535B (zh) 基于免疫优化和模糊决策的多目标厂级负荷优化方法
CN103997062A (zh) 一种微网能量管理控制方法
CN108053057A (zh) 一种基于CVaR的虚拟电厂优化调度建模方法
Aliari et al. Planning for integration of wind power capacity in power generation using stochastic optimization
Li et al. Optimal planning of Electricity–Hydrogen hybrid energy storage system considering demand response in active distribution network
Khosa et al. Performance analysis based on probabilistic modelling of Quaid-e-Azam Solar Park (QASP) Pakistan
Liu et al. A framework for uncertainty and sensitivity analysis of district energy systems considering different parameter types
Diachuk et al. Transition of Ukraine to the Renewable Energy by 2050
CN113779874A (zh) 一种离网微电网建设的多目标优化方法
Xue et al. Optimal planning for distributed energy systems with carbon capture: Towards clean, economic, independent prosumers
Khan et al. Integrative decision-making framework for techno-economic planning and sustainability assessment of renewable dominated standalone hybrid microgrids infrastructure at provincial scale of Pakistan
Betancourt-Torcat et al. Design multiperiod optimization model for the electricity sector under uncertainty–A case study of the Emirate of Abu Dhabi
Gao et al. Thermoelectric optimization of integrated energy system considering wind-photovoltaic uncertainty, two-stage power-to-gas and ladder-type carbon trading
Al-Masri et al. Optimal configuration of a large scale on-grid renewable energy systems with different design strategies

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information

Inventor after: Guo Li

Inventor after: Wang Chengshan

Inventor after: Yu Zhouzi

Inventor after: Yang Yong

Inventor after: Zhao Bo

Inventor after: Zhang Xuesong

Inventor after: Zhou Dan

Inventor before: Guo Li

Inventor before: Wang Chengshan

Inventor before: Yu Zhouzi

Inventor before: Zhao Bo

Inventor before: Zhang Xuesong

Inventor before: Zhou Dan

COR Change of bibliographic data
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20140611

WD01 Invention patent application deemed withdrawn after publication