CN108446865A - Thermo-electrically based on interval method couples multipotency streaming system power methods of risk assessment - Google Patents

Thermo-electrically based on interval method couples multipotency streaming system power methods of risk assessment Download PDF

Info

Publication number
CN108446865A
CN108446865A CN201810340625.0A CN201810340625A CN108446865A CN 108446865 A CN108446865 A CN 108446865A CN 201810340625 A CN201810340625 A CN 201810340625A CN 108446865 A CN108446865 A CN 108446865A
Authority
CN
China
Prior art keywords
power
streaming system
multipotency streaming
electric
supply network
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
CN201810340625.0A
Other languages
Chinese (zh)
Other versions
CN108446865B (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.)
Guangzhou Power Supply Bureau of Guangdong Power Grid Co Ltd
Beijing King Star Hi Tech System Control Co Ltd
Original Assignee
Guangzhou Power Supply Bureau Co Ltd
Beijing King Star Hi Tech System Control Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Guangzhou Power Supply Bureau Co Ltd, Beijing King Star Hi Tech System Control Co Ltd filed Critical Guangzhou Power Supply Bureau Co Ltd
Priority to CN201810340625.0A priority Critical patent/CN108446865B/en
Publication of CN108446865A publication Critical patent/CN108446865A/en
Application granted granted Critical
Publication of CN108446865B publication Critical patent/CN108446865B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0635Risk analysis of enterprise or organisation activities
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06313Resource planning in a project environment
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply

Landscapes

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

Abstract

The present invention relates to a kind of coupled thermomechanics multipotency streaming system power methods of risk assessment based on interval method, belongs to the risk assessment technology field of integrated energy system.The methods of risk assessment of the present invention, establish the stochastic model in the device model of multipotency streaming system, multipotency streaming system, energy storage device in consideration system, flexible unit Optimized Operation probabilistic consumption is acted on, provide the interval method for solving this stochastic problems, calculate the uncertainty of regenerative resource in following one day, the uncertainty of electric user's heat user energy consumption can caused by multipotency streaming system power-balance risk size.This method more accurately calculates the operation risk of electro thermal coupling multipotency streaming system.Be conducive to improve safety in operation, and help to reduce the operating cost of multipotency streaming system.To explore the polynary consumption of regenerative resource, multifunctional system combined dispatching is made that certain contribution.

Description

Thermo-electrically based on interval method couples multipotency streaming system power methods of risk assessment
Technical field
The present invention relates to a kind of, and the thermo-electrically based on interval method couples multipotency streaming system power methods of risk assessment, belongs to The risk assessment technology field of integrated energy system.
Background technology:
Multipotency streaming system covers that electric, hot, cold, gas etc. is multiple to flow subsystem, and each subsystem passes through cogeneration (CHP/ CCHP), the equipment such as heat pump are converted and are coupled.The rise of energy internet in recent years so that with electric, hot, cold, gas is a variety of The multipotency network of form transmission energy shows its superiority.On the one hand the foundation of multipotency streaming system has considerable benefit, separately The coupling of one side multipotency also causes the uncertain factor in each subsystem to interact, and can influence the just common energy of user, very The overall security of multipotency streaming system can extremely be influenced.Multipotency streaming system risk assessment is exactly before risk case occurs, to possible Caused by influence and the possibility of loss is assessed, to carry out early warning to the traffic control personnel of multipotency streaming system, so as to Carry out corresponding control means.Multipotency streaming system risk assessment needs to consider a variety of couplings that can be between stream, needs to consider The uncertainty of regenerative resource and multipotency stream customer charge, compared to the risk assessment of conventional electric power system, multipotency stream risk Assessment has high complexity.
One multipotency streaming system generally comprises the renewable energy power generations such as wind-power electricity generation, photovoltaic generation, the work(of these power supplys Rate changes and is affected with the random variation of weather, this may result in grid power imbalance;Electric load, thermic load itself There is randomness, unbalanced power may also can be led to by being coupled to grid side;Energy storage device, flexible unit in power grid are quickly adjusted Output power can reduce the unbalanced probability of grid power.By risk assessment to above-mentioned factor carry out comprehensive analysis after It can obtain the size of multipotency streaming system operation risk.Traffic control personnel can be instructed to implement control hand when risk is larger Section improves safety in operation;Multipotency streaming system can be instructed to take relatively radical operation when risk is smaller or no risk Mode, to reduce the operating cost of multipotency streaming system.
Invention content
Purpose of the present invention is to propose a kind of thermo-electrically coupling multipotency streaming system power risk assessment side based on interval method Method, establishes the stochastic model in the device model of multipotency streaming system, multipotency streaming system, energy storage device in consideration system, flexibly The Optimized Operation of unit acts on probabilistic consumption, provides the interval method for solving this stochastic problems, then calculates The uncertainty of regenerative resource, the uncertainty of electric user's heat user energy consumption can be flat to multipotency streaming system power in one day following Risk size caused by weighing apparatus.
Thermo-electrically coupling multipotency streaming system power methods of risk assessment proposed by the present invention based on interval method, including with Lower step:
(1) an electric-thermal coupling multipotency streaming system operation stochastic model is established:
If in electric-thermal coupling multipotency streaming system including S cogeneration units, E electric energy storage device, H hot energy storage Equipment, L1 electric load, L2 thermic load and P photovoltaic plant, electric-thermal coupling multipotency streaming system pass through interconnection and external electrical Net is connected;
The active power of the photovoltaic generation of the one day future each scheduling instance t obtained from electric-thermal coupling multipotency streaming system is pre- The interval prediction data of measured value, electric load active power predicted value and thermic load active power predicted value, establish stochastic model It is as follows:
Ppvdown(t)≤Ppv(t)≤Ppvup(t)
Pedown(t)≤Pe(t)≤Peup(t)
Phdown1(t)≤Ph,1(t)≤Phup,1(t)
Phdown2(t)≤Ph,2(t)≤Phup,2(t)
...
Phdown,L2(t)≤Ph,L2(t)≤Phup,L2(t)
Wherein, t is scheduling instance, t=1,2 ..., T, for distinguishing the following intraday T scheduling instance, Ppv(t) it is Total generated output of all P photovoltaic plants of t scheduling instances, Pe(t) it is total wattful power of all L1 electric loads of t scheduling instances Rate, Ph,1(t)、Ph,2(t)、...、Ph,L2(t) power of each thermic load, P are indicated respectivelypvdown(t)、Pedown(t)、Phdown1(t)、 Phdown2(t)、...、Phdown,L2(t) be respectively t scheduling instances photovoltaic generation, electric load, thermic load power interval prediction under Limit, Ppvup(t)、Peup(t)、Phup,1(t)、Phup,2(t)、...、Phup,L2(t) be respectively t scheduling instances photovoltaic generation power, The upper limit that electric load power, thermic load power interval are predicted, Ppv(t)、Pe(t)、Ph,1(t)、Ph,2(t)、...、Ph,L2(t) be with Machine variable;
(2) constraint of electric-thermal coupling multipotency streaming system is established, including:
(2-1) energy storage device runs constraints:
The constraints of E electric energy storage device is:
0≤Pdis(t)≤Pdmax, 0≤Pchar(t)≤Pcmax
0≤SoC(t)≤SoCmax
SoC (t)=SoC (t-1)-Pdis(t)+Pchar(t), t=2,3 ..., T
Wherein, Pdis(t) it is total discharge power of all E electric energy storage devices of t-th of scheduling instance, Pchar(t) it is t Total charge power of the electric energy storage device of all E of a scheduling instance, PdmaxFor the maximum value of total discharge power, PcmaxAlways to charge The maximum value of power, SoC (t) are all E total storing up electricity states of electric energy storage device of t-th of scheduling slot, SoCmaxFor stored energy capacitance Maximum value.Pdis(t)、Pchar(t), SoC (t) is respectively the amount for being scheduled control;
The constraints of H hot energy storage devices is:
0≤Qdis,hh(t)≤Qdmax,hh, 0≤Qchar,hh(t)≤Qcmax,hh,0≤EHhh(t)≤EHmax,hh
EHhh(t)=EHhh(t-1)+Qchar,hh(t)-Qdis,hh(t)
Hh=1,2 ..., H
Wherein, Qdis,hh(t) it is the heat release power of each hot energy storage device of t-th of scheduling instance, Qdmax,hhIt is set for each hot energy storage Standby heat release power maximum value, Qchar,hh(t) it is each hot energy storage device Endothermic power of t-th of scheduling instance, Qcmax,hhFor each heat storage The Endothermic power maximum value of energy equipment, EHhh(t) it is the heat accumulation state of each hot energy storage device of t-th of scheduling instance, EHmax,hhIt is each The maximum value of hot energy storage device heat storage capacity, Qdis,hh(t)、Qchar,hh(t)、EHhh(t) it is respectively the amount for being scheduled control;
(2-2) cogeneration units run constraints:
By the feasible zone of s-th of cogeneration units in S cogeneration units, the D of convex polygonsA vertex is denoted as (Pd,s,Qd,s), d=1,2...Ds, feasible zone is described with the linear combination on vertex, then a bit (p in feasible zones,qs) full Foot:
Wherein, the electrical power p of each cogeneration units1(t)、p2(t)、...、pS(t), each cogeneration units Thermal power q1(t)、q2(t)、...、qS(t);
Cogeneration units active power climbing constraints be:
Wherein, RAMPs downAnd RAMPs upRespectively upward, the downward climbing speed of s platforms cogeneration units active power The maximum value of rate, ps(t) and ps(t-1) it is respectively s platforms cogeneration units when t-th of scheduling instance and the t-1 are dispatched The active power at quarter;
(2-3) obtains the operation constraints of power grid and heat supply network in electric-thermal coupling multipotency streaming system:
From electric-thermal couple multipotency streaming system Energy Management System in, obtain electric-thermal coupling multipotency streaming system in power grid, The safe operation of heat supply network constrains Ψsc, including:Grid branch active power allows maximum value and grid branch active power to allow most Small value, electric-thermal, which couples heat supply network node temperature in multipotency streaming system, allows maximum value, heat supply network node temperature to allow minimum value, heat supply network Bypass flow allows maximum value and heat supply network bypass flow to allow minimum value;
(2-4) heat supply network trend constraint condition is:
Assuming that the flow in heat-net-pipeline remains unchanged, the heat supply network power flow equation for establishing matrix form is as follows:
In above-mentioned heat supply network power flow equation, A is incidence matrix, remembers in heat supply network there is N number of node and B branch, is associated in foundation MatrixWith lower incidence matrixA
In above-mentioned heat supply network power flow equation:CpFor the specific heat capacity of heat supply network heat transfer mediator, m is a branch in B branch Flow, the subscript of m indicate that the serial number of branch, M are the vector of the flow composition of B branch, TnFor each node temperature composition of heat supply network Vector, TeFor the vector that the temperature of each branch end of heat supply network forms, TaFor environment temperature, QJFor the thermic load of heat supply network node The vector of composition, C are the diagonal matrixs using the pipeline heat loss coefficient of each branch as diagonal element, and λ is that the unit of each branch of heat supply network is long Thermal conductivity is spent, the subscript of λ indicates that the serial number of branch, L are the length of each branch of heat supply network, and the subscript of L indicates the serial number of branch;
(3) object function of thermo-electrically coupling multipotency streaming system power-balance risk is established:
When calculating the power least risk value of scheduling instance t, the object function OF of power-balance risklower(t) it is:
OFlower(t)=Ptie(t)
When calculating the power greatest risk value of some t moment, the object function OF of power-balance riskupper(t) it is
OFupper(t)=- Ptie(t)
Wherein, Ptie(t) it is the power on interconnection, the power on interconnection is equal to imbalance power;
(4) the operating cost object function OF of thermo-electrically coupling multipotency streaming system is establishedinner
Wherein, S is cogeneration units number, and s is the serial number of cogeneration units, Hprices(t) it is s-th of thermoelectricity The f timesharing prices for sending out unit thermal power of coproduction unit, Eprices(t) it is that s-th cogeneration units send out unit electricity The timesharing price of power, Ctie(t) it is
Ctie(t)=Price_n (t) * Ptie1(t)+Price_ex*Ptie2(t)
Ptie(t)=Ptie1(t)+Ptie2(t)
0≤Ptie1(t)≤Ptie_limit
0≤Ptie2(t)
Wherein CtieFor interconnection purchases strategies, Ptie1For permitted imbalance power threshold value, Ptie2For more than threshold value Power, Price_n (t) be imbalance power in threshold value when interconnection on tou power price, Price_ex be more than threshold value Partial punishment electricity price, Ptie_limitFor the upper limit of the power of interconnection carrying;
(5) using the object function of the thermo-electrically of above-mentioned steps (3) coupling multipotency streaming system power-balance risk as outer layer mesh Scalar functions, above-mentioned steps (1) are used as outer layer constraints, with the thermo-electricallies of above-mentioned steps (4) couple the operation of multipotency streaming system at It is internal layer constraints that this object function, which is used as internal layer object function, above-mentioned steps (2), builds a bi-level optimal model, and Using interior point method, 2T solution is carried out, OF is obtainedlower(t) and OFupper(t), t=1,2 ... T is to get to each scheduling The upper bound of the thermo-electrically coupling multipotency streaming system power risk at moment and lower bound, realize and couple multipotency streaming system power to thermo-electrically Risk assessment.
It is proposed by the present invention based on interval method thermo-electrically coupling multipotency streaming system methods of risk assessment, feature and Advantage is:
Methods of risk assessment proposed by the present invention is established random in the device model of multipotency streaming system, multipotency streaming system Property model, energy storage device in consideration system, flexible unit Optimized Operation probabilistic consumption is acted on, provide solution this Then the interval method of stochastic problems calculates the uncertainty of regenerative resource in one day following, electric user's heat user energy consumption Uncertainty can caused by multipotency streaming system power-balance risk size.Compared to traditional Study of Risk Evaluation Analysis for Power System, The method of the present invention considers electric-thermal coupled problem, it is also considered that the Optimized Operation of energy storage device, flexible unit is not in system Deterministic consumption effect, more accurately calculates the operation risk of electric-thermal coupling multipotency streaming system.The method of the present invention is conducive to Improve safety in operation, and contribute to reduce multipotency streaming system operating cost, for explore regenerative resource it is polynary consumption, it is more It system combined can dispatch and be made that certain contribution.
Description of the drawings
Fig. 1 is the schematic diagram for the electric-thermal coupling multipotency streaming system that the method for the present invention is related to.
Fig. 2 is feasible zone convex polygon schematic diagram in the cogeneration units operation constraint that the method for the present invention is related to.
Specific implementation mode
Thermo-electrically coupling multipotency streaming system power methods of risk assessment proposed by the present invention based on interval method, including with Lower step:
(1) an electric-thermal coupling multipotency streaming system operation stochastic model is established:
If in electric-thermal coupling multipotency streaming system including S cogeneration units, E electric energy storage device, H hot energy storage Equipment, L1 electric load, L2 thermic load and P photovoltaic plant, electric-thermal coupling multipotency streaming system pass through interconnection and external electrical Net is connected;
Electric-thermal involved in the method for the present invention couples multipotency streaming system as shown in Figure 1, being in multipotency stream garden in its center Portion, each component part can merge that be integrated into energy source station, renewable energy source station (photovoltaic plant), all kinds of cool and thermal powers negative in garden Lotus;The part that garden carries out Power Exchange with external electrical network is interconnection.Energy source station is interior including cogeneration units (CHP) and cold Thermoelectricity energy storage device.For a multipotency stream garden, the electricity of the cogeneration units in garden, which is contributed, to be tended not to meet garden Electric load, in order to meet energy balance require garden need by the Tie line Power with external electrical network.Thus multipotency The risk size of stream garden power-balance can be equivalent to the risk size of dominant eigenvalues overrate.
The active power of the photovoltaic generation of the one day future each scheduling instance t obtained from electric-thermal coupling multipotency streaming system is pre- The interval prediction data of measured value, electric load active power predicted value and thermic load active power predicted value, establish stochastic model It is as follows:
Ppvdown(t)≤Ppv(t)≤Ppvup(t)
Pedown(t)≤Pe(t)≤Peup(t)
Phdown1(t)≤Ph,1(t)≤Phup,1(t)
Phdown2(t)≤Ph,2(t)≤Phup,2(t)
...
Phdown,L2(t)≤Ph,L2(t)≤Phup,L2(t)
Wherein, t is scheduling instance, t=1,2 ..., T, for distinguishing the following intraday T scheduling instance, Ppv(t) it is Total generated output of all P photovoltaic plants of t scheduling instances, Pe(t) it is total wattful power of all L1 electric loads of t scheduling instances Rate, Ph,1(t)、Ph,2(t)、...、Ph,L2(t) power of each thermic load, P are indicated respectivelypvdown(t)、Pedown(t)、Phdown1(t)、 Phdown2(t)、...、Phdown,L2(t) be respectively t scheduling instances photovoltaic generation, electric load, thermic load power interval prediction under Limit, Ppvup(t)、Peup(t)、Phup,1(t)、Phup,2(t)、...、Phup,L2(t) be respectively t scheduling instances photovoltaic generation power, The upper limit that electric load power, thermic load power interval are predicted, Ppv(t)、Pe(t)、Ph,1(t)、Ph,2(t)、...、Ph,L2(t) be with Machine variable;
(2) constraint of electric-thermal coupling multipotency streaming system is established, including:
(2-1) energy storage device runs constraints:
The constraints of E electric energy storage device is:
0≤Pdis(t)≤Pdmax, 0≤Pchar(t)≤Pcmax
0≤SoC(t)≤SoCmax
SoC (t)=SoC (t-1)-Pdis(t)+Pchar(t), t=2,3 ..., T
Wherein, Pdis(t) it is total discharge power of all E electric energy storage devices of t-th of scheduling instance, Pchar(t) it is t Total charge power of the electric energy storage device of all E of a scheduling instance, PdmaxFor the maximum value of total discharge power, PcmaxAlways to charge The maximum value of power, SoC (t) are all E total storing up electricity states of electric energy storage device of t-th of scheduling slot, SoCmaxFor stored energy capacitance Maximum value.Pdis(t)、Pchar(t), SoC (t) is respectively the amount for being scheduled control;
The constraints of H hot energy storage devices is:
0≤Qdis,hh(t)≤Qdmax,hh, 0≤Qchar,hh(t)≤Qcmax,hh,0≤EHhh(t)≤EHmax,hh
EHhh(t)=EHhh(t-1)+Qchar,hh(t)-Qdis,hh(t)
Hh=1,2 ..., H
Wherein, Qdis,hh(t) it is the heat release power of each hot energy storage device of t-th of scheduling instance, Qdmax,hhIt is set for each hot energy storage Standby heat release power maximum value, Qchar,hh(t) it is each hot energy storage device Endothermic power of t-th of scheduling instance, Qcmax,hhFor each heat storage The Endothermic power maximum value of energy equipment, EHhh(t) it is the heat accumulation state of each hot energy storage device of t-th of scheduling instance, EHmax,hhIt is each The maximum value of hot energy storage device heat storage capacity, Qdis,hh(t)、Qchar,hh(t)、EHhh(t) it is respectively the amount for being scheduled control;
(2-2) cogeneration units run constraints:
There are certain relationship, the i.e. feasible zone of point (p, q) between the generated output p and heat supply output q of cogeneration units There are boundary, cogeneration units to be divided into back pressure type unit and extraction condensing type unit, the feasible zone of (p, q) of two kinds of units is presented Convex polygon, by s-th of cogeneration of heat and power machine in S cogeneration units, the cogeneration units that the method for the present invention is related to are transported Row constraint convex polygon schematic diagram is as shown in Figure 2, wherein Fig. 2 (a) is the feasible zone of back pressure type unit, and Fig. 2 (b) is extraction condensing type The feasible zone of unit, the D of convex polygonsA vertex is denoted as (Pd,s,Qd,s), d=1,2...Ds, retouched with the linear combination on vertex Feasible zone is stated, then a bit (p in feasible zones,qs) meet:
Wherein, the electrical power p of each cogeneration units1(t)、p2(t)、...、pS(t), each cogeneration units Thermal power q1(t)、q2(t)、...、qS(t);
Cogeneration units active power climbing constraints be:
Wherein, RAMPs downAnd RAMPs upRespectively upward, the downward climbing speed of s platforms cogeneration units active power The maximum value of rate, ps(t) and ps(t-1) it is respectively s platforms cogeneration units when t-th of scheduling instance and the t-1 are dispatched The active power at quarter;
(2-3) obtains the operation constraints of power grid and heat supply network in electric-thermal coupling multipotency streaming system:
From electric-thermal couple multipotency streaming system Energy Management System in, obtain electric-thermal coupling multipotency streaming system in power grid, The safe operation of heat supply network constrains Ψsc, including:Grid branch active power allows maximum value and grid branch active power to allow most Small value, electric-thermal, which couples heat supply network node temperature in multipotency streaming system, allows maximum value, heat supply network node temperature to allow minimum value, heat supply network Bypass flow allows maximum value and heat supply network bypass flow to allow minimum value;
(2-4) heat supply network trend constraint condition is:
Assuming that the flow in heat-net-pipeline remains unchanged, the heat supply network power flow equation for establishing matrix form is as follows:
In above-mentioned heat supply network power flow equation, A is incidence matrix, in order to describe the thermodynamic behaviour of heating network, needs pipeline Topological relation, remember in heat supply network there is a N number of node and B branch, incidence matrix in foundationWith lower incidence matrixA
In above-mentioned heat supply network power flow equation:CpFor the specific heat capacity of heat supply network heat transfer mediator, m is a branch in B branch Flow, the subscript of m indicate that the serial number of branch, M are the vector of the flow composition of B branch, TnFor each node temperature composition of heat supply network Vector, TeFor the vector that the temperature of each branch end of heat supply network forms, TaFor environment temperature, QJFor the thermic load of heat supply network node The vector of composition, C are the diagonal matrixs using the pipeline heat loss coefficient of each branch as diagonal element, and λ is that the unit of each branch of heat supply network is long Thermal conductivity is spent, the subscript of λ indicates that the serial number of branch, L are the length of each branch of heat supply network, and the subscript of L indicates the serial number of branch;
The vector M of above-mentioned each bypass flow composition of heat supply network, the vector T of each node temperature composition of heat supply networknWith each branch of heat supply network The vector T of last temperature compositioneFor the variable of control can be scheduled.
(3) object function of thermo-electrically coupling multipotency streaming system power-balance risk is established:
The part of thermo-electrically coupling multipotency streaming system unbalanced power is provided by interconnection, therefore the size of power risk can be with It is indicated with dominant eigenvalues.
When calculating the power least risk value of scheduling instance t, the object function OF of power-balance risklower(t) it is:
OFlower(t)=Ptie(t)
When calculating the power greatest risk value of some t moment, the object function OF of power-balance riskupper(t) it is
OFupper(t)=- Ptie(t)
Wherein, Ptie(t) it is the power on interconnection, the power on interconnection is equal to imbalance power;
(4) the operating cost object function OF of thermo-electrically coupling multipotency streaming system is establishedinner
It is all uncommon when thermo-electrically coupling multipotency streaming system does not have power-balance risk and when emergent power imbalance risk Hope operating cost as low as possible, it is assumed that unbalanced power is in certain threshold value on thermo-electrically coupling multipotency streaming system interconnection When be considered as there is no power-balance risk,
Wherein, S is cogeneration units number, and s is the serial number of cogeneration units, Hprices(t) it is s-th of thermoelectricity The f timesharing prices for sending out unit thermal power of coproduction unit, Eprices(t) it is that s-th cogeneration units send out unit electricity The timesharing price of power, Ctie(t) it is
Ctie(t)=Price_n (t) * Ptie1(t)+Price_ex*Ptie2(t)
Ptie(t)=Ptie1(t)+Ptie2(t)
0≤Ptie1(t)≤Ptie_limit
0≤Ptie2(t)
Wherein CtieFor interconnection purchases strategies, Ptie1For permitted imbalance power threshold value, Ptie2For more than threshold value Power, Price_n (t) be imbalance power in threshold value when interconnection on tou power price, Price_ex be more than threshold value Partial punishment electricity price, Ptie_limitFor the upper limit of the power of interconnection carrying;
(5) using the object function of the thermo-electrically of above-mentioned steps (3) coupling multipotency streaming system power-balance risk as outer layer mesh Scalar functions, above-mentioned steps (1) are used as outer layer constraints, with the thermo-electricallies of above-mentioned steps (4) couple the operation of multipotency streaming system at It is internal layer constraints that this object function, which is used as internal layer object function, above-mentioned steps (2), builds a bi-level optimal model, and Using interior point method, 2T solution is carried out, OF is obtainedlower(t) and OFupper(t), t=1,2 ... T is to get to each scheduling The upper bound of the thermo-electrically coupling multipotency streaming system power risk at moment and lower bound, realize and couple multipotency streaming system power to thermo-electrically Risk assessment.
For t moment, all it is unfavorable in the photovoltaic generation power at 96 moment of whole day, electric load power, thermic load power Under the scene for cutting down t moment dominant eigenvalues, (each scheduling mode all wraps the optimal scheduling mode for obtaining for t moment Operation plan containing 96 moment), the dominant eigenvalues under this scheduling mode of this scene are exactly the upper bound of dominant eigenvalues;It is similar , be all conducive to cut down t moment interconnection in the photovoltaic generation power at 96 moment of whole day, electric load power, thermic load power Under the scene of power, the optimal scheduling mode for t moment is obtained, the dominant eigenvalues under this scheduling mode of this scene are just It is the lower bound of dominant eigenvalues.The bound for having obtained dominant eigenvalues has also just obtained the thermo-electrically coupling multipotency stream system of t moment System power risk.

Claims (1)

1. a kind of thermo-electrically based on interval method couples multipotency streaming system power methods of risk assessment, it is characterised in that this method Include the following steps:
(1) an electric-thermal coupling multipotency streaming system operation stochastic model is established:
If include in electric-thermal coupling multipotency streaming system S cogeneration units, E electricity energy storage device, a hot energy storage devices of H, L1 electric load, L2 thermic load and P photovoltaic plant, electric-thermal coupling multipotency streaming system pass through interconnection and external electrical network phase Even;
The active power predicted value of the photovoltaic generation of the one day future each scheduling instance t obtained from electric-thermal coupling multipotency streaming system, The interval prediction data of electric load active power predicted value and thermic load active power predicted value, it is as follows to establish stochastic model:
Ppvdown(t)≤Ppv(t)≤Ppvup(t)
Pedown(t)≤Pe(t)≤Peup(t)
Phdown1(t)≤Ph,1(t)≤Phup,1(t)
Phdown2(t)≤Ph,2(t)≤Phup,2(t)
...
Phdown,L2(t)≤Ph,L2(t)≤Phup,L2(t)
Wherein, t is scheduling instance, t=1,2 ..., T, for distinguishing the following intraday T scheduling instance, Ppv(t) it is t tune Spend total generated output of moment all P photovoltaic plants, Pe(t) it is total active power of all L1 electric loads of t scheduling instances, Ph,1(t)、Ph,2(t)、...、Ph,L2(t) power of each thermic load, P are indicated respectivelypvdown(t)、Pedown(t)、Phdown1(t)、 Phdown2(t)、...、Phdown,L2(t) be respectively t scheduling instances photovoltaic generation, electric load, thermic load power interval prediction under Limit, Ppvup(t)、Peup(t)、Phup,1(t)、Phup,2(t)、...、Phup,L2(t) be respectively t scheduling instances photovoltaic generation power, The upper limit that electric load power, thermic load power interval are predicted, Ppv(t)、Pe(t)、Ph,1(t)、Ph,2(t)、...、Ph,L2(t) be with Machine variable;
(2) constraint of electric-thermal coupling multipotency streaming system is established, including:
(2-1) energy storage device runs constraints:
The constraints of E electric energy storage device is:
0≤Pdis(t)≤Pdmax, 0≤Pchar(t)≤Pcmax
0≤SoC(t)≤SoCmax
SoC (t)=SoC (t-1)-Pdis(t)+Pchar(t), t=2,3 ..., T
Wherein, Pdis(t) it is total discharge power of all E electric energy storage devices of t-th of scheduling instance, Pchar(t) it is t-th of scheduling Total charge power of all E of moment electric energy storage device, PdmaxFor the maximum value of total discharge power, PcmaxFor total charge power Maximum value, SoC (t) are all E total storing up electricity states of electric energy storage device of t-th of scheduling slot, SoCmaxFor the maximum of stored energy capacitance Value.Pdis(t)、Pchar(t), SoC (t) is respectively the amount for being scheduled control;
The constraints of H hot energy storage devices is:
0≤Qdis,hh(t)≤Qdmax,hh, 0≤Qchar,hh(t)≤Qcmax,hh,0≤EHhh(t)≤EHmax,hh
EHhh(t)=EHhh(t-1)+Qchar,hh(t)-Qdis,hh(t)
Hh=1,2 ..., H
Wherein, Qdis,hh(t) it is the heat release power of each hot energy storage device of t-th of scheduling instance, Qdmax,hhFor each hot energy storage device Heat release power maximum value, Qchar,hh(t) it is each hot energy storage device Endothermic power of t-th of scheduling instance, Qcmax,hhIt is set for each hot energy storage Standby Endothermic power maximum value, EHhh(t) it is the heat accumulation state of each hot energy storage device of t-th of scheduling instance, EHmax,hhFor each heat storage The maximum value of energy equipment heat storage capacity, Qdis,hh(t)、Qchar,hh(t)、EHhh(t) it is respectively the amount for being scheduled control;
(2-2) cogeneration units run constraints:
By the feasible zone of s-th of cogeneration units in S cogeneration units, the D of convex polygonsA vertex is denoted as (Pd,s, Qd,s), d=1,2...Ds, feasible zone is described with the linear combination on vertex, then a bit (p in feasible zones,qs) meet:
Wherein, the electrical power p of each cogeneration units1(t)、p2(t)、...、pS(t), the hot merit of each cogeneration units Rate q1(t)、q2(t)、...、qS(t);
Cogeneration units active power climbing constraints be:
Wherein, RAMPs downAnd RAMPs upRespectively upward, the downward creep speed of s platforms cogeneration units active power Maximum value, ps(t) and ps(t-1) it is respectively s platforms cogeneration units in t-th of scheduling instance and the t-1 scheduling instance Active power;
(2-3) obtains the operation constraints of power grid and heat supply network in electric-thermal coupling multipotency streaming system:
From the Energy Management System that electric-thermal couples multipotency streaming system, power grid, heat supply network in electric-thermal coupling multipotency streaming system are obtained Safe operation constrain Ψsc, including:It is minimum that grid branch active power allows maximum value and grid branch active power to allow Value, electric-thermal, which couples heat supply network node temperature in multipotency streaming system, allows maximum value, heat supply network node temperature to allow minimum value, heat supply network branch Road flow allows maximum value and heat supply network bypass flow to allow minimum value;
(2-4) heat supply network trend constraint condition is:
Assuming that the flow in heat-net-pipeline remains unchanged, the heat supply network power flow equation for establishing matrix form is as follows:
In above-mentioned heat supply network power flow equation, A is incidence matrix, remembers there is N number of node and B branch in heat supply network, incidence matrix in foundation With lower incidence matrixA
In above-mentioned heat supply network power flow equation:CpFor the specific heat capacity of heat supply network heat transfer mediator, m is the flow of a branch in B branch, m Subscript indicate the serial number of branch, M is the vector of the flow composition of B branch, TnFor the vector of heat supply network each node temperature composition, TeFor the vector that the temperature of each branch end of heat supply network forms, TaFor environment temperature, QJFor heat supply network node thermic load form to Amount, C is the diagonal matrix using the pipeline heat loss coefficient of each branch as diagonal element, and λ is the unit length thermal conductivity of each branch of heat supply network, The subscript of λ indicates that the serial number of branch, L are the length of each branch of heat supply network, and the subscript of L indicates the serial number of branch;
(3) object function of thermo-electrically coupling multipotency streaming system power-balance risk is established:
When calculating the power least risk value of scheduling instance t, the object function OF of power-balance risklower(t) it is:
OFlower(t)=Ptie(t)
When calculating the power greatest risk value of some t moment, the object function OF of power-balance riskupper(t) it is
OFupper(t)=- Ptie(t)
Wherein, Ptie(t) it is the power on interconnection, the power on interconnection is equal to imbalance power;
(4) the operating cost object function OF of thermo-electrically coupling multipotency streaming system is establishedinner
Wherein, S is cogeneration units number, and s is the serial number of cogeneration units, Hprices(t) it is s-th of cogeneration of heat and power machine The f timesharing prices for sending out unit thermal power of group, Eprices(t) it is that s-th cogeneration units send out specific electric power Timesharing price, Ctie(t) it is
Ctie(t)=Price_n (t) * Ptie1(t)+Price_ex*Ptie2(t)
Ptie(t)=Ptie1(t)+Ptie2(t)
0≤Ptie1(t)≤Ptie_limit
0≤Ptie2(t)
Wherein CtieFor interconnection purchases strategies, Ptie1For permitted imbalance power threshold value, Ptie2For the power more than threshold value, Price_n (t) be imbalance power in threshold value when interconnection on tou power price, Price_ex be more than threshold portion Punish electricity price, Ptie_limitFor the upper limit of the power of interconnection carrying;
(5) using the object function of the thermo-electrically of above-mentioned steps (3) coupling multipotency streaming system power-balance risk as outer layer target letter Number, above-mentioned steps (1) are used as outer layer constraints, and the operating cost mesh of multipotency streaming system is coupled with the thermo-electrically of above-mentioned steps (4) For scalar functions as internal layer object function, above-mentioned steps (2) are internal layer constraints, build a bi-level optimal model, and use Interior point method carries out 2T solution, obtains OFlower(t) and OFupper(t), t=1,2 ... T is to get to each scheduling instance Thermo-electrically coupling multipotency streaming system power risk the upper bound and lower bound, realize to thermo-electrically couple multipotency streaming system power wind Danger assessment.
CN201810340625.0A 2018-04-17 2018-04-17 Thermo-electrically based on interval method couples multipotency streaming system power methods of risk assessment Active CN108446865B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810340625.0A CN108446865B (en) 2018-04-17 2018-04-17 Thermo-electrically based on interval method couples multipotency streaming system power methods of risk assessment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810340625.0A CN108446865B (en) 2018-04-17 2018-04-17 Thermo-electrically based on interval method couples multipotency streaming system power methods of risk assessment

Publications (2)

Publication Number Publication Date
CN108446865A true CN108446865A (en) 2018-08-24
CN108446865B CN108446865B (en) 2019-08-16

Family

ID=63200695

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810340625.0A Active CN108446865B (en) 2018-04-17 2018-04-17 Thermo-electrically based on interval method couples multipotency streaming system power methods of risk assessment

Country Status (1)

Country Link
CN (1) CN108446865B (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109494816A (en) * 2018-12-28 2019-03-19 清华-伯克利深圳学院筹备办公室 The methods of risk assessment and device of the multipotency streaming system of electric-thermal coupling
CN109755949A (en) * 2019-01-07 2019-05-14 中国电力科学研究院有限公司 A kind of optimizing distribution method and device of combined heat and power energy-accumulating power station power
CN109816201A (en) * 2018-12-19 2019-05-28 中国农业大学 A kind of characterizing method of polyisocyanate mass-energy source homogeneity coupling
CN109991494A (en) * 2019-03-28 2019-07-09 国网山东省电力公司经济技术研究院 A kind of electric-thermal coupled system fault early warning method and prior-warning device
CN109993445A (en) * 2019-04-04 2019-07-09 国家电网有限公司 A kind of integrated energy system vulnerability assessment method considering photovoltaic prediction error
CN111047069A (en) * 2019-10-23 2020-04-21 浙江大学 Thermoelectric combined scheduling flexibility assessment method considering energy storage characteristics of heat supply network
CN111898806A (en) * 2020-07-14 2020-11-06 国网综合能源服务集团有限公司 Electric-thermal coupling source storage and load integration multi-energy flow park operation optimization method and system
WO2021159893A1 (en) * 2020-02-13 2021-08-19 清华大学 Optimal electricity-heat multi-energy flow system scheduling method based on heat supply phasor model

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140203639A1 (en) * 2013-01-21 2014-07-24 Hamilton Sundstrand Corporation Reconfigurable matrix-based power distribution architecture
CN105958480A (en) * 2016-05-27 2016-09-21 清华大学 Combined static safety analysis method for electrical-gas coupling multiple energy flow system
CN106067677A (en) * 2016-05-27 2016-11-02 清华大学 A kind of electro thermal coupling multipotency streaming system static security analysis method
CN106339794A (en) * 2016-08-16 2017-01-18 清华大学 Electric-thermal coupling multi-energy flow network node energy price calculation method
CN107590579A (en) * 2017-07-18 2018-01-16 燕山大学 A kind of heat and power system Joint economics dispatching method using heat supply network energy storage characteristic
CN107800158A (en) * 2017-10-17 2018-03-13 清华大学 A kind of electro thermal coupling multipotency streaming system Optimization Scheduling for taking into account economy and efficiency

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140203639A1 (en) * 2013-01-21 2014-07-24 Hamilton Sundstrand Corporation Reconfigurable matrix-based power distribution architecture
CN105958480A (en) * 2016-05-27 2016-09-21 清华大学 Combined static safety analysis method for electrical-gas coupling multiple energy flow system
CN106067677A (en) * 2016-05-27 2016-11-02 清华大学 A kind of electro thermal coupling multipotency streaming system static security analysis method
CN106339794A (en) * 2016-08-16 2017-01-18 清华大学 Electric-thermal coupling multi-energy flow network node energy price calculation method
CN107590579A (en) * 2017-07-18 2018-01-16 燕山大学 A kind of heat and power system Joint economics dispatching method using heat supply network energy storage characteristic
CN107800158A (en) * 2017-10-17 2018-03-13 清华大学 A kind of electro thermal coupling multipotency streaming system Optimization Scheduling for taking into account economy and efficiency

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109816201A (en) * 2018-12-19 2019-05-28 中国农业大学 A kind of characterizing method of polyisocyanate mass-energy source homogeneity coupling
CN109494816B (en) * 2018-12-28 2020-07-31 清华-伯克利深圳学院筹备办公室 Risk assessment method and device for electric-thermal coupling multi-energy flow system
CN109494816A (en) * 2018-12-28 2019-03-19 清华-伯克利深圳学院筹备办公室 The methods of risk assessment and device of the multipotency streaming system of electric-thermal coupling
CN109755949B (en) * 2019-01-07 2021-12-17 中国电力科学研究院有限公司 Optimal distribution method and device for power of combined heat and power energy storage power station
CN109755949A (en) * 2019-01-07 2019-05-14 中国电力科学研究院有限公司 A kind of optimizing distribution method and device of combined heat and power energy-accumulating power station power
CN109991494A (en) * 2019-03-28 2019-07-09 国网山东省电力公司经济技术研究院 A kind of electric-thermal coupled system fault early warning method and prior-warning device
CN109991494B (en) * 2019-03-28 2022-05-27 国网山东省电力公司经济技术研究院 Electric-thermal coupling system fault early warning method and early warning device
CN109993445A (en) * 2019-04-04 2019-07-09 国家电网有限公司 A kind of integrated energy system vulnerability assessment method considering photovoltaic prediction error
CN111047069A (en) * 2019-10-23 2020-04-21 浙江大学 Thermoelectric combined scheduling flexibility assessment method considering energy storage characteristics of heat supply network
CN111047069B (en) * 2019-10-23 2022-08-09 浙江大学 Thermoelectric combined scheduling flexibility assessment method considering energy storage characteristics of heat supply network
WO2021159893A1 (en) * 2020-02-13 2021-08-19 清华大学 Optimal electricity-heat multi-energy flow system scheduling method based on heat supply phasor model
CN111898806A (en) * 2020-07-14 2020-11-06 国网综合能源服务集团有限公司 Electric-thermal coupling source storage and load integration multi-energy flow park operation optimization method and system
CN111898806B (en) * 2020-07-14 2023-10-20 国网综合能源服务集团有限公司 Electric-thermal coupling source storage integrated multipotency flow park operation optimization method and system

Also Published As

Publication number Publication date
CN108446865B (en) 2019-08-16

Similar Documents

Publication Publication Date Title
CN108446865B (en) Thermo-electrically based on interval method couples multipotency streaming system power methods of risk assessment
CN107665377A (en) A kind of multiple source-coupled integrated energy system planing method
CN107832979A (en) A kind of factory integration energy resource system economic optimization dispatching method for considering cascaded utilization of energy
CN109784573A (en) A kind of energy internet Multipurpose Optimal Method and device
Chen et al. Optimal coordinative operation strategy of the electric–thermal–gas integrated energy system considering CSP plant
CN109034508A (en) Consider the dual probabilistic integrated energy system robust Optimization Scheduling of electric heating
CN109345012A (en) Garden energy internet running optimizatin method based on comprehensive evaluation index
CN109993345B (en) Garden-oriented dynamic economic dispatching method for multi-energy complementary system for island operation
CN105373842A (en) Micro-grid energy optimization and evaluation method based on full energy flow model
CN108596442A (en) The integrated energy system economic load dispatching method of meter and Conditional Lyapunov ExponentP
CN108039722A (en) A kind of distribution type renewable energy system Optimal Configuration Method suitable for alternating current-direct current mixing
CN111860966B (en) Energy storage-containing comprehensive energy system scheduling method based on fuzzy correlation opportunity planning
CN113036819A (en) Electric heating comprehensive energy system robust scheduling method considering source-load bilateral uncertainty
CN110443467A (en) A kind of regional complex energy resource system solar energy digestion capability appraisal procedure
CN105446163A (en) Micro-grid versatile flow modeling approach based on different types of energy
CN106056246A (en) Energy storage capacity optimization method of cooling heating and power multi-potent flow microgrid considering operation
CN105186583B (en) Energy router and its energy dispatching method based on multi-agent modeling
CN109816201A (en) A kind of characterizing method of polyisocyanate mass-energy source homogeneity coupling
CN112036652A (en) Photovoltaic-energy storage integrated energy system planning method based on opportunity constraint planning
CN112886571B (en) Decomposition, coordination and optimization operation method and device of electric heating comprehensive energy system based on boundary variable feasible region
CN109494816A (en) The methods of risk assessment and device of the multipotency streaming system of electric-thermal coupling
CN109768567A (en) A kind of Optimization Scheduling coupling multi-energy complementation system
CN112084633B (en) Electric heating coupling system optimal scheduling method considering dynamic characteristics of steam network
CN111598320B (en) Micro-energy system scheduling method
CN115081867A (en) Energy storage planning method, device and equipment for comprehensive energy system

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20231212

Address after: 510620, No. two, No. 2, Tianhe South Road, Guangzhou, Guangdong, Tianhe District

Patentee after: Guangzhou Power Supply Bureau of Guangdong Power Grid Co.,Ltd.

Address before: 815, Unit 1, 8th Floor, Building 1, No. 30 Yuzhi East Road, Huilongguan Town, Changping District, Beijing, 100085

Patentee before: BEIJING KING STAR HI-TECH SYSTEM CONTROL Co.,Ltd.

Patentee before: Guangzhou Power Supply Bureau of Guangdong Power Grid Co.,Ltd.

Effective date of registration: 20231212

Address after: 815, Unit 1, 8th Floor, Building 1, No. 30 Yuzhi East Road, Huilongguan Town, Changping District, Beijing, 100085

Patentee after: BEIJING KING STAR HI-TECH SYSTEM CONTROL Co.,Ltd.

Patentee after: Guangzhou Power Supply Bureau of Guangdong Power Grid Co.,Ltd.

Address before: 1201, Block D, Yousheng Building, No. 28 Chengfu Road, Haidian District, Beijing, 100083

Patentee before: BEIJING KING STAR HI-TECH SYSTEM CONTROL Co.,Ltd.

Patentee before: GUANGZHOU POWER SUPPLY Co.,Ltd.