[go: up one dir, main page]

CN114648250A - An Integrated Energy System Planning Approach for Parks Considering Integrated Demand Response and Carbon Emissions - Google Patents

An Integrated Energy System Planning Approach for Parks Considering Integrated Demand Response and Carbon Emissions Download PDF

Info

Publication number
CN114648250A
CN114648250A CN202210378037.2A CN202210378037A CN114648250A CN 114648250 A CN114648250 A CN 114648250A CN 202210378037 A CN202210378037 A CN 202210378037A CN 114648250 A CN114648250 A CN 114648250A
Authority
CN
China
Prior art keywords
formula
energy
demand
carbon
period
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.)
Pending
Application number
CN202210378037.2A
Other languages
Chinese (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.)
Hefei University of Technology
Original Assignee
Hefei University of Technology
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 Hefei University of Technology filed Critical Hefei University of Technology
Priority to CN202210378037.2A priority Critical patent/CN114648250A/en
Publication of CN114648250A publication Critical patent/CN114648250A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; 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/06313—Resource planning in a project environment
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/12—Simultaneous equations, e.g. systems of linear equations
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/004—Artificial life, i.e. computing arrangements simulating life
    • G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; 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/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Mathematical Physics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Analysis (AREA)
  • Health & Medical Sciences (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Operations Research (AREA)
  • Pure & Applied Mathematics (AREA)
  • Marketing (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Software Systems (AREA)
  • Educational Administration (AREA)
  • Development Economics (AREA)
  • Quality & Reliability (AREA)
  • General Health & Medical Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Game Theory and Decision Science (AREA)
  • Water Supply & Treatment (AREA)
  • Biomedical Technology (AREA)
  • Public Health (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Primary Health Care (AREA)
  • Artificial Intelligence (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Evolutionary Computation (AREA)
  • Molecular Biology (AREA)

Abstract

The invention discloses a park comprehensive energy system planning method considering comprehensive demand response and carbon emission, which comprises the following steps: 1. establishing a stepped carbon emission mechanism model and considering the influence of the stepped carbon emission mechanism model on the system operation; 2. establishing a comprehensive demand response model based on the real-time electric heating and cooling demand of the user and a demand elasticity demand model; 3. considering system actual constraints, establishing a park comprehensive energy system day-ahead scheduling optimization model by taking low-carbon, low-energy-consumption and high-user satisfaction weighted scheduling sum as an optimization target, and performing coordinated optimization on supply and demand sides; 4. and constructing a park comprehensive energy system planning model considering comprehensive demand response and carbon emission, considering a scheduling strategy of comprehensive energy, and solving by adopting a particle swarm algorithm to obtain an optimal equipment configuration planning scheme of the park comprehensive energy system. The invention optimizes the planning scheme of the comprehensive energy system under the conditions of meeting the user requirements and being as low as possible in environment, thereby improving the safety and the high efficiency of the comprehensive energy system.

Description

计及综合需求响应和碳排放的园区综合能源系统规划方法An Integrated Energy System Planning Approach for Parks Considering Integrated Demand Response and Carbon Emissions

技术领域technical field

本发明属于综合能源系统规划领域，特别涉及一种计及综合需求响应和碳排放的园区规划方法；The invention belongs to the field of comprehensive energy system planning, in particular to a park planning method that takes into account comprehensive demand response and carbon emissions;

背景技术Background technique

随着能源技术进一步的发展，冷热电气等多种能源耦合的综合能源系统(integrated energy system，IES)得到深入研究和广泛应用，人类社会的进步与发展，带来了日益严峻的化石燃料短缺与环境污染问题，可再生能源代替化石能源得到积极推动，制定了碳排放控制承诺，传统的调度方法只对系统供给侧进行优化，不能调动系统需求侧的潜力，无法对碳排放进行控制，已有的规划方案不能很好地平衡系统用户的满意度和系统运行的安全性；With the further development of energy technology, integrated energy systems (IES) coupled with various energy sources such as cooling, heating and electricity have been deeply researched and widely used. The progress and development of human society have brought about an increasingly severe shortage of fossil fuels. Due to environmental pollution problems, renewable energy has been actively promoted to replace fossil energy, and carbon emission control commitments have been formulated. The traditional scheduling method only optimizes the supply side of the system, cannot mobilize the potential of the demand side of the system, and cannot control carbon emissions. Some planning schemes cannot well balance the satisfaction of system users and the security of system operation;

发明内容SUMMARY OF THE INVENTION

本发明是为了解决上述现有技术存在的不足之处，提出一种计及综合需求响应和碳排放的园区综合能源系统规划方法，以期能在满足用户用能需求的基础上降低系统能耗和碳排放，提高可再生能源消纳能力，从而能提高系统运行的安全性，得到最优的规划建设方案。In order to solve the above-mentioned shortcomings of the prior art, the present invention proposes a comprehensive energy system planning method for a park that takes into account comprehensive demand response and carbon emissions, in order to reduce system energy consumption and energy consumption on the basis of satisfying users' energy consumption needs. Carbon emissions, improve the absorption capacity of renewable energy, so as to improve the safety of system operation and obtain the optimal planning and construction scheme.

本发明为达到上述发明目的，采用如下技术方案：The present invention adopts the following technical scheme in order to achieve the above-mentioned purpose of the invention:

本发明一种计及综合需求响应和碳排放的园区综合能源系统规划方法，所述园区综合能源系统包括供给侧和需求侧，所述供给侧包括热电联产设备CHP、燃气轮机GB和P2G机组，所述需求侧包括电热泵EHP、中央空调AC、换热器HE和制冷机AF，其特点是，所述园区综合能源系统规划方法是按如下步骤进行：The present invention is a method for planning a comprehensive energy system in a park that takes into account comprehensive demand response and carbon emissions. The comprehensive energy system in the park includes a supply side and a demand side, and the supply side includes cogeneration equipment CHP, gas turbine GB, and P2G units. The demand side includes electric heat pump EHP, central air conditioner AC, heat exchanger HE and refrigerator AF, and is characterized in that the comprehensive energy system planning method of the park is carried out according to the following steps:

步骤一、通过基准值法计算系统的初始无偿碳排放量：Step 1. Calculate the initial free carbon emissions of the system through the benchmark value method:

步骤1.1、利用式(1)-式(4)计算无偿碳排放量：Step 1.1. Calculate free carbon emissions using equations (1)-(4):

Figure BDA0003590982500000011
Figure BDA0003590982500000011

Figure BDA0003590982500000012
Figure BDA0003590982500000012

Figure BDA0003590982500000013
Figure BDA0003590982500000013

Figure BDA0003590982500000014
Figure BDA0003590982500000014

式(1)-式(4)中，

Figure BDA0003590982500000015
表示向上级能源网络申领电量的无偿初始碳排放量；
Figure BDA0003590982500000016
表示热电联产设备CHP的无偿初始碳排放量；
Figure BDA0003590982500000017
表示燃气轮机组GB的无偿初始碳排放量；
Figure BDA0003590982500000018
表示单位电量的无偿初始碳排放量；
Figure BDA0003590982500000021
表示t时段内向上级能源网络申领电量；
Figure BDA0003590982500000022
表示单位热量的无偿初始碳排放量；
Figure BDA0003590982500000023
表示热电联产设备发电量e向发热量h的折算系数；
Figure BDA0003590982500000024
表示t时段内热电联产设备CHP用于产热消耗天然气能量；
Figure BDA0003590982500000025
表示t时段内热电联产设备CHP用于产电消耗天然气能量；
Figure BDA00035909825000000214
表示t时段内燃气轮机组设备GB的消耗天然气能量；T表示调度时间的集合；In formula (1) - formula (4),
Figure BDA0003590982500000015
Represents the free initial carbon emissions from applying for electricity from the upper-level energy network;
Figure BDA0003590982500000016
Represents the unpaid initial carbon emissions of CHP equipment;
Figure BDA0003590982500000017
Represents the free initial carbon emissions of the gas turbine unit GB;
Figure BDA0003590982500000018
Represents the free initial carbon emissions per unit of electricity;
Figure BDA0003590982500000021
Represents the application for electricity from the upper-level energy network within the t period;
Figure BDA0003590982500000022
Represents the unpaid initial carbon emissions per unit of heat;
Figure BDA0003590982500000023
Represents the conversion coefficient from the power generation e to the calorific value h of the cogeneration equipment;
Figure BDA0003590982500000024
Indicates that the combined heat and power equipment CHP is used to generate heat and consume natural gas energy during the t period;
Figure BDA0003590982500000025
Indicates that the combined heat and power equipment CHP is used to generate electricity and consume natural gas energy during the t period;
Figure BDA00035909825000000214
Represents the natural gas energy consumption of the gas turbine unit equipment GB in the t period; T represents the set of dispatching time;

步骤1.2、利用式(5)-式(6)计算实际碳排放：Step 1.2. Calculate the actual carbon emissions using equations (5)-(6):

M＝Melcbuy+MCHP+MGB-MP2G (5)M=M elcbuy +M CHP +M GB -M P2G (5)

Figure BDA0003590982500000026
Figure BDA0003590982500000026

式(5)-式(6)中，Melcbuy表示向上级能源网络申领电量碳排放，MCHP表示热电联产设备CHP的碳排量，MP2G表示P2G机组的消耗碳量；βg表示P2G机组的碳捕获系数；

Figure BDA00035909825000000215
表示t时段内P2G设备消耗的电量；In equations (5)-(6), M elcbuy represents the carbon emission of electricity applied to the upper energy network, M CHP represents the carbon emission of CHP of cogeneration equipment, M P2G represents the carbon consumption of P2G units; β g represents Carbon capture factor of P2G units;
Figure BDA00035909825000000215
Indicates the power consumed by the P2G device in the t period;

步骤1.3、利用式(7)计算阶梯碳排放量Ccarbon：Step 1.3, use formula (7) to calculate the step carbon emission C carbon :

Figure BDA0003590982500000027
Figure BDA0003590982500000027

式(7)中，c表示碳排放基准系数；α为碳排放系数的增长率；d为阶梯区间的长度；In formula (7), c is the carbon emission benchmark coefficient; α is the growth rate of the carbon emission coefficient; d is the length of the step interval;

步骤二、基于用户实时电热冷需求量和需求弹性理论，建立综合需求响应模型：Step 2. Based on the user's real-time electric heating and cooling demand and demand elasticity theory, establish a comprehensive demand response model:

步骤2.1、利用式(8)-式(10)计算冷热电替代系数：Step 2.1, use formula (8)-formula (10) to calculate the replacement coefficient of cold, heat and electricity:

Figure BDA0003590982500000028
Figure BDA0003590982500000028

Figure BDA0003590982500000029
Figure BDA0003590982500000029

Figure BDA00035909825000000210
Figure BDA00035909825000000210

式(8)-式(10)中，

Figure BDA00035909825000000211
和
Figure BDA00035909825000000212
是t时段第i类用户的电、热和冷直接能源需求；
Figure BDA00035909825000000213
和
Figure BDA0003590982500000031
表示第i类用户在t时段电冷负荷、电热负荷的用能偏好；
Figure BDA0003590982500000032
分别表示热电替代系数和冷电替代系数，并由
Figure BDA0003590982500000033
和
Figure BDA0003590982500000034
进行归一化后得到；δi表示第i类型用户非刚性电负荷的比例；In formula (8) - formula (10),
Figure BDA00035909825000000211
and
Figure BDA00035909825000000212
is the direct energy demand for electricity, heat and cooling of the i-th user in period t;
Figure BDA00035909825000000213
and
Figure BDA0003590982500000031
Represents the energy use preference of the i-th user in the electric cooling load and electric heating load in the t period;
Figure BDA0003590982500000032
respectively represent the thermoelectric substitution coefficient and the cold electric substitution coefficient, and are given by
Figure BDA0003590982500000033
and
Figure BDA0003590982500000034
Obtained after normalization; δ i represents the proportion of the non-rigid electrical load of the i-th type of user;

步骤2.2、利用式(11)-式(12)计算综合需求响应量：Step 2.2, use formula (11) - formula (12) to calculate the comprehensive demand response:

Figure BDA0003590982500000035
Figure BDA0003590982500000035

Figure BDA0003590982500000036
Figure BDA0003590982500000036

式(11)-式(12)中，

Figure BDA0003590982500000037
表示第i类型用户t时段的原始电负荷；εtt′表示t与t′时段之间的需求互弹性系数；当t＝t′时，εtt′表示需求自弹性系数；
Figure BDA0003590982500000038
表示t′时段的基础电需求惩罚；
Figure BDA0003590982500000039
表示t′时段的电需求惩罚变化量；
Figure BDA00035909825000000310
和
Figure BDA00035909825000000311
表示第i类型用户在t时段的电负荷、热负荷和冷负荷响应量；In formula (11) - formula (12),
Figure BDA0003590982500000037
Represents the original electrical load of the i-th type of user in period t; ε tt' represents the demand mutual elasticity coefficient between t and t'period; when t=t', ε tt' represents the demand self-elasticity coefficient;
Figure BDA0003590982500000038
represents the basic electricity demand penalty in the t'period;
Figure BDA0003590982500000039
Represents the amount of electricity demand penalty change in the t'period;
Figure BDA00035909825000000310
and
Figure BDA00035909825000000311
Represents the electrical load, heating load and cooling load response of the i-th type of user in the t period;

步骤2.3、利用式(13)计算舒适度补偿量：Step 2.3, use formula (13) to calculate the amount of comfort compensation:

Figure BDA00035909825000000312
Figure BDA00035909825000000312

式(13)中，Ccomf表示舒适度响应补偿量，λcomf表示单位能量的舒适度补偿量；In formula (13), C comf represents the comfort response compensation amount, and λ comf represents the comfort level compensation amount per unit energy;

步骤三、建立园区综合能源系统的日前调度优化模型，考虑实际约束且以低碳低能耗高用户满意度加权调和数为优化目标函数，对供给和需求侧进行协调优化：Step 3: Establish a day-ahead scheduling optimization model for the integrated energy system of the park, consider the actual constraints and take the weighted harmonic number of low carbon, low energy consumption and high user satisfaction as the optimization objective function to coordinate and optimize the supply and demand sides:

步骤3.1、利用式(14)-式(16)计算系统低碳低能耗高用户满意度加权调和数：Step 3.1. Use equations (14)-(16) to calculate the weighted harmonic number of low carbon, low energy consumption and high user satisfaction of the system:

min C＝Cbuy+Ccarbon+Cpena+Ccomf (14)min C=C buy +C carbon +C pena +C comf (14)

Figure BDA00035909825000000313
Figure BDA00035909825000000313

Figure BDA00035909825000000314
Figure BDA00035909825000000314

式(14)-式(16)中，Cbuy是向上级能源网络申领总能量；Cpena是弃风光惩罚；

Figure BDA00035909825000000317
和
Figure BDA00035909825000000318
分别表示t时段向上级能源网络申领的电量、申领的天然气量和申领的热量；cpena是弃风光惩罚系数；
Figure BDA00035909825000000315
和
Figure BDA00035909825000000316
是风、光预测功率；Ppv,t和Pwt,t是系统t时段实际消纳的光电功率与风电功率；In Equation (14)-Equation (16), C buy is the total energy applied to the upper energy network; C pena is the penalty for abandoning the scenery;
Figure BDA00035909825000000317
and
Figure BDA00035909825000000318
Represents the amount of electricity applied for, the amount of natural gas applied for, and the amount of heat applied to the upper energy network in t period; c pena is the penalty coefficient for abandoning wind and solar;
Figure BDA00035909825000000315
and
Figure BDA00035909825000000316
are the predicted power of wind and light; P pv,t and P wt,t are the photovoltaic power and wind power actually consumed by the system during period t;

步骤3.2、利用式(17)-式(19)定义能源供给侧约束：Step 3.2. Use equations (17)-(19) to define energy supply side constraints:

Figure BDA0003590982500000041
Figure BDA0003590982500000041

Figure BDA0003590982500000042
Figure BDA0003590982500000042

Figure BDA0003590982500000043
Figure BDA0003590982500000043

式(17)-式(19)中，

Figure BDA0003590982500000044
和
Figure BDA0003590982500000045
表示系统t时段P2G机组的出力、热电联产设备CHP的电出力、热电联产设备的CHP热出力和燃气锅炉GB的出力；η表示设备效率；In formula (17) - formula (19),
Figure BDA0003590982500000044
and
Figure BDA0003590982500000045
Represents the output of P2G units, the electrical output of CHP of co-generation equipment, the CHP thermal output of co-generation equipment and the output of gas-fired boiler GB in the system t period; η represents the equipment efficiency;

步骤3.2、利用式(20)-式(24)定义能源需求侧约束：Step 3.2, using equations (20)-(24) to define energy demand side constraints:

Figure BDA00035909825000000415
Figure BDA00035909825000000415

Figure BDA00035909825000000416
Figure BDA00035909825000000416

Figure BDA0003590982500000046
Figure BDA0003590982500000046

Figure BDA0003590982500000047
Figure BDA0003590982500000047

Figure BDA0003590982500000048
Figure BDA0003590982500000048

式(20)-式(24)中，

Figure BDA00035909825000000417
和
Figure BDA00035909825000000418
表示系统t时段电热泵EHP出力、中央空调AC出力、换热器HE出力和制冷机AF出力；
Figure BDA00035909825000000419
和
Figure BDA00035909825000000420
表示响应前系统t时段电负荷需求、热负荷需求和冷负荷需求；In formula (20) - formula (24),
Figure BDA00035909825000000417
and
Figure BDA00035909825000000418
Indicates the EHP output of the electric heat pump, the AC output of the central air conditioner, the HE output of the heat exchanger and the AF output of the refrigerator in the system t period;
Figure BDA00035909825000000419
and
Figure BDA00035909825000000420
Represents the electrical load demand, heating load demand and cooling load demand of the system before the response period t;

步骤3.3、利用式(25)-式(31)定义设备运行约束：Step 3.3. Use equations (25)-(31) to define equipment operating constraints:

Figure BDA0003590982500000049
Figure BDA0003590982500000049

Figure BDA00035909825000000410
Figure BDA00035909825000000410

Figure BDA00035909825000000411
Figure BDA00035909825000000411

Figure BDA00035909825000000412
Figure BDA00035909825000000412

Figure BDA00035909825000000413
Figure BDA00035909825000000413

Figure BDA00035909825000000414
Figure BDA00035909825000000414

Figure BDA0003590982500000051
Figure BDA0003590982500000051

步骤3.4、利用式(32)-式(33)定义风光出力约束：Step 3.4, use Equation (32)-Equation (33) to define the wind and light output constraints:

Figure BDA0003590982500000052
Figure BDA0003590982500000052

Figure BDA0003590982500000053
Figure BDA0003590982500000053

步骤3.5、利用式(34)定义电需求惩罚变化量约束：Step 3.5, use formula (34) to define the electricity demand penalty variation constraint:

Figure BDA0003590982500000054
Figure BDA0003590982500000054

步骤四、构建计及综合需求响应和碳排放的园区综合能源系统规划模型：Step 4. Build a comprehensive energy system planning model for the park that takes into account comprehensive demand response and carbon emissions:

步骤4.1、利用式(35)-式(37)构建园区综合能源系统规划模型的目标函数：Step 4.1. Use equations (35)-(37) to construct the objective function of the comprehensive energy system planning model of the park:

Figure BDA0003590982500000055
Figure BDA0003590982500000055

Figure BDA0003590982500000056
Figure BDA0003590982500000056

Figure BDA0003590982500000057
Figure BDA0003590982500000057

式(35)-式(37)中，ΩY为规划年集合；ΩD为规划的设备集合；

Figure BDA0003590982500000058
为第d个设备的类型集合；ΩS为季度的集合；Cinvest为投资耗材；
Figure BDA0003590982500000059
为第y年的运行耗材；
Figure BDA00035909825000000510
为第d个设备第c种类型的投资耗材；xc,d为布尔变量，表示是否投资第d类设备第c种类型；
Figure BDA00035909825000000511
为第y年第s季度典型日的低碳低能耗高用户满意度加权调和数；ns为第s季度的天数；ρ为损耗率；In equations (35)-(37), Ω Y is the set of planning years; Ω D is the set of planned equipment;
Figure BDA0003590982500000058
is the type collection of the d-th equipment; Ω S is the collection of quarters; C invest is the investment consumables;
Figure BDA0003590982500000059
Consumables for the yth year;
Figure BDA00035909825000000510
is the investment consumables of the c-th type of the d-th equipment; x c, d is a Boolean variable, indicating whether to invest in the c-th type of the d-th equipment;
Figure BDA00035909825000000511
is the weighted harmonic number of low carbon, low energy consumption and high user satisfaction on typical days in the s quarter of the yth year; n s is the number of days in the s quarter; ρ is the loss rate;

步骤4.2、利用式(38)定义设备的投资类型约束：Step 4.2, use formula (38) to define the investment type constraint of equipment:

Figure BDA00035909825000000512
Figure BDA00035909825000000512

式(38)表示对于任意d类设备，投资的类型不超过1种；Equation (38) indicates that for any type d equipment, there is no more than one type of investment;

步骤4.3、利用式(39)计算第y年第s季度负荷预测值Pj,y,s：Step 4.3, use formula (39) to calculate the load forecast value P j,y,s in the s quarter of the yth year:

Pj,y,s＝(1+γ)tPj,0,s,j∈ΩB,y∈ΩY,s∈ΩS (39)P j,y,s =(1+γ) t P j,0,s ,j∈Ω B ,y∈Ω Y ,s∈Ω S (39)

式(39)中，γ为负荷年增长率；Pi,0,s为当前年第s季度节点j的负荷值；In formula (39), γ is the annual growth rate of load; P i,0,s is the load value of node j in the sth quarter of the current year;

步骤五、采用粒子群算法求解园区综合能源系统规划模型：Step 5. Use particle swarm algorithm to solve the planning model of the comprehensive energy system of the park:

步骤5.1、输入初始参数，包括：粒子群种群规模M、学习因子c1和c2、惯性权重w、粒子群繁殖代数Mc、蒙特卡罗模拟次数Ms、置信区间β；Step 5.1. Input initial parameters, including: particle swarm population size M, learning factors c 1 and c 2 , inertia weight w, particle swarm reproduction algebra M c , Monte Carlo simulation times M s , confidence interval β;

步骤5.2、随机生成M个初始粒子并构成粒子集合M＝{m1,m2,…,mk,…,mM}，其中，mk额第k个粒子，表示从ΩD中选择不同的设备构成的第k个规划方案，且mk＝{mk1,mk2,…,mkd,…mkD}，其中，mkd表示第k个粒子选的d类设备的容量；Step 5.2. Randomly generate M initial particles and form a particle set M={m 1 , m 2 ,...,m k ,...,m M }, where m k is the kth particle, which means selecting different particles from Ω D The kth planning scheme composed of the equipment of , and m k ={m k1 ,m k2 ,...,m kd ,... m kD }, where m kd represents the capacity of the type d equipment selected by the kth particle;

步骤5.3、根据每个粒子对应的设备容量规划方案，计算投资耗材；Step 5.3. Calculate the investment consumables according to the equipment capacity planning scheme corresponding to each particle;

步骤5.4、根据式(39)，更新第y年的负荷需求，并计算第y年第s季度典型日的运行耗材；计算规划年内总的运行耗材，将计算的投资耗材与低碳低能耗高用户满意度加权调和数之和作为每个粒子的适应度；Step 5.4. According to formula (39), update the load demand in the yth year, and calculate the operating consumables for the typical day in the yth year and the s quarter; calculate the total operating consumables in the planning year, and compare the calculated investment consumables with low-carbon, low-energy consumption and high energy consumption. The sum of user satisfaction weighted harmonic numbers is used as the fitness of each particle;

步骤5.5、更新粒子位置和速度，从而获得新的粒子；Step 5.5, update the particle position and velocity to obtain new particles;

步骤5.6、重复步骤5.3-步骤5.5，直至达到给定的粒子群繁殖代数Mc为止；Step 5.6, repeat steps 5.3-5.5 until the given particle swarm reproduction algebra M c is reached;

步骤5.7、将最好的粒子所对应的设备容量作为园区综合能源系统的最优规划方案。Step 5.7, take the equipment capacity corresponding to the best particles as the optimal planning scheme of the comprehensive energy system of the park.

与现有技术相比，本发明的有益效果在于：Compared with the prior art, the beneficial effects of the present invention are:

本发明考虑了碳排放和综合需求响应在满足用户需求和尽可能低碳环保的情况下降低系统运营商的运行耗材，提供了综合能源系统规划方案，从而控制了系统碳排放，增强了系统可再生能源消纳能力，进而提高了综合能源系统的可靠性和高效性。The present invention considers carbon emissions and comprehensive demand response to reduce the operating consumables of system operators under the condition of meeting user needs and being as low-carbon and environmentally friendly as possible, and provides a comprehensive energy system planning scheme, thereby controlling system carbon emissions and enhancing system availability. The ability to absorb renewable energy, thereby improving the reliability and efficiency of the integrated energy system.

附图说明Description of drawings

图1为园区综合能源系统架构图；Figure 1 is the architecture diagram of the integrated energy system of the park;

图2为阶梯碳排放机制示意图；Figure 2 is a schematic diagram of the stepped carbon emission mechanism;

图3为本发明方法的流程图。Figure 3 is a flow chart of the method of the present invention.

具体实施方式Detailed ways

本实施例中，一种计及综合需求响应和碳排放的园区综合能源系统规划方法应用于如图1所示的园区综合能源系统，该园区综合能源系统包括供给侧和需求侧，供给侧包括热电联产设备CHP、燃气轮机GB和P2G机组，需求侧包括电热泵EHP、中央空调AC、换热器HE和制冷机AF，该园区综合能源系统规划方法的主要步骤包括：In this embodiment, a comprehensive energy system planning method for the park that takes into account comprehensive demand response and carbon emissions is applied to the comprehensive energy system of the park as shown in FIG. 1 . The comprehensive energy system of the park includes a supply side and a demand side, and the supply side includes Combined heat and power equipment CHP, gas turbine GB and P2G units, the demand side includes electric heat pump EHP, central air conditioner AC, heat exchanger HE and refrigerator AF, the main steps of the comprehensive energy system planning method of the park include:

1)通过综合需求响应充分调动用户侧的优化潜力，在运行模型目标函数中按照图2所示考虑碳排放量，弃风光控制系统碳排放和增强可再生能源消纳能力，在运行模型的基础上建立规划模型，提出的规划方案能够完全适应低碳低能耗高用户满意度的运行方案；1) Fully mobilize the optimization potential of the user side through comprehensive demand response, consider carbon emissions in the objective function of the operating model as shown in Figure 2, abandon the carbon emissions of the wind-solar control system and enhance the ability to absorb renewable energy, on the basis of the operating model The planning model is established on the above, and the proposed planning scheme can fully adapt to the operation scheme of low carbon, low energy consumption and high user satisfaction;

2)基于系统用户的用能历史数据，得到不同类型用户在不同时段的冷热电替代系数，从而计算综合需求响应量，能够兼顾系统用户的满意度；2) Based on the historical energy consumption data of system users, the replacement coefficients of cooling, heating and electricity for different types of users in different time periods are obtained, so as to calculate the comprehensive demand response amount, which can take into account the satisfaction of system users;

3)建立了园区综合能源系统日前调度优化模型，目标函数以最小化系统低碳低能耗高用户满意度加权调和数，包括向上级能源网络申领量，碳排放量，弃风光惩罚，综合需求响应舒适度补偿量；约束条件包括综合能源平衡约束，设备运行约束，风光出力约束；3) Established a day-ahead scheduling optimization model for the integrated energy system of the park. The objective function is to minimize the system’s low-carbon, low-energy consumption, and high-user satisfaction weighted reconciliation number, including the amount of claims to the upper energy network, carbon emissions, penalty for abandoning wind and solar, and comprehensive demand. Response comfort compensation; constraints include comprehensive energy balance constraints, equipment operation constraints, and wind and solar output constraints;

4)按照图3所示采用粒子群算法求解设备容量规划模型，针对每次抽样结果调用园区综合能源系统日前调度优化模型，求出规划和运行耗材，以最好的粒子作为优化问题的最优规划方案；具体的说，是按如下步骤进行：4) As shown in Figure 3, the particle swarm algorithm is used to solve the equipment capacity planning model, and the day-ahead scheduling optimization model of the integrated energy system of the park is called for each sampling result to find the planning and operation consumables, and the best particles are used as the optimal solution for the optimization problem. The planning scheme; specifically, it is carried out according to the following steps:

步骤一、系统的初始无偿碳排放量通过基准值法计算得到：Step 1. The initial free carbon emissions of the system are calculated by the benchmark value method:

步骤1.1、利用式(1)-式(4)计算无偿碳排放量：Step 1.1. Calculate free carbon emissions using equations (1)-(4):

Figure BDA0003590982500000071
Figure BDA0003590982500000071

Figure BDA0003590982500000072
Figure BDA0003590982500000072

Figure BDA0003590982500000073
Figure BDA0003590982500000073

Figure BDA0003590982500000074
Figure BDA0003590982500000074

式(1)表示系统的总初始无偿碳排放量；式(2)、式(3)和式(4)分别表示向上级能源网络申领电量的无偿碳排放量、热电联产设备CHP的无偿碳排放量和燃气轮机GB的无偿碳排放量；Equation (1) represents the total initial gratuitous carbon emissions of the system; Equations (2), (3) and (4) represent the gratuitous carbon emissions from applying for electricity from the upper energy network, and the gratuitous carbon emissions from the combined heat and power equipment CHP, respectively. Carbon emissions and unpaid carbon emissions from gas turbine GB;

式(1)-式(4)中，

Figure BDA0003590982500000075
表示向上级能源网络申领电量的无偿初始碳排放量；
Figure BDA0003590982500000076
表示热电联产设备CHP的无偿初始碳排放量；
Figure BDA0003590982500000077
表示燃气轮机组GB的无偿初始碳排放量；
Figure BDA0003590982500000078
表示单位电量的无偿初始碳排放量；
Figure BDA00035909825000000714
表示t时段内向上级能源网络申领电量；
Figure BDA0003590982500000079
表示单位热量的无偿初始碳排放量；
Figure BDA00035909825000000710
表示热电联产设备CHP发电量向发热量的折算系数；
Figure BDA00035909825000000711
表示t时段内热电联产设备CHP用于产热消耗天然气能量；
Figure BDA00035909825000000712
表示t时段内热电联产设备CHP用于产电消耗天然气能量；
Figure BDA00035909825000000715
表示t时段内燃气轮机组GB消耗天然气能量；In formula (1) - formula (4),
Figure BDA0003590982500000075
Represents the free initial carbon emissions from applying for electricity from the upper-level energy network;
Figure BDA0003590982500000076
Represents the unpaid initial carbon emissions of CHP equipment;
Figure BDA0003590982500000077
Represents the free initial carbon emissions of the gas turbine unit GB;
Figure BDA0003590982500000078
Represents the free initial carbon emissions per unit of electricity;
Figure BDA00035909825000000714
Represents the application for electricity from the upper-level energy network within the t period;
Figure BDA0003590982500000079
Represents the unpaid initial carbon emissions per unit of heat;
Figure BDA00035909825000000710
Represents the conversion coefficient from CHP power generation to calorific value of cogeneration equipment;
Figure BDA00035909825000000711
Indicates that the combined heat and power equipment CHP is used to generate heat and consume natural gas energy during the t period;
Figure BDA00035909825000000712
Indicates that the combined heat and power equipment CHP is used to generate electricity and consume natural gas energy during the t period;
Figure BDA00035909825000000715
Represents the natural gas energy consumed by the gas turbine unit GB in the t period;

步骤1.2、利用式(5)-式(6)计算实际碳排放：Step 1.2. Calculate the actual carbon emissions using equations (5)-(6):

M＝Melcbuy+MCHP+MGB-MP2G (5)M=M elcbuy +M CHP +M GB -M P2G (5)

Figure BDA00035909825000000713
Figure BDA00035909825000000713

式(5)表示系统实际碳排放；式(6)表示P2G机组消纳碳量；Equation (5) represents the actual carbon emission of the system; Equation (6) represents the amount of carbon consumed by P2G units;

式(5)-式(6)计算计算系统实际碳排量，Melcbuy表示向上级能源网络申领电量碳排放，MCHP表示热电联产设备CHP碳排量，MP2G表示P2G机组消耗碳量；βg表示P2G机组的碳捕获系数；

Figure BDA00035909825000000814
表示t时段内P2G机组消耗的电量；Equations (5)-(6) are used to calculate the actual carbon emissions of the calculation system, M elcbuy represents the carbon emissions of electricity applied to the upper energy network, M CHP represents the CHP carbon emissions of cogeneration equipment, and M P2G represents the carbon consumption of P2G units ; β g represents the carbon capture coefficient of the P2G unit;
Figure BDA00035909825000000814
Indicates the power consumed by the P2G unit in the t period;

步骤1.3、利用式(7)计算阶梯碳排放量：Step 1.3, use formula (7) to calculate the carbon emissions of the ladder:

Figure BDA0003590982500000081
Figure BDA0003590982500000081

式(7)中c表示碳排放基准系数；α为碳排放系数的增长率；d为阶梯区间的长度；In formula (7), c represents the carbon emission benchmark coefficient; α is the growth rate of the carbon emission coefficient; d is the length of the step interval;

步骤二、基于用户实时电热冷需求量和需求弹性理论，建立综合需求响应模型：Step 2. Based on the user's real-time electric heating and cooling demand and demand elasticity theory, establish a comprehensive demand response model:

步骤2.1、利用式(8)-式(10)计算冷热电替代系数：Step 2.1, use formula (8)-formula (10) to calculate the replacement coefficient of cold, heat and electricity:

Figure BDA0003590982500000082
Figure BDA0003590982500000082

Figure BDA0003590982500000083
Figure BDA0003590982500000083

Figure BDA0003590982500000084
Figure BDA0003590982500000084

式(8)和式(9)是电冷热替代系数的计算方法；式(10)是电冷热替代系数的约束；Equation (8) and Equation (9) are the calculation methods of the electric cooling and heating substitution coefficient; Equation (10) is the constraint of the electric cooling and heating substitution coefficient;

式(8)-式(10)中，

Figure BDA0003590982500000085
和
Figure BDA0003590982500000086
是t时段第i类用户的电、热和冷直接能源需求；
Figure BDA0003590982500000087
和
Figure BDA0003590982500000088
表示第i类用户在t时段电冷负荷、电热负荷的用能偏好；
Figure BDA0003590982500000089
分别表示热电替代系数和冷电替代系数，并由
Figure BDA00035909825000000810
和
Figure BDA00035909825000000811
进行归一化后得到；δi表示第i类型用户非刚性电负荷的比例；In formula (8) - formula (10),
Figure BDA0003590982500000085
and
Figure BDA0003590982500000086
is the direct energy demand for electricity, heat and cooling of the i-th user in period t;
Figure BDA0003590982500000087
and
Figure BDA0003590982500000088
Represents the energy use preference of the i-th user in the electric cooling load and electric heating load in the t period;
Figure BDA0003590982500000089
respectively represent the thermoelectric substitution coefficient and the cold electric substitution coefficient, and are given by
Figure BDA00035909825000000810
and
Figure BDA00035909825000000811
Obtained after normalization; δ i represents the proportion of the non-rigid electrical load of the i-th type of user;

步骤2.2、利用式(11)-式(12)计算综合需求响应量：Step 2.2, use formula (11) - formula (12) to calculate the comprehensive demand response:

Figure BDA00035909825000000812
Figure BDA00035909825000000812

Figure BDA00035909825000000813
Figure BDA00035909825000000813

式(11)是根据需求弹性理论计算电负荷响应量；式(12)是电冷热负荷响应量；Equation (11) is the electric load response calculated according to the demand elasticity theory; Equation (12) is the electric cooling and heating load response;

式(11)-式(12)中，

Figure BDA0003590982500000091
表示第i类型用户t时段的原始电负荷；εtt′表示t与t′时段之间的需求互弹性系数；当t＝t′时，εtt′表示需求自弹性系数；；
Figure BDA0003590982500000092
表示t′时段的基础电需求惩罚；
Figure BDA0003590982500000093
表示t′时段的电需求惩罚变化量；
Figure BDA0003590982500000094
和
Figure BDA0003590982500000095
表示第i类型用户t时段的电、热和冷负荷响应量；In formula (11) - formula (12),
Figure BDA0003590982500000091
Represents the original electrical load of the i-th type of user in period t; ε tt' represents the demand mutual elasticity coefficient between t and t'period; when t=t', ε tt' represents the demand self-elasticity coefficient;
Figure BDA0003590982500000092
represents the basic electricity demand penalty in the t'period;
Figure BDA0003590982500000093
Represents the amount of electricity demand penalty change in the t'period;
Figure BDA0003590982500000094
and
Figure BDA0003590982500000095
Represents the response of electricity, heat and cooling loads of the i-th type of user during period t;

步骤2.3、利用式(13)计算舒适度补偿量：Step 2.3, use formula (13) to calculate the amount of comfort compensation:

Figure BDA0003590982500000096
Figure BDA0003590982500000096

式(13)用于计算系统因综合需求响应产生的舒适度补偿量；Equation (13) is used to calculate the comfort compensation amount of the system due to the comprehensive demand response;

式(13)中，Ccomf表示舒适度响应补偿量，λcomf表示单位能量的舒适度补偿量；In formula (13), C comf represents the comfort response compensation amount, and λ comf represents the comfort level compensation amount per unit energy;

步骤三、建立园区综合能源系统日前调度优化模型，考虑实际约束且以低碳低能耗高用户满意度加权调和数最小为优化目标函数，对供给和需求侧进行协调优化：Step 3: Establish a day-ahead scheduling optimization model for the integrated energy system of the park, consider the actual constraints and take the minimum weighted harmonic number of low-carbon, low-energy consumption and high user satisfaction as the optimization objective function to coordinate and optimize the supply and demand sides:

步骤3.1、利用式(14)-式(16)计算低碳低能耗高用户满意度加权调和数：Step 3.1. Use equations (14)-(16) to calculate the weighted harmonic number of low carbon, low energy consumption and high user satisfaction:

min C＝Cbuy+Ccarbon+Cpena+Ccomf (14)min C=C buy +C carbon +C pena +C comf (14)

Figure BDA0003590982500000097
Figure BDA0003590982500000097

Figure BDA0003590982500000098
Figure BDA0003590982500000098

式(14)是系统运行模型目标函数；式(15)和式(16)是向上级能源网络申领量和弃风光惩罚；Equation (14) is the objective function of the system operation model; Equation (15) and Equation (16) are the amount applied to the upper energy network and the penalty for abandoning wind and solar;

式(14)-式(16)中，Cbuy是向上级能源网络申领总能量；Cpena是弃风光惩罚；

Figure BDA00035909825000000914
和
Figure BDA00035909825000000915
分别表示t时段的向上级能源网络申领的电量、申领的天然气量和申领的热量；cpena是弃风光惩罚系数；
Figure BDA0003590982500000099
和
Figure BDA00035909825000000910
是风光预测功率；Ppv,t和Pwt,t是系统t时段实际消纳的光电功率与风电功率；In Equation (14)-Equation (16), C buy is the total energy applied to the upper energy network; C pena is the penalty for abandoning the scenery;
Figure BDA00035909825000000914
and
Figure BDA00035909825000000915
Represents the amount of electricity applied for, the amount of natural gas applied for, and the amount of heat applied to the upper energy network in period t, respectively; c pena is the penalty coefficient for abandoning wind and solar;
Figure BDA0003590982500000099
and
Figure BDA00035909825000000910
is the predicted wind power; P pv,t and P wt,t are the photovoltaic power and wind power actually consumed by the system during t period;

步骤3.2、利用式(17)-式(19)定义能源供给侧约束：Step 3.2. Use equations (17)-(19) to define energy supply side constraints:

Figure BDA00035909825000000911
Figure BDA00035909825000000911

Figure BDA00035909825000000912
Figure BDA00035909825000000912

Figure BDA00035909825000000913
Figure BDA00035909825000000913

式(17)是能源供给侧电能平衡约束；式(18)是能源供给侧热能平衡约束；式(19)是能源供给侧天然气平衡约束；Equation (17) is the energy balance constraint on the energy supply side; Equation (18) is the heat energy balance constraint on the energy supply side; Equation (19) is the natural gas balance constraint on the energy supply side;

式(17)-式(19)中，

Figure BDA00035909825000001013
和
Figure BDA0003590982500000102
表示P2G机组出力、热电联产设备CHP电出力、热电联产设备CHP热出力和燃气锅炉GB出力；η表示设备效率；In formula (17) - formula (19),
Figure BDA00035909825000001013
and
Figure BDA0003590982500000102
Represents the output of P2G units, the CHP electrical output of cogeneration equipment, the CHP thermal output of cogeneration equipment, and the GB output of gas boilers; η represents equipment efficiency;

步骤3.2、利用式(20)-式(24)定义能源需求侧约束：Step 3.2, using equations (20)-(24) to define energy demand side constraints:

Figure BDA00035909825000001014
Figure BDA00035909825000001014

Figure BDA00035909825000001015
Figure BDA00035909825000001015

Figure BDA0003590982500000103
Figure BDA0003590982500000103

Figure BDA0003590982500000104
Figure BDA0003590982500000104

Figure BDA0003590982500000105
Figure BDA0003590982500000105

式(20)和式(21)表示能源供给侧与用户侧的电能和热能平衡；式(22)表示用户侧电能平衡；式(23)和式(24)表示用户侧热能和冷能平衡；Equation (20) and Equation (21) represent the balance of electric energy and thermal energy between the energy supply side and the user side; Equation (22) represents the balance of electrical energy on the user side; Equation (23) and Equation (24) represent the balance of heat energy and cold energy on the user side;

式(20)-式(24)中，

Figure BDA00035909825000001016
和
Figure BDA00035909825000001017
表示系统t时段电热泵EHP出力、中央空调AC出力、换热器HE出力和制冷机AF出力；
Figure BDA00035909825000001018
和
Figure BDA00035909825000001019
表示响应前系统t时段电负荷需求、热负荷需求和冷负荷需求；In formula (20) - formula (24),
Figure BDA00035909825000001016
and
Figure BDA00035909825000001017
Indicates the EHP output of the electric heat pump, the AC output of the central air conditioner, the HE output of the heat exchanger and the AF output of the refrigerator in the system t period;
Figure BDA00035909825000001018
and
Figure BDA00035909825000001019
Represents the electrical load demand, heating load demand and cooling load demand of the system before the response period t;

步骤3.3、利用式(25)-式(31)定义设备运行约束：Step 3.3. Use equations (25)-(31) to define equipment operating constraints:

Figure BDA0003590982500000106
Figure BDA0003590982500000106

Figure BDA0003590982500000107
Figure BDA0003590982500000107

Figure BDA0003590982500000108
Figure BDA0003590982500000108

Figure BDA0003590982500000109
Figure BDA0003590982500000109

Figure BDA00035909825000001010
Figure BDA00035909825000001010

Figure BDA00035909825000001011
Figure BDA00035909825000001011

Figure BDA00035909825000001012
Figure BDA00035909825000001012

式(25)、式(26)和式(27)是系统向上级网络申领电、热和气的上下限约束；式(28)和式(29)是热电联产设备CHP发电容量约束和发热容量约束；式(30)和式(31)是燃气轮机GB和P2G设备的容量约束；Equation (25), Equation (26) and Equation (27) are the upper and lower limit constraints for the system to apply for electricity, heat and gas to the superior network; Equation (28) and Equation (29) are the CHP power generation capacity constraints and heat generation of the cogeneration equipment. Capacity constraints; Equations (30) and (31) are the capacity constraints of gas turbine GB and P2G equipment;

步骤3.4、利用式(32)-式(33)定义风光出力约束：Step 3.4, use Equation (32)-Equation (33) to define the wind and light output constraints:

Figure BDA0003590982500000111
Figure BDA0003590982500000111

Figure BDA0003590982500000112
Figure BDA0003590982500000112

式(32)和式(33)是光电出力限制和风电出力限制；Equations (32) and (33) are the photovoltaic output limit and the wind power output limit;

步骤3.5、利用式(34)定义电需求惩罚变化量约束：Step 3.5, use formula (34) to define the electricity demand penalty variation constraint:

Figure BDA0003590982500000113
Figure BDA0003590982500000113

式(34)是电需求惩罚变化量上下限限制；Equation (34) is the upper and lower limit of the electricity demand penalty change;

步骤四、构建计及综合需求响应和碳排放的园区综合能源系统规划模型：Step 4. Build a comprehensive energy system planning model for the park that takes into account comprehensive demand response and carbon emissions:

步骤4.1、利用式(35)-式(37)构建园区综合能源系统规划模型的目标函数：Step 4.1. Use equations (35)-(37) to construct the objective function of the comprehensive energy system planning model of the park:

Figure BDA0003590982500000114
Figure BDA0003590982500000114

Figure BDA0003590982500000115
Figure BDA0003590982500000115

Figure BDA0003590982500000116
Figure BDA0003590982500000116

式(35)是系统规划模型的目标函数；式(36)和式(37)是计算投资耗材和低碳低能耗高用户满意度加权调和数；Equation (35) is the objective function of the system planning model; Equation (36) and Equation (37) are the weighted harmonic numbers for calculating investment consumables and low carbon, low energy consumption and high user satisfaction;

式(35)-式(37)中，ΩY为规划年集合；ΩD为规划的设备集合；

Figure BDA0003590982500000117
为第d个设备的类型集合；ΩS为季度的集合；Cinvest为投资耗材；
Figure BDA0003590982500000118
为第y年的运行耗材；
Figure BDA0003590982500000119
为第d个设备第c种类型的投资耗材；xc,d为布尔变量；表示是否投资第d类设备第c种类型；
Figure BDA00035909825000001110
为第y年第s季度典型日的低碳低能耗高用户满意度加权调和数，按照式(14)计算；ns为第s季度的天数；ρ为损耗率；In equations (35)-(37), Ω Y is the set of planning years; Ω D is the set of planned equipment;
Figure BDA0003590982500000117
is the type collection of the d-th equipment; Ω S is the collection of quarters; C invest is the investment consumables;
Figure BDA0003590982500000118
Consumables for the yth year;
Figure BDA0003590982500000119
is the investment consumables of the c-th type of the d-th equipment; x c, d is a Boolean variable; it indicates whether to invest in the c-th type of the d-th equipment;
Figure BDA00035909825000001110
is the weighted harmonic number of low-carbon, low-energy consumption and high user satisfaction for typical days in the sth quarter of the yth year, calculated according to formula (14); n s is the number of days in the sth quarter; ρ is the loss rate;

步骤4.2、利用式(38)定义设备的投资类型约束：Step 4.2, use formula (38) to define the investment type constraint of equipment:

Figure BDA00035909825000001111
Figure BDA00035909825000001111

式(38)表示对于任意d类设备，投资的类型不超过1种；Equation (38) indicates that for any type d equipment, there is no more than one type of investment;

步骤4.3、利用式(39)计算第y年第s季度负荷预测值Pj,y,s：Step 4.3, use formula (39) to calculate the load forecast value P j,y,s in the s quarter of the yth year:

Pj,y,s＝(1+γ)tPj,0,s,j∈ΩB,y∈ΩY,s∈ΩS (39)P j,y,s =(1+γ) t P j,0,s ,j∈Ω B ,y∈Ω Y ,s∈Ω S (39)

式(39)中，γ为负荷年增长率；Pi,0,s为当前年第s季度节点j的负荷值；In formula (39), γ is the annual growth rate of load; P i,0,s is the load value of node j in the sth quarter of the current year;

步骤五、如图3所示，采用粒子群算法求解园区综合能源系统规划模型：Step 5. As shown in Figure 3, the particle swarm algorithm is used to solve the planning model of the comprehensive energy system of the park:

步骤5.1、输入初始参数，包括粒子群种群规模M、学习因子c1和c2、惯性权重w、粒子群繁殖代数Mc、蒙特卡罗模拟次数Ms、置信区间β；Step 5.1. Input initial parameters, including particle swarm population size M, learning factors c 1 and c 2 , inertia weight w, particle swarm reproduction algebra M c , Monte Carlo simulation times M s , confidence interval β;

步骤5.2、随机生成M个初始粒子，粒子集合M＝{m1,m2,…,mk,…,mM}，其中mk表示第k个粒子，即从ΩD中选择不同的设备构成的第k个规划方案，且mk＝{mk1,mk2,…,mkd,…mkD}，其中mkd表示第k个粒子选的d类设备的容量；Step 5.2. Randomly generate M initial particles, the particle set M={m 1 , m 2 ,...,m k ,..., m M }, where m k represents the kth particle, that is, select different devices from Ω D The kth planning scheme formed, and m k ={m k1 ,m k2 ,...,m kd ,... m kD }, where m kd represents the capacity of the type d equipment selected by the kth particle;

步骤5.3、针对每个粒子，根据其规划设备容量方案mk，计算投资耗材；Step 5.3, for each particle, calculate the investment consumables according to its planned equipment capacity scheme m k ;

步骤5.4、根据式(39)，更新第y年的负荷需求，并计算第y年第s季度典型日的运行耗材；计算规划年内总的运行耗材，将计算的投资耗材与低碳低能耗高用户满意度加权调和数之和作为每个粒子的适应度；Step 5.4. According to formula (39), update the load demand in the yth year, and calculate the operating consumables for the typical day in the yth year and the s quarter; calculate the total operating consumables in the planning year, and compare the calculated investment consumables with low-carbon, low-energy consumption and high energy consumption. The sum of user satisfaction weighted harmonic numbers is used as the fitness of each particle;

步骤5.5、更新粒子位置和速度，从而获得新的粒子；Step 5.5, update the particle position and velocity to obtain new particles;

步骤5.6、重复步骤5.3-5.5，直至达到给定的粒子群繁殖代数Mc；Step 5.6, repeat steps 5.3-5.5 until reaching the given particle swarm reproduction algebra M c ;

步骤5.7、将最好的粒子所对应的设备容量作为园区综合能源系统的最优规划方案。Step 5.7, take the equipment capacity corresponding to the best particles as the optimal planning scheme of the comprehensive energy system of the park.

Claims (1)

1. A park integrated energy system planning method considering integrated demand response and carbon emission, the park integrated energy system comprises a supply side and a demand side, the supply side comprises a combined heat and power generation device CHP, a gas turbine GB and a P2G unit, the demand side comprises an electric heat pump EHP, a central air conditioner AC, a heat exchanger HE and a refrigerator AF, and the park integrated energy system planning method is characterized by comprising the following steps:
step one, calculating the initial gratuitous carbon emission of the system through a reference value method:
step 1.1, calculating the emission of the gratuitous carbon by using a formula (1) to a formula (4):
Figure FDA0003590982490000011
Figure FDA0003590982490000012
Figure FDA0003590982490000013
Figure FDA0003590982490000014
in the formulae (1) to (4),
Figure FDA0003590982490000015
representing the uncompensated initial carbon emission of applying electric quantity to a superior energy network;
Figure FDA0003590982490000016
represents the gratuitous initial carbon emission of the cogeneration plant CHP;
Figure FDA0003590982490000017
the method comprises the steps of representing the uncompensated initial carbon emission of a gas turbine unit GB;
Figure FDA0003590982490000018
represents the amount of the uncompensated initial carbon emission per unit of electricity;
Figure FDA0003590982490000019
the method comprises the steps of representing the application of electric quantity to a superior energy network within a t period;
Figure FDA00035909824900000110
represents the amount of the initial carbon emissions per unit of heat;
Figure FDA00035909824900000111
a conversion coefficient representing the conversion from the generated energy e of the cogeneration equipment to the calorific value h;
Figure FDA00035909824900000112
indicating that the heat and power cogeneration plant CHP consumes energy of natural gas for heat production in the period t;
Figure FDA00035909824900000113
indicating that the natural gas energy is consumed when the combined heat and power generation device CHP is used for generating electricity in the period t;
Figure FDA00035909824900000114
representing the natural gas energy consumed by the gas turbine unit equipment GB within the time period t; t represents a set of scheduling times;
step 1.2, calculating the actual carbon emission by using the formula (5) to the formula (6):
M＝Melcbuy+MCHP+MGB-MP2G (5)
Figure FDA00035909824900000115
in the formulae (5) to (6), MelcbuyRepresents the application of electric quantity and carbon emission to a superior energy network, MCHPDenotes the carbon emission, M, of the CHP of the cogeneration plantP2GRepresents the carbon consumption of the P2G unit; beta is agRepresents the carbon capture coefficient of the P2G unit;
Figure FDA00035909824900000116
represents the amount of power consumed by the P2G device during the t period;
step 1.3, calculating the step carbon emission C by using the formula (7)carbon：
Figure FDA0003590982490000021
In the formula (7), c represents a carbon emission reference coefficient; alpha is the growth rate of the carbon emission coefficient; d is the length of the step interval;
step two, establishing a comprehensive demand response model based on the real-time electric heating and cooling demand of the user and a demand elasticity theory:
step 2.1, calculating the cold, heat and electricity substitution coefficient by using the formula (8) to the formula (10):
Figure FDA0003590982490000022
Figure FDA0003590982490000023
Figure FDA0003590982490000024
in the formulae (8) to (10),
Figure FDA0003590982490000025
and
Figure FDA0003590982490000026
is the electricity, heat and cold direct energy demand of the ith class of users at the time period t;
Figure FDA0003590982490000027
and
Figure FDA0003590982490000028
representing the energy utilization preference of the electric cooling load and the electric heating load of the ith class of users in the t period;
Figure FDA0003590982490000029
respectively represent the thermoelectric substitution coefficient and the cold electric substitution coefficient, and is prepared from
Figure FDA00035909824900000210
And
Figure FDA00035909824900000211
carrying out normalization to obtain; deltaiRepresenting the proportion of non-rigid electrical loads of the ith type of user;
step 2.2, calculating the comprehensive demand response quantity by using the formula (11) to the formula (12):
Figure FDA00035909824900000212
Figure FDA00035909824900000213
in the formulae (11) to (12),
Figure FDA00035909824900000214
representing the original electrical load of the ith type user t period; epsilontt′Representing the required mutual elastic coefficient between the t and t' periods; when t is t', epsilontt′Representing a required self-elastic coefficient;
Figure FDA00035909824900000215
a base electrical demand penalty representing a t' period;
Figure FDA00035909824900000216
representing electrical demand penalty variations for a period t';
Figure FDA00035909824900000217
and
Figure FDA00035909824900000218
representing the electric load, heat load and cold load response of the ith type user in the t period;
step 2.3, calculating comfort compensation quantity by using the formula (13):
Figure FDA0003590982490000031
in the formula (13), CcomfIndicating a comfort response compensation quantity, λcomfA comfort compensation quantity representing a unit of energy;
establishing a day-ahead scheduling optimization model of the park comprehensive energy system, considering practical constraints and taking the weighted scheduling sum with low carbon, low energy consumption and high user satisfaction as an optimization objective function, and performing coordinated optimization on the supply side and the demand side:
step 3.1, calculating the weighted sum of the low-carbon, low-energy-consumption and high-user satisfaction by using the formula (14) to the formula (16):
min C＝Cbuy+Ccarbon+Cpena+Ccomf (14)
Figure FDA0003590982490000032
Figure FDA0003590982490000033
in the formulae (14) to (16), CbuyThe total energy is claimed to the upper-level energy network; cpenaPunishment of wind and light abandonment;
Figure FDA0003590982490000034
and
Figure FDA0003590982490000035
respectively representing electric quantity, natural gas quantity and heat quantity of application from the t time period to the superior energy network; c. CpenaThe wind and light abandonment penalty coefficient;
Figure FDA0003590982490000036
and
Figure FDA0003590982490000037
is wind, light predicted power; ppv,tAnd Pwt,tThe photoelectric power and the wind power actually consumed by the system in the t period;
step 3.2, defining the energy supply side constraint by using the formula (17) to the formula (19):
Figure FDA0003590982490000038
Figure FDA0003590982490000039
Figure FDA00035909824900000310
in the formulae (17) to (19),
Figure FDA00035909824900000311
and
Figure FDA00035909824900000312
the output of a P2G unit, the electricity output of a CHP (combined heat and power generation) device, the CHP heat output of the CHP device and the output of a gas boiler GB in the t period of the system are represented; η represents the plant efficiency;
step 3.2, defining the energy demand side constraint by using the formula (20) to the formula (24):
Figure FDA00035909824900000313
Figure FDA00035909824900000314
Figure FDA00035909824900000315
Figure FDA0003590982490000041
Figure FDA0003590982490000042
in the formulae (20) to (24),
Figure FDA0003590982490000043
and
Figure FDA0003590982490000044
the output of an electric heat pump EHP, the output of a central air conditioner AC, the output of a heat exchanger HE and the output of a refrigerating machine AF in the t period of the system are represented;
Figure FDA0003590982490000045
and
Figure FDA0003590982490000046
representing the electric load demand, the heat load demand and the cold load demand of the system in a period t before response;
step 3.3, defining the equipment operation constraint by using the formula (25) to the formula (31):
Figure FDA0003590982490000047
Figure FDA0003590982490000048
Figure FDA0003590982490000049
Figure FDA00035909824900000410
Figure FDA00035909824900000411
Figure FDA00035909824900000412
Figure FDA00035909824900000413
step 3.4, defining the wind-solar output constraint by using the formula (32) to the formula (33):
Figure FDA00035909824900000414
Figure FDA00035909824900000415
and 3.5, defining an electrical demand penalty variation constraint by using an equation (34):
Figure FDA00035909824900000416
step four, constructing a park comprehensive energy system planning model considering comprehensive demand response and carbon emission:
step 4.1, constructing an objective function of the park comprehensive energy system planning model by using the formula (35) to the formula (37):
Figure FDA00035909824900000417
Figure FDA00035909824900000418
Figure FDA00035909824900000419
in the formula (35) to the formula (37), ΩYIs a planning year set; omegaDIs a set of planned devices;
Figure FDA00035909824900000420
is a set of types for the d device; omegaSA set of quarters; cinvestIs an investment consumable;
Figure FDA0003590982490000051
is a running consumable in the y year;
Figure FDA0003590982490000052
the investment supplies of the type c are the d equipment; x is the number ofc,dThe Boolean variable indicates whether to invest the type c of the type d equipment;
Figure FDA0003590982490000053
the low-carbon low-energy-consumption high-user satisfaction weighted sum of the numbers is obtained for the typical day of the s quarter in the y year; n issNumber of days of the s quarter; rho is the loss rate;
and 4.2, defining the investment type constraint of the equipment by using an equation (38):
Figure FDA0003590982490000054
formula (38) indicates that the type of investment does not exceed 1 for any class d equipment;
step 4.3, calculating the predicted value P of the load of the s-th quarter of the y year by using the formula (39)j,y,s：
Pj,y,s＝(1+γ)tPj,0,s,j∈ΩB,y∈ΩY,s∈ΩS (39)
In the formula (39), γ represents the annual load growth rate; pi,0,sIs the load value of node j in the s quarter of the current year;
step five, solving a park comprehensive energy system planning model by adopting a particle swarm algorithm:
step 5.1, inputting initial parameters, including: particle swarm size M and learning factor c1And c2Inertia weight w, particle swarm reproduction algebra McMonte Carlo simulation times MsConfidence interval beta;
step 5.2, randomly generating M initial stagesStarting particles and forming a set of particles M ═ M1,m2,…,mk,…,mMIn which m iskKth particle, expressed from ΩDIn the k-th planning scheme, and mk＝{mk1,mk2,…,mkd,…mkDIn which m iskdRepresenting the capacity of the class d equipment selected by the kth particle;
step 5.3, calculating investment consumables according to the equipment capacity planning scheme corresponding to each particle;
step 5.4, updating the load demand in the y year according to the formula (39), and calculating the operation consumable items in the s quarter typical day in the y year; calculating total operation consumables in a planning year, and taking the sum of the calculated investment consumables and the weighted modulation number with low carbon, low energy consumption and high user satisfaction as the fitness of each particle;
step 5.5, updating the position and the speed of the particles so as to obtain new particles;
step 5.6, repeating the step 5.3 to the step 5.5 until a given particle swarm breeding algebra M is reachedcUntil the end;
and 5.7, taking the equipment capacity corresponding to the best particles as an optimal planning scheme of the park comprehensive energy system.
CN202210378037.2A 2022-04-12 2022-04-12 An Integrated Energy System Planning Approach for Parks Considering Integrated Demand Response and Carbon Emissions Pending CN114648250A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210378037.2A CN114648250A (en) 2022-04-12 2022-04-12 An Integrated Energy System Planning Approach for Parks Considering Integrated Demand Response and Carbon Emissions

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210378037.2A CN114648250A (en) 2022-04-12 2022-04-12 An Integrated Energy System Planning Approach for Parks Considering Integrated Demand Response and Carbon Emissions

Publications (1)

Publication Number Publication Date
CN114648250A true CN114648250A (en) 2022-06-21

Family

ID=81996520

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210378037.2A Pending CN114648250A (en) 2022-04-12 2022-04-12 An Integrated Energy System Planning Approach for Parks Considering Integrated Demand Response and Carbon Emissions

Country Status (1)

Country Link
CN (1) CN114648250A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114971071A (en) * 2022-06-22 2022-08-30 合肥综合性国家科学中心能源研究院(安徽省能源实验室) Sequence planning method of park integrated energy system considering wind-solar connection and electric-heat hybrid energy storage
CN115470564A (en) * 2022-10-08 2022-12-13 江苏智慧用能低碳技术研究院有限公司 Public building energy system coordination control method and control assembly thereof
CN117035244A (en) * 2023-10-10 2023-11-10 成都市智慧蓉城研究院有限公司 Space planning information acquisition method and system based on identification analysis
CN117787809A (en) * 2024-01-12 2024-03-29 中国科学院沈阳应用生态研究所 Accounting methods and systems for energy conservation and carbon reduction through biomass multi-source and multi-path energy conversion and utilization

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2015101244A4 (en) * 2015-09-08 2015-10-15 Cooper, James MR Exporting Australia's renewable energy via international trading swaps of ammonia.
CN106022503A (en) * 2016-03-17 2016-10-12 北京睿新科技有限公司 Micro-grid capacity programming method meeting coupling type electric cold and heat demand
CN113420930A (en) * 2021-06-30 2021-09-21 中国人民解放军国防科技大学 Comprehensive energy system load side optimal scheduling method and system considering multi-energy complementation
CN113487188A (en) * 2021-07-08 2021-10-08 重庆理工大学 Comprehensive energy system optimal scheduling method considering electric and gas joint price guide mechanism

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2015101244A4 (en) * 2015-09-08 2015-10-15 Cooper, James MR Exporting Australia's renewable energy via international trading swaps of ammonia.
CN106022503A (en) * 2016-03-17 2016-10-12 北京睿新科技有限公司 Micro-grid capacity programming method meeting coupling type electric cold and heat demand
CN113420930A (en) * 2021-06-30 2021-09-21 中国人民解放军国防科技大学 Comprehensive energy system load side optimal scheduling method and system considering multi-energy complementation
CN113487188A (en) * 2021-07-08 2021-10-08 重庆理工大学 Comprehensive energy system optimal scheduling method considering electric and gas joint price guide mechanism

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
晋旭东等: "考虑用户响应不确定性的园区综合能源系统分布鲁棒低碳调度", 《电力系统自动化》, vol. 47, no. 16, 25 August 2023 (2023-08-25), pages 10 - 21 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114971071A (en) * 2022-06-22 2022-08-30 合肥综合性国家科学中心能源研究院(安徽省能源实验室) Sequence planning method of park integrated energy system considering wind-solar connection and electric-heat hybrid energy storage
CN114971071B (en) * 2022-06-22 2025-02-14 合肥综合性国家科学中心能源研究院(安徽省能源实验室) Timing planning method for integrated energy system of park considering wind-solar access and electric-thermal hybrid energy storage
CN115470564A (en) * 2022-10-08 2022-12-13 江苏智慧用能低碳技术研究院有限公司 Public building energy system coordination control method and control assembly thereof
CN117035244A (en) * 2023-10-10 2023-11-10 成都市智慧蓉城研究院有限公司 Space planning information acquisition method and system based on identification analysis
CN117035244B (en) * 2023-10-10 2024-02-02 成都市智慧蓉城研究院有限公司 Space planning information acquisition method and system based on identification analysis
CN117787809A (en) * 2024-01-12 2024-03-29 中国科学院沈阳应用生态研究所 Accounting methods and systems for energy conservation and carbon reduction through biomass multi-source and multi-path energy conversion and utilization

Similar Documents

Publication Publication Date Title
CN114648250A (en) An Integrated Energy System Planning Approach for Parks Considering Integrated Demand Response and Carbon Emissions
CN109523092B (en) Multi-energy complementary combined cooling, heating and power system and cooperative scheduling method thereof
CN112464477A (en) Multi-energy coupling comprehensive energy operation simulation method considering demand response
CN111445090A (en) Double-layer planning method for off-grid type comprehensive energy system
CN108832656A (en) Multi-objective programming method for micro-energy grid based on power-to-gas and renewable energy utilization
CN110807588A (en) An optimal scheduling method for a multi-energy coupled integrated energy system
CN108154309A (en) The energy internet economy dispatching method of meter and the more load dynamic responses of cool and thermal power
CN113592365B (en) An energy optimal scheduling method and system considering carbon emissions and green power consumption
CN110826815A (en) Regional comprehensive energy system operation optimization method considering comprehensive demand response
CN117974365B (en) Multi-objective operation optimization method and system for electric heating comprehensive energy coupling system
CN109523065A (en) A kind of micro- energy net Optimization Scheduling based on improvement quanta particle swarm optimization
CN106786509A (en) Large-scale wind power and the combined heat and power dispatching method based on many scenario simulations off the net
CN111737884A (en) A multi-objective stochastic programming method for a micro-energy network with multiple clean energy sources
CN110361969B (en) Optimized operation method of cooling, heating and power comprehensive energy system
CN112671040A (en) Day-ahead optimal scheduling method of multi-energy complementary system considering maximum new energy consumption
CN107069783A (en) Heat storage electric boiler merges energy-storage system optimal control method
CN113988471A (en) Multi-objective optimization method for micro-grid operation
CN112883630B (en) Multi-microgrid system day-ahead optimization economic dispatching method for wind power consumption
Yang et al. Optimal dispatch for a combined cooling, heating and power microgrid considering building virtual energy storage
CN117172126A (en) NSGA-II capacity configuration-based optimized operation method for electric hydrogen comprehensive energy system
CN107749645A (en) A kind of method for controlling high-voltage large-capacity thermal storage heating device
CN117293915A (en) Multi-time scale optimization operation method and device for hydrogen electric coupling system
CN116502921A (en) Park comprehensive energy system optimization management system and coordination scheduling method thereof
Guo et al. Comprehensive energy system with combined heat and power photovoltaic-thermal power stations and building phase change energy storage for island regions and its coordinated dispatch strategy
CN116432824A (en) Integrated energy system optimization method and system based on multi-objective particle swarm

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