CN107749638A - 多微电网组合的虚拟发电厂分布式随机非重叠抽样的无中心优化方法 - Google Patents
多微电网组合的虚拟发电厂分布式随机非重叠抽样的无中心优化方法 Download PDFInfo
- Publication number
- CN107749638A CN107749638A CN201710976876.3A CN201710976876A CN107749638A CN 107749638 A CN107749638 A CN 107749638A CN 201710976876 A CN201710976876 A CN 201710976876A CN 107749638 A CN107749638 A CN 107749638A
- Authority
- CN
- China
- Prior art keywords
- power source
- microgrid
- distributed power
- intermittent
- distributed
- 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.)
- Granted
Links
- 238000005070 sampling Methods 0.000 title claims abstract description 51
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000005457 optimization Methods 0.000 title claims abstract description 28
- 241000039077 Copula Species 0.000 claims abstract description 80
- 238000009826 distribution Methods 0.000 claims abstract description 54
- 238000010219 correlation analysis Methods 0.000 claims abstract description 10
- 238000004458 analytical method Methods 0.000 claims abstract description 8
- 230000005611 electricity Effects 0.000 claims description 27
- 238000004146 energy storage Methods 0.000 claims description 19
- 238000012423 maintenance Methods 0.000 claims description 12
- 239000011159 matrix material Substances 0.000 claims description 11
- 230000001419 dependent effect Effects 0.000 claims description 9
- 238000005315 distribution function Methods 0.000 claims description 8
- 230000001186 cumulative effect Effects 0.000 claims description 7
- 238000010248 power generation Methods 0.000 claims description 7
- 238000012546 transfer Methods 0.000 claims description 7
- 230000005540 biological transmission Effects 0.000 claims description 6
- 238000005553 drilling Methods 0.000 claims description 6
- 238000004891 communication Methods 0.000 claims description 4
- 238000003860 storage Methods 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims description 3
- 239000004744 fabric Substances 0.000 claims description 3
- 230000035945 sensitivity Effects 0.000 claims description 3
- 238000005303 weighing Methods 0.000 claims description 3
- 239000000446 fuel Substances 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 claims 1
- 230000005684 electric field Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000002775 capsule Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 239000000295 fuel oil Substances 0.000 description 1
- 238000002513 implantation Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 241000894007 species Species 0.000 description 1
Classifications
-
- H02J3/383—
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- H02J3/386—
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02B—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
- Y02B10/00—Integration of renewable energy sources in buildings
- Y02B10/10—Photovoltaic [PV]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/50—Photovoltaic [PV] energy
- Y02E10/56—Power conversion systems, e.g. maximum power point trackers
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Marketing (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种多微电网组合的虚拟发电厂分布式随机非重叠抽样的无中心优化方法,包括考虑间歇性分布式电源的影响,建立分布式控制的虚拟发电厂分层随机优化模型;利用Copula理论分析间歇性分布式电源两两之间的相关性;利用分布式随机非重叠抽样法对间歇性分布式电源相关性分析后的结果进行抽样,从而将虚拟发电厂分层随机优化模型转换为分层确定性优化模型来求解。该方法考虑了虚拟发电厂所含的间歇性分布式电源的影响,运用Copula理论对间歇性分布式电源进行相关性分析,将虚拟电厂分层随机优化模型转换为分层确定性优化模型来求解;并对间歇性分布式电源的历史数据进行经验分布估计,提高了预测精度。
Description
技术领域
本发明涉及分布式电源优化调度方法,特别是涉及多微电网组合的虚拟发电厂分布式随机非重叠抽样的无中心优化方法。
背景技术
近年来,随着能源的大量使用,全球化能源日趋紧张,环境破坏日益严重,且造成了温室效应。因此,人类开始反思自身发展模式,降低能源需求和减少排放的理念逐渐深入人心,并且各种以节能减排为主题的绿色行动在世界范围内开展起来。在此背景下,分布式发电(Distributed Generation,DG)得到了世界各国的重视。
然而,一方面由于分布式能源自身的可控性小,或者通讯和控制硬件不足;另一方面由于电网中存在数量庞大且出力不确定的分布式电源或微网群,使得电网优化的数学模型非常复杂,这给电网调度计划的制定带来巨大的困难。虚拟发电厂(Virtual PowerPlant,VPP)提供了一种处于分布式电源与电网中间层的能源管理形式,虚拟发电厂的提出可消除单个分布式电源并网运行的随机性与不可调度性,宏观上表现出分布式电源和用户负荷的可控性。电网将对分布式电源的协调控制权移交给控制中心,摆脱了分布式电源管理问题维数庞大、控制复杂的难题。分布式电源和用户以微网等形式参与电网统一调度,可有效解决大量分布式电源单独并网运行难以调度管理的问题,是智能电网的发展方向之一。
但是,在对虚拟发电厂进行调度的过程中,其所含的间歇性分布式电源可能会产生一些随机性、不确定性,从而使得优化调度的结果不准确。
发明内容
发明目的:针对现有技术的不足,提供一种能够提高对含随机性出力的虚拟发电厂优化控制的精度的多微电网组合的虚拟发电厂分布式随机非重叠抽样的无中心优化方法。
技术方案:一种多微电网组合的虚拟发电厂分布式随机非重叠抽样的无中心优化方法,包括以下步骤:
(1)考虑间歇性分布式电源的影响,建立分布式控制的虚拟发电厂分层随机优化模型;
(2)利用Copula理论分析间歇性分布式电源两两之间的相关性;
(3)利用分布式随机非重叠抽样法对间歇性分布式电源相关性分析后的结果进行抽样,从而将虚拟发电厂分层随机优化模型转换为分层确定性优化模型来求解。
进一步的,所述步骤(1)包括:
(11)以虚拟发电厂全局与主电网之间的交换功率总成本最小化为目标,以电能交换约束和传输线路限制为约束条件,建立上层优化模型;
上层目标为:
其中,TC表示虚拟发电厂全局与主电网之间的交换功率的总成本,IMG是虚拟发电厂分布式调度控制的微网数量,ρ是每个微网与主电网之间的能源交换价格,是虚拟发电厂内部微网i与主电网之间的交流电量,表示主电网将电能转移到微网 i,否则,表示微网i将电能转移到主电网;
上层约束为:
3)电能交换的上限约束和下限约束
其中,i=1,2,...,IMG,和分别代表主电网和微网i之间的电量交换的下限和上限;
4)传输线路限制约束
其中,m=1,2,...M,M为最大线路数,ηi-m是电力从微网注入到线路的灵敏度,Tm表示线路m的能量流限制;
(12)以虚拟发电厂内各微网发电成本最小为目标,以微网内各发电单元出力为约束,同时考虑间歇性分布式电源发电出力的不确定性,建立下层随机优化模型;
下层目标:
其中,j=1,2,...,Jmt,Ji是微网i分布式发电机的数量,是微网i的发电成本,表示微网i中的发电单元j的输出功率,表示微网i中的发电单元j的单位成本,则是微网i发电单元j的发电成本函数;
下层约束:
4)功率平衡约束
其中,是微网i中风力发电的出力,是微网i中的光伏发电的输出,是微型涡轮机的输出,和分别是储能的充电能力和放电能力,和分别是储能充电的实际利用率和放电的实际利用率,是微网i的负荷;
5)可控分布式电源的上下限输出
其中,j=1,2,...,Jmt,和分别是微网i中可控分布式电源j出力的上下限, Jmt为微网i中可控分布式电源总数;
6)储能约束
SOCi min≤SOC i≤SOCi max
其中,和是充电功率的上下限,和分别是微网i中放电功率的上下限,SOCi是微网i中电池的存储状态,和分别是电池储能的上限和下限;
4)虚拟电厂内各微电网中间歇性分布式电源发电受自然条件的影响具有较强的波动性与随机性,且同一地区地理位置两两发电单元之间的输出功率具有一定的相关性,故在此采用机会约束的形式表示其运行约束,第i个微网内间歇性分布式电源发电的概率约束表达式如下所示:
其中,为微网i中可控分布式电源k的功率出力,它满足Copula联合概率分布函数;是虚拟发电厂i中可控分布式电源k出力的上下限;Ki为微网i中可控分布式电源总数;α为机会约束条件的置信水平,该值是预先指定己知的。
优选的,所述微网内各发电单元包括可控分布式电源、间歇性分布式电源和储能,所述可控分布式电源为微型涡轮机。
更进一步的,所述步骤(12)中间歇性分布式电源包括风力发电和光伏发电,则下层目标具体为:
其中,是微网i中风力发电的出力,是风电设备的运行和维护成本,是微网i中的光伏发电的输出,是光伏发电设备的运行和维护成本,是微型涡轮机的运行和维护成本,是微型涡轮机的输出,是微型涡轮机的燃油成本,是储能设备的运行和维护成本,和分别是储能的充电能力和放电能力。
进一步的,所述步骤(2)中对各间歇性分布式电源的历史数据进行经验分布估计,得到各自的经验累积分布函数,根据平方欧式距离最小的方法选择合适的理论Copula 函数,从而建立各间歇性分布式电源两两之间的联合概率分布来描述其相关性;其相关性分析包括以下步骤:
(21)对同一地区相邻的间歇性分布式电源进行两两选择,根据某一间歇性分布式电源出力历史数据,计算不同Copula函数相应的参数估计值;
(22)分别根据每种不同的Copula函数计算其与经验Copula函数之间的平方欧氏距离,选择距离最小的Copula函数作为描述两两间歇性分布式电源出力之间相关性;
其中经验Copula函数定义如下:
式中:u,v∈[0,1],I[·]为示性函数,当F(xi)≤u时,否则当G(yi)≤v时,否则
进一步的,所述步骤(3)中在间歇性分布式电源相关性分析的基础上,利用分布式随机非重叠抽样方法对所选择的Copula函数生成满足间歇性分布式电源出力相关性的随机数矩阵进行抽样,步骤如下:
(31)对某一间歇性分布式电源进行独立抽样,建立其与其他各间歇性分布式电源的两两联合概率分布,根据平方欧氏距离最小的方法选择Copula函数的类型来作为描述间歇性分布式电源出力相关性的最优函数;
计算平方欧氏距离的表达式如下:
式中:F(x)和G(y)分别是两个间歇性分布式电源出力的边缘累积分布函数; CEi(F(x),G(y))为经验Copula函数上第i个采样值;CTi(F1(x),G(y))为Copula函数上第i个采样值;N为计算点个数;
(32)由所选的理论Copula函数的类型描述间歇性分布式电源两两之间的相关性,得出各自的相关性系数;为了避免相关性太大对出力预测精度的影响,通过分布式筛选,在众多间歇性分布式电源联合Copula函数中选择Pearson线性相关系数最小的间歇性分布式电源联合Copula函数作为抽样;
Pearson线性相关系数表示如下:
其中,X、Y为两个对象的样本,
假设某一间歇性分布式电源i与间歇性分布式电源j1、j2、.......、jn的相关系数分别为ρij1、ρij2、.......、ρijn,通过各间歇性分布式电源之间的通信连接,可以传递有链接分布式电源之间的信息;由此,为了避免相关性太大对出力预测精度的影响,可以选择与间歇性分布式电源i Pearson线性相关系数最小的另一间歇性分布式电源建立联合概率分布作为抽样的样本;
ρi抽样=min|ρij|j=j1,j2,...,jn
其中,ρi抽样对应最终所选的两个间歇性分布式电源联合概率分布对应的相关系数,ρij表示间歇性分布式电源i与其它各个间歇性分布式电源所建立的联合概率分布的相关系数;
(33)利用所选择的Copula函数生成满足间歇性分布式电源出力相关性的随机数矩阵;
假设针对每个间歇性分布式电源出力随机变量生成随机数的个数为N,利用所选择的Copula函数生成满足间歇性分布式电源出力相关性的随机数矩阵为DN×K;
其中,di,j表示节点j对节点i的影响权重,di,j∈[0,1);
(34)用Copularnd()函数对样本进行分布式随机抽样;
D=Copularnd(Copula类型,Copula相关系数,抽样次数)
(35)对抽样结果进行所需处理。
更进一步的,基于分布式随机非重叠抽样将随机约束条件转换成确定性的约束条件,步骤如下:
(a)从Copula联合概率分布式中抽样获得N组间歇性分布式电源的样本数据 (P1n,P2n,L,Pin)其中,n=1,2,...,N,i表示间歇性分布式电源的个数;
(b)基于指示函数I(t)来计算满足约束条件的样本数,即当抽取的样本值满足约束条件时,I(t)=1,否则,I(t)=0;I(t)表达式如下:
其中,P1、P2分别为两个样本概率,分别为两个样本设定的约束概率;
(c)计算所抽取的样本满足约束条件的概率,并将该概率记为PN,则:
(d)判断PN是否满足PN≥α,若满足,则继续下一步;否则,返回(a)进行重新抽样计算;
(e)至此便可将原来不确定性的概率约束条件转换为确定性的约束条件:
有益效果:与现有技术相比,本发明的方法具有以下优点:首先,在对间歇性的风电、光伏发电等分布式电源相关性分析的基础上,既考虑到了不确定性电源的分布式电特点,又避免抽样样本的重复性,利用分布式随机非重叠抽样的方法将随机约束条件转化成了确定性条件,从而将虚拟电厂分层随机优化模型转换为易于求解的分层确定性优化模型;其次,对间歇性分布式电源的历史数据进行经验分布估计,提高了预测精度;再次,考虑到了间歇性的风电、光伏发电等分布式电源之间的相关性,并运用Copula 理论对其相关性进行了分析;然后,考虑到了虚拟发电厂所含的随机性出力,并根据上下层的优化目标建立了含随机变量的分层随机优化模型。
附图说明
图1是本发明方法的流程图;
图2是本发明方法中的抽样方法流程图。
具体实施方式
下面结合附图对本发明的技术方案进行详细的描述。
如图1所示,本发明的一种多微电网组合的虚拟发电厂分布式随机非重叠抽样的无中心优化方法,包括以下步骤:
(1)考虑间歇性分布式电源的影响,建立分布式控制的虚拟发电厂分层随机优化模型。
(11)以虚拟发电厂全局与主电网之间的交换功率总成本最小化为目标,以电能交换约束和传输线路限制为约束条件,建立上层优化模型。
上层目标:
其中,TC表示虚拟发电厂全局与主电网之间的交换功率的总成本,IMG是虚拟发电厂分布式调度控制的微网数量,ρ是每个微网与主电网之间的能源交换价格,在本实施例中,ρ被认为是统一的电力市场价格,是虚拟发电厂内部微网i与主电网之间的交流电量,表示主电网将电能转移到微网i,否则,表示微网i将电能转移到主电网。
上层约束:
1)电能交换的上限约束和下限约束
其中,i=1,2,...,IMG,和分别代表电网和虚拟发电厂之间的电量交换的下限和上限。
2)传输线路限制约束
其中,m=1,2,...M,M为最大线路数,ηi-m是电力从微网注入到线路的灵敏度,Tm表示线路的能量流限制;
(12)以虚拟发电厂内各微网发电成本最小为目标,以各发电单元处理为约束,同时考虑间歇性分布式电源发电出力的不确定性,建立下层随机优化模型。
下层目标:
其中,j=1,2,...,Jmt,Ji是微网i分布式发电机的数量,是微网i的发电成本,表示微网i中的发电单元j的输出功率,表示微网i中的发电单元j的单位成本,则是微网i发电单元j的发电成本函数;是微网i中风力发电的出力,是风电设备的运行和维护成本,是微网i中的光伏发电的输出,是光伏发电设备的运行和维护成本,是微型涡轮机的运行和维护成本,是微型涡轮机的输出,是微型涡轮机的燃油成本,是储能设备的运行和维护成本,和分别是储能的充电能力和放电能力。
下层约束:
7)功率平衡约束
其中,是储能充电的实际利用率,是储能放电的实际利用率,是微网i的负荷。
8)可控分布式电源(微型涡轮机)的上下限输出
其中,j=1,2,...,Jmt,和分别是微网i中可控分布式电源j出力的上下限, Jmt为微网i中可控分布式电源(微型涡轮机)总数。
9)储能约束
SOCi min≤SOC i≤SOCi max
其中,和分别是充电功率的上下限,和分别是微网i中放电功率的上下限;SOCi是微网i中电池的存储状态,分别是电池储能的上限和下限。
4)虚拟电厂内各微电网中风力与光伏等间歇性分布式电源发电受自然条件的影响具有较强的波动性与随机性,且同一地区地理位置较近的风力发电单元与光伏发电单元、风力发电单元之间及光伏发电单元之间的输出功率具有一定的相关性,故在此采用机会约束的形式表示其运行约束,第i个微网内间歇性分布式电源发电的概率约束表达式如下所示:
其中,为微网i中可控分布式电源k的功率出力,它满足Copula联合概率分布函数;是虚拟发电厂i中可控分布式电源k出力的上下限;Ki为微网i中可控分布式电源总数;α为机会约束条件的置信水平,该值是预先指定己知的。
(2)利用Copula理论分析风电、光伏发电等间歇性分布式电源两两之间的相关性。
对风电、光伏发电的历史数据进行经验分布估计,得到各自的经验累积分布函数,根据平方欧式距离最小的方法选择合适的理论Copula函数,从而建立两者之间的联合概率分布来描述其相关性。
同一区域地理位置相近的多个风力发电单元与光伏发电单元、多个风力发电单元间与多个光伏发电单元间受同一风源、光源以及相近温度、湿度等影响,其输出功率间表现出一定的相关关系。当电网接入大批量的风力、光伏发电单元后,其调度运行将受到一定的影响。对风电场与风电场、风电场与光伏发电、光伏发电与光伏发电出力之间的相关性进行分析,可以提高对此类不确定性电源出力的预测,进而有利于虚拟发电厂的优化调度。
基于历史数据和Copula理论的风力、光伏出力相关性分析可分为以下步骤:
对同一地区相邻的风电场或光伏发电进行两两选择,根据风电场出力历史数据,计算不同Copula函数相应的参数估计值。
利用MATLAB软件统计工具箱中的copulafit函数计算不同Copula函数相应的参数估计值如下:
RHOHAT=copulafit('Gaussian',U)
[RHOHAT,nuhat]=copulafit('t',U)
paramhat=copulafit('family',U)
式中:U表示由两个风电场出力随机变量的累积分布值构成的n×2的矩阵,n表示累积分布函数值的个数,其元素的取值范围为[0,1]。当选择Normal-Copula函数或 t-Copula函数时,RHOHAT为线性相关矩阵的参数估计值,nuhat为t-Copula函数的自由度估计值;当选择阿基米德型Copula函数时,family代表Copula函数的种类,包括Gumbel-Copula函数、Clayton-Copula函数和Frank-Copula函数,paramhat参数为阿基米德型Copula函数的参数估计值。
(22)分别根据每种不同的Copula函数计算其与经验Copula函数之间的平方欧氏距离,选择距离最小的Copula函数作为描述风、光发电场出力之间相关性。
不同的Copula函数所具有的分布特征各不相同,使用不同的Copula函数对风、光发电场出力之间的相关性进行建模,所得到的联合概率分布也会有一定的差别。因而为了更好地描述风光出力之间的相关性,选择合适的Copula函数是分析风力、光伏发电出力相关性的关键。评价所选择Copula函数的准确度通常根据Copula函数与经验Copula 函数之间的平方欧氏距离的大小。平方欧氏距离法首先根据风、光发电出力的历史数据构建出风、光发电场出力之间的经验Copula函数;然后根据所得的参数估计值计算各类Copula函数与经验Copula函数之间的平方欧氏距离;最后选择平方欧氏距离最小的 Copula函数作为描述风、光发电出力相关性的最优函数。
其中经验Copula函数定义如下:
式中:u,v∈[0,1];I[·]为示性函数;当F(xi)≤u时,否则当G(yi)≤v时,否则
(3)利用分布式随机非重叠抽样法对风电、光伏发电等间歇性分布式电源相关性分析后的结果进行抽样,从而将虚拟电厂分层随机优化模型转换为分层确定性优化模型来求解。
在风光相关性分析的基础上,利用分布式随机非重叠抽样方法对所选择的Copula函数生成满足风光出力相关性的随机数矩阵进行抽样,如图2所示,主要步骤如下:
(31)对某一间歇性分布式电源进行独立抽样,建立其与其他各间歇性分布式电源的两两联合概率分布,根据平方欧氏距离最小的方法选择Copula函数的类型来作为描述风、光发电出力相关性的最优函数。
计算平方欧氏距离的表达式如下:
式中:F(x)和G(y)分别是两个风、光发电场出力的边缘累积分布函数;
CEi(F(x),G(y))为经验Copula函数上第i个采样值;CTi(F1(x),G(y))为Copula函数上第i个采样值;N为计算点个数。
(32)由所选的理论Copula函数的类型描述风光两两之间的相关性,得出各自的相关性系数。为了避免相关性太大对出力预测精度的影响,通过分布式筛选,在众多风光联合Copula函数中选择Pearson线性相关系数最小的风光联合Copula函数作为抽样;
Pearson线性相关系数表示如下:
其中,Xi、Yi分别为两个对象的样本,
假设风电场i与光伏发电j1、j2、.......、jn的相关系数分别为ρij1、ρij2、.......、ρijn,通过各分布式电源之间的通信连接,可以传递有链接分布式电源之间的信息。由此,为了避免相关性太大对出力预测精度的影响,可以选择与风电场i Pearson线性相关系数最小的光伏发电建立联合概率分布作为抽样的样本。
ρi抽样=min|ρij| j=j1,j2,...,jn
其中,ρi抽样对应最终所选的风电联合概率分布对应的相关系数,ρij表示风电厂i与各个光伏发电所建立的联合概率分布的相关系数。
(33)利用所选择的Copula函数生成满足风电出力相关性的随机数矩阵。
假设针对每个分布式电源出力随机变量生成随机数的个数为N,利用所选择的Copula函数生成满足分布式电源出力相关性的随机数矩阵为DN×K。
其中,di,j表示节点j对节点i的影响权重,di,j∈[0,1)。
(34)用Copularnd()函数对样本进行分布式随机抽样;
D=Copularnd(Copula类型,Copula相关系数,抽样次数);
(35)对抽样结果进行所需处理。
基于分布式随机非重叠抽样将随机约束条件转换成确定性的约束条件,步骤如下:
(a)从Copula联合概率分布式中抽样获得N组间歇性分布式电源的样本数据(P1n,P2n)(n=1,2,...,N)(在此以二元Copula联合概率分布式为例);
(b)基于指示函数I(t)来计算满足约束条件的样本数,即当抽取的样本值满足约束条件时,I(t)=1,否则,I(t)=0。I(t)表达式如下:
其中,P1、P2分别是两个样本概率,分别为两个样本设定的约束概率。
(c)计算所抽取的样本满足约束条件的概率,并将该概率记为PN,则
(d)判断PN是否满足PN≥α,若满足,则继续下一步;否则,返回(a)进行重新抽样计算;
(e)至此便可将原来不确定性的概率约束条件转换为确定性的约束条件:
Claims (7)
1.一种多微电网组合的虚拟发电厂分布式随机非重叠抽样的无中心优化方法,其特征在于,包括以下步骤:
(1)考虑间歇性分布式电源的影响,建立分布式控制的虚拟发电厂分层随机优化模型;
(2)利用Copula理论分析间歇性分布式电源两两之间的相关性;
(3)利用分布式随机非重叠抽样法对间歇性分布式电源相关性分析后的结果进行抽样,从而将虚拟发电厂分层随机优化模型转换为分层确定性优化模型来求解。
2.根据权利要求1所述的一种多微电网组合的虚拟发电厂分布式随机非重叠抽样的无中心优化方法,其特征在于,所述步骤(1)包括:
(11)以虚拟发电厂全局与主电网之间的交换功率总成本最小化为目标,以电能交换约束和传输线路限制为约束条件,建立上层优化模型;
上层目标为:
其中,TC表示虚拟发电厂全局与主电网之间的交换功率的总成本,IMG是虚拟发电厂分布式调度控制的微网数量,ρ是每个微网与主电网之间的能源交换价格,是虚拟发电厂内部微网i与主电网之间的交流电量,表示主电网将电能转移到微网i,否则,表示微网i将电能转移到主电网;
上层约束为:
1)电能交换的上限约束和下限约束
其中,i=1,2,...,IMG,和分别代表主电网和微网i之间的电量交换的下限和上限;
2)传输线路限制约束
其中,m=1,2,...M,M为最大线路数,ηi-m是电力从微网注入到线路的灵敏度,Tm表示线路m的能量流限制;
(12)以虚拟发电厂内各微网发电成本最小为目标,以微网内各发电单元出力为约束,同时考虑间歇性分布式电源发电出力的不确定性,建立下层随机优化模型;
下层目标:
其中,j=1,2,...,Jmt,Ji是微网i分布式发电机的数量,是微网i的发电成本,表示微网i中的发电单元j的输出功率,表示微网i中的发电单元j的单位成本,则是微网i发电单元j的发电成本函数;
下层约束:
1)功率平衡约束
其中,是微网i中风力发电的出力,是微网i中的光伏发电的输出,是微型涡轮机的输出,和分别是储能的充电能力和放电能力,和分别是储能充电的实际利用率和放电的实际利用率,是微网i的负荷;
2)可控分布式电源的上下限输出
其中,j=1,2,...,Jmt,和分别是微网i中可控分布式电源j出力的上下限,Jmt为微网i中可控分布式电源总数;
3)储能约束
SOCi min≤SOCi≤SOCi max
其中,和是充电功率的上下限,和分别是微网i中放电功率的上下限,SOCi是微网i中电池的存储状态,和分别是电池储能的上限和下限;
4)虚拟电厂内各微电网中间歇性分布式电源发电受自然条件的影响具有较强的波动性与随机性,且同一地区地理位置两两发电单元之间的输出功率具有一定的相关性,故在此采用机会约束的形式表示其运行约束,第i个微网内间歇性分布式电源发电的概率约束表达式如下所示:
其中,为微网i中可控分布式电源k的功率出力,它满足Copula联合概率分布函数;是虚拟发电厂i中可控分布式电源k出力的上下限;Ki为微网i中可控分布式电源总数;α为机会约束条件的置信水平,该值是预先指定己知的。
3.根据权利要求2所述的一种多微电网组合的虚拟发电厂分布式随机非重叠抽样的无中心优化方法,其特征在于:所述微网内各发电单元包括可控分布式电源、间歇性分布式电源和储能,所述可控分布式电源为微型涡轮机。
4.根据权利要求3所述的一种多微电网组合的虚拟发电厂分布式随机非重叠抽样的无中心优化方法,其特征在于:所述步骤(12)中间歇性分布式电源包括风力发电和光伏发电,则下层目标具体为:
其中,是微网i中风力发电的出力,是风电设备的运行和维护成本,是微网i中的光伏发电的输出,是光伏发电设备的运行和维护成本,是微型涡轮机的运行和维护成本,是微型涡轮机的输出,是微型涡轮机的燃油成本,是储能设备的运行和维护成本,和分别是储能的充电能力和放电能力。
5.根据权利要求1所述的一种多微电网组合的虚拟发电厂分布式随机非重叠抽样的无中心优化方法,其特征在于,所述步骤(2)中对各间歇性分布式电源的历史数据进行经验分布估计,得到各自的经验累积分布函数,根据平方欧式距离最小的方法选择合适的理论Copula函数,从而建立各间歇性分布式电源两两之间的联合概率分布来描述其相关性;其相关性分析包括以下步骤:
(21)对同一地区相邻的间歇性分布式电源进行两两选择,根据某一间歇性分布式电源出力历史数据,计算不同Copula函数相应的参数估计值;
(22)分别根据每种不同的Copula函数计算其与经验Copula函数之间的平方欧氏距离,选择距离最小的Copula函数作为描述两两间歇性分布式电源出力之间相关性;
其中经验Copula函数定义如下:
式中:u,v∈[0,1],I[·]为示性函数,当F(xi)≤u时,否则当G(yi)≤v时,否则
6.根据权利要求1所述的一种多微电网组合的虚拟发电厂分布式随机非重叠抽样的无中心优化方法,其特征在于,所述步骤(3)中在间歇性分布式电源相关性分析的基础上,利用分布式随机非重叠抽样方法对所选择的Copula函数生成满足间歇性分布式电源出力相关性的随机数矩阵进行抽样,步骤如下:
(31)对某一间歇性分布式电源进行独立抽样,建立其与其他各间歇性分布式电源的两两联合概率分布,根据平方欧氏距离最小的方法选择Copula函数的类型来作为描述间歇性分布式电源出力相关性的最优函数;
计算平方欧氏距离的表达式如下:
式中:F(x)和G(y)分别是两个间歇性分布式电源出力的边缘累积分布函数;CEi(F(x),G(y))为经验Copula函数上第i个采样值;CTi(F1(x),G(y))为Copula函数上第i个采样值;N为计算点个数;
(32)由所选的理论Copula函数的类型描述间歇性分布式电源两两之间的相关性,得出各自的相关性系数;为了避免相关性太大对出力预测精度的影响,通过分布式筛选,在众多间歇性分布式电源联合Copula函数中选择Pearson线性相关系数最小的间歇性分布式电源联合Copula函数作为抽样;
Pearson线性相关系数表示如下:
其中,X、Y为两个对象的样本,
假设某一间歇性分布式电源i与间歇性分布式电源j1、j2、.......、jn的相关系数分别为ρij1、ρij2、.......、ρijn,通过各间歇性分布式电源之间的通信连接,可以传递有链接分布式电源之间的信息;由此,为了避免相关性太大对出力预测精度的影响,可以选择与间歇性分布式电源i Pearson线性相关系数最小的另一间歇性分布式电源建立联合概率分布作为抽样的样本;
ρi抽样=min|ρij| j=j1,j2,...,jn
其中,ρi抽样对应最终所选的两个间歇性分布式电源联合概率分布对应的相关系数,ρij表示间歇性分布式电源i与其它各个间歇性分布式电源所建立的联合概率分布的相关系数;
(33)利用所选择的Copula函数生成满足间歇性分布式电源出力相关性的随机数矩阵;
假设针对每个间歇性分布式电源出力随机变量生成随机数的个数为N,利用所选择的Copula函数生成满足间歇性分布式电源出力相关性的随机数矩阵为DN×K;
其中,di,j表示节点j对节点i的影响权重,di,j∈[0,1);
(34)用Copularnd()函数对样本进行分布式随机抽样;
D=Copularnd(Copula类型,Copula相关系数,抽样次数)
(35)对抽样结果进行所需处理。
7.根据权利要求6所述的一种多微电网组合的虚拟发电厂分布式随机非重叠抽样的无中心优化方法,其特征在于,基于分布式随机非重叠抽样将随机约束条件转换成确定性的约束条件,步骤如下:
(a)从Copula联合概率分布式中抽样获得N组间歇性分布式电源的样本数据(P1n,P2n,L,Pin)其中,n=1,2,...,N,i表示间歇性分布式电源的个数;
(b)基于指示函数I(t)来计算满足约束条件的样本数,即当抽取的样本值满足约束条件时,I(t)=1,否则,I(t)=0;I(t)表达式如下:
其中,P1、P2分别为两个样本概率,P1 max、分别为两个样本设定的约束概率;
(c)计算所抽取的样本满足约束条件的概率,并将该概率记为PN,则:
(d)判断PN是否满足PN≥α,若满足,则继续下一步;否则,返回(a)进行重新抽样计算;
(e)至此便可将原来不确定性的概率约束条件转换为确定性的约束条件:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710976876.3A CN107749638B (zh) | 2017-10-19 | 2017-10-19 | 多微电网组合的虚拟发电厂分布式随机非重叠抽样的无中心优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710976876.3A CN107749638B (zh) | 2017-10-19 | 2017-10-19 | 多微电网组合的虚拟发电厂分布式随机非重叠抽样的无中心优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107749638A true CN107749638A (zh) | 2018-03-02 |
CN107749638B CN107749638B (zh) | 2021-02-02 |
Family
ID=61252966
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710976876.3A Active CN107749638B (zh) | 2017-10-19 | 2017-10-19 | 多微电网组合的虚拟发电厂分布式随机非重叠抽样的无中心优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107749638B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108763627A (zh) * | 2018-04-13 | 2018-11-06 | 西北工业大学 | 结构机构失效概率灵敏度分解方法、计算方法及应用 |
CN110729721A (zh) * | 2019-10-22 | 2020-01-24 | 国网江西省电力有限公司经济技术研究院 | 电力系统全局备用容量计算方法 |
CN111858827A (zh) * | 2020-07-31 | 2020-10-30 | 成都新潮传媒集团有限公司 | 一种地图点位抽稀展示方法、装置及计算机设备 |
CN112734285A (zh) * | 2021-01-22 | 2021-04-30 | 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 | 一种电站送风机超温预警方法及系统 |
CN113469412A (zh) * | 2021-06-02 | 2021-10-01 | 国核电力规划设计研究院有限公司 | 一种综合能源系统实时运行策略优化方法及系统 |
CN113962477A (zh) * | 2021-11-10 | 2022-01-21 | 国网浙江省电力有限公司杭州供电公司 | 一种产业电量关联聚集预测方法、装置、设备及存储介质 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103903066A (zh) * | 2014-04-04 | 2014-07-02 | 东南大学 | 一种虚拟电厂分层随机优化调度方法 |
CN104065060A (zh) * | 2014-06-09 | 2014-09-24 | 徐多 | 独立微电网系统双层经济调度优化方法 |
CN104779611A (zh) * | 2015-03-23 | 2015-07-15 | 南京邮电大学 | 基于集中式和分布式双层优化策略的微电网经济调度方法 |
CN105631599A (zh) * | 2015-12-30 | 2016-06-01 | 国网甘肃省电力公司电力科学研究院 | 一种虚拟电厂多目标运行调度方法 |
CN105790258A (zh) * | 2016-03-15 | 2016-07-20 | 国电南瑞科技股份有限公司 | 基于正态Copula函数的拉丁超立方抽样法概率潮流计算方法 |
CN106602584A (zh) * | 2017-02-06 | 2017-04-26 | 上海电力设计院有限公司 | 一种基于双层优化模型的多能互补微网储能优化配置方法 |
CN107038292A (zh) * | 2017-04-01 | 2017-08-11 | 三峡大学 | 一种基于自适应多变量非参数核密度估计的多风电场出力相关性建模方法 |
WO2017145664A1 (ja) * | 2016-02-26 | 2017-08-31 | 日本電気株式会社 | 最適化システム、最適化方法および最適化プログラム |
-
2017
- 2017-10-19 CN CN201710976876.3A patent/CN107749638B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103903066A (zh) * | 2014-04-04 | 2014-07-02 | 东南大学 | 一种虚拟电厂分层随机优化调度方法 |
CN104065060A (zh) * | 2014-06-09 | 2014-09-24 | 徐多 | 独立微电网系统双层经济调度优化方法 |
CN104779611A (zh) * | 2015-03-23 | 2015-07-15 | 南京邮电大学 | 基于集中式和分布式双层优化策略的微电网经济调度方法 |
CN105631599A (zh) * | 2015-12-30 | 2016-06-01 | 国网甘肃省电力公司电力科学研究院 | 一种虚拟电厂多目标运行调度方法 |
WO2017145664A1 (ja) * | 2016-02-26 | 2017-08-31 | 日本電気株式会社 | 最適化システム、最適化方法および最適化プログラム |
CN105790258A (zh) * | 2016-03-15 | 2016-07-20 | 国电南瑞科技股份有限公司 | 基于正态Copula函数的拉丁超立方抽样法概率潮流计算方法 |
CN106602584A (zh) * | 2017-02-06 | 2017-04-26 | 上海电力设计院有限公司 | 一种基于双层优化模型的多能互补微网储能优化配置方法 |
CN107038292A (zh) * | 2017-04-01 | 2017-08-11 | 三峡大学 | 一种基于自适应多变量非参数核密度估计的多风电场出力相关性建模方法 |
Non-Patent Citations (2)
Title |
---|
SUI PENG等: "Probabilistic Power Flow for AC/VSC-MTDC Hybrid Grids Considering Rank Correlation Among Diverse Uncertainty Sources", 《 IEEE TRANSACTIONS ON POWER SYSTEMS 》 * |
冯其芝: "含微电网的虚拟电厂分层随机优化调度", 《中国优秀硕士学位论文全文数据库》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108763627A (zh) * | 2018-04-13 | 2018-11-06 | 西北工业大学 | 结构机构失效概率灵敏度分解方法、计算方法及应用 |
CN110729721A (zh) * | 2019-10-22 | 2020-01-24 | 国网江西省电力有限公司经济技术研究院 | 电力系统全局备用容量计算方法 |
CN110729721B (zh) * | 2019-10-22 | 2023-07-25 | 国网江西省电力有限公司经济技术研究院 | 电力系统全局备用容量计算方法 |
CN111858827A (zh) * | 2020-07-31 | 2020-10-30 | 成都新潮传媒集团有限公司 | 一种地图点位抽稀展示方法、装置及计算机设备 |
CN111858827B (zh) * | 2020-07-31 | 2023-11-14 | 成都新潮传媒集团有限公司 | 一种地图点位抽稀展示方法、装置及计算机设备 |
CN112734285A (zh) * | 2021-01-22 | 2021-04-30 | 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 | 一种电站送风机超温预警方法及系统 |
CN112734285B (zh) * | 2021-01-22 | 2023-05-02 | 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 | 一种电站送风机超温预警方法及系统 |
CN113469412A (zh) * | 2021-06-02 | 2021-10-01 | 国核电力规划设计研究院有限公司 | 一种综合能源系统实时运行策略优化方法及系统 |
CN113469412B (zh) * | 2021-06-02 | 2024-04-09 | 国核电力规划设计研究院有限公司 | 一种综合能源系统实时运行策略优化方法及系统 |
CN113962477A (zh) * | 2021-11-10 | 2022-01-21 | 国网浙江省电力有限公司杭州供电公司 | 一种产业电量关联聚集预测方法、装置、设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN107749638B (zh) | 2021-02-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107749638A (zh) | 多微电网组合的虚拟发电厂分布式随机非重叠抽样的无中心优化方法 | |
Dinh et al. | A home energy management system with renewable energy and energy storage utilizing main grid and electricity selling | |
WO2022048127A1 (zh) | 一种电热泵-热电联合系统的优化调控方法及系统 | |
Gonzalez et al. | Environmental and cost optimal design of a biomass–Wind–PV electricity generation system | |
de Queiroz | Stochastic hydro-thermal scheduling optimization: An overview | |
CN104779611B (zh) | 基于集中式和分布式双层优化策略的微电网经济调度方法 | |
Amjady et al. | Short-term load forecast of microgrids by a new bilevel prediction strategy | |
Aluisio et al. | An optimization procedure for Microgrid day-ahead operation in the presence of CHP facilities | |
Mandal et al. | Short-term combined economic emission scheduling of hydrothermal systems with cascaded reservoirs using particle swarm optimization technique | |
Shukla et al. | Clustering based unit commitment with wind power uncertainty | |
CN108734350A (zh) | 一种含微电网的配电网独立与联合调度的求解方法 | |
JP2022550225A (ja) | 電気ヒートポンプ-熱電併給システムを最適化する調整方法及びシステム | |
Sharma et al. | A critical and comparative review of energy management strategies for microgrids | |
CN107423852A (zh) | 一种计及典型场景的光储联合电站优化运营方法 | |
CN109508499A (zh) | 多时段多场景分布式电源最优接入位置与容量研究方法 | |
CN105703368B (zh) | 一种集成主动配电网和输电网多重不确定能量流建模 | |
Kunya et al. | Review of economic dispatch in multi-area power system: State-of-the-art and future prospective | |
CN103580020A (zh) | 一种基于NSGA-II和Look-ahead的含风电场电力系统多目标动态优化调度方法 | |
Forouzandehmehr et al. | Stochastic dynamic game between hydropower plant and thermal power plant in smart grid networks | |
Chen et al. | A battery management strategy in microgrid for personalized customer requirements | |
Dagdougui et al. | Global energy management system for cooperative networked residential green buildings | |
CN103715688A (zh) | 并网型微网优化配置方法及装置 | |
CN107681655A (zh) | 一种潮汐流能发电场协调规划方法 | |
CN110400056A (zh) | 基于风光时空特性的梯级水电日前优化调度方法及装置 | |
CN108667077A (zh) | 一种风蓄联合系统优化调度方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |