Parameterized fast decoupled power flow methods for obtaining the maximum loading point of power systems. Part I. Mathematical modeling

The conventional Newton and fast decoupled power flow (FDPF) methods have been considered inadequate to obtain the maximum loading point of power systems due to ill-conditioning problems at and near this critical point. It is well known that the PV and Q-theta decoupling assumptions of the fast decoupled power flow formulation no longer hold in the vicinity of the critical point. Moreover, the Jacobian matrix of the Newton method becomes singular at this point. However, the maximum loading point can be efficiently computed through parameterization techniques of continuation methods. In this paper it is shown that by using either theta or V as a parameter, the new fast decoupled power flow versions (XB and BX) become adequate for the computation of the maximum loading point only with a few small modifications. The possible use of reactive power injection in a selected PV bus (Q(PV)) as continuation parameter (mu) for the computation of the maximum loading point is also shown. A trivial secant predictor, the modified zero-order polynomial which uses the current solution and a fixed increment in the parameter (V, theta, or mu) as an estimate for the next solution, is used in predictor step. These new versions are compared to each other with the purpose of pointing out their features, as well as the influence of reactive power and transformer tap limits. The results obtained with the new approach for the IEEE test systems (14, 30, 57 and 118 buses) are presented and discussed in the companion paper. The results show that the characteristics of the conventional method are enhanced and the region of convergence around the singular solution is enlarged. In addition, it is shown that parameters can be switched during the tracing process in order to efficiently determine all the PV curve points with few iterations. (C) 2003 Elsevier B.V. All rights reserved.

Alternative parameters for the continuation power flow method

The conventional Newton's method has been considered inadequate to obtain the maximum loading point (MLP) of power systems. It is due to the Jacobian matrix singularity at this point. However, the MLP can be efficiently computed through parameterization techniques of continuation methods. This paper presents and tests new parameterization schemes, namely the total power losses (real and reactive), the power at the slack bus (real or reactive), the reactive power at generation buses, the reactive power at shunts (capacitor or reactor), the transmission lines power losses (real and reactive), and transmission lines power (real and reactive). Besides their clear physical meaning, which makes easier the development and application of continuation methods for power systems analysis, the main advantage of some of the proposed parameters is that its not necessary to change the parameter in the vicinity of the MLP. Studies on the new parameterization schemes performed on the IEEE 118 buses system show that the ill-conditioning problems at and near the MLP are eliminated. So, the characteristics of the conventional Newton's method are not only preserved but also improved. (C) 2003 Elsevier B.V. B.V. All rights reserved.

An Improved Parameterization Technique for the Continuation Power Flow

Continuation methods have been long used in P-V curve tracing due to their efficiency in the resolution of ill-conditioned cases, with close to singular Jacobian matrices, such as the maximum loading point of power systems. Several parameterization techniques have been proposed to avoid matrix singularity and successfully solve those cases. This paper presents a simple geometric parameterization technique to overcome the singularity of the Jacobian matrix by the addition of a line equations located at the plane determined by a bus voltage magnitude and the loading factor. This technique enlarges the set of voltage variables that can be used to whole P-V curve tracing, without ill-conditioning problems and no need of parameter changes. Simulation results, obtained for large realistic Brazilian and American power systems, show that the robustness and efficiency of the conventional power flow are not only preserved but also improved.

Continuation load flow using automatically determined branch megawatt losses as parameters

The conventional Newton's method is considered to be inadequate for the computation of the maximum loading point (MLP) of power systems since: (i) it encounters difficulties in the vicinity of the MLP: and (ii) the load flow Jacobian matrix becomes singular at the MLP. It is well known that continuation methods are powerful and useful tools that are able to trace the solution PV curve without experiencing such diffculties. However, continuation methods require a parameterisation so that a modified, well conditioned set of load flow equations is obtained. In particular, the Jacobian matrix associated with this modified set of equations should not be singular at the MLP. The authors propose that the actual power losses in transmission branches (lines and transformers) are used to parameterise the approach. Specific procedures for the automatic determination of the most appropriate parameter (branch) are proposed. Such procedures include the utilisation of fast voltage-stability indices. Simulation results are presented to show that the proposed method is able to trace the whole solution PV curve very efficiently.

Study of alternative schemes for the parameterization step of the continuation power flow method based on physical parameters, part I: Mathematical modeling

The conventional power flow method is considered to be inadequate to obtain the maximum loading point because of the singularity of Jacobian matrix. Continuation methods are efficient tools for solving this kind of problem since different parameterization schemes can be used to avoid such ill-conditioning problems. This paper presents the details of new schemes for the parameterization step of the continuation power flow method. The new parameterization options are based on physical parameters, namely, the total power losses (real and reactive), the power at the slack bus (real or reactive), the reactive power at generation buses, and transmission line power losses (real and reactive). The simulation results obtained with the new approach for the IEEE test systems (14, 30, 57, and 118 buses) are presented and discussed in the companion paper. The results show that the characteristics of the conventional method are not only preserved but also improved.

An efficient geometric parameterization technique for the continuation power flow

Continuation methods have been shown as efficient tools for solving ill-conditioned cases, with close to singular Jacobian matrices, such as the maximum loading point of power systems. Some parameterization techniques have been proposed to avoid matrix singularity and successfully solve those cases. This paper presents a new geometric parameterization scheme that allows the complete tracing of the P-V curves without ill-conditioning problems. The proposed technique associates robustness to simplicity and, it is of easy understanding. The Jacobian matrix singularity is avoided by the addition of a line equation, which passes through a point in the plane determined by the total real power losses and loading factor. These two parameters have clear physical meaning. The application of this new technique to the IEEE systems (14, 30, 57, 118 and 300 buses) shows that the best characteristics of the conventional Newton's method are not only preserved but also improved. (C) 2006 Elsevier B.V. All rights reserved.

Assessment of branch outage contingencies using the continuation method

This paper provides a contribution to the contingency analysis of electric power systems under steady state conditions. An alternative methodology is presented for static contingency analyses that only use continuation methods and thus provides an accurate determination of the loading margin. Rather than starting from the base case operating point, the proposed continuation power flow obtains the post-contingency loading margins starting from the maximum loading and using a bus voltage magnitude as a parameter. The branch selected for the contingency evaluation is parameterised using a scaling factor, which allows its gradual removal and assures the continuation power flow convergence for the cases where the method would diverge for the complete transmission line or transformer removal. The applicability and effectiveness of the proposed methodology have been investigated on IEEE test systems (14, 57 and 118 buses) and compared with the continuation power flow, which obtains the post-contingency loading margin starting from the base case solution. In general, for most of the analysed contingencies, few iterations are necessary to determine the post-contingency maximum loading point. Thus, a significant reduction in the global number of iterations is achieved. Therefore, the proposed methodology can be used as an alternative technique to verify and even to obtain the list of critical contingencies supplied by the electric power systems security analysis function. © 2013 Elsevier Ltd. All rights reserved.

Topics in bifurcation theory

I. Existence and Structure of Bifurcation Branches

The problem of bifurcation is formulated as an operator equation in a Banach space, depending on relevant control parameters, say of the form G(u,λ) = 0. If dimN(G_u(u_O,λ_O)) = m the method of Lyapunov-Schmidt reduces the problem to the solution of m algebraic equations. The possible structure of these equations and the various types of solution behaviour are discussed. The equations are normally derived under the assumption that G^O_λεR(G^O_u). It is shown, however, that if G^O_λεR(G^O_u) then bifurcation still may occur and the local structure of such branches is determined. A new and compact proof of the existence of multiple bifurcation is derived. The linearized stability near simple bifurcation and "normal" limit points is then indicated.

II. Constructive Techniques for the Generation of Solution Branches

A method is described in which the dependence of the solution arc on a naturally occurring parameter is replaced by the dependence on a form of pseudo-arclength. This results in continuation procedures through regular and "normal" limit points. In the neighborhood of bifurcation points, however, the associated linear operator is nearly singular causing difficulty in the convergence of continuation methods. A study of the approach to singularity of this operator yields convergence proofs for an iterative method for determining the solution arc in the neighborhood of a simple bifurcation point. As a result of these considerations, a new constructive proof of bifurcation is determined.

Obtaining bifurcation diagrams with a thermoacoustic network model

Linear techniques can predict whether the non-oscillating (steady) state of a thermoacoustic system is stable or unstable. With a sufficiently large impulse, however, a thermoacoustic system can reach a stable oscillating state even when the steady state is also stable. A nonlinear analysis is required to predict the existence of this oscillating state. Continuation methods are often used for this but they are computationally expensive. In this paper, an acoustic network code called LOTAN is used to obtain the steady and the oscillating solutions for a horizontal Rijke tube. The heat release is modelled as a nonlinear function of the mass flow rate. Several test cases from the literature are analysed in order to investigate the effect of various nonlinear terms in the flame model. The results agree well with the literature, showing that LOTAN can be used to map the steady and oscillating solutions as a function of the control parameters. Furthermore, the nature of the bifurcation between steady and oscillating states can be predicted directly from the nonlinear terms inside the flame model. Copyright © 2012 by ASME.

Three Cuts for Accelerated Interval Propagation

This paper addresses the problem of nonlinear multivariate root finding. In an earlier paper we described a system called Newton which finds roots of systems of nonlinear equations using refinements of interval methods. The refinements are inspired by AI constraint propagation techniques. Newton is competative with continuation methods on most benchmarks and can handle a variety of cases that are infeasible for continuation methods. This paper presents three "cuts" which we believe capture the essential theoretical ideas behind the success of Newton. This paper describes the cuts in a concise and abstract manner which, we believe, makes the theoretical content of our work more apparent. Any implementation will need to adopt some heuristic control mechanism. Heuristic control of the cuts is only briefly discussed here.

Análisis de sensibilidad para un modelo de medición de percepción de marca basado en simulación computacional por agentes

Se presenta el análisis de sensibilidad de un modelo de percepción de marca y ajuste de la inversión en marketing desarrollado en el Laboratorio de Simulación de la Universidad del Rosario. Este trabajo de grado consta de una introducción al tema de análisis de sensibilidad y su complementario el análisis de incertidumbre. Se pasa a mostrar ambos análisis usando un ejemplo simple de aplicación del modelo mediante la aplicación exhaustiva y rigurosa de los pasos descritos en la primera parte. Luego se hace una discusión de la problemática de medición de magnitudes que prueba ser el factor más complejo de la aplicación del modelo en el contexto práctico y finalmente se dan conclusiones sobre los resultados de los análisis.

Método da continuação aplicado na análise de contingência de linhas de transmissão

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Esquemas alternativos para o passo de parametrização do método da continuação baseados em parâmetros físicos

O método de fluxo de carga convencional é considerado inadequado para se obter o ponto de máximo carregamento (PMC) de sistemas de potência, devido à singularidade da matriz Jacobiana neste ponto. Os métodos da continuação são ferramentas eficientes para a solução deste tipo de problema, visto que técnicas de parametrização podem ser utilizadas para evitar a singularidade da matriz Jacobiana. Neste trabalho, novas opções para a etapa de parametrização do método da continuação são apresentadas. Mostra-se que variáveis com claro significado físico podem ser utilizadas na etapa de parametrização. As seguintes variáveis foram testadas: perda total de potência ativa e reativa, potência ativa e reativa na barra de referência, potência reativa das barras de geração, e as perdas de potência ativa e reativa nas linhas de transmissão (LT). Além de facilitar a implementação computacional do método de continuação, as técnicas de parametrização apresentadas simplificam a definição matemática e o entendimento do método por parte de engenheiros de potência, visto que os métodos de continuação existentes na literatura sempre utilizam técnicas de parametrização complexas, e de interpretação puramente geométrica. Resultados obtidos com a nova metodologia para os sistemas testes do IEEE (14, 30, 57 e 118 barras) mostram que as características de convergência do método de fluxo de carga convencional são melhoradas na região do PMC. Além disso, durante o traçado das curvas PV, as diversas técnicas de parametrização podem ser comutadas entre si possibilitando o cálculo de todos os pontos da curva com um número reduzido de iterações. Diversos testes são realizados para proporcionar a comparação do desempenho dos esquemas de parametrização propostos.

Determinação do ponto de máximo carregamento de sistemas de potência utilizando o fluxo de carga desacoplado rápido parametrizado

Os métodos de fluxo de carga por Newton-Raphson e fluxo de carga desacoplado rápido convencionais são considerados inadequados para a obtenção do ponto de máximo carregamento de sistemas de potência, devido à problemas de mal-condicionamento neste ponto crítico e na sua vizinhança. Neste ponto a matriz Jacobiana do método de Newton-Raphson torna-se singular e considera-se que não são mais válidas as hipóteses de desacoplamento P-V e Q-teta utilizadas para a formulação do método fluxo de carga desacoplado rápido. No entanto, mostra-se neste trabalho, que com pequenas modificações, as versões XB e BX do fluxo de carga desacoplado rápido tornam-se adequadas para a obtenção do ponto de máximo carregamento. Estas novas versões modificadas são comparadas entre si com o intuito de explicitar suas características, assim como da influência da atuação dos limites de geração de potência reativa e de tap's de transformadores. Os resultados obtidos para os sistemas testes do IEEE (14, 30, 57 e 118 barras) mostram que as características de convergência das versões originais são preservadas. Além disso, durante o traçado das curvas PV, os diversos métodos podem ser comutados entre si possibilitando o cálculo de todos os pontos da curva com um número reduzido de iterações.

Fluxo de carga continuado utilizando estratégia de parametrização geométrica baseada em parâmetros físicos

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)