CN111950171A - Backup configuration method for gas-thermal inertia backup participation park comprehensive energy system - Google Patents

Backup configuration method for gas-thermal inertia backup participation park comprehensive energy system Download PDF

Info

Publication number
CN111950171A
CN111950171A CN202010637501.6A CN202010637501A CN111950171A CN 111950171 A CN111950171 A CN 111950171A CN 202010637501 A CN202010637501 A CN 202010637501A CN 111950171 A CN111950171 A CN 111950171A
Authority
CN
China
Prior art keywords
time
standby
heat
inertia
gas
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
CN202010637501.6A
Other languages
Chinese (zh)
Other versions
CN111950171B (en
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.)
Nanjing Dongbo Intelligent Energy Research Institute Co ltd
Original Assignee
Nanjing Dongbo Intelligent Energy Research Institute 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 Nanjing Dongbo Intelligent Energy Research Institute Co ltd filed Critical Nanjing Dongbo Intelligent Energy Research Institute Co ltd
Priority to CN202010637501.6A priority Critical patent/CN111950171B/en
Priority to LU500360A priority patent/LU500360B1/en
Priority to PCT/CN2020/107349 priority patent/WO2022000700A1/en
Publication of CN111950171A publication Critical patent/CN111950171A/en
Application granted granted Critical
Publication of CN111950171B publication Critical patent/CN111950171B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • 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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/04—Constraint-based CAD
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/14—Pipes
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • 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
    • 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—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Electricity, gas or water supply

Abstract

The invention discloses a backup configuration method of a gas-heat inertia backup park integrated energy system. The invention can fully utilize the gas-heat inertia standby in the comprehensive energy system of the park, further synthesize various standby forms to deal with the problem of system power shortage, and improve the operation economy of the system on the premise of ensuring the reliability level of the system.

Description

Backup configuration method for gas-thermal inertia backup participation park comprehensive energy system
Technical Field
The invention relates to a backup configuration method of a gas-thermal inertia backup participation park comprehensive energy system, and belongs to the technical field of comprehensive energy.
Background
As an important auxiliary service of the power system, the power system backup ensures that the system is operated safely and leaves a certain margin. In the park level comprehensive energy system, the standby is mainly used for solving the problem of system power shortage caused by uncertain factors such as new energy output prediction uncertainty, load prediction uncertainty and unit failure shutdown. Under the background of the current electric power market reform, the reliability and the economy of the operation of the comprehensive energy system are equally important, and the traditional conservative standby configuration method is necessarily abandoned, and the standby configuration of the comprehensive energy system in various standby forms is researched and integrated.
Considering the internal electrical and thermal multi-energy coupling characteristic of the comprehensive energy system, besides the traditional power generation side standby and demand side standby, the gas-thermal inertia standby can also provide power support for the system, and a new idea is provided for the standby configuration of the park-level comprehensive energy system. Wherein, gas inertia standby means that a gas transmission pipeline provides standby power for the system by releasing pipe storage and reducing pipe pressure, and thermal inertia standby means that a thermal load building provides standby power for the system by sacrificing operation comfort and reducing room temperature. The system integrates multiple standby forms, fully considers complementary configuration of gas-heat inertia standby, demand side standby and power generation side standby, can improve the operation economy of the system on the premise of ensuring the reliability level of the system, and optimizes the standby configuration of the system for coping with the power shortage of the system at lower cost.
Disclosure of Invention
In order to solve the problems, the invention provides a backup configuration method for a gas-heat inertia backup and participation park integrated energy system.
The invention adopts the following technical scheme for solving the technical problems:
a backup configuration method for a gas-thermal inertia backup participation park comprehensive energy system comprises the following steps:
(1) based on the natural gas pipeline transient model, a gas inertia standby coping system power shortage output model is established, and the method specifically comprises the following steps:
1) establishing a natural gas pipeline transient model based on a continuity equation and a momentum equation of the dynamic natural gas flow;
2) solving a pressure response model at the tail end of the natural gas pipeline based on the finite element approximation idea;
3) considering the operation constraint of a natural gas system, constructing a gas inertia standby coping system power shortage output model;
(2) considering thermal time lag, thermal loss and thermal inertia characteristics of a thermodynamic system, establishing a thermal inertia standby coping system power shortage output model, specifically comprising the following steps:
a) comprehensively considering thermal time lag, thermal loss and thermal inertia characteristics of the thermodynamic system, and establishing a thermodynamic system model;
b) based on time-frequency domain transformation, solving a thermal load building room temperature response model;
c) considering the operation constraint of a thermodynamic system, constructing a thermal inertia standby coping system power shortage output model;
(3) comprehensively considering gas heat inertia standby, power generation side standby and demand side standby, constructing a park comprehensive energy standby model by taking the minimized purchase total cost of the park comprehensive energy system standby, and performing park comprehensive energy standby configuration.
Further, the step (1) is specifically as follows:
1) the transient model of the natural gas pipeline is as follows:
Figure BDA0002568522680000021
in the formula, rho, v and P are density, flow velocity and pressure intensity of natural gas respectively, lambda, D and theta are friction coefficient and inner diameter of a pipeline and an inclination angle between the pipeline and a horizontal plane respectively, g is gravity acceleration, and x and t are time variable and space variable respectively;
2) solving a natural gas pipeline tail end pressure response model by utilizing Laplace transform:
Figure BDA0002568522680000022
wherein A, L, T is the cross-sectional area, length and temperature, R, of the natural gas pipeline, respectivelyMIs the gas constant of natural gas, Pout(t)、
Figure BDA0002568522680000023
Respectively the pressure at the end of the natural gas pipeline and the first and second derivatives, f, of the pressure as a function of time tout(t)、
Figure BDA0002568522680000024
Respectively the tail end flow of the natural gas pipeline and the first derivative thereof which change along with the time t;
3) the gas inertia standby coping system power shortage output model is as follows:
Figure BDA0002568522680000025
the constraint conditions are as follows:
Figure BDA0002568522680000026
wherein the content of the first and second substances,
Figure BDA0002568522680000027
and
Figure BDA0002568522680000028
are respectively PoutUpper and lower limits of (t), t1And t2Starting and stopping times, G, of gas inertia reserve supply, respectivelyMIs the calorific value of natural gas, f1、f2Respectively at t for the pipe1Time t2End of time flow.
Further, the step (2) is specifically as follows:
1) the thermodynamic system model is as follows:
Figure BDA0002568522680000031
in the formula, τn(t) is the thermal time lag of the transmission pipeline n corresponding to the thermal load building m along with the change of the time t, lnFor the length of the transport pipe n, vn(t) is the hot water flow rate of the transport pipe n as a function of time t,
Figure BDA0002568522680000032
for the heat loss power of the transport pipe n during heat transport,
Figure BDA0002568522680000033
for the pipe heat loss rate, T, of the transport pipe nm(t) and
Figure BDA0002568522680000034
the indoor temperature of the heat-loaded building m as a function of time t and its first derivative, Hm(t) the heating power of the heat network to the heat-loaded building m as a function of time t, Lm(t) the heat loss power of the heat-loaded building m as a function of time t, CASpecific heat capacity of indoor air for heat-loaded buildings, MAFor the indoor air quality of a heat-loaded building,
Figure BDA0002568522680000035
heat dissipation coefficient for heat-loaded buildings, Tout(t) represents the outdoor temperature of the heat load building as a function of time t, and m is 1, 2.
2) Solving a thermal load building room temperature response model by using Laplace transform:
Figure BDA0002568522680000036
3) the method comprises the following steps of (1) constructing a power model of the thermal inertia standby for coping with the system power shortage:
Figure BDA0002568522680000037
in the formula, Tm,cNormal comfortable temperature, Tm,lThe minimum comfortable temperature is set as the temperature of the air,
Figure BDA0002568522680000038
thermal inertia backup for thermal load building m, varying over time t, t3、t4Respectively the start and end times of the thermal inertia backup supply.
Further, the park comprehensive energy standby model is as follows:
Figure BDA0002568522680000039
wherein the gas-heat inertia standby mathematical model is
Figure BDA00025685226800000310
at、btRespectively the gas inertia and the thermal inertia reserve capacity price which change along with the time t,
Figure BDA00025685226800000311
reserve capacity put for gas inertia, thermal inertia, which varies with time t; the standby mathematical model at the power generation side is
Figure BDA00025685226800000312
ct、dtReserve capacity price and electricity price, R, over time t after clearing for the reserve marketS、
Figure BDA0002568522680000041
The reserve capacity of the power generation side and the actually input reserve capacity are changed along with the time t; the demand side standby mathematical model is
Figure BDA0002568522680000042
Figure BDA0002568522680000043
Spare capacity for demand-side users i, respectively, as a function of time tPrice, electricity compensation price, RD,i、
Figure BDA0002568522680000044
The demand response volume of the demand side user i and the actually input spare volume are changed along with the time t, and k is the number of the demand side users.
Compared with the prior art, the invention adopting the technical scheme has the following technical effects:
on the basis of fully considering the inertia characteristics of the comprehensive energy gas-heat system, the invention realizes the complementary configuration of gas-heat inertia standby, demand side standby and power generation side standby by utilizing the standby configuration of the gas-heat inertia standby participation park comprehensive energy system and integrating various standby forms so as to deal with the problem of system power shortage and improve the running economy of the system on the premise of ensuring the reliability level of the system.
Drawings
FIG. 1 is a general flow diagram of the process of the present invention;
FIG. 2 is a schematic view of a gas inertia system;
FIG. 3 is a schematic diagram of gas inertia backup versus system power deficit, where (a) is terminal flow, (b) is terminal gas pressure, and (c) is gas backup power;
FIG. 4 is a schematic view of a thermal inertia system;
fig. 5 is a schematic diagram of the thermal inertia backup responding to the system power shortage, wherein (a) is heat source input power, (b) is heat supply network supply power, (c) is indoor actual temperature, and (d) is thermal backup power output.
Detailed Description
The technical solutions provided by the present invention will be described in detail with reference to specific examples, and it should be understood that the following specific embodiments are only illustrative and not intended to limit the scope of the present invention.
A backup configuration method of a gas-thermal inertia backup participation park integrated energy system is shown in figure 1 and comprises the following steps:
(1) based on the natural gas pipeline transient model, a gas inertia standby coping system power shortage model is established
The natural gas tube reserve has negative feedback regulation characteristics: when the power demand of the system is increased, the flow at the tail end of the pipeline can be increased, partial pipeline storage is released to the gas turbine set, the problem of power shortage of the system is relieved, the pressure intensity of the pipeline is reduced, and the pipeline storage is reduced; when the power demand of the system is recovered to be normal, the flow at the tail end of the pipeline is recovered, the transmission pipeline stores part of the natural gas supplied by the gas source, the pressure of the pipeline rises again, and the storage of the pipeline is recovered to be normal, as shown in fig. 2. The natural gas pipe can be regarded as a dynamic standby of the system power shortage in consideration of the inertia regulation characteristic of the natural gas pipe.
1) And establishing a natural gas pipeline transient model based on a continuity equation and a momentum equation of the dynamic natural gas flow.
The natural gas pipeline transmission transient process can be characterized as follows:
Figure BDA0002568522680000051
in the formula, rho, v and P are density, flow velocity and pressure intensity of natural gas respectively, lambda, D and theta are friction coefficient and inner diameter of a pipeline and an inclination angle between the pipeline and a horizontal plane respectively, g is gravity acceleration, and x and t are time variable and space variable respectively.
2) And solving a pressure response model at the tail end of the natural gas pipeline based on the finite element approximation idea.
Based on the idea of finite element approximation, the following formula is used for simplification:
Figure BDA0002568522680000052
Figure BDA0002568522680000053
wherein f isout、finRespectively the outlet flow and the inlet flow (kg/s), Pout、PinRespectively, the outlet pressure and the inlet pressure (Pa) of the pipeline.
And (4) finishing to obtain a second order equation:
Figure BDA0002568522680000054
wherein A, L, T is the cross-sectional area, length and temperature, R, of the natural gas pipeline, respectivelyMIs the gas constant of natural gas, Pout(t)、
Figure BDA0002568522680000055
Respectively the pressure at the end of the natural gas pipeline and the first and second derivatives, f, of the pressure as a function of time tout(t)、
Figure BDA0002568522680000056
Respectively the natural gas pipeline end flow rate and the first derivative thereof as a function of time t.
Let t1The power shortage of the system is increased instantaneously at the moment, and the flow at the tail end of the pipeline is increased from f instantaneously1Rises to f2. The pressure response process at the tail end of the natural gas pipeline is a linear superposition result of step response and impulse response, the linear superposition result can be obtained by utilizing Laplace transform, when the flow at the tail end of the natural gas pipeline is instantaneously increased, the pressure at the tail end is reduced according to a negative exponential curve, and the pipeline is released along with the reduction.
3) And (4) considering the operation constraint of the natural gas system, and constructing a gas inertia standby coping system power shortage output model.
In actual operation, the gas pressure at the tail end of the natural gas pipeline must not exceed the upper and lower operating limits:
Figure BDA0002568522680000057
once the terminal gas pressure drops to t2When the time is lower than the operation lower limit, the standby power support cannot be provided for the system, the flow at the tail end of the pipeline is recovered, and the pressure at the tail end rises to a normal value according to a negative index curve.
Gas inertia standby maximum supply time
Figure BDA0002568522680000058
Can be defined as:
Figure BDA0002568522680000061
available air standby power model RG(t) the following:
Figure BDA0002568522680000062
wherein the content of the first and second substances,
Figure BDA0002568522680000063
and
Figure BDA0002568522680000064
are respectively PoutUpper and lower limits of (t), t1And t2Starting and stopping times, G, of gas inertia reserve supply, respectivelyMIs the calorific value of natural gas, f1、f2Respectively at t for the pipe1Time t2End of time flow.
In the embodiment of the invention, if the parameters are respectively v-5 m/s, RM＝519.1J/(KG·K)， T＝25℃＝273.15K，λ＝0.05，D＝0.5m，A＝0.19635m2，Pin(t)＝0.3MPa，
Figure BDA0002568522680000065
GM＝29.044MJ/kg
Let t1When the load is 1h, the system load demand is increased instantaneously, and the flow at the tail end of the pipeline is instantaneously changed from a normal value f1Increase to f2Wherein f is1＝1.2kg/s，f2＝1.3kg/s。
The data are all substituted into a second-order equation, and the solution is obtained by utilizing Laplace transform:
Pout(t)＝173000-[25385-8e-0.04827(t-1*3600)+25393e-0.0008161(t-1*3600)]
Figure BDA0002568522680000066
based on the deduction, in order to deal with the power shortage of the system, the flow at the tail end of the natural gas pipeline is increased instantaneously, the pressure at the tail end is reduced according to a negative index curve, and the stored gas is released along with the pressure; once the terminal gas pressure drops to t2When the time is lower than the operation lower limit, the standby power support cannot be continuously provided for the system, the flow at the tail end of the natural gas pipeline is recovered to be normal, and the pressure at the tail end rises to a normal value according to a negative index curve. The end flow f can be obtainedout(t), terminal pressure pout(t) gas reserve power RGThe (t) schematic diagrams are shown in fig. 3 (a) to (c).
(2) Considering thermal time lag, thermal loss and thermal inertia of a thermodynamic system, and establishing a thermal inertia standby coping system power shortage model
The heat transfer process in a known heat net is shown in fig. 4. The heat source transfers energy to the first-stage heat exchanger in the form of heating steam, and the temperature of the heat supply network is changed at the heat source without changing the network flow and the flow rate of the system. Under normal operating conditions, the building is always maintained at an optimal comfortable temperature; when the system has power shortage, the heat supply to the building is actively reduced, and the temperature of the building is continuously reduced to the lowest comfortable temperature to provide power support; when the hot standby time of the heat load building reaches the upper limit, the system immediately recovers to supply heat to the building normally, and the temperature of the building continuously rises to the optimal comfortable temperature. The inertia regulation characteristic of the heat load building is considered as a dynamic backup of the system power shortage.
a) And comprehensively considering the thermal time lag, the thermal loss and the thermal inertia characteristics of the thermodynamic system to establish a thermodynamic system model.
Figure BDA0002568522680000071
Figure BDA0002568522680000072
Figure BDA0002568522680000073
Figure BDA0002568522680000074
In the formula, the heat load building m corresponds to the transmission pipeline n, and the heat time lag tau of the transmission pipeline n along with the change of the time tn(t) length l of transport pipe nnProportional to the flow velocity v of the hot watern(t) is inversely proportional; thermal power loss in the transmission pipeline n during heat transmission
Figure BDA0002568522680000075
And heat loss rate of pipeline
Figure BDA0002568522680000076
And length of pipe lnIs in direct proportion; taking into account the thermal inertia of the building itself, Tm(t) represents the indoor temperature of the heat load building m as a function of time t,
Figure BDA0002568522680000077
is TmFirst derivative of (t), Hm(t) heating power of heat network to heat load building m, L, as a function of time tm(t) the heat loss power of the heat-loaded building m as a function of time t, CASpecific heat capacity of indoor air for heat-loaded buildings, MAFor the indoor air quality of a heat-loaded building,
Figure BDA0002568522680000078
the heat dissipation coefficient (related to building envelope heat dissipation, cold air infiltration heat dissipation and ventilation heat dissipation) of a heat load building is Tout(t) represents the outdoor temperature of the heat load building as a function of time t, and m is 1, 2.
b) And solving the heat load building room temperature response model based on time-frequency domain transformation.
Suppose an outdoor temperature T within a certain timeoutWithout change, the thermodynamic system model can be simplified to obtain:
Figure BDA0002568522680000079
t4at the moment, the energy supply power of the heat source is set to be instantaneously from the normal value H1Down to the lowest value H2Then, it can be inferred that:
Figure BDA00025685226800000710
in the formula, Tm,cNormal comfortable temperature, Tm,lIs the lowest comfort temperature.
The m indoor temperature of the heat load building is obtained as first-order step response, and the temperature is reduced according to a negative exponential curve.
c) And (4) considering the operation constraint of the thermodynamic system, and constructing a thermal inertia standby coping system power shortage output model.
Considering the practical situation, the time for the heat load building to put in the hot standby power is limited, if the maximum supply time of the hot inertia standby is
Figure BDA00025685226800000711
Then define t4At that moment, the thermal inertia standby is stopped to fill the power shortage for the system, and the thermal load building gradually rises back to the optimal comfortable temperature after a certain thermal time lag process.
Figure BDA00025685226800000712
t5＝t3+τn(t3)
In the formula, t5The time at which the power supplied by the heat supply network to the heat-loaded building changes is indicated, taking into account the thermal lag.
According to the above, a park comprehensive energy system hot standby power model R can be obtainedH(t) is represented by the formula.
Figure BDA0002568522680000081
In the formula, t3、t4Respectively the start and end times of the thermal inertia backup supply.
In the embodiment of the invention, the total length of the transmission pipeline and the flow rate of hot water are taken as ln＝2000m，vn(t)＝0.8268m/s。
The thermal time lag can be found as:
Figure BDA0002568522680000082
obtaining CAMA＝1GJ/℃，
Figure BDA0002568522680000083
Tm,c＝25℃，Tm,l＝15℃，Tout＝-5℃。
Let t3At 1h, the system is under power shortage and the heat source energy supply drops instantaneously, as shown in fig. 5 (a).
Assuming that 500 heat load buildings are all in the comprehensive energy park, the method is obtained according to the formula:
Figure BDA0002568522680000084
Figure BDA0002568522680000085
based on the above analysis and calculation, the heat supply network supplies the load thermal power Hm(T) temperature T in heat-loaded buildingm(t) Hot Standby Power
Figure BDA0002568522680000086
The schematic diagrams are shown in fig. 5 (b) to (d).
(3) And the backup configuration of the comprehensive energy of the park is carried out by cooperatively considering gas-heat inertia backup, power generation side backup and demand side backup.
The gas-heat inertia standby mathematical model is as follows:
Figure BDA0002568522680000087
in the formula, at、btRespectively the gas inertia and heat inertia reserve capacity prices which change along with the time t,
Figure BDA0002568522680000088
the spare capacity is the gas inertia and the heat inertia which change along with the time t.
The standby mathematical model at the power generation side is as follows:
Figure BDA0002568522680000089
in the formula, ct、dtReserve capacity price and electricity price, R, over time t after clearing for the reserve marketS、
Figure BDA0002568522680000091
The backup capacity of the power generation side and the backup capacity actually put in are respectively changed with time t.
The demand side standby mathematical model is:
Figure BDA0002568522680000092
in the formula (I), the compound is shown in the specification,
Figure BDA0002568522680000093
respectively, reserve capacity price, electricity compensation price, R of demand side user iD,i、
Figure BDA0002568522680000094
Respectively, demand as a function of time tThe side user i needs the transaction capacity and the actually input spare capacity, and k is the number of the required side users.
Comprehensively considering gas heat inertia standby, power generation side standby and demand side standby, constructing a park comprehensive energy standby model, and taking the lowest total cost of standby purchase in a researched time period as a target function on the premise of ensuring the reliability level of a system:
Figure BDA0002568522680000095
the technical means disclosed in the invention scheme are not limited to the technical means disclosed in the above embodiments, but also include the technical scheme formed by any combination of the above technical features. It should be noted that those skilled in the art can make various improvements and modifications without departing from the principle of the present invention, and such improvements and modifications are also considered to be within the scope of the present invention.

Claims (4)

1. A backup configuration method for a gas-thermal inertia backup participation park comprehensive energy system is characterized by comprising the following steps:
(1) based on the natural gas pipeline transient model, a gas inertia standby coping system power shortage output model is established, and the method specifically comprises the following steps:
1) establishing a natural gas pipeline transient model based on a continuity equation and a momentum equation of the dynamic natural gas flow;
2) solving a pressure response model at the tail end of the natural gas pipeline based on the finite element approximation idea;
3) considering the operation constraint of a natural gas system, constructing a gas inertia standby coping system power shortage output model;
(2) considering thermal time lag, thermal loss and thermal inertia characteristics of a thermodynamic system, establishing a thermal inertia standby coping system power shortage output model, specifically comprising the following steps:
a) comprehensively considering thermal time lag, thermal loss and thermal inertia characteristics of the thermodynamic system, and establishing a thermodynamic system model;
b) based on time-frequency domain transformation, solving a thermal load building room temperature response model;
c) considering the operation constraint of a thermodynamic system, constructing a thermal inertia standby coping system power shortage output model;
(3) comprehensively considering gas heat inertia standby, power generation side standby and demand side standby, constructing a park comprehensive energy standby model by taking the aim of minimizing the total purchase cost of the park comprehensive energy system standby, and performing park comprehensive energy standby configuration.
2. The spare configuration method for the gas-thermal inertia spare participation park integrated energy system according to claim 1, wherein the step (1) is specifically as follows:
1) the transient model of the natural gas pipeline is as follows:
Figure RE-FDA0002686440050000011
in the formula, rho, v and P are density, flow velocity and pressure intensity of natural gas respectively, lambda, D and theta are friction coefficient and inner diameter of a pipeline and an inclination angle between the pipeline and a horizontal plane respectively, g is gravity acceleration, and x and t are time variable and space variable respectively;
2) solving a natural gas pipeline tail end pressure response model by utilizing Laplace transform:
Figure RE-FDA0002686440050000012
wherein A, L, T is the cross-sectional area, length and temperature, R, of the natural gas pipeline, respectivelyMIs the gas constant of natural gas, Pout(t)、
Figure RE-FDA0002686440050000013
Respectively the pressure at the end of the natural gas pipeline and the first and second derivatives, f, of the pressure as a function of time tout(t)、
Figure RE-FDA0002686440050000014
Respectively the tail end flow of the natural gas pipeline and the first derivative thereof which change along with the time t;
3) the gas inertia standby coping system power shortage output model is as follows:
Figure RE-FDA0002686440050000021
the constraint conditions are as follows:
Figure RE-FDA0002686440050000022
wherein the content of the first and second substances,
Figure RE-FDA0002686440050000023
and
Figure RE-FDA0002686440050000024
are respectively PoutUpper and lower limits of (t), t1And t2Starting and stopping times, G, of gas inertia reserve supply, respectivelyMIs the calorific value of natural gas, f1、f2Respectively at t for the pipe1Time t2End of time flow.
3. The spare configuration method for the gas-thermal inertia spare participation park integrated energy system according to claim 1, wherein the step (2) is specifically as follows:
1) the thermodynamic system model is as follows:
Figure RE-FDA0002686440050000025
in the formula, τn(t) is the thermal time lag of the transmission pipeline n corresponding to the thermal load building m along with the change of the time t, lnFor the length of the transport pipe n, vn(t) is the hot water flow rate of the transport pipe n as a function of time t,
Figure RE-FDA0002686440050000026
for the heat loss power of the transport pipe n during heat transport,
Figure RE-FDA0002686440050000027
for the pipe heat loss rate, T, of the transport pipe nm(t) and
Figure RE-FDA0002686440050000028
the indoor temperature of the heat-loaded building m as a function of time t and its first derivative, Hm(t) the heating power of the heat network to the heat-loaded building m as a function of time t, Lm(t) the heat loss power of the heat-loaded building m as a function of time t, CASpecific heat capacity of indoor air for heat-loaded buildings, MAFor the indoor air quality of a heat-loaded building,
Figure RE-FDA0002686440050000029
heat dissipation coefficient for heat-loaded buildings, Tout(t) represents the outdoor temperature of the heat load building as a function of time t, and m is 1, 2.
2) Solving a thermal load building room temperature response model by using Laplace transform:
Figure RE-FDA00026864400500000210
3) the method comprises the following steps of (1) constructing a power model of the thermal inertia standby for coping with the system power shortage:
Figure RE-FDA00026864400500000211
in the formula, Tm,cNormal comfortable temperature, Tm,lThe minimum comfortable temperature is set as the temperature of the air,
Figure RE-FDA00026864400500000212
is changed with time tThermal inertia standby of chemical heat load building m, t3、t4Respectively the start and end times of the thermal inertia backup supply.
4. The backup configuration method for the gas-thermal inertia backup participation park integrated energy system according to claim 1, wherein the park integrated energy backup model is:
Figure RE-FDA0002686440050000031
wherein the gas-heat inertia standby mathematical model is
Figure RE-FDA0002686440050000032
at、btRespectively the gas inertia and heat inertia reserve capacity prices which change along with the time t,
Figure RE-FDA0002686440050000033
reserve capacity put for gas inertia, thermal inertia, which varies with time t; the standby mathematical model at the power generation side is
Figure RE-FDA0002686440050000034
ct、dtReserve capacity price and electricity price, R, over time t after clearing for the reserve marketS、
Figure RE-FDA0002686440050000035
The reserve capacity of the power generation side and the actually input reserve capacity are changed along with the time t; the demand side standby mathematical model is
Figure RE-FDA0002686440050000036
Figure RE-FDA0002686440050000037
The reserve capacity price of the demand-side user i, respectively, as a function of time t,Electricity compensation price, RD,i、
Figure RE-FDA0002686440050000038
The demand-side users i demand response transaction capacity and the actually-input spare capacity which change along with the time t, and k is the number of the demand-side users.
CN202010637501.6A 2020-07-03 2020-07-03 Backup configuration method for gas-thermal inertia backup participation park comprehensive energy system Active CN111950171B (en)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN202010637501.6A CN111950171B (en) 2020-07-03 2020-07-03 Backup configuration method for gas-thermal inertia backup participation park comprehensive energy system
LU500360A LU500360B1 (en) 2020-07-03 2020-08-06 Method for configuring park integrated energy system reserve with gas and thermal inertia reserves
PCT/CN2020/107349 WO2022000700A1 (en) 2020-07-03 2020-08-06 Backup configuration method for park comprehensive energy system in which gas and thermal inertia backups participate

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010637501.6A CN111950171B (en) 2020-07-03 2020-07-03 Backup configuration method for gas-thermal inertia backup participation park comprehensive energy system

Publications (2)

Publication Number Publication Date
CN111950171A true CN111950171A (en) 2020-11-17
CN111950171B CN111950171B (en) 2022-03-15

Family

ID=73339936

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010637501.6A Active CN111950171B (en) 2020-07-03 2020-07-03 Backup configuration method for gas-thermal inertia backup participation park comprehensive energy system

Country Status (3)

Country Link
CN (1) CN111950171B (en)
LU (1) LU500360B1 (en)
WO (1) WO2022000700A1 (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113256009A (en) * 2021-05-31 2021-08-13 东南大学 Two-stage multi-standby configuration method, system and device based on gas-heat inertia
CN113283107A (en) * 2021-06-10 2021-08-20 东南大学 Evaluation method and model for inertial characteristics of gas-thermal system in comprehensive energy system
CN113344357A (en) * 2021-05-31 2021-09-03 东南大学 Design method of comprehensive energy system based on frequency domain dynamic indexes
CN114881416A (en) * 2022-04-08 2022-08-09 东南大学溧阳研究院 Comprehensive energy inertial support method considering flexibility of cogeneration unit

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114048928B (en) * 2022-01-12 2022-04-19 汉谷云智(武汉)科技有限公司 Building short-term load prediction method with high migratability
CN114362124B (en) * 2022-01-14 2023-10-27 国网(苏州)城市能源研究院有限责任公司 Electric heating cooperative control method and light storage straight and flexible system
CN115081247B (en) * 2022-07-19 2023-04-07 河北农业大学 Thermodynamic system reliability evaluation method based on Bayesian network
CN116362478B (en) * 2023-02-15 2024-03-22 浙江大学 Electric-gas coupling system risk scheduling method considering comprehensive energy hub flexibility

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120232701A1 (en) * 2011-03-07 2012-09-13 Raphael Carty Systems and methods for optimizing energy and resource management for building systems
CN109376912A (en) * 2018-09-29 2019-02-22 东南大学 Cooling heating and power generation system running optimizatin method based on civil building thermal inertia
CN109474025A (en) * 2018-10-08 2019-03-15 国网能源研究院有限公司 Optimized dispatching model of park level comprehensive energy system
CN109472050A (en) * 2018-09-30 2019-03-15 东南大学 Co-generation unit incorporation time scale dispatching method based on thermal inertia
CN109978625A (en) * 2019-03-28 2019-07-05 河海大学 It is a kind of meter and electric heating gas network integrated energy system multiple target running optimizatin method
CN110503271A (en) * 2019-08-30 2019-11-26 南京工业大学 A kind of polymorphic type energy storage configuration method of integrated energy system

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9116897B2 (en) * 2009-07-20 2015-08-25 Schneider Electric It Corporation Techniques for power analysis
CN108808663B (en) * 2018-06-12 2020-04-28 浙江大学 Industrial user heat demand response method based on multi-energy complementation
CN110544017A (en) * 2019-08-12 2019-12-06 上海交通大学 Energy system reliability assessment method considering thermal inertia and energy network constraint

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120232701A1 (en) * 2011-03-07 2012-09-13 Raphael Carty Systems and methods for optimizing energy and resource management for building systems
CN109376912A (en) * 2018-09-29 2019-02-22 东南大学 Cooling heating and power generation system running optimizatin method based on civil building thermal inertia
CN109472050A (en) * 2018-09-30 2019-03-15 东南大学 Co-generation unit incorporation time scale dispatching method based on thermal inertia
CN109474025A (en) * 2018-10-08 2019-03-15 国网能源研究院有限公司 Optimized dispatching model of park level comprehensive energy system
CN109978625A (en) * 2019-03-28 2019-07-05 河海大学 It is a kind of meter and electric heating gas network integrated energy system multiple target running optimizatin method
CN110503271A (en) * 2019-08-30 2019-11-26 南京工业大学 A kind of polymorphic type energy storage configuration method of integrated energy system

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
仪忠凯等: "计及热网储热和供热区域热惯性的热电联合调度策略", 《电网技术》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113256009A (en) * 2021-05-31 2021-08-13 东南大学 Two-stage multi-standby configuration method, system and device based on gas-heat inertia
CN113344357A (en) * 2021-05-31 2021-09-03 东南大学 Design method of comprehensive energy system based on frequency domain dynamic indexes
CN113344357B (en) * 2021-05-31 2023-10-03 东南大学 Design method of comprehensive energy system based on frequency domain dynamic index
CN113256009B (en) * 2021-05-31 2023-11-28 东南大学 Two-stage multi-standby configuration method, system and device based on aerothermal inertia
CN113283107A (en) * 2021-06-10 2021-08-20 东南大学 Evaluation method and model for inertial characteristics of gas-thermal system in comprehensive energy system
CN114881416A (en) * 2022-04-08 2022-08-09 东南大学溧阳研究院 Comprehensive energy inertial support method considering flexibility of cogeneration unit

Also Published As

Publication number Publication date
LU500360B1 (en) 2022-01-10
CN111950171B (en) 2022-03-15
WO2022000700A1 (en) 2022-01-06

Similar Documents

Publication Publication Date Title
CN111950171B (en) Backup configuration method for gas-thermal inertia backup participation park comprehensive energy system
WO2018130231A1 (en) Integrated energy system optimization method based on heating grid and home thermal inertia
CN105656064B (en) Realize that thermoelectricity unit participates in the method that peak-load regulating is dispatched using heat supply time lag
CN109764573B (en) Control method of building type combined cooling heating and power system considering full utilization of waste heat
CN106593771B (en) Electric heating decoupled system and its control method in a kind of co-generation unit
CN109359839B (en) Node heat price calculation method for regional heat supply network in comprehensive energy system
CN110188492B (en) Combined cooling heating and power micro-grid optimized scheduling method considering heat supply network characteristics
CN112529244A (en) Comprehensive energy system collaborative optimization operation method considering electric load demand response
EP1978311A2 (en) System solar heating autonomous and independent of another energy source
CN116308881A (en) Multi-time scale scheduling method for comprehensive energy system utilizing heat supply pipe network for heat storage
CN113283107A (en) Evaluation method and model for inertial characteristics of gas-thermal system in comprehensive energy system
CN106918069B (en) Control method of phase-change heat storage heating system
US20210190328A1 (en) Combined heating and cooling system
CN109253492B (en) Adjusting method for solving hydraulic imbalance of heat supply system
CN106440902A (en) Heat storage adjusting system based on combined cooling heating and power and adjustment method thereof
CN114218765A (en) Comprehensive energy inertial support method for coping with power shortage of power grid
CN110781598A (en) Method for calculating running state of regional comprehensive energy system
CN110716429A (en) Control method and device of combined cooling heating and power system, computer and storage medium
JP2007042377A (en) Fuel cell power generation system
CN205481281U (en) Intelligence millet electricity heat accumulation heating system
CN112115587A (en) Generator set load instruction optimization system based on storage utilization
CN219389683U (en) Regulation and control device of cogeneration unit
CN208124422U (en) Card slot type electric heating accumulation of heat heating system
FR2902180A3 (en) Air utilizing device for e.g. individual dwelling, has ventilation extractor heating cool loop of heat pump via battery, where thermal balance of cool loop is ensured completely by external air via cooler or compressor package
CN112100828B (en) Electric heating system control method considering load quasi-dynamic characteristic of heating power network

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