基于混合Copula理论的风电场多机聚合等值方法技术

本发明专利技术公开了一种基于混合Copula理论的风电场多机聚合等值方法。通过构建能够适应复杂地形环境下的风电场风速预测模型，考虑了风机的排列布置产生的多尾流效应，同时还考虑了地面粗糙度和障碍物遮挡的影响，形成精确的风电场风速分布模型。然后运用Copula理论，提出计及尾部特性的Copula函数来表征风速相关性，构建Mixed‑Copula函数来描述两台风机风速的联合分布，继而通过检验Spearman秩相关系数，判断两变量的相关程度。相关程度高的两风机进行等值处理划分为一类，实现了风电场等值简化。

【技术实现步骤摘要】
基于混合Copula理论的风电场多机聚合等值方法
本专利技术涉及一种能够适应复杂地形环境下的风电场多机聚合等值方法。
技术介绍
由于风电场内部地形空间因素的差异，造成相同风机类型风机风速也存在差异，已有的文献对风电场内部的风机进行了多机等值或单机等值的研究，然而地形差异将风机进行单机等值的模型并不精确。为了更为接近真实的实际风速进行风电场多机群聚合等值，更好的描述相邻风机之间的变量关系并实现多机群聚合等值。传统法的等值是同调等值法，其分类思想是根据发电机的功角进行划分。但是风机与传统的发电机不同，现有根据风速进行划分，相同风速环境的情况下，风况简单，等值方便。已有研究使用K-MEANS算法，将风机的状态变量矩阵作为划分聚类的分类标准。但是传统的K-MEANS算法有着跟随状态矩阵的变化而修改划分类别的缺点，由于风速的波动性带来的是反复进行分类和划分的变更，大大增加了计算难度和运算时间。因此，在对风电场多机聚合时，选择风速切入风速作为等值参考条件，根据切入风速的标准选择等值方法。
技术实现思路
本专利技术目的是在于解决现有风电场多机聚合等值研究中未考虑到风电场内部实际地貌特征以及节点之间风速相关性的问题。为了解决上述技术问题，专利技术人采用了如下的技术方案：一种基于混合Copula理论的风电场多机聚合等值方法，包括以下步骤：通过考虑复杂地形环境下影响风速的因素，主要包括风机排列布置带来的尾流效应、地面粗糙度、障碍物遮挡，运用WAsP软件构建复杂地形的风电场模型，通过仿真获得风速数据。利用历史数据统计得出风机形状参数和尺度参数。选取能够有效反映出尾部特性的Copula函数，通过将反映不同特性的Copula函数进行组合来构成新的Mixed-Copula函数用以描述相邻两列风机风速的联合概率分布；通过非层级聚类的EM算法进行Mixed-Copula函数的参数估计，估计出每个Copula函数的反映尾部特征的权系数和反映风速节点相关性的相关系数，最终得到相邻两列节点风速的联合概率分布；获得相邻两列节点风速的联合概率分布之后，根据风速服从威布尔分布进行逆变换，然后将整个风电场的风速特性进行重新划分，得到整个风电场新的风速分布；基于构造Mixed-Copula函数对两个变量的相关性进行分析，分别计算秩相关系数矩阵，根据矩阵计算得出两变量的相关性，通过判定相关性程度，当判定视为有较强的相关性连接时，可以等值为一种类型，从而对风电场进行重新等值划分。具体步骤为：第一步，通过考虑复杂地形环境下影响风速的因素，主要包括风机排列布置带来的尾流效应、地面粗糙度、障碍物遮挡，运用WAsP软件构建复杂地形的风电场模型，通过仿真获得风速数据，利用历史数据统计得出风机形状参数和尺度参数。WAsP软件输入数据包括气象数据、地面粗糙度、障碍物遮挡数据，下面主要介绍数据来源和组成。气象数据由当地的气象站、台提供的时间序列数据。其内容主要包含风速，风向，测量点标准气压、温度以及海拔。WAsP内的风向数据被划分为12个扇区，均匀的将360°划分为每30°一个扇区，根据国际惯例，将风速划分至相应扇区。地面粗糙度数据据地形的不同情况，粗糙度可以划分为若干个等级。在一定距离内，地面情况越复杂，粗糙度等级越高，变化层次越多，粗糙度越大对风的影响越大。障碍物遮挡数据，一般将障碍物视为长度、宽度、高度分别取固定值的长方体来考虑。考虑障碍物到某点的距离和方位，障碍物的输入可直接输入障碍条件，也可以手动输入。通过计算，WAsP可以输出的数据包括：每台风机的平均风速和极限风速；风向玫瑰图、风速段风向图；威布尔分布拟合参数等，获得风速数据。第二步，选取能够有效反映出尾部特性的Copula函数，通过将反映不同特性的Copula函数进行组合来构成新的Mixed-Copula函数用以描述相邻两列风机风速的联合概率分布；通过考虑风速存在尾部特性，选用能够反映出对称尾部特性的Frank-Copula函数和能够反映出非对称尾部特性随机变量的相关关系的Gumbel-Copula和Clayton-Copula函数组成Mixed-Copula函数，其表达式如下：其中，Ci分别表示不同的Copula函数，λi分别表示Frank-Copula函数、Gumbel-Copula和Clayton-Copula函数的权重系数，表示在Mixed-Copula函数中占据的比例；θi则表示的不同Copula分布函数的相关系数，表征随机变量之间相关系数；u，v表示两个随机变量服从[0,1]之间的均匀分布；K最大取值范围为3；i取1，2，3。第三步，通过非层级聚类的EM算法进行Mixed-Copula函数的参数估计，估计出每个Copula函数的反映尾部特征的权系数和反映风速节点相关性的相关系数，最终得到相邻两列节点风速的联合概率分布；EM是一个在已知部分相关变量的情况下，估计未知变量的迭代技术，其思想是获得的参量期望最大化。通过EM法实现Mixed-Copula函数的参数估计，具体实现步骤如下：首先初始化参数，引入一个未知变量z表征位置：z＝{z1,z2,z3…zi}，i取1,2,3…k(2)其中，zi有且仅有一个为1，其余均为0，用来表征所属的Copula函数。并且在zi＝1时可以获得(3)式，p(z)表示随机变量z的概率；p(zi＝1)＝λi(3)对于观测样本x＝(y，z)，y＝(u，v)其条件概率的表达式为：其中，N为采样点数，k为分布个数，X为整个采样样本然后构造一个最大似然函数并求取其期望，分别如式(8)、(9)表示：其中，表示i采样点下的函数，表示i+1采样点的函数；yi表示采样点i的y函数；得到初始化参数后，确定初值进行EM算法的求解，主要分两步：第1步(E步)：第一步是计算期望Q，利用对隐藏变量的现有估计值，计算其最大似然估计值；第2步(M步)：求解到最大化的Q，最大化的实现是在E步上求得的最大似然值来计算参数的值。M步上找到的参数估计值被用于下一个E步计算中，这个过程不断交替进行。当满足收敛条件时，停止计算，得到此时条件下的待估参数，从而确定Mixed-Copula函数的待估计的参数。第四步获得相邻两列节点风速的联合概率分布之后，根据风速服从威布尔分布进行逆变换，然后将整个风电场的风速特性进行重新划分，将风速进行转化从而得到整个风电场新的风速分布。第五步基于构造Mixed-Copula函数对两个变量的相关性进行分析，分别计算秩相关系数矩阵，根据矩阵计算得出两变量的相关性，通过判定相关性程度，当判定视为有较强的相关性连接时，可以等值为一种类型，从而对风电场进行重新等值划分。将两个变量U和V各自排序形成R和S秩相关系数矩阵，秩的大小表征了其变量值代表的大小的次序，当U和V完全相关时，每一项D＝(Ri-Si)＝0，通过这个来度量U和V的相关程度，当D越大，表征U和V偏离程度越大，相关程度越低。因为D值可正可负，所以选择用D的平方值来观察。Spearman秩相关等级系数是衡量两个变量相关程度的重要指标，基于上述所说的基础，给出Spearman秩相关假设检验步骤：1)假设H0：U和V不相关；H1：U和V相关。2)计算假设检验的统计量rs：式中，3)判断是否给出假设检验的置信水平，当没有给出假设检验的置本文档来自技高网...

【技术保护点】
一种基于混合Copula理论的风电场多机聚合等值方法，其特征在于，包括以下的步骤：1)通过考虑复杂地形环境下影响风速的因素，运用WAsP软件构建复杂地形的风电场模型，通过仿真获得风速数据；利用历史数据统计得出风机形状参数和尺度参数；2)选取能够有效反映出尾部特性的Copula函数，通过将反映不同特性的Copula函数进行组合来构成新的Mixed‑Copula函数用以描述相邻两列风机风速的联合概率分布；3)通过非层级聚类的EM算法进行Mixed‑Copula函数的参数估计，估计出每个Copula函数的反映尾部特征的权系数和反映风速节点相关性的相关系数，最终得到相邻两列节点风速的联合概率分布；4)获得相邻两列节点风速的联合概率分布之后，根据风速服从威布尔分布进行逆变换，然后将整个风电场的风速特性进行重新划分，得到整个风电场新的风速分布；5)基于构造的Mixed‑Copula函数对两个变量的相关性进行分析，分别计算秩相关系数矩阵，根据矩阵计算得出两变量的相关性，通过判定相关性程度，当判定视为有较强的相关性连接时，等值为一种类型，从而对风电场进行重新等值划分。

【技术特征摘要】
1.一种基于混合Copula理论的风电场多机聚合等值方法，其特征在于，包括以下的步骤：1)通过考虑复杂地形环境下影响风速的因素，运用WAsP软件构建复杂地形的风电场模型，通过仿真获得风速数据；利用历史数据统计得出风机形状参数和尺度参数；2)选取能够有效反映出尾部特性的Copula函数，通过将反映不同特性的Copula函数进行组合来构成新的Mixed-Copula函数用以描述相邻两列风机风速的联合概率分布；3)通过非层级聚类的EM算法进行Mixed-Copula函数的参数估计，估计出每个Copula函数的反映尾部特征的权系数和反映风速节点相关性的相关系数，最终得到相邻两列节点风速的联合概率分布；4)获得相邻两列节点风速的联合概率分布之后，根据风速服从威布尔分布进行逆变换，然后将整个风电场的风速特性进行重新划分，得到整个风电场新的风速分布；5)基于构造的Mixed-Copula函数对两个变量的相关性进行分析，分别计算秩相关系数矩阵，根据矩阵计算得出两变量的相关性，通过判定相关性程度，当判定视为有较强的相关性连接时，等值为一种类型，从而对风电场进行重新等值划分。2.根据权利要求1所述的基于混合Copula理论的风电场多机聚合等值方法，其特征在于，所述步骤1)中，构建能够适应复杂地形环境下的风电场风速预测模型，考虑风机的排列布置产生的多尾流效应，同时还考虑了地面粗糙度和障碍物遮挡的影响，运用WAsP软件构建复杂地形的风电场模型，通过仿真获得风速数据，利用历史数据统计得出风机形状参数和尺度参数。3.根据权利要求1所述的基于混合Copula理论的风电场多机聚合等值方法，其特征在于，所述步骤2)和3)中，通过考虑风速存在尾部特性，选用能够反映出对称尾部特性的Frank-Copula函数、对下尾部变化比较敏感Gumbel-Copula函数以及对上尾部相关变化敏感Clayton-Copula函数组成Mixed-Copula函数，其表达式如下：其中，Ci分别表示不同的Copula函数，λi分别表示Frank-Copula函数、Gumbel-Copula和Clayton-Copula函数的权重系数，表示...

【专利技术属性】
技术研发人员：周前，张梦月，卜京，王俊，杨阳，汪成根，张宁宇，胡昊明，
申请(专利权)人：国网江苏省电力公司电力科学研究院，国家电网公司，江苏省电力试验研究院有限公司，
类型：发明
国别省市：江苏,32