This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Green Resource Allocation and Energy Management in Heterogeneous Small Cell Networks Powered by Hybrid Energy

Qiaoni Han, Bo Yang, Nan Song, Yuwei Li and Ping Wei Q. Han is with the Department of Automation, Tianjin University, Tianjin 300072, China; B. Yang, N. Song, Y. Li, and P. Wei are with the Key Laboratory of System Control and Information Processing, Ministry of Education of China, and the Department of Automation, Shanghai Jiao Tong University, Shanghai 200240, China.
Abstract

In heterogeneous networks (HetNets), how to improve spectrum efficiency is a crucial issue. Meanwhile increased energy consumption inspires network operators to deploy renewable energy sources as assistance to traditional electricity. Based on above aspects, we allow base stations (BSs) to share their licensed spectrum resource with each other and adjust transmission power to adapt to the renewable energy level. Considering the sharing fairness among BSs, we formulate a multi-person bargaining problem as a stochastic optimization problem. We divide the optimization problem into three parts: data rate control, resource allocation and energy management. An online dynamic control algorithm is proposed to control admission rate and resource allocation to maximize the transmission and sharing profits with the least grid energy consumption. Simulation results investigate the time-varying data control and energy management of BSs and demonstrate the effectiveness of the proposed scheme.

Index Terms:
HetNets, hybrid energy, spectrum sharing, energy management, fairness, Lyapunov optimization

I Introduction

The explosive growth of global mobile data traffic and mounting energy problems have led to higher requirement for both area spectral efficiency and energy efficiency[1]. Future cellular networks are now rapidly evolving to heterogeneous network (HetNet) architectures with small-cell base stations (SBSs) improving quality of service (QoS) and macrocell base stations (MBSs) guaranteeing global coverage area[2]. There are two kinds of spectrum sharing in HetNets: orthogonal sharing and non-orthogonal sharing[3, 4, 5]. Orthogonal sharing allows BSs to operate in the orthogonal subchannels, which lowers the interference with each other. Non-orthogonal spectrum sharing allows BSs to reuse the available spectrum resources with higher interference at the receivers[6].

Although the energy consumption of individual SBS is small, a large number of SBSs lead to higher energy consumption. In fact, the communication networks expend approximately 60 billion kWh per year[7]. High energy consumption of HetNets is an urgent problem to be solved at present, so we want to reduce energy consumption without sacrificing network performance. Energy harvesting technology is an effective method to reduce energy cost[8]. However, due to the fluctuation of energy harvesting, the renewable energy generation may not adapt to the load traffic conditions. To guarantee steady communication, the power grid provides auxiliary power supply.

Spectrum is another vital resource in wireless communications, besides energy[9]. However, existing spectrum allocation schemes often lead to a low spectrum utility, because the BSs of different mobile network operators (MNOs) only consider their own spectrum requirements instead of collective interests[10]. Due to the BSs’ self-interest, we need to design incentives to encourage BSs with idle spectrum resource to share it to improve spectrum utilization.

Therefore, the wireless resource allocation scheme should be adaptive to the finite spectrum resources and stochastic energy sources with the least grid energy consumption to ensure the QoS of mobile users.

I-A Prior Work and Motivation

Previous researches have provided an overview on resource allocation and interference management in a two-tier HetNet in recent years. Palanisamy and Nirmala [11] have presented a comprehensive survey on downlink interference management strategies in HetNets. Zhang et al. [12] studied resource optimization for the interference management and self-organization schemes in a hybrid self-organized small cell network. In [13], the authors considered the sum-rate optimization problem with power control for uplink transmission in a HetNet, then a new practical near-optimal distributed algorithm that eliminates these network overheads was proposed. Lang et al. [14] adopted the heterogeneous genetic algorithm to solve the resource allocation problem for cognitive decode-and-forward relay network. A game-theoretical scheme using energy-efficient resource allocation and interference pricing was proposed in [15] for an interference-limited environment without complete channel state information.

Apart from designing efficient resource allocation policies, there are some energy-saving studies for HetNet. The deployment of renewable energy has been considered in [16, 17, 18, 19]. Gong et al. [16] considered energy-efficient wireless resource management in cellular networks powered by renewable energy. The goal is to minimize the average grid power consumption while satisfying the users’ QoS (blocking probability) requirements. It uses statistical information for traffic intensity and renewable energy to adaptively switch BSs’ on-off state and adjust the allocated resource blocks. Zhang et al. [17] considered coordinated multi-point (CoMP) communication systems and integrated renewable energy sources (RES) into smart grid. A convex optimization problem is formulated to minimize the worst-case energy transaction cost while guaranteeing the QoS of users. Yang et al. [18] jointly considered the resource allocation and energy management of HetNet powered by hybrid energy with the consideration of spatio-temporal diversity of traffic and renewable energy. Qin et al. [19] investigated the resource allocation for orthogonal frequency-division multiple access (OFDMA) wireless networks powered by renewable energy and traditional grid. There is a trade-off between network throughput and grid energy consumption. The proposed main solution is similar to [18].

Due to the finite spectrum resources, several spectrum sharing algorithms have been proposed in literature such as [20, 21, 22, 23, 24, 25]. Specifically, spectrum sharing in the same operator is considered in [20, 21, 22]. Zhang et al. [20] proposed a distributed resource management scheme in a HetNet, which divides the problem into two sub-problems: BS cluster and subchannel allocation. Note that only the BSs in the same cluster are allowed to share the spectrum band. Lee et al. [21] utilized the cognitive radio (CR) technique to improve the performance of wireless powered communication network. Two coexisting models for spectrum sharing were proposed and the authors considered the interference between two models to protect the primary user transmission. The coexisting network including device-to-device links and cellular links was adopted and a new game model called Bayesian overlapping coalition formation game was proposed in [22]. Moreover, multi-MNO coordinate their spectrum resources sharing in [23, 24, 25]. A common spectrum pool to share spectrum was investigated in [23]. It formulates the sharing between the MNOs as a repeated game and determines rules to decide whether to participate in the game at each stage game. An underlay/overlay (hybrid) transmission mode to share the spectrum of licensed users was introduced in [24]. The objective is to maximize network throughput and eliminate interference. Zhang et al. [25] considered the spectrum sharing among multi-MNO in the unlicensed spectrum. A hierarchical game including Kalai-Smorodinsky bargaining game and Stackelberg game is proposed to reduce the interference.

Motivated by the above works, we can use renewable energy to reduce energy consumption and share spectrum to improve spectrum utilization. There are some existing works in this field. But most proposed solutions ignore the fairness between BSs and time-varying of network states. Although some studies, e.g. [26] consider a fair resource allocation, they focused mainly on the network throughput without sharing costs. However, selfish BSs emphasize on improving their own utility without considering other BSs of different MNOs. Thus, we need to design incentive mechanism to promote the BSs to share idle spectrum resource.

In this paper, we consider joint subchannel allocation and power control under cross-tier interference, fairness between SBSs and renewable energy constraints. The main contributions of this paper are as follows:

  • First, we jointly consider subchannel allocation and power control strategy for spectrum sharing in a HetNet, while considering the cross-tier interference constraint.

  • Second, a multi-person bargaining problem is modelled to measure the fairness among SBSs. An online dynamic resource allocation scheme consisting of spectrum pricing, admission flow control and resource allocation is developed by solving a stochastic optimization problem.

  • Third, simulation results are presented to show the effectiveness of the proposed approach. Our approach can substantially achieve the trade-off between transmission profits and network throughput.

The rest of this paper is organized as follows. In Section II, we illustrate the system model and constraint conditions. Section III formulates a stochastic optimization problem based on the Nash bargaining. We present an online dynamic resource allocation and energy management scheme in section IV. Numerical results are given in Section V with discussions. Conclusion and future work are presented in Section VI.

II System and Queue Model

II-A Network Model

We consider a downlink OFDMA two-tier HetNet. As shown in Fig.1, SBSs belonging to different MNOs are overlaid in the coverage of a MBS. There are two types of mobile user equipments (UEs): macro-cell UE (MUE) and small cell UE (SUE) corresponding to the service BSs. All SBSs are assumed to be close access mode that allows only the authorized SUEs to access the corresponding SBS. We assume that the SBSs occupy different spectrum bands. Thus, the co-tier interference between SBSs can be considered to be negligible.

Accordingly, let 𝒩0={0,1,2,,N}\mathcal{N}_{0}=\left\{0,1,2,...,N\right\} denote the BS set, where n=0n=0 represents MBS and 𝒩={1,2,,N}\mathcal{N}=\left\{1,2,...,N\right\} denotes SBSs. As a consequence, the set of MUEs served by MBS is 𝒰0\mathcal{U}_{0} and the set of SUEs associated with SBS nn is 𝒰n={1,2,3,,Un}\mathcal{U}_{n}=\left\{1,2,3,...,U_{n}\right\}. The OFDMA system has a bandwidth FF divided into MM subchannels, i.e., ={1,2,,M}\mathcal{M}=\left\{1,2,...,M\right\} and the bandwidth of subchannel mm is ϖm\varpi_{m}, while the spectrum band occupied by SBS nn is denoted by n\mathcal{M}_{n}\in\mathcal{M}, and n𝒩n=,ij=,ij\cup_{n\in\mathcal{N}}\mathcal{M}_{n}=\mathcal{M},\mathcal{M}_{i}\cap\mathcal{M}_{j}=\varnothing,\forall i\neq j. Moreover, since MBS has the priority to choose subchannels, it is assumed that each subchannel is always occupied by one MUE at a given time, and to guarantee the QoS of MUEs, the SBSs adjust their power to reduce the cross-tier interference.

As the subchannel conditions and energy harvesting process change dynamically, we consider the HetNet operates in a time-slotted manner where the subchannel and energy harvesting conditions keep stable in a time slot. We use 𝒯={0,1,2,,t,,T}\mathcal{T}=\left\{0,1,2,...,t,...,T\right\} to denote the time-slot set and the length of each time slot defaults to 1.

Refer to caption
Fig. 1: An example of the considered network model

At the beginning of time slot tt, each SBS observes the channel state information (CSI) and current queue state information (QSI). We use pmkM(t)p_{mk}^{M}(t) to denote the power allocated to MUE kk served by MBS on subchannel mnm\in\mathcal{M}_{n} and hnmuM(t)h_{nmu}^{M}(t) to denote the interference subchannel gain from MBS to SUE uu in SBS nn on subchannel mnm\in\mathcal{M}_{n}. The received signal-to-interference-plus-noise ratio (SINR) of the uuth SUE from SBS nn on subchannel mnm\in\mathcal{M}_{n} at the tt-th time slot is given by:

γnmu(t)=pnmu(t)hnmu(t)pmkM(t)hnmuM(t)+σ2,mn\mathit{\gamma}_{nmu}(t)=\frac{p_{nmu}(t)h_{nmu}(t)}{p_{mk}^{M}(t)h_{nmu}^{M}(t)+\sigma^{2}},\forall m\in\mathcal{M}_{n} (1)

where pnmu(t)p_{nmu}(t) represents the transmission power of SBS nn to SUE uu on subchannel mnm\in\mathcal{M}_{n}. hnmu(t)h_{nmu}(t) is the channel gain between the SUE uu and SBS nn on subchannel mnm\in\mathcal{M}_{n} including pathloss, shadowing and other factors. σ2\sigma^{2} denotes the noise power level. Compared with the cross-tier interference, the co-tier interference can be considered to be de minimis. Therefore, the SBSs are primarily interfered by the MBS.

The transmission rate of SBS nn to SUE uu on subchannel mnm\in\mathcal{M}_{n} in the non-cooperative case is

Rnmu(t)=ϖmlog2(1+γnmu(t)),mnR_{nmu}(t)=\varpi_{m}\mathrm{\log_{2}}\left(1+\mathit{\gamma}_{nmu}(t)\right),\forall m\in\mathcal{M}_{n} (2)

Then, the total throughput of SUE uu serviced by SBS nn can be computed by:

Rnu(t)=mnxnmu(t)Rnmu(t)R_{nu}(t)=\sum_{m\in\mathcal{M}_{n}}x_{nmu}(t)R_{nmu}(t) (3)

We define a binary variable xnmu(t)x_{nmu}(t) as subchannel allocation index. xnmu(t)=1x_{nmu}(t)=1 indicates the subchannel mnm\in\mathcal{M}_{n} is allocated to SUE uu of SBS nn and otherwise, xnmu(t)=0x_{nmu}(t)=0. In this transmission model, each subchannel can be occupied by at most one SUE at the given time. Limited by the transmission characteristics, we have the following limitations:

Subchannel allocation constraint: Each subchannel can be allocated to at most one SUE in each SBS at time slot tt to avoid intra-cell interference. Each SBS uses different subchannels to avoid inter-tier interference. Thus, we have

n=1Nu=1Unxnmu(t)1,m,t\sum\limits_{n=1}^{N}\sum\limits_{u=1}^{U_{n}}x_{nmu}(t)\leq 1,\forall m,t (4)
xnmu(t){0,1},n,m,u,tx_{nmu}(t)\in\left\{0,1\right\},\forall n,m,u,t (5)

SBS power constraint: Each SBS has a peak transmission power pnmaxp_{n}^{max} to avoid the overload. Each SBS’s transmission power is limited by

mnu𝒰nxnmu(t)pnmu(t)pnmax,n,t\sum\limits_{m\in\mathcal{M}_{n}}\sum\limits_{u\in\mathcal{U}_{n}}x_{nmu}(t)p_{nmu}(t)\leq p_{n}^{max},\forall n,t (6)

Cross-tier interference constraint: The interference from a SBS to a MUE occurs when the SBS transmits in the same subchannel occupied by the MUE. To protect the QoS of MUEs, we denote ImS(t)I_{m}^{S}(t) as the maximum tolerable cross-tier interference temperature in subchannel mm for the assigned MUE bb at time slot tt. Thus, we have

n=1Nu=1Unxnmu(t)pnmu(t)hnmb(t)ImS(t),m,t\sum\limits_{n=1}^{N}\sum\limits_{u=1}^{U_{n}}x_{nmu}(t)p_{nmu}(t)h_{nmb}(t)\leq I_{m}^{S}(t),\forall m,t (7)

where hnmb(t)h_{nmb}(t) is the interference channel gain on subchannel mm from SBS nn to MUE bb served by the MBS at time slot tt. The interference constraint means that SBSs are allowed to transmit signals on the same subchannel with the MBS only if the total interference is kept under a tolerable level.

In this paper, we focus on the joint power control and subchannel allocation problem to improve network performance.

II-B SUE Traffic and Data Queue Model

We assume there is a data buffer for each SUE to store the arriving packets. Let the stochastic process Anu(t)[0,Amax]A_{nu}(t)\in\left[0,A^{max}\right] denote the arrival traffic amount for SUE uu of SBS nn in time slot tt. Anu(t)A_{nu}(t) is independent and identically distributed (i.i.d) over different time slots with the average arrival rate λnu\lambda_{nu}, i.e., 𝔼[Anu(t)]=λnu\mathbb{E}[A_{nu}(t)]=\lambda_{nu}, where 𝔼[]\mathbb{E}[\cdot] represents the expectation. To achieve the normal transmission, the admitted traffic amount Dnu(t)D_{nu}(t) has the following constraint:

0Dnu(t)Anu(t)Amax,n,u,t0\leq D_{nu}(t)\leq A_{nu}(t)\leq A^{max},\forall n,u,t (8)

Thus, the finite data queue backlog Qnu(t)Q_{nu}(t) of SUE uu in small cell nn is formulated as

Qnu(t+1)=[Qnu(t)Rnu(t)]++Dnu(t)Q_{nu}(t+1)=\left[Q_{nu}(t)-R_{nu}(t)\right]^{+}+D_{nu}(t) (9)

where Rnu(t)R_{nu}(t) and Dnu(t)D_{nu}(t) are the output and input rate of data queue QnuQ_{nu} respectively and [x]+=max{x,0}\left[x\right]^{+}=max\left\{x,0\right\}. At the beginning, we assume Qnu(0)=0Q_{nu}(0)=0. Due to the time-varying subchannels and packet arrival process, the arrival and departure processes are both stochastic. Thus, the QnuQ_{nu} is varying over time slot. Let R(t)R(t) and D(t)D(t) be the output rate and input rate at time slot t for queue Q(t)Q(t), respectively. For discrete time process Q(t)Q(t) evolves as following,

Q(t+1)=[Q(t)R(t)]++D(t)Q(t+1)=\left[Q(t)-R(t)\right]^{+}+D(t) (10)

we need to guarantee the queue stability. According to [27], a queue is defined as strongly stable if

limT1Tt=0T1𝔼[Q(t)]<\underset{T\rightarrow\infty}{\lim}\frac{1}{T}\sum_{t=0}^{T-1}\mathbb{E}\left[Q(t)\right]<\infty (11)

A network is called strongly stable if all the individual queues are strongly stable. Therefore, it is a crucial prerequisite for the network to guarantee the stability of all queues.

II-C Energy Harvesting and Energy Queue Model

According to the EARTH project,we adopt a linear approximation power consumption as[28]:

Pn(t)=Pnc+Δnu=1Unm=1Mxnmu(t)pnmu(t)P_{n}(t)=P_{n}^{c}+\Delta_{n}\sum\limits_{u=1}^{U_{n}}\sum_{m=1}^{M}x_{nmu}(t)p_{nmu}(t) (12)

where PncP_{n}^{c} is the static power consumption including the cooling system, baseband processor and so on. Δn\Delta_{n} is the slope of the transmission power dependent power consumptions. According to the maximum power constraint (6), the power consumption satisfies the following constraint:

Pn(t)Pnc+Δnpnmax=ΔPnmaxP_{n}(t)\leq P_{n}^{c}+\Delta_{n}p_{n}^{max}\overset{\Delta}{=}P_{n}^{max} (13)

To avoid the transmission interruption, the MBS uses a conventional grid source and the SBSs are powered by not only renewable energy generation but also the power grid. Each SBS harvests energy from ambient energy sources such as solar and stores the energy in its battery. If the storage energy is not enough for transmission, the SBS will be supplied by the grid. We model the energy harvest process as a random process. Let En(t)E_{n}(t) denote the harvested energy in time slot tt in small cell nn. En(t)E_{n}(t) is assumed to be i.i.d. with the maximum value EnmaxE_{n}^{max}. Notice that SBS may have no a priori knowledge of energy harvest process, which is appropriate for the situation.

At time slot tt, the battery level Sn(t)S_{n}(t) is determined by previous battery level, harvested energy and SBS power consumption as follows:

Sn(t+1)=[Sn(t)Fn(t)]++Jn(t)S_{n}(t+1)=\left[S_{n}(t)-F_{n}(t)\right]^{+}+J_{n}(t) (14)

where Fn(t)F_{n}(t) and Jn(t)J_{n}(t) are the discharge energy and charge energy of the battery device in SBS nn at time slot tt respectively. The storage energy cannot exceed the capacity of the battery, i.e., Sn(t)SnmaxS_{n}(t)\leq S_{n}^{max}.

The power consumption of SBS nn consists of two parts: discharge energy from battery devices Fn(t)F_{n}(t) and power grid Gn(t)G_{n}(t). Then, we have

Pn(t)=Fn(t)+Gn(t)P_{n}(t)=F_{n}(t)+G_{n}(t) (15)

The discharge energy should be less than the battery energy level and the charge energy should consider the battery capacity. Thus, the discharge energy and charge energy should satisfy the following constraints:

0Fn(t)Sn(t)0\leq{F_{n}}\left(t\right)\leq{S_{n}}\left(t\right) (16)
0Jn(t)min{SnmaxSn(t),En(t)}0\leq J_{n}(t)\leq\min\left\{S_{n}^{max}-S_{n}(t),E_{n}(t)\right\} (17)

Let F¯n=limT1Tt=0T1𝔼{Fn(t)}\overline{F}_{n}=\underset{T\rightarrow\infty}{\lim}\frac{1}{T}\sum_{t=0}^{T-1}\mathbb{E}\left\{F_{n}(t)\right\} and J¯n=limT1Tt=0T1𝔼{Jn(t)}\overline{J}_{n}=\underset{T\rightarrow\infty}{\lim}\frac{1}{T}\sum_{t=0}^{T-1}\mathbb{E}\left\{J_{n}(t)\right\} denote the time-average discharge energy and the time-average charge energy respectively. The stability of energy queue Sn(t)S_{n}(t) should be guaranteed.

II-D Spectrum Sharing

At each time slot, the common spectrum pool formed by MNOs from reusing the spectrum of MBS is partitioned between the SBSs. The total spectrum band allocated to SBS nn at time slot tt is Bn(t)=u𝒰nmnϖmxnmu(t)B_{n}(t)=\sum\limits_{u\in\mathcal{U}_{n}}\sum\limits_{m\in\mathcal{M}_{n}}\varpi_{m}x_{nmu}(t). We assume that each SBS uses the same amount of subchannels without any payment in the initial state. However, the SBS pays for an extra payment if it uses more subchannels than it has put in the common pool. Similarly, a SBS can rent its free spectrum resource to other SBSs for benefits.

We define a utility of SUE uu of SBS nn in terms of admitted data rate to represent the satisfaction that SBS nn sends data to SUE uu, then the time-average expected profit of SBS nn is given as

C¯n=u=1UnD¯nu(t)+O¯nϕG¯n\overline{C}_{n}=\sum_{u=1}^{U_{n}}\overline{D}_{nu}(t)+\overline{O}_{n}-\phi\overline{G}_{n} (18)

where

G¯n=limT1Tt=0T1𝔼{Gn(t)}\overline{G}_{n}=\underset{T\rightarrow\infty}{\lim}\frac{1}{T}\sum_{t=0}^{T-1}\mathbb{E}\left\{G_{n}(t)\right\}\\ (19)
O¯n=limT1Tt=0T1𝔼{On(t)}\overline{O}_{n}=\underset{T\rightarrow\infty}{\lim}\frac{1}{T}\sum_{t=0}^{T-1}\mathbb{E}\left\{O_{n}(t)\right\}\\ (20)
On(t)=αn(t)(Bn0Bn(t))+βn(t)(Bn(t)Bn0)+O_{n}(t)=\alpha_{n}(t)\left(B_{n}^{0}-B_{n}(t)\right)^{+}-\beta_{n}(t)\left(B_{n}(t)-B_{n}^{0}\right)^{+} (21)

ϕ\phi is a constant to balance gains from spectrum capacity and power consumption, which is a constant. Bn0B_{n}^{0} denotes the initial spectrum band of SBS nn. αn(t)\alpha_{n}(t) is the per-unit price of spectrum transferred from SBS nn to other SBSs and βn(t)\beta_{n}(t) is the per-unit price of extra spectrum rented from other SBSs. The spectrum prices all have the maximum qnmaxq_{n}^{max}. Notice that during any time slot, the SBS nn can only choose one spectrum sharing scheme: lease or rent spectrum band.

III Problem Formulation

In this paper, our objective is to maximize the total time-average profit gains of SBSs (or MNOs) by sharing their spectrum resources with each other and minimizing the grid power consumption at the same time. Furthermore, we need to allocate the spectrum resources fairly so that the SBSs have an incentive to share their spectrum. Meanwhile the proposed resource allocation and energy management scheme should satisfy interference constraints, energy constraints and network stability. Then, based on the Nash bargaining problem studied in [29], the optimization problem can be formulated as

OP1:max𝐃,𝐗,𝐏,𝐆,𝜶,𝜷n=1Nlog(C¯nCnmin)s.t.C1:C¯nCnmin,nC2:0Dnu(t)Anu(t)Amax,n,u,tC3:m=1Mu=1Unxnmu(t)pnmu(t)pnmax,nC4:n=1Nu=1Unxnmu(t)pnmu(t)hnmu(t)ImS(t),mC5:n=1Nu=1Unxnmu(t)1,n,mC6:xnmu(t){0,1},n,m,uC7:0Fn(t)Sn(t)C8:0Jn(t)min{SnmaxSn(t),En(t)}C9:Sn(t+1)=[Sn(t)Fn(t)]++Jn(t)C10:Queues Qn,u(t)are strongly stable,n,u\begin{split}\mathrm{OP}_{1}:&\quad\max\limits_{\mathbf{D,X,P,G},\bm{\alpha,\beta}}\sum\limits_{n=1}^{N}\log(\overline{C}_{n}-{C_{n}^{min}})\\ \text{s.t.}\quad&C1:\overline{C}_{n}\geq C_{n}^{min},\forall n\\ &C2:0\leq D_{nu}(t)\leq A_{nu}(t)\leq A^{max},\forall n,u,t\\ &C3:\sum\limits_{m=1}^{M}\sum\limits_{u=1}^{U_{n}}x_{nmu}(t)p_{nmu}(t)\leq p_{n}^{max},\forall n\\ &C4:\sum\limits_{n=1}^{N}\sum\limits_{u=1}^{U_{n}}x_{nmu}(t)p_{nmu}(t)h_{nmu}(t)\leq I_{m}^{S}(t),\forall m\\ &C5:\sum\limits_{n=1}^{N}\sum\limits_{u=1}^{U_{n}}x_{nmu}(t)\leq 1,\forall n,m\\ &C6:x_{nmu}(t)\in\left\{0,1\right\},\forall n,m,u\\ &C7:0\leq{F_{n}}\left(t\right)\leq{S_{n}}\left(t\right)\\ &C8:0\leq J_{n}(t)\leq\min\left\{S_{n}^{max}-S_{n}(t),E_{n}(t)\right\}\\ &C9:S_{n}(t+1)=\left[S_{n}(t)-F_{n}(t)\right]^{+}+J_{n}(t)\\ &C10:\text{Queues }Q_{n,u}\left(t\right)\text{are strongly stable},\forall n,u\end{split} (22)

where variables 𝐃={Dnu(t)}\mathbf{D}=\left\{D_{nu}(t)\right\}, 𝐗={xnmu(t)}\mathbf{X}=\left\{x_{nmu}(t)\right\}, 𝐏={pnmu(t)}\mathbf{P}=\left\{p_{nmu}(t)\right\}, 𝐆={Gn(t)}\mathbf{G}=\left\{G_{n}(t)\right\} are admitted data rate, subchannel assignment, power control and grid energy consumption variables of the network respectively. C1C1 ensures that SBSs will benefit from sharing spectrum resources and CnminC_{n}^{min} is the minimal profit of SBS nn in the non-cooperative way, which is set to be a time-invariant non-negative constant. C2C2-C9C9 are the admitted data amount, subchannel, power and energy constraints. C10C10 is imposed to ensure the network stability.

As for the objective function in optimization problem OP1, since the loglog is a continuous strictly increasing function, maxn=1Nlog(C¯nCnmin)\max\sum_{n=1}^{N}\log(\overline{C}_{n}-{C_{n}^{min}}) is equivalent to maxn=1N(C¯nCnmin)\max\prod_{n=1}^{N}(\overline{C}_{n}-{C_{n}^{min}}), which is the Nash product. Based on the Nash Bargaining Solution (NBS), by maximizing the Nash product, or equivalently the loglog sum in optimization problem OP1, a unique and fair Pareto optimality solution can be provided [29].

According to the energy queue model, we can find the current energy management policy has a coupling relationship with future energy state. The coupling nature, which comes with stochastic renewable energy sources, makes OP1\mathrm{OP}_{1} more complicated. By taking iterated expectation and using telescoping sums in (14) over t=0,1,,T1t=0,1,...,T-1, we can get

𝔼{S(T1)}S(0)=t=0T1𝔼{Fn(t)+Jn(t)}\mathbb{E}\left\{S(T-1)\right\}-S(0)=\sum_{t=0}^{T-1}\mathbb{E}\left\{-F_{n}(t)+J_{n}(t)\right\} (23)

As the energy level of the battery devices Sn(t)S_{n}(t) is bounded by SnmaxS_{n}^{max} , we have F¯n=J¯n\overline{F}_{n}=\overline{J}_{n} after dividing both sides of (23) by TT and taking limitation of TT\rightarrow\infty.

Moreover, it is noted that, the optimization problem OP1\text{OP}_{1} aims to maximize an objective function of a time-averaged variable C¯nCnmin\overline{C}_{n}-{C_{n}^{min}}, which is difficult to be tackled. To solve OP1\text{OP}_{1}, we introduce an auxiliary variable μn\mu_{n} to reformulate problem OP1\mathrm{OP}_{1} as

OP2:max𝐃,𝐗,𝐏,𝐆,𝜶,𝜷n=1Nlog(μn¯)s.t.C1:μ¯nC¯nCnminC2:C¯nCnmin,nC3:C2C6C4:F¯n=J¯nC5:Queues Qn,u(t)are strongly stable,n,u\begin{split}\mathrm{OP}_{2}:&\quad\max\limits_{\mathbf{D,X,P,G},\bm{\alpha,\beta}}\sum\limits_{n=1}^{N}\log(\overline{\mu_{n}})\\ \text{s.t.}\quad&C1^{\prime}:\overline{\mu}_{n}\leq\overline{C}_{n}-{C_{n}^{min}}\\ &C2^{\prime}:\overline{C}_{n}\geq C_{n}^{min},\forall n\\ &C3^{\prime}:C2-C6\\ &C4^{\prime}:\overline{F}_{n}=\overline{J}_{n}\\ &C5^{\prime}:\text{Queues }Q_{n,u}\left(t\right)\text{are strongly stable},\forall n,u\end{split} (24)

where μn\mu_{n} can be understood as the lower bound of C¯nCnmin\overline{C}_{n}-{C_{n}^{min}} as shown in C1C1^{\prime}. Based on the Jensen inequality, log(μn)¯\overline{\log(\mu_{n})} is the lower bound of log(μ¯n)\log(\overline{\mu}_{n}) with the non-decreasing concave logarithmic utility function. Thus, it is feasible to change the objective function log(μ¯n)\log(\overline{\mu}_{n}) to log(μn)¯\overline{\log(\mu_{n})}, i.e., the optimization problem OP2\text{OP}_{2} is an optimization problem with a time-averaged objection function, which is beneficial for the following problem solving.

To satisfy the time-average constraints C1C1^{\prime} and C2C2^{\prime}, we introduce the concept of virtual queue technology[27], and use virtual queues of Y(t)Y(t) and Z(t)Z(t). Specifically, the virtual queue Y(t)Y(t) associated with constraint C1C1^{\prime} updates as follows:

Yn(t+1)=[Yn(t)ynout(t)]++ynin(t)Y_{n}(t+1)=\left[Y_{n}(t)-y_{n}^{out}(t)\right]^{+}+y_{n}^{in}(t) (25)

where

ynin(t)=μn(t)+Cnmin+ϕGn(t)y_{n}^{in}(t)=\mu_{n}(t)+C_{n}^{min}+\phi G_{n}(t) (26)
ynout(t)=u=1UnDnu(t)+On(t)y_{n}^{out}(t)=\sum_{u=1}^{U_{n}}{D}_{nu}(t)+O_{n}(t) (27)

The virtual queue Z(t)Z(t) associated with constraint C2C2^{\prime} updates as follows:

Zn(t+1)=[Zn(t)znout(t)]++znin(t)Z_{n}(t+1)=\left[Z_{n}(t)-z_{n}^{out}(t)\right]^{+}+z_{n}^{in}(t) (28)

where

znin(t)=Cnmin+ϕGn(t)z_{n}^{in}(t)=C_{n}^{min}+\phi G_{n}(t) (29)
znout(t)=u=1UnDnu(t)+On(t)z_{n}^{out}(t)=\sum_{u=1}^{U_{n}}{D}_{nu}(t)+O_{n}(t) (30)

If the virtual queues of Yn(t)Y_{n}(t) and Zn(t)Z_{n}(t) are strongly stable, then the constraints C1C1^{\prime} and C2C2^{\prime} are satisfied.

III-A Lyapunov Optimization

According to the queues, 𝐐={Qnu(t)}\mathbf{Q}=\left\{Q_{nu}(t)\right\}, 𝐒={Sn(t)}\mathbf{S}=\left\{S_{n}(t)\right\}, 𝐘={Yn(t)}\mathbf{Y}=\left\{Y_{n}(t)\right\}, 𝐙={Zn(t)}\mathbf{Z}=\left\{Z_{n}(t)\right\}, and Θ(t)=[𝐐,𝐒,𝐘,𝐙]\Theta(t)=[\mathbf{Q},\mathbf{S},\mathbf{Y},\mathbf{Z}]. We define the Lyapunov function as

L(Θ(t))=12[n=1Nu=1UnQnu(t)2+n=1N(Sn(t)ρn)2+n=1NYn2(t)+n=1NZn2(t)]\begin{split}L(\Theta(t))=\frac{1}{2}\left[\sum\limits_{n=1}^{N}\sum\limits_{u=1}^{U_{n}}Q_{nu}(t)^{2}+\sum\limits_{n=1}^{N}\left(S_{n}(t)-\rho_{n}\right)^{2}+\sum\limits_{n=1}^{N}Y_{n}^{2}(t)+\sum\limits_{n=1}^{N}Z_{n}^{2}(t)\right]\end{split} (31)

where ρn\rho_{n} is a perturbation factor, which ensures there is enough energy in the energy queue. Without loss of generality, we assume that all queues are empty when t=0t=0 such that L(Θ(t))=0L(\Theta(t))=0.

The Lyapunov function L(Θ(t))L(\Theta(t)) is a scalar measure of network congestion. Intuitively, if L(Θ(t))L(\Theta(t)) is small then all queues are small; and if L(Θ(t))L(\Theta(t)) is large then at least one queue is large. Thus, by minimizing the drift in the Lyapunov function (i.e., by minimizing a difference in the Lyapunov function from one time slot to the next), queues Qnu(t)Q_{nu}(t), Sn(t)S_{n}(t), Yn(t)Y_{n}(t), Zn(t)Z_{n}(t) can be stabilized. By using expression (31), the drift in the Lyapunov function (i.e., the expected change in the Lyapunov function from one time slot to the next) can be written as

(Θ(t))=𝔼{L(Θ(t+1))L(Θ(t))|Θ(t)}\bigtriangleup(\Theta(t))=\mathbb{E}\left\{L(\Theta(t+1))-L(\Theta(t))|\Theta(t)\right\} (32)

We now use the drift-plus-penalty minimization method to solve optimization problem OP2. In this method, a control policy that solves OP2 is obtained by minimizing the upper bound on the following drift-plus-penalty expression. We can obtain the following drift-plus-penalty term:

ΔV(t)=(Θ(t))Vn=1N𝔼{f(μn(t))|Θ(t)}penaltyterm\Delta_{V}(t)=\bigtriangleup(\Theta(t))\underset{\mathrm{penalty\ term}}{\underbrace{-V\sum_{n=1}^{N}\mathbb{E}\left\{f(\mu_{n}(t))|\Theta(t)\right\}}} (33)

where VV is a non-negative tunable parameter which balances the maximization of network utility and the minimization of queue length to a state level. According to [30], when VV is sufficiently large, the optimization algorithm approaches the optimal capacity.

Lemma 1. The drift-plus-penalty term is upper bounded as

ΔV(t)H+n=1Nu=1UnQnu(t)𝔼{Dnu(t)Rnu(t)|Θ(t)}\displaystyle\Delta_{V}(t)\leq H+\sum\limits_{n=1}^{N}\sum\limits_{u=1}^{U_{n}}Q_{nu}(t)\mathbb{E}\left\{D_{nu}(t)-R_{nu}(t)|\Theta(t)\right\}
+n=1NWn(t)𝔼{Cnmin+ϕGn(t)u=1UnDnu(t)On(t)|Θ(t)}\displaystyle+\sum\limits_{n=1}^{N}W_{n}(t)\mathbb{E}\Bigg{\{}C_{n}^{min}+\phi G_{n}(t)-\sum_{u=1}^{U_{n}}{D}_{nu}(t)-O_{n}(t)|\Theta(t)\Bigg{\}}
+n=1N{(Sn(t)ρn)𝔼{Jn(t)Pn(t)+Gn(t)|Θ(t)}\displaystyle+\sum\limits_{n=1}^{N}\Bigg{\{}\left(S_{n}(t)-\rho_{n}\right)\mathbb{E}\left\{J_{n}(t)-P_{n}(t)+G_{n}(t)|\Theta(t)\right\}
+Yn(t)𝔼{μn(t)|Θ(t)}}Vn=1N𝔼{f(μn(t))|Θ(t)}\displaystyle+Y_{n}(t)\mathbb{E}\left\{\mu_{n}(t)|\Theta(t)\right\}\Bigg{\}}-V\sum_{n=1}^{N}\mathbb{E}\left\{f(\mu_{n}(t))|\Theta(t)\right\} (34)

where Wn(t)=Yn(t)+Zn(t)W_{n}(t)=Y_{n}(t)+Z_{n}(t) denotes the value of virtual queues. HH is a positive constant that satisfies the following inequality constraint for all time slots:

H\displaystyle H\geq 12[n=1Nu=1UnAmax2+n=1NJnmax2+n=1NSnmax2+n=1Nu=1Un𝔼{Rnu(t)2|Θ(t)}\displaystyle\frac{1}{2}\Bigg{[}\sum\limits_{n=1}^{N}\sum\limits_{u=1}^{U_{n}}A_{max}^{2}+\sum\limits_{n=1}^{N}J_{n}^{max2}+\sum\limits_{n=1}^{N}S_{n}^{max2}+\sum\limits_{n=1}^{N}\sum\limits_{u=1}^{U_{n}}\mathbb{E}\left\{R_{nu}(t)^{2}|\Theta(t)\right\}
+n=1N𝔼{ynin(t)2+ynout(t)2|Θ(t)}+n=1N𝔼{znin(t)2+znout(t)2|Θ(t)}]\displaystyle+\sum\limits_{n=1}^{N}\mathbb{E}\left\{y_{n}^{in}(t)^{2}+y_{n}^{out}(t)^{2}|\Theta(t)\right\}+\sum\limits_{n=1}^{N}\mathbb{E}\left\{z_{n}^{in}(t)^{2}+z_{n}^{out}(t)^{2}|\Theta(t)\right\}\Bigg{]} (35)

proof: See Appendix A.

By Lemma 1, we have transformed the optimization problem OP2 into minimizing the right-side term of (73) at each time slot tt. According to [26], the control policy should be adjusted to minimize the upper bound. Thus, we will decompose the optimization problem and present an online dynamic control algorithm for the green resource allocation and energy management.

IV Online Control Algorithm for SBSs

In this section, we focus on the online algorithm design based on the previous subsection. In the multi-person bargaining problem, SBSs will compete with the total spectrum band. We propose an online dynamic control algorithm by drift-plus-penalty method in Algorithm 1. At each time slot, the Algorithm 1 is partitioned to four steps:

1) Admitted data rate control: each SBS decides the admission rate of each SUE according to the data queue length and the virtual queue length;

2) Auxiliary variable decision: each SBS decides the auxiliary variable μn(t)\mu_{n}(t);

3) Spectrum pricing and resource allocation: each SBS decides either lease or rent the spectrum band; and each SBS designs subchannel assignment and power allocation scheme;

4) Battery energy management: each SBS decides the amount of discharge energy and charge energy to reduce the grid energy consumption according to the perturbation variable.

Algorithm 1 Online Dynamic Control Algorithm for Green Resource Allocation and Energy Management
  Initial Stage:
  At the beginning of all time slot, initialize all the queues Qnu(0),Sn(0),Yn(0),Zn(0)=0Q_{nu}(0),S_{n}(0),Y_{n}(0),Z_{n}(0)=0
  Step 1–Admission rate control:
  For each SBS nn, calculate the flow rate by solving the problem in (36).
  Step 2–Auxiliary variable decision:
  For each SBS nn, calculate the variable μn(t)\mu_{n}(t) by solving the problem in (38).
  Step 3–Adaptive resource allocation:
  For each SBS nn, compute spectrum sharing variables αn(t)\alpha_{n}(t), βn(t)\beta_{n}(t) by solving the problem in (40) and resource allocation variables pnmu(t)p_{nmu}(t), xnmu(t)x_{nmu}(t) by solving the problem in (43).
  Step 4–Battery energy management:
  For each SBS nn, calculate the discharge energy, charge energy and grid energy by solving the problem in (IV-D).
  End Stage:
  Update the queue state and go to Step 1.

IV-A Admission Rate Control

By observing the second and third terms on the right-hand side of (73), we can get the rate control problem regardless of the other variables:

maxWn(t)Dnu(t)Qnu(t)Dnu(t)\displaystyle\max W_{n}(t)D_{nu}(t)-Q_{nu}(t)D_{nu}(t) (36)
s.t.0Dnu(t)Anu(t),n,u,t\displaystyle\text{s.t.}\quad 0\leq D_{nu}(t)\leq A_{nu}(t),\forall n,u,t

We can observe the monotonicity (36) and then the solution of the above problem can be solved easily,

Dnu(t)={Anu(t),Wn(t)Qnu(t)0,Wn(t)<Qnu(t)\displaystyle D_{nu}(t)=\left\{\begin{matrix}A_{nu}(t),&W_{n}(t)\geq Q_{nu}(t)\\ 0,&W_{n}(t)<Q_{nu}(t)\end{matrix}\right. (37)

IV-B Auxiliary Variable Decision

Each SBS decides the auxiliary variable μn(t)\mu_{n}(t) by

maxVlog(μn(t))Yn(t)μn(t)\displaystyle\max V\log(\mu_{n}(t))-Y_{n}(t)\mu_{n}(t) (38)
s.t.0μn(t)μnmax,n,t\displaystyle\mathrm{s.t.}\quad 0\leq\mu_{n}(t)\leq\mu_{n}^{max},\forall n,t

log(x)\log(x) is the non-decreasing concave utility function, and then we get the optimal solution. μnmax\mu_{n}^{max} is the maximal profit by sharing spectrum resources.

μn(t)={μnmax,Yn(t)Vlog(μnmax)μn,Vlog(μnmax)<Yn(t)<Vlog(t)0,Yn(t)Vlog(t)\displaystyle\mu_{n}(t)=\left\{\begin{matrix}\mu_{n}^{max},&Y_{n}(t)\leq V\log^{\prime}(\mu_{n}^{max})\\ \mu_{n}^{*},&V\log^{\prime}(\mu_{n}^{max})<Y_{n}(t)<V\log^{\prime}(t)\\ 0,&Y_{n}(t)\geq V\log^{\prime}(t)\end{matrix}\right. (39)

where μn\mu_{n}^{*} should satisfy Vlog(μn)=Yn(t)V\log^{\prime}(\mu_{n}^{*})=Y_{n}(t).

IV-C Dynamic Resource Allocation

We can formulate the pricing problem as

maxn=1NWn(t)(\displaystyle\max\sum_{n=1}^{N}W_{n}(t)\Big{(} αn(t)(Bn0Bn(t))+βn(t)(Bn(t)Bn0)+)\displaystyle\alpha_{n}(t)\left(B_{n}^{0}-B_{n}(t)\right)^{+}-\beta_{n}(t)\left(B_{n}(t)-B_{n}^{0}\right)^{+}\Big{)}
s.t.0αn(t),βn(t)qnmax\displaystyle\text{s.t.}\quad 0\leq\alpha_{n}(t),\beta_{n}(t)\leq q_{n}^{max} (40)

We can find that the subchannel price αn(t)\alpha_{n}(t) and βn(t)\beta_{n}(t) are related to the state of virtual queue Wn(t)W_{n}(t). The longer the queue length of WnW_{n}, the more utility gain the SBS nn will obtain by renting its free spectrum resources and vice versa. Moreover, it is easily seen that the pricing problem is tightly coupled with the power and subchannel allocation problem. In the following, to facilitate problem solving, we firstly analyze the resource allocation in two-SBS case, and then extend the results to the multi-SBS case.

1) Two-SBS case (N=2N=2): According to the spectrum pricing problem, the SBS with the longer queue length of WnW_{n}, which is supposed the SBS 1, will rent spectrum resources to another one (SBS 2). The price is chosen as

[α,β]={α1(t)=q1max,β1(t)=0α2(t)=0,β2(t)=q2max\left[\alpha,\beta\right]=\left\{\begin{matrix}\alpha_{1}(t)=q_{1}^{max},\beta_{1}(t)=0\\ \alpha_{2}(t)=0,\beta_{2}(t)=q_{2}^{max}\end{matrix}\right. (41)

By observing the remaining term on the right-hand side of (73), we can formulate the resource allocation problem with the given spectrum sharing scheme as follows:

max𝐩,𝐱n=1Nu=1UnQnu(t)Rnu(t)+n=1NWn(t)On(t)+n=1N(Sn(t)ρn)Pn(t)s.t.C3:m=1Mu=1Unxnmu(t)pnmu(t)pnmax,nC4:n=1Nu=1Unxnmu(t)pnmu(t)hnmu(t)ImS(t),mC5:n=1Nu=1Unxnmu(t)1,mC6:xnmu(t)[0,1],n,m,uC11:pnmu(t)0,n,m,u\begin{split}\underset{\mathbf{p,x}}{\max}\sum\limits_{n=1}^{N}\sum\limits_{u=1}^{U_{n}}&Q_{nu}(t)R_{nu}(t)+\sum\limits_{n=1}^{N}W_{n}(t)O_{n}(t)+\sum\limits_{n=1}^{N}\left(S_{n}(t)-\rho_{n}\right)P_{n}(t)\\ \text{s.t.}\ &C3:\sum\limits_{m=1}^{M}\sum\limits_{u=1}^{U_{n}}x_{nmu}(t)p_{nmu}(t)\leq p_{n}^{max},\forall n\\ &C4:\sum\limits_{n=1}^{N}\sum\limits_{u=1}^{U_{n}}x_{nmu}(t)p_{nmu}(t)h_{nmu}(t)\leq I_{m}^{S}(t),\forall m\\ &C5:\sum\limits_{n=1}^{N}\sum\limits_{u=1}^{U_{n}}x_{nmu}(t)\leq 1,\forall m\\ &C6:x_{nmu}(t)\in\left[0,1\right],\forall n,m,u\\ &C11:p_{nmu}(t)\geq 0,\forall n,m,u\\ \end{split} (42)

We relax the binary variable xnmu(t)x_{nmu}(t) to a continuous variable x^nmu(t)[0,1]\hat{x}_{nmu}(t)\in\left[0,1\right]. For notational brevity, denote the power allocated to SUE uu on subchannel mm as snmu=x^nmupnmus_{nmu}=\hat{x}_{nmu}p_{nmu}. The optimization problem can be rewritten as:

max𝐩,𝐱n=1Nm=1Mu=1Un(Qnu(t)ϖmx^nmu(t)×log2(1+snmu(t)hnmu(t)x^nmu(t)(I0mu(t)+σ2))+ηn(t)snmu(t)+θnm(t)x^nmu(t))s.t.C3¯:m=1Mu=1Unsnmu(t)pnmax,nC4¯:n=1Nu=1Unsnmu(t)hnmu(t)ImS(t),mC5¯:n=1Nu=1Unx^nmu(t)1,mC6¯:x^nmu(t)[0,1]n,m,uC11¯:snmu(t)0,n,m,u\begin{split}\underset{\mathbf{p,x}}{\max}&\sum\limits_{n=1}^{N}\sum\limits_{m=1}^{M}\sum\limits_{u=1}^{U_{n}}\Bigg{(}Q_{nu}(t)\varpi_{m}\hat{x}_{nmu}(t)\times\log_{2}\left({1+\frac{{{s_{nmu}}\left(t\right){h_{nmu}}\left(t\right)}}{{\hat{x}_{nmu}(t)\left({{I_{0mu}}\left(t\right)+{\sigma^{2}}}\right)}}}\right)+\eta_{n}(t)s_{nmu}(t)+\theta_{nm}(t)\hat{x}_{nmu}(t)\Bigg{)}\\ \text{s.t.}\ &\bar{C3}:\sum\limits_{m=1}^{M}\sum\limits_{u=1}^{U_{n}}s_{nmu}(t)\leq p_{n}^{max},\forall n\\ &\bar{C4}:\sum\limits_{n=1}^{N}\sum\limits_{u=1}^{U_{n}}s_{nmu}(t)h_{nmu}(t)\leq I_{m}^{S}(t),\forall m\\ &\bar{C5}:\sum\limits_{n=1}^{N}\sum\limits_{u=1}^{U_{n}}\hat{x}_{nmu}(t)\leq 1,\forall m\\ &\bar{C6}:\hat{x}_{nmu}(t)\in\left[0,1\right]\forall n,m,u\\ &\bar{C11}:s_{nmu}(t)\geq 0,\forall n,m,u\\ \end{split} (43)

where θnm(t)\theta_{nm}(t) is related to αn(t)\alpha_{n}(t), βn(t)\beta_{n}(t), ϖm\varpi_{m} and Wn(t)W_{n}(t), θnm(t)=qnmaxϖmWn(t){\theta_{nm}}\left(t\right)=-q_{n}^{\max}{\varpi_{m}}{W_{n}}\left(t\right) and ηn(t)=Sn(t)ρnϕWn(t)\eta_{n}(t)=S_{n}(t)-\rho_{n}-\phi W_{n}(t). I0mu=pmkM(t)hnmuM(t)I_{0mu}=p_{mk}^{M}(t)h_{nmu}^{M}(t). The subchannel and power allocation strategy in (43) can be solved by using the Lagrangian dual decomposition method. By ignoring the time variables, the partial Lagrangian function is given by

L(x\displaystyle L(x ,s,λ1,λ2,λ3)=n=1Nm=1MLnm(x,s,λ1,λ2,λ3)+n=1Nλ1,npnmax+m=1Mλ2,mImS+m=1Mλ3,m\displaystyle,s,\lambda_{1},\lambda_{2},\lambda_{3})=\sum\limits_{n=1}^{N}\sum\limits_{m=1}^{M}L_{nm}(x,s,\lambda_{1},\lambda_{2},\lambda_{3}){\rm{+}}\sum\limits_{n=1}^{N}{{\lambda_{1,n}}p_{n}^{\max}+}\sum\limits_{m=1}^{M}{{\lambda_{2,m}}I_{m}^{S}+}\sum\limits_{m=1}^{M}{{\lambda_{3,m}}} (44)

with

Lnm(x,s,λ1,λ2,λ3)=\displaystyle L_{nm}(x,s,\lambda_{1},\lambda_{2},\lambda_{3})= u=1Un[Qnu(t)ϖmx^nmu(t)×log2(1+snmu(t)hnmu(t)x^nmu(t)(I0mu(t)+σ2))\displaystyle\sum\limits_{u=1}^{U_{n}}\Bigg{[}Q_{nu}(t)\varpi_{m}\hat{x}_{nmu}(t)\times\log_{2}\left({1+\frac{{{s_{nmu}}\left(t\right){h_{nmu}}\left(t\right)}}{{\hat{x}_{nmu}(t)\left({{I_{0mu}}\left(t\right)+{\sigma^{2}}}\right)}}}\right)
+ηnsnmu(t)+θnm(t)x^nmu(t)]λ1,nu=1Unsnmu(t)\displaystyle+\eta_{n}s_{nmu}(t)+\theta_{nm}(t)\hat{x}_{nmu}(t)\Bigg{]}-\lambda_{1,n}\sum\limits_{u=1}^{U_{n}}s_{nmu}(t)
λ2,mu=1Unsnmu(t)hnmu(t)λ3,mu=1Unx^nmu(t)\displaystyle-\lambda_{2,m}\sum\limits_{u=1}^{U_{n}}s_{nmu}(t)h_{nmu}(t)-\lambda_{3,m}\sum\limits_{u=1}^{U_{n}}\hat{x}_{nmu}(t)

where λ1\lambda_{1}, λ2\lambda_{2} and λ3\lambda_{3} are the Lagrange multipliers for constraints C3¯\bar{C3}, C4¯\bar{C4}, C5¯\bar{C5} in (43) respectively. The boundary constraints C5¯\bar{C5} and C11¯\bar{C11} will be absorbed in the Karush-Kuhn-Tucker (KKT) conditions. Thus, the Lagrangian dual function is defined as:

g(λ)=maxx,sL(x,s,λ1,λ2,λ3)g(\lambda)=\max\limits_{x,s}L(x,s,\lambda_{1},\lambda_{2},\lambda_{3}) (45)

The dual problem can be expressed as:

ming(λ1,λ2,λ3)s.t.λ1,λ2,λ30\begin{split}\min\limits&\quad g(\lambda_{1},\lambda_{2},\lambda_{3})\\ \text{s.t.}&\quad\lambda_{1},\lambda_{2},\lambda_{3}\geq 0\end{split} (46)

According to the KKT conditions, the optimal solutions of the problem should satisfy the following conditions:

Lnm(t)snmu(t)\displaystyle\frac{\partial L_{nm}(t)}{\partial s_{nmu}(t)} =1ln2Qnu(t)ϖmx^nmu(t)snmu(t)+x^nmu(t)(I0mu(t)+σ2)hnmu(t)+ηnλ1,nλ2,mhnmu(t)=0\displaystyle=\frac{1}{\ln 2}\frac{Q_{nu}(t)\varpi_{m}\hat{x}_{nmu}(t)}{s_{nmu}(t)+\frac{\hat{x}_{nmu}(t)\left({{I_{0mu}}\left(t\right)+{\sigma^{2}}}\right)}{h_{nmu}(t)}}+\eta_{n}-\lambda_{1,n}-\lambda_{2,m}h_{nmu}(t)=0
x^nmu(t)[0,1],snmu(t)0\displaystyle\hat{x}_{nmu}(t)\in\left[0,1\right],s_{nmu}(t)\geq 0

Thus, we can get the optimal power allocation:

snmu(t)=[1ln2ϖmQnu(t)ωnmuηnI0mu(t)+σ2hnmu(t)]+x^nmu(t)s_{nmu}^{*}(t)=\left[\frac{1}{\ln 2}\frac{\varpi_{m}Q_{nu}(t)}{\omega_{nmu}-\eta_{n}}-\frac{I_{0mu}(t)+\sigma^{2}}{h_{nmu}(t)}\right]^{+}\hat{x}_{nmu}(t) (47)

where ωnmu=λ1,n+λ2,mhnmu(t)\omega_{nmu}=\lambda_{1,n}+\lambda_{2,m}h_{nmu}(t)

Then, we will make use of the results of power allocation for subcarrier assignment. We decompose (45) into UnU_{n} independent subproblems. Each subproblem is formulated as following:

Lnm(𝐏)=u=1UnLnmu(𝐏){L_{nm}}(\mathbf{P})=\sum\limits_{u=1}^{{U_{n}}}{{L_{nmu}}(\mathbf{P})} (48)

where

Lnmu(𝐏)=\displaystyle L_{nmu}(\mathbf{P})= Qnu(t)ϖmx^nmu(t)×log2(1+snmu(t)hnmu(t)x^nmu(t)(I0mu(t)+σ2))+θnmx^nmu(t)+ηnsnmu(t)\displaystyle Q_{nu}(t)\varpi_{m}\hat{x}_{nmu}(t)\times\log_{2}\left(1+\frac{s^{*}_{nmu}(t)h_{nmu}(t)}{\hat{x}_{nmu}(t)(I_{0mu}(t)+\sigma^{2})}\right)+\theta_{nm}\hat{x}_{nmu}(t)+\eta_{n}s^{*}_{nmu}(t)
λ1,nsnmu(t)λ2,msnmu(t)hnmu(t)λ3,mx^nmu(t)\displaystyle-\lambda_{1,n}s^{*}_{nmu}(t)-\lambda_{2,m}s^{*}_{nmu}(t)h_{nmu}(t)-\lambda_{3,m}\hat{x}_{nmu}(t) (49)

Substituting (48) into (50), the objective of subcarrier assignment is to maximize Lnm(𝐏)L_{nm}(\mathbf{P}) for all SUEs associated with SBS nn. For any subcarrier mm, it will be assigned to the SUE who has the biggest Lnmu(𝐏)L_{nmu}(\mathbf{P}). Let mum^{*}_{u} be the result of subcarrier mm’s assignment, which is given by:

mu=argmaxuLmnu and Lmunu>0m_{u}^{*}=\mathop{\arg\max}\limits_{u}{L_{mnu}}\text{ and }{L_{m_{u}^{*}nu}}>0 (50)
xnmu(t)={1ifm=mu0otherwisex_{nmu}^{*}(t)=\left\{\begin{matrix}1&\text{if}\quad m=m^{*}_{u}\\ 0&\text{otherwise}\end{matrix}\right. (51)

We use the subgradient method to update the Lagrange multipliers:

λ1,n(i+1)=[λ1,n(i)d1(i)(pnmaxm=1Mu=1Unsnmu(t))]+\lambda_{1,n}^{(i+1)}=\Bigg{[}\lambda_{1,n}^{(i)}-d_{1}^{(i)}(p_{n}^{max}-\sum\limits_{m=1}^{M}\sum\limits_{u=1}^{U_{n}}s_{nmu}(t))\Bigg{]}^{+} (52)
λ2,m(i+1)=[λ2,m(i)d2(i)(ImSn=1Nu=1Unsnmu(t)hnmu(t))]+\lambda_{2,m}^{(i+1)}=\Bigg{[}\lambda_{2,m}^{(i)}-d_{2}^{(i)}\left(I_{m}^{S}-\sum\limits_{n=1}^{N}\sum\limits_{u=1}^{U_{n}}s_{nmu}(t)h_{nmu}(t)\right)\Bigg{]}^{+} (53)
λ3,m(i+1)=[λ3,m(i)d3(i)(1n=1Nu=1U0xnmu(t))]+\lambda_{3,m}^{(i+1)}=\Bigg{[}\lambda_{3,m}^{(i)}-d_{3}^{(i)}(1-\sum\limits_{n=1}^{N}\sum\limits_{u=1}^{U_{0}}x_{nmu}(t))\Bigg{]}^{+} (54)

where d1(i),d2(i),d3(i)d_{1}^{(i)},d_{2}^{(i)},d_{3}^{(i)} are the step sizes of iteration ii.

Finally, we obtain the proposed spectrum pricing and resource allocation algorithm.

It is noted that, according to [32], the convergence of Eqs. (53)-(55) can be guaranteed by adopting diminishing step sizes dk(i),k=1,2,3d_{k}^{(i)},k=1,2,3. Moreover, since only Step 3 in Algorithm 1 involves iterations. Therefore, the convergence of the proposed algorithm is determined by Eqs. (53)-(55), and by adopting diminishing step sizes dk(i),k=1,2,3d_{k}^{(i)},k=1,2,3, the convergence of the proposed Algorithm 1 can be guaranteed. On the other hand, in practice, the subproblem (47) is solved by each SBS locally for MUnMU_{n} times during one iteration ii and the computational complexity at each SBS is O(MUn)O(MU_{n}).

2) Multi-SBS case (N>2N>2): For the spectrum sharing of multiple SBSs, we decompose the original problem into several two-SBS problems. We first group SBSs into multiple pairs and then use dynamic resource allocation in two-SBS case for each pair. Notice that if the number of SBSs is odd, we can set a virtual SBS with average indexes determined by historical data. From this, we can establish a SBS-pair problem:

maxi=1Nj=1NaijC~ijs.t.i=1Naij=1,j=1Naij=1,aij{0,1}\begin{split}\max&\quad\sum\limits_{i=1}^{N}\sum\limits_{j=1}^{N}a_{ij}\widetilde{C}_{ij}\\ \text{s.t.}&\sum_{i=1}^{N}a_{ij}=1,\sum_{j=1}^{N}a_{ij}=1,a_{ij}\in\left\{0,1\right\}\end{split} (55)

where aija_{ij} is the pair parameter. If SBS ii and SBS jj form a spectrum sharing pair, aij=1a_{ij}=1. Otherwise, aij=0a_{ij}=0. C~ij\widetilde{C}_{ij} is the relative benefit for the SBS ii sharing spectrum with SBS jj, which is compared with the payoff before sharing. The pairing issues can be solved by Hungarian method [33]. Moreover, the complexity of the Hungarian method is O(N4)O(N^{4}). Thus, the overall complexity for each iteration of the proposed online dynamic control algorithm in multi-SBS case is O(MUn+N4)O(MU_{n}+N^{4}).

IV-D Battery Energy Management

By observing the third term on the right-hand side of (73), we can get the harvested energy management problem for each SBS nn:

min\displaystyle\min\,\, (Sn(t)ρn+ϕWn(t))Gn(t)+(Sn(t)ρn)Jn(t)\displaystyle(S_{n}(t)-\rho_{n}+\phi W_{n}(t))G_{n}(t)+(S_{n}(t)-\rho_{n})J_{n}(t)
s.t. Pn(t)=Fn(t)+Gn(t),n,t\displaystyle P_{n}(t)=F_{n}(t)+G_{n}(t),\forall n,t
0Fn(t)Sn(t)\displaystyle 0\leq{F_{n}}\left(t\right)\leq{S_{n}}\left(t\right) (56)
0Jn(t)min{SnmaxSn(t),En(t)}\displaystyle 0\leq J_{n}(t)\leq\min\left\{S_{n}^{max}-S_{n}(t),E_{n}(t)\right\}

The solution of (IV-D) consists of three situations:

a) Sn(t)>ρnS_{n}(t)>\rho_{n}

The storage device does’t harvest renewable energy and provide main energy for transmission. Thus, we have the optimal energy management scheme:

{Jn(t)=0Fn(t)=min{Pn(t),Sn(t)}Gn(t)=0\displaystyle\left\{\begin{matrix}J_{n}^{*}(t)&=&0\\ F_{n}^{*}(t)&=&\min\left\{P_{n}^{*}(t),S_{n}(t)\right\}\\ G_{n}^{*}(t)&=&0\end{matrix}\right. (57)

b) ρnϕWn(t)<Sn(t)ρn\rho_{n}-\phi W_{n}(t)<S_{n}(t)\leq\rho_{n}

The SBS harvests energy to feed the battery and the battery provides main energy for transmission. The energy harvest and power supply scheme is

{Jn(t)=min{SnmaxSn(t),En(t)}Fn(t)=min{Pn(t),Sn(t)}Gn(t)=0\displaystyle\left\{\begin{matrix}J_{n}^{*}(t)&=&\min\{S_{n}^{max}-S_{n}(t),E_{n}(t)\}\\ F_{n}^{*}(t)&=&\min\left\{P_{n}^{*}(t),S_{n}(t)\right\}\\ G_{n}^{*}(t)&=&0\end{matrix}\right. (58)

c) Sn(t)ρnϕWn(t)S_{n}(t)\leq\rho_{n}-\phi W_{n}(t)

The battery level is inadequate for normal transmission then the SBS will be supplied by the grid. The energy scheduling scheme is

{Jn(t)=min{SnmaxSn(t),En(t)}Fn(t)=min{Pn(t),Sn(t)}Gn(t)=max{0,Pn(t)Fn(t)}\displaystyle\left\{\begin{matrix}J_{n}^{*}(t)&=&\min\{S_{n}^{max}-S_{n}(t),E_{n}(t)\}\\ F_{n}^{*}(t)&=&\min\left\{P_{n}^{*}(t),S_{n}(t)\right\}\\ G_{n}^{*}(t)&=&\max\left\{0,P_{n}^{*}(t)-F_{n}^{*}(t)\right\}\end{matrix}\right. (59)

IV-E Discussions on Algorithm Implementation

It is noted that, in the spectrum pricing and resource allocation for multi-SBS case (N>2N>2), the Hungarian method-based pairing may make the proposed algorithm be inefficient for a large-scale or dense LTE/5G network. Then, the zoning or clustering scheme in LTE can be adopted to decrease the complexity of the proposed scheme [20, 34, 35, 36]. Particularly, within a zone, several SBSs bargain on spectrum sharing and then perform the spectrum pricing and resource allocation.

On the other hand, based on the newly emerged software-defined networking (SDN) architecture, which can accelerate the innovations for both hardware forwarding infrastructure and software networking algorithms through control and data separation, enable efficient and adaptive sharing of network resources, and achieve maximum spectrum efficiency and enhance energy efficiency [37], we can designe a virtual green resource allocation and energy management (vGRAEM) scheme, which will reduce the communication overhead over air interface and avoid the information leaking to UEs. Mainly, by leveraging cloud computing and network virtualization, virtual UEs (vUEs) and virtual SBSs (vSBSs) can be generated in the radio access networks controller (RANC) to emulate a resource allocation and energy management solution. The implementation of vGRAEM is shown in Fig. 2. Generally, it consists of three phases. The 1st phase is the initial network measurements; in the 2nd phase, the RANC first generates vUEs and vSBSs and then simulates the dynamic control algorithm for resource allocation and energy management based on the information collected in the 1st phase; in the 3rd phase, the RANC informs individual SBSs about resource allocation and energy management decisions.

Refer to caption
Fig. 2: Implementation of vGRAEM

V Performance Analysis

In this section, we analyse the properties of the queues and proposed online dynamic control algorithm.

Theorem 1: By setting the perturbation ρn\rho_{n} as

ρn=SnmaxEnmax\displaystyle\rho_{n}=S_{n}^{max}-E_{n}^{max} (60)

the energy queue is bounded by 0Sn(t)Snmax0\leq S_{n}(t)\leq S_{n}^{max}.

proof: The proof can be found in Appendix B.

Theorem 2: a) We define f(μn(t)){f\left({\mu_{n}^{*}\left(t\right)}\right)} as our optimal decision and suppose there exist finite and positive constants ε\varepsilon and VV. The proposed online dynamic control algorithm ensures that the data queue and virtual queues have an upper bound.

limT\displaystyle\underset{T\rightarrow\infty}{\mathrm{lim}} sup1Tt=0T1𝔼[n=1Nu=1UnQnu(t)+n=1NYn(t)+n=1NZn(t)]\displaystyle\mathrm{sup}\frac{1}{T}\sum_{t=0}^{T-1}\mathbb{E}\left[{\sum\limits_{n=1}^{N}{\sum\limits_{u=1}^{{U_{n}}}{{Q_{nu}}\left(t\right)}}+\sum\limits_{n=1}^{N}{{Y_{n}}\left(t\right)+\sum\limits_{n=1}^{N}{{Z_{n}}\left(t\right)}}}\right]
H+V(f¯f)ε\displaystyle\leq\frac{{{H}+V\left({\overline{f}-{f^{*}}}\right)}}{{\varepsilon}} (61)

b) The proposed algorithm can achieve the near optimal capacity

f¯fHV\overline{f}\geq{f^{*}}-\frac{{{H}}}{V} (62)

where f¯=limsupT1Tt=0T1V𝔼{n=1Nf(μ(t))}\overline{f}\!\!=\!\!\lim\mathop{\sup}\limits_{T\to\infty}\frac{1}{T}\sum\limits_{t=0}^{T-1}{V\mathbb{E}\left\{{\sum\limits_{n=1}^{N}{f\left({\mu\left(t\right)}\right)}}\right\}}, f=n=1Nf(μn(t)){f^{*}}\!\!=\!\!{\sum\limits_{n=1}^{N}{f\left({\mu_{n}^{*}\left(t\right)}\right)}}

proof: The proof can be found in Appendix C.

VI Simulation Results

We simulate a HetNet where a MBS is underlaid with three uniformly distributed SBSs. The SUEs are randomly distributed in the coverage of their serving SBSs as shown in Fig. 3. The total bandwidth is 30 MHz and we suppose the SBSs divide the spectrum band equally in the initial state. pnmaxp_{n}^{max}=0.1 W, IthS=21010I_{th}^{S}=2*10^{-10}, σ2=BN0\sigma^{2}=BN_{0} where N0=174N_{0}=-174 dBm/Hz is the AWGN power spectral density. The capacity of the SBS storage device is 500 Wh. We simulate the channel as path loss, Rayleigh fading and shadowing effect with mean zero and deviation 10dB. The MBS transmits at its peak power Pm=40P_{m}=40dBm and transmits on each subchannel at the same power. The data arrival is subject to poisson distribution with the mean value 4 packets/slot. The mean packet size is 5000 bits/packet. The static power consumption is 3.2W and the power conversion factor is Δn=4\Delta_{n}=4. The energy harvesting process follows a stationary stochastic process. In addition, we run the simulation for T=1000T=1000 time slots in the Matlab software environment.

Refer to caption
Fig. 3: A macrocell with three coexisting SBSs
Refer to caption
Fig. 4: Dynamics of data queues and energy queues
Refer to caption
Fig. 5: Dynamics of virtual queues Y and Z

Firstly, we observe the queue stability in Fig. 4 and Fig. 5 with V=10V=10. Because all SUEs’ data queues undergo similar trends, we take one SUE of each SBS as example. Fig. 4 shows the dynamics of data queues QQ and energy queues SS. Fig. 5 shows the dynamics of virtual queues YY and ZZ. The upper figure in Fig. 4 demonstrates the results in Theorem 2.a) that the data queues are bounded. The bottom figure in Fig. 4 shows the storage device is below the capacity of storage device according to Theorem 1. Moreover, in Fig. 5, we can also see that virtual queues Yn(t)Y_{n}(t) and Zn(t)Z_{n}(t) are bounded. Hence, Fig. 4 and Fig. 5 show that all the queues are bounded, which means that the network system is stabilized and that the long-term time-average constraints in optimization problem OP2 are satisfied.

Refer to caption
Fig. 6: Average queue backlog versus parameter V

Fig. 6 shows the average network backlog, i.e., 1TUt=0T1n=1Nu=1UnQnu(t)\frac{1}{TU}\!\sum_{t=0}^{T-1}\!\sum_{n=1}^{N}\!\sum_{u=1}^{U_{n}}Q_{nu}(t), versus parameter VV. Compared with no spectrum sharing resource allocation algorithm (NSRA) and time division multiple access resource allocation algorithm (TDRAA), the average queue length in the proposed spectrum sharing algorithm is smaller than that in other algorithms for any VV, which suggests that the proposed spectrum sharing algorithm can reduce network congestion.

Refer to caption
Fig. 7: Performace comparision between proposed algorithm and other algorithms

Fig. 7 shows the network performance comparison between the proposed algorithm and the other two algorithms. Results show that the average network profits obtained by sharing spectrum resources between SBSs are greater than their selfish transmission. Furthermore, the performance of the proposed algorithm is better than other algorithms, which demonstrates its effectiveness and feasibility.

VII conclusion

In this paper, we study the spectrum sharing problem between SBSs powered by renewable energy. We allow the SBS to share free spectrum resources and cooperate with other SBSs for dynamic resource allocation. A multi-bargaining framework is modeled to measure the fairness of sharing profits of SBSs. To solve the problem, we use Lyapunov optimization to decompose the stochastic optimization problem. We develop an online dynamic control algorithm to obtain the optimal transmission power and subchannel assignment and the energy scheduling. Furthermore, the simulation results show the better performance than other algorithms.

Appendix A proof of lemma 1

For any non-negative real numbers xx, yy and zz, there holds (|xy|+z)2x2+y2+z22x(yz)\left(\left|x-y\right|+z\right)^{2}\leq x^{2}+y^{2}+z^{2}-2x\left(y-z\right). We can get the results as follow:

Qnu(t+1)2Qnu(t)2+Rnu(t)2+Dnu(t)2+2Qnu(t)(Dnu(t)Rnu(t))\displaystyle Q_{nu}(t+1)^{2}\leq Q_{nu}(t)^{2}+R_{nu}(t)^{2}+D_{nu}(t)^{2}+2Q_{nu}(t)(D_{nu}(t)-R_{nu}(t)) (63)
Yn(t+1)2Yn(t)2+ynin(t)2+ynout(t)2+2Yn(t)(ynin(t)ynout(t))\displaystyle Y_{n}(t+1)^{2}\leq Y_{n}(t)^{2}+y_{n}^{in}(t)^{2}+y_{n}^{out}(t)^{2}+2Y_{n}(t)(y_{n}^{in}(t)-y_{n}^{out}(t)) (64)
Zn(t+1)2Zn(t)2+znin(t)2+znout(t)2+2Zn(t)(znin(t)znout(t))\displaystyle Z_{n}(t+1)^{2}\leq Z_{n}(t)^{2}+z_{n}^{in}(t)^{2}+z_{n}^{out}(t)^{2}+2Z_{n}(t)(z_{n}^{in}(t)-z_{n}^{out}(t)) (65)

Specially, the energy queue follows the inequality:

(Sn\displaystyle(S_{n} (t+1)ρn)2(Sn(t)ρn)2\displaystyle(t+1)-\rho_{n})^{2}-\left(S_{n}(t)-\rho_{n}\right)^{2}
=(Sn(t+1)+Sn(t)2ρn)(Sn(t+1)Sn(t))\displaystyle=\left(S_{n}(t+1)+S_{n}(t)-2\rho_{n}\right)\left(S_{n}(t+1)-S_{n}(t)\right)
=(2Sn(t)2ρn+Jn(t)Fn(t))(Jn(t)Fn(t))\displaystyle=\left(2S_{n}(t)-2\rho_{n}+J_{n}(t)-F_{n}(t)\right)\left(J_{n}(t)-F_{n}(t)\right)
=2(Sn(t)ρn)(Jn(t)Fn(t))+(Jn(t)Fn(t))2\displaystyle=2\left(S_{n}(t)-\rho_{n}\right)\left(J_{n}(t)-F_{n}(t)\right)+\left(J_{n}(t)-F_{n}(t)\right)^{2}
2(Sn(t)ρn)(Jn(t)Fn(t))+Jn(t)2+Fn(t)2\displaystyle\leq 2\left(S_{n}(t)-\rho_{n}\right)\left(J_{n}(t)-F_{n}(t)\right)+J_{n}(t)^{2}+F_{n}(t)^{2}
2(Sn(t)ρn)(Jn(t)Fn(t))+Jnmax2+Fnmax2\displaystyle\leq 2\left(S_{n}(t)-\rho_{n}\right)\left(J_{n}(t)-F_{n}(t)\right)+J_{n}^{max2}+F_{n}^{max2} (66)

where Jnmax=max{Snmax,Enmax}J_{n}^{max}=\max\{S_{n}^{max},E_{n}^{max}\} and Fnmax=SnmaxF_{n}^{max}=S_{n}^{max}.

By employing above inequalities, we can prove Lemma 1.

Appendix B proof of theorem 1

We initialize the storage device 0Sn(0)Snmax0\leq S_{n}(0)\leq S_{n}^{max} and suppose the limit holds for time slot tt. We define ρn=SnmaxEnmax\rho_{n}=S_{n}^{max}-E_{n}^{max}. Then we will prove it holds for the next time slot t+1t+1. There are two cases to be considered:

(1) 0Sn(t)SnmaxEnmax0\leq S_{n}(t)\leq S_{n}^{max}-E_{n}^{max}

It means Sn(t)ρnS_{n}(t)\leq\rho_{n}, then

Jn(t)=min{SnmaxSn(t),En(t)}=En(t)\displaystyle J_{n}^{*}(t)=\min\{S_{n}^{max}-S_{n}(t),E_{n}(t)\}=E_{n}(t) (67)

so the storage device level at the next time slot is

Sn(t+1)\displaystyle S_{n}(t+1) =[Sn(t)Fn(t)]++Jn(t)Sn(t)+EnmaxSnmax\displaystyle=\left[S_{n}(t)-F_{n}^{*}(t)\right]^{+}+J_{n}^{*}(t)\leq S_{n}(t)+E_{n}^{max}\leq S_{n}^{max} (68)

(2) SnmaxEnmax<S(t)SnmaxS_{n}^{max}-E_{n}^{max}<S(t)\leq S_{n}^{max}

It means Sn(t)<ρnS_{n}(t)<\rho_{n}. In this case, SBS will not accept renewable energy currently, i.e., Jn(t)=0J_{n}^{*}(t)=0 and then

Sn(t+1)=Sn(t)Fn(t)Snmax\displaystyle S_{n}(t+1)=S_{n}(t)-F_{n}^{*}(t)\leq S_{n}^{max} (69)

In summary, Sn(t+1)S_{n}(t+1) is bounded if Sn(t)S_{n}(t) is bounded and this completes the proof.

Appendix C proof of theorem 2

a) Lemma 2. For arbitrary small positive real number ε1{\varepsilon_{1}} and ε2{\varepsilon_{2}}, there exists an algorithm that makes independent, stationary and randomized decisions at each time slot based only on the observed network state, which satisfies

𝔼[Dnu(t)Rnu(t)|Θ(t)]=ε1\mathbb{E}\left[{D_{nu}^{*}\left(t\right)-R_{nu}^{*}\left(t\right)|\Theta\left(t\right)}\right]=-{\varepsilon_{1}} (70)
𝔼[μn(t)+ϕGn(t)+Cnminu=1UnDnu(t)On(t)|Θ(t)]=ε2\mathbb{E}\left[{\mu_{n}^{*}\left(t\right)+\phi G_{n}^{*}\left(t\right)+C_{n}^{\min}-\sum\limits_{u=1}^{{U_{n}}}{D_{nu}^{*}\left(t\right)-{O_{n}^{*}\left(t\right)}|\Theta\left(t\right)}}\right]=-{\varepsilon_{2}} (71)
𝔼[Jn(t)|Θ(t)]=𝔼[Fn(t)|Θ(t)]\mathbb{E}\left[{J_{n}^{*}\left(t\right)|\Theta\left(t\right)}\right]=\mathbb{E}\left[{F_{n}^{*}\left(t\right)|\Theta\left(t\right)}\right] (72)

where Dnu(t),Rnu(t),μn(t),Gn(t),On(t),Jn(t),Fn(t){D_{nu}^{*}\left(t\right),{R_{nu}^{*}\left(t\right)}},{\mu_{n}^{*}\left(t\right)},{G_{n}^{*}\left(t\right)},{O_{n}^{*}\left(t\right)},{J_{n}^{*}\left(t\right)},{F_{n}^{*}\left(t\right)} are corresponding results under the stationary algorithm. The similar proof of Lemma 2. can be found in [31]

Since the algorithm is to minimize the right side of (34) under the constraints, we have

ΔV(t)H+n=1Nu=1UnQnu(t)𝔼{Dnu(t)Rnu(t)|Θ(t)}\displaystyle\Delta_{V}(t)\leq H+\sum\limits_{n=1}^{N}\sum\limits_{u=1}^{U_{n}}Q_{nu}(t)\mathbb{E}\left\{D^{*}_{nu}(t)-R^{*}_{nu}(t)|\Theta(t)\right\}
+n=1NYn(t)𝔼{μn(t)+Cnmin+ϕGn(t)u=1UnDnu(t)On(t)|Θ(t)}\displaystyle+\sum\limits_{n=1}^{N}Y_{n}(t)\mathbb{E}\Bigg{\{}\mu^{*}_{n}(t)+C_{n}^{min}+\phi G^{*}_{n}(t)-\sum_{u=1}^{U_{n}}{D^{*}}_{nu}(t)-O^{*}_{n}(t)|\Theta(t)\Bigg{\}}
+n=1NZn(t)𝔼{Cnmin+ϕGn(t)u=1UnDnu(t)On(t)|Θ(t)}\displaystyle+\sum\limits_{n=1}^{N}Z_{n}(t)\mathbb{E}\Bigg{\{}C_{n}^{min}+\phi G^{*}_{n}(t)-\sum_{u=1}^{U_{n}}{D^{*}}_{nu}(t)-O^{*}_{n}(t)|\Theta(t)\Bigg{\}}
+n=1N{(Sn(t)ρn)𝔼{Jn(t)Fn(t)|Θ(t)}\displaystyle+\sum\limits_{n=1}^{N}\Bigg{\{}\left(S_{n}(t)-\rho_{n}\right)\mathbb{E}\left\{J^{*}_{n}(t)-F^{*}_{n}(t)|\Theta(t)\right\}
Vn=1N𝔼{f(μn(t))|Θ(t)}\displaystyle-V\sum_{n=1}^{N}\mathbb{E}\left\{f(\mu^{*}_{n}(t))|\Theta(t)\right\} (73)

where Dnu(t),Rnu(t),μn(t),Gn(t),On(t),Jn(t),Fn(t){D_{nu}^{*}\left(t\right),{R_{nu}^{*}\left(t\right)}},{\mu_{n}^{*}\left(t\right)},{G_{n}^{*}\left(t\right)},{O_{n}^{*}\left(t\right)},{J_{n}^{*}\left(t\right)},{F_{n}^{*}\left(t\right)} are corresponding results under the stationary algorithm referred to Lemma 2. Substituting (71) (72) and (73) into (74). we have

Δ\displaystyle\Delta (t)VHε1n=1Nu=1UnQnu(t)ε2n=1NYn(t)(ε2+μn(t))n=1NZn(t)Vn=1Nf(μn(t)){}_{V}(t)\leq H-{\varepsilon_{1}}\sum\limits_{n=1}^{N}{\sum\limits_{u=1}^{{U_{n}}}{{Q_{nu}}\left(t\right)}}-{\varepsilon_{2}}\sum\limits_{n=1}^{N}{{Y_{n}}\left(t\right)}-\left({{\varepsilon_{2}}+{\mu^{*}_{n}}\left(t\right)}\right)\sum\limits_{n=1}^{N}{{Z_{n}}\left(t\right)}-V\sum\limits_{n=1}^{N}{f\left({\mu_{n}^{*}\left(t\right)}\right)} (74)

By taking iterated expectation and using telescoping sums over t=0,1,,T1t=0,1,...,T-1, we get

𝔼[L(T)L(0)]t=0T1V𝔼{n=1Nf(μ(t))}\displaystyle\mathbb{E}\left[L(T)-L(0)\right]-\sum\limits_{t=0}^{T-1}{V\mathbb{E}\left\{{\sum\limits_{n=1}^{N}{f\left({\mu\left(t\right)}\right)}}\right\}}
THε1t=0T1𝔼[n=1Nu=1UnQnu(t)]ε2t=0T1𝔼[n=1NYn(t)]\displaystyle\leq TH-{\varepsilon_{1}}\sum\limits_{t=0}^{T-1}{\mathbb{E}\left[{\sum\limits_{n=1}^{N}{\sum\limits_{u=1}^{{U_{n}}}{{Q_{nu}}\left(t\right)}}}\right]}-{\varepsilon_{2}}\sum\limits_{t=0}^{T-1}{\mathbb{E}\left[{\sum\limits_{n=1}^{N}{{Y_{n}}\left(t\right)}}\right]}
(ε2+μn(t))t=0T1𝔼[n=1NZn(t)]Vt=0T1n=1Nf(μn(t))\displaystyle-\left({{\varepsilon_{2}}+{\mu^{*}_{n}}\left(t\right)}\right)\sum\limits_{t=0}^{T-1}{\mathbb{E}\left[{\sum\limits_{n=1}^{N}{{Z_{n}}\left(t\right)}}\right]-V\sum\limits_{t=0}^{T-1}{\sum\limits_{n=1}^{N}{f\left({\mu_{n}^{*}\left(t\right)}\right)}}} (75)

Considering 𝔼{L(T)}>0\mathbb{E}\left\{L(T)\right\}>0. we can find ε=min{ε1,ε2}\varepsilon{\rm{=}}\min\left\{{{\varepsilon_{1}},{\varepsilon_{2}}}\right\}. Then we can get

𝔼[L(T)L(0)]t=0T1V𝔼{n=1Nf(μ(t))}\displaystyle\mathbb{E}\left[L(T)-L(0)\right]-\sum\limits_{t=0}^{T-1}{V\mathbb{E}\left\{{\sum\limits_{n=1}^{N}{f\left({\mu\left(t\right)}\right)}}\right\}}
THεt=0T1𝔼{n=1Nu=1UnQnu(t)+n=1NYn(t)+n=1NZn(t)}Vt=0T1n=1Nf(μn(t))\displaystyle\leq TH-\varepsilon\sum\limits_{t=0}^{T-1}{\mathbb{E}\left\{{\sum\limits_{n=1}^{N}{\sum\limits_{u=1}^{{U_{n}}}{{Q_{nu}}\left(t\right)}}+\sum\limits_{n=1}^{N}{{Y_{n}}\left(t\right)+\sum\limits_{n=1}^{N}{{Z_{n}}\left(t\right)}}}\right\}}-V\sum\limits_{t=0}^{T-1}{\sum\limits_{n=1}^{N}{f\left({\mu_{n}^{*}\left(t\right)}\right)}} (76)

and further get

εt=0T1𝔼{n=1Nu=1UnQnu(t)+n=1NYn(t)+n=1NZn(t)}\displaystyle\varepsilon\sum\limits_{t=0}^{T-1}{\mathbb{E}\left\{{\sum\limits_{n=1}^{N}{\sum\limits_{u=1}^{{U_{n}}}{{Q_{nu}}\left(t\right)}}+\sum\limits_{n=1}^{N}{{Y_{n}}\left(t\right)+\sum\limits_{n=1}^{N}{{Z_{n}}\left(t\right)}}}\right\}}
𝔼[L(0)]+TH+t=0T1VE{n=1Nf(μ(t))}Vt=0T1n=1Nf(μn(t))\displaystyle\leq\mathbb{E}\left[L(0)\right]+TH+\sum\limits_{t=0}^{T-1}{VE\left\{{\sum\limits_{n=1}^{N}{f\left({\mu\left(t\right)}\right)}}\right\}}-V\sum\limits_{t=0}^{T-1}{\sum\limits_{n=1}^{N}{f\left({\mu_{n}^{*}\left(t\right)}\right)}} (77)

Dividing (77) by εT{\varepsilon T} and considering 𝔼[L(0)]=0\mathbb{E}\left[{L\left(0\right)}\right]=0, we have

limT\displaystyle\underset{T\rightarrow\infty}{\mathrm{lim}} sup1Tt=0T1𝔼[n=1Nu=1UnQnu(t)+n=1NYn(t)+n=1NZn(t)]\displaystyle\mathrm{sup}\frac{1}{T}\sum_{t=0}^{T-1}\mathbb{E}\left[{\sum\limits_{n=1}^{N}{\sum\limits_{u=1}^{{U_{n}}}{{Q_{nu}}\left(t\right)}}+\sum\limits_{n=1}^{N}{{Y_{n}}\left(t\right)+\sum\limits_{n=1}^{N}{{Z_{n}}\left(t\right)}}}\right]
THεT+V[t=0T1𝔼{n=1Nf(μn(t))}t=0T1n=1Nf(μn(t))]εT\displaystyle\leq\frac{{TH}}{{\varepsilon T}}+\frac{{V\left[{\sum\limits_{t=0}^{T-1}{\mathbb{E}\left\{{\sum\limits_{n=1}^{N}{f\left({{\mu_{n}}\left(t\right)}\right)}}\right\}-\sum\limits_{t=0}^{T-1}{\sum\limits_{n=1}^{N}{f\left({\mu_{n}^{*}\left(t\right)}\right)}}}}\right]}}{{\varepsilon T}}
=ΔH+V(f¯f)ε\displaystyle\overset{\Delta}{=}\frac{{{H}+V\left({\overline{f}-{f^{*}}}\right)}}{{\varepsilon}} (78)

Hence, the proof of Theorem 2.a) is completed. Similarly, we can prove Theorem 2.b).

References

  • [1] Cisco Visual Networking Index. Global mobile data traffic forecast update, 2015–2020 white paper. link: http://goo. gl/ylTuVx, 2016.
  • [2] Harpreet S. Dhillon, Li Ying, Pavan Nuggehalli, Zhouyue Pi, and Jeffrey G. Andrews. Fundamentals of heterogeneous cellular networks with energy harvesting. IEEE Transactions on Wireless Communications, 13(5):2782–2797, 2014.
  • [3] Sumudu Samarakoon, Mehdi Bennis, Walid Saad, Mérouane Debbah, and Matti Latva-Aho. Ultra dense small cell networks: Turning density into energy efficiency. IEEE Journal on Selected Areas in Communications, 34(5):1267–1280, 2016.
  • [4] Najam Ul Hasan, Waleed Ejaz, Naveed Ejaz, Hyung Seok Kim, and Minho Jo. Network selection and channel allocation for spectrum sharing in 5g heterogeneous networks. IEEE Access, 4:980–992, 2016.
  • [5] Fazlul Kader and Soo Young Shin. Cooperative spectrum sharing with space time block coding and non-orthogonal multiple access. In Eighth International Conference on Ubiquitous & Future Networks, 2016.
  • [6] Johannes Lindblom and Erik G. Larsson. Does non-orthogonal spectrum sharing in the same cell improve the sum-rate of wireless operators? In IEEE International Workshop on Signal Processing Advances in Wireless Communications, 2012.
  • [7] Remco Litjens, Haibin Zhang, Ivo Noppen, Lei Yu, Eleftherios Karipidis, and Borner Kai. System-level assessment of non-orthogonal spectrum sharing via transmit beamforming. In Vehicular Technology Conference, 2013.
  • [8] Francesco Guidolin, Antonino Orsino, Leonardo Badia, and Michele Zorzi. Statistical analysis of non orthogonal spectrum sharing and scheduling strategies in next generation mobile networks. In Wireless Communications & Mobile Computing Conference, 2013.
  • [9] Derrick Wing Kwan Ng, Ernest S. Lo, and Robert Schober. Energy-efficient resource allocation in ofdma systems with hybrid energy harvesting base station. IEEE Transactions on Wireless Communications, 12(7):3412–3427, 2013.
  • [10] Han Tao and Nirwan Ansari. Powering mobile networks with green energy. IEEE Wireless Communications, 21(1):90–96, 2014.
  • [11] P. Palanisamy and S. Nirmala. Downlink interference management in femtocell networks - a comprehensive study and survey. In International Conference on Information Communication & Embedded Systems, 2013.
  • [12] Heli Zhang, Yongbin Wang, and Ji Hong. Resource optimization based interference management for hybrid self-organized small cell network. IEEE Transactions on Vehicular Technology, 65(2):936–946, 2016.
  • [13] Manh Ho Tai, Nguyen H. Tran, Cuong T. Do, S. M. Ahsan Kazmi, Eui Nam Huh, and Choong Seon Hong. Power control for interference management and qos guarantee in heterogeneous networks. IEEE Communications Letters, 19(8):1402–1405, 2015.
  • [14] Hung Sheng Lang, Shih Chun Lin, and Wen Hsien Fang. Subcarrier pairing and power allocation with interference management in cognitive relay networks based on genetic algorithms. IEEE Transactions on Vehicular Technology, 65(9):7051–7063, 2016.
  • [15] Shengrong Bu, F. Richard Yu, and Gamini Senarath. Interference-aware energy-efficient resource allocation for ofdma-based heterogeneous networks with incomplete channel state information. In IEEE International Conference on Communications, 2013.
  • [16] Gong Jie, John S. Thompson, Zhou Sheng, and Zhisheng Niu. Base station sleeping and resource allocation in renewable energy powered cellular networks. Communications IEEE Transactions on, 62(11):3801–3813, 2013.
  • [17] Zhang Yu, Wang Xin, Georgios B Giannakis, and Shuyan Hu. Distributed robust resource allocation for renewable powered wireless cellular networks. In IEEE International Black Sea Conference on Communications & Networking, 2015.
  • [18] Bo Yang, Yanyan Shen, Qiaoni Han, Cailian Chen, Xinping Guan, and Weidong Zhang. Energy-efficient resource allocation for time-varying ofdma relay systems with hybrid energy supplies. IEEE Systems Journal, 12(1):702–713, 2018.
  • [19] Qin Meng, Qinghai Yang, Yang Jian, Daeyoung Park, and Kyung Sup Kwak. Energy-aware resource allocation for ofdma wireless networks with hybrid energy supplies. Iet Communications, 11(11):1671–1678, 2017.
  • [20] Haibo Zhang, Dingde Jiang, Fangwei Li, Houbing Song, and Huaiyu Dai. Cluster-based resource allocation for spectrum-sharing femtocell networks. IEEE Access, 4(99):8643–8656, 2017.
  • [21] Seunghyun Lee and Rui Zhang. Cognitive wireless powered network: Spectrum sharing models and throughput maximization. IEEE Transactions on Cognitive Communications & Networking, 1(3):335–346, 2015.
  • [22] Yong Xiao, Kwang Cheng Chen, Chau Yuen, Zhu Han, and Luiz A. Dasilva. A bayesian overlapping coalition formation game for device-to-device spectrum sharing in cellular networks. IEEE Transactions on Wireless Communications, 14(7):4034–4051, 2015.
  • [23] Bikramjit Singh, Konstantinos Koufos, Olav Tirkkonen, and Randall Berry. Co-primary inter-operator spectrum sharing over a limited spectrum pool using repeated games. In IEEE International Conference on Communications, 2015.
  • [24] K. B. Shashika Manosha, Nandana Rajatheva, and Matti Latva-Aho. Overlay/underlay spectrum sharing for multi-operator environment in cognitive radio networks. In Vehicular Technology Conference, 2011.
  • [25] Huaqing Zhang, Yong Xiao, Lin X. Cai, Dusit Niyato, Lingyang Song, and Zhu Han. A hierarchical game approach for multi-operator spectrum sharing in lte unlicensed. In IEEE Global Communications Conference, 2015.
  • [26] Yashuang Guo, Qinghai Yang, Jiayi Liu, and Kyung Sup Kwak. Cross-layer rate control and resource allocation in spectrum-sharing ofdma small-cell networks with delay constraints. IEEE Transactions on Vehicular Technology, 66(5):4133–4147, 2017.
  • [27] Michael J. Neely. Stochastic network optimization with application to communication and queueing systems. Synthesis Lectures on Communication Networks, 3(1):211, 2010.
  • [28] Gunther Auer, Vito Giannini, Claude Desset, Istvan Godor, Per Skillermark, Magnus Olsson, M. A. Imran, Dario Sabella, Manuel J. Gonzalez, and Oliver Blume. How much energy is needed to run a wireless network? IEEE Wireless Communications, 18(5):40–49, 2012.
  • [29] John F. Nash, Jr. The bargaining problem. Econometrica, 18(2):155–162, 1950.
  • [30] Michael J Neely, Eytan Modiano, and Chih-Ping Li. Fairness and optimal stochastic control for heterogeneous networks. IEEE/ACM Transactions On Networking, 16(2):396–409, 2008.
  • [31] Michael J Neely. Energy optimal control for time-varying wireless networks. IEEE transactions on Information Theory, 52(7):2915–2934, 2006.
  • [32] Johansson Bjorn, Pablo Soldati, and Mikael Johansson. Mathematical decomposition techniques for distributed cross-layer optimization of data networks. IEEE Journal on Selected Areas in Communications, 24(8):1535–1547, 2006.
  • [33] Taejoon Kim and Miaomiao Dong. An iterative hungarian method to joint relay selection and resource allocation for d2d communications. Wireless Communications Letters IEEE, 3(6):625–628, 2016.
  • [34] Mingyi Hong, Ruoyu Sun, Hadi Baligh, and Zhiquan Luo. Joint Base Station Clustering and Beamformer Design for Partial Coordinated Transmission in Heterogeneous Networks. IEEE Journal on Selected Areas in Communications, 31(2):226-240, 2013.
  • [35] Ryuma Seno, Tomoaki Ohtsuki, Wenjie Jiang, and Y Takatori. Complexity reduction of pico cell clustering for interference alignment in heterogeneous networks. In IEEE-APCC, 267-271, 2015.
  • [36] Wanming Hao, Osamu Muta, Haris Gacanin, and Hiroshi Furukawa. Dynamic Small Cell Clustering and Non-Cooperative Game-Based Precoding Design for Two-Tier Heterogeneous Networks With Massive MIMO. IEEE Transactions on Communications, 66(2):675-687, 2018.
  • [37] Qiaoni Han, Bo Yang, Guowang Miao, Cailian Chen, Xiaocheng Wang, and Xinping Wang. Backhaul-Aware User Association and Resource Allocation for Energy-Constrained HetNets. IEEE Transactions on Vehicular Technology, 66(1):580–593, 2017.