Preventive-Security-Constrained Optimal Power Flow Model Considering IPFC Control Modes

<jats:p>The interline power flow controller (IPFC) is one of the most versatile integrated flexible alternating current transmission systems (FACTS) controllers and can realize power flow control for multiple transmission lines in modern power systems. However, control characteristics are ignored in conventional IPFC models, in which unreasonable assumptions about injected voltages may lead to security problems in realistic operation. Besides, preventive security constraints considering IPFC control modes are not included in optimal power flow (OPF) control of the system with IPFC, squandering IPFC control potential. To solve these problems, a preventive-security-constrained optimal power flow (PSCOPF) model considering IPFC control modes is proposed in this paper. IPFC control characteristics under different control modes are analyzed and employed as constraints of the optimization model. The iterative updates of converter output voltages for different control modes are derived respectively for power flow calculation, and the power and voltages required in the objective function and constraints of the proposed model can then be obtained. Through optimal selection of IPFC control modes and control parameters, the proposed model can better reconcile the economical and secure operation of the system. Numerical results demonstrate the efficient performance and superiority of the PSCOPF model considering IPFC control modes.</jats:p>

section_content: With the expansion of power grids and the integration of renewable energy, the problem of transmission overload is becoming increasingly serious [1]. The challenge deteriorates the economic operation of power systems and may even lead to severe threats to the safety of the system [2]. FACTS technologies provide rapid voltage or power flow regulation to overcome the challenge, thus effectively improving the operation of the system [3]. IPFC is one of the most potent and versatile integrated FACTS controllers [4]. With excellent potential and broad application prospects, IPFC provides full control for transmitted active power and reactive power flow simultaneously for multiple transmission lines [5]. IPFC has the capability to adjust main transmission parameters including voltage amplitude and angle, line impedance, and power flow based on corresponding control modes [6]. It is characterized as one of the best-featured FACTS devices, with an appropriate strategy, IPFC can improve static security of the power systems and mitigate power transmission bottlenecks [7, 8]. 

For the system with IPFC, optimal power flow (OPF) is generally used to improve the economy and security of the power system in steady state [9]. OPF has been attracting researchers due to its contributions to the steady operation of the system, which is extremely important for a system with large regulation demand due to the integration of renewable energy [10]. The optimization of security-constrained economic dispatch is studied in [11], which sets economic indicators as the objective function and safety indicators as constraints [12]. OPF in [13] is used to realize congestion management and improve the economic operation of the system. OPF in [14] sets multi-objectives, including generation cost, active power loss, and voltage stability, for the economic and secure operation of the system. The obtained optimal operating point struck a balance between minimization of operating cost, voltage profile deviation, and feeder congestion [15]. 

IPFC models are required to be constructed at first for the OPF of the system. Variableimpedance-based models were proposed in [16], and impedance was introduced as decision variables to realize an economic optimum. However, the physical characteristics and operation constraints are hard to describe in terms of impedance, and low-impedance branches could cause ill conditioning of the admittance matrix. Voltage source models and power injection models are more suitable for the OPF problem [17]. In voltage source models, converters of IPFC devices are equalized by controllable voltage sources, but the introduction of additional buses in voltage source models leads to an asymmetric admittance matrix [18]. In power injection models, injected power and symmetric admittance matrices can be obtained and used for Newton-Raphson power flow calculation [19], for which Jacobian codes can be reused to reduce redundant matrix correction and programming complexity [20]. 

However, there are still limitations in the current research on OPF with IPFC [21]. Continuity and reliability of the power supply are important in existing power systems, and the current research about OPF focuses on the steady state of the power system [22, 23]. IPFC control modes will affect power flow distribution after contingencies, but conventional IPFC and other FACTS models ignore control characteristics and regard injected voltages as unchanged before and after contingencies [24]. The unreasonable assumptions about injected voltages may lead to security problems in realistic operation. Additionally, remote contingencies are difficult to detect, so IPFC is hard to adjust accurately and rapidly after contingencies, which requires early selection of the IPFC control mode to realize preventive control [25]. 

To solve the OPF model with an IPFC, a nonlinear interior point algorithm incorporating a generalized IPFC is implemented in [26] to solve the OPF model and minimize operating costs. A sequential linear programming algorithm for the nonconvex OPF problem is proposed in [27], which captures transmission losses and reactive power, realizing more accurate locational marginal prices using LP solvers. Multi-objective OPF problems with IPFC were solved by an adaptive clonal selection algorithm and sequential quadratic programming respectively in [14, 28]. However, the iterative update schemes for different IPFC control modes are not devised and no practical method can be used to calculate the required power and voltage in the optimization, making it difficult to solve the OPF model when IPFC control modes are considered [29]. 

To overcome the drawbacks above, a PSCOPF model considering IPFC control modes and a corresponding solution method are proposed. Note that we concentrate on 15 min level optimization in the paper without consideration of dynamic characteristics, which will be studied in the future research. The main contributions of this paper include:

(1) A PSCOPF model considering IPFC control modes is established to fully utilize IPFC control potential, improving the economy and security of the system. IPFC control modes affect power flow distribution after contingencies and control characteristics under different control modes are analyzed, which are then employed as constraints of the optimization model. Economy and security margin are taken as the multi-objective function in the optimization, and the preventive-security-constrained technique is involved in realizing the optimal selection of IPFC control modes and control parameters in advance. (2) The corresponding method to solve the proposed model is deduced. Iterative schemes of converter output voltages and equivalent injected power for different IPFC control modes are derived respectively, and the power and voltages required in the proposed model can then be obtained. The calculated power and voltages in the objective function and constraints are further used to obtain the optimal result of the proposed model. (3) Advantages of the proposed model are proved in the case study through comparison with the original system, with the conventional model, and with optimization not considering IPFC control modes. Numerical results show that the optimal solution of the proposed model can decrease operation costs, improve voltage stability, and eliminate the risk of overload. 

The remainder of this paper is as follows. The PSCOPF model considering IPFC control modes is constructed in Section 2. The solution method for the proposed model is provided in Section 3. The case study and conclusion are given in Sections 4 and 5, respectively.
section_name: Introduction

section_name: PSCOPF Model Considering IPFC Control Modes

section_content: To realize the optimal operation of the system in economy and security, the multiobjective function of the PSCOPF model can then be expressed as:

where F is the overall objective function, E C is the economy function, S e is the security function, and λ is the weight coefficient to keep the economy and security indexes at the same level and measure the respective ratio. 

(1) Economy Function

In addition to IPFC regulation, some generators can also participate in optimal power flow control. The economy objective function intends to minimize the operating cost of the system under the same load distribution and can be calculated as:

where N g is the number of generators involved in optimal regulation; P gi is the active power output of generator i; and c 2i, c 1i, and c 0i are the cost coefficients of generator i. 

(2) Security Function

As renewable energy integration increases, the intermittent and stochastic characteristics of the output power have a significant impact on system operation, especially voltage stability [30]. The L-index, proposed by Kessel and Glavitsch in 1986, is used to predict voltage stability margins and is included in the security objective function [31]. For the PQ node β L and the PV node β G in the system, the corresponding relationship between voltages and currents can be expressed as:

where V L, I L, V G, and I G are voltage and current vectors of the PQ and PV nodes, respectively, and Z LL, F LG, K GL, and Y GG are corresponding matrices and can be calculated from the admittance matrix as:

Energies 2024, 17, 1660 4 of 16

For any PQ node j, the voltage stability index L j can be calculated as:

where F ji is the (j,i) element in matrix F LG and is the load participation factor, V i ∠θ i is the ith PV node voltage in polar coordinates while i∈β G, and V j ∠θ j is the jth PQ node voltage in polar coordinates while j∈β L. L j will reach 1. 0 when the voltage of bus j reaches its collapse point [32]. Therefore, if L j of a PQ node reaches 1. 0, the system is at the critical point of voltage stability. If L j of a PQ node is larger than 1. 0, voltage instability has occurred. If L j is less than 1. 0 for all PQ nodes, the system is voltage-stable. 

V si is the voltage stability index of the system under normal operation and can be expressed as:

To maintain stability before and after contingencies, the preventive-security-constrained technique is integrated to optimize the security of the system under various circumstances. After N -1 contingency, the N -1 line disconnects and admittance matrix Y and corresponding matrices in (4) change accordingly. Voltage stability index L i (c) and index of the system V si (c) after N -1 contingency can be obtained as:

where (•) (c) represents the variables under contingency set N C. As a part of the security objective function, the voltage stability index of the system V s is then given by:

The other part of the security objective function adds the number of overloaded lines N lo as a penalty function to avoid transmission overload under normal operation and after contingencies, which is expressed as:

where N loi is the number of overloaded lines under normal operation and N loi (c) is the number of overloaded lines after contingencies. 

Combined with ( 9 ) and ( 10 ), the security objective function S e is calculated as:

where σ is the penalty parameter, indicating punishment intensity with a large value.
section_name: Objective Function

section_content: The integration of IPFC will introduce new constraints to the system. To construct the constraints of the system and IPFC, an equivalent IPFC model is required first. IPFC compensates for multiple transmission lines simultaneously by insertion of series voltage source converters (VSC). The series converters can be equivalent to additional injected power to nodes based on the power injection method, and the IPFC equivalent power injection model with two series branches is shown in Figure 1. 

tives, which are referred to as the "Master" converters. The other converter is referred to as the "Slave" and controls DC voltage to keep the stability of the converters. In Figure 1, the "Master" line of the system is transmission line ij, and the "Slave" line is transmission line ik. Nodes m and n are added as additional nodes for subsequent model derivation and power flow calculation. Pinj_i, Qinj_i, Pinj_m, Qinj_m, Pinj_n, and Qinj_n are equivalent injected power generated by converters at buses i, m, and n, respectively. Xseim and Xsein are equivalent impedances of the series output voltages. Vi∠θi, Vj∠θj, and Vk∠θk are voltages of buses i, j, and k, respectively. Rij, Rik, Xij, Xik, Bcij, and Bcik are resistance, reactance, and admittance against the ground of transmission lines ij and ik, respectively. 

Equivalent IPFC model with two series branches based on the power injection model. 

(
section_name: Constraints

section_content: Combined with IPFC equivalent injected power, the constraints of power balance can be modified as:

where Nb is the total number of buses, including additional buses m and n. Nodes i and j are any nodes of the system, and i, j∈Nb. Pgi and Qgi are the active and reactive power, respectively, of the generator at bus i. Pdi and Qdi are the active and reactive loads, respectively, at bus i. Pinj_i and Qinj_i are the active and reactive injected powers, respectively, at bus i. Pinj_i and Qinj_i are corresponding values related to IPFC injected power, and Pinj_i = 0, Qinj_i = 0 when buses i and j are not directly connected with IPFC. Gij and Bij are the conductance and susceptance, respectively, of transmission line ij. θij is the phase angle difference between nodes i and j. 

Power balance after N -1 contingency should also be satisfied as:

IPFC converters are connected to series transformers respectively on the AC side, and all converters are linked together at their DC terminals. Active power can be exchanged through the common DC link of VSCs. Therefore, accurate and flexible power flow control of multiple transmission lines can be realized at the same time. To maintain stability of Converters selected by the operator are then given priority to achieve control objectives, which are referred to as the "Master" converters. The other converter is referred to as the "Slave" and controls DC voltage to keep the stability of the converters. In Figure 1, the "Master" line of the system is transmission line ij, and the "Slave" line is transmission line ik. Nodes m and n are added as additional nodes for subsequent model derivation and power flow calculation. P inj_i, Q inj_i, P inj_m, Q inj_m, P inj_n, and Q inj_n are equivalent injected power generated by converters at buses i, m, and n, respectively. X seim and X sein are equivalent impedances of the series output voltages. V i ∠θ i, V j ∠θ j, and V k ∠θ k are voltages of buses i, j, and k, respectively. R ij, R ik, X ij, X ik, B cij, and B cik are resistance, reactance, and admittance against the ground of transmission lines ij and ik, respectively. 

(1) Equality Constraints Combined with IPFC equivalent injected power, the constraints of power balance can be modified as:

where N b is the total number of buses, including additional buses m and n. Nodes i and j are any nodes of the system, and i, j∈N b. P gi and Q gi are the active and reactive power, respectively, of the generator at bus i. P di and Q di are the active and reactive loads, respectively, at bus i. P inj_i and Q inj_i are the active and reactive injected powers, respectively, at bus i. P inj_i and Q inj_i are corresponding values related to IPFC injected power, and P inj_i = 0, Q inj_i = 0 when buses i and j are not directly connected with IPFC. G ij and B ij are the conductance and susceptance, respectively, of transmission line ij. θ ij is the phase angle difference between nodes i and j. 

Power balance after N -1 contingency should also be satisfied as:

IPFC converters are connected to series transformers respectively on the AC side, and all converters are linked together at their DC terminals. Active power can be exchanged through the common DC link of VSCs. Therefore, accurate and flexible power flow control of multiple transmission lines can be realized at the same time. To maintain stability of the common DC bus voltage, the balance of active power exchange between converters of IPFC must be guaranteed and the sum of active power flowing through the DC link of converters must be zero. Ignoring the active power loss of IPFC, the active power conservation between converters can be deduced as:

For the master converter of IPFC, there are four control modes, including impedance compensation control mode (ICCM), phase regulation control mode (PRCM), voltage regulation control mode (VRCM), and constant power control mode (CPCM) [6]. ICCM can match the existing capacitive line compensation in the system and controls the equivalent impedance, PRCM can be used for phase shifting and controls the voltage phase angle difference between connected buses, VRCM changes and controls the magnitude of the bus voltage, and CPCM schedules power flow directly and controls active and reactive power of the transmission line. The control mode variable X∈{1, 2, 3, 4} is incorporated into the optimization model for the optimal control mode, and corresponding settings are obtained as the optimization results as well. 

X = 1 represents that the master converter works in ICCM. The equivalent impedances Z im and Z im (c) are intended to be equal to the set control target Z ref as follows:

X = 2 represents that the master converter works in PRCM. The voltage phase angle differences θ mi and θ mi (c) between bus m and bus i, respectively, are intended to be equal to the set control target θ ref, and the voltage amplitude at bus m is intended to stay the same as that of bus i. These are given by:

(16) X = 3 represents that the master converter works in VRCM. The voltage amplitudes V m and V m (c) at bus m are intended to be equal to the set control target V ref, and the voltage phase angle at bus m is intended to stay the same as at bus i and can be obtained as:

X = 4 represents that the master converter works in CPCM. The active power and reactive power P mj, Q mj, P mj (c), and Q mj (c) of the transmission line mj are intended to be equal to the set control targets P mjref, Q mjref, respectively, as follows:

The slave converter needs to control the DC bus voltage, and only one control dimension remains. There are merely two modes: constant active power mode and quadrature compensating voltage control mode, and the former mode is generally selected for the emphasis of active power in actual system operation [6]. The active powers P nk and P nk (c) of the transmission line nk are intended to be equal to the set control target P nkref and are given by:

Energies 2024, 17, 1660

It is intended that thermal stability be realized after PSCOPF optimization, so in addition to converters, the generator output plan is intended to remain unchanged before and after N -1 contingency. The constraint should be satisfied as:

where V gi is the voltage amplitude of the generator at PV bus i. 

(2) Inequality Constraints In the power system, voltage amplitudes of the buses and output active power of the generators are required to satisfy the upper and lower bounds:

where V i min and V i max are the minimum and maximum magnitudes, respectively, of voltage amplitude at any bus and P gi min and P gi max are the minimum and maximum magnitudes, respectively, of generator output active power. 

Additionally, IPFC converters are required to meet the corresponding inequality constraints as follows:

where V se, V se, and θ se are the converter output AC voltage, voltage amplitude, and phase angle, respectively; I iσ and I iσ are the current and current amplitude through the converter, respectively; and P dc max, V se max, and I iσ max are the maximum interactive active power, output AC voltage amplitude, and output AC current of the converter, respectively. 

Control modes of IPFC converters correspond to the control mode variable X and are analyzed in the proposed optimization model. In the optimization, the optimized objective function and control setting under each control mode are calculated, and the final optimized result of the proposed model is the optimum of all control modes in which the optimal control mode and setting values are included.
section_name: 1) Equality Constraints

section_name: Solution Method for the PSCOPF Model Considering IPFC Control Modes

section_content: The PSCOPF model above is a complex nonlinear optimization model that is hard to solve by the traditional method. The PSO algorithm is easy to implement and is applied to solve the model. Objective functions are set as the fitness function in the PSO algorithm. The optimization variables include IPFC control modes, IPFC control targets, voltage amplitudes of PV buses, and generator output active power. The equality constraints are satisfied in power flow calculation, and the inequality constraints are treated by introducing the penalty function. 

Iterative schemes of converter output voltages and equivalent injected power for different IPFC control modes are derived respectively and used in the power flow calculation. The power and voltages required in the proposed model can then be obtained by power flow calculation. The calculated power and voltages in the objective function and constraints are further used to obtain the optimal result of the proposed optimization model. Figure 2 shows the process of solving the proposed model. During the solving process, power flow iterative calculation is greatly related to IPFC control modes and is illustrated in detail in the next section.
section_name: Framework of the Solution Method

section_content: The power flow calculation method considering IPFC control modes can generally be divided into two parts. The first part transforms converter output AC voltages into injected power via the power injection method and updates the voltages based on reference values of the controlled variables and iteration values of state variables [33]. The other part calculates the power flow using general methods, such as Newton-Raphson, when iteration values of IPFC injection power are available. The iterative update of converter output AC voltages is crucial, which realizes the variable exchange between the two calculation parts. Figure 3 shows the overall flow chart for power flow calculation considering IPFC, and iterative schemes for each control mode are subsequently derived. In the IPFC equivalent power injection model in Figure 1, series converters inject power into the common bus i and additional buses m and n equivalently. Then, the power flow of lines ij and ik are adjusted. Equivalent injected power generated by the converters at buses i, m, and n can be calculated as: In the IPFC equivalent power injection model in Figure 1, series converters inject power into the common bus i and additional buses m and n equivalently. Then, the power flow of lines ij and ik are adjusted. Equivalent injected power generated by the converters at buses i, m, and n can be calculated as:

Output voltages of series converters in the injected power calculation in ( 23 )-( 28 ) depend on IPFC control modes, and the iterative updates of series output voltages during power flow calculation are derived as follows. 

(1) Iterative Updates of Series Output Voltages for the Master Converter When the master converter works in ICCM, equivalent impedance Z im and current I im between bus i and bus m can be calculated as follows:

ε is the accuracy of the convergence. After the kth iteration of power flow calculation, if the criterion of convergence Z re f -Z (k) im < ε is not met, the iteration process will continue. The converter output AC voltage can be updated as:

Then, the equivalent injected power can be calculated based on (25) and (26). After the (k + 1)th iteration of power flow calculation with the (k + 1)th injected power, values of state variables can be updated, and the (k + 1)th value of equivalent impedance Z im between bus i and bus m can be calculated by ( 29 ) and (30). In the iteration process, the equivalent impedance Z im will tend to the control target Z ref. 

When the master converter works in PRCM, the voltage phase angle difference θ mi between bus m and bus i is intended to be equal to the control target θ ref, and the voltage amplitude at bus m is intended to keep the same as bus i. After the kth iteration of power flow calculation, if the criterion of convergence θ re fθ (k) mi < ε is not met, the voltage difference between bus m and bus i can be modified to:

Energies 2024, 17, 1660

The converter output AC voltage can be updated with values of state variables can be updated as:

Calculation of the current I im between bus i and bus m is the same as (30) with kth values of state variables. Then, the iteration process of power flow calculation is consistent with ICCM. The controlled variable, the voltage phase angle difference θ mi, will tend to the control target θ ref. 

When the master converter works in VRCM, the voltage amplitude V m at bus m is intended to be equal to the control target V ref, and the voltage phase angle at bus m is intended to keep the same as that at bus i. After the kth iteration power flow calculation,

< ε is not met, the voltage difference between bus m and bus i can be modified to:

The converter output AC voltage can be updated with values of state variables after the kth iteration, and the calculation of V seim ∠θ seim is the same as (34). The iteration process is consistent with that of ICCM. The controlled variable, the voltage amplitude V m at bus m, will tend to the control target V ref. 

When the master converter works in CPCM, the active power and reactive power P mj, Q mj of the transmission line mj are intended to be equal to the control targets P mjref and Q mjref, respectively. Power of the transmission line mj can be obtained as follows:

After the kth iteration of power flow calculation, if the criteria P mjre f -

mi < ε are not met, the converter output AC voltages after PQ decomposition V seimp and Vseimq can be deduced from ( 36 ) and (37), respectively, as follows:

Therefore, the converter output AC voltage can be updated as:

Then, the (k + 1)th equivalent injected power, values of state variables and power of the transmission line can be successively calculated. In the iteration process, the controlled variables-the active power and reactive power P mj and Q mj, respectively-will tend to the control targets P mjref Q mjref, respectively. (2) Iterative Updates of Series Output Voltages for the Slave Converter For transmission line nk with the slave converter, the active power is intended to be equal to the control target P nkref in constant active power control mode. After the kth iteration of power flow calculation, if the criterion of convergence P nkre f -P (k) nk < ε is not met, active and reactive components of the converter output AC voltage amplitudes V seinq, Vseinp after the kth iteration can be obtained successively by ( 14 ) and (36), respectively:

The output AC voltage of the slave converter can be updated as:

The equivalent injected power of the master converter and the slave converter can be obtained according to the specific control modes. The PSCOPF model considers the four control modes of the IPFC master converter, and power flow calculation of the system under a specific IPFC control mode can be realized. The optimal result of the PSCOPF model is the optimal scheduling strategy of the system with IPFC.
section_name: Power Flow Calculation in the PSCOPF Model Solution Process

section_content: In order to verify the optimization results of the proposed model and corresponding solution method, the overall PSCOPF model is compared to three scenarios, including the original system, optimization not considering preventive-security-constrained control, and optimization not considering IPFC control modes. Numerical calculation and analysis are conducted based on the IEEE 39-bus New England case, as shown in Figure 4. Considering the imbalanced power flow distribution of Line 2-3 (active power is nearly 320 MW) and Line 2-1 (active power is nearly 160 MW), the master converter and slave converter of IPFC are placed on Line 2-3 and Line 2-1, respectively, to realize a better power flow distribution.
section_name: Case Studies

section_content: Equivalent impedance X seim and X sein of the IPFC series output voltages are taken as X seim = X sein = 0. 04 p. u. The contingency set N C in the security objective function includes Lines 3-4, 3-18, 4-14, 6-7, 10-13, 15-16, 16-21, 17-18, 21-22, 23-24, and 28-29. PSO algorithm is applied to solve the PSCOPF model, with parameters of IPFC and generators optimized jointly. Related coefficients to be set mainly include inertia coefficient w and acceleration coefficients c 1 and c 2. The inertia coefficient w controls the momentum of the particle by weighing the contribution of the previous velocity to ensure convergent behavior and optimal tradeoff of exploration and exploitation. Acceleration coefficients c 1 and c 2 control the stochastic influence of the cognitive and social components on the overall velocity of a particle, where c 1 expresses how much confidence a particle has in itself and c 2 expresses how much confidence a particle has in its neighbors. Referring to [34, 35], which have similar optimization problems, the inertia coefficient w is set as 0. 729, and acceleration coefficients c 1 and c 2 are both set as 1. 4962. The number of particle swarms is 50, the number of iterations is 150, weight coefficient λ in the objective function is 20,000, and the penalty parameter σ in the security objective function is 5000. 

In order to verify the optimization results of the proposed model and corresponding solution method, the overall PSCOPF model is compared to three scenarios, including the original system, optimization not considering preventive-security-constrained control, and optimization not considering IPFC control modes. Numerical calculation and analysis are conducted based on the IEEE 39-bus New England case, as shown in Figure 4. Considering the imbalanced power flow distribution of Line 2-3 (active power is nearly 320 MW) and Line 2-1 (active power is nearly 160 MW), the master converter and slave converter of IPFC are placed on Line 2-3 and Line 2-1, respectively, to realize a better power flow distribution.
section_name: Optimization Results of the PSCOPF Model Considering IPFC Control Modes

section_content: Equivalent impedance Xseim and Xsein of the IPFC series output voltages are taken as Xseim = Xsein = 0. 04 p. u. The contingency set NC in the security objective function includes Lines 3-4, 3-18, 4-14, 6-7, 10-13, 15-16, 16-21, 17-18, 21-22, 23-24, and 28-29. PSO algorithm is applied to solve the PSCOPF model, with parameters of IPFC and generators optimized jointly. Related coefficients to be set mainly include inertia coefficient w and acceleration coefficients c1 and c2. The inertia coefficient w controls the momentum of the particle by weighing the contribution of the previous velocity to ensure convergent behavior and optimal tradeoff of exploration and exploitation. Acceleration coefficients c1 and c2 control the stochastic influence of the cognitive and social components on the overall velocity of a particle, where c1 expresses how much confidence a particle has in itself and c2 expresses how much confidence a particle has in its neighbors. Referring to [34, 35], which have similar optimization problems, the inertia coefficient w is set as 0. 729, and acceleration coefficients c1 and c2 are both set as 1. 4962. The number of particle swarms is Based on the coefficients set above, the program is implemented in MATLAB R2018b. After 150 iterations, the result shows that the optimal control mode is ICCM, and the corresponding control target is Z ref = -0. 03462 + j 0. 01887p. u. The control target of the slave converter is set as P nkref = 0. 7843p. u. 

The conventional FACTS model ignores control characteristics and regards injected voltages unchanged before and after N -1 contingencies [36]. However, IPFC usually works in constant power control mode in actual operation, where the injected voltage would vary when contingencies occur. The unreasonable assumptions in the conventional model may lead to security problems. In order to reveal the potential risk of the conventional model and indicate the ability and reliability of IPFC power flow regulation with the proposed PSCOPF method, respective optimization results are compared. Table 1 shows the numerical results in the four corresponding circumstances, and the corresponding objective function values of the system without IPFC are calculated based on the original state parameters. 

Table 1 shows that the security of the system has been significantly improved through the proposed PSCOPF model, including voltage stability and the risk of overload. The conventional strategy seems to eliminate the risk of overload, but these results are obtained under the unrealistic assumption that the injected voltage remains unchanged before and after contingencies. When the conventional strategy is under realistic operation, i. e., the constant power control mode, the risk of overload is inevitable, indicating the potential security problem of the unreasonable assumption in the conventional model. Additionally, both the operation cost and voltage stability of the conventional model under realistic operation deteriorated compared to the proposed model. Figure 5 shows that there are risks of overload in the actual operation of the system without IPFC, while IPFC helps to regulate power flow distribution and improve the security of the system. Load violations can be eliminated after N -1 contingencies under the proposed PSCOPF model and unrealistic conventional model. However, ignorance of IPFC control characteristics is unreasonable. Overload appears under the realistic operation, which may even lead to cascading failure in certain extreme situations. The proposed model considers the actual operation properties of IPFC, and the system is secure under all N -1 contingencies, proving the feasibility and effectiveness of the proposed model.
section_name: Optimization Results of the PSCOPF Model Considering IPFC Control Modes

section_content: To verify the importance of considering IPFC control modes, a PSCOPF model that does not consider IPFC control modes is also compared. The objective function, constraints, contingency set NC and corresponding optimization parameters keep consistent with the proposed PSCOPF model considering IPFC control modes. The only difference is that the master converter of IPFC uses the fixed and most common control mode, CPCM. After implementation in MATLAB, the optimization results of the two models are compared, as shown in Table 2.
section_name: Comparison with a Fixed IPFC Control Mode

section_content: To verify the importance of considering IPFC control modes, a PSCOPF model that does not consider IPFC control modes is also compared. The objective function, constraints, contingency set NC and corresponding optimization parameters keep consistent with the proposed PSCOPF model considering IPFC control modes. The only difference is that the master converter of IPFC uses the fixed and most common control mode, CPCM. After implementation in MATLAB, the optimization results of the two models are compared, as shown in Table 2. 

Table 2 shows that due to the priority to satisfy the elimination of load violations, operation cost is slightly sacrificed when considering IPFC control modes. However, the PSCOPF model with the fixed CPCM control mode cannot satisfy the prerequisite to eliminate the risk of overload, and the voltage stability of the system is also worse. Therefore, it is meaningful to consider IPFC control modes in the PSCOPF optimization. In Figure 6, CPCM corresponds to the fixed IPFC control mode. The optimal control mode corresponds to the optimization results considering four control modes of the master converter and ICCM is the calculated optimal result in the tested system. ICCM reveals better performance in heavy load lines and helps to eliminate overload after N -1 contingencies. However, load violation exists after the N -1 contingency of line 21-22 under the fixed IPFC control mode, dissatisfying the security constraints in actual operation. The security problem can be addressed by changing to the optimal IPFC control mode without the need for additional equipment costs, proving the superiority of the proposed model considering IPFC control modes. 

The power and voltages required in the objective function and constraints of the proposed model are obtained through power flow calculation. Iterative updates of converter output voltages and equivalent injected power under each IPFC control mode are utilized to obtain the required parameters and further solve the optimization model. The optimal solution of the proposed model can decrease operation cost and improve the security of the system.
section_name: Comparison with a Fixed IPFC Control Mode

section_content: IPFC control modes are considered in the PSCOPF model to improve security-constrained OPF. Security is given priority in the paper, and potential overloads are avoided through the proposed model. Dynamic voltage stability as well as the transition duration of IPFC control mode switching can be considered in the optimization to avoid avalanchelike voltage decrease after emergency disturbances. After considering the dynamic process, field experiments can be conducted in hardware and then in the real power system.
section_name: Future Research

section_content: A PSCOPF model considering IPFC control modes is established in this paper. It rec- In Figure 6, CPCM corresponds to the fixed IPFC control mode. The optimal control mode corresponds to the optimization results considering four control modes of the master converter and ICCM is the calculated optimal result in the tested system. ICCM reveals better performance in heavy load lines and helps to eliminate overload after N -1 contingencies. However, load violation exists after the N -1 contingency of line 21-22 under the fixed IPFC control mode, dissatisfying the security constraints in actual operation. The security problem can be addressed by changing to the optimal IPFC control mode without the need for additional equipment costs, proving the superiority of the proposed model considering IPFC control modes. 

The power and voltages required in the objective function and constraints of the proposed model are obtained through power flow calculation. Iterative updates of converter output voltages and equivalent injected power under each IPFC control mode are utilized to obtain the required parameters and further solve the optimization model. The optimal solution of the proposed model can decrease operation cost and improve the security of the system.
section_name: Conclusions

section_content: IPFC control modes are considered in the PSCOPF model to improve security-constrained OPF. Security is given priority in the paper, and potential overloads are avoided through the proposed model. Dynamic voltage stability as well as the transition duration of IPFC control mode switching can be considered in the optimization to avoid avalanche-like voltage decrease after emergency disturbances. After considering the dynamic process, field experiments can be conducted in hardware and then in the real power system.
section_name: Future Research

section_content: A PSCOPF model considering IPFC control modes is established in this paper. It reconciles the economy and security of the system, and the preventive-security-constrained technique is involved in achieving IPFC control potential after contingencies. Constraints of the proposed model generally include power system operation constraints and IPFC operation constraints, in which IPFC control characteristics under different control modes are fully analyzed. The optimal selection of IPFC control modes and control parameters can then be obtained to improve economy and security of the system. Numerical results prove the superiority of the proposed model considering IPFC control modes, avoiding overloads under unrealistic assumption of the constant converter output voltages or fixed control modes.
section_name: Conclusions

section_content: Funding: This work was supported by Science and Technology Project of State Grid Corporation of China "Planning technologies and applications to improve the flexibility of regional power grid structures with high proportions of new energy consumption" ( 5100-202318036A-1-1-ZN ).

section_name: Data Availability Statement:

section_content: The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.
section_name: Data Availability Statement:

section_name: Conflicts of Interest:

section_content: Author Hui Cai was employed by the company State Grid Jiangsu Electric Power Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
section_name: Conflicts of Interest: