CN114862457A - 梯级水电站群参与日前现货市场的竞价电量边界确定方法 - Google Patents

梯级水电站群参与日前现货市场的竞价电量边界确定方法 Download PDF

Info

Publication number
CN114862457A
CN114862457A CN202210446037.1A CN202210446037A CN114862457A CN 114862457 A CN114862457 A CN 114862457A CN 202210446037 A CN202210446037 A CN 202210446037A CN 114862457 A CN114862457 A CN 114862457A
Authority
CN
China
Prior art keywords
hydropower station
water
water level
constraint
hydropower
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
CN202210446037.1A
Other languages
English (en)
Other versions
CN114862457B (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.)
Dalian University of Technology
China Southern Power Grid Co Ltd
Original Assignee
Dalian University of Technology
China Southern Power Grid 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 Dalian University of Technology, China Southern Power Grid Co Ltd filed Critical Dalian University of Technology
Priority to CN202210446037.1A priority Critical patent/CN114862457B/zh
Publication of CN114862457A publication Critical patent/CN114862457A/zh
Application granted granted Critical
Publication of CN114862457B publication Critical patent/CN114862457B/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/02—Marketing; Price estimation or determination; Fundraising
    • G06Q30/0283—Price estimation or determination
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/23—Clustering techniques
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • 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
    • 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/0202—Market predictions or forecasting for commercial activities
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/04—Constraint-based CAD
    • 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/14—Marketing, i.e. market research and analysis, surveying, promotions, advertising, buyer profiling, customer management or rewards

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Development Economics (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • Finance (AREA)
  • Accounting & Taxation (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Data Mining & Analysis (AREA)
  • Economics (AREA)
  • Marketing (AREA)
  • General Business, Economics & Management (AREA)
  • Evolutionary Computation (AREA)
  • Game Theory and Decision Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Primary Health Care (AREA)
  • Tourism & Hospitality (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Public Health (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Biology (AREA)
  • Water Supply & Treatment (AREA)
  • Medical Informatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Software Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Geometry (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及电力市场、水电调度及电网运行领域，公开了一种梯级水电站群参与日前现货市场的竞价电量边界确定方法。随着我国电力现货市场建设不断完善，水电参与现货市场势在必行。在电力市场中，出清效率关乎电网安全平稳运行。本发明针对梯级上下游水库复杂水力联系导致现货出清模型复杂、出清效率低，且容易导致水电参与现货市场出现弃水和交易结果难以执行等问题，进行梯级水电站日前市场化边界条件确定，提出了水电站绝对可行电量区间的概念，揭示了梯级水电站电量边界随自身可发电水量变化的演变规律，完全规避了上下游电站不匹配的问题，提高了出清模型求解效率，为市场运营者管理水电参与现货市场提供一种切实可行的思路。

Description

梯级水电站群参与日前现货市场的竞价电量边界确定方法
技术领域
本发明涉及电力市场、水电调度及电网运行领域，特别涉及一种梯级水电站群参与日前现货市场的竞价电量边界确定方法。
背景技术
随着我国电力现货市场建设不断完善，水电参与现货市场势在必行。不同于火电等其他电源，水电参与现货市场主要面临以下难题：
1)水电现货交易中上下游竞价匹配困难：即由于梯级上下游电站可能归属不同主体，水力联系却又紧密耦合，造成上游电站竞争电量多、下游电站竞争电量少导致下游电站弃水，或是上游电站竞争电量少、下游电站竞争电量多导致下游无水可发的情况。
2)复杂水力约束降低出清效率：在电力市场中，出清效率关乎电网安全及其平稳运行。目前国内外的电力市场现货出清，大多是通过混合整数线性规划(MILP)进行出清，然后通过商用求解器求解，现货出清效率跟约束数密切相关。水电中存在水位-库容、机组NHQ曲线等众多非凸、非线性约束，通常采用分段线性化方法将其加入到现货出清约束条件中，这就造成MILP出清模型约束数急剧增加，导致出清效率大幅下降，难以满足现货出清的应用要求。
在这种情况下，对于水电居主导地位或者水电占比较大的电网，如何简化梯级水电站群参与日前现货市场出清，将上下游之间的水力联系转化为电力联系，确定其竞价电量边界值得深入研究。
本发明成果依托于国家自然科学基金(91547201)，以梯级水电参与电力现货市场为研究背景，以澜沧江梯级水电站群为研究对象，结合南方电网现货市场建设，提出了一种梯级水电站群参与日前现货市场的竞价电量边界确定方法。本发明方法可以解耦梯级水电之间的水力联系，规避上下游电站竞价不匹配等问题，提高水电参与现货出清的效率，为市场运营者管理水电参与现货市场提供一种切实可行的思路，具有很强实用性和广泛推广价值。
发明内容
本发明要解决的技术问题是提供一种梯级水电站群参与日前现货市场的竞价电量边界确定方法，以期提高水电参与现货出清的效率，避免梯级水电出清因为上下游竞价不匹配导致的交易结果难执行等问题。
本发明的技术方案为：本发明揭示了一种梯级水电站群参与日前现货市场的竞价电量边界确定方法，具体包括如下步骤：
(1)提出梯级水电站群绝对可行区间概念
当上游水电站m平均出库流量为最大值
Figure BDA0003615587680000011
时，下游水电站n的发电水量变化区间为
Figure BDA0003615587680000012
其对应的发电电量变化区间为
Figure BDA0003615587680000013
式中：
Figure BDA0003615587680000014
Figure BDA0003615587680000015
当上游电站m平均出库流量为最小值Rm 时，下游水电站n的发电水量变化区间为
Figure BDA0003615587680000021
其对应的发电电量变化区间为
Figure BDA0003615587680000022
式中：
Figure BDA0003615587680000023
Figure BDA0003615587680000024
其中，
Figure BDA0003615587680000025
表示初始库容，单位m3、qn表示区间来流，单位m3/s、
Figure BDA0003615587680000026
表示末库容最小值，单位m3、
Figure BDA0003615587680000027
表示末库容最大值，单位m3、Δ为时段长，单位h。无论上游电站来水条件如何，
Figure BDA0003615587680000028
之间发电流量对应的交易电量
Figure BDA0003615587680000029
总是可行的，本发明称此区间为绝对可行区间。
(2)根据步骤(1)所提概念，建立梯级水电站群绝对可行电量区间最大模型
模型可描述为：已知日前市场中各水电站的初始水位、期末水位上下限以及预测来水，在满足水电站运行约束条件下，确定各水电站的电量上下限，使梯级水电站群在日前市场中的绝对电量可行区间最大。
模型的目标函数为：
Figure BDA00036155876800000210
Figure BDA00036155876800000211
E n,t＝W n,t/rn,t (7)
式中：
Figure BDA00036155876800000212
分别为水电站n在t时段的发电量上限、发电水量上限，单位分别为kW·h、m3；E n,t、W n,t为水电站n在t时段的发电量下限、发电水量下限，单位分别为kW·h、m3；rn,t为水电站n在t时段的平均耗水率，m3/kW·h；T为总调度时段数，h，N为梯级水电站总个数，λ为惩罚因子，Sn,t为水电站n在t时段的弃水流量，m3/s。
模型的约束条件包括线性约束和非线性约束，其中，线性约束包括水量平衡约束、始末库水位约束、出库流量约束、发电引用流量约束、出力约束等，非线性约束包括水位库容、尾水位泄量、水头耗水率等。
(3)采取估算策略确定步骤(2)模型的始末库水位约束条件
上述模型中，始末库水位约束在确定竞价电量边界时尚未确定，故需采取一定的策略估算。
对于水库初始水位，采取“以电定水”方式推算，具体方式为从当前时刻电站水位开始，计算出剩余计划电量所需水量，根据水量平衡方程估算出下一日的日初水位。
对于梯级各水库日末水位范围约束，运用模糊聚类分析方法估算出日末理想水位，结合调度人员根据工程经验给出的电站理想水位的允许偏差值，得出水库日末水位范围。可以表示为下式：
Figure BDA00036155876800000213
式中：
Figure BDA00036155876800000214
ΔZ、
Figure BDA00036155876800000215
分别为第n个水电站的期望末水位、末水位下、上偏差值。
(4)将步骤(2)所提出的模型通过线性化处理转化为线性模型
采用分段线性插值方法将水位库容、尾水位泄量、水头耗水率等非线性关系转化为线性约束，进而将模型转化为线性模型。
(5)采用单纯形法求解步骤(4)得到的线性模型，得到梯级水电站群最大可行电量区间，即为梯级水电站群参与日前现货市场的竞价电量边界。
本发明对比现有技术有如下有益效果：本发明一种梯级水电站群参与日前现货市场的竞价电量边界确定方法，提出了绝对可行区间的概念，将水电站复杂的水力联系转化成电力联系，揭示了梯级各电站电量竞价区间随自身调蓄能力变化情况的演变规律，完全规避了电站竞价过程中上下游竞价不匹配的问题，提高了含水电的现货市场出清模型的求解速度，确保了水电现货出清结果的合理性和可行性，同时也就降低了水电参与现货市场的风险，从而提高梯级水电参与市场的积极性，为市场运营者管理水电参与现货市场提供一定的决策依据。
附图说明
图1是绝对可行区间示意图；
图2是本发明的总体求解框架示意图；
图3是现货出清流程图；
图4是梯级水电站绝对可行电量区间最大模型的求解流程图；
图5是实施例中澜沧江梯级拓扑结构示意图；
图6是实施例中汛期满发计算结果图。
具体实施方式
下面结合附图和实例对本发明作进一步的描述，主要包括梯级水电参与日前现货市场的竞价电量边界的建模、求解方法及实际应用三部分。
1.日前市场化边界条件的建模方法
(1)确定电量竞价区间—绝对可行区间
为简化梯级水电站上下游紧密的水力联系，本发明提出绝对可行区间的概念。具体分析如下：
当上游水电站m平均出库流量为最大值
Figure BDA0003615587680000031
时，下游水电站n的发电水量变化区间为
Figure BDA0003615587680000032
其对应的发电电量变化区间为
Figure BDA0003615587680000033
式中：
Figure BDA0003615587680000034
Figure BDA0003615587680000035
当上游电站m平均出库流量为最小值Rm 时，下游水电站n的发电水量变化区间为
Figure BDA0003615587680000036
其对应的发电电量变化区间为
Figure BDA0003615587680000037
式中：
Figure BDA0003615587680000038
Figure BDA0003615587680000039
其中，
Figure BDA00036155876800000310
表示初始库容，m3、qn表示区间来流，m3/s、
Figure BDA00036155876800000311
表示末库容最小值，m3、
Figure BDA00036155876800000312
表示末库容最大值，m3、Δ为时段长，h。无论上游电站来水条件如何，
Figure BDA00036155876800000313
之间发电流量对应的交易电量
Figure BDA0003615587680000041
总是可行的，这一区间称为绝对可行区间。
如图1所示，无论上游电站来水条件如何之间发电流量对应的交易电量总是可行的，这一区间称为绝对可行区间。为此本发明将其作为水电站日前市场的电量竞价区间。
(2)目标函数
梯级水电因为调节性能差异，在电力市场竞价过程中，竞价空间表现不同，而其空间又跟梯级上下游电站的调节性能差异密切相关。为此在保证梯级水电弃水最少的情况，满足清洁能源消纳准则，同时兼顾水电参与市场有一定的弹性空间，从而水电站追求电量竞价空间最大。因此目标函数选取梯级水电站绝对可行电量区间之和最大作为目标函数。具体可描述为：已知日前市场中各水电站的初始水位、期末水位上下限以及预测来水，在满足水电站运行约束条件下，确定各水电站的电量上下限，使梯级水电站群在日前市场中绝对可行电量区间最大。数学表达式为：
Figure BDA0003615587680000042
Figure BDA0003615587680000043
E n,t＝W n,t/rn,t (15)
式中：
Figure BDA0003615587680000044
分别为水电站n在t时段的发电量上限、发电水量上限，kW·h、m3；En,t、W n,t为水电站n在t时段的发电量下限、发电水量下限，kW·h、m3；rn,t为水电站n在t时段的平均耗水率，m3/kW·h；T为总调度时段数，h，N为梯级水电站总个数，λ为惩罚因子，Sn,t为水电站n在t时段的弃水流量，m3/s。
(3)约束条件
a)水量平衡方程
Vn,t＝Vn,t-1+3600×(In,t-Qn,t-Sn,t)Δt (16)
Figure BDA0003615587680000045
式中：Vn,t为第n个水电站第t时段的库容，单位m3；In,t、qn,t、Qn,t、Sn,t分别是第n个水电站第t时段的入库流量、区间流量、发电流量和弃水流量，单位m3/s；mn为第n个水电站的第m个直接上游水电站；Ωn为第n个水电站的直接上游水电站集合，对于龙头水电站该集合为空集；Δt为计算时段。
b)始末库水位约束
Figure BDA0003615587680000046
Figure BDA0003615587680000047
式中：
Figure BDA0003615587680000048
ΔZ、
Figure BDA0003615587680000049
分别为第n个水电站的初始水位、期望末水位、末水位下、上偏差值。
c)发电引用流量约束
Figure BDA00036155876800000410
式中：Q n,t、
Figure BDA0003615587680000051
分别是第n个水电站第t时段的发电引用流量下限、上限，单位m3/s。
d)出库流量约束
Figure BDA0003615587680000052
式中：
Figure BDA0003615587680000053
分别是第n个水电站第t时段的出库流量下限、上限，单位m3/s。
e)出力约束
Figure BDA0003615587680000054
式中：N n,t、
Figure BDA0003615587680000055
分别是第n个水电站第t时段的出力下限、上限，单位kW。
f)水位库容约束
Zn,t＝fn,ZV(Vn,t) (23)
式中：fn,ZV(·)为水电站n的水位-库容关系函数。
g)尾水位泄量约束
Zn,t＝fn,ZQ(Qn,t) (24)
式中：fn,ZQ(·)为水电站n的尾水位-泄量关系函数。
h)水头耗水率约束
rn,t＝fn,rh(hn,t) (25)
式中：fn,rh(·)为水电站n的水头-耗水率关系函数。
2.日前市场化边界条件模型求解方法
(1)总体求解思路
梯级水电站群绝对可行电量区间优化问题包含复杂的水力非线性条件约束，例如水位库容曲线、尾水位泄量曲线。该问题是一个多维、动态、强耦合、非线性规划问题，直接求解非常困难。为保证模型求解效率及工程应用的可行性，本发明将非线性约束转化为线性约束，将非线性规划问题转化为线性规划问题，从而可以利用经典的单纯形法进行精确求解。总体求解思路见图2。
(2)非线性约束线性化策略
采用分段线性插值方法处理非线性约束。以水位库容约束为例说明。根据每一个水电站的调度期内库容变动范围
Figure BDA0003615587680000056
在区间范围内选取k-2个节点，将第n个水电站的库容离散成k个区间，如式(26)—(30)所示。
Figure BDA0003615587680000057
Zn,k＝f(Vn,k) (27)
αn,t,kVn,k-1≤Vn,t,k≤αn,t,kVn,k-1 (28)
Figure BDA0003615587680000058
式中：Vn,k、Zn,k分别为第n个水电站的第k个插值节点的库容和水位；αn,t,k为0-1变量，用于表示库容所在区间范围；
Figure BDA0003615587680000061
确保第n水电站第t时段末库容落在唯一区间。
综上，任意电站任意时段末的库容可由(30)式求解。
Figure BDA0003615587680000062
其余非线性约束处理策略与上述类似，此处不在赘述。
(3)电站日初及日末水位范围以及耗水率估算策略
实际现货市场出清流程如图3所示，其中边界条件准备位于第一步骤，即需要在运行日前一天没有结束前必须计算出下一日的梯级水电电量边界条件，然而此时下一日日初水位并未确定，为了实际工程应用需求，本发明采取以电定水的方式推算下一日日初水位。具体方式为从当前时刻电站水位开始计算，算出剩余计划电量所需水量，即可估算出下一日的日初水位。
关于日末水位，由水量平衡方程公式(16)可知，梯级每座电站日末水位范围主要取决于入库流量和出库流量最大值和最小值，其中入库流量为上游出库流量和区间流量之和。因此可以采用自身及上游水库的初始水位、区间流量、历史出力的最大值和最小值作为典型特征物理量，运用模糊聚类分析方法可以估算出日末理想水位值，然后结合调度人员的工程经验给出电站允许偏差理想水位的上偏值和下偏值，从而得出水库的日末水位范围。在此基础上，计算平均出库流量(泄量)，然后根据尾水位泄量曲线线性插值得到尾水位，接着计算发电水头，最后由水头耗水率曲线分段线性插值确定每座电站的平均耗水率。
(4)求解流程
结合上述总体求解思路、约束条件处理策略以及线性规划模型求解方法，给出梯级水电站群发电量区间优化问题的详细求解过程：
a)确定基础数据，包括上下游水力联系、水位库容关系、尾水位泄量关系、水头耗水率关系等信息，构建梯级水电站绝对可行电量区间最大模型。
b)将初始条件输入a)模型中，包括梯级每座水电站的调度期初始水位，天然来水流量。
c)采用(2)所提的分段线性插值方法对非线性约束条件线性化处理，将a)模型转化成线性规划模型。
d)采用(3)策略确定梯级每座电站的日末水位上下限以及平均耗水率。
e)设置目标函数和约束条件，构建绝对可行电量区间最大模型。具体约束条件包括水量平衡约束、库水位约束、出库流量约束、发电引用流量约束、出力约束、末水位上下限约束等。
f)将目标函数及约束条件转化成线性规划模型标准形式，引入单纯形法求出线性规划模型的初始基本可行解，及一组可行的出库流量和发电流量。
g)判断目标函数是否达到最优。
h)若目标函数达到最优，则计算结束；若不是最优解，则进行换基迭代，直至目标函数最优为止。
梯级水电站群绝对可行电量区间最大模型求解流程图如图4所示。
3.实例分析
(1)模型求解结果
现以某日澜沧江梯级水电站群发电调度为例，验证所提方法的有效性。澜沧江梯级投产水电站11座，涵盖日、季、年等不同调节类型水库，梯级电站拓扑结构见图5，水电站主要参数见表1。该日梯级各电站的预测来水、水库水位上下限均来源实际调度场景。优化计算结果见表2，各电站绝对可行日电量竞价区间见图6。
表1水电站主要参数
Figure BDA0003615587680000071
表2优化计算结果
Figure BDA0003615587680000072
根据计算结果，该日澜沧江梯级水电群日前市场电量竞价区间带宽为5799万kW·h。电站A、B为日调节水库，调节库容有限，且处于高水位运行，由于日天然来水已明显超过满发流量，水库蓄满后仍有多余水量，为避免产生弃水电量电站为满发运行，电量上下限重合，竞价区间带宽为0是合理的。电站C为季调节水库，该日仍有一定调节库容，为避免下游日调节电站D弃水，充分发挥调蓄作用，将水位抬升。电站G、J具有年调节能力，是梯级主力发电水库，承担了主要调节任务，日前电量竞价区间带宽分别为1969、1430万kW·h，分别占梯级总电量竞价区间带宽的33.95％、24.66％。
(2)方法效果检验
将通过上述方法确定的竞价电量边界代入含水电参与现货市场出清模型检验，具体以云南省参与省调平衡的122座水电站和11座火电站进行出清计算，并重点考虑与不简化水力约束的出清时间对比结果，以及A-K电站中标电量和弃水对比情况。此处出清模型采用“发电侧报量报价，用电侧不报量不报价”日前现货出清方式，以系统总购电费用最小为目标函数。结果如下表3、4所示：
表3现货清算模型计算时间对照表
Figure BDA0003615587680000081
由上表3可以看出电量边界条件的加入，使得出清算法仅需迭代1次，平均出清计算时间从165秒减少至85秒，计算效率提高了48.5％，大大提高实际应用效率。
表4电站中标电量及弃水情况对比结果
Figure BDA0003615587680000082
由上表4可以看出，当不简化梯级电站的水力约束时，由于梯级电站隶属于不同的利益主体，互相之间缺乏参与现货市场的信息交互，造成发电侧电站不合理的报量报价行为，这极容易导致上游电站中标电量多、下游电站中标电量少，使得下游C，D，F出现严重的弃水现象。这样的出清结果对于水电运行是不合理的，易造成清洁能源严重的浪费。加上梯级电站的电量竞价边界条件之后，规范了电站报量的合理范围，A电站的中标电量从989万减少至306万、B电站中标电量从435万减少至88万，抬高了水电站运行水位，充分利用水电站的蓄水能力，使下游电站均变成无弃水状态，有效解决了上下游电站中标电量不匹配而导致的弃水问题。

Claims (2)

1.一种梯级水电站群参与日前现货市场的竞价电量边界确定方法，其特征在于，具体步骤如下：
(1)提出梯级水电站群绝对可行区间概念
当上游水电站m平均出库流量为最大值
Figure FDA0003615587670000011
时，下游水电站n的发电水量变化区间为
Figure FDA0003615587670000012
其对应的发电电量变化区间为
Figure FDA0003615587670000013
式中：
Figure FDA0003615587670000014
Figure FDA0003615587670000015
当上游电站m平均出库流量为最小值Rm 时，下游水电站n的发电水量变化区间为
Figure FDA0003615587670000016
其对应的发电电量变化区间为
Figure FDA0003615587670000017
式中：
Figure FDA0003615587670000018
Figure FDA0003615587670000019
其中，
Figure FDA00036155876700000110
表示初始库容，单位m3、qn表示区间来流，单位m3/s、
Figure FDA00036155876700000111
表示末库容最小值，单位m3、
Figure FDA00036155876700000112
表示末库容最大值，单位m3、Δ为时段长，单位h；无论上游电站来水条件如何，
Figure FDA00036155876700000113
之间发电流量对应的交易电量
Figure FDA00036155876700000114
总是可行的，称此区间为绝对可行区间；
(2)根据步骤(1)所提概念，建立梯级水电站群绝对可行电量区间最大模型
模型描述为：已知日前市场中各水电站的初始水位、期末水位上下限以及预测来水，在满足水电站运行约束条件下，确定各水电站的电量上下限，使梯级水电站群在日前市场中的绝对电量可行区间最大；
模型的目标函数为：
Figure FDA00036155876700000115
Figure FDA00036155876700000116
E n,t＝W n,t/rn,t (7)
式中：
Figure FDA00036155876700000117
分别为水电站n在t时段的发电量上限、发电水量上限，单位分别为kW·h、m3；E n,t、W n,t为水电站n在t时段的发电量下限、发电水量下限，单位分别为kW·h、m3；rn,t为水电站n在t时段的平均耗水率，m3/kW·h；T为总调度时段数，h，N为梯级水电站总个数，λ为惩罚因子，Sn,t为水电站n在t时段的弃水流量，m3/s；
模型的约束条件包括线性约束和非线性约束，其中，线性约束包括水量平衡约束、始末库水位约束、出库流量约束、发电引用流量约束和出力约束，非线性约束包括水位库容、尾水位泄量和水头耗水率；
(3)采取估算策略确定步骤(2)模型的始末库水位约束条件
对于水库初始水位，采取“以电定水”方式推算，具体方式为从当前时刻电站水位开始，计算出剩余计划电量所需水量，根据水量平衡方程估算出下一日的日初水位；
对于梯级各水库日末水位范围约束，运用模糊聚类分析方法估算出日末理想水位，结合调度人员根据工程经验给出的电站理想水位的允许偏差值，得出水库日末水位范围，表示为下式：
Figure FDA0003615587670000021
式中：
Figure FDA0003615587670000022
ΔZ、
Figure FDA0003615587670000023
分别为第n个水电站的期望末水位、末水位下、上偏差值；
(4)将步骤(2)所提出的模型通过线性化处理转化为线性模型
采用分段线性插值方法非线性关系转化为线性约束，进而将模型转化为线性模型；
(5)采用单纯形法求解步骤(4)得到的线性模型，得到梯级水电站群最大可行电量区间，即为梯级水电站群参与日前现货市场的竞价电量边界。
2.根据权利要求1所述的一种梯级水电站群参与日前现货市场的竞价电量边界确定方法，其特征在于，所述的梯级水电站绝对可行电量区间最大模型的约束条件包括：
a)水量平衡方程
Vn,t＝Vn,t-1+3600×(In,t-Qn,t-Sn,t)Δt (9)
Figure FDA0003615587670000024
式中：Vn,t为第n个水电站第t时段的库容，单位m3；In,t、qn,t、Qn,t、Sn,t分别是第n个水电站第t时段的入库流量、区间流量、发电流量和弃水流量，单位m3/s；mn为第n个水电站的第m个直接上游水电站；Ωn为第n个水电站的直接上游水电站集合，对于龙头水电站该集合为空集；Δt为计算时段；
b)始末库水位约束
Figure FDA0003615587670000025
Figure FDA0003615587670000026
式中：
Figure FDA0003615587670000027
ΔZ、
Figure FDA0003615587670000028
分别为第n个水电站的初始水位、期望末水位、末水位下、上偏差值；
c)发电引用流量约束
Figure FDA0003615587670000029
式中：Q n,t、
Figure FDA00036155876700000210
分别是第n个水电站第t时段的发电引用流量下限、上限，单位m3/s；
d)出库流量约束
Figure FDA00036155876700000211
式中：
Figure FDA00036155876700000212
分别是第n个水电站第t时段的出库流量下限、上限，单位m3/s；
e)出力约束
Figure FDA00036155876700000213
式中：N n,t、
Figure FDA00036155876700000214
分别是第n个水电站第t时段的出力下限、上限，单位kW；
f)水位库容约束
Zn,t＝fn,ZV(Vn,t) (16)
式中：fn,ZV(·)为水电站n的水位-库容关系函数；
g)尾水位泄量约束
Zn,t＝fn,ZQ(Qn,t) (17)
式中：fn,ZQ(·)为水电站n的尾水位-泄量关系函数；
h)水头耗水率约束
rn,t＝fn,rh(hn,t) (18)
式中：fn,rh(·)为水电站n的水头-耗水率关系函数。
CN202210446037.1A 2022-04-26 2022-04-26 梯级水电站群参与日前现货市场的竞价电量边界确定方法 Active CN114862457B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210446037.1A CN114862457B (zh) 2022-04-26 2022-04-26 梯级水电站群参与日前现货市场的竞价电量边界确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210446037.1A CN114862457B (zh) 2022-04-26 2022-04-26 梯级水电站群参与日前现货市场的竞价电量边界确定方法

Publications (2)

Publication Number Publication Date
CN114862457A true CN114862457A (zh) 2022-08-05
CN114862457B CN114862457B (zh) 2024-08-13

Family

ID=82633460

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210446037.1A Active CN114862457B (zh) 2022-04-26 2022-04-26 梯级水电站群参与日前现货市场的竞价电量边界确定方法

Country Status (1)

Country Link
CN (1) CN114862457B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115545856A (zh) * 2022-10-06 2022-12-30 大连川禾绿能科技有限公司 一种计及电价不确定性的梯级水电日前市场组合竞价方法
CN115983733A (zh) * 2023-01-03 2023-04-18 中国南方电网有限责任公司 基于水位控制的电力市场出清数据处理方法和装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20200090285A1 (en) * 2018-03-16 2020-03-19 Dalian University Of Technology Method for short-term generation scheduling of cascade hydropower plants coupling cluster analysis and decision tree
CN111160640A (zh) * 2019-12-23 2020-05-15 大连理工大学 水电站耗水率迭代修正的电力系统日前现货市场出清方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20200090285A1 (en) * 2018-03-16 2020-03-19 Dalian University Of Technology Method for short-term generation scheduling of cascade hydropower plants coupling cluster analysis and decision tree
CN111160640A (zh) * 2019-12-23 2020-05-15 大连理工大学 水电站耗水率迭代修正的电力系统日前现货市场出清方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张粒子;刘方;许通;蔡华祥;蒋燕;徐宏;: "多运营主体梯级水电站参与的日前市场出清模型", 电力系统自动化, no. 16, 17 July 2018 (2018-07-17) *
贾泽斌等: "基于梯级电量联动控制的高比例水电电力市场日前现货出清方法", 中国电机工程学报, vol. 41, no. 23, 5 December 2021 (2021-12-05) *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115545856A (zh) * 2022-10-06 2022-12-30 大连川禾绿能科技有限公司 一种计及电价不确定性的梯级水电日前市场组合竞价方法
CN115545856B (zh) * 2022-10-06 2023-06-16 大连川禾绿能科技有限公司 一种计及电价不确定性的梯级水电日前市场组合竞价方法
CN115983733A (zh) * 2023-01-03 2023-04-18 中国南方电网有限责任公司 基于水位控制的电力市场出清数据处理方法和装置
CN115983733B (zh) * 2023-01-03 2023-10-20 中国南方电网有限责任公司 基于水位控制的电力市场出清数据处理方法和装置

Also Published As

Publication number Publication date
CN114862457B (zh) 2024-08-13

Similar Documents

Publication Publication Date Title
CN112465323B (zh) 一种耦合日电量分解及日前市场竞价的梯级水电站短期鲁棒调度方法
CN101705671B (zh) 黄河上游梯级水电站运行设计与优化调度方法及设备
CN114862457A (zh) 梯级水电站群参与日前现货市场的竞价电量边界确定方法
CN111737846B (zh) 一种水电厂参与电力现货市场出清计算的方法
CN110555590A (zh) 流域梯级上下游电站参与电力现货市场二次出清方法
CN110472824A (zh) 一种考虑调峰需求的梯级水电站短期多目标优化调度方法
Wang et al. Short-term optimal scheduling of cascade hydropower plants shaving peak load for multiple power grids
CN106127342B (zh) 一种基于scuc的跨区联络线交易优化方法
CN110310022B (zh) 一种基于电力现货市场价格的电厂月度发电计划分解方法
CN107492886A (zh) 一种区域电力市场下含风电的电网月度购电方案优化方法
CN112803494B (zh) 一种含风光水火的多目标agc协调优化方法及系统
CN107506894A (zh) 一种考虑非恒定耦合约束的水电群调度方法
CN109978331B (zh) 一种高比例水电现货市场下日电量分解方法
CN106228273B (zh) 一种水电外送交易优化模型的构建方法
Ali Development and Improvement of Renewable Energy Integrated with Energy Trading Schemes based on Advanced Optimization Approaches
CN109508813B (zh) 一种兼顾调峰和弃水的水电电量分解模型建立方法
CN110570008A (zh) 流域梯级下游电站发电量全额收购的日前市场出清方法
CN110570037A (zh) 基于scuc的电厂中长期极限交易电量优化方法及系统
CN117114877A (zh) 一种基于虚拟电厂的中长期电力交易方法及系统
CN117318179A (zh) 一种基于能量共享机制的优化调度方法、设备和介质
CN115545768B (zh) 一种考虑合同分解的大水电跨省跨区日前随机竞价方法
CN110556821A (zh) 考虑交互功率控制和双边竞价交易的多微电网双层优化调度方法
CN114936672A (zh) 一种基于纳什谈判法的多虚拟电厂联合调度方法
CN112288466A (zh) 一种电力市场环境下考虑气电耦合的发电商日前申报方法
CN113283672A (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