Using Improved Entropy-Cloud Model to Select High-Speed Railway Express Freight Train Service Sites

1School of Transportation and Logistics, Southwest Jiaotong University, Chengdu, Sichuan 610031, China 2National United Engineering Laboratory of Integrated and Intelligent Transportation, Southwest Jiaotong University, Chengdu, Sichuan 610031, China 3Traffic and Transportation Engineering Experiment and Comprehensive Innovation Center, Southwest Jiaotong University, Chengdu, Sichuan 610031, China


Introduction
With the rapid development of China's economy, the demand for Third Party Logistics (TPL) among regional centers keeps increasing, especially for high value-added goods, for example, foods and drinks, electronics, medical supplies, and textile.The goods owners would rather afford to pay higher transport costs for transport service that is convenient, time-saving, and reliable than vice versa.The high-speed railway (HSR) has a lot of advantages such as large transport capacity, low energy consumption, punctual operation, safety, and reliability, which satisfies the requirements of the TPL transportation.There is a huge potential to develop the HSR express serving for the TPL.Also, it is one of the most important strategic planning for China railway to obtain more shares in high value-added goods transportation market.
The high-speed trains provide TPL service based on the HSR express network; the network is formed with sites and edges.In this paper we focus on selecting and classifying the possible HSR stations existing in China.For each HSR station, there are three choices considered: selecting and developing as a high-speed railway express logistics center (LC), selecting and developing as a high-speed railway express logistics base (LB), or no TPL service.The cargos are gathered to the LC and LB by roadway or railway.Each selected site will serve one area so the sites must be located in central cities of each logistics area.
In this paper we formulate a new approach based on improved entropy-cloud model (IECM) to solve the HSR express sites selection problem (SSP).After analyzing the factors related to the high-speed railway express SSP, for example, the existing railway network conditions, traffic environment, express freight market demand, and policy, we try to explain the internal development law of the sites, plan the layout and location of the sites, and divide the level of the sites according to the actual TPL demand.The remainder of this paper is organized as follows.In Section 2, we give a detailed literature review on railway express SSP.Section 3 is devoted to the description of the railway express freight train SSP evaluation indicators; the data of the indicators will be applied as the inputs of the proposed approach.Section 4 is devoted to introducing the translation-corrected improved entropy-cloud model, including formulation, proof, and the SSP outputs of the approach.In Section 5 we give the comparison analysis and results discussions of the case study.Section 6 presents the major conclusions and gives an outline of future research tasks.

Literature Review
This section presents a review of the railway sites selection problem literature that provides the background and establishes a framework for this research.This section also identifies the contributions of the research presented in this paper.First, we will introduce the related literatures about the railway and intermodal hub-and-spoke network design problem.Then, we will introduce other approaches which are used to solve the railway sites selection problem.After that, we will give a conclusion about the existing literatures, as well as the contributions and objectives about the approach used in this paper.
Most of the researches formulate the railway express SSP as a hub-and-spoke network [1] or intermodal hub-andspoke network design problem [2][3][4][5][6].The hub-and-spoke network problem was originally formulated to be a quadratic integer programming model by O'Kelly [7] and then later is referred to as a -hub median problem by Campbell [8].The -hub median problem aims to select  hub facilities from a set of available candidates, to allocate spoke nodes to hubs and to simultaneously route freight flows.There are two kinds of mathematical models to formulate the hub-and-spoke network problem, including path-based mixed integer linear programming models with four-dimensional variables (e.g., [9][10][11][12]) and origin-based mixed integer linear programming models with three-dimensional variables (e.g., [13][14][15]).Meng and Wang [3], both, apply a mathematical program with equilibrium constraints model for the intermodal hub-and-spoke network design problem with multiple stakeholders and multitype containers.There are two main contributions about the model: formulating the user equilibrium (UE) behavior of intermodal operators in route choice for any given network design decision of the network planner and using a cost function that is capable of reflecting the transition from scale economies to scale diseconomies in distinct flow regimes for carriers or hub operators.A hybrid genetic algorithm (HGA) was proposed to solve the model.From a network design perspective, Alumur et al. [4] proposed a linear mixed integer programming model which considered transportation costs and travel time to solve the multimodal hub location and hub network design problem.Considering the same mixed uncertainties in both transportation cost and travel time, Yang et al. [5] developed a novel modeling framework for the intermodal hub-and-spoke network design problem to jointly minimize the expected value of total transportation costs and the maximum travel time requirement in term of critical value.A new hybrid methodology by combining fuzzy random simulation (FRS) technique and multistart simulated annealing (MSA) algorithm was designed to solve the proposed model.Serper and Alumur [6] applied a mixed integer programming model and a variable neighborhood search algorithm to determine the locations and capacities of hubs, which are transportation modes to serve at hubs, allocation of nonhub nodes to hubs, and the number of vehicles of each type to operate on the hub network to route the demand between origin-destination pairs with minimum total cost.There are four kinds of costs in the objective function including fixed costs of establishing hubs with different capacities, purchasing and operational costs of vehicles, transportation costs, and material handling costs.Jeong et al. [1] applied the hub-and-spoke network in railways for the first time, where a central planner is to find transport routes, frequency of service, length of trains to be used, and transportation volume.A linear integer programming model whose objective function includes not only the typical operational cost, but also cost due to the transit time spent by freight in the network is established to formulate the problem, and heuristic algorithms were used to solve large-scale instances occurring in European rail freight systems.In order to solve the hub-and-spoke service network design problem for express freight train running between railway logistic terminals, Hai et al. [16] applied the Principal Component Analysis (PCA) to determine the node's level of the regional logistics hub-and-spoke network.They considered influence factors including GDP, population, number of industrial enterprises, industrial added value, total retail sales of social consumer goods, highway mileage, inland navigation mileage, railway operation mileage, and so on.Sun and Li [17] used Fuzzy Cluster Analysis (FCA) to classify the railway express stations' level.FCA involves four steps: establishing evaluation indicators system, dimensionless processing, and establishing fuzzy relation and fuzzy cluster.Huang et al. [18] established a mixed integer linear programming model to minimize the total operation costs, including fixed costs, variable costs, and penalty cost, and designed a genetic algorithm to solve the model and the scale of the problem.
Moreover, some other approaches have also been applied to solve the railway SSP.Yin et al. [19][20][21] established three different mixed integer programming models to solve railway baggage and parcel logistics base and distribution sites problem.Three different algorithms were designed to solve the problems, including a heuristic algorithm combined scan and local search [19], a two-stage method combined improved scan method and iterative method [20], and a two-stage method combined scan search and tabu search algorithm [21].AHP transforms a flexible multicriteria decision-making problem into a hierarchy with respect to one or more criteria.It has been widely used for a variety of sites selection problem; for example, Mohajeri and Amin [22] applied the methods of AHP and data envelopment analysis (DEA) to find the optimum site for a railway station.There are four main criteria in that paper, including rail-related, passenger services, architecture and urbanism, and economics.There are also 26 subcriteria as well as five alternatives.The comparison matrices of AHP were used to obtain the weights and priorities of the railway station candidates and a DEA model was proposed to determine the optimum site for a railway station.Bersani et al. [23] proposed an AHP based methodological approach to evaluate the optimal locations to install a Train Conformity Check System (TCCS) on a railway section.Ballis [24] used AHP to select the airport-site location on the Island of Samothraki, Greece.The AHP method has also been used for transshipment site selection [25] and transit site selection [26].Since a few years ago, artificial intelligence and expert systems have been used in the SSP for searching the appropriate solutions, such as Geographic Information System (GIS); see Padmaja et al. [27].There are three steps for a SSP.(i) Create vulnerability map with various factors.(ii) Create P map and I map to apply the equation of Pts map.(iii) GIS software locates the best sites.Yu et al. [28] proposed an approach to deal with the hub location problem in urban public transit network design.Firstly, they introduced an indicator defined as passenger attraction to describe the attraction level of a transit node to passenger.Secondly, a mathematical model aimed at maximizing served populations per construction cost was proposed to optimize multihub location based on the selected candidate nodes.Dimitriou [29] applied hierarchical variational inequalities programming to optimal competitive freight network design problem.The key points are the particularities emerging in this case of transportation facilities strategic, operational planning, multiple game-theoretic, and equilibrium problems that are structured in cascade and in hierarchies.In order to study the integrated tactical planning of hub locations and the design of a frequency service network, Rothenbacher et al. [30] presented a pathbased model and solved it with a branch-and-price-and-cut algorithm.The constraints in the model included multiple transshipments of requests at hubs, transport time limits for requests, request splitting, and outsourcing possibilities.
From the above literatures, we can find the following.(i) Most of the researchers try to select the sites by using a transportation costs based optimization model, or travel time based optimization model, or both.They did not consider the service market demand, railway development level, integrated transportation convergence development level, and so on.(ii) Some researchers determine the node's level by considering the GDP, population, number of industrial enterprises, industrial added value, total retail sales of social consumer goods, highway mileage, inland navigation mileage, railway operation mileage, and so on, but they did not consider their influence on the related policies.(iii) Some approaches, methods, and techniques (e.g., AHP) require experience and knowledge from the evaluators and experts; sometimes the calculation results are subjective and inaccurate.In this paper we will formulate an IECM to solve the HSR express SSP.First, we will introduce the study areas and minimum research unit in China.Then we formulate the SSP evaluation indicators, including HSR express freight train service market demand, HSR development level, integrated transportation convergence development level, and the related policies.The collected data about the indicators will be used as the inputs for the evaluation.The improved entropy method will be applied to obtain the weight for each indicator.Compared with AHP, the improved entropy method is more objective and accurate.Then the cloud model will be applied to select and classify the HSR express freight service sites.

HSR Express Freight Train SSP Evaluation
Indicators: The Inputs (i) HSR Express Freight Train Service Market Demand.There are 6 subindicators of the HSR express freight train service market demand: gross domestic product (GDP), regional industrial added value, number of large enterprises, total import and export trade value in the region, total retail sales of social consumer goods, and regional trade volume.GDP reflects the result of all local resident units in a certain period of time.Different levels of areas with different GDP are suitable for the construction of different levels of logistics center, because of the close relationship between logistics demand and GDP.We use the regional industrial added value to formulate the huge role for industry and related businesses on the logistics needs and apply the number of large enterprises to measure the important impact on the operation of the logistics centers by the potential large number of customers.Because the size of customers is the main service object of the railway logistics centers, the larger customers scale is conducive for the construction of highlevel logistics centers.We use the total import and export trade value in the region to formulate the development of international trade, which is one of the most important evaluation indicators.We also use total retail sales of social consumer goods to formulate the national economy of the relevant industries to the social groups and residents supply sales of consumer goods, which reflects the final decision.Also, we apply the regional trade volume to reflect the demands of the transportation and logistics; the high-level logistics centers are suitable to establish in the regions with large-scale transport demands.

(iii) Integrated Transportation Convergence Development
Level.There are 3 subindicators: road network density, which reflects the convergence situation between the logistics centers and the road network; white goods throughput in the main water ports, which shows development level of the international logistics; and the cargo throughput in the main airports, which reflects the convergence situation between the logistics centers and the regional aviation level.
(iv) The Related Policies.There are 3 subindicators.The first is administrative level, which is evaluated on the scale of 9, 7, 5, 3, or 1.We use 9 to represent the municipality directly under the jurisdiction of the cities, 7 to represent the subprovincial cities, 5 to represent the general open cities, part of the provincial capital cities, or the border cities.We also use 3 to represent the prefecture-level cities and 1 for regional cities.The second subindicator is the level in the national logistics strategic planning.Usually, we take 3 to show the national first-level logistics center cities and to show the second-level central logistics cities, and the others take 0. The last subindicator is the level in the National Strategic Planning.If the related National Strategic Planning existed in the regional under evaluation, we take 1; otherwise, we take 0.
We collect initial data about the ,  = 1, 2, 3, . . .,  second-class indicator during the ,  = 1, 2, 3, . . .,  year; the initial data will be used as the inputs of the whole evaluation processes.For each second-class logistics area ,   , represents the collected value of the  second-class indicator during the  year.

The Translation-Corrected Improved
Entropy-Cloud Model: Formulation, Proof, and SSP Outputs

The Translation-Corrected Improved Entropy Weight
Method.In this paper, we use the translation-corrected entropy weight method to obtain the indicators' weight.The entropy weight method (EWM) is based on Shannon entropy, originally developed by Shannon and Weaver [31].Shannon entropy is a concept which is proposed as a measure of uncertainty in information and formulated in terms of probability theory.The concept of entropy is well suited to measuring the relative intensities of contrast criteria in order to represent the average intrinsic information transmitted for decision-making [32]; then the method is applied to describe the thermodynamics information systems [33].The EWM has already been applied to many research fields such as clinical neurophysiology [34], transport systems [35], environmental time series data, and conflict analysis [36,37].The uncertainty of signals in communication processes is called information entropy; the lower the information entropy, the higher the weight.The existing entropy weight method works as follows to calculate the weight.
The collected data about the subindicators mentioned in Section 2 can be formulated as a information decision matrix  = [  , ] × and then calculate the initial decision matrix as standardizing decision matrix  = [  , ] × .There are two types of standardized methods: when the indicators are benefit-type, the calculation for normalization can be expressed as formula (1); if the indicators belong to cost-type, the calculation for normalization can be expressed as formula (2).From Section 3, we can find that all the subindicators belong to benefit-type, so all subindicators will be normalized by using the following: where   , is the decision matrix after normalization, ( The information entropy for each evaluation criteria  is defined as follows: And the weight obtained from information entropy for each evaluation criteria  is expressed as follows: , ∀, , where 0 ≤  , ≤ 1 and ∑  =1  , = 1, ∀.From the above EWM calculation formulas we can find that if the initial subindicator data has small fluctuation amplitude, or the data remains unchanged during the operation periods, which means, ∀, for each , the value   , = 1, or   , = 0, and we use R to represent the condition mentioned above.If  satisfies R, then it shows this subindicator  has no impact on the total HSR express freight train selection process and the evaluation will fail.In order to solve this problem, we corrected the existing EWM by adding a translation data into formula (3), detailed as follows: , ∀, , where Δ  , is the translation value.After the translation correction in formula (7), some of the information contained in the existing data will be lost.In statistics, we can use the variance to measure the amount of information contained in a set of data.The variance of the data without translation will be larger than the corrected data.In order to maintain the integrity of the original data, the data information loss must be controlled in an acceptable range.For each evaluation criteria , we assign Information Loss Tolerance Coefficient (ILTC)  , to restrict the information loss artificially.We next give the detailed reasoning, proof, and remarks about the properties.
Assuming there is a set of nonnegative numbers,  1 ,  2 ,  3 , . . .,   , . . .,   , summation is  = ∑  =1   .If the translation value is a constant Δ, 0 < Δ < 1, then summation after translation is   = ∑  =1 (  + Δ) =  +  ⋅ Δ.After the calculation in formula (3) and formula (7) From formula (11) we can prove that the variance of the data without translation will be larger than the corrected data, which suggests that some information contained in the existing data is lost.We must control the initial data information loss in an acceptable range. , will work as follows: var According to formula (13), we can obtain the relationship among  , , , , Δ using quadratic equation root function: We establish a function to formulate the relationship among  , , , , Δ, where  , , ,  are independent variables.Δ belongs to dependent variable.The function is shown as follows: Taken the derivative of each variable, we can obtain the following three constant established formulas: The above three formulas (16) show that when the value of ,  , is fixed, the larger , the smaller Δ.When the value of ,  , is fixed, the larger , the larger Δ.When the value of ,  is fixed, the larger  , , the larger Δ.In this paper, for the evaluation criteria , we can control Δ  , in a reasonable range by setting  , in an acceptable range, artificially.After the process, the data information loss can be controlled in a reasonable range.The inherent defects of the entropy weight method can be restrained to a certain extent and make all indicators play a role in the evaluation processes.

The Cloud Theory and Cloud Model.
Cloud theory was first delimited and applied by Li et al. [38].Cloud model (CM) can synthetically describe the randomness and fuzziness of qualitative information and implement uncertain transformations between qualitative information and its quantitative instantiations based on probability statistics and fuzzy set theory.Usually the CM uses the forward cloud transformation (FCT) and the backward cloud transformation (BCT) to implement the cognitive transformations, which can reduce greatly loss of information.Now it has been used extensively in many fields such as natural language processing, data mining, decision analysis, intelligent control, and image processing.Li et al. [39] used CM to evaluate the smart distribution grid; the evaluation model included two parts: synthesized cloud and remarks cloud, and the final evaluation results are given in the form of probability distribution.Also, Fan et al. [40] applied the improved evaluation method based on CM for situation consistency within the battlefield of joint operations.The improved assessment approach is based on Cloud Gravity Centre Theory and uses Structure Entropy Weight to calculate the weight of the indexes.Furthermore, the CM has also been used for risk management, for example, Zhang et al. [41], Li et al. [42], and Zhang et al. [43], and route assessment; see Ma and Xu [44].Wu et al. [45] proposed a cloud-based decision framework to solve the waste-to-energy plant site selection problem (SSP).Deng et al. [46] combined Heuristic Gaussian cloud transformation and fuzzy forecasting method to solve the novel hybrid water quality time series prediction problem.Some of the recent papers have been also devoted to the development of CM theoretical research.Wang et al. [47] proposed a linguistic multicriteria group decision-making method based on three kinds of cloud aggregation operators, that is, cloud weighted arithmetic averaging operator, cloud-ordered weighted arithmetic averaging operator, and cloud hybrid arithmetic operator.Also, Wang et al. [48] used CM to solve the uncertain linguistic multicriteria group decision-making problem.Wang et al. [49] established a generic normal cloud model and then designed and compared some algorithms to solve the model, which included 2nd2nd-order generic normal cloud model (2nd2nd-GFCT), 2nd2nd-order generic forward normal cloud transformation algorithm (2nd2nd-GIBCT), 2nd2nd-order generic backward cloud transformation algorithm (2nd2nd-GBCT), pthpth-order generic forward normal cloud transformation algorithm (pthpth-GFCT), and the backward cloud transformation algorithm (pthpth-GBCT).Wu et al. [50] So, the evaluation process for each station based on CM works as follows.
Step 1.We use Expert Scoring Method to obtain the basic evaluation criteria cloud.
All the second-class logistics areas have three handling alternatives: renovate and develop as a LC, or LB, or the second-class logistics area has no HSR logistics business to handle (referred to as ).Each kind of alternative is formulated by a basic visualization evaluation cloud.We use ESM together with MATLAB to control the convergence speed and quality of expertise, and the final visualization evaluation clouds must meet the overall national strategic planning, such as the number limitation for each kind of station.The three basic evaluation clouds  LC ,  LB ,   will be used as criteria during the final evaluation process; please refer to formula ( 21)-( 24) to get the core calculation steps.Three basic evaluation clouds were visualized in Figure 2(b).
Step 2. For each second-class logistics area  and second-class indicator , use the obtained standardizing decision matrix  = [  , ] × to do the following calculation.Firstly, calculate the sample mean  ,  and sample variance ( , ) 2 , respectively.
Then we can obtain the Entropy Value 4.3.The SSP Outputs.For each second-class logistics area  under evaluation, to obtain the final evaluation cloud   , we use the weight  , to multiply the initial evaluation cloud  , .
Use MATLAB to obtain the visible   ,  LC ,  LB , and   in the same interface and compare the four figures: the evaluated second-class logistics area  belongs to one of the alternatives ( LC ,  LB , or   ) if   is closer to one of the figures.Also, we can use the following evaluation method to obtain the final result, adjunctively.
The detailed circuit to solve the HSR express train SSP is presented in Figure 3.

Case Study: Comparison Analysis and Results Discussions
We now give our case study by taking the Yangtze River Delta first-class logistics area in China as a background.
Yangtze River Delta is one of the most developed logistics area in China, with an area of 354,400 km 2 and a total population of 227 million.It consists of the three provinces of Anhui, Jiangsu, and Zhejiang, with Shanghai as the center of economic activities.In addition to the main city of Shanghai, there are a number of regional central cities, such as Hangzhou and Nanjing.Other cities such as Wenzhou, Ningbo, Wuxi, and Suzhou are also highly developed.With the completion of the Hefei-Nanjing Passenger Railway (the first HSR in the Yangtze River Delta) in 2008, 24 out of Outputs: 39 cities were accessible by HSR in the Yangtze River Delta by 2014 [51].The highly developed HSR networks not only reduce the passenger travel time between cities but also provide good foundation for the development of HSR express freight train services.In order to select the LCs and LBs, we divided the Yangtze River Delta into 13 second-class logistics areas, shown in Figure 1.For each second-class logistics area , we collected the initial data of the 15 subindicators from 2008 to 2015 and formed an information decision matrix  = [  , ] 8×15 .The initial matrix of the 13 second-class logistics areas will not be presented in this paper because of its huge scale; for detailed information please see China Railway Corporation and Southwest Jiaotong University [52]; all the used data was obtained from National Statistics Bureau of the People's Republic of China [53], Ministry of Transport of the People's Republic of China [54], and Ministry of Commerce of the People's Republic of China [55].Next we presented the weight for each subindicator of the second-class logistics areas based on improved entropy method and the HSR freight train service site selection results.

The Weight.
For the Yangtze River Delta logistics area, some subindicators such as administrative level and level in the national logistics strategic planning remained unchanged during 2008 to 2015, so these two subindicators need to be translation-corrected.In order to compare the difference of final weight under different ILTC  , , we took the Shanghai second-class logistics area as an example.We let  , = 0.2, 0.15, 0.1, 0.05, and 0.01, calculated the translation value according to formula (14), and obtained the value ranges.We used the maximum of Δ  , in this paper.The weight  , under each  , is presented in Table 1.
From Table 1 we can find that the smaller  , is, the closer the final weight  , is to the weight without translationcorrected ( , = 0).So, if we want to control the data information loss in a reasonable range, we should set  , in a sufficiently small range so that the inherent defects of the entropy weight method can be restrained to a certain extent and make all indicators play a role in the weight obtaining processes.Therefore, we set the ILTC  , = 0.01 for the other 12 second-class logistics areas, and the final weights for each subindicator of the 13 second-class logistics areas are as presented in Table 2.
From Table 2 we can find the following.(i) For the development of HSR logistics in the Yangtze River Delta area, the related policies are not the most important influencing factors.(ii) The key influencing factors are GDP, total retail sales of social consumer goods, HSR operational line length, and road network density, because the weight values of the four subindicators are larger than the others.the output of the 10 experts' evaluation and scoring.For each subindicator of the second-class logistics area, we can obtain an evaluation cloud.In this paper, we only presented (shown in Table 3) four second-class logistics areas' clouds (including 1 LC and 3 LBs) because of the data's huge scale.

The Selection
According to the formulas ( 26)-( 27), we multiplied the weight in Table 2 and evaluation clouds in Table 3 and obtained the final evaluation cloud   for each second-class logistics area .Next, we used MATLAB to make the clouds   visual and compared   with  LC ,  LB , and   (also we can use the formula 30).The final results are presented in Table 4 and Figure 4.
South Jiangsu (SJ) was selected as the HSRELC (red cloud in Figure 4) with Nanjing chosen as the location city.Shanghai (SH), North Zhejiang (NZ), and Central Anhui (CA) were selected as HSRELBs (yellow clouds in Figure 4), with Shanghai, Hangzhou, and Hefei chosen as the location cities, respectively.The other 9 second-class logistics areas have no HSR freight train service (blue clouds in Figure 4).
Meanwhile, we can apply the same approach to obtain the HSR freight train service sites for other 7 first-class logistics areas, including Pearl River Delta, North China, Northeast China, Bohai logistics area, Central China, Southwest China, and Northwest China.The 8 LCs, 23 LBs, and the collected HSR railway lines formed the basic freight train service network in China as shown in Figure 5.

Conclusion and Further Study Works
The High-Speed Railway (HSR) express freight train service belongs to an emerging transport product, which is the most important approach for China railway to obtain more   There are two parts for the sites selection processes.First, we obtained the weight of each subindicator based on translation-corrected improved entropy weight method.The mathematical proof and the Shanghai case study results showed that the smaller  , is, the closer final weight  , is to the weight without translation correction.So, if we want to control the data information loss in a reasonable range, we should set the ILTC in a sufficiently small range so that the inherent defects of the entropy weight method can be restrained to a certain extent and make all indicators play a role in the weight obtaining processes.Second, we applied the cloud model (CM) to obtain the evaluation cloud for each subindicator.We used the weight to multiply the initial evaluation cloud, and the resulting outputs were presented visually for easy comparison and final selection.
We next gave our case study taking the Yangtze River Delta first-class logistics area in China as a background.In order to select the LCs and LBs, we divided the Yangtze River Delta into 13 second-class logistics areas.For each second-class logistics area , we collected the initial data of the 15 subindicators from 2008 to 2015.We selected the South Jiangsu as the LC with Nanjing chosen as the location city.Shanghai, North Zhejiang, and Central Anhui were selected as LBs with Shanghai, Hangzhou, and Hefei chosen as the location cities, respectively.The other 9 second-class logistics areas have no HSR freight train service.Meanwhile, we applied the same approach to obtain the HSR freight train service sites for other 7 first-class logistics areas resulting in 8 LCs, 23 LBs, and the collective HSR railway lines forming the basic freight train service network in China.
The proposed new site selection approach is also suitable for other related problems, but the evaluation indicators may be different.So, how to apply this approach into actual engineering practice and how to determine the ILTC, as well as the value of Δ  , reasonably in the practice, are further study works.

Figure 1 :
Figure 1: Layout of 8 first-class logistics areas and second-class logistics areas (using Yangtze River Delta as an example).
proposed a new model, which combined the normal cloud model (CM) and Fruit Fly Optimization Algorithm (FOA), to solve the problem.According to the previous research results, the digital signature of clouds reflects the quantitative characteristics about the qualitative concept.For each second-class logistics area  and each second-class indicator ,  = 1, 2, . . ., ,  ∈  + , we use three kinds of values to formulate a cloud, including Expectation Value  ,  , Entropy Value  ,  , and Super Entropy Value  ,  . ,  represents the average point of the cloud drops group. ,  shows the fuzziness and randomness in the qualitative concept and represents the acceptable range of concept in the domain. ,  is the uncertainty measurement of Entropy Value  ,  or, more specifically, the entropy of the Entropy Value  ,  , which represents the condensation degree of the cloud drops group.Smaller  ,  value implies a better condensation degree of the cloud drops group, lower measures of dispersion, and thinner thickness of the cloud.The three basic values used to formulate the CM are shown in Figure 2(a).If we use  as the quantitative characteristics with exact value,  , is the qualitative concept; then  , =  (

Figure 2 :
Figure 2: Basic elements of CM applied in this paper.

Figure 3 :
Figure 3: Detailed circuit to solve the HSR express train SSP.
Results.To satisfy the basic national strategic planning of the HSR freight train service sites selection, we invited 10 experts (including 5 Professors and 1 Associate Professor from School of Transportation and Logistics, Southwest Jiaotong University, and 4 Senior Engineers from China Railway Siyuan Survey and Design Group Co. Ltd.; all experts' specialization is railway operation network planning and design) to evaluate and score for the three kinds of criteria.After transferring the scoring value into evaluation values  ,  ,  ,  , and  ,  , we used MATLAB to make the clouds visual.The final evaluation clouds criteria are  LC = (0.93,0.03, 0.006),  LB = (0.85,0.04, 0.007),   = (0.33,0.20, 0.008), shown in Figure 2(b), which is

Figure 4 :
Figure 4: The visual results of the 13 second-class logistics areas in the Yangtze River Delta area.

Figure 5 :
Figure 5: The HSR freight train service network in China.
freight train SSP evaluation indicators system, we divide the whole China into 8 first-class logistics areas, including Pearl River Delta (Guangxi, Guangdong and Hainan Province), Yangtze River Delta (Anhui, Jiangsu, Shanghai and Zhejiang Province), North China (Inner Mongolia and Shanxi Province), Northeast China (Heilongjiang and Jilin Province), Bohai logistics area (Beijing, Tianjin, Hebei, Liaoning, and Shandong Province), Central China (Henan, Hubei, Hunan, Jiangxi, and Fujian Province), Southwest China (Xizang, Sichuan, Yunan, Chongqing, and Guizhou Province), and Northwest China (Xinjiang, Qinghai, Gansu, Ningxia, and Shanxi Province).We then divide each first-class logistics area into several second-class logistics areas; see Figure1.In this paper we focus on how to select and classify the HSR express freight train service station in the second-class logistics areas scale.For each second-class logistics area , we establish the HSR express freight train SSP evaluation indicators system, which includes 4 first-class indicators and 15 second-class indicators as described below.
China has 34 provinces, with different economic development and HSR operated lines length.In order to establish the HSR express , ) max is the maximum value of   , , and (  , ) min is the minimum value of   , .If (  , ) max = (  , ) min , we let   , = 0 and then calculate the probability of each   , as follows: , ∀, , .

Table 1 :
The weight  , under each  , (the Shanghai second-class logistics area case).

Table 2 :
Weights for each subindicator of the 13 second-class logistics areas (the Yangtze River Delta case).
SH: Shanghai; NZ: North Zhejiang; for other abbreviations please see Figure1.