Disasters often cause significant loss of life and property, which can be both natural (earthquakes, hurricanes) and man-made (terrorist attacks, chemical leakages). Statistics show that more than 500 disasters are estimated to strike our planet every year, around 75,000 people are killed and more than 200 million others are impacted [1] [2] . Especially in the 21st century, disasters occur more frequently, so most countries have established several series of emergency plan to deal with such disasters. However, only relying on emergency plans, the emergency relief supplies (especially grain) cannot be delivered to the stricken area in time, because the disasters generally occurs rapidly and they are difficult to predict, the rapid response times and huge demand is required, and even more emergency facilities may be damaged.

The Chinese Government has established a set of central and local grain reserve system, in response to unexpected public emergency events. However, the current China’s grain reserves are mostly unprocessed, so it needs time and conditions to transform the raw grain into processed grain. The 2008 Wenchuan earthquake provides an excellent example that directly revealed the shortage of processed grain reserves. Therefore, China has been making the planning and construction of processed grain reserve system. Many national policies such as “National Food Security and Long-Term Planning Framework (2008-2020)” and “Food Industry Development Plan during the 12th Five-Year Period”, indicate that the number of provincial designated emergency food processing enterprise increases from 1700 to 2000, the number of food supply industry increases from 4000 to 5000, and the processed grain emergency reserves (EPI) of city should meet the consumption needs of more than 15 days.

At present, although the theorists pay more attention to emergency rescue, concerning to the study about the emergency supply of grain is very insufficient. Through analyzing the operation of emergency supply of grain, it is not difficult to find that the emergency supply system of grain can be viewed as a three-stage supply chain formed by the raw grain suppliers, core enterprise and supply outlets. Due to the short shelf life of processed grain, high storage conditions, requirement of quick rotation, there will lead to a substantial increase in the cost of processed grain reserves, if reserved by local government. Therefore some enterprises which have some production and reserve capacity are assigned and commissioned to stock up on processed grain and also receive subsidies from local government. In the absence of disaster, core enterprise may carry out normal production, marketing and other routine operations on the premise of a certain quality and quantity of processed grain reserves. Once disaster occurs, core enterprise need to distribute the processed grain reserved in advance to the supply outlets of stricken area for the first time, to gives priority to the demands from stricken area.

The current researches on emergency network optimization are mainly concentrated in emergency facility location. The traditional facility location problem can be divided into: p-median problem, p-center problem and covering problem [3] -[7] . And covering problem is the most widely used model in facility location problem, especially in the location of emergency facilities. With the further research, scholars have gradually realized that the basic hypothesis of traditional covering problem is not reasonable in many cases. When the distance between the disaster area and facilities is less than certain distance (time), that was thought to be completely covered; otherwise, the disaster area will not be covered, this assumption is called 0 - 1 covering problem. Therefore scholars put forward a series of improvement ideas, such as the maximal covering location with coverage gradually [8] , the set covering model with multi-coverage [9] [10] , multiple quantity and quality of coverage model [11] , multistage covering location model [12] -[14] , the maximal covering location with backup facilities [15] , and other facility location models with different constraints [16] -[24] .

Through the summary of many facility location models, we found that the current related research only limits the conventional emergency facilities or general emergency materials reserve site selection problem, and ignores the study of the emergency supplies of grain, a special emergency material, which has obviously different characteristics from the general emergency materials in the process of emergency stockpile. Once the disaster occurred, a huge amount of grain supply is needed from disaster area. Because of the short shelf life, grain is very strict with reserve environment and needed to rotate on a regular basis in the process of reserves. So, it is of great importance to research on grain emergency network, especially to determine the scope of reserve enterprise, location, service scope and inventory strategy.

For the number of reserve enterprise is uncertain, the established location-inventory model of grain emergency network is a NP-Hard problem. Classic algorithms for solving location-inventory problems is the Lagrangian relaxation algorithm [25] , but with more complicated constraints and planning scope expands, the algorithm need a longer computation time. So in this paper, we use the genetic algorithm to solve the location-inventory model.

1) Coding work

Choose dimensional 0 - 1 integer row vector as chromosome structure. When the value of the bit is equal to 1, the candidate enterprise is selected, or that is not selected.

2) Initialization

Set evolution algebra counter, set maximum evolution algebra, randomly generated an individual as the initial group.

3) Calculation the fitness of each individual

Step 1: Assign each supply outlet to its nearest reserve enterprise, until all retailers are assigned;

Step 2: Test on each supply outlet in turn. Delete the supply outlet from the enterprise which it was originally assigned to, and then assign this supply outlet to the other reserve enterprise. If the new scheme can reduce the total cost, update this scheme; if not, assign this supply outlet to another reserve enterprise until the total cost cannot be reduced by assigning this supply outlet to any reserve enterprise.

Step 3: If all the supply outlets get through the examination of step 2 and the original scheme has not changed, then in step 4; if not, return to step 2 to examine each supply outlet in turn.

Step 4: According to the above steps to get the optimal scheme and its corresponding location decision variables. Calculate the corresponding chromosome fitness through objective function (4), and terminate the calculation.

4) Choose operation, crossover operation and variation operation

The roulette wheel selection method is adopted for the choice of individuals. First calculate the fitness of each

individual, and then choose probability method to guarantee the excellent genetic structure by

heredity to the next generation of most likely. Using discrete restructuring to make crossover operation can produce new individual more effectively from two similar structures in the parent. Generation replacement method is used for population iteration, to keep only one when appear some same individual and avoid the GA convergence too fast, can also prevent the solution into a local optimal solution at the same time.

5) Conditions for the termination judgment

If, the one with the evolutionary process have been with the maximum fitness individual as optimal solution output, terminate the calculation. When the algorithm terminates, the best chromosome is the optimal solution of the objective function [26] [27] .

In this section, we present computational results from an experiment. We assume that there are 10 candidate reserve enterprises and 50 supply outlets, and their positions are random uniform distribution on the [0, 100] of the plane. In addition randomly generated in uniform distribution [300, 1000], randomly generated in uniform distribution [100, 3000], randomly generated in uniform distribution [1000, 3000], randomly generated in uniform distribution [1/360, 3/360], , , ,.

Table 1 shows the parameters for computational operation that we assumed.

The computational results are shown in Table 2 and Figure 1.

Table 2 shows that there are 10 candidate reserve enterprises and 50 supply outlets in the computational experiment. The optimal objective function value is 475,000 and 5 candidate reserve enterprises have been selected finally. The calculation process has elapsed for 16.5 seconds after 76 evolutions. Shorter elapsed time and less number of iterations show the effectiveness of the algorithm.

In Figure 1, there are three kinds of symbol depict the reserve enterprises and supply outlets. Among them, “o” was on behalf of the supply outlets, “+” was on behalf of the selected reserve enterprises and “△” was on behalf of the unselected reserve enterprises. We can find that 5 candidate reserve enterprises have been selected finally, in order to minimize the total cost.

In this computational examination, we have designed genetic algorithm to solve the model, which does not mean that only genetic algorithm is appropriate. We choose genetic algorithm is just because of its advantage, which has good global search ability and can quickly search out all the solution in the solution space, and won’t get into local optimal solution of the rapid decline in trap.

We have presented a facility location model for grain emergency network that incorporates inventory factors and rotation mechanism, and then designed genetic algorithm based on Matlab to solve the model. Finally we verified the feasibility and effectiveness of the algorithm by computational examples.

We assume that the storage and service capacity of reserve enterprise are unlimited in this paper, but in reality, they are key points which should be considered and cannot be ignored. The storage limit of reserve enterprise will affect the final number and location assigned to the supply outlets, as well as their own inventory control strategy. And also, the construction cost will increase with the ascension of the capacity, which is no longer a fixed value but a function about the capacity. So the grain emergency network optimization with capacity limits will be a further development direction.

This article could not have been completed without the help and support of many people. We express our deepest gratitude to the reviewers and editors, whose able guidance helped in the publication of this article. This work is supported by National Social Science Fund (Grant No. 14BGL062), the National Science & Technology Pillar Program during the 12th Five-Year Plan Period (Grant No. 2014BAL07B04, 2014BAL07B05) and the Beijing Modern Logistics Research Base Open Project (Grant No. JD2013003).