CN109684763A - It is a kind of based on individual be this model integrated energy system modeling method - Google Patents
It is a kind of based on individual be this model integrated energy system modeling method Download PDFInfo
- Publication number
- CN109684763A CN109684763A CN201910000445.2A CN201910000445A CN109684763A CN 109684763 A CN109684763 A CN 109684763A CN 201910000445 A CN201910000445 A CN 201910000445A CN 109684763 A CN109684763 A CN 109684763A
- Authority
- CN
- China
- Prior art keywords
- individual
- model
- follows
- integrated energy
- modules
- 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
Classifications
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
It is the integrated energy system modeling method of this model that the invention discloses a kind of based on individual, and steps are as follows: determining the physical model of integrated energy system, each comprising modules of integrated energy system is taken out accordingly, to modules naming number；Classify to all modules, determines basis modeling unit, and establish a body Model of basic modeling unit；A body Model of system modules is inherited out from basic modeling unit model, and establishes the connection relationship between each module；System model is established, and each intermodule input/output relation is obtained according to system model；Modules, and then solving system are solved using forward-backward sweep method.Each module plug and play in the present invention is effectively reduced the modeling of complication system and modifies the complexity of process, and has been sufficiently reserved the model characteristics of each module in system, has significant engineering practical value and wide application prospect.
Description
Technical field
The present invention relates to complex system modeling technical fields, and in particular to a kind of based on the individual comprehensive energy for this model
System modeling method.
Background technique
Integrated energy system is more complicated compared to common system, more enlarged.The difference of it and conventional system, one
It is the complexity of the heterogeneity and network in structure.Integrated energy system usually combines various energy resources, across natural gas, electricity
The different energy sources net such as energy, thermal energy, cold energy.Every kind of energy networks have different physical rules and characteristic, therefore are suitable for difference
Method Modeling, and have different method for solving.In addition, it has combined natural gas network, power grid, cold and heat supply network, each
Network is all runed by different economic entities, and in addition to necessity interaction, they seek privacy and the operation of internal information
Independence.Second, the structure of integrated energy system complexity leads to the complexity of its behavioral characteristics.In power grid, electric energy is with light
Speed is propagated, and can reach stable state within several seconds；And in cold and hot net, waterpower dynamic process is propagated with the velocity of sound, it usually needs several
It can be only achieved waterpower quasi-steady state within ten seconds.In contrast, heating power dynamic process is propagated with water velocity, and speed is slower, is needed rather
Clock reaches heating power quasi-steady state to dozens of minutes.The dynamic characteristic of each part of integrated energy system is different, needs a kind of suitable more
The model of modal dynamic analysis.
Existing modeling method is that power grid, natural gas grid, the cold and hot net in comprehensive energy net are regarded as a network to be collected
Middle modeling, this model is extremely complex, in order to guarantee the solvability of network, it will usually ignore many details, carry out line to model
Property；Meanwhile the model of this network entirety is upper extremely difficult in modification.In terms of optimization, it is assumed that they do not consider itself benefit
Benefit, using whole system interests as target, the arrangement of obedience system entirety very.In terms of dynamic simulation analysis, using system
One simulation step length emulates whole system, the too long dynamic characteristic for being difficult to get fast dynamics submodule of simulation step length,
And simulation step length is too in short-term, generates bulk redundancy calculating becoming in module fastly, wastes computer computing resource.
It is important to note that in the process of construction of future new era energy resource system, different part in system, no
Same function can not be completed by same tissue.In the case where lacking unified model, between the system of different tissues operation
It can not merge and coordinate, therefore also just cannot achieve the optimizing decision of total system.Energy resource system needs a model to consider
Various energy resource systems, different main market players, different time scales various problems requirement, and at present can there is no a unification
Modeling method.
Summary of the invention
The purpose of the present invention is to solve drawbacks described above in the prior art, and providing one based on individual is this model
Integrated energy system modeling method, this method can effectively realize the unified Modeling of the complication system with heterogeneous subsystem, fill
While code insurance stays the details characteristic of subsystem, guarantee that subsystems interaction under unified standard form forms complication system,
Modeling, optimization and dynamic simulation analysis suitable for energy resource system of new generation.
The purpose of the present invention can be reached by adopting the following technical scheme that:
It is a kind of based on individual be this model integrated energy system modeling method, steps are as follows for the modeling method:
S1, the physical model for determining integrated energy system take out each comprising modules of integrated energy system, to each accordingly
Comprising modules naming number；
S2, classify to all modules, determine basis modeling unit, and establish the individual mould of basic modeling unit
Type；
S3, a body Model that system modules are inherited out from basic modeling unit model, and establish between each module
Connection relationship；
S4, system model is established, and each intermodule input/output relation is obtained according to system model；
S5, according to known parameters, modules are gradually solved using forward-backward sweep method, and then obtain the solution of whole system.
Further, a body Model of the basic modeling unit is defined as follows:
A body Model is indicated using penton:
∑=(I, K, X, F, O)
Wherein, I is the input set of individual, and K is the intrinsic Knowledge Set of individual, and X is the state set inside individual, and F is individual
Behavior collection, O be individual output collection；
(1) individual is expressed as i in the input of t momentt
it=(i, t) | i ∈ I, t ∈ T }
Wherein, T is the time series set of individual with environmental interaction, and individual input set indicates are as follows:
I={ i1,i2,…,it}；
(2) individual is k in the intrinsic representation of knowledge of t momentt
kt=(k, t) | k ∈ K, t ∈ Tk}
Wherein, TkIt is time series set corresponding to the Knowledge Set of individual, individual Knowledge Set indicates are as follows:
K={ k1,k2,…,kt}；
(3) individual is expressed as X in the state of t momentt
Xt=(X, t) | X ∈ X, t ∈ [0, et]}
Wherein, etIt is the accurate maximum time step-length for capturing system dynamic characteristic, individual state set indicates are as follows:
X={ X1,X2,…,Xt}；
(4) state X of the individual at the t+1 momentt+1More new formula is expressed as
Xt+1=F (Xt,kt,it)
Wherein, F (,) is the behavior equation of individual, describes the specific behavior of individual and evolutionary mechanism, individual is in t
The state X at+1 momentt+1, it is a mapping of the individual in t moment state, it is related with input set, Knowledge Set, due to behavior collection
In the presence of individual can also automatically update the state of t moment in the case where no outside stimulus；
(5) individual is expressed as o in the output of t momentt
ot=(o, t) | i ∈ O, t ∈ T }
Wherein, individual output set representations are as follows:
O={ o1,o2,…,ot}。
Further, the system model is defined as follows:
M=(V, T)
Wherein, V={ Σ | Σ ∈ V }, V are sub- groups of individuals, and T is the association for describing connection relationship between sub- individual
Matrix, the constitution element T of incidence matrixijIt is defined as follows:
Tij=+1, work as ΣiOutput be ΣjInput；
Tij=-1, works as ΣjOutput be ΣiInput；
Tij=0, work as ΣiAnd ΣjBetween without direct correlation.
Further, steps are as follows for forward-backward sweep method in the step S5:
S501, since the node farthest from power supply point, it is primary against the direction of power transmission using route voltage rating
Calculate the power loss and power distribution in partition line impedance, wherein power loss Δ S calculation formula are as follows:
Wherein P ", Q " is respectively the active power and reactive power of line end, VNFor the voltage rating of route, R, X divide
Not Wei route resistance value and reactance value, then the trend S ' calculation formula of route front end are as follows:
S '=S "+Δ S
Wherein, S " is the trend of line end；
S502, the power acquired using step S501 calculate partition line along the direction of power transmission since power supply point
The voltage landing on road finds out each node voltage respectively, wherein the amplitude Δ V of voltage landingmWith the phase angle Δ V of voltage landingaPoint
It does not calculate are as follows:
Wherein V ' is the voltage of route front end, and P ', Q ' are respectively the active power and reactive power of route front end, then line
The voltage of road end is expressed as:
The present invention has the following advantages and effects with respect to the prior art:
1) it is disclosed by the invention based on individual be this model integrated energy system modeling method be a kind of unified complexity
System modeling method, in complication system individual of different nature can internal individually modeling, be sufficiently reserved the details of subsystem
Characteristic, while the interaction of unified standard can be carried out.
2) disclosed by the invention based on the individual integrated energy system modeling method for this model, it can fully consider complexity
The requirement of different main market players's require information privacies and operation independence, meets current energy reform in integrated energy system
Background.
3) present invention is used based on individual as the integrated energy system modeling method of this model, is suitble to different in complication system
The subsystem of time scale carries out simulation analysis, has greatly saved computing resource and has accelerated complication system simulation process.
Detailed description of the invention
Fig. 1 is the flow chart disclosed by the invention based on the integrated energy system modeling method that individual is this model；
Fig. 2 is integrated energy system schematic diagram in the embodiment of the present invention；
Fig. 3 is to carry out classification schematic diagram to all modules in the embodiment of the present invention；
Fig. 4 is a body Model connection relationship of modules in the embodiment of the present invention.
Specific embodiment
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is
A part of the embodiment of the present invention, instead of all the embodiments.Based on the embodiments of the present invention, those of ordinary skill in the art
Every other embodiment obtained without making creative work, shall fall within the protection scope of the present invention.
Embodiment
As shown in Fig. 2, the present embodiment illustrates this hair as simulation object using a simple 4 node power system
Bright is the integrated energy system modeling method of this model based on individual, Fig. 1 be embodiment it is disclosed based on individual be this model
The flow chart of integrated energy system modeling method, steps are as follows:
Step S1, the electric system schematic diagram of the present embodiment is as shown in Figure 2.Determine the physical model of integrated energy system,
Each comprising modules of integrated energy system are taken out accordingly, to modules naming number.To each composition of electric system in Fig. 2
Module number consecutively, including route 1, route 2, load 1, load 2, confluence 3, transformer, generator.
Step S2, classify to all modules, determine basis modeling unit, and establish of basic modeling unit
Body Model.Classify to the electric system comprising modules in step S1, fall into 5 types altogether, respectively route class, load class,
Confluence class, transformer class, class generator.Classification results as shown in figure 3, establish this 5 class basis modeling unit according to Fig. 3 respectively
A body Model.By taking route individual as an example, a body Model can be indicated are as follows:
Wherein, P, Q, V be respectively active and reactive, voltage, and R, X are the resistance of route, reactance, subscriptBinRepresent route
The input of body, subscriptBoutRepresent the output of route individual, subscriptmin、maxRespectively represent the bound of relevant parameter.
Step S3, a body Model of system modules is inherited out from basic modeling unit model, and establishes each module
Between connection relationship.The electric system 7 determined in step S1 can be inherited out according to a body Model of 5 class basis modeling units
A body Model of a comprising modules, for this sentences route 1, a body Model can be indicated are as follows:
System physical model according to Fig.2, the connection relationship of 7 modules is as shown in Figure 4 in system.
Step S4, system model is established, and each intermodule input/output relation is obtained according to system model.The individual of system
Model can indicate are as follows:
ΣEPS=(VEPS,TEPS) (5)
VEPS={ ΣLoad,1,ΣLoad,2,ΣBranch,1,ΣBranch,2,ΣIS,3,ΣTransformer,ΣGenerator} (6)
According in system modelTEPSEach intermodule input/output relation can be obtained are as follows:
Step S5, according to known parameters, modules are gradually solved, and then obtain the solution of whole system.According to each module
Relationship between input and output, brings known parameters into, can the active power of each node in solution system, idle using forward-backward sweep method
Power, the amplitude of voltage and phase angle are as shown in table 1 below:
The active power, reactive power, the amplitude of voltage and phase angle value table of each node in 1. system of table
The above embodiment is a preferred embodiment of the present invention, but embodiments of the present invention are not by above-described embodiment
Limitation, other any changes, modifications, substitutions, combinations, simplifications made without departing from the spirit and principles of the present invention,
It should be equivalent substitute mode, be included within the scope of the present invention.
Claims (4)
1. a kind of based on the individual integrated energy system modeling method for this model, which is characterized in that the modeling method step
It is rapid as follows:
S1, the physical model for determining integrated energy system take out each comprising modules of integrated energy system, to each composition accordingly
Module naming number；
S2, classify to all modules, determine basis modeling unit, and establish a body Model of basic modeling unit；
S3, a body Model that system modules are inherited out from basic modeling unit model, and establish the company between each module
Connect relationship；
S4, system model is established, and each intermodule input/output relation is obtained according to system model；
S5, according to known parameters, modules are gradually solved using forward-backward sweep method, and then obtain the solution of whole system.
2. according to claim 1 a kind of based on the individual integrated energy system modeling method for this model, feature exists
In a body Model of the basic modeling unit is defined as follows:
A body Model is indicated using penton:
∑=(I, K, X, F, O)
Wherein, I is the input set of individual, and K is the intrinsic Knowledge Set of individual, and X is the state set inside individual, and F is the row of individual
For collection, O is the output collection of individual；
(1) individual is expressed as i in the input of t momentt
it=(i, t) | i ∈ I, t ∈ T }
Wherein, T is the time series set of individual with environmental interaction, and individual input set indicates are as follows:
I={ i1,i2,…,it}；
(2) individual is k in the intrinsic representation of knowledge of t momentt
kt=(k, t) | k ∈ K, t ∈ Tk}
Wherein, TkIt is time series set corresponding to the Knowledge Set of individual, individual Knowledge Set indicates are as follows:
K={ k1,k2,…,kt}；
(3) individual is expressed as X in the state of t momentt
Xt=(X, t) | X ∈ X, t ∈ [0, et]}
Wherein, etIt is the accurate maximum time step-length for capturing system dynamic characteristic, individual state set indicates are as follows:
X={ X1,X2,…,Xt}；
(4) state X of the individual at the t+1 momentt+1More new formula is expressed as
Xt+1=F (Xt,kt,it)
Wherein, F (,) is the behavior equation of individual, describes the specific behavior of individual and evolutionary mechanism, individual is in t+1
The state X at quartert+1, it is individual in one of t moment state mapping, it is related with input set, Knowledge Set, due to the presence of behavior collection,
Individual can also automatically update the state of t moment in the case where no outside stimulus；
(5) individual is expressed as o in the output of t momentt
ot=(o, t) | i ∈ O, t ∈ T }
Wherein, individual output set representations are as follows:
O={ o1,o2,…,ot}。
3. according to claim 1 a kind of based on the individual integrated energy system modeling method for this model, feature exists
In the system model is defined as follows:
M=(V, T)
Wherein, V={ Σ | Σ ∈ V }, V are sub- groups of individuals, and T is the incidence matrix for describing connection relationship between sub- individual,
The constitution element T of incidence matrixijIt is defined as follows:
Tij=+1, work as ΣiOutput be ΣjInput；
Tij=-1, works as ΣjOutput be ΣiInput；
Tij=0, work as ΣiAnd ΣjBetween without direct correlation.
4. according to claim 1 a kind of based on the individual integrated energy system modeling method for this model, feature exists
In steps are as follows for forward-backward sweep method in the step S5:
S501, since the node farthest from power supply point, using route voltage rating, once calculated against the direction of power transmission
Separate the power loss and power distribution in line impedance, wherein power loss Δ S calculation formula are as follows:
Wherein P ", Q " is respectively the active power and reactive power of line end, VNFor the voltage rating of route, R, X are respectively line
The resistance value and reactance value on road, then the trend S ' calculation formula of route front end are as follows:
S '=S "+Δ S
Wherein, S " is the trend of line end；
S502, the power acquired using step S501 calculate partition route along the direction of power transmission since power supply point
Voltage landing finds out each node voltage respectively, wherein the amplitude Δ V of voltage landingmWith the phase angle Δ V of voltage landingaIt counts respectively
It calculates are as follows:
Wherein V ' is the voltage of route front end, and P ', Q ' are respectively the active power and reactive power of route front end, then route is last
The voltage at end is expressed as:
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910000445.2A CN109684763A (en) | 2019-01-02 | 2019-01-02 | It is a kind of based on individual be this model integrated energy system modeling method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910000445.2A CN109684763A (en) | 2019-01-02 | 2019-01-02 | It is a kind of based on individual be this model integrated energy system modeling method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109684763A true CN109684763A (en) | 2019-04-26 |
Family
ID=66190373
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910000445.2A Pending CN109684763A (en) | 2019-01-02 | 2019-01-02 | It is a kind of based on individual be this model integrated energy system modeling method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109684763A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112330493A (en) * | 2020-11-24 | 2021-02-05 | 南方电网科学研究院有限责任公司 | Energy system modeling and comprehensive analysis method, device and storage medium |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103928925A (en) * | 2014-04-17 | 2014-07-16 | 国家电网公司 | Power distribution network load flow calculation method based on forward-backward sweep |
-
2019
- 2019-01-02 CN CN201910000445.2A patent/CN109684763A/en active Pending
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103928925A (en) * | 2014-04-17 | 2014-07-16 | 国家电网公司 | Power distribution network load flow calculation method based on forward-backward sweep |
Non-Patent Citations (1)
| Title |
|---|
| 谭玉华: "个体为本的综合能源系统建模及仿真", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112330493A (en) * | 2020-11-24 | 2021-02-05 | 南方电网科学研究院有限责任公司 | Energy system modeling and comprehensive analysis method, device and storage medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Hua et al. | Optimal energy management strategies for energy Internet via deep reinforcement learning approach | |
| James et al. | Delay aware intelligent transient stability assessment system | |
| Wang et al. | Short-term wind power prediction based on improved small-world neural network | |
| CN109523060A (en) | Ratio optimization method of the high proportion renewable energy under transmission and distribution network collaboration access | |
| Zheng et al. | Distributed real-time dispatch of integrated electricity and heat systems with guaranteed feasibility | |
| Fang et al. | Deep reinforcement learning for scenario-based robust economic dispatch strategy in internet of energy | |
| CN103455948B (en) | A kind of distribution system multi-dimensional multi-resolution Modeling and the method for analysis | |
| Bahmanifirouzi et al. | Multi-objective stochastic dynamic economic emission dispatch enhancement by fuzzy adaptive modified theta particle swarm optimization | |
| Chassin et al. | Load modeling and calibration techniques for power system studies | |
| Hansen et al. | A dynamic market mechanism for markets with shiftable demand response | |
| CN106600145A (en) | Power distribution network reliability analysis method and device | |
| Jain et al. | Adaptive wavelet neural network-based fast dynamic available transfer capability determination | |
| Poluru et al. | Enhancement of ATC by optimizing TCSC configuration using adaptive moth flame optimization algorithm | |
| CN110020743A (en) | A kind of interconnection passway for transmitting electricity method for planning capacity | |
| Sanki et al. | Application of a novel PIPDF controller in an improved plug-in electric vehicle integrated power system for AGC operation | |
| CN109684763A (en) | It is a kind of based on individual be this model integrated energy system modeling method | |
| Hu et al. | Frequency prediction model combining ISFR model and LSTM network | |
| CN110120673A (en) | Distributed transmission & distribution collaboration idle work optimization method and system based on Thevenin's equivalence parameter identification | |
| Vohra et al. | End-to-end learning with multiple modalities for system-optimised renewables nowcasting | |
| KR20190088127A (en) | Method and apparatus for providing transmission loss factor based on power system component modeling and big data processing | |
| Li et al. | A multi-agent deep reinforcement learning-based “Octopus” cooperative load frequency control for an interconnected grid with various renewable units | |
| CN110829411A (en) | Flexible power distribution network operation domain observation method considering alternating current power flow | |
| CN116093995A (en) | Multi-target network reconstruction method and system for power distribution system | |
| Arefi et al. | A novel algorithm based on honey bee mating optimization for distribution harmonic state estimation including distributed generators | |
| Erol | Delay margin computation in micro grid systems with time delay by using fractional order controller |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190426 |