Backhaul-Aware Dimensioning and Planning of Millimeter-Wave Small Cell Networks

The massive deployment of Small Cells (SCs) is increasingly being adopted by mobile operators to face the exponentially growing traffic demand. Using the millimeter-wave (mmWave) band in the access and backhaul networks will be key to provide the capacity that meets such demand. However, dimensioning and planning have become complex tasks, because the capacity requirements for mmWave links can significantly vary with the SC location. In this work, we address the problem of SC planning considering the backhaul constraints, assuming that a line-of-sight (LOS) between the nodes is required to reliably support the traffic demand. Such a LOS condition reduces the set of potential site locations. Simulation results show that, under certain conditions, the proposed algorithm is effective in finding solutions and strongly efficient in computational cost when compared to exhaustive search approaches.


Introduction
The deployment of the Fifth Generation (5G) of mobile communications has recently started, and it is envisaged to have a great impact on the digital society. In addition to the speed improvement of mobile broadband services, 5G will allow the definitive launch of massive machine communications and mission-critical machine communications [1,2]. The new services may require up to 10 Gbps peak data rate, 1000 times the bandwidth per unit area, up to 100 times the devices connected simultaneously, low latency (down to 1 ms), ultra-high reliability (up to 99.999% availability) and low power consumption (up to 10-year battery life). The use of high-band spectrum, virtualization, edge computing and network slicing are the main pillars of the unprecedented technological renewal required for current access, backhaul (BH) and core networks [3,4].
To cope with the explosive growth in demand for bandwidth, network densification using millimeter-wave (mmWave) Small Cells (SCs) is considered a key enabler of 5G systems [5,6]. SCs are short-range, low-cost, low-power base stations that provide additional capacity and extended coverage to the traditional macro-only cellular network. With dense SC deployments, the whole system throughput can be enhanced by increasing the frequency reuse. However, densifying the access network affects the BH network, which should be re-designed to avoid performance bottleneck, management complexity and high costs to the infrastructure provider [7]. While connecting SCs to the BH via optical fiber may result too expensive (e.g., in sparsely located fiber access areas), mmWave communications are an attractive alternative for supporting high bandwidth BH connections [8,9]. They can be established proposed a self-planning framework to update the cellular infrastructure in response to various events such as the deployment of a new tenant's slice (i.e., an isolated logical network). The cell deployment should be accompanied with intelligent resource scheduling and allocation in the BH to provide optimal network performance. In this line, the work in [33] proposes an optimization algorithm for steering data traffic of multiple slices in the BH.
Despite the existence of relevant work on cellular planning, the consideration of key 5G enablers (e.g., mmWave, massive MIMO, wireless BH, multi-tenancy, etc.) is scattered in the literature. Most of the above-mentioned works either address the planning problem in the access network (i.e., determine the cell location) or they attack problems in the BH network, such as routing and spectrum usage. The combination of these two problems is not sufficiently addressed. Although there are some proposals in the literature [7,21], their formulation neglects key considerations such as the multi-hop capability in the BH network. In addition, the spatial traffic distribution models used for performance evaluation are unlikely in reality.
In this work, we tackle the SC planning problem under two major assumptions: (i) SCs are wirelessly backhauled via the SC and/or MC layer; and (ii) a self-planning framework enables incremental deployment on top of the existing infrastructure. Unlike previous works, the first hypothesis implies two constraints at the same time: a maximum number of backhauled SCs per node and a maximum number of hops until reaching the MC. The second hypothesis allows considering a more realistic scenario, where a HetNet infrastructure is available to support the deployment of the new SCs. The main contributions of this work are as follows. We formulate it as a mixed integer non-linear programming problem, considering several key conditions that have not previously been considered together, such as the required LOS between nodes, the maximum number of hops, the number of simultaneous BH links and the number of SCs backhauled by every node. Then, we present an iterative two-step algorithm to solve the problem with limited complexity. The consideration of computationally-efficient methods is relevant since the wireless BH solutions have significantly lower associated time-to-deploy than wired options. In a first step, the heuristic algorithm provides an initial solution for satisfying the user demand. In a second step, the algorithm refines the solution to adapt it to the BH constraints. For the evaluation, we utilized standard models from 3GPP and a practical interference model for mmWave beamforming. The evaluation was carried out for different MC densities in the area of interest to establish various BH conditions that affect the behavior of the planning algorithm. In addition, a sensitivity analysis of the BH antenna gain as a key radio parameter was performed to study its effect on spectral efficiency. Simulation results show that the proposed algorithm can provide effective solutions with reduced computational costs and comparable performance to brute-force search approaches. This low complexity can ensure its application to large-scale scenarios. Moreover, the algorithm can easily be integrated with existing commercial planning tools, which usually combine live network measurements with predictions for accurate network planning.
The remainder of this paper is organized as follows. Section 2 starts with the system model and the formulation of the BH-aware cell dimensioning and planning problem. Section 3 describes the proposed computationally-efficient method considering both user demands and BH constraints. In Section 4, after describing the simulation scenario and setup, the performance of the proposed method is evaluated. Finally, Section 5 provides some concluding remarks.

System Model
Consider a two-tier HetNet scenario consisting of a set B = B M ∪ B S of base stations, where B M are the MCs and B S are the SCs, all managed by a certain operator. The access network serves a set U of user equipments (UEs), each demanding a capacity of D u . The sub-6 GHz and mmWave bands are used for access transmissions in MCs and SCs, respectively. In the BH network, each MC has a high-speed optical fiber connection, while each SC j can connect either to another SC or directly to a MC using an mmWave LOS link with a capacity C j BH . In this way, the SCs can act as aggregator nodes to relay UE data. The required LOS between the access and BH nodes limits the area of interest to a set S b LOS of candidate sites for placing the SC b. The spectrum in the access network is organized in resource blocks (RBs) that are allocated to the nodes in a semi-persistent manner to meet current demands [21]. Each UE u is allocated a fraction W u of the spectrum allocated to its serving node. Non-overlapped mmWave spectrum in the V-band is allocated for access and BH networks. The V-band is selected because of its unlicensed or lightly licensed spectrum use. Due to the attenuation effects of oxygen absorption, this band can be effectively used to connect closely spaced SCs with short links that are daisy-chained and aggregated for transport to the network core [4]. Let W BH be the bandwidth dedicated to the BH in each node. Beamforming is employed in the BH network to enable simultaneous transmissions (beams) to different nodes using the same RBs. However, SCs acting as relays in a multi-hop BH link will use time or frequency multiplexing for transmissions with the nodes at the previous and next hops. The SCs can communicate directly with the other nodes through the X2 protocol, which is an ultra-fast broadband-related protocol that allows mobile operators to use different topologies to offload traffic. For the access network, suppose that a UE u is served by certain node b. Let γ u,r be the Signal-to-Interference-plus-Noise Ratio (SINR) of this UE when transmitting on the RB r, defined as: where p RX u,b is the received power by the UE u from the node b, L j is the cell load factor of the node j, π j,r is a function that takes the value 1 when the RB r is allocated to the node j and the value 0 otherwise, I can take the value of M or S depending on the cell type (i.e., MC or SC) and p N is the noise power measured in one RB. The cell load factor is defined as the relation between the service demand and cell capacity according to the work in [32]. The spectral efficiency SE u,r of the UE u in the RB r is calculated from the γ u,r based on the following SINR mapping [34]: where SE max is the maximum achievable spectral efficiency with link adaptation; γ min and γ max are the minimum and maximum SINR values, respectively; and ρ specifies the attenuation factor, which represents implementation losses. The capacity offered by the node b to the UE u is expressed as: where R b is the subset of RBs allocated to the serving node b. The fraction W u allocated to the UE depends on the resource scheduling scheme. For example, assuming a round-robin scheme, it is given by: where W RB is the spectral occupancy of a RB (e.g., 180 kHz assuming a subcarrier spacing of 15 kHz) and U b is the subset of UEs served by the same node b.
For the BH network, the SCs are equipped with directional antennas whose gain pattern is modeled as in [35]: where θ m is the main-lobe width and A max and A min are the main and side-lobe gains, respectively. The antenna beams in the BH are aligned, so that the effective gain on BH links is A 2 max . Each candidate SC location in the area of interest has a probability of LOS propagation from a given BH node. This probability is taken from the 3GPP models in [36], which are applicable to frequencies up to 100 GHz. Unlike what happens with the UEs in the access links, an improvement of the LOS probability is applied to the BH links due to the local planning activities. According to the model presented in [37], the LOS probability p after local planning optimization is determined as follows: where p is the LOS probability without local planning optimization calculated as in [36] and τ is a parameter that depends on the site planning searching radius (e.g., τ = 3.1 for 50-m radius). In addition, the LOS/NLOS conditions at multiple positions are also correlated to provide spatial consistency.
Following the model in [38], the correlated LOS/NLOS conditions are generated by applying an exponential spatial filter to the independent random values with a given correlation distance (e.g., 50 m for 3GPP urban macro and micro scenarios [36]): where v c x ,c y and v c x ,c y are the independent and correlated LOS visibility variables, respectively, in the location (c x , c y ), where c x and c y are the xand y-coordinates. The parameter ∆d represents the correlation distance, d E stands for the Euclidean distance and M × N is the total number of grid points that are considered in the area of interest. Since the independent LOS visibility variable is Boolean to represent LOS and NLOS states, the obtained value of the correlated variable is rounded to the nearest 0 or 1. From this information, the set S b LOS of candidate sites for placing an SC can be obtained. The same exponential spatial filter is also applied to the shadow fading in order to provide spatial consistency. The SINR experienced by the SC b when served by its BH node h b is determined as: is the channel gain (i.e., path loss and shadow fading) between the SC and its BH node; L BH h j represents the cell load factor of the BH node h ; A h j (θ b ) and A b (θ h j ) are the transmitter and receiver antenna gains in the direction determined by the nodes b and h j , respectively; and p BH N stands for the noise power measured in the bandwidth of the BH network, W BH . Note that this expression assumes a single polarized MIMO system. Although cross-polarization helps reduce the undesired radio interference, it requires the use of special antennas and is more sensitive to interference and rain [39]. The cell load factor of BH nodes is given by the traffic load of their served SCs, connected with either direct or multi-hop links. The spectral efficiency SE BH b of the SC b is obtained using Equation (2). Lastly, the corresponding capacity C BH b offered by this SC is determined as: The operator offers the network infrastructure to multiple tenants. With the aim of satisfying the highly diversified radio access needs in space and time of these tenants, we adopt the self-planning framework presented in [32]. Under this framework, the aggregated UE demand and actual network capacity are compared to trigger particular planning actions such as deploying new SCs or reallocating the available spectrum.

Problem Description
In our problem formulation, we assume an existing infrastructure formed by a set B 0 = B M,0 ∪ B S,0 of base stations that has previously been deployed in certain locations s b0 , where b 0 ∈ B 0 stands for the node. We define the variable z s SC ∈ {0,1} at location s ∈ S to indicate with value 1 that an SC is deployed there. Likewise, z s MC ∈ {0,1} indicates the locations of the MCs. The variable x i,j ∈ {0,1} specifies that the node i is backhauled via the node j and y u,j ∈ {0,1} indicates that the UE u is served by the node j. For simplicity, we consider that the cell sizes and maximum distance for BH link availability are such that any SC can be backhauled to the MC through one SC at most, which implies the multi-hop BH link to be limited to two hops. In addition, the problem is focused on the downlink because of the asymmetric nature of data traffic. Given this, the SC dimensioning and planning problem can be expressed as follows: min subject to where α and β are parameters to adjust the required user throughput in the network and N S is the maximum number of SCs that can be backhauled by a certain SC. The objective is to minimize the deployment cost (i.e., the number of newly deployed SCs) while guaranteeing the UE demands.
In classical cellular planning, the deployment of new base stations is motivated by the insufficient quality of the radio links between the UEs and the base stations. In this way, Equation (11) represents the quality of service in the access network, unifying both coverage and capacity constraints for the UEs.
The cellular planning problem with wireless BH requires the consideration of a set of constraints. In particular, Equation (12) denotes the BH capacity constraint, which implies that an SC acting as relay aggregates the traffic from its backhauled SCs to that from its served UEs. Equation (13) establishes that the multi-hop BH link is limited to a maximum number of hops. Exceeding this condition may cause higher delays due to the processing time at each node. Equation (14) limits the number of SCs backhauled by one node. This condition help avoid the congestion of links in the BH network. With Equation (15), two SCs cannot be mutually assigned as their BH nodes. Breaking this condition would create undesirable loops in the BH network. Equations (16) and (17) ensure that every SC (UE) is backhauled (served) by exactly one node. Multiple links from the same SC can be used for load balancing or redundancy purposes. However, for simplicity, such objectives are out of the scope of this paper. Finally, Equations (18) and (19) represent the initial set of deployed SCs and MCs whose locations must remain fixed. As previously stated, this represents a realistic assumption for HetNets, where the current infrastructure is insufficient to cope with the increasing demand.

Proposed Dimensioning and Planning Algorithm
The formulated problem for cellular dimensioning and planning under non-ideal BH conditions is a mixed integer non-linear program. Thus, our proposed solution is a computationally-efficient BH-aware (CEBA) algorithm that follows the heuristic approach described in Algorithms 1 and 2, where the variable h b indicates the BH node for the node b and p b RX is the received power from b.
The UE demand D u is expressed as traffic density, with d x,y representing the traffic demand in the grid point (c x , c y ). The function F(·) is a filter to smooth the traffic demand variable. Dimensioning is based on iteratively adding a new SC to the planning problem until the network (both access and BH) capacity is enough to satisfy the UE demands (Steps 2-3 in Algorithm 1). The planning problem, where the number of SCs is invariant, is divided into two stages: the first aims at optimally placing the SCs based on the UE demands (Steps 4-7 in Algorithm 1), while the second further optimizes the location of the SCs to satisfy the required capacity in the BH links (Algorithm 2). The planning based on UE demands serves to increase the SC deployment in areas where the demand is high, making the cell size smaller. Following this principle, the k-means algorithm is first applied to identify high-and low-traffic subareas, resulting in a set of centroids as initial candidate locations for the SCs. The filter F(·) serves to remove irrelevant, small peaks of the traffic demand distribution before applying the k-means algorithm. The solution is refined through the steepest-ascent method to bring the SCs closer to the high-traffic subareas, achieving better spectral efficiency. Then, the location of some SCs are adjusted to match the location of the currently deployed SCs based on the minimum distance to them. With respect to parameter settings, the transmit power is set to have γ = 9 dB at √ 3/2 of the inter-site distance [32]. The number of RBs allocated to each SC is set to be proportional to the UE demand, and the RB selection is made as in [32] to reduce inter-cell interference. if b has NLOS BH then 3: create a set C of candidate BH nodes: C ← B\{b, nodes with NLOS BH or 2-hop BH link} 4: for c ∈ C 5: given if constraint in (14) is not met for b then 13: create a set G of SCs backhauled by b 14: for g ∈ G 15: create a set C of candidate BH nodes: C ← B\{b,g, nodes with NLOS BH or 2-hop BH link} 16: for c ∈ C 17: given sort g ∈ G by r g in descending order with new index g' 22: for g' ∈ G 23: if constraint in Equation (14)  Once the access network is planned, the wireless BH links are determined based on the maximum received power from candidate BH nodes having LOS. Candidate MCs are prioritized over SCs since the former implies direct links. In addition, SCs backhauled by another SC cannot serve as BH node because of the multi-hop constraint in Equation (13). It is possible that some SCs cannot have a LOS BH link from their location or that the condition in Equation (14) is not met. In these cases, Algorithm 2 is executed to make the BH-aware optimization process. The algorithm starts by modifying the candidate SC location in case of NLOS. Specifically, Steps 1-10 replace the infeasible location by a new location having LOS to the BH node that minimizes the distance between both locations. Since there is an initial set of deployed nodes, a valid location with LOS can be found out for backhauling any SC. Then, if an SC exceeds the maximum number of backhauled SCs, Steps 11-32 aim at finding another alternative for each backhauled SC until the constraint in Equation (14) is met. The location of these SCs can be modified to have LOS to the newly assigned BH node. The selection of the SCs to satisfy the constraint in Equation (14) is based on the minimum displacement to have a LOS position. Given that the MCs are not limited by the number of backhauled SCs, the algorithm will provide a solution that satisfies the constraint in Equation (14) at any case.

Simulation Scenario and Setup
The simulated scenario covers an urban area of 1 km × 1 km with a set of deployed MCs and SCs, which are located according to the spatial variations of the traffic demand. The deployment scenario is simulated following a snapshot-based model, where each snapshot involves a random realization of the same traffic probability distribution. The evaluation of a wide range of different situations ensures reliable statistical significance analysis. Let us assume that the expansion of a certain tenant's service results in a lack of capacity with the current infrastructure. The UE demand is randomly and non-uniformly distributed over the considered area. The configured UE service demand (20 Mbps) is representative for average human mobile or nomadic applications in busy hours, provided that virtual reality and high-definition video are among the most typical 5G mobile broadband services. With respect to the BH network, the antenna main-lobe width θ m is set to 10 • and the gains A max and A min are set to 12 and −2 dBi, respectively, unless stated otherwise [35]. The cell bandwidth W BH in the BH is fixed to 1 GHz, and it does not interfere with the access spectrum. Regarding the proposed algorithm, N S = 2 is set to limit the required capacity per BH link, α = 0.7 allows partial service degradation and β = 0.8 allows small congestion in BH links. The other relevant simulation parameters are shown in Table 1. As a reference for comparison, several methods are considered in the evaluation. Since solving the formulated problem requires impractical computational time, the optimal solution is approximated through a combination of exhaustive search and heuristic strategy. In particular, an exhaustive-search BH-aware (ESBA) and an exhaustive-search ideal-BH (ESIB) dimensioning and planning algorithms are developed following the same principles as in [16,32]. Specifically, the SCs are iteratively added one by one in locations that are selected by brute-force from the set of candidate locations to maximize the spectral efficiency in the access. In ESBA, the candidate locations are those having LOS BH and satisfying the BH-related constraints in Equations (12)- (17). The ESIB algorithm, which does not consider these constraints, serves to establish an upper-bound on the performance for the non-ideal BH scenario. In addition to these methods, a random-selection BH-aware (RSBA) dimensioning and planning algorithm is developed to establish an approximation lower-bound. In this case, the SCs are iteratively added one by one in random locations that satisfy the BH-related constraints.
The experiments were conducted on a Windows 10 desktop machine featuring an Intel Core i7-6700 @ 3.40 GHz CPU and 8.00 GB RAM. Figure 1 represents the deployment scenario resulting from applying each planning algorithm to support a total traffic demand of 5 Gbps. The colored contour lines indicate the traffic demand density. The dots represent the UE locations for a certain realization of the statistical distribution of the traffic demand. In the four cases, the algorithm starts from an initial set of deployed nodes, marked with blue triangles, which were derived from previous planning stage(s). The newly planned SCs are represented by red triangles. The grey lines indicate the BH links between nodes, except for ESIB where the BH is ideal (i.e., BH constraints are not considered). For CEBA, ESBA and ESIB, it is observed that the regions with higher traffic densities are covered by a larger number of SCs. In the case of CEBA and ESBA, these SCs are optimally placed to both increase resource usage efficiency and support BH links with LOS and complying with the maximum number of backhauled SCs and multi-hops. Additionally, it is observed that ESBA exhibits a lower cell density (or larger cell service areas) in low traffic regions than CEBA. Looking at the RSBA deployment, it is noted that the SCs are randomly placed in locations where the BH link is ensured, without keeping a relationship with the traffic demand.  Figure 2 compares the performance of the four algorithms. In particular, the number of planned SCs and the average spectral efficiency in the access network are shown for different MC densities and different magnitudes of the traffic demand. Varying the MC density serves to establish different BH conditions, while varying the traffic demand allows analyzing capacity issues in the scenario. As expected, RSBA and ESIB establish lower and upper bounds, respectively, on the level of performance reachable by the methods. The optimal method, ESBA, approximates the performance of ESIB as the MC density increases, showing that the BH constraints are less relevant for higher densities. In particular, for 2 MCs/km 2 , the impact of the BH network is appreciable on both indicators by comparing ESBA with ESIB. To understand the increased number of SCs, note that the distance between MCs and SCs is limited by the maximum allowed path loss. Consequently, for a low number of MCs in the scenario, multi-hop BH links may frequently be required to reach distant areas, increasing the number of SCs. For 3 MCs/km 2 , the BH impact is appreciable on the spectral efficiency and insignificant on the number of SCs. Since in this case there are more locations satisfying the LOS condition, multi-hop BH links are less probable. Thus, ESBA approximates the same number of SCs as in the ideal case given by ESIB. However, deviations from the optimal location due to the BH constraints still impact the spectral efficiency as shown in Figure 2d. The proposed method, CEBA, is more sensitive to the BH conditions than ESBA. For 2 MCs/km 2 , CEBA increases the number of SCs, on average, by 26% in comparison with ESBA. Such an increase is accentuated for higher traffic demands. Looking at the spectral efficiency, as the traffic demand increases, CEBA approaches the ESBA's curve as a consequence of the increasing number of deployed SCs. This behavior of CEBA changes drastically for MC densities above 2 MCs/km 2 . Specifically, for 3 MCs/km 2 and low-medium traffic demands, the number of SCs is only increased by 10% compared to ESBA (and ESIB). It also deviates significantly from the random case. This is due to the smaller distances between MCs and SCs, which minimize the SC relocations made by CEBA.

Simulation Results and Discussions
For high traffic demands (above 5-6 Gbps), the slope of the CEBA's curve representing the number of SCs experiences a moderate increase. This effect is appreciable for the two values of MC densities shown in Figure 2. This is mainly because of the presence of BH links with insufficient capacity to meet the demand (e.g., due to traffic aggregation in multi-hop BH links), even though the LOS condition is satisfied. To explain this, consider the access and BH capacity conditions (Equations (11) and (12), respectively), which take part in the main loop of CEBA's Algorithm 1. It is observed (not shown for brevity) that the first iterations of this loop are triggered by breaking Equation (11), while the last iterations are triggered by breaking Equation (12). Since the bandwidth in the BH is not modified, the only solution is to increase the number of SCs to distribute the access demand among more BH links until Equation (12) is satisfied. This issue can also be solved by increasing the antenna gain as discussed in the next paragraph.
To study the effect of the directional high-gain and narrow-beam antennas on BH network performance, three different values of the antenna gain are evaluated for the scenario with 3 MCs/km 2 and CEBA planning method. Figure 3 shows the number of planned SCs and the average spectral efficiency measured in the BH network. It is observed that, for high traffic demands, the number of planned SCs is smaller for 18 and 25 dBi gains. This is due to an increase of about 1 bps/Hz in the spectral efficiency, as shown in Figure 3b. In particular, when the BH capacity condition defined in Equation (12) is checked in the main loop of CEBA's Algorithm 1, it returns a true value, indicating that the BH links are not congested and, therefore, additional SCs are not needed. The improvement in the average spectral efficiency is more pronounced in the range of low-to-medium values of the antenna gain. Accordingly, 18 dBi gain provides a good trade-off between spectral efficiency and antenna complexity.  Table 2 shows a comparison of the planning algorithms in terms of computational complexity. Specifically, it provides the simulation time of each method for different traffic demands and an MC density of 2 MCs/km 2 . In general, this time increases with the traffic demand since more iterations are required to place additional SCs. Looking at the differences between methods, RSBA gives the lowest values since it follows a very simple logic for cell planning. Consequently, most of the simulation time is spent in network performance computation. The proposed method, CEBA, increases the simulation time with respect to RSBA by a factor of 2-3. This increment is mainly determined by the execution of the k-means, the steepest-ascent and the BH-aware optimization algorithms. Nevertheless, such an increase factor reveals a good trade-off between complexity and effectiveness of the solution. On the contrary, ESBA and ESIB provide the greatest values of the simulation time due to the high computational load of the exhaustive search approach. In particular, the increase factors are 25-37, which are considerably higher than those obtained by CEBA. Contrary to expectations, ESIB achieves larger simulation times than ESBA even though the BH network is not simulated. The underlying reason is that the set of candidate solutions that are evaluated by brute-force in the case of ESBA is smaller due to the required LOS condition. The simulation time is also impacted by the scenario size and the grid size, since they determine the number of candidate locations. The required computation resources for exhaustive-search methods such as ESBA and ESIB represent a bottleneck, as their scalability is very limited when large-scale scenarios are considered. Thus, the benefits of using computationally-efficient methods such as CEBA are very valued.

Conclusions
In this work, the authors propose a computationally-efficient algorithm for the dimensioning and planning of SCs under non-ideal BH conditions. The main contributions to this field rely on the consideration of a set of constraints in the BH that limit the location of the SCs, which are deployed on top of an existing HetNet infrastructure that is available to support the deployment of the new SCs. Specifically, the algorithm defines two phases to first determine a candidate solution that satisfies the user demand and then modify the solution to meet the BH constraints. Since some SCs may not have a LOS BH link from their initial candidate location, the BH-aware optimization modifies such locations based on the minimum displacement to have a LOS position. Simulation results show that, for sufficient MC density and low-medium traffic demands, the performance of the solution in terms of number of deployed SCs approximates those obtained by exhaustive-search methods, using much less computational resource.