ABSTRACT
We address the problem of expected cost minimization of meeting the uncertain fuel demand during a time planning horizon, where supply is provided by selecting discrete shipments with lead times. Due to uncertainty and the passage of time, corrective actions can be taken such as cancellation and postponement on supply of shipments with associated costs and delays. This problem is modeled as a stochastic multistage capacitated discrete lotsizing problem with lead times. Computational experiments were performed on the resolution of various instances of the model for four information structures of uncertainty. The experimental optimal values and stochastic rating measures obtained show the validity and interest of the stochastic model, as well as the benefits that can be obtained with respect to a deterministic variant of the model that considers average demand.
Keywords:
oil and gas procurement; oil supply chain; stochastic lotsizing; multistage stochastic integer programming; postponement
1 INTRODUCTION
Distribution and storage of primary products are downstream oil supply chain activities. These involve complex logistic planning under uncertainty of product features and resources. In most nonoil producing countries, or where the production is not sufficient to cover the internal demand, it is necessary to import either crude oil, or even refined products, to cover demand. One of the most important and cheapest transportation modes is by ship. This has an important impact in the supply chain, as it is necessary to negotiate not only volume and price, but also how the delivery will be carried out; shipments are of fixed sizes, ship routes are complex and travel times are usually long. This means that supply contracts must be fixed much in advance of the actual times where the fuel will be needed as stated by Moraes and Faria^{20}20. MORAES LA & FARIA LF. 2016. A stochastic programming approach to liquified natural gas planning. Pesquisa Operacional, 36: 151165, http://www.scielo.br/scielo.php?script=sci_arttext&pid=S010174382016000100151&nrm=iso.
http://www.scielo.br/scielo.php?script=s...
. Additionally, the demands may vary significantly from the best forecasts; this stochastic component introduces an additional complexity and is a source of costs.
This paper focuses on the minimization of the expectation of costs incurred in decisions made to meet the uncertain demand of fuel over a finite discrete time planning horizon. The fuel is purchased through distinguishable optional shipments with a given size and a time of delivery delay. Shipment delivery delay is relevant in the planning horizon. Due to the passage of time and the unveil of demand uncertainty, it is possible that an acquired shipment that has not been received is no longer necessary. In this case, it could be decided to cancel its acquisition or postpone its delivery; decisions, which require minimum accomplishment times and associated costs. The motivation for this model is a real problem arising in a stateowned Uruguayan oil company that deals with fuel acquisition under contractual and logistic conditions for the energy sector. The demand that the company faces is uncertain, given that thermal electricity generation, as a complement in an electrical system, is highly dependent on renewable sources^{24}24. TESTURI CE, ZIMBERG B & FERRARI G. 2012. Modelado estocástico múltiple etapa de adquisición de combustible para la generación de electricidad bajo demanda incierta. Tech. Rep. INCO RT 1207, Instituto de Computación, Facultad de Ingeniería, Universidad de la República, https://www.fing.edu.uy/inco/pedeciba/bibliote/reptec/TR1207.pdf.
https://www.fing.edu.uy/inco/pedeciba/bi...
^{), (}^{32}32. ZIMBERG B, TESTURI CE & FERRARI G. 2019. Stochastic modeling of fuel procurement for electricity generation with contractual terms and logistics constraints. Computers and Chemical Engineering, 123: 4963, https://doi.org/10.1016/j.compchemeng.2018.12.021.
https://doi.org/https://doi.org/10.1016/...
. While this particular situation of the company focuses on the application of fuel procurement for thermal generation, the present work discusses the formulation and solution of a more general variant of the problem, which can represent all situations in which oil procurement is carried out by selecting discrete cargo options in the context of uncertain demand.
Among different approaches to assist planning decisions in the oil industry, stochastic programming has been shown to be a very successful technique, as demonstrated by Dempster et al.^{10}10. DEMPSTER MAH, HICKS PEDRÓ N N, MEDOVA EA, SCOTT JE & SEMBOS A. 2000. Planning logistics operations in the oil industry. Journal of the Operational Research Society, 51(11): 12711288, https://doi.org/10.1057/palgrave.jors.2601043.
https://doi.org/https://doi.org/10.1057/...
and Pinto et al.^{22}22. PINTO J, JOLY M & MORO L. 2000. Planning and scheduling models for refinery operations. Computers and Chemical Engineering, 24(9): 22592276, http://www.sciencedirect.com/science/article/pii/S0098135400005718.
http://www.sciencedirect.com/science/art...
. Another successful application is presented by AlOthman et al.^{3}3. ALOTHMAN WB, LABABIDI HM, ALATIQI IM & ALSHAYJI K. 2008. Supply chain optimization of petroleum organization under uncertainty in market demands and prices. European Journal of Operational Research, 189(3): 822840, http://www.sciencedirect.com/science/article/pii/S0377221707006595.
http://www.sciencedirect.com/science/art...
, who developed a multiperiod stochastic supply chain planning model under uncertain market demands and prices. Christiansen et al.^{8}8. CHRISTIANSEN M, FAGERHOLT K, NYGREEN B & RONEN D. 2013. Ship routing and scheduling in the new millennium. European Journal of Operational Research, 228(3): 467483, http://www.sciencedirect.com/science/article/pii/S0377221712009125.
http://www.sciencedirect.com/science/art...
described a comprehensive review on ship routing and scheduling. A stochastic problem of shipping of oil products among ports and the management of product storage was described by Agra et al.^{1}1. AGRA A, CHRISTIANSEN M, DELGADO A & HVATTUM LM. 2015. A maritime inventory routing problem with stochastic sailing and port times. Computers and Operations Research, 61: 1830, http://www.sciencedirect.com/science/article/pii/S0305054815000210.
http://www.sciencedirect.com/science/art...
. Mixedinteger programming models are used for scheduling refinedoil shipping as shown by Ye, Liang, and Zhu^{31}31. YE Y, LIANG S & ZHU Y. 2017. A mixedinteger linear programmingbased scheduling model for refinedoil shipping. Computers and Chemical Engineering, 99: 106116, http://www.sciencedirect.com/science/article/pii/S0098135417300315.
http://www.sciencedirect.com/science/art...
. Xu et al.^{30}30. XU J, QU H, WANG S & XU Q. 2017. A New Proactive Scheduling Methodology for FrontEnd Crude Oil and Refinery Operations under Uncertainty of Shipping Delay. Industrial and Engineering Chemistry Research, 56(28): 80418053, https://doi.org/10.1021/acs.iecr.7b01496.
https://doi.org/https://doi.org/10.1021/...
incorporated the uncertainty of crude shipping delays within crude oil planning and scheduling operations. Oliveira and Hamacher^{21}21. OLIVEIRA F & HAMACHER S. 2012. Stochastic Benders decomposition for the supply chain investment planning problem under demand uncertainty. Pesquisa Operacional, 32: 663678, http://www.scielo.br/scielo.php?script=sci_arttext&pid=S010174382012000300010&nrm=iso.
http://www.scielo.br/scielo.php?script=s...
considered an integrated fuel distribution network design and binary capacity expansion problem under a multiproduct and multiperiod setting.
The fuelsupply problem tackled in this paper can be modeled as an extension of the lotsizing formulation of Wagner and Whitin^{27}27. WAGNER HM & WHITIN TM. 1958. Dynamic Version of the Economic Lot Size Model. Management Science, 5(1): 8996, http://mansci.journal.informs.org/content/5/1/89.abstract.
http://mansci.journal.informs.org/conten...
; and particularly of the variant with variable capacity or discrete lotsizing of Fleischmann^{12}12. FLEISCHMANN B. 1990. The discrete lotsizing and scheduling problem. European Journal of Operational Research, 44(3): 337348, http://www.sciencedirect.com/science/article/pii/0377221790902457.
http://www.sciencedirect.com/science/art...
. For the case in which the parameters are known with certainty (deterministic case), the size of the lot is continuous, and without capacity or with constant capacity, the problem has efficient resolution through dynamic programming as shown by Wagner and Whitin^{27}27. WAGNER HM & WHITIN TM. 1958. Dynamic Version of the Economic Lot Size Model. Management Science, 5(1): 8996, http://mansci.journal.informs.org/content/5/1/89.abstract.
http://mansci.journal.informs.org/conten...
and Wagelmans et al.^{26}26. WAGELMANS A, HOESEL SV & KOLEN A. 1992. Economic Lot Sizing: An O(nlogn) Algorithm That Runs in Linear Time in the WagnerWhitin Case. Operations Research, 40: 145156, http://www.jstor.org/stable/3840844.
http://www.jstor.org/stable/3840844...
. Bitran and Yanasse^{6}6. BITRAN GR & YANASSE HH. 1982. Computational Complexity of the Capacitated Lot Size Problem. Management Science, 28(10): 11741186, http://www.jstor.org/stable/2630946.
http://www.jstor.org/stable/2630946...
established that the variant with discrete sizing is a generalization of the binary knapsack problem, and belongs to the NPhard complexity class.
In the case that the parameters are random (stochastic case), the problem can be formulated by stochastic programming^{5}5. BIRGE JR & LOUVEAUX F. 1997. Introduction to Stochastic Programming. Springer, New York, springer Series in Operations Research and Financial Engineering.. Stochastic programming model decisions are over time interspersed with random events. These events are assumed discrete establishing, by the time dynamics, a tree of scenarios. The objective of the model is to determine an optimal solution that provides coverage for all the scenarios. Haugen, Løkketangen, and Woodruff^{16}16. HAUGEN KK, LØKKETANGEN A & WOODRUFF DL. 2001. Progressive hedging as a metaheuristic applied to stochastic lotsizing. European Journal of Operational Research, 132(1): 116122, http://www.sciencedirect.com/science/article/pii/S0377221700001168.
http://www.sciencedirect.com/science/art...
describe the scenariobased stochastic lotsizing problem. Ahmed, King and Parija^{2}2. AHMED S, KING AJ & PARIJA G. 2003. A MultiStage Stochastic Integer Programming Approach for Capacity Expansion under Uncertainty. Journal of Global Optimization, 26(1): 324, http://dx.doi.org/10.1023/A:1023062915106.
https://doi.org/http://dx.doi.org/10.102...
establish an adjusted extended formulation of the noncapacitated problem and showed that the WagnerWhitin conditions are not satisfied for the stochastic variant. Guan et al.^{13}13. GUAN Y, AHMED S, MILLER AJ & NEMHAUSER GL. 2006. On formulations of the stochastic uncapacitated lotsizing problem. Operations Research Letters, 34(3): 241250, http://www.sciencedirect.com/science/article/pii/S016763770500057X.
http://www.sciencedirect.com/science/art...
^{), (}^{14}14. GUAN Y, AHMED S, NEMHAUSER GL & MILLER AJ. 2006. A branchandcut algorithm for the stochastic uncapacitated lotsizing problem. Mathematical Programming, 105(1): 5584, http://dx.doi.org/10.1007/s1010700505729.
https://doi.org/http://dx.doi.org/10.100...
propose valid inequalities and a branchandcut algorithm for the noncapacitated variant.
Shipping delay or delivery time of lots, due to production, transport or capacity restrictions, is modeled in some deterministic continuous nocapacitated lotsizing problems. Lee et al.^{19}19. LEE CY, C¸ ETINKAYA S & WAGELMANS APM. 2001. A Dynamic LotSizing Model with Demand Time Windows. Management Science, 47(10): 13841395, http://www.jstor.org/stable/822493.
http://www.jstor.org/stable/822493...
present a variant in which demands have a compliance interval that has efficient resolution by dynamic programming. Brahimi and DauzèrePérès et al.^{7}7. BRAHIMI N, DAUZÈREPÉRÈS S & NAJID N. 2005. Capacitated multiitem lotsizing problems with time windows. Tech. rep., Ecole des Mines de Nantes.^{), (}^{9}9. DAUZÈ REPÉRÈS S, BRAHIMI N, NAJID N & NORDLI A. 2005. Uncapacitated lotsizing problems with time windows. Tech. rep., Ecole des Mines de SaintEtienne. present two variants according to whether the lots are or are not distinguishable with respect to delivery times. These authors propose efficient algorithms of dynamic programming for the distinguishable case and for the undistinguishable case when the orderdelivery windows are not inclusive. For these variants, Wolsey^{28}28. WOLSEY AL. 2006. Lotsizing with production and delivery time windows. Mathematical Programming, 107(3): 471489, http://dx.doi.org/10.1007/s1010700506753.
https://doi.org/http://dx.doi.org/10.100...
sets tight extended formulations. Van den Heuvel and Wagelmans^{25}25. VAN DEN HEUVEL W & WAGELMANS AP. 2008. Four equivalent lotsizing models. Operations Research Letters, 36(4): 465470, http://www.sciencedirect.com/science/article/pii/S0167637708000035.
http://www.sciencedirect.com/science/art...
show the equivalence of the lotsizing problem with production time windows. For the stochastic scenariobased case, Huang and Küçükyavuz^{17}17. HUANG K & KÜÇÜKYAVUZ S. 2008. On stochastic lotsizing problems with random lead times. Operations Research Letters, 36(3): 303308, http://www.sciencedirect.com/science/article/pii/S0167637707001460.
http://www.sciencedirect.com/science/art...
establish that the problem can be efficiently solved in the scenario tree size when delivery windows do not intersect in time. Furthermore, Jiang and Guan improve the efficiency of this procedure^{18}18. JIANG R & GUAN Y. 2011. An O(N2)time algorithm for the stochastic uncapacitated lotsizing problem with random lead times. Operations Research Letters, 39(1): 7477, http://www.sciencedirect.com/science/article/pii/S0167637710001422.
http://www.sciencedirect.com/science/art...
. The peculiar aspects of the treated problem that were not found in the literature are the corrective decisions of cancellation and postponement with time delay in a stochastic modeling.
The content of the present work is described below. In Section 2 the algebraic model of the problem is presented in two subsections, first a deterministic variant is described, together with the entities of the problem, and then the stochastic variant is presented. In Section 3 experiments are established to determine validity of the model. In Section 4, conclusions and future work are discussed.
2 MODEL
The modelled problem deals with the minimization of the cost expectation incurred in decisions made to meet the uncertain demand of fuel over a finite discrete time planning horizon. The problem is reduced to a single fuel in order to simplify the proposal without loss of generality. To satisfy the demand, there are optional distinguishable shipments, denoted as cargoes, with a nonsplittable quantity of the acquired product. The cargoes have relevant delivery times in the planning horizon; so that a significant amount of time elapses between the purchase decision and when the cargo is received. Due to the passage of time and the unveil of demand uncertainty, it could happen that at any given time a cargo, that was previously acquired and not received is no longer necessary. In this case it could be decided to cancel its acquisition or postpone its delivery; decisions, which, in turn, have minimum execution times in relation to the time of delivery and associated costs. After attending the demand in a given period, the amount of remaining fuel is stored up to a certain capacity, to satisfy future demand in a flexible manner.
The mathematical formulation of the problem is presented below. First, a deterministic model of the problem is developed, where definitions and general structure are established. After defining the information structure of the uncertainty from a discrete stochastic process, the stochastic variant of the model is presented.
2.1 Deterministic model
The following is a deterministic formulation of the problem, where the data are assumed to be known with certainty. A discrete time sequential decision process is considered, in which the decisions taken at a given period depend only on the information available up to that period. Main entities of the problem are described as index sets in Table 1. The planning time is represented by a set of discrete time periods, T. The set of cargoes, C, is categorized into already acquired cargoes, A, and possible cargoes to be acquired, P. Decisions are made to purchase possible cargoes, or to cancel and postpone already acquired cargoes for different periods. Due to the nature of the problem, these decisions have relevant compliance times in the planning horizon. Each cargo has a minimum delivery time, between the time the purchase decision is made and the cargo is received. Decisions to cancel and postpone an already purchased cargo must be made before a certain minimum time, prior to the receipt of the cargo. In addition, when deciding to postpone a previously purchased cargo, the time elapsed between the original reception period and the postponement period can not be less than a given minimum postponement time. All these constraints cause some latency in the decision making process.
Parameters
Parameters are described in Table 2. As mentioned before, we consider a single product (or multiple products which are interchangeable). The demand for the product in each period, d
_{t} , is known. Due to store constraints, the inventory of the product at the end of each period is restricted between a minimum volume, s, and a maximum volume,
The period in which an already acquired cargo is received, τ^{c}, is fixed, and it is decided in previous acquisitions (i.e. previous model resolutions). Each cargo has a given volume, q ^{c} . Decisions on each cargo have latency times measured in periods. The delivery time of a cargo, γ^{c}, establishes the length of the wait time (measured in periods) between the acquisition decision and the actual arrival of the cargo. The minimum time for cancellation of a cargo, δ^{c}, establishes the minimum number of periods prior to the delivery period at which the cargo may be cancelled. The minimum postponement time for a cargo, ε^{c}, establishes the minimum number of periods after the delivery period in which the postponed cargo can be received. The achievement period of decisions on acquisition, cancellation and postponement must take place within the planning horizon.
For each cargo there are unit costs per volume associated with the decisions to acquire, ca ^{c} , cancel, cc ^{c} , and postpone, cp ^{c} , it. In addition, there is a unit cost associated with storage in each period, h ^{i} .
The already acquired volume that is scheduled to be received in each period is determined by the sum of the cargoes that are received in that period as
this is an auxiliary summary parameter.
Variables
Variables are summarized in Table 3. The variables s
_{t} represent the fuel inventory at the end of each period t. The acquisition decision of each cargo c to be acquired in period t, subject to its delivery period, is modeled by the binary variables
Objective function
Our aim is to determine a minimum cost scheme of inventory, acquisition, cancellation and postponement of fuel cargoes that satisfy demand during the planning horizon. The costs of acquisitions, cancellations and postponements accrue at the time of decision making. The objective function includes budget costs of acquisition, cancellation and postponement, and inventory costs,
where: the expression (1) represents the costs of acquiring possible cargoes, the expression (2) represents the costs of cancellation minus the budgeted acquisition costs of already acquired cargoes that are cancelled, the expression (3) represents the acquisition costs minus the costs of cancellation plus the costs of postponement of the already acquired cargoes that are postponed, and the expression (4) represents inventory costs. The postponement cost includes the subtraction of cancellation costs, since a postponement is represented by a prior cancellation, but it does not incur cancellation costs.
Constraints
The main requirement is to satisfy demand while maintaining the inventory balance with contributions of the fuel previously stored, what was already acquired, what is acquired, and what is cancelled, postponed and available in storage,
where s_{0} is the given initial storage.
The amount stored in each period is constrained between lower and upper bounds,
The amount of acquired fuel that is received in each period is determined by the sum of the cargo acquisitions in the possible range of the corresponding acquisition periods,
If a cargo is acquired, the acquisition has been decided in a single period before or equal to its possible receiving period less its delivery time γ^{c},
The already acquired volume that is cancelled in each period is determined by the cancellations of the cargoes in the possible range of the corresponding cancellation periods,
If a cargo is cancelled, the cancellation is decided in a single period before or equal to its receiving period τ^{c} minus its cancellation time δ^{c},
The postponement of a cargo is modeled by the use of cancellation, that is to say, it is only possible to postpone cargoes that are cancelled. The already acquired volume that is postponed to a certain period is determined by the postponements of the cargoes in the possible range of the corresponding postponement periods,
If a cargo is postponed, it is to be received in a single period subsequent to or equal to its original receiving period τ^{c} plus its delay time ε^{c},
A cargo can be postponed, if its original arrival decision has been cancelled,
Finally, there are domain constraints of the variables
This formulation is a generalization of the discrete lotsizing problem (DLSC) described by Wolsey^{29}29. WOLSEY LA. 2002. Solving MultiItem LotSizing Problems with an MIP Solver Using Classification and Reformulation. Management Science, 48(12): 15871602, http://www.jstor.org/stable/822525.
http://www.jstor.org/stable/822525...
. The problem DLSC is a special case of the formulation in which the acquisition variables are selected among a special ordered set of type one (cf. constraint (8)), the delivery time of each cargo is zero, and there are no decisions of cancellation or postponement. This formulation belongs to the NPhard class, since DLSC belongs to the same class^{6}6. BITRAN GR & YANASSE HH. 1982. Computational Complexity of the Capacitated Lot Size Problem. Management Science, 28(10): 11741186, http://www.jstor.org/stable/2630946.
http://www.jstor.org/stable/2630946...
.
The feasibility of this formulation is conditioned to the timely availability of fuel acquisition to satisfy its demand,
2.2 Stochastic model with uncertain demand
The following model is a stochastic extension of the deterministic model in which demand is a random parameter. Previously, the uncertain information structure was established through stochastic optimization^{5}5. BIRGE JR & LOUVEAUX F. 1997. Introduction to Stochastic Programming. Springer, New York, springer Series in Operations Research and Financial Engineering.. Subsequently, the entities and formulation of the stochastic model are reported.
Uncertain information structure
The uncertain demand is represented by a discretetime stochastic process indexed in the planning periods; in such a way that each stage of the stochastic process is associated to a period. The process is defined in a finite probability space. It is assumed that the demand of the first period is deterministic, and that the demands of the remaining periods are random with known distribution function. The decisions of a period only depend on the outcomes of the random parameters of previous periods. This process is nonanticipatory of the future decisions or the realizations of the random event outcomes. This information structure can be represented by a tree structure with H levels or stages called tree of scenarios^{23}23. RÖMISCH W & SCHULTZ R. 2001. Multistage Stochastic Integer Programs: An Introduction. In: GRÖTSCHEL M, KRUMKE SO & RAMBAU J (Eds.), Online Optimization of Large Scale Systems, Springer Berlin Heidelberg, pp. 581600, http://dx.doi.org/10.1007/9783662043318_29.
https://doi.org/http://dx.doi.org/10.100...
. This is a perfect directed tree, with the root node in period
Each node of the scenario tree describes the state of the process and is identified by a period and a scenario. An alternative abbreviated notation is to identify the nodes by a single index n in a numerable set of nodes, N. For the first period,
Stochastic model formulation
From the scenario tree notation, the formulation of the deterministic model (cf. Section 2.1) is extended into a multistage stochastic optimization scheme considering the uncertainty in demand.
New index sets are established according to Table 4. First of all, the set of nodes and their subset of leaf nodes are incorporated into the formulation. Subsets of the set of nodes are established to abbreviate the denomination of nodes where it is possible to acquire each cargo,
The random parameter and mapping operators on the nodes of the stochastic model are established in Table 5. The demand parameter, d ^{n} , which in the deterministic model depends on the periods, in the stochastic model depends on the nodes of the tree.
In the stochastic model, most decisions depend on the nodes of the tree according to Table 6. In contrast to the deterministic model, the decision to postpone a cargo, into a given period, could be taken in different nodes; this is modeled by variable
The indexes of periods in deterministic parameters or variables are referred to the temporary realization of a node n by t(n). This is the case for the already acquired volume parameter, which is indexed as a _{t(n)} .
From the previous definitions, the formulation of the multistage stochastic optimization model is
s.a.
This formulation takes into account the same structural properties of the deterministic model extended with the information structure of the scenario tree. It minimizes the expectation of acquisition costs (15), cancellation costs minus acquisition costs in case of cancellation (16), postponement costs plus acquisition costs minus postponement costs (a postponement is modeled in conjunction with a cancellation) and storage costs (18).
The constraint (19) sets the volume balance for each node. The lower and upper storage bounds at each node are determined by (20). The amount of acquired fuel that is received in each node is determined by acquisitions of cargoes in the possible range of the corresponding acquisition periods according to (21). The constraint (22) states that each cargo is acquired at a single node in each path from the root node to a node whose period coincides with the receiving period minus the delivery time of the cargo. The fuel previously acquired that is cancelled at each node is determined by the cancellations of the nodes in the path from the root node to the node, whose cancellation periods are less than the delivery period less the cancellation time, according to (23). The constraint (24) states that each cargo to be cancelled is at a single node in each path from the root node to a node whose period coincides with the receiving period minus the cancellation time of the cargo. The postponement of the cargoes is modeled in conjunction with the cancellation, i.e. only cancelled cargoes can be postponed, (27). The already acquired volume that is postponed in a node is determined by the postponements of the cargoes in the nodes in the path from the root to the node for all periods superior to the period of reception plus the delay time of the node, according to (25). The constraint (26) states that each cargo to be postponed is at a single node in each path from the root node to a node in some period greater than the receiving period plus the time of postponement of the node. The formulation does not need nonanticipativity constraints because the direct indexing of the decisions on the scenario tree implicitly establishes this condition.
The stochastic formulation does not guarantee that all solutions that satisfy constraints at initial stages are also feasible in the remaining stages; that is, it does not have complete recourse^{5}5. BIRGE JR & LOUVEAUX F. 1997. Introduction to Stochastic Programming. Springer, New York, springer Series in Operations Research and Financial Engineering.. This formulation is a generalization of the deterministic one. The deterministic formulation is a special case of the stochastic for the case of a unique scenario in which the tree of scenarios is reduced to a path. Therefore, the feasibility of the stochastic formulation is conditioned to the timely availability of fuel to satisfy its demand along the path of each node of the tree,
3 EMPYRICAL ANALYSIS
In order to generate a number of diverse test instances, four scenario tree structures were considered. Each structure, depicted in Table 7, is determined by the number of direct descendants of each node (tree arity) and the number of periods of the planning horizon. The time size of the periods is measured in weeks. For each tree structure with arity g and H periods there are g
^{H  1} escenarios and
In order to show how the stochastic process is considered, the scenario tree structure with arity
Tree structure with arity
The size of each tree structure model instance (number of equations and variables) for a given distribution of already acquired cargoes (A) and possible cargoes to be acquired (P) is shown in Table 8.
Three data instances were generated for each tree structure and cargo distribution, totaling 12 instances. Each instance has an initial storage,
The computational implementation was performed using AMPL^{4}4. AMPL. AMPL Optimization Inc. http://www.ampl.com, online; accessed July, 2018.
http://www.ampl.com...
for the algebraic coding of the stochastic model, and GUROBI 6.5^{15}15. GUROBI. Gurobi Optimization Inc. http://www.gurobi.com/, online; accessed July, 2018.
http://www.gurobi.com/...
for the resolution of the instances through its branch and cut solver. The calculations were carried out on an Intel Core i7 5960X 3.5GHz computer with 20MB cache and 64GB RAM, operating with CentOS7 Linux system.
Optimal values of the instances are shown in Table 9. The optimal value and time of resolution of the stochastic model, denoted as recourse problem, is depicted at attributes RP and Time, the optimal value of the “wait and see” approach is represented with attribute WS, and the uncertainty measure “expected value of perfect information” is shown at attribute EVPI^{5}5. BIRGE JR & LOUVEAUX F. 1997. Introduction to Stochastic Programming. Springer, New York, springer Series in Operations Research and Financial Engineering.. WS represent the lower bound of the optimum in the idealized case of knowing the future with certainty; and EVPI is defined as the difference between RP and WS. It can be seen that EVPI shows a good performance, since large values validate the use of the stochastic programming approach. Unsurprisingly, instances of larger problems take longer execution time. Particularly the instances 10 and 11 reach the time limit of 900 s. Their best feasible solutions do not change after 48 and 294 s, respectively, and both optimality MIP gap’s marginal rate reduction, at the time limit, are
The “expectation of the expected value problem” for each period t, EEV_{t} , is obtained by solving the problem while fixing the solution of its variables up to time t with the solution of a variant of the problem in which the uncertain parameter is substituted by the mean. Then the “values of the stochastic solution” measure corresponding to period t, VSS_{t} , is defined as the difference between each EEV_{t} and RP. The values of VSS_{t} for the periods of the first two instance are reported in Table 10. These values measures the importance of using the distribution of the uncertain outcomes, and they are calculated according to the proposal of Escudero et al.^{11}11. ESCUDERO LF, GARÍ N A, MERINO M & PÉREZ G. 2007. The value of the stochastic solution in multistage problems. TOP, 15(1): 4864, https://doi.org/10.1007/s1175000700054.
https://doi.org/https://doi.org/10.1007/...
. The VSS_{t} values depicted confirm the model timestep advantage with respect to the solution of the expected value problem. For the case of Instance 1 an
Table 11 shows a summary of the optimal solution of the instances. For each instance, it depicts the number of cargoes that are acquired,
4 CONCLUSIONS
In this paper, we proposed a stochastic multistage capacitated discrete lotsizing model formulation for a discrete cargo fuel supply with lead times problem. The decisions of the problem were represented in detail with their delay time, aspect that for cancellation and postponement decisions is not covered in previous literature. The structure of the uncertain information was modeled by a discrete time stochastic process with finite probability, summarized in a scenario tree. Stochastic programming methodology with entities indexed by nodes of the scenario tree was used to formulate the model. The model extends deterministic models of the literature, which implied the revision of the definitions of the variables and the restrictions to take into account the structure of the scenario tree.
Computational experiments where carried out for several instances with scenarios of a variety of sizes and characteristics. Most of these experiments were solved to optimality for the mediumsize generated instances. The experimental results have shown the validity of the model. Results shows a considerable advantage over the expected value of corresponding deterministic models. Although mediumsized instances could be solved, computational times grew significantly, and for scenarios with more stages or more decisions per stage it is important to improve computational efficiency. Therefore, for future work, it is recommended to study reformulations of the stochastic model with tight high level relaxations, and develop resolution techniques according to them.
ACKNOWLEDGEMENTS
The authors thank the anonymous reviewers for their helpful comments. This work has been supported by the Programa de Desarrollo de las Ciencias Básicas (PEDECIBA), Uruguay.
REFERENCES

^{1}AGRA A, CHRISTIANSEN M, DELGADO A & HVATTUM LM. 2015. A maritime inventory routing problem with stochastic sailing and port times. Computers and Operations Research, 61: 1830, http://www.sciencedirect.com/science/article/pii/S0305054815000210
» http://www.sciencedirect.com/science/article/pii/S0305054815000210 
^{2}AHMED S, KING AJ & PARIJA G. 2003. A MultiStage Stochastic Integer Programming Approach for Capacity Expansion under Uncertainty. Journal of Global Optimization, 26(1): 324, http://dx.doi.org/10.1023/A:1023062915106.
» https://doi.org/http://dx.doi.org/10.1023/A:1023062915106 
^{3}ALOTHMAN WB, LABABIDI HM, ALATIQI IM & ALSHAYJI K. 2008. Supply chain optimization of petroleum organization under uncertainty in market demands and prices. European Journal of Operational Research, 189(3): 822840, http://www.sciencedirect.com/science/article/pii/S0377221707006595
» http://www.sciencedirect.com/science/article/pii/S0377221707006595 
^{4}AMPL. AMPL Optimization Inc. http://www.ampl.com, online; accessed July, 2018.
» http://www.ampl.com 
^{5}BIRGE JR & LOUVEAUX F. 1997. Introduction to Stochastic Programming. Springer, New York, springer Series in Operations Research and Financial Engineering.

^{6}BITRAN GR & YANASSE HH. 1982. Computational Complexity of the Capacitated Lot Size Problem. Management Science, 28(10): 11741186, http://www.jstor.org/stable/2630946
» http://www.jstor.org/stable/2630946 
^{7}BRAHIMI N, DAUZÈREPÉRÈS S & NAJID N. 2005. Capacitated multiitem lotsizing problems with time windows. Tech. rep., Ecole des Mines de Nantes.

^{8}CHRISTIANSEN M, FAGERHOLT K, NYGREEN B & RONEN D. 2013. Ship routing and scheduling in the new millennium. European Journal of Operational Research, 228(3): 467483, http://www.sciencedirect.com/science/article/pii/S0377221712009125
» http://www.sciencedirect.com/science/article/pii/S0377221712009125 
^{9}DAUZÈ REPÉRÈS S, BRAHIMI N, NAJID N & NORDLI A. 2005. Uncapacitated lotsizing problems with time windows. Tech. rep., Ecole des Mines de SaintEtienne.

^{10}DEMPSTER MAH, HICKS PEDRÓ N N, MEDOVA EA, SCOTT JE & SEMBOS A. 2000. Planning logistics operations in the oil industry. Journal of the Operational Research Society, 51(11): 12711288, https://doi.org/10.1057/palgrave.jors.2601043.
» https://doi.org/https://doi.org/10.1057/palgrave.jors.2601043 
^{11}ESCUDERO LF, GARÍ N A, MERINO M & PÉREZ G. 2007. The value of the stochastic solution in multistage problems. TOP, 15(1): 4864, https://doi.org/10.1007/s1175000700054.
» https://doi.org/https://doi.org/10.1007/s1175000700054 
^{12}FLEISCHMANN B. 1990. The discrete lotsizing and scheduling problem. European Journal of Operational Research, 44(3): 337348, http://www.sciencedirect.com/science/article/pii/0377221790902457
» http://www.sciencedirect.com/science/article/pii/0377221790902457 
^{13}GUAN Y, AHMED S, MILLER AJ & NEMHAUSER GL. 2006. On formulations of the stochastic uncapacitated lotsizing problem. Operations Research Letters, 34(3): 241250, http://www.sciencedirect.com/science/article/pii/S016763770500057X
» http://www.sciencedirect.com/science/article/pii/S016763770500057X 
^{14}GUAN Y, AHMED S, NEMHAUSER GL & MILLER AJ. 2006. A branchandcut algorithm for the stochastic uncapacitated lotsizing problem. Mathematical Programming, 105(1): 5584, http://dx.doi.org/10.1007/s1010700505729.
» https://doi.org/http://dx.doi.org/10.1007/s1010700505729 
^{15}GUROBI. Gurobi Optimization Inc. http://www.gurobi.com/, online; accessed July, 2018.
» http://www.gurobi.com/ 
^{16}HAUGEN KK, LØKKETANGEN A & WOODRUFF DL. 2001. Progressive hedging as a metaheuristic applied to stochastic lotsizing. European Journal of Operational Research, 132(1): 116122, http://www.sciencedirect.com/science/article/pii/S0377221700001168
» http://www.sciencedirect.com/science/article/pii/S0377221700001168 
^{17}HUANG K & KÜÇÜKYAVUZ S. 2008. On stochastic lotsizing problems with random lead times. Operations Research Letters, 36(3): 303308, http://www.sciencedirect.com/science/article/pii/S0167637707001460
» http://www.sciencedirect.com/science/article/pii/S0167637707001460 
^{18}JIANG R & GUAN Y. 2011. An O(N2)time algorithm for the stochastic uncapacitated lotsizing problem with random lead times. Operations Research Letters, 39(1): 7477, http://www.sciencedirect.com/science/article/pii/S0167637710001422
» http://www.sciencedirect.com/science/article/pii/S0167637710001422 
^{19}LEE CY, C¸ ETINKAYA S & WAGELMANS APM. 2001. A Dynamic LotSizing Model with Demand Time Windows. Management Science, 47(10): 13841395, http://www.jstor.org/stable/822493
» http://www.jstor.org/stable/822493 
^{20}MORAES LA & FARIA LF. 2016. A stochastic programming approach to liquified natural gas planning. Pesquisa Operacional, 36: 151165, http://www.scielo.br/scielo.php?script=sci_arttext&pid=S010174382016000100151&nrm=iso
» http://www.scielo.br/scielo.php?script=sci_arttext&pid=S010174382016000100151&nrm=iso 
^{21}OLIVEIRA F & HAMACHER S. 2012. Stochastic Benders decomposition for the supply chain investment planning problem under demand uncertainty. Pesquisa Operacional, 32: 663678, http://www.scielo.br/scielo.php?script=sci_arttext&pid=S010174382012000300010&nrm=iso
» http://www.scielo.br/scielo.php?script=sci_arttext&pid=S010174382012000300010&nrm=iso 
^{22}PINTO J, JOLY M & MORO L. 2000. Planning and scheduling models for refinery operations. Computers and Chemical Engineering, 24(9): 22592276, http://www.sciencedirect.com/science/article/pii/S0098135400005718
» http://www.sciencedirect.com/science/article/pii/S0098135400005718 
^{23}RÖMISCH W & SCHULTZ R. 2001. Multistage Stochastic Integer Programs: An Introduction. In: GRÖTSCHEL M, KRUMKE SO & RAMBAU J (Eds.), Online Optimization of Large Scale Systems, Springer Berlin Heidelberg, pp. 581600, http://dx.doi.org/10.1007/9783662043318_29.
» https://doi.org/http://dx.doi.org/10.1007/9783662043318_29 
^{24}TESTURI CE, ZIMBERG B & FERRARI G. 2012. Modelado estocástico múltiple etapa de adquisición de combustible para la generación de electricidad bajo demanda incierta. Tech. Rep. INCO RT 1207, Instituto de Computación, Facultad de Ingeniería, Universidad de la República, https://www.fing.edu.uy/inco/pedeciba/bibliote/reptec/TR1207.pdf
» https://www.fing.edu.uy/inco/pedeciba/bibliote/reptec/TR1207.pdf 
^{25}VAN DEN HEUVEL W & WAGELMANS AP. 2008. Four equivalent lotsizing models. Operations Research Letters, 36(4): 465470, http://www.sciencedirect.com/science/article/pii/S0167637708000035
» http://www.sciencedirect.com/science/article/pii/S0167637708000035 
^{26}WAGELMANS A, HOESEL SV & KOLEN A. 1992. Economic Lot Sizing: An O(nlogn) Algorithm That Runs in Linear Time in the WagnerWhitin Case. Operations Research, 40: 145156, http://www.jstor.org/stable/3840844
» http://www.jstor.org/stable/3840844 
^{27}WAGNER HM & WHITIN TM. 1958. Dynamic Version of the Economic Lot Size Model. Management Science, 5(1): 8996, http://mansci.journal.informs.org/content/5/1/89.abstract
» http://mansci.journal.informs.org/content/5/1/89.abstract 
^{28}WOLSEY AL. 2006. Lotsizing with production and delivery time windows. Mathematical Programming, 107(3): 471489, http://dx.doi.org/10.1007/s1010700506753.
» https://doi.org/http://dx.doi.org/10.1007/s1010700506753 
^{29}WOLSEY LA. 2002. Solving MultiItem LotSizing Problems with an MIP Solver Using Classification and Reformulation. Management Science, 48(12): 15871602, http://www.jstor.org/stable/822525
» http://www.jstor.org/stable/822525 
^{30}XU J, QU H, WANG S & XU Q. 2017. A New Proactive Scheduling Methodology for FrontEnd Crude Oil and Refinery Operations under Uncertainty of Shipping Delay. Industrial and Engineering Chemistry Research, 56(28): 80418053, https://doi.org/10.1021/acs.iecr.7b01496.
» https://doi.org/https://doi.org/10.1021/acs.iecr.7b01496 
^{31}YE Y, LIANG S & ZHU Y. 2017. A mixedinteger linear programmingbased scheduling model for refinedoil shipping. Computers and Chemical Engineering, 99: 106116, http://www.sciencedirect.com/science/article/pii/S0098135417300315
» http://www.sciencedirect.com/science/article/pii/S0098135417300315 
^{32}ZIMBERG B, TESTURI CE & FERRARI G. 2019. Stochastic modeling of fuel procurement for electricity generation with contractual terms and logistics constraints. Computers and Chemical Engineering, 123: 4963, https://doi.org/10.1016/j.compchemeng.2018.12.021.
» https://doi.org/https://doi.org/10.1016/j.compchemeng.2018.12.021
Publication Dates

Publication in this collection
09 May 2019 
Date of issue
JanApr 2019
History

Received
30 Nov 2017 
Accepted
31 Oct 2018