CN108416507A - The static sensitivity analysis method of integrated energy system is coupled towards electric-gas - Google Patents
The static sensitivity analysis method of integrated energy system is coupled towards electric-gas Download PDFInfo
- Publication number
- CN108416507A CN108416507A CN201810123783.0A CN201810123783A CN108416507A CN 108416507 A CN108416507 A CN 108416507A CN 201810123783 A CN201810123783 A CN 201810123783A CN 108416507 A CN108416507 A CN 108416507A
- Authority
- CN
- China
- Prior art keywords
- gas
- node
- power
- sensitivity
- pressure
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000010206 sensitivity analysis Methods 0.000 title claims abstract description 25
- 230000003068 static effect Effects 0.000 title claims abstract description 25
- 230000035945 sensitivity Effects 0.000 claims abstract description 61
- 239000000567 combustion gas Substances 0.000 claims abstract description 40
- 230000008878 coupling Effects 0.000 claims abstract description 26
- 238000010168 coupling process Methods 0.000 claims abstract description 26
- 238000005859 coupling reaction Methods 0.000 claims abstract description 26
- 239000007789 gas Substances 0.000 claims description 117
- 239000011159 matrix material Substances 0.000 claims description 38
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 30
- 239000003345 natural gas Substances 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 13
- 230000005611 electricity Effects 0.000 claims description 11
- 230000005540 biological transmission Effects 0.000 claims description 7
- 238000010276 construction Methods 0.000 claims description 3
- 240000002853 Nelumbo nucifera Species 0.000 claims description 2
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims description 2
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims description 2
- 230000009466 transformation Effects 0.000 claims description 2
- 230000017105 transposition Effects 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 abstract description 11
- 230000000694 effects Effects 0.000 abstract description 9
- 230000003993 interaction Effects 0.000 abstract description 8
- 238000002347 injection Methods 0.000 description 11
- 239000007924 injection Substances 0.000 description 11
- 238000010586 diagram Methods 0.000 description 10
- 239000002737 fuel gas Substances 0.000 description 8
- 230000008859 change Effects 0.000 description 5
- 230000012010 growth Effects 0.000 description 5
- 239000000243 solution Substances 0.000 description 5
- 238000011161 development Methods 0.000 description 3
- 230000018109 developmental process Effects 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 241000196324 Embryophyta Species 0.000 description 2
- 230000002411 adverse Effects 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000001172 regenerating effect Effects 0.000 description 2
- 238000011144 upstream manufacturing Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- 241000158500 Platanus racemosa Species 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000004378 air conditioning Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000010534 mechanism of action Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000005057 refrigeration Methods 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- 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
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- 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—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Tourism & Hospitality (AREA)
- Educational Administration (AREA)
- General Business, Economics & Management (AREA)
- Development Economics (AREA)
- Water Supply & Treatment (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a kind of static sensitivity analysis methods towards electrical couplings integrated energy system, for analyzing the mechanism of interaction effect between electric power and combustion gas energy supplying system.First, the present invention provides the Unified Power Flow model of electrical couplings integrated energy system；On this basis, the gaseous-pressure node injecting power sensitivity of electrical couplings integrated energy system is defined；Finally, in conjunction with the integrated energy system sensitivity index under typical scene, influence of the analysis grid nodes injecting power to gaseous-pressure positions the weak link of integrated energy system.Example shows that the present invention can provide auxiliary information, the safety of effective lifting system for the safe and stable operation of regional complex energy resource system.
Description
Technical field
The present invention relates to electric-gas coupling integrated energy system field more particularly to a kind of electricity-based on Unified Power Flow model
Gas couples the static sensitivity analysis method of integrated energy system.
Background technology
With increasingly deficient and environment the continuous deterioration of fossil energy, energy transition, which becomes, realizes that economic society is sustainable
The only way of development.And break each energy supplying system individually plan, independently operated existing pattern, realize energy system march toward it is more
It can merge, the completely new energy system of integrated complementary, be the important channel for pushing above-mentioned process.For this purpose,《Energy development " 13 "
Planning》Clearly proposing will " implementation provide multiple forms of energy to complement each other integrated optimization engineering " and " make overall planning electric power, combustion gas, heating power, cooling, water supply
The infrastructure such as piping lane build integrated energy supplying system ", advance the hair that various energy resources comprehensively utilize from State-level
Exhibition[1]。
Integrated energy system (Integrated Energy Systems, IES) be known as be multiple-energy-source comprehensive utilization master
One of realization means are wanted, it can pass through the scientific management and optimization between the difference such as the cold heat/electricity/gas/hydrogen grade energy in system
Scheduling improves high efficiency of energy utilization ratio, realizes the abundant consumption of regenerative resource[2].However, IES states target in realization
Meanwhile also bringing the globality risk of system safety operation.In face of the close-coupled of multiple-energy-source, between different energy supplying systems
" pulling one hair and move the whole body ", reciprocal effect to each other should be paid close attention to.On the one hand, the failure of gas burning system will conduct extremely
Electric system, and directly cause power outage, 8.15 large-scale blackouts that such as Taiwan in 2017 occurs rise because of natural gas
6 unit off-grids caused by supply discontinuity, accident cause 6,680,000 family users to have a power failure, and impacted population is more than the total people in Taiwan
The 85% of mouth[3].Similar event also appears in the U.S., Southern California of the U.S. in 2015, and Aliso Canyon natural gas leakings cause
Plant gas natural gas supply is insufficient, seriously affects the normal operation of local electric system[4].On the other hand, electric system
Adverse effect can also be transmitted to natural gas system, and jeopardize the safe operation of natural gas system.As the whole America regenerative resource is oozed
The raising of saturating rate, the frequent adjusting of the plant gas as main peak regulation resource lead to fluctuating widely for gas ductwork pressure,
Directly affect the gas transmission safety of natural gas system[5].Meanwhile the electric power link of IES often contains distributed generation resource (Distributed
Generation, DG), the interval fluctuation contributed will be conducted by the coupling link (such as Gas Generator Set) between the two to entire
IES adversely affects whole system[6].In view of this, the unfavorable factor of each energy supplying system couples comprehensive energy for electric-gas
The influence of source system (Integrated Electricity and Gas Systems, IEGS) safe operation must cause greatly
Concern.Without loss of generality, influence of the unfavorable factor for the safe operation of system is often depending on the safety of system weakness
Nargin is run, therefore, quickly positions the problem of weak link of IEGS is at current urgent need to resolve.
Document [7] establishes the steady-state analysis model of electric-gas coupled system, and the friendship of two systems is considered from energy flow angle
Interaction；Further, document [8] considers to propose electric-gas coupled system optimal energy flow model on the basis of the network operation constrains,
Decision assistant is provided for the Optimized Operation of system；Document [9] is based on Load flow calculation, and analyzing heat supply network load using point by point method increases
Add the influence for grid nodes voltage and heat supply network node heat supply temperature；Document [10] research shows that natural gas load it is continuous
Variation can carry out great influence to combustion gas network press belt in IEGS systems；Document [11] is based on energy hub model, using mixed
In the swim under Algorithm Analysis stable state each electric-gas energy supply network interaction characteristic.
However, the studies above lays particular emphasis on the Load flow calculation of IEGS, different energy supply network state (such as networks are not fully considered
Topology, pipeline configuration, load level etc.) interaction between integrated energy system different energy sources influence, and how to judge
The weak link of IEGS.
Invention content
The present invention provides a kind of static sensitivity analysis method coupling integrated energy system towards electric-gas, the present invention
Define gaseous-pressure-node injecting power sensitivity index of electric-gas coupling integrated energy system；In conjunction under typical scene
Integrated energy system sensitivity index analyzes influence of the grid nodes injecting power to gaseous-pressure；Position integrated energy system
Weak link provides decision-making foundation for integrated energy system operation control (the especially addressing of caisson), as detailed below to retouch
It states：
A kind of static sensitivity analysis method coupling integrated energy system towards electric-gas, the static sensitivity analysis
Method includes the following steps：
Establish by gas burning system model, distribution system model and the energy coupling link model group at IEGS Unified Power Flows
Model；
According to IEGS Unified Power Flow models, the Jacobian matrix of current stable operating point and current point of operation is obtained；
Gaseous-pressure-node injecting power sensitivity is calculated according to Jacobian matrix；Work(is injected to gaseous-pressure-node
Rate sensitivity carries out size sequence, and the big node set of sensitivity is the weak link of IEGS.
The IEGS Unified Power Flows model is specially：
In formula, x=[θ, V, p] indicates the state variable of IEGS, respectively power distribution network each node phase angle in addition to balance nodes,
Electric load node voltage, Gas Load node air pressure；U=[PPDS, QPDS, LNGS] indicate the disturbance variable of IEGS, respectively node
Inject active power, reactive power and Gas Load；A indicates independent parameter；PPDSActive power is injected for node；For power grid section
Point voltage；Y is node admittance matrix, and * is complex conjugate；ANGSFor the node-branch incident matrix of combustion gas network；It is light
The electromotive power output of type gas turbine；cgeFor transformation ratio.
It is described gaseous-pressure-node injecting power sensitivity is calculated according to Jacobian matrix to be specially：
According to Jacobian matrix, going out for gaseous-pressure-Gas Load sensitivity and lightweight gas turbine is calculated separately
Power-node injecting power sensitivity；
Pass through the output of gaseous-pressure-Gas Load sensitivity and lightweight gas turbine-node injecting power sensitivity
Calculate gaseous-pressure-node injecting power sensitivity.
The gaseous-pressure-node injecting power sensitivity is specially：
Wherein, PPDSActive power is injected for node；P is combustion gas node air pressure；LNGSFor Gas Load；Represent MT
Electromotive power output；For the gas consumption of MT.
Further, the gaseous-pressure-Gas Load sensitivity SggEqual to Jacobi in natural gas system Load Flow Solution
Matrix JggThe negative inverted：
In formula, ANGSFor the node-branch incident matrix of combustion gas network；fiFor pipeline i flows；∏iFor pipeline i pressure differences；
npFor number of tubes；T is matrix transposition；Diag () indicates construction diagonal matrix.
When specific implementation, if a certain node i injecting powerIt changes, has for power distribution network：
In formula,It is right for distribution network lossLocal derviation is asked to have：
In formula, Section 2 is Incremental Transmission Loss on the right of equation, represents network loss with the variation relation between node injecting power.
The advantageous effect of technical solution provided by the invention is：
1, thxe present method defines the indexs of integrated energy system air pressure safety, located the weak ring of integrated energy system
Section；
2, this method analyze under typical scene electric system in integrated energy system due to the fluctuation of load power or
It is the fluctuation that generation of electricity by new energy is contributed, causes to run the unstable mechanism of action of air pressure；
3, this method provides corresponding guidance for the selection of integrated energy system gas storage position.
Description of the drawings
Fig. 1 is a kind of flow chart for the static sensitivity analysis method coupling integrated energy system towards electric-gas；
Fig. 2 is that electric-gas couples integrated energy system schematic diagram；
Fig. 3 is IEGS example schematic diagrames；
Fig. 4 is the S of each node of power distribution networkeeSchematic diagram；
Fig. 5 is the S of each node of fuel gas networkggSchematic diagram；
Fig. 6 is gaseous-pressure-node injecting power sensitivity schematic diagram；
Fig. 7 is the air pressure schematic diagram under different load growth level；
Fig. 8 is the air pressure schematic diagram under different schemes.
Specific implementation mode
To make the object, technical solutions and advantages of the present invention clearer, embodiment of the present invention is made below further
It is described in detail on ground.
Sensitivity analysis technique is widely used to Operation of Electric Systems optimization, wherein the static analysis based on power flow equation
Method can analyze changing trend between the two by the differential relationship between different physical quantities, to judge system weakness
The information such as node, crucial branch[12]。
For this purpose, the embodiment of the present invention it is expected, by Sensitivity Analysis Method ripe in electric system, to propose for IEGS
The static sensitivity analysis method of integrated energy system is coupled towards electric-gas.This method considers the network-like of electric power-gas burning system
State, analysis gas pipe pressure deeply excavate different energy supply networks in IEGS with the interaction rule between electrical power Injection Space
Between reciprocal effect, judge influence system safety operation weakness zone, example show this method quickly and directly obtain IEGS
Weakness zone, provide decision-making foundation for the IEGS operation controls addressing of caisson (especially).
Embodiment 1
A kind of static sensitivity analysis method coupling integrated energy system towards electric-gas, referring to Fig. 1, the static state is sensitive
Degree analysis method includes the following steps：
101：Establish by gas burning system model, distribution system model and the energy coupling link model group at IEGS it is unified
Tide model；
102：According to IEGS Unified Power Flow models, the Jacobean matrix of current stable operating point and current point of operation is obtained
Battle array；
103：According to Jacobian matrix, gaseous-pressure-Gas Load sensitivity and lightweight gas turbine are calculated separately
Output-node injecting power sensitivity calculates gaseous-pressure-node injecting power sensitivity with this；
104：Size sequence is carried out to gaseous-pressure-node injecting power sensitivity, the larger node set of sensitivity is
For the weak link of IEGS.
In conclusion 101- steps 104 consider the network-like of electric power-gas burning system to the embodiment of the present invention through the above steps
State, for analysis gas pipe pressure with the interaction rule between electrical power Injection Space, this method is quick and directly obtains IEGS's
Weakness zone provides decision-making foundation for the operation control of IEGS.
Embodiment 2
The scheme in embodiment 1 is further introduced with reference to specific calculation formula, Fig. 2, specific example,
It is described below：
201：IEGS is modeled；
The embodiment of the present invention is to include combustion gas network (Natural Gas System, NGS), power distribution network (Power
Distribution System, PDS) and the IEGS of coupling link for verify the validity of this method.
Wherein, NGS is made of air source, gas pipeline, Gas Load, compressor, valve etc..Valve exists for controlling combustion gas
Circulation in pipeline or cut-off, it is assumed that valve is only in and opens or completely close completely two states, therefore the network of NGS is opened up
It flutters and is to determine.IEGS is connected with bulk power grid by distribution transformer and assumes that IEGS has signed the electric power supply of material with upstream power grid and closed
Together.Energy coupling link is the key that each energy supplying system interacts a ring in IEGS, has the characteristics that Coupling device is various.It is such as light-duty
Gas turbine (Micro-turbines, MT) realizes coupling of the gas burning system with power grid；Turn gas technology (Power to by electricity
Gas, P2G), electric energy can be converted into natural gas energy resource.In view of MT is widely used in IEGS, the embodiment of the present invention is arranged MT and makees
For the coupling link of IEGS.The electric-gas coupling integrated energy system that the embodiment of the present invention is studied is as shown in Figure 2.
The embodiment of the present invention assumes that IEGS has signed electric power supply contract with external bulk power grid, contract provision different periods
The electricity that power distribution network is interacted with bulk power grid, therefore imbalance power caused by all kinds of unfavorable factors is held by MT in operational process
Load, i.e. balance units of the MT as power distribution network.
Consider the Typical conductive scene of IEGS unfavorable factors：A large amount of uses of the refrigeration equipments such as air-conditioning cause after into summer
Electric load surge, at this time the power shortage in system supplied by MT, cause Gas Load variation, further result in combustion gas
The variation of pressure；Or due to causing MT's to go out fluctuation as the fluctuation of the DG short-time ratings of representative using wind-power electricity generation or photovoltaic, into
And cause the variation of gas combustion gas network air pressure.In short, the variation of distribution network arbitrary node injecting power can be by coupling link
Influence the air pressure level of fuel gas network.Under certain extreme scenes, or even the operation problem of gas burning system can be caused, shaped like gas pressure
Power fluctuate widely and combustion gas air pressure is out-of-limit.
Therefore can using combustion gas network air pressure level as one of important indicator of IEGS safe operations, and think work as it is unfavorable because
Air pressure declines the weak link that faster region is IEGS in IEGS when element acts on system.The weak link of acquisition system, and
It is a kind of Control Measure for improving IEGS air pressure safety margins by installing gas storage equipment in weak link.
1) gas burning system model
The model of gas burning system is mainly the steady state flow equation of pipeline, wherein and can be specifically divided into containing compressor
Pipeline and pipeline without compressor.There is pressure drops for pipe ends, especially the mesohigh in long range large capacity
Gas transmission network needs to configure compressor to be lifted the pressure of gas pipeline.The IEGS studied in view of the embodiment of the present invention
Belong to the scope of low pressure distribution network (0-75mbar), the embodiment of the present invention puts aside the influence of compressor[13]。
The node variable of gas burning system has injection gas flow and node air pressure, copies the node-classification of electric system, root
Node can be divided into pressure known node and injection flow known node according to known variables.In gas burning system, air source is balance
Node, air pressure is known and injects that flow is unknown, the air pressure of Gas Load is unknown and demand on gas known to.If gas burning system has ng
A node, npPipeline, node 1 are air source node, remaining is Gas Load node.For connecting node i, j (i, j=1,
2,…,ng) gas pipeline k, the steady state flow f in low-pressure fuel gas pipeline (0-75mbar)k(k=1,2 ..., np) can use down
Formula describes[10,13]：
Π=- (ANGS)Tp (3)
In formula, λkFor the friction coefficient of pipeline k；ΠkRepresent the pressure difference at the both ends pipeline k；DkAnd LkRepresent the straight of pipeline k
Diameter and length；G is the relative density of combustion gas；P represents combustion gas node air pressure；ANGSIt is associated with square for node-branch of combustion gas network
Battle array.Work as Πk>When 0, sk=1；Work as Πk<When 0, sk=-1.
The first and second law of analogy kirchhoff, combustion gas is flowed in a network should meet following two conditions：Inject certain
The gas flow of one node is equal with the gas flow for flowing out the node；To any one circuit of network, combustion gas is flowing over
The sum of pressure drop in journey is zero.Therefore, the flow behavior of combustion gas in a network can be described with following formula：
ANGSF=LNGS (4)
In formula, f is pipeline steady state flow column vector；LNGSFor Gas Load.
2) distribution system model
The embodiment of the present invention is assumed effectively to have solved the problems, such as the three-phase imbalance of power distribution network by commutation[14], therefore this hair
Bright embodiment, which ignores three-phase imbalance, influences Load flow calculation.If power distribution network has neA node, node 1 are balance nodes, node 2
To node 1+npvFor PV node (node of given active-power P and voltage V), remaining is PQ nodes (given active-power P and nothing
The node of work(power Q).The model of power distribution network is to reflect node power and node voltage, the node power side of relationship between phase angle
Journey[14]：
In formula, PiFor the injection active power of node i；QiFor the injection reactive power of node i；For in node i
Generator is sent out active, reactive power；For active, the reactive power of load in node i；Vi, VjFor node i, j
Voltage；Gij, BijFor node i, admittance Y between jijReal part, imaginary part；θijIt is node i with phase angle difference (i, the j between node j
=1,2 ..., ne)。
3) energy couples Link Model
The embodiment of the present invention assumes the output adjustment of MT in second grade[15], the reciprocal effect between considering different energy supply networks
When, ignore the dynamic characteristic of MT, considers influencing each other for network from energetic interaction angle.The input of MT is combustion gas, utilizes combustion gas
The high-temperature gas expansion work generated that burns exports electric energy in turn.Therefore the model of MT is to characterize its stable state input-output characteristic
Static models, MT models are as follows：
In formula,Represent the electromotive power output of MT；Represent the gas consumption of MT；cgeRepresent the conversion of MT
Coefficient.
202：IEGS static sensitivity analysis methods；The Load Flow Solution model of IEGS is segmented into Unified Power Flow model and divides
Vertical solving model[10], the two has no influence for the derivation of sensitivity analysis, and the embodiment of the present invention is using Unified Power Flow model as base
Plinth further defines gaseous-pressure-node injecting power sensitivity index, and passes through the thin of static sensitivity analyzing and positioning system
Weak link.
The embodiment of the present invention further defines gaseous-pressure-node injection work(on the basis of IEGS Unified Power Flow models
Rate sensitivity index, and pass through the weak link of static sensitivity analyzing and positioning system.
1) IEGS Unified Power Flows model；
For IEGS, the essence of Load flow calculation is that the mistake of system stable operation point is sought under given a series of conditions
Journey, Unified Power Flow model can be described as：
In formula, x=[θ, V, p] indicates the state variable of IEGS, respectively power distribution network each node phase angle in addition to balance nodes,
Electric load node voltage, Gas Load node air pressure；U=[PPDS, QPDS, LNGS] indicate the disturbance variable of IEGS, respectively node
Inject active power, reactive power and Gas Load, Gas Load；A indicates independent parameter, such as electric system node admittance, combustion
Gas grid topology, pipe parameter etc..PPDSActive power is injected for node；For grid nodes voltage；Y is node admittance square
Battle array, this inside formula (5)-(8) it is stated that, * is the mathematical operation symbol of complex conjugate.
Unified Power Flow model is solved by Newton-Laphson method, iteration form is as follows：
In formula, Δ x(k+1)For the variable quantity of the state variable of IEGS in kth+1 this iteration；x(k+1)For in kth+1 this iteration
The state variable of IEGS；Δ F is the deviation of power flow equation group；J is Jacobian matrix, be can be expressed as[7]：
Primary condition is given according to system, simultaneous formula (10)-(12) solve power flow equation group and obtain the current stabilization of system
Operating point.
Wherein, the element meaning in Jacobian matrix J is known to those skilled in the art, the embodiment of the present invention to this not
It repeats.
2) gaseous-pressure-node injecting power sensitivity matrix
Influence of the apparent unfavorable factor to different energy supply networks, being equal to analysis IEGS, (such as node injects in unfavorable factor
Changed power) under the action of system running state variation, and unfavorable factor generally can use system disturbance variable characterize, because
This IEGS static sensitivity analyzes universal expression：
According to the variation relation of different physical quantities in different energy supply networks, various informative sensitivity can be constructed and referred to
Mark.Reciprocal effect between two energy supply networks is conducted by MT, and the balancing machine that MT is power distribution network is arranged according to Section 1 Scene,
And it is connected to the node h in fuel gas network.Consider that grid nodes are injected work(by unfavorable factor Typical conductive scene, the embodiment of the present invention
Unfavorable factor of the fluctuation of rate as IEGS, state variable of the combustion gas air pressure as IEGS define gaseous-pressure-according to formula (13)
Node injecting power sensitivity matrix SgeIt is as follows：
Wherein,It is matrix SgeIn element, represent variation (i=1,2 ..., n of power distribution network node i injecting powere) right
Combustion gas node j (j=1,2 ..., ng) air pressure influence.
Define SggFor gaseous-pressure-Gas Load sensitivity, SeeAs follows, the root for MT outputs-node injecting power sensitivity
Know according to formula (9)For ceg, then SgeFor the product of this three.
It is matrix SggIn element, represent MT demand on gas variation for combustion gas node j air pressure changes influence (j=
1,2,…,ng)；It is matrix SeeElement, represent power distribution network node i injecting electric power variation for MT contribute variation shadow
It rings.
For fuel gas network, S known to analysis mode (1)-(4)ggEqual to Jacobian matrix J in natural gas system Load Flow SolutionggIt asks
Inverse negative：
In formula, function diag () indicates construction diagonal matrix.
For power distribution network, if a certain node i injecting powerIt changes, has for power distribution network：
In formula,For distribution network loss, it is contemplated that other node injecting powers are constant, and formula (19) is to node injecting powerLocal derviation is asked to have：
In formula, Section 2 is Incremental Transmission Loss on the right of equation, represents network loss with the variation relation between node injecting power.
Intuitively, it contributes for a lossless network MT identical with node injecting power variable quantity.However, for an actual electricity
Net, network topology, load level can influence the Incremental Transmission Loss of network, be handed over to influence the energy between each energy supply network
Mutually.And network loss is determined by each node voltage and phase angle[16], further deriving Incremental Transmission Loss has：
Notice power grid Jacobian matrix J in formula (12)eeIt inverts and is represented by：
Contrast (20) and (22), it is seen that Jacobian matrix JeeCan be to solve SeeCorresponding information, simultaneous formula are provided
(20)-(22) can solve See。
203：The realization step of gaseous-pressure-node injecting power sensitivity analysis.
Value it is bigger show combustion gas node j (j=1,2 ..., ng) air pressure change to power distribution network node i (i=1,2 ...,
ne) electrical power injection variation sensitivity.It can be seen that when unfavorable factor acts on IEGS, the big node air pressure change amplitude of sensitivity is most
Greatly, combustion gas network air pressure restrained boundary will be reached at first, can be regarded as the short slab of system safety operation.SgeConsider different confessions
The state of energy network, reflects influence of the disturbance of electrical power Injection Space to natural gas air pressure level, can be operation and control
Important information is provided.
It is worth noting that, the difference of the on-positions MT will influence the trend distribution of IEGS, therefore in meter sensitivity
It is required for first giving the deployed position of MT before matrix, changes the primary condition of Load flow calculation, form different stable operating points.
Sensitivity Analysis Method step is：
1) natural gas network state information is read, including：Natural gas line parameter, tube topology information form ANGS, read
The power distribution network network information is taken, Y is formed；
2) it solves power flow equation group according to formula (10)~(12) and obtains the current stable operating points of IEGS and current operation
The Jacobian matrix J of point；
3) according to Jacobian matrix J, S is calculated separately according to formula (16), (17), (20)~(22)gg、See；
4) gaseous-pressure-node injecting power sensitivity S is calculated according to formula (14)ge；
5) to gaseous-pressure-node injecting power sensitivity SgeLevel of sensitivity sequence is carried out, sensitivity number is larger
Node (specific numerical value selection according in practical application set, and generally takes preceding the 30% of sensitivity number size,
The embodiment of the present invention is without limitation) the set as weak link of IEGS.
In conclusion 201- steps 203 consider the network-like of electric power-gas burning system to the embodiment of the present invention through the above steps
State, for analysis gas pipe pressure with the interaction rule between electrical power Injection Space, this method is quick and directly obtains IEGS's
Weakness zone provides decision-making foundation for IEGS operation controls (addressing of such as caisson).
Embodiment 3
Below the feasibility of scheme in Examples 1 and 2 is verified with specific example, in conjunction with Fig. 3-Fig. 8, it is as detailed below to retouch
It states：
To verify the validity of this method, by taking a typical IEGS as an example, emulation explanation is carried out.As shown in figure 3, this hair
The IEGS examples of bright embodiment are by IEEE-33 Node power distribution systems and 11 node combustion gas networks of modification by MT is coupled
At[13,17], EBi and GBi indicate grid nodes and combustion gas node respectively.
Power distribution network is connected by EB1 with external bulk power grid, it is assumed that power distribution network is with the electrical power that external bulk power grid obtains
3500kW, EB1 is PV node at this time.The EB2 of the GB11 and power distribution network of MT connection combustion gas networks, as balance unit, therefore EB2
For the balance nodes of power distribution network.It is assumed that the warning air pressure of load bus (GB2~GB11) safe operation is in fuel gas network
The air pressure of 20mbar, air source node GB1 are 75mbar, cgeFor 6.1kW/m3.Combustion gas network data is shown in Table A1, A2.
1 gas pipeline data of Table A
2 Gas Load data of Table A
1) Calculation of Sensitivity and weak link are analyzed first, according to the operating point of the primary condition computing system of IEGS and refined
Than matrix J.S is calculated according to formula (20)~(22)eeThe results are shown in Figure 4, and block diagram size indicates MT output variations pair in figure
In the sensitivity of grid nodes injecting power variation.Figure 4, it is seen that remoter from the electrical distance of MT, node injects work(
The fluctuation that MT caused by the variation of rate contributes is more apparent.EB2 electrical distances as where EB18, with MT in Fig. 3 farther out, therefore are worked as
After load in EB18 nodes increases the workload demand of 1 unit, MT will accordingly increase the active power output of 1.1404 units.And Fig. 3
Middle EB19 is closer from MT electrical distances, be all increase by 1 unit workload demand after, MT need to only increase by 1.0007 units it is active go out
Power.The position of power distribution network interior joint injecting power variation is different, and influencing also difference, reason for the output of MT is power distribution network
Since electrical distance distance causes network loss not of uniform size during progress active power rebalancing.
S is calculated according to formula (16), (17)ggAs shown in figure 5, block diagram size indicates that gas net air pressure change is negative to combustion gas in figure
The sensitivity of lotus variation.It should be noted that the result of calculation of Fig. 5 is the result when MT accesses the GB11 of gas burning system.
When Gas Load 1 unit of increase in demand of MT, the air pressure of GB11 accordingly declines 0.3175 unit, and close
The GB2 air pressures of air source accordingly decline 0.1155 unit.The increase with flow in gas pipeline, pipe are can be seen that according to formula (1)
The pressure difference at road both ends will increase.When the combustion gas of 1 newly-increased unit is delivered to GB11 by combustion gas network, pipeline meeting in network
Cause the increase of draught head due to the increase of flow, the increase of draught head is cumulative together, is presented as that the air pressure of GB11 declines most
Decline unobvious for the apparent GB2 air pressures for being in network upstream.
In conjunction with See、SggS is further calculated according to formula (14)geAs shown in fig. 6, in figure determine a bit (x, y,
Z), expression be gas net node y level of sensitivity z of the air pressure to grid nodes x power swings.Gas net section as can be seen from FIG. 6
8,9,11 corresponding sensitivity numbers of point are larger, and air pressure change speed caused by grid nodes power swing is fast compared with other nodes.
Take the weak node of before sensitivity number 30% node as IEGS, at this time the weak node of combustion gas network be GB8, GB9,
GB11}。
To verify the validity of this method, the variation of air pressure under different injecting power growth levels is emulated by Load flow calculation
As shown in Figure 7.Each active growth rate of load bus (comparing origin node load level) is 2%, when the load growth rate of node is
When 10%, out-of-limit (circled in figure) occurs for GB11 air pressures.The increase meeting because of load is can be seen that from the data in Fig. 6
Make MT contribute to rise, corresponding gas net air pressure declines, and the air pressure of weak node reaches the boundary of safe operation at first.Another party
Face, GB11 is in line end and heavy duty from actual physics system, Fig. 3, and Gas Load variable quantity is to system air pressure
It influences maximum.It can be seen that sensitivity analysis result meets the actual conditions of network.
2) gas storage position is preferred
It is to improve combustion gas air pressure level and ensure having for system safety operation that caisson is arranged at system weak node
Effect means.Within the system short time use gas peak period, caisson by system inject combustion gas, play stable air pressure level,
Ensure the effect of air-flow continuous and stable.But it is uneconomic to all nodes setting caisson, the embodiment of the present invention is according to spirit
Gas storage is arranged in sensitivity result of calculation at system weak node.
By taking load growth 10% as an example, one capacity of setting is 20m3Caisson.Following 3 kinds are selected in gas storage position
Scheme：
1) it is not provided with gas storage
2) gas storage is set at non-weak link such as GB4；
3) gas storage is set at system weak node GB11.
In Load flow calculation, gas storage equipment approximate can be handled using the corresponding Gas Load of reduction.Pass through trend meter
The air pressure level for calculating system under different on-positions is as shown in Figure 8.Although effectively having been raised in GB4 nodes access gas storage equipment whole
Body air pressure level, but scheme 2 and scheme 3 are compared, gas storage, which is arranged, in weak node is lifted more effectively system air pressure, is
System safe operation nargin significantly improves.Static sensitivity analysis method of the embodiment of the present invention based on Unified Power Flow model, with electricity
Based on Force system and gas burning system Unified Power Flow model, IEGS gaseous-pressures-node injecting power sensitivity matrix is had studied,
Analyze influence of the grid nodes injecting power to gaseous-pressure, the weak link of localization region integrated energy system, for comprehensive energy
Source system safe and stable operation provides auxiliary information, the safety of effective lifting system.It can be obtained by sample calculation analysis and such as be drawn a conclusion：
1) gaseous-pressure of node-node injecting power level of sensitivity can effectively reflect each load bus air pressure of fuel gas network
Safe operation nargin；
2) sensitivity index can be used for quick positioning system weak link, accurately determines air pressure in system and most easily crosses the border section
Point avoids a large amount of calculating of point by point method；3) gaseous-pressure-node injecting power sensitivity is true, effectively reflects load section
The fluctuation of point Injection Space installs gas storage equipment to the influence degree of fuel gas network air pressure level in system weakness, can be effective
Improve air pressure safe operation nargin, the safety of lifting system.
Bibliography
[1] National Energy Board energy developments " 13 " plan [C] 2016:17-18.
[2] Jia Hongjie, wangdan, Xu Xiandong wait research [J] the electric system of regional complex energy resource system some problems automatic
Change, 2015,39 (7):198-207.
JIA H,Wang D,XU X,et al.Research on some key problems related to
integrated energy systems[J].Automation of Electric Power Systems,2015,39(7):
198-207.(in Chinese).
[3] pot that the Taiwan Wang Ping has a power failure on a large scale should who carry on the back [N] People's Daily Overseas Editions, 2017-08-16 (004)
[4]Fairley P.The troubled link between gas and electricity grids
[News][J].IEEE Spectrum,2016,53(6):11-12.
[5]Chertkov M,Backhaus S,Lebedev V.Cascading of fluctuations in
interdependent energy infrastructures:Gas-grid coupling[J].Applied Energy,
2015,160:541-551.
[6]Touretzky C R,McGuffin D L,ZIEGSmer J C,et al.The effect of
distributed electricity generation using natural gas on the electric and
natural gas grids[J].Applied Energy,2016,177:500-514.
[7]Martinez-Mares A,Fuerte-Esquivel C R.A unified gas and power flow
analysis in natural gas and electricity coupled networks[J].IEEE Transactions
on Power Systems,2012,27(4):2156-2166.
[8]Zlotnik A,Roald L,Backhaus S,et al.Coordinated scheduling for
interdependent electric power and natural gas infrastructures[J].IEEE
Transactions on Power Systems,2017,32(1):600-610.
[9] Wang Yingrui, Zeng Bo, Guo Jing wait electric-thermals-gas integrated energy system multipotency flow calculation methodologies [J] power grid skills
Art, 2016,40 (10):2942-2950.
WANG Y,ZENG B,GUO J,et al.Multi-energy flow cal-culation method for
integrated energy system containing electricity,heat and gas[J].Power System
Technology,2016,40(10):2942-2951.
[10] Wang Wei is bright, wangdan, Jia Hongjie, and is waited to consider electric power-natural gas region comprehensive energy of natural gas network state
Systematic steady state analyzes [J] Proceedings of the CSEEs, 2017,37 (5):1293-1304.
Wang Weiliang,Wang Dan,Jia Hongjie,et al.Steady state analysis of
electricity-gas regional integrated energy system with consideration of NGS
network status[J]Proceedings of the CSEE,,2017,37(5):1293-1304.
[11] Xu Xiandong, Jia Hongjie, Jin little Long wait regional complexs energy resource system electricity/gas/hot mixing power flow algorithm research
[J] Proceedings of the CSEEs, 2015,35 (14):3634-3642.
Xiandong X,Hongjie J,Xiaolong J,et al.Study on hybrid heat-gas-power
flow algorithm for integrated community energy system[J].Proceedings of the
CSEE,2015,35(14):3634-3642.
[12] a kind of comprehensive multi-index static voltage stability weakness busbar identifications of Luo Qingyue, Luo Ting, Sun Liuqing, Liu Baiyang
Method [J] electric power system protection and controls, 2017, V45 (3):82-86.
LUO Qingyue,LUO Ting,SUN Liuqing,et al.A multi-index integrated
static voltage stability assessment method for weak buses identification[J]
.Power System Protection and Control,2017,V45(3):82-86.
[13]Abeysekera M,Wu J,Jenkins N,et al.Steady state analysis of gas
networks with distributed injection of alternative gas[J].Applied Energy,
2016,164:991-1002.
[14] Fang Hengfu, Sheng Wanxing, Wang Jinli wait the power distribution stations three-phase load unbalance real-time online administering method to grind
Study carefully [J] Proceedings of the CSEEs, 2015,35 (9):2185-2193.
[15]Xu X,Jia H,Chiang H D,et al.Dynamic modeling and interaction of
hybrid natural gas and electricity supply system in microgrid[J].IEEE
Transactions on Power Systems,2015,30(3):1212-1221.
[16] Kong Xiangyu, room is big-and-middle, and marquis helps Incremental Transmission Loss algorithm [J] the electric power network techniques of magnificent based on DC power flow,
2007,31(15):39-43.
KONG X,FANG D,HOU Y.An incremental transmission loss algorithm based
on DC power flow[J].Power System Technology,2007,15:012.
To the model of each device in addition to doing specified otherwise, the model of other devices is not limited the embodiment of the present invention,
As long as the device of above-mentioned function can be completed.
It will be appreciated by those skilled in the art that attached drawing is the schematic diagram of a preferred embodiment, the embodiments of the present invention
Serial number is for illustration only, can not represent the quality of embodiment.
The foregoing is merely presently preferred embodiments of the present invention, is not intended to limit the invention, it is all the present invention spirit and
Within principle, any modification, equivalent replacement, improvement and so on should all be included in the protection scope of the present invention.
Claims (6)
1. a kind of static sensitivity analysis method coupling integrated energy system towards electric-gas, which is characterized in that the static spirit
Basis of sensitivity analysis method includes the following steps：
Establish by gas burning system model, distribution system model and the energy coupling link model group at IEGS Unified Power Flow moulds
Type；
According to IEGS Unified Power Flow models, the Jacobian matrix of current stable operating point and current point of operation is obtained；
Gaseous-pressure-node injecting power sensitivity is calculated according to Jacobian matrix；To gaseous-pressure-node injecting power spirit
Sensitivity carries out size sequence, and the big node set of sensitivity is the weak link of IEGS.
2. a kind of static sensitivity analysis method coupling integrated energy system towards electric-gas according to claim 1,
It is characterized in that, the IEGS Unified Power Flows model is specially：
In formula, x=[θ, V, p] indicates that the state variable of IEGS, respectively power distribution network each node phase angle in addition to balance nodes, electricity are negative
Lotus node voltage, Gas Load node air pressure；U=[PPDS, QPDS, LNGS] indicate IEGS disturbance variable, respectively node inject
Active power, reactive power and Gas Load；A indicates independent parameter；PPDSActive power is injected for node；For grid nodes electricity
Pressure；Y is node admittance matrix, and * is complex conjugate；ANGSFor the node-branch incident matrix of combustion gas network；For light-duty combustion gas
The electromotive power output of turbine；cgeFor transformation ratio.
3. a kind of static sensitivity analysis method coupling integrated energy system towards electric-gas according to claim 1,
It is characterized in that, it is described gaseous-pressure-node injecting power sensitivity is calculated according to Jacobian matrix to be specially：
According to Jacobian matrix, output-section of gaseous-pressure-Gas Load sensitivity and lightweight gas turbine is calculated separately
Point injecting power sensitivity；
Pass through the output of gaseous-pressure-Gas Load sensitivity and lightweight gas turbine-node injecting power Calculation of Sensitivity
Go out gaseous-pressure-node injecting power sensitivity.
4. a kind of static sensitivity analysis side coupling integrated energy system towards electric-gas according to claim 1 or 3
Method, which is characterized in that the gaseous-pressure-node injecting power sensitivity is specially：
Wherein, PPDSActive power is injected for node；P is combustion gas node air pressure；LNGSFor Gas Load；Represent the defeated of MT
Go out electrical power；For the gas consumption of MT.
5. a kind of static sensitivity analysis method coupling integrated energy system towards electric-gas according to claim 3,
It is characterized in that,
The gaseous-pressure-Gas Load sensitivity SggEqual to Jacobian matrix J in natural gas system Load Flow SolutionggThat inverts is negative
Number：
In formula, ANGSFor the node-branch incident matrix of combustion gas network；fiFor the flow of pipeline i；ΠiFor the gaseous-pressure of pipeline i
Difference；npFor the number of gas pipeline；T is matrix transposition；Diag () indicates construction diagonal matrix.
6. a kind of static state coupling integrated energy system towards electric-gas according to any claim in claim 1-5
Sensitivity Analysis Method, which is characterized in that
If a certain node i injecting powerIt changes, has for power distribution network：
In formula,It is right for distribution network lossLocal derviation is asked to have：
In formula, Section 2 is Incremental Transmission Loss on the right of equation, represents network loss with the variation relation between node injecting power.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810123783.0A CN108416507B (en) | 2018-02-07 | 2018-02-07 | Static sensitivity analysis method for electricity-gas coupling comprehensive energy system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810123783.0A CN108416507B (en) | 2018-02-07 | 2018-02-07 | Static sensitivity analysis method for electricity-gas coupling comprehensive energy system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108416507A true CN108416507A (en) | 2018-08-17 |
| CN108416507B CN108416507B (en) | 2021-08-20 |
Family
ID=63127971
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810123783.0A Active CN108416507B (en) | 2018-02-07 | 2018-02-07 | Static sensitivity analysis method for electricity-gas coupling comprehensive energy system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108416507B (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110543695A (en) * | 2019-08-14 | 2019-12-06 | 天津大学 | electric-gas coupling comprehensive energy system feasible region calculation method |
| CN110601185A (en) * | 2019-09-17 | 2019-12-20 | 武汉大学 | Unified power flow model and random matrix-based comprehensive energy system weak point identification method |
| CN110647040A (en) * | 2019-10-10 | 2020-01-03 | 东北电力大学 | Safety control method and device of comprehensive energy system |
| CN111767660A (en) * | 2020-09-02 | 2020-10-13 | 国网江西省电力有限公司电力科学研究院 | Optimal arrangement method and device for gas storage devices of electric coupling energy system |
| CN112510680A (en) * | 2020-10-26 | 2021-03-16 | 华能定边新能源发电有限公司 | Collaborative planning method for distributed power supply and electric-to-gas equipment in comprehensive energy system |
| CN112964961A (en) * | 2021-02-04 | 2021-06-15 | 南方电网科学研究院有限责任公司 | Electric-gas coupling comprehensive energy system fault positioning method and system |
| CN113609651A (en) * | 2021-07-13 | 2021-11-05 | 华南理工大学 | Unified power flow model-based electric-gas interconnection system action mechanism analysis method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105186525A (en) * | 2015-10-29 | 2015-12-23 | 山东大学 | Reactive voltage control partitioning method under wind power integration |
| US20150378387A1 (en) * | 2013-05-14 | 2015-12-31 | Rensselaer Polytechnic Institute | Methods of computing steady-state voltage stability margins of power systems |
| CN105356447A (en) * | 2015-10-27 | 2016-02-24 | 河海大学 | Analysis method for steady-state security region of electric-pneumatic interconnected integrated energy system |
| CN106777708A (en) * | 2016-12-21 | 2017-05-31 | 天津大学 | Steady state analysis method of electric power-natural gas regional comprehensive energy system |
-
2018
- 2018-02-07 CN CN201810123783.0A patent/CN108416507B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150378387A1 (en) * | 2013-05-14 | 2015-12-31 | Rensselaer Polytechnic Institute | Methods of computing steady-state voltage stability margins of power systems |
| CN105356447A (en) * | 2015-10-27 | 2016-02-24 | 河海大学 | Analysis method for steady-state security region of electric-pneumatic interconnected integrated energy system |
| CN105186525A (en) * | 2015-10-29 | 2015-12-23 | 山东大学 | Reactive voltage control partitioning method under wind power integration |
| CN106777708A (en) * | 2016-12-21 | 2017-05-31 | 天津大学 | Steady state analysis method of electric power-natural gas regional comprehensive energy system |
Non-Patent Citations (2)
| Title |
|---|
| YUBIN HE等: "Decentralized Optimization of Multi-Area Electricity-Natural Gas Flows Based on Cone Reformulation", 《IEEE》 * |
| 艾小猛等: "一种考虑天然气系统动态过程的气电联合系统优化运行模型", 《电网技术》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110543695A (en) * | 2019-08-14 | 2019-12-06 | 天津大学 | electric-gas coupling comprehensive energy system feasible region calculation method |
| CN110543695B (en) * | 2019-08-14 | 2023-08-01 | 天津大学 | Electric-gas coupling comprehensive energy system feasible region calculation method |
| CN110601185A (en) * | 2019-09-17 | 2019-12-20 | 武汉大学 | Unified power flow model and random matrix-based comprehensive energy system weak point identification method |
| CN110647040A (en) * | 2019-10-10 | 2020-01-03 | 东北电力大学 | Safety control method and device of comprehensive energy system |
| CN110647040B (en) * | 2019-10-10 | 2022-10-04 | 东北电力大学 | Safety control method and device of comprehensive energy system |
| CN111767660A (en) * | 2020-09-02 | 2020-10-13 | 国网江西省电力有限公司电力科学研究院 | Optimal arrangement method and device for gas storage devices of electric coupling energy system |
| CN112510680A (en) * | 2020-10-26 | 2021-03-16 | 华能定边新能源发电有限公司 | Collaborative planning method for distributed power supply and electric-to-gas equipment in comprehensive energy system |
| CN112964961A (en) * | 2021-02-04 | 2021-06-15 | 南方电网科学研究院有限责任公司 | Electric-gas coupling comprehensive energy system fault positioning method and system |
| CN113609651A (en) * | 2021-07-13 | 2021-11-05 | 华南理工大学 | Unified power flow model-based electric-gas interconnection system action mechanism analysis method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108416507B (en) | 2021-08-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108416507A (en) | The static sensitivity analysis method of integrated energy system is coupled towards electric-gas | |
| Sedghisigarchi et al. | Dynamic and transient analysis of power distribution systems with fuel Cells-part II: control and stability enhancement | |
| Ara et al. | Modelling of Optimal Unified Power Flow Controller (OUPFC) for optimal steady-state performance of power systems | |
| Oo et al. | Transmission capacity improvement using unified power flow controller with new control strategy | |
| Majid et al. | Minimization of transmission loss in application of HVDC networks under load increase scenario | |
| Hou et al. | Optimal successive start-up strategy of two hydraulic coupling pumped storage units based on multi-objective control | |
| Aftab et al. | IEC 61850 communication based dual stage load frequency controller for isolated hybrid microgrid | |
| Antoniadou-Plytaria et al. | Chalmers campus as a testbed for intelligent grids and local energy systems | |
| CN106816879B (en) | A kind of active distribution network voltage control method for coordinating based on urgent need response | |
| Kobayashi et al. | Reactive power control method between DG using ICT for proper voltage control of utility distribution system | |
| Zhang et al. | Multi-objectives OPF of AC-DC systems considering VSC-HVDC integration | |
| Jain et al. | Studying the impact of solar PV on power system dynamics using integrated transmission and distribution network models | |
| Almassalkhi | Optimization and Model-predictive Control for Overload Mitigation in Resilient Power Systems. | |
| Li et al. | Large‐scale source‐grid‐load friendly interactive system introduction and real load shedding verification test technology | |
| Liu et al. | Feasible Operation Regions of Electricity-gas Integrated Energy Distribution System | |
| Hashemi et al. | Impact of distributed generation on unbalanced distribution networks | |
| Lin et al. | Empowering the New Intelligent Power System with Digital Technologies to Achieve the “Dual Carbon” Goals | |
| Shi et al. | Quasi-AC Optimal Power Flow for VSC-MTDC Systems | |
| Metweely et al. | Optimization of optimal power flow problems with FACTS devices using PSO technique | |
| Xu et al. | Integrated optimal power flow for distribution networks in local and urban scales | |
| Jiankun et al. | Method of calculating the maximum transmission capability with unified power flow controller | |
| Feng et al. | Research on Path Connectivity Optimization Model for Active Distribution Network Reconfiguration | |
| Haoming et al. | Engineering practice of large‐scale millisecond accurate load control system of large power grid | |
| Wang et al. | A review of the global energy internet and the suggestions to China | |
| Li et al. | A SOP planning method for distribution network based on improved genetic algorithm considering EV integration |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |