CN109038648A - 一种基于Copula函数的风光联合出力建模方法 - Google Patents
一种基于Copula函数的风光联合出力建模方法 Download PDFInfo
- Publication number
- CN109038648A CN109038648A CN201810754220.1A CN201810754220A CN109038648A CN 109038648 A CN109038648 A CN 109038648A CN 201810754220 A CN201810754220 A CN 201810754220A CN 109038648 A CN109038648 A CN 109038648A
- Authority
- CN
- China
- Prior art keywords
- power output
- copula function
- copula
- function
- parameter
- 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
- 241000039077 Copula Species 0.000 title claims abstract description 113
- 238000000034 method Methods 0.000 title claims abstract description 57
- 238000007476 Maximum Likelihood Methods 0.000 claims abstract description 36
- 238000005315 distribution function Methods 0.000 claims abstract description 32
- 230000003068 static effect Effects 0.000 claims abstract description 31
- 239000011159 matrix material Substances 0.000 claims description 15
- 241000208340 Araliaceae Species 0.000 claims description 4
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims description 4
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 4
- 235000008434 ginseng Nutrition 0.000 claims description 4
- 238000005553 drilling Methods 0.000 claims description 2
- 230000005684 electric field Effects 0.000 claims description 2
- 230000007423 decrease Effects 0.000 abstract 1
- 230000005611 electricity Effects 0.000 description 9
- 238000010586 diagram Methods 0.000 description 6
- 230000007547 defect Effects 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 241000196324 Embryophyta Species 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 241001123248 Arma Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- H02J3/383—
-
- 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
- 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
- Photovoltaic Devices (AREA)
Abstract
本发明公开了一种基于Copula函数的风光联合出力建模方法,包括:将风电场和光伏电场的输出功率转化为风电出力序列和光伏出力序列后进行核密度估计,得到风电出力和光伏出力的边缘分布函数;将风电出力和光伏出力的边缘分布函数代入Copula函数得到Copula函数的静态相关系数,利用Copula函数的静态相关系数建立风光联合出力的静态模型;将风电出力序列和光伏出力序列代入演进方程,利用极大似然估计法得到Copula函数的动态相关系数,利用Copula函数的动态相关系数建立风光联合出力的动态模型;从风光联合出力的静态模型和动态模型中获取最优模型。本发明减小了估计误差,使模型拟合优度更高。
Description
技术领域
本发明属于新能源出力不确定性建模领域,更具体地,涉及一种基于Copula函数的风光联合出力建模方法。
背景技术
近年来,我国新能源装机容量逐年增加,装机总量居世界第一,消纳总量实现了快速增长。随着新能源渗透率的不断提高,消纳矛盾也逐渐凸显,如何减少弃风、弃光,成为学术研究的焦点问题。而风电场、光伏电站等间歇式能源动态特性建模验证工作则是解决此问题的基础。
为建立风光联合出力动态Copula模型,首先需要建立风电、光伏电场出力的边缘分布,目前常见建立边缘分布的方法分为两大类:一类是先拟合风速、光强的分布,然后将风速、光强的分布转化为输出功率的分布;另一类是直接通过核密度估计方法拟合风电、光伏输出功率的分布。其中,基于利用风速、光强分布建立边缘分布的方法直观且容易理解。但是,该方法建立的边缘分布无法考虑弃风、弃光的情况,且在后续估计环节会引入误差。此外,长期连续的风速、光强数据较风电、光伏输出功率数据而言更难以获取。
Copula函数是一维边缘分布与多维联合分布之间的连接函数,是构建多维随机变量联合概率分布的有效工具。确定Copula函数分布特性的关键在于相关系数的估计,而联合分布相关系数的求取依赖于边缘分布函数的形式。根据边缘分布函数的选取方法的不同,常用的几种Copula函数相关系数估计方法有两阶段估计法(IFM)、基于经验分布函数的极大似然估计法(CML)和基于非参数核密度的极大似然估计法(MLK)。其中,MLK属于半参数估计法,不需要估计边缘分布中的未知参数,因此估计效果不依赖边缘分布的拟合优度。在风光出力这类数据总量较大、无法获得确切边缘分布函数的情况下,选用MLK进行参数估计的误差最小,效果最为理想。此外,针对不同类型的数据,适用于建模的Copula函数是不同的,采用单一Copula函数建立模型往往无法做到拟合优度最高。
静态Copula模型面向数据总体相关性分析有明显优势,但无法刻画不同时段相关性的变化,无法描述相关性的动态特性,在面对风电、光伏此类相关性随时间发生本质变化的数据时无法较好拟合。
由此可见,现有技术在描述风光联合出力时均存在缺陷,存在拟合优度不高的问题。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提供了一种基于Copula函数的风光联合出力建模方法,由此解决现有技术在描述风光联合出力时均存在缺陷,存在拟合优度不高的技术问题。
为实现上述目的,本发明提供了一种基于Copula函数的风光联合出力建模方法,包括:
(1)将同一地区的风电场和光伏电场的输出功率转化为风电出力序列和光伏出力序列后进行核密度估计,得到风电出力和光伏出力的边缘分布函数;
(2)将风电出力和光伏出力的边缘分布函数代入Copula函数,利用极大似然估计法得到Copula函数的静态相关系数,利用Copula函数的静态相关系数建立基于Copula函数的风光联合出力的静态模型;
(3)将风电出力序列和光伏出力序列代入演进方程,利用极大似然估计法得到Copula函数的动态相关系数,利用Copula函数的动态相关系数建立基于Copula函数的风光联合出力的动态模型;
(4)根据赤池信息准则和贝叶斯信息准则从基于Copula函数的风光联合出力的静态模型和基于Copula函数的风光联合出力的动态模型中获取最优模型。
进一步地,Copula函数包括:N-Copula函数、t-Copula函数、Clayton-Copula函数和SJC-Copula函数。
进一步地,步骤(1)包括:
获取同一地区的风电场和光伏电场的输出功率,以风电场的额定出力为第一标幺值对风电场的输出功率标幺化,得到风电出力序列,以光伏电场的额定出力为第二标幺值对光伏电场的输出功率标幺化,得到光伏出力序列,对风电出力序列和光伏出力序列进行核密度估计,得到风电出力和光伏出力的边缘分布函数。
进一步地,步骤(2)还包括:利用极大似然估计法得到基于Copula函数的风光联合出力的静态模型的对数极大似然估计值。
进一步地,演进方程包括:N-Copula函数的演进方程、t-Copula函数的演进方程、Clayton-Copula函数的演进方程和SJC-Copula函数的演进方程;
所述N-Copula函数和t-Copula函数的演进方程为:
Qt=R(1-α-β)+α(εt-1ε′t-1)+βQt-1
所述Clayton-Copula函数的演进方程为:
所述SJC-Copula函数的演进方程为:
其中,Rt为t时刻的N-Copula函数和t-Copula函数的动态相关系数矩阵,Qt为t时刻的演进相关系数矩阵,R、β和α分别为N-Copula函数和t-Copula函数的演进方程的第一参数、第二参数和第三参数,εt-1为(t-1)时刻风电出力和光伏出力的边缘分布函数的伪逆函数序列矩阵,ε′t-1为εt-1的转置,Qt-1为(t-1)时刻的演进相关系数矩阵,θC.t为t时刻的Clayton-Copula函数的动态相关系数,ωC、βC和αC分别为Clayton-Copula函数的演进方程的第一参数、第二参数和第三参数,wt-j为(t-j)时刻的风电出力序列,st-j为(t-j)时刻的光伏出力序列,θt-1为(t-1)时刻的Clayton-Copula函数的动态相关系数,和分别为t时刻的SJC-Copula函数的上尾相关系数和下尾相关系数,和分别为(t-1)时刻的SJC-Copula函数的上尾相关系数和下尾相关系数,ωU、βU和αU分别为SJC-Copula函数的演进方程的第一上尾参数、第二上尾参数和第三上尾参数,ωL、βL和αL分别为SJC-Copula函数的演进方程的第一下尾参数、第二下尾参数和第三下尾参数,
进一步地,步骤(3)包括:
将风电出力序列和光伏出力序列代入N-Copula函数的演进方程,利用极大似然估计法得到N-Copula函数的演进方程的第一参数、第二参数和第三参数,将风电出力序列和光伏出力序列代入t-Copula函数的演进方程,利用极大似然估计法得到t-Copula函数的演进方程的第一参数、第二参数和第三参数,将风电出力序列和光伏出力序列代入Clayton-Copula函数的演进方程,利用极大似然估计法得到Clayton-Copula函数的演进方程的第一参数、第二参数和第三参数,将风电出力序列和光伏出力序列代入SJC-Copula函数的演进方程,利用极大似然估计法得到SJC-Copula函数的演进方程的第一上尾参数、第二上尾参数和第三上尾参数,SJC-Copula函数的演进方程的第一下尾参数、第二下尾参数和第三下尾参数,然后将风电出力序列和光伏出力序列分别代入N-Copula函数的演进方程、t-Copula函数的演进方程、Clayton-Copula函数的演进方程和SJC-Copula函数的演进方程,得到Copula函数的动态相关系数,包括:N-Copula函数的动态相关系数矩阵、t-Copula函数的动态相关系数矩阵、Clayton-Copula函数的动态相关系数、SJC-Copula函数的上尾相关系数和下尾相关系数,利用Copula函数的动态相关系数建立基于Copula函数的风光联合出力的动态模型。
进一步地,步骤(3)还包括:利用极大似然估计法得到基于Copula函数的风光联合出力的动态模型的对数极大似然估计值。
进一步地,步骤(4)包括:
将基于Copula函数的风光联合出力的静态模型的对数极大似然估计值和基于Copula函数的风光联合出力的动态模型的对数极大似然估计值代入赤池信息准则和贝叶斯信息准则,得到每个模型的拟合优度,将拟合优度最大值对应的模型作为最优模型。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,能够取得下列有益效果:
(1)静态Copula模型面向数据总体相关性分析有明显优势,但无法刻画不同时段相关性的变化,无法描述相关性的动态特性,在面对风电、光伏此类相关性随时间发生本质变化的数据时无法较好拟合。本发明将风电出力序列和光伏出力序列代入演进方程,利用极大似然估计法得到Copula函数的动态相关系数,利用Copula函数的动态相关系数建立基于Copula函数的风光联合出力的动态模型;通过引入动态Copula函数,刻画了风光联合出力的动态相关性,减小了估计误差,使模型拟合优度更高。
(2)本发明直接采用风电出力序列和光伏出力序列作为演进方程的输入,能够考虑弃风、弃光的情况。利用极大似然估计法得到Copula函数的静态相关系数、Copula函数的动态相关系数以及演进方程中的参数,减小了估计误差,使模型拟合优度更高。使得本发明的方法计算量小且拟合优度高,能有效处理风光联合出力的建模问题。
附图说明
图1是本发明实施例提供的一种基于Copula函数的风光联合出力建模方法的流程图;
图2是本发明实施例1提供的风电场出力统计直方图及核密度估计曲线示意图;
图3是本发明实施例1提供的光伏电场出力统计直方图及核密度估计曲线示意图;
图4是本发明实施例1提供的风光联合出力静态N-Copula模型示意图;
图5是本发明实施例1提供的风光联合出力静态t-Copula模型示意图;
图6是本发明实施例1提供的风光联合出力静态Clayton-Copula模型示意图;
图7是本发明实施例1提供的风光联合出力静态SJC-Copula模型示意图;
图8是本发明实施例1提供的风光联合出力动态N-Copula模型相关系数示意图;
图9是本发明实施例1提供的风光联合出力动态t-Copula模型相关系数示意图;
图10是本发明实施例1提供的风光联合出力动态Clayton-Copula模型相关系数示意图;
图11是本发明实施例提供的风光联合出力动态SJC-Copula模型上尾相关系数示意图;
图12是本发明实施例提供的风光联合出力动态SJC-Copula模型下尾相关系数示意图;
图13本发明实施例提供的基于赤池信息准则进行拟合优度检验对比示意图;
图14本发明实施例提供的基于贝叶斯信息准则进行拟合优度检验对比示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
如图1所示,一种基于Copula函数的风光联合出力建模方法,包括:
(1)获取同一地区的风电场和光伏电场的输出功率,以风电场的额定出力为第一标幺值对风电场的输出功率标幺化,得到风电出力序列,以光伏电场的额定出力为第二标幺值对光伏电场的输出功率标幺化,得到光伏出力序列,对风电出力序列和光伏出力序列进行核密度估计,得到风电出力和光伏出力的边缘分布函数。
(2)将风电出力和光伏出力的边缘分布函数代入Copula函数,利用极大似然估计法得到Copula函数的静态相关系数,利用Copula函数的静态相关系数建立基于Copula函数的风光联合出力的静态模型;
(3)将风电出力序列和光伏出力序列代入演进方程,利用极大似然估计法得到Copula函数的动态相关系数,利用Copula函数的动态相关系数建立基于Copula函数的风光联合出力的动态模型;
(4)根据赤池信息准则和贝叶斯信息准则从基于Copula函数的风光联合出力的静态模型和基于Copula函数的风光联合出力的动态模型中获取最优模型。
具体地,本发明基于非参数核密度估计法拟合出力,建立风电出力和光伏出力的边缘分布函数,其中非参数核密度公式为:
其中,Kh为核函数,N为序列容量,Xi为风速样本,i=1,2,...,N,v为风电出力序列或者光伏出力序列。
其次,本发明采用分别用Copula函数建立风光联合出力模型。Copula函数专门用来刻画变量间的非线性相关性,动态Copula函数还能描述此类相关性的动态特性。
基于Copula函数的联合分布函数和联合密度函数表达式分别为:
F(x1,x2,…,xn)=C(F(x1),F(x2),…,F(xn)) (2)
其中,F(xn)和f(xn)分别为变量边缘分布的分布函数和密度函数。
Copula分布函数C(F(x1),F(x2),…,F(xn))与Copula密度函数c(f(x1),f(x2),…,f(xn))的关系为:
由式(2)(3)可知,联合分布函数即为选取的Copula函数,而联合密度函数即为选取的密度函数与各边缘分布函数的积。
其次,本发明采用N-Copula、t-Copula、Clayton-Copula、SJC-Copula函数建立多种风光联合出力动态模型。表达式如下
1)Clayton-Copula:
Clayton-Copula函数在二维情况下表达式为:
其中u,v为边缘分布函数,θC.t为动态相关系数。
2)SJC-Copula:
SJC-Copula函数由Joe-Clayton Copula(JC-Copula)函数变换而来,表征SJC-Copula函数和JC-Copula函数变换关系的函数表达式如下:
二维情况下Joe-Clayton Copula函数表达式为:
其中,u,v为边缘分布函数,分别为动态上、下尾相关系数。
3)N-Copula和t-Copula
N-Copula、t-CoDula函数结构类似,在二维情况下表达式分别为:
其中u,v为边缘分布函数,ρN,t、ρt,t为N-Copula函数和t-Copula函数的动态相关系数,kt为自由度。
其次,本发明采用演进方程使Copula函数动态化。在Copula模型中,对相关系数的估计转化为对演进方程参数的估计。利用非参数核密度估计方法,将风电出力序列wi、si代入Clayton-Copula函数的演进方程,其表达式为:
其中,ωC、βC和αC分别为Clayton-Copula函数的演进方程的第一参数、第二参数和第三参数,为logistic函数,引入此函数保证θC.t∈(-1,1)。
演进方程将静态的变为动态的,而对应的似然函数也从相关系数的函数变为演进方程系数ωC、βC和αC的函数。求解该似然函数即可求得对数极大似然估计值LogL和对应的ωC、βC和αC。将ωC、βC和αC及wi、si代入演进方程,即可求得动态相关系数序列θC.t。
Clayton-Copula模型的演进特性类似于有限制条件的ARMA(1,10)模型,参数βC控制相关系数的递进关系,参数αC控制前10组数据对相关系数的影响。Clayton-Copula函数具有不对称的尾部特征,下尾高,上尾低,对联合出力在分布下尾处的变化比较敏感,即对风光出力较小时的相关特性能进行有效刻画。
SJC-Copula函数的分布特性由上尾相关系数和下尾相关系数刻画,其演进方程分别为:
其中,ωU、βU、αU和ωL、βL、αL分别为SJC-Copula函数的上尾参数和下尾参数。是logistic函数,保证
SJC-Copula函数具有不对称的尾部特征,对联合出力上尾、下尾的变化均敏感,即对风光出力较小时和较大时的相关特性均能进行有效刻画。
N-Copula函数的分布特性由相关系数ρN,t刻画,动态t-Copula函数的分布特性由相关系数ρt,t和自由度k刻画。N-Copula和t-Copula的相关系数矩阵可依DCC(1,1)分解为
其中,Qt和的演化方程为:
Qt=R(1-α-β)+α(εt-1ε′t-1)+βQt-1 (14)
R、β和α分别为N-Copula函数和t-Copula函数的演进方程的第一参数、第二参数和第三参数,εt-1为(t-1)时刻风电出力和光伏出力的边缘分布函数的伪逆函数序列矩阵,ε′t-1为εt-1的转置,Qt-1为(t-1)时刻的演进相关系数矩阵,且必须满足0<α<1,0<β<1,0<α+β<1,同时Ωt必须是正定的,以保证Rt∈(-1,1)。
这两种Copula分布都具有对称的尾部特征,无法描述随机变量之间非对称的尾部相依关系。其中,二维动态t-Copula的尾部较厚,因此对对称的尾部关系可以进行较好的刻画。
其次,本发明采用赤池信息准则(AIC)和贝叶斯信息准则(BIC)选取最优模型。
AIC和BIC的定义式如下:
AIC=2k-2ln(L) (15)
BIC=kln(n)-2ln(L) (16)
其中,k为模型参数个数,L为极大似然估计值,n为样本数量。结合定义可知,AIC和BIC是体现模型复杂度和拟合优劣的综合指标,其数值越小表示模型拟合优度越高,越简洁。
实施例1
一种基于Copula函数的风光联合出力建模方法,包括:
(1)基于非参数核密度估计方法,对风电出力、光伏出力序列标幺化并进行核密度估计,建立风电出力和光伏出力的边缘分布函数。
具体地,在本实施例中,使用湖北随州一光伏电场和一风电场输出功率实测数据,构建风光出力的边缘分布。数据选取2016年6月全月30天的同步测量数据Wi、Si,由于光伏电站在夜间无出力,故提取上午7点至傍晚7点的出力数据进行分析,每5分钟取一个测量点,则每天取144点,30天共计4320组数据。光伏电站的额定出力为100MW,风电场的额定出力为148.5MW。以额定出力为标幺值对测量数据标幺化,得到标幺化的风电、光伏出力序列wi、si:
wi=Wi/100 (17)
si=Si/148.5 (18)
利用式(1)对标幺化的风电、光伏出力序列wi、si作非参数核密度估计,得到风电出力的边缘分布函数fw(x1)、光伏出力边缘分布函数fs(x2),如图2、图3所示,表达式为:
(2)基于静态copula理论,建立基于N-Copula、t-Copula、Clayton-Copula、SJC-Copula函数的风光联合出力的静态模型。
具体的,在本实施例中,将风电出力的边缘分布函数fw(x1)、光伏出力边缘分布函数fs(x2)分别代入Copula函数表达式式(5)、(6)、(7)、(8)、(9),得到确定Copula函数的静态相关系数,Clayton-Copula函数的相关系数表示为θC,SJC-Copula函数的上尾、下尾相关系数分别表示为τU,τL,N-Copula函数和t-Copula函数的相关系数分别表示为ρN、ρt。用极大似然估计法计算相关系数,即可获得静态Copula模型,如图4、图5、图6、图7所示。
(3)基于动态copula理论,引入演进方程,求出演进方程参数,建立基于不同Copula函数的风光联合出力的动态模型。
具体地,在本实施例中,将风电出力序列和光伏出力序列代入对应的演进方程式(10)、(11)、(12)、(13)中,重新利用极大似然估计法计算演进方程的参数,并将风电出力的边缘分布函数fw(x1)、光伏出力边缘分布函数fs(x2)代入演进方程,计算动态相关系数序列,Clayton-Copula函数的相关系数表示为θC.t,SJC-Copula函数的上尾、下尾相关系数分别表示为N-Copula函数和t-Copula函数的相关系数分别表示为ρN,t、ρt,t,如图8、图9、图10、图11、图12所示。将动态相关系数代入Copula函数表达式(5)、(6)、(7)、(8)、(9),即可获得动态Copula模型。
(4)根据赤池信息准则和贝叶斯信息准则选出了适合输入数据的最优模型。
具体地,在本实施例中,将极大似然估计法求得的似然函数LogL代入式(15)、(16),即可计算不同Copula函数对应的拟合优度。比较拟合优度,选出最优Copula模型。拟合优度对比如图13,图14。由图可知,SJC-Copula模型是最适用于本实施例的Copula模型。本发明针对现有技术中静态Copula模型面向数据总体相关性分析有明显优势,但无法刻画不同时段相关性的变化,无法描述相关性的动态特性,在面对风电、光伏此类相关性随时间发生本质变化的数据时无法较好拟合。本发明将风电出力序列和光伏出力序列代入演进方程,利用极大似然估计法得到Copula函数的动态相关系数,利用Copula函数的动态相关系数建立基于Copula函数的风光联合出力的动态模型;通过引入动态Copula函数,刻画了风光联合出力的动态相关性,减小了估计误差,使模型拟合优度更高。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种基于Copula函数的风光联合出力建模方法,其特征在于,包括:
(1)将同一地区的风电场和光伏电场的输出功率转化为风电出力序列和光伏出力序列后进行核密度估计,得到风电出力和光伏出力的边缘分布函数;
(2)将风电出力和光伏出力的边缘分布函数代入Copula函数,利用极大似然估计法得到Copula函数的静态相关系数,利用Copula函数的静态相关系数建立基于Copula函数的风光联合出力的静态模型;
(3)将风电出力序列和光伏出力序列代入演进方程,利用极大似然估计法得到Copula函数的动态相关系数,利用Copula函数的动态相关系数建立基于Copula函数的风光联合出力的动态模型;
(4)根据赤池信息准则和贝叶斯信息准则从基于Copula函数的风光联合出力的静态模型和基于Copula函数的风光联合出力的动态模型中获取最优模型。
2.如权利要求1所述的一种基于Copula函数的风光联合出力建模方法,其特征在于,所述Copula函数包括:N-Copula函数、t-Copula函数、Clayton-Copula函数和SJC-Copula函数。
3.如权利要求1或2所述的一种基于Copula函数的风光联合出力建模方法,其特征在于,所述步骤(1)包括:
获取同一地区的风电场和光伏电场的输出功率,以风电场的额定出力为第一标幺值对风电场的输出功率标幺化,得到风电出力序列,以光伏电场的额定出力为第二标幺值对光伏电场的输出功率标幺化,得到光伏出力序列,对风电出力序列和光伏出力序列进行核密度估计,得到风电出力和光伏出力的边缘分布函数。
4.如权利要求1或2所述的一种基于Copula函数的风光联合出力建模方法,其特征在于,所述步骤(2)还包括:利用极大似然估计法得到基于Copula函数的风光联合出力的静态模型的对数极大似然估计值。
5.如权利要求3所述的一种基于Copula函数的风光联合出力建模方法,其特征在于,所述演进方程包括:N-Copula函数的演进方程、t-Copula函数的演进方程、Clayton-Copula函数的演进方程和SJC-Copula函数的演进方程;
所述N-Copula函数和t-Copula函数的演进方程为:
Qt=R(1-α-β)+α(εt-1ε′t-1)+βQt-1
所述Clayton-Copula函数的演进方程为:
所述SJC-Copula函数的演进方程为:
其中,Rt为t时刻的N-Copula函数和t-Copula函数的动态相关系数矩阵,Qt为t时刻的演进相关系数矩阵,R、β和α分别为N-Copula函数和t-Copula函数的演进方程的第一参数、第二参数和第三参数,εt-1为(t-1)时刻风电出力和光伏出力的边缘分布函数的伪逆函数序列矩阵,ε′t-1为εt-1的转置,Qt-1为(t-1)时刻的演进相关系数矩阵,θC.t为t时刻的Clayton-Copula函数的动态相关系数,ωC、βC和αC分别为Clayton-Copula函数的演进方程的第一参数、第二参数和第三参数,wt-j为(t-j)时刻的风电出力序列,st-j为(t-j)时刻的光伏出力序列,θt-1为(t-1)时刻的Clayton-Copula函数的动态相关系数, 和分别为t时刻的SJC-Copula函数的上尾相关系数和下尾相关系数,和分别为(t-1)时刻的SJC-Copula函数的上尾相关系数和下尾相关系数,ωU、βU和αU分别为SJC-Copula函数的演进方程的第一上尾参数、第二上尾参数和第三上尾参数,ωL、βL和αL分别为SJC-Copula函数的演进方程的第一下尾参数、第二下尾参数和第三下尾参数,
6.如权利要求5所述的一种基于Copula函数的风光联合出力建模方法,其特征在于,所述步骤(3)包括:
将风电出力序列和光伏出力序列代入N-Copula函数的演进方程,利用极大似然估计法得到N-Copula函数的演进方程的第一参数、第二参数和第三参数,将风电出力序列和光伏出力序列代入t-Copula函数的演进方程,利用极大似然估计法得到t-Copula函数的演进方程的第一参数、第二参数和第三参数,将风电出力序列和光伏出力序列代入Clayton-Copula函数的演进方程,利用极大似然估计法得到Clayton-Copula函数的演进方程的第一参数、第二参数和第三参数,将风电出力序列和光伏出力序列代入SJC-Copula函数的演进方程,利用极大似然估计法得到SJC-Copula函数的演进方程的第一上尾参数、第二上尾参数和第三上尾参数,SJC-Copula函数的演进方程的第一下尾参数、第二下尾参数和第三下尾参数,然后将风电出力序列和光伏出力序列分别代入N-Copula函数的演进方程、t-Copula函数的演进方程、Clayton-Copula函数的演进方程和SJC-Copula函数的演进方程,得到Copula函数的动态相关系数,包括:N-Copula函数的动态相关系数矩阵、t-Copula函数的动态相关系数矩阵、Clayton-Copula函数的动态相关系数、SJC-Copula函数的上尾相关系数和下尾相关系数,利用Copula函数的动态相关系数建立基于Copula函数的风光联合出力的动态模型。
7.如权利要求1或2所述的一种基于Copula函数的风光联合出力建模方法,其特征在于,所述步骤(3)还包括:利用极大似然估计法得到基于Copula函数的风光联合出力的动态模型的对数极大似然估计值。
8.如权利要求1或2所述的一种基于Copula函数的风光联合出力建模方法,其特征在于,所述步骤(4)包括:
将基于Copula函数的风光联合出力的静态模型的对数极大似然估计值和基于Copula函数的风光联合出力的动态模型的对数极大似然估计值代入赤池信息准则和贝叶斯信息准则,得到每个模型的拟合优度,将拟合优度最大值对应的模型作为最优模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810754220.1A CN109038648B (zh) | 2018-07-10 | 2018-07-10 | 一种基于Copula函数的风光联合出力建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810754220.1A CN109038648B (zh) | 2018-07-10 | 2018-07-10 | 一种基于Copula函数的风光联合出力建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109038648A true CN109038648A (zh) | 2018-12-18 |
CN109038648B CN109038648B (zh) | 2020-11-17 |
Family
ID=64640461
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810754220.1A Expired - Fee Related CN109038648B (zh) | 2018-07-10 | 2018-07-10 | 一种基于Copula函数的风光联合出力建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109038648B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109753763A (zh) * | 2019-03-06 | 2019-05-14 | 国网江苏省电力有限公司 | 一种风光联合出力概率建模方法 |
CN110287601A (zh) * | 2019-06-27 | 2019-09-27 | 浙江农林大学 | 一种毛竹胸径年龄二元联合分布精确估算方法 |
CN110414099A (zh) * | 2019-07-15 | 2019-11-05 | 东南大学 | 一种基于Copula的相关结构模型综合选型方法 |
CN110611422A (zh) * | 2019-09-17 | 2019-12-24 | 福州大学 | 一种考虑mmc模块间相关性的子模块冗余优化配置方法及系统 |
CN110705099A (zh) * | 2019-09-30 | 2020-01-17 | 华北电力大学 | 一种检定风电场出力相关性的方法 |
CN111162566A (zh) * | 2020-01-21 | 2020-05-15 | 中国电力科学研究院有限公司 | 一种确定风光联合发电系统的联合出力特性的方法和系统 |
CN111274542A (zh) * | 2020-01-10 | 2020-06-12 | 国网江苏省电力有限公司扬州供电分公司 | 基于藤copula和混合偏移正态分布的置信容量评估方法及装置 |
CN112163695A (zh) * | 2020-09-09 | 2021-01-01 | 华南理工大学 | 基于copula函数的风力光伏发电预测方法、系统和介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102801157A (zh) * | 2012-07-24 | 2012-11-28 | 江苏省电力设计院 | 基于Copula理论的风光互补发电系统可靠性评估方法 |
WO2015085308A1 (en) * | 2013-12-07 | 2015-06-11 | Cardinal Wind, Inc. | Computer-implemented data analysis methods and systems for wind energy assessments |
CN104915788A (zh) * | 2015-06-24 | 2015-09-16 | 华南理工大学 | 一种考虑多风场相关性的电力系统动态经济调度的方法 |
CN107038292A (zh) * | 2017-04-01 | 2017-08-11 | 三峡大学 | 一种基于自适应多变量非参数核密度估计的多风电场出力相关性建模方法 |
CN107834547A (zh) * | 2017-11-15 | 2018-03-23 | 南通大学 | 一种考虑风电场输出功率关联特性的输电网规划方法 |
CN108074038A (zh) * | 2017-12-11 | 2018-05-25 | 国网江苏省电力有限公司经济技术研究院 | 一种考虑可再生能源及负荷多时空分布特性的电力生产模拟方法 |
-
2018
- 2018-07-10 CN CN201810754220.1A patent/CN109038648B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102801157A (zh) * | 2012-07-24 | 2012-11-28 | 江苏省电力设计院 | 基于Copula理论的风光互补发电系统可靠性评估方法 |
WO2015085308A1 (en) * | 2013-12-07 | 2015-06-11 | Cardinal Wind, Inc. | Computer-implemented data analysis methods and systems for wind energy assessments |
CN104915788A (zh) * | 2015-06-24 | 2015-09-16 | 华南理工大学 | 一种考虑多风场相关性的电力系统动态经济调度的方法 |
CN107038292A (zh) * | 2017-04-01 | 2017-08-11 | 三峡大学 | 一种基于自适应多变量非参数核密度估计的多风电场出力相关性建模方法 |
CN107834547A (zh) * | 2017-11-15 | 2018-03-23 | 南通大学 | 一种考虑风电场输出功率关联特性的输电网规划方法 |
CN108074038A (zh) * | 2017-12-11 | 2018-05-25 | 国网江苏省电力有限公司经济技术研究院 | 一种考虑可再生能源及负荷多时空分布特性的电力生产模拟方法 |
Non-Patent Citations (3)
Title |
---|
王小红等: "基于时变Copula函数的风电出力相关性分析", 《电力系统及其自动化学报》 * |
王璐: "《中国股市与债市溢出效应的数量研究》", 31 December 2009, 西南交通大学出版社 * |
赵继超等: "基于Copula理论的风光互补发电系统可靠性评估", 《电力自动化设备》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109753763B (zh) * | 2019-03-06 | 2022-08-23 | 国网江苏省电力有限公司 | 一种风光联合出力概率建模方法 |
CN109753763A (zh) * | 2019-03-06 | 2019-05-14 | 国网江苏省电力有限公司 | 一种风光联合出力概率建模方法 |
CN110287601A (zh) * | 2019-06-27 | 2019-09-27 | 浙江农林大学 | 一种毛竹胸径年龄二元联合分布精确估算方法 |
CN110287601B (zh) * | 2019-06-27 | 2022-11-15 | 浙江农林大学 | 一种毛竹胸径年龄二元联合分布精确估算方法 |
CN110414099A (zh) * | 2019-07-15 | 2019-11-05 | 东南大学 | 一种基于Copula的相关结构模型综合选型方法 |
CN110611422B (zh) * | 2019-09-17 | 2021-01-29 | 福州大学 | 一种考虑mmc模块间相关性的子模块冗余优化配置方法及系统 |
CN110611422A (zh) * | 2019-09-17 | 2019-12-24 | 福州大学 | 一种考虑mmc模块间相关性的子模块冗余优化配置方法及系统 |
CN110705099B (zh) * | 2019-09-30 | 2021-06-11 | 华北电力大学 | 一种检定风电场出力相关性的方法 |
CN110705099A (zh) * | 2019-09-30 | 2020-01-17 | 华北电力大学 | 一种检定风电场出力相关性的方法 |
CN111274542A (zh) * | 2020-01-10 | 2020-06-12 | 国网江苏省电力有限公司扬州供电分公司 | 基于藤copula和混合偏移正态分布的置信容量评估方法及装置 |
CN111162566B (zh) * | 2020-01-21 | 2022-07-01 | 中国电力科学研究院有限公司 | 一种确定风光联合发电系统的联合出力特性的方法和系统 |
CN111162566A (zh) * | 2020-01-21 | 2020-05-15 | 中国电力科学研究院有限公司 | 一种确定风光联合发电系统的联合出力特性的方法和系统 |
CN112163695A (zh) * | 2020-09-09 | 2021-01-01 | 华南理工大学 | 基于copula函数的风力光伏发电预测方法、系统和介质 |
Also Published As
Publication number | Publication date |
---|---|
CN109038648B (zh) | 2020-11-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109038648A (zh) | 一种基于Copula函数的风光联合出力建模方法 | |
Huang et al. | Forecasting hourly solar irradiance using hybrid wavelet transformation and Elman model in smart grid | |
CN109524993B (zh) | 用于中长期优化调度的风电光伏典型周出力场景生成方法 | |
WO2017054537A1 (zh) | 一种长时间尺度光伏出力时间序列建模方法及装置 | |
CN103020423B (zh) | 基于copula函数获取风电场出力相关特性的方法 | |
Huang et al. | Prediction of wind power by chaos and BP artificial neural networks approach based on genetic algorithm | |
Borges et al. | A model to represent correlated time series in reliability evaluation by non-sequential Monte Carlo simulation | |
Cui et al. | Short‐term city electric load forecasting with considering temperature effects: An improved ARIMAX model | |
CN104679738B (zh) | 互联网热词挖掘方法及装置 | |
CN106099932B (zh) | 一种考虑不确定性的时空相关性的日前计划潮流分析方法 | |
CN104573876A (zh) | 基于时序长记忆模型的风电场短期风速预测方法 | |
Khan et al. | Determination of Weibull parameter by four numerical methods and prediction of wind speed in Jiwani (Balochistan) | |
CN103106344A (zh) | 一种建立电力系统聚类负荷模型的方法 | |
CN106055846A (zh) | 一种基于bsnn‑arx的光伏逆变器模型辨识方法 | |
CN111797132B (zh) | 考虑时空相关性的多可再生能源电站功率场景生成方法 | |
CN105162654A (zh) | 一种基于局部社团信息的链路预测方法 | |
Raischel et al. | Uncovering wind turbine properties through two-dimensional stochastic modeling of wind dynamics | |
Morrissey et al. | An improved method for estimating the wind power density distribution function | |
CN111162537B (zh) | 基于组合Copula函数的拉丁超立方抽样法概率潮流计算方法 | |
CN110212592B (zh) | 基于分段线性表达的火电机组负荷调节最大速率估计方法及系统 | |
Sin et al. | Applications of Monte Carlo method in chemical, biochemical and environmental engineering | |
Singh et al. | Parameter estimation of four diode solar photovoltaic cell using hybrid algorithm | |
CN106650060B (zh) | 一种光伏电池内阻衰减系数预测方法 | |
CN112508445A (zh) | 一种基于气温及导线温度量测的架空导线动态增容方法及系统 | |
CN112332416A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20201117 Termination date: 20210710 |
|
CF01 | Termination of patent right due to non-payment of annual fee |