Evaluating Energy Forward Dynamics Modeled as a Subordinated Hilbert-Space Linear Functional

Show more

1. Introduction

The forward pricing dynamics of an incomplete Arrow Debreu world reveals interesting challenges to speculators and one such challenge is the stochastic nature of the return process for every investment in an underlying commodity stock. In this study, we define such return process as $X\left(t\right)$ and assume a frequently changing property of time series data in discrete space. Also, the pricing model of the return process is assumed to be impacted by the fluctuation in the price caused either by contago or normal backwardation. Hence, the evaluation of the forward dynamics of the underlying commodity while incorporating the time-change component as an incomplete energy market is very significant for this study.

Interestingly, we adopt the completeness properties of a Banach space (a special type of Hilbert space) such that the return process $X\left(t\right)$ is defined in some normed space (i.e. complete, without hole) to enable us capture all discrete moving forward rates in the corresponding forward curves of the pricing dynamics, as such adopting similar approach by [1] and [2]. This approach is considered such that the distance between two nodes defined by the daily change price on the curve is defined in norm spaces with no gap in the sequence ${X}_{t}\mathrm{,}{X}_{t+1}\mathrm{,}{X}_{t+2}\mathrm{,}\cdots \mathrm{,}{X}_{t+n}$ , where $n=1,2,3,\cdots $ with each representing a node in a forward curve of the energy forward contracts.

Furthermore, for our time-change evaluation process and according to [3], the usual subordination procedure can be used to generate a Banach space valued Lévy processes and with the dynamics from other commodities markets like power and gas we pick a motivation for this study due to the presence of strong seasonality patterns, high degree of idiosyncratic risk over different market segments and leading to spikes in the forwards contract curve. Similarly, in other to capture all forward driven prices and capture the corresponding rates, the return process will be moved from general Lévy process to a subordinated Lévy process in discrete time.

Meanwhile, the arbitrage-free forward price is derived from the buy-and-hold strategy in the underlying spot commodity and the forward price dynamics is thus implied from a given stochastic model of the spot commodity based upon similar method used by [4] and [5]. Therefore, the representation of the forward price as the conditional expected value of the spot at time of delivery is represented. As such the expectation is estimated with respect to an equivalent martingale measure Q which is only possible if the price of the spot commodity is specified by a semimartingale dynamics.

In literature ( [4] [6] [7] and [8]), it is clear that the fundamental relationship between the spot and forward is highly delicate in energy markets and it is only fair to model the forward price dynamics directly.

2. Preliminaries

2.1. Characteristics Function of a Subordinated Lévy Process

Definition 1 *A return process
${\left\{X\left(t\right)\right\}}_{t\ge 0}$ used in this study as Lévy process*,* with non-decreasing paths almost surely*,* are called subordinators and such processes can be thought of as random models of the evolution*.* *

A wide class of return processes appearing in application is obtained by subordination of Lévy processes with drift. Also, a stochastic time-change to the Lévy process amounts to stochastically altering the clock on which the Lévy process runs.

Meanwhile, for the return process $X\left(t\right)$ , the subordinated process is such that

$X\left(t\right)=L\left(\Theta \left(t\right)\right)$ (1)

where the subordinator $\Theta \left(t\right)$ and $L\left(t\right)$ are increasing Lévy processes with independent and stationary increments.

According to [9] every semi-martingale $X\left(t\right)$ can be written as a time-changed Brownian motion (Lévy process), where the random time $\Theta \left(t\right)$ is a positive and increasing semi-martingale. As a result, there exist an incremental and stationary Brownian motion $\left(B\left(u\right)\mathrm{,}u\ge 0\right)$ and a random time-change $\Theta \left(t\right)$ which forms an increasing stochastic process

$X\left(t\right)=B\left(\Theta \left(t\right)\right)$ (2)

Hence, every semi-martingale can also be written as a time-changed Lévy process such as

$X\left(t\right)={X}_{t}=L\left(\Theta \left(t\right)\right)$ (3)

The distribution of increments $\Delta L\left(\Theta \left(t\right)\right)$ is said to be subordinate to the distribution of increment $\Delta L\left(t\right)$ ; where $\Delta \left(t\right)$ is the data frequency or a clock measuring the speed of the evolution.

Meanwhile, in developing a subordinated Lévy process with varying time, the characteristic function of a subordinated process $X\left(t\right)$ is utilized which is obtained by composition of the Laplace exponent of $\Theta \left(t\right)$ with the characteristics exponent of $L\left(t\right)$ .

Suppose we let ${\left({\Theta}_{t}\right)}_{t\ge 0}$ be a subordinator that is a Lévy process whose trajectories are increasing. Since $\Theta \left(t\right)$ is a positive random variable for all t it is described by its Laplace transform rather than the Fourier transform. Let the characteristic triplet of $\Theta \left(t\right)$ be $\left(\Theta \mathrm{,}\rho \mathrm{,}\alpha \right)$ . Then the moment generating function of $\Theta \left(t\right)$ is given by;

$E\left[{\text{e}}^{u{\Theta}_{t}}\right]={\text{e}}^{t{L}^{*}\left(u\right)},\forall u\le 0$ (4)

where

${L}^{*}\left(u\right)=\alpha u+{\displaystyle {\int}_{0}^{\infty}}\left({\text{e}}^{ux}-1\right)\rho (\; d\; x\; )$

and ${L}^{*}\left(u\right)$ is the Laplace exponent of $\Theta \left(t\right)$ . Since the time-changed process ${X}_{t}={L}_{{\Theta}_{t}}=L\left({\Theta}_{t}\right)$ is a stochastic process evaluated at a stochastic time, its characteristics function involves expectation over two sources of randomness

${\Phi}_{{X}_{t}}\left(u\right)=E\left[{\text{e}}^{iu{L}_{{\Theta}_{t}}}\right]$

${\Phi}_{{X}_{t}}\left(u\right)=E\left[E\left[{\text{e}}^{iu{L}_{\u03f5}}|{\Theta}_{t}=\u03f5\right]\right]$

where the inside expectation is taken on ${L}_{{\Theta}_{t}}$ , conditional on a fixed value of ${\Theta}_{t}=\u03f5$ and the outside expectation is on all possible value of ${\Theta}_{t}$ . If the random time ${\Theta}_{t}$ is independent of ${L}_{t}$ , the randomness due to the Lévy process can be integrated out using moment generating form of ${L}_{t}$ given by

${\Phi}_{{L}_{t}}\left(u\right)=E\left[{\text{e}}^{iu{L}_{t}}\right]$

${\Phi}_{{L}_{t}}\left(u\right)=\left[{\text{e}}^{-i{\Psi}_{L}\left(u\right)}\right]$

and a characteristics function of ${X}_{t}$ , that is

${\Phi}_{{X}_{t}}\left(u\right)=E\left[{\text{e}}^{iu{L}_{{\Theta}_{t}}}\right]$

${\Phi}_{{X}_{t}}\left(u\right)=E\left[E\left[{\text{e}}^{iu{L}_{\u03f5}}|{\Theta}_{t}=\u03f5\right]\right]$

${\Phi}_{{X}_{t}}\left(u\right)=E{\text{e}}^{-i{\Psi}_{L}(\; u\; )}$

${\Phi}_{{X}_{t}}\left(u\right)={L}_{{\Theta}_{t}}^{*}\left({\Psi}_{L}(\; u\; )\right)$

Under independence, the characteristics function of ${X}_{t}={L}_{{\Theta}_{t}}$ is just the Laplace transform of ${\Theta}_{t}$ evaluated at the characteristics component of ${L}_{t}$ .

Let ${\left({L}_{t}\right)}_{t\ge 0}$ be a Lévy Process on $\mathbb{R}$ with characteristics exponent $\Psi \left(u\right)$ and triplet $\left(\sigma \mathrm{,}v\mathrm{,}\gamma \right)$ and let ${\left({\Theta}_{t}\right)}_{t\ge 0}$ be a subordinator with Laplace exponent ${L}^{\mathrm{*}}\left(u\right)$ and triplet $\left(\Theta \mathrm{,}\rho \mathrm{,}a\right)$ . Then the subordinated process $X\left(t\right)$ defined by;

$X\left(t\right)=L\left(\Theta \left(t\right)\right)\mathrm{,}$ (5)

is also a subordinated Lévy process with characteristics function

$E\left[{\text{e}}^{iu{X}_{t}}\right]={\text{e}}^{t{L}^{*}\left(\Psi \left(u\right)\right)}.$ (6)

Therefore, in defining a parametric Lévy process is to obtain an Lévy process by subordinating a Brownian motion with an independent increasing Lévy process. However, for a continuous time change process, we it is necessary to evaluate the stochastic time integral of the subordinated Lévy process linear in each argument.

2.2. Stochastic Integrals of a Subordinated Lévy Process in Hilbert Space

Since the time-changed return processes ${X}_{t}$ and ${Y}_{t}$ are independent subordinated Lévy processes evaluated at a stochastic time, and using the approach in [10], we derive some general results on the stochastic integral represented as linear functional in a subordinated Hilbert space. Suppose we define a subordinated stochastic integral Y in the form

$Y\left(t\right)={\displaystyle {\int}_{0}^{t}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\Psi \left(s\right)\text{d}L\left(s\right)\mathrm{,}$ (7)

for a Lévy process $L\left(t\right)$ with values in a separable Hilbert space U and an integral stochastic process $\Psi \mathrm{:}{\mathbb{R}}_{+}\times U\to H$ , and considering the stochastic partial differential equation

$\text{d}X\left(t\right)=\left(AX\left(t\right)+\alpha \left(t\right)\right)\text{d}t+\sigma \left(t\right)\text{d}L\left(t\right)\mathrm{,}t\in \mathbb{R}$ (8)

on the interval $\mathbb{R}=\left[\mathrm{0,}T\right]$ with H being a separable Hilbert spaces equipped with some boundary conditions where

· $A\mathrm{:}H\to H$ is a first order differential operator

· $\alpha \mathrm{:}\mathbb{R}\to H$ is a mapping with Brochner-integrable trajectories

· $L\mathrm{:}\mathbb{R}\to H$ is a square-integrable, H-valued stochastic process with covariance operator $Q\in {L}_{1}^{+}$ .

· $\sigma \mathrm{:}\mathbb{R}\to L\left({Q}^{1/2}\left(H\right)\mathrm{,}H\right)$ is an operator-valued process.

Suppose for convenience, we are only interested in one-dimensional martingales, i.e., we are interested in

$Y\left(t\right)=\Upsilon \left(Y\left(t\right)\right)\mathrm{,}$ (9)

where $\Upsilon $ is a continuous linear functional on the state space H of Y. Therefore, by replacing $Y\left(t\right)$ with $X\left(t\right)$ assuming the characteristic of the subordinated process we have;

$X\left(t\right)={\displaystyle {\int}_{0}^{t}}\left(\Upsilon \circ \Psi \left(s\right)\right)\text{d}L\left(s\right)\mathrm{.}$ (10)

Suppose we take U as finite dimensional, and some standard (real-valued) Brownian motion ${\left\{B\left(t\right)\right\}}_{t\ge 0}$ and Itô integral stochastic process $\sigma $ such that

$X\left(t\right)={\displaystyle {\int}_{0}^{t}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\sigma \left(s\right)\text{d}B\left(s\right)$ (11)

This leads us to Theorem 1 which shows that a similar representation holds if U is any separable Hilbert space.

Theorem 1. *Let
$n\in \mathbb{N}$ and H*,* *U be separable Hilbert spaces.* Let L be a square integrable and mean zero *U-valued Wiener process with coveriance
$Q\in {L}_{1}^{+}$ .* Assume that
$Dimran\left(Q\right)\ge n$ and Q is positive definite*.* Let
$\Psi \in {\aleph}_{L}^{2}$ *,
$\Upsilon \in L\left(H\mathrm{,}{\mathbb{R}}^{n}\right)$ * and define *

$X\left(t\right)=\Upsilon \left({\displaystyle {\int}_{0}^{t}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\Psi \left(s\right)\text{d}L\left(s\right)\right)\mathrm{.}$

Then there is an n-dimensional standard Brownian motion B such that

$X\left(t\right)\mathrm{=}{\displaystyle {\int}_{0}^{t}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\sigma \left(s\right)\text{d}B\left(s\right)\mathrm{,}$

where $\sigma \left(s\right)={\left(\Upsilon \Psi \left(s\right)Q\Psi {\left(s\right)}^{\mathrm{*}}{\Upsilon}^{\mathrm{*}}\right)}^{1/2}\in {\aleph}_{B}^{2}\left({\mathbb{R}}^{\mathrm{*}}\right)$ . If $\sigma \left(s\right)$ is invertible in ${\mathbb{R}}^{n\times n}$ for $\lambda \otimes P$ -almost any $s\in \mathbb{R}$ , then ${\sigma}^{-1}\in {\aleph}_{X}^{2}\left({\mathbb{R}}^{n}\right)$ and

$B\left(t\right)={\displaystyle {\int}_{0}^{t}}{\left(\sigma \left(s\right)\right)}^{-1}\text{d}X\left(s\right)\mathrm{.}$ (12)

*Proof*: Let
$\Psi $ be an elementary random variable such that there exist
$n\in \mathbb{N}$ ,
$0\le {a}_{j}\le {b}_{j}<\infty $ , and
${F}_{{a}_{j}}$ -measurable square integrable random variables
${Y}_{j}$ and
${\psi}_{j}\in L\left(U\mathrm{,}H\right)$ such that

$\Psi ={\displaystyle \underset{j=1}{\overset{n}{\sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{Y}_{j}{1}_{\left[{a}_{j},{b}_{j}\right]{\psi}_{j}}.$

By definition of L, we have $L\left(t\right)=W\left(\Theta \left(t\right)\right)$ for a U-valued Wiener process W with martingale covariance ${Q}^{W}$ and a subordinator $\Theta $ . Let $\Gamma $ be the isometric embedding given in Lemma (1). Note that $\Gamma \left(\Upsilon \Psi =\Upsilon \Gamma \left(\Psi \right)\right)$ because this holds if $\Psi $ is elementary. From Lemma (1) and Theorem (1) we have;

$\begin{array}{c}\Upsilon \left({\displaystyle {\int}_{0}^{t}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\Psi \left(s\right)\text{d}L\left(s\right)\right)={\displaystyle {\int}_{0}^{t}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\Upsilon \Psi \left(s\right)\text{d}L\left(s\right)\\ ={\displaystyle {\int}_{0}^{\Theta \left(t\right)}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\Upsilon \Gamma \left(\Psi \right)\left(s\right)\text{d}W\left(s\right)\\ ={\displaystyle {\int}_{0}^{\Theta \left(t\right)}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\sigma}_{1}\left(s\right)\text{d}B(\; s\; )\end{array}$

where B is a standard Brownian motion on ${\mathbb{R}}^{n}$ and

${\sigma}_{1}\left(t\right)=\Gamma \left({\left(\Upsilon \Psi {Q}^{W}{\Psi}^{*}{\Upsilon}^{*}\right)}^{1/2}\right)\left(t\right).$

Applying Lemma (1) again gives

$\Upsilon \left({\displaystyle {\int}_{0}^{t}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\Psi \left(s\right)\text{d}L\left(s\right)\right)={\displaystyle {\int}_{0}^{t}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\sigma \left(s\right)\text{d}N\left(s\right)$ (13)

where $N\left(t\right)=B\left(\Theta \left(t\right)\right)$ , $t\in {\mathbb{R}}_{+}$ and

$\sigma \left(t\right)=\left(\Upsilon \Psi \left(t\right){Q}^{W}\Psi {\left(t\right)}^{\mathrm{*}}{\Upsilon}^{\mathrm{*}}\right)\mathrm{.}$

here, N is an n-dimensional Lévy process with the desired properties and the martingale covariance ${Q}^{L}$ of L and ${Q}^{W}$ of W coincide based on the above theorem ( [10]).

Based on results in [3] we conclude the construction of the stochastic integral based on Equation (10);

$X\left(t\right)={\displaystyle {\int}_{0}^{t}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\Gamma \left(s\right)\text{d}L(\; s\; )$

where $L\in {M}^{2}\left(U\right)$ and $\Gamma \left(s\right)=\Upsilon \circ \Psi \left(s\right)$ are operators from U to another Hilbert space H.

However, in order to express the subordination in terms of standard Brownian motion, we make the following definition:

Definition 2. *A subordinated Brownian moion L with values in some Hilbert space U is a Lévy process such that there is a U-valued Brownian motion B and a subordinator
$\Theta $ which is independent of B such that *

$L\left(t\right)=B\left(\Theta \left(t\right)\right)\mathrm{.}$

Subordinated Brownian motion L, N are of same type if there are Brownian motions ${B}_{L}$ , ${W}_{N}$ and subordinators ${\Theta}_{L}$ , ${\Theta}_{N}$ such that ${B}_{N}$ , ${\Theta}_{N}$ are independent, ${B}_{L}$ , ${\Theta}_{L}$ are independent, ${\Theta}_{L}$ , ${\Theta}_{N}$ have the same law and for any $t\in {\mathbb{R}}_{+}$ . Subordinated Brownian motions have some similarities with Brownian motion if the subordinator has finite first moment. In particular, the set of integrands can be compared easily. First, we recall the notion of time-changed filtrations:

Definition 3. *A time-change is a right-continuous increasing family
${\left(\Theta \left(t\right)\right)}_{t\in {\mathbb{R}}_{+}}$ of stopping times with respect to some right-continuous filtration
${\left({F}_{t}\right)}_{t\in {\mathbb{R}}_{+}}$ *.* The time changed filtration is the filtration given by *

${F}_{t}^{\Theta}={F}_{\Theta \left(t\right)}={\displaystyle \underset{s>t}{\cap}{F}_{\Theta \left(s\right)}t}\in {\mathbb{R}}_{+}\mathrm{.}$

Let X be an F-adapted stochastic process and assume that the time-change is finite valued. The time-changed process by ${X}^{\Theta}\left(t\right)\mathrm{=}X\left(\Theta \left(t\right)\right)$ , $t\in {\mathbb{R}}_{+}$ . Therefore, adding the following results on stochastic integration with respect to subordinated Brownian motions.

Lemma 1. *Let B be a mean-zero Brownian motion with values in a separable Hilbert space U relative to some filtration
${\left({F}_{t}\right)}_{t\in {\mathbb{R}}_{+}}$ *,* H be another separable Hilbert space*,*
$\psi \in {\aleph}_{L}^{2}\left(G\mathrm{,}H\right)$ *,*
$\Theta $ be a non-zero subordinator with finite moment such that
$\Theta \left(t\right)$ is a stopping time for each
$t\in {\mathbb{R}}_{+}$ *.* Let
${\left({G}_{t}\right)}_{t\in {\mathbb{R}}_{+}}$ be the time-changed filtration given by
${G}_{t}={F}_{\Theta \left(t\right)}$ and
$L\left(t\right)=B\left(\Theta \left(t\right)\right)$ *,*
$t\in {\mathbb{R}}_{+}$ *.* Then L is a U-valued square integrable Lévy process and there is an isometric embedding
$\Gamma \mathrm{:}{\aleph}_{L}^{2}\left(G\mathrm{,}H\right)\to {\aleph}_{B}^{2}\left(F\mathrm{,}H\right)$ such that *

${\int}_{0}^{t}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\psi \left(s\right)\text{d}L\left(s\right)={\displaystyle {\int}_{0}^{\Theta \left(t\right)}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\Gamma \left(\psi \right)\left(s\right)\text{d}B\left(s\right)\mathrm{,}P\text{-almostsurely$ (14)

where the left stochastic integral is with respect to the filtration ${\left({G}_{t}\right)}_{t\in {\mathbb{R}}_{+}}$ , the right stochastic integral is with respect to the filtration ${\left({F}_{t}\right)}_{t\in {\mathbb{R}}_{+}}$ .

Proof: See proof in [7].

Therefore, we can see from the above construction that Equation (14) holds for elementary integrands, and thus for all integrands by a density argument in the Hilbert space. We can use results from Lemma (1) to derive a link between functionals of the infinite-dimensional stochastic integral and finite dimensional versions of it.

2.3. Evaluating Energy Forward Dynamics as a Subordinated Hilbert Space

In this section, we use the properties of Hilbert space to represent the forward and futures prices in energy markets as an element of Hilbert space of functions. Motivated by results in [10], it is observed that the various relevant forwards and futures contracts traded in energy markets, which deliver the underlying over a period rather than at a fixed time in the future, can be understood as a bounded operator on a suitable Hilbert space.

We begin by stating the following relevant assumption from the Filipovic space [11] which supports the Hilbert space appropriate for our considerations.

H1: The functions $h\in H$ are continuous and the pointwise evaluation ${\mathcal{J}}_{\tau}\left(h\right)=h\left(\tau \right)$ is a continuous linear functional on H, for all $\tau \in {R}_{+}$ .

H2: The semigroup $\left\{S\left(t\right)\mathrm{|}t\in {\mathbb{R}}_{+}\right\}$ is strongly continuous in H with infinitesimal generator denoted by A.

Therefore, we use the previously defined class of Hilbert space which follows from the Filipovic space satisfying (H1)-(H2) and which are coherent with economical reasoning about the forward curve $\tau \leftrightarrow h\left(\tau \right)$ .

Since in practice the forward curve is obtained by smoothing data points using smooth fitting methods it is reasonable to assume;

${\int}_{0}^{\infty}}\text{\hspace{0.05em}}{h}^{\prime}{\left(\tau \right)}^{2}\text{d}\tau <\infty $ (15)

defined on ${H}_{w}$ of absolutely continuous functions $h\mathrm{:}{\mathbb{R}}_{+}\to \mathbb{R}$ . Moreover, the curve flattens for large time to maturity $\tau $ . There is no reason to believe that the forward rate for an instantaneous loan that begins in 10 years differs much from one which one day later. We take this into account by penalizing irregularities of $h\left(\tau \right)$ for large $\tau $ by some increasing weighing function $w\left(\tau \right)\ge 1$ , that ${\int}_{0}^{\infty}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}w\left(\tau \right){h}^{\prime}{\left(\tau \right)}^{2}\text{d}\tau $ for a given continuous and increasing weight function $w\mathrm{:}{\mathbb{R}}_{+}\to \left[\mathrm{1,}\infty \right)$ with $w\left(0\right)=1$ . The norm of ${H}_{w}$ is ${\Vert g\Vert}_{w}^{2}=\langle g,g\rangle $ for the inner product

$\langle f,h\rangle =f\left(0\right)h\left(0\right)+{\displaystyle {\int}_{0}^{\infty}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}w\left(\tau \right){h}^{\prime}\left(\tau \right){f}^{\prime}\left(\tau \right)\text{d}\tau .$

here $f\mathrm{,}h\in {H}_{w}$ . We assume that ${\int}_{0}^{\infty}}\text{\hspace{0.05em}}{w}^{-1}\left(\tau \right)\text{d}\tau <\infty $ . The typical choice of weight function is that of an exponential function $w\left(\tau \right)\mathrm{=}{\text{e}}^{\left(\stackrel{\u02dc}{w}\tau \right)}$ for a constant $w>0$ , in which case the integrability condition on the inverse of w is trivially satisfied. From [11], we know that ${H}_{w}$ is a separable Hilbert space. As we shall see, one can realize energy forward prices as linear operators on and in fact interpret energy forward prices as stochastic processes with values in this space.

A simple motivation for the appropriate choice of is the classical model for the dynamics of energy spot prices in the so-called Schwartz dynamics [12].

Here, the spot price at time is given by

for a standardized subordinated return process defined in an Ornstein-Uhlenbeck (OU) process

(16)

with a corresponding risk-neutral process given by

(17)

driven by a Lévy process. We assume that has exponential moments, , are constants, and. It is simple to see that is continuously differentiable for every t, and

From [12], we get the forward price at time, for a contract delivering at time, is given;

with being the logarithm of the moment generating function of. Suppose we model the spot price directly under the pricing measure Q. Letting, we find (by slightly abusing the notation)

Assuming that the weight function w is such that

Then it follows that from the Cauchy-Schwartz inequality and the assumption. Hence, f is uniformly bounded in since

But then

which shows that. If L is a driftless Lévy process, the exponential moment condition on yields that has the representation

(18)

for a constant and Lévy measure. But by the monotone convergence theorem and L’Hopital’s rule, we find that

such that Equation (18) becomes

and therefore when is small. Thus, a sufficient condition for is.

Therefore, we now represent the realization in of general energy forward and futures contracts with a delivery period. Suppose that is the forward price at time t of a contract on energy delivering over the time interval, where [12]. Then we can express this price as

(19)

where is the forward price for a contract “delivering energy” at the fixed time T, and is a deterministic weight function defined by;

(20)

for the forward-style contracts and

(21)

for the futures style.

Here is the risk-free interest rate which we suppose to be constant. In the energy market on NYMEX, say, WTI oil is delivered physically at a location over a given delivery period like month or quarter. We will therefore have the same expression (19) for the oil forward prices as in the case of energy forwards. Meanwhile, futures on WTI oil index deliver the money-equivalent from the aggregate index value over a specified period. Hence, the futures can be expressed as

where is the futures price of a contract that delivers the corresponding index at the fixed delivery time, i.e., the weight of the futures contracts can be expressed by (19) with

(22)

as the weight function. Aiming at a so-called Museila representation of in (19). Define as the time until start of delivery of the swap and as the length of delivery of the swap. With the notation, we write

(23)

for the weight function defined by

(24)

where, and. Referring to the different cases of the weight function, we find that for a forward contract and for the forward-style oil swap. Slightly more interesting are the future-style energy swaps, yielding

(25)

here, we use the future-style contract as defined in Equation (21). Note that all these cases result in a weight function which is independent of time. Furthermore, the only case that depend on and y is given in Equation (25), which becomes in fact stationary in the sense that depends on.

However, when the asset price is a spot commodity, the solution for forwards price (under deterministic interest rates, futures prices equal forward prices) is given by;

(26)

where is the time to expiration of the contract, and with and hence

where is defined by the geometric Lévy process given by;

and

For a Poisson process, measures the sensitivity of the process and at

and we have

Suppose then which violates our sensitivity condition of the process, hence we select, then Equation (26) becomes

(27)

Therefore, suppose that the rate of convenience yield follows the mean-reverting process then recall from Equation (17) that SDE;

(28)

where is the speed of adjustment, is the long-run mean yield, and is the increment to a standard Lévy process. However, from [13], the, where denotes the correlation coefficient. The variance of the change in the net marginal rate of convenience yield (volatility) is representing the measure of jump sizes in the forward driven incomplete market.

Hence, we adopt [14] approach and construct a no-arbitrage portfolio that includes two futures contracts of different maturities and the spot commodity and this approach leads to the following differential equation for futures price;

(29)

with boundary condition

(30)

Since convenience yield is non-traded, the differential Equation in (29) depends on investor risk preferences embedded in the market price of risk for convenience yield,. To obtain a solution, a version of the Feynman-Kac Theorem is involved and, for tractability, the market price of risk is assumed to be constant. This is equivalent to assuming in a general equilibrium framework that the representative investor has a logarithmic utility function. In this special case, the marginal utility of wealth is independent of wealth; and the market price of risk, which is given by the covariance of the change in the convenience yield with the rate of change in the marginal utility wealth, is constant.

Therefore, the Feynman-Kac solution for futures (forwards) price is given by;

(31)

where the expectation is taken with respect to the risk-neutral processes from Equation (17);

(32)

and

(33)

Then, we define the two-factor theoretical futures and forward prices depend on the current level of the spot price, the current level of the convenience yield, , time to maturity, the parameters of the joint process, and the price of a zero-coupon bond with maturity at time T, based on similar approach by [13] such that suppose;

(34)

The standard Jump-diffusion for follows the transformation (34), Ito’s lemma, and the risk-neutral diffusion for spot prices, and

(35)

Integrating (35) on both sides with respect to;

however, from Equation (34) we take, then;

(36)

since the effect of the interest rate in the drift of the spot process cancels out the effect of the interest rate discount factor. The forward solution is therefore given by

(37)

The solution for the forward price in Equation (37) is independent of the assumption of stochastic interest rates as long as is independent of the spot interest rate, and matches the market price of the zero-coupon bond with maturity at T in a no-arbitrage interest model.

For the diffusion assumed for spot price, is independent of S and r, the distribution of is normal, and is the moment-generating function of a normal distribution. In this case, the solution for the forward price is

(38)

The expected value of is

(39)

To solve the integral in (39), the solution for the risk-neutralized stochastic differential equation for convenience yield is used which is given by

(40)

Substituting the expected value of Equation (40) into Equation (39) and integrating with respect to gives

(41)

where

and

The variance of, Equation (36), is

(42)

The variance and the covariance terms in Equation (42) are obtained by replacing with the random part of the process (40), giving

(43)

Solving for the variance and covariance in Equation (43) gives

(44)

The final result for is obtained by using, Equation (41), and, Equation (44) in Equation (38) the forward solution is given by;

(45)

where

which gives the pricing framework for estimation of the energy forward dynamics modeled as a subordinated Hilbert space of linear functional.

3. Main Results

3.1. Exploring the Forward Curve of Energy Commodity

In this section, we select the West Texas Intermediate (WTI) Oil data for analysis and to evaluate the variation in forward rates in a given time interval. The data contains the daily prices of front month WTI oil price traded by NYMEX (New York Mercantile Exchange).

The front month WTI oil price is a futures contract with the shortest duration that could be purchased in the NYMEX market.

The WTI daily oil price is normalized in Figure 1 in other to refine the time series properties of the WTI oil data and the weight functions of the forward-style and the futures-style contract to enhance better pricing model for estimation of the forward prices. Between the first 50 monthly price observation, it is clear that noise exist in the data.

However, the spikes level is benchmarked using the properties Hilbert space linear functional of jump-reduction to reduce the price movement because the stochastic price return contains high volatility with massive jumps.

This is obvious in the daily change prices for the underlying asset as seen in Figure 2 and risk averse investors tend to drive the market with demand, the daily price change seems to favour the market share capitalization for a long period as seen in the forward curve.

Figure 1. WTI crude oil price return with no-jump model.

Figure 2. Forward rates forward curve of WTI crude oil price return with jump-diffusion.

Furthermore, it follows from Equation (23) and as seen in Figure 2 that the dynamics of forward price curves can be represent as Hilbert space functional as a result of this assertion we accept that if is an -valued stochastic process, then will be a stochastic process with values in as well.

As the subordinated process, the time change stochastic process as shown in the forward curve in Figure 2 defines how the incomplete market trading options of WTI oil futures and forwards displays the highs and lows of the market situation.

Investors in the energy market are either willing to be risk neutral and trade for high profits or not but the energy forward pricing framework clearly favors the idea that with changing and the forward rates is bound to increase marginal yield over time.

3.2. Estimation of Forward Prices with and without Lévy Jump

In order to avoid false detection of jumps in the forward curves pricing framework, we employed the following model for simulation;

(46)

where is volatility, the drift component and is a Wiener proces that follows a Poisson distribution. In other to estimate the simulation model for a case of constant volatility verses stochastic volatility without Lévy jumps, we integrate both sides of Equation (46);

(47)

With a constant volatility and the standard Brownian motion following a Poisson distribution such that, then the simulated stock prices of the stochastic model is presented in Figure 3 with 1000 observations.

Figure 3 clearly indicates that the fluctuation in the stock prices without the presence of Lévy jumps is benchmarked by the volatility value of upwards. However, we notice one spike in price of the stock which is absolute with a 0.05 level of significance and assigned to the variation in price dynamics of the model.

Figure 3. Model based stock prices without Lévy jumps.

With a stochastic volatility and the standard Brownian motion following a Poisson distribution such that, then the simulated stock prices of the stochastic model is presented in Figure 4 with 1000 observations.

With the big-jump and small-jump rules, Figure 5 clearly shows the detection and arrival of jumps from the small-jumps to big-jumps as the frequency of the observation increases with an initial volatility, given a sensitivity level. Although there exist a drift part, the number of jumps with constant volatility in all 300 observations is seen to be evenly distributed by both rules as can be seen in Table 1.

The Log prices as described by Figure 5 and Figure 6 respectively shows that although there appears to be a varying changes in forward rates, with a slight change in volatility, the magnitude of price fluctuation in a stochastic volatility situation is the same with that of a constant volatility which further reinforces the impact of the presence of closed-gap normed spaces in the sequence of forward rates as introduce by the properties of the Banach space (special type of Hilbert space) (Table 2).

The forward curve without Lévy jumps as described in Figure 7 is estimated with the initial Lévy process, the average stock price for WTI, and is the drift component for which the price fluctuate with sensitivity and constant volatility of, and respectively where the jumps, are evaluated by the small-jump and big-jump rules accordingly with a convenient yield and correlation coefficient of the forward rates as and respectively.

4. Conclusion and Suggestions

Representing energy forward and futures contracts as a time-changing stochastic process in a Hilbert-space of functions shows clearly, that an arbitrage-free forward price can be derived from the buy-and hold strategy in the energy market thereby enabling investors in the market willing to be salvage from the market uncertainties as well as Arrow-Debreu situations to execute a spot or forward contracts depending on the time and place the market becomes favorable.

Figure 4. Model based stock prices without Lévy jumps.

Figure 5. Simulated constant price volatility without jumps.

Figure 6. Simulated stochastic price volatility without Jumps.

Figure 7. Estimated Forward curve for stock prices without Lévy jumps.

Table 1. Parameters of estimation for jump sizes with constant volatility.

Table 2. Parameters of estimation for jump sizes with stochastic volatility.

However, subsequent evaluation can take into consideration the impact of regime-shifts present in the time-changes as well as measure the sizes of such shifts and their corresponding impact to the price volatility.

5. Primary Contribution

· The multivariate subordinated processes are deduced as a moment generating function in order to bench mark the time-change process.

· The subordinated process in a continuous time is represented as a stochastic integral as Hilbert space linear functional.

· The energy forward is represented while incorporating the completeness properties of Banach space (a special type of Hilbert space) for an incomplete energy market.

· The pricing framework for energy forward contracts is used to evaluate the daily change price of an underlying energy contracts.

Acknowledgements

Our sincere thanks to the Pan African University Institute for Basic Science and Technology, Kenya for their support, and special thanks to Dr. Jane Aduda Akinyi for moral support guide throughout the development of this paper.

References

[1] Barndorff-Nielsen, O.E., Benth, F.E. and Veraart, A.E. (2013) Modelling Energy Spot Prices by Volatility Modulated Lévy-Driven Volterra Processes. Bernoulli, 19, 803-845.

https://doi.org/10.3150/12-BEJ476

[2] Benth, F.E., Benth, J.S. and Koekebakker, S. (2008) Stochastic Modelling of Electricity and Related Markets. World Scientific, 11, 352. https://doi.org/10.1142/6811

[3] Benth, F.E. and Kruhner, P. (2015) Derivatives Pricing in Energy Markets: an Infinite-Dimensional Approach. SIAM Journal on Financial Mathematics, 6, 825-869.

https://doi.org/10.1137/15100268X

[4] Benth, F.E. and Zdanowicz, H. (2016) Pricing and Hedging of Energy Spread Options and Volatility Modulated Volterra Processes. International Journal of Theoretical and Applied Finance, 19, No. 1. https://doi.org/10.1142/S0219024916500023

[5] Benth, F.S. and Benth, J.S. (2012) Modeling and Pricing in Financial Markets for Weather Derivatives. Technical Report, World Scientific, Singapore.

https://doi.org/10.1142/8457

[6] Benth, F.E. and Krühner, P. (2012) Subordination of Hilbert Space Valued Lévy Processes. arXiv preprint arXiv:1211.6266.

[7] Clewlow, L. and Strickland, C. (2000) Energy Derivatives: Pricing and Risk Management. Lacima Publications, London.

[8] Filipovic, D. (2001) Consistency Problems for Heath-Jarrow-Morton Interest Rate Models. Springer Science & Business Media. https://doi.org/10.1007/b76888

[9] Hilliard, J.E. and Reis, J. (1998) Valuation of Commodity Futures and Options under Stochastic Convenience Yields, Interest Rates, and Jump Diffusions in the Spot. Journal of Financial and Quantitative Analysis, 33, 61-86.

https://doi.org/10.2307/2331378

[10] Hull, J. and White, A. (1987) The Pricing of Options on Assets with Stochastic Volatilities. The Journal of Finance, 42, 281-300.

https://doi.org/10.1111/j.1540-6261.1987.tb02568.x

[11] Kallsen, J. and Kruhner, P. (2014) On a Heath-Jarrow-Morton Approach for Stock Options. Finance and Statistics, 19, 583-615.

https://doi.org/10.1007/s00780-015-0263-1

[12] Kilian, L. and Murphy, D. (2013) The Role of Inventories and Speculative Trading in the Global Market for Crude Oil. Journal of Applied Econometrics, 29, 454-478.

https://doi.org/10.1002/jae.2322

[13] Monroe, I. (1972) On Embedding Right Continuous Martingales in Brownian Motion. Annals of Mathematical Statistics, 43, 1293-1311.

[14] Peszat, S. and Zabczyk, J. (2007) Stochastic Partial Differential Equations with Lévy Noise. Cambridge University Press, Cambridge.

https://doi.org/10.1017/CBO9780511721373