CN110738377B - Day-ahead market clearing method and device containing wind power bidding, and computer equipment - Google Patents

Day-ahead market clearing method and device containing wind power bidding, and computer equipment Download PDF

Info

Publication number
CN110738377B
CN110738377B CN201911022977.2A CN201911022977A CN110738377B CN 110738377 B CN110738377 B CN 110738377B CN 201911022977 A CN201911022977 A CN 201911022977A CN 110738377 B CN110738377 B CN 110738377B
Authority
CN
China
Prior art keywords
node
wind turbine
real
market
power
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.)
Active
Application number
CN201911022977.2A
Other languages
Chinese (zh)
Other versions
CN110738377A (en
Inventor
张昆
周华锋
高红亮
朱文
顾慧杰
彭超逸
胡亚平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Southern Power Grid Co Ltd
Original Assignee
China Southern Power Grid Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Southern Power Grid Co Ltd filed Critical China Southern Power Grid Co Ltd
Priority to CN201911022977.2A priority Critical patent/CN110738377B/en
Publication of CN110738377A publication Critical patent/CN110738377A/en
Application granted granted Critical
Publication of CN110738377B publication Critical patent/CN110738377B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • 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/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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02B—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
    • Y02B10/00—Integration of renewable energy sources in buildings
    • Y02B10/30—Wind power
    • 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/10—Energy trading, including energy flowing from end-user application to grid
    • 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)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Finance (AREA)
  • Development Economics (AREA)
  • Accounting & Taxation (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Game Theory and Decision Science (AREA)
  • Health & Medical Sciences (AREA)
  • Tourism & Hospitality (AREA)
  • Quality & Reliability (AREA)
  • Data Mining & Analysis (AREA)
  • Operations Research (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The application relates to a method and a device for clearing a market in the day-ahead, a computer device and a storage medium. The method comprises the following steps: pre-establishing a day-ahead market clearing model containing wind power bidding; the day-ahead market clearing model comprises a first objective function taking the minimum market clearing total cost as a target and a corresponding first constraint condition; pre-establishing a real-time market clearing model containing wind power bidding; the real-time market clearing model comprises a second objective function taking the minimum load shedding expected cost as a target and a corresponding second constraint condition; and when the real power generation cost quoted by the conventional unit and the real wind power probability distribution function declared by the wind turbine unit, the day-ahead market clearing model and the real-time market clearing model achieve the aim of minimizing the total cost of market clearing. By applying the embodiment of the application, the income of each market member is the maximum under the condition that all market members declare respective real quotation.

Description

Day-ahead market clearing method and device containing wind power bidding, and computer equipment
Technical Field
The application relates to the field of power dispatching, in particular to a day-ahead market clearing method with wind power bidding, a day-ahead market clearing device with wind power bidding, computer equipment and a storage medium.
Background
As early as the nineties of the last century, many power markets abroad have achieved a shift from regulated monopolies to market competition. Among them, node Marginal Price (LMP) is introduced into power markets represented by the united states PJM power market, texas, california, new york, and new england. Under the node marginal electricity price mechanism, each market member declares respective quotation, and a system operator (ISO) makes market clearing decisions according to the quotation of the market member. Between 1990 and 2001, the united kingdom has imposed a consolidated power market. When bidding on the generator set in the market, strategic quotation exists, so that the benefit of the generator set is maximized, and the operation efficiency of the power market is reduced due to the strategic quotation, which is not beneficial to the efficient and healthy operation of the power market.
Most literature on electricity market bidding mechanisms assumes the convexity of the cost of electricity generation; in a large amount of literature on behavior analysis of a generator agent in a large-scale electric power market, a generator set is assumed to be priced or some numerical methods aiming at solving game theory balance, but most of the documents lack the research that a wind turbine generator directly participates in the market bidding condition in the day before.
Disclosure of Invention
Therefore, it is necessary to provide a method for clearing a day-ahead market including wind power bidding, a device for clearing a day-ahead market including wind power bidding, a computer device and a computer readable storage medium for solving the problem that the conventional market bidding mechanism, i.e., an LMP mechanism, has the problem that a generator set maximizes market strength to maximize benefits of the generator set, destroys market order, and causes economic efficiency loss, in view of the above technical problems.
A method of day-ahead marketing, comprising:
pre-establishing a day-ahead market clearing model containing wind power bidding; the day-ahead market clearing model comprises a first objective function taking the minimum market clearing total cost as a target and a corresponding first constraint condition;
pre-establishing a real-time market clearing model containing wind power bidding; the real-time market clearing model comprises a second objective function taking the minimum load shedding expected cost as a target and a corresponding second constraint condition;
and when the real power generation cost declared by the conventional unit and the real wind power probability distribution function declared by the wind turbine unit, the day-ahead market clearing model and the real-time market clearing model achieve the aim of minimizing the total cost of market clearing.
Preferably, the first and second electrodes are formed of a metal,
the first objective function is:
Figure BDA0002247812510000021
the first constraint includes:
and power balance constraint:
Figure BDA0002247812510000022
constraint of the power transmission network:
Figure BDA0002247812510000023
Figure BDA0002247812510000024
and (3) output constraint of a conventional unit:
Figure BDA0002247812510000025
output restraint of the wind turbine generator:
Figure BDA0002247812510000026
wherein n represents a node included in the power system; i represents the ith node;
Figure BDA0002247812510000027
representing the total cost of market clearing; c represents that self power generation cost quotations reported by all conventional units to a system operator ISO form a row vector;
Figure BDA0002247812510000028
representing wind power probability distribution functions predicted by quotations of all wind generation sets to form vectors; sigma fi(Pi) Representing the total power generation cost of the conventional unit;
Figure BDA0002247812510000029
representing the expected load shedding cost generated by the system operating in a real-time market due to the uncertainty of the output power of the wind turbine generator when the day-ahead market containing wind power bidding is cleared;
Figure BDA00022478125100000210
representing the output power of the wind turbine generator in the market at the node i in the day ahead; diRepresenting the load demand power of the node i; phi is equal to (phi)1,...,φn) Representing the power constraint upper limit of each node of the actual wind turbine generator in the system;
Figure BDA00022478125100000211
representing the upper limit of the output of the conventional unit; hliRepresenting a power generation transfer distribution factor of the node i to the transmission line l;
Figure BDA00022478125100000212
representing the transmission capacity of the transmission line l; pwind,iRepresenting the actual real-time power of the wind turbine at node i.
Preferably, the first and second electrodes are formed of a metal,
the second objective function is:
Figure BDA00022478125100000213
the second constraint includes:
and power balance constraint:
Figure BDA00022478125100000214
constraint of the power transmission network:
Figure BDA00022478125100000215
Figure BDA0002247812510000031
output restraint of the wind turbine generator:
Figure BDA0002247812510000032
real-time market load shedding constraint:
Figure BDA0002247812510000033
wind power curtailment of the wind turbine generator:
Figure BDA0002247812510000034
wherein the content of the first and second substances,
Figure BDA0002247812510000035
representing the wind curtailment power of the wind turbine generator at the node i;
Figure BDA0002247812510000036
representing the load shedding amount of the wind turbine system at the node i; b is the unit cost of the ISO shedding load of the system operator in the real-time market.
Preferably, if
Figure BDA0002247812510000037
In order to achieve the optimal solution,
under VCG mechanism, system operator ISO pays for conventional machine set at node i
Figure BDA0002247812510000038
Comprises the following steps:
Figure BDA0002247812510000039
under the VCG mechanism, the system operator ISO pays the cost of the wind turbine generator at the node i
Figure BDA00022478125100000310
Comprises the following steps:
Figure BDA00022478125100000311
wherein, P represents the column vector of the clear output of all the conventional generator sets of the power system in the market in the day ahead;
Figure BDA00022478125100000312
representing the optimal market clearing total cost when all market members participate in clearing;
Figure BDA00022478125100000313
the optimal market clearing total cost of the new system is shown when the conventional unit at the node i is not included;
Figure BDA00022478125100000314
representing the optimal market clearing total cost of the new system when the wind turbine set at the node i is not included;
Figure BDA00022478125100000315
and representing the probability distribution function of the wind power submitted by the other wind turbines.
Preferably, under the VCG mechanism, the incentive compatibility and the individual rationality of the conventional unit at least comprise that the clear profit produced by the conventional unit is maximum in the day before if and only if the conventional unit declares the real power generation cost;
for the conventional set at the node i, if all the wind power sets declare the probability distribution function of the real wind power, the other conventional sets except the conventional set at the node i declare the power generation cost as
Figure BDA00022478125100000316
In time, if the conventional unit at the node i reports the false power generation cost
Figure BDA00022478125100000317
The payment obtained is then:
Figure BDA00022478125100000318
wherein the content of the first and second substances,
Figure BDA00022478125100000319
the conventional unit at the node i does not participate in market clearing, and the declaration of power generation cost of all the conventional units of the new system
Figure BDA00022478125100000320
The market of the conventional unit is clear;
Figure BDA00022478125100000321
representing the total output of all conventional units of the system;
Figure BDA00022478125100000322
representing the false power generation cost declared by the conventional unit at the node i
Figure BDA00022478125100000323
The optimal market clearing;
if the conventional unit at the node i reports the real power generation cost CiThen the payment obtained is:
Figure BDA00022478125100000324
wherein, PiRepresenting the actual power generation cost C declared by the conventional unit at the node iiThe optimal market clearing;
when the conventional unit reports the false power generation cost at the node i
Figure BDA0002247812510000041
In time, the net profit is:
Figure BDA0002247812510000042
wherein the content of the first and second substances,
Figure BDA0002247812510000043
and
Figure BDA0002247812510000044
declaring false power generation cost for conventional unit at node i respectively
Figure BDA0002247812510000045
In time, the system optimizes the regular unit clearing at the scheduling node i
Figure BDA0002247812510000046
At their false generation costs, respectively
Figure BDA0002247812510000047
And true power generation cost CiThe cost of electricity generation;
Figure BDA0002247812510000048
the system is cleared for all market members to participate in, and the conventional unit at the node i declares the false power generation cost as
Figure BDA0002247812510000049
Reporting power generation cost by other conventional units except the conventional unit at the node i
Figure BDA00022478125100000410
The sum of the optimal power generation cost of other conventional units of the system;
when the conventional unit at the node i declares the real power generation cost CiIn time, the net profit is:
Figure BDA00022478125100000411
wherein f isi(Pi,Ci) Representing the reported power generation cost of other conventional units except the conventional unit at the node i
Figure BDA00022478125100000412
Actual power generation cost C reported by conventional unit at node iiIn time, the conventional unit at the market clearing node i clears PiAt its true power generation cost CiThe cost of electricity generation;
Figure BDA00022478125100000413
shows that all market members participate in clearing, and other conventional units declare the power generation cost
Figure BDA00022478125100000414
Actual power generation cost C reported by conventional unit at node iiThe sum of the optimal power generation cost of other conventional units except the conventional unit at the node i;
in the above formula for calculating net profit, the first term at the right end of the two equations is irrelevant to the power generation cost declared by the conventional unit at the node i; and P isiAnd
Figure BDA00022478125100000415
declaring real power generation cost C for conventional units at node i respectivelyiAnd other conventional units report the power generation cost
Figure BDA00022478125100000416
The optimal clearing plan of the conventional generator set of the system is as follows:
Figure BDA00022478125100000417
therefore, the conventional unit at the node i declares the real power generation cost CiThe obtained net profit is not less than the declared false power generation cost
Figure BDA00022478125100000418
Net profit, so there are:
Figure BDA00022478125100000419
when the conventional unit at the node i declares the real power generation cost, the excitation compatibility of the conventional unit;
for conventional unitsPhysical, when the conventional unit at node i declares the real generating cost CiWhen the temperature of the water is higher than the set temperature,
Figure BDA0002247812510000051
equivalent to the increase of P in the clearing of all the conventional units participating in the marketiConstraint of 0; therefore, the optimization of the daily output model can be reduced, the objective function value is not less than that of all conventional units participating in the daily market output model, so that the method comprises the following steps:
Figure BDA0002247812510000052
preferably, under the VCG mechanism, the excitation compatibility and the individual rationality of the wind turbine generator at least comprise that if and only if the wind turbine generator declares a real wind power probability distribution function, the clear profit of the wind turbine generator is maximum in the day ahead;
for the wind turbine generator set at the node i, if all conventional wind turbine generator sets claim the real power generation cost, when other wind turbine generator sets claim the wind power probability distribution function
Figure BDA0002247812510000053
If the wind turbine generator at the node i declares the probability distribution function of false wind power as
Figure BDA0002247812510000054
The payment obtained is then:
Figure BDA0002247812510000055
wherein the content of the first and second substances,
Figure BDA0002247812510000056
representing wind turbine generator based on node i
The probability distribution function of the virtual reporting wind power is
Figure BDA0002247812510000057
Reporting wind of other wind turbine generatorsElectric power probability distribution function
Figure BDA0002247812510000058
The load shedding amount of the system;
if the wind turbine generator at the node i declares the probability distribution function of the real wind power as
Figure BDA0002247812510000059
The payment obtained is then:
Figure BDA00022478125100000510
wherein the content of the first and second substances,
Figure BDA00022478125100000511
representing a probability distribution function of declaring real wind power based on a wind turbine at node i
Figure BDA00022478125100000512
Reporting wind power probability distribution function of other wind turbine generators
Figure BDA00022478125100000513
The load shedding amount of the system;
in the payment calculation formula, the first term at the right end of the equation is irrelevant to the wind power probability distribution function declared by the wind turbine at the node i; while
Figure BDA00022478125100000514
Declaring real wind power probability distribution function for wind turbine generator at node i
Figure BDA00022478125100000515
Reporting wind power probability distribution function of other wind turbine generators
Figure BDA00022478125100000516
The optimal load shedding of the system is determined by the time and the time
Figure BDA00022478125100000517
Therefore, the method comprises the following steps:
Figure BDA00022478125100000518
therefore, the wind turbine generator set at the node i declares a real wind power probability distribution function
Figure BDA00022478125100000519
The obtained payment is not less than the probability distribution function of declaring false wind power
Figure BDA00022478125100000520
Payment when, that is:
Figure BDA00022478125100000521
the wind turbine generator at the node i declares a real wind power probability distribution function, namely the excitation compatibility of the wind turbine generator;
for the individual rationality of the wind turbine,
Figure BDA0002247812510000061
equivalent to the situation that all wind turbine generators participate in market clearing, the number of the wind turbine generators is increased
Figure BDA0002247812510000062
Therefore, the optimization feasible region of the model of the day-ahead output is reduced, and the objective function value is not less than that of the model of the day-ahead market output of all wind turbines, so that the method comprises the following steps:
Figure BDA0002247812510000063
preferably, under a VCG mechanism, the conventional unit declares the real power generation cost, the wind turbine unit declares the real wind power probability distribution function, and the market clearing total cost is minimized.
A day-ahead market dispensing apparatus, comprising:
the first model establishing module is used for pre-establishing a day-ahead market clearing model containing wind power bidding; the day-ahead market clearing model comprises a first objective function taking the minimum market clearing total cost as a target and a corresponding first constraint condition;
the second model establishing module is used for pre-establishing a real-time market clearing model containing wind power bidding; the real-time market clearing model comprises a second objective function taking the minimum load shedding expected cost as a target and a corresponding second constraint condition;
and when the real power generation cost declared by the conventional unit and the real wind power probability distribution function declared by the wind turbine unit, the day-ahead market clearing model and the real-time market clearing model achieve the aim of minimizing the total cost of market clearing.
A computer device comprising a memory and a processor, the memory storing a computer program, the processor implementing the following steps when executing the computer program:
pre-establishing a day-ahead market clearing model containing wind power bidding; the day-ahead market clearing model comprises a first objective function taking the minimum market clearing total cost as a target and a corresponding first constraint condition;
pre-establishing a real-time market clearing model containing wind power bidding; the real-time market clearing model comprises a second objective function taking the minimum load shedding expected cost as a target and a corresponding second constraint condition;
and when the real power generation cost declared by the conventional unit and the real wind power probability distribution function declared by the wind turbine unit, the day-ahead market clearing model and the real-time market clearing model achieve the aim of minimizing the total cost of market clearing.
A computer device comprising a memory and a processor, the memory storing a computer program, the processor implementing the following steps when executing the computer program:
pre-establishing a day-ahead market clearing model containing wind power bidding; the day-ahead market clearing model comprises a first objective function taking the minimum market clearing total cost as a target and a corresponding first constraint condition;
pre-establishing a real-time market clearing model containing wind power bidding; the real-time market clearing model comprises a second objective function taking the minimum load shedding expected cost as a target and a corresponding second constraint condition;
and when the real power generation cost declared by the conventional unit and the real wind power probability distribution function declared by the wind turbine unit, the day-ahead market clearing model and the real-time market clearing model achieve the aim of minimizing the total cost of market clearing.
A computer-readable storage medium, on which a computer program is stored which, when executed by a processor, carries out the steps of:
pre-establishing a day-ahead market clearing model containing wind power bidding; the day-ahead market clearing model comprises a first objective function taking the minimum market clearing total cost as a target and a corresponding first constraint condition;
pre-establishing a real-time market clearing model containing wind power bidding; the real-time market clearing model comprises a second objective function taking the minimum load shedding expected cost as a target and a corresponding second constraint condition;
and when the real power generation cost declared by the conventional unit and the real wind power probability distribution function declared by the wind turbine unit, the day-ahead market clearing model and the real-time market clearing model achieve the aim of minimizing the total cost of market clearing.
According to the day-ahead market clearing method, the day-ahead market clearing device, the computer equipment and the computer readable storage medium, the day-ahead market clearing model containing the wind power bidding and the real-time market clearing model containing the wind power bidding are established in advance, and when the real power generation cost quoted by the conventional unit and the real wind power probability distribution function quoted by the wind power unit are quoted, the day-ahead market clearing model and the real-time market clearing model reach the aim of minimizing the total cost of the market clearing. By applying the embodiment of the application, under the condition that all market members declare respective real quotation, the income of each market member is the maximum; meanwhile, the overall economic benefit of the market is improved through theoretical analysis and example analysis, so that the superiority of the VCG mechanism in the aspects of market member incentive compatibility, individual rationality and economic efficiency of the market is realized.
Drawings
FIG. 1 is a schematic flow chart of a method for exporting a day-ahead market including wind power bidding according to an embodiment;
FIG. 2 is a modified IEEE14 node system of one embodiment;
FIG. 3 is a table diagram of parameters of a genset of an embodiment;
FIG. 4 is a table diagram of parameters of a load for one embodiment;
FIG. 5 is a tabular illustration of output, system payments, and net profits for a generator set declaring a true quote, according to an embodiment;
6-7 are graphs of net profit levels for a conventional unit declaring different power generation cost factors, according to one embodiment;
8-9 are diagrams of net profit levels for a wind turbine generator declaring different wind power probability distribution parameters, under an embodiment;
FIG. 10 is a block diagram illustrating a day-ahead market clearing apparatus with wind power bidding according to an embodiment;
FIG. 11 is an internal block diagram of a computer device of an embodiment.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application more apparent, the present application is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the present application and are not intended to limit the present application.
In one embodiment, as shown in fig. 1, a flow diagram of a day-ahead market clearing method with wind power bidding is provided. The method can be applied to various computer devices such as computers and notebook computers, is suitable for the clearing of the day-ahead market containing wind power bidding of various multi-node power systems, and applies the Vickrey-Clarke-Groves (VCG) market bidding mechanism to the design of the day-ahead market clearing containing wind power bidding, so that the best economic benefit is achieved when the day-ahead market clearing is beneficial to a conventional generator set (conventional generator set) and a wind generator set to submit real quotes. Specifically, the market clearing method of the present application may include the steps of:
step S110, pre-establishing a day-ahead market clearing model containing wind power bidding; the day-ahead market clearing model comprises a first objective function and a corresponding first constraint condition, wherein the first objective function takes minimum market clearing total cost as a target.
Step S120, a real-time market clearing model containing wind power bidding is established in advance; wherein the real-time market clearing model comprises a second objective function and a corresponding second constraint condition, wherein the second objective function is used for targeting the minimum load shedding expected cost.
The utility model provides a mathematical model of market clearing, including the day-ahead market clearing model and the real-time market clearing model that contain wind-powered electricity generation competitive bidding:
the day-ahead market clearing model containing wind power bidding is as follows:
a first objective function:
Figure BDA0002247812510000091
the first constraint condition is:
1) and power balance constraint:
Figure BDA0002247812510000092
2) constraint of the power transmission network:
Figure BDA0002247812510000093
Figure BDA0002247812510000094
3) and (3) output constraint of a conventional unit:
Figure BDA0002247812510000095
4) output restraint of the wind turbine generator:
Figure BDA0002247812510000096
wherein n represents a node included in the power system; i represents the ith node;
Figure BDA0002247812510000097
representing the total cost of market clearing; c represents that self power generation cost quotations reported by all conventional units to a system operator ISO form a row vector;
Figure BDA0002247812510000098
representing wind power probability distribution functions predicted by quotations of all wind generation sets to form vectors; sigma fi(Pi) Representing the total power generation cost of the conventional unit;
Figure BDA0002247812510000099
representing the expected load shedding cost generated by the system operating in a real-time market due to the uncertainty of the output power of the wind turbine generator when the day-ahead market containing wind power bidding is cleared;
Figure BDA00022478125100000910
representing the output power of the wind turbine generator in the market at the node i in the day ahead; diRepresenting the load demand power of the node i; phi is equal to (phi)1,...,φn) Representing actual wind power of each node in systemThe upper limit of the unit power constraint;
Figure BDA00022478125100000911
representing the upper limit of the output of the conventional unit; hliRepresenting a power generation transfer distribution factor of the node i to the transmission line l;
Figure BDA00022478125100000912
representing the transmission capacity of the transmission line l; pwind,iRepresenting the actual real-time power of the wind turbine at node i.
The real-time market clearing model containing wind power bidding is as follows:
a second objective function:
Figure BDA0002247812510000101
the second constraint condition is as follows:
1) and power balance constraint:
Figure BDA0002247812510000102
2) constraint of the power transmission network:
Figure BDA0002247812510000103
Figure BDA0002247812510000104
3) output restraint of the wind turbine generator:
Figure BDA0002247812510000105
4) real-time market load shedding constraint:
Figure BDA0002247812510000106
5) wind power curtailment of the wind turbine generator:
Figure BDA0002247812510000107
wherein the content of the first and second substances,
Figure BDA0002247812510000108
representing the wind curtailment power of the wind turbine generator at the node i;
Figure BDA0002247812510000109
representing the load shedding amount of the wind turbine system at the node i; b is the unit cost of the ISO shedding load of the system operator in the real-time market.
And S130, when the real power generation cost quoted by the conventional unit and the real wind power probability distribution function declared by the wind turbine unit are used, the day-ahead market clearing model and the real-time market clearing model achieve the aim of minimizing the total cost of market clearing.
Suppose that
Figure BDA00022478125100001010
For the optimal solution, under the VCG mechanism, the system operator ISO pays the cost of the conventional unit at the node i
Figure BDA00022478125100001011
Comprises the following steps:
Figure BDA00022478125100001012
under the VCG mechanism, the system operator ISO pays the cost of the wind turbine generator at the node i
Figure BDA00022478125100001013
Comprises the following steps:
Figure BDA00022478125100001014
wherein, P represents the column vector of the clear output of all the conventional generator sets of the power system in the market in the day ahead;
Figure BDA00022478125100001015
representing the optimal market clearing total cost when all market members participate in clearing;
Figure BDA00022478125100001016
the optimal market clearing total cost of the new system is shown when the conventional unit at the node i is not included;
Figure BDA00022478125100001017
representing the optimal market clearing total cost of the new system when the wind turbine set at the node i is not included;
Figure BDA0002247812510000111
and representing the probability distribution function of the wind power submitted by the other wind turbines.
The above-mentioned formulas (1) to (15) represent a conventional, market-clearing model as well as a VCG mechanism model.
In order to more clearly explain the embodiments of the present application, the following description will be made using fig. 2 to 9.
(one) adopt IEEE14 node system of fig. 2 modification, there are: 4 conventional units which are respectively positioned at the node 1, the node 2, the node 3 and the node 6; 1 wind power generator set located at node 8; and 3 load nodes are respectively positioned at the node 9, the node 13 and the node 14. The parameters of the genset and load are shown in the tables of fig. 3 and 4.
For the conventional set at the node i, assuming that all the wind power sets report a real wind power probability distribution function, when other conventional sets except the conventional set at the node i report the power generation cost
Figure BDA0002247812510000112
In time, the generator set i can choose whether to declare the true cost of power generation. If it declares false power generation cost
Figure BDA0002247812510000113
The payment obtained is then:
Figure BDA0002247812510000114
wherein the content of the first and second substances,
Figure BDA0002247812510000115
the conventional unit at the node i does not participate in market clearing, and the declaration of power generation cost of all the conventional units of the new system
Figure BDA0002247812510000116
The market of the conventional unit is clear;
Figure BDA0002247812510000117
representing the total output of all conventional units of the system;
Figure BDA0002247812510000118
representing the false power generation cost declared by the conventional unit at the node i
Figure BDA0002247812510000119
The market is best cleared.
If the conventional unit at the node i reports the real power generation cost CiThen the payment obtained is:
Figure BDA00022478125100001110
wherein, PiRepresenting the actual power generation cost C declared by the conventional unit at the node iiThe market is best cleared.
For a conventional unit, the goal of participating in market clearing is to maximize its own net profit, i.e., the payment made by the conventional unit minus the unit's own cost of generating electricity. Therefore, when the conventional unit i gives a false report of the self power generation cost
Figure BDA00022478125100001111
In time, the net profit is:
Figure BDA00022478125100001112
wherein the content of the first and second substances,
Figure BDA00022478125100001113
and
Figure BDA00022478125100001114
declaring false power generation cost for conventional unit at node i respectively
Figure BDA00022478125100001115
In time, the system optimizes the regular unit clearing at the scheduling node i
Figure BDA00022478125100001116
At their false generation costs, respectively
Figure BDA00022478125100001117
And true power generation cost CiThe cost of electricity generation;
Figure BDA00022478125100001118
the system is cleared for all market members to participate in, and the conventional unit at the node i declares the false power generation cost as
Figure BDA00022478125100001119
Reporting power generation cost by other conventional units except the conventional unit at the node i
Figure BDA0002247812510000121
The sum of the optimal power generation cost of other conventional units of the system;
when the conventional unit at the node i declares the real power generation cost CiIn time, the net profit is:
Figure BDA0002247812510000122
wherein f isi(Pi,Ci) Representing the reported power generation cost of other conventional units except the conventional unit at the node i
Figure BDA0002247812510000123
Actual power generation cost C reported by conventional unit at node iiIn time, the conventional unit at the market clearing node i clears PiAt its true power generation cost CiThe cost of electricity generation;
Figure BDA0002247812510000124
shows that all market members participate in clearing, and other conventional units declare the power generation cost
Figure BDA0002247812510000125
Actual power generation cost C reported by conventional unit at node iiAnd the sum of the optimal power generation cost of other conventional units except the conventional unit at the node i.
For the above equations (18) and (19), the first term at the right end of the equation is irrelevant to the electricity generation cost quoted by the conventional unit at the node i; and P isiAnd
Figure BDA0002247812510000126
declaring real power generation cost C for conventional units at node i respectivelyiAnd other conventional units report the power generation cost
Figure BDA0002247812510000127
The optimal clearing plan of the conventional generator set of the system is as follows:
Figure BDA0002247812510000128
therefore, the conventional unit at the node i declares the real power generation cost CiThe obtained net profit is not less than the declared false power generation cost
Figure BDA0002247812510000129
Net profit on time, i.e.:
Figure BDA00022478125100001210
based on the result of the formula (21), the conventional unit at the node i reports the real power generation cost, namely the excitation compatibility of the conventional unit. For the individual rationality of a conventional unit, according to equation (19),
Figure BDA00022478125100001211
equivalently setting the output of the conventional unit at the node i to be 0 in the model of all the conventional units participating in market clearing, namely adding P in the clearing modeliConstraint 0. Therefore, the optimization feasible domain of the output model is reduced, and the objective function value is not less than that of all conventional units participating in the output of the market at present. Thus:
Figure BDA00022478125100001212
(II) for the wind turbine generator at the node i, assuming that all conventional wind turbine generators report real power generation cost, when other wind turbine generators report a wind power probability distribution function
Figure BDA00022478125100001213
And then, the wind turbine generator at the node i can select whether to declare a real wind power probability distribution function. If the probability distribution function of the wind power of the self-body is false report of the probability distribution function of the wind power of the self-body is
Figure BDA00022478125100001214
(because system operator ISO selects wind turbine generator system to go out the market, according to the principle that the power generation power is big priority, therefore the false report of this application means that wind turbine generator system exaggerates self wind power probability distribution in order to obtain more generated power), then the VCG payment that obtains is:
Figure BDA0002247812510000131
wherein the content of the first and second substances,
Figure BDA0002247812510000132
the probability distribution function of the wind power based on the false report of the wind turbine generator at the node i is represented as
Figure BDA0002247812510000133
Reporting wind power probability distribution function of other wind turbine generators
Figure BDA0002247812510000134
The load shedding amount of the system;
if the wind turbine generator at the node i declares the probability distribution function of the real wind power as
Figure BDA0002247812510000135
Then the VCG payment obtained is:
Figure BDA0002247812510000136
wherein the content of the first and second substances,
Figure BDA0002247812510000137
representing a probability distribution function of declaring real wind power based on a wind turbine at node i
Figure BDA0002247812510000138
Reporting wind power probability distribution function of other wind turbine generators
Figure BDA0002247812510000139
The load shedding amount of the system;
for the above equations (23) and (24), the first term at the right end of the equation is irrelevant to the probability distribution function of the wind power declared by the wind turbine at the node i; while
Figure BDA00022478125100001310
Is at node iProbability distribution function for reporting real wind power of wind turbine generator
Figure BDA00022478125100001311
Reporting wind power probability distribution function of other wind turbine generators
Figure BDA00022478125100001312
The optimal load shedding of the system is determined by the time and the time
Figure BDA00022478125100001313
Therefore, the method comprises the following steps:
Figure BDA00022478125100001314
therefore, the wind turbine generator set at the node i declares a real wind power probability distribution function
Figure BDA00022478125100001315
The payment obtained is not less than the payment of the virtual strike, i.e.:
Figure BDA00022478125100001316
based on the result of the formula (26), the wind turbine at the node i declares a true wind power probability distribution function, which is the excitation compatibility of the wind turbine. For the individual rationality of the wind turbine, according to equation (24),
Figure BDA00022478125100001317
equivalently setting the wind power probability distribution function of the wind turbine generator i to be 0 in the model of all the wind turbine generators participating in market clearing, namely adding the probability distribution function to the clearing model
Figure BDA00022478125100001318
A constraint condition. Therefore, the optimization feasible region of the output model is reduced, and the objective function value is not less than that of all wind turbines participating in the output of the market in the day before. Thus:
Figure BDA00022478125100001319
fig. 5 is a schematic diagram of the output, system payment and net profit when the generator set declares the real power generation cost and the wind power probability distribution function in the modified IEEE14 node system. Fig. 6-7 are schematic diagrams of net profit levels when different power generation cost coefficients are declared by a conventional wind turbine generator, and fig. 8-9 are schematic diagrams of net profit levels when different wind power probability distribution parameters are declared by a wind turbine generator.
And (III) under the condition that the conventional generator set and the wind generating set meet excitation compatibility and individuality, the generator set voluntarily participates in the system day-ahead market, and reports the real generating cost and the wind power probability distribution function. Therefore, the minimization of the total cost of market clearing is automatically satisfied by the above-mentioned (one) and (two).
The method provides an incentive wind power competitive bidding day-ahead spot market clearing mechanism, namely under a VCG mechanism, the market capacity of all market members is eliminated, the node power price is not influenced by the decision of a single market member any more, and the income of the single market member is not influenced by the decision of other market members any more; under the condition that all members declare respective real quotation, the income of each market member is the maximum; meanwhile, the overall economic benefit of the market is improved through theoretical analysis and example analysis, so that the superiority of the VCG mechanism in the aspects of market member incentive compatibility, individual rationality and economic efficiency of the market is realized.
It should be understood that, although the steps in the flowchart of fig. 1 are shown in order as indicated by the arrows, the steps are not necessarily performed in order as indicated by the arrows. The steps are not performed in the exact order shown and described, and may be performed in other orders, unless explicitly stated otherwise. Moreover, at least a portion of the steps in fig. 1 may include multiple sub-steps or multiple stages that are not necessarily performed at the same time, but may be performed at different times, and the order of performance of the sub-steps or stages is not necessarily sequential, but may be performed in turn or alternately with other steps or at least a portion of the sub-steps or stages of other steps.
In one embodiment, as shown in fig. 10, there is provided a day-ahead market clearing device including a wind power bidding, including:
the first model establishing module 210 is used for establishing a day-ahead market clearing model containing wind power bidding in advance; the day-ahead market clearing model comprises a first objective function taking the minimum market clearing total cost as a target and a corresponding first constraint condition;
the second model establishing module 220 is used for pre-establishing a real-time market clearing model containing wind power bidding; the real-time market clearing model comprises a second objective function taking the minimum load shedding expected cost as a target and a corresponding second constraint condition;
the clearing module 230 is configured to, when a real power generation cost quoted by the conventional wind turbine generator and a real wind power probability distribution function declared by the wind turbine generator, enable the day-ahead market clearing model and the real-time market clearing model to achieve a goal of minimizing a total cost of market clearing.
Preferably, the first objective function is:
Figure BDA0002247812510000151
the first constraint includes:
and power balance constraint:
Figure BDA0002247812510000152
constraint of the power transmission network:
Figure BDA0002247812510000153
Figure BDA0002247812510000154
and (3) output constraint of a conventional unit:
Figure BDA0002247812510000155
output restraint of the wind turbine generator:
Figure BDA0002247812510000156
wherein n represents a node included in the power system; i represents the ith node;
Figure BDA0002247812510000157
representing the total cost of market clearing; c represents that self power generation cost quotations reported by all conventional units to a system operator ISO form a row vector;
Figure BDA0002247812510000158
representing wind power probability distribution functions predicted by quotations of all wind generation sets to form vectors; sigma fi(Pi) Representing the total power generation cost of the conventional unit;
Figure BDA0002247812510000159
representing the expected load shedding cost generated by the system operating in a real-time market due to the uncertainty of the output power of the wind turbine generator when the day-ahead market containing wind power bidding is cleared;
Figure BDA00022478125100001510
representing the output power of the wind turbine generator in the market at the node i in the day ahead; diRepresenting the load demand power of the node i; phi is equal to (phi)1,...,φn) Representing the power constraint upper limit of each node of the actual wind turbine generator in the system;
Figure BDA00022478125100001511
representing the upper limit of the output of the conventional unit; hliRepresenting a power generation transfer distribution factor of the node i to the transmission line l;
Figure BDA00022478125100001512
representing the transmission capacity of the transmission line l; pwind,iRepresenting the actual real-time power of the wind turbine at node i.
Preferably, the second objective function is:
Figure BDA00022478125100001513
the second constraint includes:
and power balance constraint:
Figure BDA00022478125100001514
constraint of the power transmission network:
Figure BDA00022478125100001515
Figure BDA00022478125100001516
output restraint of the wind turbine generator:
Figure BDA00022478125100001517
real-time market load shedding constraint:
Figure BDA00022478125100001518
wind power curtailment of the wind turbine generator:
Figure BDA00022478125100001519
wherein the content of the first and second substances,
Figure BDA00022478125100001520
representing the wind curtailment power of the wind turbine generator at the node i;
Figure BDA00022478125100001521
representing the load shedding amount of the wind turbine system at the node i; b is real timeIn the market, the system operator ISO cuts the unit cost of the load.
Preferably, if
Figure BDA0002247812510000161
In order to achieve the optimal solution,
under VCG mechanism, system operator ISO pays for conventional machine set at node i
Figure BDA0002247812510000162
Comprises the following steps:
Figure BDA0002247812510000163
under the VCG mechanism, the system operator ISO pays the cost of the wind turbine generator at the node i
Figure BDA0002247812510000164
Comprises the following steps:
Figure BDA0002247812510000165
wherein, P represents the column vector of the clear output of all the conventional generator sets of the power system in the market in the day ahead;
Figure BDA0002247812510000166
representing the optimal market clearing total cost when all market members participate in clearing;
Figure BDA0002247812510000167
the optimal market clearing total cost of the new system is shown when the conventional unit at the node i is not included;
Figure BDA0002247812510000168
representing the optimal market clearing total cost of the new system when the wind turbine set at the node i is not included;
Figure BDA0002247812510000169
indicating the remaining windAnd the wind power probability distribution function submitted by the generator set.
Preferably, under the VCG mechanism, the incentive compatibility and the individual rationality of the conventional unit at least comprise that the clear profit produced by the conventional unit before the day is maximum if and only if the conventional unit declares the real power generation cost quotation;
for the conventional set at the node i, if all the wind power sets declare the probability distribution function of the real wind power, the other conventional sets except the conventional set at the node i declare the power generation cost as
Figure BDA00022478125100001610
In time, if the conventional unit at the node i reports the false power generation cost
Figure BDA00022478125100001611
The payment obtained is then:
Figure BDA00022478125100001612
wherein the content of the first and second substances,
Figure BDA00022478125100001613
the conventional unit at the node i does not participate in market clearing, and the declaration of power generation cost of all the conventional units of the new system
Figure BDA00022478125100001614
The market of the conventional unit is clear;
Figure BDA00022478125100001615
representing the total output of all conventional units of the system;
Figure BDA00022478125100001616
representing the false power generation cost declared by the conventional unit at the node i
Figure BDA00022478125100001617
The optimal market clearing;
if the conventional unit at the node i declares the truthCost of electricity generation CiThen the payment obtained is:
Figure BDA00022478125100001618
wherein, PiRepresenting the actual power generation cost C declared by the conventional unit at the node iiThe optimal market clearing;
when the conventional unit reports the false power generation cost at the node i
Figure BDA00022478125100001619
In time, the net profit is:
Figure BDA00022478125100001620
wherein the content of the first and second substances,
Figure BDA0002247812510000171
and
Figure BDA0002247812510000172
declaring false power generation cost for conventional unit at node i respectively
Figure BDA0002247812510000173
In time, the system optimizes the regular unit clearing at the scheduling node i
Figure BDA0002247812510000174
At their false generation costs, respectively
Figure BDA0002247812510000175
And true power generation cost CiThe cost of electricity generation;
Figure BDA0002247812510000176
the system is cleared for all market members to participate in, and the conventional unit at the node i declares the false power generation cost as
Figure BDA0002247812510000177
Reporting power generation cost by other conventional units except the conventional unit at the node i
Figure BDA0002247812510000178
The sum of the optimal power generation cost of other conventional units of the system;
when the conventional unit at the node i declares the real power generation cost CiIn time, the net profit is:
Figure BDA0002247812510000179
wherein f isi(Pi,Ci) Representing the reported power generation cost of other conventional units except the conventional unit at the node i
Figure BDA00022478125100001710
Actual power generation cost C reported by conventional unit at node iiIn time, the conventional unit at the market clearing node i clears PiAt its true power generation cost CiThe cost of electricity generation;
Figure BDA00022478125100001711
shows that all market members participate in clearing, and other conventional units declare the power generation cost
Figure BDA00022478125100001712
Actual power generation cost C reported by conventional unit at node iiThe sum of the optimal power generation cost of other conventional units except the conventional unit at the node i;
in the above formula for calculating net profit, the first term at the right end of the two equations is irrelevant to the power generation cost declared by the conventional unit at the node i; and P isiAnd
Figure BDA00022478125100001713
declaring real power generation cost C for conventional units at node i respectivelyiAnd other conventional units report the power generation cost
Figure BDA00022478125100001714
The optimal clearing plan of the conventional generator set of the system is as follows:
Figure BDA00022478125100001715
therefore, the conventional unit at the node i declares the real power generation cost CiThe obtained net profit is not less than the declared false power generation cost
Figure BDA00022478125100001716
Net profit, so there are:
Figure BDA00022478125100001717
when the conventional unit at the node i declares the real power generation cost, the excitation compatibility of the conventional unit;
for the individual rationality of the conventional unit, when the conventional unit at the node i declares the real power generation cost CiWhen the temperature of the water is higher than the set temperature,
Figure BDA00022478125100001718
equivalent to the increase of P in the clearing of all the conventional units participating in the marketiConstraint of 0; therefore, the optimization of the daily output model can be reduced, the objective function value is not less than that of all conventional units participating in the daily market output model, so that the method comprises the following steps:
Figure BDA00022478125100001719
preferably, under the VCG mechanism, the excitation compatibility and the individual rationality of the wind turbine generator at least comprise that if and only if the wind turbine generator declares a real wind power probability distribution function, the clear profit of the wind turbine generator is maximum in the day ahead;
for the wind turbine generator set at the node i, if all conventional wind turbine generator sets claim the real power generation cost, when other wind turbine generator sets claim the wind power probability distribution functionNumber of
Figure BDA0002247812510000181
If the wind turbine generator at the node i declares the probability distribution function of false wind power as
Figure BDA0002247812510000182
The payment obtained is then:
Figure BDA0002247812510000183
wherein the content of the first and second substances,
Figure BDA0002247812510000184
representing wind turbine generator based on node i
The probability distribution function of the virtual reporting wind power is
Figure BDA0002247812510000185
Reporting wind power probability distribution function of other wind turbine generators
Figure BDA0002247812510000186
The load shedding amount of the system;
if the wind turbine generator at the node i declares the probability distribution function of the real wind power as
Figure BDA0002247812510000187
The payment obtained is then:
Figure BDA0002247812510000188
wherein the content of the first and second substances,
Figure BDA0002247812510000189
representing a probability distribution function of declaring real wind power based on a wind turbine at node i
Figure BDA00022478125100001810
Declaration of other wind turbine generatorsWind power probability distribution function
Figure BDA00022478125100001811
The load shedding amount of the system;
in the payment calculation formula, the first term at the right end of the equation is irrelevant to the wind power probability distribution function declared by the wind turbine at the node i; while
Figure BDA00022478125100001812
Declaring real wind power probability distribution function for wind turbine generator at node i
Figure BDA00022478125100001813
Reporting wind power probability distribution function of other wind turbine generators
Figure BDA00022478125100001814
The optimal load shedding of the system is determined by the time and the time
Figure BDA00022478125100001815
Therefore, the method comprises the following steps:
Figure BDA00022478125100001816
therefore, the wind turbine generator set at the node i declares a real wind power probability distribution function
Figure BDA00022478125100001817
The obtained payment is not less than the probability distribution function of declaring false wind power
Figure BDA00022478125100001818
Payment when, that is:
Figure BDA00022478125100001819
the wind turbine generator at the node i declares a real wind power probability distribution function, namely the excitation compatibility of the wind turbine generator;
for the individual rationality of the wind turbine,
Figure BDA00022478125100001820
equivalent to the situation that all wind turbine generators participate in market clearing, the number of the wind turbine generators is increased
Figure BDA00022478125100001821
Therefore, the optimization feasible region of the model of the day-ahead output is reduced, and the objective function value is not less than that of the model of the day-ahead market output of all wind turbines, so that the method comprises the following steps:
Figure BDA00022478125100001822
preferably, under a VCG mechanism, the conventional unit declares a true power generation cost quotation, the wind turbine unit declares a true wind power probability distribution function, and the market clearing total cost is minimized, and the method at least comprises the steps that under the condition that the conventional unit and the wind turbine unit meet excitation compatibility and individuality, the conventional unit declares the true power generation cost and the wind turbine unit declares the true wind power probability distribution function, and the aim of minimizing the market clearing total cost is automatically met by a day-ahead market clearing model.
For specific limitations of the day-ahead market clearing device, reference may be made to the limitations of the day-ahead market clearing method above, which are not described in detail herein. The modules in the market clearing device can be wholly or partially realized by software, hardware and a combination thereof. The modules can be embedded in a hardware form or independent from a processor in the computer device, and can also be stored in a memory in the computer device in a software form, so that the processor can call and execute operations corresponding to the modules.
The day-ahead market clearing device provided by the above can be used for executing the day-ahead market clearing method provided by any of the above embodiments, and has corresponding functions and beneficial effects.
In one embodiment, a computer device is provided, which may be a terminal, and its internal structure diagram may be as shown in fig. 11. The computer device includes a processor, a memory, a network interface, a display screen, and an input device connected by a system bus. Wherein the processor of the computer device is configured to provide computing and control capabilities. The memory of the computer device comprises a nonvolatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of an operating system and computer programs in the non-volatile storage medium. The network interface of the computer device is used for communicating with an external terminal through a network connection. The computer program is executed by a processor to implement a method of indoor positioning of an air sensor. The display screen of the computer equipment can be a liquid crystal display screen or an electronic ink display screen, and the input device of the computer equipment can be a touch layer covered on the display screen, a key, a track ball or a touch pad arranged on the shell of the computer equipment, an external keyboard, a touch pad or a mouse and the like.
Those skilled in the art will appreciate that the architecture shown in fig. 11 is merely a block diagram of some of the structures associated with the disclosed aspects and is not intended to limit the computing devices to which the disclosed aspects apply, as particular computing devices may include more or less components than those shown, or may combine certain components, or have a different arrangement of components.
In one embodiment, a computer device is provided, comprising a memory in which a computer program is stored and a processor which, when executing the computer program, implements the method of the day ahead market clearing described above.
In one embodiment, a computer-readable storage medium is provided, on which a computer program is stored, which, when being executed by a processor, carries out the method of the day-ahead market clearing described above.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by hardware instructions of a computer program, which can be stored in a non-volatile computer-readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. Any reference to memory, storage, database, or other medium used in the embodiments provided herein may include non-volatile and/or volatile memory, among others. Non-volatile memory can include read-only memory (ROM), Programmable ROM (PROM), Electrically Programmable ROM (EPROM), Electrically Erasable Programmable ROM (EEPROM), or flash memory. Volatile memory can include Random Access Memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in a variety of forms such as Static RAM (SRAM), Dynamic RAM (DRAM), Synchronous DRAM (SDRAM), Double Data Rate SDRAM (DDRSDRAM), Enhanced SDRAM (ESDRAM), Synchronous Link DRAM (SLDRAM), Rambus Direct RAM (RDRAM), direct bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).
The technical features of the above embodiments can be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the above embodiments are not described, but should be considered as the scope of the present specification as long as there is no contradiction between the combinations of the technical features.
The above-mentioned embodiments only express several embodiments of the present application, and the description thereof is more specific and detailed, but not construed as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the concept of the present application, which falls within the scope of protection of the present application. Therefore, the protection scope of the present patent shall be subject to the appended claims.

Claims (8)

1. A day-ahead market clearing method containing wind power bidding is characterized by comprising the following steps:
pre-establishing a day-ahead market clearing model containing wind power bidding; the day-ahead market clearing model comprises a first objective function taking the minimum market clearing total cost as a target and a corresponding first constraint condition; the first mentionedAn objective function is:
Figure 682441DEST_PATH_IMAGE001
(ii) a Wherein i represents the ith node;
Figure 98379DEST_PATH_IMAGE002
representing the total cost of market clearing;
Figure 834254DEST_PATH_IMAGE003
expressing self power generation cost quotations reported by all conventional units to a system operator ISO to form row vectors;
Figure 298733DEST_PATH_IMAGE004
representing wind power probability distribution functions predicted by quotations of all wind generation sets to form vectors;
Figure 311820DEST_PATH_IMAGE005
representing the total power generation cost of the conventional unit;
Figure 774025DEST_PATH_IMAGE006
representing the expected load shedding cost generated by the system operating in a real-time market due to the uncertainty of the output power of the wind turbine generator when the day-ahead market containing wind power bidding is cleared;
Figure 856251DEST_PATH_IMAGE007
representing the power constraint upper limit of each node of the actual wind turbine generator in the system; the market clearing total cost is obtained based on a row vector formed by self power generation cost quotations submitted to a system operator ISO by all conventional units and a vector formed by a wind power probability distribution function predicted by all wind turbine quotations;
pre-establishing a real-time market clearing model containing wind power bidding; the real-time market clearing model comprises a second objective function taking the minimum load shedding expected cost as a target and a corresponding second constraint condition; the second objective function is:
Figure 62104DEST_PATH_IMAGE008
(ii) a Wherein, b is the unit cost of ISO load shedding of a system operator in a real-time market;
order to
Figure 851068DEST_PATH_IMAGE009
An optimal solution for the market clearing model is generated,
under a Vickrey-Clarke-Groves market bidding mechanism, a system operator ISO pays the cost of a conventional unit at a node i
Figure 639769DEST_PATH_IMAGE010
Comprises the following steps:
Figure 350236DEST_PATH_IMAGE011
under a Vickrey-Clarke-Groves market bidding mechanism, a system operator ISO pays the cost of the wind turbine generator at the node i
Figure 953256DEST_PATH_IMAGE012
Comprises the following steps:
Figure 65568DEST_PATH_IMAGE013
wherein the content of the first and second substances,
Figure 744942DEST_PATH_IMAGE014
representing the clear output of all conventional generator sets of the power system in the market at the day before to form a column vector;
Figure 739443DEST_PATH_IMAGE015
representing the optimal market clearing total cost when all market members participate in clearing;
Figure 287099DEST_PATH_IMAGE016
the optimal market clearing total cost of the new system is shown when the conventional unit at the node i is not included;
Figure 112973DEST_PATH_IMAGE017
representing the optimal market clearing total cost of the new system when the wind turbine set at the node i is not included;
Figure 353461DEST_PATH_IMAGE018
representing the probability distribution function of the wind power submitted by other wind turbines;
under a Vickrey-Clarke-Groves market bidding mechanism and under the condition that the conventional unit meets excitation compatibility and individual rationality conditions of the conventional unit, the clear profit is maximum in the day ahead if and only if the conventional unit declares the real power generation cost; the excitation compatibility of the conventional unit represents the real power generation cost declared by the conventional unit at a node i; the individual rationality of the conventional units represents the market clearing behavior mode of a single conventional unit only according to the self state and condition;
for the conventional set at the node i, if all the wind power sets declare the function of the probability distribution of the real wind power, the other conventional sets except the conventional set at the node i declare the power generation cost as
Figure 569679DEST_PATH_IMAGE019
In time, if the conventional unit at the node i reports the false power generation cost
Figure 796392DEST_PATH_IMAGE020
Then the payment obtained is:
Figure 617717DEST_PATH_IMAGE021
wherein the content of the first and second substances,
Figure 888162DEST_PATH_IMAGE022
shows that the conventional unit at the node i does not participate in market clearing, new systemReporting power generation cost of all conventional units of the system
Figure 794938DEST_PATH_IMAGE023
The market of the conventional unit is clear;
Figure 481134DEST_PATH_IMAGE024
representing the total output of all conventional units of the system;
Figure 297912DEST_PATH_IMAGE025
representing the false power generation cost declared by the conventional unit at the node i
Figure 880203DEST_PATH_IMAGE026
The optimal market clearing;
if the conventional unit at the node i reports the real power generation cost
Figure 133329DEST_PATH_IMAGE027
Then the payment obtained is:
Figure 826479DEST_PATH_IMAGE028
wherein the content of the first and second substances,
Figure 153555DEST_PATH_IMAGE029
representing the actual power generation cost declared by the conventional unit at the node i
Figure 779184DEST_PATH_IMAGE027
The optimal market clearing;
when the conventional unit reports the false power generation cost at the node i
Figure 660552DEST_PATH_IMAGE030
In time, the net profit is:
Figure 16447DEST_PATH_IMAGE031
wherein the content of the first and second substances,
Figure 666871DEST_PATH_IMAGE032
and
Figure 122124DEST_PATH_IMAGE033
declaring false power generation cost for conventional unit at node i respectively
Figure 366154DEST_PATH_IMAGE030
In time, the system optimizes the regular unit clearing at the scheduling node i
Figure 135527DEST_PATH_IMAGE025
At their false generation costs, respectively
Figure 30671DEST_PATH_IMAGE034
And true cost of electricity generation
Figure 860087DEST_PATH_IMAGE027
The cost of electricity generation;
Figure 512785DEST_PATH_IMAGE035
the system is cleared for all market members to participate in, and the conventional unit at the node i declares the false power generation cost as
Figure 226794DEST_PATH_IMAGE036
Reporting the power generation cost of other conventional units except the conventional unit at the node i
Figure 586231DEST_PATH_IMAGE037
The sum of the optimal power generation cost of other conventional units of the system;
when the conventional unit at the node i reports the real power generation cost
Figure 711182DEST_PATH_IMAGE027
When the temperature of the water is higher than the set temperature,the net profit is:
Figure 788859DEST_PATH_IMAGE038
Figure 306560DEST_PATH_IMAGE039
wherein the content of the first and second substances,
Figure 582820DEST_PATH_IMAGE040
representing the reported power generation cost of other conventional units except the conventional unit at the node i
Figure 19618DEST_PATH_IMAGE037
And the conventional unit at the node i reports the real power generation cost
Figure 443646DEST_PATH_IMAGE027
In time, the conventional unit at the market clearing node i is cleared
Figure 624091DEST_PATH_IMAGE041
At its true cost of electricity generation
Figure 836417DEST_PATH_IMAGE027
The cost of electricity generation;
Figure 444116DEST_PATH_IMAGE042
shows that all market members participate in clearing, and other conventional units declare the power generation cost
Figure 621019DEST_PATH_IMAGE037
And the conventional unit at the node i reports the real power generation cost
Figure 136314DEST_PATH_IMAGE027
The sum of the optimal power generation cost of other conventional units except the conventional unit at the node i;
in the above formula for calculating net profit, the first term at the right end of the two equations is irrelevant to the power generation cost declared by the conventional unit at the node i; while
Figure 590429DEST_PATH_IMAGE041
And
Figure 978816DEST_PATH_IMAGE043
declaring real power generation cost for conventional units at node i respectively
Figure 518382DEST_PATH_IMAGE027
And other conventional units report the power generation cost
Figure 165264DEST_PATH_IMAGE037
The optimal clearing plan of the conventional generator set of the system is as follows:
Figure 270623DEST_PATH_IMAGE044
therefore, the conventional unit at the node i reports the real power generation cost
Figure 220125DEST_PATH_IMAGE027
The obtained net profit is not less than the declared false power generation cost
Figure 591194DEST_PATH_IMAGE030
Net profit, so there are:
Figure 917133DEST_PATH_IMAGE045
when the conventional unit at the node i declares the real power generation cost, the real power generation cost is the excitation compatibility of the conventional unit;
the individual rationality of the conventional unit comprises the following steps: when the conventional unit at the node i reports the real power generation cost
Figure 876999DEST_PATH_IMAGE027
When the temperature of the water is higher than the set temperature,
Figure 122036DEST_PATH_IMAGE046
equivalent to the fact that all conventional units participate in market clearing, the method increases
Figure 636194DEST_PATH_IMAGE047
The constraint of (2); therefore, the optimization of the daily output model can be reduced, the objective function value is not less than that of all conventional units participating in the daily market output model, so that the method comprises the following steps:
Figure 110031DEST_PATH_IMAGE048
；
under a Vickrey-Clarke-Groves market bidding mechanism, and under the condition that the wind turbine generator meets the individual physiological conditions of excitation compatibility and the wind turbine generator, the clear profit of the wind turbine generator is the largest day ahead when and only when the wind turbine generator declares a real wind power probability distribution function; the excitation compatibility of the wind turbine generator represents a real wind power probability distribution function declared by the wind turbine generator at a node i; the individual rationality of the wind turbines represents the market clearing behavior of a single wind turbine only according to the self state and conditions;
for the wind turbine at the node i, if all conventional wind turbines declare the real power generation cost, when other wind turbines declare the function of the wind power probability distribution
Figure 127666DEST_PATH_IMAGE049
If the wind turbine at the node i declares the probability distribution of false wind power as
Figure 481287DEST_PATH_IMAGE050
Then the payment obtained is:
Figure 341796DEST_PATH_IMAGE051
wherein the content of the first and second substances,
Figure 743958DEST_PATH_IMAGE052
the function of the probability distribution of the wind power is reported based on the fact that the wind turbine generator at the node i is
Figure 488536DEST_PATH_IMAGE053
And the function of reporting the wind power probability distribution of other wind turbine generators is
Figure 216320DEST_PATH_IMAGE049
The load shedding amount of the wind turbine system is calculated;
if the wind turbine generator at the node i declares the function of the probability distribution of the real wind power as
Figure 236229DEST_PATH_IMAGE054
Then the payment obtained is:
Figure 566716DEST_PATH_IMAGE055
wherein the content of the first and second substances,
Figure 558943DEST_PATH_IMAGE056
the function for declaring the probability distribution of the real wind power based on the wind turbine at the node i is
Figure 67416DEST_PATH_IMAGE054
And the function of reporting the wind power probability distribution of other wind turbine generators is
Figure 309041DEST_PATH_IMAGE049
The load shedding amount of the wind turbine system is calculated;
for the above payment calculation formula, the first term at the right end of the equation and the wind power declared by the wind turbine at the node iFunction independent of the probability distribution; while
Figure 318585DEST_PATH_IMAGE056
Function for declaring real wind power probability distribution for wind turbine generator at node i
Figure 555532DEST_PATH_IMAGE054
And other wind turbine generators declare functions of wind power probability distribution
Figure 93960DEST_PATH_IMAGE049
The optimal load shedding of the system is determined by the time and the time
Figure 167090DEST_PATH_IMAGE057
Therefore, the following are:
Figure 511483DEST_PATH_IMAGE058
therefore, the wind turbine generator at the node i declares the function of the probability distribution of the real wind power
Figure 478302DEST_PATH_IMAGE054
The obtained payment is not less than the function of declaring the probability distribution of the false wind power
Figure 312266DEST_PATH_IMAGE050
Payment when, that is:
Figure 997325DEST_PATH_IMAGE059
the wind turbine generator at the node i reports a function of the probability distribution of the real wind power, namely the excitation compatibility of the wind turbine generator;
for the individual rationality of the wind turbine,
Figure 224039DEST_PATH_IMAGE060
equivalent to the situation that all wind turbine generators participate in market clearing, the number of the wind turbine generators is increased
Figure 842102DEST_PATH_IMAGE061
Therefore, the optimization feasible region of the model of the day-ahead output is reduced, and the objective function value is not less than that of the model of the day-ahead market output of all wind turbines, so that the method comprises the following steps:
Figure 253492DEST_PATH_IMAGE062
；
when the conventional unit declares the real power generation cost and the wind turbine unit declares the real wind power probability distribution function, the day-ahead market clearing model and the real-time market clearing model achieve the aim of minimizing the total cost of market clearing, and the method comprises the following steps: under a Vickrey-Clarke-Groves market bidding mechanism, if the conventional unit declares a real power generation cost and the wind turbine unit declares a function of a real wind power probability distribution, the market clearing total cost is minimized, and the method comprises the following steps: under the condition that the conventional unit and the wind turbine set meet excitation compatibility and individual rationality, the conventional unit declares the real power generation cost and the wind turbine set declares the function of the real wind power probability distribution, and the aim of minimizing the total cost of market clearing is automatically met by a day-ahead market clearing model.
2. The method of claim 1,
the first constraint includes:
and power balance constraint:
Figure 19322DEST_PATH_IMAGE063
constraint of the power transmission network:
Figure 908781DEST_PATH_IMAGE064
Figure 646930DEST_PATH_IMAGE065
and (3) output constraint of a conventional unit:
Figure 830219DEST_PATH_IMAGE066
output restraint of the wind turbine generator:
Figure 958712DEST_PATH_IMAGE067
wherein the content of the first and second substances,
Figure 776495DEST_PATH_IMAGE068
representing nodesiThe output power of the wind turbine generator in the market before the day; n represents a node included in the power system;
Figure 369150DEST_PATH_IMAGE069
representing the load demand power of the node i;
Figure 122343DEST_PATH_IMAGE070
representing the upper limit of the output of the conventional unit;
Figure 613498DEST_PATH_IMAGE071
representing a power generation transfer distribution factor of the node i to the transmission line l;
Figure 641497DEST_PATH_IMAGE072
representing the transmission capacity of the transmission line i.
3. The method of claim 2,
the second constraint includes:
and power balance constraint:
Figure 291921DEST_PATH_IMAGE073
constraint of the power transmission network:
Figure 75069DEST_PATH_IMAGE074
Figure 443734DEST_PATH_IMAGE075
output restraint of the wind turbine generator:
Figure 275424DEST_PATH_IMAGE076
real-time market load shedding constraint:
Figure 655720DEST_PATH_IMAGE077
wind power curtailment of the wind turbine generator:
Figure 750715DEST_PATH_IMAGE078
wherein the content of the first and second substances,
Figure 465730DEST_PATH_IMAGE079
representing the load shedding amount of the wind turbine system at the node i;
Figure 38794DEST_PATH_IMAGE080
representing the wind curtailment power of the wind turbine generator at the node i;
Figure 460548DEST_PATH_IMAGE081
representing the actual real-time power of the wind turbine at node i.
4. The utility model provides a day-ahead market goes out clear device that contains wind-powered electricity generation competitive bidding which characterized in that includes:
the first model establishing module is used for pre-establishing a day-ahead market clearing model containing wind power bidding; the day-ahead market clearing model comprises a first objective function taking the minimum market clearing total cost as a target and a corresponding first constraint condition; the first objective function is：
Figure 336232DEST_PATH_IMAGE082
(ii) a Wherein i represents the ith node;
Figure 413909DEST_PATH_IMAGE083
representing the total cost of market clearing;
Figure 852981DEST_PATH_IMAGE084
expressing self power generation cost quotations reported by all conventional units to a system operator ISO to form row vectors;
Figure 457137DEST_PATH_IMAGE004
representing wind power probability distribution functions predicted by quotations of all wind generation sets to form vectors;
Figure 628356DEST_PATH_IMAGE005
representing the total power generation cost of the conventional unit;
Figure 990067DEST_PATH_IMAGE085
representing the expected load shedding cost generated by the system operating in a real-time market due to the uncertainty of the output power of the wind turbine generator when the day-ahead market containing wind power bidding is cleared;
Figure 777370DEST_PATH_IMAGE086
representing the power constraint upper limit of each node of the actual wind turbine generator in the system; the market clearing total cost is obtained based on a row vector formed by self power generation cost quotations submitted to a system operator ISO by all conventional units and a vector formed by a wind power probability distribution function predicted by all wind turbine quotations;
the second model establishing module is used for pre-establishing a real-time market clearing model containing wind power bidding; the real-time market clearing model comprises a second objective function taking the minimum load shedding expected cost as a target and a corresponding second constraint condition; the second objective function is:
Figure 111399DEST_PATH_IMAGE087
(ii) a Wherein, b is the unit cost of ISO load shedding of a system operator in a real-time market;
order to
Figure 843732DEST_PATH_IMAGE009
An optimal solution for the market clearing model is generated,
under a Vickrey-Clarke-Groves market bidding mechanism, a system operator ISO pays the cost of a conventional unit at a node i
Figure 692739DEST_PATH_IMAGE010
Comprises the following steps:
Figure 411296DEST_PATH_IMAGE011
under a Vickrey-Clarke-Groves market bidding mechanism, a system operator ISO pays the cost of the wind turbine generator at the node i
Figure 475198DEST_PATH_IMAGE012
Comprises the following steps:
Figure 50536DEST_PATH_IMAGE013
wherein the content of the first and second substances,
Figure 590102DEST_PATH_IMAGE014
representing the clear output of all conventional generator sets of the power system in the market at the day before to form a column vector;
Figure 971405DEST_PATH_IMAGE015
representing the optimal market clearing total cost when all market members participate in clearing;
Figure 280026DEST_PATH_IMAGE088
the optimal market clearing total cost of the new system is shown when the conventional unit at the node i is not included;
Figure 26266DEST_PATH_IMAGE017
representing the optimal market clearing total cost of the new system when the wind turbine set at the node i is not included;
Figure 662914DEST_PATH_IMAGE018
representing the probability distribution function of the wind power submitted by other wind turbines;
under a Vickrey-Clarke-Groves market bidding mechanism and under the condition that the conventional unit meets excitation compatibility and individual rationality conditions of the conventional unit, the clear profit is maximum in the day ahead if and only if the conventional unit declares the real power generation cost; the excitation compatibility of the conventional unit represents the real power generation cost declared by the conventional unit at a node i; the individual rationality of the conventional units represents the market clearing behavior mode of a single conventional unit only according to the self state and condition;
for the conventional set at the node i, if all the wind power sets declare the function of the probability distribution of the real wind power, the other conventional sets except the conventional set at the node i declare the power generation cost as
Figure 723274DEST_PATH_IMAGE019
In time, if the conventional unit at the node i reports the false power generation cost
Figure 11036DEST_PATH_IMAGE089
Then the payment obtained is:
Figure 928177DEST_PATH_IMAGE090
wherein the content of the first and second substances,
Figure 442334DEST_PATH_IMAGE091
indicating that the conventional unit at the node i does not participate in marketReporting power generation cost of all conventional units of clean and new system
Figure 916172DEST_PATH_IMAGE023
The market of the conventional unit is clear;
Figure 730545DEST_PATH_IMAGE024
representing the total output of all conventional units of the system;
Figure 287428DEST_PATH_IMAGE092
representing the false power generation cost declared by the conventional unit at the node i
Figure 882357DEST_PATH_IMAGE034
The optimal market clearing;
if the conventional unit at the node i reports the real power generation cost
Figure 284520DEST_PATH_IMAGE027
Then the payment obtained is:
Figure 218978DEST_PATH_IMAGE093
wherein the content of the first and second substances,
Figure 825058DEST_PATH_IMAGE029
representing the actual power generation cost declared by the conventional unit at the node i
Figure 517071DEST_PATH_IMAGE027
The optimal market clearing;
when the conventional unit reports the false power generation cost at the node i
Figure 785241DEST_PATH_IMAGE030
In time, the net profit is:
Figure 902101DEST_PATH_IMAGE031
wherein the content of the first and second substances,
Figure 269629DEST_PATH_IMAGE032
and
Figure 855462DEST_PATH_IMAGE033
declaring false power generation cost for conventional unit at node i respectively
Figure 661744DEST_PATH_IMAGE030
In time, the system optimizes the regular unit clearing at the scheduling node i
Figure 508477DEST_PATH_IMAGE025
At their false generation costs, respectively
Figure 437119DEST_PATH_IMAGE034
And true cost of electricity generation
Figure 369303DEST_PATH_IMAGE027
The cost of electricity generation;
Figure 979276DEST_PATH_IMAGE035
the system is cleared for all market members to participate in, and the conventional unit at the node i declares the false power generation cost as
Figure 555882DEST_PATH_IMAGE036
Reporting the power generation cost of other conventional units except the conventional unit at the node i
Figure 530791DEST_PATH_IMAGE037
The sum of the optimal power generation cost of other conventional units of the system;
when the conventional unit at the node i reports the real power generation cost
Figure 12588DEST_PATH_IMAGE027
In time, the net profit is:
Figure 488569DEST_PATH_IMAGE094
Figure 44315DEST_PATH_IMAGE095
wherein the content of the first and second substances,
Figure 65492DEST_PATH_IMAGE040
representing the reported power generation cost of other conventional units except the conventional unit at the node i
Figure 769006DEST_PATH_IMAGE037
And the conventional unit at the node i reports the real power generation cost
Figure 924043DEST_PATH_IMAGE027
In time, the conventional unit at the market clearing node i is cleared
Figure 458930DEST_PATH_IMAGE041
At its true cost of electricity generation
Figure 837959DEST_PATH_IMAGE027
The cost of electricity generation;
Figure 232031DEST_PATH_IMAGE042
shows that all market members participate in clearing, and other conventional units declare the power generation cost
Figure 532038DEST_PATH_IMAGE037
And the conventional unit at the node i reports the real power generation cost
Figure 327955DEST_PATH_IMAGE027
Optimal power generation cost of other conventional units except the conventional unit at node iAnd;
in the above formula for calculating net profit, the first term at the right end of the two equations is irrelevant to the power generation cost declared by the conventional unit at the node i; while
Figure 877885DEST_PATH_IMAGE041
And
Figure 618308DEST_PATH_IMAGE043
declaring real power generation cost for conventional units at node i respectively
Figure 849570DEST_PATH_IMAGE027
And other conventional units report the power generation cost
Figure 375360DEST_PATH_IMAGE037
The optimal clearing plan of the conventional generator set of the system is as follows:
Figure 96191DEST_PATH_IMAGE096
therefore, the conventional unit at the node i reports the real power generation cost
Figure 464856DEST_PATH_IMAGE027
The obtained net profit is not less than the declared false power generation cost
Figure 358862DEST_PATH_IMAGE030
Net profit, so there are:
Figure 863793DEST_PATH_IMAGE045
when the conventional unit at the node i declares the real power generation cost, the real power generation cost is the excitation compatibility of the conventional unit;
the individual rationality of the conventional unit comprises the following steps: when the conventional unit at the node i reports the real power generation cost
Figure 755526DEST_PATH_IMAGE027
When the temperature of the water is higher than the set temperature,
Figure 955694DEST_PATH_IMAGE046
equivalent to the fact that all conventional units participate in market clearing, the method increases
Figure 794337DEST_PATH_IMAGE047
The constraint of (2); therefore, the optimization of the daily output model can be reduced, the objective function value is not less than that of all conventional units participating in the daily market output model, so that the method comprises the following steps:
Figure 216091DEST_PATH_IMAGE048
；
under a Vickrey-Clarke-Groves market bidding mechanism, and under the condition that the wind turbine generator meets the individual physiological conditions of excitation compatibility and the wind turbine generator, the clear profit of the wind turbine generator is the largest day ahead when and only when the wind turbine generator declares a real wind power probability distribution function; the excitation compatibility of the wind turbine generator represents a real wind power probability distribution function declared by the wind turbine generator at a node i; the individual rationality of the wind turbines represents the market clearing behavior of a single wind turbine only according to the self state and conditions;
for the wind turbine at the node i, if all conventional wind turbines declare the real power generation cost, when other wind turbines declare the function of the wind power probability distribution
Figure 341042DEST_PATH_IMAGE049
If the wind turbine at the node i declares the probability distribution of false wind power as
Figure 418719DEST_PATH_IMAGE050
Then the payment obtained is:
Figure 592212DEST_PATH_IMAGE097
wherein the content of the first and second substances,
Figure 947101DEST_PATH_IMAGE052
the function of the probability distribution of the wind power is reported based on the fact that the wind turbine generator at the node i is
Figure 118319DEST_PATH_IMAGE098
And the function of reporting the wind power probability distribution of other wind turbine generators is
Figure 807926DEST_PATH_IMAGE049
The load shedding amount of the wind turbine system is calculated;
if the wind turbine generator at the node i declares the function of the probability distribution of the real wind power as
Figure 785110DEST_PATH_IMAGE054
Then the payment obtained is:
Figure 119139DEST_PATH_IMAGE055
wherein the content of the first and second substances,
Figure 429353DEST_PATH_IMAGE099
the function for declaring the probability distribution of the real wind power based on the wind turbine at the node i is
Figure 216043DEST_PATH_IMAGE054
And the function of reporting the wind power probability distribution of other wind turbine generators is
Figure 996917DEST_PATH_IMAGE049
The load shedding amount of the wind turbine system is calculated;
for the above payment calculation formula, the first term at the right end of the equation and the wind power declared by the wind turbine at the node iFunction independent of power probability distribution; while
Figure 575666DEST_PATH_IMAGE099
Function for declaring real wind power probability distribution for wind turbine generator at node i
Figure 88687DEST_PATH_IMAGE054
And other wind turbine generators declare functions of wind power probability distribution
Figure 159411DEST_PATH_IMAGE049
The optimal load shedding of the system is determined by the time and the time
Figure 557026DEST_PATH_IMAGE057
Therefore, the following are:
Figure 600068DEST_PATH_IMAGE100
therefore, the wind turbine generator at the node i declares the function of the probability distribution of the real wind power
Figure 674203DEST_PATH_IMAGE054
The obtained payment is not less than the function of declaring the probability distribution of the false wind power
Figure 497803DEST_PATH_IMAGE050
Payment when, that is:
Figure 558163DEST_PATH_IMAGE059
the wind turbine generator at the node i reports a function of the probability distribution of the real wind power, namely the excitation compatibility of the wind turbine generator;
for the individual rationality of the wind turbine,
Figure 596657DEST_PATH_IMAGE060
equivalent to the situation that all wind turbine generators participate in market clearing, the number of the wind turbine generators is increased
Figure 513797DEST_PATH_IMAGE061
Therefore, the optimization feasible region of the model of the day-ahead output is reduced, and the objective function value is not less than that of the model of the day-ahead market output of all wind turbines, so that the method comprises the following steps:
Figure 762376DEST_PATH_IMAGE062
；
the clearing module is used for enabling the day-ahead market clearing model and the real-time market clearing model to achieve the aim of minimizing the total clearing cost of the market when the real power generation cost is declared by a conventional unit and the real wind power probability distribution function is declared by a wind turbine unit; the method is specifically used for: under a Vickrey-Clarke-Groves market bidding mechanism, if the conventional unit declares a real power generation cost and the wind turbine unit declares a function of a real wind power probability distribution, the market clearing total cost is minimized, and the method comprises the following steps: under the condition that the conventional unit and the wind turbine set meet excitation compatibility and individual rationality, the conventional unit declares the real power generation cost and the wind turbine set declares the function of the real wind power probability distribution, and the aim of minimizing the total cost of market clearing is automatically met by a day-ahead market clearing model.
5. The apparatus of claim 4, wherein the first constraint comprises:
and power balance constraint:
Figure 485482DEST_PATH_IMAGE063
constraint of the power transmission network:
Figure 768695DEST_PATH_IMAGE064
Figure 856737DEST_PATH_IMAGE065
and (3) output constraint of a conventional unit:
Figure 467978DEST_PATH_IMAGE066
output restraint of the wind turbine generator:
Figure 870141DEST_PATH_IMAGE101
wherein the content of the first and second substances,
Figure 804599DEST_PATH_IMAGE102
representing nodesiThe output power of the wind turbine generator in the market before the day; n represents a node included in the power system;
Figure 391438DEST_PATH_IMAGE069
representing the load demand power of the node i;
Figure 349029DEST_PATH_IMAGE070
representing the upper limit of the output of the conventional unit;
Figure 692899DEST_PATH_IMAGE103
representing a power generation transfer distribution factor of the node i to the transmission line l;
Figure 481863DEST_PATH_IMAGE072
representing the transmission capacity of the transmission line i.
6. The apparatus of claim 5, wherein the second constraint comprises:
and power balance constraint:
Figure 114970DEST_PATH_IMAGE104
constraint of the power transmission network:
Figure 684491DEST_PATH_IMAGE074
Figure 694036DEST_PATH_IMAGE075
output restraint of the wind turbine generator:
Figure 337506DEST_PATH_IMAGE076
real-time market load shedding constraint:
Figure 16881DEST_PATH_IMAGE077
wind power curtailment of the wind turbine generator:
Figure 949065DEST_PATH_IMAGE078
wherein the content of the first and second substances,
Figure 559037DEST_PATH_IMAGE105
representing the load shedding amount of the wind turbine system at the node i;
Figure 384911DEST_PATH_IMAGE080
representing the wind curtailment power of the wind turbine generator at the node i;
Figure 359820DEST_PATH_IMAGE081
representing the actual real-time power of the wind turbine at node i.
7. A computer device comprising a memory and a processor, the memory storing a computer program, wherein the processor when executing the computer program implements the steps of the method of any one of claims 1 to 3.
8. A computer-readable storage medium, on which a computer program is stored which, when being executed by a processor, carries out the steps of the method for day-ahead market clearing according to any one of claims 1 to 3.
CN201911022977.2A 2019-10-25 2019-10-25 Day-ahead market clearing method and device containing wind power bidding, and computer equipment Active CN110738377B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911022977.2A CN110738377B (en) 2019-10-25 2019-10-25 Day-ahead market clearing method and device containing wind power bidding, and computer equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911022977.2A CN110738377B (en) 2019-10-25 2019-10-25 Day-ahead market clearing method and device containing wind power bidding, and computer equipment

Publications (2)

Publication Number Publication Date
CN110738377A CN110738377A (en) 2020-01-31
CN110738377B true CN110738377B (en) 2021-11-30

Family

ID=69271469

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911022977.2A Active CN110738377B (en) 2019-10-25 2019-10-25 Day-ahead market clearing method and device containing wind power bidding, and computer equipment

Country Status (1)

Country Link
CN (1) CN110738377B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111860994B (en) * 2020-07-14 2023-12-08 北京清能互联科技有限公司 Spot market transaction method and device considering that two sides participate in peak shaving
CN113112296B (en) * 2021-04-09 2022-07-01 国电南瑞南京控制系统有限公司 Typical scene library construction method and device for centralized spot market
CN115630527B (en) * 2022-11-08 2023-11-03 中国南方电网有限责任公司 Exception model detection method, device, computer equipment and storage medium

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6088689A (en) * 1995-11-29 2000-07-11 Hynomics Corporation Multiple-agent hybrid control architecture for intelligent real-time control of distributed nonlinear processes
CN108879796A (en) * 2018-08-10 2018-11-23 广东电网有限责任公司 Electric power ahead market goes out clear calculation method, system, device and readable storage medium storing program for executing
CN110021965A (en) * 2019-03-20 2019-07-16 南方电网科学研究院有限责任公司 Method for formulating wind power bidding strategy based on prediction error
CN110224398A (en) * 2019-06-12 2019-09-10 广东电网有限责任公司 Power spot market goes out clearing method, device, system and storage medium to one kind a few days ago

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6088689A (en) * 1995-11-29 2000-07-11 Hynomics Corporation Multiple-agent hybrid control architecture for intelligent real-time control of distributed nonlinear processes
CN108879796A (en) * 2018-08-10 2018-11-23 广东电网有限责任公司 Electric power ahead market goes out clear calculation method, system, device and readable storage medium storing program for executing
CN110021965A (en) * 2019-03-20 2019-07-16 南方电网科学研究院有限责任公司 Method for formulating wind power bidding strategy based on prediction error
CN110224398A (en) * 2019-06-12 2019-09-10 广东电网有限责任公司 Power spot market goes out clearing method, device, system and storage medium to one kind a few days ago

Also Published As

Publication number Publication date
CN110738377A (en) 2020-01-31

Similar Documents

Publication Publication Date Title
CN110738377B (en) Day-ahead market clearing method and device containing wind power bidding, and computer equipment
Moiseeva et al. Exercise of market power on ramp rate in wind-integrated power systems
Gilli et al. A data-driven optimization heuristic for downside risk minimization
Vuelvas et al. Limiting gaming opportunities on incentive-based demand response programs
CN103729698B (en) Requirement responding scheduling method for wind power uncertainty
Steriotis et al. Strategic and network-aware bidding policy for electric utilities through the optimal orchestration of a virtual and heterogeneous flexibility assets’ portfolio
Rocha et al. Impact of CO2 cap-and-trade programs on restructured power markets with generation capacity investments
Downward et al. Electricity retail contracting under risk-aversion
Aghaei et al. Uncertainty management in multiobjective hydro-thermal self-scheduling under emission considerations
Gil et al. Modeling competition in transmission expansion
Dimitriadis et al. Optimization-based economic analysis of energy storage technologies in a coupled electricity and natural gas market
Johnathon et al. A proposed hedge-based energy market model to manage renewable intermittency
Spiridonova Transmission capacities and competition in Western European electricity market
Neame et al. Offer stack optimization in electricity pool markets
CN109063859B (en) Power grid equipment maintenance optimization processing method and device
Haidar et al. A market framework for energy bidding decision-making strategy to provide a competitive mechanism in the context of deregulated electricity market
Chen et al. Chance-constrained multi-stage stochastic energy system expansion planning with demand satisfaction flexibility
CN110827067A (en) Day-ahead market clearing method and device with demand side bidding, and computer equipment
CN113177655A (en) Comprehensive energy system multi-main-body operation optimization method and device based on reinforcement learning
CN111754076A (en) Method and equipment suitable for electric power wholesale market mode evaluation
CN110288274A (en) Generation schedule preparation method, device, equipment and the storage medium of multicomputer
Li et al. Optimal wind power capacity decision consider commitment contracts under uncertain power supply and electricity demand in China
CN112072644A (en) Thermal power plant day-ahead market yield decision-making operation plan optimization method and system
Chamba et al. Application of hybrid heuristic optimization algorithms for solving optimal regional dispatch of energy and reserve considering the social welfare of the participating markets
Ciornei et al. Efficient hybrid optimization solution for the economic dispatch with nonsmooth cost function

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
GR01 Patent grant
GR01 Patent grant