# Valuation of Swing Options under a Regime-Switching Mean-Reverting Model

Authors: Lingjie Shao; Kaili Xiang
Journal: Mathematical Problems in Engineering
 (2019)

---

## Abstract

In this paper, we study the valuation of swing options on electricity in a model where the underlying spot price is set to be the product of a deterministic seasonal pattern and Ornstein-Uhlenbeck process with Markov-modulated parameters. Under this setting, the difficulties of pricing swing options come from the various constraints embedded in contracts, e.g., the total number of rights constraint, the refraction time constraint, the local volume constraint, and the global volume constraint. Here we propose a framework for the valuation of the swing option on the condition that all the above constraints are nontrivial. To be specific, we formulate the pricing problem as an optimal stochastic control problem, which can be solved by the trinomial forest dynamic programming approach. Besides, empirical analysis is carried out on the model. We collect historical data in Nord Pool electricity market, extract the seasonal pattern, calibrate the Ornstein-Uhlenbeck process parameters in each regime, and also get market price of risk. Finally, on the basis of calibration results, a specific numerical example concerning all typical constraints is presented to demonstrate the valuation procedure.


---

## Body

## 1. Introduction

In the long run, energy liberalisation in electricity markets improves economic benefits, since that it overcomes the overcapacity problem in regulated markets and improves efficiency in the operation of networks and transport services (see [1] for more detail). However, in the short term, it exposes participants to various risks, e.g., price risks and quantity risks, because that electricity energy is hard to store and agents usually can not control their consumption. The emergence of swing options brings participants opportunities to hedge these risks according to [2]. Swings provide the option holder with flexibilities in deciding when to exercise and choosing the trading volume of each exercise, with strike price predetermined. However, these flexibilities usually are limited for the benefit of the option writer. Firstly, the chosen exercise times should be in the range of a series of predetermined days. Secondly, the total number of exercise rights is limited. Thirdly, each time the holder chooses to exercise, they face the upper and lower bounds for the delivery volume. Fourthly, the total amount of delivery volume is limited. Fifthly, the time span between the two consecutive exercises must be longer than the refraction time. All the details of these constraints are presented in [3]. Comparing to vanilla options, all of these above features add complexity to the valuation of swings.For the valuation of option, we need to model the underling spot price. According to [4], electricity price in real markets usually has three main characteristics, including seasonal pattern, mean-reversion and spikes. We refer the reader to [5] for various models used in electricity market. Among them, jump models and regime-switching models are popular. In particular, if the real data shows that the price remained high for some time after a jump, then regime-switching models are better to depict this phenomenon [6]. For example, different types of regime-switching models were selected in [6–10]. Besides, various kinds of jump-diffusion models were applied in [11–13].Adopting a specific setting on swings, Carmona and Touzi [14] applied the Snell envelope theory to formulate the pricing of swings as an optimal multiple-stopping problem, which can be reduced to compound single stopping problem. And they analyzed this problem in the Black-Scholes framework and use Monte-Carlo simulation to construct a numerical example. To tackle the optimal multiple-stopping problem in the form of several variational inequalities, finite element method was used by M. Wilhelm and C. Winter [15], and dual method was used in [16–18]. Moreover, optimal stopping boundaries associated with swings’ pricing problems were studied in [19–21].Besides, the valuation problem can also be formulated as an optimal stochastic control problem subject to constraints. Facing this problem, many researchers employed the dynamic programming method, which lead to the Bellman’s equation for the value function. And to deal with the conditional expectations in the Bellman’s equation, different kinds of numerical methods were applied; e.g., tree methods were used in [3, 22, 23], Monte-Carlo method was used in [24, 25], and the finite difference method was used in [26, 27]. Besides the dynamic programming method, the method of discretizing the underlying probability space was used in [28] to overcome the difficulty of facing infinitely many cases.In this paper, to characterize the electricity price, we use a seasonal regime-switching mean-reverting model withn states, after adding the mean-reverting characteristic into the model adopted in [22]. The motivation to add this modification comes from empirical researches of [6], in which electricity spot prices have been verified to have three basic features, including seasonal pattern, mean-reversion, and spikes. As for the constraints of option, besides usual ones, we take refraction time constraint and a penalty function of more general form into consideration based on [22]. Under these settings, we formulate the pricing of swings as an optimal stochastic control problem subject to several constraints, including the refraction time, the local volume, and the global volume. Then a trinomial forest dynamic programming approach is adopted, and more applications of this approach can be found in [3, 22, 23], etc. Finally, based on the daily average data in Nord Pool, we calibrate the model with two regimes and construct a concrete example concerning all nontrivial constraints and a penalty function of general form.The structure of this paper is as follows. Section2 introduces the model for the underlying electricity price and characterizes the swing option. Section 3 presents the procedure of valuing swing options with the trinomial forest dynamic programming approach. Section 4 concentrates on the empirical calibration and finding the market price of risk when the number of regime is two. Section 5 gives a numerical example based on the results in Section 4 for the valuation of the swing option with all nontrivial constrains. Section 6 concludes the paper.


## 2. Electricity Price Model and Characteristics of Swing Options

### 2.1. Characteristics of Electricity Spot Prices

Our modelling is based on the electricity spot price data of the day-ahead market in Nord Pool, and we focus on the stochastic component of the deseasonalized prices. Firstly, the stochastic component has the mean-reverting feature captured by Lucia and Schwartz [29], which shows that the electricity prices tend to revert to the mean value over time. Secondly, the time series of electricity prices show spikes, which can sometimes be depicted by jump models or regime-switching models (see [30]). All these two characteristics of the deseasonalized prices are the basis of modelling in this paper.


### 2.2. Models for Electricity Spot Prices

Throughout this paper, for a fixed maturityT>0, we work on a probability space (Ω,F,{Ft}0≤t≤T,P), where {Ft}0≤t≤T is the natural filtration generated by Brownian Motion {Bt}0≤t≤T in physical measure P and an independent observable continuous-time Markov chain {αt}0≤t≤T. Assume that {αt}0≤t≤T takes values in a finite state space M:={1,2,⋯,m0} and is generated by Q={qij}, where qij≥0 for i,j∈M,j≠i, and ∑j=1m0‍qij=0 for each i∈M. The transitional probability is given by(1)PijΔ≔Pαt+Δ=j∣αt=i=qijΔ+oΔ,j≠i,1+qijΔ+oΔ,j=i,for a small time interval Δ>0. Besides, we denote by ρi the probability of being in state i at the initial time t=0, i.e., (2)ρi≔Pα0=i,i∈M.We start our modelling for the spot price from the one-factor model used in [29], which captures seasonal and mean-reversion feature. The model is(3)St=exp⁡ft+Xt,(4)dXt=-βXtdt+σdBtP.In (3), the logarithm of the electricity price St is divided into two parts, one of which is the seasonal part f(t), and the other is the deseasonalized part Xt. Equation (4) describes the dynamic of the deseasonalized part by Vasicek model, where the long-term mean is set to be zero, and the speed of reversion together with instantaneous volatility is denoted by β and σ respectively, both of which are positive constants.In this paper, our model is(5)St=ft·Xt,(6)dXt=βαtξαt-Xtdt+σαtdBtP.In (5), the modification is based on the fact that the electricity price can be negative (see [5]). And in (6), to capture the feature that electricity price show spikes, we replace (4) with a n-state regime-switching mean-reverting model, which was proposed in [8, 10]. Specifically, β(·):M→R+, ξ(·):M→R, and σ(·):M→R+ represent the speed of adjustment, the long-term mean level and the volatility, respectively. Note that, for simplicity, we denote βαt≔β(αt),ξαt≔ξ(αt),σαt≔σ(αt).


### 2.3. Swings and Pricing Problem

In this subsection, we introduce the settings on swing options and formulate the corresponding pricing problem. To begin with, we introduce mathematical notations and use them to sketch the details of a swing option from the holder’s view.Given a time interval[0,T], n points are selected from it to be the possible exercise times, which offer opportunities for the holder to exercise the swing option. We denote these days by τ0,τ1⋯τn-1, where 0<τ0<⋯<τn-1=T<∞. Let Δτl≔τl+1-τl(l=0,1,⋯,n-2), and for simplicity we set Δτl to be equal, i.e., Δτk=Δτ for l=0,1,⋯,n-2. For convenience, we use k instead of τk to be the subscript of all variables in this paper.At timeτk, the holder has one right to trade uk volume(s) of electricity with the strike price Ek, where uk usually is limited. We define(7)Nk≔∑i=0k1ui>0,Uk≔∑i=0kui,k=0,1,⋯,n-1, which denote the number of rights used and the total volume delivered respectively up to τk+. In particularly, we let N-1=0 and U-1=0, which means that the swing option has not been exercised before time τ0.Now we can introduce the main constraints in the contract of the swing option. In general, there are four kinds as follows:(1)
Total number of rights constraint. The total number of rights for holder to exercise is limited to NA∈N+, and it should be less than or equal to the number of possible exercise days, i.e., NA≤n.(2)
Refraction time constraint. Right after exercising, the holder should wait for at least the time period ΔtR(we call it refraction time) to exercise another right again. Here, we assume that ΔtR=κΔτ, where κ is some predetermined positive integer.(3)
Local volume constraint. For k=0,⋯,n-1, the delivery volume uk∈A≔{-m1,-m1+1,⋯,0,⋯,m2-1,m2}, where m1,m2∈N+.(4)
Global volume constraint. Once the holder’s total delivery volume Un-1 is beyond the global volume constraint (Vmin and Vmax), she will be penalized. And we adopt the general type penalty function introduced in [3], which is defined as(8)ψUn-1=-C1ifUn-1<Vmin,0ifVmin≤Un-1≤Vmax,-C2U-VmaxifUn-1>Vmax,whereC1 and C2 are positive constants.In this paper, we consider the valuation problem of the swing option with model (5) and (6) under all the above constraints, e.g., the total number of rights constraint, the refraction time constraint, the local volume constraint, and the global volume constraint. Now we formulate the pricing problem. Firstly, we introduce the strategy of exercising swing option in Definition 1.Definition 1.
The sequence{uk}k=0n-1 is defined as a strategy if it is Fk-adapted and satisfies the following conditions: (1)
Nk≤n, k=0,1,⋯,n-1.(2)
Ifui,uj>0, then τi-τj≥κΔτ, i≠j.(3)
uk∈A≔{-m1,-m1+1,⋯,0,⋯,m2-1,m2}, where m1,m2∈N+, k=0,1,⋯,n-1.
Moreover, the set of all admissible strategies is called admissible set and is denoted byD.Secondly, we define the value function in Definition2, which is widely used in stochastic control theory.Definition 2.
Given the risk-neutral probability measureQ, at time τj-(j∈{1,2,⋯,n-1}), the value of swing option with N¯j-1 right(s) used and U¯j-1 volume(s) totally delivered is given by(9)VjN¯j-1,U¯j-1s,i=supukk=jn-1∈DEQ∑k=jn-1e-rfk-jΔτ·ukϕSk,Ek+e-rfn-1-jΔτ·ψUn-1∣Sj=s,αj=i,Nj-1=N¯j-1,Uj-1=U¯j-1,where s∈R, i∈M, N¯j-1∈{0,1,⋯,NA}, U¯j-1∈N, rf represents the risk-free rate, and ϕ(Sk,Ek)≔max(Sk-Ek,Ek-Sk,0).It is pointed out that our final goal is to solve this primal optimal stochastic control problem; that is, we need to findV0N-1,U-1(s,i) and the optimal strategy {uk}k=0n-1. Note that if we delete the local constraint and global constraint and set NA=1, then this problem is pricing a Bermuda option. In addition, if we set NA=n, then it becomes the problem of pricing a strip of European options, and the differences of these European option are only embodied in expiry times and strike prices.We deal with the primal optimal stochastic control problem via the method of backward induction as follows:(i) IfNk-1=NA for k∈{n-1,n-2,⋯,0}, then(10)VkNk-1,Uk-1Sk,αk=e-rfn-1-kΔτ·ψUk-1;(ii) If Nk-1∈{NA-1,NA-2,⋯,0} for k=n-1, then(11)VkNk-1,Uk-1Sk,αk=supuk∈Auk·ϕSk,Ek+ψUk-1+uk;(iii) If Nk-1∈{NA-1,NA-2,⋯,0} for k∈{n-2,⋯,n-κ}, then(12)VkNk-1,Uk-1Sk,αk=supuk∈Auk·ϕSk,Ek+1uk=0·e-rfΔτEQVk+1Nk-1,Uk-1Sk+1,αk+1∣Fk+1uk>0·e-rfn-1-kΔτψUk-1+uk; (iv) If Nk-1∈{NA-1,NA-2,⋯,0} for k∈{n-κ-1,⋯,0}, then(13)VkNk-1,Uk-1Sk,αk=supuk∈Auk·ϕSk,Ek+e-rfk~-kΔτEQVk~Nk-1+1uk>0,Uk-1+ukSk~,αk~∣Fk,where (14)k~=k+1ifuk=0,k+κifuk>0. With (10)-(13), we can adopt tree method to fulfill our goal of pricing. Therefore, our next step is to build trinomial tree for St.


## 2.1. Characteristics of Electricity Spot Prices

Our modelling is based on the electricity spot price data of the day-ahead market in Nord Pool, and we focus on the stochastic component of the deseasonalized prices. Firstly, the stochastic component has the mean-reverting feature captured by Lucia and Schwartz [29], which shows that the electricity prices tend to revert to the mean value over time. Secondly, the time series of electricity prices show spikes, which can sometimes be depicted by jump models or regime-switching models (see [30]). All these two characteristics of the deseasonalized prices are the basis of modelling in this paper.


## 2.2. Models for Electricity Spot Prices

Throughout this paper, for a fixed maturityT>0, we work on a probability space (Ω,F,{Ft}0≤t≤T,P), where {Ft}0≤t≤T is the natural filtration generated by Brownian Motion {Bt}0≤t≤T in physical measure P and an independent observable continuous-time Markov chain {αt}0≤t≤T. Assume that {αt}0≤t≤T takes values in a finite state space M:={1,2,⋯,m0} and is generated by Q={qij}, where qij≥0 for i,j∈M,j≠i, and ∑j=1m0‍qij=0 for each i∈M. The transitional probability is given by(1)PijΔ≔Pαt+Δ=j∣αt=i=qijΔ+oΔ,j≠i,1+qijΔ+oΔ,j=i,for a small time interval Δ>0. Besides, we denote by ρi the probability of being in state i at the initial time t=0, i.e., (2)ρi≔Pα0=i,i∈M.We start our modelling for the spot price from the one-factor model used in [29], which captures seasonal and mean-reversion feature. The model is(3)St=exp⁡ft+Xt,(4)dXt=-βXtdt+σdBtP.In (3), the logarithm of the electricity price St is divided into two parts, one of which is the seasonal part f(t), and the other is the deseasonalized part Xt. Equation (4) describes the dynamic of the deseasonalized part by Vasicek model, where the long-term mean is set to be zero, and the speed of reversion together with instantaneous volatility is denoted by β and σ respectively, both of which are positive constants.In this paper, our model is(5)St=ft·Xt,(6)dXt=βαtξαt-Xtdt+σαtdBtP.In (5), the modification is based on the fact that the electricity price can be negative (see [5]). And in (6), to capture the feature that electricity price show spikes, we replace (4) with a n-state regime-switching mean-reverting model, which was proposed in [8, 10]. Specifically, β(·):M→R+, ξ(·):M→R, and σ(·):M→R+ represent the speed of adjustment, the long-term mean level and the volatility, respectively. Note that, for simplicity, we denote βαt≔β(αt),ξαt≔ξ(αt),σαt≔σ(αt).


## 2.3. Swings and Pricing Problem

In this subsection, we introduce the settings on swing options and formulate the corresponding pricing problem. To begin with, we introduce mathematical notations and use them to sketch the details of a swing option from the holder’s view.Given a time interval[0,T], n points are selected from it to be the possible exercise times, which offer opportunities for the holder to exercise the swing option. We denote these days by τ0,τ1⋯τn-1, where 0<τ0<⋯<τn-1=T<∞. Let Δτl≔τl+1-τl(l=0,1,⋯,n-2), and for simplicity we set Δτl to be equal, i.e., Δτk=Δτ for l=0,1,⋯,n-2. For convenience, we use k instead of τk to be the subscript of all variables in this paper.At timeτk, the holder has one right to trade uk volume(s) of electricity with the strike price Ek, where uk usually is limited. We define(7)Nk≔∑i=0k1ui>0,Uk≔∑i=0kui,k=0,1,⋯,n-1, which denote the number of rights used and the total volume delivered respectively up to τk+. In particularly, we let N-1=0 and U-1=0, which means that the swing option has not been exercised before time τ0.Now we can introduce the main constraints in the contract of the swing option. In general, there are four kinds as follows:(1)
Total number of rights constraint. The total number of rights for holder to exercise is limited to NA∈N+, and it should be less than or equal to the number of possible exercise days, i.e., NA≤n.(2)
Refraction time constraint. Right after exercising, the holder should wait for at least the time period ΔtR(we call it refraction time) to exercise another right again. Here, we assume that ΔtR=κΔτ, where κ is some predetermined positive integer.(3)
Local volume constraint. For k=0,⋯,n-1, the delivery volume uk∈A≔{-m1,-m1+1,⋯,0,⋯,m2-1,m2}, where m1,m2∈N+.(4)
Global volume constraint. Once the holder’s total delivery volume Un-1 is beyond the global volume constraint (Vmin and Vmax), she will be penalized. And we adopt the general type penalty function introduced in [3], which is defined as(8)ψUn-1=-C1ifUn-1<Vmin,0ifVmin≤Un-1≤Vmax,-C2U-VmaxifUn-1>Vmax,whereC1 and C2 are positive constants.In this paper, we consider the valuation problem of the swing option with model (5) and (6) under all the above constraints, e.g., the total number of rights constraint, the refraction time constraint, the local volume constraint, and the global volume constraint. Now we formulate the pricing problem. Firstly, we introduce the strategy of exercising swing option in Definition 1.Definition 1.
The sequence{uk}k=0n-1 is defined as a strategy if it is Fk-adapted and satisfies the following conditions: (1)
Nk≤n, k=0,1,⋯,n-1.(2)
Ifui,uj>0, then τi-τj≥κΔτ, i≠j.(3)
uk∈A≔{-m1,-m1+1,⋯,0,⋯,m2-1,m2}, where m1,m2∈N+, k=0,1,⋯,n-1.
Moreover, the set of all admissible strategies is called admissible set and is denoted byD.Secondly, we define the value function in Definition2, which is widely used in stochastic control theory.Definition 2.
Given the risk-neutral probability measureQ, at time τj-(j∈{1,2,⋯,n-1}), the value of swing option with N¯j-1 right(s) used and U¯j-1 volume(s) totally delivered is given by(9)VjN¯j-1,U¯j-1s,i=supukk=jn-1∈DEQ∑k=jn-1e-rfk-jΔτ·ukϕSk,Ek+e-rfn-1-jΔτ·ψUn-1∣Sj=s,αj=i,Nj-1=N¯j-1,Uj-1=U¯j-1,where s∈R, i∈M, N¯j-1∈{0,1,⋯,NA}, U¯j-1∈N, rf represents the risk-free rate, and ϕ(Sk,Ek)≔max(Sk-Ek,Ek-Sk,0).It is pointed out that our final goal is to solve this primal optimal stochastic control problem; that is, we need to findV0N-1,U-1(s,i) and the optimal strategy {uk}k=0n-1. Note that if we delete the local constraint and global constraint and set NA=1, then this problem is pricing a Bermuda option. In addition, if we set NA=n, then it becomes the problem of pricing a strip of European options, and the differences of these European option are only embodied in expiry times and strike prices.We deal with the primal optimal stochastic control problem via the method of backward induction as follows:(i) IfNk-1=NA for k∈{n-1,n-2,⋯,0}, then(10)VkNk-1,Uk-1Sk,αk=e-rfn-1-kΔτ·ψUk-1;(ii) If Nk-1∈{NA-1,NA-2,⋯,0} for k=n-1, then(11)VkNk-1,Uk-1Sk,αk=supuk∈Auk·ϕSk,Ek+ψUk-1+uk;(iii) If Nk-1∈{NA-1,NA-2,⋯,0} for k∈{n-2,⋯,n-κ}, then(12)VkNk-1,Uk-1Sk,αk=supuk∈Auk·ϕSk,Ek+1uk=0·e-rfΔτEQVk+1Nk-1,Uk-1Sk+1,αk+1∣Fk+1uk>0·e-rfn-1-kΔτψUk-1+uk; (iv) If Nk-1∈{NA-1,NA-2,⋯,0} for k∈{n-κ-1,⋯,0}, then(13)VkNk-1,Uk-1Sk,αk=supuk∈Auk·ϕSk,Ek+e-rfk~-kΔτEQVk~Nk-1+1uk>0,Uk-1+ukSk~,αk~∣Fk,where (14)k~=k+1ifuk=0,k+κifuk>0. With (10)-(13), we can adopt tree method to fulfill our goal of pricing. Therefore, our next step is to build trinomial tree for St.


## 3. Valuation of Swing Options in Trinomial Forest

Since that the trinomial forest method for pricing strong path-depend options is simple and fast even when the number of regime state is large (see [31]), we adopt it to calculate the conditional expectation in inductive equations. Note that in (10)-(13), expectations are taken under risk-neutral measure Q. Therefore, we start from finding the risk-neutral measures in each regime. Let(15)BtQ=BtP+μσit,i∈M,t∈0,T,where μ is an unknown constant, which can be calibrated on market data. By Girsanov Theorem, we know that BtQ is Brownian Motion under the risk-neutral measure Q, which satisfies (16)dQdP=exp-∫0TμσidBsQ-12∫0Tμσi2ds,i∈M. Then we can rewrite SDE (6) as(17)dXt=βαtξαt∗-Xtdt+σαtdBtQ, with X0=x0 and ξαt∗≔ξαt-μ/βαt.


### 3.1. Tree-Building Procedure for Deseasonalized Prices

In order to build a tree forXt based on the SDE (17), we discretize the time and the space. Given time interval [0,T], we select a time sequence 0<t0<t1<⋯<tn-1=T with equal step Δt=(T-t0)/n~. Similarly, we let the space size be positive constant σ¯>0. Both Δt and σ¯ are to be chosen later.For each node in the tree, we need to know the value and regime state of this node. Hence, we introduce the pair(Xk,αk)≔(Xtk, αtk), k=0,1,⋯,n-1 to convey the necessary information on all of the nodes in the tree.We now show the tree construction, which is initially designed by Liu [32]. For the node with (Xk,αk)=(x,i), we consider how it evolves in the next time. Firstly, the Markov chain evolves from αk=i to αk+1=j basing on the transitional probability defined in (1). Secondly, Xk+1 may take three values denoted by xk+1t, xk+1m, and xk+1b, where the superscripts t, m, and b are short for "top", "middle", and "bottom". Hence, there are three branches emanating from Xk=x. Besides, Liu [32] designed three structures of the branches, and in different structures xk+1t, xk+1m, and xk+1b are also different. Choosing which structure depends on the comparison result between x and the two threshold values xi,1 and xi,2 as follows:Structure(a): If xi,1≤x≤xi,2, then the possible values of Xk+1 in regime αk+1 are up to x+liσ¯Δt,remaining at x, and down to x-liσ¯Δt;Structure(b): If x<xi,1, then the possible values of Xk+1 in regime αk+1 are up to x+2liσ¯Δt, up to x+liσ¯Δt, and remaining at x;Structure(c): If x>xi,2, then the possible values of Xk+1 in regime αk+1 are remaining at x, down to x-liσ¯Δt, and down to x-2liσ¯Δt.Hereli depending on current state i is constant to be chosen later. And xi,1,xi,2 are defined as(18)xi,1=ξi∗-liσ¯-liσ¯2-σi2βiΔt,xi,2=ξi∗+liσ¯-liσ¯2-σi2βiΔt,i∈M,where the parameters are as before. It is emphasized that xi,1 and xi,2 change with the state of regime i. More details on this mechanism of the tree designing can be found in [32].In different structures, we introduce the notationsPx,it,j,Px,im,j, and Px,ib,j, which denote the conditional probabilities transforming from node (Xk,αk)=(x,i) to the nodes (Xk+1,αk+1)=(xk+1t,j), (Xk+1,αk+1)=(xk+1m,j), and (Xk+1,αk+1)=(xk+1b,j), respectively. We can calculate them by matching the mean E[Xk+1-Xk∣Xk=x] and the variance E[(Xk+1-Xk)2∣Xk=x] derived from tree with their counterpart derived from SDE (6). The explicit formula for Px,it,j,Px,im,j, and Px,ib,j in different structures is as follows, and for simplicity we write Pij=Pij(Δt).Structure (a) (19) P x , i t , j = P i j · σ i 2 + λ i x l i σ ¯ Δ t + λ i 2 x Δ t 2 l i σ ¯ 2 , P x , i m , j = P i j · 1 - σ i 2 + λ i 2 x Δ t l i σ ¯ 2 , P x , i b , j = P i j - P x , i t , j - P x , i m , j ,Structure (b) (20) P x , i t , j = P i j · σ i 2 - λ i x l i σ ¯ Δ t + λ i 2 x 2 l i σ ¯ 2 , P x , i m , j = P i j · - σ i 2 - 2 λ i x l i σ ¯ Δ t + λ i 2 x Δ t l i σ ¯ 2 , P x , i b , j = P i j - P x , i t , j - P x , i m , j ,Structure (c) (21) P x , i t , j = P i j · 1 + σ i 2 + 3 λ i x l i σ ¯ Δ t + λ i 2 x Δ t 2 l i σ ¯ 2 , P x , i m , j = P i j · - σ i 2 + 2 λ i x l i σ ¯ Δ t + λ i 2 x Δ t l i σ ¯ 2 , P x , i b , j = P i j - P x , i t , j - P x , i m , j ,where λi(x)=βi(ξi∗-x). To keep these probabilities positive, it is required that the parameters Δt,σ¯,li(i∈M) satisfy the following inequalities:(22)2σi3≤liσ¯≤2σi,i∈M,(23)Δt≤min1≤i≤m021-σi/liσ¯2βi.


### 3.2. The Valuation Formula of the Swing Option in Trinomial Forest

Note that{τk}k=0n-1 and {tk}k=0n-1 are two partitions of [0,T], but for simplicity, we let tk=τk for k=0,⋯,n-1 and assume that Δt=Δτ, and Δτ satisfies (22) and (23).Now we can write the recursive value functions with the trinomial tree model forXt as follows:(i) IfNk-1=NA for k∈{n-1,n-2,⋯,0}, then(24)VkNk-1,Uk-1Sk,αk=e-rfn-1-kΔτ·ψUk-1;(ii) If Nk-1∈{NA-1,NA-2,⋯,0} for k=n-1, then(25)VkNk-1,Uk-1Sk,αk=supuk∈Auk·ϕSk,Ek+ψUk-1+uk;(iii) If Nk-1∈{NA-1,NA-2,⋯,0} for k∈{n-2,⋯,n-κ}, then(26)VkNk-1,Uk-1Sk,αk=supuk∈Auk·ϕSk,Ek+1uk=0·e-rfΔτ∑γk=t,m,b∑αk+1=1m0Pxk,αkγk,αk+1·Vk+1Nk-1,Uk-1S¯k+1,αk+1+1uk>0·e-rfn-1-kΔτψUk-1+uk; (iv) If Nk-1∈{NA-1,NA-2,⋯,0} for k∈{n-κ-1,⋯,0}, then(27)VkNk-1,Uk-1Sk,αk=supuk∈Auk·ϕSk,Ek+e-rfk~-kΔτ∑γk=t,m,b∑αk+1=1m0⋯∑γk~-1=t,m,b∑αk~=1m0∏ω=kk~-1Pxωγω-1,αω-1,αωγω,αω+1·Vk~Nk-1+1uk>0,Uk-1+ukS¯k~,αk~, where (28)k~=k+1ifuk=0,k+κifuk>0.We need to explainS¯k+1 in (26), xωγω-1,αω-1forω∈{k,⋯,k~-1} in (27), and S¯k~ in (27). Note that S¯k+1 is a special case of S¯k~; hence we firstly consider xωγω-1,αω-1forω∈{k,⋯,k~-1}. When ω=k, we define xωγω-1,αω-1≔Xk, which is known at time τk. When ω>k, we define xωγω-1,αω-1 as follows:(i) ifxαω-1,1≤xω-1γω-2,αω-2≤xαω-1,2, then (29)xωγω-1,αω-1≔xω-1γω-2,αω-2+lαω-1σ¯Δτ,γω-1=t,xω-1γω-2,αω-2,γω-1=m,xω-1γω-2,αω-2-lαω-1σ¯Δτ,γω-1=b, (ii) if xω-1γω-2,αω-2<xαω-1,1, then (30)xωγω-1,αω-1≔xω-1γω-2,αω-2+2lαω-1σ¯Δτ,γω-1=t,xω-1γω-2,αω-2+lαω-1σ¯Δτ,γω-1=m,xω-1γω-2,αω-2,γω-1=b, (iii) if xω-1γω-2,αω-2>xαω-1,2, then (31)xωγω-1,αω-1≔xω-1γω-2,αω-2,γω-1=t,xω-1γω-2,αω-2-lαω-1σ¯Δτ,γω-1=m,xω-1γω-2,αω-2-2lαω-1σ¯Δτ,γω-1=b. Secondly, we explain S¯k~ as follows:(i) ifxαk~-1,1≤xk~-1γk~-2,αk~-2≤xαk~-1,2, then (32)S¯k~≔xk~-1γk~-2,αk~-2+lαk~-1σ¯Δτfτk~,γk~-1=t,xk~-1γk~-2,αk~-2fτk~,γk~-1=m,xk~-1γk~-2,αk~-2-lαk~-1σ¯Δτfτk~,γk~-1=b, (ii) if xk~-1γk~-2αk~-2<xαk~-1,1, then (33)S¯k~≔xk~-1γk~-2,αk~-2+2lαk~-1σ¯Δτfτk~,γk~-1=t,xk~-1γk~-2,αk~-2+lαk~-1σ¯Δτfτk~,γk~-1=m,xk~-1γk~-2,αk~-2fτk~,γk~-1=b, (iii) if xk~-1γk~-2αk~-2>xαk~-1,2, then (34)S¯k~≔xk~-1γk~-2αk~-2fτk~,γk~-1=t,xk~-1γk~-2,αk~-2-lαk~-1σ¯Δτfτk~,γk~-1=m,xk~-1γk~-2,αk~-2-2lαk~-1σ¯Δτfτk~,γk~-1=b.In Figure4, we illustrate the pricing procedure depicted by formulas (24)-(27).


## 3.1. Tree-Building Procedure for Deseasonalized Prices

In order to build a tree forXt based on the SDE (17), we discretize the time and the space. Given time interval [0,T], we select a time sequence 0<t0<t1<⋯<tn-1=T with equal step Δt=(T-t0)/n~. Similarly, we let the space size be positive constant σ¯>0. Both Δt and σ¯ are to be chosen later.For each node in the tree, we need to know the value and regime state of this node. Hence, we introduce the pair(Xk,αk)≔(Xtk, αtk), k=0,1,⋯,n-1 to convey the necessary information on all of the nodes in the tree.We now show the tree construction, which is initially designed by Liu [32]. For the node with (Xk,αk)=(x,i), we consider how it evolves in the next time. Firstly, the Markov chain evolves from αk=i to αk+1=j basing on the transitional probability defined in (1). Secondly, Xk+1 may take three values denoted by xk+1t, xk+1m, and xk+1b, where the superscripts t, m, and b are short for "top", "middle", and "bottom". Hence, there are three branches emanating from Xk=x. Besides, Liu [32] designed three structures of the branches, and in different structures xk+1t, xk+1m, and xk+1b are also different. Choosing which structure depends on the comparison result between x and the two threshold values xi,1 and xi,2 as follows:Structure(a): If xi,1≤x≤xi,2, then the possible values of Xk+1 in regime αk+1 are up to x+liσ¯Δt,remaining at x, and down to x-liσ¯Δt;Structure(b): If x<xi,1, then the possible values of Xk+1 in regime αk+1 are up to x+2liσ¯Δt, up to x+liσ¯Δt, and remaining at x;Structure(c): If x>xi,2, then the possible values of Xk+1 in regime αk+1 are remaining at x, down to x-liσ¯Δt, and down to x-2liσ¯Δt.Hereli depending on current state i is constant to be chosen later. And xi,1,xi,2 are defined as(18)xi,1=ξi∗-liσ¯-liσ¯2-σi2βiΔt,xi,2=ξi∗+liσ¯-liσ¯2-σi2βiΔt,i∈M,where the parameters are as before. It is emphasized that xi,1 and xi,2 change with the state of regime i. More details on this mechanism of the tree designing can be found in [32].In different structures, we introduce the notationsPx,it,j,Px,im,j, and Px,ib,j, which denote the conditional probabilities transforming from node (Xk,αk)=(x,i) to the nodes (Xk+1,αk+1)=(xk+1t,j), (Xk+1,αk+1)=(xk+1m,j), and (Xk+1,αk+1)=(xk+1b,j), respectively. We can calculate them by matching the mean E[Xk+1-Xk∣Xk=x] and the variance E[(Xk+1-Xk)2∣Xk=x] derived from tree with their counterpart derived from SDE (6). The explicit formula for Px,it,j,Px,im,j, and Px,ib,j in different structures is as follows, and for simplicity we write Pij=Pij(Δt).Structure (a) (19) P x , i t , j = P i j · σ i 2 + λ i x l i σ ¯ Δ t + λ i 2 x Δ t 2 l i σ ¯ 2 , P x , i m , j = P i j · 1 - σ i 2 + λ i 2 x Δ t l i σ ¯ 2 , P x , i b , j = P i j - P x , i t , j - P x , i m , j ,Structure (b) (20) P x , i t , j = P i j · σ i 2 - λ i x l i σ ¯ Δ t + λ i 2 x 2 l i σ ¯ 2 , P x , i m , j = P i j · - σ i 2 - 2 λ i x l i σ ¯ Δ t + λ i 2 x Δ t l i σ ¯ 2 , P x , i b , j = P i j - P x , i t , j - P x , i m , j ,Structure (c) (21) P x , i t , j = P i j · 1 + σ i 2 + 3 λ i x l i σ ¯ Δ t + λ i 2 x Δ t 2 l i σ ¯ 2 , P x , i m , j = P i j · - σ i 2 + 2 λ i x l i σ ¯ Δ t + λ i 2 x Δ t l i σ ¯ 2 , P x , i b , j = P i j - P x , i t , j - P x , i m , j ,where λi(x)=βi(ξi∗-x). To keep these probabilities positive, it is required that the parameters Δt,σ¯,li(i∈M) satisfy the following inequalities:(22)2σi3≤liσ¯≤2σi,i∈M,(23)Δt≤min1≤i≤m021-σi/liσ¯2βi.


## 3.2. The Valuation Formula of the Swing Option in Trinomial Forest

Note that{τk}k=0n-1 and {tk}k=0n-1 are two partitions of [0,T], but for simplicity, we let tk=τk for k=0,⋯,n-1 and assume that Δt=Δτ, and Δτ satisfies (22) and (23).Now we can write the recursive value functions with the trinomial tree model forXt as follows:(i) IfNk-1=NA for k∈{n-1,n-2,⋯,0}, then(24)VkNk-1,Uk-1Sk,αk=e-rfn-1-kΔτ·ψUk-1;(ii) If Nk-1∈{NA-1,NA-2,⋯,0} for k=n-1, then(25)VkNk-1,Uk-1Sk,αk=supuk∈Auk·ϕSk,Ek+ψUk-1+uk;(iii) If Nk-1∈{NA-1,NA-2,⋯,0} for k∈{n-2,⋯,n-κ}, then(26)VkNk-1,Uk-1Sk,αk=supuk∈Auk·ϕSk,Ek+1uk=0·e-rfΔτ∑γk=t,m,b∑αk+1=1m0Pxk,αkγk,αk+1·Vk+1Nk-1,Uk-1S¯k+1,αk+1+1uk>0·e-rfn-1-kΔτψUk-1+uk; (iv) If Nk-1∈{NA-1,NA-2,⋯,0} for k∈{n-κ-1,⋯,0}, then(27)VkNk-1,Uk-1Sk,αk=supuk∈Auk·ϕSk,Ek+e-rfk~-kΔτ∑γk=t,m,b∑αk+1=1m0⋯∑γk~-1=t,m,b∑αk~=1m0∏ω=kk~-1Pxωγω-1,αω-1,αωγω,αω+1·Vk~Nk-1+1uk>0,Uk-1+ukS¯k~,αk~, where (28)k~=k+1ifuk=0,k+κifuk>0.We need to explainS¯k+1 in (26), xωγω-1,αω-1forω∈{k,⋯,k~-1} in (27), and S¯k~ in (27). Note that S¯k+1 is a special case of S¯k~; hence we firstly consider xωγω-1,αω-1forω∈{k,⋯,k~-1}. When ω=k, we define xωγω-1,αω-1≔Xk, which is known at time τk. When ω>k, we define xωγω-1,αω-1 as follows:(i) ifxαω-1,1≤xω-1γω-2,αω-2≤xαω-1,2, then (29)xωγω-1,αω-1≔xω-1γω-2,αω-2+lαω-1σ¯Δτ,γω-1=t,xω-1γω-2,αω-2,γω-1=m,xω-1γω-2,αω-2-lαω-1σ¯Δτ,γω-1=b, (ii) if xω-1γω-2,αω-2<xαω-1,1, then (30)xωγω-1,αω-1≔xω-1γω-2,αω-2+2lαω-1σ¯Δτ,γω-1=t,xω-1γω-2,αω-2+lαω-1σ¯Δτ,γω-1=m,xω-1γω-2,αω-2,γω-1=b, (iii) if xω-1γω-2,αω-2>xαω-1,2, then (31)xωγω-1,αω-1≔xω-1γω-2,αω-2,γω-1=t,xω-1γω-2,αω-2-lαω-1σ¯Δτ,γω-1=m,xω-1γω-2,αω-2-2lαω-1σ¯Δτ,γω-1=b. Secondly, we explain S¯k~ as follows:(i) ifxαk~-1,1≤xk~-1γk~-2,αk~-2≤xαk~-1,2, then (32)S¯k~≔xk~-1γk~-2,αk~-2+lαk~-1σ¯Δτfτk~,γk~-1=t,xk~-1γk~-2,αk~-2fτk~,γk~-1=m,xk~-1γk~-2,αk~-2-lαk~-1σ¯Δτfτk~,γk~-1=b, (ii) if xk~-1γk~-2αk~-2<xαk~-1,1, then (33)S¯k~≔xk~-1γk~-2,αk~-2+2lαk~-1σ¯Δτfτk~,γk~-1=t,xk~-1γk~-2,αk~-2+lαk~-1σ¯Δτfτk~,γk~-1=m,xk~-1γk~-2,αk~-2fτk~,γk~-1=b, (iii) if xk~-1γk~-2αk~-2>xαk~-1,2, then (34)S¯k~≔xk~-1γk~-2αk~-2fτk~,γk~-1=t,xk~-1γk~-2,αk~-2-lαk~-1σ¯Δτfτk~,γk~-1=m,xk~-1γk~-2,αk~-2-2lαk~-1σ¯Δτfτk~,γk~-1=b.In Figure4, we illustrate the pricing procedure depicted by formulas (24)-(27).


## 4. Empirical Calibration and the Market Price of Risk

Before implementing numerical examples of pricing swing options with the above trinomial forest approach, we do calibration for model (5) and (6). Our dataset covers the period from 2008/1/1 to 2017/11/16; therefore in this section, we let t∈{0,1,⋯,3607}, where 0 denotes 2008/1/1, 3597 denotes 2017/11/16, and the time step 1 denotes one day. We set the states of the regime to be 2, and what we need to calibrate are f(t) and θ≔(βi,ξi,σi,Pii,ρi),i=1,2. After getting calibration results, we compare the forecasting and actual data. In addition, we also estimate the market price of risk.


### 4.1. Empirical Calibration

In this subsection we carry out calibration. Firstly, we obtain the system daily average electricity price in Nord Pool market from the Bloomberg service, which is presented in Figure1. Then we deal with the seasonal pattern f(t) by using the moving average and least squares fit techniques (see [4] for more details). We assume that f(t) has the weekly part fwe(t) together with the annually part fan(t) and that f(t)=fwe(t)·fan(t). Then we apply moving average technique to fwe(t); that is, we let fwe(t)=fwe(t-7) for t>7. And the annual cycle can be estimated with the sinusoidal form fant=Asin(2π(t+B)/365)+Ct+D via least squares fit. The results of fwe(t) and fan(t) are presented in Table 1 and (35), respectively.(35)fant=-4.7728sin2π365t+275.7260-0.0067t+47.8687Table 1
f w e ( t ) of one loop cycle.
t 0 1 2 3 4 5 6 fwe(t) 1.0405 1.0334 1.0140 0.9351 0.9018 1.0350 1.0402Figure 1
Daily average system electricity price.After removing the seasonal partf(t), we focus on the estimation of θ=(βi,ξi,σi,Pii,ρi) for i=1,2. To obtain θ^=(β^i,ξ^i,σ^i,P^ii,ρ^i),i=1,2, we discretize (6) to have that(36)Xt+1=Xt+βαtξαt-Xt+σαtϵt,where ϵt~N(0,1). Then by employing Expectation-Maximization (EM) algorithm (see [33–35]), we conduct estimation. And values of θ^ are shown in Table 2.Table 2
Calibration results for two-state regime-switching model.
Regime β ^ i ξ ^ i σ ^ i P ^ i i P ^ ( α 0 = i ) 1 0.1244 1.0460 0.1875 0.7739 1 2 0.0181 1.0195 0.0517 0.9574 0


### 4.2. Assessment of Forecasts

In this part, we assess the predicted yields by calculating the differences between the one-step ahead forecasts and the actual data, as conducted in [36]. In particular, we choose to use the mean square error (MSE), the mean absolute Error (MAE), the mean absolute percent error (MAPE) and the median relative absolute percentage error (MdAPE) as indicators. Firstly, the one-step ahead predicted Xt+1 are calculated by (37)EPXt+1∣Ft=EPβ^αtξ^αt+1-β^αtXt+σ^αtϵt∣Ft=∑i=12P^αt=i·β^iξ^i+1-β^iXt. With the calibrated parameters in Tables 1 and 2 and (35), we can calculate MSE, MAE, MAPE, and MdAPE, and results of these errors are presented in Table 3. Moreover, the actual data and the one-step ahead forecast are plotted in Figure 2, and we find that curves of them are very close to each other.Table 3
Errors for one-step ahead forecast.
Forecast error measures MSE MAE MAPE MdAPE Values 0.0081 0.0560 0.0617 0.0378Figure 2
Actual data and one-step ahead forecast ofXt.


### 4.3. The Market Price of Risk

Note that the fair price of the swing option in Definition2 is the discounted expected future payoff under a martingale measure Q. In order to find Q, we apply the risk premium approach proposed in [37, 38]. Define the risk premium as(38)RPT^≔EPST^∣F0-F0T^,where EP is the expectation under the physical measure P and F0T^ is the price of future contract at present time 0 with the delivery time T^. With an equation connecting the risk premium RP(T^) and the market price of risk μ emerging in (15), we can get an estimation of μ. All the calculation details are presented in Appendix A. Here we just show the result(39)μ^=-0.0032.


## 4.1. Empirical Calibration

In this subsection we carry out calibration. Firstly, we obtain the system daily average electricity price in Nord Pool market from the Bloomberg service, which is presented in Figure1. Then we deal with the seasonal pattern f(t) by using the moving average and least squares fit techniques (see [4] for more details). We assume that f(t) has the weekly part fwe(t) together with the annually part fan(t) and that f(t)=fwe(t)·fan(t). Then we apply moving average technique to fwe(t); that is, we let fwe(t)=fwe(t-7) for t>7. And the annual cycle can be estimated with the sinusoidal form fant=Asin(2π(t+B)/365)+Ct+D via least squares fit. The results of fwe(t) and fan(t) are presented in Table 1 and (35), respectively.(35)fant=-4.7728sin2π365t+275.7260-0.0067t+47.8687Table 1
f w e ( t ) of one loop cycle.
t 0 1 2 3 4 5 6 fwe(t) 1.0405 1.0334 1.0140 0.9351 0.9018 1.0350 1.0402Figure 1
Daily average system electricity price.After removing the seasonal partf(t), we focus on the estimation of θ=(βi,ξi,σi,Pii,ρi) for i=1,2. To obtain θ^=(β^i,ξ^i,σ^i,P^ii,ρ^i),i=1,2, we discretize (6) to have that(36)Xt+1=Xt+βαtξαt-Xt+σαtϵt,where ϵt~N(0,1). Then by employing Expectation-Maximization (EM) algorithm (see [33–35]), we conduct estimation. And values of θ^ are shown in Table 2.Table 2
Calibration results for two-state regime-switching model.
Regime β ^ i ξ ^ i σ ^ i P ^ i i P ^ ( α 0 = i ) 1 0.1244 1.0460 0.1875 0.7739 1 2 0.0181 1.0195 0.0517 0.9574 0


## 4.2. Assessment of Forecasts

In this part, we assess the predicted yields by calculating the differences between the one-step ahead forecasts and the actual data, as conducted in [36]. In particular, we choose to use the mean square error (MSE), the mean absolute Error (MAE), the mean absolute percent error (MAPE) and the median relative absolute percentage error (MdAPE) as indicators. Firstly, the one-step ahead predicted Xt+1 are calculated by (37)EPXt+1∣Ft=EPβ^αtξ^αt+1-β^αtXt+σ^αtϵt∣Ft=∑i=12P^αt=i·β^iξ^i+1-β^iXt. With the calibrated parameters in Tables 1 and 2 and (35), we can calculate MSE, MAE, MAPE, and MdAPE, and results of these errors are presented in Table 3. Moreover, the actual data and the one-step ahead forecast are plotted in Figure 2, and we find that curves of them are very close to each other.Table 3
Errors for one-step ahead forecast.
Forecast error measures MSE MAE MAPE MdAPE Values 0.0081 0.0560 0.0617 0.0378Figure 2
Actual data and one-step ahead forecast ofXt.


## 4.3. The Market Price of Risk

Note that the fair price of the swing option in Definition2 is the discounted expected future payoff under a martingale measure Q. In order to find Q, we apply the risk premium approach proposed in [37, 38]. Define the risk premium as(38)RPT^≔EPST^∣F0-F0T^,where EP is the expectation under the physical measure P and F0T^ is the price of future contract at present time 0 with the delivery time T^. With an equation connecting the risk premium RP(T^) and the market price of risk μ emerging in (15), we can get an estimation of μ. All the calculation details are presented in Appendix A. Here we just show the result(39)μ^=-0.0032.


## 5. Numerical Example: All Constraints Are Nontrivial

In this section, with all the results obtained in Section4, we implement numerical examples in the case that all constrains are nontrivial.The settings on swing option are as follows. The holding period is[0,T], where 0 represents the date 2008/1/1 and T=1 represents 2008/12/31. Set Δτ=1/3 (120 days), which implies that τ0, τ1, τ2, and τ3 representing 2008/1/6, 2008/5/5, 2008/9/2, and 2008/12/31, respectively. The refraction time (2Δτ) is 240 days, i.e., κ=2. The corresponding strike prices Ek(k=0,⋯n-1) are chosen to be 45 euros. The possible delivery volume(s) each time are 0 MWh, 1 MWh, or 2 MWh, respectively. The total number of rights NA is set to be 2. As for the penalty function, we adopt a specific form of function (8) as(40)ψUn-1=0if0≤Un-1≤3,-3Un-1-3ifUn-1>3.For the tree-building part, we setσ¯=0.1, l1=3, and l2=1 and use the information of Table 2. Firstly, by the method introduced in Section 3.1, we build a tree for Xt, which is presented in Figure 3(a). Secondly, we calculate the values of f(t) at 4 possible exercise days, which are approximately equal to 54.42, 46.05, 45.88, and 51.85, respectively. Thirdly, we use the relationship(41)St=Xt·ftto get a trinomial tree for St, which is shown in Figure 3(b).Figure 3
Trinomial tree for underlying variable.
(a)
Trinomial Tree forXt (b)
Trinomial Tree forStFigure 4
Connection Between the level withN rights exercised and the level with N+1 rights exercised in the Swing forest.Besides, we setrf=0.1. Under all these above settings and with the information listed in Table 2, we carry out a numerical example. The results are as follows.Firstly, we present the values of the swing option under different cases. To be specific, we present the case where the swing option has been used all its rights (2 rights) in Table 4, the case where the swing option has been used 1 right in Table 5 and the case where the swing option has not been used any right in Table 6. Secondly, we see from Table 6 the values of the swing option at time τ0 in two regimes, respectively, which are 19.2 under state 1 and 15.2 under state 2. Thirdly, we show details of calculating the bolded numbers in these tables. In Table 4, (42)-2.7=24e-rf·3Δτ·ψ4=e-0.1·3·1/3·-3·4-3. In Table 5, (43)60.2=25maxψ1,1·75.1-45+ψ1+1,2·75.1-45+ψ1+2,37.2=26maxe-rfΔτ0.13×57.2+0.47×39.2+0.18×21.2+0.04×57.2+0.14×39.2+0.05×21.2,1·max58.5-45,0+e-rfΔτ·ψ2+1,2·max58.5-45,0+e-rfΔτ·ψ2+2=max37.2,13.5,24.1. In Table 6, (44)19.2=27maxe-rfΔτ0.16×30.9+0.47×9.7+0.14×2.5+0.05×31.4+0.14×6.1+0.04×0.1,1·max50.5-45,0+e-rf·2Δτ∑i=136aibici,2·max50.5-45,0+e-rf·2Δτ∑i=136dieifi,=max11.9,15.6,19.2,where sequences {ai}, {bi}, {ci}, {di}, {ei}, and {fi} are presented in Appendix B.Table 4
Option values when2 rights are used (no rights left).
When2 volume delivered When3 volume delivered When4 volume delivered τ0 τ 1 τ 2 τ 3 τ0 τ 1 τ 2 τ 3 τ0 τ 1 τ 2 τ 3  0.0,0.0 0.0,0.0 0.0,0.0 0.0,0.0 0.0,0.0 0.0,0.0 0.0,0.0 0.0,0.0 -2.7,-2.7 -2.8,-2.8 - 2.9 , - 2.9 -3.0,-3.0 0.0,0.0 0.0,0.0 0.0,0.0  0.0,0.0 0.0,0.0 0.0,0.0  -2.8,-2.8 - 2.9 , - 2.9 -3.0,-3.0 0.0,0.0 0.0,0.0 0.0,0.0  0.0,0.0 0.0,0.0 0.0,0.0  -2.8,-2.8 - 2.9 , - 2.9 -3.0,-3.0 0.0,0.0 0.0,0.0 0.0,0.0  0.0,0.0 0.0,0.0 0.0,0.0  -2.8,-2.8 - 2.9 , - 2.9 -3.0,-3.0 0.0,0.0 0.0,0.0 0.0,0.0  0.0,0.0 0.0,0.0 0.0,0.0  -2.8,-2.8 - 2.9 , - 2.9 -3.0,-3.0 0.0,0.0 0.0,0.0   0.0,0.0 0.0,0.0   - 2.9 , - 2.9 -3.0,-3.0 0.0,0.0 0.0,0.0   0.0,0.0 0.0,0.0   - 2.9 , - 2.9 -3.0,-3.0 0.0,0.0 0.0,0.0   0.0,0.0 0.0,0.0   - 2.9 , - 2.9 -3.0,-3.0 0.0,0.0 0.0,0.0   0.0,0.0 0.0,0.0   - 2.9 , - 2.9 -3.0,-3.0 0.0,0.0 0.0,0.0   0.0,0.0 0.0,0.0   - 2.9 , - 2.9 -3.0,-3.0 0.0,0.0 0.0,0.0   0.0,0.0 0.0,0.0   - 2.9 , - 2.9 -3.0,-3.0 0.0,0.0    0.0,0.0    - 3.0 , - 3.0 0.0,0.0    0.0,0.0    - 3.0 , - 3.0 0.0,0.0    0.0,0.0    - 3.0 , - 3.0 0.0,0.0    0.0,0.0    -3.0,-3.0 0.0,0.0    0.0,0.0    - 3.0 , - 3.0 0.0,0.0    0.0,0.0    - 3.0 , - 3.0Table 5
Option values when1 right is used (1 right left).
When1 volume delivered When2 volume delivered τ0 τ 1 τ 2 τ 3 τ0 τ 1 τ 2 τ 3  11.0,11.0 22.2,22.8 40.2,40.5 60.2 , 60.2 8.3,8.3 19.5,20.0 37.2 , 37.6 57.2,57.2 13.6,11.6 28.7,29.0 48.2,48.2  11.4,8.8 25.8,26.1 45.2,45.2 9.7,6.1 22.9,23.5 42.2,42.2  7.6,3.9 20.0,20.6 39.2,39.2 5.8,1.8 17.5,17.7 36.2,36.2  4.8,1.0 15.1,14.8 33.2,33.2 2.4,0.1 13.2,12.0 30.2,30.2  1.6,0.0 10.8,9.0 27.2,27.2 8.7,6.1 24.2,24.2   6.3,3.6 21.2,21.2 4.1,1.1 18.2,18.2   3.4,0.6 15.2,15.2 2.6,0.1 12.4,12.4   2.0,0.8 9.4,9.4 1.4,0.0 6.4,6.4   0.8,0.0 3.4,3.4 0.1,0.0 0.4,0.4   0.0,0.0 0.2,0.2 0.0,0.0 0.0,0.0   0.0,0.0 0.0,0.0 0.0,0.0    0.0,0.0 0.0,0.0    0.0,0.0 0.0,0.0    0.0,0.0 0.0,0.0    0.0,0.0 0.0,0.0    0.0,0.0 0.0,0.0    0.0,0.0Table 6
Option values when0 rights are used (2 rights left).
When0 volume delivered τ 0 τ 1 τ 2 τ 3  19.2,15.2 30.9,31.4 40.2,40.5 60.2,60.2 14.0,12.0 28.7,29.0 48.2,48.2 9.7,6.1 22.9,23.5 42.2,42.2 5.8,1.8 17.5,17.7 36.2,36.2 2.5,0.1 13.2,12.0 30.2,30.2 8.7,6.1 24.2,24.2 4.1,1.1 18.2,18.2 2.6,0.1 12.4,12.4 1.4,0.0 6.4,6.4 0.1,0.0 0.4,0.4 0.0,0.0 0.0,0.0 0.0,0.0 0.0,0.0 0.0,0.0 0.0,0.0 0.0,0.0 0.0,0.0


## 6. Conclusions

Under the settings that the underlying asset follows a seasonal mean-reverting regime-switching model withn regimes and that some typical constraints are nontrivial, we have proposed a trinomial forest framework for the valuation of the swing option, which is a multiple-layer tree extension of the traditional trinomial tree approach. It is known that tree approach is based on the risk-neutral pricing theory. Thus we have introduced the risk premium to transform the physical measure into the risk-neutral one. Then the problem of pricing swing options has been stated as an optimal stochastic control problem, where complex constraints on options are inherited by admissible control. In light of stochastic control theory, we have figured out that finding the option price is equivalent to obtaining the associated value function, which satisfies a Bellman equation. On the basis of the above work, we have succeeded in stating pricing formula for swing options via backward induction.To show the details of pricing, we have conducted a numerical example based on the calibration results. Comparing to previous work, in our example we have added the refraction time constraint and a more general form global constraint to make the valuation more complicated but closer to the real trading of swing options. This example has also shown that the trinomial forest method does work when both the underlying asset process and the constraints embedded in swing contract are complex.



---

*Source: 5796921-2019-01-09.xml*