CN112308329A - 兼容中长期实物合同的日前市场出清模型构建方法 - Google Patents

兼容中长期实物合同的日前市场出清模型构建方法 Download PDF

Info

Publication number
CN112308329A
CN112308329A CN202011237229.9A CN202011237229A CN112308329A CN 112308329 A CN112308329 A CN 112308329A CN 202011237229 A CN202011237229 A CN 202011237229A CN 112308329 A CN112308329 A CN 112308329A
Authority
CN
China
Prior art keywords
unit
period
power
day
market
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.)
Granted
Application number
CN202011237229.9A
Other languages
English (en)
Other versions
CN112308329B (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.)
North China Electric Power University
State Grid Fujian Electric Power Co Ltd
Trading Center of State Grid Fujian Electric Power Co Ltd
Original Assignee
North China Electric Power University
State Grid Fujian Electric Power Co Ltd
Trading Center of State Grid Fujian Electric Power 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 North China Electric Power University, State Grid Fujian Electric Power Co Ltd, Trading Center of State Grid Fujian Electric Power Co Ltd filed Critical North China Electric Power University
Priority to CN202011237229.9A priority Critical patent/CN112308329B/zh
Publication of CN112308329A publication Critical patent/CN112308329A/zh
Application granted granted Critical
Publication of CN112308329B publication Critical patent/CN112308329B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • 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
    • G06Q30/00—Commerce
    • G06Q30/02—Marketing; Price estimation or determination; Fundraising
    • G06Q30/0201—Market modelling; Market analysis; Collecting market data
    • G06Q30/0206—Price or cost determination based on market factors
    • 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
    • G06Q30/00—Commerce
    • G06Q30/02—Marketing; Price estimation or determination; Fundraising
    • G06Q30/0283—Price estimation or determination
    • 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
    • G06Q30/00—Commerce
    • G06Q30/06—Buying, selling or leasing transactions
    • G06Q30/0601—Electronic shopping [e-shopping]
    • G06Q30/0611—Request for offers or quotes
    • 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

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Strategic Management (AREA)
  • Development Economics (AREA)
  • Economics (AREA)
  • Finance (AREA)
  • Accounting & Taxation (AREA)
  • Marketing (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Human Resources & Organizations (AREA)
  • Game Theory and Decision Science (AREA)
  • Health & Medical Sciences (AREA)
  • Tourism & Hospitality (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Water Supply & Treatment (AREA)
  • Public Health (AREA)
  • Data Mining & Analysis (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种兼容中长期实物合同的日前市场出清模型构建方法。本发明构建模型符合我国的电力交易与调管理体制，能够满足不同市场运行场景的需要；该模型有效提高日前市场出清结果的可执行性，减少实时平衡机制的调整成本，提高日前市场和实时平衡两阶段的整体社会福利。

Description

兼容中长期实物合同的日前市场出清模型构建方法
技术领域
本发明涉及一种兼容中长期实物合同的日前市场出清模型构建方法。
背景技术
全球的电力市场化改革兴起于20世纪90年代，经过三十多年的发展，欧洲、北美、澳大利亚、南美等许多国家或地区已经建立了相对完善的电力批发市场，尽管各个国家或地区的电力市场设计存在明显差异，但其共同特征是建立了分时定价的电力现货市场。2015年3月，我国开启了新一轮的电力体制改革，电改配套文件《关于推进电力市场建设的实施意见》提出在具备条件的地区开展电力现货市场试点，并将电力市场模式分为分散式和集中式两种。
分散式和集中式市场模式最大的区别在于，前者允许市场主体签订中长期的实物合同，偏离中长期实物合同的发用电需求自愿参与日前、日内等现货交易；而后者则要求市场主体必须全电量参与现货交易，中长期交易只允许签订差价合同、电力期货期权等金融性的合同。两种市场模式并没有绝对的优劣之分，需要根据各地区的电源结构、电网特点、市场结构、管理体制等因素进行选择。国内首批八个现货试点省份中，蒙西和福建选择采用分散式市场模式外，广东、浙江等其他省份选择采用集中式市场模式，各试点地区的现货市场目前仍处于模拟运行和规则完善阶段，尚未正式投入运行。
现货市场出清模型研究方面，目前国内外研究文献主要是基于传统的集中式市场模式现货出清模型，考虑风电、储能、电动汽车、需求响应、虚拟电厂等新兴市场主体的技术特点，对市场出清模型进行相应改进。关于分散式市场模式现货出清模型的研究较少，主要欧洲分散式市场的运行机制为主，现有关于分散式市场模式的研究文献没有考虑国外中长期实物合同交割方式与我国电力交易和调度管理体制的适应性及其对现货市场出清模型设计的影响。本申请在总结分析欧洲分散式市场模式特点和日前市场出清模型现状的基础上，根据我国电力交易和调度管理运行体制与欧洲的差异，设计了一种日前现货交易与中长期实物合同交割、电网阻塞管理联合优化的市场出清模型，并通过算例的对比分析，验证了本申请模型及机制的有效性。
发明内容
本发明的目的在于提供一种兼容中长期实物合同的日前市场出清模型构建方法，考虑机组爬坡和最大最小技术出力等机组运行约束，当机组的中长期实物合同不满足上述运行约束时，将按日前市场出清价被动购入或售出相应电力电量以满足约束；而对于风电、光伏等间歇性电源占比较高的电力系统，建立滚动交易的日内市场以促进可再生能源消纳，使得机组有机会主动参与日内交易以满足运行约束。
为实现上述目的，本发明的技术方案是：一种兼容中长期实物合同的日前市场出清模型构建方法，包括如下步骤：
步骤S1、构建兼容中长期实物合同的日前市场出清模型，该模型以社会福利最大化为优化目标，可对日前市场的电能与旋转备用进行联合优化出清，目标函数的具体数学表达式如下：
Figure BDA0002768201870000021
式中：△t为时段间隔；T为时段总数；NG、ND分别为机组和负荷的总数；
Figure BDA0002768201870000022
分别为t时段第i个机组和第j个负荷买电报价的数目；
Figure BDA0002768201870000023
为t时段第i个机组第k个买电报价的价格；
Figure BDA0002768201870000024
为t时段第i个机组第k个买电报价的中标电力；
Figure BDA0002768201870000025
为t时段第j个负荷第w个买电报价的价格；
Figure BDA0002768201870000026
为t时段第j个负荷第w个买电报价的中标电力；
Figure BDA0002768201870000027
分别为为t时段第i个机组和第j个负荷卖电报价的数目；
Figure BDA0002768201870000028
为t时段第i个机组第h个卖电报价的价格；
Figure BDA0002768201870000029
为t时段第i个机组第h个卖电报价的中标电力；
Figure BDA00027682018700000210
为t时段第j个负荷第v个卖电报价的价格；
Figure BDA00027682018700000211
为t时段第j个负荷第v个卖电报价的中标电力；
Figure BDA00027682018700000212
为t时段第i个机组中标的旋转备用容量；
Figure BDA00027682018700000213
为t时段第i个机组上旋转备用的报价；yi,t为启机变量，StartUpi为机组i的启机成本；
Figure BDA00027682018700000214
为第m个需求侧灵活报价在t时段的调用状态0-1变量，
Figure BDA00027682018700000215
为可中断负荷的报价；λ为供需平衡松弛变量的惩罚因子，δt为t时段供需不平衡电量；
步骤S2、设置约束条件：
1)机组和负荷中标电力的上下限约束
Figure BDA00027682018700000216
式中：
Figure BDA00027682018700000217
为t时段第i个机组第k个买电报价的容量上限；
Figure BDA00027682018700000218
为t时段第i个机组第h个卖电报价的容量上限；
Figure BDA00027682018700000219
为t时段第j个负荷第w个买电报价的容量上限；
Figure BDA00027682018700000220
为t时段第j个负荷第v个卖电报价的容量上限；
2)机组上旋转备用中标容量的上下限约束
Figure BDA00027682018700000221
式中，
Figure BDA0002768201870000031
为t时段机组i可提供的旋转备用容量上限；
3)机组出力上下限约束
Figure BDA0002768201870000032
式中：
Figure BDA0002768201870000033
为t时段机组i的实物合同位置；
Figure BDA0002768201870000034
为机组已经t时段机组i通过中长期合同已售出的旋备容量；ui,t为t时段第i个机组的启停机状态0-1变量，对于自主确定启停机的机组，ui,t设为常量1或0，不参与状态变量优化；
Figure BDA0002768201870000035
为机组i的出力上限；
Figure BDA0002768201870000036
为机组i的出力下限；
4)机组出力爬坡约束
Figure BDA0002768201870000037
式中，RUi为机组i在△t时段内的上爬坡速率；RDi为机组i在△t时段内的下爬坡速率；yi,t、zi,t分别为机组启停控制0-1变量，yi,t＝1、zi,t＝1分别表示机组i在第t时段执行启动操作和执行停机操作，否则为0；Pstart,i、Pshut,i分别为启动功率速度限制和停机功率速度限制；
5)启停机逻辑约束
Figure BDA0002768201870000038
当机组运行状态ui,t发生变化时，式(6)可以确保yi,t、zi,t取正确的值；
6)机组最小启停机时间约束
Figure BDA0002768201870000039
式中：e为t时段启停机约束校验的起始时段；Gi为机组i在初始时刻仍需连续运行的时段数；Li为机组i在初始时刻仍需连续停机的时段数；
Figure BDA00027682018700000310
为机组i的最小运行时间；
Figure BDA00027682018700000311
为机组i的最小停机时间；
7)需求侧灵活报价运行约束
Figure BDA0002768201870000041
Tm为第m个需求侧灵活报价的可调节时长；
8)系统电力平衡约束
Figure BDA0002768201870000042
9)系统旋转备用需求约束
Figure BDA0002768201870000043
式中：
Figure BDA0002768201870000044
为t时段系统总负荷需求预测；
Figure BDA0002768201870000045
为系统在t时段的旋转备用总需求；旋转备用需求可以由市场运营机构以中长期合同方式购买，不足部分可以在日前市场中购买，RSt为系统在t时段需要从日前市场中购买的旋转备用需求；
10)网络传输约束
Figure BDA0002768201870000046
式中：
Figure BDA0002768201870000047
为第l条线路的潮流传输极限；模型中潮流的求解采用直流潮流算法，Se,l为第e个节点注入功率对第l条线路潮流的灵敏度；Pe,net为第e个节点的净注入功率；
Figure BDA0002768201870000048
相较于现有技术，本发明具有以下有益效果：本发明构建模型符合我国的电力交易与调管理体制，能够满足不同市场运行场景的需要；该模型有效提高日前市场出清结果的可执行性，减少实时平衡机制的调整成本，提高日前市场和实时平衡两阶段的整体社会福利。
附图说明
图1为IEEE 30节点系统图。
图2为各场景下的日前市场SCUC优化结果。
图3为各场景下的发电侧各机组中标结果。
图4为各场景下的各节点边际出清价格。
图5为双节点示例系统。
具体实施方式
下面结合附图，对本发明的技术方案进行具体说明。
本发明提供了一种兼容中长期实物合同的日前市场出清模型构建方法，包括如下步骤：
步骤S1、构建兼容中长期实物合同的日前市场出清模型，该模型以社会福利最大化为优化目标，可对日前市场的电能与旋转备用进行联合优化出清，目标函数的具体数学表达式如下：
Figure BDA0002768201870000051
式中：△t为时段间隔；T为时段总数；NG、ND分别为机组和负荷的总数；
Figure BDA0002768201870000052
分别为t时段第i个机组和第j个负荷买电报价的数目；
Figure BDA0002768201870000053
为t时段第i个机组第k个买电报价的价格；
Figure BDA0002768201870000054
为t时段第i个机组第k个买电报价的中标电力；
Figure BDA0002768201870000055
为t时段第j个负荷第w个买电报价的价格；
Figure BDA0002768201870000056
为t时段第j个负荷第w个买电报价的中标电力；
Figure BDA0002768201870000057
分别为为t时段第i个机组和第j个负荷卖电报价的数目；
Figure BDA0002768201870000058
为t时段第i个机组第h个卖电报价的价格；
Figure BDA0002768201870000059
为t时段第i个机组第h个卖电报价的中标电力；
Figure BDA00027682018700000510
为t时段第j个负荷第v个卖电报价的价格；
Figure BDA00027682018700000511
为t时段第j个负荷第v个卖电报价的中标电力；
Figure BDA00027682018700000512
为t时段第i个机组中标的旋转备用容量；
Figure BDA00027682018700000513
为t时段第i个机组上旋转备用的报价；yi,t为启机变量，StartUpi为机组i的启机成本；
Figure BDA00027682018700000514
为第m个需求侧灵活报价在t时段的调用状态0-1变量，
Figure BDA00027682018700000515
为可中断负荷的报价；λ为供需平衡松弛变量的惩罚因子，δt为t时段供需不平衡电量；
步骤S2、设置约束条件：
1)机组和负荷中标电力的上下限约束
Figure BDA00027682018700000516
式中：
Figure BDA00027682018700000517
为t时段第i个机组第k个买电报价的容量上限；
Figure BDA00027682018700000518
为t时段第i个机组第h个卖电报价的容量上限；
Figure BDA00027682018700000519
为t时段第j个负荷第w个买电报价的容量上限；
Figure BDA00027682018700000520
为t时段第j个负荷第v个卖电报价的容量上限；
2)机组上旋转备用中标容量的上下限约束
Figure BDA00027682018700000521
式中，
Figure BDA0002768201870000061
为t时段机组i可提供的旋转备用容量上限；
3)机组出力上下限约束
Figure BDA0002768201870000062
式中：
Figure BDA0002768201870000063
为t时段机组i的实物合同位置；
Figure BDA0002768201870000064
为机组已经t时段机组i通过中长期合同已售出的旋备容量；ui,t为t时段第i个机组的启停机状态0-1变量，对于自主确定启停机的机组，ui,t设为常量1或0，不参与状态变量优化；
Figure BDA0002768201870000065
为机组i的出力上限；
Figure BDA0002768201870000066
为机组i的出力下限；
4)机组出力爬坡约束
Figure BDA0002768201870000067
式中，RUi为机组i在△t时段内的上爬坡速率；RDi为机组i在△t时段内的下爬坡速率；yi,t、zi,t分别为机组启停控制0-1变量，yi,t＝1、zi,t＝1分别表示机组i在第t时段执行启动操作和执行停机操作，否则为0；Pstart,i、Pshut,i分别为启动功率速度限制和停机功率速度限制；
5)启停机逻辑约束
Figure BDA0002768201870000068
当机组运行状态ui,t发生变化时，式(6)可以确保yi,t、zi,t取正确的值；
6)机组最小启停机时间约束
Figure BDA0002768201870000069
式中：e为t时段启停机约束校验的起始时段；Gi为机组i在初始时刻仍需连续运行的时段数；Li为机组i在初始时刻仍需连续停机的时段数；
Figure BDA00027682018700000610
为机组i的最小运行时间；
Figure BDA00027682018700000611
为机组i的最小停机时间；
8)需求侧灵活报价运行约束
Figure BDA0002768201870000071
Tm为第m个需求侧灵活报价的可调节时长；
8)系统电力平衡约束
Figure BDA0002768201870000072
9)系统旋转备用需求约束
Figure BDA0002768201870000073
式中：
Figure BDA0002768201870000074
为t时段系统总负荷需求预测；
Figure BDA0002768201870000075
为系统在t时段的旋转备用总需求；旋转备用需求可以由市场运营机构以中长期合同方式购买，不足部分可以在日前市场中购买，RSt为系统在t时段需要从日前市场中购买的旋转备用需求；
10)网络传输约束
Figure BDA0002768201870000076
式中：
Figure BDA0002768201870000077
为第l条线路的潮流传输极限；模型中潮流的求解采用直流潮流算法，Se,l为第e个节点注入功率对第l条线路潮流的灵敏度；Pe,net为第e个节点的净注入功率；
Figure BDA0002768201870000078
以下为本发明的具体实现过程。
1、兼容中长期实物合同的日前市场出清模型
1.1、报价及出清机制
本发明的日前市场出清模型考虑机组爬坡和最大最小技术出力等机组运行约束，当机组的中长期实物合同不满足上述运行约束时，将按日前市场出清价被动购入或售出相应电力电量以满足约束。对于风电、光伏等间歇性电源占比较高的电力系统，通常需要建立滚动交易的日内市场以促进可再生能源消纳，机组有机会主动参与日内交易以满足运行约束，此时在日前市场出清模型中可以不再考虑机组爬坡和最小技术出力的约束。
(1)机组自主决定启停机计划，日前市场出清优化仅涉及增量开机、不涉及停机
日前市场开市前，签有中长期实物合同的机组可自主选择在某些时段停机运行，不参与相应时段的出清优化，对应实物合同电量按日前市场出清价格购入；否则，市场出清系统默认签有实物合同的机组在各时段开机运行，参与日前市场出清优化。没有实物合同且处于停机状态的机组可申报启机报价及分段卖电报价，由市场出清模型决定其是否启机运行。
启机机组的实物合同位置在某些时段低于其最小出力时，市场出清系统自动将其实物合同位置调整为最小出力；相邻时段实物合同位置不满足机组爬坡约束时，市场出清系统自动将下一时段实物合同位置调整为满足最大爬坡约束；实物合同位置调整引起的偏差电量，按日前市场出清价格进行结算。
(2)实物合同位置作为日前市场出清的边界条件，开机机组需要同时申报反向交易报价
对于开机机组，其各时段实物合同位置由系统设为市场出清计算的边界条件，剩余发电容量可申报分段卖电报价；开机机组需要同时申报各时段反向交易报价(即买电报价)，报价容量不超过其相应时段实物合同位置与其最小技术出力之间的容量差值。签订实物合同的用户可自愿申报反向交易报价(即卖电报价)，报价容量不超过其相应时段的实物合同位置。
通过反向交易报价机制，发电侧可以实现低边际成本机组对高边际成本机组的替代发电，用户侧可以降低进入实时市场的偏差电量，而调度机构则可以实现电网的阻塞管理。
(3)引入需求侧灵活报价
与用户侧的反向交易报价不同，需求侧灵活报价只申报负荷可下调容量、价格和可调节时长，但并不指定具体投标时段，而是由日前市场出清模型根据社会福利最大化原则确定其最优中标时段。需求侧灵活报价可以满足可调节负荷参与日前市场的需要。
(4)设置市场供需平衡松弛变量
日前市场出清模型中，当自主开机的机组最小出力之和超过低谷某些时段的系统负荷需求时，引入市场供需平衡松弛变量保证模型求解。
日前市场出清后、实时市场运行前，由调度机构根据系统平衡需要决定是否购买停机服务，引起的费用可以从相应时段的日前市场电费结算盈余中支付，不足部分再由日前市场相应时段存在正偏差电量的机组和实时市场相应时段存在负偏差电量的用户按偏差电量占比共同承担。
1.2、目标函数
模型以社会福利最大化为优化目标，可以对日前市场的电能与旋转备用进行联合优化出清，目标函数的具体数学表达式如下：
Figure BDA0002768201870000091
式中：△t为时段间隔；T为时段总数；NG、ND分别为机组和负荷的总数；
Figure BDA0002768201870000092
分别为t时段第i个机组和第j个负荷买电报价的数目；
Figure BDA0002768201870000093
为t时段第i个机组第k个买电报价的价格；
Figure BDA0002768201870000094
为t时段第i个机组第k个买电报价的中标电力；
Figure BDA0002768201870000095
为t时段第j个负荷第w个买电报价的价格；
Figure BDA0002768201870000096
为t时段第j个负荷第w个买电报价的中标电力；
Figure BDA0002768201870000097
分别为为t时段第i个机组和第j个负荷卖电报价的数目；
Figure BDA0002768201870000098
为t时段第i个机组第h个卖电报价的价格；
Figure BDA0002768201870000099
为t时段第i个机组第h个卖电报价的中标电力；
Figure BDA00027682018700000910
为t时段第j个负荷第v个卖电报价的价格；
Figure BDA00027682018700000911
为t时段第j个负荷第v个卖电报价的中标电力；
Figure BDA00027682018700000912
为t时段第i个机组中标的旋转备用容量；
Figure BDA00027682018700000913
为t时段第i个机组上旋转备用的报价；yi,t为启机变量，StartUpi为机组i的启机成本；
Figure BDA00027682018700000914
为第m个需求侧灵活报价在t时段的调用状态0-1变量，
Figure BDA00027682018700000915
为可中断负荷的报价；λ为供需平衡松弛变量的惩罚因子，δt为t时段供需不平衡电量；
1.3、约束条件
1)机组和负荷中标电力的上下限约束
Figure BDA00027682018700000916
式中：
Figure BDA00027682018700000917
为t时段第i个机组第k个买电报价的容量上限；
Figure BDA00027682018700000918
为t时段第i个机组第h个卖电报价的容量上限；
Figure BDA00027682018700000919
为t时段第j个负荷第w个买电报价的容量上限；
Figure BDA00027682018700000920
为t时段第j个负荷第v个卖电报价的容量上限；
2)机组上旋转备用中标容量的上下限约束
Figure BDA00027682018700000921
式中，
Figure BDA00027682018700000922
为t时段机组i可提供的旋转备用容量上限；
3)机组出力上下限约束
Figure BDA0002768201870000101
式中：
Figure BDA0002768201870000102
为t时段机组i的实物合同位置；
Figure BDA0002768201870000103
为机组已经t时段机组i通过中长期合同已售出的旋备容量；ui,t为t时段第i个机组的启停机状态0-1变量，对于自主确定启停机的机组，ui,t设为常量1或0，不参与状态变量优化；
Figure BDA0002768201870000104
为机组i的出力上限；
Figure BDA0002768201870000105
为机组i的出力下限；
4)机组出力爬坡约束
Figure BDA0002768201870000106
式中，RUi为机组i在△t时段内的上爬坡速率；RDi为机组i在△t时段内的下爬坡速率；yi,t、zi,t分别为机组启停控制0-1变量，yi,t＝1、zi,t＝1分别表示机组i在第t时段执行启动操作和执行停机操作，否则为0；Pstart,i、Pshut,i分别为启动功率速度限制和停机功率速度限制；
5)启停机逻辑约束
Figure BDA0002768201870000107
当机组运行状态ui,t发生变化时，式(6)可以确保yi,t、zi,t取正确的值；
6)机组最小启停机时间约束
Figure BDA0002768201870000108
式中：e为t时段启停机约束校验的起始时段；Gi为机组i在初始时刻仍需连续运行的时段数；Li为机组i在初始时刻仍需连续停机的时段数；
Figure BDA0002768201870000109
为机组i的最小运行时间；
Figure BDA00027682018700001010
为机组i的最小停机时间；
7)需求侧灵活报价运行约束
Figure BDA00027682018700001011
Tm为第m个需求侧灵活报价的可调节时长；
8)系统电力平衡约束
Figure BDA0002768201870000111
9)系统旋转备用需求约束
Figure BDA0002768201870000117
式中：
Figure BDA0002768201870000112
为t时段系统总负荷需求预测；
Figure BDA0002768201870000113
为系统在t时段的旋转备用总需求；旋转备用需求可以由市场运营机构以中长期合同方式购买，不足部分可以在日前市场中购买，RSt为系统在t时段需要从日前市场中购买的旋转备用需求；
10)网络传输约束
Figure BDA0002768201870000114
式中：
Figure BDA0002768201870000115
为第l条线路的潮流传输极限；模型中潮流的求解采用直流潮流算法，Se,l为第e个节点注入功率对第l条线路潮流的灵敏度；Pe,net为第e个节点的净注入功率；
Figure BDA0002768201870000116
2、算例仿真
为验证本发明机制及模型的有效性，采用图1所示的IEEE-30节点标准测试系统进行仿真分析，模型采用GAMS编程求解。各机组、各节点负荷24时段的中长期实物合同数据及日前市场报价数据详见附表。
日前市场的定价机制一般分为系统边际电价、分区边际电价和节点边际电价三种，适用不同阻塞情况的电网结构。为便于分析，本发明算例不对电网模型进行简化，直接采用节点边际电价进行分析。
基于图1所示的仿真系统，对以下四种场景下的日前市场出清结果进行分析：(1)网络无输电容量约束，6台机组全部决定开机运行，即场景1。(2)网络无输电容量约束，机组5决定停机运行，即场景2。(3)网络存在输电容量约束，6台机组全部决定开机运行，即场景3。(4)网络存在输电容量约束，机组5决定停机运行，即场景4。
2.1、仿真结果
四种场景中各节点的负荷报价数据、开机机组的报价数据、双边合同数据、负荷预测数据等保持不变，各场景下的日前市场SCUC优化结果、发电侧各机组中标结果及各节点边际出清价格分别如图2、图3和图4所示。
如图4所示，场景1和场景2均无输电容量的约束，同一时段各节点之间的市场出清价格相同，实质上为系统边际电价。相比场景1，场景2中由于5号机组选择停机运行，机组5在各时段的双边合同电量需要由其他机组进行替代发电，导致部分时段(如时段7、时段11、时段14、时段19)的边际机组转移至更高报价的机组或容量区间，相应时段的市场出清价格要高于场景1。
场景1中，时段3和时段4的市场出清价格为日前市场最低限价-300元/MWh，这是由于场景1中六台机组均选择开机运行，开机机组在负荷低谷的时段3和时段4最小出力之和超过了系统负荷需求，模型求解引入了供需平衡松弛变量，市场出清价格默认设定为价格下限。而场景2中，由于机组5选择停机运行，系统负荷需求能满足开机机组最小运行出力的需要，时段3和时段4没有出现场景1中的负市场出清价。
场景3和场景4中的部分时段存在输电线路容量越限情况，受约束影响的节点市场出清价格出现差异，且价格高于场景1、场景2。此外，场景4下由于2号机组选择停机运行，导致部分时段边际机组转移至更高报价的机组或容量区间，相应时段部分节点的市场出清价格要高于场景3(例如，场景3中节点2在时段14的市场出清价格为240元/MWh，而在场景4的市场出清价格则为266.30元/MWh)。
2.2、结果分析
算例仿真表明，本发明日前市场模型能够满足不同运行场景的需要，实现基于机组运行约束的中长期实物合同调整、基于阻塞管理的中长期实物合同调整以及基于帕累托改进的发电权转让交易，具有很好的现实可行性和适应性。
(1)基于机组运行约束的中长期实物合同调整。以场景2中的时段4为例，六台机组的最小运行出力要求分别为50MW、20MW、15MW、10MW、10MW、12MW，各机组在时段4的双边合同电力分别为60MW、12MW、10MW、5MW、4MW、8MW。除机组1外，其他机组均无法满足最小运行出力需要。日前市场出清后，负荷侧中标净值12MW，机组5在时段4中标-4MW实现停机运行，机组2、机组3、机组4、机组6在时段4分别中标8MW、5MW、5MW、4MW满足最小出力运行，而机组1中标-6MW满足系统平衡需要。时段4的市场出清价格为180元/MW，为机组1在50MW至60MW容量段的反向交易报价，也即机组1的报价决定了市场出清价格，其他机组因运行约束而强制中标的电量只能作为价格接受者。
(2)基于帕累托改进的发电权转让交易。以场景2中的时段8为例，六台机组在时段4的双边合同电力分别为120MW、30MW、25MW、12MW、10MW、18MW，均能满足最小运行出力需要。其中，机组1在120MW-170MW、170-190MW、190MW-200MW容量区间的卖电报价分别为200元/MWh、230元/MWh、270元/MWh，机组2在20MW-30MW区间的买电报价为205元/MWh，机组3在15MW-25MW区间的买电报价为210元/MWh，机组4在10MW-12MW区间的买电报价为245元/MWh，机组6在12MW-18MW区间的买电报价为235元/MWh。机组2、机组3、机组4和机组6的买电报价高于机组1在第一段容量区间的卖电报价，在经济性上可以进行替代发电实现帕累托改进。日前市场出清后，负荷侧在时段8中标净值13MW，机组5中标-10MW实现停机运行，机组2标-9MW运行出力点降至21MW，机组3、机组4、机组6分别中标-10MW、-2MW、-6MW降至各自最小出力运行，而机组1中标50MW满足系统供需平衡，系统出清价格为205元/MWh。
(3)基于阻塞管理的中长期实物合同调整。当日前市场发生基于阻塞管理的实物合同调整时，系统必然处于线路传输受限的情形，市场出清出现节点电价。由于基于阻塞管理的实物合同调整与基于帕累托改进的发电权交易在日前市场出清优化时往往同时发生，因而在图1所示的多节点环网仿真系统中难以单独分析其机理，本发明采用如图5所示的双节点系统进行阐释，线路传输容量为100MW。
图5中，机组G1和机组G2与负荷L在T时段分别签订了120MW、60MW的双边实物合同。负荷L预测该时段的用电需求预测值为190MW，日前市场报买电报价10MW、300元/MWh；机组G1在日前市场的卖电报价为30MW(120MW-150MW区间)、250元/MWh，买电报价为70MW(50MW-120MW区间)、180元/MWh；机组G2在日前市场的卖电报价为40MW(60MW-100MW区间)、280元/MWh，买电报价为10MW(50MW-60MW区间)、220元/MWh。日前市场考虑线路输电约束出清后，G1中标-20MW，所在节点出清价格为180元/MWh；G2中标30MW，所在节点出清价格为280元/MWh，该时段的日前市场电费结算盈余为-2000元，即阻塞管理成本为2000元。
2.3、电费结算
日前市场出清后，中长期实物合同对应电力电量仍按照合同约定价格结算，中长期实物合同之外的电量按照日前市场形成的各时段节点边际电价进行结算。如表1所示，四种场景下日前市场的电费结算盈亏分别为3300元、0元、10983.67元、8995.54元。
表1日前市场电费结算
Figure BDA0002768201870000131
Figure BDA0002768201870000141
场景1中，由于时段3和时段4的拟开机机组最小出力之和超过了系统负荷需求，市场出清价格设为价格下限(-300元/MWh)，最小出力高于双边合同电力值的机组(机组2、机组3、机组4和机组5)，其溢出电量以负价格卖出，导致电费结算盈余为正值。
场景2中，机组5决定停机运行，各时段均能实现供需平衡，而且买卖双方均按系统边际电价进行结算，电费结算盈余为零。
场景3和场景4中，部分时段发生阻塞，节点之间存在价差，结算后产生阻塞盈余。场景3和场景4中的日前市场电费结算盈余包括因双边合同阻塞管理引起的成本，当阻塞管理成本较高时，会出现电费结算盈余为负的情形。
2.4、模型对比
以IEEE30系统为例，选取场景4，从市场出清结果和社会福利两方面对比分析本文模型与欧洲现行模型的差异。模型对比基于以下共同前提：(1)日前市场只开展电能量交易；(2)市场主体在两种模型中的报价策略保持不变；(3)负荷的实时用电需求与日前预测相同；(4)实时平衡机制以调整成本最小化为目标，按报价结算平衡服务电量。两种模型的出清结果分别见表2和表3，表中数据为24时段的累计数据，实时调整成本为实时平衡服务的结算资金亏空。
表2欧洲现行模型出清结果
Figure BDA0002768201870000142
表3本发明模型出清结果
Figure BDA0002768201870000143
Figure BDA0002768201870000151
如表2和表3所示，两种模型的市场出清结果存在较大差异。现行模型的日前市场中标电量主要由价格较低的机组1获得(中标1026MWh)；而本发明模型中，由于考虑了中长期实物合同交割对电网运行的影响，线路传输能力限制了机组1的中标电量(仅中标663.15MWh)，机组2和机组3的中标电量则相应增加。日前市场中标电量的不同，导致两种模型形成的日前发电计划也存在显著差异。当负荷侧实时用电需求与日前预测相同时，现行模型形成的日前发电计划在实时运行阶段需要购买较多的平衡服务以满足电网运行约束，而本发明模型形成的日前发电计划无需调用平衡服务。从最终实时发电量的角度，在市场主体报价策略一致的前提下，两种模型下各机组最终的实时发电量相差很小，引起差别的原因主要在于日前发电计划的不同导致了平衡服务报价的容量分段不完全一致。
表4社会福利比较
Figure BDA0002768201870000152
如表4所示，本发明日前市场模型中，报价较低的机组受电网运行约束的限制不一定能够完全中标，日前市场的社会福利小于欧洲现行模型。但是，在实时运行阶段，欧洲现行模型形成的日前发电计划产生了较高的调整成本，无论由发电侧还是用户侧承担这部分成本均会造成社会福利的损失。综合考虑日前和实时两阶段整体的社会福利，本发明模型要优于欧洲现行模型。
综上可见，本发明设计的日前市场出清模型，可以实现日前现货交易与中长期实物合同交割、电网阻塞管理的联合优化，与我国的电力交易与调管理体制具有良好的衔接性。多场景的仿真算例表明本文模型能够满足不同市场运行场景的需要，具有较强的可行性和灵活性。与欧洲分散式市场模式现行日前市场出清模型的算例对比表明，本发明模型能够有效提高日前市场出清结果的可执行性，减少实时平衡机制的调整成本，提高日前和实时两阶段整体的社会福利。需要强调的是，对于电网阻塞比较严重的地区，在组织中长期实物交易时应考虑关键的输电约束，降低日前市场的阻塞管理成本。
以上是本发明的较佳实施例，凡依本发明技术方案所作的改变，所产生的功能作用未超出本发明技术方案的范围时，均属于本发明的保护范围。

Claims (1)

1.一种兼容中长期实物合同的日前市场出清模型构建方法，其特征在于，包括如下步骤：
步骤S1、构建兼容中长期实物合同的日前市场出清模型，该模型以社会福利最大化为优化目标，可对日前市场的电能与旋转备用进行联合优化出清，目标函数的具体数学表达式如下：
Figure FDA0002768201860000011
式中：△t为时段间隔；T为时段总数；NG、ND分别为机组和负荷的总数；
Figure FDA0002768201860000012
分别为t时段第i个机组和第j个负荷买电报价的数目；
Figure FDA0002768201860000013
为t时段第i个机组第k个买电报价的价格；
Figure FDA0002768201860000014
为t时段第i个机组第k个买电报价的中标电力；
Figure FDA0002768201860000015
为t时段第j个负荷第w个买电报价的价格；
Figure FDA0002768201860000016
为t时段第j个负荷第w个买电报价的中标电力；
Figure FDA0002768201860000017
分别为为t时段第i个机组和第j个负荷卖电报价的数目；
Figure FDA0002768201860000018
为t时段第i个机组第h个卖电报价的价格；
Figure FDA0002768201860000019
为t时段第i个机组第h个卖电报价的中标电力；
Figure FDA00027682018600000110
为t时段第j个负荷第v个卖电报价的价格；
Figure FDA00027682018600000111
为t时段第j个负荷第v个卖电报价的中标电力；
Figure FDA00027682018600000112
为t时段第i个机组中标的旋转备用容量；
Figure FDA00027682018600000113
为t时段第i个机组上旋转备用的报价；yi,t为启机变量，StartUpi为机组i的启机成本；
Figure FDA00027682018600000114
为第m个需求侧灵活报价在t时段的调用状态0-1变量，
Figure FDA00027682018600000115
为可中断负荷的报价；λ为供需平衡松弛变量的惩罚因子，δt为t时段供需不平衡电量；
步骤S2、设置约束条件：
1)机组和负荷中标电力的上下限约束
Figure FDA00027682018600000116
式中：
Figure FDA00027682018600000117
为t时段第i个机组第k个买电报价的容量上限；
Figure FDA00027682018600000118
为t时段第i个机组第h个卖电报价的容量上限；
Figure FDA00027682018600000119
为t时段第j个负荷第w个买电报价的容量上限；
Figure FDA00027682018600000120
为t时段第j个负荷第v个卖电报价的容量上限；
2)机组上旋转备用中标容量的上下限约束
Figure FDA00027682018600000121
式中，
Figure FDA00027682018600000122
为t时段机组i可提供的旋转备用容量上限；
3)机组出力上下限约束
Figure FDA0002768201860000021
式中：
Figure FDA0002768201860000022
为t时段机组i的实物合同位置；
Figure FDA0002768201860000023
为机组已经t时段机组i通过中长期合同已售出的旋备容量；ui,t为t时段第i个机组的启停机状态0-1变量，对于自主确定启停机的机组，ui,t设为常量1或0，不参与状态变量优化；
Figure FDA0002768201860000024
为机组i的出力上限；
Figure FDA0002768201860000025
为机组i的出力下限；
4)机组出力爬坡约束
Figure FDA0002768201860000026
式中，RUi为机组i在△t时段内的上爬坡速率；RDi为机组i在△t时段内的下爬坡速率；yi,t、zi,t分别为机组启停控制0-1变量，yi,t＝1、zi,t＝1分别表示机组i在第t时段执行启动操作和执行停机操作，否则为0；Pstart,i、Pshut,i分别为启动功率速度限制和停机功率速度限制；
5)启停机逻辑约束
Figure FDA0002768201860000027
当机组运行状态ui,t发生变化时，式(6)可以确保yi,t、zi,t取正确的值；
6)机组最小启停机时间约束
Figure FDA0002768201860000028
式中：e为t时段启停机约束校验的起始时段；Gi为机组i在初始时刻仍需连续运行的时段数；Li为机组i在初始时刻仍需连续停机的时段数；
Figure FDA0002768201860000029
为机组i的最小运行时间；
Figure FDA00027682018600000210
为机组i的最小停机时间；
7)需求侧灵活报价运行约束
Figure FDA00027682018600000211
Tm为第m个需求侧灵活报价的可调节时长；
8)系统电力平衡约束
Figure FDA0002768201860000031
9)系统旋转备用需求约束
Figure FDA0002768201860000032
式中：
Figure FDA0002768201860000033
为t时段系统总负荷需求预测；
Figure FDA0002768201860000034
为系统在t时段的旋转备用总需求；旋转备用需求可以由市场运营机构以中长期合同方式购买，不足部分可以在日前市场中购买，RSt为系统在t时段需要从日前市场中购买的旋转备用需求；
10)网络传输约束
Figure FDA0002768201860000035
式中：
Figure FDA0002768201860000036
为第l条线路的潮流传输极限；模型中潮流的求解采用直流潮流算法，Se,l为第e个节点注入功率对第l条线路潮流的灵敏度；Pe,net为第e个节点的净注入功率；
Figure FDA0002768201860000037
CN202011237229.9A 2020-11-09 2020-11-09 兼容中长期实物合同的日前市场出清模型构建方法 Active CN112308329B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011237229.9A CN112308329B (zh) 2020-11-09 2020-11-09 兼容中长期实物合同的日前市场出清模型构建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011237229.9A CN112308329B (zh) 2020-11-09 2020-11-09 兼容中长期实物合同的日前市场出清模型构建方法

Publications (2)

Publication Number Publication Date
CN112308329A true CN112308329A (zh) 2021-02-02
CN112308329B CN112308329B (zh) 2022-08-05

Family

ID=74326487

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011237229.9A Active CN112308329B (zh) 2020-11-09 2020-11-09 兼容中长期实物合同的日前市场出清模型构建方法

Country Status (1)

Country Link
CN (1) CN112308329B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114117326A (zh) * 2021-12-06 2022-03-01 国网河南省电力公司经济技术研究院 基于系统安全运行约束的微电网市场双阶段交易优化机制

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050114255A1 (en) * 1999-10-20 2005-05-26 Accenture Global Services Gmbh Spot market clearing
CN105447775A (zh) * 2016-01-11 2016-03-30 广东电网有限责任公司 内嵌电网安全约束的电力日前市场撮合交易方法
CN108711890A (zh) * 2018-06-27 2018-10-26 广东电网有限责任公司 日前市场出清方法、系统、装置及计算机可读存储介质
CN110110984A (zh) * 2019-04-26 2019-08-09 国网福建省电力有限公司 一种考虑用户需求响应的电力现货市场出清方法
CN110852535A (zh) * 2018-07-24 2020-02-28 华北电力大学 一种考虑中长期交易和风电不确定性的日前市场出清模型
CN111832824A (zh) * 2020-07-14 2020-10-27 国网电力科学研究院有限公司 电力调频市场交易出清及结算方法、装置及系统
CN111860994A (zh) * 2020-07-14 2020-10-30 北京清能互联科技有限公司 一种考虑发用两侧参与调峰的现货市场交易方法及装置

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050114255A1 (en) * 1999-10-20 2005-05-26 Accenture Global Services Gmbh Spot market clearing
CN105447775A (zh) * 2016-01-11 2016-03-30 广东电网有限责任公司 内嵌电网安全约束的电力日前市场撮合交易方法
CN108711890A (zh) * 2018-06-27 2018-10-26 广东电网有限责任公司 日前市场出清方法、系统、装置及计算机可读存储介质
CN110852535A (zh) * 2018-07-24 2020-02-28 华北电力大学 一种考虑中长期交易和风电不确定性的日前市场出清模型
CN110110984A (zh) * 2019-04-26 2019-08-09 国网福建省电力有限公司 一种考虑用户需求响应的电力现货市场出清方法
CN111832824A (zh) * 2020-07-14 2020-10-27 国网电力科学研究院有限公司 电力调频市场交易出清及结算方法、装置及系统
CN111860994A (zh) * 2020-07-14 2020-10-30 北京清能互联科技有限公司 一种考虑发用两侧参与调峰的现货市场交易方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张粒子 等: "兼容中长期实物合同的日前市场出清模型", 《电力系统自动化》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114117326A (zh) * 2021-12-06 2022-03-01 国网河南省电力公司经济技术研究院 基于系统安全运行约束的微电网市场双阶段交易优化机制
CN114117326B (zh) * 2021-12-06 2024-07-09 国网河南省电力公司经济技术研究院 基于系统安全运行约束的微电网市场双阶段交易优化机制

Also Published As

Publication number Publication date
CN112308329B (zh) 2022-08-05

Similar Documents

Publication Publication Date Title
CN106845807B (zh) 基于调峰辅助服务的结算方法及装置
Luo et al. Distributed peer-to-peer energy trading based on game theory in a community microgrid considering ownership complexity of distributed energy resources
CN110782281B (zh) 一种多业主梯级电站流域电量转移的日前市场出清方法
Zhang et al. Energy trading with demand response in a community-based P2P energy market
CN110378729B (zh) 一种基于动态能源价格策略的综合需求响应方法
CN108460482A (zh) 考虑可中断负荷和日前小时电价的风电消纳方法
CN112488744A (zh) 一种交易驱动的虚拟电厂定制化构建方法
Fan et al. Research on an optimization model for wind power and thermal power participating in two-level power market transactions
CN113888209A (zh) 虚拟电厂参与电力市场和碳交易市场的协同竞价方法
CN116231655A (zh) 一种计及源-荷多类型备用的虚拟电厂双层优化调度方法
CN115062835A (zh) 基于虚拟发电厂的主动配电网分布式资源优化调度方法
Hongliang et al. Virtual power plant participates in the two-level decision-making optimization of internal purchase and sale of electricity and external multi-market
CN112308329B (zh) 兼容中长期实物合同的日前市场出清模型构建方法
CN114912943A (zh) 一种考虑碳排放交易的虚拟电厂中长期联合调度方法
Ghaemi et al. Local energy communities with strategic behavior of multi-energy players for peer-to-peer trading: A techno-economic assessment
Aminlou et al. ADMM-based fully decentralized peer to peer energy trading considering a shared CAES in a local community
CN113690875A (zh) 一种微网实时互动等值模型建立方法
CN117391718A (zh) 一种基于动态减排因子的绿电-ccer互认交易系统
CN110556821B (zh) 考虑交互功率控制和双边竞价交易的多微电网双层优化调度方法
CN117114877A (zh) 一种基于虚拟电厂的中长期电力交易方法及系统
CN108418210A (zh) 提高风电消纳能力的可中断负荷机制设计优化方法
Yang et al. The internal resource aggregation strategy of virtual power plants based on different basic configurations
CN112381399A (zh) 一种能源系统日前-实时出清方法及系统
Ma et al. Optimal peer-to-peer energy transaction for distributed prosumers
Hui et al. Analysis and Research on the Key Issues of Electricity Market Structure and Transmission Congestion Based on Financial Transmission Rights

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