CN109858711B - Wind power consumption day-ahead scheduling method considering price type demand response and participation of CSP power station - Google Patents

Wind power consumption day-ahead scheduling method considering price type demand response and participation of CSP power station Download PDF

Info

Publication number
CN109858711B
CN109858711B CN201910172439.5A CN201910172439A CN109858711B CN 109858711 B CN109858711 B CN 109858711B CN 201910172439 A CN201910172439 A CN 201910172439A CN 109858711 B CN109858711 B CN 109858711B
Authority
CN
China
Prior art keywords
power
csp
load
wind power
wind
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
CN201910172439.5A
Other languages
Chinese (zh)
Other versions
CN109858711A (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.)
Northeast Electric Power University
Original Assignee
Northeast Dianli University
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 Northeast Dianli University filed Critical Northeast Dianli University
Priority to CN201910172439.5A priority Critical patent/CN109858711B/en
Publication of CN109858711A publication Critical patent/CN109858711A/en
Application granted granted Critical
Publication of CN109858711B publication Critical patent/CN109858711B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention relates to a day-ahead wind power consumption scheduling method considering price type demand response and CSP power station participation, aiming at the serious problem of wind abandon caused by weak power grid regulation capacity and single price mechanism of a power market, the invention considers the uncertainty of wind power prediction and price type demand response and system safety constraint on the premise of considering comprehensive cost, constructs a scheduling model participating in wind power consumption based on the CSP power station and the price type demand response, combines the CSP power station and the price type demand response with wind power optimization scheduling, and improves the wind power consumption capacity of the system through coordinated scheduling on two sides of source load.

Description

Wind power consumption day-ahead scheduling method considering price type demand response and participation of CSP power station
Technical Field
The invention relates to a day-ahead wind power consumption scheduling method considering price type demand response and participation of a CSP power station.
Background
The renewable energy sources such as wind, light and the like in northwest areas of China are rich, and multi-source complementation becomes a main characteristic of a local new energy power system. However, the wind abandoning phenomenon is increasingly serious due to the comprehensive influence of technical and policy factors such as strong wind wave mobility, weak power grid regulation capability, low load level, single price mechanism of the current power market and the like. How to consider technical and market factors, reduce the amount of abandoned wind power and improve the wind power consumption level is a practical problem which needs to be researched and solved urgently in a large-scale wind power system.
By strengthening a demand side management and user response system (price-based demand response, PDR), the wind power consumption capacity can be improved; meanwhile, renewable energy sources with regulation capacity can be fully utilized for generating electricity, such as: solar photo-thermal power (CSP) realizes output power stabilization and complementation to fluctuating power supplies such as wind power and the like, and further improves the internet surfing space of the wind power. In fact, renewable energy sources with strong scheduling can be utilized on the power supply side, such as: the CSP power station makes up the deficiency of the wind power output anti-peak-shaving characteristic, and improves the wind power consumption together through source-load coordination in cooperation with peak shaving and valley filling of the load side price type demand response.
Disclosure of Invention
The invention aims to solve the technical problem of providing a wind power consumption day-ahead scheduling method which is scientific, reasonable, high in applicability and good in effect and takes price type demand response and CSP power station participation into consideration. The model combines the CSP power station and the price type demand response with the wind power optimization scheduling, and improves the wind power consumption capacity of the system through the coordinated scheduling of the source load and the load.
The technical scheme adopted for solving the technical problem is as follows: a wind power consumption day-ahead scheduling method considering price type demand response and CSP power station participation is characterized in that wind power prediction and price type demand response are considered together with uncertainty and system safety constraint on the premise of considering comprehensive cost, and a scheduling model based on CSP power station and price type demand response participation wind power consumption is constructed; the model combines a CSP power station and a price type demand response with wind power optimization scheduling, improves the capacity of a system for absorbing wind power through coordinated scheduling of two sides of a source load, and specifically comprises the following steps:
1) uncertainty analysis of price type demand response and wind power forecast
(a) Uncertainty of price type demand response
In price type demand response, according to the psychological principle of consumers, when a dispatching center makes a day-ahead dispatching plan, the electricity utilization habits of users are changed to a certain extent by means of real-time electricity prices with time scales of hour levels, and the purposes of peak clipping, valley filling and wind power internet space increasing are achieved;
in the price type demand response, the influence of the electricity price change rate on the load response rate is expressed by using an elasticity coefficient, and the load response rate for a T period is modeled as an expression (1):
Figure BDA0001988372500000021
wherein: phi is aΔq,tFor the response rate of the load for the period t,
Figure BDA0001988372500000022
the change rate of the electricity price of the load in the time interval T, and T is equal to (1: T); e is a price type demand elasticity matrix expressed as formula (2):
Figure BDA0001988372500000023
wherein: epsiloniiCoefficient of self-elasticity,. epsilonijFor cross elasticity coefficients, subscripts i and j denote the ith and jth scheduling periods, respectively;
the uncertainty degree of the load response rate is summarized by adopting a triangular membership function, and a fuzzy expression of the load response rate and a membership parameter have the following relationship and are calculated as a formula (3):
φΔq,t=(φΔql,t,φΔq2,t,φΔq3,t) (3)
Figure BDA0001988372500000024
wherein:
Figure BDA0001988372500000025
is phiΔq,tThe fuzzy expression of (1);
φΔq1,t,φΔq2,t,φΔq3,tmembership degree parameter of load response rate in t time period;
δtfor the maximum error level of the predicted load response rate at time t, it is calculated as (5):
Figure BDA0001988372500000026
wherein: lambda [ alpha ]1A proportionality coefficient before the electricity price factor is dominant;
Figure BDA0001988372500000027
is the inflection point of the rate of change of electricity price;
Figure BDA0001988372500000028
the load change rate corresponding to the electricity price change rate inflection point;
Figure BDA0001988372500000029
the maximum/minimum value of the rate of change of the electricity price;
(b) uncertainty of wind power prediction
Wind power is a power supply which has volatility in output power and is difficult to control, and is a typical intermittent and low-power-density power supply, and wind power prediction often has a certain prediction error; the uncertainty of the prediction error is described by a normal distribution, which is expressed as formula (6):
Figure BDA00019883725000000210
wherein: omegatWind power prediction error at the time t; n represents that the wind power prediction error obeys normal distribution;
Figure BDA0001988372500000031
predicting power for the wind power at the t moment; wmThe installed capacity of the wind power plant;
the actual output of the wind power plant is the sum of the predicted output of the wind power plant and the prediction error, and meanwhile, the actual output of the wind power plant does not exceed the loading capacity of the wind power plant, and the calculation is as the formula (7):
Figure BDA0001988372500000032
wherein:
Figure BDA0001988372500000033
real output of the wind power plant at the moment t is obtained;
2) establishment of CSP power station heat-electricity conversion model
The solar mirror field collects heat energy generated by solar radiation to the heat collection tower, and the heat energy collected by the heat collection tower is calculated as the formula (8):
Figure BDA0001988372500000034
wherein:
Figure BDA0001988372500000035
the heat energy collected by the heat collection tower at the moment t is obtained; etad-thTo the light-to-heat conversion efficiency; sSFIs the area of a solar mirror field;
Dtrepresents the solar direct radiation index (DNI) at time t;
the heat energy collected by the heat collecting tower can be stored in the heat storage device through the heat transfer fluid and also can be directly supplied to a power generation system for power generation, and the heat energy collected by the heat collecting tower has the possibility of being abandoned in order to ensure the stable operation of the CSP power station; meanwhile, the heat energy used for power generation of the power generation system can also be supplemented by the heat storage system, and is expressed as the formula (9):
Figure BDA0001988372500000036
wherein:
Figure BDA0001988372500000037
representing the heat storage device heat release efficiency;
Figure BDA0001988372500000038
represents the heat energy rejected by the heat collection tower at the moment t; eSURepresents the thermal energy consumed by the power generation system during startup; u. oftA start variable indicating the power generation system, 1 indicating that the power generation system is started at that time;
the heat storage amount of the heat storage device at each time is expressed by the following formula (10):
Figure BDA0001988372500000039
wherein:
Figure BDA00019883725000000310
respectively representing the heat storage amount of the heat storage device at the moment t and the moment t-1; etalossIndicating a rate of heat loss from the heat storage device;
Figure BDA00019883725000000311
represents the charging efficiency of the heat storage device;
Figure BDA00019883725000000312
representing the charging power of the heat storage device;
Figure BDA00019883725000000313
representing the heat release power of the heat storage device; Δ t represents a time interval;
the thermodynamic dynamic differential equation is not considered in the static model of the CSP power station, and the power generation power and the heat energy consumption power of the CSP power station can be regarded as a linear relation and are calculated as the formula (11):
Figure BDA00019883725000000314
wherein:
Figure BDA00019883725000000315
representing the generated power of the CSP power station for supplying load at the time t; etaRs,uThe positive rotation standby coefficient of the CSP power station;
Figure BDA00019883725000000316
providing positive rotation standby for the CSP power station at the time t; kappaCSPThe heat-electricity conversion efficiency of the CSP power station power generation system is shown;
Figure BDA00019883725000000317
representing the heat energy consumed by the power generation system at time t;
3) establishment of day-ahead scheduling model considering PDR and CSP power station participating in wind power consumption
a) Objective function
The method comprises the following steps of comprehensively considering the power generation cost of a conventional thermal power generating unit, the operation and maintenance cost of a wind power and CSP power station, the renewable energy environmental benefit brought by grid-connected consumption, the cost of the thermal power generating unit and the CSP power station participating in load standby and wind power prediction error standby, establishing a combined scheduling model containing the wind power, the CSP power station and the thermal power station and based on optimal cost, and calculating a target function as a formula (12):
min F=D1+D2+D3+D4-D5 (12)
wherein: f is the comprehensive cost of the system when thermal power, wind power and CSP power stations participate in the optimized dispatching; d1The power generation cost when the load is supplied to the thermal power generating unit; d2The operation and maintenance cost of the wind power is reduced; d3The operation and maintenance cost for supplying power to the load by the CSP power station; d4The cost for the thermal power generating unit and the CSP power station to participate in load and wind power prediction error standby is saved; d5Environmental benefits brought to the power generation grid-connected consumption of wind power and CSP power stations;
the output of the thermal power generating unit is flexible and controllable, and the stable operation of a power grid can be ensured through reasonable optimized scheduling; in order to meet the load requirement in the scheduling process, the output is often required to be adjusted, even the unit is scheduled to be started and stopped, the power generation cost at the moment mainly comprises the coal consumption cost of the thermal power unit and the start and stop cost of the unit, and the calculation is as the formulas (13) and (14):
Figure BDA0001988372500000041
fi(PGi,t)=aiPGi,t 2+biPGi,t+ci (14)
wherein: t is the total duration; n is a radical ofGThe number of the thermal power generating units is; f. ofiThe coal consumption cost of the thermal power generating unit i is obtained; pGi,tThe output of the thermal power generating unit i in the time period t is obtained; v isi,tV and vi,t-1Respectively the running states of the thermal power generating unit i at t and t-1, if vi,tV denotes unit operation as 1i,tWhen the unit is stopped, 0 represents the unit is stopped; siThe start-stop cost of the unit i is calculated; a isi,bi,ciThe coal consumption cost parameter is the coal consumption cost parameter of the thermal power generating unit i;
wind power generation belongs to renewable energy power generation, coal is not consumed in the power generation process, but certain operation and maintenance cost can be generated in the power generation process of a fan due to uncertainty of wind speed; the operation maintenance cost of the wind power and the wind power output power can be approximately regarded as a linear relation, and the linear relation is calculated as an expression (15):
Figure BDA0001988372500000042
wherein: kWOperating and maintaining costs for the wind farm; pW,tGenerating power of the wind power plant at the time t;
the CSP power station operation and maintenance cost is approximately regarded as a linear function of the generated power, and is calculated as the formula (16):
Figure BDA0001988372500000043
wherein: kCSPAs a sheet of CSP power stationA bit run maintenance cost;
in order to deal with uncertainty and emergency of load and wind power prediction, a certain rotation reserve capacity needs to be reserved, and the reserve cost of a thermal power unit and a CSP power station is calculated as a formula (17):
Figure BDA0001988372500000044
wherein:
Figure BDA0001988372500000045
indicating the participation of the CSP power station in the load standby,
Figure BDA0001988372500000046
Representing the cost coefficient of the CSP power station participating in wind power standby;
Figure BDA0001988372500000051
indicating a positive spinning reserve provided by the CSP station for the load at time t,
Figure BDA0001988372500000052
Representing positive rotation standby power provided by the CSP power station for wind power at the time t;
Figure BDA0001988372500000053
representing the cost coefficient of the thermal power generating unit participating in the load standby,
Figure BDA0001988372500000054
Representing the cost coefficient of the thermal power generating unit participating in wind power standby;
Figure BDA0001988372500000055
representing positive rotation standby provided for the load by the thermal power generating unit i at the moment t;
Figure BDA0001988372500000056
representing positive power provided by thermal power generating unit i for wind power at time tRotating the standby power;
wind power, CSP power station renewable energy are incorporated into the power networks and are consumed and have reduced thermal power unit's the power generation on-line volume, effectively reduce the emission of pollutants such as sulfur and sodium sulfate, can bring environmental benefit, calculate (18) formula:
wherein: pW,tPower is consumed for the grid connection of the wind power plant at the time t; rhoCSPEnvironmental benefit coefficient brought to CSP power station grid-connected consumption; rhoWEnvironmental benefit coefficient brought to wind power grid connection consumption;
b) system constraints
When the network transmission loss is not counted, the sum of the output power of each unit is equal to the value after the power grid load response change, and is expressed as a formula (19):
Figure BDA0001988372500000057
wherein: l istThe load power value before the system needs to respond at the time t;
Figure BDA0001988372500000058
the load response rate triangular fuzzy number is converted into a deterministic variable, and then the expected value of the load response quantity at the time t is expressed as a formula (20):
Figure BDA0001988372500000059
the transmission capacity constraint of the transmission line is expressed as formula (21):
Figure BDA00019883725000000510
wherein: pij,maxIs the maximum transmission capacity of the transmission line between nodes i and j; b isijIs the susceptance between nodes ij; thetai,tIs the voltage phase angle of node i; thetaj,tIs the voltage phase angle of node j;
the opportunity constraint is adopted to determine the wind power reserve capacity, the wind power reserve cost is reduced on the basis of ensuring the safety of a power grid, and the concrete constraint is expressed as a formula (22):
Figure BDA00019883725000000511
wherein:
Figure BDA00019883725000000512
the positive rotation reserve capacity of the thermal power generating unit i at the moment t is obtained;
Figure BDA00019883725000000513
the negative rotation reserve capacity of the thermal power generating unit i at the moment t is obtained;
Figure BDA00019883725000000514
positive rotation reserve capacity provided for CSP power station at time t,
Figure BDA00019883725000000515
The negative spinning reserve capacity provided by the CSP plant at time t,
Figure BDA00019883725000000516
μLis the load reserve factor; alpha and beta are respectively confidence coefficients meeting the positive and negative rotation standby constraints;
the heat storage device of the CSP power station has to be within the rated limit range, and the heat storage device cannot be simultaneously charged/discharged in each scheduling period, and the specific constraint is expressed as the following formula (23) and formula (24):
Figure BDA0001988372500000061
Figure BDA0001988372500000062
wherein:
Figure BDA0001988372500000063
the maximum charging power of the heat storage system;
Figure BDA0001988372500000064
the maximum heat release power of the heat storage system;
the heat storage quantity of the CSP power station heat storage device is expressed as a formula (25):
Figure BDA0001988372500000065
wherein:
Figure BDA0001988372500000066
the minimum heat storage quantity of the heat storage device is obtained; xiTSThe maximum heat storage capacity of the heat storage device is expressed by taking FLH as a unit;
the specific constraint of the maximum/minimum output of the thermal power generating unit is expressed as a formula (26):
UGi,tPi,min≤PGi,t≤UGi,tPi,max (26)
wherein: u shapeGi,tThe method comprises the following steps that 1 represents operation, and 0 represents shutdown, wherein the thermal power generating unit i is in an operating state; pi,minThe minimum output of the thermal power generating unit i is obtained; p isi,maxThe maximum output of the thermal power generating unit i is obtained;
the ramp rate constraint expression of the thermal power generating unit is as follows (27):
-rdi≤PGi,t-PGi,t-1≤rui (27)
during day-ahead scheduling, the wind power on-line power cannot exceed the predicted value, and the specific expression is (28):
Figure BDA0001988372500000067
wherein:
Figure BDA0001988372500000068
for wind farms at time tThe predicted output of (2);
the sum of the expected load response values after the demand response is zero in the whole scheduling period, and is expressed as formula (29):
Figure BDA0001988372500000069
considering the benefit of the user, the load variation needs to be limited, and the satisfaction of the user power utilization mode and the satisfaction of the power utilization expense are used as measurement indexes and expressed as formulas (30) and (31):
Figure BDA00019883725000000610
Figure BDA00019883725000000611
wherein:
Figure BDA00019883725000000612
the minimum value of the satisfaction degree of the user power utilization mode is set;
Figure BDA00019883725000000613
the minimum value of the satisfaction degree of electricity cost expenditure;
Figure BDA00019883725000000614
the load value after the demand response at the time t is calculated as (32):
Figure BDA0001988372500000071
the load value after the price type demand response is between the upper limit and the lower limit of the load value before the response, and the specific constraint expression is (33):
LminminLt≤E(Δqt)≤LmaxmaxLt (33)
wherein: l ismaxIs the maximum value, L, of the original load curveminIs the minimum value of the original load curve; etamax、ηminTaking the value eta as the peak-to-valley difference coefficient of the demand response loadmax≥1,ηmin≤1。
The effect of the day-ahead wind power consumption scheduling method considering the price type demand response and the participation of the CSP power station is represented as follows: on the premise of considering the comprehensive cost, the uncertainty of wind power prediction and price type demand response and system safety constraints are considered, a scheduling model based on the CSP power station and the price type demand response participating in wind power consumption is constructed, and the capacity of the system for consuming wind power is improved through coordinated scheduling of the source load and the load. The method is scientific and reasonable, and has strong applicability and good effect.
Drawings
FIG. 1 is a 30-node wiring diagram of the modified IEEE;
FIG. 2CSP power plant components and energy flow diagrams;
FIG. 3 illustrates the pre-and post-response loads of the price type demand and the performance diagrams of various types of units;
FIG. 4 is a real-time electricity price graph for each dispatch period after a price type demand response;
FIG. 5 illustrates the heat storage device for heat and charge conditions and the amount of heat stored during each scheduled time period;
FIG. 6CSP power station thermal power consumption per period.
Detailed Description
The present invention will be further explained with reference to the drawings and embodiments, in which a method for scheduling wind power generation in consideration of the price type demand response and participation of CSP power stations in the day ahead is described.
The invention relates to a day-ahead wind power consumption scheduling method considering price type demand response and participation of a CSP power station, which is characterized in that uncertainty and system safety constraint of wind power prediction and price type demand response are considered on the premise of considering comprehensive cost, and a scheduling model based on CSP power station and price type demand response participation wind power consumption is constructed; the model combines a CSP power station and a price type demand response with wind power optimization scheduling, improves the wind power consumption capability of the system through coordinated scheduling of source load two sides, and specifically comprises the following steps:
1) uncertainty analysis of price type demand response and wind power forecast
(a) Uncertainty of price type demand response
In price type demand response, according to the psychological principle of consumers, when a dispatching center makes a day-ahead dispatching plan, the electricity utilization habits of users are changed to a certain extent by means of real-time electricity prices with time scales of hour levels, and the purposes of peak clipping, valley filling and wind power internet space increasing are achieved;
in the price type demand response, the influence of the electricity price change rate on the load response rate is expressed by using an elasticity coefficient, and the load response rate for a T period is modeled as an expression (1):
Figure BDA0001988372500000081
wherein: phi is aΔq,tFor the response rate of the load for the period t,
Figure BDA0001988372500000082
the change rate of the electricity price of the load in the time interval T, and T is equal to (1: T);
e is a price type demand elasticity matrix expressed as formula (2):
Figure BDA0001988372500000083
wherein: epsiloniiCoefficient of self-elasticity,. epsilonijFor cross elasticity coefficients, subscripts i and j denote the ith and jth scheduling periods, respectively;
the uncertainty degree of the load response rate is summarized by adopting a triangular membership function, and a fuzzy expression of the load response rate and a membership parameter have the following relationship and are calculated as a formula (3):
φΔq,t=(φΔq1,t,φΔq2,t,φΔq3,t) (3)
Figure BDA0001988372500000084
wherein:
Figure BDA0001988372500000085
is phiΔq,tThe fuzzy expression of (1);
φΔq1,t,φΔq2,t,φΔq3,tmembership parameter of load response rate in t time period;
δtfor the maximum error level of the predicted load response rate at time t, it is calculated as (5):
Figure BDA0001988372500000086
wherein: lambda [ alpha ]1A proportionality coefficient before the electricity price factor is dominant;
Figure BDA0001988372500000087
is the inflection point of the rate of change of electricity price;
Figure BDA0001988372500000088
the load change rate corresponding to the electricity price change rate inflection point;
Figure BDA0001988372500000089
the maximum/minimum value of the rate of change of the electricity price;
(b) uncertainty of wind power prediction
Wind power is a power supply which has volatility in output power and is difficult to control, and is a typical intermittent and low-power-density power supply, and wind power prediction often has a certain prediction error; the uncertainty of the prediction error is described by a normal distribution, which is expressed as formula (6):
Figure BDA00019883725000000810
wherein: omegatWind power prediction error at the time t; n represents that the wind power prediction error obeys normal distribution;
Figure BDA00019883725000000811
predicting the power of the wind power at the time t; wmIs the installed capacity of the wind farm;
the actual output of the wind power plant is the sum of the predicted output of the wind power plant and the prediction error, and meanwhile, the actual output of the wind power plant does not exceed the loading capacity of the wind power plant, and the calculation is as the formula (7):
Figure BDA0001988372500000091
wherein:
Figure BDA0001988372500000092
real output of the wind power plant at the time t is obtained;
2) establishment of CSP power station heat-electricity conversion model
The solar mirror field collects heat energy generated by solar radiation to the heat collection tower, and the heat energy collected by the heat collection tower is calculated as the formula (8):
Figure BDA0001988372500000093
wherein:
Figure BDA0001988372500000094
the heat energy collected by the heat collection tower at the moment t is obtained; etad-thTo the light-to-heat conversion efficiency; sSFIs the area of a solar mirror field;
Dtrepresents the solar direct radiation index (DNI) at time t;
the heat energy collected by the heat collecting tower can be stored in the heat storage device through the heat transfer fluid and also can be directly supplied to a power generation system for power generation, and the heat energy collected by the heat collecting tower has the possibility of being abandoned in order to ensure the stable operation of the CSP power station; meanwhile, the heat energy used for power generation of the power generation system can also be supplemented by the heat storage system, and is expressed as the formula (9):
Figure BDA0001988372500000095
wherein:
Figure BDA0001988372500000096
represents the heat-storage device heat-release efficiency;
Figure BDA0001988372500000097
represents the heat energy rejected by the heat collection tower at the moment t; eSURepresents the thermal energy consumed by the power generation system during startup; u. oftA start variable indicating the power generation system, 1 indicating that the power generation system is started at that time;
the heat storage amount of the heat storage device at each time is expressed by the following formula (10):
Figure BDA0001988372500000098
wherein:
Figure BDA0001988372500000099
respectively representing the heat storage amount of the heat storage device at the moment t and the moment t-1; etalossIndicating a rate of heat loss from the heat storage device;
Figure BDA00019883725000000910
represents the charging efficiency of the heat storage device;
Figure BDA00019883725000000911
representing the charging power of the heat storage device;
Figure BDA00019883725000000912
representing the heat release power of the heat storage device; Δ t represents a time interval;
the thermodynamic dynamic differential equation is not considered in the static model of the CSP power station, and the power generation power and the heat energy consumption power of the CSP power station can be regarded as a linear relation and are calculated as the formula (11):
Figure BDA00019883725000000913
wherein:
Figure BDA00019883725000000914
representing the generated power of the CSP power station for supplying load at the time t; etaRs,uThe positive rotation standby coefficient of the CSP power station;
Figure BDA00019883725000000915
providing positive rotation standby for the CSP power station at the time t; kappa typeCSPThe heat-electricity conversion efficiency of the CSP power station power generation system is shown;
Figure BDA00019883725000000916
representing the heat energy consumed by the power generation system at time t;
3) establishment of day-ahead scheduling model considering participation of PDR (plant data Rate) and CSP (compact strip production) power station in wind power consumption
a) Objective function
The method comprises the following steps of comprehensively considering the power generation cost of a conventional thermal power generating unit, the operation and maintenance cost of a wind power and CSP power station, the renewable energy environmental benefit brought by grid-connected consumption, the cost of the thermal power generating unit and the CSP power station participating in load standby and wind power prediction error standby, establishing a combined scheduling model containing the wind power, the CSP power station and the thermal power station and based on optimal cost, and calculating a target function as a formula (12):
min F=D1+D2+D3+D4-D5 (12)
wherein: f is the comprehensive cost of the system when the thermal power, wind power and CSP power stations participate in the optimized scheduling; d1The power generation cost when the load is supplied to the thermal power generating unit; d2The operation and maintenance cost of the wind power is reduced; d3The operation and maintenance cost for supplying power to the load by the CSP power station; d4The cost for the thermal power generating unit and the CSP power station to participate in load and wind power prediction error standby is saved; d5Environmental benefits brought to the power generation grid-connected consumption of wind power and CSP power stations;
the output of the thermal power generating unit is flexible and controllable, and the stable operation of a power grid can be ensured through reasonable optimized scheduling; in order to meet the load requirement in the scheduling process, the output is often required to be adjusted, even the unit is scheduled to be started and stopped, the power generation cost at the moment mainly comprises the coal consumption cost of the thermal power unit and the start and stop cost of the unit, and the calculation is as the formulas (13) and (14):
Figure BDA0001988372500000101
fi(PGi,t)=aiPGi,t 2+biPGi,t+ci (14)
wherein: t is the total duration; n is a radical ofGThe number of the thermal power generating units is; f. ofiThe coal consumption cost of the thermal power generating unit i is obtained; pGi,tThe output of the thermal power generating unit i in the time period t is obtained; v isi,tV and vi,t-1Respectively the running states of the thermal power generating unit i at t and t-1, if vi,tV denotes unit operation as 1i,tWhen the unit is stopped, 0 represents the unit is stopped; siThe start-stop cost of the unit i is calculated; a isi,bi,ciThe coal consumption cost parameter is the coal consumption cost parameter of the thermal power generating unit i;
wind power generation belongs to renewable energy power generation, coal is not consumed in the power generation process, but certain operation and maintenance cost can be generated in the power generation process of a fan due to uncertainty of wind speed; the operation and maintenance cost of the wind power and the wind power output power can be approximately regarded as a linear relation, and are calculated as an expression (15):
Figure BDA0001988372500000102
wherein: k isWOperating and maintaining costs for the wind farm; p isW,tGenerating power of the wind power plant at the time t;
the operation and maintenance cost of the CSP power station is approximately regarded as a linear function of the generated power, and is calculated as the formula (16):
Figure BDA0001988372500000103
wherein: kCSPThe unit operation and maintenance cost of the CSP power station is saved;
in order to deal with uncertainty and emergency of load and wind power prediction, a certain rotation reserve capacity needs to be reserved, and the reserve cost of a thermal power unit and a CSP power station is calculated as a formula (17):
Figure BDA0001988372500000104
wherein:
Figure BDA0001988372500000105
indicating the participation of the CSP power station in the load standby,
Figure BDA0001988372500000106
Representing the cost coefficient of the CSP power station participating in wind power standby;
Figure BDA0001988372500000107
indicating a positive spinning reserve provided by the CSP plant for the load at time t,
Figure BDA0001988372500000108
Representing positive rotation standby power provided by the CSP power station for wind power at the moment t;
Figure BDA0001988372500000111
representing the cost coefficient of the thermal power generating unit participating in the load standby,
Figure BDA0001988372500000112
Representing the cost coefficient of the thermal power generating unit participating in wind power standby;
Figure BDA0001988372500000113
indicating that thermal power unit i is providing load at time tRotating positively for standby;
Figure BDA0001988372500000114
representing positive rotation standby power provided by the thermal power generating unit i for wind power at the time t;
wind power, CSP power station renewable energy are incorporated into the power networks and are consumed and have reduced thermal power unit's the power generation on-line volume, effectively reduce the emission of pollutants such as sulfur and sodium sulfate, can bring environmental benefit, calculate (18) formula:
wherein: pW,tPower is consumed for the grid connection of the wind power plant at the time t; rhoCSPEnvironmental benefit coefficient brought to CSP power station grid-connected consumption; rhoWEnvironmental benefit coefficient brought to wind power grid connection consumption;
b) system constraints
When the network transmission loss is not counted, the sum of the output power of each unit is equal to the value after the power grid load response change, and is expressed as a formula (19):
Figure BDA0001988372500000115
wherein: l istThe load power value before the system needs to respond at the time t;
Figure BDA0001988372500000116
the load response rate triangular fuzzy number is converted into a deterministic variable, and then the expected value of the load response quantity at the time t is expressed as a formula (20):
Figure BDA0001988372500000117
the transmission capacity constraint of the transmission line is expressed as formula (21):
Figure BDA0001988372500000118
wherein: pij,maxIs the maximum of the transmission line between nodes i and jA transmission capacity; b isijIs the susceptance between nodes ij; thetai,tIs the voltage phase angle of node i; thetaj,tIs the voltage phase angle of node j;
the opportunity constraint is adopted to determine the wind power reserve capacity, the wind power reserve cost is reduced on the basis of ensuring the safety of a power grid, and the concrete constraint is expressed as a formula (22):
Figure BDA0001988372500000119
wherein:
Figure BDA00019883725000001110
the positive rotation reserve capacity of the thermal power generating unit i at the moment t is obtained;
Figure BDA00019883725000001111
the negative rotation reserve capacity of the thermal power generating unit i at the moment t is obtained;
Figure BDA00019883725000001112
positive rotation reserve capacity provided for CSP power station at time t,
Figure BDA00019883725000001113
The negative spinning reserve capacity provided by the CSP plant at time t,
Figure BDA00019883725000001114
μLis the load reserve factor; alpha and beta are respectively confidence coefficients meeting the positive and negative rotation standby constraints;
the heat storage device of the CSP power station has to be within the rated limit range, and the heat storage device cannot be simultaneously charged/discharged in each scheduling period, and the specific constraint is expressed as the following formula (23) and formula (24):
Figure BDA0001988372500000121
Figure BDA0001988372500000122
wherein:
Figure BDA0001988372500000123
the maximum charging power of the heat storage system;
Figure BDA0001988372500000124
the maximum heat release power of the heat storage system;
the heat storage quantity of the CSP power station heat storage device is expressed as a formula (25):
Figure BDA0001988372500000125
wherein:
Figure BDA0001988372500000126
the minimum heat storage quantity of the heat storage device is obtained; xi shapeTSThe maximum heat storage capacity of the heat storage device is expressed by taking FLH as a unit;
the specific constraint of the maximum/minimum output of the thermal power generating unit is expressed as a formula (26):
UGi,tPi,min≤PGi,t≤UGi,tPi,max (26)
wherein: u shapeGi,tThe method comprises the following steps that 1 represents operation, and 0 represents shutdown, wherein the thermal power generating unit i is in an operating state; pi,minThe minimum output of the thermal power generating unit i is obtained; pi,maxThe maximum output of the thermal power generating unit i is obtained;
the ramp rate constraint expression of the thermal power generating unit is as follows (27):
-rdi≤PGi,t-PGi,t-1≤rui (27)
during day-ahead scheduling, the wind power on-line power cannot exceed the predicted value, and the specific expression is (28):
Figure BDA0001988372500000127
wherein:
Figure BDA0001988372500000128
the predicted output of the wind power plant at the time t is obtained;
the sum of the expected load response values after the demand response is zero in the whole scheduling period, and is expressed as formula (29):
Figure BDA0001988372500000129
considering the benefit of the user, the load variation needs to be limited, and the satisfaction of the user power utilization mode and the satisfaction of the power utilization expense are used as measurement indexes and expressed as formulas (30) and (31):
Figure BDA00019883725000001210
Figure BDA00019883725000001211
wherein:
Figure BDA00019883725000001212
the minimum value of the satisfaction degree of the power utilization mode of the user is set;
Figure BDA00019883725000001213
the minimum value of the satisfaction degree of electricity cost expenditure;
Figure BDA00019883725000001214
the load value after the demand response at the time t is calculated as (32):
Figure BDA00019883725000001215
the load value after the price type demand response is between the upper limit and the lower limit of the load value before the response, and the specific constraint expression is (33):
LminminLt≤E(Δqt)≤LmaxmaxLt (33)
wherein: l ismaxIs the maximum value, L, of the original load curveminIs the minimum value of the original load curve; etamax、ηminTaking the value eta as the peak-to-valley difference coefficient of the demand response loadmax≥1,ηmin≤1。
The embodiment takes an improved IEEE-30 node system as an example, a CSP power station and a price type demand response are calculated by combining with wind power optimization scheduling, the system comprises 6 conventional thermal power units, a wind power plant and a CSP power station, and a system wiring diagram is shown in figure 1; the data are from measured data of wind power plants and CSP power plants, and the data can be obtained by using a data acquisition device of a commercial product familiar to the technical personnel in the field. The CSP plant main parameters are shown in table 1.
Example the calculation conditions are illustrated below:
1) the electricity price before price type demand response is 400 yuan/MWh;
2) coefficient of proportionality λ1=0.5;
3) Inflection point of electricity price change rate
Figure BDA0001988372500000131
The value is +/-0.3, and the maximum/minimum value of the electricity price change rate is +/-0.5;
4) load reserve factor μLThe value is 0.1;
5) cost of CSP plant participation in load and wind backup
Figure BDA0001988372500000132
6) Cost coefficient K of CSP power station power generation operation maintenanceCSP40 yuan/MW;
7) wind power operation maintenance cost KWIs 30 yuan/MW;
8) environmental benefit coefficient rho generated by grid-connected consumption of wind power and CSP power stationCSP=ρW120 yuan/WM;
9) the cost coefficient pi G L-pi G W-120 yuan/MW of the thermal power generating unit participating in load and wind power backup;
10) minimum value of satisfaction degree of user power utilization mode
Figure BDA0001988372500000133
The value is 0.9, and the minimum value of the satisfaction degree of electricity cost expenditure
Figure BDA0001988372500000134
The value is 1;
10) peak to valley difference coefficient ηmax、ηminThe values are 1.05 and 0.95 respectively;
11) the confidence alpha is 95% that satisfies the positive and negative rotation standby constraints.
TABLE 1CSP Power station parameters
Figure BDA0001988372500000135
Figure BDA0001988372500000141
Under the above calculation conditions, the results of the day-ahead scheduling calculation of the price type demand response and the wind power consumption participated in by the CSP power station by applying the method of the invention are as follows:
1. uncertainty analysis of price type demand response and wind power forecast
In price type demand response, according to the psychological principle of consumers, when a scheduling center makes a day-ahead scheduling plan, the electricity utilization habits of users can be changed to a certain extent by means of real-time electricity prices with time scales of small levels, and therefore the purposes of cutting peaks, filling valleys and increasing wind power internet space are achieved, load response has uncertainty, and the influence of the uncertainty on scheduling results is analyzed in scheduling.
Electricity is a power supply which has volatility in output power and is not easy to control, and is a typical intermittent and low-power-density power supply, wind power prediction often has a certain prediction error, and uncertainty of the prediction error is described by normal distribution.
Establishment of CSP power station thermal-electric conversion model
The CSP power station components and energy flow diagram are shown in FIG. 2, the CSP power station containing the heat storage device mainly comprises a heat collection system, namely: the solar energy mirror field and the heat collecting tower, the heat storage device and the power generation system. The CSP electric station has good energy time shifting characteristics due to the arrangement of the heat storage device, and the defects of intermittency and fluctuation of solar energy are effectively overcome, so that the CSP electric station can be guaranteed to provide stable and controllable power supply. The CSP power station utilizes the heat collected by the solar mirror field to generate power through a certain energy flow exchange process and finally through Rankine cycle.
3. Establishment and solution of day-ahead scheduling model considering price type demand response and participation of CSP power station in wind power consumption
The method comprises the steps of comprehensively considering the power generation cost of a conventional thermal power generating unit, the operation and maintenance cost of a wind power and CSP power station, the renewable energy environmental benefit brought by grid-connected consumption, and the cost of the thermal power generating unit and the CSP power station participating in load standby and wind power prediction error standby, and establishing a combined scheduling model containing the wind power, the CSP power station and the thermal power station and based on optimal cost.
The price type demand response guides the user to change the electricity utilization behavior by changing the electricity price according to the economic principle, and has the functions of peak clipping and valley filling. The load before and after the demand response and the output situation of each type of unit are shown in fig. 3, and the real-time electricity price after the response is shown in fig. 4. As can be seen from fig. 3 and 4, the wind power has a back-peak regulation property, the price type demand response changes the electricity price at each moment, and the electricity usage habits of the users are influenced by the electricity price, so that the peak clipping and valley filling effects are achieved, and the grid-connected consumption of the wind power is increased to a certain extent. However, the thermal power generating unit has certain minimum output and still generates partial abandoned wind in the overlapping time period of the load valley and the wind power high-power generation due to the limitation of the minimum output and the negative rotation standby constraint of the thermal power generating unit.
The output of the CSP power station provided with the heat storage device is flexible and controllable, the output of a thermal power generating unit can be reduced while the wind power anti-peaking property is compensated, and further the environmental benefit is increased. The heat charging and discharging conditions of the CSP power station heat storage device and the heat storage quantity in each scheduling period are shown in FIG. 5.
As can be seen from fig. 5, the CSP power station has the characteristic of energy time shifting, the heat storage device releases the heat stored in the last period of the previous scheduling day in the morning for the CSP power station to output, the heat storage is performed in the period of the day when the solar radiation resources are sufficient, and the stored heat is released at night, thereby ensuring the stable output of the CSP power station.
The output of the CSP power station is stable and controllable, the CSP power station is used for providing standby for wind power and load prediction errors, and the CSP power station makes up the defect of prediction precision through renewable energy sources and ensures the stable operation of a power grid. The specific thermal power consumption of the CSP power station on the scheduled day is shown in FIG. 6.
The terms of calculation, illustration and the like in the embodiments of the present invention are used for further description, are not exhaustive, and do not limit the scope of the claims, and those skilled in the art can conceive other substantially equivalent alternatives without creative efforts based on the teachings of the embodiments of the present invention, which are within the scope of the present invention.

Claims (1)

1. A day-ahead wind power consumption scheduling method considering price type demand response and participation of a CSP power station is characterized in that uncertainty and system safety constraint of wind power prediction and price type demand response are considered under the premise of considering comprehensive cost, and a scheduling model based on the CSP power station and the price type demand response and participating in wind power consumption is constructed; the model combines a CSP power station and a price type demand response with wind power optimization scheduling, improves the capacity of a system for absorbing wind power through coordinated scheduling of two sides of a source load, and specifically comprises the following steps:
1) uncertainty analysis of price type demand response and wind power forecast
(a) Uncertainty of price type demand response
In price type demand response, according to the psychological principle of consumers, when a scheduling center makes a day-ahead scheduling plan, the electricity utilization habits of users can be changed by means of real-time electricity prices with time scales of hour levels, and the purposes of peak clipping, valley filling and wind power online space increasing are achieved;
in the price type demand response, the influence of the electricity price change rate on the load response rate is expressed by using an elasticity coefficient, and the load response rate for a T period is modeled as an expression (1):
Figure FDA0003529538990000011
wherein: phi is a unit ofΔq,tFor the response rate of the load for the period t,
Figure FDA0003529538990000012
the change rate of the electricity price of the load in the time interval T, and T is equal to (1: T);
e is a price type demand elasticity matrix expressed as formula (2):
Figure FDA0003529538990000013
wherein: epsiloniiCoefficient of self-elasticity,. epsilonijFor cross elasticity coefficients, subscripts i and j denote the ith and jth scheduling periods, respectively;
the uncertainty degree of the load response rate is summarized by adopting a triangular membership function, and a fuzzy expression of the load response rate and a membership parameter have the following relationship and are calculated as a formula (3):
Figure FDA0003529538990000014
Figure FDA0003529538990000015
wherein:
Figure FDA0003529538990000016
is phiΔq,tThe fuzzy expression of (1);
φΔq1,t,φΔq2,t,φΔq3,tmembership degree parameter of load response rate in t time period;
δtfor the maximum error level of the predicted load response rate at time t, it is calculated as (5):
Figure FDA0003529538990000021
wherein: lambda1A proportionality coefficient before the electricity price factor is dominant;
Figure FDA0003529538990000022
is the inflection point of the rate of change of electricity price;
Figure FDA0003529538990000023
the load change rate is corresponding to the inflection point of the electricity price change rate;
Figure FDA0003529538990000024
the maximum/minimum value of the rate of change of the electricity price;
(b) uncertainty of wind power prediction
Wind power is a power supply which has volatility in output power and is difficult to control, and is a typical intermittent and low-power-density power supply, and prediction errors exist in wind power prediction; the uncertainty of the prediction error is described by a normal distribution, which is expressed as formula (6):
Figure FDA0003529538990000025
wherein: omegatWind power prediction error at the time t; n represents that the wind power prediction error obeys normal distribution;
Figure FDA0003529538990000026
predicting the power of the wind power at the time t; wmFor installations in wind farmsMachine capacity;
the actual output of the wind power plant is the sum of the predicted output of the wind power plant and the prediction error, and meanwhile, the actual output of the wind power plant does not exceed the installed capacity, and the calculation is as the formula (7):
Figure FDA0003529538990000027
wherein:
Figure FDA0003529538990000028
real output of the wind power plant at the time t is obtained;
2) establishment of CSP power station heat-electricity conversion model
The solar mirror field collects heat energy generated by solar radiation to the heat collection tower, and the heat energy collected by the heat collection tower is calculated as the formula (8):
Figure FDA0003529538990000029
wherein:
Figure FDA00035295389900000210
the heat energy collected by the heat collection tower at the moment t is obtained; etad-thTo the light-to-heat conversion efficiency; sSFIs the area of a solar mirror field; dtRepresenting the direct solar radiation index at the moment t;
the heat energy collected by the heat collecting tower can be stored in the heat storage device through the heat transfer fluid and also can be directly supplied to the power generation system for power generation, so that the heat energy collected by the heat collecting tower has the possibility of being discarded in order to ensure the stable operation of the CSP power station; meanwhile, the heat energy used for generating electricity by the power generation system can also be supplemented by the heat storage system, and is expressed as the formula (9):
Figure FDA00035295389900000211
wherein:
Figure FDA00035295389900000212
representing the heat storage device heat release efficiency;
Figure FDA00035295389900000213
represents the heat energy rejected by the heat collection tower at the moment t; eSURepresents the thermal energy consumed by the power generation system during startup; u. oftA start variable indicating the power generation system, 1 indicating that the power generation system is started at that time;
the heat storage amount of the heat storage device at each time is expressed by the following formula (10):
Figure FDA00035295389900000214
wherein:
Figure FDA00035295389900000215
respectively representing the heat storage amount of the heat storage device at the moment t and the moment t-1; etalossRepresenting the rate of heat loss from the heat storage device;
Figure FDA00035295389900000216
represents the charging efficiency of the heat storage device;
Figure FDA00035295389900000217
representing the charging power of the heat storage device;
Figure FDA00035295389900000218
representing the heat release power of the heat storage device; Δ t represents a time interval;
the thermodynamic dynamic differential equation is not considered in the static model of the CSP power station, and the power generation power and the heat energy consumption power of the CSP power station can be regarded as a linear relation and are calculated as the formula (11):
Figure FDA0003529538990000031
wherein:
Figure FDA0003529538990000032
representing the generated power of the CSP power station for supplying load at the time t; etaRs,uThe positive rotation standby coefficient of the CSP power station;
Figure FDA0003529538990000033
providing positive rotation standby for the CSP power station at the time t; kappaCSPRepresenting the heat-electricity conversion efficiency of the CSP power station power generation system;
Figure FDA0003529538990000034
representing the heat energy consumed by the power generation system at time t;
3) establishment of day-ahead scheduling model considering PDR and CSP power station participating in wind power consumption
a) Objective function
The method comprises the following steps of comprehensively considering the power generation cost of a conventional thermal power generating unit, the operation and maintenance cost of a wind power and CSP power station, the renewable energy environmental benefit brought by grid-connected consumption, the cost of the thermal power generating unit and the CSP power station participating in load standby and wind power prediction error standby, establishing a combined scheduling model containing the wind power, the CSP power station and the thermal power station and based on optimal cost, and calculating a target function as a formula (12):
min F=D1+D2+D3+D4-D5 (12)
wherein: f is the comprehensive cost of the system when thermal power, wind power and CSP power stations participate in the optimized dispatching; d1The power generation cost when supplying load to the thermal power generating unit; d2The operation and maintenance cost of the wind power is reduced; d3The operation and maintenance cost for supplying power to the load by the CSP power station; d4The cost for the thermal power generating unit and the CSP power station to participate in load and wind power prediction error standby is saved; d5Environmental benefits brought to the power generation grid-connected consumption of wind power and CSP power stations;
the output of the thermal power generating unit is flexible and controllable, and the stable operation of a power grid can be ensured through reasonable optimized scheduling; in order to meet the load requirement in the scheduling process, output needs to be adjusted often, and even the unit needs to be scheduled for starting and stopping, and the power generation cost at this time includes the coal consumption cost of the thermal power unit and the start and stop cost of the unit, and is calculated as an equation (13) and an equation (14):
Figure FDA0003529538990000035
fi(PGi,t)=aiPGi,t 2+biPGi,t+ci (14)
wherein: t is the total duration; n is a radical ofGThe number of the thermal power generating units is; f. ofiThe coal consumption cost of the thermal power generating unit i is obtained; pGi,tThe output of the thermal power generating unit i in the time period t is obtained; v isi,tV and vi,t-1Respectively the running states of the thermal power generating unit i at t and t-1, if vi,tV denotes unit operation as 1i,tWhen the unit is stopped, 0 represents the unit is stopped; siThe start-stop cost of the unit i is calculated; a isi,bi,ciThe coal consumption cost parameter is the coal consumption cost parameter of the thermal power generating unit i;
wind power generation belongs to renewable energy power generation, coal is not consumed in the power generation process, but due to uncertainty of wind speed, operation and maintenance cost can be generated in the power generation process of a fan; the operation and maintenance cost of the wind power and the wind power output power can be approximately regarded as a linear relation, and are calculated as an expression (15):
Figure FDA0003529538990000036
wherein: kWOperating and maintaining costs for the wind farm; pW,tGenerating power of the wind power plant at the time t;
the CSP power station operation and maintenance cost is approximately regarded as a linear function of the generated power, and is calculated as the formula (16):
Figure FDA0003529538990000041
wherein: kCSPThe unit operation and maintenance cost of the CSP power station is saved;
in order to deal with uncertainty and emergency of load and wind power prediction, the reserve capacity of rotation needs to be reserved, and the reserve cost of a thermal power unit and a CSP power station is calculated as a formula (17):
Figure FDA0003529538990000042
wherein:
Figure FDA0003529538990000043
indicating the participation of the CSP power station in the load standby,
Figure FDA0003529538990000044
Representing the cost coefficient of the CSP power station participating in wind power standby;
Figure FDA0003529538990000045
indicating a positive spinning reserve provided by the CSP station for the load at time t,
Figure FDA0003529538990000046
Representing positive rotation standby power provided by the CSP power station for wind power at the time t;
Figure FDA0003529538990000047
representing the cost coefficient of the thermal power generating unit participating in the load standby,
Figure FDA0003529538990000048
Representing the cost coefficient of the thermal power generating unit participating in wind power standby;
Figure FDA0003529538990000049
representing positive rotation standby provided by the thermal power generating unit i for the load at the time t;
Figure FDA00035295389900000410
representing positive rotation standby power provided by the thermal power generating unit i for wind power at the time t;
wind power, CSP power station renewable energy are incorporated into the power networks and are consumed and have reduced thermal power unit's the net volume of generating electricity, effectively reduce the emission of sulfur and nitrate pollutant, can bring environmental benefit, calculate as (18) formula:
Figure FDA00035295389900000411
wherein: pW,tPower is consumed for the grid connection of the wind power plant at the time t; rhoCSPEnvironmental benefit coefficient brought to CSP power station grid-connected consumption; rhoWEnvironmental benefit coefficient brought to wind power grid connection consumption;
b) system constraints
When the network transmission loss is not counted, the sum of the output power of each unit is equal to the value after the power grid load response change, and is expressed as a formula (19):
Figure FDA00035295389900000412
wherein: l istThe load power value before the system needs to respond at the time t;
Figure FDA00035295389900000413
the load response rate triangular fuzzy number is converted into a deterministic variable, and then the expected value of the load response quantity at the time t is expressed as a formula (20):
Figure FDA00035295389900000414
the transmission capacity constraint of the transmission line is expressed as formula (21):
-Pij,max≤Biji,tj,t)≤Pij,max (21)
wherein: pij,maxIs the maximum transmission capacity of the transmission line between nodes i and j; b isijIs the susceptance between nodes ij; thetai,tIs the voltage phase angle of node i; thetaj,tIs the voltage phase angle of node j;
the opportunity constraint is adopted to determine the wind power reserve capacity, the wind power reserve cost is reduced on the basis of ensuring the safety of a power grid, and the concrete constraint is expressed as a formula (22):
Figure FDA0003529538990000051
wherein: cr { } is a confidence expression;
Figure FDA0003529538990000052
the positive rotation reserve capacity of the thermal power generating unit i at the moment t is obtained;
Figure FDA0003529538990000053
the negative rotation reserve capacity of the thermal power generating unit i at the moment t is obtained;
Figure FDA0003529538990000054
positive rotation reserve capacity provided for CSP power station at time t,
Figure FDA0003529538990000055
The negative spinning reserve capacity provided for the CSP plant at time t,
Figure FDA0003529538990000056
μLis the load reserve factor; alpha and beta are respectively confidence coefficients meeting the positive and negative rotation standby constraints;
the heat storage device of the CSP power station has to be within the rated limit range, and the heat storage device cannot be simultaneously charged/discharged in each scheduling period, and the specific constraint is expressed as the following formula (23) and formula (24):
Figure FDA0003529538990000057
Pt chaPt dis=0 (24)
wherein:
Figure FDA0003529538990000058
the maximum charging power of the heat storage system;
Figure FDA0003529538990000059
the maximum heat release power of the heat storage system;
the heat storage quantity of the CSP power station heat storage device is expressed as a formula (25):
Figure FDA00035295389900000510
wherein:
Figure FDA00035295389900000511
the minimum heat storage quantity of the heat storage device is obtained; xiTSThe maximum heat storage capacity of the heat storage device is expressed by the number of hours of full load;
Figure FDA00035295389900000512
the maximum output power of the CSP power station;
the specific constraint of the maximum/minimum output of the thermal power generating unit is expressed as a formula (26):
UGi,tPi,min≤PGi,t≤UGi,tPi,max (26)
wherein: u shapeGi,tThe method comprises the following steps that 1 represents operation, and 0 represents shutdown, wherein the thermal power generating unit i is in an operating state; pi,minThe minimum output of the thermal power generating unit i is obtained; pi,maxThe maximum output of the thermal power generating unit i is obtained;
the ramp rate constraint expression of the thermal power generating unit is as follows (27):
-rdi≤PGi,t-PGi,t-1≤rui (27)
wherein: r isuiThe upward climbing rate of the unit i is obtained; r isdiThe downward climbing speed of the unit i is obtained;
during day-ahead scheduling, the wind power on-line power cannot exceed the predicted value, and the specific expression is (28):
Figure FDA0003529538990000061
wherein:
Figure FDA0003529538990000062
the predicted output of the wind power plant at the time t is obtained;
the sum of the expected load response values after the demand response is zero in the whole scheduling period, and is expressed as formula (29):
Figure FDA0003529538990000063
considering the benefit of the user, the load variation needs to be limited, and the satisfaction degree of the power utilization mode and the satisfaction degree of the power utilization expense expenditure of the user are used as measurement indexes and are expressed as formulas (30) and (31):
Figure FDA0003529538990000064
Figure FDA0003529538990000065
wherein:
Figure FDA0003529538990000066
using electricity for userA minimum value of satisfaction;
Figure FDA0003529538990000067
the minimum value of the satisfaction degree of electricity cost expenditure;
Figure FDA0003529538990000068
the load value after the demand response at the time t is calculated as (32):
Figure FDA0003529538990000069
the load value after the price type demand response is between the upper limit and the lower limit of the load value before the response, and the specific constraint expression is (33):
Figure FDA00035295389900000610
wherein: l ismaxIs the maximum value of the original load curve, LminIs the minimum value of the original load curve; etamax、ηminTaking the value eta for the peak-to-valley difference coefficient of the demand response loadmax≥1,ηmin≤1。
CN201910172439.5A 2019-03-07 2019-03-07 Wind power consumption day-ahead scheduling method considering price type demand response and participation of CSP power station Active CN109858711B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910172439.5A CN109858711B (en) 2019-03-07 2019-03-07 Wind power consumption day-ahead scheduling method considering price type demand response and participation of CSP power station

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910172439.5A CN109858711B (en) 2019-03-07 2019-03-07 Wind power consumption day-ahead scheduling method considering price type demand response and participation of CSP power station

Publications (2)

Publication Number Publication Date
CN109858711A CN109858711A (en) 2019-06-07
CN109858711B true CN109858711B (en) 2022-05-24

Family

ID=66900166

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910172439.5A Active CN109858711B (en) 2019-03-07 2019-03-07 Wind power consumption day-ahead scheduling method considering price type demand response and participation of CSP power station

Country Status (1)

Country Link
CN (1) CN109858711B (en)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110224398B (en) * 2019-06-12 2023-04-18 广东电网有限责任公司 Method, device and system for clearing current market of electric power day by day and storage medium
CN110224399B (en) * 2019-06-12 2023-06-30 广东电网有限责任公司 Electric power spot market clearing method and device based on wind power plant and storage medium
CN110247408B (en) * 2019-06-28 2021-03-02 广东电网有限责任公司 Three-phase unbalance management method for transformer area based on demand response
CN110535185B (en) * 2019-08-09 2022-03-18 东北电力大学 Multi-source optimized operation method of photothermal power station with consideration of demand response
CN111047345B (en) * 2019-10-11 2022-04-12 国网甘肃省电力公司电力科学研究院 Electricity price data adjusting method considering new energy consumption rate and user satisfaction degree
CN110930073B (en) * 2019-12-14 2022-09-20 东北电力大学 Day-ahead scheduling method for wind-light-photo-thermal combined power generation system considering price type demand response
CN111277005B (en) * 2020-02-19 2022-03-18 东北电力大学 Multi-source power system multi-time scale scheduling method considering source-load coordination optimization
CN111740413B (en) * 2020-06-28 2021-11-30 东北电力大学 Optimal scheduling method of thermal power peak regulation initiative-considering energy storage-containing power system
CN111950807B (en) * 2020-08-26 2022-03-25 华北电力大学(保定) Comprehensive energy system optimization operation method considering uncertainty and demand response
CN113394817B (en) * 2021-06-28 2023-03-28 国网甘肃省电力公司电力科学研究院 Multi-energy capacity optimal configuration method of wind, light, water and fire storage system
CN113947427A (en) * 2021-10-14 2022-01-18 广西大学 Long-term price guiding method for dynamic differential control of flexible energy hybrid network

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105634024A (en) * 2016-02-23 2016-06-01 国网江苏省电力公司电力科学研究院 Price demand response-based intraday economic scheduling model and linear solving method
CN106655246A (en) * 2016-10-18 2017-05-10 国网黑龙江省电力有限公司哈尔滨供电公司 Method of solving robust two-layer optimization model based on wind power prediction and demand response
CN106651473A (en) * 2017-01-23 2017-05-10 国网福建省电力有限公司 Method for promoting wind power acceptance level by considering recent electricity price per hour and multiple demand responses
CN106961124A (en) * 2017-02-27 2017-07-18 南阳理工学院 A kind of economic load dispatching model responded based on price demand
CN107679776A (en) * 2017-11-14 2018-02-09 上海泛智能源装备有限公司 A kind of energy scheduling method and distributed energy resource system, distributed energy network system
CN108460482A (en) * 2018-01-31 2018-08-28 南京理工大学 Consider interruptible load and a few days ago the wind electricity digestion method of hour electricity price
CN108923472A (en) * 2018-08-27 2018-11-30 东北电力大学 Combine power output dispatching method with fired power generating unit based on the photo-thermal power station of Optimum cost

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105634024A (en) * 2016-02-23 2016-06-01 国网江苏省电力公司电力科学研究院 Price demand response-based intraday economic scheduling model and linear solving method
CN106655246A (en) * 2016-10-18 2017-05-10 国网黑龙江省电力有限公司哈尔滨供电公司 Method of solving robust two-layer optimization model based on wind power prediction and demand response
CN106651473A (en) * 2017-01-23 2017-05-10 国网福建省电力有限公司 Method for promoting wind power acceptance level by considering recent electricity price per hour and multiple demand responses
CN106961124A (en) * 2017-02-27 2017-07-18 南阳理工学院 A kind of economic load dispatching model responded based on price demand
CN107679776A (en) * 2017-11-14 2018-02-09 上海泛智能源装备有限公司 A kind of energy scheduling method and distributed energy resource system, distributed energy network system
CN108460482A (en) * 2018-01-31 2018-08-28 南京理工大学 Consider interruptible load and a few days ago the wind electricity digestion method of hour electricity price
CN108923472A (en) * 2018-08-27 2018-11-30 东北电力大学 Combine power output dispatching method with fired power generating unit based on the photo-thermal power station of Optimum cost

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
多种需求响应和日前小时电价优化促进风电接纳研究;张艺渊 等;《电工电能新技术》;20170721;第第37卷卷(第01期);57-65 *
考虑风电消纳的电热联合系统源荷协调优化调度;崔雪 等;《电力自动化设备》;20180703;第第38卷卷(第07期);74-81 *
计及需求响应的含风电电力系统日前经济调度;曹慧秋;《中国优秀硕士学位论文全文数据库工程科技II辑》;20190115(第12期);C042-969 *

Also Published As

Publication number Publication date
CN109858711A (en) 2019-06-07

Similar Documents

Publication Publication Date Title
CN109858711B (en) Wind power consumption day-ahead scheduling method considering price type demand response and participation of CSP power station
CN111738502B (en) Multi-energy complementary system demand response operation optimization method for promoting surplus wind power consumption
CN109523052B (en) Virtual power plant optimal scheduling method considering demand response and carbon transaction
CN109617142B (en) CCHP type micro-grid multi-time scale optimization scheduling method and system
CN112990523B (en) Hierarchical optimization operation method for regional comprehensive energy system
Guo et al. Energy storage technology and its typical application in new energy grid connection system
CN113541195B (en) Method for consuming high-proportion renewable energy in future power system
CN112418537B (en) Optimized scheduling method for multi-energy cloud energy storage system
CN110676846B (en) Power peak regulation method, dispatching center and power peak regulation system
Xu et al. Optimal intraday rolling operation strategy of integrated energy system with multi-storage
CN112531785A (en) Photo-thermal-heat-storage active power distribution network main and auxiliary combined scheduling method under multiple time scales
CN111882452A (en) Marginal cost calculation method for integrated energy system participating in demand side response
Zhai et al. Low carbon economic dispatch of regional integrated energy system considering load uncertainty
CN116805192A (en) Comprehensive energy system double-layer planning optimization method considering optimal energy rejection rate and application thereof
Hu et al. Low carbon economic dispatch of multi-energy combined system considering carbon trading
CN115912419A (en) Heat storage-electricity storage cooperative scheduling method based on load aggregation quotient
CN115936336A (en) Virtual power plant capacity configuration and regulation operation optimization method
Li et al. Multi-energy coordinated operation optimization model for wind-solar-hydro-thermal-energy storage system considering the complementary characteristics of different power resources
CN114759614A (en) Multi-energy system robust optimization scheduling strategy based on game theory
CN113627720A (en) Comprehensive park wind-solar energy storage system configuration method based on electric energy substitution
Mei et al. Multi-objective coordinated optimal scheduling of virtual power plants based on demand side response
Wang et al. Quantitative and apportionment method of auxiliary service cost of new energy subject
CN115563816B (en) Low-carbon-oriented photovoltaic and wind power generation grid connection and energy storage optimization method and device
Huan et al. Optimization of Integrated Energy Service Provider Considering Carbon Trading
Xiaofei et al. Research on the quantitative evaluation method of wind power and photovoltaic power curtailment in real-time market transactions

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