Calculation process and acceleration method considering far-end cooperative adjustment of multiple power electronic devices

Abstract

The invention discloses a computing process and acceleration method considering remote collaborative adjustment of various power electronic devices, which comprises the following steps: 1) acquiring input power system data, and establishing physical models of various power electronic devices; 2) calculating a modularized local sensitivity matrix based on a physical model of the power electronic equipment; 3) calculating a target physical quantity y; 4) on the basis of the modular local sensitivity matrix and the target physical quantity y, obtaining the minimum variable quantity delta u (k) of the control parameter of the power electronic equipment during the kth alternating iteration through pseudo-inverse calculation, and updating the control parameter u (k + 1) = u (k) + delta u (k); and 4) judging whether the active power of the tie line and the alternating iteration power flow meet the convergence criterion during the k-th iteration, if so, outputting the control parameters of the power electronic equipment, and otherwise, returning to the step 2). According to the method, the calculation result of the alternating iteration method can be reused, local sensitivity calculation can be directly obtained, and the solving efficiency of a far-end cooperative adjustment strategy is improved.

Description

Technical Field

[0001] The present invention relates to the field of power electronic device adjustment, and in particular, to a calculation process and a speed-up method for remotely collaborative adjustment of multiple power electronic devices. Background Art

[0002] With the continuous expansion of the scale of the power system, the technical challenges brought about by the increase in the grid-connected capacity of renewable energy and the growth in the number of deployed power electronic devices (PEDs) have become a research hotspot. The large-scale access of renewable energy has led to significant stochastic characteristics of the power system, which has further increased the probability of faults in key areas, and thus triggered cascading power outages. The rapid progress of materials science and control technology has promoted the wide application of power electronic devices in the power system. The application of power electronic devices in the power system mainly focuses on two major systems: high voltage direct current (HVDC) and flexible AC transmission systems (FACTS). As the core device of HVDC, the voltage source converter (VSC) can independently adjust active and reactive power and adjust the current direction while keeping the voltage polarity unchanged. In recent years, the research on AC-DC hybrid systems based on VSC-HVDC has been continuously increasing. The flexible AC transmission system covers a variety of devices such as unified power flow controller (UPFC), generalized unified power flow controller (GUPFC), static synchronous compensator (STATCOM), and static synchronous series compensator (SSSC), and has comprehensive functions such as voltage regulation, series compensation, power control, and phase adjustment [7]. These devices have opened up a new path for ensuring the safe and stable operation of modern power systems by realizing the refined control and dynamic optimization of electrical parameters [8].

[0003] Under normal circumstances, the adjustment of power electronic devices mainly focuses on adjusting the physical quantities at their locations, i.e., local adjustment. However, for the large-scale installation of power electronic devices, relying solely on local adjustment may lead to the inability of the physical quantities in the areas without installed devices to meet the requirements, and even the risk of over-limit may occur. Therefore, it is particularly important to study the coordinated adjustment strategy of power electronic devices. It is necessary to consider not only the comprehensive influence of each power electronic device on the controlled physical quantity, but also their overall adjustment, so as to improve the ability of power electronic devices to respond to system requirements under different operating conditions (including load changes, new energy grid connection, system faults, etc.), and enhance the flexibility and security of the operation of the entire power system.

[0004] Under normal circumstances, the adjustment of power electronic devices mainly focuses on adjusting the physical quantities at their locations, i.e., local adjustment. However, for the large-scale installation of power electronic devices, relying solely on local adjustment may lead to the inability of the physical quantities in the areas without installed devices to meet the requirements, and even the risk of over-limit may occur. Therefore, it is particularly important to study the coordinated adjustment strategy of power electronic devices. It is necessary to consider not only the comprehensive influence of each power electronic device on the controlled physical quantity, but also their overall adjustment, so as to improve the ability of power electronic devices to respond to system requirements under different operating conditions (including load changes, new energy grid connection, system faults, etc.), and enhance the flexibility and security of the operation of the entire power system.

[0005] However, the solution process of traditional optimization methods is extremely time-consuming. Considering the processing of the optimization model to improve the solution efficiency will make the accuracy control of the solution results less transparent and lack credibility. Summary of the Invention

[0006] The object of the present invention is to provide a calculation process and speed-up method for considering the remote coordinated adjustment of multiple power electronic devices, including the following steps:

[0007] 1) Obtain the input power system data and establish the physical models of various power electronic devices;

[0008] 2) Calculate the modular local sensitivity matrix based on the physical models of the power electronic devices;

[0009] 3) Calculate the target physical quantity y;

[0010] 4) Based on the modular local sensitivity matrix and the target physical quantity y, obtain the minimum change amount Δu of the control parameters of the power electronic devices at the kth alternating iteration through pseudo-inverse calculation (k) , and update the control parameter u (k+1) = u (k) + Δu (k) ;

[0011] 4) Determine whether the active power of the tie line and the alternative iterative power flow at the k-th iteration both satisfy the convergence criterion. If so, output the control parameters of the power electronic device; otherwise, return to step 2).

[0012] Furthermore, the basic data of the power system includes the topological structure of the power system, the magnitudes and locations of the sources and loads, and the types and quantities of the power electronic devices.

[0013] Furthermore, the physical models of various power electronic devices include the steady-state model of the power electronic device, the AC system model, and the DC system model;

[0014] The steady-state model of the power electronic device includes the unbalanced equations of the VSC converter station, GUPFC, STATCOM, and SSSC.

[0015] Furthermore, the unbalanced equation of the VSC converter station is as follows:

[0016]

[0017] In the formula, P c.dc is the active power flowing from the converter station to the DC system; ΔP c.dc is the unbalance of the active power flowing from the converter station to the DC system; P loss is the active power loss; P c is the active power of node c; ΔP f , ΔQ f are the unbalances of the active power and reactive power output from node f; P cf , Q cf are the active power and reactive power flowing from node f to node c; P fp , Q fp are the active power and reactive power flowing from node f to node p; Q f is the reactive power output from node f; Q p is the reactive power flowing from the converter station to node p; is the reactive power reference value; ΔQ p is the unbalance of the reactive power flowing to node p;

[0018] The unbalanced equation of GUPFC is as follows:

[0019]

[0020] In the formula, V p is the voltage amplitude of node p; is the reference voltage amplitude of node p; ΔV p is the unbalance of the voltage amplitude of node p; Q sh is the reactive power flowing from node p to the shunt side of GUPFC; is the reactive power reference value flowing from node p to the shunt side of the GUPFC; ΔQ sh is the reactive power imbalance flowing from node p to the shunt side of the GUPFC; P nm is the active power flowing from node n to node m; is the active and reactive power; ΔP nm , ΔQ nm are the active and reactive power imbalances flowing from node n to node m; PEsh = Re(Vsh·Ish * ) is the real part of the active power flowing from node p to the shunt side of the GUPFC; PEse k = Re(Vse k ·I *nm,k ) is the real part of the active power; the power loss on the DC side P dc = 0; ΔPEx is the real part of the active power imbalance;

[0021] The imbalance equation of the STATCOM is as follows:

[0022]

[0023] In the formula, Q sta , are the reactive power and reference reactive power of the STATCOM; ΔQ sta is the reactive power imbalance of the STATCOM; PEsta = Re(Vsta·Ista * ) is the real part of the active power of the STATCOM; ΔPE is the active power imbalance of the STATCOM;

[0024] The imbalance equation of the SSSC is as follows:

[0025]

[0026] In the formula, ΔPEse is the active power imbalance of the SSSC; PEse = Re(Vsssc·Isssc * ) is the real part of the active power of the SSSC.

[0027] Furthermore, the steady-state model of the power electronic device is simplified as follows:

[0028]

[0029] Among them, the vector F VSC = [ΔP c.dc , ΔP f , ΔQ f , ΔQ p T , the vector F GUPFC = [ΔV p ​ / ΔQsh,ΔP nm ,ΔQ nm ,ΔPEx] T , vector F STATCOM =[ΔV p / ΔQsh,ΔPEsta] T , vector F SSSC =[ΔP nm / ΔQ nm ,ΔPEsta] T ; x VSC , x GUPFC , x STATCOM and x SSSC respectively represent the bus voltages of the power electronic device, i.e., and f VSC , f GUPFC , f STATCOM , f SSSC is a multivariate function of x VSC , x GUPFC , x STATCOM and x SSSC .

[0030] Furthermore, the AC system model is as follows:

[0031]

[0032] where is the voltage of buses i and j; Y ij is the element in the i-th row and j-th column of the AC bus admittance matrix; P i , Q i are the active and reactive power injections on bus i respectively; ΔP i and ΔQ i are the active and reactive power imbalances on bus i;

[0033] The simplified form of the AC system model is as follows:

[0034] F AC =f AC (x AC )=0 (7)

[0035] where F AC =[ΔP,ΔQ] T , x AC is the bus voltage of the AC system; f AC is a function of voltage V AC .

[0036] The DC system model is as follows:

[0037] ΔP dci= P Gi -P Li +P c.dc,i -V dc,i Σ j V dcj Y dc.ij = 0 (8)

[0038] Wherein, V dc,i is the voltage of the DC bus i; Y dc.ij is the element at the i-th row and j-th column in the conductance matrix; P Gdc,i , P Ldc,i are respectively the active power generation and active load power on the DC bus i; ΔP dci is the active power imbalance of the DC system node; V dcj is the voltage of node j; P c.dc,i is the active power of the DC bus i;

[0039] The simplified form of the DC system model is as follows:

[0040] F dc = f dc (x dc ) = 0 (9)

[0041] Wherein, F dc = [ΔP dc1 , ΔP dc2 , ΔP dc3 , …] T , x dc represents the DC system bus voltage; f DC is a function of V DC .

[0042] Furthermore, in step 2), the steps of calculating the modular local sensitivity matrix include:

[0043] 2.1) Import the current control parameters of the power electronic devices, and based on the Newton-Raphson method, solve the steady-state models of the power electronic devices in the physical models of various power electronic devices to calculate the power flows of each module of the power electronic devices;

[0044] 2.2) Based on the power flows of each module of the power electronic devices, use the Newton-Raphson method to solve the AC system model in the physical models of various power electronic devices to obtain the AC system power flow and update the bus voltages of the AC system;

[0045] 2.3) Based on the power flows of each module of the power electronic devices, use the Newton-Raphson method to solve the DC system model in the physical models of various power electronic devices to obtain the DC system power flow and update the bus voltages of the DC system;

[0046] 2.4) Calculate the modular local sensitivity matrix of the target output y with respect to the control parameters u of the power electronic device based on the bus voltages of the AC and DC systems. That is:

[0047]

[0048] where v represents the bus voltage of the power electronic device, and f PEDs is the unbalanced equation of the power electronic device.

[0049] Furthermore, the target value y is as follows:

[0050]

[0051] where represents the current active power value of tie line i k -j k ; n t is the number of tie lines; is the set value of the active power of tie line i k -j k .

[0052] Furthermore, the minimum change amount Δu of the control parameters of the power electronic device during the k-th alternating iteration is as follows: (k) as follows:

[0053]

[0054] Furthermore, the convergence criterion satisfied by the active power of the tie line during the k-th iteration is as follows:

[0055] Δy (k) <ε AIM-in (13)

[0056] The convergence criterion satisfied by the alternating iteration power flow during the k-th iteration is as follows:

[0057] max|χ (k) -χ (k-1) |<ε AIM (14)

[0058] where χ represents all state variables in the system, including the bus voltages of the AC and DC systems and the power electronic device; the superscript (k) represents the k-th iteration; ε AIM is the convergence coefficient of the deterministic power flow.

[0059] The technical effect of the present invention is beyond doubt. The present invention proposes a remote coordination and regulation strategy for power electronic devices. This strategy first uses the Alternate Iteration Method (AIM) for preliminary calculations, and optimally adjusts power electronic devices based on local sensitivity to ensure that the convergence criteria of deterministic power flow calculations and system adjustments are met. Even if the power electronic devices in these areas are not adjacent in spatial distribution, this strategy can still accurately regulate the target values of key areas (such as the active power of tie lines).

[0060] The present invention proposes a modular local sensitivity calculation method, which does not require matrix inversion calculations, can reuse the calculation results of the alternate iteration method, and can directly obtain local sensitivity calculations, improving the solution efficiency of the remote collaborative adjustment strategy. Brief Description of the Drawings

[0061] Figure 1 Equivalent circuit diagrams for power electronic devices: (a) VSC converter station; (b) GUPFC (or UPFC); (c) STATCOM; (d) SSSC;

[0062] Figure 2 For the overall flowchart;

[0063] Figure 3 For the test case introduction;

[0064] Figure 4 For the active power of the tie line and the absolute value of the difference from the target value (|Δy|) when only considering control through GUPFC.

[0065] Figure 5 For the change of the active power of the tie line with the adjustment iteration times when using different methods in Case b: (a) AIM-out; (b) AIM-in; (c) UIM-in. Detailed Embodiment

[0066] The present invention will be further described below in conjunction with embodiments, but it should not be understood that the above-mentioned subject scope of the present invention is limited to the following embodiments. Without departing from the above technical ideas of the present invention, various substitutions and changes should be included in the protection scope of the present invention according to the common general knowledge and conventional means in the art.

[0067] Embodiment 1:

[0068] See Figures 1 to 5 , a calculation process and speed-up method for remote collaborative adjustment considering multiple power electronic devices, including the following steps:

[0069] 1) Obtain the input power system data and establish physical models of various power electronic devices;

[0070] 2) Calculate the modular local sensitivity matrix based on the physical model of the power electronic equipment;

[0071] 3) Calculate the target physical quantity y;

[0072] 4) Based on the modular local sensitivity matrix and the target physical quantity y, calculate the minimum change Δu of the control parameters of the power electronic equipment at the k - th alternating iteration through pseudo - inverse (k) , and update the control parameter u (k+1) = u (k) +Δu (k) ;

[0073] 4) Judge whether the active power of the tie - line and the alternating iteration power flow at the k - th iteration both meet the convergence criterion. If so, output the control parameters of the power electronic equipment; otherwise, return to step 2).

[0074] The basic data of the power system include the topological structure of the power system, the magnitudes and positions of the power sources and loads, and the types and quantities of the power electronic equipment.

[0075] The physical models of various power electronic equipment include the steady - state model of the power electronic equipment, the AC system model, and the DC system model;

[0076] The steady - state model of the power electronic equipment includes the unbalanced equations of the VSC converter station, GUPFC, STATCOM, and SSSC.

[0077] The unbalanced equation of the VSC converter station is as follows:

[0078]

[0079] In the formula, P c.dc is the active power flowing from the converter station to the DC system; ΔP c.dc is the unbalance of the active power flowing from the converter station to the DC system; P loss is the active power loss; P c is the active power of node c; ΔP f , ΔQ f are the unbalances of the active power and reactive power output from node f; P cf , Q cf are the active power and reactive power flowing from node f to node c; P fp , Q fp are the active power and reactive power flowing from node f to node p; Q f is the reactive power output from node f; Q p is the reactive power flowing from the converter station to node p; is the reactive power reference value; ΔQ pis the reactive power imbalance flowing into node p;

[0080] The imbalance equation of GUPFC is as follows:

[0081]

[0082] In the formula, V p is the voltage amplitude of node p; is the reference voltage amplitude of node p; ΔV p is the voltage amplitude imbalance of node p; Q sh is the reactive power flowing from node p to the shunt side of GUPFC; is the reference value of the reactive power flowing from node p to the shunt side of GUPFC; ΔQ sh is the reactive power imbalance flowing from node p to the shunt side of GUPFC; P nm is the active power flowing from node n to node m; are the active and reactive powers; ΔP nm , ΔQ nm are the active and reactive power imbalances flowing from node n to node m; PEsh = Re(Vsh·Ish * ) is the real part of the active power flowing from node p to the shunt side of GUPFC; PEse k = Re(Vse k ·I *nm,k ) is the real part of the active power; the power loss P dc = 0; ΔPEx is the real part imbalance of the active power;

[0083] The imbalance equation of STATCOM is as follows:

[0084]

[0085] In the formula, Q sta , are the reactive power and reference reactive power of STATCOM; ΔQ sta is the reactive power imbalance of STATCOM; PEsta = Re(Vsta·Ista * ) is the real part of the active power of STATCOM; ΔPE is the active power imbalance of STATCOM;

[0086] The imbalance equation of SSSC is as follows:

[0087]

[0088] In the formula, ΔPEse is the active power imbalance of SSSC; PEse = Re(Vsssc·Isssc *) is the real part of the active power of the SSSC.

[0089] The steady-state model of the power electronic device is simplified as follows:

[0090]

[0091] Among them, the vector F VSC = [ΔP c.dc , ΔP f , ΔQ f , ΔQ p T , the vector F GUPFC = [ΔV p / ΔQsh, ΔP nm , ΔQ nm , ΔPEx] T , the vector F STATCOM = [ΔV p / ΔQsh, ΔPEsta] T , the vector F SSSC = [ΔP nm / ΔQ nm , ΔPEsta] T ; x VSC , x GUPFC , x STATCOM and x SSSC respectively represent the bus voltages of the power electronic device, that is and f VSC , f GUPFC , f STATCOM , f SSSC are multivariate functions of x VSC , x GUPFC , x STATCOM and x SSSC .

[0092] The AC system model is as follows:

[0093]

[0094] Among them, is the voltage of buses i and j; Y ij is the element in the i-th row and j-th column of the AC bus admittance matrix; P i , Q i are the active and reactive power injections on bus i respectively; ΔP i and ΔQ i are the active and reactive power imbalances on bus i;

[0095] The simplified form of the AC system model is as follows:

[0096] FAC = f AC (x AC ) = 0 (7)

[0097] where F AC = [ΔP, ΔQ] T , x AC is the bus voltage of the AC system; f AC is a function of voltage V AC .

[0098] The DC system model is as follows:

[0099] ΔP dci = P Gi - P Li + P c.dc,i - V dc,i ∑ j V dcj Y dc.ij = 0 (8)

[0100] where V dc,i is the voltage of DC bus i; Y dc.ij is the element in the i-th row and j-th column of the conductance matrix; P Gdc,i , P Ldc,i are the active power generation and active load power on DC bus i respectively; ΔP dci is the active power unbalance of the DC system node; V dcj is the voltage of node j; P c.dc,i is the active power of DC bus i.

[0101] The simplified form of the DC system model is as follows:

[0102] F dc = f dc (x dc ) = 0 (9)

[0103] where F dc = [ΔP dc1 , ΔP dc2 , ΔP dc3 ,…] T , x dc represents the DC system bus voltage; f DC is a function of V DC .

[0104] In step 2), the steps of calculating the modular local sensitivity matrix include:

[0105] 2.1) Import the current control parameters of the power electronic equipment, and based on the Newton-Raphson method, solve the steady-state model of the power electronic equipment in the physical models of various power electronic devices to calculate the power flow of each module of the power electronic equipment;

[0106] 2.2) Based on the power flow of each module of the power electronic equipment, use the Newton-Raphson method to solve the AC system model in the physical models of various power electronic devices to obtain the AC system power flow and update the bus voltages of the AC system;

[0107] 2.3) Based on the power flow of each module of the power electronic equipment, use the Newton-Raphson method to solve the DC system model in the physical models of various power electronic devices to obtain the DC system power flow and update the bus voltages of the DC system;

[0108] 2.4) Based on the bus voltages of the AC and DC systems, calculate the modular local sensitivity matrix of the target output y with respect to the control parameters u of the power electronic equipment That is:

[0109]

[0110] Among them, v represents the bus voltage of the power electronic equipment, f PEDs is the unbalanced equation of the power electronic equipment,

[0111] The target value y is as follows:

[0112]

[0113] Among them, represents the current active power value of the tie line i k -j k ; n t is the number of tie lines; is the set value of the active power of the tie line i k -j k .

[0114] The minimum change amount Δu of the control parameters of the power electronic equipment during the k-th alternating iteration (k) is as follows:

[0115]

[0116] The convergence criterion satisfied by the active power of the tie line during the k-th iteration is as follows:

[0117] Δy (k) <ε AIM-in (13)

[0118] The convergence criterion satisfied by the alternating iteration power flow during the k-th iteration is as follows:

[0119] max|χ (k) -χ (k-1) |<ε AIM (14)

[0120] wherein, χ represents all state variables in the system, including the bus voltages of AC / DC systems and power electronic devices; the superscript (k) represents the k-th iteration; ε AIM is the convergence coefficient of the deterministic power flow.

[0121] Embodiment 2:

[0122] A calculation process and acceleration method for considering the remote collaborative adjustment of multiple power electronic devices, comprising the following steps:

[0123] 1) Obtain the input power system data and establish the physical models of various power electronic devices;

[0124] 2) Calculate the modular local sensitivity matrix based on the physical models of the power electronic devices;

[0125] 3) Calculate the target physical quantity y;

[0126] 4) Based on the modular local sensitivity matrix and the target physical quantity y, obtain the minimum change amount Δu of the control parameters of the power electronic devices at the k-th alternating iteration through pseudo-inverse calculation (k) and update the control parameter u (k+1) =u (k) +Δu (k) ;

[0127] 4) Determine whether the active power of the tie line and the alternating iteration power flow at the k-th iteration both satisfy the convergence criterion. If so, output the control parameters of the power electronic devices; otherwise, return to step 2).

[0128] Embodiment 3:

[0129] A calculation process and acceleration method for considering the remote collaborative adjustment of multiple power electronic devices, the technical content is the same as that of Embodiment 2. Further, the basic data of the power system includes the topological structure of the power system, the magnitude and location of the source and load, and the type and quantity of the power electronic devices.

[0130] Embodiment 4:

[0131] A calculation process and acceleration method for considering the remote collaborative adjustment of multiple power electronic devices, the technical content is the same as any one of Embodiments 2-3. Further, the physical models of various power electronic devices include the steady-state model of the power electronic device, the AC system model, and the DC system model;

[0132] The steady-state model of the power electronic device includes the unbalanced equations of the VSC converter station, GUPFC, STATCOM, and SSSC.

[0133] Example 5:

[0134] A calculation process and acceleration method for considering the remote collaborative adjustment of multiple power electronic devices, the technical content is the same as any one of Examples 2-4. Further, the unbalanced equation of the VSC converter station is as follows:

[0135]

[0136] In the formula, P c.dc is the active power flowing from the converter station to the DC system; ΔP c.dc is the unbalance of the active power flowing from the converter station to the DC system; P loss is the active power loss; P c is the active power of node c; ΔP f , ΔQ f are the unbalances of the active power and reactive power output by node f; P cf , Q cf are the active power and reactive power flowing from node f to node c; P fp , Q fp are the active power and reactive power flowing from node f to node p; Q f is the reactive power output by node f; Q p is the reactive power flowing from the converter station to node p; is the reactive power reference value; ΔQ p is the unbalance of the reactive power flowing to node p;

[0137] The unbalanced equation of GUPFC is as follows:

[0138]

[0139] In the formula, V p is the voltage amplitude of node p; is the reference voltage amplitude of node p; ΔV p is the unbalance of the voltage amplitude of node p; Q sh is the reactive power flowing from node p to the shunt side of GUPFC; is the reactive power reference value flowing from node p to the shunt side of GUPFC; ΔQ sh is the unbalance of the reactive power flowing from node p to the shunt side of GUPFC; P nm is the active power flowing from node n to node m; is the active and reactive power; ΔP nm , ΔQ nmis the active and reactive power imbalance from node n to node m; \(P_{Esh} = Re(V_{sh} \cdot I_{sh})\) * ) is the real part of the active power flowing from node p to the parallel side of GUPFC; \(P_{Ese}\) k \( = Re(V_{se}\) k \(\cdot I\) *nm,k ) is the real part of the active power; the power loss \(P\) on the DC side dc \( = 0\); \(\Delta P_{Ex}\) is the active power real part imbalance;

[0140] The imbalance equation of STATCOM is as follows:

[0141]

[0142] In the formula, \(Q\) sta , are the reactive power and reference reactive power of STATCOM; \(\Delta Q\) sta is the reactive power imbalance of STATCOM; \(P_{Esta} = Re(V_{sta} \cdot I_{sta}\) * ) is the real part of the active power of STATCOM; \(\Delta P_{E}\) is the active power imbalance of STATCOM;

[0143] The imbalance equation of SSSC is as follows:

[0144]

[0145] In the formula, \(\Delta P_{Ese}\) is the active power imbalance of SSSC; \(P_{Ese} = Re(V_{sssc} \cdot I_{sssc}\) * ) is the real part of the active power of SSSC.

[0146] Example 6:

[0147] A calculation process and acceleration method considering the remote collaborative adjustment of multiple power electronic devices, the technical content is the same as any one of Examples 2 - 5. Further, the steady-state model of the power electronic device is simplified as follows:

[0148]

[0149] Among them, the vector \(\mathbf{F}\) VSC \( = [\Delta P\) c.dc , \(\Delta P\) f , \(\Delta Q\) f , \(\Delta Q\) p T , the vector \(\mathbf{F}\) GUPFC \( = [\Delta V\) p / \(\Delta Q_{sh}\), \(\Delta P\) nm , \(\Delta Q\) nm , \(\Delta P_{Ex}]\) T , the vector \(\mathbf{F}\) STATCOM ​= [ΔV p / ΔQsh, ΔPEsta] T , vector F SSSC = [ΔP nm / ΔQ nm , ΔPEsta] T ; x VSC , x GUPFC , x STATCOM and x SSSC respectively represent the bus voltages of the power electronic device, i.e., and f VSC , f GUPFC , f STATCOM , f SSSC is a multivariate function of x VSC , x GUPFC , x STATCOM and x SSSC .

[0150] Example 7:

[0151] A calculation process and acceleration method considering remote collaborative adjustment of multiple power electronic devices, the technical content is the same as any one of Examples 2 - 6. Further, the AC system model is as follows:

[0152]

[0153] Among them, is the voltage of buses i and j; Y ij is the element in the i-th row and j-th column of the AC bus admittance matrix; P i , Q i are the active and reactive power injections on bus i respectively; ΔP i and ΔQ i are the active and reactive power imbalances on bus i;

[0154] The simplified form of the AC system model is as follows:

[0155] F AC = f AC (x AC ) = 0 (7)

[0156] Among them, F AC = [ΔP, ΔQ] T , x AC is the bus voltage of the AC system; f AC is a function of voltage V AC .

[0157] The DC system model is as follows:

[0158] ΔP dci= P Gi -P Li +P c.dc,i -V dc,i ∑ j V dcj Y dc.ij = 0 (8)

[0159] Among them, V dc,i is the voltage of the DC bus i; Y dc.ij is the element in the i-th row and j-th column of the conductance matrix; P Gdc,i , P Ldc,i are respectively the active power generation and active load power on the DC bus i; ΔP dci is the active power imbalance of the DC system node; V dcj is the voltage of node j; P c.dc,i is the active power of the DC bus i;

[0160] The simplified form of the DC system model is as follows:

[0161] F dc = f dc (x dc ) = 0 (9)

[0162] Among them, F dc = [ΔP dc1 , ΔP dc2 , ΔP dc3 ,…] T , x dc represents the DC system bus voltage; f DC is a function of V DC .

[0163] Example 8:

[0164] A calculation process and acceleration method considering the remote collaborative adjustment of various power electronic devices. The technical content is the same as any one of Examples 2-7. Further, in step 2), the steps of calculating the modular local sensitivity matrix include:

[0165] 2.1) Import the current control parameters of the power electronic device, and based on the Newton-Raphson method, solve the steady-state model of the power electronic device in the physical model of various power electronic devices, and calculate the power flow of each module of the power electronic device;

[0166] 2.2) Based on the power flow of each module of the power electronic device, use the Newton-Raphson method to solve the AC system model in the physical model of various power electronic devices, obtain the AC system power flow, and update the bus voltage of the AC system;

[0167] 2.3) Based on the power flow of each module of the power electronic device, use the Newton-Raphson method to solve the DC system model in the physical models of various power electronic devices, obtain the DC system power flow, and update the bus voltages of the DC system;

[0168] 2.4) Based on the bus voltages of the AC and DC systems, calculate the modular local sensitivity matrix of the target output y with respect to the control parameters u of the power electronic device That is:

[0169]

[0170] where v represents the bus voltage of the power electronic device, and f PEDs is the unbalanced equation of the power electronic device,

[0171] Example 9:

[0172] A calculation process and acceleration method considering the remote collaborative adjustment of multiple power electronic devices, the technical content is the same as any one of Examples 2-8. Further, the target value y is as follows:

[0173]

[0174] where, represents the current value of the active power of the tie line i k -j k ; n t is the number of tie lines; is the set value of the active power of the tie line i k -j k .

[0175] Example 10:

[0176] A calculation process and acceleration method considering the remote collaborative adjustment of multiple power electronic devices, the technical content is the same as any one of Examples 2-9. Further, the minimum change amount Δu of the control parameters of the power electronic device during the k-th alternating iteration (k) is as follows:

[0177]

[0178] Example 11:

[0179] A calculation process and acceleration method considering the remote collaborative adjustment of multiple power electronic devices, the technical content is the same as any one of Examples 2-10. Further, the convergence criterion satisfied by the active power of the tie line during the k-th iteration is as follows:

[0180] Δy (k) <ε AIM-in (13)

[0181] The convergence criterion satisfied by the alternating iterative power flow at the k-th iteration is as follows:

[0182] max|χ (k) -χ (k-1) |<ε AIM (14)

[0183] Wherein, χ represents all state variables in the system, including the bus voltages of AC / DC systems and power electronic devices; the superscript (k) represents the k-th iteration; ε AIM is the convergence coefficient of the deterministic power flow.

[0184] Embodiment 12:

[0185] A calculation process and speed-up method for considering the coordinated adjustment of multiple power electronic devices at a remote end, the steps include:

[0186] Step S1: Input power system data and establish physical models of various power electronic devices;

[0187] Step S2: Perform a loop of alternating iterative calculations and calculate the local sensitivity matrix based on modular local sensitivity;

[0188] Step S3: Set the set value of the target physical quantity and calculate the target value y.

[0189] Step S4: Obtain the minimum norm solution through pseudo-inverse calculation, and then the minimum change amount Δu (k) of the control parameters of the power electronic device at the k-th alternating iteration can be obtained, and the parameters are updated.

[0190] Step S5: Judge whether the active power of the tie line and the alternating iterative calculation at the k-th iteration satisfy the convergence criterion. If so, output the control parameters of the power electronic device.

[0191] The characteristics of Step S1 are as follows:

[0192] S11: Steady-state model of power electronic device

[0193] The equivalent circuits of VSC converter stations, GUPFC (or UPFC), STATCOM, and SSSC are as Figure 1 shown. Note that the steady-state model of the VSC converter station takes into account the losses related to power exchange, while other models, such as GUPFC (or UPFC), SSSC, and STATCOM, ignore these power losses. Considering Figure 1 the positive direction in, the local regulation and unbalance equations of the power electronic device are as follows:

[0194] (1) VSC converter station

[0195] The detailed control parameters for the local regulation of the VSC converter station are given in Table 1.

[0196] Table 1 Control Parameters and Types of VSC Converter Station

[0197]

[0198] For the power station, since the AC bus voltage and the complex power S p = P p + jQ p are known, the voltage of the power station and can be calculated by the following formula:

[0199] For the voltage station, the unbalanced equation is:

[0200]

[0201] where P c.dc is known and Q s ref is the set value.

[0202] (2) GUPFC (or UPFC)

[0203] For the shunt side, GUPFC (or UPFC) can control the magnitude of the AC bus voltage or the injected reactive power For the series side, GUPFC (or UPFC) can control the power flow The unbalanced equation of GUPFC is:

[0204]

[0205] where PEsh = Re(Vsh · Ish * ) and PEse k = Re(Vse k · I *nm,k ), and the "*" in the superscript is the conjugate of the complex number. In addition, the power loss on the DC side is ignored, i.e., P dc ≈ 0.

[0206] (3) STATCOM

[0207] STATCOM can control the magnitude of the AC bus voltage or the injected reactive power The unbalanced equation of GUPFC is as follows:

[0208]

[0209] where P_Esta = Re(V_sta · I_sta * ), and P dc ≈ 0.

[0210] (4) SSSC

[0211] The SSSC can control the active power or the reactive power The unbalanced equation of the SSSC is:

[0212]

[0213] where P_Ese = Re(V_sssc · I_sssc * ), and P dc ≈ 0.

[0214] The unbalanced equations of the power electronic devices in Eqs. (1)-(4) can be simplified to:

[0215]

[0216] where F VSC = [ΔP c.dc , ΔP f , ΔQ f , ΔQ p T , F GUPFC = [ΔV p / ΔQ_sh, ΔP nm , ΔQ nm , ΔP_Ex] T , F STATCOM = [ΔV p / ΔQ_sh, ΔP_Esta] T , F SSSC = [ΔP nm / ΔQ nm , ΔP_Esta] T , x VSC , x GUPFC , x STATCOM and x SSSC respectively represent the bus voltages of the power electronic devices, i.e., and

[0217] S12: AC System Model

[0218] The power unbalanced equation of bus i in the AC system is:

[0219]

[0220] where,​ is the voltage of bus i / j; Y ij is the element in the i-th row and j-th column of the AC bus admittance matrix; P i , Q i are the active and reactive power injections (loads, generation, and power injections from power electronic devices) at bus i respectively. And Equation (6) can be simplified to:

[0221] F AC = f AC (x AC ) = 0 (7)

[0222] where F AC = [ΔP, ΔQ] T , and x AC is the bus voltage of the AC system.

[0223] S13: DC system model

[0224] The unbalanced equation in the DC system is:

[0225] ΔP dci = P Gi - P Li + P c.dc,i - V dc,i Σ j V dcj Y dc.ij = 0 (8)

[0226] where V dc,i is the voltage of DC bus i; Y dc.ij is the element in the i-th row and j-th column of the conductance matrix; P Gdc,i , P Ldc,i are the active power generation and active power load at DC bus i respectively. And Equation (8) can be simplified to:

[0227] F dc = f dc (x dc ) = 0 (9)

[0228] where F dc = [ΔP dc1 , ΔP dc2 , ΔP dc3 ,…] T , and x dc represents the DC system bus voltage.

[0229] 3. According to the calculation process and acceleration method for considering the remote collaborative adjustment of multiple power electronic devices, the characteristics of step S2 are as follows:

[0230] There are mainly two methods for remote control, namely Π-out and Π-in (proposed in this paper), where Π represents any other iterative solution method such as the Alternating Iterative Method (AIM), the Unified Iterative Method (UIM), or DPF. The process of Π-out is: first, execute DPF (meeting the convergence criterion of DPF), and then adjust the flexibility setting (meeting the convergence criterion of adjustment); the process of Π-in is: first, perform one iteration in DPF, and then adjust the power electronic device to meet the convergence criterion of DPF and the adjustment. This patent proposes AIM-in.

[0231] S21: Perform one iterative calculation using the Alternating Iterative Method

[0232] Import the current control parameters of the power electronic device, and then calculate the power flow equations of each module of the power electronic device based on the Newton-Raphson method. It should be noted that each module is independent and can be theoretically calculated in parallel; then calculate the power injected from the power electronic device into the AC system and solve the AC system power flow equation based on the Newton-Raphson method, and then determine the power injected from the VSC converter station into the DC system and solve the DC system power flow equation based on the Newton-Raphson method; and update the bus voltages of the AC and DC systems.

[0233] S22: Solution method for local sensitivity

[0234] To achieve remote coordinated regulation through the power electronic device, it is necessary to calculate the local sensitivity matrix of the target output y (such as the active power of the tie line) with respect to the control parameters u of the power electronic device

[0235] Traditional local sensitivity

[0236] The local sensitivity matrix obtained through the traditional solution method Focuses on the coupling relationship of the entire system. The target output y is affected by the state variables χ and u, so y can be expressed as y = y(χ, u). The power flow equation of the entire system is:

[0237] f w (χ, u) = 0 (10)

[0238] where, f w is the unbalanced equation of the entire system.

[0239] Can be obtained by the following formula:

[0240]

[0241] Modular local sensitivity

[0242] The local sensitivity matrix can also be obtained through modular local sensitivity It is calculated. The modular local sensitivity depends on intermediate variables, i.e., the active and reactive power injections s (from loads, power generation equipment, and power electronic equipment) of each bus, which are used as a bridge to indirectly consider the linear relationship between modules. The relationship between y and u can be expressed as:

[0243]

[0244] where v represents the bus voltage of the power electronic equipment, and f PEDs is the unbalanced equation of the power electronic equipment (see Equation (5)),

[0245] Actually, calculating the local sensitivity matrix (including the traditional local sensitivity matrix and the modular local sensitivity matrix ) mainly depends on the calculation of the inverse matrix. When using the Unified Iteration Method (UIM), since the inverse matrix in has been calculated, the traditional local sensitivity method is more applicable. When using AIM, since the inverse matrix of and of have been calculated separately in each module, this indicates that the modular local sensitivity method is more suitable.

[0246] 4. According to the calculation process and acceleration method for considering the coordinated adjustment of multiple power electronic devices at the far end, the characteristics of step S3 are as follows:

[0247] Compared with local regulation, in far-end regulation, the physical quantity to be controlled in the system is not adjacent to the power electronic equipment, and a regulation strategy for the power electronic equipment needs to be designed to regulate the active power of the tie line. If the set value of the active power of the tie line i k -j k is determined and the target value y is determined by the following formula.

[0248]

[0249] where represents the current value of the active power of the tie line i k -j k ; n t is the number of tie lines. Since the target value y = 0, the difference from the target value y is Δy = 0 - y = -y. Obviously, y > 0, and the smaller y is, the closer the actual active power of the tie line i k -j k is to the set value.

[0250] 5. According to the AC-DC power grid remote coordination regulation strategy considering power electronic devices, the features of step S4 are as follows:

[0251] Δy and Δu satisfy Since the dimension of u is always greater than or equal to that of y, the minimum norm solution is obtained through pseudo-inverse calculation, and then the minimum change amount Δu of the control parameters of the power electronic device at the k-th alternating iteration can be obtained. (k) .

[0252]

[0253] Furthermore, the control parameter u of the power electronic device is updated. (k) : u (k+1) = u (k) + Δu (k) .

[0254] 6. According to the calculation process and speed-up method for remote collaborative adjustment of multiple power electronic devices, the features of step S5 are as follows:

[0255] Judge whether the active power of the tie line at the k-th iteration meets the convergence criterion for adjustment:

[0256] Δy (k) < ε AIM-in (15)

[0257] And judge whether two iterations (the k-th iteration and the (k - 1)-th iteration) meet the convergence criterion of deterministic power flow through the following formula.

[0258] max|χ (k) - χ (k-1) | < ε AIM (16)

[0259] Where χ represents all state variables in the system, including the bus voltages of the AC-DC system and power electronic devices; the superscript (k) represents the k-th iteration; ε AIM is the convergence coefficient of deterministic power flow. Usually, ε AIM = 10 -6 . If the convergence criterion in formula (16) is not met, start the iteration again from the power flow calculation step; otherwise, output the final calculation result and terminate the calculation process. ε AIM-in is the convergence coefficient for adjustment. In this patent, ε AIM-in = 10 -4 . If any one of formula (16) or formula (15) is not met, repeat steps S2 to S5. Otherwise, calculate other required data.

[0260] Example 13:

[0261] Verification of a calculation process and speed-up method considering remote collaborative adjustment of multiple power electronic devices is as follows:

[0262] The example in this embodiment aims to illustrate the remote coordination regulation characteristics and local regulation characteristics of power electronic devices. Figure 3 Among them, lines 23-24, 30-38, and 33-37 are tie lines between two regions, and their rated active powers are 1.1 p.u., 2.0 p.u., and 0.2 p.u., respectively. Assume The basic parameters of power electronic devices (such as equivalent impedance, control coefficient, etc.) are from " Ayan, Optimizing reactive power flow of HVDC systems using geneticalgorithm[J].International Journal of Electrical Power&Energy Systems,2014,55:1-12", "Z.Fan, Z.Yang, K.Xie, etal.General Steady-State Modeling andLinearization of Power Electronic Devices in AC-DC Hybrid Grid[J].IEEETransactions on Power Systems,2021,36(6):5746-5755", and the detailed information is shown in Table 2. It should be noted that the control parameter P37-33 of GUPFC can directly adjust the active power of tie line 33-37, and in example a, only the P 37-33 parameters (i.e., local regulation) of GUPFC are considered. Example b further considers the remote coordination regulation of power electronic devices.

[0263] The example parameters are from MATPOWER, and all tests are carried out on a computer with an Intel(R) Core(TM) i7-13700K CPU @ 3.40 GHz and 32 GB 6400 MHz DDR5 using MATLAB 2023a.

[0264] Table 2 Parameters of the example

[0265]

[0266] (1) Local regulation of power electronic devices

[0267] To illustrate the limitations of local regulation, Case a was designed. This case only considered the control parameter P of the GUPFC 37-33 , which can directly adjust the active power of the tie line 33 - 37 (i.e., local regulation). Since the influence of the P 37-33 parameter of the GUPFC on the active power of the other two tie lines is not clear, the P 37-33 parameter of the GUPFC was adjusted near the opposite number of (the power direction of P 37-33 is opposite to that of ), that is, the P 37-33 parameter of the GUPFC was changed from 0.11 to 0.2. The active power of the tie line and the absolute value of the difference from the target value (|Δy|) are shown as Figure 4 .

[0268] According to Figure 4 , the minimum value of |Δy| was set to 0.2756 for the P 37-33 of the GUPFC (not satisfying Equation (15)), which indicates that it is difficult to make the active power of the tie line reach the target value only through the local regulation of the P 37-33 of the GUPFC. In addition, as the P 37-33 of the GUPFC changes, it can be determined that the active power of the tie line 33 - 37 is almost equal to the opposite number of the P 37-33 of the GUPFC, but no intuitive relationship between the active power of the tie lines 23 - 24 and 30 - 38 and the P37 - 33 of the GUPFC can be obtained. Therefore, it is unwise to change the active power of the tie lines where non - adjacent (power electronic devices) are located through local regulation. Unless the control parameters of the power electronic devices are continuously adjusted manually after each DPF calculation, the influence of the regulation of the power electronic devices on the expected target is not clear. Thus, it is necessary to study the remote coordinated regulation strategy to give full play to the synergistic effect of various power electronic devices, so as to effectively achieve the regulation goals of all tie lines.

[0269] (2) Remote Coordinated Regulation of Power Electronic Devices

[0270] To verify the performance of the AIM combined with modular local sensitivity in achieving the remote coordinated regulation of power electronic devices, Case b was designed. Three methods were adopted in Case b, namely AIM - out, AIM - in (the method proposed in this patent), and UIM - in. The variation of the active power with the number of regulation iterations under different methods is shown in the figure.

[0271] As can be seen from the figure, the results show that all three methods can achieve remote coordinated regulation of power electronic devices. In addition, as can be seen from the figure, AIM-in requires a total of 22 adjustment iterations, but it is the most efficient (1.11 seconds) in achieving remote coordinated regulation. Moreover, regardless of the method used, the green area in the figure indicates that under the same convergence criterion (here ε AIM-in = 10 -4 ), some fine-tuning is required to reach the target. And if ε AIM-in is appropriately increased, the number of adjustment iterations can be reduced.

Claims

1. A calculation process and speed-up method considering remote coordinated adjustment of multiple power electronic devices, characterized in that: The following steps are involved: 1) Obtain input power system data and establish physical models of various power electronic equipment; 2) Calculate the modular local sensitivity matrix based on the physical model of the power electronics equipment; 3) Calculate the target physical quantity y; 4) Based on the modular local sensitivity matrix and the target physical quantity y, the minimum change Δu of the control parameters of the power electronic equipment at the kth alternating iteration is obtained by pseudo-inverse calculation (k) , and update the control parameter u (k+1) =u (k) +Δu (k) ; 4) Determine whether the active power of the tie line and the alternating iterative power flow meet the convergence criteria at the kth iteration. If so, output the control parameters of the power electronic equipment, otherwise return to step 2).

2. According to claim 1, a calculation process and speed-up method considering remote coordinated adjustment of multiple power electronic devices is characterized in that: The basic data of the power system include the topological structure of the power system, the size and location of sources and loads, and the type and quantity of power electronic equipment.

3. According to claim 1, a calculation process and speed-up method considering remote coordinated adjustment of multiple power electronic devices is characterized in that: The physical models of various power electronic equipment include steady-state models of power electronic equipment, AC system models, and DC system models; The power electronic equipment steady-state model includes unbalanced equations of VSC converter station, GUPFC, STATCOM and SSSC.

4. A calculation process and speed-up method considering remote coordinated adjustment of multiple power electronic devices according to claim 3, characterized in that: The unbalance equation of the VSC converter station is as follows: Where P c.dc is the active power flowing from the converter station to the DC system; ΔP c.dc P is the unbalanced amount of active power flowing from the converter station to the DC system; loss is the active power loss; P c is the active power of node c; ΔP f , ΔQ f P is the unbalanced amount of active power and reactive power output by node f; cf , Q cf is the active power and reactive power flowing from node f to node c; P fp , Q fp is the active power and reactive power flowing from node f to node p; Q f is the reactive power output from node f; Q p is the reactive power flowing from the converter station to the node p; is the reactive power reference value; ΔQ p is the unbalanced reactive power flowing to node p; The unbalanced equation of GUPFC is as follows: Where V p is the voltage amplitude of node p; is the reference voltage amplitude of node p; ΔV p is the voltage amplitude imbalance at node p; Q sh is the reactive power flowing from node p to the parallel side of GUPFC; is the reactive power reference value flowing from node p to the parallel side of GUPFC; ΔQ sh is the unbalanced reactive power flowing from node p to the parallel side of GUPFC; P nm is the active power flowing from node n to node m; is the active and reactive power; ΔP nm , ΔQ nm is the unbalanced amount of active and reactive power flowing from node n to node m; PEsh = Re (Vsh·Ish * ) is the real part of the active power flowing from node p to the parallel side of GUPFC; PEse k =Re(Vse k I *nm,k ) is the real part of active power; the power loss on the DC side is P dc =0; ΔPEx is the real part unbalance of active power; The unbalance equation of STATCOM is as follows: In the formula, Q sta , is the reactive power of STATCOM and the reference reactive power; ΔQ sta is the reactive power unbalance of STATCOM; PEsta=Re(Vsta·Ista * ) is the real part of the active power of STATCOM; ΔPE is the active power unbalance of STATCOM; The unbalanced equation for SSSC is shown below: Where ΔPEse is the unbalanced active power of SSSC; PEse = Re(Vsssc·Isssc * ) is the real part of the SSSC active power.

5. A calculation process and speed-up method considering remote coordinated adjustment of multiple power electronic devices according to claim 4, characterized in that: The steady-state model of power electronic equipment is simplified as follows: Among them, the vector F VSC =[ΔP c.dc ,ΔP f ,ΔQ f ,ΔQ p ] T , vector F GUPFC =[ΔV p / ΔQsh,ΔP nm ,ΔQ nm ,ΔPEx] T , vector F STATCOM =[ΔV p / ΔQsh,ΔPEsta] T , vector F SSSC =[ΔP nm / ΔQ nm ,ΔPEsta] T ;x VSC , x GUPFC , x STATCOM and x SSSC They represent the bus voltage of the power electronic equipment, namely and f VSC 、f GUPFC 、f STATCOM 、f SSSC For x VSC , x GUPFC , x STATCOM and x SSSC A multivariate function of .

6. The calculation process and speed-up method considering remote coordinated adjustment of multiple power electronic devices according to claim 2, characterized in that: The AC system model is shown below: in, is the voltage of busbars i and j; Y ij is the element in the i-th row and j-th column of the AC bus admittance matrix; P i , Q i are the active and reactive power injection on bus i respectively; ΔP i and ΔQ i is the unbalanced amount of active and reactive power on bus i; The simplified form of the AC system model is as follows: F AC =f AC (x AC )=0 (7) Among them, F AC =[ΔP,ΔQ] T , x AC is the bus voltage of the AC system; f AC is about voltage V AC Function of The DC system model is shown below: ΔP dci =P Gi -P Li +P c.dc,i -V dc,i ∑ j V dcj Y dc.ij =0 (8) Among them, V dc,i is the voltage of DC bus i; Y dc.ij is the element in the i-th row and j-th column of the conductivity matrix; P Gdc,i , P Ldc,i are respectively the active generation power and active load power on DC bus i; ΔP dci is the unbalanced active power of the DC system node; V dcj is the voltage at node j; P c.dc,i is the active power of DC bus i; The simplified form of the DC system model is as follows: F dc =f dc (x dc )=0 (9) Among them, F dc =[ΔP dc1 ,ΔP dc2 ,ΔP dc3 ,…] T , x dc represents the DC system bus voltage; f DC It's about V DC function.

7. The calculation process and speed-up method considering remote coordinated adjustment of multiple power electronic devices according to claim 1, characterized in that: In step 2), the step of calculating the modular local sensitivity matrix includes: 2.1) Import the current control parameters of the power electronic equipment, and solve the steady-state model of the power electronic equipment in the physical model of various power electronic equipment based on the Newton-Ray method, and calculate the power flow of each module of the power electronic equipment; 2.2) Based on the power flow of each module of the power electronic equipment, the Newton-Ray method is used to solve the AC system model in the physical model of various power electronic equipment, obtain the AC system power flow, and update the bus voltage of the AC system; 2.3) Based on the power flow of each module of the power electronic equipment, the DC system model in the physical model of various power electronic equipment is solved by using the Newton-Ray method to obtain the DC system power flow and update the bus voltage of the DC system; 2.4) Based on the bus voltage of the AC and DC systems, the modular local sensitivity matrix of the target output y relative to the control parameter u of the power electronic equipment is calculated Right now: Where, v represents the bus voltage of the power electronic equipment, f PEDs is the unbalanced equation for power electronic equipment, 8. The calculation process and speed-up method considering remote coordinated adjustment of multiple power electronic devices according to claim 1, characterized in that: The target value y is as follows: in, Represents contact line i k -j k The current value of active power; n t is the number of tie lines; For contact line i k -j k The active power setting value.

9. The calculation process and speed-up method considering remote coordinated adjustment of multiple power electronic devices according to claim 1, characterized in that: The minimum change Δu of the control parameters of the power electronic equipment during the kth alternating iteration (k) As shown below:

10. The calculation process and speed-up method considering remote coordinated adjustment of multiple power electronic devices according to claim 1, characterized in that: The convergence criterion satisfied by the active power of the tie line at the kth iteration is as follows: Δy (k) <e AIM-in (13) The convergence criterion satisfied by the alternating iterative power flow at the kth iteration is as follows: max|x (k) -x (k-1) |<e AIM (14) in, χ represents all state variables in the system, including the bus voltages of the AC and DC systems and power electronic equipment; the superscript (k) indicates the kth iteration; ε AIM is the convergence coefficient of the deterministic power flow.

Images