CN108984894A - A kind of EHS pneumoelectric based on CCP purchases optimization method - Google Patents
A kind of EHS pneumoelectric based on CCP purchases optimization method Download PDFInfo
- Publication number
- CN108984894A CN108984894A CN201810747877.5A CN201810747877A CN108984894A CN 108984894 A CN108984894 A CN 108984894A CN 201810747877 A CN201810747877 A CN 201810747877A CN 108984894 A CN108984894 A CN 108984894A
- Authority
- CN
- China
- Prior art keywords
- energy
- power
- ehs
- pneumoelectric
- formula
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/06—Electricity, gas or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
Abstract
The present invention relates to a kind of, and the EHS pneumoelectric based on CCP purchases optimization method.Present invention aims at making up under renewable energy access, single energy hub is difficult to take into account the status of power supply reliability and economy problems, proposes that pneumoelectric purchases model outside the energy cluster system based on chance constrained programming.For the uncertainty of wind-powered electricity generation in renewable energy and photovoltaic, multimode modeling is carried out to stochastic variable using wind-powered electricity generation-photovoltaic probability density function, connecting each other between multiple energy hubs is established by Coupling method, by one to more, form energy cluster system, and constraint planning of improving the occasion is to power tie line and natural gas line power transmission upper limit proper treatment, both the randomness of renewable energy energy supply power had been fully considered, energy supply reliability, safety and economy are taken into account again, develop multi-energy system for economy of large scale and advantageous reference is provided.
Description
Technical field
The present invention relates to power system automatic fields, and in particular to a kind of EHS pneumoelectric based on CCP purchases optimization side
Method.
Background technique
With penetrating deep into for global energy, using energy source develops towards multipotency coordination, the direction provided multiple forms of energy to complement each other, and the energy is mutual
Networking becomes the integrated developing direction of various energy resources net.Traditional electric power, heating power and mutually independent operation mould of natural gas network
Formula can not adapt to current production of energy and Land use systems.Various energy resources are supplied and use energy demand high abstraction by energy hub
And classify, realize the function of energy conversion, storage, distribution.General cell, hospital, school etc. can be abstracted as single energy collection
Line device, but single energy hub model flexibility is weak, the schedulable energy is few, is unfavorable for the shared complementation of the energy.Meanwhile by can
Inhomogeneities, energy line concentration is presented in the limitation of randomness and fluctuation possessed by regenerated resources itself, renewable energy access
Device often needs to buy a certain amount of electric energy and natural gas to external electrical network, gas online shopping to meet electrothermal load demand, and it is bought
Amount is often limited to the power transmission ceiling restriction of power tie line between system and external electrical network, gas net, natural gas line.Cause
This, when energy hub workload demand can not be satisfied, system is to keep energy supply reliability, it will usually which investment construction power passes
The higher power tie line road of defeated upper limit value or natural gas line meet workload demand, so increase system investments construction at
This.Further, since lacking the clear model for considering multiple energy hub interconnections, existing research is mostly with single energy hub
Basis, such as consider that installing relevant device additional reaches energy supply requirement only for single energy hub.Under background above, existing skill
Art and model are difficult to take into account the economy of system in the case where guarantee system energizes reliability, meets workload demand.
Summary of the invention
Technical problem to be solved by the invention is to provide a kind of, and the EHS pneumoelectric based on CCP purchases optimization method, and it is difficult to solve
To take into account the economy problems of system in the case where guarantee system energizes reliability, meets workload demand.
The technical scheme to solve the above technical problems is that a kind of EHS pneumoelectric based on CCP purchases optimization side
Method, comprising the following steps:
S1, the mathematical model for establishing EH energy conversion devices and EH energy storage device, and show that EH's outputs and inputs relationship；
S2, EH renewable energy power output is calculated；
S3, by the mathematical model of EH energy conversion devices and EH energy storage device, EH output and input relationship and EH can
Renewable sources of energy power output establishes the EHS pneumoelectric based on CCP and purchases Optimized model, and purchases optimization mould using the EHS pneumoelectric based on CCP
Type is purchased EHS pneumoelectric and is optimized.
Based on the above technical solution, the present invention can also be improved as follows.
Further, the step S1 specifically includes the following steps:
S11, the mathematical model for establishing EH energy conversion devices:
Calculate the power supply and heat supply power output of CHP unit, calculation formula are as follows:
In formula (1) and formula (2),WithRespectively the t moment power generation in the d days of CHP unit and heating function
Rate,For the amount of natural gas of unit time CHP unit consumption, ηCHP,eleAnd ηCHP,thThe respectively gentle heat of pneumoelectric transformation efficiency
Transformation efficiency；
Calculate the heats power of GF, calculation formula are as follows:
In formula (3),WithThe respectively amount of natural gas of GF heats power and unit time consumption, ηGFFor GF
Gas thermal transition efficiency；
Calculate the heat production power of GSHP, calculation formula are as follows:
In formula (4),WithThe electric energy of the heat production power of respectively GSHP and unit time consumption,
For the performance parameter of GSHP；
S12, the mathematical model for establishing EH energy storage device:
Charging and discharging state limitation are as follows:
The charge and discharge energy work rate limit are as follows:
Capacity limit are as follows:
In above formula, x is energy, indicates electricity, thermal energy with es and hs respectively,WithFor energy storage device
Energy state and exoergic state are filled, is the variable of 0-1,WithRespectively energy-storage system
Charge and discharge energy work rate bound,WithEnergy exoergic power is respectively filled,For energy-storage system energy storage capacity, δxFor
Energy-storage system rate of energy loss, ηx,charge、ηx,dischargeIt respectively fills, exergic efficiency,WithRespectively minimum energy storage
Amount and maximum energy storage capacity,WithThe respectively energy storage capacity of d days finish times and initial time, Δ t are time step
1h；
S13, outputting and inputting for single EH is obtained by the mathematical model of energy conversion devices each in EH and energy storage device
Relationship are as follows:
In formula (11),For electrical energy demands,For thermal demand, ηTransFor transformer efficiency,For from electricity
Online shopping electricity,The electricity generated for honourable resource through blower and photovoltaic cell group.
Further, the step S2 specifically includes the following steps:
S21, photovoltaic generating system active power output P is calculatedPV, calculation formula are as follows:
In formula (13), r is radiancy, and A is the gross area of solar cell module, M, AmRespectively solar battery
The area of the battery pack number of packages of square matrix, m-th battery component, η, ηmRespectively total photoelectric conversion effect of solar cell module
Rate, m-th of photoelectric conversion efficiency, the probability density function of r are as follows:
In formula (14), Γ () is gamma function, rmaxFor greatest irradiation degree, α, β are Bata profile shape parameter；
S22, wind power output P is calculatedWT, calculation formula are as follows:
In formula (15), v is the wind speed of blower, vc、vfAnd vsWind speed, cut-out wind speed and rated wind speed are respectively cut,
PRRated power a, b for blower are constant, the probability density function of v are as follows:
In formula (16), K is form parameter, reflects the characteristic distributions of wind speed v, and C is scale parameter, and reflection this area is flat
The size of equal wind speed；
S23, renewable energy power output P is calculatedre, calculation formula are as follows:
Pre=PPV+PWT (17)。
Further, the specific steps of the step S3 are as follows:
S31, objective function is established:
TC=min (CGas+CGrid) (18)
In formula (18), TCFor the totle drilling cost that pneumoelectric is purchased, CGasIt always purchases gas expense for EHS in dispatching cycle to use, CGridFor
Total power purchase expense, wherein
In above formula, n be EHS in EH number, T be one day 24 period, t=1,2, T, D be total activation week
Phase, i=1,2,3, it is energy hub number in EHS,For power purchase price,For the unit time when EHS electric power is insufficient
When, from power grid purchase of electricity,For when EHS power surplus, to power grid electricity sales amount,For sale of electricity price；
S32, setting bound for objective function, constraint condition includes formula (1)-formula (12), additionally include with
Lower constraint:
Electric equilibrium constraint:
Thermal balance constraint:
Amount of natural gas Constraints of Equilibrium:
0-1 constraint:
CHP unit electric heating units limits:
The hot units limits of GF unit:
The hot units limits of GSHP unit:
The pneumoelectric amount of purchasing chance constraint:
In above formula,WithElectricity is injected when electric power deficiency between respectively i-th of EH and adjacent EH or major network
Electricity is sent outside when with power surplus,Qigong rate is purchased for the EHS unit time,WithFor 0-1 decision variable,For
Power purchase power maximum value,Absolute value be sale of electricity power maximum value, Pi maxFor i-th of EH injecting electric power maximum value, Pi min
Absolute value be send electrical power maximum value outside,AndRespectively CHP unit electricity goes out
The upper lower limit value of power and heat power output,WithRespectively GF unit heat power output upper lower limit value,WithRespectively
For the maximin of GSHP unit heat power output, Pr { } is the probability that event is set up,It is purchased in qigong rate for the EHS unit time
Limit value, αgridAnd αgasRespectively correspond to the confidence level of power tie line, natural gas line chance constrained programming condition；
S33, the pneumoelectric amount of purchasing chance constraint is converted to certainty constraint:
In above formula,Maximum transmission power is designed for power tie line,Maximum transmitted is designed for natural gas line
Power,WithRespectively 0-1 decision variable；
S34, Optimized model is purchased using the EHS pneumoelectric based on CCP EHS pneumoelectric is purchased and optimize.
The beneficial effects of the present invention are: present invention aims at making up under renewable energy access, single energy hub
It is difficult to take into account the status of power supply reliability and economy problems, propose outside the energy cluster system based on chance constrained programming
Pneumoelectric purchases model.For the uncertainty of wind-powered electricity generation in renewable energy and photovoltaic, wind-powered electricity generation-photovoltaic probability density letter is utilized
Several pairs of stochastic variables carry out multimode modeling, connecting each other between multiple energy hubs are established by Coupling method, by one
To more, formation energy cluster system, and constraint planning of improving the occasion is on power tie line and natural gas line power transmission
Proper treatment is limited, had not only fully considered the randomness of renewable energy energy supply power, but also takes into account energy supply reliability, safety and warp
Ji property develops multi-energy system for economy of large scale and provides advantageous reference.
Detailed description of the invention
Fig. 1 is the total flow chart of steps of the present invention；
Fig. 2 is single energy hub schematic diagram in the embodiment of the present invention；
Fig. 3 is energy cluster system schematic diagram in the embodiment of the present invention；
Fig. 4, which is that situation 1 is lower in the embodiment of the present invention, purchases gas power purchase power CDF curve；
Fig. 5, which is that situation 2 is lower in the embodiment of the present invention, purchases gas power purchase power CDF curve.
Specific embodiment
The principle and features of the present invention will be described below with reference to the accompanying drawings, and the given examples are served only to explain the present invention, and
It is non-to be used to limit the scope of the invention.
Energy hub (energy hub, EH) be defined as in a kind of description multi-energy system the energy, load, network it
Between exchange, coupled relation input-output port model.EH can play the multiple types energy and mutually convert, distributes and store up
The effect deposited has very strong flexibility.EH can adjust the power output and operation of distinct device according to model optimization target difference
State, and then meet the needs of corresponding load.Have studied the EH based on cogeneration herein, meter and three kinds of electricity, air and heat energy
Form, single EH using chance constrained programming (chance constrained programming, CCP) as shown in Fig. 2, built
Mould can describe to improve robustness and reliability by the not true property of stochastic variable bring well.
Energy conversion devices specifically include that photovoltaic cell (photovoltaic, PV), blower (wind in single EH
Turbine, WT), cogeneration units (combined heat and power, CHP), earth source heat pump (ground source
Heat pump, GSHP), gas fired-boiler (gas furnace, GF), power transformer (transformer, Trans)；Meanwhile
Using storage (electric storage, ES), heat accumulation (heat storage, HS) system can realize the energy storage and it is more when
The transfer of section.As shown in Figure 2, EH is the multiport network of a broad sense, and port, which is divided into, outputs and inputs 2 kinds.Input side packet
Containing the electricity P bought through power tie line from power gridgridEH is flowed by transformer, honourable resource is produced through blower and photovoltaic cell group
Raw electricity Pre, the amount of natural gas P that enters through natural gas line from gas online shoppinggas, can be by a m dimensional vector VinIt indicates, it is similar
Seemingly, outlet side includes electrical energy demands LeleWith thermal demand Lth, with a n-dimensional vector LoutIt indicates.Coupling matrix C defines Vin
And LoutTransformational relation within a certain period of time.It can thus be concluded that the input of single EH (uses VinIndicate) it exports and (uses LoutIndicate) relationship
Are as follows:
Lout=CVin
As shown in Figure 1, a kind of EHS pneumoelectric based on CCP purchases optimization method, comprising the following steps:
S1, the mathematical model for establishing EH energy conversion devices and EH energy storage device, and show that EH's outputs and inputs relationship；
S2, EH renewable energy power output is calculated；
S3, by the mathematical model of EH energy conversion devices and EH energy storage device, EH output and input relationship and EH can
Renewable sources of energy power output establishes the EHS pneumoelectric based on CCP and purchases Optimized model, and purchases optimization mould using the EHS pneumoelectric based on CCP
Type is purchased EHS pneumoelectric and is optimized.
Step S1 specifically includes the following steps:
S11, the mathematical model for establishing EH energy conversion devices:
Calculate the power supply and heat supply power output of CHP unit, calculation formula are as follows:
In formula (1) and formula (2),WithRespectively the t moment power generation in the d days of CHP unit and heating function
Rate,For the amount of natural gas of unit time CHP unit consumption, ηCHP,eleAnd ηCHP,thThe respectively gentle heat of pneumoelectric transformation efficiency
Transformation efficiency；
Calculate the heats power of GF, calculation formula are as follows:
In formula (3),WithThe respectively amount of natural gas of GF heats power and unit time consumption, ηGFFor GF
Gas thermal transition efficiency；
Calculate the heat production power of GSHP, calculation formula are as follows:
In formula (4),WithThe electric energy of the heat production power of respectively GSHP and unit time consumption,
For the performance parameter of GSHP；
S12, the mathematical model for establishing EH energy storage device:
Charging and discharging state limitation are as follows:
The charge and discharge energy work rate limit are as follows:
Capacity limit are as follows:
In above formula, x is energy, indicates electricity, thermal energy with es and hs respectively,WithEnergy is filled for energy storage device
The variable (1 indicates working condition, and 0 indicates off working state) that state and exoergic state are 0-1,
WithRespectively energy-storage system charge and discharge energy work rate bound,WithEnergy exoergic power is respectively filled,For the d days t moment energy-storage system energy storage capacities, δxFor energy-storage system rate of energy loss, ηx,charge、ηx,dischargeRespectively
It fills, exergic efficiency,WithRespectively minimum energy storage capacity and maximum energy storage capacity,WithRespectively at the end of d days
It carves and the energy storage capacity of initial time, Δ t is time step 1h；
S13, outputting and inputting for single EH is obtained by the mathematical model of energy conversion devices each in EH and energy storage device
Relationship are as follows:
In formula (11),For electrical energy demands,For thermal demand, ηTransFor transformer efficiency,For from electricity
Online shopping electricity,The electricity generated for honourable resource through blower and photovoltaic cell group.
In embodiments of the present invention, step S2 specifically includes the following steps:
S21, photovoltaic generating system active power output P is calculatedPV, calculation formula are as follows:
In formula (13), r is radiancy, and A is the gross area of solar cell module, M, AmRespectively solar battery
The area of the battery pack number of packages of square matrix, m-th battery component, η, ηmRespectively total photoelectric conversion effect of solar cell module
Rate, m-th of photoelectric conversion efficiency, the probability density function of r are as follows:
In formula (14), Γ () is gamma function, rmaxFor greatest irradiation degree, α, β are Bata profile shape parameter；
S22, wind power output P is calculatedWT, calculation formula are as follows:
In formula (15), v is the wind speed of blower, vc、vfAnd vsWind speed, cut-out wind speed and rated wind speed are respectively cut,
PRFor the rated power of blower, a, b are constant, the probability density function of v are as follows:
In formula (16), K is form parameter, reflects the characteristic distributions of wind speed v, and C is scale parameter, and reflection this area is flat
The size of equal wind speed；
S23, renewable energy power output P is calculatedre, calculation formula are as follows:
Pre=PPV+PWT (17)。
In embodiments of the present invention, the specific steps of step S3 are as follows:
Optimization object is improved to as shown in figure 3, as an example by the EHS of the EH composition of 3 connections, each EH can be with table
Show a community in daily life, such as: cell, school etc..N in Fig. 3 is natural air-air source, by natural gas line and gas
Net is connected.Entire EHS deficiency pneumoelectric amount is supplied by natural gas line and power tie line by gas net power grid.
S31, objective function is established:
TC=min (CGas+CGrid) (18)
In formula (18), TCFor the totle drilling cost that pneumoelectric is purchased, CGasIt always purchases gas expense for EHS in dispatching cycle to use, CGridFor
Total power purchase expense, wherein
In above formula, n be EHS in EH number, T be one day 24 period, t=1,2, T, D be total activation week
Phase, i=1,2,3 be energy hub number in EHS,For power purchase price,It is the unit time when EHS electric power deficiency,
From power grid purchase of electricity (for positive value),For when EHS power surplus, to power grid electricity sales amount (for negative value),For sale of electricity valence
Lattice；
S32, setting bound for objective function, constraint condition includes formula (1)-formula (12), additionally include with
Lower constraint:
Electric equilibrium constraint:
Thermal balance constraint:
Amount of natural gas Constraints of Equilibrium:
0-1 constraint:
CHP unit electric heating units limits:
The hot units limits of GF unit:
The hot units limits of GSHP unit:
The pneumoelectric amount of purchasing chance constraint:
In above formula,WithElectricity is injected when electric power deficiency between respectively i-th of EH and adjacent EH or major network
Electricity is sent outside when with power surplus,Qigong rate is purchased for the EHS unit time,WithFor 0-1 decision variable,(for
Positive value) it is power purchase power maximum value,The absolute value of (for negative value) is sale of electricity power maximum value, Pi max(for positive value) is i-th
A EH injecting electric power maximum value, Pi minThe absolute value of (for negative value) is to send electrical power maximum value outside,AndThe respectively upper lower limit value of CHP unit electricity power output and heat power output,WithRespectively GF unit heat
Power output upper lower limit value,WithThe respectively maximin of GSHP unit heat power output, Pr { } are what event was set up
Probability,Qigong rate upper limit value, α are purchased for the EHS unit timegridAnd αgasRespectively correspond to power tie line, natural gas line
The confidence level of chance constrained programming condition；
S33, the pneumoelectric amount of purchasing chance constraint is converted to certainty constraint:
In above formula,Maximum transmission power is designed for power tie line,Maximum transmitted is designed for natural gas line
Power,WithRespectively 0-1 decision variable；
S34, Optimized model is purchased using the EHS pneumoelectric based on CCP EHS pneumoelectric is purchased and optimize.
It is optimization Simulation object that the present invention, which chooses the system that Fig. 3 includes 3 energy hubs, probe into natural gas line and
Power tie line transimission power limits the influence to the EHS pneumoelectric amount of purchasing optimization performance using chance constrained programming.It is dispatched in text
Period is D=30d, T=24h, time step 1h, and it is completely the same to set each device parameter in different energy sources hub.
Influence for verifying model to optimum results, point 2 kinds of different situations are probed into using chance constrained programming to not respectively
Chance constraint is carried out with object and limits the influence purchased to pneumoelectric, while meter and electrical load requirement respond the shadow purchased to pneumoelectric
It rings.√ is indicated under this kind of situation comprising the equipment or the factor in table, >=indicate the chance constraint that object is limited under this kind of situation
Planning confidence level is taken as 80%, 85%, 90%, 95%, 100% respectively.When confidence level is taken as 100%, namely to this
Limitation object is constrained using traditional conservative, rather than chance constraint.Detail classification is shown in Table 1.
Table 1
The purchase gas power purchase power of each period is generated into CDF curve, is convenient for comparative analysis.And consider the safety of system
Property, natural gas line design maximum transmission power is taken asConventional constraint purchases qigong rate upper limit valuePower tie line design maximum transmission power is taken asConventional constraint power purchase function
Rate upper limit value
1st kind of situation chance constraint only limits the purchase qigong rate of natural gas line between EHS and gas net, other constraints use
Conventional constraint.The confidence level upper limit takes 0.8, and applies different confidence levels 80%, 85%, 90%, 95% and 100%
Study lower influence of the chance constraint to optimum results of situation 1.The EHS pneumoelectric amount of purchasing CDF curve is as shown in Figure 4 under situation 1.And
And from when can be seen that purchase qigong rate is about changed between 14800kW from 0kW in Fig. 4 (a), cumulative probability becomes from 0
Change to 1.With the rising of confidence level, qigong rate is purchased in 0kW to 7400kW (including 7400kW) interval range, i.e. conventional constraint
The distribution of (when confidence level α=100%, all purchase gas performance numbers are distributed in the interval range) is more in interval range；
And in 7400kW to 14800kW interval range, i.e. distribution outside conventional constraint interval range is less.This is because renewable energy
The uncertainty of power output causes when renewable energy undercapacity in EHS, and it is negative to meet electric heating that EHS can purchase gas by power purchase
The demand of lotus, and lower confidence level can make the elasticity of natural gas system more preferable at this time, so in order to guarantee whole system
Economical operation, under the premise of meeting natural gas line transimission power chance constraint, EHS meets system by increasing purchase qigong rate
The demand of electrothermal load.When power purchase power is about changed between 4500kW from 0kW in Fig. 4 (b), power purchase power accumulation is general
Rate changes to 100% from 0, and under different confidence levels the shape and variation tendency of CDF curve have the characteristics that it is substantially similar.
2nd kind of situation chance constraint only limits the flow of power of power tie line between EHS and power grid, emulates remaining condition
It is completely the same with situation 1.The EHS pneumoelectric amount of purchasing CDF curve is as shown in Figure 5 under situation 2.
From Fig. 5 (a) it is found that cumulative probability becomes from 0 when purchase qigong rate is about changed between 7400kW from 0kW
Change to 1, and under different confidence levels the shape and variation tendency of CDF curve have the characteristics that it is substantially similar.With confidence level
Rising, purchasing distribution of the gas performance number in 7400kW increased significantly, and this is mainly due to apply chance constraint to power tie line
Afterwards, lower confidence level can make the elasticity of electric system stronger；When renewable energy undercapacity, electrothermal load can not be by
When meeting, due to purchasing the limitation of the gas upper limit of the power, system can deploy more power resources to meet electrothermal load demand, simultaneously
Also the number that purchase gas performance number reaches its upper limit can be reduced, and it can be seen that power purchase power about becomes from 0kW from Fig. 5 (b)
When change between 9000kW, cumulative probability changes to 1 from 0.With the rising of confidence level, power purchase power is arrived in 0kW
In 4500kW (not including 4500kW) interval range, power purchase power Distribution value is roughly the same, its power purchase performance number at 4500kW
It is distributed more；And in 4500kW to 9000kW interval range, power purchase power Distribution value is less.This is because renewable energy goes out
The uncertainty of power causes when renewable energy undercapacity in EHS, and EHS can purchase gas by power purchase to meet electrothermal load
Demand, and lower confidence level can make the elasticity of electric system more preferable at this time, so the economy in order to guarantee whole system
Operation, under the premise of meeting power tie line interaction power chance constraint, EHS meets system electric heating by increasing power purchase power
The demand of load.
It is as shown in table 2 that pneumoelectric acquisition expenses in dispatching cycle is acquired under every kind of situation difference confidence level.
Table 2
From 2 across comparison of table it can be seen that a little higher than situation 2 of the pneumoelectric acquisition expenses of situation 1, this is primarily due in EHS
Have GSHP, it can be achieved that electric energy to thermal energy unidirectional conversion, when natural gas system and electric system are all respectively provided with elasticity, consumption
The economy specific consumption natural gas that electric energy meets electrothermal load demand is high；From the point of view of longitudinal comparison, 2 kinds of situations are showed with machine
The rising of confidence level can be constrained, the characteristics of pneumoelectric acquisition expenses rises with it, this is primarily due to set with chance constraint
Believe horizontal rising, the elasticity of system is compressed, it is caused to purchase in certain period pneumoelectrics due to being limited to conventional constraint, can not
It realizes that most economical pneumoelectric purchases mode, thus its operating cost is caused to rise.
The foregoing is merely presently preferred embodiments of the present invention, is not intended to limit the invention, it is all in spirit of the invention and
Within principle, any modification, equivalent replacement, improvement and so on be should all be included in the protection scope of the present invention.
Claims (4)
1. a kind of EHS pneumoelectric based on CCP purchases optimization method, which comprises the following steps:
S1, the mathematical model for establishing EH energy conversion devices and EH energy storage device, and show that EH's outputs and inputs relationship；
S2, EH renewable energy power output is calculated；
S3, by the mathematical model of EH energy conversion devices and EH energy storage device, EH output and input relationship and EH is renewable
Energy power output establishes the EHS pneumoelectric based on CCP and purchases Optimized model, and purchases Optimized model pair using the EHS pneumoelectric based on CCP
EHS pneumoelectric, which is purchased, to be optimized.
2. the EHS pneumoelectric according to claim 1 based on CCP purchases optimization method, which is characterized in that the step S1 tool
Body the following steps are included:
S11, the mathematical model for establishing EH energy conversion devices:
Calculate the power supply and heat supply power output of CHP unit, calculation formula are as follows:
In formula (1) and formula (2),WithRespectively the t moment power generation in the d days of CHP unit and heats power,For the amount of natural gas of unit time CHP unit consumption, ηCHP,eleAnd ηCHP,thThe respectively gentle heat of pneumoelectric transformation efficiency turns
Change efficiency；
Calculate the heats power of GF, calculation formula are as follows:
In formula (3),WithThe respectively amount of natural gas of GF heats power and unit time consumption, ηGFFor GF gas heat
Transformation efficiency；
Calculate the heat production power of GSHP, calculation formula are as follows:
In formula (4),WithThe electric energy of the heat production power of respectively GSHP and unit time consumption,For
The performance parameter of GSHP；
S12, the mathematical model for establishing EH energy storage device:
Charging and discharging state limitation are as follows:
The charge and discharge energy work rate limit are as follows:
Capacity limit are as follows:
In above formula, x is energy, indicates electricity, thermal energy with es and hs respectively,WithEnergy is filled for energy storage device
State and exoergic state are the variable of 0-1,WithRespectively energy-storage system charge and discharge
Energy work rate bound,WithEnergy exoergic power is respectively filled,For t moment energy-storage system energy storage in the d days
Amount, δxFor energy-storage system rate of energy loss, ηx,charge、ηx,dischargeIt respectively fills, exergic efficiency,WithRespectively most
Small energy storage capacity and maximum energy storage capacity,WithThe respectively energy storage capacity of d days finish times and initial time, when Δ t is
Between step-length 1h；
S13, show that single EH's outputs and inputs relationship by the mathematical model of energy conversion devices each in EH and energy storage device
Are as follows:
In formula (11),For electrical energy demands,For thermal demand, ηTransFor transformer efficiency,To be purchased from power grid
Electricity,The electricity generated for honourable resource through blower and photovoltaic cell group.
3. the EHS pneumoelectric according to claim 1 based on CCP purchases optimization method, which is characterized in that the step S2 tool
Body the following steps are included:
S21, photovoltaic generating system active power output P is calculatedPV, calculation formula are as follows:
In formula (13), r is radiancy, and A is the gross area of solar cell module, M, AmRespectively solar cell array
Battery pack number of packages, m-th battery component area, η, ηmRespectively total photoelectric conversion efficiency of solar cell module, m
A photoelectric conversion efficiency, the probability density function of r are as follows:
In formula (14), Γ () is gamma function, rmaxFor greatest irradiation degree, α, β are Bata profile shape parameter；
S22, wind power output P is calculatedWT, calculation formula are as follows:
In formula (15), v is the wind speed of blower, vc、vfAnd vsRespectively cut wind speed, cut-out wind speed and rated wind speed, PRFor
The rated power of blower, a, b are constant, the probability density function of v are as follows:
In formula (16), K is form parameter, reflects the characteristic distributions of wind speed v, and C is scale parameter, reflects this area's average wind
The size of speed；
S23, renewable energy power output P is calculatedre, calculation formula are as follows:
Pre=PPV+PWT (17)。
4. the EHS pneumoelectric according to claim 1 based on CCP purchases optimization method, which is characterized in that the step S3's
Specific steps are as follows:
S31, objective function is established:
TC=min (CGas+CGrid) (18)
In formula (18), TCFor the totle drilling cost that pneumoelectric is purchased, CGasIt always purchases gas expense for EHS in dispatching cycle to use, CGridAlways to purchase
The electricity charge are used, wherein
In above formula, n be EHS in EH number, T be one day 24 period, t=1,2, T, D be the total activation period, i=
1,2,3 is energy hub number in EHS,For power purchase price,It is the unit time when EHS electric power deficiency, from power grid
Purchase of electricity,For when EHS power surplus, to power grid electricity sales amount,For sale of electricity price；
S32, setting bound for objective function, constraint condition includes formula (1)-formula (12), additionally include it is following about
Beam:
Electric equilibrium constraint:
Thermal balance constraint:
Amount of natural gas Constraints of Equilibrium:
0-1 constraint:
CHP unit electric heating units limits:
The hot units limits of GF unit:
The hot units limits of GSHP unit:
The pneumoelectric amount of purchasing chance constraint:
In above formula,WithElectricity and electric power are injected between respectively i-th of EH and adjacent EH or major network when electric power deficiency
Send electricity when superfluous outside,Qigong rate is purchased for the EHS unit time,WithFor 0-1 decision variable,For power purchase function
Rate maximum value,Absolute value be sale of electricity power maximum value, Pi maxFor i-th of EH injecting electric power maximum value, Pi minIt is absolute
Value is to send electrical power maximum value outside,AndRespectively CHP unit electricity power output and
The upper lower limit value of heat power output,WithRespectively GF unit heat power output upper lower limit value,WithRespectively
The maximin of GSHP unit heat power output, Pr { } are the probability that event is set up,The gas upper limit of the power is purchased for the EHS unit time
Value, αgridAnd αgasRespectively correspond to the confidence level of power tie line, natural gas line chance constrained programming condition；
S33, the pneumoelectric amount of purchasing chance constraint is converted to certainty constraint:
In above formula,Maximum transmission power is designed for power tie line,Maximum transmitted function is designed for natural gas line
Rate,WithRespectively 0-1 decision variable；
S34, Optimized model is purchased using the EHS pneumoelectric based on CCP EHS pneumoelectric is purchased and optimize.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810747877.5A CN108984894A (en) | 2018-07-09 | 2018-07-09 | A kind of EHS pneumoelectric based on CCP purchases optimization method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810747877.5A CN108984894A (en) | 2018-07-09 | 2018-07-09 | A kind of EHS pneumoelectric based on CCP purchases optimization method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108984894A true CN108984894A (en) | 2018-12-11 |
Family
ID=64536618
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810747877.5A Pending CN108984894A (en) | 2018-07-09 | 2018-07-09 | A kind of EHS pneumoelectric based on CCP purchases optimization method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108984894A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111769553A (en) * | 2020-07-06 | 2020-10-13 | 四川大学 | Optimal operation method of gas-electricity interconnection system considering wind electricity uncertainty |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090013697A1 (en) * | 2003-09-04 | 2009-01-15 | David Charles Landry | Simultaneous Underground Cavern Development and Fluid Storage |
| CN104392286A (en) * | 2014-12-02 | 2015-03-04 | 山东大学 | Microgrid operation optimizing method by considering combined supply of cooling, heating and power with stored energy operation strategy |
| CN107565696A (en) * | 2017-09-04 | 2018-01-09 | 四川大学 | A kind of energy accumulation capacity configuration based on Legnedre polynomial energy state function |
-
2018
- 2018-07-09 CN CN201810747877.5A patent/CN108984894A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090013697A1 (en) * | 2003-09-04 | 2009-01-15 | David Charles Landry | Simultaneous Underground Cavern Development and Fluid Storage |
| CN104392286A (en) * | 2014-12-02 | 2015-03-04 | 山东大学 | Microgrid operation optimizing method by considering combined supply of cooling, heating and power with stored energy operation strategy |
| CN107565696A (en) * | 2017-09-04 | 2018-01-09 | 四川大学 | A kind of energy accumulation capacity configuration based on Legnedre polynomial energy state function |
Non-Patent Citations (1)
| Title |
|---|
| 倪伟 等: "基于机会约束规划的能源集线器系统气电购置优化建模", 《电网技术》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111769553A (en) * | 2020-07-06 | 2020-10-13 | 四川大学 | Optimal operation method of gas-electricity interconnection system considering wind electricity uncertainty |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Jurasz et al. | Integrating a wind-and solar-powered hybrid to the power system by coupling it with a hydroelectric power station with pumping installation | |
| CN111445090B (en) | Double-layer planning method for off-grid type comprehensive energy system | |
| CN103490410B (en) | Micro-grid planning and capacity allocation method based on multi-objective optimization | |
| CN106849190B (en) | A kind of microgrid real-time scheduling method of providing multiple forms of energy to complement each other based on Rollout algorithm | |
| CN108229025B (en) | Economic optimization scheduling method for cooling, heating and power combined supply type multi-microgrid active power distribution system | |
| CN110689189B (en) | Combined cooling, heating and power supply and demand balance optimization scheduling method considering energy supply side and demand side | |
| CN105305419B (en) | A kind of independent micro-grid capacity configuration optimizing method containing compressed-air energy storage | |
| CN106786753B (en) | The system and its adjusting method of the Regional Energy internet of multi-user | |
| CN106026152A (en) | Charging and discharging scheduling method for electric vehicles connected to micro-grid | |
| CN106372742A (en) | Power-to-gas multi-source energy storage type microgrid day-ahead optimal economic dispatching method | |
| CN102623989A (en) | Method for optimization and configuration of intermittent distributed generation (DG) | |
| CN104537443A (en) | Co-generation type micro-grid economy coordination and optimization dispatching method | |
| CN111898806B (en) | Electric-thermal coupling source storage integrated multipotency flow park operation optimization method and system | |
| CN111934360B (en) | Virtual power plant-energy storage system energy collaborative optimization regulation and control method based on model predictive control | |
| CN110365013A (en) | A kind of capacity optimization method of photo-thermal-photovoltaic-wind-powered electricity generation combined generating system | |
| CN114330827B (en) | Distributed robust self-scheduling optimization method for multi-energy flow virtual power plant and application thereof | |
| CN102593855B (en) | Method for stabilizing fluctuation of output power of renewable energy power supply in power system | |
| CN106786610A (en) | A kind of photovoltaic high permeability line voltage optimization method based on energy-storage battery | |
| CN105186583B (en) | Energy router and its energy dispatching method based on multi-agent modeling | |
| CN114595868A (en) | Source network and storage collaborative planning method and system for comprehensive energy system | |
| CN103971183B (en) | The Optimizing Site Selection and capacity collocation method of a kind of photovoltaic plant | |
| CN115330021A (en) | Comprehensive energy operation optimization system and method considering methane electric heat utilization ratio | |
| CN102684228A (en) | Method for optimizing configuration of intermittent distribution type power supply based on complementary | |
| CN108491977A (en) | A kind of weak robust Optimization Scheduling of micro- energy net | |
| CN108984894A (en) | A kind of EHS pneumoelectric based on CCP purchases optimization method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181211 |
|
| RJ01 | Rejection of invention patent application after publication |