CN110119524A - A kind of energy resource system Multiple Time Scales interaction optimizing operation method and system - Google Patents

A kind of energy resource system Multiple Time Scales interaction optimizing operation method and system Download PDF

Info

Publication number
CN110119524A
CN110119524A CN201910144273.6A CN201910144273A CN110119524A CN 110119524 A CN110119524 A CN 110119524A CN 201910144273 A CN201910144273 A CN 201910144273A CN 110119524 A CN110119524 A CN 110119524A
Authority
CN
China
Prior art keywords
user
operator
model
load
moment
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201910144273.6A
Other languages
Chinese (zh)
Inventor
刘念
周琳洁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
North China Electric Power University
Original Assignee
North China Electric Power University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by North China Electric Power University filed Critical North China Electric Power University
Priority to CN201910144273.6A priority Critical patent/CN110119524A/en
Publication of CN110119524A publication Critical patent/CN110119524A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06Electricity, gas or water supply

Abstract

The present invention discloses a kind of energy resource system Multiple Time Scales interaction optimizing operation method and system.The optimizing operation method includes: the heating network established according to described for dsc data in practical heating system for thermal model;According to the delay time for determining transmission heat ratio same period transmission electricity for thermal model;Establish the electric load model and thermic load model of load side user;The model of operator is established according to the operating status of energy resource system;The benefit function for determining user determines the benefit function of the operator according to operator's model, obtains user's benefit function and operator's benefit function;Betting model is established using theory of games according to user's benefit function, operator's benefit function and delay time, the user and the corresponding game strategies of the maximum benefit of the operator are obtained according to betting model.Establishing betting model by theory of games can guarantee that operator and user respectively under the premise of maximum revenue, reach game equilibrium.

Description

A kind of energy resource system Multiple Time Scales interaction optimizing operation method and system
Technical field
The present invention relates to energy field, more particularly to a kind of energy resource system Multiple Time Scales interact optimizing operation method and System.
Background technique
Integrated energy system is able to achieve the scientific dispatch of various energy resources in system, realizes the efficient utilization of the energy, effectively slow The energy problem that solution faces at present, alleviating energy crisis.But electric, heat energy subsystem dynamic process time scale in system It differs larger, cannot achieve and optimize scheduling in same time scale.If only thermal energy transmission is idealized as without heat Delay, then there are gaps with actual conditions.In addition, industrial park includes a variety of energy main bodys, workload demand amount is big, but energy benefit With low efficiency, cause to the electricity shortage with energy main body.It is therefore desirable to interact to the integrated energy system of industrial park excellent Change, the enthusiasm of each demand response for participating in main body can not only be promoted, additionally it is possible to provide better energy supply clothes for user Business.
Summary of the invention
The object of the present invention is to provide a kind of more time rulers of the energy resource system that can be realized same time scale Optimized Operation Degree interaction optimizing operation method and system.
To achieve the above object, the present invention provides following schemes:
A kind of energy resource system Multiple Time Scales interaction optimizing operation method, the optimizing operation method include:
Acquire energy resource system for dsc data and power data;
According to the heating network in the practical heating system of dsc data foundation for thermal model;
According to the delay time for determining transmission heat ratio same period transmission electricity for thermal model;
According to the electric load model and thermic load model for establishing load side user for dsc data and the power data;
The model of operator is established according to the operating status of the energy resource system, obtains operator's model;
The benefit function that the user is determined according to the electric load model of the user, the thermic load model, according to institute The benefit function that operator's model determines the operator is stated, user's benefit function and operator's benefit function are obtained;
It is established according to user's benefit function, operator's benefit function and the delay time using theory of games Betting model obtains the user and the corresponding game strategies of the maximum benefit of the operator according to the betting model.
Optionally, the heating network established according to described for dsc data in practical heating system is specific for thermal model Include:
The temperature difference for obtaining carrier side and user side, obtains pressure and user's lateral pressure at the main pipe of the heating tube;
According to pressure at the main pipe and user's lateral pressure specific frictional resistance;
Wherein, PopeFor pressure at the main pipe, PuserFor user's lateral pressure, L is pipe range, αjFor jet chimney part The ratio of resistance and on-way resistance;
It tables look-up to obtain at steam outlet and user entry according to the pressure value of the carrier side and the user side Compare volume;
Steam averag density is calculated than volume according to described:
Wherein, ρopeFor density at main pipe, ρuserFor user side density;
The specific frictional resistance is standardized, standard specific frictional resistance is obtained
It is tabled look-up determining pipe diameter d according to the standard specific frictional resistance;
Steam mass flow:
Wherein, u is heating agent flow velocity, m/s;ρ is heating agent density, kg/m3;D is pipe diameter;
Temperature drop of the steam from exit to from user entry:
Wherein, qsplFor unit length heat supply pipeline heat waste, L is pipe range, and c is heating agent thermal capacitance;
Wherein, twFor insulating tube hull-skin temperature, tsFor tube circumference soil environment temperature, λgFor the thermal coefficient of soil, Hl For pipeline equivalent buried depth, Disl,outFor insulating layer outer diameter,
Wherein, twFor insulating tube hull-skin temperature, tsFor tube circumference soil environment temperature, t0For instrumentation tubes hull-skin temperature, λgFor the thermal coefficient of soil, λislFor the thermal coefficient for keeping the temperature layer material, D0For the pipe outside diameter that works, HlIt (m) is pipeline equivalent buried depth;
The temperature drop
Optionally, described to determine that the delay time of transmission heat ratio same period transmission electricity specifically wraps for thermal model according to described It includes:
The heat lag time
Wherein, KdelayFor heat lag coefficient;
Wherein, λislFor insulating layer thermal conductivity, λprtFor protective layer thermal conductivity, Disl,inFor insulating layer internal diameter, DprtFor protection Layer outer diameter;
Optionally, it is described according to it is described for dsc data and the power data establish load side user electric load model and Thermic load model specifically includes:
Electric load model
Calculate the sum of the net load that system net load of all users at the h moment is all users:
Wherein,Load for the user i in moment h, Pi hFor the photovoltaic predicted value at i-th family, n is that user is total Number,Net load for the user i in moment h, NLiFor the user i all the period of time at each moment net load set;
It whether is in the peak clipping period according to electricity consumption, electric load model of the user i in moment h is defined as follows:
Wherein, hex is the set of out-of-limit period,For the user i moment h fixation electric load,It is described User i moment h translatable load incrementss,Be the user i in moment h can reduction plans amount.
Thermic load model
Overall heat duty of the user i in moment h:
Wherein,For the user i the moment h fixation thermic load,It is the user i the moment h's Load can be increased.
Optionally, the operating status according to the energy resource system establishes the model of operator, obtains operator's model It specifically includes:
It is equal to the energy resource system in h-T in moment h total heat dutiesdelayWhen total quantity of heat production
The generated energy of the energy resource system of the moment h:
Wherein, ηchpFor generating efficiency, ηlossFor heat loss factor, δheatFor hot coefficient, θ is hotspot stress, and L is that natural gas is low Calorific value.
Optionally, the benefit that the user is determined according to electric load model, the thermic load model of the user Function determines that the benefit function of the operator specifically includes according to operator's model:
The utility models of user are the net profit model of user:
For with can utility function
Wherein, kiFor preference coefficient;
Energy cost is used in moment h for user i
Wherein,For the sale of electricity electricity price of operator, γ is unit heating power price;
The economic compensation obtained for user i in the participation peak clipping at h moment
Wherein,Indicate that user's peak clipping compensation electricity price for participating in interaction is given by operator, d indicates that operator gives user's Thermic load discount;
lsiFor translatable load, lciFor can reduction plans, Δ hiFor thermic load can be increased;
Operator's utility function model:
proM=EM-CM
For operator's income
Wherein,For moment h sale of electricity valence from operator to power grid,For moment h power purchase from operator to power grid Valence, γ are unit heating power price,For the sale of electricity electricity price of operator, LhFor in user's total load of moment h,For the energy Total power generation of the source system in moment h;
CMFor operator's income
Wherein, pgasFor Gas Prices,For translatable load of out-of-limit period,For the out-of-limit period can reduction plans,For peak clipping making up price,For more thermic load can be increased in limited time,Electricity price is compensated for peak clipping.
Optionally, described to be used according to user's benefit function, operator's benefit function and the delay time Theory of games establishes betting model, described to obtain the user and the maximum benefit pair of the operator according to the betting model The game strategies answered specifically include:
Using the operator as leader, the user is used as with follower;
The game strategies of the operator are that peak clipping compensates electricity price, the game strategies of the user be translatable electric load, It directly cuts down electric load and increases thermic load;
Initial peak clipping compensation electricity price is sent to the user by the operator;
The user formulates user's game strategies according to the initial peak clipping compensation electricity price and the benefit function, and User's game strategies are sent to the operator;
The operator calculates operator's benefit value according to user's game strategies and the benefit function;
Judge whether operator's benefit value and user's benefit value reach maximum value, if so, game reaches equal Weighing apparatus, the initial peak clipping compensation electricity price is optimal game strategies;Otherwise, the operator reformulates peak clipping compensation electricity price.
To achieve the goals above, the present invention also provides following schemes:
A kind of energy resource system Multiple Time Scales interaction optimization operating system, the operating system include:
Data acquisition module, for acquire energy resource system for dsc data and power data;
Heat supply model building module, for according to the confession for the heating network in the practical heating system of dsc data foundation Thermal model;
Delay time determining module, for according to the delay for determining transmission heat ratio same period transmission electricity for thermal model Time;
Load side model building module, for establishing load side user's for dsc data and the power data according to described Electric load model and thermic load model;
Operator's model building module is obtained for establishing the model of operator according to the operating status of the energy resource system Obtain operator's model;
Benefit function determining module, for according to the electric load model of the user, thermic load model determination The benefit function of user, determines the benefit function of the operator according to operator's model, obtain user's benefit function and Operator's benefit function;
Game strategies obtain module, for according to user's benefit function, operator's benefit function and described prolonging The slow time establishes betting model using theory of games, obtains the user according to the betting model and the operator is maximum The corresponding game strategies of benefit.
The specific embodiment provided according to the present invention, the invention discloses following technical effects: the invention discloses one kind Energy resource system Multiple Time Scales interact optimizing operation method and system, by calculating the heat lag time, and are dissolved into operation In the foundation of quotient's basic model and benefit function model, establish the model of Multiple Time Scales energy transportation process, realize by Heat lag takes into account in the operation interactive process of the optimization between user and operator.The solution of present invention application leader-followers games frame contains The interaction of heat lag optimizes operation problem, can guarantee that under the premise of maximum revenue, it is equal respectively to reach game by operator and user Weighing apparatus, acquires equilibrium solution, in gambling process, constantly interacts between operator and user, and dynamic generates the current generation most Dominant strategy, so that Optimization Solution is more efficient, feasibility is stronger, and convergence result more meets reality.
Detailed description of the invention
It in order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, below will be to institute in embodiment Attached drawing to be used is needed to be briefly described, it should be apparent that, the accompanying drawings in the following description is only some implementations of the invention Example, for those of ordinary skill in the art, without any creative labor, can also be according to these attached drawings Obtain other attached drawings.
Fig. 1 is the flow chart that energy resource system Multiple Time Scales provided by the invention interact optimizing operation method;
Fig. 2 is the composition block diagram of energy resource system Multiple Time Scales provided by the invention interaction optimization operating system.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other Embodiment shall fall within the protection scope of the present invention.
The object of the present invention is to provide a kind of more time rulers of the energy resource system that can be realized same time scale Optimized Operation Degree interaction optimizing operation method and system.
In order to make the foregoing objectives, features and advantages of the present invention clearer and more comprehensible, with reference to the accompanying drawing and specific real Applying mode, the present invention is described in further detail.
As shown in Figure 1, the present invention provides a kind of energy resource system Multiple Time Scales to interact optimizing operation method, the optimization Operation method includes:
Step 100: acquire energy resource system for dsc data and power data;
Step 200: according to the heating network in the practical heating system of dsc data foundation for thermal model;
Step 300: according to the delay time for determining transmission heat ratio same period transmission electricity for thermal model;
Step 400: according to the electric load model and heat for establishing load side user for dsc data and the power data Load model;
Step 500: establishing the model of operator according to the operating status of the energy resource system, obtain operator's model;
Step 600: the benefit letter of the user is determined according to the electric load model of the user, the thermic load model Number determines the benefit function of the operator according to operator's model, obtains user's benefit function and operator's benefit letter Number;
Step 700: game is used according to user's benefit function, operator's benefit function and the delay time Theory establishes betting model, obtains the user and the corresponding game of the maximum benefit of the operator according to the betting model Strategy.
The step 200: according to having for thermal model for the heating network established in practical heating system for dsc data Body includes:
The temperature difference for obtaining carrier side and user side, obtains pressure and user's lateral pressure at the main pipe of the heating tube;
According to pressure at the main pipe and user's lateral pressure specific frictional resistance;
Wherein, PopeFor pressure at the main pipe, PuserFor user's lateral pressure, L is pipe range, αjFor jet chimney part The ratio of resistance and on-way resistance;
It tables look-up to obtain at steam outlet and user entry according to the pressure value of the carrier side and the user side Compare volume;
Steam averag density is calculated than volume according to described:
Wherein, ρopeFor density at main pipe, ρuserFor user side density;
The specific frictional resistance is standardized, standard specific frictional resistance is obtained
It is tabled look-up determining pipe diameter d according to the standard specific frictional resistance;
Steam mass flow:
Wherein, u is heating agent flow velocity, m/s;ρ is heating agent density, kg/m3;D is pipe diameter;
Temperature drop of the steam from exit to from user entry:
Wherein, qsplFor unit length heat supply pipeline heat waste, L is pipe range, and c is heating agent thermal capacitance;
Wherein, twFor insulating tube hull-skin temperature, tsFor tube circumference soil environment temperature, λgFor the thermal coefficient of soil, Hl For pipeline equivalent buried depth, Disl,outFor insulating layer outer diameter,
Wherein, twFor insulating tube hull-skin temperature, tsFor tube circumference soil environment temperature, t0For instrumentation tubes hull-skin temperature, λgFor the thermal coefficient of soil, λislFor the thermal coefficient for keeping the temperature layer material, D0For the pipe outside diameter that works, HlIt (m) is pipeline equivalent buried depth;
The temperature drop
Described 300: determine that the delay time of transmission heat ratio same period transmission electricity specifically includes for thermal model according to described:
The heat lag time
Wherein, KdelayFor heat lag coefficient;
Wherein, λislFor insulating layer thermal conductivity, λprtFor protective layer thermal conductivity, Disl,inFor insulating layer internal diameter, DprtFor protection Layer outer diameter;
The step 400: according to the electric load model for establishing load side user for dsc data and the power data It is specifically included with thermic load model:
Electric load model
Calculate the sum of the net load that system net load of all users at the h moment is all users:
Wherein,Load for the user i in moment h, Pi hFor the photovoltaic predicted value at i-th family, n is that user is total Number,Net load for the user i in moment h, NLiFor the user i all the period of time at each moment net load set;
It whether is in the peak clipping period according to electricity consumption, electric load model of the user i in moment h is defined as follows:
Wherein, hex is the set of out-of-limit period,For the user i moment h fixation electric load,It is described User i moment h translatable load incrementss,Be the user i in moment h can reduction plans amount;
Thermic load model
Overall heat duty of the user i in moment h:
Wherein,For the user i the moment h fixation thermic load,It is the user i in the moment h Increase load.
The step 500: establishing the model of operator according to the operating status of the energy resource system, obtains operator's model It specifically includes:
It is equal to the energy resource system in h-T in moment h total heat dutiesdelayWhen total quantity of heat production
The generated energy of the energy resource system of the moment h:
Wherein, ηchpFor generating efficiency, ηlossFor heat loss factor, δheatFor hot coefficient, θ is hotspot stress, and L is that natural gas is low Calorific value.
The step 600: the benefit of the user is determined according to the electric load model of the user, the thermic load model Function determines that the benefit function of the operator specifically includes according to operator's model:
The utility models of user are the net profit model of user:
For with can utility function
Wherein, kiFor preference coefficient;
Energy cost is used in moment h for user i
Wherein,For the sale of electricity electricity price of operator, γ is unit heating power price;
The economic compensation obtained for user i in the participation peak clipping at h moment
Wherein,Indicate that user's peak clipping compensation electricity price for participating in interaction is given by operator, d indicates that operator gives user's Thermic load discount;
lsiFor translatable load, lciFor can reduction plans, Δ hiFor thermic load can be increased;
Operator's utility function model:
proM=EM-CM
For operator's income
Wherein,For moment h sale of electricity valence from operator to power grid,For moment h power purchase from operator to power grid Valence, γ are unit heating power price,For the sale of electricity electricity price of operator, LhFor in user's total load of moment h,For the energy Total power generation of the source system in moment h;
CMFor operator's income
Wherein, pgasFor Gas Prices,For translatable load of out-of-limit period,For the out-of-limit period can reduction plans,For peak clipping making up price,For more thermic load can be increased in limited time,Electricity price is compensated for peak clipping.
The step 700: it is used according to user's benefit function, operator's benefit function and the delay time Theory of games establishes betting model, described to obtain the user and the maximum benefit pair of the operator according to the betting model The game strategies answered specifically include:
Using the operator as leader, the user is used as with follower;
The game strategies of the operator are that peak clipping compensates electricity price, the game strategies of the user be translatable electric load, It directly cuts down electric load and increases thermic load;
Initial peak clipping compensation electricity price is sent to the user by the operator;
The user formulates user's game strategies according to the initial peak clipping compensation electricity price and the benefit function, and User's game strategies are sent to the operator;
The operator calculates operator's benefit value according to user's game strategies and the benefit function;
Judge whether operator's benefit value and user's benefit value reach maximum value, if so, game reaches equal Weighing apparatus, the initial peak clipping compensation electricity price is optimal game strategies;Otherwise, the operator reformulates peak clipping compensation electricity price.
As shown in Fig. 2, the present invention also provides a kind of energy resource system Multiple Time Scales to interact optimization operating system, the fortune Row system includes:
Data acquisition module 1, for acquire energy resource system for dsc data and power data;
Heat supply model building module 2, for according to the heating network in the practical heating system of dsc data foundation For thermal model;
Delay time determining module 3, for determining prolonging for transmission heat ratio same period transmission electricity for thermal model according to described The slow time;
Load side model building module 4, for establishing load side user for dsc data and the power data according to described Electric load model and thermic load model;
Operator's model building module 5 is obtained for establishing the model of operator according to the operating status of the energy resource system Obtain operator's model;
Benefit function determining module 6, for according to the electric load model of the user, thermic load model determination The benefit function of user, determines the benefit function of the operator according to operator's model, obtain user's benefit function and Operator's benefit function;
Game strategies obtain module 7, for according to user's benefit function, operator's benefit function and described prolonging The slow time establishes betting model using theory of games, obtains the user according to the betting model and the operator is maximum The corresponding game strategies of benefit.
Beneficial effects of the present invention:
By using interaction optimizing operation method of the present invention, the result after being optimized is applied to actual industrial Garden integrated energy system considers in the Optimized Operation a few days ago of heat lag effect.Based on interaction optimization method in the present invention Basic load data include electric load, thermic load (steam load), meet requirement of the industrial park for load type;It participates in excellent The interaction both sides of change process are garden operator and user, meet the development actual conditions of industrial park;Consider that thermal energy is transmitting In the process since the characteristic of heat transfer medium causes transmission time to be slower than the electric energy of same time tranfer, thus caused by heat lag imitate It answers, i.e. thermal energy and electric energy transmission is asynchronous, and optimum results more meet actual production requirement.The target of interaction optimization is for one day In, customer charge is limited beyond the peak load of power grid defined, and user receives the peak clipping compensation of operator, by adjusting that can put down Load is moved, it can the reduction of reduction plans and increase thermic load realization for out-of-limit period peak load.By using institute of the present invention Optimization method is stated, optimum results are applied to the industrial park integrated energy system Optimized Operation a few days ago for considering heat lag effect In, ENERGY PLANNING, reasonable arrangement can be carried out in advance according to demand history data under the premise of guaranteeing that meeting customer charge requires User's energy, and guarantee that load within the KB limit that power grid is set, evades the risk of load moon line in advance, is conducive to electricity Net is safely and reliably run.
Each embodiment in this specification is described in a progressive manner, the highlights of each of the examples are with other The difference of embodiment, the same or similar parts in each embodiment may refer to each other.For system disclosed in embodiment For, since it is corresponded to the methods disclosed in the examples, so being described relatively simple, related place is said referring to method part It is bright.
Used herein a specific example illustrates the principle and implementation of the invention, and above embodiments are said It is bright to be merely used to help understand method and its core concept of the invention;At the same time, for those skilled in the art, foundation Thought of the invention, there will be changes in the specific implementation manner and application range.In conclusion the content of the present specification is not It is interpreted as limitation of the present invention.

Claims (8)

1. a kind of energy resource system Multiple Time Scales interact optimizing operation method, which is characterized in that the optimizing operation method includes:
Acquire energy resource system for dsc data and power data;
According to the heating network in the practical heating system of dsc data foundation for thermal model;
According to the delay time for determining transmission heat ratio same period transmission electricity for thermal model;
According to the electric load model and thermic load model for establishing load side user for dsc data and the power data;
The model of operator is established according to the operating status of the energy resource system, obtains operator's model;
The benefit function that the user is determined according to the electric load model of the user, the thermic load model, according to the fortune Battalion's quotient module type determines the benefit function of the operator, obtains user's benefit function and operator's benefit function;
Game is established using theory of games according to user's benefit function, operator's benefit function and the delay time Model obtains the user and the corresponding game strategies of the maximum benefit of the operator according to the betting model.
2. a kind of energy resource system Multiple Time Scales according to claim 1 interact optimizing operation method, which is characterized in that institute State specifically including for thermal model according to the heating network established in practical heating system for dsc data:
The temperature difference for obtaining carrier side and user side, obtains pressure and user's lateral pressure at the main pipe of the heating tube;
According to pressure at the main pipe and user's lateral pressure specific frictional resistance;
Wherein, PopeFor pressure at the main pipe, PuserFor user's lateral pressure, L is pipe range, αjFor jet chimney local resistance With the ratio of on-way resistance;
Table look-up to obtain the ratio body at steam outlet and user entry according to the pressure value of the carrier side and the user side Product;
Steam averag density is calculated than volume according to described:
Wherein, ρopeFor density at main pipe, ρuserFor user side density;
The specific frictional resistance is standardized, standard specific frictional resistance is obtained
It is tabled look-up determining pipe diameter d according to the standard specific frictional resistance;
Steam mass flow:
Wherein, u is heating agent flow velocity, m/s;ρ is heating agent density, kg/m3;D is pipe diameter;
Temperature drop of the steam from exit to from user entry:
Wherein, qsplFor unit length heat supply pipeline heat waste, L is pipe range, and c is heating agent thermal capacitance;
Wherein, twFor insulating tube hull-skin temperature, tsFor tube circumference soil environment temperature, λgFor the thermal coefficient of soil, HlFor pipe Road equivalent buried depth, Disl,outFor insulating layer outer diameter,
Wherein, twFor insulating tube hull-skin temperature, tsFor tube circumference soil environment temperature, t0For instrumentation tubes hull-skin temperature, λgFor The thermal coefficient of soil, λislFor the thermal coefficient for keeping the temperature layer material, D0For the pipe outside diameter that works, HlIt (m) is pipeline equivalent buried depth;
The temperature drop
3. a kind of energy resource system Multiple Time Scales according to claim 1 interact optimizing operation method, which is characterized in that institute It states and determines that the delay time of transmission heat ratio same period transmission electricity specifically includes for thermal model according to described:
The heat lag time
Wherein, KdelayFor heat lag coefficient;
Wherein, λislFor insulating layer thermal conductivity, λprtFor protective layer thermal conductivity, Disl,inFor insulating layer internal diameter, DprtOutside for protective layer Diameter;
4. a kind of energy resource system Multiple Time Scales according to claim 1 interact optimizing operation method, which is characterized in that institute It states and establishes the electric load model of load side user for dsc data and the power data and thermic load model specifically wraps according to described It includes:
Electric load model
Calculate the sum of the net load that system net load of all users at the h moment is all users:
Wherein,Load for the user i in moment h, Pi hFor the photovoltaic predicted value at i-th family, n is total number of users, Net load for the user i in moment h, NLiFor the user i all the period of time at each moment net load set;
It whether is in the peak clipping period according to electricity consumption, electric load model of the user i in moment h is defined as follows:
Wherein, hex is the set of out-of-limit period,For the user i moment h fixation electric load,For the user i In the translatable load incrementss of moment h,Be the user i in moment h can reduction plans amount.
Thermic load model
Overall heat duty of the user i in moment h:
Wherein,For the user i the moment h fixation thermic load,For user i the increasing in the moment h Application of load.
5. a kind of energy resource system Multiple Time Scales according to claim 1 interact optimizing operation method, which is characterized in that institute The model for establishing operator according to the operating status of the energy resource system is stated, operator's model is obtained and specifically includes:
It is equal to the energy resource system in h-T in moment h total heat dutiesdelayWhen total quantity of heat production
The generated energy of the energy resource system of the moment h:
Wherein, ηchpFor generating efficiency, ηlossFor heat loss factor, δheatFor hot coefficient, θ is hotspot stress, and L is natural gas low-heat Value.
6. a kind of energy resource system Multiple Time Scales according to claim 1 interact optimizing operation method, which is characterized in that institute State the electric load model according to the user, the thermic load model determines the benefit function of the user, according to the operation Quotient module type determines that the benefit function of the operator specifically includes:
The utility models of user are the net profit model of user:
For with can utility function
Wherein, kiFor preference coefficient;
Energy cost is used in moment h for user i
Wherein,For the sale of electricity electricity price of operator, γ is unit heating power price;
The economic compensation obtained for user i in the participation peak clipping at h moment
Wherein,Indicate that user's peak clipping compensation electricity price for participating in interaction is given by operator, the heat that user gives in d expression operator is negative Lotus discount;
lsiFor translatable load, lciFor can reduction plans, Δ hiFor thermic load can be increased;
Operator's utility function model:
proM=EM-CM
For operator's income
Wherein,For moment h sale of electricity valence from operator to power grid,For the power purchase valence in moment h operator to power grid, γ For unit heating power price,For the sale of electricity electricity price of operator, LhFor in user's total load of moment h,For the energy system It unites in the total power generation of moment h;
CMFor operator's income
Wherein, pgasFor Gas Prices,For translatable load of out-of-limit period,For the out-of-limit period can reduction plans, For peak clipping making up price,For more thermic load can be increased in limited time,Electricity price is compensated for peak clipping.
7. a kind of energy resource system Multiple Time Scales according to claim 1 interact optimizing operation method, which is characterized in that institute It states and game mould is established using theory of games according to user's benefit function, operator's benefit function and the delay time Type, it is described specifically to be wrapped according to the betting model acquisition user and the corresponding game strategies of the maximum benefit of the operator It includes:
Using the operator as leader, the user is used as with follower;
The game strategies of the operator are that peak clipping compensates electricity price, and the game strategies of the user are translatable electric load, directly It cuts down electric load and increases thermic load;
Initial peak clipping compensation electricity price is sent to the user by the operator;
The user formulates user's game strategies according to the initial peak clipping compensation electricity price and the benefit function, and by institute It states user's game strategies and is sent to the operator;
The operator calculates operator's benefit value according to user's game strategies and the benefit function;
Judge whether operator's benefit value and user's benefit value reach maximum value, if so, game reaches balanced, institute Stating initial peak clipping compensation electricity price is optimal game strategies;Otherwise, the operator reformulates peak clipping compensation electricity price.
8. a kind of energy resource system Multiple Time Scales interaction optimization operating system, which is characterized in that the operating system includes:
Data acquisition module, for acquire energy resource system for dsc data and power data;
Heat supply model building module, for according to the heating network established for dsc data in practical heating system for hot-die Type;
Delay time determining module, when for according to the delay for determining transmission heat ratio same period transmission electricity for thermal model Between;
Load side model building module, for negative according to the electricity for establishing load side user for dsc data and the power data Lotus model and thermic load model;
Operator's model building module is transported for establishing the model of operator according to the operating status of the energy resource system Seek quotient module type;
Benefit function determining module, for determining the user according to electric load model, the thermic load model of the user Benefit function, determine the benefit function of the operator according to operator's model, obtain user's benefit function and operation Quotient's benefit function;
Game strategies obtain module, when for according to user's benefit function, operator's benefit function and the delay Between betting model established using theory of games, the user and the maximum benefit of the operator are obtained according to the betting model Corresponding game strategies.
CN201910144273.6A 2019-02-27 2019-02-27 A kind of energy resource system Multiple Time Scales interaction optimizing operation method and system Pending CN110119524A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910144273.6A CN110119524A (en) 2019-02-27 2019-02-27 A kind of energy resource system Multiple Time Scales interaction optimizing operation method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910144273.6A CN110119524A (en) 2019-02-27 2019-02-27 A kind of energy resource system Multiple Time Scales interaction optimizing operation method and system

Publications (1)

Publication Number Publication Date
CN110119524A true CN110119524A (en) 2019-08-13

Family

ID=67520556

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910144273.6A Pending CN110119524A (en) 2019-02-27 2019-02-27 A kind of energy resource system Multiple Time Scales interaction optimizing operation method and system

Country Status (1)

Country Link
CN (1) CN110119524A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111626587A (en) * 2020-05-21 2020-09-04 浙江大学 Comprehensive energy system topology optimization method considering energy flow delay characteristics

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070247331A1 (en) * 2006-03-31 2007-10-25 Bruce Angelis Integrated data collection, anomaly detection and investigation, such as integrated mobile utility meter reading, theft detection and investigation system
ES2444634T3 (en) * 2011-10-12 2014-02-26 Etogas Gmbh Procedure for manufacturing a methane-rich product gas, as well as proper installation for it
CN104408563A (en) * 2014-11-27 2015-03-11 国网宁夏电力公司 Regional planning method for wind power heating
US20150142368A1 (en) * 2013-05-31 2015-05-21 Patrick Andrew Shiel Method for determining mechanical heat-up lag (MHL) of a building from the building's natural thermal lag (NTL)
CN204648717U (en) * 2015-04-27 2015-09-16 东南大学 Anti-freeze solar water supply installation
CN106683001A (en) * 2016-12-06 2017-05-17 华北电力大学 Thermal power plant set identification data selection method based on historical operation data
CN108628176A (en) * 2018-06-29 2018-10-09 国网山东省电力公司电力科学研究院 It is a kind of meter and heat supply network combined heat and power Optimization Scheduling flexible
CN108764734A (en) * 2018-05-30 2018-11-06 南方电网科学研究院有限责任公司 A kind of integrated energy system interactive simulation method and device based on leader-followers games
US20180356105A1 (en) * 2017-04-28 2018-12-13 Southeast University Modeling Method of Combined Heat and Power Optimal Dispatching Model
CN109240366A (en) * 2018-08-10 2019-01-18 北京科技大学 A kind of hot activation building system equivalent outdoor temperature forecast Control Algorithm

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070247331A1 (en) * 2006-03-31 2007-10-25 Bruce Angelis Integrated data collection, anomaly detection and investigation, such as integrated mobile utility meter reading, theft detection and investigation system
ES2444634T3 (en) * 2011-10-12 2014-02-26 Etogas Gmbh Procedure for manufacturing a methane-rich product gas, as well as proper installation for it
US20150142368A1 (en) * 2013-05-31 2015-05-21 Patrick Andrew Shiel Method for determining mechanical heat-up lag (MHL) of a building from the building's natural thermal lag (NTL)
CN104408563A (en) * 2014-11-27 2015-03-11 国网宁夏电力公司 Regional planning method for wind power heating
CN204648717U (en) * 2015-04-27 2015-09-16 东南大学 Anti-freeze solar water supply installation
CN106683001A (en) * 2016-12-06 2017-05-17 华北电力大学 Thermal power plant set identification data selection method based on historical operation data
US20180356105A1 (en) * 2017-04-28 2018-12-13 Southeast University Modeling Method of Combined Heat and Power Optimal Dispatching Model
CN108764734A (en) * 2018-05-30 2018-11-06 南方电网科学研究院有限责任公司 A kind of integrated energy system interactive simulation method and device based on leader-followers games
CN108628176A (en) * 2018-06-29 2018-10-09 国网山东省电力公司电力科学研究院 It is a kind of meter and heat supply network combined heat and power Optimization Scheduling flexible
CN109240366A (en) * 2018-08-10 2019-01-18 北京科技大学 A kind of hot activation building system equivalent outdoor temperature forecast Control Algorithm

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
LIU NIAN .ECT: "Energy Sharing Provider for PV Prosumer Clusters A Hybrid Approach Using Stochastic Programming and Stackelberg Game", 《IEEE TRANSACTIONS ON INDUSTRIAL ELECTRONICS》 *
王炳坤: "《城市规划中的工程规划》", 31 December 1994, 天津大学出版社 *
荣爽: "促进供暖期风电消纳的多热源容量规划与协调调度策略", 《中国博士论文全文数据库-工程科技II辑》 *
马丽 等: "基于主从博弈策略的社区能源互联网分布式能量管理", 《电网技术》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111626587A (en) * 2020-05-21 2020-09-04 浙江大学 Comprehensive energy system topology optimization method considering energy flow delay characteristics
CN111626587B (en) * 2020-05-21 2022-09-06 浙江大学 Comprehensive energy system topology optimization method considering energy flow delay characteristics

Similar Documents

Publication Publication Date Title
WO2020155515A1 (en) Blockchain-based dual-source energy internet transaction method and device
CN105576710B (en) A kind of configuration method of distributed generation resource in integrated energy system
CN109727158A (en) A kind of electric heating integrated energy system dispatching method based on the weak robust optimization of improvement
CN110378729B (en) Comprehensive demand response method based on dynamic energy price strategy
CN110163411A (en) A kind of regional complex energy system operation optimization method
CN110163443A (en) Consider the micro- energy net Optimization Scheduling in the natural gas pressure regulating station of electric-gas integration requirement response
CN109740824A (en) A kind of multi-energy system Method for optimized planning considering heating network and thermic load
CN109787231A (en) A kind of integrated energy system distributed energy optimization method and system
CN109861302A (en) A kind of energy internet based on leader-followers games optimal control method a few days ago
CN114358432A (en) Multi-energy system optimization scheduling method and device considering demand response and carbon transaction
Yu et al. Optimal operation of the combined heat and power system equipped with power‐to‐heat devices for the improvement of wind energy utilization
CN113435095A (en) Method and system for optimizing scheduling of comprehensive energy system
Sun et al. How to implement real-time pricing in China? A solution based on power credit mechanism
Wang et al. Stackelberg game-based optimal scheduling of integrated energy systems considering differences in heat demand across multi-functional areas
CN110119524A (en) A kind of energy resource system Multiple Time Scales interaction optimizing operation method and system
Luo et al. Optimal operation and cost–benefit allocation for multi‐participant cooperation of integrated energy system
CN108491973A (en) A kind of electric-thermal association system Optimization Scheduling considering Demand Side Response
CN110245819A (en) A kind of energy efficiency evaluating method of multiple-energy-source interacted system
Wang et al. Flexibility-based improved model of multi-energy hubs using linear weighted sum algorithm
Zhao et al. Demand side management: A green way to power Beijing
CN108960564A (en) Supply of cooling, heating and electrical powers type microgrid distribution system dispatching method based on centralization modeling
CN109886584A (en) A kind of heat supply network dispatching method and device based on intelligent body
Yu et al. Energy hub model for energy management in energy internet
CN109002963A (en) Electric heating energy market equilibrium method for solving and device
CN115034551A (en) Energy interaction optimization method and device for multi-park comprehensive energy system

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