_{1}

Security constrained multi area multi objective dynamic economic dispatch (SCMAMODED) with renewable energy (RE) and all the possible MTDC stability constraints is formulated for the first time. The stability merits of multi terminal DC (MTDC) tie lines as compared to the traditional HVAC forms the main objective of this paper. Probabilistic load flow (PLF) is applied to determine the system parameters while the uncertainties are modelled using Scenario Based Method (SBM). The simulation results reveal that with the use of MTDC tie lines, the frequency and voltage stability in the MAMODED with renewable energy sources (RES) are enhanced while keeping the MTDC power exchange interface nodes at secure levels.

Multi area stability is the continuance of intact operation of an interconnected system following a disturbance. Security evaluation demands stability analysis. This paper is concerned with stability studies in a multi area system with increased renewable energy (RE) penetration levels. Three recent works have considered multi area economic dispatch (MAED) with renewable energy (RE). In [

Paper organization: This paper is outlined as follows; Section II highlights the effects of HVDC tie lines on stability, Section III is the problem formulation, Section IV is the methodology, Section V contains results and analysis and finally Section VI is the conclusion.

When two dynamic areas are interconnected by an HVAC tie line, it is necessary that there should be sufficiently close frequency control on each of the areas. Such an interconnection is called synchronous interconnection or synchronous tie. Thus HVAC tie line provides a rigid connection between any two areas. Even with coordinated control of the AC areas, the operation of the HVAC tie line suffers from three technical problems [

The use of HVDC tie lines provides several security advantages to the interconnected power system. HVDC interconnections or HVDC tie lines provide a loose coupling between the two interconnected areas. This is called non-synchronous (asynchronous) tie line. This interconnection has four merits [

A multi terminal DC (MTDC) has more than two converter stations and interconnecting DC transmission lines. The Radial MTDC technology of HVDC is the choice for asynchronous tie lines in this paper because of the following reasons [

Contributions: Five objective SCMAMODED with RE and MTDC tie lines has been considered for the first time. More accurate cubic cost functions have been used for the thermal and emissions functions. MTDC stability constraints have been formulated alongside security constraints.

MAMODED with RE and MTDC tie line was formulated in [

min f = ∑ m = 1 M { W 1 F + ( 1 − W 1 ) E } (1)

where M is the total number of areas and W_{1} is the weighting factor between the fuel cost and emissions. The fuel cost function is defined as

F = ∑ t = 1 T [ ∑ i = 1 N F ( P i , t , s ) + F ( O C t ) ] + ∑ j = 1 W F ( P j , t , s ) + ∑ k = 1 S F ( P k , t , s ) + F ( MTDC ) (2)

where T is the operation time, N, W and S are the total thermal, wind and solar units respectively and F ( P i , t , s ) , F ( P j , t , s ) and F ( P k , t , s ) are the corresponding fuel cost functions at time t and scenario s .

It is worth to note that F ( O C t ) is the outage cost defined by

F ( O C t ) = ∑ s = 0 S [ b Δ P L , t , s + c Δ Q L , t , s ] (3)

where Δ P L , t , s and Δ Q L , t , s are the amounts of real and reactive load-shedding in scene s at time t, b and c are the constants of real and reactive load-shedding costs. Unlike the operating cost, the outage cost determined based on the risk in each scene.

The transmission cost in MAMODED for power transfer between areas using MTDC tie lines can be expressed as [

F ( MTDC ) = ∑ m = 1 M − 1 ∑ k = m + 1 M f m k T m k (4a)

where T m k is the MTDC tie line from area m to k, f m k is the transmission cost coefficient relevant to T m k and T is the vector of real power transmission between areas.

The emissions objective function, E, is formulated as

E = ∑ l = 1 G E ( P i , t , s , P j , t , s , P k , t , s ) ; G = N + W + S (4b)

where G is the total number of generators in the system

The SCMAMODED is solved subject to the following system, security and stability constraints

1) Power Balance Constraints (PBC)

g b ( P b , Q b ) = 0 (5)

・ Real power balance (RPB)

∑ m = 1 M { ∑ i = 1 N P i , t , s + ∑ j = 1 W n w , t P j , t , s + ∑ k = 1 S n s , t P k , t , s + Δ P L , t , s } = ∑ t = 1 T P D , t (6)

where n w , t ∈ [ 0 , 1 ] and n s , t ∈ [ 0 , 1 ] are control variables used in adjusting the power output of the wind and solar units respectively, P D , t is the interconnected system load at time t.

・ Reactive Power Balance(RPB)

∑ m = 1 M { ∑ i = 1 N Q i , t , s + ∑ j = 1 W n w , t Q j , t , s + ∑ k = 1 S n s , t Q k , t , s + Δ Q L , t , s } = ∑ t = 1 T Q D , t (7)

2) Generation Units Constraints (GUC)

g U ( P i , t , s , P j , t , s , P k , t , s ) ≤ 0 (8)

・ Real power generation constraints

P i , t , s , m min ≤ P i , t , s ≤ P i , t , s , m max (9a)

P j , t , s , m min ≤ P j , t , s ≤ P j , t , s , m max (9b)

P k , t , s , m min ≤ P k , t , s ≤ P k , t , s , m max (9c)

・ Reactive power generation constraints

Q i , t , s , m min ≤ Q i , t , s ≤ Q i , t , s , m max (10a)

Q j , t , s , m min ≤ Q j , t , s ≤ Q j , t , s , m max (10b)

Q k , t , s , m min ≤ Q k , t , s ≤ Q k , t , s , m max (10c)

3) System Security Constraints (SSC)

g S ( P i , t , s , P j , t , s , P k , t , s , S Y ) ≤ 0 (11)

・ Load shedding limits

P D , t ≥ Δ P L , t , s , m ≥ 0 (12a)

Q D , t ≥ Δ Q L , t , s , m ≥ 0 (12b)

・ Multi Scene power constraints

P i , t , s , m − P i , t , s − 1 , m ≥ 0 (13)

・ Ramping speed of thermal generating units constraints

P i , max d o w n ≤ P i , t , s − 1 , m − P i , t − 1 , s , m ≤ P i , max u p (14a)

P i , max d o w n ≤ P i , t , s , m − P i , t − 1 , s − 1 , m ≤ P i , max u p (14b)

4) System Stability Constraints (SSC)

| g s s ( h v d c ) | ≤ 0 (15a)

・ Line Flow Constraint (MTLFC)

| P t | ≤ P t , max , t = 1 , 2 , 3 , ⋯ , N t

where N t the number of is tie lines and P t is the active power flow in the tie line t

・ Converter Tap Ratio Constraint (MCTRC)

T min ≤ T ≤ T max (17)

・ Converter Ignition Angle Constraint (MCIAC): Facilitates fast and reliable control of power flows

α min ≤ α ≤ α max (18)

・ Converter extinction Angle Constraint (MCEAC): Facilitates rapid control to increase transient stability limit

γ min ≤ γ ≤ γ max (19)

・ Current Constraint (HCC)

I d c , min ≤ I d c ≤ I d c , max (20)

・ HVDC Voltage Constraint (HVC)

V d c , min ≤ V d c ≤ V d c , max (21)

・ MTDC Transmission Capacity Limits (MTCL): The transfer including both generation and reserve from area m to k should not exceed the MTDC tie line transfer capacities for stability and security considerations. This can be expressed as

T m k , min ≤ T m k + R m k ≤ T m k , max (22)

where T m k , min and T m k , max represents the MTDC tie line transmission capability

In this section, we consider the probalistic load flow(PLF) for the MAMODED with RES, uncertainty and variability modelling of the RE sources using Scenario-based method (SBM) and a summary of the solution process of the SCMAMODED.

Load flow methods are classified into two; Deterministic and Probabilistic methods. In this paper, probabilistic load flow (PLF) is applied since it is capable of handling the uncertain (intermittent) RE. PLF requires inputs with probability density function (PDF) or cumulative density function (CDF) which are the distribution functions for wind speed and solar radiation intensity respectively [

AC PLF Model for RES

The PLF model is given by

w = f ( x ) , Δ x = J 0 − 1 Δ w = S 0 Δ w (23a)

z = g ( x ) , Δ z = G 0 J 0 − 1 Δ w = T 0 Δ w (23b)

The following steps are involved [

Step 1: Read the interconnected system, wind, and Solar PV data

Step 2: Run the deterministic load flow (DLF) using newton rap son (NR) method. This will help in obtaining expected values of nodal voltages, line flows, S 0 and T 0 .

Step 3: Determine the cumulative values of thermal and RES generation and loading according to their respective PDFs where applicable.

Step 4: By computation, determine the cumulative of the generated active power, absorbed reactive power, power injections and the state variables (∆x and ∆z).

Step 5: Obtain the PDF and CDF of the respective ∆x and ∆z.

For a multivariate general function, y = F ( X ) where X is a vector containing the uncertain input values, the SBM uncertainty modelling is a method for finding the expected value of y. A set of scenarios, Ω s is generated for describing the probable values of X such that [

y = ∑ s ∈ Ω s π s F ( X s ) (24)

where π s is the probability of state s.

In an interconnected system, uncertainties and variability in RES power generation, reserve market prices and the load profile of the system, results into a probabilistic MAMODED which is formulated in this paper. The total cost of energy procurement is given by

C T = ∑ S , t , m π s P P s , m ( t ) λ s ( t ) + ∑ i , t , m C i m ( P i m ( t ) ) (25)

where C T is the total cost paid by the retailer, π s is the probability of scenario s, λ s ( t ) is the purchased power from reserve market in time t, scenario s and area m, λ s ( t ) is the price of energy purchased from the power reserve in time t, scenario s ($/MWh) and C i m ( P i m ( t ) ) is the production cost of the i^{th} thermal unit located in area m in time t. The objective function of a rational retailer that is to be maximized is defined by

F = ∑ S , t , m π s P D , s , m ( t ) λ c ( t ) − C T (26)

where λ c ( t ) the price of energy is sold to customers in time t, and P D , s , m ( t ) is the load demand at time, scenario s and area m.

Step 1: Get the prediction data for the RE power, load and other known conditions using PLF

Step 2: Set the range of the RE power fluctuation and the risk constant(s) of the load shedding

Step 3: Build the SCMAMODED model

Step 4: Invoke the power system stability checks (MTDC constraints etc.)

Step 5: Solve the problem using Modified PSO (MPSO)

There are six thermal sources in each area with cubic cost coefficients and emission coefficients (for SO_{2}, NO_{X} and CO_{2}). The necessary wind and solar power parameters, the renewable power that is available at different scenarios of the study are as in [_{1} is a central area with a thermal source and a RE reserve, it is connected to the other four areas through tie lines. The MTDC tie-lines’ flow limits are as shown in

Area | Objectives | % RE |
---|---|---|

m_{1} | Thermal(T), Reserve(R) | 16 |

m_{2} | Thermal(T), Wind(W) | 33 |

m_{3} | Thermal(T), Solar(S) | 33 |

m_{4} | Thermal(T), Wind(W), Solar(S) | 67 |

m_{5} | Thermal(T) | 0 |

Overall Penetration | 30 |

MTDC Tie Lines | Power (MW) |
---|---|

T_{1-3} | 2350 |

T_{1-}_{4} | 2300 |

T_{1-}_{5} | 2350 |

T_{2}_{-3} | 2300 |

T_{3}_{-}_{5} | 2350 |

T_{2}_{-}_{4} | 2350 |

Voltages observed at the chosen node for the MTDC and HVAC are as shown in

adjacent area(s) results in a corresponding frequency change in that area also. The regulating power is obtained from the neighboring area(s), regardless of whether thermal or RES are involved. The result is that RES and thermal generators in that area would respond, like in traditional case of HVAC interconnection. However, the vital difference is that the MTDC tie lines do not put the other area(s) at risk and therefore limits the effect of the disturbance on the other areas and the entire system to an acceptable level [

A frequency control curve is set such that the MTDC tie lines do not function for all low level contingencies (i.e. the dead band zone). In case of severe faults, actions of the asynchronous MTDC are in direct correlation with frequency values. If the frequency in the affected area approaches a first under-frequency RES-load shedding value, MTDC tie lines offers all the possible help. This is made possible by increasing the flexible tie line power flows to the maximum levels that are still safe for the neighboring areas with intermittent power sources (RES) and RES storage [

The simplified MTDC control characteristic is shown in

HVDC and HVAC interconnections have their own separate advantages and all the merits of both interconnections are combined in MTDC tie lines and preserved. The multi area system operate in an asynchronous manner and

during the disturbances, and thus, the areas help each other. If the risk of severe faults and their consequences are constrained in a practical manner to one area only then this is practically possible. When the frequencies interconnected areas are restored to the preset values, then the active power exchanged over the MTDC lines returns to its original set values as frequency oscillation in the areas can be quickly damped. This means that real multi area system separation should never happen in a practical real time system.

The only scenario when power exchange in between any two areas is fully cut off in HVAC interconnection is that when the areas are simultaneously in extremely severe contingency state so that none of the two areas is able to help each other either by import or export of power. However in this case, the presence of RES and the flexible HVDC (Radial MTDC) will help alleviate this situation to a great extent.

In this paper, two scenarios of power transfer at the MTDC interface have been considered. During the healthy functioning of the MTDC interconnection, the power flows are monitored and controlled from a centralized point. In the event that one MTDC is tripped, the stressed multi area multi objective system responds to compensate for the power loss by increasing in a flexible manner the power flow in the remaining healthy MTDCs so as to maintain the preset power levels unaltered .The scenario 1 has been simulated and results are presented in

From

MTDC increases the transient stability of the interconnected system without any increase in the installed capacity. Thus the overall stability of the MAMODED with RE is improved.

Scenario 2 considers MTDC tie lines with interconnected area control algorithms (IACA) [

In this paper SCMAMODED with all the possible security and stability constraints

has been formulated with RES put into consideration. The uncertainties have been modelled using scenario based method (SBM) and a probalistic load flow (PLF) used in parameter determination. From the simulated results, it is clear that MTDC interconnection is preferred to HVAC due to their stability merits even in presence of RES. Further with control strategies the merits of both HVDC and HVAC interconnections can be realized using MTDC tie lines. This paper is concerned only about multi area stability (continuance of intact operation following a disturbance); the security aspects of the interconnected system (degree of risk in the ability to survive imminent disturbance or contingencies) are proposed as an area of further research. Further, environmental emission levels need to be included in the security-stability analysis.

The author declares no conflicts of interest regarding the publication of this paper.

Musau, M.P. (2018) Stability Analysis of a Multi Area System with Renewable Energy and Multi-Terminal DC Tie Lines. Journal of Power and Energy Engineering, 6, 85-96. https://doi.org/10.4236/jpee.2018.69009