一种基于概率统计的平均风向计算方法与流程

本发明涉及风力发电领域，特别涉及一种基于概率统计的平均风向计算方法。
背景技术：
：一、名词解释⑴风向：风吹来方向。测风塔实测风向数据记录间隔一般为10分钟，即间隔10分钟保存一次数据，单位一般采用(°)。⑵方位：方向位置。将360°圆周以22.5°为间隔划分为16个风向扇区，每个扇区对应一个方位，分别为：N(正北)、NNE(东北偏北)、NE(东北)、ENE(东北偏东)、E(正东)、ESE(东南偏东)、SE(东南)、SSE(东南偏南)、S(正南)、SSW(西南偏南)、SW(西南)、WSW(西南偏西)、W(正西)、WNW(西北偏西)、NW(西北)、(NNW西北偏北)。二、背景说明风电是目前最具有发展前景的可再生清洁能源之一，相比其他形式的可再生能源，具有技术成熟、成本相对较低、对环境影响较小等优势。随着并网与消纳等“下游因素”影响日益凸显，我国风电“版图”正在发生变化，即发展热点由风能资源丰富区逐步转向并网条件相对较好的低风速地区，主要集中在我国中南和西南地区。我国低风速地区的显著特点为：资源条件较差、地形条件复杂、气候条件恶劣等。因此，如何保证低风速风电场风能资源评估工作的精细化、科学化、准确化对风电场规划设计、风功率预测、后评估等工作起着至关重要的作用。风电场场址范围内代表性测风塔的实测数据可真实客观地反映该区域的风能资源情况，是进行风能资源评估的基础。《风电场风能资源评估方法》(GB/T18710-2002)明确规定：进行风能资源评估时，应取得至少连续一年完整的风场逐小时数据，并且实测数据应具有长期代表性。但实际工程中，测风塔数据质量难以保证，主要存在以下几方面问题：①高海拔山地风电场，受冰冻及其他因素影响，测风设备会出现故障或停测现象，导致测风数据存在大量缺测或不合理数据，必须进行缺测、不合理数据的插补替换。②代表性测风塔测风年一般不是平风年(即为大风年或小风年)，实测数据不能反映该地区长期的平均水平，必须进行代表年订正。③风电场开发业主实际需要，项目前期工作开展阶段测风塔实际收集数据未满一年，需要将实测数据延长至满一年。针对以上问题，通常的解决办法是根据临近测风塔、当地参证气象站或中尺度数据进行测风塔数据的插补、替换、延长以及代表年订正。当目标测风塔周围没有其他参照测风塔时，以上工作只能依据参证气象站或中尺度数据进行，而参证气象站或中尺度数据记录值均为小时平均值，因此需要首先统一两者的时间尺度，然后建立相关关系。小时平均风速一般直接采用算数平均方法直接计算得到，而小时平均风向的计算方法则有多种。目前，最常用的方法有：算术平均法、首风向代表法以及利用专业测风数据分析软件。其中首向代表法是直接取该小时内的第一个10分钟风向代表本小时段的平均风向；算术平均法是取该小时内的6个10分钟平均风向值得算术平均值作为本小时段的平均风向。实际工程设计中，设计单位通常采用北京木联能软件技术有限公司开发的测风数据验证和评估软件直接进行小时平均风向计算，该软件小时平均风向的计算原理为矢量平均法，即首先将风向按单位长度投影到X轴和Y轴，同时将加上对应风速权重，再对这些风向在X轴和Y轴上的投影分别求平均值，最后在进行矢量求和。2012年，韩爽等人提出经过多次数据筛选，剔除偏离主风向较远的个别数据，形成优化矢量平均法。目前，风能资源评估中涉及到小时平均风向的处理手段主要由以下两种：⑴设计人员手动处理设计人员手动处理主要是根据实际工程经验，采用首风向代表法和算术平均法两种计算小时平均风向。两种方法小时平均风向计算结果受人为因素的影响均较大，且计算结果的质量难以保证，难以满足可研、初设等设计深度较高阶段的要求。主要弊端体现如下：①当10分钟风向变化范围较大或者风向变化频繁时，首向代表法计算的小时平均风向可能在一定程度上偏离实际主风向。②当相邻两个10分钟风向之差大于180°时，算术平均方法计算的小时平均风向明显错误，与实际风向反向。⑵商业软件处理目前，国内从事风电场工程设计的设计单位(企业)普遍采用由北京木联能软件技术有限公司开发的测风数据验证和评估软件进行风能资源评估，该软件计算小时平均风向的原理为矢量平均法。另外，该软件属商业软件，须付费取得电子狗(按个计)后方可使用，且需定期更新。采用测风数据验证和评估软件计算小时平均风向可以满足风电场全设计阶段的要求，但受软件自身商业诉求的制约，主要存在以下弊端：①通常情况下，受电子狗个数的限制，实际工程应用过程中只能单人单任务进行操作，严重影响工作效率。另外，若设计人员从事国外项目，受限于电子狗，难以实现现场办公。②严格意义上，矢量平均法只是从风向物理意义角度(即风向是一个矢量)对其进行小时平均，并未考虑气温、气压、地形、设备等其他影响风向测量值的外在因素，这会导致计算得到的小时平均风向与实际风向存在一定的偏差。综上所述，无论是采用手动处理还是软件处理，小时平均风向计算均会产生大量不合理数据，并且两种处理手段都是将小时平均风向计算与不合理数据检验两个过程分开进行，即顺序过程(先计算，后验证)，这样耗时费力，工作量大。技术实现要素：本发明要解决的技术问题是：提供一种基于概率统计的平均风向计算方法，使工作人员快速、准确的确定某段时间的平均风向。为解决上述问题，本发明采用的技术方案是：将一段时间内实测的风向按16个风向扇区进行概率统计，以发生概率最大的风向扇区对应的中心风向值作为该段时间的平均风向，并对不合理的平均风向值进行修正。进一步的，本发明的具体步骤如下：步骤一、导入一段时间内的风向数据；步骤二、将风向数据按16个风向扇区划分标准由风向角度值转化为扇区方位；步骤三、统计该段时间内各扇区方位的频率，以发生概率最大的风向扇区对应的中心风向值作为该段时间的平均风向；步骤四、检验所得的平均风向是否合理，若不合理，则对平均风向进行修正。进一步的，16个风向扇区划分标准为：将360°圆周以22.5°为间隔划分为16个风向扇区。进一步的，当出现多个风向频率相等的情况，则筛选并剔除该时段内明显偏离主风向的方位，重新计算风向频率，以发生概率最大的风向扇区对应的中心风向值作为该段时间的平均风向；若两次剔除后，仍出现多个风向频率相等的情况，则以该段时间内最大风速对应的方位值作为该段时间的平均风向。进一步的,检验平均风向是否合理的方法为：令待检平均风向为p，令待检平均风向后面时间的第一个平均风向为p2，令待检平均风向后面时间的第二个平均风向为p2；检验p与p1的差值是否小于等于45°，若是，则说明p合理；否则继续检验p1与p2的差值是否小于等于45°，若是，则说明p不合理。进一步的,对平均风向值进行修正的方法包括：a.令待修正平均风向所在时间段为t，将t内最大风速对应的风向值作为t的平均风向；b.按所述检验平均风向是否合理的方法，对步骤a中确定的t的平均风向进行检验，判断是否合理，若是，则将所述最大风速对应的风向值作为修正后的值，修正结束；否则进入步骤c；c.按所述检验平均风向是否合理的方法，检验当t内除最大风速对应的风向值以外的其他风向值作为t内的平均风向时，是否有合理的平均风向，若有，则将其中一个合理的平均风向作为修正后的值，修正结束；若没有，则取t的前一时间段或后一时间段的平均风向值作为t的平均风向，修正结束。本发明的有益效果是：本发明不但可以减少设计人员对商业软件的过度依赖，为自用软件开发奠定良好基础；还将外在因素对风向的影响考虑在内，真实反映风向客观变化规律；同时将不合理数据验证与小时平均风向计算同步进行，即并行过程，提高了小时平均风向的数据质量，有效降低了工作量，利于提高工作效率。附图说明图1是本发明的流程图。图2是实测风向值风向扇区划分示意图。具体实施方式如图1所示，本发明步骤如下：步骤一、导入一段时间内的风向数据；步骤二、将风向数据按16个风向扇区划分标准由风向角度值转化为扇区方位；步骤三、统计该段时间内各扇区方位的频率，以发生概率最大的风向扇区对应的中心风向值作为该段时间的平均风向；步骤四、检验所得的平均风向是否合理，若不合理，则对平均风向进行修正。步骤五、导出平均风向。以下以计算小时平均风向为例，对本发明做进一步说明。1、基础数据导入基础数据导入主要包括：①风塔实逐10分钟时间序列；②对应的逐10分钟实测风向值，数据格式如表1所示。表1基础数据格式示意表2、风向标准化风向标准化主要包括：①实测风向值按16个风向扇区进行划分；②实测风向值转化为相对应的扇区方位。第一步，实测风向值按16个风向扇区进行划分，划分标准如表2所示：将360°圆周以22.5°为间隔划分为16个风向扇区，划分后的直观图如图2所示；表2实测风向值扇区划分标准表2中，以N为第1扇区，顺时针依次NNE为第2扇区，NE为第3扇区，……，NNW为第16扇区。第二步，将风向数据按16个风向扇区划分标准由风向角度值转化为扇区方位。将表1中风向角度值转化为扇区方位后如表3所示。表3风向标准化格式示意表时间(Date&TimeStamp)风向(Vg)扇区2016-11-19:00SSW102016-11-19:10SSW102016-11-19:20SW112016-11-19:30SW112016-11-19:40SW112016-11-19:50WSW122016-11-110:00SW112016-11-110:10SW112016-11-110:20WSW122016-11-110:30SW112016-11-110:40SW112016-11-110:50SW11………3、风向频率统计风向频率统计主要包括：①计算逐小时内各实测风向频率；②取实测风向频率最大对应的风向扇区方位作为小时平均风向扇区方位；③小时平均风向扇区方位转化为相对应的扇区中心方向角度值。第一步，计算逐小时内各实测风向频率，计算公式如下，并以风向频率最大对应的风向扇区方位作为小时平均风向扇区方位，数据格式详见表4。式中：逐小时内风向方位的总观测次数应为6；逐小时内各风向方位出现的次数应不大于6。表4小时平均风向扇区方位格式示意表例：2016-11-19:00至2016-11-19:50时段内：P(SSW)＝2/6*100％＝33.3％；P(SW)＝3/6*100％＝50.0％(最大)；P(WSW)＝1/6*100％＝16.7％。当出现多个风向频率相等的情况，说明该时段内风向变化频繁，此时需要筛选并剔除该时段内明显偏离主风向的方位，重新计算风向频率，直到某风向频率出现最大值。若两次剔除后，仍出现多个风向频率相等的情况，说明采用概率统计法计算本小时的平均风向意义不大，可用该小时内最大风速对应的方位作为本小时的平均风向方位。第二步，将小时平均风向扇区方位转化为对应扇区中心方向角度值，转化的标准表如表5所示。小时平均风向方位统一为对应扇区中心方向角度值的数据格式见表6。表5小时平均风向扇区方位转化为对应扇区中心方向角度值参照表表6小时平均风向方位统一为对应扇区中心方向角度值的数据格式示意表时间(Date&TimeStamp)小时平均风向方位小时平均风向值(Vsg)2016-11-19:00SSW202.52016-11-110:00SW225………4、数据检验及修正将得出的小时平均风向值进行数据验证并进行不合理数据的处理，并且该过程是与步骤3中第二步并行实施，参照《风电场风能资源评估方法》(GB/T18710-2002)规范中对不合理数据的检验标准。具体如下：第一步，检验平均风向是否合理检验相邻两个小时平均风向差值，一般是待测小时(以9:00为例)和待测小时后一小时(10:00)，判断向差值是否小于等于45°，若满足，则说明9:00的平均风向数据合理；若不满足，继续检验后延一小时相邻两个小时(即10:00和11:00)平均风向差值是否小于等于45°，若满足，则说明9:00的平均风向数据不合理，反之说明10:00平均数据不合理。第二步，对平均风向值进行修正对不合理小时平均风向数据进行处理，考虑到风向的随机性以及瞬间不可突变性，返回该小时内原始6个10分钟实测风向值，按照步骤4的第一步检验该小时内最大风速对应的风向是否满足要求，满足则取最大风速对应的风向值修正为本小时的平均风向，若不满足，重复步骤4的第一步检验其他5个风向值，直到满足要求。若6个风向值均不满足要求，则直接采用前一小时或后一小时平均风向值作为该时刻的小时平均风向。5、目标数据导出将检验及修正后的数据导出，目标数据导出后，增加小时平均风速数据，便可导入风电场风能资源评估专业软件(WT、WAsP等)进行风电场区域风能资源分布模拟。实际验证效果显示，本发明准确性高、应用性强。以上描述了本发明的基本原理和主要的特征，说明书的描述只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。当前第1页1 2 3