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 PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000013461 design Methods 0.000 title claims abstract description 15
- 238000005457 optimization Methods 0.000 title claims description 15
- 230000005611 electricity Effects 0.000 claims abstract description 20
- 230000005540 biological transmission Effects 0.000 claims abstract description 6
- 230000006870 function Effects 0.000 claims description 85
- 239000007789 gas Substances 0.000 claims description 29
- 238000004364 calculation method Methods 0.000 claims description 14
- 238000001816 cooling Methods 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 10
- 239000002918 waste heat Substances 0.000 claims description 10
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 8
- 238000010521 absorption reaction Methods 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 6
- 238000005265 energy consumption Methods 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 4
- 239000003345 natural gas Substances 0.000 claims description 4
- 238000005057 refrigeration Methods 0.000 claims description 3
- 238000013459 approach Methods 0.000 claims description 2
- 230000007423 decrease Effects 0.000 claims description 2
- 230000010354 integration Effects 0.000 claims description 2
- 230000008901 benefit Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 238000010438 heat treatment Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000002485 combustion reaction Methods 0.000 description 2
- AMXOYNBUYSYVKV-UHFFFAOYSA-M lithium bromide Chemical compound [Li+].[Br-] AMXOYNBUYSYVKV-UHFFFAOYSA-M 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000003546 flue gas Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008450 motivation Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06313—Resource planning in a project environment
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/003—Load forecast, e.g. methods or systems for forecasting future load demand
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-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
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:
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:
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:
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:
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:
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:
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.
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:
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:
wherein r is the annual capital investment rate, niThe service life of the equipment.
The power grid electricity purchasing cost is calculated as
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:
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:
in the formula (I), the compound is shown in the specification,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:
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.
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 ifx2∈I,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
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
As constraints on energy prices, namely:
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:
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:
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:
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:
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.
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)
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 |
-
2021
- 2021-11-26 CN CN202111418001.4A patent/CN114021385A/en active Pending
Cited By (1)
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 |