CN109918795A - A kind of energy hinge canonical matrix Steady state modeling method for taking into account physical property and reliability - Google Patents
A kind of energy hinge canonical matrix Steady state modeling method for taking into account physical property and reliability Download PDFInfo
- Publication number
- CN109918795A CN109918795A CN201910178690.2A CN201910178690A CN109918795A CN 109918795 A CN109918795 A CN 109918795A CN 201910178690 A CN201910178690 A CN 201910178690A CN 109918795 A CN109918795 A CN 109918795A
- Authority
- CN
- China
- Prior art keywords
- energy
- hinge
- coupling
- energy hinge
- reliability
- 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
The present invention discloses a kind of energy hinge canonical matrix Steady state modeling method for taking into account physical property and reliability, it is related to integrated energy system safe and stable operation field, the following steps are included: it is uncertain to describe scene first with stochastic response surface, consider integrated energy system energy hinge coupling matrix and its failure state transfer matrix, establishes energy hinge canonical matrix steady-state model；Then electricity, the air and heat impedance operator for considering the coupling of energy hinge, establish the energy hinge steady state energy function for describing its characteristic, and the Bifurcation Characteristics of analysing energy hinge nonlinear balance equation group solve the maximum load capability of energy hinge stable state；In conjunction with the multipotency stream transfer factor of electric power system tide transfer factor solving system, the coupling mechanism of transmission and reliability that integrated energy system describes overload type failure are described using the multipotency stream transfer factor acquired.The present invention is adapted to energy source interconnection development trend, most important to the safe and stable operation safety analysis of integrated energy system.
Description
Technical field
The present invention relates to a kind of energy hinge canonical matrix Steady state modeling methods for taking into account physical property and reliability, belong to
Integrated energy system safe and stable operation field.
Background technique
Currently, the transition of China's energy be faced with demand slows down, traditional excess capacity, environmental problem are prominent, whole efficiency compared with
The problems such as low.In this context, integrated energy system comes into being, and theory breaches the technology of traditional energy system, market
It is the comprehensive energy resource system of all kinds of energy unified plannings such as electricity, air and heat, cold, United Dispatching with management barrier, for pushing
Energy resource structure transition, promotes China's energy revolution to be of great significance.Integrated energy system is bound to that future source of energy system will be become
Main Morphology.
Currently, the concept of integrated energy system, physical structure and correlation model are had been carried out both at home and abroad more deep
Research, but existing research mostly be for it is a certain it is specific/assume Regional Energy system modeled, to comprehensive energy system
All kinds of typical physical equipment and its mathematical model do not carry out systematic combing and summary in system.
The present invention is in the case where considering honourable various uncertain factors, by the various analysis methods of conventional electric power system
Integrated energy system energy hinge coupling matrix Steady state modeling is promoted the use of, to analysis integrated energy resource system energy hinge coupling
The transient stability of conjunction is of great significance.
Summary of the invention
In view of the deficiencies of the prior art, " a kind of energy hinge canonical matrix for taking into account physical property and reliability of the invention
Steady state modeling method " establishes the energy hinge stable state canonical matrix model for considering the coupling of multipotency stream and its failure；Consider energy
The electricity of hinge coupling, air and heat impedance operator, establish the energy hinge steady state energy function for describing its characteristic, analysing energy hinge
The Bifurcation Characteristics of nonlinear balance equation group solve the maximum load capability of energy hinge stable state；Turn in conjunction with electric power system tide
The multipotency stream transfer factor for moving factor solving system describes integrated energy system using the multipotency stream transfer factor acquired and describes
The coupling mechanism of transmission and reliability of overload type failure.
The present invention adopts the following technical scheme: a kind of energy hinge canonical matrix stable state for taking into account physical property and reliability
Modeling method, this method comprises the following steps:
Step 1: establishing the energy hinge stable state canonical matrix model for considering comprehensive energy coupling and its failure；
Step 2: considering electricity, the air and heat impedance operator of the coupling of energy hinge, the energy hinge that foundation describes its characteristic is steady
State energy function；
Step 3: considering electricity, the air and heat impedance operator of the coupling of energy hinge, analysing energy hinge nonlinear balance equation
The Bifurcation Characteristics of group；
Step 4: solving the maximum load capability of energy hinge stable state；
Step 5: in conjunction with the multipotency stream transfer factor of electric power system tide transfer factor solving system；
Step 6: describing integrated energy system using the multipotency stream transfer factor that step (5) are acquired and describe overload type failure
The coupling mechanism of transmission and reliability.
Detailed description of the invention
Fig. 1 is inventive algorithm flow chart.
Specific embodiment
The present invention the following steps are included:
Step 1: establishing the energy hinge stable state canonical matrix model for considering comprehensive energy coupling and its failure；
(1) scene uncertainty is described with stochastic response surface
Consider that the probabilistic probabilistic loadflow model of scene is represented by F, state variable to be asked is the defeated of model F
Y=(Y is responded out1,Y2…,Ym)T, each scene uncertain factor is the input stochastic variable X=(X of model F1,X2…Xn)T,
Model F may be expressed as:
Y=F (X)=F (X1,X2…Xn) (1)
Scene uncertainty, the i.e. probability distribution of the response Y of analysis model formula (1) are described with response phase method immediately.
(2) consider energy hinge coupling matrix and its failure state transfer matrix
Various elements are mutually coupled in integrated energy system first, and connection relationship is complicated, and problem to be solved is exactly
How equivalent inner element couples transforming relationship, establishes the external port model of mapping electricity, air and heat input and output.In order to solve this
A problem describes multiple-energy-source input/output relation with a coupling matrix using energy hinge model.That is:
P in formulajAnd LiThe different form energy for respectively inputting, exporting, i ∈ { 1,2 ..., n }, j ∈ { 1,2 ... m }；cij
For coupling factor.
Then the failure state transfer matrix for considering integrated energy system energy hinge, finally can be obtained energy hinge stable state
Canonical matrix model.Its state-transition matrix is as follows:
In formulaFor the probability of happening of primary fault；For derivative
The probability that failure occurs.
Step 2: considering electricity, the air and heat impedance operator of the coupling of energy hinge, the energy hinge that foundation describes its characteristic is steady
State energy function；
Energy function method is commonly used in analysis transient stability analysis of power system, i.e., by including net in construction electric system
The energy function of the different elements such as network, generator, load carrys out the transient stability of analysis system.By the impedance in electric system
Characteristic expansion is applied to integrated energy system, considers electricity, the air and heat impedance operator of the coupling of energy hinge, establishes and describe its characteristic
Energy hinge steady state energy Functional Analysis integrated energy system transient stability.
Step 3: considering electricity, the air and heat impedance operator of the coupling of energy hinge, analysing energy hinge nonlinear balance equation
The Bifurcation Characteristics of group；
For so-called saddle i.e. when bifurcation occurs in nonlinear system, the equalization point of system becomes one in fork
The phenomenon that a saddle point and a node.Electricity, the air and heat impedance operator of the coupling of energy hinge are considered first, then analysing energy pivot
The Bifurcation Characteristics of knob nonlinear balance equation group carry out structural machine to integrated energy system energy hinge using bifurcation theory
Reason analysis.
Step 4: obtaining energy hinge stable state maximum load capability；
Electricity, air and heat impedance operator and the comprehensive energy of the coupling of energy hinge are considered on the basis of step (2) and step (3)
Source system capacity hinge saddle phenomenon obtains energy hinge stable state maximum load capability.
Step 5: in conjunction with the multipotency stream transfer factor of electric power system tide transfer factor solving system；
The transfer of conventional electric power system load flow usually occurs after the topological structure of network changes, when certain in system is propped up
After road is removed, trend thereon transfers to the branch that other in system operate normally and may cause certain normal branches
Overload.
When certain branch excision event occurs, the transfer factor column vector relative to other branches in network is as follows:
Flow transferring relativity factor is the only proportionality coefficient with network topology structure and relating to parameters, can be cut branch occurs
Except event is calculated before.Therefore the flow transferring relativity factor in electric system can be expanded and be applied to comprehensive energy by we
In system, transfer factor can be flowed by obtaining integrated energy system multipotency stream.
Throughput in gas distributing system, the water flow in heat distribution pipe network are regarded as comprehensive energy by multipotency stream transfer factor
Heterogeneous energy stream in the system of source, the original energy of removed branch or pipeline after system to that branch or pipeline excision event occur flow
The behavior being transferred on other branches or pipeline is referred to as multipotency circulation and moves, and above-mentioned energy stream is to be transferred to it according to a certain percentage
On his branch or pipeline, this ratio is defined as multipotency stream transfer factor.
Step 6: the coupling of integrated energy system overload type failure is described using the multipotency stream transfer factor that step (5) are acquired
Close the mechanism of transmission and reliability.
After branch or pipeline excision event occurs according to system, it is removed branch or original can circulate of pipeline moves to other branch
Behavior on road or pipeline can analyze the coupling mechanism of transmission of integrated energy system overload type failure；Then according to multipotency stream
The certain proportion being transferred on other branches or pipeline describes the reliability of integrated energy system.
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 institute
There is equivalent technical solution to also belong to protection category of the invention.
Claims (5)
1. a kind of energy hinge canonical matrix Steady state modeling method for taking into account physical property and reliability, it is characterised in that:
Step (1) establishes the energy hinge stable state canonical matrix model for considering comprehensive energy coupling and its failure；
Step (2) considers the electricity of energy hinge coupling, air and heat impedance operator, establishes the energy hinge stable state energy for describing its characteristic
Flow function；
Step (3) considers the electricity of energy hinge coupling, air and heat impedance operator, analysing energy hinge nonlinear balance equation group
Bifurcation Characteristics；
The maximum load capability of step (4) solution energy hinge stable state；
Step (5) combines the multipotency stream transfer factor of electric power system tide transfer factor solving system；
Step (6) is passed using the coupling that the multipotency stream transfer factor that step (5) are acquired describes integrated energy system overload type failure
Broadcast mechanism and reliability.
2. a kind of energy hinge canonical matrix Steady state modeling side for taking into account physical property and reliability according to claim 1
Method, which is characterized in that the energy hinge stable state canonical matrix for considering comprehensive energy coupling and its failure is established in the step 1)
Model.Scene uncertainty is described with stochastic response surface first, then considers integrated energy system energy hinge coupling matrix
And its failure state transfer matrix, establish energy hinge canonical matrix steady-state model.
3. a kind of energy hinge canonical matrix Steady state modeling side for taking into account physical property and reliability according to claim 1
Method, which is characterized in that the electricity of analysing energy hinge coupling, air and heat impedance operator in the step 2), by electric system energy
Function obtains comprehensive energy energy function.The impedance operator in electric system is promoted the use of into integrated energy system analysis first
The electricity of energy hinge coupling, air and heat impedance operator；Finally integrated energy system energy pivot is constructed by electric system energy function
The energy function of knob coupling carrys out the transient stability of analysis integrated energy resource system.
4. a kind of energy hinge canonical matrix Steady state modeling side for taking into account physical property and reliability according to claim 1
Method, which is characterized in that energy hinge stable state maximum load capability is obtained in the step 4), on the basis of step 2) and step 3)
Consider that electricity, air and heat impedance operator and the integrated energy system energy hinge saddle phenomenon of the coupling of energy hinge obtain energy
Hinge stable state maximum load capability.
5. a kind of energy hinge canonical matrix Steady state modeling side for taking into account physical property and reliability according to claim 1
Method, which is characterized in that the multipotency stream transfer factor that the step 5) acquires, is the concept based on electric power system tide transfer factor
It acquires, i.e., the heterogeneous energy being regarded as the throughput in gas distributing system, the water flow in heat distribution pipe network in integrated energy system
Stream, original can circulate of removed branch or pipeline after system to that branch or pipeline excision event occur move to other branches or pipe
Behavior on road is referred to as multipotency circulation and moves, and above-mentioned energy stream is to be transferred on other branches or pipeline according to a certain percentage, by this
One ratio is defined as multipotency stream transfer factor.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910178690.2A CN109918795A (en) | 2019-03-11 | 2019-03-11 | A kind of energy hinge canonical matrix Steady state modeling method for taking into account physical property and reliability |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910178690.2A CN109918795A (en) | 2019-03-11 | 2019-03-11 | A kind of energy hinge canonical matrix Steady state modeling method for taking into account physical property and reliability |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109918795A true CN109918795A (en) | 2019-06-21 |
Family
ID=66964013
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910178690.2A Pending CN109918795A (en) | 2019-03-11 | 2019-03-11 | A kind of energy hinge canonical matrix Steady state modeling method for taking into account physical property and reliability |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109918795A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110266004A (en) * | 2019-06-28 | 2019-09-20 | 西安交通大学 | A kind of standardization construction method of integrated energy system energy hinge model |
| CN113327044A (en) * | 2021-06-10 | 2021-08-31 | 国家电网有限公司 | Collaborative operation quality analysis system of energy hub |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20040225486A1 (en) * | 2003-05-08 | 2004-11-11 | Mullis Vance Condary | System for generating a resultant model for a power system |
| 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 |
| CN106058942A (en) * | 2016-08-03 | 2016-10-26 | 长沙理工大学 | Energy hub optimizing model taking wind power nondeterminacy into consideration and including power to gas and CCHP |
| CN106777708A (en) * | 2016-12-21 | 2017-05-31 | 天津大学 | Steady state analysis method of electric power-natural gas regional comprehensive energy system |
| CN107992963A (en) * | 2017-11-23 | 2018-05-04 | 东南大学 | Coordinate the optimization method of bilateral interests based on multipotency streaming system Bi-level Programming Models |
| CN108594658A (en) * | 2018-04-24 | 2018-09-28 | 长沙理工大学 | Multi-objective optimization model for probability maximum load margin of electric-gas coupling system and solving method thereof |
-
2019
- 2019-03-11 CN CN201910178690.2A patent/CN109918795A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20040225486A1 (en) * | 2003-05-08 | 2004-11-11 | Mullis Vance Condary | System for generating a resultant model for a power system |
| 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 |
| CN106058942A (en) * | 2016-08-03 | 2016-10-26 | 长沙理工大学 | Energy hub optimizing model taking wind power nondeterminacy into consideration and including power to gas and CCHP |
| CN106777708A (en) * | 2016-12-21 | 2017-05-31 | 天津大学 | Steady state analysis method of electric power-natural gas regional comprehensive energy system |
| CN107992963A (en) * | 2017-11-23 | 2018-05-04 | 东南大学 | Coordinate the optimization method of bilateral interests based on multipotency streaming system Bi-level Programming Models |
| CN108594658A (en) * | 2018-04-24 | 2018-09-28 | 长沙理工大学 | Multi-objective optimization model for probability maximum load margin of electric-gas coupling system and solving method thereof |
Non-Patent Citations (2)
| Title |
|---|
| YI WANG: "Standardized Matrix Modeling of Multiple Energy Systems", 《IEEE》 * |
| 鲍玉川等: "基于综合能源系统的最优潮流解耦算法", 《电网与清洁能源》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110266004A (en) * | 2019-06-28 | 2019-09-20 | 西安交通大学 | A kind of standardization construction method of integrated energy system energy hinge model |
| CN113327044A (en) * | 2021-06-10 | 2021-08-31 | 国家电网有限公司 | Collaborative operation quality analysis system of energy hub |
| CN113327044B (en) * | 2021-06-10 | 2023-11-21 | 国家电网有限公司 | Collaborative operation quality analysis system of energy hub |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106067677B (en) | A kind of electric-thermal couples multipotency streaming system static security analysis method | |
| CN108921727A (en) | Consider the regional complex energy resource system reliability estimation method of thermic load dynamic characteristic | |
| CN105576657B (en) | Random fuzzy continuous power flow modeling method for active power distribution network and power transmission network integrated system | |
| CN110134833A (en) | Diagram data modeling and method towards power network topology management | |
| CN112398946A (en) | Source network load storage scheduling control system and configuration method | |
| CN107330056A (en) | Wind power plant SCADA system and its operation method based on big data cloud computing platform | |
| CN101141064A (en) | Method for distributed tidal current analyzing by exchange boundary node state and net damage information | |
| CN109145428A (en) | A kind of information physical emerging system vulnerability assessment method under cascading failure mode | |
| CN105376156A (en) | Multi-attribute decision-making based power backbone transmission network route planning method | |
| CN104810853B (en) | Digraph-based islanding method of distribution network including distributed generations | |
| CN109409730A (en) | A kind of energy microgrid site selecting method based on complex network characteristic evaluation | |
| CN109918795A (en) | A kind of energy hinge canonical matrix Steady state modeling method for taking into account physical property and reliability | |
| CN107292481A (en) | A kind of power network key node appraisal procedure based on pitch point importance | |
| CN112491034B (en) | Electric power information physical system modeling method based on alternating current power flow model | |
| CN114529123A (en) | Urban intelligent energy network hierarchical planning method | |
| CN105335900A (en) | Integrated data processing method based on power equipment state monitoring system | |
| CN105790256B (en) | Power distribution network access distributed generation unit critical path recognition methods based on multi-agent technology | |
| WO2021203480A1 (en) | Power distribution network automation system comprehensive planning method based on reliability constraint | |
| CN109768543B (en) | Elastic bottom-preserving net rack search modeling method based on mixed integer linear programming | |
| CN109617132B (en) | Resource allocation and network reconstruction optimization method for improving resilience of elastic power distribution network | |
| CN108599974B (en) | Micro-grid distributed communication topology design method based on graph theory connectivity | |
| Changchao et al. | Research on the frequency synchronization control strategy for power system | |
| CN112287559A (en) | Power distribution network power supply capacity analysis method and device for cold and heat pipe network virtual energy storage | |
| CN109934486A (en) | A method of extracting the differentiation of multipotency stream load nargin and its interaction feature in faults coupling communication process | |
| CN115906605A (en) | Distributed information energy system multi-agent model construction method and application |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190621 |