CN111680387A - Decomposition method and device of real-time simulation model of active power distribution network - Google Patents
Decomposition method and device of real-time simulation model of active power distribution network Download PDFInfo
- Publication number
- CN111680387A CN111680387A CN202010316298.2A CN202010316298A CN111680387A CN 111680387 A CN111680387 A CN 111680387A CN 202010316298 A CN202010316298 A CN 202010316298A CN 111680387 A CN111680387 A CN 111680387A
- Authority
- CN
- China
- Prior art keywords
- real
- sub
- network
- time
- power distribution
- 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
- 238000009826 distribution Methods 0.000 title claims abstract description 118
- 238000004088 simulation Methods 0.000 title claims abstract description 110
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 96
- 238000000034 method Methods 0.000 title claims abstract description 37
- 239000011159 matrix material Substances 0.000 claims description 24
- 238000004364 calculation method Methods 0.000 claims description 14
- 230000008569 process Effects 0.000 claims description 8
- 238000010200 validation analysis Methods 0.000 claims description 6
- 230000001052 transient effect Effects 0.000 abstract description 4
- 238000004458 analytical method Methods 0.000 abstract description 2
- 238000012827 research and development Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 19
- 238000010248 power generation Methods 0.000 description 11
- 238000004146 energy storage Methods 0.000 description 9
- 238000004590 computer program Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 238000003860 storage Methods 0.000 description 5
- 238000006467 substitution reaction Methods 0.000 description 4
- 238000004891 communication Methods 0.000 description 2
- 239000000446 fuel Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000002699 waste material Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/04—Power grid distribution networks
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Networks & Wireless Communication (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a decomposition method and a decomposition device of a real-time simulation model of an active power distribution network, which are used for determining the number of pre-decomposed networks and pre-decomposed nodes; performing pre-decomposition on the real-time simulation model of the active power distribution network to obtain a plurality of sub-networks, and additionally installing interfaces at pre-decomposition nodes; the sub-networks are distributed to the real-time simulators, and the nodes in the sub-networks are adjusted, so that the real-time simulation model of the active power distribution network can be objectively decomposed, the decomposition time is greatly shortened, the average resource utilization rate of the real-time simulators is considered, resources are saved, the sub-networks obtained by pre-decomposition are adjusted, and the decomposition accuracy is improved; the invention can fully utilize the resource utilization rate of the real-time simulator, improve the real-time simulation speed of the power distribution network to the greatest extent, enlarge the real-time simulation scale of the power distribution network, improve the electromagnetic transient simulation capability of the power distribution network, and provide technical support for operation analysis, equipment research and development, scheduling control and the like of the power distribution network.
Description
Technical Field
The invention relates to the technical field of power system simulation, in particular to a decomposition method and a decomposition device of a real-time simulation model of an active power distribution network.
Background
With the promotion of the construction and the transformation of the power distribution network of a national power grid company, a large number of novel devices such as distributed power supplies, electric vehicles and the like are connected into the power distribution network, and the characteristics of large scale, numerous nodes, complex devices and various operation modes of the power distribution network are increasingly prominent. On one hand, the power distribution network covers various devices such as loads, transformers, circuit breakers, ring main units and the like, and simultaneously comprises a large number of power electronic devices such as distributed power supplies, energy storage elements and static var compensators; on the other hand, the power distribution network is developed into a double-power-supply ring network and multi-power-supply latticed wiring mode from the original single-power-supply radial structure and hand-pull wiring mode, so that the scale of the power distribution network is continuously enlarged, the number of nodes is greatly increased, the complexity degree is increased in an exponential form, and higher requirements and challenges are provided for a real-time simulation technology of a complex power distribution network.
The decomposition method of the real-time simulation model of the active power distribution network is one of effective means for solving the problem of real-time simulation of the power distribution network. The model decomposition is based on the idea of grouping and grouping, the power distribution network is subjected to model decomposition processing, the power distribution network is decomposed into a plurality of sub-networks, one state space system is decomposed into two or more state space groups, and each state space group solves a corresponding state space matrix. The model decomposition and the parallel computation are combined, the computation burden of one processor can be effectively reduced through the multi-processor parallel computation based on the decomposition model, the simulation scale is improved, and the rapid simulation of the complex power distribution network is realized. In the prior art, the real-time simulation model of the active power distribution network with a large scale or containing a large number of power electronic devices is generally decomposed in a manual mode, specifically, the real-time simulation model of the active power distribution network is decomposed by manually additionally installing an interface, and the method has the advantages of strong subjectivity, long consumed time, serious resource waste and inaccurate decomposition result.
Disclosure of Invention
In order to overcome the defects of strong subjectivity, long time consumption, serious resource waste and inaccurate decomposition result in the prior art, the invention provides a decomposition method of a real-time simulation model of an active power distribution network, which comprises the following steps:
determining the number of pre-decomposed networks and pre-decomposed nodes based on the number of real-time simulators and the number of switching devices in the real-time simulation model of the active power distribution network;
pre-decomposing the real-time simulation model of the active power distribution network to obtain a plurality of sub-networks based on the pre-decomposed network number and the pre-decomposition nodes, and additionally installing interfaces at the pre-decomposition nodes;
the sub-network is assigned to each real-time simulator via the interface, and nodes within the sub-network are adjusted based on an average resource utilization of all real-time simulators.
The method for determining the number of pre-decomposed networks and pre-decomposed nodes based on the number of real-time simulators and the number of switching devices in the real-time simulation model of the active power distribution network comprises the following steps:
determining the number of pre-decomposed networks based on the number of distributed power sources in the active power distribution network real-time simulation model and the number of real-time simulators;
determining the number of system matrixes of the active power distribution network according to the number of switching devices in the real-time simulation model of the active power distribution network;
and selecting the number of distributed power supplies contained in each sub-network based on the number of the system matrixes, and selecting a pre-decomposition node based on the number of the distributed power supplies and a preset decomposition principle.
The preset decomposition principle comprises the following steps:
dividing the distributed power supplies of the adjacent nodes into a sub-network;
if the distance between a certain node in the sub-network after the pre-decomposition and other nodes exceeds the preset maximum node distance, dividing the node into other sub-networks;
the maximum node interval is set based on a topological structure of a real-time simulation model of the active power distribution network.
The adjusting the sub-network based on the average resource utilization of all real-time simulators comprises:
determining the resource utilization rate of each real-time simulator according to the process of solving the state space matrix of each sub-network by the real-time simulators;
calculating the average resource utilization rate of all real-time simulators based on the resource utilization rate of each real-time simulator;
and dividing the nodes containing the distributed power supplies in the sub-networks corresponding to the sub-networks with the resource utilization rate larger than the average resource utilization rate into the sub-networks corresponding to the sub-networks with the resource utilization rate smaller than the average resource utilization rate, so that the resource utilization rate of the real-time simulator is as close as possible to the average resource utilization rate.
The average resource utilization of all the real-time simulators is calculated according to the following formula:
in the formula, eavgRepresenting the average resource utilization of all real-time simulators; resource represents the total memory space of all real-time simulators; the validation represents the total computational workload of all real-time simulators,
n represents the number of real-time simulators; utimA memory space representing an mth real-time simulator; e.g. of the typemRepresenting the resource utilization of the mth real-time simulator.
The real-time simulator solving the state space matrix of each sub-network comprises the following steps:
1) each real-time simulator initializes the sub-networks distributed by the real-time simulator to obtain the initial state of each sub-network;
2) each real-time simulator solves the state space matrix of each sub-network based on the initial state of each sub-network to obtain the voltage signal and the current signal of each sub-network;
3) exchanging voltage signals and current signals of adjacent sub-networks based on interface relations and interface types between different sub-networks;
4) and each sub-network updates the state space matrix of the sub-network based on the voltage signal and the current signal obtained by exchange, and returns to 2) to continue solving until the simulation is finished.
The interface comprises a voltage-mode ITM interface or a current-mode ITM interface.
Each real-time simulator initializes the sub-networks allocated to the real-time simulator, and the method comprises the following steps:
if a voltage type ITM interface is adopted, acquiring the initial value of the three-phase current of the opposite side sub-network from the opposite side sub-network, and transmitting the three-phase voltage of the local side sub-network into the opposite side sub-network;
and if a current type ITM interface is adopted, acquiring the initial value of the three-phase voltage of the opposite side sub-network from the opposite side sub-network, and transmitting the three-phase current of the local side sub-network into the opposite side sub-network.
Based on the same invention concept, the invention also provides a decomposition device of the real-time simulation model of the active power distribution network, which comprises the following steps:
the determining module is used for determining the number of pre-decomposed networks and pre-decomposed nodes based on the number of real-time simulators and the number of switching devices in the real-time simulation model of the active power distribution network;
the decomposition module is used for pre-decomposing the real-time simulation model of the active power distribution network to obtain a plurality of sub-networks based on the pre-decomposed network number and the pre-decomposition nodes, and additionally installing interfaces at the pre-decomposition nodes;
and the adjusting module is used for distributing the sub-network to each real-time simulator through the interface and adjusting the nodes in the sub-network based on the average resource utilization rate of all the real-time simulators.
The determining module is specifically configured to:
determining the number of pre-decomposed networks based on the number of distributed power sources in the active power distribution network real-time simulation model and the number of real-time simulators;
determining the number of system matrixes of the active power distribution network according to the number of switching devices in the real-time simulation model of the active power distribution network;
and selecting the number of distributed power supplies contained in each sub-network based on the number of the system matrixes, and selecting a pre-decomposition node based on the number of the distributed power supplies and a preset decomposition principle.
The preset decomposition principle comprises the following steps:
dividing the distributed power supplies of the adjacent nodes into a sub-network;
if the distance between a certain node in the sub-network after the pre-decomposition and other nodes exceeds the preset maximum node distance, dividing the node into other sub-networks;
the maximum node interval is set based on a topological structure of a real-time simulation model of the active power distribution network.
The adjustment module includes:
the determining unit is used for determining the resource utilization rate of each real-time simulator according to the process of solving the state space matrix of each sub-network by the real-time simulator;
the calculating unit is used for calculating the average resource utilization rate of all the real-time simulators based on the resource utilization rate of each real-time simulator;
and the adjusting unit is used for dividing the nodes containing the distributed power supplies in the sub-networks corresponding to the condition that the resource utilization rate of the real-time simulator is greater than the average resource utilization rate into the sub-networks corresponding to the condition that the resource utilization rate of the real-time simulator is less than the average resource utilization rate, so that the resource utilization rate of the real-time simulator is as close as possible to the average resource utilization rate.
The calculating unit calculates the average resource utilization rate of the real-time simulator according to the following formula:
in the formula, eavgRepresenting an average resource utilization of the real-time simulator; resource represents the total memory space of all real-time simulators; the validation represents the total computational workload of all real-time simulators,
n represents the number of real-time simulators; utimA memory space representing an mth real-time simulator; e.g. of the typemRepresenting the resource utilization of the mth real-time simulator.
The computing unit is specifically configured to:
1) each real-time simulator initializes the sub-networks distributed by the real-time simulator to obtain the initial state of each sub-network;
2) each real-time simulator solves the state space matrix of each sub-network based on the initial state of each sub-network to obtain the voltage signal and the current signal of each sub-network;
3) exchanging voltage signals and current signals of adjacent sub-networks based on interface relations and interface types between different sub-networks;
4) and each sub-network updates the state space matrix of the sub-network based on the voltage signal and the current signal obtained by exchange, and returns to 2) to continue solving until the simulation is finished.
The interface comprises a voltage-mode ITM interface or a current-mode ITM interface.
The computing unit is specifically configured to:
if a voltage type ITM interface is adopted, acquiring the initial value of the three-phase current of the opposite side sub-network from the opposite side sub-network, and transmitting the three-phase voltage of the local side sub-network into the opposite side sub-network;
and if a current type ITM interface is adopted, acquiring the initial value of the three-phase voltage of the opposite side sub-network from the opposite side sub-network, and transmitting the three-phase current of the local side sub-network into the opposite side sub-network.
The technical scheme provided by the invention has the following beneficial effects:
in the decomposition method of the real-time simulation model of the active power distribution network, the number of pre-decomposed networks and pre-decomposed nodes are determined based on the number of real-time simulators and the number of switching devices in the real-time simulation model of the active power distribution network; pre-decomposing the real-time simulation model of the active power distribution network to obtain a plurality of sub-networks based on the pre-decomposed network number and the pre-decomposition nodes, and additionally installing interfaces at the pre-decomposition nodes; the sub-networks are distributed to the real-time simulators through the interfaces, the nodes in the sub-networks are adjusted based on the average resource utilization rate of all the real-time simulators, the real-time simulation model of the active power distribution network can be objectively decomposed, the decomposition time is greatly shortened, the average resource utilization rate of the real-time simulators is considered, resources are saved, the sub-networks obtained through pre-decomposition are adjusted, and the decomposition accuracy is improved;
according to the technical scheme provided by the invention, the number of the node admittance matrixes is determined according to the number of the switching devices in the real-time simulation model of the active power distribution network, and the real-time simulation model of the active power distribution network is intelligently decomposed into a plurality of sub-networks, namely, a state space system is automatically decomposed into a plurality of state space groups, and each state space group solves the corresponding state space matrix, so that the calculation efficiency is greatly improved;
according to the technical scheme provided by the invention, each real-time simulator corresponds to one sub-network, so that parallel calculation of a plurality of real-time simulators is realized, the resource utilization rate of the real-time simulators can be fully utilized, the real-time simulation speed of the power distribution network is increased to the greatest extent, and the real-time simulation scale of the power distribution network is enlarged;
the technical scheme provided by the invention improves the electromagnetic transient simulation capability of the power distribution network and provides technical support for operation analysis, equipment research and development, scheduling control and the like of the power distribution network.
Drawings
FIG. 1 is a flowchart of a decomposition method of a real-time simulation model of an active power distribution network according to an embodiment of the present invention;
FIG. 2 is a simplified topology block diagram of a controllable voltage source inverter in an embodiment of the present invention;
FIG. 3 is a diagram of a photovoltaic grid-connected power generation unit in the embodiment of the invention;
FIG. 4 is a schematic diagram of a parallel simulation of multiple real-time simulators in an embodiment of the invention;
FIG. 5 is a diagram of a 10kV voltage class distribution network according to an embodiment of the present invention;
FIG. 6 is a circuit diagram of a real-time simulation model of an active power distribution network according to an embodiment of the present invention;
FIG. 7 is a circuit diagram of an exploded real-time simulation model of an active power distribution network according to an embodiment of the present invention;
FIG. 8 is a diagram of a circuit structure of a differentiated real-time simulation model of an active power distribution network in an embodiment of the present invention;
FIG. 9 is a schematic diagram of a voltage mode ITM interface in an embodiment of the present invention;
fig. 10 is a schematic diagram of a current mode ITM interface in an embodiment of the present invention.
Detailed Description
The present invention will be described in further detail with reference to the accompanying drawings.
Example 1
The embodiment 1 of the invention provides a decomposition method of a real-time simulation model of an active power distribution network, a specific flow chart is shown in fig. 1, and the specific process is as follows:
s101: determining the number of pre-decomposed networks and pre-decomposed nodes based on the number of real-time simulators and the number of switching devices in the real-time simulation model of the active power distribution network;
s102: performing pre-decomposition on the real-time simulation model of the active power distribution network to obtain a plurality of sub-networks based on the number of pre-decomposed networks and pre-decomposition nodes, and additionally installing interfaces at the pre-decomposition nodes;
s103: the subnetworks are assigned to respective real-time simulators through the interfaces, and nodes within the subnetworks are adjusted based on an average resource utilization of all real-time simulators.
The method for determining the number of pre-decomposed networks and pre-decomposed nodes based on the number of real-time simulators and the number of switching devices in the real-time simulation model of the active power distribution network comprises the following steps:
1) determining the number of pre-decomposed networks based on the number of distributed power sources in the active power distribution network real-time simulation model and the number of real-time simulators;
distributed power supplies such as photovoltaic power supplies, fans, fuel cells and micro gas turbines need to be connected into a power grid through power electronic converters so as to solve the problems of energy transfer among different voltage levels, frequencies and alternating current and direct current systems. The demand of the transient model of the equipment elements comprising the power electronic converter on the calculation amount of the real-time simulator is an important factor influencing the number of the pre-decomposition networks.
according to the topology shown in FIG. 2, the following inverter AC side line voltage and DC side voltage U can be obtaineddcThe relationship of (1):
in the above formula, Sa、Sb、ScApparent power, u, of A, B, C phases respectivelyabFor the AC side AB phase line voltage of the inverter, ubcFor the AC side of the inverter, the voltage of the BC phase line ucaIs the AC phase line voltage of the AC side of the inverter.
In the embodiment 1 of the present invention, taking the photovoltaic grid-connected power generation unit shown in fig. 3 as an example, the direct current generated by the photovoltaic array usually needs to be converted into alternating current by the power electronic converter before being connected to the power grid. The photovoltaic grid-connected power generation unit is composed of a photovoltaic array, a power electronic converter, a maximum power controller and a grid-connected controller, and is shown in the attached figure 3. In FIG. 3, PmppIrradiance is 1kW/m at a certain temperature2A reference value of the output power of the photovoltaic array; t is the current temperature, IrrIs the current irradiance, FTFor the current output power coefficient, EFFFor the current efficiency coefficient, PPVIs the output power of the photovoltaic array, SPVBeing photovoltaic arraysApparent power, PnAnd QnRespectively the active power and the reactive power output by the inverter.
The quantity of switching devices of the power electronic converter is large, the switching states are mutually coupled, a large amount of pre-calculation needs a large storage space and high calculation capacity of a real-time simulator, and difficulty is brought to real-time of a model.
Transient simulation of the power distribution network is to solve a differential-algebraic equation set, the resolving scale is large, and the resolving speed is low. In order to realize real-time simulation, a multi-real-time simulator simulation mechanism is adopted, as shown in fig. 4, and in fig. 4, the RTDS is a real-time simulator. The method comprises the steps that a network to be simulated is decomposed into a plurality of sub-networks, the sub-networks are connected through communication lines, each sub-network is distributed to one real-time simulator, the resolving results of the real-time simulators are shared through the communication lines, the resolving scale of the large-scale complex power distribution network is reduced, and system-level parallel simulation is achieved.
A structure diagram of a 10kV voltage class power distribution network is shown in fig. 5, UPOC is a reactive compensator, a 10kV voltage class power distribution network is composed of a power distribution network and a plurality of distributed power generation units (photovoltaic, wind power, gas turbine, fuel cell, etc.), energy storage and load connected to the network, and the power distribution network and a large number of distributed power generation (energy storage) units should be resolved in consideration of the large resolving scale of a single distributed power generation (energy storage) unit. When the number of the real-time simulators is not less than 4, the number n of the pre-decomposed networks is taken as 4; when the number of the real-time simulators is less than 4, the number n of the pre-decomposed networks is taken as the number of the real-time simulators.
2) Determining the number of system matrixes of the active power distribution network according to the number of switching devices in the real-time simulation model of the active power distribution network;
3) and selecting the number of distributed power supplies contained in each sub-network based on the number of the system matrixes, and selecting a pre-decomposition node based on the number of the distributed power supplies and a preset decomposition principle.
The decomposition uniformity and the decomposition node number of the real-time simulation model of the active power distribution network directly influence the speed of parallel computation, and are important factors influencing the simulation parallel performance. And decomposing the real-time simulation model of the active power distribution network into n sub-networks according to the pre-decomposition nodes.
When the real-time simulation is carried out, the real-time simulator can pre-calculate each switch state and store the node admittance matrix obtained by the pre-calculation. For the distributed power generation (energy storage) unit using the voltage-type inverter shown in fig. 2, the number of switching devices included in each unit is 6, and the number of corresponding node admittance matrixes is 2664.
Taking a power distribution network with a certain voltage class of 10kV as shown in fig. 5 as an example, the total number Δ of distributed power generation (energy storage) units in the real-time simulation model of the active power distribution network is 8. When the active power distribution network real-time simulation model is not decomposed and is used as a state space system for resolving, x (t + delta t) ═ Akx(t)+Bku (t + Δ t), k 1,2, where x, u are the state variables and the input vector, respectively, ak，BkRespectively, the state matrixes corresponding to the kth switching sequence. The value of the number k of the node admittance matrixes obtained by the pre-calculation of the real-time simulator reaches 248Therefore, a huge storage space and a large amount of calculation are required for the real-time simulator, which may cause that the real-time simulator cannot complete all calculations within one step length, and finally cannot perform real-time simulation on the real-time simulation model of the active power distribution network.
And decomposing the active power distribution network real-time simulation model, wherein the number of distributed power generation (energy storage) units contained in each sub-network is 8/4-2, and the decomposition result is shown as a dotted line box in fig. 5. Decomposing an active power distribution network real-time simulation model into 4 sub-networks, and performing parallel computation by using 4 real-time simulators:
x1(t+Δt)＝Ak1x1(t)+Bk1u1(t+Δt),k1＝1,2,...,212
x2(t+Δt)＝Ak2x2(t)+Bk2u2(t+Δt),k2＝1,2,...,212
x3(t+Δt)＝Ak3x3(t)+Bk3u3(t+Δt),k3＝1,2,...,212
x4(t+Δt)＝Ak4x4(t)+Bk4u4(t+Δt),k4＝1,2,...,212
4 sub-networks are respectively calculated on different real-time simulators, and the number of node admittance matrixes obtained by pre-calculation is 212The burden of the real-time simulator during operation is reduced, and the simulation scale is enlarged.
And counting the total number of the distributed power generation (energy storage) units in the power distribution network to be delta. The method comprises the following steps of performing pre-decomposition on a real-time simulation model of the active power distribution network according to a preset decomposition principle, setting the number of distributed power generation (energy storage) units contained in each sub-network to be delta/n, wherein the preset decomposition principle comprises the following steps: 1) dividing the distributed power supplies of the adjacent nodes into a sub-network; 2) if the distance between a certain node in the sub-network after the pre-decomposition and other nodes exceeds the preset maximum node distance, dividing the node into other sub-networks; the maximum node interval is set based on the topological structure of the real-time simulation model of the active power distribution network.
The interface additionally arranged at the pre-decomposition node comprises a voltage type ITM interface or a current type ITM interface, the interface in the embodiment of the invention adopts an Ideal Transformer Model (ITM) interface, the ITM is based on a substitution theorem, a controlled voltage source and a controlled current source are used as signal receiving devices, and voltage or current signals on the opposite side of the interface are received, so that the decomposition of the real-time simulation model of the active power distribution network is realized.
The embodiment of the invention carries out pre-decomposition on the real-time simulation model of the active power distribution network shown in fig. 6, the circuit structure diagram of the decomposed real-time simulation model of the active power distribution network is shown in fig. 7, and the real-time simulation model of the active power distribution network is decomposed into a sub-network 1 and a sub-network 2 at the dotted line. Based on the substitution theorem, a controlled current source is equivalent to the sub-network 2 in the sub-network 1, and the current of the controlled current source is equal to the measured line current i of the sub-network 2; sub-network 1 is equivalent to sub-network 2 by a controlled voltage source with a voltage equal to the measured interface voltage u of sub-network 1. Only 2 variables of the sub-network are required with ITM: (1) the interface voltage u of the sub-network 1 is input as a control signal of a controlled voltage source in the sub-network 2; (2) the line current i of sub-network 2 is input as a control signal for the controlled current source in sub-network 1. Since the controlled current source cannot be directly connected in series with the inductive element and the controlled voltage source cannot be directly connected in parallel with the capacitive element, it needs to be differentially processed.
The difference equation of the inductance volt-ampere characteristic obtained by adopting the implicit trapezoidal integration method is as follows:
wherein,
similarly, the difference equation of the capacitance volt-ampere characteristic can be obtained as follows:
wherein,
the differentiated post-ITM interface is shown in fig. 8. The ITM interface is divided into a voltage-mode ITM interface (shown in fig. 9) and a current-mode ITM interface (shown in fig. 10), E, according to the type of the interfacel(t)、E2(t) Thevenin equivalent voltages, Z, for sub-network 1 and sub-network 2, respectivelyl、Z2The thevenin equivalent impedances of sub-network 1 and sub-network 2, respectively. For voltage-mode ITM interfaces, the interface voltage u of the sub-network 11After a step length delay, the interface current i is transmitted to the controlled voltage source control end of the sub-network 2, and the interface current i of the sub-network 22Directly transmitting the current to a controlled current source control end of the sub-network 1; for current mode ITM interfaces, the interface voltage u of the sub-network 22Directly transmitted to the controlled voltage source control terminal of the sub-network 1, the interface current i of the sub-network 11And transmitting the delayed signal to the current source control end of the sub-network 2 after one step length.
Adjusting the sub-network based on the average resource utilization of all real-time simulators, comprising:
determining the resource utilization rate of each real-time simulator according to the process of solving the state space matrix of each sub-network by the real-time simulators;
calculating the average resource utilization rate of all real-time simulators based on the resource utilization rate of each real-time simulator;
and dividing the nodes containing the distributed power supplies in the sub-networks corresponding to the sub-networks with the resource utilization rate larger than the average resource utilization rate into the sub-networks corresponding to the sub-networks with the resource utilization rate smaller than the average resource utilization rate, so that the resource utilization rate of the real-time simulator is close to the average resource utilization rate as much as possible and is kept in a proper range, and the resource utilization rate is neither too high nor too low.
The average resource utilization of all real-time simulators is calculated as follows:
in the formula, eavgRepresenting an average resource utilization of the real-time simulator; resource represents the total memory space of all real-time emulators,
the validation represents the total computational workload of all real-time simulators,
n represents the number of real-time simulators; utimA memory space representing an mth real-time simulator; e.g. of the typemRepresenting the resource utilization of the mth real-time simulator.
The real-time simulator solves the state space matrix of each sub-network, and the method comprises the following steps:
1) each real-time simulator initializes the sub-networks distributed by the real-time simulator to obtain the initial state of each sub-network;
2) each real-time simulator solves the state space matrix of each sub-network based on the initial state of each sub-network to obtain the voltage signal and the current signal of each sub-network;
3) exchanging voltage signals and current signals of adjacent sub-networks based on interface relations and interface types between different sub-networks;
4) and each sub-network updates the state space matrix of the sub-network based on the voltage signal and the current signal obtained by exchange, and returns to 2) to continue solving until the simulation is finished.
Each real-time simulator initializes the sub-network allocated to itself, and the method comprises the following steps:
if a voltage type ITM interface is adopted, acquiring the initial value of the three-phase current of the opposite side sub-network from the opposite side sub-network, and transmitting the three-phase voltage of the local side sub-network into the opposite side sub-network;
and if a current type ITM interface is adopted, acquiring the initial value of the three-phase voltage of the opposite side sub-network from the opposite side sub-network, and transmitting the three-phase current of the local side sub-network into the opposite side sub-network.
Example 2
Based on the same inventive concept, embodiment 2 of the present invention further provides a decomposition apparatus for a real-time simulation model of an active power distribution network, including:
the determining module is used for determining the number of pre-decomposed networks and pre-decomposed nodes based on the number of real-time simulators and the number of switching devices in the real-time simulation model of the active power distribution network;
the decomposition module is used for pre-decomposing the real-time simulation model of the active power distribution network to obtain a plurality of sub-networks based on the pre-decomposed network number and pre-decomposition nodes, and additionally installing interfaces at the pre-decomposition nodes;
and the adjusting module is used for distributing the sub-networks to each real-time simulator through the interfaces and adjusting the nodes in the sub-networks based on the average resource utilization rate of all the real-time simulators.
The determination module is specifically configured to:
determining the number of pre-decomposed networks based on the number of distributed power sources in the active power distribution network real-time simulation model and the number of real-time simulators;
determining the number of system matrixes of the active power distribution network according to the number of switching devices in the real-time simulation model of the active power distribution network;
and selecting the number of distributed power supplies contained in each sub-network based on the number of the system matrixes, and selecting a pre-decomposition node based on the number of the distributed power supplies and a preset decomposition principle.
The preset decomposition principle comprises the following steps:
dividing the distributed power supplies of the adjacent nodes into a sub-network;
if the distance between a certain node in the sub-network after the pre-decomposition and other nodes exceeds the preset maximum node distance, dividing the node into other sub-networks;
the maximum node interval is set based on the topological structure of the real-time simulation model of the active power distribution network.
The adjustment module includes:
the determining unit is used for determining the resource utilization rate of each real-time simulator according to the process of solving the state space matrix of each sub-network by the real-time simulator;
the calculating unit is used for calculating the average resource utilization rate of all the real-time simulators based on the resource utilization rate of each real-time simulator;
and the adjusting unit is used for dividing the nodes containing the distributed power supplies in the sub-networks corresponding to the condition that the resource utilization rate of the real-time simulator is greater than the average resource utilization rate into the sub-networks corresponding to the condition that the resource utilization rate of the real-time simulator is less than the average resource utilization rate, so that the resource utilization rate of the real-time simulator is as close as possible to the average resource utilization rate.
The calculating unit calculates the average resource utilization rate of all real-time simulators according to the following formula:
in the formula, eavgRepresenting the average resource utilization of all real-time simulators; resource represents the total memory space of all real-time simulators; the validation represents the total computational workload of all real-time simulators,
n represents the number of real-time simulators; utimA memory space representing an mth real-time simulator; e.g. of the typemRepresenting resources of the mth real-time simulatorSource utilization.
The calculation unit is specifically configured to:
1) each real-time simulator initializes the sub-networks distributed by the real-time simulator to obtain the initial state of each sub-network;
2) each real-time simulator solves the state space matrix of each sub-network based on the initial state of each sub-network to obtain the voltage signal and the current signal of each sub-network;
3) exchanging voltage signals and current signals of adjacent sub-networks based on interface relations and interface types between different sub-networks;
4) and each sub-network updates the state space matrix of the sub-network based on the voltage signal and the current signal obtained by exchange, and returns to 2) to continue solving until the simulation is finished.
The interface comprises a voltage-mode ITM interface or a current-mode ITM interface.
The calculation unit is specifically configured to:
if a voltage type ITM interface is adopted, acquiring the initial value of the three-phase current of the opposite side sub-network from the opposite side sub-network, and transmitting the three-phase voltage of the local side sub-network into the opposite side sub-network;
and if a current type ITM interface is adopted, acquiring the initial value of the three-phase voltage of the opposite side sub-network from the opposite side sub-network, and transmitting the three-phase current of the local side sub-network into the opposite side sub-network.
For convenience of description, each part of the above apparatus is separately described as being functionally divided into various modules or units. Of course, the functionality of the various modules or units may be implemented in the same one or more pieces of software or hardware when implementing the present application.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
Finally, it should be noted that: the above embodiments are only intended to illustrate the technical solution of the present invention and not to limit the same, and a person of ordinary skill in the art can make modifications or equivalent substitutions to the specific embodiments of the present invention with reference to the above embodiments, and any modifications or equivalent substitutions which do not depart from the spirit and scope of the present invention are within the protection scope of the present invention as claimed in the appended claims.
Claims (16)
1. A decomposition method of a real-time simulation model of an active power distribution network is characterized by comprising the following steps:
determining the number of pre-decomposed networks and pre-decomposed nodes based on the number of real-time simulators and the number of switching devices in the real-time simulation model of the active power distribution network;
pre-decomposing the real-time simulation model of the active power distribution network to obtain a plurality of sub-networks based on the pre-decomposed network number and the pre-decomposition nodes, and additionally installing interfaces at the pre-decomposition nodes;
the sub-network is assigned to each real-time simulator via the interface, and nodes within the sub-network are adjusted based on an average resource utilization of all real-time simulators.
2. The decomposition method for the real-time simulation model of the active power distribution network according to claim 1, wherein the determining the number of pre-decomposed networks and the number of pre-decomposition nodes based on the number of real-time simulators and the number of switching devices in the real-time simulation model of the active power distribution network comprises:
determining the number of pre-decomposed networks based on the number of distributed power sources in the active power distribution network real-time simulation model and the number of real-time simulators;
determining the number of system matrixes of the active power distribution network according to the number of switching devices in the real-time simulation model of the active power distribution network;
and selecting the number of distributed power supplies contained in each sub-network based on the number of the system matrixes, and selecting a pre-decomposition node based on the number of the distributed power supplies and a preset decomposition principle.
3. The decomposition method of the real-time simulation model of the active power distribution network according to claim 2, wherein the decomposition principle comprises:
dividing the distributed power supplies of the adjacent nodes into a sub-network;
if the distance between a certain node in the sub-network after the pre-decomposition and other nodes exceeds the preset maximum node distance, dividing the node into other sub-networks;
the maximum node interval is set based on a topological structure of a real-time simulation model of the active power distribution network.
4. The decomposition method for the real-time simulation model of the active power distribution network according to claim 1, wherein the adjusting the nodes in the sub-network based on the average resource utilization of all real-time simulators comprises:
determining the resource utilization rate of each real-time simulator according to the process of solving the state space matrix of each sub-network by the real-time simulators;
calculating the average resource utilization rate of all real-time simulators based on the resource utilization rate of each real-time simulator;
and dividing the nodes containing the distributed power supplies in the sub-networks corresponding to the sub-networks with the resource utilization rate larger than the average resource utilization rate into the sub-networks corresponding to the sub-networks with the resource utilization rate smaller than the average resource utilization rate, so that the resource utilization rate of the real-time simulator is as close as possible to the average resource utilization rate.
5. The decomposition method for the real-time simulation model of the active power distribution network according to claim 4, wherein the average resource utilization of all the real-time simulators is calculated according to the following formula:
in the formula, eavgRepresenting the average resource utilization of all real-time simulators; resource represents the total memory space of all real-time simulators; the validation represents the total computational workload of all real-time simulators,
n represents the number of real-time simulators; utimA memory space representing an mth real-time simulator; e.g. of the typemRepresenting the resource utilization of the mth real-time simulator.
6. The decomposition method for the real-time simulation model of the active power distribution network according to claim 4, wherein the solving of the state space matrix of each sub-network by the real-time simulator comprises:
1) each real-time simulator initializes the sub-networks distributed by the real-time simulator to obtain the initial state of each sub-network;
2) each real-time simulator solves the state space matrix of each sub-network based on the initial state of each sub-network to obtain the voltage signal and the current signal of each sub-network;
3) exchanging voltage signals and current signals of adjacent sub-networks based on interface relations and interface types between different sub-networks;
4) and each sub-network updates the state space matrix of the sub-network based on the voltage signal and the current signal obtained by exchange, and returns to 2) to continue solving until the simulation is finished.
7. The decomposition method for the real-time simulation model of the active power distribution network according to claim 6, wherein the interface comprises a voltage-mode ITM interface or a current-mode ITM interface.
8. The decomposition method for the real-time simulation model of the active power distribution network according to claim 7, wherein each real-time simulator initializes its own assigned sub-network, comprising:
if a voltage type ITM interface is adopted, acquiring the initial value of the three-phase current of the opposite side sub-network from the opposite side sub-network, and transmitting the three-phase voltage of the local side sub-network into the opposite side sub-network;
and if a current type ITM interface is adopted, acquiring the initial value of the three-phase voltage of the opposite side sub-network from the opposite side sub-network, and transmitting the three-phase current of the local side sub-network into the opposite side sub-network.
9. The utility model provides a decomposition device of real-time simulation model of active power distribution network which characterized in that includes:
the determining module is used for determining the number of pre-decomposed networks and pre-decomposed nodes based on the number of real-time simulators and the number of switching devices in the real-time simulation model of the active power distribution network;
the decomposition module is used for pre-decomposing the real-time simulation model of the active power distribution network to obtain a plurality of sub-networks based on the pre-decomposed network number and the pre-decomposition nodes, and additionally installing interfaces at the pre-decomposition nodes;
and the adjusting module is used for distributing the sub-network to each real-time simulator through the interface and adjusting the nodes in the sub-network based on the average resource utilization rate of all the real-time simulators.
10. The decomposition device for the real-time simulation model of the active power distribution network according to claim 10, wherein the determination module is specifically configured to:
determining the number of pre-decomposed networks based on the number of distributed power sources in the active power distribution network real-time simulation model and the number of real-time simulators;
determining the number of system matrixes of the active power distribution network according to the number of switching devices in the real-time simulation model of the active power distribution network;
and selecting the number of distributed power supplies contained in each sub-network based on the number of the system matrixes, and selecting a pre-decomposition node based on the number of the distributed power supplies and a preset decomposition principle.
11. The decomposition device for the real-time simulation model of the active power distribution network according to claim 10, wherein the preset decomposition principle comprises:
dividing the distributed power supplies of the adjacent nodes into a sub-network;
if the distance between a certain node in the sub-network after the pre-decomposition and other nodes exceeds the preset maximum node distance, dividing the node into other sub-networks;
the maximum node interval is set based on a topological structure of a real-time simulation model of the active power distribution network.
12. The decomposition device for the real-time simulation model of the active power distribution network according to claim 9, wherein the adjusting module comprises:
the determining unit is used for determining the resource utilization rate of each real-time simulator according to the process of solving the state space matrix of each sub-network by the real-time simulator;
the calculating unit is used for calculating the average resource utilization rate of all the real-time simulators based on the resource utilization rate of each real-time simulator;
and the adjusting unit is used for dividing the nodes containing the distributed power supplies in the sub-networks corresponding to the condition that the resource utilization rate of the real-time simulator is greater than the average resource utilization rate into the sub-networks corresponding to the condition that the resource utilization rate of the real-time simulator is less than the average resource utilization rate, so that the resource utilization rate of the real-time simulator is as close as possible to the average resource utilization rate.
13. The decomposition method for the real-time simulation model of the active power distribution network according to claim 12, wherein the calculating unit calculates the average resource utilization rate of the real-time simulator according to the following formula:
in the formula, eavgRepresenting an average resource utilization of the real-time simulator; resource represents the total memory space of all real-time simulators; the validation represents the total computational workload of all real-time simulators,
n represents the number of real-time simulators; utimA memory space representing an mth real-time simulator; e.g. of the typemRepresenting the resource utilization of the mth real-time simulator.
14. The decomposition device for the real-time simulation model of the active power distribution network according to claim 12, wherein the calculation unit is specifically configured to:
1) each real-time simulator initializes the sub-networks distributed by the real-time simulator to obtain the initial state of each sub-network;
2) each real-time simulator solves the state space matrix of each sub-network based on the initial state of each sub-network to obtain the voltage signal and the current signal of each sub-network;
3) exchanging voltage signals and current signals of adjacent sub-networks based on interface relations and interface types between different sub-networks;
4) and each sub-network updates the state space matrix of the sub-network based on the voltage signal and the current signal obtained by exchange, and returns to 2) to continue solving until the simulation is finished.
15. The decomposition device for the real-time simulation model of the active power distribution network according to claim 14, wherein the interface comprises a voltage-mode ITM interface or a current-mode ITM interface.
16. The decomposition device for the real-time simulation model of the active power distribution network according to claim 15, wherein the calculation unit is specifically configured to:
if a voltage type ITM interface is adopted, acquiring the initial value of the three-phase current of the opposite side sub-network from the opposite side sub-network, and transmitting the three-phase voltage of the local side sub-network into the opposite side sub-network;
and if a current type ITM interface is adopted, acquiring the initial value of the three-phase voltage of the opposite side sub-network from the opposite side sub-network, and transmitting the three-phase current of the local side sub-network into the opposite side sub-network.
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010316298.2A CN111680387A (en) | 2020-04-21 | 2020-04-21 | Decomposition method and device of real-time simulation model of active power distribution network |
| PCT/CN2020/113692 WO2021212724A1 (en) | 2020-04-21 | 2020-09-07 | Decomposition method and device for real-time simulation model of active power distribution network |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010316298.2A CN111680387A (en) | 2020-04-21 | 2020-04-21 | Decomposition method and device of real-time simulation model of active power distribution network |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111680387A true CN111680387A (en) | 2020-09-18 |
Family
ID=72451734
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010316298.2A Pending CN111680387A (en) | 2020-04-21 | 2020-04-21 | Decomposition method and device of real-time simulation model of active power distribution network |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN111680387A (en) |
| WO (1) | WO2021212724A1 (en) |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7321834B2 (en) * | 2005-07-15 | 2008-01-22 | Chang Gung University | Method for calculating power flow solution of a power transmission network that includes interline power flow controller (IPFC) |
| CN106777636B (en) * | 2016-12-07 | 2020-08-07 | 天津大学 | Block parallel method for micro-grid electromagnetic transient real-time simulation |
| CN108063442A (en) * | 2017-12-30 | 2018-05-22 | 贵州大学 | A kind of electric power system alternating current power grid real-time simulation apparatus and its emulation mode |
| CN108256217A (en) * | 2018-01-17 | 2018-07-06 | 广东电网有限责任公司电力科学研究院 | Numerical model analysis simulated power interface algorithm and device based on ideal transformer method |
| CN109829178B (en) * | 2018-10-26 | 2020-04-14 | 贵州电网有限责任公司 | Transformer substation real-time simulation method based on model segmentation method |
| CN110471309A (en) * | 2019-07-26 | 2019-11-19 | 西北工业大学 | For power hardware loop system emulation mode |
-
2020
- 2020-04-21 CN CN202010316298.2A patent/CN111680387A/en active Pending
- 2020-09-07 WO PCT/CN2020/113692 patent/WO2021212724A1/en active Application Filing
Also Published As
| Publication number | Publication date |
|---|---|
| WO2021212724A1 (en) | 2021-10-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Morstyn et al. | Model predictive control for distributed microgrid battery energy storage systems | |
| CN107171341B (en) | Integrated reactive power optimization method for power transmission and distribution network based on distributed computation | |
| Li et al. | Dynamic equivalent modeling of two-staged photovoltaic power station clusters based on dynamic affinity propagation clustering algorithm | |
| CN103777525A (en) | User-defined interface of wind power station simulation machine and RTDS (real time digital simulator) | |
| Prokhorov et al. | Real time hybrid simulation of electrical power systems: concept, tools, field experience and smart grid challenges | |
| CN108347057A (en) | A kind of AC and DC power system tidal current computing method of the power grid of mixed DC containing LCC-MMC | |
| CN103605293B (en) | A kind of substation relay protection emulation test system and emulation test method thereof | |
| Weiwei et al. | An implementation technology of electromagnetic transient real-time simulation for large-scale grid based on HYPERSIM | |
| CN115481856A (en) | Comprehensive energy system multi-scale scheduling method and system considering comprehensive demand response | |
| CN106021754B (en) | Consider the serial-parallel power grid Probabilistic Load Flow algorithm of VSC reactive power constraints adjustable strategies | |
| Lu et al. | Autonomous power management and load sharing in isolated micro-grids by consensus-based droop control of power converters | |
| CN112072692A (en) | Impedance equivalence method and device for new energy power generation station | |
| CN111934307B (en) | Flat operation control method and system for direct current power distribution network | |
| CN109617086A (en) | A kind of power grid dynamic passive compensation equipment coordination control system and method | |
| CN109245096A (en) | A kind of active distribution network net capability calculation method | |
| CN110676836B (en) | Parallel dynamic compensation method and system of BUCK converter based on disturbance observer | |
| CN107508318B (en) | Active control method and system based on voltage sensitivity partition | |
| CN111680387A (en) | Decomposition method and device of real-time simulation model of active power distribution network | |
| CN112783002A (en) | Digital-analog hybrid simulation method and system for direct-current power distribution network | |
| Wang et al. | Performance evaluation of distributed energy resource management algorithm in large distribution networks | |
| Eissa et al. | A novel approach for optimum allocation of Flexible AC Transmission Systems using Harmony Search technique | |
| Yousefi et al. | Voltage control approach based on PCPM distributed algorithm in the presence of high PV penetration: a stochastic modeling | |
| Hao et al. | Research of distribution network reconfiguration with renewable energy power generation unit | |
| Yu et al. | Reactive Voltage Control in Photovoltaic Power Stations Considering Source-Load Uncertainty | |
| Xu et al. | Research on Electromechanical-electromagnetic Hybrid Simulation Algorithm for Large Scale Power System based on Boundary Nodes Grouping and Decoupling |
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 |