Interior-Point Methods Applied to the Predispatch Problem of a Hydroelectric System with Scheduled Line Manipulations

Transmission line manipulations in a power system are necessary for the execution of preventative or corrective maintenance in a network, thus ensuring the stability of the system. In this study, primal-dual interior-point methods are used to minimize costs and losses in the generation and transmission of the predispatch active power flow in a hydroelectric system with previously scheduled line manipulations for preventative maintenance, over a period of twenty-four hours. The matrix structure of this problem and the modification that it imposes on the system is also broached in this study. From the computational standpoint, the effort required to solve a problem with or without line manipulations is similar, and the reasons for this are also discussed in this study. Computational results sustain our findings.


Introduction
In a power system, continual improvements in operational conditions are always sought; however, to achieve such improvements, constant maintenance and improvement work in the network is required, and scheduled shut-downs are sometimes required.Maintenance is necessary to avoid electrical short circuits and overloads.Thus, the objective is to ensure the continual supply of electric power, with the fewest interruptions of the shortest duration possible, thus maintaining the service levels required by the legislation.
For instance, in Brazil the System Operation Center (COS) is responsible for executing, authorizing and supervising scheduled and emergency line manipulations in the power transmission system, as well as for monitoring the system and reestablishing the power grid in the event of isolated or generalized contingencies.These activities, executed in real time, require situational awareness and management of the execution of the necessary line manipulations, with the aim of ensuring the integrity of personnel and installations, whilst guaranteeing the reliability of the system and the continuity and quality of supply.
In this study, the predispatch will be modeled on a twenty-four hour period, representing the dispatch of a single day, and line manipulations will be an input data point; in other words, they will already have been scheduled, and thus the line manipulations considered herein will be of a preventative nature.It is important to emphasize that there are other types of manipulations, such as generator shut downs, tap changing procedures and flexible alternating current transmission system (FACTS) adjustments, but these are not considered in this study.
Furthermore, by making use of the speed and robustness of interior-point methods [1][2][3], the intention is to obtain more efficient implementations for predispatch by means of exploiting the matrix structure of the resulting system.

The Predispatch Problem
Predispatch is a short-term operational problem, in this case, short-term refers to one-week, or even one-day, planning.The intention is to fulfill demands and satisfy the energy targets that have been defined in long-term planning.
Flow restrictions can be divided into blocks that repeat themselves over a certain time interval, representing the power system in these intervals, and thus we obtain an independent formulation of the Kirchhoff laws, where power flows are represented in order to enable a direct consideration of the transmission limits as restrictions and transmission losses as a performance criterion [4].
A predispatch system with m buses, n lines and g generators, where line manipulations are considered, can be modeled in the following manner: , where,  represents the active power flow;  represents the active power generation; represents the quadratic component of the generation cost;    and  objective function weight constants; represents the power generation target established in long-term planning.For this model, the two objective function components are quadratic with separable variables.The first component represents the value of the transmission losses.The second component characterizes the generation cost of the plants [5].
Problem (1) may be simplified using changes in variables [6] and adding slack variables, thus we get the primal problem in its standard form: , Whose respective dual problem may be expressed as: Matrix B, formed by the juxtaposed rows of the incidence and reactance matrices, is no longer constant throughout the time intervals t, and may be partioned as: In a more detailed form: With this partitioning, the columns of incidence matrix A are divided so that T contains the edges of a spanning tree and N is formed by the remaining edges, which belong to the co-tree [7], and the reactance matrix X is partitioned in a similar manner.
Matrices B and E vary in accordance with the time intervals (B k and E k ), reflecting the modifications to the networks and buses by line manipulations carried out throughout the study window.The reason that these matrices vary according to the time intervals is that the network is no longer constant throughout these t-intervals.Every time a line manipulation occurs, matrix B formed by the juxtaposed rows of the incidence and reactance matrix, and matrix E of order m × g should vary in accordance with the changes imposed on the network.
The primal-dual interior-point methods consist of the application of Newton's method to the optimality conditions [8]; therefore, the following linear system is obtained: when several variable substitutions are made, we get: whose direct solution requires great computational effort, because, has the dimension of the number of generators.

Scheduled Line Manipulations
Line manipulations are executed with the intention of adapting the transmission network to the load variation anipulation azilian power system our to six line manipulations are executed per day.The changes considered herein are other words, Matrix B is formed by the juxtaposed rows of the inonstant thro e that a line the number of rows and dy a has the dimension of (variation of the power demand) throughout the day.In most time intervals, line m s are not executed in the Br , which means that the transmission network rarely alters from one time interval to another, normally f called preventative line manipulations; in alterations to the network so that maintenance can be carried out and power blackouts avoided.

Study of the Matrix Structure for the Problem with Line Manipulations
We assume that in this system, i previously scheduled line manipulations take place over t time intervals, where each time interval corresponds to the period of one hour, and i line manipulations are few with the average value ranging from zero to six, which is typical of the Brazilian system.
cidence and reactance matrices, and is no longer c ughout t time intervals, because every tim manipulation is executed, we can infer that a row and a column of matrix B are removed (inserted), in the event that more than one line manipulation occurs in the same time interval, more rows and columns of matrix B are removed (inserted).Thus, when we consider a system with line manipulations at different time intervals, we will use the following representation: where, 1, 2, , 24 k   Every time a line manipulation occurs during interval k, matrix B will have its size modified, because it is formed by the incidence matrix, which represents the network topology, and by the reactance matrix.Therefore, if a line is disconnected from the grid at time k, then the incidence matrix should represent the network at that exact instant in time; in other words, a column of A is ed, therefore the reactance associated with that row should also be extracted from the matrix.eliminat As the size of matrix B may be modified with each line manipulation, we must adjust the system to these changes; in other words, in order to obtain the product and the sum total with B, the size of some of the matrices involved in the system should also be altered.
The initial objective is to solve the linear system ( 4), but for this we realize that too much computational effort would be necessary to solve it directly, because in the first equation, we have the following matrix: whose size is the number of rows at instant t, while the size of vector dy a is the number of generators.A more efficient solution follows these steps [6]: Step 1) Consider a matrix k B  , consisting of matrix B k and canonical vector e j , thus (note that this matrix is square and nonsingu r).
Step 2) A row and a column is added to matrix la   , in order to adjust its size for multiplication of the matrices.
Step 3) In matrix D k , remove the j th th row and the j co ingly, we are removing a generation bus from the matrix and its size is not altered, rather only replaced with zero in the oved j th row and lumn, with j varying from 1, 2, •••, m, where its intersections cannot be zero.Accord rem column.
which has the size (n + 1) × (n + 1), and the vector Step 4) We define: this system will be solved in two stages [9 We will first solve the following linear system: ]: Accordingly, the "inversions" of and m is important to emphasize that we are supposing that in t time intervals, i-scheduled line In order to solve (5), we have dy which is difficulties, using, for example, the LU factor of difficulties, using, for example, the LU factor of found without h does vary throughout the iterations.In algerbraic form found without h does vary throughout the iterations.In algerbraic form whic whic not , we write: not , we write: The Sherman-Morrison-Woodbury formula [10] is used to solve the linear system (5), in this case the computation of invertible matrix where U and V are matrices of size p × q and S is of size q × q, and therefore suitable for our problem  S is a diagonal matrix of size g × g, whose elements are those that belong to matrix k Matrix its columns removed from the identity matrix can be decomposed once only before the iterative process, in an identical fashion to what can be done with matrix B, in the case of the problem without line manipulations.

Implementation Details for the Developed Methods
In order to devinterior-points methods for a system with line manipulations, the network s to be adap d.With each line manipulation carried out, a row an a co ved) from m d), r It is worth em sizing that as the spanning tree is not t.
te d need lumn can be inserted (remo atrix B. The spanning tree is represented by matrix T and must not be modified; in other words, only the branches belonging to matrix N, formed by the additional branches of the spanning tree, can be connected (disconnecte thus facilitating implementation and resulting in greate computational efficiency.pha unique and the line manipulations are previously known, any branch may be manipulated, and all this would require is the construction of another spanning tree, as long as the branch to be manipulated did not belong to i The tests undertaken in this study use the starting point shown in the Equation ( 6).This was defined as in [12], which presented satisfactory results in previous experiments.

Computational Results
The networks in which the tests were carried out include the IEEE30 bus system, representing the Midwest of the United States of America, and also South-Southeastern-Central-West region wi (SSECO1732) and the Brazilian interconnected system, consisting of 1993 buses.
In the computation experiments carried out, the primal-dual interior-point method was used and tests were n with the number of line manipulations varying from which is what normally stem.
ipulated in the grids were chose ere to di the Brazilian th 1732 buses ru zero to six over a 24-hour period, happens in practice in the Brazilian power sy The implementation was developed using Matlab 7.0, with a precision of 10 -3 , in order to satisfy the optimal conditions of the problem.The computer has an Intel Centrino Core 2 Duo processor, with 4 GB of RAM and a speed of 2.13 GHz.
The branches to be man n randomly, and we emphasize that several tests were carried out and drastic differences from one test to another were noted, because the number of iterations and the computational time can vary significantly, depending on the branch that is being manipulated.If we w sconnect a branch where the flow was high, the grid would have to find a new way of fulfilling the demand, thus reflecting on the number of iterations.However, other reasons may be responsible for the increase in the computational effort, because depending on the branches m branch with a high flo anipulated, the network resulting from these manipulations may not adapt very efficiently.
The first column of the Tables 1, 2 and 3 pertains to the number of manipulations and they are accumulative; for example, in the event of three manipulations, it consists of the branches previously manipulated (1 and 2) plus the new third branch.
Note that in all the tests carried out, the line manipulations represent a greater computational effort, but one that is acceptable due to the amount of support they provide to the grid.The number of iterations necessary for convergence may vary significantly, depending on the branch that is being manipulated.If a w is disconnected, the grid will have to find a new

Conclusions
In this study, interior-point methods are used to solve a predispatch problem in a hydroelectric system.This research contributes to solving this problem with transmission line manipulations.When these manipulations take place, the topology of the network is altered.The characteristics of this problem and its importance to the Brazilian power system have been the drivers behind this research.
In practice, the number of mod work is small, vary period.In this study network are known in advance.Thus, the matrices associated with the problem can be analyzed and decomposed before the interior-point methods are applied.An identical approach has been taken by other authors for a scenario without manipulations.
Line manipulations represent the disconnection and/or return to operational status, of certain transmission lines.For each occ e incidence matrix and a row in the reactance matrix are removed.The methodology used in the development of this study is the primal-dual interior-point method, because this presents satisfactory results for optimum power flow problems.
From the computational standpoint, the effort to solve a problem with or without modification in network topology is similar.Even with oblem matrix, the number of linear systems that need to be solved is still the same, compared to the problem without manipulations.Furthermore, the number of iterations necessary for the convergence of the interior-point method depends on how important the manipulated branches are to the grid.
The new method does not result in greater costs, although it is known that anches that are being disconnected, because if the load of a manipulated branch is very high, the grid may not be able to readapt, thus resulting in convergence challenges isting approaches that consider schedu nowledgements .Innorta and R. Ricci, "The Flexibility of wer Flow Algorithms Facing s," Electrical Power E Energy in this study may be proven not only by means of convergence, but also for non-convergence of some of the tests carried out.Accordingly, it is possible to prove that the implementation actually optimizes the problems covered.
There are other ex led transmission lines manipulation.For instance, line manipulation can be simulated.In this case, there is not network topology modification as oppose of the adopted approach.However, it presents stability difficulties and higher computational cost.
Thus, the solved model gets much closer to the actual problem with an insignificant additional computational effort, and without presenting any issues regarding numerical stability.


E represents a matrix of order m × g, with each column containing exactly one element equal to 1 and all other elements equal to 0; of this number of line manipulations