CN111160640A - Method for clearing current market of hydropower station water consumption rate iterative correction electric power system before day - Google Patents

Method for clearing current market of hydropower station water consumption rate iterative correction electric power system before day Download PDF

Info

Publication number
CN111160640A
CN111160640A CN201911338730.1A CN201911338730A CN111160640A CN 111160640 A CN111160640 A CN 111160640A CN 201911338730 A CN201911338730 A CN 201911338730A CN 111160640 A CN111160640 A CN 111160640A
Authority
CN
China
Prior art keywords
power
unit
hydropower
water consumption
consumption rate
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.)
Granted
Application number
CN201911338730.1A
Other languages
Chinese (zh)
Other versions
CN111160640B (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.)
Dalian University of Technology
Original Assignee
Dalian University of Technology
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 Dalian University of Technology filed Critical Dalian University of Technology
Priority to CN201911338730.1A priority Critical patent/CN111160640B/en
Publication of CN111160640A publication Critical patent/CN111160640A/en
Application granted granted Critical
Publication of CN111160640B publication Critical patent/CN111160640B/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
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06315—Needs-based resource requirements planning or analysis
    • 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/0601—Electronic shopping [e-shopping]
    • G06Q30/0611—Request for offers or quotes
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
    • 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)
  • Human Resources & Organizations (AREA)
  • Development Economics (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Marketing (AREA)
  • Finance (AREA)
  • Accounting & Taxation (AREA)
  • Theoretical Computer Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Business, Economics & Management (AREA)
  • Game Theory and Decision Science (AREA)
  • Tourism & Hospitality (AREA)
  • Quality & Reliability (AREA)
  • Operations Research (AREA)
  • Health & Medical Sciences (AREA)
  • Educational Administration (AREA)
  • Data Mining & Analysis (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 invention belongs to the field of spot markets with participation of large-scale hydropower, and relates to a method for clearing the spot market of a hydropower station in the day ahead of an electric power system by iteratively correcting the water consumption rate. The water consumption rate is used for replacing complex NHQ curve power generation amount optimization calculation, the cascade power station space-time hydraulic power and power complex relation is coupled into an SCUC/SCED model, the water consumption rate is iteratively corrected through hydropower fine calculation and is converged within the accuracy of 0.01; the complete day-ahead spot shipment clearing solution method which comprises the steps of beforehand information management, declaration, SCUC/SCED and hydropower fine calculation is constructed. The method effectively solves the problem that the generated energy of the power stations on the upstream and the downstream of the cascade is not matched when large-scale hydropower participates in the daily market clearing, and simultaneously avoids the solving difficulty brought by complex constraint of hydropower to the clear optimization model. The method has important significance for the promotion and construction of the electric power spot market in the hydropower high-proportion area and the cascade power station with tight space-time coupling and the upstream and downstream power stations belonging to different benefit subjects.

Description

Method for clearing current market of hydropower station water consumption rate iterative correction electric power system before day
Technical Field
The invention belongs to the field of large-scale hydropower-participated power system spot market in the day ahead, relates to a hydropower station water consumption rate iterative correction power system spot market clearing method, and is a novel method for solving the problem that the upstream and downstream generated energy of a cascade power station is not matched when the spot market participated by a multi-operation main cascade power station is cleared.
Background
The method is characterized in that the clear nature of the existing market of the power system at present is a safety constraint unit combination problem of electric energy clear and reserve, mature clear algorithms exist, for the power system containing water and electricity, particularly the power system with high water-electricity ratio and large water-electricity scale, no mature clear algorithms exist at home and abroad so far, the situation is the situation of Yunnan and Sichuan power grids in China, until 2018, 6666 ten thousand kW and 2699 hundred million kWh of the Yunnan hydropower installation respectively account for 71.52% and 83.21% of the whole province, 7823.9 ten thousand kW and 2982.2 kWh of generated power are respectively 79.57% and 79.30% of the generated power at the end of the generation, the existing market at present day of the Yunnan and Sichuan power grids in which the water and electricity are dominant, the method is greatly different from the existing power market at home and other power sources mainly in the world, ①, the tight coupling of the hydropower station and the tight coupling of the hydropower station, the tight hydraulic power station and the tight coupling between the upstream and downstream ladder power grids in which the same thermal power supply is mainly used, the same, the linear and the linear coordination relationship between the same and the linear calculation of the straight line, the adjustment of the straight line, the adjustment of the straight line, the adjustment of the straight line, the straight line of the adjustment of the straight line, the adjustment of the straight line of the adjustment of the straight line, the straight line.
Aiming at the problems, the invention provides a method for clearing the current spot of an electric power system day ahead based on iterative correction of the water consumption rate of a hydropower station by depending on a major plan support project (91547201) of the national science fund, the method comprises the steps of taking the hydropower station as a bidding unit and taking a thermal power unit as a bidding unit to participate in SCUC/SCED calculation, replacing complex NHQ curve power generation amount optimization calculation with the water consumption rate, coupling the cascade power station space-time hydraulic power and power complex relation into an SCUC/SCED model, and iteratively correcting the water consumption rate by combining the clearing result of the power station through fine calculation of hydropower, so that the problem that the power generation amounts of upstream and downstream power stations of the cascade power station are not matched when the current market is cleared is effectively solved, and meanwhile, the solving difficulty brought by the clearing optimization model due to nonlinear constraint is avoided.
Disclosure of Invention
The invention aims to provide an efficient and practical technical method for optimizing clearing of the existing electric power market containing the step hydropower, fully considers the space-time coupling relation of the hydropower of the step hydropower and provides an effective technical path for the implementation of the existing electric power market containing the step hydropower.
The technical scheme of the invention is as follows:
a method for clearing a spot market of a hydropower station before day by iteratively correcting water consumption rate of a hydropower station comprises the following steps of beforehand information management, declaration, an SCUC/SCED clearing model and hydropower refined calculation:
(1) determining SCUC/SCED clear model constraint boundary conditions
Step1.1: extracting the prediction information of the power utilization side in the day from the load prediction system, wherein the prediction information comprises the load prediction of a whole network system, the bus load prediction, the solar output process of wind power and photovoltaic power, reliability information and standby constraint information;
step1.2: the grid-connected generator set provides power generation side information which comprises minimum continuous start-up and shut-down time of the generator set, minimum/maximum technical output of the generator set, up-down climbing rate of the generator set, start-up and shut-down state of the generator set one day before, maintenance plan of the generator set and state constraint of the generator set; the hydroelectric generating set needs to provide hydraulic parameters and the flow of leaving the reservoir in the previous day;
step1.3: the power dispatching mechanism provides power grid technical parameters including node pair section sensitivity information, power station pair section sensitivity information, power generation node and power station corresponding relation, section limit and direction, power grid topological model and unit start-stop plan.
(2) The power generation enterprise market declares according to the unit of competing for price, wherein thermal power uses the unit as the unit of bidding for price, and water and electricity uses the power station as the unit of competing for price, and specific default declares information contains: the starting cost of the unit is one incremental power generation amount price curve in 96 time periods throughout the day, and the maximum power generation amount price curve does not exceed five sections; the generated electricity price curve is composed of declared energy price (yuan/MWh) and corresponding generated electricity capacity section (MW).
(3) The method is characterized in that the lowest electricity purchasing cost of the power grid side is used as a target function of the SCUC/SCED output clear optimization, a mixed integer programming model is adopted in the SCUC/SCED model, the water consumption rate is used for replacing the complex NHQ curve generated energy optimizing calculation, and under the condition that the water consumption rate is assumed to be fixed, the hydraulic power and electric power space-time coupling constraints of the cascade power station are nested into the output clear model, and the method specifically comprises the following steps:
step3.1: hydraulic parameters reported according to a hydroelectric bidding unit, comprisingThe current reservoir capacity and the average water consumption rate of the previous day, the water consumption rate of the coming date is assumed to be the average water consumption rate r of the previous dayrThe output p of the hydropower bidding unitr(t) is represented by pr(t)＝ur(t)/rr*3.6，ur(t) is the average warehousing flow of the r bidding units in the t time period;
step3.2: the hydraulic coupling constraint related to the r hydropower bidding unit is coupled to the outbound model and comprises water balance constraint, power generation flow upper and lower limit constraint, discharge upper and lower limit constraint, reservoir capacity constraint, water abandoning constraint and initial reservoir capacity constraint;
step3.3: for different step watershed power stations, steps Step3.1 and Step3.2 need to be repeated to nest the step hydraulic constraints into the supernatant model for the supernatant calculation.
(4) And (3) refining the clear result of the hydropower bidding unit by using an electricity water-ordering method to update the daily average water consumption rate, which comprises the following specific steps:
step4.1: determining the water consumption rate r of each hydropower bidding unit in each period by adopting a method of determining water consumption by electricityr(t) and calculating the average daily water consumption
Figure BDA0002331674810000021
Wherein T is a set of time periods;
step4.2: updating the water consumption rate by adopting a dichotomy, and reselecting the water consumption rate for the r hydropower bidding unit
Figure BDA0002331674810000031
Performing clearing calculation;
step4.3: if it is
Figure BDA0002331674810000032
Stopping the iteration water consumption rate and outputting a clear result.
The invention has the following beneficial effects: according to the method, the water consumption rate of the hydropower station is used for replacing the complex NHQ curve power generation amount optimization calculation, the water consumption rate is used for coupling the hydraulic power and power space-time coupling relation of the cascade power station to the clearing model, the problem that the power generation amounts of the upstream power station and the downstream power station of the cascade power station are not matched when the power spot market is cleared is effectively solved, and an efficient solving method is provided for the power spot clearing containing the large-scale cascade water power space-time coupling problem.
Drawings
FIG. 1 is a flow chart of an electric power spot inventory solution based on iterative correction of hydropower station water consumption rates;
FIG. 2 is a graphical illustration of a bid amount curve declared by a bidding unit;
fig. 3 is a diagram of a full network load balancing process.
Detailed Description
The following further describes a specific embodiment of the present invention with reference to the drawings and technical solutions.
The energy distribution of China has strong regional characteristics, the water and electricity resources in the southwest region are rich, and the high-water-electricity-ratio regions have certain particularity when the construction of the spot market is carried out due to the unique natural characteristics of the regions. The cascade power station has close time-space hydraulic power connection, the discharge of an upstream power station directly influences the power generation condition of a downstream power station, and meanwhile, the characteristics of time lag constraint, uncertainty of water and power output and the like exist between the upstream power station and the downstream power station, so that great challenges are brought to the establishment and solution of a export model of a current stock market of a power system. If the hydropower and power connection of the cascade power station is not considered in the clearing model, the problem that the upstream and downstream power generation quantities of the cascade power station are not matched is caused, namely the upstream power station has too high quotation, the water drainage is less, the downstream power station has too low quotation, but the medium-winning power quantity cannot be finished due to poor self-regulation performance; it is also possible that the upstream quote is low, but the bid amount cannot be met due to less water coming from itself. Based on the problems, the invention provides a method for clearing the spot market of the power system before the day based on iterative correction of the water consumption rate of the hydropower station.
And selecting the lowest power purchase cost of the power grid side as an objective function of SCUC/SCED calculation.
Figure BDA0002331674810000033
In the formula:
Figure BDA0002331674810000034
operating costs for bidding units;
Figure BDA0002331674810000035
a boot cost for a bid unit;
Figure BDA0002331674810000036
a shutdown cost for a bid unit; t and j are respectively the indexes of the time interval and the bidding unit; t, J are collections of time periods and bid units, respectively. The terms in the objective function are defined as follows:
bid unit operating cost
Figure BDA0002331674810000037
The definition is as follows:
Figure BDA0002331674810000038
Figure BDA0002331674810000041
δ1(j,t)≤T1,j-P j,
Figure BDA0002331674810000042
δl(j,t)≤Tl,j-Tl-1,j,
Figure BDA0002331674810000043
Figure BDA0002331674810000044
δl(j,t)≥0,
Figure BDA0002331674810000045
Aj＝F1,j×p j
in the formula: p is a radical ofj(t)Representing the output of the j bidding unit at the time t; v. ofj(t) is a binary variable, the j bidding unit is started to be 1 at the time of t and is shut down to be 0; fl,jBidding for the l section of the j bidding unit; deltal(j, t) represents the I segment output of the j bidding unit at the time t; NLjThe number of quote segments representing j bid units;P jrepresenting the minimum output of the j bidding unit;p jand the starting point of the contribution interval of the j bidding unit quoted in the 1 st segment is shown.
The starting and shutdown costs of the bidding unit
Figure BDA0002331674810000046
The self characteristics of hydropower are considered, and the starting and shutdown costs are generally set to be 0; for a bidding unit of the thermal power generating unit:
Figure BDA0002331674810000047
Figure BDA0002331674810000048
Figure BDA0002331674810000049
Figure BDA00023316748100000410
in the formula:
Figure BDA00023316748100000411
representing the cost of the kth interval of the j stepped starting cost function of the thermal power generating unit; cjRepresenting the shutdown cost of the thermal power generating unit j; NDjRepresenting the number of the starting cost segments of the thermal power generating unit j; and R represents a hydroelectric generating set.
The above objective function needs to satisfy the following constraints:
and (4) system constraint:
1) power balance constraint
Figure BDA00023316748100000412
In the formula: d (t) represents the total load at time t.
2) System backup constraints
Figure BDA00023316748100000413
In the formula:
Figure BDA00023316748100000414
representing the possible maximum output of the j unit at the time t, α is the load reserve rate.
3) Safety restraint
Figure BDA0002331674810000051
Figure BDA0002331674810000052
In the formula: pfl,tThe active power of the line l passing through at the moment t;P fl、
Figure BDA0002331674810000053
upper and lower limits for the active transmission allowed by line l; di(t) is the load power of node i; sfl-j、Sfl-iThe transfer factors of the bidding unit j and the node load i to the line l are constants related to the line impedance parameter.
The bidding unit constraint:
1) upper and lower bound constraints for bid unit operation and standby
Figure BDA0002331674810000054
Figure BDA0002331674810000055
In the formula:
Figure BDA0002331674810000056
an upper limit is imposed on bid unit j.
2) Slope rate constraint of up-down climbing of bidding unit
Figure BDA0002331674810000057
Figure BDA0002331674810000058
Figure BDA0002331674810000059
In the formula: RU (RU)jThe rate of ascent for j bid units; RDjA downward ramp rate for j bid units; SUjThe starting ramp rate for the j bid unit; SDjThe shutdown ramp rate for the j bid unit.
3) Minimum continuous boot time constraint of bidding unit
Figure BDA00023316748100000510
Figure BDA00023316748100000511
Figure BDA00023316748100000512
Figure BDA0002331674810000061
Figure BDA0002331674810000062
Figure BDA0002331674810000063
In the formula: UT (unified device)jMinimum continuous run time for j bidding units;
Figure BDA0002331674810000064
bid for j the number that the unit has run before participating in the calculated 1 st period.
4) Minimum continuous off time constraint of bidding unit
Figure BDA0002331674810000065
Figure BDA0002331674810000066
Figure BDA0002331674810000067
Figure BDA0002331674810000068
Figure BDA0002331674810000069
Figure BDA00023316748100000610
In the formula: DTjMinimum downtime to bid on units for j;
Figure BDA00023316748100000611
the bid for j unit has been off for the number of periods before participating in the calculated 1 st period.
Hydropower bidding unit specific constraints:
1) water balance constraint
vrr(t+1)＝vrr(t)+3600×(Ir(t)-ur(t))Δt
Figure BDA00023316748100000612
In the formula: vr isr(t)，vrr(t +1) represents the storage capacity of the reservoir r at time t and time t +1, respectively, Ir(t) represents the warehousing flow rate at time t, ur(t) represents the delivery flow at time t, including the generated flow qr-1(t-T) and reject flow Rate
Figure BDA00023316748100000613
Tau represents the water flow time lag between the power station (r-1) and the power station r.
2) Restraint of output
qr(t)＝pr(t)*rr/3.6
In the formula: r isrRepresenting the average water consumption coefficient of the plant r.
3) Upper and lower limit constraints of generated flow
Figure BDA00023316748100000614
In the formula:
Figure BDA0002331674810000071
Q rthe upper limit and the lower limit of the generated flow are provided.
4) Waste water flow restriction
Figure BDA0002331674810000072
5) Outbound flow constraint
Figure BDA0002331674810000073
In the formula:
Figure BDA0002331674810000074
U rthe upper limit and the lower limit of the ex-warehouse flow.
6) Capacity constraint
Figure BDA0002331674810000075
In the formula:
Figure BDA0002331674810000076
represents the upper and lower limits of the storage capacity.
7) Initial storage capacity constraint
Figure BDA0002331674810000077
In the formula:
Figure BDA0002331674810000078
a capacity constraint is initiated for the scheduling period.
Based on the method, the invention is completely applied once, and the specific implementation scheme is realized according to the following steps (1) to (4), and the flow chart of the specific method is shown in figure 1:
(1) determining SCUC/SCED clear model constraint boundary conditions, wherein the detailed steps are as follows:
step1.1, extracting daily power utilization side prediction information from a load prediction system, wherein the daily power utilization side prediction information comprises the load prediction of a whole network system, the bus load prediction, the daily output process of wind power and photovoltaic power, reliability information and standby constraint information;
step1.2, providing power generation side information by the grid-connected generator set, wherein the power generation side information comprises minimum continuous start-up and shut-down time of the generator set, minimum/maximum technical output of the generator set, up-down climbing rate of the generator set, start-up and shut-down states of the generator set one day before, maintenance plans of the generator set and state constraints of the generator set; the hydroelectric generating set needs to provide hydraulic parameters and the flow of leaving the reservoir in the previous day;
and Step1.3, the power dispatching mechanism provides power grid technical parameters including node pair section sensitivity information, power station section sensitivity information, power generation node and power station corresponding relation, section limit and direction, a power grid topological model and a unit stop plan.
(2) The power generation enterprise market declares according to the unit of competing for price, wherein thermal power uses the unit as the unit of bidding for price, and water and electricity uses the power station as the unit of competing for price, and specific default declares information contains: the starting cost of the unit and a curve of increasing power generation amount price in a 96-period whole day are shown in figure 2, and the maximum is not more than five sections. The generated electricity price curve is composed of declared energy price (yuan/MWh) and corresponding generated electricity capacity section (MW).
(3) The method is characterized in that the lowest electricity purchasing cost of the power grid side is used as a target function of the SCUC/SCED output clear optimization, a mixed integer programming model is adopted in the SCUC/SCED model, the water consumption rate is used for replacing the complex NHQ curve generated energy optimizing calculation, the hydraulic power and electric power space-time coupling constraint of the cascade power station is nested into the output clear model, and the method specifically comprises the following steps:
step3.1: according to the hydraulic parameters reported by the hydropower bidding unit, the current reservoir capacity and the average water consumption rate of the previous day are included, and the water consumption rate of the Qing day is assumed to be the average water consumption rate r of the previous dayrThe output p of the hydropower bidding unitr(t) may be represented by pr(t)＝ur(t)/rr*3.6，ur(t) is the average warehousing flow of the r bidding units in the t time period;
step3.2: hydraulic coupling constraints related to the r hydropower bidding unit are also coupled into the clear model and respectively comprise water balance constraints, power generation flow upper and lower limit constraints, discharge upper and lower limit constraints, reservoir capacity constraints, water abandoning constraints and initial reservoir capacity constraints;
step3.3: for different Step watershed power stations, steps 1 and 2 are repeated to nest the Step hydraulic constraints into the supernatant model for the supernatant calculation.
(4) And (3) refining the clear result of the hydropower bidding unit by using an electricity water-ordering method to update the daily average water consumption rate, which comprises the following specific steps:
step4.1: determining the water consumption rate r of each hydropower bidding unit in each period by adopting a method of determining water consumption by electricityr(t) and calculating the average daily water consumption
Figure BDA0002331674810000081
Step4.2: updating the water consumption rate by adopting a dichotomy, and reselecting the water consumption rate for the r hydropower bidding unit
Figure BDA0002331674810000082
Performing clearing calculation;
step4.3: if it is
Figure BDA0002331674810000083
Stopping the iteration water consumption rate and outputting a clear result.
The method is used for simulating that the power system in the south West province is cleared in stock before a certain day before a flood, the water and electricity saving resources are rich, the installation scale of the water and electricity exceeds 6600 ten thousand kW by 2018, the annual generated energy of the water and electricity is close to 2700 hundred million kWh, the province is a typical power grid in China regardless of installation or generated energy, and the installation scale of the system is large; a plurality of large cascade watersheds are positioned in the province, the hydraulic power connection and the electric power connection between cascade power stations are closely coupled in time and space, the hydraulic power connection and the electric power connection exist between a time period and power stations on upstream and downstream, the same power station is merged into different connecting lines, and the power stations on the upstream and downstream of the same watersheds and the different watersheds are merged into the complex electric power connection of the same connecting lines; the water-saving power stations have a large number and different adjusting performances, the functions of the water-saving power stations in an electric power system and market competition are different, how to simplify and optimize the water-saving power stations is realized, and the efficiency of the shipment of the spot market is directly related; meanwhile, the upstream and downstream power stations of the provincial part of the watershed have complex interest relationship, and mutual-benefit bidding strategies are difficult to form in the market environment. Therefore, the province practical situation is used as an engineering background, and the method provided by the patent is applied to check the effectiveness and the feasibility of the implementation. When the spot shipment is cleared before the province is simulated, the thermal power unit serves as a bidding unit, the hydropower station serves as a bidding unit, each bidding unit declares a pricing relation in a form shown in fig. 2, the wind power and photovoltaic power output is equivalent and simplified into direct load deduction, and the grid-frame topological structure of the power grid is calculated according to the actual condition of the province. In order to verify the feasibility of the method, the cascade hydraulic constraint under the fixed water consumption rate is added to the cascade power station in the province A watershed, and the daily average water consumption rate is iteratively corrected through the fine calculation of the water and electricity, namely, the method disclosed by the patent ensures the balance of electric power and electric quantity in each time period as shown in figure 3; meanwhile, the daily average water consumption rate of each power station in the basin A gradually becomes stable through finite iterations, and as shown in table 1, the balance matching relationship of the winning power quantities between the upstream power station and the downstream power station of the cascade power station is maintained, so that the winning power quantities in the cascade power station have performability, and the reliability of a calculation result is ensured. The method effectively solves the problem that the generated energy of the power stations at the upstream and downstream of the cascade power station is not matched when the power spot market is cleared, and provides an efficient and feasible solving way for clearing the power spot in the area with high water-electricity ratio.
TABLE 1 convergence of average water consumption rate in basin A day
Figure BDA0002331674810000091

Claims (2)

1. A method for solving the existing market clearance of a hydropower station in the day front by iterative correction of water consumption rate of a hydropower station is characterized by comprising the following steps:
(1) determining SCUC/SCED clear model constraint boundary conditions
Step1.1, extracting daily power utilization side prediction information from a load prediction system, wherein the daily power utilization side prediction information comprises the load prediction of a whole network system, the bus load prediction, the daily output process of wind power and photovoltaic power, reliability information and standby constraint information;
step1.2, providing power generation side information by the grid-connected generator set, wherein the power generation side information comprises minimum continuous start-up and shut-down time of the generator set, minimum/maximum technical output of the generator set, up-down climbing rate of the generator set, start-up and shut-down states of the generator set one day before, maintenance plans of the generator set and state constraints of the generator set; the hydroelectric generating set needs to provide hydraulic parameters and the flow of leaving the reservoir in the previous day;
step1.3, the power dispatching mechanism provides power grid technical parameters, including node pair section sensitivity information, power station pair section sensitivity information, power generation node and power station corresponding relation, section limit and direction, power grid topology model and unit must start and stop plan;
(2) the power generation enterprise market declares according to the unit of competing for price, wherein thermal power uses the unit as the unit of bidding for price, and water and electricity uses the power station as the unit of competing for price, and specific default declares information contains: the starting cost of the unit is one incremental power generation amount price curve in 96 time periods throughout the day, and the maximum power generation amount price curve does not exceed five sections; the power generation volume price curve is composed of declared energy price and corresponding power generation capacity section;
(3) the method is characterized in that the lowest electricity purchasing cost of the power grid side is used as a target function of the SCUC/SCED output clear optimization, a mixed integer programming model is adopted in the SCUC/SCED model, the water consumption rate is used for replacing the complex NHQ curve generated energy optimizing calculation, under the condition that the water consumption rate is assumed to be fixed, the hydraulic power space-time coupling constraint of the cascade power station is nested into the output clear model, and the method specifically comprises the following steps:
step3.1: according to the hydraulic parameters reported by the hydropower bidding unit, the current reservoir capacity and the average water consumption rate of the previous day are included, and the water consumption rate of the Qing day is assumed to be the average water consumption rate r of the previous dayrThe output p of the hydropower bidding unitr(t) is represented by pr(t)＝ur(t)/rr*3.6，ur(t) is the average warehousing flow of the r bidding units in the t time period;
step3.2: the hydraulic coupling constraint related to the r hydropower bidding unit is coupled to the outbound model and comprises water balance constraint, power generation flow upper and lower limit constraint, discharge upper and lower limit constraint, reservoir capacity constraint, water abandoning constraint and initial reservoir capacity constraint;
step3.3: for different step watershed power stations, the steps Step3.1 and Step3.2 are required to be repeated to nest the step hydraulic power constraint into a supernatant model for calculating the supernatant;
(4) and (3) refining the clear result of the hydropower bidding unit by using an electricity water-ordering method to update the daily average water consumption rate, which comprises the following specific steps:
step4.1: determining the water consumption rate r of each hydropower bidding unit in each period by adopting a method of determining water consumption by electricityr(t) and calculating the average daily water consumption
Figure FDA0002331674800000011
Wherein T is a set of time periods;
step4.2: updating the water consumption rate by adopting a dichotomy, and reselecting the water consumption rate for the r hydropower bidding unit
Figure FDA0002331674800000012
Performing clearing calculation;
step4.3: if it is
Figure FDA0002331674800000021
Stopping the iteration water consumption rate and outputting a clear result.
2. The method according to claim 1, wherein the objective function is:
Figure FDA0002331674800000022
in the formula:
Figure FDA0002331674800000023
operating costs for bidding units;
Figure FDA0002331674800000024
a boot cost for a bid unit;
Figure FDA0002331674800000025
a shutdown cost for a bid unit; t and j are respectively the indexes of the time interval and the bidding unit; t, J are sets of periods and bid units, respectively; the terms in the objective function are defined as follows:
bid unit operating cost
Figure FDA0002331674800000026
The definition is as follows:
Figure FDA0002331674800000027
Figure FDA0002331674800000028
Figure FDA0002331674800000029
Figure FDA00023316748000000210
Figure FDA00023316748000000211
Figure FDA00023316748000000212
Aj＝F1,j×p j
in the formula: p is a radical ofj(t) represents the output of the j bidding unit at the time t; v. ofj(t) is a binary variable, the j bidding unit is started to be 1 at the time of t and is shut down to be 0; fl,jBidding for the l section of the j bidding unit; deltal(j, t) represents the I segment output of the j bidding unit at the time t; NLjThe number of quote segments representing j bid units;P jrepresenting the minimum output of the j bidding unit;p jrepresenting the starting point of the output interval of the j bidding unit in the 1 st section;
the starting and shutdown costs of the bidding unit
Figure FDA00023316748000000213
The method is defined as follows, the self characteristics of the hydropower are considered, and the starting and shutdown costs are set to be 0; for a bidding unit of the thermal power generating unit:
Figure FDA00023316748000000214
Figure FDA00023316748000000215
Figure FDA00023316748000000216
Figure FDA0002331674800000031
in the formula:
Figure FDA0002331674800000032
representing the cost of the kth interval of the j stepped starting cost function of the thermal power generating unit; cjRepresenting the shutdown cost of the thermal power generating unit j; NDjRepresenting the number of the starting cost segments of the thermal power generating unit j; and R represents a hydroelectric generating set.
CN201911338730.1A 2019-12-23 2019-12-23 Method for clearing current market of hydropower station water consumption rate iterative correction electric power system before day Active CN111160640B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911338730.1A CN111160640B (en) 2019-12-23 2019-12-23 Method for clearing current market of hydropower station water consumption rate iterative correction electric power system before day

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911338730.1A CN111160640B (en) 2019-12-23 2019-12-23 Method for clearing current market of hydropower station water consumption rate iterative correction electric power system before day

Publications (2)

Publication Number Publication Date
CN111160640A true CN111160640A (en) 2020-05-15
CN111160640B CN111160640B (en) 2023-04-07

Family

ID=70557825

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911338730.1A Active CN111160640B (en) 2019-12-23 2019-12-23 Method for clearing current market of hydropower station water consumption rate iterative correction electric power system before day

Country Status (1)

Country Link
CN (1) CN111160640B (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111861162A (en) * 2020-07-06 2020-10-30 大连理工大学 Method for clearing daily spot goods of high-proportion hydroelectric system under linkage control of step electric quantity
CN111950801A (en) * 2020-08-21 2020-11-17 中国电力科学研究院有限公司 Cross-section interactive day-ahead market clearing method, system, equipment and storage medium
CN111967986A (en) * 2020-05-27 2020-11-20 中国南方电网有限责任公司 Electric power market clearing method, system, device and medium containing cascade hydropower
CN112103943A (en) * 2020-08-13 2020-12-18 云南电网有限责任公司 Safety checking method and device for delivery of electric power spot market in the day-ahead and storage medium
CN112767018A (en) * 2021-01-11 2021-05-07 中国电力科学研究院有限公司 Provincial and intra-provincial power market coupling clearing method, system, device and medium
CN112907078A (en) * 2021-02-20 2021-06-04 中国电力科学研究院有限公司 Market clearing method, system, equipment and readable storage medium for water-electricity coupling
CN113452079A (en) * 2021-06-07 2021-09-28 云南电网有限责任公司 Power generation self-driving system and method considering water and power basin step scheduling
CN113468755A (en) * 2021-07-12 2021-10-01 云南电网有限责任公司 Market iteration clearing method and device based on refined hydraulic connection check
CN113869975A (en) * 2021-09-23 2021-12-31 大连理工大学 Multi-stage clearing method for day-ahead spot market by coupling water and electricity energy consumption
CN115423508A (en) * 2022-08-29 2022-12-02 大连川禾绿能科技有限公司 Strategic bidding method of cascade hydropower in uncertain carbon-electricity coupling market
CN111967986B (en) * 2020-05-27 2024-05-28 中国南方电网有限责任公司 Electric power market clearing method, system, device and medium containing step hydropower

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107563586A (en) * 2017-07-19 2018-01-09 北京微肯佛莱科技有限公司 A kind of ahead market based on step power station coupled relation goes out clear mechanism system
CN110555590A (en) * 2019-07-31 2019-12-10 云南电网有限责任公司 secondary clearing method for watershed cascade upstream and downstream power stations participating in electric power spot market
CN110570008A (en) * 2019-07-31 2019-12-13 云南电网有限责任公司 Day-ahead market clearing method for full acquisition of power generation of down-stream power station in drainage basin cascade

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107563586A (en) * 2017-07-19 2018-01-09 北京微肯佛莱科技有限公司 A kind of ahead market based on step power station coupled relation goes out clear mechanism system
CN110555590A (en) * 2019-07-31 2019-12-10 云南电网有限责任公司 secondary clearing method for watershed cascade upstream and downstream power stations participating in electric power spot market
CN110570008A (en) * 2019-07-31 2019-12-13 云南电网有限责任公司 Day-ahead market clearing method for full acquisition of power generation of down-stream power station in drainage basin cascade

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张粒子等: "多运营主体梯级水电站参与的日前市场出清模型", 《电力系统自动化》 *

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111967986A (en) * 2020-05-27 2020-11-20 中国南方电网有限责任公司 Electric power market clearing method, system, device and medium containing cascade hydropower
CN111967986B (en) * 2020-05-27 2024-05-28 中国南方电网有限责任公司 Electric power market clearing method, system, device and medium containing step hydropower
CN111861162B (en) * 2020-07-06 2023-12-19 大连理工大学 Method for clearing daily current goods of high-proportion hydroelectric system under step electric quantity linkage control
CN111861162A (en) * 2020-07-06 2020-10-30 大连理工大学 Method for clearing daily spot goods of high-proportion hydroelectric system under linkage control of step electric quantity
CN112103943A (en) * 2020-08-13 2020-12-18 云南电网有限责任公司 Safety checking method and device for delivery of electric power spot market in the day-ahead and storage medium
CN111950801A (en) * 2020-08-21 2020-11-17 中国电力科学研究院有限公司 Cross-section interactive day-ahead market clearing method, system, equipment and storage medium
CN112767018A (en) * 2021-01-11 2021-05-07 中国电力科学研究院有限公司 Provincial and intra-provincial power market coupling clearing method, system, device and medium
CN112907078A (en) * 2021-02-20 2021-06-04 中国电力科学研究院有限公司 Market clearing method, system, equipment and readable storage medium for water-electricity coupling
CN113452079A (en) * 2021-06-07 2021-09-28 云南电网有限责任公司 Power generation self-driving system and method considering water and power basin step scheduling
CN113452079B (en) * 2021-06-07 2022-11-29 云南电网有限责任公司 Power generation self-driving system and method considering water and power basin step scheduling
CN113468755A (en) * 2021-07-12 2021-10-01 云南电网有限责任公司 Market iteration clearing method and device based on refined hydraulic connection check
CN113468755B (en) * 2021-07-12 2022-10-04 云南电网有限责任公司 Market iteration clearing method and device based on refined hydraulic connection check
CN113869975A (en) * 2021-09-23 2021-12-31 大连理工大学 Multi-stage clearing method for day-ahead spot market by coupling water and electricity energy consumption
CN113869975B (en) * 2021-09-23 2024-04-16 大连理工大学 Multi-stage clearing method for day-ahead spot market coupled with water and electricity abandonment energy consumption
CN115423508B (en) * 2022-08-29 2023-07-18 大连川禾绿能科技有限公司 Strategy bidding method for cascade hydropower in uncertain carbon-electricity coupling market
CN115423508A (en) * 2022-08-29 2022-12-02 大连川禾绿能科技有限公司 Strategic bidding method of cascade hydropower in uncertain carbon-electricity coupling market

Also Published As

Publication number Publication date
CN111160640B (en) 2023-04-07

Similar Documents

Publication Publication Date Title
CN111160640B (en) Method for clearing current market of hydropower station water consumption rate iterative correction electric power system before day
CN109858759B (en) Industrial park comprehensive energy balance scheduling method
CN108764519B (en) Optimal configuration method for capacity of park energy Internet energy equipment
CN110310173B (en) Electric quantity distribution method for renewable energy to participate in medium-and-long-term electric power transaction
CN109636671B (en) Wind power consumption optimization strategy considering multi-level market linkage type
Yuansheng et al. A two-stage planning and optimization model for water-hydrogen integrated energy system with isolated grid
CN105006844A (en) Electric power system day-ahead robust scheduling system on intermittent power generation grid connected condition
CN109586284B (en) Random production simulation method of transmitting-end power system considering energy curtailment and application
CN110232583B (en) Electric power market marginal price planning method considering carbon emission right
CN115423260A (en) Quantitative analysis method for new energy utilization of electric power market and policy service
Dong et al. China's solar subsidy policy: Government funding yields to open markets
CN109978331B (en) Method for decomposing daily electric quantity in high-proportion water-electricity spot market
Hu et al. A two-stage dispatching method for wind-hydropower-pumped storage integrated power systems
Rizvi et al. Optimal scheduling of virtual power plants utilizing wind power and electric vehicles
Pruckner et al. A study on the impact of packet loss and latency on real-time demand response in smart grid
CN110119850B (en) Heat storage amount two-stage optimization scheduling method based on photo-thermal power generation adjustment
CN114417625B (en) Seasonal energy storage solution method considering wind-gas complementary characteristics
He et al. Research on intraday rolling optimal dispatch including pumped storage power station
CN108022021A (en) A kind of cascade hydropower participate in inside the province with transferring electricity from the west to the east marketing-mix trading strategies optimization method
CN114465226A (en) Method for establishing multi-level standby acquisition joint optimization model of power system
CN114243691A (en) Low-carbon optimized scheduling method for power system
CN114362255A (en) Multi-target day-ahead scheduling optimization method and system for source-network charge storage power system
CN113572165A (en) Electric power spot market clearing mechanism considering water-light storage complementary power generation
CN113394820A (en) Optimized scheduling method for new energy grid-connected power system
Zhou et al. Research on peak shaving auxiliary service market considering the participation of pumped storage power stations

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