下载:
-
风能资源是气候资源的重要组成部分，相对于不可再生的化石燃料，风能是最具潜力的替代能源，也是永久性的、清洁的可再生能源[1-4].进行大型风能资源的开发时，前期风资源评估和风电场选址十分关键，通常复杂地形应根据实际情况布设测风塔进行观测，评估后再进行开发[5].在面积辽阔的四川省，盆地内气候特点为少风多云，通常被认为风能资源不佳[4, 6-8]；在部分盆地周围的高山和高原地区，气候和地形与盆地截然不同，仍然存在部分风能资源较好地带[8-10]，且地广人稀的特点利于修建风电场.以往四川省风能资源评估的研究，大多基于四川省基本气象站10 m高度观测数据，采用插值方法来研究空间分布情况.传统方法不能满足现有风能资源开发的要求：气象站本身数量分布有限，随着城市发展，现多位于城区附近，地理位置不满足开展风电机安装的要求[11]；传统气象站测风高度不能满足现有风电发电机要求；风能资源具有很强的局地性，传统插值方法无法满足现有开发要求[12-15].对于风能资源实地开发，还需实地竖立测风塔，进行为期一年以上的风能数据观测，才能更好地反应局地风能资源是否满足开发的情况.目前，关于四川省“三州”地区风能资源的研究较少，关于布拖县风能资源目前尚无探寻[16-19].本研究基于位于四川省凉山州南部的布拖县测风塔的实测数据和气象观测站长期数据，对布拖县局地风资源情况予以分析和评估，以此揭示布拖县该区域内的风资源特点，为该区域的风能资源评估提供方法和资料，为风电场选址、规划和建站等风能资源开发工作提供相应参考材料.
全文HTML
-
本研究所用的测风塔风速资料来源于四川省发展改革委员会风能资源数据库，为四川省凉山州彝族自治区布拖县两座测风塔(A，B)2012年为期一年的逐时实测风速资料；订正数据来自四川省气象局，为布拖气象观测站1981-2010年共计30年的风速逐时数据；本研究所用分析数据为A，B测风塔经距离最近的布拖县气象局风速资料，采用长年代法[5]订正推算后的代表年数据.测风塔A，B高70 m，地处海拔3 300 m左右的较为平坦的高山山脊，与布拖气象站距离在25 km以内；设备的安装无塔影效应，测风仪器为美国NRG系统，测风仪均经过标定，风向标安装位置根据当地磁偏角修正.
-
空气密度直接影响风能的大小，在同等风速条件下，空气密度越大风能越大.空气密度计算公式如下：
式中：ρ为气密度(kg/m3)；P为气压(hPa)；t为气温(℃)；e为水汽压(hPa).
根据测风塔的实测气温、气压、相对湿度观测数据，计算测风塔空气密度值.
-
式中：VE为平均风速；n为计算时段内(年、月)平均风速序列的个数；Vt为观测的风速序列[9].
空气密度计算公式为：
式中：ρ为空气密度(kg/m3)；R为气体常数；T为年平均开氏温标绝对温度；P为年平均气压(Pa)[14].
用(2)式计算所得实际观测高度处的空气密度结果，依据以下公式推算轮毂高度处的空气密度：
式中：ρz为订正轮毂高度处z的空气密度；ρh为实地观测高度处的空气密度[20].
年平均风功率密度为：
式中：νi为风速观测序列(m/s)；D为年平均风功率密度(W/m2)；ρ为空气密度(kg/m3)；vi为逐小时风速(m/s)；n为统计时段小时数(h)[11].
-
由于目前测风塔观测时间一般比较短，难以代表当地长年平均风况特征，为满足运行期长达20年的风电场风能资源评估需要，利用布拖气象观测站30年风速资料，对测风塔一年的风速实测资料进行订正，对拟建风电场区域的风能资源进行长年代评估.计算各参证站观测年度年平均风速相对本站近20年累年平均风速距平百分率，计算公式为：
式中：η为累年平均风速距平百分率；V为参证站在现场观测时段的平均风速；V为参证站累年平均风速.根据算出的参证站风速距平百分率η和测风塔观测年度的风速值V，得出测风塔长年代风速订正公式：
长年代风功率密度的计算公式为：
式中：D为测风塔观测时段的风功率密度.
-
某处风向分布特征通常用风向在16方位出现的频率来描述，这里的风向频率指在一定时段内，风出现在各个方位的次数占风向在全方位出现的总次数比值.
风能密度计算公式为：
式中：DWE为一定时段内风能密度(W·h/m2)；vi为第i个记录的风速(m/s)值；n为一定时段内记录总数；ρ为空气密度，由空气密度公式(2)和(3)给出.
风能的密度分布情况通常由一定时间段内，风能密度在各方向上与全方位的总风能密度比值来描述.
-
风向频率通过统计观测时间段内所观测的风向资料计算. 16个方位各风向在统计时间段内(年或月)小时数量与总观测小时数的比值就是该点风向频率.
风能方向频率的计算：根据一定时间段内风速和风向的逐时观测资料，计算16个方位统计需计算各方位具有的风能能量和总能量的比值，该值即为此方向的风能频率[9, 13].年风能方向频率计算公式如下：
式中：F东为一年内东方位的风所具有能量与该处总能量的比值；i为东风风向出现风的小时数量.
风速频率指的是每1 m/s即为一个风速区间，在要统计计算年份的测风序列中，各风速区间内风速出现的频率.每个风速区间取的中间值，比如风速区间6.6 m/s到7.5 m/s[13]则取中间值7 m/s.
-
通常用平均风速和标准差来估算K和A威布尔参数.
平均风速ν用μ来估算，σ用准差Sv来估算：
式中：Vi为风速观测数据序列；n为特定时段内的风速序列个数[20].
Weibull两参数K，A通过以下两式计算：
式中：
${\mathit{\Gamma }\left( {1 + \frac{1}{k}} \right)}$ 为伽马函数.采用威布尔分布曲线来与风速频率曲线拟合，其中概率密度函数由下式来计算：
式中：f(x)为概率密度函数；A和K分别为尺度参数和为形状参数[20].
-
底层大气的层结状态和地表的粗糙度极大地决定了近地层风速垂直分布情况.若大气层结为中性，无论采用幂指数还是对数方程来描述风速的垂直廓线，都可以获得较好的效果.我国《建筑结构设计规范》中推荐采用幂指数公式，而根据相关研究所得实地测量数据，经多次检验结果表明，相对于对数公式，幂指数公式能够更加精确地拟合风速的垂直廓线，其表达式如下：
式中：V1和V2分别为高度Z1和Z2处的风速(m/s)；此处Z1取10 m高度；α为风切变指数，其值的大小可以体现风速垂直切变的强度[13].
1.1. 数据来源
1.2. 计算方法
1.2.1. 空气密度
1.2.2. 年平均风速和年平均风功率密度
1.2.3. 风能资源长年代订正
1.2.4. 风向和风能密度分布
1.2.5. 风向频率、风能方向频率和风速平率计算
1.2.6. 风频曲线及威布尔参数计算
1.2.7. 风速垂直切变
-
据布拖县气象站数据所得，布拖县1981-2010年风速为2.1 m/s，其平均风速变化情况如图 1所示.近30年来，布拖县年平均风速呈下降趋势，2010年较1981年相比，年均风速降低41%，与四川省大多数气象站近年来的变化趋势一致.这可能与气象站周围建筑物逐年增多，城市化进程加快有关[12].故在新建风电场的局地区域设立测风塔进行观测十分必要，因为长年代气象站所测风速并不能很好地反应局地风资源情况.
由布拖县1981-2010年平均风速月变化情况可知，1-4月风速较大，3月最大，风速为3.2 m/s；7-9月风速较小，8月最小，为1.4 m/s.
-
图 2所示为布拖县风向频率分布情况.根据布拖县气象站多年相关气象资料统计计算，该地区主要风向出现在ENE，出现的频率为11%；在WSW方向次数也较多，出现频率都为9%，该处静风频率为27%.
-
A测风塔平均风速和风功率密度月变化情况如图 3所示，B测风塔平均风速和风功率密度数值和月变化情况与A塔一致.测风塔A塔70 m和10 m高度年平均风速分别为7.3，6.7 m/s，年均风功率密度分别为338，271 W/m2；B塔70 m和10 m高度年平均风速分别为8.8，9.1 m/s，年均风功率密度分别为511，586 W/m2. A和B测风塔在10 m高度和70 m高度的平均风速和风功率密度月变化情况高度相似，在4-6月为低值月份，在11-次年3月为高值月份，最大值都为2月，A塔70 m高度该月平均风速和风功率密度为10.6 m/s和804 W/m2，B塔为13 m/s和1121 W/m2. A塔和B塔的平均风速和风功率密度均呈现出冬季较大，夏季较小的规律.
-
由图 4可知A测风塔平均风速和风功率密度日变化情况，B测风塔平均风速和风功率密度数值和日变化情况和A塔一致. A测风塔70 m和10 m高度平均风速和风功率密度日变化规律高度相似，均呈现出了白天较小，夜间较大的规律；最小值出现在12时，此时的70 m平均风速和风功率密度分别为6.2 m/s和208 W/m2，10 m平均风速和风功率密度分别为5.7 m/s和171 W/m2；最大值出现在晚上20时，此时的70 m平均风速和风功率密度分别为8.9 m/s和489 W/m2，10 m平均风速和风功率密度分别为8.2 m/s和400 W/m2. B测风塔的平均风速和风功率密度日变化情况不大，较为稳定，但风速和风功率密度间日变化规律一致性较差，且无明显规律.平均风速最大值出现在20时和21时，70 m和10 m分别为9.1，9.4 m/s；平均风功率密度最大值出现在9时，70 m和10 m分别为561，650 W/m2.
-
风电机启动风速至破坏风速间的风速成为有效风速，范围通常为3~25 m/s.测风塔等级风速频率分布如图 5所示.由图可知，A塔70 m和10 m的有效风速出现频率之和分别达到了86%和84%，风速主要分布区间在1~16 m/s，70 m和10 m的风速频率出现最高的速段分别为5，4 m/s；B塔70 m和10 m的有效风速出现频率分别达到89%，风速均主要分布区间在1~18 m/s区间，为70 m风速频率出现最高的速段，10 m风速频率最高速段为11 m/s.
-
A和B风塔70 m高度风向和风能密度方向分布如图 6所示. A，B两塔风向和风能密度都分布在西南方向.其中，A塔70 m和10 m主导风向均为SSW，频率均为26%；风能密度方向均为SW，频率分别为40%和66%；B塔70 m和10 m主导风向为SW，在该风向出现的频率均为41%，风能密度主要方向与主导风向同向为SW，在改方向出现的频率分别均为66%.不同高度层的风向和风能密度方向都较为集中，利于风能资源的有效开发和利用.
-
A和B测风塔70 m高度的威布尔尺度和形状参数如表 1所示，风速威布尔分布曲线如图 7所示，A和B测风塔10 m高度分布曲线与70 m高度一致.可以看出A塔风速威布尔分布中，出现频率最高风速段为6 m/s；B塔出现频率最高风速段则为7 m/s. A，B两塔就与威布尔函数分布曲线吻合度而言，70 m高度比10 m高度好.
-
据幂指数方法，计算得A和B测风塔风切变指数如图 8所示. A塔模拟风速随着高度增加而增大，而B塔则出现了负切变，这可能与测风塔所在的地形复杂有关.
2.1. 布拖气象站风参数背景分析
2.1.1. 风速年际和年内变化分析
2.1.2. 风向频率
2.2. 观测塔长年代风能参数分析
2.2.1. 风速和风功率密度年内变化
2.2.2. 风速和风功率密度小时变化
2.2.3. 风速频率分布
2.2.4. 风向频率和风能密度比例分布
2.2.5. 风频曲线及威布尔参数
2.2.6. 风速垂直切变
-
1) 经分析，本研究中布拖县气象站的多年观测值结果与布拖县测风塔实际测量的平均风速和风向值差异较大.结果反映了布拖县风能资源较强的局地性，这与其他地区实测风能资源评估结果类似.研究结果也说明，在布拖县地区进行风能资源开发，于修建风电场之前，必须先设立测风塔，采用实地观测资料.此举可较大避免风机布设和修建场址的不合理性[20].
2) 布拖县的A塔10 m、B塔的10 m和70 m高度风能资源有很好的开发潜力.根据“风电场风能资源评估方法”[9]中的标准：布拖县A塔在10 m高度风能资源应用于并网风能发电达到“好”的等级，B塔在70 m高度和10 m高度风能资源分别达到“好”和“很好”.且A和B测风塔区域风能在时间和方向的分布上都较为稳定和集中，利于风能资源开发.
3) 布拖县属亚热带湿润季风气候，夏季高温多雨，冬季温暖湿润.年平均气温10.1 ℃，1月均温1.3 ℃，7月均温17.3 ℃.本研究区域位于地形较为平坦的高海拔山脊地带，在该处设立风电场，不仅附近公路顺畅，便于运输修建材料，还不占用耕地，无需增加赔偿成本；冬无严寒，夏无酷暑，气候温和，无极端气象灾害，便于设备安装和风电并网，对修建风电场和建成后的运行维护有利.挑选具体风电场址修建时，可在风能资源较好的地区，在复杂山地地形中挑选较为平坦的山脊地段沿线布设，充分利用该地区的风能资源，正常稳定运行风电机，减少风机间的相互影响，提高风电场的电量指标[20].
4) 本研究受资料所限，仅对位于凉山州中东部的布拖县测风点区域进行分析.凉山州还存在大量地形和公路条件较好，且风能资源有潜力的区域，可继续对这些区域进行详实的实地调查实验，竖塔观测，结合气象资料，对四川省风能潜力进行更加深入的挖掘.