Dynamic Base Station Operation in Large-Scale Green Cellular Networks Yue Ling Che, Lingjie Duan, and Rui Zhang
arXiv:1512.07469v1 [cs.IT] 23 Dec 2015
Abstract In this paper, to minimize the on-grid energy cost in a large-scale green cellular network, we jointly design the optimal BS on/off operation policy and the on-grid energy purchase policy from a network-level perspective. We consider the BSs are aggregated as a microgrid with hybrid energy supplies and an associated central energy storage (CES), which can store the harvested renewable energy and the purchased on-grid energy over time. Due to the fluctuations of the on-grid energy prices, the harvested renewable energy, and the network traffic loads over time, as well as the BS coordination to hand over the traffic offloaded from the inactive BSs to the active BSs, it is generally NP-hard to find a network-level optimal adaptation policy that can minimize the on-grid energy cost over a long-term and yet assures the downlink transmission quality at the same time. Aiming at the network-level dynamic system design, we jointly apply stochastic geometry (Geo) for large-scale green cellular network analysis and dynamic programming (DP) for adaptive BS on/off operation design and on-grid energy purchase design, and thus propose a new Geo-DP design approach. By this approach, we obtain the optimal BS on/off policy, which shows that the optimal BSs’ active operation probability in each horizon is just sufficient to assure the required downlink transmission quality with time-varying load in the large-scale cellular network. However, due to the curse of dimensionality of the DP, it is of high complexity to obtain the optimal on-grid energy purchase policy. We thus propose a suboptimal on-grid energy purchase policy with low-complexity, where the low-price on-grid energy is over-purchased in the current horizon only when the current storage level and the future renewable energy level are both low. Simulation results show that the suboptimal on-grid energy purchase can achieve near-optimal performance. We also compare the proposed policy with the existing schemes, to show that our proposed policy can more efficiently save the on-grid energy cost over time. Index Terms Base station on/off operation, hybrid energy supplies, on-grid energy cost minimization, energy storage management, stochastic geometry, dynamic programming.
I. I NTRODUCTION The dramatically increased mobile traffic can easily lead to an unacceptable total energy consumption level in the future cellular network, and energy-efficient cellular network design has become one of the critical issues for developing the 5G wireless communication systems [1]. Due to the resulting high energy cost and the CO2 Y. L. Che and L. Duan are with the Engineering Systems and Design Pillar, Singapore University of Technology and Design (e-mail: yueling
[email protected]; lingjie
[email protected]). R. Zhang is with the Department of Electrical and Computer Engineering, National University of Singapore (e-mail:
[email protected]). He is also with the Institute for Infocomm Research, A*STAR, Singapore.
2
emissions, it is important to improve the cellular network energy efficiency (EE) by greening the cellular network [2]. It has been noted that the operations of base stations (BSs) contribute to most of the energy consumptions in the cellular networks, which are estimated as 60-80 percent of the total network energy consumptions [3]. Hence, reducing the BSs’ energy consumptions is crucial for developing energy-efficient or green cellular networks. In traditional cellular networks, BSs are usually densely deployed to meet the peak-hour traffic load, which causes unnecessary energy consumption for many lightly-loaded BSs during the low traffic load period (e.g., early morning). It is envisioned that in the future energy-efficient 5G wireless communication system, the cellular network should dynamically shut down some BSs or adapt the BS transmissions according to the time-varying traffic load. Although appealing, due to the significant spatial and temporal fluctuations of the traffic load in the cellular networks, optimal BS on/off operation design is a challenging problem [4]. For example, the authors in [5] studied optimal energy saving for the cellular networks by reducing the number of active cells, where simplified cellular network model with hexagonal cells and uniformly distributed traffic load is assumed. In [8], the authors proposed adaptive cell zooming scheme according to the traffic load fluctuations, but without considering downlink transmission quality explicitly. In [6], by studying the impact of each BS’s active operation on the network performance, the authors proposed a suboptimal BS on/off scheme that can be implemented in a distributed manner to save the BSs’ energy consumption. However, the impact of traffic fluctuation on the BS’s on/off decisions as well as the downlink transmission quality were not explicitly considered. The issue of downlink transmission quality was considered in [7], where the BS’s transmit power and cell range adaptations to the traffic load fluctuation were proposed to improve the EE of the cellular network, but [7] only considered a single-cell scenario. In [9], an energy group-buying scheme was proposed to let the network operators make the day-ahead and real-time energy purchase as a single group, where the BSs of the network operators share the wireless traffic to maximally turn lightly-loaded BSs into sleep mode. To more effectively reduce the carbon footprint and improve the EE, the energy harvesting (EH) enabled BSs, which are able to harvest clean and free renewable energy (e.g., solar and wind energy) from the surrounding environment for utilization, have been developed for the green cellular networks recently. For example, in [10], by assuming a renewable energy storage, the authors applied the dynamic programming (DP) technique to study the energy allocations for a point-to-point transmission, to adapt to the time-selective channel fading. In [11], by using a birth-death process to model the harvested renewable energy level at each BS for a large-scale K -tier green cellular network, the authors characterized an energy availability region of the BSs based on tools from stochastic geometry [12]. Various issues of the EH-enabled wireless communications are also discussed in [13]-
3
[16]. However, it is noted that the available renewable energy in nature is limited and intermittent, and thus cannot always guarantee sufficient energy supply to power up the BSs. By combining the merits of the reliable on-grid energy and the cheap renewable energy, hybrid energy supplied BS design has become a promising method to power up the BSs [17]. One of the main design issues is to properly exploit the available renewable energy to minimize the on-grid energy cost or maximize the EE of the system. We note that some initial work has been proposed to address this issues via, e.g., the optimal packet scheduling [18], the optimal BS resource allocation [19], and low-complexity online algorithm design based on the Lyapunov optimization [20]. In [21]-[23], to mitigate the fluctuations of the renewable energy harvested by each BS, the authors also considered energy cooperation and sharing between the BSs. Moreover, by applying the DP technique to adapt to the fluctuations of both the traffic loads and the harvested renewable energy, hybrid energy supplied BS on/off operation design has also been studied in the literature. For example, in [24], by assuming each BS in the network is supplied by either the renewable energy or the on-grid energy, the authors studied the BS sleep control problem to save the energy consumptions of all the BSs. In [25], by assuming each BS is jointly supplied by renewable energy and on-grid energy, the authors studied the joint optimization of the BS on/off operation and resource allocation under a user blocking probability constraint at each BS. Due to the traffic offloaded from the inactive BSs to the active BSs, all the BSs in the network are essentially coordinated to support to all the users in both [24] and [25]. However, due to the dimensionality curse of the DP [26], it is generally NP-hard to conduct a network-level BS coordination to decide their optimal on/off status. As a result, in [24] and [25], the authors proposed a cluster-based BS coordination scheme to decide the BSs’ on/off status in a large-scale network, where the BSs form different clusters, and the BSs in each cluster apply the DP technique to decide their on/off status. However, such a cluster-based BS coordination scheme cannot achieve network-level optimal BS on/off decisions in general. In this paper, we consider a large-scale cellular network, where all the BSs are aggregated as a microgrid with hybrid energy supplies and an associated central energy storage (CES). The microgrid works in island mode if it has enough local supply from the renewable farm; otherwise, it connects to the main grid for purchasing the on-grid energy [27]. By adapting to the fluctuations of the harvested renewable energy, the network traffic loads, as well as the on-grid energy prices over time, the microgrid applies the DP technique to decide the amount of on-grid energy to purchase to satisfy all the BSs’ energy demand in the network, as well as the probability that each BS stays active, which is referred to as the BSs’ active operation probability, at each time, so as to minimize the total on-grid energy cost over time in the network. Unlike [24] and [25], we aim at a network-level optimal BS
4
on/off operation policy design, such that the on-grid energy cost in the cellular network can be efficiently saved. The key contributions of this paper are summarized as follows: •
Novel large-scale network model with hybrid-energy supplied BSs: Based on Poisson point processes (PPPs), we first model the large-scale network by jointly considering the mobility of mobile terminals (MTs), traffic offloading from an inactive BS to its neighboring active BSs, as well as BS frequency reuse for reducing the downlink network interference level. Since all the BSs are supplied by the hybrid energy stored in the CES, we then model the hybrid energy management at the CES, where the on-grid energy is purchased with time-varying price to ensure that the energy demands of all the BSs in the network can be satisfied.
•
New Geo-DP approach for dynamic network optimization: To pursue a network-level design, we jointly apply stochastic geometry (Geo) for large-scale green cellular network analysis and the DP technique for adaptive BS on/off operation and on-grid energy purchase, and thus propose a new Geo-DP approach. By this approach, the microgrid centrally decides the optimal BSs’ active operation probability and the purchased on-grid energy in each time horizon, so as to minimize the total on-grid energy cost over all time horizons, under the BSs’ downlink transmission quality constraint and the BSs’ energy demand constraint.
•
Network-level optimal policy design: We aim at optimal offline policy design, which updates the BSs’ on/off status and on-grid energy purchase decisions in a predictable future period. By studying the impact of the BSs’ active operation probability on the storage level at the CES, we find an optimal BS on/off policy, which shows that the optimal BSs’ active operation probability in each horizon is just sufficient to assure the required downlink transmission quality with time-varying load in the large-scale cellular network. However, due to the curse of dimensionality of the DP, it is of high complexity to obtain the optimal on-grid energy purchase policy. We thus propose a suboptimal on-grid energy purchase policy with low-complexity, where the lowprice on-grid energy is over-purchased in the current horizon only when the current storage level and the future renewable energy level are both low.
•
Efficient network-level BS coordination to minimize the on-grid energy cost: By conducting extensive simulations in Section VI, we show that the proposed suboptimal on-grid energy purchase policy can achieve near-optimal performance. We also show that our proposed policy is robust to the prediction errors of the on-grid energy prices, the harvested renewable energy, and the network traffic loads, and thus can be properly implemented in practice. Moreover, we conduct a large-scale network simulation to compare our proposed policy with the existing schemes, and show that our proposed policy can more efficiently save the on-grid energy cost in the network.
5
Fig. 1.
Illustration of the cellular network supplied by both renewable energy and on-grid energy.
As a powerful tool, stochastic geometry has been widely applied to model the EH enabled wireless communication networks. For example, besides [11], we also noticed that the EH-based ad hoc networks, cognitive radio networks, cooperative communication networks have been studied in [28], [29], and [30], respectively, by adopting the PPP-based network modeling. In our previous work [31], we also applied tools from stochastic geometry to analyze the spatial throughput of a wireless powered communication network, where the users exploit their harvested radio frequency (RF) energy to power up their communications. Although the network-level system performance can be characterized in these existing studies, the adaptation to the fluctuations of the network parameters (e.g., the harvested renewable energy, the network traffic loads, and the on-grid energy prices) has not been properly addressed. To our best knowledge, our newly proposed Geo-DP approach that can lead to a network-level optimal adaptation policy design has not been studied in the literature. II. N ETWORK O PERATION M ODEL We consider a large-scale cellular network, where the traffic loads, the on-grid energy prices, the available renewable energy, as well as the wireless fading channels from each BS to the MTs are all varying over time in the network. The BSs are able to switch between active (i.e., on status) and inactive (i.e., off status) operation modes to save the total on-grid energy cost over time. Supported by smart grid development, as illustrated in Fig. 1, we aggregate all the BSs as a microgrid with hybrid energy supplies and an associated CES to defend against the supply/load variation [27]. For efficient energy management, the cellular network deploys the CES in the microgrid, which connects to both the local renewable farm (e.g., wind turbines or solar panels) and the power grid, and thus can store the harvested renewable energy (e.g., wind or solar energy) from the environment and the
6
Fig. 2.
Dynamic network operations over two different time scales: horizons and slots.
purchased on-grid energy from the power grid at the same time1 . All the BSs are connected with the CES and are powered up by its stored hybrid energy. Similar to [7], it is assumed that the CES is able to charge and discharge at the same time. In the following subsections, we first model each BS’s operations over time, and then present the large-scale cellular network model based on homogeneous PPPs. The hybrid energy management model at the CES will be elaborated later in Section III. A. Time Scales for Network Operations As compared to the variations of the downlink wireless fading channels from each BS to the MTs, the renewable energy arrival rates, the network traffic loads, as well as the on-grid energy prices are all varying slowly over time [11]. We thus consider the network operates over two different time scales [15]: one is the long time scale and is referred to as horizons, and the other is the short time scale and is referred to as slots, as shown in Fig. 2. We define a horizon as a reasonably large time period where the renewable energy arrival rate, the network traffic load, and the on-grid energy price all remain unchanged. Each horizon consists of N slots, and the wireless fading channels vary over slots. We focus on in total T horizons, 1 ≤ T < ∞, and assume each horizon has a fixed duration τ > 0. For simplicity, we assume τ = 1 hour in the sequel, and thus interchangeably use Wh and W to measure the energy consumption in each horizon. By adapting to the variations of the on-grid energy prices, the renewable energy arrival rates, the network traffic loads, as well as the wireless fading channels, the microgrid centrally decides the amount of on-grid energy to purchase and the active operation probability of the BSs in the network for each horizon t ∈ {1, ..., T }, which are denoted by G(t) > 0 and ρ(t), respectively, so as to minimize the total on-grid energy cost over all T horizons. 1
A CES is easier and cheaper to manage as compared to distributed storage at each BS. One can also view this CES as a collection of
distributed storage with perfect energy exchange links.
7
Similar to [13], we consider the optimal offline policy design, and assume the fluctuations of the renewable energy arrival rates, the network traffic loads, as well as the on-grid energy prices in all T horizons can be accurately predicted [32]. We will show that even given reasonable prediction errors, our approach works quite well in Section VI by simulation. In each horizon t, based on ρ(t) determined by the microgrid, each BS independently decides to be active or inactive. If a BS decides to be active, it transmits to its associated MTs with a fixed power level [33], denoted by PB > 0. If a BS decides to be inactive, it keeps silent and its originally associated MTs are handed over to its neighboring active BSs. B. PPP-based Network Model This subsection models the large-scale cellular network based on stochastic geometry. We assume the BSs are deployed randomly according to a homogeneous PPP, denoted by Φ(λB ), of BS density λB > 0. For a given active operation probability ρ(t) in horizon t, due to the independent on/off decisions at the BSs, the point process formed by the active BSs in each horizon t ∈ {1, ..., T } is still a homogeneous PPP, denoted by Φ(λB ρ(t)), with active BS density given by λB ρ(t). We assume the MTs are also distributed as a homogeneous PPP, denoted by Φ(λm (t)), of MT density λm (t) > 0 in horizon t. All Φ(λm (t))’s, ∀t ∈ {1, ..., T }, and Φ(λB ) are mutually independent. The MTs independently move in the network over slots in each horizon based on the random walk model [34]. Each MT is associated with its nearest active BS from Φ(λB ρ(t)) in each slot of horizon t, for achieving good communication quality. When an active BS becomes inactive, its originally associated MTs will be transferred to their nearest BSs that are still active as in [5]-[8] and [18]-[25]. From [35], the resulting coverage areas of the active BSs in each slot comprise a Voronoi tessellation on the plane R2 . We refer to traffic load of an active BS as the number of MTs that are associated with it, i.e., those are located within the Voronoi cell of the active BS, as in [11]. The following proposition studies the average traffic load of an active BS. Proposition 2.1: Given the BSs’ active operation probability ρ(t) in horizon t, ∀t ∈ {1, ..., T }, the average traffic load for an active BS from Φ(λB ρ(t)) is given by D(t) =
λm (t) λB ρ(t) .
The proof of Proposition 2.1 is obtained by noticing that D(t) equals the product of the MT density λm (t) and the average coverage area of each BS, where the latter term is
1 λB ρ(t)
from Lemma 1 in [11], and thus is
omitted here for brevity. From Proposition 2.1, the average traffic load D(t) at each BS is increasing over the MT density λm (t), and decreasing over the BS density λB as well as the BSs’ active operation probability ρ(t), as expected. Moreover, since the inactive BSs do not transmit information to the MTs and thus have zero traffic load, the average traffic load for each BS from Φ(λB ) is also obtained as D(t) for each horizon t ∈ {1, ..., T }.
8
C. Frequency Reuse Model To effectively control the potentially large downlink interference, which is mainly caused by the nearby BSs that are operating over the same frequency band, this subsection presents a random frequency reuse model to assign different frequency bands to the neighboring BSs. Specifically, we assume the total bandwidth available for the network is W Hz. To properly support the traffic load D(t), we assume the total bandwidth is equally divided into δ(t) bands in each horizon t, and each BS is randomly assigned with one out of δ(t) bands for operation as in [35]. We consider orthogonal multiple access for the MTs in the Voronoi cell of each active BS as in, e.g., [11] and [35], and each MT is assigned with a channel with fixed bandwidth B Hz for communication, where 0 < B < W . Since each active BS is only assigned with one band, by letting each BS’s required total bandwidth
to support its associated MTs be equal to the average bandwidth per band, i.e., D(t)B = W/δ(t), we obtain δ(t) =
by substituting D(t) =
λm (t) λB ρ(t)
W 1 W λB ρ(t) = , B D(t) B λm (t)
(1)
from Proposition 2.1. It is observed from (1) that the number of bands δ(t) decreases
as the traffic load D(t) increases, as expected. For each horizon t ∈ {1, ..., T } with a given δ(t), due to the BSs’ independent operation band selection as well as the independent on/off decisions over horizons, the point process formed by the co-band active BSs is still a homogeneous PPP, denoted by Φ(λcB (t)), with co-band active BS 1 density λcB (t) = λB ρ(t) δ(t) . As a result, for each MT, it only receives interference from all other active BSs that
operate over the same band as its associated BS with probability 1/δ(t). In the next subsection, based on the PPP-based network model and the frequency reuse mode, we present the downlink transmission model from each BS to its associated MT. D. Downlink Transmission Model We measure the downlink transmission quality between the BS and its associated MTs based on the received signal-to-interference-plus-noise-ratio (SINR) at the MTs. We assume Rayleigh flat fading channels with path-loss as in [11], [28], [31], and [35], where the channel gains vary over slots as shown in Fig. 2. Let Φ(λB ) = {X}, X ∈ R2 , denote the coordinates of the BSs. In slot n ∈ {1, ..., N } of horizon t ∈ {1, ..., T }, let hX (t, n) represent
an exponentially distributed random variable with unit mean to model Rayleigh fading from BS X to the origin o = (0, 0) in the plane R2 . We assume hX (t, n)’s are mutually independent over any time slot n ∈ {1, ..., N }, any
horizon t ∈ {1, ..., T }, and any BS X ∈ Φ(λB ). Suppose a typical MT is located at the origin and is associated with the active BS X ∈ Φ(λB ρ(t)) in slot n of horizon t. We express the desired signal strength received at the typical MT as |X|−α hX (t, n), where |X| is the distance from the associated active BS X to the origin, α > 2
9
is the path-loss exponent. Under the random frequency reuse model, since the received interference at the typical MT is caused by the co-band BSs from Φ(λcB (t)), we express the received interference at the typical MT as P −α h (t, n). Thus, the received SINR at the typical MT from BS X in slot n of horizon Y Y ∈Φ(λcB (t)),Y 6=X PB |Y | t, denoted by SINRo (t, n), is obtained as
SINRo (t, n) = P
PB |X|−α hX (t, n) , ∀t ∈ {1, ..., T }, ∀n ∈ {1, ..., N }, −α h (t, n) + σ 2 Y Y ∈Φ(λcB (t)),Y 6=X PB |Y |
(2)
where σ 2 is the noise power at the typical MT. It is also noted that |X| in (2) also represents the minimum distance from the typical MT to all the active BSs from Φ(λB ρ(t)), regardless of their operating frequency bands, i.e., no other active BSs can be closer than |X|. By using the null probability of a PPP, the probability density 2
function (pdf) of |X| in (2) is then obtained as f|X| (r) = 2πλB ρ(t)re−λB ρ(t)πr for r ≥ 0 [35]. We say a downlink transmission is successful if the received SINR at the MT is no smaller than a predefined threshold β > 0. III. H YBRID E NERGY M ANAGEMENT M ODEL This section presents the hybrid energy management model at the CES, by modeling the on-grid and renewable energy consumptions of all the BSs in the network. In the following, we first model the energy consumption at both active and inactive BSs. We then detail the hybrid energy management at the CES so as to meet the energy demands of the BSs. A. Energy Consumption at the BS According to the BS energy consumption breakdown provided by the practical studies in [36], various components of the BS, e.g., the power amplifier, base-band operation, and power supply for power grid and/or back-haul connections, all contribute to the total BS energy consumption. It is also shown in [36] that the BS energy consumption can be well approximated by a linear power model, where the total energy consumption of an active transmitting BS is affine/linear to its traffic load, with the corresponding slope determined by the power amplifier efficiency. Such a linear power model has been widely adopted in the literature (see, e.g., [10] and [25]). We thus also adopt it to model the energy consumption for each BS in this paper. Specifically, for a given horizon t, although the traffic load of each active BS generally varies over space, by applying our PPP-based network model in Section II, we can obtain the average energy consumption of each active BS over space as Pon (t) = Pa +
PB D(t), t ∈ {1, ..., T }, µ
(3)
where the average traffic load D(t) of the active BS is given in Proposition 2.1, µ ∈ (0, 1) is the power amplifier efficiency, and Pa > 0 represents the active BS’s average energy consumption for supporting its basic operations
10
and is taken as a constant [36]. For an inactive BS in horizon t, since it keeps silent in the horizon, we model its average energy consumption as a constant, which is given by Pof f (t) = Ps , t ∈ {1, ..., T },
(4)
where similar to Pa in (3), the constant Ps > 0 represents the average energy consumption of an inactive BS for supplying its basic operations. We assume Pa > Ps [36]. Based on the energy consumption model for each active or inactive BS, we present the hybrid energy management model at the CES to meet all the BSs’ energy demand. B. Energy Supplies at the CES As shown in Fig. 1, the CES can store both the renewable energy that is harvested from the local renewable farm and the on-grid energy that is purchased from the power grid. In each horizon, the microgrid centrally decides the amount of on-grid energy to purchase from the power grid and the BSs’ active operation probability ρ(t), by jointly considering the BSs’ downlink transmission quality, the energy demands of all the BSs in each horizon, and the hybrid energy storage level at the CES. We focus on the average energy consumption of all the BSs normalized by spatial area. For a given active operation probability ρ(t) in horizon t, by considering the BS density λB and its average traffic load D(t) in each horizon t, we can obtain the average energy consumption of all the BSs over the BSs’ on/off status and the network area in each horizon as E(t) = λB ρ(t)Pon (t) + λB (1 − ρ(t))Pof f (t) (a)
= λB ρ(t)Pgap + λB Ps +
PB λm (t) , t ∈ {1, ..., T }. µ
(5) (6)
where equality (a) follows by substituting (3) and (4) into (5), and Pgap = Pa − Ps > 0. It is easy to find that the unit of E(t) is given by W/unit-area. It is also observed from (6) that E(t) increases monotonically over the BSs’ active operation probability ρ(t) in each horizon t. Denote the storage level of the CES at the beginning of horizon t as B(t) ≥ 0, the amount of on-grid energy to purchase in horizon t as G(t), and the renewable energy arrival rate in horizon t as λe (t). All B(t), G(t), and λe (t) have the same unit as E(t), i.e., W/unit-area. According to Fig. 1, to satisfy the demanded energy of all the
BSs in horizon t, we obtain the following energy demand constraint: B(t) + G(t) + λe (t) ≥ E(t), t ∈ {1, ..., T }.
(7)
If the purchased on-grid energy G(t) and/or harvested renewable energy λe (t) cannot be completely consumed by the BSs in the current horizon, their residual portions are both stored in the CES for future use. The dynamics of
11
the storage level B(t) is obtained as B(t + 1) = min (B(t) + λe (t) + G(t) − E(t), C) ,
(8)
where 0 < C < ∞ is the given storage capacity of the CES. For simplicity, we assume the initial storage level at the beginning of horizon t = 1 is B(1) = 0. IV. P ERFORMANCE M ETRICS
AND
P ROBLEM F ORMULATION
In this section, based on the network operation model and the hybrid energy management model in Sections II and III, respectively, we propose a new Geo-DP approach to study the cellular network and formulate the on-grid energy cost minimization problem from the network-level perspective. In particular, we first use stochastic geometry to characterize the downlink transmission performance from the BSs to their associated MTs. We then further apply DP to formulate the on-grid energy cost minimization problem under the BSs’ successful downlink transmission probability constraint. A. Successful Downlink Transmission Probability As discussed in Section II, if a BS is active in a particular horizon, it transmits to its associated MTs in each slot within this horizon, while the channel fadings from the BSs to the MTs as well as the locations of the MTs vary over slots. In this subsection, by applying tools from stochastic geometry, we study the active BSs’ successful downlink transmission probability to their associated MTs in each slot. Specifically, from Section II-B, due to the stationarity of the considered homogeneous PPPs for the active BSs and the MTs in each horizon as well as the independent moving of the MTs over slots, we focus on a typical MT in each slot of horizon t ∈ {1, ..., T }, which is assumed to be located at the origin of the plane R2 , without loss of generality. Denote the successful downlink transmission probability for the typical MT in slot n of horizon t as Psuc (t, n), ∀t ∈ {1, ..., T }, ∀n ∈ {1, ..., N }. According to Section II-D, for a targeted SINR threshold β , we define Psuc (t, n) as Psuc (t, n) = P (SINRo (t, n) ≥ β) , ∀t ∈ {1, ..., T }, ∀n ∈ {1, ..., N },
(9)
with SINRo (t, n) for a typical MT given in (2). By applying the probability generating functional (PGFL) of a PPP [12], we explicitly express Psuc (t, n) in the following proposition. Proposition 4.1: In slot n of horizon t, ∀t ∈ {1, ..., T }, ∀n ∈ {1, ..., N }, the successful downlink transmission probability for the typical MT is Psuc (t, n) = πλB ρ(t)
Z
∞ 0
e−ax e−bx
α/2
dx,
(10)
12 (t)B with a = πλB ρ(t) + π λmW v , where B and W are the fixed channel bandwidth for each MT and the total R∞ bandwidth available in the network, respectively, as introduced in Section II-C, and v = β 2/α β −2/α 1+u1α/2 du,
and b = βσ 2 /PB . When α = 4, (10) admits a closed-form expression with s Υ(t)2 πPB Psuc (t, n) = πλB ρ(t) exp Q(Υ(t)), βσ 2 2
where Υ(t) =
λm (t)B v) (πλB ρ(t)+π √ W
√
2βσ2
PB
with v =
is the standard Gaussian tail probability.
√ √ β π/2 − arctan 1/ β , and Q(x) =
(11) √1 2π
R∞ x
2 exp − u2 du
Proposition 4.1 is proved using a method similar to that for proving Theorem 2 in [35], and thus is omitted for brevity. By using a similar method as in our previous work [31], one can also easily validate Proposition 4.1 by simulation. Moreover, from Proposition 4.1, due to the independent moving of the MTs as introduced in Section II-B, the expressions of Psuc (t, n) are observed to be identical for all slots within each horizon t ∈ {1, ..., T ]}. Thus, for notational simplicity, we omit the slot index n and use Psuc (t) to represent the typical BS’s successful downlink transmission probability in any slot of horizon t in the sequel. B. Problem Formulation Denote the on-grid energy price in each horizon t as a(t) ($/W). The total on-grid energy cost across all the BSs P over all T horizons is thus obtained as Tt=1 a(t)G(t), with unit $/unit-area. To minimize the total on-grid energy cost under the energy demand constraint given in (7), the microgrid can take advantage of the on-grid energy’s
price variations by purchasing more on-grid energy when its price is low and store it for future use, and purchasing less (or even zero) on-grid energy when its price is high. To ensure the QoS for each BS’s downlink transmissions, we also apply a successful downlink transmission probability constraint such that Psuc (t) ≥ 1 − ǫ with ǫ ≪ 1 and t ∈ {1, ..., T } [35]. As a result, by jointly considering the energy demand constraint and the successful downlink
transmission probability constraint, we optimize the BSs’ active operation probability ρ(t) and the purchased on-grid energy G(t) in each horizon, and formulate the on-grid energy cost minimization problem as (P1) :
min.
ρ(t),G(t)
T X
a(t)G(t)
t=1
s.t. Psuc (t) ≥ 1 − ǫ, ∀t ∈ {1, ..., T }
(12)
B(t) + G(t) + λe (t) ≥ E(t), ∀t ∈ {1, ..., T } G(t) ≥ 0, ∀t ∈ {1, ..., T }
(13)
0 ≤ ρ(t) ≤ 1, ∀t ∈ {1, ..., T },
(14)
B(t + 1) = min (B(t) + λe (t) + G(t) − E(t), C) , ∀t ∈ {1, ..., T − 1},
13
where E(t) is determined by ρ(t) and is given in (6). However, for a general α, due to the complicated expression of Psuc (t) in (10), it is generally difficult to analyze the effects of ρ(t) on the successful downlink transmission probability constraint. The following proposition simplifies the successful downlink transmission probability constraint and obtains an equivalent constraint of ρ(t). Proposition 4.2: When α = 4, for each horizon t, there exists a minimum required BSs’ active operation √ g0 βσ2 1−ǫ v , probability ρmin (t) = max λλmB(t)B as a non-decreasing function of time-varying MT density √ 3 W ǫ π 2 λB PB 2 g0 g0 , such λm (t), where v is given in Proposition 4.1 and g0 is the unique solution to g0 Q 2π = (1 − ǫ) exp − 4π that Psuc (t) ≥ 1 − ǫ is equivalent to ρ(t) ≥ ρmin (t).
Proposition 4.2 is proved by using the same method as in our previous work [31], which is omitted here for brevity. It is noted that when the MT density and thus the required total bandwidth λm (t)B of the MTs increase, due to the decreased number of bands δ(t) for frequency reuse, as given in (1), the downlink interference level increases. In this case, it is observed from Proposition 4.2 that the minimum required BSs’ active operation probability ρmin (t) also increases, such that more BSs need to be active to increase the desired signal strength at the MT, so as to assure a high Psuc (t). Due to the variations of MT density λm (t), it is also observed that ρmin (t) is also varying over time in general. Moreover, similar to [31], the noise power σ 2 6= 0 in the expression
of ρmin (t) provides a valid minimum required active operation probability for each BS, which is important to assure a sufficiently large Psuc in a noise-dominant network. Therefore, for the ease of analysis, we focus on the case of α = 4 in the sequel.2 By applying Proposition 4.2, we find an equivalent problem to problem (P1) as follows (P2) :
min.
ρ(t),G(t)
T X
a(t)G(t)
(15)
t=1
s.t. G(t) ≥ max(E(t) − B(t) − λe (t), 0), ∀t ∈ {1, ..., T }, ρmin (t) ≤ ρ(t) ≤ 1, ∀t ∈ {1, ..., T },
(16) (17)
B(t + 1) = min (B(t) + λe (t) + G(t) − E(t), C) , ∀t ∈ {1, ..., T − 1},
where the new constraint (16) is obtained by combining the energy demand constraint B(t) + G(t) + λe (t) ≥ E(t) and the non-zero on-grid energy constraint G(t) ≥ 0 in problem (P1), and the other new constraint (17) is obtained by combining the constraints (12) and (14) in problem (P1), as well as Proposition 4.2. To avoid the trivial case where the load is too high to meet even when all BSs are activated, we assume ρmin (t) ≤ 1 in the sequel. In the 2
The value of α does not affect the main results of this paper. Moreover, for other cases with α 6= 4, it is easy to verify from (10) that
there always exist feasible regions for ρ(t) such that Psuc ≥ 1 − ǫ is guaranteed, though it is more difficult to exactly define these feasible regions. Still, our later analysis for studying the on-grid energy cost minimization problem applies.
14
next section, we focus on solving problem (P2) by using the DP technique. V. P OLICY D ESIGN
FOR
O N - GRID E NERGY C OST M INIMIZATION
This section studies the optimal on-grid energy cost minimization policy for solving problem (P2). The optimal on-grid energy minimization policy consists of two parts. One is the optimal BS on/off policy, which decides the optimal BSs’ active operation probability ρ(t) in each horizon t. The other is the optimal on-grid energy purchase policy, which decides the amount of on-grid energy G(t) to purchase in each horizon t. In the following, by studying the cost-to-go function that is defined based on the Bellman’s equations [26], we first investigate the optimal BS on/off policy, and then focus on designing the optimal on-grid energy purchase policy. A. Cost-To-Go Function In this subsection, we define the cost-to-go function for problem (P2). Due to the fluctuations of the on-grid energy prices, the MT density, as well as the renewable energy arrival rates over time, the decision of G(t) and ρ(t) in each horizon t are related to both current and future values of the on-grid energy prices, the MT density,
and the renewable energy arrival rates, as well as the current storage level B(t) at the CES. We thus define the system state in horizon t as Λ(t) = (B(t), (λe (t), ..., λe (T )) , (a(t), ..., a(T )) , (λm (t), ..., λm (T ))). Based on the Bellman’s equations [26], for a given system state Λ(t) in horizon t, we then define the cost-to-go function Jt (Λ(t)) for problem (P2) as min ρmin (T )≤ρ(T )≤1, a(T )G(T ), if t = T, G(T )≥max(E(T )−B(T )−λe (T ),0) Jt (Λ(t)) = min ρmin (t)≤ρ(t)≤1, a(t)G(t) + Jt+1 (Λ(t + 1)|Λ(t)), if 1 ≤ t < T.
(18)
G(t)≥max(E(t)−B(t)−λe (t),0)
When 1 ≤ t < T , the first term a(t)G(t) in (18) represents the instantaneous on-grid energy cost in the current
horizon t, and the second term Jt+1 (Λ(t + 1)|Λ(t)) gives the minimized on-grid energy cost over all the future horizons from t + 1 to T , where for a given Λ(t) and thus B(t), B(t + 1) in Λ(t + 1) is updated according to (8). When t = T , only the instantaneous on-grid energy cost is considered in (18). By finding the cost-to-go function Jt (Λ(t)) in each horizon t, one can equivalently obtain the optimal BS on/off policy and the optimal on-grid
energy purchase policy, which gives the optimal ρ∗ (t) and G∗ (t) in each horizon t for problem (P2), respectively. Moreover, unlike [37], due to the update of the storage level in (8) that correlates the optimal G∗ (t) and ρ∗ (t) over time, a joint design of the optimal BS on/off policy and the optimal on-grid energy purchase policy is generally required. In the next two subsections, we first study the optimal BS on/off policy, under a given on-grid energy purchase policy. Then by substituting the obtained optimal BS on/off policy into the cost-to-go function (18), we study the optimal on-grid energy purchase policy.
15
B. Optimal BS On/Off Policy This subsection studies the optimal BS on/off policy for problem (P2), by supposing that an arbitrary on-grid energy purchase policy is given, i.e., the microgrid knows the amount of on-grid energy G(t) to purchase for a given system state Λ(t) in each horizon t. We first study the impact of storage level B(t) on the cost-to-go function Jt (Λ(t)) in the following lemma. Lemma 5.1: The cost-to-go function Jt (Λ(t)) in horizon t is non-increasing over the current storage level B(t). Proof: Lemma 5.1 is proved based on the mathematical induction method. First, at the last horizon t = T , the optimal G∗ (T ) is given as G∗ (T ) = max(E(T ) − B(T ) − λe (T ), 0). As a result, G∗ (T ) and thus JT (Λ(T )) = a(T )G∗ (T ) is non-increasing over B(T ). Next, suppose at horizon t + 1, Jt+1 (Λ(t + 1)) is non-increasing over B(t + 1). Then, at horizon t, when B(t) increases, since B(t + 1) is non-decreasing from (8), the future cost Jt+1 (Λ(t + 1)|Λ(t)) is non-increasing. Moreover, similar to the case in horizon T , the current cost a(t)G(t) is
also non-increasing over B(t). Therefore, by combing the current and future cost, we find that Jt (Λ(t)) is also non-increasing over B(t). Lemma 5.1 thus follows. Next, by studying the impact of ρ(t) in the current horizon on the storage level B(t + 1) in the next horizon, we obtain the following lemma. Lemma 5.2: For any given on-grid energy purchase policy, the cost-to-go function Jt (Λ(t)) in horizon t is non-decreasing over ρ(t). The proof of Lemma 5.2 can be easily obtained by applying Lemma 5.1 and the fact that for any given storage level B(t) and purchase amount G(t) in the current horizon, the storage level B(t + 1) in the next horizon is non-increasing over ρ(t) in the current horizon according to (6) and (8), and thus is omitted here for brevity. Finally, as a straightforward result from Lemma 5.2, to minimize the overall on-grid energy cost in all T horizons, we obtain the optimal BS on/off policy as follows. Proposition 5.1 (Optimal BS on/off policy): The optimal ρ∗ (t) in each horizon t for problem (P2) is given by ρ∗ (t) = ρmin (t), ∀t ∈ {1, ..., T },
(19)
where ρmin (t) is given in Proposition 4.2 and depends on time-varying MT density λm (t). Remark 5.1: From Proposition 5.1, since a small ρ(t) yields a low energy demand E(t) of the BSs, and thus a low on-grid energy cost in general, the optimal ρ∗ (t) in each horizon t is always selected as the minimum required ρmin (t) which is just sufficient to assure the required downlink transmission quality. Based on ρ∗ (t), each BS
then independently decides its on/off status in each horizon t. It is worth noting that since the traffic loads of the inactive BSs are handed over to their neighboring active BSs, as described in Section II, although the exact
16
on/off status of the BSs are determined independently, the network-level downlink transmission quality can still be assured with ρ∗ (t) = ρmin (t). It is also noted from Proposition 4.2 that the optimal ρ∗ (t) = ρmin (t) generally varies over time to adapt to the variations of the MT density λm (t). C. Optimal On-Grid Energy Purchase Policy This subsection studies the on-grid energy purchase policy for problem (P2). By substituting the optimal ρ∗ (t) into (18), the cost-to-go function for problem (P2) is simplified as ( minG(T )≥max(Emin (T )−B(T )−λe (T ),0) a(T )G(T ), if t = T, Jt (Λ(t)) = minG(t)≥max(Emin (t)−B(t)−λe (t),0) a(t)G(t) + Jt+1 (Λ(t + 1)|Λ(t)), if 1 ≤ t < T.
(20)
where by substituting ρ∗ (t) = ρmin (t) into (6), we obtain the minimum value of E(t), denoted by Emin (t) = E(t)|ρ(t)=ρmin (t) , as Emin (t) = λB ρmin (t)Pgap + λB Ps +
PB λm (t) , µ
(21)
and accordingly, the storage level B(t + 1) in system state Λ(t + 1) is updated from B(t) in system state Λ(t) as B(t + 1) = min (B(t) + λe (t) + G(t) − Emin (t), C) .
(22)
However, unlike the design of optimal BS on/off policy, due to the fluctuations of the on-grid energy price a(t) over time, the cost-to-go function Jt (Λ(t)) in (20) is generally neither non-increasing nor non-decreasing over G(t) in each horizon t. Although based on (20), one can always find the optimal G∗ (t), the complexity of finding G∗ (t) in each horizon increases exponentially over the total number of horizons T , due to the curse of dimensionality for the DP problem. As a result, in the following, to obtain tractable and insightful on-grid energy purchase policy for problem (P2), we focus on designing a suboptimal on-grid energy purchase policy, by studying the optimal solutions of G∗ (t) for the last horizon T and the horizon T − 1. We also assume the storage capacity C is sufficiently large with C ≥ max (Emin (1), ..., Emin (T )), such that the BSs’ demanded energy for each horizon can be properly stored. First, we look at the last horizon with t = T , and assume the system state Λ(T ) = (B(T ), λe (T ), a(T ), λm (T )) in horizon T is given. From (20), it is easy to find that the optimal G∗ (T ) in the last horizon T is given by G∗ (T ) = max(Emin (T ) − B(T ) − λe (T ), 0).
(23)
It is also noted that the optimal G∗ (T ) is a myopic solution, which minimizes the instantaneous on-grid energy cost in the current horizon. Moreover, when G∗ (T ) is myopic, it is noted that the optimal ρ∗ (T ) that minimizes E(T ) and thus the instantaneous on-grid energy cost a(T )G(T ) is also a myopic solution.
17
Next, we consider the horizon with t = T − 1, and assume the system state in horizon T − 1, i.e., Λ(T − 1) = (B(T − 1), (λe (T − 1), λe (T )) , (a(T − 1), a(T )) , (λm (T − 1), λm (T ))), is given. We also assume the harvested
renewable energy in horizon T − 1 is reasonably small with λe (T − 1) < C + Emin (T − 1) − B(T − 1), as compared to the storage capacity and the energy demand of all the BSs. Based on (20), we jointly consider the impact of G(T − 1) on the instantaneous on-grid energy cost in the current horizon T − 1 and the future on-grid energy cost JT (Λ(T )|Λ(T − 1)), where the storage level B(T ) in the next horizon T is updated from B(T − 1) according to (22), and obtain the optimal G∗ (T − 1) in the following proposition. Proposition 5.2: Given the system state Λ(T − 1) in horizon T − 1, the optimal G∗ (T − 1) for problem (P2) is given as follows: •
Large current storage regime with B(T − 1) ≥ Emin (T ) + Emin (T − 1) − λe (T − 1) − λe (T ): we have G∗ (T − 1) = max(Emin (T − 1) − B(T − 1) − λe (T − 1), 0).
•
(24)
Low current storage regime with B(T − 1) < Emin (T ) + Emin (T − 1) − λe (T − 1) − λe (T ): – High future renewable energy regime with λe (T ) ≥ Emin (T ): G∗ (T ) is also given by (24). – Low future renewable energy regime with λe (T ) < Emin (T ): we have ( max(Emin (T − 1) − B(T − 1) − λe (T − 1), 0), if a(T − 1) ≥ a(T ), ∗ G (T − 1) = Emin (T ) + Emin (T − 1) − λe (T − 1) − λe (T ) − B(T − 1), if a(T − 1) < a(T ). (25) Proof: Please refer to Appendix A.
Observation 1: From Proposition 5.2, the decision of G∗ (T − 1) is jointly determined by the fluctuations of MT density, on-grid energy prices, the renewable energy arrival rates in both current and future horizons. Moreover, the on-grid energy G∗ (T − 1) is over-purchased only when all the following three conditions are satisfied: 1) The current storage level B(T − 1) is within the low current storage regime; 2) The future renewable energy λe (T ) alone cannot to satisfy the demanded energy Emin (T ) in the future; 3) The current on-grid energy price is lower than the future with a(T − 1) < a(T ). When all the three conditions are satisfied, the microgrid purchases more energy than the demanded in the current horizon, so as to ensure both current and future energy demands can be satisfied. In addition, it is also observed that if the on-grid energy prices a(T − 1) = a(T ), G∗ (T − 1) is not affected by the fluctuations of the on-grid energy prices and becomes a myopic solution in all cases of Proposition 5.2, just as G∗ (T ) in (23). Finally, we extend the above observations for the optimal G∗ (t) obtained in horizon T and horizon T − 1 to the on-grid energy purchase decision in an arbitrary horizon t ∈ {1, ..., T }, and propose a suboptimal on-grid energy purchase policy for problem (P2) in the following, by jointly considering the impact of storage capacity C .
18
Proposition 5.3 (Suboptimal on-grid energy purchase policy): Denote G⋆ (t) as the suboptimal solution of the purchased on-grid energy to problem (P2) in each horizon t ∈ {1, ..., T }. The microgrid over-purchases the on-grid energy in horizon t if all of the three following conditions are satisfied at the same time: P P 1) Low current storage level: min (B(t) − Emin (t) + λe (t), C) < Tk=t+1 Emin (k) − Tk=t+1 λe (k); P P 2) Low future renewable energy: Tk=t+1 Emin (k) − Tk=t+1 λe (k) ≥ 0;
3) Low current on-grid energy price: a(t) < min (a(t + 1), ..., a(T )). P PT T ⋆ Let ω(t) = min k=t+1 Emin (k) − k=t+1 λe (k), C . Given the system state Λ(t) in horizon t, G (t) is
given by
G⋆ (t) =
(
max (ω(t) + Emin (t) − λe (t) − B(t), 0) , if conditions 1)-3) hold,
max(Emin (t) − B(t) − λe (t), 0),
otherwise.
(26)
From Proposition 5.3, the microgrid over-purchases the low-price energy in the current horizon for future usage only when the current storage and future renewable supply are both low. It is also easy to verify that G⋆ (t) = G∗ (t) when t = T or t = T − 1. Moreover, when a(t) = a(t′ ), ∀t, t′ ∈ {1, ..., T } and t 6= t′ , G⋆ (t) becomes the myopic solution max(Emin (t)− B(t)− λe (t), 0) for all horizons, which is similar to the case with t = T − 1 in Proposition 5.2. Moreover, it is also easy to find that as compared to the DP-based approach to search the optimal G∗ (t) in each horizon by applying (20), the proposed suboptimal on-grid energy purchase policy is of low complexity to implement. We will validate the performance of the proposed suboptimal on-grid energy purchase policy for multi-horizon in Section VI. To conclude, by combining the optimal suboptimal BS on/off policy and the suboptimal on-grid energy purchase policy, given in Proposition 5.1 and Proposition 5.3, respectively, the on-grid energy cost minimization problem (P2) is solved. VI. S IMULATION R ESULTS This section studies the proposed network-level on-grid energy cost minimization policy by providing extensive numerical and simulation results. We first illustrate the adaptation of ρ(t) and G(t) over horizons by applying the proposed policy. We then validate the performance of the proposed suboptimal on-grid energy purchase policy. After that, we study the effects of the prediction errors in terms of the network profiles λe (t)’s, λm (t)’s, and a(t)’s. At last, we conduct a large-scale network-level simulation to compare the performance of the proposed policy with two reference schemes. Without specified otherwise, we consider the following simulation parameters in this section. We set the basic energy consumptions of each BS in the active and inactive modes as Pa = 130W and Ps = 75W, respectively, the
19
Active MT ratio
θ(t)
0.18 0.12 0.06
2
λe(t) (W/m )
0
0.25 0.2 0.15 0.1 0.05 0
2
4
6
8
a(t) ($/W)
Fig. 3.
12 14 Hour
16
18
20
22
24
18
20
22
24
18
20
22
24
Solar energy arrival rate
2
4
6
8
−5
6 5 4 3 2 1 0
10
12 14 Hour
16
On−grid energy price
x 10
2
10
4
6
8
10
12 14 Hour
16
Network profiles on the variations of λm (t), λe (t), and a(t) for a single day.
transmit power of each active BS as PB = 20W , and the power amplifier parameter as µ = 21.3% [36]. We also set the noise power σ 2 = 10−9 W, the targeted SINR level β = 2, the path-loss exponent α = 4, and the maximum outage ǫ = 0.05 to ensure a sufficiently high downlink successful transmission probability Psuc (t) in each horizon. According to the real data measured in [36], [38], and [39], we consider the variations of the MT density λm (t), the solar energy arrival rate λe (t), as well as the on-grid energy price a(t) over a single day, and show them in Fig. 3. In particular, instead of λm (t), Fig. 3 shows the active MT ratio θ(t) for each horizon t. By denoting the all total MT density in the network as λall m > 0, we can obtain the MT density in each horizon as λm (t) = λm θ(t). −3 2 −4 2 We set λall m = 8 × 10 /m , and the BS density as λB = 5 × 10 /m . We also set the storage capacity C at the
CES, given in (22), as C = 0.2 (W/m2 ). A. Illustration for Dynamic Network Operation In this subsection, we illustrate the dynamic adaptation of the BSs’ active operation probability and the amount of on-grid energy to purchase over horizons, by applying the optimal BS on/off policy and the suboptimal on-grid energy purchase policy, given in Proposition 5.1 and Proposition 5.3, respectively. By applying the network profiles given in Fig. 3, Fig. 4 shows the variations of the storage level B(t), as well as the adaptations of the BSs’ active operation probability and the amount of purchased on-grid energy to these network profiles over T = 24 hours. From Proposition 5.1, the optimal ρ∗ (t) = ρmin (t) in each horizon, where ρmin (t) given in Proposition 4.2 varies monotonically over the MT density λe (t). It is thus observed from Fig. 4
20
Storage level at the CES 2
B(t) (W/m )
0.2 0.15 0.1 0.05 0
2
4
6
8
10
2
4
2
4
12 14 16 18 Hour BSs’ active operation probability
6
8
10
6
8
20
22
24
18
20
22
24
18
20
22
24
0.8 ρ (t)
0.6 0.4 0.2
2
G (t) (W/m )
0
Fig. 4.
0.3 0.24 0.18 0.12 0.06 0
12 14 16 Hour On−grid energy purchase
10
12 14 Hour
16
Adaptations of ρ∗ (t) and G⋆ (t) and the variation of B(t) for a single day.
that ρ∗ (t) follows the exact variation trend of λe (t) in Fig. 3. It is also observed that when the solar energy arrival rate is substantially large from t = 7 to t = 19 in Fig. 3, the storage level B(t) in Fig. 4 is also quite high, and achieves the capacity of the CES in most horizons from t = 7 to t = 19. However, when the solar energy arrival rate is approaching to zero from t = 1 to t = 6 and from t = 20 to t = 24 in Fig. 3, the storage level B(t) in Fig. 4 is also almost zero in each of these horizons. Moreover, from t = 1 to t = 5, since the storage level is almost zero in Fig. 4, but the MT density is not low and the on-grid energy price is non-decreasing in Fig. 3, it is observed from Fig. 4 that the purchased on-grid energy follows the variation trends of the MT density and thus is myopic solution in each of these horizons. At horizon t = 6, it is observed from Fig. 3 that a(6) at t = 6 is the lowest value among all the on-grid energy prices over T horizons. Thus, from the suboptimal on-grid energy purchase policy given in Proposition 5.3, it is noted that the purchased on-grid energy G(6) at t = 6 achieves the highest value among all the T horizons. From t = 7 to t = 19, due to the large storage level B(t), the purchased on-grid energy is decreasing in these horizons and becomes zero from t = 10 to t = 18 in Fig. 4. At last, from t = 19 to t = 24, due to the almost zero storage level in Fig. 4 as well as the large MT density in Fig. 3, the
purchased on-grid energy is largely increased in these horizons in Fig. 4.
21
B. Suboptimal Policy Performance Validation
C OMPARISON OF
TABLE I PT TOTAL ON - GRID ENERGY COST t=1 a(t)G(t) UNDER OPTIMAL AND
Number of horizons −6
Optimal (×10
2
$/m )
Suboptimal (×10−6 $/m2 ) −6
Absolute error (×10
2
$/m )
SUBOPTIMAL POLICES .
T =1
T =2
T =3
T =4
T =5
2.1004935
3.6441487
5.005297
6.1624714
7.0429205
2.1004935
3.6441487
5.005329
6.1624757
7.0429250
4.3 × 10
4.5 × 10−5
0
0
3.2 × 10
−5
−5
Since the proposed BS on/off policy in Proposition 5.1 is optimal, the suboptimality of the proposed on-grid energy minimization policy is caused by the suboptimality of the proposed on-grid energy purchase policy in Proposition 5.3. Thus, by applying the optimal BS on/off policy, we compare the performance of the suboptimal on-grid energy purchase policy in Proposition 5.3 with the optimal on-grid energy purchase policy that is obtained by exhaustive search. To search the optimal G∗ (t), we first find the search range of G(t) in each horizon t. From (20), the minimum required G(t) is given by Gmin (t) = max(Emin (t)−B(t)−λe (t), 0). It is also easy to obtain that P P to save the on-grid energy cost, the maximum required G(t) is given by Gmax (t) = Tk=t Emin (k) − Tk=t λe (k),
which can just satisfy the energy demand of all the BSs from the current to all future horizons. Then, we apply the DP technique from [26] to search the optimal G∗ (t) ∈ [Gmin (t), Gmax (t)] in each horizon t, where the range of [Gmin (t), Gmax (t)] is quantized with a sufficiently small step-size 1.0 × 10−8 . In Table I, we compare PT PT ∗ ⋆ t=1 a(t)G (t) that is obtained by the optimal search with t=1 a(t)G (t) that is obtained by applying the
proposed suboptimal on-grid energy purchase policy in Proposition 5.3. Due to the curse of dimensionality to apply the DP technique for searching the optimal G∗ (t), we only consider the total number of horizons T ∈ {1, .., 5} in Table I, using network profiles starting from t = 2 in Fig. 3 as an example. It is observed from Table I that
both the optimal and suboptimal total on-grid energy cost increases over T . It is also observed that as the number of horizons T increases, since the number of optimal G∗ (t) and suboptimal G⋆ (t) with different values in each P P horizon t also increases in general, the absolute error, given by Tt=1 a(t)G⋆ (t) − Tt=1 a(t)G∗ (t), increases over T accordingly. However, it is observed that the absolute error at each T is quite small and is within 4.5 × 10−11 .
As a result, the proposed suboptimal on-grid energy purchase policy in Proposition 5.3 and thus the proposed suboptimal on-grid energy minimization policy in Section V can achieve near-optimal performance. C. Effects of Prediction Errors In this subsection, by applying our proposed on-grid energy cost minimization policy in Section V, we study the effects of the prediction errors for λe (t)’s, λm (t)’s, and a(t)’s. In particular, we focus on the effects of the prediction errors for the solar energy arrival rate λe (t)’s, and assume λm (t)’s and a(t)’s in all horizons are still
22
−5
5
x 10
4.5 W/ prediction error W/o prediction error
2
Toal on−grid energy cost ($/m )
4 3.5 3 2.5 2 1.5 1 0.5 0
Fig. 5.
2
4
6
8
10
12 14 Hour
16
18
20
22
24
Effects of the predictions errors on the network performance.
perfectly known. We also observe similar simulation results for the cases where prediction errors exist in λm (t) and/or a(t), which are thus omitted here for brevity. Specifically, we assume a prediction error ∆(t) exists in each horizon t. We also assume ∆(t)’s are identical and independently distributed zero-mean Gaussian variables over different horizons, where the standard deviation is denoted by η . Since the prediction error ∆(t) can be large due to the large T = 24 under consideration, based on the value of λe (t) ∈ [0.01 W/m2 , 0.22 W/m2 ] shown in Fig. 3, we reasonably consider a large standard deviation for the prediction errors with η = 0.1 in each horizon. As a result, under the consideration of prediction errors, the microgrid adopts λe (t) + ∆(t) as the solar energy arrival rate in each horizon t to decide ρ∗ (t) and G⋆ (t). However, since the storage level B(t) in each horizon can be perfectly known by the microgrid, the storage level B(t + 1) in horizon t + 1 is still updated according to the actual solar energy arrival rate λe (t) in horizon t. P Fig. 5 compares the total on-grid energy cost Tt=1 a(t)G(t) over the number of horizons T under both cases
with and without the prediction errors. For the case with prediction errors, at each T in Fig. 5, we consider 1000 realizations of the T horizons with predictions errors, and take the average on-grid energy cost in T horizons over the 1000 realizations as the total on-grid energy cost for comparison. It is observed from Fig. 5 that since we focus on the average energy consumption, as discussed in III-B, the network performance under the case with
the prediction errors is always very close to that under the case without errors at each T . Hence, our proposed policy is robust to the prediction errors. It is also observed from Fig. 5 that the total on-grid energy cost generally
23
increases with the number of horizons T , as expected. D. Geo-DP based Scheme versus Reference Schemes To show the performance of our proposed scheme based on the Geo-DP approach, where the BSs are coordinated from a network-level to serve the MTs, this subsection conducts Monte Carlo simulations to compare the proposed scheme with two reference schemes, which are specified as follows. First, we consider a simple reference scheme without BSs’ mutual coordination. In this scheme, each BS individually decides whether to be active or not in each horizon: if there exist MTs locating in a BS’s coverage, it keeps active, or becomes inactive, otherwise. We assume the BS density λB is sufficiently large, such that when all the BSs are active, the downlink transmission quality is assured in the network. Thus, since only BSs with no associated MTs turn to be inactive, the downlink transmission quality can be assured in this scheme. The total energy demand of the BSs is calculated according to (3) and (4), depending on whether a BS is active or not, respectively. Second, we consider a reference scheme with cluster-based BS coordination, which has been proposed in the existing literature to coordinate the BSs’ on/off operations in a large-scale network (see, e.g., [24] and [25]). However, it is noted that due to the considered different performance metrics, the existing cluster-based BS coordination policy cannot be applied directly. For example, in [24] and [25], the authors considered a blocking probability constraint for the new MT’s access as the QoS constraint, while we use the SINR-based downlink transmission quality as the QoS constraint. As a result, we consider the following cluster-based BS coordination scheme for comparison. Specifically, we assume each cluster includes at most two BSs, and any two BSs that are located within L > 0 meters form a cluster. Each BS can belong to at most one cluster. If a BS does not belong to any cluster, it decides its on/off status as that in the scheme with no BS coordination. For any two BSs that belong to one cluster, the BS that has less associated MTs becomes inactive, while the other becomes active and also helps to serve the MTs of the inactive BS. The total energy consumption of all the BSs is also calculated as the same as that in the scheme with no BS coordination. It is easy to verify that the downlink transmission quality can be assured with a small L, where the MTs are always connected to the active BSs that are located closely for good communication quality. To obtain the total on-grid energy cost under our proposed scheme and the two reference schemes, we conduct a large-scale network simulation. We generate Φ(λB ) for BSs and Φ(λm (t)) for MTs in a square of [0m, 1000m]× [0m, 1000m], according to the method described in [12]. Due to the computation complexity that
increases substantially over both network size as well as the number of horizons T , we set T = 3 and the number
24
4
9
x 10
Total on−grid energy cost ($)
8.5
Proposed Geo−DP Cluster−based No coordination
8
7.5
7
6.5
6
5.5 2
2.5
3 3.5 4 all 2 Total MT density λ (/m ) m
Fig. 6.
4.5
5 −3
x 10
Network performance comparison with two reference schemes.
of slots within each horizon as N = 1 in this simulation, and take network profiles of λe (t) and θ(t) from Fig. 3 for t ∈ {7, 8, 9}. For the cluster-based scheme, we use a small L with L = 100m to assure the downlink transmission
quality. We consider in total 20000 network realizations, and take the average value over the 20000 realizations as the network performance for each considered scheme. For simplicity, we assume identical on-grid energy price with a(t) = 1 ($/W) in each horizon t. In this case, it is easy to find that the on-grid energy purchase decision for all three considered schemes is always myopic in each horizon from Section V. Similar performance can also be observed for all three schemes with non-identical on-grid energy prices, which are thus omitted here for brevity. Fig. 6 compares the total on-grid energy cost across T horizons over the total MT density λall m . It is observed all from Fig. 6 that as the total MT density λall m increases, since the actual MT density λm (t) = λm θ(t) also increases,
the total on-grid energy cost of all three schemes increases accordingly, so as to properly serve the increased MTs in the network. It is also observed that at each value of λall m , the total on-grid energy cost under the proposed network-level BS coordination using the Geo-DP approach is always the minimum, that under the cluster-based scheme is always the medium, and that under the scheme with no BS coordination is always the largest. Moreover, it is also worth noting that although we only show the performance of the cluster-based scheme with two BSs in one cluster, when the number of BSs in each cluster increases, the performance of the cluster-based scheme is expected to gradually approach our proposed network-level BS coordination scheme. However, unlike our proposed closed-
25
form optimal BS on/off operation policy, which has low-complexity to coordinate all the BSs’ on/off operations in the network, the complexity of coordinating the increased number of BSs in each cluster can be largely increased for the cluster-based scheme to assure the downlink transmission quality [24], [25]. As a result, our proposed network-level BS coordination scheme based on the Geo-DP approach can more efficiently solve the on-grid energy cost minimization problem for a large-scale network. VII. C ONCLUSION In this paper, we considered a large-scale green cellular network, where the BSs are aggregated as a microgrid with hybrid energy supplies and an associated CES. We proposed a new Geo-DP approach to jointly apply stochastic geometry and the DP technique to conduct network-level dynamic system design. By optimally selecting the BSs’ active operation probability as well as the amount of on-grid energy to purchase in each horizon, we studied the on-grid energy cost minimization problem. We first studied the optimal BS on/off policy, which shows that the optimal BSs’ active operation probability in each horizon is just sufficient to assure the required downlink transmission quality with time-varying load in a large-scale network. We then proposed a suboptimal on-grid energy purchase policy, where the low on-grid energy is over-purchased in the current horizon only when both the current storage level and future renewable energy level are both low. We also compared our proposed policy with the existing schemes, and showed that our proposed policy can efficiently save the on-grid energy cost in the large-scale network. It is also of our interest to minimize the network-level on-grid energy cost by jointly considering the uplink and downlink traffic loads in our future work. A PPENDIX A P ROOF
TO
P ROPOSITION 5.2
According to (20), to find the optimal G∗ (T − 1), we need to first find the optimal G∗ (T ) in horizon T to obtain the exact expression of JT (Λ(T )|Λ(T − 1)) for a given system state Λ(T − 1) in horizon T − 1. From (22), for a given B(T − 1) in horizon T − 1, the storage level of B(T ) is obtained as B(T ) = min (B(T − 1) + λe (T − 1) + G(T − 1) − Emin (T − 1), C) .
(27)
From (23), it is observed that when G(T − 1) = 0, if B(T ) in (27) can be assured to be no smaller than Emin (T ) − λe (T ), we always have G∗ (T ) = 0. Moreover, under the assumption with λe (T − 1) < C + Emin (T − 1) − B(T − 1) and C ≥ max (Emin (1), ..., Emin (T )), it is easy to verify that (27) can be equivalently reduced to B(T ) = B(T − 1) + λe (T − 1) + G(T − 1) − Emin (T − 1), and thus B(T ) ≥ Emin (T ) − λe (T ) is equivalent to B(T −1) ≥ Emin (T )+Emin (T −1)−λe (T −1)−λe (T ). We thus consider two cases to find the optimal G∗ (T −1).
26
One is the case with large current storage regime B(T − 1) ≥ Emin (T ) + Emin (T − 1) − λe (T − 1) − λe (T ), and the other is the case with small current storage regime B(T − 1) < Emin (T ) + Emin (T − 1) − λe (T − 1) − λe (T ). First, for case with large current storage regime, since G∗ (T − 1) = 0 and thus JT (Λ(T )|Λ(T − 1)) = 0, according to (20), we obtain JT −1 (Λ(T − 1)) = minG(T −1)≥max(Emin (T −1)−B(T −1)−λe (T −1),0) a(T − 1)G(T − 1). As a result, the optimal G∗ (T − 1) in this case is the myopic solution with G∗ (T − 1) = max(Emin (T − 1) − B(T − 1) − λe (T − 1), 0).
Next, for the case with small current storage regime, since G∗ (T ) is not assured to be 0 or Emin (T ) − B(T ) − λe (T ), B(T ) cannot be reduced as in the previous case without considering the storage capacity C . As this case
is more complicated than the previous one, based on the possible values of G∗ (T ), we consider and compare the following two options. 1) Option 1: In this option, we assume G∗ (T ) = 0. According to (23) and (27), G∗ (T ) = 0 holds if and only if min (B(T − 1) + λe (T − 1) + G(T − 1) − Emin (T − 1), C) ≥ Emin (T ) − λe (T ). Then, according to (20), we solve the following optimization problem to find G∗ (T − 1). (P3) :
min.
G(T −1)
a(T − 1)G(T − 1)
s.t. G(T − 1) ≥ max(Emin (T − 1) − B(T − 1) − λe (T − 1), 0), min (B(T − 1) + λe (T − 1) + G(T − 1) − Emin (T − 1), C) ≥ Emin (T ) − λe (T ),
(28)
where G∗ (T − 1) is given as follows: •
Case with λe (T ) ≥ Emin (T ): It is easy to verify that (28) always holds when G(T − 1) = max(Emin (T − 1)−B(T −1)−λe (T −1), 0) under the small current storage regime. As a result, G∗ (T −1) = max(Emin (T − 1) − B(T − 1) − λe (T − 1), 0).
•
Case with λe (T ) < Emin (T ): (28) holds in this case under two scenarios. One is that G(T − 1) ≥ C + Emin (T −1)−λe (T −1)−B(T −1), such that min (B(T − 1) + λe (T − 1) + G(T − 1) − Emin (T − 1), C) = C > Emin (T ) − λe (T ). The other is that G(T − 1) < C + Emin (T − 1) − λe (T − 1) − B(T − 1), such that min (B(T − 1) + λe (T − 1) + G(T − 1) − Emin (T − 1), C) = B(T −1)+λe (T −1)+G(T −1)−Emin (T − 1), which renders (28) to be equivalent to G(T −1) ≥ Emin (T )+Emin (T −1)−λe (T −1)−λe (T )−B(T −1).
It is easy to verify that the second scenario yields a smaller G∗ (T − 1) = Emin (T ) + Emin (T − 1) − λe (T − 1) − λe (T ) − B(T − 1). As a result, we obtain the cost-to-go function in this case as JT −1 (Λ(T − 1)) = a(T − 1) (Emin (T ) + Emin (T − 1) − λe (T − 1) − λe (T ) − B(T − 1)).
2) Option 2: In this option, we assume G∗ (T ) = Emin (T )−B(T )−λe (T ). According to (23) and (27), G∗ (T ) = Emin (T ) − B(T ) − λe (T ) holds if and only if min (B(T − 1) + λe (T − 1) + G(T − 1) − Emin (T − 1), C)
Emin (T ) − λe (T ). We let U = a(T ) (Emin (T ) + Emin (T − 1) − λe (T − 1) − λe (T ) − B(T − 1)). According to (20) and (27), we solve the following optimization problem to find G∗ (T − 1) in this option. (P4) :
min.
G(T −1)
(a(T − 1) − a(T ))G(T − 1) + U
s.t. G(T − 1) ≥ max(Emin (T − 1) − B(T − 1) − λe (T − 1), 0), B(T − 1) + λe (T − 1) + G(T − 1) − Emin (T − 1) < Emin (T ) − λe (T ).
(29) (30)
It is easy to find G∗ (T − 1) as follows: •
Case with λe (T ) ≥ Emin (T ): It is observed that if Emin (T − 1) − B(T − 1) − λe (T − 1) ≥ 0, (29) is reduced to G(T − 1) ≥ Emin (T − 1) − B(T − 1) − λe (T − 1). However, (30) requires G(T − 1) < Emin (T ) − λe (T ) + Emin (T − 1) − B(T − 1) − λe (T − 1), where the right-hand side is smaller than Emin (T −1)−B(T −1)−λe (T −1). Similar case is also observed when Emin (T −1)−B(T −1)−λe (T −1) < 0.
As a result, no feasible solution of G(T − 1) exists in this case for Option 2, which implies that under this case the optimal G∗ (T − 1) is obtained by using Option 1. •
Case with λe (T ) < Emin (T ): Clearly, if a(T −1) ≥ a(T ), G∗ (T −1) = max(Emin (T −1)−B(T −1)−λe (T − 1), 0), with the cost-to-go function given by JT −1 (Λ(T − 1)) = (a(T − 1)− a(T )) max(Emin (T − 1)− B(T − 1) − λe (T − 1), 0) + a(T ) (Emin (T ) + Emin (T − 1) − λe (T − 1) − λe (T ) − B(T − 1)). If a(T − 1) < a(T ), G∗ (T − 1) = Emin (T ) + Emin (T − 1) − λe (T − 1) − λe (T ) − B(T − 1), with the cost-to-go function given
by JT −1 (Λ(T − 1)) = a(T − 1) (Emin (T ) + Emin (T − 1) − λe (T − 1) − λe (T ) − B(T − 1)). Finally, by comparing the cost-to-go functions that are obtained under both options, it is easy to verify that the cost-to-go functions that are obtained under Option 2 for both cases with a(T − 1) ≥ a(T ) and a(T − 1) < a(T ) are always no larger than that under Option 1. As a result, for the case with small current storage regime, the optimal G∗ (T − 1) is given by that obtained under Option 2. Proposition 5.2 thus follow. R EFERENCES [1] K. Davaslioglu and E. Ayanoglu, “Quantifying potential energy efficiency gain in green cellular wireless networks,” IEEE Commun. Surv. & Tut., vol. 16, no. 4, pp. 2065-2091, 4th Quart., 2014. [2] J. Wu, S. Rangan, and H. Zhang, Green communications: Theoretical fundamentals, algorithms and applications. Boca Raton, FL, USA: CRC Press, 2012. [3] M. A. Marsan et al.,“Optimal energy savings in cellular access networks,” in Proc. IEEE Int. Conf. Commun. (ICC), Jun. 2009. [4] L. M. Correia et al., “Challenges and enabling technologies for energy aware mobile radio networks,” IEEE Commun. Mag., vol. 48, no. 11, pp. 66-72, Nov. 2010. [5] M. Marsan, L. Chiaraviglio, D. Ciullo, and M. Meo, “Optimal energy savings in cellular access networks,” in Proc. IEEE Int. Conf. Commun. Workshops, Jun. 2009.
28
[6] E. Oh, K. Son, and B. Krishnamachari, “Dynamic base station switching-on/off strategies for green cellular networks,” IEEE Trans. Wireless Commun., vol. 12, no. 5, May 2013. [7] S. Luo, R. Zhang, and T. J. Lim, “Optimal power and range adaptation for green broadcasting,” IEEE Trans. Wireless Commun., vol. 12, no. 9, pp. 4592-4603, Sep. 2013. [8] S. Zhou, J. Gong, Z. Yang, Z,. Niu, and P. Yang, “Green mobile access network with dynamic base station energy saving,” in Proc. ACM MobiCom, Sep. 2009. [9] J. Xu, L. Duan, and R. Zhang, “Energy group-buying with loading sharing for green cellular networks,” IEEE J. Sel. Areas Commun., to appear. Available [online] at http://arxiv.org/pdf/1508.06093v1.pdf. [10] C. K. Ho and R. Zhang, “Optimal energy allocation for wireless communications with energy harvesting constraints,” IEEE Trans. Sig. Process., vol. 60, no. 9, pp. 4808-4818, Sep. 2012. [11] H. S. Dhillon, Y. Li, P. Nuggehalli, Z. Pi, and J. G. Andrews, “Fundamentals of heterogeneous cellular networks with energy harvesting,” IEEE Trans. Wireless Commun., vol. 13, no. 5, May 2014. [12] D. Stoyan, W. S. Kendall, and J. Mecke, Stochastic geometry and its applications, 2nd edition. John Wiley and Sons, 1995. [13] D. Gunduz, K. Stamatiou, N. Michelusi and M. Zorzi, “Designing intelligent energy harvesting communication systems,” IEEE Commun. Mag., vol. 52, no. 1, pp. 210-216, Jan. 2014. [14] S. Ulukus, A. Yener, E. Erkip, O. Simeone, M. Zorzi, P. Grover, and K. Huang, “Energy harvesting wireless communications: A review of recent advances,” IEEE J. Sel. Areas Commun., vol. 33, no. 3, pp. 360-381, Mar. 2015. [15] H. Li, J. Xu, R. Zhang, and S. Cui, “A general utility optimization framework for energy harvesting based wireless communications,” IEEE Commun. Mag., vol. 53, no. 4, pp.79-85, April 2015. [16] J. Xu, L. Duan, and R. Zhang, “Cost-aware green cellular networks with energy and communication cooperation,” IEEE Commun. Mag., vol. 53, no. 5, pp. 257-263, May 2015. [17] T. Han and N. Ansari, “On optimizing green energy utilization for cellular networks with hybrid energy supplies,” IEEE Trans. Wireless Commun., vol. 12, no. 8, pp. 3872-3882, Aug. 2013. [18] Y. Mao, G. Yu, and Z. Zhang,“On the optimal transmission policy in hybrid energy supply wireless communication systems,” IEEE Trans. Wireless Commun., vol. 13, no. 1, pp. 6422-6430, Nov. 2014. [19] D. W. K. Ng, E. S. Lo, and R. Schober, “Energy-efficient resource allocation in OFDMA system with hybrid energy harvesting base station,” IEEE Trans. Wireless Commun., vol. 12, no. 7, pp. 3412-3427, Jul. 2013. [20] Y. Mao, J. Zhang, and K. B. Letaief, “A lyapunov optimization approach for green cellular networks with hybrid energy supplies,” IEEE J. Sel. Areas Commun., accepted. Available [online] at http://arxiv.org/pdf/1509.01886v2.pdf. [21] J. Xu and R. Zhang, “CoMP meets smart grid: a new communication and energy cooperation paradigm,” IEEE Trans. Veh. Technol., vol. 64, no. 6, pp. 2476-2488, Jun. 2015. [22] Y. K. Chia, S. Sun, and R. Zhang, “Energy cooperation in cellular networks with renewable powered base stations,” IEEE Trans. Wireless Commun., vol. 13, no. 12, pp. 6996-7010, Dec. 2014. [23] Y. Guo, J. Xu, L. Duan, and R. Zhang, “Joint energy and spectrum cooperation for cellular communication systems,” IEEE Trans. Commun., vol. 62, no. 10, pp. 3678-3691, Oct. 2014. [24] S. Zhou, J. Gong, and Z. Niu, “Sleep control for base stations powered by heterogeneous energy sources,” in Proc. 2013 Intl. Conf. ICT Convergence (ICTC), Oct. 2013. [25] J. Gong, J. S. Thompson, S. Zhou, and Z. Niu, “Base station sleeping and resource allocation in renewable energy powered cellular networks,” IEEE Trans. Commun., vol. 62, no.11, pp. 3801-3813, Nov. 2014. [26] D. Bertsekas, Dynamic Programming and Optimal Control, Athena Scientific, vol. 1, 2005. [27] X. Fang, S. Misra, G. Xue, and D. Yang, “Smart grid-The new and improved power grid: A survey,” IEEE Commun. Surveys Tuts., vol. 14, no. 4, pp. 944-980, 4th Quart., 2012. [28] K. Huang, “Spatial throughput of mobile ad hoc networks with energy harvesting,” IEEE Trans. Inf. Theory, vol. 59, no. 11, pp. 7597-7612, Nov. 2013. [29] S. Lee, R. Zhang, and K. Huang, “Opportunistic wireless energy harvesting in cognitive radio networks,” IEEE Trans. Wireless Commun., vol. 12, no. 9, pp. 4788-4799, Sep. 2013. [30] T. Khan, P. Orlik, and K. J. Kim, “A stochastic geometry analysis of cooperative wireless networks powered by energy harvesting,” in Proc. IEEE Int. Conf. Commun. (ICC), Jun. 2015. [31] Y. L. Che, L. Duan, and R. Zhang, “Spatial throughput maximization of wireless powered communication networks,” IEEE J. Sel. Areas Commun., vol. 33, no. 8, pp. 1534-1548, Aug. 2015.
29
[32] C. Wong, S. Sen, S. Ha, and M. Chiang, “Optimized day-ahead pricing for smart grids with device-specific scheduling flexibility,” IEEE J. Sel. Areas Commun., vol. 30, no. 6, pp. 1075-1085, Jul. 2012. [33] S. Sesia, M. Baker, and I. Toufik, LTE - the UMTS long term evolution: from theory to practice, John Wiley & Sons, 2009. [34] F. Baccelli and B. Blaszczyszyn, Stochastic Geometry and Wireless Networks, Volume I: Theory. NOW: Foundations and Trends in Networking, 2009. [35] J. G. Andrews, F. Baccelli, and R. K. Ganti, “A tractable approach to coverage and rate in cellular networks,” IEEE Trans. Commun., vol. 59, no. 11, pp. 3122-3134, Nov. 2011. [36] Energy efficiency analysis of the reference systems, areas of improvements and target breakdown. Available [online] at http://cordis.europa.eu/docs/projects/cnect/3/247733/080/deliverables/001-EARTHWP2D23v2.pdf. [37] A. Fu, E. Modiano, and J. Tsitsiklis, “Optimal dnergy allocation for delay-constrained data transmission over a time-varying channel,” in Proc. IEEE INFOCOM 2003, Apr. 2003. [38] National Renewable Energy Laboratory. Available [online] at http://www.nrel.gov/midc/unlv/. [39] Hourly Ontario Energy Price. Available [online] at http://www.ieso.ca/Pages/Power-Data/price.aspx.