A Lattice Method for Jump-Diffusion Process Applied to Transmission Expansion

Total Page:16

File Type:pdf, Size:1020Kb

A Lattice Method for Jump-Diffusion Process Applied to Transmission Expansion A Lattice Method for Jump-Diffusion Process Applied to Transmission Expansion Investments under Demand and Distributed Generation Uncertainties Fikri Kucuksayacigil a, K. Jo Min b a University of Arkansas, Fayetteville, AR, [email protected] b Industrial and Manufacturing Systems Engineering Department, Iowa State University, Ames, IA Corresponding author: Fikri Kucuksayacigil, 4179 Bell Engineering Center, University of Arkansas, Fayetteville, AR, [email protected], 515-203-9808 Abstract Strategic decision making for rare events is considered to be an arduous course of actions as their impacts cannot accurately be modeled. Transmission of electrical power is a major challenge in this day and age due to the prevalent of rare events. Decision makers of generation units, distribution utilities, and transmission networks do not often cooperate, which creates severe uncertainties for all players. Growth of demand for electricity and installation or removal of distributed generation (DG), considered to be a rare event, are among uncertainties encountered by transmission network owners. Expansion decisions for transmission lines should be strategically executed because installations of DGs may create a stranded cost for transmission owners. In this study, we propose a real options framework that quantifies the values of transmission investments under demand and DG uncertainties to guide decision makers of transmission companies regarding how to adapt their expansions decisions. We model demand uncertainty as a geometric Brownian motion (GBM) process and DG uncertainty as a Poisson arrival process. We devise a computationally-efficient lattice diagram to discretize both processes in a single grid, which can also be employed to model other types of rare events in various sectors. The proposed framework is demonstrated on a realistic transmission network. Numerical study shows that our proposed lattice model is computationally superior over existing diagrams. As for managerial insights, it shows that depending on the locations of DG installations, DG penetration may not reduce the value of a transmission network. If the center at which a DG is installed has a large-capacity local generator, the installation may not undervalue the transmission network. Keywords: Rare events, distributed generations, jump-diffusion process, lattice framework, real options, transmission expansion investments 1 Introduction After deregulation in the electricity market of the United States, decision makers of transmission companies / transmission owners confront critical uncertainties when they make investments. The reason is that transmission owners may not have prior information of decisions made by generation and distribution companies or communities. The trend in electricity demand is one of the severe uncertainties because it may exhibit large fluctuations [1]. Installation of distributed generations (DGs) creates another uncertainty because they have potential to defer large-scale transmission expansion investments. When transmission owners are not informed of DG installations, transmission investments tend to be redundant. Experts in the electricity market have initiated discussions toward evaluating the impact of DGs on costs and benefits of transmission investments. They point out that transmission investments could be better-planned if the rate of future adoption of DGs could be correctly estimated [2]. Interested readers can also see [3] for a discussion concerning transmission investment redundancy associated with a badly-designed distribution network. DGs have been installed to meet local demand for electricity during the previous decade in various sizes ranging 2 from only a couple of megawatts to tens of megawatts. Reference [4] shows a summary data listing various DG technologies preferred by utilities and societies as well as capacities installed in each year from 2006 to 2015. Transmission investments usually require huge capital costs, causing strategic decision making to be indispensable, and presence of uncertainties turns this problem to even more intractable form. A modeling framework to quantify the values of transmission investments in the presence of such uncertainties is a vital requirement. By values of transmission investment, we mean their monetary value, which is estimated with their future profits subject to an uncertain decision environment. This description is actually the definition of real options value of investments. In this study, we propose and demonstrate a real options approach to show how transmission investments can be assessed under continuous (infinitesimal changes in infinitesimal time interval) and discrete uncertainties (random shocks or rare events), namely uncertainties of growth in electricity demand and DG installations and removals. While a wide array of studies has been conducted for evaluating transmission investments under electricity demand uncertainty, the evaluation of DG penetration is a fairly new area of research. Most studies model DG installation uncertainty with discrete penetration scenarios, but our model stands unique because it uses a Poisson arrival process. It has a significant advantage over defining discrete scenarios because it is sufficient to determine only one parameter, arrival rate. As it will be seen in mathematical formulations, the only parameter as to uncertainty of DG penetration is arrival rate. As soon as we determine this parameter, we capture the corresponding uncertainty. Poisson distribution is an approximation of a counting process which would normally be modeled with a binomial distribution. First-order approximation enables us to suffice with arrival rate. In discrete penetration scenarios, one needs to know what scenarios are and what probabilities they have. It is challenging to estimate both 3 simultaneously with an acceptable accuracy. If we had perfect information about the distribution of uncertainties, we would be able to generate all relevant probability distributions. Investment evaluation problems modeled with real options methodology can be solved through three alternative approaches. While analytical techniques are worthwhile because managerial insights do not depend on numerical values of parameters, dealing with an analytically tractable model often requires advancing many unrealistic and restrictive assumptions. Monte Carlo simulation, proposed as one alternative for American options, provides researchers with modeling flexibilities such as ability of handling with jump and diffusion processes without enforcing a sequence between processes. On the other hand, it has a significant drawback from a computational perspective. For example, [5] ran through 50,000 paths to obtain an average value of American options for each of the stochastic processes underlying stock prices they were interested in modeling. Finally, with respect to lattice methods, it is often claimed that they require much less computation time compared to Monte Carlo simulation while returning more accurate results [6]. With lattice methods, however, one may not be sure about stability of results. In this study, we take advantage of lattice frameworks to model demand growth. As will be seen in Sect. 3, estimating the revenue of a transmission network requires solving an optimization problem. This problem is impractical to embed in an analytical model, which is the fundamental reason why we do not use an analytical approach We also want to avoid simulation due to its expensive computation time. There are numerous studies in the energy literature describing discretizing continuous evolutions with lattice or tree diagrams to solve problems without sacrificing their core properties [7 - 9]. Merit of our lattice model is that quantification of transmission network values is quickly accomplished. Although our framework finds the basis in a dynamic programming tree, often challenged because of its curse of dimensionality, we mitigate 4 this problem using approximation techniques, proposed in this study. As will be seen in next sections, our lattice model discretizes both diffusion and jump processes. Whereas diffusion process in our case is random smooth growth of demand for electricity, jump process is random installation of DGs in a power network. This study is structured as follows: In Sect. 2, we present the literature on work in close proximity to our paper. We then identify a discrete counterpart of the underlying uncertainty path in Sect. 3 and introduce a method to reduce its computational complexity. We also show how we quantify transmission network values. This is followed by a numerical example in Sect. 4 in which we demonstrate our framework on a realistic transmission network and show computational superiority of our proposed lattice model over existing ones. In Sect. 5, we discuss some assumptions made throughout the paper and present generalizations of our framework. Sect. 6concludes the paper by summarizing key points and important managerial insights. 2 Literature review A rich body of literature on analyzing transmission investments has been presented. Since we approach this problem by considering DG uncertainties, we restrict our attention to studies that incorporate DG penetrations with transmission investments. The following is a quick overview of those studies. Hejeejo and Qiu [10] examine transmission investments and power networks with intermittent DG resources. Investments are put in place to meet peak electricity demand while maintaining network
Recommended publications
  • Numerical Solution of Jump-Diffusion LIBOR Market Models
    Finance Stochast. 7, 1–27 (2003) c Springer-Verlag 2003 Numerical solution of jump-diffusion LIBOR market models Paul Glasserman1, Nicolas Merener2 1 403 Uris Hall, Graduate School of Business, Columbia University, New York, NY 10027, USA (e-mail: [email protected]) 2 Department of Applied Physics and Applied Mathematics, Columbia University, New York, NY 10027, USA (e-mail: [email protected]) Abstract. This paper develops, analyzes, and tests computational procedures for the numerical solution of LIBOR market models with jumps. We consider, in par- ticular, a class of models in which jumps are driven by marked point processes with intensities that depend on the LIBOR rates themselves. While this formulation of- fers some attractive modeling features, it presents a challenge for computational work. As a first step, we therefore show how to reformulate a term structure model driven by marked point processes with suitably bounded state-dependent intensities into one driven by a Poisson random measure. This facilitates the development of discretization schemes because the Poisson random measure can be simulated with- out discretization error. Jumps in LIBOR rates are then thinned from the Poisson random measure using state-dependent thinning probabilities. Because of discon- tinuities inherent to the thinning process, this procedure falls outside the scope of existing convergence results; we provide some theoretical support for our method through a result establishing first and second order convergence of schemes that accommodates thinning but imposes stronger conditions on other problem data. The bias and computational efficiency of various schemes are compared through numerical experiments. Key words: Interest rate models, Monte Carlo simulation, market models, marked point processes JEL Classification: G13, E43 Mathematics Subject Classification (1991): 60G55, 60J75, 65C05, 90A09 The authors thank Professor Steve Kou for helpful comments and discussions.
    [Show full text]
  • Superprocesses and Mckean-Vlasov Equations with Creation of Mass
    Sup erpro cesses and McKean-Vlasov equations with creation of mass L. Overb eck Department of Statistics, University of California, Berkeley, 367, Evans Hall Berkeley, CA 94720, y U.S.A. Abstract Weak solutions of McKean-Vlasov equations with creation of mass are given in terms of sup erpro cesses. The solutions can b e approxi- mated by a sequence of non-interacting sup erpro cesses or by the mean- eld of multityp e sup erpro cesses with mean- eld interaction. The lat- ter approximation is asso ciated with a propagation of chaos statement for weakly interacting multityp e sup erpro cesses. Running title: Sup erpro cesses and McKean-Vlasov equations . 1 Intro duction Sup erpro cesses are useful in solving nonlinear partial di erential equation of 1+ the typ e f = f , 2 0; 1], cf. [Dy]. Wenowchange the p oint of view and showhowtheyprovide sto chastic solutions of nonlinear partial di erential Supp orted byanFellowship of the Deutsche Forschungsgemeinschaft. y On leave from the Universitat Bonn, Institut fur Angewandte Mathematik, Wegelerstr. 6, 53115 Bonn, Germany. 1 equation of McKean-Vlasovtyp e, i.e. wewant to nd weak solutions of d d 2 X X @ @ @ + d x; + bx; : 1.1 = a x; t i t t t t t ij t @t @x @x @x i j i i=1 i;j =1 d Aweak solution = 2 C [0;T];MIR satis es s Z 2 t X X @ @ a f = f + f + d f + b f ds: s ij s t 0 i s s @x @x @x 0 i j i Equation 1.1 generalizes McKean-Vlasov equations of twodi erenttyp es.
    [Show full text]
  • A Switching Self-Exciting Jump Diffusion Process for Stock Prices Donatien Hainaut, Franck Moraux
    A switching self-exciting jump diffusion process for stock prices Donatien Hainaut, Franck Moraux To cite this version: Donatien Hainaut, Franck Moraux. A switching self-exciting jump diffusion process for stock prices. Annals of Finance, Springer Verlag, 2019, 15 (2), pp.267-306. 10.1007/s10436-018-0340-5. halshs- 01909772 HAL Id: halshs-01909772 https://halshs.archives-ouvertes.fr/halshs-01909772 Submitted on 31 Oct 2018 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. A switching self-exciting jump diusion process for stock prices Donatien Hainaut∗ ISBA, Université Catholique de Louvain, Belgium Franck Morauxy Univ. Rennes 1 and CREM May 25, 2018 Abstract This study proposes a new Markov switching process with clustering eects. In this approach, a hidden Markov chain with a nite number of states modulates the parameters of a self-excited jump process combined to a geometric Brownian motion. Each regime corresponds to a particular economic cycle determining the expected return, the diusion coecient and the long-run frequency of clustered jumps. We study rst the theoretical properties of this process and we propose a sequential Monte-Carlo method to lter the hidden state variables.
    [Show full text]
  • MERTON JUMP-DIFFUSION MODEL VERSUS the BLACK and SCHOLES APPROACH for the LOG-RETURNS and VOLATILITY SMILE FITTING Nicola Gugole
    International Journal of Pure and Applied Mathematics Volume 109 No. 3 2016, 719-736 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu AP doi: 10.12732/ijpam.v109i3.19 ijpam.eu MERTON JUMP-DIFFUSION MODEL VERSUS THE BLACK AND SCHOLES APPROACH FOR THE LOG-RETURNS AND VOLATILITY SMILE FITTING Nicola Gugole Department of Computer Science University of Verona Strada le Grazie, 15-37134, Verona, ITALY Abstract: In the present paper we perform a comparison between the standard Black and Scholes model and the Merton jump-diffusion one, from the point of view of the study of the leptokurtic feature of log-returns and also concerning the volatility smile fitting. Provided results are obtained by calibrating on market data and by mean of numerical simulations which clearly show how the jump-diffusion model outperforms the classical geometric Brownian motion approach. AMS Subject Classification: 60H15, 60H35, 91B60, 91G20, 91G60 Key Words: Black and Scholes model, Merton model, stochastic differential equations, log-returns, volatility smile 1. Introduction In the early 1970’s the world of option pricing experienced a great contribu- tion given by the work of Fischer Black and Myron Scholes. They developed a new mathematical model to treat certain financial quantities publishing re- lated results in the article The Pricing of Options and Corporate Liabilities, see [1]. The latter work became soon a reference point in the financial scenario. Received: August 3, 2016 c 2016 Academic Publications, Ltd. Revised: September 16, 2016 url: www.acadpubl.eu Published: September 30, 2016 720 N. Gugole Nowadays, many traders still use the Black and Scholes (BS) model to price as well as to hedge various types of contingent claims.
    [Show full text]
  • Sensitivity Estimation of Sabr Model Via Derivative of Random Variables
    Proceedings of the 2011 Winter Simulation Conference S. Jain, R. R. Creasey, J. Himmelspach, K. P. White, and M. Fu, eds. SENSITIVITY ESTIMATION OF SABR MODEL VIA DERIVATIVE OF RANDOM VARIABLES NAN CHEN YANCHU LIU The Chinese University of Hong Kong 709A William Mong Engineering Building Shatin, N. T., HONG KONG ABSTRACT We derive Monte Carlo simulation estimators to compute option price sensitivities under the SABR stochastic volatility model. As a companion to the exact simulation method developed by Cai, Chen and Song (2011), this paper uses the sensitivity of “vol of vol” as a showcase to demonstrate how to use the pathwise method to obtain unbiased estimators to the price sensitivities under SABR. By appropriately conditioning on the path generated by the volatility, the evolution of the forward price can be represented as noncentral chi-square random variables with stochastic parameters. Combined with the technique of derivative of random variables, we can obtain fast and accurate unbiased estimators for the sensitivities. 1 INTRODUCTION The ubiquitous existence of the volatility smile and skew poses a great challenge to the practice of risk management in fixed and foreign exchange trading desks. In foreign currency option markets, the implied volatility is often relatively lower for at-the-money options and becomes gradually higher as the strike price moves either into the money or out of the money. Traders refer to this stylized pattern as the volatility smile. In the equity and interest rate markets, a typical aspect of the implied volatility, which is also known as the volatility skew, is that it decreases as the strike price increases.
    [Show full text]
  • Option Pricing with Jump Diffusion Models
    UNIVERSITY OF PIRAEUS DEPARTMENT OF BANKING AND FINANCIAL MANAGEMENT M. Sc in FINANCIAL ANALYSIS FOR EXECUTIVES Option pricing with jump diffusion models MASTER DISSERTATION BY: SIDERI KALLIOPI: MXAN 1134 SUPERVISOR: SKIADOPOULOS GEORGE EXAMINATION COMMITTEE: CHRISTINA CHRISTOU PITTIS NIKITAS PIRAEUS, FEBRUARY 2013 2 TABLE OF CONTENTS ς ABSTRACT .............................................................................................................................. 3 SECTION 1 INTRODUCTION ................................................................................................ 4 ώ SECTION 2 LITERATURE REVIEW: .................................................................................... 9 ι SECTION 3 DATASET .......................................................................................................... 14 SECTION 4 .............................................................................................................................α 16 4.1 PRICING AND HEDGING IN INCOMPLETE MARKETS ....................................... 16 4.2 CHANGE OF MEASURE ............................................................................................ρ 17 SECTION 5 MERTON’S MODEL .........................................................................................ι 21 5.1 THE FORMULA ........................................................................................................... 21 5.2 PDE APPROACH BY HEDGING ................................................................ε ............... 24 5.3 EFFECT
    [Show full text]
  • Particle Filtering-Based Maximum Likelihood Estimation for Financial Parameter Estimation
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Spiral - Imperial College Digital Repository Particle Filtering-based Maximum Likelihood Estimation for Financial Parameter Estimation Jinzhe Yang Binghuan Lin Wayne Luk Terence Nahar Imperial College & SWIP Techila Technologies Ltd Imperial College SWIP London & Edinburgh, UK Tampere University of Technology London, UK London, UK [email protected] Tampere, Finland [email protected] [email protected] [email protected] Abstract—This paper presents a novel method for estimating result, which cannot meet runtime requirement in practice. In parameters of financial models with jump diffusions. It is a addition, high performance computing (HPC) technologies are Particle Filter based Maximum Likelihood Estimation process, still new to the finance industry. We provide a PF based MLE which uses particle streams to enable efficient evaluation of con- straints and weights. We also provide a CPU-FPGA collaborative solution to estimate parameters of jump diffusion models, and design for parameter estimation of Stochastic Volatility with we utilise HPC technology to speedup the estimation scheme Correlated and Contemporaneous Jumps model as a case study. so that it would be acceptable for practical use. The result is evaluated by comparing with a CPU and a cloud This paper presents the algorithm development, as well computing platform. We show 14 times speed up for the FPGA as the pipeline design solution in FPGA technology, of PF design compared with the CPU, and similar speedup but better convergence compared with an alternative parallelisation scheme based MLE for parameter estimation of Stochastic Volatility using Techila Middleware on a multi-CPU environment.
    [Show full text]
  • Lévy Finance *[0.5Cm] Models and Results
    Stochastic Calculus for L´evyProcesses L´evy-Process Driven Financial Market Models Jump-Diffusion Models General L´evyModels European Style Options Stochastic Calculus for L´evy Processes L´evy-Process Driven Financial Market Models Jump-Diffusion Models Merton-Model Kou-Model General L´evy Models Variance-Gamma model CGMY model GH models Variance-mean mixtures European Style Options Equivalent Martingale Measure Jump-Diffusion Models Variance-Gamma Model NIG Model Professor Dr. R¨udigerKiesel L´evyFinance Stochastic Calculus for L´evyProcesses L´evy-Process Driven Financial Market Models Jump-Diffusion Models General L´evyModels European Style Options Stochastic Integral for L´evyProcesses Let (Xt ) be a L´evy process with L´evy-Khintchine triplet (α, σ, ν(dx)). By the L´evy-It´odecomposition we know X = X (1) + X (2) + X (3), where the X (i) are independent L´evyprocesses. X (1) is a Brownian motion with drift, X (2) is a compound Poisson process with jump (3) distributed concentrated on R/(−1, 1) and X is a square-integrable martingale (which can be viewed as a limit of compensated compound Poisson processes with small jumps). We know how to define the stochastic integral with respect to any of these processes! Professor Dr. R¨udigerKiesel L´evyFinance Stochastic Calculus for L´evyProcesses L´evy-Process Driven Financial Market Models Jump-Diffusion Models General L´evyModels European Style Options Canonical Decomposition From the L´evy-It´odecomposition we deduce the canonical decomposition (useful for applying the general semi-martingale theory) Z t Z X (t) = αt + σW (t) + x µX − νX (ds, dx), 0 R where Z t Z X xµX (ds, dx) = ∆X (s) 0 R 0<s≤t and Z t Z Z t Z Z X X E xµ (ds, dx) = xν (ds, dx) = t xν(dx).
    [Show full text]
  • Estimating Single Factor Jump Diffusion Interest Rate Models
    ESTIMATING SINGLE FACTOR JUMP DIFFUSION INTEREST RATE MODELS Ghulam Sorwar1 University of Nottingham Nottingham University Business School Jubilee Campus Wallaton Road, Nottingham NG8 1BB Email: [email protected] Tele: +44(0)115 9515487 ABSTRACT Recent empirical studies have demonstrated that behaviour of interest rate processes can be better explained if standard diffusion processes are augmented with jumps in the interest rate process. In this paper we examine the performance of both linear and non- linear one factor CKLS model in the presence of jumps. We conclude that empirical features of interest rates not captured by standard diffusion processes are captured by models with jumps and that the linear CKLS model provides sufficient explanation of the data. Keywords: term structure, jumps, Bayesian, MCMC JEL: C11, C13, C15, C32 1 Corresponding author. Email [email protected] 1. INTRODUCTION Accurate valuation of fixed income derivative securities is dependent on the correct specification of the underlying interest rate process driving it. This underlying interest rate is modelled as stochastic process comprising of two components. The first component is the drift associated with the interest rates. This drift incorporates the mean reversion that is observed in interest rates. The second component consists of a gaussian term. To incorporate heteroskedasticity observed in interest rates, the gaussian term is multiplied with the short term interest rate raised to a particular power. Different values of this power leads to different interest rate models. Setting this power to zero yields the Vasicek model (1977), setting this power to half yields the Cox-Ingersoll-Ross model (1985b), setting this power to one and half yields the model proposed by Ahn and Gao (1999).
    [Show full text]
  • An Infinite Hidden Markov Model with Similarity-Biased Transitions
    An Infinite Hidden Markov Model With Similarity-Biased Transitions Colin Reimer Dawson 1 Chaofan Huang 1 Clayton T. Morrison 2 Abstract ing a “stick” off of the remaining weight according to a Beta (1; γ) distribution. The parameter γ > 0 is known We describe a generalization of the Hierarchical as the concentration parameter and governs how quickly Dirichlet Process Hidden Markov Model (HDP- the weights tend to decay, with large γ corresponding to HMM) which is able to encode prior informa- slow decay, and hence more weights needed before a given tion that state transitions are more likely be- cumulative weight is reached. This stick-breaking process tween “nearby” states. This is accomplished is denoted by GEM (Ewens, 1990; Sethuraman, 1994) for by defining a similarity function on the state Griffiths, Engen and McCloskey. We thus have a discrete space and scaling transition probabilities by pair- probability measure, G , with weights β at locations θ , wise similarities, thereby inducing correlations 0 j j j = 1; 2;::: , defined by among the transition distributions. We present an augmented data representation of the model i.i.d. θj ∼ H β ∼ GEM(γ): (1) as a Markov Jump Process in which: (1) some jump attempts fail, and (2) the probability of suc- G0 drawn in this way is a Dirichlet Process (DP) random cess is proportional to the similarity between the measure with concentration γ and base measure H. source and destination states. This augmenta- The actual transition distribution, πj, from state j, is drawn tion restores conditional conjugacy and admits a from another DP with concentration α and base measure simple Gibbs sampler.
    [Show full text]
  • Final Report (PDF)
    Foundation of Stochastic Analysis Krzysztof Burdzy (University of Washington) Zhen-Qing Chen (University of Washington) Takashi Kumagai (Kyoto University) September 18-23, 2011 1 Scientific agenda of the conference Over the years, the foundations of stochastic analysis included various specific topics, such as the general theory of Markov processes, the general theory of stochastic integration, the theory of martingales, Malli- avin calculus, the martingale-problem approach to Markov processes, and the Dirichlet form approach to Markov processes. To create some focus for the very broad topic of the conference, we chose a few areas of concentration, including • Dirichlet forms • Analysis on fractals and percolation clusters • Jump type processes • Stochastic partial differential equations and measure-valued processes Dirichlet form theory provides a powerful tool that connects the probabilistic potential theory and ana- lytic potential theory. Recently Dirichlet forms found its use in effective study of fine properties of Markov processes on spaces with minimal smoothness, such as reflecting Brownian motion on non-smooth domains, Brownian motion and jump type processes on Euclidean spaces and fractals, and Markov processes on trees and graphs. It has been shown that Dirichlet form theory is an important tool in study of various invariance principles, such as the invariance principle for reflected Brownian motion on domains with non necessarily smooth boundaries and the invariance principle for Metropolis algorithm. Dirichlet form theory can also be used to study a certain type of SPDEs. Fractals are used as an approximation of disordered media. The analysis on fractals is motivated by the desire to understand properties of natural phenomena such as polymers, and growth of molds and crystals.
    [Show full text]
  • Jump-Diffusion Modeling of a Markov Jump Linear System with Applications in Microgrids Gilles Mpembele, Member, IEEE, and Jonathan W
    1 Jump-Diffusion Modeling of a Markov Jump Linear System with Applications in Microgrids Gilles Mpembele, Member, IEEE, and Jonathan W. Kimball, Senior Member, IEEE Abstract—This paper proposes a jump-diffusion model for the the derivation of the ordinary differential equations that govern analysis of a microgrid operating in grid-tied and standalone the evolution of the system statistics. modes. The model framework combines a linearized power The application of stochastic models to the analysis of the system dynamic model with a continuous-time Markov chain (CTMC). The power system has a stochastic input modeled with dynamics of power systems is relatively new [4], [8]. This a one-dimensional Wiener process and multiplicative diffusion representation is particularly relevant for a class of stochastic coefficient. The CTMC gives rise to a compound Poisson pro- processes called Stochastic Hybrid Systems (SHSs). In this cess that represents jumps between different operating modes. framework, the linearized power system dynamic model is Starting from the resulting stochastic differential equation (SDE), combined with discrete transitions of the system variables the proposed approach uses Itoˆ calculus and Dynkin’s formula to derive a system of ordinary differential equations (ODEs) triggered by stochastic events [9], [10], [11], [4], [5], [7]. The that describe the evolution of the moments of the system. The linearized model is represented by a system of differential approach is validated with the IEEE 37-bus system modified to algebraic equations (DAEs) that includes the evolution of form a microgrid. The results match a Monte Carlo simulation the dynamic states, and of the output variables expressed as but with far lower computational complexity.
    [Show full text]