Characterizing the Path Coverage of Random Wireless Sensor Networks

Abstract

Wireless sensor networks are widely used in security monitoring applications to sense and report specific activities in a field. In path coverage, for example, the network is in charge of monitoring a path and discovering any intruder trying to cross it. In this paper, we investigate the path coverage properties of a randomly deployed wireless sensor network when the number of sensors and also the length of the path are finite. As a consequence, Boolean model, which has been widely used previously, is not applicable. Using results from geometric probability, we determine the probability of full path coverage, distribution of the number of uncovered gaps over the path, and the probability of having no uncovered gaps larger than a specific size. We also find the cumulative distribution function (cdf) of the covered part of the path. Based on our results on the probability of full path coverage, we derive a tight upper bound for the number of nodes guaranteeing the full path coverage with a desired reliability. Through computer simulations, it is verified that for networks with nonasymptotic size, our analysis is accurate where the Boolean model can be inaccurate.

1. Introduction

Wireless sensor networks (WSNs) have many applications in security monitoring. In such applications, since it is essential to keep track of all activities within the field, network coverage is of great importance and must be considered in the network design stage.

Path coverage is one of the monitoring examples, where WSNs are deployed to sense a specific path and report possible efforts made by intruders to cross it. In a manual network deployment, the desired level of the path coverage can be achieved by proper placement of the sensors over the area. When it is not possible to deploy the network manually, random deployment, for example, dropping sensors from an aircraft, is used. Due to the randomness of the sensors location, network coverage expresses a stochastic behavior and the desired (full) path coverage is not guaranteed. Thus, a detailed analysis of the random network coverage can be ultimately useful in the network design stage to determine the node density for achieving the desired area/path coverage.

Path coverage by a random network (or barrier coverage which is a relaxed version of the path coverage) has been the focus of some previous work [1–6]. In [1], assuming that a random network is deployed over an infinite area with nodes following a Poisson distribution, authors investigate the path coverage of the network. They first study the path coverage over an infinite straight line when the sensor has a random sensing range. Then, they show that in the asymptotic situation, where the sensing range of the sensors tends to 0 and the node density approaches infinity, the results are extendible to finite linear and curvilinear paths. Further, a path coverage analysis is proposed for a high-density Poisson-distributed network in [2] where sensors have a fixed sensing range. The path coverage analysis of [1, 2] is based on the Boolean model of [7], where a Poisson point process is justified.

Kumar et al. study -barrier coverage provided by a random WSN in [3]. To this end, they develop a theoretical model revealing the behavior of the network coverage over a long narrow belt. It is assumed that the sensors are spread over the belt according to a Poisson distribution. The authors propose an algorithm determining whether an area is -barrier covered or not. Also, they introduce the concepts of weak and strong barrier coverage over the belt and derive the condition on the sensors density guaranteeing the weak barrier coverage.

The focus of [4] is on the strong barrier coverage. First, authors present a condition insuring the strong barrier coverage over a strip where the sensors locations follow a Poisson point process. Then, by considering asymptotic situation (on the network size and number of nodes) and using Percolation theory [8], they determine, with a probability approaching 1, whether the network has a strong barrier coverage or not. Then, they use their analysis to devise a distributed algorithm to build strong barrier coverage over the strip.

In this work, unlike most existing studies which focus on asymptotic setups, we study the path coverage of a finite random network (in terms of both network size and the number of nodes). As a result, the Boolean model is not accurate. Alternatively, the methodology of this work is based on some results from geometric probability. Our focus is on the path coverage for a circle, but extension to other path shapes is briefly discussed.

In the ideal case, all sensors are located exactly on the path. This, however, is not a practical assumption for randomly deployed networks. To consider the inaccuracy of the sensors locations, we assume that sensors are inside a ring containing the circular path. As a result, the portion of the path covered by any given sensor is not deterministic. Moreover, other factors may affect the sensing range of a sensor. Thus, our analysis is not based on a fixed sensing range. Indeed, we first develop a random model for the covered segment of the path by each sensor. Then, we study the distribution of the number of uncovered gaps on the path. The full path coverage is a special case where the number of gaps is zero. This is used to determine a tight bound on the number of active sensors assuring the full path coverage with a desired reliability. Also, we find the probability of having all possible gaps smaller than a given size. This probability reflects the reliability of detecting an intruding object with a known size.

In addition to studying the number of gaps, we present a simplified analysis for deriving the cumulative distribution function (cdf) of the covered part of the path. This simplified analysis is based on using the expected value of the covered part of the path by a sensor instead of considering the precise random model. We observe that the simplified analysis can provide a fairly accurate approximation of the path coverage.

Since our analysis studies the effect of the number of nodes on the path coverage of a finite size network, it can readily be used in the design of practical networks. In fact, using our results, one can determine the number of nodes in the network to satisfy a desired level of coverage. An example is provided.

The paper is organized as follows. Section 2 introduces the network model and defines the problem. Our coverage analysis is presented in Section 3. Section 4 includes computer simulations verifying our analysis. Finally, Section 5 concludes the paper.

2. Preliminary

2.1. Network Model

We consider sensors monitoring a circular path with unit circumference, called . In an ideal case, the sensors are precisely located on the circular path, but this is not usually true for a randomly deployed network. In order to take this fact into account, we assume that sensors are randomly spread over a ring containing (See Figure 1). We assume a symmetric distribution for sensors, that is, the sensor density does not depend on the polar angle and is determined only by the distance from the center. It is generally desired to have more sensors in the vicinity of Thus, distributions with larger values close to are preferred. When no effort is made to put the sensor as close as possible to the path ( sensors are spread totally randomly), the uniform distribution is obtained. Hence, in the sense of placement efforts, uniform distribution reflects the worst case. We consider uniform distribution to verify our analysis by computer simulation in Section 4. Our analysis, however, is presented for any given symmetric distribution. Also, notice that since the number of sensors is finite and known, Poisson distribution, which has been the focus of existing asymptotic analysis, is not applicable.

Figure 1

Network structure.

We also assume that sensors sensing range may vary from to . Obviously, for a fixed sensing range, model . Without loss of generality, it is assumed that the width of the ring is smaller than or equal to and the desired circular path is located at the middle of the ring. Since the sensors farther than to the path do not contribute to the path coverage, our assumption on the ring width does not hurt the generality of the analysis.

2.2. Motivation and Problem Definition

Our goal in this work is to investigate the quality of the path coverage by the network. To this end, we study the following features of the path coverage which provide us with a concrete insight to the performance of the network.

• (i) Distribution of the number of gaps. Due to the randomness of the network implementation, sensors may not cover the whole path. In this case, one or more gaps appear. Assume that represents the number of gaps on . We are interested to find the probability of having gaps, shown by .

• (ii) Full path coverage. It is desired to provide a complete coverage of the path. Since the full path coverage is identical to having no gaps, one can equivalently find . This can simply be found from the derived distribution of .

• (iii) Reliability of the network in detecting objects. It is important to investigate whether the network is able to detect an object, while the path is not fully covered and there may exist some gaps. Basically, we need to consider the size of the gaps in addition to their number. If one knows the size of the intruders beforehand, it is not necessary to provide the full path coverage. Instead, it is possible to deploy a network such that while the path is not fully covered, the size of the possible gaps is smaller than the intruders. Clearly, implementing a network with possible small gaps requires fewer number of nodes and consequently is less expensive. To this end, we find the probability of having all gaps smaller than a given length , denoted by .

• (iv) Distribution of the covered part of the path. The covered part of the path, , has a stochastic nature and its distribution provides a general view of the entire path coverage. In fact, the covered part of the network reflects the combined effect of the number of gaps and their sizes. We derive the cdf of , .

3. Path Coverage Analysis

In this section, we present our analysis of the path coverage. For this purpose, we take advantage of existing results in geometric probability and extend them to our case. After the exact coverage analysis, a less complex approximate analysis is also presented.

An arbitrary point on is covered if it is within the sensing range of at least one sensor. Here, we assume that the sensing area of sensor is a circle denoted by . The covered part of the path by each is its intersection with which is an arc, called . Thus, the total covered part of the path is

(1)

Notice that the length of 's depends on the location of the sensor within the ring-shaped network area and its sensing range. Considering an arbitrary point as the origin on and choosing the clockwise direction as the positive direction, each starts from and continues (clockwise) until , (Figure 2). In other words, determines the most left point of the arc and specifies the most right point of the arc. There are two noteworthy issues here. First, the size of 's and their positions are random because of the random placement of the sensors over the ring. Second, is not necessarily connected and there may exist several uncovered gaps on . The number of uncovered gaps on and their size can reflect the possible opportunities for intruders to pass without being detected by the sensors. If is fully covered.

Figure 2

The random arcs placed clockwise on .

The problem of covering a circle with random arcs has been studied in geometric probability [9–15]. In some cases, it is assumed that the arcs have a fixed length [9, 12, 13, 15] or the analysis is conducted in the asymptotic situation [10, 14]. Asymptotic analysis is suitable for the situation where the sizes of 's are significantly smaller than . In the following, we initially use the result of [11] on the coverage of a circle with random arcs of random sizes. This helps us to provide an exact explanation of the path coverage. Then, we use the mean value of 's to provide a simplified approximate analysis based on the fixed-length random arc placement over the circle [12].

3.1. Exact Analysis

We apply the following theorem from [11] to find the exact distribution of the number of gaps on .

Theorem 1.

Assume that arcs are distributed independently with a uniform distribution over a circle of circumference 1. If shows the cdf of the arc length over , the distribution of the number of uncovered gaps on the circle, , is

(2)

where

(3)

Proof.

See [11].

To apply Theorem 1 for finding the number of uncovered gaps on , we first prove the uniformity of the arc distribution over in the following lemma.

Lemma 1.

For a symmetric distribution of the sensors over the path, the location of the intersection of the sensors sensing range and is uniformly distributed over .

Proof.

We equivalently show that the center points of the arcs are uniformly distributed over the circle. For this purpose, consider a small element with length on . Then, we build a sector of the ring based on this length element whose left and right sides pass the left and right ends of the length element (Figure 3). The center point of the arcs, resulted from the intersection of the sensing area of the nodes within the sector and falls within . Due to the independence of the sensors distribution from the polar angle, all elements with length on the circle have the same chance to include an arc center point. Therefore, the distribution of the arc centers, and consequently arc locations, is uniform on .

Figure 3

Distribution of the arc position over the circle.

Following Lemma 1, in order to find the distribution of the number of gaps on we need or in our case the cdf of 's. Notice that 's are independent and identically distributed (i.i.d) random variables. We find in the appendix for arbitrary distributions of sensor location and sensing range.

As a result of Theorem 1 and Lemma 1, we have the following corollary.

Corollary 1.

The probability of the full path coverage, , is

(4)

Furthermore, one can show that the expected number of gaps on is [11]

(5)

where is the mean of 's.

can be used to find an upper bound on the number of nodes in the network guaranteeing the full path coverage with a given reliability. This is presented below.

Corollary 2.

To guarantee a full path coverage with probability , the following relation holds

(6)

Proof.

Recall Markov's inequality for a positive random variable

(7)

where . If we let be the random variable of the number of gaps , and put , we have

(8)

Combining (5) and (8) results in (6).

Using (6), it is straightforward to find an upper bound on guaranteeing a desired level of coverage, . Later, our simulations show that this bound is in fact very tight.

Another feature of the path coverage that we like to study is the quality of the coverage in terms of the size of the gaps on . Assume that we like to guarantee detecting any object bigger than a particular size, say . To assure detecting such objects, all of the gaps have to be smaller than . Hence, we like to find the probability of having no gaps larger than , where is the length of the largest gap on .

Corollary 3.

The probability of having no gaps larger than is

(9)

Proof.

Consider a realization of the random placement of arcs on the path. Now, one can consider a scenario where the length of each arc is increased by . If there exists a gap smaller than in the first scenario, this gap will be covered in the second scenario since the arcs are longer. On the other hand, a gap with any size in the second scenario will be a gap with length more than in the first scenario. Notice that the above discussion is valid for any realization of the network. Thus, instead of investigating the probability of having no gaps longer than in the first scenario, we look for the probability of the full coverage in the second scenario. Denoting the length of the arcs in the second scenario by , one can think of them as being drawn from the distribution or equivalently

(10)

This completes the proof.

Using the same approach taken for finding the upper bound on in (6), one can derive an upper bound on the number of nodes to guarantee having all gaps smaller than .

We also like to investigate that is, the portion of which is covered by the nodes. To find , we first reorder the arcs based on their starting points, 's. Thus . Now, we divide to arcs , where is an arc starting from and ending at . Finally, starts from and ends at . Since we have random arcs intersecting with , there exist of such spacings on the circle. These spacings may or may not be covered by the network. Adding the covered parts of the path together, we have

(11)

where . Also, 's denotes the length of and

(12)

Notice that in (12) we assume rotational indices for 's. It means that if we replace the index with . In (12), is the length of the connected part of starting from and continuing clockwise. When , the whole spacing is covered and function should return . When , a portion of remains uncovered and there exists a gap at the right side of . Thus, function returns . It is noteworthy that because of the problem symmetry, 's are identically distributed random variables. Thus, we use a single random variable to refer to them.

The distribution of can be well approximated by a Gaussian distribution using central limit theorem (CLT) where the mean value of , , is . Here, denotes the mean value of . Also, where and represent the variance of and , respectively. In reality, one can safely simplify (12) to

(13)

This is because 's are i.i.d. and thus it is very unlikely that, for example, .

3.2. Approximate Analysis

In the following, we present an approximate analysis simplifying our path coverage study. The idea of this approximate analysis is to consider a model where a set of fixed-length arcs are spread randomly over instead of using the actual random-sized arcs. The length of these fixed arcs is equal to the mean value of the random-sized arcs in the original case. We denote the mean value of these random arcs with . In this case, it can be shown that the number of uncovered gaps on is distributed as follows [12]:

(14)

where . The same technique as before is applicable to find the probability of having no gaps larger than on . For this purpose, we just need to use instead of in (14). In addition, the distribution of can be derived when the arc size is fixed [12]. In this case, we have

(15)

where is the cdf of .

One can also calculate the expected value of . To this aim, we first consider the uncovered part of the path, , and find its expected value, called . Then can be found using the fact .

An arbitrary point on remains uncovered when there is no covering it. This is equivalent to having none of 's within an arc with length whose right end point is . There are sensors in the network, hence, the probability of having uncovered, , is

(16)

Consequently,

(17)

3.3. Some Remarks

Remark 1.

Our path coverage analysis is applicable to any closed path, for example, ellipse, with finite length when the location of the path segment covered by an arbitrary sensor is uniformly distributed over the path. For this purpose, we just need to have the distribution of the intersection of sensors sensing range and the path. Also, the analysis is applicable to linear path coverage. In fact, the problem of covering a circle with random arcs can be transformed to the problem of covering a piece of line, say the interval , with random intervals. In this case, sensors are deployed randomly over a strip surrounding the linear path. It is notable that in the linear case, Torus convention [7] is applied. In Torus convention, it is assumed that if a part of the random interval goes out of the line segment, it comes in from the other side of the line piece. However, when the length of the random intervals is small compared to the line piece, one can remove the Torus convention and the analysis remains quite accurate.

Remark 2.

In many WSNs, the number of active sensors in the network changes with time. This can be due to, for example, sleep scheduling or death of some nodes. Since our analysis is provided for arbitrary it can accommodate such situations, simply by replacing with in relevant equations. Consequently, the coverage can be studied as a function of time.

4. Simulation

In this section, we demonstrate the accuracy of our analysis via computer simulations. We have inspected two scenarios for the sensors sensing range. In the first scenario, we assume a network with sensors all having a fixed sensing range equal to . The sensors are uniformly deployed inside a ring around the circular path, where has unit circumference. In the second scenario, the sensors sensing range is also uniformly distributed between and . A zero sensing range can represent a dead sensor.

We evaluate random properties such as the full coverage, number of uncovered gaps, tightness of the bound presented in (6), the intruder detectability, and the portion of the covered path using simulation, and compare the results with our theoretical analysis.

4.1. Uncovered Gaps

Probability density of the number of uncovered gaps on the path, , was derived in Section 3.1. Figure 4 shows the probability mass function (pmf) of the number of uncovered gaps via simulation for . Here, we have assumed that the sensing range of all sensor nodes is fixed and is equal to . The theoretical results using (2) have also been sketched for comparison. It can be concluded that the formulation derived in Section 3.1 quite accurately describes the pmf of the number of uncovered gaps on the path. The third curve in Figure 4 is the result of approximation analysis in Section 3.2. Parameter in (14) is set to be the expected value of random variable , derived in the appendix.

Figure 4

Probability of the number of uncovered gaps when and .

It is clear from Figure 4 that the results from the approximate analysis are fairly close to the exact analysis and the simulations. Due to the complexity of the evaluation of exact analysis, we compare the rest of our simulation results with the approximate analysis presented in Section 3.2 to characterize the coverage properties of the network.

In the case of fixed sensing range, as the width of the ring becomes smaller, the variance of decreases and the arc lengths become closer to , making the approximate analysis more accurate. To study the worst case, in our analysis, we assume the ring width is equal to . Notice that any node outside this ring does not contribute to the path coverage. For random sensing range, we choose . Notice that since , there will be nodes in the ring that will not contribute to the path coverage.

Figures 5 and 6 demonstrate the probability of full coverage versus number of sensors deployed in the region. Figure 5 shows the results for the fixed sensing range scenario. is estimated through simulation for different values of sensing range, , and is compared with the theoretical results using (14). As seen from the curves, our theoretical formulation can effectively predict probability of full path coverage. Figure 6 represents the results of variable sensing range scenario. The sensing range of sensor nodes is randomly selected from a uniform distribution between and . We have run the simulation based on three uniform distributions, , and compared with the theoretical results. For theoretical calculations, we have computed the average arc length for the case of random sensing range using (A.22) in the appendix, and then substituted the resulting into the approximate formula (14). From Figure 6, we can see that the theoretical analysis in Section 3.2 can closely describe the probability of full coverage for random sensing range scenario.

Figure 5

Probability of full coverage, that is, , versus the number of sensor nodes for fixed sensing range.

Figure 6

Probability of full coverage, that is, , versus the number of sensor nodes for random sensing range.

We have also tested our analysis for full coverage of a straight line segment instead of a circle. Figure 7 depicts the coverage of a straight line segment of length when sensing range is fixed and is equal to . The solid line is the result of Poisson assumption in [2] and the dashed line is the result of our formulation. It can be seen that specially for smaller number of sensor nodes, the Boolean model is not well applicable to describe the coverage of small networks.

Figure 7

Probability of full coverage of a straight linear path of length when .

The expected number of uncovered gaps, , after deploying sensors in the ring is given by (5). In Figure 8, has been calculated versus for three values of fixed sensing range, , using simulation as well as the analysis. The expected number of gaps for variable sensing ranges is shown in Figure 9. The sensors sensing range has been taken from the three uniform distributions used previously.

Figure 8

Expected value of the number of gaps versus number of sensor nodes for fixed sensing range.

Figure 9

Expected value of the number of gaps versus number of sensor nodes for random sensing range.

In Section 3.1, we used Markov's inequality to find a relation between the number of nodes and probability of full coverage over the path, presented in (6). The smallest number of nodes that satisfies (6) can efficiently be found by conducting a simple binary search. We denote the value of found via search by . Table 1 shows the value of calculated for probability of full coverage being equal to . The results found by inequality (6) and simulation are shown for comparison. For probability of full coverage closer to one (the region of interest in practice), gets even closer to the value of satisfying the desired reliability found via simulation. For example, for probability of full coverage equal to , we found for various values of It can be inferred that (6) provides a tight upper bound on the number of nodes needed for full path coverage.

Table 1 Upper bound on the number of sensors for guaranteeing full coverage with probability .

4.2. Detectability of the Object

As discussed in Section 3.1, in real applications of sensor networks, we might not be interested in the full coverage of a path, yet we need to make sure that there are no uncovered gaps on the path, larger than a certain maximum length . The probability of this kind of coverage, , was given by (9). We use simulation to find for values of equal to , when . Again, comparing simulation results with theoretical ones in Figure 10 verifies our formulation. Our study on the size of the gaps is useful for decreasing the cost of the network implementation. In fact, if we know the size of the intruders, instead of providing a full path coverage, we can design the network with less number of nodes to have all gaps smaller than .

Figure 10

Probability of having no gaps larger than when .

4.3. Covered Part of the Path

The covered portion of the path, , is another important metric for path coverage in a WSN. Indeed, is a random variable whose cdf is approximated in Section 3.2. Figure 11 shows the cdf of , for and . As it can be seen, our path coverage analysis is more accurate for larger values of .

Figure 11

The cdf of the covered part of the path when .

The formulation for expected covered part, , is derived in Section 3.2. Figure 12 shows simulation and theoretical results for versus , when .

Figure 12

Expected size of the covered portion of the path versus .

Figure 13

Covered part of the path by a single sensor.

5. Conclusion

In this paper, we studied the path coverage of a random WSN when neither the area size nor the number of network nodes were infinite. Hence, the widely used Boolean model was no longer valid. Moreover, due to the randomness of the sensors placement over the area, network coverage was nondeterministic. Thus, a probabilistic solution was taken for determining the network coverage features. Our analysis considered the number of gaps, probability of full path coverage, probability of having all uncovered gaps smaller than a specific size, and the cdf of the covered length of the path. All these characteristics were found as a function of the number of sensors We also proposed a tight upper bound on required for full coverage. Through computer simulations, we verified the accuracy of our approach. Since our study was performed for finite , using our results on various features of path coverage, one can find the necessary number of sensors for a certain quality of coverage.

Appendix

In the following, we find the cdf of the intersection between the sensing area of the sensors and , called . First, we study the situation where sensors have a fixed sensing range and they are uniformly distributed over the ring. Then, we investigate the general case where sensors can have a random sensing range varying from to and have any symmetric distribution over the ring.

Let us first discuss the case where the sensors have a fixed sensing range. Figure 13 shows the ring-shaped network containing . As mentioned previously, the circumference of is 1, hence, the radius of is . It is also assumed that the ring width is and , where is the sensing radius of the sensors. Notice that in Figure 13 shows the distance of a sensor from the center of the ring. Since the sensors are uniformly distributed over the area, it can be easily shown that the cdf of , is as follows:

(A1)

We use to derive .In Figure 13, the intersection of the sensing area of an arbitrary sensor with is denoted by . By forming a triangle whose vertices are the center of , sensor location, and one of the points where the sensing circle of the sensor meets , one can write

(A2)

On the other hand, we have

(A3)

Replacing with in (A.2) results in

(A4)

Solving (A.4), we have

(A5)

Equivalently,

(A6)

Now having the cdf of and using the relation between and in (A.5) and (A.6), we will derive . To this end, one can state

(A7)

Thus,

(A8)

Replacing in (A.8) using (A.1), we obtain

(A9)

where

(A10)

Moreover, when and when . We use (A.9) for our exact analysis in order to characterize the path coverage features of the network. Notice that when is small, can be approximated as follows:

(A11)

In addition to the cdf of the arc length, we use the mean value of for our approximate analysis. Recall that for an arbitrary random variable distributed over ,

(A12)

where is the mean value of and is the cdf of . Using (A.12), can be found as follows:

(A13)

Notice that when .Now assume that both sensing range and sensor location are random and we like to find . Sensing range of the sensors, , varies over with probability density function (pdf) . Also, such that , because sensors located farther than from the path do not contribute in the path coverage. It is noteworthy that where

(A14)

This can simply be justified using (A.6).To find , we partition the problem to two separate cases. In the first case, sensing area of the sensor does not intersect with , that is, . This happens when or . If , this never happens and sensing area of the sensor always intersects with and consequently . If , we have

(A15)

To evaluate two terms in the right side of the above equation, we use the joint distribution of and . Notice that in the case where sensors sensing range is independent from their location, , where is the pdf of over . To evaluate and , we simply have to integrate from over the area where or . It can be shown that

(A16)

Notice that states that the pdf of , has a Dirac delta function at .When sensing area of a sensor intersects with path, . To find in this case, we first find . For this purpose, we apply Jacobian transformation to derive , the joint distribution of and , from . Using (A.6) and Jacobian transformation, one can show that

(A17)

where

(A18)

Having is found by integration over . To integrate over , the region of integration has to be determined carefully. For any arbitrary value of , there exist an infinite number of pairs satisfying (A.6); however, to guarantee an intersection between the sensing range of the sensor and , should fall within where

(A19)

In fact, and are the desired integral bounds. Thus,

(A20)

Consequently,

(A21)

The mean value of , used in our approximate analysis, is also derived as follows:

(A22)