Power system load flow calculation method based on recursive decomposition and related device
By transforming the current calculation equation of the power system into a state variable expression form that combines constants and nonlinear functions, and processing it based on recursive decomposition, the existing methods solve the problem of computational failure and error when solving high-dimensional nonlinear equation systems, and efficient and accurate current calculation is achieved.
Patent Information
- Application Number
- CN202510477229.2
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2025-04-16
- Publication Date
- 2025-08-01
AI Technical Summary
The current power system current calculation methods are prone to failure and have large errors when solving high-dimensional nonlinear equation systems. Especially when Newton-Ravson's method is sensitive to initial values and cannot be calculated when the Jacobian matrix is singular, the linearization method leads to large calculation errors.
The recursive decomposition method is used to transform the power system flow calculation equation into a state variable expression form that combines constant and nonlinear functions, and the estimation method of the constant part and the iterative calculation formula of the nonlinear function part are determined through the recursive decomposition method, and the iterative calculation formula of the nonlinear function part is iteratively calculated until the termination condition is reached.
Without computing the Jacobian matrix, it is possible to perform calculations when the Jacobian matrix is singular, reducing computational complexity and improving computational accuracy, and is suitable for different types of power systems and trend computing scenarios.
Smart Images
Figure CN120414554A_ABST
Abstract
Description
Technical Field
[0001] This application belongs to a power flow calculation method, and specifically relates to a power flow calculation method and related devices for a power system based on recursive decomposition. Background Art
[0002] In power flow calculation of a power system, it is necessary to solve a high-dimensional non-linear equation set, and the calculation complexity is relatively high, usually difficult to solve directly. In conventional power flow calculation of a power system, the Newton-Raphson method is mainly used for calculation. However, this method requires calculating the Jacobian matrix. When the Jacobian matrix is singular, this method will completely fail. Moreover, this method is sensitive to the initial value. When the initial value is poorly selected, the convergence speed of the power flow calculation is difficult to meet the expectation. In addition to the Newton-Raphson method, some researchers have also proposed linearization methods. However, when dealing with non-linear problems through linearization methods, large calculation errors will be caused due to the lack of non-linear characteristics. Summary of the Invention
[0003] This application aims at the technical problems of high-dimensional non-linear equation set solving methods in current power flow calculation of a power system, such as easy calculation failure and large errors, and provides a power flow calculation method and related devices for a power system based on recursive decomposition.
[0004] To achieve the above object, this application is implemented by adopting the following technical solutions: In the first aspect, this application proposes a power flow calculation method for a power system based on recursive decomposition, including: Obtain the composition information of the power system network model; Combined with the composition information of the power system network model, transform the power flow calculation equation of the power system to be solved into a state variable expression form composed of constants and non-linear functions; For the transformed power flow calculation equation of the power system, based on the recursive decomposition method, determine the estimation method of the constant part and the iterative calculation formula of the non-linear function part; Perform iterative calculation on the non-linear function part. When the iterative calculation termination condition is reached, integrate the estimated value of the constant part and the iterative updated value of the non-linear function part to obtain the calculation result of the state variable of the power flow calculation equation of the power system.
[0005] Furthermore, the composition information of the power system network model includes at least one of node data, node type, and admittance matrix information.
[0006] Furthermore, the transformation into the state variable expression form composed of constants and non-linear functions includes: For the i th load node, the state variable expression form is:
[0007]
[0008] Among them, n is the number of nodes in the power system network model, is the voltage of node i ; is the voltage of node k ; is the reactive power of node i ; is the active power of node i ; is the magnitude of the self-admittance of node i ; is the phase of the self-admittance of node i ; is the phase of the mutual admittance between node i and node k ; is the magnitude of the mutual admittance between node i and node k ; is the voltage phase of node i ; is the voltage phase of node k ; For the i th generating node, the expression form of the state variable is: .
[0009] Furthermore, the estimation method of the constant part includes: For the i th load node, the estimation formula of the constant part is:
[0010]
[0011] Among them, is the constant part of , that is, the initial value of ; is the constant part of , that is, the initial value of ; is the active power of node i ; For the i th generating node, the estimation formula of the constant part is:
[0012] Among them, is the nodei The magnitude of the self - admittance, is the node i The phase of the self - admittance.
[0013] Furthermore, the iterative calculation formula of the non - linear function part includes: For the i th load node, the iterative calculation formula of the non - linear function part is:
[0014]
[0015] Wherein, is the voltage phase of the node at the (m + 1)-th iteration i , is the node i and the node k The magnitude of the mutual - admittance between them, is the node i and the node k The phase of the mutual - admittance between them, is the voltage of the node at the (m + 1)-th iteration i , is the node i Reactive power of the node, is the node i The magnitude of the self - admittance, is the node i The phase of the self - admittance, is the node i Voltage of the node, is the th iteration, is the th iteration node i Voltage of the node, is the th iteration node k Voltage of the node, is the th iteration node i Voltage phase of the node, is the th iteration node k Voltage phase of the node; For the i th generation node, the iterative calculation formula of the non - linear function part is: .
[0016] Furthermore, the calculation results of the state variables of the power - flow calculation equation of the power system include: For the i th load node, the calculation results of the state variables include:
[0017] Among them, is the voltage phase of node i . is the k -th iteration voltage phase of node i . is the voltage of node i , is the k -th iteration voltage amplitude of node i ; For the i -th power generation node, the calculation result of the state variable includes: .
[0018] On the second hand, the present application proposes a power flow calculation system based on recursive decomposition, including: An information acquisition module, configured to acquire the composition information of the power system network model; An equation transformation module, configured to combine the composition information of the power system network model to transform the power flow calculation equation to be solved into a state variable expression form combined with constants and non-linear functions; A decomposition module, configured to determine the estimation method of the constant part and the iterative calculation formula of the non-linear function part for the transformed power flow calculation equation based on the recursive decomposition method; A calculation module, configured to perform iterative calculation on the non-linear function part, and when the iterative calculation termination condition is reached, integrate the estimated value of the constant part and the iterative updated value of the non-linear function part to obtain the calculation result of the state variable of the power flow calculation equation.
[0019] Further, in the equation transformation module, transforming into a state variable expression form combined with constants and non-linear functions includes: For the i -th load node, the state variable expression form is:
[0020]
[0021] Among them, n is the number of nodes in the power system network model, is the voltage of node i , is the voltage of node k , is the reactive power of node i , is the i The active power is the node i the magnitude of the self - admittance is the node i the phase of the self - admittance is the node i and the node k the phase of the mutual - admittance between them is the node i and the node k the magnitude of the mutual - admittance between them is the voltage phase of the node i is the voltage phase of the node is the node k ; For the i th generation node, the expression form of the state variable is: .
[0022] Furthermore, in the decomposition module, the estimation method of the constant part includes: For the i th load node, the estimation formula of the constant part is:
[0023]
[0024] Wherein, is the constant part of, that is, the initial value of , is the constant part of, that is, the initial value of , is the active power of the node i ; For the i th generation node, the estimation formula of the constant part is:
[0025] Wherein, is the magnitude of the self - admittance of the node i , is the node i the phase of the self - admittance.
[0026] Furthermore, in the decomposition module, the iterative calculation formula of the non - linear function part includes: For the i th load node, the iterative calculation formula of the non - linear function part is:
[0027]
[0028] Among them, is the voltage phase of the (m + 1)-th iteration node i ; is the node i and the node k the magnitude of the mutual admittance between them; is the node i and the node k the phase of the mutual admittance between them; is the voltage of the (m + 1)-th iteration node i ; is the node i the reactive power of it; is the node i the magnitude of the self-admittance; is the node i the phase of the self-admittance; is the node i the voltage of it; is the -th iteration; is the -th iteration node i the voltage of it; is the -th iteration node k the voltage of it; is the -th iteration node i the voltage phase of it; is the -th iteration node k the voltage phase of it; For the i -th power generation node, the iterative calculation formula for the non-linear function part is: .
[0029] Furthermore, in the said calculation module, the calculation results of the state variables of the power flow calculation equation of the power system include: For the i -th load node, the calculation results of the state variables include:
[0030] Among them, is the voltage phase of the node i ; is the k -th iteration node i the voltage phase of it; is the node i the voltage of it; is the k -th iteration node iThe voltage amplitude; For the i th power generation node, the calculation results of the state variables include: .
[0031] In a third aspect, the present application proposes an electronic device, including: A memory for storing a computer program; A processor for implementing the steps of the above-mentioned power flow calculation method based on recursive decomposition when executing the computer program.
[0032] In a third aspect, the present application proposes a computer-readable storage medium, in which a computer program is stored, and when the computer program is executed by a processor, the steps of the above-mentioned power flow calculation method based on recursive decomposition are implemented.
[0033] Compared with the prior art, the present application has the following beneficial effects: The present application proposes a power flow calculation method based on recursive decomposition, transforms the power flow calculation equation to be solved, and respectively determines the estimation method of the constant part and the iterative calculation formula of the non-linear function part based on the recursive decomposition method. After iteration, the estimated value of the constant part is integrated to obtain the calculation result of the state variable. The power flow calculation method proposed by the present application, compared with the prior art, does not need to calculate the Jacobian matrix. When the Newton-Raphson method cannot perform power flow calculation due to the singularity of the Jacobian matrix, the method of the present application can still perform the calculation. In addition, the present application performs power flow calculation through iteration, and the calculation complexity is relatively low.
[0034] The present application also proposes a power flow calculation system based on recursive decomposition, an electronic device and a computer storage medium, which have all the advantages of the above-mentioned power flow calculation method. Drawings
[0035] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the following will briefly introduce the drawings required in the embodiments. It should be understood that the following drawings only show some embodiments of the present application, and therefore should not be regarded as limiting the scope. For those of ordinary skill in the art, other related drawings can be obtained based on these drawings without creative efforts.
[0036] Figure 1 It is the first flow schematic diagram of the power flow calculation method based on recursive decomposition of the present application; Figure 2 It is the second flow schematic diagram of the power flow calculation method based on recursive decomposition of the present application; Figure 3It is the system structure and line parameter diagram of a 3-node system in the embodiments of this application; Figure 4 It is a schematic diagram of a power system power flow calculation system based on recursive decomposition in this application. Specific implementation manners
[0037] To make the objectives, technical solutions and advantages of the embodiments of this application clearer, the technical solutions in the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings in the embodiments of this application. Obviously, the described embodiments are some but not all of the embodiments of this application. Generally, the components of the embodiments of this application described and illustrated in the accompanying drawings here can be arranged and designed in various different configurations.
[0038] Therefore, the detailed description of the embodiments of this application provided in the accompanying drawings is not intended to limit the scope of this application claimed, but merely represents selected embodiments of this application. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of this application without creative efforts shall fall within the scope of protection of this application.
[0039] It should be noted that: like reference numerals and letters denote like items in the following drawings. Therefore, once an item is defined in one drawing, it does not need to be further defined and explained in subsequent drawings.
[0040] In the description of the embodiments of this application, it should be noted that if terms such as "upper", "lower", "horizontal", "inner", etc. indicate an orientation or positional relationship based on the orientation or positional relationship shown in the drawings, or the orientation or positional relationship in which the inventive product is customarily placed during use, it is only for the convenience of describing this application and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and thus should not be construed as a limitation of this application. In addition, terms such as "first", "second", etc. are only used for descriptive distinction and cannot be construed as indicating or implying relative importance.
[0041] In addition, if the term "horizontal" appears, it does not mean that the component is required to be absolutely horizontal, but it can be slightly inclined. For example, "horizontal" only means that its direction is more horizontal relative to "vertical", and does not mean that the structure must be completely horizontal, but it can be slightly inclined.
[0042] In the description of the embodiments of the present application, it should also be noted that unless otherwise clearly specified and limited, if the terms "set", "install", "connect", and "link" appear, they should be understood in a broad sense. For example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be directly connected or indirectly connected through an intermediate medium, and it can be the communication inside two components. For those of ordinary skill in the art, the specific meanings of the above terms in the present application can be understood according to specific situations.
[0043] Power flow calculation of a power system is a basic electrical calculation for studying the steady-state operation of the power system. The task is to determine the operating state of the entire system according to the given operating conditions and network structure, such as the voltages (magnitudes and phase angles) of each bus, the power distribution in the network, and the power losses, etc. Power flow calculation is one of the most basic and important calculations in power system analysis and an important part of power system steady-state calculation. The core of power flow calculation is to establish a set of non-linear equations based on the nodal admittance matrix. These equations describe the complex relationships among the node voltages, currents, and powers in the power system. In practical applications, the basic equations of power flow calculation are usually established in real number form of power equations in rectangular coordinates or polar coordinates. Given the network connection and the parameters of each branch, the nodal admittance matrix in power flow calculation can be formed.
[0044] Since power flow calculation involves high-dimensional non-linear equations, it is usually difficult to solve directly. Therefore, various effective calculation algorithms are adopted in the field of power system analysis: (1) Newton-Raphson method: Taylor-expand the node power equations of the power system at a certain initial point and retain the first-order terms to obtain a set of linear equations. By iteratively solving this set of linear equations, the power flow solution of the power system can be gradually approximated. The Newton-Raphson method has the advantages of fast convergence speed and high calculation accuracy in power flow calculation of power systems, but it also has limitations such as the singularity of the Jacobian matrix and sensitivity to the initial value. Due to the complexity of the power system, the Jacobian matrix may be singular or ill-conditioned, and other methods need to be used for processing at this time.
[0045] (2) In addition to the Newton-Raphson method, some researchers have proposed linearization methods. However, when dealing with nonlinear problems through linearization methods, large computational errors will be caused due to the lack of nonlinear characteristics. In recent years, the decomposition calculation of nonlinear algebraic equations has received wide attention. This method decomposes nonlinear problems into linear and nonlinear parts, and expands the nonlinear terms into a series of polynomials that are easy to handle through means such as the Galerkin decomposition method and the recursive decomposition method. This method no longer requires linearization, can retain the essential characteristics of nonlinear problems, and can quickly converge to the exact solution through iterative calculation and step-by-step approximation, showing good application prospects in the field of solving nonlinear algebraic equations.
[0046] The power flow calculation of the power system essentially involves solving a system of nonlinear algebraic equations. Based on this, the present application proposes a power flow calculation method and related device for the power system. This method transforms the power flow equation of the power system, expresses it in the form of constants and nonlinear functions, and further designs calculation formulas for the initial value and iterative update value of the state variables based on the transformed power flow equation. Finally, the results of the state variables in the power flow calculation are obtained by accumulating the initial value and iterative update value of the state variables.
[0047] Based on the above situation, the present application proposes a power flow calculation method and related device for the power system based on recursive decomposition. The following will describe the present application in detail with reference to the embodiments and the drawings.
[0048] As Figure 1 shown, the following is the first flow schematic diagram of the power flow calculation method for the power system based on recursive decomposition in the present application, which may include: S101, obtaining the composition information of the power system network model.
[0049] It should be noted that the composition information of the power system network model may include the structure of the power grid (such as the connection relationship of equipment such as transmission lines and transformers), component parameters (such as resistance, reactance, admittance, etc.), and generation and load parameters (such as active power, reactive power, etc.). After obtaining this information, a mathematical model of the power system can be constructed to provide a basis for subsequent power flow calculations.
[0050] S102, combining the composition information of the power system network model, transforming the power flow calculation equation to be solved for the power system into a state variable expression form combined with constants and nonlinear functions.
[0051] In practical applications, the power flow calculation equation of the power system is usually a set of complex nonlinear equations, involving state variables such as the voltage and power of nodes. These nonlinear equations can be transformed into a state variable expression form combined with constants and nonlinear functions, which is helpful for subsequent recursive decomposition and iterative calculation. The specific transformed equation form can be similar to: F(X)=C+G(X) Among them, F(X) is the power flow calculation equation to be solved, X is the state variable, C is the constant part, and G(X) is the non-linear function part.
[0052] S103. For the power flow calculation equation of the transformed power system, based on the recursive decomposition method, determine the estimation method of the constant part and the iterative calculation formula of the non-linear function part.
[0053] It should be noted that the recursive decomposition method is a technique that decomposes complex problems into simpler problems, and can decompose non-linear equations into a constant part and a non-linear function part. For the constant part, an estimation method needs to be determined, such as direct assignment, prediction based on historical data, etc. For the non-linear function part, an iterative calculation formula needs to be determined for iterative calculation in subsequent steps.
[0054] S104. Perform iterative calculation on the non-linear function part. When the iterative calculation termination condition is reached, integrate the estimated value of the constant part and the iterative updated value of the non-linear function part to obtain the calculation result of the state variable of the power flow calculation equation of the power system.
[0055] It should be noted that the iterative calculation process is a process of repeatedly solving non-linear equations. By continuously adjusting the values of the state variables, it gradually approaches the true solution. During the iterative calculation process, an iterative calculation termination condition can be set, such as the change amount of the state variable is less than a certain threshold, the number of iterations reaches a certain upper limit, etc. When the iterative calculation termination condition is met, integrate the estimated value of the constant part and the iterative updated value of the non-linear function part to obtain the calculation result of the state variable of the power flow calculation equation of the power system. These results include the voltage amplitude and phase angle of each bus node in the power grid, the power distribution of each branch, the power loss of the network, etc.
[0056] This application transforms the power flow calculation equation of the power system to be solved into a state variable expression form combined with constants and non-linear functions, and processes it based on the recursive decomposition method, which can solve non-linear equations more accurately. This method can retain the essential characteristics of non-linear problems and avoid the calculation errors that may be brought by linearization methods. The recursive decomposition method decomposes complex non-linear equations into a constant part and a non-linear function part, and processes them separately, which helps to simplify the calculation process and improve the calculation efficiency. In addition, the way of gradually approaching the true solution through iterative calculation can reduce the calculation amount while ensuring the calculation accuracy. The method of this application does not require excessive simplification or assumptions for the original problem, so it has strong generality and applicability, and can be applied to different types of power systems and different power flow calculation scenarios. The method of this application can provide support for power system analysis and optimization. For example, it can be used to evaluate the system performance under different operating modes, discover the weak links in the power grid, optimize the power supply capacity and access points, etc.
[0057] As Figure 2 shown, the second flow schematic diagram of the power system power flow calculation method based on recursive decomposition in this application may include: S201. Determine the basic information for power system power flow calculation, including node data, node types, and admittance matrix information. Among them, the node types include PQ and PV.
[0058] It should be noted that the PQ node represents the load node in the power grid, and its active power (P) and reactive power (Q) are given known quantities. The voltage amplitude and phase angle of the PQ node are quantities to be solved and need to be obtained through power flow calculation. PQ nodes account for the vast majority in the power system and are usually connected to a large number of loads. The PV node represents the generator node or the node equipped with reactive power compensation equipment, and its active power (P) and voltage amplitude (V) are given known quantities. The reactive power and phase angle of the PV node are quantities to be solved and need to be obtained through power flow calculation. PV nodes are usually fewer in number in the power system but play an important role in maintaining the voltage stability of the system.
[0059] S202. Transform the power system power flow calculation equation and construct it into a form expressed by state variables combined with constants and non-linear functions.
[0060] For the i th PQ node, its state variables are V i , δ i , and its expression form is as follows:
[0061]
[0062] Among them, n is the number of nodes in the power system network model, is the voltage of node i , is the voltage of node k , is the reactive power of node i , is the active power of node i , is the magnitude of the self-admittance of node i , is the phase of the self-admittance of node i , is the phase of the mutual admittance between node i and node k , is the phase of the mutual admittance between node i and node kThe magnitude of the mutual admittance, is the voltage phase of node i , is the voltage phase of node k .
[0063] For the i th PV node, its state variables are δ i , and its expression form is as follows: .
[0064] S203. Based on the reformed power flow equation form of the power system and the recursive decomposition method, an estimation method for the constant part of the state variables and an iterative calculation formula for the non-linear function part are designed.
[0065] Specifically, the following methods and formulas can be adopted: For the i th load node, its state variables are V i , δ i , and the estimation formula for its constant part is as follows:
[0066]
[0067] Among them, is the constant part of , that is, the initial value of , is the constant part of , that is, the initial value of .
[0068] For the i th generating node, its state variables are δ i , and the estimation formula for its constant part is as follows: .
[0069] For the i th load node, its state variables are V i , δ i , and the iterative calculation formula for its non-linear function part is:
[0070]
[0071] Among them, is the th iteration, is the voltage of the node at the i th iteration, is the voltage of the node at the k th iteration, is the voltage phase of the node at the i th iteration, is the voltage phase of the node at the k th iteration.
[0072] For the i th power generation node, its state variable is δ i , and the iterative formula for the non - linear function part is as follows: .
[0073] S204. When the iterative update value reaches the design threshold, integrate the initial value of the constant part and the iterative update value of the non - linear function part, and obtain the calculation result of the state variable of the power flow of the power system after accumulation.
[0074] The specific calculation formula can be: For the i th load node, its state variables are V i , δ i , and the calculation formula is as follows:
[0075] For the i th power generation node, its state variable is δ i , and the calculation formula is as follows: .
[0076] For the power flow calculation method of the present application, on the one hand, there is no need to calculate the Jacobian matrix. When the Newton - Raphson method cannot perform power flow calculation due to the singularity of the Jacobian matrix, the method of the present application can still be used for calculation. On the other hand, the proposed method performs power flow calculation through iteration, and the calculation complexity is relatively low.
[0077] The present application is further illustrated by an example as follows: As Figure 3 shown, it is a diagram of a 3 - node system structure and line parameters. Taking Figure 3Taking the test system as an example, assume that Node 1 in the system is the balancing node, with the voltage magnitude and phase being 1.05 p.u. and 0 rad respectively; Node 2 is the PQ node, with the active power and reactive power being 4.0 p.u. and 2.5 p.u. respectively; Node 3 is the PV node, with the voltage magnitude being 1.04 p.u. and the active power being 2.0 p.u. Figure 3 In it, the line admittances of lines 1-2, 1-3, and 2-3 are 0.02 + j0.04, 0.01 + j0.03, and 0.0125 + j0.025 respectively. G represents the generator. According to the line parameters of the system, the admittance matrix of the system is calculated as follows:
[0078] According to the power flow calculation equation, the state variables to be solved for the system are the voltage magnitude V 2 and phase angle δ 2 of Node 2, and the voltage phase angle δ 3 of Node 3. According to the formula for the constant part, the estimation results of the constant part are as follows:
[0079] According to the formula for the non-linear function part, when m = 1, 2, 3, the calculation results of each state variable are shown in the following table:
[0080] When m = 30, the iterative process ends. The initial value and the results of the non-linear function part are accumulated, and the final calculation results are as follows:
[0081] When using the Newton-Raphson method to calculate the state variables in the power flow of the system, the calculation results are respectively , verifying the accuracy of the proposed method.
[0082] As Figure 4 shown, it is a schematic diagram of a power system power flow calculation system based on recursive decomposition according to the present application, which may include: An information acquisition module, used to acquire the composition information of the power system network model; An equation transformation module, used to transform the power system power flow calculation equation to be solved into a state variable expression form combined with constants and non-linear functions; A decomposition module, used to determine the estimation method of the constant part and the iterative calculation formula of the non-linear function part for the transformed power system power flow calculation equation based on the recursive decomposition method; A calculation module is configured to perform iterative calculations on the non-linear function part. When the termination condition of the iterative calculation is reached, the estimated value of the constant part and the iteratively updated value of the non-linear function part are integrated to obtain the calculation result of the state variables of the power flow calculation equation for the power system.
[0083] In some embodiments of the power flow calculation system for a power system based on recursive decomposition according to the present application, in the equation transformation module, in the expression form of the state variable transformed into a combination of a constant and a non-linear function, for the i th load node and the i th generation node, the expression form of the state variable may include: For the i th load node, the expression form of the state variable is:
[0084]
[0085] Wherein, n is the number of nodes in the power system network model, is the voltage of node i , is the voltage of node k , is the reactive power of node i , is the active power of node i , is the magnitude of the self-admittance of node i , is the phase of the self-admittance of node i , is the phase of the mutual admittance between node i and node k , is the magnitude of the mutual admittance between node i and node k , is the voltage phase of node i , is the voltage phase of node k ; For the i th generation node, the expression form of the state variable is: .
[0086] In some embodiments of the power flow calculation system for a power system based on recursive decomposition according to the present application, in the decomposition module, the estimation method for the constant part may include: For the i th load node, the estimation formula for the constant part is:
[0087]
[0088] Among them, is the constant part of the initial value of is the constant part of the initial value of; For the i th power generation node, the estimation formula for the constant part is: .
[0089] In some embodiments of the power system power flow calculation system based on recursive decomposition in the present application, in the decomposition module, the iterative calculation formula for the non-linear function part may include: For the i th load node, the iterative calculation formula for the non-linear function part is:
[0090]
[0091] Among them, is the th iteration, is the voltage of node at the i th iteration, is the voltage of node at the k th iteration, is the voltage phase of node at the i th iteration, is the voltage phase of node at the k th iteration; For the i th power generation node, the iterative calculation formula for the non-linear function part is: .
[0092] In some embodiments of the power system power flow calculation system based on recursive decomposition in the present application, in the calculation module, the calculation result of the state variables of the power system power flow calculation equation may include: For the i th load node, the calculation result of the state variables includes:
[0093] For the iA power generation node, and the calculation results of the state variables include: 。
[0094] It should be noted that in several embodiments provided in the present application, it should be understood that the disclosed devices and methods can be implemented in other ways. For example, the system embodiments described above are merely illustrative. For example, the division of each module is only a logical function division. In actual implementation, there may be other division methods. For example, multiple modules can be combined or integrated into another device, or some features can be ignored or not executed. The modules described as separate components may or may not be physically separated. The components shown as modules can be one physical unit or multiple physical units, that is, they can be located in one place, or they can be distributed to multiple different places. Some or all of the modules can be selected according to actual needs to achieve the purpose of the solution of this embodiment.
[0095] In addition, in each embodiment of the present invention, each module can be integrated in a processing unit, or each module can exist physically alone, or two or more modules can be integrated in one unit. The above integrated unit can be implemented in the form of hardware or in the form of a software functional unit.
[0096] The embodiments of the present application further provide an electronic device, which may include one or more processors, a memory, and a communication interface.
[0097] Among them, the memory and the communication interface are coupled to the processor. For example, the memory and the communication interface can be coupled together through a bus.
[0098] Among them, the communication interface is used for data transmission with other devices. The memory stores computer program code. The computer program code includes computer instructions. When the computer instructions are executed by the processor, the electronic device executes the steps of the above power flow calculation method based on recursive decomposition.
[0099] Among them, the processor can be a processor or a controller. For example, it can be a Central Processing Unit (CPU), a general-purpose processor, a Digital Signal Processor (DSP), an Application-Specific Integrated Circuit (ASIC), a Field Programmable Gate Array (FPGA), or other programmable logic devices, transistor logic devices, hardware components, or any combination thereof. It can implement or execute various exemplary logic blocks, modules, and circuits described in conjunction with the present disclosure. The processor can also be a combination that realizes computing functions, such as a combination of one or more microprocessors, a combination of a DSP and a microprocessor, and so on. The processor can be used to support the electronic device in executing the method steps provided in the above embodiments.
[0100] Among them, the bus can be a Peripheral Component Interconnect (PCI) bus or an Extended Industry Standard Architecture (EISA) bus, etc. The above bus can be divided into an address bus, a data bus, a control bus, etc.
[0101] A computer-readable storage medium provided by an embodiment of the present application stores a computer program, and when the computer program is executed by a processor, it implements the steps of the above-mentioned power flow calculation method based on recursive decomposition.
[0102] The computer-readable storage medium involved in the present application includes a random access memory (RAM), an internal memory, a read-only memory (ROM), an electrically programmable ROM, an electrically erasable programmable ROM, a register, a hard disk, a removable disk, a CD ROM, or any other form of storage medium well-known in the technical field.
[0103] The above are only the preferred embodiments of the present application and are not used to limit the present application. For those skilled in the art, the present application can have various changes and modifications. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the protection scope of the present application.
Claims
1. A power flow calculation method for a power system based on recursive decomposition, characterized in that, Including: Obtaining the composition information of the power system network model; Combining the composition information of the power system network model, transforming the power flow calculation equation of the power system to be solved into a state variable expression form composed of constants and non-linear functions; For the transformed power flow calculation equation of the power system, based on the recursive decomposition method, determining the estimation method of the constant part and the iterative calculation formula of the non-linear function part; Performing iterative calculation on the non-linear function part. When the iterative calculation termination condition is reached, integrating the estimated value of the constant part and the iterative updated value of the non-linear function part to obtain the calculation result of the state variables of the power flow calculation equation of the power system.
2. The power flow calculation method based on recursive decomposition according to claim 1, wherein The composition information of the power system network model includes at least one of node data, node type, and admittance matrix information.
3. The power flow calculation method based on recursive decomposition according to claim 1, wherein The transformation into the state variable expression form composed of constants and non-linear functions includes: For the i th load node, the expression form of the state variable is: Among them, n is the number of nodes in the power system network model, is the voltage of node i ; is the voltage of node k ; is the reactive power of node i ; is the active power of node i ; is the magnitude of the self-admittance of node i ; is the phase of the self-admittance of node i ; is the phase of the mutual admittance between node i and node k ; is the magnitude of the mutual admittance between node i and node k ; is the voltage phase of node i ; is the voltage phase of node k ; For the i th power generation node, the expression form of the state variable is: 。 4. The power flow calculation method based on recursive decomposition according to claim 1, wherein The estimation method of the constant part includes: For the i th load node, the estimation formula for the constant part is: Among them, is the constant part of, that is the initial value of is the constant part of, that is the initial value of is the active power of node i ; For the i th power generation node, the estimation formula for the constant part is: Among them, is the magnitude of the self-admittance of node i and is the phase of the self-admittance of node i .
5. The power flow calculation method based on recursive decomposition according to claim 1, wherein The iterative calculation formula of the non-linear function part includes: For the i th load node, the iterative calculation formula for the non-linear function part is: Wherein, is the voltage phase of the (m + 1)-th iteration node i , is the node i and the node k is the magnitude of the mutual admittance therebetween, is the node i and the node k is the phase of the mutual admittance therebetween, is the voltage of the (m + 1)-th iteration node i , is the node i 's reactive power, is the node i 's self-admittance magnitude, is the node i 's self-admittance phase, is the node i 's voltage, is the -th iteration, is the -th iteration node i 's voltage, is the -th iteration node k 's voltage, is the -th iteration node i 's voltage phase, is the -th iteration node k 's voltage phase; For the i th power generation node, the iterative calculation formula for the non-linear function part is as follows: 。 6. The power flow calculation method based on recursive decomposition according to claim 1, characterized in that The calculation result of the state variables of the power flow calculation equation of the power system includes: For the i th load node, the calculation results of the state variables include: Among them, is the voltage phase of the node i . is the voltage phase of the k -th iteration node i . is the voltage of the node i , is the voltage amplitude of the k -th iteration node i . For the i th power generation node, the calculation results of the state variables include: 。 7. A power flow calculation system for a power system based on recursive decomposition, characterized in that, Including: An information acquisition module for obtaining the composition information of the power system network model; An equation transformation module for combining the composition information of the power system network model and transforming the power flow calculation equation of the power system to be solved into a state variable expression form composed of constants and non-linear functions; A decomposition module for, based on the recursive decomposition method, determining the estimation method of the constant part and the iterative calculation formula of the non-linear function part for the transformed power flow calculation equation of the power system; A calculation module for performing iterative calculation on the non-linear function part. When the iterative calculation termination condition is reached, integrating the estimated value of the constant part and the iterative updated value of the non-linear function part to obtain the calculation result of the state variables of the power flow calculation equation of the power system.
8. The power flow calculation system based on recursive decomposition according to claim 7, wherein In the equation transformation module, the transformation into the state variable expression form composed of constants and non-linear functions includes: For the i th load node, the state variable expression is: Among them, n is the number of nodes in the power system network model, is the voltage of node i ; is the voltage of node k ; is the reactive power of node i ; is the active power of node i ; is the magnitude of the self - admittance of node i ; is the phase of the self - admittance of node i ; is the phase of the mutual - admittance between node i and node k ; is the magnitude of the mutual - admittance between node i and node k ; is the voltage phase of node i ; is the voltage phase of node k ; For the i th power generation node, the expression form of the state variable is: 。 9. The power flow calculation system based on recursive decomposition according to claim 7, characterized in that, In the decomposition module, the estimation method of the constant part includes: For the i th load node, the estimation formula for the constant part is: Among them, is the constant part of, that is the initial value of, is the constant part of, that is the initial value of, is the active power of node i ; For the i th power generation node, the estimation formula for the constant part is: Among them, is the magnitude of the self-admittance of node i and is the phase of the self-admittance of node i .
10. The power flow calculation system based on recursive decomposition according to claim 7, characterized in that, In the decomposition module, the iterative calculation formula of the non-linear function part includes: For the i th load node, the iterative calculation formula for the non-linear function part is: Among them, is the voltage phase of the (m + 1)-th iteration node i , is the magnitude of the mutual admittance between node i and node k , is the phase of the mutual admittance between node i and node k , is the voltage of the (m + 1)-th iteration node i , is the reactive power of node i , is the magnitude of the self-admittance of node i , is the phase of the self-admittance of node i , is the voltage of node i , is the -th iteration, is the -th iteration voltage of node i , is the -th iteration voltage of node k , is the -th iteration voltage phase of node i , is the -th iteration voltage phase of node k ; For the i th power generation node, the iterative calculation formula for the non-linear function part is: 。 11. The power flow calculation system based on recursive decomposition according to claim 7, characterized in that, In the calculation module, the calculation result of the state variables of the power flow calculation equation of the power system includes: For the i th load node, the calculation results of the state variables include: wherein, is the voltage phase of the node i ; is the voltage phase of the node k at the i th iteration; is the voltage of the node i ; is the voltage amplitude of the node k at the i th iteration. For the i th power generation node, the calculation results of the state variables include: 。 12. An electronic device, characterized in that, Including: A memory for storing computer programs; A processor for implementing the steps of the power system power flow calculation method based on recursive decomposition according to any one of claims 1 to 6 when executing the computer program.
13. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program, and when the computer program is executed by the processor, it implements the steps of the power system power flow calculation method based on recursive decomposition according to any one of claims 1 to 6.