CN114021385A - Master-slave game-based optimization design method and device for regional comprehensive energy system - Google Patents

Master-slave game-based optimization design method and device for regional comprehensive energy system Download PDF

Info

Publication number
CN114021385A
CN114021385A CN202111418001.4A CN202111418001A CN114021385A CN 114021385 A CN114021385 A CN 114021385A CN 202111418001 A CN202111418001 A CN 202111418001A CN 114021385 A CN114021385 A CN 114021385A
Authority
CN
China
Prior art keywords
energy
user
load
function
game
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
CN202111418001.4A
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.)
State Grid Shanghai Electric Power Co Ltd
East China Power Test and Research Institute Co Ltd
Original Assignee
State Grid Shanghai Electric Power Co Ltd
East China Power Test and Research Institute Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Shanghai Electric Power Co Ltd, East China Power Test and Research Institute Co Ltd filed Critical State Grid Shanghai Electric Power Co Ltd
Priority to CN202111418001.4A priority Critical patent/CN114021385A/en
Publication of CN114021385A publication Critical patent/CN114021385A/en
Pending legal-status Critical Current

Links

Images

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
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06312Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
    • 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/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06313Resource planning in a project environment
    • 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/003Load forecast, e.g. methods or systems for forecasting future load demand
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/04Constraint-based CAD

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Marketing (AREA)
  • Tourism & Hospitality (AREA)
  • Game Theory and Decision Science (AREA)
  • Development Economics (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Educational Administration (AREA)
  • Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Power Engineering (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention relates to an optimal design method and device of a master-slave game-based regional comprehensive energy system, wherein the method comprises the following steps: constructing an energy transmission model taking a distributed cold, heat and electricity triple power supply system as a core, and setting energy balance conditions including electric balance conditions, heat balance conditions, cold balance conditions and prime mover output characteristic constraints; establishing a master-slave game model between a distributed energy station and energy users, wherein a game party comprises an energy supply side and a strategy set thereof, and a user side and a strategy set thereof, and a utility function comprises a utility function of the energy supply side and a utility function of the user side; the utility function of the energy supply side comprises the income of the energy station and the cost of the energy station; the utility function at the user side includes a user satisfaction function and a user's energy payment fee. Compared with the prior art, the method obtains the load and price configuration result of the comprehensively considered regional comprehensive energy system, and realizes comprehensive optimal configuration of the regional comprehensive energy system.

Description

Master-slave game-based optimization design method and device for regional comprehensive energy system
Technical Field
The invention relates to the field of regional integrated energy systems, in particular to a master-slave game-based optimization design method and device for a regional integrated energy system.
Background
The regional comprehensive energy system is a multi-element and three-dimensional energy network system formed by connecting links of energy production, transmission and distribution, storage, use and the like by taking distributed energy producers and consumers as nodes and taking electric, cold and hot energy networks as links. The regional integrated energy system can be regarded as the coupling and integration of a traditional single distributed energy supply system and a regional energy system which takes a regional heating and cooling system as a main body. The fundamental motivation of the method is to solve the problem of matching and balancing of multiple energy sources existing in the energy supply side and the energy utilization side in an area.
As a multi-energy combined supply system, complicated system coordination and energy regulation follow, so that the operation management of the regional comprehensive energy system is very important for realizing the comprehensive benefits of the system. The main purpose of the operation management research of the regional comprehensive energy system is to maximally mine the system income and the value on the premise of meeting the mutual matching of system energy supply and user requirements, so that the energy and the requirements need to be comprehensively considered, the overall optimization is realized, and the comprehensive consideration and the overall optimization cannot be realized by the conventional method.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provides a master-slave game-based regional comprehensive energy system optimization design method and device comprehensively considering energy and requirements
The purpose of the invention can be realized by the following technical scheme:
an optimal design method of a master-slave game-based regional comprehensive energy system comprises the following steps:
constructing an energy transmission model taking a distributed cold, heat and electricity triple power supply system as a core, and setting energy balance conditions including electric balance conditions, heat balance conditions, cold balance conditions and prime mover output characteristic constraints;
establishing a master-slave game model between a distributed energy station and an energy user, wherein a game party of the master-slave game model comprises an energy supply side and a strategy set thereof, and a user side and a strategy set thereof, and utility functions of the master-slave game model comprise utility functions of the energy supply side and utility functions of the user side;
the utility function of the energy supply side comprises the income of the energy station and the cost of the energy station;
the utility function of the user side comprises a user satisfaction function and the energy payment cost of the user.
Further, the solving process of the master-slave game model comprises the following steps:
s1: firstly, initializing energy prices of cold, heat and electricity, and substituting the energy prices into a utility function at a user side;
s2: obtaining load configuration under the condition that the utility function at the user side is optimal, substituting the optimal load configuration into the utility function at the energy supply side, and obtaining the optimal energy price at the moment;
s3: and repeating the step S2 until the deviation of the optimal energy price after two adjacent iterations is smaller than a preset deviation threshold value.
Further, in the load configuration solving process, a user-side actual load adjustment constraint is set, and a calculation expression of the user-side actual load adjustment constraint is as follows:
(dK,N,mh-lK,N,mh≤vKdK,N,mh)
in the formula (d)K,N,mhFor the hourly demand load of the user in the energy category K,/K,N,mhFor the actual energy load under energy category K, N is the user, m is the month, N is the hour, vKThe load adjustment ratio for the energy type K is adjusted.
Further, the calculation expression of the electrical balance condition is as follows:
Figure BDA0003376311860000021
wherein E is electric energy, N is user, m is month, N is hour, lB,N,m,hFor actual electrical load, Qgrid,m,hPurchasing electric power, Q, for the electric networkPM,m,h、QEC,m,hThe power generated by the prime mover and the power consumed by the electric refrigerator are respectively.
Further, the calculation expression of the thermal equilibrium condition is as follows:
Figure BDA0003376311860000022
in the formula IH,N,m,hFor actual heat load, Hh,PM,m,hFor heat supply from waste heat of prime mover GGB,m,hIs the gas consumption of the gas boiler, etaHB、ηGBThe efficiency of the heat exchanger and the efficiency of the gas boiler are respectively, and omega is the low calorific value of natural gas.
Further, the calculation expression of the cold balance condition is as follows:
Figure BDA0003376311860000023
in the formula IC,N,m,hIs made ofIntercalary cold load, HcPM,|m,hRefrigeration capacity of waste heat of prime mover, COPAC、COPECThe coefficients of performance of an absorption chiller and an electric chiller, respectively.
Further, the user satisfaction function satisfies:
when the actual load is smaller than the user demand load, the function takes a positive value to indicate that the user is not satisfied, and as the actual load approaches the demand load, the function value is reduced, and the dissatisfaction degree is reduced; when the actual load is larger than the demand load, the user satisfaction function is a negative value and represents user satisfaction, and as the actual load increases, the satisfaction function continuously decreases, the user satisfaction continuously increases, but the change amplitude gradually slows down and finally tends to be stable.
Further, the calculation expression of the user satisfaction function is as follows:
Figure BDA0003376311860000031
in the formula of UEU,NAs a function of user satisfaction, E is power, N is user, m is month, N is hour, lK,N,mhFor the actual energy load in energy class K, dK,N,mhFor the hourly demand load of the user in the energy category K, αKThe elastic parameter is the energy consumption of the user and is a negative value; beta is aKIs the energy price under normal load.
Further, the calculation expression of the cost function of the energy station is as follows:
CDES=Cinv+Cgrid+Cgas+Com
in the formula, CDESAs a function of the cost of the energy station, CinvConversion of annual capital costs for the plant, CgridPurchase of electricity for the grid, CgasFor gas charge, ComThe cost of operating and maintaining the equipment.
The invention also provides a master-slave game-based optimization design device for the regional comprehensive energy system, which comprises a memory and a processor, wherein the memory stores a computer program, and the processor calls the computer program to execute the steps of the method.
Compared with the prior art, the invention has the following advantages:
in the aspect of an energy system of a regional comprehensive energy system, an energy transmission model taking a distributed cold, heat and electricity triple power supply system as a core is constructed, and electric balance conditions, heat balance conditions, cold balance conditions and prime mover output characteristic constraints are considered; and then, an optimization solution is carried out by constructing a master-slave game model between the distributed energy station and the energy users, the aspects of energy supply and user requirements are comprehensively considered, and the load and price configuration result of the comprehensively considered regional comprehensive energy system is obtained by sequentially carrying out game optimization, so that the comprehensive optimization configuration of the regional comprehensive energy system is realized.
Drawings
Fig. 1 is a schematic structural diagram of a power supply and heat supply model provided in an embodiment of the present invention;
fig. 2 is a schematic structural diagram of a power supply and cooling model provided in an embodiment of the present invention;
FIG. 3 is a schematic structural diagram of a heating and cooling model according to an embodiment of the present invention;
fig. 4 is a schematic structural diagram of a cold, hot, and electricity triple power supply system of an energy station provided in an embodiment of the present invention;
fig. 5 is a schematic diagram of a two-side game on demand in a regional integrated energy system according to an embodiment of the present invention;
FIG. 6 is a schematic diagram of a user satisfaction function provided in an embodiment of the present invention;
fig. 7 is a flow chart of game model solving provided in the embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. The components of embodiments of the present invention generally described and illustrated in the figures herein may be arranged and designed in a wide variety of different configurations.
Thus, the following detailed description of the embodiments of the present invention, presented in the figures, is not intended to limit the scope of the invention, as claimed, but is merely representative of selected embodiments of the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
It should be noted that: like reference numbers and letters refer to like items in the following figures, and thus, once an item is defined in one figure, it need not be further defined and explained in subsequent figures.
Example 1
The embodiment provides an optimal design method of a master-slave game-based regional comprehensive energy system, which comprises the following steps:
constructing an energy transmission model taking a distributed cold, heat and electricity triple power supply system as a core, and setting energy balance conditions including electric balance conditions, heat balance conditions, cold balance conditions and prime mover output characteristic constraints;
establishing a master-slave game model between a distributed energy station and an energy user, wherein a game party of the master-slave game model comprises an energy supply side and a strategy set thereof, and a user side and a strategy set thereof, and utility functions of the master-slave game model comprise utility functions of the energy supply side and utility functions of the user side;
the utility function of the energy supply side comprises the income of the energy station and the cost of the energy station;
the utility function at the user side includes a user satisfaction function and a user's energy payment fee.
The process of the present invention is described in detail below
Physical modeling of regional integrated energy system
1) Power supply and heat supply model
The power supply and heat supply model mainly considers the areas with high heat load, generally adopts a gas turbine or
The gas internal combustion engine is a prime mover, and the auxiliary equipment is a waste heat boiler and a heat conversion device, as shown in fig. 1.
2) Power supply and cold supply model
The cooling load is generally provided by various refrigerators in the cooling model according to the demand of regional energy users, as shown in fig. 2.
3) Heat and cold supply model
The heating-cooling model generally uses a heat converter and a refrigerator as the output ends of a heat load and a cooling load, and can be equipped with a water source heat pump or an air source heat pump as an auxiliary energy supply device according to requirements, as shown in fig. 3.
In order to realize the cold-heat-electricity demand in the regional comprehensive energy system, the distributed energy system on the supply side region in the section takes a distributed cold-heat-electricity triple supply system (figure 4) as a core and mainly comprises a prime mover (such as a gas internal combustion engine, a gas turbine and the like), a gas boiler, an absorption refrigerator, an electric refrigerator and other equipment.
The prime motor generates electricity to provide electricity demand for energy users, and the insufficient electricity is supplemented by power purchasing of a power grid; the heat of the power generation waste heat of the prime motor (high-temperature flue gas, cylinder sleeve water waste heat and the like) is exchanged by the heat exchanger to supply heat (including heating and hot water requirements) for users, and the gas boiler is reserved; the lithium bromide absorption refrigerator and the electric refrigerator meet the cold load required by the user. Based on the energy supply structure, on the premise of constructing each unit equipment model and meeting the requirements of cold, hot and electric state loads of end users, the cold, hot and electric energy balance relation of the comprehensive energy system is established by taking the multi-energy flow as a link.
Energy balance
In the constructed physical model of the energy station, the energy demand of a user side needs to be met, and the balance of cold, heat and electric energy sources is established. The model established in the study divides a day into 24 periods, and each period is used as a basic unit for reasonable optimization. In addition, in order to increase the applicability range of the study, the operation optimization study was performed in the study of the present embodiment with a time span of one year.
1) Electric balance
In the whole system, assuming that the power supply is charged by the energy station, when the power station cannot generate power to meet the user demand, the power station purchases power from the power grid. The following equation gives the actual electrical load to the user on the left and the electrical load provided by the energy station on the right. The electric balance calculation formula is as follows:
Figure BDA0003376311860000051
wherein E is electric energy, N is user, m is month, N is hour, lE,N,m,hFor actual electrical load, Qgrid,m,hPurchasing electric power, Q, for the electric networkPM,m,h、QEC,m,hThe power generated by the prime mover and the power consumed by the electric refrigerator are respectively.
2) Heat balance
The following equation gives the actual heat load for the user on the left and the heat load provided by the energy station on the right. The heat load is provided by the gas boiler and the waste heat of the prime motor through the conversion of the heat exchanger. The heat balance calculation formula is:
Figure BDA0003376311860000061
in the formula IH,N,m,hFor actual heat load, Hh,PM,,m,hFor heat supply from waste heat of prime mover GGB,m,hIs the gas consumption of the gas boiler, etaHE、ηGBThe efficiency of the heat exchanger and the efficiency of the gas boiler are respectively, and omega is the low calorific value of natural gas.
3) Cold balance
The actual load of the cooling load is shown on the left side of the equation, and the cooling load output from the energy station is shown on the right side. The cooling load is provided by both the electric refrigerator and the absorption refrigerator.
Figure BDA0003376311860000062
In the formula IC,N,m,hFor actual cooling load, HcPM,|m,hRefrigeration capacity of waste heat of prime mover, COPAC、COPECAre respectively of absorption typeRefrigerator and electric refrigerator performance coefficients.
4) Prime mover output characteristics
The equipment capacity of the prime mover is required to be within the allowable range of the equipment, and in addition, the relationship between the heat output and the electric output is shown as follows:
0≤QPM,m,h≤QPM,max
HPM,m,h+HcPM,m,h=QPM,m,h·λ
in the formula, QPM,maxλ is the prime mover thermoelectric ratio, which is the rated capacity of the prime mover.
5) Other device output characteristics
The operation of other devices must be within an allowable range, and the respective output force is restricted by an upper limit and a lower limit:
0≤QEC,m,h≤QEC,max
0≤HcPM,m,h≤QAC,max
0≤ωGB,m,h≤QGB,max
in the formula, QEC,maxRated capacity, Q, of an electric refrigeratorAC,maxRated capacity, Q, for absorption chillersGB,maxIs rated capacity of the gas boiler.
Third, design and operation optimization model based on principal and subordinate game
3.1 base framework for Game optimization
The present example assumes that the regional integrated energy system under study consists of one distributed energy station and N individual users, thus forming a 1-N type energy supply architecture. And the supply and demand sides establish an interaction mechanism of energy flow and value flow based on the multi-energy demands of cold, heat, electricity and the like and corresponding prices, and finally establish the optimal design and operation scheme of the regional comprehensive energy system through game coordination. In consideration of the obvious primary and secondary relations at the supply and demand sides, a 1-N type primary and secondary game model between the distributed energy station and the energy users is established, and the multi-energy balance transaction problem between the energy station and the energy users is analyzed. Fig. 5 is a schematic diagram of a supply and demand two-side game in a regional integrated energy system. In fig. 5, the energy station is used as a leader of the game, the user is a follower, and the price of the cold, hot and electric multiple energy and the corresponding load are decision variables. Generally, a game includes several main factors, namely, a game party policy and a utility function, and the principal and subordinate game basic model proposed in this embodiment can be expressed as:
Γ={S∪Z,{WDES},{WEU}}
in the formula, S is a strategy set of an energy supply side and mainly comprises electricity, heat and cold prices, and Z is a strategy set of a user and mainly is the actual load of the user. WDESAs utility function of the energy supply side, WEUIs a utility function of the user side; by(s)*,z*) Indicating master and slave game balance points,(s)*,z*)∈S×Z。
As a leader of a master-slave game, the distributed energy station makes unit energy price on the basis of maximizing income, and the energy users respond to the strategy of the leader to realize utility maximization on the premise of considering the satisfaction degree of energy consumption. The final solution to the game model would be a Stackerlberg balance. In the equilibrium solution, the energy station side makes w through a strategy set SDESFunction maximization, and solving an objective function W in a strategy set Z by an energy userEUThe optimal solution of (1).
3.2 participant utility function
The concept of the utility function is initially applied to a function of the relationship between the utility obtained by a consumer in the consumption process and a consumed commodity in the economics so as to measure the satisfaction degree of the user in the consumption process, and with the popularization of energy marketization, the utility function is gradually used in the energy market to represent the benefit of two trading parties in the energy consumption process. In the game process, each participant has different satisfaction degrees on the game result, and different profits or payments are obtained, and the profits or payments adopt the concept of effectiveness. Higher utility means a higher degree of satisfaction achieved by the participants.
1) Utility function of supply side energy station
In the regional integrated energy system, a supply side distributed energy station sells cold, heat, electricity and other energy sources to make up investment and operation cost and obtain certain income, a utility function of the system mainly comprises two parts, and income of the energy station and cost of the energy station are specifically as shown in the following formula:
Max WDES=FDES-CDES
in the formula: w is aDESAs utility function of the energy station, FDESEnergy sales revenue for energy stations, CDESAs a function of the cost of the energy station.
(1) Revenue of energy station
Energy sales revenue F for energy stationDESThe expression is shown in the following formula and consists of total energy purchase payments of N users. Since the time span of this study is one year, the energy station revenue is expressed here as the product of the revenue per month and the number of days per month, as shown in the following equation.
FDES=∑∑∑(pC·lC,N,m,h+pH·lH,N,m,h+pE·lE,N,m,h)·mom
In the formula, pC、pH、pERespectively unit cold, hot, electric price, momDays corresponding to m months.
(2) Cost of energy station
The computational expression of the cost function of the energy station is as follows:
CDES=Ginv+Cgrid+Cgas+Com
in the formula, CDESAs a function of the cost of the energy station, CinvConversion of annual capital costs for the plant, CgridPurchase of electricity for the grid, CgasFor gas charge, ComThe cost of operating and maintaining the equipment.
The annual investment conversion cost of the equipment is as follows:
Figure BDA0003376311860000081
in the formula, Qi,maxIs the rated capacity of the device i, δiCRF, unit capacity of equipmentiFor the return on investment coefficient, the calculation formula is as follows:
Figure BDA0003376311860000082
wherein r is the annual capital investment rate, niThe service life of the equipment.
The power grid electricity purchasing cost is calculated as
Figure BDA0003376311860000083
In the formula, pgrid,mAnd purchasing electricity price for the power grid.
The gas cost is composed of two parts of natural gas consumed by a prime mover and a gas boiler:
Figure BDA0003376311860000084
in the formula etaPM,HFor the efficiency of the prime movergas1、pgas2The price of the gas for power generation and heat supply is respectively.
The calculation formula of the operation and maintenance cost of the equipment is as follows:
Figure BDA0003376311860000085
in the formula (I), the compound is shown in the specification,
Figure BDA0003376311860000086
for operating maintenance factor, Q, of the device ii,m,hIs the time-wise output power of device i.
2) Demand side user utility function
In order to further optimize the user experience, the present embodiment associates the degree of the actual load of the user deviating from the theoretical demand load with the satisfaction function of the user, and constructs the user satisfaction function at the user side as one of the important bases of the user decision.
The utility function of the regional energy user is the sum of the utility functions of all users, and the calculation formula is as follows:
Figure BDA0003376311860000091
in the formula, WEUFor the sum of all user utility functions of the region, WEU,NIs the utility function of user N.
The utility function of the energy user is the sum of the satisfaction function of the user and the energy purchase payment cost. Taking user N as an example, since the user utility function is required to take the maximum value, a negative sign is added before the formula, that is, the formula
WEU,N=-(UEU,N+FDES,N)
(1) Function of user satisfaction
UEU,NRepresents the satisfaction function of the user, which should have the following requirements: when the actual load is smaller than the user demand load, the function takes a positive value to indicate that the user is not satisfied, and as the actual load approaches the demand load, the function value is reduced, and the dissatisfaction degree is reduced; when the actual load is larger than the demand load, the user satisfaction function is a negative value and represents user satisfaction, and as the actual load increases, the satisfaction function continuously decreases, the user satisfaction continuously increases, but the change amplitude gradually slows down and finally tends to be stable.
As shown in fig. 6: when the actual load is smaller than the demand load, namely 1 < d, the user is dissatisfied, as shown in the left side of fig. 6, the satisfaction function of the user is a positive value, and the satisfaction function of the user is reduced as the actual load is closer to the demand load, which indicates that the dissatisfaction condition of the user is reduced. And secondly, when the actual load of the user is greater than the demand load, the satisfaction function of the user is a negative value, which indicates that the user is satisfied, and as the actual load increases, the satisfaction function of the user continues to decrease, and the satisfaction condition of the user increases, but the satisfaction function of the user does not always decrease, and tends to a fixed value when reaching a certain condition.
In this embodiment, the satisfaction function of the user is the sum of the satisfaction of the user for purchasing electricity, heat and cold, i.e.
Figure BDA0003376311860000092
In the formula of UEU,NAs a function of user satisfaction, E is power, N is user, m is month, N is hour, lK,N,mhFor the actual energy load in energy class K, dK,N,mhFor the hourly demand load of the user in the energy category K, αKThe elastic parameter is the energy consumption of the user and is a negative value; beta is aKIs the energy price under normal load.
Here, the energy demand flexibility refers to a degree of reaction of a relative change in energy demand to a relative change in price over a certain period of time. Research results show that the demand elasticity of the commodities for the user is related to the importance degree of the commodities, and the more critical commodities are demanded, the smaller the demand elasticity is, and the larger the demand elasticity is.
(2) Energy payment for a user
The energy payment cost of the user is the same as the energy selling cost of the energy station.
Fourth, model solution
4.1 Master-Slave Game Balancing Presence Attestation
In the principal and subordinate game theory, the solution of solving the equilibrium of the principal and subordinate game (Stackerlberg) is the final purpose of the game and is also the theoretical verification of the research result. The balance means that when all participants obtain a balance strategy, any participant cannot improve self-interest by only changing own strategy, namely, the balance strategy is a strategy that benefits the most for each rational participant under a certain environment.
When Nash equilibrium exists in the game model, the definition of Nash equilibrium indicates that(s)*,z*) The method is a balanced solution of a game model, and at the moment, the benefit of the energy supply side strategy S and the benefit of the user side strategy Z can reach the optimal value of Nash balance meaning. Let S, Z be an immediate subset of the metric space and Z be a non-empty convex set of the metric space, it can be seen that Stackelberg equilibrium exists when the master-slave game simultaneously satisfies the following conditions:
1)WEUis a continuity function with respect to policy set Z;
2)WEUis aboutK,N,m,hA quasi-convex function of (a);
3)wDESis a continuity function with respect to the set of policies S.
Firstly, proving the conditions (1) and (3), since the utility function of the energy station is calculated according to the difference between the energy selling yield and the cost of the energy station, the utility function of the user side is calculated in the above way, and the continuity of the two functions with respect to each variable is obvious. The presence of Stackerlberg equilibrium is proved, with an emphasis on whether the above condition (2) is satisfied. By definition "function f (x) is defined in the open interval I if
Figure BDA0003376311860000101
x2∈I,
Figure BDA0003376311860000102
Comprises the following steps:
f[rx1+(1-r)x2]≤rf(x1)+(1-r)f(x2) The f (x) interval I is called a convex function or a downward convex function. "by deductive calculation, WEUSatisfy the convex function characteristic in the definition, are aboutK,N,m,hA pseudo-concave function of (a).
In addition, as can be seen from the analysis, the energy users have unique optimal response to the energy price, and all strategy sets participating in the game are non-empty and compact, so that a balanced solution exists in the master-slave game of the energy transaction between the energy station and the users.
4.1 solving of Game model
The principal and subordinate game model constructed by the embodiment needs to reach a balanced result after multiple rounds of game play of game participants. Firstly, initializing the prices of cold, heat and electric energy sources at a supply side and substituting the prices into a utility function at a user side; then, load under the optimal condition of the user utility function is obtained, the optimal load is substituted into the utility function at the energy station side, and the optimal energy price at the moment is obtained; this process is repeated until a balanced solution for the optimal price is obtained. The solving flow chart is shown in fig. 7.
The load adjustment of the user must meet normal load requirements. In addition, in order to ensure that the user does not influence the reasonable energy demand of the user due to the low or high energy price, the user-side actual load adjustment is constrained as follows:
(dK,N,m,h-lK,N,m,h)≤vKdK,N,m,h
in the formula, vKThe load is adjusted in proportion, and the values are different according to different energy loads.
And the energy supply side determines energy pricing according to the load demand adjusted by the user and optimizes the income of the energy supply side. The user side objective function wEUObtaining the optimal load expression as
Figure BDA0003376311860000111
And substituting the optimal load into an energy station side utility function for optimizing to obtain the optimal energy price. Obtaining the relation between the energy price and the actual load at the same time
Figure BDA0003376311860000112
As constraints on energy prices, namely:
Figure BDA0003376311860000113
in the formula:
Figure BDA0003376311860000114
obtaining optimal load for user side
Figure BDA0003376311860000115
The optimal energy price.
The solving software applied in the embodiment is a Yalmip toolbox in Matlab, and is used for solving the interactive game between the supply and demand parties.
The embodiment also provides an optimal design device of the master-slave game-based regional integrated energy system, which comprises a memory and a processor, wherein the memory stores a computer program, and the processor calls the computer program to execute the steps of the optimal design method of the regional integrated energy system based on the master-slave game.
The foregoing detailed description of the preferred embodiments of the invention has been presented. It should be understood that numerous modifications and variations could be devised by those skilled in the art in light of the present teachings without departing from the inventive concepts. Therefore, the technical solutions available to those skilled in the art through logic analysis, reasoning and limited experiments based on the prior art according to the concept of the present invention should be within the scope of protection defined by the claims.

Claims (10)

1. An optimal design method of a master-slave game-based regional comprehensive energy system is characterized by comprising the following steps:
constructing an energy transmission model taking a distributed cold, heat and electricity triple power supply system as a core, and setting energy balance conditions including electric balance conditions, heat balance conditions, cold balance conditions and prime mover output characteristic constraints;
establishing a master-slave game model between a distributed energy station and an energy user, wherein a game party of the master-slave game model comprises an energy supply side and a strategy set thereof, and a user side and a strategy set thereof, and utility functions of the master-slave game model comprise utility functions of the energy supply side and utility functions of the user side;
the utility function of the energy supply side comprises the income of the energy station and the cost of the energy station;
the utility function of the user side comprises a user satisfaction function and the energy payment cost of the user.
2. The method for optimally designing the regional integrated energy system based on the master-slave game as claimed in claim 1, wherein the solving process of the master-slave game model comprises the following steps:
s1: firstly, initializing energy prices of cold, heat and electricity, and substituting the energy prices into a utility function at a user side;
s2: obtaining load configuration under the condition that the utility function at the user side is optimal, substituting the optimal load configuration into the utility function at the energy supply side, and obtaining the optimal energy price at the moment;
s3: and repeating the step S2 until the deviation of the optimal energy price after two adjacent iterations is smaller than a preset deviation threshold value.
3. The optimization design method of the regional integrated energy system based on the principal and subordinate game as claimed in claim 2, wherein in the load configuration solving process, the method further comprises setting a user-side actual load adjustment constraint, and a calculation expression of the user-side actual load adjustment constraint is as follows:
(dK,N,mh-lK,N,mh≤vKdK,N,mh)
in the formula (d)K,N,mhFor the hourly demand load of the user in the energy category K,/K,N,mhFor the actual energy load under energy category K, N is the user, m is the month, N is the hour, vKThe load adjustment ratio for the energy type K is adjusted.
4. The method for optimally designing the regional integrated energy system based on the principal and subordinate game as claimed in claim 1, wherein the computational expression of the electrical balance condition is as follows:
Figure FDA0003376311850000011
wherein E is electric energy, N is user, m is month, N is hour, lB,N,m,hFor actual electrical load, Qgrid,m,hPurchasing electric power, Q, for the electric networkPM,m,h、QEC,m,hThe power generated by the prime mover and the power consumed by the electric refrigerator are respectively.
5. The method for optimally designing the regional integrated energy system based on the principal and subordinate game as claimed in claim 1, wherein the computational expression of the thermal equilibrium condition is as follows:
Figure FDA0003376311850000021
in the formula IH,N,m,hFor actual heat load, Hh,PM,,m,hFor heat supply from waste heat of prime mover GGB,m,hIs the gas consumption of the gas boiler, etaHB、ηGBThe efficiency of the heat exchanger and the efficiency of the gas boiler are respectively, and omega is the low calorific value of natural gas.
6. The method for optimally designing the regional integrated energy system based on the principal and subordinate game as claimed in claim 1, wherein the calculation expression of the cold balance condition is as follows:
Figure FDA0003376311850000022
in the formula IC,N,m,hFor actual cooling load, HcPM,|m,hRefrigeration capacity of waste heat of prime mover, COPAC、COPECThe coefficients of performance of an absorption chiller and an electric chiller, respectively.
7. The method for optimally designing the regional integrated energy system based on the principal and subordinate game as claimed in claim 1, wherein the user satisfaction function satisfies the following conditions:
when the actual load is smaller than the user demand load, the function takes a positive value to indicate that the user is not satisfied, and as the actual load approaches the demand load, the function value is reduced, and the dissatisfaction degree is reduced; when the actual load is larger than the demand load, the user satisfaction function is a negative value and represents user satisfaction, and as the actual load increases, the satisfaction function continuously decreases, the user satisfaction continuously increases, but the change amplitude gradually slows down and finally tends to be stable.
8. The method for optimally designing the regional integrated energy system based on the principal and subordinate game as claimed in claim 1, wherein the computational expression of the user satisfaction function is as follows:
Figure FDA0003376311850000023
in the formula of UEU,NAs a function of user satisfaction, E is power, N is user, m is month, N is hour, lK,N,mhFor the actual energy load in energy class K, dK,N,mhFor the hourly demand load of the user in the energy category K, αKThe elastic parameter is the energy consumption of the user and is a negative value; beta is aKIs the energy price under normal load.
9. The method for optimally designing the regional integrated energy system based on the principal and subordinate game as claimed in claim 1, wherein the computational expression of the cost function of the energy station is as follows:
CDES=Cinv+Cgrid+Cgas+Com
in the formula, CDESAs a function of the cost of the energy station, CinvConversion of annual capital costs for the plant, CgridPurchase of electricity for the grid, CgasFor gas charge, ComThe cost of operating and maintaining the equipment.
10. An optimal design device of a master-slave game-based regional energy integration system, which is characterized by comprising a memory and a processor, wherein the memory stores a computer program, and the processor calls the computer program to execute the steps of the method according to any one of claims 1 to 9.
CN202111418001.4A 2021-11-26 2021-11-26 Master-slave game-based optimization design method and device for regional comprehensive energy system Pending CN114021385A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111418001.4A CN114021385A (en) 2021-11-26 2021-11-26 Master-slave game-based optimization design method and device for regional comprehensive energy system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111418001.4A CN114021385A (en) 2021-11-26 2021-11-26 Master-slave game-based optimization design method and device for regional comprehensive energy system

Publications (1)

Publication Number Publication Date
CN114021385A true CN114021385A (en) 2022-02-08

Family

ID=80066420

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111418001.4A Pending CN114021385A (en) 2021-11-26 2021-11-26 Master-slave game-based optimization design method and device for regional comprehensive energy system

Country Status (1)

Country Link
CN (1) CN114021385A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116227838A (en) * 2022-12-30 2023-06-06 国网江苏省电力有限公司电力科学研究院 Regional distributed resource collaborative optimization scheduling method and system

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116227838A (en) * 2022-12-30 2023-06-06 国网江苏省电力有限公司电力科学研究院 Regional distributed resource collaborative optimization scheduling method and system

Similar Documents

Publication Publication Date Title
Gu et al. Bi-level optimal low-carbon economic dispatch for an industrial park with consideration of multi-energy price incentives
Luo et al. Distributed peer-to-peer energy trading based on game theory in a community microgrid considering ownership complexity of distributed energy resources
Wang et al. Game theory-based multi-agent capacity optimization for integrated energy systems with compressed air energy storage
Qin et al. Optimal operation of integrated energy systems subject to coupled demand constraints of electricity and natural gas
Chen et al. Optimal day‐ahead scheduling of multiple integrated energy systems considering integrated demand response, cooperative game and virtual energy storage
Zhang et al. Low-carbon optimal operation of distributed energy systems in the context of electricity supply restriction and carbon tax policy: A fully decentralized energy dispatch strategy
CN111507529A (en) Economic-environment and energy dynamic pricing based double-layer game micro-energy network capacity optimization planning method
Zhang et al. Parametric optimization-based peer-to-peer energy trading among commercial buildings considering multiple energy conversion
CN112600253B (en) Park comprehensive energy collaborative optimization method and equipment based on optimal energy utilization efficiency
CN110489915B (en) Electric-heat combined scheduling method and system considering comprehensive demand response
CN115186902A (en) Regulating and controlling method, device, terminal and storage medium of greenhouse comprehensive energy system
CN112016825A (en) Centralized transaction optimization decision method of regional comprehensive energy system
CN115829142A (en) Industrial enterprise comprehensive energy system optimization planning method
Cao et al. Optimal operation of cold–heat–electricity multi-energy collaborative system based on price demand response
CN115081891A (en) Value-added service decision method for power consumer
CN112907274A (en) Park comprehensive energy trading price determination method based on Stackelberg double-layer game model
CN115204562A (en) Interconnected micro energy network distributed collaborative optimization scheduling method and system considering multi-energy sharing
CN116822697A (en) Comprehensive energy system low-carbon economic optimization method considering master-slave playing and demand response
CN117689234A (en) Multi-main-body double-layer game-based park comprehensive energy system scheduling method and system
CN114021385A (en) Master-slave game-based optimization design method and device for regional comprehensive energy system
Li et al. Multi-objective optimization of regional integrated energy system matrix modeling considering exergy analysis and user satisfaction
CN117829538A (en) Optimization method of energy hub network system considering fraud chess and multiple scenes
Wang et al. Capacity and operation joint optimization for integrated energy system based on Nash bargaining game
CN117522435A (en) Electric power market decision double-layer clearing model considering grading electricity price
CN111724026A (en) Optimization method for coupling operation of multi-energy network and water distribution network

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination