CN109829742A - 一种基于申报价结算的火电机组最优报价计算方法及系统 - Google Patents

一种基于申报价结算的火电机组最优报价计算方法及系统 Download PDF

Info

Publication number
CN109829742A
CN109829742A CN201811322519.6A CN201811322519A CN109829742A CN 109829742 A CN109829742 A CN 109829742A CN 201811322519 A CN201811322519 A CN 201811322519A CN 109829742 A CN109829742 A CN 109829742A
Authority
CN
China
Prior art keywords
generating unit
power generating
fired power
peak regulation
quotation
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
CN201811322519.6A
Other languages
English (en)
Other versions
CN109829742B (zh
Inventor
燕京华
孙振
韩彬
胡超凡
杨军峰
罗治强
戴赛
崔晖
丁强
邹精
郑晓雨
葛睿
陈龙翔
李伟刚
李博
许丹
张传成
蔡帜
闫翠会
黄国栋
胡晨旭
张加力
门德月
刘芳
屈富敏
胡晓静
李宇轩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
State Grid Jiangxi Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
State Grid Jiangxi 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 State Grid Corp of China SGCC, China Electric Power Research Institute Co Ltd CEPRI, State Grid Jiangxi Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201811322519.6A priority Critical patent/CN109829742B/zh
Publication of CN109829742A publication Critical patent/CN109829742A/zh
Application granted granted Critical
Publication of CN109829742B publication Critical patent/CN109829742B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70Smart grids as climate change mitigation technology in the energy generation sector
    • YGENERAL 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
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS 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/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
    • YGENERAL 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
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS 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/00Market activities related to the operation of systems integrating technologies related to power network operation or related to communication or information technologies
    • Y04S50/14Marketing, i.e. market research and analysis, surveying, promotions, advertising, buyer profiling, customer management or rewards

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

一种基于申报价结算的火电机组最优报价计算方法,包括获取火电机组的参数;基于所述火电机组的参数,计算火电机组的调峰总成本;基于所述火电机组的调峰总成本和预先建立的以次日市场中标最高价格为随机变量的概率密度函数,计算火电机组最优报价。本发明提供的技术方案能够解决火电企业参与调峰辅助服务市场的竞价难题,提高火电机组参与市场的水平,提升市场运行的效率,增加火电企业的收益。

Description

一种基于申报价结算的火电机组最优报价计算方法及系统
技术领域
本发明涉及电力市场技术领域,具体涉及一种基于申报价结算的火电机组最优报价计算方法及系统。
背景技术
随着风电等清洁能源的快速发展,其装机容量在各省的总装机容量中占比越来越大。由于风电等清洁能源的反调峰特性明显,火电机组大多不愿意降出力太多,因为这会导致发电平均成本快速上升,调度中心在低谷时段的调峰压力十分巨大,调峰辅助服务市场成为了各省份解决调峰难题的重要手段。
对于调峰辅助服务市场,火电机组在最小技术出力到调峰基准功率之间的范围是该机组的有偿调峰区域。在报价时通常有两种规则,一种是已经规定了各火电机组的报价档,只需要申报价格,另一种是允许火电机组在有偿调峰区域内,自由分三档报价。两种报价模式都会要求报价单调递增。对于不同的报价模式,采取的竞价策略有所区别。对于机组参与调峰市场影响的“三公”电量进度,会在非调峰时段逐步平衡。
火电企业参与调峰市场,其收益受到多方面的影响,首先,是对自身机组参与调峰市场的成本;其次,是该企业对次日最高中标价格的预测准确度;最后,是该企业对机组进行的报价。为了在调峰市场中尽可能的提高自身的收益,增强竞争力,需要建立完善的竞价策略。
对于火电企业来说,随着机组性能的逐步老化,不断调整自身参与调峰的成本模型,在市场运营过程中,不断提高预测次日最高中标价的准确率水平,这样才能提升此报价策略的使用效果,在参与到调峰辅助服市场服务于电网安全的同时,最大化自身的收益。
发明内容
为了解决现有技术中所存在的火电企业参与调峰辅助服务市场的竞价难题,提高火电机组参与市场的水平,提升市场运行的效率,增加火电企业的收益,本发明提提供一种基于申报价结算的火电机组最优报价计算方法及系统。
本发明提供的技术方案是:
一种基于申报价结算的火电机组最优报价计算方法,其改进之处在于,所述方法包括:
获取火电机组的参数;
基于所述火电机组的参数,计算火电机组的调峰总成本;
基于所述火电机组的调峰总成本和预先建立的以次日市场中标最高价格为随机变量的概率密度函数,计算火电机组最优报价。
优选地,所述基于所述火电机组的参数,计算火电机组的调峰总成本包括:
获取该机组在最小技术出力到额定容量之间的出力点的出力值和平均煤耗,建立平均煤耗函数;
基于所述平均煤耗函数和电煤价格,建立火电机组的发电平均成本函数;
基于所述火电机组的发电平均成本函数计算火电机组的调峰总成本。
优选地,所述火电机组的发电平均成本函数如下式所示:
B=z(P)
式中,B:火电机组的发电平均成本;P:火电机组的发电出力。
优选地,所述火电机组的调峰总成本按下式进行计算:
cs(Δp)=(Fa+FΔp)-Fsta
式中,cs:火电机组的调峰总成本计算函数;Δp:火电机组调峰的中标量;Fa:火电机组调峰时段的发电成本;FΔp:火电机组参与调峰市场而少发的电量的补发成本;Fsta:火电机组在有偿调峰基准功率下发相同电量的成本。
优选地,所述火电机组调峰时段的发电成本Fa按下式进行计算:
Fa=z(Psta-Δp)(Psta-Δp)Δt;
火电机组参与调峰市场而少发的电量的补发成本FΔp按下式进行计算:
式中,Pi:本月已有的计划时刻火电机组的出力功率中去掉其中低于调峰基准功率的数据;F:去掉其中低于调峰基准功率的数据后本月已有的计划时刻火电机组的出力功率中的剩余数据的个数;
火电机组的有偿调峰基准功率Fsta按下式进行计算:
Fsta=z(Psta)PstaΔt。
优选地,所述预先建立以次日市场中标最高价格为随机变量的概率密度函数包括:
根据历史预测数据,计算每天实际中标最高价相对于预测值的偏移比例;
计算所述偏移比例的期望值和标准差,判断概率分布类型,得到以所述偏移比例为随机变量的概率密度函数;
依据次日市场中标最高价格和所述每天实际中标最高价相对于预测值的偏移比例的关系,基于所述以偏移比例为随机变量的概率密度函数计算得到以次日市场中标最高价格为随机变量的概率密度函数。
优选地,所述次日市场中标最高价格和所述每天实际中标最高价相对于预测值的偏移比例的关系如下式所示:
Mi=αiVi+Vi
式中,Mi次日市场中标最高价格:αi:每天实际中标最高价相对于预测值的偏移比例;Vi:预测次日市场中标最高价。
优选地,所述基于所述火电机组的调峰总成本和预先建立的以次日市场中标最高价格为随机变量的概率密度函数,计算火电机组最优报价包括:
市场规则中在有偿调峰区域内有多档报价;
判断市场规则中的报价模式,若报价模式规定了火电机组的报价档,则基于预先建立的以次日市场中标最高价格为随机变量的概率密度函数,计算基于所述报价档的火电机组最优报价;
若报价模式是允许火电机组在有偿调峰区域内自行分档报价,则基于所述火电机组的调峰总成本和预先建立的以次日市场中标最高价格为随机变量的概率密度函数,计算各档报价时火电机组最优报价。
优选地,所述基于所述报价档的火电机组最优报价的计算式如下:
式中,Dk:第k档的最优报价;Rk:第k档报价的期望收益;Δt:调峰的时间段时长;Ck:第k档的调峰平均成本;w:规定的报价档中每档占火电机组额定功率的比例;f(m):以次日市场中标最高价格为随机变量的概率密度函数;
其中,第k档的调峰平均成本计算按下式进行计算:
Ck=(cs(wkPcap)-cs(w(k-1)Pcap))/wPcap
式中,cs:火电机组的调峰总成本;Pcap:火电机组额定功率。
优选地,所述市场规则中在有偿调峰区域内的多档包括:第一档、第二档和第三档;
所述各档报价时火电机组最优报价按下式进行计算:
式中,R1:第一档报价的期望收益;R2:第二档报价的期望收益;R3:第三档报价的期望收益;E1:第一档功率上限;E2:第二档功率上限;C1:第一档的调峰平均成本;C2:第二档的调峰平均成本;C3:第三档的调峰平均成本;cs:火电机组的调峰总成本;D1:第一档最优报价;D2:第二档最优报价;D3:第三档最优报价;Δt:调峰的时间段时长;f(m):以次日市场中标最高价格为随机变量的概率密度函数。
一种基于申报价结算的火电机组最优报价计算系统,包括:
获取模块:用于获取火电机组的参数;
第一计算模块:用于基于所述火电机组的参数,计算火电机组的调峰总成本;
第二计算模块:用于基于所述火电机组的调峰总成本和预先建立的概率密度函数,计算火电机组最优报价。
与现有技术相比,本发明的有益效果为:
本发明提供的技术方案能够根据火电机组的发电煤耗的实际情况,计算出该机组参与调峰辅助服务市场时不同调峰量的总成本,使得火电企业能够准确掌握自身的成本;同时,能够适应不同的市场规则,提供了只申报价格和由火电企业自行在有偿调峰区间内分三段报价的两种报价模式下的最优价格计算方法,使火电企业参与市场能够获得最大收益。
本发明提供的技术方案将预测的次日最高中标价作为随机变量建立概率密度函数,减小了预测值可能带来的误差的影响,进而求出火电机组参与市场能够获得最大收益期望值的最优报价。
附图说明
图1为本发明基于申报价结算的火电机组最优报价计算方法的示意图;
图2为本发明基于申报价结算的火电机组最优报价计算方法的流程图;
图3为本发明基于申报价结算的火电机组最优报价计算系统示意图。
具体实施方式
为了更好地理解本发明,下面结合说明书附图和实例对本发明的内容做进一步的说明。
实施例1
一种基于申报价结算的火电机组最优报价计算方法,如图1所示,所述方法包括:
步骤1:获取火电机组的参数;
步骤2:基于所述火电机组的参数,计算火电机组的调峰总成本;
步骤3:基于所述火电机组的调峰总成本和预先建立的以次日市场中标最高价格为随机变量的概率密度函数,计算火电机组最优报价。
步骤1:获取火电机组的参数,包括最小技术出力到额定容量之间多个出力点的平均煤耗、有偿调峰基准功率、本月已有的计划时刻该机组的出力功率。
火电机组最优报价计算流程如图2所示。
步骤2:基于所述火电机组的参数,计算火电机组的调峰总成本,包括:
获取该机组在最小技术出力到额定容量之间的出力点的出力值和平均煤耗,建立平均煤耗函数;
基于所述平均煤耗函数和当前电煤价格,建立火电机组的发电平均成本函数;
基于所述火电机组的发电平均成本函数计算火电机组的调峰总成本。
具体地,所述火电机组的发电平均成本函数如下式所示:
B=z(P)
式中,B:火电机组的发电平均成本;P:火电机组的发电出力。
首先,获取火电机组在最小技术出力到额定容量之间多个出力点的平均煤耗。然后,根据出力值和对应的平均煤耗拟合出平均煤耗与出力点的函数关系,通常二次函数。最后,基于最新的电煤价格,建立火电机组的发电平均成本与发电出力的函数。
具体地,所述火电机组的调峰总成本按下式进行计算:
cs(Δp)=(Fa+FΔp)-Fsta
式中,cs:火电机组的调峰总成本计算函数;Δp:火电机组的调峰中标量;Fa:火电机组调峰时段的发电成本;FΔp:火电机组参与调峰市场而少发的电量的补发成本;Fsta:火电机组在有偿调峰基准功率下发相同电量的成本。
具体地,所述火电机组调峰时段的发电成本Fa按下式进行计算:
Fa=z(Psta-Δp)(Psta-Δp)Δt;
火电机组参与调峰市场而少发的电量的补发成本FΔp按下式进行计算:
式中,Pi:本月已有的计划时刻火电机组的出力功率中去掉其中低于调峰基准功率的数据;F:去掉其中低于调峰基准功率的数据后本月已有的计划时刻火电机组的出力功率中的剩余数据的个数;
火电机组的有偿调峰基准功率Fsta按下式进行计算:
Fsta=z(Psta)PstaΔt。
具体地,所述预先建立以次日市场中标最高价格为随机变量的概率密度函数包括:
根据历史预测数据,计算每天实际中标最高价相对于预测值的偏移比例;
计算所述偏移比例的期望值和标准差,判断概率分布类型,得到以所述偏移比例为随机变量的概率密度函数;
依据次日市场中标最高价格和所述每天实际中标最高价相对于预测值的偏移比例的关系,基于所述以偏移比例为随机变量的概率密度函数计算得到以次日市场中标最高价格为随机变量的概率密度函数。
具体地,所述次日市场中标最高价格和所述每天实际中标最高价相对于预测值的偏移比例的关系如下式所示:
Mi=αiVi+Vi
式中,Mi次日市场中标最高价格:αi:每天实际中标最高价相对于预测值的偏移比例;Vi:预测次日市场中标最高价。
步骤3:基于所述火电机组的调峰总成本和预先建立的以次日市场中标最高价格为随机变量的概率密度函数,计算火电机组最优报价包括:
判断市场规则中的报价模式,若报价模式规定了火电机组的报价档,则基于预先建立的以次日市场中标最高价格为随机变量的概率密度函数,计算只需要申报价格时火电机组最优报价;
具体地,所述只需要申报价格时火电机组最优报价按下式进行计算:
式中,Dk:第k档的最优报价;Rk:第k档报价的期望收益;Δt:调峰的时间段时长;Ck:第k档的调峰平均成本,w:规定的报价档中每档占火电机组额定功率的比例;Pcap:为火电机组额定功率;f(m):以次日市场中标最高价格为随机变量的概率密度函数;
其中,第k档的调峰平均成本计算按下式进行计算:
Ck=(cs(wkPcap)-cs(w(k-1)Pcap))/wPcap
式中,cs:火电机组的调峰总成本;Pcap:火电机组额定功率。
市场规则规定每wPcap为一档,求出在有偿调峰区域内该机组一共有N档报价,第k档的对应的申报调峰量范围是(wPcap(k-1),wPcapk]。其中Pcap为该机组额定功率。
若报价模式是允许火电机组在有偿调峰区域内自行分三档报价,则基于所述火电机组的调峰总成本和预先建立的以次日市场中标最高价格为随机变量的概率密度函数,计算自行分三档报价时火电机组最优报价。
设所分的三档分别为(0,E1],(E1,E2],(E2,Psta-Pmin],对应的报价分别为D1,D2,D3
具体地,所述各档报价时火电机组最优报价按下式进行计算:
式中,R1:第一档报价的期望收益;R2:第二档报价的期望收益;R3:第三档报价的期望收益;E1:第一档功率上限;E2:第二档功率上限;C1:第一档的调峰平均成本;C2:第二档的调峰平均成本;C3:第三档的调峰平均成本;cs:火电机组的调峰总成本;Psta:机组的有偿调峰基准功率;Pmin:火电机组最小技术出力;D1:第一档最优报价;D2:第二档最优报价;D3:第三档最优报价;Δt:调峰的时间段时长;f(m):以次日市场中标最高价格为随机变量的概率密度函数。
实施例2:
基于同一种发明构思,本发明还提供了一种基于申报价结算的火电机组最优报价计算系统,包括:
获取模块:用于获取火电机组的参数;
第一计算模块:用于基于所述火电机组的参数,计算火电机组的调峰总成本;第二计算模块:用于基于所述火电机组的调峰总成本和预先建立的概率密度函数,计算火电机组最优报价。
第一计算模块中,基于所述火电机组的参数,计算火电机组的调峰总成本,包括:
获取该机组在最小技术出力到额定容量之间的出力点的出力值和平均煤耗,建立平均煤耗函数;
基于所述平均煤耗函数和当前电煤价格,建立火电机组的发电平均成本函数;
基于所述火电机组的发电平均成本函数计算火电机组的调峰总成本。
具体地,所述火电机组的发电平均成本函数如下式所示:
B=z(P)
式中,B:火电机组的发电平均成本;P:火电机组的发电出力。
首先,获取火电机组在最小技术出力到额定容量之间多个出力点的平均煤耗。然后,根据出力值和对应的平均煤耗拟合出平均煤耗与出力点的函数关系,通常二次函数。最后,基于最新的电煤价格,建立火电机组的发电平均成本与发电出力的函数。
具体地,所述火电机组的调峰总成本按下式进行计算:
cs(Δp)=(Fa+FΔp)-Fsta
式中,cs:火电机组的调峰总成本计算函数;Δp:火电机组的调峰中标量;Fa:火电机组调峰时段的发电成本;FΔp:火电机组参与调峰市场而少发的电量的补发成本;Fsta:火电机组在有偿调峰基准功率下发相同电量的成本。
具体地,所述火电机组调峰时段的发电成本Fa按下式进行计算:
Fa=z(Psta-Δp)(Psta-Δp)Δt;
火电机组参与调峰市场而少发的电量的补发成本FΔp按下式进行计算:
式中,Pi:本月已有的计划时刻火电机组的出力功率中去掉其中低于调峰基准功率的数据;F:去掉其中低于调峰基准功率的数据后本月已有的计划时刻火电机组的出力功率中的剩余数据的个数;
火电机组的有偿调峰基准功率Fsta按下式进行计算:
Fsta=z(Psta)PstaΔt。
具体地,所述预先建立以次日市场中标最高价格为随机变量的概率密度函数包括:
根据历史预测数据,计算每天实际中标最高价相对于预测值的偏移比例;
计算所述偏移比例的期望值和标准差,判断概率分布类型,得到以所述偏移比例为随机变量的概率密度函数;
依据次日市场中标最高价格和所述每天实际中标最高价相对于预测值的偏移比例的关系,基于所述以偏移比例为随机变量概率密度函数计算得到以次日市场中标最高价格为随机变量的概率密度函数。
具体地,所述次日市场中标最高价格和所述每天实际中标最高价相对于预测值的偏移比例的关系如下式所示:
Mi=αiVi+Vi
式中,Mi次日市场中标最高价格:αi:每天实际中标最高价相对于预测值的偏移比例;Vi:预测次日市场中标最高价。
第二计算模块中,基于所述火电机组的调峰总成本和预先建立的以次日市场中标最高价格为随机变量的概率密度函数,计算火电机组最优报价包括:
判断市场规则中的报价模式,若报价模式规定了火电机组的报价档,则基于预先建立的以次日市场中标最高价格为随机变量的概率密度函数,计算只需要申报价格时火电机组最优报价;
具体地,所述只需要申报价格时火电机组最优报价按下式进行计算:
式中,Dk:第k档的最优报价;Rk:第k档报价的期望收益;Δt:调峰的时间段时长;Ck:第k档的调峰平均成本;w:规定的报价档中每档占火电机组额定功率的比例;Pcap:为火电机组额定功率;f(m):以次日市场中标最高价格为随机变量的概率密度函数;
其中,第k档的调峰平均成本计算按下式进行计算:
Ck=(cs(wkPcap)-cs(w(k-1)Pcap))/wPcap
式中,cs:火电机组的调峰总成本;Pcap:火电机组额定功率。
市场规则规定每wPcap为一档,求出在有偿调峰区域内该机组一共有N档报价,第k档的对应的申报调峰量范围是(wPcap(k-1),wPcapk]。其中Pcap为该机组额定功率。
若报价模式是允许火电机组在有偿调峰区域内自行分三档报价,则基于所述火电机组的调峰总成本和预先建立的以次日市场中标最高价格为随机变量的概率密度函数,计算自行分三档报价时火电机组最优报价。
设所分的三档分别为(0,E1],(E1,E2],(E2,Psta-Pmin],对应的报价分别为D1,D2,D3
具体地,所述各档报价时火电机组最优报价按下式进行计算:
式中,R1:第一档报价的期望收益;R2:第二档报价的期望收益;R3:第三档报价的期望收益;E1:第一档功率上限;E2:第二档功率上限;C1:第一档的调峰平均成本;C2:第二档的调峰平均成本;C3:第三档的调峰平均成本;cs:火电机组的调峰总成本;Psta:机组的有偿调峰基准功率;Pmin:火电机组最小技术出力;D1:第一档最优报价;D2:第二档最优报价;D3:第三档最优报价;Δt:调峰的时间段时长;f(m):以次日市场中标最高价格为随机变量的概率密度函数。
上述事实说明,本发明提供的技术方案能够根据火电机组的发电煤耗的实际情况,计算出该机组参与调峰辅助服务市场时不同调峰量的总成本,使得火电企业能够准确掌握自身的成本,为参与市场打好了基础。
本发明提供的技术方案将次日最高中标价预测值以随机变量进行建模,减小了误差,进而求出火电机组参与市场能够获得最大收益期望值的报价模型。
本发明提供的技术方案能够应对不同的报价规则,一种是规定了报价分档,只申报价格的情况,另一种是由火电企业自行在有偿调峰区间内分三段报价,使得此竞价策略能够更好的应对不同的市场规则。
显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上仅为本发明的实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均包含在申请待批的本发明的权利要求范围之内。

Claims (11)

1.一种基于申报价结算的火电机组最优报价计算方法,其特征在于,所述方法包括:
获取火电机组的参数;
基于所述火电机组的参数,计算火电机组的调峰总成本;
基于所述火电机组的调峰总成本和预先建立的以次日市场中标最高价格为随机变量的概率密度函数,计算火电机组最优报价。
2.如权利要求1所述的方法,其特征在于,所述基于所述火电机组的参数,计算火电机组的调峰总成本包括:
获取该机组在最小技术出力到额定容量之间的出力点的出力值和平均煤耗,建立平均煤耗函数;
基于所述平均煤耗函数和电煤价格,建立火电机组的发电平均成本函数;
基于所述火电机组的发电平均成本函数计算火电机组的调峰总成本。
3.如权利要求2所述的方法,其特征在于,所述火电机组的发电平均成本函数如下式所示:
B=z(P)
式中,B:火电机组的发电平均成本;P:火电机组的发电出力。
4.如权利要求1所述的方法,其特征在于,所述火电机组的调峰总成本按下式进行计算:
cs(Δp)=(Fa+FΔp)-Fsta
式中,cs:火电机组的调峰总成本计算函数;Δp:火电机组的调峰中标量;Fa:火电机组调峰时段的发电成本;FΔp:火电机组参与调峰市场而少发的电量的补发成本;Fsta:火电机组在有偿调峰基准功率下发相同电量的成本。
5.如权利要求4所述的方法,其特征在于,所述火电机组调峰时段的发电成本Fa按下式进行计算:
Fa=z(Psta-Δp)(Psta-Δp)Δt;
火电机组参与调峰市场而少发的电量的补发成本FΔp按下式进行计算:
式中,Pi:本月已有的计划时刻火电机组的出力功率中去掉其中低于调峰基准功率的数据;F:去掉其中低于调峰基准功率的数据后本月已有的计划时刻火电机组的出力功率中的剩余数据的个数;
火电机组的有偿调峰基准功率Fsta按下式进行计算:
Fsta=z(Psta)PstaΔt。
6.如权利要求1所述的方法,其特征在于,所述预先建立以次日市场中标最高价格为随机变量的概率密度函数包括:
根据历史预测数据,计算每天实际中标最高价相对于预测值的偏移比例;
计算所述偏移比例的期望值和标准差,判断概率分布类型,得到以所述偏移比例为随机变量的概率密度函数;
依据次日市场中标最高价格和所述每天实际中标最高价相对于预测值的偏移比例的关系,基于所述以偏移比例为随机变量的概率密度函数计算得到以次日市场中标最高价格为随机变量的概率密度函数。
7.如权利要求6所述的方法,其特征在于,所述次日市场中标最高价格和所述每天实际中标最高价相对于预测值的偏移比例的关系如下式所示:
Mi=αiVi+Vi
式中,Mi次日市场中标最高价格:αi:每天实际中标最高价相对于预测值的偏移比例;Vi:预测次日市场中标最高价。
8.如权利要求1所述的方法,其特征在于,所述基于所述火电机组的调峰总成本和预先建立的以次日市场中标最高价格为随机变量的概率密度函数,计算火电机组最优报价包括:
市场规则中在有偿调峰区域内有多档报价;
判断市场规则中的报价模式,若报价模式规定了火电机组的报价档,则基于预先建立的以次日市场中标最高价格为随机变量的概率密度函数,计算基于所述报价档的火电机组最优报价;
若报价模式是允许火电机组在有偿调峰区域内自行分档报价,则基于所述火电机组的调峰总成本和预先建立的以次日市场中标最高价格为随机变量的概率密度函数,计算各档报价时火电机组最优报价。
9.如权利要求8所述的方法,其特征在于,所述基于所述报价档的火电机组最优报价的计算式如下:
式中,Dk:第k档的最优报价;Rk:第k档报价的期望收益;Δt:调峰的时间段时长;Ck:第k档的调峰平均成本;w:规定的报价档中每档占火电机组额定功率的比例;f(m):以次日市场中标最高价格为随机变量的概率密度函数;
其中,第k档的调峰平均成本计算按下式进行计算:
Ck=(cs(wkPcap)-cs(w(k-1)Pcap))/wPcap
式中,cs:火电机组的调峰总成本;Pcap:火电机组额定功率。
10.如权利要求9所述的方法,其特征在于,所述市场规则中在有偿调峰区域内的多档包括:第一档、第二档和第三档;
所述各档报价时火电机组最优报价按下式进行计算:
式中,R1:第一档报价的期望收益;R2:第二档报价的期望收益;R3:第三档报价的期望收益;E1:第一档功率上限;E2:第二档功率上限;C1:第一档的调峰平均成本;C2:第二档的调峰平均成本;C3:第三档的调峰平均成本;cs:火电机组的调峰总成本;D1:第一档最优报价;D2:第二档最优报价;D3:第三档最优报价;Δt:调峰的时间段时长;f(m):以次日市场中标最高价格为随机变量的概率密度函数。
11.一种基于申报价结算的火电机组最优报价计算系统,其特征在于,所述系统包括:
获取模块:用于获取火电机组的参数;
第一计算模块:用于基于所述火电机组的参数,计算火电机组的调峰总成本;
第二计算模块:用于基于所述火电机组的调峰总成本和预先建立的概率密度函数,计算火电机组最优报价。
CN201811322519.6A 2018-11-08 2018-11-08 一种基于申报价结算的火电机组最优报价计算方法及系统 Active CN109829742B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811322519.6A CN109829742B (zh) 2018-11-08 2018-11-08 一种基于申报价结算的火电机组最优报价计算方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811322519.6A CN109829742B (zh) 2018-11-08 2018-11-08 一种基于申报价结算的火电机组最优报价计算方法及系统

Publications (2)

Publication Number Publication Date
CN109829742A true CN109829742A (zh) 2019-05-31
CN109829742B CN109829742B (zh) 2021-12-03

Family

ID=66859796

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811322519.6A Active CN109829742B (zh) 2018-11-08 2018-11-08 一种基于申报价结算的火电机组最优报价计算方法及系统

Country Status (1)

Country Link
CN (1) CN109829742B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111027802A (zh) * 2019-06-06 2020-04-17 国网辽宁省电力有限公司 一种基于竞价辅助服务市场的自动发电控制方法及装置
CN111222975A (zh) * 2019-11-18 2020-06-02 中国银行股份有限公司 一种本币利率互换套算方法及系统
CN111740413A (zh) * 2020-06-28 2020-10-02 东北电力大学 一种考虑火电调峰主动性与需求响应的含储能电力系统优化调度方法
CN112465534A (zh) * 2020-10-20 2021-03-09 华电电力科学研究院有限公司 一种热电机组参与调峰辅助服务竞价计算方法
CN117196683A (zh) * 2023-09-26 2023-12-08 国网宁夏电力有限公司 基于lstm的火电企业月内滚动撮合交易竞价系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010081165A2 (en) * 2009-01-12 2010-07-15 Battelle Memorial Institute Nested, hierarchical resource allocation schema for management and control of an electric power grid
CN106844916A (zh) * 2017-01-09 2017-06-13 北京清能互联科技有限公司 一种基于调峰辅助服务市场的发电和用电组织方法及装置
CN106845807A (zh) * 2017-01-09 2017-06-13 北京清能互联科技有限公司 基于调峰辅助服务的结算方法及装置
CN108123492A (zh) * 2017-12-26 2018-06-05 国网辽宁省电力有限公司 一种考虑火电电蓄热联合调峰的日前发电计划优化方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010081165A2 (en) * 2009-01-12 2010-07-15 Battelle Memorial Institute Nested, hierarchical resource allocation schema for management and control of an electric power grid
CN106844916A (zh) * 2017-01-09 2017-06-13 北京清能互联科技有限公司 一种基于调峰辅助服务市场的发电和用电组织方法及装置
CN106845807A (zh) * 2017-01-09 2017-06-13 北京清能互联科技有限公司 基于调峰辅助服务的结算方法及装置
CN108123492A (zh) * 2017-12-26 2018-06-05 国网辽宁省电力有限公司 一种考虑火电电蓄热联合调峰的日前发电计划优化方法

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111027802A (zh) * 2019-06-06 2020-04-17 国网辽宁省电力有限公司 一种基于竞价辅助服务市场的自动发电控制方法及装置
CN111027802B (zh) * 2019-06-06 2023-07-25 国网辽宁省电力有限公司 一种基于竞价辅助服务市场的自动发电控制方法及装置
CN111222975A (zh) * 2019-11-18 2020-06-02 中国银行股份有限公司 一种本币利率互换套算方法及系统
CN111740413A (zh) * 2020-06-28 2020-10-02 东北电力大学 一种考虑火电调峰主动性与需求响应的含储能电力系统优化调度方法
CN111740413B (zh) * 2020-06-28 2021-11-30 东北电力大学 一种考虑火电调峰主动性的含储能电力系统优化调度方法
CN112465534A (zh) * 2020-10-20 2021-03-09 华电电力科学研究院有限公司 一种热电机组参与调峰辅助服务竞价计算方法
CN117196683A (zh) * 2023-09-26 2023-12-08 国网宁夏电力有限公司 基于lstm的火电企业月内滚动撮合交易竞价系统

Also Published As

Publication number Publication date
CN109829742B (zh) 2021-12-03

Similar Documents

Publication Publication Date Title
CN109829742A (zh) 一种基于申报价结算的火电机组最优报价计算方法及系统
US20210344225A1 (en) Energy storage modeling and control
Zhai et al. Modeling and analysis of lithium battery operations in spot and frequency regulation service markets in Australia electricity market
KR101490611B1 (ko) 풍력 발전 시스템, 풍력 발전 제어 장치 및 풍력 발전 제어 방법
CN105046395B (zh) 一种含多类型新能源的电力系统日内滚动计划编制方法
Luickx et al. Considerations on the backup of wind power: Operational backup
CN102567815A (zh) 基于实际电网运行数据的事后理想计划分析方法
CN108281989B (zh) 一种风电电力系统经济调度方法及装置
US20220029424A1 (en) Hybrid power plant
Taylor et al. Consolidated dynamic pricing of power system regulation
CN103887813A (zh) 基于风功率预测不确定度的风电系统运行的控制方法
CN103280845B (zh) 一种水火电系统跨省调峰负荷分配方法
Mégel et al. Stochastic dual dynamic programming to schedule energy storage units providing multiple services
JP2021141778A (ja) 計画作成装置、計画作成方法およびコンピュータプログラム
Hoffman et al. Analysis tools for sizing and placement of energy storage in grid applications
CN111008463A (zh) 一种考虑发电侧储能的容量配置优化方法、装置及设备
CN102904248A (zh) 基于风电出力不确定性集合的电力系统调度方法
Park et al. Optimal VPP Operation Considering Network Constraint Uncertainty of DSO
CN103490421B (zh) 区域电网直调抽水蓄能电站群短期多电网负荷分配方法
Simpkins et al. Optimizing battery sizing and dispatching to maximize economic return
CN107833101A (zh) 一种基于节点电价机制的可再生能源交易费用的计算方法
CN110880777A (zh) 一种电网侧储能应急调峰备用能力的评估方法和装置
Steber Integration of decentralized battery energy storage systems into the german electrical power system
Shang et al. Grid-Side Energy Storage System Day-Ahead Bidding Strategy Based on Two-Level Decision in Spot Market
CN113011030B (zh) 基于cps1的调频容量分配方法、装置及存储介质

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