CN109886567B - 一种考虑体感温度和辐射强度的短期负荷预测方法 - Google Patents

一种考虑体感温度和辐射强度的短期负荷预测方法 Download PDF

Info

Publication number
CN109886567B
CN109886567B CN201910101802.4A CN201910101802A CN109886567B CN 109886567 B CN109886567 B CN 109886567B CN 201910101802 A CN201910101802 A CN 201910101802A CN 109886567 B CN109886567 B CN 109886567B
Authority
CN
China
Prior art keywords
load
day
predicted
daily
data
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
Application number
CN201910101802.4A
Other languages
English (en)
Other versions
CN109886567A (zh
Inventor
喻乐
张珂珩
张晶
施磊
谢旭
涂孟夫
史佩然
耿琳
宁健
沈茂亚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
NARI Group Corp
North China Grid Co Ltd
Nari Technology Co Ltd
NARI Nanjing Control System Co Ltd
Original Assignee
State Grid Corp of China SGCC
NARI Group Corp
North China Grid Co Ltd
Nari Technology Co Ltd
NARI Nanjing Control System Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, NARI Group Corp, North China Grid Co Ltd, Nari Technology Co Ltd, NARI Nanjing Control System Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201910101802.4A priority Critical patent/CN109886567B/zh
Publication of CN109886567A publication Critical patent/CN109886567A/zh
Application granted granted Critical
Publication of CN109886567B publication Critical patent/CN109886567B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明公开了一种考虑体感温度和辐射强度的短期负荷预测方法,包括如下步骤:1)查询历史负荷和气象等样本数据;2)计算历史体感温度数据和日负荷水平;3)从历史样本数据集中,以待预测日的日类型信息和气象数据为依据,选择最优“模式相似日”,最终计算得到归一化负荷曲线;4)建立考虑体感温度和日照强度的神经网络预测模型,得到待预测日的负荷水平;5)通过归一化曲线和负荷水平计算待预测日的负荷数据。本发明充分考虑影响负荷的体感温度和影响分布式光伏发电的日照强度的影响,充分考虑历史负荷的自身变化规律,将负荷水平和负荷模式分开预测,减少了神经网络输入维度,降低了网络训练负担,提高了计算的效率。

Description

一种考虑体感温度和辐射强度的短期负荷预测方法
技术领域
本发明涉及一种考虑体感温度和辐射强度的短期负荷预测方法,属于电力系统负荷预测技术领域。
背景技术
目前,短期负荷预测是电网公司编制日前发电计划、实时运行控制工作开展的基础,通过准确的负荷预测,可以经济安排发电机组的启停、合理安排发电机组的计划检修,在保证电网稳定、可靠供电的同时,又能有效地降低发电成本、提高社会整体经济效益。
随着居民生活水平和生活质量的不断提高,居民生活用电负荷在电网负荷中所占的比重也越来越大,最明显的是夏季空调负荷和冬季电采暖负荷,也就是气象条件变化导致人体舒适感觉变化引起的用电负荷;另外,近几年随着国家政策影响,分布式光伏装机呈高速增长态势。这两者与电网负荷的相关性非常高,影响了短期负荷预测工作的开展。
长期以来,电网的调度工作者及科研工作者在这一课题上进行了卓有成效的研究。但是,如何综合考虑温度、湿度等气象因素和分布式光伏对短期负荷预测的影响,一直未得到有效地解决。
发明内容
目的:为了克服现有技术中存在的不足,本发明提供一种考虑体感温度和辐射强度的短期负荷预测方法。
技术方案:为解决上述技术问题,本发明采用的技术方案为:
本发明公开了一种考虑体感温度和辐射强度的短期负荷预测方法,包括以下步骤:
步骤1:查询历史的负荷、气象和日类型信息作为样本数据。样本数据包括:近5年的日96点负荷数据、日24小时气象数据和日类型信息。
步骤2:计算历史体感温度数据和日负荷水平。利用温度、湿度、风速计算相应时间的历史体感温度;对日96点负荷进行统计,分别得到每日的最高负荷和最低负荷。
步骤3:建立考虑历史体感温度和日照辐射强度的神经网络(Neural Network)预测模型,得到待预测日的负荷水平。根据近N天的历史日最高负荷值、日最低负荷值、历史体感温度、日照辐射强度建立神经网络的短期负荷水平预测模型,经过参数优化后,通过该预测模型得到待预测日的日最高负荷值、日最低负荷值。
步骤4:从历史样本数据集中,通过凝聚层次聚类算法,对近5年的日负荷曲线进行聚类分析;根据待预测日的日类型,选择特定类的实际日,对该类中实际日的气象数据和待预测日的预测气象数据做偏差分析,选择气象数据相近的作为待预测日的“模式相似日”;分别计算实际日归一化负荷系数,然后平均作为待预测日的归一化负荷系数。
步骤5:通过待预测日归一化负荷系统和日最高负荷值、日最低负荷值计算待预测日96点负荷数据。
作为优选方案,步骤1中,日24小时气象数据至少包括24小时的温度、湿度、风速、风向和日照辐射强度。日类型指某个日期为工作日还是休息日,若为工作日为,则确定为星期几;若为休息日,则确定普通周末还是特殊节假日,包括:元旦、春节、元宵、五一、端午、十一。
作为优选方案,步骤2中,历史体感温度Tg计算公式如下:
Figure BDA0001964333180000031
其中,Ta为常规气温,h为相对湿度,v为风速大小。
步骤2中,日最高负荷和日最低负荷的计算公式如下:
Lmax=max(L)
Lmin=min(L)
其中,Lmax为日最高负荷,Lmin为日最低负荷,L为日96点负荷值的数组。
作为优选方案,步骤3中,考虑体感温度和日照强度的神经网络(Neural Network)预测模型包括:
输入层:输入变量包括近N天的日24小时体感温度、辐射强度、近N天的日最高负荷值、日最低负荷值;待预测日的日24小时体感温度、辐射强度;
输出层:日最高负荷值、日最低负荷值。
作为优选方案,步骤3中,神经网络模型,求解器使用随机梯度下降(StochasticGradient Descent,SGD)算法。其激活函数采用sigmoid函数,即:
g(x)=1/(1+e-x)
式中,x是预测输入向量,g为激活函数输出值。
作为优选方案,步骤4中,凝聚层次聚类时,通过计算每个个体样本之间的欧氏距离来达到最终的聚类归并效果。
计算任意两个历史负荷之间的欧氏距离方法如下:设n维样本空间S中任意两个数据序列X、Y分别为X={x1,x2,......,xn-1,xn},Y={y1,y2,......,yn-1,yn},其欧氏距离d(X,Y)为:
Figure BDA0001964333180000041
设定负荷曲线间距离阈值为d。在样本集S中有N条样本负荷数据序列S={L1,L2,……,Ln-1,LN},任意取出负荷数据Lk(k=1,2,...,N),并计算该负荷数据与其他负荷数据序列的欧氏距离。该算法开始时,每一个负荷数据均为一个“簇”,通过上述距离公式,对相近的负荷数据进行“簇”的合并,为同一个“簇”;反复上述合并过程,直到达到预定的分类目标即可停止。分类完成后,根据负荷的自身变化特性,分类包括:工作日、周末、三天的节假日(元旦、五一、端午)、七天的节假日(十一和春节)。
作为优选方案,步骤4中,根据待预测日的日期,选择特定的分类。该特定分类集合中实际日气象数据和待预测日的气象数据偏差分析计算公式如下:
Ew=∑∝i*|Wr,i-Wf,i|
EW为集合类中实际日气象和待预测日预测气象的偏差;∝i为第i项气象因素的系数;Wr,i为实际日第i项气象因素的实际值,Wf,i为待预测日第i项气象因素的预测值。对选择范围内所有实际日的气象偏差EW按升序排序,选择偏差最小的M(M一般选择5)天作为“模式相似日”,参与归一化负荷系数计算。
步骤4中,“模式相似日”的归一化负荷系数计算的公式如下:
Ln,i=(Li-Lmin)/(Lmax-Lmin),i=1,2,..,96
Ln,i为某实际日第i点的归一化负荷系数;Li某实际日第i点的负荷值,Lmax和Lmin分别为该日最高负荷和最低负荷。
步骤4中,计算预测日的归一化负荷系数,依次对“模式相似日”中实际日第i(i=1,2,..,96)点进行如下计算:
Figure BDA0001964333180000051
Ln,i,j为“模式相似日”中第j(j=1,2,…,M)日的第i点的归一化负荷系数,Lfn,i为预测日的第i点的归一化负荷系数。
步骤5中,待预测日96点负荷值计算公式,具体如下:
Lfi=Lmin+Lfn,i*(Lmax-Lmin),i=1,2,..,96
Lfi待预测日第i点的负荷值;Lfn,i为预测日的第i点的归一化负荷系数,Lmax和Lmin分别为预测日最高负荷值和最低负荷值。
有益效果:本发明提供的一种考虑体感温度和辐射强度的短期负荷预测方法,采用体感温度和日照辐射强度作为输入因素,综合了考虑湿度、温度和分布式光伏发电对电网负荷水平的影响;同时基于层次聚类算法,保证了负荷曲线变化趋势的合理性。具有以下优点:
1、通过选用综合考虑温度、湿度和辐射的体感温度指标,综合考虑气象因素对负荷的影响;通过单独考虑辐射实现对分布式光伏的考虑。
2、通过聚类分析计算和气象偏差过滤,合理考虑了负荷自身的变化趋势。
3、通过对日最高负荷和日最低负荷分别建模,降低了神经网络维度,提高了模型训练效率。
附图说明
图1为本发明的方法流程示意图;
图2为本发明神经网络模型结构图。
具体实施方式
下面结合附图对本发明作更进一步的说明。
如图1所示,本发明的一个实施案例,包含了采用本发明方法,在进行短期负荷预测处理过程中,它的特征、目的和优点可以从实施例的步骤中看出。
一种考虑体感温度和辐射强度的短期负荷预测方法,具体包括如下步骤:
第一步:查询历史日期的日负荷和气象和日类型信息作为样本数据。样本数据主要包括:近5年的日96点负荷数据、日24小时气象数据和日类型信息。
第二步:计算历史日期的24小时历史体感温度数据、日最高负荷和日最低负荷。历史体感温度Tg计算公式如下:
Figure BDA0001964333180000061
其中,Ta为温度,h为湿度,v为风速。
第三步:根据待预测日期,选择合理的样本集;分别建立最大负荷预测训练模型和最小负荷预测训练模型。其输入量包括:输入变量包括近N天的历史体感温度、辐射强度、近N天的日最高负荷值、日最低负荷值;待预测日的历史体感温度、辐射强度;其输出为:日最高负荷、日最低负荷值。如图2所示。经过模型训练后,分别预测得到待预测日的日最高负荷和日最低负荷。
步骤4:从历史样本数据集中,通过凝聚层次聚类算法,对近5年的日负荷曲线进行聚类分析;根据待预测日的日类型,选择特定类的实际日,对该类中实际日的气象数据和待预测日的预测气象数据做偏差分析,选择气象数据相近的作为待预测日的“模式相似日”;分别计算实际日归一化负荷系数,然后平均作为待预测日的归一化负荷系数。
步骤5:通过待预测日归一化负荷系统和日最高负荷值、日最低负荷值计算待预测日96点负荷数据。
本发明实际的应用效果:
本发明技术方案在某省级电网短期系统负荷预测中得到应用,应用效果符合预期。实际应用表明,本发明能够综合考虑体感温度和日照辐射强度对电网负荷的影响,在减少计算量的基础上,有效降低预测误差,提高预测精度。
以上显示和描述了本发明的基本原理、主要特征及优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (8)

1.一种考虑体感温度和辐射强度的短期负荷预测方法,其特征在于:包括如下步骤:
步骤1:查询历史的负荷、气象和日类型信息作为样本数据;
步骤2:计算历史体感温度数据和日负荷水平;利用温度、湿度、风速计算相应时间的历史体感温度;对日96点负荷进行统计,分别得到每日的最高负荷和最低负荷;
步骤3:建立考虑历史体感温度和日照辐射强度的神经网络预测模型,得到待预测日的负荷水平;根据近N天的历史日最高负荷值、日最低负荷值、历史体感温度、日照辐射强度建立神经网络的短期负荷水平预测模型,经过参数优化后,通过该预测模型得到待预测日的日最高负荷值、日最低负荷值;
步骤4:从历史样本数据集中,通过凝聚层次聚类算法,对近N年的日负荷曲线进行聚类分析;根据待预测日的日类型,选择特定类的实际日,对特定类中实际日的气象数据和待预测日的预测气象数据做偏差分析,选择气象数据相近的作为待预测日的“模式相似日”;分别计算实际日归一化负荷系数,然后平均作为待预测日的归一化负荷系数;
步骤5:通过待预测日归一化负荷系统和日最高负荷值、日最低负荷值计算待预测日96点负荷数据;
所述步骤2中,历史体感温度Tg计算公式如下:
Figure FDA0003716917770000011
其中,Ta为常规气温,h为相对湿度,v为风速大小;
日最高负荷和日最低负荷的计算公式如下:
Lmax=max(L)
Lmin=min(L)
其中,Lmax为日最高负荷,Lmin为日最低负荷,L为日96点负荷值的数组。
2.根据权利要求1所述的一种考虑体感温度和辐射强度的短期负荷预测方法,其特征在于:步骤1中,日24小时气象数据至少包括24小时的温度、湿度、风速、风向和日照辐射强度;日类型指某个日期为工作日还是休息日,若为工作日,则确定为星期几;若为休息日,则确定普通周末还是特殊节假日,包括:元旦、春节、元宵、五一、端午、十一。
3.根据权利要求1所述的一种考虑体感温度和辐射强度的短期负荷预测方法,其特征在于:所述步骤3中,考虑体感温度和日照强度的神经网络预测模型包括:
输入层:输入变量包括近N天的日24小时体感温度、辐射强度、近N天的日最高负荷值、日最低负荷值;待预测日的日24小时体感温度、辐射强度;
输出层:日最高负荷值、日最低负荷值;
所述神经网络预测模型的求解器使用随机梯度下降算法,其激活函数采用sigmoid函数,即:
g(x)=1/(1+e-x)
式中,x是预测输入向量,g为激活函数输出值。
4.根据权利要求1所述的一种考虑体感温度和辐射强度的短期负荷预测方法,其特征在于:所述步骤4中,凝聚层次聚类时,通过计算每个个体样本之间的欧氏距离来达到最终的聚类归并效果;
计算任意两个历史负荷之间的欧氏距离方法如下:设n维样本空间S中任意两个数据序列X、Y分别为X={x1,x2,……,xn-1,xn},Y={y1,y2,……,yn-1,yn},其欧氏距离d(X,Y)为:
Figure FDA0003716917770000031
设定负荷曲线间距离阈值为d;在样本集S中有N条样本负荷数据序列S={L1,L2,……,Ln-1,LN},任意取出负荷数据Lk,k=1,2,…,N,并计算该负荷数据与其他负荷数据序列的欧氏距离;分类算法开始时,每一个负荷数据均为一个“簇”,通过欧氏距离公式,对相近的负荷数据进行“簇”的合并,为同一个“簇”;反复上述合并的过程,直到达到预定的分类目标即可停止;分类完成后,根据负荷的自身变化特性,分类包括:工作日、周末、三天的节假日、七天的节假日。
5.根据权利要求1所述的一种考虑体感温度和辐射强度的短期负荷预测方法,其特征在于:所述步骤4中,根据待预测日的日期,选择特定的分类;特定分类集合中实际日气象数据和待预测日的气象数据偏差分析计算公式如下:
EW=∑∝i*|Wr,i-Wf,i|
EW为集合类中实际日气象和待预测日预测气象的偏差;∝i为第i项气象因素的系数;Wr,i为实际日第i项气象因素的实际值,Wf,i为待预测日第i项气象因素的预测值;对选择范围内所有实际日的气象偏差EW按升序排序,选择偏差最小的M天作为“模式相似日”,参与归一化负荷系数计算。
6.根据权利要求5所述的一种考虑体感温度和辐射强度的短期负荷预测方法,其特征在于:所述M设置为5。
7.根据权利要求1所述的一种考虑体感温度和辐射强度的短期负荷预测方法,其特征在于:所述步骤4中,“模式相似日”的归一化负荷系数计算的公式如下:
Ln,i=(Li-Lmin)/(Lmax-Lmin), i=1,2,..,96
Ln,i为某实际日第i点的归一化负荷系数;Li某实际日第i点的负荷值,Lmax和Lmin分别为该日最高负荷和最低负荷;
计算预测日的归一化负荷系数,依次对“模式相似日”中实际日第i点进行如下计算:
Figure FDA0003716917770000041
Ln,i,j为“模式相似日”中第j,j=1,2,…,M日的第i点的归一化负荷系数,Lfn,i为预测日的第i点的归一化负荷系数。
8.根据权利要求1所述的一种考虑体感温度和辐射强度的短期负荷预测方法,其特征在于:所述步骤5中,待预测日96点负荷值计算公式,具体如下:
Lfi=Lmin+Lfn,i*(Lmax-Lmin),i=1,2,..,96
Lfi为待预测日第i点的负荷值;Lfn,i为预测日的第i点的归一化负荷系数,Lmax和Lmin分别为预测日最高负荷值和最低负荷值。
CN201910101802.4A 2019-01-31 2019-01-31 一种考虑体感温度和辐射强度的短期负荷预测方法 Active CN109886567B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910101802.4A CN109886567B (zh) 2019-01-31 2019-01-31 一种考虑体感温度和辐射强度的短期负荷预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910101802.4A CN109886567B (zh) 2019-01-31 2019-01-31 一种考虑体感温度和辐射强度的短期负荷预测方法

Publications (2)

Publication Number Publication Date
CN109886567A CN109886567A (zh) 2019-06-14
CN109886567B true CN109886567B (zh) 2022-09-30

Family

ID=66927792

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910101802.4A Active CN109886567B (zh) 2019-01-31 2019-01-31 一种考虑体感温度和辐射强度的短期负荷预测方法

Country Status (1)

Country Link
CN (1) CN109886567B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110782094B (zh) * 2019-10-28 2022-06-21 国网江苏省电力有限公司苏州供电分公司 一种居民用户细粒度用电行为的分项负荷预测方法
CN111665879A (zh) * 2020-05-26 2020-09-15 深圳市真元天成科技有限公司 室内体感温度控制装置及控制系统、智能床垫
CN111695733B (zh) * 2020-06-11 2022-09-09 国网天津市电力公司 一种基于改进模糊聚类算法的“煤改电”负荷预测方法
CN113011045B (zh) * 2021-04-13 2022-09-30 河海大学 夏季用电负荷中空调功率获取方法及计算机可读介质
CN113536695B (zh) * 2021-08-06 2022-12-23 广东电网有限责任公司 一种基于时段划分的日负荷预测方法及相关装置
CN114742263B (zh) * 2022-03-02 2024-03-01 北京百度网讯科技有限公司 负荷预测方法、装置、电子设备及存储介质
CN116258280B (zh) * 2023-05-12 2023-08-11 国网湖北省电力有限公司经济技术研究院 一种基于时间序列聚类的短期负荷预测方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105069525A (zh) * 2015-07-30 2015-11-18 广西大学 全天候96点日负荷曲线预测及优化修正系统
CN106779129A (zh) * 2015-11-19 2017-05-31 华北电力大学(保定) 一种考虑气象因素的短期电力负荷预测方法
CN107423836A (zh) * 2017-04-07 2017-12-01 山东大学 基于体感温度的短期负荷预测方法
CN107578122A (zh) * 2017-08-21 2018-01-12 国网山东省电力公司济南供电公司 一种基于体感温度和日期类型的电网负荷预测方法与系统
CN108229754A (zh) * 2018-01-31 2018-06-29 杭州电子科技大学 基于相似日分段和lm-bp网络的短期负荷预测方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105069525A (zh) * 2015-07-30 2015-11-18 广西大学 全天候96点日负荷曲线预测及优化修正系统
CN106779129A (zh) * 2015-11-19 2017-05-31 华北电力大学(保定) 一种考虑气象因素的短期电力负荷预测方法
CN107423836A (zh) * 2017-04-07 2017-12-01 山东大学 基于体感温度的短期负荷预测方法
CN107578122A (zh) * 2017-08-21 2018-01-12 国网山东省电力公司济南供电公司 一种基于体感温度和日期类型的电网负荷预测方法与系统
CN108229754A (zh) * 2018-01-31 2018-06-29 杭州电子科技大学 基于相似日分段和lm-bp网络的短期负荷预测方法

Also Published As

Publication number Publication date
CN109886567A (zh) 2019-06-14

Similar Documents

Publication Publication Date Title
CN109886567B (zh) 一种考虑体感温度和辐射强度的短期负荷预测方法
CN109858673A (zh) 一种光伏发电系统功率预测方法
CN109934395B (zh) 一种基于分时分区气象数据的多小水电地区短期电力负荷预测方法
CN110796307B (zh) 一种综合能源系统分布式负荷预测方法及系统
CN107423836B (zh) 基于体感温度的短期负荷预测方法
CN112039056A (zh) 一种新能源两阶段优化调度方法
CN112510703B (zh) 一种考虑煤耗曲线修正的多能源接入电网优化调度方法
CN111092451B (zh) 一种基于配电网网格的光伏消纳预警方法
CN105701570A (zh) 一种基于全过程技术提升的短期电力需求分析方法
CN111612244B (zh) 基于qra-lstm的日前光伏功率非参数概率预测方法
Qiuyu et al. Short-term load forecasting based on load decomposition and numerical weather forecast
CN116187601B (zh) 一种基于负荷预测的综合能源系统运行优化方法
CN111525552B (zh) 一种基于特征信息的三阶段短期风电场群功率预测方法
CN116436008A (zh) 一种用于光储充电站的电力调度方法及终端
CN115689146A (zh) 一种多能源智能楼宇负荷调度方法及系统
CN108346009A (zh) 一种基于用户模型自学习的电力生产配置方法和装置
Sun et al. The short-term load forecasting using an artificial neural network approach with periodic and nonperiodic factors: a case study of Tai'an, Shandong Province, China
CN110147908A (zh) 一种基于三维最优相似度与改进布谷鸟算法的风电功率预测方法
CN110135617A (zh) 一种气象因素影响下基于负荷特性的改进短期负荷的预测方法
Liu et al. Short-term load prediction method for power distributing method based on back-propagation neural network
CN116454928A (zh) 一种考虑多时间尺度的多类型储能协同调度方法
CN115912419A (zh) 一种基于负荷聚合商的储热-储电协同调度方法
CN112491049B (zh) 一种考虑在线煤耗曲线的多能源接入电网优化调度方法
CN115169839A (zh) 一种基于数据-物理-知识联合驱动的供暖负荷调度方法
CN109149644B (zh) 一种基于大数据分析的光储一体在线策略匹配和协同优化方法

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