LU500360B1 - Method for configuring park integrated energy system reserve with gas and thermal inertia reserves - Google Patents
Method for configuring park integrated energy system reserve with gas and thermal inertia reserves Download PDFInfo
- Publication number
- LU500360B1 LU500360B1 LU500360A LU500360A LU500360B1 LU 500360 B1 LU500360 B1 LU 500360B1 LU 500360 A LU500360 A LU 500360A LU 500360 A LU500360 A LU 500360A LU 500360 B1 LU500360 B1 LU 500360B1
- Authority
- LU
- Luxembourg
- Prior art keywords
- reserve
- thermal
- time
- refer
- inertia
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000010248 power generation Methods 0.000 claims abstract description 14
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 102
- 239000007789 gas Substances 0.000 claims description 61
- 239000003345 natural gas Substances 0.000 claims description 51
- 230000005540 biological transmission Effects 0.000 claims description 20
- 230000004044 response Effects 0.000 claims description 17
- 238000013178 mathematical model Methods 0.000 claims description 9
- 230000009466 transformation Effects 0.000 claims description 9
- 230000001052 transient effect Effects 0.000 claims description 9
- 230000017525 heat dissipation Effects 0.000 claims description 6
- 238000010438 heat treatment Methods 0.000 claims description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 5
- 230000001133 acceleration Effects 0.000 claims description 3
- 101150067055 minC gene Proteins 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 241001125929 Trisopterus luscus Species 0.000 description 2
- 230000033228 biological regulation Effects 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000004458 analytical method Methods 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
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 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
- 238000005457 optimization Methods 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
- 238000009423 ventilation Methods 0.000 description 1
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)
- Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Development Economics (AREA)
- Computer Hardware Design (AREA)
- Game Theory and Decision Science (AREA)
- Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Educational Administration (AREA)
- Mathematical Analysis (AREA)
- Mathematical Physics (AREA)
- Fluid Mechanics (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
Disclosed is a method for configuring a park integrated energy system reserve with gas and thermal inertia reserves. The method comprises the following steps: firstly building an output model for a gas inertia reserve to cope with system power vacancy, then building an output model for a thermal inertia reserve to cope with the system power vacancy, and finally configuring the park integrated energy reserve by cooperatively considering the gas and thermal inertia reserves, a power generation side reserve and a demand side reserve. The present invention can fully utilize the gas and thermal inertia reserves in the park integrated energy system, and integrate multiple forms of reserves to cope with the system power vacancy, so that the system operation economy is improved on the premise that the system reliability level is assured.
Description
[0001] The present invention relates to a method for configuring a park integrated energy system reserve with gas and thermal inertia reserves, and belongs to the technical field of integrated energy.
[0002] As an important auxiliary service of a power system, power system reserves ensure the safe operation of the system and leave a certain adequacy. In a park-level integrated energy system, reserves are mainly used to cope with the system power vacancy caused by uncertain factors such as new energy output forecast uncertainty, load forecast uncertainty and unit failure outage. In the context of current power market reform, the reliability and economy of the integrated energy system are equally important, and it is necessary to abandon a traditional conservative reserve configuration method and study integrated energy system reserve configuration that integrates multiple forms of reserves.
[0003] In view of the internal electrical and thermal multi-energy coupling characteristics of the integrated energy system, in addition to traditional power generation side reserve and demand side reserve, the gas and thermal inertial reserves can also provide power support for the system, thereby providing a new idea for the reserve configuration of the park-level integrated energy system. The gas inertia reserve refers to the provision of reserve power for the system by a gas transmission pipeline through releasing pipe storage and reducing pipe pressure, and the thermal inertia reserve refers to the provision of reserve power for the system by thermal load buildings through sacrificing the operating comfort and lowering the room temperature. Multiple forms of reserves are integrated, and complementary configuration of the gas and thermal inertia reserves, the demand side reserve and the power generation side reserve is fully considered, so that the system operation economy can be improved under the premise of ensuring the system reliability level, and the reserve configuration for the system to cope with the system power vacancy 1s optimized with the relatively low cost. 7500360
[0004] In order to solve the above problem, the present invention provides a method for configuring a park integrated energy system reserve with gas and thermal inertia reserves. According to the present invention, the gas and thermal inertia reserves in a park integrated energy system can be fully utilized for reserve optimization configuration to cope with the system power vacancy, so that the system operation economy is improved on the premise that the system reliability level is assured.
[0005] The present invention adopts the following technical scheme for solving the above technical problem:
[0006] A method for configuring a park integrated energy system reserve with gas and thermal inertia reserves, comprising the following steps:
[0007] (1) building an output model for a gas inertia reserve to cope with system power vacancy based on a natural gas pipeline transient model, specifically comprising:
[0008] 1) building the natural gas pipeline transient model based on a continuity equation and a momentum equation of dynamic natural gas flow;
[0009] 2) solving a natural gas pipeline tail end pressure response model based on the thought of finite element approximation;
[0010] 3) building the output model for the gas inertia reserve to cope with the system power vacancy by considering the operation constraints of a natural gas system;
[0011] (2) building an output model for a thermal inertia reserve to cope with the system power vacancy by considering the thermal time lag, thermal loss and thermal inertia characteristics of a thermodynamic system, specifically comprising:
[0012] a) building a thermodynamic system model by comprehensively considering the thermal time lag, thermal loss and thermal inertia characteristics of the thermodynamic system;
[0013] b) solving a thermal load building room temperature response model based on time-frequency domain transformation;
[0014] c) building the output model for the thermal inertia reserve to cope with the system power vacancy by considering the operation constraints of the thermodynamic system;
[0015] and (3) building a park integrated energy reserve model and configuring the park 0980 integrated energy reserve by comprehensively considering the gas and thermal inertia reserves, a power generation side reserve and a demand side reserve to minimize the total park integrated energy system reserve purchase cost.
[0016] Further, the step (1) specifically comprises the following steps:
[0017] 1) the natural gas pipeline transient model is as follows: | = + = =0 Op | OpvT | OP | Av” | eg sin 0=0
[0018] ot Ox ox 2D
[0019] wherein, #, Y and P refer to density, flow velocity and pressure of natural gas respectively, À, D and © refer to friction coefficient of a pipeline, inner diameter of the pipeline and inclination angle between the pipeline and the horizontal plane respectively, & refers to gravitational acceleration, and * and ! refer to time variable and space variable respectively;
[0020] 2) solving the natural gas pipeline tail end pressure response model by utilizing Laplace transformation: 2 LOS PO + P= PO] F0 Zt
[0021] M Ar
[0022] wherein, 4, L and 7 refer to the cross sectional area, the length and the temperature of the natural gas pipeline respectively, Rır refers to a natural gas gas constant, Fou UN Pau) and Pau (1) refer to the natural gas pipeline tail end pressure intensity changed along with the time ?, the first-order derivative and the second-order derivative of the natural gas pipeline tail end pressure intensity respectively, and Four (1) and Fous) refer to the natural gas pipeline tail end flow changed along with the time ! and the first-order derivative of the natural gas pipeline tail end flow respectively;
[0023] and 3) the output model for the gas inertia reserve to cope with the system power vacancy 1s as follows:
[0024]
0 O<t<t,/>t RS (f) — | G ( 1 2)
[0025] ml I (t <t<t,)
[0026] the constraint conditions are as follows: ‘ou = Ful) = Fy
[0027] wherein, Fou and Fou refer to an upper limit and a lower limit of Fu) respectively, i and © refer to the start time and the cut-off time of gas inertia reserve supply respectively, Gy is a natural gas thermal value, and Ji and ha refer to the tail end flow of the pipeline at the time 1 and the time ” respectively.
[0028] Further, the step (2) specifically comprises the following steps:
[0029] 1) the thermodynamic system model is as follows: / Ta (0) =— v, (0) loss! = ul x1, 1,0 =[H,0-L, OV CM,
[0030] L,()=e;" lz, 0-1, 0)
[0031] wherein, %» () refers to the thermal time lag, changed along with the time ‘ of a transmission pipeline n corresponding to a thermal load building m, L, refers to the length of the transmission pipeline n, va (1) refers to the hot water flow velocity, changed with the time 7, of
H the transmission pipeline n, loss, refers to the thermal loss power of the transmission pipeline n
H in a thermal transmission process, Fin refers to the pipeline thermal loss rate of the transmission y T (1) TL, . . 1 pipeline n, ” and 7m refer to the indoor temperature, changed with the time /, of the thermal load building m and the first-order derivative of the indoor temperature respectively, HC) refers to the heating power, changed with the time /, of a thermal network to the thermal load building m, Lu) refers to the heat loss power, changed with the time /, of the thermal load building m, Ca refers to the indoor air specific heat capacity of the thermal load building,
Joss LU500360 M, refers to the indoor air quality of the thermal load building, n refers to a heat dissipation coefficient of the thermal load building, Lou) represents the outdoor temperature, changed with the time ?, of the thermal load building, m=12,..z , and 7 refer to the total number of thermal load buildings in the park integrated energy system; 5 [0032] 2) solving the thermal load building room temperature response model by utilizing Laplace transformation:
[0033] CMT, (0) +e,"T, (1) =H, (04e Tu :
[0034] and 3) building the output model for the thermal inertia reserve to cope with the system power vacancy: 0 (O<t<t, t>1,) R° () = > > 0SS S RIM)=D €, (ne Tan) (t; <t<1)
[0035] m=l m=l
[0036] wherein, ”“ refers to the normal comfortable temperature, ” is the lowest R(t) . oo . . Î comfortable temperature, ” is thermal inertia reserve, changed along with the time ’, of the thermal load building m, and f and # refer to the start time and the cut-off time of thermal inertia reserve supply respectively.
[0037] Further, the park integrated energy reserve model is as follows: i T minC (6) = H[C,) + C,() + C,(0]
[0038] I
[0039] wherein, a gas and thermal inertia reserve mathematical model is _ G H GO) =aR’ +b. , 3 and b, refer to a gas inertia reserve capacity price and a thermal
RÉ RI inertia reserve capacity price changed along with the time ! respectively, and and refer to the gas inertia input reserve capacity and the thermal inertia input reserve capacity changed along with the time ! respectively; a power generation side reserve mathematical model . =cR +dR° LL oo is GO=q dR , % and d, refer to a reserve capacity price and an electric quantity price R° R° changed along with the time ’ after a reserve market is cleared, and ‘* refer to the transaction capacity and the actually input reserve capacity, changed along with the time / Fr 00980 power generation side reserve respectively; and a demand side reserve mathematical model is k OS era) 4 4 i=l , ‘ and ‘ refer to a reserve capacity price and an electric quantity
RPA compensation price, changed along with the time /, of a demand side user i respectively, and Rr refer to the demand response transaction capacity and the actually input reserve capacity, changed along with the time ’, of the demand side user i respectively, and k is the number of demand side users.
[0040] Compared with the prior art, the above technical scheme has the following technical effects:
[0041] according to the present invention, the gas and thermal inertia reserves are utilized for the reserve configuration of the park integrated energy system on the basis of fully considering the gas and thermal inertia characteristics of the integrated energy system, and multiple forms of reserves are integrated for complementary configuration of the gas and thermal inertia reserves, the demand side reserve and the power generation side reserve to cope with the system power vacancy, so that the system operation economy is improved on the premise that the system reliability level is assured.
[0042] Fig. 1 is a general flow chart of the method of the present invention;
[0043] Fig. 2 is a schematic diagram of a gas inertia system;
[0044] Fig. 3 is a schematic diagram for the gas inertia reserve to cope with the system power vacancy, wherein, (a) is the tail end flow, (b) is the tail end gas pressure, and (c) is the gas reserve power;
[0045] Fig. 4 is a schematic diagram of a thermal inertial system.
[0046] Fig. 5 is a schematic diagram for the thermal inertia reserve to cope with the system power vacancy, wherein, (a) is the input power of a heat source, (b) is the supply power of a thermal network, (c) is the actual indoor temperature, and (d) is the power output of the thermal reserve.
[0047] The technical scheme provided by the present invention is described in detail in combination with the embodiments below, and it should be understood that the following embodiments are only used to illustrate the present invention and not to limit the scope of the present invention.
[0048] A method for configuring a park integrated energy system reserve with gas and thermal inertia reserves, as shown in the Fig. 1, comprising the following steps:
[0049] (1) building an output model for a gas inertia reserve to cope with system power vacancy based on a natural gas pipeline transient model,
[0050] wherein a natural gas pipeline storage reserve has negative feedback regulation characteristics: when the system power demand is increased, the tail end flow of a pipeline can be increased, part of pipeline storage is released to a gas turbine unit, the system power vacancy is alleviated, the pressure of the pipeline is lowered, and the pipeline storage is reduced; when the system power demand returns to normal, the tail end flow of the pipeline is restored, a transmission pipeline stores part of natural gas supplied by a gas source, the pressure of the pipeline rises again, and the pipeline storage returns to the normal value, as shown in Fig. 2; in view of the inertia regulation characteristics of the natural gas pipeline storage, the gas inertia reserve can be regarded as a dynamic reserve for the system power vacancy;
[0051] 1) building the natural gas pipeline transient model based on a continuity equation and a momentum equation of dynamic natural gas flow,
[0052] wherein a natural gas pipeline transient process can be characterized as follows: = + = =0 êpy | Op” OP | Any in 0=0
[0053] ot ox ox 2D
[0054] wherein, ©, ” and P refer to density, flow velocity and pressure of natural gas respectively, 4 D and 9 refer to friction coefficient of a pipeline, inner diameter of the pipeline and inclination angle between the pipeline and the horizontal plane respectively, & refers to gravitational acceleration, and * and ’ refer to time variable and space variable respectively; 7500360
[0055] 2) solving a natural gas pipeline tail end pressure response model based on the thought of finite element approximation,
[0056] wherein simplification is performed by using the following formula based on the thought of finite element approximation: F Sort Sin
[0057] Ox L oP Pat Prin
[0058] aL
[0059] wherein, Fou and En refer to the outlet flow and the inlet flow (kg/s) of the pipeline respectively, Pau and Pr refer to the outlet pressure and the inlet pressure (Pa) of the pipeline respectively;
[0060] a second-order equation can be obtained by arrangement: 2 LOS PO + P= PO] F0 Zt
[0061] M M
[0062] wherein, 4, L and 7 refer to the cross sectional area, the length and the temperature of the natural gas pipeline respectively, Rır refers to a natural gas gas constant, Fou) ; Po) and Pau (1) refer to the natural gas pipeline tail end pressure changed along with the time 7, the first-order derivative and the second-order derivative of the natural gas pipeline tail end pressure respectively, and four) and fou () refer to the natural gas pipeline tail end flow changed along with the time * and the first-order derivative of the natural gas pipeline tail end flow respectively;
[0063] when the system power vacancy 1s set to be increased instantaneously at the time f, , the tail end flow of the pipeline is increased from h wo /o instantaneously; a natural gas pipeline tail end pressure response process is a linear superposition result of step response and impulse response, which can be solved by using Laplace transformation, and when the tail end flow of the natural gas pipeline is increased instantaneously, the tail end gas pressure is reduced according to a negative exponential curve and the pipeline storage is released accordingly;
[0064] 3) building the output model for the gas inertia reserve to cope with the system a . . LU500360 power vacancy by considering the operation constraints of a natural gas system,
[0065] wherein the tail end gas pressure of the natural gas pipeline shall not exceed the operating upper and lower limits in actual operation:
[0066] Do = Pout (1) s Pout
[0067] once the tail end gas pressure drops to be lower than the operating lower limit at the time th) reserve power support cannot be provided for the system, the tail end flow of the pipeline is restored, and the tail end gas pressure rises to the normal value according to the negative exponential curve;
[0068] the longest supply time “mx of the gas inertia reserve can be defined as follows: Ton = —1 10069] PS R°() ;
[0070] the gas reserve power model can be obtained as follows: 0 O<t<t,/>t RS (f) — | G ( 1 2)
[0071] u 1) (t <t<t,)
[0072] wherein, Fou and Tow refer to an upper limit and a lower limit of respectively, hand © refer to the start time and the stop time of gas inertia reserve supply respectively, Gyr 1s a natural gas thermal value, and Ji and À refer to the tail end flow of the pipeline at the time t 1 4 I and the time ? respectively;
[0073] in the embodiments, if Y=5m/s, Ru =519.19 /(KG-K) T =25C=273.15K, A=005 | D=05m _ A=0.19635m* P,()=03MPa Pac =0.145MPa and Gy, =29.044MJ /kg in parameters,
[0074] and when 4 =1h , the system load demand is increased, the tail end flow of the pipeline is increased from the normal value hwo fo , wherein, h=12K8/8 and f, =13kg/s.
[0075] the above data are all substituted into the second-order equation, and the following formulas can be obtained by using the Laplace transformation:
10076] P,, (f) = 173000 - [25385 _ Be 0.048276 153600) + 25393 e70-0008161 0113600) | OMW (O<t<t,/>t,) R°(0)= |
[0077] 2.9044MW (t, <t<t,)
[0078] based on the above deduction, in order to cope with the system power vacancy, the tail end flow of the natural gas pipeline is increased instantaneously, the tail end gas pressure is reduced according to the negative exponential curve, and the pipeline storage is released accordingly; once the tail end gas pressure drops to be lower than the operating lower limit at the time th) the reserve power support cannot be provided for the system, the tail end flow of the pipeline is restored, and the tail end gas pressure rises to the normal value according to the negative exponential curve; the tail end flow Fou 0) the tail end gas pressure Pow () and the gas reserve power R° 0) can be obtained and are as shown in (a) to (c) in the Fig. 3;
[0079] (2) building an output model for a thermal inertia reserve to cope with the system power vacancy by considering the thermal time lag, thermal loss and thermal inertia characteristics of a thermodynamic system,
[0080] wherein a heat transfer process in the known thermal network 1s as shown in Fig. 4; a quality adjustment method is assumed to be adopted, a heat source transfers energy to a primary heat exchanger in a form of heating steam, only the water supply temperature of the thermal network is changed at the heat source, and the network flow and flow velocity of the system are not changed; under normal operating conditions, a building is always maintained at the optimal comfortable temperature; when the system has power vacancy, heating of the building is actively reduced, and the building temperature continuously drops to the lowest comfortable temperature to provide power support, when the input thermal reserve time of the thermal load building reaches the upper limit, the system immediately restores normal heating of the building, and the building temperature continuously rises to the optimal comfortable temperature; in view of the inertia regulation characteristics of the thermal load building, the thermal inertia reserve can be regarded as a dynamic reserve for the system power vacancy;
[0081] a) building a thermodynamic system model by comprehensively considering the thermal time lag, thermal loss and thermal inertia characteristics of the thermodynamic system,
I LU500360 7,0) = v(t) loss! = ul" x1, 1,0) =|H,0-L,olc.M,
[0082] L(t) =r, ()-T,,(0]
[0083] wherein, a transmission pipeline corresponding to the thermal load building m is a pipeline n, and the thermal time lag, changed with the time /, of the transmission pipeline n is in direct proportion to the length Ly of the transmission pipeline n and is in inverse proportion to y Va (D. loss” .__ , the hot water flow velocity ”* ” ; the thermal loss power » of the transmission pipeline n
H in the thermal transmission process is in direct proportion to the thermal loss rate “ and the length ', of the pipeline; in view of thermal inertia of the building, 1,0) represents the indoor temperature, changed with the time ‘, of the thermal load building m, 1,0) is the first-order derivative of 7m 0) HC) refers to the heating power, changed along with the time 7, of the thermal network to the thermal load building m, Lu) refers to the heat loss power, changed along with the time /, of the thermal load building m, Ca refers to the indoor air specific heat capacity of the thermal load building, Ma refers to the indoor air quality of the thermal load loss building, é» refers to a heat dissipation coefficient (related to heat dissipation of a building envelope, heat dissipation of cold air infiltration and heat dissipation of ventilation) of the thermal load building, Lou (1) represents the outdoor temperature, changed along with the time /, of the thermal load building, ™ 7 1,2,..., 7 ‚and Z refers to the total number of thermal load buildings of the park integrated energy system;
[0084] b) solving a thermal load building room temperature response model based on time-frequency domain transformation,
[0085] wherein if the outdoor temperature Tow is assumed to be unchanged within a certain time, the thermodynamic system model can be obtained by simplification:
[0086] CMT, (0)+e,T,(1)=H,(0)+&,"T,,
[0087] when the energy supply power of the heat source is set to be reduced to the 00980 minimum value > from the normal value Hı instantaneously at the time [a the following formula can be deduced:
[0088] CM Th) + ERTL) = (ne ~T, DÉC) + ET
[0089] wherein, Loe refers to the normal comfortable temperature, and Tons refers to the lowest comfortable temperature;
[0090] the indoor temperature of the thermal load building m is first-order step response, and the temperature is reduced according to the negative exponential curve;
[0091] c) building the output model for the thermal inertia reserve to cope with the system power vacancy by considering the operation constraints of the thermodynamic system,
[0092] wherein in view of the limited input thermal reserve power time of the thermal load building under the actual condition, if the longest thermal inertia reserve supply time is Tas ‚the time %# is defined, the thermal inertia reserve stops to fill the system with the power vacancy, and the temperature of the thermal load building gradually rises to the optimal comfortable temperature after a certain thermal time lag process; t, = +T
[0093] f,=t +T,(t,)
[0094]
[0095] wherein, fs represents the time for changing the power supplied to the thermal load building by the thermal network in view of the thermal time lag;
[0096] according to the preceding text, the thermal reserve power model R° (1) of the park integrated energy system can be obtained and is the following formula: 0 (O<t<ty t>1,) H — zZ z #6) | SRINO=>er(T,.-T,) (t,<t<H)
[0097] mal mal
[0098] wherein, 5 and #4 refer to the start time and the cut-off time of thermal inertia reserve supply respectively;
[0099] in the embodiments, the total length and the hot water flow velocity of the se TE I, =2000m v (1) = 0.8268m/s alee transmission pipeline are ” and ” respectively;
[0100] the thermal time lag can be calculated out and is as follows: [ 2000 T,(1)=—— = —— = 2418.96s = 0.672h
[0101] v,(t) 0.8268
[0102] CM, =1GJ/C gl = 0.5075 kJ /sC Las = 25°C La =15C Tout =-5C.
[0103] when #7 1h the system has the power vacancy, and the energy supply of the heat source is reduced instantaneously, as shown in (a) of the Fig. 5;
[0104] when a total of 500 thermal load buildings in an integrated energy park is assumed, the following formulas can be obtained according to the previous formula:
0.5075 41 672*3600) T,(t)=-101-e 100 +25
[0105] OMW O<t<t,,/>t Ro | 2.5375MW | oY
[0106] (t, <t<t,)
[0107] according to the above analysis and calculation, the thermal power HC) of supply load of the thermal network, the indoor temperature 7.0) of the thermal load building
H and the thermal reserve power Ry (1) are as shown in (b) to (d) of the Fig. 5;
[0108] and (3) performing park integrated energy reserve configuration by comprehensively considering the gas and thermal inertia reserves, a power generation side reserve and a demand side reserve,
[0109] wherein a gas and thermal inertia reserve mathematical model is as follows:
[0110] C,(t)=a, RE +b,R]
[0111] wherein, 4 and b, refer to a gas inertia reserve capacity price and a thermal . . . R° RŸ inertia reserve capacity price changed along with the time ? respectively, and and refer to the gas inertia input reserve capacity and the thermal inertia input reserve capacity changed along with the time ! respectively;
[0112] a power generation side reserve mathematical model is as follows:
[0113] C,(1)=cR+d,R}
[0114] wherein, “ and d, refer to a reserve capacity price and an electric quantity price . . . R° R° changed along with the time ’ after a reserve market is cleared, and and ‘* refer to the transaction capacity and the actually input reserve capacity, changed along with the time ’, of power generation side reserve respectively;
[0115] a demand side reserve mathematical model is as follows: k C(t) =} (eRP" +f/ RE")
[0116] =
[0117] wherein, €; and f; refer to a reserve capacity price and an electric quantity
RPA compensation price, changed along with the time /, of a demand side user i respectively, and Le refer to the demand response transaction capacity and the actually input reserve capacity, changed along with the time ’, of the demand side user i respectively, and k is the number of demand side users:
[0118] and the park integrated energy reserve model is built by comprehensively considering the gas and thermal inertia reserves, the power generation side reserve and the demand side reserve to minimize the total reserve purchase cost in a research time period under the premise of ensuring the system reliability level:
T minC(t) = > [C,(t) + C, (1) + C5 (1)]
[0119] I ;
[0120] The technical means disclosed in the scheme of the present invention not only are not limited to the technical means disclosed in the above embodiments, but also include the technical scheme composed of any combination of the above technical features. It should be pointed out that several improvements and modifications can be made without departing from the principle of the present invention for those of ordinary skill in the technical field, and these improvements and modifications are also regarded as the protective scope of the present invention.
Claims (5)
1. A method for configuring a park integrated energy system reserve with gas and thermal inertia reserves, characterized by comprising the following steps: (1) building an output model for a gas inertia reserve to cope with system power vacancy based on a natural gas pipeline transient model, specifically comprising: 1) building the natural gas pipeline transient model based on a continuity equation and a momentum equation of dynamic natural gas flow; 2) solving a natural gas pipeline tail end pressure response model based on the thought of finite element approximation, 3) building the output model for the gas inertia reserve to cope with the system power vacancy by considering the operation constraints of a natural gas system; (2) building an output model for a thermal inertia reserve to cope with the system power vacancy by considering the thermal time lag, thermal loss and thermal inertia characteristics of a thermodynamic system, specifically comprising: a) building a thermodynamic system model by comprehensively considering the thermal time lag, thermal loss and thermal inertia characteristics of the thermodynamic system; b) solving a thermal load building room temperature response model based on time-frequency domain transformation; c) building the output model for the thermal inertia reserve to cope with the system power vacancy by considering the operation constraints of the thermodynamic system; and (3) building a park integrated energy reserve model and configuring the park integrated energy reserve by comprehensively considering the gas and thermal inertia reserves, a power generation side reserve and a demand side reserve to minimize the total park integrated energy system reserve purchase cost.
2. The method for configuring the park integrated energy system reserve with the gas and thermal inertia reserves according to claim 1, characterized in that the step (1) specifically comprises the following steps: 1) the natural gas pipeline transient model is as follows:
op pv =0 ot Ox 2 2 opv | Dev „OP Av + pgsin0 =0 ot ox ox 2D wherein, #, Y and P refer to density, flow velocity and pressure of natural gas respectively, À, D and refer to friction coefficient of a pipeline, inner diameter of the pipeline and inclination angle between the pipeline and the horizontal plane respectively, & refers to gravitational acceleration, and Ÿ and ’ refer to time variable and space variable respectively; 2) solving the natural gas pipeline tail end pressure response model by utilizing Laplace transformation: AL , A’ AL ; A ; Adv bo (+= Pour (0 PO A, (OS Fou (0) R,T 2DR,,T L 2D wherein, 4, L and 7 refer to the cross sectional area, the length and the temperature of the natural gas pipeline respectively, Rır refers to a natural gas gas constant, Four 0. Pour (7) and Fou () refer to the natural gas pipeline tail end pressure changed along with the time 7, the first-order derivative and the second-order derivative of the natural gas pipeline tail end pressure . 0) f(t) SATE respectively, and “#%/ and “7 refer to the natural gas pipeline tail end flow changed along with the time ! and the first-order derivative of the natural gas pipeline tail end flow respectively; and 3) the output model for the gas inertia reserve to cope with the system power vacancy is as follows: 0 O<t<t,/>t RE | (0<t<t, 151) Gulf 7h) (t, <t<t,) the constraint conditions are as follows: ‘ou = OE wherein, Fou and Fou refer to an upper limit and a lower limit of Fou (1) respectively, hi and © refer to the start time and the cut-off time of gas inertia reserve supply respectively, Gyr is a natural gas thermal value, and Ji and À refer to the tail end flow of the pipeline at the time % and the time © respectively.
3. The method for configuring the park integrated energy system reserve with the gas and thermal inertia reserves according to claim 1, characterized in that the step (2) specifically comprises the following steps: 1) the thermodynamic system model is as follows: / Ta (0) =— v, (1) loss! = ul x1, 1,0 =[H,0-L, OV CM, L,0=e;>[, 01,0) wherein, 7.) refers to the thermal time lag, changed along with the time ’ of a transmission pipeline n corresponding to a thermal load building m, ', refers to the length of the transmission pipeline n, vu (1) refers to the hot water flow velocity, changed with the time 7, of
H the transmission pipeline n, loss, refers to the thermal loss power of the transmission pipeline n
H in a thermal transmission process, “» refers to the pipeline thermal loss rate of the transmission pipeline n, 1,0) and In () refer to the indoor temperature, changed with the time /, of the thermal load building m and the first-order derivative of the indoor temperature respectively, HC) refers to the heating power, changed with the time /, of a thermal network to the thermal load building m, Lu) refers to the heat loss power, changed with the time ?, of the thermal load building m, C4 refers to the indoor air specific heat capacity of the thermal load building, loss M refers to the indoor air quality of the thermal load building, “= refers to a heat dissipation coefficient of the thermal load building, Lou (1) represents the outdoor temperature, changed with the time /, of the thermal load building, ” =L 2,02 and Z refer to the total number of thermal load buildings in the park integrated energy system; 2) solving the thermal load building room temperature response model by utilizing Laplace
. LU500360 transformation: CM, (1) + eT, (A) = H, (A) + eT, . and 3) building the output model for the thermal inertia reserve to cope with the system power vacancy: 0 (O<t<t, t>1,) RY (f) = = = loss S RIM)=D €, (ne Tan) (t; <t<1) m=l m=l = Tue T,, ; wherein, ”“ refers to the normal comfortable temperature, ‘” is the lowest comfortable
H temperature, Ry (1) is thermal inertia reserve, changed along with the time 7, of the thermal load building m, and Land #4 refer to the start time and the cut-off time of thermal inertia reserve supply respectively.
4. The method for configuring the park integrated energy system reserve with the gas and thermal inertia reserves according to claim 1, characterized in that the park integrated energy system reserve model is as follows:
T minC (6) = H[C,) + C,() + C,(0] t=1 _ G H wherein, a gas and thermal inertia reserve mathematical model is Cl) =aR’ +b,R , a and b, refer to a gas inertia reserve capacity price and a thermal inertia reserve capacity price . . . RC RÏ Ce changed along with the time ! respectively, and * and refer to the gas inertia input reserve capacity and the thermal inertia input reserve capacity changed along with the time / . LL ; . C()=cR+dR respectively; a power generation side reserve mathematical model 1s (1) =¢, A , and d, refer to a reserve capacity price and an electric quantity price changed along with the time / . R° RS . .
after a reserve market is cleared, and ** refer to the transaction capacity and the actually input reserve capacity, changed along with the time ‘, of power generation side reserve Cyn) =D (e,R” +R) respectively; and a demand side reserve mathematical model is i=l , 0 and f refer to a reserve capacity price and an electric quantity compensation price, changed . . . RP RP} along with the time /, of a demand side user i respectively, and refer to the demand response transaction capacity and the actually input reserve capacity, changed along with the time ! _ of the demand side user i respectively, and k is the number of demand side users.
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 (1)
| Publication Number | Publication Date |
|---|---|
| LU500360B1 true LU500360B1 (en) | 2022-01-10 |
Family
ID=73339936
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| LU500360A LU500360B1 (en) | 2020-07-03 | 2020-08-06 | Method for configuring park integrated energy system reserve with gas and thermal inertia reserves |
Country Status (3)
| Country | Link |
|---|---|
| CN (1) | CN111950171B (en) |
| LU (1) | LU500360B1 (en) |
| WO (1) | WO2022000700A1 (en) |
Families Citing this family (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113256009B (en) * | 2021-05-31 | 2023-11-28 | 东南大学 | Two-stage multi-standby configuration method, system and device based on aerothermal inertia |
| CN113344357B (en) * | 2021-05-31 | 2023-10-03 | 东南大学 | Design method of comprehensive energy system based on frequency domain dynamic index |
| CN113283107A (en) * | 2021-06-10 | 2021-08-20 | 东南大学 | Evaluation method and model for inertial characteristics of gas-thermal system in comprehensive energy system |
| CN114218765A (en) * | 2021-11-26 | 2022-03-22 | 国网江苏省电力有限公司 | Comprehensive energy inertial support method for coping with power shortage of power grid |
| 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 |
| CN114723221A (en) * | 2022-03-11 | 2022-07-08 | 上海电力大学 | Thermoelectric combined optimization scheduling method for integral centralized heat supply and demand response |
| CN114881416B (en) * | 2022-04-08 | 2023-05-30 | 东南大学溧阳研究院 | Comprehensive energy inertial support method considering flexibility of cogeneration unit |
| CN114741988A (en) * | 2022-04-26 | 2022-07-12 | 东南大学 | Frequency modulation method considering gas-heat inertia of comprehensive energy system |
| CN114925501B (en) * | 2022-04-26 | 2024-10-01 | 东南大学 | Comprehensive energy system reliability assessment method considering gas-heat inertia |
| 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 |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU2010276368B2 (en) * | 2009-07-20 | 2016-02-25 | Schneider Electric It Corporation | Techniques for power analysis |
| US9244444B2 (en) * | 2011-03-07 | 2016-01-26 | Callida Energy Llc | Systems and methods for optimizing energy and resource management for building systems |
| CN108808663B (en) * | 2018-06-12 | 2020-04-28 | 浙江大学 | Industrial user heat demand response method based on multi-energy complementation |
| CN109376912B (en) * | 2018-09-29 | 2021-07-02 | 东南大学 | Operation optimization method of combined cooling heating and power system based on thermal inertia of civil building |
| CN109472050B (en) * | 2018-09-30 | 2023-10-24 | 东南大学 | Thermal inertia-based hybrid time scale scheduling method for cogeneration system |
| CN109474025B (en) * | 2018-10-08 | 2022-05-13 | 国网能源研究院有限公司 | Optimized dispatching model of park level comprehensive energy system |
| CN109978625B (en) * | 2019-03-28 | 2022-07-29 | 河海大学 | Comprehensive energy system multi-target operation optimization method considering electric and thermal gas network |
| CN110544017A (en) * | 2019-08-12 | 2019-12-06 | 上海交通大学 | Energy system reliability assessment method considering thermal inertia and energy network constraint |
| CN110503271B (en) * | 2019-08-30 | 2023-05-26 | 南京工业大学 | Multi-type energy storage configuration method of comprehensive energy system |
-
2020
- 2020-07-03 CN CN202010637501.6A patent/CN111950171B/en active Active
- 2020-08-06 WO PCT/CN2020/107349 patent/WO2022000700A1/en active Application Filing
- 2020-08-06 LU LU500360A patent/LU500360B1/en active IP Right Grant
Also Published As
| Publication number | Publication date |
|---|---|
| CN111950171A (en) | 2020-11-17 |
| WO2022000700A1 (en) | 2022-01-06 |
| CN111950171B (en) | 2022-03-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| LU500360B1 (en) | Method for configuring park integrated energy system reserve with gas and thermal inertia reserves | |
| WO2018130231A1 (en) | Integrated energy system optimization method based on heating grid and home thermal inertia | |
| CN107808216B (en) | Method for constructing comprehensive minimum optimization model for wind curtailment and electrical thermal load reduction of electricity-gas-heat interconnection system | |
| CN104813131A (en) | Thermal energy storage system comprising a combined heating and cooling machine and a method for using the thermal energy storage system | |
| Xu et al. | Dynamic simulation of thermal energy storage system of Badaling 1 MW solar power tower plant | |
| CN112613671B (en) | Microgrid distribution robust optimization scheduling method considering phase change water tank and building energy storage | |
| WO2022193422A1 (en) | Multi-device site selection method for comprehensive energy virtual power plant | |
| CN114218765A (en) | Comprehensive energy inertial support method for coping with power shortage of power grid | |
| Urbaneck et al. | Solar district heating in East Germany–transformation in a cogeneration dominated city | |
| Zhang et al. | Operation strategy and dynamic performance study of integrated solar combined-cycle system | |
| CN114294708B (en) | Method for adjusting heat storage of long-distance heat supply pipe network | |
| Wang et al. | Performance comparison of different combined heat and compressed air energy storage systems integrated with organic Rankine cycle | |
| CN102878611B (en) | Heating pipe network fine regulation system and method | |
| Li et al. | Study on the dynamic characteristics of a concentrated solar power plant with the supercritical CO2 Brayton cycle coupled with different thermal energy storage methods | |
| Dai et al. | Modeling and analysis of electrical heating system based on entransy dissipation-based thermal resistance theory | |
| CN110829470A (en) | Fused salt heat storage-based power grid peak regulation system | |
| GB2575858A (en) | Balanced energy network | |
| JP2023149312A (en) | Controller, underground heat utilization system, control method and program | |
| Chen et al. | Optimally coordinated operation of a multi-energy microgrid with coupled electrical and heat networks | |
| WO2020015241A1 (en) | Method for calculating interval power flow of thermal network taking into consideration heat transfer coefficient interval | |
| Kulakowski et al. | Design of a packed bed thermal storage unit for a solar system | |
| CN220038651U (en) | Ground source heat pump system capable of supplementing cold and heat | |
| CN109446714A (en) | Cool and thermal power transmitting energy transmission specific frictional resistance in energy internet determines method | |
| Miao et al. | Frequency Response Strategies Considering Gas-thermal Inertia in the Integrated Energy System | |
| CN109447450A (en) | The CHP system dispatching method of building comprehensive heating raising wind electricity digestion |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FG | Patent granted |
Effective date: 20220110 |