A Reliability-Oriented Local-Area Model for Large-Scale Wireless Sensor Networks

Large-scale wireless sensor networks (WSNs) have demonstrated some complex features which are similar to those of other types of complex networks, such as social networks. Based on these complex features, evolution process and characteristic of WSNs, we represent a WSN topologically by building a suitable model, which is named as the reliability-oriented local-area model (ROLM) and aimed at improving the performance of WSNs. For analyzing the performance of the ROLM, we define the reliability as the probability of that the relative error between the measurement and the true value is equal to or less than ε (ε ≥ 0) and proposed a parameter η to measure the reliability of the network. Based on them, we use η to analyze the influence of network structure on the reliability, and compared the reliabilities of the ROLM and the existing WSNs. Experiment results prove that the large-scale WSN follows a power-law distribution, and it has scale-free characteristic and small world characteristic. And it also shows that, comparing with existingmodel, ROLMnot only balances energy consumption by limiting the connectivity of each node to prolong the lifetime of the network, but also improves the reliability substantially. And the ROLM can be used to express the topology of reliability-oriented WSNs and analyze the structure preferably.


Introduction
Complex networks are currently being studied in many scientific fields [1] and many systems actually can be described through complex networks such as biological networks [2], metabolic networks [3], social networks [4], software network [5], scientific collaboration networks [6], and World Wide Web [7,8].Wireless sensor networks (WSNs), as a particular type of complex networks, are different from other types of complex networks .A large-scale WSN consists of a large number of distributed sensor nodes, and therefore a large-scale WSN has some complex features of complex networks [9].Due to short communication range and energy constraints, these sensor nodes autonomously establish connections through wireless communications with the other nodes that are within their local-area.Here the concept of local-area can be seen as a domain surrounding these nodes where their signals can be reached.In fact, there is also localarea in other networks, for example, a community structure in social networks, a regional cooperative group in economic networks, and a domain in the internet, all of which are localareas.
In order to model a WSN, most researches consider that the WSN consists of the same types of nodes and links.Reference [10] shows that it is possible to model a WSN as a small world network by using multiple cluster nodes that can emulate the long edge required in small world networks.In [11], a novel evolving network model is based on random walk to study the fault tolerance of the WSN due to node failure and discuss the spreading dynamic behavior of viruses in the evolution model.Considering the diversities of nodes and links in a real WSN, [9] proposes a local-world heterogeneous model for the WSN.In [9], nodes are divided into two classes: sensor nodes (SNs) and cluster-head nodes (CHs), and there are two kinds of links: the bidirectional edge between CHs (CH↔CH) and the directional link pointing to CHs from SNs (SN → CH).The model in [9] is different from the aforementioned models with homogeneous nodes and links, and it balances energy consumption by limiting the connectivity of CHs to prolong the lifetime of the network.

Mathematical Problems in Engineering
However, recent advancements in wireless communications and sensor technologies have enabled most SNs to relay data packets [12]; that is, there are always SN↔SN in many practical applications.Reference [9] only considers two kinds of links and neglects the links among SNs, which makes the model in [9] be not suitable for most current practical applications.
On the other hand, in practical WSN applications, SNs are usually deployed in some inaccessible and dangerous environments to gather information from the physical world to a sink node (base station); for example, WSNs can be used to monitor environmental changes (such as weather, gas in coal mine) [13,14], monitor habitats [15,16], track objects [17], manage disasters [18], and so on.Most recent researches about WSNs try to improve energy efficiency and prolong the lifetime of the network.Several international research projects dedicated to energy-efficient (EE) wireless communications have been carried out.In [19], a linkadaptive transmission scheme for MIMO-OFDM systems is proposed, which maximizes EE in terms of bits-per-Joule using dynamic power allocation based on the channel state as well as the circuit power consumption.Based on the complex network theory, a new topological evolving model is proposed in.In the evolution of the topology of sensor networks, the energy-aware mechanism is taken into account, and the phenomenon of change of the link and node in the network is discussed.Although they are energy-efficient and have long lifetime, they are not suitable for some applications with the requirement of high reliability (this is critically important in some study of the WSN [20], and the concept of reliability will be given in the next section).Reference [21] shows that, in WSNs, due to environment noise, reliable data communications cannot be definitely achieved.Even under ideal conditions, the packet loss rate of a WSN may be above 1% or close to 1% due to packet collisions.When the communication environment becomes hostile, the packet loss rate will definitely go up, and a successful data transmission over links in the WSN can be guaranteed with a certain probability (less than 100%).As most applications of WSNs need to aggregate sensed data from environment [20], and if we want to aggregate the sensed data in the network, just a part of sensed data from SNs will be sent to the sink node, the aggregation result provided by the WSN cannot be ensured to be accurate, and there is a need to study this issue for WSN's some aggregation operators.As we define the reliability to describe this issue in our research, we call this issue as the reliability of WSN.Even though some researches such as [9] can prolong the lifetime of WSNs through balancing energy consumption, the reliability of the WSN should also be taken into account for reliable data communications.
For the above reasons, a reliability-oriented network model which is consistent with the actual network is very important to the WSN application development.In this paper, we propose a reliability-oriented local-area model for the WSN with the help of growth and preferential attachment mechanism.The influence of the local-area scale , the scale of WSNs, the strategy of evolution on the reliability of WSNs, and topological properties are investigated.Firstly, to measure the reliability of the WSN, we present the concept of reliability and verify its mathematical rationality and find that the reliability  of a large-scale WSN can be represented as a function of the maximum value  of cluster layer of SNs.Secondly, we show that this model has intrinsic characteristics, the assortative connectivity correlation, and characteristic of power law.The energy efficiency of this model is better than that of [9], and the lifetime is identical and it has higher reliability compared with the model in [9].Moreover, from the perspective of the reliability of the WSN, we reveal that the reliability of the WSN evolving from ROLM is much better than that of [9], and when the network scale is larger than 1000, the reliability of the WSN evolving from [9] is far less than that of the ROLM.

Model of WSNs
There are two kinds of nodes (SNs and CHs, which will be identified by IEEE address [22]) in the WSN: the SNs will connect with a CH or other SNs, and the CHs can connect with the sink node or other CHs.These two kinds of nodes perform different functions in the WSN evolving from ROLM, such that SN is responsible for sensing information from geographical environment and sending its sensed data to a CH which it belongs to, while a CH collects sensed data from its cluster members.After processing the whole data it will retransmit the aggregation result to the next hop CH.Considering the complex futures of WSNs, ROLM is designed based on the existing complex system modeling method.Take BA model, for example, [23].There is one new-incoming node entering a preexisting network at every step and choosing one or some nodes from the preexisting network to connect with a certain probability.After a node connects to the network, there will be another new-incoming node entering the network in the next step until the number of nodes in this network reaches the preset network scale.
As a special kind of complex networks, WSNs have some specific characteristics.Therefore, we should consider some factors such as the node transmission range, the hopconstraints between nodes, the reliability of network, and transmission delay.In WSNs, because of the constraint of node transmission range, each node in the network can only communicate with those nodes located within its coverage, which is named as local-area connections.On the other hand, when the monitoring area has the similar condition, every wireless link (the link between two nodes which has one hop between them) in the network has the same packet loss probability, which means that the sensed data will be sent successfully over one link with a probability  ( ≤ 1) [24].And apparently, if the sensed data has to be delivered to a node two hops away, then it will be over two links and the successful transmission probability of this data is  2 .As the probability  is equal to or less than 1, the more hops the sensed data needs to be delivered, the lower the probability that the sensed data will be sent successfully to the sink node.Therefore, just a part of sensed data will be sent to the sink node when the WSN performing an aggregation operator and the reliability of the WSN for this aggregation operator will decrease consequently.ROLM tries to reduce the number of links between two connected nodes to improve the reliability.Based on the local-area and fitness models in complex networks theory, this paper studies the evolution model while considering the constraints of transmission range and the reliability of the WSN in a novel approach.

Problem Definition.
To facilitate researching on ROLM and the reliability of the WSN evolving from ROLM in this paper, we will briefly introduce some related definitions as follows in this part.It includes sink layer, cluster layer, round, -estimate, reliability, , unbiased estimate, and round.
Definition 1 (sink layer).For node  (CH or SN), if its father node is the sink node, then we define that its sink layer is 1 and express it as   = 1; if its father node 's   = 1, then its   = 2, and so on.Definition 2 (cluster layer).For a SN i , if its father node is a CH, then we define that its cluster layer is 1 and express it as    = 1; if its father node 's    = 1, then its    = 2, and so on.There are many types of aggregation result (Sum, Average, Max, and so on) which we can get from the WSNs [25].In order to specifically study the reliability of the WSN evolving from ROLM, we take the SUM aggregation operator as an example and propose the concept of reliability and  parameter to measure the reliability of the WSN.For the other aggregation operator, their processing methods are the same as that of SUM aggregation operator, such that we can get the Average result from the process of calculating the Sum result (the sum value in every cluster divided by the number of members of this cluster, and then the CH will transfer their quotient values to the sink node).
At time , if the WSN can only aggregate a part of sensed data from its nodes, it will provide an approximated SUM, expressed in Sum(  )  .If it aggregates all sensed data in the WSN, it will get an accurate SUM, called Sum(  ).
Definition 7 (unbiased estimate).Sum(  )  is an unbiased estimate of Sum(  ) if the mathematical expectation of Sum(  )  is equal to Sum(  ); that is, (Sum(  )  ) = Sum(  ); otherwise Sum(  )  is a biased estimator of Sum(  ).Definition 8 (round).The amount of time from the sink node initiates a data collection command to all sensed data from SNs being sent to the sink node, no matter how long it will last.This amount of time is called one round.

Reliability-Oriented Local-Area Model.
In this section, we model a WSN as a network with growth and preferential attachment and propose an evolution model ROLM with high reliability.Different from [9], the WSN evolving from ROLM includes three kinds of links.They are the bidirectional link between CHs (CH↔CH), the directional link from SN to CH (SN → CH), and the bidirectional link between SNs (SN↔ SN).For CH i , let    be the current degree that is the total number of links pointing to CH i from other nodes, and let    be its initial energy value.And for a SN i , provide that    is the degree that is the total number of links pointing to SN i from other SNs and The basic settings for ROLM can be described as shown in Algorithm 1.
The growth in the WSN evolving from ROLM is starting with an existing network with a small number of  0 nodes (all of these  0 nodes are connected with each other.As SNs connect to the sink node by CHs, there must be at least one CH among these  0 nodes).Referring to the BA modeling method in complex networks, ROLM can be described as follows. ( where    and    are the residual energy of CH i and SN i . Considering the fact that every link connected to node  will consume some energy of node , we simply define where  is a constant and can be regarded as the energy consumed by every link.Formulas ( 4) and ( 5) mean that the closer the degree of node  is to   max  or   max  , the less the residual energy of the node  will be left.Replacing    and    with ( 4) and ( 5) in formulas ( 2) and ( 3), respectively, we can get (c) After step (b), when a new-incoming SN enters into the preexisting network, and if the degree of an existing SN i is    <   max  , return to step (b); otherwise remove SN i from the local-area Φ; that is to say, a new node cannot connect to this SN i even though SN i is within the local-area Φ.When a new-incoming node (CH or SN) enters into the pre-existing network, if the degree of CH i is    < (  max  − 1) (it could avoid all these   max  nodes connected to one CH i are SNs, and the new-incoming CH cannot joint into the pre-existing network as all CHs within the localarea Φ have reached their saturated degrees), and the new-incoming node is a SN, then return to step (b), otherwise remove CH i from the local-area Φ; if the degree of CH i is    <   max  , and the new-incoming node is a CH, then return to step (b), otherwise remove CH i from the local-area Φ.
The description of ROLM with pseudocode is shown as in Pseudocode 1.
For a WSN with a certain number of nodes, through the above three steps, the nodes in the WSN can be connected with each other.When the number of nodes is 5000,  max  = 7,  max  = 20, the topology graph is shown in Figure 1(a), and Figure 1(b) is an example for detailed topology of one cluster in Figure 1(a).The three steps of ROLM show that the probability of a new-incoming node connecting to a preexisting CH i or SN i not only depends on its current state:    (or    ) and    (or    ) but also has the constraint of    <   max  ≤  max  (or    <   max  ≤  max  ), and    (   ) depends on the node 's initial energy.ROLM has considered the diversities of nodes and links and balanced the energy consumption globally by limiting the links to a CH and the links to a father SN.Furthermore, it also improves the reliability of WSNs for an aggregation operator (the smaller the    is, the higher the reliability of WSNs will be, and it will be proved in the next section for SUM aggregation operator) without shortening the lifetime of the network.As high-reliability and energy efficiency are both critically important for most WSNs applications [26,27], the WSN evolving from ROLM can perform more efficiently under such environment.

Mathematic Principle of 𝜂
Since this paper is dedicated to present an evolution model ROLM to improve the reliability of WSNs, in the following sections we analytically calculate the reliability of the WSN evolving from ROLM for SUM aggregation operator.Section 2.2 shows that, when the number of network nodes is given, through the three steps in Section 2.2, we can connect these nodes to form a complete WSN.Here, we provide that As the WSN for practical applications is not ideal deterministic network model, a data transmission over the wireless link of the real WSN is successfully conducted with a certain probability [28].Here, when the sink node initiating one time of data collecting at time , we provide that a link SN → CH or SN↔SN is connected successfully with probability .In our model, CHs' initial energy is larger; we define that the timeout retransmission mechanism is used among CHs and sink node, and therefore data communications over CH↔CH links can be successfully conducted with probability 1.
According to the above analysis, at time , the sensed data of SNs with    = 1 will be sent successfully to the sink node with a certain probability .Similarly, the sensed data from SNs with    = 2, SNs with    = 3, . .., and SNs with    =  will be sent successfully to the sink node with the probabilities  2 ,  3 , . . .,   , respectively.To facilitate our research on data analysis, let  () ,  ( 2 ) ,  ( 3 ) , . . .,  (  ) be the set of sensed data that are sent successfully to the sink node from the SNs with    = 1, SNs with    = 2, SNs with    = 3, . .., and SNs with    =  at time  respectively.In order to deliver the sensed data successfully to the sink node, the medium access schedule of our network model is similar to that in [29].
In order to specifically study the reliability of the WSN, we take the SUM aggregation operator as an example (the other operator also can be got with the same processing method) and propose the concept of reliability and  parameter to measure the reliability of the WSN.Here, the accurate SUM of the WSN at time  is defined as Sum(  ) = ∑ =1   , .In the next section, we will analyze the mathematic principle of  and study whether the approximated SUM can replace the accurate SUM, which is one of the basic requirements of researching on the reliability of the WSN.It also determines whether we need to measure the reliability of the WSN.After that, we will research on how to calculate the value of  and use it to measure the reliability of the WSN.

Estimator of
,2 According to the definition of unbiased estimate in Section 2.1, the following Theorem 9 will show that Sum(  )  is the unbiased estimator of Sum(  ).Theorem 9. Let (Sum (  )  ) be the expectation of Sum (  )  , and let Var(Sum (  )  ) be the variance.Then, Theorem 9 will be proved in Appendix A. It shows that the mathematic estimator of Sum (  )  is the unbiased estimator of Sum (  ) and the upper bound of Var(Sum (  )  ) is inversely proportional to .That is to say, with the increase of , the upper bound of Var(Sum (  )  ) can be arbitrarily small.Reference [30] shows that, with  increasing, the relative error between Sum (  )  and Sum (  ) gradually decreases, and if  is sufficiently large, this relative error can be arbitrarily small.For any , let the variable  V (1 ≤ V ≤ ) denote the following equations: . . .
There is Sum( . Firstly, we need to prove that Sum(  )  follows normal distribution.In view of the linear combination of  independent normal distribution variables still follow normal distribution, through proving that the sum of sensed data of SNs with    = 1, the sum of sensed data of SNs with    = 2, . .., and the sum of sensed data of SNs with    =  all follow normal distribution to prove that Sum(  )  follows normal distribution.Reference [31] shows that, if sensed data from SNs with    = V (1 ≤ V ≤ ) is conformed to Lyapunov condition, the sum of sensed data from SNs with    = V will meet the application conditions of central limit theorem; that is, the sum of sensed data from SNs with    = V will follow normal distribution.And Theorem 10 proves that sensed data from SNs with    = 1, SNs with    = 2, . .., and SNs with    =  are all conformed to Lyapunov condition, respectively.Theorem 10.These  groups of sequence of random variables  V (1 ≤ V ≤ ) satisfy Lyapunov condition; that is, ∃ V > 0 satisfies the following: where 1 ≤ V ≤ ,  V is the number of sensed data from SNs with cluster layer V at time , and =1  V , and for all  (1 ≤  ≤  V ) there are  V = ( V ) and  V = Var( V ).
The proof for Theorem 10 will be given in Appendix B [31].Theorem 10 shows that these  groups of random variable sequences  V (1 ≤ V ≤ ) satisfy Lyapunov conditions.That is, the sum of sensed data from SNs with =1  V follows normal distribution.As the successful transmissions of the sensed data from SNs with    = V to the sink node are independent of other SNs with    , Sum(  )  is the sum of these  independent variable normal distributions Sum( ,V )  .Thus Sum(  )  follows normal distribution.For a given relative error limit , the following Theorem 11 describes how to calculate the  of the WSN.Theorem 11.Define  = 1 − , and  /2 is the /2 quantile of standardized normal distribution, if  /2 satisfies Then, the probability that the relative error between Sum(  )  and Sum(  ) satisfies the given error limit  will be equal to or greater than ; that is, We will prove Theorem 11 in Appendix C [32].In Theorem 11,  is an important parameter which could measure the reliability of the WSN in our research.Combining with the above analysis, we can conclude that the problem of the calculation of  can be described as follows.

The Degree Distribution of Nodes
To achieve a better understanding of the considered model, in this section, ROLM will be analyzed in comparison with the model in [9].When we apply the preferential attachment strategy of [9] to the network evolution, the new-incoming node (CH or SN) will connect to CH i in the local-area Φ according to the probability ∏   , and SN i in the local-area Φ will be connected by a new-incoming SN according to the probability ∏   : Here  is a constant and can be regarded as the energy consumed by every link.Considering that a small number of CHs would make lots of SNs connected to them and large number of CHs will waste the network resource, this paper will control the number of CHs in a low proportion.In order to compare these two evolution models fairly, in our experiment, all parameters are set in the same way in [9], and set the number of WSN nodes to be 5000,  max  = 7,  max  = 20,  = 0.1, (  ) ∼ [0.7,1], (  ) ∼ [0.3,0.7].Then we will use the evolution strategy of ROLM and the evolution strategy of [9] to connect the nodes and define that there are two kinds of nodes and three kinds of links in this WSN.Next, we will focus on the cumulative degree distribution log 10 (()) of nodes in WSNs evolving from ROLM and [9].To find the cumulative degree distribution log 10 (()), we first need to do degree statistics (i.e., get the probability of a node with k links).The statistical result is shown as in Figure 2.
Figure 2 shows the statistical results of  = 25,  = 35, and  = 45 for WSNs evolving from ROLM and [9], respectively.The left and right three figures in Figure 2 indicate that the change of local-area scale  does not cause the cumulative degree distribution to change greatly.For the WSN evolving from ROLM, the slopes of fitting curves for these three  are within [−2.5, −2].And for [9], the slopes of fitting curves are within [−3.6, −3.2].That is, in these two evolution models, the scale of local-area  has a weak effect on the cumulative degree distribution.
On the other hand, the six constant variables of the above six fitting curves exhibit that the power-law fitting of the WSN evolving from ROLM is better than that of [9], the WSN evolving from ROLM follows a power-law distribution that can be found generally in the so-called scale-free networks, and the power-law features of the WSN evolving from ROLM are more obvious than that of [9].The six figures in Figure 2 could also reveal the proportion of the nodes (SNs who are father nodes and CHs) with high degree and low degree.There are a low proportion of nodes with high degree and the proportion of nodes with low degree is large.And the degree of most of the nodes    (or    ) is far smaller than their maximum degree  max  (or  max  ), which could prevent the CHs and father SNs from consuming energy too quickly; then the lifetime of the WSN would be prolonged.

The Discussion of Network Performances
In this section, we will discuss the performance of WSNs evolving from ROLM and [9].As our purpose is to give an evolution model which is reliability-oriented, we will first analyze the reliability of WSNs evolving from ROLM and [9].Considering that the advantage of [9] is prolonging the lifetime of WSN, we also analyze the lifetime of WSNs evolving from ROLM and [9] in this section.Finally, we study the average length of shortest path of these two WSNs.As the other parameters of formula (11) are decided by the actual use of the WSN, and when we compare the reliability of the WSN evolving from ROLM with that of [9], these parameters have no effects on the comparative results if we set them in the same way.Therefore, we just choose an example setting for these experiments in this part.And for different network scales and , we can get corresponding WSNs evolving from ROLM and [9] and the maximum value of cluster layer of SNs for these networks.Then, in combination with Section 3, we can calculate the reliability of the corresponding networks with different network scale and .The results are shown in Figure 3 (for  = 0.85) and Figure 4 (for  = 0.9).Firstly, the four figures in Figures 3 and 4 show that, with the increasing of network scale, the reliability of network presents a downward trend (not monotonously decreasing).As the CHs enter into the preexisting network in probability  and the new-incoming nodes will link to CH i or SN i in the local-area Φ in accordance with the probability ∏   or ∏   in the WSN evolving from ROLM and ∏   or ∏   in the WSN evolving from [9], WSNs evolving from ROLM or [9] are not unique.That is, every evolution may end with a different network, and this leads to the fluctuating in Figures 3 and 4. Secondly, for the same  and network scale, the WSN evolving from ROLM will get a higher reliability than that of [9], such that the reliability of the WSN evolving from ROLM is 0.857 when   = 1000,  = 0.85, while the reliability of [9] is only 0.0011.This is due to the difference between the preferential attachment principles of ROLM and [9].Reference [9] only considers balancing energy consumption globally to prolong the lifetime of the WSN, while ROLM not only considers balancing energy consumption globally but also tries to improve the reliability of data aggregation.The experiment results shown in Figures 3 and 4 indicate that ROLM can improve the reliability of WSN significantly.Lastly, experiment results in Figures 3(a) and 4(a) show that ROLM can guarantee the WSN a higher reliability for smallscale WSN and large-scale WSN.Even under the conditions where   = 20000,  = 0.85, the reliability of the WSN evolving from ROLM can maintain about 0.6399.

Network Lifetime of ROLM.
After discussing the reliability of the WSN evolving from ROLM, we will analyze the network lifetime of the WSN evolving from ROLM in this part.Same as Section 5.1, here we will compare the lifetime of the WSN evolving from ROLM with that of [9].According to the definition of round in Section 2.1, in order to facilitate the statistics of round, we define that all nodes will aggregate its own sensed data and the sensed data from other nodes with fixed-length  bits and then send it to its father node.According to [34], for one SN, the energy expenditure of sending 1-bit data is equal to that of executing 1000 commands.Therefore, here we just calculate the communication energy expenditure and ignore that of executing overhead.In this group's experiments, assuming a node expends  energy to send a  bits packet, we set  =  max  /500.That is, the energy of each node decreases by  in every round.If a node has less than  energy, then we consider that the node and all its son nodes are dead (the sensed data of its son nodes need it to relay their sensed data).Except from the   = 5000, the rest of other variables are set as in Section 5.1.Under the above assumptions, the number of nodes alive will decrease with the increase of network running rounds, and the comparing results of lifetimes of WSNs evolving from ROLM and [9] are shown in Figure 5.
There are three groups of number in Figure 5; they are network lifetimes in different network scale of 1000, 3000, and 5000.From Figures 5(a), 5(c), and 5(e), we find that, under the same assumptions, WSNs evolving from ROLM and [9] have the same lifetime, the network lifetime of both WSNs evolving from ROLM and [9] is 438, and their initial time of attenuation is 125.That is, some nodes begin to die from the 125th round, and all nodes are dead in the 438th round.Both of them have set the same  max  and  max  to prolong the lifetime of network. Figures ) are detailed graphs of one part from Figures 5(a), 5(c), and 5(e), respectively.These three figures reveal that the number of nodes alive in the WSN evolving from [9] decreased quickly and the number of death nodes in the 125th round is much larger than that of ROLM.For the network with scale of 1000, the number of alive nodes in the WSN evolving from ROLM is larger than or equal to that of [9] among the rounds within [1,175].And for the network with scale of 5000, the number of alive nodes in the WSN evolving from ROLM is larger than or equal to that of [9] among the rounds within [1,200].That is, the WSN evolving from ROLM is more energy-efficient than that of [9].This is because [9] only considers the energy consumption of nodes and did not control the maximum value of nodes' cluster layer, which leads to the result that the new-incoming nodes have a tendency to link to the SNs that are far away from the sink node which can increase the hops among any two nodes and energy consumption for processing data and finally decrease the energy efficient.
The above two groups of experiments about the network reliability and the network lifetime prove that ROLM can ensure a higher reliability which is much larger than that of [9], ROLM can be used in larger-scale network where [9] could not, and the energy efficiency of the WSN evolving from ROLM among some rounds is better than that of [9].

⟨𝐿⟩-
The Average Length of Shortest Path.In the WSN, SNs collect and send data to their father nodes (SNs or CHs) and then forward them to the sink node via multihops, so all SNs or CHs who act as a father node play an important role in data  The scale of network N t = [0, 1000], q = 0.9 (b) The changing trend of [9] reliability Figure 4: The changing trend of reliability with the increase of network scale ( = 0.9).collecting and transmitting.We use ⟨⟩-the average length of shortest path to measure the average hops among SNs and CHs for data processing; its calculation can be written as [23] where   is the length of shortest path between any two nodes (SNs or CHs)  and ,  is the number of nodes in the network, and   = +∞ if there is no path between nodes  and .Different from [9], considering that every path in the network can reflect the efficiency of the network, we use the   between every two nodes in the network to calculate the value of ⟨⟩.The smaller ⟨⟩ means fewer hops and less energy consumption for data processing.On the other hand, [9] just analyzed the relationships among , , and ⟨⟩, and these can reflect the influences of  and  on ⟨⟩ but cannot reflect the efficiency of the evolution models.Therefore, in this part, we study the relationship among the scale of network , ⟨⟩, and  ( = 25 and   6).
From Figure 6, we can find that, for each , ⟨⟩ of the WSN evolving from ROLM is smaller than that of [9] no matter whether  = 25 (shown as in Figure 6(a)) or  = 45 (shown as in Figure 6(b)), such that when  = 2000 and  = 25, there are ⟨⟩ = 9.587719786 for the WSN evolving from ROLM and ⟨⟩ = 42.89694797 for the WSN evolving from [9].When  = 2000 and  = 45, there are ⟨⟩ = 9.441435432 for the WSN evolving from ROLM and ⟨⟩ = 59.72642678 for the WSN evolving from [9].That means the number of hops among two nodes in the WSN evolving from ROLM is less than that in [9] which can ensure a better synchronization, and the energy consumption for data processing in the WSN evolving from ROLM is less than that of [9].It also explains the results shown in Figure 5 in which the WSN evolving from ROLM has better energy efficiency than that of [9].That is because we try to reduce the cluster layer in the WSN evolving from ROLM while [9] does not consider it.
The two blue lines in Figure 6 indicate that the slopes of blue lines decrease with the increase of  (when  > 1000), and the red lines do not have this phenomenon.That is, the WSN evolving from ROLM has more obvious small world characteristic than that of [9].For different values of C, the change extent of ⟨⟩ in the WSN evolving from ROLM is weaker than that in the WSN evolving from [9], which means that  has a weak effect on the WSN evolving from ROLM.Because [9] only considers the energy of nodes to evolve the WSN, for big , the new-incoming node will tend to connect to the preexisting node whose residual energy is larger (no matter how many its cluster layers), and the change of  affects the ⟨⟩ of the WSN evolving from [9].
The above analysis and Section 4 show that the WSN evolving from ROLM has steadier complex network features than that of [9].

Conclusions
Consider large-scale WSNs have showed many complex features and an appropriate evolution model is very necessary for improving the performances of WSNs.In this paper, we are trying to build an evolution model which could evolve a WSN with high reliability and find out a way to measure the reliability of network.Describing the WSN as a model with two kinds of nodes and three kinds of links, we proposed a reliability-oriented large-scale model (ROLM) to improve the reliability of the WSN for some aggregation operators and presented the reliability parameter  to measure the reliability of network for SUM aggregation operator.After proving the mathematical rationality of parameter  theoretically and deriving the calculating method of , we use  as one of performance metrics of WSNs.
The analysis of degree distribution of nodes shows that both of WSNs evolving from ROLM and [9] follow a powerlaw distribution and ROLM has a more obvious scale-free characteristic.From the experiments, we can find that, for a specific network scale, the reliability of the WSN evolving from ROLM is higher than that of [9]; for example, when  = 0.85, the evolution model in [9] is unusable when the network scale is larger than 1000 nodes, while ROLM can ensure the reliability of the WSN to be 0.6399 even with 20000 nodes in the network.On the other hand, the experiment Mathematical Problems in Engineering results also show that ROLM balanced energy consumption by limiting the connectivity of nodes to prolong the life of the network, and the WSN evolving from ROLM achieves the same long lifetime with that of [9].Within many rounds, the energy efficiency of the WSN evolving from ROLM is better than that of [9].Lastly, we also give experiments to analyze the average length of shortest path of WSNs evolving from ROLM and [9].The results indicate that the hops among two nodes in the WSN evolving from ROLM are fewer than those in [9] which can ensure a better synchronization, and the energy consumption for data processing in the WSN evolving from ROLM is less than that of [9].The experiments about the degree distribution and the average length of shortest path show that the WSN evolving from ROLM has steadier complex network features than that of [9].And the experiments about the network reliabilities and lifetimes prove that the energy efficiency of the WSN evolving from ROLM within some rounds is better than that of [9], ROLM can ensure a higher reliability which is much larger than that of [9], and ROLM can be used in larger-scale network where [9] could not be.
Without considering the analysis of other    (their researching methods are the same as V = 1), we next describe the analysis of SNs with    = ; that is, V = .

C. The Proof for Theorem 11
From formula (11)

Figure 1 :
Figure1: The topology graph and subgraph of topology of the WSN, where nodes' sizes represent their current degree, the nodes' color reflects their community structure, nodes with the same color means that they have the same modularity, and the communication among them will be more frequent.

3. 2 .
Calculation of .The steps of calculating the value of  are (a) proving that Sum(  )  follows normal distribution; (b) transforming the normal distribution into standard normal distribution; (c) utilizing characteristics of standard normal distribution to calculate the value of .

Figure 3 :
Figure 3: The changing trend of reliability with the increase of network scale ( = 0.85).
The concepts of   and    show that CH i has only value of   ; and SN j has both   and    , and the values of   and    for SN j are different.
represents the initial energy value of SN i .Define that    is a random variable following the uniform distribution (  ) in the interval [ min  ,  max  ] and    is a random variable following the uniform distribution (  ) in the interval [ min  ,  max  ].Set    as much larger than    .Here,  max  ≥  min  ≥  max  ≥  min  ,  max  is the maximum initial energy of CHs, and a CH with  maxc could connect with less than  max  nodes.Similarly,  max  is the maximum initial energy of SNs, and the SN with  maxs could connect with less than  max  SNs.Let   max  and   max  be the saturated degree of CH i and SN i respectively; then we define //starting with an existing network with a small number of  0 nodes and at least one CH node among them //using the Function ROLM Model to construct the WSN Function ROLM Model Begain For  = 1 to NumberOf Networks do { //Firstly, we mark  nodes randomly from the pre-existing network as the local-area Φ local-area Φ:  nodes,  ≥ max{ max  ,  max  }   is the residual energy of SN  ,   and   are  be the number of SNs in the WSN at time .Let   (1 ≤  ≤   ) be the sensed data of SN i at time , and let   = { 1 ,  2 , . . .} be the set of all sensed data in the WSN at time .Since the value of sensed data is bounded, we use sup(  ) to denote the upper bound of all sensed data, and the metric of data depends on the category of sensors.From the analysis in Section 2.2, except the sink node and the links between CHs and sink node, we can construct a WSN consisting of two kinds of nodes (CHs and SNs) and three kinds of links (CH↔CH, SN → CH, and SN↔SN).The WSN is divided into a number of clusters which are disjoined with each other, and a cluster consists of one CH and many SNs.Here, we define that the WSN is divided into  clusters.Let  = max{   , 1 ≤  ≤ (1 − )  } be the maximum value of cluster layer of SNs at time , and let  1 ,  2 ,  3 , . . .,   be the numbers of SNs whose   } be the set of sensed data from SNs with    = 1 at time , and let  ,2 ,  ,3 , . . .,  , be the sets of sensed data from SNs with     ,  ,1 ,  ,2 ,  ,3 , . . .,  , can be described as  ,1 ⊆   ,  ,2 ⊆   ,  ,3 ⊆   . ..,  , ⊆   ,   = { ,1 ,  ,2 ,  ,3 . . .,  , }.