Optimization of Container Shipping Network Reconfiguration under RCEP

: Due to its advantages of large transportation volume and low transportation cost, container shipping has become an important transportation mode in current international trade. The recovery of the shipping industry in 2021 and the signing of RCEP make the reconfiguration and optimization of the container shipping network a very important task at present. The network service capability is an important factor affecting the container shipping network. Based on the complex network theory, the coefficients of port location, the importance of distance and route among ports are taken as the service capacity reconfiguration coefficients of the global container shipping network. A max-min mixed integer model is established for global container shipping network re-configuration. A sort of communication-reducing conjugate gradient method based on Krylov Sub-space (CR-CG-KS) is proposed to reduce the reconfiguration computation. The results show that the global container shipping network does not need large-scale reconfiguration but requires small-scale changes to optimize the network feature vector centrality and make the network more balanced and stable. This study is beneficial for business managers to proactively respond to the future development of the shipping network and improve the operational efficiency of the global container shipping network.


Introduction
Because of the advantages of high throughput, low costs and fewer access restrictions, maritime transportation accounts for more than 70% of global trade [1].Container shipping has become an important mode of transportation mode with many advantages such as high security, high loading and unloading and high energy efficiency.Liner shipping operation is generally employed, that is, sailing according to a predetermined schedule and fixed routes, stopping at fixed ports, and carrying out international shipping at fixed freight rates.Affected by economic, political, policy and other factors, the container shipping network is constantly developing and changing.
Regional Comprehensive Economic Partnership (RCEP) is an agreement initiated by ASEAN (Association of Southeast Asian Nations) in 2012 and developed by 15 members, including China, Japan, Korea, Australia, New Zealand and ten ASEAN countries.The signing of RCEP marks the official departure of the most populous, largest economic and trade scale and most promising free trade area.On 1 January 2022, RCEP officially entered into force.As the largest free trade zone in effect in the world, RCEP can enhance the foreign trade competitiveness of member states, especially to promote trade among countries that have not signed bilateral free trade agreements.The signing of RCEP has also had a significant impact on the maritime industry.It will significantly promote the increase in container shipping between Chinese coastal ports and Japan, South Korea, ASEAN, Australia and New Zealand, reduce export costs, and further increase the proportion of near-ocean container routes [2].Facing more and more frequent trade exchanges, reconfiguring and optimizing the global container shipping network is also a growing concern.
With the development of complex network theory, a wide range of transportation networks have been studied within the framework of complex networks, which are able to capture the interactions between the entities [3].They demonstrate structure properties that are not evident when each component is considered independently [4].For this reason, many scholars have studied different transportation modes from the perspective of complex networks in recent years.In aviation networks, Jia T. and Jiang B. [5] explored the U.S. airport model using complex network theory and classified airports according to various network metrics.Wang et al. [6] used a complex network approach to study the network structure and node centrality of cities in the Air Transportation Network of China (ATNC),and confirmed that centrality captures an important aspect of locational advantage in the ATNC and is important in shaping the spatial pattern of economic activity.LIN et al. [7] applied characteristics of nodes and edges to analyze and compare the relationship of urban traffic growth mechanisms, but there are no significant small-world characteristics in urban networks [8], which differs from shipping networks.Container shipping networks are large-scale networks consisting of many ports, routes, ships and other elements.They can also be understood as a complex structure, similarly to other transportation networks, consisting of a set of vertices (ports) connected by edges (routes) [9].
Seaports and shipping routes form the topology of the container shipping network [10].The importance of topology in transportation modes has prompted many researchers to investigate and describe the structure and dynamics that make transportation a complex network [11].Current research on maritime transportation networks has focused on the topology of the world's maritime transportation networks, including degree distributions, clustering coefficients, and average shortest path lengths [12].Meanwhile, the topological characteristics (such as number of lines, throughput, average degree, clustering coefficient, average distance, power law distribution index, and community) of the world ocean transportation network are compared with those of the sub-networks of container transportation, dry cargo transportation and oil transportation, respectively.There are also many studies on the weight distribution and seaport distance distribution of the world ocean transport network.
Connectivity is a basic concept in topology.In recent years, the interconnection levels of different infrastructure networks have been analyzed in various transportation networks such as railways, highways and shipping [13].Container shipping is a typical liner shipping method, and the United Nations Conference on Trade and Development (UNCTAD) proposes the Liner Shipping Connectivity Index (LSCI) as the basis for evaluating the trade cost and competitiveness of various countries, which covers the following aspects: the number of ships expected to berth in a country per week, the designed throughput, the number of scheduled liner routes to and from the country, the number of liner companies providing services to and from the country, the average size of ships deployed on scheduled routes and the number of countries that trade directly with the country by sea via transshipment [14].Therefore, the connectivity of the container shipping network can not only measure the liner connection index but also reflect the service capacity of the ports [15].
Since the introduction of the LSCI index, there has been an increasing amount of research on the connectivity of shipping networks.Among them, Pan et al. [16] showed that the global shipping network is represented by an unsigned Laplacian matrix, which could be decomposed to produce its eigenvectors and corresponding eigenvalues.The maximum gap among eigenvalues was then used to determine the optimal number of communities within the network.This efficient method identified major port communities and analyzed the network connectivity of the global shipping network according to the community structure.
However, understanding these structures is inherently more difficult due to their dynamic complexity [17], the evolution of the network [18], the diversity of connections [19] or the diversity of nodes.There is little research to enhance the properties of the ports and their links [20].With the rapid development of the global economy, commodity characteristics and environmental factors will also affect shipping decisions [21].We need to consider the status quo of the global economy to pay attention to the sustainable development of famous international organizations [22], carry out research on the connectivity and reconfiguration of container shipping networks, predict the next developments in maritime transport and make recommendations to the global maritime industry, and propose predictive implementation solutions for port decision-making.
In this paper, the impact of RCEP is focused on, and a new global container shipping network is established based on the complex network theory, which is more realistic and current.A new concept of future port expectations is proposed based on the existing expectations of ports, and the increase or decrease in port link expectations is calculated for the reconfiguration and optimization of global container ports.The applied complex network analysis provides insight into the operational and geographical dynamics of the ports participating in the container shipping network, helping lines and ports understand and measure their competitive position in the network [23].

Reconfiguration Model of Container Shipping Network
In the container shipping network, nodes correspond to ports with different specifics.Some ports have the advantages of being able to accommodate large ships, better geographical location, and higher route accessibility.We classify such ports as ports with high service capacity.Some ports are more remote with few direct connections, which are classified as ports with lower service capacity.However, after data analysis and literature research, it is found that the service capacity of many ports is significantly higher than the current demand.At the same time, due to the surge in trade demand, the service capacity of some ports cannot meet the transportation demand, and congestion often occurs.Therefore, the reconfiguration and optimization of the port service capability in the global container shipping network is a problem that needs to be considered.
Based on throughput, port location, distance among ports and route importance, this paper focuses on the port connection and port centrality to establish a large-scale sparse matrix linear system to realize the reconfiguration optimization of global container shipping network service capacity.A floating-point arithmetic model is established, and the optimization result of the model is measured by the largest eigenvalue of the weighted adjacency matrix.However, it would take a long time to solve the large sparse matrix, which is the bottleneck in the calculation of the port service capacity of the entire shipping network.A communication-reducing conjugate gradient method based on Krylov Subspace (CR-CG-KS) is put forward to quickly calculate the reconfiguration of shipping network service capabilities.At the same time, the service capacity reconfiguration coefficient of the global container shipping network is designed to facilitate the decision and selection for the service capacity reconfiguration optimization of global container ports, and the eigenvector centrality is used to measure the performance of reconfiguration optimization.

Container Shipping Network Model Construction
The container shipping network is a relatively complex stochastic system, and there are many factors affecting the network operation efficiency.These influencing factors will interact with each other, thus having an extremely far-reaching impact on the shipping network.The following assumptions are made for the model: (1) The service capacity and efficiency of different ports are the same, and the price competition between ports is not considered.
(2) Only the conventional container transportation network is considered, and the interference of factors such as random events is not considered.(3) The container vessels deployed on the route are all regular operating vessels of the world fleet at this stage, and the number of vessels is sufficient.(4) The distance calculation only considers the straight-line distance between two ports, without considering some passing terminals.(5) The fuel cost does not consider gaps among them and is calculated only on the global average fuel price.
Consider a container transportation network, in which  is the set of ports, E is the set of links, and W is a vector of link weights representing the link frequency. is the weighted adjacency matrix of the network, and each element of  can be expressed as: where  =  , and both  and  are the sailing frequency between ports  and .In this paper, if there is a container route from port  to port , the link is a directional oneway link.That means the container shipping network cannot be seen as a symmetrical net.
Let  be the degree matrix： is the Laplacian for the container shipping network : In this model, the largest eigenvalue of the frequency-weighted adjacency matrix is the node-weighted average link frequency.Thus, the largest eigenvalue can be used as a measure of network connectivity with the same units as the elements of the adjacency matrix.The eigenvector, which corresponds to the largest eigenvalue of a non-negative adjacency matrix with a unit criterion, is denoted as the "principal eigenvector" in this paper.In addition, this study does not require the adjacency matrix to be symmetric, as the world container transport network is highly asymmetric.
It is well known that connecting two nodes bidirectionally can ensure strong connectivity of the network.In this paper, we aim to improve the network connectivity by adding links frequency among hub ports or among ports located in RCEP or among hub ports and ports located in the RCEP.The key notations are shown in Table 1.In addition to the impact of port throughput, we also consider both the hub ports and the ports covered by the RCEP:

.2. Cost-Based Importance Evaluation of Container Ports
The fuel cost is the major component of the ship operating cost.For the transportation cost, we mainly consider the fuel cost, which is proportional to the oceanic distance.Distances between nodes are calculated in Equation ( 5) by latitude and longitude coordinates: ( ,  ) and ( ,  ), where r = 6371.004km.

Importance Evaluation of Container Shipping Routes
Suppose that the container shipping network consists of x routes and y ports on the target transport route.Then, the importance of each link on that route can be calculated as the number of routes including this link: The Borda Count is a single-winner election method.In the Borda Count method, voters rank candidates in descending order of preference, giving each candidate a certain score based on the total number of candidates involved, and ultimately each candidate is selected as the winner based on the total score received from all voters.Because the method tends to select the option that is widely accepted by voters rather than the one preferred by the majority, Borda Count is often described as a consensus-based voting system rather than a majority-minority voting system.The method provides a reasonable way to integrate the evaluation results of the various weighting indicators by avoiding result bias caused by decision makers' personal preferences.In the importance ranking of container ports, the various indicators can be treated as voters and the ports as candidates.Then the ranking of the ports obtained according to the different centrality indicators can be calculated for each port based on the following equation.
where n denotes the total number of ports in the sample network. (),  () and  () represent the scores of ports based on corresponding ranking, respectively.Based on these scores, we can derive the combined centrality score of port i as follows: where  () represents the composite score of port i.  (),  ()   () denote the ranking of port i in the sample network based on throughput capacity, transportation cost and route importance, respectively.The weights are determined by the throughput capacity, transportation cost and route importance, respectively.

Mathematical Model
From Section 2.2, A is an n × n Hermitian matrix.The Rayleigh quotient RA is: Let Uk denote a subspace of dimension n − k + 1.By checking the dimensionality, it can be known that the intersection of span{U1, ..., Uk} and the subspace Un-k+1 is not zero.Therefore, there exists a vector v at this intersection: Its Rayleigh quotient is Hence, Since this is true for all subspace U, we can conclude that Therefore, the largest eigenvalue can be regarded as the node-weighted average link frequency.Eigenvector centrality, i.e., the largest eigenvalue and its associated eigenvectors, is our chosen criterion to measure the connectivity of the container shipping network.
In the whole process of service capability reconfiguration analysis, A is the weighted adjacency matrix of the global container shipping network, and the weight is set to the frequency of anchoring and service capability.For the operation of the matrix, we use the standard floating-point operation method.
where ε is the unit rounding of the machine.Most liner optimization models of container shipping are mixed-integer models [24,25].By analyzing the frequency-weighted adjacency matrix, the largest eigenvalues and their corresponding eigenvectors, this paper proposes a max-min integer optimization model based on a Rayleigh quotient for better connectivity in the shipping network.
(, ) = 0, ∀(, ) ∉  (25) −  = 0, ∀(, ) ∈ ∈ {0,1} × (28) The objective is to maximize the largest eigenvalue by changing the service capability of ports.Constraint (17) ensures that the link frequency cannot be greater than the maximum link frequency for each node.Constraints (18) and (19) ensure that Zij = 0 if Iij = 1 for each candidate pair of nodes i and j, thus A(i,j)≥1 as Iij is a binary variable according to constraint (27).Constraints ( 20)-( 23) ensure that containers can only be routed via an open hub port.Constraints ( 24)-( 26) ensure that Bij, Iij, Zij are set as zero for each non-candidate node pair.Constraints (27) ensure that the in-coming and out-going shipping activities between ports i and j are mutual.Constraints ( 28) and ( 29) define the domain of the decision variables.This max-min model is NP-hard as it contains a fractional objective function and integer constraints.In this paper, a communication-reducing conjugate gradient method based on Krylov Subspace is proposed to find an optimal or near-optimal solution.

The Reconfiguration Optimization Algorithm
From Section 2, the global container shipping network can be seen as a large sparse matrix.The Krylov subspace methods (KSMs) are a class of iterative methods that can maintain the sparsity of the coefficient matrix and require less computation at each iteration step.Therefore, the Krylov subspace method is the preferred method for solving this type of problem.However, in the Krylov subspace iteration, the classical formulation needs to move the data in each iteration, which causes a performance bottleneck and thus increases the running time.
In this section, after the introduction of KSMs, the Krylov subspace parallel computing strategy is detailed, which can improve the calculation efficiency and maintain numerical stability.Then a modified algorithm, CR-CG-KS, is proposed to improve the stability and accuracy of reconfiguration optimization.

Krylov Subspace Methods
In the n-th iteration, the solution xn and the residual rn are updated to where xn and rn have different rounding patterns at finite precision.That is, the expression of xn does not depend on rn, while the expression of rn does not depend on xn.Therefore, computational errors that occur in xn cannot be self-corrected.These errors accumulate over many iterations, leading to a deviation between the true residual and the updated residual.Writing the true residual as  −  =  + ( −  −  ), we can constrain its criterion by || −  −  || is the magnitude of the deviation between the true residual and the updated residual.If this deviation becomes large, it will limit the maximum achievable accuracy.
In the S-step Krylov subspace algorithm, the iterative loop is divided into an outer loop and an inner loop, and the inner loop computes S steps of the iterative process for each outer loop.Such a formulation has been re-derived many times in the last decades, mostly with the aim of increasing parallelism [26] and avoiding data movement between two memory levels on sequential machines and between processors in parallel environments [27].

Krylov Subspace Parallel Computing Strategy
Since Krylov subspace computation requires frequent access to main memory, for large sparse matrix operations, data transmission in each step of computation will occupy most of the time, which is a bottleneck for efficient computation.For this problem, the operations in iterations need to be rearranged to make full use of the locally saved data.In the calculation, some data are used by each computer node, such as parallel calculation of sparse matrix-vector multiplication and inner product calculation requiring global communication; however, some data are stored in the local memory of the computer and do not need to be communicated.If they run separately, the amount of calculation data will be very large, and the simultaneous calculation will block most of the communication time.Many scholars have adopted the S-step approach to increase parallelism [26,28] and avoid data moving between two memory levels on sequential machines and between processors in parallel environments [29], which has been derived many times and proven to be important steps to enhance the parallel nature of the Krylov subspace method.
Based on the S-step method, this paper reorders the calculation sequence of the iterative method, which delays the solution correction by one iteration step.That is, the solution correction does not have to wait for the completion of the inner product calculation, thereby realizing the overlapping of communication and calculation.This method not only improves the calculation efficiency but also maintains the numerical stability of the iterative method [30].

Conjugate Gradient Algorithm (CG) for Communication Reduction
CR-CG-KS consists of an outer loop indexed by k and an inner loop iterating from j = 1 to s.Its objective is to determine the information needed to compute the conjugate gradient vectors psk+j, rsk+j and xsk+j for 1 ≤ j ≤ s and s >0.From the properties of CG, it can be known that  ,  ∈  (,  ) +  (,  )  −  ∈  (,  ) +  (,  ) The j-th Krylov subspace with starting vector v generated from the initial matrix A is The conjugate gradient vector from iteration sk + 1 to sk + s for the subspace with the starting vector Rsk can be calculated by Equations ( 34) and ( 35):  = { () , ⋯ ,  () } ( ) =  (,  ) (34)  = { () , ⋯ ,  () } ( ) =  (,  ) (35) where  () is a polynomial of degree j that satisfies the three-term recurrence formula: In the inner loop, we just update its coordinates in ,: It is set that  =  ⋅  in the inner loop computation, which greatly reduces the memory cost in each iteration.The matrix  can be computed by the update reduction in each outer loop, and its asymptotic delay cost is the same as that of a single inner product computation.Since the dimensions of  and  are both (2s + 1) × (2s + 1),  and  can be computed locally in the inner loop. where Assuming the algorithm terminates at step sk+j, return  and ̂ .The deviation between the true residual  −  and the updated residual  −  are limited as follows: Equation ( 46) can be further simplified as: It is expected that the true residuals and updated residuals of iteration sk+j in CR-CG-KS are the same size as those in CG, which play an important role in determining the accuracy of CR-CG-KS relative to CG. to ensure that there is a non-linear increase in error since the last substitution step to avoid unnecessary substitutions.Then, the algorithm must start a new outer loop that generates a new S-step base relative to the residuals being replaced.
Next, to perform an inertial update, we must be able to estimate the norm of A easily, preferably only once.We discuss how to obtain the remaining norms without increasing the asymptotic computation or communication cost of CR-CG-KS.
Hence, we can finally obtain

Numerical Processing and Preliminary Network Generation
According to Alphaliner's statistics in 2021, the three major shipping alliances account for about 79% of the global container shipping market share (2M accounts for 34%, Ocean Alliance accounts for 28%, THE Alliance accounts for 17%), and the remaining shipping companies only account for 21% of shipping market share.The top 10 global container shipping companies and their market shares in 2021 are shown in Table 2. To construct a global container shipping network under RCEP, we use the information published on the websites of Clarkson, Maersk, CMA-CGM, Evergreen Marine, Hapag-Lloyd, COSCO Container, and CSCL in 2021.The schedules show the paths of ships calling at ports during their voyages, including the names of ports and shipping routes.In our network model, links are established only between consecutive ports along the shipping routes [31,32].In total, 237 major ports in the world and 1434 undirected links among them are included in the sample for the global container shipping network.The weight of each link is represented by the service frequency.Thus, the frequency-weighted adjacency matrix A provides a realistic representation of the global shipping network.The global container shipping network is constructed based on the construction principle of the L space network [33].

Numerical Example
A 256-order sparse matrix is used as the test matrix to evaluate the residual replacement strategy of this paper.The error is set to 10 , and  ' =  =  =  = 1 [30].CG and the S-step Krylov subspace conjugate gradient algorithm (S-CG) are used to conduct comparative experiments.The result comparisons when S = 4, 8 and 12 are shown in Figure 1a-c From Figure 1, CG reaches convergence at the 33rd time when S = 4, and the convergence speed is the fastest, but the residual values are larger than the other two methods, and the calculation accuracy is not high.When S = 8, the number of iterations is 41, and the residual value is also more than 100 times different from those when S = 4.When S = 12, the number of iterations is 217.Although the residual meets the convergence requirements, the calculation time is long, and the calculation memory is too much.It can be seen that the number of inner loops has a great influence on the calculation results obtained by CG.
In contrast, S-CG sequentially selects the orthogonal S-direction vectors of A to form an orthogonal orientation matrix, while minimizing the residual vectors in these S-directions) [28].When S = 4 and S = 8, the variation of iteration number and residual distance is similar to those of CG.The residual value reaches the optimal value when S = 8, but the convergence requirement cannot be satisfied when S = 12.
CR-CG-KS proposed in this paper has the same iteration number and residual values when S = 4, 8 and 12.The residual values are the minimum values among those of these three algorithms.It can be seen that CR-CG-KS has better stability and accuracy.

The Major Container Shipping Network
According to the data collection of the global container shipping network with 237 major ports and considering the practical application, the data are discrete and analyzed to obtain the optimal solution and scheme.Figure 2a-c show the result comparisons among the three methods.From Figure 2a-c, the convergence speed of CG is always the fastest.However, its final convergence accuracy is worse than that of CR-CG-KS.During the early iteration period, the errors of S-CG and CR-CG-KS are almost the same regardless of whether the inner loop step is 4, 8 or 12.However, the convergence speed of S-CG becomes much lower than that of CR-CG-KS, and its final error is significantly larger.As S increases, the convergence speed of CG and S-CG becomes slower, and their convergence accuracy becomes worse.The change of S has little effect on the convergence speed and accuracy of CR-CG-KS.The errors of CR-CG-KS always reach their minimum in the 126th iteration, and the minimum errors of CR-CG-KS are always 1.26 × 10 −14 .It can be seen that the improvement of CR-CG-KS for the global container shipping network reconfiguration is significant.
Specifically, we compare the expected reference value of the global container shipping network with the planned reconfiguration coefficients.The global container shipping network reconfiguration results are shown in Figure 3a, and their fitting comparison is shown in Figure 3b.If P is greater than or equal to 10%, it is needed to enhance the corresponding service capability; if P is less than or equal to −10%, the service capability is redundant, which can be reduced.
where  () is the preliminary calculation value of the global container port service capacity,  () is the global container port service capacity calculated by CR-CG-KS.From Figure 3a, there is a certain difference between the calculation result and the predetermined reconfiguration result, which indicates that some container ports need reconfiguration optimization to a certain extent.Figure 3b shows the expecting reconfiguration reference value and the first-order linear fitting of the calculated reconfiguration results.The slope and the intercept of the fitted curve are 0.920 and 0.005, respectively.It can be seen that the change in the reconfiguration value is not significant, which is highly similar to the preset reconfiguration value.
In particular, ports with sufficient service capacity and underserved capacity are listed in Tables 3 and 4, respectively.Proportions of global container transport network reconfiguration with different P ranges are shown in Table 5.From Tables 3-5, there are three ports with P below -10%, accounting for 1% of all ports, namely Tampa, Rijeka, Ensenada, which have redundant service capabilities and need to be reduced; eight ports with more than 10%, accounting for 3% of all ports, namely Taipei, Miami, Denmark, Dublin, Kaliningrad, Morocco, Bangkok, Yantai, which need to increase their service capabilities.In addition, 122 ports with −10% to 0, accounting for 52% of all ports; 108 ports with 0-10%, accounting for 44% of all ports.These 230 ports do not need to optimize the port service capacity.

The Global Container Shipping Network Considering Non-Major Ports
Based on these 237 more important ports, we continued to collect some other relatively less important ports.Finally, we collected a total of 625 ports.Similar to the reconfiguration and optimization of the major container shipping networks in the world, we perform a similar simulation analysis on the global container shipping networks of these 625 ports.The result comparisons among three methods with different S are shown in Figure 4a-c.
The expected reference value of the global container shipping network was calculated and compared with the calculated planned reconfiguration coefficients.The global container shipping network reconfiguration results are shown in Figure 5a, and their fitting comparison is shown in Figure 5b.From Figure 5a, the calculation results of the reconfiguration and optimization of the global container shipping network of 625 ports are somewhat different from the reconfiguration results present in this paper.It indicates that the global container shipping network has a certain difference, and some ports need to be refactored and optimized.Figure 5b gives the first-order linear fit of the expected reconfiguration reference value and the calculated reconfiguration result.From Figure 5b, compared with the major container shipping network in Section 4.3.1,there are a few ports with a larger difference; but at the same time, the slope and intercept of the fitted curve are 0.9539 and 0.0015, respectively, which are larger than the first order of the previous one.The fitting effect is better.It can be seen that the method in this paper is also applicable to larger calculation examples.Similarly, the reconfiguration proportions of the global container shipping network under different P ranges are shown in Table 6.It can be seen that 25 ports have redundant service capabilities and need to be reduced, and 23 ports need to increase their service capabilities.Their proportion in all ports are 4.00% and 3.68%, respectively.Clarksons noted that key "hot spots" for container network congestion this year include the United States, China and Northern Europe.The negative impact of emergencies such as the new crown epidemic on the maritime industry will be very long-term and profound.Even if the epidemic is over and the port congestion problem is solved, the maritime market will need sufficient time to restore its former order and vitality.
Based on the eigenvector centrality, this paper scores the centrality of the global container shipping network and ranks the top 10 centrality score tables in Table 7.It can be From the analysis in Section 4.2, it can be seen that CR-CG-KS is not very sensitive to S. Therefore, in this section, we will only use the inner loop parameter of S = 4 for the sensitivity simulation of the shipping network.Since Shanghai Port is a core port, the throughput of its directly connected ports decreases, which has a great impact on the entire network.In the container shipping network of 237 ports, the optimal solution is obtained in the 126th iteration, and the residual 2-norm is 1.26 × 10 −14 ; in the container shipping network of 625 ports, the optimal solution is obtained in the 587th iteration, the residual 2-norm is 1.345 ×10 −14 .Different from the normal network in Sections 4.3.1 and 4.3.2, the container shipping network with 237 ports needs to iterate 125 times, and the residual 2-norm is 6.158 ×10 −15 ; the container shipping network with 625 ports needs to iterate 583 times, and the residual 2-norm is 7.454 ×10 −15 .It can be seen that the number of iterations of the algorithm in this paper increases a lot with the scale of the network.Although the throughput of some ports has changed greatly, even if this port is a very important port, it has little effect on the number of iterations of the algorithm.However, when the throughput of the core port is reduced, the residual 2-norm is also reduced.

Network Connectivity Analysis
The reconfiguration will change the centralities of ports, the links and the topology of container shipping network.For example, there are obvious changes in the central ranking of the 11 ports that need to be optimized, which is shown in Table 8; and the major container shipping network with 237 ports after reconfiguration has a significant difference to that before reconfiguration, which can be seen in Figure 7a,b.This will lead to changes in network performance.Network efficiency is an important property of network performance research.In the domain of complex networks, Latora and Marchiori demonstrated the relationship between network efficiency and the shortest paths [34].In the research of container transportation networks, network efficiency such as E[A] in Equation (52) has been widely used as a measure of connectivity.
As seen in Equation ( 52), the higher the link weight, the lower the efficiency, which is not applicable to the case of service capacity considered in this paper.Wang [35] used eigenvector centrality to determine the accessibility of ports in liner shipping networks.Their results validated the superiority of this approach and demonstrated that the principal eigenvectors were also an important indicator of the nature of ports, and eigenvector centrality could provide a good measure of network connectivity.
To verify the connectivity after reconfiguration of the global container shipping network, this paper takes the product of the network efficiency and the largest eigenvalue as the connectivity measure coefficient and calculates the connectivity improvement proportional coefficient based on the model in Section 2.3: where ( ) and ( ) represent the network efficiency before and after network reconfiguration, respectively, and  and  represent the maximum eigenvalues before and after the reconfiguration of the global container shipping network, respectively.From the calculation, the connectivity coefficient of the major container shipping network with 237 ports before reconfiguration is 50.72, while the eigenvalue after reconfiguration is 52.88, and the overall property is improved by 4.26%.The connectivity coefficient of the global container shipping network with 625 ports is 183.29, while the eigenvalue after optimization is 190.53, and the overall property is improved by 3.95%.It can be seen that after the service capability reconfiguration of the global container shipping network, its connectivity has also been improved.

Conclusions
In this paper, a new mixed integer model of global containerized shipping network reconfiguration is developed based on the combination of connectivity frequency, throughput, port distance and route importance.A communication-reducing conjugate gradient method based on Krylov Subspace is proposed to quickly calculate the shipping network reconfiguration values.The comprehensive reconfiguration performance is measured by the maximum eigenvalue of the weighted adjacency matrix.Due to the rapid development of the maritime industry and the increasing emphasis on containerized shipping, this study is beneficial for business managers to proactively respond to the future development of the maritime network and to improve the operational efficiency of the global containerized shipping network.
In future work, we will comprehensively consider the reconfigurability of ports, routes, and container flows to further improve the practicability of models and algorithms.In addition, we will investigate the possibility of applying the proposed approach to other transport networks.

Figure 7 .
Figure 7.The global container shipping network comparison under service capability configuration: (a) before configuration (b) after configuration.
NR Number of ports in RCEP  A binary variable that equals 1 if port k is port in RCEP, and 0 otherwise  A binary variable that equals 1 if port k is port in set H, and 0 otherwise qij The number of containers to be transported from origin port i to destination port j Oij The cost for transporting one TEU (twenty-foot equivalent unit) container between ports i and j.DISij The distance between ports (or dummy nodes) i and j M Set of links in the container shipping network C Set of container shipping routes  If there is a link between ports i and j, and the link belongs to both sets M and C, it is 1; otherwise, it 2.2.Container Port Service Capability 2.2.1.Importance Evaluation of Container Ports Based on Location The optimization model maximizing the algebraic connectivity is as follows:

Table 2 .
Top 10 container shipping companies in the world as of 2021.

Table 3 .
Excellent service capacity ports.

Table 4 .
Service capacity needs to be optimized for these ports.

Table 5 .
Proportion of global container transport network reconfiguration.

Table 6 .
Proportion of global container transport network reconfiguration.