Fuel Consumption Estimation System and Method with Lower Cost

This study proposes a fuel consumption estimation system and method with lower cost. On-board units can report vehicle speed, and user devices can send fuel information to a data analysis server. Then the data analysis server can use the proposed fuel consumption estimation method to estimate the fuel consumption based on driver behaviours without fuel sensors for cost savings. The proposed fuel consumption estimation method is designed based on a genetic algorithm which can generate gene sequences and use crossover and mutation for retrieving an adaptable gene sequence. The adaptable gene sequence can be applied as the set of fuel consumption in accordance with the pattern of driver behaviour. The practical experimental results indicated that the accuracy of the proposed fuel consumption estimation method was about 95.87%.


Introduction
With the development of the economic environment and the evolution of mobile communications, the intelligent transportation system has been more and more popular for obtaining fleet management services for the logistics industries and bus carriers [1].A logistics company may have hundreds or thousands of trucks to provide freight services.However, the fuel cost of these industries is the most important challenge of fleet management services.For instance, Taiwan Institute of Economic Research (TIER) reported that the fuel costs of logistics industries and bus carriers were about 35.8 billion dollars [2] and 3.4 billion dollars [3], respectively, in Taiwan in 2015.Therefore, monitoring and saving fuel consumption efficiently can improve the profits for the logistics industries and bus carriers and reduce the air pollution from carbon dioxide (i.e., CO 2 ) for city governance [4,5].
For the measurement of fuel consumption, some studies used fuel sensors to detect the remaining quantity of fuel and calculated the differences among the remaining quantities.Furthermore, the data from the on-board diagnostics (OBD) could also be retrieved to obtain the fuel system status, vehicle speed, and engine revolutions per minute (RPM) [6][7][8][9].Although these methods can measure the fuel consumption, fuel sensors and OBD devices should be equipped to report the fuel quantity data periodically, albeit with higher cost.For instance, the cost of a fuel level sensor is about 150,000 dollars, and the cost of data communications is about 13,000 dollars per month in Taiwan for a fleet of 1000 vehicles.Moreover, these methods cannot support estimating fuel consumption in accordance with driver behaviours.For the estimation of fuel consumption based on driver behaviours, some studies have proposed using gravity sensors, accelerometers, and OBD devices to collect and analyse the data of azimuth, acceleration, movement records, and fuel quantities [10][11][12][13].Some studies used a genetic algorithm (GA) and a neural network to analyse fuel consumption and classify driver behaviours [14][15][16].Although the relation between fuel consumption and driver behaviour can be estimated, sensors and OBD devices are required in these methods.Furthermore, the measurement errors of sensors and OBD devices have not been discussed, and signal interference may lead to large estimation errors.
Therefore, a lower cost solution for a system is proposed and implemented to estimate the fuel consumption for the logistics industries.In this system, On-board units (OBUs) can send the information of movement to a data analysis server, and users can input and send the information of fuel quantity through user devices.Then the data analysis server can analyse the movement information and the fuel quantity information to estimate the fuel consumption based on driver behaviours without fuel sensors for saving costs.The proposed fuel consumption estimation method is designed based on a GA [17][18][19] which can generate gene sequences and use crossover and mutation to retrieve an adaptable gene sequence.The adaptable gene sequence can be applied as the set of fuel consumption in accordance with the pattern of driver behaviour.
In the next section, the architecture of the proposed consumption estimation system is described in detail.Section 3 proposes a consumption estimation method based on GA.In Section 4, practical experiments are designed to evaluate the proposed methods, and the results of these experiments are also analyzed and discussed in this section.The conclusions and future work of this paper are presented in Section 5.

Fuel Consumption Estimation System
The proposed fuel consumption estimation system includes OBUs, user devices, a data analysis server, and a database server (shown in Figure 1).The OBU can send the movement information which includes the timestamp, location (i.e., longitude and latitude), and vehicle speed to the data analysis server.Furthermore, a user can use his user device to input the fuel quantity after refuelling.The data analysis server can store these data in the database server and perform a fuel consumption estimation method to estimate fuel consumption in accordance with the pattern of driver behaviour.
month in Taiwan for a fleet of 1000 vehicles.Moreover, these methods cannot support estimating fuel consumption in accordance with driver behaviours.
For the estimation of fuel consumption based on driver behaviours, some studies have proposed using gravity sensors, accelerometers, and OBD devices to collect and analyse the data of azimuth, acceleration, movement records, and fuel quantities [10][11][12][13].Some studies used a genetic algorithm (GA) and a neural network to analyse fuel consumption and classify driver behaviours [14][15][16].Although the relation between fuel consumption and driver behaviour can be estimated, sensors and OBD devices are required in these methods.Furthermore, the measurement errors of sensors and OBD devices have not been discussed, and signal interference may lead to large estimation errors.
Therefore, a lower cost solution for a system is proposed and implemented to estimate the fuel consumption for the logistics industries.In this system, On-board units (OBUs) can send the information of movement to a data analysis server, and users can input and send the information of fuel quantity through user devices.Then the data analysis server can analyse the movement information and the fuel quantity information to estimate the fuel consumption based on driver behaviours without fuel sensors for saving costs.The proposed fuel consumption estimation method is designed based on a GA [17][18][19] which can generate gene sequences and use crossover and mutation to retrieve an adaptable gene sequence.The adaptable gene sequence can be applied as the set of fuel consumption in accordance with the pattern of driver behaviour.
In the next section, the architecture of the proposed consumption estimation system is described in detail.Section 3 proposes a consumption estimation method based on GA.In Section 4, practical experiments are designed to evaluate the proposed methods, and the results of these experiments are also analyzed and discussed in this section.The conclusions and future work of this paper are presented in Section 5.

Fuel Consumption Estimation System
The proposed fuel consumption estimation system includes OBUs, user devices, a data analysis server, and a database server (shown in Figure 1).The OBU can send the movement information which includes the timestamp, location (i.e., longitude and latitude), and vehicle speed to the data analysis server.Furthermore, a user can use his user device to input the fuel quantity after refuelling.The data analysis server can store these data in the database server and perform a fuel consumption estimation method to estimate fuel consumption in accordance with the pattern of driver behaviour.

On-Board Unit
The OBU includes a position module, a middleware module, and a communication module as follows.

1.
The position module can support a global positioning system (GPS) to receive and analyze satellite signals for estimating the location (i.e., longitude and latitude) and speed of the vehicle.

2.
The communication module can support the techniques of long term evolution (LTE), and the OBU can connect with the data analysis server through the communication module and cellular networks.

3.
The middleware module can support hypertext transfer protocol (HTTP) and representational state transfer (REST), and the OBU can periodically call application program interfaces (APIs) and send the movement information (e.g., OBU ID, car type, driver ID, timestamp, longitude, latitude, and vehicle speed) to the data analysis server.
In this study, an OBU stores an OBU ID, a car type, and a driver ID.There are C N OBUs, T N car types, and D N drivers.The OBU can send the movement information to the data analysis server every 30 s.For instance, Driver 1 drove the car which was equipped with OBU 1 on 1 January 2015, and the type of this car was Type 1.The position module of OBU 1 was used to estimate the location (i.e., 120.5423383 • E and 24.09490167 • N) and speed (i.e., 44 km/h) at 06:00:00, and the middleware module was used to call REST APIs of the data analysis server for the transmission of the movement information (shown in Table 1).

User Device
The user device includes a user interface, a middleware module, and a communication module as follows.

1.
The user interface can be used to input the OBU ID, timestamp, and fuel quantity after refuelling.

2.
The communication module can support the techniques of LTE, and the connection between a user device and the data analysis server can be built through the communication module.

3.
The middleware module can support the techniques of HTTP and REST, and the user device can send the information such as OBD ID, timestamp, and fuel quantity to the data analysis server through the middleware module.
In this study, a user can input the information of fuel quantity (e.g., OBU ID, timestamp, and fuel quantity) through the user interface of the user device after refuelling, and the fuel quantity information can be sent to the data analysis server through the middleware module.For instance, the car which was equipped with OBU 1 was refueled with 43.04 L of gas at 18:51:00 on 5 January 2015.Then a user inputted the OBU ID (i.e., OBU 1), timestamp (i.e., 5 January 2015 18:51:00), and fuel quantity (i.e., 43.03 L) through the user interface of the user device, and the middleware module was Symmetry 2017, 9, 105 4 of 15 used to send the inputted data (i.e., the fuel quantity information) to the data analysis server (shown in Table 2).

Data Analysis Server
The data analysis server includes a middleware module, communication module, and an optimized composition module as follows.

1.
The middleware module can obtain several REST APIs to receive the movement information (e.g., OBU ID, car type, driver ID, timestamp, longitude, latitude, and vehicle speed) and the fuel quantity information (e.g., OBU ID, timestamp, and fuel quantity) from the OBUs and user devices through HTTP.These data can be stored in a database server.

2.
The communication module can support Ethernet and build the connections among the data analysis server and other devices (e.g., OBUs, user devices, and a database server).

3.
The optimized composition module can use the proposed fuel consumption estimation method to collect and analyse the movement information and fuel quantity information for generating the estimated results of fuel consumption.
In this study, the data analysis server can request Google Maps [20] or Chunghwa Telecom GeoWeb [21] to find the corresponding road type of the location through the middleware module.For instance, the road type of the location which is positioned at 120.5423383 • E and 24.09490167 • N is an urban road.The data analysis server can add the column of road type into the movement information (shown in Table 3), and the modified movement information can be stored in a database server.

Database Server
The database server includes a storage module, a computation module, and a communication module as follows.The communication module can support Ethernet, and the connection between the database server and the data analysis server can be built by this module.

2.
The computation module can receive the requests from the data analysis server through the communication module and access the storage module in accordance with the requests.

3.
The storage module can perform the operations of creation, update, deletion, and query.

Fuel Consumption Estimation Method
This study proposes a fuel consumption estimation method which includes a movement information collection method, a fuel information collection method, and an optimized composition method.The details of each method are illustrated in the following subsections.

Movement Information Collection Method
The process of the movement information collection method includes: (1) receiving the movement information from the OBUs; (2) analyzing and storing the movement information; and (3) calculating the amount of each vehicle speed interval (i.e., driver behaviour in this study) for each OBU and each driver during a time interval.
In this study, the movement information collection method can be performed by the optimized composition module of the data analysis server to retrieve the modified movement information (shown in Table 3).The vehicle speed (v) of each movement record can be converted into a vehicle speed interval, and the amount of each vehicle speed interval for each OBU and each driver during a time interval can be calculated.This study chooses 10 km/h as a vehicle speed interval and one year as a time interval for the estimation of driver behaviour.For instance, Driver 1 drove a car which was equipped with OBU 1 during 2015; c C,D 1,1,1 records which include idle speed (i.e., the value of v is zero) are reported by OBU 1 during 2015; c C,D 1,1,2 records which include the speed between 0 km/h and 10 km/h are reported by OBU 1 during 2015; consequently, c C,D 1,1,14 records which include the speed higher than 120 km/h are reported by OBU 1 during 2015.Furthermore, Driver D N drove a car which was equipped OBU C N during 2015; c C,D C N ,D N ,1 records which include idle speed are reported by OBU C N during 2015; c C,D C N ,D N ,2 records which include the speed between 0 km/h and 10 km/h are reported by OBU C N during 2015; consequently, c C,D C N,D N ,14 records which include the speed higher than 120 km/h are reported by OBU C N during 2015 (shown in Table 4).
Table 4.The movement information of each OBU and each driver during 2015 (The unit of v is km/h).

Fuel Information Collection Method
The process of the fuel information collection method includes: (1) receiving the fuel quantity information from the user devices; (2) analyzing and storing the fuel quantity information; and (3) calculating the amount of fuel quantities for each OBU and each driver during a time interval.
In this study, the fuel information collection method can retrieve the fuel quantity information (shown in Table 2) and analyse the amount of fuel quantities for each OBU and each driver during a time interval.This study chooses one month or one year as a time interval.For instance, Driver 1 drove a car which was equipped with OBU 1 during 2015; the summary of fuel quantities of OBU 1 during January 2015 is ).Furthermore, Driver D N drove a car which was equipped with OBU C N during 2015; the summary of fuel quantities of OBU C N during January 2015 is

Optimized Composition Method
In this subsection, the design of the optimized composition method is presented in Section 3.3.1,and a case study of this method is given in Section 3.3.2.

The Process of the Method
The process of the optimized composition method includes: (1) receiving the patterns of driver behaviours from the movement information collection method; (2) receiving the patterns of fuel consumption from the fuel information collection method; and (3) performing a GA (shown in Figure 2) to analyse the set of fuel consumption in accordance with the pattern of driver behaviour.The steps of the optimized composition method are illustrated as follows.

1.
The values of the parameters including the amount of initial maternal DNA (deoxyribonucleic acid) sequences (count g ), the number of evolution times (count c ), the maximum number of iterations (count i ), crossover rate (α), and mutation rate (β) are initially given in this step.

2.
The model of fitness function is designed for finding the cost of each DNA sequence which includes several chromosomes.

3.
The process of the initial population can generate count g maternal DNA sequences.

4.
Each DNA sequence can be adopted into the model of fitness function for the cost calculation.

5.
The process of the convergence check can be performed to check the values of the number of evolution times (count c ) and the maximum number of iterations (count i ).If the number of evolution times (count c ) is equal to the maximum number of iterations (count i ), the adaptable DNA sequence with the lowest cost is outputted as the estimated results of the fuel consumption based on the driver behaviour; otherwise the number of evolution times (count c ) is increased by one.6.
The process of gene selection can select two of the maternal DNA sequences for crossover and mutation.7.
The process of gene crossover can generate a child's DNA sequences in the first generation in accordance with the crossover rate (α) and the maternal DNA sequences in the first generation.8.
The process of gene mutation can generate a child's DNA sequences in the second generation in accordance with the mutation rate (β) and the maternal DNA sequences in the second generation.9.
The process of gene reproduction can support two new generated child's DNA sequences being substituted for two original maternal DNA sequences which have the highest cost.1.The values of the parameters including the amount of initial maternal DNA (deoxyribonucleic acid) sequences ( g count ), the number of evolution times ( c count ), the maximum number of iterations ( i count ), crossover rate (α), and mutation rate (β) are initially given in this step.

A Case Study
The parameters of GA in this study are adopted as follows: the amount of initial maternal DNA sequences (count g ) is 14; the initial number of evolution times (count c ) is 0; the maximum number of iterations (count i ) is 1000; the crossover rate (α) is 100%; the mutation rate (β) is 7%.A DNA sequence q C,D j,i = q C,D j,i,1 , q C,D j,i,2 , ..., q C,D j,i,14 includes 14 chromosomes (i.e., q C,D j,i = 14) for the j-th OBU driven by the i-th driver.Each chromosome is encoded as a float.For instance, the parameter q C,D j,i,1 can be used to estimate the quantity of fuel consumption at idle speed during each 30 s period for the j-th OBU driven by the i-th driver.Furthermore, the model of fitness function s = 14 (shown in Figure 3) is adopted to estimate the cost of each DNA sequence.The unit of cost is a liter in this model.The best DNA sequence has the lowest cost (i.e., s = 14 For the initial population, 14 maternal DNA sequences are randomly generated, and each DNA sequence includes 14 chromosomes.For instance, the first maternal DNA sequence is q 1 14 for the j-th OBU driven by the i-th driver; the second maternal DNA sequence is q 2 C,D j,i = q 2 C,D j,i,1 , q 2 C,D j,i,2 , ..., q 2 C,D j,i,14 for the j-th OBU driven by the i-th driver; consequently, the fourteenth maternal DNA sequence is q 14 C,D j,i = q 14 C,D j,i,1 , q 14 C,D j,i,2 , ..., q 14 C,D j,i,14 for the j-th OBU driven by the i-th driver (shown in Table 5).
,..., , = for the j-th OBU driven by the i-th driver; the second maternal DNA sequence is ,..., , = for the j-th OBU driven by the i-th driver; consequently, the fourteenth maternal DNA sequence is ,..., , = for the j-th OBU driven by the i-th driver (shown in Table 5). ) This study gives a case study of OBU 1 driven by Driver 1 to explain the process of the optimized composition method.The set of movement information of OBU 1 driven by Driver 1

Chromosome
Chromosome 1 Chromosome 2 . . .Chromosome 14 . . .DNA Sequence 14 (q 14 This study gives a case study of OBU 1 driven by Driver 1 to explain the process of the optimized composition method.The set of movement information of OBU 1 driven by Driver 1 during 2015 1,1,2 , ..., c C,D 1,1,14 = {103100, 66752, ..., 4}, and the fuel quantity of OBU 1 driven by Driver 1 during 2015 is recorded as 10921.364L (i.e., Q C,D 1,1 = 10921.364).Furthermore, the 14 maternal DNA sequences are randomly generated as shown in Table 6.For example, the first maternal DNA sequence is q C,D 11,1 = {0.013249146, 0.018487159, ..., 0.551971137}; the second maternal DNA sequence is q C,D 21,1 = {0.016574516, 0.02331678, ..., 0.553625064}; consequently, the fourteenth maternal DNA sequence is q C,D 141,1 = {0.01539256, 0.021892833, ..., 0.555117159}.For finding the cost of each DNA sequence, the set of chromosomes in a maternal DNA sequence is adopted into the fitness function.In this study, the cost of the h-th maternal DNA sequence is defined for the j-th OBU driven by the i-th driver.For instance, the cost of the first maternal DNA sequence is calculated by Equation ( 1) for OBU 1 driven by Driver 1; the cost of the second maternal DNA sequence is calculated by Equation ( 2) for OBU 1 driven by Driver 1; consequently, the cost of the fourteenth maternal DNA sequence is calculated by Equation ( 3) for OBU 1 driven by Driver 1.For the convergence check, the adaptable DNA sequence with the lowest cost is output as the estimated results of the fuel consumption based on driver behaviour if the number of evolution times (count c ) is equal to the maximum number of iterations (count i ); otherwise the number of evolution times (count c ) is increased by one.
For gene selection, this study uses the roulette wheel selection method to select two of the maternal DNA sequences.For instance, DNA Sequence 1 (q C,D 11,1 = {0.013249146, 0.018487159, ..., 0.551971137}) and DNA Sequence 2 (q C,D 21,1 = {0.016574516, 0.02331678, ..., 0.553625064}) (shown in Table 7) are selected as the maternal DNA sequences in the first generation for OBU 1 driven by Driver 1.For gene crossover, the 1-point crossover method is performed in accordance with the crossover rate (α).For instance, the value of the crossover point is randomly determined as 2. Then the two child's DNA sequences (shown in Table 8) are generated as Equations ( 4) and ( 5) according to the two maternal DNA sequences (shown in Table 7).The line in Equations ( 4) and ( 5) is the crossover point.
= {0.013249146,|0.02331678 , ..., 0.553625064} DNA Sequence 2 q C,D = {0.016574516,|0.018487159 , ..., 0.551971137} (5) For gene mutation, the set of binary vectors (η = {η 1 , η 2 , ..., η 14 }) is randomly generated in accordance with the mutation rate (β).If the value of η n is equal to one, the value of the n-th chromosome will be changed after the process of gene mutation.For instance, the child's DNA sequences in the first generation (shown in Table 8) are selected as the maternal DNA sequences in the second generation (shown in Table 9).Furthermore, the set of binary vectors is randomly generated as η = {1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}.Then the two child's DNA sequences (shown in Table 10) are determined according to the two maternal DNA sequences (shown in Table 9).For gene reproduction, the two new generated child's DNA sequences are substituted for two of the original maternal DNA sequences which have the highest cost.For instance, the costs of DNA Sequence 2 (q 2 C,D j,i ) and DNA Sequence 14 (q 14 C,D j,i ) are higher than the costs of the other maternal DNA sequences in Table 6 for OBU 1 driven by Driver 1.Therefore, the two new generated child's DNA sequences in Table 10 are substituted for DNA Sequence 2 (q 2 C,D j,i ) and DNA Sequence 14 (q 14 C,D j,i ) (shown in Table 11).Furthermore, the costs of these two reproduced DNA sequences are determined by Equations ( 6) and (7)  After gene reproduction, the convergence check is performed repeatedly.The adaptable DNA sequence with the lowest cost is output as the estimated results of the fuel consumption based on the driver behaviour when the number of evolution times (count c ) is equal to the maximum number of iterations (count i ).In the case of OBU 1 driven by Driver 1 in this study, the adaptable DNA sequence with the lowest cost, DNA Sequence 14 (q 14 C,D j,i ), is output as the estimated results of the fuel consumption based on the driver behaviour (shown in Equation ( 8)).Therefore, the parameter q C,D 1,1,1 (i.e., 0.0125003 L) is used to estimate the quantity of fuel consumption at idle speed during each 30 s Symmetry 2017, 9, 105 11 of 15 period; consequently, the quantity of fuel consumption at higher speed (>120 km/h) during each 30 s period is estimated as q C,D 1,1,14 (i.e., 0.551971137 L).

Practical Experimental Results and Discussions
This section describes the practical experimental environments and data, and some experimental designs are given to evaluate the proposed method in the following subsections.

Experimental Environments
The practical experimental data including the movement information and fuel quantity of fifteen trucks during November and December 2016 were collected from an eFMS (e-Fleet Management Service) system which was built by Chunghwa Telecom [22] for the evaluation of the proposed method.This study used the Package 'GA' [23] to implement the GA algorithm for fuel consumption estimation.The movement information and fuel quantity information during November 2016 were used as training data, and the movement information and fuel quantity information during December 2016 were used as testing data to evaluate the performance of the proposed fuel consumption estimation system and method.In the training data, 420,673 movement records were retrieved, and 105,569 km were driven by the trucks; in the testing data, 414,798 movement records were retrieved, and 107,651 km were driven by the trucks (shown in Table 12).Table 13 shows that the fuel quantity was about 31,687.786L from 280 refuelling times in the training data, and the fuel quantity was about 33,164.136L from 286 refuelling times in the testing data.

Experimental Designs
Three cases were designed to evaluate the performances of the fuel consumption estimation with or without road types.Furthermore, this study compared the performances of the proposed method and the fuel economy guide [24].Each case is illustrated as follows.

•
Case 1: The fuel consumption estimation method considered the different road types (e.g., urban road and highway).In this case, each DNA sequence included 28 chromosomes.

•
Case 2: The fuel consumption estimation method did not consider the different road types.In this case, each DNA sequence included 14 chromosomes.

•
Case 3: The significant differences between the results of Case 1 and Case 2 were evaluated.
In this study, the mean absolute percentage error (MAPE) was used to measure the accuracy of the fuel consumption estimation method by Equation (9).Furthermore, the t-test and F-test were adopted to evaluate the significant differences between the results of Case 1 and Case 2.

Experimental Results and Discussions
This subsection used the records during November 2016 (i.e., training data) with the proposed method to retrieve the adaptable DNA sequence as the estimated fuel consumption based on the driver behaviour.Then the records during December 2016 (i.e., testing data) were used to test the performance of the proposed method.In Case 1 (i.e., 28 chromosomes), the testing results indicated that the accuracies of the proposed method and fuel economy guide were 95.54% and 24.06% (shown in Table 14), respectively.Furthermore, Table 15 shows that the accuracy of the proposed method was 95.87% in Case 2 (i.e., 14 chromosomes).The Bureau of Energy did not consider driver behaviour, traffic condition, and weather when it performed its experiments, so larger errors of fuel consumption estimation were generated in accordance with the fuel economy guide [24].Therefore, the proposed method can provide precise fuel consumption estimation in practical environments for the logistics industries.The practical data were collected from 15 trucks to evaluate the accuracy of the proposed method in Cases 1 and 2. In Case 3, a two-tailed t-test was performed to determine the differences between the mean accuracy of the proposed method in Case 1 and Case 2. The p-value of this two-tailed t-test was measured as 0.6396 which is higher than the alpha level of 0.05, so the H0 (i.e., null hypothesis) was accepted under 14 degrees of freedom (i.e., 14 = 15 − 1).Furthermore, this study performed a two-tailed F-test to determine the differences between the variance of the accuracy of the proposed method in Case 1 and Case 2. The p-value of this two-tailed F-test was measured as 0.1160 which is higher than the alpha level of 0.05, so the H0 (i.e., null hypothesis) was accepted under 14 degrees of freedom (i.e., 14 = 15 − 1).There were no significant differences between the results of Case 1 and Case 2. Therefore, the fuel consumption estimation without road types can be adopted for reducing the computation costs in the future.

Conclusions and Future Work
As monitoring and reducing fuel costs is the greatest challenge for the logistics industries and bus carriers, a fuel consumption estimation system and method based on a GA is proposed to collect the movement information and fuel quantity information for analyzing the relation of fuel consumption and driver behaviour.This study models the pattern of driver behaviour as a gene sequence, and the GA including crossover and mutation processes is used to retrieve an adaptable gene sequence.The adaptable gene sequence can be applied as the set of fuel consumption in accordance with the pattern of driver behaviour.In practical experimental environments, 835,471 movement records from 15 trucks were collected for the evaluation of the proposed method.The practical results indicated that the accuracy of the proposed fuel consumption estimation method was 95.87%, which was higher than the fuel economy guide from the Bureau of Energy in Taiwan.
In the future, the evaluation model of the driver performance can be designed based on the proposed fuel consumption estimation method.Furthermore, some suggestions and guidelines can be made based on the driver behaviours for the driver's reference.

Figure 1 .
Figure 1.The architecture of the fuel consumption estimation system.Figure 1.The architecture of the fuel consumption estimation system.

Figure 1 .
Figure 1.The architecture of the fuel consumption estimation system.Figure 1.The architecture of the fuel consumption estimation system.

10 .Figure 2 .
Figure 2. The process of the fuel consumption estimation method.

2 . 6 . 9 .
The model of fitness function is designed for finding the cost of each DNA sequence which includes several chromosomes.3. The process of the initial population can generate g count maternal DNA sequences.4. Each DNA sequence can be adopted into the model of fitness function for the cost calculation.5.The process of the convergence check can be performed to check the values of the number of evolution times ( c count ) and the maximum number of iterations ( i count ).If the number of evolution times ( c count ) is equal to the maximum number of iterations ( i count ), the adaptable DNA sequence with the lowest cost is outputted as the estimated results of the fuel consumption based on the driver behaviour; otherwise the number of evolution times ( The process of gene selection can select two of the maternal DNA sequences for crossover and mutation.7. The process of gene crossover can generate a child's DNA sequences in the first generation in accordance with the crossover rate (α) and the maternal DNA sequences in the first generation.8.The process of gene mutation can generate a child's DNA sequences in the second generation in accordance with the mutation rate (β) and the maternal DNA sequences in the second generation.The process of gene reproduction can support two new generated child's DNA sequences being substituted for two original maternal DNA sequences which have the highest cost.10.The costs of the two reproduced DNA sequences can be measured by using the model of fitness function, and the GA is performed repeatedly.

Figure 2 .
Figure 2. The process of the fuel consumption estimation method.

Figure 3 .
Figure 3.The model of the fitness function.

Figure 3 .
Figure 3.The model of the fitness function.

Table 2 .
Fuel quantity information.

Table 3 .
The modified movement information.

Table 6 .
The maternal DNA sequence for OBU 1 driven by Driver 1.

Table 7 .
The maternal DNA sequences in the first generation.

Table 8 .
The child's DNA sequences in the first generation.

Table 9 .
The maternal DNA sequences in the second generation.

Table 10 .
The child's DNA sequences in the second generation.

Table 11 .
The maternal DNA sequences after the first iteration.

Table 12 .
The movement information in practical experimental environments.

Table 13 .
The fuel quantity information in practical experimental environments.

Table 14 .
The accuracy of fuel consumption estimation with road types (i.e., Case 1).

Table 15 .
The accuracy of fuel consumption estimation without road types (i.e., Case 2).