## Services on Demand

## Article

## Indicators

- Cited by SciELO
- Access statistics

## Related links

- Cited by Google
- Similars in SciELO
- Similars in Google

## Share

## Pesquisa Operacional

##
*Print version* ISSN 0101-7438

### Pesqui. Oper. vol.32 no.1 Rio de Janeiro Jan./Apr. 2012 Epub Apr 03, 2012

#### http://dx.doi.org/10.1590/S0101-74382012005000009

**Long-term planning of a container terminal under demand uncertainty and economies of scale**

**Antonio G.N. Novaes ^{I,}^{*}; Bernd Scholz-Reiter^{II}; Vanina M. Durski Silva^{I}; Hobed Rosa^{I}**

^{I}Federal University of Santa Catarina (UFSC), Florianópolis, SC. E-mails: novaes@deps.ufsc.br / vaninadurski@gmail.com / obbdy@yahoo.com

^{II}BIBA, University of Bremen, Bremen, Germany. E-mail: bsr@biba.uni-bremen.de

**ABSTRACT**

Governments of many developing countries are presently granting concessions to private organizations to operate container terminals. A key element in this process is the technical and economic evaluation of the project. In practice, the demand levels over time are volatile. The first part of the paper deals with a capacity expansion problem with stochastic demand. Next, an optimization model employs a dynamic programming routine to find the best epochs to install new berths. Ship waiting times are computed with queuing models and are introduced into the main model as upper restrictions based on international practice. Finally, a real options analysis is performed with the objective of evaluating the effects of an eventual abandonment of the project at a pre-defined time horizon. The model is then applied to the container terminal of the Port of Rio Grande, Brazil.

**Keywords:** container terminals, capacity expansion, dynamic programming, abandon option.

**1 INTRODUCTION **

Governments in many developing countries traditionally regarded ports as national strategic assets which should remain under public control. However, in some cases container ports occasionally suffer from several critical problems that contribute to high logistics costs throughout the economy. These problems include equipment obsolescence, inefficiencies in labor development and labor allocation, lack of harbor capacity, and inefficiencies in the port administration model. Waiting times for berthing, on the other hand, are usually too high and yard space is often insufficient. In order to alleviate public investments normally scarce, it became common in developing countries the granting of concessions to private companies to operate port terminals. To give support to this operation, a crescent need for robust computer models intended to analyze port developments and their economic feasibility has been reflected in the literature.

The optimization of logistics activities within maritime container terminals has attracted the attention of researchers and practitioners for more than forty years. Novaes & Frankel (1966) analyzed the generation process of unitized cargo (container, pallet, roll-on, roll-off) in a port and its relation to liner shipping, with the aid of a bulk queuing model. Taborga (1969), in his pioneering work on port planning, applied dynamic programming and queuing models to define the optimal epochs to add new berths to the terminal, and at same time choosing the best port technology with the objective of maximizing seaport returns. The economic appraisal of port investments was analyzed by Goss (1967). Over the past recent years, the typical operating scheme of a container terminal has changed significantly. Larger container ships that call a modern port are often "classified ships" (Daganzo, 1990; Huang *et al.*, 2007), in the sense that an appropriate berth slot is allocated to the vessel prior to its arrival. A non-classified ship, on the other hand, may have to wait until a berth of appropriate length becomes available to accommodate it. This means that the classical queuing models, that assume a *first-in, first-out *servicing discipline, are no long adequate to represent this kind of port operation. New empirical approximate queuing approaches have been proposed to partially bypass this limitation (Huang *et al.*, 2007; Morrison, 2007). In parallel, a number of researchers have relied on simulation techniques in order to plan resource utilization at a container terminal (Legato & Mazza, 2001; Kia *et al.*, 2002; Shabayek & Yeung, 2002), thus bypassing the queuing formulation drawback. But one serious limitation of the simulation approach in the context of developing countries is the lack of extensive and detailed data to be used in the models.

The use of containers for intercontinental maritime transport has dramatically increased around the world in the last decade (Notteboom, 2004). From about 150 million twenty-feet equivalent units (TEU) in 1996, the world container turnover has reached 496 million TEU in 2008. A further increase is expected in the upcoming years. Although the flow of containers in Brazil is approximately 1% of the world movement, the participation of containers in the country's foreign trade has increased more than 11% a year in the period 1996-2008 (The World Bank, 2008). Brazil is following a trajectory of economic growth. Part of this growth is fueled by international trade, which is expected to continue to increase in line with Brazilian economic expansion. In such a scenario, port container terminals play an important role and represent a key asset of Brazil's logistics system. In fact, the country has one of the longest coastlines in the world and the existence of container terminals in points of the Atlantic coast is an important advantage in international trade. However, in spite of this location advantage, the Brazilian port system suffers from some critical problems that contribute to high logistics costs throughout the economy. Among other limitations, the lack of harbor capacity, excessive ship waiting times, and inefficiencies in port administration are worth mentioning. Thus, improvements in port efficiency would have a major impact on reducing overall logistic costs.

These factors, among others, indicate that the use of simple, yet robust operations research models in container port development analysis, in the context of developing countries, could help in the economic evaluation of terminal concession initiatives. But this is not an easy task. Forecasting future scenarios depends on a number of intervening factors, such as technology evolution, liner shipping strategies and operating schemes, regional port competition, and government policies (Daschkovska & Scholz-Reiter, 2009; Günther & Kim, 2005; Kim & Günther, 2007). As a consequence, the forecasted demand levels tend to show a relatively high degree of uncertainty. The lack of reliable operational data, on the other hand, precludes in many cases the use of more detailed methods of analysis such as simulation techniques, for example.

The objective of this paper is to present a robust, yet simple approach, to develop a long-term container terminal planning and evaluation model intended to analyze a port concession project in a developing country like Brazil. One characteristic of the methodology is the use of the basic data available in most developing countries. This paper is the improvement of the contents of two previous papers presented at academic conferences, namely Novaes *et al. *(2010a, 2010b). Although the analytical techniques used in this work are not new (capacity expansion modelling, queuing models to estimate waiting times, dynamic programming, real options), the contribution of the paper is essentially an integration of the techniques in one unique model framework, with the objective to solve practical problems of container port planning in developing countries, where the availability of data is often limited and the economic background is usually volatile.

The paper is structured as follows: a mathematical formulation which allows for the introduction of demand uncertainty in the capacity expansion modelling is described in Section 2. A critical review of queuing model developments intended to estimate ship waiting times are analyzed in Section 3. The basic real option concepts and the application of the Black-Scholes equation (Smit & Trigeorgis, 2004) to an abandonment option are described in Section 4. The *Tecon *container terminal of the Port of Rio Grande, Brazil, object of our application, is analyzed next in Section 5. The structure of the capacity expansion model is discussed in Section 6, including the description of a dynamic programming model used to minimize the net present value of berth investments, and its application to the case under analysis. In Section 7, the value of an abandonment option to be exercised at a pre-defined horizon during the project lifetime is estimated with the objective of improving its economical feasibility. Finally, conclusions are drawn in Section 8.

**2 CAPACITY EXPANSION WITH DEMAND VOLATILITY **

Capacity expansion is the process of adding facilities over time in order to satisfy rising demand (Manne, 1961; 1967). Capacity expansion decisions generally add up to a massive commitment of capital. The efficient investment of capital depends on making appropriate decisions in expansion undertakings, in such a way the demand remains satisfied over an extended time period, with a minimum discounted lifespan cost. The literature shows, along the last forty years, an evolving sequence of capacity expansion models, starting with a simple deterministic and linear demand formulation (Manne, 1961), and extending it to non-linear, probabilistic demand situations (Manne, 1967; Freidenfelds, 1980, 1981; Higle & Corrado, 1992; Bean *et al.*, 1992; Novaes & Souza, 2005; Novaes *et al.*, 2010a).

**2.1 The cost function **

The basic way to economically minimize a project lifespan cost is to compute its Net Present Value (*NPV *), where investments and costs are discounted using a continuous interest rate *r*.

The project will increase its own value by adopting the investment plan that minimizes *NPV*. The classical *NPV *calculation presupposes that unknown-risk future monetary values are summarized by their expected values. Assuming no backlogs in demand are admitted, the net present value of investments, costs, and salvage value is

where *I*_{0} is the investment at time *t *= 0, *I _{i} *is the value of the money inversion of order

*i*, (1

__<__

*i*

__<__

*n*),

*r*is the annual continuous interest rate,

*t*is the time instant when the investment

_{i}*I*is made,

_{i}*C*(

*t*) is the system operating cost at time

*t, T*is the duration of the project, and is the salvage value of the project at time

*T*.

The cost functions associated with port operations reported in the literature vary slightly among authors (Huang *et al.*, 2007). Usually these equations incorporate ship costs and berth operating costs in different forms. With regard to berth operating cost, we assume that it depends only on the demand level (throughput). In our approach, demand is an external variable, not depending on port servicing prices and service levels. Additionally, no demand backlog is admitted to occur, *i.e. *it is assumed that the port will provide adequate facilities before full congestion may occur. Consequently, the integral in the right hand side of expression (1) is constant in our capacity expansion model, and therefore does not affect the minimization of *NPV *. Thus, the capacity expansion model adopted in this application assumes only port investments and its salvage values, which reduces expression (1) to

Later, when analyzing cash flow and the abandon option, terminal operating costs and servicing prices will be introduced (see Section 7).

**2.2 Economies of scale **

An important characteristic of most capacity problems is the recognition of economies of scale, *i.e. *large installations usually cost less per produced unit than small ones. But, if the demand level is continuously rising on the long run and demand backlogging is not permitted, excess capacity will occur (Fig. 1). Thus, there is a trade-off between scale economies and excess capacity cost, leading to a compromised optimal solution.

The term "learning curve" effect (Couto & Teixeira, 2005) states that the more time a task is performed, the less time it will be necessary to workers do each subsequent iteration. Learning curve theory states that as the quantity of units produced doubles, labor costs will decrease at a predictable rate. The experience curve, on the other hand, is broader in scope, since it encompasses far more than just labor cost. Now, each time cumulative volume doubles, value added costs (including construction, administration, logistics, etc.) fall by a constant and predictable percentage. Mathematically the experience curve is described by a power law function

where *I *^{(1) }is the value of the building cost of the first unit, *I *^{(m) }is the value of the investment on the *m*^{th }unit, and θ is the elasticity of building cost with regard to the number of units built in sequence. A experience curve that depicts a 25% cost reduction for every doubling of the number of built units is called a "75% experience curve" indicating that unit cost drops to 75% of its original level when installing the second item of the series, and so on. Let δ be the experience factor, obtained via (3) as

Applying logarithms to (4) one has

and the total investment related to the installation of *m *units in series is

Let *m _{i} *be the number of berths to be installed simultaneously at time τ

*, and the corresponding investment. To incorporate returns to scale in the analysis, we substitute variable for*

_{i}*I*in equation (1), thus adding the decision variable

_{i}*m*into the model.

_{i}**2.3 Capacity expansion modelling**

**2.3.1 The linear expansion model**

Manne, in his seminal approach to capacity expansion problems (Manne, 1961) departed from a deterministic demand formulation with a constant linear growth rate. Under this approach it is assumed that the facility has an infinite economic life and, whenever demand catches up with the existing capacity, *x *additional unities of same capacity are installed. No backlog in demand is admitted. The objective is to determine the optimal capacity addition *x*, and the consequent time interval Δ*T *between increments. The time origin is the point *t*_{0}, followed by regeneration points at *t*_{0} + Δ*T , t*_{0} + 2Δ*T , t*_{0} + 3Δ*T *, that represent instants at which the previously existing excess capacity has just been wiped out. Note that when the system has reached *t*_{0} + ΔT , the future looks identical with the way it appeared at time *t*_{0}. The investment necessary to add a capacity increment of size *x*, is assumed to be given by a relationship in the form of a power function, with economies of scale

The demand, in this case, varies linearly over time and is expressed *d _{t} *= µ

*t*, with µ > 0. The objective function to be minimized is of type (2). Let

*NPV*(

*x*) be a continuous function of

*x*that represents the sum of all discounted future investments looking forward from a point of regeneration

*t*. We may write down the following recursive equation (Manne, 1961)

The first term on the right hand side of (8) represents the investment cost incurred at the beginning of the current cycle. Since it is assumed that the facilities have an infinite economic life, the salvage value of the project in this case is nil. The second term represents the sum of all installation costs incurred in subsequent cycles, discounted from the next point of regeneration back to instant *t*_{0}, with a continuous interest rate *r*. It follows directly that

But, since *t = x*/µ, expression (9) becomes

The optimal expansion size is obtained by minimizing *NPV *(*x*) with respect to *x*, and the optimal expansion time interval Δ*t *is given by Δ*t *= /µ (Manne, 1961, 1967). Manne (1961) extended his deterministic model to a stochastic demand formulation, with an underlying linear growing rate µ and standard variation σ (Hull, 1997; Novaes & Souza, 2005). Such a model is based on Brownian motion theory, specifically the Bachelier-Wiener diffusion process in continuous time. One is interested in determining the instant *t *where demand first exceeds the installed capacity. At that time additional facilities will be supplied. In stochastic process terminology, one seeks the first passage time at point *x*, with probability *u*(*t, x*). The following difference equation applies to the problem (Manne, 1961; Novaes & Souza, 2005)

where *p *is the probability that one incremental step Δ*x *is directed to increasing *x*, and *q *= 1 – *p *is the probability that the step Δ*x *is directed to the opposite direction. Expanding expression (11) according to Taylor's theorem up to terms of second order, and dividing bothsides by Δ*t*, one gets

The Laplace transform of *u*(*t, x*) is

where *r*, in this formulation, is both the transformer parameter and the discount interest rate (Manne, 1961). Applying the Laplace transform to each side of (12), and recalling the boundary condition *u*(0, *x*) = 0, one gets a second-order linear differential equation with respect to *x*. Its characteristic equation has two real roots

After some transformations and simplifications (Manne, 1961), one gets

Note that λ_{2} is a function of µ, σ , and *r*, being independent of the quantity *x*. Since *u*(*t, x*)*dt* represents the probability with which exactly *t *years have elapsed between two successive points of regeneration, the probability analogue of (8) is (Manne, 1961)

Just as in the deterministic case, the first term on the right hand side of (16) is the present cost of installing a facility of capacity *x*. The second term of (16) is the probability that the next point of regeneration will occur in *t *units of time, discounting the corresponding cost back to the present, and integrating over all *t*. Departing from (16), and recalling the basic definition of Laplace transform (13), one has (Manne, 1961)

Substituting (15) for into (17), one has

From (14) it follows that λ_{2} < 0 for µ > 0, and *r *> 0. Comparing (10) and (18), it is seen that the following equivalence holds (Bean *et al.*, 1992)

meaning the linear stochastic model can be solved through an equivalent deterministic model in which the original discount rate *r *is substituted by the modified rate *r** (Srinivasan, 1967; Freidenfelds, 1980; Bean *et al.*, 1992; Higle & Corrado, 1992). Manne (1961) indicated that, in the deterministic model, the optimal expansion size increases as the discount rate *r *is lowered. For the probabilistic model, Manne (1961) showed that the greater the variance σ^{2}, the lower will be the value of *r**, and the higher will be the optimal value .

**2.3.2 Generalization to non-linear demand cases **

Assume now that *d*(*t*), the demand for product or service, is a semi-Markov process {*d*(*t*), *t *__>__ 0} (Bean *et al.*, 1992). One is interested in selecting a sequence of capacity expansions to satisfy the demand *d*(*t*) over an infinite horizon at minimum expected discounted cost using a continuous interest rate *r*. A continuous semi-Markov process {*d*(*t*), *t *__>__ 0} is said to be transformed Brownian motion with underlying rate µ and standard variation σ if there exists a non-negative increasing deterministic transformation *h *such that *d*(*t*) = *h*(*d*(*t*)), where {*d*(*t*), *t *__>__ 0} is (linear) Brownian motion with drift µ > 0 and volatility σ > 0. The function *h *is referred to as the transforming function. Bean *et al. *(1992) demonstrated that if {*h*(*d*(*t*)), *t *__>__ 0} is a transformed Brownian motion as defined above, the interest rate *r** to be used in the equivalent deterministic problem is given by the same expression (19) obtained by Manne (1961) for demand following an ordinary (linear) Brownian motion. In other words, *r** is independent of the transformation *h*. In more general terms, Bean *et al. *(1992) demonstrated that a stochastic capacity expansion problem may be solved via a deterministic problem formulation in which (a) the random expansion epochs τ(*x*) are replaced by their expected values; and (b) the original continuous discount factor *r *is replaced by its equivalent value *r** given by (19). In general terms, capacity expansion models with general cost structures and with demand based on a transformed Brownian motion can be solved with the equivalent deterministic problem and the equivalent discount factor (Bean *et al.*, 1992).

**2.3.3 First passage time**

The Laplace transform of the probability distribution *u*(*t, x*), given by (15), allows to computethe expected value and the variance of the first passage time to the level *x*. One has (Novaes & Souza, 2005)

and

For the general case with *d*(*t*) = *h*(*d*(*t*)) one has (Novaes & Souza, 2005)

**3 QUEUING APPROXIMATIONS FOR CONTAINER PORT PLANNING **

Researchers and practitioners have applied queuing models to port planning for more than forty years (Novaes & Frankel, 1966; Taborga, 1969; Noritake, 1978; Noritake & Kimura, 1983; Huang & Wu, 2005; Huang *et al.*, 2007). The first papers dealt with general cargo operations and no product specialization, with the berthing process represented by a *M*/*M*/*c *queue in Kendall's notation (Page, 1972), and with a *first-come, first-served *discipline. Novaes & Frankel (1966) utilized bulk queues in their research work, but the objective of the analysis was not to compute ship waiting times, but to analyze the cargo generating process at the port. Along the years, the handling of cargo at ports became more and more specialized. As a consequence, the queue models that best fit the ship waiting time process have changed, with the Erlang* *and the *M*/*Ek*/*c *Kendall types being the most commonly used in container terminal planning (Huang *et al.*, 1995; Kozan, 1997; Huang *et al.*, 2007). The Erlang distribution has the following probability density function

where *x *is the random variable under analysis, λ is the mean arrival rate, and *k *is the order of the distribution. The Erlang distribution can be viewed as being composed by *k *phases, each phase formed by a negative exponential distribution with the same average length 1/λ. The parameter *k *is an integer, and for *k *= 1 the distribution repeats the negative exponential distribution. The mean and the standard variation for the Erlang distribution (23) are *E*[*x*]= *k*/λ and σ* _{x} *= /λ, respectively. In general, the Erlangian queue is represented by

*Ek*

_{1}/

*Ek*

_{2}/

*c*, where the interarrival times and service times are distributed according to Erlang distributions of order

*k*

_{1}and

*k*

_{2}, respectively. Other queue models can also be treated mathematically. The reader is referred to Fogliatti & Mattos (2007).

**3.1 Approximate methods**

For more complex applications, involving *GI*/*G*/*c *queues with general independent interarrival times and a general distribution of service time, it is customary to use numerical tables or approximate formulas to estimate mean waiting times. Page (1972) proposed an interpolation formula using the numerical results of the four extreme cases represented by the queues *M*/*M*/*c, D*/*M*/*c, M*/*D*/*c *and *D*/*D*/*c*. The resulting queue length for the *D*/*D*/*c *case is nil (Page, 1972). Expressing the mean waiting time as a fraction of the mean service time, Page's interpolation formula is

where * * is the relative waiting time for a * * queue, and *q*^{(D/M/c)}, *q*^{(M/D/c) }and *q*^{(M/M/c) }are the relative waiting times for the queues *D*/*M*/*c, M*/*D*/*c *and *M*/*M*/*c*, respectively. In (24), v* _{a} *and v

*are coefficients of variation of the inter-arrival time and of the servicetime distribution functions respectively, given by*

_{s}where *t _{a} *and

*t*are the ship inter-arrival time and the service time, respectively. For the Erlang queue

_{s}*E*

_{k}_{1}/

*E*

_{k}_{2}/

*c*, one has v

*= 1*

_{a}^{ }and v

*= 1. Let λ be the mean ship arrival rate at the terminal and µ the mean service rate. Page (1972) presented numerical results in table format, which requires interpolation for intermediate values of ρ, the average berth utilization factor given by ρ = λ/µ*

_{s}*c*, which must be less than the unit to let the queue attain a steady state solution. Furthermore, Page's tables are restricted to

*c*

__<__10.

Another approximation, valid for the queue *M*/*E _{k}*/

*c*, is (Cosmetatos, 1976; Huang

*et al.*, 2007)

where π_{0} is the probability of having no elements in the line waiting, according to a *M*/*M*/*c* queue. Recalling that ρ = (λ/µ*c*), *E*[*t _{s}*] = 1/

*t*and

_{s}*q = W*

_{1}/

*E*[

*t*], substituting these expressions in (26), and making the necessary simplifications, one has the Cosmetatos' formula

_{s}where *q *is the mean waiting time expressed as a fraction of the mean service time. Since most of the main container terminals in the world show a *M*/*Ek*/*c *queuing system, the above formula is extensively used in their planning process. Another approximate formula (Sakasegawa, 1976; Morrison, 2007), which is applied to the general queue *GI*/*G*/*c*, is defined as

A numerical comparison was made among the approximations due to Page (1972), Cosmetatos (1976), and Sakasegawa (1976), for the queue *M*/*E*_{3}/*c*, which is largely used in container terminal planning. For *c *in the range 1 to 8, the results of the latter two methods agreed quite well. On the other hand, Page tables are limited to *c *__<__ 10, which, together with the table format of the results, restrict the applicability of his method. Cosmetatos formula (27) was selected to be usedin the application described in Section 6.

**4 THE REAL OPTIONS APPROACH TO PROJECT EVALUATION **

In basic terms, the classical discounted cash flow technique assumes that the organization will accept an investing alternative if the net present value of the project cash flow, computed with an appropriate discount rate, is positive. The calculation of the of the *NPV *project presupposes that unknown-risk future cash flows be summarized by their expected values. It is given by

where *R*(*t*) is the revenue of the project at time *t *and *E*(*t*) is the sum of all expenditures at time *t*, with the other variables having the same meaning as in (1). The project is accepted as economically feasible when *NPV' *> 0. In (29) ship waiting-time costs are not considered as expenditures since vessels belong to external organizations. Instead, constraints limiting ship waiting times are added to the model. These constraints are based on best maritime practices among leading container terminals of the world, in a benchmarking process.

**4.1 A real options overview**

It is widely recognized that the traditional discounted-cash flow approach to the appraisal of capital investment projects, such as the above *NPV *rule, cannot properly capture management's flexibility to adapt and revise later decisions in response to unexpected market developments (Trigeorgis, 1996). The value of an investment with options, on the other hand, is usually greater than the one obtained with the traditional approach because an option truncates possible losses at some point in time. In order to reduce possible negative effects, some measures are taken beforehand in order to limit losses to pre-defined values (Trigeorgis, 1996; Hull, 1997). Thus, one is able to withstand riskier situations, limiting losses to a bearable, pre-established level. In general, a firm's ability to anticipate and to respond flexibly to such changes reduces its expected downside risk and thus lowers costs, and/or raises its revenues. In fact, flexibility adds value to most financial undertakings (Trigeorgis, 1996).

In general terms, an option is the right, but not an obligation, to take an action in the future. Options are valuable when there is uncertainty. For example, a *call option *contract traded on the financial exchanges gives the buyer the right to buy an asset at a certain date and will be exercised only if the price of the asset on that date exceeds the previously specified price. When the transaction involves the right to buy an underlying asset at a fixed price on a future date, one has a *call option*. Conversely, a *put option *occurs when it involves the right to sell the underlying asset at a fixed price. A European option is exercised at a predefined date; conversely, an American option may be exercised at any date prior to the predetermined deadline (Hull, 1997).

In addition to financial options (stock options, forwards, future contracts, etc.), the term *real option *refers to the extension of financial option theory to real, physical assets (Mun, 2002; Smit & Trigeorgis, 2004). The main categories of real options are growth options, the option to expand, the option to wait, the option to switch inputs, outputs, or processes, the option to reduce scale (downsizing), and the abandonment or exit option (Smit & Trigeorgis, 2004). Abandonment options are important in capital-intensive industries where management want to have the flexibility to capture some resale value for assets if their in-use value to the company falls (Novaes & Souza, 2005). This is the kind of real option that will be applied and discussed in this paper.

**4.2 Valuing an abandonment real option**

Basically, there are three general methods to solve real options problems: the *partial differential equation *approach, the *dynamic programming *approach, and the simulation approach (Mun, 2002). The method employed in this paper, the Black and Scholes equation, belongs to the first group and is based on the Wiener-Ito process (Trigeorgis, 1996; Mun, 2002; Smit & Trigeorgis, 2004). Growth models of money-related variables generally presuppose that deviations are not linear, being expressed as a proportion (a percentage) of the asset value. Let *S *be a money-related variable for which we want to study the variation *dS*/*S*. One writes the corresponding Ito process equation (Hull, 1997)

where µ*S *is the expected instantaneous growth rate, σ^{2 }is the volatility, and *dz *is a Wiener process (Hull, 1997). When (30) is observed, the value *S *of the asset follows a log-normal process (Hull, 1997). Departing from equation (30), and considering a generic function *G *of *S *and *t*, Ito derived a lemma that bears his name. Let *S *be the value of an asset such that the variable *dS*/*S *follows the Ito process. Assume that *f *is the price of an option or another derivative of the asset under consideration, with its value depending on *S *and *t*. Departing from Ito's lemma and making substitutions and simplifications, the Black and Scholes differential equation is obtained (Hull, 1997)

where:

*S *= the present value of the asset;

*f *= the price of the derivative, which is a function of *S, t, r *and σ;

*r _{f} *= the risk-free interest rate.

Assume that demand grows stochastically and exponentially with an underlying growing rate µ and standard variation η. The integration of (30), making η = 0, yields the expected growth function of *S *

where *S*_{0} is the value of *S *at *t *= 0. Considering the process evolution from instant *t *– 1 to instant *t*, one has

leading to

Taking a time series containing *n *values of *u _{t}*, one estimates its mean and standard deviation. The time

*t*is usually measured in days, weeks, or months, so we have to compute the equivalent annual standard deviation in order to evaluate the volatility σ to be used in the calculations (Hull, 1997). It is interesting to observe that the Black-Scholes equation separates the risk-free return on the investment, represented by

*r*, from the risk generated by the market variability, expressed by its volatility σ. For a European put option the boundary condition to solve (31) is (Hull, 1997)

_{f}where *X *is the exercise price and *T _{A} *is the time to expiration. The integration of (31), subject to (35), yields

where *N *{·} is the cumulative unit normal probability distribution, with

The value of a European call option is similarly determined assuming the appropriate boundary condition (Hull, 1997)

The Black-Scholes equation gives a way of computing the economic impact of a real option as, for example, the exiting or abandonment option. One way of obtaining the value of an abandonment option is to make the strike price *X *equal to the salvage value of the abandoning project at time *T _{A}*, and the value of the underlying asset

*S*equal to the present value of the project cash flow. The corresponding

*put*option is calculated with (36) and its value is added to the net present value of the project, improving its economic performance. It must be pointed out that the Black-Scholes equation formulation assumes normal distributions to represent stochastic variables. But random variables associated with demand in real options problems may show high kurtosis levels, contrary to the normal distribution which is not skewed. In the application to be described forth, the available data do not allow for kurtosis calculation, and thus the Black-Scholes equation gives, in fact, an approximate solution.

An alternative model to evaluate an abandonment option was developed by Margrabe (1978). Suppose a firm is evaluating an investment proposal at time zero, covering a time horizon of *T* years. At time *T _{A}*(

*T*<

_{A}*T*), the firm has the option to abandon or retain the project. If the project is abandoned, the cash flow realized by the firm is

*X*. On the other hand, if the project is retained, the expected value to be realized by the firm is

*S*. Admitting that both

*X*and

*S*are lognormally distributed (Schnabel, 1992), Margrabe's exchange option pricing model is (Margrabe, 1978)

where *X *is the discounted value of the project asset when abandonment occurs and *S *is the discounted value of the project's value-in-use (Schnabel, 1992). Here *f *is a call option on the abandonment value, and a put option on the project's value-in-use (Margrabe, 1978; Schnabel, 1992). Margrabe (1978) shows that the option (39) is worth more alive than dead, which implies that its holder will not exercise it until the last possible moment. In other words, *f *is a European type option.

Let σ* _{p} *be the volatility of the abandonment process. Margrabe (1978) decomposes into two components, and , according to the following formula

where σ* _{X} *is the volatility associated with the abandonment of the project, σ

*is the volatility associated with the value-in-use of the project, and ρ is the correlation between*

_{S}*X*and

*S*(Margrabe, 1978; Schnabel, 1992). If σ

*= 0, Margrabe's model reduces itself to the Black and Scholes original abandonment option formulation. Since in practical applications the abandonment value is not deterministic (*

_{X}*i.e.*, σ

*> 0), the use of Margrabe's model is of advantage when analyzing abandonment options. The value of the option*

_{X}*f*, if positive, is added to the net present value

*NPV'*, given by (29), improving the economic performance of the project. If

*f*

__<__0, the abandonment scheme is not economically feasible, and therefore must be discarded or modified.

The employment of the Black-Scholes equation to compute the economic impact of an abandonment option such as the one analyzed in this paper is not free from criticism.

**5 THE TECON CONTAINER TERMINAL **

**5.1 An overview of the terminal**

In order to illustrate the practical application of the methodology to a container port terminal located in a developing country, we have used operational information obtained from the Port of Rio Grande, located in the south of Brazil. Its container terminal – *Tecon *– is operated on a 24-hour/day basis by a private company, Wilson and Sons, since 1997, under a 25-year concession period, renewable for an equal time, and granted by the Brazilian Government. The concession is of the "landlord port" type, in which the domain of the port authority is restricted to the provision of the infrastructure, while investment in the superstructure and terminal operation is the responsibility of the licensed private company (Cullinane *et al.*, 2005; Cullinane & Song, 2002).

Brazil is following a trajectory of higher economic growth, with rates that are expected to surpass 5% a year. Part of the economic growth is fueled by international trade, which is expected to continue to go up in line with Brazilian economic expansion. Brazil's exports by value and by volume have grown significantly during the last thirty years. More importantly, the data indicates that export value has increased more rapidly than export volume. For example, the index of manufactured goods by value, increased almost 23-fold over the period examined. The volume of these exports, however, increased only about 12-fold over the same period of time. Overall, the value of total exports grew 50% more than their volume. Thus, the opportunity cost of marginal time spent in logistics activities has become higher. In this context, time, reliability, just-in-time delivery, and security are critical considerations now. Therefore, the ability to trade with other countries and to grow competitively will be increasingly intertwined with the country logistics system (The World Bank, 2008).

In this context, ports are a key asset in Brazil's logistics context, serving the entire coastline, one of the longest in the world. However, in spite of the locational advantage, the port system suffers from several critical problems that impair its development and contribute to high logistics costs throughout the economy. These include equipment obsolescence, inefficiencies in labor development and labor allocation, lack of harbor capacity, and inefficiencies in port administration. Waiting times for berthing are too high and yard space is insufficient in several ports (The World Bank, 2008). Thus, the employment of robust and up-to-date techniques to analyze terminal expansions and its economic evaluation can lead to better port performance, promoting higher efficiency and competition.

The main products exported through the *Tecon *are tobacco, thermoplastic resins, shoes, frozen poultry, furniture, leather, and auto-parts. Presently, the terminal has three berths, with a total extension of 850.0 m, and 12.5 m depth. Container ship waiting times are greater than the ones observed in other international leading ports of similar capacity. For instance, 33% of container ships that arrived in Rio Grande during 2008, waited more than 12 hours, and 14% waited more than 24 hours. These figures are similar to the results shown by the Port of Santos, the largest in the country, where 35% of ships waited more than 12 hours and more than half of these waited more than 24 hours in the same time period (The World Bank, 2008). In 2008, the *Tecon *operated 356,575 containers, equivalent to 598,196 TEUs.

**5.2 Data analysis**

Although employing a good deal of real data in the analysis, obtained directly from the port administration, we have added some hypothetical variables in order to simulate the future development of the terminal. Table 1 shows the distribution of containers (both load and unload), and expressed in TEU among the various ship capacity classes, at the *Tecon *in 2008. About 45% of the container flow was performed by ships in the capacity range 1500-3000 TEU. On the other hand, about 37% of the movement was done by Post-Panamax ships in the capacity range 4000-5500 TEU. In 2008, the average number of containers loaded/unloaded per ship call was approximately 695 TEU.

It can observed in Table 1 that the number of containers operated per ship call, in absolute terms, and considering both loading and unloading, tends to increase with ship size (column *d*).

However, when one analyzes the ratio between the number of TEUs operated per ship call, and the corresponding ship capacity, this index drops expressively with ship size. For instance, for ships in the range 5000-5500 TEU, the number of containers operated per ship call is only 21% of the ship capacity, while for vessels in the 500-1000 TEU range, the ratio is more than double that figure.

Another important index, considered when analyzing terminal productivity, is the average number of containers loaded/unloaded per ship-hour. Let δ represent the fraction of forty-foot containers on the total number of boxes. The berth productivity index, expressed as a number of containers transferred per hour, per ship (load/unload), can be estimated as

where *n _{TEU} *is the mean number of TEUs transferred per ship call (loaded/unloaded) and

*E*[

*t*] is the average berth occupancy time. The variable did not vary significantly in the period 2001-2009, with a trend toward δ = 0.70, value that has been adopted in the application.

_{s}On the other hand, *E*[*t _{s}*] can be estimated as a function of

*n*, as depicted in Figure 2, and expressed approximately as

_{TEU}

It can been seen that there is a decreasing return to scale in terms of the time necessary to load/unload a ship, as a function of *n _{TEU}*. As a general rule, the container ships travelling south-bound in the Atlantic East Coast of South America stop in Santos and proceed directly to Buenos Aires, Argentina. In the way back, they stop in Montevideo (Uruguay), and next in Rio Grande. Thus, the vessels arrive at

*Tecon*with a light load, generally carrying only the containers embarked in Buenos Aires and Montevideo, plus the units to be discharged at the

*Tecon*. One important characteristic of the Rio Grande's

*Tecon*operation observed in the last years is that the growth of demand is satisfied mostly by increasing the number of ship calls, instead of augmenting the number of containers per stop. As a result, the

*Tecon*depth of 12.5 m is not a restriction today and probably not in the near future too, but it could be a severe constraint if the shipping lines substantially change their operating strategy such as, for instance, the establishment of a hub pole in Rio Grande (SCP, 2006).

Figure 3 shows the monthly evolution of the *Tecon *throughput from January 1999 to December 2009. Up to December 2003 the throughput (expressed in TEU) grew at an average rate of 21% a year. Due to the recent international crisis, the *Tecon *throughput evolution showed a downward inflection, but the trend is already changing, as seen in Figure 3. The volatility of the monthly series of throughput, from January 1999 to December 2009, is now computed. To do this (Hull, 1997) one converts these absolute monthly values *D*(*t*) = *t *= 0, 1, 2,... into relative returns *u*(*t*) = *D*(*t*)/*D*(*t *– 1), *t *= 0.1, 2,..., and then taking the natural logarithms, thus forming a series

The coefficient of variation of *u*(*t*) in the series is 0.201. Converting it to annual basis, one gets the demand volatility *cv* = σ/µ = 0.201 = 0.70, which will be used in our analysis.

Ship inter-arrival times at the *Tecon *are considered next. The analysis of 846 container ship arrivals in 2008 led to a mean inter-arrival time of 10.3 hours. Figure 4 shows the fitting result of a negative exponential distribution to the data. Ship service time, on the other hand, comprises the period of time the vessel effectively occupies a berth. Analyzing the same 2008 ship sample, the mean service time was 13.7 hours. An Erlang distribution with *k *= 6 was fitted to the data, as shown in Figure 5. Although similar studies reported in the literature have indicated Erlang distributions of order 3 (Huang *et al.*, 1995, 2007), the *Tecon *result can be justified by the shorter spread of the number of containers loaded/unloaded per ship call, as compared to other containerports around the world.

In 2008, the mean ship waiting time, computed for the sample of 846 container ship arrivals, was 12.3 hours. Cosmetatos (1976) equation (27) was applied to the observed 2008 data. The number of berths was *c *= 3, the average berth utilization factor was ρ = 0.45, *k *= 6, and π_{0} = 0.250, yielding *q *= 0.071. Thus, the expected ship waiting time is *T _{q} = q *×

*E*[

*t*] = 0.071 × 13.7 = 0.97 hours. This means that ship waiting time at the

_{s}*Tecon*was twelve times the standard figure in 2008. This result is perhaps questionable since the third berth was put into service only in October 2008, but a report from the World Bank confirms the excessive ship waiting time at the

*Tecon*(The World Bank, 2008).

**6 TECON CAPACITY EXPANSION PLANNING **

**6.1 Model structure **

As seen in Section 2, the objective of the capacity expansion modelling is to minimize the net present value of the sum of investments, with no demand backlog and ship waiting times respecting upper limits imposed by best maritime practices observed internationally. In Subsection 2.1, the following relation was defined

In addition to the objective function (44) to be minimized, constraints limiting ship waiting time are added to the model. These constraints are based on best maritime practices among leading container terminals of the world, in a benchmarking process, but since the resulting improvement will not be instantaneous, their upper values are assumed to continuously decrease along time. Let (*t*) be the maximum tolerable ship waiting time at instant *t*. Thus, the estimated ship waiting time at time *t *must respect the constraint

The berth loading/unloading productivity index is represented by the mean number of containers that are loaded/unloaded per ship, per hour. Considering a project lifetime of *T *= 25 years, we divided it into five-year periods, as indicated in Table 2, column (a). Column (b), of Table 2, shows the maximum ship expected capacity to be observed in each period. Column (c) shows the required berth length to receive ships with the capacities listed in column (b). Column (d) indicates the mean ship capacity in each epoch. Column (e) shows the expected number *n _{TEU}* of TEU operated per ship call in each time horizon. The mean berth occupancy time

*E*[

*t*] is exhibited in column (f). Finally, the berth productivity index, represented by the mean number of containers loaded/unloaded per ship-hour, are shown in column (g). These figures are the expected coefficients to prevail during the lifetime of the project. The growth of the variable

_{s}*n*at the

_{TEU}*Tecon*was assumed to follow a polynomial expression with

*t*

with *t *expressed in years. The values of *n _{TEU} *computed via (46) for the time horizons of the project are exhibited in column (e), of Table 2. For each

*t*, the value of

*n*(

_{TEU}*t*) is computed. Thus, with (42), the variable

*E*[

*t*] is estimated. Finally, the values of

_{s}

*Tecon*in2003 was 21 containers/ship/hour (SCP, 2006), which has improved substantially during the last years, reaching about 30 containers/ship/hour in 2008. Thus, the productivity levels indicated in Table 2 represent, in fact, feasible targets, and were adopted in the optimization model.

**6.2 The equivalent discount rate**

First, the project basic discount rate was computed based on the *Capital Asset Pricing Model (CAPM)*, developed by Sharpe and Lintner (Campbell *et al.*, 1997). The *CAPM *was developed, at least in part, to explain the differences in risk premium across assets. Let *r*'_{f} be the discrete-time, annual rate of return on the risk-free asset. Consider, next, a portfolio involving all types of assets traded in the local economy. Let be the discrete average rate of return of such a portfolio. The risky assets contemplated by the investor have returns that are not known with certainty at the time the investments are made. The *CAPM *asserts that the correct measure of risk of the asset is a coefficient called *beta*. The basic equation is

where *r*' is the expected rate of return, *r*'_{f} is the risk-free rate of return, and β is a measure of the relative risk of the asset with reference to the whole portfolio (Campbell *et al.*, 1997).

The risk-free return is the interest rate offered by entities that are entirely creditworthy during the period of a loan, such as the rates paid by the US Treasury and by some European government bonds. In the *Tecon *application we have considered Brazilian figures, assuming *r*'_{f }= 0.036 or 3.6% a year, * * = 0.10 and β = 1.2, leading to *r*' = 0.1128. The corresponding continuous-time discount rate is *r *= ln(1 + *r*') = 0.1069. The equivalent discount rate to be used in the model is given by (19), putting σ/µ = *cv* = 0.70, which is the volatility of the throughput series (Manne, 1961; Freidenfelds, 1981; Novaes & Souza, 2005). Applying (19) with *r *= 0.1069, one has *r** = 0.1048, value to be used in our analysis. As seen in Subsection 2.3.2, the impact of the stochastic behaviour of demand on the results is totally reflected by the discount rate reduction from *r *to *r** (Higle & Corrado, 1992; Bean *et al.*, 1992).

**6.3 Investments on berth expansions**

In 2008, a 250-meter berth was put to service at *Tecon*, with a cost of US$ 50 million. Assuming an additional of US$ 10 million to cover the acquisition of container cranes and other equipments, and dividing the total investment by the berth length, one has a unit cost of US$ 240,000 per meter. Further expansion investments were assumed to be proportional to berth length, given in column (c), Table 2.

It is assumed, in the model, that the mean ship waiting time constraint (45) follows an exponentially decaying function, with = 12 hours at time * **t *= 0, reaching = 2 hours at * **t *= 10 years, and = 1 at * **t *= 14 years, with = 1 hour afterwards. In fact, due to budget restrictions, inefficiencies in port administration, lack of coordination between the terminal and shipping lines and cargo forwarders, among other factors, it is expected that * **Tecon *performance will not improve in a very short time period.

It was admitted that a berth investment is linearly depreciated in 30 years, with nil value at the end of its lifespan. Suppose a berth is installed at time *t*, with an investment *I _{t}*. The salvage value of the berth at the termination time

*T*of the project is (Fig. 6)

with *v* = 0 otherwise, where *T _{B} *is the berth depreciation time.

**6.4 The Dynamic Programming Algorithm**

The dynamic programming concept has been broadly used in the literature to solve sequential decision problems of many kinds. Specifically, it can be used to optimize decision processes over time, which is the case of the problem treated in this paper. Dynamic Programming is, in essence, a general approach to optimization, rather than a mathematical algorithm as, for instance, the Simplex method used to solve linear programming problems. One essential feature of the dynamic programming approach is the reduction of a multivariate problem to a succession of single variable cases. At each stage of the optimization process the analysis of the resulting values of these single variables, with the appropriate decision rules, defines a *policy*. The procedures of dynamic programming are based on the *Principle of Optimality *due to Richard Bellman: an optimal policy has the property that, whatever the initial state and the initial decision are, the remaining decisions must constitute an optimal policy with regard to the state resulting from the first decision (Bellman & Dreyfus, 1962).

In the dynamic programming algorithm developed for this application the *stages *represent the instants *t *(0 __<__ *t *__<__ *T *) for which the model analyses the possibility of introducing additional berths to the container terminal. The *states *of the process, on the other hand, represent the number of berths that will be available at stage *t*, including the ones installed at previous stages plus the berths to be added at time *t*. At each stage of the process the minimum number of required berths *M*_{min} is externally imposed in order to guarantee ρ(*t*)< 1, and *W _{q}*(

*t*)

__<__

*t*). These restrictions help in reducing computational time, since the model investigates only feasible states. Thus, at each stage, the model first computes the minimum required number of berths. Let

*c*(

_{e}*t*) > 0 be the number of berths already in operation at the beginning of stage

*t*, and let

*m*(

*t*) be the quantity of berths to be added at that stage. The following procedure yields the value of

*M*

_{min}:

Step 0: Initialize *c _{e}*(

*t*) and make

*m*(

*t*) ← 0.

Step 1: Make *c*(*t*) → *c _{e}*(

*t*) +

*m*(

*t*) and compute ρ(

*t*) = λ(

*t*)/[µ(

*t*)

*c*(

*t*)] and

*W*(

_{q}*t*) =

*q*·

*E*[

*t*].

_{s}Step 2: If ρ(*t*) __>__ 1 or *W _{q}*(

*t*) > (

*t*), then make

*m*(

*t*) ←

*m*(

*t*) + 1 and go back to Step 1.

Step 3: Return *M*_{min} ← *c*(*t*).

Assuming the two operating restrictions ρ(*t*) < 1 and *W _{q}*(

*t*)

__<__

*t*) are respected, the maximum number of berths to be considered in the analysis is imposed by external considerations, mainly budget constraints and physical expansion limitations. In the former procedure the average ship arrival rate is defined as λ(

*t*) =

*D*(

*t*)/(360 · 24 ·

*n*) and the berth average servicing rate is µ(

_{TEU}*t*) = 1/

*E*[

*t*], where

_{s}*n*is given by (46) and

_{TEU}*E*[

*t*] by (42). All states that satisfy ρ(

_{s}*t*)< 1 and

*W*(

_{q}*t*)

__<__

*t*) are technically feasible, and therefore, can be evaluated according to the objective function (44).

The dynamic programming transitions from the generic stage *t*_{1} to the generic stage *t*_{2} are depicted in Figure 7. In accordance to this evolving structure, the problem is solved in backward steps adopting the following algorithm (Hastings, 1973):

Step 0: Initialize the values of *c _{e}*(0),

*M*

_{max}, and other parameters of the problem.

Step 1: Assign value zero to the terminal states.

Step 2: Repeat Steps 3, 4, 5 e 6 for each instant *t*.

Step 3: Compute *M*_{min} and repeat Steps 4, 5 e 6 for each feasible state.

Step 4: Repeat Steps 5 e 6 for each capacity expansion decision *m*(*t*), until the maximum level *M*_{max} is reached.

Step 5: Compute the objective function value for decision *m*(*t*), at the present stage *t*.

Step 6: If the present decision is better than the previous one, label it as temporary optimal.

Step 7: Stop. Return the optimal capacity expansion sequence.

In our application the maximum number of berths considered in the analysis is ten, meaning a total number of decision variables equal to twenty: (a) the time *t *of capacity expansion in number of ten, and (b) the number of berths *m*(*t*) to be added simultaneously at time *t*, also equally limited to ten. Of course, in terms of dynamic programming this is a small size problem. When dealing with large problems it is possible that the classical dynamic programming approach be limited by dimensional restrictions. In such a case, approximate dynamic programming algorithms (Powell, 2007) can be used.

**6.5 Capacity expansion results**

It is assumed that the throughput *D*(*t*), expressed in TEU, will grow at an annual discrete rate of 10% a year during the first five years, 8% a year during the next ten years, and 6% a year thereafter. The terminal is already in operation at time *t *= 0, with three berths and an annual throughput of *D*_{0} = 650, 000 TEU/year. The project lifespan is *T *= 25 years, and berths and transferring equipments (gantry cranes) are assumed to depreciate in 30 years. The expected investment required to install a new berth, with the necessary cranes, is approximately US$ 240,000 per meter. Depending on the epoch to build a new berth, the expected ship size increases in accordance to column (b), Table 2, and the required berth length will vary accordingly, as shown in column (c) of the same table. It was assumed a 75% experience curve (see Subsection 2.2). In order to estimate average ship waiting times, the Cosmetatos formula (27) was employed, assuming a *M*/*E*_{6}/*c *queue type, with a threshold value of one hour to take into account ship maneuvering, customs, and other formal activities when ship arrives and departs.

With the above assumptions, the resulting optimal expansion plan is exhibited in Table 3. Apart from the three berths already in operation at time *t *= 0, it would be necessary to add four new berths, at different epochs during the project lifespan, as shown in Table 3. The capacity expansion at year 9, month 8, has benefitted from scale effect: instead of one additional berth, the model pointed out an expansion of two berths, with the second costing 75% of the first one. Figure 8 shows the evolution of the number of berths required along the project lifespan.

**7 ECONOMIC FEASIBILITY ANALYSIS **

**7.1 The basic scenario **

Suppose a private firm is negotiating, with the public port agency, the concession to operate the *Tecon *terminal for a period of 25 years. At the end of the concession period the government agency agrees to fully reimburse the private operator with a sum equal to the non-depreciated part of the investments made on berths and equipments, following the scheme indicated in Subsection 6.3. The port agency also requires that the proponent pay a lump sum of 65 million dollars to the government at the signature of the contract.

The economic evaluation to be described next was performed considering the point of view of the proponent firm. Of course, a similar analysis should be performed by the public port agency, considering not only its own economic and financial objectives, but the social impacts as well. Such analysis is not done here due to lack of space and because it would involve considerations of quite different nature. In addition, if the private operator will decide to abandon the project at a certain date, it would mean that the market value of the project at that moment would perhaps not cover the non-depreciated part of the investments. Then, the government agency might not accept to reimburse the private operator to cover the non-depreciated part of the investments fully. But contractual clauses are binding and depend on other factors as, for example, public incentives to attract private investments, political strategies, etc. Of course, other contractual terms can be tested with this methodology, which is not done here due to space limitations.

In order to compute the *NPV' *value of the project, the prospective operator estimates a net margin of 11 dollars per transferred *TEU*. Thus, the value of *R*(*t*) – *E*(*t*) in relation (29) is simply 11 × *D*(*t*), where *D*(*t*) is the throughput at time *t*, expressed in TEUs. The volatility σ* _{p} *of the project is obtained with relation (40). The values of σ

*and σ*

_{S}*can be estimated bysimulation. In this application, however, we assume hypothetical values due to space limitations. Let σ*

_{X}*= 0.70 be the volatility of the project's value-in-use. With regard to σ*

_{S}*, the abandonment option is composed by two parts: the present value of the cash flow accrued by the project during the period 0 <*

_{X}*t*

__<__

*T*and the present value of the salvage value to be accrued at time

_{A}*T*. The latter has a deterministic value in this application, since, under the contractual terms, the government accepts to pay a pre-defined amount of money to the private operator at the end of the deal. Thus, one has σ

_{A}*< σ*

_{X}*, say σ*

_{S}*= 0.40. The correlation between σ*

_{X}*and σ*

_{X}*was assumed to be ρ = 0.4, leading to σ*

_{S}*= 0.72.*

_{p}Making *r*' = 11.28% a year (discrete), or equivalently *r *= 0.1069 (continuous), and considering the investments depicted in Table 3 plus the down payment of *I*_{0} = 65 million, relation (29) yields *NPV *=–0.43 million dollars, meaning the project is economically unfeasible.

In the sequel, the prospective operator analyzes the possibility of proposing an abandonment option to the port agency. The discounted value of the project asset abandonment value is

where 0 __<__ φ_{X} __<__ 1 is a reducing factor to be applied to the salvage value, and *i _{A} *indicates the order of the last investment such that

*t*

_{i}A__<__

*T*. The prospective operator will try to increase the value of φ, while the port agency will argue in the opposite direction. In fact, the state agency wants to discourage an early abandonment of the project by reducing the value of φ

_{A}*as much as possible. At the end of the discussion, both parties agree to set ϕ*

_{X}*= 0.8. The discounted value of the project's value-in-use, on the other hand, is*

_{X}Now, the whole salvage value of the project will be reimbursed to the private operator, as agreed, covering the non-depreciated part of the investments at time *T *.

The prospective operator tried different values of *T _{A}*, but the project remained infeasible for

*T*< 17 years. For

_{A}*T*= 17 years and

_{A}*r*= 0.05827 (equivalent to

_{f }*r*'

_{f }= 0.06 discrete discount rate), the Margrabe model produced

*S*= 91.21 million and

*X*= 30.29 million. Applying relation (39) with σ

*= 0.72 and the above values of*

_{p}*S*and

*X*, one gets

*f*= 17.21 million, which is the value of the put/call option according to Margrabe's method. Finally, the modifiednet present value of the project is equal to

*NPV'*= –0.43 + 17.21 = 16.77 million (Fig. 9). With

*T*= 17 years the project is feasible now, and the project should be submitted to a sensitivity analysis before it is presented to the port agency for further discussion and approval.

_{A}

**7.2 Sensitivity analysis**

A sensitivity analysis was performed over the basic scenario described in Subsection 7.1. Three experience factors δ (Subsection 2.1) were assumed: 0.75, 0.80 and 0.85. On the other hand, three demand forecasting curves were analyzed, as indicated in Table 4. The basic scenario corresponds to δ = 0.75 (Subsection 2.2) and to "intermediate" demand evolution pattern.

The model produced returns to scale for δ = 0.75 when adopting the "intermediate" and the "optimistic" demand growing pattern. The second berth expansion event, for both cases, required two berths implemented simultaneously, as depicted in Table 3 for the "intermediate" situation. For the "pessimistic" growth pattern no return to scale was registered. The same financing criteria adopted for the "intermediate" demand growing pattern has been applied to the "pessimistic" and "optimistic" cases. As expected, the latter produced better *NPV' * result, but the "pessimistic" case showed a negative *NPV' * value. If the lump sum to be paid to the port agency would be reduced from 65 million dollars to 60 million, with all the other conditions unchanged, the project would be feasible for the three demand pattern scenarios. The economic result obtained for δ = 0.80, on the other hand, although smaller when compared with the δ = 0.75 scenario due to the inexistence of returns to scale, repeated the same scheme, with the "pessimistic" pattern requiring a 5 million dollars reduction in the lump sum payment. Changing from δ = 0.80 to δ = 0.85, the results were the same since no returns to scale were incorporated into the model. Thus, in order to have a positive result for the three alternative scenarios depicted in Table 4, the prospective operator could argue to get a reduction of 5 million dollars in the down payment, maintaining all the other factors defined in Subsection 7.1 unchanged.

**8 CONCLUSIONS **

We have analysed here, in considerable detail, the long-term planning of a container terminal in a developing country, with a model that incorporates the berth capacity expansion optimization, the ship waiting times, and the economic feasibility analysis of the investment, where the prospective operator makes use of an abandonment option to add value to the project and turn it economically feasible. Most part of the data was obtained from the *Tecon *container terminal, Port of Rio Grande, Brazil. In developing countries like Brazil, the port operational data is usually incomplete, and therefore the planning of future operations and marketing strategies are limited. On the other hand, the allocation of funds for expansions and maintenance heavily depends on political factors. With the growing trend to transfer at least part of port operations to private organizations, it becomes necessary to evaluate in more detail the economic perspectives of such endeavours, even considering data limitations and other restrictions. Thus, a model framework as the one presented in this paper, even with its limitations, can help the parties in the process of getting a satisfactory agreement.

**ACKNOWLEDGMENTS **

This research has been supported by the Brazilian CNPq/Fapesc/Pronex Project no. 1.0810-00684 and by Capes/German Research Foundation (DFG), Bragecrim Project no. 2009-2.

**REFERENCES **

[1] BEAN JC, HIGLE JL & SMITH RL. 1992. Capacity expansion under stochastic demands. *Operations Research*, **40**(2): S210–S216. [ Links ]

[2] BELLMAN R & DREYFUS SE. 1962. *Applied Dynamic Programming. *Princeton University Press, Princeton, New Jersey. [ Links ]

[3] CAMPBELL JH, LO AW & MACKINLAY AC. 1997. *The Econometrics of Financial Markets. *Princeton University Press, Princeton, New Jersey. [ Links ]

[4] COSMETATOS GP. 1976. Some approximate equilibrium results for the multi-server queue (M/G/r).* Operational Research Quartely*, **27**: 615–620. [ Links ]

[5] COUTO JP & TEIXEIRA JC. 2005. Using linear model for learning curve effect on highrise floor construction. *Construction Management and Economics*, **23**: 355–364. [ Links ]

[6] CULLINANE K & SONG DW. 2002. Port privatization policy and practice. *Transport Reviews*, **22**(1): 55–75. [ Links ]

[7] CULLINANE K, JI P & WANG T. 2005. The relationship between privatization and DEA estimates of efficiency in the container port industry. *Journal of Economics and Business*, **57**: 433–462. [ Links ]

[8] DAGANZO CF. 1990. The productivity of multipurpose seaport terminals. *Transportation Science*, **24**(3): 205–216. [ Links ]

[9] DASCHKOVSKA K & SCHOLZ-REITER B. 2009. How can electronic seals contribute to the efficiency of global container system? *Proceedings Second International Conference on Dynamics in Logistics **– LDIC 2009*, 144-53, Bremen, Germany. [ Links ]

[10] FOGLIATTI MC & MATTOS NMC. 2007. *Queuing Theory. *Editora Interciencia, Rio de Janeiro (in Portuguese). [ Links ]

[11] FREIDENFELDS J. 1980. Capacity expansion when demand is a birth-death process. *Operations Research*, **28**: 712–721. [ Links ]

[12] FREIDENFELDS J. 1981. *Capacity Expansion: Analysis of Simple Models with Applications. *North-Holland, New York. [ Links ]

[13] GOSS RO. 1967. Towards an economic appraisal of port investments. *Journal of Transport Economics and Policy*, September, 1967, 249–271. [ Links ]

[14] GÜNTHER HO & KIM KH. 2005. *Container Terminals and Automated Transport Systems. *Springer, Berlin – Heidelberg. [ Links ]

[15] HASTINGS NAJ. 1973. *Dynamic Programming with Management Applications. *London: Butterworths, London. [ Links ]

[16] HIGLE JL & CORRADO CJ. 1992. Economic investment times for capacity expansion problems.* European Journal of Operational Research*, **59**: 288–293. [ Links ]

[17] HUANG WC, CHISAKI T & LI G. 1995. A study on the container port planning by using cost function with IND based on system. *Journal of the Eastern Asia Society for Transportation Studies*, **1**: 263–276. [ Links ]

[18] HUANG WC & WU SC. 2005. The estimation of the initial number of berths in a port system based on cost-function. *Journal of Marine Science and Technology*, **13**(1): 35–45. [ Links ]

[19] HUANG WC, KUO TC & WU SC. 2007. A comparison of analytical methods and simulation for container terminal planning. *Journal of Chinese Institute of Industrial Engineers*, **24**(3): 200–209. [ Links ]

[20] HULL JC. 1997. *Options, Futures, and Other Derivatives. *Prentice-Hall, Upper Saddle River, New Jersey. [ Links ]

[21] KIA M, SHAYAN E & GHOTB F. 2002. Investigation of port capacity under a new approach by computer simulation. *Computers & Industrial Engineering*, **42**: 533–540. [ Links ]

[22] KIM K & GÜNTHER HO. 2007. *Container Terminals and Cargo Systems. *Springer, Berlin-Heidelberg. [ Links ]

[23] KOZAN E. 1997. Comparison of analytical and simulation planning models of seaport container terminals. *Transportation Planning and Technology*, **20**: 235–248. [ Links ]

[24] LEGATO P & MAZZA RM. 2001. Berth planning and resources optimization at a container terminal via discrete event simulation. *European Journal of Operational Research*, **133**: 537–547. [ Links ]

[25] MANNE AS. 1961. Capacity Expansion and Probabilistic Growth. *Econometrica*, **29**: 632–649. [ Links ]

[26] MANNE AS. 1967. Investments for Capacity Expansion: Size, Location and Time-Phasing. Allen and Unwin, London. [ Links ]

[27] MARGRABE W. 1978. The value of an option to exchange one asset for another. *The Journal of Finance*, **XXXIII**(1): 177–186. [ Links ]

[28] MORRISON JR. 2007. Practical extensions to cycle time approximations for the G/G/m-queue with applications. *IEEE Transactions on Automation Science and Engineering*, **4**(4): 525–532. [ Links ]

[29] MUN J. 2002. *Real options analysis. *Wiley, Hoboken, New Jersey. [ Links ]

[30] NORITAKE M. 1978. A study on optimum number of berths in public wharf. *Proceedings of the Japanese Society of Civil Engineers*, **278**: 113–122. [ Links ]

[31] NORITAKE M & KIMURA S. 1983. Optimum number and capacity of seaport berths. *Journal of Waterway, Port, Coastal and Ocean Engineering*, **109**: 323–329. [ Links ]

[32] NOTTEBOOM TE. 2004. Container shipping and ports: an overview. *Review of Network Economics*,** 3**(2): 86–106. [ Links ]

[33] NOVAES AG & FRANKEL EG. 1966. A queuing model for unitized cargo generation. *Operations Research*, **14**(1): 100–132. [ Links ]

[34] NOVAES AG & SOUZA JC. 2005. A real options approach to a classical capacity expansion problem. *Pesquisa Operacional*, **25**(2): 159–181. [ Links ]

[35] NOVAES AG, SCHOLZ-REITER B, SILVA VMD & ROSA H. 2010a. Long term planning of a container port terminal under demand uncertainty and economies of scale. *12*^{th }*World Conference on Transport Research*, Lisbon, July 11-15, 2010. [ Links ]

[36] NOVAES AG, FRAZZON EM, SILVA VMD & ROSA H. 2010b. Dynamic investment planning for container port ventures. Proceedings of ICIL 2010, *International Conference on Industrial Logistics*, Rio de Janeiro, 21–28, March 2010. [ Links ]

[37] PAGE E. 1972. *Queueing Theory in OR. *Butterworths, London. [ Links ]

[38] POWELL WB. 2007. *Approximate Dynamic Programming – Solving the Curses of Dimensionality. *Wiley, New Jersey. [ Links ]

[39] SAKASEGAWA H. 1976. An approximation formula *L _{q} *

__~__αρ

^{β}/(1 – ρ).

*Annals of the Institute of Statistical Mathematics*,

**29**(A): 67–75. [ Links ]

[40] SCHNABEL JA. 1992. Uncertainty and the abandonment option. *The Engineering Economist*, **37**(2): 172–177. [ Links ]

[41] SCP – COORDINATION AND PLANNING SECRETARIAT, GOVERNMENT OF THE STATE OF RIO GRANDE DO SUL, BRAZIL. 2006. A study of regional development and logistics for the State of Rio Grande do Sul (in Portuguese), Consortium Booz Allen – Fipe – HLC, Porto Alegre, RS. [ Links ]

[42] SHABAYEK AA & YEUNG WW. 2002. A simulation model for the Kwai Chung container terminals in Hong Kong. *European Journal of Operational Research*, **140**: 1–11. [ Links ]

[43] SMIT HTJ & TRIGEORGIS L. 2004. Strategic Investments – Real Options and Games. Princeton University Press, Princeton. [ Links ]

[44] SRINIVASAN TN. 1967. Geometric rate of growth of demand. In: MANNE AS. (Ed.), *Investments for capacity expansion*, Allen & Unwin, Sydney. [ Links ]

[45] TABORGA PN. 1969. Determination of an optimal policy for seaport growth and development.* Department of Civil Engineering, Massachusetts Institute of Technology*, Publication R69-43, Cambridge, Mass. [ Links ]

[46] THE WORLD BANK. 2008. *How to decrease freight logistics costs in Brazil*. Report no. 40020-BR, Washington, DC. [ Links ]

[47] TRIGEORGIS L. 1996. *Real Options – Managerial Flexibility and Strategy in Resource Allocation*. The MIT Press, Cambridge, Mass. [ Links ]

Received November 3, 2010 / Accepted May 31, 2011

* Corresponding author