CN105356447B - A kind of electrical interconnection integrated energy system Analysis of Steady-state Security Region method - Google Patents
A kind of electrical interconnection integrated energy system Analysis of Steady-state Security Region method Download PDFInfo
- Publication number
- CN105356447B CN105356447B CN201510707513.0A CN201510707513A CN105356447B CN 105356447 B CN105356447 B CN 105356447B CN 201510707513 A CN201510707513 A CN 201510707513A CN 105356447 B CN105356447 B CN 105356447B
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- mtd
- gas
- msubsup
- 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
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- 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
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Abstract
The invention discloses a kind of electrical interconnection integrated energy system Analysis of Steady-state Security Region method.Electrical interconnection integrated energy system steady state energy flow model is initially set up；Definition and model referring next to power system static security domain, propose the concept of electrical interconnection integrated energy system Steady State Security Region；Then by repeating energy stream calculation, one group of border operating point on security domain boundaries is obtained；Linear Hyperplane fit is finally based on, the Steady State Security Region of electrical interconnection integrated energy system can be obtained.The present invention carries security domain models and assesses and provide theoretical foundation the real time execution scheduling for electrical interconnection integrated energy system and safety on line.
Description
Technical field
The present invention relates to a kind of electric-gas to interconnect integrated energy system Analysis of Steady-state Security Region method, belongs to comprehensive energy system
Safety on line of uniting is analyzed, assessment technology field.
Background technology
Compared to fired power generating unit, gas turbine group carbon emission amount is few, dynamic response is quick, the construction period is short, therefore in recent years
Carry out gas turbine power generation proportion to be stepped up.Correspondingly, power system and the coupling of natural gas system also constantly increase.The opposing party
Face, electricity turn advantage of the gas technology in terms of extensive energy storage, are expected to realize the two-way flow of power network and gas network energy stream.Power train
The coupling interconnection of system and natural gas system is expected to promote in terms of two the large-scale grid connection of new energy：1) gas turbine group is quick
Dynamic response characteristic can be used for the fluctuation for stabilizing new energy output；2) new energy that can not be dissolved for power network, can be turned by electricity
This portion of energy is converted into natural gas by gas technology, is stored in natural gas network.Therefore, electric-gas interconnection integrated energy system
It is the only way which must be passed that human society changes to low-carbon, sustainable energy system.
However, the safe operation that the depth between power network and gas net is coupled to the two brings potential risks.Such as
The U.S. in 2011 stops the supple of gas or steam greatly (electricity) event, and as shown in Figure 2, the event has its source in two aspects to its Evolution Mechanism：1) pole
Hold weather；2) power network couples with the depth of gas net.Contract can be interrupted by typically being signed due to gas turbine, when gas transmission resistance occurs for gas net
When plug or failure, gas turbine group will be removed first；In this case, when gas turbine power generation proportion is larger, great Liang electricity
The missing in source will cause large-scale power outage so that natural gas pressurizing point loses firm power supply, causes natural gas grid
Network gas transmission ability further reduces, i.e. accident spreads to gas net from power network in turn after gas net spreads to power network, causes thing
Therefore grow in intensity.
In fact, when power network couples with gas net depth, the generating proportion of gas turbine, which depends not only on itself installation, to be held
Amount, the security constraint of natural gas network is more dependent upon, and traditional power system static security domain is then difficult to embody this point, i.e., it is electric
Force system static security field result is excessively optimistic.
The content of the invention：
Goal of the invention：The technical problems to be solved by the invention are that solve existing power system static security domain to be difficult to embody
The problem of safe field result of the security constraint of natural gas network, i.e. power system static is excessively optimistic.
Technical scheme：The present invention to achieve the above object, adopts the following technical scheme that：
A kind of electric-gas interconnects integrated energy system Analysis of Steady-state Security Region method, comprises the following steps：
1) power system mesomeric state energy flow model is established：
For the branch road l of connecting node i, j, its branch power is represented by：
In formula：PlWith QlRespectively branch road active power and reactive power；V and θ are respectively node voltage amplitude and phase angle,
θij=θi-θj；gijWith bijRespectively branch road conductance and susceptance；gsh,iWith bsh,iRespectively conductance and susceptance over the ground；
Meanwhile each node flows into power and is necessarily equal to flow out power：
TGPG-TLPL=AePl；
TGQG-TLQL=AeQl；
In formula：TGFor grid nodes-generator incidence matrix, TLFor grid nodes-load incidence matrix, AeFor power network section
Point-branch road incidence matrix；PG、QGActive power, reactive power, P are injected for generatorL、QLFor the active power of load absorption,
Reactive power；
2) natural gas system steady state energy flow model is established：
For connecting node m and n pipeline l ', the flow for flowing through the pipeline is
In formula：Fl′For pipeline l ' flows；π is node pressure；kl′For the constant related to the pipeline l ' pressure losses；
Based on the balance of each node flow of natural gas system, can obtain：
TSFS-TDFD=AgFl′；
In formula：TSFor gas net node-source of the gas incidence matrix, TDFor gas net node-load incidence matrix, AgFor gas net node-
Pipeline incidence matrix；FSFlow, F are injected for source of the gasDFor gas load flow；
3) gas turbine couples
The gas quantity of gas turbine consumption is in following non-linear relation with its active output：
In formula：GgFor gas turbine set；It is active for gas turbine output,The air-flow absorbed for gas turbine
Amount；
4) power system static security domain
Power system static security domain is the one group of power note for meeting direction of energy constraint and network static security constraint
Enter collection：
In formula：ye=(PG,PL,QL), xe=(VL,θ)；Re(V under being constrained for static networkL, θ) and feasible zone in space；fe
For non-linear trend function；
5) electric-gas interconnection integrated energy system Steady State Security Region
With reference to the definition of power system static security domain, electric-gas interconnection integrated energy system Steady State Security Region may be defined as：
Meet that one group of energy stream injection of power network and the constraint of gas network energy stream and Static Security Constraints is gathered：
In formula：Y=(PG,PL,QL,FS,FD) injected for energy stream；X=(VL,θ,πD)；R is power network and gas net static security
The feasible zone in the lower x spaces of constraint；F is nonlinear energy flow equation；
6) operating point of security domain boundaries
To keep known quantity equal with unknown quantity quantity, it is known that amountIn some elements (it is assumed that) it is set to unknown
Amount, correspondingly, power system energy stream is equations turned to be：
In formula：flimitConstrained corresponding to the crucial of equality constraint is converted into；
Vector is injected for revised power；
By solving revised energy flow equation, an operating point on security domain boundaries can be obtained；By adjusting P 'G's
Size, one group of border operating point can be obtained；
7) Hyperplane fit
It is assumed that security domain boundaries are hyperplane, then power system security domain border is represented by：
In formula：CBeFor power system key constraint set；ηj,iFor j-th of crucial constraint centralized node i bill of power system
The hyperplane coefficient of machine, PG, iContributed for node i generated power；
Similarly, natural gas system security domain boundaries are：
In formula：CBgFor natural gas system key constraint set；τj,iFor j-th of crucial constraint centralized node i of natural gas system
The hyperplane coefficient of source of the gas；
Finally, electric-gas interconnection integrated energy system Steady State Security Region Ω is represented by：
As optimization, in the step 3) gas turbine in power network equivalent to power supply, while equivalent to gas in gas net
Load；Therefore gas turbine is connected to power network and gas net.
As optimization, the border of security domain is converted into equation about corresponding to the crucial constraint of some inequality in the step 6)
Beam, the equation constraint can be added in energy flow equation.
As optimization, the hyperplane coefficient η on electric-gas interconnection integrated energy system Steady State Security Region border in the step 7)
It can be obtained with τ using least square fitting.
Beneficial effect：The present invention is compared with prior art：When power network couples with gas net depth, the generating ratio of gas turbine
Itself installed capacity is depended not only on again, is more dependent upon the security constraint of natural gas network, and traditional power system static is pacified
Universe is then difficult to embody this point, i.e. the safe field result of power system static is excessively optimistic.The present invention proposes the comprehensive energy of electric-gas interconnection
Source static system security domain models, while count and the Static Security Constraints of power network and gas net, so as to be the comprehensive energy of electric-gas interconnection
Source system safety on line is analyzed, and the prevention and control before accident and the Corrective control after accident are provided fundamental basis.
Brief description of the drawings：
Fig. 1 is the schematic flow sheet of the present invention；
Fig. 2 is that the U.S. in 2011 stops the supple of gas or steam greatly (electricity) event Evolution Mechanism；
Fig. 3 is NGS5 node system schematic diagrames；
Fig. 4 is Steady State Security Region border schematic diagram；
Fig. 5 is that Steady State Security Region visualizes schematic diagram.
Embodiment：
The techniqueflow of 1~5 pair of invention is described in detail below in conjunction with the accompanying drawings：
Establish electric-gas interconnection integrated energy system steady state energy flow model
Establish power system mesomeric state energy flow model：
It is assumed that power system contains ne+ 1 node, beBar circuit；Definition node 0 is balance nodes, G={ 1,2, L, negBe
PV node, L={ neg+1,neg+2,L,neIt is PQ nodes, BeFor sets of lines.
For the branch road l of connecting node i, j, its branch power is represented by：
In formula：PlWith QlRespectively branch road active power and reactive power；V and θ are respectively node voltage amplitude and phase angle,
θij=θi-θj；gijWith bijRespectively branch road conductance and susceptance；gsh,iWith bsh,iRespectively conductance and susceptance over the ground.
Meanwhile each node flows into power and is necessarily equal to flow out power：
TGPG-TLPL=AePl；
TGQG-TLQL=AeQl；
In formula：TGFor grid nodes-generator incidence matrix, TLFor grid nodes-load incidence matrix, AeFor power network section
Point-branch road incidence matrix.PG、QGActive power, reactive power, P are injected for generatorL、QLFor the active power of load absorption,
Reactive power.
Establish natural gas system steady state energy flow model
Firstly the need of explanation, natural gas system can be stated by transient Model and steady-state model.Wherein transient Model is more
To be specific, but it is quite big by distributed constant and time-varying state variable description, computation complexity；Comparatively speaking, steady-state model is complicated
Spend low, and precision is in tolerance interval.Therefore steady state energy flow model is selected in the modeling to natural gas system herein.
It is assumed that natural gas system contains ngIndividual node and bgBar circuit；Define W={ 1,2, L, ngsIt is pressure known node collection,
H={ ngs+1,ngs+2,L,ngIt is that flow injects known node collection；For set H, S={ n are definedgs+1,ngs+2,L,ngs1Be
Source of the gas set of node, D={ ngs1+1,ngs1+2,L,ngsIt is load bus collection, H=S ∪ D.
For connecting node m and n pipeline l ', the flow for flowing through the pipeline is
In formula：Fl′For pipeline l ' flows；π is node pressure；kl′For the constant related to the pipeline l ' pressure losses.
Based on the balance of each node flow of natural gas system, can obtain：
TSFS-TDFD=AgFl′；
In formula：TSFor gas net node-source of the gas incidence matrix, TDFor gas net node-load incidence matrix, AgFor gas net node-
Pipeline incidence matrix；FSFlow, F are injected for source of the gasDFor gas load flow.
Calculate gas turbine coupling
Gas turbine in power network equivalent to power supply, while equivalent to gas load in gas net.Therefore gas turbine connection
Power network and gas net.The gas quantity of gas turbine consumption is in following non-linear relation with its active output：
In formula：GgFor gas turbine set；It is active for gas turbine output,The air-flow absorbed for gas turbine
Amount.
Calculate electric-gas interconnection integrated energy system Steady State Security Region
Calculate power system static security domain：
Power system static security domain is defined as one group of work(for meeting direction of energy constraint and network static security constraint
Rate injection collection.
Under power system static security constraint, (VL, θ) and the feasible zone in space is：
In formula：RV、Rl、AndRespectively voltage magnitude constraint, tributary capacity constraint, generated power with it is idle
(V under constraintL, θ) space feasible zone；ReFor (the V under all power system static security constraintsL, θ) space feasible zone；Sl
For branch road complex power,
Then, power system static security domain is represented by：
In formula：ye=(PG,PL,QL), xe=(VL,θ)；feFor non-linear trend function.
Calculate electric-gas interconnection integrated energy system Steady State Security Region
With reference to the definition of power system static security domain, electric-gas interconnection integrated energy system Steady State Security Region may be defined as：
Meet that one group of energy stream injection of power network and the constraint of gas network energy stream and Static Security Constraints is gathered.
Under natural gas system Static Security Constraints, πDThe feasible zone in space is
In formula：RFWith Rππ respectively under source of the gas capacity-constrained, node pressure constraintDThe feasible zone in space；RgIt is all
The lower π of natural gas system constraintDThe feasible zone in space.
With reference to power system static security domain models, electric-gas interconnection integrated energy system Steady State Security Region is represented by：
In formula：Y=(PG,PL,QL,FS,FD) injected for energy stream；X=(VL,θ,πD)；R=Re×Rg, it is power network and gas net
The feasible zone in x spaces under Static Security Constraints；F is nonlinear energy flow equation.
Calculate the security domain of decision space
The task of system coordinator is under given system structure and load, by the output for adjusting generator and source of the gas
Reaching ensures the purpose of the safe efficient operation of system.Under this application background, pass through the absorption of fixed electric load and gas load
Energy, the security domain of electric-gas interconnection integrated energy system can be reduced to security domain in decision space：
Calculate the operating point of security domain boundaries
The border of security domain corresponds to the crucial constraint of some inequality and is converted into equality constraint, and the equation constraint can be added to
In energy flow equation.By taking power system as an example, to keep known quantity equal with unknown quantity quantity, it is known that amountIn some member
Element (it is assumed that) unknown quantity is set to, correspondingly, power system energy stream is equations turned to be：
In formula：flimitConstrained corresponding to the crucial of equality constraint is converted into；
Vector is injected for revised power.
By solving revised energy flow equation, an operating point on security domain boundaries can be obtained.By adjusting P 'G's
Size, one group of border operating point can be obtained.
Carry out Hyperplane fit：
It is assumed that security domain boundaries are hyperplane, then power system security domain border is represented by：
In formula：CBeFor power system key constraint set；η is CBeCorresponding hyperplane coefficient.
Similarly, natural gas system security domain boundaries are：
In formula：CBgFor natural gas system key constraint set；τ is CBgCorresponding hyperplane coefficient.
It should be noted that above-mentioned natural gas security domain boundaries meter and the natural gas of gas turbine group consumptionEnter one
Step, by the active output of gas turbine groupInstead of its gas consumptionThen natural gas system Steady State Security Region border
It is converted into：
Meanwhile the hyperplane coefficient η and τ on electric-gas interconnection integrated energy system Steady State Security Region border can use a most young waiter in a wineshop or an inn
Multiplication is fitted to obtain.
Border based on power system Yu natural gas system Steady State Security Region, the static peace of electric-gas interconnection integrated energy system
Universe Ω is represented by：
Sample calculation analysis
The example that the present invention tests is to be saved by power system example IEEE39 node systems and natural gas system example NGS5
Dot system (as shown in annex Fig. 3) is formed.It is assumed that IEEE39 nodes interior joint 39 and 30 connected generator of node are combustion gas wheel
Machine, and be connected respectively with the node 2 in NGS5 nodes and node 5.
Based on least square fitting, IEEE39 node systems interior joint 37 and 30 running fire motors and NGS5 nodes system
System interior joint 4 connects the Steady State Security Region border of air accumulator as shown in annex Fig. 4.Define least square fitting error ε：
ε=| | (PG,FS)EF- (PG,FS)LS||/||(PG,FS)EF| | × 100%
In formula：(PG,FS)EFFor the border operating point being calculated according to repeated power flow；(PG,FS)LSIntend for least square method
Close obtained border operating point.
The maximum and average value of error ε are respectively 3.83% and 0.82%, and the error precision disclosure satisfy that Practical
Demand.
Based on the common factor of each Steady State Security Region border (hyperplane), Steady State Security Region is visualized as shown in annex Fig. 5.
In addition, hyperplane coefficient η embodies sensitivity of the crucial constraint to generator output.Pressed for natural gas node 5
Force constraint, the hyperplane coefficient of generator 30 are approximately 3 times of generator 37.Therefore, when substantial amounts of electric energy is needed by generator
During 30 and 37 supply, from natural gas system safety perspective, it should consider to be powered by generator 37 first.
Claims (4)
1. a kind of electric-gas interconnects integrated energy system Analysis of Steady-state Security Region method, it is characterised in that：Comprise the following steps：
1) power system mesomeric state energy flow model is established：
For the branch road l of connecting node i, j, its branch power is represented by：
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mi>l</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>g</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>g</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msubsup>
<mi>V</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msub>
<mi>g</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>V</mi>
<mi>i</mi>
</msub>
<msub>
<mi>V</mi>
<mi>j</mi>
</msub>
<msub>
<mi>cos&theta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>b</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>V</mi>
<mi>i</mi>
</msub>
<msub>
<mi>V</mi>
<mi>j</mi>
</msub>
<msub>
<mi>sin&theta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mi>l</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>b</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msubsup>
<mi>V</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>b</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>V</mi>
<mi>i</mi>
</msub>
<msub>
<mi>V</mi>
<mi>j</mi>
</msub>
<msub>
<mi>cos&theta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>g</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>V</mi>
<mi>i</mi>
</msub>
<msub>
<mi>V</mi>
<mi>j</mi>
</msub>
<msub>
<mi>sin&theta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mo>&ForAll;</mo>
<mi>l</mi>
<mo>&Element;</mo>
<msub>
<mi>B</mi>
<mi>e</mi>
</msub>
<mo>;</mo>
</mrow>
In formula：PlWith QlRespectively branch road active power and reactive power；V and θ are respectively node voltage amplitude and phase angle, θij=
θi-θj；gijWith bijRespectively branch road conductance and susceptance；gsh,iWith bsh,iRespectively conductance and susceptance over the ground, BeFor sets of lines；
Meanwhile each node flows into power and is necessarily equal to flow out power：
TGPG-TLPL=AePl；
TGQG-TLQL=AeQl；
In formula：TGFor grid nodes-generator incidence matrix, TLFor grid nodes-load incidence matrix, AeFor grid nodes-
Road incidence matrix；PG、QGActive power, reactive power, P are injected for generatorL、QLFor the active power of load absorption, idle work(
Rate；
2) natural gas system steady state energy flow model is established：
For connecting node m and n pipeline l ', the flow for flowing through the pipeline is
<mrow>
<msub>
<mi>F</mi>
<msup>
<mi>l</mi>
<mo>&prime;</mo>
</msup>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<msup>
<mi>l</mi>
<mo>&prime;</mo>
</msup>
</msub>
<msub>
<mi>s</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<msqrt>
<mrow>
<msubsup>
<mi>&pi;</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>&pi;</mi>
<mi>n</mi>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>s</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mn>1</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&pi;</mi>
<mi>m</mi>
</msub>
<mo>&GreaterEqual;</mo>
<msub>
<mi>&pi;</mi>
<mi>n</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&pi;</mi>
<mi>m</mi>
</msub>
<mo><</mo>
<msub>
<mi>&pi;</mi>
<mi>n</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
In formula：Fl′For pipeline l ' flows；π is node pressure；kl′For the constant related to the pipeline l ' pressure losses；
Based on the balance of each node flow of natural gas system, can obtain：
TSFS-TDFD=AgFl′；
In formula：TSFor gas net node-source of the gas incidence matrix, TDFor gas net node-load incidence matrix, AgFor gas net node-pipeline
Incidence matrix；FSFlow, F are injected for source of the gasDFor gas load flow；
3) gas turbine couples
The gas quantity of gas turbine consumption is in following non-linear relation with its active output：
<mrow>
<msub>
<mi>F</mi>
<mrow>
<msub>
<mi>G</mi>
<mi>g</mi>
</msub>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
<msubsup>
<mi>P</mi>
<mrow>
<msub>
<mi>G</mi>
<mi>g</mi>
</msub>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>&beta;</mi>
<mi>i</mi>
</msub>
<msub>
<mi>P</mi>
<mrow>
<msub>
<mi>G</mi>
<mi>g</mi>
</msub>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&gamma;</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<mi>i</mi>
<mo>&Element;</mo>
<msub>
<mi>G</mi>
<mi>g</mi>
</msub>
<mo>;</mo>
</mrow>
In formula：GgFor gas turbine set；It is active for gas turbine output,The throughput absorbed for gas turbine；
4) power system static security domain
Power system static security domain is the one group of power injection collection for meeting direction of energy constraint and network static security constraint：
<mrow>
<msub>
<mi>&Omega;</mi>
<mi>e</mi>
</msub>
<mo>:</mo>
<mo>=</mo>
<mo>{</mo>
<msub>
<mi>y</mi>
<mi>e</mi>
</msub>
<mo>:</mo>
<mo>&Exists;</mo>
<msub>
<mi>x</mi>
<mi>e</mi>
</msub>
<mo>&Element;</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>&DoubleRightArrow;</mo>
<msub>
<mi>f</mi>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>e</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>y</mi>
<mi>e</mi>
</msub>
<mo>}</mo>
<mo>;</mo>
</mrow>
In formula：ye=(PG,PL,QL), xe=(VL,θ)；Re(V under being constrained for static networkL, θ) and feasible zone in space；feTo be non-
Linear trend function；
5) electric-gas interconnection integrated energy system Steady State Security Region
With reference to the definition of power system static security domain, electric-gas interconnection integrated energy system Steady State Security Region may be defined as：Meet
Power network is constrained with gas network energy stream and the injection of one group of energy stream of Static Security Constraints is gathered：
<mrow>
<mi>&Omega;</mi>
<mo>:</mo>
<mo>=</mo>
<mo>{</mo>
<mi>y</mi>
<mo>:</mo>
<mo>&Exists;</mo>
<mi>x</mi>
<mo>&Element;</mo>
<mi>R</mi>
<mo>&DoubleRightArrow;</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>y</mi>
<mo>}</mo>
<mo>;</mo>
</mrow>
In formula：Y=(PG,PL,QL,FS,FD) injected for energy stream；X=(VL,θ,πD)；R is power network and gas net Static Security Constraints
The feasible zone in lower x spaces；F is nonlinear energy flow equation；
6) operating point of security domain boundaries
To keep known quantity equal with unknown quantity quantity, it is known that amountIn some elements (it is assumed that) unknown quantity is set to,
Correspondingly, power system energy stream is equations turned is：
<mrow>
<mo>{</mo>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>y</mi>
<mi>e</mi>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>f</mi>
<mi>e</mi>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>e</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>y</mi>
<mi>e</mi>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>lim</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>,</mo>
<msubsup>
<mi>P</mi>
<msub>
<mi>G</mi>
<mi>g</mi>
</msub>
<mo>&prime;</mo>
</msubsup>
<mo>,</mo>
<msub>
<mi>P</mi>
<mi>L</mi>
</msub>
<mo>,</mo>
<msub>
<mi>Q</mi>
<mi>L</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>;</mo>
</mrow>
</mrow>
In formula：flimitConstrained corresponding to the crucial of equality constraint is converted into；To repair
Power injection vector after just；
By solving revised energy flow equation, an operating point on security domain boundaries can be obtained；By adjusting P 'GSize,
One group of border operating point can be obtained；
7) Hyperplane fit
It is assumed that security domain boundaries are hyperplane, then power system security domain border is represented by：
<mrow>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>&Element;</mo>
<mi>G</mi>
</mrow>
</munder>
<msub>
<mi>&eta;</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mi>j</mi>
<mo>&Element;</mo>
<msub>
<mi>CB</mi>
<mi>e</mi>
</msub>
<mo>;</mo>
</mrow>
In formula：CBeFor power system key constraint set；ηj,iFor the super of j-th of crucial constraint centralized node i bill machine of power system
Floor coefficient, PG, iContributed for node i generated power；
Similarly, natural gas system security domain boundaries are：
<mrow>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>&Element;</mo>
<mi>S</mi>
</mrow>
</munder>
<msub>
<mi>&tau;</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>F</mi>
<mrow>
<mi>S</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>&Element;</mo>
<msub>
<mi>G</mi>
<mi>g</mi>
</msub>
</mrow>
</munder>
<msub>
<mi>&eta;</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>F</mi>
<mrow>
<msub>
<mi>G</mi>
<mi>g</mi>
</msub>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mi>j</mi>
<mo>&Element;</mo>
<msub>
<mi>CB</mi>
<mi>g</mi>
</msub>
<mo>;</mo>
</mrow>
In formula：CBgFor natural gas system key constraint set；τj,iFor j-th of crucial constraint centralized node i source of the gas of natural gas system
Hyperplane coefficient；
Finally, electric-gas interconnection integrated energy system Steady State Security Region Ω is represented by：
<mrow>
<mi>&Omega;</mi>
<mo>:</mo>
<mo>=</mo>
<mo>{</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mi>G</mi>
</msub>
<mo>,</mo>
<msub>
<mi>F</mi>
<mi>S</mi>
</msub>
<mo>)</mo>
</mrow>
<mfenced open = "|" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>&Element;</mo>
<msub>
<mi>G</mi>
<mi>e</mi>
</msub>
</mrow>
</munder>
<msub>
<mi>&eta;</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>&le;</mo>
<mn>1</mn>
<mo>,</mo>
<mi>j</mi>
<mo>&Element;</mo>
<msub>
<mi>CB</mi>
<mi>e</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>&Element;</mo>
<msub>
<mi>G</mi>
<mi>g</mi>
</msub>
</mrow>
</munder>
<msub>
<mi>&eta;</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
<msubsup>
<mi>P</mi>
<mrow>
<msub>
<mi>G</mi>
<mi>g</mi>
</msub>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>&beta;</mi>
<mi>i</mi>
</msub>
<msub>
<mi>P</mi>
<mrow>
<msub>
<mi>G</mi>
<mi>g</mi>
</msub>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&gamma;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>&Element;</mo>
<mi>S</mi>
</mrow>
</munder>
<msub>
<mi>&tau;</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>F</mi>
<mrow>
<mi>S</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>&le;</mo>
<mn>1</mn>
<mo>,</mo>
<mi>j</mi>
<mo>&Element;</mo>
<msub>
<mi>CB</mi>
<mi>g</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>P</mi>
<mi>G</mi>
<mi>min</mi>
</msubsup>
<mo>&le;</mo>
<msub>
<mi>P</mi>
<mi>G</mi>
</msub>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mi>G</mi>
<mi>max</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>F</mi>
<mi>S</mi>
<mi>min</mi>
</msubsup>
<mo>&le;</mo>
<msub>
<mi>F</mi>
<mi>S</mi>
</msub>
<mo>&le;</mo>
<msubsup>
<mi>F</mi>
<mi>S</mi>
<mi>max</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>}</mo>
<mo>.</mo>
</mrow>
2. electric-gas according to claim 1 interconnects integrated energy system Analysis of Steady-state Security Region method, it is characterised in that：
In the step 3) gas turbine in power network equivalent to power supply, while equivalent to gas load in gas net；Therefore gas turbine
It is connected to power network and gas net.
3. electric-gas according to claim 1 interconnects integrated energy system Analysis of Steady-state Security Region method, it is characterised in that：
The border of security domain is converted into equality constraint corresponding to the crucial constraint of some inequality in the step 6), and the equation constraint can add
Enter into energy flow equation.
4. electric-gas according to claim 1 interconnects integrated energy system Analysis of Steady-state Security Region method, it is characterised in that：
The hyperplane coefficient η and τ on electric-gas interconnection integrated energy system Steady State Security Region border can use least square in the step 7)
Method is fitted to obtain.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510707513.0A CN105356447B (en) | 2015-10-27 | 2015-10-27 | A kind of electrical interconnection integrated energy system Analysis of Steady-state Security Region method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510707513.0A CN105356447B (en) | 2015-10-27 | 2015-10-27 | A kind of electrical interconnection integrated energy system Analysis of Steady-state Security Region method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105356447A CN105356447A (en) | 2016-02-24 |
| CN105356447B true CN105356447B (en) | 2018-01-19 |
Family
ID=55332356
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510707513.0A Active CN105356447B (en) | 2015-10-27 | 2015-10-27 | A kind of electrical interconnection integrated energy system Analysis of Steady-state Security Region method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105356447B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108964041A (en) * | 2018-07-23 | 2018-12-07 | 燕山大学 | A kind of control method of electricity-interconnection integrated energy system peak load deliverability |
Families Citing this family (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106067677B (en) * | 2016-05-27 | 2018-07-06 | 清华大学 | A kind of electric-thermal couples multipotency streaming system static security analysis method |
| CN105939029B (en) * | 2016-06-21 | 2018-06-19 | 中国南方电网有限责任公司电网技术研究中心 | Integrated energy system programme acquisition methods and system |
| CN106208157B (en) * | 2016-07-19 | 2018-07-24 | 河海大学 | The electric-gas for turning gas based on electricity interconnects integrated energy system peak load shifting method |
| CN105978037B (en) * | 2016-08-03 | 2018-08-28 | 河海大学 | Multi-period electrical interconnection system optimal load flow computational methods based on forecasting wind speed |
| CN106777708B (en) * | 2016-12-21 | 2021-02-05 | 天津大学 | Steady state analysis method of electric power-natural gas regional comprehensive energy system |
| CN107194055B (en) * | 2017-05-17 | 2020-04-21 | 重庆大学 | Electric-gas interconnection system reliability modeling considering electric gas conversion device and evaluation method thereof |
| CN107276068B (en) * | 2017-05-24 | 2019-07-19 | 国家电网公司 | A kind of power grid multiple failure Steady State Security Region appraisal procedure based on the state addition method |
| CN107832959B (en) * | 2017-11-16 | 2021-12-03 | 广东电网有限责任公司电力调度控制中心 | Voltage stability evaluation method considering load characteristics and power supply constraints |
| CN108416507B (en) * | 2018-02-07 | 2021-08-20 | 天津大学 | Static sensitivity analysis method for electricity-gas coupling comprehensive energy system |
| CN109031952B (en) * | 2018-07-18 | 2020-06-12 | 河海大学 | Hybrid control method for electricity-gas interconnection comprehensive energy system |
| CN108964039A (en) * | 2018-07-23 | 2018-12-07 | 长沙理工大学 | A kind of electric-gas coupled system steady stability of wind power plant containing DFIG domain tangent plane method for solving |
| CN109217299A (en) * | 2018-09-29 | 2019-01-15 | 河海大学 | A method of electrical interconnection integrated energy system optimal energy stream is solved based on second order cone optimization algorithm |
| CN109742757A (en) * | 2019-03-11 | 2019-05-10 | 长沙理工大学 | A kind of generation method of integrated energy system load security domain |
| CN110543695B (en) * | 2019-08-14 | 2023-08-01 | 天津大学 | Electric-gas coupling comprehensive energy system feasible region calculation method |
| CN111416340B (en) * | 2020-03-10 | 2022-07-19 | 天津大学 | Regional comprehensive energy system optimization control method based on practical security domain |
| CN112434402B (en) * | 2020-10-22 | 2022-08-12 | 天津大学 | Interval practical safety domain modeling method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6988025B2 (en) * | 2000-11-28 | 2006-01-17 | Power Measurement Ltd. | System and method for implementing XML on an energy management device |
| CN101534010A (en) * | 2009-04-14 | 2009-09-16 | 天津大学 | Method for solving the local boundary of cut-set voltage stability regions on the basis of perturbation |
| CN102623988A (en) * | 2012-03-27 | 2012-08-01 | 清华大学 | Automatic generating method for day-ahead plan power flow in power grid |
| CN104734155A (en) * | 2015-04-14 | 2015-06-24 | 河海大学 | Method of acquiring available transfer capacity of electricity-gas interconnected energy system |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101489533B1 (en) * | 2010-11-29 | 2015-02-03 | 엘에스산전 주식회사 | Method for selecting an available transfer capability |
-
2015
- 2015-10-27 CN CN201510707513.0A patent/CN105356447B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6988025B2 (en) * | 2000-11-28 | 2006-01-17 | Power Measurement Ltd. | System and method for implementing XML on an energy management device |
| CN101534010A (en) * | 2009-04-14 | 2009-09-16 | 天津大学 | Method for solving the local boundary of cut-set voltage stability regions on the basis of perturbation |
| CN102623988A (en) * | 2012-03-27 | 2012-08-01 | 清华大学 | Automatic generating method for day-ahead plan power flow in power grid |
| CN104734155A (en) * | 2015-04-14 | 2015-06-24 | 河海大学 | Method of acquiring available transfer capacity of electricity-gas interconnected energy system |
Non-Patent Citations (1)
| Title |
|---|
| 含多类型分布式电源的主动配电网分布式三相状态估计;卫志农等;《电力系统自动化》;20150510;第39卷(第9期);第68-74页 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108964041A (en) * | 2018-07-23 | 2018-12-07 | 燕山大学 | A kind of control method of electricity-interconnection integrated energy system peak load deliverability |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105356447A (en) | 2016-02-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105356447B (en) | A kind of electrical interconnection integrated energy system Analysis of Steady-state Security Region method | |
| CN104734155B (en) | A kind of method obtaining electrical interconnection energy resource system available transmission capacity | |
| CN104734147B (en) | A kind of integrated energy system likelihood energy stream analysis method | |
| CN105046369B (en) | A kind of modeling of electrical combined hybrid system and Optimization Scheduling based on energy centre | |
| Beerten et al. | VSC MTDC systems with a distributed DC voltage control-A power flow approach | |
| CN106096269B (en) | The Interval Power Flow calculation method of natural gas grid in a kind of electric-gas coupled system | |
| CN109031952B (en) | Hybrid control method for electricity-gas interconnection comprehensive energy system | |
| Li et al. | Control strategy of battery state of charge for wind/battery hybrid power system | |
| CN107808216B (en) | Method for constructing comprehensive minimum optimization model for wind curtailment and electrical thermal load reduction of electricity-gas-heat interconnection system | |
| CN105005940A (en) | Correlation-considered GEPOPF calculation method | |
| CN105745809B (en) | Obtain the symmetry approach of the non-linear effective power flow of multiterminal element power network | |
| CN109376428B (en) | Reliability evaluation method, device, equipment and storage medium of comprehensive energy system | |
| CN105576657A (en) | Random fuzzy continuation power flow modeling and algorithm of active distribution network and transmission network integrated system under energy source interconnection | |
| CN106096777A (en) | A kind of Optimization Scheduling of electrical couplings multipotency streaming system | |
| CN107947245B (en) | Equivalent optimal power flow model construction method considering natural gas system constraint | |
| CN106026169A (en) | Decomposition-coordination optimization method based on multi-microgrid merging into power distribution network | |
| Chen et al. | Steady-state security regions of electricity-gas integrated energy systems | |
| CN104318326A (en) | Net rack optimizing model for improving renewable energy source acceptance ability | |
| US20230369852A1 (en) | Method for calculating optimal energy flow of integrated electricity-gas system based on sequential convex programming | |
| CN107122924A (en) | A kind of intelligent distribution system bearing capacity evaluation method | |
| CN113256152A (en) | Reliability evaluation method and system for gas-electricity coupling system based on Bayesian network | |
| CN106655194A (en) | UPFC-considered ATC calculation method of electric power system including wind power | |
| CN109742757A (en) | A kind of generation method of integrated energy system load security domain | |
| CN104899414B (en) | A kind of electrical interconnection energy resource system Calculation of Available Transfer Capability method | |
| CN103761680A (en) | Grid and provincial dispatching method and device for AC/DC interconnected large power grid with wind farm |
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 |