CN113610394A - 一种基于区域互联的能源市场双边竞价出清方法 - Google Patents

一种基于区域互联的能源市场双边竞价出清方法 Download PDF

Info

Publication number
CN113610394A
CN113610394A CN202110900365.XA CN202110900365A CN113610394A CN 113610394 A CN113610394 A CN 113610394A CN 202110900365 A CN202110900365 A CN 202110900365A CN 113610394 A CN113610394 A CN 113610394A
Authority
CN
China
Prior art keywords
energy
market
power
clearing
region
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
CN202110900365.XA
Other languages
English (en)
Other versions
CN113610394B (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.)
Tianjin University
Nari Technology Co Ltd
Original Assignee
Tianjin University
Nari Technology 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 Tianjin University, Nari Technology Co Ltd filed Critical Tianjin University
Priority to CN202110900365.XA priority Critical patent/CN113610394B/zh
Publication of CN113610394A publication Critical patent/CN113610394A/zh
Application granted granted Critical
Publication of CN113610394B publication Critical patent/CN113610394B/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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06315—Needs-based resource requirements planning or analysis
    • 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
    • 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
    • G06F2113/00—Details relating to the application field
    • G06F2113/04—Power grid distribution networks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/06—Wind turbines or wind farms
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08—Thermal analysis or thermal optimisation

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Physics & Mathematics (AREA)
  • Development Economics (AREA)
  • Theoretical Computer Science (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Physics & Mathematics (AREA)
  • Finance (AREA)
  • Accounting & Taxation (AREA)
  • Game Theory and Decision Science (AREA)
  • Marketing (AREA)
  • General Business, Economics & Management (AREA)
  • Operations Research (AREA)
  • Tourism & Hospitality (AREA)
  • Quality & Reliability (AREA)
  • Educational Administration (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开一种基于区域互联的能源市场双边竞价出清方法，包括基于零售侧能源市场，构建互联区域电热联合市场框架；以各区域总社会剩余最大化为目标，提出互联区域能源市场双边竞价出清策略；采用内点法对上述模型进行求解，根据所得到市场出清结果对分布式能源站供能、用户用能行为及区域间传输能量进行调节。该方法为能源互联区域提供了能源市场的交易机制，可实现资源的跨区域优化分配，同时提升各区域市场参与主体的利益。

Description

一种基于区域互联的能源市场双边竞价出清方法
技术领域
本发明涉及综合能源系统能源市场领域，尤其涉及一种基于区域互联的能源市场双边竞价出清方法。
背景技术
随着国民经济的迅速发展，以化石能源为主的能源消费体系使得我国的能源短缺和环境问题日益凸显。因此，综合能源系统(integrated energy system,IES)的概念得以提出和推广，特指在规划、建设和运行等过程中，通过对能源的产生/传输/分配(能源供应网络)、转换、存储、消费等环节进行有机协调与优化后形成的能源产供销一体化系统，被认为是未来人类社会能源的主要承载形式。IES中丰富的可调度资源可满足用户多种能源梯级利用，提升社会供能的安全可靠性。同时，随着电力体制改革的推进和相关支撑技术的发展，能源市场对于资源的优化调节作用日益凸显，以电能为主的综合能源交易迅速发展。目前对于能源市场的研究大多局限于电力市场，随着IES的普及与应用，单一能源服务方式已无法满足用户的多能源需求。同时，随着能源耦合更加紧密，区域间的能源互联成为趋势，研究互联区域的能源市场出清策略对于资源优化分配，提升系统运行可靠性与经济性，促进供需双方利益提升具有重要意义。
互联区域能源市场出清策略研究的目标是使不同区域的供需双方参与统一的能源市场，通过能源市场中的竞价与分析，得到使所有参与主体共同利益提升的出清结果，确定此时的出清价格与能源数量，并根据此结果调节分布式能源站与用户的供用能行为。基于双边竞价的出清策略是能源市场顺利交易的基础，是各利益主体持续参与市场并保持积极性的关键。
因此，有必要基于零售侧能源市场，研究互联区域的能源市场双边竞价出清策略。如何对互联分布式能源站运行机理进行建模，如何在更高的层面上实现跨区域的资源优化及各区域市场主体共同利益的最大化，是亟待解决的问题。
发明内容
本发明的目的是为了克服现有技术中的不足，提供一种基于区域互联的能源市场双边竞价出清方法，该方法可以有效的解决综合能源系统中跨区域的资源优化分配问题，实现能源市场交易后各区域总社会剩余的最大化。
本发明的目的是通过以下技术方案实现的：
一种基于区域互联的能源市场双边竞价出清方法，包括以下步骤：
基于零售侧能源市场，构建互联区域电热联合市场框架；
以市场交易后各区域的总社会剩余最大化为目标，建立互联区域能源市场双边竞价出清模型；
采用内点法对互联区域能源市场双边竞价出清模型进行求解得到市场出清结果，根据所得到市场出清结果对分布式能源站供能、用户用能行为及区域间传输能量进行调节；
对于分布式能源站，调节站内各种能源耦合设备的出力以响应市场出清结果，同时通过控制管道阀门与各控制装置调节区域间所传输的热能与电能；对于多能用户，根据市场出清结果与自身需求调节负荷情况。
进一步的，所述互联区域电热联合市场框架具体为：
分布式能源站与用户分别为能源的生产者与消费者，通过市场交易实现资源的优化分配；在需求侧，引入负荷多代理系统，即通过电、热负荷代理帮助用户进行竞价，负责收集相应区域内的市场与用户信息计算竞价价格，并聚合形成各自的需求曲线，提交至电热联合能源市场进行竞价；而供给侧则由综合能源服务商负责接收上级能源网络的能源价格信息及分布式能源站的设备信息，对设备供能边际成本进行排序并提交至市场竞价；各区域能源市场根据供需竞价信息，以各区域总社会剩余最大为目标进行出清，决定此时的能源出清数量与价格以及区域间传输的能源数量；在市场出清后，由综合能源服务商根据出清结果发出控制信号，调节分布式能源站供能与用户用能行为。
进一步的，建立互联区域能源市场双边竞价出清模型的具体步骤如下：
(1)分段讨论市场出清情况；
(2)设定市场出清模型的目标函数为最大化社会剩余；
(3)设定分布式能源站中各设备的市场出清的约束条件。
进一步的，具体的：
市场出清模型的目标函数如下
Figure BDA0003199571140000021
其中，i表示第i个互联区域，n为互联区域数量，
Figure BDA0003199571140000022
和
Figure BDA0003199571140000023
分别为电力和热力社会剩余，
Figure BDA0003199571140000024
和
Figure BDA0003199571140000025
分别为第i个分布式能源站中供应电能、热能功率，
Figure BDA0003199571140000026
和
Figure BDA0003199571140000027
分别为市场出清电、热负荷需求，
Figure BDA0003199571140000031
和
Figure BDA0003199571140000032
分别为市场出清电、热能源价格，
Figure BDA0003199571140000033
和
Figure BDA0003199571140000034
分别为第i个DES中第k种和第l种设备产能的边际成本；
Figure BDA0003199571140000035
和
Figure BDA0003199571140000036
为第i个DES中第k种和第l种设备出力；
Figure BDA0003199571140000037
和
Figure BDA0003199571140000038
为描述第i个DES中第k种和第l种设备选取情况的二元变量，数值为1表示选择该设备产能，数值为0则表示不使用该设备；
Figure BDA0003199571140000039
和
Figure BDA00031995711400000310
分别为电、热负荷功率下限，
Figure BDA00031995711400000311
和
Figure BDA00031995711400000312
分别为电力；WT表示风机、HP表示热泵、CHP表示热电联产、GF表示燃气锅炉、HE表示热交换器；
市场出清模型的约束条件为互联区域分布式能源站的能量平衡约束和网络潮流及设备出力的不等式约束；互联区域分布式能源站能量平衡约束为：
Figure BDA00031995711400000313
Figure BDA00031995711400000314
其中，Le和Lh分别为电力和热力负荷功率；
Figure BDA00031995711400000315
和
Figure BDA00031995711400000316
分别为热电联产利用天然气生产电能和热能的效率；ηGF为燃气锅炉利用天然气生产热能的效率；ηPEC和ηHE分别为电力电子转换器和热交换器的能量转换效率；ηHP为热泵利用电能生产热能的效率；Pe、Pg、Ph和PWT分别为上级电力、天然气、热力网络及风机的功率；λ1为分配系数，表示流经热电联产天然气占输入总量的比例；λ2为分配系数，表示流经电力电子转换器电能占输入总量的比例；
不等式约束为：
Figure BDA00031995711400000317
Figure BDA00031995711400000318
Pg_min≤Pg≤Pg_max
PHP_min≤PHP≤PHP_max
Le_min≤Le≤Le_max
Lh_min≤Lh≤Lh_max
其中，Vi p为节点i的p相电压，Vi p_max和Vi p_min代表其上、下限；Sij为线路ij上所传输的视在功率，
Figure BDA0003199571140000041
代表其上限；
Figure BDA0003199571140000042
为管道中的水流量，
Figure BDA0003199571140000043
和
Figure BDA0003199571140000044
代表其上、下限；Pg_min和Pg_max分别为天然气输入功率Pg的下限及上限；PHP为热泵设备输出功率，PHP_min和PHP_max分别为其下限及上限；Le_max和Le_min分别为电负荷功率Le的上、下限，Lh_max和Lh_min分别为热负荷功率Lh的上、下限。
与现有技术相比，本发明的技术方案所带来的有益效果是：
1、本发明构建了有助于提升互联区域能源供需双方利益的能源市场出清策略，可提升市场主体参与积极性，充分发挥市场自发调节作用；
2、通过互联区域能源市场出清，可确定每个时刻的能源价格与供需数量，从而实现跨区域的资源优化分配，有助于提升系统的经济性。
3、基于能源市场交易所实现的资源优化分配，可有效地实现综合需求响应，进而促进可再生能源消纳率，有助于能源的清洁化使用。
附图说明
图1为互联分布式能源站运行示意图；
图2为电热水器电热耦合示意图；
图3a和图3b为电热用户竞价曲线；
图4为电热联合市场框架图；
图5为能源市场社会剩余示意图；
图6a和图6b为独立区域能源市场出清过程示意图；
图7a和图7b为互联区域能源市场出清过程示意图；
图8为能源市场出清算法流程图。
具体实施方式
以下结合附图和具体实施例对本发明作进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
为了实现各区域资源优化分配、提升市场参与主体利益，本发明实施例提供了一种基于区域互联的能源市场双边竞价出清方法，详见下文描述：
综上所述，本发明实施例有助于提升市场参与主体利益，可以指导能源市场制定最有利于资源优化分配的出清策略。
1.基于零售侧能源市场，构建互联区域电热联合市场框架；
首先简要介绍能源供应者分布式能源站。综合能源系统中，分布式能源站(DES)是一种分布在用户侧，将用户与各能源网络紧密相连的能源耦合形式，是满足用户多样化用能需求并构成区域多能互补能源系统的重要载体。本发明以A型与B型两种典型DES为例进行分析，其主要设备包括：风机(wind turbine,WT)、电力电子转换器(power electronicconverter,PEC)、热泵(heat pump,HP)、热电联产(combined heat and power，CHP)、燃气锅炉(gas furnace，GF)、热交换器(heat exchanger,HE)。
在互联DES的运行过程中，通过在每个DES的输出端设置控制装置和阀门来控制能源的传输，根据能否满足该区域用户的能源需求及实现各区域利益提升来决定各区域能源的优化分配，其运行示意图如图1所示。互联A型与B型DES的能量平衡方程分别为：
Figure BDA0003199571140000051
Figure BDA0003199571140000052
其中，Le和Lh分别为电力和热力负荷功率；
Figure BDA0003199571140000053
和
Figure BDA0003199571140000054
分别为CHP利用天然气生产电能和热能的效率；ηGF为GF利用天然气生产热能的效率；ηPEC和ηHE分别为PEC和HE的能量转换效率；ηHP为HP利用电能生产热能的效率；Pe、Pg、Ph和PWT分别为上级电力、天然气、热力网络及风机的功率；λ1为分配系数，表示流经CHP天然气占输入总量的比例；λ2为分配系数，表示流经PEC电能占输入总量的比例。
接着介绍能源消费者用户负荷模型。本实施例所研究用户能源需求主要包括电能与热能，所研究负荷分别为电热水器与空间供暖负荷。对于电热水器，采用等值热力学参数法建立电热水器设备数学模型，模拟其运行过程。其电热耦合关系如图2所示，其水温的动态变化可用下式表示：
1)当电热水器工作Δt时间后，水箱中的水温可通过下式计算：
Figure BDA0003199571140000061
2)当电热水器关闭Δt时间后，水箱中的水温可通过下式计算：
Figure BDA0003199571140000062
3)当电热水器处于关闭状态且用户使用热水时，水箱中水温
Figure BDA0003199571140000063
计算如下：
Figure BDA0003199571140000064
式中：t0表示当前时间；Ttank为水箱内实际水温；Tam为外界环境温度；R为电热水器等值热阻；C为热容；Qewh为电热水器操作热功率；M为水箱中的水量；mu为使用的水量。
对于空间供暖负荷，其运行状态受室外温度，建筑材料，太阳辐射和电气设备加热等方面的影响。假设单栋楼宇中所有用户具有相同的室内供热温度和温度设定值，居民楼宇的热平衡方程可表示如下式：
Figure BDA0003199571140000065
式中：ρ为空气密度；Ca为空气的比热容；V为居民楼宇有效供热体积；
Figure BDA0003199571140000066
为第i栋楼宇的室内供热温度；ΔQi为楼宇i的热增益，若其值为正，表示楼宇获得热量，温度上升；若为负，表示楼宇损失热量，温度降低。
考虑到维护结构的热量损耗以及太阳辐射带来的热增益，居民楼宇热平衡方程也可用下式表示：
Figure BDA0003199571140000067
式中：
Figure BDA0003199571140000068
和
Figure BDA0003199571140000069
为第i栋楼宇围护结构(墙体、窗户和地板)的热量消耗；
Figure BDA00031995711400000610
为太阳辐射得热量；QDES为分布式能源站输入热量。
围护结构(墙体、窗户、地板)的耗热量可以通过下式获得：
Figure BDA00031995711400000611
Figure BDA00031995711400000612
Figure BDA00031995711400000613
式中：β为热传输系数；S为维护结构面积；下标walls、win、floor分别表示墙体、窗户及地板；
Figure BDA00031995711400000614
为第i栋楼宇的室内温度设定值；
Figure BDA00031995711400000615
为第i栋楼宇的室内温度设计值；γ和δ为建筑物的温差修正系数和方向修正系数，用于表示居民楼宇方位对热量损失的影响。
太阳辐射得热量
Figure BDA0003199571140000071
可由下式计算：
Figure BDA0003199571140000072
式中：αwin为玻璃透射系数，SA为窗户遮阳系数；Swin为楼宇每个墙面窗户总面积；SQ为总太阳辐射量。
接着介绍用户在双边竞价市场中的竞价策略，其示意图如图3a和图3b所示。在市场中引入负荷多代理系统，由电、热负荷代理收集相应区域内的市场与用户信息计算竞价价格，并聚合形成各自的需求曲线，提交至电热联合能源市场。对于电力用户，通过市场预测电能价格的平均值和标准差及收集区域内负荷状态(由当前水温和所需水温之间的差值定义)，计算电力用户竞价价格，其与实际水温Ttank呈线性负相关关系，具体如下所示：
Figure BDA0003199571140000073
其中，pbid为电力用户竞价价格，pa和psd为未来两小时预测电价平均值与标准差，c为用户预定义舒适度，Tset为温度设定值；Tmax为最高水温。
对于热力用户，与电力用户竞价价格类似，热力用户竞价价格也与市场预测热价等因素有关，并且与热功率输入QDES呈线性负相关关系，具体可表示为：
Figure BDA0003199571140000074
Figure BDA0003199571140000075
其中，ph_high和ph_low分别为热力竞价价格的上下限，ph_a和ph_sd分别为未来两小时预测热价的平均值与标准差，ch为用户预定义热舒适度，
Figure BDA0003199571140000076
为第i栋楼宇用户期望热量输入，
Figure BDA0003199571140000077
和
Figure BDA0003199571140000078
分别为各建筑物最大和最小的热输入功率，可利用下式计算：
Figure BDA0003199571140000079
Figure BDA0003199571140000081
以所建立用户电热负荷模型(即公式(3)-(11))与竞价模型(即公式(12)-(16))为基础，构建综合能源系统中互联区域的电热联合零售能源市场框架，其示意图如图4所示，具体为：DES与用户分别为能源的生产者与消费者，通过市场交易实现资源的优化分配。在需求侧，引入负荷多代理系统，由电、热负荷代理收集相应区域内的市场与用户信息计算竞价价格，并聚合形成各自的需求曲线，提交至电热联合能源市场进行竞价；而供给侧则由综合能源服务商负责接收上级能源网络的能源价格信息及DES的设备信息，对设备供能边际成本进行排序并提交至市场竞价。各区域能源市场根据供需竞价信息，以各区域总社会剩余最大为目标进行出清，决定此时的能源出清数量与价格以及区域间传输的能源数量。在市场出清后，由综合能源服务商根据出清结果发出控制信号，调节DES供能与用户用能行为。
2.提出互联区域能源市场双边竞价出清模型，指导能源市场以各区域总社会剩余最大为目标进行出清：
双边竞价市场的出清结果通常由供需曲线决定。对于需求曲线，认为其由一段代表不可控负荷的水平线段与代表可控负荷的递减曲线组成；而供应曲线，则认为其是由DES内各设备的供能边际成本所决定的阶梯型线段。供需曲线中价格与功率的关系可由下式表示：
Figure BDA0003199571140000082
Figure BDA0003199571140000083
其中，ps为供应价格，pmax和pmin分别为供应价格的上下限，qs为供应能源数量，qw为最低供能成本设备供应能源数量，qg为次低供能成本设备供应能源数量，pd为需求侧竞价价格，d为需求曲线斜率，qd为需求负荷量，qun为不可控负荷需求量，qr为可调节负荷量，pms为需求侧最高竞价价格。
上式可进一步整理为：
Ax＝b (19)
其中，A为供需矩阵，向量x为供需能源价格及数量，b为能源价格与数量的解向量，具体如下：
Figure BDA0003199571140000091
其中，Δq为传输能源，当Δq>0时，代表向其他区域传输能源，反之为接受能源；而当Δp>0时，其代表价格补贴，反之则代表能源税。
本发明通过社会剩余来表示市场参与者在交易过程中获得的收益，在几何上表示为供需曲线所围部分了，其示意图如图5所示。
图6a至土7b中展示了两区域在两种不同情况下双边竞价能源市场的出清过程。图6a和图6b中，DES仅为独立区域供能，市场在供需曲线的交点处出清，此时供应能源与价格都与需求侧相等，区域1的出清价格比区域2的低。出清结果可表示为：
x1＝A-1b0＝[qc1 pc1 qc1 pc1]T (21)
x2＝A-1b0＝[qc2 pc2 qc2 pc2]T (22)
其中，x1和x2分别表示区域1和区域2的供需能源价格及数量，qc1表示区域1的能源供应量，此时与负荷需求量相等；qc2表示区域2的负荷需求量与能源需求量，此时与负荷需求量相等；pc1表示区域1的竞价价格，此时与供应价格相等；pc2表示区域2的竞价价格，此时与供应价格相等。并且此时的解向量为
Figure BDA0003199571140000092
图7a和图7b中，DES为互联区域供能，若不考虑传输线约束，区域1DES多余的能源通过传输线路及管道供应区域2中额外负荷，此时两区域能源市场出清点均沿需求曲线发生移动，直到两个区域的出清价格相等。此时，区域1的出清价格升高，用户能源需求降低，而区域2的出清价格降低，用户能源需求升高。由于出清点的移动，图形中出现供需曲线间的三角形部分，其表示互联所增加的社会剩余。对于区域1，DES向其他区域传输能量，市场出清点上升，这意味着该DES可以更高的价格出售能源，此时，增加的社会剩余为生产者剩余；对于区域2，DES接受从其他区域传输来的能源，市场出清点下降，这意味着该区域的用户可以较低的价格购买能源，此时，社会剩余的增加为消费者剩余。市场出清结果可表示为：
x1＝A-1b1＝[qc3 pct qc1 pct]T (24)
x2＝A-1b2＝[qc4 pct qc2 pct]T (25)
其中，pct为区域统一出清价格，qc3、qc3分别代表区域1与区域2的负荷需求量，且此时区域1与区域2的解向量分别为：
Figure BDA0003199571140000101
基于以上分析，本发明以各区域总社会剩余最大化为目标，对互联区域内的供需资源进行优化分配。具体目标函数如下式所示：
Figure BDA0003199571140000102
其中，上标i表示第i个互联区域，n为互联区域数量，
Figure BDA0003199571140000103
和
Figure BDA0003199571140000104
分别为电力和热力社会剩余，
Figure BDA0003199571140000105
和
Figure BDA0003199571140000106
分别为第i个DES中供应电能、热能功率，
Figure BDA0003199571140000107
和
Figure BDA0003199571140000108
分别为市场出清电、热负荷需求，
Figure BDA0003199571140000109
和
Figure BDA00031995711400001010
分别为市场出清电、热能源价格，
Figure BDA00031995711400001011
和
Figure BDA00031995711400001012
为第i个DES中第k种和第l种设备产能的边际成本；
Figure BDA00031995711400001013
和
Figure BDA00031995711400001014
为第i个DES中第k种和第l种设备出力；
Figure BDA00031995711400001015
和
Figure BDA00031995711400001016
为描述第i个DES中第k种和第l种设备选取情况的二元变量，数值为1表示选择该设备产能，数值为0则表示不使用该设备；
Figure BDA00031995711400001017
和
Figure BDA00031995711400001018
分别为电、热负荷功率下限，
Figure BDA00031995711400001019
和
Figure BDA00031995711400001020
分别为电力、热力竞价价格上限。
市场出清模型的约束条件为互联区域分布式能源站的能量平衡约束和网络潮流及设备出力等不等式约束；
互联区域分布式能源站能量平衡约束为：
Figure BDA00031995711400001021
Figure BDA0003199571140000111
不等式约束为：
Pe_min≤Pe≤Pe_max (30)
Ph_min≤Ph≤Ph_max (31)
Pg_min≤Pg≤Pg_max (32)
PHP_min≤PHP≤PHP_max (33)
Le_min≤Le≤Le_max (34)
Lh_min≤Lh≤Lh_max (35)
其中，Pe_min、Ph_min和Pg_min分别为上级网络电能、热能和天然气的输入功率下限；Pe_max、Ph_max和Pg_max分别为上级网络电能、热能和天然气的输入功率上限；PHP为热泵设备输出功率，PHP_min和PHP_max分别为其下限及上限；Le_max和Le_min分别为电负荷功率上、下限，Lh_max和Lh_min分别为热负荷功率上、下限。
3.根据以上能源市场架构与双边竞价市场出清模型，可利用内点法求得市场出清结果，算法流程如图8所示。对出清结果的调节过程具体为：根据所得到市场出清结果，决定此时能源价格与供需能源数量，由综合能源服务商向分布式能源站与用户传递控制信号，调整能源站内各设备出力，调节能源站输出端控制装置及阀门以控制区域间传输能量，并通过对负荷的直接控制调节用户能源需求。
综上所述，本发明实施例提供了一种基于区域互联的能源市场双边竞价出清方法，有助于提升互联区域能源供需双方利益的能源市场出清策略，可提升市场主体参与积极性，充分发挥市场自发调节作用，实现跨区域的资源优化分配，可直接应用于互联区域的多能源市场中，辅助能源市场顺利交易。
本发明实施例对各器件的型号除做特殊说明的以外，其他器件的型号不做限制，只要能完成上述功能的器件均可。
本领域技术人员可以理解附图只是一个优选实施例的示意图，上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。
本发明并不限于上文描述的实施方式。以上对具体实施方式的描述旨在描述和说明本发明的技术方案，上述的具体实施方式仅仅是示意性的，并不是限制性的。在不脱离本发明宗旨和权利要求所保护的范围情况下，本领域的普通技术人员在本发明的启示下还可做出很多形式的具体变换，这些均属于本发明的保护范围之内。

Claims (4)

1.一种基于区域互联的能源市场双边竞价出清方法，其特征在于，包括以下步骤：
基于零售侧能源市场，构建互联区域电热联合市场框架；
以市场交易后各区域的总社会剩余最大化为目标，建立互联区域能源市场双边竞价出清模型；
采用内点法对互联区域能源市场双边竞价出清模型进行求解得到市场出清结果，根据所得到市场出清结果对分布式能源站供能、用户用能行为及区域间传输能量进行调节；
对于分布式能源站，调节站内各种能源耦合设备的出力以响应市场出清结果，同时通过控制管道阀门与各控制装置调节区域间所传输的热能与电能；对于多能用户，根据市场出清结果与自身需求调节负荷情况。
2.根据权利要求1所述的一种基于区域互联的能源市场双边竞价出清方法，其特征在于，所述互联区域电热联合市场框架具体为：
分布式能源站与用户分别为能源的生产者与消费者，通过市场交易实现资源的优化分配；在需求侧，引入负荷多代理系统，即通过电、热负荷代理帮助用户进行竞价，负责收集相应区域内的市场与用户信息计算竞价价格，并聚合形成各自的需求曲线，提交至电热联合能源市场进行竞价；而供给侧则由综合能源服务商负责接收上级能源网络的能源价格信息及分布式能源站的设备信息，对设备供能边际成本进行排序并提交至市场竞价；各区域能源市场根据供需竞价信息，以各区域总社会剩余最大为目标进行出清，决定此时的能源出清数量与价格以及区域间传输的能源数量；在市场出清后，由综合能源服务商根据出清结果发出控制信号，调节分布式能源站供能与用户用能行为。
3.根据权利要求1所述的一种基于区域互联的能源市场双边竞价出清方法，其特征在于，建立互联区域能源市场双边竞价出清模型的具体步骤如下：
(1)分段讨论市场出清情况；
(2)设定市场出清模型的目标函数为最大化社会剩余；
(3)设定分布式能源站中各设备的市场出清的约束条件。
4.根据权利要求3所述的一种基于区域互联的能源市场双边竞价出清方法，其特征在于，具体的：
市场出清模型的目标函数如下
Figure FDA0003199571130000011
其中，i表示第i个互联区域，n为互联区域数量，
Figure FDA0003199571130000021
和
Figure FDA0003199571130000022
分别为电力和热力社会剩余，
Figure FDA0003199571130000023
和
Figure FDA0003199571130000024
分别为第i个分布式能源站中供应电能、热能功率，
Figure FDA0003199571130000025
和
Figure FDA0003199571130000026
分别为市场出清电、热负荷需求，
Figure FDA0003199571130000027
和
Figure FDA0003199571130000028
分别为市场出清电、热能源价格，
Figure FDA0003199571130000029
和
Figure FDA00031995711300000210
分别为第i个DES中第k种和第l种设备产能的边际成本；
Figure FDA00031995711300000211
和
Figure FDA00031995711300000212
为第i个DES中第k种和第l种设备出力；
Figure FDA00031995711300000213
和
Figure FDA00031995711300000214
为描述第i个DES中第k种和第l种设备选取情况的二元变量，数值为1表示选择该设备产能，数值为0则表示不使用该设备；
Figure FDA00031995711300000215
和
Figure FDA00031995711300000216
分别为电、热负荷功率下限，
Figure FDA00031995711300000217
和
Figure FDA00031995711300000218
分别为电力；WT表示风机、HP表示热泵、CHP表示热电联产、GF表示燃气锅炉、HE表示热交换器；
市场出清模型的约束条件为互联区域分布式能源站的能量平衡约束和网络潮流及设备出力的不等式约束；互联区域分布式能源站能量平衡约束为：
Figure FDA00031995711300000219
Figure FDA00031995711300000220
其中，Le和Lh分别为电力和热力负荷功率；
Figure FDA00031995711300000221
和
Figure FDA00031995711300000222
分别为热电联产利用天然气生产电能和热能的效率；ηGF为燃气锅炉利用天然气生产热能的效率；ηPEC和ηHE分别为电力电子转换器和热交换器的能量转换效率；ηHP为热泵利用电能生产热能的效率；Pe、Pg、Ph和PWT分别为上级电力、天然气、热力网络及风机的功率；λ1为分配系数，表示流经热电联产天然气占输入总量的比例；λ2为分配系数，表示流经电力电子转换器电能占输入总量的比例；
不等式约束为：
Figure FDA00031995711300000223
Figure FDA00031995711300000224
Pg_min≤Pg≤Pg_max
PHP_min≤PHP≤PHP_max
Le_min≤Le≤Le_max
Lh_min≤Lh≤Lh_max
其中，
Figure FDA0003199571130000031
为节点i的p相电压，
Figure FDA0003199571130000036
和
Figure FDA0003199571130000032
代表其上、下限；Sij为线路ij上所传输的视在功率，
Figure FDA0003199571130000033
代表其上限；
Figure FDA0003199571130000037
为管道中的水流量，
Figure FDA0003199571130000034
和
Figure FDA0003199571130000035
代表其上、下限；Pg_min和Pg_max分别为天然气输入功率Pg的下限及上限；PHP为热泵设备输出功率，PHP_min和PHP_max分别为其下限及上限；Le_max和Le_min分别为电负荷功率Le的上、下限，Lh_max和Lh_min分别为热负荷功率Lh的上、下限。
CN202110900365.XA 2021-08-06 2021-08-06 一种基于区域互联的能源市场双边竞价出清方法 Active CN113610394B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110900365.XA CN113610394B (zh) 2021-08-06 2021-08-06 一种基于区域互联的能源市场双边竞价出清方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110900365.XA CN113610394B (zh) 2021-08-06 2021-08-06 一种基于区域互联的能源市场双边竞价出清方法

Publications (2)

Publication Number Publication Date
CN113610394A true CN113610394A (zh) 2021-11-05
CN113610394B CN113610394B (zh) 2023-06-13

Family

ID=78339586

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110900365.XA Active CN113610394B (zh) 2021-08-06 2021-08-06 一种基于区域互联的能源市场双边竞价出清方法

Country Status (1)

Country Link
CN (1) CN113610394B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116934115A (zh) * 2023-07-18 2023-10-24 天津大学 一种考虑时变虚拟储能的实时端对端能量交易方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111784188A (zh) * 2020-07-14 2020-10-16 中国南方电网有限责任公司 基于电力市场的综合能源服务商市场出清机方法
CN111882111A (zh) * 2020-06-30 2020-11-03 华南理工大学 一种源网荷储协同互动的电力现货市场出清方法
CN112149914A (zh) * 2020-09-30 2020-12-29 四川大学 一种多约束条件下电力市场资源优化配置的方法
CN112381473A (zh) * 2021-01-14 2021-02-19 南方电网数字电网研究院有限公司 一种区域综合能源系统集群划分及参与日前市场出清方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111882111A (zh) * 2020-06-30 2020-11-03 华南理工大学 一种源网荷储协同互动的电力现货市场出清方法
CN111784188A (zh) * 2020-07-14 2020-10-16 中国南方电网有限责任公司 基于电力市场的综合能源服务商市场出清机方法
CN112149914A (zh) * 2020-09-30 2020-12-29 四川大学 一种多约束条件下电力市场资源优化配置的方法
CN112381473A (zh) * 2021-01-14 2021-02-19 南方电网数字电网研究院有限公司 一种区域综合能源系统集群划分及参与日前市场出清方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王丹 等: ""基于电-热联合市场出清的综合需求响应建模及策略"", 《电力系统自动化》, pages 1 - 3 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116934115A (zh) * 2023-07-18 2023-10-24 天津大学 一种考虑时变虚拟储能的实时端对端能量交易方法
CN116934115B (zh) * 2023-07-18 2024-04-09 天津大学 一种考虑时变虚拟储能的实时端对端能量交易方法

Also Published As

Publication number Publication date
CN113610394B (zh) 2023-06-13

Similar Documents

Publication Publication Date Title
CN111738502B (zh) 促进富余风电消纳的多能互补系统需求响应运行优化方法
Wu et al. Bi-level optimization model for integrated energy system considering the thermal comfort of heat customers
CN114362168B (zh) 一种能源互联系统的设备选型方法
Gu et al. Residential CCHP microgrid with load aggregator: Operation mode, pricing strategy, and optimal dispatch
Gabbar et al. Optimal scheduling of interconnected micro energy grids with multiple fuel options
CN113780776B (zh) 基于需求侧的电力系统碳运行调度方法、装置及设备
Chen et al. A robust optimization framework for energy management of CCHP users with integrated demand response in electricity market
CN102509175B (zh) 分布式供电系统可靠性优化方法
Jordehi Two-stage stochastic programming for risk-aware scheduling of energy hubs participating in day-ahead and real-time electricity markets
Tamashiro et al. Optimal components capacity based multi-objective optimization and optimal scheduling based MPC-optimization algorithm in smart apartment buildings
Zhi et al. Scenario-based multi-objective optimization strategy for rural PV-battery systems
Kou et al. A distributed energy management approach for residential demand response
CN114971372A (zh) 基于需求响应的微能网多时空尺度优化调度方法
Hayati et al. A two-stage stochastic optimization scheduling approach for integrating renewable energy sources and deferrable demand in the spinning reserve market
CN113610394B (zh) 一种基于区域互联的能源市场双边竞价出清方法
CN110661255B (zh) 一种多能源系统的热电优化运行方法、装置及设备
CN112365034A (zh) 一种电热综合能源系统调度方法及系统
Li et al. Aggregator’s scheduling and offering strategy for renewable integration based on information gap decision theory
Yu et al. Transactive control of a residential community with solar photovoltaic and battery storage systems
CN111523697A (zh) 一种综合能源服务成本分摊与定价计算方法
Denysiuk et al. Multiagent system for community energy management
CN112215687B (zh) 一种基于云边协同的家庭微网自主竞价系统及方法
CN115906456A (zh) 考虑需求侧响应不确定性的含氢储能ies调度优化模型
CN112257951B (zh) 一种基于合作博弈的综合能源系统与配电公司的优化运行方法
CN114282708A (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