Network Planning Method for Capacitated Metro-Based Underground Logistics System

Underground logistics system (ULS) tends to alleviate traffic congestion, increase city logistics efficiency, mitigate the negative effects of traditional logistics processes, and improve the sustainability of urban areas. However, the relatively high cost and risk of underground construction are serious obstacles to implementing ULS. Integrating ULS into modern metro system (M-ULS) is considered to be feasible and efficient to solve this problem. ,is paper aims at developing a metro system-based ULS network planning method. First, an evaluation model of underground freight volume was proposed considering service capacity, freight flow, and regional accessibility. Second, a set of mixed integer programming model was developed to solve the problem of optimal nodes’ location-allocation (LAP) in the network. ,en, a hybrid algorithm was designed with a combination of E-TOPSIS, exact algorithm, and heuristic algorithm. Finally, two lines of Nanjing Metro were selected as a case to validate the proposed planning method. ,e results showed that the new system can significantly reduce the construction costs of ULS and alleviate traffic congestion. Moreover, the potential of metro stations and underground tunnels can be fully exploited to achieve higher logistics benefits.


Introduction
Traffic congestion is one of the "dilemmas" encountered by metropolises around the world.
e fifth type of transportation and supply system, "underground logistics system" (ULS), has been developed to address the energy consumption, environmental pollution, and safety issues [1].Various forms of ULS with different locomotives or driving methods have been developed and applied, such as CargoCap underground freight transportation system in Germany [2], OLS-ASH project in Netherland [3], and Pipe §net vacuum freight capsule in Italy [4][5][6][7].ULS is widely recognized as a great innovation in the logistics industry because of its huge benefits to the city.It can quickly transport goods through underground pipes or tunnels within cities or between cities [8,9].However, it is undeniable that the high cost, long construction period, and high-risk of underground projects [10][11][12] are serious obstacles to ULS development.Reducing the construction cost is critical to improve the uptake of ULS [13].
Dampier and Marinov proposed a new conception of "integrating underground freight transport into modern metro system" (metro-based ULS, designated as M-ULS), which offered a feasible way with limited budget [14].As a special form and promising alternative of ULS network, M-ULS has become one of the hottest research topics.Especially, setting up the intermodal network aims at maximizing the efficiency of cargo transport without interfering with the passenger transport system.e organizational process of this new system is complicated.It involves the linkage and optimization of the dual objectives of the passenger and cargo transportation system.e cargo needs to be sent to the demand points via reasonable metro stations and the distribution system on the ground.erefore, the global optimization of the M-ULS efficiency can be classified as a special type of capacitated multistage facility location problem (CMFLP).e aim of this paper is thus to design a multilevel node location-allocation solution for M-ULS through mathematical modeling and simulation.
Research on ULS has been ongoing for nearly 30 years, but M-ULS is still an emerging area.Existing studies verified the feasibility of M-ULS and proposed concept designs.M-ULS is defined as an intermodal mechanism of mixed passengers and freight, in which batch driving or attachment carriages are enabled by utilizing metro railways [15].eir advantage not only lies on the reservation of traditional ULS benefits in alleviating traffic congestion, improving urban environment and freight efficiency, but also contributes to a huge saving of the underground space, reducing construction costs and period.
e passenger transport efficiency would not be interfered through rational technology design and operation management [15].In Tokyo, new subway-integrated ULS received public endorsement in a pilot project [16].In Newcastle upon Tyne, metro network was used for the collection of small-sized to medium-sized parcels, low-density high-value goods, and recyclable material [17].In Korea, the freight volume of Seoul Metro 50 platforms was predicted to validate the reduction of social cost and environmental problem [18].Toronto subway has also been evaluated on how such a modal shift appealed to cargo shippers [19].All in all, real-world simulations were widely developed in recent years [20][21][22].
However, few articles with the topics on the planning method of M-ULS network have been published.Hörl et al. and Masson et al. conducted qualitative analyses of node location problem of single-line metro freight transport [23,24].Fangting et al. proposed a location-routing model of multiline transfer problem of M-ULS [25], which aimed at minimizing delivery time to improve the efficiency of metro freight.
e optimization of the single-layer network topology was carried out; however, the influence of urban freight flow distribution on the system service level was not involved.Ambrosino and Scutella [26] proposed a four-level complex network distribution structure including facility location, warehousing, transshipment, and strategic decisions.
e findings have contribution for analyzing the CMFLP of M-ULS, although it does not directly discuss M-ULS.M-ULS itself is a complex system.Functional classification of different types of facilities is necessary to meet the needs of intermodal transport.In order to pursue higher logistics benefits, a systematic evaluation of the urban freight transport demand is required due to strict capacity constraints of M-ULS.e freight flow direction and node location are two key indicators of M-ULS network layout.
ereby, this paper proposed a hierarchical mixed integer programming model (MIPM) combined with comprehensive evaluation that could reflect the different relationships of nodes.Compared to the modeling benchmark of general facility location problem [27][28][29], three additional important factors had been considered: (i) Constraint of metro network.Compared to the traditional logistics LAP model, M-ULS distribution path, which is restricted by the existing metro network, cannot be freely selected and adjusted immediately.(ii) New system service area.e definition of new system service objects has to be integrated into the optimization process because maximizing the limited system service capability is a key to M-ULS planning.(iii) Mixed transport organization is complicated.e design principle of existing metro network is passenger oriented, and the passenger transportation efficiency cannot be influenced by the new system.Improving the service performance of the new system requires redesigning an organizational process of passenger and freight mixed transportation, and the global optimization of freight transportation process.
e contributions of this study mainly lie in three aspects: (1) an effective organizational planning method for M-ULS was proposed; (2) a mixed integer programming model with multiple constraints was developed to solve the capacitated metro-based underground logistics system location-allocation problem (CM-ULSLAP); (3) a hybrid heuristic algorithm was developed for the simulation of realworld.Two lines of Nanjing Metro were selected as a case to validate the proposed planning method.e findings can promote the implementation of ULS.In addition, this study offers new ideas and quantitative optimization approaches for linear or network-based complex engineering projects.
e supplied case study provides additional insights for research and practices.

Problem Description
e CM-ULSLAP can be described as follows: there are m operating metro lines, n urban distribution centers (designated as UDCs), and k demand points.Firstly, cargo from the open UDCs which provide ULS services is imported into underground network through the nearby metro stations (designated as ULS terminal).en, dedicated cargo trains travel between ULS terminals and different underground depots (designated as UDs).ose UDs refer to the retrofitted metro stations, in which UD-oriented trains are bypassed to the parallel freight platforms before entering the passenger platforms.In this process, the frequency of freight fleet keeps interlaced with metro, and existing transfer points can be activated as underground hubs (designated as UHs) to support underground transshipment.Finally, freight units floating up from UDs are allocated to the corresponding ground terminals (designated as GTs) by ground electric vehicles.Cargo is delivered to the customers who are within the service scope of M-ULS after automated sorting and packaging at GTs. e reversed procedures are presented in mirror image.Figure 1 demonstrates the organizational breakdown structures with planning issues of M-ULS.
e initial data that can be used for the metro freight planning usually focus on several incomplete customers characteristics, such as location, regular delivery route, and handover period record.Other important input data are the demand for travel, which is commonly formulated as the origin-destination (OD) matrix [30].e freight OD matrix in this study is extracted based on the counts of travels multiplied by order capacity from one area to another.On this basis, the alternative solution aims to improve the current customer service with low cost and high efficiency.
e high efficiency of M-ULS depends on the coordination of its own technical integration and freight flow distribution.In order to get this, the joint freight function should not lead 2 Advances in Civil Engineering to any degradation or disruption of passenger service [15].As a result, goods and passengers must be separated.Meanwhile, the efficiency loss caused by intermodal transport can also be compensated by a highly integrated urban logistics benchmark.Taking into account the regional development density and spatial impact on the surrounding environment, in M-ULS, freight flows cannot be detained too much after reaching corresponding underground destinations.Instead, they should be transferred directly to the newly built ground terminals for processing.By setting such a transit layer between customers and underground distribution network, the transport efficiency and system service capacity can be largely improved by dividing the M-ULS from the traditional logistics systems.Based on the planning principles discussed above, the node facilities are divided into the following four layers:

ULS Terminal.
e ULS terminal plays a role in docking UDC with underground network.e freight units which have been well sorted and packed in UDC will be simply coded at the ULS terminal to each exact ground terminal or underground depot.In this paper, the location of the ULS terminal was determined by calculating the weighted distance to nearby urban distribution centers.

Underground Depot (UD).
e UD is defined as an enhanced metro station which has the function of cargo handling and connects the underground logistics network with upper ground terminals.e number and distribution of UD directly affect the system service.In this paper, a MILP method with corresponding heuristic algorithm was introduced to search the optimal location-allocation of UD-GT.

Ground Terminal (GT).
GT is one of the optimization targets, which is affiliated with UD.GT is also a key part of the entire M-ULS supply chain because it is directly oriented to customers whose delivery requirements are time efficient.
e cargo handling capacity of GT is the major handicap in the performance of system service.In this paper, the location of ground terminal was transformed into the capacitated set covering problem (CSCP).And a set coverage model for GT location is established.

Demand Point.
Demand point is also one of the optimization targets.Because of the limited transportation capacity of the metro system which is far from meeting cargo supply for all customers, a screening evaluation of demand points is essential for maximizing logistics benefits and service scope.
In order to simplify the study, the following assumptions for CM-ULSLAP are proposed: (1) suburban UDC reasonably allocates goods into the underground network, and the freight handling capacity of ULS terminal is unlimited; (2) each candidate GT can only correspond to one UD, and the source of freight flows is unique; (3) taking a full carriage of dedicated trains as a freight unit, the less than truckload (LTL) transportation are not allowed in the underground distribution network; (4) the processing time and transfer cost at underground hub are not considered, and all UHs are activated by default; (5) intermodal transport does not involve the urban road network layout and additional expenses caused by traffic congestion; (6) M-ULS does not reschedule the original timetable of metro system, but underground freight flow frequency can be adjusted based on the traffic volume of passengers; (7) construction and storage condition at UD are not considered; and (8) only part of customers' attributes are available.

E-TOPSIS for Evaluation of Underground
Freight Volume e entropy weighted TOPSIS method (E-TOPSIS) is applied in the weighting process of the evaluation indicator Advances in Civil Engineering system [31].e information entropy calculated in the indicator matrix is used to weight the dispersion degree of demand characteristics, all that can avoid the influence of human subjective factors [32].e final points of the evaluation objects in multiobjective decision-making are gotten by the TOPSIS method which can solve the higher requirement for the sample.
is paper used unity of freight flows (UFFs), regional accessibility (RA), and source distance (SD) as three dimensions to reflect UDC's preference for the characteristics of customers' demand.erein, variables UFF and RA were obtained from Nathanail et al. [33], who maintained that the urban traffic alleviation and the development of integrated transport to reduce the scattered freight delivery were two of the most remarkable characteristics of the sustainable city logistics.In our study, we took these two quantitative features as the functional requirements of the metro freight system, where UFF was formulated as the accumulated OD from different directions according to Masson et al. [24] and RA was formulated as the accumulated delivery period from different directions divided by average path length according to Montoya-Torres et al. [34].SD was formulated as a variable that supplemented UFF with distance influence, which involved the data of spatial distance and freight OD from different directions.ese 3 variables share the initial data mentioned above.
(i) Unity of freight flows.Since UD does not have the function of parcel disassembly, logistics benefits are promoted along with the growth of underground mileages [35].erefore, M-ULS should provide centralized supplies as much as possible to areas with huge demand: where Q is the total amount of freight picked up and delivered per day at demand point j.
(ii) Regional accessibility.e efficiency of urban logistics is largely influenced by traffic congestion [36][37][38].By removing goods from above to underground, the freight mobility of some regions with poor accessibility can be effectively improved [39].e index of regional accessibility can be characterized by the time-road response function: where t ij is the average delivery time form UDC i to demand point j, S is the length of the delivery route, and τ ij ∈ (0, 1) is the time sensitive factor.
(iii) Source distance.Underground logistics can generate considerable internal and external benefits [40].M-ULS should take full advantage of the length of metro lines.Putting goods which are far from the source (ULS terminal) but in huge demand into the system will lead to the overall increase of logistics values.e source distance index can be portrayed as where ‖d ij ‖ is the Euclidean distance between UDC i and demand point j, c i (x j ) is the freight correction factor of the sample, which represents the proportion of the freight volume from direction i, and α j is an amplification coefficient.An evaluation model for the priority level of M-ULS service is constructed.
(1) e initial decision matrix.Assuming there are m evaluation indicators and evaluation indicator sets have n subsets, the evaluation value of indicator i in subset j is x ij .
e decision matrix of all subsets is (2) Standardization of the decision matrix.In order to solve the uniform of indicators' units, this paper makes the standardization on all indicators.Take the best value of indicator i as x best i , the worst as x worst i .e efficacy coefficient α is introduced to obtain the normalized data matrix B: (3) Determine the indicators weight.e indicator weights are determined by the expert evaluation method in the TOPSIS.e results are quite subjective.erefore, this paper adopts the information entropy method to determine the indicators weight.
e i denotes the entropy of indicator i in the standardized decision matrix: e dispersity of the evaluation value of indicator i can be represented as follows: (5) Calculating the relative distance.e selection scheme of candidate GT for this phase is calculated with the relative proximity ρ j of positive and negative ideal solutions under weighted matrix K: where d − j denotes the Euclidean distance between k ij and negative ideal solution: For the candidate location scheme, ρ j is generally between 0 and 1. e closer to 1, the more appropriate the corresponding scheme.Finally, demand points in the newly built GT are selected to accept the M-ULS service according to the value of ρ j listing in descending order, until one of the nodes or underground sections is overloaded.Underground freight OD can be uniformly formulated as

Parameter Definition
(1) Define V C , V S , and V P as sets of metro station, candidate UD, and GT, respectively.
e cost parameters: c 1 , integrated unit price for underground transport; c 2 � β • c 1 , comprehensive price of ground transport, β is the freight price coefficient, which reflects the pricing gap between two modes; and c 3 , construction cost of UD (depreciated to daily).
e distance parameters: d ij , the distance between UDC j and metro station i; r ik , the shortest path for goods routed to the UD k, in which underground transfer is included; and Z x , andZ i , coordinates of GT and metro station.
e freight parameters: Q ij , the amount of goods sent by UDC j to certain metro line where the UD i is located; Q k , the amount of goods picked up by the UD k per day; u j max , capacity of metro line j.
M i : 1, if metro station i is selected as UD; 0, otherwise.P xk : 1, if GT x is visited by the vehicles from UD k; 0, otherwise.
H ij : 1, if metro station i is selected as underground access of goods from UDC j; 0, otherwise.

Set Coverage Model for Ground Terminal.
Sets A(j) of all demand points covered by ground terminal a j and sets B(i) of all GTs that cover demand point b i are given.Let Y ij be the distribution coefficient of cargo throughput for a j to b i and X j be the decision of whether to locate ground terminal at a j : min  j,a j ∈V * R X j . ( subject to j,a j ∈B(i) e objective ( 14) is to search the minimum number of GT under the current demand scenario; constraints (15) and ( 16) indicate the maximum service scope and cargo handling capacity of the ground terminal; and constraint (17) specifies the sum of freight demand b i shared by each GT should be equal to the amount of OD obtained in (13).
For unified presentation, the above model is written in the form of matrix inequality: K is a 2n * n(n + 1) dimensional matrix.In the first n line, elements of column 1 to n equal 0; n + 1 to 2n is a diagonal matrix with principal diagonal D 11 , . . ., D n1 , and so forth.In the line n + 1, column 1 to n principal diagonal is −Q j max , while column n + 1 to 2n 2 has a single line element R at the same time compose the feasible space of the GT location under the maximum service capability.

Location-Allocation Model for Underground
subject to Formula ( 20) is the objective function, which consists of transportation cost and node construction cost for the intermodal logistics system.Constraint ( 21) reflects the capacity restrictions of underground freight handling at UD; constraint (22) ensures that the capacity of metro lines are not violated; constraint (23) stipulates that the ULS terminal on a single metro line is unique, but the UDC can be assigned to different ULS terminals based on reality; constraint (24) guarantees that UD is connected to at least one ground terminal of which ownership is unique; constraint (25) specifies that the number of metro stations that are selected as underground depots equals y; and constraints ( 26)-( 28) are the range of relevant decision variables.

Solution Method
CM-ULSLAP is an NP-hard problem in which the customer layer (demand point), middle layer (GT), underground logistics layers (UD and UH), and supply layer (UDC) are integrated into a complex intermodal network. is problem is computationally intractable due to the following reasons.First, there is not a single layer of node information that can be fully invoked.Different intermodal strategies result in the conflict between customer demands and economic considerations.e evaluation procedure is hard to carry out.On the contrary, the possible paths for each OD pairs are so huge that it is computationally expensive and unsustainable to enumerate all of them for each configuration involving transport mode choice and corresponding routing transfer.Once the supply relationship and the generation strategy of candidate GT are determined, CM-ULSLAP reduces to the bi-level capacitated intermodal network design problem (BCIND).For calculation efficiency and stability, the access routes from UDCs to ULS terminals were determined in advance to cut off the supply layer.en the entropy weighted TOPSIS method (E-TOPSIS) was applied to evaluate the demand characteristics in model (P1).Next, an exact algorithm for set coverage 6 Advances in Civil Engineering model is served for obtaining feasible solutions in model (P2), thereby cutting off the customer layer by replacing abundant demand points with candidate GT sets.Immune Clone Selection Algorithm (ICSA) was applied to optimize the location of UD and allocation of deterministic GT which following the minimum objective cost in model (P3).An adaptive mechanism of underground route navigation was also embedded in ICSA where fleet flows always follow the shortest path that is currently available until the certain section is blocked.e ICSA was inspired by the principle of clonal selection of antibodies in the biological immune system [41].It was used to explain the basic characteristics of adaptive immune response to antigen stimulation.e ICSA has been widely used in solving problems such as combinatorial optimization, intelligent optimization, and production scheduling due to its powerful data search capabilities [42][43][44].However, the convergence rate of original ICSA is slower, immune probability and cloning probability are relatively fixed, and the degree of change is relatively low when solving complex problems.In this paper, we introduced a self-adaptive mutation operation based on normal distribution to achieve a uniform and dynamic global high-level ambiguity for each antibody that meets the mutation rate for the characteristics of multiple UD feasible solutions and frequent information feedback at each stage.Probability variation enhances the randomness and stability of the search and can effectively prevent the solution from falling into local optimum.e main operation of the heuristic algorithm is shown in Figure 3.
e solution steps of GT location are listed in Algorithm 1.
Step 1: Generate Initial Population.e initial antibodies were obtained from the memory cells that were generated randomly from the feasible solution space.
e coding of antibodies includes the UD open strategies and the location information of newly built GT which were fed back by E-TOPSIS.e initial antibody population was recorded as P 0 (RC), which was composed by RC number randomly generated antibodies.
Step 2: Diversity Evaluation.Select k individuals with the highest fitness value and t individuals with the lowest affinity from the parent population P n (RC) to compose the solution vector, where k � RC × 20% and t � RC × 10%.Affinity was used to indicate the matching degree between antibody and antigen, namely, the optimal satisfaction between the candidate GT sets and open UD in objective (19).Lehmer mean [45] was recommended to the affinity determination.e affinity of the uth solution vector can be described as Step 3: Clone Operation.e clone ratio was calculated by the evaluation of affinity and similarity between antigen and other antibodies.
e cloned population Z n (N C ) is produced by replicating k + t ± d number of selected antibodies.Specifically, the total number of clones generated for all selected antibodies is Step 4: Gene Mutation.Select a group of clones from Z n (N RC ) and perform Gaussian mutation on them to obtain population S n .e mutation rate adopts an adaptive strategy [46] which is related to the adaptation value f k of the antibody.
is mutation operation can be described as m j � normrnd (m j , σ, 1, 1), where m j is the jth attribute of the clone, and normrnd is a normal-distribution random number with a mean of m j and a standard deviation of σ. e σ domain of the antibody is adaptively adjusted to σ � ω • I u /f k following its adaptation value and affinity.
Step 5: Immune Selection.e memory population M is consisted of a batch of antibodies with the highest f k from the set S n ∪ Z n .en select 30% k number of antibodies with the highest fitness value from M to update the equal number of antibodies with the lowest fitness value in the population P n to form the progeny population P n+1 .
Step 6: Node Search Reset.After reaching the maximum number of iterations N, check I (n+1) . en update the objective cost J b (i).
Let y � y + 1, returning to step 2; if not exists, step forward.
Step 7. e algorithm terminates when reaching iteration number N. Extract the most-affinity antibodies with their attributions as the preferred configuration of M-ULS.

Background and Parameter Setting.
e Nanjing city is located at latitude N32 ∘ 02 ′ and longitude E118 ∘ 46 ′ in East China.is city is situated on the Yangtze river bank.Based on the 2016 census, Nanjing had a population of 8.27 million.Express delivery volume of Nanjing had reached 630 million parcels in 2017 (ASKCI consulting).e total amount of social logistics is nearly 3 trillion RMB, with a 7.96% annual growth rate (Nanjing Municipal Bureau of Commerce).ree major transportation hubs, Nanjing Railway Station, Nanjing South Railway Station, and Lukou Airport, run through the city of which freight throughput capacity breaks 800,000 tons per year.Nowadays, there are 9 metro lines and 164 stations in Nanjing.e total length of metro lines is 347 kilometers, the 7th longest in the world.
For the sake of studying the performance of metro freight in the real-world case, Nanjing Metro was taken as a case to schedule the optimal configuration layout under minimum objective cost.A scenario of 2 open UDC, Line 1 and Line 2 with an activated UH and 31 stations, was Advances in Civil Engineering involved based on the real-world location.As shown in Figure 4, there are 274 demand points distributed in a range of 201 square kilometers from Xinjiekou to Nanjing South Railway Station.Relevant data and parameter values were obtained by the Monte Carlo simulation method, as shown in Table 1.e underground freight volume of each UDC was limited to 10,000 tons per day to be in consistent with the transport capacity of underground tunnels.Cargos from two metro lines can be transferred at Xinjiekou Station where the function of UD is invalid.Advances in Civil Engineering end if (8) end for (9) for j ∈ A(j) do (10) assign to j′ subsequently according to length B(i) (11) while Q j ′ max � 0 or A(j′) � ∅ do (12) end while (13) if true then (23) return to line 4 (24) end if (25)  Advances in Civil Engineering e proposed solution algorithm was coded in MATLAB R2016, and all experiments were run on a desktop computer with Windows 10, Intel Core I7-7700K 4.2 GHz processor, and 32 GB of RAM. e initial population size p � 50 and maximum iteration number G max � 100.

Simulation Results.
e E-TOPSIS method was applied to illuminate the range of system services at first.e number of qualified demand points is reduced to 228, in which the relative proximity of the most urgent area maxρ j equals 0.1916, and the critical proximity ρ j � 0.1812.Metro Line 1 saturates faster than the Line 2 with the cargo throughput of 9,747 tons per day and 7,758 tons per day, respectively.Algorithm 1 gave a set of feasible GT location demonstrated in Figure 5. e number of ground terminal that satisfies the constraints is determined as fifteen.en the optimization results from ICSA indicate that the minimum cost of the M-ULS-integrated logistics system is 185,900 RMB per day, while the transport cost of traditional ground logistics is 198,700 RMB per day.UDs were selected to locate on the following six stations: Xinmofanmalu, Andemen, Hanzhongmen, Daxinggong, Minggugong, and Jiqingmendajie.Figure 6 depicts the optimal open strategy of UD and corresponding allocation.Minggugong Station captured the busiest UD, with the daily freight volume of 3914.9 tons, whereas Andemen Station had the lowest freight volume of 1483.9 tons with single affiliated GT.
e specific node information is shown in Table 2.
Table 3 depicts the comparison of traffic volume from 2 UDCs to 15 GTs under traditional point-to-point distribution mode and M-ULS, where traffic volume (TV) is defined as freight OD quantity multiplied by the transport distance [47].It can be seen that the ground traffic volume (GTV) of GT-4 has been reduced by 69.32%, which benefits the most from underground logistics.Most GTVs got a different degree of alleviation between 30% and 65%.Note that GT-15 in Figure 5 is relatively close to the freight sources, and the GTV has increased by 36.71% after adopting M-ULS.erefore, the traditional point-to-point model is more economical for nodes near such UDCs.
When it came to the underground traffic volume (UTV) in M-ULS, all GTs exceeded 50% of the underground logistics ratio, where the maximum reaching 75.71% at GT-5.
e average above and underground mileage at Xinmofanmalu accounted for 34.03% and 65.97%, respectively.Table 4 demonstrates the overall variation of traffic volume.Benefit from the higher utilization of underground distribution, less than half of the ground delivery is required.For this part of goods that need to be delivered to the downtown areas, the original long-distance highway travel is decomposed into several subtransport processes above and underground.As a result, fewer freight vehicles are needed.
e mobility on urban traffic loops will be significantly improved and some congested road sections will be effectively alleviated from traffic pressure.10 Advances in Civil Engineering UD-opened n c .Note that the total cost of the intermodal system dropped slightly along with the increase of β and then reached the same level with traditional ground logistics when β equaled 2.52.Accordingly, there is a wide range of pricing options for underground freight transport, but the ratio of the discreet unit price of M-ULS to traditional ground logistics should be less than 1 : 2.5 in order to achieve economic advantages.On the contrary, the total cost of intermodal system ascended along with the increase of open UD, while the transport cost dropped due to the advance in network connectivity.Both trends were relatively stable.erefore, in M-ULS-integrated logistics system, optimal cost is not sensitive to the number of UD opened.It is considerable to set up additional underground depots so that the capacity and robustness system service can be further augmented.

Conclusions
is paper developed a multistage model with a combined algorithm based on hierarchical optimization, which aimed to deal with optimal nodes location-allocation problem in the cooperative operation of M-ULS.e applicability of the model was proven by selecting the appropriate metro lines and stations Advances in Civil Engineering for problem solving.e critical issue of developing a M-ULS network is to establish an operational framework and optimization method for multiple nodes.ree types of nodes were defined as underground depot (UD), underground hub (UH), and ground terminal (GT), in accordance with the current state and planning of the city, logistics, and metro network.e M-ULS is in limited capacity, which has a multiple network hierarchy and multiple sources of freight.erefore, it is necessary to screen the demand points and integrate them into the model for global optimization.A set of mixed integer programming model combined with comprehensive evaluation was developed as a more rigorous and convenient analysis approach to solve CM-ULSLAP.e model was embedded with three evaluation variables, namely, UFF, RA, and SD.A hybrid algorithm composed of E-TOPSIS, exact algorithm, and heuristic algorithm was designed to obtain the solution of the model.A real-world simulation based on Nanjing Metro was carried out to verify the effectiveness of proposed models and algorithms.
It was proved that the potential of metro stations and underground tunnels can be fully exploited to achieve higher logistics benefits.According to the results, optimized M-ULS contributed to reducing the overall ground traffic volume by 49.09% and the daily comprehensive cost by 6.44% within the service scope.Served by the new logistics system, 62.47% of the total freight was completed by ULS, and the remaining 37.53% was delivered to all demand points which were delivered by road distribution.Furthermore, compared to traditional pointto-point logistics mode, the price gap of M-ULS optimal cost has mostly remained stable, regardless of the number of UDs. e realization of these advantages neither affects the normal operation of the metro system nor further increases the amount of underground space development.e other comprehensive benefits for urban transport, environment, and society are obvious.is paper not only proposes a well-designed locationallocation model for joint operation of ULS and metro but provides a possible method for the collaborative research of different urban infrastructures.Moreover, it can provide novel ideas for the sustainable urban development.
Despite the above contributions, this study has several noteworthy limitations.First, the proposed model was exclusively applied to cross-shaped or #-shaped radial metro layout.
e loop network will have some special problems that cannot be ignored, for example, the trade-off among transshipment efficiency, routing selection, and budget.Second, the proposed deterministic system in which the UHs were always opened was actually not economical.Although the new problem always starts from certainty, future research needs to explore the uncertainty problems.Robust optimization approach can be used to study the UH scheduling strategy under conditions of stochastic disruption and demand uncertainty.Additionally, research on network planning or underground space can be further extended by considering the operational conflicts between M-ULS and metro system. is study was supported by the National Natural Science Foundation of China (71631007 and 71601095).12 Advances in Civil Engineering

Figure 1 :
Figure 1: Operational modules and issues of M-ULS-integrated urban logistics.
Finally, the solution space of combined models has been reduced from five dimensions M * N * K * |V C | * |V P | to three dimensions M * |V C | * |V ″ P | which contributes a lot to the computational efficiency.e relationships and decomposition of models are shown in Figure 2.

end while ALGORITHM 1 :Figure 4 :
Figure 4: Sketch map of metro lines and demand areas.
Figure 7  indicates the variation of optimal cost with freight price coefficient β and number of

15 Figure 5 :
Figure 5: Service scope of M-ULS and GT location.

Figure 7 :
Figure 7: Variation of objective cost against the change of β and n c .
bigger and indicator i is more important.If T ij is relatively concentrated, indicator i is less important.If all T ij are equal and the distribution is absolutely concentrated, indicator i is invalid.

Table 1 :
Parameters of metro freight.

Table 2 :
Summary of optimal node information.

Table 3 :
Comparison of traffic volume between point-to-point mode and M-ULS.

Table 4 :
Overall variation of traffic volume.