1. Introduction: MAED with Renewable Energy
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  , a multi area multi objective dynamic economic dispatch (MAMODED) for a retailer, a model taking into account hydrothermal, wind and a power reserve market is proposed. Scenario-based method (SBM) is adopted for uncertainty modeling while optimality condition decomposition (OCD) is the proposed solution tool. The proposed approach demonstrates the real-time scheduling of joint thermal systems with undispatchable wind energy in a three-area and five-area multi area system. In  , a penalty function-hybrid direct search method (PF-HDSM) is formulated to handle a multi-area problem with wind and thermal sources. In these two works however, security and stability aspects, and solar cost functions are ignored. In  , security and control areas of MAMODED with wind and solar REs were addressed. Co-operative optimization strategy with RE source (RES) curtailment, dispatch-redispatch strategies and RES storage operation are considered. However, the stability analysis was not done. Further, quadratic cost functions and HVAC tie lines are considered.
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.
2. Effects of HVDC Tie Lines on Stability
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  : 1) Since the interconnection is a synchronous tie, frequency disturbance in one area is transferred to the other area(s). 2) Power swing in one area affects the other area(s) in the interconnected system since large power swings may result in frequency tripping. Further, major faults in any of the areas network may lead to complete failure of the whole multi area system. 3) There is increased fault level if an existing AC area is interconnected with another AC area(s) by HVAC tie line. This is because the additional parallel lines reduce the equivalent reactance of the interconnected system.
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  : 1) Since the HVDC tie lines between the areas is asynchronous, the two areas can either operate at the same frequency or may have different frequencies. This facilitates interconnection of AC systems at different frequencies to enable them operate independently and to maintain their frequency standards 2) HVDC link provides fast and reliable control of magnitude and direction of load flow by controlling the firing angle of the converters. The rapid control of power flow increases the limit of transient stability in the multi area system. 3) Further, power swings in the interconnected AC networks can be damped rapidly by modulating the power flowing through the HVDC line. Thus the stability of the multi area system is increased. 4) Fault level in each dynamic area remains unchanged.
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   : 1) It is more flexible and economical 2) Frequency oscillations in the areas can be damped quickly 3) There is inherent overload capability. Hence transient stability of the interconnected system is increased without an increase in installed capacity. Thus, the overall stability of the multi area system is improved. 4) Heavily loaded AC area networks can be reinforced by using MTDC systems
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.
3. Scmamoded Formulation
MAMODED with RE and MTDC tie line was formulated in  . In this paper, the stability aspects of the same problem are investigated. The security constrained MAMODED (SCMAMODED) is formulated as
where M is the total number of areas and W1 is the weighting factor between the fuel cost and emissions. The fuel cost function is defined as
where T is the operation time, N, W and S are the total thermal, wind and solar units respectively and , and are the corresponding fuel cost functions at time t and scenario s .
It is worth to note that is the outage cost defined by
where and 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 
where is the MTDC tie line from area m to k, is the transmission cost coefficient relevant to and T is the vector of real power transmission between areas.
The emissions objective function, E, is formulated as
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)
・ Real power balance (RPB)
where and are control variables used in adjusting the power output of the wind and solar units respectively, is the interconnected system load at time t.
・ Reactive Power Balance(RPB)
2) Generation Units Constraints (GUC)
・ Real power generation constraints
・ Reactive power generation constraints
3) System Security Constraints (SSC)
・ Load shedding limits
・ Multi Scene power constraints
・ Ramping speed of thermal generating units constraints
4) System Stability Constraints (SSC)
・ Line Flow Constraint (MTLFC)
where the number of is tie lines and is the active power flow in the tie line t
・ Converter Tap Ratio Constraint (MCTRC)
・ Converter Ignition Angle Constraint (MCIAC): Facilitates fast and reliable control of power flows
・ Converter extinction Angle Constraint (MCEAC): Facilitates rapid control to increase transient stability limit
・ Current Constraint (HCC)
・ HVDC Voltage Constraint (HVC)
・ 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
where and represents the MTDC tie line transmission capability
4. Proposed Methodology
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.
4.1. Probabilistic Load Flow (PLF)
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    . The PLF involves buses such as PQ (negative load), PV (voltage controlled), PX and RX. Numerical and Analytical methods are applied in solving the PLF.
AC PLF Model for RES
The PLF model is given by
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, and .
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.
4.2. Scenario-Based Method (SBM)
For a multivariate general function, 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, is generated for describing the probable values of X such that  ;
where 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
where is the total cost paid by the retailer, is the probability of scenario s, is the purchased power from reserve market in time t, scenario s and area m, is the price of energy purchased from the power reserve in time t, scenario s ($/MWh) and is the production cost of the ith thermal unit located in area m in time t. The objective function of a rational retailer that is to be maximized is defined by
where the price of energy is sold to customers in time t, and is the load demand at time, scenario s and area m.
4.3. The Solution Process
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)
5. Results and Analysis
There are six thermal sources in each area with cubic cost coefficients and emission coefficients (for SO2, NOX and CO2). The necessary wind and solar power parameters, the renewable power that is available at different scenarios of the study are as in  . The area generation capacity for wind and solar 300 MW and 250 MW respectively. The test case consists of five interconnected areas. The one line diagram of the 5-area MAMODED is as shown in  except that the power sources in each area are as shown in Table 1. The area m1 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 Table 2.
Table 1. Power sources in the areas.
Table 2. MTDC tie line flow limits.
5.1. Frequency Control
Voltages observed at the chosen node for the MTDC and HVAC are as shown in Figure 1. The voltage angles for two chosen nodes for the MTDC tie line flow limits are as depicted in Figure 2. In order to fully exploit advantages of MTDC interconnection technology, MTDC links monitoring a five-area system frequency at both ends of the link are considered. This is as shown in Figure 3 the real and reactive power flows are varied as per the frequency deviation between the interconnecting nodes. This implies that MTDC tie line ends responds to thermal and RES loss in one area as in the other areas, leading to a change in power flow from the faulted point. This change in exchanged power from the
Figure 1. Voltage at the HVAC and MTDC links.
Figure 2. Voltage angles (dig) at the observed MTDC node.
Figure 3. Frequency difference at both ends of MTDC link.
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  . Further the MTDC interconnection ensures a fast damping of the frequency oscillations in between the areas   .
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  . Technically speaking, the action of the Converter Ignition Angle (CIA) which facilitates fast and reliable control of power flows and the Converter Extinction Angle (CEA) which facilitates rapid control of power flow so as to increase transient stability limit makes frequency control possible.
The simplified MTDC control characteristic is shown in Figure 4. From this characteristic, it can be seen that ±2.0 pu represents maximum/minimum power ratings of standard MTDC links. In order to provide protection to the healthy interconnected system, if the frequency in the adjacent area falls below a certain value, for example 49.50 or 49.65 Hz, then further control support actions of MTDC links by the CIA and CEA must be suspended at the present values. The vital advantage of the MTDC link control is that all the control processes are localized in individual MTDC tie line controllers. The control actions are very fast since there are no complex mathematical computations involved. In addition, there is no need of additional information exchange system, complex system control centers or sensitive network data sharing topology  .
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
Figure 4. MTDC control droop characteristic  .
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.
5.2. Power Transfer at the Interface
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 Figure 5.
From Figure 5, it can been that immediately the MTDC tie line transferring 2500 MW is faulted and thus tripped, the power flow control center adjusts to compensate for the deficiency in power wheeling and increases the probabilistic load flows (PLF) on the remaining MTDC lines. Due to the flexible nature of the MTDCs and the presence of RES and storage a total recovery of the deficit is realized, unlike the case in the traditional HVAC interconnection. However, one very important merit of MTDC technology comes to light here, that is, the MTDC tie lines can be overloaded. The inherent overload capability of the
Figure 5. Active power transfer at Scenario 1.
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)  . The results are presented in Figure 6 and Figure 7. From these results, it can be seen that as soon as the frequency in one area of the interconnected system start falling as a result of generation loss, all MTDC links automatically increases their power wheeling capacities (import/export) so as to offer maximum support to the endangered interconnected system even it mean MTDC overload. Since two of these MTDC lines are operating close to their power limits, there is no increase of the power exchange to maximum, for there a practical possibility of overloading the flexible HVDC interconnection. The MTDC damps the frequency and power swings in the areas rapidly by modulating its power flows. The preset values are more accurately obtained by easily tuning the MTDC control and stability parameters.
In this paper SCMAMODED with all the possible security and stability constraints
Figure 6. Active power transfer at Scenario 2.
Figure 7. System frequency at Scenario 2.
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.
 Soroudi, A. and Rabiee, A. (2013) Optimal Multi-Area Generation Schedule Considering Renewable Resources Mix: A Real-Time Approach. Generation, Transmission & Distribution, IET, 7, 1011-1026.
 Chen, C.-L., Chen, Z.-Y. and Lee, T.-Y. (2014) Multi-Area Economic Generation and Reserve Dispatch Considering Large-Scale Integration of Wind Power. International Journal of Electrical Power & Energy Systems, 55, 171-178.
 Kahl, M., Freye, C. and Leibfried, T. (2015) A Cooperative Multi-Area Optimization with Renewable Generation and Storage Devices. Power Systems, IEEE Transactions on, 30, 2386-2395.
 Musau, M.P., Odero, A.N. and Wekesa, C.W. (2016) Multi Area Multi Objective Dynamic Economic Dispatch with Renewable Energy and Multi Terminal DC Tie Lines. 6th IEEE International Conference on Power Systems (ICPS 2016), New Delhi, 4-6 March 2016.
 Paensuwan, N. and Yokoyama, A. (2009) Risk-Based TTC Calculation of a Power System with Renewable Energy Resources. IEEE Power Tech Conference, Bucharest, 28 June-2 July 2009, 1-8.
 UCTE, IPS/UPS (2008) Synchronous Interconnection of the IPS/UPS with the UCTE: Summary of Investigations and Conclusions.
 Kim, C., Sood, V.K., Jang, G., Lim, S. and Lee, S. (2009) HVDC Transmission, Power Conversion Applications in Power Systems. John Wiley & Sons (Asia) Pte Ltd, Asia.