Research Article  Open Access
Stendley Busan, Muhammad Murtadha Othman, Ismail Musirin, Azah Mohamed, Aini Hussain, "A New Algorithm for the Available Transfer Capability Determination", Mathematical Problems in Engineering, vol. 2010, Article ID 795376, 30 pages, 2010. https://doi.org/10.1155/2010/795376
A New Algorithm for the Available Transfer Capability Determination
Abstract
This paper presents a fast and accurate method to determine the available transfer capability. Ralston's method is used to predict the two trajectory points of voltage magnitude, power flow, and maximum generator rotor angle difference. Then, the cubic spline interpolation technique is used to accurately trace the PV, PS, or P curves between two points of trajectory. The PV, PS and P curves represent as the variations of voltage magnitude, power, flow and maximum generator rotor angle difference due to the increase of power transfer. The actual available transfer capability value is determined at the intersection point between the curve and the constraints limit. The effectiveness of the proposed method is verified by referring to the results of ATC for a case study of 2737Polish system and 39New England bus system. The proposed method gives satisfactorily accurate and fast computation of ATC as compared to recursive AC power flow method.
1. Introduction
Transferring an electric power from one place to another is an alternative way to provide effective electric power required by the demand. This may assist towards reduction in a system operational cost. Nowadays, the power trade activity which is involved in the wholesale power market requires accurate information of power transfer between areas. Such vital information can help power marketers, sellers, and buyers in planning, operation, and reserving transmission services [1]. There are two significant indices in the transfer capability assessment, namely, the total transfer capability (TTC) and the available transfer capability (ATC). TTC represents as the maximum amount of power that can be transferred over the interconnected transmission network in a reliable manner while meeting all of a specific set of defined pre and postcontingency system conditions [2]. On the other hand, ATC is a measure of the additional amount of power that flows across the interface, over and above the base case flows without jeopardizing power system security [3].
The determination of ATC for a large and complex power system usually utilizes excessive amount of computational time. This instigates to a new development of a fast and accurate method in determining the ATC value. Various approaches have been proposed to determine ATC such as using the methods of DC power flow [1], AC power flow [4], optimal power flow [5], sensitivity [6], curve fittingbased cubic spline interpolation technique [7], and artificial neural network [8]. The method based on linear DC power flow considering distribution factors is considered fast but less accurate for transfer capability analysis because the DC network model does not require the voltage magnitude and reactive power component in the power flow calculation. Therefore, the computation is based on the linear DC power flow resulting in an inaccurate ATC value, especially for the heavily stressed system that is caused by critical contingencies. The AC power flow method gives an accurate solution in determining the ATC because it considers the effects of reactive power flows and voltage limits. However, transfer capability evaluation using repetitive AC power flows is time consuming because it requires a load flow solution at every transfer step size. To avoid many repetitive AC power flow solutions, curve fitting technique such as cubic spline interpolation technique has been used [7]. There are various curve fitting techniques that are used for voltage stability analysis such as the least square fit of secondorder polynomial [9], cubic spline interpolation [10] and quadratic approximation [11].
This paper proposes a new approach to determine fast and accurate value of ATC by using Ralston's method incorporating with cubic spline interpolation technique. The Ralston's method is categorized under the secondorder RungeKutta method and this is similar to Heun's and midpoint methods. However, the Ralston's method is superior to Heun's and midpoint methods in terms of providing a minimum bound of truncation error in extrapolation [12]. The Ralston's method is used to determine the two trajectory points of voltage magnitude, power flow, or maximum generator rotor angle difference. Then, the cubicspline interpolation technique is used to accurately trace the PV, PS, or P curves between the two trajectory points of voltage magnitude , power flow or maximum generator rotor angle difference (), respectively. The PV, PS and P curves represent as the variations of voltage magnitude, power flow, and maximum generator rotor angle difference due to the increase of power transfer, respectively. The ATC is then determined at a point when the voltage magnitude limit, power flow limit or generator rotor angle difference limit intersects the curve. In the ATC estimation, the transmission line and voltage magnitude limits are considered as the steadystate security constraints whereas the maximum generator angle difference limit is referred to as the transient stability constraint. The transfer capability of a system is analyzed under two different sets of transfer, which are the areatoarea ATC and pointtopoint ATC. Areatoarea ATC is the additional amount of power transferred from the selling area to the buying area without jeopardizing the system security. On the other hand, pointtopoint ATC is the additional amount of power transferred from the selling bus to the buying bus without violating the system security. The effectiveness of the proposed method in estimating fast and accurate computation of ATC is verified on the case studies of 2737bus Polish system [13] and 39New England bus system. The proposed method can be used in the probabilistic assessment of transfer capability. This is due to the fact that the proposed method is able to accurately determine the ATC in a less computational time for every system operating condition. The system operating conditions such as the transmission line failures are usually generated by the Monte Carlo simulation technique or parametric bootstrap technique. The ATCs are then used in the probabilistic based risk or uncertain assessment of transfer capability.
2. Problem Formulation
The first section describes the problem definition of ATC followed by the explanation of Ralston’s method that is used to determine the two trajectory points of voltage magnitude, power flow, and maximum generator rotor angle difference. The last section provides a detailed explanation of cubic spline formulation that is used for tracing the PV, PS, and P curves between the two trajectory points for accurate ATC determination.
2.1. ATC Problem Definition
ATC is defined as the TTC less than the transmission reliability margin (TRM), less than the sum of existing transmission commitments (ETCs), and capacity benefit margin (CBM) [2, 14, 15]. The TRM is the amount of transmission capability necessary to ensure that the interconnected system is secure under a reasonable range of uncertainties in system conditions. The CBM is the amount of transmission transfer capability reserved by load serving entities to ensure access to generation from interconnected systems to meet generation reliability requirements. The ETC is the normal transmission flows included in the given case. The methods to determine the TRM, CBM, and ETC margins may vary among regions, power pools, individual system, and loadserving entities.
ATC must satisfy certain principles balancing both technical and commercial issues, so that the interconnected transmission network is performed based on the commercial requirements associated with transmission service requests. The following principles identify the requirements for the calculation and application of ATC.(a)Electricity demand and supply cannot be treated independently of one another. All system conditions must be considered to accurately access the capabilities of the transmission network.(b)Electric power flows resulting from each power transfer use the entire network and are not governed by the commercial terms of the transfer. (c)ATC calculations must use a regional or widearea approach to capture the interactions of electric power flows among individual, regional, subregional, and multiregional systems.(d)The determination of ATC must accommodate reasonable uncertainties in system conditions and provide operating flexibility to ensure a secure operation of the interconnected network.
In the determination of ATC, the transmission lines flow and voltage magnitudes limits have to be taken into account in the calculation. All these limits can be handled by the AC load flow power system model. Limits due to transient or oscillations are not often addressed in the ATC determination because these limits are crudely approximated by flow limits [16]. However, the large disturbance such as system faults, loss of generator, or equipment outages could lead to undesirable behavior that affects the stability of a system. The undesirable behavior is associated with the transient stability which could lead to great losses and costly to the utilities. Therefore, it is necessary to consider the transient stability constraints within the ATC calculation.
2.2. Formulation of Ralston’s Method
A generic ATC computation is performed by solving recursive AC power flow calculations due to the increased amount of power transfers between areas or buses. The ATC is then determined by referring to the increase amount of power transfer caused to the violation of a system constraint such as the voltage magnitude limit, transmission line limit or generator rotor angle difference limit. In the ATC computation using the recursive AC power flow solution, the variations of voltage magnitude , MVA power flow and maximum generator rotor angle difference () due to the increase of MW power transfers can be described in terms of PV, PS, and P curves, respectively. By considering the PV, PS and P curves as the quadratic polynomial form, Ralston's method can be used to approximate the two trajectory points of voltage magnitude, power flow, and maximum generator rotor angle difference. Then, the cubicspline interpolation technique is used to accurately trace the PV, PS, or P curves between the two specific points of trajectory. The proposed methodology is able to provide accurate value of ATC by taking into account the steadystate and transient stability constraints.
The Ralston’s method is used to approximate the two trajectory points of voltage magnitude, power flow, and maximum generator rotor angle difference [12], and it is derived from a basic extrapolation equation that is given by (2.1) Equation (2.1) represents that the is used to extrapolate from an old value of to a new value of over a distance of . The second order of (2.1) gives where The , , and are the unknown constants used to satisfy the three conditions By referring to (2.2a), the value of is assumed to be 2/3 thus resulting in the values of and This yields to a Ralston’s method given by (2.3). where, Note that is equivalent to the power transfer, . It is worth mentioning that the Ralston’s method given in (2.3) is used to extrapolate the secondorder polynomial curvature. The secondorder polynomial is represented by The first order of (2.5) yields (2.6a) and it is representing as , Thus, in (2.4a) is given by By applying (2.6a) into (2.4b), The value of constants and can be determined by using the least square method [12]. The step size, h, of power transfer is determined as below, where is the number of incremental steps for power transfer. In this case study, is specified as 4 in which it is reasonable enough to provide accurate approximation of the two trajectory points. is the initial power transfer that is 1 MW. is the lookahead power transfer which may cause the violation of voltage, transmission line, or generator rotor angle difference limits. The methodology of lookahead power transfer is explained elaborately in Section 3.
The next power transfer, for each th incremental step is determined by using.The values obtained from (2.6b), (2.6c), (2.6d), and (2.7a) are used in (2.3) so that the Ralston’s method could perform the extrapolation in order to obtain the two points of trajectory. The value of y that is obtained from (2.3) represents as the voltage magnitude, power flow, or maximum generator rotor angle difference. The is increased at each th incremental step by using (2.7a) until the Ralston’s method in (2.3) gives y value that violates the system constraint. The last two values of y represent as the two trajectory points of voltage magnitude, power flow, or maximum generator rotor angle difference. It is then used in the cubicspline interpolation technique to accurately trace the PV, PS or P curves.
It is obvious that the step size, of power transfer given in (2.6d) is highly dependent on the number of incremental steps, . A reasonable number of incremental steps, , need to be specified so that the Ralston's method in (2.3) could perform the extrapolation with a minimum computational time. The number of incremental steps, , is specified under two categories. First, a large number of incremental steps, , may cause a lengthy computational time in the extrapolation process due to a small step size, . Second, a small number of incremental steps, , yield a fast computational time in the extrapolation process due to a large step size, . In this case study, four incremental steps are used, and it is reasonable enough for the Ralston's method to provide fast and accurate approximation of the two trajectory points. Furthermore, the cubicspline interpolation technique is used to accurately trace the curve between the two trajectory points with a large value of step size, , and this will be discussed in the next subsection.
2.3. CubicSpline Interpolation Technique
The methodology of cubicspline interpolation technique is basically based on determining the four known points and then fitting appropriate curves to the four points. In the cubicspline technique [7], tracing the curves , and begins with finding the value for parameters and which are given as The values for parameters and are used in the cubicspline equations in order to obtain the curve functions of ), and , which are given as,
In the PV curve fitting, the parameters and can be described as , and , respectively. On the other hand, the parameters and can also be described as , and , respectively for the case of PS curve fitting. For the case of P curve fitting, the parameters and can be described as , and , respectively. is the cubicspline function that is used for tracing the curves of voltage magnitude, , MVA power flow, , and maximum generator rotor angle difference, (). is the increase of power transfer by 1 MW between and . is the number of three incremental steps, that is, 1, 2, and 3. Specifically, is used for tracing the curves between the four points of with respect to the increase of k_{l} by 1 MW from to . Whereby, f () represents as the four points of voltage magnitude, , MVA power flow, , or maximum generator rotor angle difference, , which are obtained from the AC power flow solutions. The four points of real power transfer, can also be described as , where = 1, 2, 3, and 4. For an example, the curve from point to point is traced by using with the increase of k_{1} by 1 MW from x_{1} = 1 MW to x_{2} = 300 MW.
2.4. Determination of PV, PS, and P Curves Using CubicSpline Interpolation Technique
Generally, there are two main procedures involved in the PV curve fitting using the cubicspline interpolation technique. First, the voltage at each point of real power transfer, , is obtained by solving the AC power flow solution. Second, the cubicspline interpolation technique is used for tracing the voltage curves, , based on the four voltage points, which are obtained from the previous four AC power flow solutions and this is shown in Figure 1. Particularly, the curve from point to point is traced by using with the increase of by 1 MW from to . Then, the next curve from point to point is traced by using with the increase of by 1 MW from to . Finally, the last curve from point to point is traced by using with the increase of by 1 MW from to .
Similarly, the cubicspline interpolation technique is used in tracing the PS curves representing as the MVA power flow variations with respect to the increase of power transfer. The procedures that are used in tracing the curves , and between the four points of MVA power flow are similar to those described for tracing the voltage curves, except that the voltage variables in (2.8) to (2.11)are replaced by the MVA power flow variables. This is similar to a case whereby the cubicspline interpolation technique is used for tracing the P curves.
2.5. Transient Stability Constraint
In this study, the transient stability is obtained by analyzing the “first swing” of each generator. Transient stability is referred to as the generator rotor angle that is returning to its synchronism state after the fault is cleared. A classical model of a synchronous generator is used in this case of study, and the details can be accessed in [17]. The transient stability, based rotor angle is measured by referring to the difference between relative rotor angle, with respect to the center of inertia (COI) [18–20]. The transient stability limit should be less or equal to 180°, and it is given by where is rotor angle of th generator. is rotor angle difference of gth generator. is generator inertia constant in seconds for gth generator.
The maximum rotor angle difference for unstable condition is taken at the end of time simulation. This is due to the fact that the relative rotor angle is monotonically increasing if the generator losing its synchronism. On the other hand, for a stable condition, the maximum relative rotor angle is taken within the simulation time interval. This is because the increased relative rotor angle is returning back to its synchronism state after the fault is cleared.
3. Procedure of ATC Evaluation Using Ralston’s Method Incorporating CubicSpline Interpolation Technique
Generally, the main steps involved in the transfer capability computation are the definition of a base case, determination of network response, and finding the maximum transfer or ATC. Determination of the areatoarea and pointtopoint ATCs using the Ralston’s method incorporating cubicspline interpolation technique is described as follows.
(a)Establish a solved base case AC power flow solution. (b)Specify the area or point of transfers. For the pointtopoint transfer, a generator is considered as a selling bus and a load is a buying bus. However, the areatoarea transfer considers participation of all generators in the specified selling area and all loads in the specified buying area. (c)Simultaneously, increase the power generation and load at the selected buses or areas at three incremental steps in order to obtain the variations of voltage , the power flow and maximum generator rotor angle difference (). Where, is the bus number, is the transmission line connected between bus and bus and is the number of generator. The AC power flow solution should be performed for each incremental step of and The amount of power generation, is equivalent to the amount of power transfer, Then, the sensitivity method is used to identify the sensitive bus, transmission line or generator that has the highest potential to be violated due to the increase amount of power transfer [21]. Equations (3.1), (3.2), (3.3) and (3.4) are the sensitivity methods that are used to approximate the amount of power transfer, P, corresponding to each bus, transmission line and generator. Then, the sensitive line, bus or generator is selected based on the minimum amount of power transfer: where and are the linear estimation of power transfer based on the violations of lower voltage limit, upper voltage limit, thermal limit, and generator rotor angle difference limit, respectively. and are the lower and upper voltage limits which are 0.9 p.u. and 1.1 p.u., respectively. is the transmission line limit. is the generator rotor angle difference limit which is 180°. is the power transfer for every th incremental step. is incremental steps. , and are the voltage magnitude at each bus, i, power flow at each transmission line, ij, and maximum rotor angle difference at each generator bus, g, respectively.The sensitive bus, i, transmission line, ij, or generator, g, is selected based on the minimum value of power transfer amongst or . represents as the minimum value of power transfer given by. (d)Determine the lookahead power transfer based on the sensitive bus, transmission line or generator. The methodology that is used to determine the lookahead power transfer is initially derived from the formulation of secondorder polynomial that is given in [22]. Further derivation of the firstorder quadratic formulation in (2.6a) yields to, Equation (3.7) is obtained by substituting (3.6) into (2.5), By expanding (3.7), Equation (3.8) is derived to become Whereby, By deriving (3.9), In this case, is the system parameter constraint such as the transmission line rating, lower limit of voltage magnitude that is 0.9 p.u., upper limit of voltage magnitude which is specified at 1.1 p.u., or generator rotor angle difference limit specified as 180°.The lookahead power transfer, , is determined by using (3.12) which is derived from (2.6a). In (3.12), the is determined by using (3.11) It is obvious that the lookahead power transfer formulation in (3.12) is similar to (3.6). The y_{o} and which are obtained by using (3.10) and (3.11), respectively, are used in (3.12) to determine the . The values of , and are calculated by using the least square method [12] that utilizes the V_{n}, S_{n} or _{n} at three incremental steps of power transfer, P_{n}. The values of V_{n}, S_{n} or _{n} are determined by referring to the sensitive bus, transmission line or generator obtained from procedure (c). This shows that the is determined by referring to the sensitive bus, transmission line or generator.(e)Use the Ralston’s method in (2.3) to determine the two trajectory points of voltage magnitude, power flow or maximum generator rotor angle difference. This refers to the sensitive bus, transmission line or generator obtained from procedure (c). Initially, the AC power flow solution is performed at three incremental steps of with = . This is performed to obtain the variation of V_{n}, S_{n} or at the sensitive bus, transmission line or generator. Then, the V_{n}, S_{n} or _{n}, and P_{n} are used in the least square method [12] to determine the new values of , , and . The is specified as the last power transfer of P_{n} so that , , and are determined at stable system condition. Hence, accurate estimation of two trajectory points could be obtained by using the Ralston’s method that takes into account the , and . Specifically, the , , and are used in (2.6b), (2.6c), (2.6d) and (2.7a) so that the Ralston’s method in (2.3) is able to determine the two trajectory points of voltage magnitude, power flow or maximum generator rotor angle difference. The is not an optimum value of ATC. Therefore, the Ralston’s method is used to extrapolate at the two trajectory points for optimum or accurate determination of ATC.(f)Use the cubicspline to trace the PV, PS or P curve between the two points of trajectory. The PV, PS or P curve is determined based on the sensitive bus, transmission line or generator.(g)Determine areatoarea or pointtopoint ATCs which are the maximum power transfer obtained when the voltage limit, the MVA line rating, or generator rotor angle difference limit intersects the PV, PS or P curve at two trajectory points, respectively.
The above procedures are summarized in terms of flowchart shown in Figure 2. The Ralston’s method incorporating with cubicspline interpolation technique gives a faster ATC computation which implies less AC power flow solutions as compared to the ATC computation method using the recursive AC Newton Raphson power flow solutions [4, 23].
4. Results and Discussion
The performance of the Ralston's method incorporating with cubicspline interpolation technique that used in the determination of ATC is verified in terms of accuracy and computation speed. CPU timing for the transfer capability analysis was obtained using 2.4 GHz, Intel Core 2 Duo with 1 GB of memory. The 2737bus Polish power system is used as a test case to illustrate the determination of ATC using the proposed technique. The system is comprising of 6 areas namely, area 1, area 2, area 3, area 4, area 5, and area 6. The 2737bus system is modelled with 193 generation units, 2544 load units and 3506 lines. In this study, the upper and lower voltage limits are assumed to be 1.1 p.u. and 0.9 p.u., respectively. The thermal limit is also used as a system constraint in the ATC computation. However, the generator rotor angle difference limit of 180° is not considered in the ATC computation for the 2737bus Polish power system. This is due to the fact that detailed information of generating unit is not available in order to compute ATC by considering the generator rotor angle difference limit.
Nevertheless, the generator rotor angle difference limit as well as the transmission line limit and voltage magnitude limit is considered as the constraints of ATC computation for a case study of 39New England bus system. The system is consisting of 10 generation units, 29 load units and 46 transmission lines [24]. The system data is given in Tables 7, 8, and 9. The transmission line limit information is taken from [25]. The system is comprised into three areas namely area 1, area 2 and area 3 as illustrated in Figure 3.
4.1. Faulted Bus and Tripping Line Selection
The transient stability analysis is performed to ensure that the system is operating in a secure manner without violating the generator rotor angle limit during the occurrence of fault. In this case study of transfer capability assessment, it is assumed that a threephase fault is occurring at a particular transmission line. The faulted bus is referring to as the nearest bus which is connected to a faulty line [17]. Therefore, the faulty line should be tripped in order to clear the fault so that a stable generator rotor angle could be obtained during power transfer.
4.2. Fault Critical Clearing Time and Final Simulation Time Selection
In a transient stability analysis, the faulty line should be tripped at a certain fault critical clearing time and this criterion does affect the stability of generator rotor angle. Therefore, a set of relays and protecting circuits should operate within the fault critical clearing time so that the fault is cleared without causing any loss of synchronism for the generators [26, 27]. However, the determination of fault critical clearing time is not considered in this case study of transfer capability assessment. Nevertheless, the fault critical clearing time is set at a typical maximum allowable time so that the generator rotor angle is stable during the occurrence of fault. In [28], an analysis to estimate the fault critical clearing time has been conducted on the 39New England bus system. In conjunction to this system, the estimated fault critical clearing time for the 10 machines is best to be within the range of 0.13 to 0.24 second. Therefore, in this case study, the fault critical clearing time of 0.15 second is chosen to clear the fault. It is selected based on the fact that the fault is expected to be cleared before reaching the end of fault critical clearing time. Thus, the fault critical clearing time of 0.15 second is viable enough to ensure the stability of generator rotor angle.
The rotor angles with respect to the COI reference frame of all generators are initially increasing or decreasing until a peak value is reached. Then, the rotor angle starts returning to its stable equilibrium point and it is said to be the first swing stable. On the other hand, a system is said to be the first swing unstable if the postfault angle is increasing or decreasing monotonically for at least one of the machines [20, 28]. In this case study, duration for the simulation time is within the range of until seconds, and it is chosen as to analyze the stability of first swing generator rotor angle difference [20, 28].
4.3. ATC Results Using Ralston’s Method Incorporating CubicSpline Interpolation Technique: Case Study of 2737Bus Polish Power System
Prior to the ATC determination, the cubicspline interpolation technique is used to trace the PS or PV curves between the two points of trajectory determined by the Ralston’s method. The ATC is then determined by referring to the maximum power transfer that causes the limiting levels of MVA power flow or voltage magnitude intersects the PS or PV curves, respectively. Detailed information of generating unit is not available for the computation of ATC computation considering the generator rotor angle difference limit. In this case study, the load bus 2737 is chosen to describe the determination of ATC using the proposed method. The load bus 2737 is a sensitive bus that limits the MW transfer from area 2 to area 6.
By referring to Figure 4, the sensitivity method given in (3.1) linearly determines the amount of ATC, = 1243 MW, that violates the lower limit of voltage magnitude. Linear approximation of power transfer considered in the sensitivity method usually gives inaccurate value of ATC especially for a large system. In the AC power flow solution, the amount of = 1243 MW could be very large, which may cause instability to the system condition. Hence, accurate nonlinear estimation of ATC could be obtained by considering the quadratic form of PV or PS curves. In this case study, the sensitivity method is used to determine the sensitive bus or transmission line which is based on the minimum amount of power transfer as given in (3.5). The sensitive bus or transmission line has the potential to cause the first violation of voltage magnitude or thermal limits due to the increase of power transfer, respectively. Hence, fast and accurate estimation of ATC could be determined by referring only to the sensitive bus or transmission line. In this case study, the sensitive load bus 2737 limits the power transfer between area 2 to area 6, and it is shown in Figure 4. In Figure 4, it is observed that the = 524 MW does not exceed the power transfer at the second trajectory point of voltage magnitude. Hence, the system condition is stable when = 524 MW is below the second trajectory point. The is not an optimum or accurate value of ATC, and it is determined by using (3.12). However, the unstable system conditions may occur when power transfer exceeds the second trajectory point. Therefore, the value of 524 MW is used in the Ralston's method to accurately extrapolate at the two trajectory points for optimum or accurate estimation of ATC. Figure 4 shows four incremental steps performed of power transfer by using the Ralston's method.
The cubicspline interpolation technique is then used to fit the PV curve between the two trajectory points at the sensitive bus 2737. The two trajectory points are at the power transfer of 462 MW and 720 MW. The power transfer at the two trajectory points is then used in the cubicspline interpolation technique for PV curve fitting. In Figure 5, it is shown that the cubicspline interpolation technique traces the PV curve based on the four points of voltage magnitude. The four points of voltage magnitude are obtained from the AC power flow solutions considering the power transfers of 462 MW, 548 MW, 634 MW and 720 MW. A point is noted where voltage limit of 0.9 p.u. intersects the PV curve. This point yields an actual ATC value of 699 MW.
4.4. Results of AreaToArea ATC and PointToPoint ATC: Case Study of 2737Bus Polish Power System
Tables 1 and 2 represent the results of the areatoarea ATC and the pointtopoint ATC, respectively. The ATCs obtained from the Ralston’s method incorporating with cubicspline interpolation technique are compared with the ATCs obtained from the recursive AC power flow method. The comparisons are made in terms of accuracy and time taken in computing the ATC.


The results shown in Tables 1 and 2 indicate that the ATC obtained is due to the overloaded line. For instance, by referring to Table 1, the ATC is 914 MW for transfer case between areas 1 and 2, and it is obtained due to the overloaded line 713–449. The overloaded line that limits the increase of power transfer occurs for the transfer case between buses. For an example, by referring to Table 2, the ATC is 140 MW for the transfer case between buses 98 and 2636, and it is obtained due to the overloaded line 2507–2513. Simulations that have been carried out on the test system indicate that the ATCs are determined not due to the violation of voltage limit.
The recursive AC power flow method is a basic approach to determine accurate value of ATC. This means that the proposed method is able to compute accurate value of ATC since it is similar to the result determined by the recursive AC power flow method, and it is shown in Tables 1 and 2. Both methods are able to calculate accurate value of ATC since they considers nonlinear condition of reactive power flows and voltage magnitudes. The proposed method is better than the linear DC power flow method in terms of accuracy for the ATC computation. In terms of computational time, it is noted that the proposed method computes a much faster ATC value as compared to the recursive AC power flow method. This is due to the fact that the proposed method does not perform many recursive load flow solutions in the ATC determination. Finally, the results have shown that the proposed method is able to provide accurate value of ATC with less computational time.
4.5. Results of AreaToArea ATC and PointToPoint ATC: Case Study of 39New England Bus System
Tables 3 and 4 present the result of areatoarea ATCs and pointtopoint ATCs for a test system of 39New England buses, respectively. The ATCs are determined by considering the systems constraints of transmission line limit, voltage stability limit, and generator rotor angle difference limit. Whereby, the rotor angle difference limit is referred to as the transient stability limit. In this case study, it is assumed that a single (N1) contingency type of threephase fault happened at bus 10. A threephase fault yields the most severe fault current compared to the other types of unsymmetrical fault. The computation of ATCs due to transient stability limit is performed by considering the tripping of line 10–13 for clearing the fault. The fault critical clearing time of second, is specified for the tripping at line 10–13. The selection of fault critical clearing time, second has been explained elaborately in Section 4.2. The transient response of generator rotor angles is monitored for 1.5 seconds in the case study of ATC that takes into account the transient stability limit. Since, the tripping of line 10–13 is taken into account in the determination of ATC considering the transient stability limit the tripping of line 10–13 is considered as a single (N1) contingency of three phase fault. On the other side, a double (N2) contingency may sometimes occur, which would be critical to the system operating conditions especially during the power transfer. Indeed, this is an intriguing issue that needs be considered for further analysis on the impact of N2 contingencies on the transfer capability assessmentbased transient stability limit. In the transient stability analysis, a small time step size of 0.03 second is used in the timedomain of rotor angle. The small time step size is used to discard the higherorder terms of Taylor series expansion used in the Euler's method. Therefore, the error of rotor angle approximation is decreased for every successive point of time domain [17].


In Table 3, the minimum interarea ATC of 64 MW is obtained for the transfer case from area 1 to area 2 and the transfer case from area 2 to area 3 yields a maximum interarea ATC value of 350 MW. On the other hand, by referring to Table 4, the minimum pointtopoint ATC of 41 MW is obtained for the power transfer from bus 32 to bus 24. The power transfer case from bus 34 to bus 26 yields a maximum pointtopoint ATC value of 353 MW. For both cases of power transfer, the ATCs are obtained based on the violation of transmission line limit. It is obvious that the proposed method provides relatively similar results of ATC compared to the recursive AC power flow method. In terms of computational time, it is obvious that the proposed technique gives a fast ATC calculation compared to the recursive AC power flow method. This is because the Ralston's method incorporated with the cubicspline interpolation technique executes fewer numbers of power flow solution in the ATC determination in comparison to the recursive AC power flow method.
4.6. Performance Comparison of Ralston’s Method Incorporating with CubicSpline Interpolation Technique at Various Numbers of Steps in Power Transfer: Case Study of 39New England Bus System
Table 5 presents the performance of areatoarea ATC computation based on several numbers of steps, n, applied in the proposed method. The same case study as discussed in Section 4.5 is used in this analysis. Basically, the power transfer step size, h, is specified depending on the number of steps, n, applied in (2.6d). The performance of the proposed technique is evaluated based on four difference numbers of steps, n. In addition, the percentage of relative error in applying the proposed method to determine the ATC is also computed for every transaction in the case study. By comparing the ATC values obtained in Tables 3 and 5, the interarea transactions computed by using the proposed technique with gives the range of relative error percentage within 0.36% to 3.85%. Most of these transactions are 2 MW higher than the ATC values obtained by using the recursive AC power flow solution. This is referring to the case study of interarea power transfer from selling area 1 to buying area 2, selling area 2 to buying area 1, and selling area 3 to buying area 1. The highest percentage of relative error of 3.85% is obtained from selling area 1 to buying area 3. Furthermore, the performance of the proposed method in ATC computation is evaluated based on and incremental steps of power transfer. It is observed that the proposed method with and gives relatively similar results as compared to the ATCs obtained by using the recursive AC power flow solution as shown in Table 3. This indicates that it is important to use and in the Ralston’s method for accurate extrapolation at the two trajectory points. Therefore, it will give relatively accurate value of ATC compared to the result determined by the Ralston’s method with . Besides that, the Ralston's method with the number of steps specified at and also gives less percentage of relative error at 1.28% and 0.29% from the selling area 1 to buying area 3 and selling area 2 to buying area 3, respectively. In Table 5, it is observed that the proposed method with gives similar results compared to the ATCs obtained by using the recursive AC power flow solution as shown in Table 3. This is due to the fact that as the number of steps, , is increased, the power transfer step size, , will become narrow. This may assist towards a better result of curve fitting performed by the cubicspline interpolation technique. Thus, accurate estimation of ATC result can be obtained.





In terms of computational time, the Ralston’s method with, , and is utilizing the same computational time of 1.8 seconds in determining the ATC. However, it is much faster than the Ralston’s method with that takes 0.05 minutes in determining the value of ATC. It is obvious that and improves the performance of Ralston’s method which gives similar results of ATC with less computational time. By referring to the abovementioned detail explanations, the proposed method with is used in the other power transfer case study for a fast and accurate estimation of ATC.
4.7. Determination of ATC due to the Transient Stability Limit
This section discusses the variations of maximum generator rotor angle difference that is used in (3.4) to determine the sensitive generator. The sensitivity method given in (3.4) is used to identify the generator that has the highest potential to be violated due to the increased amount of power transfer [21]. The sensitive generator is important for the proposed method that assists towards fast computation of ATC. This is because the increase of power transfer that varies the maximum amount of generator rotor angle difference is observed only at the sensitive generator.
In order to observe the variations of maximum generator rotor angle differences, therefore, the power transfer is increased at 3 stages, which are 1 MW, 150.5 MW, and 300 MW for the transfer case from area 2 to area 3 and it is shown in Figure 6. It is observed that generator G5 gives rapid changes of maximum generator rotor angle difference due to the increase of power transfer from 1 MW to 300 MW. The maximum generator angle difference angle is 64.94°, 71.45°, and 78.5° for the increase of power transfer of 1 MW, 150.5 MW and 300 MW, respectively. This shows that G5 is the sensitive generator. On the other hand, the increase of the power transfer does not cause rapid changes of maximum generator rotor angle difference at generator G3. The maximum rotor angle difference of generator G3 is recorded at 74.14°, 73.93°, and 73.72° for 1 MW, 150.5 MW, and 300 MW of power transfer, respectively. Therefore, G3 is not a sensitive generator. Hence, a fast and accurate ATC computation could be obtained by considering the variations of maximum generator rotor angle difference only at the sensitive generator G5.
(a) Power transfer of 1 MW
(b) Power transfer of 150.5 MW
(c) Power transfer of 300 MW
The maximum generator rotor angle difference of G5 is then used in (3.12) to determine the lookahead power transfer, , of 850.67 MW. The is then used in the Ralston’s method to extrapolate the P curve at two trajectory points and it is shown in Figure 7. By referring to Figure 7, the system is actually operating in a stable condition when = 850.67 MW is less than the power transfer at the second trajectory point. Therefore, the that is used in the Ralston’s method is able to accurately extrapolate at the two trajectory points for optimum or accurate estimation of ATC. In Figure 7, the Ralston’s method performs four incremental steps of power transfer and it is reasonable enough to provide fast and accurate approximation of the two trajectory points. The two trajectory points are located at the power transfer of 850.67 MW and 1275.50 MW. It is believed that the actual value of ATC that causes the violation of generator rotor angle difference limit falls within these two points of trajectory and it can be observed in Figure 7.
Then, the AC power flow solution is performed at four incremental steps of power transfer, which are 850.67 MW, 992.28 MW, 1133.89 MW and 1275.50 MW. The four power transfers are obtained by equal division between the range of ATC = 850.67 MW and ATC = 1275.50 MW located at the two trajectory points. Figure 8 shows the variations of maximum generator rotor angle difference for the increase of four power transfers. It is obvious that the sensitive generator is still referring to G5. This is due to the fact that the increase of power transfers from 850.67 MW to 1275.50 MW causing the rapid changes of maximum generator rotor angle difference at G5. The maximum generator angle difference angle is 164.82°, 282.98°, 396.69°, and 529.41° for the increase of power transfer of 850.67 MW, 992.28 MW, 1133.89 MW, and 1275.50 MW, respectively.
(a) Power transfer of 850.67 MW
(b) Power transfer of 992.28 MW
(c) Power transfer of 1133.89 MW
(d) Power transfer of 1275.50 MW
Finally, the increase of four power transfers and the four maximum rotor angle differences are used in the cubicspline interpolation technique to fit the P curve between the two trajectory points of the sensitive generator G5. In Figure 9, it is shown that the cubicspline interpolation technique traces the P curve based on the four points of power transfer and maximum generator rotor angle difference. A point is noted where the generator rotor angle difference limit of 180° intersects the P curve. This point yields an actual ATC value of 870 MW.
4.8. Results of AreaToArea ATC due to Transient Stability Limit: Case Study of 39New England Bus System
Table 5 represents the results of ATC obtained only by considering the violation of generator rotor angle difference limit. This is to prove that the proposed method is also robust in determining the ATC based on the violation of transient stability limit. Whereby, the rotor angle difference limit is referred to as the transient stability limit. In this case study, it is assumed that a threephase fault happened at bus 10. Therefore, the computation of ATCsbased transient stability limit is performed by considering the tripping of line 10–13 for clearing the fault. The fault critical clearing time of second is specified for the tripping at line 10–13. The transient response of generator rotor angle difference is monitored for 1.5 seconds in the case study of ATC that takes into account the transient stability limit. The proposed method gives the ATC value of 448 MW for the interarea transfer from area 1 to area 2, and this is referring to the violation of the rotor angle difference limit of generator G3 which is located at bus 32. The proposed method accurately determines the ATC that is relatively similar to the ATC value of 449 MW determined by the recursive AC power flow method. This is similar to the other cases of interarea power transfer. By comparing the ATC values obtained in Tables 3 and 5, it can be concluded that the ATCs are actually obtained due to the violation of transmission line limit. It also shows that the violation of generator rotor angle difference limit gives higher value of ATC in comparison to the violation of transmission line limit. However, in a real system operating condition, the first violation may occur at the transmission line limit that yields the amount of power transfer. On the other hand, it is noted that the proposed method performs a fast ATC computation as compared to the recursive AC power flow method. The results have shown that the proposed method is able to provide fast and accurate value of ATC although it refers to the limit of transient stability. The results of pointtopoint ATC are not shown in this section. This is because any further increase of power transfer between buses does not violate the limit of generator rotor angle difference. This is due to the fact that any increase of power transfer from a particular generator does not adversely affect the synchronism of generating system.
5. Conclusion
This paper presents a new approach that is used to evaluate the areatoarea and pointtopoint ATCs. The proposed method is based on the Ralston method incorporating with cubicspline interpolation technique. The Ralston's method is used to determine the two trajectory points of voltage magnitude, power flow, and maximum generator rotor angle difference. Then, the cubicspline interpolation technique is used to trace the PV, PS, or P curves between the two trajectory points. The curve fitting procedure is performed to reduce the time in ATC computation. By referring to the PV, PS, or P curves between the two trajectory points, the ATC is determined at the intersecting point of voltage, MVA power flow or generator rotor angle difference limits. The effectiveness of the proposed method in determining the ATC is verified on a case study of 2737bus Polish and 39New England bus power systems. It is proven that the Ralston's method incorporating with cubicspline interpolation technique is a fast and accurate method for ATC evaluation as compared to the ATC method using recursive AC power flow method. The proposed method is an effective way to speed up the ATC computation. The proposed method is useful for the utilities in a deregulated electricity market in which the ATCs are required to be posted in a realtime market signal so that all transmission users have the same chance to access transmission information.
References
 G. Hamoud, “Assessment of available transfer capability of transmission systems,” IEEE Transactions on Power Systems, vol. 15, no. 1, pp. 27–32, 2000. View at: Publisher Site  Google Scholar
 North America Electric Reliability Council (NERC), “Available transfer capability definitions and determination,” June 1996, http://www.nerc.com/docs/docs/pubs/atcfinal.pdf. View at: Google Scholar
 P. W. Sauer, “Technical challenges of computing available transfer capability (ATC) in electric power systems,” in Proceedings of the 30th Annual Hawaii International Conference on System Sciences, vol. 5, pp. 589–593, Wailea, Hawaii, USA, 1997. View at: Google Scholar
 M. Shaaban, Y. Ni, H. Dai, and F. F. Wu, “Considerations in calculating total transfer capability,” in Proceedings of the International Conference Power System Technology (POWERCON '98), vol. 2, pp. 1356–1360, 1998. View at: Google Scholar
 M. Shaaban, Y. Ni, and F. F. Wu, “Transfer capability computations in deregulated power systems,” in Proceedings of the 33rd Annual Hawaii International Conference System Sciences, pp. 1–5, 2000. View at: Google Scholar
 R. D. Christie, B. F. Wollenberg, and I. Wangensteen, “Transmission management in the deregulated environment,” Proceedings of the IEEE, vol. 88, no. 2, pp. 170–195, 2000. View at: Publisher Site  Google Scholar
 M. M. Othman, A. Mohamed, and A. Hussain, “Fast evaluation of available transfer capability using cubicspline interpolation technique,” Electric Power Systems Research, vol. 73, no. 3, pp. 335–342, 2005. View at: Publisher Site  Google Scholar
 M. M. Othman, A. Mohamed, and A. Hussain, “A neural network based ATC assessment incorporating novel feature selection and extraction methods,” Electric Power Components and Systems, vol. 32, no. 11, pp. 1121–1136, 2004. View at: Publisher Site  Google Scholar
 G. C. Ejebe, G. D. Irisarri, S. Mokhtari, O. Obadina, P. Ristanovic, and J. Tong, “Methods for contingency screening and ranking for voltage stability analysis of power systems,” IEEE Transactions on Power Systems, vol. 11, no. 1, pp. 350–356, 1996. View at: Publisher Site  Google Scholar
 H.D. Chiang, C.S. Wang, and A. J. Flueck, “Lookahead voltage and load margin contingency selection functions for largescale power systems,” IEEE Transactions on Power Systems, vol. 12, no. 1, pp. 173–180, 1997. View at: Publisher Site  Google Scholar
 O. B. Fosso, N. Flatabo, B. Hakavik, and A. T. Holen, “Comparison of methods for calculation of margins to voltage instability,” in Proceedings of the Athens Power Technology Conference, Greece, 1993. View at: Publisher Site  Google Scholar
 S. C. Chapra and R. P. Canale, Numerical Methods for Engineers: With Software and Programming Application, McGrawHill, New York, NY, USA, 4th edition, 2003.
 R. D. Zimmerman and D. Gan, MATPOWER: version 3.2. School of Electrical Engineering, Cornell University, 2007.
 M. M. Othman, A. Mohamed, and A. Hussain, “Available transfer capability assessment using evolutionary programming based capacity benefit margin,” International Journal of Electrical Power and Energy Systems, vol. 28, no. 3, pp. 166–176, 2006. View at: Publisher Site  Google Scholar
 M. M. Othman, A. Mohamed, and A. Hussain, “Determination of transmission reliability margin using parametric bootstrap technique,” IEEE Transactions on Power Systems, vol. 23, no. 4, pp. 1689–1700, 2008. View at: Publisher Site  Google Scholar
 I. Dobson, S. Greene, and R. Rajaraman, “Electric power transfer capability: concepts, applications, sensitivity and uncertainty,” http://www.pserc.cornell.edu/tcc/info.html. View at: Google Scholar
 H. Saadat, Power System Analysis, McGrawHill, New York, NY, USA, 2nd edition, 2004.
 Y. Xia, K. W. Chan, M. Liu, and J. Wu, “Calculation of available transfer capability with transient stability constraints,” in Proceedings of the IEEE International Conference on Electric Utility Deregulation, Restructuring and Power Technologies (DRPT '04), vol. 1, pp. 128–132, Hong Kong, 2004. View at: Google Scholar
 Y. Yuan, J. Kubokawa, T. Nagata, and H. Sasaki, “A solution of dynamic available transfer capability by means of stability constrained optimal power flow,” in Proceedings of the IEEE Bologna Power Technology Conference, Italy, 2003. View at: Google Scholar
 M. H. Haque, “Novel method of finding the first swing stability margin of a power system from time domain simulation,” IEE Proceedings: Generation, Transmission and Distribution, vol. 143, no. 5, pp. 413–419, 1996. View at: Publisher Site  Google Scholar
 Y.K. Wu, “A novel algorithm for ATC calculations and applications in deregulated electricity markets,” International Journal of Electrical Power and Energy Systems, vol. 29, no. 10, pp. 810–821, 2007. View at: Publisher Site  Google Scholar
 M. ElHawary, Electrical Power Systems: Design and Analysis, Power Systems Engineering Series, WileyInterscience, IEEE Press, New York, NY, USA, 1995.
 Y. Ou and C. Singh, “Assessment of available transfer capability and margins,” IEEE Transactions on Power Systems, vol. 17, no. 2, pp. 463–468, 2002. View at: Publisher Site  Google Scholar
 M. A. Pai, Energy function analysis for power system stability, Kluwer Academic Publishers, Boston, 1989.
 C. K. Babulal and P. S. Kannan, “A novel approach for ATC computation in deregulated environment,” Journal of Electrical Systems, pp. 146–161, 2006. View at: Google Scholar
 P. Kundur, Power System Stability and Control, McGrawHill, New York, NY, USA, 1994.
 M. Pavella and P. G. Murthy, Transient Stability of Power Systems: Theory and Practice, John Wiley & Sons, New York, NY, USA, 1994.
 M. H. Haque and A. H. M. A. Rahim, “Determination of first swing stability limit of multimachine power systems through Taylor series expansions,” IEE Proceedings, vol. 136, pp. 373–379, 1989. View at: Google Scholar
Copyright
Copyright © 2010 Stendley Busan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.