<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">AM</journal-id><journal-title-group><journal-title>Applied Mathematics</journal-title></journal-title-group><issn pub-type="epub">2152-7385</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/am.2012.37113</article-id><article-id pub-id-type="publisher-id">AM-19995</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Energy Portfolio Management with Entry Decisions over an Infinite Horizon
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>hen</surname><given-names>Liu</given-names></name><xref ref-type="aff" rid="aff1"><sub>1</sub></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><label>1</label><addr-line>Department of Engineering Management &amp;amp; System Engineering Missouri University of Science &amp;amp; Technology Rolla, MO 65409, USA</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>zliu@mst.edu</email></corresp></author-notes><pub-date pub-type="epub"><day>30</day><month>07</month><year>2012</year></pub-date><volume>03</volume><issue>07</issue><fpage>760</fpage><lpage>764</lpage><history><date date-type="received"><day>May</day>	<month>5,</month>	<year>2012</year></date><date date-type="rev-recd"><day>June</day>	<month>5,</month>	<year>2012</year>	</date><date date-type="accepted"><day>June</day>	<month>13,</month>	<year>2012</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  We study a firm that has a conventional plant and considers introducing a new plant as an alternative to generate electricity. The firm’s decision includes the optimal entry time for the new plant, and the optimal dispatch between the existing plant and the new plant after it has been constructed to maximize the expected profit over an infinite time horizon. Under geometric Brownian motion, we formulate the problems as non-regular mixed optimal stopping/control problem. Due to the intractability of the mixed problem, we decompose it into two auxiliary problems, and characterize the optimal strategies in closed-form by standard value-matching and smooth-pasting conditions. Our numerical example confirms our theoretical results.
 
</p></abstract><kwd-group><kwd>Energy Portfolio Management; Geometric Brownian Motion; Optimal Stopping</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Climate change is now recognized as the major environmental problem facing the world. The factor of most concern that causes climate change is the increase in carbon dioxide levels due to emissions from fossil fuel combustion. Therefore, construction of an alternative power plant is crucial to reducing carbon dioxide emission. By introducing an alternative plant, the firm therefore has an energy portfolio. A decision on when a new plant should be built must consider this uncertainty due to fluctuations in electricity prices. We study this problem with the optimal entry decision for the new plant, given fixed capital investment and the optimal dispatch decision for the conventional plant and the alternative plant with the objective of maximizing the long-term expected profit under geometric Brownian motion for the electricity prices over an infinite horizon.</p><p>Investments in and operations of power plants have been widely studied. Deng et al. [<xref ref-type="bibr" rid="scirp.19995-ref1">1</xref>] valued electricity derivatives by futures-based replication due to the nonstorable nature of electricity. Tseng and Barz [<xref ref-type="bibr" rid="scirp.19995-ref2">2</xref>] evaluated a power plant in the short-term with unit commitment constraints using a real-options approach. By the same approach, Tseng and Lin [<xref ref-type="bibr" rid="scirp.19995-ref3">3</xref>] evaluated a power plant involving processes of electricity and fuel prices. Thompson et al. [<xref ref-type="bibr" rid="scirp.19995-ref4">4</xref>] studied the valuation and optimal operations of hydroelectric and thermal power generators through an optimal control and partial-integral-differential-equations (PIDEs) approach. Takashima et al. [<xref ref-type="bibr" rid="scirp.19995-ref5">5</xref>] analyzed the optimal entry strategy of two firms under price uncertainty and competition with real options and game theory. Recently, Deng et al. [<xref ref-type="bibr" rid="scirp.19995-ref6">6</xref>] studied the same problem over a finite time horizon, and solved the resulting partial differential equation (PDE) by finite difference method. Liu [<xref ref-type="bibr" rid="scirp.19995-ref7">7</xref>] studied the optimal time to abandon a plant of a firm with a portfolio of two plants over an infinite horizon by stochastic control approach, and the optimal policy is obtained in closed-form.</p><p>This paper studies the timing that a firm invests in a new technology. We assume the firm owns a plant, and considers adding a new plant while maximizing the expected long-term profit. Since the firm can generate electricity by a portfolio of two plants, the optimal dispatch of these two plants needs to be determined after the new plant is constructed. (Admittedly, it is possible to include more plants with various generating methods; we discuss this extension in Section 5. Under the geometric Brownian motion of long-term electricity prices [<xref ref-type="bibr" rid="scirp.19995-ref8">8</xref>], we formulate the decision problem as a mixed stochastic control problem. Due to the intractability of the mixed problem, we decompose it into two auxiliary problems: one is a regular stochastic control problem, the other one is an optimal stopping problem. The solution to the auxiliary problems is equivalent to the original control problem, and is obtained in closed-form.</p><p>Our contribution is two-fold: First, we formulate the problem as a mixed stochastic control problem. Therefore, the optimal entry decision and optimal dispatch are addressed accordingly. To the best of our knowledge, the mixed control problem cannot be solved directly. Second, we decompose the intractable mixed problem into two auxiliary problems, and the corresponding value functions satisfy a standard Hamilton-Jacobi-Bellman (HJB) equation or variational inequality (VI). We obtain the closed-form solutions to the value functions.</p><p>The rest of the paper is organized as follows. We formulate the decision problem as a non-standard stochastic control problem in Section 2. In Section 3, we write the equivalent form of the value function to the control problem. By standard arguments, we obtain the closedform of the value function. In Section 4, we provide a numerical example to confirm our results and sensitivity analysis. Finally, conclusions and future research directions are presented in Section 5.</p></sec><sec id="s2"><title>2. Problem Formulation</title><p>We introduce the following notation to formulate the problem.</p><p>• <img src="13-7400834\ec57e60d-a99b-4ae2-ade6-9a0c767d7331.jpg" />: Electricity price [$/MWh]</p><p>• <img src="13-7400834\4599aa54-91d7-447c-8b58-ac8d045a5b4d.jpg" />: the risk-adjusted discount rate</p><p>• <img src="13-7400834\49ebd4da-27e6-4200-8a14-866c064b2032.jpg" />: maximum proportion of total wealth invested in alternative method (AL)</p><p>• <img src="13-7400834\34dccad0-f20c-454b-84d1-19951472fbc8.jpg" />: proportion of total wealth invested in AL [Decision variable]</p><p>• <img src="13-7400834\c0dafecf-1ddf-4d0b-8f99-f3165696456d.jpg" />: Production rate of the conventional method (CON)</p><p>• <img src="13-7400834\34c7a426-1341-4e50-be15-cd0d13df32db.jpg" />: Production rate of AL</p><p>• <img src="13-7400834\3dd111b1-29db-45bb-b07e-dd504fd0919b.jpg" />: Total cost of generating <img src="13-7400834\a0673347-5f22-4f75-82ae-3fa46570b0d8.jpg" /> units of electricity from CON</p><p>• <img src="13-7400834\398df7ff-d968-4575-8aaa-d93d9bda9347.jpg" />: Total cost of generating <img src="13-7400834\9c26653c-2346-4088-b185-5683acb48400.jpg" /> units of electricity from AL</p><p>• <img src="13-7400834\746fa403-0d3c-48bf-9eb1-25690d698b44.jpg" />: time to construct AL [Decision variable]</p><p>• <img src="13-7400834\64e9cf9b-2492-4495-b145-d91f72faaf2a.jpg" />: capital investment for constructing AL [\$]</p><p>We assume the long-term electricity price follows the standard geometric Brownian motion [<xref ref-type="bibr" rid="scirp.19995-ref8">8</xref>]:</p><disp-formula id="scirp.19995-formula32409"><label>(1)</label><graphic position="anchor" xlink:href="13-7400834\cf83b04f-43ea-4d1c-bd96-6dae09101f1f.jpg"  xlink:type="simple"/></disp-formula><p>where μ and σ are expected growth rate and volatility of the electricity price respectively, and {Bt} is a Wiener processes. Assume x is the initial position of electricity price. That is, X<sub>0</sub> = x.</p><p>The objective of the firm is to choose an optimal stopping time τ to construct AL, and an optimal proportion α after τ in order to maximize the expected profit. If we define the expected discounted profit functional J given initial electricity price x, proportional investment in AL α, and the time to construct AL τ as</p><disp-formula id="scirp.19995-formula32410"><label>(2)</label><graphic position="anchor" xlink:href="13-7400834\4823fa3a-36dc-4dea-a70d-fe1fc3706cc0.jpg"  xlink:type="simple"/></disp-formula><p>subject to (1), where E<sub>x</sub> is the expectation with respect to x and αt is obviously a function of time t, then the value function u is defined as</p><disp-formula id="scirp.19995-formula32411"><label>(3)</label><graphic position="anchor" xlink:href="13-7400834\24b22ee0-0461-46bb-9433-98199f05a139.jpg"  xlink:type="simple"/></disp-formula><p>where Γ is the set of stopping times.</p></sec><sec id="s3"><title>3. Solution Methods</title><p>In order to solve the non-standard stochastic control problem (3), we proceed as follows: By introducing an auxiliary function x, we obtain the equivalent function w to the value function u, which solves an optimal stopping problem. We then get the closed-form solution of w by value-matching and smooth-pasting conditions.</p><sec id="s3_1"><title>3.1. Equivalent Problem to (3)</title><p>We define the auxiliary function v as the expected profit from the portfolio assuming AL is constructed:</p><disp-formula id="scirp.19995-formula32412"><label>(4)</label><graphic position="anchor" xlink:href="13-7400834\ff2aad63-0aaf-495e-b17e-7216eef686c7.jpg"  xlink:type="simple"/></disp-formula><p>It is easy to get</p><disp-formula id="scirp.19995-formula32413"><label>(5)</label><graphic position="anchor" xlink:href="13-7400834\a3c694d0-f1e8-4e23-a180-bd58b564b340.jpg"  xlink:type="simple"/></disp-formula><p>by the fact that<img src="13-7400834\d780477e-f460-444b-bf3b-4dea6046ddae.jpg" />, where</p><disp-formula id="scirp.19995-formula32414"><label>(6)</label><graphic position="anchor" xlink:href="13-7400834\5f2d638e-fc0d-4c6d-9f4d-98c015f5bd95.jpg"  xlink:type="simple"/></disp-formula><p>If we define the value function w as</p><disp-formula id="scirp.19995-formula32415"><label>(7)</label><graphic position="anchor" xlink:href="13-7400834\75c75f9b-217b-452b-9375-2b962d4ec501.jpg"  xlink:type="simple"/></disp-formula><p>we can formally prove that w is equivalent to u in (3) (See [<xref ref-type="bibr" rid="scirp.19995-ref9">9</xref>]).</p><p>The value function w satisfies the combination of variational inequality of optimal stopping problem as follows</p><disp-formula id="scirp.19995-formula32416"><label>(8)</label><graphic position="anchor" xlink:href="13-7400834\d29f68fb-5fe9-4e4f-9c3a-624252ef6d27.jpg"  xlink:type="simple"/></disp-formula><p>where the generator L is defined as</p><disp-formula id="scirp.19995-formula32417"><label>(9)</label><graphic position="anchor" xlink:href="13-7400834\e19fc6ae-ef38-46cb-9c66-873d6637878f.jpg"  xlink:type="simple"/></disp-formula></sec><sec id="s3_2"><title>3.2. Closed-Form Solution to w</title><p>First we consider the solution to the following ordinary differential equation (ODE):</p><disp-formula id="scirp.19995-formula32418"><label>(10)</label><graphic position="anchor" xlink:href="13-7400834\1d97d875-daaa-496d-ab66-fbab9c046ffb.jpg"  xlink:type="simple"/></disp-formula><p>A special solution <img src="13-7400834\301794c2-95ca-4852-9ef7-1473f19a1e05.jpg" /> to (10) can be easily identified as</p><disp-formula id="scirp.19995-formula32419"><label>(11)</label><graphic position="anchor" xlink:href="13-7400834\617c0a2f-8b90-4144-98e0-f7fd8a994ad1.jpg"  xlink:type="simple"/></disp-formula><p>If we try a function w of the form</p><p><img src="13-7400834\2ea37528-8965-40f1-9128-9d3492ad4b57.jpg" />, for some constant β,&#160;&#160;&#160;&#160; &#160;(12)</p><p>we get</p><disp-formula id="scirp.19995-formula32420"><label>(13)</label><graphic position="anchor" xlink:href="13-7400834\18af774a-78c8-433d-a893-07f5c699ec15.jpg"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.19995-formula32421"><label>(14)</label><graphic position="anchor" xlink:href="13-7400834\490e7980-fd98-4d73-a67a-b57305bdc5ea.jpg"  xlink:type="simple"/></disp-formula><p>Note that</p><disp-formula id="scirp.19995-formula32422"><label>(15)</label><graphic position="anchor" xlink:href="13-7400834\3e7a9346-e62b-41fc-8761-27371aba59ff.jpg"  xlink:type="simple"/></disp-formula><p>Therefore if we assume<img src="13-7400834\a4c5a1a0-e7c5-4d68-a29d-9ff1dd0bebcc.jpg" />, then we can get there exists <img src="13-7400834\4e75489b-598e-4db0-b742-7a99a3929bfd.jpg" /> such that <img src="13-7400834\78a86313-2143-45c3-8ae0-d2841b79a692.jpg" /></p><p>Next we solve for the explicit form of solution to Problem (7). With the value of<img src="13-7400834\82faba81-51fc-4be8-a805-6735a77e9c72.jpg" />, we put</p><disp-formula id="scirp.19995-formula32423"><label>(16)</label><graphic position="anchor" xlink:href="13-7400834\168731a6-2732-46dd-a899-3e980bec23a9.jpg"  xlink:type="simple"/></disp-formula><p>for constants C and <img src="13-7400834\13ad12ba-99c1-4e92-9d1d-c6823caddbae.jpg" /> to be determined.</p><p>By value matching condition [<xref ref-type="bibr" rid="scirp.19995-ref10">10</xref>] at<img src="13-7400834\192a2b42-eaac-4e77-aaf8-db92cfda095c.jpg" />, we have</p><disp-formula id="scirp.19995-formula32424"><label>(17)</label><graphic position="anchor" xlink:href="13-7400834\fad5d9ac-f525-4a9e-9165-3bf137ae4dbb.jpg"  xlink:type="simple"/></disp-formula><p>and smooth pasting condition [<xref ref-type="bibr" rid="scirp.19995-ref10">10</xref>] at<img src="13-7400834\f5224c23-33c7-4b5b-9699-96dc4075481e.jpg" />, we have</p><disp-formula id="scirp.19995-formula32425"><label>(18)</label><graphic position="anchor" xlink:href="13-7400834\ffae1011-3e59-481f-a07f-fd3c1b467b05.jpg"  xlink:type="simple"/></disp-formula><p>It is easy to see that</p><disp-formula id="scirp.19995-formula32426"><label>(19)</label><graphic position="anchor" xlink:href="13-7400834\aca1bb8c-fadf-48dc-9d8c-9e6dbe99803e.jpg"  xlink:type="simple"/></disp-formula><p>and</p><disp-formula id="scirp.19995-formula32427"><label>(20)</label><graphic position="anchor" xlink:href="13-7400834\d315ad11-654f-44ac-9d37-e3a1bd04d42e.jpg"  xlink:type="simple"/></disp-formula><p>(20) requires that</p><disp-formula id="scirp.19995-formula32428"><label>(21)</label><graphic position="anchor" xlink:href="13-7400834\96817ab6-0f85-48a0-b8be-99a81440ac45.jpg"  xlink:type="simple"/></disp-formula><p>Therefore we obtain</p><disp-formula id="scirp.19995-formula32429"><label>(22)</label><graphic position="anchor" xlink:href="13-7400834\808ea56e-c912-498b-b84f-d8ef7afaf308.jpg"  xlink:type="simple"/></disp-formula><p>(20) also requires that</p><disp-formula id="scirp.19995-formula32430"><label>(23)</label><graphic position="anchor" xlink:href="13-7400834\da84dd4f-85ca-4411-8a14-44f2620429ed.jpg"  xlink:type="simple"/></disp-formula><p>Plugging (22) into (20) yields</p><disp-formula id="scirp.19995-formula32431"><label>(24)</label><graphic position="anchor" xlink:href="13-7400834\5ebbe3a3-1e78-413e-8cc7-c392af2ebf69.jpg"  xlink:type="simple"/></disp-formula><p>In summary, the value function u is equivalent to w, which has the closed-form (16). The unknown constants are determined by (22) and (24) with constraints (23).</p></sec></sec><sec id="s4"><title>4. Numerical Example and Sensitivity Analysis</title><p>A numerical example with the following parameters</p><p><img src="13-7400834\fc65a0d7-3821-4fe1-925a-182e11d205a2.jpg" /></p><p>The results are as follows</p><p><img src="13-7400834\66a87fec-e0c8-4b01-b729-daa52e87285b.jpg" /></p><p>They indicate that when electricity price is higher than 91.97, it is optimal to construct the alternative plant. Otherwise, the decision-maker needs to keep the conventional method.</p><p>Next we carry about sensitivity analysis by changing one of the given parameters to see how it affects the threshold level <img src="13-7400834\67ed0098-ed1b-47e7-a30c-f9f6eb3020bd.jpg" /> as in (22).</p><p>First we change the capital investment K from $400 to $600 with other parameters fixed. <xref ref-type="fig" rid="fig1">Figure 1</xref> shows the threshold level <img src="13-7400834\b773b439-50bf-4319-9058-e68e3f916706.jpg" /> increases linearly from $81.76 to $102.2 as K increases. This confirms the result of (22) and is consistent with the intuition that higher capital investment discourages the firm from constructing the new plant.</p><p>Second, if the electricity price growth rate μ changes from 1% to 3% with other parameters fixed, then the threshold level <img src="13-7400834\2496442d-1d3a-4650-86fb-f1bc47a5a873.jpg" /> decreases from $84.38 to $77.47 as shown in <xref ref-type="fig" rid="fig2">Figure 2</xref>. This result explains the fact: the higher the electricity price, the earlier a firm tends to construct the new plant. <xref ref-type="fig" rid="fig2">Figure 2</xref> also shows a nonlinear relationship between the growth rate μ and the threshold level <img src="13-7400834\204e99cb-4ed7-4212-845d-0503925e9e46.jpg" /> as in (14) and (22).</p><p>Third, if the electricity price volatility σ changes from 10% to 200% with other parameters fixed, then the threshold level <img src="13-7400834\eba659ef-542f-4eb6-8b29-79a2d300e912.jpg" /> increases drastically from $84.38 to $2372.1 as shown in <xref ref-type="fig" rid="fig3">Figure 3</xref>. This impact of price volatility can be explained analytically from (14) and (22): as volatility σ is the leading term of quadratic Equation (14), it has huge impact on the solution β and the threshold level <img src="13-7400834\9ab15e84-9a65-4938-adcc-8de156ee05a4.jpg" /> through (22).</p><p>We omit the sensitivity analysis with other parameters as it is straightforward to carry out given (22).</p></sec><sec id="s5"><title>5. Conclusions</title><p>We study the optimal entry decision for alternative plant given a fixed capital investment, and the optimal dispatch decision between the conventional plant and the alternative plant. By introducing two auxiliary problems, we solved the mixed stochastic control problem in closed form.</p><p>This paper can be generalized in the following ways. First, the analysis in this paper is just based on one stochastic process (the electricity price); we admit that there are other stochastic processes that can affect the decisions, such as the cost of the carbon dioxide emission. We would have more stochastic processes in addition to electricity prices process, and we need to solve the multidimensional optimal control problem. Second, switching costs will be incurred when we abandon the conventional method, and singular control technique would be employed to study this problem.</p></sec><sec id="s6"><title>REFERENCES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.19995-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">S. J. Deng, B. Johnson and A. Sogomonian, “Exotic Electricity Options and the Valuation of Electricity Generation and Transmission Assets,” Decision Support Systems, Vol. 30, No. 3, 2001, pp. 383-392. 
doi:10.1016/S0167-9236(00)00112-3</mixed-citation></ref><ref id="scirp.19995-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">C. L. Tseng and G. Barz, “Short-Term Generation Asset Valuation: A Real Options Approach,” Operations Research, Vol. 50, No. 2, 2002, pp. 297-310. 
doi:10.1287/opre.50.2.297.429</mixed-citation></ref><ref id="scirp.19995-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">C. L. Tseng and K. Lin, “Framework Using Two-Factor Price Lattices for Generation Asset Valuation,” Operations Research, Vol. 55, No. 2, 2007, pp. 234-251. 
doi:10.1287/opre.1060.0355</mixed-citation></ref><ref id="scirp.19995-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">M. Thompson, M. Davison and H. Rasmussen, “Valuation and Optimal Operation of Electric Power Plants in Competitive Markets,” Operations Research, Vol. 52, No. 4, 2004, pp. 546-562. doi:10.1287/opre.1040.0117</mixed-citation></ref><ref id="scirp.19995-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">R. Takashima, M. Goto, H. Kimura and H. Madarame, “Entry into the Electricity Market: Uncertainty, Competition, and Mothballing Options,” Energy Economics, Vol. 30, No. 4, 2008, pp. 1809-1830.  
doi:10.1016/j.eneco.2007.05.002</mixed-citation></ref><ref id="scirp.19995-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">J. Deng, Z. Liu and S. E. Grasman, “Energy Portfolio Investment with Delayed Entry Decisions,” under revision. </mixed-citation></ref><ref id="scirp.19995-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Z. Liu, “Energy Portfolio Management with Abandonment Option over an Infinite Horizon,” Proceedings of the Production and Operations Management Society (POMS), Chicago, 20-23 April 2012. </mixed-citation></ref><ref id="scirp.19995-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">R. Pindyck, “The Long-Run Evolution of Energy Prices,” The Energy Journal, Vol. 20, No. 2, 1999, pp. 1-27. 
doi:10.5547/ISSN0195-6574-EJ-Vol20-No2-1</mixed-citation></ref><ref id="scirp.19995-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">X. Guo and H. Pham, “Optimal Partially Reversible Investment with Entry Decision and General Production Function,” Stochastic Processes and Their Applications, Vol. 115, No. 5, 2005, pp. 705-736. 
doi:10.1016/j.spa.2004.12.002</mixed-citation></ref><ref id="scirp.19995-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">A. K. Dixit and R. S. Pindyck, “Investment under Uncertainty,” Princeton University Press, Princeton, 1994.</mixed-citation></ref></ref-list></back></article>