CN107608008B - 一种基于广义大气浑浊度的晴空时段的检测方法 - Google Patents
一种基于广义大气浑浊度的晴空时段的检测方法 Download PDFInfo
- Publication number
- CN107608008B CN107608008B CN201710834253.2A CN201710834253A CN107608008B CN 107608008 B CN107608008 B CN 107608008B CN 201710834253 A CN201710834253 A CN 201710834253A CN 107608008 B CN107608008 B CN 107608008B
- Authority
- CN
- China
- Prior art keywords
- turbidity
- clear
- atmospheric
- current
- time
- 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.)
- Active
Links
Images
Landscapes
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
Description
技术领域
本发明涉及的是一种检测方法,具体地说是晴空状态的检测方法。
背景技术
随着大规模光伏(PV)发电***的逐渐并入电网运行,光伏发电对电网的影响愈加增强。光伏发电***依靠光伏阵列将太阳能转化为电能,光伏发电***瞬时输出电功率的大小主要与光伏阵列所在地理位置的瞬时太阳辐照度相关。这里所谓的太阳辐照度是指太阳辐射经过大气传输过程中的吸收、散射、反射等作用后,到达固体地球表面上单位面积的辐射能量。
太阳辐射在从太阳到地表的传输过程中,会受到各种传输介质的影响,其某些特定的波长成分将会衰减,其中影响衰减的主要介质为大气成分和云层。大气中各种成分都或多或少地会对太阳辐射存在一定的衰减和散射作用。通常,在晴空条件下可以只采用大气浑浊度来描述大气各种成分对于太阳辐照度的综合衰减程度。这里所说的晴空状态被定义为:在整个可视天空范围内不存在或极少地存在云层的天空状态。而在非晴空条件下,除了大气浑浊度外,通常还要采用晴空因子来描述云层等非大气成分太阳辐照度的衰减作用。一般来说,若在一天中,某天空区域的大气浑浊度的变化较小,则表示该天空区域的大气环境在这一天中,基本处于一个近乎平稳的状态。
因为相比于大气中的各类成分,云层对辐照度的衰减程度更为严重,所以,云层是晴空因子主要影响因素。有研究表明:由于云层的遮挡作用,严重时,甚至可以使PV发电***的瞬时输出功率下降80%以上。此外,因云层的消长和移动等变化,亦会加剧作用于光伏阵列的太阳辐照度的波动性。因此,光伏发电***相对于风电等其它可再生能源发电***来说,其输出功率会呈现出更为明显的波动性和短时突变性。非晴空条件下,若透过云层到达光伏板表面的辐照度过低或辐照度波动过大,都将严重降低PV发电***的瞬时出力功率,进而对PV注入电网的有功功率和电网频率造成不利影响。对于这种情况,现行的电力网***应对的策略大多是采取热备用容量紧急投入的处理方式;若这种情况发生较为频繁,并联电网甚至会采取暂时切除PV的调度策略。
相比非晴空状态,当A地T时段处于晴空状态时,光伏阵列将接收到最大的太阳辐照,而且辐照度的变化也会较为平缓。与此对应,A地T时段内,处于晴空状态时的PV发电***相对于非晴空情况,其输出电能的电量也更多,而且波动次数更少。这些都有利于维持PV***持续而稳定地并网供电。
由上述可见:在光伏发电站前期选址和规划设计阶段,不仅需要利用该地区的历史辐照度数据及相关地理数据,年平均晴空时长和晴空辐照度时序分布数据也将是极为重要的寻址和规划设计依据。晴空(时段)检测则可以帮助设计人员在几个待选PV电站位置中,计算每个待选站点的晴空时长,确定其中太阳能资源最优的地点,进而为光伏电站选址的提供重要的依据参数。因此,晴空阶段的检测对PV发电***(孤立型或并网型)具有重要的意义。晴空(时段)检测还有其它的应用场合,例如为光伏阵列重组(重构)提供重组启动信号。晴空检测对于光伏发电***优化运行(包括PV阵列重构)、提高***运行稳定性具有重要的意义,可为优化光伏发电***前期规划和后期运行提供了帮助。
发明内容
本发明的目的在于提供将云层遮挡的影响纳入到广义大气浑浊度中,利用当前时刻的直射辐照度推算广义大气浑浊度的一种基于广义大气浑浊度的晴空时段的检测方法。
本发明的目的是这样实现的:
本发明一种基于广义大气浑浊度的晴空时段的检测方法,其特征是:
a通过直射辐照度传感器检测当前时刻的直射辐照度I(tk),检测的时间间隔为Δt;
c计算晴空大气浑浊度;
d判断当前是否为晴空状态。
本发明还可以包括:
(1)将时刻tk和数据采集地点的地理信息,包括海拔h、经度、纬度输入到SPA算法中,得到日地相对距离d和日地平均距离d0的比值、太阳天顶角Z(tk)和太阳方位角aZ,计算地外辐照度I0:
I0=Ie·(d0/d)2
其中,Ie为太阳辐照度常数,这里为1361.2W/m2;
(2)由上一步中所得的太阳天顶角Z(tk),求相对光学大气质量m(tk):
(3)根据测量数据地点的海拔h,确定用于计算所测地点的海拔修正参数b:
(1)检测针对当天,在当前已推算出的大气浑浊度序列是否存在连续的方差即小于1且数据量超过10的数据段。其中,方差的计算公式如下:
(2)当已计算出的大气浑浊度序列中存在数量大于10且方差小于1的子序列时,这部分子序列即为或近似为晴空大气浑浊度;在上述情况下,通过分段拟合差值法来获得当前时刻的晴空大气浑浊度具体做法为:利用前一天的晴空大气浑浊度序列和当日已获得的小段晴空大气浑浊度数据在相应时间段进行拟合。处理的详细过程是:将之前检测出的方差小于1的大气浑浊度数据段与前一天中相应时间段的晴空大气浑浊度相减,得到两者相应时间段的差值数据序列;再利用拟合法对与的差值数据进行拟合,即可得到与该时间段所有数据点的偏差总和最小的一条拟合直线;
直射辐照度的检测时间为1min,该天气状态即为这1min的天气状态,日晴空时长判断即为每天判定为晴空的时间长度相加,晴空辐照度时序分布即依据晴空检测结果对检测辐照度进行依次标记得出。
本发明的优势在于:在光伏发电***规划与设计中,希望选取太阳能资源储量最大且遮挡较少的地区为光伏电站待建厂址。本发明将云层遮挡的影响纳入到广义大气浑浊度中,利用当前时刻的直射辐照度推算广义大气浑浊度,并其与晴空大气浑浊度进行比较,最终可得出当前的天空状态(是否晴空)。本发明提出的新方法既可以减少晴空检测所需要的数据量,又提高了检测的时效性。同时,运用本方法可以减少因增设检测设备所需的晴空检测***的建设成本。
附图说明
图1为本发明的流程图;
图2为相邻两天晴空大气浑浊度序列及拟合差值区域图;
图3为分段修正法原理图。
具体实施方式
下面结合附图举例对本发明做更详细地描述:
结合图1-3,大气浑浊度原本仅表示在理想晴空条件下的大气浑浊情况(后文简称为“狭义大气浑浊度”或“晴空大气浑浊度”),本发明为了同时考虑云层和大气成分对太阳光辐射的衰减作用,将云层对地面辐照度的影响也并入大气浑浊度来表征,这里定义为“广义大气浑浊度”(下文将简称为大气浑浊度,并记为)。在本发明中,将原“狭义大气浑浊度”称之为“晴空大气浑浊度”,并记为在晴空条件下,与的值相同,均表示大气成分对辐照度的衰减;在非晴空(如多云、阴天等)天气条件下,的数值体现的是大气成分、云层及其他影响因素共同作用下的辐照度衰减程度,而则表示除云层外的其他因素对辐照度的影响程度。
本发明解决该技术问题所采用的技术方案简单表述为:
a、读取当前tk时刻的直射辐照度数据I(tk)
(1)将时刻t和数据采集地点的地理信息(如海拔、经度、纬度等)输入到太阳位置算法(Solar Position Algorithm,SPA)中,得到日地相对距离d和日地平均距离d0的比值、太阳天顶角Z和太阳方位角aZ,并计算地外辐照度I0;
(2)由上一步中得到的太阳天顶角,求取相对光学大气质量m;
(3)根据测量数据地点的海拔h,计算海拔修正参数b;
因为与大气质量m及其他参数的相关性均较小,所以只要Ineichen模型中的参量的量值不变,晴空大气浑浊度就基本维持在一个稳定的状态(数值)。但是,前后两天之间的晴空大气浑浊度可能会存在一定的差距。这是因为:一天中的高空大气环境较为稳定,但经历夜晚后,高空大气环境会发生一定的变化。因此,相邻两天的大气浑浊度曲线会存在一定的差距(如图2)。因此,在本发明中,将基于前一天的晴空大气浑浊度序列,选取分段线性拟合差值的方法或采用历史数据拟合机构预测数据的方法来得出当天的晴空大气浑浊度序列。
d、判断当前时刻是否为晴空状态。
将当前大气浑浊度与晴空大气浑浊度进行对比分析,可得出当前天空状态是否为晴空。
本发明是一种针对光伏发电***的规划和运行过程中需要获取晴空参数和晴空时长的现实需求,基于广义大气浑浊度序列,通过计算获取晴空时长和判定当前实时天空状态的方法。其具体实施包括如下步骤:
a、通过直射辐照度传感器检测当前时刻的直射辐照度I(tk)。检测的时间间隔Δt可为几分钟,这里取1min。
按照如下步骤进行广义大气浑浊度的计算:
(1)将时刻tk和数据采集地点的地理信息,如海拔h、经度、纬度等输入到SPA算法中。可以得到日地相对距离d和日地平均距离d0的比值、太阳天顶角Z(tk)和太阳方位角aZ。然后,按照(2)式计算地外辐照度I0。
I0=Ie·(d0/d)2 (2)
其中,Ie为太阳辐照度常数,这里为1361.2W/m2.
(2)由上一步中所得的太阳天顶角Z(tk),求相对光学大气质量m(tk)
(3)根据测量数据地点的海拔h,确定用于计算所测地点的海拔修正参数b:
c、晴空大气浑浊度计算
晴空大气浑浊度主要受高空中的大气成分变化的影响。因为依靠单个地点的检测手段对高空气体进行高质量的检测的难度较大,所以可以依靠专业的测试机构发布的测试数据,例如SoDa Service***。SoDa Service***依托其在世界各地建立的天气数据服务站所上传的测量数据,对世界各地的晴空大气浑浊度进行记录并协调预测,并定时发布其记录和预测结果。但是这些专业的测试机构所预测并公开的晴空大气浑浊度的往往两点之间的时间间隔较大,对大部分地区来说,往往仅能提供该地区的每天的晴空大气浑浊度的日平均值。
因此本发明采用如附图3中的分段修正法基于前一天的晴空大气浑浊度序列,通过拟合法获取当前及当日的晴空大气浑浊度而在当日已计算出的的数据量小于10或者根据当前已检测出的数据仍无法利用分段修正法推测出当前时刻的的情况下,采用专业的测试机构,如SoDa Service***,在前一天所发布的对本日的的预测数据。
(1)检测针对当天,在当前已推算出的大气浑浊度序列是否存在连续的方差(公式(6))小于D(这里D取1)且数据量超过10的数据段。如不存在,代表此时无法采用分段修正法计算晴空大气浑浊度,此时将采用采用专业的测试机构预测数据,则依据公式(7)对前一天的日照时间内的晴空大气浑浊度序列取平均值(t1~tm为位于日照时间内的数据采集时间点)。
求出与采用专业的测试机构对晴空大气浑浊度的预测数据之间的差值,并将其与前一天日照时间内的晴空大气浑浊度序列的每个值相加,即可得到基于采用专业的测试机构数据的当前时刻的晴空大气浑浊度其计算公式如式(8)(t1~tm为位于日照时间内的数据采集时间点):
(2)如图3,当已计算出的大气浑浊度序列中存在数量大于10且方差小于1的子序列时,这部分子序列即为或近似为晴空大气浑浊度。在上述情况下,通过分段拟合差值法来获得当前时刻的晴空大气浑浊度具体做法为:利用前一天的晴空大气浑浊度序列和当日已获得的小段晴空大气浑浊度数据在相应时间段进行拟合,利用拟合法得出的数据作为当前时刻的晴空大气浑浊度的预测值。具体过程如图3,将之前检测出的方差小于1的大气浑浊度数据段与前一天中相应时间段的晴空大气浑浊度相减,得到两者相应时间段的差值数据序列。相关示意图如图3所示,图中以零点为原点,x轴为时间轴,y轴为大气浑浊度。利用拟合法(如最小二乘法)对与的差值数据进行拟合即可得到与该段所有数据点的偏差和最小的一条拟合直线。
(3)步骤(2)中得到的拟合曲线代表了图3中三角形阴影区域的边界,整个阴影区域表示前一天的晴空大气浑浊度和当日晴空大气浑浊度的差将当前时刻tk带入到步骤(2)中拟合出的直线中可得出当前时刻两大气浑浊度的拟合差值所构成的区域。将拟合差值与前一天中对应当前时刻的晴空大气浑浊度相加,即可得出当前时刻的晴空大气浑浊度
d、判断当前是否为晴空状态。
因为直射辐照度的检测时间为1min,所以该天气状态即为这1min的天气状态,日晴空时长判断即为每天判定为晴空的时间长度相加,晴空辐照度时序分布即依据晴空检测结果对检测辐照度进行依次标记得出。
Claims (1)
1.一种基于广义大气浑浊度的晴空时段的检测方法,其特征是,包括以下步骤:
a通过直射辐照度传感器检测当前时刻的直射辐照度I(tk),检测的时间间隔为Δt;
d判断当前是否为晴空状态;
(Ⅰ)将时刻tk和测量数据地点的地理信息,包括海拔h、经度、纬度输入到太阳位置算法中,得到日地相对距离d和日地平均距离d0的比值、太阳天顶角Z(tk)和太阳方位角aZ,由日地相对距离d和日地平均距离d0的比值计算地外辐照度I0:
I0=Ie·(d0/d)2
其中,Ie为太阳辐照度常数,这里为1361.2W/m2;
(Ⅱ)由上一步中所得的太阳天顶角Z(tk),求相对光学大气质量m(tk):
(Ⅲ)根据测量数据地点的海拔h,确定用于计算所测地点的海拔修正参数b:
(2)当已计算出的大气浑浊度序列中存在数量大于10且方差小于1的子序列时,这部分子序列即为可以表征晴空大气浑浊度的数据段,在上述情况下,通过分段拟合差值法来获得当前时刻的晴空大气浑浊度具体做法为:利用前一天的晴空大气浑浊度序列和当日已获得的小段晴空大气浑浊度数据在相应时间段进行拟合,处理的详细过程是:将之前检测出的方差小于1的大气浑浊度数据段与前一天中相应时间段的晴空大气浑浊度相减,得到两者相应时间段的差值数据序列;利用拟合法对与的差值数据进行拟合即可得到与该时间段所有数据点的偏差总和最小的一条拟合直线;
步骤d具体为:
直射辐照度的检测时间为1min,该天气状态即为这1min的天气状态,日晴空时长判断即为每天判定为晴空的时间长度相加,晴空辐照度时序分布即依据晴空检测结果对检测辐照度进行依次标记得出。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710834253.2A CN107608008B (zh) | 2017-09-15 | 2017-09-15 | 一种基于广义大气浑浊度的晴空时段的检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710834253.2A CN107608008B (zh) | 2017-09-15 | 2017-09-15 | 一种基于广义大气浑浊度的晴空时段的检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107608008A CN107608008A (zh) | 2018-01-19 |
CN107608008B true CN107608008B (zh) | 2020-09-25 |
Family
ID=61060300
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710834253.2A Active CN107608008B (zh) | 2017-09-15 | 2017-09-15 | 一种基于广义大气浑浊度的晴空时段的检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107608008B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109344491A (zh) * | 2018-09-27 | 2019-02-15 | 河北工业大学 | 一种考虑天气状态变化和云层遮挡的太阳辐照度建模方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4463388B2 (ja) * | 2000-06-29 | 2010-05-19 | 小糸工業株式会社 | 視認状況測定装置 |
CN104181612B (zh) * | 2014-08-13 | 2016-03-30 | 中国人民解放军理工大学 | 一种红外与激光联合的地基测云方法 |
CN105005937B (zh) * | 2015-04-28 | 2019-02-15 | 国家电网公司 | 一种基于晴空指数的光伏电站出力时序模拟方法 |
-
2017
- 2017-09-15 CN CN201710834253.2A patent/CN107608008B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN107608008A (zh) | 2018-01-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Cebecauer et al. | Uncertainty sources in satellite-derived direct normal irradiance: how can prediction accuracy be improved globally | |
Marques Filho et al. | Global, diffuse and direct solar radiation at the surface in the city of Rio de Janeiro: Observational characterization and empirical modeling | |
Zhao et al. | Optimal PV panel tilt angle based on solar radiation prediction | |
JP2013529051A (ja) | 太陽光発電予測システム並びに方法 | |
CN105335560B (zh) | 光伏发电功率波动性及其自动发电控制备用需求计算方法 | |
CN109543323B (zh) | 一种基于逐时气象数据的光伏***发电量的检测方法 | |
Kam et al. | Comparative Weibull distribution methods for reliable global solar irradiance assessment in France areas | |
CN114019579A (zh) | 高时空分辨率近地表空气温度重构方法、***、设备 | |
CN110781458B (zh) | 一种基于混合回归模型预测地表太阳辐照度的方法 | |
Lyden et al. | Modelling, parameter estimation and assessment of partial shading conditions of photovoltaic modules | |
Marzo et al. | Solar extinction map in Chile for applications in solar power tower plants, comparison with other places from sunbelt and impact on LCOE | |
Malvoni et al. | Improvements in the predictions for the photovoltaic system performance of the Mediterranean regions | |
Mahmud et al. | Solar energy resource assessment of the geba catchment, Northern Ethiopia | |
Huang et al. | Predicting intra-hour variability of solar irradiance using hourly local weather forecasts | |
KR20160078920A (ko) | 목표 지점의 일사량 계산 방법 및 장치 | |
CN109599899B (zh) | 一种新能源运行模拟边界条件的设定方法 | |
Benatiallah et al. | An empirical model for estimating solar radiation in the Algerian Sahara | |
CN107608008B (zh) | 一种基于广义大气浑浊度的晴空时段的检测方法 | |
Yao et al. | A new modified method of all-sky radiance distribution based on the principle of photothermal integration | |
Jayaraman et al. | Temporal and spatial variations of the solar radiation observed in Singapore | |
CN116247663A (zh) | 光伏功率突变预测方法 | |
Wang et al. | Method for short-term photovoltaic generation power prediction base on weather patterns | |
Eddine et al. | Estimation model for global solar radiation in Lebanon | |
Rincón et al. | Assessment of short-term irradiance forecasting based on post-processing tools applied on WRF meteorological simulations | |
Anadol et al. | The effect of snowfall and icing on the sustainability of the power output of a grid-connected photovoltaic system in Konya, Turkey |
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 |