一种基于广义大气浑浊度的晴空时段的检测方法与流程

本发明涉及的是一种检测方法，具体地说是晴空状态的检测方法。

背景技术：

随着大规模光伏(pv)发电系统的逐渐并入电网运行，光伏发电对电网的影响愈加增强。光伏发电系统依靠光伏阵列将太阳能转化为电能，光伏发电系统瞬时输出电功率的大小主要与光伏阵列所在地理位置的瞬时太阳辐照度相关。这里所谓的太阳辐照度是指太阳辐射经过大气传输过程中的吸收、散射、反射等作用后，到达固体地球表面上单位面积的辐射能量。

太阳辐射在从太阳到地表的传输过程中，会受到各种传输介质的影响，其某些特定的波长成分将会衰减，其中影响衰减的主要介质为大气成分和云层。大气中各种成分都或多或少地会对太阳辐射存在一定的衰减和散射作用。通常，在晴空条件下可以只采用大气浑浊度来描述大气各种成分对于太阳辐照度的综合衰减程度。这里所说的晴空状态被定义为：在整个可视天空范围内不存在或极少地存在云层的天空状态。而在非晴空条件下，除了大气浑浊度外，通常还要采用晴空因子来描述云层等非大气成分太阳辐照度的衰减作用。一般来说，若在一天中，某天空区域的大气浑浊度的变化较小，则表示该天空区域的大气环境在这一天中，基本处于一个近乎平稳的状态。

因为相比于大气中的各类成分，云层对辐照度的衰减程度更为严重，所以，云层是晴空因子主要影响因素。有研究表明：由于云层的遮挡作用，严重时，甚至可以使pv发电系统的瞬时输出功率下降80％以上。此外，因云层的消长和移动等变化，亦会加剧作用于光伏阵列的太阳辐照度的波动性。因此，光伏发电系统相对于风电等其它可再生能源发电系统来说，其输出功率会呈现出更为明显的波动性和短时突变性。非晴空条件下，若透过云层到达光伏板表面的辐照度过低或辐照度波动过大，都将严重降低pv发电系统的瞬时出力功率，进而对pv注入电网的有功功率和电网频率造成不利影响。对于这种情况，现行的电力网系统应对的策略大多是采取热备用容量紧急投入的处理方式；若这种情况发生较为频繁，并联电网甚至会采取暂时切除pv的调度策略。

相比非晴空状态，当a地t时段处于晴空状态时，光伏阵列将接收到最大的太阳辐照，而且辐照度的变化也会较为平缓。与此对应，a地t时段内，处于晴空状态时的pv发电系统相对于非晴空情况，其输出电能的电量也更多，而且波动次数更少。这些都有利于维持pv系统持续而稳定地并网供电。

由上述可见：在光伏发电站前期选址和规划设计阶段，不仅需要利用该地区的历史辐照度数据及相关地理数据，年平均晴空时长和晴空辐照度时序分布数据也将是极为重要的寻址和规划设计依据。晴空(时段)检测则可以帮助设计人员在几个待选pv电站位置中，计算每个待选站点的晴空时长，确定其中太阳能资源最优的地点，进而为光伏电站选址的提供重要的依据参数。因此，晴空阶段的检测对pv发电系统(孤立型或并网型)具有重要的意义。晴空(时段)检测还有其它的应用场合，例如为光伏阵列重组(重构)提供重组启动信号。晴空检测对于光伏发电系统优化运行(包括pv阵列重构)、提高系统运行稳定性具有重要的意义，可为优化光伏发电系统前期规划和后期运行提供了帮助。

技术实现要素：

本发明的目的在于提供将云层遮挡的影响纳入到广义大气浑浊度中，利用当前时刻的直射辐照度推算广义大气浑浊度的一种基于广义大气浑浊度的晴空时段的检测方法。

本发明的目的是这样实现的：

本发明一种基于广义大气浑浊度的晴空时段的检测方法，其特征是：

(1)通过直射辐照度传感器检测当前时刻的直射辐照度i(tk)，检测的时间间隔为δt；

(2)计算广义大气浑浊度

(3)计算晴空大气浑浊度；

(4)判断当前是否为晴空状态。

本发明还可以包括：

1、按照如下步骤进行广义大气浑浊度的计算：

(1)将时刻tk和数据采集地点的地理信息，包括海拔h、经度、纬度输入到spa算法中，得到日地相对距离d和日地平均距离d0的比值、太阳天顶角z(tk)和太阳方位角az，计算地外辐照度i0：

i0＝ie·(d0/d)²

其中，ie为太阳辐照度常数，这里为1361.2w/m²；

(2)由上一步中所得的太阳天顶角z(tk)，求相对光学大气质量m(tk)：

(3)根据测量数据地点的海拔h，确定用于计算所测地点的海拔修正参数b：

(4)根据以上所求结果和当前时刻tk的直射辐照度i(tk)，计算当前时刻tk的广义大气浑浊度

2、当前时刻tk的晴空大气浑浊度的计算按照以下步骤实施：

(1)检测针对当天，在当前已推算出的大气浑浊度序列是否存在连续的方差即小于d且数据量超过10的数据段：

如不存在，代表此时无法采用分段修正法计算晴空大气浑浊度，此时采用测试机构的预测数据，依据下式对前一天的日照时间内的晴空大气浑浊度序列取平均值t1～tm为位于日照时间内的数据采集时间点：

求出与测试机构对晴空大气浑浊度的预测数据之间的差值，并将其与前一天日照时间内的晴空大气混浊度序列的每个值相加，得到基于采用测试机构数据的当前时刻的晴空大气浑浊度其计算公式如下式，

(2)当已计算出的大气混浊度序列中存在数量大于10且方差小于1的子序列时，这部分子序列即为或近似为晴空大气浑浊度；在上述情况下，通过分段拟合差值法来获得当前时刻的晴空大气浑浊度具体做法为：利用前一天的晴空大气浑浊度序列和当日已获得的小段晴空大气浑浊度数据在相应时间段进行拟合，利用拟合法得出的数据作为当前时刻的晴空大气浑浊度的预测值；

(3)将当前时刻tk带入到步骤(2)中拟合出的直线中得出当前时刻两大气浑浊度的拟合差值所构成的区域，将拟合差值与前一天中对应当前时刻的晴空大气浑浊度相加，得出当前时刻的晴空大气浑浊度

(4)将检测为晴空状态的数据点的大气浑浊度纳入到下一时刻的对进行直线拟合的过程中，对拟合直线进行修正，准确化所需的拟合值，指导达到预设标准。

3、将当前广义大气浑浊度与步骤c中所得的晴空大气浑浊度进行比较，并采用如下分段函数对当前晴空与否进行判定：

直射辐照度的检测时间为1min，该天空状态即为这1min的天空状态，日晴空时长判断即为每天判定为晴空的时间长度相加，晴空辐照度时序分布即依据晴空检测结果对检测辐照度进行依次标记得出。

本发明的优势在于：在光伏发电系统规划与设计中，希望选取太阳能资源储量最大且遮挡较少的地区为光伏电站待建厂址。本发明将云层遮挡的影响纳入到广义大气浑浊度中，利用当前时刻的直射辐照度推算广义大气浑浊度，并其与晴空大气浑浊度进行比较，最终可得出当前的天空状态(是否晴空)。本发明提出的新方法既可以减少晴空检测所需要的数据量，又提高了检测的时效性。同时，运用本方法可以减少因增设检测设备所需的晴空检测系统的建设成本。

附图说明

图1为本发明的流程图；

图2为相邻两天晴空大气浑浊度序列及拟合差值区域图；

图3为分段修正法原理图。

具体实施方式

下面结合附图举例对本发明做更详细地描述：

结合图1-3，大气浑浊度原本仅表示在理想晴空条件下的大气浑浊情况(后文简称为“狭义大气浑浊度”或“晴空大气浑浊度”)，本发明为了同时考虑云层和大气成分对太阳光辐射的衰减作用，将云层对地面辐照度的影响也并入大气浑浊度来表征，这里定义为“广义大气浑浊度”(下文将简称为大气浑浊度，并记为)。在本发明中，将原“狭义大气浑浊度”称之为“晴空大气浑浊度”，并记为在晴空条件下，与的值相同，均表示大气成分对辐照度的衰减；在非晴空(如多云、阴天等)天气条件下，的数值体现的是大气成分、云层及其他影响因素共同作用下的辐照度衰减程度，而则表示除云层外的其他因素对辐照度的影响程度。

本发明提出一种基于广义大气浑浊度的晴空时段检测方法。需将当前广义大气浑浊度与晴空大气混浊度进行比较，从而判断当前天空是否处于晴空状态。

本发明解决该技术问题所采用的技术方案简单表述为：

a、读取当前tk时刻的直射辐照度数据i(tk)

b、通过晴空模型计算大气浑浊度所选晴空模型为(1)式所示的ineichen晴空模型，并按照如下步骤计算大气浑浊度：

(1)将时刻t和数据采集地点的地理信息(如海拔、经度、纬度等)输入到太阳位置算法(solarpositionalgorithm,spa)中，得到日地相对距离d和日地平均距离d0的比值、太阳天顶角z和太阳方位角az，并计算地外辐照度i0；

(2)由上一步中得到的太阳天顶角，求取相对光学大气质量m；

(3)根据测量数据地点的海拔h，计算海拔修正参数b；

(4)根据以上所求结果和当前时刻tk的直射辐照度i(tk)，计算当前时刻tk的广义大气浑浊度

c、晴空大气浑浊度的计算

因为与大气质量m及其他参数的相关性均较小，所以只要ineichen模型中的参量的量值不变，晴空大气浑浊度就基本维持在一个稳定的状态(数值)。但是，前后两天之间的晴空大气浑浊度可能会存在一定的差距。这是因为：一天中的高空大气环境较为稳定，但经历夜晚后，高空大气环境会发生一定的变化。因此，相邻两天的大气浑浊度曲线会存在一定的差距(如图2)。因此，在本发明中，将基于前一天的晴空大气浑浊度序列，选取分段线性拟合差值的方法或采用历史数据拟合机构预测数据的方法来得出当天的晴空大气浑浊度序列。

d、判断当前时刻是否为晴空状态。

将当前大气浑浊度与晴空大气浑浊度进行对比分析，可得出当前天空状态是否为晴空。

本发明是一种针对光伏发电系统的规划和运行过程中需要获取晴空参数和晴空时长的现实需求，基于广义大气浑浊度序列，通过计算获取晴空时长和判定当前实时天空状态的方法。其具体实施包括如下步骤：

a、通过直射辐照度传感器检测当前时刻的直射辐照度i(tk)。检测的时间间隔δt可为几分钟，这里取1min。

b、通过式(1)的晴空模型计算广义大气浑浊度

按照如下步骤进行广义大气浑浊度的计算：

(1)将时刻tk和数据采集地点的地理信息，如海拔h、经度、纬度等输入到spa算法中。可以得到日地相对距离d和日地平均距离d0的比值、太阳天顶角z(tk)和太阳方位角az。然后，按照(2)式计算地外辐照度i0。

i0＝ie·(d0/d)²(2)

其中，ie为太阳辐照度常数，这里为1361.2w/m².

(2)由上一步中所得的太阳天顶角z(tk)，求相对光学大气质量m(tk)

(3)根据测量数据地点的海拔h，确定用于计算所测地点的海拔修正参数b：

(4)根据以上所求结果和当前时刻tk的直射辐照度i(tk)，计算当前时刻tk的广义大气浑浊度

c、晴空大气浑浊度计算

晴空大气浑浊度主要受高空中的大气成分变化的影响。因为依靠单个地点的检测手段对高空气体进行高质量的检测的难度较大，所以可以依靠专业的测试机构发布的测试数据，例如sodaservice系统。sodaservice系统依托其在世界各地建立的天气数据服务站所上传的测量数据，对世界各地的晴空大气混浊度进行记录并协调预测，并定时发布其记录和预测结果。但是这些专业的测试机构所预测并公开的晴空大气浑浊度的往往两点之间的时间间隔较大，对大部分地区来说，往往仅能提供该地区的每天的晴空大气混浊度的日平均值。

因此本发明采用如附图3中的分段修正法基于前一天的晴空大气浑浊度序列，通过拟合法获取当前及当日的晴空大气浑浊度而在当日已计算出的的数据量小于10或者根据当前已检测出的数据仍无法利用分段修正法推测出当前时刻的的情况下，采用专业的测试机构，如sodaservice系统，在前一天所发布的对本日的的预测数据。

当前时刻tk的晴空大气浑浊度的计算可按照以下步骤实施：

(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)中拟合出的直线中可得出当前时刻两大气浑浊度的拟合差值所构成的区域。将拟合差值与前一天中对应当前时刻的晴空大气浑浊度相加，即可得出当前时刻的晴空大气浑浊度

(4)将检测为晴空状态的数据点的大气浑浊度也纳入到下一时刻的对进行直线拟合的过程中，对拟合直线进行不断地修正，不断准确化所需的拟合值。

d、判断当前是否为晴空状态。

将当前广义大气浑浊度与步骤c中所得的晴空大气浑浊度进行比较，并采用如下分段函数对当前晴空与否进行判定：

因为直射辐照度的检测时间为1min，所以该天空状态即为这1min的天空状态，日晴空时长判断即为每天判定为晴空的时间长度相加，晴空辐照度时序分布即依据晴空检测结果对检测辐照度进行依次标记得出。