POWER FLOW CALCULATION METHOD AND DEVICE FOR AC-DC INTERCONNECTED POWER SYSTEM, STORAGE MEDIUM AND TERMINAL

Abstract

Provided are a power flow calculation method and device for an alternating current (AC)-direct current (DC) interconnected power system, a storage medium and terminal. A state of a DC network is calculated according to a control mode of a converter station, and a connection point between the DC network and an AC network is equivalent to a power node; and power flow calculation is performed through a Newton-Raphson method.

Claims

1. A power flow calculation method for an alternating current (AC)-direct current (DC) interconnected system, comprising: solving a conductance matrix for a DC network of the AC-DC interconnected power system, and acquiring a resistance between any two of converters in the DC network, or acquiring a resistance between any two of connection points of the DC network of hierarchical structures; acquiring a DC voltage and an active power of a node corresponding to each of the converters according to a structure of the DC network; acquiring a control mode of the each of the converters; calculating a reactive power injection amount of the each of the converters into an AC power grid according to the DC voltage, the active power, and the control mode of the each of the converters; and performing power flow calculation through a Newton-Raphson method according to the resistance between the any two of the converters or the resistance between the any two of the connection points, the DC voltage, the active power, the reactive power injection amount and the control mode of the each of the converters.

2. The method of claim 1, wherein the acquiring the DC voltage and the active power of the node corresponding to the each of the converters according to the structure of the DC network comprises: acquiring a node parameter corresponding to the each of the converters according to the structure of the DC network; constructing, according to the node parameter, an equation set: $\{\begin{array}{c}{I}_{\mathrm{dk}}={.Math.}_{j=1}^{n_c}.Math.G_{\mathrm{kj}}.Math.V_{\mathrm{dj}}\\ P_{\mathrm{dk}}=I_{\mathrm{dk}}.Math.V_{\mathrm{dk}}\end{array},$ wherein V.sub.dk is the DC voltage of the node corresponding to the each of the converters, P.sub.dk is the active power of the node corresponding to the each of the converters, and I.sub.dk is a DC current flowing into a converter station k, G.sub.kj is an admittance matrix element between a node k corresponding to the converter station k and a node j, V.sub.dj is a voltage of a DC bus connected to a converter j, and n.sub.c is a number of converters in the DC network; and calculating the DC voltage and the active power of the node corresponding to the each of the converters according to the equation set.

3. The method of claim 2, wherein in condition that the DC network of the AC-DC interconnected power system comprises a hierarchical structure, active power outputted by a single converter on a series side of the DC network is proportional to a voltage ratio of the single converter.

4. The method of claim 3, wherein in condition that the DC network of the AC-DC interconnected power system comprises the hierarchical structure, the following relationship is satisfied: $\{\begin{array}{c}{I}_{\mathrm{di}.Math..Math.1}=I_{\mathrm{di}.Math..Math.2}=I_d\\ I_d=\frac{V_{\mathrm{dr}}-\left(V_{\mathrm{di}.Math..Math.1}+V_{\mathrm{di}.Math..Math.2}\right)}{R_d}\end{array};$ wherein I.sub.di1 is a current flowing through a high-voltage converter of a converter station in the hierarchical structure and I.sub.di2 is a current flowing through a low-voltage converter of the converter station in the hierarchical structure, I.sub.d denotes a current flowing through the converter station, V.sub.dr is a sending end voltage of the DC network, V.sub.di1 denotes a DC voltage of the high-voltage converter in the hierarchical structure and V.sub.di2 denotes a DC voltage of the low-voltage converter in the hierarchical structure, and R.sub.d is a resistance of a DC line; and wherein the calculating the DC voltage and the active power of the node corresponding to the each of the converters according to the equation set specifically comprises: calculating the DC voltage and the active power according to the following equation set: $\{\begin{array}{c}{P}_{\mathrm{idk}}=k_{\mathrm{idk}}.Math.P_d\\ V_{\mathrm{idk}}=k_{\mathrm{idk}}.Math.V_d\end{array};$ wherein k.sub.idk is a voltage ratio of a converter k in the hierarchical structure, P.sub.d is active power of the DC network injected into a converter station, and V.sub.d is a DC voltage of a node connected to the converter station, P.sub.idk is active power outputted by the converter k in the hierarchical structure, and V.sub.idk is a DC voltage applied across the converter k in the hierarchical structure.

5. The method of claim 1, wherein the control mode of the each of the converters comprises a first type control mode and a second type control mode; wherein the first type control mode comprises constant active power control mode, a constant DC voltage control mode, and a constant DC current control mode; and the second type control mode comprises a constant transformation ratio control mode and a constant overlap angle control mode.

6. The method of claim 5, wherein in condition that the control mode of the each of the converters is the first type control mode, the calculating the reactive power injection amount of the each of the converters into the AC power grid according to the DC voltage, the active power, and the control mode of the each of the converters comprises calculating the reactive power injection amount according to the following equation set: $\{\begin{array}{c}{P}_{\mathrm{dk}}=I_{\mathrm{dk}}.Math.V_{\mathrm{dk}}\\ Q_{\mathrm{dk}}=V_{\mathrm{dk}}.Math.I_{\mathrm{dk}}.Math.\tan .Math..Math.{}_k\end{array},$ wherein I.sub.dk is a DC current flowing into a converter k, P.sub.dk is the active power, V.sub.dk is the DC voltage, .sub.k is a power factor of the converter k, and Q.sub.dk is the reactive power injection amount.

7. The method of claim 5, wherein in condition that the control mode of the each of the converters is the second type control mode and the second type control mode is the constant overlap angle control mode, the calculating the reactive power injection amount of the each of the converters into the AC power grid according to the DC voltage, the active power, and the control mode of the each of the converters specifically comprises calculating the reactive power injection amount according to the following equation: $Q_{\mathrm{dk}}=\frac{k_y.Math.P_{\mathrm{dk}}\left(V_{\mathrm{dk}}^2+P_{\mathrm{idk}}.Math.X_c\right).Math.\sqrt{1-\frac{V_{\mathrm{dk}}^4.Math.{\cos }^2\left({}_d\right)}{{k_y^2\left(V_{\mathrm{dk}}^2+P_{\mathrm{idk}}.Math.X_c\right)}^2}}}{V_{\mathrm{dk}}^2},$ wherein V.sub.dk is a DC power transmission voltage, P.sub.dk is active power flowing into a converter k, P.sub.idk is active power of the DC network injected into an AC node i, .sub.d is a control angle of the converter, and X.sub.c is an overlap resistance, k.sub.y is a converter constant, and Q.sub.dk is the reactive power injection amount.

8. The method of claim 5, wherein in condition that the control mode of the each of the converters is the second type control mode and the second type control mode is the constant transformation ratio mode, the calculating the reactive power injection amount of the each of the converters into the AC power grid according to the DC voltage, the active power, and the control mode of the each of the converters specifically comprises calculating the reactive power injection amount according to the following equation: $Q_{\mathrm{dk}}=\frac{P_{\mathrm{dk}}.Math.\sqrt{-V_{\mathrm{dk}}+k_y^2.Math.k_T^2.Math.V_a^2}}{V_{\mathrm{dk}}}$ wherein V.sub.dk is a DC power transmission voltage, P.sub.dk is a active power flowing into a converter k, V.sub.a is a voltage amplitude of a node connected to the converter, k.sub.T is a transformation ratio, and k.sub.y is a converter constant.

9. The method of claim 5, wherein the performing the power flow calculation through the Newton-Raphson method according to the resistance between the any two of the converters or the resistance between the any two of the connection points, the DC voltage, the active power, the reactive power injection amount and the control mode of the each of the converters comprises: acquiring an unbalance amount of the active power and an unbalance amount of the reactive power injection amount in the power flow calculation according to the resistance between the any two of the converters or the resistance between the any two of the connection points, the DC voltage, the active power, the reactive power injection amount and the control mode of the each of the converters; establishing a Jacobian matrix for the power flow calculation according to the control mode of the each of the converters, the unbalance amount of the active power and the unbalance amount of the reactive power injection amount, wherein in condition that the control mode of the each of the converters is the second type control mode and the second type control mode is the constant overlap angle control mode, a Jacobian matrix parameter of the node corresponding to the each of the converters is determined merely by an AC network parameter, and in condition that the control mode of the each of the converters is the second type control mode and the second type control mode is the constant transformation ratio control mode, the Jacobian matrix parameter of the node corresponding to the each of the converters is corrected after being determined by the AC network parameter; and performing the power flow calculation through the Newton-Raphson method according to the Jacobian matrix.

10. The method of claim 9, wherein in condition that the control mode of the each of the converters is the second type control mode and the second type control mode is the constant transformation ratio control mode, the performing the power flow calculation through the Newton-Raphson method according to the Jacobian matrix further comprises: correcting an element Lii of the Jacobian matrix as follows: $L_{\mathrm{ii}}=-V_i.Math.\underset{ji,ji}{.Math.}.Math.V_j\left(G_{\mathrm{ij}}.Math.\sin .Math..Math.{}_{\mathrm{ij}}-B_{\mathrm{ij}}.Math.\cos .Math..Math.{}_{\mathrm{ij}}\right)+2.Math.V_i^2.Math.B_{\mathrm{ii}}-\frac{k_y^2.Math.k_T^2.Math.P_{\mathrm{dk}}.Math.V_a}{V_{\mathrm{dk}}.Math.\sqrt{-V_{\mathrm{dk}}^2+k_y^2.Math.k_T^2.Math.V_i^2}}$ wherein i is a node of an AC network connected to the each of the converters, V.sub.i is a voltage amplitude of the node i, G.sub.ij and B.sub.ij are respectively a real part and an imaginary part of an admittance matrix, V.sub.a is a voltage amplitude of a node connected to the each of the converters, V.sub.dk is a DC power transmission voltage, P.sub.dk is active power flowing into a converter k, k.sub.T is a transformation ratio, k.sub.y is a converter constant, .sub.ij is a control angle of the node i, H, N and L are block matrices of the Jacobian matrix, P is the unbalance amount of the active power, Q is the unbalance amount of the reactive power injection amount, and and V are correction amounts of variables in an iterative process.

11. The method of claim 10, wherein the performing the power flow calculation through the Newton-Raphson method further comprises: determining whether a calculation result of the power flow satisfies a convergence condition; in condition that the calculation result of the power flow satisfies the convergence condition, completing the power flow calculation; and in condition that the calculation result of the power flow does not satisfy the convergence condition, acquiring the unbalance amount of the active power and the unbalance amount of the reactive power injection amount in the power flow calculation again.

12. A terminal, comprising a display screen, a memory, a processor and computer programs stored in the memory and executable by the processor, wherein the processor is configured to, when executing the computer programs which, implement the power flow calculation method for the alternating current (AC)-direct current (DC) interconnected power system of any one of claim 1.

13. A power flow calculation method for an alternating current (AC)-direct current (DC) interconnected power system having DC hierarchical structures, comprising: solving a conductance matrix for a DC network, and acquiring a resistance between any two of converters, or acquiring a resistance between any two of connection points of a DC power grid of hierarchical structures; analyzing a control mode of each of the converters, and calculating a voltage and active power of each node; calculating a reactive power injection amount into an AC power grid according to the control mode, and the voltage and the active power of the each node; and performing power flow calculation through a Newton-Raphson method to obtain a calculation result.

14. The method of claim 13, comprising: determining part of parameters of the each node according to the control mode of the each of the converters, and constructing an equation set: $\begin{array}{c}{I}_{\mathrm{dk}}=\underset{j=1}{\overset{n_c}{.Math.}}.Math.G_{\mathrm{kj}}.Math.V_{\mathrm{dj}}\\ P_{\mathrm{dk}}=I_{\mathrm{dk}}.Math.V_{\mathrm{dk}}\end{array}\},$ wherein I.sub.dk is a DC current flowing into a converter station k, and G.sub.kj is an admittance matrix element between a node k corresponding to the converter station k and a node j; and calculating the voltage and the active power of the each node according to a specific equation set: $\begin{array}{c}{I}_{\mathrm{di}.Math..Math.1}=I_{\mathrm{di}.Math..Math.2}=I_d\\ \frac{V_{\mathrm{dr}}-\left(V_{\mathrm{di}.Math..Math.1}+V_{\mathrm{di}.Math..Math.2}\right)}{R_d}=I_d\end{array}\},$ wherein I.sub.di1 is a current flowing through a high-voltage converter of the converter station in the hierarchical structures and I.sub.di2 is a current flowing through a low-voltage converter of the converter station in the hierarchical structures; I.sub.d denotes a current flowing through the converter station, V.sub.dr is a sending end voltage of the DC network; V.sub.di1 denotes a DC voltage of the high-voltage converter in the hierarchical structures and V.sub.di2 denotes a DC voltage of the low-voltage converter in the hierarchical structures; and R.sub.d is a resistance of a DC line.

15. The method of claim 13, wherein an active power outputted by a single converter on a series side is proportional to a voltage ratio of the single converter.

16. The method of claim 13, comprising: calculating the reactive power injection amount into the AC power grid according to a control mode of a converter at a node, which is a control mode for a constant overlap angle, by using the following equation: $Q_{\mathrm{idc}}=\frac{k_y.Math.P_{\mathrm{idc}}\left(V_{d.Math..Math.c}^2+P_{\mathrm{idc}}.Math.X_c\right)\sqrt{1-\frac{V_{\mathrm{dk}}^4.Math.{\cos }^2\left({}_d\right)}{{k_y^2\left(V_{\mathrm{dk}}^2+P_{\mathrm{idc}}.Math.X_c\right)}^2}}}{V_{\mathrm{dk}}^2},$ calculating the reactive power injection amount into the AC power grid according to the control mode of the converter at the node, which is a control mode for a constant transformation ratio, by using the following equation: $Q_{d.Math..Math.c}=\frac{P_{d.Math..Math.c.Math.}.Math.\sqrt{-V_{d.Math..Math.c}+k_y^2.Math.k_T^2.Math.V_a^2}}{V_{d.Math..Math.c}},$ and calculating an active power injection amount by calculating the state of the DC network; wherein V.sub.dc is a voltage of a DC network node connected to the converter station, .sub.d is a control angle of the converter which comprises a gating delay angle of a rectifier and an extinction advance angle of an inverter; k.sub.T is a transformation ratio; X.sub.c is an overlap resistance, a variable k.sub.y is introduced considering an effect of an overlap angle, .sub.i is a power factor angle corresponding to active power and reactive power absorbed by the converter from an AC system; and V.sub.a is a voltage amplitude of an AC network connected to the converter.

17. The method of claim 13, comprising: when the hierarchical structures are involved, completing equivalence of the power node according to a voltage ratio of a converter in the hierarchical structures.

18. The method of claim 13, wherein a specific process for calculating the reactive power injection amount into the AC power grid comprises: in condition that a control mode of a converter corresponding to a node is a constant overlap angle control mode, calculating the reactive power injection amount by using the following equation: $Q_{\mathrm{idc}}=\frac{k_y.Math.P_{\mathrm{idc}}\left(V_{d.Math..Math.c}^2+P_{\mathrm{idc}}.Math.X_c\right).Math.\sqrt{1-\frac{V_{\mathrm{dk}}^4.Math.{\cos }^2\left({}_d\right)}{{k_y^2\left(V_{d.Math..Math.c}^2+P_{\mathrm{idc}}.Math.X_c\right)}^2}}.Math.\sec \left({}_d\right)}{V_{d.Math..Math.c}^2},$ wherein V.sub.dc is a voltage of a DC network node connected to the converter station, .sub.d is a control angle of the converter which comprises a gating delay angle of a rectifier and an extinction advance angle of an inverter, k.sub.t is a transformation ratio, X.sub.c is an overlap resistance, a variable k.sub.y is introduced considering an effect of an overlap angle, and .sub.i is a power factor angle corresponding to active power and reactive power absorbed by the converter from an AC system; or in condition that a control mode of a converter corresponding to a node is a constant transformation ratio control mode, calculating the reactive power injection amount by using the following equation: $Q_{d.Math..Math.c}=\frac{P_{d.Math..Math.c}.Math.\sqrt{-V_{d.Math..Math.c}+k_y^2.Math.k_T^2.Math.V_a^2}}{V_{d.Math..Math.c}}$ and calculating a derivative of the reactive power injection amount with respect to an AC voltage corresponding to the reactive power injection amount; wherein V.sub.a is a voltage amplitude of an AC network connected to the converter; and when the hierarchical structures are involved, calculating a power effect of each of the hierarchical structures on a connection point with the AC power grid according to the following equation set: $\{\begin{array}{c}{P}_{\mathrm{idk}}=k_{\mathrm{idk}}.Math.P_d\\ V_{\mathrm{idk}}=k_{\mathrm{idk}}.Math.V_d\end{array}$ wherein k.sub.idk is a voltage ratio of a converter k in the hierarchical structures, P.sub.d is active power of the DC network injected into the converter station, and V.sub.d is a DC voltage of a node connected to the converter station, P.sub.idk is active power outputted by the converter k in the hierarchical structures, and V.sub.idk is a DC voltage applied across the converter k in the hierarchical structures.

19. The method of claim 13, wherein a specific process for performing the power flow calculation through the Newton-Raphson method comprises: setting an initial value of the AC network and calculating an unbalance amount of a power flow power equation; constructing a Jacobian matrix, in condition that the control mode of the each of the converters is a constant overlap angle control mode, a Jacobian matrix parameter of a node corresponding to the each of the converters is determined merely by an AC network parameter; and in condition that the control mode of the each of the converters is a constant transformation ratio control mode, the Jacobian matrix parameter of the node corresponding to the each of the converters is corrected after being determined by the AC network parameter; and correcting the AC network parameter, determining whether the calculation result of the power flow satisfies a convergence condition, and in condition that the calculation result of the power flow satisfies the convergence condition, ending an iteration; in condition that the calculation result of the power flow is not satisfied with the convergence condition, repeating the specific process

20. The method of claim 19, wherein a correction mode is: $L_{\mathrm{ii}}=-V_i.Math.\underset{ji,ji}{.Math.}.Math.V_j\left(G_{\mathrm{ij}}.Math.\sin .Math..Math.{}_{\mathrm{ij}}-B_{\mathrm{ij}}.Math.\cos .Math..Math.{}_{\mathrm{ij}}\right)+2.Math.V_i^2.Math.B_{\mathrm{ii}}-\frac{k_y^2.Math.k_T^2.Math.P_{d.Math..Math.c}.Math.V_a}{V_{d.Math..Math.c}.Math.\sqrt{-V_{d.Math..Math.c}^2+k_y^2.Math.k_T^2.Math.V_i^2}}$ wherein V.sub.i ji, jiV.sub.j(G.sub.ij sin .sub.ijB.sub.ij cos .sub.ij)+2V.sub.i.sup.2 B.sub.ii is a calculation formula of a Jacobian matrix element L in traditional pure AC power flow calculation, V.sub.i is a voltage amplitude of a node i; G.sub.ij and B.sub.ij are respectively a real part and an imaginary part of an admittance matrix; and V.sub.a is a voltage amplitude of a node connected to the each of the converters and is numerically consistent with V.sub.i.

Description

BRIEF DESCRIPTION OF DRAWINGS

[0078] FIG. 1 is a flowchart of a power flow calculation method for an AC-DC interconnected power system according to an embodiment of the present disclosure;

[0079] FIG. 2 is a schematic diagram illustrating a circuit structure of a DC network of hierarchical structures according to an embodiment of the present disclosure;

[0080] FIG. 3 is a flowchart of a method for acquiring a DC voltage and active power according to an embodiment of the present disclosure;

[0081] FIG. 4 is a flowchart of a method for preforming power flow calculation through a Newton-Raphson method according to an embodiment of the present disclosure;

[0082] FIG. 5 is a flowchart of another method for preforming power flow calculation through a Newton-Raphson method according to an embodiment of the present disclosure;

[0083] FIG. 6 is a block diagram of a power flow calculation device for an AC-DC interconnected power system according to an embodiment of the present disclosure; and

[0084] FIG. 7 is a structural diagram of a terminal according to an embodiment of the present disclosure.

DETAILED DESCRIPTION

[0085] Hereinafter the present disclosure will be further described in detail in conjunction with the drawings and embodiments. It may be understood that the specific embodiments set forth below are intended to illustrate and not to limit the present disclosure. Additionally, it is to be noted that, for ease of description, merely part, not all, of the structures related to the present disclosure are illustrated in the drawings.

[0086] An embodiment of the present disclosure provides a power flow calculation method for an AC-DC interconnected power system, which is applicable to power flow calculation in an AC-DC interconnected power system. The AC-DC interconnected power flow calculation method according to the embodiment of the present disclosure may be executed by a power flow calculation device for the AC-DC interconnected power system which may be implemented by software and/or hardware. FIG. 1 is a flowchart of a power flow calculation method for the AC-DC interconnected power system according to an embodiment of the present disclosure. As shown in FIG. 1, the power flow calculation method for the AC-DC interconnected power system includes steps described below.

[0087] In S110, a conductance matrix is solved for a DC network of the AC-DC interconnected power system, and a resistance between any two of converters in the DC network is acquired, or a resistance between any two of connection points of the DC network of hierarchical structures is acquired.

[0088] Specifically, the AC-DC interconnected power system includes the DC network and an AC power grid. The DC network and an AC network are connected to each other through a converter. The converter may convert an AC signal in the AC power grid into a DC signal to be inputted into the DC network, and the converter may also convert the DC signal in the DC network into the AC signal to be inputted into the AC network, so that the converter of the AC-DC interconnected power system has important contribution to a stable operation of the AC-DC interconnected power system. For a multi-terminal DC network, the resistance between the any two of the converters in the DC network may be obtained by solving the conductance matrix of the DC network. Alternatively, in condition that hierarchical structures exist the DC network of the AC-DC interconnected power system, the resistance between the any two of the connection points of the DC network of the hierarchical structures may be obtained by solving the conductance matrix of the DC network.

[0089] In S120, a DC voltage and an active power of a node corresponding to each of the converters are acquired according to a structure of the DC network. Specifically, a parameter of the DC network generally includes capacitance, inductance and the like. However, when a research and an analysis are conducted on the stable operation of the AC-DC interconnected power system, merely a resistance characteristic of the DC network is generally considered, and the DC network is represented by an admittance matrix G.sub.d of nodes of the DC network:

[00012]$G_d=\left[\begin{array}{ccc}{G}_{11}& \mathrm{.Math.}& G_{1.Math..Math.n}\\ \mathrm{.Math.}& & \mathrm{.Math.}\\ G_{n.Math..Math.1}& \mathrm{.Math.}& G_{\mathrm{nn}}\end{array}\right].$

[0090] An injection current I.sub.d of a node of the DC network may be expressed as: I.sub.d=G.sub.dV.sub.d. Where I.sub.d is the injection current of the node of the DC network, and V.sub.d is an injection voltage of the node of the DC network. Accordingly, the DC voltage I.sub.dk and the active power P.sub.dk of the node corresponding to the each of the converters may be calculated by using the following equation set:

[00013]$\{\begin{array}{c}{i}_{\mathrm{dk}}={.Math.}_{j=1}^{n_c}.Math.G_{\mathrm{kj}}.Math.V_{\mathrm{dj}}\\ P_{\mathrm{dk}}=I_{\mathrm{dk}}.Math.V_{\mathrm{dk}}\end{array}.$

[0091] Where V.sub.dk is the DC voltage of the node corresponding to the each of the converters, P.sub.dk is the active power of the node corresponding to the each of the converters, and I.sub.dk is a DC current flowing into a converter station k, G.sub.kj is an admittance matrix element between a node k corresponding to the converter station k and a node j, V.sub.dj is a voltage of a DC bus connected to a converter j, and n.sub.c is a number of converters in the DC network.

[0092] In addition, when the DC network of the AC-DC interconnected power system includes the hierarchical structure, active power outputted by a single converter on a series side of the DC network is proportional to a voltage ratio of the single converter. FIG. 2 is a schematic diagram illustrating a circuit structure of a DC network of a hierarchical structure according to an embodiment of the present disclosure. As shown in FIG. 2, when the DC network of the AC-DC interconnected power system includes the hierarchical structures, the injection current I.sub.d of the node of the DC network should also satisfy the following relationship:

[00014]$\{\begin{array}{c}{I}_{\mathrm{di}.Math..Math.1}=I_{\mathrm{di}.Math..Math.2}=I_d\\ I_d=\frac{V_{\mathrm{dr}}-\left(V_{\mathrm{di}.Math..Math.1}+V_{\mathrm{di}.Math..Math.2}\right)}{R_d}.\end{array}$

[0093] Where I.sub.di1 is a current flowing through a high-voltage converter of a converter station in the hierarchical structures and I.sub.di2 is a current flowing through a low-voltage converter of the converter station in the hierarchical structures; I.sub.d denotes a current flowing through the converter station; V.sub.dr is a sending end voltage of the DC network; V.sub.di1 denotes a DC voltage of the high-voltage converter in the hierarchical structures and V.sub.di2 denotes a DC voltage of the low-voltage converter in the hierarchical structures; and R.sub.d is a resistance of a DC line. Accordingly, a DC voltage V.sub.idk and active power P.sub.idk of a node corresponding to a converter i are specifically calculated as follows:

[00015]$\{\begin{array}{c}{P}_{\mathrm{idk}}=k_{\mathrm{idk}}.Math.P_d\\ V_{\mathrm{idk}}=k_{\mathrm{idk}}.Math.V_d\end{array}.$

[0094] Where k.sub.idk is a voltage ratio of a converter k in the hierarchical structures, P.sub.d is active power of the DC network injected into a converter station, and V.sub.d is a DC voltage of a node connected to the converter station, P.sub.idk is active power outputted by the converter k in the hierarchical structures, and V.sub.idk is a DC voltage applied across the converter k in the hierarchical structures.

[0095] In S130, a control mode of the each of the converters is acquired.

[0096] Specifically, for a traditional commutated converter, each converter has two independent control variables. Assuming that a transformer tap related to the converter i may be adjusted seamlessly, a turn ratio k.sub.ti of a transformer may be linearly controlled. Therefore, active power P.sub.dci, a DC voltage V.sub.dci and a DC current I.sub.dci of a DC bus connected to the converter i may be considered as control variables in a first type control mode, which may be defined as a D-axis control mode. A transformation ratio (turn ratio) k.sub.ti of the transformer related to the converter i and a control angle .sub.i of the converter i may be considered as control variables in a second type control mode, which may be defined as an E-axis control mode. Accordingly, the control mode of the each of the converters may be divided into the first type control mode and the second type control mode.

[0097] In S140, a reactive power injection amount of the each of the converters into the AC power grid is calculated according to the DC voltage, the active power, and the control mode of the each of the converters.

[0098] Specifically, the converters in the AC-DC interconnected power system have different control modes, and the reactive power injection amounts of the converters into the AC power grid are calculated in different manners. The control mode of the each of the converters may be divided into the first type control mode and the second type control mode. The first type control mode includes a constant active power control mode, a constant DC voltage control mode, a constant DC current control mode and the like. The second type control mode includes a constant transformation ratio control mode, a constant overlap angle control mode and the like. The overlap angle here is a determined value of the control angle of the converter.

[0099] When the control mode of the each of the converters is the first type control mode, the reactive power injection amount of the each of the converters into the AC power grid is calculated according to the DC voltage, the active power, and the control mode of the each of the converters, which is specifically preformed according to the following equation:

[00016]$\{\begin{array}{c}{P}_{\mathrm{dk}}=I_{\mathrm{dk}}.Math.V_{\mathrm{dk}}\\ Q_{\mathrm{dk}}=V_{\mathrm{dk}}.Math.I_{\mathrm{dk}}.Math.\tan .Math..Math.{}_k\end{array}.$

[0100] Where I.sub.dk is a DC current flowing into the converter k, P.sub.dk is the active power, V.sub.dk is the DC voltage, .sub.k is a power factor of the converter, and Q.sub.dk is the reactive power injection amount.

[0101] When the control mode of the converter is the second type control mode, the equation set for calculating the reactive injection quantity needs to be combined with basic equations of the converter to calculate the reactive power injection amount. The basic equations of the converter are listed as follows:

[00017]$\begin{array}{c}{V}_{d^{*}}=k_{T^{*}}{V}_{s^{*}}\cos .Math..Math..Math..Math.d-X_{C^{*}}{I}_{d^{*}}\\ V_{d^{*}}=k_y{k}_{T^{*}}{V}_{s^{*}}\cos .Math..Math.\\ I_{s^{*}}=k_y{k}_{T^{*}}{I}_{d^{*}}\end{array}\}.$

[0102] Where V.sub.d* is a per unit of a DC power transmission voltage, I.sub.d* is a per unit of a DC power transmission current, V.sub.s* is a per unit of a line voltage of an AC bus, I.sub.s* is a base frequency AC current injected to the converter, k.sub.T* is the transformation ratio, .sub.d is the control angle of the converter which includes a gating delay angle of a rectifier and an extinction advance angle of an inverter, is a power factor angle corresponding to active power (absorbed by the rectifier and emitted by the inverter) and reactive power absorbed by the converter from an AC system, X.sub.c* is a per unit of an overlap resistance, and k.sub.y is a converter constant which approximates 0.995 by a simplified analysis with consideration of an effect of the overlap angle.

[0103] When the control mode of the each of the converters is the second type control mode, and the second type control mode is the constant overlap angle control mode, the reactive power injection amount of the each of the converters into the AC power grid is calculated according to the DC voltage, the active power, and the control mode of the each of the converters, and the following equation is derived from the above equations:

[00018]$Q_{\mathrm{dk}}=\frac{k_y.Math.P_{\mathrm{dk}}\left(V_{\mathrm{dk}}^2+P_{\mathrm{idk}}.Math.X_C\right)\sqrt{1-\frac{V_{\mathrm{dk}}^4.Math.{\cos }^2\left({}_d\right)}{{k_y^2\left(V_{\mathrm{dk}}^2+P_{\mathrm{idk}}.Math.X_C\right)}^2}}}{V_{\mathrm{dk}}^2}$

[0104] Where V.sub.dk is the DC power transmission voltage, P.sub.dk is active power flowing into the converter k, P.sub.idk is active power of the DC network injected into an AC node i, .sub.d is the control angle of the converter, and X.sub.c is the overlap resistance, k.sub.y is the converter constant, and Q.sub.dk is the reactive power injection amount.

[0105] When the control mode of the each of the converters is the second type control mode, and the second type control mode is the constant transformation ratio control mode, the reactive power injection amount of the each of the converters into the AC power grid is calculated according to the DC voltage, the active power, and the control mode of the each of the converters, and the following equation is derived from the above equations:

[00019]$Q_{\mathrm{dk}}=\frac{P_{\mathrm{dk}}.Math.\sqrt{-V_{\mathrm{dk}}+k_y^2.Math.k_T^2.Math.V_a^2}}{V_{\mathrm{dk}}}.$

[0106] Where V.sub.dk is the DC power transmission voltage, P.sub.dk is the active power flowing into the converter k, V.sub.a is a voltage amplitude of a node connected to the converter, k.sub.T is the transformation ratio, and k.sub.y is the converter constant.

[0107] In S150, power flow calculation is performed through a Newton-Raphson method according to the resistance between the any two of the converters or the resistance between the connection points, the DC voltage, the active power, the reactive power injection amount and the control mode of the each of the converters.

[0108] Specifically, when the existing Newton-Raphson method is used for performing the power flow calculation, an initial value is generally selected for performing iterative calculation, and a power flow calculation result is related to the selection of the initial value so that a large selection scale of the initial value leads to increased iteration times, and a Jacobian matrix is large in scale and is not advantageous to calculation. The power flow calculation in the AC-DC interconnected power system is performed through the Newton-Raphson method by acquiring the resistance between the any two of the converters or the resistance between the connection points and calculating the DC voltage, the active power and the reactive power injection amount according to the control mode of the each of the converters.

[0109] In the embodiments of the present disclosure, the corresponding parameters are obtained by analyzing the control mode of each converter in the AC-DC interconnected power system, and the power flow calculation is performed through the Newton-Raphson method, thereby avoiding the problem that the calculation is not facilitated due to a large selection scale of the initial value and a large scale of the Jacobian matrix. Therefore, the AC-DC interconnected power flow calculation has better convergence, reduced calculation complexity, an improved calculation rate, and reduced costs.

[0110] Optionally, based on the preceding embodiment, a specific method for acquiring the DC voltage and the active power is optimized. FIG. 3 is a flowchart of a method for acquiring a DC voltage and active power according to an embodiment of the present disclosure. As shown in FIG. 3, the step in which the DC voltage and the active power of the node corresponding to the each of the converters according to the structure of the DC network specifically includes steps described below.

[0111] In S131, a node parameter corresponding to the each of the converters is acquired according to the structure of the DC network.

[0112] In S132, an equation set is constructed according to the node parameter as follows:

[00020]$\{\begin{array}{c}{I}_{\mathrm{dk}}={.Math.}_{j=1}^{n_c}.Math.G_{\mathrm{kj}}.Math.V_{\mathrm{dj}}\\ P_{\mathrm{dk}}=I_{\mathrm{dk}}.Math.V_{\mathrm{dk}}\end{array}.$

[0113] Where V.sub.dk is the DC voltage of the node corresponding to the each of the converters, P.sub.dk is the active power of the node corresponding to the each of the converters, and I.sub.dk is the DC current flowing into the converter station k, G.sub.kj is the admittance matrix element between the node k corresponding to the converter station k and the node j, V.sub.dj is the voltage of the DC bus connected to the converter j, and n.sub.c is the number of converters in the DC network.

[0114] In S133, the DC voltage and the active power of the node corresponding to the each of the converters are calculated according to the equation set.

[0115] Specifically, in condition that the DC network in the AC-DC interconnected power system has different structures, the DC voltage and the active power of the node corresponding to the each of the converters are calculated in different manners. For a converter in a general DC network, an equation set may be constructed according to the node parameter to calculate the DC voltage and the active power of the node corresponding to the converter:

[00021]$\{\begin{array}{c}{I}_{\mathrm{dk}}={.Math.}_{j=1}^{n_c}.Math.G_{\mathrm{kj}}.Math.V_{\mathrm{dj}}\\ P_{\mathrm{dk}}=I_{\mathrm{dk}}.Math.V_{\mathrm{dk}}\end{array}.$

[0116] Therefore, the DC voltage and the active power are calculated according to the above equation set.

[0117] For the DC network of the hierarchical structures, the active power outputted by the single converter on the series side of the DC network is proportional to the voltage ratio of the single converter.

[0118] The DC voltage and the active power may be calculated according to the above equation set:

[00022]$\{\begin{array}{c}{P}_{\mathrm{idk}}=k_{\mathrm{idk}}.Math.P_d\\ V_{\mathrm{idk}}=k_{\mathrm{idk}}.Math.V_d\end{array}.$

[0119] Optionally, based on the preceding embodiment, a method for preforming the power flow calculation through the Newton-Raphson method is optimized. FIG. 4 is a flowchart of a method for preforming power flow calculation through a Newton-Raphson method according to an embodiment of the present disclosure. As shown in FIG. 4, the step in which the power flow calculation is performed through the Newton-Raphson method according to the resistance between the any two of the converters or the resistance between the connection points, the DC voltage, the active power, the reactive power injection amount and the control mode of the each of the converters specifically includes the steps described below.

[0120] In S1511, an unbalance amount of the active power and an unbalance amount of the reactive power injection amount in the power flow calculation are acquired according to the resistance between the any two of the converters or the resistance between the connection points, the DC voltage, the active power, the reactive power injection amount and the control mode of the each of the converters.

[0121] In S1512, a Jacobian matrix for the power flow calculation is established according to the control mode of the each of the converters, the unbalance amount of the active power and the unbalance amount of the reactive power injection amount. When the control mode of the each of the converters is the second type control mode, and the second control mode is the constant overlap angle control mode, a Jacobian matrix parameter of the node corresponding to the each of the converters is determined merely by an AC network parameter. When the control mode of the each of the converters is the second type control mode, and the second type control mode is the constant transformation ratio control mode, the Jacobian matrix parameter of the node corresponding to the each of the converters is corrected after being determined by the AC network parameter.

[0122] In S1513, the power flow calculation is performed through the Newton-Raphson method according to the Jacobian matrix.

[0123] Specifically, different converters in the AC-DC interconnected power system have different control modes, and the active power and the reactive power injection amount of a node corresponding to one converter is related to the control mode of the one converter. For example, in a general iterative process of the power flow calculation, an unbalance equation is as follows:

[00023]$\begin{array}{c}.Math..Math.P_{\mathrm{idk}}=P_{\mathrm{is}}-V_i.Math.\underset{j=i}{\overset{n_c}{.Math.}}.Math..Math.V_j\left(G_{\mathrm{ij}}.Math.\cos .Math..Math.{}_{\mathrm{ij}}-B_{\mathrm{ij}}.Math.\sin .Math..Math.{}_{\mathrm{ij}}\right)-\left(P_{\mathrm{idk}}\right)=0\\ .Math..Math.Q_{\mathrm{idk}}=Q_{\mathrm{is}}-V_i.Math.\underset{j=i}{\overset{n_c}{.Math.}}.Math..Math.V_j\left(G_{\mathrm{ij}}.Math.\sin .Math..Math.{}_{\mathrm{ij}}-B_{\mathrm{ij}}.Math.\cos .Math..Math.{}_{\mathrm{ij}}\right)-\left(Q_{\mathrm{idk}}\right)=0\\ i=n_c+k\\ k=1,2,\mathrm{.Math.}.Math.,n_c\end{array}\}.$

[0124] Where P.sub.idk and Q.sub.idk are both scalars and positive values, and signs of P.sub.idk and Q.sub.idk are selected according to a rule that a positive sign is selected for a rectification side and a negative sign is selected for an inversion side, P.sub.is and Q.sub.is are total injection power of a system generator and a load node; .sub.ij is a difference between phase angles of the node i and the node j, and G.sub.ij and B.sub.ij are respectively the real part and the imaginary part of the admittance matrix element.

[0125] The Jacobian matrix is constructed as follows:

[00024]$\{\begin{array}{c}\left[\begin{array}{c}.Math..Math.P\\ .Math..Math.Q\end{array}\right]=J\left[\begin{array}{c}\\ V\end{array}\right]\\ J=\left[\begin{array}{cc}H& N\\ J& L\end{array}\right]\end{array},$

where H, N and L are block matrices of the Jacobian matrix, P is the unbalance amount of the active power, Q is the unbalance amount of the reactive power injection amount, and and V are correction amounts of variables in an iterative process.

[0126] When the control mode of the converter corresponding to the node is the second type control mode, and the second type control mode is the constant transformation ratio control mode, the reactive power injection amount is calculated by using the following calculation formula:

[00025]$Q_{\mathrm{dk}}=\frac{P_{\mathrm{dk}}.Math.\sqrt{-V_{\mathrm{dk}}+k_y^2+k_T^2.Math.V_a^2}}{V_{\mathrm{dk}}}$

[0127] At this time, the unbalance amount of the reactive power injection amount should be calculated by using the following equation:

[00026]$\frac{{\mathrm{dQ}}_{\mathrm{dk}}}{V_a}=\frac{k_y^2.Math.k_T^2.Math.P_{\mathrm{dk}}.Math.V_a}{V_{\mathrm{dk}}.Math.\sqrt{-V_{\mathrm{dk}}+k_y^2.Math.k_T^2.Math.V_a^2}}$

[0128] The calculated reactive power injection amount is substituted into the Jacobian matrix, and the power flow calculation in the AC-DC interconnected power system is implemented through the Newton-Raphson method.

[0129] Optionally, based on the preceding embodiment, the method for preforming the power flow calculation through the Newton-Raphson method may be further optimized. FIG. 5 is a flowchart of another method for preforming power flow calculation through a Newton-Raphson method according to an embodiment of the present disclosure. As shown in FIG. 5, performing the power flow calculation through the Newton-Raphson method includes steps described below.

[0130] In S1521, the unbalance amount of the active power and the unbalance amount of the reactive power injection amount in the power flow calculation are calculated according to the resistance, the voltage of the DC bus, the active power and the reactive power injection amount.

[0131] In S1522, the Jacobian matrix for the power flow calculation is established according to the control mode of the each of the converters, the unbalance amount of the active power and the unbalance amount of the reactive power injection amount.

[0132] In S1523, when the control mode of the each of the converters is the second type control mode, and the second type control mode is the constant transformation ratio control mode, an element Lii of the Jacobian matrix is corrected.

[0133] In S1524, whether an calculation result of the power flow satisfies a convergence condition is determined. If yes, proceed to S1525. If not, proceed to S1521, and the unbalance amount of the active power and the unbalance amount of the reactive power injection amount in the power flow calculation are acquired again.

[0134] In S1525, the power flow calculation is implemented.

[0135] Specifically, the element Lii of the Jacobian matrix is calculated by using the following calculation formula:

[00027]$L_{\mathrm{ii}}=-V_i.Math.\underset{ji,ji}{.Math.}.Math..Math.V_j\left(G_{\mathrm{ij}}.Math.\sin .Math..Math.{}_{\mathrm{ij}}-B_{\mathrm{ij}}.Math.\cos .Math..Math.{}_{\mathrm{ij}}\right)+2.Math.V_i^2.Math.B_{\mathrm{ii}}.$

[0136] When the control mode of the each of the converters is the second type control mode, and the second type control mode is the constant transformation ratio control mode, the element Lii of the Jacobian matrix also needs to be corrected as follows:

[00028]$L_{\mathrm{ii}}=-V_i.Math.\underset{ji,ji}{.Math.}.Math.V_j\left(G_{\mathrm{ij}}.Math.\sin .Math..Math.{}_{\mathrm{ij}}-B_{\mathrm{ij}}.Math.\cos .Math..Math.{}_{\mathrm{ij}}\right)+2.Math.V_i^2.Math.B_{\mathrm{ii}}-\frac{k_y^2.Math.k_T^2.Math.P_{\mathrm{dk}}.Math.V_a}{V_{\mathrm{dk}}.Math.\sqrt{-V_{\mathrm{dk}}^2+k_y^2.Math.k_T^2.Math.V_i^2}},$

where i is a node of an AC network connected to the each of the converters; V.sub.i is a voltage amplitude of the node i, G.sub.ij and B.sub.ij are respectively a real part and an imaginary part of an admittance matrix, V.sub.a is a voltage amplitude of a node connected to the each of the converters, V.sub.dk is the DC power transmission voltage; P.sub.dk is the active power flowing into the converter k, k.sub.T is the transformation ratio, k.sub.y is the converter constant, .sub.ij is a control angle of the node i, H, N and L are the block matrices of the Jacobian matrix, P is the unbalance amount of the active power, Q is the unbalance amount of the reactive power injection amount, and and V are the correction amounts of the variables in the iterative process.

[0137] After the correction, the convergence of the calculation result of the power flow needs to be verified. When the calculation result of the power flow satisfies the convergence condition, the iteration process of the power flow calculation is ended and a corresponding result is outputted. When the calculation result of the power flow does not satisfy the convergence condition, the unbalance amount of the active power and the unbalance amount of the reactive power injection amount in the power flow calculation need to be calculated again until the calculation result of the power flow satisfies the convergence condition.

[0138] In the embodiments of the present disclosure, the corresponding parameters are obtained by analyzing the control mode of each converter in the AC-DC interconnected power system, and the power flow calculation is performed through the Newton-Raphson method, thereby avoiding the problem that the calculation is not facilitated due to a large selection scale of the initial value and a large scale of the Jacobian matrix. Therefore, the AC-DC interconnected power flow calculation has better convergence, reduced calculation complexity, an improved calculation rate, and reduced costs.

[0139] An embodiment of the present disclosure further provides a power flow calculation device for an AC-DC interconnected power system, which is applicable to power flow calculation in an AC-DC interconnected power system. The AC-DC interconnected power flow calculation device in this embodiment may be implemented by software and/or hardware. FIG. 6 is a block diagram of a power flow calculation device for an AC-DC interconnected power system according to an embodiment of the present disclosure. As shown in FIG. 6, the power flow calculation device for the AC-DC interconnected power system includes a resistance acquisition module 61, a control mode acquisition module 62, a DC voltage and active power acquisition module 63, a reactive power injection amount calculation module 64 and a power flow calculation module 65.

[0140] The resistance acquisition module 61 is configured to solve a conductance matrix for a DC network of the AC-DC interconnected power system, and acquire a resistance between any two of converters in the DC network, or acquire a resistance between any two of connection points of each hierarchical structures of the DC network.

[0141] The control mode acquisition module 62 is configured to acquire a control mode of each of the converters.

[0142] The DC voltage and active power acquisition module 63 is configured to acquire a DC voltage and active power of a node corresponding to the each of the converters according to the control mode of the each of the converters.

[0143] The reactive power injection amount calculation module 64 is configured to calculate a reactive power injection amount of the each of the converters into an AC power grid according to the DC voltage, the active power, and the control mode of the each of the converters.

[0144] The power flow calculation module 65 is configured to perform power flow calculation through a Newton-Raphson method according to the resistance between the any two of the converters or the resistance between the connection points, a voltage of a DC bus, the active power, the reactive power injection amount and the control mode of the each of the converters.

[0145] In the embodiments of the present disclosure, the corresponding parameters are obtained by analyzing the control mode of each converter in the AC-DC interconnected power system, and the power flow calculation is performed through the Newton-Raphson method, thereby avoiding the problem that the calculation is not facilitated due to a large selection scale of the initial value and a large scale of the Jacobian matrix. Therefore, the power flow calculation for the AC-DC interconnected power system has better convergence, reduced calculation complexity, an improved calculation rate, and reduced costs.

[0146] An embodiment of the present disclosure further provides a storage medium including computer-executable instructions and configured to store computer programs, which when executed by a processor, are used for implementing the power flow calculation method for the AC-DC interconnected power system according to the embodiments of the present disclosure. The method includes steps described below.

[0147] A conductance matrix is solved for a DC network of an AC-DC interconnected power system, and a resistance between any two of converters in the DC network is acquired, or a resistance between any two of connection points of each hierarchical structures of the DC network is acquired.

[0148] A DC voltage and active power of a node corresponding to each of the converters are acquired according to a structure of the DC network.

[0149] A control mode of the each of the converters is acquired.

[0150] A reactive power injection amount of the each of the converters into an AC power grid is calculated according to the DC voltage, the active power, and the control mode of the each of the converters.

[0151] Power flow calculation is performed through a Newton-Raphson method according to the resistance between the any two of the converters or the resistance between the connection points, the DC voltage, the active power, the reactive power injection amount and the control mode of the each of the converters.

[0152] The storage medium is any one of various types of memory apparatus or storage apparatus. The term storage medium is intended to include a mounting medium such as a compact disc read-only memory (CD-ROM), a floppy disk or a magnetic tape device; a computer system memory or a random access memory (RAM) such as a dynamic random access memory (DRAM), a double data rate (DDR) RAM, a static random access memory (SRAM), an extended data output (EDO) RAM, or a Rambus RAM; a non-volatile memory such as a flash memory or a magnetic medium (such as a hard disk or an optical storage device); a register or other similar types of memory elements, etc. The storage medium may also include other types of memory or combinations thereof. In addition, the storage medium may be located in a first computer system in which programs are executed, or may be located in a different second computer system connected to the first computer system through a network such as the Internet. The second computer system may provide program instructions to a first computer for execution. The term storage medium may include two or more storage media which can reside at different positions, such as in different computer systems connected through a network. The storage medium may store program instructions (e.g., embodied as computer programs) which are executable by one or more processors.

[0153] Of course, in the storage medium including the computer-executable instructions according to the embodiment of the present disclosure, the computer-executable instructions implement not only the operations of the power flow calculation method for the AC-DC interconnected power system described above but also operations related to the power flow calculation method for the AC-DC interconnected power system according to any embodiment of the present disclosure.

[0154] An embodiment of the present disclosure further provides a terminal which may integrate the power flow calculation device for the AC-DC interconnected power system according to the embodiment of the present disclosure. FIG. 7 is a structural diagram of a terminal according to an embodiment of the present disclosure. As shown in FIG. 7, the terminal may include a display (not shown), a memory 101, a central processing unit (CPU) 102 (also referred to as a processor), a circuit board (not shown) and a power circuit (not shown). The CPU 102 and the memory 101 are disposed on the circuit board. The power circuit is configured to supply each circuit or component of the terminal with power. The memory 101 is configured to store computer programs. The CPU 102 reads and executes the computer programs stored in the memory 101. When executing the computer programs, the CPU 102 implements steps described below. A conductance matrix is solved for a DC network of an AC-DC interconnected power system, and a resistance between any two of converters in the DC network is acquired, or a resistance between any two of connection points of the DC network of hierarchical structures is acquired. A DC voltage and an active power of a node corresponding to each of the converters are acquired according to a structure of the DC network. A control mode of the each of the converters is acquired. A reactive power injection amount of the each of the converters into an AC power grid is calculated according to the DC voltage, the active power, and the control mode of the each of the converters. Power flow calculation is performed through a Newton-Raphson method according to the resistance between the any two of the converters or the resistance between the connection points, the DC voltage, the active power, the reactive power injection amount and the control mode of the each of the converters.

[0155] It should be understood that the illustrated terminal 100 is merely one example of the terminal, and that the terminal 100 may include more or fewer components than the components shown in the figure, may combine two or more components, or may have a different configuration of components. The various components shown in the figure may be implemented in hardware, software, or a combination of hardware and software, which includes one or more signal processing and/or application-specific integrated circuits. The terminal 100 may be, for example, a computer.

[0156] The terminal according to the embodiment of the present disclosure implements and performs the operations of the power flow calculation method for the AC-DC interconnected power system according to any embodiment of the present disclosure, and may effectively perform the power flow calculation for the AC-DC interconnected power system.

[0157] The power flow calculation device for the AC-DC interconnected power system, the storage medium and the terminal according to the above-mentioned embodiments can execute the power flow calculation method for the AC-DC interconnected power system according to any embodiment of the present disclosure and have function modules and beneficial effects corresponding to this method. For technical details not described in detail in the above-mentioned embodiments, reference can be made to the power flow calculation method for the AC-DC interconnected power system according to any embodiment of the present disclosure.

[0158] An embodiment of the present disclosure further provides a power flow calculation method for an AC-DC interconnected power system, which is applied in DC hierarchical structures. The method includes steps described below.

[0159] In step 1, a conductance matrix is solved for a DC network, and a resistance between each two of converters is acquired, or a resistance between any two of connection points of a DC power grid of hierarchical structures is acquired.

[0160] In step 2, a control mode of each of the converters is analyzed, and a voltage and active power of each node are calculated by using the following equation:

[00029]$\begin{array}{cc}\begin{array}{c}{I}_{\mathrm{dk}}={.Math.}_{j=1}^{n_c}.Math.G_{\mathrm{kj}}.Math.V_{\mathrm{dj}}\\ P_{\mathrm{dk}}=I_{\mathrm{dk}}.Math.V_{\mathrm{dk}}\end{array}\},& \left(6\right)\end{array}$

where I.sub.dk is a DC current flowing into a converter station k, and G.sub.kj is an admittance matrix element between a node k corresponding to the converter station k and a node j.

[0161] In step 3, a reactive power injection amount into an AC power grid is calculated according to the control mode, and the voltage and the active power of the each node obtained in step 2.

[0162] Step 3 includes steps described below.

[0163] In step 3.1, in condition that a control mode of a converter corresponding to a node is a constant overlap angle control mode, the reactive power injection amount is calculated by using an equation (2).

[0164] Then proceed to step 3.3; otherwise proceed to step 3.2.

[0165] In step 3.2, in condition that the control mode of the converter corresponding to the node is a constant transformation ratio control mode, the reactive power injection amount is calculated by using an equation (3), and a derivative of the reactive power injection amount with respect to an AC voltage corresponding to the reactive power injection amount is calculated by using the following equation:

[00030]$\begin{array}{cc}\frac{{\mathrm{dQ}}_{d.Math..Math.c}}{V_a}=\frac{k_{}^2.Math.K_t^2.Math.P_{d.Math..Math.c}.Math.V_a}{V_{d.Math..Math.c}.Math.\sqrt{-V_{d.Math..Math.c}^2+k_{}^2.Math.K_t^2.Math.V_a^2}}.& \left(7\right)\end{array}$

[0166] In step 3.3, with the hierarchical structures involved, a power effect of each layer on a connection point with the AC power grid is calculated according to an equation set (4).

[0167] In step 4, power flow calculation is performed through a Newton-Raphson method.

[0168] Step 4 includes steps described below.

[0169] In step 4.1, an initial value of the AC network is set and an unbalance amount of a power flow power equation is calculated.

[0170] In step 4.2, a Jacobian matrix is constructed as the following equation set:

[00031]$\begin{array}{cc}\{\begin{array}{c}\left[\begin{array}{c}.Math..Math.P\\ .Math..Math.Q\end{array}\right]=J\left[\begin{array}{c}.Math..Math.\\ .Math..Math.V\end{array}\right]\\ J=\left[\begin{array}{cc}H& N\\ J& L\end{array}\right]\end{array},& \left(8\right)\end{array}$

where J is the Jacobian matrix, P and Q are both the unbalance amount of the power equation, and and V are correction amounts of variables in an iterative process. When the control mode of the converter is the constant overlap angle, L.sub.ii corresponding to a node corresponding to the converter is determined merely by an AC network parameter. When the control mode of the converter is the constant transformation ratio, L.sub.ii corresponding to the node corresponding to the converter needs to be further corrected according to an equation (5) after being determined by the AC network parameter.

[0171] In step 4.3, the correction amounts are calculated and the AC network parameter is corrected, a convergence condition is checked, and an iteration is ended when the convergence condition is met; otherwise, proceed to step 4.1.

[0172] In step 5, a result is outputted.

[0173] Specifically, DC network modeling is described below.

[0174] (1) Traditional DC Network Model

[0175] A DC line parameter includes capacitance, inductance and the like. However, a stable operation is considered for the power flow calculation, and thus the DC line as a whole exhibits a resistance characteristic. The DC network is represented by an admittance matrix G.sub.d of nodes as:

[00032]$G_d=\left[\begin{array}{cccc}{G}_{11}& G_{12}& \mathrm{.Math.}& G_{1.Math..Math.n}\\ G_{12}& G_{22}& \mathrm{.Math.}& G_{2.Math.n}\\ \mathrm{.Math.}& \mathrm{.Math.}& & \mathrm{.Math.}\\ G_{n.Math..Math.1}& G_{n.Math..Math.2}& \mathrm{.Math.}& G_{\mathrm{nn}}\end{array}\right].$

[0176] An injection current of a node may be expressed as: I.sub.d=G.sub.dV.sub.d Where I.sub.d is an injection current of a DC node, and V.sub.d is a DC voltage.

[0177] Basic equations of the converter are listed as follows:

[00033]$\begin{array}{cc}\begin{array}{c}{V}_{d^{*}}=k_{T^{*}}.Math.V_{s^{*}}.Math.\cos .Math..Math.{}_d-X_{c^{*}}.Math.I_{d^{*}}\\ V_{d^{*}}=k_{}.Math.k_{T^{*}}.Math.V_{s^{*}}.Math.\cos .Math..Math.\\ I_{s^{*}}=k_{}.Math.k_{T^{*}}.Math.I_{d^{*}}\end{array}\}& \left(9\right)\end{array}$

where a mark * denotes a per unit value, V.sub.dci is a DC power transmission voltage, I.sub.dci is a DC power transmission current, V.sub.i.sub.si is a vector of a line voltage of an AC bus, I.sub.ci is a base frequency AC current injected to the converter, n.sub.ii is a number of bridges included in the converter, k.sub.Ti is a transformation ratio, .sub.i is a control angle of the converter which includes a gating delay angle of a rectifier and an extinction advance angle of an inverter, X.sub.ci is an overlap resistance, k.sub.=0.995 by a simplified analysis with consideration of an effect of an overlap angle, is a power factor angle corresponding to active power (absorbed by the rectifier and emitted by the inverter) and reactive power absorbed by the converter from an AC system.

[0178] (2) DC Hierarchical Structures Involved

[0179] A simple DC power transmission structure in a layered access manner is shown in FIG. 2.

[0180] Series coupling exists between DC power transmission nodes in the layered access manner. A DC node 1 and a DC node 2 shown in FIG. 2 satisfy the following relationship:

[00034]$\begin{array}{c}{I}_{\mathrm{di}.Math..Math.1}=I_{\mathrm{di}.Math..Math.2}=I_d\\ \frac{V_{\mathrm{dr}}-\left(V_{\mathrm{di}.Math..Math.1}+V_{\mathrm{di}.Math..Math.2}\right)}{R_d}=I_d\end{array}\}.$

[0181] Meanings of variables are consistent with those in an equation set (1).

[0182] The active power outputted by a single converter on a series side is proportional to a voltage ratio of the single converter, that is,

[00035]$\{\begin{array}{c}{P}_{\mathrm{idk}}=k_{\mathrm{idk}}.Math..Math..Math.P_d\\ V_{\mathrm{idk}}=k_{\mathrm{idk}}.Math..Math..Math.V_d\end{array}$

where k.sub.idk is a voltage ratio of a converter k in the hierarchical structures, P.sub.d is active power of the DC network injected into a converter station, and V.sub.d is a DC voltage of a node connected to the converter station, P.sub.idk is active power outputted by the converter k in the hierarchical structures, and V.sub.idk is a DC voltage applied across the converter k in the hierarchical structures.

[0183] (3) Control Strategy of the Converter Station

[0184] For a traditional commutated converter, each converter has two independent control variables. Assuming that a transformer tap may be adjusted seamlessly, a turn ratio k.sub.T may be linearly controlled. Therefore, active power P.sub.dc, a DC voltage V.sub.dc and a DC current I.sub.dc of a DC bus are defined as D-axis control variables; and the transformation ratio k.sub.T and the control angle of the converter are referred to as E-axis control variables.

TABLE-US-00001 TABLE 1 Control strategy of the converter D-axis Control E-axis Control Constant P.sub.dc Constant k.sub.T Constant V.sub.dc Constant Constant I.sub.dc

[0185] As for the D-axis control, in the DC network, the D-axis control of the converter at one end must be a voltage control mode, and no matter whether the D-axis control of the converters at other ends is constant P.sub.dc or constant I.sub.dc, G.sub.kj is obtained with a resistance of the DC network known, and then the voltage and the active power of the converter at each end are calculated according to an equation set (6).

[0186] The E-axis control includes two types of control.

[0187] {circle around (1)} The constant overlap angle is selected for the converter.

[0188] Output power may be expressed as:

[00036]$\begin{array}{cc}\begin{array}{c}{P}_{\mathrm{idc}}=V_{\mathrm{dk}}.Math.I_{\mathrm{dk}}\\ Q_{\mathrm{idc}}=V_{\mathrm{dk}}.Math.I_{\mathrm{dk}}.Math.\mathrm{tg}.Math..Math.{}_k\end{array}\},& \left(10\right)\end{array}$

where .sub.k is a power factor of the converter, and V.sub.dk and I.sub.dk are respectively a voltage and a current of a DC node connected to the converter.

[0189] An equation set (10) is combined with an equation set (9) to obtain the following equation:

[00037]$Q_{\mathrm{idc}}=\frac{k_{}.Math.P_{\mathrm{idc}}\left(V_{\mathrm{dk}}^2+P_{\mathrm{idc}}.Math.X_c\right).Math.\sqrt{1-\frac{V_{\mathrm{dk}}^4.Math.{\cos }^2\left({}_d\right)}{{k_{}^2\left(V_{\mathrm{dk}}^2+P_{\mathrm{idc}}.Math.X_c\right)}^2}}.Math.\sec \left({}_d\right)}{V_{\mathrm{dk}}^2}.$

[0190] {circle around (2)} The constant transformation ratio is selected for the converter.

[0191] The equation set (10) is combined with the equation set (9) to obtain the following equation:

[00038]$Q_{d.Math..Math.c}=\frac{P_{d.Math..Math.c}.Math.\sqrt{-V_{d.Math..Math.c}^2+k_{}^2.Math.K_t^2.Math.V_a^2}}{V_{d.Math..Math.c}}.$

[0192] 3. Power Flow Calculation

[0193] In the iterative process of the power flow calculation, an unbalance equation is as follows:

[00039]$\begin{array}{cc}\begin{array}{c}.Math..Math.P_{\mathrm{idc}}=P_{\mathrm{is}}-V_i.Math.\underset{j=1}{\overset{n}{.Math.}}.Math.V_j\left(G_{\mathrm{ij}}.Math.\cos .Math..Math.{}_{\mathrm{ij}}+B_{\mathrm{ij}}.Math.\sin .Math..Math.{}_{\mathrm{ij}}\right)-\left(P_{\mathrm{idc}}\right)=0\\ .Math..Math.Q_{\mathrm{idc}}=Q_{\mathrm{is}}-V_i.Math.\underset{j=1}{\overset{n}{.Math.}}.Math.V_j\left(G_{\mathrm{ij}}.Math.\sin .Math..Math.{}_{\mathrm{ij}}-B_{\mathrm{ij}}.Math.\cos .Math..Math.{}_{\mathrm{ij}}\right)-\left(Q_{\mathrm{idc}}\right)=0\\ i=n_a+k\\ k=1,2,\mathrm{.Math.}.Math.,n_c\end{array}\}& \left(11\right)\end{array}$

where P.sub.idc and Q.sub.idc are both scalars and positive values, and signs of P.sub.idc and Q.sub.idc are selected according to a rule that a positive sign is selected for a rectification side and a negative sign is selected for an inversion side; P.sub.is and Q.sub.is are total injection power of a system generator and a load node; .sub.ij is a difference between phase angles of the node i and the node j; and G.sub.ij and B.sub.ij are respectively a real part and an imaginary part of an admittance matrix element.

[0194] (1) When the constant overlap angle is selected for the E-axis control of the converter, no change is needed for the original Jacobian matrix used for the power flow calculation.

[0195] (2) When the constant transformation ratio is selected for the E-axis control of the converter,

[00040]$\frac{{\mathrm{dQ}}_{d.Math..Math.c}}{V_a}=\frac{k_{}^2.Math.K_t^2.Math.P_{d.Math..Math.c}.Math.V_a}{V_{d.Math..Math.c}.Math.\sqrt{-V_{d.Math..Math.c}^2+k_{}^2.Math.K_t^2.Math.V_a^2}}.$

[0196] In the iterative process of the power flow calculation, the unbalance amounts are consistent with those in an equation set (11).

[0197] The Jacobian matrix is modified as follows:

[00041]$L_{\mathrm{ii}}=-V_i.Math.\underset{ji,ji}{.Math.}.Math.V_j\left(G_{\mathrm{ij}}.Math.\sin .Math..Math.{}_{\mathrm{ij}}-B_{\mathrm{ij}}.Math.\cos .Math..Math.{}_{\mathrm{ij}}\right)+2.Math.V_i^2.Math.B_{\mathrm{ii}}-\frac{k_{}^2.Math.K_t^2.Math.P_{d.Math..Math.c}.Math.V_a}{V_{d.Math..Math.c}.Math.\sqrt{-V_{d.Math..Math.c}^2+k_{}^2.Math.K_t^2.Math.V_i^2}},$

where i is an AC node connected to the converter,

[00042]$-V_i.Math.\underset{ji,ji}{.Math.}.Math.V_j\left(G_{\mathrm{ij}}.Math.\sin .Math..Math.{}_{\mathrm{ij}}-B_{\mathrm{ij}}.Math.\cos .Math..Math.{}_{\mathrm{ij}}\right)+2.Math.V_i^2.Math.B_{\mathrm{ii}}$

is a calculation formula of a Jacobian matrix element L in traditional pure AC power flow calculation; V.sub.i is a voltage amplitude of the node; G.sub.ij and B.sub.ij are respectively a real part and an imaginary part of an admittance matrix; and V.sub.a is a voltage amplitude of a node connected to the converter and is numerically consistent with V.sub.i.

[0198] It is to be noted that the above are merely preferred embodiments of the present disclosure and the technical principles used therein. It will be understood by those skilled in the art that the present disclosure is not limited to the specific embodiments described herein. Those skilled in the art can make various apparent modifications, adaptations, combinations and substitutions without departing from the scope of the present disclosure. Therefore, while the present disclosure has been described in detail through the above-mentioned embodiments, the present disclosure is not limited to the above-mentioned embodiments and may include more other equivalent embodiments without departing from the concept of the present disclosure. The scope of the present disclosure is determined by the scope of the appended claims.