CN117610780A - 一种多维多时间尺度能源间互补性评价方法 - Google Patents
一种多维多时间尺度能源间互补性评价方法 Download PDFInfo
- Publication number
- CN117610780A CN117610780A CN202311579636.1A CN202311579636A CN117610780A CN 117610780 A CN117610780 A CN 117610780A CN 202311579636 A CN202311579636 A CN 202311579636A CN 117610780 A CN117610780 A CN 117610780A
- Authority
- CN
- China
- Prior art keywords
- wind
- water
- complementarity
- light
- photovoltaic
- 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.)
- Pending
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 14
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 52
- 238000004364 calculation method Methods 0.000 claims abstract description 32
- 238000000034 method Methods 0.000 claims abstract description 24
- 238000010606 normalization Methods 0.000 claims abstract description 17
- 230000000295 complement effect Effects 0.000 claims abstract description 14
- 238000004458 analytical method Methods 0.000 claims abstract description 6
- 238000012545 processing Methods 0.000 claims abstract description 6
- 230000005611 electricity Effects 0.000 claims description 7
- 238000010248 power generation Methods 0.000 claims description 3
- 230000003313 weakening effect Effects 0.000 claims description 3
- -1 light Substances 0.000 claims description 2
- 238000011161 development Methods 0.000 abstract description 2
- 230000007547 defect Effects 0.000 description 5
- 241000345998 Calamus manan Species 0.000 description 2
- 230000002776 aggregation Effects 0.000 description 2
- 238000004220 aggregation Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 235000012950 rattan cane Nutrition 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000009194 climbing Effects 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Data Mining & Analysis (AREA)
- Tourism & Hospitality (AREA)
- Entrepreneurship & Innovation (AREA)
- General Business, Economics & Management (AREA)
- Operations Research (AREA)
- Marketing (AREA)
- Health & Medical Sciences (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Educational Administration (AREA)
- Software Systems (AREA)
- Wind Motors (AREA)
- Primary Health Care (AREA)
- Public Health (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Algebra (AREA)
- Quality & Reliability (AREA)
- Databases & Information Systems (AREA)
- General Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Water Supply & Treatment (AREA)
Abstract
本发明公开一种多维多时间尺度能源间互补性评价方法,包括:初始计算条件准备:包括新能源历史实际出力数据和水电的历史实际出力数据;对风电、水电和光伏出力数据进行Max‑Min线性归一化处理;计算风电、水电和光伏各自的装机容量占比;根据水风光装机容量的相对差值,计算折减系数;计算i时刻风‑水光,水‑风光和光‑水风三种组合的联合出力;计算联合波动量和单独波动量;计算在相应时间尺度T内风电,水电和光伏出力的互补性指数;本发明有效提高了互补性量化指标的真实性和可靠性,为水风光一体化基地开发提供了互补分析支撑。
Description
技术领域
本发明涉及能源间的互补特性评价技术领域,具体地指一种多维多时间尺度能源间互补性评价方法。
背景技术
可变可再生能源(VREs,主要包括风电、光伏和水电)发电与气象因素紧密耦合而呈现强堆积波动性,导致以新能源为主体的新型电力系统面临巨大的灵活性需求挑战。得益于水电的快速爬坡能力和水库储能,水风光多能互补运行已成为电网应对风光不确定性,保障可再生能源大规模消纳的重要手段。有效量化评价水风光多能源互补能力是开展水风光多能互补规划和调度运行的重要基础。在电力系统领域中,互补性是对风电、光伏、水电等多种VREs组合出力曲线经聚合之后的波动性相比聚合前较单一的VREs组合的波动性变化这一现象的描述。现阶段对于互补能力的评价受到很多关键因素的制约,现有的评价指标/方法
(1)仅适用于评估水光互补、水光互补和风光互补等两种能源之间的互补性;
(2)不能有效捕捉待分析时间序列之间的负相关,包括整体波动趋势和相邻时间段的波动;
(3)很大程度上忽略了装机容量对互补性的影响;
(4)只关注单一时间尺度上的互补性。
发明内容
本发明的目的在于克服上述不足,提供一种多维多时间尺度能源间互补性评价方法,以解决背景技术中提出的问题。
本发明为解决上述技术问题,所采用的技术方案是:一种多维多时间尺度能源间互补性评价方法,它包括如下步骤:
S1、初始计算条件准备:包括新能源历史实际出力数据和水电的历史实际出力数据;
S2、对风电、水电和光伏出力数据进行Max-Min线性归一化处理;
S3、计算风电、水电和光伏各自的装机容量占比;
S4、根据水风光装机容量的相对差值,计算折减系数,在装机容量大小差距过大时,对互补性进行削弱;
S5、计算i时刻风-水光,水-风光和光-水风三种组合的联合出力;
S6、计算i→j时刻的联合波动量和单独波动量;
S7、计算在相应时间尺度T内风电,水电和光伏出力的互补性指数,用C表示,C的取值范围为[0,1],越接近1代表互补性越好,反之越接近0代表互补性越差;
S8、在给定风水光装机容量之和时,得到不同组合下的风水光的互补特性,进而引入权重,为不同时间尺度的互补性指数赋予相应的重要性,实现考虑了多时间尺度下最优风水光装机配比的计算,为联合发电系统的规划提供有效辅助。
进一步地,所述步骤S1中,新能源历史实际出力数据包括风电的历史实际出力数据和光伏的历史实际出力数据;
所述步骤S2中,对风电、水电和光伏出力数据进行Max-Min线性归一化处理的计算公式如下:
式中:P'(i)为i时刻风电、水电或光伏的出力;Pmax,Pmin分别为风电、水电或光伏时间序列的最大值和最小值;P(i)为i时刻归一化后的风电、水电或光伏的出力。
进一步地,所述步骤S3中,计算风电、水电和光伏各自的装机容量占比的计算公式如下:
式中:ICW为风电的装机容量;ICH为水电的装机容量;ICS为光伏的装机容量;α为风电的装机容量占比;β为水电的装机容量占比;γ为光伏的装机容量占比。
进一步地,所述步骤S4中,折减系数的计算公式如下:
式中:a为风-水光组合时的折减系数;b为水-风光组合时的折减系数;c为光-水风组合时的折减系数。
进一步地,所述步骤S5中,计算公式如下:
式中:Pw(i),Ph(i),Ps(i)分别为i时刻风电,水电和光伏的归一化后的理论出力;Pw→h,s(i),Ph→w,s(i),Ps→w,h(i)分别为i时刻归一化后的风-水光组合,水-风光组合和光-水风组合下以装机容量为权重的联合出力。
进一步地,所述步骤S6中,计算i→j时刻的联合波动量和单独波动量,其中j>i,计算公式如下:
式中:ΔPw(i,j),ΔPh(i,j),ΔPs(i,j)为归一化后的风电,水电和光伏的单独出力波动量;ΔPh,s(i,j),ΔPw,s(i,j),ΔPw,h(i,j)为归一化后的水光,风光和风水的联合出力波动量;ΔPw→h,s(i,j),ΔPh→w,s(i,j),ΔPs→w,h(i,j)为归一化后的风-水光,水-风光和光-水风组合下的以装机容量为权重的联合出力波动量。
进一步地,所述步骤S7中,具体计算公式如下:
式中:N为分析时间尺度内的数据量。
进一步地,所述步骤S8中,具体计算公式如下:
C=pyearCyear+pmonthCmonth+pdayCday+phourChour
式中:Cyear,Cmonth,Cday,Chour分别为年、月、日和小时尺度下互补性指数;pyear,pmonth,pday,phour则为其对应的权重。
本发明的有益效果:
1、本发明方法从肯达尔系数的秩相关分析和藤状的图形结构思想出发,以波动性互补评价为基础,能够统筹兼顾局部和整体互补性,在不受振动的影响对三种以上能源间的互补性进行评价;
2、本发明通过藤状的物理结构克服了传统指标无法计算三种以及更高维度的资源间互补性的缺陷;改进后的指标同时考虑了联合的波动性和水风光单独的波动性,消除了风光波动性对互补性的影响;且在计算中以秩相关计算的思想,打破了时段的连续性,有效捕捉待分析时间序列整体波动趋势和相邻时间段的波动的影响,将所有时间尺度纳入统一的评价框架,消除了传统指标无法综合各种时间尺度的缺陷;
3、本发明以波动性互补评价为基础,能够统筹兼顾局部和整体互补性,在不受振动的影响对三种以上能源间的互补性进行评价,有效提高了互补性量化指标的真实性和可靠性,为水风光一体化基地开发提供了互补分析支撑。
附图说明
图1为本发明方法流程图;
图2为本发明方法涉及的指标物理结构图;
图3为本发明针对的四种特殊情况示意图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步的详细描述。
如图1所示,一种多维多时间尺度能源间互补性评价方法,它包括如下步骤:
S1、初始计算条件准备,包括新能源历史实际出力数据和水电的历史实际出力数据。所述的新能源包括:风电、光伏;
S2、对风电、水电和光伏出力数据进行Max-Min线性归一化处理,计算公式如下:
式中:P'(i)为i时刻风电、水电或光伏的出力;Pmax,Pmin分别为风电、水电或光伏时间序列的最大值和最小值;P(i)为i时刻归一化后的风电、水电或光伏的出力。
S3、计算风电、水电和光伏各自的装机容量占比,计算公式如下:
式中:ICW为风电的装机容量;ICH为水电的装机容量;ICS为光伏的装机容量;α为风电的装机容量占比;β为水电的装机容量占比;γ为光伏的装机容量占比。
S4、根据水风光装机容量的相对差值,计算折减系数,在装机容量大小差距过大时,对互补性进行削弱,计算公式如下:
式中:a为风-水光组合时的折减系数;b为水-风光组合时的折减系数;c为光-水风组合时的折减系数。
S5、计算i时刻风-水光,水-风光和光-水风三种组合的联合出力,计算公式如下
式中:Pw(i),Ph(i),Ps(i)分别为i时刻风电,水电和光伏的归一化后的理论出力;Pw→h,s(i),Ph→w,s(i),Ps→w,h(i)分别为i时刻归一化后的风-水光组合,水-风光组合和光-水风组合下以装机容量为权重的联合出力。
S6、计算i→j时刻的联合波动量和单独波动量,其中j>i,计算公式如下:
式中:ΔPw(i,j),ΔPh(i,j),ΔPs(i,j)为归一化后的风电,水电和光伏的单独出力波动量;ΔPh,s(i,j),ΔPw,s(i,j),ΔPw,h(i,j)为归一化后的水光,风光和风水的联合出力波动量;ΔPw→h,s(i,j),ΔPh→w,s(i,j),ΔPs→w,h(i,j)为归一化后的风-水光,水-风光和光-水风组合下的以装机容量为权重的联合出力波动量。
S7、计算在相应时间尺度T内风电,水电和光伏出力的互补性指数,用C表示,C的取值范围为[0,1],越接近1代表互补性越好,反之越接近0代表互补性越差,具体计算公式如下:
式中:N为分析时间尺度内的数据量。
S8、在给定风水光装机容量之和时,可以得到不同组合下的风水光的互补特性,进而引入权重,为不同时间尺度的互补性指数赋予相应的重要性,能够实现考虑了多时间尺度下最优风水光装机配比的计算,为联合发电系统的规划提供有效辅助,计算公式如下:
C=pyearCyear+pmonthCmonth+pdayCday+phourChour
式中:Cyear,Cmonth,Cday,Chour分别为年、月、日和小时尺度下互补性指数;pyear,pmonth,pday,phour则为其对应的权重。
本实施例通过藤的物理结构克服了传统指标无法计算三种以及更高维度的资源间互补性的缺陷,指标物理结构如图2所示。
改进后的指标同时考虑了联合的波动性和水风光单独的波动性,消除了风光波动性对互补性的影响。且在计算中以秩相关计算的思想,打破了时段的连续性,有效捕捉待分析时间序列整体波动趋势和相邻时间段的波动的影响,将所有时间尺度纳入统一的评价框架,消除了传统指标无法综合各种时间尺度的缺陷。
本实施例方法具体实施后的效果:列举四种特殊情况如图3所示。
在图3中的情况1时,即y_1以极大幅度下降,y_2以极小幅度上升的情况,互补性较差;在图3中的情况2下,后半段风光呈现负相关,存在互补性,互补性相对情况1明显更高;在图3中的情况3,即风光同频同向振动,但整体走向相反和图3中的情况4,即风光同频反向振动,整体走向相同两种情况下,从整体趋势的角度上来看情况3具有较好的互补性,而情况4互补性较差。分别应用斯皮尔曼、肯达尔和基于波动性的互补性指标三种传统指标对互补性进行计算,对比本发明提出的指标和上述分析的互补性真实情况,结果如表1所示。可以看出的是,本发明所提出的指标在各种特殊情况下均能反映出真实的互补性。
表1特殊情况下指数比较分析
上述的实施例仅为本发明的优选技术方案,而不应视为对于本发明的限制,本申请中的实施例及实施例中的特征在不冲突的情况下,可以相互任意组合。本发明的保护范围应以权利要求记载的技术方案,包括权利要求记载的技术方案中技术特征的等同替换方案为保护范围。即在此范围内的等同替换改进,也在本发明的保护范围之内。
Claims (8)
1.一种多维多时间尺度能源间互补性评价方法,其特征在于:它包括如下步骤:
S1、初始计算条件准备:包括新能源历史实际出力数据和水电的历史实际出力数据;
S2、对风电、水电和光伏出力数据进行Max-Min线性归一化处理;
S3、计算风电、水电和光伏各自的装机容量占比;
S4、根据水风光装机容量的相对差值,计算折减系数,在装机容量大小差距过大时,对互补性进行削弱;
S5、计算i时刻风-水光,水-风光和光-水风三种组合的联合出力;
S6、计算i→j时刻的联合波动量和单独波动量;
S7、计算在相应时间尺度T内风电,水电和光伏出力的互补性指数,用C表示,C的取值范围为[0,1],越接近1代表互补性越好,反之越接近0代表互补性越差;
S8、在给定风水光装机容量之和时,得到不同组合下的风水光的互补特性,进而引入权重,为不同时间尺度的互补性指数赋予相应的重要性,实现考虑了多时间尺度下最优风水光装机配比的计算,为联合发电系统的规划提供有效辅助。
2.根据权利要求1所述的一种多维多时间尺度能源间互补性评价方法,其特征在于:所述步骤S1中,新能源历史实际出力数据包括风电的历史实际出力数据和光伏的历史实际出力数据;
所述步骤S2中,对风电、水电和光伏出力数据进行Max-Min线性归一化处理的计算公式如下:
式中:P'(i)为i时刻风电、水电或光伏的出力;Pmax,Pmin分别为风电、水电或光伏时间序列的最大值和最小值;P(i)为i时刻归一化后的风电、水电或光伏的出力。
3.根据权利要求1所述的一种多维多时间尺度能源间互补性评价方法,其特征在于:所述步骤S3中,计算风电、水电和光伏各自的装机容量占比的计算公式如下:
式中:ICW为风电的装机容量;ICH为水电的装机容量;ICS为光伏的装机容量;α为风电的装机容量占比;β为水电的装机容量占比;γ为光伏的装机容量占比。
4.根据权利要求1所述的一种多维多时间尺度能源间互补性评价方法,其特征在于:所述步骤S4中,折减系数的计算公式如下:
式中:a为风-水光组合时的折减系数;b为水-风光组合时的折减系数;c为光-水风组合时的折减系数。
5.根据权利要求1所述的一种多维多时间尺度能源间互补性评价方法,其特征在于:所述步骤S5中,计算公式如下:
式中:Pw(i),Ph(i),Ps(i)分别为i时刻风电,水电和光伏的归一化后的理论出力;Pw→h,s(i),Ph→w,s(i),Ps→w,h(i)分别为i时刻归一化后的风-水光组合,水-风光组合和光-水风组合下以装机容量为权重的联合出力。
6.根据权利要求1所述的一种多维多时间尺度能源间互补性评价方法,其特征在于:所述步骤S6中,计算i→j时刻的联合波动量和单独波动量,其中j>i,计算公式如下:
式中:ΔPw(i,j),ΔPh(i,j),ΔPs(i,j)为归一化后的风电,水电和光伏的单独出力波动量;ΔPh,s(i,j),ΔPw,s(i,j),ΔPw,h(i,j)为归一化后的水光,风光和风水的联合出力波动量;ΔPw→h,s(i,j),ΔPh→w,s(i,j),ΔPs→w,h(i,j)为归一化后的风-水光,水-风光和光-水风组合下的以装机容量为权重的联合出力波动量。
7.根据权利要求1所述的一种多维多时间尺度能源间互补性评价方法,其特征在于:所述步骤S7中,具体计算公式如下:
式中:N为分析时间尺度内的数据量。
8.根据权利要求1所述的一种多维多时间尺度能源间互补性评价方法,其特征在于:所述步骤S8中,具体计算公式如下:
C=pyearCyear+pmonthCmonth+pdayCday+phourChour
式中:Cyear,Cmonth,Cday,Chour分别为年、月、日和小时尺度下互补性指数;pyear,pmonth,pday,phour则为其对应的权重。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311579636.1A CN117610780A (zh) | 2023-11-22 | 2023-11-22 | 一种多维多时间尺度能源间互补性评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311579636.1A CN117610780A (zh) | 2023-11-22 | 2023-11-22 | 一种多维多时间尺度能源间互补性评价方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117610780A true CN117610780A (zh) | 2024-02-27 |
Family
ID=89957456
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311579636.1A Pending CN117610780A (zh) | 2023-11-22 | 2023-11-22 | 一种多维多时间尺度能源间互补性评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117610780A (zh) |
-
2023
- 2023-11-22 CN CN202311579636.1A patent/CN117610780A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Eltamaly et al. | Load management as a smart grid concept for sizing and designing of hybrid renewable energy systems | |
CN109256799B (zh) | 一种基于样本熵的新能源电力系统优化调度方法 | |
CN107276122B (zh) | 适应大规模可再生能源并网的调峰资源调用决策方法 | |
Ghiassi-Farrokhfal et al. | Solar power shaping: An analytical approach | |
CN105281362A (zh) | 电力系统消纳风电能力的评估方法 | |
CN103366030A (zh) | 风电与抽水蓄能电站联合承担系统基荷运行模拟方法 | |
Tahir et al. | Optimal ESS size calculation for ramp rate control of grid-connected microgrid based on the selection of accurate representative days | |
CN110046780B (zh) | 基于分时电价的孤岛微电网需求响应经济调度方法及系统 | |
CN105262146A (zh) | 含风电的电力系统备用容量计算方法和系统 | |
Mégel et al. | Reducing the computational effort of stochastic multi-period dc optimal power flow with storage | |
CN109543993B (zh) | 分析光伏电站的方法、计算机存储介质及计算机设备 | |
CN114301089A (zh) | 一种风光联合发电系统储能容量配置优化方法 | |
CN112736899A (zh) | 一种微电网规划方案评价指标计算方法以及装置 | |
Ding et al. | A dynamic period partition method for time-of-use pricing with high-penetration renewable energy | |
CN110571862A (zh) | 一种光伏电站与电力负荷的时序匹配度分析方法及系统 | |
CN111882228A (zh) | 一种含分布式电源的配电网可靠性评估方法 | |
CN117610780A (zh) | 一种多维多时间尺度能源间互补性评价方法 | |
Hyrzyński et al. | Complementarity of wind and photovoltaic power generation in conditions similar to Polish | |
CN115600838A (zh) | 一种水电站调节能力评价方法、设备和介质 | |
Hajiah et al. | Optimal sizing of wind power systems in three high wind potential zones in kuwait for remote housing electrification | |
Born et al. | Development and demonstration of a renewable energy based demand/supply decision support tool for the building design profession | |
Lin et al. | Scenarios analysis method based on wind and photovoltaic power output correlation | |
Arshad et al. | Optimized Sizing of Stand-Alone PV Systems for Off-Grid Electrification | |
CN111126772A (zh) | 一种模拟电力现货市场平衡的方法及系统 | |
CN110619439A (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 |