CN111931977B - 一种考虑电—热能传输价值投入的虚拟电厂扩展规划模型构建方法 - Google Patents

一种考虑电—热能传输价值投入的虚拟电厂扩展规划模型构建方法 Download PDF

Info

Publication number
CN111931977B
CN111931977B CN202010582336.9A CN202010582336A CN111931977B CN 111931977 B CN111931977 B CN 111931977B CN 202010582336 A CN202010582336 A CN 202010582336A CN 111931977 B CN111931977 B CN 111931977B
Authority
CN
China
Prior art keywords
representing
equipment
load
node
layer
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
CN202010582336.9A
Other languages
English (en)
Other versions
CN111931977A (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.)
Southeast University
State Grid Liaoning Electric Power Co Ltd
Economic and Technological Research Institute of State Grid Jiangsu Electric Power Co Ltd
Original Assignee
Southeast University
State Grid Liaoning Electric Power Co Ltd
Economic and Technological Research Institute of State Grid Jiangsu 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 Southeast University, State Grid Liaoning Electric Power Co Ltd, Economic and Technological Research Institute of State Grid Jiangsu Electric Power Co Ltd filed Critical Southeast University
Priority to CN202010582336.9A priority Critical patent/CN111931977B/zh
Publication of CN111931977A publication Critical patent/CN111931977A/zh
Application granted granted Critical
Publication of CN111931977B publication Critical patent/CN111931977B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • 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

Landscapes

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

Abstract

本发明公开了一种考虑电—热能传输价值投入的虚拟电厂扩展规划模型构建方法,建立了包括上层设备定容层和下层设备选址层的虚拟电厂双层协调扩展规划模型,上层考虑虚拟电厂总投资运行投入最低,并同时考虑了规划方案对下层运行投入的影响,以设备建设容量、变压器及水泵的待建设状态、与外界交互的电能、热能量为决策变量,并将决策变量传递给下层;在下层的设备选址问题中,以产生的能源传输投入、燃料投入最低为目标,并根据上层的宏观规划方案,确定下层的决策变量,即各个设备的安装地址,并将下层的决策变量与目标函数统一返回至上层,如此反复迭代,直至达到上下层都满意的协调解,从而确定虚拟电厂的扩展规划方案。

Description

一种考虑电—热能传输价值投入的虚拟电厂扩展规划模型构 建方法
技术领域
本发明涉及电力系统规划领域,特别是涉及一种考虑电—热能传输价值投入的虚拟电厂扩展规划模型构建方法。
背景技术
近年来,分布式能源的广泛应用以及互联网技术的深度融合是能源互联网的主要特征之一。大量分布式资源的发展使得未来电网资源呈现出分布式小容量、资源种类多样、数量规模庞大、地域分散等特征,难以通过电网调度中心集中调度的方式对这些资源进行调控,将对电网的安全可靠运行造成不可忽视的影响。同时,各类分布式能源具有一定的调节能力,又为电网提供了可以利用的资源。基于云计算、物联网及先进的ICT技术,虚拟电厂是整合管理分布式资源的一种有效模式。虚拟电厂既可以应用在用户侧, 包括居民用户、工商业用户等不同类型,也可以应用在发电侧,包括冷热电联产、小型风场、小型水电站等。
目前,针对虚拟电厂的研究内容,多集中于优化调度与交易竞价方面,而在虚拟电厂优化配置方面所做的工作相对较少。面对虚拟电厂内部日益增长的负荷需求,虚拟电厂也面临扩展规划的必要,研究各类型分布式资源的运行特性,统筹考虑多种能源的需求,实现对各类资源的优化配置,既是构建多能虚拟电厂需要解决的关键技术问题之一,也可为其优化运行提供基础和指导。
目前在虚拟电厂的扩展规划研究领域已有一定的成果,但是,现有的关于虚拟电厂扩展规划的研究大多未考虑能源传输价值投入,事实上,由于虚拟电厂的聚合将占用一定的公共资源,能源传输价值投入是不可忽视的部分,在虚拟电厂规划时忽略能源传输价值投入将使得整体规划方案出现偏差。
发明内容
发明目的:为了应对虚拟电厂内部日益增加的能源需求,本发明考虑总投资运行投入,建立虚拟电厂扩展规划模型,提供了一种考虑电—热能传输价值投入的虚拟电厂扩展规划模型构建方法。
技术方案:一种考虑电—热能传输价值投入的虚拟电厂扩展规划模型构建方法,包括以下步骤:
步骤1、根据虚拟电厂区域内各个负荷节点的地理位置、以及用能情况,基于各节点初始负荷信息,预测各节点新增负荷量,从而确定总负荷需求,其中,负荷节点表示虚拟电厂区域内的用电用热终端,这些负荷节点处可以安装其它能源设备来供能;
步骤2、以总投资建设投入和运行投入最低为目标,建立上层设备定容层模型,所述上层设备定容层模型以设备建设容量、变压器及水泵的待建设状态、与外界交互的电能、热能量为决策变量;
步骤3、以总运行投入最低目标,建立下层设备选址层模型,所述下层设备选址层模型以需安装设备的负荷节点位置为决策变量;
步骤4、上下层通过各自的决策变量,采用交替迭代求解策略实现上下层模型的协调运行;
步骤5、确定虚拟电厂扩展规划方案,即虚拟电厂的设备建设容量、变压器及水泵的待建设状态、与外界交互的电能、热能量、需安装设备的负荷节点:
min C=Cinv1+Cinv2+Cmain+Cinter+Copr (1)
式(1)中,Cinv1表示CHP机组,燃气机组WT,燃气锅炉GB三种设备的投资建设投入,Cinv2表示变压器和水泵的投资建设投入,Cmain表示三种设备的年维修投入,Cinter表示与外界交互的电、热能量的投入,Copr表示下层目标函数;
Figure BDA0002552805110000021
式(2)中,设备序号1、2、3分别代表CHP机组,燃气机组WT,燃气锅炉GB;
Figure BDA0002552805110000022
表示第j种设备的建设容量;
Figure BDA0002552805110000023
表示第j种设备的单位容量建设投入;
Figure BDA0002552805110000024
表示第 j种设备建设投入的等年值系数;r表示折现率,Tj表示第j种设备的生命周期;
Figure BDA0002552805110000025
式(3)中,Costsub表示变压器的总建设投入;xsub表示变压器待建设状态的决策变量,为0-1变量,若该值为1,则代表新建变压器,若该值为0,则代表不需新建变压器;Costbump表示水泵的总建设投入;xbump表示水泵待建设状态的决策变量,若该值为1,则代表新建水泵,若该值为0,则代表不需新建水泵;r表示折现率,Ts和Tb分别表示变压器和水泵两种设备的生命周期;
Figure BDA0002552805110000031
式(4)中,
Figure BDA0002552805110000032
表示第j种设备的单位容量年维护投入;
Cinter=ωρ(kePgrid+khHheat) (5)
式(5)中,ke、kh分别表示与外界交互的电、热价值转换系数;Pgrid、Hheat分别表示与外界交互的电量与热量;ρ表示每年的平均运行小时数;ω表示平均负荷率,即一年内平均负荷功率与最大负荷功率的比值;
步骤2中所满足的约束条件如下:
(21)设备安装容量约束:
Figure BDA0002552805110000033
式(6)中,
Figure BDA0002552805110000034
分别表示第j种设备最小建设容量和最大建设容量;
(22)电负荷需求约束:
Figure BDA0002552805110000035
式(7)中,η1_e表示CHP机组的发电效率;P1 fix表示CHP机组的建设容量;η2表示燃气轮机的发电效率;
Figure BDA0002552805110000036
表示燃气轮机的建设容量;Pgrid表示与外界交互的电量;N为区域内负荷节点总数;Pi表示各个节点初始电负荷;Pnew,i表示各个节点新增的电负荷最大功率;
Figure BDA0002552805110000037
表示该区域内某一时刻下所有节点的最大电负荷需求,同时包括初始负荷与新增负荷;
(23)热负荷需求约束:
Figure BDA0002552805110000038
式(8)中,η2_e表示CHP机组的发热效率;P1 fix表示CHP机组的建设容量;η3表示燃气锅炉的发热效率;
Figure BDA0002552805110000041
表示燃气锅炉的建设容量;Hi表示各个节点初始热负荷;Hheat表示与外界交互的电量;Hnew,i表示各个节点新增的热负荷最大功率,
Figure BDA0002552805110000042
表示该区域内所有节点的最大热负荷需求,同时包括初始负荷与新增负荷;
(24)购电量约束:
Pgrid≤xsubPsub+Psub0 (9)
式(9)中,Pgrid表示与外界交互的电量;xsub表示变压器待建设状态的决策变量;Psub表示新增变压器的容量;Psub0表示原始变压器的容量;
(25)购热量约束:
HHeat≤xbumpHbump+Hbumo0 (10)
式(10)中,Hheat表示与外界交互的热量;xbump表示水泵待建设状态的决策变量;Hbump表示新增水泵由于水流量的限制导致的热功率最大传输量;Hbumo0表示原水泵最大可承受的热功率。
步骤3中,下层优化模型的目标函数为:
minCopr=ωρ(CP,trans+CH,trans+Cgas) (11)
式(11)中,CP,trans表示电能的传输价值投入;CH,trans表示热能传输价值投入;Cgas表示天然气的燃料投入;Copr表示年总运行投入;
Figure BDA0002552805110000043
式(12)中,CP,trans表示电能的传输价值投入;Nl表示该区域配电网的总支路数量;N为区域内负荷节点总数;l表示支路序号;PTDFPi l表示节点i处的电负荷对线路l的功率传输分布因子;由于在某个负荷节点安装设备后,该设备所发出的电功率将抵消一部分该节点的电负荷,因此,ΔPi表示在节点i处安装设备后该节点处的净电负荷;λP表示传输单位电能经过单位距离所消耗的传输价值,λP与PTDFPi l均为事先核定值,其中,节点表示电力网络和热力网络中的负荷节点;
Figure BDA0002552805110000051
式(13)中,η1_e表示CHP机组的发电效率;η2表示燃气轮机的发电效率;c1,i表示第1种设备,即CHP设备在节点i处的安装状态0-1变量,若该值为1,则代表CHP设备安装在了节点i处;若该值为0,则代表CHP设备没有安装在了节点i处;c2,i与之相同;
Figure BDA0002552805110000052
式(14)中,CH,trans表示热能传输价值投入;Np表示该区域热力管网的总管道数量,p表示管道序号;
Figure BDA0002552805110000053
表示节点i处的热负荷对管道p的流量传输分布因子;由于在某个负荷节点安装设备后,该设备所发出的热功率将抵消一部分该节点的热负荷,因此,ΔHi表示在节点i安装设备后该节点处的净热负荷;λH表示传输单位热能经过单位距离所消耗的传输价值,λH
Figure BDA0002552805110000054
均为事先核定值;
Figure BDA0002552805110000055
式(15)中,η2_e表示CHP机组的发热效率;η3表示燃气锅炉的发热效率;c3,i表示第3种设备,即燃气锅炉设备在节点i处的安装状态0-1变量,若该值为1,则代表燃气锅炉安装在节点i处,若该值为0,则代表燃气锅炉没有安装在节点i处;
Figure BDA0002552805110000056
式(16)中,kgas表示天然气价值转换系数;
步骤3中所满足的约束条件如下:
(31)设备安装约束:
Figure BDA0002552805110000057
该式表示第j种设备只能安装在一个节点处;
(32)支路功率传输约束:
Figure BDA0002552805110000058
式(18)中,Pl为线路l的最大传输容量;
(33)管道流量传输约束:
Figure BDA0002552805110000061
式(19)中,Hp为管道p的最大传输容量。
有益效果:与现有技术相比,本发明公开了一种考虑电—热能传输价值投入的虚拟电厂扩展规划模型构建方法,建立了包括上层设备定容层和下层设备选址层的虚拟电厂双层协调扩展规划模型,上层考虑的虚拟电厂总投资运行投入最低并考虑到该区域内原始负荷与新增负荷的总负荷需求,以设备建设状态、变压器及水泵的建设计划、外界购能量为决策变量,确定上层的宏观规划方案,并将决策变量传递给下层;在下层的设备选址问题中,以产生的能源传输价值投入、燃料投入最低为目标,并根据上层的宏观规划方案,确定下层的决策变量,即各个设备的安装地址,并将下层的决策变量与目标函数统一返回至上层,如此反复修正迭代,直至达到上下层都满意的协调解,从而确定虚拟电厂的最低扩展规划方案。
附图说明
图1是本发明方法流程图。
具体实施方式
下面结合实施例和说明书附图对本发明作进一步详细描述。
本发明的考虑电—热能传输价值投入的虚拟电厂扩展规划模型构建方法,包括以下步骤:
(1)针对虚拟电厂的扩展规划问题同时涉及设备的定容与选址两个层面的特性,建立了包括上层设备定容层和下层设备选址层的虚拟电厂双层协调扩展规划模型,上下两层决策通过各自的决策变量相互影响、相互作用。
(2)以总投资建设投入和运行投入最低为目标,建立上层设备定容层模型,所述上层优化模型以设备建设容量、变压器及水泵的待建设状态、与外界交互的电能、热能量为决策变量,其目标函数为:
min C=Cinv1+Cinv2+Cmain+Cinter+Copr (1)
式(1)中,Cinv1表示CHP机组,燃气机组WT,燃气锅炉GB三种设备的投资建设投入,Cinv2表示变压器和水泵的投资建设投入,Cmain表示三种设备的年维修投入,Cinter表示与外界交互的电、热能量的投入,Copr表示下层目标函数;
Figure BDA0002552805110000071
式(2)中,设备序号1、2、3分别代表CHP机组,燃气机组WT,燃气锅炉GB;
Figure BDA0002552805110000072
表示第j种设备的建设容量;
Figure BDA0002552805110000073
表示第j种设备的单位容量建设投入;
Figure BDA0002552805110000074
表示第 j种设备建设投入的等年值系数;r表示折现率,Tj表示第j种设备的生命周期;
Figure BDA0002552805110000075
式(3)中,Costsub表示变压器的总建设投入;xsub表示变压器待建设状态的决策变量,为0-1变量,若该值为1,则代表新建变压器,若该值为0,则代表不需新建变压器;Costbump表示水泵的总建设投入;xbump表示水泵待建设状态的决策变量,若该值为1,则代表新建水泵,若该值为0,则代表不需新建水泵;r表示折现率,Ts和Tb分别表示变压器和水泵两种设备的生命周期;
Figure BDA0002552805110000076
式(4)中,
Figure BDA0002552805110000077
表示第j种设备的单位容量年维护投入;
Cinter=ωρ(kePgrid+khHheat) (5)
式(5)中,ke、kh分别表示与外界交互的电、热价值转换系数;Pgrid、Hheat分别表示与外界交互的电量与热量;ρ表示每年的平均运行小时数;ω表示平均负荷率,即一年内平均负荷功率与最大负荷功率的比值;
步骤2中所满足的约束条件如下:
(21)设备安装容量约束:
Figure BDA0002552805110000078
式(6)中,
Figure BDA0002552805110000081
分别表示第j种设备最小建设容量和最大建设容量;
(22)电负荷需求约束:
Figure BDA0002552805110000082
式(7)中,η1_e表示CHP机组的发电效率;P1 fix表示CHP机组的建设容量;η2表示燃气轮机的发电效率;
Figure BDA0002552805110000083
表示燃气轮机的建设容量;Pgrid表示与外界交互的电量;N为区域内负荷节点总数;Pi表示各个节点初始电负荷;Pnew,i表示各个节点新增的电负荷最大功率;
Figure BDA0002552805110000084
表示该区域内某一时刻下所有节点的最大电负荷需求,同时包括初始负荷与新增负荷;
(23)热负荷需求约束:
Figure BDA0002552805110000085
式(8)中,η2_e表示CHP机组的发热效率;P1 fix表示CHP机组的建设容量;η3表示燃气锅炉的发热效率;
Figure BDA0002552805110000086
表示燃气锅炉的建设容量;Hi表示各个节点初始热负荷;Hheat表示与外界交互的电量;Hnew,i表示各个节点新增的热负荷最大功率,
Figure BDA0002552805110000087
表示该区域内所有节点的最大热负荷需求,同时包括初始负荷与新增负荷;
(24)购电量约束:
Pgrid≤xsubPsub+Psub0 (9)
式(9)中,Pgrid表示与外界交互的电量;xsub表示变压器待建设状态的决策变量;Psub表示新增变压器的容量;Psub0表示原始变压器的容量;
(25)购热量约束:
HHeat≤xbumpHbump+Hbumo0 (10)
式(10)中,Hheat表示与外界交互的热量;xbump表示水泵待建设状态的决策变量;Hbump表示新增水泵由于水流量的限制导致的热功率最大传输量;Hbumo0表示原水泵最大可承受的热功率。
步骤3中,下层优化模型的目标函数为:
minCopr=ωρ(CP,trans+CH,trans+Cgas) (11)
式(11)中,CP,trans表示电能的传输价值投入;CH,trans表示热能传输价值投入;Cgas表示天然气的燃料投入;Copr表示年总运行投入;
Figure BDA0002552805110000091
式(12)中,CP,trans表示电能的传输价值投入;Nl表示该区域配电网的总支路数量;N为区域内负荷节点总数;l表示支路序号;PTDFPi l表示节点i处的电负荷对线路l的功率传输分布因子;由于在某个负荷节点安装设备后,该设备所发出的电功率将抵消一部分该节点的电负荷,因此,ΔPi表示在节点i处安装设备后该节点处的净电负荷;λP表示传输单位电能经过单位距离所消耗的传输价值,λP与PTDFPi l均为事先核定值,其中,节点表示电力网络和热力网络中的负荷节点;
Figure BDA0002552805110000092
式(13)中,η1_e表示CHP机组的发电效率;η2表示燃气轮机的发电效率;c1,i表示第1种设备,即CHP设备在节点i处的安装状态0-1变量,若该值为1,则代表CHP设备安装在了节点i处;若该值为0,则代表CHP设备没有安装在了节点i处;c2,i与之相同;
Figure BDA0002552805110000093
式(14)中,CH,trans表示热能传输价值投入;Np表示该区域热力管网的总管道数量,p表示管道序号;
Figure BDA0002552805110000094
表示节点i处的热负荷对管道p的流量传输分布因子;由于在某个负荷节点安装设备后,该设备所发出的热功率将抵消一部分该节点的热负荷,因此,ΔHi表示在节点i安装设备后该节点处的净热负荷;λH表示传输单位热能经过单位距离所消耗的传输价值,λH
Figure BDA0002552805110000095
均为事先核定值;
Figure BDA0002552805110000101
式(15)中,η2_e表示CHP机组的发热效率;η3表示燃气锅炉的发热效率;c3,i表示第3种设备,即燃气锅炉设备在节点i处的安装状态0-1变量,若该值为1,则代表燃气锅炉安装在节点i处,若该值为0,则代表燃气锅炉没有安装在节点i处;
Figure BDA0002552805110000102
式(16)中,kgas表示天然气价值转换系数;
步骤3中所满足的约束条件如下:
(31)设备安装约束:
Figure BDA0002552805110000103
该式表示第j种设备只能安装在一个节点处;
(32)支路功率传输约束:
Figure BDA0002552805110000104
式(18)中,Pl为线路l的最大传输容量;
(33)管道流量传输约束:
Figure BDA0002552805110000105
式(19)中,Hp为管道p的最大传输容量。
(4)上下层通过各自的决策变量,采用交替迭代求解策略实现上下层模型的协调运行;上层首先根据自身的约束条件生产一系列满足约束条件的初始解,并输入下层;下层通过自身的求解,会将初始解对应的下层决策变量及目标函数返还给上层,上层接收到下层的返回值后,会根据自身的优化方法更新自身的决策变量;上下两层如此反复迭代,直到符合迭代终止条件或达到最大进化代数,最终获得双层规划模型的协调优化解。
(5)以步骤4中的优化解为参考,得到虚拟电厂扩展规划方案。
上述实施例仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和等同替换,这些对本发明权利要求进行改进和等同替换后的技术方案,均落入本发明的保护范围。

Claims (4)

1.一种考虑电—热能传输价值投入的虚拟电厂扩展规划模型构建方法,其特征在于,包括以下步骤:
步骤1、根据虚拟电厂区域内各个负荷节点的地理位置、以及用能情况,基于各节点初始负荷信息,预测各节点新增负荷量,从而确定总负荷需求,其中,负荷节点表示虚拟电厂区域内的用电用热终端,这些负荷节点处安装其它能源设备来供能;
步骤2、以总投资建设投入和运行投入最低为目标,建立上层设备定容层模型,所述上层设备定容层模型以设备建设容量、变压器及水泵的待建设状态、与外界交互的电能、热能量为决策变量;
步骤3、以总运行投入最低目标,建立下层设备选址层模型,所述下层设备选址层模型以需安装设备的负荷节点位置为决策变量;
步骤4、上下层通过各自的决策变量,采用交替迭代求解策略实现上下层模型的协调运行;
步骤5、确定虚拟电厂扩展规划方案,即虚拟电厂的设备建设容量、变压器及水泵的待建设状态、与外界交互的电能、热能量、需安装设备的负荷节点;
步骤2中,上层优化模型的目标函数为:
min C=Cinv1+Cinv2+Cmain+Cinter+Copr (1)
式(1)中,Cinv1表示CHP机组,燃气机组WT,燃气锅炉GB三种设备的投资建设投入,Cinv2表示变压器和水泵的投资建设投入,Cmain表示三种设备的年维修投入,Cinter表示与外界交互的电、热能量的投入,Copr表示下层目标函数;
Figure FDA0003679168260000011
式(2)中,设备序号1、2、3分别代表CHP机组,燃气机组WT,燃气锅炉GB;
Figure FDA0003679168260000012
表示第j种设备的建设容量;
Figure FDA0003679168260000013
表示第j种设备的单位容量建设投入;
Figure FDA0003679168260000014
表示第j种设备建设投入的等年值系数;r表示折现率,Tj表示第j种设备的生命周期;
Figure FDA0003679168260000015
式(3)中,Costsub表示变压器的总建设投入;xsub表示变压器待建设状态的决策变量,为0-1变量,若决策变量值为1,则代表新建变压器,若决策变量值为0,则代表不需新建变压器;Costbump表示水泵的总建设投入;xbump表示水泵待建设状态的决策变量,若决策变量值为1,则代表新建水泵,若决策变量值为0,则代表不需新建水泵;r表示折现率,Ts和Tb分别表示变压器和水泵两种设备的生命周期;
Figure FDA0003679168260000021
式(4)中,
Figure FDA0003679168260000022
表示第j种设备的单位容量年维护投入;
Cinter=ωρ(kePgrid+khHheat) (5)
式(5)中,ke、kh分别表示与外界交互的电、热价值转换系数;Pgrid、Hheat分别表示与外界交互的电量与热量;ρ表示每年的平均运行小时数;ω表示平均负荷率,即一年内平均负荷功率与最大负荷功率的比值;
步骤2中所满足的约束条件如下:
(21)设备安装容量约束:
Figure FDA0003679168260000023
式(6)中,
Figure FDA0003679168260000024
分别表示第j种设备最小建设容量和最大建设容量;
(22)电负荷需求约束:
Figure FDA0003679168260000025
式(7)中,η1_e表示CHP机组的发电效率;P1 fix表示CHP机组的建设容量;η2表示燃气轮机的发电效率;P2 fix表示燃气轮机的建设容量;Pgrid表示与外界交互的电量;N为区域内负荷节点总数;Pi表示各个节点初始电负荷;Pnew,i表示各个节点新增的电负荷最大功率;
Figure FDA0003679168260000026
表示该区域内某一时刻下所有节点的最大电负荷需求,同时包括初始负荷与新增负荷;
(23)热负荷需求约束:
Figure FDA0003679168260000031
式(8)中,η2_e表示CHP机组的发热效率;P1 fix表示CHP机组的建设容量;η3表示燃气锅炉的发热效率;P3 fix表示燃气锅炉的建设容量;Hi表示各个节点初始热负荷;Hheat表示与外界交互的电量;Hnew,i表示各个节点新增的热负荷最大功率,
Figure FDA0003679168260000032
表示该区域内所有节点的最大热负荷需求,同时包括初始负荷与新增负荷;
(24)购电量约束:
Pgrid≤xsubPsub+Psub0 (9)
式(9)中,Pgrid表示与外界交互的电量;xsub表示变压器待建设状态的决策变量;Psub表示新增变压器的容量;Psub0表示原始变压器的容量;
(25)购热量约束:
HHeat≤xbumpHbump+Hbumo0 (10)
式(10)中,Hheat表示与外界交互的热量;xbump表示水泵待建设状态的决策变量;Hbump表示新增水泵由于水流量的限制导致的热功率最大传输量;Hbumo0表示原水泵最大可承受的热功率。
2.根据权利要求1所述的一种考虑电—热能传输价值投入的虚拟电厂扩展规划模型构建方法,其特征在于:
步骤2中建立的所述上层设备定容层模型以总投资建设投入和运行投入最低为目标,且约束条件包括设备安装容量约束、与外界交互的电、热功率约束和电、热负荷需求约束;
步骤3中建立的所述下层设备选址层模型以总运行投入最低为目标,其中包含电、热传输价值投入、且约束条件包括设备安装数量约束、支路传输功率约束和管道传输容量约束。
3.根据权利要求1所述的一种考虑电—热能传输价值投入的虚拟电厂扩展规划模型构建方法,其特征在于,上层采用遗传算法求解,求解待建设设备的建设容量、变压器及水泵的待建设状态、与外界交互的电能、热能量,下层采用Cplex求解器求解,求解需要按照设备的负荷节点位置,上层首先产生一系列满足约束条件的初始个体,并输入下层;对于上层传入的每一个个体,下层均会将该个体对应的下层决策变量并返还给上层,从而计算考虑下层运行的个体适应度,并进行遗传算法的交叉、变异操作,更新得到新一代的个体,并重复上述操作;上下两层如此反复迭代,直到符合迭代终止条件或达到最大进化代数,最终获得最低的虚拟电厂扩展规划方案,即虚拟电厂的设备建设容量、变压器及水泵的待建设状态、与外界交互的电能、热能量、需安装设备的负荷节点。
4.根据权利要求1所述的一种考虑电—热能传输价值投入的虚拟电厂扩展规划模型构建方法,其特征在于,步骤3中,下层优化模型的目标函数为:
min Copr=ωρ(CP,trans+CH,trans+Cgas) (11)
式(11)中,CP,trans表示电能的传输价值投入;CH,trans表示热能传输价值投入;Cgas表示天然气的燃料投入;Copr表示年总运行投入;
Figure FDA0003679168260000041
式(12)中,CP,trans表示电能的传输价值投入;Nl表示该区域配电网的总支路数量;N为区域内负荷节点总数;l表示支路序号;PTDFPi l表示节点i处的电负荷对线路l的功率传输分布因子;由于在某个负荷节点安装设备后,该设备所发出的电功率将抵消一部分该节点的电负荷,因此,ΔPi表示在节点i处安装设备后该节点处的净电负荷;λP表示传输单位电能经过单位距离所消耗的传输价值,λP与PTDFPi l均为事先核定值,其中,节点表示电力网络和热力网络中的负荷节点;
ΔPi=Pi+Pnew,i-c1,iη1_eP1 fix-c2,iη2P2 fix (13)
式(13)中,η1_e表示CHP机组的发电效率;η2表示燃气轮机的发电效率;c1,i表示第1种设备,即CHP设备在节点i处的安装状态0-1变量,若该值为1,则代表CHP设备安装在了节点i处;若该值为0,则代表CHP设备没有安装在了节点i处;c2,i与之相同;
Figure FDA0003679168260000042
式(14)中,CH,trans表示热能传输价值投入;Np表示该区域热力管网的总管道数量,p表示管道序号;
Figure FDA0003679168260000043
表示节点i处的热负荷对管道p的流量传输分布因子;由于在某个负荷节点安装设备后,该设备所发出的热功率将抵消一部分该节点的热负荷,因此,ΔHi表示在节点i安装设备后该节点处的净热负荷;λH表示传输单位热能经过单位距离所消耗的传输价值,λH
Figure FDA0003679168260000051
均为事先核定值;
ΔHi=Hnew,i+Hi-c1,iη2_eP1 fix-c3,iη3P3 fix (15)
式(15)中,η2_e表示CHP机组的发热效率;η3表示燃气锅炉的发热效率;c3,i表示第3种设备,即燃气锅炉设备在节点i处的安装状态0-1变量,若该值为1,则代表燃气锅炉安装在节点i处,若该值为0,则代表燃气锅炉没有安装在节点i处;
Figure FDA0003679168260000052
式(16)中,kgas表示天然气价值转换系数;
步骤3中所满足的约束条件如下:
(31)设备安装约束:
Figure FDA0003679168260000053
该式表示第j种设备只能安装在一个节点处;
(32)支路功率传输约束:
Figure FDA0003679168260000054
式(18)中,Pl为线路l的最大传输容量;
(33)管道流量传输约束:
Figure FDA0003679168260000055
式(19)中,Hp为管道p的最大传输容量。
CN202010582336.9A 2020-06-23 2020-06-23 一种考虑电—热能传输价值投入的虚拟电厂扩展规划模型构建方法 Active CN111931977B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010582336.9A CN111931977B (zh) 2020-06-23 2020-06-23 一种考虑电—热能传输价值投入的虚拟电厂扩展规划模型构建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010582336.9A CN111931977B (zh) 2020-06-23 2020-06-23 一种考虑电—热能传输价值投入的虚拟电厂扩展规划模型构建方法

Publications (2)

Publication Number Publication Date
CN111931977A CN111931977A (zh) 2020-11-13
CN111931977B true CN111931977B (zh) 2022-11-18

Family

ID=73316985

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010582336.9A Active CN111931977B (zh) 2020-06-23 2020-06-23 一种考虑电—热能传输价值投入的虚拟电厂扩展规划模型构建方法

Country Status (1)

Country Link
CN (1) CN111931977B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112949940A (zh) * 2021-03-17 2021-06-11 东南大学 一种综合能源虚拟电厂多设备选址方法
CN117035202B (zh) * 2023-10-10 2024-01-23 国网山西省电力公司电力科学研究院 考虑需求响应的电热综合能源系统双层协同扩展规划方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108446809A (zh) * 2018-04-09 2018-08-24 国网河南省电力公司经济技术研究院 一种区域综合能源设备及网络双层优化配置方法
CN109063992A (zh) * 2018-07-18 2018-12-21 国网重庆市电力公司经济技术研究院 考虑区域综合能源系统优化运行的配电网扩展规划方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108446809A (zh) * 2018-04-09 2018-08-24 国网河南省电力公司经济技术研究院 一种区域综合能源设备及网络双层优化配置方法
CN109063992A (zh) * 2018-07-18 2018-12-21 国网重庆市电力公司经济技术研究院 考虑区域综合能源系统优化运行的配电网扩展规划方法

Also Published As

Publication number Publication date
CN111931977A (zh) 2020-11-13

Similar Documents

Publication Publication Date Title
Lv et al. Model predictive control based robust scheduling of community integrated energy system with operational flexibility
Zhang et al. Optimal operation of integrated electricity and heat system: A review of modeling and solution methods
Liu et al. Coordinated operation of multi-integrated energy system based on linear weighted sum and grasshopper optimization algorithm
Lin et al. Decentralized solution for combined heat and power dispatch through benders decomposition
Lu et al. Coordinated dispatch of multi-energy system with district heating network: Modeling and solution strategy
CN109858759A (zh) 一种工业园区综合能源平衡调度方法
CN110659830A (zh) 面向综合能源系统的多能源微网规划方法
CN109711601A (zh) 电-气-热综合能源系统分布式优化调度方法和装置
CN103903073A (zh) 一种含分布式电源及储能的微电网优化规划方法及系统
CN111931977B (zh) 一种考虑电—热能传输价值投入的虚拟电厂扩展规划模型构建方法
CN105955931A (zh) 面向高密度分布式光伏消纳的区域能源网络优化调度方法
Fan et al. Review on coordinated planning of source-network-load-storage for integrated energy systems
Wei et al. Distributed scheduling of smart buildings to smooth power fluctuations considering load rebound
CN112182887A (zh) 一种综合能源系统规划优化仿真方法
CN114595868A (zh) 一种综合能源系统源网荷储协同规划方法及系统
Liu et al. Two-stage optimal economic scheduling for commercial building multi-energy system through internet of things
CN116341881B (zh) 计及热网灵活性的电-热系统鲁棒超前调度方法及系统
Cheng et al. Energy management of CHP-based microgrid with thermal storage for reducing wind curtailment
Li et al. Intraday multi-objective hierarchical coordinated operation of a multi-energy system
JP2017200311A (ja) マイクログリッドの運転計画装置および方法、並びにマイクログリッドの運転計画装置で使用される地域エネルギー管理装置およびエネルギー管理装置
Hu et al. A decentralized market framework for procurement of operating reserves from district energy systems
Liu et al. Cooperative operation for integrated multi-energy system considering transmission losses
CN112214883B (zh) 一种微能源网综合能源系统运行优化方法和系统
An et al. Real-time optimal operation control of micro energy grid coupling with electricity-thermal-gas considering prosumer characteristics
CN113011655A (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