A Joint Time-Series Simulation Method and System for Cyber-Physical Systems of Distribution Networks

Abstract

The present invention provides a method and system for joint time-series simulation of a distribution network cyber-physical system, including: based on the acquired time-series data of the distribution network, performing physical simulation using a pre-established electrical physical system model of the distribution network in the physical sampling period; when the sampling period of the communication network is reached, performing information communication simulation using a pre-established information communication model of the distribution network based on the acquired communication data; splicing the physical simulation and information communication simulation results to obtain the joint simulation result of the distribution network cyber-physical system. This method and system implement a joint solution calculation framework for the distribution network cyber-physical system, reduce the coupling synchronization error between information communication and electrical physical simulation, and solve the problem of online and time-series joint simulation in non-real-time scenarios in the two simulation domains of the distribution network cyber-physical system.

Description

Technical Field

[0001] The present invention belongs to the technical field of cyber-physical systems for distribution networks, and particularly relates to a method and system for joint time-series simulation of cyber-physical systems for distribution networks. Background Art

[0002] In the cyber-physical system for distribution networks, the deep integration and real-time interaction between the information and communication system and the physical system present new characteristics such as complex dynamic process evolution, cross-space fault chain propagation, and sharp increase in risk probability, and face a series of new problems such as information attacks, communication link blockages, and compound fault evolution. Traditional distribution network simulations focus on the primary electrical system of the distribution network, describe the continuous-time response model of the electrical system with a set of differential-algebraic equations, and use numerical integration algorithms such as explicit integration or implicit integration for solution. By discrete time steps, a relatively accurate estimation of the current state of the electrical system is obtained, so as to obtain the time-domain response waveform of the electrical system. The information and communication system belongs to a discrete system, and uses a discrete state model to describe network discrete parameters such as data queue length, data packets, and discrete events, transforms the complex communication process into a specific event queue, and simulates it through a discrete event simulation tool. It can be seen that traditional distribution network simulation and information and communication simulation are basically separated in theory and method. The distribution physical system is a time-varying system, and its electrical quantities flow through power nodes and branches in the form of voltage, current, power, etc. The information and communication system is a discrete system, and its information changes are triggered by discrete events. The existing simulation methods are mainly based on the joint simulation scheme between mature commercial power real-time simulators RTLAB, RTDS, and communication simulation software OPNET. There is currently no solution for online and time-series simulation scenarios for distribution networks under non-real-time conditions. Summary of the Invention

[0003] To overcome the deficiency of the lack of online and time-series simulation solutions in the above-mentioned prior art under non-real-time scenarios, the present invention proposes a method and system for joint time-series simulation of cyber-physical systems for distribution networks. The method and system propose a solution calculation framework for cyber-physical systems for distribution networks based on a joint solution mechanism, and design a simulation synchronization interaction mechanism based on trigger synchronization. By adopting the combination of distribution network electrical physical simulation algorithms and information and communication discrete event simulation, and designing a simulation synchronization mechanism and switching method, the integrated calculation of the information flow of the distribution network communication system and the electrical physical energy flow is realized, solving the problem of the lack of joint simulation methods for cyber-physical systems for distribution networks, and realizing the online and time-series simulation solution of cyber-physical systems for distribution networks under steady-state conditions.

[0004] The solution adopted to achieve the above object is as follows:

[0005] A method for joint time-series simulation of a cyber-physical system for a distribution network, which is improved in that it includes:

[0006] Based on the obtained time-series data of the distribution network, physical simulation is performed using a pre-established electrical physical system model of the distribution network during the physical sampling period;

[0007] When the sampling period of the communication network is reached, based on the obtained communication data, information communication simulation is performed using a pre-established information communication model of the distribution network;

[0008] The physical simulation result and the information communication simulation result are spliced to obtain the joint simulation result of the distribution network cyber-physical system.

[0009] The first preferred technical solution provided by the present invention is improved in that the physical simulation using a pre-established electrical physical system model of the distribution network during the physical sampling period based on the obtained time-series data of the distribution network includes:

[0010] S101: Collect the time-series data of the distribution network at the same time section from the distribution automation system, the power consumption information acquisition system, and the production management system based on the time scale;

[0011] S102: Based on the time-series data and the pre-established electrical physical system model of the distribution network, perform power flow calculation of the distribution network electrical physical system;

[0012] S103: Determine whether the power flow calculation converges: if so, execute S104; otherwise, adjust the line impedance and load parameters of the distribution network electrical physical system model and execute S102 again;

[0013] S104 records the physical simulation result on the event axis including time series in the form of a physical simulation data packet.

[0014] Among them, the time-series data includes: the state data of switches in the distribution network, the active and reactive power data of loads, and the node voltage and current data of the measuring points.

[0015] The second preferred technical solution provided by the present invention is improved in that the determination of whether the power flow calculation converges includes:

[0016] Select a node in the distribution network electrical physical system model, and determine whether the difference between the voltages of the node in two consecutive calculations is within a preset range within the set number of power flow iteration calculations: if so, the node voltage calculation converges, otherwise the node voltage calculation does not converge;

[0017] Determine whether there is a node in the distribution network electrical physical system model whose node voltage calculation does not converge: if so, the power flow calculation does not converge; otherwise the power flow calculation converges.

[0018] The third preferred technical solution provided by the present invention is improved in that the establishment of the distribution network electrical physical system model includes:

[0019] Based on the distribution network topology, a distribution network electrical physical system model is established by using nodal equations and power balance equations;

[0020] The nodal equation is shown as follows:

[0021] I B = Y B U B

[0022] In the formula, I B represents the branch current vector, Y B represents the admittance matrix, and U B represents the nodal voltage vector;

[0023] The power balance equation is shown as follows:

[0024]

[0025]

[0026] In the formula, P i represents the active power injected into electrical physical node i, Q i represents the reactive power injected into electrical physical node i; U i represents the voltage of electrical physical node i, and U j represents the voltage of electrical physical node j; G ij represents the conductance between electrical physical nodes i and j; B ij represents the susceptance between electrical physical nodes i and j; δ ij represents the phase angle difference between electrical physical nodes i and j.

[0027] The fourth preferred technical solution provided by the present invention is improved in that, based on the acquired communication data, information communication simulation is performed by using a pre-established distribution network information communication model, including:

[0028] Bring the acquired communication data into the pre-established distribution network information communication model, and perform information communication simulation on the information communication network with a finite automaton;

[0029] Record the information communication simulation results on the event axis in the form of communication simulation data packets;

[0030] The communication data includes: the status data of switches in the distribution network, the node voltage and current data of the measurement points, and the transformer status data.

[0031] The fifth preferred technical solution provided by the present invention is improved in that, after performing information communication simulation by using the pre-established distribution network information communication model and before splicing the physical simulation and information communication simulation results based on the simulation time, it further includes:

[0032] Determine whether the information communication simulation result contains a control signal:

[0033] If so, update the distribution network topology and operation data according to the control signal, and perform physical simulation again according to the updated distribution network topology and operation data;

[0034] Otherwise, transfer to splicing the physical simulation and the information communication simulation results based on the simulation time.

[0035] The sixth preferred technical solution provided by the present invention is improved in that the establishment of the distribution network information communication model includes:

[0036] Based on the communication network of the distribution network, establish a distribution network information communication model by using a node model and a network model;

[0037] The node model is shown as follows:

[0038] NL lm =(W lm ,Rb lm ,RB lm ), l∈[1,n], m∈[1,n]

[0039] In the formula, n represents the total number of communication network nodes, NL lm is a tuple, W lm represents the bandwidth between communication network nodes l and m, Rb lm represents the transmission rate between communication network nodes l and m, RB lm represents the baud rate between communication network nodes l and m. When l = m, NL lm is a communication node. When l ≠ m, NL lm is a communication link;

[0040] The network model is shown as follows:

[0041]

[0042] A is an n×n adjacency matrix composed of NL lm (l, m∈[1,n]). When NL lm ≠0, it means that communication network nodes l and m are directly connected, and NL lm represents the performance of link lm; when NL lm =0, it means that there is no direct connection between lm.

[0043] The seventh preferred technical solution provided by the present invention is improved in that the splicing of the physical simulation and the information communication simulation results to obtain the joint simulation result of the distribution network cyber-physical system includes:

[0044] Read the physical simulation data packet and the communication simulation data packet from the event axis according to time respectively;

[0045] Stitch the physical simulation data packet and the communication simulation data packet to obtain the joint simulation result of the power distribution network cyber-physical system.

[0046] An improved joint time-sequence simulation system for a power distribution network cyber-physical system, comprising: an electrical physical simulation module, an information communication simulation module, and a data synthesis module;

[0047] The electrical physical simulation module is used to perform physical simulation on the basis of the obtained time-sequence data of the power distribution network by using a pre-established power distribution network electrical physical system model in a physical sampling period;

[0048] The information communication simulation module is used to perform information communication simulation by using a pre-established power distribution network information communication model based on the obtained communication data when the sampling period of the communication network arrives;

[0049] The data synthesis module is used to stitch the physical simulation and information communication simulation results to obtain the joint simulation result of the power distribution network cyber-physical system.

[0050] The improved eighth preferred technical solution provided by the present invention is that the electrical physical simulation module includes: a time-sequence data acquisition unit, a power flow calculation unit, a convergence judgment unit, and a physical data recording unit;

[0051] The time-sequence data acquisition unit is used to collect the time-sequence data of the power distribution network at the same time section from the distribution automation system, the power consumption information acquisition system, and the production management system based on the time scale;

[0052] The power flow calculation unit is used to perform power flow calculation on the power distribution network electrical physical system based on the time-sequence data and a pre-established power distribution network electrical physical system model;

[0053] The convergence judgment unit is used to judge whether the power flow calculation converges: if so, call the physical data recording unit; otherwise, call the power flow calculation unit;

[0054] The physical data recording unit is used to record the physical simulation result on the event axis including time sequence in the form of a physical simulation data packet;

[0055] Wherein, the time-sequence data includes: the state data of switches in the power distribution network, the active and reactive power data of loads, and the node voltage and current data of the measuring points.

[0056] The ninth preferred technical solution provided by the present invention is improved in that the information communication simulation module includes: a finite automaton simulation unit and a communication data recording unit;

[0057] The finite automaton simulation unit is used to bring the obtained communication data into the pre-established distribution network information communication model, and perform information communication simulation on the information communication network with a finite automaton;

[0058] The communication data recording unit is used to record the information communication simulation result on the event axis in the form of a communication simulation data packet;

[0059] The communication data includes: the state data of switches in the distribution network, the node voltage and current data of the measurement points, and the transformer state data.

[0060] The tenth preferred technical solution provided by the present invention is improved in that the data synthesis module includes: a data reading unit and a data splicing unit;

[0061] The data reading unit is used to respectively read the physical simulation data packet and the communication simulation data packet from the event axis according to time;

[0062] The data splicing unit is used to splice the physical simulation data packet and the communication simulation data packet to obtain the joint simulation result of the distribution network information physical system.

[0063] Compared with the closest prior art, the beneficial effects of the present invention are as follows:

[0064] A joint time-series simulation method and system for a distribution network information physical system proposed by the present invention, based on the obtained time-series data of the distribution network, performs physical simulation using the pre-established distribution network electrical physical system model during the physical sampling period; when the sampling period of the communication network arrives, based on the obtained communication data, performs information communication simulation using the pre-established distribution network information communication model; finally, splices the physical simulation and information communication simulation results to obtain the joint simulation result of the distribution network information physical system, realizing the joint solution calculation framework of the distribution network information physical system, reducing the coupling synchronization error between information communication and electrical physical simulation, and solving the online and time-series joint simulation problems in the non-real-time scenarios of the two simulation domains of the distribution network information physical system. Description of the Drawings

[0065] Figure 1 It is a schematic flow chart of a joint time-series simulation method for a distribution network information physical system provided by the present invention;

[0066] Figure 2 It is a schematic diagram of the basic calculation framework of a joint time-series simulation method for a distribution network information physical system provided by the present invention;

[0067] Figure 3 Schematic diagram of the specific process of a joint time-sequential simulation method for a cyber-physical system of a distribution network provided by the present invention;

[0068] Figure 4 Schematic diagram of a synchronization method based on trigger synchronization in a joint time-sequential simulation method for a cyber-physical system of a distribution network provided by the present invention;

[0069] Figure 5 Schematic diagram of the basic structure of a joint time-sequential simulation system for a cyber-physical system of a distribution network provided by the present invention;

[0070] Figure 6 Schematic diagram of the detailed structure of a joint time-sequential simulation system for a cyber-physical system of a distribution network provided by the present invention. Specific embodiments

[0071] The following further elaborates on the specific embodiments of the present invention with reference to the accompanying drawings.

[0072] Example 1:

[0073] A schematic diagram of the process of a joint time-sequential simulation method for a cyber-physical system of a distribution network provided by the present invention is as Figure 1 shown, including:

[0074] 1. A joint time-sequential simulation method for a cyber-physical system of a distribution network, characterized by including:

[0075] Step 1: Based on the obtained time-sequential data of the distribution network, physical simulation is performed using a pre-established electrical physical system model of the distribution network in the physical sampling period;

[0076] Step 2: When the sampling period of the communication network is reached, information communication simulation is performed using a pre-established information communication model of the distribution network based on the obtained communication data;

[0077] Step 3: The physical simulation and information communication simulation results are spliced to obtain the joint simulation result of the cyber-physical system of the distribution network.

[0078] Specifically, the present invention provides a joint time-series simulation method for a distribution network cyber-physical system. First, a model of the distribution network electrical-physical system and node and network models of information communication are established. Secondly, the time-series data of the distribution network at time t is read, and physical simulation is performed on the physical network of the distribution network, that is, power flow calculation is carried out. After the calculation converges, it is judged whether it reaches the sampling period T of the communication network. If it reaches the communication acquisition period T, it will transfer to the information communication simulation subsystem for simulation calculation. Information communication time simulation is carried out through a finite automaton, and a simulation synchronization method based on event time series is adopted between the distribution network information communication and the electrical-physical. After both calculations are completed, the data is interactively spliced to achieve the full-network simulation of the distribution network electrical-physical and information communication. Then, the time-series data of the next moment t + 1 is read for calculation, and through circulation, the joint time-series simulation of the distribution network cyber-physical system is realized.

[0079] The basic calculation framework of the joint time-series simulation method for the distribution network cyber-physical system in this invention patent is as Figure 2 shown. By decomposing the distribution network cyber-physical system into two subsystems, namely the distribution network electrical-physical system and the distribution network information communication, modeling is carried out, and then the electrical system solver and the information communication solver are respectively used for joint solution. A synchronization method based on trigger synchronization is adopted for interactive synchronization in the middle (as Figure 4 shown), and finally the data of the two is synthesized, which is the full-network simulation result of the distribution network cyber-physical system. The specific steps are as Figure 3 shown and are described in detail as follows:

[0080] Step 101: Establish a model of the distribution network cyber-physical system, mainly including the distribution network electrical-physical model and the information communication model, that is, the physical subsystem and the communication subsystem. Among them, the distribution network electrical-physical model is modeled using node equations and power balance equations:

[0081]

[0082] In the formula, I B is the branch current, Y B is the admittance matrix, U B is the node voltage; P i , Q i are the injected powers of the electrical-physical nodes i and j respectively, U i , U j are the voltages of the electrical-physical nodes i and j respectively; G ij represents the conductance between the electrical-physical nodes i and j; B ij represents the susceptance between the electrical-physical nodes i and j; δ ij is the phase angle difference between the electrical-physical nodes i and j.

[0083] The distribution network information communication model is modeled and characterized using node models and network models:

[0084]

[0085] Wherein, NL lm is a tuple; W lm represents the bandwidth between communication network nodes l and m; Rb lm represents the transmission rate between communication network nodes l and m; RB lm represents the baud rate between communication network nodes l and m. When l = m, NL lm is a communication node; when l ≠ m, NL lm is a communication link. A is an n×n adjacency matrix composed of NL lm (l, m ∈ [1, n]), where n is the number of nodes in the communication network system. When NL lm ≠0, it means that communication network nodes l and m are directly connected, and NL lm represents the performance of link lm; when NL lm =0, it means that there is no direct connection between lm.

[0086] Step 102: Read the time-series data at time t, including the online data from the distribution automation system, the power consumption information acquisition system, and the production management system PMS. Ensure the data is at the same time section through time stamps. Among them, the time-series data includes: the status data of switches in the distribution network, the active and reactive power data of loads, and the node voltage and current data of measurement points.

[0087] Step 103: Perform power flow calculation on the distribution network electrical physical system, that is, perform simulation based on the distribution network electrical physical model. The power flow calculation of the distribution network is carried out using the traditional Newton-Raphson algorithm.

[0088] Step 104: Determine whether it converges. If it does not converge, modify the parameters and re-enter Step 3 for calculation; if it converges, transfer to Step 105. Among them, the modified parameters include the line impedance and load parameters of the distribution network electrical physical system model.

[0089] The determination condition for convergence is:

[0090]

[0091] Wherein, is the voltage value of any node i in the distribution network at the current calculation time k, U i (k-1) is the voltage value of node i in the distribution network at the (k - 1)-th calculation, and ε is the set convergence value (such as 0.00001). That is, when the difference between the voltage values of the previous and the next time is small enough, it is determined to converge. Among them, k should be less than the maximum value kmax of the set iteration times.

[0092] The specific convergence determination condition is as follows: Arbitrarily select a node in the electrical physical system model of the distribution network, and judge whether the difference between the voltages of the node in two consecutive calculations is within the preset range ε within the set number of power flow iteration calculations kmax. If so, the voltage calculation of the node converges; otherwise, the voltage calculation of the node does not converge.

[0093] Judge whether there is a node in the electrical physical system model of the distribution network whose voltage calculation does not converge. If so, the power flow calculation does not converge; otherwise, the power flow calculation converges. That is, only when the voltage calculations of all nodes converge is the power flow calculation considered to converge.

[0094] Step 105: Judge whether it reaches the sampling period of communication. If it reaches the sampling period T, go to Step 106 for information communication simulation. If it does not reach the sampling period T, output the simulation result data, read the data at the next moment t + 1, and repeat Step 2.

[0095] Step 106: When it reaches the sampling period T of the information communication system, trigger the communication simulation solver. The communication simulation solver uses a classic finite automaton. Trigger the synchronization mechanism as Figure 4 shown. After the data packet is generated after the physical simulation ends one time step, it is recorded on the virtual event axis and waits for the sampling period T. After reaching the sampling period, trigger the information communication simulation, and realize the simulation calculation of the data on the same section of information communication and electrical physics through the event axis with a sequential relationship, so as to ensure the interaction synchronization between the two simulation solvers.

[0096] Step 107: Receive and parse the data packet, and use a finite automaton to perform simulation calculations on the information communication network.

[0097] Among them, the data packet includes: the state data of switches in the distribution network, the voltage and current data of nodes with measurement points, and the state data of transformers.

[0098] Step 108: After obtaining the information communication simulation results, judge whether these results contain control signals. If they contain control signals, in Step 103, update the network topology and operation data and perform loop calculations again.

[0099] Step 109: After both the electrical physics and information communication solvers complete the calculations, read the simulation results on the event axis. The complete calculation result data of the distribution network information communication system can be obtained through the concatenated results, and the simulation results at the complete moment t can be obtained through data concatenation.

[0100] Step 110: Read the data at the next moment t + 1, go to Step 102, and perform loop calculations to obtain the simulation results of the distribution network cyber-physical system at the moment t + 1. Through the loop of time series, the time series simulation of the distribution network cyber-physical system is realized.

[0101] Embodiment 2:

[0102] Based on the same inventive concept, the present invention also provides a joint time-sequence simulation system for a distribution network cyber-physical system. Since the principles of these devices for solving technical problems are similar to those of the joint time-sequence simulation method for a distribution network cyber-physical system, the repeated parts will not be elaborated here.

[0103] The basic structure of this system is as Figure 5 shown, including: an electrical-physical simulation module, an information-communication simulation module, and a data synthesis module;

[0104] Among them, the electrical-physical simulation module is used to perform physical simulation on the pre-established distribution network electrical-physical system model in the physical sampling period based on the acquired time-sequence data of the distribution network;

[0105] The information-communication simulation module is used to perform information-communication simulation on the pre-established distribution network information-communication model based on the acquired communication data when the sampling period of the communication network arrives;

[0106] The data synthesis module is used to splice the physical simulation and information-communication simulation results to obtain the joint simulation result of the distribution network cyber-physical system.

[0107] The detailed structure of the joint time-sequence simulation system for a distribution network cyber-physical system is as Figure 6 shown.

[0108] Among them, the electrical-physical simulation module includes: a time-sequence data acquisition unit, a power flow calculation unit, a convergence judgment unit, and a physical data recording unit;

[0109] The time-sequence data acquisition unit is used to acquire the time-sequence data of the distribution network at the same time section from the distribution automation system, the power consumption information acquisition system, and the production management system based on the time stamp;

[0110] The power flow calculation unit is used to perform power flow calculation on the distribution network electrical-physical system based on the time-sequence data and the pre-established distribution network electrical-physical system model;

[0111] The convergence judgment unit is used to judge whether the power flow calculation converges: if so, call the physical data recording unit; otherwise, call the power flow calculation unit;

[0112] The physical data recording unit is used to record the physical simulation result in the form of a physical simulation data packet on the event axis including time sequence;

[0113] Among them, the time-sequence data includes: the state data of switches in the distribution network, the active and reactive power data of loads, and the node voltage and current data of the measuring points.

[0114] Among them, the convergence judgment unit includes a node judgment subunit and a whole network judgment subunit;

[0115] The node judgment subunit is used to randomly select nodes in the electrical physical system model of the distribution network, and judge whether the difference between the node voltages in two consecutive calculations is within a preset range within the set number of power flow iteration calculations: if so, the node voltage calculation converges; otherwise, the node voltage calculation does not converge;

[0116] The whole network judgment subunit is used to judge whether there are nodes with non-convergent node voltage calculations in the electrical physical system model of the distribution network: if so, the power flow calculation does not converge; otherwise, the power flow calculation converges.

[0117] Among them, the joint time-sequence simulation system of the distribution network cyber-physical system further includes a physical modeling module for establishing an electrical physical system model of the distribution network;

[0118] The physical modeling module establishes an electrical physical system model of the distribution network based on the distribution network topology structure, using node equations and power balance equations.

[0119] Among them, the information communication simulation module includes: a finite automaton simulation unit and a communication data recording unit;

[0120] The finite automaton simulation unit is used to bring the obtained communication data into the pre-established distribution network information communication model, and perform information communication simulation on the information communication network with a finite automaton;

[0121] The communication data recording unit is used to record the information communication simulation results on the event axis in the form of communication simulation data packets;

[0122] The communication data includes: the state data of switches in the distribution network, the node voltage and current data of nodes with measurement points, and the transformer state data.

[0123] Among them, the joint time-sequence simulation system of the distribution network cyber-physical system further includes a control signal processing module;

[0124] The control signal processing module is used to judge whether the information communication simulation results contain control signals:

[0125] If so, update the distribution network topology structure and operation data according to the control signal, and perform physical simulation again according to the updated distribution network topology structure and operation data;

[0126] Otherwise, call the data synthesis module.

[0127] Among them, the joint time-sequence simulation system of the distribution network cyber-physical system further includes a communication modeling module for establishing a distribution network information communication model;

[0128] The communication modeling module, based on the communication network of the distribution network, establishes the information communication model of the distribution network by using the node model and the network model.

[0129] Among them, the data synthesis module includes: a data reading unit and a data splicing unit;

[0130] The data reading unit is used to read the physical simulation data packet and the communication simulation data packet respectively according to time on the event axis;

[0131] The data splicing unit is used to splice the physical simulation data packet and the communication simulation data packet to obtain the joint simulation result of the information physical system of the distribution network.

[0132] Those skilled in the art should understand that the embodiments of the present application can be provided as methods, systems, or computer program products. Therefore, the present application can adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present application can adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0133] The present application is described with reference to the flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or block in the flowchart and / or block diagram can be realized by computer program instructions, and the combination of the processes and / or blocks in the flowchart and / or block diagram can also be realized by computer program instructions. These computer program instructions can be provided to the processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing devices to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing devices generate a device for realizing the functions specified in Figure 1 one process or multiple processes and / or blocks Figure 1 one block or multiple blocks.

[0134] These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer-readable memory generate a manufactured product including an instruction device, and the instruction device realizes the functions specified in Figure 1 one process or multiple processes and / or blocks Figure 1 one block or multiple blocks.

[0135] These computer program instructions can also be loaded onto a computer or other programmable data processing device, so that a series of operation steps are executed on the computer or other programmable device to generate a computer-implemented process, and thus the instructions executed on the computer or other programmable device provide for realizing the processFigure 1 one process or multiple processes and / or blocks Figure 1 steps of the functions specified in one block or multiple blocks.

[0136] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of this application rather than limit the scope of its protection. Although this application has been described in detail with reference to the above embodiments, those of ordinary skill in the art should understand that after reading this application, various changes, modifications or equivalent substitutions can still be made to the specific implementation manners of the application. However, these changes, modifications or equivalent substitutions are all within the scope of the claims pending approval of this application.

Claims

1. A joint time-series simulation method for a cyber-physical system of a distribution network, characterized in that, Including: Based on the obtained time-series data of the distribution network, physical simulation is carried out using a pre-established electrical physical system model of the distribution network in the physical sampling period; When the sampling period of the communication network arrives, based on the obtained communication data, information communication simulation is carried out using a pre-established information communication model of the distribution network; Stitch the physical simulation and information communication simulation results to obtain the joint simulation result of the distribution network cyber-physical system; The physical simulation carried out using a pre-established electrical physical system model of the distribution network based on the obtained time-series data of the distribution network in the physical sampling period includes: S101: Collect the time-series data of the distribution network at the same time section from the distribution automation system, power consumption information collection system, and production management system based on the time scale; S102: Based on the time-series data and the pre-established electrical physical system model of the distribution network, perform power flow calculation of the electrical physical system of the distribution network; S103: Determine whether the power flow calculation converges: if so, execute S104; otherwise, adjust the line impedance and load parameters of the electrical physical system model of the distribution network, and execute S102 again; S104 Record the physical simulation result in the form of a physical simulation data packet on the event axis containing time series; Wherein, the time-series data includes: the state data of switches in the distribution network, the active and reactive power data of loads, and the node voltage and current data of nodes with measurement points; The establishment of the electrical physical system model of the distribution network includes: Based on the distribution network topology, establish an electrical physical system model of the distribution network using node equations and power balance equations; The node equation is shown as follows: I B = Y B U B where, I B represents the branch current vector, Y B represents the admittance matrix, and U B represents the node voltage vector; The power balance equation is shown as follows: Wherein, P i represents the active power injected by the electrical physical node i, and Q i represents the reactive power injected by the electrical physical node i; U i represents the voltage of the electrical physical node i, and U j represents the voltage of the electrical physical node j; G ij represents the conductance between the electrical physical nodes i and j; B ij represents the susceptance between the electrical physical nodes i and j; δ ij represents the phase angle difference between the electrical physical nodes i and j; The information communication simulation carried out using a pre-established information communication model of the distribution network based on the obtained communication data includes: Bring the obtained communication data into the pre-established information communication model of the distribution network, and perform information communication simulation on the information communication network using a finite automaton; Record the information communication simulation result in the form of a communication simulation data packet on the event axis; The communication data includes: the state data of switches in the distribution network, the node voltage and current data of nodes with measurement points, and the transformer state data; The establishment of the information communication model of the distribution network includes: Based on the communication network of the distribution network, establish an information communication model of the distribution network using node models and network models; The node model is shown as follows: NL lm =(W lm , Rb lm , RB lm ), l ∈ [1, n], m ∈ [1, n] where n represents the total number of communication network nodes, NL lm is a tuple, and W lm represents the bandwidth between communication network nodes l and m, and Rb lm represents the transmission rate between communication network nodes l and m, and RB lm represents the baud rate between communication network nodes l and m. When l = m, NL lm is a communication node. When l ≠ m, NL lm is a communication link; The network model is shown as follows: A is composed of NL lm (where l, m ∈ [1, n]) to form an n×n adjacency matrix. When NL lm ≠0, it indicates that communication network nodes l and m are directly connected, and NL lm represents the performance of link lm; when NL lm =0, it indicates that there is no direct connection between l and m.

2. The method according to claim 1, characterized in that The determination of whether the power flow calculation converges includes: Optionally select a node in the electrical physical system model of the distribution network, and determine whether the difference between the node voltages in two consecutive calculations is within a preset range within the set number of power flow iteration calculations: if so, the node voltage calculation converges, otherwise the node voltage calculation does not converge; Determine whether there are nodes with non-convergent node voltage calculations in the electrical physical system model of the distribution network: if so, the power flow calculation does not converge; otherwise the power flow calculation converges.

3. The method according to claim 1, characterized in that, After the information communication simulation is carried out using a pre-established information communication model of the distribution network and before the physical simulation and information communication simulation results are stitched based on the simulation time, it further includes: Determine whether the information communication simulation result contains a control signal: If so, update the distribution network topology structure and operation data according to the control signal, and perform physical simulation again based on the updated distribution network topology structure and operation data; Otherwise, transfer to splicing the physical simulation and information communication simulation results based on the simulation time.

4. The method according to claim 1, characterized in that Splicing the physical simulation and information communication simulation results to obtain the joint simulation result of the distribution network cyber-physical system, including: Read the physical simulation data packet and communication simulation data packet from the event axis according to time respectively; Splice the physical simulation data packet and the communication simulation data packet to obtain the joint simulation result of the distribution network cyber-physical system.

5. A joint time-series simulation system for a cyber-physical system of a distribution network, characterized in that, Including: An electrical physical simulation module, an information communication simulation module, and a data synthesis module; The electrical physical simulation module is used to perform physical simulation on the pre-established distribution network electrical physical system model in the physical sampling period based on the obtained time series data of the distribution network; The information communication simulation module is used to perform information communication simulation on the pre-established distribution network information communication model based on the obtained communication data when the sampling period of the communication network arrives; The data synthesis module is used to splice the physical simulation and information communication simulation results to obtain the joint simulation result of the distribution network cyber-physical system; The electrical physical simulation module includes: a time series data acquisition unit, a power flow calculation unit, a convergence judgment unit, and a physical data recording unit; The time series data acquisition unit is used to collect the time series data of the distribution network at the same time section from the distribution automation system, the power consumption information acquisition system, and the production management system based on the time scale; The power flow calculation unit is used to perform power flow calculation on the distribution network electrical physical system based on the time series data and the pre-established distribution network electrical physical system model; The convergence judgment unit is used to judge whether the power flow calculation converges: if so, call the physical data recording unit; otherwise, call the power flow calculation unit; The physical data recording unit is used to record the physical simulation result in the form of a physical simulation data packet on the event axis including time series; Wherein, the time series data includes: the state data of switches in the distribution network, the active and reactive power data of loads, and the node voltage and current data of measurement points; The establishment of the distribution network electrical physical system model includes: Based on the distribution network topology structure, establish a distribution network electrical physical system model using node equations and power balance equations; The node equation is shown as follows: I B = Y B U B where I B represents the branch current vector, Y B represents the admittance matrix, and U B represents the node voltage vector; The power balance equation is shown as follows: where P i represents the active power injected by the electrical physical node i, and Q i represents the reactive power injected by the electrical physical node i; U i represents the voltage of the electrical physical node i, and U j represents the voltage of the electrical physical node j; G ij represents the conductance between the electrical physical nodes i and j; B ij represents the susceptance between the electrical physical nodes i and j; δ ij represents the phase angle difference between the electrical physical nodes i and j. The information communication simulation module includes: a finite automaton simulation unit and a communication data recording unit; The finite automaton simulation unit is used to bring the obtained communication data into the pre-established distribution network information communication model and perform information communication simulation on the information communication network with a finite automaton; The communication data recording unit is used to record the information communication simulation result in the form of a communication simulation data packet on the event axis; The communication data includes: the state data of switches in the distribution network, the node voltage and current data of measurement points, and the transformer state data; The establishment of the distribution network information communication model includes: Based on the communication network of the distribution network, establish a distribution network information communication model using node models and network models; The node model is shown as follows: NL lm =(W lm , Rb lm , RB lm ), l ∈ [1, n], m ∈ [1, n] where n represents the total number of communication network nodes, NL lm is a tuple, W lm represents the bandwidth between communication network nodes l and m, Rb lm represents the transmission rate between communication network nodes l and m, RB lm represents the baud rate between communication network nodes l and m. When l = m, NL lm is a communication node. When l ≠ m, NL lm is a communication link; The network model is shown as follows: A is composed of NL lm (l, m ∈ [1, n]) to form an n×n adjacency matrix. When NL lm ≠0, it indicates that communication network nodes l and m are directly connected, and NL lm represents the performance of link lm; when NL lm =0, it indicates that there is no direct connection between l and m.

6. The system according to claim 5, wherein The data synthesis module includes: a data reading unit and a data splicing unit; The data reading unit is used to respectively read physical simulation data packets and communication simulation data packets according to time on the event axis; The data splicing unit is used to splice the physical simulation data packets and the communication simulation data packets to obtain the joint simulation result of the power distribution network cyber-physical system.

Images