Next Article in Journal / Special Issue
Effect of Additional Order in Two-Stage Supply Chain Contract under the Demand Uncertainty
Previous Article in Journal
Hybrid Mamdani Fuzzy Rules and Convolutional Neural Networks for Analysis and Identification of Animal Images
Previous Article in Special Issue
Application of the Generalized Laplace Homotopy Perturbation Method to the Time-Fractional Black–Scholes Equations Based on the Katugampola Fractional Derivative in Caputo Type
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Application of Optimal Control to Sugarcane Harvesting in Thailand

by
Wisanlaya Pornprakun
1,2,†,
Surattana Sungnul
1,2,*,†,
Chanakarn Kiataramkul
1,2,† and
Elvin James Moore
1,2,†
1
Department of Mathematics, King Mongkut’s University of Technology North Bangkok, Pracharat 1 Road, Bangkok 10800, Thailand
2
Centre of Excellence in Mathematics, CHE, Si Ayutthaya Road, Bangkok 10400, Thailand
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 8 February 2021 / Revised: 16 March 2021 / Accepted: 16 March 2021 / Published: 19 March 2021

Abstract

:
The sugar industry is of great importance to the Thai economy. In general, the government sets sugarcane prices at the beginning of each harvesting season based on type (fresh or fired), sweetness (sugar content) and gross weight. The main aim of the present research is to use optimal control to find optimal sugarcane harvesting policies for fresh and fired sugarcane for the four sugarcane producing regions of Thailand, namely North, Central, East and North-east, for harvesting seasons 2012/13, 2013/14, 2014/15, 2017/18 and 2018/19. The optimality problem is to determine the harvesting policy which gives maximum profit to the farmers subject to constraints on the maximum amount that can be cut in each day, where a harvesting policy is defined as the amount of each type of sugarcane harvested and delivered to the sugar factories during each day of a harvesting season. The results from the optimal control methods are also compared with results from three optimization methods, namely bi-objective, linear programming and quasi-Newton. The results suggest that discrete optimal control is the most effective of the five methods considered. The data used in this paper were obtained from the Ministry of Industry and the Ministry of Agriculture and Co-operatives of the Royal Thai government.

1. Introduction

The main aims of this paper are: (1) to use two optimal control methods to determine optimal harvesting policies for sugarcane in Thailand, where an optimal harvesting policy is defined as the amount of sugarcane harvested each day during a crop year that gives maximum profit to farmers; (2) To compare the efficiency and results of the two optimal control methods with three optimization methods. The two optimal control methods we consider are based on the discrete and continuous Pontryagin maximum principles (see, e.g., [1,2,3]). The three optimization methods we consider are the bi-objective ε -constraints method (see, e.g., [4]), a quasi-Newton optimization method (see, e.g., [5,6,7]) and linear programming (see, e.g., [8]). We use the five optimal control and optimization methods to determine optimal harvesting policies for the two types of sugarcane (fresh and fired) in the four regions (North, Central, East and North-east) of Thailand for crop years 2012/13, 2013/14, 2014/15, 2017/18 and 2018/19. We use data obtained from the Ministry of Industry and the Ministry of Agriculture and Co-operatives of the Royal Thai government [9,10,11,12,13]. The present paper is an extension of a previous paper [14] in which bi-objective and quasi-Newton methods were used to find optimal harvesting policies for the crop years 2012/13, 2013/14 and 2014/15. For these two optimization methods, it was necessary to divide a harvesting season into 15-day periods and to assume that there was zero growth of sugarcane during a season.
Although there have been few, if any, applications of optimal control to sugarcane harvesting, there have been applications of optimal control to similar types of problems in other areas which could be adapted to sugarcane farming [3]. For example, in 2011, Kiataramkul et al. [15] studied an optimal control problem for optimal nutritional intake for fetal growth in sheep. In 2014, Dine, Lenhart and Behncke [16] investigated discrete-time optimal harvesting of fish populations with age structure with the two objectives of maximizing the profit of fishing by finding the optimal harvesting strategy for each age class and of finding the optimal net size. In 2015, Kiataramkul and Matkhao [17] studied the optimal control problem of food intake of swine during the post weaning period with the objective of minimizing the amount of food fed to the swine and achieving an optimal weight at time of sale. In 2016, Puengpo et al. [18] applied continuous and discrete optimal control models to find optimal methods for sheep and swine feeding problems.
In contrast to optimal control, there have been a number of applications of optimization methods to sugar production. Some examples of applications of optimization methods to sugarcane harvesting and sugar production include the following. In 2012, Gomes [19] studied a bi-objective mathematical model for choosing sugarcane varieties in Brazil which have a harvest biomass residual that could be used in electricity generation. The bi-objective optimization model aimed to minimize the cost of harvesting the residual biomass and to maximize the revenue from sale of the generated electricity. In 2016, Sungnul et al. [20] studied a bi-objective optimization model to find an optimal time of harvesting for sugarcane growers in the North-east region of Thailand. The aim of this work was to help farmers to find the optimal harvesting time in order to maximize the revenue and to minimize the cost. Sungnul et al. used the ε -constraints method [4] to find the optimal cutting pattern by using the revenue as the objective function and the costs as constraints. In 2017, Sungnul et al. [21] extended the work in [20] to find the optimal harvesting times for all of the four regions of Thailand. Quasi-Newton optimization methods are well-known methods of optimization that have been used for many years to find optimal solutions for problems in many areas of science, finance and industry (see, e.g., [5,6,7]). In 2019, Pornprakun et al. [14] used the bi-objective and quasi-Newton optimization methods (see, e.g., [5,6,7]) to determine optimal harvesting policies for sugarcane in Thailand in order to maximize the total revenue and minimize harvesting cost for crop years 2012/13, 2013/14 and 2014/15.

2. Mathematical Models

We used the following mathematical model as the basic model for maximizing profit from sugarcane harvesting.
max u J [ u ] = t 0 t f P ( t , x ( t ) , u ( t ) ) d t , subject   to d x ( t ) d t = r x ( t ) 1 x ( t ) K u ( t ) , 0 u ( t ) u m a x , x ( t ) 0 , x ( t 0 ) = a , x ( t f ) = 0 ,
where J [ u ] is the profit functional to be maximized by selecting the control variable u ( t ) , which is the rate (tonnes/day) of cutting sugarcane at time t. Also, P ( t , x ( t ) , u ( t ) ) is the profit for sugarcane cut at time t (baht/day),
  • P ( t , x ( t ) , u ( t ) ) is the profit for sugarcane cut at time t (baht/day),
  • x ( t ) is the total amount of sugarcane (tonnes) on the farms in a region at time t,
  • u m a x is the maximum rate of cutting sugarcane (tonnes/day),
  • t 0 is the initial time at start of harvesting season (day),
  • t f is the final time at end of harvesting season (day),
  • a is the total amount of sugarcane on farms at the initial time t 0 (tonnes),
  • r x ( t ) 1 x ( t ) K is a logistic growth function for the rate of increase in weight of sugarcane (tonnes per day) on the farms in a region at time t, where r is a specific growth rate (1/day), and K is a constant which represents the carrying capacity of the farms in the region (tonnes).
We also considered the following discrete version of (1).
max u J [ u ] = k = 0 N 1 P ( t k , x k , u k ) , subject   to x k + 1 = x k + r x k 1 x k K u k , 0 u k u m a x , x k 0 , x 0 = x ( t 0 ) = a , x N = x ( t N ) = 0 ,
where x k = x ( t k ) and u k = u ( t k ) and t k , k = 0 , 1 , 2 , , N 1 are N periods of a cutting season and t N is the end of the cutting season. For the optimal control methods, we used 1-day periods, for the bi-objective and quasi-Newton methods we found it necessary to assume 15-day periods, and for linear programming we considered both 1-day and 15-day periods.

3. Optimization and Optimal Control Methods

In this section, we give details of the three optimization and two optimal control methods that we used to obtain optimal policies for sugarcane harvesting in Thailand.

3.1. Bi-Objective Optimization

For the bi-objective optimization problem, we separated the profit P ( t k , x k , u k ) per 15-day period into a revenue term R ( t k , x k , u k ) and a cost term C ( t k , x k , u k ) and assumed that r = 0 , that is, no growth during the harvesting season. We then used the ε -constraints method [4,14] to solve the bi-objective optimization problem
max u R ( u ) = k = 0 N 1 R ( t k , x k , u k ) , subject   to C ( u ) = k = 0 N 1 C ( t k , u k , x k ) ε r , x k + 1 = x k u k , 0 u k u m a x , x k 0 , x 0 = a , x N = 0 ,
where the ε r are a set of values of costs between a minimum and a maximum value of the cost for a feasible cutting pattern.

3.2. Quasi-Newton Optimization

As in the previous paper [14], we used the constrained optimization function fmincon in Matlab with the “active-set” algorithm based on quasi-Newton method [5,6,7] to find the optimal cutting patterns to maximize the profit in (2) for 15-day periods and we assumed that there was no growth during the harvesting season.

3.3. Linear Programming

We used the Matlab program linprog with the “dual-simplex” algorithm based on linear programming [8] to solve the optimization problem in (2) for both 15-day periods and daily periods. We again assumed that there was no growth during the harvesting season because growth would make the optimization problem a nonlinear problem.
The results of three optimization methods for the years 2012/13, 2013/14 and 2014/15 have been published in [14] and the results for 2017/18, 2018/19 are included in this paper.

3.4. Pontryagin Maximum Principle for Continuous Optimal Control

For a given type, region and harvesting season, the continuous optimal control problem is of the form given in (1). We used the continuous Pontryagin maximum principle to solve this continuous optimal control problem (see, e.g., [1,2,3]). It should be noted that the optimal control problem for sugar harvesting has state variable constraints on the amount that can be cut each day and on the total amount that is available for cutting.
The first step is to define a Hamiltonian
H ( t , x ( t ) , u ( t ) , λ ( t ) ) = P ( t , x ( t ) , u ( t ) ) + λ ( t ) T f ( t , x ( t ) , u ( t ) ) ,
where f ( t , x ( t ) , u ( t ) ) = r x ( t ) 1 x ( t ) K u ( t ) and the λ ( t ) are called costate variables. Then, the equations for the state variables x ( t ) and costate variables λ ( t ) are given by
State : d x d t = f ( t , x ( t ) , u ( t ) ) = H ( t , x ( t ) , u ( t ) , λ ( t ) ) λ ( t ) ,
Costate : d λ d t = H ( t , x ( t ) , u ( t ) , λ ( t ) ) x ( t ) = P ( t , x ( t ) , u ( t ) ) x ( t ) λ ( t ) T f ( t , x ( t ) , u ( t ) ) x ( t ) .
Since the boundary conditions on the state variables in (1) are given values at fixed initial and final times, the appropriate boundary conditions on the state and costate equations are
State : x ( t 0 ) = a , x ( t f ) = 0 given . Costate : λ ( 0 ) , λ ( t f ) free .
Then, the Pontryagin maximum principle states that the optimal control u * ( t ) is obtained by finding the maximum of H ( t , x * ( t ) , u ( t ) , λ * ( t ) ) with respect to u ( t ) , where the x * ( t ) , λ * ( t ) are the solutions of the state and costate equations for u * ( t ) .
In general, the maximum of the Hamiltonian can occur at the minimum or maximum values of the cutting constraints or at an internal point. For the problem in (1) the state equation is a linear function of the control u. In this case, there is no internal optimal value of u and the control is called a “bang-bang” control and the optimal values u * ( t ) are either the minimum u ( t ) = 0 or maximum u ( t ) = u max . However, as noted above, the sugar harvesting problem has the extra complication that the constraints on the control u * ( t ) also depend on the values of the state variables x * ( t ) because of the conditions x * ( t ) 0 , x * ( t f ) = 0 . We have included these conditions as explained in step 4 of the following algorithm.

Continuous Optimal Control Algorithm

The algorithm that we used for the solution of the continuous optimal control problem is as follows.
  • Select an initial feasible cutting pattern u ( t k ) , k = 0 , 1 , , N 1 such that x ( t ) 0 and x ( t f ) = 0 .
  • Solve the state and costate equations for this pattern to obtain values of x ( t k ) and λ ( t k ) . For reasons of numerical stability, the state equations must be integrated forwards in time and the costate equations must be integrated backwards in time.
  • Compute the Hamiltonian values H ( t k , x ( t k ) , u ( t k ) , λ ( t k ) ) and the partial derivative values H ( t k , x ( t k ) , u ( t k ) , λ ( t k ) ) u ( t k ) .
  • In our examples, we have found that all partial derivatives are positive. Therefore, if there were no constraints on the state variables, the optimal cutting pattern would be to cut the maximum amount available for cutting each day. In order to satisfy the state variable constraints x ( t ) 0 and x ( t f ) = 0 , we chose the optimal u ( t k ) values as follows.
    (a)
    Sort the derivatives H ( t k , x ( t k ) , u ( t k ) , λ ( t k ) ) u ( t k ) in decreasing order.
    (b)
    Find the number of days M required to cut all available cane at the maximum rate u max and select the times t k corresponding to the top M values of the derivatives.
    (c)
    If M u max = a , set the new cutting pattern with u ( t k ) = u max at these M values of t k and u ( t k ) = 0 at all other times. However, if M u max > a , set u ( t k ) = u max for the top M 1 derivative values and
    u ( t k ) = a ( M 1 ) u max where t k is time of derivative M.
    Note: A slight modification is required if growth can occur ( r > 0 ).
  • Then, stop if the new cutting pattern is within a selected tolerance of the old cutting pattern or return to step 2 and repeat if the new pattern is not within the selected tolerance of the old.

3.5. Pontryagin Maximum Principle for Discrete Optimal Control

The discrete version of the continuous optimal control problem in (1) is as follows.
max u J [ u ] = k = 0 N 1 P ( t k , x k , u k ) , subject   to x k + 1 = x k + r x k 1 x k K u k , 0 u k u m a x , x 0 = x ( t 0 ) = a , x N = x ( t N ) = 0 ,
where x k = x ( t k ) and u k = u ( t k ) .
As for the continuous case, the first step is to define a Hamiltonian
H ( k , x k , u k , λ k + 1 ) = P ( t k , x k , u k ) + λ k + 1 T f ( k , x k , u k ) ,
where f ( k , x k , u k ) = x k + r x k 1 x k K u k .
Then, the equations for the state variables x k and costate variables λ k are given by
State : x k + 1 = f ( k , x k , u k ) = H ( k , x k , u k , λ k + 1 ) λ k + 1 ,
Costate : λ k = H ( k , x k , u k , λ k + 1 ) x k = P ( k , x k , u k ) x k + λ k + 1 T f ( k , x k , u k ) x k .
Since the boundary conditions on the state variables in (8) are given values at fixed initial and final times, the appropriate boundary conditions on the state and costate equations are
State : x 0 = a , x N = 0 given . Costate : λ 0 , λ N free .
As for the continuous case, the optimal cutting pattern will be “bang-bang”. The algorithm for finding the optimal cutting pattern for the continuous case can also be used for the discrete case.

4. The Commercial Cane Sugar System (CCS) in Thailand

As details of the CCS system have been given in [14,22], we will only give a brief review here. Thai farmers harvest two types of sugarcane, namely fresh and fired, where fired sugarcane is burnt to remove leaves before it is cut so that it can be cut manually by workers, whereas fresh sugarcane is usually cut by machines which can remove the leaves as the cane is cut. In the CCS system, the Royal Thai government sets the price of fresh and fired sugarcane for each of the four sugarcane producing regions in Thailand, namely North, Central, East and North-east. The price of the sugarcane is based on the two main factors of weight and sweetness, where sweetness or CCS is the percentage of sugar in the sugarcane. All data in this paper have been obtained from the Office of the Cane and Sugar Board (OCSB), Ministry of Industry and the Ministry of Agriculture and Co-operatives of the Royal Thai government [9,10,11,12,13]. The OCSB reports data for 15-day periods of a harvesting season which typically starts around 15 November and ends around 15 June of the following year.
In the following, we will use the notation, i = A for fresh and i = B for fired sugarcane. For the four sugarcane-growing regions of Thailand, we will use the notation j = 1 for North, j = 2 for Central, j = 3 for East, and j = 4 for North-east.

4.1. Price of Sugarcane

As stated above, the price of sugarcane is set on the basis of type, weight, and sweetness (CCS) [14,20].
  • Price based on weight: The basic price per weight of sugarcane (baht/tonne) is fixed by the Royal Thai government for each crop year. This basic price is the same for all regions. However, the actual price paid to farmers for fired sugarcane is 30 baht/tonne less than the basic price. Then, at the end of harvesting for the year, factories in each region will share the total amount of money deducted from fired sugarcane sales in that region to farmers who sold fresh sugarcane at a rate not exceeding 70 baht/tonne of fresh sugarcane delivered. We will use the notation P w ( i , j ) to denote actual price based on weight (baht/tonne) for type i in region j.
    The actual prices per tonne based on weight for sugarcane for a given crop year are:
    Fresh : P w ( A , j ) = P w + 30 a ( B , j ) a ( A , j ) P w + 70 , Fired : P w ( B ) = P w 30 ,
    where P w is the basic price, a ( A , j ) is the total amount of fresh sugarcane (tonnes) harvested in region j and a ( B , j ) is the total amount of fired sugarcane (tonnes) harvested in region j for a given year.
  • Price based on sweetness (CCS): Each year, the Royal Thai government sets a basic price per tonne for sugarcane with 10 CCS, where CCS is the percentage by weight of sugar in sugarcane. This price is the same for fresh and fired sugarcane and for all regions. The actual price per tonne received by farmers is then adjusted if the CCS is different from 10. For sugarcane harvested in period k in region j the actual price per tonne for a given year is
    P c ( k , j ) = P c ( 1 + 0.06 y ( k , j ) ) ,
    where P c is the basic price per tonne for sugarcane with 10 CCS set by the government for the year, and y ( k , j ) = C C S ( k , j ) 10 , where C C S ( k , j ) is the average CCS from sugarcane harvested in period k in region j and the factor 0.06 is the rate of change of price per 1 CCS from the base level of 10.
Therefore, the total prices to farmers (baht/tonne) for a given crop year are:
Fresh : P ( k , A , j ) = P w ( A , j ) + P c ( k , j ) , Fired : P ( k , B , j ) = P w ( B ) + P c ( k , j ) .
The basic prices per weight P w and for sweetness P c set by the Royal Thai government and the total amounts delivered to the factories are shown in Table 1. It can be seen from this data that the amount of fresh sugarcane is appreciably less than the amount of fired in all regions in all crop years and that the amount of sugarcane was increasing from 2012/13 to 2017/18 and then decreased in 2018/19.
The CCS values are shown in Table 2 and Table 3 for 15-day periods for all regions for crop years 2017/18 and 2018/19. Similar tables for the crop years 2012/13, 2013/14 and 2014/15 have been published previously in [14]. It can be seen that in 2017/18 the CCS value in all regions reached a maximum in April and then decreased slightly, whereas in 2018/19 the CCS value increased in all regions during the harvesting season.
The total prices for the fresh and fired sugarcane are shown in Table 4 and Table 5 for 15-day periods for all regions for crop years 2017/18 and 2018/19. It can be seen that the prices in all regions increase rapidly at the beginning of a crop year and then slowly at the end due to the changes in CCS values. It can also be seen that the fresh prices are appreciably higher than the fired prices due to the price adjustments discussed in Section 4.1.

4.2. Costs of Production 

The costs of production (baht/tonne) can be separated into cutting costs, transport costs and maintenance costs. In this paper, we assume that the cutting and transport costs are fixed costs that depend only on the amount cut in a given period, whereas the maintenance costs (baht/tonne) are variable costs that depend on the amount of uncut sugarcane remaining on the farm after cutting in previous periods.
We assume that the total cost of sugarcane production (baht) of type i in period k in region j in a given year is
C ( i , j ) = k = 0 N 1 C h ( i , j ) u k ( i , j ) + C t ( i , j ) u k ( i , j ) + C m ( i , j ) x k ( i , j ) ,
where u k ( i , j ) is the weight of sugarcane (tonnes) cut in period k, x k ( i , j ) is the weight of uncut sugarcane remaining on the farm in period k, C h ( i , j ) is harvesting cost of sugarcane (baht/tonne), C t ( i , j ) is transport cost for delivering sugarcane to the mills (baht/tonne) and C m ( i , j ) is maintenance cost (baht/tonne) for sugarcane remaining on the farms.
The harvesting, transport and maintenance costs obtained from the OCSB are shown in Table 6 for each region for each crop year. It can be seen that the harvesting costs are lower in all regions for the years 2017/18 and 2018/19 for both fresh and fired cane, whereas there is no clear pattern for the transport and maintenance costs.

4.3. Total Revenue and Total Profit

For each type i and region j, the total revenue for a given year is
R ( u ( i , j ) ) = k = 0 N 1 P ( k , i , j ) u k ( i , j ) ,
and the total profit is
J ( u ( i , j ) ) = k = 0 N 1 P t ( k , i , j ) C h ( i , j ) C t ( i , j ) u k ( i , j ) k = 0 N 1 C m ( i , j ) x k ( i , j ) .

4.4. Cutting Constraints

We assume that there are constraints u m a x ( i , j ) on the amount of sugarcane of type i that can be cut each day in a given region j. These constraints could be due, for example, to the number of machines available for cutting fresh sugarcane or the number of workers avaiable for cutting fired sugarcane. In addition, there can be constraints due to the number of trucks available for transporting the cut cane to the mill and the cutting capacities of the mills. We also add the constraints that 0 u k ( i , j ) x k ( i , j ) and that all sugarcane must be cut, that is, x N ( i , j ) = 0 at the end of the harvesting season.

5. Results

5.1. Profits

Examples of the optimal profits computed from the three optimization and two optimal control methods are shown in Table 7 for fresh sugarcane for the four regions of Thailand for the crop year 2017/18. It can be seen that all optimization and optimal control methods give the same results. A comparison of the theoretical results with the actual profit show that the best agreement is for mcf = 1, that is, full maintenance cost with an upper limit of approximately half of the total amount available cut each 15-day period. We have obtained similar results for fresh and fired sugarcane for the four regions of Thailand for the crop years 2012/13, 2013/14, 2014/15 and 2017/18.

5.2. Optimal Daily Cutting Patterns

In this section, we show the optimal daily cutting patterns that we obtained from discrete optimal control. We found that continuous optimal control gave the same optimal cutting patterns. Also, for the special case of zero growth, we found that linear programming gave the same optimal daily cutting patterns. For comparison, optimal cutting patterns computed from the bi-objective and quasi-Newton methods have been published in [14] for 15-day cutting periods for the crop years 2012/13, 2013/14 and 2014/15.
The results in Figure 1 and Figure 2 show the optimal cutting patterns for the crop year 2017/18 if prices, costs, and profits and an upper bound on the amount cut per day are the main factors involved in determining the optimal profits of the farmers. We have computed optimal daily cutting patterns for the crop years 2012/13, 2013/14, 2014/15, 2018/19 and obtained similar results. For each type, region and year, we examined the effect of changing the values of the upper bound on the maximum cutting per day u m a x (tonnes per day) and the effect of reducing the maintenance costs by a factor mcf ( 0 m c f 1 ) of the actual maintenance cost. It can be seen that the optimal cutting patterns are sensitive to changes in the maintenance cost as they range from cutting all sugarcane as early as possible for mcf = 1 to cutting all sugarcane as late as possible for mcf = 0.
Some examples of actual cutting patterns obtained from OCSB data [12,13] are shown in Figure 3 and Figure 4 for comparison with the theoretical optimal cutting patterns. It can be seen that the actual cutting patterns for both fresh and fired sugarcane in 2017/18 are similar to the cutting patterns for mcf = 0.01 in Figure 1 and Figure 2 with a maximum cutting of between 0.15–0.2 total available sugarcane in each region. As far as the authors are aware, the actual cutting patterns are due to constraints on cutting on the farms and to the capacities of the sugar mills.

5.3. Linear Progamming Dual Variables and Hamiltonian Derivatives

In the optimal control algorithms in Section 3, we used the maximum values of the Hamiltonian derivatives H ( t k , x k , u k , λ k ) u k to find the optimal cutting patterns. It is also well known in linear programming that nonzero values of dual variables λ correspond to active constraints and that the values of dual variables correspond to the extra profit (baht) if a cutting constraint can be increased by one unit (tonne). It is therefore expected that maximum Hamiltonian derivative values should correspond to maximum λ values. As an example, Figure 5 shows a plot of the Hamiltonian derivative values for fresh sugarcane for the North-east region in 2017/18 and Figure 6 shows a similar plot for the dual variables from linear programming.

5.4. Computation Times

As we have seen, the three optimization methods and the two optimal control methods all give similar optimal cutting patterns and profits. However, there is a big difference in computation times between the different methods. As shown in Table 8, the linear programming, discrete and continuous optimal control methods all give very fast computation times in the order discrete optimal control, linear programming, continuous optimal control. However, the optimal control methods have the advantage that the objective functions and constraints can be nonlinear. As shown in Table 9, the linear programming again gives a fast computation time whereas the bi-objective and quasi-Newton methods are much slower. However, both the bi-objective and quasi-Newton methods again have the advantage that the objective functions and constaints can be nonlinear.

6. Discussion and Conclusions

We have compared five different optimization and optimal control methods for optimization of sugarcane harvesting in the four sugarcane growing areas of Thailand, namely, North, Central, East and North-east regions for fresh and fired sugarcane for the crop years 2012/13, 2013/14, 2014/15, 2017/18 and 2018/19. To build the mathematical models, we have used price and cost data from the Office of the Cane and Sugar Board [9,10,11,12,13].
For the bi-objective and quasi-Newton methods, we have assumed that there is no growth in sugarcane during the crop year and have divided the crop year into 15-day periods. For the discrete and continuous optimal control methods, we have allowed for the possibility of growth during a crop year through a logistic function and have divided the crop year into 1-day periods. We have also used linear programming to find optimal cutting patterns for both 1-day and 15-day periods for a no-growth model.
As shown in Table 7, the profits obtained from the five different optimization methods have been the same to at least 4-digit accuracy for both fresh and fired for all regions for the crop year 2017/18. We have obtained similar results for both types of sugarcane for all four regions of Thailand for the crop years 2012/13, 2013/14, 2014/15, 2017/18 and 2018/19. It can also be seen from the table that the actual profits are within the range of computed profits for mcf = 1 (full maintenance cost). It is also clear from the results that changes in the maintenance cost are one of the main factors affecting the optimal cutting patterns, with high maintenance costs suggesting cutting as early as possible and low maintenance costs suggesting cutting later.
Although the bi-objective, quasi-Newton, linear programming and optimal control methods give the same optimal cutting patterns and profits for the case of zero growth, there is a big difference in computation times as shown in Table 8 and Table 9 with the linear programming and optimal control methods being orders of magnitude faster than the bi-objective and quasi-Newton methods. Further, for the bi-objective, quasi-Newton methods, we have considered 15-day cutting periods for a total of approximately 12 decision variables, whereas for the linear programming and optimal control methods we have considered daily cutting periods for a total of approximately 180 decision variables. However, linear programming can only be used if the growth is zero and the profit function is a linear function of the cutting patterns and the constraints. Finally, we have found that the programming for the discrete optimal control is much simpler than for the continuous optimal control and that the computation times are also shorter.
There are a number of important factors that have not been considered in this paper. For example, there are likely to be changes in the harvesting and transport costs per tonne as a function of amount cut per day and also changes in the costs of harvesting and transport during the cutting season. Including these effects would make the profit a nonlinear function of control and change the solution from the bang-bang control observed in the models in this paper. Unfortunately, we have not been able to obtain the cost data that would be required to include these effects.

Author Contributions

Conceptualization, W.P., S.S., C.K. and E.J.M.; methodology, W.P., S.S., C.K. and E.J.M.; software, W.P., S.S. and E.J.M.; validation, S.S., C.K. and E.J.M.; formal analysis, W.P., S.S. and E.J.M.; writing—original draft preparation, W.P.; writing—review and editing, S.S., C.K. and E.J.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by King Mongkut’s University of Technology North Bangkok (Contract No. KMUTNB-60-GOV-071) and Centre of Excellence in Mathematics, the Commission on Higher Education, Thailand.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data on sugarcane in this paper was obtained from the Ministry of Industry and the Ministry of Agriculture and Co-operatives of the Royal Thai government. These data can be found here: [http://www.ocsb.go.th] (accessed on 18 March 2021).

Acknowledgments

This research was partially supported by the Graduate College, King Mongkut’s University of Technology North Bangkok. The authors would like to express their thanks to the anonymous referees for their time and helpful comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pontryagin, L.S. Mathematical Theory of Optimal Processes; Routledge: Abingdon-on-Thames, UK, 2018. [Google Scholar]
  2. Luenberger, D.G. Introduction to Dynamic Systems: Theory, Models & Applications; John Wiley & Sons: New York, NY, USA, 1979. [Google Scholar]
  3. Lenhart, S.; Workman, J.T. Optimal Control Applied to Biological Models; Chapman & Hall/CRC: Boca Raton, FL, USA, 2007. [Google Scholar]
  4. Deb, K. Multi-Objective Optimization Using Evolutionary Algorithms; John Wiley & Sons: New York, NY, USA, 2005. [Google Scholar]
  5. Gill, P.E.; Murray, W.; Wright, M.H. Practical Optimization; Wiley-Interscience: Chichester, UK, 1981. [Google Scholar]
  6. Rao, S.S. Engineering Optimization: Theory and Practice, 4th ed.; Wiley-Interscience: New York, NY, USA, 2009. [Google Scholar]
  7. Fletcher, R. Practical Methods of Optimization, 2nd ed.; John Wiley & Sons: New York, NY, USA, 2013. [Google Scholar]
  8. Chvátal, V. Linear Programming; W.H. Freeman: San Francisco, CA, USA, 1983. [Google Scholar]
  9. The Cane and Sugar Board, Ministry of Industry of Thailand. Data of Sugarcane. 2013. Available online: http://www.ocsb.go.th (accessed on 18 March 2021).
  10. The Cane and Sugar Board, Ministry of Industry of Thailand. Data of Sugarcane. 2014. Available online: http://www.ocsb.go.th (accessed on 18 March 2021).
  11. The Cane and Sugar Board, Ministry of Industry of Thailand. Data of Sugarcane. 2015. Available online: http://www.ocsb.go.th (accessed on 18 March 2021).
  12. The Cane and Sugar Board, Ministry of Industry of Thailand. Data of Sugarcane. 2018. Available online: http://www.ocsb.go.th (accessed on 18 March 2021).
  13. The Cane and Sugar Board, Ministry of Industry of Thailand. Data of Sugarcane. 2019. Available online: http://www.ocsb.go.th (accessed on 18 March 2021).
  14. Pornprakun, W.; Sungnul, S.; Kiataramkul, C.; Moore, E.J. Determining optimal policies for sugarcane harvesting in Thailand using bi-objective and quasi-Newton optimization methods. Adv. Differ. Equations 2019, 2019, 1–15. [Google Scholar] [CrossRef]
  15. Kiataramkul, C.; Wake, G.; Ben-Tal, A.; Lenbury, Y. Optimal nutritional intake for fetal growth. Math. Biosci. Eng. 2011, 8, 723–732. [Google Scholar] [PubMed]
  16. Dine, W.; Lenhart, S.; Behncke, H. Discrete time optimal harvesting of fish populations with age structure. Lett. Biomath. 2014, 1, 193–207. [Google Scholar] [CrossRef]
  17. Kiataramkul, C.; Matkhao, P. Optimal Control Problem of Food Intake of Swine During Post Weaning Period. Kasetsart J. (Nat. Sci.) 2015, 49, 983–989. [Google Scholar]
  18. Puengpao, S.; Kiataramkul, C.; Moore, E.J. Application of discrete optimal control to feeding of farm animals. In Proceedings of the International Conference in Mathematics and Applications (ICMA-MU2016), Bangkok, Thailand, 17–19 December 2016; pp. 151–160. [Google Scholar]
  19. Gomes, F.R.A. Bi-objective mathematical model for choosing sugarcane varieties with harvest residual biomass in energy cogeneration. Int. J. Agric. Biol. Eng. 2012, 5, 50–58. [Google Scholar]
  20. Sungnul, S.; Pornprakun, W.; Prasattong, S.; Baitiang, C. Optimal Time of Sugarcane Harvesting for Sugar Factories in Thailand. In Proceedings of the International Conference in Mathematics and Applications (ICMA-MU2016), Bangkok, Thailand, 17–19 December 2016; pp. 185–194. [Google Scholar]
  21. Sungnul, S.; Pornprakun, W.; Prasattong, S.; Baitiang, C. Multi-Objective Mathematical Model for the Optimal Time to Harvest Sugarcane. Appl. Math. 2017, 8, 329. [Google Scholar] [CrossRef] [Green Version]
  22. Naranong, V.; Tansakun, S.; Chantaprom, T. A Study of Reform of the Structure of Thailand’s Sugar and Cane Industry. Thail. Dev. Res. Inst. 2013. Available online: https://dric.nrct.go.th/index.php?/Search/SearchDetail/278291 (accessed on 18 March 2021).
Figure 1. Optimum cutting of fresh sugarcane in four regions of Thailand in year 2017/18 showing effects of changing maintenance cost by fraction mcf of actual cost. North: (ac), Central: (df), East: (gi), North-east: (jl).
Figure 1. Optimum cutting of fresh sugarcane in four regions of Thailand in year 2017/18 showing effects of changing maintenance cost by fraction mcf of actual cost. North: (ac), Central: (df), East: (gi), North-east: (jl).
Computation 09 00036 g001
Figure 2. Optimum cutting of fired sugarcane in four regions of Thailand in year 2017/18 showing effects of changing maintenance cost by fraction mcf of actual cost. North: (ac), Central: (df), East: (gi), North-east: (jl).
Figure 2. Optimum cutting of fired sugarcane in four regions of Thailand in year 2017/18 showing effects of changing maintenance cost by fraction mcf of actual cost. North: (ac), Central: (df), East: (gi), North-east: (jl).
Computation 09 00036 g002
Figure 3. Actual cutting patterns for fresh and fired sugarcane for 2017/18.
Figure 3. Actual cutting patterns for fresh and fired sugarcane for 2017/18.
Computation 09 00036 g003
Figure 4. Actual cutting patterns for fresh and fired sugarcane for 2018/19.
Figure 4. Actual cutting patterns for fresh and fired sugarcane for 2018/19.
Computation 09 00036 g004
Figure 5. Example of H ( t k , x k , u k , λ k ) u k for fresh sugarcane in North-east region in 2017/18. Compare optimal cutting plots in Figure 1d–f.
Figure 5. Example of H ( t k , x k , u k , λ k ) u k for fresh sugarcane in North-east region in 2017/18. Compare optimal cutting plots in Figure 1d–f.
Computation 09 00036 g005
Figure 6. Example of dual variables from linear programming for fresh sugarcane in North-east region in 2017/18. Compare Hamiltonian derivative plots in Figure 5.
Figure 6. Example of dual variables from linear programming for fresh sugarcane in North-east region in 2017/18. Compare Hamiltonian derivative plots in Figure 5.
Computation 09 00036 g006
Table 1. The basic prices (Baht/tonne) set by the Royal Thai Government and the total amounts ( 10 6 tonnes) of fresh and fired sugarcane delivered to the factories in each region in each crop year.
Table 1. The basic prices (Baht/tonne) set by the Royal Thai Government and the total amounts ( 10 6 tonnes) of fresh and fired sugarcane delivered to the factories in each region in each crop year.
TypeYearBasic Price (baht)Total Amount (tonnes)
P w P c NorthCentralEastNorth-East
Fresh2012/13198.47999.207.0369.3090.95516.909
2013/14194.67958.317.0769.6001.23520.013
2014/15197.41900.007.3849.2351.26419.028
2017/18141.00880.0011.39911.7351.74120.623
2018/19141.00700.0013.94411.0691.70524.221
Fired2012/13140.00999.2017.56221.1893.73123.311
2013/14140.00958.3117.13420.4793.23724.892
2014/15140.00900.0017.99918.7663.94628.338
2017/18140.00880.0021.64124.8584.94237.989
2018/19140.00700.0019.16320.4684.17436.226
Table 2. The Commercial Cane Sugar System (CCS) value in each region for crop year 2017/18.
Table 2. The Commercial Cane Sugar System (CCS) value in each region for crop year 2017/18.
PeriodNorthCentralEastNorth-East
(1) 1–15 December 201710.3879.93210.48211.655
(2) 16–31 December 201710.72610.26210.80712.009
(3) 1–15 January 201811.03210.48510.99412.313
(4) 16–31 January 201811.37810.85611.36712.596
(5) 1–14 February 201811.65411.13411.67012.840
(6) 15–28 February 201811.85611.32811.91113.021
(7) 1–15 March 201812.00511.46312.07513.174
(8) 16–31 March 201812.12011.58612.21113.291
(9) 1–15 April 201812.14711.60312.25313.333
(10) 16–30 April 201812.13911.59512.26513.337
(11) 1–15 May 201812.11211.57212.23513.311
(12) 16–31 May 201812.10111.56412.18913.298
(13) 1–15 June 201812.10111.56112.18713.298
Table 3. The Commercial Cane Sugar System (CCS) value in each region for crop year 2018/19.
Table 3. The Commercial Cane Sugar System (CCS) value in each region for crop year 2018/19.
PeriodNorthCentralEastNorth-East
(1) 16–30 November 20189.7609.752-11.530
(2) 1–15 December 201810.17510.11411.31611.909
(3) 16–31 December 201810.48610.32211.64112.163
(4) 1–15 January 201910.82210.62211.85512.507
(5) 16–31 January 201911.18210.96112.17612.819
(6) 1–14 February 201911.44411.21312.40413.036
(7) 15–28 February 201911.63911.41912.59913.179
(8) 1–15 March 201911.83411.59612.79013.311
(10) 16–31 March 201911.98611.65412.85813.427
(11) 1–15 April 201912.01511.65812.85813.464
(12) 16–30 April 201912.01511.658-13.481
(13) 1–15 May 201912.01511.658-13.480
Note: The-means that no sugarcane was delivered to the mills in that period.
Table 4. The price (baht/tonne) of fresh and fired sugarcane in each region for crop year 2017/18.
Table 4. The price (baht/tonne) of fresh and fired sugarcane in each region for crop year 2017/18.
PeriodNorthCentralEastNorth-East
FreshFiredFreshFiredFreshFiredFreshFired
(1) 1–15 December 20171041.45989.761017.42966.541046.47994.601108.381054.40
(2) 16–31 December 20171059.321007.021034.82983.351063.591011.141127.101072.48
(3) 1–15 January 20181075.491022.631046.63994.761073.471020.681143.121087.96
(4) 16–31 January 20181093.771040.291066.191013.651093.161039.701158.091102.42
(5) 1–14 February 20181108.351054.371080.861027.821109.151055.151170.951114.83
(6) 15–28 February 20181118.991064.651091.121037.731121.901067.461180.481124.05
(7) 1–15 March 20181126.881072.271098.261044.621130.551075.821188.591131.87
(8) 16–31 March 20181132.941078.121104.741050.881137.741082.761194.751137.83
(9) 1–15 April 20181134.361079.491105.621051.741139.971084.911197.001140.00
(10) 16–30 April 20181133.941079.091105.231051.361140.581085.501197.201140.19
(11) 1–15 May 20181132.521077.721104.021050.191139.021084.001195.821138.86
(12) 16–31 May 20181131.961077.171103.561049.751136.591081.651195.141138.21
(13) 1–15 June 2018--1103.431049.621136.461081.531195.14-
Note: The-means that no sugarcane was delivered to the mills in that period.
Table 5. The price (baht/tonne) of fresh and fired sugarcane in each region for crop year 2018/19.
Table 5. The price (baht/tonne) of fresh and fired sugarcane in each region for crop year 2018/19.
PeriodNorthCentralEastNorth-East
FreshFiredFreshFiredFreshFiredFreshFired
(1) 16–30 November 2018830.94780.37830.57780.02--905.28851.52
(2) 1–15 December 2018848.33797.02845.77794.57896.29842.92921.17866.73
(3) 16–31 December 2018861.41809.54854.53802.95909.92855.97931.83876.94
(4) 1–15 January 2019875.51823.03867.13815.01918.92864.58946.27890.76
(5) 16–31 January 2019890.66837.53881.37828.64932.40877.48959.39903.32
(6) 1–14 February 2019901.66848.06891.93838.75941.96886.63968.51912.05
(7) 15–28 February 2019909.83855.88900.58847.03950.17894.49974.52917.80
(8) 1–15 March 2019918.04863.74908.04854.16958.17902.15980.06923.10
(10) 16–31 March 2019924.40869.82910.46856.48961.05904.91984.95927.78
(11) 1–15 April 2019925.61870.99910.66856.67961.05904.90986.50929.26
(12) 16–30 April 2019--910.66856.67--987.18929.92
(13) 1–15 May 2019------987.18929.91
Note: The-means that no sugarcane was delivered to the mills in that period.
Table 6. The average costs of sugarcane production (baht/tonne) in each region for each crop year.
Table 6. The average costs of sugarcane production (baht/tonne) in each region for each crop year.
YearRegionHarvesting Cost C h ( i , j )
(baht/tonne)
Transport Cost C t ( i , j )
(baht/tonne)
Maintenance Cost C m ( i , j )
(baht/tonne)
2012/13North924.28136.5152.70
Central872.95147.6694.24
East836.25148.8992.51
North-east765.18141.9756.30
2013/14North815.89149.4666.07
Central781.41147.0486.59
East912.61165.66106.01
North-east875.84151.6459.55
2014/15North1061.42182.3679.38
Central954.16155.3082.32
East1024.71194.8994.30
North-east987.00140.1291.71
2017/18North771.27143.6284.94
Central771.91137.2284.83
East848.51146.21102.43
North-east871.78150.7466.80
2018/19North738.86138.40104.87
Central754.86138.2483.26
East768.63157.6875.48
North-east775.09149.5859.69
Table 7. The total profit ( × 10 10 Baht) of fresh sugarcane in each region in crop year 2017/18.
Table 7. The total profit ( × 10 10 Baht) of fresh sugarcane in each region in crop year 2017/18.
RegionmcfP.O.C.Actual ProfitTotal Profit ( 10 10 Baht)
Bi-objLin prog.qnDisc.Cont.
North00.20.117760.928080.928080.928080.928060.92806
0.40.929060.929060.929060.929060.92906
0.60.929290.929290.929290.929320.92932
0.80.929350.929350.929350.929390.92939
1.00.929400.929400.929400.929420.92942
0.250.20.386050.386090.386090.385980.38598
0.40.596450.596620.596620.605160.60516
0.60.666480.666570.666570.678160.67816
0.80.701630.701670.701670.714620.71462
1.00.736770.736770.736770.736520.73652
10.2−1.06910−1.06910−1.06910−1.06922−1.06920
0.4−0.16008−0.16008−0.16008−0.12257−0.12257
0.60.142700.142700.142700.193090.19309
0.80.294210.294210.294210.350620.35062
1.00.445730.445730.445730.445480.44548
Central00.20.096940.933540.933540.933530.933530.93353
0.40.934470.934470.934450.934450.93445
0.60.934710.934710.934660.934660.93466
0.80.934750.934750.934720.934720.93472
1.00.934800.934800.934750.934750.93475
0.250.20.374820.374820.374630.374630.37463
0.40.594080.594180.602800.602800.60280
0.60.667860.668020.680110.680110.68011
0.80.705060.705350.718780.718780.71878
1.00.742690.742690.741860.741860.74186
10.2−1.13210−1.13210−1.13228−1.13228−1.13230
0.4−0.18941−0.18941−0.15079−0.15079−0.15079
0.60.125530.125530.177800.177800.17780
0.80.283420.283420.341850.341850.34185
1.00.441310.441310.440480.440480.44048
East00.20.000560.149250.149250.149250.149240.14924
0.40.149460.149460.149460.149460.14946
0.60.149500.149500.149500.149510.14951
0.80.149510.149510.149510.149520.14952
1.00.149520.149520.149520.149530.14953
0.250.20.051080.051080.051080.051050.05105
0.40.090440.090460.090460.092010.09201
0.60.103760.103790.103790.105960.10596
0.80.110470.110520.110520.112940.11294
1.00.117250.117250.117250.117100.11710
10.2−0.21515−0.21515−0.21515−0.21518−0.21518
0.4−0.04798−0.04798−0.04798−0.04113−0.04113
0.60.007940.007940.007940.017220.01722
0.80.035980.035980.035980.046350.04635
1.00.064010.064010.064010.063860.06386
North-east00.20.144161.812211.812201.812201.812101.81210
0.41.813601.813601.813601.813701.81370
0.61.814061.814101.814101.814201.81420
0.81.814351.814401.814401.814401.81440
1.01.814641.814601.814601.814501.81450
0.250.20.719140.719130.719130.718920.71892
0.41.160301.160501.160501.178551.17850
0.61.305901.306201.306201.330301.33030
0.81.378131.378701.378701.405661.40570
1.01.451241.451201.451201.450901.45090
10.2−2.26570−2.26570−2.26570−2.26595−2.26600
0.4−0.39165−0.39165−0.39165−0.31416−0.31416
0.60.231660.231660.231660.335340.33534
0.80.542960.542960.542960.659040.65904
1.00.854270.854270.854270.853930.85393
Note: P.O.C. means the percentage of total available sugarcane that can be cut in a 15-day period.
Table 8. CPU time (seconds) of fresh and fired sugarcane cutting in each day in crop year 2017/18 and 2018/19.
Table 8. CPU time (seconds) of fresh and fired sugarcane cutting in each day in crop year 2017/18 and 2018/19.
Crop YearMethodmcfType of Sugarcane
FreshFired
2017/18discrete00.0110.017
continuous0.2130.152
linear programming0.0370.041
discrete10.0100.013
continuous0.1670.157
linear programming0.0410.037
2018/19discrete00.0100.010
continuous0.1710.162
linear programming0.0310.034
discrete10.0110.010
continuous0.1400.133
linear programming0.0360.036
Table 9. CPU time (seconds) of fresh and fired sugarcane cutting in each 15-day period in crop year 2017/18 and 2018/19.
Table 9. CPU time (seconds) of fresh and fired sugarcane cutting in each 15-day period in crop year 2017/18 and 2018/19.
Crop YearMethodmcfType of Sugarcane
FreshFired
2017/18bi-objective00.0240.031
linear programming0.0190.016
quasi-Newton83.74372.457
bi-objective15.2265.385
linear programming0.0380.031
quasi-Newton12.90515.426
2018/19bi-objective00.0230.023
linear programming0.0400.035
quasi-Newton98.37483.416
bi-objective15.0104.926
linear programming0.0330.036
quasi-Newton13.04215.668
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pornprakun, W.; Sungnul, S.; Kiataramkul, C.; Moore, E.J. An Application of Optimal Control to Sugarcane Harvesting in Thailand. Computation 2021, 9, 36. https://0-doi-org.brum.beds.ac.uk/10.3390/computation9030036

AMA Style

Pornprakun W, Sungnul S, Kiataramkul C, Moore EJ. An Application of Optimal Control to Sugarcane Harvesting in Thailand. Computation. 2021; 9(3):36. https://0-doi-org.brum.beds.ac.uk/10.3390/computation9030036

Chicago/Turabian Style

Pornprakun, Wisanlaya, Surattana Sungnul, Chanakarn Kiataramkul, and Elvin James Moore. 2021. "An Application of Optimal Control to Sugarcane Harvesting in Thailand" Computation 9, no. 3: 36. https://0-doi-org.brum.beds.ac.uk/10.3390/computation9030036

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop