Power flow convergence method and system based on direct current power flow preprocessing and staged decoupling

Abstract

The invention provides a power flow convergence method and system based on direct current power flow preprocessing and staged decoupling, and belongs to the field of power system analysis and control. Comprising the following steps: performing multi-stage verification and parameter adjustment on original power flow data, executing active balance preliminary verification of a power system, and generating a node active power set initial value set; setting an initial value set based on the active power of the nodes to carry out direct-current power flow calculation, obtaining an initial phase angle of each node, identifying a heavy-load line with a power angle difference exceeding a preset threshold value, and carrying out adjustment to obtain initial phase angle distribution; in combination with a preset initial voltage amplitude and the initial phase angle distribution, performing staging decoupling load flow calculation and adjustment to obtain load flow parameters; and taking the power flow parameter as an initial value, and carrying out full-coupling power flow calculation to obtain a final power flow solution meeting the convergence precision requirement. According to the method, the physical root of the convergence problem can be intelligently diagnosed and positioned, and the convergence and robustness of load flow calculation under complex working conditions are effectively improved through small-range targeted adjustment.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of power system analysis and control, and particularly relates to a power flow convergence method and system based on DC power flow preprocessing and stage-by-stage decoupling. BACKGROUND

[0002] The statements in this section merely provide background information related to the application and do not necessarily constitute prior art.

[0003] With the deepening of the "double carbon" goal, the global energy structure is undergoing a deep transformation, and the new power system presents the core characteristics of large-scale and high-proportion new energy grid connection. Wind power, solar power and other new energy power generation technologies have become the main body of incremental power systems due to their clean and low-carbon advantages. However, the output of such power sources is significantly affected by natural conditions, and has strong volatility, intermittency and spatio-temporal uncertainty, which is quite different from the stable output characteristics of traditional synchronous generators.

[0004] This characteristic leads to extremely frequent and intense dynamic adjustment of the power grid operation mode. In the scenario of large fluctuations in new energy output and increasing cross-regional power exchange, the active and reactive power distribution of the power system is often in an extreme spatio-temporal imbalance, making the solution of the power flow equation fall into a sick interval, the condition number of the Jacobian matrix increases, and the numerical stability decreases, posing a serious challenge to traditional power flow calculation methods. In the current power grid operation, extreme conditions such as overloaded section power, weak local voltage support and limited system balance constraints frequently occur, directly leading to the numerical power flow calculation deviating from the convergence domain, causing frequent oscillation of numerical iteration or even complete divergence. This problem not only affects the efficiency of core businesses such as power grid planning, design and operation, but also delays fault handling, situation awareness and other critical decisions due to the inability to effectively solve the power flow, threatening the safe and stable operation of the power system.

[0005] Existing power flow convergence enhancement techniques have significant limitations: on the one hand, most methods only optimize algorithm steps from the numerical iteration level, lack deep exploration of the physical nature of power flow divergence under fluctuating operation modes, and are difficult to accurately locate the weak link causing calculation divergence; on the other hand, the full-coupling solution mechanism of the traditional Newton-Raphson method is prone to mutual interference of correction amounts in each dimension due to the ill-conditioned Jacobian matrix in extreme conditions with high penetration of new energy, resulting in numerical coupling oscillation problems; while the traditional PQ decomposition method has decoupling, but its alternating iteration mechanism often causes P-θ and Q-U dimensional residual alternating mapping, leading to convergence path detours and even divergence problems. SUMMARY

[0006] To overcome the above deficiencies of the prior art, the present application provides a power flow convergence method and system based on DC power flow preprocessing and stage-by-stage decoupling, aiming to solve the problems of the traditional power flow calculation method in the face of complex operating conditions, such as the nonlinear coupling degree of power flow equations being aggravated, the ill-conditioning of the Jacobian matrix, and the iterative value being easily out of the convergence domain and the numerical solution being difficult to converge stably, due to the large-scale access of new energy and its uncertainty in output, resulting in the dramatic fluctuation of power grid operating conditions.

[0007] To achieve the above object, one or more embodiments of the present application provide the following technical solutions: The present application provides a power flow convergence method based on DC power flow preprocessing and stage-by-stage decoupling in the first aspect; The power flow convergence method based on DC power flow preprocessing and stage-by-stage decoupling comprises: The original power flow data is subjected to multi-level checking and parameter adjustment, and a preliminary check of system active power balance is performed to generate a node active power setting initial value set satisfying basic physical feasibility; Based on the node active power setting initial value set, DC power flow calculation is performed to obtain initial phase angles of each node, identify overloaded lines with power angle difference exceeding a preset threshold, and make all line power angle differences satisfy the constraint through local active resource adjustment to obtain an initial phase angle distribution; Combined with a preset initial voltage amplitude and the initial phase angle distribution, stage-by-stage decoupled power flow calculation and adjustment are performed to obtain power flow parameters; Based on the obtained power flow parameters as initial values, full-coupled power flow calculation is performed to obtain a final power flow solution satisfying the convergence accuracy requirement.

[0008] As a further technical solution, the multi-level checking and parameter adjustment of the original power flow data, and the preliminary check of system active power balance to generate a node active power setting initial value set satisfying basic physical feasibility, comprises: The original power flow data is subjected to network topology checking and device parameter checking, the network topology checking comprises identifying and removing isolated buses, checking the configuration of balancing machines and the consistency of circuit breaker states; the device parameter checking comprises checking the non-negativity and reasonableness of generator, transformer and line parameters, and the reasonableness of node power and voltage setting values; The initial total active imbalance of the system is calculated, and it is judged whether the power to be borne by the balancing machine exceeds its preset upper and lower limits; if the balancing machine power exceeds the limit, the excess power is proportionally distributed to other controllable generators in the system except the balancing machine to obtain a node active power setting initial value set satisfying basic physical feasibility.

[0009] As a further technical solution, the DC power flow calculation based on the node active power setting initial value set to obtain initial phase angles of each node comprises: Substituting the active power injection power of all nodes in the initial set of active power settings into the DC power flow equation matrix, as shown below:

[0010] in, This represents the active power injection amount for all unbalanced nodes. The phase angles to be solved for all non-equilibrium nodes; Here is the admittance matrix; Solving the linear equation yields the initial phase angle of each non-equilibrium node. Setting the phase angle of the equilibrium node to 0 gives the initial phase angle of all nodes in the network.

[0011] As a further technical solution, the identification of heavily loaded lines with a power angle difference exceeding a preset threshold includes: Based on the DC power flow calculation results, calculate the absolute value of the power angle difference between the two ends of all lines; The absolute value of the power angle difference is compared with a preset static stability safety threshold to identify lines that exceed the limit; Using the sending and receiving busbars of each over-limit line as the center, a local topology adjustment area is delineated using a breadth-first search algorithm; According to the preset priority order, the distribution of active resources in the local topology area is adjusted, and the active power output of the opposite area is adjusted synchronously to maintain the overall active power balance of the system. Recalculate the DC power flow until the absolute value of the power angle difference of all lines is less than the static stability safety threshold.

[0012] As a further technical solution, the phased decoupling of power flow calculation and adjustment includes: Phase 1: Fixed voltage amplitude, only... P-θ The active power flow is solved iteratively. If the iteration fails to converge, the weak bus is located based on the maximum phase angle correction, and topology search and active power resource adjustment are performed on its associated region until... P-θ Iterative convergence; Second stage: Based on the phase angle distribution obtained in the first stage, fix the voltage phase angle and only perform... Q-U The reactive power flow is solved iteratively. If the iteration fails to converge, the weak bus is located based on the maximum voltage correction, and a topology search and reactive power resource adjustment are performed on its associated region until... Q-U Iterative convergence.

[0013] As a further technical solution, the method involves locating the weak busbar based on the maximum phase angle correction, and then performing topology search and active power resource adjustment on its associated region until... P-θ Iterative convergence includes: when P-θIf the iteration reaches the preset maximum number of iterations and still fails to converge, extract the maximum phase angle correction amount within the iteration cycle that caused the non-convergence and its corresponding bus node. Starting from the aforementioned busbar, a breadth-first search with a finite number of steps is used to determine the key adjustment area; Perform at least one of the following operations on the critical adjustment area: verify and correct line and transformer parameters, adjust the active power output of generators in the area according to the principle of sensitivity or proportionality, and reduce non-critical loads in the area.

[0014] As a further technical solution, the method involves locating weak buses based on the maximum voltage correction amount and performing topology search and reactive power resource adjustment on their associated areas until... Q-U Iterative convergence includes: when Q-U If the iteration fails to converge after reaching the preset maximum number of iterations, extract the maximum voltage correction amount within the iteration period that caused the non-convergence and its corresponding bus node. Starting from the aforementioned busbar, a breadth-first search with a finite number of steps is used to determine the key adjustment area; The reactive power status of the region is determined by the sign of the maximum voltage correction. If the maximum voltage correction is positive, capacitive reactive power equipment is put into operation or inductive reactive power equipment is cut off according to priority. If the maximum voltage correction is negative, inductive reactive power equipment is put into operation or capacitive reactive power equipment is cut off according to priority.

[0015] As a further technical solution, the step of performing fully coupled power flow calculations based on the obtained power flow parameters as initial values ​​to obtain a final power flow solution that meets the convergence accuracy requirements includes: Using the power flow parameters obtained after phased decoupling and adjustment as the initial values ​​for iteration, the active power imbalance and reactive power imbalance of each node are calculated in the k-th iteration. Construct the Jacobian matrix or coefficient matrix at the current iteration point, solve the correction equation to obtain the phase angle correction and voltage magnitude correction for each node, and correct the voltage magnitude and phase angle of each node respectively; Determine whether the active power imbalance and reactive power imbalance at each node satisfy the convergence condition. If both conditions are met, output the final power flow solution; otherwise, set... k = k +1, continue iterating until the convergence accuracy requirement is met.

[0016] The second aspect of the present invention provides a power flow convergence system based on DC power flow preprocessing and staged decoupling.

[0017] A power flow convergence system based on DC power flow preprocessing and staged decoupling includes: The data verification module is configured to perform multi-level verification and parameter adjustment on the raw power flow data, and to perform preliminary verification of the system active power balance to generate a set of initial values ​​for node active power settings that meet basic physical feasibility. The DC power flow preprocessing module is configured to: perform DC power flow calculation based on the set of initial values ​​of the active power of the nodes, obtain the initial phase angle of each node, identify the heavily loaded lines whose power angle difference exceeds a preset threshold, adjust the local active power resources to make the power angle difference of all lines meet the constraints, and obtain the initial phase angle distribution. The phased decoupled power flow calculation and adjustment module is configured to: combine the preset initial voltage amplitude with the initial phase angle distribution to perform phased decoupled power flow calculation and adjustment to obtain power flow parameters; The fully coupled power flow verification module is configured to perform fully coupled power flow calculations based on the obtained power flow parameters as initial values, and obtain the final power flow solution that meets the convergence accuracy requirements.

[0018] A third aspect of the present invention provides a computer-readable storage medium having a program stored thereon, which, when executed by a processor, implements the steps in the power flow convergence method based on DC power flow preprocessing and phased decoupling as described in the first aspect of the present invention.

[0019] A fourth aspect of the present invention provides an electronic device, including a memory, a processor, and a program stored in the memory and executable on the processor, wherein the processor executes the program to implement the steps in the power flow convergence method based on DC power flow preprocessing and phased decoupling as described in the first aspect of the present invention.

[0020] The above one or more technical solutions have the following beneficial effects: (1) This invention eliminates parameter configuration errors through data verification, then quickly locates and adjusts heavily loaded lines by utilizing the linearization characteristics of DC power flow, and finally precisely pinpoints specific weak points through the correction amount analysis in the phased decoupling iteration. This diagnostic path from macro to micro and from parameters to physics completely changes the traditional method's reliance on manual experience for trial and adjustment, and realizes intelligent diagnosis and precise location of the root cause of the problem.

[0021] (2) High-quality initial phase angle values ​​are obtained through DC power flow preprocessing, providing a good starting point for subsequent nonlinear iterations. A fully phased decoupling strategy of "global active power first, then local reactive power" is adopted, effectively cutting off the traditional PQ decomposition method. P-θ and Q-UThe iterative process involves cross-coupling oscillation paths. When any iteration fails to converge, a regional adjustment mechanism based on the maximum correction amount can implement precise intervention within a minimal scope, avoiding blind disturbances to the entire system. Finally, the high-precision approximate solution obtained in the preceding stages is used as the initial value, enabling the fully coupled power flow calculation to converge rapidly within a very small number of iterations. The entire scheme forms a complete closed loop from preprocessing and decoupling adjustment to global verification, significantly improving the convergence and numerical robustness of power flow calculations under complex conditions such as renewable energy fluctuations and heavy load operation.

[0022] (3) In the DC preprocessing stage, this invention adjusts active power resources only in local areas of lines with abnormal power angle differences; in the phased decoupling stage, the adjustment area with the smallest impact range is precisely delineated through a finite number of topology searches based on the maximum phase angle or voltage correction. The adjustment methods are prioritized according to equipment response speed and adjustment cost, ensuring the economy and speed of adjustment. This mechanism can solve the convergence problem while maintaining the stability of the original system operating state to the greatest extent, providing reliable and efficient decision support for dispatchers.

[0023] Advantages of additional aspects of the invention will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. Attached Figure Description

[0024] The accompanying drawings, which form part of this invention, are used to provide a further understanding of the invention. The illustrative embodiments of the invention and their descriptions are used to explain the invention and do not constitute an improper limitation of the invention.

[0025] Figure 1 This is a flowchart of the method in the first embodiment.

[0026] Figure 2 This is a schematic diagram of the DC power flow preprocessing flow of the first embodiment.

[0027] Figure 3 This is a schematic diagram of the secondary topology search in the first embodiment.

[0028] Figure 4 For the first embodiment P-θ A schematic diagram of the iterative active power flow adjustment process.

[0029] Figure 5 For the first embodiment Q-U A schematic diagram of the iterative reactive power flow adjustment process.

[0030] Figure 6 This is a schematic diagram of the power flow path and convergence region in the first embodiment.

[0031] Figure 7This is a system structure diagram of the second embodiment. Detailed Implementation

[0032] It should be noted that the following detailed descriptions are exemplary and intended to provide further illustration of the invention. Unless otherwise specified, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention pertains.

[0033] It should be noted that the terminology used herein is for the purpose of describing particular implementations only and is not intended to limit the exemplary implementations of the present invention.

[0034] Where there is no conflict, the embodiments and features in the embodiments of the present invention can be combined with each other.

[0035] The overall approach proposed in this invention is as follows: First, through multi-level data verification and preliminary active power balance check, initial value deviations caused by unreasonable original parameter configurations and active power supply-demand imbalances are eliminated; second, high-quality initial phase angle values ​​are quickly constructed using linearized solutions of the DC power flow model, and the identification and adjustment of overloaded sections are achieved based on the power angle difference threshold warning; subsequently, following the physical logic of "global active power first, then local reactive power," the following steps are performed respectively... P-θ and Q-U The process involves independent iterations, and the system's weak points are located using the maximum phase angle correction and maximum voltage correction. A regional power flow convergence improvement strategy is implemented using breadth-first search to guide the power flow back to the convergence domain. Finally, the high-quality approximate solution obtained after phased decoupling and adjustment is used as the starting point for iteration, and the final verification is performed through fully coupled power flow calculation.

[0036] Example 1 This embodiment discloses a power flow convergence method based on DC power flow preprocessing and staged decoupling; like Figure 1 As shown, the power flow convergence method based on DC power flow preprocessing and staged decoupling includes: Step S1 involves performing multi-level verification and parameter adjustment on the original power flow data, and conducting a preliminary check of the system's active power balance to generate a set of initial node active power settings that meet basic physical feasibility requirements.

[0037] First, network topology verification is carried out. A power grid topology traversal algorithm is used to identify and remove isolated buses without connected branches to ensure that all operating nodes are connected to the main grid. Verification and assurance are made that each regional sub-grid is equipped with at least one balancing machine. If there is a case where no balancing machine is configured, the synchronous generator with the largest capacity in the region is automatically designated as the balancing machine. The status of the circuit breakers at both ends of each line is compared. If there is an inconsistency where one end is closed and the other end is open, it is corrected to a unified status according to the power grid operation specifications.

[0038] Then, equipment parameter verification was performed, checking the non-negativity of the impedance parameters of all generators and transformers, and simultaneously verifying the lines and transformers. X / R Check whether the parameters are within the reasonable range of the project; check the sign of the reactive power compensation equipment capacity and correct capacitive and inductive configuration errors; verify the active and reactive power configurations of all PQ nodes to ensure that their power factors are within the reasonable range and avoid extreme configurations; verify the voltage settings of all PV nodes to ensure that the voltage settings of PV nodes are within the acceptable range allowed by the power grid operation standards.

[0039] Furthermore, a preliminary verification of the system's active power balance is conducted, including: active power balance check, power flow problem localization, and strategies to improve power flow convergence. Specifically: When dealing with power flow problems in large power grids, it is necessary to first verify the active power balance constraints of the system to ensure that the total active power supply and demand of the entire system are roughly balanced. If there is a huge difference between the total generation and the total load in the initial stage, the balancing machine will not be able to handle this excessively large unbalanced power, which will directly lead to the failure of subsequent power flow calculations.

[0040] First, the initial unbalanced power of the system is calculated, as shown in the following formula: (1) In the formula, This refers to the power imbalance of the system. The total load power of the system. For system line loss, Given the total power generation of the system, the system line loss formula is as follows: (2) in, μ This is a proportionality coefficient, typically around 1%-5%. Then, calculate the total power that the system balancing machine should handle: (3) In the formula, Set the power for the balancing machine. The total power required for the balancing machine.

[0041] Secondly, in the process of locating the trend problem, all balancing mechanisms are traversed to determine the responsibilities they need to bear. Are the values ​​within the preset maximum and minimum upper and lower limits of the balancing machine? If the power limit of the balancing machine is outside the range, it means that the initial active power configuration for the current operating condition is physically infeasible and there is a risk of exceeding the limit. Therefore, the correction process must be initiated.

[0042] Once the balancing machine's power is determined to be exceeding the limit, a power flow convergence improvement strategy is adopted, initiating an adaptive adjustment mechanism based on minimum intervention and proportional allocation to reduce the balancing machine's power exceeding the limit. The active power is distributed reasonably to other controllable generators within the system to achieve a new active power balance.

[0043] (4) By selecting generators within the region, excluding balancing machines, whose current active power output has not reached its power limit, an adjustable generator set is formed. The power exceeding the limit will be located. The power is allocated to the adjustable generator set according to the original generator node power ratio. Each generator in the system. The allocation formula is as follows: (5) in, This represents the adjusted active power of the k-th generator. Let be the original active power of the k-th generator; This represents the active power adjustment of the k-th generator.

[0044] The final output is a set of initial values ​​for the nodal active power settings that satisfy the basic active power balance. It is used for subsequent DC power flow preprocessing.

[0045] Step S2: Based on the initial set of active power values ​​for the nodes, perform DC power flow calculations to obtain the initial phase angle of each node. Identify heavily loaded lines whose power angle difference exceeds a preset threshold. Adjust local active power resources to ensure that the power angle difference of all lines meets the constraints, thus obtaining the initial phase angle distribution, as shown below. Figure 2 As shown.

[0046] Step S21: After completing the preliminary balance check of the system's active power and obtaining the set of active power setpoints that satisfy the constraints. Then, the rapid computational advantage of the DC power flow model is used to determine the initial phase angle distribution of the system. This model is based on simplifying assumptions and linearizes the complex nonlinear power flow equations.

[0047] The following assumptions are made: (1) The voltage at each node of the power system is usually near the rated voltage, and is approximately considered to be .

[0048] (2) The voltage phase angle difference between the two ends of the line is very small, so it is assumed that... Therefore, there is , .

[0049] (3) In ultra-high voltage power grids, the line resistance is much smaller than the reactance, therefore the resistance is negligible. .

[0050] Based on the above assumptions, nodes in the traditional AC power flow equation i To the node j active power transmitted It is linearized and depends only on the phase angle difference between the two ends and the line reactance. Its calculation formula is as follows: (6) In the formula, For nodes i To the node j The active power transmitted, The difference in phase angle between the two ends, For line reactance, This is the line susceptance.

[0051] Will Active power injection of all nodes Substituting the values, we obtain the matrix form of the DC power flow equations for the entire system: (7) In the formula, This represents the active power injection amount for all unbalanced nodes. The phase angles to be solved for all non-equilibrium nodes are given, and the phase angles for the equilibrium nodes are given. Set the reference value to 0. Let be the admittance matrix, with its diagonal and off-diagonal elements as follows: (8) in, The self-susceptance between nodes.

[0052] By solving this simplified linear equation, the initial phase angle of each node can be calculated very quickly.

[0053] Step S22: Locate the power flow problem and identify heavily loaded lines whose power angle difference exceeds a preset threshold.

[0054] Calculate the absolute value of the phase angle difference between the two busbars at both ends of all transmission lines in the DC power flow solution. : (9) In DC power flow diagnosis, the power angle difference between the two busbars at both ends of the line is considered. Compared with the preset static stability threshold A comparison is made. This threshold is usually set in the range of 30° to 45°, which is significantly less than the theoretical static stability limit, and serves as a safety warning line for power flow operation.

[0055] If the absolute value of the phase angle difference of a certain transmission line Exceeding the preset threshold The system immediately determines that the line is under severe overload or that there is a fundamental active power imbalance in its region. This large power angle difference means that its static stability margin is insufficient, and the power flow equations in this region are highly nonlinear, which may lead to convergence failure in subsequent staged decoupling calculations. The system will immediately initiate a local active power adjustment strategy.

[0056] Step S23: Taking the two busbars at both ends of the AC transmission line with abnormal power angle as the center, a breadth-first search (BFS) algorithm with a finite number of steps (two or three levels of neighborhood) is used, such as... Figure 3 As shown. Quickly define a local topology region that includes all adjacent power supply and load nodes of the sending node. and the local topology region containing all adjacent nodes of the receiving node. .

[0057] In local topological regions as well as Internally, the system systematically identifies all controllable active power devices, including synchronous generators, energy storage, wind power, and photovoltaics, and obtains their current active power output in real time. and its effective upper and lower limits The adjustability margins of various types of equipment are calculated.

[0058] After completing resource identification and margin assessment of all adjustable active power equipment within the local topology area, power flow convergence improvement measures are initiated. The core objective of these measures is to adjust the active power balance within the area through the principles of minimal intervention and fastest response speed, thereby eliminating overload on lines with abnormal phase angles and... Reduced to a safe threshold γ To avoid concentrating the effects of adjustments on the balancing bus, if active power injection is increased at the receiving end of an abnormal line, the system should correspondingly reduce generator output at the sending end of the transmission line to maintain a balance between total system power generation and total load demand.

[0059] To adjust the receiving end Taking the region's merits as an example, the focus of the adjustment will first be on the region. , The power supply within the system. Given the rapid and bidirectional regulation capabilities of energy storage, the first priority is to dispatch energy storage devices in local areas to fully utilize the system's energy storage resources.

[0060] When energy storage regulation capacity is insufficient, the second stage begins, adjusting the output of the synchronous generator, which serves as the primary power source. The adjustment amount is constrained within its adjustable margin range and allocated according to the current generator output ratio.

[0061] If, after the above two-stage power flow convergence improvement measures, the power angle difference of the abnormal line still has not been reduced to below the safety threshold, the output of wind power and photovoltaic units will be adjusted according to the current output ratio.

[0062] If the line is still overloaded, then identify the local area at the receiving end. Standby generator sets that are out of service but not under maintenance should be set to start status.

[0063] If all power regulation measures have been implemented, but the abnormal line is still overloaded, then the final safeguard measure should be activated, and the receiving end should be disconnected. Partial load on the region.

[0064] Adjusting the receiving end After the region's active power output is completed, the sending end should be adjusted accordingly based on the above priority order. The active power output of the region. After implementing any stage of power flow convergence improvement measures, this invention will recalculate the DC power flow and check the power angle difference of all branches. This adjustment-calculation cycle will continue until the power flow converges and the power angle difference of all lines meets the threshold constraint. This concludes the preprocessing based on DC power flow.

[0065] The final phase angles of each node obtained from the DC power flow preprocessing Vector and preset voltage amplitude This combination provides a high-quality initial power flow value with a reasonable phase angle distribution for subsequent phased decoupling iterations. .

[0066] Step S3: Combine the preset initial voltage amplitude with the initial phase angle distribution to perform phased decoupled power flow calculation and adjustment to obtain power flow parameters.

[0067] Unlike the standard PQ decomposition method P-θ and Q-U The mechanism of alternating correction within a single large cycle employs a decoupling strategy with absolute temporal independence in this embodiment. The system blocks the numerical interaction between the two by controlling logic parameters, only when the voltage amplitude is locked. P-θ Phase locking is only enabled after the iteration converges globally. Q-U Iteration. This mechanism prevents unconverged residuals from repeatedly coupling and oscillating between the two iterative dimensions due to alternating corrections.

[0068] The phased decoupling power flow convergence improvement described in this embodiment is achieved through logical reconstruction of the standard PQ decomposition method, realizing a three-mode switching operation mechanism based on variable instructions. By introducing two core parameters, 'Decoupled P Solution Logic' and 'Decoupled Q Solution Logic', the system can flexibly switch algorithm modes according to the calculation stage. When 'Decoupled P Solution Logic' is set to True and 'Decoupled Q Solution Logic' is set to False, the solver is a single... P-θ Iterative active power flow calculation; when 'Decoupled P Solution Logic' is False and 'Decoupled QSolution Logic' is True, the solver is a single... Q-U Iterative reactive power flow calculation; when both are set to True, the algorithm logic is the traditional PQ decomposition method.

[0069] To clarify the computational principles of each mode, the power correction equations for the nodes are given below: (10) In the formula, This refers to the imbalance of active power. This refers to the reactive power imbalance. This is the sensitivity matrix of active power to voltage phase angle; This is the sensitivity matrix of active power to voltage amplitude; This is the sensitivity matrix of reactive power to voltage phase angle; This is the sensitivity matrix of reactive power to voltage amplitude; This is the phase angle correction amount; This is the voltage correction amount; This represents the voltage amplitude.

[0070] The elements of the Jacobian matrix are expressed as follows: (11) (12) (13) (14) (15) (16) (17) (18) in, for H The off-diagonal elements of the matrix; for N The off-diagonal elements of the matrix; for J The off-diagonal elements of the matrix; for L The off-diagonal elements of the matrix; for H The diagonal elements of a matrix; These are the diagonal elements of matrix N; for J The diagonal elements of a matrix; for L The diagonal elements of a matrix; For connecting nodes i With nodes j The negative values ​​of the conductance of the branches between them; for i Nodes and j The phase angle difference between nodes.

[0071] Because high-voltage transmission networks typically have high... X / R Due to the characteristics of the ratio, the branch susceptance is much greater than the conductance. Based on the expressions for each element in the Jacobian matrix in Equation 11-18, the correction matrix for active power on voltage amplitude can be obtained. N The correction matrix for voltage amplitude is much lower than that for reactive power. L Similarly, the reactive power correction matrix for voltage phase angle... J The ones are also much smaller than H Matrix. Therefore P-U and Q-θ The nodes are weakly coupled, and the PQ decomposition method simplifies the node voltage correction equation to the following formula.

[0072] (19) The traditional PQ decomposition method will 2( n The system of linear equations of order -1 becomes two n The system of linear equations with coefficients of -1 is simplified to a symmetric matrix that remains invariant during iteration. and .

[0073] (20) (twenty one) The nodal power imbalance expressed in polar coordinates is: (twenty two) (twenty three) In the formula, For the first i Active power imbalance at nodes; For the first i The given active power of the node; For the first i The reactive power imbalance at the node; For the first i The given reactive power of the node.

[0074] In the iterative process of the traditional PQ decomposition method, the phase angle correction amount... and voltage correction amount The updates continue to alternate within the same outer loop. While this alternating update mechanism is faster than the fully coupled Newton-Raphson method, it may fail under heavy load conditions if... P-θ It has not yet subsided, and what followed... Q-U Iteration generates new voltage correction values The newly revised voltage is obtained through Equation 22. This can have an impact, causing the convergence path to become circuitous or to oscillate at the convergence edge.

[0075] To avoid the disturbance to convergence caused by this alternating update mechanism, in this embodiment, the power flow calculation is decomposed into two stages, wherein the first stage is... P-θ Iterative phase, the second phase is Q-U The iterative phase consists of two completely independent yet continuous phases.

[0076] Since active power balance involves all generators and loads in the system, its imbalance... Adjustments need to be made through the active power sources across the entire network, exhibiting a global nature. Reactive power, on the other hand, is primarily used to support voltage, and its imbalance... Typically, regulation is achieved through locally configured reactive power sources, exhibiting significant locality. Based on these characteristics, a strategy is adopted to address global issues first, followed by local issues, starting with... P-θ Iterative improvement of active power flow convergence, followed by further improvement based on... Q-U The convergence of reactive power flow is improved through iteration.

[0077] Step S31, based on P-θ Iterative active power flow adjustment, such as Figure 4 As shown.

[0078] Step S311: In the phased decoupling power flow convergence improvement, in the first stage, the solver parameters 'Decoupled P Solution Logic' are set to True and 'Decoupled Q Solution Logic' to False. The voltage magnitude remains constant during the iteration process, and only the active power equation is considered for iterative calculation of the active power flow. The calculation steps are as follows: (1) Substitute the initial power flow values ​​obtained from DC power flow preprocessing Set, keeping the voltage amplitude constant during the iteration process, let the number of iterations be... n = 0; (2) Substitute into Equation 22 to calculate the active power imbalance at each node. and find ; (3) Solve Formula 20 to obtain the voltage phase angle correction for each node. ; (4) Correct the phase angle of the voltage at each node. : (twenty four) (5) Judgment n If the maximum number of iterations has been reached, proceed to step (7); otherwise, proceed to step (6). (6) Determine the active power imbalance at each node. Check if the convergence condition is met. If it is met, proceed to step (7); if not, set the iteration count to zero. n = n +1, and return to step (2) for iteration; (7) Output the final power flow data set.

[0079] Step S312: Locating the power flow problem where the active power does not converge.

[0080] when P-θ The maximum active power imbalance within a preset maximum number of iterations. Still unable to meet the preset convergence accuracy When this occurs, it indicates that the system operating conditions corresponding to the current scheduling settings may have approached or exceeded the static stability limit, and the calculation model cannot find a solution under the current settings.

[0081] In this embodiment, if the power flow calculation fails to converge after reaching the maximum number of iterations, a correction based on the maximum phase angle is performed. Problem identification and regional adjustment strategies. The essence of active power transmission in a network is driven by the phase angle difference between nodes. During the iteration process, the active power imbalance... The power residual only reflects the current node and is a numerical indicator characterizing the balance gap; while the phase angle correction... This reflects the degree of phase angle correction the system attempts to eliminate the difference, and is an indicator of its balancing ability. Phase angle correction amount The larger the value, the greater the active power balancing capacity required by this node, and the easier it is to reach its node balancing capacity limit. Compared to , This reflects the fundamental reason why the system's active power does not converge.

[0082] locking The bus in question is actually located at the physical weak point in the power grid where active power dispatching resistance is greatest. Using this bus as the starting point, the adjustment area can be delineated using BFS, which can directly locate the cause of iterative divergence, thereby guiding the power flow calculation from the edge of the physical limit back to the effective convergence region.

[0083] In each iteration of the calculation, the system records the active power imbalance of all nodes in the network in real time. With phase angle correction Simultaneously extract and output the maximum phase angle correction within the current iteration cycle. and its corresponding bus nodes. Focus on observing the period approaching the maximum iteration count. Bus node. If the maximum phase angle correction is always located at a specific bus in multiple consecutive iterations, it indicates that the active power balance capacity at that node has reached the physical limit, and the bus can be determined as the physical weak point that causes the overall power flow calculation to diverge.

[0084] Step S313: Strategies to improve the convergence of power flow that does not converge with active power.

[0085] First, a finite-level BFS topology search program (typically a second or third-level neighborhood search) is run based on the magnitude of the correction to locate the adjacent regions of the bus node. Second, the model parameters of this region are validated, with a focus on checking the impedance parameters of the lines and transformers within that region. If a branch reactance is identified... X If the value far exceeds the physical norm, a standard parameter library will be called for comparison and correction to eliminate active power obstruction caused by model parameter errors. Finally, it will be determined whether there is a line overload or near transmission limit; if so, power flow will be adjusted by increasing the power output at the sending end. Generator units with high sensitivity to bottleneck sections within the region will be identified, and their active power output at the sending end will be increased proportionally to the current output, ensuring that the unit output does not exceed the limit. If adjusting the generator output at the sending end still cannot improve the power flow convergence, load reduction control will be implemented in the receiving end area of ​​the line. Priority will be given to cutting off the non-critical loads that have the most significant impact on the bottleneck section within this area, reducing the power demand at the receiving end.

[0086] After all adjustments are completed, the system will synchronously transmit the corrected generator output, load level, and verified branch parameters back to the initial values ​​for power flow calculation. The system will then restart. P-θ Iteratively, power flow calculations are performed on the adjusted data baseline. If the calculations meet the convergence accuracy requirements... ε p Then proceed sequentially to the next step. Q-U The process involves iterative steps and reactive power adjustment. If convergence is still not achieved, the number of steps in the first-level BFS search is increased, and the above diagnostic logic is iterated until a set of physically feasible power flow data for the entire network is obtained. Finally, the converged power flow data with improved active power flow convergence is output.

[0087] Step S32, based on Q-U Iterative reactive power flow adjustment, such as Figure 5 As shown.

[0088] Step S321: After completing the first stage of active power flow convergence and locking the global phase angle distribution, the second stage of phased decoupling to improve power flow convergence is performed. The solver parameters 'Decoupled P Solution Logic' are set to False and 'Decoupled Q Solution Logic' to True. In this stage, the voltage phase angle remains constant during the iteration process; only the voltage amplitude is iteratively solved for the reactive power balance equation. The calculation steps are as follows: a) Substitute into the first stage P-θ The system power flow parameter set obtained after iterative convergence, keeping the voltage phase angle constant during the iteration process, is set to the number of iterations. m = 0.

[0089] b) Substitute into Equation 23 to calculate the reactive power imbalance at each node. and find .

[0090] c) Solve Equation 21 to obtain the voltage amplitude correction for each node. .

[0091] d) Correct the magnitude of the voltage at each node. : (27) e) Determine the number of iterations m Has the maximum number of iterations been reached? If so, proceed to step g; otherwise, proceed to step f.

[0092] f) Determine the reactive power imbalance at each node. Does it meet the preset convergence accuracy? ε q If the condition is met, proceed to step g; otherwise, set the iteration count to zero.m = m +1, and return to step b for iteration.

[0093] g) Output the final power flow data set Step S322: Locating the power flow problem of non-converging reactive power.

[0094] when Q-U Iterations within a preset maximum number of iterations Within, its maximum reactive power imbalance Δ Q Still unable to meet the preset convergence accuracy This indicates a severe imbalance between the current reactive power resource allocation and the voltage control target, causing the computational model to be unable to obtain a feasible solution under the current physical constraints.

[0095] When the reactive power flow calculation fails to converge, apply a correction based on the maximum voltage amplitude. Problem localization and regional adjustment strategies. Reactive power balance is highly coupled with node voltage amplitude. During the iteration process, reactive power imbalance... It only reflects the reactive power difference at the node, while the voltage amplitude correction is not included. It is the response of the reactive power difference after being mapped by the Jacobian matrix, reflecting the voltage adjustment range required by the system to eliminate the reactive power difference. The larger the value, the higher the reactive power regulation capability required for the node to maintain the target voltage, and the closer it is to its reactive power regulation limit. Compared to , This approach more fundamentally reflects the root cause of the system's reactive voltage non-convergence. A breadth-first search is then performed starting from this bus to delineate the critical regions associated with it.

[0096] In each Q-U During iterative calculations, the system records the reactive power imbalance of all nodes in the network in real time. With voltage amplitude correction Simultaneously extract the maximum voltage correction amount within the current iteration cycle. and its corresponding bus nodes. Focus on observing the period approaching the maximum iteration count. Bus node. If in multiple consecutive iterations, If a node is always located on a specific bus, it indicates that the voltage balance capability at that node has reached its physical limit, and the bus can be identified as a weak point in the power grid that causes reactive power flow divergence.

[0097] Step S323: Strategy to improve the convergence of power flow when reactive power fails to converge.

[0098] When reactive power flow fails to converge, a finite-level BFS topology search program is first run based on the magnitude of the correction, typically searching a second or third-level neighborhood. This bus voltage correction is then used to determine the final bus voltage level. The sign of the reactive power is used to determine whether the busbar is in a state of insufficient or excessive reactive power. For the type of insufficient reactive power, first increase the reactive power output of the Static Var Compensator (SVC) and Static Var Generator (SVG) with the fastest response speed in the area. If the reactive power deficit exceeds its dynamic response capability, further adjust the discrete reactive power equipment, connect parallel capacitors, and disconnect the parallel reactors in the area. If it is still not balanced, further disconnect the inductive load. For the type of excessive reactive power, first adjust the Static Var Compensator (SVC) and Static Var Generator (SVG) with the fastest response speed to reduce their reactive power output. If the reactive power deficit exceeds its dynamic response capability, connect parallel reactors and disconnect the parallel capacitors in the area. If the reactive power non-convergence is caused by the reciprocating jump of the DC converter transformer tap changer during the iteration process, the tap changer step size will be appropriately increased to avoid the algorithm triggering continuous step adjustment of discrete parameters under small reactive power fluctuations, thereby cutting off the numerical oscillation source and stabilizing the reactive power in the area.

[0099] After the reactive power flow convergence improvement is completed, the system will synchronously transmit the corrected grid parameters back to the system. Q-U The initial values ​​for the iterative calculation. The system restarts. Q-U Iteratively, power flow calculations are performed on the adjusted data baseline. If the calculations meet the convergence accuracy requirements... If convergence is still not achieved, the subsequent fully coupled power flow calculation will proceed sequentially. If convergence is still not achieved, the number of steps in the first-level BFS search will be increased and the above diagnostic logic will be iterated until a set of physically feasible power flow data for the entire network is obtained, and finally, the complete power flow data after the reactive power flow convergence is improved will be output.

[0100] Step S4: Based on the obtained power flow parameters as initial values, perform fully coupled power flow calculations to obtain the final power flow solution that meets the convergence accuracy requirements.

[0101] After completion P-θ and Q-U After two phases of decoupling and adjustment, the system has obtained high-precision power flow parameters. While this decoupling is efficient, it deliberately ignores... P-U and Q-θ The cross-coupling between them means that the result can only be considered as a high-precision approximate solution.

[0102] Therefore, a fully coupled power flow verification is introduced in the final stage. Since the preceding decoupling adjustments have pushed the system power flow parameters into the minimum neighborhood of the final solution, forming extremely high-quality initial values ​​for iteration, the verification procedure can flexibly use the Newton-Raphson method (NR method) or the PQ decomposition method for global closed-loop solution. If the PQ decomposition method is used, both the 'Decoupled Power Logic' and 'Decoupled Q Solution Logic' parameters are set to True.

[0103] The calculation steps are as follows: 1) Load the system power flow parameter set obtained after phased decoupling power flow calculation, and set the iteration number. k =0.

[0104] 2) Calculate the first according to equation 22-23. k Active power imbalance at each node in the next iteration and reactive power imbalance .

[0105] 3) If the NR method is used, the complete Jacobian matrix is ​​constructed according to Equations 11-18 to solve Equation 10 and obtain the phase angle correction for each node. and voltage amplitude correction amount If the PQ decomposition method is used, a coefficient matrix is ​​constructed, and the coefficients of each node are solved using formulas 20-21. as well as .

[0106] 4) Correct the amplitude of the voltage at each node. : (28) 5) Correct the phase angle of the voltage at each node. : (29) 6) Determine the active power imbalance at each node. and reactive power imbalance Check if the convergence condition is met. If it is met, proceed to step 7); if not, set the iteration count to zero. k = k +1, return to step 2) and iterate.

[0107] 7) Calculate the power of the slack node and the total power of all lines. Because the aforementioned phased adjustments provide high-quality approximate initial values, the fully coupled power flow calculation stage requires only a very small number of iterations to quickly and accurately locate a feasible power flow solution that meets the system's operational constraints and convergence accuracy standards. For example... Figure 6 As shown, through the power flow convergence improvement scheme proposed in this invention, non-convergent power flow conditions are gradually guided back to the effective convergence domain.

[0108] Example 2 This embodiment discloses a power flow convergence system based on DC power flow preprocessing and staged decoupling; like Figure 7 As shown, the power flow convergence system based on DC power flow preprocessing and staged decoupling includes: The data verification module is configured to perform multi-level verification and parameter adjustment on the raw power flow data, and to perform preliminary verification of the system active power balance to generate a set of initial values ​​for node active power settings that meet basic physical feasibility. The DC power flow preprocessing module is configured to: perform DC power flow calculation based on the set of initial values ​​of the active power of the nodes, obtain the initial phase angle of each node, identify the heavily loaded lines whose power angle difference exceeds a preset threshold, adjust the local active power resources to make the power angle difference of all lines meet the constraints, and obtain the initial phase angle distribution. The phased decoupled power flow calculation and adjustment module is configured to: combine the preset initial voltage amplitude with the initial phase angle distribution to perform phased decoupled power flow calculation and adjustment to obtain power flow parameters; The fully coupled power flow verification module is configured to perform fully coupled power flow calculations based on the obtained power flow parameters as initial values, and obtain the final power flow solution that meets the convergence accuracy requirements.

[0109] Example 3 The purpose of this embodiment is to provide a computer-readable storage medium.

[0110] A computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps in the power flow convergence method based on DC power flow preprocessing and phased decoupling as described in Example 1.

[0111] Example 4 The purpose of this embodiment is to provide an electronic device.

[0112] An electronic device includes a memory, a processor, and a program stored in the memory and executable on the processor. When the processor executes the program, it implements the steps in the power flow convergence method based on DC power flow preprocessing and phased decoupling as described in Embodiment 1.

[0113] The steps and methods involved in the apparatuses of Embodiments 2, 3, and 4 above correspond to those in Embodiment 1. For specific implementation details, please refer to the relevant description section of Embodiment 1. The term "computer-readable storage medium" should be understood as a single medium or multiple media including one or more instruction sets; it should also be understood as including any medium capable of storing, encoding, or carrying an instruction set for execution by a processor and enabling the processor to perform any of the methods in this invention.

[0114] Those skilled in the art will understand that the modules or steps of the present invention described above can be implemented using general-purpose computer devices. Optionally, they can be implemented using computer-executable program code, thereby allowing them to be stored in a storage device for execution by a computer device, or they can be fabricated as separate integrated circuit modules, or multiple modules or steps can be fabricated as a single integrated circuit module. The present invention is not limited to any particular combination of hardware and software.

[0115] While the specific embodiments of the present invention have been described above in conjunction with the accompanying drawings, this is not intended to limit the scope of protection of the present invention. Those skilled in the art should understand that various modifications or variations that can be made by those skilled in the art without creative effort based on the technical solutions of the present invention are still within the scope of protection of the present invention.

Claims

1. A power flow convergence method based on DC power flow preprocessing and decoupled power flow in stages, characterized in that, include: The original power flow data is subjected to multi-level verification and parameter adjustment, and the initial verification of the system active power balance is performed to generate a set of initial values ​​for node active power settings that meet basic physical feasibility. DC power flow calculation is performed based on the set of initial values ​​of the active power of the nodes, the initial phase angle of each node is obtained, the heavily loaded lines whose power angle difference exceeds the preset threshold are identified, and the power angle difference of all lines is satisfied by adjusting the local active power resources to obtain the initial phase angle distribution. By combining the preset initial voltage amplitude with the initial phase angle distribution, staged decoupled power flow calculation and adjustment are performed to obtain power flow parameters; Using the obtained power flow parameters as initial values, a fully coupled power flow calculation is performed to obtain the final power flow solution that meets the convergence accuracy requirements.

2. The DC power flow pre-processing and phase-decoupled power flow convergence method of claim 1, wherein, The process involves multi-level verification and parameter adjustment of the original power flow data, and preliminary verification of the system's active power balance, generating a set of initial node active power settings that meet basic physical feasibility, including: The original power flow data is subjected to network topology verification and device parameter verification. The network topology verification includes identifying and removing isolated buses, verifying the consistency of balancing machine configuration and circuit breaker status. The device parameter verification includes verifying the non-negativity and rationality of generator, transformer, and line parameters, as well as the rationality of node power and voltage setting values. The initial total active power imbalance of the system is calculated, and it is determined whether the power to be borne by the balancing machine exceeds its preset upper and lower power limits. If the power of the balancing machine exceeds the limit, the excess power is distributed proportionally to other controllable generators in the system other than the balancing machine, so as to obtain the initial set of node active power settings that meet the basic physical feasibility.

3. The DC power flow pre-processing and phase-decoupled power flow convergence method of claim 1, wherein, The process of calculating DC power flow based on the initial set of active power values ​​at the nodes to obtain the initial phase angle of each node includes: Substituting the active power injection power of all nodes in the initial set of active power settings into the DC power flow equation matrix, as shown below: wherein, is the active injection of all non-balance nodes; is the set of unknown phase angles of all non-balance nodes; is the admittance matrix; Solving the linear equation yields the initial phase angles of each non-equilibrium node. Setting the phase angles of the equilibrium nodes to 0 gives the initial phase angles of all nodes in the network.

4. The DC power flow pre-processing and phase-decoupled power flow convergence method of claim 1, wherein, The identification of overloaded lines with a power angle difference exceeding a preset threshold includes: Based on the DC power flow calculation results, calculate the absolute value of the power angle difference between the two ends of all lines; The absolute value of the power angle difference is compared with a preset static stability safety threshold to identify lines that exceed the limit; Using the sending and receiving busbars of each over-limit line as the center, a local topology adjustment area is delineated using a breadth-first search algorithm; According to the preset priority order, the distribution of active resources in the local topology area is adjusted, and the active power output of the opposite area is adjusted synchronously to maintain the overall active power balance of the system. Recalculate the DC power flow until the absolute value of the power angle difference of all lines is less than the static stability safety threshold.

5. The DC power flow pre-processing and phase-decoupled power flow convergence method of claim 1, wherein, The phased decoupling power flow calculation and adjustment includes: First stage: fixed voltage amplitude, only P-θ iterative solution of active power flow, if the iteration does not converge, the weak bus is located based on the maximum phase angle correction amount, and the topological search and active power resource adjustment are performed on its associated region until P-θ iteration converges; Second stage: On the basis of locking the phase angle distribution obtained in the first stage, the voltage phase angle is fixed, and only QU iterative solution of reactive power flow; if the iteration does not converge, the weak bus is located based on the maximum voltage correction amount, and the topological search and reactive power resource adjustment are performed on the associated region until QU iteration converges.

6. The DC power flow pre-processing and phase-decoupled power flow convergence method of claim 5, wherein, The weak bus is positioned based on the maximum phase angle correction amount, and a topological search and active resource adjustment are performed on the associated region until P-θ iterative convergence, comprising: When P-θ When the iteration reaches a preset maximum number of times and still does not converge, the maximum phase angle correction amount in the iteration period leading to non-convergence and its corresponding bus node are extracted; Starting from the aforementioned busbar, a breadth-first search with a finite number of steps is used to determine the key adjustment area; Perform at least one of the following operations on the critical adjustment area: verify and correct line and transformer parameters, adjust the active power output of generators in the area according to the principle of sensitivity or proportionality, and reduce non-critical loads in the area.

7. The DC power flow pre-processing and phase-decoupled power flow convergence method of claim 5, wherein, The weak bus is positioned based on the maximum voltage correction amount, and a topology search and reactive power resource adjustment are performed on the associated area until QU Iterative convergence includes: When QU When the iteration reaches a preset maximum number of times and still does not converge, the maximum voltage correction amount in the iteration period leading to non-convergence and its corresponding bus node are extracted; A key adjustment region is determined by breadth-first search with a limited number of steps, starting from the bus; The reactive state of the region is determined according to the positive and negative nature of the maximum voltage correction amount; if the maximum voltage correction amount is positive, capacitive reactive power is put in or inductive reactive power equipment is cut off according to the priority; if the maximum voltage correction amount is negative, inductive reactive power is put in or capacitive reactive power equipment is cut off according to the priority.

8. A power flow convergence system based on DC power flow preprocessing and phase-wise decoupling, characterized in that, It comprises: A data verification module configured to perform multi-level verification and parameter adjustment on original power flow data, and perform preliminary checking of system active balance to generate a set of initial values of node active power settings meeting basic physical feasibility; A DC power flow preprocessing module configured to perform DC power flow calculation based on the set of initial values of node active power settings, to obtain initial phase angles of each node, to identify overloaded lines with power angle differences exceeding a preset threshold, and to make all line power angle differences meet the constraints through local active resource adjustment to obtain an initial phase angle distribution; A phased decoupling power flow calculation and adjustment module configured to perform phased decoupling power flow calculation and adjustment in combination with a preset initial voltage amplitude and the initial phase angle distribution to obtain power flow parameters; A full-coupling power flow verification module configured to perform full-coupling power flow calculation based on the obtained power flow parameters as initial values to obtain a final power flow solution meeting convergence accuracy requirements.

9. A computer-readable storage medium having stored thereon a program, characterized in that, The program is executed by the processor to implement the steps in the power flow convergence method based on DC power flow preprocessing and phased decoupling according to any one of claims 1-7.

10. An electronic device comprising a memory, a processor, and a program stored on the memory and executable on the processor, characterized in that, The processor executes the program to implement the steps in the power flow convergence method based on DC power flow preprocessing and phased decoupling according to any one of claims 1-7.

Images