CN115840740A - 一种用于光伏电站的太阳能资源缺测数据插补方法 - Google Patents

Abstract

本发明公开了一种用于光伏电站的太阳能资源缺测数据插补方法，结合《太阳能资源评估方法》(GB/T 37256—2019)规定，并充分考虑与缺测时刻最近的实测时刻数据及其太阳总辐照度、大气温度和大气湿度等关键因素，提出了一种系统的、合理的太阳能资源缺测数据插补方法，该方法计算出的数据更贴近实际情况，同时极大程度地提高了数据插补的效率、合理性和准确性；本发明还提出了太阳能资源数据监测系统，通过相应的数据采集卡、各类传感器、上位机程序和数据库，可以实时采集实际工程所需监测的太阳能资源数据，并将实测数据写入数据库，不仅可以用于第三方平台数据共享，而且方便后期工程人员进行数据分析处理。

Description

一种用于光伏电站的太阳能资源缺测数据插补方法

技术领域

本发明涉及太阳能资源数据监测及数据插补技术领域，具体涉及一种用于光伏电站的太阳能资源缺测数据插补方法。

背景技术

太阳能资源数据广泛应用于能源、建筑、气象、电力、农业等领域，尤其对于光伏电站而言，其太阳能资源数据的完整性和准确性对项目选址和后期发电量的计算十分重要。2013年年底国家能源局出台了《光伏电站暂行管理办法》，其中里面第12条明确规定光伏项目开发前期，必须要有现场实测数据，或者有代表性的实测数据。目前国内光伏电站在前期可行性研究报告的编写过程中使用的太阳能资源数据主要来源为当地气象部门或其它商业软件提供的数据。

在实际工程项目开发前期，应对项目当地的太阳能资源数据进行数据采集和存储，以方便工程人员进行后期数据分析。但许多时候由于监测设备的故障、测量精度或其它外部的影响，所监测的短期太阳能资源实测数据会存在不合理值和缺测值。《太阳能资源评估方法》(GB/T37256—2019)明确规定，需要从实测数据系列中选择与缺测时刻最近、天气现象相同、有实测数据的时刻进行数据插补，实际工程项目操作中缺少对缺测数据的插补运算方法。

专利号为“CN201410652759”的“缺测数据补全方法以及数据收集装置”，该方法在发生缺测时，对能量使用量的变动式样与所述发生缺测的测量位置的能量使用量的变动式样类似的测量位置的、与所述缺测的期间相当的期间的每个测量时刻间隔的能量使用量比率进行计算，并使用该计算出的能量使用量比率，计算所述缺测的期间的每个测量时刻间隔的能量使用量，补全所述缺测数据；该方法虽然能补全缺测数据，但是由于仅参考历史数据，没有结合缺测时刻的实际因素，所以计算出来的数据不够客观和准确。

发明内容

本发明的目的在于提供一种用于光伏电站的太阳能资源缺测数据插补方法，结合《太阳能资源评估方法》(GB/T37256—2019)规定并充分考虑与缺测时刻最近的实测时刻数据及天气现象等关键因素，提出了一种系统的、合理的太阳能资源缺测数据插补方法，极大程度地提高了数据插补的效率、合理性和准确性。

为解决上述技术问题，本发明采用如下的技术方案：一种用于光伏电站的太阳能资源缺测数据插补方法，包括以下步骤：

S110、确定数据总条目，具体为：从数据库中取出需要检测太阳能资源数据，依据设置的步长Ts计算太阳能资源数据的数据条目总量为n；

S120、对n条太阳能资源数据进行合理性检查，记录不合理的数据条目后将其中不合理的数据值置为空值，太阳能资源数据中有置空项该条数据即为缺测数据；

S130、对n条太阳能资源数据进行完整性检查，计算得到有效数据完整率r_sd，通过有效数据完整率r_sd判断是否需要数据插补计算，如果需要插补计算则继续判断是否满足插补计算要求，满足插补计算要求则进行下一步；

S140、对满足插补计算要求的n条太阳能资源数据进行缺测数据插补计算，具体步骤为：

a、定义数据时刻；

b、搜索第一个实测数据并记录第一个实测数据时刻为t_p；

c、搜索第一个缺测数据并记录第一个缺测数据时刻为t_f；

d、由第t_f向前搜索，离t_f最近的第一个实测时刻为t_i，同时由第t_f向后搜索，离t_f最近的第一个实测时刻为t_j；

e、依据t_f、t_i和t_j计算得出t_f时刻的数据，并将该数据插补到t_f时刻太阳能资源数据中；

f、计算插补完整率r_p，并根据r_p判断是否继续进行缺测数据插补计算，如r_p<100％，则依次进行c、d、e和f步骤，如r_p＝100％，则完成本次缺测数据插补计算。

前述用于光伏电站的太阳能资源缺测数据插补方法，其中所述步骤S110中太阳能资源数据的数据条目总量为n，其中每条太阳能资源数据包括的字段有太阳总辐照度、太阳直接辐照度、太阳散射辐照度、风速风向、大气压强、大气温度和大气湿度，结合规范可知数据写入步长Ts应该小于1小时，依据监测系统所设定，10min或者20min均可。

前述用于光伏电站的太阳能资源缺测数据插补方法，其中所述步骤S120中对太阳能资源数据进行合理性检查，是根据《太阳能资源评估方法》(GB/T37256—2019)规定，从气候学界限值、内部一致性、变化范围3个方面对n条太阳能资源数据的合理性进行检查，记录不合理的数据条目，将其中不合理的数据值置为空值，太阳能资源数据中有置空项该条数据即为缺测数据。

前述用于光伏电站的太阳能资源缺测数据插补方法，其中所述步骤S130具体为：

S131、数据完整性检查：进一步检查n条太阳能资源数据的完整性，根据《太阳能资源评估方法》(GB/T37256—2019)规定，太阳能资源数据的实测数据记录应包括至少1年的连续、完整数据，数据记录应至少包括小时值，小时值的数据完整率应不低于95％；即n条太阳能资源数据的有效数据完整率应不低于95％，有效数据完整率r_sd的定义如下：

式中，N₀为预期记录的数据总条目，这里为n，N₁为没有太阳直接辐照度的数据条目，N₂为不合理的数据条目，N₁+N₂＝m；

具体为，N₀是太阳能资源数据所预期记录的总数据条目；譬如太阳直接辐射度预期记录3年，Ts为1h，则N₀＝1x24x3x365。

S132、依据数据完整率判断是否满足数据插补条件：假定所缺测的太阳能资源要素为太阳直接辐照度，DHI代表太阳直接辐照度；

首先判断是否满足数据插补要求，若太阳直接辐照度数据的完整率r_sd＝100％，则无需进行插补；若95％≤r_sd<100％，则需要对其进行数据插补；若r_sd<95％，则判断为实测数据异常不满足数据插补要求。

前述用于光伏电站的太阳能资源缺测数据插补方法，其中所述步骤S140具体为：

S141、定义数据时刻：将存在太阳直接辐照度的第一条太阳能资源数据对应起始时刻为t₁，则最后一条，即第n条数据对应结束时刻为t_n，两条相邻数据之间的时间间隔为步长Ts，n条数据里至少应包括1条实测数据；

S142、搜索第一个实测数据，对n条太阳能资源数据自前向后进行搜索，查询并记录第一条太阳能直接辐照度不为空的太阳能资源数据的时刻，即第一条太阳直接辐照度的实测数据时刻，见下式，

DHI(t_p)≠NULL t₁≤t_p≤t_n

式中，t_p为第一个实测数据时刻；

S143、搜索第一个缺测数据，对记录有太阳能直接辐照度的n条数据自前向后进行搜索，查询并记录第一条太阳能直接辐照度为空的太阳能资源数据的时刻，即第一条太阳直接辐照度的缺测数据时刻，

DHI(t_f)＝NULL t₁≤t_f≤t_n

式中，t_f为第一个缺测数据时刻；

S144、由步骤S143获取第一个缺测数据时刻t_f后，若向前搜索存在太阳能直接辐照度不为空的时刻，则记为t_i，若向后搜索存在太阳能直接辐照度不为空的时刻，则记为t_j，需要满足：

S145、依据t_f、t_i和t_j计算得出t_f时刻的太阳直接辐照度数据，并将该数据插补到t_f时刻太阳能资源数据的太阳直接辐照度数据项中；

S146、在步骤S145的数据插补成功后，将n条数据自前向后进行检索，当太阳直接辐照度数据为空值时进行记录，检索完成后所记录的空值条目总数即为m，计算所插补数据的插补完整率r_p，如r_p<100％，则依次进行S143、S144、S145和S146步骤，如r_p＝100％，则完成本次缺测数据插补计算；

插补完整率r_p的计算公式为：

式中，M₀为太阳能资源数据条目总量，这里为n条；M₁为太阳直接辐照度数据为空的条目总量，这里为m条。

前述用于光伏电站的太阳能资源缺测数据插补方法，其中所述步骤S145中计算缺测的太阳直接辐照度数据：

S1451、t_i存在，t_j不存在，根据t_i时刻的地外水平面太阳辐照度、t_i时刻的太阳直接辐照度和t_f时刻的地外水平面太阳辐照度，计算t_f时刻的太阳直接辐照度，具体计算如下：

上式中，EHI为地外水平面太阳辐照度，DHI为太阳直接辐照度，

将计算t_f时刻的太阳直接辐照度插补到t_f时刻太阳能资源数据对应数据项中；

S1452、t_i不存在，t_j存在，根据t_j时刻的地外水平面太阳辐照度、t_j时刻的太阳直接辐照度和t_f时刻的地外水平面太阳辐照度，计算t_f时刻的太阳直接辐照度，具体计算如下：

上式中，EHI为地外水平面太阳辐照度，DHI为太阳直接辐照度，

将计算t_f时刻的太阳直接辐照度插补到t_f时刻太阳能资源数据对应数据项中；

S1453、若t_i和t_j都存在，因为是对列数据自前向后搜索，则会出现三种情况分别为：

a、t_i时刻离t_f时刻更近，即t_f–t_i<t_j–t_f；此时按照步骤S1451中的计算公式对t_f时刻进行数据插补；

b、t_j时刻离t_f时刻更近，即t_f–t_i>t_j–t_f；此时按照步骤S1452中的计算公式对t_f时刻进行数据插补；

c、t_i时刻和t_j时刻离t_f时刻一样近，即t_f–t_i＝t_j–t_f，此时需要考虑天气现象，介于本发明所监测的太阳能资源数据，以太阳总辐照度、大气温度和大气湿度作为天气差异情况的衡量指标，根据相关工程经验并结合各指标对光伏电站的影响情况，设定其指标权重分别占比，太阳总辐照度占比为0.6、大气温度占比为0.1，大气湿度占比为0.3，分别计算t_i时刻和t_j时刻相较于t_j时刻的天气现象差异情况，具体计算如下：

式中，GHI为太阳总辐照度，T为大气温度，H为大气湿度，w代表天气现象差异情况，w数值越小表明两个时刻的天气现象更加相近；

进一步地，根据天气现象差异情况对t_f时刻进行数据插补，具体计算如下：

前述用于光伏电站的太阳能资源缺测数据插补方法，还包括有太阳能资源数据监测系统，所述太阳能资源数据监测系统包括光伏电站、数据采集装置和数据处理服务器，光伏电站内布置有数据采集装置，数据采集装置电性连接数据处理服务器。

前述用于光伏电站的太阳能资源缺测数据插补方法，其中所述数据采集装置包括总辐射传感器、直接辐射传感器、散射辐射传感器、风速风向传感器、气压传感器、温度传感器、湿度传感器和数据采集卡，总辐射传感器、直接辐射传感器、散射辐射传感器、风速风向传感器、气压传感器、温度传感器和湿度传感器均电性连接于数据采集卡；

数据采集卡采集光伏电站的重要太阳能资源数据，所述重要太阳能资源数据包括太阳总辐照度、太阳直接辐照度、太阳散射辐照度、风速风向、大气压强、大气温度和大气湿度。

前述用于光伏电站的太阳能资源缺测数据插补方法，其中所述数据处理服务器中配置有通讯接口、上位机LabVIEW程序和数据库，所述数据采集卡电性连接通讯接口，所述上位机LabVIEW程序将通讯接口接收到的重要太阳能资源数据进行处理，并将处理后的太阳能资源数据存储到数据库。

前述用于光伏电站的太阳能资源缺测数据插补方法，其中所述上位机LabVIEW程序包括用户登录模块、数据采集模块、数据存储模块、状态推送模块和报表打印模块；

所述用户登录模块具体为：用户登录模块基于数据库绑定用户表和角色表的动态信息权限，角色为电站管理员的用户赋予所有功能权限，角色为监测员的用户可查看LabVIEW前面板并实时监测光伏电站系统的太阳能资源数据，角色为工程师的用户可查看后台程序框图并对其功能进行修改完善；

所述数据采集模块具体为：采用“生产者——消费者”模型作为程序设计的主体结构，进一步设定数据采集卡的各项工作参数，包括物理通道，采样率、采样数和数据写入时间步长，对光伏电站所需监测的太阳能资源数据进行采集；

所述数据存储模块具体为：根据所设定的数据写入步长，将光伏电站的各项太阳能资源数据在LabVIEW前面板进行实时显示并自动写入数据库做后续分析；

所述状态推送模块具体为：通过对数据库中的历史数据进一步分析后可对当地太阳能资源进行评估，然后将评估结果实时推送给电站管理员；

所述报表打印模块：在用户需要查看资源报表时，该模块可以按逐时、逐日、逐月、逐年等条件将光伏电站的太阳能资源数据进行筛选统计并打印。

与现有技术相比，本发明的有益之处在于结合《太阳能资源评估方法》(GB/T37256—2019)规定，并充分考虑与缺测时刻最近的实测时刻数据及其太阳总辐照度、大气温度和大气湿度等关键因素，提出了一种系统的、合理的太阳能资源缺测数据插补方法，该方法计算出的数据更贴近实际情况，同时极大程度地提高了数据插补的效率、合理性和准确性；本发明还提出了太阳能资源数据监测系统，通过相应的数据采集卡、各类传感器、上位机程序和数据库，可以实时采集实际工程所需监测的太阳能资源数据，并将实测数据写入数据库，不仅可以用于第三方平台数据共享，而且方便后期工程人员进行数据分析处理。

附图说明

图1是本发明的缺测数据插补方法流程图；

图2是本发明中太阳能资源数据监测系统结构示意图；

下面结合附图和具体实施方式对本发明作进一步的说明。

具体实施方式

本发明的实施例1：一种用于光伏电站的太阳能资源缺测数据插补方法，包括以下步骤：

S110、确定数据总条目，具体为：从数据库中取出需要检测太阳能资源数据，依据设置的步长Ts计算太阳能资源数据的数据条目总量为n；

S120、对n条太阳能资源数据进行合理性检查，记录不合理的数据条目后将其中不合理的数据值置为空值，太阳能资源数据中有置空项该条数据即为缺测数据；

S130、对n条太阳能资源数据进行完整性检查，计算得到有效数据完整率r_sd，通过有效数据完整率r_sd判断是否需要数据插补计算，如果需要插补计算则继续判断是否满足插补计算要求，满足插补计算要求则进行下一步；

S140、对满足插补计算要求的n条太阳能资源数据进行缺测数据插补计算，具体步骤为：

a、定义数据时刻；

b、搜索第一个实测数据并记录第一个实测数据时刻为t_p；

c、搜索第一个缺测数据并记录第一个缺测数据时刻为t_f；

d、由第t_f向前搜索，离t_f最近的第一个实测时刻为t_i，同时由第t_f向后搜索，离t_f最近的第一个实测时刻为t_j；

e、依据t_f、t_i和t_j计算得出t_f时刻的数据，并将该数据插补到t_f时刻太阳能资源数据中；

f、计算插补完整率r_p，并根据r_p判断是否继续进行缺测数据插补计算，如r_p<100％，则依次进行c、d、e和f步骤，如r_p＝100％，则完成本次缺测数据插补计算。

本发明的实施例2：一种用于光伏电站的太阳能资源缺测数据插补方法，其特征在于，包括以下步骤：

S110、确定数据总条目，具体为：从数据库中取出需要检测太阳能资源数据，依据步长Ts计算太阳能资源数据的数据条目总量为n，其中每条太阳能资源数据包括的字段有太阳总辐照度、太阳直接辐照度、太阳散射辐照度、风速风向、大气压强、大气温度和大气湿度，结合规范可知数据写入步长Ts应该小于1小时，依据监测系统所设定，10min或者20min均可；

S120、对n条太阳能资源数据进行合理性检查，是根据《太阳能资源评估方法》(GB/T37256—2019)规定，从气候学界限值、内部一致性、变化范围3个方面对n条太阳能资源数据的合理性进行检查，记录不合理的数据条目，将其中不合理的数据值置为空值，太阳能资源数据中有置空项该条数据即为缺测数据；

S131、数据完整性检查：进一步检查n条太阳能资源数据的完整性，根据《太阳能资源评估方法》(GB/T37256—2019)规定，太阳能资源数据的实测数据记录应包括至少1年的连续、完整数据，数据记录应至少包括小时值，小时值的数据完整率应不低于95％；即n条太阳能资源数据的有效数据完整率应不低于95％，有效数据完整率r_sd的定义如下：

式中，N₀为预期记录的数据总条目，这里为n，N₁为没有太阳直接辐照度的数据条目，N₂为不合理的数据条目，N₁+N₂＝m；

具体为，N₀是太阳能资源数据所预期记录的总数据条目；譬如太阳直接辐射度预期记录3年，Ts为1h，则N₀＝1x24x3x365。

S132、依据数据完整率判断是否满足数据插补条件：假定所缺测的太阳能资源要素为太阳直接辐照度，DHI代表太阳直接辐照度；

首先判断是否满足数据插补要求，若太阳直接辐照度数据的完整率r_sd＝100％，则无需进行插补；若95％≤r_sd<100％，则需要对其进行数据插补；若r_sd<95％，则判断为实测数据异常不满足数据插补要求；

S140、对满足插补计算要求的n条太阳能资源数据进行缺测数据插补计算，具体步骤为：

S141、定义数据时刻：将存在太阳直接辐照度的第一条太阳能资源数据对应起始时刻为t₁，则最后一条，即第n条数据对应结束时刻为t_n，两条相邻数据之间的时间间隔为步长Ts，n条数据里至少应包括1条实测数据；

S142、搜索第一个实测数据，对n条太阳能资源数据自前向后进行搜索，查询并记录第一条太阳能直接辐照度不为空的太阳能资源数据的时刻，即第一条太阳直接辐照度的实测数据时刻，见下式，

DHI(t_p)≠NULL t₁≤t_p≤t_n

式中，t_p为第一个实测数据时刻；

S143、搜索第一个缺测数据，对记录有太阳能直接辐照度的n条数据自前向后进行搜索，查询并记录第一条太阳能直接辐照度为空的太阳能资源数据的时刻，即第一条太阳直接辐照度的缺测数据时刻，

DHI(t_f)＝NULL t₁≤t_f≤t_n

式中，t_f为第一个缺测数据时刻；

S144、由步骤S143获取第一个缺测数据时刻t_f后，若向前搜索存在太阳能直接辐照度不为空的时刻，则记为t_i，若向后搜索存在太阳能直接辐照度不为空的时刻，则记为t_j，需要满足：

S145、依据t_f、t_i和t_j计算得出t_f时刻的太阳直接辐照度数据，并将该数据插补到t_f时刻太阳能资源数据的太阳直接辐照度数据项中；

具体为：

S1451、t_i存在，t_j不存在，根据t_i时刻的地外水平面太阳辐照度、t_i时刻的太阳直接辐照度和t_f时刻的地外水平面太阳辐照度，计算t_f时刻的太阳直接辐照度，具体计算如下：

上式中，EHI为地外水平面太阳辐照度，DHI为太阳直接辐照度，

将计算t_f时刻的太阳直接辐照度插补到t_f时刻太阳能资源数据对应数据项中；

S1452、t_i不存在，t_j存在，根据t_j时刻的地外水平面太阳辐照度、t_j时刻的太阳直接辐照度和t_f时刻的地外水平面太阳辐照度，计算t_f时刻的太阳直接辐照度，具体计算如下：

上式中，EHI为地外水平面太阳辐照度，DHI为太阳直接辐照度，

将计算t_f时刻的太阳直接辐照度插补到t_f时刻太阳能资源数据对应数据项中；

S1453、若t_i和t_j都存在，因为是对列数据自前向后搜索，则会出现三种情况分别为：

a、t_i时刻离t_f时刻更近，即t_f–t_i<t_j–t_f；此时按照步骤S1451中的计算公式对t_f时刻进行数据插补；

b、t_j时刻离t_f时刻更近，即t_f–t_i>t_j–t_f；此时按照步骤S1452中的计算公式对t_f时刻进行数据插补；

c、t_i时刻和t_j时刻离t_f时刻一样近，即t_f–t_i＝t_j–t_f，此时需要考虑天气现象，介于本发明所监测的太阳能资源数据，以太阳总辐照度、大气温度和大气湿度作为天气差异情况的衡量指标，根据相关工程经验并结合各指标对光伏电站的影响情况，设定其指标权重分别占比，太阳总辐照度占比为0.6、大气温度占比为0.1，大气湿度占比为0.3，分别计算t_i时刻和t_j时刻相较于t_j时刻的天气现象差异情况，具体计算如下：

式中，GHI为太阳总辐照度，T为大气温度，H为大气湿度，w代表天气现象差异情况，w数值越小表明两个时刻的天气现象更加相近；

进一步地，根据天气现象差异情况对t_f时刻进行数据插补，具体计算如下：

S146、在步骤S145的数据插补成功后，将n条数据自前向后进行检索，当太阳直接辐照度数据为空值时进行记录，检索完成后所记录的空值条目总数即为m，计算所插补数据的插补完整率r_p，如r_p<100％，则依次进行S143、S144、S145和S146步骤，如r_p＝100％，则完成本次缺测数据插补计算；

插补完整率r_p的计算公式为：

式中，M₀为太阳能资源数据条目总量，这里为n条；M₁为太阳直接辐照度数据为空的条目总量，这里为m条。

本发明的实施例3：一种用于光伏电站的太阳能资源缺测数据插补方法，包括以下步骤：

步骤1：根据附图2所搭建的光伏电站太阳能资源数据监测系统，结合总辐射传感器、直接辐射传感器、散射辐射传感器、风速风向传感器、气压传感器、温度传感器和湿度传感器，实时采集某地光伏电站的太阳能资源数据，设定数据写入步长Ts为10min，即每隔10min将所采集数据做算术平均并写入预设数据库。

本实例选取了当地光伏电站某两天的实测数据，步长Ts为10min，每小时实测数据6条，数据总条数为2x24x6＝288条，其中缺测段实测数据如表1所示，缺测数据条目总量为8。

表1某地光伏电站某天的部分实测数据

注：太阳总辐照度(GHI)和太阳直接辐照度(DHI)单位均为W/m²，“---”表示缺测数据。

步骤2：数据合理性检查。根据《GB/T37256—2019太阳能资源评估方法》规定，从气候学界限值、内部一致性、变化范围3个方面对288条实测数据的合理性进行检查可知，不存在不合理数据。

步骤3：根据式(1)计算数据完整率r_sd。

步骤4：由步骤3计算结果可知r_sd≈2.7％，故需要对缺测段数据表1中太阳能直接辐照度的缺测数据进行插补运算。

步骤5：定义数据时刻，本次太阳能资源数据条目总量n为288条，所设定的数据写入步长Ts为10min，故定义数据时刻为t₁到t₂₈₈，缺测段数据表1中数据时刻为t_q1到t_q11。

步骤6：搜索第一个实测数据。并记录第一个实测数据时刻为t₁。

步骤7：搜索第一个缺测数据。通过表1搜索并记录第一个缺测数据时刻为t_q2，且DHI(t_q2)＝NULL。

步骤8：由步骤7获取第一个太阳能直接辐照度缺测数据时刻t_q2后，分别向前和向后同时搜索离t_q2最近的第一个实测时刻。t_q2时刻向上搜索最近的实测时刻为t_q1，向下搜索最近的实测时刻为t_q9。

步骤9：由步骤8可知缺测时刻t_q2向上和向下均存在实测时刻，但向上的实测时刻t_q1离缺测时刻t_q2更近，故根据计算式：

计算t_q2时刻的太阳直接辐照度，即：

其中地外水平面太阳辐照度(EHI)的定义如下：

上式中，φ表示当地纬度，单位为度；δ表示太阳赤纬，单位为度；w_t表示时角，单位为度；EDNI表示地外法向太阳辐照度，单位为W/m²；E₀为太阳常数，取1366.1W/m²；n_y表示积日，即日期在一年的序数。根据式(2)和式(3)计算可知DHI(t_q2)＝2.4862，根据该插补运算结果更新数据库，即将t_q2时刻的太阳能直接辐照度由空值插补为2.4862。

步骤10：再次计算插补完整率r_p，公式为：

计算缺测完整率发现r_p小于100％，故重新搜索新数据库中的缺测数据进行插补；

当t_q2到t_q7时刻的数据均插补完成后，按照自前向后的搜索原则，接下来需要对t_q8时刻的空值进行插补。

步骤11：对于t_q8时刻来说，此时t_q7时刻和t_q9时刻离t_q8时刻均相差一个数据步长且不为空值，故需要进一步考虑天气因素，根据式(4)计算t_q7时刻和t_q9时刻相较于t_q8时刻的天气差异情况结果如下：

查询表1代入相应数值计算可得w(t_q7-q8)＝0.4153、w(t_q9-q8)＝2.6457，根据所发明的方法，采用t_q7时刻的数据对t_q8时刻的空值进行插补，即：

根据上式计算可知DHI(t_q8)＝1.8090。

步骤12：完成上述步骤后，表1中还剩t_q10时刻的数据为空值，同理，对于t_q10时刻来说，此时t_q9时刻和t_q11时刻离t_q10时刻均相差一个数据步长且不为空值，其插补方法同上，通过式(4)计算可知w(t_q9-q10)＝0.9443、w(t_q11-q10)＝0.7826，故需要采用t_q11时刻的数据对t_q10时刻的空值进行插补，本实例不再赘述。综上，根据本发明提出的缺测数据插补方法对表1的缺测数据进行插补后的结果如表2所示。

表2缺测段数据插补运算结果

本发明的一种实施例的工作原理：首先确定太阳能资源数据总条目，从数据库中取出需要检测太阳能资源数据，依据步长Ts计算太阳能资源数据的数据条目总量为n；对n条太阳能资源数据进行合理性检查，记录不合理的数据条目后将其中不合理的数据值置为空值；对n条太阳能资源数据进行完整性检查，计算得到有效数据完整率r_sd，通过有效数据完整率r_sd判断是否需要数据插补计算，如果需要插补计算则继续判断是否满足插补计算要求；对满足插补计算要求的n条太阳能资源数据进行缺测数据插补计算，具体步骤为：

a、将存在太阳直接辐照度的第一条太阳能资源数据对应起始时刻为t₁，则最后一条，即第n条数据对应结束时刻为t_n。两条相邻数据之间的时间间隔为数据写入步长为Ts；

b、搜索第一个实测数据，对n条太阳能资源数据自前向后进行搜索，查询并记录第一条太阳能直接辐照度不为空的太阳能资源数据的时刻，即第一条太阳直接辐照度的实测数据时刻，

DHI(t_p)≠NULL t₁≤t_p≤t_n

式中，t_p为第一个实测数据时刻；

c、搜索第一个缺测数据，对记录有太阳能直接辐照度的n条数据自前向后进行搜索，查询并记录第一条太阳能直接辐照度为空的太阳能资源数据的时刻，即第一条太阳直接辐照度的缺测数据时刻，

DHI(t_f)＝NULL t₁≤t_f≤t_n

式中，t_f为第一个缺测数据时刻；

d、由步骤c获取第一个缺测数据时刻后，若向前搜索存在太阳能直接辐照度不为空的时刻，则记为t_i，若向后搜索存在太阳能直接辐照度不为空的时刻，则记为t_j，需要满足：

e、依据t_f、t_i和t_j计算得出t_f时刻的数据，并将该数据插补到t_f时刻太阳能资源数据中；

计算缺测的太阳直接辐照度数据过程为：

t_i存在，t_j不存在，根据t_i时刻的地外水平面太阳辐照度、t_i时刻的太阳直接辐照度和t_f时刻的地外水平面太阳辐照度，计算t_f时刻的太阳直接辐照度，具体计算如下：

上式中，EHI为地外水平面太阳辐照度，DHI为太阳直接辐照度，

将计算t_f时刻的太阳直接辐照度插补到t_f时刻太阳能资源数据中；

t_i不存在，t_j存在，根据t_j时刻的地外水平面太阳辐照度、t_j时刻的太阳直接辐照度和t_f时刻的地外水平面太阳辐照度，计算t_f时刻的太阳直接辐照度，具体计算如下：

上式中，EHI为地外水平面太阳辐照度，DHI为太阳直接辐照度，

将计算t_f时刻的太阳直接辐照度插补到t_f时刻太阳能资源数据中；

若t_i和t_j都存在，因为是对列数据自前向后搜索，则会出现三种情况分别为：

情况一、t_i时刻离t_f时刻更近，即t_f–t_i<t_j–t_f；此时按照t_i存在，t_j不存在的情况对t_f时刻进行数据插补；

情况二、t_j时刻离t_f时刻更近，即t_f–t_i>t_j–t_f；此时按照t_i不存在，t_j存在的情况对t_f时刻进行数据插补；

情况三、t_i时刻和t_j时刻离t_f时刻一样近，即t_f–t_i＝t_j–t_f，此时需要考虑天气现象，介于本发明所监测的太阳能资源数据，以太阳总辐照度、大气温度和大气湿度作为天气差异情况的衡量指标，根据相关工程经验并结合各指标对光伏电站的影响情况，设定其指标权重分别占比，太阳总辐照度占比为0.6、大气温度占比为0.1，大气湿度占比为0.3，分别计算t_i时刻和t_j时刻相较于t_j时刻的天气现象差异情况，具体计算如下：

式中，GHI为太阳总辐照度，T为大气温度，H为大气湿度，w代表天气现象差异情况，w数值越小表明两个时刻的天气现象更加相近；

进一步地，根据天气现象差异情况对t_f时刻进行数据插补，具体计算如下：

f、计算插补完整率r_p，并根据r_p判断是否继续进行缺测数据插补计算，如r_p<100％，则依次进行c、d、e和f步骤，如r_p＝100％，则完成本次缺测数据插补计算。

Claims (10)

1.一种用于光伏电站的太阳能资源缺测数据插补方法，其特征在于，包括以下步骤：

S110、确定数据总条目，具体为：从数据库中取出需要检测太阳能资源数据，依据设置的步长Ts计算太阳能资源数据的数据条目总量为n；

S120、对n条太阳能资源数据进行合理性检查，记录不合理的数据条目后将其中不合理的数据值置为空值，太阳能资源数据中有置空项该条数据即为缺测数据；

S130、对n条太阳能资源数据进行完整性检查，计算得到有效数据完整率r_sd，通过有效数据完整率r_sd判断是否需要数据插补计算，如果需要插补计算则继续判断是否满足插补计算要求，满足插补计算要求则进行下一步；

S140、对满足插补计算要求的n条太阳能资源数据进行缺测数据插补计算，具体步骤为：

a、定义数据时刻；

b、搜索第一个实测数据并记录第一个实测数据时刻为t_p；

c、搜索第一个缺测数据并记录第一个缺测数据时刻为t_f；

d、由第t_f向前搜索，离t_f最近的第一个实测时刻为t_i，同时由第t_f向后搜索，离t_f最近的第一个实测时刻为t_j；

e、依据t_f、t_i和t_j计算得出t_f时刻的数据，并将该数据插补到t_f时刻太阳能资源数据中；

f、计算插补完整率r_p，并根据r_p判断是否继续进行缺测数据插补计算，如r_p<100％，则依次进行c、d、e和f步骤，如r_p＝100％，则完成本次缺测数据插补计算。

2.根据权利要求1所述的用于光伏电站的太阳能资源缺测数据插补方法，其特征在于，所述步骤S110中太阳能资源数据的数据条目总量为n，其中每条太阳能资源数据包括的字段有太阳总辐照度、太阳直接辐照度、太阳散射辐照度、风速风向、大气压强、大气温度和大气湿度，结合规范可知数据写入步长Ts应该小于1小时。

3.根据权利要求1所述的用于光伏电站的太阳能资源缺测数据插补方法，其特征在于，所述步骤S120中对太阳能资源数据进行合理性检查，是根据《太阳能资源评估方法》(GB/T37256—2019)规定，从气候学界限值、内部一致性、变化范围3个方面对n条太阳能资源数据的合理性进行检查，记录不合理的数据条目，将其中不合理的数据值置为空值，太阳能资源数据中有置空项该条数据即为缺测数据。

4.根据权利要求1所述的用于光伏电站的太阳能资源缺测数据插补方法，其特征在于，所述步骤S130具体为：

S131、数据完整性检查：进一步检查n条太阳能资源数据的完整性，根据《太阳能资源评估方法》(GB/T 37256—2019)规定，太阳能资源数据的实测数据记录应包括至少1年的连续、完整数据，数据记录应至少包括小时值，小时值的数据完整率应不低于95％；即n条太阳能资源数据的有效数据完整率应不低于95％，有效数据完整率r_sd的定义如下：

式中，N₀为预期记录的数据总条目，N₁为没有太阳直接辐照度的数据条目，N₂为不合理的数据条目；

S132、依据数据完整率判断是否满足数据插补条件：假定所缺测的太阳能资源要素为太阳直接辐照度，DHI代表太阳直接辐照度；

首先判断是否满足数据插补要求，若太阳直接辐照度数据的完整率r_sd＝100％，则无需进行插补；若95％≤r_sd<100％，则需要对其进行数据插补；若r_sd<95％，则判断为实测数据异常不满足数据插补要求。

5.根据权利要求1所述的用于光伏电站的太阳能资源缺测数据插补方法，其特征在于，所述步骤S140具体为：

S141、定义数据时刻：将存在太阳直接辐照度的第一条太阳能资源数据对应起始时刻为t₁，则最后一条，即第n条数据对应结束时刻为t_n，两条相邻数据之间的时间间隔为步长Ts，n条数据里至少应包括1条实测数据；

S142、搜索第一个实测数据，对n条太阳能资源数据自前向后进行搜索，查询并记录第一条太阳能直接辐照度不为空的太阳能资源数据的时刻，即第一条太阳直接辐照度的实测数据时刻，见下式，

DHI(t_p)≠NULL t₁≤t_p≤t_n

式中，t_p为第一个实测数据时刻；

S143、搜索第一个缺测数据，对记录有太阳能直接辐照度的n条数据自前向后进行搜索，查询并记录第一条太阳能直接辐照度为空的太阳能资源数据的时刻，即第一条太阳直接辐照度的缺测数据时刻，

DHI(t_f)＝NULL t₁≤t_f≤t_n

式中，t_f为第一个缺测数据时刻；

S144、由步骤S143获取第一个缺测数据时刻t_f后，若向前搜索存在太阳能直接辐照度不为空的时刻，则记为t_i，若向后搜索存在太阳能直接辐照度不为空的时刻，则记为t_j，需要满足：

S145、依据t_f、t_i和t_j计算得出t_f时刻的太阳直接辐照度数据，并将该数据插补到t_f时刻太阳能资源数据的太阳直接辐照度数据项中；

S146、在步骤S145的数据插补成功后，将n条数据自前向后进行检索，当太阳直接辐照度数据为空值时进行记录，检索完成后所记录的空值条目总数即为m，计算所插补数据的插补完整率r_p，如r_p<100％，则依次进行S143、S144、S145和S146步骤，如r_p＝100％，则完成本次缺测数据插补计算；

插补完整率r_p的计算公式为：

式中，M₀为太阳能资源数据条目总量，这里为n条；M₁为太阳直接辐照度数据为空的条目总量，这里为m条。

6.根据权利要求5所述的用于光伏电站的太阳能资源缺测数据插补方法，其特征在于，所述步骤S145中计算缺测的太阳直接辐照度数据：

S1451、t_i存在，t_j不存在，根据t_i时刻的地外水平面太阳辐照度、t_i时刻的太阳直接辐照度和t_f时刻的地外水平面太阳辐照度，计算t_f时刻的太阳直接辐照度，具体计算如下：

上式中，EHI为地外水平面太阳辐照度，DHI为太阳直接辐照度，

将计算t_f时刻的太阳直接辐照度插补到t_f时刻太阳能资源数据对应数据项中；

S1452、t_i不存在，t_j存在，根据t_j时刻的地外水平面太阳辐照度、t_j时刻的太阳直接辐照度和t_f时刻的地外水平面太阳辐照度，计算t_f时刻的太阳直接辐照度，具体计算如下：

上式中，EHI为地外水平面太阳辐照度，DHI为太阳直接辐照度，

将计算t_f时刻的太阳直接辐照度插补到t_f时刻太阳能资源数据对应数据项中；

S1453、若t_i和t_j都存在，因为是对列数据自前向后搜索，则会出现三种情况分别为：

a、t_i时刻离t_f时刻更近，即t_f–t_i<t_j–t_f；此时按照步骤S1451中的计算公式对t_f时刻进行数据插补；

b、t_j时刻离t_f时刻更近，即t_f–t_i>t_j–t_f；此时按照步骤S1452中的计算公式对t_f时刻进行数据插补；

c、t_i时刻和t_j时刻离t_f时刻一样近，即t_f–t_i＝t_j–t_f，此时需要考虑天气现象，介于本发明所监测的太阳能资源数据，以太阳总辐照度、大气温度和大气湿度作为天气差异情况的衡量指标，根据相关工程经验并结合各指标对光伏电站的影响情况，设定其指标权重分别占比，太阳总辐照度占比为0.6、大气温度占比为0.1，大气湿度占比为0.3，分别计算t_i时刻和t_j时刻相较于t_j时刻的天气现象差异情况，具体计算如下：

式中，GHI为太阳总辐照度，T为大气温度，H为大气湿度，w代表天气现象差异情况，w数值越小表明两个时刻的天气现象更加相近；

进一步地，根据天气现象差异情况对t_f时刻进行数据插补，具体计算如下：

7.根据权利要求1所述的用于光伏电站的太阳能资源缺测数据插补方法，其特征在于，还包括有太阳能资源数据监测系统，所述太阳能资源数据监测系统包括光伏电站、数据采集装置和数据处理服务器，光伏电站内布置有数据采集装置，数据采集装置电性连接数据处理服务器。

8.根据权利要求7所述的用于光伏电站的太阳能资源缺测数据插补方法，其特征在于，所述数据采集装置包括总辐射传感器、直接辐射传感器、散射辐射传感器、风速风向传感器、气压传感器、温度传感器、湿度传感器和数据采集卡，总辐射传感器、直接辐射传感器、散射辐射传感器、风速风向传感器、气压传感器、温度传感器和湿度传感器均电性连接于数据采集卡；

数据采集卡采集光伏电站的重要太阳能资源数据，所述重要太阳能资源数据包括太阳总辐照度、太阳直接辐照度、太阳散射辐照度、风速风向、大气压强、大气温度和大气湿度。

9.根据权利要求8所述的用于光伏电站的太阳能资源缺测数据插补方法，其特征在于，所述数据处理服务器中配置有通讯接口、上位机LabVIEW程序和数据库，所述数据采集卡电性连接通讯接口，所述上位机LabVIEW程序将通讯接口接收到的重要太阳能资源数据进行处理，并将处理后的太阳能资源数据存储到数据库。

10.根据权利要求9所述的用于光伏电站的太阳能资源缺测数据插补方法，其特征在于，所述上位机LabVIEW程序包括用户登录模块、数据采集模块、数据存储模块、状态推送模块和报表打印模块；

所述用户登录模块具体为：用户登录模块基于数据库绑定用户表和角色表的动态信息权限，角色为电站管理员的用户赋予所有功能权限，角色为监测员的用户可查看LabVIEW前面板并实时监测光伏电站系统的太阳能资源数据，角色为工程师的用户可查看后台程序框图并对其功能进行修改完善；

所述数据采集模块具体为：采用“生产者——消费者”模型作为程序设计的主体结构，进一步设定数据采集卡的各项工作参数，包括物理通道，采样率、采样数和数据写入时间步长，对光伏电站所需监测的太阳能资源数据进行采集；

所述数据存储模块具体为：根据所设定的数据写入步长，将光伏电站的各项太阳能资源数据在LabVIEW前面板进行实时显示并自动写入数据库做后续分析；

所述状态推送模块具体为：通过对数据库中的历史数据进一步分析后可对当地太阳能资源进行评估，然后将评估结果实时推送给电站管理员；

所述报表打印模块：在用户需要查看资源报表时，该模块可以按逐时、逐日、逐月、逐年等条件将光伏电站的太阳能资源数据进行筛选统计并打印。

Images