Research Article  Open Access
A Planning and Optimization Framework for Ultra Dense Cellular Deployments
Abstract
To accommodate the everexpanding wireless data traffic volumes, mobile network operators are complementing their macrocellular networks by deploying lowpower base stations (or small cells) to offload traffic from congested macrocells and to reuse spectrum. To that end, Ultra Dense Network (UDN) deployments provide means to aggressively reuse spectrum, thus providing significant enhancements in terms of system capacity. However, these deployments entail several challenges, including the increased complexity in network planning and optimization. In this paper, we propose a versatile optimization framework for planning UDN deployments. The planning and optimization framework is underpinned by metrics that consider scalability in terms of number of users, cost of densification, and fairness. The proposed methodology is evaluated using a realworld UDN planning case. The numerical results expose a number of interesting insights, including the impact of different bandwidth allocation strategies and spatial service demand distribution on the performance of various network topologies. Specifically, we provide a performance comparison of the optimized UDN topologies versus random (unplanned), regular grid, and heuristically derived UDN topologies. This comparison further underlines the need for flexible network planning and optimization frameworks as different operator performance metrics of interest may require different radio access networks configurations.
1. Introduction
Mobile network operators face the continuous challenge of upgrading their networks which are rapidly expanding traffic volumes. This trend is mostly attributed to the increased adoption of smart devices (e.g., smartphones). Recent projections for global mobile traffic growth anticipate a tenfold increase in average monthly data consumption from the current 2–5 GB/month to 20–50 GB/month by 2020 [1]. Moreover, the average yearonyear subscriber growths of 5%–15% are expected to continue well into the next decade, notably with most of this demand from emerging markets [1, 2]. At the same time, user expectation on service quality also continues to increase, with highspeed connectivity becoming the baseline requirement for most users, regardless of their location or network load conditions.
To accommodate those projected traffic growths and meet user needs, mobile network operators are densifying their networks through heterogeneous deployment of lowpower base stations (BSs) or small cells (typically less than 10 W transmit power) to complement the existing highpower (20 W or higher) macrocells (umbrella coverage) [3–6]. Nowadays, small cell is a term to refer to compact lowpower BSs (e.g., microcells, picocells, or femtocells) and other macrocellular network extensions (e.g., relays, remote radio heads) that are deployed to enhance coverage and capacity in homes, enterprise environments, underserved areas, and other indoor and outdoor traffic hotspots [6]. This exponential data traffic growth is already setting an imperative requirement for Ultra Dense Networks (UDNs), identified as a key enabler for the 5th generation (5G) [1, 5, 7–9]. Indeed, 5G targets the operation in higher frequencies together with smaller cell sizes to achieve the envisioned extreme Mobile Broadband (xMBB) [1, 7]. To that end, UDNs are characterized by small cell deployments with intersite distances (ISD) of a few tens of meters for outdoor deployments (even shorter distances for indoor) and site density exceeding 100 sites/km^{2} in dense urban scenarios. This is in contrast to legacy 4th generation (4G) heterogeneous network deployments with typical site densities of less than 10 sites/km^{2} and ISD of a few hundred meters [1, 5, 7]. There is currently no commonly accepted definition on what network deployment constitutes a UDN. The definitions provided in different scientific literature have typically attempted to define UDNs in terms of cell density or the cell density relative to active user density (see [10] and references quoted therein). Other UDN definitions promulgated by industry include that of UDN being network with a small cell deployed outdoor on every lamp post or indoor with spacing of less than 10 m [1]. In this study, we adopt the pragmatic viewpoint from [1], whereby UDNs are considered to be an evolution from legacy dense networks and small cell network deployments with ISD of less than 100 m and are projected to become more prominent after year 2020.
However, UDNs are creating new and significant challenges for mobile network operators, with network planning and optimization notably becoming increasingly complex with denser network deployments [1, 7, 8]. In this context, planning refers to the process of determining the number, location, and configuration of base stations (e.g., small cells) to provide wireless access to users (and things) guaranteeing a certain targeted Quality of Service (QoS). In this process, dimensioning is the initial step used to solve the problem of estimating the required number of base stations needed to meet the capacity needs of a given service demand volume [11]. Thereafter, more precise network planning is carried out to evaluate cell site locations and initial cell parameters, and eventually optimization procedures in live networks are also used to continuously adjust cell’s parameters to further optimize both coverage and capacity. However, in practical scenarios, dimensioning and site positioning are nontrivial problems because the services are heterogeneous, that is, different QoS requirements, and the spatiotemporal distribution of the service demand is both nonuniform and dynamic. Furthermore, the challenges of site acquisition naturally scale with increased densification, thus obliging operators to consider leveraging base stations sites (mostly small cells) available in unplanned (suboptimal) locations [1, 4, 12]. Additionally, the ongoing evolution of radio access technologies, together with the new radio access concepts and paradigms expected for 5G, is blurring the traditional boundary between planning and optimization tasks. Indeed, as per discussion presented in [11], planning and optimization are iterative tasks that should be increasingly intertwined. These thoughts are echoed by the authors of [13], who also highlight the need for a rethink of planning and optimization in the context of dense heterogeneous networks, emphasizing that effective planning attains an even distribution of the load among cells, a goal that in the opinion of the authors of [11] (and corroborated by the authors of this paper) is a valid way to enhance system performance.
This paper addresses the aforementioned challenges by proposing an optimization framework for planning of UDN, which is suitable for realworld deployments. The corresponding research problem can be stated as follows.
Research Problem. Determine the set of network topologies with a certain number of access points (within an interval of interest, i.e., minimum and maximum node density) that is best compatible with a given spatial distribution of the service demand distribution (in statistical terms) and a certain performance metric.
Thus, the contribution of this paper, associated with the previous research problem, can be summarized as follows.
Main Contribution. A single and multiobjective optimization framework for planning of UDN deployments: The optimization allows obtaining network topologies which can be optimized for any arbitrary spatial traffic distribution (STD) (hereafter, the terms “spatial traffic distribution” and “spatial service demand distribution” are used interchangeably) and performance metrics, such as spectral efficiency or celledge performance.
Additional Contributions. In addition, several other minor contributions include:(1)The comparative analysis of several bandwidth allocation policies in the context of network planning.(2)A simple heuristic for planning of UDNs.
The numerical results from a realworld planning case (evaluated under a variety of conditions) reveal a number of interesting insights:(i)Bandwidth allocation strategies can facilitate the identification of optimized UDN topologies that may enable an operator to flexibly prioritize either system capacity or celledge performance.(ii)The results from the benchmarking clearly indicate that, in case of nonuniform STD, optimization is mandatory as the performance of regular and userdeployed (random) topologies is poor, while quasioptimal performances accompanied with significant gains can be attained through the use of heuristic planning and optimization.
The rest of the paper is organized as follows: the next section presents the system model. The performance metrics and proposed optimization formulations are introduced in Section 3. In Section 4, a background and description of the planning case study are presented together with the description of the spatial service demand distributions, benchmarks, and parameters and assumptions used in numerical evaluations. Section 5 provides a concise analysis of the numerical results. Finally, the concluding discussions and overview of potential research directions are provided in Section 6.
2. System Model
As indicated previously, the goal is to plan an ultra dense cellular network composed of lowpower BSs for a target service area . The service area is divided into small area elements or pixels (in this paper, the terms “area elements” and “pixels” will be used interchangeably) in which the average received power can be assumed to be constant.
In this study, the downlink of an Orthogonal Frequency Division Multiple Access (OFDMA) based cellular network with system bandwidth is considered. To carry out the planning, it is assumed that a set of candidate locations have been previously defined in the target service area. In each of these locations, a BS could be placed, and a maximum transmit power is assumed.
The radio propagation, that is, the network geometry, is captured by the matrix that indicates the average channel gain between each BS and area element. The vectors and , both , correspond to the transmit power of each BS in Reference Signals (RS) and data channels, respectively. The average RS received power can be calculated by means of the following expression:The operator denotes Hadamard (pointwise) operations. The binary vector indicates the allocation of a BS in the candidate locations, and hence, is referred to as “network topology” as it determines the number and location of BSs. Therefore, is the planning (optimization) variable. Hereafter, all the dependencies on are omitted for the sake of clarity. For instance, in (1). gives the average RS received power in the th pixel from the th BS.
Cell selection, the association of each pixel to a serving BS, is based on the average RS received power. Therefore, the th pixel (the th row in ) is served by cell ifThe coverage pattern associated with each network topology is represented by the binary coverage matrices and , both in . If the th area element is served by , then . is the binary complement of . It is assumed that each area element is either served by one cell or outofcoverage. It is considered that the th area element is outofcoverage if at least one of the following three conditions is not fulfilled: (i)The RS received power is larger than a minimum value: .(ii)The Signal to Interference plus Noise Ratio (SINR) is larger than a threshold: .(iii)The average channel gain between the area element and its serving BS is larger than .
The outage associated with a network topology is captured by the vector . If the th area element is outofcoverage, then , and 0 otherwise.
A certain knowledge of the spatial distribution of the service demand is assumed. In practice, this is known by operators in statistical terms [14]. This information is stored in the vector . can be regarded as a Probability Density Function (PDF) in two dimensions, and hence, it indicates the probability, in the event of a new user, that the th pixel has the user on it. Thus, .
In this work, full load is assumed to model the intercell interference, which is a reasonable assumption for planning purposes. Other models, such as loadcoupling [15], can be easily incorporated in the model, if needed. Thus, the vector representing the average SINR at each area element is given bywhere is the noise power. The operators and denote Hadamard (pointwise) operations. It is customary to define link performance as a nondecreasing function of the SINR. In this work, Shannon’s bound is considered, and hence, the resulting spectral efficiency is stored in the vector , and its elements are calculated as follows:In (4), the idea is to discard the contribution of the pixels that are outofcoverage (by means of ), thus penalizing network topologies with significant coverage holes in the optimization procedure.
The list of symbols is provided in Basic Notation in Notation for convenience.
3. Performance Metrics and Optimization
3.1. Performance Metrics for Radio Access Network Planning
Generally speaking, planning is about determining the number and location of BSs in the service area. Evidently, the network deployment should be done such that the maximum benefit is obtained; that is, network capacity is maximized (more users) with minimal cost (less infrastructure deployment), while guaranteeing a certain level of coverage, QoS, and fairness. In order to address this problem by means of optimization, several metrics (and constraints) are required. In this work, the following objectives are considered.
(i) Number of BSs (f_{1}). In principle, the deployment should be done with the minimum possible number of BSs to minimize both the Capital Expenditure (CAPEX) and the energy consumption that is part of the Operational Expenditure (OPEX).
(ii) Network Capacity (f_{2}). This metric captures the average aggregate rate the network is able to deliver. Thus, represents a systemoriented performance indicator.
(iii) CellEdge Performance (f_{3}). This metric captures the performance in the weakest zones of the service area. Thus, is a useroriented performance indicator and promotes fairness.
The definition of the previous metrics is given next. The number of BSs () in a network topology is simply the number of ‘1’s in the corresponding , and hence,
From a planning point of view, it is important to consider the spatial distribution of the service demand. In other words, the planning should favor network topologies that provides more capacity to the zones of the service area where the traffic is more likely to appear. Given that this information is contained in the vector , it can be used to weight the different pixels according to their importance; that is, pixels with more traffic are more important. Thus, the weighted spectral efficiency vector is defined as follows: . Note that, indeed, the scalar represents the expected spectral efficiency at area element level because is a PDF.
In cellular networks, a very important aspect is the frequency reuse; that is, the system bandwidth can be reutilized at each cell. This is the most distinctive aspect of cellular networks that allows these systems to provide radio access to a large amount of users (and things). The way in which the bandwidth is allocated to the users largely determines the resulting system capacity and/or users’ satisfaction. In this sense, celledge performance [16] is a wellknown, yet important, problem in OFDMAbased cellular networks that can affect (negatively) user’s experience. Broadly speaking, users at celledges are relatively more expensive in terms of radio resources, as their SINR is typically very low. In this work, this aspect is considered from the planning point of view, and consequently, two different bandwidth allocation strategies are considered and integrated in the performance metrics.
(1) Uniform Bandwidth Allocation (UBA). The objective is to evaluate the resulting aggregate capacity assuming that the bandwidth of each cell is equally distributed over its coverage area (pixels).
(2) Proportional Bandwidth Allocation (PBA). The objective is to evaluate the resulting aggregate capacity assuming that the bandwidth of each cell is distributed over its coverage area (pixels) proportionally to the SINR or, equivalently, the spectral efficiency of the pixels.
The vectors and , both in , indicate the bandwidth that would be allocated to each area element under the uniform and proportional bandwidth allocation, respectively. They are defined as follows:where the vector contains the inverse of the number of pixels associated with each BS. The proportional allocation is as follows:Equation (7) divides the bandwidth of each cell proportionally to the spectral efficiency of the area elements in the cell.
Thus, the network capacity metric is defined as follows: Hereafter, superscripts “” and “” are used to indicate UBA and PBA, respectively, as follows: and .
Celledge performance is defined, for planning purposes herein, as the aggregate rate of the worst of the service area. Given the vectors and , where the sorting is in ascending order, then the metric representing the celledge performance would be given by where . Since and can be evaluated for both UBA and PBA, then this study is able to utilize four possible objective functions for comparison purposes, namely, , , , and .
3.2. Optimization Problem Formulation
In this work, two different optimization formulations are considered. They can be used depending on the network planning strategy of the operator. On the one hand, if the network operator’s target is to maximize the network aggregate capacity (), a multiobjective problem is proposed as this metric is in conflict with ; that is, generally speaking, the denser the network, the higher the capacity due to the more aggressive frequency reuse. On the other hand, if the operator’s target is to provide a more homogeneous coverage, that is, less variability at pixel level, a singleobjective problem is proposed with as objective function, and the required number of cells as an input.
3.2.1. Multiobjective Optimization
Multiobjective optimization [17] can be used when multiple conflicting objectives need to be simultaneously optimized (a brief introduction to multiobjective and evolutionary optimization is presented in Appendix A). This is the case of and in the planning framework presented herein. Thus, in order to obtain network topologies featuring the best tradeoff between the number of BSs () and network capacity (), the following multiobjective optimization problem is proposed:
In problem (10a), (10b), (10c), and (10d), constraint (10b) guarantees that a minimum fraction () of the area elements has coverage. Constraint (10c) defines the search space, that is, the domain of the variable . In practice, and due to the nature of the environments in which UDNs are envisioned to be deployed, network operators usually have an estimate of the number of BSs that is required/feasible, and hence, the optimization can be further localized. This is accomplished by means of constraint (10d), where these limits are set.
Problem (10a), (10b), (10c), and (10d) is a combinatorial problem belonging to the class NPcomplete. The search space defined by the optimization variable (the total number of network topologies) is a set of size , where , as indicated, is the number of candidate locations. Even for a small set of candidate locations, say , the number of network topologies would be larger than , which makes it infeasible to compare all possible topologies by means of timeconsuming and computationally heavy system level simulations. For this reason, the proposed planning approach is a convenient strategy. The objective space (or image) is defined by the possible values of the objective functions. Due to the mathematical structure of and , the objective space is highly nonlinear, nonconvex, and full of discontinuities and local optima [18]. Thus, a multiobjective evolutionary algorithm (MOEA) [19], the Nondominated Sorting Genetic Algorithm II (NSGAII) [20], is used to address (10a), (10b), (10c), and (10d). A brief description is provided in Appendix A.
3.2.2. SingleObjective Optimization
Singleobjective optimization is proposed if planning needs to be carried out following a maxmin approach, such as the maximization of aggregate rate in the area elements with weak coverage. Thus, the problem of maximizing celledge performance (), for planning, can be written as follows:
Problem (11a), (11b), (11c), and (11d) and its constraints are similar to (10a), (10b), (10c), and (10d), except that it contains only one objective function. Constraints (11b) and (11c) are equal to (10b) and (10c), respectively. Constraint (11d) indicates that only solutions with BSs are accepted. This is so because, in general, the metric is proportional to due to the reduced level of intercell interference in network topologies with less BSs. However, the fact that the planning framework presented herein can consider both types of optimization (i.e., singleobjective and multiobjective) is an indication of the versatility of the framework. Indeed, more than two objectives can also be considered, but at expense of significant increase in complexity. Moreover, in the opinion of the authors, the simultaneous consideration of more than two metrics would add complexity in the meaningful interpretation of the results.
4. Planning Case Study
4.1. Deployment Scenario
The network densification as planned by operators is both difficult and highly constrained in certain scenarios. This includes the fast expanding highdensity urban and periurban settlements in emerging market areas. Indeed, 90% of the urban population growth by 2050 is expected to be concentrated in Asia and Africa. These settlements already have populations densities typically in the range of 40,000–200,000 people/km^{2} [21]. Mobile broadband networks continue to be the primary means for wireless connectivity in these densely populated areas [2, 22], which makes them a highly compelling target for the deployment of UDN. Unfortunately, some challenges related to UDN deployment are further exacerbated in these areas due to the limited availability of legacy infrastructure for small cell backhaul, energy scarcity, difficulties in site acquisition, need for securing network assets at sites, and relatively low Average Revenue Per User (ARPU) compared to more developed economies [22].
One of the interesting approaches is to leverage thirdparty nonoperator entities, such as individual end users, households, microenterprises, and public venue owners, to deploy sharedaccess small cells that will provide service as an integral part of the operator’s network. An example is the neighborhood small cell concept by Qualcomm promoting the use of privately deployed residential small cells as sharedaccess points [12]. A key distinction between thirdparty deployments and operatorled deployments is that, in the former case, the small cells deployments are unplanned; that is, the location in which small cells are deployed is not originally defined by the operator’s network planning procedures. However, although small cells are deployed autonomously by thirdparties, the operator retains remote management via core network and the use of SelfOrganizing Networks (SON) [23, 24].
Therefore, to contextualize the proposed UDN planning and optimization framework in a realistic setting, we consider a case study for UDN in a highdensity urban settlement. To that end, we use the Hanna Nassif ward in Dar es Salaam, Tanzania, as the planning study case. Hanna Nassif has an estimated population density of 40000 people/km^{2}. The approximately 1 km^{2} Hanna Nassif area includes around 3000 buildings (mostly 3–6 m tall) and is located on a terrain with a topographical difference of 19 m. A threedimensional (3D) representation of this scenario is shown in Figure 1. We assume that all candidate locations (indicated in whiteblue dots) are outdoor at rooftop level. Rooftop deployed sharedaccess small cells provide improved outdoor coverage compared to indoor deployed small cells and enable lineofsight (LOS) or near LOS (nLOS) conditions for the implementation of highcapacity wireless backhauling [25]. Moreover, the rooftop is also a convenient location for offgrid operation of the small cells through energy harvesting from ambient renewable energy sources (solar, wind, etc.) [26].
4.2. Parameters and Assumptions
The radio coverage estimations are based on realistic 3D building vectors and topographical data (see Figure 1, for the simulation area) and are evaluated using the deterministic dominant path model implemented in the WinProp propagation modeling tool [27]. The simulation parameters and assumptions follow the IMTAdvanced guidelines [28] and are listed in Table 1. The building penetration losses are modeled explicitly to account for the outdoortoindoor propagation and it is assumed that all buildings outer walls are based on one material (10 cm brick) and inbuilding losses (e.g., due to internal walls, doors) are approximated by an exponential decay model. Furthermore, as noted previously, the case study area includes a large number of densely packed small houses or buildings of relatively low height built on a land with significant topographical differences over short distances. The deployment of rooftop smalls in this environment creates significant multipath propagation, waveguiding effects, and outdoortoindoor propagation that can be suitably captured by a raytracing tool like WinProp (see, e.g., [29]). Indeed, it is noted that ray tracing provides significant accuracy improvements compared to classical models but with reduced generality in terms of scenario selection [30]. The precision (pixel resolution) and inclusion of timevarying nonstationary objects (e.g., cars, pedestrians) may further influence accuracy and computational effort of different raytracing computations. To that end, the pixel resolution in our study is , which provides the right tradeoff between modeling accuracy and feasible computational time. However, in WinProp environment we do not include the scattering effects of nonstationary objects for the sake of simplicity in capturing the scenario. Moreover, the usage of same raytracing path loss results for all topologies considered in the study ensures that the absolute accuracy of the ray tracing does not influence the performance comparisons between the topologies. The parameters for the calibration of the NSGAII algorithm are also shown in Table 1. Calibration parameters and complexity aspects are provided in Section 5.4.

4.3. Spatial Traffic Distributions
Two different spatial traffic distributions have been considered for the numerical evaluations presented in the next section: uniform and nonuniform. Uniform STD implies that the service demand is uniformly distributed in the coverage area. Nonuniform STD implies that the traffic is more likely to appear in certain areas (hotspots), which is representative of how traffic is commonly distributed in practice. The system model and optimization framework presented herein are able to consider any arbitrary STD by means of the vector , as it is explained in Section 3.1. Figure 2 shows a representation of the nonuniform STD used in numerical evaluations. A 2dimensional representation (map) is provided in Figure 2(a), where red areas correspond to hotspots. Figure 2(b) shows the CDF of the probability at pixel level in order to provide a perspective of the level of irregularity of the nonuniform STD. The uniform STD is also indicated by the vertical dashed line. This input, as it will be seen shortly, has a profound impact on the optimization process and leads to significantly different network topologies.
(a) Nonuniform service demand distribution map ()
(b) CDF of
4.4. Benchmarks
In order to clarify the merit of the proposed optimization framework, several benchmarks have been considered.
(1) Random Deployments. As explained earlier, openaccess outdoor small cells (randomly deployed by users) would provide a sustainable (cost efficient) approach to network densification. Thus, random topologies (located at candidate locations) are considered.
(2) Regular Deployments. Regularly deployed access points, that is, geometrically regular cells, are the natural choice for uniform STD and are often considered as a baseline. Two regular topologies ( and ) are considered and are shown in Figure 3.

(3) Heuristic Network Planning. A basic heuristic network planning for small cells is also considered. The corresponding pseudocode (Algorithm 1) and its explanation are provided in Appendix B. The resulting network topologies ( and ) from the heuristic algorithm are also shown in Figure 3 (both for uniform and nonuniform STD).
5. Numerical Results
5.1. Optimization for Planning
The solution of multiobjective optimization problems, such as (10a), (10b), (10c), and (10d), is a set of network topologies, all them featuring Pareto efficiency with respect to the objective functions, in this case and . The image of this set is called Pareto Front (see Appendix A), and it is shown in Figure 4 for both uniform and nonuniform spatial traffic distribution, using different target ranges for . The  and axes indicate the values of and , respectively. Both uniform and proportional bandwidth allocation are considered in each case. The conflicting nature of and is evident; that is, the larger the number of BSs, the larger the network capacity, which is expected due to the higher frequency reuse. The result shows that uniform bandwidth allocation provides significantly better network capacity compared with proportional bandwidth allocation; see (6) and (7). In the figure, the difference (≈60%) in terms of between optimized topologies with 180 BSs (for the case of uniform STD) is shown as a reference. This result, although expected, is important because it is wellaligned with system level simulations [31], and more precisely with different scheduling policies (Round Robin or Proportional Fair), where it is wellknown that fairness is always traded by spectral efficiency. Thus, the performance metrics proposed herein for planning capture the same behavior.
Analogously, problem (11a), (11b), (11c), and (11d) can be solved for both uniform and proportional bandwidth allocation; however, only one solution can be obtained in each case (since singleobjective optimization can be regarded as a particular case of multiobjective optimization, (11a), (11b), (11c), and (11d) was also solved using NSGAII). In order to compare optimized topologies, solutions of problem (10a), (10b), (10c), and (10d) with have been selected, as this value was also used in (11a), (11b), (11c), and (11d), where . Hereafter, the notation used to describe the four optimized network topologies is indicated in Notation Used to Refer to the Optimized Topologies in Notation.
The comparison is shown in Figure 5. The figure shows the (normalized) performance of each network topology and objective function, that is, , , , and . As it can be seen, each of these topologies maximizes one metric, and hence, they are all Pareto efficient. Both uniform and nonuniform spatial traffic distribution are considered. It is worth noting that, in case of nonuniform spatial traffic distribution, the performance of each optimized topology with respect to the metrics for which it is not optimized becomes notably poor. Since uniform STD is not common in practice, the result confirms that optimization is highly needed.
A pictorial representation of these optimized topologies is shown in Figure 6. It is clear that they are different, as each of them is optimized for a certain metric. Therefore, it can be concluded that, in general, each performance indicator (representing different operator interests) requires an independent optimization. The spatial correlation between the density of BSs of the optimized topologies (for the case of nonuniform STD) and the nonuniform pattern considered for numerical evaluations should be noted (see Figure 2(a)).
It is worth noting that, in the study case considered herein ( and ), there are more than possible network topologies (the number of combinations of candidate locations with BSs).
In order to justify and explain the previous differences in performance, Figure 7 provides an alternative point of view. The results correspond to the case of uniform STD, but the same analysis applies for the nonuniform case. The figure shows the statistics at pixel level of the resulting spectral efficiency (Figure 7(a)) and average rate (Figure 7(b)), both normalized. As it is clear, there is not major difference in terms of the statistic of the resulting spectral efficiency; only presents a slightly better distribution. This is somehow expected considering their dense nature and the relatively homogeneous propagation characteristics of the study case; see Figure 1. However, the impact of the different bandwidth allocation policies is noticeable, as it can be seen in Figure 7(b). The topologies and clearly enhance not only the fairness among area elements, but also celledge performance (despite that, for instance, is a metric measuring network capacity). The results allow us to conclude that the bandwidth allocation policies, introduced herein in the context of planning, facilitate the identification of network topologies to enhance either system capacity or celledge performance (better fairness), thus, demonstrating the effectiveness and flexibility of the proposed planning framework. So far, the optimized topologies have been compared among themselves; next, benchmarking is presented in order to provide a quantitative perspective of the merit of the proposed planning framework.
(a) Spectral efficiency
(b) Achievable rate
5.2. Benchmarking: Random Deployments
Optimized topologies are compared with random deployments. In order to make the comparison fair, 1000 random deployments with BSs (located in candidate locations) have been considered and evaluated in terms of , , , and . As additional measure of fairness, Jain’s Index (Jain’s index is usually used to determine whether users or applications are receiving a fair share of system resources) [32] has also been considered taking the resulting rate at pixel level as input. Thus, if the resulting rate in the th pixel is , Jain’s index is given bywhere , being and the worst and ideal case, respectively. It is worth mentioning that the same set of random topologies is used in each case.
The comparisons for uniform and nonuniform STD are shown in Figures 8 and 9, respectively. Figures 8(a), 8(b), 9(a), and 9(b) show the empirical Cumulative Density Function (CDF) of and for the random deployments (thinsolid lines). Figures 8(c), 8(d), 9(c), and 9(d) show the CDF of Jain’s index for the random deployments (thinsolid lines). In all the cases, the performance of the corresponding optimized topology is also indicated (vertical thickdashed lines). Since the analysis of Figures 8 and 9 is similar; the focus is placed initially on Figure 8. Next, relevant particularities of Figure 9 are discussed.
(a) CDF of
(b) CDF of
(c) CDF of
(d) CDF of
(e) Legend
(a) CDF of
(b) CDF of
(c) CDF of
(d) CDF of
(e) Legend
From Figure 8(a), it can be seen that the gain achieved by the optimized topologies with respect to the random deployments is, on average, 14% and 19%, for uniform and proportional bandwidth allocation, respectively. In practice, these gains can be even larger, as in random topologies (e.g., deployed by customers), BSs can potentially be anywhere, without respecting a minimum distance as it has been set for the candidate locations used herein. The difference (≈60%) between the optimized topologies is also shown, indicated as and . The result makes the significant benefit that can be achieved by means of proper UDN planning evident. The same behavior has been verified for celledge performance (), although with greater gains, as it can be seen in Figure 8(b). Thus, it can be concluded that if the network operator target is maximizing the celledge performance, the need for planning is mandatory. Gains with respect to random deployments are, on average, and , while the celledge performance gain by using the proportional bandwidth allocation is (between optimized topologies). In this manner, and for the case of study, it would be possible to trade a gain of in terms of by a reduction of in terms of , or vice versa.
Figures 8(c) and 8(d) validate that, by means of the proportional bandwidth allocation policy, the overall system fairness (measured in terms of Jain’s index) is significantly enhanced, for both and . The high gains in terms of Jain’s index ( and ) and the fact that, in most of the cases, the optimized topologies provide better fairness than random deployments although Jain’s index is not the objective function under optimization should be remarked.
As mentioned, the analysis of Figure 9 is similar as well as the obtained results, qualitatively speaking. However, optimization gains when nonuniform service demand distribution is considered are notably larger. This is expected because uniform STD is indeed the worst case from the optimization point of view, which is wellaligned with the information theoretical intuition that uniform distribution maximizes the entropy. Thus, both situations are presented herein, making (even in the worst case) the merit and potential benefit of the proposed optimization framework clear.
5.3. Benchmarking: Regular Deployments and Heuristic Planning
In order to provide another quantitative perspective of the merit of the proposed scheme, a comparative assessment with regular deployments and network topologies obtained through basic heuristic planning, in terms of network capacity (), have also been carried out. The results are shown in Figure 10. Figures 10(a) and 10(b) correspond to uniform STD and Figures 10(c) and 10(d) correspond to nonuniform STD. For each case, both bandwidth allocation policies (uniform and proportional) have been taken into account. Similarly to the benchmarking with random topologies, larger gains (between and ) are consistently obtained when nonuniform STD is considered (the case of practical interest). This is expected as regular deployments are not likely to be effective for nonuniform STDs. However, the gains even in the worst case (uniform STD) range from to . It is worth mentioning the dense nature of the deployment considered for this comparison: BSs in the target area (resulting in ≈514 access points per km^{2}). Thus, the proposed optimization formulation succeeds in finding neartooptimal and Pareto efficient network topologies under different bandwidth allocation strategies and spatial traffic distribution conditions. The effectiveness of the optimization comes from the fact that the objective functions take into account the spatial traffic distribution (prioritizing the topologies providing better capacity where the traffic is more likely to appear) and the radio propagation characteristics obtained from every single candidate location (favoring better wireless links with serving access points and more isolation among interfering cells). As mentioned earlier, other performance metrics (e.g., Jain’s index) can also be considered depending on operator’s needs/interest. Finally, it is worth mentioning the notorious performance of the heuristic planning (Algorithm 1), which can also be used if optimization tools are not available; and in that sense, it can also be considered as a part of the UDN planning framework presented herein.
(a) Uniform bandwidth allocation:
(b) Proportional bandwidth allocation:
(c) Uniform bandwidth allocation:
(d) Proportional bandwidth allocation:
(e) Legend
All in all, the previous results show (and confirm) that effective planning is strongly recommended when it comes to UDN. The evaluation and benchmarking presented provide evidence of the merit, flexibility, and effectiveness of the proposed framework under several conditions.
5.4. Complexity and Calibration Aspects
To close this section, complexity and calibration aspects are discussed. According to [20], the complexity of NSGAII is , where and correspond to the population size and the number of objective functions, respectively. In our case, and can be set depending on the scale of the problem. However, there is a consensus about the size of the population when using genetic algorithms, such as NSGAII, and it is considered that during calibration populations of 20 up to 100 individuals can be used. Values greater than 100 hardly achieve significant gains and the same global convergence is obtained [33]. In evolutionary algorithms, a termination criterion is usually defined/needed. One metric used to measure the level of convergence is the hypervolume indicator [19]. It reflects the size of volume dominated by the estimated Pareto Front. In this work, the search is terminated if the improvement in the hypervolume is smaller than a threshold (0.001%) after a certain number of generations (in this study, 20). Finally, crossover and mutation probabilities are set to and (one mutation per solution, on average), respectively, as indicated in Table 1.
6. Conclusion and Future Work
Network planning is the process of determining the amount and location of BSs that are required to provide wireless access to a given service area. In this article, a novel planning framework for UDN has been presented to address the general planning problem: that of finding network topologies which are compatible with a certain spatial traffic distribution while maximizing the benefit in terms of a given performance metric. To that end, a flexible multiobjective framework (system model and optimization formulation) has been presented. One interesting feature of the proposed scheme is the use of different bandwidth allocation strategies, which have been shown to have a profound impact on the ability to identify network topologies delivering the best performances in terms of the metrics of interest (network capacity, fairness, or celledge performance).
In addition, by benchmarking with regular and random deployments as well as network topologies obtained though simple heuristic planning, the notable benefit that can be achieved by means of the proposed optimization framework was quantified. Based on these results, the need for methodological planning means (such as the one provided in this article) becomes strongly recommended, specially in scenarios where the service demand is not uniformly distributed (the case commonly found in practice).
In the light of the results, it can be concluded that planning is still a quite valid tool for mobile operators, both nowadays and in future 5G systems, where small cells are of utmost importance. The framework studied herein is rich and admits several future research directions. Considering additional objective functions is definitely a study item, as well as other optimization formulations that could be specific for certain use cases, such as Downlink Uplink Decoupling. In addition, upgrading fixed/existing deployments is of great practical interest, as well as adaptation of coverage patterns, and power optimization. Finally, planning studies for UDN operating at higher frequency bands is also on our roadmap.
Appendix
A. Multiobjective and Evolutionary Optimization
Multiobjective optimization is the discipline that focuses on the resolution of problems involving simultaneous optimization of several conflicting objectives. The target is to find a subset of good solutions from a set (the domain of an optimization variable ) according to a set of criteria , with cardinality greater than one. Since the objectives are in conflict, improving one of them implies worsening another. Consequently, in the context of multiobjective optimization, there is no one single optimal solution but an optimum set . A central concept in multiobjective optimization is Pareto efficiency. A solution has Pareto efficiency if and only if there does not exist a solution , such that dominates . A solution is preferred to (dominates) another solution , (), if is better than in at least one criterion and not worse than any of the remaining ones. The set , composed of all the solutions featuring Pareto efficiency, is called optimal nondominated set, and its image is known as the Optimal Pareto Front (OPF). Figure 11 illustrates this idea, where it can be seen that the nondominated solutions (in blue) are the ones not dominated by any other, and hence, they provide a tradeoff that cannot be improved.
In multiobjective optimization, it is unusual to obtain the OPF due to problem complexity; instead, a nearoptimal or estimated Pareto Front (PF) is found. Interested readers are referred to [17] for an indepth discussion. In general, solving multiobjective (combinatorial) problems, such as (10a), (10b), (10c), and (10d), is very difficult [17]. Indeed, optimal solutions for problems belonging to the class NPcomplete cannot be found in polynomial time. For this reason, heuristicbased algorithms are commonly used to address this type of problems, but unfortunately, heuristic solutions are typically problemspecific, and hence, their use is limited. Thus, a more general class of solutions, the socalled “metaheuristics,” have become very popular and an active research field [18]. Metaheuristics can be used to solve very general kind of multiobjective optimization problems. These methods allow finding good solutions by efficiently exploring the search space and operating efficiently with multiple criteria and a large number of design variables. In addition, no requirements on the mathematical structure of the objective functions are required (e.g., convexity or continuity).
Multiobjective evolutionary algorithms (MOEAs) [19] are metaheuristics that fulfill the previous goals. MOEAs are populationbased metaheuristics that simulate the process of natural evolution and they are convenient due to their generalpurpose nature. One MOEA, the Nondominated Sorting Genetic Algorithm II (NSGAII) [20], is wellrecognized as a reference in the field of evolutionary optimization as it has desirable features, such as elitism (the ability to preserve good solutions), and mechanisms to flexibly improve convergence and distribution. Interested readers are referred to [17, 19], and the references therein, for an indepth treatment of multiobjective and evolutionary optimization.
B. Heuristic Network Planning for Small Cell Deployments
Heuristic solutions have been considered for solving complex optimization problems in many fields, including engineering. Since planning and topology optimization for UDN are NPcomplete problems, heuristics are naturally an option. In this appendix, a simple heuristic for planning of UDN is presented. The idea is adapted to planning from the Minimum Distance Algorithm used in Cell SwitchOff [34]. The pseudocode is shown in Algorithm 1.
The basic idea is to start identifying the best single access point (a network with only one cell) for a given metric (), Line . This would be the network with cell. Then the loop (Lines to ) iteratively find the best new candidate location to create the network with cells (and the previously selected locations). Thus, the loop in the th level has to search for candidate locations. The loop is repeated until locations are selected. In Line , corresponds to the Hamming distance between two different network topologies, which guarantees that the previously found cells are fixed in the th iteration. In Line , corresponds to the performance of the network topology in terms of the metric .
Notation
Basic Notation:  System bandwidth 
:  Number of candidate locations 
:  Number of area elements 
:  Channel gain matrix 
:  Target received power 
:  Minimum channel gain 
:  Spatial demand distribution 
:  Spectral efficiency vector 
:  Average SINR vector 
:  Inverse of cell’s size 
:  Maximum transmit power per cell 
:  Set of area elements in the target area 
:  Received power matrix 
:  Power vectors: pilots and data channels 
:  Network topology 
:  Coverage of the th BS 
:  Coverage matrices 
:  Minimum SINR 
:  Coverage vector 
:  Coverage threshold 
:  Optimal topology with respect to using uniform bandwidth allocation (UBA) 
:  Optimal topology with respect to using proportional bandwidth allocation (PBA) 
:  Optimal topology with respect to using uniform bandwidth allocation (UBA) 
:  Optimal topology with respect to using proportional bandwidth allocation (PBA). 
Competing Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
Hanna Nassif GIS data was kindly provided by Professor R. Sliuzas of ITCFaculty of GeoInformation Science & Earth Observation, University of Twente. This work is partially supported by the Academy of Finland under Grants 287249 and 284634. Furthermore, this work received partial support from the National Bureau of Science, Technology and Innovation of Panama (SENACYT) through the project RAPIDO (ITE15021).
References
 Nokia, Ultra dense network (UDN) white paper. Nokia Solutions and Neworks Oy, June 2016, http://resources.alcatellucent.com/asset/200295.
 Ericsson, Ericsson mobility report: on the pulse of the networked society, June 2016, https://www.ericsson.com/res/docs/2016/ericssonmobilityreport2016.pdf.
 J. Zander and P. Mähönen, “Riding the data tsunami in the cloud: myths and challenges in future wireless access,” IEEE Communications Magazine, vol. 51, no. 3, pp. 145–151, 2013. View at: Publisher Site  Google Scholar
 I. Hwang, B. Song, and S. Soliman, “A holistic view on hyperdense heterogeneous and small cell networks,” IEEE Communications Magazine, vol. 51, no. 6, pp. 20–27, 2013. View at: Publisher Site  Google Scholar
 D. LópezPérez, M. Ding, H. Claussen, and A. H. Jafari, “Towards 1 Gbps/UE in cellular systems: understanding ultradense small cell deployments,” IEEE Communications Surveys and Tutorials, vol. 17, no. 4, pp. 2078–2101, 2015. View at: Publisher Site  Google Scholar
 Small Cell Forum (SCF), “Small cells—what is the big idea? Small Cell Forum,” February, 2012 http://scf.io/. View at: Google Scholar
 X. Ge, S. Tu, G. Mao, C.X. Wang, and T. Han, “5G ultradense cellular networks,” IEEE Wireless Communications, vol. 23, no. 1, pp. 72–79, 2016. View at: Publisher Site  Google Scholar
 H. Peng, Y. Xiao, Y.N. Ruyue, and Y. Yifei, “Ultra dense network: challenges, enabling technologies and new trends,” China Communications, vol. 13, no. 2, Article ID 7405723, pp. 30–40, 2016. View at: Publisher Site  Google Scholar
 W. Yu, H. Xu, H. Zhang, D. Griffith, and N. Golmie, “Ultradense networks: survey of state of the art and future directions,” in Proceedings of the 25th International Conference on Computer Communication and Networks (ICCCN '16), pp. 1–10, Waikoloa, Hawaii, USA, August 2016. View at: Publisher Site  Google Scholar
 M. Kamel, W. Hamouda, and A. Youssef, “Ultradense networks: a survey,” IEEE Communications Surveys & Tutorials, vol. 18, no. 4, pp. 2522–2545, 2016. View at: Publisher Site  Google Scholar
 M. Jaber, Z. Dawy, N. Akl, and E. Yaacoub, “Tutorial on LTE/LTEA cellular network dimensioning using iterative statistical analysis,” IEEE Communications Surveys and Tutorials, vol. 18, no. 2, pp. 1355–1383, 2016. View at: Publisher Site  Google Scholar
 Qualcomm Research, “Neighborhood small cells for hyperdense deployments: taking hetnets to the next level,” Tech. Rep., Qualcomm, San Diego, Calif, USA, 2013, https://www.qualcomm.com/. View at: Google Scholar
 S. Wang and C. Ran, “Rethinking cellular network planning and optimization,” IEEE Wireless Communications, vol. 23, no. 2, pp. 118–125, 2016. View at: Publisher Site  Google Scholar
 X. Zhou, Z. Zhao, R. Li, Y. Zhou, and H. Zhang, “The predictability of cellular networks traffic,” in Proceedings of the International Symposium on Communications and Information Technologies (ISCIT '12), pp. 973–978, October 2012. View at: Publisher Site  Google Scholar
 I. Siomina and D. Yuan, “Analysis of cell load coupling for LTE network planning and optimization,” IEEE Transactions on Wireless Communications, vol. 11, no. 6, pp. 2287–2297, 2012. View at: Publisher Site  Google Scholar
 G. D. Gonzalez, M. GarciaLozano, S. Ruiz Boque, and D. S. Lee, “Optimization of soft frequency reuse for irregular LTE macrocellular networks,” IEEE Transactions on Wireless Communications, vol. 12, no. 5, pp. 2410–2423, 2013. View at: Publisher Site  Google Scholar
 Y. Sawaragi, H. Nakayama, and T. Tanino, Theory of Multiobjective Optimization, vol. 176 of Mathematics in Science and Engineering, Academic Press, Orlando, Fla, USA, 1st edition, 1985. View at: MathSciNet
 T. Weise, Global Optimization Algorithms&Theory and Application, 2nd edition, 2009, http://www.itweise.de/.
 C. A. Coello Coello, G. B. Lamont, and D. A. Van Veldhuizen, Evolutionary algorithms for solving multiobjective problems, Genetic and Evolutionary Computation Series, Springer, NY, USA, Second edition, 2007. View at: MathSciNet
 K. Deb, A. Pratap, S. Agarwal, and T. Meyarivan, “A fast and elitist multiobjective genetic algorithm: NSGAII,” IEEE Transactions on Evolutionary Computation, vol. 6, no. 2, pp. 182–197, 2002. View at: Publisher Site  Google Scholar
 UN, “World urbanization prospects: the 2014 revision, highlights,” UN Report ST/ESA/SER.A/352, United Nations, Department of Economic and Social Affairs, Population Division, New York, NY, USA, 2014, https://esa.un.org/unpd/wup/Publications/Files/WUP2014Highlights.pdf. View at: Google Scholar
 GSMA, “The mobile economy Africa 2016,” 2016, http://www.gsma.com/mobileeconomy/africa/. View at: Google Scholar
 H. Klessig, D. Öhmann, A. I. Reppas et al., “From immune cells to selforganizing ultradense small cell networks,” IEEE Journal on Selected Areas in Communications, vol. 34, no. 4, pp. 800–811, 2016. View at: Publisher Site  Google Scholar
 3GPP, “Evolved Universal Terrestrial Radio Access Network (EUTRAN); selfconfiguring and selfoptimizing network (SON) use cases and solutions,” 3GPP Standard TR 36.902, 2011, http://www.3gpp.org/dynareport/36902.htm. View at: Google Scholar
 P. Amin, N. S. Kibret, E. Mutafungwa, B. B. Haile, J. Hamalainen, and J. K. Nurminen, “Performance study for offgrid selfbackhauled small cells in dense informal settlements,” in Proceedings of the 25th IEEE Annual International Symposium on Personal, Indoor, and Mobile Radio Communication (IEEE PIMRC '14), pp. 1652–1657, September 2014. View at: Publisher Site  Google Scholar
 Y. Mao, Y. Luo, J. Zhang, and K. B. Letaief, “Energy harvesting small cell networks: Feasibility, deployment, and operation,” IEEE Communications Magazine, vol. 53, no. 6, pp. 94–101, 2015. View at: Publisher Site  Google Scholar
 Altair, Winprop overview, http://www.altairhyperworks.com/product/FEKO/WinProp.
 ITUR, “Guidelines for evaluation of radio interface technologies for IMTAdvanced,” ITU Std. ITUR, M.2135, 2008. View at: Google Scholar
 L. Simic, J. Riihijarvi, and P. Mahonen, “Can statistical propagation models be saved by real 3D city data?: A Regionalized Study of Radio Coverage in New York City,” in Proceedings of the IEEE International Symposium on Dynamic Spectrum Access Networks (DySPAN '15), pp. 285–288, October 2015. View at: Publisher Site  Google Scholar
 P. Ferrand, M. Amara, S. Valentin, and M. Guillaud, “Trends and challenges in wireless channel modeling for evolving radio access,” IEEE Communications Magazine, vol. 54, no. 7, pp. 93–99, 2016. View at: Publisher Site  Google Scholar
 D. González, S. Ruiz, M. GarcíaLozano, J. Olmos, and A. Serra, “System level evaluation of LTE networks with semidistributed intercell interference coordination,” in Proceedings of the IEEE 20th International Symposium on Personal, Indoor and Mobile Radio Communications, pp. 1497–1501, IEEE, Tokyo, Japan, September 2009. View at: Publisher Site  Google Scholar
 R. Jain, The Art of Computer Systems Performance Analysis: Techniques for Experimental Design, Measurement, Simulation, and Modeling, John Wiley & Sons, New York, NY, USA, 1990.
 J. C. Spall, Introduction to Stochastic Search and Optimization: Estimation, Simulation, and Control, vol. 65, John Wiley & Sons, 2005.
 D. G. Gonzalez, J. Hamalainen, H. Yanikomeroglu, M. GarciaLozano, and G. Senarath, “A novel multiobjective cell switchoff framework for cellular networks,” IEEE Access, vol. 4, pp. 7883–7898, 2016. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2017 David González González et al. 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.