Discrete Dynamics in Nature and Society

Volume 2013 (2013), Article ID 391709, 8 pages

http://dx.doi.org/10.1155/2013/391709

## Robustness Measure of China’s Railway Network Topology Using Relative Entropy

School of Traffic and Transportation Engineering, Central South University, Changsha 410075, China

Received 4 July 2013; Accepted 5 September 2013

Academic Editor: Wuhong Wang

Copyright © 2013 Fenling Feng and Lei Wang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

This study focused on China’s railway network topology issue from a robustness measure perspective. Relative entropy is used in the study as a measurement of robustness of railway network topology. It is found that the entropy-based measure provides more informative analyses compared to a traditional graph measure. The results indicate that the railway network in the 12th five-year plan has improved robustness when compare to 2008 with respect to deliberate and random attacks.

#### 1. Introduction

Since China’s first railway was built in 1865, there have been tremendous achievements in its development over 145 years. Many new stations and lines were constructed; a huge network has been formed. The planning of a railway network should not only be evaluated from external indexes, but also analyzed on its overall performance. The research on railway network topology structure takes the railway network as a complicated system and studies its topology structure using statistical methods, to provide a theoretical basis and foundation for network planning, construction, and scientific research. The robustness of network topology is one of the most important and basic features of complicated systems.

Research on robustness is based primarily on graph theory, which includes average shortest path, node degree, aggregation coefficient, and the number of interfaces. But these measures of robustness do not offer adequate amounts of information. This article uses relative entropy as the robustness measure of China’s railway network topology. It also analyzes the relationship of relative entropy and ordering employing information theory methods from a systematic view. The robustness of China’s railway network in 2008 and in the 12th five-year plan was analyzed.

At present, in robustness research of network topology, scholars usually take definitions from graph theory as measures of robustness. Albert et al. [1] compared the connectivity of random networks and scale-free networks (a typical complex network) and further pointed out that the robustness of the two networks is significantly different under random attacks and malicious attacks. Watts [2] and Strogatz [3] put forward the concept of small world networks. Cohen et al. [4] found the robustness of an Internet network under random attacks and malicious attacks is based on Percolation theory. Bollobas and Riordan [5] made a mathematical analysis for robustness complex networks based on random graph theory. And Bars et al. [6] proposed a robustness measure in accordance with the Gastinel-Kahan theorem investigating stability conditions of discrete perturbed closed-loop systems. Rossi [7] presented new results for assessing the robustness of a configuration for multipurpose machines. The robustness measure of a configuration is returned by assessing the minimum magnitude of disturbances affecting the forecast demand that may lead to breaking the deadline provided by the decision maker. Border et al. [8] found that the macroscopic structure of the Web is considerably more intricate than suggested by earlier experiments on a smaller scale.

The following research used graph theory and measured the invulnerability of the network topological structure from the angle of connectivity. Holme et al. [9] used the overall efficiency and max connected subgraph to measure the network performance after an attack. Moreno et al. [10] found the test of strength for complex networks such as the Internet was more stringent than others recent tests, like the random removal of nodes. Xiao and Dong [11] defined the connectivity of , and it represents average communication paths between all nodes of the group, which is to measure and evaluate the survivability of network topology. Lu and Dong’s [12] domain is the topology of the network invulnerability as the total number of paths which may establish within a network group divided by the pathways between the groups. Criado et al. [13] considered it necessary for the network survivability index to have some important characteristics of the definition of network survivability. Then, they gave the invulnerability index a function definition and two kinds of survivability functions.

Railway networks have been hot topic system studies and measure methods. Banik and Dasgupta [14] deal with the use of Petri nets in modeling railway networks and designing appropriate control logic to avoid collisions. Feng et al. [15] used a spline interpolation method, a numerical differential five-point formula and the method of least squares to solve for a synergistic coefficient. Tomoeda et al. [16] presented a rescheduling method of homogenization to alleviate congestion of crowded train. Khemakhem et al. [17] provided an accurate estimate of schedule robustness, and they introduced surrogate measures. Santiago et al. [18] found through empirically analyzed that the SDH network operated by Telefónica in Spain shares remarkable topological properties with other real complex networks, such as the Internet network. Wang et al. [19, 20] analyzed the methods for the safety of traffic systems.

For years, scholars focused on using entropy to measure complex systems. Xu and Hu [21] proposed a new degree dependence entropy (DDE) descriptor to describe the degree dependence relationship and corresponding characteristics. Anand and Bianconi [22] defined the Shannon entropy of a network ensemble and proposed that it relates to the Gibbs and von Neumann entropies of network ensembles.

The paper will concentrate on the domain of proposing a new robustness approach for measuring China’s railway network. A model based on a simulated Chinas railway network will be described. Results from an implementation with a random attack network will be reported.

#### 2. Model of Railway Geography Network

From a geographical point of view, the main elements which constitute a railway network are the railway station and the railway line. Therefore, the actual railway network can be abstracted as a geography network by viewing the railway station as nodes and the railway line between each node as a side. With regard to the construction of the railway geography network, the following assumptions are needed:(1)the railway line is bidirectional; the railway geography network is undirected as directions of its sides are not taken into account;(2)the importance of the railway station and the railway line is relative; that is to say, the railway geography network is an unweighted network;(3)the carrying capacity of the railway line will not be taken into consideration;(4)multiple sites in the same city are regarded as a unified one (e.g., Beijing railway station, Beijing west railway station, Beijing south railway station, and Beijing north railway station in Beijing, are all seen as “Beijing”).

According to this assumption, we can build a railway geography network. The network is an undirected and unweighted diagram connected by the railway station and the railway line, which can be described as , of which represents a collection of various railway stations and represents a collection of railway lines connecting each railway station.

The following takes Urumqi station and its nearby stations as an example to illustrate this railway geographic network model. Figure 1 is the actual railway network. The railway geographic network can be worked out after it is abstracted, and a sketch can be derived. That is Figure 2, where . From Figure 2, we can further obtain its adjacency matrix, shown in Figure 3.

As there are five thousand stations of different levels or so across the country, it is hard to obtain so much data. This paper researches some of the main sites according to China’s railway line diagrams published by the Ministry of Railways in December 2008 and in the 12th five-year plan released in October 2012. Some main sites are to be taken to establish the 2008 railway geography network and the 12th five-year plan railway geography network.

#### 3. Measure Model of Railway Geography Network

##### 3.1. Definition of Relative Entropy

Entropy is a measure of the uncertainty of random variables, and it is also a measure of the amount of information required to describe random variables in an ordinary sense.

Relative entropy is a measure of the distance between two random distributions. In statistics, it corresponds to the logarithm expectation of likelihood ratios. The relative entropy or Kullback-Leibler distance of two probability density functions and is defined as follows:

In accordance with the intrinsic meaning of relative entropy, relative entropy can be viewed as the ordering distance between the railway network and a standard network. Therefore, as long as the completely disordered network is determined, the ordering distance between the railway network in 2008 and a completely disordered network can be obtained by relative entropy. Then we can get quantitative data on the railway network robustness. Therefore, the relative entropy can be seen as a measure of the ordering of the railway network in 2008. The robustness of the railway network can be obtained by analyzing the variation of relative entropy when the railway network is attacked.

##### 3.2. Properties of Relative Entropy

In line with convention, the relative entropy possesses properties as follows:

Therefore, if there is a character to satisfy , , then .

In addition, relative entropy has nonsymmetry; namely, . The proof is as follows.

Provided that , has two distributions, , , , , then the following expressions are established:

Suppose that , , the following expressions can be obtained by calculating:

So it is proved.

The nonsymmetry of relative entropy refers to the distance from to , and is different from the distance from to . So, it is a must to determine what the network and represent, respectively. Typically, represents true distribution, while represents theoretical distribution. So, will be deemed the railway network distribution being studied, and will be treated as the railway network with a completely disordered distribution.

##### 3.3. Railway Network with a Completely Disordered Distribution—Maximum Information Entropy

Next, the completely disordered distribution of the railway network will be determined. Information theory only gives the solve function for the maximum entropy. This paper will take China’s railway network as a complex system. Topology information entropy is at its maximum when the railway network is in a uniform structure at the moment the ordering of the system is worst; when the railway network is configured with a star structure, topology information entropy is at its minimum, and here the ordering of the system is best. The construction of the railway network is a process that moves from disorder to order. Therefore, the condition when ordering is at its worst is the worst case for the railway network. This paper selects maximum information entropy as the complete railway network with a disordered distribution.

The relative entropy is the ordering distance between the railway network being studied and the railway network with a completely disordered distribution. This means that the larger relative entropy is, the more ordered the railway network is.

When the railway network is completely disordered, all nodes have the same degree distribution ; that is to say, there is for any node , . At the moment the maximum information entropy is the following:

##### 3.4. Measure of Railway Geography Network Robustness

Accepting the definition of relative entropy as a measure of robustness requires a distribution as a basis. This paper is intended to employ node degree distribution.

Degree is defined as the number of adjacent sides of the nodes; namely , which represents the total number of the sides connected to points. Degree distribution function represents the probability that any node’s degree is . In another words, degree distribution means the proportion of nodes with degree of among the total number of nodes in the network. As an important geometrical property of the network, degree distribution could describe some characteristics of the network.

If the degree distribution density function of the railway network is expressed as , then the robustness measure can be obtained from the expression (1) as follows:

#### 4. Proof of Measure Testability of China Railway Geography Network

To define a new measure, it is necessary to analyze whether or not it is a measurable function. Next, analysis is to be conducted on whether the relative entropy function is a measurable function or not by using the Lebsegue measure theory in the function of read variable.

##### 4.1. Definition and Theorem

First, a few basic definitions and theorems.

*Definition 1 (measurable set). *Set ; if any collection , there is

Then is Lebsegue measurable or measurable in short, is the outer measure, is supplementary set of , and is the whole measurable set.

*Definition 2 (measurable function). *Set is a measurable set, and set is a function from to . If any section , there is

Then is Lebsegue measurable, or measurable in short.

*Definition 3 (measure). *Measure is a set of functions , defined in the collection of sets constituted by the subsets of given set , valued in the expansion of real functions and meeting certain conditions. These conditions should be consistent with the following properties.(1)Normative: .(2)Nonnegativity: for all , .(3)Countable additively: for any non overlapping row in , there is .

Theorem 4. *Set is a measurable space, and set is a function defined in the measurable space. Then the fact that is a measurable function and the fact that for any , is equal.*

##### 4.2. Construction of Measurable Sets

*Step 1. *Since the relative entropy is actually a probability distribution function of a set of random variables, a probability space is constructed first. When it comes to the Chinese railway network topology, the probability space can be built based on the node degree distribution. Assuming the selection of degree is the event , is maximum degree of nodes, sample space is as follows:

The value of degree is derived from a function on the sample space. From this, we can use the subsets of the sample space to make a -domain . Proportion of degree is the probability that a node’s degree is . For this reason, set , ; if , then . So a probability space is set for random phenomena as the degree distribution of Chinese railway network topology.

*Step 2. *According to the probability distribution and the definition of probability measure, a new probability space is induced in the above random variables, of which ( refers to Borel collection); there is . Function is a probability distribution of random variable and distribution in short. According to the real variable theorem, this set function is a probability measure on the set functions .

Here, it is required to know that and are both the set of events, and the difference between them is that refers to setting a few values for node degrees, while refers to when the value of several node degrees are real numbers. In other words, is a simple event and belongs to the set of real numbers. However, all real numbers are basic events, and all Borel sets are events; hence, is also an event.

Thus, a new probability space is induced from the probability space . represents a probability measure of event .

*Step 3. *As the relative entropy is the probability distribution function of the random variable, the relative entropy function can be written as . According to , we can figure out that is a composite function.

##### 4.3. Proof of Testability

Next, we will prove that the relative entropy is a measurable function.

Because the expression is composed of simple functions, is a continuous function in . Because , is an open set in for any real number and an open set is a Borel set, is Borel measurable. Since is measurable, according to Theorem 4, there is . So for any , , composite function , is a measurable function in , according to Theorem 4.

At this point, it is proven that the relative entropy is a measurable function.

#### 5. Analysis of Railway Network Robustness in “The 12th Five-Year Plan”

To test and verify the utility of the new robustness measure, two networks were selected. The first one is a circuit diagram of China’s railway network in 2008, in which 313 nodes and 472 sides were selected; the other is a circuit diagram of China’s railway network in “the 12th five-year plan”, in which there are 420 nodes and 753 sides selected (Figure 4). By comparing the relative entropy changes of railway network in “the 12th five-year Plan” and the railway network in 2008 under two different types of attacks, we can get to know each configuration’s robustness and their pros and cons.

The superiority of the relative entropy measure can be shown by comparing relative entropy with the clustering coefficient of graph theory. Clustering coefficient is defined as follows: if the node is collected to other nodes by sides, then there are sides at most between the sides. So, the clustering coefficient is a ratio of the actual number of sides between the sides and the maximum number of possible sides . can be described as .

##### 5.1. Initial State

We can determine that the relative entropy of China’s railway network in 2008 is 0.6501 and in “the 12th five-year plan” is 0.6931 by inputting the initial data of the railway geography network in 2008 and the “12th five-year plan” to the expression (6). The calculated average clustering coefficient in 2008 is 0.0699 and in “the 12th five-year plan” is 0.0918. From this, ordering and aggregation degrees of node in “the 12th five-year plan” are higher than in 2008. It is observed that the conclusion derived from relative entropy and average clustering coefficient is consistent.

##### 5.2. Deliberate Attack

In this part, we’ll analyze the robustness of the railway geography network in 2008 and “the 12th five-year plan” under deliberate attack. Deliberate attack refers to those attacks on important sites on purpose. It shows up in the deletion of node degrees from largest to smallest on the geography network model. The result is shown in Figure 5.

From the point of view of relative entropy, similarities of the two railway networks under deliberate attack are as follows. Along with increasing attack strength, the relative entropy of railway network declines, and the ordering gets lower and lower. The fact that relative entropy fluctuates is because that deliberate attack is from larger node degree to smaller node degree. When the nodes of large value are deleted, the railway network will be divided into several smaller networks. The relative entropy is a superposition of ordering of these small networks, so there will be a rally. When deliberate attacks increase to a certain extent, the railway will be completely disordered. And when the node degree is 1, each node will be of equal status, and then the network will remain in a completely disordered state until the attack is over.

Compared to the railway network in 2008, the railway network in “the 12th five-year plan” has the following advantages. During the early stage when a deliberate attack is about 10 percent, the relative entropy is decreasing slower. But the railway network in “the 12th five-year plan" still has higher ordering and stronger robustness when the deliberate attack is relatively small. The change in relative entropy of the railway network in “the 12th five-year plan” is more stable than that in 2008. The railway network in “the 12th five-year plan” becomes quickly disordered when the deliberate attack is at 70 percent. But the railway network in 2008 becomes quickly disordered when the deliberate attack is at 50 percent. The railway network in “the 12th five-year plan” becomes completely disordered when the deliberate attack is at 90 percent. But the railway network in 2008 becomes completely disordered when the deliberate attack is 80 percent.

From the point of view of the average degree of aggregation, the average degree of aggregation in 2008 and “the 12th five-year plan” is declining under deliberate attack. When attacks reach a certain extent, the clustering coefficient is 0. The decreasing trend of average degree of aggregation in “the 12th five-year plan” is more stable than that in 2008. Clustering coefficient is more durable than that in 2008.

##### 5.3. Random Attack

After a description of a deliberate attack, analysis will be given on the robustness of railway geography network in 2008 and “the 12th five-year plan” under random attack. Random attack refers to aimless attacks on sites. In the geography network model, a random attack deletes node degrees randomly. The result is shown in Figure 6.

From the point of view of relative entropy, both the railway networks in 2008 and “the 12th five-year plan” have common downward trends and volatility changes and end in complete disorder, just the same as with those under deliberate attack.

In the meantime, the railway network in “the 12th five-year plan” under random attacks still has the following advantages. When a random attack is less than 10 percent, the ordering of the railway network in “the 12th five-year plan” is higher than that in 2008. The railway network in “the 12th five-year plan” is more stable than that in 2008 with stronger robustness. The railway network in “the 12th five-year plan” becomes quickly disordered when random attacks reach 90 percent. But the railway network in 2008 becomes quickly disordered when random attacks are at 80 percent or so. The railway network in “the 12th five-year plan” arrives at complete disorder when the random attack is 98 percent. But the railway network in 2008 becomes completely disordered when the random attack is at 95 percent.

From the point of view of the average degree of aggregation, the average degree of aggregation in 2008 and “the 12th five-year plan” declines under random attack. The clustering coefficient will be 0 in the end, which is the same as that under deliberate attack.

The decreasing trend of average degree of aggregation in “the 12th five-year plan” is more stable than that in 2008. The clustering coefficient becomes 0 later than that in 2008.

By comparing the conclusions above, we can see that on the one hand, as a robustness measure, relative entropy could be used to get similar conclusions. On the other hand, the relative entropy measure can provide this information as follows. The overall trend of networks under attack. The rate of measure change at the initial attack stage. When networks will become in rapidly disorder at the initial attack stage. When networks will be in a state of almost complete disorder at the last stage. The common graph theory measure can only show the first and the fourth information. So, taking relative entropy into account can provide more accurate figures and greater amount of information.

Furthermore, representing the order of networks, the relative entropy can be used to analyze the robustness for a view of systems, which means that relative entropy has a measure function beyond graph theory.

#### 6. Conclusions

This paper deals with the robustness measure approach to railway network. First, the entropy-based robustness measure approach is established. It was found that the entropy-based method provides greater amounts of information than the graph theory measure. Based on the method proposed in this paper, stronger robustness was found for the 12th five-year plan railway network as compared to that of 2008 network, with respect to deliberate attack. Furthermore, the 12th five-year plan railway network has greater stability and ordering than the 2008 network.

Further study can be conducted to compare the relative entropy measure to other graph theory and existing research on the robustness of railway traffic networks and transport networks in China.

Network robustness is not unique to railway networks, the relative entropy measure can also be used with other types of complex networks, such as biological networks, internet networks, and social networks. This method to assess the relative entropy is also universal.

#### Acknowledgments

The authors would like to thank Matt Rolph and center for communication practices at RPI for the discussions and comments on the paper. Project supported by the Fundamental Research Funds for the Central Universities (Grant no. 2010QZZD021) and the Ministry of Railway Science and Technology Research Development Program (Grant no. 2012X012-A), China.

#### References

- R. Albert, H. Jeong, and A. L. Barabasi, “Error and attack tolerance of complex networks,”
*Nature*, vol. 406, no. 6794, pp. 378–382, 2000. View at Google Scholar - D. J. Watts, “The “new” science of networks,”
*Annual Review of Sociology*, vol. 30, pp. 243–270, 2004. View at Publisher · View at Google Scholar · View at Scopus - S. H. Strogatz, “Exploring complex networks,”
*Nature*, vol. 410, no. 6825, pp. 268–276, 2001. View at Publisher · View at Google Scholar · View at Scopus - R. Cohen, K. Erez, D. Ben-Avraham, and S. Havlin, “Resilience of the Internet to random breakdowns,”
*Physical Review Letters*, vol. 85, no. 21, pp. 4626–4628, 2000. View at Publisher · View at Google Scholar · View at Scopus - B. Bollobas and O. Riordan, “Robustness and vulnerability of scale-free random graphs,”
*Internet Mathematics*, vol. 1, no. 1, pp. 1–35, 2003. View at Google Scholar - R. Bars, I. Jaworska, and S. Tzafestas, “A mimo discrete-system robustness measure for internal model control structure,”
*Systems Science*, vol. 23, no. 1, pp. 25–35, 1997. View at Google Scholar · View at Scopus - A. Rossi, “A robustness measure of the configuration of multi-purpose machines,”
*International Journal of Production Research*, vol. 48, no. 4, pp. 1013–1033, 2010. View at Publisher · View at Google Scholar · View at Scopus - A. Broder, R. Kumar, F. Maghoul et al., “Graph structure in the Web,”
*Computer Networks*, vol. 33, no. 1–6, pp. 309–320, 2000. View at Publisher · View at Google Scholar · View at Scopus - P. Holme, B. J. Kim, C. N. Yoon, and S. K. Han, “Attack vulnerability of complex networks,”
*Physical Review E*, vol. 65, no. 5, Article ID 056109, 14 pages, 2002. View at Publisher · View at Google Scholar · View at Scopus - Y. Moreno, J. B. Gómez, and A. F. Pacheco, “Instability of scale-free networks under node-breaking avalanches,”
*Europhysics Letters*, vol. 58, no. 4, pp. 630–636, 2002. View at Publisher · View at Google Scholar · View at Scopus - Z. Xiao and Z. Dong, “Improved GIB synchronization method for OFDM system,”
*IEEE Telecommunications*, vol. 2, no. 8, pp. 1417–1421, 2003. View at Google Scholar - H. Lu and Z. Dong, “Carrier frequency offset estimation of DAB receiver based on phase reference symbol,”
*IEEE Transactions on Consumer Electronics*, vol. 46, no. 1, pp. 127–130, 2000. View at Publisher · View at Google Scholar · View at Scopus - R. Criado, J. Flores, B. Hernández-Bermejo, J. Pello, and M. Romance, “Effective measurement of network vulnerability under random and intentional attacks,”
*Journal of Mathematical Modelling and Algorithms*, vol. 4, no. 3, pp. 307–316, 2005. View at Publisher · View at Google Scholar · View at Scopus - M. Banik and R. Dasgupta, “Modeling of railway network using petri nets,” in
*Proceedings of the International Conference on Modeling, Simulation & Visualization Methods*, pp. 175–179, 2011. - F. Feng, D. Lan, and L. Yang, “Analysis on the synergy evolutionary development of the collecting, distributing, and transporting system of railway heavy haul ransportation,”
*Discrete Dynamics in Nature and Society*, vol. 2012, Article ID 329120, 22 pages, 2012. View at Publisher · View at Google Scholar - A. Tomoeda, M. Komatsu, I. Y. Yoo et al., “Real-time railway network simulation and alleviating congestion of crowded trains,” in
*Proceedings of the ICROS-SICE International Joint Conference (ICCAS-SICE '09)*, pp. 4932–4936, August 2009. View at Scopus - M. A. Khemakhem and H. Chtourou, “Efficient robustness measures for the resource-constrained project scheduling problem,”
*International Journal of Industrial and Systems Engineering*, vol. 14, no. 2, pp. 245–267, 2013. View at Google Scholar - A. Santiago, J. P. Cardenas, M. L. Mouronte, V. Feliu, and R. M. Benito, “Modeling the topology of SDH Networks,”
*International Journal of Modern Physics C*, vol. 19, no. 12, pp. 1809–1820, 2008. View at Publisher · View at Google Scholar · View at Scopus - W. Wang, W. Zhang, H. Guo, H. Bubb, and K. Ikeuchi, “A safety-based approaching behavioural model with various driving characteristics,”
*Transportation Research C*, vol. 19, no. 6, pp. 1202–1214, 2011. View at Publisher · View at Google Scholar · View at Scopus - W. Wang, X. Jiang, S. Xia, and Q. Cao, “Incident tree model and incident tree analysis method for quantified risk assessment: an in-depth accident study in traffic operation,”
*Safety Science*, vol. 48, no. 10, pp. 1248–1262, 2010. View at Publisher · View at Google Scholar · View at Scopus - X. Xu and X. Hu, “Degree dependence entropy: a new descriptor for complex networks,” in
*Communications in Computer and Information Science*, vol. 3 of*AsiaSim*, pp. 19–26, 2012. View at Google Scholar - K. Anand and G. Bianconi, “Entropy measures for networks: toward an information theory of complex topologies,”
*Physical Review E*, vol. 80, no. 4, Article ID 045102, 2009. View at Publisher · View at Google Scholar · View at Scopus