Allocation To Minimize Fuel Consumptio
Allocation To Minimize Fuel Consumptio
218
2019 IFAC MMM
D.M. Bajany et al. / IFAC PapersOnLine 52-14 (2019) 207–212
Stellenbosch, South Africa, August 28-30, 2019 209
and S is the number of shovels. For a complete shift, the idle time Ij of the j-th shovel is
calculated as:
N
U
The loading time tyl,j of a truck y with the j-th shovel is Ij = sh − y
Zlo,jβ tyl,j , (11)
determined by: y=1 β=1
Cy
tyl,j = , (5) where sh is the shift duration, and, the fuel Fidle,j con-
Cj sumed by the j-th shovel during its idling period is given
where Cy is the loading capacity of the truck y and Cj is as:
the capacity of the shovel j. Fidle,j = Ij fj,idle , (12)
The fuel consumption of a dump truck y per single cycle with fj,idle , the fuel consumption rate of shovel j during
is given as: idling time. For a complete shift, the fuel Fu,j consumed
by the j-th shovel during periods when it was used to load
y y
trucks is calculated as:
Fijβ = fe,ij tye,ij + flo,jβ
y
tylo,jβ N U
y
y
+ fidle (tyl,j + tyu + tyw,j + tyw,β ), (6) Fu,j = fj Zlo,jβ tyl,j , (13)
y=1 β=1
y where fj is the fuel consumption of shovel j during the
where fe,ijis the fuel consumption of the empty truck y
moving from the unloading point i to the shovel j, flo,jβ y time it is used to load a truck.
is the fuel consumption of the loaded truck y moving from The total fuel Fj consumed by the j-th shovel during a
the shovel j to the unloading point β, fidle y
is the fuel shift is determined by:
consumption of the truck y during engine idling time, tyu Fj = Fidle,j + Fu,j . (14)
is the unloading time of the truck y, tyw,j is the waiting The total fuel F s consumed by all the shovels for the whole
time of the truck y at the shovel j, tyw,β is the waiting time shift duration is given as:
S
of the truck y at the dumping point β. Note that i = β
implies that the truck returns to the original unloading Fs = Fj . (15)
point where it started its cycle. j=1
The number of times xyj a y truck is loaded by the j-th The total fuel F consumed for whole haulage operations
shovel during a shift and the number of times a y truck during a shift is calculated as:
has dumped his load at the dumping β during a shift are F = F t + F s. (16)
calculated by equations (7) and (8) respectively: There are many variables that can be used to evaluate
U
U
the haul truck energy efficiency (Krzyzanowska (2007)).
xyj = y
Ze,ij or xyj = y
Zlo,jβ , (7) In this study, the liters of fuel per tonne moved is adopted
i=1 β=1 as a performance indicator for the truck-shovel dispatching
S S
system. The liters of fuel per tonne moved at the end of a
shift is evaluated by:
xyβ = y
Ze,ij or xyβ = y
Zlo,jβ , (8)
F
j=1 j=1 LT ts = N U S . (17)
with U , the number of unloading points, and S, the y
Zlo,jβ Cy
number total of shovels.
y=1 β=1 j=1
For xyj cycles that a y truck has accomplished to the j-
As previously mentioned, the objective of this study is to
th shovel during a shift, the fuel consumption F tyj of that
minimize the fuel consumption of both trucks and shovels.
truck is calculated as: The model is built in such a way that all the requirements
U
U
of dump sites are met and the optimal number of trips
F tyj = y
Ze,ij y
fe,ij tye,ij + y
Zlo,jβ y
flo,jβ tylo,jβ of each truck on each route of the pit is determined. The
i=1 β=1 technical specifications, such as payload of trucks, loaded
U
capacity of shovels, and fuel consumption in function of the
y y
+ Zlo,jβ fidle tyl,j + tyu + tyw,j + tyw,β , (9) operating conditions of equipment are directly considered
β=1 in the mathematical model. The objective function is given
y in the following:
where Ze,ij is the number of journeys that the truck y has
min F. (18)
travelled from the i-th unloading point to the j-th shovel
y The operating constraints of the problem include the fol-
during a shift and Zlo,jβ is the number of journeys that
lowing:
the truck y has travelled from the j-th shovel to the β-th - The material transported from all shovels should be
unloading point during a shift. greater than the handling demand of each unloading point
From the above equation, for a given shift, the total (of ore and waste). This constraint ensures that trucks are
amount of fuel F t consumed by all trucks is calculated dispatched so that the production target at the mine is
as: satisfied. The way this requirement is taken into account
S N
Ft = F tyj , (10) is shown in equation (19):
j=1 y=1 N
S
y
with N , the number of trucks used in the haulage opera- Zlo,jβ Cy ≥ Dβ sh, ∀β ∈ {i, w}, (19)
tions. y=1 j=1
219
2019 IFAC MMM
210
Stellenbosch, D.M. Bajany et al. / IFAC PapersOnLine 52-14 (2019) 207–212
South Africa, August 28-30, 2019
with Dβ , the handling demand per hour of the unloading times that a truck travels empty from that unloading point
point i and sh, the shift duration. to different shovels can not exceed one. This constraint is
- Material transported by trucks during the shift duration written as follows:
from each loading point is less than the shovel capacity S S
y y
allocated to that pit. Ze,ij − Zlo,jβ ≤ 1 f or i = β, ∀j, ∀y. (27)
j=1 j=1
N
U
y - The difference between the number of times an empty
Zlo,jβ Cy ≤ Cj sh, ∀j. (20)
truck y travels from an unloading point i to different
y=1 β=1
shovels and the number of time this truck dumps its load at
- The ore transported by trucks during the shift duration that unloading point can not exceed one. This constraint
from a loading point j to the unloading points of ore is is written as follows:
less or equal to their available quantity at the considered S S
loading point. y
Zlo,jβ − y
Ze,ij ≤ 1 f or i = β, ∀j, ∀y. (28)
N
y j=1 j=1
Zlo,jβ C y ≤ Oj , (21)
y=1
- The last constraint (29) ensures that the number of trips
of trucks are positive and integer.
∀j, ∀sh, β ∈ i with i = {i1 , i2 , ..., iq }, with Oj , the y y
available quantity of ore at the loading point j and q, the Ze,ij ∈ N + ; Zlo,jβ ∈ N + ; i = 1, 2, ..., U ; β = 1, 2, ..., U ;
total number of ore unloading points. j = 1, 2, ..., S; y = 1, 2, ..., N. (29)
- The waste transported by trucks during the shift duration
Constraints (26), (27) and (28) ensure that the continuity
from a loading point j to the unloading points of waste is
of loading and transportation is maintained.
less or equal to their available quantity at the considered
loading point.
N 5. SIMULATIONS
y
Zlo,jβ Cy ≤ Wj , (22)
y=1 A case study of a hypothetical downgrade open-pit mine
with two unloading points (one for waste materials and
∀j, ∀sh, β ∈ w with w = {w1 , w2 , ..., wr }, with Wj , the
one for ores) and three shovels located at three pits is
available quantity of waste at the loading point j and r,
considered for optimization and simulation of the model.
the total number of waste unloading points.
Among these shovels, two have a capacity of 1600 t/h and
- For the sh-th shift, the utilization time of each shovel
one, a capacity of 2000 t/h. The length of each road in the
is less or equal to the shift duration. This constraint is
haulage network is given in Table 1. The total resistance
written as follows:
of the haul road is assumed to be equal to 4%. Operating
N
speeds of loaded trucks are estimated using the Rimpull-
xyj tyl,j ≤ sh, ∀j. (23) Speed-Grade curves given in Caterpillar.Inc (2007) 1 . The
y=1
mine topography and operating parameters of the loading
and haulage equipment used as input parameters of the
From equation (7), the constraint (23) can also be written
model are shown in Table 1. The three possible shovel
as:
allocations in different pits are determined by equation
3!
N
U (2) as Pj = 2!1! = 3, they are displayed on Fig. 1.
y
Zlo,jβ tyl,j ≤ sh, ∀j. (24)
y=1 β=1 6. SIMULATION RESULT
- The sum of cycle times of all cycles performed by a truck
y during the shift is less or equal to the shift duration. In order to show the effectiveness of the developed method
and the impact of the waiting time on the fuel consumed
S
U S
U
in the haulage operations, three dispatch scenarios with
different waiting time at loading and unloading points have
F tyj = y
Ze,ij y
fe,ij tye,ij + y
Zlo,jβ y
flo,jβ tylo,jβ been considered in the simulations. The waiting times con-
j=1 i=1 j=1 β=1
sidered at loading and dumping points were as follows: 0.30
S
U min for the first scenario, 0.48 min for the second scenario
y y
+ Zlo,jβ fidle tyl,j + tyu + tyw,j + tyw,β ≤ sh, and 0.78 min for the third scenario. The shift duration of
j=1 β=1 each scenario was 8 hours. Solving the developed model
∀y. (25) with a mixed integer nonlinear programming algorithm
- The number of trips that a truck y makes to a shovel for the three scenarios, the following results were found.
equals to the number of trips that the same truck leaves For the same amount of material transported (37248 t),
that shovel. Equation (26) shows how this requirement is the first scenario, with the shovels allocate following the
considered: allocation (a), results in 130.99 milliliters per tonne and
4879.80 liters of fuel consumed for the whole haulage
U
U
operation; with the allocation (b), it results in 136,90
y y milliliters per tonne and 5099.30 liters of fuel consumed for
Ze,ij = Zlo,jβ , ∀j, ∀y. (26)
i=1 β=1 the whole haulage operation and with the allocation (c),
- The difference between the number of times a truck y its results in 132,92 milliliters per tonne and 4851,12 liters
dumps its load at an unloading point i and the number of 1 Caterpilar. Inc, 777D Off-Highway Truck, USA.2007.
220
2019 IFAC MMM
D.M. Bajany et al. / IFAC PapersOnLine 52-14 (2019) 207–212
Stellenbosch, South Africa, August 28-30, 2019 211
Chart Title
Table 1. Mine topography and resources 3500
3000
Litres
1500
Rated payload 96 t
1000
Gross vehicle weight 161 030 kg
500
Speed of empty truks 50 km/h
0
Speed of loaded trucks 36 km/h Allocation (a) Allocation (b) Allocation (c)
Scenario 1: 0.30 min 2325,13 2554,64 2530,2
Fuel consumption (idle time) 22.38 L/h
Scenario 2: 0.48 min 2498,4 2741,5 2682,35
Fuel consumption empty truck 44.76 L/h Senario 3: 0.78 min 2525,23 2945,33 2792,16
Fuel consumption loaded truck 78.33 L/h
Shovel capacity 1600 t/h
2000 t/h
Fig. 2. Fuel consumed in the haulage by trucks and shovels
per scenario Chart Title
Fuel consumption of shovel
- during idle time 1600 t/h - 6.6 L/h 5600
2000 t/h -9.5 L/h 5400
- during working time 1600 t/h - 117 L/h 5200
Liters
2000 t/h -130 L/h 5000
Mine topography Downgrade mine 4800
Fuel consumption loaded truck 78.33 L/h 4600
Gradient 1:14 40 4400
Distance between dumping Allocation (a) Allocation (b) Allocation (c)
and loading point i1 − j1 1.5 km Scenario 1: 0.30 min 4879,8 5099,3 4951,12
Scenario 2: 0.48 min 4956,8 5203,37 5052,16
i1 − j2 2.5 km
Scenario 3: 0.78 min 5259,14 5520,5 5338,45
i1 − j3 3 km
i2 − j1 3.5 km
i2 − j2 2.5 km Fig. 3. FuelChart
consumed
Title by trucks Chart Title
i2 − j3 1.5 km
Demand of ore at unloading
point i1 2375 t/h 32% 30%
35%
44%
of fuel consumed for the whole haulage operation. The 24% 35%
second scenario, with the allocation (a), results in 133,08
Shovel j1=1600t/h
milliliters per tonne and 4956 liters of fuel consumed for Shovel j2=1600t/h Chart Title
Shovel j1=1600t/h
Shovel j2=2000t/h
the whole haulage operation; with the allocation (b), it Allocation (a)
Shovel j3=2000t/h Allocation (b) Shovel j3=1600t/h
milliliters per tonne and 5259,14 liters of fuel consumed Shovel j1=2000t/h
for the whole haulage operation; with the allocation (b), Shovel j2=1600t/h
Allocation(c)
it results in 148,21 milliliters per tonne and 5520.50 liters Shovel j3=1600t/h
tonne.
34%
The applicability of the developed model to determine 43%
the best shovel allocation was also demonstrated in this
study. Indeed, from the obtained results it is clear that 23%
the second allocation is the least efficient one compared to Shovel j1=2000t/h
the other two. Indeed, for the first scenario, by allocating Shovel j2=1600t/h
the shovel with the highest capacity at the second loading Shovel j3=1600t/h
221
2019 IFAC MMM
212
Stellenbosch, D.M. Bajany et al. / IFAC PapersOnLine 52-14 (2019) 207–212
South Africa, August 28-30, 2019
Chart Title
Shovel j2=1600t/h Shovel j2=2000t/h Bahrami, M. (2015). Theoretical and experimental
Allocation (a) Shovel j3=2000t/h Shovel j3=1600t/h
investigation into anti-idling a/c system for trucks.
Energy Conversion and Management, 98, 173–183.
34%
Bajany, D., Xia, X., and Zhang, L. (2017a). A milp model
44% for truck-shovel scheduling to minimize fuel consump-
tion. Energy Procedia, 105, 2739–2745.
22%
Bajany, M.D. et al. (2017b). A mixed integer linear pro-
Shovel j1=2000t/h gramming model for truck-shovel scheduling to minimize
Allocation (c)
Shovel j2=1600t/h
Shovel j3=1600t/h
fuel consumption. Master’s thesis, University of Preto-
ria.
Fig. 6. Contribution of each shovel to the total handling Caterpillar.Inc (2007). 777D Off-Highway Truck. USA.
demand (Waiting time: 0.78 min) Da Cunha Rodovalho, E., Lima, H.M., and de Tomi, G.
(2016). New approach for reduction of diesel consump-
of this observation is that, with the allocation (b), the tion by comparing different mining haulage configura-
total travel distances of loaded and empty trucks are much tions. Journal of Environmental Management, 172, 177–
longer than those obtained with the other two allocations. 185.
Hence, trucks consume more fuel with the allocation (b). Jochens, P. (1980). The energy requirements of the mining
This is shown in the Fig. 3. Another explanation of this and metallurgical industry in south africa: presidential
observation is that the shovel of capacity 2000 t/h is not address. Journal of the Southern African Institute of
properly utilized with the allocation (b). This can be seen Mining and Metallurgy, 80(9), 331–343.
by comparing the contributions of this shovel to the total Kecojevic, V. and Komljenovic, D. (2010). Haul truck
handling demand displayed in Fig. 4, Fig. 5 and Fig. 6. To fuel consumption and co 2 emission under various engine
be specific, only 35% of materials moved in the haulage load conditions. Mining Engineering, 62(12), 44–48.
are loaded by the shovel of capacity 2000 t/h; whereas Krzyzanowska, J. (2007). The impact of mixed fleet haul-
for the two other allocations, 44% (allocation (a)) and ing on mining operations at venetia mine. Journal of
39% (allocation (c)) of materials are loaded by this shovel. The Southern African Institute of Mining and Metal-
With the 2000 t/h shovel being the most fuel efficient, less lurgy, 107(4), 215–224.
utilization of it resulted in more fuel consumption for the Lashgari, A., Johnson, C., Kecojevic, V., Lusk, B., and
allocation (b). Hoffman, J. (2014). Nox emission of equipment and
blasting agents in surface mining. In Proceedings of
The results of this work are important as the proposed the 143rd Annual Meeting and Exhibition of the Society
method considers the waiting time of trucks at shovels of Mining, Metallurgy and Exploration, Salt Lake City,
and dumping points and the specific characteristic of each UT.
truck and shovel used in the haulage, unlike previously May, M.A. (2013). Applications of queuing theory for open-
reported optimization work. pit truck/shovel haulage systems. Ph.D. thesis, Virginia
Tech.
7. CONCLUSION Peralta, S., Sasmito, A.P., and Kumral, M. (2016). Re-
liability effect on energy consumption and greenhouse
A model for minimization of fuel consumption of dump gas emissions of mining hauling fleet towards sustainable
trucks and shovels in open-pit mine was developed. To mining. Journal of Sustainable Mining, 15(3), 85–94.
evaluate the developed model, a case study of an open- Rahman, S.A., Masjuki, H., Kalam, M., Abedin, M.,
pit mine with two unloading points and three shovels was Sanjid, A., and Sajjad, H. (2013). Impact of idling on
considered. In the case of a mixed fleet of shovels, for each fuel consumption and exhaust emissions and available
possible shovel allocation, the optimal number of trips idle-reduction technologies for diesel vehicles–a review.
that each truck should make on each route of the mine Energy Conversion and Management, 74, 171–182.
during a shift was realized. It was possible to determine Sahoo, L.K., Bandyopadhyay, S., and Banerjee, R. (2014).
the best shovel allocation which is the allocation with Benchmarking energy consumption for dump trucks in
the lowest fuel consumption in the haulage operation and mines. Applied Energy, 113, 1382–1396.
therefore with the lower diesel emission. The results of this Suzdaleva, E. and Nagy, I. (2014). Data-based speed-limit-
work serve the foundation for future research that could respecting eco-driving system. Transportation Research
include the fact that, in practice, the waiting time varies Part C: Emerging Technologies, 44, 253–264.
stochastically. Zhang, L. and Xia, X. (2015). An integer programming
approach for truck-shovel dispatching problem in open-
REFERENCES pit mines. Energy Procedia, 75, 1779–1784.
Abdelaziz, E., Saidur, R., and Mekhilef, S. (2011). A
review on energy saving strategies in industrial sector.
Renewable and Sustainable Energy Reviews, 15(1), 150–
168.
222