CN107833101B - 一种基于节点电价机制的可再生能源交易费用的计算方法 - Google Patents

一种基于节点电价机制的可再生能源交易费用的计算方法 Download PDF

Info

Publication number
CN107833101B
CN107833101B CN201711107239.9A CN201711107239A CN107833101B CN 107833101 B CN107833101 B CN 107833101B CN 201711107239 A CN201711107239 A CN 201711107239A CN 107833101 B CN107833101 B CN 107833101B
Authority
CN
China
Prior art keywords
renewable energy
time
power
output
real
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
CN201711107239.9A
Other languages
English (en)
Other versions
CN107833101A (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.)
Northwest Branch Of State Grid Power Grid Co
Xian Jiaotong University
Original Assignee
Northwest Branch Of State Grid Power Grid Co
Xian Jiaotong University
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 Northwest Branch Of State Grid Power Grid Co, Xian Jiaotong University filed Critical Northwest Branch Of State Grid Power Grid Co
Priority to CN201711107239.9A priority Critical patent/CN107833101B/zh
Publication of CN107833101A publication Critical patent/CN107833101A/zh
Application granted granted Critical
Publication of CN107833101B publication Critical patent/CN107833101B/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
    • 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
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S50/00—Market activities related to the operation of systems integrating technologies related to power network operation or related to communication or information technologies
    • Y04S50/10—Energy trading, including energy flowing from end-user application to grid

Landscapes

  • Business, Economics & Management (AREA)
  • Economics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Finance (AREA)
  • Accounting & Taxation (AREA)
  • Health & Medical Sciences (AREA)
  • Marketing (AREA)
  • General Physics & Mathematics (AREA)
  • Strategic Management (AREA)
  • General Business, Economics & Management (AREA)
  • Public Health (AREA)
  • Tourism & Hospitality (AREA)
  • Primary Health Care (AREA)
  • Human Resources & Organizations (AREA)
  • General Health & Medical Sciences (AREA)
  • Water Supply & Treatment (AREA)
  • Development Economics (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种基于节点电价机制的可再生能源交易费用的计算方法，将可再生能源出力分阶段引入最优直流潮流模型中，建立基于节点电价机制的最优直流潮流计算模型，通过计算考虑可再生能源日前预测出力和可再生能源实时出力偏差的最优直流潮流，得到可再生能源的实际交易费用。所得结果能够考虑可再生能源出力预测误差的影响，能为高可再生能源渗透率下电网的运行方式调整提供优化方案，尤其适用于大规模、多类型可再生能源共存的电力交易费用计算。

Description

一种基于节点电价机制的可再生能源交易费用的计算方法
技术领域
本发明属于电力系统领域，具体涉及一种基于节点电价机制的可再生能源交易费用的计算方法。
背景技术
随着我国可再生能源发展力度的加大，可再生能源发电在电力系统中的渗透率不断增长，其对电力系统运行造成的不利影响也越来越严重。该影响不只体现在对电网安全性和可靠性的影响，也体现在对电力市场的影响。在含可再生能源参与的电力市场中，如何结算具有不确定性特征的可再生能源出力的交易费用，成为亟待解决的重要问题之一。
我国目前的可再生能源交易费用计算方法通常采用固定的可再生能源标杆电价对可再生能源的实际出力直接进行结算。现有方法存在一些缺陷，其并未考虑可再生能源发电机组对电力系统运行安全性的影响，对所有可再生能源出力按照统一固定价格进行结算，难以引导可再生能源发电方提高自身出力预测精度并改善出力质量。实际上，由于电力系统运行本身受到输电容量、调峰容量等技术约束限制，当大规模的可再生能源参与交易时，如果可再生能源出力的预测精度较低使得可再生能源实际出力偏差较大，往往会导致上述技术约束的越限，给电力系统的安全稳定运行带来了严重的影响。
发明内容
为了解决上述问题，本发明的目的在于提供一种能够引导可再生能源发电方降低自身出力预测误差的基于节点电价机制的可再生能源交易费用的计算方法。
为达到上述目的，本发明采用的技术方案是，一种基于节点电价机制的可再生能源交易费用的计算方法，包括以下步骤：
1)获取电力系统基本技术数据、电力系统运行约束基础数据、电力系统实时运行数据、非可再生能源发电企业报价数据以及可再生能源发电企业报价数据；
电力系统基本技术数据为节点数据、输电线路数据、变压器数据和负荷数据；
电力系统运行约束基础数据为各非可再生能源发电机组有功出力的上限值和下限值、输电线路的传输容量限值、可再生能源发电机组日前预测出力值；
电力系统实时运行数据为负荷实时偏差值、可再生能源发电机组实时预测出力值、各非可再生能源发电机组可提供的实时有功调节能力的上限值和下限值；
非可再生能源发电企业报价数据为非可再生能源发电企业对其所属非可再生能源发电机组出力费用的日前报价值，以及非可再生能源发电机组参与系统实时调节的实时报价值；
可再生能源发电企业报价数据为可再生发电企业对其所属可再生能源机组出力费用的日前报价值，以及可再生能源出力实时正偏移量费用的实时报价值；
2)以电力系统发电成本最小为目标函数构建最优直流潮流优化目标；然后构建约束条件，完成面向日前节点电价结算的最优直流潮流模型的构建；其中，约束条件包括系统功率平衡约束、直流潮流约束以及非可再生能源发电机组出力约束；
3)设定面向日前节点电价结算的最优直流潮流模型中的收敛准则参数，收敛准则参数包括优化模型收敛精度和最大迭代次数；将电力系统基本技术数据、电力系统运行约束基础数据、非可再生能源发电企业报价数据、可再生能源发电企业报价数据以及收敛准则参数输入到面向日前节点电价结算的最优直流潮流模型，对计及可再生能源日前预测出力的日前场景进行最优直流潮流分析求解，求解结果包括各节点的节点电价和各非可再生能源发电机组的有功出力；
4)以电力系统综合运行成本最小为目标函数构建最优直流潮流优化目标；然后构建约束条件，完成面向实时调节费用修正的最优直流潮流模型的构建；其中，约束条件包括系统功率平衡约束、直流潮流约束、可再生能源发电机组出力约束，非可再生能源发电机组调节能力约束；
5)设定面向实时调节费用修正的最优直流潮流模型中的收敛准则参数，收敛准则参数包括优化模型收敛精度和最大迭代次数；将电力系统基本技术数据、电力系统运行约束基础数据、电力系统实时运行数据、非可再生能源发电企业报价数据、可再生能源发电企业报价数据、日前场景各非可再生能源发电机组的有功出力结果以及收敛准则参数输入到面向实时调节费用修正的最优直流潮流模型，对计及可再生能源实时出力偏移量的实时场景进行最优直流潮流分析求解，求解结果包括各非可再生能源发电机组的有功调节量、各可再生能源发电机组有功出力的实时偏移量；
6)根据步骤1)得到的日前场景下的可再生能源发电机组的日前预测出力和步骤3)得到的节点电价结果，得到可再生能源交易的日前费用；获取实时场景下可再生能源出力的偏移量，以及非可再生能源发电机组的有功调节结果，得到可再生能源交易的实时调节费用；可再生能源交易的日前费用以及实时调节费用之和即为可再生能源交易的最终结算费用；
7)将步骤6)得到的可再生能源交易的最终结算费用反馈至可再生能源发电企业，可再生能源发电企业按照上述求解结果获得交易费用；将计及可再生能源出力实时变化量的实时场景最优直流潮流结果反馈给电网企业，电网企业参考最优直流潮流结果进行全系统实时运行方式调整。
所述的步骤2)中最优直流潮流优化目标是由公式
Figure BDA0001464663290000031
获得的，其中，
Figure BDA0001464663290000032
为非可再生能源发电机组Ci在t时刻的有功出力，
Figure BDA0001464663290000033
为非可再生能源发电机组Ci在t时刻的报价参数，SGC为所有非可再生能源发电机组集合，
Figure BDA0001464663290000034
为可再生能源发电机组Ri在t时刻的日前预测有功出力，
Figure BDA0001464663290000035
为可再生能源发电机组Ri在t时刻的报价参数，SGR为所有可再生能源发电机组集合，ST为所有运行时段的集合。
所述的步骤2)中系统功率平衡约束为：
Figure BDA0001464663290000036
其中，
Figure BDA0001464663290000037
为t时刻的负荷节点Di的有功负荷，SD为所有负荷节点的集合；
直流潮流约束为
Figure BDA0001464663290000038
其中，
Figure BDA0001464663290000039
为位于节点Bi的非可再生能源发电机组Ci在t时刻的有功出力，
Figure BDA00014646632900000310
为位于节点Bi的可再生能源发电机组Ri在t时刻的日前预测有功出力，
Figure BDA00014646632900000311
为位于节点Bi的t时刻的有功负荷；
Figure BDA00014646632900000312
为功率传输分布因子矩阵第l行第Bi列元素，
Figure BDA00014646632900000313
为支路l的传输功率上限值，Fl 为支路l的传输功率下限值，SL为所有支路的集合；
非可再生能源发电机组出力约束为
Figure BDA00014646632900000314
其中
Figure BDA00014646632900000315
为非可再生能源发电机组Ci的有功出力上限值，
Figure BDA00014646632900000316
为非可再生能源发电机组Ci的有功出力下限值。
所述的步骤4)中最优直流潮流优化目标如下公式所示：
Figure BDA0001464663290000041
其中，
Figure BDA0001464663290000042
为可再生能源发电机组Ri在t时刻实时有功出力的负偏移值，当获得可再生能源发电机组Ri在t时刻的实时预测出力
Figure BDA0001464663290000043
时，负偏移值
Figure BDA0001464663290000044
获得；
Figure BDA0001464663290000045
为可再生能源发电机组Ri在t时刻的实时正偏移有功出力，
Figure BDA0001464663290000046
为可再生能源发电机组Ri在t时刻对其实时正偏移有功出力的报价参数；
Figure BDA0001464663290000047
为日前场景下非可再生能源发电机组Ci在t时刻的有功出力结果，
Figure BDA0001464663290000048
和
Figure BDA0001464663290000049
分别为非可再生能源发电机组Ci在t时刻的负有功调节出力和正有功调节出力，
Figure BDA00014646632900000410
和
Figure BDA00014646632900000411
分别为非可再生能源发电机组Ci在t时刻对其负有功调节出力和正有功调节出力的报价参数。
所述的步骤4)中系统功率平衡约束为：
Figure BDA00014646632900000412
其中，
Figure BDA00014646632900000413
为t时刻负荷节点Di的有功负荷实时偏差值；
直流潮流约束为:
Figure BDA00014646632900000414
其中，
Figure BDA00014646632900000415
和
Figure BDA00014646632900000416
分别为位于节点Bi的非可再生能源发电机组Ci在t时刻的负有功调节出力和正有功调节出力；
Figure BDA00014646632900000417
为位于节点Bi的可再生能源发电机组Ri在t时刻实时有功出力的负偏移值，
Figure BDA00014646632900000418
为位于节点Bi的可再生能源发电机组Ri在t时刻的实时正偏移有功出力；
Figure BDA00014646632900000419
为t时刻位于节点Bi的有功负荷的实时偏差值；
可再生能源发电机组出力约束为
Figure BDA00014646632900000420
其中
Figure BDA00014646632900000421
为可再生能源发电机组Ri在t时刻实时有功出力的正偏移值，当获得可再生能源发电机组Ri在t时刻的实时预测出力
Figure BDA00014646632900000422
时，正偏移值
Figure BDA00014646632900000423
由
Figure BDA00014646632900000424
获得；
非可再生能源发电机组调节能力约束为
Figure BDA00014646632900000425
其中
Figure BDA00014646632900000426
和
Figure BDA00014646632900000427
分别为非可再生能源发电机组Ci在t时刻正有功调节出力和负有功调节出力的上限值。
所述的步骤6)中对应可再生能源发电机组Ri的日前费用
Figure BDA0001464663290000051
为：
Figure BDA0001464663290000052
其中，
Figure BDA0001464663290000053
为节点Bi在t时刻的节点电价，单位为元/kWh，节点电价
Figure BDA0001464663290000054
由公式
Figure BDA0001464663290000055
获得，μt对应于t时刻系统功率平衡约束
Figure BDA0001464663290000056
的对偶乘子取值，λlt对应于t时刻支路l直流潮流约束
Figure BDA0001464663290000057
的对偶乘子取值；
对应可再生能源发电机组Ri的实时调节费用
Figure BDA0001464663290000058
为：
Figure BDA0001464663290000059
其中
Figure BDA00014646632900000510
为可再生能源发电机组Ri于t时刻的系统调节费用的分摊系数，
Figure BDA00014646632900000511
由公式
Figure BDA00014646632900000512
获得；
对应可再生能源发电机组Ri的交易费用
Figure BDA00014646632900000513
为：
Figure BDA00014646632900000514
所述非可再生能源发电为火力发电。
与现有技术相比，本发明的有益效果在于：
从系统安全运行方面来讲，本发明在考虑可再生能源交易时利用了最优潮流模型的求解结果，其求解结果能够为电网企业调整运行方式提供参考，同时考虑系统运行约束，将各输电线路的传输功率以及系统调峰容量等限制在允许范围内，能够保证在大规模可再生能源参与市场交易时电力系统安全可靠的运行；从电网发展角度来讲，本发明基于节点电价机制计算可再生能源的日前结算费用，以市场化手段合理的处理了电力系统中存在的网络阻塞问题，能够体现出位于不同地理位置时各可再生能源发电方的实际电能价值，引导可再生能源发电方的合理规划；同时，本发明通过引入系统调节费用修正可再生能源的结算费用，考虑了可再生能源实际出力偏差对系统实时调度的影响，以市场化手段实现了对可再生能源出力预测误差的考核机制，有利于引导可再生能源发电方提高自身出力预测精度；从经济效益方面来讲，本发明通过引入系统调节费用，在保证安全运行前提下合理的促进了可再生能源的消纳，由于可再生能源的低成本性，实现了电力系统运行成本的降低，提高了其经济性；从环保角度来讲，本发明间接实现了可再生能源和传统能源之间的发电权交易，减少传统火电机组的发电量，降低了环保成本，实现了环保效益的提高，符合可持续发展的战略要求。
附图说明
图1为本发明的流程图。
具体实施方式
下面结合附图对本发明作进一步说明，但本发明的内容不仅仅局限于此。
参见图1，本发明包括基于节点电价机制的可再生能源交易费用的计算方法包括以下步骤：
(一)应用本发明所提模型时，需要获取计算模型输入数据包括：电力系统基本技术数据、电力系统运行约束基础数据、电力系统实时运行数据、非可再生能源发电企业报价数据以及可再生能源发电企业报价数据；
其中，电力系统基本技术数据为节点数据、输电线路数据、变压器数据、负荷数据；节点数据如节点编号Bi和系统平衡节点的编号Bs；输电线路数据如支路l的电抗值Xl；变压器数据如变压器TRi的阻抗值
Figure BDA0001464663290000061
负荷数据如负荷节点Di在t时刻的负荷值
Figure BDA0001464663290000062
电力系统运行约束基础数据包括非可再生能源发电机组约束数据如机组Ci有功出力的上限值
Figure BDA0001464663290000063
和下限值
Figure BDA0001464663290000064
输电线路约束数据如支路l传输容量的上限值
Figure BDA0001464663290000065
和下限值Fl ；可再生能源发电机组约束数据如机组Ri在t时刻的日前预测保障出力值
Figure BDA0001464663290000066
电力系统实时运行数据包括有功负荷的实时偏差数据如t时刻负荷节点Di的有功负荷偏差值
Figure BDA0001464663290000067
可再生能源发电机组实时预测出力数据如t时刻机组Ri的实时预测有功出力值
Figure BDA0001464663290000068
各非可再生能源发电机组调节能力数据如t时刻机组Ci可提供实时有功调节能力的上限值
Figure BDA0001464663290000069
和下限值
Figure BDA00014646632900000610
非可再生能源发电企业报价数据包括非可再生能源发电机组出力的日前报价数据如t时刻机组Ci对其出力
Figure BDA00014646632900000611
的日前报价值
Figure BDA00014646632900000612
非可再生能源发电机组参与实时调节的报价数据如t时刻机组Ci的对其负有功调节出力
Figure BDA00014646632900000613
和正有功调节出力
Figure BDA00014646632900000614
的报价值
Figure BDA00014646632900000615
和
Figure BDA00014646632900000616
可再生能源发电企业报价数据包括可再生能源机组出力的日前报价数据如t时刻机组Ri对其日前预测出力
Figure BDA00014646632900000617
的日前报价值
Figure BDA00014646632900000618
可再生能源出力实时偏移量的报价数据如t时刻机组Ri对其正偏移有功出力
Figure BDA00014646632900000619
的报价值
Figure BDA00014646632900000620
其中，l∈Sl，TRi∈STR，Bi∈SB，Ti∈ST，Di∈SD，Ci∈SGC，Ri∈SGR；SB表示电力系统所有节点集合，STR表示所有变压器集合，Sl表示所有支路集合，ST表示所有时段集合，SD表示所有负荷节点集合，SGC表示所有非可再生能源发电机组集合，SGR表示所有可再生能源发电机组集合，下同。
(二)从相关部门获取上述信息后，按照下述步骤进行可再生能源交易费用的计算：
第1步：以电力系统发电成本最小为目标函数建立优化目标，所述优化目标由下式获得：
Figure BDA0001464663290000071
其中，
Figure BDA0001464663290000072
为非可再生能源发电机组Ci在t时刻的有功出力，
Figure BDA0001464663290000073
为非可再生能源发电机组Ci在t时刻的报价参数，
Figure BDA0001464663290000074
为可再生能源发电机组Ri在t时刻的日前预测有功出力，
Figure BDA0001464663290000075
为可再生能源发电机组Ri在t时刻的报价参数。
第2步：构建系统运行约束条件完成面向日前节点电价结算最优潮流模型的构建；其中，约束条件包括系统功率平衡约束、直流潮流约束以及非可再生能源发电机组出力约束；
系统功率平衡约束中整个系统的总发电量和总负荷量必须满足平衡等式约束为：
Figure BDA0001464663290000076
式中，
Figure BDA0001464663290000077
为非可再生能源发电机组Ci在t时刻的有功出力，
Figure BDA0001464663290000078
为可再生能源发电机组Ri在t时刻的日前预测有功出力，
Figure BDA0001464663290000079
为t时刻负荷节点Di的有功负荷；
直流潮流约束为
Figure BDA00014646632900000710
式中，
Figure BDA00014646632900000711
为位于节点Bi的非可再生能源发电机组Ci在t时刻的有功出力，
Figure BDA00014646632900000712
为位于节点Bi的可再生能源发电机组Ri在t时刻的日前预测有功出力，
Figure BDA00014646632900000713
为t时刻位于节点Bi的有功负荷；
Figure BDA00014646632900000714
为功率传输分布因子矩阵第l行第Bi列元素，
Figure BDA00014646632900000715
为支路l的传输功率上限值，Fl为支路l的传输功率下限值；
非可再生能源发电机组出力约束为
Figure BDA00014646632900000716
其中
Figure BDA00014646632900000717
为非可再生能源发电机组Ci的有功出力上限值，
Figure BDA00014646632900000718
为非可再生能源发电机组Ci的有功出力下限值。
第3步：设定面向日前节点电价结算的最优直流潮流模型中的收敛准则参数，收敛准则参数包括优化模型收敛精度和最大迭代次数；将电力系统基本技术数据、电力系统运行约束基础数据、非可再生能源发电企业报价数据、可再生能源发电企业报价数据以及收敛准则参数输入到面向日前节点电价结算的最优直流潮流模型。
对计及可再生能源日前预测出力的日前场景进行最优直流潮流分析求解，求解t时刻节点Bi的节点电价结果
Figure BDA0001464663290000081
以及发电机组Ci的有功出力结果
Figure BDA0001464663290000082
其中，节点电价
Figure BDA0001464663290000083
的单位为元/kWh，其由公式
Figure BDA0001464663290000084
获得，μt对应于t时刻系统功率平衡约束
Figure BDA0001464663290000085
的对偶乘子取值，λlt对应于t时刻支路l直流潮流约束
Figure BDA0001464663290000086
的对偶乘子取值。
第4步：以电力系统综合运行成本最小为目标函数建立优化目标，所述优化目标由下式获得：
Figure BDA0001464663290000087
其中，
Figure BDA0001464663290000088
为可再生能源发电机组Ri在t时刻有功出力的实时负偏移值，当获得可再生能源发电机组Ri在t时刻的实时预测出力
Figure BDA0001464663290000089
时，负偏移值
Figure BDA00014646632900000810
由
Figure BDA00014646632900000811
获得，
Figure BDA00014646632900000812
为可再生能源发电机组Ri在t时刻的实时正偏移有功出力，
Figure BDA00014646632900000813
为可再生能源发电机组Ri在t时刻对其实时正偏移有功出力的报价参数；
Figure BDA00014646632900000814
为日前场景下非可再生能源发电机组Ci在t时刻的有功出力结果，
Figure BDA00014646632900000815
和
Figure BDA00014646632900000816
分别为非可再生能源发电机组Ci在t时刻的负有功调节出力和正有功调节出力，
Figure BDA00014646632900000817
和
Figure BDA00014646632900000818
分别为非可再生能源发电机组Ci在t时刻对其负有功调节出力和正有功调节出力的报价参数。
第5步：构建系统运行约束条件完成面向实时调节费用修正的最优直流潮流模型的构建；其中，约束条件包括系统功率平衡约束、直流潮流约束、可再生能源发电机组出力约束，非可再生能源发电机组调节能力约束。
系统功率平衡约束中整个系统的总发电量、总负荷量、总偏移量以及总调节量必须满足平衡等式约束为：
Figure BDA0001464663290000091
其中，
Figure BDA0001464663290000092
为t时刻有功负荷的实时偏差值；
直流潮流约束为:
Figure BDA0001464663290000093
其中，
Figure BDA0001464663290000094
和
Figure BDA0001464663290000095
分别为位于节点Bi的非可再生能源发电机组Ci在t时刻的负有功调节出力和正有功调节出力；
Figure BDA0001464663290000096
为位于节点Bi的可再生能源发电机组Ri在t时刻有功出力的实时负偏移值，
Figure BDA0001464663290000097
为位于节点Bi的可再生能源发电机组Ri在t时刻的实时正偏移有功出力；
Figure BDA0001464663290000098
为t时刻位于节点Bi的有功负荷的实时偏差值；
可再生能源发电机组出力约束为
Figure BDA0001464663290000099
其中
Figure BDA00014646632900000910
为可再生能源发电机组Ri在t时刻有功出力的实时正偏移值，当获得可再生能源发电机组Ri在t时刻的实时预测出力
Figure BDA00014646632900000911
时，正偏移值
Figure BDA00014646632900000912
由
Figure BDA00014646632900000913
获得；
非可再生能源发电机组调节能力约束为
Figure BDA00014646632900000914
其中
Figure BDA00014646632900000915
和
Figure BDA00014646632900000916
分别为非可再生能源发电机组Ci在t时刻正有功调节出力和负有功调节出力的上限值。
第6步：设定面向实时调节费用修正的最优直流潮流模型中的收敛准则参数，收敛准则参数包括优化模型收敛精度和最大迭代次数；将电力系统基本技术数据、电力系统运行约束基础数据、电力系统实时运行数据、非可再生能源发电企业报价数据、可再生能源发电企业报价数据、日前场景各非可再生能源发电机组的有功出力结果以及收敛准则参数输入到面向实时调节费用修正的最优直流潮流模型。
对计及可再生能源实时出力偏移量的实时场景进行最优直流潮流分析求解，求解t时刻非可再生能源发电机组Ci的向上有功调节量
Figure BDA00014646632900000917
和向下有功调节量
Figure BDA00014646632900000918
以及可再生能源发电机组Ri的实时正偏移有功出力
Figure BDA00014646632900000919
第7步：获取可再生能源发电机组的日前预测出力和节点电价结果，得到可再生能源交易的日前费用如可再生能源机组Ri的日前费用
Figure BDA00014646632900000920
获取实时场景下可再生能源出力的偏移量，以及非可再生能源发电机组的有功调节量结果，得到可再生能源交易的实时调节费用如可再生能源机组Ri的调节费用
Figure BDA0001464663290000101
Figure BDA0001464663290000102
其中
Figure BDA0001464663290000103
为可再生能源发电机组Ri于t时刻的系统调节费用的分摊系数，
Figure BDA0001464663290000104
由公式
Figure BDA0001464663290000105
获得；
可再生能源交易的日前费用以及实时调节费用之和即为可再生能源交易的最终结算费用如可再生能源机组Ri的最终结算费用
Figure BDA0001464663290000106
第8步：将得到的可再生能源交易费用反馈至可再生能源发电企业，可再生能源发电企业按照上述求解结果获得交易费用；将计及可再生能源出力实时变化量的实时场景最优直流潮流结果反馈给电网企业，电网企业参考求解结果进行全系统实时运行方式调整。

Claims (7)

1.一种基于节点电价机制的可再生能源交易费用的计算方法，其特征在于，包括以下步骤：
1)获取电力系统基本技术数据、电力系统运行约束基础数据、电力系统实时运行数据、非可再生能源发电企业报价数据以及可再生能源发电企业报价数据；
电力系统基本技术数据为节点数据、输电线路数据、变压器数据和负荷数据；
电力系统运行约束基础数据为各非可再生能源发电机组有功出力的上限值和下限值、输电线路的传输容量限值、可再生能源发电机组日前预测出力值；
电力系统实时运行数据为负荷实时偏差值、可再生能源发电机组实时预测出力值、各非可再生能源发电机组可提供的实时有功调节能力的上限值和下限值；
非可再生能源发电企业报价数据为非可再生能源发电企业对其所属非可再生能源发电机组出力费用的日前报价值，以及非可再生能源发电机组参与系统实时调节的实时报价值；
可再生能源发电企业报价数据为可再生发电企业对其所属可再生能源机组出力费用的日前报价值，以及可再生能源出力实时正偏移量费用的实时报价值；
2)以电力系统发电成本最小为目标函数构建最优直流潮流优化目标；然后构建约束条件，完成面向日前节点电价结算的最优直流潮流模型的构建；其中，约束条件包括系统功率平衡约束、直流潮流约束以及非可再生能源发电机组出力约束；
3)设定面向日前节点电价结算的最优直流潮流模型中的收敛准则参数，收敛准则参数包括优化模型收敛精度和最大迭代次数；将电力系统基本技术数据、电力系统运行约束基础数据、非可再生能源发电企业报价数据、可再生能源发电企业报价数据以及收敛准则参数输入到面向日前节点电价结算的最优直流潮流模型，对计及可再生能源日前预测出力的日前场景进行最优直流潮流分析求解，求解结果包括各节点的节点电价和各非可再生能源发电机组的有功出力；
4)以电力系统综合运行成本最小为目标函数构建最优直流潮流优化目标；然后构建约束条件，完成面向实时调节费用修正的最优直流潮流模型的构建；其中，约束条件包括系统功率平衡约束、直流潮流约束、可再生能源发电机组出力约束，非可再生能源发电机组调节能力约束；
5)设定面向实时调节费用修正的最优直流潮流模型中的收敛准则参数，收敛准则参数包括优化模型收敛精度和最大迭代次数；将电力系统基本技术数据、电力系统运行约束基础数据、电力系统实时运行数据、非可再生能源发电企业报价数据、可再生能源发电企业报价数据、日前场景各非可再生能源发电机组的有功出力结果以及收敛准则参数输入到面向实时调节费用修正的最优直流潮流模型，对计及可再生能源实时出力偏移量的实时场景进行最优直流潮流分析求解，求解结果包括各非可再生能源发电机组的有功调节量、各可再生能源发电机组有功出力的实时偏移量；
6)根据步骤1)得到的日前场景下的可再生能源发电机组的日前预测出力和步骤3)得到的节点电价结果，得到可再生能源交易的日前费用；获取实时场景下可再生能源出力的偏移量，以及非可再生能源发电机组的有功调节结果，得到可再生能源交易的实时调节费用；可再生能源交易的日前费用以及实时调节费用之和即为可再生能源交易的最终结算费用；
7)将步骤6)得到的可再生能源交易的最终结算费用反馈至可再生能源发电企业，可再生能源发电企业按照上述求解结果获得交易费用；将计及可再生能源出力实时变化量的实时场景最优直流潮流结果反馈给电网企业，电网企业参考最优直流潮流结果进行全系统实时运行方式调整。
2.根据权利要求1所述的基于节点电价机制的可再生能源交易费用计算方法，其特征在于：所述的步骤2)中最优直流潮流优化目标是由公式
Figure FDA0003169495930000021
Figure FDA0003169495930000022
获得的，其中，
Figure FDA0003169495930000023
为非可再生能源发电机组Ci在t时刻的有功出力，
Figure FDA0003169495930000024
为非可再生能源发电机组Ci在t时刻的报价参数，SGC为所有非可再生能源发电机组集合，
Figure FDA0003169495930000025
为可再生能源发电机组Ri在t时刻的日前预测有功出力，
Figure FDA0003169495930000026
为可再生能源发电机组Ri在t时刻的报价参数，SGR为所有可再生能源发电机组集合，ST为所有运行时段的集合。
3.根据权利要求1所述的基于节点电价机制的可再生能源交易费用计算方法，其特征在于，所述的步骤2)中系统功率平衡约束为：
Figure FDA0003169495930000027
其中，
Figure FDA0003169495930000028
为t时刻的负荷节点Di的有功负荷，SD为所有负荷节点的集合，SGC为所有非可再生能源发电机组集合，SGR为所有可再生能源发电机组集合，ST为所有运行时段的集合；
直流潮流约束为
Figure FDA0003169495930000031
其中，
Figure FDA0003169495930000032
为位于节点Bi的非可再生能源发电机组Ci在t时刻的有功出力，
Figure FDA0003169495930000033
为位于节点Bi的可再生能源发电机组Ri在t时刻的日前预测有功出力，
Figure FDA0003169495930000034
为位于节点Bi的t时刻的有功负荷；
Figure FDA0003169495930000035
为功率传输分布因子矩阵第l行第Bi列元素，
Figure FDA0003169495930000036
为支路l的传输功率上限值，Fl 为支路l的传输功率下限值，SL为所有支路的集合；
非可再生能源发电机组出力约束为
Figure FDA0003169495930000037
其中
Figure FDA0003169495930000038
为非可再生能源发电机组Ci的有功出力上限值，
Figure FDA0003169495930000039
为非可再生能源发电机组Ci的有功出力下限值。
4.根据权利要求1所述的基于节点电价机制的可再生能源交易费用计算方法，其特征在于，所述的步骤4)中最优直流潮流优化目标如下公式所示：
Figure FDA00031694959300000310
其中，
Figure FDA00031694959300000311
为可再生能源发电机组Ri在t时刻实时有功出力的负偏移值，当获得可再生能源发电机组Ri在t时刻的实时预测出力
Figure FDA00031694959300000312
时，负偏移值
Figure FDA00031694959300000313
由
Figure FDA00031694959300000314
获得；
Figure FDA00031694959300000315
为可再生能源发电机组Ri在t时刻的实时正偏移有功出力，
Figure FDA00031694959300000316
为可再生能源发电机组Ri在t时刻对其实时正偏移有功出力的报价参数；
Figure FDA00031694959300000317
为日前场景下非可再生能源发电机组Ci在t时刻的有功出力结果，
Figure FDA00031694959300000318
和
Figure FDA00031694959300000319
分别为非可再生能源发电机组Ci在t时刻的负有功调节出力和正有功调节出力，
Figure FDA00031694959300000320
和
Figure FDA00031694959300000321
分别为非可再生能源发电机组Ci在t时刻对其负有功调节出力和正有功调节出力的报价参数，SGC为所有非可再生能源发电机组集合，SGR为所有可再生能源发电机组集合，ST为所有运行时段的集合，
Figure FDA00031694959300000322
为非可再生能源发电机组Ci在t时刻的报价参数，
Figure FDA00031694959300000323
为可再生能源发电机组Ri在t时刻的日前预测有功出力，
Figure FDA00031694959300000324
为可再生能源发电机组Ri在t时刻的报价参数。
5.根据权利要求1所述的基于节点电价机制的可再生能源交易费用计算方法，其特征在于，所述的步骤4)中系统功率平衡约束为：
Figure FDA0003169495930000041
Figure FDA0003169495930000042
其中，
Figure FDA0003169495930000043
为t时刻负荷节点Di的有功负荷实时偏差值，
Figure FDA0003169495930000044
为日前场景下非可再生能源发电机组Ci在t时刻的有功出力结果，
Figure FDA0003169495930000045
为可再生能源发电机组Ri在t时刻的日前预测有功出力，
Figure FDA0003169495930000046
为t时刻的负荷节点Di的有功负荷，SD为所有负荷节点的集合，SGC为所有非可再生能源发电机组集合，SGR为所有可再生能源发电机组集合，ST为所有运行时段的集合；
直流潮流约束为:
Figure FDA0003169495930000047
其中，
Figure FDA0003169495930000048
为非可再生能源发电机组Ci在t时刻的有功出力，
Figure FDA0003169495930000049
和
Figure FDA00031694959300000410
分别为位于节点Bi的非可再生能源发电机组Ci在t时刻的负有功调节出力和正有功调节出力；
Figure FDA00031694959300000411
为位于节点Bi的可再生能源发电机组Ri在t时刻实时有功出力的负偏移值，
Figure FDA00031694959300000412
为位于节点Bi的可再生能源发电机组Ri在t时刻的实时正偏移有功出力；
Figure FDA00031694959300000413
为t时刻位于节点Bi的有功负荷的实时偏差值；
可再生能源发电机组出力约束为
Figure FDA00031694959300000414
其中
Figure FDA00031694959300000415
为可再生能源发电机组Ri在t时刻实时有功出力的正偏移值，当获得可再生能源发电机组Ri在t时刻的实时预测出力
Figure FDA00031694959300000416
时，正偏移值
Figure FDA00031694959300000417
由
Figure FDA00031694959300000418
获得；
非可再生能源发电机组调节能力约束为
Figure FDA00031694959300000419
其中
Figure FDA00031694959300000420
和
Figure FDA00031694959300000421
分别为非可再生能源发电机组Ci在t时刻正有功调节出力和负有功调节出力的上限值。
6.根据权利要求1所述的基于节点电价机制的可再生能源交易费用计算方法，其特征在于，所述的步骤6)中对应可再生能源发电机组Ri的日前费用
Figure FDA00031694959300000422
为：
Figure FDA00031694959300000423
其中，ST为所有运行时段的集合，
Figure FDA00031694959300000424
为可再生能源发电机组Ri在t时刻的日前预测有功出力，
Figure FDA00031694959300000425
为节点Bi在t时刻的节点电价，单位为元/kWh，节点电价
Figure FDA00031694959300000426
由公式
Figure FDA00031694959300000427
获得，μt为对应于t时刻系统平衡约束
Figure FDA00031694959300000428
的对偶乘子取值，λlt为对应于t时刻支路l直流潮流约束
Figure FDA0003169495930000051
的对偶乘子取值，SL为所有支路的集合，
Figure FDA0003169495930000052
为功率传输分布因子矩阵第l行第Bi列元素；
对应可再生能源发电机组Ri的实时调节费用
Figure FDA0003169495930000053
为：
Figure FDA0003169495930000054
其中ST为所有运行时段的集合，SGC为所有非可再生能源发电机组集合，SGR为所有可再生能源发电机组集合，
Figure FDA0003169495930000055
为可再生能源发电机组Ri在t时刻的实时正偏移有功出力，
Figure FDA0003169495930000056
为可再生能源发电机组Ri在t时刻实时有功出力的负偏移值，
Figure FDA0003169495930000057
为节点Bi在t时刻的节点电价，
Figure FDA0003169495930000058
为可再生能源发电机组Ri于t时刻的系统调节费用的分摊系数，
Figure FDA0003169495930000059
由公式
Figure FDA00031694959300000510
获得SD为所有负荷节点的集合，
Figure FDA00031694959300000511
和
Figure FDA00031694959300000512
分别为非可再生能源发电机组Ci在t时刻的负有功调节出力和正有功调节出力，
Figure FDA00031694959300000513
和
Figure FDA00031694959300000514
分别为非可再生能源发电机组Ci在t时刻对其负有功调节出力和正有功调节出力的报价参数，
Figure FDA00031694959300000515
为t时刻位于节点Bi的有功负荷的实时偏差值；
对应可再生能源发电机组Ri的交易费用
Figure FDA00031694959300000516
为：
Figure FDA00031694959300000517
7.根据权利要求1所述的基于节点电价机制的可再生能源交易费用计算方法，其特征在于，所述非可再生能源发电为火力发电。
CN201711107239.9A 2017-11-10 2017-11-10 一种基于节点电价机制的可再生能源交易费用的计算方法 Active CN107833101B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711107239.9A CN107833101B (zh) 2017-11-10 2017-11-10 一种基于节点电价机制的可再生能源交易费用的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711107239.9A CN107833101B (zh) 2017-11-10 2017-11-10 一种基于节点电价机制的可再生能源交易费用的计算方法

Publications (2)

Publication Number Publication Date
CN107833101A CN107833101A (zh) 2018-03-23
CN107833101B true CN107833101B (zh) 2021-08-27

Family

ID=61654123

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711107239.9A Active CN107833101B (zh) 2017-11-10 2017-11-10 一种基于节点电价机制的可再生能源交易费用的计算方法

Country Status (1)

Country Link
CN (1) CN107833101B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109102155B (zh) * 2018-07-09 2022-03-15 中国南方电网有限责任公司 一种超短期节点边际电价概率预测方法及系统
CN109523309B (zh) * 2018-10-26 2021-04-30 武汉华中思能科技有限公司 一种基于潮流分析计算的电力报价仿真方法及系统
CN110232583B (zh) * 2019-01-31 2021-01-05 广东电力交易中心有限责任公司 一种考虑碳排放权的电力市场边际价格规划方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007028720A (ja) * 2005-07-13 2007-02-01 Hitachi Ltd 電力供給調整時損失予測システム、電力供給調整時損失予測方法、電力供給調整時損失予測プログラム
CN104036414A (zh) * 2014-06-19 2014-09-10 西安交通大学 一种基于最优潮流的转运交易费用的计算方法
CN104156581A (zh) * 2014-08-01 2014-11-19 重庆大学 计及系统可靠性与购电风险的峰谷分时电价确定方法
CN106339794A (zh) * 2016-08-16 2017-01-18 清华大学 一种电‑热耦合多能流网络节点能价计算方法
CN106886842A (zh) * 2017-03-22 2017-06-23 华北电力大学 一种分布式供电系统的部署方法、装置和计算设备
CN107301537A (zh) * 2017-06-16 2017-10-27 昆明电力交易中心有限责任公司 一种区域电力市场购电主体结算系统和方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007028720A (ja) * 2005-07-13 2007-02-01 Hitachi Ltd 電力供給調整時損失予測システム、電力供給調整時損失予測方法、電力供給調整時損失予測プログラム
CN104036414A (zh) * 2014-06-19 2014-09-10 西安交通大学 一种基于最优潮流的转运交易费用的计算方法
CN104156581A (zh) * 2014-08-01 2014-11-19 重庆大学 计及系统可靠性与购电风险的峰谷分时电价确定方法
CN106339794A (zh) * 2016-08-16 2017-01-18 清华大学 一种电‑热耦合多能流网络节点能价计算方法
CN106886842A (zh) * 2017-03-22 2017-06-23 华北电力大学 一种分布式供电系统的部署方法、装置和计算设备
CN107301537A (zh) * 2017-06-16 2017-10-27 昆明电力交易中心有限责任公司 一种区域电力市场购电主体结算系统和方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
An Extension of Reduced Disjunctive Model for Multi-Stage Security-Constrained Transmission Expansion Planning;Yao Zhang等;《IEEE Transactions on Power Systems》;20170523;全文 *
Flexible transmission network expansion planning under uncertainty based on a self-adaptive clustering technique;Yunhao Li等;《International Conference on Renewable Power Generation》;20160407;全文 *
区域电力市场中的分区备用模型;王建学等;《中国电机工程学报》;20060930;全文 *
电力市场中区域电价综合模型;丁晓莺等;《电力系统自动化》;20040525;全文 *

Also Published As

Publication number Publication date
CN107833101A (zh) 2018-03-23

Similar Documents

Publication Publication Date Title
CN107833101B (zh) 一种基于节点电价机制的可再生能源交易费用的计算方法
CN101777764B (zh) 基于区域控制偏差和机组发电偏差的调度反馈控制系统
CN104951899A (zh) 含大规模可再生能源的配电公司多时间尺度优化调度方法
CN101281637A (zh) 基于超平面形式安全域边界的电力系统优化潮流和实时定价方法
CN101431240A (zh) 一种计及梯级电站的省地一体化avc系统控制方法
CN112529304B (zh) 一种基于配额制考虑风险的两级电力市场优化运行方法
Hou et al. A review on China's current situation and prospects of poverty alleviation with photovoltaic power generation
CN102593839A (zh) 计及电网全运行方式的发电机励磁系统调差系数整定方法
CN103927590B (zh) 信息不对称情况下的电力系统优化控制及费用计算方法
CN107492886A (zh) 一种区域电力市场下含风电的电网月度购电方案优化方法
CN105515010A (zh) 一种基于合作对策的二级电压协调控制方法及系统
CN116436039A (zh) 一种虚拟电厂辅助电网二次调频优化调度方法及装置
CN115358519A (zh) 一种虚拟电厂优化调度方法及装置
CN115481856A (zh) 考虑综合需求响应的综合能源系统多尺度调度方法及系统
CN105610200A (zh) 一种基于同步协调控制的火电厂全功率控制方法
Tang et al. Unit maintenance strategy considering the uncertainty of energy intensive load and wind power under the carbon peak and carbon neutral target
CN108197788A (zh) 一种对等控制模式下微电网电压频率偏差估计方法
Zhang et al. The economic impact of energy storage co-deployment on renewable energy in China
CN112531703A (zh) 一种多能源虚拟电厂提供多市场与本地服务的优化方法
CN108108871A (zh) 风电场群输电设备的选型方法
CN110826778A (zh) 一种主动适应新能源发展的负荷特性优化计算方法
CN116683529A (zh) 风光储工业园区优化运行的实时控制方法、系统及设备
Liu et al. An Overview of ancillary service development in China: status, problems and countermeasures
CN112508365B (zh) 一种在线煤耗曲线滚动修正方法及系统
CN109004689A (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