CN116228461A - P2P transaction-based multi-virtual power plant random optimization scheduling method - Google Patents

P2P transaction-based multi-virtual power plant random optimization scheduling method Download PDF

Info

Publication number
CN116228461A
CN116228461A CN202211551953.8A CN202211551953A CN116228461A CN 116228461 A CN116228461 A CN 116228461A CN 202211551953 A CN202211551953 A CN 202211551953A CN 116228461 A CN116228461 A CN 116228461A
Authority
CN
China
Prior art keywords
representing
virtual power
power plant
user
transaction
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
CN202211551953.8A
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.)
Yancheng Power Supply Co of State Grid Jiangsu Electric Power Co Ltd
Original Assignee
Yancheng Power Supply Co of State Grid Jiangsu Electric Power 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 Yancheng Power Supply Co of State Grid Jiangsu Electric Power Co Ltd filed Critical Yancheng Power Supply Co of State Grid Jiangsu Electric Power Co Ltd
Priority to CN202211551953.8A priority Critical patent/CN116228461A/en
Publication of CN116228461A publication Critical patent/CN116228461A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • 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
    • 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
    • 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/06315Needs-based resource requirements planning or analysis
    • 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/007Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources
    • H02J3/0075Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources for providing alternative feeding paths between load and source according to economic or energy efficiency considerations, e.g. economic dispatch
    • 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/008Circuit arrangements for ac mains or ac distribution networks involving trading of energy or energy transmission rights
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Marketing (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • Health & Medical Sciences (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Power Engineering (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • Primary Health Care (AREA)
  • General Health & Medical Sciences (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses a P2P transaction-based multi-virtual power plant random optimization scheduling method, which belongs to the field of power system energy optimization scheduling and comprises the steps of establishing a virtual power plant aggregation model comprising photovoltaics, energy storage, a central air conditioning system, a gas turbine and user loads; establishing a random optimization scheduling model considering photovoltaic uncertainty; adopting a Shapley value method to distribute the cooperation surplus of the multiple virtual power plants; and finally solving the optimal scheduling model by adopting GAMS software and outputting a result. The invention is beneficial to the active participation of the P2P transaction by the multiple virtual power plants, and can promote the overall operation income more than the traditional centralized transaction, thereby reasonably carrying out the P2P transaction and promoting the economic benefit and the environmental protection benefit of the arrangement.

Description

P2P transaction-based multi-virtual power plant random optimization scheduling method
Technical Field
The invention belongs to the technical field of power, relates to the aspect of power system power supply optimal scheduling, and discloses a multi-virtual power plant random optimal scheduling method based on P2P transaction.
Background
With the problems of large-scale access of distributed power sources, lack of timely information interaction among the distributed power sources and difficulty in scheduling among the distributed power sources, virtual power plants (virtual power plant, VPP) aggregate units such as Photovoltaic (PV), energy Storage (ES), central air conditioner (central air-conditioning system) and flexible load through advanced control, metering, communication and other technologies, and realize coordinated and optimized operation of a plurality of distributed energy sources through a higher level.
At present, the forms of electric energy sharing and distributed transaction between virtual power plants can be basically divided into the following two types: (1) Energy sharing and transaction mechanisms based on cooperative game mode; (2) Energy sharing and transaction mechanisms based on non-cooperative gaming. In the former, the main form is that each virtual power plant integrally participates in market trade in a alliance form, and the internal surplus electric energy is preferentially consumed in an energy sharing form. Domestic scholars put forward a cooperative game model based on single-to-many and multiple-to-multiple transaction conditions, and adopt a nucleolus method to realize cooperative surplus redistribution, and improve the nucleolus method by methods of clustering, calculating excess and the like, so that the method is still applicable to huge virtual power plant groups. Still scholars propose a day-ahead optimal scheduling strategy for building groups based on energy sharing, and redistribute the overall benefits after scheduling is completed by adopting a shape method. Unlike cooperative gaming, non-cooperative gaming-based power trading schemes focus on balancing under various principal interest conflicts, and can be categorized into platform-dominant and participant-dominant based trading schemes. At present, a market transaction model based on platform leading exists in China, wherein a Stackelberg game relationship is formed between a leader market operator and a follower, electric energy sharing is guided through electricity price, the near consumption of clean energy is promoted, and the transmission loss in the energy transaction process is reduced. The domestic scholars research a non-cooperative game model based on participant leading, establish a non-cooperative competition relationship among a plurality of sellers, and obtain an optimal bidding strategy of each electricity seller based on Nash equilibrium solution, wherein the game relationship comprises a master-slave game relationship between an agent and a virtual power plant and an evolution game relationship among a plurality of virtual power plants.
Disclosure of Invention
The invention provides a multi-virtual power plant random optimization scheduling method based on P2P transaction, which can better provide an optimal strategy scheme for a decision maker.
The invention particularly relates to a multi-virtual power plant random optimization scheduling method based on P2P transaction, which comprises the following steps:
step 1, constructing a virtual power plant aggregation model comprising photovoltaic, energy storage, a central air conditioning system, a gas turbine and user loads;
step 2, establishing a random optimization scheduling model considering photovoltaic uncertainty;
step 3, constructing a multi-virtual power plant cooperation residual distribution model based on a Shapley value method;
and step 4, solving the optimal scheduling model by adopting GAMS software and outputting a result.
The step 1 specifically comprises the following steps:
user load model:
the user consumes certain electric energy to bring certain profit, so that the profit generated by consuming energy by the user is:
U i (d i,t )=k i,t ln(1+d i,t ) (1)
wherein k is i,t Representing a combination of a weight coefficient and a preference parameter, reflecting the weight of the electricity utilization effect of the user and the energy consumption preference of the user, d i,t Representing the energy consumed by the user.
When there is a shortage between the load required by the user and the load that can be provided, the user will tend to be dissatisfied, so the dissatisfaction of the user is described by the cost of dissatisfaction of the user, as shown in the following formula:
Figure SMS_1
wherein beta is i Representing the priority coefficient of each virtual power plant, and beta i > 0. With a larger beta i Meaning that it is more concerned about comfort than economy. Omega t Is a weight factor for t-period dissatisfaction costs.
Figure SMS_2
Representing the ideal electricity demand of the user in the virtual power plant i during period t. As can be seen from the above, when the actual load is equal to the desired electricity demand, i.e.
Figure SMS_3
When (I)>
Figure SMS_4
The value of (2) is 0; when->
Figure SMS_5
When it is indicated that the user load demand is not met, and as d i,t Is reduced (1)>
Figure SMS_6
Is growing faster and faster; when->
Figure SMS_7
When it is indicated that the load demand has been met, and as the actual electricity load d is used i,t Increase of->
Figure SMS_8
The speed of the decrease of (1) starts to slow down, i.e. the user satisfaction reaches saturation.
And (3) energy storage model:
the battery capacity in the operation of the energy storage system may be expressed as:
Figure SMS_9
in the method, in the process of the invention,
Figure SMS_10
representing the state of charge of the stored energy, u c 、u d Charge and discharge power respectively representing energy storage capacity, +.>
Figure SMS_11
And->
Figure SMS_12
Respectively representing the maximum charge-discharge efficiency of the stored energy.
Photovoltaic model:
the predicted output power of the photovoltaic array, based on the solar radiation intensity, is:
Figure SMS_13
wherein eta is ct For the energy conversion efficiency of the photovoltaic array,
Figure SMS_14
is the total output of the photovoltaic unit at the moment t, S CA Is the photovoltaic array area. G t The predicted solar radiation intensity at a certain point in time t is indicated.
Central air conditioning model:
the invention adopts a modeling method based on cold and hot load calculation to describe a heat storage model of a public building, for the public building, a central air conditioning unit provides refrigerating capacity for the indoor during working, and meanwhile, the indoor and outdoor heat sources of the building have the heat release function and the heat storage function of the inner wall of the building so that the room temperature is continuously increased, and the specific expression is as follows:
Figure SMS_15
Figure SMS_16
/>
Figure SMS_17
Figure SMS_18
in which Q i,s,t Total cooling capacity provided for a central air conditioning system;
Figure SMS_19
the power consumption of the refrigerating unit in the refrigerating process and the power consumption of the cold storage tank in the cold releasing process are respectively shown. />
Figure SMS_20
The indoor temperature is represented, alpha, beta and gamma are respectively characteristic parameters of the intelligent building and an air conditioning system thereof, Δt is the time interval, < >>
Figure SMS_21
Is the output of the central air conditioning system, < >>
Figure SMS_22
To the capacity of the cold accumulation tank at the moment t, eta i,st 、η i,re The cold accumulation/release efficiency of the cold accumulation groove is respectively.
The step 2 specifically comprises the following steps:
the invention aims to establish a random optimization scheduling model considering uncertainty of the photovoltaic output, wherein randomness is expressed as uncertainty of the photovoltaic output, so that a scene s is adopted to represent actual output of the photovoltaic, and the maximum social benefit generated by virtual power plant operation is established. Thus, the objective function is:
Figure SMS_23
wherein u is P2P The objective function of the multiple virtual power plants when participating in the P2P transaction and the objective function of the multiple virtual power plants when not participating in the P2P transaction. Wherein N is i (i=1, 2, 3) represents the number of virtual power plants, s represents the photovoltaic scene, ρ s The probability of representing the scene s is indicated,
Figure SMS_24
represents profit of the user from consuming energy, < +.>
Figure SMS_25
Representing costs incurred by a virtual power plant in trade with a superordinate grid,/->
Figure SMS_26
Representing the cost of energy storage operation>
Figure SMS_27
Representing user dissatisfaction costs, +.>
Figure SMS_28
Representing the costs incurred by the operation of the gas turbine, +.>
Figure SMS_29
Representing the cost of the photovoltaic array generation. The specific calculation formula of each part is as follows:
1) Profit generated by consuming energy by users
U i (d i,t )=k i,t ln(1+d i,t ) (10)
Wherein k is i,t Representing a combination of a weight coefficient and a preference parameter, reflecting the weight of the electricity utilization effect of the user and the energy consumption preference of the user, d i,t Representing the energy consumed by the user.
2) Cost of electric market trade
Figure SMS_30
In the method, in the process of the invention,
Figure SMS_31
representing the price of electricity purchase->
Figure SMS_32
Representing the price of electricity selling, ->
Figure SMS_33
Representing the electricity purchase amount @ and @ of>
Figure SMS_34
Indicating the sales power.
3) Energy storage charge and discharge costs
Figure SMS_35
In the method, in the process of the invention,
Figure SMS_36
and->
Figure SMS_37
Respectively charging and discharging power of energy storage, +.>
Figure SMS_38
And->
Figure SMS_39
The cost coefficients of the charging and discharging of the energy storage are respectively.
4) User dissatisfaction cost
Figure SMS_40
Wherein beta is i Representing the priority coefficient of each virtual power plant, and beta i And is more than or equal to 0. With a larger beta i Meaning that it is more concerned about comfort than economy. Omega t Is a weight factor for t-period dissatisfaction costs.
Figure SMS_41
Representing the ideal electricity demand of the user in the virtual power plant i during period t. As can be seen from the above, when the actual load is equal to the desired electricity demand, i.e.
Figure SMS_42
When (I)>
Figure SMS_43
The value of (2) is 0; when->
Figure SMS_44
When it is indicated that the user load demand is not met, and as d i,t Is reduced (1)>
Figure SMS_45
And grow faster and faster. When->
Figure SMS_46
When it is indicated that the load demand has been met, and as the actual electricity load d is used i,t Increase of->
Figure SMS_47
The speed of the decrease of (1) starts to slow down, i.e. the user satisfaction reaches saturation.
5) Cost of operation of gas turbine
Figure SMS_48
Wherein a is MT 、b MT 、c MT A correlation coefficient representing the power generation cost of the gas turbine,
Figure SMS_49
indicating the gas turbine output.
6) Cost of photovoltaic output
The PV aggregated inside the virtual power plant is mainly a PV panel or roof PV, and its basic operation maintenance cost is:
Figure SMS_50
in the method, in the process of the invention,
Figure SMS_51
the power generated by the PV in the producer i at the moment t; />
Figure SMS_52
Maintenance costs for PV unit power generation in producer i.
7) P2P transaction cost
Figure SMS_53
Wherein, c ij Representing the trading coefficient between virtual power plant i and virtual power plant j.
The constraint conditions related to the invention are as follows:
1) Dissatisfaction function constraint
From the dissatisfaction function, the dissatisfaction function
Figure SMS_54
Should be relative to d i,t Is a monotonically decreasing function of (a). But->
Figure SMS_55
At this time, the dissatisfaction function begins to increment during this period. Thus (S)>
Figure SMS_56
Should be less than or equal to 2. Furthermore, in order to guarantee the basic and comfort requirements of the electrical load, the actual electrical load d i,t Should not be less than the ideal load +.>
Figure SMS_57
Half of (a) is provided. Thus (S)>
Figure SMS_58
The values of (2) must satisfy the following relationship:
Figure SMS_59
wherein d i,t Representing the load of the user and,
Figure SMS_60
indicating the desired power demand by the user.
2) Energy storage constraint
Figure SMS_61
Figure SMS_62
Figure SMS_63
Figure SMS_64
Figure SMS_65
/>
In the method, in the process of the invention,
Figure SMS_66
indicating the state of charge of the stored energy->
Figure SMS_67
Respectively representing upper and lower limits of the energy storage capacity;
Figure SMS_68
respectively representing the maximum charge and discharge power of the stored energy, < + >>
Figure SMS_69
Respectively representing the charge and discharge power of the stored energy at the time t.
3) Gas turbine constraints
Figure SMS_70
Figure SMS_71
Where ramp represents the rate of climb of the gas turbine,
Figure SMS_72
representing maximum and minimum output of the gas turbine, respectively,/->
Figure SMS_73
The output of the gas turbine at time t is shown.
4) Central air conditioning load constraint
Figure SMS_74
Figure SMS_75
Figure SMS_76
Figure SMS_77
In the method, in the process of the invention,
Figure SMS_78
respectively represents the power consumption of the refrigerating machine in the refrigerating process and the power consumption of the cold accumulation tank in the cold accumulation and release process, Q i ch,max Represents the maximum refrigerating capacity of the refrigerating unit, Q i st,max 、Q i,remax The maximum cold accumulation amount and the cold release amount of the cold accumulation groove are respectively +.>
Figure SMS_79
The capacity at the time t of the cold accumulation tank and the upper limit thereof are respectively set.
5) Power balance constraint
Figure SMS_80
In the method, in the process of the invention,
Figure SMS_83
representing photovoltaic output, ">
Figure SMS_84
Indicating gas turbine output,/->
Figure SMS_86
Representing the purchase power, < >>
Figure SMS_82
Represents the energy storage discharge power, P i,j,t Power representing P2P transactions, +.>
Figure SMS_85
Indicating the power of electricity selling, < >>
Figure SMS_87
Representing the stored charge power>
Figure SMS_88
Indicating central air conditioning power,/->
Figure SMS_81
Representing the load power.
The step 3 specifically comprises the following steps:
in fact, for the stability of cooperation among multiple virtual power plants, it is necessary to reasonably distribute the cooperation residuals of the multiple virtual power plants by using a certain profit sharing means. The invention adopts a Shapley value method as a residual distribution method for cooperation of multiple virtual power plants. The Shapley method has two main characteristics: 1. emphasis on fairness, i.e., the more virtual power plants that cooperate to produce utility contributions should distribute more revenue; 2. the importance of certain virtual power plants to the collaborative formation may be reflected. Without a particular virtual power plant, it is difficult to form a collaboration. Thus, it can help build up a collaborative virtual power plant that is allocated more revenue.
The benefit of the centralized transaction among the virtual power plants is shown as a formula (30), under the condition of sharing the energy of multiple virtual power plants, the overall operation benefit is increased, the excessive benefit becomes the emerging benefit, the actual benefit after the i-th virtual power plant is cooperated is shown as a formula (32),
f 0 =f 1 0 +f 2 0 +···+f i 0 (30)
f=f 1 +f 2 +···+f i (31)
Figure SMS_89
wherein f 0 Representing a virtual power plant total prior to participation in a P2P transactionIncome (E); f represents the total income of the virtual power plant after participating in P2P transaction; f (f) i 0 Representing the revenue of the virtual power plant i before participating in the P2P transaction; f (f) i Representing the benefits of the virtual power plant i for Shapley profit allocation after participation in the P2P transaction; s represents the number of virtual power plants in the set S; n represents the number of total virtual power plants; f (f) S\{i} Represents federated profits without virtual power plant i, (|S| -1) |! (n- |S|) is-! Representing the weighting factors.
The step 4 specifically comprises the following steps:
and solving an optimal scheduling model by adopting GAMS software, outputting a result, analyzing the advantage of P2P transaction relative to centralized transaction, and distributing profits by adopting a Shapley value method according to the contribution degree of different virtual power plants to P2P.
Drawings
FIG. 1 is a basic flow chart of the method of the present invention;
FIG. 2 is a sales of electricity to an upper grid by the virtual power plant 1;
FIG. 3 is a diagram illustrating the power trade of a virtual power plant with a superior grid under a centralized trade and a P2P trade;
FIG. 4 is a diagram illustrating power interactions between virtual power plants;
FIG. 5 is a diagram of the stored energy power situation in a virtual power plant;
FIG. 6 is a graph showing the energy of each unit in the virtual power plant 1 under a centralized transaction;
FIG. 7 is a graph showing the energy of each unit in the virtual power plant 1 under P2P transactions;
FIG. 8 is a diagram of profit sharing for each virtual power plant.
Detailed Description
Embodiments of the present invention will be described below with reference to the drawings.
As shown in fig. 1, the invention relates to a multi-virtual power plant random optimization scheduling method based on P2P transaction, which specifically comprises the following steps:
step 1, constructing a virtual power plant aggregation model comprising photovoltaic, energy storage, a central air conditioning system, a gas turbine and user loads;
step 2, establishing a random optimization scheduling model considering photovoltaic uncertainty;
step 3, constructing a multi-virtual power plant cooperation residual distribution model based on a Shapley value method;
and step 4, solving the optimal scheduling model by adopting GAMS software and outputting a result.
In order to verify that the energy sharing performed in the P2P transaction process can improve the social benefit, the effectiveness of the method provided by the simulation example verification model formed by three virtual power plants is designed. The virtual power plant 1 comprises a micro gas turbine, a central air conditioning system, photovoltaics, energy storage and a load, and the virtual power plants 2 and 3 comprise the central air conditioning system, the photovoltaics, the energy storage and the load. The electricity purchase price of the virtual power plant is according to the time-sharing electricity price issued by Jiangsu 2021 in 1 month, and at the time of peak (11:00-17:00, 22:00-24:00): 1.89 yuan/kWh, peak time (8:00-11:00, 17:00-22:00): 2.97 yuan/kWh, gu Shi (0:00-8:00): 0.85 yuan/kWh, and according to the current distributed power residual point internet surfing policy, the electricity purchasing price is kept unchanged to 0.85 yuan/kWh throughout the day. For better analysis of P2P trade conditions between virtual power plants, only 9:00-19:00 time periods are considered herein. In order to compare the overall operation modes and operation costs of multiple virtual power plants under different transaction mechanisms, the following comparative calculation examples are set: case1, P2P energy transaction is not carried out among the virtual power plants, and each virtual power plant only carries out centralized energy transaction with an upper power grid; case2, P2P energy sharing is carried out among the virtual power plants. The profit and cost of trading with three virtual power plants are shown in table 1, and it can be seen that the cost of trading with the higher grid for virtual power plants in case2 is reduced by 846.978 yuan compared to case1, however, because the trade between virtual power plants requires additional payment of service fees, the cost of dispersion increases by 24.327 yuan. Overall, the P2P transaction mechanism may reduce the operating cost by about 4.99% compared to the centralized transaction set.
Table 1 virtual Power plant operating amount under two examples
Figure SMS_90
Figure SMS_91
The electric energy trade situation between the virtual power plant 1 and the upper power grid under the two calculation examples is shown in fig. 2, and it can be seen from the figure that under the condition of centralized trade, the electric energy selling quantity of the virtual power plant 1 to the upper power grid is higher than that of the P2P trade situation. It can be seen that when P2P transactions are performed between virtual power plants, the multi-electric virtual power plant can flexibly select the transaction object, and enable clean energy and residual electric energy to be consumed as much as possible inside the alliance.
The market transaction electricity of the virtual power plant and the upper power grid under two examples is shown in fig. 3. The power interactions between the various virtual power plants are shown in FIG. 4, and each of the energy storage conditions in the virtual power plants, as well as the overall energy storage conditions, are shown in FIG. 5. By comprehensively analyzing the several graphs, it can be known that at t=10 to 11 hours, the electric quantity of the three virtual power plants is expressed as follows: the virtual power plant 1 is surplus in electric quantity, and the virtual power plants 2 and 3 are insufficient in electric quantity, and at this time, the virtual power plant 1 transmits electric energy to the virtual power plants 2 and 3 and stores the remaining electric energy in the energy storage. At this time, as shown in the figure, the energy storage in the virtual power plant 1 increases. At t=12 h, the three virtual power plant capacities appear as: the power of the virtual power plant 1 and the virtual power plant 3 are surplus and the power of the virtual power plant 2 is insufficient, at this time, the virtual power plant 1 and the virtual power plant 3 transmit electric power to the virtual power plant 2, and store the remaining energy in the energy storage in the units of the virtual power plant 1 and the virtual power plant 3, at this time, the energy storage capacities in the virtual power plant 1 and the virtual power plant 3 are increased to some extent as shown in fig. 5. At t=13 h, the power of the three virtual power plants appears as: the virtual power plant 1 and the virtual power plant 2 are full of power and the virtual power plant 3 is insufficient of power. At this time, the virtual power plant 1 and the virtual power plant 2 transmit surplus energy to the virtual power plant 3, but at this time, the electric energy shortage of the virtual power plant 3 cannot be completely balanced, and thus, the energy storage unit of the virtual power plant 3 releases part of the energy to ensure energy balance. At t=15-16 h, the power of the three virtual power plants appears as: the virtual power plant 1 is surplus and the virtual power plants 2 and 3 are deficient. At this time, the virtual power plant 1 and the virtual power plant 2 transmit electric power into the virtual power plant 3 and the electric power shortage of the virtual power plant 2 and the virtual power plant 3 cannot be completely balanced, and at this time, the electric power shortage is replenished by the energy storage released by the virtual power plant 2 and the virtual power plant 3. At t=17 to 18 hours, the electric quantity of the three virtual power plants represents surplus electric energy of the virtual power plant 1 and the virtual power plant 3, the electric energy of the virtual power plant 2 is insufficient, and the electric energy deficiency of the virtual power plant 2 is supplemented by the virtual power plant 3. At this time, the energy is released by the influence of the market price, and surplus electric energy of the virtual power plant 1 and the virtual power plant 3 is sold to the market to obtain the maximum profit.
As can be seen from fig. 6 and fig. 7, the energy management policies of the virtual power plant aggregation units in the two scenarios are basically consistent, but the energy storage scheduling policies of the virtual power plants in different transaction modes are slightly different, for example, the energy storage of the virtual power plant 1 is larger in the case of case2 in charge-discharge power, and smaller in the case of case1, which is mainly because in the case of case2, power shortage is invoked for balancing the power balance of the P2P transaction to perform energy balance management.
The operational benefits of each virtual power plant in case of case1 and case2 are shown in FIG. 8. Compared with a centralized transaction mode, under the condition that the virtual power plants participate in the P2P transaction, the respective operation profits can be improved to a certain extent under the condition that the Shapley value is carried out for profit allocation, profit allocation can be carried out according to the contribution value of each virtual power plant, fairness is guaranteed, and therefore each virtual power plant can actively participate in the P2P transaction.
It is noted that relational terms such as first and second, and the like are used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Moreover, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising one … …" does not exclude the presence of other like elements in a process, method, article, or apparatus that comprises the element.
Although embodiments of the present invention have been shown and described, it will be understood by those skilled in the art that various changes, modifications, substitutions and alterations can be made therein without departing from the principles and spirit of the invention, the scope of which is defined in the appended claims and their equivalents.

Claims (5)

1. A multi-virtual power plant random optimization scheduling method based on P2P transaction is characterized by comprising the following steps:
step 1, constructing a virtual power plant aggregation model comprising photovoltaic, energy storage, a central air conditioning system, a gas turbine and user loads;
step 2, establishing a random optimization scheduling model considering photovoltaic uncertainty;
step 3, constructing a multi-virtual power plant cooperation residual distribution model based on a Shapley value method;
and step 4, solving the optimal scheduling model by adopting GAMS software and outputting a result.
2. The P2P transaction-based multi-virtual power plant random optimization scheduling method according to claim 1, wherein the step 1 specifically includes:
constructing a user load model:
the user consumes certain electric energy to bring certain profit, so that the profit generated by consuming energy by the user is:
U i (d i,t )=k i,t ln(1+d i,t ) (1)
wherein k is i,t Representing a combination of a weight coefficient and a preference parameter, reflecting the weight of the electricity utilization effect of the user and the energy consumption preference of the user, d i,t Representing energy consumed by a user;
when there is a shortage between the load required by the user and the load that can be provided, the user will tend to be dissatisfied, so the dissatisfaction of the user is described by the cost of dissatisfaction of the user, as shown in the following formula:
Figure QLYQS_1
wherein beta is i Representing the priority coefficient of each virtual power plant, and beta i > 0, with larger beta i Meaning that it is more concerned about comfort than economy; omega t Is a weight factor for t period dissatisfaction costs;
Figure QLYQS_2
representing the ideal electricity demand of the user in the virtual power plant i in the period t; as can be seen from the above, when the actual load is equal to the desired electricity demand, i.e.
Figure QLYQS_3
When (I)>
Figure QLYQS_4
The value of (2) is 0; when->
Figure QLYQS_5
When it is indicated that the user load demand is not met, and as d i,t Is reduced (1)>
Figure QLYQS_6
Is growing faster and faster; when->
Figure QLYQS_7
When it is indicated that the load demand has been met, and as the actual electricity load d is used i,t Increase of->
Figure QLYQS_8
The speed of the decrease of (1) starts to slow down, i.e. the user satisfaction reaches saturation;
constructing an energy storage model:
the battery capacity in the operation of the energy storage system may be expressed as:
Figure QLYQS_9
in the method, in the process of the invention,
Figure QLYQS_10
representing the state of charge of the stored energy, u c 、u d Charge and discharge power respectively representing energy storage capacity, +.>
Figure QLYQS_11
And->
Figure QLYQS_12
Respectively representing the maximum charge and discharge efficiency of the stored energy;
building a photovoltaic model:
the predicted output power of the photovoltaic array, based on the solar radiation intensity, is:
P t pre =η ct ·S CA ·G t (4)
wherein eta is ct For photovoltaic array energy conversion efficiency, P t pre Is the total output of the photovoltaic unit at the moment t, S CA Is the photovoltaic array area; g t The predicted solar radiation intensity at a certain point in time t is indicated.
Constructing a central air conditioner model:
the regulation and control of the central air conditioning system should meet the indoor temperature requirement of building users, a modeling method based on cold and hot load calculation is adopted to describe a heat storage model of the public building, for the public building, the central air conditioning unit provides refrigerating capacity for the indoor space during the working period, meanwhile, the indoor and outdoor heat sources of the building have the heat release effect and the heat storage effect of the building inner wall so that the indoor temperature is continuously increased, and the specific expression is as follows:
Figure QLYQS_13
Figure QLYQS_14
Figure QLYQS_15
Figure QLYQS_16
in which Q i,s,t Total cooling capacity provided for a central air conditioning system;
Figure QLYQS_17
respectively representing the power consumption of the refrigerating unit in the refrigerating process and the power consumption of the cold storage tank in the cold releasing process; />
Figure QLYQS_18
The indoor temperature is represented, alpha, beta and gamma are respectively characteristic parameters of the intelligent building and an air conditioning system thereof, Δt is the time interval, < >>
Figure QLYQS_19
Is the output of the central air conditioning system, < >>
Figure QLYQS_20
To the capacity of the cold accumulation tank at the moment t, eta i,st 、η i,re The cold accumulation/release efficiency of the cold accumulation groove is respectively.
3. The P2P transaction-based multi-virtual power plant random optimization scheduling method according to claim 1, wherein the step 2 specifically comprises:
a random optimization scheduling model considering the uncertainty of the photovoltaic output is established, the randomness is expressed as the uncertainty of the photovoltaic output, so that the scene S is adopted to represent the actual output of the photovoltaic, the maximization of the social benefit generated by the operation of the virtual power plant is established, and the objective function is as follows:
Figure QLYQS_21
wherein u is P2P An objective function when the multiple virtual power plants participate in the P2P transaction and an objective function when the multiple virtual power plants do not participate in the P2P transaction; wherein N is i (i=1, 2, 3) represents the number of virtual power plants, s represents the photovoltaic scene, ρ s Representing the probability of scene s, U i (d i t ) Represents the profit of the user from consuming energy,
Figure QLYQS_22
representing the costs incurred by the virtual power plant in trade with the upper grid,
Figure QLYQS_23
representing the cost of energy storage operation>
Figure QLYQS_24
Representing user dissatisfaction costs, +.>
Figure QLYQS_25
Representing the costs incurred by the operation of the gas turbine, +.>
Figure QLYQS_26
Representing the cost of the photovoltaic array generation; the specific calculation formula of each part is as follows:
1) Profit generated by consuming energy by users
U i (d i,t )=k i,t ln(1+d i,t ) (10)
Wherein k is i,t Representing a combination of a weight coefficient and a preference parameter, reflecting the weight of the electricity utilization effect of the user and the energy consumption preference of the user, d i,t Representing energy consumed by a user;
2) Cost of electric market trade
Figure QLYQS_27
In the method, in the process of the invention,
Figure QLYQS_28
representing the price of electricity purchase->
Figure QLYQS_29
Representing the price of electricity selling, ->
Figure QLYQS_30
Representing the electricity purchase amount @ and @ of>
Figure QLYQS_31
Representing the sales power;
3) Energy storage charge and discharge costs
Figure QLYQS_32
In the method, in the process of the invention,
Figure QLYQS_33
and->
Figure QLYQS_34
Respectively charging and discharging power of energy storage, +.>
Figure QLYQS_35
And->
Figure QLYQS_36
Respectively the cost coefficients of energy storage charge and discharge; />
4) User dissatisfaction cost
Figure QLYQS_37
Wherein beta is i Representing the priority coefficient of each virtual power plant, and beta i Not less than 0, has larger beta i Virtual plant user intent of (3)Meaning that it is more concerned with comfort than economy; omega t Is a weight factor for t period dissatisfaction costs;
Figure QLYQS_38
representing the ideal electricity demand of the user in the virtual power plant i in the period t; as can be seen from the above, when the actual load is equal to the desired electricity demand, i.e.
Figure QLYQS_39
When (I)>
Figure QLYQS_40
The value of (2) is 0; when->
Figure QLYQS_41
When it is indicated that the user load demand is not met, and as d i,t Is reduced (1)>
Figure QLYQS_42
Is growing faster and faster; when->
Figure QLYQS_43
When it is indicated that the load demand has been met, and as the actual electricity load d is used i,t Increase of->
Figure QLYQS_44
The speed of the decrease of (1) starts to slow down, i.e. the user satisfaction reaches saturation;
5) Cost of operation of gas turbine
Figure QLYQS_45
Wherein a is MT 、b MT 、c MT A correlation coefficient representing the power generation cost of the gas turbine,
Figure QLYQS_46
representing the output of the gas turbine;
6) Cost of photovoltaic output
The PV aggregated inside the virtual power plant is mainly a PV panel or roof PV, and its basic operation maintenance cost is:
Figure QLYQS_47
in the method, in the process of the invention,
Figure QLYQS_48
the power generated by the PV in the producer i at the moment t; />
Figure QLYQS_49
Maintenance cost for the PV unit power generation in producer i;
7) P2P transaction cost
Figure QLYQS_50
Wherein, c ij Representing a trading coefficient between virtual power plant i and virtual power plant j;
the constraints involved are as follows:
1) Dissatisfaction function constraint
From the dissatisfaction function, the dissatisfaction function
Figure QLYQS_51
Is relative to d i,t But is monotonically decreasing in function of
Figure QLYQS_52
At this time, the dissatisfaction function begins to increment during this period; thus (S)>
Figure QLYQS_53
Should be less than or equal to 2; furthermore, in order to guarantee the basic and comfort requirements of the electrical load, the actual electrical load d i,t Should not be less than the ideal load +.>
Figure QLYQS_54
Half of (2); thus (S)>
Figure QLYQS_55
The values of (2) must satisfy the following relationship:
Figure QLYQS_56
wherein d i,t Representing the load of the user and,
Figure QLYQS_57
representing the ideal required power of the user;
2) Energy storage constraint
Figure QLYQS_58
Figure QLYQS_59
Figure QLYQS_60
Figure QLYQS_61
Figure QLYQS_62
In the method, in the process of the invention,
Figure QLYQS_63
indicating the state of charge of the stored energy->
Figure QLYQS_64
Respectively representing upper and lower limits of the energy storage capacity; />
Figure QLYQS_65
Respectively representing the maximum charge and discharge power of the stored energy, < + >>
Figure QLYQS_66
Respectively representing the charge and discharge power of the stored energy at the time t;
3) Gas turbine constraints
Figure QLYQS_67
Figure QLYQS_68
Where ramp represents the rate of climb of the gas turbine,
Figure QLYQS_69
representing the maximum and minimum output of the gas turbine respectively,
Figure QLYQS_70
the output of the gas turbine at the time t is shown;
4) Central air conditioning load constraint
Figure QLYQS_71
Figure QLYQS_72
Figure QLYQS_73
Figure QLYQS_74
In the method, in the process of the invention,
Figure QLYQS_75
respectively represents the power consumption of the refrigerating machine in the refrigerating process and the power consumption of the cold accumulation tank in the cold accumulation and release process, Q ich,max Represents the maximum refrigerating capacity of the refrigerating unit, Q ist,max 、Q i,remax The maximum cold accumulation amount and the cold release amount of the cold accumulation groove are respectively +.>
Figure QLYQS_76
The capacity at the moment t of the cold accumulation tank and the upper limit thereof are respectively;
5) Power balance constraint
Figure QLYQS_77
In the method, in the process of the invention,
Figure QLYQS_80
representing photovoltaic output, ">
Figure QLYQS_82
Indicating gas turbine output,/->
Figure QLYQS_84
Representing the purchase power, < >>
Figure QLYQS_79
Represents the energy storage discharge power, P i,j,t Power representing P2P transactions, +.>
Figure QLYQS_81
Indicating the power of electricity selling, < >>
Figure QLYQS_83
Representing the stored charge power>
Figure QLYQS_85
Indicating central air conditioning power,/->
Figure QLYQS_78
Representing the load power.
4. The P2P transaction-based multi-virtual power plant random optimization scheduling method according to claim 1, wherein the step 3 specifically comprises:
adopting a Shapley value method as a multi-virtual power plant cooperation residual distribution method; the benefit of the centralized transaction among the virtual power plants is shown as a formula (30), under the condition of sharing the energy of multiple virtual power plants, the overall operation benefit is increased, the excessive benefit becomes the emerging benefit, the actual benefit after the i-th virtual power plant is cooperated is shown as a formula (32),
Figure QLYQS_86
f=f 1 +f 2 +···+f i (31)
Figure QLYQS_87
wherein f 0 Representing the total revenue of the virtual power plant before participating in the P2P transaction; f represents the total income of the virtual power plant after participating in P2P transaction; f (f) i 0 Representing the revenue of the virtual power plant i before participating in the P2P transaction; f (f) i Representing the benefits of the virtual power plant i for Shapley profit allocation after participation in the P2P transaction; s represents the number of virtual power plants in the set S; n represents the number of total virtual power plants; f (f) S\{i} Represents federated profits without virtual power plant i, (|S| -1) |! (n- |S|) is-! Representing the weighting factors.
5. The P2P transaction-based multi-virtual power plant random optimization scheduling method according to claim 1, wherein the step 4 specifically comprises: and solving an optimal scheduling model by adopting GAMS software, outputting a result, and analyzing the advantage of the P2P transaction relative to the centralized transaction.
CN202211551953.8A 2022-12-05 2022-12-05 P2P transaction-based multi-virtual power plant random optimization scheduling method Pending CN116228461A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211551953.8A CN116228461A (en) 2022-12-05 2022-12-05 P2P transaction-based multi-virtual power plant random optimization scheduling method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211551953.8A CN116228461A (en) 2022-12-05 2022-12-05 P2P transaction-based multi-virtual power plant random optimization scheduling method

Publications (1)

Publication Number Publication Date
CN116228461A true CN116228461A (en) 2023-06-06

Family

ID=86581287

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211551953.8A Pending CN116228461A (en) 2022-12-05 2022-12-05 P2P transaction-based multi-virtual power plant random optimization scheduling method

Country Status (1)

Country Link
CN (1) CN116228461A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117151398A (en) * 2023-09-01 2023-12-01 深圳市科中云技术有限公司 Central air conditioner regulation and control method and system based on virtual power plant

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117151398A (en) * 2023-09-01 2023-12-01 深圳市科中云技术有限公司 Central air conditioner regulation and control method and system based on virtual power plant

Similar Documents

Publication Publication Date Title
CN109523052B (en) Virtual power plant optimal scheduling method considering demand response and carbon transaction
CN111008739B (en) Optimal regulation and control and income distribution method and system for cogeneration virtual power plant
CN109727158A (en) A kind of electric heating integrated energy system dispatching method based on the weak robust optimization of improvement
CN112308334A (en) Master-slave cooperation game-based multi-virtual power plant joint optimization scheduling method
CN108009693A (en) Grid-connected micro-capacitance sensor dual blank-holder based on two-stage demand response
CN114050571B (en) Comprehensive energy system energy hub control method considering carbon flow
CN113361875B (en) Optimization scheduling method for multi-microgrid comprehensive energy system considering demand side response and shared energy storage
CN109884888B (en) Multi-building micro-grid model prediction regulation and control method based on non-cooperative game
CN112398164A (en) Micro-energy-source network group optimization operation and cost distribution method containing shared energy storage system
CN112564102B (en) Multi-microgrid load optimization scheduling method and system
He et al. A new cooperation framework with a fair clearing scheme for energy storage sharing
CN113705906A (en) Energy coordination optimization operation method and system for comprehensive energy park
CN110245810A (en) A kind of Regional Energy internet coordination optimization operation method
CN115829142A (en) Industrial enterprise comprehensive energy system optimization planning method
CN116579560A (en) Virtual power plant energy distribution method considering multi-type load demand response
CN113746089A (en) Multi-user-oriented multi-time-scale power package and family energy optimization method
Hayati et al. A two-stage stochastic optimization scheduling approach for integrating renewable energy sources and deferrable demand in the spinning reserve market
CN116128201A (en) Multi-virtual power plant point-to-point energy trading method based on non-cooperative game
CN115423260A (en) Quantitative analysis method for new energy utilization of electric power market and policy service
CN116228461A (en) P2P transaction-based multi-virtual power plant random optimization scheduling method
CN114970962A (en) Optimization method of electric and thermal comprehensive energy system
CN112531785B (en) Photo-thermal-heat-storage active power distribution network main and auxiliary combined scheduling method under multiple time scales
CN115730747A (en) Multi-subject benefit distribution method of comprehensive energy system and application thereof
CN116780627A (en) Micro-grid regulation and control method in building park
CN112257951A (en) Optimized operation method of comprehensive energy system and power distribution company based on cooperative game

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