CN112529622A - Virtual power plant-based method for clearing multiple micro-micro main bodies participating in spot market - Google Patents
Virtual power plant-based method for clearing multiple micro-micro main bodies participating in spot market Download PDFInfo
- Publication number
- CN112529622A CN112529622A CN202011443256.1A CN202011443256A CN112529622A CN 112529622 A CN112529622 A CN 112529622A CN 202011443256 A CN202011443256 A CN 202011443256A CN 112529622 A CN112529622 A CN 112529622A
- Authority
- CN
- China
- Prior art keywords
- market
- time
- vpp
- representing
- day
- 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 41
- 238000004146 energy storage Methods 0.000 claims abstract description 61
- 230000005611 electricity Effects 0.000 claims abstract description 49
- 239000000126 substance Substances 0.000 claims description 33
- 238000010248 power generation Methods 0.000 claims description 24
- 210000002500 microbody Anatomy 0.000 claims description 11
- 230000002860 competitive effect Effects 0.000 claims description 10
- 230000006870 function Effects 0.000 claims description 10
- 230000005540 biological transmission Effects 0.000 claims description 9
- 238000005457 optimization Methods 0.000 claims description 8
- 230000009467 reduction Effects 0.000 claims description 7
- 238000003860 storage Methods 0.000 claims description 7
- 238000007599 discharging Methods 0.000 claims description 6
- 238000007781 pre-processing Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000008901 benefit Effects 0.000 abstract description 9
- 230000008569 process Effects 0.000 abstract description 8
- 230000000007 visual effect Effects 0.000 abstract description 7
- 238000004458 analytical method Methods 0.000 abstract description 5
- 239000003795 chemical substances by application Substances 0.000 description 12
- 238000010586 diagram Methods 0.000 description 10
- 230000007246 mechanism Effects 0.000 description 5
- 230000004044 response Effects 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000035515 penetration Effects 0.000 description 2
- 241001123248 Arma Species 0.000 description 1
- 238000012614 Monte-Carlo sampling Methods 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 230000009194 climbing Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 230000003862 health status Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000010348 incorporation Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000001502 supplementing effect Effects 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Images
Classifications
-
- 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
- G06Q30/00—Commerce
- G06Q30/02—Marketing; Price estimation or determination; Fundraising
- G06Q30/0201—Market modelling; Market analysis; Collecting market data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/90—Details of database functions independent of the retrieved data types
- G06F16/903—Querying
- G06F16/9038—Presentation of query results
-
- 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
- G06Q30/00—Commerce
- G06Q30/02—Marketing; Price estimation or determination; Fundraising
- G06Q30/0201—Market modelling; Market analysis; Collecting market data
- G06Q30/0206—Price or cost determination based on market factors
-
- 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
- G06Q30/00—Commerce
- G06Q30/02—Marketing; Price estimation or determination; Fundraising
- G06Q30/0283—Price estimation or determination
-
- 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
- G06Q30/00—Commerce
- G06Q30/06—Buying, selling or leasing transactions
- G06Q30/08—Auctions
-
- 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
-
- Y—GENERAL 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS 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
- Y04S50/00—Market activities related to the operation of systems integrating technologies related to power network operation or related to communication or information technologies
- Y04S50/14—Marketing, i.e. market research and analysis, surveying, promotions, advertising, buyer profiling, customer management or rewards
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Accounting & Taxation (AREA)
- Finance (AREA)
- Development Economics (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Game Theory and Decision Science (AREA)
- Data Mining & Analysis (AREA)
- Health & Medical Sciences (AREA)
- Databases & Information Systems (AREA)
- Water Supply & Treatment (AREA)
- Primary Health Care (AREA)
- Tourism & Hospitality (AREA)
- Human Resources & Organizations (AREA)
- General Health & Medical Sciences (AREA)
- Computational Linguistics (AREA)
- Public Health (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention provides a virtual power plant-based clearing method for a multi-element small micro-subject participating spot market, which solves the risk problem brought to VPP income by the uncertainty of a competitor in a trading decision game process. The method comprises the following steps: firstly, constructing a double-layer clearing model considering the uncertainty of the electricity price and a constraint condition of the double-layer clearing model based on the VPP parameters of the power grid side and the information of the electric energy buyer and the electric energy seller; secondly, solving the double-layer output model based on the constraint conditions of the double-layer output model to obtain the optimal bidding result and the optimal income of the day-ahead energy and duration frequency modulation market; and finally, obtaining the successful bidding condition, the optimal quotation and the clearing price of each main body of the market according to the optimal bidding result and the optimal income. According to the invention, through the two-stage double-layer clearing model and the two-stage double-layer constraint condition established by the clearing method, the visual benefit analysis and bidding decision of the peak shaving market in which the power grid side energy storage participates can be simply and conveniently realized by considering the uncertainty of the electricity price.
Description
Technical Field
The invention relates to the technical field of data processing, in particular to a virtual power plant-based method for clearing a plurality of small micro-subjects participating in spot market.
Background
Environmental issues have led to an increase in the penetration of renewable energy sources in several power systems around the world over the last few years. However, the integration of renewable energy power generation presents many challenges to the operation of power systems, particularly due to the limited predictability and high variability of the primary "fuels" (i.e., wind, solar radiation), which make renewable energy power plants non-dispatchable.
To date, a great deal of research has been conducted to find the best way to incorporate renewable energy power plants into existing power systems within the concept of smart grids. Among them, Demand Response (DR) and distributed energy storage (DS) have proven to be the most important key factors for smart grid deployment so far, as their incorporation into the demand side of the grid can alleviate the challenges of large-scale demand response penetration and improve system reliability. The use of Advanced Metering Infrastructure (AMI) can provide consumers with real-time electricity prices, can increase their active participation, and help to efficiently integrate intermittent resources into the grid. In the same direction, the fast response characteristic of distributed energy storage makes it a valuable regulatory resource, and can balance the negative effects of intermittently regenerated resources.
At present, a Virtual Power Plant (VPP) is one of the most effective and potential technologies recognized in the world for dealing with large-scale renewable energy grid connection. With the continuous promotion of the progress of the electric power market in China, the advantages of VPPs (virtual private Power Point) as independent bodies participating in the energy market and frequency modulation are widely accepted, but a complete market transaction mechanism and an effective model are not provided for exciting and standardizing the exertion of the energy storage peak regulation capability. The VPP is an important frequency modulation means as a high-quality bidirectional modulation resource with mature large-scale application technical conditions, can relieve the influence of the increase of the frequency modulation demand on the conventional unit, and reduces the loss of the conventional unit caused by frequent output reduction. However, even if the VPP is added to the fm auxiliary service, the mechanism problem of the fm auxiliary service market still cannot be solved, and the commercial application thereof is more dependent on the solution of the market mechanism problem.
In the market mechanism problem, the electricity price is the basic and core factor of the electric power market transaction, the electricity price is determined by the market, the rationalization and the validation of the electricity price in the transaction link can be guaranteed, and the electricity price can be guaranteed to correctly reflect the actual value of the electric power in the production and the transaction. The day-ahead electricity price prediction can provide decision basis for the bidding of power generators and electricity selling companies on spot market, and the accurate day-ahead electricity price prediction can greatly reduce the price risk of market participants and bring stable benefits to the market participants. In the current electric power market, the electricity purchase price is obtained by the trading center by integrating the quotation of the power generation party, the market demand and the condition of the electric power system. The electricity purchase price is determined by the value of the power commodity and is influenced by the supply and demand relation, unlike the past power control period. Therefore, the price of electricity should contain abundant information, such as the cost of the network unit, the supply and demand status of the market, and the health status of the system (including the transmission network and the generator). However, the price of electricity in the electricity market is uncertain due to its dynamic changes, resulting in higher bidding risk and failure to obtain stable profit. Therefore, the above factors need to be taken into account when formulating a storage bidding strategy. There is currently a few literature studying the bidding strategy of VPPs participating in the spot-market, and the influence of the randomness of self-quote and the uncertainty of competitor strategy in the bidding decision process on the bidding strategy.
Disclosure of Invention
Aiming at the technical problems that the risk brought by the uncertainty of a competitor in the process of introducing a multi-scenario description transaction decision game to the VPP income cannot ensure that the VPP reasonably participates in market-oriented transaction, the invention provides a clearing method of a virtual power plant-based multi-element small micro-subject participating spot market from the participation of the multi-element small micro-subject in the spot market, and provides a two-stage double-layer combined bidding model; firstly, in the first stage, a plurality of small micro bodies are aggregated into a VPP according to the clustering idea and are scheduled by a unified operation mechanism, and secondly, the whole VPP is preprocessed according to the analysis of the output characteristics of the plurality of small micro bodies, so that the upper limit and the lower limit of the whole VPP for charging and discharging can be determined; the second stage is that the VPP participates in a master-slave game model of a joint bidding strategy of the day-ahead energy and real-time balance market, and the goal of the master-slave game model is that the day-ahead profit maximization is combined with the minimization of the expected real-time production and consumption unbalance cost; the two-stage double-layer clearing model and the two-stage double-layer constraint condition established by the clearing method can simply and conveniently realize the visual benefit analysis and bidding decision of the peak shaving market in which the power grid side energy storage participates, considering the uncertainty of the electricity price.
The technical scheme of the invention is realized as follows:
a virtual power plant-based method for clearing a plurality of small micro-subjects participating in spot market includes the following steps:
the method comprises the following steps: constructing a double-layer clearing model considering the uncertainty of the electricity price based on the VPP parameters of the power grid side and the information of the electric energy buyer and the electric energy seller;
step two: constructing a constraint condition of a double-layer clearing model based on the VPP parameters of the power grid side, the information of both electric energy buyers and sellers and the double-layer clearing model;
step three: solving the double-layer clearing model based on the constraint conditions of the double-layer clearing model to obtain the optimal bidding result and the optimal profit of the day-ahead frequency modulation market of the day-ahead energy and duration;
step four: and obtaining the successful bidding condition, the optimal quotation and the clearing price of each main body of the market according to the optimal bidding result and the optimal income.
The double-layer clearing model comprises a first-stage model, a second-stage inner-layer model and a second-stage outer-layer model, wherein the first-stage model is a VPP coordination optimization scheduling model, the second-stage inner-layer model is a VPP participation day-ahead market income maximization model, and the second-stage outer-layer model is a full social welfare maximization model; the constraint conditions of the double-layer output model comprise the constraint conditions of a first-stage model, the constraint conditions of a second-stage inner-layer model and the constraint conditions of a second-stage outer-layer model; the constraint conditions of the first section model comprise VPP internal power balance constraint, wind power constraint, photovoltaic constraint, distributed energy storage constraint and adjustable load constraint; the constraint conditions of the second-stage inner layer model comprise VPP competitive bidding electric quantity constraint, VPP up-down capacity constraint and VPP up-down mileage constraint; and the constraint conditions of the second-stage outer layer model comprise energy market power balance constraint, network safety constraint, system frequency modulation capacity constraint, system frequency modulation mileage constraint and conventional unit power constraint.
The VPP coordination optimization scheduling model is as follows:
wherein the content of the first and second substances,representing the amount of power VPP delivered to the grid at time t,and (4) predicting the load quantity in the power grid in the period t.
The VPP participation day-ahead market profit maximization model is as follows:
where R denotes the scenario of a competitor's offer, probrRepresenting the probability of occurrence of a scenario R, I representing the investment portfolio parameter set of the VPP, D representing the internal load parameter set of the virtual power plant, T representing the time set of 24h,the system marginal liquidity price representing the energy market at day-ahead in scenario R and time t,representing the amount of electricity reported by the ith unit in the market trading day before and accepted by the ISO in the scenario R and time t,representing the load declaration curve that the d-th load is reported in the market before the day and is accepted by ISO in the scene R and the time t,represents the clearing price of the day-ahead service market capacity in scenario R and time t,representing the up-modulation capacity reported by the ith unit in the service market at the day before and accepted by ISO in the scenario R and time t,represents the lower frequency modulation capacity reported by the ith unit in the service market at the day before and accepted by ISO in the scene and time t,represents the clearing price of the service market fm mileage at scene R and time t,representing the frequency-up mileage reported by the ith unit in the service market in the day-ahead and accepted by ISO in the scene R and the time t,and (3) indicating the frequency-down mileage which is reported by the ith unit in the service market at the day before and accepted by ISO in the scene R and the time t.
The social welfare maximization model is as follows:
wherein the content of the first and second substances,represents the quotation of the discharge of a plurality of small micro-bodies j in a virtual power plant in the scene R and the time t in the energy market,at scene R andsimulating the price of charging the multielement micro-body j in the power plant in the energy market in the interval t,representing the quoted frequency modulation capacity of a plurality of small micro-agents j in a virtual power plant in the scene R and the time t in the auxiliary service market,represents the quotation of the multiple small micro-main bodies in the virtual power plant for the day-ahead frequency modulation market mileage in the scene R and the time t,representing the quoted capacity of the conventional unit g in the auxiliary service market for the scenario R and time t,representing the offer of mileage in the fm market for the conventional unit g in the scenario R and time t,representing the amount of electricity reported by the multivariate micro agent j in the market trading day before and accepted by the ISO in the scenario R and time t,representing the load declaration curve that the load is reported in the market for the day-ahead trading and accepted by the ISO in the scenario R and the time t,the up-modulation capacity reported by the multivariate micro-agent j in the service market at the day before and accepted by ISO in the scene R and time t,the down-tuning capacity reported by the multivariate micro-agent j in the service market in the day ahead and accepted by ISO in the scenario R and time t,representing the up-frequency mileage reported by the multivariate micro-agent j in the service market in the day before and accepted by ISO in the scene R and the time t,represents the down-frequency mileage reported by the multivariate micro-agent j in the service market in the day before and accepted by ISO in the scene R and the time t,represents the lower frequency-modulation capacity reported by the conventional unit g in the service market in the day ahead and accepted by the ISO in the scene and time t,represents the lower frequency-modulation capacity reported by the conventional unit g in the service market in the day ahead and accepted by the ISO in the scene and time t,representing the up-tuned mileage reported by the conventional unit g in the day-ahead service market and accepted by ISO within the scenario R and time t,represents the lower frequency-modulation mileage reported by the conventional unit g in the service market in the day ahead and accepted by the ISO in the scene R and the time t.
The VPP internal power balance constraint is:
wherein the content of the first and second substances,the output of the wind generating set in the virtual power plant in the time t is shown,photovoltaic power generation in virtual power plant within representation time tThe magnitude of the electric output force is large,representing the amount of stored energy in the virtual power plant during time t,represents the amount of stored energy and discharge in a virtual power plant within the time t, PL,shRepresenting the maximum load transferable during time t, LtRepresenting the amount of load in the grid over time t,representing the electric quantity of wind power, photovoltaic and stored energy to the VPP internal load within the time t;
the wind power constraint is as follows:
wherein the content of the first and second substances,representing the electric quantity delivered to the power grid by the wind generating set at the time of the VPP at t,representing the amount of power supplied by the wind turbines to the loads in the virtual power plant during time t,representing the amount of power supplied by the wind turbine to the energy storage unit over time t,representing the maximum power generation of the wind turbine generator set within time t;
the photovoltaic constraints are:
wherein the content of the first and second substances,the amount of power delivered to the grid by the photovoltaic power generation in the VPP at time t,in time t, the electric quantity supplied to the load in the virtual power plant by the photovoltaic power generation,representing the amount of power supplied by the photovoltaic power generation to the energy storage unit during time t,representing the maximum power generation of the wind turbine generator set over time t;
the distributed energy storage constraints are:
wherein, PBSS,dRepresenting the amount of stored energy and discharge in the virtual power plant at time t, PBSS,cRepresenting the amount of stored energy in the virtual power plant at time t,representing the size of the energy storage capacity in the virtual power plant at the moment t,the energy storage capacity in the virtual power plant at the time t-1 is shown, delta t represents the time interval of charging and discharging,are all binary variables, and are provided with the following functions,indicating that the energy storage unit is in a charging state,indicating that the energy storage unit is not in a charging state,indicating that the energy storage unit is in a discharge state,indicating that the energy storage unit is not in the discharge state, etaBSSIndicating the charging efficiency of the energy storage unit, EminIndicating minimum storage capacity of energy storage plant, EmaxRepresenting the maximum storage capacity of the energy storage power station;
the adjustable load constraints are:
wherein the content of the first and second substances,representing the power of the internal load of the virtual power plant at time t,representing the amount of power purchased by the load in the VPP to the grid at time t,representing the amount of load reduction during the time t,which indicates the amount of increase of the load during the time t,are all binary variables, and are provided with the following functions,representing a reduction in the virtual plant load over time t,representing a virtual plant load increase over time t, otherwisePL,sh,totRepresenting the maximum load allowed to be transferred during the entire scheduling.
The VPP bidding electric quantity constraint is as follows:
wherein the content of the first and second substances,for the amount of power sold by the VPP during the day-ahead energy market time t,the method comprises the following steps that the electricity purchasing quantity of a VPP in the day-ahead energy market time t is shown, and G is the maximum competitive bidding electricity quantity of an electricity market;
the VPP drop-down capacity constraint is:
wherein the content of the first and second substances,representing the upper frequency modulation capacity reported by the ith unit in the service market at the moment in the scene R and the time t,the lower frequency modulation capacity reported by the ith unit in the service market at the moment in the scene R and the time t is shown,representing the electricity sales amount reported by the ith unit in the energy market in the scene R and the time t,the power purchasing quantity reported by the ith unit in the energy market in the scene R and the time t,represents the maximum capacity of the ith unit in the VPP,represents the maximum turndown capacity of the ith unit in the VPP,indicating the maximum charging power of the VPP market at the day-ahead,represents the maximum discharge power of the VPP on the market at the day-ahead;
the VPP up-down mileage constraint is as follows:
wherein the content of the first and second substances,showing the frequency-up mileage reported by the ith unit in the service market in the time t and the scene R,showing the lower frequency-modulation mileage, epsilon, reported by the ith unit in the service market in the time t and in the scene RgA mileage multiplier representing a multivariate micro-body g.
The energy market power balance constraint is:
in the formula: pL,tExpressed as the magnitude of the system load over time t;
the network security constraints are:
pl,min≤pl,t≤pl,max
wherein p isl,tFor the flow of the line in time t, pl,maxFor maximum transmission capacity of the line, pl,minRespectively, the minimum transmission capacity of the line;
the system frequency modulation capacity constraint is as follows:
wherein the content of the first and second substances,representing the up-modulation capacity of the system over time t,representing the down-modulation capacity of the system over time t,indicating the up-conversion capacity in VPP during time t,indicating the down-scaled capacity in VPP over time t,represents the upper frequency modulation capacity of the conventional unit in time t,expressing the lower frequency-reducing capacity of the bid in the conventional unit within time t;
the system frequency modulation mileage constraint is as follows:
wherein the content of the first and second substances,represents the up-modulated mileage of the system over time t,represents the lower frequency-modulated mileage of the system over time t,represents the up-frequency-modulated mileage in the VPP within the time t,represents the lower frequency-reduction mileage in VPP within time t,represents the bidding mileage of the conventional unit within the time t,expressing the lower frequency-regulating mileage of the conventional unit within time t;
the conventional unit power constraint is as follows:
wherein, Pg,t,rRepresenting the output value, P, of the conventional unit g over time tg,maxIs the maximum output of the conventional unit g, Pg,minThe minimum output of the conventional unit g.
The method for obtaining the optimal bidding result and the optimal profit of the day-ahead energy and duration day-ahead frequency modulation market comprises the following steps:
step 3.1, carrying out data preprocessing on the VPP after the multiple small micro-main bodies are gathered, and generating multiple bidding strategies k by taking the maximum profit of the VPP as an objective function;
step 3.2, generating a typical quotation scene r of the competitor according to the marginal cost of the competitor multiplied by a curve random scale factor [0.75,1.25 ];
step 3.3, the energy market and the frequency modulation market are cleared in a combined mode in the day ahead;
and 3.4, calculating the probability weighted average of all power price prediction scenes to obtain clearing prices in each market time t in the day, and calculating the total profit of the VPP under each strategy.
The bid-winning condition of each main body in the market comprises a bid-winning condition of each main body in the day-ahead energy market of the market, a bid-winning condition of the day-ahead auxiliary market frequency modulation capacity of each main body in the market and a bid-winning condition of day-ahead auxiliary market frequency modulation mileage of each main body in the market; the optimal quotation comprises the optimal quotation of the VPP in each market and the optimal report of the VPP in each market; the clearing price refers to the clearing price of a day-ahead energy market and a day-ahead frequency modulation market.
The beneficial effect that this technical scheme can produce:
1. the invention can establish a power price bidding model fully considering the uncertainty of the price of electricity in the day-ahead market and day-ahead frequency modulation, and provides a reference basis for VPP to carry out market bidding, operation decision and energy storage economic feasibility assessment.
2. The method can construct a clearing model of the day-ahead energy and day-ahead frequency modulation market, and enables the VPP on the power grid side to be used as a buyer and seller to participate in the unified clearing process of the day-ahead energy and day-ahead frequency modulation market. Therefore, the VPP can be effectively stimulated to participate in frequency modulation auxiliary service, more benefits are obtained, the loss caused by the output power of the thermal power generating unit is reduced frequently, the electricity purchasing cost of a demand party is reduced, and the benefit requirements of supply and demand parties are met.
3. The method can realize the related algorithm through programming, obtain the output result by reading the input data, and visually display the output result in a data and graph mode by carrying out visual interface design on each module interface, thereby more simply and conveniently realizing the participation of the VPP in the frequency modulation auxiliary service market transaction process.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a diagram of a game model architecture for VPP participation in spot market bidding in accordance with the present invention;
FIG. 3 is a flow chart of model solution of the present invention;
FIG. 4 is a schematic diagram of visualized charge and discharge power of the energy storage power station in each time period according to the present invention;
FIG. 5 is a diagram illustrating the power consumption after VPP internal load optimization according to the present invention;
FIG. 6 is a diagram of system load power prediction in accordance with the present invention;
FIG. 7 is a chart of bid winning status in the day-ahead energy market of each subject of the market according to the present invention;
FIG. 8 is a plot of the bid amount of the FM capacity of the day-ahead auxiliary market for each market subject of the present invention;
FIG. 9 is a chart of bid-winning status of the frequency-modulated mileage of each market subject in the day-ahead auxiliary market of the present invention;
FIG. 10 is a diagram of the best bid in each market for the VPP of the present invention;
FIG. 11 is a graph of the optimal volume of the VPP of the present invention in various markets;
FIG. 12 is a plot of the day-ahead energy market consolidated clearing price profile and total load curve of the present invention;
FIG. 13 is a plot of the day-ahead auxiliary market FM capacity distribution and FM mileage clearing price curves of the present invention.
Detailed Description
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 only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be obtained by a person skilled in the art without inventive effort based on the embodiments of the present invention, are within the scope of the present invention.
As shown in fig. 1, an embodiment of the present invention provides a virtual power plant-based method for clearing a plurality of small micro-subjects participating in spot market, which includes the following specific steps:
the method comprises the following steps: constructing a double-layer clearing model considering the uncertainty of the electricity price based on the VPP parameters of the power grid side and the information of the electric energy buyer and the electric energy seller;
in this embodiment, a day-ahead energy and day-ahead auxiliary service market trading platform may be provided, and based on an energy and frequency modulation market trading mode of day-ahead centralized bidding trading provided in the trading platform, the electric energy buyer and seller may claim related information, so that the trading platform may form a day-ahead energy and frequency modulation market for centralized bidding and unified clearing of electricity prices.
Generally, the information of the buyer and the seller of the electric energy includes price information and electric quantity information of each buyer and the seller, such as parameters of different quotes of the conventional generator set and the VPP in different bidding periods, the price of the outgoing clear electricity in the day, the profit of the outgoing clear electricity in the day, the electric quantity sold and the electric quantity purchased.
It should be noted that, the conventional generator set is an electric energy seller, and the VPP may be used as an electric energy seller in the market or an electric energy buyer in the market. Which can be profitable by low-buying high-selling through the two-stage double-deck pull-out model described herein according to time-of-use pricing.
Meanwhile, the VPP parameter on the grid side generally includes a parameter related to the energy storage unit and a parameter related to the renewable energy unit. For example, the bid amount of the energy storage unit and the renewable energy source, the number of the energy storage unit participating in frequency modulation and the demand side response user, the cyclic loss cost of one day of energy storage, the number of charging and discharging cycles of the energy storage full life cycle, the energy storage charging and discharging depth and the like.
Generally, in a mature market, the price quoted for electricity generation should be able to represent the cost of electricity generation. Therefore, in the present embodiment, the compensation fee for the settlement of the quoted price can be approximately regarded as the power generation cost of the power generation manufacturer. If the power generation manufacturer quotes based on the power generation cost or the energy storage cost, a clearing model can be constructed based on the uncertainty of the electricity price according to the VPP parameters of the power grid side and the information of the electric energy buyers and sellers.
As shown in fig. 2, the double-layer export model includes a first-stage model, a second-stage inner-layer model and a second-stage outer-layer model, wherein the first-stage model is a VPP coordinated optimization scheduling model, the second-stage inner-layer model is a VPP participation day-ahead market profit maximization model, and the second-stage outer-layer model is a social welfare maximization model.
The VPP coordinated optimization scheduling model coordinates and schedules each small micro main body in the main body by taking the maximum target of the goodness of fit between the output characteristic curve of each main body and the load curve of the power grid. And taking the optimized output curve as the basis of the next-stage system operation optimization, wherein the target function is as follows:
wherein the content of the first and second substances,representing the amount of power VPP delivered to the grid at time t,and (4) predicting the load quantity in the power grid in the period t.
The VPP participation day-ahead market profit maximization model is as follows:
where R denotes the scenario of a competitor's offer, probrRepresenting the probability of occurrence of a scenario R, I representing the investment portfolio parameter set of the VPP, D representing the internal load parameter set of the virtual power plant, T representing the time set of 24h,the system marginal liquidity price representing the energy market at day-ahead in scenario R and time t,representing the amount of electricity reported by the ith unit in the market trading day before and accepted by the ISO in the scenario R and time t,representing the load declaration curve that the d-th load is reported in the market before the day and is accepted by ISO in the scene R and the time t,represents the clearing price of the day-ahead service market capacity in scenario R and time t,representing the up-modulation capacity reported by the ith unit in the service market at the day before and accepted by ISO in the scenario R and time t,represents the lower frequency modulation capacity reported by the ith unit in the service market at the day before and accepted by ISO in the scene and time t,represents the clearing price of the service market fm mileage at scene R and time t,representing the frequency-up mileage reported by the ith unit in the service market in the day-ahead and accepted by ISO in the scene R and the time t,and (3) indicating the frequency-down mileage which is reported by the ith unit in the service market at the day before and accepted by ISO in the scene R and the time t.
The social welfare maximization model is as follows:
wherein the content of the first and second substances,represents the quotation of the discharge of a plurality of small micro-bodies j in a virtual power plant in the scene R and the time t in the energy market,simulating the price of charging the multi-element micro main body j in the power plant in the scene R and the time t in the energy market,representing the quoted frequency modulation capacity of a plurality of small micro-agents j in a virtual power plant in the scene R and the time t in the auxiliary service market,represents the quotation of the multiple small micro-main bodies in the virtual power plant for the day-ahead frequency modulation market mileage in the scene R and the time t,representing the quoted capacity of the conventional unit g in the auxiliary service market for the scenario R and time t,representing the offer of mileage in the fm market for the conventional unit g in the scenario R and time t,representing the amount of electricity reported by the multivariate micro agent j in the market trading day before and accepted by the ISO in the scenario R and time t,representing the load declaration curve that the load is reported in the market for the day-ahead trading and accepted by the ISO in the scenario R and the time t,at scene R and timethe up-modulation capacity reported by the multivariate micro-agent j in the service market at the day before and accepted by ISO within t,the down-tuning capacity reported by the multivariate micro-agent j in the service market in the day ahead and accepted by ISO in the scenario R and time t,representing the up-frequency mileage reported by the multivariate micro-agent j in the service market in the day before and accepted by ISO in the scene R and the time t,represents the down-frequency mileage reported by the multivariate micro-agent j in the service market in the day before and accepted by ISO in the scene R and the time t,represents the lower frequency-modulation capacity reported by the conventional unit g in the service market in the day ahead and accepted by the ISO in the scene and time t,represents the lower frequency-modulation capacity reported by the conventional unit g in the service market in the day ahead and accepted by the ISO in the scene and time t,representing the up-tuned mileage reported by the conventional unit g in the day-ahead service market and accepted by ISO within the scenario R and time t,represents the lower frequency-modulation mileage reported by the conventional unit g in the service market in the day ahead and accepted by the ISO in the scene R and the time t.
Step two: constructing a constraint condition of a double-layer clearing model based on the VPP parameters of the power grid side, the information of both electric energy buyers and sellers and the double-layer clearing model;
the constraint conditions of the double-layer output model comprise the constraint conditions of a first-stage model, the constraint conditions of a second-stage inner-layer model and the constraint conditions of a second-stage outer-layer model; the constraint conditions of the first section model comprise VPP internal power balance constraint, wind power constraint, photovoltaic constraint, distributed energy storage constraint and adjustable load constraint; the constraint conditions of the second-stage inner layer model comprise VPP competitive bidding electric quantity constraint, VPP up-down capacity constraint and VPP up-down mileage constraint; and the constraint conditions of the second-stage outer layer model comprise energy market power balance constraint, network safety constraint, system frequency modulation capacity constraint, system frequency modulation mileage constraint and conventional unit power constraint.
The VPP internal power balance constraint comprises absorption quantities of a plurality of small micro bodies with output equal to load and stored energy, and the expression is as follows:
wherein the content of the first and second substances,the output of the wind generating set in the virtual power plant in the time t is shown,represents the magnitude of the photovoltaic power generation output in the virtual power plant within the time t,representing the amount of stored energy in the virtual power plant during time t,represents the amount of stored energy and discharge in a virtual power plant within the time t, PL,shRepresenting the maximum load transferable during time t, LtRepresenting the amount of load in the grid over time t,electric quantity for expressing wind power, photovoltaic and energy storage to supply internal load of VPP within time t。
The wind power constraint is that the output power of the wind power is equal to the sum of the internal load given by the wind power and the internal load of the virtual power plant and is constrained by the maximum power, and the expression is as follows:
wherein the content of the first and second substances,representing the electric quantity delivered to the power grid by the wind generating set at the time of the VPP at t,representing the amount of power supplied by the wind turbines to the loads in the virtual power plant during time t,representing the amount of power supplied by the wind turbine to the energy storage unit over time t,representing the maximum power production of the wind park over time t.
The photovoltaic constraint is that the output power of the photovoltaic is equal to the sum of the internal load given by the photovoltaic and the internal load of the virtual power plant and is constrained by the maximum power, and the expression is as follows:
wherein the content of the first and second substances,the amount of power delivered to the grid by the photovoltaic power generation in the VPP at time t,in time t, the electric quantity supplied to the load in the virtual power plant by the photovoltaic power generation,representing the amount of power supplied by the photovoltaic power generation to the energy storage unit during time t,representing the maximum power generation of the wind park over time t.
The distributed energy storage constraint comprises charge-discharge constraint, capacity constraint and state constraint of all energy storage power stations, and the expression is as follows:
wherein, PBSS,dRepresenting the amount of stored energy and discharge in the virtual power plant at time t, PBSS,cRepresenting the amount of stored energy in the virtual power plant at time t,representing the size of the energy storage capacity in the virtual power plant at the moment t,the energy storage capacity in the virtual power plant at the time t-1 is shown, delta t represents the time interval of charging and discharging,are all binary variables, and are provided with the following functions,indicating that the energy storage unit is in a charging state,indicating that the energy storage unit is not in a charging state,indicating that the energy storage unit is in a discharge state,indicating that the energy storage unit is not in the discharge state, etaBSSIndicating the charging efficiency of the energy storage unit, EminIndicating minimum storage capacity of energy storage plant, EmaxRepresenting the maximum storage capacity of the energy storage power station;
the adjustable load constraints comprise all adjustable load up-down capacity constraints and state constraints, and the expression is as follows:
wherein the content of the first and second substances,representing the power of the internal load of the virtual power plant at time t,representing the amount of power purchased by the load in the VPP to the grid at time t,representing the amount of load reduction during the time t,which indicates the amount of increase of the load during the time t,are all binary variables, and are provided with the following functions,representing a reduction in the virtual plant load over time t,representing a virtual plant load increase over time t, otherwisePL,sh,totThe VP representing the maximum load allowed to be transferred during the entire schedulingThe P bidding electric quantity constraint comprises VPP and main network transmission power limitation, and the bidding electric quantity of the VPP in the electric power market needs to meet the following constraint formula:
wherein the content of the first and second substances,for the amount of power sold by the VPP during the day-ahead energy market time t,the electricity quantity purchased by the VPP in the energy market time t in the day ahead, and G is the maximum competitive bidding electricity quantity of the electricity market.
The VPP up-down capacity constraint comprises a maximum up-down capacity constraint and a maximum charge-discharge constraint, and the expression is as follows:
wherein the content of the first and second substances,representing the upper frequency modulation capacity reported by the ith unit in the service market at the moment in the scene R and the time t,the lower frequency modulation capacity reported by the ith unit in the service market at the moment in the scene R and the time t is shown,representing the electricity sales amount reported by the ith unit in the energy market in the scene R and the time t,the power purchasing quantity reported by the ith unit in the energy market in the scene R and the time t,represents the maximum capacity of the ith unit in the VPP,represents the maximum turndown capacity of the ith unit in the VPP,indicating the maximum charging power of the VPP market at the day-ahead,indicating the maximum discharge power of the VPP market at the day-ahead.
The VPP up-down mileage constraint comprises a maximum up-down history constraint, and the expression is as follows:
wherein the content of the first and second substances,showing the frequency-up mileage reported by the ith unit in the service market in the time t and the scene R,showing the lower frequency-modulation mileage, epsilon, reported by the ith unit in the service market in the time t and in the scene RgA mileage multiplier representing a multivariate micro-body g.
The energy market power balance constraint comprises that the discharge power of the virtual power plant and the discharge power of the conventional unit are equal to the total load required by the system, and the expression is as follows:
in the formula: pL,tExpressed as the magnitude of the system load over time t.
The network security constraint comprises maximum and minimum transmission capacity constraints of a line, and the expression is as follows:
pl,min≤pl,t≤pl,max
wherein p isl,tFor the flow of the line in time t, pl,maxFor maximum transmission capacity of the line, pl,minRespectively the minimum transmission capacity of the line.
The system frequency modulation capacity constraint comprises that the total up-down regulation capacity of the system is equal to the sum of the up-down regulation capacity of the virtual power plant and the up-down regulation capacity of the unit, and the expression is as follows:
wherein the content of the first and second substances,representing the up-modulation capacity of the system over time t,representing the down-modulation capacity of the system over time t,indicating the up-conversion capacity in VPP during time t,indicating the down-scaled capacity in VPP over time t,represents the upper frequency modulation capacity of the conventional unit in time t,and (4) indicating the lower frequency-reducing capacity of the conventional unit during the time t.
The system frequency modulation mileage constraint comprises that the total upper and lower mileage of the system is equal to the sum of the upper and lower mileage of the virtual power plant and the upper and lower mileage of the unit, and the expression is as follows:
wherein the content of the first and second substances,represents the up-modulated mileage of the system over time t,represents the lower frequency-modulated mileage of the system over time t,represents the up-frequency-modulated mileage in the VPP within the time t,represents the lower frequency-reduction mileage in VPP within time t,represents the bidding mileage of the conventional unit within the time t,and (4) indicating the lower frequency-regulating mileage in the conventional unit within the time t.
The conventional unit power constraint comprises maximum up-down regulation power constraint, and the expression is as follows:
wherein, Pg,t,rRepresenting the output value, P, of the conventional unit g over time tg,maxIs the maximum output of the conventional unit g, Pg,minThe minimum output of the conventional unit g. The conventional unit constraint also comprises other constraints such as climbing constraint, declared frequency modulation capacity and mileage constraint and the like. Wherein, the declared frequency modulation capacity and mileage constraint are the same as the energy storage form. The restriction of the bid amount includes that the sum of the bid amounts in each time interval of the conventional unit does not exceed the sum of the bid amounts in each time interval of the conventional unitThe sum of the rated capacity and the bid amount in each energy storage period does not exceed the rated capacity and the rated power.
Step three: solving the double-layer clearing model based on the constraint conditions of the double-layer clearing model to obtain the optimal bidding result and the optimal profit of the day-ahead frequency modulation market of the day-ahead energy and duration;
as shown in fig. 3, the method for obtaining the optimal bidding result and the optimal profit of the day-ahead energy and duration day-ahead frequency modulation market includes:
step 3.1, carrying out data preprocessing on the VPP after the multiple small micro-main bodies are gathered, and generating multiple bidding strategies k by taking the maximum profit of the VPP as an objective function; in step 3.1, due to uncertainty and randomness of output of renewable energy sources (wind power, photovoltaic and the like) and distribution universality, scheduling of the renewable energy sources is difficult. With the development of the power internet of things and the data communication technology, the multiple small micro-bodies can be aggregated into one VPP and scheduled by a unified operator. After a VPP is aggregated, an operator needs to perform data preprocessing on the day-ahead energy and frequency modulation market quotation information input by each main body based on methods such as abnormal value screening and supplementing, feature selection, optimized parameter configuration and the like to generate a quotation and report curve for 24h in the future, so that various bidding strategies of the day-ahead energy and frequency modulation market are regenerated.
Step 3.2, generating a typical quotation scene r of the competitor according to the marginal cost of the competitor multiplied by a curve random scale factor [0.75,1.25 ]; in step 3.2, since the multi-step bidding strategy is generated in step 3.1, but the bidding strategies of other competitors, such as the bidding strategies of the traditional unit and the bidding strategies of other VPPs, are also considered, the ARMA is adopted to fit the data, a series of data is obtained through Monte Carlo sampling, and the scene is cut down by adopting the rapid antecedent method, so that a typical target scene set of the competitors is obtained.
Step 3.3, the energy market and the frequency modulation market are cleared in a combined mode in the day ahead; in step 3.3, the VPP is used as a leader in a competitive game to ensure that the energy of the VPP and the income of the frequency modulation market are maximized in the day ahead, the ISO is used as a follower in the competitive game to listen to the quoted price and the quoted volume of each unit, the maximization of the social welfare is taken as a target, namely the minimum electricity purchasing cost is taken as a target, the combined clearing of the energy market and the frequency modulation market in the day ahead is carried out, and the competitive electric quantity and the competitive price of each unit for 24 hours are obtained.
And 3.4, calculating the probability weighted average of all power price prediction scenes to obtain clearing prices in each market time t in the day, and calculating the total profit of the VPP under each strategy. In step 3.4, considering the uncertainty of the scenes, each scene corresponds to a probability, the successful bidding electric quantity and the bidding price of each main body in each scene are obtained, and the clearing price of each market in the day-ahead time t period in each scene is obtained through the weighted average. And obtaining the best quotation strategy of the VPP participating in the day-ahead energy market and the day-ahead frequency modulation market by comparing the total profit of the VPP in each scene.
Through the steps 3.1 to 3.4, the optimal bidding result and the optimal income of the seller in the market can be obtained according to the two-stage double-layer clearing model and the two-stage double-layer constraint condition.
Step four: and obtaining the successful bidding condition, the optimal quotation and the clearing price of each main body of the market according to the optimal bidding result and the optimal income. The bid-winning condition of each main body in the market comprises a bid-winning condition of each main body in the day-ahead energy market of the market, a bid-winning condition of the day-ahead auxiliary market frequency modulation capacity of each main body in the market and a bid-winning condition of day-ahead auxiliary market frequency modulation mileage of each main body in the market; the optimal quotation comprises the optimal quotation of the VPP in each market and the optimal report of the VPP in each market; the clearing price refers to the clearing price of a day-ahead energy market and a day-ahead frequency modulation market. The embodiment outputs a bidding situation diagram of each main body day-ahead energy market, a bidding capacity bidding situation diagram of each main body day-ahead auxiliary market, a bidding mileage bidding situation diagram of each main body day-ahead auxiliary market, a different market main body profit comparison table, an optimal quotation diagram of VPP in each market, and a clearing price diagram of the day-ahead energy market and the day-ahead frequency modulation market based on the regional dimension and the time dimension in a visual mode.
By encapsulating the related algorithms involved in the steps one to three, the process of clearing can be displayed through a simple and convenient visual interface. For example, a Graphical User Interface (GUI) module in MATLAB may be employed to package the relevant algorithms.
In the embodiment, an example is constructed by a VPP formed by a wind turbine generator, a photovoltaic generator, an energy storage generator and 4000 families, and the output parameters of each small micro main body in the VPP are shown in figures 4 and 5; the total load situation of the system is shown in figure 6.
As can be seen from fig. 7, in the energy market at day time, the conventional units G1, G2 provide the vast majority of energy, whereas VPP is mainly 1: the energy is provided during the periods of 00-6: 00 and 13: 00-15: 00, because the wind power provides larger output at night, the sunlight has stronger illumination capability at noon, and the photovoltaic provides larger output at the moment. In the following step 9: 00-11: 00 and 17: 00-19: during 00 the VPP needs to purchase power from the grid due to the large power usage by the load within the VPP (see fig. 4).
Comparing fig. 8 and 9, the VPP assumes the main fm capacity and fm mileage tasks of the system. This is because when the same modulation capacity is provided as compared with the modulation output of the conventional unit, the VPP unit, which includes the distributed energy storage device, responds to the modulation signal more frequently in the upward and downward directions, and thus can provide more frequency modulation mileage than the conventional unit, and is preferentially called by the system in the frequency modulation market. The revenue situation of different market entities in each market is shown in table 1.
TABLE 1 comparison table of different market subject profits
As can be seen from table 1, the main gains of the conventional units G1, G2 and G6 mainly come from the day-ahead energy market, and the main gains of the VPP mainly include the day-ahead energy market gains and the frequency modulation capacity gains.
As can be seen from FIG. 10, the VPP discharge price is always higher than the charge price because the VPP earns profit from the price difference. The VPP FM capacity quote, although fluctuating, remains entirely on a straight line, whereas the VPP FM mileage quote is a straight line.
As can be seen more intuitively from fig. 11, the VPP charges its internal energy storage during periods of low load, and discharges it when the load is high. The VPP frequency modulation capability is higher, so the frequency modulation capacity and the medium-rate rating of the frequency modulation mileage are higher.
As can be seen from fig. 12, the clear electricity price of the energy market is similar to the load trend, and shows a variation trend of low night and high day, because the conventional units quote according to the power generation cost, the quoted price is linearly related to the power generation amount, and the quoted price trend is similar to the load trend.
As can be seen from fig. 13, the fm capacity and the odds price of fm mileage of the aftermarket fluctuate over time. They all follow a law: 4: 00-6: price of 00 out is lowest, 12: 00-14: 00. 20: 00-23: the 00 clearing price is the highest, but the overall clearing price is not greatly fluctuated, and is maintained at an approximate level.
According to the invention, by establishing a clearing model and constraint conditions, visual benefit analysis and bidding decision analysis can be carried out on the current energy and frequency modulation market in which VPP participates on the basis of considering the uncertainty of the electricity price. The invention can establish the electric power price bidding model fully considering the uncertainty of the price of electricity in the market in the day ahead, and provides a reference basis for market bidding, operation decision and energy storage economic feasibility evaluation of VPPs. The method can construct a clearing model of the day-ahead energy and frequency modulation market, and enables the VPP on the power grid side to be used as a buyer and a seller to participate in the unified clearing process of the day-ahead energy and frequency modulation market. Therefore, the VPP can be effectively stimulated to participate in the peak shaving auxiliary service, more benefits are obtained, the loss caused by the output power of the thermal power generating unit is reduced frequently, the electricity purchasing cost of a demand party is reduced, and the benefit requirements of supply and demand parties are met. The invention can realize the related algorithm through programming, obtain the output result through reading the input data, carry on the visual interface design to each module interface so as to show the output result through the way of data and figure visualization, realize VPP to the participation of the energy and frequency modulation market transaction process before the day more simply, conveniently.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents, improvements and the like that fall within the spirit and principle of the present invention are intended to be included therein.
Claims (10)
1. A virtual power plant-based method for clearing a plurality of small micro-subjects participating in spot market is characterized by comprising the following steps:
the method comprises the following steps: constructing a double-layer clearing model considering the uncertainty of the electricity price based on the VPP parameters of the power grid side and the information of the electric energy buyer and the electric energy seller;
step two: constructing a constraint condition of a double-layer clearing model based on the VPP parameters of the power grid side, the information of both electric energy buyers and sellers and the double-layer clearing model;
step three: solving the double-layer clearing model based on the constraint conditions of the double-layer clearing model to obtain the optimal bidding result and the optimal profit of the day-ahead frequency modulation market of the day-ahead energy and duration;
step four: and obtaining the successful bidding condition, the optimal quotation and the clearing price of each main body of the market according to the optimal bidding result and the optimal income.
2. The virtual power plant-based multi-element micro-subject spot market clearing method according to claim 1, wherein the double-layer clearing model comprises a first-stage model, a second-stage inner-layer model and a second-stage outer-layer model, wherein the first-stage model is a VPP coordination optimization scheduling model, the second-stage inner-layer model is a VPP participation day-ahead market profit maximization model, and the second-stage outer-layer model is a full social welfare maximization model; the constraint conditions of the double-layer output model comprise the constraint conditions of a first-stage model, the constraint conditions of a second-stage inner-layer model and the constraint conditions of a second-stage outer-layer model; the constraint conditions of the first section model comprise VPP internal power balance constraint, wind power constraint, photovoltaic constraint, distributed energy storage constraint and adjustable load constraint; the constraint conditions of the second-stage inner layer model comprise VPP competitive bidding electric quantity constraint, VPP up-down capacity constraint and VPP up-down mileage constraint; and the constraint conditions of the second-stage outer layer model comprise energy market power balance constraint, network safety constraint, system frequency modulation capacity constraint, system frequency modulation mileage constraint and conventional unit power constraint.
3. The virtual power plant based multi-element micro-agent spot market clearing method according to claim 2, wherein the VPP coordinated optimization scheduling model is:
4. The virtual power plant based export method of multivariate micro-agents participating in spot market of claim 2, wherein the VPP participation day-ahead market profit maximization model is:
where R denotes the scenario of a competitor's offer, probrRepresenting the probability of occurrence of a scenario R, I representing the investment portfolio parameter set of the VPP, D representing the internal load parameter set of the virtual power plant, T representing the time set of 24h,the system marginal liquidity price representing the energy market at day-ahead in scenario R and time t,representing the amount of electricity reported by the ith unit in the market trading day before and accepted by the ISO in the scenario R and time t,representing the load declaration curve that the d-th load is reported in the market before the day and is accepted by ISO in the scene R and the time t,represents the clearing price of the day-ahead service market capacity in scenario R and time t,representing the up-modulation capacity reported by the ith unit in the service market at the day before and accepted by ISO in the scenario R and time t,represents the lower frequency modulation capacity reported by the ith unit in the service market at the day before and accepted by ISO in the scene and time t,represents the clearing price of the service market fm mileage at scene R and time t,representing the frequency-up mileage reported by the ith unit in the service market in the day-ahead and accepted by ISO in the scene R and the time t,and (3) indicating the frequency-down mileage which is reported by the ith unit in the service market at the day before and accepted by ISO in the scene R and the time t.
5. The virtual power plant based multi-element micro-agent spot market clearing method according to claim 4, wherein the social welfare maximization model is:
wherein the content of the first and second substances,represents the quotation of the discharge of a plurality of small micro-bodies j in a virtual power plant in the scene R and the time t in the energy market,simulating the price of charging the multi-element micro main body j in the power plant in the scene R and the time t in the energy market,representing the quoted frequency modulation capacity of a plurality of small micro-agents j in a virtual power plant in the scene R and the time t in the auxiliary service market,represents the quotation of the multiple small micro-main bodies in the virtual power plant for the day-ahead frequency modulation market mileage in the scene R and the time t,representing the quoted capacity of the conventional unit g in the auxiliary service market for the scenario R and time t,representing the offer of mileage in the fm market for the conventional unit g in the scenario R and time t,showing that multiple micro-subjects j are reported in the market of the day-ahead trade in the scene R and the time tAnd the amount of power accepted by the ISO,representing the load declaration curve that the load is reported in the market for the day-ahead trading and accepted by the ISO in the scenario R and the time t,the up-modulation capacity reported by the multivariate micro-agent j in the service market at the day before and accepted by ISO in the scene R and time t,the down-tuning capacity reported by the multivariate micro-agent j in the service market in the day ahead and accepted by ISO in the scenario R and time t,representing the up-frequency mileage reported by the multivariate micro-agent j in the service market in the day before and accepted by ISO in the scene R and the time t,represents the down-frequency mileage reported by the multivariate micro-agent j in the service market in the day before and accepted by ISO in the scene R and the time t,represents the lower frequency-modulation capacity reported by the conventional unit g in the service market in the day ahead and accepted by the ISO in the scene and time t,represents the lower frequency-modulation capacity reported by the conventional unit g in the service market in the day ahead and accepted by the ISO in the scene and time t,shows that the conventional unit g reports in the service market at the day before the scene R and the time tThe up-tuned mileage accepted by ISO,represents the lower frequency-modulation mileage reported by the conventional unit g in the service market in the day ahead and accepted by the ISO in the scene R and the time t.
6. The virtual power plant based multivariate small micro-agent spot market clearing method according to claim 2 or 3, wherein the VPP internal power balance constraint is:
wherein the content of the first and second substances,the output of the wind generating set in the virtual power plant in the time t is shown,represents the magnitude of the photovoltaic power generation output in the virtual power plant within the time t,representing the amount of stored energy in the virtual power plant during time t,represents the amount of stored energy and discharge in a virtual power plant within the time t, PL,shRepresenting the maximum load transferable during time t, LtRepresenting the amount of load in the grid over time t,representing the electric quantity of wind power, photovoltaic and stored energy to the VPP internal load within the time t;
the wind power constraint is as follows:
wherein the content of the first and second substances,representing the electric quantity delivered to the power grid by the wind generating set at the time of the VPP at t,representing the amount of power supplied by the wind turbines to the loads in the virtual power plant during time t,representing the amount of power supplied by the wind turbine to the energy storage unit over time t,representing the maximum power generation of the wind turbine generator set within time t;
the photovoltaic constraints are:
wherein the content of the first and second substances,the amount of power delivered to the grid by the photovoltaic generation in the VPP at the time,in time t, the electric quantity supplied to the load in the virtual power plant by the photovoltaic power generation,representing the amount of power supplied by the photovoltaic power generation to the energy storage unit during time t,representing the maximum power generation of the wind turbine generator set over time t;
the distributed energy storage constraints are:
wherein, PBSS,dRepresenting the amount of stored energy and discharge in the virtual power plant at time t, PBSS,cRepresenting the amount of stored energy in the virtual power plant at time t,representing the size of the energy storage capacity in the virtual power plant at the moment t,the energy storage capacity in the virtual power plant at the time t-1 is shown, delta t represents the time interval of charging and discharging,are all binary variables, and are provided with the following functions,indicating that the energy storage unit is in a charging state,indicating that the energy storage unit is not in a charging state,indicating that the energy storage unit is in a discharge state,indicating that the energy storage unit is not in the discharge state, etaBSSIndicating the charging efficiency of the energy storage unit, EminIndicating minimum storage capacity of energy storage plant, EmaxRepresenting the maximum storage capacity of the energy storage power station;
the adjustable load constraints are:
wherein the content of the first and second substances,representing the power of the internal load of the virtual power plant at time t,representing the amount of power purchased by the load in the VPP to the grid at time t,representing the amount of load reduction during the time t,which indicates the amount of increase of the load during the time t,are all binary variables, and are provided with the following functions,representing a reduction in the virtual plant load over time t,representing a virtual plant load increase over time t, otherwisePL,sh,totRepresenting the maximum load allowed to be transferred during the entire scheduling.
7. The virtual power plant-based multi-element micro-agent spot market clearing method according to claim 2 or 4, wherein the VPP bidding power constraint is as follows:
wherein the content of the first and second substances,for the amount of power sold by the VPP during the day-ahead energy market time t,the method comprises the following steps that the electricity purchasing quantity of a VPP in the day-ahead energy market time t is shown, and G is the maximum competitive bidding electricity quantity of an electricity market;
the VPP drop-down capacity constraint is:
wherein the content of the first and second substances,representing the upper frequency modulation capacity reported by the ith unit in the service market at the moment in the scene R and the time t,the lower frequency modulation capacity reported by the ith unit in the service market at the moment in the scene R and the time t is shown,representing the electricity sales amount reported by the ith unit in the energy market in the scene R and the time t,the power purchasing quantity reported by the ith unit in the energy market in the scene R and the time t,represents the maximum capacity of the ith unit in the VPP,represents the maximum turndown capacity of the ith unit in the VPP,indicating the maximum charging power of the VPP market at the day-ahead,represents the maximum discharge power of the VPP on the market at the day-ahead;
the VPP up-down mileage constraint is as follows:
wherein the content of the first and second substances,showing the frequency-up mileage reported by the ith unit in the service market in the time t and the scene R,showing the lower frequency-modulation mileage, epsilon, reported by the ith unit in the service market in the time t and in the scene RgA mileage multiplier representing a multivariate micro-body g.
8. The virtual power plant based multi-element micro-agent spot market clearing method according to claim 7, wherein the energy market power balance constraint is:
in the formula: pL,tExpressed as the magnitude of the system load over time t;
the network security constraints are:
pl,min≤pl,t≤pl,max
wherein p isl,tFor the flow of the line in time t, pl,maxFor maximum transmission capacity of the line, pl,minRespectively, the minimum transmission capacity of the line;
the system frequency modulation capacity constraint is as follows:
wherein the content of the first and second substances,representing the up-modulation capacity of the system over time t,representing the down-modulation capacity of the system over time t,indicating the up-conversion capacity in VPP during time t,indicating the down-scaled capacity in VPP over time t,represents the upper frequency modulation capacity of the conventional unit in time t,expressing the lower frequency-reducing capacity of the bid in the conventional unit within time t;
the system frequency modulation mileage constraint is as follows:
wherein the content of the first and second substances,represents the up-modulated mileage of the system over time t,represents the lower frequency-modulated mileage of the system over time t,represents the up-frequency-modulated mileage in the VPP within the time t,represents the lower frequency-reduction mileage in VPP within time t,represents the bidding mileage of the conventional unit within the time t,expressing the lower frequency-regulating mileage of the conventional unit within time t;
the conventional unit power constraint is as follows:
wherein, Pg,t,rRepresenting the output value, P, of the conventional unit g over time tg,maxIs the maximum output of the conventional unit g, Pg,minThe minimum output of the conventional unit g.
9. The virtual power plant-based ex-clearing method for multivariate micro-subjects participating in spot market according to claim 6, wherein the method for obtaining optimal bidding results and optimal profits of day-ahead energy and duration day-ahead frequency-modulated market is as follows:
step 3.1, carrying out data preprocessing on the VPP after the multiple small micro-main bodies are gathered, and generating multiple bidding strategies k by taking the maximum profit of the VPP as an objective function;
step 3.2, generating a typical quotation scene r of the competitor according to the marginal cost of the competitor multiplied by a curve random scale factor [0.75,1.25 ];
step 3.3, the energy market and the frequency modulation market are cleared in a combined mode in the day ahead;
and 3.4, calculating the probability weighted average of all power price prediction scenes to obtain clearing prices in each market time t in the day, and calculating the total profit of the VPP under each strategy.
10. The virtual power plant-based method for coming out of the spot market by participation of multiple small micro-agents in the market according to claim 1, wherein the bid-winning condition of each agent in the market comprises a day-ahead energy market bid-winning condition of each agent in the market, a day-ahead auxiliary market frequency modulation capacity bid-winning condition of each agent in the market, and a day-ahead auxiliary market frequency modulation mileage bid-winning condition of each agent in the market; the optimal quotation comprises the optimal quotation of the VPP in each market and the optimal report of the VPP in each market; the clearing price refers to the clearing price of a day-ahead energy market and a day-ahead frequency modulation market.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011443256.1A CN112529622A (en) | 2020-12-08 | 2020-12-08 | Virtual power plant-based method for clearing multiple micro-micro main bodies participating in spot market |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011443256.1A CN112529622A (en) | 2020-12-08 | 2020-12-08 | Virtual power plant-based method for clearing multiple micro-micro main bodies participating in spot market |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112529622A true CN112529622A (en) | 2021-03-19 |
Family
ID=75000132
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011443256.1A Pending CN112529622A (en) | 2020-12-08 | 2020-12-08 | Virtual power plant-based method for clearing multiple micro-micro main bodies participating in spot market |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112529622A (en) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112952918A (en) * | 2021-03-23 | 2021-06-11 | 国网山西省电力公司 | Electric power peak regulation and frequency modulation combined clearing method based on marginal benefit |
CN113098040A (en) * | 2021-04-09 | 2021-07-09 | 国网新疆电力有限公司经济技术研究院 | Power grid side energy storage capacity optimal configuration method for obtaining multi-scene benefits |
CN113255957A (en) * | 2021-04-02 | 2021-08-13 | 国网山东省电力公司济南供电公司 | Quantitative optimization analysis method and system for uncertain factors of comprehensive service station |
CN113538066A (en) * | 2021-07-30 | 2021-10-22 | 国网上海市电力公司 | Control method, system, equipment and medium for virtual power plant to participate in power market |
CN113610435A (en) * | 2021-08-24 | 2021-11-05 | 内蒙古电力(集团)有限责任公司电力调度控制分公司 | Peak regulation method and device for power system based on thermal power, self-contained power plant and energy storage |
CN113627762A (en) * | 2021-07-30 | 2021-11-09 | 国网山西省电力公司电力科学研究院 | Virtual power plant peak regulation method based on excitation electricity price |
CN113837444A (en) * | 2021-08-26 | 2021-12-24 | 国网河南省电力公司经济技术研究院 | Trading clearing optimization method for virtual power plant participating in multiple markets |
CN114362227A (en) * | 2022-01-12 | 2022-04-15 | 华北电力大学 | Electric automobile frequency modulation reporting method considering user intention |
CN115204443A (en) * | 2022-05-12 | 2022-10-18 | 电力规划总院有限公司 | Wind, light, fire and storage combined participation market bidding optimal strategy model |
CN116823332A (en) * | 2023-06-29 | 2023-09-29 | 广东电网有限责任公司广州供电局 | Quantitative analysis system for virtual power plant operation benefits considering distributed resources |
CN113610435B (en) * | 2021-08-24 | 2024-06-07 | 内蒙古电力(集团)有限责任公司电力调度控制分公司 | Power system peak regulation method and device based on thermal power, self-contained power plant and energy storage |
-
2020
- 2020-12-08 CN CN202011443256.1A patent/CN112529622A/en active Pending
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112952918A (en) * | 2021-03-23 | 2021-06-11 | 国网山西省电力公司 | Electric power peak regulation and frequency modulation combined clearing method based on marginal benefit |
CN112952918B (en) * | 2021-03-23 | 2023-04-14 | 国网山西省电力公司 | Electric power peak regulation and frequency modulation combined clearing method based on marginal benefit |
CN113255957A (en) * | 2021-04-02 | 2021-08-13 | 国网山东省电力公司济南供电公司 | Quantitative optimization analysis method and system for uncertain factors of comprehensive service station |
CN113098040A (en) * | 2021-04-09 | 2021-07-09 | 国网新疆电力有限公司经济技术研究院 | Power grid side energy storage capacity optimal configuration method for obtaining multi-scene benefits |
CN113098040B (en) * | 2021-04-09 | 2023-07-11 | 国网新疆电力有限公司经济技术研究院 | Power grid side energy storage capacity optimal configuration method for obtaining multi-scene benefits |
CN113627762B (en) * | 2021-07-30 | 2024-03-22 | 国网山西省电力公司电力科学研究院 | Virtual power plant peak shaving method based on excitation electricity price |
CN113627762A (en) * | 2021-07-30 | 2021-11-09 | 国网山西省电力公司电力科学研究院 | Virtual power plant peak regulation method based on excitation electricity price |
CN113538066A (en) * | 2021-07-30 | 2021-10-22 | 国网上海市电力公司 | Control method, system, equipment and medium for virtual power plant to participate in power market |
CN113538066B (en) * | 2021-07-30 | 2024-02-27 | 国网上海市电力公司 | Control method, system, equipment and medium for virtual power plant to participate in power market |
CN113610435B (en) * | 2021-08-24 | 2024-06-07 | 内蒙古电力(集团)有限责任公司电力调度控制分公司 | Power system peak regulation method and device based on thermal power, self-contained power plant and energy storage |
CN113610435A (en) * | 2021-08-24 | 2021-11-05 | 内蒙古电力(集团)有限责任公司电力调度控制分公司 | Peak regulation method and device for power system based on thermal power, self-contained power plant and energy storage |
CN113837444A (en) * | 2021-08-26 | 2021-12-24 | 国网河南省电力公司经济技术研究院 | Trading clearing optimization method for virtual power plant participating in multiple markets |
CN113837444B (en) * | 2021-08-26 | 2023-12-26 | 国网河南省电力公司经济技术研究院 | Transaction clearing optimization method for virtual power plant participating in multiple markets |
CN114362227B (en) * | 2022-01-12 | 2022-09-09 | 华北电力大学 | Electric automobile frequency modulation reporting method considering user intention |
CN114362227A (en) * | 2022-01-12 | 2022-04-15 | 华北电力大学 | Electric automobile frequency modulation reporting method considering user intention |
CN115204443B (en) * | 2022-05-12 | 2023-07-18 | 电力规划总院有限公司 | Wind-solar-fire-storage combined market bidding optimal clear scheduling method |
CN115204443A (en) * | 2022-05-12 | 2022-10-18 | 电力规划总院有限公司 | Wind, light, fire and storage combined participation market bidding optimal strategy model |
CN116823332A (en) * | 2023-06-29 | 2023-09-29 | 广东电网有限责任公司广州供电局 | Quantitative analysis system for virtual power plant operation benefits considering distributed resources |
CN116823332B (en) * | 2023-06-29 | 2024-05-07 | 广东电网有限责任公司广州供电局 | Quantitative analysis system for virtual power plant operation benefits considering distributed resources |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112529622A (en) | Virtual power plant-based method for clearing multiple micro-micro main bodies participating in spot market | |
Mansour-Saatloo et al. | A risk-averse hybrid approach for optimal participation of power-to-hydrogen technology-based multi-energy microgrid in multi-energy markets | |
CN110188950B (en) | Multi-agent technology-based optimal scheduling modeling method for power supply side and demand side of virtual power plant | |
Damisa et al. | Microgrid energy and reserve management incorporating prosumer behind‐the‐meter resources | |
Luo et al. | Distributed peer-to-peer energy trading based on game theory in a community microgrid considering ownership complexity of distributed energy resources | |
Ju et al. | Robust purchase and sale transactions optimization strategy for electricity retailers with energy storage system considering two-stage demand response | |
Chen et al. | Research on day-ahead transactions between multi-microgrid based on cooperative game model | |
Parastegari et al. | Joint operation of wind farms and pump-storage units in the electricity markets: Modeling, simulation and evaluation | |
CN111311012A (en) | Multi-agent-based micro-grid power market double-layer bidding optimization method | |
Tavakoli et al. | Optimal probabilistic operation of energy hub with various energy converters and electrical storage based on electricity, heat, natural gas, and biomass by proposing innovative uncertainty modeling methods | |
Bagheri et al. | Stochastic optimization and scenario generation for peak load shaving in Smart District microgrid: sizing and operation | |
Chang et al. | Day-ahead bidding strategy of cloud energy storage serving multiple heterogeneous microgrids in the electricity market | |
Jin et al. | An overview of virtual power plant development from the perspective of market participation | |
Golla et al. | Evaluating the impact of regulation on the path of electrification in Citizen Energy Communities with prosumer investment | |
Peng et al. | Review on bidding strategies for renewable energy power producers participating in electricity spot markets | |
Liu et al. | Impact of industrial virtual power plant on renewable energy integration | |
Datta et al. | Energy management of multi-microgrids with renewables and electric vehicles considering price-elasticity based demand response: A bi-level hybrid optimization approach | |
Han et al. | A low-carbon dispatch strategy for power systems considering flexible demand response and energy storage | |
Kong et al. | Independence enhancement of distributed generation systems by integrating shared energy storage system and energy community with internal market | |
CN113765105A (en) | Micro-grid group energy management method based on dynamic random model | |
Alikhani et al. | Optimal implementation of consumer demand response program with consideration of uncertain generation in a microgrid | |
CN112580850A (en) | Clearing method and system for electric power peak regulation market | |
CN117114877A (en) | Medium-and-long-term power transaction method and system based on virtual power plant | |
CN117332937A (en) | Multi-energy complementary virtual power plant economic dispatching method considering demand response | |
Jiang et al. | Monthly electricity purchase and decomposition optimization considering wind power accommodation and day-ahead schedule |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210319 |
|
RJ01 | Rejection of invention patent application after publication |