^{1}

^{2}

^{*}

^{1}

Cold-chain demand of fresh agricultural products is increasing in China , while network layout of cold-chain logistics is in disorder and its cost is huge. To address this problem, this paper casts an optimal model of cold-chain logistics network and tackles it with genetic algorithms. This optimal model takes running total cost of logistics network as the objective, and embeds a nonlinear mixed integer programming including two assignment issues. The model determines optimal layout and logistics management for pre-cooling stations and logistics center for fresh agricultural products. Our main contribution is to consider construction cost and operation cost of cold chain logistics simultaneously. Case study illustrates the effectiveness of the model.

Guangdong has various characteristic fruits, such as lychee, banana, strawberry, dragon fruit, longan, etc., and its output of fruit is huge. The output of fruits such as cherry tomato, cucumber, cabbage and mustard, is also at the top in the list. In addition, there are many sea areas adjacent to Guangdong with various aquatic products, which provide a broad development prospects for cold-chain logistics. However, the cold-chain food is still self-producing and self-selling in most areas of Guangdong, and the refrigeration circulation rate is low comparing to some other provinces [

Some references have studied in similar issue. Xie and Zhao [_{2} emissions were considered simultaneously. Zheng, Li and Wu [

Our main contribution in this paper is that, the construction cost and operation cost of cold chain logistics are also considered simultaneously, to establish network optimization model for cold chain logistics. Unlikely with the other research, we study the double layer structure of cold chain logistics network, that is, layer one is from places of origins to pre-cooling station while layer two is from pre-cooling station to distribution centers.

The planning of cold-chain network is to determine the layout and quantity of logistics distribution nodes, including producing area of agricultural products, pre-cooling stations, refrigeration house, distribution center, etc.

Different fresh agricultural products have different refrigerated transportation conditions, such as precooling mode, refrigerated transportation mode, temperature and humidity, etc. Aquatic products and vegetables are generally pre-cooled in the production area, and the pre-cooling method is relatively backward. Road and railway transportation are main methods of trans-regional transportation for vegetables. Generally, meat is pre-cooled in processing plants, and transport under frozen and chilled conditions. Aquatic products have two type of transportation after pre-cooling. The long-distance transportation of quick-frozen aquatic products is generally carried by railway, while the long-distance transportation of fresh aquatic products is mainly carried by air.

In view of different fresh agricultural products are not identical to the requirement of storage and transportation environment, so for a specific category (including aquatic products, vegetables, meat, fruits, etc.), this paper studies and constructs the logistics network of fresh agricultural products of this category, that is, to determine the transportation network between the origin, pre-cooling station and fresh distribution center.

The determination of the logistics network is mainly based on the cost of establishing precooling station, distribution center, storage cost and transportation cost of fresh agricultural products. At the same time, when constructing the logistics network of fresh agricultural products, we should take into account the transportation demand of producing area, the product demand of market, the storage capacity of precooling station and distribution center, etc. The determination of the logistics network takes establishment cost of pre-cooling station and distribution center, storage cost and transportation cost as the main objectives. Meanwhile, in the construction of fresh agricultural products logistics network, we should take into account the transportation needs of the origin, market demand for products, storage capacity of pre-cooling station and distribution center, and other factors.

Before establishing an optimization model, the following assumptions are made for the fresh agricultural product logistics network:

1) Due to the huge investment for establishing a new cold chain logistics distribution center, we assume that the cold chain logistics distribution center has been completed and put into use, and there is no major change in the short term, that is, the location problem of distribution center of cold chain logistics is transformed into the selection problem of existing distribution center.

2) Assume that fresh agricultural produce at the production site will not be directly transported to the cold chain logistics distribution center even if the daily average supply is small, but will be sent to the pre-cooling station for pre-cooling before being transported to the cold chain logistics distribution center. For pre-cooling stations, the number of fresh agricultural products to be pre-cooled per day is called the average daily purchase volume.

3) Because fresh agricultural products have characteristic of diversity, seasonality and the dynamic of logistics operations, in order to accurately reflect the actual cost of unit supply, we takes 2 months as the cost accounting cycle.

4) Since the pre-cooling station has the function of temporarily storing for cold-chain food, and the procurement cycle of pre-cooling station or refrigeration storage is short, the centralized procurement cycle is set within one week in consideration of the capacity limitation.

The optimization model of the cold chain logistics network takes the total cost as objective function. The total cost accounting includes the rental cost, inventory cost, operation cost of logistics activities, and the transportation cost from the place of origin to the pre-cooling station and from pre-cooling station to distribution center. In terms of constraints, each place of origin will only select one pre-cooling station; each pre-cooling station also corresponds to a cold chain logistics distribution center; the distance from place of origin to pre-cooling station and the distance from pre-cooling station to cold chain logistics distribution center should not exceed its maximum limit; the quantity of fresh agricultural products received by pre-cooling station and cold chain logistics distribution center shall not exceed its maximum processing capacity; the quantity of fresh agricultural products processed by all pre-cooling stations is equal to the number of cold chain logistics distribution center. The optimization model of cold chain logistics network is established as follows:

Objective Function:

min [ 60 ∑ t C t L Z T + 60 ∑ j C j H ∑ i q i s X i j + q j t d j t ( C i j + C j t ) 60 P j × ∑ t Z t ∑ j Y j + q j t C t S 60 P j × ∑ t Z t ∑ j Y j + q j t C t J 60 P j × ∑ t Z t ∑ j Y j ] s.t.

∑ j X i j = 1 ,

∑ t X j t = 1 ,

d i j X i j ≤ L ( i , j ) max ,

d j t X j t ≤ D ( i , j ) max ,

∑ j q i j ≤ Q j H ,

∑ j q j t ≤ Z t Q t ,

∑ j q i S X i j P j = ∑ t q j t .

where,

i ―the number of origin of agricultural products, i ∈ [ I ] ;

j ―the number of alternative pre-cooling stations, j ∈ [ J ] ;

t ―the number of alternative distribution centers, t ∈ [ T ] ;

C t L denotes daily average rental cost cold-chain logistics distribution centers, (unit: RMB);

C t S denotes unit inventory cost of the cold-chain logistics distribution centers (unit: RMB);

C t J denotes the operating cost per unit logistics activity of cold chain logistics distribution center t (unit: RMB);

q i S denotes the average daily supply of production site i (unit: ton/day);

C j H denotes the unit processing cost of the pre-cooling station j (unit: RMB);

Q j H denotes the maximum processing capacity of pre-cooling station j (unit: tons);

Q t denotes the maximum inventory capacity of the cold chain logistics distribution center t (unit: tons);

d i j denotes the distance between production site i and pre-cooling station j (unit: km);

d j t denotes the distance between pre-cooling station j and cold chain logistics distribution center t (unit: km);

D ( i , j ) max denotes the maximum distance restriction (unit: km) from the place of production i to the pre-cooling station j;

D ( j , k ) max denotes the maximum distance restriction (unit: km) from the pre-cooling station j to the cold-chain logistics distribution center t;

C i j denotes the unit transportation cost (unit: RMB/ton・km) from the place of production i to the pre-cooling station j;

C j t denotes the unit transportation cost (unit: RMB/ton・km) from the pre-cooling station j to the cold chain logistics distribution center t;

Q min p denotes the minimum number of pre-cooling stations;

Q min D denotes the minimum number of cold-chain logistics distribution centers;

Q max D denotes the maximum number of cold-chain logistics distribution centers;

q i j denotes the supply quantity (unit: ton) from the production site to the pre-cooling station j;

q j t denotes the periodic supply (unit: ton) from the pre-cooling station j to the cold-chain logistics distribution center t;

P j denotes the centralized procurement cycle of the pre-cooling station j, namely, the time span from the beginning of procurement to car-loading (long distance);

X i j , X j t , Y j , Z t are 0 - 1 variables. If they are equal to 1, it means that a function works, while if they are equal to 0, it means that this function is unavailable.

Genetic algorithms (GA) is a mathematical algorithm that simulates the evolutionary process of organisms according to genetic mechanisms. The algorithm has the advantages of self-adaptability and global probabilistic search and its essence is an efficient, parallel and global search method, which can automatically acquire and accumulate information about the search space in the process of searching. The search process is controlled adaptively to obtain the optimal solution. In each generation of the genetic algorithm, individual selection is based on the fitness values of the individuals in the problem domain and the reconstruction method borrowed from natural genetics to generate a new approximate solution, which promotes the evolution of individuals in the population. The new individuals are more capable of adapting to the environment than the original individuals [

Genetic algorithm is a process of global optimization by using the idea of survival of the fittest and natural selection. In the process of genetics, the continuous replication, crossover and mutation of genes promote the continuous reproduction and evolution of the population. Finally, the optimal individual converges, and the process is shown in

As can be seen from

1) The objective function of the problem to be optimized is understood as the adaptability of a biological population to the environment;

2) The variables to be optimized are corresponding to the individuals of the biological population;

3) The algorithm for solving the optimization model is compared with the evolutionary process of the biological population.

The genetic algorithm’s operational flow consists of seven steps [

Step 1: Initialization. Firstly, the gene is encoded to form a chromosome, and then multiple individuals encoded by the chromosome are randomly selected to generate initial population;

Step 2: Evaluate individual’s superiority and inferiority, that is, to determine the fitness function, to calculate the fitness value of the individuals, and to judge whether the optimization criterion is satisfied or not. If satisfied, it means that the final generation is reached, and the optimal solution can be decoded and output; if not, then turn to step three;

Step3: Determine the selection probability according to the fitness value, and execute the selection operator to enter the next generation from the preferred part of the existing population. The higher fitness it is, the more likely the individuals are to be copied;

Step 4: Determine the crossover probability, perform a crossover operator, and preferably select some individuals to enter the next generation;

Step 5: Determine the mutation probability, execute the mutation operator, and select new individuals to enter the new population;

Step 6: If the pre-set evolution times is reached in advance or the optimal individual is obtained in advance, or the optimal individual has not been improved in successive generations, then execute step seven, otherwise, return to step two;

Step 7: Output the individual with the optimal fitness in the population as the optimal solution or satisfactory solution to the model.

The basic operation process of the genetic algorithms is as follows [

1) Initialization: Set the evolution number t = 0, set the maximum evolution number T, and randomly generate M individuals as the initial population P(0).

2) Individual evaluation: Calculate the fitness of each individual in the population P(t).

3) Selection operator: The selection operator is applied to the group. The purpose of the selection is to directly inherit the optimized individuals to the next generation or to generate a new individual through the pairing to regenerate to the next generation. The selection operator is based on the fitness assessment of the individuals in the group.

4) Crossover operator: The crossover operator is applied to the group. The so-called crossover refers to the operation of replacing the partial structure of two parent individuals to generate a new individual. The core function of the genetic algorithms is the crossover operator.

5) Mutation operator: The mutation operator is applied to the population, that is, changes the gene values on certain genomes of individuals in the population. The next generation population P(t + 1) was obtained after selection, crossover and mutation from population P(t).

6) Termination condition judgement: If t = T, the individual with the greatest fitness, which is obtained during the evolution process, is outputed as the optimal solution, and the calculation is terminated.

According to food production in Guangdong Province, the 8 palaces, Shenzhen, zhuhai, Huizhou, Foshan, Zhanjiang, Shantou, Dongguan and Shaoguan, aremain places of agricultural products, which will be coded as 1, 2, 3, 4, 5, 6, 7, and 8, respectively. Suppose that six alternative pre-cooling stations are given, which are coded as I, II, III, IV, V, VI, respectively. Besides the agricultural products logistics distribution center of Huizhou and Zhanjiang, another alternative distribution center of cold chain logistics is given, and they are coded as A, B, and C, respectively. The distances of the eight origins of Shenzhen, Zhuhai, Huizhou, Shunde, Zhanjiang, Shantou, Dongguan, Shaoguan to Guangzhou are shown in

In this example, the distance from 8 origins to 6 alternative pre-cooling stations is shown in

Using the MATLAB R2012 genetic algorithms toolbox to solve the model, through the multiple debugging of the parameters [

Shenzhen | zhuhai | Huizhou | Shunde | Shantou | Zhanjiang | Dongguan | Shaoguan | |
---|---|---|---|---|---|---|---|---|

Guanzhou | 139 | 116 | 140 | 15 | 527 | 466 | 82 | 227 |

i/j | I | II | III | IV | V | VI | Average daily supply |
---|---|---|---|---|---|---|---|

1 | 15 | 22 | 45 | 10 | 39 | 70 | 54 |

2 | 36 | 41 | 50 | 73 | 53 | 92 | 72 |

3 | 20 | 39 | 17 | 33 | 28 | 53 | 87 |

4 | 22 | 34 | 87 | 54 | 30 | 102 | 61 |

5 | 322 | 245 | 507 | 112 | 23 | 210 | 91 |

6 | 32 | 155 | 466 | 76 | 51 | 409 | 60 |

7 | 40 | 72 | 64 | 37 | 145 | 80 | 50 |

8 | 104 | 250 | 77 | 91 | 46 | 126 | 78 |

No. | Place of origins | Pre-cooling station | Cold chain logistics distribution center |
---|---|---|---|

1 | 1, 7 | IV | A |

2 | 2, 3, 4, 5, 6, 8 | V | B |

The problem of cold chain logistics network layout optimization is a complex but significant problem. The fundamental goal is to solve the current network layout chaos and high cost-consuming for domestic cold chain logistics whose demand is growing. Therefore, this paper uses genetic algorithm to optimize the cold chain logistics network, so as to determine the layout of the pre-cooling station and the cold chain logistics distribution center in the cold chain logistics network, and solve transportation management issues from the production area to the pre-cooling station, and from the pre-cooling station to cold chain logistics distribution center in logistics network. The example shows that the established optimization model for cold chain logistics network can effectively solve the network layout problem and transportation management problem, and reduce the operating cost of cold chain logistics. The model also has good adaptability and reference to other logistics network layout problems. The future work is to collect more data in cold chain logistics in Guangdong and other provinces, and to apply this model to tackle these network layout problems.

This research was partly supported by the NSF of Guangdong under grant no. 2018A0303130136, the Soft Science Project of Guangdong under grant no. 2015A070704059 and 2015A030402008, Special Fund Project of Guangdong Provincial Rural Revitalization Strategy under grant no. Yuecainong [

The authors declare no conflicts of interest regarding the publication of this paper.

Liang, K.H., Zhang, W.F. and Zhang, M. (2019) Optimization Model of Cold-Chain Logistics Network for Fresh Agricultural Products. Journal of Applied Mathematics and Physics, 7, 476-485. https://doi.org/10.4236/jamp.2019.73034