CN112510716B - Power flow calculation method and device of power supply system, storage medium and electronic device - Google Patents

Power flow calculation method and device of power supply system, storage medium and electronic device Download PDF

Info

Publication number
CN112510716B
CN112510716B CN202011425402.8A CN202011425402A CN112510716B CN 112510716 B CN112510716 B CN 112510716B CN 202011425402 A CN202011425402 A CN 202011425402A CN 112510716 B CN112510716 B CN 112510716B
Authority
CN
China
Prior art keywords
target
voltage
train
node
supply system
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.)
Active
Application number
CN202011425402.8A
Other languages
Chinese (zh)
Other versions
CN112510716A (en
Inventor
席嫣娜
王方敏
李占赫
李伟
宋宝同
沈卓轩
李笑倩
宝海龙
李笑彤
魏应冬
陆超
王思涵
刘新萌
张利
刘兆燕
王婷婷
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua University
State Grid Corp of China SGCC
State Grid Beijing Electric Power Co Ltd
Economic and Technological Research Institute of State Grid Beijing Electric Power Co Ltd
Original Assignee
Tsinghua University
State Grid Corp of China SGCC
State Grid Beijing Electric Power Co Ltd
Economic and Technological Research Institute of State Grid Beijing Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tsinghua University, State Grid Corp of China SGCC, State Grid Beijing Electric Power Co Ltd, Economic and Technological Research Institute of State Grid Beijing Electric Power Co Ltd filed Critical Tsinghua University
Priority to CN202011425402.8A priority Critical patent/CN112510716B/en
Publication of CN112510716A publication Critical patent/CN112510716A/en
Application granted granted Critical
Publication of CN112510716B publication Critical patent/CN112510716B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/04Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
    • H02J3/06Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/10Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The embodiment of the invention provides a power flow calculation method, a power flow calculation device, a storage medium and an electronic device of a power supply system, wherein the method comprises the following steps: initializing the node voltage in the power supply system according to the target time; after the initialization is completed, a target operation is executed, the target operation comprising the processing steps of: repeatedly performing predetermined processing on the power supply system until an absolute value of a difference value of a first voltage and a second voltage of a first target node included in the power supply system is smaller than a first maximum allowable error; under the condition that the absolute value of the difference value of the first voltage and the second voltage is smaller than the first maximum allowable error, judging whether the target time reaches the end time or not; in the case where it is determined that the target time reaches the end time, a result of performing the predetermined processing is determined as a power flow calculation result. The invention solves the problem that the power supply system can not be accurately analyzed and calculated in the related technology.

Description

Power flow calculation method and device of power supply system, storage medium and electronic device
Technical Field
The embodiment of the invention relates to the technical field of electrical automation, in particular to a power flow calculation method and device of a power supply system, a storage medium and an electronic device.
Background
In the field of direct-current traction power supply, in order to reduce impact on a power grid and directly absorb regenerative braking energy of a locomotive on a contact network as much as possible, system-level cooperative control needs to be introduced, meanwhile, in order to increase the power supply range and capacity of each traction substation under the condition of meeting the requirement of power quality, a converter based on a full-control device can be adopted to replace the original uncontrolled diode rectification scheme, and a flexible direct-current traction power supply system can meet the requirement, and has the advantages of energy conservation, high power quality and the like. The load flow calculation is used as a steady-state analysis method of the power system, and the voltage loss, network loss, locomotive regenerative braking energy consumption and other information of the whole system can be calculated efficiently and conveniently. Therefore, it is necessary to apply load flow calculation to perform steady-state analysis and calculation on the flexible direct-current traction power supply system.
However, the current load flow calculation method of the direct current traction power supply system mainly aims at a mode of uncontrolled rectification of a diode or a mode with an inversion feedback device, and lacks calculation and analysis of the flexible direct current traction power supply system. In addition, the load flow calculation in the traction power supply system needs to know the locomotive load, the locomotive load depends on the voltage of a vehicle-end contact network, and the voltage of the contact network needs to be determined through the load flow calculation, so that the load flow calculation and the traction calculation in the traction power supply system have a coupling relation in nature. The original method firstly completes traction calculation, then utilizes locomotive position, power or current information obtained by the traction calculation to perform load flow calculation, and ignores the coupling relation between the locomotive traction calculation and the load flow calculation.
Aiming at the problem that the coupling relation between traction calculation and load flow calculation is neglected in the related technology, so that accurate analysis and calculation cannot be carried out on a power supply system, an effective solution is not provided at present.
Disclosure of Invention
The embodiment of the invention provides a power flow calculation method, a power flow calculation device, a storage medium and an electronic device of a power supply system, which are used for at least solving the problem that the power supply system cannot be accurately analyzed and calculated due to the fact that the coupling relation between traction calculation and power flow calculation cannot be considered in the related technology.
According to an embodiment of the present invention, there is provided a power flow calculation method for a power supply system, including: initializing the node voltage in the power supply system according to the target time; after the initialization is completed, executing a target operation, wherein the target operation comprises the following processing steps: repeatedly executing predetermined processing on the power supply system until the absolute value of the difference between a first voltage and a second voltage of a first target node included in the power supply system is smaller than a first maximum allowable error, wherein the predetermined processing comprises sequentially executing traction calculation, load flow calculation and judging whether the absolute value of the difference between the first voltage and the second voltage is smaller than the maximum allowable error; under the condition that the absolute value of the difference value of the first voltage and the second voltage is determined to be smaller than the first maximum allowable error, judging whether the target time reaches the end time or not; in a case where it is determined that the target time reaches the end time, a result of performing the predetermined processing is determined as a power flow calculation result.
In one exemplary embodiment, after determining whether the target time reaches the end time, the method further comprises: under the condition that the target time is determined not to reach the end time, updating the target time according to a preset updating rule, and re-determining the updated target time as the target time; and performing the initialization on the node voltage in the power supply system by using the target time, and executing the target operation after the initialization is completed.
In an exemplary embodiment, updating the target time according to a predetermined update rule includes: and updating the target time by adding the target time and a preset time step, wherein the value obtained after the addition is the updated target time.
In one exemplary embodiment, performing the power flow calculation on the power supply system includes: repeatedly executing the load flow calculation until the absolute value of the difference value between the third voltage of the second target node calculated in the current load flow calculation and the fourth voltage of the second target node calculated in the previous load flow calculation is smaller than a second maximum allowable error; in the repeated execution of the power flow calculation, the voltage value in each traction used in the next power flow calculation and the resistance value of the target resistor are updated according to a preset updating mode.
In an exemplary embodiment, the predetermined update mode comprises a system level control mode.
According to another embodiment of the present invention, there is also provided a power flow calculation apparatus of a power supply system, including: the initialization module is used for initializing the node voltage in the power supply system according to the target time; a processing module, configured to execute a target operation after the initialization is completed, where the target operation includes the following processing steps: repeatedly executing predetermined processing on the power supply system until the absolute value of the difference value between a first voltage and a second voltage of a first target node included in the power supply system is smaller than a first maximum allowable error, wherein the predetermined processing includes traction calculation, power flow calculation and judgment of whether the absolute value of the difference value between the first voltage and the second voltage is smaller than the maximum allowable error, which are sequentially executed; under the condition that the absolute value of the difference value of the first voltage and the second voltage is determined to be smaller than the first maximum allowable error, judging whether the target time reaches the end time or not; and determining a result obtained by executing the predetermined processing as a power flow calculation result in the case that the target time is determined to reach the end time.
In one exemplary embodiment, the processing module further comprises: the updating unit is used for updating the target time according to a preset updating rule and re-determining the updated target time as the target time under the condition that the target time is determined not to reach the end time; and the processing unit is used for carrying out the initialization on the node voltage in the power supply system by using the target time and executing the target operation after the initialization is completed.
In one exemplary embodiment, the update unit includes: and the updating subunit is used for updating the target time by adding the target time and a preset time step, wherein the value obtained after the addition is the updated target time.
According to a further embodiment of the present invention, there is also provided a computer-readable storage medium having a computer program stored thereon, wherein the computer program is arranged to perform the steps of any of the above method embodiments when executed.
According to yet another embodiment of the present invention, there is also provided an electronic device, comprising a memory in which a computer program is stored and a processor configured to run the computer program to perform the steps of any of the method embodiments described above.
According to the invention, the load flow calculation and the traction calculation are coupled with each other, and are subjected to iterative calculation at each moment, so that the steady-state analysis and calculation of the power supply system can be better performed. The problem that the coupling relation between traction calculation and load flow calculation cannot be considered in the related technology, so that accurate analysis and calculation cannot be carried out on a power supply system is solved, and the effect of more accurately analyzing and calculating the power supply system is achieved.
Drawings
Fig. 1 is a block diagram of a hardware configuration of a mobile terminal of a power flow calculation method of a power supply system according to an embodiment of the present invention;
fig. 2 is a flowchart of a power flow calculation method of a power supply system according to an embodiment of the present invention;
fig. 3 is a flowchart of a power flow calculation method of a preferred power supply system according to an embodiment of the present invention;
FIG. 4 is a flow chart of a traction calculation for a power supply system according to an embodiment of the present invention;
FIG. 5 is a schematic diagram of the internal circuitry of a preferred traction device according to an embodiment of the present invention;
FIG. 6 is a schematic diagram of a preferred power flow calculation node arrangement according to an embodiment of the invention;
FIG. 7 is a flow chart of a preferred power flow calculation according to an embodiment of the invention;
fig. 8 is a block diagram of a power flow calculation device of a power supply system according to an embodiment of the present invention.
Detailed Description
Embodiments of the present invention will be described in detail below with reference to the accompanying drawings in conjunction with embodiments.
It should be noted that the terms "first," "second," and the like in the description and claims of the present invention and in the drawings described above are used for distinguishing between similar elements and not necessarily for describing a particular sequential or chronological order.
The method embodiments provided in the embodiments of the present application may be executed in a mobile terminal, a computer terminal, or a similar computing device. Taking an example of the operation on a mobile terminal, fig. 1 is a hardware structure block diagram of the mobile terminal of a power flow calculation method of a power supply system according to an embodiment of the present invention. As shown in fig. 1, the mobile terminal may include one or more processors 102 (only one is shown in fig. 1) (the processor 102 may include, but is not limited to, a processing device such as a microprocessor MCU or a programmable logic device FPGA, etc.) and a memory 104 for storing data, wherein the mobile terminal may further include a transmission device 106 for communication functions and an input-output device 108. It will be understood by those of ordinary skill in the art that the structure shown in fig. 1 is only an illustration and is not intended to limit the structure of the mobile terminal. For example, the mobile terminal may also include more or fewer components than shown in FIG. 1, or have a different configuration than shown in FIG. 1.
The memory 104 may be used to store a computer program, for example, a software program and a module of application software, such as a computer program corresponding to the power flow calculation method of the power supply system in the embodiment of the present invention, and the processor 102 executes various functional applications and data processing by running the computer program stored in the memory 104, so as to implement the above method. The memory 104 may include high speed random access memory, and may also include non-volatile memory, such as one or more magnetic storage devices, flash memory, or other non-volatile solid-state memory. In some examples, the memory 104 may further include memory located remotely from the processor 102, which may be connected to the mobile terminal over a network. Examples of such networks include, but are not limited to, the internet, intranets, local area networks, mobile communication networks, and combinations thereof.
The transmission device 106 is used to receive or transmit data via a network. Specific examples of the network described above may include a wireless network provided by a communication provider of the mobile terminal. In one example, the transmission device 106 includes a Network adapter (NIC), which can be connected to other Network devices through a base station so as to communicate with the internet. In one example, the transmission device 106 may be a Radio Frequency (RF) module, which is used to communicate with the internet in a wireless manner.
In the present embodiment, a power flow calculation method of a power supply system is provided, and fig. 2 is a flowchart of the power flow calculation method of the power supply system according to the embodiment of the present invention, as shown in fig. 2, the flowchart includes the following steps:
step S202, initializing the node voltage in the power supply system according to target time;
step S204, after the initialization is completed, executing target operation, wherein the target operation comprises the following processing steps:
s2042, repeatedly executing predetermined processing on the power supply system until an absolute value of a difference between a first voltage and a second voltage of a first target node included in the power supply system is smaller than a first maximum allowable error, wherein the predetermined processing includes sequentially executed traction calculation, power flow calculation, and determination of whether the absolute value of the difference between the first voltage and the second voltage is smaller than the maximum allowable error;
s2044, under the condition that the absolute value of the difference value between the first voltage and the second voltage is smaller than the first maximum allowable error, judging whether the target time reaches the end time or not;
s2046, in the case where it is determined that the target time reaches the end time, determining a result of performing the predetermined processing as a power flow calculation result.
Through the steps, the load flow calculation and the traction calculation are coupled with each other, and are subjected to iterative calculation at each moment, so that the steady-state analysis and calculation of the power supply system can be better performed. The problem that the coupling relation between traction calculation and load flow calculation cannot be considered in the related technology, so that accurate analysis and calculation cannot be carried out on a power supply system is solved, and the effect of more accurately analyzing and calculating the power supply system is achieved.
The main body of the above steps may be a computing device in the traction power supply system, or a computer terminal device, or a processor with human-computer interaction capability configured on a storage device, or a processing device or a processing unit with similar processing capability, but is not limited thereto. The following description will be given by taking the operation performed by the computing device in the traction power supply system as an example (which is only an exemplary description, and other devices or modules may also perform the operation in actual operation):
in the above embodiment, the computing device may initialize the node voltage of the power supply system at a certain time (e.g. t), and after the initialization is completed, sequentially perform the traction calculation and the power flow calculation, and the traction calculation and the power flow calculation are both performed cyclically (for example, the cycle is an outer cycle), and the condition of the end of the cycle is that the absolute value of the difference of the voltages of a certain node (e.g. the node Train of the locomotive) in two consecutive cycles (e.g. the j (th) and the j (1) (th) cycles) in the power supply system is smaller than a specified maximum allowable error (e.g. the error is epsilon) t ) I.e. | U t Train,j -U t Train,j-1 |<ε t The loop calculation is ended, and after the loop calculation is ended, a preset end time (e.g., t) is reached at the time t end ) And obtaining a load flow calculation result. The method provided by the embodiment fully considers the coupling relationship between the traction calculation and the load flow calculation, and performs iterative calculation on the traction calculation and the load flow calculation at each moment, so that the steady-state analysis and calculation can be better performed on the power supply system, and the problem that the coupling relationship between the traction calculation and the load flow calculation cannot be considered in the related technology is solved.
In an optional embodiment, after determining whether the target time reaches the end time, the method further comprises: under the condition that the target time is determined not to reach the end time, updating the target time according to a preset updating rule, and re-determining the updated target time as the target time; performing the initialization on the node voltage in the power supply system using the target time, and performing the target operation after the initialization is completed. In the present embodiment, in the case where it is determined that the target time does not reach the end time (i.e., t < t) end In this case), the target time is updated according to a predetermined update rule, and the updated target time is determined as the target time again, for example, x seconds is added on the basis of the target time to obtain a new target time, optionally, in practical applications, y seconds is decreased on the basis of the target time to obtain a new target time, after the new target time is obtained, the node voltage in the power supply system is initialized, and the operation corresponding to the foregoing step S204 is performed.
In an optional embodiment, updating the target time according to a predetermined update rule includes: and updating the target time by adding the target time and a preset time step, wherein the value obtained after the addition is the updated target time. In this embodiment, the target time may be updated by adding a certain time step Δ t to the target, in practical applications, Δ t may be set as needed, and the updated target time is t + Δ t.
In an optional embodiment, performing the power flow calculation on the power supply system includes: repeatedly executing the power flow calculation until the absolute value of the difference between the third voltage of the second target node calculated in the current power flow calculation and the fourth voltage of the second target node calculated in the previous power flow calculation is smaller than a second maximum allowable error; in the process of repeatedly executing the power flow calculation, the voltage value used in each traction in the next power flow calculation and the resistance value of the target resistor are updated according to a preset updating mode. In this embodiment, based on the result of the aforementioned traction calculation (e.g. the jth outer-loop traction calculation), a loop (for example, the j-th outer-loop traction calculation is distinguished, and the loop is an inner loop) to perform the power flow calculation, and the condition that the inner loop ends is that the absolute value of the difference between the voltage value of a certain node in the power supply system in two consecutive loop calculations (e.g. the k-th and k-1-th loops) is smaller than a prescribed maximum allowable error (e.g. the error is epsilon) p ) I.e. U t k,j -U t k-1,j |<ε p The inner loop calculation is ended, and the voltage value and the resistance value related to the traction and the resistance value of the target resistor in the inner loop calculation process can be different in the k-th and k-1-th inner loop calculations, namely, the voltage value and the resistance value calculated in the k-th inner loop calculation can be updated on the basis of the voltage value and the resistance value calculated in the k-1-th inner loop calculation.
In an optional embodiment, the predetermined update mode comprises a system level control mode. In this embodiment, the voltage value and the resistance value of the target resistor related to each traction can be updated in a system-level control manner in different power flow inner loop calculations.
It is to be understood that the above-described embodiments are only a few, and not all, embodiments of the present invention.
The present invention will be described in detail with reference to the following examples:
fig. 3 is a flowchart of a power flow calculation method of a preferred power supply system according to an embodiment of the invention, and as shown in fig. 3, the flowchart includes the following steps:
s302, initialization of basic parameters, such as t =1,j =1,k =1;
s304, inputting information such as system parameters, network structures, locomotive characteristics and the like;
s306, initializing the node voltage in the system;
s308, carrying out traction calculation;
s310, carrying out load flow calculation;
s312, judging the first type, namely judging the | U t Train,j -U t Train,j-1 |<ε t Whether the condition is satisfied, namely when the time is t, the node voltage vector U of the locomotive after the jth iteration t Train,j Node voltage vector U of locomotive after j-1 iteration t Train,j-1 Whether the absolute value of the difference is less than the maximum allowable error epsilon of the outer loop t
S314, under the condition that the first type judgment condition is not satisfied, adding 1 to the number of outer loops, namely preparing to execute the (j + 1) th iterative computation, and setting k to be 1;
s316, under the condition that the first type judgment condition is satisfied, the second type judgment is carried out again, and whether the time t is equal to the preset end time t or not is judged end
S318, if the second type of determination condition is not satisfied, returning the parameters j and k to the initial value 1, increasing a step Δ t for the time t to obtain a new target time t + Δ t, and repeating the operations corresponding to the steps S306 to S316;
s320, when the second type of determination condition is satisfied, that is, the time t = t end And when the flow is finished, the load flow calculation result is obtained.
In the above embodiment, the traction calculation may be described in conjunction with a traction calculation flowchart of the power supply system according to the embodiment of the present invention, as shown in fig. 4, the traction calculation flowchart includes the following steps:
s402, initializing the weight of the train, such as the position χ, the speed v, the acceleration a, the voltage U of a contact network at the train end, the current of the train and the running state;
s404, determining which running state the train is in (such as traction/inertia/braking/stopping);
s406, determining I according to U, the running state and v according to the traction characteristics of the electric locomotive;
s408, calculating a according to Newton' S second law by using the v, I and the resistance information;
s410, calculating a new velocity value v from v = v + α x Δ t;
s412, according to χ = χ + v · Δ t +0.5 · α · Δ t 2 Calculating a new position x;
and S414, finishing the traction calculation and outputting train operation information.
In the above embodiment, the locomotive is modeled as a current source, and at each time, the magnitude I of the locomotive current needs to be updated according to the speed v, the end contact network voltage U and the operating state.
Alternatively, the method can be used in a flexible dc traction power supply system, and droop control is generally adopted in the research of the flexible dc traction power supply system, or a method for improving droop control is adopted. The basic form of droop control is:
U ref =U 0 -k 0 ·I (1)
in the formula of U ref Representing the converter output voltage reference, U 0 Intercept, k, representing sag characteristic 0 The slope of the droop characteristic is shown, I represents the output current of the converter, namely the output current of a traction substation in the flexible direct current traction power supply system. Because the inner ring controller of the converter has higher control speed, the output voltage of the converter can quickly track the reference value, and in the steady-state analysis, the formula (1) can be simplified into the formula:
U=U 0 -k 0 ·I (2)
in the formula, U represents the output voltage of the converter, namely the output node voltage of the traction substation in the flexible direct current traction power supply system.
FIG. 5 is a schematic diagram of the internal circuitry of a preferred traction device according to an embodiment of the present invention, as shown in FIG. 5, wherein R k The size of (2) and k 0 Equal, U 0 Is equal to U in the formula (2) 0 Are equal.
In order to more clearly illustrate the above embodiments, the preferred power flow calculation node setting diagram according to the embodiments of the present invention is specifically described, a node setting diagram corresponding to fig. 5 is shown in fig. 6, three types of nodes, namely TSSU, TSSI, and Train, are shared in fig. 6, and the nth locomotive is represented by a Train n node. The mth TSS is represented by two nodes of TSSU m and TSSI m in power flow calculation, wherein the node voltage of TSSU m is a known quantity U 0 The node injection current of TSSI m is 0. The resistance between two nodes of TSSU m and TSSI m is R k
In an alternative embodiment, the above method is not only applicable to droop control situations, but also by dynamically adjusting U 0 And R k Can also be applied to various system-level control situations, as shown in fig. 7, a preferred power flow calculation flowchart according to an embodiment of the present invention is shown, and the flowchart includes the following steps:
s702, the method comprises the steps of calculating a vehicle position, a system parameter, a network structure,R of each traction station k Calculating a node admittance matrix Y;
s704, updating I according to the traction calculation result Train Is shown by TSSI Setting zero and U TSSU
S706, calculating U according to YU = I Train,k And U TSSI,k
S708, calculating U TSSI,k =U TSSI,k ·α+U TSSI,k-1 ·(1-α),U Train,k =U Train,k ·α+U Train,k-1 ·(1-α);
S710, calculating I according to YU = I TSSU,k
S712, judge | U t k,j -U t k-1,j |<ε p Whether the condition is satisfied;
s714, if the above-mentioned determination condition is not satisfied, adding 1 to the number k of internal cycles;
s716, the R of each traction station can be updated in a system-level control mode k And U TSSU
After the above steps are performed, step S312 and the subsequent operations in the foregoing embodiments may also be performed.
And alpha is a damping factor, and under the condition that the tide is not easy to converge, the algorithm convergence can be assisted by adjusting alpha. U shape Train,k 、U TSSI,k And I TSSU,k Respectively representing the node voltage vector of the Train node, the node voltage vector of the TSSI and the injection current of the TSSU node when the iteration number of the inner ring is k p The maximum allowable error of the inner ring.
Through the description of the foregoing embodiments, it is clear to those skilled in the art that the method according to the foregoing embodiments may be implemented by software plus a necessary general hardware platform, and certainly may also be implemented by hardware, but the former is a better implementation mode in many cases. Based on such understanding, the technical solutions of the present invention may be embodied in the form of a software product, which is stored in a storage medium (e.g., ROM/RAM, magnetic disk, optical disk) and includes instructions for enabling a terminal device (e.g., a mobile phone, a computer, a server, or a network device) to execute the method according to the embodiments of the present invention.
In this embodiment, there is further provided a power flow calculation apparatus of a power supply system, and fig. 8 is a block diagram of a structure of the power flow calculation apparatus of the power supply system according to the embodiment of the present invention, as shown in fig. 8, the apparatus includes:
an initialization module 802, configured to initialize a node voltage in a power supply system according to a target time;
a processing module 804, configured to execute a target operation after the initialization is completed, where the target operation includes the following processing steps: repeatedly executing predetermined processing on the power supply system until the absolute value of the difference value between a first voltage and a second voltage of a first target node included in the power supply system is smaller than a first maximum allowable error, wherein the predetermined processing includes traction calculation, power flow calculation and judgment of whether the absolute value of the difference value between the first voltage and the second voltage is smaller than the maximum allowable error, which are sequentially executed; under the condition that the absolute value of the difference value of the first voltage and the second voltage is determined to be smaller than the first maximum allowable error, judging whether the target time reaches the end time or not; in a case where it is determined that the target time reaches the end time, a result of performing the predetermined processing is determined as a power flow calculation result.
In an alternative embodiment, the processing module 804 includes: the updating unit is used for updating the target time according to a preset updating rule under the condition that the target time is determined not to reach the end time, and re-determining the updated target time as the target time; and the processing unit is used for carrying out the initialization on the node voltage in the power supply system by using the target time and executing the target operation after the initialization is completed.
In an optional embodiment, the updating unit includes: and the updating subunit is used for updating the target time by adding the target time and a preset time step, wherein the value obtained after the addition is the updated target time.
In an alternative embodiment, the processing module 804 includes: a power flow calculation unit for performing the power flow calculation on the power supply system by: repeatedly executing the power flow calculation until the absolute value of the difference between the third voltage of the second target node calculated in the current power flow calculation and the fourth voltage of the second target node calculated in the previous power flow calculation is smaller than a second maximum allowable error; in the repeated execution of the power flow calculation, the voltage value in each traction used in the next power flow calculation and the resistance value of the target resistor are updated according to a preset updating mode.
In an optional embodiment, the predetermined update mode comprises a system level control mode.
It should be noted that the above modules may be implemented by software or hardware, and for the latter, the following may be implemented, but not limited to: the modules are all positioned in the same processor; alternatively, the modules are located in different processors in any combination.
An embodiment of the present invention further provides a computer-readable storage medium, in which a computer program is stored, wherein the computer program is configured to perform the steps in any of the above method embodiments when executed.
In an exemplary embodiment, the computer-readable storage medium may include, but is not limited to: various media capable of storing computer programs, such as a usb disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a removable hard disk, a magnetic disk, or an optical disk.
Embodiments of the present invention further provide an electronic device, comprising a memory in which a computer program is stored and a processor configured to execute the computer program to perform the steps in any of the above method embodiments.
In an exemplary embodiment, the electronic apparatus may further include a transmission device and an input/output device, wherein the transmission device is connected to the processor, and the input/output device is connected to the processor.
For specific examples in this embodiment, reference may be made to the examples described in the above embodiments and exemplary embodiments, and details of this embodiment are not repeated herein.
According to the power supply system load flow calculation method, the load flow calculation and the traction calculation are coupled with each other, and are subjected to iterative calculation at each moment, so that the steady-state analysis and calculation can be better performed on the power supply system. The problem that the coupling relation between traction calculation and load flow calculation cannot be considered in the related technology is solved, and the effect of more accurately analyzing and calculating the power supply system is achieved.
It will be apparent to those skilled in the art that the modules or steps of the present invention described above may be implemented in a general purpose computing device, they may be centralized in a single computing device or distributed across a network of multiple computing devices, and they may be implemented in program code that is executable by a computing device, such that they may be stored in a memory device and executed by a computing device, and in some cases, the steps shown or described may be executed in an order different from that shown or described herein, or they may be separately fabricated into individual integrated circuit modules, or multiple modules or steps therein may be fabricated into a single integrated circuit module. Thus, the present invention is not limited to any specific combination of hardware and software.
The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the principle of the present invention should be included in the protection scope of the present invention.

Claims (10)

1. A power flow calculation method of a power supply system is characterized by comprising the following steps:
initializing the node voltage in the power supply system according to the target time;
after the initialization is completed, executing a target operation, wherein the target operation comprises the following processing steps:
repeatedly executing predetermined processing on the power supply system until the absolute value of the difference value between a first voltage and a second voltage of a first target node included in the power supply system is smaller than a first maximum allowable error, wherein the predetermined processing includes traction calculation, power flow calculation and judgment of whether the absolute value of the difference value between the first voltage and the second voltage is smaller than the maximum allowable error, which are sequentially executed;
under the condition that the absolute value of the difference value of the first voltage and the second voltage is determined to be smaller than the first maximum allowable error, judging whether the target time reaches the end time or not;
determining a result obtained by executing the predetermined processing as a power flow calculation result in the case where it is determined that the target time reaches the end time;
wherein the content of the first and second substances,
the traction calculation includes the steps of:
initializing parameters of the train to obtain the following parameters: the method comprises the following steps of (1) initial weight, initial position, initial speed, initial acceleration, initial train end catenary voltage, initial train current and initial running state; determining a target running state of the train at the target time; determining a target current of the train at the target time according to the initial train end contact network voltage, the target running state and the initial speed; calculating a target acceleration based on the initial speed, the target current and the resistance information; calculating a new speed according to the initial speed, the target acceleration and a preset time step length; calculating a new position according to the initial position, the new speed, the preset time step length and the target acceleration; outputting a traction calculation result of the train, wherein the traction calculation result comprises the target current and the new position of the train;
the power flow calculation comprises the following steps:
based on the new position, system parameters, network structure, resistance R of each traction station k Calculating a node admittance matrix Y; updating the current I of the train based on the traction calculation result Train Injecting a current I into the node of the second target node TSSI Set to zero and set the node voltage U of the third target node TSSU (ii) a Calculating a voltage vector U of the first target node according to YU = I Train,k And a voltage vector U of the second target node TSSI,k (ii) a Calculate U TSSI,k =U TSSI,k ·α+U TSSI,k-1 ·(1-α),U Train,k =U Train,k ·α+U Train,k-1 1- α; calculating an injection current vector I of the third target node according to YU = I TSSU,k (ii) a Determine | U t k,j -U t k-1,j |<ε p Whether the condition is satisfied; under the condition that the judgment condition is not satisfied, adding 1 to the internal cycle times k; wherein R of each of the tractors k And U TSSU Can be updated by a system-level control mode, alpha is a damping factor, and alpha is allowed to be adjusted, U Train,k 、U TSSI,k And I TSSU,k Respectively representing the node voltage vector of the first target node, the voltage vector of the second target node and the injection current of the third target node when the iteration number of the inner ring is k p The maximum allowable error of the inner ring.
2. The method of claim 1, wherein after determining whether the target time reaches an end time, the method further comprises:
under the condition that the target time is determined not to reach the end time, updating the target time according to a preset updating rule, and re-determining the updated target time as the target time;
and performing the initialization on the node voltage in the power supply system by using the target time, and executing the target operation after the initialization is completed.
3. The method of claim 2, wherein updating the target time according to a predetermined update rule comprises:
and updating the target time by adding the target time and a preset time step, wherein the value obtained after the addition is the updated target time.
4. The method of claim 1, wherein performing the power flow calculation on the power supply system comprises:
repeatedly executing the power flow calculation until the absolute value of the difference between the third voltage of the second target node calculated in the current power flow calculation and the fourth voltage of the second target node calculated in the previous power flow calculation is smaller than a second maximum allowable error;
in the repeated execution of the power flow calculation, the voltage value in each traction used in the next power flow calculation and the resistance value of the target resistor are updated according to a preset updating mode.
5. The method of claim 4, wherein the predetermined update mode comprises a system level control mode.
6. A power flow calculation apparatus for a power supply system, comprising:
the initialization module is used for initializing the node voltage in the power supply system according to the target time;
a processing module, configured to execute a target operation after the initialization is completed, where the target operation includes the following processing steps:
repeatedly executing predetermined processing on the power supply system until the absolute value of the difference value between a first voltage and a second voltage of a first target node included in the power supply system is smaller than a first maximum allowable error, wherein the predetermined processing includes traction calculation, power flow calculation and judgment of whether the absolute value of the difference value between the first voltage and the second voltage is smaller than the maximum allowable error, which are sequentially executed;
under the condition that the absolute value of the difference value of the first voltage and the second voltage is determined to be smaller than the first maximum allowable error, judging whether the target time reaches the end time or not;
determining a result obtained by executing the predetermined processing as a power flow calculation result in the case where it is determined that the target time reaches the end time;
wherein the processing module may perform the traction calculation as follows:
initializing parameters of the train to obtain the following parameters: the method comprises the following steps of (1) obtaining initial weight, initial position, initial speed, initial acceleration, initial train end catenary voltage, initial train current and initial running state;
determining a target running state of the train at the target time;
determining a target current of the train at the target time according to the initial train end contact network voltage, the target running state and the initial speed;
calculating a target acceleration based on the initial speed, the target current and the resistance information;
calculating a new speed according to the initial speed, the target acceleration and a preset time step length;
calculating a new position according to the initial position, the new speed, the preset time step length and the target acceleration;
outputting a traction calculation result of the train, wherein the traction calculation result comprises the target current and the new position of the train;
the processing module may also perform the load flow calculation as follows:
based on the new position, system parameters, network structure, resistance R of each traction station k Calculating a node admittance matrix Y;
updating the current I of the train based on the traction calculation result Train Injecting a current I into the node of the second target node TSSI Set to zero and set the node voltage U of the third target node TSSU
Calculating a voltage vector U of the first target node according to YU = I Train,k And a voltage vector U of the second target node TSSI,k
Calculate U TSSI,k =U TSSI,k ·α+U TSSI,k-1 ·(1-α),U Train,k =U Train,k ·α+U Train,k-1 ·(1-α);
Calculating an injection current vector I of the third target node according to YU = I TSSU,k
Determine | U t k,j -U t k-1,j |<ε p Whether the condition is satisfied;
under the condition that the judgment condition is not satisfied, adding 1 to the internal cycle times k;
wherein R of each of the tractors k And U TSSU Can be updated by a system-level control mode, alpha is a damping factor, and alpha is allowed to be adjusted, U Train,k 、U TSSI,k And I TSSU,k Respectively representing the node voltage vector of the first target node, the voltage vector of the second target node and the injection current of the third target node when the iteration number of the inner ring is k p The maximum allowable error of the inner ring.
7. The apparatus of claim 6, wherein the processing module further comprises:
the updating unit is used for updating the target time according to a preset updating rule under the condition that the target time is determined not to reach the end time, and re-determining the updated target time as the target time;
and the processing unit is used for carrying out the initialization on the node voltage in the power supply system by using the target time and executing the target operation after the initialization is finished.
8. The apparatus of claim 7, wherein the updating unit comprises:
and the updating subunit is used for updating the target time by adding the target time and a preset time step, wherein a value obtained after the addition is the updated target time.
9. A computer-readable storage medium, characterized in that a computer program is stored in the computer-readable storage medium, wherein the computer program, when being executed by a processor, carries out the steps of the method as claimed in any one of the claims 1 to 5.
10. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the steps of the method as claimed in any of claims 1 to 5 are implemented when the computer program is executed by the processor.
CN202011425402.8A 2020-12-08 2020-12-08 Power flow calculation method and device of power supply system, storage medium and electronic device Active CN112510716B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011425402.8A CN112510716B (en) 2020-12-08 2020-12-08 Power flow calculation method and device of power supply system, storage medium and electronic device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011425402.8A CN112510716B (en) 2020-12-08 2020-12-08 Power flow calculation method and device of power supply system, storage medium and electronic device

Publications (2)

Publication Number Publication Date
CN112510716A CN112510716A (en) 2021-03-16
CN112510716B true CN112510716B (en) 2022-11-04

Family

ID=74971573

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011425402.8A Active CN112510716B (en) 2020-12-08 2020-12-08 Power flow calculation method and device of power supply system, storage medium and electronic device

Country Status (1)

Country Link
CN (1) CN112510716B (en)

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111725812B (en) * 2020-06-04 2023-04-11 海仿(上海)科技有限公司 Load flow calculation method, device, equipment and storage medium of large-scale power distribution system

Also Published As

Publication number Publication date
CN112510716A (en) 2021-03-16

Similar Documents

Publication Publication Date Title
US20210133536A1 (en) Load prediction method and apparatus based on neural network
CN111638646B (en) Training method and device for walking controller of quadruped robot, terminal and storage medium
CN109541961A (en) For the discrete state event-driven simulation method of power electronics hybrid system emulation
CN112818588B (en) Optimal power flow calculation method, device and storage medium of power system
CN109710404B (en) Task scheduling method in distributed system
CN113113928B (en) Flexible-direct system direct-current bus voltage control method and device based on deep reinforcement learning
CN103885026A (en) Electric energy meter calibrating method
CN107742885A (en) Power distribution network voltage power sensitivity estimation method based on regular matching pursuit
CN109638883B (en) Power grid voltage control method and device, computer equipment and storage medium
CN103810646A (en) Improved projection integral algorithm based active power distribution system dynamic simulation method
JP6894724B2 (en) Power equipment planning method, power equipment manufacturing method, power equipment planning equipment and programs
Li et al. Learning-based predictive control via real-time aggregate flexibility
CN112510716B (en) Power flow calculation method and device of power supply system, storage medium and electronic device
CN104090496A (en) Smart grid control operation continuous analog simulation method
CN105404529A (en) Popup method, device and system of operation system popup for reinstallation
CN106054665B (en) A kind of large-scale photovoltaic inverter system divides group&#39;s equivalent modeling method
CN113033027B (en) Gait control method and device of intelligent robot, storage medium and electronic device
CN106789289B (en) Method and apparatus for virtual network mapping
CN115562299A (en) Navigation method and device of mobile robot, mobile robot and medium
CN113992520A (en) Virtual network resource deployment method and system
Puin et al. Performance assessment of an optimization strategy proposed for power systems
CN109861304A (en) A kind of micro-capacitance sensor economic load dispatching method considering communication Time-varying time-delays
CN107729644A (en) It is a kind of that the method for improving available time and realizing minimum dominant set is calculated using film
CN115129093B (en) Temperature control method, temperature control device and storage medium for power assembly
CN110059326A (en) Simulation method and simulation system for debugging microwave device

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
GR01 Patent grant
GR01 Patent grant