CN105046369B - A kind of modeling of electrical combined hybrid system and Optimization Scheduling based on energy centre - Google Patents
A kind of modeling of electrical combined hybrid system and Optimization Scheduling based on energy centre Download PDFInfo
- Publication number
- CN105046369B CN105046369B CN201510498181.XA CN201510498181A CN105046369B CN 105046369 B CN105046369 B CN 105046369B CN 201510498181 A CN201510498181 A CN 201510498181A CN 105046369 B CN105046369 B CN 105046369B
- Authority
- CN
- China
- Prior art keywords
- node
- energy
- gas
- variable
- compressor
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Abstract
The modeling of electrical combined hybrid system and Optimization Scheduling that the invention discloses a kind of based on energy centre, the present invention has initially set up electric power networks, natural gas network and energy centre model, and electric power networks and natural gas network couple to form electrical combined hybrid system by energy centre.Then using total energy cost as objective function, consider various constraint conditions and establish electrical combined hybrid system Optimized Operation mathematical model.It is solved with primal dual interior point method, slack variable and barrier parameter are successively introduced in solution procedure becomes model on the model of only equality constraint, it is re-introduced into Lagrange multiplier and obtains Lagrangian, the Nonlinear System of Equations formed later with its KKT condition of Newton Algorithm.The Simulation Example of construction is the result shows that the present invention is better than single optimization effect to the effect of optimization of electrical combined hybrid system.
Description
Technical field
The modeling of electrical combined hybrid system and Optimization Scheduling that the present invention relates to a kind of based on energy centre, belong to multiple-energy-source
Comprehensively control utilizes field.
Technical background
Being widely used of fossil energy leads to that environmental problem becomes increasingly conspicuous and fossil energy is petered out, traditional economy
It is difficult to continue with social development mode, needs to establish more efficient, environmental protection and Sustainable type using energy source mode, energy centre
The it is proposed of (energy hub, EH) makes it possible to establish new using energy source mode.Energy centre is defined as being turned by the energy
Change equipment and energy storage device constitutes, can be realized the pseudo-entity that various energy resources are mutually converted and stored.Energy centre is as more
The intermediary of kind energy form, promotes the connection of various energy resources form and energy resource system even closer, and consider economy
It influences.Energy centre input terminal connects the supply of various energy resources, and output end connects the consumption of various energy resources, and different energy sources form exists
It mutually converts therebetween, this Optimized Operation just for various energy resources form provides possibility.
Compared to other non-renewable energy, influence of the natural gas to environment be smaller, rich reserves and is easy to store;Gas turbine
With the increasingly development of combined-cycle power plant, shortens and regulate and control more so that less, the more efficient, construction period is invested in natural gas power
Add flexible and convenient;In addition, deepening constantly for " shale gas revolution ", will cause Gas Prices to be remarkably decreased.It is contemplated that naturally
The ratio of gas power generation will step up, natural gas also can by its many advantages in future source of energy using being occupied importantly in mode
Position.Natural gas network and electric power networks contact more and more closely, it may be said that following using energy source mode will be natural gas grid
The product that network and electric power networks are highly coupled.
Forefathers have enough research to the separate network of electric power networks and natural gas network, such as optimal load flow (optimal
Power flow, OPF) etc. classical problems, but it is rare consider electric power networks and natural gas network synthesis optimizing and scheduling problem.Closely
Nian Lai, since status of the natural gas in using energy source is more and more prominent, the comprehensive study of electric power networks and natural gas network is drawn
The extensive concern of domestic and foreign scholars is played.The research of electrical combined hybrid system is broadly divided into two ranks of coordinated planning and coordinated operation
Section, problem are concentrated mainly on mutually converting and influencing, the association of supervision department for the processing of uncertain factor, electric load and electric load
Adjust planning and management and the dynamic behaviour of combined hybrid system and reliability etc..
Electrical combined hybrid system Optimized Operation is substantially nonlinear programming problem.The present invention is using total energy cost as mesh
Scalar functions, it is contemplated that active and reactive power Constraints of Equilibrium in electric power networks, balance nodes phase angle Constraints of Equilibrium, generated power,
Idle units limits, voltage constraint;Flow equilibrium constrains in natural gas network, compressor quantity of state Constraints of Equilibrium, and gas source point is defeated
Constraint and node pressure constraint out;The input and output Constraints of Equilibrium and scheduling factor of energy centre constrain.The present invention is using former-right
Even interior point method (primal-dual interior point method, PDIPM) solves the electrical mixed connection system based on energy centre
System Optimal Scheduling, because PDIPM has, convergence is good, and computational efficiency is high, strong robustness, and initial value selection is insensitive, does not have
Identify the advantages that constraint set is difficult.Emulation finally is programmed to the example of construction, interpretation of result demonstrates the method for the present invention
Advantage.
Summary of the invention
Goal of the invention: the present invention is a kind of based on the new of energy centre for providing the technical issues of solution needed for the prior art
The modeling of type energy supply form, that is, electrical combined hybrid system and Optimized Operation.
Technical solution: the present invention to achieve the above object, adopts the following technical scheme that
The present invention is a kind of based on the electrical combined hybrid system modeling of energy centre concept and Optimization Scheduling, it is characterised in that
The method is successively realized according to the following steps:
1) parameter information of electric power networks is obtained, comprising: the headend node and endpoint node of transmission line of electricity are numbered, branch π
The resistance of type equivalent circuit, reactance, over the ground shunt conductance, susceptance, transformer voltage ratio and impedance, each node load and generator
Export active and reactive constraint, each node voltage constraint;
For each electric power networks interior joint i:
Wherein: Pi, QiRespectively electric power networks interior joint i active power and reactive power;ei, fiRespectively node i voltage
The real and imaginary parts of vector;Gij, BijThe respectively real and imaginary parts of the i-th row of node admittance matrix jth column element.
2) parameter information of natural gas network is obtained, comprising: the headend node and endpoint node of gas pipeline are numbered, pipeline
Physical characteristics, the quantity of state of compressor, each node gas load and the gas source point natural gas such as length, internal diameter and efficiency of transmission it is defeated
It constrains out, each node pressure constraint;
Introduce includes pipeline flow equation, compressor flowrate consumption equation and flux balance equations;
Under ideal conditions, pipeline k can be indicated from node i to the flow value of node j with following equation:
Wherein: fkijFor pipeline flow value;FkFor coefficient of pipe friction;DkFor internal diameter of the pipeline;G
For gas-gravity coefficient;LkFor duct length;πiFor node i pressure value;πjFor node j pressure value;π0For reference pressure value;T0For
Normal temperature value;TkaFor mean gas temperature;ZaFor the average gas compressed coefficient;
For the complete turbulent form of high voltage network, flow equation can be further simplified are as follows:
Wherein:ε is pipeline efficiency;
Under ideal gas conditions, the energy consumption equation of compressor can be indicated are as follows:
Wherein:fCkTo pass through the gas flow of compressor;πiCompressor is injected for gas
Pressure;πjCompressor pressure is exported for gas;ZkiFor the gas compressibility factor of suction port of compressor;TkiFor compressor output temperature;α
For heat gain coefficient;ηkFor compressor efficiency;
It is converted into the flow value of consumption:
τk=αTk+βTkHkij+γTkHk 2 ij;
The flux balance equations of each node can be indicated with following matrix form:
(A+U) τ=0 f+w-T;
Wherein: A, U and T is respectively that branch, compressor and compressor consume direction matrix related with node, and f is branch
Flow value vector;W is that the gas of each node injects vector;τ is each compressor consumed flow value vector;
3) parameter information of energy centre is obtained, comprising: each energy centre input and output and electric power networks and natural gas grid
The connection of network, each energy centre load and the constraint of each energy centre scheduling factor;
Multiple-input and multiple-output energy centre model can be described with matrix equation:
Wherein: P and L is respectively input and output vector;Matrix C is referred to as coupling matrix;
List energy centre input and output equilibrium equation:
Wherein: ηTIndication transformer efficiency;ηGTeIndicate that gas turbine gas turns the efficiency of electricity;ηGThIndicate that gas turbine gas turns
The efficiency of heat;ηFIndicate the efficiency of gas fired-boiler;
4) using total energy cost as objective function f (x), control variable u and state variable are chosenAccording to it is various about
Beam condition column write equality constraint and inequality constraints, establish electrical combined hybrid system nonlinear programming mathematics model and are asked with PDIPM
Solution;
Establish electrical combined hybrid system Optimal Operation Model:
min f(x);
S.t. h (x)=0;
gmin≤g(x)≤gmax;
Wherein: h indicates equality constraint;G indicates inequality constraints;gmax, gminRespectively indicate the upper and lower of inequality constraints
Limit;
5) broad sense slack variable l and u is introduced, converts equality constraint for inequality constraints in model;
Introduce broad sense slack variable l=[l1,…,lr]T, u=[u1,…,ur]TEquation is converted by broad sense inequality:
g(x)+u-gmax=0;
g(x)-l-gmin=0;
It should meet: u > 0, l > 0;
6) in order to guarantee that objective function obtains minimum value in feasible zone, log-barrier parameter μ is introduced to objective function, it will
Model conversation is the model for containing only equality constraint;
Slack variable and barrier parameter are introduced by model conversation into the model of only equality constraint:
S.t. h (x)=0;
g(x)+u-gmax=0;
g(x)-l-gmin=0;
Wherein: Discontinuous Factors (barrier parameter) μ > 0;
It introduces Lagrange multiplier and obtains Lagrangian:
7) above-mentioned equality constraint model is solved with method of Lagrange multipliers, introduces Lagrange multiplier z, w and y, is drawn
Ge Lang function, according to Lagrangian extreme value, there are necessary condition (KKT condition) column to write Nonlinear System of Equations;
One group of Nonlinear System of Equations is obtained according to Lagrangian extreme value existence condition (KKT condition), with newton-La Fuxunfa
It solves, single order update equation are as follows:
Wherein: Lx, Ly, Lz, Lw, changed for KKT equation last time
The residual error in generation;The Hessian matrix of respectively h (x) and g (x);
8) duality gap Gap=l is definedTz-uTW and Discontinuous FactorsWherein σ=[0,1] is center parameter;
9) the above-mentioned Nonlinear System of Equations of Newton Algorithm is used, Jacobian matrix, Hessian matrix and constant term etc. are calculated, is solved
Update equation group obtains the correction amount of each former variable and dual variable, is multiplied by step-length and is modified to variable, until duality gap
Gap is less than convergence precision (ε=10-6), otherwise do not restrain;
It solves update equation and obtains correction amount, calculate former variable and dual variable step-length:
Variable is corrected as the following formula:
10) Optimized Operation result is exported, comprising: electric power networks generated power, idle power output and each node pressure, naturally
The output of gas network gas source point natural gas and each node pressure value, each scheduling factor of energy centre.
The utility model has the advantages that the present invention is in terms of existing technologies: the present invention has initially set up electric power networks, natural gas network
With energy centre model, electric power networks and natural gas network couple to form electrical combined hybrid system by energy centre.Then with total
Energy cost be objective function, consider various constraint conditions and establish electrical combined hybrid system Optimized Operation mathematical model.With original-
Dual interior point solves, and successively introducing slack variable and barrier parameter in solution procedure becomes model the mould of only equality constraint
Type is re-introduced into Lagrange multiplier and obtains Lagrangian, the non-linear side formed later with its KKT condition of Newton Algorithm
Journey group.The Simulation Example of construction is the result shows that the present invention is better than single optimization effect to the effect of optimization of electrical combined hybrid system.
Detailed description of the invention
Fig. 1 is the method for the present invention flow chart;
Fig. 2 is energy centre model schematic;
Fig. 3 is energy centre example schematic diagram;
Fig. 4 is the different energy sources center input results comparison diagram of numerical results analysis.
Specific embodiment
The techniqueflow of invention is described in detail with reference to the accompanying drawing:
Power network model embodiment
The topological structure of electric power networks can be described with admittance matrix Y.The state variable of electric power networks has node voltage, section
Point power and branch power etc., here with the Cartesian form of voltage to egress active power and reactive power:
Wherein: Pi, QiRespectively electric power networks interior joint i active power and reactive power;ei, fiRespectively node i voltage
The real and imaginary parts of vector;Gij, BijThe respectively real and imaginary parts of the i-th row of node admittance matrix jth column element.
Natural gas network model embodiment
Typical natural gas network includes one or more gas source points (natural gas production and storage point), one or more
Load (power plant, other networks or storage point), pipeline, compressor and other equipment.Base for natural gas network modelling
There are three essentialities: 1) pipeline;2) compressor;3) decanting point (node or bus).Two most important parts of natural gas model
It is pipeline flow model and compressor model.
Natural gas line flow equation describes natural gas flow magnitude and pipe ends pressure value and relevant gaseousness
Relationship between matter, pipe characteristic and operating status.It forms flow equation to need to consider several factors, none flow equation
Suitable for all situations.Under ideal conditions, pipeline k can be indicated from node i to the flow value of node j with following equation:
Wherein: fkijFor pipeline flow value;FkFor coefficient of pipe friction;DkFor internal diameter of the pipeline;G
For gas-gravity coefficient;LkFor duct length;πiFor node i pressure value;πjFor node j pressure value;π0For reference pressure value;T0For
Normal temperature value;TkaFor mean gas temperature;ZaFor the average gas compressed coefficient.
For the complete turbulent form of high voltage network, flow equation can be further simplified are as follows:
Wherein:ε is pipeline efficiency.
Compressor is mounted in gas network to transmit gas and the compensation downstream pressure as caused by frictional resistance etc.
Loss.The operation of compressor needs to consume a large amount of energy, in large-scale natural gas network, drives main compressor most economical
Energy source is the natural gas by compressor.Equation for the network analysis key of compressor is energy consumption, it is to flow through
The gas flow of compressor and injection, flow out compressor gas pressure function.Under ideal gas conditions, in addition experience
Quantization, the energy consumption equation of adiabatic compressor can indicate are as follows:
Wherein:fCkTo pass through the gas flow of compressor;πiCompressor is injected for gas
Pressure;πjCompressor pressure is exported for gas;ZkiFor the gas compressibility factor of suction port of compressor;TkiFor compressor output temperature;α
For heat gain coefficient;ηkFor compressor efficiency.
It is converted into the flow value of consumption:
According to Kirchhoff's first law: all inflows are zero with the flow algebraical sum for flowing out certain network node.Therefore it is each
The flux balance equations of a node can be indicated with following matrix form:
(A+U) τ=0 f+w-T;
Wherein: A, U and T is respectively that branch, compressor and compressor consume direction matrix related with node, and f is branch
Flow value vector;W is that the gas of each node injects vector;τ is each compressor consumed flow value vector.
Energy centre model embodiment
Energy centre, which is defined as being made of energy conversion equipment and energy storage device, can be realized various energy resources mutually converts
With the pseudo-entity of storage.It is connected with each other between energy centre by energy source equipment, constitutes multiple-energy-source network system.In energy
Source central interior, the energy may be consumed by load, or be converted into other forms.
Energy centre is there are many energy form α, β ... ∈ ε, each energy form can be energy centre input or
Output.These input energies (output energy) are defined as Pα,Pβ..., Pω(Lα,Lβ,…Lω).As shown in Figure 2.For how defeated
The energy centre for entering multi output can be described with following matrix equation:
Wherein: P and L is respectively input and output vector;Matrix C is referred to as coupling matrix.Coupling matrix is mathematically described
Energy is input to from energy centre the distribution of output.Element in coupling matrix is coupling factor.Each coupling factor table
Show that is specifically input to a specific output.
Fig. 3 gives a kind of example of energy centre.Energy centre input is provided by electric power and natural gas, and the energy is in the energy
Central interior supplies electric load and thermic load by transformer, gas turbine and gas fired-boiler conversion, output.Wherein ν is scheduling
The factor.Energy centre input and output equilibrium equation can be listed accordingly:
Le=ηTPe+νηGTePg;
Lh=ν ηGThPg+(1-ν)ηFPg;
Write as matrix form:
Wherein: ηTIndication transformer efficiency;ηGTeIndicate that gas turbine gas turns the efficiency of electricity;ηGThIndicate that gas turbine gas turns
The efficiency of heat;ηFIndicate the efficiency of gas fired-boiler.
Electrical combined hybrid system Optimal Operation Model embodiment
According to above-mentioned model, electric power networks are coupled with natural gas network by energy centre, form electrical combined hybrid system.With
Total energy cost is objective function, considers the various constraints of electric power networks, natural gas network and energy centre, is established electrical mixed
Connection system Optimal Operation Model:
min f(x);
S.t. h (x)=0;
gmin≤g(x)≤gmax;
Wherein: h indicates equality constraint;G indicates inequality constraints;gmax, gminRespectively indicate the upper and lower of inequality constraints
Limit.
Selection variable is carried out, variable is controlled:
1) electric power networks control variable: generated power power output PG;Generator reactive power output QR;
2) natural gas network-control variable: gas source point deliverability of gas NG;
3) energy centre controls variable: scheduling factor ν.
The then control variable that the present invention chooses are as follows:
U=[PG1,...,PGng,QG1,...,QGng,NG1,…NGnGS,ν1,…,νnEH]T;
Wherein: ngFor electric power networks generator number of units;nGSFor natural gas network gas source number;nEHFor energy centre number.
State variable:
1) power networks state variable: the real part e and imaginary part f of node voltage vector;
2) natural gas network state variables: the pressure π of the node and flow f by compressorC;
3) energy centre state variable: the input quantity P of electric energyeWith the input quantity P of natural gasg。
The then state variable that the present invention chooses are as follows:
Wherein: nbFor the node number of electric power networks;nGBFor natural gas network node number;nGCTo be pressed in natural gas network
The number of contracting machine.
Total variable-definition are as follows:
The present invention is using the total energy cost of system as objective function:
Wherein: a2i, a1i, a0iFor energy loss parameter of curve;PiIt here include generator in power network for generalized energy
The active power and natural gas network gas source point of sending provide the energy of gas.
Equality constraint
1) electric power networks equality constraint:
ΔPi=PGi-PDi-Pei-Pi=0;
ΔQi=QRi-QDi-Qi=0;
Wherein: Δ Pi, Δ QiFor each node active and reactive power amount of unbalance;Δ angle is balance nodes angle restriction;
tanθbalFor the angle tangent value of balance nodes;PGi, QRiThe respectively active and reactive power output of generator i;PDi, QDiRespectively save
The active and reactive load of point i;PeiIt is inputted for electric energy of the node i to energy centre.
2) natural gas network equality constraint:
ΔWi=NGi+wi-Fi-Pgi=0;
Wherein: Δ WiFor the amount of unbalance of node-flow magnitude each in natural gas network;NGiIt is gas source point to the gas of node i
Injection rate;wi, FiThe respectively gas injection rate and load of node i;PgiIt is inputted for natural gas of the node i to energy centre;
ΔPCkFor the amount of unbalance of compressor k-state amount in natural gas network;RkRespectively the given value of compressor k-state amount and
Calculated value.
3) energy centre equality constraint:
ΔLei=Lei-ηTPei-νiηGTePgi;
ΔLhi=Lhi-νiηGThPgi-(1-νi)ηFPgi;
Wherein: Δ Lei, Δ LhiThe respectively amount of unbalance of the power budget and thermic load of energy centre i;Lei, LhiRespectively
For the power budget and thermic load of energy centre i;Pei, PgiRespectively the electric energy of energy centre i and natural gas input;νiFor the energy
The scheduling factor of center i.
Inequality constraints
1) electric power networks inequality constraints:
PGimin≤PGi≤PGimax;
QRimin≤QRi≤QRimax;
Wherein: PGimin, PGimaxThe lower and upper limit of active power are issued by generator;QRimin, QRimaxFor generator institute
Send out the lower and upper limit of reactive power;For the lower and upper limit of node voltage amplitude square.
2) natural gas network inequality constraints:
0≤NGi≤NGimax;
πimin≤πi≤πimax;
Wherein: NGimaxFor the upper limit of gas source point gas supply each in natural gas network;πiminWith πimaxRespectively each section
Press the lower and upper limit of force value.
3) energy centre inequality constraints:
0≤νi≤1;
It is solved by primal dual interior point method, is firstly introduced into broad sense slack variable l=[l1,…,lr]T, u=[u1,…,ur]T
Equation is converted by broad sense inequality:
g(x)+u-gmax=0;
g(x)-l-gmin=0;
It should meet: u > 0, l > 0.
Then, in order to guarantee that objective function obtains minimum value in feasible zone, log-barrier parameter is introduced to objective function
μ works as u, and when l is close to restrained boundary, objective function tends to be infinitely great.Optimal Operation Model is the mould of only equality constraint at this time
Type:
S.t. h (x)=0;
g(x)+u-gmax=0;
g(x)-l-gmin=0;
Wherein: Discontinuous Factors (barrier parameter) μ > 0.
At this moment it can directly be solved by method of Lagrange multipliers.Introduce Lagrange multiplier y=[y1,...,ym]T, z=
[z1,…,zr]T, w=[w1,...,wm]TObtain Lagrangian:
According to necessary condition existing for Lagrangian extreme value (KKT condition) available Nonlinear System of Equations:
Wherein: e is that each element is 1rDimensional vector;L=diag (l);U=diag (u);Z=diag (z);W=diag
It (w) is diagonal matrix.
It calculates available:
Define complementary gap (duality gap) are as follows: Gap=lTz-uTw.Practice have shown that the Discontinuous Factors μ in objective function is pressed
Following formula computational convergence is preferable:
Wherein: σ=[0,1] is center parameter.When σ=1, be conducive to the feasibility for improving algorithm, but between reducing antithesis
Gap is without effect;σ=0 reduces the duality gap of algorithm, is conducive to the optimality of solution.It is mentioned as far as possible on the basis of guaranteeing feasibility
The optimality of high solution, σ generally take 0.1.
Above-mentioned Nonlinear System of Equations, the matrix form of single order update equation group can be solved with newton-La Fuxunfa are as follows:
Wherein: Lx, Ly, Lz, Lw, changed for KKT equation last time
The residual error in generation;The Hessian matrix of respectively h (x) and g (x).
It solves update equation and obtains correction amount, calculate former variable and dual variable step-length:
Variable is corrected as the following formula:
Successively iteration, until duality gap Gap is less than convergence precision.
Embodiment
To verify effect of the present invention, IEEE14 node power network in Matpower and certain 14 node natural gas network are led to
Cross one example of energy centre coupling configurations in Fig. 3.Wherein electric power networks and natural gas network are joined by 5 energy centres
5 load points that natural gas network Central Plains is come are connected to energy centre by system, and the natural gas as energy centre inputs;It chooses
Heavier 5 load points of duty ratio are connect with energy centre in electric power networks, and the electric energy as energy centre inputs, energy centre
Load side considers power budget and thermal energy output, and considers the influence of energy centre scheduling factor simultaneously.Each energy centre and electricity
Power network and the specific connection of natural gas network and energy centre output situation are shown in Table 1.
It is programmed and is emulated with example of the Matlab to construction, simulation result is analyzed.
It is limited in 500~1500psia when by node pressure value each in natural gas network, takes natural gas damage curve system
Number a2=0, a1=2, a0=0, when energy centre output takes data in table 1, the minimum cost of combined hybrid system complex optimum is $
5338.5239.The minimum cost of power network single optimization is 8081.526, and the minimum cost of natural gas network single optimization is
a1·NGGHV/baseMVA=89.98, therefore when the totle drilling cost of power network and natural gas single optimization is 8081.526+
89.98=$ 8171.506.It can be seen that the complex optimum cost of combined hybrid system is significantly lower than single optimization cost, power network is illustrated
The complex optimum of network and natural gas grid combined hybrid system has significant economic benefit.
In order to more fully illustrate combined hybrid system complex optimum cost lower than single optimization cost, pressure same as above is chosen
Power range and natural gas damage curve coefficient emulate different energy centre power output value and thermal energy output valve, optimization knot
Fruit is as shown in Figure 4.From comparison diagram it is found that exporting for different energy centres, combined hybrid system complex optimum cost is below solely
The cost of vertical optimization, further illustrates the advantage of combined hybrid system complex optimum.
The connection of 1 combined hybrid system node of table is exported with energy centre
Claims (3)
1. a kind of modeling of electrical combined hybrid system and Optimization Scheduling based on energy centre, it is characterised in that: including following step
It is rapid:
1) parameter information of electric power networks is obtained, comprising: the headend node and endpoint node of transmission line of electricity are numbered, branch π type etc.
Resistance, the reactance of circuit are imitated, over the ground shunt conductance, susceptance, transformer voltage ratio and impedance, each node load and generator export
Active and reactive constraint, each node voltage constraint;
2) parameter information of natural gas network is obtained, comprising: the headend node and endpoint node of gas pipeline are numbered, the length of pipeline
Degree, internal diameter and efficiency of transmission, the quantity of state of compressor, each node gas load and gas source point natural gas output constraint, each node
Pressure confines;
3) multiple-input and multiple-output energy centre model can be described with matrix equation:
Wherein: P and L is respectively input and output vector, and factor indicates the energy of input energy sources in P, and factor indicates the output energy in L
Energy;Matrix C is referred to as coupling matrix, in C factor indicate output the energy and input energy sources coupled relation, referred to as coupling because
Son;
List energy centre input and output equilibrium equation:
Wherein: ηTIndication transformer efficiency;ηGTeIndicate that gas turbine gas turns the efficiency of electricity;ηGThIndicate that gas turbine gas turns heat
Efficiency;ηFIndicate the efficiency of gas fired-boiler;ν indicates scheduling factor;
4) electrical combined hybrid system Optimal Operation Model is established:
min f(x);
S.t.h (x)=0;
gmin≤g(x)≤gmax;
Wherein: f (x) indicates objective function;H indicates equality constraint;G indicates inequality constraints;gmax, gminRespectively indicate inequality
The upper and lower limit of constraint;
5) slack variable and barrier parameter are introduced by model conversation into the model of only equality constraint:
S.t.h (x)=0;
g(x)+u-gmax=0;
g(x)-l-gmin=0;
Wherein: lj、ujIndicate the element in slack variable l, u;K indicates variable label;The dimension of r expression slack variable;Obstacle ginseng
Number μ > 0, barrier parameter are also referred to as Discontinuous Factors;
6) it introduces Lagrange multiplier and obtains Lagrangian:
Wherein: yT、zT、wTIndicate the transposition of Lagrange multiplier y, z, w;L, u indicates slack variable;
7) one group of Nonlinear System of Equations is obtained according to Lagrangian extreme value existence condition, is solved with newton-La Fuxunfa, single order
Update equation are as follows:
Wherein: Lx, Ly, Lz, Lw,For KKT equation last iteration
Residual error;WithThe Hessian matrix of respectively h (x) and g (x);T indicates transposition;I indicates unit matrix;y,z,
W indicates Lagrange multiplier;X indicates variable;Δ z, Δ l, Δ w, Δ u, Δ x, Δ y indicate the variable quantity of relevant variable;
8) it solves update equation and obtains correction amount, calculate former variable and dual variable step-length:
Wherein: αpIndicate former variable step-length;αdIndicate dual variable step-length;K indicates variable label;
9) variable is corrected as the following formula:
10) press above-mentioned steps iteration, until duality gap be less than convergence precision, obtain the optimal of electrical combined hybrid system Optimized Operation
Solution.
2. the modeling of electrical combined hybrid system and Optimization Scheduling according to claim 1 based on energy centre, feature
It is: in the step 1), for electric power networks interior joint i:
Wherein: Pi, QiRespectively electric power networks interior joint i active power and reactive power;ei, fiRespectively node i voltage vector
Real and imaginary parts;Gij, BijConductance and susceptance respectively between node i and node j.
3. the modeling of electrical combined hybrid system and Optimization Scheduling according to claim 1 based on energy centre, feature
Be: in the step 2), introducing includes pipeline flow equation, compressor flowrate consumption equation and flux balance equations;
Under ideal conditions, pipeline K can be indicated from node m to the flow value of node n with following equation:
Wherein: fKmnFor pipeline flow value;FKFor coefficient of pipe friction;DKFor internal diameter of the pipeline;G is
Gas-gravity coefficient;LKFor duct length;πmFor node m pressure value;πnFor node n pressure value;π0For reference pressure value;T0For mark
Quasi- temperature value;TKaFor mean gas temperature;ZaFor the average gas compressed coefficient;
For the complete turbulent form of high voltage network, flow equation can be further simplified are as follows:
Wherein:ε is pipeline efficiency;
Under ideal gas conditions, the energy consumption equation of compressor can be indicated are as follows:
Wherein:fCKTo pass through the gas flow of compressor;πmNode m is injected for compressor gas
Pressure;πnFor the pressure of compressor gas output node n;ZKmFor the gas compressibility factor of suction port of compressor;TKnFor compressor
Output temperature;α is heat gain coefficient;ηKFor compressor efficiency;
It is converted into the flow value of consumption:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510498181.XA CN105046369B (en) | 2015-08-13 | 2015-08-13 | A kind of modeling of electrical combined hybrid system and Optimization Scheduling based on energy centre |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510498181.XA CN105046369B (en) | 2015-08-13 | 2015-08-13 | A kind of modeling of electrical combined hybrid system and Optimization Scheduling based on energy centre |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105046369A CN105046369A (en) | 2015-11-11 |
CN105046369B true CN105046369B (en) | 2019-04-30 |
Family
ID=54452896
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510498181.XA Active CN105046369B (en) | 2015-08-13 | 2015-08-13 | A kind of modeling of electrical combined hybrid system and Optimization Scheduling based on energy centre |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105046369B (en) |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105958531B (en) * | 2016-05-27 | 2018-10-19 | 清华大学 | A kind of electric-gas coupling multipotency flow network method for estimating state |
CN106096777B (en) * | 2016-06-12 | 2019-05-21 | 清华大学 | A kind of Optimization Scheduling of electric-gas coupling multipotency streaming system |
CN106056251B (en) * | 2016-06-12 | 2019-06-18 | 清华大学 | A kind of Optimization Scheduling of electric-thermal coupling multipotency streaming system |
CN106056478B (en) * | 2016-06-12 | 2019-05-21 | 清华大学 | The Interval Power Flow calculation method of heat supply network in a kind of electric-thermal coupled system |
CN106447152A (en) * | 2016-06-21 | 2017-02-22 | 中国南方电网有限责任公司电网技术研究中心 | Energy center-based power-gas-heat system co-scheduling method and system |
CN108510299B (en) * | 2017-02-28 | 2021-12-10 | 腾讯科技(北京)有限公司 | Method and device for determining multimedia resource delivery strategy |
CN107292456B (en) * | 2017-08-01 | 2020-11-24 | 重庆大学 | Electric-gas energy flow distribution type collaborative optimization calculation method based on alternating direction multiplier method |
CN107394828A (en) * | 2017-08-24 | 2017-11-24 | 河海大学 | A kind of electrical interconnection integrated energy system Optimal Operation Analysis method based on Probabilistic Load Flow |
CN107404118A (en) * | 2017-09-06 | 2017-11-28 | 河海大学 | Electrical interconnection system probability optimal load flow computational methods based on stochastic response surface |
CN108062599B (en) * | 2017-12-12 | 2020-08-18 | 西安交通大学 | Modeling method for cooperative operation of electric power system and natural gas system |
CN108258679B (en) * | 2017-12-25 | 2020-05-12 | 国网浙江省电力有限公司经济技术研究院 | Electric-thermal comprehensive energy system optimization scheduling method considering heat storage characteristics of heat supply network |
CN108846507B (en) * | 2018-05-29 | 2020-11-10 | 西安交通大学 | Day-ahead economic dispatching method of electric-gas coupling system based on mixed integer second-order cone programming |
CN109086474A (en) * | 2018-06-12 | 2018-12-25 | 河海大学 | Natural gas system optimal load flow method for solving based on second order cone optimization algorithm |
CN108596525B (en) * | 2018-06-29 | 2022-07-22 | 国家电网有限公司 | Robust optimization scheduling method for micro-energy network with complementary cold-heat-electricity multi-energy |
CN109190849B (en) * | 2018-10-18 | 2022-05-31 | 清华-伯克利深圳学院筹备办公室 | Energy price prediction method and device, pricing equipment and storage medium |
CN110428185B (en) * | 2019-08-08 | 2020-11-03 | 河海大学 | Pseudo measurement model-based robust state estimation method for electricity-heat interconnection comprehensive energy system |
CN112039750A (en) * | 2020-09-01 | 2020-12-04 | 四川大学 | Electrical combined system network state estimation method based on comprehensive energy system |
CN112069638A (en) * | 2020-09-08 | 2020-12-11 | 广东电网有限责任公司电力科学研究院 | Low-voltage distribution network energy storage equipment access point location method and related device |
CN112818492B (en) * | 2021-02-02 | 2022-12-02 | 山东大学 | Electric-gas coupling network energy flow solving method considering state variables of energy source station |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102522746A (en) * | 2011-12-13 | 2012-06-27 | 河海大学 | VSC-HVDC AC/DC optimal power flow method based on primal-dual interior point algorithm |
CN103366097A (en) * | 2013-07-24 | 2013-10-23 | 国家电网公司 | Calculation method of optimal power flow based on class extension variable interior point method |
CN103793755A (en) * | 2014-01-14 | 2014-05-14 | 河海大学 | Method for obtaining optimal power flow of mixed direct-current transmission system based on prime-dual interior point method |
CN104392285A (en) * | 2014-11-28 | 2015-03-04 | 河海大学 | Optimal power flow acquiring method for electrical power system including hybrid DC power transmission |
CN104393592A (en) * | 2014-11-27 | 2015-03-04 | 河海大学 | Temperature influence-considering optimal power flow algorithm of power system |
-
2015
- 2015-08-13 CN CN201510498181.XA patent/CN105046369B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102522746A (en) * | 2011-12-13 | 2012-06-27 | 河海大学 | VSC-HVDC AC/DC optimal power flow method based on primal-dual interior point algorithm |
CN103366097A (en) * | 2013-07-24 | 2013-10-23 | 国家电网公司 | Calculation method of optimal power flow based on class extension variable interior point method |
CN103793755A (en) * | 2014-01-14 | 2014-05-14 | 河海大学 | Method for obtaining optimal power flow of mixed direct-current transmission system based on prime-dual interior point method |
CN104393592A (en) * | 2014-11-27 | 2015-03-04 | 河海大学 | Temperature influence-considering optimal power flow algorithm of power system |
CN104392285A (en) * | 2014-11-28 | 2015-03-04 | 河海大学 | Optimal power flow acquiring method for electrical power system including hybrid DC power transmission |
Non-Patent Citations (1)
Title |
---|
大规模风电场并网的优化调度建模及算法研究;孙小璞;《中国优秀硕士学位论文全文数据库》;20131215(第S2期);全文 |
Also Published As
Publication number | Publication date |
---|---|
CN105046369A (en) | 2015-11-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105046369B (en) | A kind of modeling of electrical combined hybrid system and Optimization Scheduling based on energy centre | |
CN103490410B (en) | Micro-grid planning and capacity allocation method based on multi-objective optimization | |
CN104734168B (en) | Microgrid running optimization system and method based on power and heat combined dispatching | |
CN104734155B (en) | A kind of method obtaining electrical interconnection energy resource system available transmission capacity | |
CN110866627B (en) | Multi-zone electricity-gas coupling comprehensive energy system optimal scheduling method considering step gas price | |
CN105939029B (en) | Integrated energy system programme acquisition methods and system | |
CN105576657B (en) | Random fuzzy continuous power flow modeling method for active power distribution network and power transmission network integrated system | |
CN108846507A (en) | Electric-gas coupled system based on MIXED INTEGER Second-order cone programming economic load dispatching method a few days ago | |
CN103280821B (en) | Multi-period dynamic reactive power optimization method of intelligent power distribution system | |
CN105119275B (en) | A kind of algorithm of meter and the Electrical Power System Dynamic optimal load flow of THE UPFC | |
CN108023364A (en) | Power distribution network distributed generation resource maximum access capability computational methods based on convex difference planning | |
CN107196315A (en) | The extendable power-less optimized controlling method of the power distribution network containing light-preserved system | |
CN107910871A (en) | A kind of energy internet is provided multiple forms of energy to complement each other control method | |
CN107732982A (en) | Consider the integrated energy system Multiple Time Scales dispatching method of Model Predictive Control | |
CN112035984B (en) | Collaborative planning method for comprehensive energy system of electricity-gas-storage area | |
CN106655226A (en) | Active power distribution network asymmetric operation optimization method based on intelligent soft open point | |
Fan et al. | Review on coordinated planning of source-network-load-storage for integrated energy systems | |
CN102163845B (en) | Optimal configuration method of distributed generations (DG) based on power moment algorithm | |
CN104659812B (en) | A kind of many microgrid control method for coordinating based on PREDICTIVE CONTROL | |
CN110445186B (en) | Self-synchronizing microgrid control system and secondary frequency modulation control method | |
CN106253288A (en) | A kind of optimal load flow algorithm containing THE UPFC based on automatic differential | |
CN105226649A (en) | A kind of provincial power network power generation dispatching optimization method improved based on bus load prediction | |
Sun et al. | Multi-energy flow calculation method for we-energy based energy internet | |
CN104466979A (en) | Transmission grid reactive power optimization regulation and control system and method | |
CN110323779B (en) | Method and system for dynamically aggregating power of distributed power generation and energy storage device |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |