专利分类
专利分类

Mixed-Copula函数的风电场旋转备用容量优化方法专利

专利号:201610395800.7

销售价
1500
Mixed-Copula函数的风电场旋转备用容量优化方法专利二维码
  • 累计销量0
  • 浏览次数31
  • 累计评论0
首页

专利名称:Mixed-Copula函数的风电场旋转备用容量优化方法

技术领域:风电场,容量优化

IPC主分类号:H02J3/38

申请号:CN201610395800.7

公开日:2018-09-25

说明书

基于Mixed-Copula函数的风电场旋转备用容量优化方法

技术领域

[0001] 本发明属于新能源技术以及风电场参与调频容量优化领域,具体涉及一种考虑风电场内部的每台风机实际地貌特征,以及计及相邻列之间的风速相关性的特征,利用Mixed-Copula函数得到风电场出力进而实现优化旋转备用容量的方法。

背景技术

[0002] 基于目前传统备用容量的提供方式,风电场存在很大的弃风量的现象,造成能源的大量损失,现已有考虑将转子惯性动能纳入电网频率调节的调频方式。在此基础上,不仅可以有效降低弃风量,并且可以更大限度发挥风机调频的裕量,因此确定一种合理的旋转备用容量的优化方法是具有重要意义的。
[0003] 现有的旋转备用容量优化方法主要是针对单台风力机进行的旋转备用优化研究,未涉及到风电场内部联合功率分布的内容研究;有文章针对考虑转子惯性动能纳入到风机参与调频容量进行了研究,但是风电场内部特性只考虑了风机间的尾流效应,将风电场内部的同一迎风面的单排风机等价于一台,这种假设并没有考虑到风电场内部节点之间有关风速相关性的影响,同时也没有考虑实际风电场内部地貌特征对风速的影响,这些都会对风电场内部的实际风速分布造成影响,从而对于风电场的实际有功出力的计算产生误差。因此在进行风电场旋转备用容量优化的时候,需要考虑风电场内部风机节点相关性和地貌特征造成的影响。

发明内容

[0004] 本发明目的是在于解决现有风电场的旋转备用容量优化研究中未考虑到风电场内部实际地貌特征以及节点之间风速相关性的问题。
[0005] 为了解决上述技术问题,发明人采用了如下的技术方案:
[0006] 一种基于Mixed-Copula函数的风电场旋转备用容量优化方法,包括以下步骤:
[0007] 选取能够有效反映出尾部特性的Copula函数,通过将反映不同特性的Copula函数进行组合来构成新的Mixed-Copula函数用以描述相邻两列风机风速的联合概率分布;通过非层级聚类的EM算法进行Mixed-Copula函数的参数估计,估计出每个Copula函数的反映尾部特征的权系数和反映风速节点相关性的相关系数,最终得到相邻两列节点风速的联合概率分布;获得相邻两列节点风速的联合概率分布之后,根据风速服从威布尔分布进行逆变换,然后将整个风电场的风速特性进行重新划分,再结合每个节点所处的实际地貌特征,计及相对风电场最低点的海拔高度,将风速进行转化得到整个风电场新的风速分布;然后根据风机有功出力和风速关系的计算公式获得风电场实际有功出力;根据获取最大旋转备用容量为优化的目标函数,在可调节的桨距角和转子转速的区间内确定最优转速和桨距角,以实现风电场的经济效益最大化;最终求解该非线性约束目标优化的问题。
[0008] 具体步骤为:
[0009] 第一步,选取能够有效反映出尾部特性的Copula函数,通过将反映不同特性的Copula函数进行组合来构成新的Mixed-Copula函数用以描述相邻两列风机风速的联合概率分布;
[0010] 通过考虑风速存在尾部特性,选用能够反映出对称尾部特性的Frank-Copula函数和能够反映出非对称尾部特性随机变量的相关关系的Gumbel-Copula和Clayton-Copula函数组成Mixed-Copula函数,其表达式如下:
[0011]
[0012] 其中,Ci分别表示不同的Copula函数,λi分别表示Frank-Copula函数、Gumbel-Copula和Clayton-Copula函数的权重系数,表示在Mixed-Copula函数中占据的比例;θi则表示的不同Copula分布函数的相关系数,表征随机变量之间相关系数;u,v表示两个随机变量服从[0,1]之间的均匀分布。
[0013] 第二步,通过非层级聚类的EM算法进行Mixed-Copula函数的参数估计,估计出每个Copula函数的反映尾部特征的权系数和反映风速节点相关性的相关系数,最终得到相邻两列节点风速的联合概率分布;
[0014] EM是一个在已知部分相关变量的情况下,估计未知变量的迭代技术,其思想是获得的参量期望最大化。通过EM法实现Mixed-Copula函数的参数估计,具体实现步骤如下:
[0015] 首先初始化参数,引入一个未知变量z表征位置:
[0016] z={z1,z2,z3L zi},i取1,2,3…k                (2)
[0017] 其中,zi有且仅有一个为1,其余均为0,用来表征所属的Copula函数。并且在zi=1时可以获得(3)式,p(z)表示随机变量z的概率;
[0018] p(zi=1)=λi                (3)
[0019]
[0020] 对于观测样本x=(y,z),y=(u,v)其条件概率的表达式为:
[0021]
[0022]
[0023]
[0024] 其中, N为采样点数,k为分布个数,X为整个采样样本
[0025] 然后构造一个最大似然函数并求取其期望,分别如式(8)、(9)表示:
[0026]
[0027]
[0028] 其中,
[0029] 得到初始化参数后,确定初值进行EM算法的求解,主要分两步:
[0030] E步骤:第一步是计算期望Q,利用对隐藏变量的现有估计值,计算其最大似然估计值;
[0031] M步骤:求解到最大化的Q,最大化的实现是在E步上求得的最大似然值来计算参数的值。M步上找到的参数估计值被用于下一个E步计算中,这个过程不断交替进行。当满足收敛条件时,停止计算,得到此时条件下的待估参数,从而确定Mixed-Copula函数的待估计的参数。
[0032] 第三步,在获得相邻两列节点风速的联合概率分布之后,根据风速服从威布尔分布进行逆变换,然后将整个风电场的风速特性进行重新划分,再结合每个节点所处的实际地貌特征,计及相对风电场最低点的海拔高度,将风速进行转化得到整个风电场新的风速分布;
[0033] 确定Mixed-Copula函数的参数后,便可以通过(1)式获得表征相邻两列节点风速的联合概率分布,进而通过威布尔逆变换得到新的风速,同理全部风电场进行风速重新划分,得到新的风电场风速分布。
[0034] 又因为实际风电场地貌特征也会对实际风电场风速的分布产生影响,所以每个节点除了要考虑节点的相关性,还应该根据地形对风速进行根据经验公式转化,经验公式如下(10)式:
[0035]
[0036] 其中,V0表示假设整个风电场都位于海平面高度时候的风速值,H表示风机的轮毂高度,h为每个节点实际地形高度距离风电场海平面的海拔高度,α为经验数值。
[0037] 第四步,然后根据风机有功出力和风速的计算公式获得风电场实际有功出力;基于风电场重新划分的风速分布,根据(11)式风机有功出力和风速的三次方成比例的计算公式获得风电场实际有功出力。
[0038]
[0039] 式中,Pwf是风电场基于实时风速计算的风机有功出力,M、N分别表示风电场按照M×N排列分布的行和列,ρ为空气密度,fN为电网额定频率,βij、ωij为相应节点风机的桨距角和转子转速,G为齿轮箱变比,p为风力发电机的极对数,Cp为风力机特性函数,vij为考虑重新规划后的当前节点的风速。
[0040] 第五步,根据获取最大旋转备用容量为优化的目标函数,在可调节的桨距角和转子转速的区间内确定最优转速和桨距角,以实现风电场的经济效益最大化;建立以最大旋转备用容量为目标函数,满足相应的约束条件,以实现风电场经济最大化,同时可以减小弃风量的目的。具体目标函数和约束条件如下:
[0041] min(PMPPT-Pwf)              (12)
[0042]
[0043] 式中,Pmppt(i,j)表示每个节点的风机工作在最大功率点时的有功功率,H为发电机时间常数,ω_optij为当风机运行在最大功率点时候的最优转速,Tres为一次调频时间,Pre为参考功率,ωmax和βmax分别代表最大的转速限制和桨距角限制,Pg是风机实际有功出力,ΔE为转子惯性动能,σ为置信水平。
[0044] 由于(13)式中存在非线性等式约束条件,因此直接求解存在困难,选择用内点法可以解决这一问题,最终确定满足约束条件的最优解,实现风电场能够提供的旋转备用容量最大化,以及增大经济效益。
[0045] 第六步,通过内点法求解该非线性约束目标优化的问题。

附图说明

[0046] 图1是本发明的机遇Mixed-Copula函数的风电场旋转备用容量优化方法的总体实施流程图。
[0047] 图2是本发明中EM算法流程框图。
[0048] 图3是本发明中Mixed-Copula函数联合分布概率密度示意图。
[0049] 图4是本发明中风电场风速重新分布图。

具体实施方式

[0050] 本发明提出了一种基于Mixed-Copula函数的风电场旋转备用容量优化方法。该方法通过考虑风电场内部的单台风机实际地貌特征即风机安装节点的海拔高度,以及计及相邻列之间的风速相关性的特征基础上,将更为接近真实的实际风电功率的输出进行风电场的旋转备用容量优化,以达到最小弃风量的目标,同时优化的结果可以给调频提供更多的有功支持。因为对于M×N排列的风电场内部的每台风机所处节点的风速实际并不完全相同,在关于风电场旋转备用容量优化的研究中,现有文献将每个迎风面的一排风机等效为一个风机的假设,并不能够完全反映风电场的实际风功率分布效果。基于上述原因,本发明首先通过构造一种Mixed-Copula函数,利用其联合概率分布求取相邻两列风机风速的联合概率分布,再通过威布尔逆变换得到两列的计及相关性的风速值,并且根据每个节点的地貌特征进行风速转化,按照此方法可以重新划分整个风电场的风速分布;然后根据风机有功出力与风速关系的计算公式获得风电场实际有功出力,根据获得最大化的旋转备用容量的目标,确定在此目标下的优化指令,即风力机组的转速和桨距角;最后求解该非线性约束的优化目标,实现风电场获得更大的经济效益。
[0051] 结合图1-图4所示,根据本发明的实施例,具体优化方法的实现步骤如下:
[0052] 第一步,选取能够有效反映出尾部特性的Copula函数,通过将反映不同特性的Copula函数进行组合来构成新的Mixed-Copula函数用以描述相邻两列风机风速的联合概率分布;
[0053] 通过考虑风速存在尾部特性,选用能够反映出对称尾部特性的Frank-Copula函数和能够反映出非对称尾部特性随机变量的相关关系的Gumbel-Copula和Clayton-Copula函数组成Mixed-Copula函数,其表达式如下:
[0054]
[0055] 其中,Ci分别表示不同的Copula函数,λi分别表示Frank-Copula函数、Gumbel-Copula和Clayton-Copula函数的权重系数,表示在Mixed-Copula函数中占据的比例;θi则表示的不同Copula分布函数的相关系数,表征随机变量之间相关系数;u,v表示两个随机变量服从[0,1]之间的均匀分布。
[0056] Gumbel-Copula函数:
[0057]
[0058] Frank-copula函数:
[0059]
[0060] Clayton-copula函数:
[0061]
[0062] λi和θi是待估参数。
[0063] 第二步,通过非层级聚类的EM算法进行Mixed-Copula函数的参数估计,估计出每个Copula函数的反映尾部特征的权系数和反映风速节点相关性的相关系数,最终得到相邻两列节点风速的联合概率分布;
[0064] EM是一个在已知部分相关变量的情况下,估计未知变量的迭代技术,其思想是获得的参量期望最大化。通过EM法实现Mixed-Copula函数的参数估计,具体实现步骤如下:
[0065] 首先初始化参数,引入一个未知变量z表征位置:
[0066] z={z1,z2,z3L zi},i取1,2,3…k          (5)
[0067] 其中,zi有且仅有一个为1,其余均为0,用来表征所属的Copula函数。并且在zi=1时可以获得(6)式,p(z)表示随机变量z的概率;
[0068] p(zi=1)=λi                       (6)
[0069]
[0070] 对于观测样本x=(y,z),y=(u,v)其条件概率的表达式为:
[0071]
[0072]
[0073]
[0074] 其中, N为采样点数,k为分布个数,X为整个采样样本
[0075] 然后构造一个最大似然函数并求取其期望,分别如式(11)、(12)表示:
[0076]
[0077]
[0078] 其中,
[0079] 得到初始化参数后,确定初值进行EM算法的求解,主要分两步:
[0080] 第一步(E步):第一步是计算期望Q,利用对隐藏变量的现有估计值,计算其最大似然估计值;
[0081] 第二步(M步):求解到最大化的Q,最大化的实现是在E步上求得的最大似然值来计算参数的值。M步上找到的参数估计值被用于下一个E步计算中,这个过程不断交替进行。当满足收敛条件时,停止计算,得到此时条件下的待估参数,从而确定Mixed-Copula函数的待估计的参数。求解得到的Mixed-Copula函数表征的计及风速和地貌特征的联合风速概率密度分布。
[0082] 第三步,获得相邻两列节点风速的联合概率分布之后,根据风速服从威布尔分布进行逆变换,然后将整个风电场的风速特性进行重新划分,再结合每个节点所处的实际地貌特征,计及相对风电场最低点的海拔高度,将风速进行转化得到整个风电场新的风速分布;
[0083] 确定Mixed-Copula函数的参数后,便可以通过(1)式获得表征相邻两列节点风速的联合概率分布,进而通过威布尔逆变换得到相邻两列的计及相关性的风速,同理全部风电场进行风速重新划分,得到新的风电场风速分布。
[0084] 又因为实际风电场地貌特征也会对实际风电场风速的分布产生影响,所以每个节点除了要考虑节点的相关性,还应该根据地形对风速进行根据经验公式转化,经验公式如下(13)式:
[0085]
[0086] 其中,V0表示假设整个风电场都位于海平面高度时候的风速值,H表示风机的轮毂高度,h为每个节点实际地形高度距离风电场海平面的海拔高度,α表示经验参数。
[0087] 第四步,然后根据风机有功出力和风速的三次方成比例的计算公式获得风电场实际有功出力;基于风电场重新划分的风速分布。
[0088] 根据(14)式风机有功出力和风速的计算公式获得风电场实际有功出力。
[0089]式中,Pwf是风电场基于实时风速计算的风机有功出力,M、N分别表示风电场按照M×N排列分布的行和列,ρ为空气密度,fN为电网额定频率,βij、ωij为相应节点风机的桨距角和转子转速,G为齿轮箱变比,p为风力发电机的极对数,Cp为风力机特性函数,vij为考虑重新规划后的当前节点的风速。
[0090] 第五步,根据获取最大旋转备用容量为优化的目标函数,在可调节的桨距角和转子转速的区间内确定最优转速和桨距角,以实现风电场的经济效益最大化;建立以最大旋转备用容量为目标函数,满足相应的约束条件,以实现风电场经济最大化,同时可以减小弃风量的目的。具体目标函数和约束条件如下:
[0091] min(PMPPT-Pwf)                 (15)
[0092]
[0093] 式中,Pmppt(i,j)表示每个节点的风机工作在最大功率点时的有功功率,H为发电机时间常数,ω_optij为当风机运行在最大功率点时候的最优转速,Tres为一次调频时间,Pre为参考功率,ωmax和βmax分别代表最大的转速限制和桨距角限制,Pg是风机实际有功出力,ΔE为转子惯性动能,σ为置信水平。
[0094] 目标函数表示考虑调频裕量时风力机应该减少的出力值,因为考虑风电场的内部风速相关性以及实际节点的海报高度后,其每个节点的功率分布可以按照独立性直接累加,从而求得整个电场的功率分布。
[0095] 由于(16)式中存在非线性等式约束条件,因此直接求解存在困难,选择用内点法可以解决这一问题,最终确定满足约束条件的最优解,实现风电场能够提供的旋转备用容量最大化,以及增大经济效益。
[0096] 第六步,通过内点法求解该非线性约束目标优化的问题。
[0097] 选取一个M×N排列的风电场,其中M值为5,N值为6,按照上面的步骤先选取Copula函数,然后通过EM求解函数的权重系数和相关系数。EM计算的初始权重参数为[1/3,1/3,1/3],初始的相关参数为[6,1,4];经过EM求解计算后得到新的权重系数[0.30,0.31,0.39]和相关参数[8.84,4.23,8.72],得到的Mixed-Copula函数表征的计及风速和地貌特征的联合风速概率密度分布如图(3)所示,能够反映出风速的尾部特性,即能够对风速出现极大值和极小值得情况进行体现,图中两个尖峰值分别反映的是上下尾部特性的影响。
[0098] 重新划分的风速分布如图(4)所示。两个风速进行比较可以看出重新规划的风速也能够拟合威布尔分布,同时又包含了风速相关性和风电场实际地形中海拔高度的不同影响。
[0099] 然后,再根据步骤四、五、六进行旋转备用容量的优化,最终求解出优化方案,在可调节的桨距角和转子转速的区间内确定最优转速和桨距角,以实现风电场的经济效益最大化。
[0100] 本发明所属技术领域中具有通常知识者,在不脱离本发明的精神和范围内,当可作各种的改动与润饰。因此,本发明的保护范围当视权利要求书所界定者为准。

Mixed-Copula函数的风电场旋转备用容量优化方法委托购买说明

填写需求表单支付预付款

平台根据需求优化购买方案

确认购买方案支付尾款

平台办理变更等待成功通知

Mixed-Copula函数的风电场旋转备用容量优化方法购买流程说明

发起委托,需要先支付100元预付款,委托不成功,全额退返预付款;

平台收到需求后,会在第一时间联系您,给到您最佳购买方案;

您在确认购买方案后,需支付全额专利购买费,预付款可抵扣购买费,专利购买费具体参见下方表格;

平台确认收款后,将帮您办理专利购买、专利过户等全流程手续;

平台代购专利失败,将全额退返专利购买费,包括预付款;

Mixed-Copula函数的风电场旋转备用容量优化方法专利购买费用

授权未缴费=专利裸价+著录项变更(200元)+登办费(当年年费+5元印花税)+恢复权利请求费1000元(按实收)+委托服务费(200元)+税金(专利裸价+委托服务费)x6%

已下证=专利裸价+著录项变更(200元)+滞纳金(按实收)+恢复权利请求费1000元(按实收)+委托服务费(200元)+税金(专利裸价+委托服务费)x6%

Mixed-Copula函数的风电场旋转备用容量优化方法购买费用说明

专利转让费用

专利买卖交易资料

Q:办理专利转让的流程及所需资料

A:专利权人变更需要办理著录项目变更手续,有代理机构的,变更手续应当由代理机构办理。

1:专利变更应当使用专利局统一制作的“著录项目变更申报书”提出。

2:按规定缴纳著录项目变更手续费。

3:同时提交相关证明文件原件。

4:专利权转移的,变更后的专利权人委托新专利代理机构的,应当提交变更后的全体专利申请人签字或者盖章的委托书。更多

Q:专利著录项目变更费用如何缴交

A:(1)直接到国家知识产权局受理大厅收费窗口缴纳,(2)通过代办处缴纳,(3)通过邮局或者银行汇款,更多缴纳方式

Q:专利转让变更,最快多久能出结果

A:著录项目变更请求书递交后,一般1-2个月左右就会收到通知,国家知识产权局会下达《转让手续合格通知书》。

更多专利转让常见问题

动态评分

0.0

没有评分数据
没有评论数据
 
X 顶部大图