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 PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 22
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 94
- 239000003345 natural gas Substances 0.000 claims description 47
- 239000007789 gas Substances 0.000 claims description 37
- 230000010485 coping Effects 0.000 claims description 17
- 230000004044 response Effects 0.000 claims description 16
- 238000010248 power generation Methods 0.000 claims description 13
- 230000008859 change Effects 0.000 claims description 11
- 230000005540 biological transmission Effects 0.000 claims description 9
- 238000013178 mathematical model Methods 0.000 claims description 9
- 230000001052 transient effect Effects 0.000 claims description 9
- 230000005611 electricity Effects 0.000 claims description 6
- 230000017525 heat dissipation Effects 0.000 claims description 6
- 238000010438 heat treatment Methods 0.000 claims description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000033228 biological regulation Effects 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- PKOMXLRKGNITKG-UHFFFAOYSA-L calcium;hydroxy(methyl)arsinate Chemical compound [Ca+2].C[As](O)([O-])=O.C[As](O)([O-])=O PKOMXLRKGNITKG-UHFFFAOYSA-L 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000006735 deficit Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 230000009125 negative feedback regulation Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000009423 ventilation Methods 0.000 description 1
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—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)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Primary Health Care (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Fluid Mechanics (AREA)
- Computing Systems (AREA)
- General Health & Medical Sciences (AREA)
- Educational Administration (AREA)
- Algebra (AREA)
- Water Supply & Treatment (AREA)
- Public Health (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
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
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:
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:
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)、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)、 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:
wherein,andare 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:
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,for the heat loss power of the transport pipe n during heat transport,for the pipe heat loss rate, T, of the transport pipe nm(t) andthe 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,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:
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:
in the formula, Tm,cNormal comfortable temperature, Tm,lThe minimum comfortable temperature is set as the temperature of the air,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:
wherein the gas-heat inertia standby mathematical model isat、btRespectively the gas inertia and the thermal inertia reserve capacity price which change along with the time t,reserve capacity put for gas inertia, thermal inertia, which varies with time t; the standby mathematical model at the power generation side isct、dtReserve capacity price and electricity price, R, over time t after clearing for the reserve marketS、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 Spare capacity for demand-side users i, respectively, as a function of time tPrice, electricity compensation price, RD,i、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:
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:
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:
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)、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)、 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:
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.
available air standby power model RG(t) the following:
wherein,andare 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, 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)]
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.
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 transmissionAnd heat loss rate of pipelineAnd 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,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,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:
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:
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 isThen 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.
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.
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:
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:
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 PowerThe 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:
in the formula, at、btRespectively the gas inertia and heat inertia reserve capacity prices which change along with the time t,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:
in the formula, ct、dtReserve capacity price and electricity price, R, over time t after clearing for the reserve marketS、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:
in the formula,respectively, reserve capacity price, electricity compensation price, R of demand side user iD,i、 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:
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:
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:
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)、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)、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:
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:
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,for the heat loss power of the transport pipe n during heat transport,for the pipe heat loss rate, T, of the transport pipe nm(t) andthe 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,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:
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:
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:
wherein the gas-heat inertia standby mathematical model isat、btRespectively the gas inertia and heat inertia reserve capacity prices which change along with the time t,reserve capacity put for gas inertia, thermal inertia, which varies with time t; the standby mathematical model at the power generation side isct、dtReserve capacity price and electricity price, R, over time t after clearing for the reserve marketS、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 The reserve capacity price of the demand-side user i, respectively, as a function of time t,Electricity compensation price, RD,i、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.
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 |
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 |
LU500360A LU500360B1 (en) | 2020-07-03 | 2020-08-06 | Method for configuring park integrated energy system reserve with gas and thermal inertia reserves |
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)
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 (8)
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 |
CN114741988A (en) * | 2022-04-26 | 2022-07-12 | 东南大学 | Frequency modulation method considering gas-heat inertia of comprehensive energy system |
CN115081247B (en) * | 2022-07-19 | 2023-04-07 | 河北农业大学 | Thermodynamic system reliability evaluation method based on Bayesian network |
CN116681545B (en) * | 2023-01-30 | 2024-06-21 | 兰州理工大学 | Facility agriculture park near-zero carbon implementation method considering biomass-P2G coupling |
CN116362478B (en) * | 2023-02-15 | 2024-03-22 | 浙江大学 | Electric-gas coupling system risk scheduling method considering comprehensive energy hub flexibility |
CN117669982B (en) * | 2023-12-18 | 2024-05-17 | 山东建筑大学 | Park comprehensive energy system capacity configuration method considering user multi-comfort |
CN117993896B (en) * | 2024-04-07 | 2024-07-12 | 浙江大学 | Comprehensive energy system toughness improving method considering thermal inertia under extreme ice and snow disasters |
Citations (6)
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 | 南京工业大学 | Multi-type energy storage configuration method of comprehensive energy system |
Family Cites Families (3)
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 |
-
2020
- 2020-07-03 CN CN202010637501.6A patent/CN111950171B/en active Active
- 2020-08-06 LU LU500360A patent/LU500360B1/en active IP Right Grant
- 2020-08-06 WO PCT/CN2020/107349 patent/WO2022000700A1/en active Application Filing
Patent Citations (6)
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 | 南京工业大学 | Multi-type energy storage configuration method of comprehensive energy system |
Non-Patent Citations (1)
Title |
---|
仪忠凯等: "计及热网储热和供热区域热惯性的热电联合调度策略", 《电网技术》 * |
Cited By (6)
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 |
WO2022000700A1 (en) | 2022-01-06 |
CN111950171B (en) | 2022-03-15 |
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 | |
CN109190785B (en) | Operation optimization method for electric-thermal coupling comprehensive energy system | |
CN112529244B (en) | Comprehensive energy system collaborative optimization operation method considering electric load demand response | |
CN109447323B (en) | Comprehensive energy system two-stage capacity configuration method considering node heat price | |
CN109764573B (en) | Control method of building type combined cooling heating and power system considering full utilization of waste heat | |
CN110188492B (en) | Combined cooling heating and power micro-grid optimized scheduling method considering heat supply network characteristics | |
CN109359839B (en) | Node heat price calculation method for regional heat supply network in comprehensive energy system | |
CN114218765A (en) | Comprehensive energy inertial support method for coping with power shortage of power grid | |
CN110474335B (en) | Comprehensive energy system operation method based on evaluation criterion | |
CN113283107A (en) | Evaluation method and model for inertial characteristics of gas-thermal system in comprehensive energy system | |
CN116308881A (en) | Multi-time scale scheduling method for comprehensive energy system utilizing heat supply pipe network for heat storage | |
CN114962222B (en) | Advanced adiabatic compressed air energy storage energy hub equipment and modeling method thereof | |
CN106440902B (en) | A kind of heat accumulation regulating system and its adjusting method based on cogeneration | |
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 | |
CN210717770U (en) | Heat storage water tank system for peak regulation of heat supply load | |
CN113761727A (en) | Method for constructing optimal scheduling model of thermoelectric combined system with distributed electric heat pumps | |
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 | |
CN112100828B (en) | Electric heating system control method considering load quasi-dynamic characteristic of heating power network | |
EP3892925B1 (en) | A heating system, a method, a computer program, a computer-readable medium and a control device | |
CN112050273B (en) | Hot water pipe network virtual energy storage control method based on quasi-dynamic time delay characteristic |
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 |