CN109871647A - A kind of multipotency stream load nargin method for solving based in faults coupling communication process - Google Patents
A kind of multipotency stream load nargin method for solving based in faults coupling communication process Download PDFInfo
- Publication number
- CN109871647A CN109871647A CN201910178725.2A CN201910178725A CN109871647A CN 109871647 A CN109871647 A CN 109871647A CN 201910178725 A CN201910178725 A CN 201910178725A CN 109871647 A CN109871647 A CN 109871647A
- Authority
- CN
- China
- Prior art keywords
- energy
- load
- multipotency
- stream
- multipotency stream
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
Invention describes a kind of multipotency stream load nargin method for solving based in faults coupling communication process, are related to integrated energy system safety analysis field, comprising the following steps: prepare the initial data of multipotency streaming system；It is analyzed by multipotency stream load nargin, obtains multipotency stream load security domain；Based on multipotency stream load security domain and its boundary characteristic, current point of operation is obtained to the distance of corresponding security domain boundaries, load margin under stable situation is obtained, obtains the load margin of energy hinge.The present invention is adapted to energy source interconnection development trend, proposes the static security method towards multipotency streaming system, and the security and stability of multipotency streaming system operation can be improved.
Description
Technical field
The present invention relates to a kind of multipotency stream load nargin method for solving based in faults coupling communication process, belong to synthesis
Energy resource system safety analysis field.
Background technique
Energy internet is the product of internet and energy depth integration, is the research hotspot of current energy field with before
The important morphological of edge and future source of energy system.It is ground under national science foundation of the US subsidy in the research center FREEDM in 2008
Study carefully energy internet, the same year Germany E-Energy carries out 6 demonstrative projects and practices energy internet.Jeremy Rifkin works in 2011
" the third time industrial revolution " is using energy internet as the important composition of the third time industrial revolution.Since 2012, the country is also opened
Relevant discussion and research are opened up.Energy revolution, power system reform, internet+etc. national strategies driving under, energy internet
A large amount of concerns are caused in academia and industrial circle.
In recent years, with integrated energy system interconnection scale constantly expand, high permeability renewable energy large-scale grid connection,
Heterogeneous energy stream intercouples, power electronic equipment is widely applied, and in enhancing integrated energy system conveying capacity, improves comprehensive energy
While system operational flexibility and economic effect property, the problems such as also bringing power to inject fluctuation, trend distribution randomness, aggravation
Integrated energy system operation risk, can flow stability to system and bring the hidden danger for being difficult to estimate.
Load margin can flow the performance indicator of maintenance level as measurement integrated energy system, reacted system and born load
And when fault disturbance, maintenance can flow stable ability.The uncertain factor of influence load margin computational accuracy and efficiency is increasingly
It is more.How under energy internet and renewable energy large-scale grid connection background, the fast of integrated energy system load margin is realized
Speed, accurate calculating, are of great significance to multiple-energy-source coupled system steady stability online evaluation.
So far, different researchers proposes a variety of load margin indexs from different angles, generally falls into
Two classes: margin index and state index.Margin index good linearity, explicit physical meaning, but it is related to seeking for critical point, it calculates
Measure larger, calculation method has direct method and two kinds of Continuation Method.For the calculation amount for reducing margin index, many documents are proposed
A variety of different practical approaches for calculating air extracts.State index includes all kinds of sensitivity indexs, characteristic value
With Singularity Value Index etc., margin index is mainly load margin.Many stability index, application on site can be flowed although proposing
Disadvantage undesirable in calculating speed and efficiency is still remained, and although has proposed many sides calculated for load margin
Method, but there are still calculate the disadvantages of complicated.
The load margin of integrated energy system how is quickly and accurately calculated, realizes quickly commenting for static energy stream stability
Estimating is still research emphasis.
Summary of the invention
For technological deficiency, the present invention is the load margin for quickly and accurately calculating integrated energy system, realizes comprehensive energy
The rapid evaluation of source power system steady-state stability is proposed based on the direct method for calculating integrated energy system load margin based on event
Multipotency stream load nargin method for solving in barrier coupling communication process, to promote the computation rate of solving system load margin.
The present invention adopts the following technical scheme: a kind of multipotency stream load nargin based in faults coupling communication process solves
Method, this method comprises the following steps:
Step 1: preparing the initial data of integrated energy system, including electric system, natural gas system and therrmodynamic system
Each node, branch initial data estimate that equipment or online energy Load flow calculation are analyzed by energy stream mode, obtain comprehensive energy system
The current operating conditions of system；
Step 2: calculating the associated static of analysis integrated energy resource system node, branch, coupling propagation path and energy hinge
Stability index can be flowed, these static energy stream stability indexs include: the relevant risk index of natural gas air pressure, flow；Therrmodynamic system
Temperature, the relevant risk index of flow；The relevant risk index of air pressure, temperature, the voltage of energy hinge etc.；
Step 3: analysis being calculated by step 2 and determines that the crucial branch, key node, crucial coupling of integrated energy system pass
Broadcast path and key energy hinge；
Step 4: forming extended power flow equations, the load margin in integrated energy system is calculated；
Step 5: being analyzed by multipotency stream load nargin, obtain multipotency stream load security domain；
Step 6: being based on multipotency stream load security domain and its boundary characteristic, obtain current point of operation to corresponding security domain boundaries
Distance, obtain load margin under stable situation, obtain the load margin of energy hinge.
Detailed description of the invention
It in order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, below will be to specific embodiment party
Attached drawing needed in formula is briefly described, it should be apparent that, the accompanying drawings in the following description is only of the invention one
A little embodiments for those of ordinary skill in the art without creative efforts, can also be according to these
Attached drawing obtains other attached drawings.In the accompanying drawings: attached drawing 1 is that the present invention is a kind of negative based on the multipotency stream in faults coupling communication process
The whole implementation flow chart of lotus nargin method for solving.
Specific embodiment
To make the object, technical solutions and advantages of the present invention clearer, embodiment of the present invention is made below further
Ground detailed description.
The present invention gives specific implementation steps:
Step 1: as the load in integrated energy system increases, the branch of system can flow stability index and change, first
First establish branch model in integrated energy system:
Step 2: being by the stability index that can flow that further analytical calculation obtains out branch k-m
Generally during system loading increase, necessarily there is the static energy stream stability index an of branch to level off to 0, therefore
It can use static energy stream stability index of the static voltage stability index as system of the branch, the corresponding branch is exactly system
Crucial branch, i.e., when system is close to when can flow stability boundaris, the stability index that can flow of crucial branch levels off to 0 at first, so being
The static energy stream stability index of system can indicate are as follows:
Step 3: with the increase of system loading, the branch of system can flow stability index and change, and be write formula (2) as increasing
Amount form is as follows:
A in formula1=2XVk；b1=-4Z2；a2=2RVk；b2=-4Z2
Branch power linearisation can be expressed as
In formula, λ is load level, herein referred as load coefficient；KP、KQFor the active and idle growth ratio of node.
Step 4: formula (5) substitution formula (4) is obtained branch k-n'sIndex reaches load coefficient when 0 are as follows:
Step 5: taking system loading nargin of the minimum value as prediction in all branches, i.e.,
The load margin of integrated energy system and its corresponding can be predicted by the sensitivity using integrated energy system
Static energy is lost steady branch, it can determines the crucial branch in integrated energy system.
Step 6: the pluriderivative based on load parameter carries out initial value and estimates, and obtains critical point initial value a0, zero eigenvalue pair
The right feature vector initial value b answered0, and multifunctional system load margin initial value λ in steps of 50；
Step 7: constructing non-linear power flow equation group and corresponding update equation group based on above-mentioned parameter, pass through
Update equation group is split into several groups of systems of linear equations with coefficient low-dimensional matrix by depression of order transformation；Line based on Riccati transformation
The tree-like multi-body system characteristic value method for solving of property, establishes multipotency stream input and output Riccati transfer matrix recurrence relation；From tree-like
Each input terminal of system starts the Riccati transfer matrix that each energy hinge is successively acquired along transmission path, and then realizes comprehensive energy
The rapid solving of the load margin of energy hinge in the system of source.
Embodiments above is merely to illustrate the present invention, and not limitation of the present invention, in relation to the common of technical field
Technical staff can also make a variety of changes and modification without departing from the spirit and scope of the present invention, therefore all
Equivalent technical solution also belongs to protection category of the invention.
Claims (3)
1. a kind of multipotency stream load nargin method for solving based in faults coupling communication process, it is characterised in that: this method packet
Include following steps:
Step 1: preparing the initial data of integrated energy system, each section including electric system, natural gas system and therrmodynamic system
Point, branch initial data estimate that equipment or online energy Load flow calculation are analyzed by energy stream mode, obtain integrated energy system
Current operating conditions；
Step 2: the static energy stream for calculating analysis integrated energy resource system node, branch, coupling propagation path and energy hinge is stablized
Index, these static energy stream stability indexs include: the relevant risk index of natural gas air pressure, flow；Therrmodynamic system temperature, flow
Relevant risk index；The relevant risk index of air pressure, temperature, the voltage of energy hinge etc.；
Step 3: analysis being calculated by step 2 and determines that the crucial branch, key node, crucial coupling of integrated energy system propagate road
Diameter and key energy hinge；
Step 4: forming extended power flow equations, calculated load nargin；
Step 5: being analyzed by multipotency stream load nargin, obtain multipotency stream load security domain；
Step 6: be based on multipotency stream load security domain and its boundary characteristic, obtain current point of operation to corresponding security domain boundaries away from
From obtaining load margin under stable situation, obtain the load margin of energy hinge.
2. a kind of multipotency stream load nargin method for solving based in faults coupling communication process according to claim 1,
It is characterized in that, in step 3 key link of integrated energy system determination, be the static energy stream stability index by step 2
Calculated result analysis determination.
3. requiring the formation extended power flow equations, calculated load nargin according to 1 step 4 of right, it is characterised in that based on negative
The pluriderivative of lotus parameter carries out initial value and estimates, and obtains critical point initial value a0, the corresponding right feature vector initial value b of zero eigenvalue0,
And multifunctional system load margin initial value λ in step 30；Constructed based on above-mentioned parameter non-linear power flow equation group and
Corresponding update equation group is converted by depression of order update equation group splitting into several groups of linear equations with coefficient low-dimensional matrix
Group；Based on the linear tree-like multi-body system characteristic value method for solving of Riccati transformation, multipotency stream input and output Riccati is established
Transfer matrix recurrence relation；The Riccati of each energy hinge is successively acquired along transmission path since each input terminal of tree system
Transfer matrix, and then realize the rapid solving of the load margin of energy hinge in integrated energy system.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910178725.2A CN109871647A (en) | 2019-03-11 | 2019-03-11 | A kind of multipotency stream load nargin method for solving based in faults coupling communication process |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910178725.2A CN109871647A (en) | 2019-03-11 | 2019-03-11 | A kind of multipotency stream load nargin method for solving based in faults coupling communication process |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109871647A true CN109871647A (en) | 2019-06-11 |
Family
ID=66920121
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910178725.2A Pending CN109871647A (en) | 2019-03-11 | 2019-03-11 | A kind of multipotency stream load nargin method for solving based in faults coupling communication process |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109871647A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111416340A (en) * | 2020-03-10 | 2020-07-14 | 天津大学 | Regional comprehensive energy system optimization control method based on practical security domain |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130066480A1 (en) * | 2011-09-07 | 2013-03-14 | Quanta Associates, L.P | Real-time monitoring of electric power system voltage stability margins |
| CN105576657A (en) * | 2016-02-06 | 2016-05-11 | 马瑞 | Random fuzzy continuation power flow modeling and algorithm of active distribution network and transmission network integrated system under energy source interconnection |
| CN107423867A (en) * | 2017-04-05 | 2017-12-01 | 长沙理工大学 | A kind of Static Voltage Security of AC and DC power system containing wind power plant domain method |
| CN108521125A (en) * | 2018-04-27 | 2018-09-11 | 长沙理工大学 | A kind of modeling of electric-gas coupled system probability continuous tide and its method for solving |
| CN108804386A (en) * | 2018-07-09 | 2018-11-13 | 东北电力大学 | A kind of parallelization computational methods of power system load nargin |
| CN109359389A (en) * | 2018-10-18 | 2019-02-19 | 东北大学 | City electric car charging decision method based on typical load dynamic game |
-
2019
- 2019-03-11 CN CN201910178725.2A patent/CN109871647A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130066480A1 (en) * | 2011-09-07 | 2013-03-14 | Quanta Associates, L.P | Real-time monitoring of electric power system voltage stability margins |
| CN105576657A (en) * | 2016-02-06 | 2016-05-11 | 马瑞 | Random fuzzy continuation power flow modeling and algorithm of active distribution network and transmission network integrated system under energy source interconnection |
| CN107423867A (en) * | 2017-04-05 | 2017-12-01 | 长沙理工大学 | A kind of Static Voltage Security of AC and DC power system containing wind power plant domain method |
| CN108521125A (en) * | 2018-04-27 | 2018-09-11 | 长沙理工大学 | A kind of modeling of electric-gas coupled system probability continuous tide and its method for solving |
| CN108804386A (en) * | 2018-07-09 | 2018-11-13 | 东北电力大学 | A kind of parallelization computational methods of power system load nargin |
| CN109359389A (en) * | 2018-10-18 | 2019-02-19 | 东北大学 | City electric car charging decision method based on typical load dynamic game |
Non-Patent Citations (3)
| Title |
|---|
| 傅旭等: "考虑N-1静态电压稳定约束的电力系统负荷裕度新算法", 《陕西电力》 * |
| 秦佳倩: "电-气耦合系统概率 Pareto 最大负荷裕度及灵敏度分析", 《电网技术》 * |
| 顾俊杰等: "基于Riccati传递矩阵法的线性树形多体系统特征值求解", 《南京理工大学学报》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111416340A (en) * | 2020-03-10 | 2020-07-14 | 天津大学 | Regional comprehensive energy system optimization control method based on practical security domain |
| CN111416340B (en) * | 2020-03-10 | 2022-07-19 | 天津大学 | Regional comprehensive energy system optimization control method based on practical security domain |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103476051B (en) | A kind of communication net node importance evaluation method | |
| CN103401238B (en) | A kind of power load modelling approach based on Measurement-based approach | |
| CN104778367B (en) | Wide area Thevenin's equivalence parameter on-line calculation method based on a single state section | |
| CN108448568B (en) | Power distribution network hybrid state estimation method based on multiple time period measurement data | |
| CN101958543B (en) | Method for estimating three-phase impedance-free nonlinear multi-source state of transformer substation | |
| Yuan et al. | Research on dynamic load modelling based on power quality monitoring system | |
| CN107038292A (en) | A kind of many output of wind electric field correlation modeling methods based on adaptive multivariable nonparametric probability | |
| CN106198218A (en) | A kind of method of the monitoring core level pipeline fatigue using strain transducer | |
| CN102495288A (en) | Method for identifying and estimating power transmission line parameters based on SCADA (Supervisory Control and Data Acquisition) measurement data | |
| CN103595040A (en) | Power distribution network comprising direct-current distributed power supply and analysis method of comprehensive load characteristics of power distribution network | |
| CN103532137A (en) | Method for estimating state of three-phase four-wire low-voltage distribution network | |
| CN106786608B (en) | A kind of uncertain harmonic flow calculation method suitable for distributed generation resource access | |
| CN104052053B (en) | Based on the power distribution network double line terminal section impedance analysis method of Linear Regression Model in One Unknown | |
| CN108054757A (en) | A kind of embedded idle and voltage N-1 Close loop security check methods | |
| CN105811404B (en) | A kind of distribution static voltage stability situation monitoring method of transmission ＆ distribution collaboration | |
| CN102508086A (en) | Parameter estimation method for electric power circuit in T-type connection | |
| CN104636829A (en) | Decoupling algorithm for increasing temperature optimal power flow (OPF) calculation efficiency of electric power system | |
| CN106372440B (en) | A kind of adaptive robust state estimation method of the power distribution network of parallel computation and device | |
| CN102340139A (en) | System and method for dynamically measuring and calculating transmission capacity of transmission line | |
| CN105787604B (en) | A kind of distribution situation method for quickly sensing of transmission & distribution collaboration | |
| CN109871647A (en) | A kind of multipotency stream load nargin method for solving based in faults coupling communication process | |
| CN101534010B (en) | Method for solving the local boundary of cut-set voltage stability regions on the basis of perturbation | |
| CN102723711B (en) | Cable power distribution network estimation method based on credibility | |
| CN107958120A (en) | A kind of system Thevenin's equivalence calculation method of parameters based on power series expansion | |
| CN102142682B (en) | Method for calculating sensitivity of branch breakage based on direct-current flow model |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190611 |
|
| WD01 | Invention patent application deemed withdrawn after publication |